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