FAIR SimulationAnalysis Framework Fair Root Mohammad AlTurany ITGSI

FAIR Simulation&Analysis Framework Fair. Root Mohammad Al-Turany (IT-GSI) Denis Bertini (IT-GSI) Florian Uhlig (IT-GSI) Ilse König (Hades-GSI) 28. 03. 2007 Root workshop

Overview • Motivations • Fair. Root Features – Geometry Interface – Runtime Database and Parameter Handling – Integrated Track follower (Geane) • FAIR experiments design studies – CBM – PANDA – HADES • Summary 28. 03. 2007 Root workshop 2

Fair@Gsi Improvements: SIS 100/SIS 300 § § Primary beam intensity 100 – 1000 x § Ion energy: Secondary beam intensity for radioactive nucl. : up to 10 000 x 20 x New features: § fast pulsed supraconducting Magnet GSI Accelerator Facilities CBM HESR UNILAC § „cooled“ Antiproton beam up to 15 Ge. V § Specific: intense „cooled“ beam radioactive ‚exotic nucleus Super- FRS Plasma Physics PANDA Atomic Physics FLAIR CR NESR 28. 03. 2007 Root workshop 3

Motivations • Which simulation engine to choose? – Need to move to modern and maintained MC: GEANT 4 – Need for • Working fast ! ( LOI, TDR deadlines … ) • Making reliable simulation – Usually: better knowledge of “old” MC’s: GEANT 3, FLUKA … – A cross-check of simulation results between different MC is needed • Better understanding of GEANT 4 ( intrinsic cuts / physics list …) • Preparing for full simulation • Use of VMC (Virtual Monte Carlo ) : an interface between MCs – With the same code, the user can switch between different MCs 28. 03. 2007 Root workshop 4

Features • The same framework can be used for Simulation and Analysis • Fully ROOT based: – – VMC and VGM for simulation IO scheme (TChain, friend TTrees, TFolders ) for persistency TTask to organize the analysis data flow Use of TGeo. Manager for Simulation and Reconstruction • Completely configurable via ROOT macros • Easy to maintain (only ROOT standard services are used) • Use of a Geometry Interface. – G 3 Native geometry – Geometry Modeller (TGeo. Manager) – Different geometry input format 28. 03. 2007 Root workshop 5

CBM Simulation Geant 3 Geometry Manager Virtual MC Geant 4 ROOT Magnet FLUKA Target PIPE Cave Geo. Interface Module TRD Primary Generator Magnetic Field STS Detector IO Manager Run Manager Tasks Run. Time Data. Base EVGEN TOF RICH Field Map ECAL 28. 03. 2007 Particle Generator ASCII Root workshop Oracle Configuration, Parameters, Geometry Root files Configuration, Parameters, Geometry Pluto Urqmd 6

CBM Analysis ROOT Root files MCPoints, Hits, Digits, Tracks Geometry Manager Geo. Interface Run Manager Module Primary Generator Magnetic Field Detector Tasks Delta 28. 03. 2007 IO Manager digitizers Run. Time Data. Base EVGEN Oracle Configuration, Parameters, Geometry Root files Configuration, Parameters, Geometry Tracking Root workshop 7

Panda Simulation 28. 03. 2007 Root workshop 8

Geometry Interface • Advantage: – more flexibility : different inputs can be used. – closer to technical drawings and analysis coordinate systems • Oracle interface – Hades geometry table design reusable 28. 03. 2007 Root workshop 9

Material & Geometry Interface Oracle DB Cbm. Geo. Interface ASCII G 3 Geometry/Material G 3 Builder Root files 28. 03. 2007 TGeo Geometry/Material TGeo. Volume TGeo. Node TGeo. Material Root. Builder XYZ Geometry/Material XYZ Root workshop 10

Runtime Database The Runtime Database is the manager class for all Parameter containers: Creation, Initialization, Output Runtime Database Container Factories read() 2 Inputs List of Parameter Containers 1 Output List of Runs write() ASCII File ROOT File Oracle 28. 03. 2007 Created in the init() function of the tasks via the container factories or in the macro Filled during initialization IO defined in the macro Root workshop 11

Version management in Oracle • For time dependent information a version management is needed which fulfills the following requirements: – It must be possible to get a consistent set of information for any date (e. g. the start time of a certain run). – To preserve the history, no information - even if wrong - should be overwritten without trace, which means that only inserts should be made, no deletes nor updates. – It must be possible to get an answer to the question: 'Which parameters were used when analyzing this run X years ago? ' (The calibration might have been optimized several times since this date. Maybe some bugs have been detected and corrected in the mean time. ) 28. 03. 2007 Root workshop 12

Version management in Oracle Time dependant entries have a time stamp (date + time with the precision of one second) in form of three columns (Format: DATE): • valid_since : First date when the entry is valid. • valid_until : Last date when the entry is still valid • invalid_since : Date when the entry is replaced by a correct entry or a better version in case of e. g. calibration parameters and therefore gets invalid. 28. 03. 2007 Root workshop 13

