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