[u/mrichter/AliRoot.git] / PWG4 / PartCorrDep / AliAnaPhoton.h

#ifndef ALIANAPHOTON_H
#define ALIANAPHOTON_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice     */
/* $Id: AliAnaPhoton.h 27413 2008-07-18 13:28:12Z gconesab $ */

//_________________________________________________________________________
//
// Class for the photon identification.
// Clusters from calorimeters are identified as photons
// and kept in the AOD. Few histograms produced.
// Produces input for other analysis classes like AliAnaPi0, 
// AliAnaParticleHadronCorrelation ... 
//

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

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

// --- ANALYSIS system ---
#include "AliAnaPartCorrBaseClass.h"

class AliAnaPhoton : public AliAnaPartCorrBaseClass {

 public: 
  AliAnaPhoton() ;              // default ctor
  virtual ~AliAnaPhoton() { ; } // virtual dtor
 private:
  AliAnaPhoton(const AliAnaPhoton & g) ;               // cpy ctor
  AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment

 public:
	
  //---------------------------------------
  // General analysis frame methods
  //---------------------------------------
  
  TObjString * GetAnalysisCuts();
  
  TList      * GetCreateOutputObjects();
  
  void         Init();

  void         InitParameters();

  void         MakeAnalysisFillAOD()  ;

  void         MakeAnalysisFillHistograms() ; 
  
  void         Print(const Option_t * opt)const;
    
  
  // Analysis methods
  
  Bool_t       ClusterSelected(AliVCluster* cl, TLorentzVector mom) ;
  
  void         FillAcceptanceHistograms();

  void         FillShowerShapeHistograms( AliVCluster* cluster, const Int_t mcTag) ;
  
  void         SwitchOnFillShowerShapeHistograms()    { fFillSSHistograms = kTRUE  ; }
  void         SwitchOffFillShowerShapeHistograms()   { fFillSSHistograms = kFALSE ; }  
  
  
  // Analysis parameters setters getters
  
  TString      GetCalorimeter()                 const { return fCalorimeter        ; }
  void         SetCalorimeter(TString  & det)         { fCalorimeter = det         ; }
    
  // ** Cluster selection methods **
  
  void         SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {
                fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; }

