ECAL in CBMROOT2

1
ECAL in CBMROOT2

• Mikhail Prokudin
• ITEP, MOSCOW

2
Linux flavours tested

• We tested ECAL in CBMROOT2 with
• Debian 3.0 (GSI)
• GCC 3.4.4
• GLibC 2.2.5
• SLC 3.0.6 (ITEP)
• GCC 3.2.3
• GLibC 2.3.2
• Slackware 10.2(my notebook)
• GCC 3.3.6
• GLibC 2.3.5

Ok
Ok
Compilation is Ok. But program fails during ECAL
creation.
3
Whats new

• We can correctly(!!!) construct ECAL in CBMROOT2
  (thanks to Denis)
• CbmEcalInf geometry container class for ECAL.
• CbmEcal?ell ECAL cell class, which knows about
  neighbor cells.
• CbmEcalModule ECAL module class, container for
  ECAL cells.
• ?bmEcalStructure ECAL container class.

4
UML diagram for new classes
CbmEcalModule
CbmEcalCell GetClusterEnergy(Int_t
num) IvansProcedure(Double_t clusterE3x3,
Double_t clusterE2x2, Double_t shape)
CbmEcalStructure (ECAL container) Construct() Rese
tModules() GetCell(Double_t x, Double_t
y) AddEnergy(Double_t x, Double_t y, Double_t
energy, Boot_t isPS)
CbmEcalInf (Geometry container) GetXPos() GetYPos(
) GetZPos() (And a lot of other getters here)
For more complex example see CbmEcalAnalysisDS
class
5
CbmEcalCell and CbmEcalModule

• CbmEcalCell (reco)
• Geometrical properties (position, size, etc.
  accessed via GetCenterX(), GetCenterY())
• Neighbors information (accessed via
  GetNeighborsList(list_num))
• Energy deposition in PS and ECAL
• Total, by all tracks (via GetEnergy(),
  GetPSEnergy(), GetTotalEnergy())
• By any track (via GetTrackEnergy(track_number),)
• Simple cluster reconstruction routines
• Total energy deposition (from all tracks)
• GetClusterEnergy(0) --- energy in 3x3 cluster
• IvansProcedure(E3x3, E2x2, shape)), where E3x3
  --- energy in 3x3 cluster, E2x2 --- energy in 2x2
  subcluster of 3x3 with max. energy,
  shapeE2x2/E3x3
• By specific track (GetTrackClusterEnergy(track_num
  ,0), TrackIvansProcedure(track_num,))
• CbmEcalModule (reco)
• Container for cells, accessed via FindCell(x,y),
  At(x,y)

6
CbmEcalStructure and CbmEcalInf

• CbmEcalStructure (reco)
• Container for ECAL modules and cells
  (GetCell(x,y), GetModule(x,y))
• Constructor for ECAL structure. Creates cells and
  modules with correct neighbors information via
  Construct()
• And some other useful routines
• AddEnergy(x,y,energy, isPreShower) --- find cell
  at x,y, and add energy to it
• ResetModules() --- set energy deposition in all
  modules and cell to 0.
• CbmEcalInf (sim reco)
• Contains all information about ECAL properties
• Geometry (accessed via GetXPos, GetYPos() )
• Granularity (accessed via GetType(Int_t x, Int_t
  y))
• Geant cuts (accessed via GetElectronCut() and
  GetHadronCut())
• Can load such information from text file
  (ecal_v1.geo by default, see next slides)

