]>
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 | |
0852acd0 | 26 | virtual ~AliAnaOmegaToPi0Gamma() ;//virtual dtor |
edff0964 | 27 | |
28 | private: | |
1e3689fa | 29 | AliAnaOmegaToPi0Gamma(const AliAnaOmegaToPi0Gamma & ex) ; // cpy ctor |
edff0964 | 30 | AliAnaOmegaToPi0Gamma & operator = (const AliAnaOmegaToPi0Gamma & ex) ;//cpy assignment |
31 | ||
32 | public: | |
0852acd0 | 33 | TList * GetCreateOutputObjects(); |
34 | void Print(const Option_t * opt) const; | |
35 | ||
36 | void InitParameters(); | |
37 | void MakeAnalysisFillHistograms(); | |
38 | void Terminate(TList * outList); | |
39 | ||
40 | TString GetInputAODPhotonName() const {return fInputAODGammaName;} | |
edff0964 | 41 | void SetInputAODPhotonName(TString name) { fInputAODGammaName = name; } |
0852acd0 | 42 | Bool_t IsBadRun(Int_t /*iRun*/) const {return kFALSE;} //Tests if this run bad according to private list |
43 | ||
4007d0a2 | 44 | void SetNEventsMixed(Int_t nevents) { fNmaxMixEv=nevents;} //events to be mixed |
d5b2a3f1 | 45 | |
46 | void SetNPID(Int_t pid) {fNpid=pid;} //N pid cut | |
47 | void SetNVtxZ(Int_t vtx){fNVtxZBin=vtx;} //N vertex Z cut | |
0852acd0 | 48 | void SetPi0MassPeakWidthCut(Double_t win){fPi0MassWindow=win;} |
d5b2a3f1 | 49 | |
4007d0a2 | 50 | void SetPi0OverOmegaPtCut(Double_t cut){fPi0OverOmegaPtCut=cut;} |
51 | void SetGammaOverOmegaPtCut(Double_t cut){fGammaOverOmegaPtCut=cut;} | |
0852acd0 | 52 | void ReadHistograms(TList * outputList); |
53 | ||
54 | private: | |
55 | ||
edff0964 | 56 | TClonesArray * fInputAODPi0; //Input AOD pi0 array |
57 | TString fInputAODGammaName; //Input AOD gamma name | |
58 | TList ** fEventsList; //event list for mixing | |
d5b2a3f1 | 59 | |
edff0964 | 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 | |
d5b2a3f1 | 66 | |
67 | Double_t *fVtxZCut; //![fNVtxZBin] vtertx z cut | |
68 | Double_t *fCent; //![fNCentBin] centrality cut | |
69 | Double_t *fRp; //![fNRpBin] reaction plane cut | |
d5b2a3f1 | 70 | |
edff0964 | 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 |