PWGGA/CaloTrackCorrelations/AliAnaParticleHadronCorrelation.h

   1 #ifndef ALIANAPARTICLEHADRONCORRELATION_H
   2 #define ALIANAPARTICLEHADRONCORRELATION_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice     */
   5
   6 //_________________________________________________________________________
   7 // Class that contains the algorithm for the analysis of
   8 // particle - hadron correlations
   9 // Particle (for example direct gamma) must be found in a previous analysis
  10 //
  11 //-- Author: Gustavo Conesa (LNF-INFN) (LPSC-IN2P3-CNRS)
  12 //           Yaxian Mao (LPSC-IN2P3-CNRS) and (CNWU) first usable implementation.
  13 //           Xiangrong Zhu (CNWU), implementtion of own mixing.
  14 //
  15
  16 // --- Analysis system ---
  17
  18 #include "AliAnaCaloTrackCorrBaseClass.h"
  19 class AliAODPWG4ParticleCorrelation ;
  20
  21 class AliAnaParticleHadronCorrelation : public AliAnaCaloTrackCorrBaseClass {
  22
  23  public:
  24
  25   AliAnaParticleHadronCorrelation() ;          // default ctor
  26   virtual ~AliAnaParticleHadronCorrelation() ; // virtual dtor
  27
  28   // General methods
  29
  30   TObjString * GetAnalysisCuts();
  31
  32   TList      * GetCreateOutputObjects();
  33
  34   void         Init();
  35
  36   void         InitParameters();
  37
  38   void         FillEventMixPool() ;
  39
  40   void         MakeAnalysisFillHistograms() ;
  41
  42   void         Print(const Option_t * opt) const;
  43
  44   // Main analysis methods
  45
  46   Bool_t       FindLeadingOppositeHadronInWindow(AliAODPWG4ParticleCorrelation * particle);
  47
  48   Bool_t       GetDecayPhotonMomentum   (Int_t indexPhoton1, Int_t indexPhoton2, Int_t idetector);
  49
  50   void         MakeChargedCorrelation   (AliAODPWG4ParticleCorrelation * particle) ;
  51
  52   void         MakeNeutralCorrelation   (AliAODPWG4ParticleCorrelation * particle) ;
  53
  54   void         MakeMCChargedCorrelation (Int_t triggerMCLable, Int_t histoIndex, Bool_t lostDecayPair) ;
  55
  56   void         MakeChargedMixCorrelation(AliAODPWG4ParticleCorrelation * particle) ;
  57
  58   // Filling histogram methods
  59
  60   void         FillChargedAngularCorrelationHistograms  (Float_t ptAssoc,  Float_t ptTrig,      Int_t   assocBin,
  61                                                          Float_t phiAssoc, Float_t phiTrig,     Float_t deltaPhi,
  62                                                          Float_t etaAssoc, Float_t etaTrig,
  63                                                          Int_t   decayTag, Float_t hmpidSignal, Int_t outTOF,
  64                                                          Int_t   cenbin,   Int_t   mcTag);
  65
  66   void         FillChargedEventMixPool();
  67
  68   Bool_t       FillChargedMCCorrelationHistograms       (Float_t mcAssocPt, Float_t mcAssocPhi, Float_t mcAssocEta,
  69                                                          Float_t mcTrigPt,  Float_t mcTrigPhi,  Float_t mcTrigEta,
  70                                                          Int_t histoIndex,  Bool_t  lostDecayPair);
  71
  72   void         FillChargedMomentumImbalanceHistograms   (Float_t ptTrig,   Float_t ptAssoc,
  73                                                          Float_t deltaPhi, Int_t cenbin, Int_t charge,
  74                                                          Int_t   assocBin, Int_t decayTag,
  75                                                          Int_t   outTOF,   Int_t mcTag );
  76
  77   void         FillChargedUnderlyingEventHistograms     (Float_t ptTrig,   Float_t ptAssoc,
  78                                                          Float_t deltaPhi, Int_t cenbin, Int_t outTOF);
  79
  80   void         FillChargedUnderlyingEventSidesHistograms(Float_t ptTrig,   Float_t ptAssoc,
  81                                                          Float_t deltaPhi);
  82
  83   void         FillDecayPhotonCorrelationHistograms     (Float_t ptAssoc,  Float_t phiAssoc, Bool_t bChargedOrNeutral);
  84
  85   void         FillNeutralEventMixPool();
  86
  87
  88   void         FillNeutralUnderlyingEventSidesHistograms(Float_t ptTrig,   Float_t ptAssoc,
  89                                                          Float_t zT,       Float_t hbpZT,
  90                                                          Float_t deltaPhi);
  91
  92   Int_t        GetMCTagHistogramIndex(Int_t tag);
  93   static const Int_t fgkNmcTypes = 10;
  94
  95   Bool_t       IsTriggerTheEventLeadingParticle();
  96
  97   // Parameter setter and getter
  98
  99   Float_t      GetMinimumTriggerPt()       const { return GetMinPt()             ; }
 100   Float_t      GetMaximumTriggerPt()       const { return GetMaxPt()             ; }
 101   void         SetTriggerPtRange(Float_t min, Float_t max)
 102                { SetMinPt(min), SetMaxPt(max)                                    ; }
 103
 104
 105   Float_t      GetMaximumAssociatedPt()    const { return fMaxAssocPt            ; }
 106   Float_t      GetMinimumAssociatedPt()    const { return fMinAssocPt            ; }
 107   void         SetAssociatedPtRange(Float_t min, Float_t max)
 108                { fMaxAssocPt   = max ;           fMinAssocPt  = min              ; }
 109
 110   Double_t     GetDeltaPhiMaxCut()         const { return fDeltaPhiMaxCut        ; }
 111   Double_t     GetDeltaPhiMinCut()         const { return fDeltaPhiMinCut        ; }
 112   void         SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
 113                { fDeltaPhiMaxCut   = phimax ;    fDeltaPhiMinCut   = phimin      ; }
 114
 115   // Leading Hadron
 116   Double_t     GetLeadHadronPhiMaxCut()    const { return fMaxLeadHadPhi         ; }
 117   Double_t     GetLeadHadronPhiMinCut()    const { return fMinLeadHadPhi         ; }
 118   void         SetLeadHadronPhiCut(Float_t min, Float_t max)
 119                { fMaxLeadHadPhi = max ;          fMinLeadHadPhi  = min           ; }
 120
 121   Double_t     GetLeadHadronPtMinCut()     const { return fMinLeadHadPt          ; }
 122   Double_t     GetLeadHadronPtMaxCut()     const { return fMaxLeadHadPt          ; }
 123   void         SetLeadHadronPtCut(Float_t min, Float_t max)
 124                { fMaxLeadHadPt  = max ;           fMinLeadHadPt  = min           ; }
 125
 126   Bool_t       IsLeadHadronCutOn()        const { return fSelectLeadingHadronAngle   ; }
 127   void         SwitchOnLeadHadronSelection()    { fSelectLeadingHadronAngle = kTRUE  ; }
 128   void         SwitchOffLeadHadronSelection()   { fSelectLeadingHadronAngle = kFALSE ; }
 129
 130   void         SwitchOnFillLeadHadronHistograms() { fFillLeadHadOppositeHisto = kTRUE  ; }
 131   void         SwitchOffFillLeadHadronHistograms(){ fFillLeadHadOppositeHisto = kFALSE ; }
 132
 133   // UE
 134
 135   Double_t     GetUeDeltaPhiMaxCut()       const { return fUeDeltaPhiMaxCut      ; }
 136   Double_t     GetUeDeltaPhiMinCut()       const { return fUeDeltaPhiMinCut      ; }
 137
 138   void         SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax)
 139                   { fUeDeltaPhiMaxCut = uephimax ;  fUeDeltaPhiMinCut = uephimin ; }
 140
 141   Bool_t       IsSeveralUEOn()             const { return fMakeSeveralUE         ; }
 142   void         SwitchOnSeveralUECalculation()    { fMakeSeveralUE      = kTRUE   ; }
 143   void         SwitchOffSeveralUECalculation()   { fMakeSeveralUE      = kFALSE  ; }
 144
 145   // Do trigger-neutral correlation
 146   Bool_t       DoNeutralCorr()             const { return fNeutralCorr           ; }
 147   void         SwitchOnNeutralCorr()             { fNeutralCorr      = kTRUE     ; }
 148   void         SwitchOffNeutralCorr()            { fNeutralCorr      = kFALSE    ; }
 149
 150   // Taking the absolute leading as the trigger or not
 151   Bool_t       DoAbsoluteLeading()         const { return fMakeAbsoluteLeading   ; }
 152   void         SwitchOnAbsoluteLeading()         { fMakeAbsoluteLeading = kTRUE  ; }
 153   void         SwitchOffAbsoluteLeading()        { fMakeAbsoluteLeading = kFALSE ; }
 154
 155   // Taking the near side leading as the trigger or not
 156   Bool_t       DoNearSideLeading()         const { return fMakeNearSideLeading   ; }
 157   void         SwitchOnNearSideLeading()         { fMakeNearSideLeading = kTRUE  ; }
 158   void         SwitchOffNearSideLeading()        { fMakeNearSideLeading = kFALSE ; }
 159
 160   // Do decay-hadron correlation if it is pi0 trigger
 161   Bool_t       IsPi0Trigger()              const { return fPi0Trigger            ; }
 162   void         SwitchOnPi0TriggerDecayCorr()     { fPi0Trigger          = kTRUE  ; }
 163   void         SwitchOffPi0TriggerDecayCorr()    { fPi0Trigger          = kFALSE ; }
 164
 165   Bool_t       IsDecayTrigger()            const { return fDecayTrigger          ; }
 166   void         SwitchOnDecayTriggerDecayCorr()   { fDecayTrigger        = kTRUE  ; }
 167   void         SwitchOffDecayTriggerDecayCorr()  { fDecayTrigger        = kFALSE ; }
 168   void         SetNDecayBits(Int_t n)            { fNDecayBits = n               ; }
 169   void         SetDecayBits(Int_t i, UInt_t bit) { if(i < 4) fDecayBits[i] = bit ; }
 170
 171   Bool_t       IsHMPIDCorrelation()        const { return fHMPIDCorrelation      ; }
 172   void         SwitchOnHMPIDCorrelation()        { fHMPIDCorrelation    = kTRUE  ; }
 173   void         SwitchOffHMPIDCorrelation()       { fHMPIDCorrelation    = kFALSE ; }
 174
 175   void         SwitchOnFillBradHistograms()      { fFillBradHisto       = kTRUE  ; }
 176   void         SwitchOffFillBradHistograms()     { fFillBradHisto       = kFALSE ; }
 177
 178   Bool_t       OnlyIsolated()              const { return fSelectIsolated        ; }
 179   void         SelectIsolated(Bool_t s)          { fSelectIsolated   = s         ; }
 180
 181   void         SetPi0AODBranchName(TString n)    { fPi0AODBranchName = n         ; }
 182
 183   void         SetNAssocPtBins(Int_t n) ;
 184   void         SetAssocPtBinLimit(Int_t ibin, Float_t pt) ;
 185
 186   Bool_t       IsMixStoredInReaderOn()     const { return fUseMixStoredInReader  ; }
 187   void         SwitchOnUseMixStoredInReader()    { fUseMixStoredInReader = kTRUE ; }
 188   void         SwitchOffUseMixStoredInReader()   { fUseMixStoredInReader = kFALSE; }
 189
 190   void         SwitchOnFillNeutralInMixedEvent() { fFillNeutralEventMixPool = kTRUE  ; }
 191   void         SwitchOffFillNeutralInMixedEvent(){ fFillNeutralEventMixPool = kFALSE ; }
 192
 193   void         SetM02Cut(Float_t min=0, Float_t max=10)  { fM02MinCut   = min ; fM02MaxCut  = max ; }
 194
 195   void         SwitchOnCorrelationVzBin()        { fCorrelVzBin          = kTRUE  ; }
 196   void         SwitchOffCorrelationVzBin()       { fCorrelVzBin          = kFALSE ; }
 197
 198   void         SwitchOnFillTriggerAODWithReferences()   { fFillAODWithReferences = kTRUE  ; }
 199   void         SwitchOffFillTriggerAODWithReferences()  { fFillAODWithReferences = kFALSE ; }
 200
 201   void         SwitchOnCheckNeutralClustersForLeading() { fCheckLeadingWithNeutralClusters = kTRUE  ; }
 202   void         SwitchOffCheckNeutralClustersForLeading(){ fCheckLeadingWithNeutralClusters = kFALSE ; }
 203
 204   void         SwitchOnFillEtaGapHistograms()    { fFillEtaGapsHisto    = kTRUE  ; }
 205   void         SwitchOffFillEtaGapHistograms()   { fFillEtaGapsHisto    = kFALSE ; }
 206
 207   void         SwitchOnFillPtImbalancePerPtABinHistograms()  { fFillMomImbalancePtAssocBinsHisto = kTRUE  ; }
 208   void         SwitchOffFillPtImbalancePerPtABinHistograms() { fFillMomImbalancePtAssocBinsHisto = kFALSE ; }
 209
 210   void         SetMCGenType(Int_t min = 0, Int_t max = 6) { if(min >= 0 && min < fgkNmcTypes) fMCGenTypeMin = min ;
 211                                                             if(max >= 0 && max < fgkNmcTypes) fMCGenTypeMax = max ; }
 212
 213  private:
 214
 215   Bool_t       fFillAODWithReferences;         // Add to the trigger particle AOD the reference to the tracks or neutrals in correlation.
 216   Bool_t       fCheckLeadingWithNeutralClusters;// Compare the trigger candidate to Leading pT with the clusters pT, by default only charged
 217   Float_t      fMaxAssocPt ;                   // Maximum associated hadron pt
 218   Float_t      fMinAssocPt ;                   // Minimum associated hadron pt
 219   Double_t     fDeltaPhiMaxCut ;               // Minimum Delta Phi Gamma-Hadron
 220   Double_t     fDeltaPhiMinCut ;               // Maximum Delta Phi Gamma-Hadron
 221   Bool_t       fSelectIsolated ;               // Select only trigger particles isolated
 222   Bool_t       fMakeSeveralUE ;                // Do analysis for several underlying events contribution
 223   Double_t     fUeDeltaPhiMaxCut ;             // Minimum Delta Phi Gamma-Underlying Hadron
 224   Double_t     fUeDeltaPhiMinCut ;             // Maximum Delta Phi Gamma-Underlying Hadron
 225   TString      fPi0AODBranchName;              // Name of AOD branch with pi0, not trigger
 226   Bool_t       fNeutralCorr ;                  // switch the analysis with neutral particles
 227   Bool_t       fPi0Trigger ;                   // switch the analysis with decay photon from pi0 trigger
 228   Bool_t       fDecayTrigger ;                 // switch the analysis with decay photon from photon trigger
 229   Int_t        fNDecayBits ;                   // in case of study of decay triggers, select the decay bit
 230   UInt_t       fDecayBits[4] ;                 // in case of study of decay triggers, select the decay bit
 231   Bool_t       fMakeAbsoluteLeading ;          // requesting absolute leading triggers
 232   Bool_t       fMakeNearSideLeading ;          // requesting near side leading (+-90º from trigger particle) triggers
 233   Int_t        fLeadingTriggerIndex ;          // Store here per event the trigger index, to avoid too many loops
 234   Bool_t       fHMPIDCorrelation    ;          // Correlate with particles on HMPID or its acceptance
 235   Bool_t       fFillBradHisto ;                // DPhi histograms calculated differently
 236   Int_t        fNAssocPtBins ;                 // Number of associated pT bins under study
 237   Float_t      fAssocPtBinLimit[20] ;          // Associated pT under study
 238   Bool_t       fCorrelVzBin ;                  // Fill one histogram per vz bin
 239
 240   TList **     fListMixTrackEvents ;           //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for tracks in stored events for mixing
 241   TList **     fListMixCaloEvents ;            //![GetNCentrBin()*GetNZvertBin()*GetNRPBin()] Containers for calo clusters in stored events for mixing
 242
 243   Bool_t       fUseMixStoredInReader;          // Signal if in the current event the pool was filled
 244   Bool_t       fFillNeutralEventMixPool;       // Add clusters to pool if requested
 245
 246   Float_t      fM02MaxCut   ;                  // Study photon clusters with l0 smaller than cut
 247   Float_t      fM02MinCut   ;                  // Study photon clusters with l0 larger than cut
 248
 249   Bool_t       fSelectLeadingHadronAngle;      // Select events with leading particle within a range
 250   Bool_t       fFillLeadHadOppositeHisto;      // Fill histograms for leading hadrons in opposite side of trigger
 251
 252   Float_t      fMinLeadHadPhi;                 // Minimum angle between the trigger and leading hadron
 253   Float_t      fMaxLeadHadPhi;                 // Maximum ange between the trigger and leading hadron
 254   Float_t      fMinLeadHadPt;                  // Minimum pT of leading hadron
 255   Float_t      fMaxLeadHadPt;                  // Maximum pT of leading hadron
 256
 257   Bool_t       fFillEtaGapsHisto;              // Fill azimuthal correlation histograms in 2 eta gaps, |eta|>0.8 and |eta|<0.01
 258   Bool_t       fFillMomImbalancePtAssocBinsHisto; // momentum imbalance histograms in bins of pT associated
 259
 260   Int_t        fMCGenTypeMin;                  // Of the fgkNmcTypes possible types, select those between fMCGenTypeMin and fMCGenTypeMax
 261   Int_t        fMCGenTypeMax;                  // Of the fgkNmcTypes possible types, select those between fMCGenTypeMin and fMCGenTypeMax
 262
 263   TVector3       fTrackVector;                 //! track momentum vector
 264   TLorentzVector fMomentum;                    //! trigger momentum
 265   TLorentzVector fDecayMom1;                   //! decay particle momentum
 266   TLorentzVector fDecayMom2;                   //! decay particle momentum
 267
 268   //Histograms
 269
 270   //trigger particles
 271   TH1F *       fhPtTriggerInput;               //! pT distribution of trigger particles before selection
 272   TH1F *       fhPtTriggerSSCut;               //! pT distribution of trigger particles after shower shape selection
 273   TH1F *       fhPtTriggerIsoCut;              //! pT distribution of trigger particles after isolation cut selection
 274   TH1F *       fhPtTriggerFidCut;              //! pT distribution of trigger particles after fiducial selection
 275   TH1F *       fhPtTrigger;                    //! pT distribution of trigger particles
 276   TH1F *       fhPtTriggerVtxBC0;              //! pT distribution of trigger particles
 277   TH1F *       fhPtTriggerPileUp[7];           //! pT distribution of trigger particles
 278   TH2F *       fhPtTriggerVzBin;               //! pT distribution of trigger particles vs vz bin
 279   TH2F *       fhPtTriggerBin;                 //! pT distribution of trigger particles, vs mixing bin
 280   TH2F *       fhPhiTrigger;                   //! phi distribution vs pT of trigger particles
 281   TH2F *       fhEtaTrigger;                   //! eta distribution vs pT of trigger particles
 282
 283   TH1F *       fhPtTriggerMC[fgkNmcTypes];     //! pT distribution of trigger particles, check the origin of the cluster : "Photon","Pi0","Pi0Decay","EtaDecay","OtherDecay","Electron","Hadron"
 284
 285   TH1F *       fhPtDecayTrigger[4];            //! pT distribution of trigger particles, tagged as decay
 286   TH1F *       fhPtDecayTriggerMC[4][fgkNmcTypes];//! pT distribution of trigger particles, tagged as decay, check the origin of the cluster
 287
 288   TH2F *       fhPtTriggerCentrality;          //! pT distribution of trigger particles vs centrality
 289   TH2F *       fhPtTriggerEventPlane;          //! pT distribution of trigger particles vs centrality
 290   TH2F *       fhTriggerEventPlaneCentrality;  //! event plane vs centrality for trigger particles
 291
 292   TH1F *       fhPtTriggerMixed;               //! pT distribution of trigger particles, used in mixing
 293   TH2F *       fhPtTriggerMixedVzBin;          //! pT distribution of trigger particles, used in mixing, vs vz bin
 294   TH2F *       fhPtTriggerMixedBin;            //! pT distribution of trigger particles vs mixing bin
 295   TH2F *       fhPhiTriggerMixed;              //! phi distribution vs pT of trigger particles, used in mixing
 296   TH2F *       fhEtaTriggerMixed;              //! eta distribution vs pT of trigger particles, used in mixing
 297
 298   // Leading hadron in the opposite side of the trigger
 299   TH2F *       fhPtLeadingOppositeHadron;        //! pT trigger : pT distribution of leading hadron oposite to trigger
 300   TH2F *       fhPtDiffPhiLeadingOppositeHadron; //! pT trigger : difference phi distribution of leading hadron oposite and trigger
 301   TH2F *       fhPtDiffEtaLeadingOppositeHadron; //! pT trigger: difference eta distribution of leading hadron oposite and trigger
 302   TH1F *       fhPtNoLeadingOppositeHadron;      //! pT trigger for events without opposite hadrons
 303   TH2F *       fhEtaPhiNoLeadingOppositeHadron;  //! location of trigger when no hadron is found on the opposite side
 304
 305   //trigger-charged histograms
 306   TH2F *       fhDeltaPhiDeltaEtaCharged ;     //! differences of eta and phi between trigger and charged hadrons
 307   TH2F *       fhPhiCharged  ;                 //! Phi distribution of charged particles
 308   TH2F *       fhEtaCharged  ;                 //! Eta distribution of charged particles
 309   TH2F *       fhDeltaPhiCharged  ;            //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 310   TH2F *       fhDeltaEtaCharged  ;            //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT
 311   TH2F *       fhDeltaPhiChargedPt  ;          //! Difference of charged particle phi and trigger particle  phi as function of charged particle pT
 312   TH2F *       fhDeltaPhiUeChargedPt ;         //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT
 313   TH1F *       fhUePart;                       //! UE particles distribution vs pt trig
 314   TH2F *       fhXECharged  ;                  //! Trigger particle -charged hadron momentum imbalance histogram
 315   TH2F *       fhXECharged_Cone2  ;            //! Trigger particle -charged hadron momentum imbalance histogram in cone2 (5pi/6-7pi/6)
 316   TH2F *       fhXEUeCharged  ;                //! Trigger particle -underlying charged hadron momentum imbalance histogram
 317   TH2F *       fhXEPosCharged  ;               //! Trigger particle -positive charged hadron momentum imbalance histogram
 318   TH2F *       fhXENegCharged  ;               //! Trigger particle -negative charged hadron momentum imbalance histogram
 319   TH2F *       fhPtHbpXECharged  ;             //! Trigger particle -charged hadron momentum HBP histogram
 320   TH2F *       fhPtHbpXECharged_Cone2  ;       //! Trigger particle -charged hadron momentum HBP histogram in cone2 (5pi/6-7pi/6)
 321   TH2F *       fhPtHbpXEUeCharged  ;           //! Trigger particle -underlying charged hadron momentum HBP histogram
 322   TH2F *       fhZTCharged  ;                  //! Trigger particle -charged hadron momentum imbalance histogram
 323   TH2F *       fhZTUeCharged  ;                //! Trigger particle -underlying charged hadron momentum imbalance histogram
 324   TH2F *       fhZTPosCharged  ;               //! Trigger particle -positive charged hadron momentum imbalance histogram
 325   TH2F *       fhZTNegCharged  ;               //! Trigger particle -negative charged hadron momentum imbalance histogram
 326   TH2F *       fhPtHbpZTCharged  ;             //! Trigger particle -charged hadron momentum HBP histogram
 327   TH2F *       fhPtHbpZTUeCharged  ;           //! Trigger particle -underlying charged hadron momentum HBP histogram
 328
 329   TH2F *       fhXEChargedMC[fgkNmcTypes]  ;   //! Trigger particle -charged hadron momentum imbalance histogram, check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0), hadron, rest of photons (prompt, FSR, ISR)
 330   TH2F *       fhDeltaPhiChargedMC[fgkNmcTypes];//! Trigger particle -charged hadron delta phi histogram, check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0), hadron, rest of photons (prompt, FSR, ISR)
 331
 332   TH2F *       fhDeltaPhiDeltaEtaChargedPtA3GeV;//! differences of eta and phi between trigger and charged hadrons, pTa > 3 GeV
 333   TH2F *       fhDeltaPhiChargedPtA3GeV  ;      //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, pTa > 3 GeV
 334   TH2F *       fhDeltaEtaChargedPtA3GeV  ;      //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT, pTa > 3 GeV
 335
 336   // Events tagged as pileup by SDD,EMCal, or combination
 337   TH2F *       fhDeltaPhiChargedPileUp[7]  ;    //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 338   TH2F *       fhDeltaEtaChargedPileUp[7]  ;    //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT
 339   TH2F *       fhDeltaPhiChargedPtA3GeVPileUp[7] ; //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, pTa > 3 GeV
 340   TH2F *       fhDeltaEtaChargedPtA3GeVPileUp[7] ; //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT, pTa > 3 GeV
 341   TH2F *       fhXEChargedPileUp[7]  ;          //! Trigger particle -charged hadron momentum imbalance histogram
 342   TH2F *       fhXEUeChargedPileUp[7]  ;        //! Trigger particle -charged hadron momentum imbalance histogram
 343   TH2F *       fhZTChargedPileUp[7]  ;          //! Trigger particle -charged hadron momentum imbalance histogram
 344   TH2F *       fhZTUeChargedPileUp[7]  ;        //! Trigger particle -charged hadron momentum imbalance histogram
 345   TH2F *       fhPtTrigChargedPileUp[7] ;       //! trigger and correlated particl pt, to be used for mean value for kt
 346
 347   TH2F *       fhDeltaPhiChargedOtherBC  ;       //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 348   TH2F *       fhDeltaPhiChargedPtA3GeVOtherBC ; //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, pTa > 3 GeV
 349   TH2F *       fhXEChargedOtherBC  ;             //! Trigger particle -charged hadron momentum imbalance histogram
 350   TH2F *       fhXEUeChargedOtherBC  ;           //! Trigger particle -charged hadron momentum imbalance histogram
 351   TH2F *       fhZTChargedOtherBC  ;             //! Trigger particle -charged hadron momentum imbalance histogram
 352   TH2F *       fhZTUeChargedOtherBC  ;           //! Trigger particle -charged hadron momentum imbalance histogram
 353   TH2F *       fhPtTrigChargedOtherBC ;          //! trigger and correlated particl pt, to be used for mean value for kt
 354
 355   TH2F *       fhDeltaPhiChargedBC0  ;           //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 356   TH2F *       fhDeltaPhiChargedPtA3GeVBC0 ;     //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, pTa > 3 GeV
 357   TH2F *       fhXEChargedBC0  ;                 //! Trigger particle -charged hadron momentum imbalance histogram
 358   TH2F *       fhXEUeChargedBC0  ;               //! Trigger particle -charged hadron momentum imbalance histogram
 359   TH2F *       fhZTChargedBC0  ;                 //! Trigger particle -charged hadron momentum imbalance histogram
 360   TH2F *       fhZTUeChargedBC0  ;               //! Trigger particle -charged hadron momentum imbalance histogram
 361   TH2F *       fhPtTrigChargedBC0 ;              //! trigger and correlated particl pt, to be used for mean value for kt
 362
 363   TH2F *       fhDeltaPhiChargedVtxBC0  ;        //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 364   TH2F *       fhDeltaPhiChargedPtA3GeVVtxBC0 ;  //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, pTa > 3 GeV
 365   TH2F *       fhXEChargedVtxBC0  ;              //! Trigger particle -charged hadron momentum imbalance histogram
 366   TH2F *       fhXEUeChargedVtxBC0  ;            //! Trigger particle -charged hadron momentum imbalance histogram
 367   TH2F *       fhZTChargedVtxBC0  ;              //! Trigger particle -charged hadron momentum imbalance histogram
 368   TH2F *       fhZTUeChargedVtxBC0  ;            //! Trigger particle -charged hadron momentum imbalance histogram
 369   TH2F *       fhPtTrigChargedVtxBC0 ;           //! trigger and correlated particl pt, to be used for mean value for kt
 370
 371   //if several UE calculation is on, most useful for jet-jet events contribution
 372   TH2F *       fhDeltaPhiUeLeftCharged  ;      //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT
 373   TH2F *       fhDeltaPhiUeLeftUpCharged;      //! Difference of charged particle from underlying events phi and trigger particle  phi
 374   TH2F *       fhDeltaPhiUeRightUpCharged;     //! Difference of charged particle from underlying events phi and trigger particle  phi
 375   TH2F *       fhDeltaPhiUeLeftDownCharged;    //! Difference of charged particle from underlying events phi and trigger particle  phi
 376   TH2F *       fhDeltaPhiUeRightDownCharged;   //! Difference of charged particle from underlying events phi and trigger particle  phi
 377   TH2F *       fhXEUeLeftCharged  ;            //! Trigger particle -underlying charged hadron momentum imbalance histogram
 378   TH2F *       fhXEUeLeftUpCharged  ;          //! Trigger particle -underlying charged hadron momentum imbalance histogram
 379   TH2F *       fhXEUeRightUpCharged ;          //! Trigger particle -underlying charged hadron momentum imbalance histogram
 380   TH2F *       fhXEUeLeftDownCharged  ;        //! Trigger particle -underlying charged hadron momentum imbalance histogram
 381   TH2F *       fhXEUeRightDownCharged ;        //! Trigger particle -underlying charged hadron momentum imbalance histogram
 382   TH2F *       fhPtHbpXEUeLeftCharged  ;       //! Trigger particle -underlying charged hadron momentum HBP histogram
 383   TH2F *       fhZTUeLeftCharged  ;            //! Trigger particle -underlying charged hadron momentum imbalance histogram
 384   TH2F *       fhPtHbpZTUeLeftCharged  ;       //! Trigger particle -underlying charged hadron momentum HBP histogram
 385
 386   //for pout and kt extraction
 387   TH2F *       fhPtTrigPout  ;                 //! Pout =associated pt*sin(delta phi) distribution vs trigger pt
 388   TH2F *       fhPtTrigCharged ;               //! trigger and correlated particl pt, to be used for mean value for kt
 389
 390   //if different multiplicity analysis asked
 391   TH2F **      fhDeltaPhiChargedMult ;         //![GetNCentrBin()] differences of phi between trigger and charged hadrons: multiplicity bin
 392   TH2F **      fhDeltaEtaChargedMult ;         //![GetNCentrBin()] differences of eta between trigger and charged hadrons: multiplicity bin
 393   TH2F **      fhXEMult  ;                     //![GetNCentrBin()] Trigger particle -charged hadron momentum imbalance histogram: multiplicity bin
 394   TH2F **      fhXEUeMult  ;                   //![GetNCentrBin()] Trigger particle -UE charged hadron momentum imbalance histogram: multiplicity bin
 395   TH2F **      fhZTMult  ;                     //![GetNCentrBin()] Trigger particle -charged hadron momentum imbalance histogram: multiplicity bin
 396   TH2F **      fhZTUeMult  ;                   //![GetNCentrBin()] Trigger particle -UE charged hadron momentum imbalance histogram: multiplicity bin
 397
 398   TH2F *       fhAssocPtBkg;                   //! Trigger pT vs associated pT for background
 399   TH2F **      fhDeltaPhiDeltaEtaAssocPtBin;   //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle  phi as function eta difference, for different associated bins
 400   TH2F **      fhDeltaPhiAssocPtBin;           //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins
 401   TH2F **      fhDeltaPhiAssocPtBinDEta08;     //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins for Delta eta > 0.8
 402   TH2F **      fhDeltaPhiAssocPtBinDEta0 ;     //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins for Delta eta = 0
 403   TH2F **      fhDeltaPhiAssocPtBinHMPID;      //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins, track with HMPID
 404   TH2F **      fhDeltaPhiAssocPtBinHMPIDAcc;   //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi for different associated pt and vz bins, track with HMPIDAcc
 405   TH2F **      fhDeltaPhiBradAssocPtBin;       //![fNAssocPtBins*GetNZvertBin()] Trigger pT vs dPhi Brad (?) for different associated pt bins
 406   TH2F *       fhDeltaPhiBrad;                 //! Trigger pT vs dPhi Brad (?) for different associated pt bins
 407   TH2F **      fhXEAssocPtBin ;                //![fNAssocPtBins] Trigger pT vs xE for different associated pt bins
 408   TH2F **      fhZTAssocPtBin ;                //![fNAssocPtBins] Trigger pT vs zT for different associated pt bins
 409   TH2F **      fhXEVZ ;                        //![GetNZvertBin()] Trigger pT vs xE for different vz bins
 410   TH2F **      fhZTVZ ;                        //![GetNZvertBin()] Trigger pT vs zT for different vz bins
 411
 412   //trigger-neutral histograms
 413   TH2F *       fhDeltaPhiDeltaEtaNeutral ;     //! differences of eta and phi between trigger and neutral hadrons (pi0)
 414   TH2F *       fhPhiNeutral   ;                //! Phi distribution of neutral particles
 415   TH2F *       fhEtaNeutral   ;                //! Eta distribution of neutral particles
 416   TH2F *       fhDeltaPhiNeutral   ;           //! Difference of neutral particle phi and trigger particle  phi as function of  trigger particle pT
 417   TH2F *       fhDeltaEtaNeutral  ;            //! Difference of neutral particle eta and trigger particle  eta as function of  trigger particle pT
 418   TH2F *       fhDeltaPhiNeutralPt  ;          //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT
 419   TH2F *       fhDeltaPhiUeNeutralPt ;         //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT
 420   TH2F *       fhXENeutral  ;                  //! Trigger particle - neutral hadron momentum imbalance histogram
 421   TH2F *       fhXEUeNeutral  ;                //! Trigger particle - neutral hadron momentum imbalance histogram
 422   TH2F *       fhPtHbpXENeutral  ;             //! Trigger particle - neutral particle momentum HBP histogram
 423   TH2F *       fhPtHbpXEUeNeutral  ;           //! Trigger particle - underlying neutral hadron momentum HBP histogram
 424   TH2F *       fhZTNeutral  ;                  //! Trigger particle - neutral hadron momentum imbalance histogram
 425   TH2F *       fhZTUeNeutral  ;                //! Trigger particle - neutral hadron momentum imbalance histogram
 426   TH2F *       fhPtHbpZTNeutral  ;             //! Trigger particle - neutral particle momentum HBP histogram
 427   TH2F *       fhPtHbpZTUeNeutral  ;           //! Trigger particle - underlying neutral hadron momentum HBP histogram
 428
 429   // If several UE calculation is on,
 430   TH2F *       fhDeltaPhiUeLeftNeutral  ;      //! Difference of charged particle from underlying events phi and trigger particle  phi as function of neutral particle pT
 431   TH2F *       fhXEUeLeftNeutral  ;            //! Trigger particle -underlying neutral hadron momentum imbalance histogram
 432   TH2F *       fhPtHbpXEUeLeftNeutral  ;       //! Trigger particle -underlying neutral hadron momentum HBP histogram
 433   TH2F *       fhZTUeLeftNeutral  ;            //! Trigger particle -underlying neutral hadron momentum imbalance histogram
 434   TH2F *       fhPtHbpZTUeLeftNeutral  ;       //! Trigger particle -underlying neutral hadron momentum HBP histogram
 435
 436   // Pi0/Eta trigger correlation, recover input photons
 437   TH2F *       fhPtPi0DecayRatio ;             //! for pi0 trigger pt and ratio of decay photon pt
 438   TH2F *       fhDeltaPhiPi0DecayCharged  ;    //! Difference of charged particle phi and decay photon from pi0/eta trigger
 439   TH2F *       fhXEPi0DecayCharged ;           //! Trigger particle (decay from pi0/eta trigger)-charged hadron momentum imbalance histogram
 440   TH2F *       fhZTPi0DecayCharged ;           //! Trigger particle (decay from pi0/eta trigger)-charged hadron momentum imbalance histogram
 441
 442   TH2F *       fhDeltaPhiPi0DecayNeutral  ;    //! Difference of neutral particle phi and decay photon from pi0/eta trigger
 443   TH2F *       fhXEPi0DecayNeutral ;           //! Trigger particle (decay from pi0/eta trigger)-neutral hadron momentum imbalance histogram
 444   TH2F *       fhZTPi0DecayNeutral ;           //! Trigger particle (decay from pi0/eta trigger)-neutral hadron momentum imbalance histogram
 445
 446   // Decay photon trigger correlation
 447   TH2F *       fhDeltaPhiDecayCharged[4]  ;    //! Difference of charged particle phi and photon decay trigger
 448   TH2F *       fhXEDecayCharged[4] ;           //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram
 449   TH2F *       fhZTDecayCharged[4] ;           //! Trigger particle (decay from pi0)-charged hadron momentum imbalance histogram
 450   TH2F **      fhDeltaPhiDecayChargedAssocPtBin;//![fNAssocPtBins*GetNZvertBin()] Tagged as decay (fDecayBits[0]) Trigger pT vs dPhi for different associated pt bins
 451
 452   // If the data is MC, correlation with generated particles
 453   // check the origin of the cluster : decay photon (pi0, eta, other), merged photon (pi0),
 454   // hadron, rest of photons (prompt, FSR, ISR)
 455   TH1F *       fhMCPtTrigger[fgkNmcTypes];               //! MC pure pT distribution of trigger particles
 456   TH2F *       fhMCPhiTrigger[fgkNmcTypes];              //! MC pure Phi distribution of trigger particles
 457   TH2F *       fhMCEtaTrigger[fgkNmcTypes];              //! MC pure Eta distribution of trigger particles
 458   TH1F *       fhMCPtTriggerNotLeading[fgkNmcTypes];     //! MC pure pT distribution of trigger not leading particles
 459   TH2F *       fhMCPhiTriggerNotLeading[fgkNmcTypes];    //! MC pure Phi distribution of trigger not leading particles
 460   TH2F *       fhMCEtaTriggerNotLeading[fgkNmcTypes];    //! MC pure Eta distribution of trigger not leading particles
 461   TH2F *       fhMCEtaCharged[fgkNmcTypes];              //! MC pure particles charged primary pt vs eta (both associated)
 462   TH2F *       fhMCPhiCharged[fgkNmcTypes];              //! MC pure particles charged primary pt vs phi (both associated)
 463   TH2F *       fhMCDeltaEtaCharged[fgkNmcTypes];         //! MC pure particles charged trigger primary pt vs delta eta (associated-trigger)
 464   TH2F *       fhMCDeltaPhiCharged[fgkNmcTypes];         //! MC pure particles charged trigger primary pt vs delta phi (associated-trigger)
 465   TH2F *       fhMCDeltaPhiDeltaEtaCharged[fgkNmcTypes]; //! MC pure particles charged associated primary pt vs delta phi (associated-trigger), in away side
 466   TH2F *       fhMCDeltaPhiChargedPt[fgkNmcTypes];       //! MC pure particles charged delta phi vs delta eta (associated-trigger)
 467   TH2F *       fhMCPtXECharged[fgkNmcTypes];             //! MC pure particles charged trigger primary pt vs xE
 468   TH2F *       fhMCPtXEUeCharged[fgkNmcTypes];           //! MC pure particles charged trigger primary pt vs xE (underlying event)
 469   TH2F *       fhMCPtXEUeLeftCharged[fgkNmcTypes];       //! MC pure particles charged trigger primary pt vs xE (underlying event,left cone)
 470   TH2F *       fhMCPtHbpXECharged[fgkNmcTypes];          //! MC pure particles charged trigger primary pt vs ln(1/xE)
 471   TH2F *       fhMCPtHbpXEUeCharged[fgkNmcTypes];        //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event)
 472   TH2F *       fhMCPtHbpXEUeLeftCharged[fgkNmcTypes];    //! MC pure particles charged trigger primary pt vs ln(1/xE) (underlying event, left cone)
 473   TH1F *       fhMCUePart[fgkNmcTypes];                  //! MC pure UE particles distribution vs pt trig
 474   TH2F *       fhMCPtZTCharged[fgkNmcTypes];             //! MC pure particles charged trigger primary pt vs zT
 475   TH2F *       fhMCPtZTUeCharged[fgkNmcTypes];           //! MC pure particles charged trigger primary pt vs zT (underlying event)
 476   TH2F *       fhMCPtZTUeLeftCharged[fgkNmcTypes];       //! MC pure particles charged trigger primary pt vs zT (underlying event, left cone)
 477   TH2F *       fhMCPtHbpZTCharged[fgkNmcTypes];          //! MC pure particles charged trigger primary pt vs ln(1/zT)
 478   TH2F *       fhMCPtHbpZTUeCharged[fgkNmcTypes];        //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event)
 479   TH2F *       fhMCPtHbpZTUeLeftCharged[fgkNmcTypes];    //! MC pure particles charged trigger primary pt vs ln(1/zT) (underlying event, left cone)
 480   TH2F *       fhMCPtTrigPout[fgkNmcTypes];              //! MC pure particles charged trigger primary pt vs pOut
 481   TH2F *       fhMCPtAssocDeltaPhi[fgkNmcTypes];         //! MC pure particles charged associated primary pt vs delta phi (associated-trigger)
 482
 483   // Mixing
 484   TH1I *       fhNEventsTrigger;               //! number of analyzed triggered events
 485   TH2F *       fhNtracksMB;                    //! total number of tracks in MB events
 486   TH2F *       fhNclustersMB;                  //! total number of clusters in MB events
 487   TH2F *       fhMixDeltaPhiCharged  ;         //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
 488   TH2F *       fhMixDeltaPhiDeltaEtaCharged  ; //! Difference of charged particle phi and trigger particle  phi as function eta difference
 489   TH2F *       fhMixXECharged;                 //! xE for mixed event
 490   TH2F *       fhMixXEUeCharged;               //! xE for mixed event in Ue region
 491   TH2F *       fhMixHbpXECharged;              //! ln(1/xE) for mixed event
 492   TH2F **      fhMixDeltaPhiChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, for different associated bins
 493   TH2F **      fhMixDeltaPhiChargedAssocPtBinDEta08;   //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, for different associated bins, delta eta > 0.8
 494   TH2F **      fhMixDeltaPhiChargedAssocPtBinDEta0;    //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT, for different associated bins, delta eta = 0
 495   TH2F **      fhMixDeltaPhiDeltaEtaChargedAssocPtBin; //![fNAssocPtBins*GetNZvertBin()] Difference of charged particle phi and trigger particle  phi as function eta difference, for different associated bins
 496
 497   TH1I *       fhEventBin;                     //! Number of triggers in a particular event bin (cen,vz,rp)
 498   TH1I *       fhEventMixBin;                  //! Number of triggers mixed in a particular bin (cen,vz,rp)
 499   TH1I *       fhEventMBBin;                   //! Number of MB events in a particular bin (cen,vz,rp)
 500
 501   AliAnaParticleHadronCorrelation(              const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor
 502   AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ; // cpy assignment
 503
 504   ClassDef(AliAnaParticleHadronCorrelation,35)
 505 } ;
 506
 507
 508 #endif //ALIANAPARTICLEHADRONCORRELATION_H
 509
 510
 511