1 // YS Subatech Mai 2002
2 // YK Subatech 6 Aug 2002
4 // EMCAL Reconstruction chain:
5 // Hits -> SDigits -> Digits -> RecPoints -> TrackSegments -> RecParticles
11 #include "STEER/AliRun.h"
12 #include "EMCAL/AliEMCALSDigitizer.h"
13 #include "EMCAL/AliEMCALDigitizer.h"
14 #include "EMCAL/AliEMCALClusterizerv1.h"
15 // #include "EMCAL/AliEMCALTrackSegmentMakerv1.h"
16 // #include "EMCAL/AliEMCALPIDv1.h"
18 void EMCALHits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {
21 // 1. write SDigits in the same file as Hits --------------- (OK)
22 //root [0] .L Reconstruct.C++
23 //root [1] SDigits2Digits()
24 // 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
25 //root [0] .L Reconstruct.C++
26 //root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)
31 AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ;
34 sdp->ExecuteTask("deb") ;
39 //________________________________________________________________________
40 void EMCALSDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
43 // 1. write SDigits in the same file as SDigits --------------- (OK)
44 //root [0] .L Reconstruct.C++
45 //root [1] SDigits2Digits()
46 // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
47 //root [0] .L Reconstruct.C++
48 //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=EMCAL, EMCAL)
54 AliEMCALDigitizer * dp = 0 ;
57 dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
58 dp->SetSplitFile() ; }
60 dp = new AliEMCALDigitizer(fileName) ;
62 dp->ExecuteTask("deb") ;
67 //________________________________________________________________________
68 void EMCALDigits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
71 // 1. write RecPoints in the same file as Digits --------------- OK
72 //root [0] .L Reconstruct.C++
73 //root [1] Digits2RecPoints()
74 // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK
75 //root [0] .L Reconstruct.C++
76 //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=EMCAL, EMCAL)
81 AliEMCALClusterizer * cp = 0 ;
84 cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ;
85 cp->SetSplitFile() ; }
87 cp = new AliEMCALClusterizerv1(fileName) ;
89 cp->ExecuteTask("deb") ;
94 // //________________________________________________________________________
95 // void EMCALRecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
98 // // 1. write TrackSegments in the same file as RecPoints --------------- (OK)
99 // //root [0] .L Reconstruct.C++
100 // //root [1] RecPoints2TrackSegments()
101 // // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed)
102 // //root [0] .L Reconstruct.C++
103 // //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
108 // AliEMCALTrackSegmentMaker * tmp = 0 ;
111 // tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ;
113 // tmp = new AliEMCALTrackSegmentMakerv1(fileName) ;
115 // tmp->ExecuteTask("deb") ;
120 // //________________________________________________________________________
121 // void EMCALTrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
124 // // 1. write RecParticles in the same file as TrackSegments --------------- (OK)
125 // //root [0] .L Reconstruct.C++
126 // //root [1] TrackSegments2RecParticles()
127 // // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed)
128 // //root [0] .L Reconstruct.C++
129 // //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
134 // AliEMCALPID * pp = 0 ;
137 // pp = new AliEMCALPIDv1("EMCAL.RecData.root") ;
139 // pp = new AliEMCALPIDv1(fileName) ;
141 // pp->ExecuteTask("deb") ;
146 // //________________________________________________________________________
147 // void EMCALDigits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
150 // // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
151 // //root [0] .L Reconstruct.C++
152 // //root [1] Digits2RecParticles()
153 // // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
154 // //root [0] .L Reconstruct.C++
155 // //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
162 // AliEMCALClusterizer * cp = 0 ;
165 // cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ;
166 // cp->SetSplitFile() ; }
168 // cp = new AliEMCALClusterizerv1(fileName) ;
170 // cp->ExecuteTask("deb") ;
175 // AliEMCALTrackSegmentMaker * tmp = 0 ;
178 // tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ;
180 // tmp = new AliEMCALTrackSegmentMakerv1(fileName) ;
182 // tmp->ExecuteTask("deb") ;
184 // AliEMCALPID * pp = 0 ;
187 // pp = new AliEMCALPIDv1("EMCAL.RecData.root") ;
189 // pp = new AliEMCALPIDv1(fileName) ;
191 // pp->ExecuteTask("deb") ;
197 //________________________________________________________________________
198 void EMCALHits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
200 // 1. write (S)Digits in the same file as Hits --------------- (OK)
201 //root [0] .L Reconstruct.C++
202 //root [1] Hits2Digits()
203 // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
204 //root [0] .L Reconstruct.C++
205 //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)
206 // Digits saved in [DET}.Digits.root (DET=EMCAL, EMCAL)
212 AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ;
214 sdp->SetSplitFile() ;
215 sdp->ExecuteTask("deb") ;
220 AliEMCALDigitizer * dp = 0 ;
223 dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ;
224 dp->SetSplitFile() ; }
226 dp = new AliEMCALDigitizer(fileName) ;
228 dp->ExecuteTask("deb") ;