1 #ifndef ALIANALYSISET_H
2 #define ALIANALYSISET_H
3 //_________________________________________________________________________
4 // Utility Class for transverse energy studies
5 // Base class for ESD & MC analysis
6 // - reconstruction and MonteCarlo output
8 //*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
9 //_________________________________________________________________________
18 class AliESDtrackCuts;
21 class AliAnalysisEtCuts;
28 virtual ~AliAnalysisEt();
32 /** Analyse the event! */
33 virtual Int_t AnalyseEvent(AliVEvent *event);
35 /** Fill the objects you want to output, classes which add new histograms should overload this. */
36 virtual void FillOutputList(TList* list);
38 /** Initialise the analysis, must be overloaded. */
42 * Creates the histograms, must be overloaded if you want to add your own.
43 * Uses the fHistogramNameSuffix to create proper histogram names
45 virtual void CreateHistograms();
46 virtual void CreateTree();
48 /** Fills the histograms, must be overloaded if you want to add your own */
49 virtual void FillHistograms();
51 /** Reset event specific values (Et etc.) */
52 virtual void ResetEventValues();
54 /** Set Particle codes/mass */
55 virtual void SetParticleCodes();
58 AliAnalysisEtCuts * GetCuts() const { return fCuts; }
59 virtual void SetCuts(const AliAnalysisEtCuts *cuts)
60 { fCuts = (AliAnalysisEtCuts *) cuts; }
62 /** Total Et in the event (without acceptance cuts) */
63 Double_t GetTotEt() const { return fTotEt; }
65 /** Total Et in the event within the acceptance cuts */
66 Double_t GetTotEtAcc() const { return fTotEtAcc; }
68 /** Total neutral Et in the event (without acceptance cuts) */
69 Double_t GetTotNeutralEt() const { return fTotNeutralEt; }
71 /** Total neutral Et in the event within the acceptance cuts */
72 Double_t GetTotNeutralEtAcc() const { return fTotNeutralEtAcc; }
74 /** Total charged Et in the event (without acceptance cuts) */
75 Double_t GetTotChargedEt() const { return fTotChargedEt; }
77 /** Total charged Et in the event within the acceptance cuts */
78 Double_t GetTotChargedEtAcc() const { return fTotChargedEtAcc; }
80 void SetTPCOnlyTrackCuts(const AliESDtrackCuts *cuts){ fEsdtrackCutsTPC = (AliESDtrackCuts *) cuts;}
84 TString fHistogramNameSuffix; /** The suffix for the histogram names */
86 AliAnalysisEtCuts *fCuts; // keeper of basic cuts
89 TDatabasePDG *fPdgDB;//data base used for looking up pdg codes
90 //these codes are stored as variables because otherwise there were issues using this with the plugin
91 Int_t fPiPlusCode;//pdg pi plus code
92 Int_t fPiMinusCode;//pdg pi minus code
93 Int_t fKPlusCode;// pdg k plus code
94 Int_t fKMinusCode;//pdg k minus code
95 Int_t fProtonCode;//pdg proton code
96 Int_t fAntiProtonCode;//pdg antiproton code
97 Int_t fLambdaCode;// pdg lambda code
98 Int_t fAntiLambdaCode;//pdg antilambda code
99 Int_t fK0SCode;//pdg k0 short code
100 Int_t fOmegaCode;//pdg omega code
101 Int_t fAntiOmegaCode;//pdg anti-omega code
102 Int_t fXi0Code;//pdg xi-0 code
103 Int_t fAntiXi0Code;//pdg anti-xi0 code
104 Int_t fXiCode;//pdg xi code
105 Int_t fAntiXiCode;//pdg anti-xi code
106 Int_t fSigmaCode;//pdg sigma code
107 Int_t fAntiSigmaCode;//pdg anti-sigma code
108 Int_t fK0LCode;//pdg k0 long code
109 Int_t fNeutronCode;//pdg neutron code
110 Int_t fAntiNeutronCode;//pdg anti-neutron code
111 Int_t fEPlusCode;//pdg positron code
112 Int_t fEMinusCode;//pdg electron code
113 Float_t fPionMass;//pdg pion mass
115 Double_t fTotEt;/** Total Et in the event (without acceptance cuts) */
116 Double_t fTotEtAcc;/** Total Et in the event within the acceptance cuts */
118 Double_t fTotNeutralEt;/** Total neutral Et in the event */
119 Double_t fTotNeutralEtAcc;/** Total neutral Et in the event within the acceptance cuts */
120 Double_t fTotChargedEt;/** Total charged Et in the event */
121 Double_t fTotChargedEtAcc;/** Total charged Et in the event within the acceptance cuts */
123 Int_t fMultiplicity;/** Multiplicity of particles in the event */
124 Int_t fChargedMultiplicity;/** Multiplicity of charged particles in the event */
125 Int_t fNeutralMultiplicity; /** Multiplicity of neutral particles in the event */
127 Double_t fBaryonEt; /** Et of identified baryons */
128 Double_t fAntiBaryonEt; /** Et of identified anti-baryons */
129 Double_t fMesonEt; /** Et of identified mesons */
131 Double_t fBaryonEtAcc; /** Et of identified baryons in calorimeter acceptance */
132 Double_t fAntiBaryonEtAcc; /** Et of identified anti-baryons in calorimeter acceptance */
133 Double_t fMesonEtAcc; /** Et of identified mesons in calorimeter acceptance */
135 Double_t fProtonEt; /** Et of identified protons */
136 Double_t fChargedKaonEt; /** Et of identified charged kaons */
137 Double_t fMuonEt; /** Et of identified muons */
138 Double_t fElectronEt; /** Et of identified electrons */
140 Double_t fProtonEtAcc; /** Et of identified protons in calorimeter acceptance */
141 Double_t fChargedKaonEtAcc; /** Et of identified charged kaons in calorimeter acceptance */
142 Double_t fMuonEtAcc; /** Et of identified muons in calorimeter acceptance */
143 Double_t fElectronEtAcc; /** Et of identified electrons in calorimeter acceptance */
145 Double_t fEtaCut;/** Cut in eta (standard |eta| < 0.5 )*/
147 Double_t fEtaCutAcc;/** Eta cut for our acceptance */
148 Double_t fPhiCutAccMin; /** Min phi cut for our acceptance in radians */
149 Double_t fPhiCutAccMax; /** Max phi cut for our acceptance in radians */
150 Double_t fDetectorRadius; /** Detector radius */
152 Double_t fClusterEnergyCut; /** Cut on the cluster energy */
153 Double_t fSingleCellEnergyCut; /** Minimum energy to cut on single cell cluster */
155 // Declare the histograms
157 /** The full Et spectrum measured */
158 TH1F *fHistEt; //Et spectrum
160 /** The full charged Et spectrum measured */
161 TH1F *fHistChargedEt; //Charged Et spectrum
163 /** The full neutral Et spectrum measured */
164 TH1F *fHistNeutralEt; //Neutral Et spectrum
166 /** The Et spectrum within the calorimeter acceptance */
167 TH1F *fHistEtAcc; //Et in acceptance
169 /** The charged Et spectrum within the calorimeter acceptance */
170 TH1F *fHistChargedEtAcc; //Charged Et in acceptance
172 /** The neutral Et spectrum within the calorimeter acceptance */
173 TH1F *fHistNeutralEtAcc; //Et in acceptance
175 /** Multiplicity of particles in the events */
176 TH1F *fHistMult; //Multiplicity
178 /** Charged multiplicity of particles in the events */
179 TH1F *fHistChargedMult; //Charged multiplicity
181 /** Neutral multiplicity of particles in the events */
182 TH1F *fHistNeutralMult; //Neutral multiplicity
184 /* Acceptance plots */
185 TH2F *fHistPhivsPtPos; //phi vs pT plot for positive tracks
186 TH2F *fHistPhivsPtNeg; //phi vs pT plot for negative tracks
189 TH1F *fHistBaryonEt; /** Et of identified baryons */
190 TH1F *fHistAntiBaryonEt; /** Et of identified anti-baryons */
191 TH1F *fHistMesonEt; /** Et of identified mesons */
193 TH1F *fHistBaryonEtAcc; /** Et of identified baryons in calorimeter acceptance */
194 TH1F *fHistAntiBaryonEtAcc; /** Et of identified anti-baryons in calorimeter acceptance */
195 TH1F *fHistMesonEtAcc; /** Et of identified mesons in calorimeter acceptance */
197 TH1F *fHistProtonEt; /** Et of identified protons */
198 TH1F *fHistChargedKaonEt; /** Et of identified charged kaons */
199 TH1F *fHistMuonEt; /** Et of identified muons */
200 TH1F *fHistElectronEt; /** Et of identified electrons */
202 TH1F *fHistProtonEtAcc; /** Et of identified protons in calorimeter acceptance */
203 TH1F *fHistChargedKaonEtAcc; /** Et of identified charged kaons in calorimeter acceptance */
204 TH1F *fHistMuonEtAcc; /** Et of identified muons in calorimeter acceptance */
205 TH1F *fHistElectronEtAcc; /** Et of identified electrons in calorimeter acceptance */
207 /* Correction plots */
208 TH2F *fHistEtRecvsEtMC; //Reconstructed Et versus MC Et
210 TH1F *fHistTMDeltaR; /* Track matching plots */
212 TTree *fTree; // optional TTree
215 AliESDtrackCuts* fEsdtrackCutsTPC;//esd track cuts for TPC tracks (which may also contain ITS hits)
218 //Declare private to avoid compilation warning
219 AliAnalysisEt & operator = (const AliAnalysisEt & g) ;//cpy assignment
220 AliAnalysisEt(const AliAnalysisEt & g) ; // cpy ctor
222 ClassDef(AliAnalysisEt, 0);
225 #endif // ALIANALYSISET_H