  void         SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min; 
                                                        fTimeCutMax = max          ; }
  Double_t     GetTimeCutMin()                  const { return fTimeCutMin         ; }
  Double_t     GetTimeCutMax()                  const { return fTimeCutMax         ; }	
	
  void         SetNCellCut(Int_t n)                   { fNCellsCut = n             ; }
  Double_t     GetNCellCut()                    const { return fNCellsCut          ; }
  
  Bool_t       IsTrackMatchRejectionOn()        const { return fRejectTrackMatch   ; }
  void         SwitchOnTrackMatchRejection()          { fRejectTrackMatch = kTRUE  ; }
  void         SwitchOffTrackMatchRejection()         { fRejectTrackMatch = kFALSE ; }  
  	  
  void         FillNOriginHistograms(Int_t n)         { fNOriginHistograms = n ; 
    if(n > 14) fNOriginHistograms = 14; }
  void         FillNPrimaryHistograms(Int_t n)        { fNPrimaryHistograms= n ;
    if(n > 7)  fNPrimaryHistograms = 7; }

  // For histograms in arrays, index in the array, corresponding to a particle
  enum mcTypes    { mcPhoton = 0,        mcPi0Decay = 1,       mcOtherDecay = 2,  
                    mcPi0 = 3,           mcEta = 4,            mcElectron = 5,       
                    mcConversion = 6,    mcOther = 7,          mcAntiNeutron = 8,    
                    mcAntiProton = 9,    mcPrompt = 10,        mcFragmentation = 11, 
                    mcISR = 12,          mcString = 13                               };  

  enum mcPTypes   { mcPPhoton = 0,       mcPPi0Decay = 1,       mcPOtherDecay = 2,  mcPOther = 3,
                    mcPPrompt = 4,       mcPFragmentation = 5,  mcPISR = 6           };  
  
  enum mcssTypes  { mcssPhoton = 0,      mcssOther = 1,       mcssPi0 = 2,         
                    mcssEta = 3,         mcssConversion = 4,  mcssElectron = 5       };  
  
  private:
 
  TString  fCalorimeter ;                // Calorimeter where the gamma is searched;
  Float_t  fMinDist ;                    // Minimal distance to bad channel to accept cluster
  Float_t  fMinDist2;                    // Cuts on Minimal distance to study acceptance evaluation
  Float_t  fMinDist3;                    // One more cut on distance used for acceptance-efficiency study
  Bool_t   fRejectTrackMatch ;           // If PID on, reject clusters which have an associated TPC track
  Double_t fTimeCutMin  ;                // Remove clusters/cells with time smaller than this value, in ns
  Double_t fTimeCutMax  ;                // Remove clusters/cells with time larger than this value, in ns
  Int_t    fNCellsCut ;                  // Accept for the analysis clusters with more than fNCellsCut cells
  Bool_t   fFillSSHistograms ;           // Fill shower shape histograms
  Int_t    fNOriginHistograms;           // Fill only NOriginHistograms of the 14 defined types
  Int_t    fNPrimaryHistograms;          // Fill only NPrimaryHistograms of the 7 defined types

  //Histograms 
  TH2F * fhNCellsE;                      //! number of cells in cluster vs E 
  TH2F * fhMaxCellDiffClusterE;          //! Fraction of energy carried by cell with maximum energy
  TH2F * fhTimeE;                        //! time of cluster vs E 

  TH1F * fhEPhoton    ;                  //! Number of identified photon vs energy
  TH1F * fhPtPhoton   ;                  //! Number of identified photon vs transerse momentum 
  TH2F * fhPhiPhoton  ;                  //! Azimuthal angle of identified  photon vs transerse momentum 
  TH2F * fhEtaPhoton  ;                  //! Pseudorapidity of identified  photon vs transerse momentum 
  TH2F * fhEtaPhiPhoton  ;               //! Pseudorapidity vs Phi of identified  photon for transerse momentum > 0.5
  TH2F * fhEtaPhi05Photon  ;             //! Pseudorapidity vs Phi of identified  photon for transerse momentum < 0.5
  
  //Shower shape
  
  TH2F * fhDispE;                         //! cluster dispersion vs E
  TH2F * fhLam0E;                         //! cluster lambda0 vs  E
  TH2F * fhLam1E;                         //! cluster lambda1 vs  E  

  TH2F * fhDispETRD;                      //! cluster dispersion vs E, SM covered by TRD
  TH2F * fhLam0ETRD;                      //! cluster lambda0 vs  E, SM covered by TRD
  TH2F * fhLam1ETRD;                      //! cluster lambda1 vs  E, SM covered by TRD 

  TH2F * fhNCellsLam0LowE;                //! number of cells in cluster vs lambda0
  TH2F * fhNCellsLam1LowE;                //! number of cells in cluster vs lambda1
  TH2F * fhNCellsDispLowE;                //! number of cells in cluster vs dispersion
  TH2F * fhNCellsLam0HighE;               //! number of cells in cluster vs lambda0, E>2
  TH2F * fhNCellsLam1HighE;               //! number of cells in cluster vs lambda1, E>2
  TH2F * fhNCellsDispHighE;               //! number of cells in cluster vs dispersion, E>2
  
  TH2F * fhEtaLam0LowE;                   //! cluster eta vs lambda0, E<2
  TH2F * fhPhiLam0LowE;                   //! cluster phi vs lambda0, E<2
  TH2F * fhEtaLam0HighE;                  //! cluster eta vs lambda0, E>2
  TH2F * fhPhiLam0HighE;                  //! cluster phi vs lambda0, E>2
  TH2F * fhLam0DispLowE;                  //! cluster lambda0 vs dispersion, E<2
  TH2F * fhLam0DispHighE;                 //! cluster lambda0 vs dispersion, E>2
  TH2F * fhLam1Lam0LowE;                  //! cluster lambda1 vs lambda0, E<2
  TH2F * fhLam1Lam0HighE;                 //! cluster lambda1 vs lambda0, E>2
  TH2F * fhDispLam1LowE;                  //! cluster disp vs lambda1, E<2
  TH2F * fhDispLam1HighE;                 //! cluster disp vs lambda1, E>2
    
  //Fill MC dependent histograms, Origin of this cluster is ...

  TH2F * fhMCDeltaE[14]  ;                      //! MC-Reco E distribution coming from MC particle     
  TH2F * fhMCDeltaPt[14] ;                      //! MC-Reco pT distribution coming from MC particle
  TH2F * fhMC2E[14]  ;                          //! E distribution, Reco vs MC coming from MC particle
  TH2F * fhMC2Pt[14] ;                          //! pT distribution, Reco vs MC coming from MC particle
  
  TH1F * fhMCE[14];                             //! Number of identified photon vs cluster energy coming from MC particle
  TH1F * fhMCPt[14];                            //! Number of identified photon vs cluster pT     coming from MC particle
  TH2F * fhMCPhi[14];                           //! Phi of identified photon coming from MC particle
  TH2F * fhMCEta[14];                           //! eta of identified photon coming from MC particle

  TH1F * fhEPrimMC[7];                          //! Number of generated photon vs energy
  TH1F * fhPtPrimMC[7];                         //! Number of generated photon vs pT   
  TH2F * fhPhiPrimMC[7];                        //! Phi of generted photon
  TH2F * fhYPrimMC[7];                          //! Rapidity of generated photon 
  
  TH1F * fhEPrimMCAcc[7];                       //! Number of generated photon vs energy, in calorimeter acceptance
  TH1F * fhPtPrimMCAcc[7];                      //! Number of generated photon vs pT, in calorimeter acceptance   
  TH2F * fhPhiPrimMCAcc[7];                     //! Phi of generted photon, in calorimeter acceptance
  TH2F * fhYPrimMCAcc[7];                       //! Rapidity of generated photon, in calorimeter acceptance   
  
  // Shower Shape MC

  TH2F * fhMCELambda0[6] ;                      //! E vs Lambda0     from MC particle
  TH2F * fhMCELambda1[6] ;                      //! E vs Lambda1     from MC particle
  TH2F * fhMCEDispersion[6] ;                   //! E vs Dispersion  from MC particle
  
  TH2F * fhMCPhotonELambda0NoOverlap ;          //! E vs Lambda0     from MC photons, no overlap
  TH2F * fhMCPhotonELambda0TwoOverlap ;         //! E vs Lambda0     from MC photons, 2 particles overlap
  TH2F * fhMCPhotonELambda0NOverlap ;           //! E vs Lambda0     from MC photons, N particles overlap
  
  TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6];  //! Lambda0 vs fraction of energy of max cell for E < 2 GeV
  TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6];  //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV
  TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6];  //! Lambda0 vs fraction of energy of max cell for E > 6 GeV
  TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6];   //! NCells  vs fraction of energy of max cell for E < 2
  TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6];   //! NCells  vs fraction of energy of max cell for 2 < E < 6 GeV
  TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6];   //! NCells  vs fraction of energy of max cell for E > 6
  TH2F * fhMCNCellsE[6];                        //! NCells per cluster vs energy
  TH2F * fhMCMaxCellDiffClusterE[6];            //! Fraction of energy carried by cell with maximum energy

  //Embedding
  TH2F * fhEmbeddedSignalFractionEnergy ;       //! Fraction of photon energy of embedded signal vs cluster energy
  
  TH2F * fhEmbedPhotonELambda0FullSignal ;      //!  Lambda0 vs E for embedded photons with more than 90% of the cluster energy
  TH2F * fhEmbedPhotonELambda0MostlySignal ;    //!  Lambda0 vs E for embedded photons with 90%<fraction<50% 
  TH2F * fhEmbedPhotonELambda0MostlyBkg ;       //!  Lambda0 vs E for embedded photons with 50%<fraction<10% 
  TH2F * fhEmbedPhotonELambda0FullBkg ;         //!  Lambda0 vs E for embedded photons with less than 10% of the cluster energy
  
  TH2F * fhEmbedPi0ELambda0FullSignal ;         //!  Lambda0 vs E for embedded photons with more than 90% of the cluster energy
  TH2F * fhEmbedPi0ELambda0MostlySignal ;       //!  Lambda0 vs E for embedded photons with 90%<fraction<50% 
  TH2F * fhEmbedPi0ELambda0MostlyBkg ;          //!  Lambda0 vs E for embedded photons with 50%<fraction<10% 
  TH2F * fhEmbedPi0ELambda0FullBkg ;            //!  Lambda0 vs E for embedded photons with less than 10% of the cluster energy
  
   ClassDef(AliAnaPhoton,19)

} ;
 
