[u/mrichter/AliRoot.git] / PWGGA / CaloTrackCorrelations / AliAnaParticleIsolation.h

#ifndef ALIANAPARTICLEISOLATION_H
#define ALIANAPARTICLEISOLATION_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */

//_________________________________________________________________________

// Class for the analysis of particle isolation
// Input is selected particles put in AOD branch (AliAODPWG4ParticleCorrelation)
//
//  Class created from old AliPHOSGammaJet
//  (see AliRoot versions previous Release 4-09)

//-- Author: Gustavo Conesa (INFN-LNF)

// --- ROOT system ---
class TH2F;
class TList ;
class TObjString;

// --- ANALYSIS system ---
#include "AliAnaCaloTrackCorrBaseClass.h"
class AliAODPWG4Particle;
class AliAODPWG4ParticleCorrelation ;


class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {

 public:   
  AliAnaParticleIsolation() ; // default ctor
  virtual ~AliAnaParticleIsolation() { ; } //virtual dtor

  // Main general methods
    
  
  TObjString * GetAnalysisCuts() ;
  
  TList      * GetCreateOutputObjects() ;
  
  void         Init() ;

  void         InitParameters() ;
  
  void         MakeAnalysisFillAOD()  ;
  
  void         MakeAnalysisFillHistograms() ; 
  
  void         Print( const Option_t * opt ) const ;
 
  //Analysis specific methods 
  
  void         FillTrackMatchingShowerShapeControlHistograms(const Int_t clusterID, 
                                                             const Int_t nLocMax,
                                                             const Int_t mcTag ) ;
  
  void         MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ; 
  
  // Analysis Setters and Getters
  
  TString      GetCalorimeter()                const { return fCalorimeter       ; }
  Int_t        GetNCones()                     const { return fNCones            ; }
  Int_t        GetNPtThresFrac()               const { return fNPtThresFrac      ; }
  Float_t      GetConeSizes(Int_t i)           const { return fConeSizes[i]      ; }
  Float_t      GetPtThresholds(Int_t i)        const { return fPtThresholds[i]   ; }
  Float_t      GetPtFractions(Int_t i)         const { return fPtFractions[i]    ; }
  
  void         SetCalorimeter(TString & det)         { fCalorimeter     = det    ; }
  void         SetNCones(Int_t ncs)                  { fNCones          = ncs    ; }
  void         SetNPtThresFrac(Int_t npt)            { fNPtThresFrac    = npt    ; }
  void         SetConeSizes(Int_t i, Float_t r)      { fConeSizes[i]    = r      ; }
  void         SetPtThresholds(Int_t i, Float_t pt)  { fPtThresholds[i] = pt     ; }
  void         SetPtFractions(Int_t i, Float_t pt)   { fPtFractions[i]  = pt     ; } 
  
  Bool_t       IsReIsolationOn()               const { return fReMakeIC          ; }
  void         SwitchOnReIsolation()                 { fReMakeIC      = kTRUE    ; }
  void         SwitchOffReIsolation()                { fReMakeIC      = kFALSE   ; }
  
  Bool_t       IsSeveralIsolationOn()          const { return fMakeSeveralIC     ; }
  void         SwitchOnSeveralIsolation()            { fMakeSeveralIC = kTRUE    ; }
  void         SwitchOffSeveralIsolation()           { fMakeSeveralIC = kFALSE   ; }

  void         SwitchOnTMHistoFill()                 { fFillTMHisto   = kTRUE    ; }
  void         SwitchOffTMHistoFill()                { fFillTMHisto   = kFALSE   ; }
  
  void         SwitchOnSSHistoFill()                 { fFillSSHisto   = kTRUE    ; }
  void         SwitchOffSSHistoFill()                { fFillSSHisto   = kFALSE   ; }

  //Histogrammes setters and getters
  
