[u/mrichter/AliRoot.git] / PWG4 / PartCorr / 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")

class TH3D ;
class AliAODEvent ;
class AliESDEvent ;

#include "AliAnalysisTaskSE.h"
 
class AliAnaPi0 : public AliAnalysisTaskSE {
       
  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
              
    //obligatory methods
    virtual void UserCreateOutputObjects();
    virtual void Init() ;
    virtual void LocalInit() { Init() ; }
    virtual void UserExec(Option_t * opt = "") ;
    virtual void Terminate(Option_t * opt = "") ;
 
    void InitParameters();
    void Print(const Option_t * opt) const;

    void SetBadRunsList(){} ;     //Set list of runs which can be used for this analysis
    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
    void SetPtMin(Float_t pt = 0.2) {fPtMin=pt ;} //Cut on minimal pt of photon

 private:
    void FillHistograms() ; //fill invariant mass distributions
    Bool_t FillFromAOD(AliAODEvent * aod) ;
    Bool_t FillFromESD(AliESDEvent * esd) ;

    Bool_t IsBadRun(Int_t /*iRun*/) {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

       Int_t fCurCentrBin ; //! bin for current event (centrality)
       Int_t fCurZvertBin ; //! bin for current event (z position)
       Int_t fCurRPBin ;    //! bin for current event (reaction plain)

       Float_t fPtMin ;     //Cut on minimum photon Pt
       Float_t fZvtxCut ;   //Cut on vertex position
       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
       Float_t fDispCut;    //Cut on dispersion, used in PID evaluation
       Float_t fTOFCut;     //Cut on TOF, used in PID evaluation
       Float_t fPhotPID;    //Baesian PID for photon

       Double_t fVert[3] ;   //Currecnt vertex position
       TList ** fEventsList ;        //! containers for photons in stored events
       TClonesArray * fCurrentEvent ; //!Container  with latest event

       TList * fOutputList ; //list of output objects

       //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

       //Calo Cells
       ClassDef(AliAnaPi0,2)
 } ;


#endif //ALIANAPI0_H
Commit	Line	Data
f3050824	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
	13	class TH3D ;
	14	class AliAODEvent ;
	15	class AliESDEvent ;
	16
	17	#include "AliAnalysisTaskSE.h"
	18
	19	class AliAnaPi0 : public AliAnalysisTaskSE {
	20
	21	public:
	22
	23	AliAnaPi0() ; // default ctor
	24	AliAnaPi0(const char *name) ; // default ctor
	25	AliAnaPi0(const AliAnaPi0 & g) ; // cpy ctor
	26	AliAnaPi0 & operator = (const AliAnaPi0 & api0) ;//cpy assignment
	27	virtual ~AliAnaPi0() ;//virtual dtor
	28
	29	//obligatory methods
	30	virtual void UserCreateOutputObjects();
	31	virtual void Init() ;
	32	virtual void LocalInit() { Init() ; }
	33	virtual void UserExec(Option_t * opt = "") ;
	34	virtual void Terminate(Option_t * opt = "") ;
	35
	36	void InitParameters();
	37	void Print(const Option_t * opt) const;
	38
	39	void SetBadRunsList(){} ; //Set list of runs which can be used for this analysis
	40	void SetEtalonHisto(TH3D * h);//Provide etalon of binning for histograms
	41
	42	//Setters for parameters of event buffers
	43	void SetNCentrBin(Int_t n=5){fNCentrBin=n ;} //number of bins in centrality
	44	void SetNZvertBin(Int_t n=5){fNZvertBin=n ;} //number of bins for vertex position
	45	void SetNRPBin(Int_t n=6) {fNrpBin=n ;} //number of bins in reaction plain
	46	void SetNMaxEvMix(Int_t n=20){fNmaxMixEv=n ;}//Maximal number of events for mixing
	47
	48	//Setters for event selection
	49	void SetZvertexCut(Float_t zcut=40.){fZvtxCut=zcut ;} //cut on vertex position
	50	void SetPtMin(Float_t pt = 0.2) {fPtMin=pt ;} //Cut on minimal pt of photon
	51
	52	private:
	53	void FillHistograms() ; //fill invariant mass distributions
	54	Bool_t FillFromAOD(AliAODEvent * aod) ;
	55	Bool_t FillFromESD(AliESDEvent * esd) ;
	56
	57	Bool_t IsBadRun(Int_t /iRun/) {return kFALSE;} //Tests if this run bad according to private list
	58
	59	private:
	60
	61	Int_t fNCentrBin ; //number of bins in event container for centrality
	62	Int_t fNZvertBin ; //number of bins in event container for vertex position
	63	Int_t fNrpBin ; //number of bins in event container for reaction plain
	64	Int_t fNPID ; //Number of possible PID combinations
65	Int_t fNmaxMixEv ; //Maximal number of events stored in buffer for mixing
66
67	Int_t fCurCentrBin ; //! bin for current event (centrality)
68	Int_t fCurZvertBin ; //! bin for current event (z position)
69	Int_t fCurRPBin ; //! bin for current event (reaction plain)
70
71	Float_t fPtMin ; //Cut on minimum photon Pt
72	Float_t fZvtxCut ; //Cut on vertex position
73	Float_t fMinDist ; //Minimal distance to bad channel to accept cluster
74	Float_t fMinDist2; //Cuts on Minimal distance to study acceptance evaluation
75	Float_t fMinDist3; //One more cut on distance used for acceptance-efficiency study
76	Float_t fDispCut; //Cut on dispersion, used in PID evaluation
77	Float_t fTOFCut; //Cut on TOF, used in PID evaluation
78	Float_t fPhotPID; //Baesian PID for photon
79
80	Double_t fVert[3] ; //Currecnt vertex position
81	TList ** fEventsList ; //! containers for photons in stored events
82	TClonesArray * fCurrentEvent ; //!Container with latest event
83
84	TList * fOutputList ; //list of output objects
85
86	//Histograms
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	//Calo Cells
98	ClassDef(AliAnaPi0,2)
99	} ;
100
101
102	#endif //ALIANAPI0_H
103
104
105