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)