PWG4/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         // total gamma ET
 182         // *******************
 183         TH1F *fHistTotGammatotET;//total ET distribution
 184
 185         // *******************
 186         // total electromagnetic ET
 187         // *******************
 188         TH1F *fHistTotEMtotET;//total ET distribution
 189
 190         // non-primary electromagnetic ET
 191         TH2F *fHistNPPElectronEtaEET;//ET - Eta vs E
 192         TH2F *fHistNPPElectronEtaPtET;//ET - Eta vs pt
 193         TH2F *fHistNPPElectronEtaET;// ET - Eta
 194         TH2F *fHistNPPElectronEtaE;// multiplicity - Eta vs E
 195         TH2F *fHistNPPElectronEtaPt;// multiplicity - Eta vs pt
 196         TH1F *fHistNPPElectrontotET;// total ET distribution
 197
 198         TH2F *fHistNPPGammaEtaEET;//ET - Eta vs E
 199         TH2F *fHistNPPGammaEtaPtET;//ET - Eta vs pt
 200         TH2F *fHistNPPGammaEtaET;// ET - Eta
 201         TH2F *fHistNPPGammaEtaE;// multiplicity - Eta vs E
 202         TH2F *fHistNPPGammaEtaPt;// multiplicity - Eta vs pt
 203         TH1F *fHistNPPGammatotET;// total ET distribution
 204
 205         TH1F *fHistTotNPPEMtotET;//total ET distribution
 206
 207         TH2F *fHistNPPPi0GammaEtaEET;//ET - Eta vs E
 208         TH2F *fHistNPPPi0GammaEtaPtET;//ET - Eta vs pt
 209         TH2F *fHistNPPPi0GammaEtaET;// ET - Eta
 210         TH2F *fHistNPPPi0GammaEtaE;// multiplicity - Eta vs E
 211         TH2F *fHistNPPPi0GammaEtaPt;// multiplicity - Eta vs pt
 212         TH1F *fHistNPPPi0GammatotET;// total ET distribution
 213
 214         // *******************
 215         // electron ET inside EMCal acceptance
 216         // *******************
 217         TH2F *fHistElectronAccEtaEET;//ET - Eta vs E
 218         TH2F *fHistElectronAccEtaPtET;//ET - Eta vs pt
 219         TH2F *fHistElectronAccEtaET;// ET - Eta
 220         TH2F *fHistElectronAccEtaE;// multiplicity - Eta vs E
 221         TH2F *fHistElectronAccEtaPt;// multiplicity - Eta vs pt
 222         TH1F *fHistElectronAcctotET;// total ET distribution
 223
 224         TH2F *fHistConvElectronAccEtaEET;//ET - Eta vs E
 225         TH2F *fHistConvElectronAccEtaPtET;//ET - Eta vs pt
 226         TH2F *fHistConvElectronAccEtaET;//  ET - Eta
 227         TH2F *fHistConvElectronAccEtaE;//  multiplicity - Eta vs E
 228         TH2F *fHistConvElectronAccEtaPt;//  multiplicity - Eta vs pt
 229         TH1F *fHistConvElectronAcctotET;//  total ET distribution
 230
 231         TH2F *fHistScatElectronAccEtaEET;//ET - Eta vs E
 232         TH2F *fHistScatElectronAccEtaPtET;//ET - Eta vs pt
 233         TH2F *fHistScatElectronAccEtaET;//  ET - Eta
 234         TH2F *fHistScatElectronAccEtaE;//  multiplicity - Eta vs E
 235         TH2F *fHistScatElectronAccEtaPt;//  multiplicity - Eta vs pt
 236         TH1F *fHistScatElectronAcctotET;//  total ET distribution
 237
 238         // *******************
 239         // total electron ET inside EMCal acceptance
 240         // *******************
 241         TH1F *fHistTotElectronAcctotET;//total ET distribution
 242
 243         // *******************
 244         // gamma ET inside EMCal acceptance
 245         // *******************
 246         TH2F *fHistGammaAccEtaEET;//ET - Eta vs E
 247         TH2F *fHistGammaAccEtaPtET;//ET - Eta vs pt
 248         TH2F *fHistGammaAccEtaET;//  ET - Eta
 249         TH2F *fHistGammaAccEtaE;//  multiplicity - Eta vs E
 250         TH2F *fHistGammaAccEtaPt;//  multiplicity - Eta vs pt
 251         TH1F *fHistGammaAcctotET;//  total ET distribution
 252
 253         TH2F *fHistAnnihGammaAccEtaEET;//ET - Eta vs E
 254         TH2F *fHistAnnihGammaAccEtaPtET;//ET - Eta vs pt
 255         TH2F *fHistAnnihGammaAccEtaET;//  ET - Eta
 256         TH2F *fHistAnnihGammaAccEtaE;//  multiplicity - Eta vs E
 257         TH2F *fHistAnnihGammaAccEtaPt;//  multiplicity - Eta vs pt
 258         TH1F *fHistAnnihGammaAcctotET;//  total ET distribution
 259
 260         TH2F *fHistScatGammaAccEtaEET;//ET - Eta vs E
 261         TH2F *fHistScatGammaAccEtaPtET;//ET - Eta vs pt
 262         TH2F *fHistScatGammaAccEtaET;//  ET - Eta
 263         TH2F *fHistScatGammaAccEtaE;//  multiplicity - Eta vs E
 264         TH2F *fHistScatGammaAccEtaPt;//  multiplicity - Eta vs pt
 265         TH1F *fHistScatGammaAcctotET;//  total ET distribution
 266
 267         TH2F *fHistConvGammaAccEtaEET;//ET - Eta vs E
 268         TH2F *fHistConvGammaAccEtaPtET;//ET - Eta vs pt
 269         TH2F *fHistConvGammaAccEtaET;//  ET - Eta
 270         TH2F *fHistConvGammaAccEtaE;//  multiplicity - Eta vs E
 271         TH2F *fHistConvGammaAccEtaPt;//  multiplicity - Eta vs pt
 272         TH1F *fHistConvGammaAcctotET;//  total ET distribution
 273
 274         TH2F *fHistNonConvGammaAccEtaEET;//ET - Eta vs E
 275         TH2F *fHistNonConvGammaAccEtaPtET;//ET - Eta vs pt
 276         TH2F *fHistNonConvGammaAccEtaET;//  ET - Eta
 277         TH2F *fHistNonConvGammaAccEtaE;//  multiplicity - Eta vs E
 278         TH2F *fHistNonConvGammaAccEtaPt;//  multiplicity - Eta vs pt
 279         TH1F *fHistNonConvGammaAcctotET;//  total ET distribution
 280
 281         // *******************
 282         // total gamma ET inside EMCal acceptance
 283         // *******************
 284         TH1F *fHistTotGammaAcctotET;//total ET distribution
 285
 286         // *******************
 287         // total electromagnetic ET inside EMCal acceptance
 288         // *******************
 289         TH1F *fHistTotEMAcctotET;//total ET distribution
 290
 291         // non-primary electromagnetic ET
 292         TH2F *fHistNPPElectronAccEtaEET;//ET - Eta vs E
 293         TH2F *fHistNPPElectronAccEtaPtET;//ET - Eta vs pt
 294         TH2F *fHistNPPElectronAccEtaE;// multiplicity - Eta vs E
 295         TH2F *fHistNPPElectronAccEtaPt;// multiplicity - Eta vs pt
 296
 297         TH2F *fHistNPPGammaAccEtaEET;//ET - Eta vs E
 298         TH2F *fHistNPPGammaAccEtaPtET;//ET - Eta vs pt
 299         TH2F *fHistNPPGammaAccEtaE;// multiplicity - Eta vs E
 300         TH2F *fHistNPPGammaAccEtaPt;//  multiplicity - Eta vs pt
 301
 302         // *******************
 303         // electron ET reconstructed in EMCal
 304         // *******************
 305         TH2F *fHistElectronRecEtaEET;//ET - Eta vs E
 306         TH2F *fHistElectronRecEtaPtET;//ET - Eta vs pt
 307         TH2F *fHistElectronRecEtaET;// ET - Eta
 308         TH2F *fHistElectronRecEtaE;// multiplicity - Eta vs E
 309         TH2F *fHistElectronRecEtaPt;// multiplicity - Eta vs pt
 310         TH1F *fHistElectronRectotET;// total ET distribution
 311
 312         TH2F *fHistConvElectronRecEtaEET;//ET - Eta vs E
 313         TH2F *fHistConvElectronRecEtaPtET;//ET - Eta vs pt
 314         TH2F *fHistConvElectronRecEtaET;//  ET - Eta
 315         TH2F *fHistConvElectronRecEtaE;//  multiplicity - Eta vs E
 316         TH2F *fHistConvElectronRecEtaPt;//  multiplicity - Eta vs pt
 317         TH1F *fHistConvElectronRectotET;//  total ET distribution
 318
 319         TH2F *fHistScatElectronRecEtaEET;//ET - Eta vs E
 320         TH2F *fHistScatElectronRecEtaPtET;//ET - Eta vs pt
 321         TH2F *fHistScatElectronRecEtaET;//  ET - Eta
 322         TH2F *fHistScatElectronRecEtaE;//  multiplicity - Eta vs E
 323         TH2F *fHistScatElectronRecEtaPt;//  multiplicity - Eta vs pt
 324         TH1F *fHistScatElectronRectotET;//  total ET distribution
 325
 326         // *******************
 327         // total Electron ET reconstructed in EMCal
 328         // *******************
 329         TH1F *fHistTotElectronRectotET;//total ET distribution
 330
 331         // *******************
 332         // gamma ET reconstructed in EMCal
 333         // *******************
 334         TH2F *fHistGammaRecEtaEET;//ET - Eta vs E
 335         TH2F *fHistGammaRecEtaPtET;//ET - Eta vs pt
 336         TH2F *fHistGammaRecEtaET;//  ET - Eta
 337         TH2F *fHistGammaRecEtaE;//  multiplicity - Eta vs E
 338         TH2F *fHistGammaRecEtaPt;//  multiplicity - Eta vs pt
 339         TH1F *fHistGammaRectotET;//  total ET distribution
 340
 341         TH2F *fHistAnnihGammaRecEtaEET;//ET - Eta vs E
 342         TH2F *fHistAnnihGammaRecEtaPtET;//ET - Eta vs pt
 343         TH2F *fHistAnnihGammaRecEtaET;//  ET - Eta
 344         TH2F *fHistAnnihGammaRecEtaE;//  multiplicity - Eta vs E
 345         TH2F *fHistAnnihGammaRecEtaPt;//  multiplicity - Eta vs pt
 346         TH1F *fHistAnnihGammaRectotET;//  total ET distribution
 347
 348         TH2F *fHistScatGammaRecEtaEET;//ET - Eta vs E
 349         TH2F *fHistScatGammaRecEtaPtET;//ET - Eta vs pt
 350         TH2F *fHistScatGammaRecEtaET;//  ET - Eta
 351         TH2F *fHistScatGammaRecEtaE;//  multiplicity - Eta vs E
 352         TH2F *fHistScatGammaRecEtaPt;//  multiplicity - Eta vs pt
 353         TH1F *fHistScatGammaRectotET;//  total ET distribution
 354
 355         // *******************
 356         // total gamma ET reconstructed in EMCal
 357         // *******************
 358         TH1F *fHistTotGammaRectotET;//total ET distribution
 359
 360         // *******************
 361         // total EM ET reconstructed in EMCal
 362         // *******************
 363         TH1F *fHistTotEMRectotET;//total ET distribution
 364
 365         // non-primary electromagnetic ET
 366         TH2F *fHistNPPElectronRecEtaEET;//ET - Eta vs E
 367         TH2F *fHistNPPElectronRecEtaPtET;//ET - Eta vs pt
 368         TH2F *fHistNPPElectronRecEtaET;// ET - Eta
 369         TH2F *fHistNPPElectronRecEtaE;// multiplicity - Eta vs E
 370         TH2F *fHistNPPElectronRecEtaPt;// multiplicity - Eta vs pt
 371         TH1F *fHistNPPElectronRectotET;// total ET distribution
 372
 373         TH2F *fHistNPPGammaRecEtaEET;//ET - Eta vs E
 374         TH2F *fHistNPPGammaRecEtaPtET;//ET - Eta vs pt
 375         TH2F *fHistNPPGammaRecEtaET;// ET - Eta
 376         TH2F *fHistNPPGammaRecEtaE;// multiplicity - Eta vs E
 377         TH2F *fHistNPPGammaRecEtaPt;// multiplicity - Eta vs pt
 378         TH1F *fHistNPPGammaRectotET;// total ET distribution
 379
 380         TH1F *fHistTotNPPEMRectotET;//total ET distribution
 381
 382         TH2F *fHistNPPPi0GammaRecEtaEET;//ET - Eta vs E
 383         TH2F *fHistNPPPi0GammaRecEtaPtET;//ET - Eta vs pt
 384         TH2F *fHistNPPPi0GammaRecEtaET;// ET - Eta
 385         TH2F *fHistNPPPi0GammaRecEtaE;// multiplicity - Eta vs E
 386         TH2F *fHistNPPPi0GammaRecEtaPt;// multiplicity - Eta vs pt
 387         TH1F *fHistNPPPi0GammaRectotET;// total ET distribution
 388
 389         // *******************
 390         // muon ET (+ and -)
 391         // *******************
 392         TH2F *fHistMuonEtaEET;//ET - Eta vs E
 393         TH2F *fHistMuonAccEtaEET;//ET - Eta vs E
 394         TH2F *fHistMuonRecEtaEET;//ET - Eta vs E
 395         TH2F *fHistMuonMatchEtaEET;//ET - Eta vs E
 396
 397         TH2F *fHistMuonEtaPtET;// ET - Eta vs pt
 398         TH2F *fHistMuonAccEtaPtET;// ET - Eta vs pt
 399         TH2F *fHistMuonRecEtaPtET;// ET - Eta vs pt
 400         TH2F *fHistMuonMatchEtaPtET;// ET - Eta vs pt
 401
 402         TH2F *fHistMuonEtaET;// ET - Eta
 403         TH2F *fHistMuonAccEtaET;// ET - Eta
 404         TH2F *fHistMuonRecEtaET;// ET - Eta
 405         TH2F *fHistMuonMatchEtaET;// ET - Eta
 406
 407         TH2F *fHistMuonEtaE;// multiplicity - Eta vs E
 408         TH2F *fHistMuonAccEtaE;// multiplicity - Eta vs E
 409         TH2F *fHistMuonRecEtaE;// multiplicity - Eta vs E
 410         TH2F *fHistMuonMatchEtaE;// multiplicity - Eta vs E
 411
 412         TH2F *fHistMuonEtaPt;// multiplicity - Eta vs pt
 413         TH2F *fHistMuonAccEtaPt;// multiplicity - Eta vs pt
 414         TH2F *fHistMuonRecEtaPt;// multiplicity - Eta vs pt
 415         TH2F *fHistMuonMatchEtaPt;// multiplicity - Eta vs pt
 416
 417         TH1F *fHistMuontotET;// total ET distribution
 418         TH1F *fHistMuonAcctotET;// total ET distribution
 419         TH1F *fHistMuonRectotET;// total ET distribution
 420         TH1F *fHistMuonMatchtotET;// total ET distribution
 421
 422         TH1F *fHistMuonRectotETDep;//total deposited ET distribution
 423         TH1F *fHistMuonMatchtotETDep;// total deposited ET distribution
 424
 425         TH2F *fHistMuonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 426         TH2F *fHistMuonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
 427
 428         TH2F *fHistMuonRecEtaPtETDep;// ET deposited - Eta vs pt
 429         TH2F *fHistMuonMatchEtaPtETDep;// ET deposited - Eta vs pt
 430
 431         TH2F *fHistMuonRecEtaETDep;// ET deposited - Eta
 432         TH2F *fHistMuonMatchEtaETDep;// ET deposited - Eta
 433
 434         TH2F *fHistMuonRecResEET;// ET - track matching residual vs E
 435         TH2F *fHistMuonRecResPtET;// ET - track matching residual vs pt
 436         TH2F *fHistMuonRecResE;// multiplicity - track matching residual vs E
 437         TH2F *fHistMuonRecResPt;// multiplicity - track matching residual vs pt
 438         TH2F *fHistMuonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
 439         TH2F *fHistMuonRecResPtETDep;// ET deposited - track matching residual vs pt
 440
 441         // *******************
 442         // pion ET (+ and -)
 443         // *******************
 444         TH2F *fHistPionEtaEET;//ET - Eta vs E
 445         TH2F *fHistPionAccEtaEET;//ET - Eta vs E
 446         TH2F *fHistPionRecEtaEET;//ET - Eta vs E
 447         TH2F *fHistPionMatchEtaEET;//ET - Eta vs E
 448
 449         TH2F *fHistPionEtaPtET;// ET - Eta vs pt
 450         TH2F *fHistPionAccEtaPtET;// ET - Eta vs pt
 451         TH2F *fHistPionRecEtaPtET;// ET - Eta vs pt
 452         TH2F *fHistPionMatchEtaPtET;// ET - Eta vs pt
 453
 454         TH2F *fHistPionEtaET;// ET - Eta
 455         TH2F *fHistPionAccEtaET;// ET - Eta
 456         TH2F *fHistPionRecEtaET;// ET - Eta
 457         TH2F *fHistPionMatchEtaET;// ET - Eta
 458
 459         TH2F *fHistPionEtaE;// multiplicity - Eta vs E
 460         TH2F *fHistPionAccEtaE;// multiplicity - Eta vs E
 461         TH2F *fHistPionRecEtaE;// multiplicity - Eta vs E
 462         TH2F *fHistPionMatchEtaE;// multiplicity - Eta vs E
 463
 464         TH2F *fHistPionEtaPt;// multiplicity - Eta vs pt
 465         TH2F *fHistPionAccEtaPt;// multiplicity - Eta vs pt
 466         TH2F *fHistPionRecEtaPt;// multiplicity - Eta vs pt
 467         TH2F *fHistPionMatchEtaPt;// multiplicity - Eta vs pt
 468
 469         TH1F *fHistPiontotET;// total ET distribution
 470         TH1F *fHistPionAcctotET;// total ET distribution
 471         TH1F *fHistPionRectotET;// total ET distribution
 472         TH1F *fHistPionMatchtotET;// total ET distribution
 473
 474         TH1F *fHistPionRectotETDep;// total deposited ET distribution
 475         TH1F *fHistPionMatchtotETDep;// total deposited ET distribution
 476
 477         TH2F *fHistPionRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 478         TH2F *fHistPionMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
 479
 480         TH2F *fHistPionRecEtaPtETDep;// ET deposited - Eta vs pt
 481         TH2F *fHistPionMatchEtaPtETDep;// ET deposited - Eta vs pt
 482
 483         TH2F *fHistPionRecEtaETDep;// ET deposited - Eta
 484         TH2F *fHistPionMatchEtaETDep;// ET deposited - Eta
 485
 486         TH2F *fHistPionRecResEET;// ET - track matching residual vs E
 487         TH2F *fHistPionRecResPtET;// ET - track matching residual vs pt
 488         TH2F *fHistPionRecResE;// multiplicity - track matching residual vs E
 489         TH2F *fHistPionRecResPt;// multiplicity - track matching residual vs pt
 490         TH2F *fHistPionRecResEDepETDep;// ET deposited - track matching residual vs E deposited
 491         TH2F *fHistPionRecResPtETDep;// ET deposited - track matching residual vs pt
 492
 493         // *******************
 494         // charged kaon (+ and -) ET
 495         // *******************
 496         TH2F *fHistKaonEtaEET;//ET - Eta vs E
 497         TH2F *fHistKaonAccEtaEET;//ET - Eta vs E
 498         TH2F *fHistKaonRecEtaEET;//ET - Eta vs E
 499         TH2F *fHistKaonMatchEtaEET;//ET - Eta vs E
 500
 501         TH2F *fHistKaonEtaPtET;// ET - Eta vs pt
 502         TH2F *fHistKaonAccEtaPtET;// ET - Eta vs pt
 503         TH2F *fHistKaonRecEtaPtET;// ET - Eta vs pt
 504         TH2F *fHistKaonMatchEtaPtET;// ET - Eta vs pt
 505
 506         TH2F *fHistKaonEtaET;// ET - Eta
 507         TH2F *fHistKaonAccEtaET;// ET - Eta
 508         TH2F *fHistKaonRecEtaET;// ET - Eta
 509         TH2F *fHistKaonMatchEtaET;// ET - Eta
 510
 511         TH2F *fHistKaonEtaE;// multiplicity - Eta vs E
 512         TH2F *fHistKaonAccEtaE;// multiplicity - Eta vs E
 513         TH2F *fHistKaonRecEtaE;// multiplicity - Eta vs E
 514         TH2F *fHistKaonMatchEtaE;// multiplicity - Eta vs E
 515
 516         TH2F *fHistKaonEtaPt;// multiplicity - Eta vs pt
 517         TH2F *fHistKaonAccEtaPt;// multiplicity - Eta vs pt
 518         TH2F *fHistKaonRecEtaPt;// multiplicity - Eta vs pt
 519         TH2F *fHistKaonMatchEtaPt;// multiplicity - Eta vs pt
 520
 521         TH1F *fHistKaontotET;// total ET distribution
 522         TH1F *fHistKaonAcctotET;// total ET distribution
 523         TH1F *fHistKaonRectotET;// total ET distribution
 524         TH1F *fHistKaonMatchtotET;// total ET distribution
 525
 526         TH1F *fHistKaonRectotETDep;// total deposited ET distribution
 527         TH1F *fHistKaonMatchtotETDep;// total deposited ET distribution
 528
 529         TH2F *fHistKaonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 530         TH2F *fHistKaonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
 531
 532         TH2F *fHistKaonRecEtaPtETDep;// ET deposited - Eta vs pt
 533         TH2F *fHistKaonMatchEtaPtETDep;// ET deposited - Eta vs pt
 534
 535         TH2F *fHistKaonRecEtaETDep;// ET deposited - Eta
 536         TH2F *fHistKaonMatchEtaETDep;// ET deposited - Eta
 537
 538         TH2F *fHistKaonRecResEET;// ET - track matching residual vs E
 539         TH2F *fHistKaonRecResPtET;// ET - track matching residual vs pt
 540         TH2F *fHistKaonRecResE;// multiplicity - track matching residual vs E
 541         TH2F *fHistKaonRecResPt;// multiplicity - track matching residual vs pt
 542         TH2F *fHistKaonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
 543         TH2F *fHistKaonRecResPtETDep;// ET deposited - track matching residual vs pt
 544
 545         // *******************
 546         // proton (anti) ET
 547         // *******************
 548         TH2F *fHistProtonEtaEET;//ET - Eta vs E
 549         TH2F *fHistProtonAccEtaEET;//ET - Eta vs E
 550         TH2F *fHistProtonRecEtaEET;//ET - Eta vs E
 551         TH2F *fHistProtonMatchEtaEET;//ET - Eta vs E
 552
 553         TH2F *fHistProtonEtaPtET;// ET - Eta vs pt
 554         TH2F *fHistProtonAccEtaPtET;// ET - Eta vs pt
 555         TH2F *fHistProtonRecEtaPtET;// ET - Eta vs pt
 556         TH2F *fHistProtonMatchEtaPtET;// ET - Eta vs pt
 557
 558         TH2F *fHistProtonEtaET;// ET - Eta
 559         TH2F *fHistProtonAccEtaET;// ET - Eta
 560         TH2F *fHistProtonRecEtaET;// ET - Eta
 561         TH2F *fHistProtonMatchEtaET;// ET - Eta
 562
 563         TH2F *fHistProtonEtaE;// multiplicity - Eta vs E
 564         TH2F *fHistProtonAccEtaE;// multiplicity - Eta vs E
 565         TH2F *fHistProtonRecEtaE;// multiplicity - Eta vs E
 566         TH2F *fHistProtonMatchEtaE;// multiplicity - Eta vs E
 567
 568         TH2F *fHistProtonEtaPt;// multiplicity - Eta vs pt
 569         TH2F *fHistProtonAccEtaPt;// multiplicity - Eta vs pt
 570         TH2F *fHistProtonRecEtaPt;// multiplicity - Eta vs pt
 571         TH2F *fHistProtonMatchEtaPt;// multiplicity - Eta vs pt
 572
 573         TH1F *fHistProtontotET;// total ET distribution
 574         TH1F *fHistProtonAcctotET;// total ET distribution
 575         TH1F *fHistProtonRectotET;// total ET distribution
 576         TH1F *fHistProtonMatchtotET;// total ET distribution
 577
 578         TH1F *fHistProtonRectotETDep;// total deposited ET distribution
 579         TH1F *fHistProtonMatchtotETDep;// total deposited ET distribution
 580
 581         TH2F *fHistProtonRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 582         TH2F *fHistProtonMatchEtaEDepETDep;// ET deposited - Eta vs E deposited
 583
 584         TH2F *fHistProtonRecEtaPtETDep;// ET deposited - Eta vs pt
 585         TH2F *fHistProtonMatchEtaPtETDep;// ET deposited - Eta vs pt
 586
 587         TH2F *fHistProtonRecEtaETDep;// ET deposited - Eta
 588         TH2F *fHistProtonMatchEtaETDep;// ET deposited - Eta
 589
 590         TH2F *fHistProtonRecResEET;// ET - track matching residual vs E
 591         TH2F *fHistProtonRecResPtET;// ET - track matching residual vs pt
 592         TH2F *fHistProtonRecResE;// multiplicity - track matching residual vs E
 593         TH2F *fHistProtonRecResPt;// multiplicity - track matching residual vs pt
 594         TH2F *fHistProtonRecResEDepETDep;// ET deposited - track matching residual vs E deposited
 595         TH2F *fHistProtonRecResPtETDep;// ET deposited - track matching residual vs pt
 596
 597         // *******************
 598         // total charged ET
 599         // *******************
 600         TH1F *fHistTotChargedtotET;//total ET distribution
 601         TH1F *fHistTotChargedAcctotET;//total ET distribution
 602         TH1F *fHistTotChargedRectotET;//total ET distribution
 603         TH1F *fHistTotChargedRectotETDep;//total deposited ET distribution
 604         TH1F *fHistTotChargedMatchtotET;//total ET distribution
 605         TH1F *fHistTotChargedMatchtotETDep;//total deposited ET distribution
 606
 607         // *******************
 608         // neutron (anti) ET
 609         // *******************
 610         TH2F *fHistNeutronEtaEET;//ET - Eta vs E
 611         TH2F *fHistNeutronAccEtaEET;//ET - Eta vs E
 612         TH2F *fHistNeutronRecEtaEET;//ET - Eta vs E
 613
 614         TH2F *fHistNeutronEtaPtET;// ET - Eta vs pt
 615         TH2F *fHistNeutronAccEtaPtET;// ET - Eta vs pt
 616         TH2F *fHistNeutronRecEtaPtET;// ET - Eta vs pt
 617
 618         TH2F *fHistNeutronEtaET;// ET - Eta
 619         TH2F *fHistNeutronAccEtaET;// ET - Eta
 620         TH2F *fHistNeutronRecEtaET;// ET - Eta
 621
 622         TH2F *fHistNeutronEtaE;// multiplicity - Eta vs E
 623         TH2F *fHistNeutronAccEtaE;// multiplicity - Eta vs E
 624         TH2F *fHistNeutronRecEtaE;// multiplicity - Eta vs E
 625
 626         TH2F *fHistNeutronEtaPt;// multiplicity - Eta vs pt
 627         TH2F *fHistNeutronAccEtaPt;// multiplicity - Eta vs pt
 628         TH2F *fHistNeutronRecEtaPt;// multiplicity - Eta vs pt
 629
 630         TH1F *fHistNeutrontotET;// total ET distribution
 631         TH1F *fHistNeutronAcctotET;// total ET distribution
 632         TH1F *fHistNeutronRectotET;// total ET distribution
 633         TH1F *fHistNeutronRectotETDep;// total deposited ET distribution
 634
 635         TH2F *fHistNeutronRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 636         TH2F *fHistNeutronRecEtaETDep;// ET deposited - Eta
 637
 638         TH2F *fHistNeutronRecEtaPtETDep;// ET deposited - Eta vs pt
 639
 640         // *******************
 641         // neutral kaon ET
 642         // *******************
 643         TH2F *fHistK0EtaEET;//ET - Eta vs E
 644         TH2F *fHistK0RecEtaEET;//ET - Eta vs E
 645
 646         TH2F *fHistK0EtaPtET;// ET - Eta vs pt
 647         TH2F *fHistK0RecEtaPtET;// ET - Eta vs pt
 648
 649         TH2F *fHistK0EtaET;// ET - Eta
 650         TH2F *fHistK0RecEtaET;// ET - Eta
 651
 652         TH2F *fHistK0EtaE;// multiplicity - Eta vs E
 653         TH2F *fHistK0RecEtaE;// multiplicity - Eta vs E
 654
 655         TH2F *fHistK0EtaPt;// multiplicity - Eta vs pt
 656         TH2F *fHistK0RecEtaPt;// multiplicity - Eta vs pt
 657
 658         TH1F *fHistK0totET;// total ET distribution
 659         TH1F *fHistK0RectotET;// total ET distribution
 660
 661         TH1F *fHistK0RectotETDep;// total deposited ET distribution
 662
 663         TH2F *fHistK0RecEtaEDepETDep;// ET deposited - Eta vs E deposited
 664         TH2F *fHistK0RecEtaETDep;// ET deposited - Eta
 665
 666         TH2F *fHistK0RecEtaPtETDep;// ET deposited - Eta vs pt
 667
 668         // *******************
 669         // Lambda(anti) ET
 670         // *******************
 671         TH2F *fHistLambdaEtaEET;//ET - Eta vs E
 672         TH2F *fHistLambdaRecEtaEET;//ET - Eta vs E
 673
 674         TH2F *fHistLambdaEtaPtET;// ET - Eta vs pt
 675         TH2F *fHistLambdaRecEtaPtET;// ET - Eta vs pt
 676
 677         TH2F *fHistLambdaEtaET;// ET - Eta
 678         TH2F *fHistLambdaRecEtaET;// ET - Eta
 679
 680         TH2F *fHistLambdaEtaE;// multiplicity - Eta vs E
 681         TH2F *fHistLambdaRecEtaE;// multiplicity - Eta vs E
 682
 683         TH2F *fHistLambdaEtaPt;// multiplicity - Eta vs pt
 684         TH2F *fHistLambdaRecEtaPt;// multiplicity - Eta vs pt
 685
 686         TH1F *fHistLambdatotET;// total ET distribution
 687         TH1F *fHistLambdaRectotET;// total ET distribution
 688
 689         TH1F *fHistLambdaRectotETDep;// total deposited ET distribution
 690
 691         TH2F *fHistLambdaRecEtaEDepETDep;// ET deposited - Eta vs E deposited
 692         TH2F *fHistLambdaRecEtaETDep;// ET deposited - Eta
 693
 694         TH2F *fHistLambdaRecEtaPtETDep;// ET deposited - Eta vs pt
 695
 696         // *******************
 697         // total neutral ET
 698         // *******************
 699         TH1F *fHistTotNeutraltotET;//total ET distribution
 700         TH1F *fHistTotNeutralRectotET;//total ET distribution
 701         TH1F *fHistTotNeutralRectotETDep;//total deposited ET distribution
 702
 703         // *******************
 704         // total ET
 705         // *******************
 706         TH1F *fHistTotaltotET;//total ET distribution
 707         TH1F *fHistTotalAcctotET;//total ET distribution
 708         TH1F *fHistTotalRectotET;//total ET distribution
 709         TH1F *fHistTotalRectotETDep;//total deposited ET distribution
 710
 711         // *******************
 712         // some checks
 713         // *******************
 714
 715         // check produced electrons
 716         TH1F *fHistElectronFirstMother;// first mother ID
 717         TH2F *fHistElectronFirstMotherXY;// first mother XY position
 718         TH1F *fHistElectronNDaughters;// number of daughters
 719         TH1F *fHistElectronDaughters;// daughters ID
 720         TH2F *fHistElectronDaughtersXY;// daughters XY position
 721
 722         TH1F *fHistElectronFirstMotherAcc;// first mother ID
 723         TH2F *fHistElectronFirstMotherXYAcc;// first mother XY position
 724         TH1F *fHistElectronNDaughtersAcc;// number of daughters
 725         TH1F *fHistElectronDaughtersAcc;// daughters ID
 726         TH2F *fHistElectronDaughtersXYAcc;// daughters XY position
 727
 728         TH1F *fHistElectronFirstMotherRec;// first mother ID
 729         TH2F *fHistElectronFirstMotherXYRec;//  first mother XY position
 730         TH1F *fHistElectronNDaughtersRec;// number of daughters
 731         TH1F *fHistElectronDaughtersRec;// daughters ID
 732         TH2F *fHistElectronDaughtersXYRec;// daughters XY position
 733
 734         TH1F *fHistNPPElectronFirstMother;//  first mother ID
 735         TH2F *fHistNPPElectronFirstMotherXY;// first mother XY position
 736         TH1F *fHistNPPElectronNDaughters;// number of daughters
 737         TH1F *fHistNPPElectronDaughters;// daughters ID
 738         TH2F *fHistNPPElectronDaughtersXY;// daughters XY position
 739
 740         TH1F *fHistNPPElectronFirstMotherAcc;//  first mother ID
 741         TH2F *fHistNPPElectronFirstMotherXYAcc;//  first mother XY position
 742         TH1F *fHistNPPElectronNDaughtersAcc;// number of daughters
 743         TH1F *fHistNPPElectronDaughtersAcc;// daughters ID
 744         TH2F *fHistNPPElectronDaughtersXYAcc;// daughters XY position
 745
 746         TH1F *fHistNPPElectronFirstMotherRec;// first mother ID
 747         TH2F *fHistNPPElectronFirstMotherXYRec;//  first mother XY position
 748         TH1F *fHistNPPElectronNDaughtersRec;// number of daughters
 749         TH1F *fHistNPPElectronDaughtersRec;// daughters ID
 750         TH2F *fHistNPPElectronDaughtersXYRec;// daughters XY position
 751
 752         // check produced gammas
 753         TH1F *fHistGammaFirstMother;// first mother ID
 754         TH2F *fHistGammaFirstMotherXY;// first mother XY position
 755         TH1F *fHistGammaNDaughters;// number of daughters
 756         TH1F *fHistGammaDaughters;// daughters ID
 757         TH2F *fHistGammaDaughtersXY;// daughters XY position
 758         TH2F *fHistConvGammaDaughtersXY;// daughters XY position
 759         TH2F *fHistNonConvGammaDaughtersXY;// daughters XY position
 760
 761         TH1F *fHistGammaFirstMotherAcc;// first mother ID
 762         TH2F *fHistGammaFirstMotherXYAcc;//  first mother XY position
 763         TH1F *fHistGammaNDaughtersAcc;// number of daughters
 764         TH1F *fHistGammaDaughtersAcc;// daughters ID
 765         TH2F *fHistGammaDaughtersXYAcc;// daughters XY position
 766         TH2F *fHistConvGammaDaughtersXYAcc;// daughters XY position
 767         TH2F *fHistNonConvGammaDaughtersXYAcc;// daughters XY position
 768
 769         TH1F *fHistGammaFirstMotherRec;// first mother ID
 770         TH2F *fHistGammaFirstMotherXYRec;//  first mother XY position
 771         TH1F *fHistGammaNDaughtersRec;// number of daughters
 772         TH1F *fHistGammaDaughtersRec;// daughters ID
 773         TH2F *fHistGammaDaughtersXYRec;// daughters XY position
 774         TH2F *fHistConvGammaDaughtersXYRec;// daughters XY position
 775         TH2F *fHistNonConvGammaDaughtersXYRec;// daughters XY position
 776
 777         TH1F *fHistNPPGammaFirstMother;// first mother ID
 778         TH2F *fHistNPPGammaFirstMotherXY;// first mother XY position
 779         TH1F *fHistNPPGammaNDaughters;// number of daughters
 780         TH1F *fHistNPPGammaDaughters;// daughters ID
 781         TH2F *fHistNPPGammaDaughtersXY;// daughters XY position
 782
 783         TH1F *fHistNPPGammaFirstMotherAcc;// first mother ID
 784         TH2F *fHistNPPGammaFirstMotherXYAcc;//  first mother XY position
 785         TH1F *fHistNPPGammaNDaughtersAcc;// number of daughters
 786         TH1F *fHistNPPGammaDaughtersAcc;// daughters ID
 787         TH2F *fHistNPPGammaDaughtersXYAcc;// daughters XY position
 788
 789         TH1F *fHistNPPGammaFirstMotherRec;// first mother ID
 790         TH2F *fHistNPPGammaFirstMotherXYRec;//  first mother XY position
 791         TH1F *fHistNPPGammaNDaughtersRec;// number of daughters
 792         TH1F *fHistNPPGammaDaughtersRec;// daughters ID
 793         TH2F *fHistNPPGammaDaughtersXYRec;// daughters XY position
 794
 795         //check projections
 796         TH2F *fHistAllERecEMC;// E reconstructed vs E MC
 797         TH2F *fHistAllPtRecPtMC;// pt reconstructed vs pt MC
 798         TH2F *fHistElectronERecEMC;// E reconstructed vs E MC
 799         TH2F *fHistGammaERecEMC;// E reconstructed vs E MC
 800
 801         TH2F *fHistChargedRes;// charged particle track matching residual
 802         TH2F *fHistChargedRes2;// charged particle track matching residual
 803         TH2F *fHistChargedRes3;// charged particle track matching residual
 804         TH2F *fHistNeutralRes;// neutral particle track matching residual
 805         TH2F *fHistElectronRes;// electron track matching residual
 806         TH2F *fHistGammaRes;// gamma track matching residual
 807
 808         TH2F *fHistIsInAcc;// EMCal acceptance check
 809
 810  private:
 811
 812   //Declare it private to avoid compilation warning
 813     AliAnalysisEmEtMonteCarlo & operator = (const AliAnalysisEmEtMonteCarlo & g) ;//cpy assignment
 814     AliAnalysisEmEtMonteCarlo(const AliAnalysisEmEtMonteCarlo & g) ; // cpy ctor
 815     ClassDef(AliAnalysisEmEtMonteCarlo, 1);
 816 };
 817
 818 #endif //ALIANALYSISEMETMONTECARLO_H