Example: Cbm. Generic. Par. Set Base class for most parameter containers Advantage: only a few lines of code to be implemented all I/O interfaces exist already Allows to store various types of parameters (handled by Cbm. Param. List): int, float, double, strings, … arrays TObjects (classes, histograms, …) in derived class implemented 28. 03. 2007 class Cbm. Par. Generic. Set : public Cbm. Par. Set { public: Cbm. Par. Generic. Set(const char* name, const char* title, const char* context) : Cbm. Par. Set(name, title, context) {} virtual ~Cbm. Par. Generic. Set() {} virtual Bool_t init(Cbm. Par. Io*); virtual Int_t write(Cbm. Par. Io*); virtual void put. Params(Cbm. Param. List*)=0; virtual Bool_t get. Params(Cbm. Param. List*)=0; virtual void print. Params(); Root workshop 14

Initialisation scheme (Analysis) Data File=1 Cbm. Task Parameters Cbm. Par. Io Run. Id 1 Sim. Data Cbm. Task: : Set. Containers() Cbm. Task: : init() Par. Cont. Cbm. Task: : Exec() Cbm. Par. Io File=2 Run. Id 2 Cbm. Task: : Reinit() Run. Id 1 Run. Id 2 Par. Cont Sim. Data Cbm. Task: : Exec() 28. 03. 2007 Root workshop 15

Combined Chain& Friend • Make use of our Cbm. Root. Manager (IO): – When the TChain switch to new file: • Clear the global list of friends • Add the correct next friend to the list • Update the corresponding pointers (TTree, Friend Tree. . ) • Combining Chain with Friend: { //… Cbm. Run *f. Run = new Cbm. Run(); f. Run->Set. Input. File(“rich_hit 1. root”); f. Run->Add. Friend(“Rguidance 1. root”); f. Run->Add. File(“rich_hit 2. root”); f. Run->Add. Friend(“Rguidance 2. root”); //… f. Run->Init(); f. Run->Run(); … } 28. 03. 2007 Root workshop 16

Geane Integration in Fair. Root • The integration into the VMC (TGeant 3) is done • In Fair. Root: – Geane can be used in the analysis or from macro – Propagation to – Length – Plane – Volume (Enter or Exit point) – Cbm. Points and/or Cbm. Track. Param can be used as input for propagation – See talk of Andrea Fontana 28. 03. 2007 Root workshop 17

Muon Absorber in CBM 28. 03. 2007 Root workshop 18

CBM experiment@GSI • tracking, vertex reconstruction: radiation hard silicon pixel/strip detectors (STS) in a magnetic dipole field • electron ID: RICH 1 & TRD (& ECAL) p suppression 104 ECAL • hadron ID: TOF (& RICH 2) • photons, p 0, m: ECAL RICHs • high speed DAQ and trigger magnet beam TOF TRDs target 28. 03. 2007 STS Root workshop 19

Analysis: π0 e+e- 28. 03. 2007 Root workshop 20

CBM: KO electrons production in the Au target (setup in the air or in vacuum B=1 T) Geant 3 Geant 4 7. 9 e/HI @ B=1 T 28. 03. 2007 Root workshop 21

Hades@GSI Coils Tof MDC 1, 2 Rich MDC 3 MDC 4 28. 03. 2007 Root workshop Shower 22

Hades example • Need to simulate heavy system at High energy – Need external stack for Geant 3: internal stack capacity reached) – Check data with geant 4 • Easy reuse of Fair. Root framework services • Interesting case ! – Gives us the opportunity to tune Geant 4 (cuts/physics list …) and compare with real Data ! – Realistic test of the framework 28. 03. 2007 Root workshop 23

HADES: MDC Hit distribution Geant 3 Geant 4 Cross-check with Hades data to be done ! 28. 03. 2007 Root workshop 24

The Panda experiment Multi purpose detector at FAIR pp, p. A collisions 1. 5 15 Ge. V/c (p momentum) 28. 03. 2007 Physics program Ø Charmonium (cc) spectroscopy Ø Open charm spectroscopy Ø Search for gluonic excitations (hybrids - glueballs) Ø Charmed hadrons in nuclei Ø Single and double Hypernuclei Ø Other options (EFF, GPD, …) Root workshop 25

PANDA Detector implementation: proposed geometry 28. 03. 2007 Root workshop 26

Detector implementation: state of art view from geometry manager COILS (dipole) Muon Detector DIRC (Cherenkov) TPC/STT Micro Vertex EMC (Fwd End. Cup) 28. 03. 2007 COILS (solenoid) EMC (barrel/Bkw End. Cup) Root workshop 27

Reconstruction example: c in EMC barrel c at rest in lab frame crystal only barrel not yet full coverage cluster c ° two clusters Clusterization Jan Zhong Dima Melnichuk 28. 03. 2007 c Root workshop 28

Summary • A VMC based framework for FAIR has been implemented – First released in March 2004 – Used for FAIR experiments design studies • Ongoing work – Check of Geant 4 VMC interface – Check of TFluka VMC interface – Use more intensively the TGeo. Manager services for the reconstruction 28. 03. 2007 Root workshop 29

Availability • Tested on – – – – – Red Hat 9. 0 (gcc 3. 2. 2 and icc 8. 1) 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 ) Fedora Core 5 (gcc 4. 1. 0) Gentoo (gcc 4. 1. 0) SL 3 -4 64 bit architectures • Binaries are also available for these platforms 28. 03. 2007 Root workshop 30