PWGLF/totEt/AliAnalysisEtMonteCarlo.h

   1 #ifndef ALIANALYSISETMONTECARLO_H
   2 #define ALIANALYSISETMONTECARLO_H
   3 //_________________________________________________________________________
   4 //  Utility Class for transverse energy studies
   5 //  Base class for MC analysis
   6 //  - MC output
   7 //
   8 //*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
   9 //_________________________________________________________________________
  10 #include "AliAnalysisEt.h"
  11 class TParticle;
  12 class TH3F;
  13 class TH2I;
  14 class AliPHOSGeoUtils;
  15 class AliPHOSGeometry;
  16 class AliStack;
  17 //class AliMCEvent;
  18 //class AliESDEvent;
  19
  20 class AliAnalysisEtMonteCarlo : public AliAnalysisEt
  21 {
  22
  23 public:
  24
  25     AliAnalysisEtMonteCarlo();
  26     virtual ~AliAnalysisEtMonteCarlo();
  27
  28     virtual Int_t AnalyseEvent(AliVEvent* event);
  29     virtual Int_t AnalyseEvent(AliVEvent* event, AliVEvent* event2);
  30     //virtual Int_t AnalyseEvent(AliMCEvent* event, AliESDEvent* event2);
  31
  32     virtual void Init();
  33     virtual void ResetEventValues();
  34     virtual void CreateHistograms();
  35     virtual void FillOutputList(TList* list);
  36
  37     virtual void FillHistograms();
  38
  39     void CalcTrackMatchVsMult(){fCalcTrackMatchVsMult = kTRUE;}
  40     void CalcForKaonCorrection(){fCalcForKaonCorrection = kTRUE;}
  41     void IsData(){fIsMC = kFALSE;}
  42
  43     void SetNumberOfChargedHadronsMatched(Float_t val){nChargedHadronsMeasured = val;}
  44     void SetTotalNumberOfChargedHadrons(Float_t val){nChargedHadronsTotal = val;}
  45     Float_t GetNumberOfChargedHadronsMatched(){return nChargedHadronsMeasured;}
  46     Float_t GetTotalNumberOfChargedHadrons(){return nChargedHadronsTotal;}
  47
  48 protected:
  49
  50     virtual bool TrackHitsCalorimeter(TParticle *part, Double_t magField=0.5);
  51
  52
  53     Int_t GetPrimMother(Int_t partIdx, AliStack *stack);
  54
  55     Int_t GetK0InFamily(Int_t partIdx, AliStack *stack);
  56
  57     Int_t PrintFamilyTree(Int_t partIdx, AliStack *stack);
  58     Int_t PrintMothers(Int_t partIdx, AliStack *stack, Int_t gen);
  59     Int_t PrintFamilyTreeShort(Int_t partIdx, AliStack *stack);
  60     Int_t PrintMothersShort(Int_t partIdx, AliStack *stack, Int_t gen);
  61
  62
  63
  64 protected:
  65
  66     Float_t nChargedHadronsMeasured;
  67     Float_t nChargedHadronsTotal;
  68
  69     Bool_t fIsMC;//if we are running over data, we still need this object to exist but we don't want to do anything.
  70
  71     Double_t fImpactParameter; // b(fm), for Hijing; 0 otherwise
  72     Int_t fNcoll; // Ncoll, for Hijing; 1 otherwise
  73     Int_t fNpart; // Ncoll, for Hijing; 2 otherwise
  74
  75     TTree *fPrimaryTree; // Tree holding info on primaries
  76
  77     Double_t fTotEtWithSecondaryRemoved; // enter comment here
  78     Double_t fTotEtSecondaryFromEmEtPrimary; // enter comment here
  79     Double_t fTotEtSecondary; // enter comment here
  80
  81     Int_t fPrimaryCode; // enter comment here
  82     Int_t fPrimaryCharge; // enter comment here
  83
  84     Double_t fPrimaryE; // enter comment here
  85     Double_t fPrimaryEt; // enter comment here
  86
  87     Double_t fPrimaryPx; // enter comment here
  88     Double_t fPrimaryPy; // enter comment here
  89     Double_t fPrimaryPz; // enter comment here
  90
  91     Double_t fPrimaryVx; // enter comment here
  92     Double_t fPrimaryVy; // enter comment here
  93     Double_t fPrimaryVz; // enter comment here
  94
  95     Bool_t fPrimaryAccepted; // enter comment here
  96     Bool_t fPrimaryMatched;
  97     Int_t fDepositedCode; // enter comment here Double_t fDepositedEt; // enter comment here
  98     Double_t fDepositedE;
  99     Double_t fDepositedEt;
 100     Int_t fDepositedCharge; // enter comment here
 101
 102     Double_t fDepositedVx; // enter comment here
 103     Double_t fDepositedVy; // enter comment here
 104     Double_t fDepositedVz; // enter comment here
 105
 106     Bool_t fSecondary;
 107
 108     Double_t fReconstructedE;
 109     Double_t fReconstructedEt;
 110
 111     Double_t fTotPx;
 112     Double_t fTotPy;
 113     Double_t fTotPz;
 114
 115
 116     Int_t fClusterMult;
 117
 118     TH3F *fHistDecayVertexNonRemovedCharged; // Decay vertex for non-removed charged particles
 119     TH3F *fHistDecayVertexRemovedCharged; // Decay vertex for non-removed charged particles
 120     TH3F *fHistDecayVertexNonRemovedNeutral; // Decay vertex for non-removed charged particles
 121     TH3F *fHistDecayVertexRemovedNeutral; // Decay vertex for non-removed charged particles
 122
 123     TH2F *fHistRemovedOrNot; // If charged/neutral particles were removed or not
 124
 125     TH2F *fHistEtNonRemovedProtons; // enter comment here
 126     TH2F *fHistEtNonRemovedAntiProtons; // enter comment here
 127     TH2F *fHistEtNonRemovedPiPlus; // enter comment here
 128     TH2F *fHistEtNonRemovedPiMinus; // enter comment here
 129     TH2F *fHistEtNonRemovedKaonPlus; // enter comment here
 130     TH2F *fHistEtNonRemovedKaonMinus; // enter comment here
 131     TH2F *fHistEtNonRemovedK0s; // enter comment here
 132     TH2F *fHistEtNonRemovedK0L; // enter comment here
 133     TH2F *fHistEtNonRemovedLambdas; // enter comment here
 134     TH2F *fHistEtNonRemovedElectrons; // enter comment here
 135     TH2F *fHistEtNonRemovedPositrons; // enter comment here
 136     TH2F *fHistEtNonRemovedMuPlus; // enter comment here
 137     TH2F *fHistEtNonRemovedMuMinus; // enter comment here
 138     TH2F *fHistEtNonRemovedNeutrons; // enter comment here
 139     TH2F *fHistEtNonRemovedAntiNeutrons; // enter comment here
 140     TH2F *fHistEtNonRemovedGammas; // enter comment here
 141     TH2F *fHistEtNonRemovedGammasFromPi0; // enter comment here
 142
 143     TH2F *fHistEtRemovedGammas; // enter comment here
 144     TH2F *fHistEtRemovedNeutrons; // enter comment here
 145     TH2F *fHistEtRemovedAntiNeutrons; // enter comment here
 146
 147     TH2F *fHistEtRemovedCharged; // enter comment here
 148     TH2F *fHistEtRemovedNeutrals; // enter comment here
 149
 150     TH2F *fHistEtNonRemovedCharged; // enter comment here
 151     TH2F *fHistEtNonRemovedNeutrals; // enter comment here
 152
 153     TH2F *fHistMultNonRemovedProtons; // enter comment here
 154     TH2F *fHistMultNonRemovedAntiProtons; // enter comment here
 155     TH2F *fHistMultNonRemovedPiPlus; // enter comment here
 156     TH2F *fHistMultNonRemovedPiMinus; // enter comment here
 157     TH2F *fHistMultNonRemovedKaonPlus; // enter comment here
 158     TH2F *fHistMultNonRemovedKaonMinus; // enter comment here
 159     TH2F *fHistMultNonRemovedK0s; // enter comment here
 160     TH2F *fHistMultNonRemovedK0L; // enter comment here
 161     TH2F *fHistMultNonRemovedLambdas; // enter comment here
 162     TH2F *fHistMultNonRemovedElectrons; // enter comment here
 163     TH2F *fHistMultNonRemovedPositrons; // enter comment here
 164     TH2F *fHistMultNonRemovedMuPlus; // enter comment here
 165     TH2F *fHistMultNonRemovedMuMinus; // enter comment here
 166     TH2F *fHistMultNonRemovedNeutrons; // enter comment here
 167     TH2F *fHistMultNonRemovedAntiNeutrons; // enter comment here
 168     TH2F *fHistMultNonRemovedGammas; // enter comment here
 169
 170     TH2F *fHistMultRemovedGammas; // enter comment here
 171     TH2F *fHistMultRemovedNeutrons; // enter comment here
 172     TH2F *fHistMultRemovedAntiNeutrons; // enter comment here
 173
 174     TH2F *fHistMultRemovedCharged; // enter comment here
 175     TH2F *fHistMultRemovedNeutrals; // enter comment here
 176
 177     TH2F *fHistMultNonRemovedCharged; // enter comment here
 178     TH2F *fHistMultNonRemovedNeutrals; // enter comment here
 179
 180     TH2F *fHistTrackMultvsNonRemovedCharged; // enter comment here
 181     TH2F *fHistTrackMultvsNonRemovedNeutral; // enter comment here
 182     TH2F *fHistTrackMultvsRemovedGamma; // enter comment here
 183
 184     TH2F *fHistClusterMultvsNonRemovedCharged; // enter comment here
 185     TH2F *fHistClusterMultvsNonRemovedNeutral; // enter comment here
 186     TH2F *fHistClusterMultvsRemovedGamma; // enter comment here
 187
 188     TH2F *fHistMultvsNonRemovedChargedE; // enter comment here
 189     TH2F *fHistMultvsNonRemovedNeutralE; // enter comment here
 190     TH2F *fHistMultvsRemovedGammaE; // enter comment here
 191
 192     Bool_t fCalcForKaonCorrection;//turns on and off creation of kaon correction histograms
 193     TH3F *fHistK0EDepositsVsPtInAcceptance; // enter comment here
 194     TH3F *fHistK0EGammaVsPtInAcceptance; // enter comment here
 195     TH3F *fHistK0EDepositsVsPtOutOfAcceptance; // enter comment here
 196     TH3F *fHistK0EGammaVsPtOutOfAcceptance; // enter comment here
 197     TH1F *fHistSimKaonsInAcceptance;// enter comment here
 198     TH1F *fHistSimK0SInAcceptance;// enter comment here
 199     TH1F *fHistSimKPlusInAcceptance;// enter comment here
 200     TH1F *fHistSimKMinusInAcceptance;// enter comment here
 201     TH1F *fHistSimK0LInAcceptance;// enter comment here
 202     TH1F *fHistSimKaonsOutOfAcceptance;// enter comment here
 203     TH1F *fHistSimKaonsInAcceptanceWithDepositsPrimaries;// enter comment here
 204     TH1F *fHistSimKaonsOutOfAcceptanceWithDepositsSecondaries;// enter comment here
 205     TH1F *fHistSimKaonsOutOfAcceptanceWithDepositsPrimaries;// enter comment here
 206
 207     Float_t fEtNonRemovedProtons; // enter comment here
 208     Float_t fEtNonRemovedAntiProtons; // enter comment here
 209     Float_t fEtNonRemovedPiPlus; // enter comment here
 210     Float_t fEtNonRemovedPiMinus; // enter comment here
 211     Float_t fEtNonRemovedKaonPlus; // enter comment here
 212     Float_t fEtNonRemovedKaonMinus; // enter comment here
 213     Float_t fEtNonRemovedK0S; // enter comment here
 214     Float_t fEtNonRemovedK0L; // enter comment here
 215     Float_t fEtNonRemovedLambdas; // enter comment here
 216     Float_t fEtNonRemovedElectrons; // enter comment here
 217     Float_t fEtNonRemovedPositrons; // enter comment here
 218     Float_t fEtNonRemovedMuMinus; // enter comment here
 219     Float_t fEtNonRemovedMuPlus; // enter comment here
 220     Float_t fEtNonRemovedGammas; // enter comment here
 221     Float_t fEtNonRemovedGammasFromPi0; // enter comment here
 222     Float_t fEtNonRemovedNeutrons; // enter comment here
 223     Float_t fEtNonRemovedAntiNeutrons; // enter comment here
 224
 225     Float_t fEtRemovedProtons; // enter comment here
 226     Float_t fEtRemovedAntiProtons; // enter comment here
 227     Float_t fEtRemovedPiPlus; // enter comment here
 228     Float_t fEtRemovedPiMinus; // enter comment here
 229     Float_t fEtRemovedKaonPlus; // enter comment here
 230     Float_t fEtRemovedKaonMinus; // enter comment here
 231     Float_t fEtRemovedK0s; // enter comment here
 232     Float_t fEtRemovedK0L; // enter comment here
 233     Float_t fEtRemovedLambdas; // enter comment here
 234     Float_t fEtRemovedElectrons; // enter comment here
 235     Float_t fEtRemovedPositrons; // enter comment here
 236     Float_t fEtRemovedMuMinus; // enter comment here
 237     Float_t fEtRemovedMuPlus; // enter comment here
 238
 239     Float_t fEtRemovedGammasFromPi0; // enter comment here
 240     Float_t fEtRemovedGammas; // enter comment here
 241     Float_t fEtRemovedNeutrons; // enter comment here
 242     Float_t fEtRemovedAntiNeutrons; // enter comment here
 243
 244     Int_t fMultNonRemovedProtons; // enter comment here
 245     Int_t fMultNonRemovedAntiProtons; // enter comment here
 246     Int_t fMultNonRemovedPiPlus; // enter comment here
 247     Int_t fMultNonRemovedPiMinus; // enter comment here
 248     Int_t fMultNonRemovedKaonPlus; // enter comment here
 249     Int_t fMultNonRemovedKaonMinus; // enter comment here
 250     Int_t fMultNonRemovedK0s; // enter comment here
 251     Int_t fMultNonRemovedK0L; // enter comment here
 252     Int_t fMultNonRemovedLambdas; // enter comment here
 253     Int_t fMultNonRemovedElectrons; // enter comment here
 254     Int_t fMultNonRemovedPositrons; // enter comment here
 255     Int_t fMultNonRemovedMuMinus; // enter comment here
 256     Int_t fMultNonRemovedMuPlus; // enter comment here
 257     Int_t fMultNonRemovedGammas; // enter comment here
 258     Int_t fMultNonRemovedNeutrons; // enter comment here
 259     Int_t fMultNonRemovedAntiNeutrons; // enter comment here
 260
 261     Int_t fMultRemovedProtons; // enter comment here
 262     Int_t fMultRemovedAntiProtons; // enter comment here
 263     Int_t fMultRemovedPiPlus; // enter comment here
 264     Int_t fMultRemovedPiMinus; // enter comment here
 265     Int_t fMultRemovedKaonPlus; // enter comment here
 266     Int_t fMultRemovedKaonMinus; // enter comment here
 267     Int_t fMultRemovedK0s; // enter comment here
 268     Int_t fMultRemovedK0L; // enter comment here
 269
 270     Int_t fMultRemovedLambdas; // enter comment here
 271     Int_t fMultRemovedElectrons; // enter comment here
 272     Int_t fMultRemovedPositrons; // enter comment here
 273     Int_t fMultRemovedMuMinus; // enter comment here
 274     Int_t fMultRemovedMuPlus; // enter comment here
 275
 276     Int_t fMultRemovedGammas; // enter comment here
 277     Int_t fMultRemovedNeutrons; // enter comment here
 278     Int_t fMultRemovedAntiNeutrons; // enter comment here
 279
 280     Int_t fTrackMultInAcc; // enter comment here
 281
 282
 283     TH2F *fHistDxDzNonRemovedCharged; // enter comment here
 284     TH2F *fHistDxDzRemovedCharged; // enter comment here
 285     TH2F *fHistDxDzNonRemovedNeutral; // enter comment here
 286     TH2F *fHistDxDzRemovedNeutral; // enter comment here
 287
 288     TH1F *fHistPiPlusMult; // enter comment here
 289     TH1F *fHistPiMinusMult; // enter comment here
 290     TH1F *fHistPiZeroMult; // enter comment here
 291
 292     TH1F *fHistPiPlusMultAcc; // enter comment here
 293     TH1F *fHistPiMinusMultAcc; // enter comment here
 294     TH1F *fHistPiZeroMultAcc; // enter comment here
 295
 296    // Int_t fPiPlusMult; // enter comment here
 297    // Int_t fPiMinusMult; // enter comment here
 298
 299     Int_t fPiZeroMult; // enter comment here
 300
 301     Int_t fPiPlusMultAcc; // enter comment here
 302     Int_t fPiMinusMultAcc; // enter comment here
 303     Int_t fPiZeroMultAcc; // enter comment here
 304
 305
 306     Int_t fNeutralRemoved; // number of neutral particles that where removed by track matching
 307     Int_t fChargedRemoved; // number of charged particles that where removed by track matching
 308     Int_t fChargedNotRemoved; // number of charged particles that were not removed
 309     Int_t fNeutralNotRemoved; // number of neutral particles that were not removed
 310     Int_t fGammaRemoved; // number of gammas removed
 311
 312     Int_t fSecondaryNotRemoved;
 313
 314     Double_t fEnergyNeutralRemoved; // energy of neutral particles that where removed by track matching
 315     Double_t fEnergyChargedRemoved; // energy of charged particles that where removed by track matching
 316     Double_t fEnergyChargedNotRemoved; // energy of charged particles that were not removed
 317     Double_t fEnergyNeutralNotRemoved; // energy of neutral particles that were not removed
 318     Double_t fEnergyGammaRemoved; // energy of gammas that were removed
 319
 320     Int_t fNClusters; // Number of clusters in event
 321
 322     Double_t fTotNeutralEtAfterMinEnergyCut; // enter comment here
 323
 324     Bool_t fCalcTrackMatchVsMult;
 325     TH1F *fHistGammasFound;
 326     TH1F *fHistGammasGenerated;
 327     TH2F *fHistGammasFoundCent;
 328     TH2F *fHistGammasFoundOutOfAccCent;
 329     TH2F *fHistGammasFoundAltCent;
 330     TH2F *fHistGammasFoundOutOfAccAltCent;
 331     TH2F *fHistGammasGeneratedCent;
 332     TH2F *fHistGammasFoundRecoEnergyCent;
 333     TH2F *fHistGammasFoundOutOfAccRecoEnergyCent;
 334     TH1F *fHistChargedTracksCut;
 335     TH1F *fHistChargedTracksAccepted;
 336     TH1F *fHistGammasCut;
 337     TH1F *fHistGammasAccepted;
 338     TH2F *fHistChargedTrackDepositsAcceptedVsPt;
 339     TH2F *fHistChargedTrackDepositsAllVsPt;
 340     TH2F *fHistChargedTrackDepositsAcceptedVsPtEffCorr;
 341     TH2F *fHistChargedTrackDepositsAllVsPtEffCorr;
 342     TH2F *fHistChargedTracksCutMult;
 343     TH2F *fHistChargedTracksAcceptedMult;
 344     TH2F *fHistChargedTracksAcceptedLowPtCentEffCorr;
 345     TH2F *fHistChargedTracksAcceptedLowPtCent;
 346     TH2F *fHistChargedTracksAcceptedLowPtCent500MeV;
 347     TH2F *fHistChargedTracksAcceptedLowPtCentNoAntiProtons;
 348     TH2F *fHistChargedTracksAcceptedLowPtCentAntiProtons;
 349     TH2F *fHistGammasCutMult;
 350     TH2F *fHistGammasAcceptedMult;
 351     TH1F *fHistBadTrackMatches;
 352     TH2F *fHistMatchedTracksEvspTBkgd;
 353     TH2F *fHistMatchedTracksEvspTSignal;
 354     TH2F *fHistMatchedTracksEvspTBkgdPeripheral;
 355     TH2F *fHistMatchedTracksEvspTSignalPeripheral;
 356     TH3F *fHistMatchedTracksEvspTBkgdvsCent;
 357     TH3F *fHistMatchedTracksEvspTSignalvsCent;
 358     TH3F *fHistMatchedTracksEvspTBkgdvsCentEffCorr;
 359     TH3F *fHistMatchedTracksEvspTSignalvsCentEffCorr;
 360     TH1F *fHistChargedTracksCutPeripheral;
 361     TH1F *fHistChargedTracksAcceptedPeripheral;
 362     TH1F *fHistGammasCutPeripheral;
 363     TH1F *fHistGammasAcceptedPeripheral;
 364     TH2F *fHistBadTrackMatchesdPhidEta;
 365     TH2F *fHistGoodTrackMatchesdPhidEta;
 366     TH1F *fHistHadronDepositsAll;
 367     TH1F *fHistHadronDepositsReco;
 368     TH2F *fHistHadronDepositsAllCent;
 369     TH2F *fHistHadronDepositsAllCent500MeV;
 370     TH2F *fHistHadronDepositsRecoCent;
 371     TH2F *fHistHadronDepositsAllvsECent;
 372     TH2F *fHistHadronDepositsRecovsECent;
 373     TH2F *fHistHadronsAllCent;
 374     TH3F *fHistMultChVsSignalVsMult;
 375     TH2F *fHistNeutralRemovedSecondaryEtVsCent;
 376     TH2F *fHistChargedRemovedSecondaryEtVsCent;
 377     TH2F *fHistNeutralNotRemovedSecondaryEtVsCent;
 378     TH2F *fHistChargedNotRemovedSecondaryEtVsCent;
 379     TH2F *fHistNeutralRemovedSecondaryNumVsNCluster;
 380     TH2F *fHistChargedRemovedSecondaryNumVsNCluster;
 381     TH2F *fHistNeutralNotRemovedSecondaryNumVsNCluster;
 382     TH2F *fHistChargedNotRemovedSecondaryNumVsNCluster;
 383     TH2F *fHistNeutralRemovedSecondaryNumVsCent;
 384     TH2F *fHistChargedRemovedSecondaryNumVsCent;
 385     TH2F *fHistNeutralNotRemovedSecondaryNumVsCent;
 386     TH2F *fHistChargedNotRemovedSecondaryNumVsCent;
 387     TH2F *fHistNeutronsEtVsCent;
 388     TH2F *fHistNeutronsNumVsCent;
 389     TH2F *fHistNotNeutronsNumVsCent;
 390     TH2F *fHistPiKPDepositedVsNch;
 391     TH2F *fHistPiKPNotTrackMatchedDepositedVsNch;
 392
 393     TH2F *fHistNeutronsDepositedVsNch;
 394     TH2F *fHistAntiNeutronsDepositedVsNch;
 395     TH2F *fHistProtonsDepositedVsNch;
 396     TH2F *fHistAntiProtonsDepositedVsNch;
 397     TH2F *fHistProtonsNotTrackMatchedDepositedVsNch;
 398     TH2F *fHistAntiProtonsNotTrackMatchedDepositedVsNch;//
 399     TH2F *fHistNeutronsDepositedVsNcl;
 400     TH2F *fHistAntiNeutronsDepositedVsNcl;
 401     TH2F *fHistProtonsDepositedVsNcl;
 402     TH2F *fHistAntiProtonsDepositedVsNcl;
 403     TH2F *fHistProtonsNotTrackMatchedDepositedVsNcl;
 404     TH2F *fHistAntiProtonsNotTrackMatchedDepositedVsNcl;
 405     TH2F *fHistSecondariesVsNch;
 406     TH2F *fHistSecondariesVsNcl;
 407     TH2F *fHistSecondariesEffCorrVsNch;
 408     TH2F *fHistSecondariesEffCorrVsNcl;
 409     TH2F *fHistSecondariesOutOfAccEffCorrVsNch;
 410     TH2F *fHistSecondariesDetectorCoverEffCorrVsNch;
 411
 412
 413     TH2F *fHistNeutronsDepositedVsNchNoEffCorr;//filled
 414     TH2F *fHistAntiNeutronsDepositedVsNchNoEffCorr;//filled
 415     TH2F *fHistProtonsDepositedVsNchNoEffCorr;//filled
 416     TH2F *fHistAntiProtonsDepositedVsNchNoEffCorr;//filled
 417     TH2F *fHistProtonsNotTrackMatchedDepositedVsNchNoEffCorr;//filled
 418     TH2F *fHistAntiProtonsNotTrackMatchedDepositedVsNchNoEffCorr;//filled
 419     TH2F *fHistNeutronsDepositedVsNclNoEffCorr;//filled
 420     TH2F *fHistAntiNeutronsDepositedVsNclNoEffCorr;//filled
 421     TH2F *fHistProtonsDepositedVsNclNoEffCorr;//filled
 422     TH2F *fHistAntiProtonsDepositedVsNclNoEffCorr;//filled
 423     TH2F *fHistProtonsNotTrackMatchedDepositedVsNclNoEffCorr;//filled
 424     TH2F *fHistAntiProtonsNotTrackMatchedDepositedVsNclNoEffCorr;//filled
 425
 426     TH3F *fHistCentVsNchVsNcl;
 427     TH3F *fHistSecondaryPositionInDetector;
 428     //TH2F *fHistSecondaryPositionInDetector2D;
 429     TH2F *fClusterPositionWeird;
 430     //TH3F *fHistSecondaryPositionInDetectorMultiple;
 431     TH1F *fSecondaryClusterEnergy; // Distribution of cluster energies
 432
 433     TH2F *fHistGammaCrossCheck;
 434     TH2F *fHistGammaCrossCheckAlt;//includes gammas that did not hit as gammas
 435     TH2F *fHistGammaEnergyCrossCheck;
 436     TH2F *fHistGammaEnergyCrossCheckAlt;//includes gammas that did not hit as gammas
 437     TH2F *fHistNeutronCrossCheck;
 438     TH2F *fHistSecondaryCrossCheck;
 439     TH2F *fHistHadronCrossCheck;
 440     TH2F *fHistKaonCrossCheck;
 441
 442     TH1F *fHistAllEnergy;//energy of all clusters passing cuts vs centrality
 443     TH1F *fHistSignalEnergy;//signal of signal clusters passing cuts vs centrality
 444     TH1F *fHistNeutronEnergy;//signal of neutron clusters passing cuts vs centrality
 445     TH1F *fHistKaonEnergy;//signal of kaon clusters passing cuts vs centrality
 446     TH1F *fHistHadronEnergy;//signal of hadron clusters passing cuts vs centrality
 447     TH1F *fHistSecondaryEnergy;//signal of secondary clusters passing cuts vs centrality
 448     TH1F *fHistSecondaryChargedEnergy;//signal of secondary clusters passing cuts vs centrality
 449     TH1F *fHistSecondaryNeutronEnergy;//signal of secondary clusters passing cuts vs centrality
 450     TH1F *fHistSecondaryGammaEnergy;//signal of secondary clusters passing cuts vs centrality
 451     TH1F *fHistSecondaryElectronEnergy;//signal of secondary clusters passing cuts vs centrality
 452     TH1F *fHistSecondaryOtherEnergy;//signal of secondary clusters passing cuts vs centrality
 453     TH1F *fHistSimulatedGammaEnergy;//signal of signal clusters passing cuts vs centrality
 454     TH1F *fHistReconstructedGammaEnergy;//signal of signal clusters passing cuts vs centrality
 455     TH1F *fHistSimulatedGammaEnergyAboveThreshold;//signal of signal clusters passing cuts vs centrality
 456     TH1F *fHistReconstructedSignalEnergy;//signal of signal clusters passing cuts vs centrality
 457
 458     TH2F *fHistFracSignalVsNClusters;//Fraction of signal vs number of clusters
 459     TH2F *fHistFracHadronsVsNClusters;//Fraction of hadrons vs number of clusters
 460     TH2F *fHistFracNeutronsVsNClusters;//Fraction of signal vs number of clusters
 461     TH2F *fHistFracKaonsVsNClusters;//Fraction of signal vs number of clusters
 462     TH2F *fHistFracSecondariesVsNClusters;//Fraction of signal vs number of clusters
 463     TH2F *fHistFracSignalVsNMultiplicity;//Fraction of signal vs charged track multiplicity
 464     TH2F *fHistFracHadronsVsNMultiplicity;//Fraction of hadrons vs charged track multiplicity
 465     TH2F *fHistFracNeutronsVsNMultiplicity;//Fraction of signal vs charged track multiplicity
 466     TH2F *fHistFracKaonsVsNMultiplicity;//Fraction of signal vs charged track multiplicity
 467     TH2F *fHistFracSecondariesVsNMultiplicity;//Fraction of signal vs charged track multiplicity
 468     TH2F *fHistFracSignalVsNMatchedTracks;//Fraction of signal vs number of matched tracks
 469     TH2F *fHistFracHadronsVsNMatchedTracks;//Fraction of hadrons vs number of matched tracks
 470     TH2F *fHistFracNeutronsVsNMatchedTracks;//Fraction of signal vs number of matched tracks
 471     TH2F *fHistFracKaonsVsNMatchedTracks;//Fraction of signal vs number of matched tracks
 472     TH2F *fHistFracSecondariesVsNMatchedTracks;//Fraction of signal vs number of matched tracks
 473     TH2F *fHistFracSignalVsNTotalTracks;//Fraction of signal vs number of total tracks
 474     TH2F *fHistFracHadronsVsNTotalTracks;//Fraction of hadrons vs number of total tracks
 475     TH2F *fHistFracNeutronsVsNTotalTracks;//Fraction of signal vs number of total tracks
 476     TH2F *fHistFracKaonsVsNTotalTracks;//Fraction of signal vs number of total tracks
 477     TH2F *fHistFracSecondariesVsNTotalTracks;//Fraction of signal vs number of total tracks
 478
 479
 480 private:
 481
 482     //Declare it private to avoid compilation warning
 483     AliAnalysisEtMonteCarlo & operator = (const AliAnalysisEtMonteCarlo & g) ;//cpy assignment
 484     AliAnalysisEtMonteCarlo(const AliAnalysisEtMonteCarlo & g) ; // cpy ctor
 485     ClassDef(AliAnalysisEtMonteCarlo, 2);
 486 };
 487
 488 #endif // ALIANALYSISETMONTECARLO_H