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

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

//_________________________________________________________________________
// Class to fill two-photon invariant mass hisograms 
// to be used to extract pi0 raw yield.
//
//-- Author: Dmitri Peressounko (RRC "KI")
//-- Adapted to PartCorr frame by Lamia Benhabib (SUBATECH)
//-- and Gustavo Conesa (INFN-Frascati)

//Root
class TList;
class TH3D ;

//Analysis
class AliAODEvent ;
class AliESDEvent ;
#include "AliAnaPartCorrBaseClass.h"

#ifdef __PHOSUTIL__
  class AliPHOSGeoUtils;
#endif

#ifdef __EMCALUTIL__
  class AliEMCALGeoUtils;
#endif

class AliAnaPi0 : public AliAnaPartCorrBaseClass {
  
  public: 
  
  AliAnaPi0() ; // default ctor
  AliAnaPi0(const char *name) ; // default ctor
  AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
  AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
  virtual ~AliAnaPi0() ;//virtual dtor
  
  TList * GetCreateOutputObjects(); 
  
  void Print(const Option_t * opt) const;
  
  void Init();
  void InitParameters();
  
  //void MakeAnalysisFillAOD() {;} //Not needed
  void MakeAnalysisFillHistograms();
  
  //    void SetBadRunsList(){;} ;     //Set list of runs which can be used for this analysis
  //To be defined in future.
  void SetEtalonHisto(TH3D * h);//Provide etalon of binning for histograms
  
  //Setters for parameters of event buffers
  void SetNCentrBin(Int_t n=5){fNCentrBin=n ;} //number of bins in centrality 
  void SetNZvertBin(Int_t n=5){fNZvertBin=n ;} //number of bins for vertex position
  void SetNRPBin(Int_t n=6)   {fNrpBin=n ;}    //number of bins in reaction plain
  void SetNMaxEvMix(Int_t n=20){fNmaxMixEv=n ;}//Maximal number of events for mixing
  
  //Setters for event selection
  void SetZvertexCut(Float_t zcut=40.){fZvtxCut=zcut ;} //cut on vertex position
  
  TString GetCalorimeter()   const {return fCalorimeter ; }
  void SetCalorimeter(TString det)    {fCalorimeter = det ; }
  
  void Terminate(TList* outputList);
  void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
	
  private:
  Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list
  
  private:
  Int_t fNCentrBin ;		//number of bins in event container for centrality
  Int_t fNZvertBin ;		//number of bins in event container for vertex position
  Int_t fNrpBin ;			//number of bins in event container for reaction plain
  Int_t fNPID ;			//Number of possible PID combinations
  Int_t fNmaxMixEv ;		//Maximal number of events stored in buffer for mixing
  Float_t fZvtxCut ;		//Cut on vertex position
  TString fCalorimeter ;     //Select Calorimeter for IM
  
  TList ** fEventsList ;	//! containers for photons in stored events
  
  //Histograms
  
  TH3D * fhEtalon ; //Etalon histo, all distributions will have same binning as this one
  
  TH3D ** fhRe1 ;  //!REAL two-photon invariant mass distribution for different centralities and PID 
  TH3D ** fhMi1 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  TH3D ** fhRe2 ;  //!REAL two-photon invariant mass distribution for different centralities and PID 
  TH3D ** fhMi2 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  TH3D ** fhRe3 ;  //!REAL two-photon invariant mass distribution for different centralities and PID 
  TH3D ** fhMi3 ;  //!MIXED two-photon invariant mass distribution for different centralities and PID
  TH3D * fhEvents;  //!Number of events per centrality, RP, zbin
  
  //Acceptance
  TH1D * fhPrimPt ;    //! Spectrum of Primary 
  TH1D * fhPrimAccPt ; //! Spectrum of primary with accepted daughters 
  TH1D * fhPrimY ;     //! Rapidity distribution of primary particles
  TH1D * fhPrimAccY ;  //! Rapidity distribution of primary with accepted daughters
  TH1D * fhPrimPhi ;   //! Azimutal distribution of primary particles
  TH1D * fhPrimAccPhi; //! Azimutal distribution of primary with accepted daughters	
  
#ifdef __PHOSUTIL__
  AliPHOSGeoUtils * fPHOSGeo ; //! PHOS geometry pointer
#endif	
  
#ifdef __EMCALUTIL__
  AliEMCALGeoUtils * fEMCALGeo ; //! EMCAL geometry pointer
#endif	

  ClassDef(AliAnaPi0,4)
} ;


#endif //ALIANAPI0_H
Commit	Line	Data
1c5acb87	1	#ifndef ALIANAPI0_H
	2	#define ALIANAPI0_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5	/* $Id: $ */
	6
	7	//_________________________________________________________________________
	8	// Class to fill two-photon invariant mass hisograms
	9	// to be used to extract pi0 raw yield.
	10	//
	11	//-- Author: Dmitri Peressounko (RRC "KI")
	12	//-- Adapted to PartCorr frame by Lamia Benhabib (SUBATECH)
	13	//-- and Gustavo Conesa (INFN-Frascati)
	14
	15	//Root
	16	class TList;
	17	class TH3D ;
	18
	19	//Analysis
	20	class AliAODEvent ;
	21	class AliESDEvent ;
	22	#include "AliAnaPartCorrBaseClass.h"
	23