#endif//ALIANAPHOTON_H
Commit	Line	Data
1c5acb87	1	#ifndef ALIANAPHOTON_H
	2	#define ALIANAPHOTON_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5	/* $Id: AliAnaPhoton.h 27413 2008-07-18 13:28:12Z gconesab $ */
	6
	7	//_________________________________________________________________________
	8	//
	9	// Class for the photon identification.
	10	// Clusters from calorimeters are identified as photons
	11	// and kept in the AOD. Few histograms produced.
6175da48	12	// Produces input for other analysis classes like AliAnaPi0,
6175da48	13	// AliAnaParticleHadronCorrelation ...
1c5acb87	14	//
	15
	16	//-- Author: Gustavo Conesa (INFN-LNF)
	17
	18	// --- ROOT system ---
	19	class TH2F ;
123fc3bd	20	class TH1F;
1c5acb87	21	class TString ;
0c1383b5	22	class TObjString;
5812a064	23	class TList ;
1c5acb87	24
	25	// --- ANALYSIS system ---
	26	#include "AliAnaPartCorrBaseClass.h"
1c5acb87	27
	28	class AliAnaPhoton : public AliAnaPartCorrBaseClass {
	29
78219bac	30	public:
5812a064	31	AliAnaPhoton() ; // default ctor
5812a064	32	virtual ~AliAnaPhoton() { ; } // virtual dtor
78219bac	33	private:
5812a064	34	AliAnaPhoton(const AliAnaPhoton & g) ; // cpy ctor
5812a064	35	AliAnaPhoton & operator = (const AliAnaPhoton & g) ; // cpy assignment
78219bac	36
78219bac	37	public:
0c1383b5	38
6175da48	39	//---------------------------------------
	40	// General analysis frame methods
	41	//---------------------------------------
c4a7d28a	42
0c1383b5	43	TObjString * GetAnalysisCuts();
6175da48	44
0c1383b5	45	TList * GetCreateOutputObjects();
c4a7d28a	46
6175da48	47	void Init();
6639984f	48
6175da48	49	void InitParameters();
	50
	51	void MakeAnalysisFillAOD() ;
	52
	53	void MakeAnalysisFillHistograms() ;
1c5acb87	54
6175da48	55	void Print(const Option_t * opt)const;
521636d2	56
3d5d5078	57
	58	// Analysis methods
	59
521636d2	60	Bool_t ClusterSelected(AliVCluster* cl, TLorentzVector mom) ;
1c5acb87	61
3d5d5078	62	void FillAcceptanceHistograms();
3d5d5078	63
3d5d5078	64	void FillShowerShapeHistograms( AliVCluster* cluster, const Int_t mcTag) ;
	65
	66	void SwitchOnFillShowerShapeHistograms() { fFillSSHistograms = kTRUE ; }
	67	void SwitchOffFillShowerShapeHistograms() { fFillSSHistograms = kFALSE ; }
	68
	69
6175da48	70	// Analysis parameters setters getters
c4a7d28a	71
521636d2	72	TString GetCalorimeter() const { return fCalorimeter ; }
	73	void SetCalorimeter(TString & det) { fCalorimeter = det ; }
	74
6175da48	75	// Cluster selection methods
6175da48	76
c4a7d28a	77	void SetMinDistanceToBadChannel(Float_t m1, Float_t m2, Float_t m3) {
521636d2	78	fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3; }
6175da48	79
c4a7d28a	80	void SetTimeCut(Double_t min, Double_t max) { fTimeCutMin = min;
521636d2	81	fTimeCutMax = max ; }
	82	Double_t GetTimeCutMin() const { return fTimeCutMin ; }
	83	Double_t GetTimeCutMax() const { return fTimeCutMax ; }
1e86c71e	84
521636d2	85	void SetNCellCut(Int_t n) { fNCellsCut = n ; }
521636d2	86	Double_t GetNCellCut() const { return fNCellsCut ; }
c4a7d28a	87
	88	Bool_t IsTrackMatchRejectionOn() const { return fRejectTrackMatch ; }
	89	void SwitchOnTrackMatchRejection() { fRejectTrackMatch = kTRUE ; }
	90	void SwitchOffTrackMatchRejection() { fRejectTrackMatch = kFALSE ; }
5812a064	91
f66d95af	92	void FillNOriginHistograms(Int_t n) { fNOriginHistograms = n ;
	93	if(n > 14) fNOriginHistograms = 14; }
	94	void FillNPrimaryHistograms(Int_t n) { fNPrimaryHistograms= n ;
	95	if(n > 7) fNPrimaryHistograms = 7; }
	96
3d5d5078	97	// For histograms in arrays, index in the array, corresponding to a particle
f66d95af	98	enum mcTypes { mcPhoton = 0, mcPi0Decay = 1, mcOtherDecay = 2,
	99	mcPi0 = 3, mcEta = 4, mcElectron = 5,
	100	mcConversion = 6, mcOther = 7, mcAntiNeutron = 8,
	101	mcAntiProton = 9, mcPrompt = 10, mcFragmentation = 11,
	102	mcISR = 12, mcString = 13 };
41121cfe	103
f66d95af	104	enum mcPTypes { mcPPhoton = 0, mcPPi0Decay = 1, mcPOtherDecay = 2, mcPOther = 3,
	105	mcPPrompt = 4, mcPFragmentation = 5, mcPISR = 6 };
	106
3f5990d6	107	enum mcssTypes { mcssPhoton = 0, mcssOther = 1, mcssPi0 = 2,
f66d95af	108	mcssEta = 3, mcssConversion = 4, mcssElectron = 5 };
3d5d5078	109
1c5acb87	110	private:
1c5acb87	111
6175da48	112	TString fCalorimeter ; // Calorimeter where the gamma is searched;
	113	Float_t fMinDist ; // Minimal distance to bad channel to accept cluster
	114	Float_t fMinDist2; // Cuts on Minimal distance to study acceptance evaluation
	115	Float_t fMinDist3; // One more cut on distance used for acceptance-efficiency study
	116	Bool_t fRejectTrackMatch ; // If PID on, reject clusters which have an associated TPC track
6175da48	117	Double_t fTimeCutMin ; // Remove clusters/cells with time smaller than this value, in ns
	118	Double_t fTimeCutMax ; // Remove clusters/cells with time larger than this value, in ns
	119	Int_t fNCellsCut ; // Accept for the analysis clusters with more than fNCellsCut cells
c4a7d28a	120	Bool_t fFillSSHistograms ; // Fill shower shape histograms
f66d95af	121	Int_t fNOriginHistograms; // Fill only NOriginHistograms of the 14 defined types
f66d95af	122	Int_t fNPrimaryHistograms; // Fill only NPrimaryHistograms of the 7 defined types
521636d2	123
2244659d	124	//Histograms
c4a7d28a	125	TH2F * fhNCellsE; //! number of cells in cluster vs E
f66d95af	126	TH2F * fhMaxCellDiffClusterE; //! Fraction of energy carried by cell with maximum energy
f15c25da	127	TH2F * fhTimeE; //! time of cluster vs E
f15c25da	128
20218aea	129	TH1F * fhEPhoton ; //! Number of identified photon vs energy
6175da48	130	TH1F * fhPtPhoton ; //! Number of identified photon vs transerse momentum
	131	TH2F * fhPhiPhoton ; //! Azimuthal angle of identified photon vs transerse momentum
	132	TH2F * fhEtaPhoton ; //! Pseudorapidity of identified photon vs transerse momentum
	133	TH2F * fhEtaPhiPhoton ; //! Pseudorapidity vs Phi of identified photon for transerse momentum > 0.5
	134	TH2F * fhEtaPhi05Photon ; //! Pseudorapidity vs Phi of identified photon for transerse momentum < 0.5
123fc3bd	135
521636d2	136	//Shower shape
f66d95af	137
521636d2	138	TH2F * fhDispE; //! cluster dispersion vs E
	139	TH2F * fhLam0E; //! cluster lambda0 vs E
	140	TH2F * fhLam1E; //! cluster lambda1 vs E
7c65ad18	141
521636d2	142	TH2F * fhDispETRD; //! cluster dispersion vs E, SM covered by TRD
	143	TH2F * fhLam0ETRD; //! cluster lambda0 vs E, SM covered by TRD
	144	TH2F * fhLam1ETRD; //! cluster lambda1 vs E, SM covered by TRD
7c65ad18	145
521636d2	146	TH2F * fhNCellsLam0LowE; //! number of cells in cluster vs lambda0
	147	TH2F * fhNCellsLam1LowE; //! number of cells in cluster vs lambda1
	148	TH2F * fhNCellsDispLowE; //! number of cells in cluster vs dispersion
	149	TH2F * fhNCellsLam0HighE; //! number of cells in cluster vs lambda0, E>2
	150	TH2F * fhNCellsLam1HighE; //! number of cells in cluster vs lambda1, E>2
	151	TH2F * fhNCellsDispHighE; //! number of cells in cluster vs dispersion, E>2
	152
521636d2	153	TH2F * fhEtaLam0LowE; //! cluster eta vs lambda0, E<2
	154	TH2F * fhPhiLam0LowE; //! cluster phi vs lambda0, E<2
	155	TH2F * fhEtaLam0HighE; //! cluster eta vs lambda0, E>2
	156	TH2F * fhPhiLam0HighE; //! cluster phi vs lambda0, E>2
	157	TH2F * fhLam0DispLowE; //! cluster lambda0 vs dispersion, E<2
	158	TH2F * fhLam0DispHighE; //! cluster lambda0 vs dispersion, E>2
	159	TH2F * fhLam1Lam0LowE; //! cluster lambda1 vs lambda0, E<2
	160	TH2F * fhLam1Lam0HighE; //! cluster lambda1 vs lambda0, E>2
	161	TH2F * fhDispLam1LowE; //! cluster disp vs lambda1, E<2
	162	TH2F * fhDispLam1HighE; //! cluster disp vs lambda1, E>2
7c65ad18	163
4c8f7c2e	164	//Fill MC dependent histograms, Origin of this cluster is ...
4c8f7c2e	165
5812a064	166	TH2F * fhMCDeltaE[14] ; //! MC-Reco E distribution coming from MC particle
	167	TH2F * fhMCDeltaPt[14] ; //! MC-Reco pT distribution coming from MC particle
	168	TH2F * fhMC2E[14] ; //! E distribution, Reco vs MC coming from MC particle
	169	TH2F * fhMC2Pt[14] ; //! pT distribution, Reco vs MC coming from MC particle
4c8f7c2e	170
5812a064	171	TH1F * fhMCE[14]; //! Number of identified photon vs cluster energy coming from MC particle
	172	TH1F * fhMCPt[14]; //! Number of identified photon vs cluster pT coming from MC particle
	173	TH2F * fhMCPhi[14]; //! Phi of identified photon coming from MC particle
	174	TH2F * fhMCEta[14]; //! eta of identified photon coming from MC particle
3d5d5078	175
5812a064	176	TH1F * fhEPrimMC[7]; //! Number of generated photon vs energy
	177	TH1F * fhPtPrimMC[7]; //! Number of generated photon vs pT
	178	TH2F * fhPhiPrimMC[7]; //! Phi of generted photon
	179	TH2F * fhYPrimMC[7]; //! Rapidity of generated photon
3d5d5078	180
5812a064	181	TH1F * fhEPrimMCAcc[7]; //! Number of generated photon vs energy, in calorimeter acceptance
	182	TH1F * fhPtPrimMCAcc[7]; //! Number of generated photon vs pT, in calorimeter acceptance
	183	TH2F * fhPhiPrimMCAcc[7]; //! Phi of generted photon, in calorimeter acceptance
	184	TH2F * fhYPrimMCAcc[7]; //! Rapidity of generated photon, in calorimeter acceptance
f66d95af	185
521636d2	186	// Shower Shape MC
521636d2	187
5812a064	188	TH2F * fhMCELambda0[6] ; //! E vs Lambda0 from MC particle
	189	TH2F * fhMCELambda1[6] ; //! E vs Lambda1 from MC particle
	190	TH2F * fhMCEDispersion[6] ; //! E vs Dispersion from MC particle
f66d95af	191
5812a064	192	TH2F * fhMCPhotonELambda0NoOverlap ; //! E vs Lambda0 from MC photons, no overlap
	193	TH2F * fhMCPhotonELambda0TwoOverlap ; //! E vs Lambda0 from MC photons, 2 particles overlap
	194	TH2F * fhMCPhotonELambda0NOverlap ; //! E vs Lambda0 from MC photons, N particles overlap
f66d95af	195
	196	TH2F * fhMCLambda0vsClusterMaxCellDiffE0[6]; //! Lambda0 vs fraction of energy of max cell for E < 2 GeV
	197	TH2F * fhMCLambda0vsClusterMaxCellDiffE2[6]; //! Lambda0 vs fraction of energy of max cell for 2< E < 6 GeV
	198	TH2F * fhMCLambda0vsClusterMaxCellDiffE6[6]; //! Lambda0 vs fraction of energy of max cell for E > 6 GeV
	199	TH2F * fhMCNCellsvsClusterMaxCellDiffE0[6]; //! NCells vs fraction of energy of max cell for E < 2
	200	TH2F * fhMCNCellsvsClusterMaxCellDiffE2[6]; //! NCells vs fraction of energy of max cell for 2 < E < 6 GeV
	201	TH2F * fhMCNCellsvsClusterMaxCellDiffE6[6]; //! NCells vs fraction of energy of max cell for E > 6
	202	TH2F * fhMCNCellsE[6]; //! NCells per cluster vs energy
	203	TH2F * fhMCMaxCellDiffClusterE[6]; //! Fraction of energy carried by cell with maximum energy
	204
3d5d5078	205	//Embedding
5812a064	206	TH2F * fhEmbeddedSignalFractionEnergy ; //! Fraction of photon energy of embedded signal vs cluster energy
3d5d5078	207
5812a064	208	TH2F * fhEmbedPhotonELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
	209	TH2F * fhEmbedPhotonELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
	210	TH2F * fhEmbedPhotonELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
	211	TH2F * fhEmbedPhotonELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
3d5d5078	212
5812a064	213	TH2F * fhEmbedPi0ELambda0FullSignal ; //! Lambda0 vs E for embedded photons with more than 90% of the cluster energy
	214	TH2F * fhEmbedPi0ELambda0MostlySignal ; //! Lambda0 vs E for embedded photons with 90%<fraction<50%
	215	TH2F * fhEmbedPi0ELambda0MostlyBkg ; //! Lambda0 vs E for embedded photons with 50%<fraction<10%
	216	TH2F * fhEmbedPi0ELambda0FullBkg ; //! Lambda0 vs E for embedded photons with less than 10% of the cluster energy
3d5d5078	217
f15c25da	218	ClassDef(AliAnaPhoton,19)
6639984f	219
1c5acb87	220	} ;
1c5acb87	221
1c5acb87	222	#endif//ALIANAPHOTON_H
	223
	224
	225