]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/totEt/AliAnalysisEmEtMonteCarlo.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Making corrections second order in cluster multiplicity
[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEmEtMonteCarlo.cxx
CommitLineData
2e12941b 1//_________________________________________________________________________
2// Utility Class for transverse energy studies
3// Base class for MC analysis
4// - MC output
5// implementation file
6//
7//*-- Author: Marcelo G. Munhoz (USP)
112b1d0c 8//_________________________________________________________________________
2e12941b 9
10#include "AliAnalysisEmEtMonteCarlo.h"
11#include "AliAnalysisEtCuts.h"
112b1d0c 12#include "AliAnalysisEtSelectorEmcal.h"
2e12941b 13#include "AliESDtrack.h"
14#include "AliStack.h"
15#include "AliVEvent.h"
16#include "AliMCEvent.h"
17#include "AliMCParticle.h"
18#include "AliTrackReference.h"
19#include "AliESDEvent.h"
20#include "TH2F.h"
21#include "TParticle.h"
22#include "AliGenHijingEventHeader.h"
23#include "AliGenPythiaEventHeader.h"
24#include "TList.h"
25#include "AliESDCaloCluster.h"
26#include "TGeoGlobalMagField.h"
27#include "AliMagF.h"
28#include "AliEMCALTrack.h"
29#include "AliESDtrackCuts.h"
30#include "AliEMCALGeometry.h"
31#include "AliExternalTrackParam.h"
32#include "AliTrackerBase.h"
33#include "TGeoManager.h"
2c70dc4d 34#include "AliCentrality.h"
2e12941b 35
112b1d0c 36using namespace std;
2e12941b 37
38ClassImp(AliAnalysisEmEtMonteCarlo);
39
40
41// ctor
42AliAnalysisEmEtMonteCarlo::AliAnalysisEmEtMonteCarlo():AliAnalysisEtMonteCarlo()
112b1d0c 43,fResCut(0),fPrimtotET(0), fPrimAcctotET(0), fPrimRectotET(0), fPrimRectotETDep(0)
44,fElectrontotET(0), fElectronAcctotET(0), fElectronRectotET(0)
45,fConvElectrontotET(0), fConvElectronAcctotET(0), fConvElectronRectotET(0), fScatElectrontotET(0), fScatElectronAcctotET(0), fScatElectronRectotET(0)
46,fTotElectrontotET(0), fTotElectronAcctotET(0), fTotElectronRectotET(0)
47
48,fGammatotET(0), fGammaAcctotET(0), fGammaRectotET(0)
49,fAnnihGammatotET(0), fAnnihGammaAcctotET(0), fAnnihGammaRectotET(0), fScatGammatotET(0), fScatGammaAcctotET(0), fScatGammaRectotET(0)
50,fTotGammatotET(0), fTotGammaAcctotET(0), fTotGammaRectotET(0)
51,fConvGammatotET(0),fNonConvGammatotET(0),fConvGammaAcctotET(0),fNonConvGammaAcctotET(0), fNPPPi0GammatotET(0), fNPPPi0GammaRectotET(0)
52
53,fTotEMtotET(0), fTotEMAcctotET(0), fTotEMRectotET(0)
54
55,fNPPElectrontotET(0), fNPPElectronRectotET(0), fNPPGammatotET(0), fNPPGammaRectotET(0)
56,fTotNPPEMtotET(0), fTotNPPEMRectotET(0)
57
58,fMuontotET(0), fPiontotET(0), fKaontotET(0), fProtontotET(0)
59,fMuonAcctotET(0), fPionAcctotET(0), fKaonAcctotET(0), fProtonAcctotET(0)
60,fMuonRectotET(0), fMuonRectotETDep(0), fPionRectotET(0), fPionRectotETDep(0), fKaonRectotET(0), fKaonRectotETDep(0), fProtonRectotET(0), fProtonRectotETDep(0)
61,fMuonMatchtotET(0), fMuonMatchtotETDep(0), fPionMatchtotET(0), fPionMatchtotETDep(0), fKaonMatchtotET(0), fKaonMatchtotETDep(0), fProtonMatchtotET(0), fProtonMatchtotETDep(0)
62,fTotChargedtotET(0), fTotChargedAcctotET(0), fTotChargedRectotET(0), fTotChargedRectotETDep(0), fTotChargedMatchtotET(0), fTotChargedMatchtotETDep(0)
63
64,fNeutrontotET(0), fNeutronAcctotET(0), fNeutronRectotET(0), fNeutronRectotETDep(0)
65,fK0totET(0), fK0RectotET(0), fK0RectotETDep(0), fLambdatotET(0), fLambdaRectotET(0), fLambdaRectotETDep(0)
66,fTotNeutraltotET(0), fTotNeutralRectotET(0), fTotNeutralRectotETDep(0)
67
68,fTotaltotET(0), fTotalAcctotET(0), fTotalRectotET(0), fTotalRectotETDep(0)
69
70,fGeoUt(0)
71
72,fHistPrimEtaEET(0)
73,fHistPrimEtaPtET(0)
74,fHistPrimEtaET(0)
75,fHistPrimtotET(0)
76
77,fHistPrimAccEtaEET(0)
78,fHistPrimAccEtaPtET(0)
79,fHistPrimAccEtaET(0)
80,fHistPrimAcctotET(0)
81
82,fHistPrimRecEtaEET(0)
83,fHistPrimRecEtaPtET(0)
84,fHistPrimRecEtaET(0)
85,fHistPrimRectotET(0)
86
87,fHistPrimRecEtaEDepETDep(0)
88,fHistPrimRecEtaPtETDep(0)
89,fHistPrimRecEtaETDep(0)
90,fHistPrimRectotETDep(0)
91
92,fHistElectronEtaEET(0)
93,fHistElectronEtaPtET(0)
94,fHistElectronEtaET(0)
95,fHistElectronEtaE(0)
96,fHistElectronEtaPt(0)
97,fHistElectrontotET(0)
98
99,fHistConvElectronEtaEET(0)
100,fHistConvElectronEtaPtET(0)
101,fHistConvElectronEtaET(0)
102,fHistConvElectronEtaE(0)
103,fHistConvElectronEtaPt(0)
104,fHistConvElectrontotET(0)
105
106,fHistScatElectronEtaEET(0)
107,fHistScatElectronEtaPtET(0)
108,fHistScatElectronEtaET(0)
109,fHistScatElectronEtaE(0)
110,fHistScatElectronEtaPt(0)
111,fHistScatElectrontotET(0)
112
113,fHistTotElectrontotET(0)
114
115,fHistGammaEtaEET(0)
116,fHistGammaEtaPtET(0)
117,fHistGammaEtaET(0)
118,fHistGammaEtaE(0)
119,fHistGammaEtaPt(0)
120,fHistGammatotET(0)
121
122,fHistAnnihGammaEtaEET(0)
123,fHistAnnihGammaEtaPtET(0)
124,fHistAnnihGammaEtaET(0)
125,fHistAnnihGammaEtaE(0)
126,fHistAnnihGammaEtaPt(0)
127,fHistAnnihGammatotET(0)
128
129,fHistScatGammaEtaEET(0)
130,fHistScatGammaEtaPtET(0)
131,fHistScatGammaEtaET(0)
132,fHistScatGammaEtaE(0)
133,fHistScatGammaEtaPt(0)
134,fHistScatGammatotET(0)
135
136,fHistConvGammaEtaEET(0)
137,fHistConvGammaEtaPtET(0)
138,fHistConvGammaEtaET(0)
139,fHistConvGammaEtaE(0)
140,fHistConvGammaEtaPt(0)
141,fHistConvGammatotET(0)
142
143,fHistNonConvGammaEtaEET(0)
144,fHistNonConvGammaEtaPtET(0)
145,fHistNonConvGammaEtaET(0)
146,fHistNonConvGammaEtaE(0)
147,fHistNonConvGammaEtaPt(0)
148,fHistNonConvGammatotET(0)
149
150,fHistTotGammatotET(0)
151
152,fHistTotEMtotET(0)
153
154,fHistNPPElectronEtaEET(0)
155,fHistNPPElectronEtaPtET(0)
156,fHistNPPElectronEtaET(0)
157,fHistNPPElectronEtaE(0)
158,fHistNPPElectronEtaPt(0)
159,fHistNPPElectrontotET(0)
160
161,fHistNPPGammaEtaEET(0)
162,fHistNPPGammaEtaPtET(0)
163,fHistNPPGammaEtaET(0)
164,fHistNPPGammaEtaE(0)
165,fHistNPPGammaEtaPt(0)
166,fHistNPPGammatotET(0)
167
168,fHistTotNPPEMtotET(0)
169
170,fHistNPPPi0GammaEtaEET(0)
171,fHistNPPPi0GammaEtaPtET(0)
172,fHistNPPPi0GammaEtaET(0)
173,fHistNPPPi0GammaEtaE(0)
174,fHistNPPPi0GammaEtaPt(0)
175,fHistNPPPi0GammatotET(0)
176
177,fHistElectronAccEtaEET(0)
178,fHistElectronAccEtaPtET(0)
179,fHistElectronAccEtaET(0)
180,fHistElectronAccEtaE(0)
181,fHistElectronAccEtaPt(0)
182,fHistElectronAcctotET(0)
183
184,fHistConvElectronAccEtaEET(0)
185,fHistConvElectronAccEtaPtET(0)
186,fHistConvElectronAccEtaET(0)
187,fHistConvElectronAccEtaE(0)
188,fHistConvElectronAccEtaPt(0)
189,fHistConvElectronAcctotET(0)
190
191,fHistScatElectronAccEtaEET(0)
192,fHistScatElectronAccEtaPtET(0)
193,fHistScatElectronAccEtaET(0)
194,fHistScatElectronAccEtaE(0)
195,fHistScatElectronAccEtaPt(0)
196,fHistScatElectronAcctotET(0)
197
198,fHistTotElectronAcctotET(0)
199
200,fHistGammaAccEtaEET(0)
201,fHistGammaAccEtaPtET(0)
202,fHistGammaAccEtaET(0)
203,fHistGammaAccEtaE(0)
204,fHistGammaAccEtaPt(0)
205,fHistGammaAcctotET(0)
206
207,fHistAnnihGammaAccEtaEET(0)
208,fHistAnnihGammaAccEtaPtET(0)
209,fHistAnnihGammaAccEtaET(0)
210,fHistAnnihGammaAccEtaE(0)
211,fHistAnnihGammaAccEtaPt(0)
212,fHistAnnihGammaAcctotET(0)
213
214,fHistScatGammaAccEtaEET(0)
215,fHistScatGammaAccEtaPtET(0)
216,fHistScatGammaAccEtaET(0)
217,fHistScatGammaAccEtaE(0)
218,fHistScatGammaAccEtaPt(0)
219,fHistScatGammaAcctotET(0)
220
221,fHistConvGammaAccEtaEET(0)
222,fHistConvGammaAccEtaPtET(0)
223,fHistConvGammaAccEtaET(0)
224,fHistConvGammaAccEtaE(0)
225,fHistConvGammaAccEtaPt(0)
226,fHistConvGammaAcctotET(0)
227
228,fHistNonConvGammaAccEtaEET(0)
229,fHistNonConvGammaAccEtaPtET(0)
230,fHistNonConvGammaAccEtaET(0)
231,fHistNonConvGammaAccEtaE(0)
232,fHistNonConvGammaAccEtaPt(0)
233,fHistNonConvGammaAcctotET(0)
234
235,fHistTotGammaAcctotET(0)
236
237,fHistTotEMAcctotET(0)
238
239,fHistNPPElectronAccEtaEET(0)
240,fHistNPPElectronAccEtaPtET(0)
241,fHistNPPElectronAccEtaE(0)
242,fHistNPPElectronAccEtaPt(0)
243
244,fHistNPPGammaAccEtaEET(0)
245,fHistNPPGammaAccEtaPtET(0)
246,fHistNPPGammaAccEtaE(0)
247,fHistNPPGammaAccEtaPt(0)
248
249,fHistElectronRecEtaEET(0)
250,fHistElectronRecEtaPtET(0)
251,fHistElectronRecEtaET(0)
252,fHistElectronRecEtaE(0)
253,fHistElectronRecEtaPt(0)
254,fHistElectronRectotET(0)
255
256,fHistConvElectronRecEtaEET(0)
257,fHistConvElectronRecEtaPtET(0)
258,fHistConvElectronRecEtaET(0)
259,fHistConvElectronRecEtaE(0)
260,fHistConvElectronRecEtaPt(0)
261,fHistConvElectronRectotET(0)
262
263,fHistScatElectronRecEtaEET(0)
264,fHistScatElectronRecEtaPtET(0)
265,fHistScatElectronRecEtaET(0)
266,fHistScatElectronRecEtaE(0)
267,fHistScatElectronRecEtaPt(0)
268,fHistScatElectronRectotET(0)
269
270,fHistTotElectronRectotET(0)
271
272,fHistGammaRecEtaEET(0)
273,fHistGammaRecEtaPtET(0)
274,fHistGammaRecEtaET(0)
275,fHistGammaRecEtaE(0)
276,fHistGammaRecEtaPt(0)
277,fHistGammaRectotET(0)
278
279,fHistGammaRecResEET(0) //Added Caio
280,fHistGammaRecResPtET(0)//Added Caio
281,fHistGammaRecResE(0) // Added Caio
282,fHistGammaRecResPt(0) //Added Caio
283,fHistGammaRecResEDepETDep(0) //Added Caio
284,fHistGammaRecResPtETDep(0) //Added Caio
285
286,fHistAnnihGammaRecEtaEET(0)
287,fHistAnnihGammaRecEtaPtET(0)
288,fHistAnnihGammaRecEtaET(0)
289,fHistAnnihGammaRecEtaE(0)
290,fHistAnnihGammaRecEtaPt(0)
291,fHistAnnihGammaRectotET(0)
292
293,fHistScatGammaRecEtaEET(0)
294,fHistScatGammaRecEtaPtET(0)
295,fHistScatGammaRecEtaET(0)
296,fHistScatGammaRecEtaE(0)
297,fHistScatGammaRecEtaPt(0)
298,fHistScatGammaRectotET(0)
299
300,fHistTotGammaRectotET(0)
301
302,fHistTotEMRectotET(0)
303
304,fHistNPPElectronRecEtaEET(0)
305,fHistNPPElectronRecEtaPtET(0)
306,fHistNPPElectronRecEtaET(0)
307,fHistNPPElectronRecEtaE(0)
308,fHistNPPElectronRecEtaPt(0)
309,fHistNPPElectronRectotET(0)
310
311,fHistNPPGammaRecEtaEET(0)
312,fHistNPPGammaRecEtaPtET(0)
313,fHistNPPGammaRecEtaET(0)
314,fHistNPPGammaRecEtaE(0)
315,fHistNPPGammaRecEtaPt(0)
316,fHistNPPGammaRectotET(0)
317
318,fHistTotNPPEMRectotET(0)
319
320,fHistNPPPi0GammaRecEtaEET(0)
321,fHistNPPPi0GammaRecEtaPtET(0)
322,fHistNPPPi0GammaRecEtaET(0)
323,fHistNPPPi0GammaRecEtaE(0)
324,fHistNPPPi0GammaRecEtaPt(0)
325,fHistNPPPi0GammaRectotET(0)
326
327,fHistMuonEtaEET(0)
328,fHistMuonAccEtaEET(0)
329,fHistMuonRecEtaEET(0)
330,fHistMuonMatchEtaEET(0)
331
332,fHistMuonEtaPtET(0)
333,fHistMuonAccEtaPtET(0)
334,fHistMuonRecEtaPtET(0)
335,fHistMuonMatchEtaPtET(0)
336
337,fHistMuonEtaET(0)
338,fHistMuonAccEtaET(0)
339,fHistMuonRecEtaET(0)
340,fHistMuonMatchEtaET(0)
341
342,fHistMuonEtaE(0)
343,fHistMuonAccEtaE(0)
344,fHistMuonRecEtaE(0)
345,fHistMuonMatchEtaE(0)
346
347,fHistMuonEtaPt(0)
348,fHistMuonAccEtaPt(0)
349,fHistMuonRecEtaPt(0)
350,fHistMuonMatchEtaPt(0)
351
352,fHistMuontotET(0)
353,fHistMuonAcctotET(0)
354,fHistMuonRectotET(0)
355,fHistMuonMatchtotET(0)
356
357,fHistMuonRectotETDep(0)
358,fHistMuonMatchtotETDep(0)
359
360,fHistMuonRecEtaEDepETDep(0)
361,fHistMuonMatchEtaEDepETDep(0)
362
363,fHistMuonRecEtaPtETDep(0)
364,fHistMuonMatchEtaPtETDep(0)
365
366,fHistMuonRecEtaETDep(0)
367,fHistMuonMatchEtaETDep(0)
368
369,fHistMuonRecResEET(0)
370,fHistMuonRecResPtET(0)
371,fHistMuonRecResE(0)
372,fHistMuonRecResPt(0)
373
374,fHistMuonRecResEDepETDep(0)
375,fHistMuonRecResPtETDep(0)
376
377,fHistPionEtaEET(0)
378,fHistPionAccEtaEET(0)
379,fHistPionRecEtaEET(0)
380,fHistPionMatchEtaEET(0)
381
382,fHistPionEtaPtET(0)
383,fHistPionAccEtaPtET(0)
384,fHistPionRecEtaPtET(0)
385,fHistPionMatchEtaPtET(0)
386
387,fHistPionEtaET(0)
388,fHistPionAccEtaET(0)
389,fHistPionRecEtaET(0)
390,fHistPionMatchEtaET(0)
391
392,fHistPionEtaE(0)
393,fHistPionAccEtaE(0)
394,fHistPionRecEtaE(0)
395,fHistPionMatchEtaE(0)
396
397,fHistPionEtaPt(0)
398,fHistPionAccEtaPt(0)
399,fHistPionRecEtaPt(0)
400,fHistPionMatchEtaPt(0)
401
402,fHistPiontotET(0)
403,fHistPionAcctotET(0)
404,fHistPionRectotET(0)
405,fHistPionMatchtotET(0)
406
407,fHistPionRectotETDep(0)
408,fHistPionMatchtotETDep(0)
409
410,fHistPionRecEtaEDepETDep(0)
411,fHistPionMatchEtaEDepETDep(0)
412
413,fHistPionRecEtaPtETDep(0)
414,fHistPionMatchEtaPtETDep(0)
415
416,fHistPionRecEtaETDep(0)
417,fHistPionMatchEtaETDep(0)
418
419,fHistPionRecResEET(0)
420,fHistPionRecResPtET(0)
421,fHistPionRecResE(0)
422,fHistPionRecResPt(0)
423,fHistPionRecResEDepETDep(0)
424,fHistPionRecResPtETDep(0)
425
426,fHistKaonEtaEET(0)
427,fHistKaonAccEtaEET(0)
428,fHistKaonRecEtaEET(0)
429,fHistKaonMatchEtaEET(0)
430
431,fHistKaonEtaPtET(0)
432,fHistKaonAccEtaPtET(0)
433,fHistKaonRecEtaPtET(0)
434,fHistKaonMatchEtaPtET(0)
435
436,fHistKaonEtaET(0)
437,fHistKaonAccEtaET(0)
438,fHistKaonRecEtaET(0)
439,fHistKaonMatchEtaET(0)
440
441,fHistKaonEtaE(0)
442,fHistKaonAccEtaE(0)
443,fHistKaonRecEtaE(0)
444,fHistKaonMatchEtaE(0)
445
446,fHistKaonEtaPt(0)
447,fHistKaonAccEtaPt(0)
448,fHistKaonRecEtaPt(0)
449,fHistKaonMatchEtaPt(0)
450
451,fHistKaontotET(0)
452,fHistKaonAcctotET(0)
453,fHistKaonRectotET(0)
454,fHistKaonMatchtotET(0)
455
456,fHistKaonRectotETDep(0)
457,fHistKaonMatchtotETDep(0)
458
459,fHistKaonRecEtaEDepETDep(0)
460,fHistKaonMatchEtaEDepETDep(0)
461
462,fHistKaonRecEtaPtETDep(0)
463,fHistKaonMatchEtaPtETDep(0)
464
465,fHistKaonRecEtaETDep(0)
466,fHistKaonMatchEtaETDep(0)
467
468,fHistKaonRecResEET(0)
469,fHistKaonRecResPtET(0)
470,fHistKaonRecResE(0)
471,fHistKaonRecResPt(0)
472
473,fHistKaonRecResEDepETDep(0)
474,fHistKaonRecResPtETDep(0)
475
476,fHistProtonEtaEET(0)
477,fHistProtonAccEtaEET(0)
478,fHistProtonRecEtaEET(0)
479,fHistProtonMatchEtaEET(0)
480
481,fHistProtonEtaPtET(0)
482,fHistProtonAccEtaPtET(0)
483,fHistProtonRecEtaPtET(0)
484,fHistProtonMatchEtaPtET(0)
2e12941b 485
112b1d0c 486,fHistProtonEtaET(0)
487,fHistProtonAccEtaET(0)
488,fHistProtonRecEtaET(0)
489,fHistProtonMatchEtaET(0)
2e12941b 490
112b1d0c 491,fHistProtonEtaE(0)
492,fHistProtonAccEtaE(0)
493,fHistProtonRecEtaE(0)
494,fHistProtonMatchEtaE(0)
2e12941b 495
112b1d0c 496,fHistProtonEtaPt(0)
497,fHistProtonAccEtaPt(0)
498,fHistProtonRecEtaPt(0)
499,fHistProtonMatchEtaPt(0)
2e12941b 500
112b1d0c 501,fHistProtontotET(0)
502,fHistProtonAcctotET(0)
503,fHistProtonRectotET(0)
504,fHistProtonMatchtotET(0)
2e12941b 505
112b1d0c 506,fHistProtonRectotETDep(0)
507,fHistProtonMatchtotETDep(0)
2e12941b 508
112b1d0c 509,fHistProtonRecEtaEDepETDep(0)
510,fHistProtonMatchEtaEDepETDep(0)
2e12941b 511
112b1d0c 512,fHistProtonRecEtaPtETDep(0)
513,fHistProtonMatchEtaPtETDep(0)
2e12941b 514
112b1d0c 515,fHistProtonRecEtaETDep(0)
516,fHistProtonMatchEtaETDep(0)
2e12941b 517
112b1d0c 518,fHistProtonRecResEET(0)
519,fHistProtonRecResPtET(0)
520,fHistProtonRecResE(0)
521,fHistProtonRecResPt(0)
2e12941b 522
112b1d0c 523,fHistProtonRecResEDepETDep(0)
524,fHistProtonRecResPtETDep(0)
2e12941b 525
112b1d0c 526,fHistTotChargedtotET(0)
527,fHistTotChargedAcctotET(0)
528,fHistTotChargedRectotET(0)
529,fHistTotChargedRectotETDep(0)
530,fHistTotChargedMatchtotET(0)
531,fHistTotChargedMatchtotETDep(0)
2e12941b 532
112b1d0c 533,fHistNeutronEtaEET(0)
534,fHistNeutronAccEtaEET(0)
535,fHistNeutronRecEtaEET(0)
2e12941b 536
112b1d0c 537,fHistNeutronEtaPtET(0)
538,fHistNeutronAccEtaPtET(0)
539,fHistNeutronRecEtaPtET(0)
2e12941b 540
112b1d0c 541,fHistNeutronEtaET(0)
542,fHistNeutronAccEtaET(0)
543,fHistNeutronRecEtaET(0)
2e12941b 544
112b1d0c 545,fHistNeutronEtaE(0)
546,fHistNeutronAccEtaE(0)
547,fHistNeutronRecEtaE(0)
2e12941b 548
112b1d0c 549,fHistNeutronEtaPt(0)
550,fHistNeutronAccEtaPt(0)
551,fHistNeutronRecEtaPt(0)
2e12941b 552
112b1d0c 553,fHistNeutrontotET(0)
554,fHistNeutronAcctotET(0)
555,fHistNeutronRectotET(0)
2e12941b 556
112b1d0c 557,fHistNeutronRectotETDep(0)
2e12941b 558
112b1d0c 559,fHistNeutronRecEtaEDepETDep(0)
560,fHistNeutronRecEtaETDep(0)
561,fHistNeutronRecEtaPtETDep(0)
562
563,fHistK0EtaEET(0)
564,fHistK0RecEtaEET(0)
565
566,fHistK0EtaPtET(0)
567,fHistK0RecEtaPtET(0)
568
569,fHistK0EtaET(0)
570,fHistK0RecEtaET(0)
571
572,fHistK0EtaE(0)
573,fHistK0RecEtaE(0)
574
575,fHistK0EtaPt(0)
576,fHistK0RecEtaPt(0)
577
578,fHistK0totET(0)
579,fHistK0RectotET(0)
580,fHistK0RectotETDep(0)
581
582,fHistK0RecEtaEDepETDep(0)
583,fHistK0RecEtaETDep(0)
584
585,fHistK0RecEtaPtETDep(0)
586
587,fHistLambdaEtaEET(0)
588,fHistLambdaRecEtaEET(0)
589
590,fHistLambdaEtaPtET(0)
591,fHistLambdaRecEtaPtET(0)
592
593,fHistLambdaEtaET(0)
594,fHistLambdaRecEtaET(0)
595
596,fHistLambdaEtaE(0)
597,fHistLambdaRecEtaE(0)
598
599,fHistLambdaEtaPt(0)
600,fHistLambdaRecEtaPt(0)
601
602,fHistLambdatotET(0)
603,fHistLambdaRectotET(0)
604,fHistLambdaRectotETDep(0)
605
606,fHistLambdaRecEtaEDepETDep(0)
607,fHistLambdaRecEtaETDep(0)
608
609,fHistLambdaRecEtaPtETDep(0)
610
611,fHistTotNeutraltotET(0)
612,fHistTotNeutralRectotET(0)
613,fHistTotNeutralRectotETDep(0)
614
615,fHistTotaltotET(0)
616,fHistTotalAcctotET(0)
617,fHistTotalRectotET(0)
618,fHistTotalRectotETDep(0)
619
620,fHistElectronFirstMother(0)
621,fHistElectronFirstMotherXY(0)
622,fHistElectronNDaughters(0)
623,fHistElectronDaughters(0)
624,fHistElectronDaughtersXY(0)
625
626,fHistElectronFirstMotherAcc(0)
627,fHistElectronFirstMotherXYAcc(0)
628,fHistElectronNDaughtersAcc(0)
629,fHistElectronDaughtersAcc(0)
630,fHistElectronDaughtersXYAcc(0)
631
632,fHistElectronFirstMotherRec(0)
633,fHistElectronFirstMotherXYRec(0)
634,fHistElectronNDaughtersRec(0)
635,fHistElectronDaughtersRec(0)
636,fHistElectronDaughtersXYRec(0)
637
638,fHistNPPElectronFirstMother(0)
639,fHistNPPElectronFirstMotherXY(0)
640,fHistNPPElectronNDaughters(0)
641,fHistNPPElectronDaughters(0)
642,fHistNPPElectronDaughtersXY(0)
643
644,fHistNPPElectronFirstMotherAcc(0)
645,fHistNPPElectronFirstMotherXYAcc(0)
646,fHistNPPElectronNDaughtersAcc(0)
647,fHistNPPElectronDaughtersAcc(0)
648,fHistNPPElectronDaughtersXYAcc(0)
649
650,fHistNPPElectronFirstMotherRec(0)
651,fHistNPPElectronFirstMotherXYRec(0)
652,fHistNPPElectronNDaughtersRec(0)
653,fHistNPPElectronDaughtersRec(0)
654,fHistNPPElectronDaughtersXYRec(0)
655
656,fHistGammaFirstMother(0)
657,fHistGammaFirstMotherXY(0)
658,fHistGammaNDaughters(0)
659,fHistGammaDaughters(0)
660,fHistGammaDaughtersXY(0)
661,fHistConvGammaDaughtersXY(0)
662,fHistNonConvGammaDaughtersXY(0)
663
664,fHistGammaFirstMotherAcc(0)
665,fHistGammaFirstMotherXYAcc(0)
666,fHistGammaNDaughtersAcc(0)
667,fHistGammaDaughtersAcc(0)
668,fHistGammaDaughtersXYAcc(0)
669,fHistConvGammaDaughtersXYAcc(0)
670,fHistNonConvGammaDaughtersXYAcc(0)
671
672,fHistGammaFirstMotherRec(0)
673,fHistGammaFirstMotherXYRec(0)
674,fHistGammaNDaughtersRec(0)
675,fHistGammaDaughtersRec(0)
676,fHistGammaDaughtersXYRec(0)
677,fHistConvGammaDaughtersXYRec(0)
678,fHistNonConvGammaDaughtersXYRec(0)
679
680,fHistNPPGammaFirstMother(0)
681,fHistNPPGammaFirstMotherXY(0)
682,fHistNPPGammaNDaughters(0)
683,fHistNPPGammaDaughters(0)
684,fHistNPPGammaDaughtersXY(0)
685
686,fHistNPPGammaFirstMotherAcc(0)
687,fHistNPPGammaFirstMotherXYAcc(0)
688,fHistNPPGammaNDaughtersAcc(0)
689,fHistNPPGammaDaughtersAcc(0)
690,fHistNPPGammaDaughtersXYAcc(0)
691
692,fHistNPPGammaFirstMotherRec(0)
693,fHistNPPGammaFirstMotherXYRec(0)
694,fHistNPPGammaNDaughtersRec(0)
695,fHistNPPGammaDaughtersRec(0)
696,fHistNPPGammaDaughtersXYRec(0)
697
698,fHistAllERecEMC(0)
699,fHistAllPtRecPtMC(0)
700,fHistElectronERecEMC(0)
701,fHistGammaERecEMC(0)
702
703,fHistChargedRes(0)
704,fHistChargedRes2(0)
705,fHistChargedRes3(0)
706,fHistNeutralRes(0)
707,fHistElectronRes(0)
708,fHistGammaRes(0)
709
710,fHistIsInAcc(0)
2e12941b 711{//constructor
112b1d0c 712 fHistogramNameSuffix = TString("EmcalMC");
2e12941b 713
112b1d0c 714 fResCut = 0.02;
715 //fResCut = fEmcalTrackDistanceCut;
2e12941b 716
112b1d0c 717 TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
718 //TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 1., 1., AliMagF::k5kG));
719 TGeoManager::Import("geometry.root");
720 //fGeoUt = new AliEMCALGeometry("EMCAL_FIRSTYEAR","EMCAL");
2e12941b 721}
722
723// dtor
724AliAnalysisEmEtMonteCarlo::~AliAnalysisEmEtMonteCarlo()
725{//destructor
112b1d0c 726 //Marcelo, have you really cleaned up all memory here? What about the histos?
727 delete fGeoUt;
728 delete fHistPrimEtaEET;
729 delete fHistPrimEtaPtET;
730 delete fHistPrimEtaET;
731 delete fHistPrimtotET;
732
733 delete fHistPrimAccEtaEET;
734 delete fHistPrimAccEtaPtET;
735 delete fHistPrimAccEtaET;
736 delete fHistPrimAcctotET;
737
738 delete fHistPrimRecEtaEET;
739 delete fHistPrimRecEtaPtET;
740 delete fHistPrimRecEtaET;
741 delete fHistPrimRectotET;
742
743 delete fHistPrimRecEtaEDepETDep;
744 delete fHistPrimRecEtaPtETDep;
745 delete fHistPrimRecEtaETDep;
746 delete fHistPrimRectotETDep;
747
748 // *******************
749 // electron ET
750 // *******************
751 delete fHistElectronEtaEET;
752 delete fHistElectronEtaPtET;
753 delete fHistElectronEtaET;
754 delete fHistElectronEtaE;
755 delete fHistElectronEtaPt;
756 delete fHistElectrontotET;
757
758 delete fHistConvElectronEtaEET;
759 delete fHistConvElectronEtaPtET;
760 delete fHistConvElectronEtaET;
761 delete fHistConvElectronEtaE;
762 delete fHistConvElectronEtaPt;
763 delete fHistConvElectrontotET;
764
765 delete fHistScatElectronEtaEET;
766 delete fHistScatElectronEtaPtET;
767 delete fHistScatElectronEtaET;
768 delete fHistScatElectronEtaE;
769 delete fHistScatElectronEtaPt;
770 delete fHistScatElectrontotET;
771
772 // *******************
773 // total electron ET
774 // *******************
775 delete fHistTotElectrontotET;
776
777 // *******************
778 // gamma ET
779 // *******************
780 delete fHistGammaEtaEET;
781 delete fHistGammaEtaPtET;
782 delete fHistGammaEtaET;
783 delete fHistGammaEtaE;
784 delete fHistGammaEtaPt;
785 delete fHistGammatotET;
786
787 delete fHistAnnihGammaEtaEET;
788 delete fHistAnnihGammaEtaPtET;
789 delete fHistAnnihGammaEtaET;
790 delete fHistAnnihGammaEtaE;
791 delete fHistAnnihGammaEtaPt;
792 delete fHistAnnihGammatotET;
793
794 delete fHistScatGammaEtaEET;
795 delete fHistScatGammaEtaPtET;
796 delete fHistScatGammaEtaET;
797 delete fHistScatGammaEtaE;
798 delete fHistScatGammaEtaPt;
799 delete fHistScatGammatotET;
800
801 delete fHistConvGammaEtaEET;
802 delete fHistConvGammaEtaPtET;
803 delete fHistConvGammaEtaET;
804 delete fHistConvGammaEtaE;
805 delete fHistConvGammaEtaPt;
806 delete fHistConvGammatotET;
807
808 delete fHistNonConvGammaEtaEET;
809 delete fHistNonConvGammaEtaPtET;
810 delete fHistNonConvGammaEtaET;
811 delete fHistNonConvGammaEtaE;
812 delete fHistNonConvGammaEtaPt;
813 delete fHistNonConvGammatotET;
814
815 // *******************
816 // total gamma ET
817 // *******************
818 delete fHistTotGammatotET;
819
820 // *******************
821 // total electromagnetic ET
822 // *******************
823 delete fHistTotEMtotET;
824
825 // non-primary electromagnetic ET
826 delete fHistNPPElectronEtaEET;
827 delete fHistNPPElectronEtaPtET;
828 delete fHistNPPElectronEtaET;
829 delete fHistNPPElectronEtaE;
830 delete fHistNPPElectronEtaPt;
831 delete fHistNPPElectrontotET;
832
833 delete fHistNPPGammaEtaEET;
834 delete fHistNPPGammaEtaPtET;
835 delete fHistNPPGammaEtaET;
836 delete fHistNPPGammaEtaE;
837 delete fHistNPPGammaEtaPt;
838 delete fHistNPPGammatotET;
839
840 delete fHistTotNPPEMtotET;
841
842 delete fHistNPPPi0GammaEtaEET;
843 delete fHistNPPPi0GammaEtaPtET;
844 delete fHistNPPPi0GammaEtaET;
845 delete fHistNPPPi0GammaEtaE;
846 delete fHistNPPPi0GammaEtaPt;
847 delete fHistNPPPi0GammatotET;
848
849 // *******************
850 // electron ET inside EMCal acceptance
851 // *******************
852 delete fHistElectronAccEtaEET;
853 delete fHistElectronAccEtaPtET;
854 delete fHistElectronAccEtaET;
855 delete fHistElectronAccEtaE;
856 delete fHistElectronAccEtaPt;
857 delete fHistElectronAcctotET;
858
859 delete fHistConvElectronAccEtaEET;
860 delete fHistConvElectronAccEtaPtET;
861 delete fHistConvElectronAccEtaET;
862 delete fHistConvElectronAccEtaE;
863 delete fHistConvElectronAccEtaPt;
864 delete fHistConvElectronAcctotET;
865
866 delete fHistScatElectronAccEtaEET;
867 delete fHistScatElectronAccEtaPtET;
868 delete fHistScatElectronAccEtaET;
869 delete fHistScatElectronAccEtaE;
870 delete fHistScatElectronAccEtaPt;
871 delete fHistScatElectronAcctotET;
872
873 // *******************
874 // total electron ET inside EMCal acceptance
875 // *******************
876 delete fHistTotElectronAcctotET;
877
878 // *******************
879 // gamma ET inside EMCal acceptance
880 // *******************
881 delete fHistGammaAccEtaEET;
882 delete fHistGammaAccEtaPtET;
883 delete fHistGammaAccEtaET;
884 delete fHistGammaAccEtaE;
885 delete fHistGammaAccEtaPt;
886 delete fHistGammaAcctotET;
887
888 delete fHistAnnihGammaAccEtaEET;
889 delete fHistAnnihGammaAccEtaPtET;
890 delete fHistAnnihGammaAccEtaET;
891 delete fHistAnnihGammaAccEtaE;
892 delete fHistAnnihGammaAccEtaPt;
893 delete fHistAnnihGammaAcctotET;
894
895 delete fHistScatGammaAccEtaEET;
896 delete fHistScatGammaAccEtaPtET;
897 delete fHistScatGammaAccEtaET;
898 delete fHistScatGammaAccEtaE;
899 delete fHistScatGammaAccEtaPt;
900 delete fHistScatGammaAcctotET;
901
902 delete fHistConvGammaAccEtaEET;
903 delete fHistConvGammaAccEtaPtET;
904 delete fHistConvGammaAccEtaET;
905 delete fHistConvGammaAccEtaE;
906 delete fHistConvGammaAccEtaPt;
907 delete fHistConvGammaAcctotET;
908
909 delete fHistNonConvGammaAccEtaEET;
910 delete fHistNonConvGammaAccEtaPtET;
911 delete fHistNonConvGammaAccEtaET;
912 delete fHistNonConvGammaAccEtaE;
913 delete fHistNonConvGammaAccEtaPt;
914 delete fHistNonConvGammaAcctotET;
915
916 // *******************
917 // total gamma ET inside EMCal acceptance
918 // *******************
919 delete fHistTotGammaAcctotET;
920
921 // *******************
922 // total electromagnetic ET inside EMCal acceptance
923 // *******************
924 delete fHistTotEMAcctotET;
925
926 // non-primary electromagnetic ET
927 delete fHistNPPElectronAccEtaEET;
928 delete fHistNPPElectronAccEtaPtET;
929 delete fHistNPPElectronAccEtaE;
930 delete fHistNPPElectronAccEtaPt;
931
932 delete fHistNPPGammaAccEtaEET;
933 delete fHistNPPGammaAccEtaPtET;
934 delete fHistNPPGammaAccEtaE;
935 delete fHistNPPGammaAccEtaPt;
936
937 // *******************
938 // electron ET reconstructed in EMCal
939 // *******************
940 delete fHistElectronRecEtaEET;
941 delete fHistElectronRecEtaPtET;
942 delete fHistElectronRecEtaET;
943 delete fHistElectronRecEtaE;
944 delete fHistElectronRecEtaPt;
945 delete fHistElectronRectotET;
946
947 delete fHistConvElectronRecEtaEET;
948 delete fHistConvElectronRecEtaPtET;
949 delete fHistConvElectronRecEtaET;
950 delete fHistConvElectronRecEtaE;
951 delete fHistConvElectronRecEtaPt;
952 delete fHistConvElectronRectotET;
953
954 delete fHistScatElectronRecEtaEET;
955 delete fHistScatElectronRecEtaPtET;
956 delete fHistScatElectronRecEtaET;
957 delete fHistScatElectronRecEtaE;
958 delete fHistScatElectronRecEtaPt;
959 delete fHistScatElectronRectotET;
960
961 // *******************
962 // total Electron ET reconstructed in EMCal
963 // *******************
964 delete fHistTotElectronRectotET;
965
966 // *******************
967 // gamma ET reconstructed in EMCal
968 // *******************
969 delete fHistGammaRecEtaEET;
970 delete fHistGammaRecEtaPtET;
971 delete fHistGammaRecEtaET;
972 delete fHistGammaRecEtaE;
973 delete fHistGammaRecEtaPt;
974 delete fHistGammaRectotET;
975
976 delete fHistGammaRecResEET;//Add Caio
977 delete fHistGammaRecResPtET;//Add Caio
978 delete fHistGammaRecResE;//Add Caio
979 delete fHistGammaRecResPt;//Add Caio
980 delete fHistGammaRecResEDepETDep;//Add Caio
981 delete fHistGammaRecResPtETDep; // Add Caio
982
983 delete fHistAnnihGammaRecEtaEET;
984 delete fHistAnnihGammaRecEtaPtET;
985 delete fHistAnnihGammaRecEtaET;
986 delete fHistAnnihGammaRecEtaE;
987 delete fHistAnnihGammaRecEtaPt;
988 delete fHistAnnihGammaRectotET;
989
990 delete fHistScatGammaRecEtaEET;
991 delete fHistScatGammaRecEtaPtET;
992 delete fHistScatGammaRecEtaET;
993 delete fHistScatGammaRecEtaE;
994 delete fHistScatGammaRecEtaPt;
995 delete fHistScatGammaRectotET;
996
997 // *******************
998 // total gamma ET reconstructed in EMCal
999 // *******************
1000 delete fHistTotGammaRectotET;
1001
1002 // *******************
1003 // total EM ET reconstructed in EMCal
1004 // *******************
1005 delete fHistTotEMRectotET;
1006
1007 // non-primary electromagnetic ET
1008 delete fHistNPPElectronRecEtaEET;
1009 delete fHistNPPElectronRecEtaPtET;
1010 delete fHistNPPElectronRecEtaET;
1011 delete fHistNPPElectronRecEtaE;
1012 delete fHistNPPElectronRecEtaPt;
1013 delete fHistNPPElectronRectotET;
1014
1015 delete fHistNPPGammaRecEtaEET;
1016 delete fHistNPPGammaRecEtaPtET;
1017 delete fHistNPPGammaRecEtaET;
1018 delete fHistNPPGammaRecEtaE;
1019 delete fHistNPPGammaRecEtaPt;
1020 delete fHistNPPGammaRectotET;
1021
1022 delete fHistTotNPPEMRectotET;
1023
1024 delete fHistNPPPi0GammaRecEtaEET;
1025 delete fHistNPPPi0GammaRecEtaPtET;
1026 delete fHistNPPPi0GammaRecEtaET;
1027 delete fHistNPPPi0GammaRecEtaE;
1028 delete fHistNPPPi0GammaRecEtaPt;
1029 delete fHistNPPPi0GammaRectotET;
1030
1031 // *******************
1032 // muon ET (+ and -)
1033 // *******************
1034 delete fHistMuonEtaEET;
1035 delete fHistMuonAccEtaEET;
1036 delete fHistMuonRecEtaEET;
1037 delete fHistMuonMatchEtaEET;
1038
1039 delete fHistMuonEtaPtET;
1040 delete fHistMuonAccEtaPtET;
1041 delete fHistMuonRecEtaPtET;
1042 delete fHistMuonMatchEtaPtET;
1043
1044 delete fHistMuonEtaET;
1045 delete fHistMuonAccEtaET;
1046 delete fHistMuonRecEtaET;
1047 delete fHistMuonMatchEtaET;
1048
1049 delete fHistMuonEtaE;
1050 delete fHistMuonAccEtaE;
1051 delete fHistMuonRecEtaE;
1052 delete fHistMuonMatchEtaE;
1053
1054 delete fHistMuonEtaPt;
1055 delete fHistMuonAccEtaPt;
1056 delete fHistMuonRecEtaPt;
1057 delete fHistMuonMatchEtaPt;
1058
1059 delete fHistMuontotET;
1060 delete fHistMuonAcctotET;
1061 delete fHistMuonRectotET;
1062 delete fHistMuonMatchtotET;
1063
1064 delete fHistMuonRectotETDep;
1065 delete fHistMuonMatchtotETDep;
1066
1067 delete fHistMuonRecEtaEDepETDep;
1068 delete fHistMuonMatchEtaEDepETDep;
1069
1070 delete fHistMuonRecEtaPtETDep;
1071 delete fHistMuonMatchEtaPtETDep;
1072
1073 delete fHistMuonRecEtaETDep;
1074 delete fHistMuonMatchEtaETDep;
1075
1076 delete fHistMuonRecResEET;
1077 delete fHistMuonRecResPtET;
1078 delete fHistMuonRecResE;
1079 delete fHistMuonRecResPt;
1080 delete fHistMuonRecResEDepETDep;
1081 delete fHistMuonRecResPtETDep;
1082
1083 // *******************
1084 // pion ET (+ and -)
1085 // *******************
1086 delete fHistPionEtaEET;
1087 delete fHistPionAccEtaEET;
1088 delete fHistPionRecEtaEET;
1089 delete fHistPionMatchEtaEET;
1090
1091 delete fHistPionEtaPtET;
1092 delete fHistPionAccEtaPtET;
1093 delete fHistPionRecEtaPtET;
1094 delete fHistPionMatchEtaPtET;
1095
1096 delete fHistPionEtaET;
1097 delete fHistPionAccEtaET;
1098 delete fHistPionRecEtaET;
1099 delete fHistPionMatchEtaET;
1100
1101 delete fHistPionEtaE;
1102 delete fHistPionAccEtaE;
1103 delete fHistPionRecEtaE;
1104 delete fHistPionMatchEtaE;
1105
1106 delete fHistPionEtaPt;
1107 delete fHistPionAccEtaPt;
1108 delete fHistPionRecEtaPt;
1109 delete fHistPionMatchEtaPt;
1110
1111 delete fHistPiontotET;
1112 delete fHistPionAcctotET;
1113 delete fHistPionRectotET;
1114 delete fHistPionMatchtotET;
1115
1116 delete fHistPionRectotETDep;
1117 delete fHistPionMatchtotETDep;
1118
1119 delete fHistPionRecEtaEDepETDep;
1120 delete fHistPionMatchEtaEDepETDep;
1121
1122 delete fHistPionRecEtaPtETDep;
1123 delete fHistPionMatchEtaPtETDep;
1124
1125 delete fHistPionRecEtaETDep;
1126 delete fHistPionMatchEtaETDep;
1127
1128 delete fHistPionRecResEET;
1129 delete fHistPionRecResPtET;
1130 delete fHistPionRecResE;
1131 delete fHistPionRecResPt;
1132 delete fHistPionRecResEDepETDep;
1133 delete fHistPionRecResPtETDep;
1134
1135 // *******************
1136 // charged kaon (+ and -) ET
1137 // *******************
1138 delete fHistKaonEtaEET;
1139 delete fHistKaonAccEtaEET;
1140 delete fHistKaonRecEtaEET;
1141 delete fHistKaonMatchEtaEET;
1142
1143 delete fHistKaonEtaPtET;
1144 delete fHistKaonAccEtaPtET;
1145 delete fHistKaonRecEtaPtET;
1146 delete fHistKaonMatchEtaPtET;
1147
1148 delete fHistKaonEtaET;
1149 delete fHistKaonAccEtaET;
1150 delete fHistKaonRecEtaET;
1151 delete fHistKaonMatchEtaET;
1152
1153 delete fHistKaonEtaE;
1154 delete fHistKaonAccEtaE;
1155 delete fHistKaonRecEtaE;
1156 delete fHistKaonMatchEtaE;
1157
1158 delete fHistKaonEtaPt;
1159 delete fHistKaonAccEtaPt;
1160 delete fHistKaonRecEtaPt;
1161 delete fHistKaonMatchEtaPt;
1162
1163 delete fHistKaontotET;
1164 delete fHistKaonAcctotET;
1165 delete fHistKaonRectotET;
1166 delete fHistKaonMatchtotET;
1167
1168 delete fHistKaonRectotETDep;
1169 delete fHistKaonMatchtotETDep;
1170
1171 delete fHistKaonRecEtaEDepETDep;
1172 delete fHistKaonMatchEtaEDepETDep;
1173
1174 delete fHistKaonRecEtaPtETDep;
1175 delete fHistKaonMatchEtaPtETDep;
1176
1177 delete fHistKaonRecEtaETDep;
1178 delete fHistKaonMatchEtaETDep;
1179
1180 delete fHistKaonRecResEET;
1181 delete fHistKaonRecResPtET;
1182 delete fHistKaonRecResE;
1183 delete fHistKaonRecResPt;
1184 delete fHistKaonRecResEDepETDep;
1185 delete fHistKaonRecResPtETDep;
1186
1187 // *******************
1188 // proton (anti) ET
1189 // *******************
1190 delete fHistProtonEtaEET;
1191 delete fHistProtonAccEtaEET;
1192 delete fHistProtonRecEtaEET;
1193 delete fHistProtonMatchEtaEET;
1194
1195 delete fHistProtonEtaPtET;
1196 delete fHistProtonAccEtaPtET;
1197 delete fHistProtonRecEtaPtET;
1198 delete fHistProtonMatchEtaPtET;
1199
1200 delete fHistProtonEtaET;
1201 delete fHistProtonAccEtaET;
1202 delete fHistProtonRecEtaET;
1203 delete fHistProtonMatchEtaET;
1204
1205 delete fHistProtonEtaE;
1206 delete fHistProtonAccEtaE;
1207 delete fHistProtonRecEtaE;
1208 delete fHistProtonMatchEtaE;
1209
1210 delete fHistProtonEtaPt;
1211 delete fHistProtonAccEtaPt;
1212 delete fHistProtonRecEtaPt;
1213 delete fHistProtonMatchEtaPt;
1214
1215 delete fHistProtontotET;
1216 delete fHistProtonAcctotET;
1217 delete fHistProtonRectotET;
1218 delete fHistProtonMatchtotET;
1219
1220 delete fHistProtonRectotETDep;
1221 delete fHistProtonMatchtotETDep;
1222
1223 delete fHistProtonRecEtaEDepETDep;
1224 delete fHistProtonMatchEtaEDepETDep;
1225
1226 delete fHistProtonRecEtaPtETDep;
1227 delete fHistProtonMatchEtaPtETDep;
1228
1229 delete fHistProtonRecEtaETDep;
1230 delete fHistProtonMatchEtaETDep;
1231
1232 delete fHistProtonRecResEET;
1233 delete fHistProtonRecResPtET;
1234 delete fHistProtonRecResE;
1235 delete fHistProtonRecResPt;
1236 delete fHistProtonRecResEDepETDep;
1237 delete fHistProtonRecResPtETDep;
1238
1239 // *******************
1240 // total charged ET
1241 // *******************
1242 delete fHistTotChargedtotET;
1243 delete fHistTotChargedAcctotET;
1244 delete fHistTotChargedRectotET;
1245 delete fHistTotChargedRectotETDep;
1246 delete fHistTotChargedMatchtotET;
1247 delete fHistTotChargedMatchtotETDep;
1248
1249 // *******************
1250 // neutron (anti) ET
1251 // *******************
1252 delete fHistNeutronEtaEET;
1253 delete fHistNeutronAccEtaEET;
1254 delete fHistNeutronRecEtaEET;
1255
1256 delete fHistNeutronEtaPtET;
1257 delete fHistNeutronAccEtaPtET;
1258 delete fHistNeutronRecEtaPtET;
1259
1260 delete fHistNeutronEtaET;
1261 delete fHistNeutronAccEtaET;
1262 delete fHistNeutronRecEtaET;
1263
1264 delete fHistNeutronEtaE;
1265 delete fHistNeutronAccEtaE;
1266 delete fHistNeutronRecEtaE;
1267
1268 delete fHistNeutronEtaPt;
1269 delete fHistNeutronAccEtaPt;
1270 delete fHistNeutronRecEtaPt;
1271
1272 delete fHistNeutrontotET;
1273 delete fHistNeutronAcctotET;
1274 delete fHistNeutronRectotET;
1275 delete fHistNeutronRectotETDep;
1276
1277 delete fHistNeutronRecEtaEDepETDep;
1278 delete fHistNeutronRecEtaETDep;
1279
1280 delete fHistNeutronRecEtaPtETDep;
1281
1282 // *******************
1283 // neutral kaon ET
1284 // *******************
1285 delete fHistK0EtaEET;
1286 delete fHistK0RecEtaEET;
1287
1288 delete fHistK0EtaPtET;
1289 delete fHistK0RecEtaPtET;
1290
1291 delete fHistK0EtaET;
1292 delete fHistK0RecEtaET;
1293
1294 delete fHistK0EtaE;
1295 delete fHistK0RecEtaE;
1296
1297 delete fHistK0EtaPt;
1298 delete fHistK0RecEtaPt;
1299
1300 delete fHistK0totET;
1301 delete fHistK0RectotET;
1302
1303 delete fHistK0RectotETDep;
1304
1305 delete fHistK0RecEtaEDepETDep;
1306 delete fHistK0RecEtaETDep;
1307
1308 delete fHistK0RecEtaPtETDep;
1309
1310 // *******************
1311 // Lambda(anti) ET
1312 // *******************
1313 delete fHistLambdaEtaEET;
1314 delete fHistLambdaRecEtaEET;
1315
1316 delete fHistLambdaEtaPtET;
1317 delete fHistLambdaRecEtaPtET;
1318
1319 delete fHistLambdaEtaET;
1320 delete fHistLambdaRecEtaET;
1321
1322 delete fHistLambdaEtaE;
1323 delete fHistLambdaRecEtaE;
1324
1325 delete fHistLambdaEtaPt;
1326 delete fHistLambdaRecEtaPt;
1327
1328 delete fHistLambdatotET;
1329 delete fHistLambdaRectotET;
1330
1331 delete fHistLambdaRectotETDep;
1332
1333 delete fHistLambdaRecEtaEDepETDep;
1334 delete fHistLambdaRecEtaETDep;
1335
1336 delete fHistLambdaRecEtaPtETDep;
1337
1338 // *******************
1339 // total neutral ET
1340 // *******************
1341 delete fHistTotNeutraltotET;
1342 delete fHistTotNeutralRectotET;
1343 delete fHistTotNeutralRectotETDep;
1344
1345 // *******************
1346 // total ET
1347 // *******************
1348 delete fHistTotaltotET;
1349 delete fHistTotalAcctotET;
1350 delete fHistTotalRectotET;
1351 delete fHistTotalRectotETDep;
1352
1353 // *******************
1354 // some checks
1355 // *******************
1356
1357 // check produced electrons
1358 delete fHistElectronFirstMother;
1359 delete fHistElectronFirstMotherXY;
1360 delete fHistElectronNDaughters;
1361 delete fHistElectronDaughters;
1362 delete fHistElectronDaughtersXY;
1363
1364 delete fHistElectronFirstMotherAcc;
1365 delete fHistElectronFirstMotherXYAcc;
1366 delete fHistElectronNDaughtersAcc;
1367 delete fHistElectronDaughtersAcc;
1368 delete fHistElectronDaughtersXYAcc;
1369
1370 delete fHistElectronFirstMotherRec;
1371 delete fHistElectronFirstMotherXYRec;
1372 delete fHistElectronNDaughtersRec;
1373 delete fHistElectronDaughtersRec;
1374 delete fHistElectronDaughtersXYRec;
1375
1376 delete fHistNPPElectronFirstMother;
1377 delete fHistNPPElectronFirstMotherXY;
1378 delete fHistNPPElectronNDaughters;
1379 delete fHistNPPElectronDaughters;
1380 delete fHistNPPElectronDaughtersXY;
1381
1382 delete fHistNPPElectronFirstMotherAcc;
1383 delete fHistNPPElectronFirstMotherXYAcc;
1384 delete fHistNPPElectronNDaughtersAcc;
1385 delete fHistNPPElectronDaughtersAcc;
1386 delete fHistNPPElectronDaughtersXYAcc;
1387
1388 delete fHistNPPElectronFirstMotherRec;
1389 delete fHistNPPElectronFirstMotherXYRec;
1390 delete fHistNPPElectronNDaughtersRec;
1391 delete fHistNPPElectronDaughtersRec;
1392 delete fHistNPPElectronDaughtersXYRec;
1393
1394 // check produced gammas
1395 delete fHistGammaFirstMother;
1396 delete fHistGammaFirstMotherXY;
1397 delete fHistGammaNDaughters;
1398 delete fHistGammaDaughters;
1399 delete fHistGammaDaughtersXY;
1400 delete fHistConvGammaDaughtersXY;
1401 delete fHistNonConvGammaDaughtersXY;
1402
1403 delete fHistGammaFirstMotherAcc;
1404 delete fHistGammaFirstMotherXYAcc;
1405 delete fHistGammaNDaughtersAcc;
1406 delete fHistGammaDaughtersAcc;
1407 delete fHistGammaDaughtersXYAcc;
1408 delete fHistConvGammaDaughtersXYAcc;
1409 delete fHistNonConvGammaDaughtersXYAcc;
1410
1411 delete fHistGammaFirstMotherRec;
1412 delete fHistGammaFirstMotherXYRec;
1413 delete fHistGammaNDaughtersRec;
1414 delete fHistGammaDaughtersRec;
1415 delete fHistGammaDaughtersXYRec;
1416 delete fHistConvGammaDaughtersXYRec;
1417 delete fHistNonConvGammaDaughtersXYRec;
1418
1419 delete fHistNPPGammaFirstMother;
1420 delete fHistNPPGammaFirstMotherXY;
1421 delete fHistNPPGammaNDaughters;
1422 delete fHistNPPGammaDaughters;
1423 delete fHistNPPGammaDaughtersXY;
1424
1425 delete fHistNPPGammaFirstMotherAcc;
1426 delete fHistNPPGammaFirstMotherXYAcc;
1427 delete fHistNPPGammaNDaughtersAcc;
1428 delete fHistNPPGammaDaughtersAcc;
1429 delete fHistNPPGammaDaughtersXYAcc;
1430
1431 delete fHistNPPGammaFirstMotherRec;
1432 delete fHistNPPGammaFirstMotherXYRec;
1433 delete fHistNPPGammaNDaughtersRec;
1434 delete fHistNPPGammaDaughtersRec;
1435 delete fHistNPPGammaDaughtersXYRec;
1436
1437 //check projections
1438 delete fHistAllERecEMC;
1439 delete fHistAllPtRecPtMC;
1440 delete fHistElectronERecEMC;
1441 delete fHistGammaERecEMC;
1442
1443 delete fHistChargedRes;
1444 delete fHistChargedRes2;
1445 delete fHistChargedRes3;
1446 delete fHistNeutralRes;
1447 delete fHistElectronRes;
1448 delete fHistGammaRes;
1449
1450 delete fHistIsInAcc;
1451 //delete TH2F * yyyyyy
2e12941b 1452}
1453
1454Int_t AliAnalysisEmEtMonteCarlo::AnalyseEvent(AliVEvent* ev)
1455{ // analyse MC event
112b1d0c 1456 //ResetEventValues();
2e12941b 1457
112b1d0c 1458 // Get us an mc event
1459 if(!ev)
f2adae2b 1460 {
4d376d01 1461 Printf("ERROR: ev does not exist");
112b1d0c 1462 return 0;
2e12941b 1463 }
112b1d0c 1464 AliMCEvent *event = dynamic_cast<AliMCEvent*>(ev);
4d376d01 1465 if(!event)
1466 {
1467 Printf("ERROR: event does not exist");
1468 return 0;
1469 }
112b1d0c 1470
1471 // Hijing header
1472 AliGenEventHeader* genHeader = event->GenEventHeader();
1473 if(!genHeader){
1474 Printf("ERROR: Event generation header does not exist");
1475 return 0;
1476 }
1477
1478 // Let's play with the stack!
1479 AliStack *stack = event->Stack();
1480
1481 if (!stack)
f2adae2b 1482 {
112b1d0c 1483 Printf("ERROR: Could not get stack");
1484 return 0;
f2adae2b 1485 }
2e12941b 1486
112b1d0c 1487 //Int_t nStackTracks = stack->GetNtrack();
1488 Int_t nStackTracks = event->GetNumberOfTracks();
2e12941b 1489
112b1d0c 1490 for (Int_t iPart = 0; iPart < nStackTracks; iPart++){
1491 AliMCParticle* aliPart = (AliMCParticle*)event->GetTrack(iPart);
1492
1493 //TParticle *part = stack->Particle(iPart);
1494 TParticle *part = aliPart->Particle();
1495 TParticle *partMom = 0;
1496 TParticle *partDaughter = 0;
2e12941b 1497
112b1d0c 1498 if (!part){
1499 Printf("ERROR: Could not get particle %d", iPart);
1500 continue;
1501 }
1502
1503 Int_t iPartMom = part->GetMother(0);
1504 Int_t iPartDaughter = 0;
1505 Int_t nPartDaughters = part->GetNDaughters();
1506
1507 TParticlePDG *pdg = part->GetPDG(0);
1508 TParticlePDG *pdgMom = 0;
1509 TParticlePDG *pdgDaugther = 0;
1510
1511 if (!pdg){
1512 Printf("ERROR-1: Could not get particle PDG %d", iPart);
1513 continue;
1514 }
2e12941b 1515
112b1d0c 1516 //create an external track param for projection
1517 AliExternalTrackParam* extParam = CreateExternalTrackParam(part);
1518
1519 if ((iPartMom>=0) && (iPartMom < nStackTracks))
1520 {
1521 partMom = stack->Particle(iPartMom);
1522 pdgMom = partMom->GetPDG(0);
1523 }
2e12941b 1524
112b1d0c 1525 // Check if it is a primary particle
2e12941b 1526
112b1d0c 1527 // Check for reasonable (for now neutral and singly charged) charge on the particle
1528 //TODO:Maybe not only singly charged?
1529 if (TMath::Abs(TMath::Abs(pdg->Charge()) - fCuts->GetMonteCarloSingleChargedParticle())<1e-3 && TMath::Abs(TMath::Abs(pdg->Charge()) - fCuts->GetMonteCarloNeutralParticle())<1e-3) continue;
2e12941b 1530
112b1d0c 1531 fMultiplicity++;
2e12941b 1532
112b1d0c 1533 // Inside ALICE central barrel acceptance
2e12941b 1534
112b1d0c 1535 if (TMath::Abs(part->Eta()) < fCuts->GetCommonEtaCut())
1536 {
1537 Double_t et = CalcET(part,pdg);
2e12941b 1538
112b1d0c 1539 if (et < 0) continue;
1540
1541 if (IsPrimary(stack,iPart,pdg,iPartMom,pdgMom))
1542 {
1543 if (stack->IsPhysicalPrimary(iPart))
1544 {
1545 fHistPrimEtaEET->Fill(part->Energy(),part->Eta(),et);
1546 fHistPrimEtaPtET->Fill(part->Pt(),part->Eta(),et);
1547 fHistPrimEtaET->Fill(et,part->Eta());
1548 fPrimtotET += et;
2e12941b 1549
112b1d0c 1550 //if (IsInAcceptance(part,pdg,extParam))
1551 if (IsInAcceptance(aliPart))
1552 {
1553 fHistPrimAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1554 fHistPrimAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1555 fHistPrimAccEtaET->Fill(et,part->Eta());
1556 fPrimAcctotET += et;
1557 }
1558 }
1559
1560 // Fill up total E_T counters for each particle species
1561 if (pdg->PdgCode() == fgProtonCode || pdg->PdgCode() == fgAntiProtonCode)
1562 {
1563 fProtontotET += et;
1564 fHistProtonEtaEET->Fill(part->Energy(),part->Eta(),et);
1565 fHistProtonEtaPtET->Fill(part->Pt(),part->Eta(),et);
1566 fHistProtonEtaET->Fill(et,part->Eta());
1567 fHistProtonEtaE->Fill(part->Energy(),part->Eta());
1568 fHistProtonEtaPt->Fill(part->Pt(),part->Eta());
2e12941b 1569
112b1d0c 1570 // inside EMCal acceptance
1571 //if (IsInAcceptance(part,pdg,extParam))
1572 if (IsInAcceptance(aliPart))
1573 {
1574 fProtonAcctotET += et;
1575 fHistProtonAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1576 fHistProtonAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1577 fHistProtonAccEtaET->Fill(et,part->Eta());
1578 fHistProtonAccEtaE->Fill(part->Energy(),part->Eta());
1579 fHistProtonAccEtaPt->Fill(part->Pt(),part->Eta());
1580 }
1581 }
1582 if (pdg->PdgCode() == fgPiPlusCode || pdg->PdgCode() == fgPiMinusCode)
1583 {
1584 fPiontotET += et;
1585 fHistPionEtaEET->Fill(part->Energy(),part->Eta(),et);
1586 fHistPionEtaPtET->Fill(part->Pt(),part->Eta(),et);
1587 fHistPionEtaET->Fill(et,part->Eta());
1588 fHistPionEtaE->Fill(part->Energy(),part->Eta());
1589 fHistPionEtaPt->Fill(part->Pt(),part->Eta());
1590 // inside EMCal acceptance
1591 //if (IsInAcceptance(part,pdg,extParam))
1592 if (IsInAcceptance(aliPart))
1593 {
1594 fPionAcctotET += et;
1595 fHistPionAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1596 fHistPionAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1597 fHistPionAccEtaET->Fill(et,part->Eta());
1598 fHistPionAccEtaE->Fill(part->Energy(),part->Eta());
1599 fHistPionAccEtaPt->Fill(part->Pt(),part->Eta());
1600 }
1601 }
1602 if (pdg->PdgCode() == fgKPlusCode || pdg->PdgCode() == fgKMinusCode)
1603 {
1604 fKaontotET += et;
1605 fHistKaonEtaEET->Fill(part->Energy(),part->Eta(),et);
1606 fHistKaonEtaPtET->Fill(part->Pt(),part->Eta(),et);
1607 fHistKaonEtaE->Fill(part->Energy(),part->Eta());
1608 fHistKaonEtaET->Fill(et,part->Eta());
1609 fHistKaonEtaPt->Fill(part->Pt(),part->Eta());
1610 // inside EMCal acceptance
1611 //if (IsInAcceptance(part,pdg,extParam))
1612 if (IsInAcceptance(aliPart))
1613 {
1614 fKaonAcctotET += et;
1615 fHistKaonAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1616 fHistKaonAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1617 fHistKaonAccEtaET->Fill(et,part->Eta());
1618 fHistKaonAccEtaE->Fill(part->Energy(),part->Eta());
1619 fHistKaonAccEtaPt->Fill(part->Pt(),part->Eta());
1620 }
1621 }
1622 if (pdg->PdgCode() == fgMuPlusCode || pdg->PdgCode() == fgMuMinusCode)
1623 {
1624 fMuontotET += et;
1625 fHistMuonEtaEET->Fill(part->Energy(),part->Eta(),et);
1626 fHistMuonEtaPtET->Fill(part->Pt(),part->Eta(),et);
1627 fHistMuonEtaET->Fill(et,part->Eta());
1628 fHistMuonEtaE->Fill(part->Energy(),part->Eta());
1629 fHistMuonEtaPt->Fill(part->Pt(),part->Eta());
1630 // inside EMCal acceptance
1631 //if (IsInAcceptance(part,pdg,extParam))
1632 if (IsInAcceptance(aliPart))
1633 {
1634 fMuonAcctotET += et;
1635 fHistMuonAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1636 fHistMuonAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1637 fHistMuonAccEtaET->Fill(et,part->Eta());
1638 fHistMuonAccEtaE->Fill(part->Energy(),part->Eta());
1639 fHistMuonAccEtaPt->Fill(part->Pt(),part->Eta());
1640 }
1641 }
1642 if (pdg->PdgCode() == fgEPlusCode || pdg->PdgCode() == fgEMinusCode)
1643 {
1644 if (stack->IsPhysicalPrimary(iPart))
1645 {//Marcelo - isn't this redundant? Isn't this inside an if statement already?
1646 fElectrontotET += et;
1647 fHistElectronEtaEET->Fill(part->Energy(),part->Eta(),et);
1648 fHistElectronEtaPtET->Fill(part->Pt(),part->Eta(),et);
1649 fHistElectronEtaET->Fill(et,part->Eta());
1650 fHistElectronEtaE->Fill(part->Energy(),part->Eta());
1651 fHistElectronEtaPt->Fill(part->Pt(),part->Eta());
2e12941b 1652
112b1d0c 1653 // inside EMCal acceptance
1654 //if (IsInAcceptance(part,pdg,extParam))
1655 if (IsInAcceptance(aliPart))
1656 {
1657 fElectronAcctotET += et;
1658 fHistElectronAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1659 fHistElectronAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1660 fHistElectronAccEtaET->Fill(et,part->Eta());
1661 fHistElectronAccEtaE->Fill(part->Energy(),part->Eta());
1662 fHistElectronAccEtaPt->Fill(part->Pt(),part->Eta());
1663 }
1664 }
1665 else if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz()))
1666 {//Marcelo - are we sure we know what this is doing? How sensitive is this to geometry?
1667 if (IsMotherPrimaryGamma(stack,iPartMom,pdgMom))
1668 {
1669 fHistConvElectronEtaEET->Fill(part->Energy(),part->Eta(),et);
1670 fHistConvElectronEtaPtET->Fill(part->Pt(),part->Eta(),et);
1671 fHistConvElectronEtaET->Fill(et,part->Eta());
1672 fHistConvElectronEtaE->Fill(part->Energy(),part->Eta());
1673 fHistConvElectronEtaPt->Fill(part->Pt(),part->Eta());
1674 fConvElectrontotET += et;
2e12941b 1675
112b1d0c 1676 // gamma mother is inside EMCal acceptance
1677 //if (IsInAcceptance(partMom,pdgMom))
1678 if (IsInAcceptance(aliPart))
1679 {
1680 fHistConvElectronAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1681 fHistConvElectronAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1682 fHistConvElectronAccEtaET->Fill(et,part->Eta());
1683 fHistConvElectronAccEtaE->Fill(part->Energy(),part->Eta());
1684 fHistConvElectronAccEtaPt->Fill(part->Pt(),part->Eta());
1685 fConvElectronAcctotET += et;
1686 }
1687 }
1688 else if (IsMotherPrimaryElectron(stack,iPartMom,pdgMom))
1689 {
1690 fHistScatElectronEtaEET->Fill(part->Energy(),part->Eta(),et);
1691 fHistScatElectronEtaPtET->Fill(part->Pt(),part->Eta(),et);
1692 fHistScatElectronEtaET->Fill(et,part->Eta());
1693 fHistScatElectronEtaE->Fill(part->Energy(),part->Eta());
1694 fHistScatElectronEtaPt->Fill(part->Pt(),part->Eta());
1695 fScatElectrontotET += et;
1696
1697 // inside EMCal acceptance - does it work?
1698 //if (IsInAcceptance(part,pdg,extParam))
1699 if (IsInAcceptance(aliPart))
1700 {
1701 fHistScatElectronAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1702 fHistScatElectronAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1703 fHistScatElectronAccEtaET->Fill(et,part->Eta());
1704 fHistScatElectronAccEtaE->Fill(part->Energy(),part->Eta());
1705 fHistScatElectronAccEtaPt->Fill(part->Pt(),part->Eta());
1706 fScatElectronAcctotET += et;
1707 }
1708 }
1709 }
1710
1711 // few checks
1712 if ((stack->IsPhysicalPrimary(iPart)) || (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())))
1713 {//Marcelo - ...isn't this redundant?
1714 if (pdgMom)
1715 fHistElectronFirstMother->Fill(pdgMom->PdgCode());
1716 fHistElectronFirstMotherXY->Fill(part->Vx(),part->Vy());
1717 fHistElectronNDaughters->Fill(nPartDaughters);
2e12941b 1718
112b1d0c 1719 iPartDaughter = part->GetLastDaughter();
1720 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
1721 {
1722 partDaughter = stack->Particle(iPartDaughter);
1723 if (partDaughter)
1724 {
1725 pdgDaugther = partDaughter->GetPDG(0);
1726 if (pdgDaugther) {
1727 fHistElectronDaughters->Fill(pdgDaugther->PdgCode());
1728 fHistElectronDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
1729 }
1730 }
1731 }
2e12941b 1732
112b1d0c 1733 // inside EMCal acceptance
1734 //if (IsInAcceptance(part,pdg,extParam))
1735 if (IsInAcceptance(aliPart))
1736 {//Marcelo - should all three of the lines below be in the if statement?
1737 if (pdgMom) fHistElectronFirstMotherAcc->Fill(pdgMom->PdgCode());
1738 fHistElectronFirstMotherXYAcc->Fill(part->Vx(),part->Vy());
1739 fHistElectronNDaughtersAcc->Fill(nPartDaughters);
2e12941b 1740
112b1d0c 1741 iPartDaughter = part->GetLastDaughter();
1742 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
1743 {
1744 partDaughter = stack->Particle(iPartDaughter);
1745 if (partDaughter)
1746 {
1747 pdgDaugther = partDaughter->GetPDG(0);
1748 if (pdgDaugther) {
1749 fHistElectronDaughtersAcc->Fill(pdgDaugther->PdgCode());
1750 fHistElectronDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
1751 }
1752 }
1753 }
1754 }
1755 }
1756 } // end of if electron
1757
1758 // some neutrals also
1759 if (pdg->PdgCode() == fgNeutronCode || pdg->PdgCode() == fgAntiNeutronCode)
1760 {
1761 fHistNeutronEtaEET->Fill(part->Energy(),part->Eta(),et);
1762 fHistNeutronEtaPtET->Fill(part->Pt(),part->Eta(),et);
1763 fHistNeutronEtaET->Fill(et,part->Eta());
1764 fHistNeutronEtaE->Fill(part->Energy(),part->Eta());
1765 fHistNeutronEtaPt->Fill(part->Pt(),part->Eta());
1766 fNeutrontotET += et;
f2adae2b 1767
112b1d0c 1768 // inside EMCal acceptance
1769 //if (IsInAcceptance(part,pdg))
1770 if (IsInAcceptance(aliPart))
1771 {
1772 fHistNeutronAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1773 fHistNeutronAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1774 fHistNeutronAccEtaET->Fill(et,part->Eta());
1775 fHistNeutronAccEtaE->Fill(part->Energy(),part->Eta());
1776 fHistNeutronAccEtaPt->Fill(part->Pt(),part->Eta());
1777 fNeutronAcctotET += et;
1778 }
1779 if(pdg->PdgCode() == fgNeutronCode)
1780 {
1781 fNeutronEt += et;
1782 }
1783 if(pdg->PdgCode() == fgAntiNeutronCode)
1784 {
1785 fAntiNeutronEt += et;
1786 }
1787 }
2e12941b 1788
112b1d0c 1789 if(pdg->PdgCode() == fgGammaCode)
1790 {
1791 if (stack->IsPhysicalPrimary(iPart))
1792 {
1793 fHistGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
1794 fHistGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
1795 fHistGammaEtaET->Fill(et,part->Eta());
1796 fHistGammaEtaE->Fill(part->Energy(),part->Eta());
1797 fHistGammaEtaPt->Fill(part->Pt(),part->Eta());
1798 fGammatotET += et;
2e12941b 1799
112b1d0c 1800 if (IsGammaConversion(stack, part, pdg))
1801 {
1802 fHistConvGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
1803 fHistConvGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
1804 fHistConvGammaEtaET->Fill(et,part->Eta());
1805 fHistConvGammaEtaE->Fill(part->Energy(),part->Eta());
1806 fHistConvGammaEtaPt->Fill(part->Pt(),part->Eta());
1807 fConvGammatotET += et;
1808 }
1809 else
1810 {
1811 fHistNonConvGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
1812 fHistNonConvGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
1813 fHistNonConvGammaEtaET->Fill(et,part->Eta());
1814 fHistNonConvGammaEtaE->Fill(part->Energy(),part->Eta());
1815 fHistNonConvGammaEtaPt->Fill(part->Pt(),part->Eta());
1816 fNonConvGammatotET += et;
1817 }
2e12941b 1818
112b1d0c 1819 Bool_t inAcc=kFALSE;
1820 // inside EMCal acceptance
1821 //if (IsInAcceptance(part,pdg))
1822 if (IsInAcceptance(aliPart))
1823 {
1824 //Printf("phi(1) = %f, eta(1) = %f",part->Phi(),part->Eta());
1825 inAcc = kTRUE;
2e12941b 1826
112b1d0c 1827 fHistGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1828 fHistGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1829 fHistGammaAccEtaET->Fill(et,part->Eta());
1830 fHistGammaAccEtaE->Fill(part->Energy(),part->Eta());
1831 fHistGammaAccEtaPt->Fill(part->Pt(),part->Eta());
1832 fGammaAcctotET += et;
1833 }
2e12941b 1834
112b1d0c 1835 if (IsInAcceptance(part,pdg))
1836 {
1837 if (IsGammaConversion(stack, part, pdg))
1838 {
1839 if (inAcc)
1840 Printf("phi(1) = %f, eta(1) = %f",part->Phi(),part->Eta());
1841
1842 fHistConvGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1843 fHistConvGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1844 fHistConvGammaAccEtaET->Fill(et,part->Eta());
1845 fHistConvGammaAccEtaE->Fill(part->Energy(),part->Eta());
1846 fHistConvGammaAccEtaPt->Fill(part->Pt(),part->Eta());
1847 fConvGammaAcctotET += et;
1848 }
1849 else
1850 {
1851 /*
1852 if (!inAcc)
1853 {
1854 Printf("phi(2) = %f, eta(2) = %f, pt = %f",TMath::RadToDeg()*part->Phi(),part->Eta(),part->Pt());
1855
1856 Int_t iPartDaughter = part->GetLastDaughter();
1857
1858 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
1859 {
1860 TParticle *partDaughter = stack->Particle(iPartDaughter);
1861 if (partDaughter)
1862 {
1863 TParticlePDG *pdgDaugther = partDaughter->GetPDG(0);
1864 if (pdgDaugther)
1865 {
1866 Double_t decayR = sqrt(pow(partDaughter->Vx(),2)+pow(partDaughter->Vy(),2));
1867 Printf("radius = %f, daughter pid = %d",decayR,pdgDaugther->PdgCode());
1868 }
1869 }
1870 }
1871
1872 for (int i=0;i<aliPart->GetNumberOfTrackReferences();i++)
1873 {
1874 AliTrackReference* aliTrkRef = aliPart->GetTrackReference(i);
1875
1876 if (aliTrkRef)
1877 {
1878 Printf("det id = %d, x=%f, y=%f, z=%f", aliTrkRef->DetectorId(),aliTrkRef->X(),aliTrkRef->Y(),aliTrkRef->Z());
1879 }
1880 }
1881
1882 }
1883 */
2e12941b 1884
112b1d0c 1885 fHistNonConvGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1886 fHistNonConvGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1887 fHistNonConvGammaAccEtaET->Fill(et,part->Eta());
1888 fHistNonConvGammaAccEtaE->Fill(part->Energy(),part->Eta());
1889 fHistNonConvGammaAccEtaPt->Fill(part->Pt(),part->Eta());
1890 fNonConvGammaAcctotET += et;
1891 }
1892 }
1893 }
1894 else if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz()))
1895 {
1896 if (IsMotherPrimaryElectron(stack,iPartMom,pdgMom))
1897 {
1898 fHistAnnihGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
1899 fHistAnnihGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
1900 fHistAnnihGammaEtaET->Fill(et,part->Eta());
1901 fHistAnnihGammaEtaE->Fill(part->Energy(),part->Eta());
1902 fHistAnnihGammaEtaPt->Fill(part->Pt(),part->Eta());
1903 fAnnihGammatotET += et;
2e12941b 1904
112b1d0c 1905 // inside EMCal acceptance
1906 //if (IsInAcceptance(part,pdg))
1907 if (IsInAcceptance(aliPart))
1908 {
1909 fHistAnnihGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1910 fHistAnnihGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1911 fHistAnnihGammaAccEtaET->Fill(et,part->Eta());
1912 fHistAnnihGammaAccEtaE->Fill(part->Energy(),part->Eta());
1913 fHistAnnihGammaAccEtaPt->Fill(part->Pt(),part->Eta());
1914 fAnnihGammaAcctotET += et;
1915 }
1916 }
1917 else if (IsMotherPrimaryGamma(stack,iPartMom,pdgMom))
1918 {
1919 fHistScatGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
1920 fHistScatGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
1921 fHistScatGammaEtaET->Fill(et,part->Eta());
1922 fHistScatGammaEtaE->Fill(part->Energy(),part->Eta());
1923 fHistScatGammaEtaPt->Fill(part->Pt(),part->Eta());
1924 fScatGammatotET += et;
2e12941b 1925
112b1d0c 1926 // inside EMCal acceptance
1927 //if (IsInAcceptance(part,pdg))
1928 if (IsInAcceptance(aliPart))
1929 {
1930 fHistScatGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
1931 fHistScatGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
1932 fHistScatGammaAccEtaET->Fill(et,part->Eta());
1933 fHistScatGammaAccEtaE->Fill(part->Energy(),part->Eta());
1934 fHistScatGammaAccEtaPt->Fill(part->Pt(),part->Eta());
1935 fScatGammaAcctotET += et;
1936 }
1937 }
1938 }
1939
1940 // few checks
1941 if ((stack->IsPhysicalPrimary(iPart)) || (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())))
1942 {
1943 if (pdgMom)
1944 fHistGammaFirstMother->Fill(pdgMom->PdgCode());
1945 fHistGammaFirstMotherXY->Fill(part->Vx(),part->Vy());
1946 fHistGammaNDaughters->Fill(nPartDaughters);
2e12941b 1947
112b1d0c 1948 iPartDaughter = part->GetLastDaughter();
1949 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
1950 {
1951 partDaughter = stack->Particle(iPartDaughter);
1952 if (partDaughter)
1953 {
1954 pdgDaugther = partDaughter->GetPDG(0);
1955 if (pdgDaugther) {
1956 fHistGammaDaughters->Fill(pdgDaugther->PdgCode());
1957 fHistGammaDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
2e12941b 1958
112b1d0c 1959 if (stack->IsPhysicalPrimary(iPart))
1960 {
1961 if (IsGammaConversion(stack, part, pdg))
1962 {
1963 fHistConvGammaDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
1964 }
1965 else
1966 {
1967 fHistNonConvGammaDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
1968 }
1969 }
1970 }
1971 }
1972 }
2e12941b 1973
112b1d0c 1974 // inside EMCal acceptance
1975 //if (IsInAcceptance(part,pdg))
1976 if (IsInAcceptance(aliPart))
1977 {
1978 if (pdgMom)//Marcelo - again, should this be all three?
1979 fHistGammaFirstMotherAcc->Fill(pdgMom->PdgCode());
1980 fHistGammaFirstMotherXYAcc->Fill(part->Vx(),part->Vy());
1981 fHistGammaNDaughtersAcc->Fill(nPartDaughters);
2e12941b 1982
112b1d0c 1983 iPartDaughter = part->GetLastDaughter();
1984 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
1985 {
1986 partDaughter = stack->Particle(iPartDaughter);
1987 if (partDaughter)
1988 {
1989 pdgDaugther = partDaughter->GetPDG(0);
1990 if (pdgDaugther) {
1991 fHistGammaDaughtersAcc->Fill(pdgDaugther->PdgCode());
1992 fHistGammaDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
2e12941b 1993
112b1d0c 1994 if (stack->IsPhysicalPrimary(iPart))
1995 {
1996 if (IsGammaConversion(stack, part, pdg))
1997 {
1998 fHistConvGammaDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
1999 }
2000 else
2001 {
2002 fHistNonConvGammaDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
2003 }
2004 }
f2adae2b 2005
112b1d0c 2006 }
2007 }
2008 }
2009 }
2010 }
2011 } // end of if gamma
2012
2013 // Neutral particles
2014 if (TMath::Abs(pdg->Charge() - fCuts->GetMonteCarloNeutralParticle()) <1e-3 )
2015 {
2016 //fNeutralMultiplicity++;
2017 fTotNeutralEt += et;
2e12941b 2018
112b1d0c 2019 // inside EMCal acceptance
2020 //if (IsInAcceptance(part,pdg))
2021 if (IsInAcceptance(aliPart))
2022 {
2023 fTotNeutralEtAcc += et;
2024 //fTotEtAcc += et;
2025 }
2026 } // end of neutral particles block
2027 //Charged particles
2028 else if (TMath::Abs( pdg->Charge() - fCuts->GetMonteCarloNeutralParticle())>1e-3 )
2029 {
2030 //fChargedMultiplicity++;
2031 fTotChargedEt += et;
2e12941b 2032
112b1d0c 2033 // inside EMCal acceptance
2034 //if (IsInAcceptance(part,pdg,extParam))
2035 if (IsInAcceptance(aliPart))
2036 {
2037 fTotChargedEtAcc += et;
2038 //fTotEtAcc += et;
2039 } // inside EMCal acceptance
2e12941b 2040
112b1d0c 2041 //if (TrackHitsCalo(extParam)) // magnetic field info not filled?
2042 //{
2043 // if (pdg->Charge() > 0) fHistPhivsPtPos->Fill(part->Phi(),part->Pt());
2044 // else if (pdg->Charge() < 0) fHistPhivsPtNeg->Fill(part->Phi(), part->Pt());
2045 //}
2046 } // end of charged particles block
2047 } // end of is primary
2048 else // not a primary
2049 {
2050 if (pdgMom)
2051 {
2052 if (pdgMom->PdgCode() == fgK0SCode)
2053 {
2054 fHistK0EtaEET->Fill(part->Energy(),part->Eta(),et);
2055 fHistK0EtaPtET->Fill(part->Pt(),part->Eta(),et);
2056 fHistK0EtaET->Fill(et,part->Eta());
2057 fHistK0EtaE->Fill(part->Energy(),part->Eta());
2058 fHistK0EtaPt->Fill(part->Pt(),part->Eta());
2059 fK0totET += et;
2060 }
f2adae2b 2061
112b1d0c 2062 if (pdgMom->PdgCode() == fgLambdaCode || pdgMom->PdgCode() == fgAntiLambdaCode)
2063 {
2064 fHistLambdaEtaEET->Fill(part->Energy(),part->Eta(),et);
2065 fHistLambdaEtaPtET->Fill(part->Pt(),part->Eta(),et);
2066 fHistLambdaEtaET->Fill(et,part->Eta());
2067 fHistLambdaEtaE->Fill(part->Energy(),part->Eta());
2068 fHistLambdaEtaPt->Fill(part->Pt(),part->Eta());
2069 fLambdatotET += et;
2070 }
2071 }
2072
2073 if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())) // exclude secondaries from interactions inside the EMCal
2074 {
2075 if (pdg->PdgCode() == fgEPlusCode || pdg->PdgCode() == fgEMinusCode)
2076 {
2077 fHistNPPElectronEtaEET->Fill(part->Energy(),part->Eta(),et);
2078 fHistNPPElectronEtaPtET->Fill(part->Pt(),part->Eta(),et);
2079 fHistNPPElectronEtaET->Fill(et,part->Eta());
2080 fHistNPPElectronEtaE->Fill(part->Energy(),part->Eta());
2081 fHistNPPElectronEtaPt->Fill(part->Pt(),part->Eta());
2082 fNPPElectrontotET += et;
2e12941b 2083
112b1d0c 2084 // inside EMCal acceptance
2085 //if (IsInAcceptance(part,pdg,extParam))
2086 if (IsInAcceptance(aliPart))
2087 {
2088 fHistNPPElectronAccEtaEET->Fill(part->Energy(),part->Eta(),et);
2089 fHistNPPElectronAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
2090 fHistNPPElectronAccEtaE->Fill(part->Energy(),part->Eta());
2091 fHistNPPElectronAccEtaPt->Fill(part->Pt(),part->Eta());
2092 }
2e12941b 2093
112b1d0c 2094 // few checks
77063fe4 2095 if(pdgMom) fHistNPPElectronFirstMother->Fill(pdgMom->PdgCode());
112b1d0c 2096 fHistNPPElectronFirstMotherXY->Fill(part->Vx(),part->Vy());
2097 fHistNPPElectronNDaughters->Fill(nPartDaughters);
2e12941b 2098
112b1d0c 2099 iPartDaughter = part->GetLastDaughter();
2100 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2101 {
2102 partDaughter = stack->Particle(iPartDaughter);
2103 if (partDaughter)
2104 {
2105 pdgDaugther = partDaughter->GetPDG(0);
2106 if (pdgDaugther) {
2107 fHistNPPElectronDaughters->Fill(pdgDaugther->PdgCode());
2108 fHistNPPElectronDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
2109 }
2110 }
2111 }
2e12941b 2112
112b1d0c 2113 // inside EMCal acceptance
2114 //if (IsInAcceptance(part,pdg,extParam))
2115 if (IsInAcceptance(aliPart))
2116 {
77063fe4 2117 if(pdgMom) fHistNPPElectronFirstMotherAcc->Fill(pdgMom->PdgCode());
112b1d0c 2118 fHistNPPElectronFirstMotherXYAcc->Fill(part->Vx(),part->Vy());
2119 fHistNPPElectronNDaughtersAcc->Fill(nPartDaughters);
2120
2121 iPartDaughter = part->GetLastDaughter();
2122 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2123 {
2124 partDaughter = stack->Particle(iPartDaughter);
2125 if (partDaughter)
2126 {
2127 pdgDaugther = partDaughter->GetPDG(0);
2128 if (pdgDaugther) {
2129 fHistNPPElectronDaughtersAcc->Fill(pdgDaugther->PdgCode());
2130 fHistNPPElectronDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
2131 }
2132 }
2133 }
2134 }
2e12941b 2135
112b1d0c 2136 } // end of if electron
2e12941b 2137
112b1d0c 2138 if(pdg->PdgCode() == fgGammaCode)
2139 {
2140 fHistNPPGammaEtaEET->Fill(part->Energy(),part->Eta(),et);
2141 fHistNPPGammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
2142 fHistNPPGammaEtaET->Fill(et,part->Eta());
2143 fHistNPPGammaEtaE->Fill(part->Energy(),part->Eta());
2144 fHistNPPGammaEtaPt->Fill(part->Pt(),part->Eta());
2145 fNPPGammatotET += et;
2e12941b 2146
112b1d0c 2147 if (pdgMom)
2148 {
2149 if (pdgMom->PdgCode() == fgPi0Code)
2150 {
2151 fHistNPPPi0GammaEtaEET->Fill(part->Energy(),part->Eta(),et);
2152 fHistNPPPi0GammaEtaPtET->Fill(part->Pt(),part->Eta(),et);
2153 fHistNPPPi0GammaEtaET->Fill(et,part->Eta());
2154 fHistNPPPi0GammaEtaE->Fill(part->Energy(),part->Eta());
2155 fHistNPPPi0GammaEtaPt->Fill(part->Pt(),part->Eta());
2156 fNPPPi0GammatotET += et;
2157 }
2158 }
2e12941b 2159
112b1d0c 2160 // inside EMCal acceptance
2161 //if (IsInAcceptance(part,pdg))
2162 if (IsInAcceptance(aliPart))
2163 {
2164 fHistNPPGammaAccEtaEET->Fill(part->Energy(),part->Eta(),et);
2165 fHistNPPGammaAccEtaPtET->Fill(part->Pt(),part->Eta(),et);
2166 fHistNPPGammaAccEtaE->Fill(part->Energy(),part->Eta());
2167 fHistNPPGammaAccEtaPt->Fill(part->Pt(),part->Eta());
2168 }
2e12941b 2169
112b1d0c 2170 // few checks
2171 if (pdgMom)
2172 fHistNPPGammaFirstMother->Fill(pdgMom->PdgCode());
2173 fHistNPPGammaFirstMotherXY->Fill(part->Vx(),part->Vy());
2174 fHistNPPGammaNDaughters->Fill(nPartDaughters);
2e12941b 2175
112b1d0c 2176 iPartDaughter = part->GetLastDaughter();
2177 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2178 {
2179 partDaughter = stack->Particle(iPartDaughter);
2180 if (partDaughter)
2181 {
2182 pdgDaugther = partDaughter->GetPDG(0);
2183 if (pdgDaugther) {
2184 fHistNPPGammaDaughters->Fill(pdgDaugther->PdgCode());
2185 fHistNPPGammaDaughtersXY->Fill(partDaughter->Vx(),partDaughter->Vy());
2186 }
2187 }
2188 }
2e12941b 2189
112b1d0c 2190 // inside EMCal acceptance
2191 //if (IsInAcceptance(part,pdg))
2192 if (IsInAcceptance(aliPart))
2193 {
2194 fHistNPPGammaFirstMotherAcc->Fill(pdgMom->PdgCode());
2195 fHistNPPGammaFirstMotherXYAcc->Fill(part->Vx(),part->Vy());
2196 fHistNPPGammaNDaughtersAcc->Fill(nPartDaughters);
2197
2198 iPartDaughter = part->GetLastDaughter();
2199 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2200 {
2201 partDaughter = stack->Particle(iPartDaughter);
2202 if (partDaughter)
2203 {
2204 pdgDaugther = partDaughter->GetPDG(0);
2205 if (pdgDaugther) {
2206 fHistNPPGammaDaughtersAcc->Fill(pdgDaugther->PdgCode());
2207 fHistNPPGammaDaughtersXYAcc->Fill(partDaughter->Vx(),partDaughter->Vy());
2208 }
2209 }
2210 }
2211 }
f2adae2b 2212
112b1d0c 2213 } // end of gamma
2214 }
2215 } // end of NOT a primary
2216 } // end of eta cut (Inside ALICE central barrel acceptance)
2e12941b 2217
112b1d0c 2218 if (extParam)
2219 delete extParam;
2e12941b 2220
112b1d0c 2221 }// end of loop over TParticles
2222 fTotEt = fTotChargedEt + fTotNeutralEt;
2223 fTotEtAcc = fTotChargedEtAcc + fTotNeutralEtAcc;
2224
2225 fTotElectrontotET = fElectrontotET + fConvElectrontotET + fScatElectrontotET;
2226 fTotElectronAcctotET = fElectronAcctotET + fConvElectronAcctotET + fScatElectronAcctotET;
2227 fTotGammatotET = fGammatotET + fAnnihGammatotET + fScatGammatotET;
2228 fTotGammaAcctotET = fGammaAcctotET + fAnnihGammaAcctotET + fScatGammaAcctotET;
2229 fTotEMtotET = fTotElectrontotET + fTotGammatotET;
2230 fTotEMAcctotET = fTotElectronAcctotET + fTotGammaAcctotET;
2231 fTotNPPEMtotET = fNPPElectrontotET + fNPPGammatotET;
2232 fTotChargedtotET = fMuontotET + fPiontotET + fKaontotET + fProtontotET;
2233 fTotChargedAcctotET = fMuonAcctotET + fPionAcctotET + fKaonAcctotET + fProtonAcctotET;
2234 fTotNeutraltotET = fNeutrontotET + fK0totET + fLambdatotET;
2235 fTotaltotET = fTotEMtotET + fTotNPPEMtotET + fTotChargedtotET + fTotNeutraltotET;
2236 fTotalAcctotET = fTotEMAcctotET + fTotChargedAcctotET;
2237
2238 //FillHistograms();
2239
2240 fHistPrimtotET->Fill(fPrimtotET);
2241 fHistPrimAcctotET->Fill(fPrimAcctotET);
2242
2243 fHistElectrontotET->Fill(fElectrontotET);
2244 fHistElectronAcctotET->Fill(fElectronAcctotET);
2245 fHistConvElectrontotET->Fill(fConvElectrontotET);
2246 fHistConvElectronAcctotET->Fill(fConvElectronAcctotET);
2247 fHistScatElectrontotET->Fill(fScatElectrontotET);
2248 fHistScatElectronAcctotET->Fill(fScatElectronAcctotET);
2249
2250 fHistTotElectrontotET->Fill(fTotElectrontotET);
2251 fHistTotElectronAcctotET->Fill(fTotElectronAcctotET);
2252
2253 fHistGammatotET->Fill(fGammatotET);
2254 fHistGammaAcctotET->Fill(fGammaAcctotET);
2255 fHistAnnihGammatotET->Fill(fAnnihGammatotET);
2256 fHistAnnihGammaAcctotET->Fill(fAnnihGammaAcctotET);
2257 fHistScatGammatotET->Fill(fScatGammatotET);
2258 fHistScatGammaAcctotET->Fill(fScatGammaAcctotET);
2259
2260 fHistTotGammatotET->Fill(fTotGammatotET);
2261 fHistTotGammaAcctotET->Fill(fTotGammaAcctotET);
2262
2263 fHistTotEMtotET->Fill(fTotEMtotET);
2264 fHistTotEMAcctotET->Fill(fTotEMAcctotET);
2265
2266 fHistConvGammatotET->Fill(fConvGammatotET);
2267 fHistNonConvGammatotET->Fill(fNonConvGammatotET);
2268 fHistConvGammaAcctotET->Fill(fConvGammaAcctotET);
2269 fHistNonConvGammaAcctotET->Fill(fNonConvGammaAcctotET);
2270
2271 fHistNPPElectrontotET->Fill(fNPPElectrontotET);
2272 fHistNPPGammatotET->Fill(fNPPGammatotET);
2273
2274 fHistTotNPPEMtotET->Fill(fTotNPPEMtotET);
2275
2276 fHistNPPPi0GammatotET->Fill(fNPPPi0GammatotET);
2277
2278 fHistMuontotET->Fill(fMuontotET);
2279 fHistMuonAcctotET->Fill(fMuonAcctotET);
2280 fHistPiontotET->Fill(fPiontotET);
2281 fHistPionAcctotET->Fill(fPionAcctotET);
2282 fHistKaontotET->Fill(fKaontotET);
2283 fHistKaonAcctotET->Fill(fKaonAcctotET);
2284 fHistProtontotET->Fill(fProtontotET);
2285 fHistProtonAcctotET->Fill(fProtonAcctotET);
2286
2287 fHistTotChargedtotET->Fill(fTotChargedtotET);
2288 fHistTotChargedAcctotET->Fill(fTotChargedAcctotET);
2289
2290 fHistNeutrontotET->Fill(fNeutrontotET);
2291 fHistNeutronAcctotET->Fill(fNeutronAcctotET);
2292 fHistK0totET->Fill(fK0totET);
2293 fHistLambdatotET->Fill(fNeutrontotET);
2294
2295 fHistTotNeutraltotET->Fill(fTotNeutraltotET);
2296
2297 fHistTotaltotET->Fill(fTotaltotET);
2298 fHistTotalAcctotET->Fill(fTotalAcctotET);
2299
2300 return 0;
2e12941b 2301}
2302
2303Int_t AliAnalysisEmEtMonteCarlo::AnalyseEvent(AliVEvent* ev,AliVEvent* ev2)
2304{ // analyse MC and real event info
112b1d0c 2305 if(!ev || !ev2){//Marcelo - should use AliError
2306 Printf("ERROR: Event does not exist");
2307 return 0;
2308 }
2309
2310 AliMCEvent *mcEvent = dynamic_cast<AliMCEvent*>(ev);
2311 AliESDEvent *realEvent = dynamic_cast<AliESDEvent*>(ev2);
77063fe4 2312 if(!mcEvent || !realEvent){//Marcelo - should use AliError
2313 Printf("ERROR: Event does not exist");
2314 return 0;
2315 }
112b1d0c 2316 fCentBin= -1;
2317 if(fDataSet==20100){//If this is Pb+Pb
2318 AliCentrality *centrality = realEvent->GetCentrality();
2319 if(fNCentBins<21) fCentBin= centrality->GetCentralityClass10(fCentralityMethod);
2320 else{ fCentBin= centrality->GetCentralityClass5(fCentralityMethod);}
2321 }
2322
2323 if(!fGeoUt){
2324 fGeoUt = AliEMCALGeometry::GetInstance("EMCAL_FIRSTYEARV1");//new AliEMCALGeometry("EMCAL_FIRSTYEAR","EMCAL");
2325 AliInfo("Creating new AliEMCALGeometry");
2326 }
2327 //fGeoUt = new AliEMCALGeometry("EMCAL_FIRSTYEAR","EMCAL");
2328 //fGeoUt->SetMisalMatrix(realEvent->GetEMCALMatrix(0),0);
2329
2330 ResetEventValues();
2331 AnalyseEvent(ev);
2332
2333 AliStack *stack = mcEvent->Stack();
2334 if (!stack)
f2adae2b 2335 {
112b1d0c 2336 Printf("ERROR: Could not get stack");
2337 return 0;
f2adae2b 2338 }
2e12941b 2339
112b1d0c 2340 Int_t nStackTracks = stack->GetNtrack();
2e12941b 2341
112b1d0c 2342 // get all emcal clusters
2343 TRefArray* caloClusters = new TRefArray();
2344 realEvent->GetEMCALClusters( caloClusters );
2e12941b 2345
112b1d0c 2346 Int_t nCluster = caloClusters->GetEntries();
2e12941b 2347
112b1d0c 2348 Float_t pos[3] = {0};
2349 TVector3 caloPos(0,0,0);
2350 TVector3 trackPos(0,0,0);
2e12941b 2351
112b1d0c 2352 // loop the clusters
2353 for (int iCluster = 0; iCluster < nCluster; iCluster++ )
f2adae2b 2354 {
112b1d0c 2355 AliESDCaloCluster* caloCluster = ( AliESDCaloCluster* )caloClusters->At( iCluster );
2356 Float_t caloE = caloCluster->E();
2357 caloCluster->GetPosition(pos);
2358 caloPos.SetXYZ(pos[0],pos[1],pos[2]);
2e12941b 2359
112b1d0c 2360 UInt_t iPart = (UInt_t)TMath::Abs(caloCluster->GetLabel());
2361 TParticle *part = stack->Particle(iPart);
2e12941b 2362
112b1d0c 2363 if (!part)
f2adae2b 2364 {//Marcelo -- use AliError
112b1d0c 2365 Printf("No MC particle %d", iCluster);
2366 continue;
2e12941b 2367 }
2368
112b1d0c 2369 TParticlePDG *pdg = part->GetPDG(0);
2e12941b 2370
112b1d0c 2371 TParticle *partMom = 0;
2372 TParticlePDG *pdgMom = 0;
2e12941b 2373
112b1d0c 2374 Int_t nPartDaughters = part->GetNDaughters();
2375 TParticle *partDaughter = 0;
2376 TParticlePDG *pdgDaugther = 0;
2e12941b 2377
112b1d0c 2378 if (!pdg)
f2adae2b 2379 {//Marcelo -- use AliError
112b1d0c 2380 Printf("ERROR-2: Could not get particle PDG %d", iPart);
2381 continue;
2e12941b 2382 }
2383
112b1d0c 2384 Int_t iPartMom = part->GetMother(0);
2385 Int_t iPartDaughter = 0;
2e12941b 2386
112b1d0c 2387 if ((iPartMom>=0) && (iPartMom < nStackTracks))
2388 {
2389 partMom = stack->Particle(iPartMom);
2390 pdgMom = partMom->GetPDG(0);
2391 }
2e12941b 2392
112b1d0c 2393 // find the track associated to this MC particle
2394 TObjArray* list = fEsdtrackCutsITSTPC->GetAcceptedTracks(realEvent);
2395 Int_t nGoodTracks = list->GetEntries();
2396 Bool_t trackFound = kFALSE;
2397 Bool_t trackProjected = kFALSE;
2398 Float_t res = 0, resMin=9999;
2399 AliESDtrack *track = 0;
2400 AliEMCALTrack *emcTrack = 0;
2401 AliEMCALTrack *emcTrack2 = 0;
2402 AliExternalTrackParam* extParamTPart = 0;
2403 AliESDtrack *esdTPart = 0;
2404 AliEMCALTrack *emcTPart = 0;
2405
2406 // find corresponding track
2407 for (Int_t iTrack = 0; iTrack < nGoodTracks; iTrack++)
2408 {
2409 track = dynamic_cast<AliESDtrack*> (list->At(iTrack));
2410 if (!track)
2411 {//Marcelo -use AliError
2412 Printf("ERROR: Could not get track %d", iTrack);
2413 continue;
2414 }
2415 else
2416 {
2417 UInt_t label = (UInt_t)TMath::Abs(track->GetLabel());
2418 if (label == iPart)
2419 {
2420 trackFound = kTRUE;
2421 emcTrack = new AliEMCALTrack(*track);
2422 fHistAllPtRecPtMC->Fill(part->Pt(),track->Pt());
2e12941b 2423
112b1d0c 2424 if (GetTrackProjection(emcTrack,trackPos,caloPos))
2425 {
2426 trackProjected = kTRUE;
2427 res = sqrt(pow(trackPos.Phi()-caloPos.Phi(),2)+pow(trackPos.Eta()-caloPos.Eta(),2));
2428 }
2429 else
2430 res = -1.;
2e12941b 2431
112b1d0c 2432 break;
2433 }
2434 }
2435 }
2436
2437 if (!trackFound)
2438 {
2439 track = 0;
2440 emcTrack = 0;
2441 res = -2.;
2442 }
2443
2444 //create an external track param for projection
2445 extParamTPart = CreateExternalTrackParam(part);
2e12941b 2446
112b1d0c 2447 // create esd and emcal tracks out of TParticle (used for projection)
2448 esdTPart = new AliESDtrack(part);
2e12941b 2449
112b1d0c 2450 if (esdTPart && extParamTPart)
2451 {
2452 esdTPart->SetOuterParam(extParamTPart,0);
2453 emcTPart = new AliEMCALTrack(*esdTPart);
2454 }
2e12941b 2455
112b1d0c 2456 // few checks
2457 // compare MC and Rec energies for all particles
2458 fHistAllERecEMC->Fill(part->Energy(),caloE);
2459 //Marcelo - doesn't it make sense to change this so that we use the tracks matched by the official code?
2460 if (TMath::Abs( pdg->Charge() - fCuts->GetMonteCarloNeutralParticle()) > 1e-3)
2461 {
2462 //Printf("calo.Phi = %f, calo.Eta = %f \n", caloPos.Phi(), caloPos.Eta());
2e12941b 2463
112b1d0c 2464 if (trackProjected)
2465 {
2466 //Printf("good track.Phi = %f, track.Eta = %f \n", trackPos.Phi(), trackPos.Eta());
2467 fHistChargedRes->Fill(trackPos.Phi()-caloPos.Phi(),trackPos.Eta()-caloPos.Eta());
2468 }
2e12941b 2469
112b1d0c 2470 if (GetTrackProjection(emcTPart,trackPos,caloPos))
2471 {
2472 fHistChargedRes2->Fill(trackPos.Phi()-caloPos.Phi(),trackPos.Eta()-caloPos.Eta());
2473 //Printf("track.Phi = %f, track.Eta = %f \n", trackPos.Phi(), trackPos.Eta());
2474 }
2e12941b 2475
112b1d0c 2476 if (GetTrackProjection(extParamTPart,trackPos))
2477 {
2478 fHistChargedRes3->Fill(trackPos.Phi()-caloPos.Phi(),trackPos.Eta()-caloPos.Eta());
2479 //Printf("track.Phi = %f, track.Eta = %f \n", trackPos.Phi(), trackPos.Eta());
2480 }
2e12941b 2481
112b1d0c 2482 }
2483 else if (TMath::Abs(pdg->Charge() - fCuts->GetMonteCarloNeutralParticle()) < 1e-3 )
2484 {
2485 fHistNeutralRes->Fill(part->Phi()-caloPos.Phi(),part->Eta()-caloPos.Eta());
2486 }
2e12941b 2487
112b1d0c 2488 if(pdg->PdgCode() == fgGammaCode)
2489 {
2490 // compare MC and Rec energies for gammas
2491 fHistGammaERecEMC->Fill(part->Energy(),caloE);
2492 fHistGammaRes->Fill(part->Phi()-caloPos.Phi(),part->Eta()-caloPos.Eta());
2493 }
2e12941b 2494
112b1d0c 2495 if (pdg->PdgCode() == fgEPlusCode || pdg->PdgCode() == fgEMinusCode)
2496 {
2497 // compare MC and Rec energies for electrons
2498 fHistElectronERecEMC->Fill(part->Energy(),caloE);
2499 if (GetTrackProjection(extParamTPart,trackPos))
2500 {
2501 fHistElectronRes->Fill(trackPos.Phi()-caloPos.Phi(),trackPos.Eta()-caloPos.Eta());
2502 }
2503 }
2e12941b 2504
112b1d0c 2505 // calculate ET
2506 Double_t et = CalcET(part,pdg);
2507 Double_t etDep = CalcETDep(caloE,part,pdg);
2e12941b 2508
112b1d0c 2509 // Check if it is a primary particle
2510 if (IsPrimary(stack,iPart,pdg,iPartMom,pdgMom))
2511 {
2512 if (TMath::Abs(TMath::Abs(pdg->Charge()) - fCuts->GetMonteCarloSingleChargedParticle())<1e-3 && TMath::Abs(TMath::Abs(pdg->Charge()) - fCuts->GetMonteCarloNeutralParticle())<1e-3) continue;
2513
2514 if (stack->IsPhysicalPrimary(iPart))
2515 {
2516 fHistPrimRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2517 fHistPrimRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2518 fHistPrimRecEtaET->Fill(et,part->Eta());
2519 fPrimRectotET += et;
2520
2521 fHistPrimRecEtaEDepETDep->Fill(part->Energy(),part->Eta(),etDep);
2522 fHistPrimRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2523 fHistPrimRecEtaETDep->Fill(etDep,part->Eta());
2524 fPrimRectotETDep += etDep;
2525 }
2526
2527 if(pdg->PdgCode() == fgGammaCode)
2528 {
2529 if (stack->IsPhysicalPrimary(iPart))
2530 {
2531 fHistGammaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2532 fHistGammaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2533 fHistGammaRecEtaET->Fill(et,part->Eta());
2534 fHistGammaRecEtaE->Fill(part->Energy(),part->Eta());
2535 fHistGammaRecEtaPt->Fill(part->Pt(),part->Eta());
2536
2537 resMin=9999;
2538 for (Int_t iTrack = 0; iTrack < nGoodTracks; iTrack++)
2539 {
2540 track = dynamic_cast<AliESDtrack*> (list->At(iTrack));
2541 if (!track)
2542 {//Marcelo -use AliError
2543 Printf("ERROR: Could not get track %d", iTrack);
2544 continue;
2545 }
2546 else
2547 {
2548 emcTrack2 = new AliEMCALTrack(*track);
2549
2550 if (GetTrackProjection(emcTrack2,trackPos,caloPos))
2551 {
2552 res = sqrt(pow(trackPos.Phi()-caloPos.Phi(),2)+pow(trackPos.Eta()-caloPos.Eta(),2));
2553 }
2554 if (res<resMin)
2555 resMin=res;
2556
2557 delete emcTrack2;
2558 emcTrack2=0;
2559 }
2560 }
2561
2562 fHistGammaRecResEET->Fill(part->Energy(),resMin,et);// Add Caio
2563 fHistGammaRecResPtET->Fill(part->Pt(),resMin,et);// Add Caio
2564 fHistGammaRecResE->Fill(part->Energy(),resMin);// Add Caio
2565 fHistGammaRecResPt->Fill(part->Pt(),resMin);// Add Caio
2566 fHistGammaRecResEDepETDep->Fill(caloE,resMin,etDep);// Add Caio
2567 fHistGammaRecResPtETDep->Fill(part->Pt(),resMin,etDep);// Add Caio
2568
2569
2570 fGammaRectotET += et;
2571 }
2572 else if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz()))
2573 {
2574 if (IsMotherPrimaryElectron(stack,iPartMom,pdgMom))
2575 {
2576 fHistAnnihGammaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2577 fHistAnnihGammaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2578 fHistAnnihGammaRecEtaET->Fill(et,part->Eta());
2579 fHistAnnihGammaRecEtaE->Fill(part->Energy(),part->Eta());
2580 fHistAnnihGammaRecEtaPt->Fill(part->Pt(),part->Eta());
2581 fAnnihGammaRectotET += et;
2582 }
2583 else if (IsMotherPrimaryGamma(stack,iPartMom,pdgMom))
2584 {
2585 fHistScatGammaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2586 fHistScatGammaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2587 fHistScatGammaRecEtaET->Fill(et,part->Eta());
2588 fHistScatGammaRecEtaE->Fill(part->Energy(),part->Eta());
2589 fHistScatGammaRecEtaPt->Fill(part->Pt(),part->Eta());
2590 fScatGammaRectotET += et;
2591 }
2592 }
2e12941b 2593
112b1d0c 2594 // few checks
2595 if (pdgMom)
2596 fHistGammaFirstMotherRec->Fill(pdgMom->PdgCode());
2597 fHistGammaFirstMotherXYRec->Fill(part->Vx(),part->Vy());
2598 fHistGammaNDaughtersRec->Fill(nPartDaughters);
2e12941b 2599
112b1d0c 2600 iPartDaughter = part->GetLastDaughter();
2601 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2602 {
2603 partDaughter = stack->Particle(iPartDaughter);
2604 if (partDaughter)
2605 {
2606 pdgDaugther = partDaughter->GetPDG(0);
2607 if (pdgDaugther) {
2608 fHistGammaDaughtersRec->Fill(pdgDaugther->PdgCode());
2609 fHistGammaDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2e12941b 2610
112b1d0c 2611 if (stack->IsPhysicalPrimary(iPart))
2612 {
2613 if (IsGammaConversion(stack, part, pdg))
2614 {
2615 fHistConvGammaDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2616 }
2617 else
2618 {
2619 fHistNonConvGammaDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2620 }
2621 }
2622 }
2623 }
2624 }
2625 } // gamma
2e12941b 2626
112b1d0c 2627 if (pdg->PdgCode() == fgEPlusCode || pdg->PdgCode() == fgEMinusCode)
2628 {
2629 if (stack->IsPhysicalPrimary(iPart))
2630 {
2631 fHistElectronRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2632 fHistElectronRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2633 fHistElectronRecEtaET->Fill(et,part->Eta());
2634 fHistElectronRecEtaE->Fill(part->Energy(),part->Eta());
2635 fHistElectronRecEtaPt->Fill(part->Pt(),part->Eta());
2636 fElectronRectotET += et;
2637 }
2638 else if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz()))
2639 {
2640 if (IsMotherPrimaryGamma(stack,iPartMom,pdgMom))
2641 {
2642 fHistConvElectronRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2643 fHistConvElectronRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2644 fHistConvElectronRecEtaET->Fill(et,part->Eta());
2645 fHistConvElectronRecEtaE->Fill(part->Energy(),part->Eta());
2646 fHistConvElectronRecEtaPt->Fill(part->Pt(),part->Eta());
2647 fConvElectronRectotET += et;
2648 }
2649 else if (IsMotherPrimaryElectron(stack,iPartMom,pdgMom))
2650 {
2651 fHistScatElectronRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2652 fHistScatElectronRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2653 fHistScatElectronRecEtaET->Fill(et,part->Eta());
2654 fHistScatElectronRecEtaE->Fill(part->Energy(),part->Eta());
2655 fHistScatElectronRecEtaPt->Fill(part->Pt(),part->Eta());
2656 fScatElectronRectotET += et;
2657 }
2658 }
2e12941b 2659
112b1d0c 2660 // few checks
2661 if (pdgMom)
2662 fHistElectronFirstMotherRec->Fill(pdgMom->PdgCode());
2663 fHistElectronFirstMotherXYRec->Fill(part->Vx(),part->Vy());
2664 fHistElectronNDaughtersRec->Fill(nPartDaughters);
2e12941b 2665
112b1d0c 2666 iPartDaughter = part->GetLastDaughter();
2667 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2668 {
2669 partDaughter = stack->Particle(iPartDaughter);
2670 if (partDaughter)
2671 {
2672 pdgDaugther = partDaughter->GetPDG(0);
2673 if (pdgDaugther) {
2674 fHistElectronDaughtersRec->Fill(pdgDaugther->PdgCode());
2675 fHistElectronDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2676 }
2677 }
2678 }
2679 } // electrons
2e12941b 2680
112b1d0c 2681 if (pdg->PdgCode() == fgMuPlusCode || pdg->PdgCode() == fgMuMinusCode)
2682 {
2683 fHistMuonRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2684 fHistMuonRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2685 fHistMuonRecEtaET->Fill(et,part->Eta());
2686 fHistMuonRecEtaE->Fill(part->Energy(),part->Eta());
2687 fHistMuonRecEtaPt->Fill(part->Pt(),part->Eta());
2688 fMuonRectotET += et;
2e12941b 2689
112b1d0c 2690 fHistMuonRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2691 fHistMuonRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2692 fHistMuonRecEtaETDep->Fill(etDep,part->Eta());
2693 fMuonRectotETDep += etDep;
2694
2695 if (trackProjected)
2696 {
2697 fHistMuonRecResEET->Fill(part->Energy(),res,et);
2698 fHistMuonRecResPtET->Fill(part->Pt(),res,et);
2699 fHistMuonRecResE->Fill(part->Energy(),res);
2700 fHistMuonRecResPt->Fill(part->Pt(),res);
2701 fHistMuonRecResEDepETDep->Fill(caloE,res,etDep);
2702 fHistMuonRecResPtETDep->Fill(part->Pt(),res,etDep);
2e12941b 2703
112b1d0c 2704 if ((res>0.) && (res<fResCut))
2705 {
2706 fHistMuonMatchEtaEET->Fill(part->Energy(),part->Eta(),et);
2707 fHistMuonMatchEtaPtET->Fill(part->Pt(),part->Eta(),et);
2708 fHistMuonMatchEtaET->Fill(et,part->Eta());
2709 fHistMuonMatchEtaE->Fill(part->Energy(),part->Eta());
2710 fHistMuonMatchEtaPt->Fill(part->Pt(),part->Eta());
2711 fMuonMatchtotET += et;
2e12941b 2712
112b1d0c 2713 fHistMuonMatchEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2714 fHistMuonMatchEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2715 fHistMuonMatchEtaETDep->Fill(etDep,part->Eta());
2716 fMuonMatchtotETDep += etDep;
2717 }
2718 }
2719 }
2e12941b 2720
112b1d0c 2721 if (pdg->PdgCode() == fgPiPlusCode || pdg->PdgCode() == fgPiMinusCode)
2722 {
2723 fHistPionRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2724 fHistPionRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2725 fHistPionRecEtaET->Fill(et,part->Eta());
2726 fHistPionRecEtaE->Fill(part->Energy(),part->Eta());
2727 fHistPionRecEtaPt->Fill(part->Pt(),part->Eta());
2728 fPionRectotET += et;
2729
2730 fHistPionRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2731 fHistPionRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2732 fHistPionRecEtaETDep->Fill(etDep,part->Eta());
2733 fPionRectotETDep += etDep;
2734
2735 if (trackProjected)
2736 {
2737 fHistPionRecResEET->Fill(part->Energy(),res,et);
2738 fHistPionRecResPtET->Fill(part->Pt(),res,et);
2739 fHistPionRecResE->Fill(part->Energy(),res);
2740 fHistPionRecResPt->Fill(part->Pt(),res);
2741 fHistPionRecResEDepETDep->Fill(caloE,res,etDep);
2742 fHistPionRecResPtETDep->Fill(part->Pt(),res,etDep);
2743
2744 if ((res>0.) && (res<fResCut))
2745 {
2746 fHistPionMatchEtaEET->Fill(part->Energy(),part->Eta(),et);
2747 fHistPionMatchEtaPtET->Fill(part->Pt(),part->Eta(),et);
2748 fHistPionMatchEtaET->Fill(et,part->Eta());
2749 fHistPionMatchEtaE->Fill(part->Energy(),part->Eta());
2750 fHistPionMatchEtaPt->Fill(part->Pt(),part->Eta());
2751 fPionMatchtotET += et;
2e12941b 2752
112b1d0c 2753 fHistPionMatchEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2754 fHistPionMatchEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2755 fHistPionMatchEtaETDep->Fill(etDep,part->Eta());
2756 fPionMatchtotETDep += etDep;
2757 }
2758 }
2759 }
2e12941b 2760
112b1d0c 2761 if (pdg->PdgCode() == fgKPlusCode || pdg->PdgCode() == fgKMinusCode)
2762 {
2763 fHistKaonRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2764 fHistKaonRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2765 fHistKaonRecEtaET->Fill(et,part->Eta());
2766 fHistKaonRecEtaE->Fill(part->Energy(),part->Eta());
2767 fHistKaonRecEtaPt->Fill(part->Pt(),part->Eta());
2768 fKaonRectotET += et;
2769
2770 fHistKaonRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2771 fHistKaonRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2772 fHistKaonRecEtaETDep->Fill(etDep,part->Eta());
2773 fKaonRectotETDep += etDep;
2774
2775 if (trackProjected)
2776 {
2777 fHistKaonRecResEET->Fill(part->Energy(),res,et);
2778 fHistKaonRecResPtET->Fill(part->Pt(),res,et);
2779 fHistKaonRecResE->Fill(part->Energy(),res);
2780 fHistKaonRecResPt->Fill(part->Pt(),res);
2781
2782 fHistKaonRecResEDepETDep->Fill(caloE,res,etDep);
2783 fHistKaonRecResPtETDep->Fill(part->Pt(),res,etDep);
2784
2785 if ((res>0.) && (res<fResCut))
2786 {
2787 fHistKaonMatchEtaEET->Fill(part->Energy(),part->Eta(),et);
2788 fHistKaonMatchEtaPtET->Fill(part->Pt(),part->Eta(),et);
2789 fHistKaonMatchEtaET->Fill(et,part->Eta());
2790 fHistKaonMatchEtaE->Fill(part->Energy(),part->Eta());
2791 fHistKaonMatchEtaPt->Fill(part->Pt(),part->Eta());
2792 fKaonMatchtotET += et;
2e12941b 2793
112b1d0c 2794 fHistKaonMatchEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2795 fHistKaonMatchEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2796 fHistKaonMatchEtaETDep->Fill(etDep,part->Eta());
2797 fKaonMatchtotETDep += etDep;
2798 }
2799 }
2800 }
2801
2802 if (pdg->PdgCode() == fgProtonCode || pdg->PdgCode() == fgAntiProtonCode)
2803 {
2804 fHistProtonRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2805 fHistProtonRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2806 fHistProtonRecEtaET->Fill(et,part->Eta());
2807 fHistProtonRecEtaE->Fill(part->Energy(),part->Eta());
2808 fHistProtonRecEtaPt->Fill(part->Pt(),part->Eta());
2809 fProtonRectotET += et;
2810
2811 fHistProtonRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2812 fHistProtonRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2813 fHistProtonRecEtaETDep->Fill(etDep,part->Eta());
2814 fProtonRectotETDep += etDep;
2815
2816 if (trackProjected)
2817 {
2818 fHistProtonRecResEET->Fill(part->Energy(),res,et);
2819 fHistProtonRecResPtET->Fill(part->Pt(),res,et);
2820 fHistProtonRecResE->Fill(part->Energy(),res);
2821 fHistProtonRecResPt->Fill(part->Pt(),res);
2822
2823 fHistProtonRecResEDepETDep->Fill(caloE,res,etDep);
2824 fHistProtonRecResPtETDep->Fill(part->Pt(),res,etDep);
2825
2826 if ((res>0.) && (res<fResCut))
2827 {
2828 fHistProtonMatchEtaEET->Fill(part->Energy(),part->Eta(),et);
2829 fHistProtonMatchEtaPtET->Fill(part->Pt(),part->Eta(),et);
2830 fHistProtonMatchEtaET->Fill(et,part->Eta());
2831 fHistProtonMatchEtaE->Fill(part->Energy(),part->Eta());
2832 fHistProtonMatchEtaPt->Fill(part->Pt(),part->Eta());
2833 fProtonMatchtotET += et;
2e12941b 2834
112b1d0c 2835 fHistProtonMatchEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2836 fHistProtonMatchEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2837 fHistProtonMatchEtaETDep->Fill(etDep,part->Eta());
2838 fProtonMatchtotETDep += etDep;
2839 }
2840 }
2841 }
2e12941b 2842
112b1d0c 2843 if (pdg->PdgCode() == fgNeutronCode || pdg->PdgCode() == fgAntiNeutronCode)
2844 {
2845 fHistNeutronRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2846 fHistNeutronRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2847 fHistNeutronRecEtaET->Fill(et,part->Eta());
2848 fHistNeutronRecEtaE->Fill(part->Energy(),part->Eta());
2849 fHistNeutronRecEtaPt->Fill(part->Pt(),part->Eta());
2850 fNeutronRectotET += et;
2851
2852 fHistNeutronRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2853 fHistNeutronRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2854 fHistNeutronRecEtaETDep->Fill(etDep,part->Eta());
2855 fNeutronRectotETDep += etDep;
2856 }
2e12941b 2857
112b1d0c 2858 if (emcTrack)
2859 delete emcTrack;
2860 if (esdTPart)
2861 delete esdTPart;
2862 if (emcTPart)
2863 delete emcTPart;
2864 if (extParamTPart)
2865 delete extParamTPart;
2866 } // end of primary tracks
2867 else // not a primary
2868 {
2869 if (pdgMom)
2870 {
2871 if (pdgMom->PdgCode() == fgK0SCode)
2872 {
2873 fHistK0RecEtaEET->Fill(part->Energy(),part->Eta(),et);
2874 fHistK0RecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2875 fHistK0RecEtaET->Fill(et,part->Eta());
2876 fHistK0RecEtaE->Fill(part->Energy(),part->Eta());
2877 fHistK0RecEtaPt->Fill(part->Pt(),part->Eta());
2878 fK0RectotET += et;
2e12941b 2879
112b1d0c 2880 fHistK0RecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2881 fHistK0RecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2882 fHistK0RecEtaETDep->Fill(etDep,part->Eta());
2883 fK0RectotETDep += etDep;
2884 }
2e12941b 2885
112b1d0c 2886 if (pdgMom->PdgCode() == fgLambdaCode || pdgMom->PdgCode() == fgAntiLambdaCode)
2887 {
2888 fHistLambdaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2889 fHistLambdaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2890 fHistLambdaRecEtaET->Fill(et,part->Eta());
2891 fHistLambdaRecEtaE->Fill(part->Energy(),part->Eta());
2892 fHistLambdaRecEtaPt->Fill(part->Pt(),part->Eta());
2893 fLambdaRectotET += et;
2894
2895 fHistLambdaRecEtaEDepETDep->Fill(caloE,part->Eta(),etDep);
2896 fHistLambdaRecEtaPtETDep->Fill(part->Pt(),part->Eta(),etDep);
2897 fHistLambdaRecEtaETDep->Fill(etDep,part->Eta());
2898 fLambdaRectotETDep += etDep;
2899 }
2900 }
2e12941b 2901
112b1d0c 2902 if (!fGeoUt->IsInEMCAL(part->Vx(),part->Vy(),part->Vz())) // exclude secondaries from interactions inside the EMCal
2903 {
2904 if (pdg->PdgCode() == fgEPlusCode || pdg->PdgCode() == fgEMinusCode)
2905 {
2906 fHistNPPElectronRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2907 fHistNPPElectronRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2908 fHistNPPElectronRecEtaET->Fill(et,part->Eta());
2909 fHistNPPElectronRecEtaE->Fill(part->Energy(),part->Eta());
2910 fHistNPPElectronRecEtaPt->Fill(part->Pt(),part->Eta());
2911 fNPPElectronRectotET += et;
2e12941b 2912
112b1d0c 2913 // few checks
77063fe4 2914 if(pdgMom) fHistNPPElectronFirstMotherRec->Fill(pdgMom->PdgCode());
112b1d0c 2915 fHistNPPElectronFirstMotherXYRec->Fill(part->Vx(),part->Vy());
2916 fHistNPPElectronNDaughtersRec->Fill(nPartDaughters);
2e12941b 2917
112b1d0c 2918 iPartDaughter = part->GetLastDaughter();
2919 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2920 {
2921 partDaughter = stack->Particle(iPartDaughter);
2922 if (partDaughter)
2923 {
2924 pdgDaugther = partDaughter->GetPDG(0);
2925 if (pdgDaugther) {
2926 fHistNPPElectronDaughtersRec->Fill(pdgDaugther->PdgCode());
2927 fHistNPPElectronDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2928 }
2929 }
2930 }
2931 } // end of if electron
2e12941b 2932
112b1d0c 2933 if(pdg->PdgCode() == fgGammaCode)
2934 {
2935 fHistNPPGammaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2936 fHistNPPGammaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2937 fHistNPPGammaRecEtaET->Fill(et,part->Eta());
2938 fHistNPPGammaRecEtaE->Fill(part->Energy(),part->Eta());
2939 fHistNPPGammaRecEtaPt->Fill(part->Pt(),part->Eta());
2940 fNPPGammaRectotET += et;
2e12941b 2941
112b1d0c 2942 if (pdgMom)
2943 {
2944 if (pdgMom->PdgCode() == fgPi0Code)
2945 {
2946 fHistNPPPi0GammaRecEtaEET->Fill(part->Energy(),part->Eta(),et);
2947 fHistNPPPi0GammaRecEtaPtET->Fill(part->Pt(),part->Eta(),et);
2948 fHistNPPPi0GammaRecEtaET->Fill(et,part->Eta());
2949 fHistNPPPi0GammaRecEtaE->Fill(part->Energy(),part->Eta());
2950 fHistNPPPi0GammaRecEtaPt->Fill(part->Pt(),part->Eta());
2951 fNPPPi0GammaRectotET += et;
2952 }
2953 }
2e12941b 2954
112b1d0c 2955 // few checks
2956 if (pdgMom)
2957 fHistNPPGammaFirstMotherRec->Fill(pdgMom->PdgCode());
2958 fHistNPPGammaFirstMotherXYRec->Fill(part->Vx(),part->Vy());
2959 fHistNPPGammaNDaughtersRec->Fill(nPartDaughters);
2e12941b 2960
112b1d0c 2961 iPartDaughter = part->GetLastDaughter();
2962 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
2963 {
2964 partDaughter = stack->Particle(iPartDaughter);
2965 if (partDaughter)
2966 {
2967 pdgDaugther = partDaughter->GetPDG(0);
2968 if (pdgDaugther) {
2969 fHistNPPGammaDaughtersRec->Fill(pdgDaugther->PdgCode());
2970 fHistNPPGammaDaughtersXYRec->Fill(partDaughter->Vx(),partDaughter->Vy());
2971 }
2972 }
2973 }
2974 } // end of gamma
2975 }
2976 } // end of NOT a primary
f2adae2b 2977 } // end of loop over clusters
112b1d0c 2978
2979 fTotElectronRectotET = fElectronRectotET + fConvElectronRectotET + fScatElectronRectotET;
2980 fTotGammaRectotET = fGammaRectotET + fAnnihGammaRectotET + fScatElectronRectotET;
2981 fTotEMRectotET = fTotElectronRectotET + fTotGammaRectotET;
2982 fTotNPPEMRectotET = fNPPElectronRectotET + fNPPGammaRectotET;
2983 fTotChargedRectotET = fMuonRectotET + fPionRectotET + fKaonRectotET + fProtonRectotET;
2984 fTotChargedRectotETDep = fMuonRectotETDep + fPionRectotETDep + fKaonRectotETDep + fProtonRectotETDep;
2985 fTotChargedMatchtotET = fMuonMatchtotET + fPionMatchtotET + fKaonMatchtotET + fProtonMatchtotET;
2986 fTotChargedMatchtotETDep = fMuonMatchtotETDep + fPionMatchtotETDep + fKaonMatchtotETDep + fProtonMatchtotETDep;
2987 fTotNeutralRectotET = fNeutronRectotET + fK0RectotET + fLambdaRectotET;
2988 fTotNeutralRectotETDep = fNeutronRectotETDep + fK0RectotETDep + fLambdaRectotETDep;
2989 fTotalRectotET = fTotEMRectotET + fTotNPPEMRectotET + fTotChargedRectotET + fTotNeutralRectotET;
2990 fTotalRectotETDep = fTotEMRectotET + fTotNPPEMRectotET + fTotChargedRectotETDep + fTotNeutralRectotETDep;
2991
2992 fHistPrimRectotET->Fill(fPrimRectotET);
2993 fHistPrimRectotETDep->Fill(fPrimRectotETDep);
2994
2995 fHistElectronRectotET->Fill(fElectronRectotET);
2996 fHistConvElectronRectotET->Fill(fConvElectronRectotET);
2997 fHistScatElectronRectotET->Fill(fScatElectronRectotET);
2998 fHistTotElectronRectotET->Fill(fTotElectronRectotET);
2999
3000 fHistGammaRectotET->Fill(fGammaRectotET);
3001 fHistAnnihGammaRectotET->Fill(fAnnihGammaRectotET);
3002 fHistScatGammaRectotET->Fill(fScatGammaRectotET);
3003 fHistTotGammaRectotET->Fill(fTotGammaRectotET);
3004
3005 fHistTotEMRectotET->Fill(fTotEMRectotET);
3006
3007 fHistNPPElectronRectotET->Fill(fNPPElectronRectotET);
3008 fHistNPPGammaRectotET->Fill(fNPPGammaRectotET);
3009 fHistTotNPPEMRectotET->Fill(fTotNPPEMRectotET);
3010
3011 fHistNPPPi0GammaRectotET->Fill(fNPPPi0GammaRectotET);
3012
3013 fHistMuonRectotET->Fill(fMuonRectotET);
3014 fHistMuonRectotETDep->Fill(fMuonRectotETDep);
3015 fHistMuonMatchtotET->Fill(fMuonMatchtotET);
3016 fHistMuonMatchtotETDep->Fill(fMuonMatchtotETDep);
3017 fHistPionRectotET->Fill(fPionRectotET);
3018 fHistPionRectotETDep->Fill(fPionRectotETDep);
3019 fHistPionMatchtotET->Fill(fPionMatchtotET);
3020 fHistPionMatchtotETDep->Fill(fPionMatchtotETDep);
3021 fHistKaonRectotET->Fill(fKaonRectotET);
3022 fHistKaonRectotETDep->Fill(fKaonRectotETDep);
3023 fHistKaonMatchtotET->Fill(fKaonMatchtotET);
3024 fHistKaonMatchtotETDep->Fill(fKaonMatchtotETDep);
3025 fHistProtonRectotET->Fill(fProtonRectotET);
3026 fHistProtonRectotETDep->Fill(fProtonRectotETDep);
3027 fHistProtonMatchtotET->Fill(fProtonMatchtotET);
3028 fHistProtonMatchtotETDep->Fill(fProtonMatchtotETDep);
3029 fHistTotChargedRectotET->Fill(fTotChargedRectotET);
3030 fHistTotChargedRectotETDep->Fill(fTotChargedRectotETDep);
3031 fHistTotChargedMatchtotET->Fill(fTotChargedMatchtotET);
3032 fHistTotChargedMatchtotETDep->Fill(fTotChargedMatchtotETDep);
3033
3034 fHistNeutronRectotET->Fill(fNeutronRectotET);
3035 fHistNeutronRectotETDep->Fill(fNeutronRectotETDep);
3036 fHistK0RectotET->Fill(fK0RectotET);
3037 fHistK0RectotETDep->Fill(fK0RectotETDep);
3038 fHistLambdaRectotET->Fill(fLambdaRectotET);
3039 fHistLambdaRectotETDep->Fill(fLambdaRectotETDep);
3040 fHistTotNeutralRectotET->Fill(fTotNeutralRectotET);
3041 fHistTotNeutralRectotETDep->Fill(fTotNeutralRectotETDep);
3042
3043 fHistTotalRectotET->Fill(fTotalRectotET);
3044 fHistTotalRectotETDep->Fill(fTotalRectotETDep);
3045
3046 //delete fGeoUt;
3047 delete caloClusters;//Marcelo - Christine - make as object & don't create new one each event
3048
3049 return 0;
2e12941b 3050}
3051
3052void AliAnalysisEmEtMonteCarlo::Init()
3053{ // init
112b1d0c 3054 AliAnalysisEt::Init();
2e12941b 3055
112b1d0c 3056 fDetectorRadius = fCuts->GetGeometryEmcalDetectorRadius();
3057 fSingleCellEnergyCut = fCuts->GetReconstructedEmcalSingleCellEnergyCut();
2e12941b 3058}
3059
3060void AliAnalysisEmEtMonteCarlo::ResetEventValues()
3061{ // reset event values
112b1d0c 3062 AliAnalysisEt::ResetEventValues();
3063
3064 fPrimtotET = 0; fPrimAcctotET = 0; fPrimRectotET = 0; fPrimRectotETDep = 0;
3065
3066 fElectrontotET = 0; fElectronAcctotET = 0; fElectronRectotET = 0;
3067 fConvElectrontotET = 0; fConvElectronAcctotET = 0; fConvElectronRectotET = 0; fScatElectrontotET = 0; fScatElectronAcctotET = 0; fScatElectronRectotET = 0;
3068 fTotElectrontotET = 0, fTotElectronAcctotET = 0, fTotElectronRectotET = 0;
3069
3070 fGammatotET = 0; fGammaAcctotET = 0; fGammaRectotET = 0;
3071 fAnnihGammatotET = 0; fAnnihGammaAcctotET = 0; fAnnihGammaRectotET = 0; fScatGammatotET = 0; fScatGammaAcctotET = 0; fScatGammaRectotET = 0;
3072 fTotGammatotET = 0, fTotGammaAcctotET = 0, fTotGammaRectotET = 0;
3073
3074 fTotEMtotET = 0, fTotEMAcctotET = 0, fTotEMRectotET = 0;
3075
3076 fConvGammatotET = 0; fNonConvGammatotET = 0; fConvGammaAcctotET = 0; fNonConvGammaAcctotET = 0; fNPPPi0GammatotET = 0; fNPPPi0GammaRectotET = 0;
3077
3078 fNPPElectrontotET = 0; fNPPElectronRectotET = 0; fNPPGammatotET = 0; fNPPGammaRectotET = 0;
3079 fTotNPPEMtotET = 0, fTotNPPEMRectotET = 0;
3080
3081 fMuontotET = 0; fPiontotET = 0; fKaontotET = 0; fProtontotET = 0;
3082 fMuonAcctotET = 0; fPionAcctotET = 0; fKaonAcctotET = 0; fProtonAcctotET = 0;
3083 fMuonRectotET = 0; fMuonRectotETDep = 0; fPionRectotET = 0; fPionRectotETDep = 0; fKaonRectotET = 0; fKaonRectotETDep = 0; fProtonRectotET = 0; fProtonRectotETDep = 0;
3084 fMuonMatchtotET = 0; fMuonMatchtotETDep = 0; fPionMatchtotET = 0; fPionMatchtotETDep = 0; fKaonMatchtotET = 0; fKaonMatchtotETDep = 0; fProtonMatchtotET = 0; fProtonMatchtotETDep = 0;
3085 fTotChargedtotET = 0, fTotChargedAcctotET = 0, fTotChargedRectotET = 0, fTotChargedRectotETDep = 0, fTotChargedMatchtotET = 0, fTotChargedMatchtotETDep = 0;
3086
3087 fNeutrontotET = 0; fNeutronAcctotET = 0; fNeutronRectotET = 0; fNeutronRectotETDep = 0;
3088 fK0totET = 0; fK0RectotET = 0; fK0RectotETDep = 0; fLambdatotET = 0; fLambdaRectotET = 0; fLambdaRectotETDep = 0;
3089 fTotNeutraltotET = 0, fTotNeutralRectotET = 0, fTotNeutralRectotETDep = 0;
3090
3091 fTotaltotET = 0, fTotalAcctotET = 0, fTotalRectotET = 0, fTotalRectotETDep = 0;
2e12941b 3092}
3093
3094void AliAnalysisEmEtMonteCarlo::CreateHistograms()
112b1d0c 3095{ // histogram related Additions
3096 //AliAnalysisEt::CreateHistograms();
3097 fSelector = new AliAnalysisEtSelectorEmcal(fCuts);
3098 fHistPrimEtaEET = CreateEtaEHisto2D("fHistPrimEtaEET_","MC E_{T}, primary particles","E_{T}(GeV)");
3099 fHistPrimEtaPtET = CreateEtaPtHisto2D("fHistPrimEtaPtET_","MC E_{T}, primary particles","E_{T}(GeV)");
3100 fHistPrimEtaET = CreateEtaEtHisto2D("fHistPrimEtaET_","MC primary particles","#");
3101 TString histname = "fHistPrimtotET_" + fHistogramNameSuffix;
3102 fHistPrimtotET = new TH1F(histname.Data(),"total ET, primary particles",fgNumOfEBins, fgEAxis);
3103
3104 fHistPrimAccEtaEET = CreateEtaEHisto2D("fHistPrimAccEtaEET_","MC E_{T}, primary particles","E_{T}(GeV)");
3105 fHistPrimAccEtaPtET = CreateEtaPtHisto2D("fHistPrimAccEtaPtET_","MC E_{T}, primary particles","E_{T}(GeV)");
3106 fHistPrimAccEtaET = CreateEtaEtHisto2D("fHistPrimAccEtaET_","MC primary particles","#");
3107 histname = "fHistPrimAcctotET_" + fHistogramNameSuffix;
3108 fHistPrimAcctotET = new TH1F(histname.Data(),"total ET, primary particles",fgNumOfEBins, fgEAxis);
3109
3110 fHistPrimRecEtaEET = CreateEtaEHisto2D("fHistPrimRecEtaEET_","MC E_{T}, primary particles","E_{T}(GeV)");
3111 fHistPrimRecEtaPtET = CreateEtaPtHisto2D("fHistPrimRecEtaPtET_","MC E_{T}, primary particles","E_{T}(GeV)");
3112 fHistPrimRecEtaET = CreateEtaEtHisto2D("fHistPrimRecEtaET_","MC primary particles","#");
3113 histname = "fHistPrimRectotET_" + fHistogramNameSuffix;
3114 fHistPrimRectotET = new TH1F(histname.Data(),"total ET, primary particles",fgNumOfEBins, fgEAxis);
3115
3116 fHistPrimRecEtaEDepETDep = CreateEtaEHisto2D("fHistPrimRecEtaEDepETDep_","MC E_{T}, primary particles","E_{T}(GeV)");
3117 fHistPrimRecEtaPtETDep = CreateEtaPtHisto2D("fHistPrimRecEtaPtETDep_","MC E_{T}, primary particles","E_{T}(GeV)");
3118 fHistPrimRecEtaETDep = CreateEtaEtHisto2D("fHistPrimRecEtaETDep_","MC primary particles","#");
3119 histname = "fHistPrimRectotETDep_" + fHistogramNameSuffix;
3120 fHistPrimRectotETDep = new TH1F(histname.Data(),"total ET, primary particles",fgNumOfEBins, fgEAxis);
3121
3122 fHistElectronEtaEET = CreateEtaEHisto2D("fHistElectronEtaEET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3123 fHistElectronEtaPtET = CreateEtaPtHisto2D("fHistElectronEtaPtET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3124 fHistElectronEtaET = CreateEtaEtHisto2D("fHistElectronEtaET_","MC primary electrons","#");
3125 fHistElectronEtaE = CreateEtaEHisto2D("fHistElectronEtaE_","MC primary electrons","#");
3126 fHistElectronEtaPt = CreateEtaPtHisto2D("fHistElectronEtaPt_","MC primary electrons","#");
3127 histname = "fHistElectrontotET_" + fHistogramNameSuffix;
3128 fHistElectrontotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3129
3130 fHistConvElectronEtaEET = CreateEtaEHisto2D("fHistConvElectronEtaEET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3131 fHistConvElectronEtaPtET = CreateEtaPtHisto2D("fHistConvElectronEtaPtET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3132 fHistConvElectronEtaET = CreateEtaEtHisto2D("fHistConvElectronEtaET_","MC electrons from conversion","#");
3133 fHistConvElectronEtaE = CreateEtaEHisto2D("fHistConvElectronEtaE_","MC electrons from conversion","#");
3134 fHistConvElectronEtaPt = CreateEtaPtHisto2D("fHistConvElectronEtaPt_","MC electrons from conversion","#");
3135 histname = "fHistConvElectrontotET_" + fHistogramNameSuffix;
3136 fHistConvElectrontotET = new TH1F(histname.Data(),"total ET, MC electrons from conversion",fgNumOfEBins, fgEAxis);
3137
3138 fHistScatElectronEtaEET = CreateEtaEHisto2D("fHistScatElectronEtaEET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3139 fHistScatElectronEtaPtET = CreateEtaPtHisto2D("fHistScatElectronEtaPtET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3140 fHistScatElectronEtaET = CreateEtaEtHisto2D("fHistScatElectronEtaET_","MC electrons from Scattering","#");
3141 fHistScatElectronEtaE = CreateEtaEHisto2D("fHistScatElectronEtaE_","MC electrons from Scattering","#");
3142 fHistScatElectronEtaPt = CreateEtaPtHisto2D("fHistScatElectronEtaPt_","MC electrons from Scattering","#");
3143 histname = "fHistScatElectrontotET_" + fHistogramNameSuffix;
3144 fHistScatElectrontotET = new TH1F(histname.Data(),"total ET, MC electrons from Scattering",fgNumOfEBins, fgEAxis);
3145
3146 histname = "fHistTotElectrontotET_" + fHistogramNameSuffix;
3147 fHistTotElectrontotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3148
3149 fHistGammaEtaEET = CreateEtaEHisto2D("fHistGammaEtaEET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3150 fHistGammaEtaPtET = CreateEtaPtHisto2D("fHistGammaEtaPtET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3151 fHistGammaEtaET = CreateEtaEtHisto2D("fHistGammaEtaET_","MC primary gammas","#");
3152 fHistGammaEtaE = CreateEtaEHisto2D("fHistGammaEtaE_","MC primary gammas","#");
3153 fHistGammaEtaPt = CreateEtaPtHisto2D("fHistGammaEtaPt_","MC primary gammas","#");
3154 histname = "fHistGammatotET_" + fHistogramNameSuffix;
3155 fHistGammatotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3156
3157 fHistAnnihGammaEtaEET = CreateEtaEHisto2D("fHistAnnihGammaEtaEET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3158 fHistAnnihGammaEtaPtET = CreateEtaPtHisto2D("fHistAnnihGammaEtaPtET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3159 fHistAnnihGammaEtaET = CreateEtaEtHisto2D("fHistAnnihGammaEtaET_","MC Annihilation gammas","#");
3160 fHistAnnihGammaEtaE = CreateEtaEHisto2D("fHistAnnihGammaEtaE_","MC Annihilation gammas","#");
3161 fHistAnnihGammaEtaPt = CreateEtaPtHisto2D("fHistAnnihGammaEtaPt_","MC Annihilation gammas","#");
3162 histname = "fHistAnnihGammatotET_" + fHistogramNameSuffix;
3163 fHistAnnihGammatotET = new TH1F(histname.Data(),"total ET, MC Annihilation gammas",fgNumOfEBins, fgEAxis);
3164
3165 fHistScatGammaEtaEET = CreateEtaEHisto2D("fHistScatGammaEtaEET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3166 fHistScatGammaEtaPtET = CreateEtaPtHisto2D("fHistScatGammaEtaPtET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3167 fHistScatGammaEtaET = CreateEtaEtHisto2D("fHistScatGammaEtaET_","MC Scattering gammas","#");
3168 fHistScatGammaEtaE = CreateEtaEHisto2D("fHistScatGammaEtaE_","MC Scattering gammas","#");
3169 fHistScatGammaEtaPt = CreateEtaPtHisto2D("fHistScatGammaEtaPt_","MC Scattering gammas","#");
3170 histname = "fHistScatGammatotET_" + fHistogramNameSuffix;
3171 fHistScatGammatotET = new TH1F(histname.Data(),"total ET, MC Scattering gammas",fgNumOfEBins, fgEAxis);
3172
3173 fHistConvGammaEtaEET = CreateEtaEHisto2D("fHistConvGammaEtaEET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3174 fHistConvGammaEtaPtET = CreateEtaPtHisto2D("fHistConvGammaEtaPtET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3175 fHistConvGammaEtaET = CreateEtaEtHisto2D("fHistConvGammaEtaET_","MC non conversion primary gammas","#");
3176 fHistConvGammaEtaE = CreateEtaEHisto2D("fHistConvGammaEtaE_","MC non conversion primary gammas","#");
3177 fHistConvGammaEtaPt = CreateEtaPtHisto2D("fHistConvGammaEtaPt_","MC non conversion primary gammas","#");
3178 histname = "fHistConvGammatotET_" + fHistogramNameSuffix;
3179 fHistConvGammatotET = new TH1F(histname.Data(),"total ET, MC non conversion primary gammas",fgNumOfEBins, fgEAxis);
3180
3181 fHistNonConvGammaEtaEET = CreateEtaEHisto2D("fHistNonConvGammaEtaEET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3182 fHistNonConvGammaEtaPtET = CreateEtaPtHisto2D("fHistNonConvGammaEtaPtET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3183 fHistNonConvGammaEtaET = CreateEtaEtHisto2D("fHistNonConvGammaEtaET_","MC non conversion primary gammas","#");
3184 fHistNonConvGammaEtaE = CreateEtaEHisto2D("fHistNonConvGammaEtaE_","MC non conversion primary gammas","#");
3185 fHistNonConvGammaEtaPt = CreateEtaPtHisto2D("fHistNonConvGammaEtaPt_","MC non conversion primary gammas","#");
3186 histname = "fHistNonConvGammatotET_" + fHistogramNameSuffix;
3187 fHistNonConvGammatotET = new TH1F(histname.Data(),"total ET, MC non conversion primary gammas",fgNumOfEBins, fgEAxis);
3188
3189 histname = "fHistTotGammatotET_" + fHistogramNameSuffix;
3190 fHistTotGammatotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3191
3192 histname = "fHistTotEMtotET_" + fHistogramNameSuffix;
3193 fHistTotEMtotET = new TH1F(histname.Data(),"total electromagnetic ET",fgNumOfEBins, fgEAxis);
3194
3195 fHistNPPElectronEtaEET = CreateEtaEHisto2D("fHistNPPElectronEtaEET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3196 fHistNPPElectronEtaPtET = CreateEtaPtHisto2D("fHistNPPElectronEtaPtET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3197 fHistNPPElectronEtaET = CreateEtaEtHisto2D("fHistNPPElectronEtaET_","MC non-primary electrons","#");
3198 fHistNPPElectronEtaE = CreateEtaEHisto2D("fHistNPPElectronEtaE_","MC non-primary electrons","#");
3199 fHistNPPElectronEtaPt = CreateEtaPtHisto2D("fHistNPPElectronEtaPt_","MC non-primary electrons","#");
3200 histname = "fHistNPPElectrontotET_" + fHistogramNameSuffix;
3201 fHistNPPElectrontotET = new TH1F(histname.Data(),"total ET, MC non-primary electrons",fgNumOfEBins, fgEAxis);
3202
3203 fHistNPPGammaEtaEET = CreateEtaEHisto2D("fHistNPPGammaEtaEET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3204 fHistNPPGammaEtaPtET = CreateEtaPtHisto2D("fHistNPPGammaEtaPtET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3205 fHistNPPGammaEtaET = CreateEtaEtHisto2D("fHistNPPGammaEtaET_","MC non-primary gammas","#");
3206 fHistNPPGammaEtaE = CreateEtaEHisto2D("fHistNPPGammaEtaE_","MC non-primary gammas","#");
3207 fHistNPPGammaEtaPt = CreateEtaPtHisto2D("fHistNPPGammaEtaPt_","MC non-primary gammas","#");
3208 histname = "fHistNPPGammatotET_" + fHistogramNameSuffix;
3209 fHistNPPGammatotET = new TH1F(histname.Data(),"total ET, MC non-primary gammas",fgNumOfEBins, fgEAxis);
3210
3211 histname = "fHistTotNPPEMtotET_" + fHistogramNameSuffix;
3212 fHistTotNPPEMtotET = new TH1F(histname.Data(),"total ET, MC non-primary electromagnetic",fgNumOfEBins, fgEAxis);
3213
3214 fHistNPPPi0GammaEtaEET = CreateEtaEHisto2D("fHistNPPPi0GammaEtaEET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3215 fHistNPPPi0GammaEtaPtET = CreateEtaPtHisto2D("fHistNPPPi0GammaEtaPtET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3216 fHistNPPPi0GammaEtaET = CreateEtaEtHisto2D("fHistNPPPi0GammaEtaET_","MC non-primary gammas","#");
3217 fHistNPPPi0GammaEtaE = CreateEtaEHisto2D("fHistNPPPi0GammaEtaE_","MC non-primary gammas","#");
3218 fHistNPPPi0GammaEtaPt = CreateEtaPtHisto2D("fHistNPPPi0GammaEtaPt_","MC non-primary gammas","#");
3219 histname = "fHistNPPPi0GammatotET_" + fHistogramNameSuffix;
3220 fHistNPPPi0GammatotET = new TH1F(histname.Data(),"total ET, MC non-primary gammas",fgNumOfEBins, fgEAxis);
3221
3222 fHistElectronAccEtaEET = CreateEtaEHisto2D("fHistElectronAccEtaEET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3223 fHistElectronAccEtaPtET = CreateEtaPtHisto2D("fHistElectronAccEtaPtET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3224 fHistElectronAccEtaET = CreateEtaEtHisto2D("fHistElectronAccEtaET_","MC primary electrons","#");
3225 fHistElectronAccEtaE = CreateEtaEHisto2D("fHistElectronAccEtaE_","MC primary electrons","#");
3226 fHistElectronAccEtaPt = CreateEtaPtHisto2D("fHistElectronAccEtaPt_","MC primary electrons","#");
3227 histname = "fHistElectronAcctotET_" + fHistogramNameSuffix;
3228 fHistElectronAcctotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3229
3230 fHistConvElectronAccEtaEET = CreateEtaEHisto2D("fHistConvElectronAccEtaEET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3231 fHistConvElectronAccEtaPtET = CreateEtaPtHisto2D("fHistConvElectronAccEtaPtET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3232 fHistConvElectronAccEtaET = CreateEtaEtHisto2D("fHistConvElectronAccEtaET_","MC electrons from conversion","#");
3233 fHistConvElectronAccEtaE = CreateEtaEHisto2D("fHistConvElectronAccEtaE_","MC electrons from conversion","#");
3234 fHistConvElectronAccEtaPt = CreateEtaPtHisto2D("fHistConvElectronAccEtaPt_","MC electrons from conversion","#");
3235 histname = "fHistConvElectronAcctotET_" + fHistogramNameSuffix;
3236 fHistConvElectronAcctotET = new TH1F(histname.Data(),"total ET, MC electrons from conversion",fgNumOfEBins, fgEAxis);
3237
3238 fHistScatElectronAccEtaEET = CreateEtaEHisto2D("fHistScatElectronAccEtaEET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3239 fHistScatElectronAccEtaPtET = CreateEtaPtHisto2D("fHistScatElectronAccEtaPtET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3240 fHistScatElectronAccEtaET = CreateEtaEtHisto2D("fHistScatElectronAccEtaET_","MC electrons from Scattering","#");
3241 fHistScatElectronAccEtaE = CreateEtaEHisto2D("fHistScatElectronAccEtaE_","MC electrons from Scattering","#");
3242 fHistScatElectronAccEtaPt = CreateEtaPtHisto2D("fHistScatElectronAccEtaPt_","MC electrons from Scattering","#");
3243 histname = "fHistScatElectronAcctotET_" + fHistogramNameSuffix;
3244 fHistScatElectronAcctotET = new TH1F(histname.Data(),"total ET, MC electrons from Scattering",fgNumOfEBins, fgEAxis);
3245
3246 histname = "fHistTotElectronAcctotET_" + fHistogramNameSuffix;
3247 fHistTotElectronAcctotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3248
3249 fHistGammaAccEtaEET = CreateEtaEHisto2D("fHistGammaAccEtaEET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3250 fHistGammaAccEtaPtET = CreateEtaPtHisto2D("fHistGammaAccEtaPtET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3251 fHistGammaAccEtaET = CreateEtaEtHisto2D("fHistGammaAccEtaET_","MC primary gammas","#");
3252 fHistGammaAccEtaE = CreateEtaEHisto2D("fHistGammaAccEtaE_","MC primary gammas","#");
3253 fHistGammaAccEtaPt = CreateEtaPtHisto2D("fHistGammaAccEtaPt_","MC primary gammas","#");
3254 histname = "fHistGammaAcctotET_" + fHistogramNameSuffix;
3255 fHistGammaAcctotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3256
3257 fHistAnnihGammaAccEtaEET = CreateEtaEHisto2D("fHistAnnihGammaAccEtaEET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3258 fHistAnnihGammaAccEtaPtET = CreateEtaPtHisto2D("fHistAnnihGammaAccEtaPtET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3259 fHistAnnihGammaAccEtaET = CreateEtaEtHisto2D("fHistAnnihGammaAccEtaET_","MC Annihilation gammas","#");
3260 fHistAnnihGammaAccEtaE = CreateEtaEHisto2D("fHistAnnihGammaAccEtaE_","MC Annihilation gammas","#");
3261 fHistAnnihGammaAccEtaPt = CreateEtaPtHisto2D("fHistAnnihGammaAccEtaPt_","MC Annihilation gammas","#");
3262 histname = "fHistAnnihGammaAcctotET_" + fHistogramNameSuffix;
3263 fHistAnnihGammaAcctotET = new TH1F(histname.Data(),"total ET, MC Annihilation gammas",fgNumOfEBins, fgEAxis);
3264
3265 fHistScatGammaAccEtaEET = CreateEtaEHisto2D("fHistScatGammaAccEtaEET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3266 fHistScatGammaAccEtaPtET = CreateEtaPtHisto2D("fHistScatGammaAccEtaPtET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3267 fHistScatGammaAccEtaET = CreateEtaEtHisto2D("fHistScatGammaAccEtaET_","MC Scattering gammas","#");
3268 fHistScatGammaAccEtaE = CreateEtaEHisto2D("fHistScatGammaAccEtaE_","MC Scattering gammas","#");
3269 fHistScatGammaAccEtaPt = CreateEtaPtHisto2D("fHistScatGammaAccEtaPt_","MC Scattering gammas","#");
3270 histname = "fHistScatGammaAcctotET_" + fHistogramNameSuffix;
3271 fHistScatGammaAcctotET = new TH1F(histname.Data(),"total ET, MC Scattering gammas",fgNumOfEBins, fgEAxis);
3272
3273 fHistConvGammaAccEtaEET = CreateEtaEHisto2D("fHistConvGammaAccEtaEET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3274 fHistConvGammaAccEtaPtET = CreateEtaPtHisto2D("fHistConvGammaAccEtaPtET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3275 fHistConvGammaAccEtaET = CreateEtaEtHisto2D("fHistConvGammaAccEtaET_","MC non conversion primary gammas","#");
3276 fHistConvGammaAccEtaE = CreateEtaEHisto2D("fHistConvGammaAccEtaE_","MC non conversion primary gammas","#");
3277 fHistConvGammaAccEtaPt = CreateEtaPtHisto2D("fHistConvGammaAccEtaPt_","MC non conversion primary gammas","#");
3278 histname = "fHistConvGammaAcctotET_" + fHistogramNameSuffix;
3279 fHistConvGammaAcctotET = new TH1F(histname.Data(),"total ET, MC non conversion primary gammas",fgNumOfEBins, fgEAxis);
3280
3281 fHistNonConvGammaAccEtaEET = CreateEtaEHisto2D("fHistNonConvGammaAccEtaEET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3282 fHistNonConvGammaAccEtaPtET = CreateEtaPtHisto2D("fHistNonConvGammaAccEtaPtET_","MC E_{T}, non conversion primary gammas","E_{T}(GeV)");
3283 fHistNonConvGammaAccEtaET = CreateEtaEtHisto2D("fHistNonConvGammaAccEtaET_","MC non conversion primary gammas","#");
3284 fHistNonConvGammaAccEtaE = CreateEtaEHisto2D("fHistNonConvGammaAccEtaE_","MC non conversion primary gammas","#");
3285 fHistNonConvGammaAccEtaPt = CreateEtaPtHisto2D("fHistNonConvGammaAccEtaPt_","MC non conversion primary gammas","#");
3286 histname = "fHistNonConvGammaAcctotET_" + fHistogramNameSuffix;
3287 fHistNonConvGammaAcctotET = new TH1F(histname.Data(),"total ET, MC non conversion primary gammas",fgNumOfEBins, fgEAxis);
3288
3289 histname = "fHistTotGammaAcctotET_" + fHistogramNameSuffix;
3290 fHistTotGammaAcctotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3291
3292 histname = "fHistTotEMAcctotET_" + fHistogramNameSuffix;
3293 fHistTotEMAcctotET = new TH1F(histname.Data(),"total electromagnetic ET",fgNumOfEBins, fgEAxis);
3294
3295 fHistNPPElectronAccEtaEET = CreateEtaEHisto2D("fHistNPPElectronAccEtaEET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3296 fHistNPPElectronAccEtaPtET = CreateEtaPtHisto2D("fHistNPPElectronAccEtaPtET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3297 fHistNPPElectronAccEtaE = CreateEtaEHisto2D("fHistNPPElectronAccEtaE_","MC non-primary electrons","#");
3298 fHistNPPElectronAccEtaPt = CreateEtaPtHisto2D("fHistNPPElectronAccEtaPt_","MC non-primary electrons","#");
3299
3300 fHistNPPGammaAccEtaEET = CreateEtaEHisto2D("fHistNPPGammaAccEtaEET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3301 fHistNPPGammaAccEtaPtET = CreateEtaPtHisto2D("fHistNPPGammaAccEtaPtET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3302 fHistNPPGammaAccEtaE = CreateEtaEHisto2D("fHistNPPGammaAccEtaE_","MC non-primary gammas","#");
3303 fHistNPPGammaAccEtaPt = CreateEtaPtHisto2D("fHistNPPGammaAccEtaPt_","MC non-primary gammas","#");
3304
3305 fHistElectronRecEtaEET = CreateEtaEHisto2D("fHistElectronRecEtaEET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3306 fHistElectronRecEtaPtET = CreateEtaPtHisto2D("fHistElectronRecEtaPtET_","MC E_{T}, primary electrons","E_{T}(GeV)");
3307 fHistElectronRecEtaET = CreateEtaEtHisto2D("fHistElectronRecEtaET_","MC primary electrons","#");
3308 fHistElectronRecEtaE = CreateEtaEHisto2D("fHistElectronRecEtaE_","MC primary electrons","#");
3309 fHistElectronRecEtaPt = CreateEtaPtHisto2D("fHistElectronRecEtaPt_","MC primary electrons","#");
3310 histname = "fHistElectronRectotET_" + fHistogramNameSuffix;
3311 fHistElectronRectotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3312
3313 fHistConvElectronRecEtaEET = CreateEtaEHisto2D("fHistConvElectronRecEtaEET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3314 fHistConvElectronRecEtaPtET = CreateEtaPtHisto2D("fHistConvElectronRecEtaPtET_","MC E_{T}, electrons from conversion","E_{T}(GeV)");
3315 fHistConvElectronRecEtaET = CreateEtaEtHisto2D("fHistConvElectronRecEtaET_","MC electrons from conversion","#");
3316 fHistConvElectronRecEtaE = CreateEtaEHisto2D("fHistConvElectronRecEtaE_","MC electrons from conversion","#");
3317 fHistConvElectronRecEtaPt = CreateEtaPtHisto2D("fHistConvElectronRecEtaPt_","MC electrons from conversion","#");
3318 histname = "fHistConvElectronRectotET_" + fHistogramNameSuffix;
3319 fHistConvElectronRectotET = new TH1F(histname.Data(),"total ET, MC electrons from conversion",fgNumOfEBins, fgEAxis);
3320
3321 fHistScatElectronRecEtaEET = CreateEtaEHisto2D("fHistScatElectronRecEtaEET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3322 fHistScatElectronRecEtaPtET = CreateEtaPtHisto2D("fHistScatElectronRecEtaPtET_","MC E_{T}, electrons from Scattering","E_{T}(GeV)");
3323 fHistScatElectronRecEtaET = CreateEtaEtHisto2D("fHistScatElectronRecEtaET_","MC electrons from Scattering","#");
3324 fHistScatElectronRecEtaE = CreateEtaEHisto2D("fHistScatElectronRecEtaE_","MC electrons from Scattering","#");
3325 fHistScatElectronRecEtaPt = CreateEtaPtHisto2D("fHistScatElectronRecEtaPt_","MC electrons from Scattering","#");
3326 histname = "fHistScatElectronRectotET_" + fHistogramNameSuffix;
3327 fHistScatElectronRectotET = new TH1F(histname.Data(),"total ET, MC electrons from Scattering",fgNumOfEBins, fgEAxis);
3328
3329 histname = "fHistTotElectronRectotET_" + fHistogramNameSuffix;
3330 fHistTotElectronRectotET = new TH1F(histname.Data(),"total ET, MC primary electrons",fgNumOfEBins, fgEAxis);
3331
3332 fHistGammaRecEtaEET = CreateEtaEHisto2D("fHistGammaRecEtaEET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3333 fHistGammaRecEtaPtET = CreateEtaPtHisto2D("fHistGammaRecEtaPtET_","MC E_{T}, primary gammas","E_{T}(GeV)");
3334 fHistGammaRecEtaET = CreateEtaEtHisto2D("fHistGammaRecEtaET_","MC primary gammas","#");
3335 fHistGammaRecEtaE = CreateEtaEHisto2D("fHistGammaRecEtaE_","MC primary gammas","#");
3336 fHistGammaRecEtaPt = CreateEtaPtHisto2D("fHistGammaRecEtaPt_","MC primary gammas","#");
3337 histname = "fHistGammaRectotET_" + fHistogramNameSuffix;
3338 fHistGammaRectotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3339
3340 fHistGammaRecResEET = CreateResEHisto2D("fHistGammaRecResEET_","MC E_{T}, primary Gammas","E_{T}(GeV)");// Add Caio
3341 fHistGammaRecResPtET = CreateResPtHisto2D("fHistGammaRecResPtET_","MC E_{T}, primary Gammas","E_{T}(GeV)");// Add Caio
3342 fHistGammaRecResE = CreateResEHisto2D("fHistGammaRecResE_","MC primary Gammas","#");// Add Caio
3343 fHistGammaRecResPt = CreateResPtHisto2D("fHistGammaRecResPt_","MC primary Gammas","#");// Add Caio
3344 fHistGammaRecResEDepETDep = CreateResEHisto2D("fHistGammaRecResEDepETDep_","MC E_{T}, Gammas","E_{T} dep (GeV)"); // Add Caio
3345 fHistGammaRecResPtETDep = CreateResPtHisto2D("fHistGammaRecResPtETDep_","MC E_{T}, Gammas","E_{T} dep (GeV)");// Add Caio
3346
3347 fHistAnnihGammaRecEtaEET = CreateEtaEHisto2D("fHistAnnihGammaRecEtaEET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3348 fHistAnnihGammaRecEtaPtET = CreateEtaPtHisto2D("fHistAnnihGammaRecEtaPtET_","MC E_{T}, Annihilation gammas","E_{T}(GeV)");
3349 fHistAnnihGammaRecEtaET = CreateEtaEtHisto2D("fHistAnnihGammaRecEtaET_","MC Annihilation gammas","#");
3350 fHistAnnihGammaRecEtaE = CreateEtaEHisto2D("fHistAnnihGammaRecEtaE_","MC Annihilation gammas","#");
3351 fHistAnnihGammaRecEtaPt = CreateEtaPtHisto2D("fHistAnnihGammaRecEtaPt_","MC Annihilation gammas","#");
3352 histname = "fHistAnnihGammaRectotET_" + fHistogramNameSuffix;
3353 fHistAnnihGammaRectotET = new TH1F(histname.Data(),"total ET, MC Annihilation gammas",fgNumOfEBins, fgEAxis);
3354
3355 fHistScatGammaRecEtaEET = CreateEtaEHisto2D("fHistScatGammaRecEtaEET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3356 fHistScatGammaRecEtaPtET = CreateEtaPtHisto2D("fHistScatGammaRecEtaPtET_","MC E_{T}, Scattering gammas","E_{T}(GeV)");
3357 fHistScatGammaRecEtaET = CreateEtaEtHisto2D("fHistScatGammaRecEtaET_","MC Scattering gammas","#");
3358 fHistScatGammaRecEtaE = CreateEtaEHisto2D("fHistScatGammaRecEtaE_","MC Scattering gammas","#");
3359 fHistScatGammaRecEtaPt = CreateEtaPtHisto2D("fHistScatGammaRecEtaPt_","MC Scattering gammas","#");
3360 histname = "fHistScatGammaRectotET_" + fHistogramNameSuffix;
3361 fHistScatGammaRectotET = new TH1F(histname.Data(),"total ET, MC Scattering gammas",fgNumOfEBins, fgEAxis);
3362
3363 histname = "fHistTotGammaRectotET_" + fHistogramNameSuffix;
3364 fHistTotGammaRectotET = new TH1F(histname.Data(),"total ET, MC primary gammas",fgNumOfEBins, fgEAxis);
3365
3366 histname = "fHistTotEMRectotET_" + fHistogramNameSuffix;
3367 fHistTotEMRectotET = new TH1F(histname.Data(),"total electromagnetic ET",fgNumOfEBins, fgEAxis);
3368
3369 fHistNPPElectronRecEtaEET = CreateEtaEHisto2D("fHistNPPElectronRecEtaEET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3370 fHistNPPElectronRecEtaPtET = CreateEtaPtHisto2D("fHistNPPElectronRecEtaPtET_","MC E_{T}, non-primary electrons","E_{T}(GeV)");
3371 fHistNPPElectronRecEtaET = CreateEtaEtHisto2D("fHistNPPElectronRecEtaET_","MC non-primary electrons","#");
3372 fHistNPPElectronRecEtaE = CreateEtaEHisto2D("fHistNPPElectronRecEtaE_","MC non-primary electrons","#");
3373 fHistNPPElectronRecEtaPt = CreateEtaPtHisto2D("fHistNPPElectronRecEtaPt_","MC non-primary electrons","#");
3374 histname = "fHistNPPElectronRectotET_" + fHistogramNameSuffix;
3375 fHistNPPElectronRectotET = new TH1F(histname.Data(),"total ET, MC non-primary electrons",fgNumOfEBins, fgEAxis);
3376
3377 fHistNPPGammaRecEtaEET = CreateEtaEHisto2D("fHistNPPGammaRecEtaEET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3378 fHistNPPGammaRecEtaPtET = CreateEtaPtHisto2D("fHistNPPGammaRecEtaPtET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3379 fHistNPPGammaRecEtaET = CreateEtaEtHisto2D("fHistNPPGammaRecEtaET_","MC non-primary gammas","#");
3380 fHistNPPGammaRecEtaE = CreateEtaEHisto2D("fHistNPPGammaRecEtaE_","MC non-primary gammas","#");
3381 fHistNPPGammaRecEtaPt = CreateEtaPtHisto2D("fHistNPPGammaRecEtaPt_","MC non-primary gammas","#");
3382 histname = "fHistNPPGammaRectotET_" + fHistogramNameSuffix;
3383 fHistNPPGammaRectotET = new TH1F(histname.Data(),"total ET, MC non-primary gammas",fgNumOfEBins, fgEAxis);
3384
3385 histname = "fHistTotNPPEMRectotET_" + fHistogramNameSuffix;
3386 fHistTotNPPEMRectotET = new TH1F(histname.Data(),"total ET, MC non-primary electromagnetic",fgNumOfEBins, fgEAxis);
3387
3388 fHistNPPPi0GammaRecEtaEET = CreateEtaEHisto2D("fHistNPPPi0GammaRecEtaEET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3389 fHistNPPPi0GammaRecEtaPtET = CreateEtaPtHisto2D("fHistNPPPi0GammaRecEtaPtET_","MC E_{T}, non-primary gammas","E_{T}(GeV)");
3390 fHistNPPPi0GammaRecEtaET = CreateEtaEtHisto2D("fHistNPPPi0GammaRecEtaET_","MC non-primary gammas","#");
3391 fHistNPPPi0GammaRecEtaE = CreateEtaEHisto2D("fHistNPPPi0GammaRecEtaE_","MC non-primary gammas","#");
3392 fHistNPPPi0GammaRecEtaPt = CreateEtaPtHisto2D("fHistNPPPi0GammaRecEtaPt_","MC non-primary gammas","#");
3393 histname = "fHistNPPPi0GammaRectotET_" + fHistogramNameSuffix;
3394 fHistNPPPi0GammaRectotET = new TH1F(histname.Data(),"total ET, MC non-primary gammas",fgNumOfEBins, fgEAxis);
3395
3396 fHistMuonEtaEET = CreateEtaEHisto2D("fHistMuonEtaEET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3397 fHistMuonAccEtaEET = CreateEtaEHisto2D("fHistMuonAccEtaEET_","MC E_{T}, primary Muons, inside EMCal acceptance","E_{T}(GeV)");
3398 fHistMuonRecEtaEET = CreateEtaEHisto2D("fHistMuonRecEtaEET_","MC E_{T}, primary Muons, reconstructed","E_{T}(GeV)");
3399 fHistMuonMatchEtaEET = CreateEtaEHisto2D("fHistMuonMatchEtaEET_","MC E_{T}, primary Muons, tracking matched","E_{T}(GeV)");
3400
3401 fHistMuonEtaPtET = CreateEtaPtHisto2D("fHistMuonEtaPtET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3402 fHistMuonAccEtaPtET = CreateEtaPtHisto2D("fHistMuonAccEtaPtET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3403 fHistMuonRecEtaPtET = CreateEtaPtHisto2D("fHistMuonRecEtaPtET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3404 fHistMuonMatchEtaPtET = CreateEtaPtHisto2D("fHistMuonMatchEtaPtET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3405
3406 fHistMuonEtaET = CreateEtaEtHisto2D("fHistMuonEtaET_","MC primary Muons","#");
3407 fHistMuonAccEtaET = CreateEtaEtHisto2D("fHistMuonAccEtaET_","MC primary Muons","#");
3408 fHistMuonRecEtaET = CreateEtaEtHisto2D("fHistMuonRecEtaET_","MC primary Muons","#");
3409 fHistMuonMatchEtaET = CreateEtaEtHisto2D("fHistMuonMatchEtaET_","MC primary Muons","#");
3410
3411 fHistMuonEtaE = CreateEtaEHisto2D("fHistMuonEtaE_","MC primary Muons","#");
3412 fHistMuonAccEtaE = CreateEtaEHisto2D("fHistMuonAccEtaE_","MC primary Muons","#");
3413 fHistMuonRecEtaE = CreateEtaEHisto2D("fHistMuonRecEtaE_","MC primary Muons","#");
3414 fHistMuonMatchEtaE = CreateEtaEHisto2D("fHistMuonMatchEtaE_","MC primary Muons","#");
3415
3416 fHistMuonEtaPt = CreateEtaPtHisto2D("fHistMuonEtaPt_","MC primary Muons","#");
3417 fHistMuonAccEtaPt = CreateEtaPtHisto2D("fHistMuonAccEtaPt_","MC primary Muons","#");
3418 fHistMuonRecEtaPt = CreateEtaPtHisto2D("fHistMuonRecEtaPt_","MC primary Muons","#");
3419 fHistMuonMatchEtaPt = CreateEtaPtHisto2D("fHistMuonMatchEtaPt_","MC primary Muons","#");
3420
3421 histname = "fHistMuontotET_" + fHistogramNameSuffix;
3422 fHistMuontotET = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3423 histname = "fHistMuonAcctotET_" + fHistogramNameSuffix;
3424 fHistMuonAcctotET = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3425 histname = "fHistMuonRectotET_" + fHistogramNameSuffix;
3426 fHistMuonRectotET = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3427 histname = "fHistMuonMatchtotET_" + fHistogramNameSuffix;
3428 fHistMuonMatchtotET = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3429
3430 histname = "fHistMuonRectotETDep_" + fHistogramNameSuffix;
3431 fHistMuonRectotETDep = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3432 histname = "fHistMuonMatchtotETDep_" + fHistogramNameSuffix;
3433 fHistMuonMatchtotETDep = new TH1F(histname.Data(),"total ET, MC primary Muons",fgNumOfEBins, fgEAxis);
3434
3435 fHistMuonRecEtaEDepETDep = CreateEtaEHisto2D("fHistMuonRecEtaEDepETDep_","MC E_{T}, primary Muons, reconstructed","E_{T} dep (GeV)");
3436 fHistMuonMatchEtaEDepETDep = CreateEtaEHisto2D("fHistMuonMatchEtaEDepETDep_","MC E_{T}, primary Muons, tracking matched","E_{T} dep (GeV)");
3437
3438 fHistMuonRecEtaPtETDep = CreateEtaPtHisto2D("fHistMuonRecEtaPtETDep_","MC E_{T}, primary Muons","E_{T} dep (GeV)");
3439 fHistMuonMatchEtaPtETDep = CreateEtaPtHisto2D("fHistMuonMatchEtaPtETDep_","MC E_{T}, primary Muons","E_{T} dep(GeV)");
3440
3441 fHistMuonRecEtaETDep = CreateEtaEtHisto2D("fHistMuonRecEtaETDep_","MC primary Muons","#");
3442 fHistMuonMatchEtaETDep = CreateEtaEtHisto2D("fHistMuonMatchEtaETDep_","MC primary Muons","#");
3443
3444 fHistMuonRecResEET = CreateResEHisto2D("fHistMuonRecResEET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3445 fHistMuonRecResPtET = CreateResPtHisto2D("fHistMuonRecResPtET_","MC E_{T}, primary Muons","E_{T}(GeV)");
3446 fHistMuonRecResE = CreateResEHisto2D("fHistMuonRecResE_","MC primary Muons","#");
3447 fHistMuonRecResPt = CreateResPtHisto2D("fHistMuonRecResPt_","MC primary Muons","#");
3448 fHistMuonRecResEDepETDep = CreateResEHisto2D("fHistMuonRecResEDepETDep_","MC E_{T}, primary Muons","E_{T} dep (GeV)");
3449 fHistMuonRecResPtETDep = CreateResPtHisto2D("fHistMuonRecResPtETDep_","MC E_{T}, primary Muons","E_{T} dep (GeV)");
3450
3451 fHistPionEtaEET = CreateEtaEHisto2D("fHistPionEtaEET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3452 fHistPionAccEtaEET = CreateEtaEHisto2D("fHistPionAccEtaEET_","MC E_{T}, primary Pions, inside EMCal acceptance","E_{T}(GeV)");
3453 fHistPionRecEtaEET = CreateEtaEHisto2D("fHistPionRecEtaEET_","MC E_{T}, primary Pions, reconstructed","E_{T}(GeV)");
3454 fHistPionMatchEtaEET = CreateEtaEHisto2D("fHistPionMatchEtaEET_","MC E_{T}, primary Pions, tracking matched","E_{T}(GeV)");
3455
3456 fHistPionEtaPtET = CreateEtaPtHisto2D("fHistPionEtaPtET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3457 fHistPionAccEtaPtET = CreateEtaPtHisto2D("fHistPionAccEtaPtET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3458 fHistPionRecEtaPtET = CreateEtaPtHisto2D("fHistPionRecEtaPtET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3459 fHistPionMatchEtaPtET = CreateEtaPtHisto2D("fHistPionMatchEtaPtET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3460
3461 fHistPionEtaET = CreateEtaEtHisto2D("fHistPionEtaET_","MC primary Pions","#");
3462 fHistPionAccEtaET = CreateEtaEtHisto2D("fHistPionAccEtaET_","MC primary Pions","#");
3463 fHistPionRecEtaET = CreateEtaEtHisto2D("fHistPionRecEtaET_","MC primary Pions","#");
3464 fHistPionMatchEtaET = CreateEtaEtHisto2D("fHistPionMatchEtaET_","MC primary Pions","#");
3465
3466 fHistPionEtaE = CreateEtaEHisto2D("fHistPionEtaE_","MC primary Pions","#");
3467 fHistPionAccEtaE = CreateEtaEHisto2D("fHistPionAccEtaE_","MC primary Pions","#");
3468 fHistPionRecEtaE = CreateEtaEHisto2D("fHistPionRecEtaE_","MC primary Pions","#");
3469 fHistPionMatchEtaE = CreateEtaEHisto2D("fHistPionMatchEtaE_","MC primary Pions","#");
3470
3471 fHistPionEtaPt = CreateEtaPtHisto2D("fHistPionEtaPt_","MC primary Pions","#");
3472 fHistPionAccEtaPt = CreateEtaPtHisto2D("fHistPionAccEtaPt_","MC primary Pions","#");
3473 fHistPionRecEtaPt = CreateEtaPtHisto2D("fHistPionRecEtaPt_","MC primary Pions","#");
3474 fHistPionMatchEtaPt = CreateEtaPtHisto2D("fHistPionMatchEtaPt_","MC primary Pions","#");
3475
3476 histname = "fHistPiontotET_" + fHistogramNameSuffix;
3477 fHistPiontotET = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3478 histname = "fHistPionAcctotET_" + fHistogramNameSuffix;
3479 fHistPionAcctotET = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3480 histname = "fHistPionRectotET_" + fHistogramNameSuffix;
3481 fHistPionRectotET = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3482 histname = "fHistPionMatchtotET_" + fHistogramNameSuffix;
3483 fHistPionMatchtotET = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3484
3485 histname = "fHistPionRectotETDep_" + fHistogramNameSuffix;
3486 fHistPionRectotETDep = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3487 histname = "fHistPionMatchtotETDep_" + fHistogramNameSuffix;
3488 fHistPionMatchtotETDep = new TH1F(histname.Data(),"total ET, MC primary Pions",fgNumOfEBins, fgEAxis);
3489
3490 fHistPionRecEtaEDepETDep = CreateEtaEHisto2D("fHistPionRecEtaEDepETDep_","MC E_{T}, primary Pions, reconstructed","E_{T} dep (GeV)");
3491 fHistPionMatchEtaEDepETDep = CreateEtaEHisto2D("fHistPionMatchEtaEDepETDep_","MC E_{T}, primary Pions, tracking matched","E_{T} dep (GeV)");
3492
3493 fHistPionRecEtaPtETDep = CreateEtaPtHisto2D("fHistPionRecEtaPtETDep_","MC E_{T}, primary Pions","E_{T} dep (GeV)");
3494 fHistPionMatchEtaPtETDep = CreateEtaPtHisto2D("fHistPionMatchEtaPtETDep_","MC E_{T}, primary Pions","E_{T} dep(GeV)");
3495
3496 fHistPionRecEtaETDep = CreateEtaEtHisto2D("fHistPionRecEtaETDep_","MC primary Pions","#");
3497 fHistPionMatchEtaETDep = CreateEtaEtHisto2D("fHistPionMatchEtaETDep_","MC primary Pions","#");
3498
3499 fHistPionRecResEET = CreateResEHisto2D("fHistPionRecResEET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3500 fHistPionRecResPtET = CreateResPtHisto2D("fHistPionRecResPtET_","MC E_{T}, primary Pions","E_{T}(GeV)");
3501 fHistPionRecResE = CreateResEHisto2D("fHistPionRecResE_","MC primary Pions","#");
3502 fHistPionRecResPt = CreateResPtHisto2D("fHistPionRecResPt_","MC primary Pions","#");
3503 fHistPionRecResEDepETDep = CreateResEHisto2D("fHistPionRecResEDepETDep_","MC E_{T}, primary Pions","E_{T} dep (GeV)");
3504 fHistPionRecResPtETDep = CreateResPtHisto2D("fHistPionRecResPtETDep_","MC E_{T}, primary Pions","E_{T} dep (GeV)");
3505
3506 fHistKaonEtaEET = CreateEtaEHisto2D("fHistKaonEtaEET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3507 fHistKaonAccEtaEET = CreateEtaEHisto2D("fHistKaonAccEtaEET_","MC E_{T}, primary Kaons, inside EMCal acceptance","E_{T}(GeV)");
3508 fHistKaonRecEtaEET = CreateEtaEHisto2D("fHistKaonRecEtaEET_","MC E_{T}, primary Kaons, reconstructed","E_{T}(GeV)");
3509 fHistKaonMatchEtaEET = CreateEtaEHisto2D("fHistKaonMatchEtaEET_","MC E_{T}, primary Kaons, tracking matched","E_{T}(GeV)");
3510
3511 fHistKaonEtaPtET = CreateEtaPtHisto2D("fHistKaonEtaPtET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3512 fHistKaonAccEtaPtET = CreateEtaPtHisto2D("fHistKaonAccEtaPtET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3513 fHistKaonRecEtaPtET = CreateEtaPtHisto2D("fHistKaonRecEtaPtET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3514 fHistKaonMatchEtaPtET = CreateEtaPtHisto2D("fHistKaonMatchEtaPtET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3515
3516 fHistKaonEtaET = CreateEtaEtHisto2D("fHistKaonEtaET_","MC primary Kaons","#");
3517 fHistKaonAccEtaET = CreateEtaEtHisto2D("fHistKaonAccEtaET_","MC primary Kaons","#");
3518 fHistKaonRecEtaET = CreateEtaEtHisto2D("fHistKaonRecEtaET_","MC primary Kaons","#");
3519 fHistKaonMatchEtaET = CreateEtaEtHisto2D("fHistKaonMatchEtaET_","MC primary Kaons","#");
3520
3521 fHistKaonEtaE = CreateEtaEHisto2D("fHistKaonEtaE_","MC primary Kaons","#");
3522 fHistKaonAccEtaE = CreateEtaEHisto2D("fHistKaonAccEtaE_","MC primary Kaons","#");
3523 fHistKaonRecEtaE = CreateEtaEHisto2D("fHistKaonRecEtaE_","MC primary Kaons","#");
3524 fHistKaonMatchEtaE = CreateEtaEHisto2D("fHistKaonMatchEtaE_","MC primary Kaons","#");
3525
3526 fHistKaonEtaPt = CreateEtaPtHisto2D("fHistKaonEtaPt_","MC primary Kaons","#");
3527 fHistKaonAccEtaPt = CreateEtaPtHisto2D("fHistKaonAccEtaPt_","MC primary Kaons","#");
3528 fHistKaonRecEtaPt = CreateEtaPtHisto2D("fHistKaonRecEtaPt_","MC primary Kaons","#");
3529 fHistKaonMatchEtaPt = CreateEtaPtHisto2D("fHistKaonMatchEtaPt_","MC primary Kaons","#");
3530
3531 histname = "fHistKaontotET_" + fHistogramNameSuffix;
3532 fHistKaontotET = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3533 histname = "fHistKaonAcctotET_" + fHistogramNameSuffix;
3534 fHistKaonAcctotET = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3535 histname = "fHistKaonRectotET_" + fHistogramNameSuffix;
3536 fHistKaonRectotET = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3537 histname = "fHistKaonMatchtotET_" + fHistogramNameSuffix;
3538 fHistKaonMatchtotET = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3539
3540 histname = "fHistKaonRectotETDep_" + fHistogramNameSuffix;
3541 fHistKaonRectotETDep = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3542 histname = "fHistKaonMatchtotETDep_" + fHistogramNameSuffix;
3543 fHistKaonMatchtotETDep = new TH1F(histname.Data(),"total ET, MC primary Kaons",fgNumOfEBins, fgEAxis);
3544
3545 fHistKaonRecEtaEDepETDep = CreateEtaEHisto2D("fHistKaonRecEtaEDepETDep_","MC E_{T}, primary Kaons, reconstructed","E_{T} dep (GeV)");
3546 fHistKaonMatchEtaEDepETDep = CreateEtaEHisto2D("fHistKaonMatchEtaEDepETDep_","MC E_{T}, primary Kaons, tracking matched","E_{T} dep (GeV)");
3547
3548 fHistKaonRecEtaPtETDep = CreateEtaPtHisto2D("fHistKaonRecEtaPtETDep_","MC E_{T}, primary Kaons","E_{T} dep (GeV)");
3549 fHistKaonMatchEtaPtETDep = CreateEtaPtHisto2D("fHistKaonMatchEtaPtETDep_","MC E_{T}, primary Kaons","E_{T} dep(GeV)");
3550
3551 fHistKaonRecEtaETDep = CreateEtaEtHisto2D("fHistKaonRecEtaETDep_","MC primary Kaons","#");
3552 fHistKaonMatchEtaETDep = CreateEtaEtHisto2D("fHistKaonMatchEtaETDep_","MC primary Kaons","#");
3553
3554 fHistKaonRecResEET = CreateResEHisto2D("fHistKaonRecResEET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3555 fHistKaonRecResPtET = CreateResPtHisto2D("fHistKaonRecResPtET_","MC E_{T}, primary Kaons","E_{T}(GeV)");
3556 fHistKaonRecResE = CreateResEHisto2D("fHistKaonRecResE_","MC primary Kaons","#");
3557 fHistKaonRecResPt = CreateResPtHisto2D("fHistKaonRecResPt_","MC primary Kaons","#");
3558 fHistKaonRecResEDepETDep = CreateResEHisto2D("fHistKaonRecResEDepETDep_","MC E_{T}, primary Kaons","E_{T} dep (GeV)");
3559 fHistKaonRecResPtETDep = CreateResPtHisto2D("fHistKaonRecResPtETDep_","MC E_{T}, primary Kaons","E_{T} dep (GeV)");
3560
3561 fHistProtonEtaEET = CreateEtaEHisto2D("fHistProtonEtaEET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3562 fHistProtonAccEtaEET = CreateEtaEHisto2D("fHistProtonAccEtaEET_","MC E_{T}, primary Protons, inside EMCal acceptance","E_{T}(GeV)");
3563 fHistProtonRecEtaEET = CreateEtaEHisto2D("fHistProtonRecEtaEET_","MC E_{T}, primary Protons, reconstructed","E_{T}(GeV)");
3564 fHistProtonMatchEtaEET = CreateEtaEHisto2D("fHistProtonMatchEtaEET_","MC E_{T}, primary Protons, tracking matched","E_{T}(GeV)");
3565
3566 fHistProtonEtaPtET = CreateEtaPtHisto2D("fHistProtonEtaPtET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3567 fHistProtonAccEtaPtET = CreateEtaPtHisto2D("fHistProtonAccEtaPtET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3568 fHistProtonRecEtaPtET = CreateEtaPtHisto2D("fHistProtonRecEtaPtET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3569 fHistProtonMatchEtaPtET = CreateEtaPtHisto2D("fHistProtonMatchEtaPtET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3570
3571 fHistProtonEtaET = CreateEtaEtHisto2D("fHistProtonEtaET_","MC primary Protons","#");
3572 fHistProtonAccEtaET = CreateEtaEtHisto2D("fHistProtonAccEtaET_","MC primary Protons","#");
3573 fHistProtonRecEtaET = CreateEtaEtHisto2D("fHistProtonRecEtaET_","MC primary Protons","#");
3574 fHistProtonMatchEtaET = CreateEtaEtHisto2D("fHistProtonMatchEtaET_","MC primary Protons","#");
3575
3576 fHistProtonEtaE = CreateEtaEHisto2D("fHistProtonEtaE_","MC primary Protons","#");
3577 fHistProtonAccEtaE = CreateEtaEHisto2D("fHistProtonAccEtaE_","MC primary Protons","#");
3578 fHistProtonRecEtaE = CreateEtaEHisto2D("fHistProtonRecEtaE_","MC primary Protons","#");
3579 fHistProtonMatchEtaE = CreateEtaEHisto2D("fHistProtonMatchEtaE_","MC primary Protons","#");
3580
3581 fHistProtonEtaPt = CreateEtaPtHisto2D("fHistProtonEtaPt_","MC primary Protons","#");
3582 fHistProtonAccEtaPt = CreateEtaPtHisto2D("fHistProtonAccEtaPt_","MC primary Protons","#");
3583 fHistProtonRecEtaPt = CreateEtaPtHisto2D("fHistProtonRecEtaPt_","MC primary Protons","#");
3584 fHistProtonMatchEtaPt = CreateEtaPtHisto2D("fHistProtonMatchEtaPt_","MC primary Protons","#");
3585
3586 histname = "fHistProtontotET_" + fHistogramNameSuffix;
3587 fHistProtontotET = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3588 histname = "fHistProtonAcctotET_" + fHistogramNameSuffix;
3589 fHistProtonAcctotET = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3590 histname = "fHistProtonRectotET_" + fHistogramNameSuffix;
3591 fHistProtonRectotET = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3592 histname = "fHistProtonMatchtotET_" + fHistogramNameSuffix;
3593 fHistProtonMatchtotET = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3594
3595 histname = "fHistProtonRectotETDep_" + fHistogramNameSuffix;
3596 fHistProtonRectotETDep = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3597 histname = "fHistProtonMatchtotETDep_" + fHistogramNameSuffix;
3598 fHistProtonMatchtotETDep = new TH1F(histname.Data(),"total ET, MC primary Protons",fgNumOfEBins, fgEAxis);
3599
3600 fHistProtonRecEtaEDepETDep = CreateEtaEHisto2D("fHistProtonRecEtaEDepETDep_","MC E_{T}, primary Protons, reconstructed","E_{T} dep (GeV)");
3601 fHistProtonMatchEtaEDepETDep = CreateEtaEHisto2D("fHistProtonMatchEtaEDepETDep_","MC E_{T}, primary Protons, tracking matched","E_{T} dep (GeV)");
3602
3603 fHistProtonRecEtaPtETDep = CreateEtaPtHisto2D("fHistProtonRecEtaPtETDep_","MC E_{T}, primary Protons","E_{T} dep (GeV)");
3604 fHistProtonMatchEtaPtETDep = CreateEtaPtHisto2D("fHistProtonMatchEtaPtETDep_","MC E_{T}, primary Protons","E_{T} dep(GeV)");
3605
3606 fHistProtonRecEtaETDep = CreateEtaEtHisto2D("fHistProtonRecEtaETDep_","MC primary Protons","#");
3607 fHistProtonMatchEtaETDep = CreateEtaEtHisto2D("fHistProtonMatchEtaETDep_","MC primary Protons","#");
3608
3609 fHistProtonRecResEET = CreateResEHisto2D("fHistProtonRecResEET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3610 fHistProtonRecResPtET = CreateResPtHisto2D("fHistProtonRecResPtET_","MC E_{T}, primary Protons","E_{T}(GeV)");
3611 fHistProtonRecResE = CreateResEHisto2D("fHistProtonRecResE_","MC primary Protons","#");
3612 fHistProtonRecResPt = CreateResPtHisto2D("fHistProtonRecResPt_","MC primary Protons","#");
3613 fHistProtonRecResEDepETDep = CreateResEHisto2D("fHistProtonRecResEDepETDep_","MC E_{T}, primary Protons","E_{T} dep (GeV)");
3614 fHistProtonRecResPtETDep = CreateResPtHisto2D("fHistProtonRecResPtETDep_","MC E_{T}, primary Protons","E_{T} dep (GeV)");
3615
3616 histname = "fHistTotChargedtotET_" + fHistogramNameSuffix;
3617 fHistTotChargedtotET = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3618 histname = "fHistTotChargedAcctotET_" + fHistogramNameSuffix;
3619 fHistTotChargedAcctotET = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3620 histname = "fHistTotChargedRectotET_" + fHistogramNameSuffix;
3621 fHistTotChargedRectotET = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3622 histname = "fHistTotChargedMatchtotET_" + fHistogramNameSuffix;
3623 fHistTotChargedMatchtotET = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3624
3625 histname = "fHistTotChargedRectotETDep_" + fHistogramNameSuffix;
3626 fHistTotChargedRectotETDep = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3627 histname = "fHistTotChargedMatchtotETDep_" + fHistogramNameSuffix;
3628 fHistTotChargedMatchtotETDep = new TH1F(histname.Data(),"total ET, MC primary charged particles",fgNumOfEBins, fgEAxis);
3629
3630 fHistNeutronEtaEET = CreateEtaEHisto2D("fHistNeutronEtaEET_","MC E_{T}, primary Neutrons","E_{T}(GeV)");
3631 fHistNeutronAccEtaEET = CreateEtaEHisto2D("fHistNeutronAccEtaEET_","MC E_{T}, primary Neutrons, inside EMCal acceptance","E_{T}(GeV)");
3632 fHistNeutronRecEtaEET = CreateEtaEHisto2D("fHistNeutronRecEtaEET_","MC E_{T}, primary Neutrons, reconstructed","E_{T}(GeV)");
3633
3634 fHistNeutronEtaPtET = CreateEtaPtHisto2D("fHistNeutronEtaPtET_","MC E_{T}, primary Neutrons","E_{T}(GeV)");
3635 fHistNeutronAccEtaPtET = CreateEtaPtHisto2D("fHistNeutronAccEtaPtET_","MC E_{T}, primary Neutrons","E_{T}(GeV)");
3636 fHistNeutronRecEtaPtET = CreateEtaPtHisto2D("fHistNeutronRecEtaPtET_","MC E_{T}, primary Neutrons","E_{T}(GeV)");
3637
3638 fHistNeutronEtaET = CreateEtaEtHisto2D("fHistNeutronEtaET_","MC primary Neutrons","#");
3639 fHistNeutronAccEtaET = CreateEtaEtHisto2D("fHistNeutronAccEtaET_","MC primary Neutrons","#");
3640 fHistNeutronRecEtaET = CreateEtaEtHisto2D("fHistNeutronRecEtaET_","MC primary Neutrons","#");
3641
3642 fHistNeutronEtaE = CreateEtaEHisto2D("fHistNeutronEtaE_","MC primary Neutrons","#");
3643 fHistNeutronAccEtaE = CreateEtaEHisto2D("fHistNeutronAccEtaE_","MC primary Neutrons","#");
3644 fHistNeutronRecEtaE = CreateEtaEHisto2D("fHistNeutronRecEtaE_","MC primary Neutrons","#");
3645
3646 fHistNeutronEtaPt = CreateEtaPtHisto2D("fHistNeutronEtaPt_","MC primary Neutrons","#");
3647 fHistNeutronAccEtaPt = CreateEtaPtHisto2D("fHistNeutronAccEtaPt_","MC primary Neutrons","#");
3648 fHistNeutronRecEtaPt = CreateEtaPtHisto2D("fHistNeutronRecEtaPt_","MC primary Neutrons","#");
3649
3650 histname = "fHistNeutrontotET_" + fHistogramNameSuffix;
3651 fHistNeutrontotET = new TH1F(histname.Data(),"total ET, MC primary Neutrons",fgNumOfEBins, fgEAxis);
3652 histname = "fHistNeutronAcctotET_" + fHistogramNameSuffix;
3653 fHistNeutronAcctotET = new TH1F(histname.Data(),"total ET, MC primary Neutrons",fgNumOfEBins, fgEAxis);
3654 histname = "fHistNeutronRectotET_" + fHistogramNameSuffix;
3655 fHistNeutronRectotET = new TH1F(histname.Data(),"total ET, MC primary Neutrons",fgNumOfEBins, fgEAxis);
3656 histname = "fHistNeutronRectotETDep_" + fHistogramNameSuffix;
3657 fHistNeutronRectotETDep = new TH1F(histname.Data(),"total ET, MC primary Neutrons",fgNumOfEBins, fgEAxis);
3658
3659 fHistNeutronRecEtaEDepETDep = CreateEtaEHisto2D("fHistNeutronRecEtaEDepETDep_","MC E_{T}, primary Neutrons, reconstructed","E_{T} dep (GeV)");
3660 fHistNeutronRecEtaETDep = CreateEtaEtHisto2D("fHistNeutronRecEtaETDep_","MC primary Neutrons","#");
3661
3662 fHistNeutronRecEtaPtETDep = CreateEtaPtHisto2D("fHistNeutronRecEtaPtETDep_","MC E_{T}, primary Neutrons","E_{T} dep (GeV)");
3663
3664 fHistK0EtaEET = CreateEtaEHisto2D("fHistK0EtaEET_","MC E_{T}, K0S daughters","E_{T}(GeV)");
3665 fHistK0RecEtaEET = CreateEtaEHisto2D("fHistK0RecEtaEET_","MC E_{T}, K0S daughters, reconstructed","E_{T}(GeV)");
3666
3667 fHistK0EtaPtET = CreateEtaPtHisto2D("fHistK0EtaPtET_","MC E_{T}, K0S daughters","E_{T}(GeV)");
3668 fHistK0RecEtaPtET = CreateEtaPtHisto2D("fHistK0RecEtaPtET_","MC E_{T}, K0S daughters","E_{T}(GeV)");
3669
3670 fHistK0EtaET = CreateEtaEtHisto2D("fHistK0EtaET_","MC K0S daughters","#");
3671 fHistK0RecEtaET = CreateEtaEtHisto2D("fHistK0RecEtaET_","MC K0S daughters","#");
3672
3673 fHistK0EtaE = CreateEtaEHisto2D("fHistK0EtaE_","MC K0S daughters","#");
3674 fHistK0RecEtaE = CreateEtaEHisto2D("fHistK0RecEtaE_","MC K0S daughters","#");
3675
3676 fHistK0EtaPt = CreateEtaPtHisto2D("fHistK0EtaPt_","MC K0S daughters","#");
3677 fHistK0RecEtaPt = CreateEtaPtHisto2D("fHistK0RecEtaPt_","MC K0S daughters","#");
3678
3679 histname = "fHistK0totET_" + fHistogramNameSuffix;
3680 fHistK0totET = new TH1F(histname.Data(),"total ET, MC K0s daughters",fgNumOfEBins, fgEAxis);
3681 histname = "fHistK0RectotET_" + fHistogramNameSuffix;
3682 fHistK0RectotET = new TH1F(histname.Data(),"total ET, MC K0s daughters",fgNumOfEBins, fgEAxis);
3683 histname = "fHistK0RectotETDep_" + fHistogramNameSuffix;
3684 fHistK0RectotETDep = new TH1F(histname.Data(),"total ET, MC K0s daughters",fgNumOfEBins, fgEAxis);
3685
3686 fHistK0RecEtaEDepETDep = CreateEtaEHisto2D("fHistK0RecEtaEDepETDep_","MC E_{T}, MC K0s daughters, reconstructed","E_{T} dep (GeV)");
3687 fHistK0RecEtaETDep = CreateEtaEtHisto2D("fHistK0RecEtaETDep_","MC K0s daughters","#");
3688
3689 fHistK0RecEtaPtETDep = CreateEtaPtHisto2D("fHistK0RecEtaPtETDep_","MC E_{T}, MC K0s daughters","E_{T} dep (GeV)");
3690
3691 fHistLambdaEtaEET = CreateEtaEHisto2D("fHistLambdaEtaEET_","MC E_{T}, Lambda daughters","E_{T}(GeV)");
3692 fHistLambdaRecEtaEET = CreateEtaEHisto2D("fHistLambdaRecEtaEET_","MC E_{T}, Lambda daughters, reconstructed","E_{T}(GeV)");
3693
3694 fHistLambdaEtaPtET = CreateEtaPtHisto2D("fHistLambdaEtaPtET_","MC E_{T}, Lambda daughters","E_{T}(GeV)");
3695 fHistLambdaRecEtaPtET = CreateEtaPtHisto2D("fHistLambdaRecEtaPtET_","MC E_{T}, Lambda daughters","E_{T}(GeV)");
3696
3697 fHistLambdaEtaET = CreateEtaEtHisto2D("fHistLambdaEtaET_","MC Lambda daughters","#");
3698 fHistLambdaRecEtaET = CreateEtaEtHisto2D("fHistLambdaRecEtaET_","MC Lambda daughters","#");
3699
3700 fHistLambdaEtaE = CreateEtaEHisto2D("fHistLambdaEtaE_","MC Lambda daughters","#");
3701 fHistLambdaRecEtaE = CreateEtaEHisto2D("fHistLambdaRecEtaE_","MC Lambda daughters","#");
3702
3703 fHistLambdaEtaPt = CreateEtaPtHisto2D("fHistLambdaEtaPt_","MC Lambda daughters","#");
3704 fHistLambdaRecEtaPt = CreateEtaPtHisto2D("fHistLambdaRecEtaPt_","MC Lambda daughters","#");
3705
3706 histname = "fHistLambdatotET_" + fHistogramNameSuffix;
3707 fHistLambdatotET = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3708 histname = "fHistLambdaRectotET_" + fHistogramNameSuffix;
3709 fHistLambdaRectotET = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3710 histname = "fHistLambdaRectotETDep_" + fHistogramNameSuffix;
3711 fHistLambdaRectotETDep = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3712
3713 fHistLambdaRecEtaEDepETDep = CreateEtaEHisto2D("fHistLambdaRecEtaEDepETDep_","MC E_{T}, MC Lambdas daughters, reconstructed","E_{T} dep (GeV)");
3714 fHistLambdaRecEtaETDep = CreateEtaEtHisto2D("fHistLambdaRecEtaETDep_","MC Lambdas daughters","#");
3715
3716 fHistLambdaRecEtaPtETDep = CreateEtaPtHisto2D("fHistLambdaRecEtaPtETDep_","MC E_{T}, MC Lambdas daughters","E_{T} dep (GeV)");
3717
3718 histname = "fHistTotNeutraltotET_" + fHistogramNameSuffix;
3719 fHistTotNeutraltotET = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3720 histname = "fHistTotNeutralRectotET_" + fHistogramNameSuffix;
3721 fHistTotNeutralRectotET = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3722 histname = "fHistTotNeutralRectotETDep_" + fHistogramNameSuffix;
3723 fHistTotNeutralRectotETDep = new TH1F(histname.Data(),"total ET, MC Lambdas daughters",fgNumOfEBins, fgEAxis);
3724
3725 histname = "fHistTotaltotET_" + fHistogramNameSuffix;
3726 fHistTotaltotET = new TH1F(histname.Data(),"total ET, all particles",fgNumOfEBins, fgEAxis);
3727 histname = "fHistTotalAcctotET_" + fHistogramNameSuffix;
3728 fHistTotalAcctotET = new TH1F(histname.Data(),"total ET, all particles",fgNumOfEBins, fgEAxis);
3729 histname = "fHistTotalRectotET_" + fHistogramNameSuffix;
3730 fHistTotalRectotET = new TH1F(histname.Data(),"total ET, all particles",fgNumOfEBins, fgEAxis);
3731 histname = "fHistTotalRectotETDep_" + fHistogramNameSuffix;
3732 fHistTotalRectotETDep = new TH1F(histname.Data(),"total ET, all particles",fgNumOfEBins, fgEAxis);
3733
3734 histname = "fHistAll_ERecvsMC_" + fHistogramNameSuffix;
3735 fHistAllERecEMC = new TH2F(histname.Data(),"E cluster Rec vs MC, all particles",fgNumOfEBins, fgEAxis,fgNumOfEBins, fgEAxis);
3736 fHistAllERecEMC->SetXTitle("E_{MC}(GeV)");
3737 fHistAllERecEMC->SetYTitle("E_{Rec}(GeV)");
3738
3739 histname = "fHistElectron_ERecvsMC_" + fHistogramNameSuffix;
3740 fHistElectronERecEMC = new TH2F(histname.Data(),"E cluster Rec vs MC, Electrons",fgNumOfEBins, fgEAxis,fgNumOfEBins, fgEAxis);
3741 fHistElectronERecEMC->SetXTitle("E_{MC}(GeV)");
3742 fHistElectronERecEMC->SetYTitle("E_{Rec}(GeV)");
3743
3744 histname = "fHistGamma_ERecvsMC_" + fHistogramNameSuffix;
3745 fHistGammaERecEMC = new TH2F(histname.Data(),"E cluster Rec vs MC, Gammas",fgNumOfEBins, fgEAxis,fgNumOfEBins, fgEAxis);
3746 fHistGammaERecEMC->SetXTitle("E_{MC}(GeV)");
3747 fHistGammaERecEMC->SetYTitle("E_{Rec}(GeV)");
3748
3749 histname = "fHistAllPtRecPtMC_" + fHistogramNameSuffix;
3750 fHistAllPtRecPtMC = new TH2F(histname.Data(),"pt track Rec vs MC, all particles",fgNumOfEBins, fgEAxis,fgNumOfEBins, fgEAxis);
3751 fHistAllPtRecPtMC->SetXTitle("p_{T}^{MC}(GeV/c)");
3752 fHistAllPtRecPtMC->SetYTitle("p_{T}^{Rec}(GeV/c)");
3753
3754 histname = "fHistChargedRes_" + fHistogramNameSuffix;
3755 fHistChargedRes = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position), charged particles",200,-0.1,0.1,200,-0.1,0.1);
3756 fHistChargedRes->SetXTitle("#Delta#phi");
3757 fHistChargedRes->SetYTitle("#Delta#eta");
3758
3759 histname = "fHistChargedRes2_" + fHistogramNameSuffix;
3760 fHistChargedRes2 = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position), charged particles",200,-0.1,0.1,200,-0.1,0.1);
3761 fHistChargedRes2->SetXTitle("#Delta#phi");
3762 fHistChargedRes2->SetYTitle("#Delta#eta");
3763
3764 histname = "fHistChargedRes3_" + fHistogramNameSuffix;
3765 fHistChargedRes3 = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position), charged particles",200,-0.1,0.1,200,-0.1,0.1);
3766 fHistChargedRes3->SetXTitle("#Delta#phi");
3767 fHistChargedRes3->SetYTitle("#Delta#eta");
3768
3769 histname = "fHistNeutralRes_" + fHistogramNameSuffix;
3770 fHistNeutralRes = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position), neutral particles",200,-0.1,0.1,200,-0.1,0.1);
3771 fHistNeutralRes->SetXTitle("#Delta#phi");
3772 fHistNeutralRes->SetYTitle("#Delta#eta");
3773
3774 histname = "fHistElectronRes_" + fHistogramNameSuffix;
3775 fHistElectronRes = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position, Electrons",200,-0.1,0.1,200,-0.1,0.1);
3776 fHistElectronRes->SetXTitle("#Delta#phi");
3777 fHistElectronRes->SetYTitle("#Delta#eta");
3778
3779 histname = "fHistGammaRes_" + fHistogramNameSuffix;
3780 fHistGammaRes = new TH2F(histname.Data(),"#Delta#phi vs #Delta#eta (track projection - cluster position, Gammas",200,-0.1,0.1,200,-0.1,0.1);
3781 fHistGammaRes->SetXTitle("#Delta#phi");
3782 fHistGammaRes->SetYTitle("#Delta#eta");
3783
3784 histname = "fHistIsInAcc_" + fHistogramNameSuffix;
3785 //fHistIsInAcc = new TH2F(histname.Data(),"X,Y position of particle projection inside EMCal",1201,-600.5,600.5,1201,-600.5,600.5);
3786 //fHistIsInAcc->SetXTitle("X (cm)");
3787 //fHistIsInAcc->SetYTitle("Y (cm)");
3788 fHistIsInAcc = new TH2F(histname.Data(),"#phhi, #eta position of particle projection inside EMCal",360,0.,360.,200,-1.,1.);
3789 fHistIsInAcc->SetXTitle("#phi");
3790 fHistIsInAcc->SetYTitle("#eta");
3791
3792 histname = "fHistElectronFirstMother_" + fHistogramNameSuffix;
3793 fHistElectronFirstMother = new TH1F(histname.Data(),"Electron First Mother PDG Code",1201,-600.5,600.5);
3794 histname = "fHistElectronFirstMotherXY_" + fHistogramNameSuffix;
3795 fHistElectronFirstMotherXY = new TH2F(histname.Data(),"Electron Mother X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3796 histname = "fHistElectronNDaughters_" + fHistogramNameSuffix;
3797 fHistElectronNDaughters = new TH1F(histname.Data(),"Number of Electron Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3798 histname = "fHistElectronDaughters_" + fHistogramNameSuffix;
3799 fHistElectronDaughters = new TH1F(histname.Data(),"Electron Daugther PDG Code",1201,-600.5,600.5);
3800 histname = "fHistElectronDaughtersXY_" + fHistogramNameSuffix;
3801 fHistElectronDaughtersXY = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3802
3803 histname = "fHistElectronFirstMotherAcc_" + fHistogramNameSuffix;
3804 fHistElectronFirstMotherAcc = new TH1F(histname.Data(),"Electron First Mother PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3805 histname = "fHistElectronFirstMotherXYAcc_" + fHistogramNameSuffix;
3806 fHistElectronFirstMotherXYAcc = new TH2F(histname.Data(),"Electron Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3807 histname = "fHistElectronNDaughtersAcc_" + fHistogramNameSuffix;
3808 fHistElectronNDaughtersAcc = new TH1F(histname.Data(),"Number of Electron Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3809 histname = "fHistElectronDaughtersAcc_" + fHistogramNameSuffix;
3810 fHistElectronDaughtersAcc = new TH1F(histname.Data(),"Electron Daugther PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3811 histname = "fHistElectronDaughtersXYAcc_" + fHistogramNameSuffix;
3812 fHistElectronDaughtersXYAcc = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3813
3814 histname = "fHistElectronFirstMotherRec_" + fHistogramNameSuffix;
3815 fHistElectronFirstMotherRec = new TH1F(histname.Data(),"Reconstructed Electron First Mother PDG Code",1201,-600.5,600.5);
3816 histname = "fHistElectronFirstMotherXYRec_" + fHistogramNameSuffix;
3817 fHistElectronFirstMotherXYRec = new TH2F(histname.Data(),"Electron Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3818 histname = "fHistElectronNDaughtersRec_" + fHistogramNameSuffix;
3819 fHistElectronNDaughtersRec = new TH1F(histname.Data(),"Number of Electron Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3820 histname = "fHistElectronDaughtersRec_" + fHistogramNameSuffix;
3821 fHistElectronDaughtersRec = new TH1F(histname.Data(),"Electron Daugther PDG Code",1201,-600.5,600.5);
3822 histname = "fHistElectronDaughtersXYRec_" + fHistogramNameSuffix;
3823 fHistElectronDaughtersXYRec = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3824
3825 histname = "fHistNPPElectronFirstMother_" + fHistogramNameSuffix;
3826 fHistNPPElectronFirstMother = new TH1F(histname.Data(),"Electron First Mother PDG Code",1201,-600.5,600.5);
3827 histname = "fHistNPPElectronFirstMotherXY_" + fHistogramNameSuffix;
3828 fHistNPPElectronFirstMotherXY = new TH2F(histname.Data(),"Electron Mother X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3829 histname = "fHistNPPElectronNDaughters_" + fHistogramNameSuffix;
3830 fHistNPPElectronNDaughters = new TH1F(histname.Data(),"Number of Electron Daugthers",11,-0.5,10.5);
3831 histname = "fHistNPPElectronDaughters_" + fHistogramNameSuffix;
3832 fHistNPPElectronDaughters = new TH1F(histname.Data(),"Electron Daugther PDG Code",1201,-600.5,600.5);
3833 histname = "fHistNPPElectronDaughtersXY_" + fHistogramNameSuffix;
3834 fHistNPPElectronDaughtersXY = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3835
3836 histname = "fHistNPPElectronFirstMotherAcc_" + fHistogramNameSuffix;
3837 fHistNPPElectronFirstMotherAcc = new TH1F(histname.Data(),"Electron First Mother PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3838 histname = "fHistNPPElectronFirstMotherXYAcc_" + fHistogramNameSuffix;
3839 fHistNPPElectronFirstMotherXYAcc = new TH2F(histname.Data(),"Electron Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3840 histname = "fHistNPPElectronNDaughtersAcc_" + fHistogramNameSuffix;
3841 fHistNPPElectronNDaughtersAcc = new TH1F(histname.Data(),"Number of Electron Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3842 histname = "fHistNPPElectronDaughtersAcc_" + fHistogramNameSuffix;
3843 fHistNPPElectronDaughtersAcc = new TH1F(histname.Data(),"Electron Daugther PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3844 histname = "fHistNPPElectronDaughtersXYAcc_" + fHistogramNameSuffix;
3845 fHistNPPElectronDaughtersXYAcc = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3846
3847 histname = "fHistNPPElectronFirstMotherRec_" + fHistogramNameSuffix;
3848 fHistNPPElectronFirstMotherRec = new TH1F(histname.Data(),"Reconstructed Electron First Mother PDG Code",1201,-600.5,600.5);
3849 histname = "fHistNPPElectronFirstMotherXYRec_" + fHistogramNameSuffix;
3850 fHistNPPElectronFirstMotherXYRec = new TH2F(histname.Data(),"Electron Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3851 histname = "fHistNPPElectronNDaughtersRec_" + fHistogramNameSuffix;
3852 fHistNPPElectronNDaughtersRec = new TH1F(histname.Data(),"Number of Electron Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3853 histname = "fHistNPPElectronDaughtersRec_" + fHistogramNameSuffix;
3854 fHistNPPElectronDaughtersRec = new TH1F(histname.Data(),"Electron Daugther PDG Code",1201,-600.5,600.5);
3855 histname = "fHistNPPElectronDaughtersXYRec_" + fHistogramNameSuffix;
3856 fHistNPPElectronDaughtersXYRec = new TH2F(histname.Data(),"Electron Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3857
3858 histname = "fHistGammaFirstMother_" + fHistogramNameSuffix;
3859 fHistGammaFirstMother = new TH1F(histname.Data(),"Gamma First Mother PDG Code",1201,-600.5,600.5);
3860 histname = "fHistGammaFirstMotherXY_" + fHistogramNameSuffix;
3861 fHistGammaFirstMotherXY = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3862 histname = "fHistGammaNDaughters_" + fHistogramNameSuffix;
3863 fHistGammaNDaughters = new TH1F(histname.Data(),"Number of Gamma Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3864 histname = "fHistGammaDaughters_" + fHistogramNameSuffix;
3865 fHistGammaDaughters = new TH1F(histname.Data(),"Gamma Daugther PDG Code",1201,-600.5,600.5);
3866 histname = "fHistGammaDaughtersXY_" + fHistogramNameSuffix;
3867 fHistGammaDaughtersXY = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3868 histname = "fHistConvGammaDaughtersXY_" + fHistogramNameSuffix;
3869 fHistConvGammaDaughtersXY = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3870 histname = "fHistNonConvGammaDaughtersXY_" + fHistogramNameSuffix;
3871 fHistNonConvGammaDaughtersXY = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3872
3873 histname = "fHistGammaFirstMotherAcc_" + fHistogramNameSuffix;
3874 fHistGammaFirstMotherAcc = new TH1F(histname.Data(),"Gamma First Mother PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3875 histname = "fHistGammaFirstMotherXYAcc_" + fHistogramNameSuffix;
3876 fHistGammaFirstMotherXYAcc = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3877 histname = "fHistGammaNDaughtersAcc_" + fHistogramNameSuffix;
3878 fHistGammaNDaughtersAcc = new TH1F(histname.Data(),"Number of Gamma Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3879 histname = "fHistGammaDaughtersAcc_" + fHistogramNameSuffix;
3880 fHistGammaDaughtersAcc = new TH1F(histname.Data(),"Gamma Daugther PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3881 histname = "fHistGammaDaughtersXYAcc_" + fHistogramNameSuffix;
3882 fHistGammaDaughtersXYAcc = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3883 histname = "fHistConvGammaDaughtersXYAcc_" + fHistogramNameSuffix;
3884 fHistConvGammaDaughtersXYAcc = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3885 histname = "fHistNonConvGammaDaughtersXYAcc_" + fHistogramNameSuffix;
3886 fHistNonConvGammaDaughtersXYAcc = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3887
3888 histname = "fHistGammaFirstMotherRec_" + fHistogramNameSuffix;
3889 fHistGammaFirstMotherRec = new TH1F(histname.Data(),"Reconstructed Gamma First Mother PDG Code",1201,-600.5,600.5);
3890 histname = "fHistGammaFirstMotherXYRec_" + fHistogramNameSuffix;
3891 fHistGammaFirstMotherXYRec = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3892 histname = "fHistGammaNDaughtersRec_" + fHistogramNameSuffix;
3893 fHistGammaNDaughtersRec = new TH1F(histname.Data(),"Number of Gamma Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3894 histname = "fHistGammaDaughtersRec_" + fHistogramNameSuffix;
3895 fHistGammaDaughtersRec = new TH1F(histname.Data(),"Gamma Daugther PDG Code",1201,-600.5,600.5);
3896 histname = "fHistGammaDaughtersXYRec_" + fHistogramNameSuffix;
3897 fHistGammaDaughtersXYRec = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3898 histname = "fHistConvGammaDaughtersXYRec_" + fHistogramNameSuffix;
3899 fHistConvGammaDaughtersXYRec = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3900 histname = "fHistNonConvGammaDaughtersXYRec_" + fHistogramNameSuffix;
3901 fHistNonConvGammaDaughtersXYRec = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3902
3903 histname = "fHistNPPGammaFirstMother_" + fHistogramNameSuffix;
3904 fHistNPPGammaFirstMother = new TH1F(histname.Data(),"Gamma First Mother PDG Code",1201,-600.5,600.5);
3905 histname = "fHistNPPGammaFirstMotherXY_" + fHistogramNameSuffix;
3906 fHistNPPGammaFirstMotherXY = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3907 histname = "fHistNPPGammaNDaughters_" + fHistogramNameSuffix;
3908 fHistNPPGammaNDaughters = new TH1F(histname.Data(),"Number of Gamma Daugthers",11,-0.5,10.5);
3909 histname = "fHistNPPGammaDaughters_" + fHistogramNameSuffix;
3910 fHistNPPGammaDaughters = new TH1F(histname.Data(),"Gamma Daugther PDG Code",1201,-600.5,600.5);
3911 histname = "fHistNPPGammaDaughtersXY_" + fHistogramNameSuffix;
3912 fHistNPPGammaDaughtersXY = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
3913
3914 histname = "fHistNPPGammaFirstMotherAcc_" + fHistogramNameSuffix;
3915 fHistNPPGammaFirstMotherAcc = new TH1F(histname.Data(),"Gamma First Mother PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3916 histname = "fHistNPPGammaFirstMotherXYAcc_" + fHistogramNameSuffix;
3917 fHistNPPGammaFirstMotherXYAcc = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3918 histname = "fHistNPPGammaNDaughtersAcc_" + fHistogramNameSuffix;
3919 fHistNPPGammaNDaughtersAcc = new TH1F(histname.Data(),"Number of Gamma Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3920 histname = "fHistNPPGammaDaughtersAcc_" + fHistogramNameSuffix;
3921 fHistNPPGammaDaughtersAcc = new TH1F(histname.Data(),"Gamma Daugther PDG Code, inside EMCal acceptance",1201,-600.5,600.5);
3922 histname = "fHistNPPGammaDaughtersXYAcc_" + fHistogramNameSuffix;
3923 fHistNPPGammaDaughtersXYAcc = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3924
3925 histname = "fHistNPPGammaFirstMotherRec_" + fHistogramNameSuffix;
3926 fHistNPPGammaFirstMotherRec = new TH1F(histname.Data(),"Reconstructed Gamma First Mother PDG Code",1201,-600.5,600.5);
3927 histname = "fHistNPPGammaFirstMotherXYRec_" + fHistogramNameSuffix;
3928 fHistNPPGammaFirstMotherXYRec = new TH2F(histname.Data(),"Gamma Mother X,Y vertex position, inside EMCal acceptance",1201,-600.5,600.5,1201,-600.5,600.5);
3929 histname = "fHistNPPGammaNDaughtersRec_" + fHistogramNameSuffix;
3930 fHistNPPGammaNDaughtersRec = new TH1F(histname.Data(),"Number of Gamma Daugthers, inside EMCal acceptance",11,-0.5,10.5);
3931 histname = "fHistNPPGammaDaughtersRec_" + fHistogramNameSuffix;
3932 fHistNPPGammaDaughtersRec = new TH1F(histname.Data(),"Gamma Daugther PDG Code",1201,-600.5,600.5);
3933 histname = "fHistNPPGammaDaughtersXYRec_" + fHistogramNameSuffix;
3934 fHistNPPGammaDaughtersXYRec = new TH2F(histname.Data(),"Gamma Daugther X,Y vertex position",1201,-600.5,600.5,1201,-600.5,600.5);
2e12941b 3935}
3936
3937void AliAnalysisEmEtMonteCarlo::FillOutputList(TList *list)
3938{//fill the output tlist
112b1d0c 3939 //AliAnalysisEt::FillOutputList(list);
3940
3941 list->Add(fHistPrimEtaEET);
3942 list->Add(fHistPrimEtaPtET);
3943 list->Add(fHistPrimEtaET);
3944 list->Add(fHistPrimtotET);
3945
3946 list->Add(fHistPrimAccEtaEET);
3947 list->Add(fHistPrimAccEtaPtET);
3948 list->Add(fHistPrimAccEtaET);
3949 list->Add(fHistPrimAcctotET);
3950
3951 list->Add(fHistPrimRecEtaEET);
3952 list->Add(fHistPrimRecEtaPtET);
3953 list->Add(fHistPrimRecEtaET);
3954 list->Add(fHistPrimRectotET);
3955
3956 list->Add(fHistPrimRecEtaEDepETDep);
3957 list->Add(fHistPrimRecEtaPtETDep);
3958 list->Add(fHistPrimRecEtaETDep);
3959 list->Add(fHistPrimRectotETDep);
3960
3961 list->Add(fHistElectronEtaEET);
3962 list->Add(fHistElectronEtaPtET);
3963 list->Add(fHistElectronEtaET);
3964 list->Add(fHistElectronEtaE);
3965 list->Add(fHistElectronEtaPt);
3966 list->Add(fHistElectrontotET);
3967
3968 list->Add(fHistConvElectronEtaEET);
3969 list->Add(fHistConvElectronEtaPtET);
3970 list->Add(fHistConvElectronEtaET);
3971 list->Add(fHistConvElectronEtaE);
3972 list->Add(fHistConvElectronEtaPt);
3973 list->Add(fHistConvElectrontotET);
3974
3975 list->Add(fHistScatElectronEtaEET);
3976 list->Add(fHistScatElectronEtaPtET);
3977 list->Add(fHistScatElectronEtaET);
3978 list->Add(fHistScatElectronEtaE);
3979 list->Add(fHistScatElectronEtaPt);
3980 list->Add(fHistScatElectrontotET);
3981
3982 list->Add(fHistTotElectrontotET);
3983
3984 list->Add(fHistGammaEtaEET);
3985 list->Add(fHistGammaEtaPtET);
3986 list->Add(fHistGammaEtaET);
3987 list->Add(fHistGammaEtaE);
3988 list->Add(fHistGammaEtaPt);
3989 list->Add(fHistGammatotET);
3990
3991 list->Add(fHistAnnihGammaEtaEET);
3992 list->Add(fHistAnnihGammaEtaPtET);
3993 list->Add(fHistAnnihGammaEtaET);
3994 list->Add(fHistAnnihGammaEtaE);
3995 list->Add(fHistAnnihGammaEtaPt);
3996 list->Add(fHistAnnihGammatotET);
3997
3998 list->Add(fHistScatGammaEtaEET);
3999 list->Add(fHistScatGammaEtaPtET);
4000 list->Add(fHistScatGammaEtaET);
4001 list->Add(fHistScatGammaEtaE);
4002 list->Add(fHistScatGammaEtaPt);
4003 list->Add(fHistScatGammatotET);
4004
4005 list->Add(fHistConvGammaEtaEET);
4006 list->Add(fHistConvGammaEtaPtET);
4007 list->Add(fHistConvGammaEtaET);
4008 list->Add(fHistConvGammaEtaE);
4009 list->Add(fHistConvGammaEtaPt);
4010 list->Add(fHistConvGammatotET);
4011
4012 list->Add(fHistNonConvGammaEtaEET);
4013 list->Add(fHistNonConvGammaEtaPtET);
4014 list->Add(fHistNonConvGammaEtaET);
4015 list->Add(fHistNonConvGammaEtaE);
4016 list->Add(fHistNonConvGammaEtaPt);
4017 list->Add(fHistNonConvGammatotET);
4018
4019 list->Add(fHistTotGammatotET);
4020
4021 list->Add(fHistTotEMtotET);
4022
4023 list->Add(fHistNPPElectronEtaEET);
4024 list->Add(fHistNPPElectronEtaPtET);
4025 list->Add(fHistNPPElectronEtaET);
4026 list->Add(fHistNPPElectronEtaE);
4027 list->Add(fHistNPPElectronEtaPt);
4028 list->Add(fHistNPPElectrontotET);
4029
4030 list->Add(fHistNPPGammaEtaEET);
4031 list->Add(fHistNPPGammaEtaPtET);
4032 list->Add(fHistNPPGammaEtaET);
4033 list->Add(fHistNPPGammaEtaE);
4034 list->Add(fHistNPPGammaEtaPt);
4035 list->Add(fHistNPPGammatotET);
4036
4037 list->Add(fHistTotNPPEMtotET);
4038
4039 list->Add(fHistNPPPi0GammaEtaEET);
4040 list->Add(fHistNPPPi0GammaEtaPtET);
4041 list->Add(fHistNPPPi0GammaEtaET);
4042 list->Add(fHistNPPPi0GammaEtaE);
4043 list->Add(fHistNPPPi0GammaEtaPt);
4044 list->Add(fHistNPPPi0GammatotET);
4045
4046 list->Add(fHistElectronAccEtaEET);
4047 list->Add(fHistElectronAccEtaPtET);
4048 list->Add(fHistElectronAccEtaET);
4049 list->Add(fHistElectronAccEtaE);
4050 list->Add(fHistElectronAccEtaPt);
4051 list->Add(fHistElectronAcctotET);
4052
4053 list->Add(fHistConvElectronAccEtaEET);
4054 list->Add(fHistConvElectronAccEtaPtET);
4055 list->Add(fHistConvElectronAccEtaET);
4056 list->Add(fHistConvElectronAccEtaE);
4057 list->Add(fHistConvElectronAccEtaPt);
4058 list->Add(fHistConvElectronAcctotET);
4059
4060 list->Add(fHistScatElectronAccEtaEET);
4061 list->Add(fHistScatElectronAccEtaPtET);
4062 list->Add(fHistScatElectronAccEtaET);
4063 list->Add(fHistScatElectronAccEtaE);
4064 list->Add(fHistScatElectronAccEtaPt);
4065 list->Add(fHistScatElectronAcctotET);
4066
4067 list->Add(fHistTotElectronAcctotET);
4068
4069 list->Add(fHistGammaAccEtaEET);
4070 list->Add(fHistGammaAccEtaPtET);
4071 list->Add(fHistGammaAccEtaET);
4072 list->Add(fHistGammaAccEtaE);
4073 list->Add(fHistGammaAccEtaPt);
4074 list->Add(fHistGammaAcctotET);
4075
4076 list->Add(fHistConvGammaAccEtaEET);
4077 list->Add(fHistConvGammaAccEtaPtET);
4078 list->Add(fHistConvGammaAccEtaET);
4079 list->Add(fHistConvGammaAccEtaE);
4080 list->Add(fHistConvGammaAccEtaPt);
4081 list->Add(fHistConvGammaAcctotET);
4082
4083 list->Add(fHistNonConvGammaAccEtaEET);
4084 list->Add(fHistNonConvGammaAccEtaPtET);
4085 list->Add(fHistNonConvGammaAccEtaET);
4086 list->Add(fHistNonConvGammaAccEtaE);
4087 list->Add(fHistNonConvGammaAccEtaPt);
4088 list->Add(fHistNonConvGammaAcctotET);
4089
4090 list->Add(fHistAnnihGammaAccEtaEET);
4091 list->Add(fHistAnnihGammaAccEtaPtET);
4092 list->Add(fHistAnnihGammaAccEtaET);
4093 list->Add(fHistAnnihGammaAccEtaE);
4094 list->Add(fHistAnnihGammaAccEtaPt);
4095 list->Add(fHistAnnihGammaAcctotET);
4096
4097 list->Add(fHistScatGammaAccEtaEET);
4098 list->Add(fHistScatGammaAccEtaPtET);
4099 list->Add(fHistScatGammaAccEtaET);
4100 list->Add(fHistScatGammaAccEtaE);
4101 list->Add(fHistScatGammaAccEtaPt);
4102 list->Add(fHistScatGammaAcctotET);
4103
4104 list->Add(fHistTotGammaAcctotET);
4105
4106 list->Add(fHistTotEMAcctotET);
4107
4108 list->Add(fHistNPPElectronAccEtaEET);
4109 list->Add(fHistNPPElectronAccEtaPtET);
4110 list->Add(fHistNPPElectronAccEtaE);
4111 list->Add(fHistNPPElectronAccEtaPt);
4112
4113 list->Add(fHistNPPGammaAccEtaEET);
4114 list->Add(fHistNPPGammaAccEtaPtET);
4115 list->Add(fHistNPPGammaAccEtaE);
4116 list->Add(fHistNPPGammaAccEtaPt);
4117
4118 list->Add(fHistElectronRecEtaEET);
4119 list->Add(fHistElectronRecEtaPtET);
4120 list->Add(fHistElectronRecEtaET);
4121 list->Add(fHistElectronRecEtaE);
4122 list->Add(fHistElectronRecEtaPt);
4123 list->Add(fHistElectronRectotET);
4124
4125 list->Add(fHistConvElectronRecEtaEET);
4126 list->Add(fHistConvElectronRecEtaPtET);
4127 list->Add(fHistConvElectronRecEtaET);
4128 list->Add(fHistConvElectronRecEtaE);
4129 list->Add(fHistConvElectronRecEtaPt);
4130 list->Add(fHistConvElectronRectotET);
4131
4132 list->Add(fHistScatElectronRecEtaEET);
4133 list->Add(fHistScatElectronRecEtaPtET);
4134 list->Add(fHistScatElectronRecEtaET);
4135 list->Add(fHistScatElectronRecEtaE);
4136 list->Add(fHistScatElectronRecEtaPt);
4137 list->Add(fHistScatElectronRectotET);
4138
4139 list->Add(fHistTotElectronRectotET);
4140
4141 list->Add(fHistGammaRecEtaEET);
4142 list->Add(fHistGammaRecEtaPtET);
4143 list->Add(fHistGammaRecEtaET);
4144 list->Add(fHistGammaRecEtaE);
4145 list->Add(fHistGammaRecEtaPt);
4146 list->Add(fHistGammaRectotET);
4147
4148
4149 list->Add(fHistGammaRecResEET);//Add Caio
4150 list->Add(fHistGammaRecResPtET);//Add CAio
4151 list->Add(fHistGammaRecResE);//Add Caio
4152 list->Add(fHistGammaRecResPt);//Add Caio
4153 list->Add(fHistGammaRecResEDepETDep);//Add Caio
4154 list->Add(fHistGammaRecResPtETDep); // Add Caio
4155
4156 list->Add(fHistAnnihGammaRecEtaEET);
4157 list->Add(fHistAnnihGammaRecEtaPtET);
4158 list->Add(fHistAnnihGammaRecEtaET);
4159 list->Add(fHistAnnihGammaRecEtaE);
4160 list->Add(fHistAnnihGammaRecEtaPt);
4161 list->Add(fHistAnnihGammaRectotET);
4162
4163 list->Add(fHistScatGammaRecEtaEET);
4164 list->Add(fHistScatGammaRecEtaPtET);
4165 list->Add(fHistScatGammaRecEtaET);
4166 list->Add(fHistScatGammaRecEtaE);
4167 list->Add(fHistScatGammaRecEtaPt);
4168 list->Add(fHistScatGammaRectotET);
4169
4170 list->Add(fHistTotGammaRectotET);
4171
4172 list->Add(fHistTotEMRectotET);
4173
4174 list->Add(fHistNPPElectronRecEtaEET);
4175 list->Add(fHistNPPElectronRecEtaPtET);
4176 list->Add(fHistNPPElectronRecEtaET);
4177 list->Add(fHistNPPElectronRecEtaE);
4178 list->Add(fHistNPPElectronRecEtaPt);
4179 list->Add(fHistNPPElectronRectotET);
4180
4181 list->Add(fHistNPPGammaRecEtaEET);
4182 list->Add(fHistNPPGammaRecEtaPtET);
4183 list->Add(fHistNPPGammaRecEtaET);
4184 list->Add(fHistNPPGammaRecEtaE);
4185 list->Add(fHistNPPGammaRecEtaPt);
4186 list->Add(fHistNPPGammaRectotET);
4187
4188 list->Add(fHistTotNPPEMRectotET);
4189
4190 list->Add(fHistNPPPi0GammaRecEtaEET);
4191 list->Add(fHistNPPPi0GammaRecEtaPtET);
4192 list->Add(fHistNPPPi0GammaRecEtaET);
4193 list->Add(fHistNPPPi0GammaRecEtaE);
4194 list->Add(fHistNPPPi0GammaRecEtaPt);
4195 list->Add(fHistNPPPi0GammaRectotET);
4196
4197 list->Add(fHistMuonEtaEET);
4198 list->Add(fHistMuonAccEtaEET);
4199 list->Add(fHistMuonRecEtaEET);
4200 list->Add(fHistMuonMatchEtaEET);
4201
4202 list->Add(fHistMuonEtaPtET);
4203 list->Add(fHistMuonAccEtaPtET);
4204 list->Add(fHistMuonRecEtaPtET);
4205 list->Add(fHistMuonMatchEtaPtET);
4206
4207 list->Add(fHistMuonEtaET);
4208 list->Add(fHistMuonAccEtaET);
4209 list->Add(fHistMuonRecEtaET);
4210 list->Add(fHistMuonMatchEtaET);
4211
4212 list->Add(fHistMuonEtaE);
4213 list->Add(fHistMuonAccEtaE);
4214 list->Add(fHistMuonRecEtaE);
4215 list->Add(fHistMuonMatchEtaE);
4216
4217 list->Add(fHistMuonEtaPt);
4218 list->Add(fHistMuonAccEtaPt);
4219 list->Add(fHistMuonRecEtaPt);
4220 list->Add(fHistMuonMatchEtaPt);
4221
4222 list->Add(fHistMuontotET);
4223 list->Add(fHistMuonAcctotET);
4224 list->Add(fHistMuonRectotET);
4225 list->Add(fHistMuonMatchtotET);
4226
4227 list->Add(fHistMuonRectotETDep);
4228 list->Add(fHistMuonMatchtotETDep);
4229
4230 list->Add(fHistMuonRecEtaEDepETDep);
4231 list->Add(fHistMuonMatchEtaEDepETDep);
4232
4233 list->Add(fHistMuonRecEtaPtETDep);
4234 list->Add(fHistMuonMatchEtaPtETDep);
4235
4236 list->Add(fHistMuonRecEtaETDep);
4237 list->Add(fHistMuonMatchEtaETDep);
4238
4239 list->Add(fHistMuonRecResEET);
4240 list->Add(fHistMuonRecResPtET);
4241 list->Add(fHistMuonRecResE);
4242 list->Add(fHistMuonRecResPt);
4243 list->Add(fHistMuonRecResEDepETDep);
4244 list->Add(fHistMuonRecResPtETDep);
4245
4246 list->Add(fHistPionEtaEET);
4247 list->Add(fHistPionAccEtaEET);
4248 list->Add(fHistPionRecEtaEET);
4249 list->Add(fHistPionMatchEtaEET);
4250
4251 list->Add(fHistPionEtaPtET);
4252 list->Add(fHistPionAccEtaPtET);
4253 list->Add(fHistPionRecEtaPtET);
4254 list->Add(fHistPionMatchEtaPtET);
4255
4256 list->Add(fHistPionEtaET);
4257 list->Add(fHistPionAccEtaET);
4258 list->Add(fHistPionRecEtaET);
4259 list->Add(fHistPionMatchEtaET);
4260
4261 list->Add(fHistPionEtaE);
4262 list->Add(fHistPionAccEtaE);
4263 list->Add(fHistPionRecEtaE);
4264 list->Add(fHistPionMatchEtaE);
4265
4266 list->Add(fHistPionEtaPt);
4267 list->Add(fHistPionAccEtaPt);
4268 list->Add(fHistPionRecEtaPt);
4269 list->Add(fHistPionMatchEtaPt);
4270
4271 list->Add(fHistPiontotET);
4272 list->Add(fHistPionAcctotET);
4273 list->Add(fHistPionRectotET);
4274 list->Add(fHistPionMatchtotET);
4275
4276 list->Add(fHistPionRectotETDep);
4277 list->Add(fHistPionMatchtotETDep);
4278
4279 list->Add(fHistPionRecEtaEDepETDep);
4280 list->Add(fHistPionMatchEtaEDepETDep);
4281
4282 list->Add(fHistPionRecEtaPtETDep);
4283 list->Add(fHistPionMatchEtaPtETDep);
4284
4285 list->Add(fHistPionRecEtaETDep);
4286 list->Add(fHistPionMatchEtaETDep);
4287
4288 list->Add(fHistPionRecResEET);
4289 list->Add(fHistPionRecResPtET);
4290 list->Add(fHistPionRecResE);
4291 list->Add(fHistPionRecResPt);
4292 list->Add(fHistPionRecResEDepETDep);
4293 list->Add(fHistPionRecResPtETDep);
4294
4295 list->Add(fHistKaonEtaEET);
4296 list->Add(fHistKaonAccEtaEET);
4297 list->Add(fHistKaonRecEtaEET);
4298 list->Add(fHistKaonMatchEtaEET);
4299
4300 list->Add(fHistKaonEtaPtET);
4301 list->Add(fHistKaonAccEtaPtET);
4302 list->Add(fHistKaonRecEtaPtET);
4303 list->Add(fHistKaonMatchEtaPtET);
4304
4305 list->Add(fHistKaonEtaET);
4306 list->Add(fHistKaonAccEtaET);
4307 list->Add(fHistKaonRecEtaET);
4308 list->Add(fHistKaonMatchEtaET);
4309
4310 list->Add(fHistKaonEtaE);
4311 list->Add(fHistKaonAccEtaE);
4312 list->Add(fHistKaonRecEtaE);
4313 list->Add(fHistKaonMatchEtaE);
4314
4315 list->Add(fHistKaonEtaPt);
4316 list->Add(fHistKaonAccEtaPt);
4317 list->Add(fHistKaonRecEtaPt);
4318 list->Add(fHistKaonMatchEtaPt);
4319
4320 list->Add(fHistKaontotET);
4321 list->Add(fHistKaonAcctotET);
4322 list->Add(fHistKaonRectotET);
4323 list->Add(fHistKaonMatchtotET);
4324
4325 list->Add(fHistKaonRectotETDep);
4326 list->Add(fHistKaonMatchtotETDep);
4327
4328 list->Add(fHistKaonRecEtaEDepETDep);
4329 list->Add(fHistKaonMatchEtaEDepETDep);
4330
4331 list->Add(fHistKaonRecEtaPtETDep);
4332 list->Add(fHistKaonMatchEtaPtETDep);
4333
4334 list->Add(fHistKaonRecEtaETDep);
4335 list->Add(fHistKaonMatchEtaETDep);
4336
4337 list->Add(fHistKaonRecResEET);
4338 list->Add(fHistKaonRecResPtET);
4339 list->Add(fHistKaonRecResE);
4340 list->Add(fHistKaonRecResPt);
4341 list->Add(fHistKaonRecResEDepETDep);
4342 list->Add(fHistKaonRecResPtETDep);
4343
4344 list->Add(fHistProtonEtaEET);
4345 list->Add(fHistProtonAccEtaEET);
4346 list->Add(fHistProtonRecEtaEET);
4347 list->Add(fHistProtonMatchEtaEET);
4348
4349 list->Add(fHistProtonEtaPtET);
4350 list->Add(fHistProtonAccEtaPtET);
4351 list->Add(fHistProtonRecEtaPtET);
4352 list->Add(fHistProtonMatchEtaPtET);
4353
4354 list->Add(fHistProtonEtaET);
4355 list->Add(fHistProtonAccEtaET);
4356 list->Add(fHistProtonRecEtaET);
4357 list->Add(fHistProtonMatchEtaET);
4358
4359 list->Add(fHistProtonEtaE);
4360 list->Add(fHistProtonAccEtaE);
4361 list->Add(fHistProtonRecEtaE);
4362 list->Add(fHistProtonMatchEtaE);
4363
4364 list->Add(fHistProtonEtaPt);
4365 list->Add(fHistProtonAccEtaPt);
4366 list->Add(fHistProtonRecEtaPt);
4367 list->Add(fHistProtonMatchEtaPt);
4368
4369 list->Add(fHistProtontotET);
4370 list->Add(fHistProtonAcctotET);
4371 list->Add(fHistProtonRectotET);
4372 list->Add(fHistProtonMatchtotET);
4373
4374 list->Add(fHistProtonRectotETDep);
4375 list->Add(fHistProtonMatchtotETDep);
4376
4377 list->Add(fHistProtonRecEtaEDepETDep);
4378 list->Add(fHistProtonMatchEtaEDepETDep);
4379
4380 list->Add(fHistProtonRecEtaPtETDep);
4381 list->Add(fHistProtonMatchEtaPtETDep);
4382
4383 list->Add(fHistProtonRecEtaETDep);
4384 list->Add(fHistProtonMatchEtaETDep);
4385
4386 list->Add(fHistProtonRecResEET);
4387 list->Add(fHistProtonRecResPtET);
4388 list->Add(fHistProtonRecResE);
4389 list->Add(fHistProtonRecResPt);
4390 list->Add(fHistProtonRecResEDepETDep);
4391 list->Add(fHistProtonRecResPtETDep);
4392
4393 list->Add(fHistTotChargedtotET);
4394 list->Add(fHistTotChargedAcctotET);
4395 list->Add(fHistTotChargedRectotET);
4396 list->Add(fHistTotChargedMatchtotET);
4397
4398 list->Add(fHistTotChargedRectotETDep);
4399 list->Add(fHistTotChargedMatchtotETDep);
4400
4401 list->Add(fHistNeutronEtaEET);
4402 list->Add(fHistNeutronAccEtaEET);
4403 list->Add(fHistNeutronRecEtaEET);
4404
4405 list->Add(fHistNeutronEtaPtET);
4406 list->Add(fHistNeutronAccEtaPtET);
4407 list->Add(fHistNeutronRecEtaPtET);
4408
4409 list->Add(fHistNeutronEtaET);
4410 list->Add(fHistNeutronAccEtaET);
4411 list->Add(fHistNeutronRecEtaET);
4412
4413 list->Add(fHistNeutronEtaE);
4414 list->Add(fHistNeutronAccEtaE);
4415 list->Add(fHistNeutronRecEtaE);
4416
4417 list->Add(fHistNeutronEtaPt);
4418 list->Add(fHistNeutronAccEtaPt);
4419 list->Add(fHistNeutronRecEtaPt);
4420
4421 list->Add(fHistNeutrontotET);
4422 list->Add(fHistNeutronAcctotET);
4423 list->Add(fHistNeutronRectotET);
4424 list->Add(fHistNeutronRectotETDep);
4425
4426 list->Add(fHistNeutronRecEtaEDepETDep);
4427 list->Add(fHistNeutronRecEtaETDep);
4428
4429 list->Add(fHistNeutronRecEtaPtETDep);
4430
4431 list->Add(fHistK0EtaEET);
4432 list->Add(fHistK0RecEtaEET);
4433
4434 list->Add(fHistK0EtaPtET);
4435 list->Add(fHistK0RecEtaPtET);
4436
4437 list->Add(fHistK0EtaET);
4438 list->Add(fHistK0RecEtaET);
4439
4440 list->Add(fHistK0EtaE);
4441 list->Add(fHistK0RecEtaE);
4442
4443 list->Add(fHistK0EtaPt);
4444 list->Add(fHistK0RecEtaPt);
4445
4446 list->Add(fHistK0totET);
4447 list->Add(fHistK0RectotET);
4448 list->Add(fHistK0RectotETDep);
4449
4450 list->Add(fHistK0RecEtaEDepETDep);
4451 list->Add(fHistK0RecEtaETDep);
4452
4453 list->Add(fHistK0RecEtaPtETDep);
4454
4455 list->Add(fHistLambdaEtaEET);
4456 list->Add(fHistLambdaRecEtaEET);
4457
4458 list->Add(fHistLambdaEtaPtET);
4459 list->Add(fHistLambdaRecEtaPtET);
4460
4461 list->Add(fHistLambdaEtaET);
4462 list->Add(fHistLambdaRecEtaET);
4463
4464 list->Add(fHistLambdaEtaE);
4465 list->Add(fHistLambdaRecEtaE);
4466
4467 list->Add(fHistLambdaEtaPt);
4468 list->Add(fHistLambdaRecEtaPt);
4469
4470 list->Add(fHistLambdatotET);
4471 list->Add(fHistLambdaRectotET);
4472 list->Add(fHistLambdaRectotETDep);
4473
4474 list->Add(fHistLambdaRecEtaEDepETDep);
4475 list->Add(fHistLambdaRecEtaETDep);
4476
4477 list->Add(fHistLambdaRecEtaPtETDep);
4478
4479 list->Add(fHistTotNeutraltotET);
4480 list->Add(fHistTotNeutralRectotET);
4481 list->Add(fHistTotNeutralRectotETDep);
4482
4483 list->Add(fHistTotaltotET);
4484 list->Add(fHistTotalAcctotET);
4485 list->Add(fHistTotalRectotET);
4486 list->Add(fHistTotalRectotETDep);
4487
4488 list->Add(fHistElectronFirstMother);
4489 list->Add(fHistElectronFirstMotherXY);
4490 list->Add(fHistElectronNDaughters);
4491 list->Add(fHistElectronDaughters);
4492 list->Add(fHistElectronDaughtersXY);
4493
4494 list->Add(fHistElectronFirstMotherAcc);
4495 list->Add(fHistElectronFirstMotherXYAcc);
4496 list->Add(fHistElectronNDaughtersAcc);
4497 list->Add(fHistElectronDaughtersAcc);
4498 list->Add(fHistElectronDaughtersXYAcc);
4499
4500 list->Add(fHistElectronFirstMotherRec);
4501 list->Add(fHistElectronFirstMotherXYRec);
4502 list->Add(fHistElectronNDaughtersRec);
4503 list->Add(fHistElectronDaughtersRec);
4504 list->Add(fHistElectronDaughtersXYRec);
4505
4506 list->Add(fHistNPPElectronFirstMother);
4507 list->Add(fHistNPPElectronFirstMotherXY);
4508 list->Add(fHistNPPElectronNDaughters);
4509 list->Add(fHistNPPElectronDaughters);
4510 list->Add(fHistNPPElectronDaughtersXY);
4511
4512 list->Add(fHistNPPElectronFirstMotherAcc);
4513 list->Add(fHistNPPElectronFirstMotherXYAcc);
4514 list->Add(fHistNPPElectronNDaughtersAcc);
4515 list->Add(fHistNPPElectronDaughtersAcc);
4516 list->Add(fHistNPPElectronDaughtersXYAcc);
4517
4518 list->Add(fHistNPPElectronFirstMotherRec);
4519 list->Add(fHistNPPElectronFirstMotherXYRec);
4520 list->Add(fHistNPPElectronNDaughtersRec);
4521 list->Add(fHistNPPElectronDaughtersRec);
4522 list->Add(fHistNPPElectronDaughtersXYRec);
4523
4524 list->Add(fHistGammaFirstMother);
4525 list->Add(fHistGammaFirstMotherXY);
4526 list->Add(fHistGammaNDaughters);
4527 list->Add(fHistGammaDaughters);
4528 list->Add(fHistGammaDaughtersXY);
4529 list->Add(fHistConvGammaDaughtersXY);
4530 list->Add(fHistNonConvGammaDaughtersXY);
4531
4532 list->Add(fHistGammaFirstMotherAcc);
4533 list->Add(fHistGammaFirstMotherXYAcc);
4534 list->Add(fHistGammaNDaughtersAcc);
4535 list->Add(fHistGammaDaughtersAcc);
4536 list->Add(fHistGammaDaughtersXYAcc);
4537 list->Add(fHistConvGammaDaughtersXYAcc);
4538 list->Add(fHistNonConvGammaDaughtersXYAcc);
4539
4540 list->Add(fHistGammaFirstMotherRec);
4541 list->Add(fHistGammaFirstMotherXYRec);
4542 list->Add(fHistGammaNDaughtersRec);
4543 list->Add(fHistGammaDaughtersRec);
4544 list->Add(fHistGammaDaughtersXYRec);
4545 list->Add(fHistConvGammaDaughtersXYRec);
4546 list->Add(fHistNonConvGammaDaughtersXYRec);
4547
4548 list->Add(fHistNPPGammaFirstMother);
4549 list->Add(fHistNPPGammaFirstMotherXY);
4550 list->Add(fHistNPPGammaNDaughters);
4551 list->Add(fHistNPPGammaDaughters);
4552 list->Add(fHistNPPGammaDaughtersXY);
4553
4554 list->Add(fHistNPPGammaFirstMotherAcc);
4555 list->Add(fHistNPPGammaFirstMotherXYAcc);
4556 list->Add(fHistNPPGammaNDaughtersAcc);
4557 list->Add(fHistNPPGammaDaughtersAcc);
4558 list->Add(fHistNPPGammaDaughtersXYAcc);
4559
4560 list->Add(fHistNPPGammaFirstMotherRec);
4561 list->Add(fHistNPPGammaFirstMotherXYRec);
4562 list->Add(fHistNPPGammaNDaughtersRec);
4563 list->Add(fHistNPPGammaDaughtersRec);
4564 list->Add(fHistNPPGammaDaughtersXYRec);
4565
4566 list->Add(fHistAllERecEMC);
4567 list->Add(fHistAllPtRecPtMC);
4568 list->Add(fHistElectronERecEMC);
4569 list->Add(fHistGammaERecEMC);
4570
4571 list->Add(fHistChargedRes);
4572 list->Add(fHistChargedRes2);
4573 list->Add(fHistChargedRes3);
4574 list->Add(fHistNeutralRes);
4575 list->Add(fHistElectronRes);
4576 list->Add(fHistGammaRes);
4577
4578 list->Add(fHistIsInAcc);
2e12941b 4579}
4580
4581//________________________________________________________________________
4582Bool_t AliAnalysisEmEtMonteCarlo::TrackHitsCalo(AliExternalTrackParam* extParam)
4583{//Does the track hit the calorimeter?
112b1d0c 4584 TVector3 pos(0,0,0);
2e12941b 4585
112b1d0c 4586 if (extParam)
f2adae2b 4587 {
112b1d0c 4588 if (GetTrackProjection(extParam,pos))
4589 {
4590 Bool_t inAcc = fGeoUt->IsInEMCAL(pos.X(),pos.Y(),pos.Z());
2e12941b 4591
112b1d0c 4592 //if (inAcc)
4593 // fHistIsInAcc->Fill(pos.X(),pos.Y());
2e12941b 4594
112b1d0c 4595 return inAcc;
4596 }
f2adae2b 4597 }
112b1d0c 4598 return kFALSE;
2e12941b 4599}
4600
4601//________________________________________________________________________
f2adae2b 4602//project to a EMCal radius
2e12941b 4603Bool_t AliAnalysisEmEtMonteCarlo::GetTrackProjection(AliExternalTrackParam *trackParam, TVector3 &trackPos)
4604{//Get the track projection
112b1d0c 4605 Bool_t proj = kFALSE;
4606 Double_t emcalR = fGeoUt->GetEMCGeometry()->GetIPDistance();
2e12941b 4607
112b1d0c 4608 if (trackParam) //it is constructed from TParticle
2e12941b 4609 {
112b1d0c 4610 Double_t trkPos[3] = {0};
2e12941b 4611
112b1d0c 4612 //Assume the track is a pion with mass 0.139GeV/c^2
4613 //Extrapolation step is 1cm
4614 if(!AliTrackerBase::PropagateTrackToBxByBz(trackParam, emcalR, 0.139, 1, kTRUE, 0.8) ) return proj;
2e12941b 4615
112b1d0c 4616 trackParam->GetXYZ(trkPos);
2e12941b 4617
112b1d0c 4618 trackPos.SetXYZ(trkPos[0],trkPos[1],trkPos[2]);
2e12941b 4619
112b1d0c 4620 proj = kTRUE;
2e12941b 4621 }
4622
112b1d0c 4623 return proj;
2e12941b 4624}
4625
4626//________________________________________________________________________
f2adae2b 4627//project to a cluster position
2e12941b 4628Bool_t AliAnalysisEmEtMonteCarlo::GetTrackProjection(AliEMCALTrack* emcTrack, TVector3 &trackPos, TVector3 clusPos)
4629{//get the track projection
112b1d0c 4630 Bool_t proj = kFALSE;
2e12941b 4631
112b1d0c 4632 if (emcTrack)
f2adae2b 4633 {
112b1d0c 4634 Double_t trkPos[3] = {0};
2e12941b 4635
112b1d0c 4636 emcTrack->PropagateToGlobal(clusPos.X(),clusPos.Y(),clusPos.Z(),0.,0.);
4637 emcTrack->GetXYZ(trkPos);
2e12941b 4638
112b1d0c 4639 trackPos.SetXYZ(trkPos[0],trkPos[1],trkPos[2]);
2e12941b 4640
112b1d0c 4641 proj = kTRUE;
f2adae2b 4642 }
2e12941b 4643
112b1d0c 4644 return proj;
2e12941b 4645}
4646
4647//________________________________________________________________________
4648Bool_t AliAnalysisEmEtMonteCarlo::IsInAcceptance(TParticle *part, TParticlePDG *pdg, AliExternalTrackParam* extParam)
4649{//is the track in the acceptance of the emcal?
112b1d0c 4650 if ((part) && (pdg))
f2adae2b 4651 {
112b1d0c 4652 if (TMath::Abs(pdg->Charge() - fCuts->GetMonteCarloNeutralParticle()) <1e-3 )
4653 {
4654 return fSelector->CutGeometricalAcceptance(*part);
4655 }
4656 else
4657 {
4658 return TrackHitsCalo(extParam);
4659 }
f2adae2b 4660 }
112b1d0c 4661
4662 return kFALSE;
2e12941b 4663}
4664
4665//________________________________________________________________________
4666Bool_t AliAnalysisEmEtMonteCarlo::IsInAcceptance(AliMCParticle *part)
4667{// is the track in the acceptance of the emcal?
112b1d0c 4668 if (part)
f2adae2b 4669 {
112b1d0c 4670 for (int i=0;i<part->GetNumberOfTrackReferences();i++)
4671 {
4672 AliTrackReference* aliTrkRef = part->GetTrackReference(i);
2e12941b 4673
112b1d0c 4674 if (aliTrkRef)
4675 {
4676 //if (aliTrkRef->DetectorId() == AliTrackReference::kEMCAL)
4677 // return kTRUE;
4678 if ( (aliTrkRef->DetectorId() == AliTrackReference::kEMCAL) || (fGeoUt->IsInEMCAL(aliTrkRef->X(),aliTrkRef->Y(),aliTrkRef->Z())) )
4679 return kTRUE;
4680 }
4681 }
f2adae2b 4682 }
2e12941b 4683
112b1d0c 4684 return kFALSE;
2e12941b 4685}
4686
4687//________________________________________________________________________
4688Bool_t AliAnalysisEmEtMonteCarlo::IsPrimary(AliStack *stack, Int_t iPart, TParticlePDG *pdg, Int_t iPartMom, TParticlePDG *pdgMom)
4689{//Is the track a primary track?
112b1d0c 4690 if (stack->IsPhysicalPrimary(iPart))
f2adae2b 4691 {
112b1d0c 4692 return kTRUE;
f2adae2b 4693 }
112b1d0c 4694 else if (pdg)
f2adae2b 4695 {
112b1d0c 4696 if (((pdg->PdgCode() == fgEPlusCode) || (pdg->PdgCode() == fgEMinusCode) || (pdg->PdgCode() == fgGammaCode)) &&
4697 ((IsMotherPrimaryGamma(stack,iPartMom,pdgMom)) || (IsMotherPrimaryElectron(stack,iPartMom,pdgMom))) )
4698 {
4699 return kTRUE;
4700 }
f2adae2b 4701 }
112b1d0c 4702
4703 return kFALSE;
2e12941b 4704}
4705
4706//________________________________________________________________________
4707Bool_t AliAnalysisEmEtMonteCarlo::IsMotherPrimaryGamma(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom)
4708{//Is the mother a primary gamma?
112b1d0c 4709 Int_t nStackTracks = stack->GetNtrack();
4710
4711 if (pdgMom)
f2adae2b 4712 {
112b1d0c 4713 if ((pdgMom->PdgCode() == fgGammaCode) && (iPartMom>=0) && (iPartMom < nStackTracks))
4714 {
4715 if (stack->IsPhysicalPrimary(iPartMom))
4716 return kTRUE;
4717 }
f2adae2b 4718 }
112b1d0c 4719 return kFALSE;
2e12941b 4720}
4721
4722//________________________________________________________________________
4723Bool_t AliAnalysisEmEtMonteCarlo::IsMotherPrimaryElectron(AliStack *stack, Int_t iPartMom, TParticlePDG *pdgMom)
4724{//is the mother a primary electron?
112b1d0c 4725 Int_t nStackTracks = stack->GetNtrack();
4726
4727 if (pdgMom)
f2adae2b 4728 {
112b1d0c 4729 if ((pdgMom->PdgCode() == fgEPlusCode || pdgMom->PdgCode() == fgEMinusCode) && (iPartMom>=0) && (iPartMom < nStackTracks))
4730 {
4731 if (stack->IsPhysicalPrimary(iPartMom))
4732 return kTRUE;
4733 }
f2adae2b 4734 }
112b1d0c 4735 return kFALSE;
2e12941b 4736}
4737
4738//________________________________________________________________________
4739Bool_t AliAnalysisEmEtMonteCarlo::IsGammaConversion(AliStack *stack, TParticle *part, TParticlePDG *pdg)
4740{//is the gamma converted?
112b1d0c 4741 if ((part) && (pdg))
f2adae2b 4742 {
112b1d0c 4743 if(pdg->PdgCode() == fgGammaCode)
4744 {
4745 Int_t nStackTracks = stack->GetNtrack();
4746 Int_t iPartDaughter = part->GetLastDaughter();
2e12941b 4747
112b1d0c 4748 if ((iPartDaughter>=0) && (iPartDaughter < nStackTracks))
4749 {
4750 TParticle *partDaughter = stack->Particle(iPartDaughter);
4751 if (partDaughter)
4752 {
4753 TParticlePDG *pdgDaugther = partDaughter->GetPDG(0);
4754 if (pdgDaugther)
4755 {
4756 if ( ((pdgDaugther->PdgCode() == fgEPlusCode) || (pdgDaugther->PdgCode() == fgEMinusCode)) && (!fGeoUt->IsInEMCAL(partDaughter->Vx(),partDaughter->Vy(),partDaughter->Vz())) )
4757 {
4758 //Double_t emcalR = fGeoUt->GetEMCGeometry()->GetIPDistance();
4759 //Double_t decayR = sqrt(pow(partDaughter->Vx(),2)+pow(partDaughter->Vy(),2));
2e12941b 4760
112b1d0c 4761 //if (decayR<emcalR)
4762 return kTRUE;
4763 }
4764 }
4765 }
4766 }
4767 }
f2adae2b 4768 }
112b1d0c 4769 return kFALSE;
2e12941b 4770}
112b1d0c 4771
2e12941b 4772//________________________________________________________________________
4773AliExternalTrackParam* AliAnalysisEmEtMonteCarlo::CreateExternalTrackParam(TParticle *part)
4774{//create external track param
112b1d0c 4775 // Calculate the AliExternalTrackParam content
4776 Double_t xref;
4777 Double_t alpha;
4778 Double_t param[5];
4779 Double_t covar[15];
4780
4781 // Calculate alpha: the rotation angle of the corresponding local system (TPC sector)
4782 alpha = part->Phi()*180./TMath::Pi();
4783 if (alpha<0) alpha+= 360.;
4784 if (alpha>360) alpha -= 360.;
4785
4786 Int_t sector = (Int_t)(alpha/20.);
4787 alpha = 10. + 20.*sector;
4788 alpha /= 180;
4789 alpha *= TMath::Pi();
4790
4791 // Covariance matrix: no errors, the parameters are exact
4792 for (int i=0; i<15; i++) covar[i]=0.;
4793
4794 // Get the vertex of origin and the momentum
4795 TVector3 ver(part->Vx(),part->Vy(),part->Vz());
4796 TVector3 mom(part->Px(),part->Py(),part->Pz());
4797
4798 // Rotate to the local coordinate system (TPC sector)
4799 ver.RotateZ(-alpha);
4800 mom.RotateZ(-alpha);
4801
4802 // X of the referense plane
4803 xref = ver.X();
4804
4805 Double_t charge;
4806 if (part->GetPDG(0))
4807 charge = part->GetPDG(0)->Charge();
4808 else
4809 return 0;
4810
4811 if (mom.Pt()>0)
f2adae2b 4812 {
112b1d0c 4813 param[0] = ver.Y();
4814 param[1] = ver.Z();
4815 param[2] = TMath::Sin(mom.Phi());
4816 param[3] = mom.Pz()/mom.Pt();
4817 param[4] = TMath::Sign(1/mom.Pt(),charge);
f2adae2b 4818 }
112b1d0c 4819 else
4820 return 0;
2e12941b 4821
112b1d0c 4822 // Set AliExternalTrackParam
4823 AliExternalTrackParam* extTrkParam = new AliExternalTrackParam(xref, alpha, param, covar);
2e12941b 4824
112b1d0c 4825 return extTrkParam;
2e12941b 4826}
4827
4828//________________________________________________________________________
4829Double_t AliAnalysisEmEtMonteCarlo::CalcET(TParticle *part, TParticlePDG *pdg)
4830{//Calculate Et
112b1d0c 4831 //***************
f2adae2b 4832 // calculate E_T
4833 //***************
4834 Double_t particleMassPart = 0; //The mass part in the Et calculation for this particle
112b1d0c 4835 Double_t protonMass = fgProtonMass;
4836
4837 if (pdg)
f2adae2b 4838 {
112b1d0c 4839 if (
4840 TMath::Abs(pdg->PdgCode()) == fgProtonCode ||
4841 TMath::Abs(pdg->PdgCode()) == fgNeutronCode ||
4842 TMath::Abs(pdg->PdgCode()) == fgLambdaCode ||
4843 TMath::Abs(pdg->PdgCode()) == fgXiCode ||
4844 TMath::Abs(pdg->PdgCode()) == fgXi0Code ||
4845 TMath::Abs(pdg->PdgCode()) == fgOmegaCode
4846 )
4847 {
4848 if (pdg->PdgCode() > 0) { particleMassPart = - protonMass;}
4849 if (pdg->PdgCode() < 0) { particleMassPart = protonMass;}
4850 }
4851 Double_t et = part->Energy() * TMath::Sin(part->Theta()) + particleMassPart;
4852 return et;
f2adae2b 4853 }
112b1d0c 4854 else
4855 return -1.;
2e12941b 4856}
112b1d0c 4857
2e12941b 4858//________________________________________________________________________
4859Double_t AliAnalysisEmEtMonteCarlo::CalcETDep(Double_t caloE, TParticle *part, TParticlePDG *pdg)
4860{//calculate et dependence
112b1d0c 4861 //***************
f2adae2b 4862 // calculate E_T
4863 //***************
4864 Double_t particleMassPart = 0; //The mass part in the Et calculation for this particle
112b1d0c 4865 Double_t protonMass = fgProtonMass;
2e12941b 4866
112b1d0c 4867 if (pdg)
f2adae2b 4868 {
112b1d0c 4869 if (
4870 TMath::Abs(pdg->PdgCode()) == fgProtonCode ||
4871 TMath::Abs(pdg->PdgCode()) == fgNeutronCode ||
4872 TMath::Abs(pdg->PdgCode()) == fgLambdaCode ||
4873 TMath::Abs(pdg->PdgCode()) == fgXiCode ||
4874 TMath::Abs(pdg->PdgCode()) == fgXi0Code ||
4875 TMath::Abs(pdg->PdgCode()) == fgOmegaCode
4876 )
4877 {
4878 if (pdg->PdgCode() > 0) { particleMassPart = - protonMass;}
4879 if (pdg->PdgCode() < 0) { particleMassPart = protonMass;}
4880 }
4881 Double_t et = caloE * TMath::Sin(part->Theta()) + particleMassPart;
4882 return et;
f2adae2b 4883 }
112b1d0c 4884 else
4885 return -1.;
2e12941b 4886}
4887
4888
112b1d0c 4889