1 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
2 * See cxx source for full Copyright notice */
5 //_________________________________________________________________________
6 // class to extract omega(782)->pi0+gamma->3gamma
8 //-- Author: Renzhuo Wan (IOPP-Wuhan, China)
9 //_________________________________________________________________________
10 #ifndef ALIANAOMEGATOPI0GAMMA_H
11 #define ALIANAOMEGATOPI0GAMMA_H
17 #include "AliAnaPartCorrBaseClass.h"
20 class AliAnaOmegaToPi0Gamma : public AliAnaPartCorrBaseClass {
24 AliAnaOmegaToPi0Gamma() ; // default ctor
25 AliAnaOmegaToPi0Gamma(const char *name) ; // default ctor
26 virtual ~AliAnaOmegaToPi0Gamma() ;//virtual dtor
29 AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & ex) ; // cpy ctor
30 AliAnaOmegaToPi0Gamma & operator = (const AliAnaOmegaToPi0Gamma & ex) ;//cpy assignment
33 TList * GetCreateOutputObjects();
34 void Print(const Option_t * opt) const;
36 void InitParameters();
37 void MakeAnalysisFillHistograms();
38 void Terminate(TList * outList);
40 TString GetInputAODPhotonName() const {return fInputAODGammaName;}
41 void SetInputAODPhotonName(TString & name) { fInputAODGammaName = name; }
42 Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list
44 void SetNEventsMixed(Int_t nevents) { fNmaxMixEv=nevents;} //events to be mixed
46 void SetNPID(Int_t pid) {fNpid=pid;} //N pid cut
47 void SetNVtxZ(Int_t vtx){fNVtxZBin=vtx;} //N vertex Z cut
48 void SetPi0MassPeakWidthCut(Double_t win){fPi0MassWindow=win;}
50 void SetPi0OverOmegaPtCut(Double_t cut){fPi0OverOmegaPtCut=cut;}
51 void SetGammaOverOmegaPtCut(Double_t cut){fGammaOverOmegaPtCut=cut;}
52 void ReadHistograms(TList * outputList);
56 TClonesArray * fInputAODPi0; //Input AOD pi0 array
57 TString fInputAODGammaName; //Input AOD gamma name
58 TList ** fEventsList; //event list for mixing
60 Int_t fNVtxZBin; //Number of vertex z cut
61 Int_t fNCentBin; //Number of centrality cut
62 Int_t fNRpBin; //Number of reaction plane cut
63 Int_t fNBadChDistBin; //Number of bad channel dist cut
64 Int_t fNpid; //Number of PID cut
65 Int_t fNmaxMixEv; //buffer size events to be mixed
67 Double_t *fVtxZCut; //![fNVtxZBin] vtertx z cut
68 Double_t *fCent; //![fNCentBin] centrality cut
69 Double_t *fRp; //![fNRpBin] reaction plane cut
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
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
91 ClassDef(AliAnaOmegaToPi0Gamma,2)