1 / 33

FAIR Simulation & Analysis Framework FairRoot

FAIR Simulation & Analysis Framework FairRoot . M. Al-Turany, D. Bertini, F. Uhlig GSI-IT. FairRoot FairRoot new features Integrated Track follower (Geane) Geant4 configuration classes New reader (ROOT converted CAD Step Format ) Fast Simulation CMake/CTest

wednesday
Download Presentation

FAIR Simulation & Analysis Framework FairRoot

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. FAIR Simulation & Analysis FrameworkFairRoot M. Al-Turany, D. Bertini, F. Uhlig GSI-IT

  2. FairRoot FairRoot new features Integrated Track follower (Geane) Geant4 configuration classes New reader (ROOT converted CAD Step Format ) Fast Simulation CMake/CTest FAIR experiments design studies CBM PANDA Summary Overview IT-Palaver

  3. FairRoot (Former CbmRoot) has started end of 2003 First released in March 2004 Oct 04 release was used to produce data for the CBM technical report June 05 release ( Hades initialization scheme adapted ) Sept. 06 PANDA collaboration decided to use CbmRoot as simulation and analysis framework Oct. 06 CbmRoot was renamed to FairRoot FairRoot IT-Palaver

  4. The integration into the VMC (TGeant3) is done In FairRoot: Geane can be used in the analysis or from macro Propagation to Length Plane Volume (Enter or Exit point) CbmPoints and/or CbmTrackPar can be used as input for propagation Geane Integration in FairRoot IT-Palaver

  5. Package to calculate the average trajectories of particles through dense materials and to calculate the transport matrix as well as the propagated errors covariance matrix in a given track representation. Geane is a tool to calculate extrapolated track parameters and propagated errors through dense materials. With VMC it is straight forward to use it. What is GEANE? (1) IT-Palaver

  6. It is a track follower: it predicts the trajectory of a charged particle in terms of mean values and errors both in forward and in backward direction. Three effects are taken into account: energy loss (affects mean values and errors) Coulomb multiple scattering (affects errors only) magnetic field (affects mean values only) Geometry and magnetic fields are handled by Geant3 (In VMC applications TGeoManager handles the geometry) What is Geane? (2) IT-Palaver

  7. Geane: Muon Absorber in CBM IT-Palaver

  8. Geane : Panda detector IT-Palaver

  9. Pulls for the whole Panda detector s=1.41 s=1.13 s=1.19 s=1.05 s=0.96 IT-Palaver

  10. Module 1 layer 1(D1S1) To Module 4 layer 6 (D4S6) Propagate MC points (No input errors) Geane::PropagateToVolume is used: Helix track representation Internal representation for Geane (can be transformed) Pull distributions where not calculated : Error are in the volume frame (have to be transformed to lab) Geane Extrapolation (Hades MDC): IT-Palaver

  11. Geane: Hades example • Hades Simulation • Geometry: 28 02 2003 • Field map • Electrons: • 0.05-0.7 GeV • Polar angle range (20., 85.) degree IT-Palaver

  12. Hades Simulation in FairRoot IT-Palaver

  13. Geane (Red) vs. MC (Black) IT-Palaver

  14. Geane (Red) vs. MC (Black) IT-Palaver

  15. A standalone STEP to ROOT geometry converter has been implemented by Tobias Stockmanns (PANDA collaboration) The output ROOT file has the geometry but not the full material properties needed by the simulation engines. A reader has been introduced to the framework that read this geometry and replace the media definition by a proper one. New Geometry Reader IT-Palaver

  16. STEP to ROOT Geometry PANDA MVD detector 9345 volumes IT-Palaver

  17. gconfig/g3Config.C gconfig/g4Config.C geant3->SetPAIR(1); geant3->SetCOMP(1); geant3->SetPHOT(1); geant3->SetPFIS(0); geant3->SetDRAY(1); geant3->SetANNI(1); geant3->SetBREM(1); geant3->SetHADR(3); geant3->SetMUNU(1); geant3->SetDCAY(1); geant3->SetLOSS(1); geant3->SetMULS(1); geant3->SetCKOV(1); geant3->SetRAYL(1); … … geant3->SetCUTS(…) geant4->SetProcess("PAIR",1); /** pair production*/ geant4->SetProcess("COMP",1); /**Compton scattering*/ geant4->SetProcess("PHOT",1); /** photo electric effect */ geant4->SetProcess("PFIS",0); /**photofission*/ geant4->SetProcess("DRAY",1); /**delta-ray*/ geant4->SetProcess("ANNI",1); /**annihilation*/ geant4->SetProcess("BREM",1); /**bremsstrahlung*/ geant4->SetProcess("HADR",1); /**hadronic process*/ geant4->SetProcess("MUNU",1); /**muon nuclear interaction*/ geant4->SetProcess("DCAY",0); /**decay*/ geant4->SetProcess("LOSS",1); /**energy loss*/ geant4->SetProcess("MULS",1); /**multiple scattering*/ geant4->SetProcess("CKOV",1); /**Cerenkov photon generation*/ geant4->SetProcess("RAYL",1); /**Rayleigh scattering*/ ... … geant4->SetCut("CUTGAM",cut1); /** gammas (GeV)*/ … Geant4/Geant3 Configuration IT-Palaver

  18. Full Simulation-Analysis Chain Determine particle properties at target vertex Event Generator Transport particles through the detector material Simulation Transport SIM Digitizer Determine detector response RAW Determine physical space point parameters from detector hits Storage Levels Hit Finder Determine momentum vector and PID for all tracks Analysis Reconstruction Physics Analysis Calculate physics observables IT-Palaver

  19. Fast Simulation-Analysis Chain Determine particle properties at target vertex Event Generator No Transport Put the events on the Stack SIM Parameterized detector Response Storage Levels Physics Analysis Calculate physics observables IT-Palaver

  20. All generators available in FairRoot can be used Events from different generators can be mixed Detector response is implemented in Tasks Give a unified output format for all different generators (CbmStack) Fast analysis can be done with TTree::Draw() The same code can be used to make fast simulation using the full simulation transport files (read only primary particles from stack) The same physics analysis code can be used for fast and full simulation Fast and full simulation results can be easily compared Fast Simulation features IT-Palaver

  21. // Load basic libraries gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C"); basiclibs(); // Load this example libraries gSystem->Load("libGeoBase"); gSystem->Load("libParBase"); gSystem->Load("libBase"); gSystem->Load("libMCStack"); gSystem->Load("libGen"); gSystem->Load("libPassive"); gSystem->Load("libPGen"); CbmRunSim *fRun = new CbmRunSim(); fRun->SetOutputFile("sim_fast.root"); Fast Simulation: Macro IT-Palaver

  22. // Create and Set Event Generator // CbmPrimaryGenerator* primGen = new CbmPrimaryGenerator(); fRun->SetGenerator(primGen); PndDpmGenerator *PndDpm = new PndDpmGenerator("../../input/dpmevt_noelastic_36755.root"); primGen->AddGenerator(PndDpm); /**switch off the transport of particles*/ primGen->DoTracking(kFALSE); fRun->Init(); fRun->Run(1000); Fast Simulation: Macro IT-Palaver

  23. Supports complex, large build environments. CMake has been proved in large projects. (KDE 4) Has powerful commands: include the ability to locate include files, libraries, executables; include external CMake files that encapsulate standard functionality; interfaces to testing systems; supports recursive directory traversal with variable inheritance; can run external programs; supports conditional builds; .... CMake & CTest IT-Palaver

  24. Already in SVN (under testing) Can run in parallel to Automake/Autoconf Need to test: Nightly Builds E-mail to the user who has committed his code Reports See for an example of the web interface: http://www.na-mic.org:8081/Insight/Dashboard/ http://lxg1417.gsi.de:8081 CMake & CTest IT-Palaver

  25. CMake/CTest : Dashboard IT-Palaver

  26. Dashboard : CBM IT-Palaver

  27. Multi purpose detector at FAIR Physics program pp, pA collisions 1.515 GeV/c (p momentum) • Charmonium (cc) spectroscopy • Open charm spectroscopy • Search for gluonic excitations (hybrids - glueballs) • Charmed hadrons in nuclei • Single and double Hypernuclei • Other options (EFF, GPD, …) The Panda experiment IT-Palaver

  28. PANDA Detector implementation: proposed geometry IT-Palaver

  29. Detector implementation: state of art view from geometry manager COILS (dipole) Muon Detector DIRC (Cherenkov) TPC/STT Micro Vertex EMC (Fwd EndCup) COILS (solenoid) EMC (barrel/Bkw EndCup) IT-Palaver

  30. crystal crystal cluster two clusters Reconstruction example: c in EMC barrel c at rest in lab frame c only barrel not yet full coverage c °  Clusterization Jan Zhong Dima Melnichuk IT-Palaver

  31. Task example - EMC reconstruction Full reconstruction scheme Migrated from Babar-like framework for PandaRoot Dima Melnychuk (Warsaw) IT-Palaver

  32. New Geane Interface is implemented and under testing Check of Geant4 VMC interface (Physics list) New G4 physics list for PANDA is needed Energy cuts can be used directly, but the processes has to be adapted to the new list Check of TFluka VMC interface Use more intensively the TGeoManager services for the reconstruction: First tests in PANDA MVD and CBM TOF has been made, and they look very promising Ongoing work IT-Palaver

  33. The framework is available via SVN https://subversion.gsi.de/fairroot/ Tested on Red Hat 9.0 (gcc 3.2.2) Suse 9.0 (gcc 3.3.1) Suse 10.1 (gcc 4.1.0) Debian (gcc 3.2.3) Fedora Core 2 (gcc 3.3.3) Fedora Core 4 (gcc 4.0.0 ) Fedora Core 5 (gcc 4.1.0) Gentoo (gcc 4.1.0) SL 3 & 4 64 bit architectures (Debian, Suse) Availability IT-Palaver

More Related