  void         SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n)       {
    fHistoNPtSumBins = n ;    fHistoPtSumMax = max ;    fHistoPtSumMin = min ;     }
  
  Int_t        GetHistoNPtSumBins()            const { return fHistoNPtSumBins   ; }
  Float_t      GetHistoPtSumMin()              const { return fHistoPtSumMin     ; }
  Float_t      GetHistoPtSumMax()              const { return fHistoPtSumMax     ; }
  
  void         SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n)    {
    fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min  ; }
  
  Int_t        GetHistoNPtInConeBins()         const { return fHistoNPtInConeBins; }
  Float_t      GetHistoPtInConeMin()           const { return fHistoPtInConeMin  ; }
  Float_t      GetHistoPtInConeMax()           const { return fHistoPtInConeMax  ; }
  
 private:
  
  TString  fCalorimeter ;                         // Calorimeter where neutral particles in cone for isolation are;
  Bool_t   fReMakeIC ;                            // Do isolation analysis
  Bool_t   fMakeSeveralIC ;                       // Do analysis for different IC
  Bool_t   fFillTMHisto;                          // Fill track matching plots
  Bool_t   fFillSSHisto;                          // Fill Shower shape plots

  // Analysis data members for multiple cones and pt thresholds 
  Int_t    fNCones ;                              //! Number of cone sizes to test
  Int_t    fNPtThresFrac ;                        //! Number of ptThres and ptFrac to test
  
  Float_t  fConeSizes[5] ;                        //! Array with cones to test
  Float_t  fPtThresholds[5] ;                     //! Array with pt thresholds to test
  Float_t  fPtFractions[5] ;                      //! Array with pt thresholds to test
  
  TH1F*    fhPtThresIsolated[5][5] ;              //! Isolated particle with pt threshold 
  TH1F*    fhPtFracIsolated[5][5] ;               //! Isolated particle with pt threshold 
  TH2F*    fhPtSumIsolated[5] ;                   //! Isolated particle with threshold on cone pt sum
  
  //Histograms  
  
  TH1F *   fhEIso ;                               //! Number of isolated particles
  TH1F *   fhPtIso ;                              //! Number of isolated particles
  TH2F *   fhPhiIso ;                             //! Phi of isolated particles
  TH2F *   fhEtaIso ;                             //! eta of isolated particles
  TH2F *   fhEtaPhiIso ;                          //! eta vs phi of isolated particles
  TH2F *   fhEtaPhiNoIso ;                        //! eta vs phi of not isolated leading particles
  TH1F *   fhPtNoIso ;                            //! Number of not isolated leading particles
  TH1F *   fhPtDecayIso ;                         //! Number of isolated Pi0 decay particles (invariant mass tag)
  TH1F *   fhPtDecayNoIso ;                       //! Number of not isolated Pi0 decay leading particles (invariant mass tag)
  TH2F *   fhEtaPhiDecayIso ;                     //! eta vs phi of isolated Pi0 decay particles
  TH2F *   fhEtaPhiDecayNoIso ;                   //! eta vs phi of not isolated leading Pi0 decay particles
  TH2F *   fhConeSumPt ;                          //! Sum Pt in the cone
  TH2F *   fhPtInCone ;                           //! Particle Pt in the cone
  TH2F *   fhFRConeSumPt ;                        //! Sum Pt in the forward region cone (phi +90)
  TH2F *   fhPtInFRCone ;                         //! Particle Pt in the forward region cone (phi +90 ) 
  
  
  //MC
  TH1F *   fhPtIsoPrompt;                         //! Number of isolated prompt gamma 
  TH2F *   fhPhiIsoPrompt;                        //! Phi of isolated prompt gamma
  TH2F *   fhEtaIsoPrompt;                        //! eta of isolated prompt gamma
  TH1F *   fhPtThresIsolatedPrompt[5][5];         //! Isolated prompt gamma with pt threshold 
  TH1F *   fhPtFracIsolatedPrompt[5][5];          //! Isolated prompt gamma with pt frac
  TH2F *   fhPtSumIsolatedPrompt[5];              //! Isolated prompt gamma with threshold on cone pt sume
  TH1F *   fhPtIsoFragmentation;                  //! Number of isolated fragmentation gamma 
  TH2F *   fhPhiIsoFragmentation;                 //! Phi of isolated fragmentation gamma
  TH2F *   fhEtaIsoFragmentation;                 //! eta of isolated fragmentation gamma
  TH1F *   fhPtThresIsolatedFragmentation[5][5];  //! Isolated fragmentation gamma with pt threshold 
  TH1F *   fhPtFracIsolatedFragmentation[5][5];   //! Isolated fragmentation gamma with pt frac
  TH2F *   fhPtSumIsolatedFragmentation[5];       //! Isolated fragmentation gamma with threshold on cone pt sume
  TH1F *   fhPtIsoPi0Decay;                       //! Number of isolated pi0 decay gamma 
  TH2F *   fhPhiIsoPi0Decay;                      //! Phi of isolated pi0 decay gamma
  TH2F *   fhEtaIsoPi0Decay;                      //! eta of isolated pi0 decay gamma
  TH1F *   fhPtThresIsolatedPi0Decay[5][5];       //! Isolated pi0 decay gamma with pt threshold 
  TH1F *   fhPtFracIsolatedPi0Decay[5][5];        //! Isolated pi0 decay gamma with pt frac
  TH2F *   fhPtSumIsolatedPi0Decay[5];            //! Isolated pi0 decay gamma with threshold on cone pt sume
  TH1F *   fhPtIsoEtaDecay;                       //! Number of isolated eta decay gamma 
  TH2F *   fhPhiIsoEtaDecay;                      //! Phi of isolated eta decay gamma
  TH2F *   fhEtaIsoEtaDecay;                      //! eta of isolated eta decay gamma
  TH1F *   fhPtThresIsolatedEtaDecay[5][5];       //! Isolated eta decay gamma with pt threshold 
  TH1F *   fhPtFracIsolatedEtaDecay[5][5];        //! Isolated eta decay gamma with pt frac
  TH2F *   fhPtSumIsolatedEtaDecay[5];            //! Isolated eta fecay gamma with threshold on cone pt sume  
  TH1F *   fhPtIsoOtherDecay;                     //! Number of isolated other decay gamma 
  TH2F *   fhPhiIsoOtherDecay;                    //! Phi of isolated other decay gamma
  TH2F *   fhEtaIsoOtherDecay;                    //! eta of isolated other decay gamma
  TH1F *   fhPtThresIsolatedOtherDecay[5][5];     //! Isolated OtherDecay gamma with pt threshold 
  TH1F *   fhPtFracIsolatedOtherDecay[5][5];      //! Isolated OtherDecay gamma with pt frac
  TH2F *   fhPtSumIsolatedOtherDecay[5];          //! Isolated OtherDecay gamma with threshold on cone pt sume	
  TH1F *   fhPtIsoConversion;                     //! Number of isolated Conversion gamma 
  TH2F *   fhPhiIsoConversion;                    //! Phi of isolated Conversion gamma
  TH2F *   fhEtaIsoConversion;                    //! eta of isolated Conversion gamma
  TH1F *   fhPtThresIsolatedConversion[5][5];     //! Isolated Conversion gamma with pt threshold 
  TH1F *   fhPtFracIsolatedConversion[5][5];      //! Isolated Conversion gamma with pt frac
  TH2F *   fhPtSumIsolatedConversion[5];          //! Isolated Conversion gamma with threshold on cone pt sume
  TH1F *   fhPtIsoUnknown;                        //! Number of isolated Unknown 
  TH2F *   fhPhiIsoUnknown;                       //! Phi of isolated Unknown
  TH2F *   fhEtaIsoUnknown;                       //! eta of isolated Unknown
  TH1F *   fhPtThresIsolatedUnknown[5][5];        //! Isolated Unknown gamma with pt threshold 
  TH1F *   fhPtFracIsolatedUnknown[5][5];         //! Isolated Unknown gamma with pt frac
  TH2F *   fhPtSumIsolatedUnknown[5];             //! Isolated Unknown gamma with threshold on cone pt sume

  TH1F *   fhPtNoIsoPi0Decay;                     //! Number of not isolated leading pi0 decay gamma 
  TH1F *   fhPtNoIsoEtaDecay;                     //! Number of not isolated leading eta decay gamma 
  TH1F *   fhPtNoIsoOtherDecay;                   //! Number of not isolated leading other decay gamma 
  TH1F *   fhPtNoIsoPrompt;                       //! Number of not isolated leading prompt gamma 
  TH1F *   fhPtIsoMCPhoton;                       //! Number of isolated leading gamma 
  TH1F *   fhPtNoIsoMCPhoton;                     //! Number of not isolated leading gamma 
  TH1F *   fhPtNoIsoConversion;                   //! Number of not isolated leading conversion gamma 
  TH1F *   fhPtNoIsoFragmentation;                //! Number of not isolated leading fragmentation gamma 
  TH1F *   fhPtNoIsoUnknown;                      //! Number of not isolated leading hadrons 

  // Track matching studies
  TH2F *   fhTrackMatchedDEta     ;               //! Eta distance between track and cluster vs cluster E
  TH2F *   fhTrackMatchedDPhi     ;               //! Phi distance between track and cluster vs cluster E
  TH2F *   fhTrackMatchedDEtaDPhi ;               //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
  TH2F *   fhdEdx  ;                              //! matched track dEdx vs cluster E 
  TH2F *   fhEOverP;                              //! matched track E cluster over P track vs cluster E, after dEdx cut 
  TH2F *   fhTrackMatchedMCParticle;              //! Trace origin of matched particle

  // Shower Shape histograms
  TH2F *   fhELambda0;                            //! Shower shape of isolated photons
  TH2F *   fhELambda1;                            //! Shower shape of isolated photons
  TH2F *   fhELambda0TRD;                         //! Shower shape of isolated photons, SM behind TRD
  TH2F *   fhELambda1TRD;                         //! Shower shape of isolated photons, SM behind TRD
  
  // Local maxima
  TH2F * fhNLocMax;                               //! number of maxima in selected clusters
  TH2F * fhELambda0LocMax1 ;                      //! E vs lambda0 of selected cluster, 1 local maxima in cluster 
  TH2F * fhELambda1LocMax1 ;                      //! E vs lambda1 of selected cluster, 1 local maxima in cluster 
  TH2F * fhELambda0LocMax2 ;                      //! E vs lambda0 of selected cluster, 2 local maxima in cluster 
  TH2F * fhELambda1LocMax2 ;                      //! E vs lambda1 of selected cluster, 2 local maxima in cluster
  TH2F * fhELambda0LocMaxN ;                      //! E vs lambda0 of selected cluster, N>2 local maxima in cluster 
  TH2F * fhELambda1LocMaxN ;                      //! E vs lambda1 of selected cluster, N>2 local maxima in cluster 
  
  //Histograms settings
  Int_t    fHistoNPtSumBins;                      // Number of bins in PtSum histograms
  Float_t  fHistoPtSumMax;                        // PtSum maximum in histogram
  Float_t  fHistoPtSumMin;	                      // PtSum minimum in histogram
  Int_t    fHistoNPtInConeBins;                   // Number of bins in PtInCone histogram
  Float_t  fHistoPtInConeMax;                     // PtInCone maximum in histogram
  Float_t  fHistoPtInConeMin;                     // PtInCone maximum in histogram 

  AliAnaParticleIsolation(              const AliAnaParticleIsolation & iso) ; // cpy ctor
  AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
  
  ClassDef(AliAnaParticleIsolation,12)
} ;