7
Example of ecal_v1.geo
CbmEcal text file By //Dr.Sys XPos0
Position of ECal center cm YPos0
Position of ECal center cm ZPos1200
Position of ECal start cm NLayers140
Number of layers ModuleSize12 Module size
cm Lead.1 Lead thickness in layer
cm Scin.1 Scintillator thickness in
layercm Tyveec.006 Tyveec thickness in layer
cm PSLead0.86 Lead thickness in preshover
cm PSScin1.00 Scintillator thickness in
layercm PSGap2.0 Gap between ecal and PS
cm CellSize1.0 Cell size in ECAL simulation
cm ECAL will be divided in cells with size
CellSize x CellSize. CbmEcalCell responsible
for real cells assembly from simulated cells.
This procedure allows to construct different
granularities at analysis step from one Geant
simulated input. CbmEcalCells could contain any
(integer) number of CellSize. EcalZSize120.0 Z
size of ECAL container cm At the moment next
two lines don't change anything ECut100e-6
Geant cuts CUTGAM CUTELE BCUTE BCUTM DCUTE
GeV HCut300e-6 Geant cuts CUTNEU CUTHAD
CUTMUO DCUTM PPCUTM GeV
8
Example of ecal_v1.geo. Structure part.
structure Zero for non instrumented area (only
rectangular holes could be constructed). Number
means number of divisions of ECAL module in X and
Y. For example 2 means module divided in 2x24
cells. 1 2 3 4
5 6 7 8 9
234567890123456789012345678901234567890123456789
01234567890123456789012345678901234567890123456789
1111111111111111111111111222222222222222222244444
44444422222222222222222221111111111111111111111111
27 111111111111111111111111122222222222222244444
44444444444444222222222222222111111111111111111111
1111 28 11111111111111111111111122222222222222244
44444444444444444442222222222222221111111111111111
11111111 29 1111111111111111111111112222222222222
44444444444444444444444442222222222222111111111111
111111111111 30 111111111111111111111112222222222
22244444444444444444444444444422222222222221111111
1111111111111111 31 11111111111111111111111222222
22222244444444444444444444444444444222222222222111
11111111111111111111 32 1111111111111111111111122
22222222224444444444444444444444444444422222222222
211111111111111111111111 33 111111111111111111111
12222222222222444444444444444444444444444442222222
2222221111111111111111111111 34
11111111111111111111112222222222224444444444444444
4444444444444442222222222221111111111111111111111
35 1111111111111111111111222222222222444444444444
44444444444444444442222222222221111111111111111111
111 36 111111111111111111111222222222222444444444
44440000000444444444444422222222222211111111111111
1111111 37 11111111111111111111122222222222244444
44444444000000044444444444442222222222221111111111
11111111111 38 1111111111111111111112222222222224
44444444444400000004444444444444222222222222111111
111111111111111 39 111111111111111111111222222222
22244444444444440000000444444444444422222222222211
1111111111111111111 40 11111111111111111111122222
22222224444444444444000000044444444444442222222222
22111111111111111111111 41 1111111111111111111112
22222222222444444444444400000004444444444444222222
222222111111111111111111111 42
11111111111111111111122222222222244444444444440000
0004444444444444222222222222111111111111111111111
43 1111111111111111111111222222222222444444444444
44444444444444444442222222222221111111111111111111
111 44 111111111111111111111122222222222244444444
44444444444444444444444222222222222111111111111111
1111111 45 11111111111111111111112222222222222444
44444444444444444444444444222222222222211111111111
11111111111 46 1111111111111111111111122222222222
24444444444444444444444444444422222222222211111111
111111111111111 47 111111111111111111111112222222
22222444444444444444444444444444442222222222221111
1111111111111111111 48 11111111111111111111111222
22222222224444444444444444444444444442222222222222
11111111111111111111111 49 1111111111111111111111
11222222222222244444444444444444444444442222222222
222111111111111111111111111 50
11111111111111111111111122222222222222244444444444
4444444444222222222222222111111111111111111111111
51 1111111111111111111111111222222222222222444444
44444444444442222222222222221111111111111111111111
111 52 111111111111111111111111122222222222222222
22444444444442222222222222222222111111111111111111
1111111 53
9
System requirements for full MC

• Using 100KeV cuts for e/? and 300KeV for
  hadrons/muons
• Memory requirements for Geant simulation of UrQMD
  central events
• 800? ??? ?BMROOT2
• 750? ??? CBMROOT
• CPU costs (2.8 GHz P4, 1Gb ram) seconds/event

No ECAL ECAL no hit processing Full ECAL
CBMROOT 45 550 600
CBMROOT2 30 250 350
10
ECAL resolution for single photon inCBMROOT and
CBMROOT2
11
Conclusions

• ECAL now can be correctly constructed in CBMROOT2
  (many thanks to Denis)
• A lot of new functionality added and bug
  correction performed
• Need documentation for new classes
• This version of ECAL software will be used for MC
  feasibility studies. We are going to produce at
  least 10K events in ITEP, would be great to run
  some production at GSI.