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

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

//_________________________________________________________________________
// class to extract omega(782)->pi0+gamma->3gamma
//
//-- Author: Renzhuo Wan (IOPP-Wuhan, China)
//_________________________________________________________________________
#ifndef ALIANAOMEGATOPI0GAMMA_H
#define ALIANAOMEGATOPI0GAMMA_H
//Root
class TList;
class TH2F ;
class TLorentzVector;
//Analysis
#include "AliAnaPartCorrBaseClass.h"
class TParticle;

class AliAnaOmegaToPi0Gamma : public AliAnaPartCorrBaseClass {
  
  public: 
  
  AliAnaOmegaToPi0Gamma() ; // default ctor
  AliAnaOmegaToPi0Gamma(const char *name) ; // default ctor
  AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & g) ; // cpy ctor
  AliAnaOmegaToPi0Gamma & operator = (const AliAnaOmegaToPi0Gamma & api0) ;//cpy assignment
  virtual ~AliAnaOmegaToPi0Gamma() ;//virtual dtor
 
  TList * GetCreateOutputObjects(); 
  void Print(const Option_t * opt) const;
  
  void InitParameters();
  void MakeAnalysisFillHistograms();
  void Terminate(TList * outList);

  TString GetInputAODPhotonName() const {return fInputAODGammaName;}
  void SetInputAODPhotonName(TString name) {fInputAODGammaName=name;}
  Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list

  void SetNEventsMixed(Int_t nevents) { fNmaxMixEv=nevents;} //events to be mixed

  void SetNPID(Int_t pid) {fNpid=pid;} //N pid cut 
  void SetNVtxZ(Int_t vtx){fNVtxZBin=vtx;} //N vertex Z cut
  void SetNBadChDist(Int_t dist){fNBadChDistBin=dist;} //N dist cut
  void SetPi0MassPeakWidthCut(Double_t win){fPi0MassWindow=win;} 

  void SetPi0OverOmegaPtCut(Double_t cut){fPi0OverOmegaPtCut=cut;}
  void SetGammaOverOmegaPtCut(Double_t cut){fGammaOverOmegaPtCut=cut;}
  void ReadHistograms(TList * outputList);

  private:

  TClonesArray * fInputAODGamma; //!Input AOD gamma array 
  TClonesArray * fInputAODPi0;   //!Input AOD pi0 array
  TString fInputAODGammaName;    //!Input AOD gamma name
  TList ** fEventsList;          //!event list for mixing 
 
  Int_t fNVtxZBin;               //!Number of vertex z cut
  Int_t fNCentBin;               //!Number of centrality cut
  Int_t fNRpBin;                 //!Number of reaction plane cut
  Int_t fNBadChDistBin;          //!Number of bad channel dist cut
  Int_t fNpid;                   //!Number of PID cut
  Int_t fNmaxMixEv;              //!buffer size events to be mixed

  Double_t *fVtxZCut;            //![fNVtxZBin] vtertx z cut
  Double_t *fCent;               //![fNCentBin] centrality cut
  Double_t *fRp;                 //![fNRpBin] reaction plane cut
  Int_t *fBadChDist;             //![fNBadChDistBin] bad channel dist
  
  Double_t fPi0Mass;             //!nominal pi0 mass
  Double_t fPi0MassWindow;       //!pi0 mass windows
  Double_t fPi0OverOmegaPtCut;   //!pi0 Pt over omega pt cut
  Double_t fGammaOverOmegaPtCut; //!gamma pt over omega pt cut

  TH2F * fhEtalon;               //!an etalon of 3D histograms
  TH2F **fRealOmega0;             //!real omega IVM(asy, pt, m), with Asy_pi0<1 
  TH2F **fMixAOmega0;             //!mixA omega IVM(asy, pt, m) 
  TH2F **fMixBOmega0;             //!mixB omega IVM(asy, pt, m) 
  TH2F **fMixCOmega0;             //!mixC omega IVM(asy, pt, m) 
  TH2F **fRealOmega1;            //!real omega IVM(asy, pt, m), with Asy_pi0<0.7
  TH2F **fMixAOmega1;            //!mixA omega IVM(asy, pt, m)
  TH2F **fMixBOmega1;            //!mixB omega IVM(asy, pt, m)
  TH2F **fMixCOmega1;            //!mixC omega IVM(asy, pt, m)
  TH2F **fRealOmega2;            //!real omega IVM(asy, pt, m), with Asy_pi0<0.8
  TH2F **fMixAOmega2;            //!mixA omega IVM(asy, pt, m)
  TH2F **fMixBOmega2;            //!mixB omega IVM(asy, pt, m)
  TH2F **fMixCOmega2;            //!mixC omega IVM(asy, pt, m)
  TH1F *fhOmegaPriPt;            //!MC primary omega pt in 2pi and |y|<0.5

  ClassDef(AliAnaOmegaToPi0Gamma,2)
} ;

