AliPhysics  vAN-20150822 (d56cf94)
 All Classes Namespaces Files Functions Variables Enumerations Enumerator Macros
AliAnaOmegaToPi0Gamma.h
Go to the documentation of this file.
1 #ifndef ALIANAOMEGATOPI0GAMMA_H
2 #define ALIANAOMEGATOPI0GAMMA_H
3 
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  * See cxx source for full Copyright notice */
6 
7 //_________________________________________________________________________
17 //_________________________________________________________________________
18 
19 //Root
20 class TList;
21 class TH2F ;
22 class TLorentzVector;
23 
24 //Analysis
26 class TParticle;
27 
29 
30  public:
31 
33 
34  AliAnaOmegaToPi0Gamma(const char *name) ;
35 
36  virtual ~AliAnaOmegaToPi0Gamma() ;
37 
38  TList * GetCreateOutputObjects();
39 
40  void Print(const Option_t * opt) const;
41 
42  void InitParameters();
43 
45 
46  void Terminate(TList * outList);
47 
48  TString GetInputAODPhotonName() const { return fInputAODGammaName ; }
49 
50  void SetInputAODPhotonName(TString & name) { fInputAODGammaName = name ; }
51 
53  Bool_t IsBadRun(Int_t /*iRun*/) const { return kFALSE;}
54 
55  void SetNCentBin(Int_t nbin) { fNCentBin = nbin ; }
56 
57  void SetNPID(Int_t pid) { fNpid = pid ; }
58 
59  void SetNVtxZ(Int_t vtx) { fNVtxZBin = vtx ; }
60 
61  void SetNDistToBadChannel(Int_t ndist) { fNBadChDistBin = ndist; }
62 
63  void SetPi0MassPeakWidthCut(Double_t win) { fPi0MassWindow = win ; }
64 
65  void SetPi0OverOmegaPtCut(Double_t cut) { fPi0OverOmegaPtCut = cut ; }
66 
67  void SetGammaOverOmegaPtCut(Double_t cut) { fGammaOverOmegaPtCut = cut ; }
68 
69  void SetEOverlapCluster(Double_t e) { fEOverlapCluster = e ; }
70 
71  void ReadHistograms(TList * outputList);
72 
73  private:
74 
75  TClonesArray * fInputAODPi0;
77  TList ** fEventsList;
78 
79  Int_t fNVtxZBin;
80  Int_t fNCentBin;
81  Int_t fNRpBin;
83  Int_t fNpid;
84 
86  Double_t * fVtxZCut;
87 
89  Double_t * fCent;
90 
92  Double_t * fRp;
93 
94  Double_t fPi0Mass;
95  Double_t fPi0MassWindow;
96  Double_t fPi0OverOmegaPtCut;
98  Double_t fEOverlapCluster;
99 
100  TH2F * fhEtalon;
101  TH2F **fRealOmega0;
102  TH2F **fMixAOmega0;
103  TH2F **fMixBOmega0;
104  TH2F **fMixCOmega0;
105  TH2F **fRealOmega1;
106  TH2F **fMixAOmega1;
107  TH2F **fMixBOmega1;
108  TH2F **fMixCOmega1;
109  TH2F **fRealOmega2;
110  TH2F **fMixAOmega2;
111  TH2F **fMixBOmega2;
112  TH2F **fMixCOmega2;
113 
114  TH2F **fhFakeOmega;
115 
116  TH1F * fhOmegaPriPt;
117 
120 
123 
125  ClassDef(AliAnaOmegaToPi0Gamma,2) ;
127 
128 } ;
129 
130 #endif //ALIANAOMEGATOPI0GAMMA_H
131 
132 
133 
TH2F ** fRealOmega1
! Real omega IVM(asy, pt, m), with Asy_pi0<0.7
TH2F ** fMixCOmega0
! MixC omega IVM(asy, pt, m)
TList * GetCreateOutputObjects()
Create histograms to be saved in output file.
Int_t fNRpBin
Number of reaction plane cut.
TH2F ** fhFakeOmega
! High pt clusters assumed as pi0 + another gamma
void Terminate(TList *outList)
Do some calculations and plots from the final histograms.
Double_t fGammaOverOmegaPtCut
gamma pt over omega pt cut
TString fInputAODGammaName
Input AOD gamma name.
TH2F ** fMixCOmega2
! MixC omega IVM(asy, pt, m)
AliAnaOmegaToPi0Gamma & operator=(const AliAnaOmegaToPi0Gamma &ex)
Assignment operator not implemented.
Int_t fNBadChDistBin
Number of bad channel dist cut.
void SetNCentBin(Int_t nbin)
TH2F ** fRealOmega0
! Real omega IVM(asy, pt, m), with Asy_pi0<1
TH2F ** fMixBOmega1
! MixB omega IVM(asy, pt, m)
TH2F ** fMixCOmega1
! MixC omega IVM(asy, pt, m)
omega(782)->pi0+gamma->3gamma
void SetInputAODPhotonName(TString &name)
Base class for CaloTrackCorr analysis algorithms.
Int_t fNVtxZBin
Number of vertex z cut.
Double_t * fCent
[fNVtxZBin]
virtual ~AliAnaOmegaToPi0Gamma()
Destructor.
TH2F ** fMixAOmega1
! MixA omega IVM(asy, pt, m)
AliAnaOmegaToPi0Gamma()
Default constructor.
void SetNDistToBadChannel(Int_t ndist)
Bool_t IsBadRun(Int_t) const
Tests if this run bad according to private list.
void SetPi0OverOmegaPtCut(Double_t cut)
void SetPi0MassPeakWidthCut(Double_t win)
TH2F ** fMixBOmega0
! MixB 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)
void SetGammaOverOmegaPtCut(Double_t cut)
Double_t fPi0MassWindow
pi0 mass windows
TH2F ** fMixBOmega2
! MixB omega IVM(asy, pt, m)
void MakeAnalysisFillHistograms()
Do analysis, fill histograms.
Int_t fNCentBin
Number of centrality cut.
TString GetInputAODPhotonName() const
TClonesArray * fInputAODPi0
Input AOD pi0 array.
void Print(const Option_t *opt) const
Print some relevant parameters set in the analysis.
TH1F * fhOmegaPriPt
! MC primary omega pt in 2pi and |y|<0.5
Double_t fPi0Mass
[fNRpBin]
Double_t fPi0OverOmegaPtCut
pi0 Pt over omega pt cut
void SetEOverlapCluster(Double_t e)
TH2F * fhEtalon
! An etalon of 3D histograms
Double_t fEOverlapCluster
the pt when the two photons overlapped
TH2F ** fMixAOmega0
! MixA omega IVM(asy, pt, m)
Double_t * fVtxZCut
Vertex z cut.
Int_t fNpid
Number of PID cut.
Double_t * fRp
[fNCentBin]
TList ** fEventsList
event list for mixing
void ReadHistograms(TList *outputList)