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