FAIR Simulation & Analysis FAIR Simulation & Analysis Framework Framework FairRoot FairRoot M. Al-Turany, D. Bertini, F. Uhlig GSI-IT
Jan 02, 2016
FAIR Simulation & Analysis FrameworkFAIR Simulation & Analysis FrameworkFairRoot FairRoot
M. Al-Turany, D. Bertini, F. Uhlig
GSI-IT
12.07.2007 IT-Palaver 2
OverviewOverview
• 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
12.07.2007 IT-Palaver 3
FairRootFairRoot
• 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
12.07.2007 IT-Palaver 4
Geane Integration in FairRoot
• 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
12.07.2007 IT-Palaver 5
What is GEANE? (1)
• 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.
12.07.2007 IT-Palaver 6
What is Geane? (2)What is Geane? (2)
• 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)
12.07.2007 IT-Palaver 7
Geane: Muon Absorber in CBMGeane: Muon Absorber in CBM
12.07.2007 IT-Palaver 8
Geane : Panda detector
12.07.2007 IT-Palaver 9
=1.41
=1.13
=1.19
=0.96
=1.05
Pulls for the whole Panda detector
12.07.2007 IT-Palaver 10
Geane Extrapolation (Hades MDC):Geane Extrapolation (Hades MDC):
• 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)
12.07.2007 IT-Palaver 11
Geane: Hades exampleGeane: Hades example
Hades Simulation
• Geometry: 28 02 2003
• Field map
• Electrons: – 0.05-0.7 GeV– Polar angle range (20., 85.) degree
12.07.2007 IT-Palaver 12
Hades Simulation in FairRootHades Simulation in FairRoot
12.07.2007 IT-Palaver 13
Geane (Red) vs. MC (Black)Geane (Red) vs. MC (Black)
12.07.2007 IT-Palaver 14
Geane (Red) vs. MC (Black)Geane (Red) vs. MC (Black)
12.07.2007 IT-Palaver 15
New Geometry ReaderNew Geometry Reader
• 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.
12.07.2007 IT-Palaver 16
STEP to ROOT STEP to ROOT GeometryGeometry
PANDA MVD detector 9345 volumes
12.07.2007 IT-Palaver 17
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(…)
gconfig/g3Config.C
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)*/ …
gconfig/g4Config.C
Geant4/Geant3 ConfigurationGeant4/Geant3 Configuration
12.07.2007 IT-Palaver 18
Full Simulation-Analysis ChainFull Simulation-Analysis Chain
Event Generator
Transport
Digitizer
Hit Finder
Reconstruction
Physics Analysis
Sim
ula
tio
nA
nal
ysis
Determine particle properties
at target vertex
Transport particles through
the detector material
Determine detector response
Determine physical space point
parameters from detector hits
Determine momentum vector
and PID for all tracks
Calculate physics observables
Storage
Levels
SIM
RAW
12.07.2007 IT-Palaver 19
Fast Simulation-Analysis ChainFast Simulation-Analysis Chain
Event Generator
Put the events on the Stack
Parameterized detector
Response
Physics Analysis
Determine particle properties
at target vertex
No Transport
Calculate physics observables
Storage
Levels
SIM
12.07.2007 IT-Palaver 20
Fast Simulation featuresFast Simulation features
• 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
12.07.2007 IT-Palaver 21
Fast Simulation: MacroFast Simulation: Macro
// 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");
12.07.2007 IT-Palaver 22
Fast Simulation: MacroFast Simulation: Macro
// 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);
12.07.2007 IT-Palaver 23
CMake & CTestCMake & CTest
• 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;
....
12.07.2007 IT-Palaver 24
CMake & CTestCMake & CTest
• 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
12.07.2007 IT-Palaver 25
CMake/CTest : Dashboard CMake/CTest : Dashboard
12.07.2007 IT-Palaver 26
Dashboard : CBMDashboard : CBM
12.07.2007 IT-Palaver 27
The Panda experiment
Multi purpose detector at FAIR
Charmonium (cc)
spectroscopy Open charm spectroscopy Search for gluonic
excitations
(hybrids -
glueballs) Charmed hadrons in nuclei Single and double
Hypernuclei Other options (EFF, GPD,
…)
Physics program
pp, pA collisions1.515 GeV/c (p momentum)
12.07.2007 IT-Palaver 28
PANDA Detector implementation: proposed geometry
12.07.2007 IT-Palaver 29
Detector implementation: state of art
Micro VertexMicro Vertex
TPC/STTTPC/STT
EMC(barrel/Bkw EndCup)
EMC(barrel/Bkw EndCup)
Muon DetectorMuon Detector
EMC(Fwd EndCup)
EMC(Fwd EndCup)
DIRC(Cherenkov)
DIRC(Cherenkov)
COILS(dipole)COILS
(dipole)
COILS(solenoid)
COILS(solenoid)
view from geometry manager
12.07.2007 IT-Palaver 30
Reconstruction example: c in EMC barrel
,C
ClusterizationJan Zhong
Dima Melnichuk
crystal crystal
cluster two clusters
c
°
c at restin lab frame
only barrelnot yet fullcoverage
c
12.07.2007 IT-Palaver 31
Task example - EMC reconstructionFull reconstruction scheme
Full reconstruction scheme
Migrated from Babar-like framework for PandaRoot
Migrated from Babar-like framework for PandaRoot
Dima Melnychuk (Warsaw)
12.07.2007 IT-Palaver 32
Ongoing workOngoing work
• 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
12.07.2007 IT-Palaver 33
AvailabilityAvailability
• 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)