#endif
Commit	Line	Data
0852acd0	1	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	2	* See cxx source for full Copyright notice */
	3	/* $Id: $ */
	4
	5	//_________________________________________________________________________
	6	// class to extract omega(782)->pi0+gamma->3gamma
	7	//
	8	//-- Author: Renzhuo Wan (IOPP-Wuhan, China)
	9	//_________________________________________________________________________
	10	#ifndef ALIANAOMEGATOPI0GAMMA_H
	11	#define ALIANAOMEGATOPI0GAMMA_H
	12	//Root
	13	class TList;
	14	class TH2F ;
	15	class TLorentzVector;
	16	//Analysis
	17	#include "AliAnaPartCorrBaseClass.h"
	18	class TParticle;
	19
	20	class AliAnaOmegaToPi0Gamma : public AliAnaPartCorrBaseClass {
	21
	22	public:
	23
	24	AliAnaOmegaToPi0Gamma() ; // default ctor
	25	AliAnaOmegaToPi0Gamma(const char *name) ; // default ctor
	26	AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & g) ; // cpy ctor
	27	AliAnaOmegaToPi0Gamma & operator = (const AliAnaOmegaToPi0Gamma & api0) ;//cpy assignment
	28	virtual ~AliAnaOmegaToPi0Gamma() ;//virtual dtor
	29
	30	TList * GetCreateOutputObjects();
	31	void Print(const Option_t * opt) const;
	32
	33	void InitParameters();
	34	void MakeAnalysisFillHistograms();
	35	void Terminate(TList * outList);
	36
	37	TString GetInputAODPhotonName() const {return fInputAODGammaName;}
	38	void SetInputAODPhotonName(TString name) {fInputAODGammaName=name;}
	39	Bool_t IsBadRun(Int_t /iRun/) const {return kFALSE;} //Tests if this run bad according to private list
	40
4007d0a2	41	void SetNEventsMixed(Int_t nevents) { fNmaxMixEv=nevents;} //events to be mixed
d5b2a3f1	42
	43	void SetNPID(Int_t pid) {fNpid=pid;} //N pid cut
	44	void SetNVtxZ(Int_t vtx){fNVtxZBin=vtx;} //N vertex Z cut
	45	void SetNBadChDist(Int_t dist){fNBadChDistBin=dist;} //N dist cut
0852acd0	46	void SetPi0MassPeakWidthCut(Double_t win){fPi0MassWindow=win;}
d5b2a3f1	47
4007d0a2	48	void SetPi0OverOmegaPtCut(Double_t cut){fPi0OverOmegaPtCut=cut;}
4007d0a2	49	void SetGammaOverOmegaPtCut(Double_t cut){fGammaOverOmegaPtCut=cut;}
0852acd0	50	void ReadHistograms(TList * outputList);
	51
	52	private:
	53
d5b2a3f1	54	TClonesArray * fInputAODGamma; //!Input AOD gamma array
	55	TClonesArray * fInputAODPi0; //!Input AOD pi0 array
	56	TString fInputAODGammaName; //!Input AOD gamma name
	57	TList ** fEventsList; //!event list for mixing
	58
	59	Int_t fNVtxZBin; //!Number of vertex z cut
	60	Int_t fNCentBin; //!Number of centrality cut
	61	Int_t fNRpBin; //!Number of reaction plane cut
	62	Int_t fNBadChDistBin; //!Number of bad channel dist cut
	63	Int_t fNpid; //!Number of PID cut
	64	Int_t fNmaxMixEv; //!buffer size events to be mixed
	65
	66	Double_t *fVtxZCut; //![fNVtxZBin] vtertx z cut
	67	Double_t *fCent; //![fNCentBin] centrality cut
	68	Double_t *fRp; //![fNRpBin] reaction plane cut
	69	Int_t *fBadChDist; //![fNBadChDistBin] bad channel dist
	70
	71	Double_t fPi0Mass; //!nominal pi0 mass
	72	Double_t fPi0MassWindow; //!pi0 mass windows
	73	Double_t fPi0OverOmegaPtCut; //!pi0 Pt over omega pt cut
	74	Double_t fGammaOverOmegaPtCut; //!gamma pt over omega pt cut
	75
	76	TH2F * fhEtalon; //!an etalon of 3D histograms
	77	TH2F **fRealOmega0; //!real omega IVM(asy, pt, m), with Asy_pi0<1
	78	TH2F **fMixAOmega0; //!mixA omega IVM(asy, pt, m)
	79	TH2F **fMixBOmega0; //!mixB omega IVM(asy, pt, m)
	80	TH2F **fMixCOmega0; //!mixC omega IVM(asy, pt, m)
	81	TH2F **fRealOmega1; //!real omega IVM(asy, pt, m), with Asy_pi0<0.7
	82	TH2F **fMixAOmega1; //!mixA omega IVM(asy, pt, m)
	83	TH2F **fMixBOmega1; //!mixB omega IVM(asy, pt, m)
	84	TH2F **fMixCOmega1; //!mixC omega IVM(asy, pt, m)
	85	TH2F **fRealOmega2; //!real omega IVM(asy, pt, m), with Asy_pi0<0.8
	86	TH2F **fMixAOmega2; //!mixA omega IVM(asy, pt, m)
	87	TH2F **fMixBOmega2; //!mixB omega IVM(asy, pt, m)
	88	TH2F **fMixCOmega2; //!mixC omega IVM(asy, pt, m)
	89	TH1F *fhOmegaPriPt; //!MC primary omega pt in 2pi and \|y\|<0.5
0852acd0	90
4007d0a2	91	ClassDef(AliAnaOmegaToPi0Gamma,2)
0852acd0	92	} ;
	93
	94	#endif
	95
	96
	97