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