#endif //ALIANAPARTICLEISOLATION_H
Commit	Line	Data
1a31a9ab	1	#ifndef ALIANAPARTICLEISOLATION_H
	2	#define ALIANAPARTICLEISOLATION_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
1a31a9ab	5
	6	//_________________________________________________________________________
	7
	8	// Class for the analysis of particle isolation
	9	// Input is selected particles put in AOD branch (AliAODPWG4ParticleCorrelation)
	10	//
	11	// Class created from old AliPHOSGammaJet
	12	// (see AliRoot versions previous Release 4-09)
	13
	14	//-- Author: Gustavo Conesa (INFN-LNF)
	15
	16	// --- ROOT system ---
	17	class TH2F;
	18	class TList ;
	19	class TObjString;
	20
	21	// --- ANALYSIS system ---
745913ae	22	#include "AliAnaCaloTrackCorrBaseClass.h"
1a31a9ab	23	class AliAODPWG4Particle;
	24	class AliAODPWG4ParticleCorrelation ;
	25
	26
745913ae	27	class AliAnaParticleIsolation : public AliAnaCaloTrackCorrBaseClass {
1a31a9ab	28
	29	public:
	30	AliAnaParticleIsolation() ; // default ctor
803d06a8	31	virtual ~AliAnaParticleIsolation() { ; } //virtual dtor
1a31a9ab	32
1a31a9ab	33	// Main general methods
803d06a8	34
1a31a9ab	35
b5dbb99b	36	TObjString * GetAnalysisCuts() ;
1a31a9ab	37
b5dbb99b	38	TList * GetCreateOutputObjects() ;
b5dbb99b	39
03bae431	40	void Init() ;
03bae431	41
b5dbb99b	42	void InitParameters() ;
803d06a8	43
1a31a9ab	44	void MakeAnalysisFillAOD() ;
	45
	46	void MakeAnalysisFillHistograms() ;
	47
b5dbb99b	48	void Print( const Option_t * opt ) const ;
1a31a9ab	49
1a31a9ab	50	//Analysis specific methods
b5dbb99b	51
5c46c992	52	void FillTrackMatchingShowerShapeControlHistograms(const Int_t clusterID,
	53	const Int_t nLocMax,
	54	const Int_t mcTag ) ;
b5dbb99b	55
b5dbb99b	56	void MakeSeveralICAnalysis( AliAODPWG4ParticleCorrelation * ph ) ;
1a31a9ab	57
1a31a9ab	58	// Analysis Setters and Getters
1a31a9ab	59
b5dbb99b	60	TString GetCalorimeter() const { return fCalorimeter ; }
	61	Int_t GetNCones() const { return fNCones ; }
	62	Int_t GetNPtThresFrac() const { return fNPtThresFrac ; }
	63	Float_t GetConeSizes(Int_t i) const { return fConeSizes[i] ; }
	64	Float_t GetPtThresholds(Int_t i) const { return fPtThresholds[i] ; }
	65	Float_t GetPtFractions(Int_t i) const { return fPtFractions[i] ; }
1a31a9ab	66
b5dbb99b	67	void SetCalorimeter(TString & det) { fCalorimeter = det ; }
	68	void SetNCones(Int_t ncs) { fNCones = ncs ; }
	69	void SetNPtThresFrac(Int_t npt) { fNPtThresFrac = npt ; }
	70	void SetConeSizes(Int_t i, Float_t r) { fConeSizes[i] = r ; }
	71	void SetPtThresholds(Int_t i, Float_t pt) { fPtThresholds[i] = pt ; }
	72	void SetPtFractions(Int_t i, Float_t pt) { fPtFractions[i] = pt ; }
1a31a9ab	73
b5dbb99b	74	Bool_t IsReIsolationOn() const { return fReMakeIC ; }
	75	void SwitchOnReIsolation() { fReMakeIC = kTRUE ; }
	76	void SwitchOffReIsolation() { fReMakeIC = kFALSE ; }
1a31a9ab	77
b5dbb99b	78	Bool_t IsSeveralIsolationOn() const { return fMakeSeveralIC ; }
	79	void SwitchOnSeveralIsolation() { fMakeSeveralIC = kTRUE ; }
	80	void SwitchOffSeveralIsolation() { fMakeSeveralIC = kFALSE ; }
09273901	81
b5dbb99b	82	void SwitchOnTMHistoFill() { fFillTMHisto = kTRUE ; }
b5dbb99b	83	void SwitchOffTMHistoFill() { fFillTMHisto = kFALSE ; }
09273901	84
b5dbb99b	85	void SwitchOnSSHistoFill() { fFillSSHisto = kTRUE ; }
b5dbb99b	86	void SwitchOffSSHistoFill() { fFillSSHisto = kFALSE ; }
09273901	87
1a31a9ab	88	//Histogrammes setters and getters
1a31a9ab	89
b5dbb99b	90	void SetHistoPtSumRangeAndNBins(Float_t min, Float_t max, Int_t n) {
b5dbb99b	91	fHistoNPtSumBins = n ; fHistoPtSumMax = max ; fHistoPtSumMin = min ; }
1a31a9ab	92
b5dbb99b	93	Int_t GetHistoNPtSumBins() const { return fHistoNPtSumBins ; }
	94	Float_t GetHistoPtSumMin() const { return fHistoPtSumMin ; }
	95	Float_t GetHistoPtSumMax() const { return fHistoPtSumMax ; }
1a31a9ab	96
b5dbb99b	97	void SetHistoPtInConeRangeAndNBins(Float_t min, Float_t max, Int_t n) {
b5dbb99b	98	fHistoNPtInConeBins = n ; fHistoPtInConeMax = max ; fHistoPtInConeMin = min ; }
1a31a9ab	99
b5dbb99b	100	Int_t GetHistoNPtInConeBins() const { return fHistoNPtInConeBins; }
	101	Float_t GetHistoPtInConeMin() const { return fHistoPtInConeMin ; }
	102	Float_t GetHistoPtInConeMax() const { return fHistoPtInConeMax ; }
1a31a9ab	103
	104	private:
	105
	106	TString fCalorimeter ; // Calorimeter where neutral particles in cone for isolation are;
	107	Bool_t fReMakeIC ; // Do isolation analysis
	108	Bool_t fMakeSeveralIC ; // Do analysis for different IC
09273901	109	Bool_t fFillTMHisto; // Fill track matching plots
	110	Bool_t fFillSSHisto; // Fill Shower shape plots
	111
803d06a8	112	// Analysis data members for multiple cones and pt thresholds
	113	Int_t fNCones ; //! Number of cone sizes to test
	114	Int_t fNPtThresFrac ; //! Number of ptThres and ptFrac to test
	115
	116	Float_t fConeSizes[5] ; //! Array with cones to test
	117	Float_t fPtThresholds[5] ; //! Array with pt thresholds to test
	118	Float_t fPtFractions[5] ; //! Array with pt thresholds to test
	119
	120	TH1F* fhPtThresIsolated[5][5] ; //! Isolated particle with pt threshold
	121	TH1F* fhPtFracIsolated[5][5] ; //! Isolated particle with pt threshold
	122	TH2F* fhPtSumIsolated[5] ; //! Isolated particle with threshold on cone pt sum
1a31a9ab	123
	124	//Histograms
	125
0fb69ade	126	TH1F * fhEIso ; //! Number of isolated particles
1a31a9ab	127	TH1F * fhPtIso ; //! Number of isolated particles
	128	TH2F * fhPhiIso ; //! Phi of isolated particles
	129	TH2F * fhEtaIso ; //! eta of isolated particles
0fb69ade	130	TH2F * fhEtaPhiIso ; //! eta vs phi of isolated particles
d0a4f937	131	TH2F * fhEtaPhiNoIso ; //! eta vs phi of not isolated leading particles
1a31a9ab	132	TH1F * fhPtNoIso ; //! Number of not isolated leading particles
803d06a8	133	TH1F * fhPtDecayIso ; //! Number of isolated Pi0 decay particles (invariant mass tag)
803d06a8	134	TH1F * fhPtDecayNoIso ; //! Number of not isolated Pi0 decay leading particles (invariant mass tag)
d0a4f937	135	TH2F * fhEtaPhiDecayIso ; //! eta vs phi of isolated Pi0 decay particles
d0a4f937	136	TH2F * fhEtaPhiDecayNoIso ; //! eta vs phi of not isolated leading Pi0 decay particles
1a31a9ab	137	TH2F * fhConeSumPt ; //! Sum Pt in the cone
	138	TH2F * fhPtInCone ; //! Particle Pt in the cone
	139	TH2F * fhFRConeSumPt ; //! Sum Pt in the forward region cone (phi +90)
	140	TH2F * fhPtInFRCone ; //! Particle Pt in the forward region cone (phi +90 )
	141
1a31a9ab	142
	143	//MC
	144	TH1F * fhPtIsoPrompt; //! Number of isolated prompt gamma
	145	TH2F * fhPhiIsoPrompt; //! Phi of isolated prompt gamma
	146	TH2F * fhEtaIsoPrompt; //! eta of isolated prompt gamma
	147	TH1F * fhPtThresIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt threshold
	148	TH1F * fhPtFracIsolatedPrompt[5][5]; //! Isolated prompt gamma with pt frac
	149	TH2F * fhPtSumIsolatedPrompt[5]; //! Isolated prompt gamma with threshold on cone pt sume
	150	TH1F * fhPtIsoFragmentation; //! Number of isolated fragmentation gamma
	151	TH2F * fhPhiIsoFragmentation; //! Phi of isolated fragmentation gamma
	152	TH2F * fhEtaIsoFragmentation; //! eta of isolated fragmentation gamma
	153	TH1F * fhPtThresIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt threshold
	154	TH1F * fhPtFracIsolatedFragmentation[5][5]; //! Isolated fragmentation gamma with pt frac
	155	TH2F * fhPtSumIsolatedFragmentation[5]; //! Isolated fragmentation gamma with threshold on cone pt sume
803d06a8	156	TH1F * fhPtIsoPi0Decay; //! Number of isolated pi0 decay gamma
	157	TH2F * fhPhiIsoPi0Decay; //! Phi of isolated pi0 decay gamma
	158	TH2F * fhEtaIsoPi0Decay; //! eta of isolated pi0 decay gamma
	159	TH1F * fhPtThresIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt threshold
	160	TH1F * fhPtFracIsolatedPi0Decay[5][5]; //! Isolated pi0 decay gamma with pt frac
	161	TH2F * fhPtSumIsolatedPi0Decay[5]; //! Isolated pi0 decay gamma with threshold on cone pt sume
	162	TH1F * fhPtIsoEtaDecay; //! Number of isolated eta decay gamma
	163	TH2F * fhPhiIsoEtaDecay; //! Phi of isolated eta decay gamma
	164	TH2F * fhEtaIsoEtaDecay; //! eta of isolated eta decay gamma
	165	TH1F * fhPtThresIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt threshold
	166	TH1F * fhPtFracIsolatedEtaDecay[5][5]; //! Isolated eta decay gamma with pt frac
	167	TH2F * fhPtSumIsolatedEtaDecay[5]; //! Isolated eta fecay gamma with threshold on cone pt sume
	168	TH1F * fhPtIsoOtherDecay; //! Number of isolated other decay gamma
	169	TH2F * fhPhiIsoOtherDecay; //! Phi of isolated other decay gamma
	170	TH2F * fhEtaIsoOtherDecay; //! eta of isolated other decay gamma
1a31a9ab	171	TH1F * fhPtThresIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt threshold
	172	TH1F * fhPtFracIsolatedOtherDecay[5][5]; //! Isolated OtherDecay gamma with pt frac
	173	TH2F * fhPtSumIsolatedOtherDecay[5]; //! Isolated OtherDecay gamma with threshold on cone pt sume
	174	TH1F * fhPtIsoConversion; //! Number of isolated Conversion gamma
	175	TH2F * fhPhiIsoConversion; //! Phi of isolated Conversion gamma
	176	TH2F * fhEtaIsoConversion; //! eta of isolated Conversion gamma
	177	TH1F * fhPtThresIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt threshold
	178	TH1F * fhPtFracIsolatedConversion[5][5]; //! Isolated Conversion gamma with pt frac
	179	TH2F * fhPtSumIsolatedConversion[5]; //! Isolated Conversion gamma with threshold on cone pt sume
	180	TH1F * fhPtIsoUnknown; //! Number of isolated Unknown
	181	TH2F * fhPhiIsoUnknown; //! Phi of isolated Unknown
	182	TH2F * fhEtaIsoUnknown; //! eta of isolated Unknown
	183	TH1F * fhPtThresIsolatedUnknown[5][5]; //! Isolated Unknown gamma with pt threshold
	184	TH1F * fhPtFracIsolatedUnknown[5][5]; //! Isolated Unknown gamma with pt frac
	185	TH2F * fhPtSumIsolatedUnknown[5]; //! Isolated Unknown gamma with threshold on cone pt sume
	186
803d06a8	187	TH1F * fhPtNoIsoPi0Decay; //! Number of not isolated leading pi0 decay gamma
	188	TH1F * fhPtNoIsoEtaDecay; //! Number of not isolated leading eta decay gamma
	189	TH1F * fhPtNoIsoOtherDecay; //! Number of not isolated leading other decay gamma
1a31a9ab	190	TH1F * fhPtNoIsoPrompt; //! Number of not isolated leading prompt gamma
	191	TH1F * fhPtIsoMCPhoton; //! Number of isolated leading gamma
	192	TH1F * fhPtNoIsoMCPhoton; //! Number of not isolated leading gamma
0fb69ade	193	TH1F * fhPtNoIsoConversion; //! Number of not isolated leading conversion gamma
	194	TH1F * fhPtNoIsoFragmentation; //! Number of not isolated leading fragmentation gamma
	195	TH1F * fhPtNoIsoUnknown; //! Number of not isolated leading hadrons
1a31a9ab	196
09273901	197	// Track matching studies
	198	TH2F * fhTrackMatchedDEta ; //! Eta distance between track and cluster vs cluster E
	199	TH2F * fhTrackMatchedDPhi ; //! Phi distance between track and cluster vs cluster E
	200	TH2F * fhTrackMatchedDEtaDPhi ; //! Eta vs Phi distance between track and cluster, E cluster > 0.5 GeV
31ae6d59	201	TH2F * fhdEdx ; //! matched track dEdx vs cluster E
	202	TH2F * fhEOverP; //! matched track E cluster over P track vs cluster E, after dEdx cut
	203	TH2F * fhTrackMatchedMCParticle; //! Trace origin of matched particle
	204
09273901	205	// Shower Shape histograms
	206	TH2F * fhELambda0; //! Shower shape of isolated photons
	207	TH2F * fhELambda1; //! Shower shape of isolated photons
b5dbb99b	208	TH2F * fhELambda0TRD; //! Shower shape of isolated photons, SM behind TRD
	209	TH2F * fhELambda1TRD; //! Shower shape of isolated photons, SM behind TRD
	210
5c46c992	211	// Local maxima
	212	TH2F * fhNLocMax; //! number of maxima in selected clusters
	213	TH2F * fhELambda0LocMax1 ; //! E vs lambda0 of selected cluster, 1 local maxima in cluster
	214	TH2F * fhELambda1LocMax1 ; //! E vs lambda1 of selected cluster, 1 local maxima in cluster
	215	TH2F * fhELambda0LocMax2 ; //! E vs lambda0 of selected cluster, 2 local maxima in cluster
	216	TH2F * fhELambda1LocMax2 ; //! E vs lambda1 of selected cluster, 2 local maxima in cluster
	217	TH2F * fhELambda0LocMaxN ; //! E vs lambda0 of selected cluster, N>2 local maxima in cluster
	218	TH2F * fhELambda1LocMaxN ; //! E vs lambda1 of selected cluster, N>2 local maxima in cluster
	219
1a31a9ab	220	//Histograms settings
	221	Int_t fHistoNPtSumBins; // Number of bins in PtSum histograms
	222	Float_t fHistoPtSumMax; // PtSum maximum in histogram
	223	Float_t fHistoPtSumMin; // PtSum minimum in histogram
	224	Int_t fHistoNPtInConeBins; // Number of bins in PtInCone histogram
	225	Float_t fHistoPtInConeMax; // PtInCone maximum in histogram
	226	Float_t fHistoPtInConeMin; // PtInCone maximum in histogram
31ae6d59	227
09273901	228	AliAnaParticleIsolation( const AliAnaParticleIsolation & iso) ; // cpy ctor
09273901	229	AliAnaParticleIsolation & operator = (const AliAnaParticleIsolation & iso) ; // cpy assignment
c5693f62	230
d0a4f937	231	ClassDef(AliAnaParticleIsolation,12)
1a31a9ab	232	} ;
	233
	234
	235	#endif //ALIANAPARTICLEISOLATION_H
	236
	237
	238