1 #ifndef ALIANALYSISEMETMONTECARLO_H
2 #define ALIANALYSISEMETMONTECARLO_H
3 //_________________________________________________________________________
4 // Utility Class for transverse energy studies
5 // Base class for MC analysis
8 //*-- Author: Marcelo G. Munhoz (USP)
9 //_________________________________________________________________________
11 #include "AliAnalysisEtMonteCarlo.h"
18 class AliEMCALGeometry;
19 class AliExternalTrackParam;
22 class AliAnalysisEmEtMonteCarlo : public AliAnalysisEtMonteCarlo
27 AliAnalysisEmEtMonteCarlo();
28 virtual ~AliAnalysisEmEtMonteCarlo();
30 virtual Int_t AnalyseEvent(AliVEvent* event);
31 virtual Int_t AnalyseEvent(AliVEvent* event, AliVEvent* event2);
34 virtual void ResetEventValues();
35 virtual void CreateHistograms();
36 virtual void FillOutputList(TList* list);
40 virtual Bool_t IsPrimary(AliStack *stack, Int_t part, TParticlePDG *pdg, Int_t partMom, TParticlePDG *pdgMom);
41 virtual Bool_t IsMotherPrimaryGamma(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom);
42 virtual Bool_t IsMotherPrimaryElectron(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom);
43 virtual Bool_t IsGammaConversion(AliStack *stack, TParticle *part, TParticlePDG *pdg);
44 virtual Bool_t IsInAcceptance(TParticle *part=0, TParticlePDG *pdg=0, AliExternalTrackParam* extParam=0);
45 virtual Bool_t IsInAcceptance(AliMCParticle *part=0);
47 virtual Bool_t TrackHitsCalo(AliExternalTrackParam *extParam);
49 virtual Bool_t GetTrackProjection(AliExternalTrackParam *trackParam, TVector3 &trackPos); // project to a radius
50 virtual Bool_t GetTrackProjection(AliEMCALTrack* emcTrack, TVector3 &trackPos, TVector3 clusPos); // project to a point
52 AliExternalTrackParam* CreateExternalTrackParam(TParticle *part);
54 virtual Double_t CalcET(TParticle *part, TParticlePDG *pdg);
55 virtual Double_t CalcETDep(Double_t caloE, TParticle *part, TParticlePDG *pdg);
59 Double_t fResCut; // b(fm), for Hijing; 0 otherwise
60 Double_t fPrimtotET, fPrimAcctotET, fPrimRectotET, fPrimRectotETDep;//Primary particles ET
62 Double_t fElectrontotET, fElectronAcctotET, fElectronRectotET;//Electron ET
63 Double_t fConvElectrontotET, fConvElectronAcctotET, fConvElectronRectotET, fScatElectrontotET, fScatElectronAcctotET, fScatElectronRectotET;//Secondary electrons ET
64 Double_t fTotElectrontotET, fTotElectronAcctotET, fTotElectronRectotET;//Total Electrons ET
66 Double_t fGammatotET, fGammaAcctotET, fGammaRectotET;//Gamma event ET
67 Double_t fAnnihGammatotET, fAnnihGammaAcctotET, fAnnihGammaRectotET, fScatGammatotET, fScatGammaAcctotET, fScatGammaRectotET;//Secondary gamma ET
68 Double_t fTotGammatotET, fTotGammaAcctotET, fTotGammaRectotET;//Total gamma ET
69 Double_t fConvGammatotET, fNonConvGammatotET, fConvGammaAcctotET, fNonConvGammaAcctotET, fNPPPi0GammatotET, fNPPPi0GammaRectotET;//Conversion and pi0 gamma ET
71 Double_t fTotEMtotET, fTotEMAcctotET, fTotEMRectotET;//Total EM ET
73 Double_t fNPPElectrontotET, fNPPElectronRectotET, fNPPGammatotET, fNPPGammaRectotET;//Non-physical primary electron ET
74 Double_t fTotNPPEMtotET, fTotNPPEMRectotET;//Total Non-physical primary electron ET
76 Double_t fMuontotET, fPiontotET, fKaontotET, fProtontotET;//Charged particles ET
77 Double_t fMuonAcctotET, fPionAcctotET, fKaonAcctotET, fProtonAcctotET;//Charged particles acceptance ET
78 Double_t fMuonRectotET, fMuonRectotETDep, fPionRectotET, fPionRectotETDep, fKaonRectotET, fKaonRectotETDep, fProtonRectotET, fProtonRectotETDep;//Charged particles reconstructed ET
79 Double_t fMuonMatchtotET, fMuonMatchtotETDep, fPionMatchtotET, fPionMatchtotETDep, fKaonMatchtotET, fKaonMatchtotETDep, fProtonMatchtotET, fProtonMatchtotETDep;//Charged particles track matched ET
80 Double_t fTotChargedtotET, fTotChargedAcctotET, fTotChargedRectotET, fTotChargedRectotETDep, fTotChargedMatchtotET, fTotChargedMatchtotETDep;//Total charged particles ET
82 Double_t fNeutrontotET, fNeutronAcctotET, fNeutronRectotET, fNeutronRectotETDep;//Neutrons ET
83 Double_t fK0totET, fK0RectotET, fK0RectotETDep, fLambdatotET, fLambdaRectotET, fLambdaRectotETDep;//K0 and Lambda ET
84 Double_t fTotNeutraltotET, fTotNeutralRectotET, fTotNeutralRectotETDep;//Neutral particles ET
86 Double_t fTotaltotET, fTotalAcctotET, fTotalRectotET, fTotalRectotETDep;//Total ET
88 AliEMCALGeometry *fGeoUt;//EMCal geometry object
90 // *******************
92 // *******************
93 TH2F *fHistPrimEtaEET;//total ET - Eta vs E
94 TH2F *fHistPrimEtaPtET;//total ET - Eta vs pt
95 TH2F *fHistPrimEtaET;//total ET - Eta
96 TH1F *fHistPrimtotET;//total ET distribution
98 TH2F *fHistPrimAccEtaEET;//acceptance ET - Eta vs E
99 TH2F *fHistPrimAccEtaPtET;//acceptance ET - Eta vs pt
100 TH2F *fHistPrimAccEtaET;//acceptance ET - Eta
101 TH1F *fHistPrimAcctotET;//acceptance ET distribution
103 TH2F *fHistPrimRecEtaEET;//reconstructed ET - Eta vs E
104 TH2F *fHistPrimRecEtaPtET;//reconstructed ET - Eta vs pt
105 TH2F *fHistPrimRecEtaET;//reconstructed ET - Eta
106 TH1F *fHistPrimRectotET;//reconstructed ET distribution
108 TH2F *fHistPrimRecEtaEDepETDep;//deposited ET - Eta vs E deposited
109 TH2F *fHistPrimRecEtaPtETDep;//deposited ET - Eta vs pt
110 TH2F *fHistPrimRecEtaETDep;//deposited ET - Eta
111 TH1F *fHistPrimRectotETDep;//deposited ET distribution
113 // *******************
115 // *******************
116 TH2F *fHistElectronEtaEET;// ET - Eta vs E
117 TH2F *fHistElectronEtaPtET;//ET - Eta vs pt
118 TH2F *fHistElectronEtaET;// ET - Eta
119 TH2F *fHistElectronEtaE;// multiplicity - Eta vs E
120 TH2F *fHistElectronEtaPt;// multiplicity - Eta vs pt
121 TH1F *fHistElectrontotET;// total ET distribution
123 TH2F *fHistConvElectronEtaEET;//ET - Eta vs E
124 TH2F *fHistConvElectronEtaPtET;//ET - Eta vs pt
125 TH2F *fHistConvElectronEtaET;// ET - Eta
126 TH2F *fHistConvElectronEtaE;// multiplicity - Eta vs E
127 TH2F *fHistConvElectronEtaPt;// multiplicity - Eta vs pt
128 TH1F *fHistConvElectrontotET;// total ET distribution
130 TH2F *fHistScatElectronEtaEET;//ET - Eta vs E
131 TH2F *fHistScatElectronEtaPtET;//ET - Eta vs pt
132 TH2F *fHistScatElectronEtaET;// ET - Eta
133 TH2F *fHistScatElectronEtaE;// multiplicity - Eta vs E
134 TH2F *fHistScatElectronEtaPt;// multiplicity - Eta vs pt
135 TH1F *fHistScatElectrontotET;// total ET distribution
137 // *******************
139 // *******************
140 TH1F *fHistTotElectrontotET;//total ET distribution
142 // *******************
144 // *******************
145 TH2F *fHistGammaEtaEET;//ET - Eta vs E
146 TH2F *fHistGammaEtaPtET;//ET - Eta vs pt
147 TH2F *fHistGammaEtaET;// ET - Eta
148 TH2F *fHistGammaEtaE;// multiplicity - Eta vs E
149 TH2F *fHistGammaEtaPt;// multiplicity - Eta vs pt
150 TH1F *fHistGammatotET;// total ET distribution
152 TH2F *fHistAnnihGammaEtaEET;//ET - Eta vs E
153 TH2F *fHistAnnihGammaEtaPtET;//ET - Eta vs pt
154 TH2F *fHistAnnihGammaEtaET;// ET - Eta
155 TH2F *fHistAnnihGammaEtaE;// multiplicity - Eta vs E
156 TH2F *fHistAnnihGammaEtaPt;// multiplicity - Eta vs pt
157 TH1F *fHistAnnihGammatotET;// total ET distribution
159 TH2F *fHistScatGammaEtaEET;//ET - Eta vs E
160 TH2F *fHistScatGammaEtaPtET;//ET - Eta vs pt
161 TH2F *fHistScatGammaEtaET;// ET - Eta
162 TH2F *fHistScatGammaEtaE;// multiplicity - Eta vs E
163 TH2F *fHistScatGammaEtaPt;// multiplicity - Eta vs pt
164 TH1F *fHistScatGammatotET;// total ET distribution
166 TH2F *fHistConvGammaEtaEET;//ET - Eta vs E
167 TH2F *fHistConvGammaEtaPtET;//ET - Eta vs pt
168 TH2F *fHistConvGammaEtaET;// ET - Eta
169 TH2F *fHistConvGammaEtaE;// multiplicity - Eta vs E
170 TH2F *fHistConvGammaEtaPt;// multiplicity - Eta vs pt
171 TH1F *fHistConvGammatotET;// total ET distribution
173 TH2F *fHistNonConvGammaEtaEET;//ET - Eta vs E
174 TH2F *fHistNonConvGammaEtaPtET;//ET - Eta vs pt
175 TH2F *fHistNonConvGammaEtaET;// ET - Eta
176 TH2F *fHistNonConvGammaEtaE;// multiplicity - Eta vs E
177 TH2F *fHistNonConvGammaEtaPt;// multiplicity - Eta vs pt
178 TH1F *fHistNonConvGammatotET;// total ET distribution
180 // *******************
182 // *******************
183 TH1F *fHistTotGammatotET;//total ET distribution
185 // *******************
186 // total electromagnetic ET
187 // *******************
188 TH1F *fHistTotEMtotET;//total ET distribution
190 // non-primary electromagnetic ET
191 TH2F *fHistNPPElectronEtaEET;//ET - Eta vs E
192 TH2F *fHistNPPElectronEtaPtET;//ET - Eta vs pt
193 TH2F *fHistNPPElectronEtaET;// ET - Eta
194 TH2F *fHistNPPElectronEtaE;// multiplicity - Eta vs E
195 TH2F *fHistNPPElectronEtaPt;// multiplicity - Eta vs pt
196 TH1F *fHistNPPElectrontotET;// total ET distribution
198 TH2F *fHistNPPGammaEtaEET;//ET - Eta vs E
199 TH2F *fHistNPPGammaEtaPtET;//ET - Eta vs pt
200 TH2F *fHistNPPGammaEtaET;// ET - Eta
201 TH2F *fHistNPPGammaEtaE;// multiplicity - Eta vs E
202 TH2F *fHistNPPGammaEtaPt;// multiplicity - Eta vs pt
203 TH1F *fHistNPPGammatotET;// total ET distribution
205 TH1F *fHistTotNPPEMtotET;//total ET distribution
207 TH2F *fHistNPPPi0GammaEtaEET;//ET - Eta vs E
208 TH2F *fHistNPPPi0GammaEtaPtET;//ET - Eta vs pt
209 TH2F *fHistNPPPi0GammaEtaET;// ET - Eta
210 TH2F *fHistNPPPi0GammaEtaE;// multiplicity - Eta vs E
211 TH2F *fHistNPPPi0GammaEtaPt;// multiplicity - Eta vs pt
212 TH1F *fHistNPPPi0GammatotET;// total ET distribution
214 // *******************
215 // electron ET inside EMCal acceptance
216 // *******************
217 TH2F *fHistElectronAccEtaEET;//ET - Eta vs E
218 TH2F *fHistElectronAccEtaPtET;//ET - Eta vs pt
219 TH2F *fHistElectronAccEtaET;// ET - Eta
220 TH2F *fHistElectronAccEtaE;// multiplicity - Eta vs E
221 TH2F *fHistElectronAccEtaPt;// multiplicity - Eta vs pt
222 TH1F *fHistElectronAcctotET;// total ET distribution
224 TH2F *fHistConvElectronAccEtaEET;//ET - Eta vs E
225 TH2F *fHistConvElectronAccEtaPtET;//ET - Eta vs pt
226 TH2F *fHistConvElectronAccEtaET;// ET - Eta
227 TH2F *fHistConvElectronAccEtaE;// multiplicity - Eta vs E
228 TH2F *fHistConvElectronAccEtaPt;// multiplicity - Eta vs pt
229 TH1F *fHistConvElectronAcctotET;// total ET distribution
231 TH2F *fHistScatElectronAccEtaEET;//ET - Eta vs E
232 TH2F *fHistScatElectronAccEtaPtET;//ET - Eta vs pt
233 TH2F *fHistScatElectronAccEtaET;// ET - Eta
234 TH2F *fHistScatElectronAccEtaE;// multiplicity - Eta vs E
235 TH2F *fHistScatElectronAccEtaPt;// multiplicity - Eta vs pt
236 TH1F *fHistScatElectronAcctotET;// total ET distribution
238 // *******************
239 // total electron ET inside EMCal acceptance
240 // *******************
241 TH1F *fHistTotElectronAcctotET;//total ET distribution
243 // *******************
244 // gamma ET inside EMCal acceptance
245 // *******************
246 TH2F *fHistGammaAccEtaEET;//ET - Eta vs E
247 TH2F *fHistGammaAccEtaPtET;//ET - Eta vs pt
248 TH2F *fHistGammaAccEtaET;// ET - Eta
249 TH2F *fHistGammaAccEtaE;// multiplicity - Eta vs E
250 TH2F *fHistGammaAccEtaPt;// multiplicity - Eta vs pt
251 TH1F *fHistGammaAcctotET;// total ET distribution
253 TH2F *fHistAnnihGammaAccEtaEET;//ET - Eta vs E
254 TH2F *fHistAnnihGammaAccEtaPtET;//ET - Eta vs pt
255 TH2F *fHistAnnihGammaAccEtaET;// ET - Eta
256 TH2F *fHistAnnihGammaAccEtaE;// multiplicity - Eta vs E
257 TH2F *fHistAnnihGammaAccEtaPt;// multiplicity - Eta vs pt
258 TH1F *fHistAnnihGammaAcctotET;// total ET distribution
260 TH2F *fHistScatGammaAccEtaEET;//ET - Eta vs E
261 TH2F *fHistScatGammaAccEtaPtET;//ET - Eta vs pt
262 TH2F *fHistScatGammaAccEtaET;// ET - Eta
263 TH2F *fHistScatGammaAccEtaE;// multiplicity - Eta vs E
264 TH2F *fHistScatGammaAccEtaPt;// multiplicity - Eta vs pt
265 TH1F *fHistScatGammaAcctotET;// total ET distribution
267 TH2F *fHistConvGammaAccEtaEET;//ET - Eta vs E
268 TH2F *fHistConvGammaAccEtaPtET;//ET - Eta vs pt
269 TH2F *fHistConvGammaAccEtaET;// ET - Eta
270 TH2F *fHistConvGammaAccEtaE;// multiplicity - Eta vs E
271 TH2F *fHistConvGammaAccEtaPt;// multiplicity - Eta vs pt
272 TH1F *fHistConvGammaAcctotET;// total ET distribution
274 TH2F *fHistNonConvGammaAccEtaEET;//ET - Eta vs E
275 TH2F *fHistNonConvGammaAccEtaPtET;//ET - Eta vs pt
276 TH2F *fHistNonConvGammaAccEtaET;// ET - Eta
277 TH2F *fHistNonConvGammaAccEtaE;// multiplicity - Eta vs E
278 TH2F *fHistNonConvGammaAccEtaPt;// multiplicity - Eta vs pt
279 TH1F *fHistNonConvGammaAcctotET;// total ET distribution
281 // *******************
282 // total gamma ET inside EMCal acceptance
283 // *******************
284 TH1F *fHistTotGammaAcctotET;//total ET distribution
286 // *******************
287 // total electromagnetic ET inside EMCal acceptance
288 // *******************
289 TH1F *fHistTotEMAcctotET;//total ET distribution
291 // non-primary electromagnetic ET
292 TH2F *fHistNPPElectronAccEtaEET;//ET - Eta vs E
293 TH2F *fHistNPPElectronAccEtaPtET;//ET - Eta vs pt
294 TH2F *fHistNPPElectronAccEtaE;// multiplicity - Eta vs E
295 TH2F *fHistNPPElectronAccEtaPt;// multiplicity - Eta vs pt
297 TH2F *fHistNPPGammaAccEtaEET;//ET - Eta vs E
298 TH2F *fHistNPPGammaAccEtaPtET;//ET - Eta vs pt
299 TH2F *fHistNPPGammaAccEtaE;// multiplicity - Eta vs E
300 TH2F *fHistNPPGammaAccEtaPt;// multiplicity - Eta vs pt
302 // *******************
303 // electron ET reconstructed in EMCal
304 // *******************
305 TH2F *fHistElectronRecEtaEET;//ET - Eta vs E
306 TH2F *fHistElectronRecEtaPtET;//ET - Eta vs pt
307 TH2F *fHistElectronRecEtaET;// ET - Eta
308 TH2F *fHistElectronRecEtaE;// multiplicity - Eta vs E
309 TH2F *fHistElectronRecEtaPt;// multiplicity - Eta vs pt
310 TH1F *fHistElectronRectotET;// total ET distribution
312 TH2F *fHistConvElectronRecEtaEET;//ET - Eta vs E
313 TH2F *fHistConvElectronRecEtaPtET;//ET - Eta vs pt
314 TH2F *fHistConvElectronRecEtaET;// ET - Eta
315 TH2F *fHistConvElectronRecEtaE;// multiplicity - Eta vs E
316 TH2F *fHistConvElectronRecEtaPt;// multiplicity - Eta vs pt
317 TH1F *fHistConvElectronRectotET;// total ET distribution
319 TH2F *fHistScatElectronRecEtaEET;//ET - Eta vs E
320 TH2F *fHistScatElectronRecEtaPtET;//ET - Eta vs pt
321 TH2F *fHistScatElectronRecEtaET;// ET - Eta
322 TH2F *fHistScatElectronRecEtaE;// multiplicity - Eta vs E
323 TH2F *fHistScatElectronRecEtaPt;// multiplicity - Eta vs pt
324 TH1F *fHistScatElectronRectotET;// total ET distribution
326 // *******************
327 // total Electron ET reconstructed in EMCal
328 // *******************
329 TH1F *fHistTotElectronRectotET;//total ET distribution
331 // *******************
332 // gamma ET reconstructed in EMCal
333 // *******************
334 TH2F *fHistGammaRecEtaEET;//ET - Eta vs E
335 TH2F *fHistGammaRecEtaPtET;//ET - Eta vs pt
336 TH2F *fHistGammaRecEtaET;// ET - Eta
337 TH2F *fHistGammaRecEtaE;// multiplicity - Eta vs E
338 TH2F *fHistGammaRecEtaPt;// multiplicity - Eta vs pt
339 TH1F *fHistGammaRectotET;// total ET distribution
341 TH2F *fHistAnnihGammaRecEtaEET;//ET - Eta vs E
342 TH2F *fHistAnnihGammaRecEtaPtET;//ET - Eta vs pt
343 TH2F *fHistAnnihGammaRecEtaET;// ET - Eta
344 TH2F *fHistAnnihGammaRecEtaE;// multiplicity - Eta vs E
345 TH2F *fHistAnnihGammaRecEtaPt;// multiplicity - Eta vs pt
346 TH1F *fHistAnnihGammaRectotET;// total ET distribution
348 TH2F *fHistScatGammaRecEtaEET;//ET - Eta vs E
349 TH2F *fHistScatGammaRecEtaPtET;//ET - Eta vs pt
350 TH2F *fHistScatGammaRecEtaET;// ET - Eta
351 TH2F *fHistScatGammaRecEtaE;// multiplicity - Eta vs E
352 TH2F *fHistScatGammaRecEtaPt;// multiplicity - Eta vs pt
353 TH1F *fHistScatGammaRectotET;// total ET distribution
355 // *******************
356 // total gamma ET reconstructed in EMCal
357 // *******************
358 TH1F *fHistTotGammaRectotET;//total ET distribution
360 // *******************
361 // total EM ET reconstructed in EMCal
362 // *******************
363 TH1F *fHistTotEMRectotET;//total ET distribution
365 // non-primary electromagnetic ET
366 TH2F *fHistNPPElectronRecEtaEET;//ET - Eta vs E
367 TH2F *fHistNPPElectronRecEtaPtET;//ET - Eta vs pt
368 TH2F *fHistNPPElectronRecEtaET;// ET - Eta
369 TH2F *fHistNPPElectronRecEtaE;// multiplicity - Eta vs E
370 TH2F *fHistNPPElectronRecEtaPt;// multiplicity - Eta vs pt
371 TH1F *fHistNPPElectronRectotET;// total ET distribution
373 TH2F *fHistNPPGammaRecEtaEET;//ET - Eta vs E
374 TH2F *fHistNPPGammaRecEtaPtET;//ET - Eta vs pt
375 TH2F *fHistNPPGammaRecEtaET;// ET - Eta
376 TH2F *fHistNPPGammaRecEtaE;// multiplicity - Eta vs E
377 TH2F *fHistNPPGammaRecEtaPt;// multiplicity - Eta vs pt
378 TH1F *fHistNPPGammaRectotET;// total ET distribution
380 TH1F *fHistTotNPPEMRectotET;//total ET distribution
382 TH2F *fHistNPPPi0GammaRecEtaEET;//ET - Eta vs E
383 TH2F *fHistNPPPi0GammaRecEtaPtET;//ET - Eta vs pt
384 TH2F *fHistNPPPi0GammaRecEtaET;// ET - Eta
385 TH2F *fHistNPPPi0GammaRecEtaE;// multiplicity - Eta vs E
386 TH2F *fHistNPPPi0GammaRecEtaPt;// multiplicity - Eta vs pt
387 TH1F *fHistNPPPi0GammaRectotET;// total ET distribution
389 // *******************
391 // *******************
392 TH2F *fHistMuonEtaEET;//ET - Eta vs E
393 TH2F *fHistMuonAccEtaEET;//ET - Eta vs E
394 TH2F *fHistMuonRecEtaEET;//ET - Eta vs E
395 TH2F *fHistMuonMatchEtaEET;//ET - Eta vs E
397 TH2F *fHistMuonEtaPtET;// ET - Eta vs pt
398 TH2F *fHistMuonAccEtaPtET;// ET - Eta vs pt
399 TH2F *fHistMuonRecEtaPtET;// ET - Eta vs pt
400 TH2F *fHistMuonMatchEtaPtET;// ET - Eta vs pt
402 TH2F *fHistMuonEtaET;// ET - Eta
403 TH2F *fHistMuonAccEtaET;// ET - Eta
404 TH2F *fHistMuonRecEtaET;// ET - Eta
405 TH2F *fHistMuonMatchEtaET;// ET - Eta
407 TH2F *fHistMuonEtaE;// multiplicity - Eta vs E
408 TH2F *fHistMuonAccEtaE;// multiplicity - Eta vs E
409 TH2F *fHistMuonRecEtaE;// multiplicity - Eta vs E
410 TH2F *fHistMuonMatchEtaE;// multiplicity - Eta vs E
412 TH2F *fHistMuonEtaPt;// multiplicity - Eta vs pt
413 TH2F *fHistMuonAccEtaPt;// multiplicity - Eta vs pt
414 TH2F *fHistMuonRecEtaPt;// multiplicity - Eta vs pt
415 TH2F *fHistMuonMatchEtaPt;// multiplicity - Eta vs pt
417 TH1F *fHistMuontotET;// total ET distribution
418 TH1F *fHistMuonAcctotET;// total ET distribution
419 TH1F *fHistMuonRectotET;// total ET distribution
420 TH1F *fHistMuonMatchtotET;// total ET distribution
422 TH1F *fHistMuonRectotETDep;//total deposited ET distribution
423 TH1F *fHistMuonMatchtotETDep;// total deposited ET distribution
425 TH2F *fHistMuonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
426 TH2F *fHistMuonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
428 TH2F *fHistMuonRecEtaPtETDep;// ET deposited - Eta vs pt
429 TH2F *fHistMuonMatchEtaPtETDep;// ET deposited - Eta vs pt
431 TH2F *fHistMuonRecEtaETDep;// ET deposited - Eta
432 TH2F *fHistMuonMatchEtaETDep;// ET deposited - Eta
434 TH2F *fHistMuonRecResEET;// ET - track matching residual vs E
435 TH2F *fHistMuonRecResPtET;// ET - track matching residual vs pt
436 TH2F *fHistMuonRecResE;// multiplicity - track matching residual vs E
437 TH2F *fHistMuonRecResPt;// multiplicity - track matching residual vs pt
438 TH2F *fHistMuonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
439 TH2F *fHistMuonRecResPtETDep;// ET deposited - track matching residual vs pt
441 // *******************
443 // *******************
444 TH2F *fHistPionEtaEET;//ET - Eta vs E
445 TH2F *fHistPionAccEtaEET;//ET - Eta vs E
446 TH2F *fHistPionRecEtaEET;//ET - Eta vs E
447 TH2F *fHistPionMatchEtaEET;//ET - Eta vs E
449 TH2F *fHistPionEtaPtET;// ET - Eta vs pt
450 TH2F *fHistPionAccEtaPtET;// ET - Eta vs pt
451 TH2F *fHistPionRecEtaPtET;// ET - Eta vs pt
452 TH2F *fHistPionMatchEtaPtET;// ET - Eta vs pt
454 TH2F *fHistPionEtaET;// ET - Eta
455 TH2F *fHistPionAccEtaET;// ET - Eta
456 TH2F *fHistPionRecEtaET;// ET - Eta
457 TH2F *fHistPionMatchEtaET;// ET - Eta
459 TH2F *fHistPionEtaE;// multiplicity - Eta vs E
460 TH2F *fHistPionAccEtaE;// multiplicity - Eta vs E
461 TH2F *fHistPionRecEtaE;// multiplicity - Eta vs E
462 TH2F *fHistPionMatchEtaE;// multiplicity - Eta vs E
464 TH2F *fHistPionEtaPt;// multiplicity - Eta vs pt
465 TH2F *fHistPionAccEtaPt;// multiplicity - Eta vs pt
466 TH2F *fHistPionRecEtaPt;// multiplicity - Eta vs pt
467 TH2F *fHistPionMatchEtaPt;// multiplicity - Eta vs pt
469 TH1F *fHistPiontotET;// total ET distribution
470 TH1F *fHistPionAcctotET;// total ET distribution
471 TH1F *fHistPionRectotET;// total ET distribution
472 TH1F *fHistPionMatchtotET;// total ET distribution
474 TH1F *fHistPionRectotETDep;// total deposited ET distribution
475 TH1F *fHistPionMatchtotETDep;// total deposited ET distribution
477 TH2F *fHistPionRecEtaEDepETDep;// ET deposited - Eta vs E deposited
478 TH2F *fHistPionMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
480 TH2F *fHistPionRecEtaPtETDep;// ET deposited - Eta vs pt
481 TH2F *fHistPionMatchEtaPtETDep;// ET deposited - Eta vs pt
483 TH2F *fHistPionRecEtaETDep;// ET deposited - Eta
484 TH2F *fHistPionMatchEtaETDep;// ET deposited - Eta
486 TH2F *fHistPionRecResEET;// ET - track matching residual vs E
487 TH2F *fHistPionRecResPtET;// ET - track matching residual vs pt
488 TH2F *fHistPionRecResE;// multiplicity - track matching residual vs E
489 TH2F *fHistPionRecResPt;// multiplicity - track matching residual vs pt
490 TH2F *fHistPionRecResEDepETDep;// ET deposited - track matching residual vs E deposited
491 TH2F *fHistPionRecResPtETDep;// ET deposited - track matching residual vs pt
493 // *******************
494 // charged kaon (+ and -) ET
495 // *******************
496 TH2F *fHistKaonEtaEET;//ET - Eta vs E
497 TH2F *fHistKaonAccEtaEET;//ET - Eta vs E
498 TH2F *fHistKaonRecEtaEET;//ET - Eta vs E
499 TH2F *fHistKaonMatchEtaEET;//ET - Eta vs E
501 TH2F *fHistKaonEtaPtET;// ET - Eta vs pt
502 TH2F *fHistKaonAccEtaPtET;// ET - Eta vs pt
503 TH2F *fHistKaonRecEtaPtET;// ET - Eta vs pt
504 TH2F *fHistKaonMatchEtaPtET;// ET - Eta vs pt
506 TH2F *fHistKaonEtaET;// ET - Eta
507 TH2F *fHistKaonAccEtaET;// ET - Eta
508 TH2F *fHistKaonRecEtaET;// ET - Eta
509 TH2F *fHistKaonMatchEtaET;// ET - Eta
511 TH2F *fHistKaonEtaE;// multiplicity - Eta vs E
512 TH2F *fHistKaonAccEtaE;// multiplicity - Eta vs E
513 TH2F *fHistKaonRecEtaE;// multiplicity - Eta vs E
514 TH2F *fHistKaonMatchEtaE;// multiplicity - Eta vs E
516 TH2F *fHistKaonEtaPt;// multiplicity - Eta vs pt
517 TH2F *fHistKaonAccEtaPt;// multiplicity - Eta vs pt
518 TH2F *fHistKaonRecEtaPt;// multiplicity - Eta vs pt
519 TH2F *fHistKaonMatchEtaPt;// multiplicity - Eta vs pt
521 TH1F *fHistKaontotET;// total ET distribution
522 TH1F *fHistKaonAcctotET;// total ET distribution
523 TH1F *fHistKaonRectotET;// total ET distribution
524 TH1F *fHistKaonMatchtotET;// total ET distribution
526 TH1F *fHistKaonRectotETDep;// total deposited ET distribution
527 TH1F *fHistKaonMatchtotETDep;// total deposited ET distribution
529 TH2F *fHistKaonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
530 TH2F *fHistKaonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
532 TH2F *fHistKaonRecEtaPtETDep;// ET deposited - Eta vs pt
533 TH2F *fHistKaonMatchEtaPtETDep;// ET deposited - Eta vs pt
535 TH2F *fHistKaonRecEtaETDep;// ET deposited - Eta
536 TH2F *fHistKaonMatchEtaETDep;// ET deposited - Eta
538 TH2F *fHistKaonRecResEET;// ET - track matching residual vs E
539 TH2F *fHistKaonRecResPtET;// ET - track matching residual vs pt
540 TH2F *fHistKaonRecResE;// multiplicity - track matching residual vs E
541 TH2F *fHistKaonRecResPt;// multiplicity - track matching residual vs pt
542 TH2F *fHistKaonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
543 TH2F *fHistKaonRecResPtETDep;// ET deposited - track matching residual vs pt
545 // *******************
547 // *******************
548 TH2F *fHistProtonEtaEET;//ET - Eta vs E
549 TH2F *fHistProtonAccEtaEET;//ET - Eta vs E
550 TH2F *fHistProtonRecEtaEET;//ET - Eta vs E
551 TH2F *fHistProtonMatchEtaEET;//ET - Eta vs E
553 TH2F *fHistProtonEtaPtET;// ET - Eta vs pt
554 TH2F *fHistProtonAccEtaPtET;// ET - Eta vs pt
555 TH2F *fHistProtonRecEtaPtET;// ET - Eta vs pt
556 TH2F *fHistProtonMatchEtaPtET;// ET - Eta vs pt
558 TH2F *fHistProtonEtaET;// ET - Eta
559 TH2F *fHistProtonAccEtaET;// ET - Eta
560 TH2F *fHistProtonRecEtaET;// ET - Eta
561 TH2F *fHistProtonMatchEtaET;// ET - Eta
563 TH2F *fHistProtonEtaE;// multiplicity - Eta vs E
564 TH2F *fHistProtonAccEtaE;// multiplicity - Eta vs E
565 TH2F *fHistProtonRecEtaE;// multiplicity - Eta vs E
566 TH2F *fHistProtonMatchEtaE;// multiplicity - Eta vs E
568 TH2F *fHistProtonEtaPt;// multiplicity - Eta vs pt
569 TH2F *fHistProtonAccEtaPt;// multiplicity - Eta vs pt
570 TH2F *fHistProtonRecEtaPt;// multiplicity - Eta vs pt
571 TH2F *fHistProtonMatchEtaPt;// multiplicity - Eta vs pt
573 TH1F *fHistProtontotET;// total ET distribution
574 TH1F *fHistProtonAcctotET;// total ET distribution
575 TH1F *fHistProtonRectotET;// total ET distribution
576 TH1F *fHistProtonMatchtotET;// total ET distribution
578 TH1F *fHistProtonRectotETDep;// total deposited ET distribution
579 TH1F *fHistProtonMatchtotETDep;// total deposited ET distribution
581 TH2F *fHistProtonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
582 TH2F *fHistProtonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
584 TH2F *fHistProtonRecEtaPtETDep;// ET deposited - Eta vs pt
585 TH2F *fHistProtonMatchEtaPtETDep;// ET deposited - Eta vs pt
587 TH2F *fHistProtonRecEtaETDep;// ET deposited - Eta
588 TH2F *fHistProtonMatchEtaETDep;// ET deposited - Eta
590 TH2F *fHistProtonRecResEET;// ET - track matching residual vs E
591 TH2F *fHistProtonRecResPtET;// ET - track matching residual vs pt
592 TH2F *fHistProtonRecResE;// multiplicity - track matching residual vs E
593 TH2F *fHistProtonRecResPt;// multiplicity - track matching residual vs pt
594 TH2F *fHistProtonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
595 TH2F *fHistProtonRecResPtETDep;// ET deposited - track matching residual vs pt
597 // *******************
599 // *******************
600 TH1F *fHistTotChargedtotET;//total ET distribution
601 TH1F *fHistTotChargedAcctotET;//total ET distribution
602 TH1F *fHistTotChargedRectotET;//total ET distribution
603 TH1F *fHistTotChargedRectotETDep;//total deposited ET distribution
604 TH1F *fHistTotChargedMatchtotET;//total ET distribution
605 TH1F *fHistTotChargedMatchtotETDep;//total deposited ET distribution
607 // *******************
609 // *******************
610 TH2F *fHistNeutronEtaEET;//ET - Eta vs E
611 TH2F *fHistNeutronAccEtaEET;//ET - Eta vs E
612 TH2F *fHistNeutronRecEtaEET;//ET - Eta vs E
614 TH2F *fHistNeutronEtaPtET;// ET - Eta vs pt
615 TH2F *fHistNeutronAccEtaPtET;// ET - Eta vs pt
616 TH2F *fHistNeutronRecEtaPtET;// ET - Eta vs pt
618 TH2F *fHistNeutronEtaET;// ET - Eta
619 TH2F *fHistNeutronAccEtaET;// ET - Eta
620 TH2F *fHistNeutronRecEtaET;// ET - Eta
622 TH2F *fHistNeutronEtaE;// multiplicity - Eta vs E
623 TH2F *fHistNeutronAccEtaE;// multiplicity - Eta vs E
624 TH2F *fHistNeutronRecEtaE;// multiplicity - Eta vs E
626 TH2F *fHistNeutronEtaPt;// multiplicity - Eta vs pt
627 TH2F *fHistNeutronAccEtaPt;// multiplicity - Eta vs pt
628 TH2F *fHistNeutronRecEtaPt;// multiplicity - Eta vs pt
630 TH1F *fHistNeutrontotET;// total ET distribution
631 TH1F *fHistNeutronAcctotET;// total ET distribution
632 TH1F *fHistNeutronRectotET;// total ET distribution
633 TH1F *fHistNeutronRectotETDep;// total deposited ET distribution
635 TH2F *fHistNeutronRecEtaEDepETDep;// ET deposited - Eta vs E deposited
636 TH2F *fHistNeutronRecEtaETDep;// ET deposited - Eta
638 TH2F *fHistNeutronRecEtaPtETDep;// ET deposited - Eta vs pt
640 // *******************
642 // *******************
643 TH2F *fHistK0EtaEET;//ET - Eta vs E
644 TH2F *fHistK0RecEtaEET;//ET - Eta vs E
646 TH2F *fHistK0EtaPtET;// ET - Eta vs pt
647 TH2F *fHistK0RecEtaPtET;// ET - Eta vs pt
649 TH2F *fHistK0EtaET;// ET - Eta
650 TH2F *fHistK0RecEtaET;// ET - Eta
652 TH2F *fHistK0EtaE;// multiplicity - Eta vs E
653 TH2F *fHistK0RecEtaE;// multiplicity - Eta vs E
655 TH2F *fHistK0EtaPt;// multiplicity - Eta vs pt
656 TH2F *fHistK0RecEtaPt;// multiplicity - Eta vs pt
658 TH1F *fHistK0totET;// total ET distribution
659 TH1F *fHistK0RectotET;// total ET distribution
661 TH1F *fHistK0RectotETDep;// total deposited ET distribution
663 TH2F *fHistK0RecEtaEDepETDep;// ET deposited - Eta vs E deposited
664 TH2F *fHistK0RecEtaETDep;// ET deposited - Eta
666 TH2F *fHistK0RecEtaPtETDep;// ET deposited - Eta vs pt
668 // *******************
670 // *******************
671 TH2F *fHistLambdaEtaEET;//ET - Eta vs E
672 TH2F *fHistLambdaRecEtaEET;//ET - Eta vs E
674 TH2F *fHistLambdaEtaPtET;// ET - Eta vs pt
675 TH2F *fHistLambdaRecEtaPtET;// ET - Eta vs pt
677 TH2F *fHistLambdaEtaET;// ET - Eta
678 TH2F *fHistLambdaRecEtaET;// ET - Eta
680 TH2F *fHistLambdaEtaE;// multiplicity - Eta vs E
681 TH2F *fHistLambdaRecEtaE;// multiplicity - Eta vs E
683 TH2F *fHistLambdaEtaPt;// multiplicity - Eta vs pt
684 TH2F *fHistLambdaRecEtaPt;// multiplicity - Eta vs pt
686 TH1F *fHistLambdatotET;// total ET distribution
687 TH1F *fHistLambdaRectotET;// total ET distribution
689 TH1F *fHistLambdaRectotETDep;// total deposited ET distribution
691 TH2F *fHistLambdaRecEtaEDepETDep;// ET deposited - Eta vs E deposited
692 TH2F *fHistLambdaRecEtaETDep;// ET deposited - Eta
694 TH2F *fHistLambdaRecEtaPtETDep;// ET deposited - Eta vs pt
696 // *******************
698 // *******************
699 TH1F *fHistTotNeutraltotET;//total ET distribution
700 TH1F *fHistTotNeutralRectotET;//total ET distribution
701 TH1F *fHistTotNeutralRectotETDep;//total deposited ET distribution
703 // *******************
705 // *******************
706 TH1F *fHistTotaltotET;//total ET distribution
707 TH1F *fHistTotalAcctotET;//total ET distribution
708 TH1F *fHistTotalRectotET;//total ET distribution
709 TH1F *fHistTotalRectotETDep;//total deposited ET distribution
711 // *******************
713 // *******************
715 // check produced electrons
716 TH1F *fHistElectronFirstMother;// first mother ID
717 TH2F *fHistElectronFirstMotherXY;// first mother XY position
718 TH1F *fHistElectronNDaughters;// number of daughters
719 TH1F *fHistElectronDaughters;// daughters ID
720 TH2F *fHistElectronDaughtersXY;// daughters XY position
722 TH1F *fHistElectronFirstMotherAcc;// first mother ID
723 TH2F *fHistElectronFirstMotherXYAcc;// first mother XY position
724 TH1F *fHistElectronNDaughtersAcc;// number of daughters
725 TH1F *fHistElectronDaughtersAcc;// daughters ID
726 TH2F *fHistElectronDaughtersXYAcc;// daughters XY position
728 TH1F *fHistElectronFirstMotherRec;// first mother ID
729 TH2F *fHistElectronFirstMotherXYRec;// first mother XY position
730 TH1F *fHistElectronNDaughtersRec;// number of daughters
731 TH1F *fHistElectronDaughtersRec;// daughters ID
732 TH2F *fHistElectronDaughtersXYRec;// daughters XY position
734 TH1F *fHistNPPElectronFirstMother;// first mother ID
735 TH2F *fHistNPPElectronFirstMotherXY;// first mother XY position
736 TH1F *fHistNPPElectronNDaughters;// number of daughters
737 TH1F *fHistNPPElectronDaughters;// daughters ID
738 TH2F *fHistNPPElectronDaughtersXY;// daughters XY position
740 TH1F *fHistNPPElectronFirstMotherAcc;// first mother ID
741 TH2F *fHistNPPElectronFirstMotherXYAcc;// first mother XY position
742 TH1F *fHistNPPElectronNDaughtersAcc;// number of daughters
743 TH1F *fHistNPPElectronDaughtersAcc;// daughters ID
744 TH2F *fHistNPPElectronDaughtersXYAcc;// daughters XY position
746 TH1F *fHistNPPElectronFirstMotherRec;// first mother ID
747 TH2F *fHistNPPElectronFirstMotherXYRec;// first mother XY position
748 TH1F *fHistNPPElectronNDaughtersRec;// number of daughters
749 TH1F *fHistNPPElectronDaughtersRec;// daughters ID
750 TH2F *fHistNPPElectronDaughtersXYRec;// daughters XY position
752 // check produced gammas
753 TH1F *fHistGammaFirstMother;// first mother ID
754 TH2F *fHistGammaFirstMotherXY;// first mother XY position
755 TH1F *fHistGammaNDaughters;// number of daughters
756 TH1F *fHistGammaDaughters;// daughters ID
757 TH2F *fHistGammaDaughtersXY;// daughters XY position
758 TH2F *fHistConvGammaDaughtersXY;// daughters XY position
759 TH2F *fHistNonConvGammaDaughtersXY;// daughters XY position
761 TH1F *fHistGammaFirstMotherAcc;// first mother ID
762 TH2F *fHistGammaFirstMotherXYAcc;// first mother XY position
763 TH1F *fHistGammaNDaughtersAcc;// number of daughters
764 TH1F *fHistGammaDaughtersAcc;// daughters ID
765 TH2F *fHistGammaDaughtersXYAcc;// daughters XY position
766 TH2F *fHistConvGammaDaughtersXYAcc;// daughters XY position
767 TH2F *fHistNonConvGammaDaughtersXYAcc;// daughters XY position
769 TH1F *fHistGammaFirstMotherRec;// first mother ID
770 TH2F *fHistGammaFirstMotherXYRec;// first mother XY position
771 TH1F *fHistGammaNDaughtersRec;// number of daughters
772 TH1F *fHistGammaDaughtersRec;// daughters ID
773 TH2F *fHistGammaDaughtersXYRec;// daughters XY position
774 TH2F *fHistConvGammaDaughtersXYRec;// daughters XY position
775 TH2F *fHistNonConvGammaDaughtersXYRec;// daughters XY position
777 TH1F *fHistNPPGammaFirstMother;// first mother ID
778 TH2F *fHistNPPGammaFirstMotherXY;// first mother XY position
779 TH1F *fHistNPPGammaNDaughters;// number of daughters
780 TH1F *fHistNPPGammaDaughters;// daughters ID
781 TH2F *fHistNPPGammaDaughtersXY;// daughters XY position
783 TH1F *fHistNPPGammaFirstMotherAcc;// first mother ID
784 TH2F *fHistNPPGammaFirstMotherXYAcc;// first mother XY position
785 TH1F *fHistNPPGammaNDaughtersAcc;// number of daughters
786 TH1F *fHistNPPGammaDaughtersAcc;// daughters ID
787 TH2F *fHistNPPGammaDaughtersXYAcc;// daughters XY position
789 TH1F *fHistNPPGammaFirstMotherRec;// first mother ID
790 TH2F *fHistNPPGammaFirstMotherXYRec;// first mother XY position
791 TH1F *fHistNPPGammaNDaughtersRec;// number of daughters
792 TH1F *fHistNPPGammaDaughtersRec;// daughters ID
793 TH2F *fHistNPPGammaDaughtersXYRec;// daughters XY position
796 TH2F *fHistAllERecEMC;// E reconstructed vs E MC
797 TH2F *fHistAllPtRecPtMC;// pt reconstructed vs pt MC
798 TH2F *fHistElectronERecEMC;// E reconstructed vs E MC
799 TH2F *fHistGammaERecEMC;// E reconstructed vs E MC
801 TH2F *fHistChargedRes;// charged particle track matching residual
802 TH2F *fHistChargedRes2;// charged particle track matching residual
803 TH2F *fHistChargedRes3;// charged particle track matching residual
804 TH2F *fHistNeutralRes;// neutral particle track matching residual
805 TH2F *fHistElectronRes;// electron track matching residual
806 TH2F *fHistGammaRes;// gamma track matching residual
808 TH2F *fHistIsInAcc;// EMCal acceptance check
812 //Declare it private to avoid compilation warning
813 AliAnalysisEmEtMonteCarlo & operator = (const AliAnalysisEmEtMonteCarlo & g) ;//cpy assignment
814 AliAnalysisEmEtMonteCarlo(const AliAnalysisEmEtMonteCarlo & g) ; // cpy ctor
815 ClassDef(AliAnalysisEmEtMonteCarlo, 1);
818 #endif //ALIANALYSISEMETMONTECARLO_H