Fixing problems in the generation of HTML documentation
[u/mrichter/AliRoot.git] / EMCAL / AliEMCALReconstruct.C
1 // YS Subatech Mai 2002
2 // YK Subatech 6 Aug 2002
3
4 // EMCAL Reconstruction chain:
5 // Hits -> SDigits -> Digits -> RecPoints -> TrackSegments -> RecParticles
6
7 //Root
8 #if !defined(__CINT__) || defined(__MAKECINT__)
9 #include "TString.h"
10
11 //AliRoot
12 #include "STEER/AliRun.h"
13 #include "EMCAL/AliEMCALSDigitizer.h"
14 #include "EMCAL/AliEMCALDigitizer.h"
15 #include "EMCAL/AliEMCALClusterizerv1.h"
16 // #include "EMCAL/AliEMCALTrackSegmentMakerv1.h"
17 // #include "EMCAL/AliEMCALPIDv1.h"
18 #endif
19
20 void EMCALHits2SDigits( Bool_t split=kFALSE, TString fileName = "galice.root") {
21
22   // usage : 
23   // 1. write SDigits in the same file as Hits --------------- (OK)
24   //root [0] .L Reconstruct.C++
25   //root [1] SDigits2Digits()
26   // 2. write SDigits in a separate file, one per detector, from Hits --------------- (OK)
27   //root [0] .L Reconstruct.C++
28   //root [1] SDigits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)
29
30   delete gAlice ; 
31   gAlice = 0 ; 
32   
33   AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ; 
34   if (split) 
35     sdp->SetSplitFile() ;
36   sdp->ExecuteTask("deb") ; 
37
38   delete sdp ;
39 }
40
41 //________________________________________________________________________
42 void EMCALSDigits2Digits( Bool_t split=kFALSE, TString fileName = "galice.root") {
43   
44  // usage : 
45   // 1. write SDigits in the same file as SDigits --------------- (OK)
46   //root [0] .L Reconstruct.C++
47   //root [1] SDigits2Digits()
48   // 2. write SDigits in a separate file, one per detector, from SDigits --------------- (OK)
49   //root [0] .L Reconstruct.C++
50   //root [1] SDigitsDigits(kTRUE) // Digits saved in [DET}.Digits.root (DET=EMCAL, EMCAL)
51
52   delete gAlice ; 
53   gAlice = 0 ; 
54   
55   // EMCAL
56   AliEMCALDigitizer * dp = 0 ; 
57  
58   if (split) {
59     dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ; 
60     dp->SetSplitFile() ; } 
61   else 
62     dp = new AliEMCALDigitizer(fileName) ; 
63   
64   dp->ExecuteTask("deb") ; 
65   
66   delete dp ;
67 }
68
69 //________________________________________________________________________
70 void EMCALDigits2RecPoints( Bool_t split=kFALSE, TString fileName = "galice.root") {
71   
72  // usage : 
73   // 1. write RecPoints in the same file as Digits --------------- OK 
74   //root [0] .L Reconstruct.C++
75   //root [1] Digits2RecPoints()
76   // 2. write RecPoints in a separate file, one per detector, from Digits --------------- OK 
77   //root [0] .L Reconstruct.C++
78   //root [1] Digits2RecPoints(kTRUE) // RecPoints saved in [DET}.RecPoints.root (DET=EMCAL, EMCAL)
79
80   delete gAlice ; 
81   gAlice = 0 ; 
82  
83   AliEMCALClusterizer * cp = 0 ; 
84  
85   if (split) {
86     cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ; 
87     cp->SetSplitFile() ; } 
88   else 
89     cp = new AliEMCALClusterizerv1(fileName) ; 
90   
91   cp->ExecuteTask("deb") ; 
92   
93   delete cp ;
94 }
95
96 // //________________________________________________________________________
97 // void EMCALRecPoints2TrackSegments( Bool_t split=kFALSE, TString fileName = "galice.root") {
98   
99 //  // usage : 
100 //   // 1. write TrackSegments in the same file as RecPoints --------------- (OK) 
101 //   //root [0] .L Reconstruct.C++
102 //   //root [1] RecPoints2TrackSegments()
103 //   // 2. write TrackSegments in a separate file, one per detector, from RecPoints --------------- (Not needed) 
104 //   //root [0] .L Reconstruct.C++
105 //   //root [1] RecPoints2TrackSegments(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
106
107 //   delete gAlice ; 
108 //   gAlice = 0 ; 
109   
110 //   AliEMCALTrackSegmentMaker * tmp = 0 ; 
111  
112 //   if (split)
113 //     tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ; 
114 //   else 
115 //     tmp = new AliEMCALTrackSegmentMakerv1(fileName) ; 
116   
117 //   tmp->ExecuteTask("deb") ; 
118   
119 //   delete tmp ;
120 // }
121
122 // //________________________________________________________________________
123 // void EMCALTrackSegments2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
124   
125 //  // usage : 
126 //   // 1. write RecParticles in the same file as TrackSegments ---------------  (OK)
127 //   //root [0] .L Reconstruct.C++
128 //   //root [1] TrackSegments2RecParticles()
129 //   // 2. write RecParticles in a separate file, one per detector, from TrackSegments --------------- (Not needed) 
130 //   //root [0] .L Reconstruct.C++
131 //   //root [1] TrackSegments2RecParticles(kTRUE) // RecParticles saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
132
133 //   delete gAlice ; 
134 //   gAlice = 0 ; 
135   
136 //   AliEMCALPID * pp = 0 ; 
137  
138 //   if (split) 
139 //     pp = new AliEMCALPIDv1("EMCAL.RecData.root") ; 
140 //   else 
141 //     pp = new AliEMCALPIDv1(fileName) ; 
142   
143 //   pp->ExecuteTask("deb") ; 
144   
145 //   delete pp ;
146 // }
147
148 // //________________________________________________________________________
149 // void EMCALDigits2RecParticles( Bool_t split=kFALSE, TString fileName = "galice.root") {
150   
151 //  // usage : 
152 //   // 1. write RecPoints, TrackSegments and RecParticles in the same file as Digits --------------- (OK)
153 //   //root [0] .L Reconstruct.C++
154 //   //root [1] Digits2RecParticles()
155 //   // 2. write RecPoints , TrackSegments and RecParticles in a separate file, one per detector, from Digits --------------- (OK)
156 //   //root [0] .L Reconstruct.C++
157 //   //root [1] Digits2RecParticles(kTRUE) // TrackSegments saved in [DET}.RecData.root (DET=EMCAL, EMCAL)
158
159  
160 //   delete gAlice ; 
161 //   gAlice = 0 ; 
162  
163 //   // EMCAL
164 //   AliEMCALClusterizer * cp = 0 ; 
165  
166 //   if (split) {
167 //     cp = new AliEMCALClusterizerv1("EMCAL.Digits.root") ; 
168 //     cp->SetSplitFile() ; } 
169 //   else 
170 //     cp = new AliEMCALClusterizerv1(fileName) ; 
171   
172 //   cp->ExecuteTask("deb") ; 
173
174 //   if (split) 
175 //     delete cp ;
176   
177 //   AliEMCALTrackSegmentMaker * tmp = 0 ; 
178   
179 //   if (split) 
180 //     tmp = new AliEMCALTrackSegmentMakerv1("EMCAL.RecData.root") ; 
181 //   else 
182 //     tmp = new AliEMCALTrackSegmentMakerv1(fileName) ; 
183   
184 //   tmp->ExecuteTask("deb") ; 
185   
186 //   AliEMCALPID * pp = 0 ; 
187  
188 //   if (split) 
189 //     pp = new AliEMCALPIDv1("EMCAL.RecData.root") ; 
190 //   else 
191 //     pp = new AliEMCALPIDv1(fileName) ; 
192   
193 //   pp->ExecuteTask("deb") ; 
194   
195 //   delete tmp; 
196 //   delete pp ; 
197 // }
198
199 //________________________________________________________________________
200 void EMCALHits2Digits (Bool_t split=kFALSE, TString fileName = "galice.root") {
201   // usage : 
202   // 1. write (S)Digits in the same file as Hits --------------- (OK)
203   //root [0] .L Reconstruct.C++
204   //root [1] Hits2Digits()
205   // 2. write (S)Digits in a separate file, one per detector, from Hits --------------- (OK)
206   //root [0] .L Reconstruct.C++
207   //root [1] Hits2Digits(kTRUE) // SDigits saved in [DET}.SDigits.root (DET=EMCAL, EMCAL)
208                                 // Digits  saved in [DET}.Digits.root  (DET=EMCAL, EMCAL)
209
210   delete gAlice ; 
211   gAlice = 0 ; 
212   
213   //EMCAL
214   AliEMCALSDigitizer * sdp = new AliEMCALSDigitizer(fileName) ; 
215   if (split) 
216     sdp->SetSplitFile() ;
217   sdp->ExecuteTask("deb") ; 
218
219   if (split) 
220     delete sdp ; 
221
222   AliEMCALDigitizer * dp = 0 ; 
223  
224   if (split) {
225     dp = new AliEMCALDigitizer("EMCAL.SDigits.root") ; 
226     dp->SetSplitFile() ; } 
227   else 
228     dp = new AliEMCALDigitizer(fileName) ; 
229   
230   dp->ExecuteTask("deb") ; 
231
232   if (split) 
233     delete dp ; 
234
235   if (!split) { 
236     delete sdp ; 
237     delete dp ; 
238   }
239 }
240
241
242