591cc579	24	#ifdef __PHOSUTIL__
6639984f	25	class AliPHOSGeoUtils;
1c5acb87	26	#endif
1c5acb87	27
591cc579	28	#ifdef __EMCALUTIL__
	29	class AliEMCALGeoUtils;
	30	#endif
	31
1c5acb87	32	class AliAnaPi0 : public AliAnaPartCorrBaseClass {
6639984f	33
	34	public:
	35
	36	AliAnaPi0() ; // default ctor
	37	AliAnaPi0(const char *name) ; // default ctor
	38	AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
	39	AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
	40	virtual ~AliAnaPi0() ;//virtual dtor
	41
	42	TList * GetCreateOutputObjects();
	43
	44	void Print(const Option_t * opt) const;
	45
	46	void Init();
	47	void InitParameters();
	48
	49	//void MakeAnalysisFillAOD() {;} //Not needed
	50	void MakeAnalysisFillHistograms();
	51
	52	// void SetBadRunsList(){;} ; //Set list of runs which can be used for this analysis
	53	//To be defined in future.
	54	void SetEtalonHisto(TH3D * h);//Provide etalon of binning for histograms
	55
	56	//Setters for parameters of event buffers
	57	void SetNCentrBin(Int_t n=5){fNCentrBin=n ;} //number of bins in centrality
	58	void SetNZvertBin(Int_t n=5){fNZvertBin=n ;} //number of bins for vertex position
	59	void SetNRPBin(Int_t n=6) {fNrpBin=n ;} //number of bins in reaction plain
	60	void SetNMaxEvMix(Int_t n=20){fNmaxMixEv=n ;}//Maximal number of events for mixing
	61
	62	//Setters for event selection
	63	void SetZvertexCut(Float_t zcut=40.){fZvtxCut=zcut ;} //cut on vertex position
	64
	65	TString GetCalorimeter() const {return fCalorimeter ; }
	66	void SetCalorimeter(TString det) {fCalorimeter = det ; }
	67
a5cc4f03	68	void Terminate(TList* outputList);
	69	void ReadHistograms(TList * outputList); //Fill histograms with histograms in ouput list, needed in Terminate.
	70
6639984f	71	private:
	72	Bool_t IsBadRun(Int_t /iRun/) const {return kFALSE;} //Tests if this run bad according to private list
	73
	74	private:
	75	Int_t fNCentrBin ; //number of bins in event container for centrality
	76	Int_t fNZvertBin ; //number of bins in event container for vertex position
	77	Int_t fNrpBin ; //number of bins in event container for reaction plain
	78	Int_t fNPID ; //Number of possible PID combinations
	79	Int_t fNmaxMixEv ; //Maximal number of events stored in buffer for mixing
	80	Float_t fZvtxCut ; //Cut on vertex position
	81	TString fCalorimeter ; //Select Calorimeter for IM
	82
	83	TList ** fEventsList ; //! containers for photons in stored events
	84
	85	//Histograms
	86
	87	TH3D * fhEtalon ; //Etalon histo, all distributions will have same binning as this one
	88
	89	TH3D ** fhRe1 ; //!REAL two-photon invariant mass distribution for different centralities and PID
	90	TH3D ** fhMi1 ; //!MIXED two-photon invariant mass distribution for different centralities and PID
	91	TH3D ** fhRe2 ; //!REAL two-photon invariant mass distribution for different centralities and PID
	92	TH3D ** fhMi2 ; //!MIXED two-photon invariant mass distribution for different centralities and PID
	93	TH3D ** fhRe3 ; //!REAL two-photon invariant mass distribution for different centralities and PID
	94	TH3D ** fhMi3 ; //!MIXED two-photon invariant mass distribution for different centralities and PID
	95	TH3D * fhEvents; //!Number of events per centrality, RP, zbin
	96
	97	//Acceptance
	98	TH1D * fhPrimPt ; //! Spectrum of Primary
	99	TH1D * fhPrimAccPt ; //! Spectrum of primary with accepted daughters
	100	TH1D * fhPrimY ; //! Rapidity distribution of primary particles
	101	TH1D * fhPrimAccY ; //! Rapidity distribution of primary with accepted daughters
	102	TH1D * fhPrimPhi ; //! Azimutal distribution of primary particles
	103	TH1D * fhPrimAccPhi; //! Azimutal distribution of primary with accepted daughters
	104
591cc579	105	#ifdef __PHOSUTIL__
6639984f	106	AliPHOSGeoUtils * fPHOSGeo ; //! PHOS geometry pointer
1c5acb87	107	#endif
6639984f	108
591cc579	109	#ifdef __EMCALUTIL__
	110	AliEMCALGeoUtils * fEMCALGeo ; //! EMCAL geometry pointer
	111	#endif
	112
	113	ClassDef(AliAnaPi0,4)
1c5acb87	114	} ;
	115
	116
	117	#endif //ALIANAPI0_H
	118
	119
	120