340 likes | 559 Views
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
E N D
FAIR Simulation & Analysis FrameworkFairRoot 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 FAIR experiments design studies CBM PANDA Summary Overview IT-Palaver
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
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
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
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
Geane: Muon Absorber in CBM IT-Palaver
Geane : Panda detector IT-Palaver
Pulls for the whole Panda detector s=1.41 s=1.13 s=1.19 s=1.05 s=0.96 IT-Palaver
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
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
Hades Simulation in FairRoot IT-Palaver
Geane (Red) vs. MC (Black) IT-Palaver
Geane (Red) vs. MC (Black) IT-Palaver
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
STEP to ROOT Geometry PANDA MVD detector 9345 volumes IT-Palaver
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
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
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
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
// 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
// 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
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
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
CMake/CTest : Dashboard IT-Palaver
Dashboard : CBM IT-Palaver
Multi purpose detector at FAIR Physics program pp, pA collisions 1.515 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
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
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
Task example - EMC reconstruction Full reconstruction scheme Migrated from Babar-like framework for PandaRoot Dima Melnychuk (Warsaw) IT-Palaver
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
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