]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG4/PartCorrDep/AliAnaParticleHadronCorrelation.h
Selection of absolute leading for triggering, MC histograms added, possibility of...
[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaParticleHadronCorrelation.h
index 0ea7b603b98f00e52fa6b772c94870264f46f3dc..668b313ea7b92925dfe9b676b2eb2fde71f55a76 100755 (executable)
-#ifndef ALIANAPARTICLEHADRONCORRELATION_H\r
-#define ALIANAPARTICLEHADRONCORRELATION_H\r
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
- * See cxx source for full Copyright notice     */\r
-/* $Id:  $ */\r
-\r
-//_________________________________________________________________________\r
-// Class that contains the algorithm for the analysis of particle - hadron correlations\r
-// Particle (for example direct gamma) must be found in a previous analysis \r
-//-- Author: Gustavo Conesa (INFN-LNF)\r
-\r
-//  Modified by Yaxian Mao:\r
-// 1. add the UE subtraction for corrlation study\r
-// 2. change the correlation variable\r
-// 3. Only use leading particle(cluster/track) as trigger for correlation (2010/07/02)\r
-// 4. Make decay photon-hadron correlations where decay contribute pi0 mass (2010/09/09)\r
-// 5. fill the pout to extract kt at the end, also to study charge asymmetry(2010/10/06) \r
-// 6. Add the possibality for event selection analysis based on vertex and multiplicity bins (10/10/2010)\r
-// 7. change the way of delta phi cut for UE study due to memory issue (reduce histograms)\r
-\r
-// --- ROOT system ---\r
-//class TH3D;\r
-\r
-// --- Analysis system ---\r
-#include "AliAnaPartCorrBaseClass.h"\r
-class AliAODPWG4ParticleCorrelation ;\r
-\r
-class AliAnaParticleHadronCorrelation : public AliAnaPartCorrBaseClass {\r
-  \r
- public: \r
-  AliAnaParticleHadronCorrelation() ; // default ctor\r
-  virtual ~AliAnaParticleHadronCorrelation() {;} //virtual dtor\r
- private:  \r
-  AliAnaParticleHadronCorrelation(const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor\r
-  AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ;//cpy assignment\r
-\r
- public:\r
-  \r
-  TObjString * GetAnalysisCuts();\r
-  TList * GetCreateOutputObjects();\r
-  \r
-  Double_t GetDeltaPhiMaxCut() const {return fDeltaPhiMaxCut ; }\r
-  Double_t GetDeltaPhiMinCut() const {return fDeltaPhiMinCut ; }\r
-  void SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)\r
-  {fDeltaPhiMaxCut =phimax;  fDeltaPhiMinCut =phimin;}\r
\r
-  Double_t GetUeDeltaPhiMaxCut() const {return fUeDeltaPhiMaxCut ; }\r
-  Double_t GetUeDeltaPhiMinCut() const {return fUeDeltaPhiMinCut ; }\r
-  void SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax)\r
-  {fUeDeltaPhiMaxCut =uephimax;  fUeDeltaPhiMinCut =uephimin;}\r
-  Bool_t IsSeveralUEOn() const {return fMakeSeveralUE ; }\r
-  void SwitchOnSeveralUECalculation()  { fMakeSeveralUE = kTRUE;}\r
-  void SwitchOffSeveralUECalculation() { fMakeSeveralUE = kFALSE;}\r
-\r
-  // Do trigger-neutral correlation\r
-  Bool_t DoNeutralCorr() const {return fNeutralCorr ; }\r
-  void SwitchOnNeutralCorr()  { fNeutralCorr = kTRUE;}\r
-  void SwitchOffNeutralCorr() { fNeutralCorr = kFALSE;}  \r
-  \r
-  // Do decay-hadron correlation if it is pi0 trigger\r
-  Bool_t IsPi0Trigger() const {return fPi0Trigger ; }\r
-  void SwitchOnDecayCorr()  { fPi0Trigger = kTRUE;}\r
-  void SwitchOffDecayCorr() { fPi0Trigger = kFALSE;}  \r
-  \r
-  Bool_t OnlyIsolated() const {return fSelectIsolated ; }\r
-  void SelectIsolated(Bool_t select) {fSelectIsolated = select ; }\r
-\r
-  void InitParameters();\r
-  \r
-  void Print(const Option_t * opt) const;\r
-  \r
-  void MakeChargedCorrelation(AliAODPWG4ParticleCorrelation * aodParticle,TObjArray* const pl,   const Bool_t bFillHisto) ;\r
-  void MakeNeutralCorrelation(AliAODPWG4ParticleCorrelation * aodParticle,TObjArray* const pl, const Bool_t bFillHisto) ;\r
-\r
-  //void MakeNeutralCorrelationFillAOD(AliAODPWG4ParticleCorrelation* const aodParticle, TObjArray* const pl, TString detector)  ;\r
-  //void MakeNeutralCorrelationFillHistograms(AliAODPWG4ParticleCorrelation* const aodParticle)  ;\r
-       \r
-  void MakeAnalysisFillAOD()  ;\r
-  void MakeAnalysisFillHistograms() ; \r
-  \r
-  //Bool_t SelectCluster(AliVCluster * calo, Double_t *vertex, TLorentzVector & mom, Int_t & pdg) ;\r
-  \r
-  void SetPi0AODBranchName(TString pi0list) {fPi0AODBranchName = pi0list;}\r
-  \r
- private:\r
-  \r
-  Double_t   fDeltaPhiMaxCut ;      // Minimum Delta Phi Gamma-Hadron\r
-  Double_t   fDeltaPhiMinCut ;      // Maximum Delta Phi Gamma-Hadron\r
-  Bool_t     fSelectIsolated ;      // Select only trigger particles isolated\r
-  Bool_t     fMakeSeveralUE ;       // Do analysis for several underlying events contribution\r
-  Double_t   fUeDeltaPhiMaxCut ;    // Minimum Delta Phi Gamma-Underlying Hadron\r
-  Double_t   fUeDeltaPhiMinCut ;    // Maximum Delta Phi Gamma-Underlying Hadron\r
-  TString    fPi0AODBranchName;     // Name of AOD branch with pi0, not trigger\r
-  Bool_t     fNeutralCorr ;          // switch the analysis with neutral particles\r
-  Bool_t     fPi0Trigger ;          // switch the analysis with decay photon from pi0 trigger\r
-  \r
-  //Histograms\r
-//  TH2F * fhNclustersNtracks; //charge and cluster multiplicity distribution\r
-  //leading particles \r
-  TH1F * fhPtLeading;         //! pT distribution of leading particles\r
-  TH2F * fhPhiLeading;        //! phi distribution vs pT of leading particles\r
-  TH2F * fhEtaLeading;        //! eta distribution vs pT of leading particles\r
-  \r
-  //trigger-charged histograms\r
-  TH2F * fhDeltaPhiDeltaEtaCharged ; //! differences of eta and phi between trigger and charged hadrons\r
-  TH2F * fhPhiCharged  ;         //! Phi distribution of charged particles\r
-  TH2F * fhEtaCharged  ;         //! Eta distribution of charged particles\r
-  TH2F * fhDeltaPhiCharged  ;    //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT\r
-  TH2F * fhDeltaEtaCharged  ;    //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT\r
-  TH2F * fhDeltaPhiChargedPt  ;  //! Difference of charged particle phi and trigger particle  phi as function of charged particle pT\r
-  TH2F * fhDeltaPhiUeChargedPt ; //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT\r
-  TH2F * fhPtImbalanceCharged  ;   //! Trigger particle -charged hadron momentim imbalance histogram\r
-  TH2F * fhPtImbalanceUeCharged  ; //! Trigger particle -underlying charged hadron momentim imbalance histogram  \r
-  TH2F * fhPtImbalancePosCharged  ;   //! Trigger particle -positive charged hadron momentim imbalance histogram\r
-  TH2F * fhPtImbalanceNegCharged  ; //! Trigger particle -negative charged hadron momentim imbalance histogram \r
-  //with different imblance varible defination HBP distribution\r
-  TH2F * fhPtHbpCharged  ;   //! Trigger particle -charged hadron momentim HBP histogram\r
-  TH2F * fhPtHbpUeCharged  ; //! Trigger particle -underlying charged hadron momentim HBP histogram  \r
-  \r
-  //if several UE calculation is on, most useful for jet-jet events contribution\r
-  TH2F * fhDeltaPhiUeLeftCharged  ;    //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT\r
-  TH2F * fhDeltaPhiUeRightCharged  ;   //! Difference of charged particle from underlying events phi and trigger particle  phi \r
-  TH2F * fhPtImbalanceUeLeftCharged  ; //! Trigger particle -underlying charged hadron momentim imbalance histogram \r
-  TH2F * fhPtImbalanceUeRightCharged ; //! Trigger particle -underlying charged hadron momentim imbalance histogram  \r
-  TH2F * fhPtHbpUeLeftCharged  ;       //! Trigger particle -underlying charged hadron momentim HBP histogram \r
-  TH2F * fhPtHbpUeRightCharged  ;      //! Trigger particle -underlying charged hadron momentim HBP histogram  \r
-\r
-  //for pout and kt extraction\r
-  TH2F * fhPtTrigPout  ; // Pout =associated pt*sin(delta phi) distribution vs trigger pt vs associated pt\r
-  TH2F * fhPtAssocDeltaPhi  ; // Pout =associated pt*sin(delta phi) distribution\r
- // TH3D * fhUePoutPtTrigPtAssoc  ; // UE Pout =associated pt*sin(delta phi) distribution vs trigger pt vs associated pt\r
-  TH2F * fhPtTrigCharged ; //trigger and correlated particl pt, to be used for mean value for kt       \r
-  \r
-  //if different multiplicity analysis asked\r
-  TH2F ** fhTrigDeltaPhiCharged ; //! differences of phi between trigger and charged hadrons\r
-  TH2F ** fhTrigDeltaEtaCharged ; //! differences of eta between trigger and charged hadrons\r
-  TH2F ** fhTrigCorr  ;    //! Trigger particle -charged hadron momentim imbalance histogram\r
-  TH2F ** fhTrigUeCorr  ;    //! Trigger particle -UE charged hadron momentim imbalance histogram\r
-  \r
-  //trigger-neutral histograms\r
-  TH2F * fhDeltaPhiDeltaEtaNeutral ; //! differences of eta and phi between trigger and neutral hadrons (pi0)\r
-  TH2F * fhPhiNeutral   ;        //! Phi distribution of neutral particles  \r
-  TH2F * fhEtaNeutral   ;        //! Eta distribution of neutral particles\r
-  TH2F * fhDeltaPhiNeutral   ;   //! Difference of neutral particle phi and trigger particle  phi as function of  trigger particle pT\r
-  TH2F * fhDeltaEtaNeutral  ;    //! Difference of neutral particle eta and trigger particle  eta as function of  trigger particle pT\r
-  TH2F * fhDeltaPhiNeutralPt  ;  //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT\r
-  TH2F * fhDeltaPhiUeNeutralPt ; //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT  \r
-  TH2F * fhPtImbalanceNeutral  ;   //! Trigger particle - neutral hadron momentum imbalance histogram \r
-  TH2F * fhPtImbalanceUeNeutral  ; //! Trigger particle - neutral hadron momentum imbalance histogram \r
-  //with different imblance varible defination HBP distribution  \r
-  TH2F * fhPtHbpNeutral  ;   //! Trigger particle -neutral particle momentim HBP histogram\r
-  TH2F * fhPtHbpUeNeutral  ; //! Trigger particle -underlying neutral hadron momentim HBP histogram  \r
-\r
-  //if several UE calculation is on, most useful for jet-jet events contribution\r
-  TH2F * fhDeltaPhiUeLeftNeutral  ;    //! Difference of charged particle from underlying events phi and trigger particle  phi as function of neutral particle pT\r
-  TH2F * fhDeltaPhiUeRightNeutral  ;   //! Difference of charged particle from underlying events phi and trigger particle  phi \r
-  TH2F * fhPtImbalanceUeLeftNeutral  ; //! Trigger particle -underlying neutral hadron momentim imbalance histogram \r
-  TH2F * fhPtImbalanceUeRightNeutral ; //! Trigger particle -underlying neutral hadron momentim imbalance histogram \r
-  TH2F * fhPtHbpUeLeftNeutral  ;       //! Trigger particle -underlying neutral hadron momentim HBP histogram \r
-  TH2F * fhPtHbpUeRightNeutral  ;      //! Trigger particle -underlying neutral hadron momentim HBP histogram \r
-  \r
-  //for decay photon trigger correlation\r
-  TH2F * fhPtPi0DecayRatio ;          //! for pi0 pt and ratio of decay photon pt\r
-  TH2F * fhDeltaPhiDecayCharged  ;   //! Difference of charged particle phi and decay trigger\r
-  TH2F * fhPtImbalanceDecayCharged ; //! Trigger particle (decay from pi0)-charged hadron momentim imbalance histogram    \r
-  TH2F * fhDeltaPhiDecayNeutral  ;   //! Difference of neutral particle phi and decay trigger\r
-  TH2F * fhPtImbalanceDecayNeutral ; //! Trigger particle (decay from pi0)-neutral hadron momentim imbalance histogram  \r
-  \r
-       \r
-  ClassDef(AliAnaParticleHadronCorrelation,6)\r
-} ;\r
\r
-\r
-#endif //ALIANAPARTICLEHADRONCORRELATION_H\r
-\r
-\r
-\r
+#ifndef ALIANAPARTICLEHADRONCORRELATION_H
+#define ALIANAPARTICLEHADRONCORRELATION_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice     */
+/* $Id:  $ */
+
+//_________________________________________________________________________
+// Class that contains the algorithm for the analysis of particle - hadron correlations
+// Particle (for example direct gamma) must be found in a previous analysis 
+//-- Author: Gustavo Conesa (INFN-LNF)
+
+//  Modified by Yaxian Mao:
+// 1. add the UE subtraction for corrlation study
+// 2. change the correlation variable
+// 3. Only use leading particle(cluster/track) as trigger for correlation (2010/07/02)
+// 4. Make decay photon-hadron correlations where decay contribute pi0 mass (2010/09/09)
+// 5. fill the pout to extract kt at the end, also to study charge asymmetry(2010/10/06) 
+// 6. Add the possibility for event selection analysis based on vertex and multiplicity bins (10/10/2010)
+// 7. change the way of delta phi cut for UE study due to memory issue (reduce histograms)
+// 8. Add the possibility to request the absolute leading particle at the near side or not, set trigger bins, general clean-up (08/2011)
+
+// --- ROOT system ---
+//class TH3D;
+
+// --- Analysis system ---
+#include "AliAnaPartCorrBaseClass.h"
+class AliAODPWG4ParticleCorrelation ;
+
+class AliAnaParticleHadronCorrelation : public AliAnaPartCorrBaseClass {
+  
+ public: 
+  AliAnaParticleHadronCorrelation() ; // default ctor
+  virtual ~AliAnaParticleHadronCorrelation() {;} //virtual dtor
+ private:  
+  AliAnaParticleHadronCorrelation(const AliAnaParticleHadronCorrelation & ph) ; // cpy ctor
+  AliAnaParticleHadronCorrelation & operator = (const AliAnaParticleHadronCorrelation & ph) ;//cpy assignment
+
+ public:
+  
+  // General methods
+  
+  TObjString * GetAnalysisCuts();
+  
+  TList      * GetCreateOutputObjects();
+  
+  void         InitParameters();
+  
+  void         MakeAnalysisFillAOD()  ;
+  
+  void         MakeAnalysisFillHistograms() ; 
+  
+  void         Print(const Option_t * opt) const;
+  
+  // Main analysis methods
+  
+  Bool_t       MakeChargedCorrelation  (AliAODPWG4ParticleCorrelation * aodParticle, TObjArray* pl, const Bool_t bFillHisto) ;
+  
+  Bool_t       MakeNeutralCorrelation  (AliAODPWG4ParticleCorrelation * aodParticle, TObjArray* pl, const Bool_t bFillHisto) ;
+  
+  void         MakeMCChargedCorrelation(AliAODPWG4ParticleCorrelation *aodParticle);
+  
+  
+  // Parameter setter and getter
+  
+  Double_t     GetDeltaPhiMaxCut()       const { return fDeltaPhiMaxCut   ; }
+  Double_t     GetDeltaPhiMinCut()       const { return fDeltaPhiMinCut   ; }
+  Double_t     GetUeDeltaPhiMaxCut()     const { return fUeDeltaPhiMaxCut ; }
+  Double_t     GetUeDeltaPhiMinCut()     const { return fUeDeltaPhiMinCut ; }
+  
+  void         SetDeltaPhiCutRange(Double_t phimin, Double_t phimax)
+    { fDeltaPhiMaxCut   = phimax ;   fDeltaPhiMinCut   = phimin   ; }
+  void         SetUeDeltaPhiCutRange(Double_t uephimin, Double_t uephimax)
+    { fUeDeltaPhiMaxCut = uephimax;  fUeDeltaPhiMinCut = uephimin ; }
+  
+  Bool_t       IsSeveralUEOn()           const { return fMakeSeveralUE    ; }
+  void         SwitchOnSeveralUECalculation()  { fMakeSeveralUE = kTRUE   ; }
+  void         SwitchOffSeveralUECalculation() { fMakeSeveralUE = kFALSE  ; }
+
+  // Do trigger-neutral correlation
+  Bool_t       DoNeutralCorr()           const { return fNeutralCorr      ; }
+  void         SwitchOnNeutralCorr()           { fNeutralCorr = kTRUE     ; }
+  void         SwitchOffNeutralCorr()          { fNeutralCorr = kFALSE    ; }  
+  
+  // Taking the absolute leading as the trigger or not
+  Bool_t       DoAbsoluteLeading()       const { return fMakeAbsoluteLeading   ; }
+  void         SwitchOnAbsoluteLeading()       { fMakeAbsoluteLeading = kTRUE  ; }
+  void         SwitchOffAbsoluteLeading()      { fMakeAbsoluteLeading = kFALSE ; }
+  
+  // Do decay-hadron correlation if it is pi0 trigger
+  Bool_t       IsPi0Trigger()            const { return fPi0Trigger       ; }
+  void         SwitchOnDecayCorr()             { fPi0Trigger = kTRUE      ; }
+  void         SwitchOffDecayCorr()            { fPi0Trigger = kFALSE     ; }  
+  
+  Bool_t       OnlyIsolated()            const { return fSelectIsolated   ; }
+  void         SelectIsolated(Bool_t s)        { fSelectIsolated   = s    ; }
+
+  void         SetPi0AODBranchName(TString n)  { fPi0AODBranchName = n    ; }
+  
+  void         SetNTriggerPtBins(Int_t n) ;     
+  void         SetTriggerPtBinLimit(Int_t ibin, Float_t pt) ;
+                
+ private:
+  
+  Double_t     fDeltaPhiMaxCut ;               // Minimum Delta Phi Gamma-Hadron
+  Double_t     fDeltaPhiMinCut ;               // Maximum Delta Phi Gamma-Hadron
+  Bool_t       fSelectIsolated ;               // Select only trigger particles isolated
+  Bool_t       fMakeSeveralUE ;                // Do analysis for several underlying events contribution
+  Double_t     fUeDeltaPhiMaxCut ;             // Minimum Delta Phi Gamma-Underlying Hadron
+  Double_t     fUeDeltaPhiMinCut ;             // Maximum Delta Phi Gamma-Underlying Hadron
+  TString      fPi0AODBranchName;              // Name of AOD branch with pi0, not trigger
+  Bool_t       fNeutralCorr ;                  // switch the analysis with neutral particles
+  Bool_t       fPi0Trigger ;                   // switch the analysis with decay photon from pi0 trigger
+  Bool_t       fMakeAbsoluteLeading ;          // requesting absolute leading while it is cluster triggers
+  Int_t        fLeadingTriggerIndex ;          // Store here per event the trigger index, to avoid too many loops
+  
+  Int_t        fNTriggerPtBins ;               // Number of trigger bins under study
+  Float_t      fTriggerPtBinLimit[10] ;        // Trigger pt limtis
+  
+  //Histograms
+
+  //leading particles 
+  TH1F *       fhPtLeading;                    //! pT distribution of leading particles
+  TH2F *       fhPhiLeading;                   //! phi distribution vs pT of leading particles
+  TH2F *       fhEtaLeading;                   //! eta distribution vs pT of leading particles
+  
+  //trigger-charged histograms
+  TH2F *       fhDeltaPhiDeltaEtaCharged ;     //! differences of eta and phi between trigger and charged hadrons
+  TH2F *       fhPhiCharged  ;                 //! Phi distribution of charged particles
+  TH2F *       fhEtaCharged  ;                 //! Eta distribution of charged particles
+  TH2F *       fhDeltaPhiCharged  ;            //! Difference of charged particle phi and trigger particle  phi as function of  trigger particle pT
+  TH2F *       fhDeltaEtaCharged  ;            //! Difference of charged particle eta and trigger particle  eta as function of  trigger particle pT
+  TH2F *       fhDeltaPhiChargedPt  ;          //! Difference of charged particle phi and trigger particle  phi as function of charged particle pT
+  TH2F *       fhDeltaPhiUeChargedPt ;         //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT
+  TH2F *       fhPtImbalanceCharged  ;         //! Trigger particle -charged hadron momentim imbalance histogram
+  TH2F *       fhPtImbalanceUeCharged  ;       //! Trigger particle -underlying charged hadron momentum imbalance histogram  
+  TH2F *       fhPtImbalancePosCharged  ;      //! Trigger particle -positive charged hadron momentum imbalance histogram
+  TH2F *       fhPtImbalanceNegCharged  ;      //! Trigger particle -negative charged hadron momentum imbalance histogram 
+  
+  //with different imblance varible defination HBP distribution
+  TH2F *       fhPtHbpCharged  ;               //! Trigger particle -charged hadron momentim HBP histogram
+  TH2F *       fhPtHbpUeCharged  ;             //! Trigger particle -underlying charged hadron momentim HBP histogram  
+  
+  //if several UE calculation is on, most useful for jet-jet events contribution
+  TH2F *       fhDeltaPhiUeLeftCharged  ;      //! Difference of charged particle from underlying events phi and trigger particle  phi as function of charged particle pT
+  TH2F *       fhDeltaPhiUeRightCharged  ;     //! Difference of charged particle from underlying events phi and trigger particle  phi 
+  TH2F *       fhPtImbalanceUeLeftCharged  ;   //! Trigger particle -underlying charged hadron momentim imbalance histogram 
+  TH2F *       fhPtImbalanceUeRightCharged ;   //! Trigger particle -underlying charged hadron momentim imbalance histogram  
+  TH2F *       fhPtHbpUeLeftCharged  ;         //! Trigger particle -underlying charged hadron momentim HBP histogram 
+  TH2F *       fhPtHbpUeRightCharged  ;        //! Trigger particle -underlying charged hadron momentim HBP histogram  
+
+  //for pout and kt extraction
+  TH2F *       fhPtTrigPout  ;                 //! Pout =associated pt*sin(delta phi) distribution vs trigger pt vs associated pt
+  TH2F *       fhPtAssocDeltaPhi  ;            //! Pout =associated pt*sin(delta phi) distribution
+  TH2F *       fhPtTrigCharged ;               //! trigger and correlated particl pt, to be used for mean value for kt 
+  
+  //if different multiplicity analysis asked
+  TH2F **      fhTrigDeltaPhiCharged ;         //! differences of phi between trigger and charged hadrons
+  TH2F **      fhTrigDeltaEtaCharged ;         //! differences of eta between trigger and charged hadrons
+  TH2F **      fhTrigCorr  ;                   //! Trigger particle -charged hadron momentim imbalance histogram
+  TH2F **      fhTrigUeCorr  ;                 //! Trigger particle -UE charged hadron momentim imbalance histogram
+    
+  TH1F *       fhTrigPt[10] ;                  //! add explanation
+  TH1F *       fhDphiTrigPtAssocPt[10];        //! add explanation
+  TH1F *       fhAssocPtBkg[10];               //! add explanation
+  TH1F *       fhXE[10] ;                      //! add explanation
+  
+  //trigger-neutral histograms
+  TH2F *       fhDeltaPhiDeltaEtaNeutral ;     //! differences of eta and phi between trigger and neutral hadrons (pi0)
+  TH2F *       fhPhiNeutral   ;                //! Phi distribution of neutral particles  
+  TH2F *       fhEtaNeutral   ;                //! Eta distribution of neutral particles
+  TH2F *       fhDeltaPhiNeutral   ;           //! Difference of neutral particle phi and trigger particle  phi as function of  trigger particle pT
+  TH2F *       fhDeltaEtaNeutral  ;            //! Difference of neutral particle eta and trigger particle  eta as function of  trigger particle pT
+  TH2F *       fhDeltaPhiNeutralPt  ;          //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT
+  TH2F *       fhDeltaPhiUeNeutralPt ;         //! Difference of neutral particle phi and trigger particle  phi as function of neutral particle particle pT  
+  TH2F *       fhPtImbalanceNeutral  ;         //! Trigger particle - neutral hadron momentum imbalance histogram 
+  TH2F *       fhPtImbalanceUeNeutral  ;       //! Trigger particle - neutral hadron momentum imbalance histogram 
+  
+  //with different imblance varible defination HBP distribution  
+  TH2F *       fhPtHbpNeutral  ;               //! Trigger particle -neutral particle momentim HBP histogram
+  TH2F *       fhPtHbpUeNeutral  ;             //! Trigger particle -underlying neutral hadron momentim HBP histogram  
+
+  //if several UE calculation is on, most useful for jet-jet events contribution
+  TH2F *       fhDeltaPhiUeLeftNeutral  ;      //! Difference of charged particle from underlying events phi and trigger particle  phi as function of neutral particle pT
+  TH2F *       fhDeltaPhiUeRightNeutral  ;     //! Difference of charged particle from underlying events phi and trigger particle  phi 
+  TH2F *       fhPtImbalanceUeLeftNeutral  ;   //! Trigger particle -underlying neutral hadron momentim imbalance histogram 
+  TH2F *       fhPtImbalanceUeRightNeutral ;   //! Trigger particle -underlying neutral hadron momentim imbalance histogram 
+  TH2F *       fhPtHbpUeLeftNeutral  ;         //! Trigger particle -underlying neutral hadron momentim HBP histogram 
+  TH2F *       fhPtHbpUeRightNeutral  ;        //! Trigger particle -underlying neutral hadron momentim HBP histogram 
+  
+  //for decay photon trigger correlation
+  TH2F *       fhPtPi0DecayRatio ;             //! for pi0 pt and ratio of decay photon pt
+  TH2F *       fhDeltaPhiDecayCharged  ;       //! Difference of charged particle phi and decay trigger
+  TH2F *       fhPtImbalanceDecayCharged ;     //! Trigger particle (decay from pi0)-charged hadron momentim imbalance histogram    
+  TH2F *       fhDeltaPhiDecayNeutral  ;       //! Difference of neutral particle phi and decay trigger
+  TH2F *       fhPtImbalanceDecayNeutral ;     //! Trigger particle (decay from pi0)-neutral hadron momentim imbalance histogram  
+  
+  //if the data is MC, fill MC information
+  TH2F *       fh2phiLeadingParticle;          //! #phi resolution for triggers
+  TH1F *       fhMCLeadingCount;               //! add explanation
+  TH2F *       fhMCEtaCharged;                 //! add explanation
+  TH2F *       fhMCPhiCharged;                 //! add explanation
+  TH2F *       fhMCDeltaEtaCharged;            //! add explanation
+  TH2F *       fhMCDeltaPhiCharged;            //! add explanation
+  TH2F *       fhMCDeltaPhiDeltaEtaCharged;    //! add explanation
+  TH2F *       fhMCDeltaPhiChargedPt;          //! add explanation
+  TH2F *       fhMCPtImbalanceCharged;         //! add explanation
+  TH2F *       fhMCPtHbpCharged;               //! add explanation
+  TH2F *       fhMCPtTrigPout ;                //! add explanation
+  TH2F *       fhMCPtAssocDeltaPhi  ;          //! Pout =associated pt*sin(delta phi) distribution
+
+       
+  ClassDef(AliAnaParticleHadronCorrelation,7)
+} ;
+
+#endif //ALIANAPARTICLEHADRONCORRELATION_H
+
+
+