]>
Commit | Line | Data |
---|---|---|
bcb6ffdb | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | // +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ | |
17 | // | |
18 | // Modified version of AliAnalysisTaskCheckCascade.h | |
19 | // Used bits of code from AliAnalysisTaskCheckPerformanceStrange | |
20 | // | |
21 | // --- David Dobrigkeit Chinellato | |
22 | // | |
23 | // +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+ | |
24 | ||
25 | #ifndef ALIANALYSISTASKEXTRACTPERFORMANCEV0_H | |
26 | #define ALIANALYSISTASKEXTRACTPERFORMANCEV0_H | |
27 | ||
28 | class TList; | |
29 | class TH1F; | |
30 | class TH2F; | |
31 | class TH3F; | |
32 | class TVector3; | |
33 | class THnSparse; | |
34 | ||
35 | class AliESDpid; | |
36 | class AliESDtrackCuts; | |
37 | class AliESDEvent; | |
38 | class AliPhysicsSelection; | |
39 | class AliCFContainer; | |
40 | ||
41 | //#include "TString.h" | |
42 | //#include "AliESDtrackCuts.h" | |
43 | #include "AliAnalysisTaskSE.h" | |
44 | ||
45 | class AliAnalysisTaskExtractPerformanceV0 : public AliAnalysisTaskSE { | |
46 | public: | |
47 | AliAnalysisTaskExtractPerformanceV0(); | |
48 | AliAnalysisTaskExtractPerformanceV0(const char *name); | |
49 | virtual ~AliAnalysisTaskExtractPerformanceV0(); | |
50 | ||
51 | virtual void UserCreateOutputObjects(); | |
52 | virtual void UserExec(Option_t *option); | |
53 | virtual void Terminate(Option_t *); | |
54 | Double_t MyRapidity(Double_t rE, Double_t rPz) const; | |
4f8ed367 | 55 | void CheckChargeV0(AliESDv0 *thisv0); |
548bd092 | 56 | |
57 | void SetIsNuclear (Bool_t lIsNuclear = kTRUE ) { fkIsNuclear = lIsNuclear; } | |
048f4f8f | 58 | void SetINT7Trigger (Bool_t lSwitchINT7 = kTRUE ) { fkSwitchINT7 = lSwitchINT7; } |
548bd092 | 59 | void SetUseOnTheFly (Bool_t lUseOnTheFly = kTRUE ) { fkUseOnTheFly = lUseOnTheFly; } |
7f4db1cd | 60 | void SetTakeAllTracks (Bool_t lTakeAllTracks = kTRUE ) { fkTakeAllTracks = lTakeAllTracks; } |
388b34d3 | 61 | void SetpARapidityShift (Double_t lRapShift = 0.465 ) { fpArapidityShift = lRapShift; } |
e039e6b3 | 62 | void SetCentralityEstimator (TString lCentralityEstimator = "V0M" ) { fCentralityEstimator = lCentralityEstimator; } |
c92e2e56 | 63 | void SetLightWeightAnalysis (Bool_t lLightWeight = kTRUE) {fkLightWeight = lLightWeight; } |
bcb6ffdb | 64 | |
65 | private: | |
66 | // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ... | |
67 | // your data member object is created on the worker nodes and streaming is not needed. | |
68 | // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14 | |
69 | TList *fListHistV0; //! List of Cascade histograms | |
c228a5af | 70 | TTree *fTree; //! Output Tree, V0 |
bcb6ffdb | 71 | |
548bd092 | 72 | AliPIDResponse *fPIDResponse; // PID response object |
76029adc | 73 | AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition |
9083ce79 | 74 | |
548bd092 | 75 | |
76 | //Objects Controlling Task Behaviour | |
77 | ||
78 | Bool_t fkIsNuclear; //if true, replace multiplicity est. by centrality (default FALSE) | |
048f4f8f | 79 | Bool_t fkSwitchINT7 ; //if true, skip FASTOnly (default FALSE) |
548bd092 | 80 | Bool_t fkUseOnTheFly; //if true, will use On-the-fly V0s instead of Offline V0s (default FALSE) |
7f4db1cd | 81 | Bool_t fkTakeAllTracks; // if true, no TPC crossed rows and ratio cut |
9083ce79 | 82 | Double_t fpArapidityShift; //pA rapidity shift (should be 0.465, usually) |
e039e6b3 | 83 | TString fCentralityEstimator; //Centrality Estimator String value (default V0M) |
c92e2e56 | 84 | |
85 | Bool_t fkLightWeight; //if true, skip a number of debugging information branches in TTree | |
86 | //(to make resulting tree output significantly smaller! | |
bcb6ffdb | 87 | |
c228a5af | 88 | //=========================================================================================== |
89 | // Variables for Tree, V0s | |
90 | //=========================================================================================== | |
bcb6ffdb | 91 | Int_t fTreeVariablePrimaryStatus; //! |
92 | Int_t fTreeVariablePrimaryStatusMother; //! | |
93 | Float_t fTreeVariableChi2V0; //! | |
94 | Float_t fTreeVariableDcaV0Daughters; //! | |
95 | Float_t fTreeVariableDcaV0ToPrimVertex; //! | |
96 | Float_t fTreeVariableDcaPosToPrimVertex; //! | |
97 | Float_t fTreeVariableDcaNegToPrimVertex; //! | |
98 | Float_t fTreeVariableV0CosineOfPointingAngle; //! | |
99 | Float_t fTreeVariableV0Radius; //! | |
100 | Float_t fTreeVariablePt; //! | |
101 | Float_t fTreeVariablePtMC; //! | |
102 | Float_t fTreeVariableRapK0Short; //! | |
103 | Float_t fTreeVariableRapLambda; //! | |
104 | Float_t fTreeVariableRapMC; //! | |
105 | Float_t fTreeVariableInvMassK0s; //! | |
106 | Float_t fTreeVariableInvMassLambda; //! | |
107 | Float_t fTreeVariableInvMassAntiLambda; //! | |
108 | Float_t fTreeVariableAlphaV0; //! | |
109 | Float_t fTreeVariablePtArmV0;//! | |
110 | Float_t fTreeVariableNegTotMomentum; //! | |
111 | Float_t fTreeVariablePosTotMomentum; //! | |
112 | Float_t fTreeVariableNegTransvMomentum; //! | |
113 | Float_t fTreeVariablePosTransvMomentum; //! | |
114 | Float_t fTreeVariableNegTransvMomentumMC; //! | |
115 | Float_t fTreeVariablePosTransvMomentumMC; //! | |
116 | ||
117 | Float_t fTreeVariableNSigmasPosProton; //! | |
118 | Float_t fTreeVariableNSigmasPosPion; //! | |
119 | Float_t fTreeVariableNSigmasNegProton; //! | |
120 | Float_t fTreeVariableNSigmasNegPion; //! | |
121 | ||
122 | Float_t fTreeVariablePtMother; //! | |
123 | Float_t fTreeVariableV0CreationRadius; //! | |
124 | Int_t fTreeVariablePID; //! | |
125 | Int_t fTreeVariablePIDPositive; //! | |
126 | Int_t fTreeVariablePIDNegative; //! | |
127 | Int_t fTreeVariablePIDMother; //! | |
128 | Int_t fTreeVariableIndexStatus; //! | |
129 | Int_t fTreeVariableIndexStatusMother; //! | |
130 | ||
0eebce93 | 131 | Int_t fTreeVariableRunNumber; //! |
132 | ULong64_t fTreeVariableEventNumber; //! | |
133 | ||
bcb6ffdb | 134 | //Note: TDistOverTotMom needs a mass hypothesis to be converted to proper decaylength. |
135 | Float_t fTreeVariableDistOverTotMom;//! | |
136 | ||
137 | Float_t fTreeVariablePosEta; //! | |
138 | Float_t fTreeVariableNegEta; //! | |
139 | ||
a18e0a39 | 140 | Float_t fTreeVariableVertexZ; //! |
141 | ||
bcb6ffdb | 142 | Int_t fTreeVariableLeastNbrCrossedRows;//! |
143 | Float_t fTreeVariableLeastRatioCrossedRowsOverFindable;//! | |
5d43de53 | 144 | |
145 | Int_t fTreeVariableMultiplicity;//! | |
146 | Int_t fTreeVariableMultiplicityV0A;//! | |
147 | Int_t fTreeVariableMultiplicityZNA;//! | |
148 | Int_t fTreeVariableMultiplicityTRK;//! | |
149 | Int_t fTreeVariableMultiplicitySPD;//! | |
150 | ||
a18e0a39 | 151 | Int_t fTreeVariableMultiplicityMC;//! |
bcb6ffdb | 152 | |
11960ef4 | 153 | Float_t fTreeVariableV0x; //! |
154 | Float_t fTreeVariableV0y; //! | |
155 | Float_t fTreeVariableV0z; //! | |
156 | ||
157 | Float_t fTreeVariableV0Px; //! | |
158 | Float_t fTreeVariableV0Py; //! | |
159 | Float_t fTreeVariableV0Pz; //! | |
160 | ||
161 | Float_t fTreeVariableMCV0x; //! | |
162 | Float_t fTreeVariableMCV0y; //! | |
163 | Float_t fTreeVariableMCV0z; //! | |
164 | ||
165 | Float_t fTreeVariableMCV0Px; //! | |
166 | Float_t fTreeVariableMCV0Py; //! | |
167 | Float_t fTreeVariableMCV0Pz; //! | |
168 | ||
653265df | 169 | Float_t fTreeVariablePVx; //! |
170 | Float_t fTreeVariablePVy; //! | |
171 | Float_t fTreeVariablePVz; //! | |
172 | ||
173 | Float_t fTreeVariableMCPVx; //! | |
174 | Float_t fTreeVariableMCPVy; //! | |
175 | Float_t fTreeVariableMCPVz; //! | |
176 | ||
df661cb7 | 177 | Bool_t fTreeVariableIsNonInjected; //! |
f068724a | 178 | //Decay Length issue debugging: ULong_t with track status |
179 | ULong64_t fTreeVariableNegTrackStatus; | |
180 | ULong64_t fTreeVariablePosTrackStatus; | |
df661cb7 | 181 | |
c228a5af | 182 | //=========================================================================================== |
183 | // Histograms | |
184 | //=========================================================================================== | |
185 | ||
bcb6ffdb | 186 | TH1F *fHistV0MultiplicityBeforeTrigSel; //! V0 multiplicity distribution |
187 | TH1F *fHistV0MultiplicityForTrigEvt; //! V0 multiplicity distribution | |
188 | TH1F *fHistV0MultiplicityForSelEvt; //! V0 multiplicity distribution | |
189 | TH1F *fHistV0MultiplicityForSelEvtNoTPCOnly; //! V0 multiplicity distribution | |
190 | TH1F *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution | |
191 | ||
5d43de53 | 192 | //Default V0M Centrality |
bcb6ffdb | 193 | TH1F *fHistMultiplicityBeforeTrigSel; //! multiplicity distribution |
194 | TH1F *fHistMultiplicityForTrigEvt; //! multiplicity distribution | |
195 | TH1F *fHistMultiplicity; //! multiplicity distribution | |
196 | TH1F *fHistMultiplicityNoTPCOnly; //! multiplicity distribution | |
197 | TH1F *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution | |
198 | ||
5d43de53 | 199 | //V0A Centrality |
200 | TH1F *fHistMultiplicityV0ABeforeTrigSel; //! multiplicity distribution | |
201 | TH1F *fHistMultiplicityV0AForTrigEvt; //! multiplicity distribution | |
202 | TH1F *fHistMultiplicityV0A; //! multiplicity distribution | |
203 | TH1F *fHistMultiplicityV0ANoTPCOnly; //! multiplicity distribution | |
204 | TH1F *fHistMultiplicityV0ANoTPCOnlyNoPileup; //! multiplicity distribution | |
205 | ||
206 | //ZNA Centrality | |
207 | TH1F *fHistMultiplicityZNABeforeTrigSel; //! multiplicity distribution | |
208 | TH1F *fHistMultiplicityZNAForTrigEvt; //! multiplicity distribution | |
209 | TH1F *fHistMultiplicityZNA; //! multiplicity distribution | |
210 | TH1F *fHistMultiplicityZNANoTPCOnly; //! multiplicity distribution | |
211 | TH1F *fHistMultiplicityZNANoTPCOnlyNoPileup; //! multiplicity distribution | |
212 | ||
213 | //TRK Centrality | |
214 | TH1F *fHistMultiplicityTRKBeforeTrigSel; //! multiplicity distribution | |
215 | TH1F *fHistMultiplicityTRKForTrigEvt; //! multiplicity distribution | |
216 | TH1F *fHistMultiplicityTRK; //! multiplicity distribution | |
217 | TH1F *fHistMultiplicityTRKNoTPCOnly; //! multiplicity distribution | |
218 | TH1F *fHistMultiplicityTRKNoTPCOnlyNoPileup; //! multiplicity distribution | |
219 | ||
220 | //SPD Centrality | |
221 | TH1F *fHistMultiplicitySPDBeforeTrigSel; //! multiplicity distribution | |
222 | TH1F *fHistMultiplicitySPDForTrigEvt; //! multiplicity distribution | |
223 | TH1F *fHistMultiplicitySPD; //! multiplicity distribution | |
224 | TH1F *fHistMultiplicitySPDNoTPCOnly; //! multiplicity distribution | |
225 | TH1F *fHistMultiplicitySPDNoTPCOnlyNoPileup; //! multiplicity distribution | |
226 | ||
1d94e33b | 227 | //Raw Data for J/Psi paper Technique |
228 | TH2F *f2dHistMultiplicityVsTrueBeforeTrigSel; //! multiplicity distribution | |
229 | TH2F *f2dHistMultiplicityVsTrueForTrigEvt; //! multiplicity distribution | |
230 | TH2F *f2dHistMultiplicityVsTrue; //! multiplicity distribution | |
231 | TH2F *f2dHistMultiplicityVsTrueNoTPCOnly; //! multiplicity distribution | |
232 | TH2F *f2dHistMultiplicityVsTrueNoTPCOnlyNoPileup; //! multiplicity distribution | |
233 | ||
234 | //Raw Data for Vertex Z position estimator change | |
235 | TH2F *f2dHistMultiplicityVsVertexZBeforeTrigSel; //! multiplicity distribution | |
236 | TH2F *f2dHistMultiplicityVsVertexZForTrigEvt; //! multiplicity distribution | |
237 | TH2F *f2dHistMultiplicityVsVertexZ; //! multiplicity distribution | |
238 | TH2F *f2dHistMultiplicityVsVertexZNoTPCOnly; //! multiplicity distribution | |
239 | TH2F *f2dHistMultiplicityVsVertexZNoTPCOnlyNoPileup; //! multiplicity distribution | |
240 | ||
241 | TH1F *fHistGenVertexZBeforeTrigSel; //! multiplicity distribution | |
242 | TH1F *fHistGenVertexZForTrigEvt; //! multiplicity distribution | |
243 | TH1F *fHistGenVertexZ; //! multiplicity distribution | |
244 | TH1F *fHistGenVertexZNoTPCOnly; //! multiplicity distribution | |
245 | TH1F *fHistGenVertexZNoTPCOnlyNoPileup; //! multiplicity distribution | |
246 | ||
bcb6ffdb | 247 | //---> Filled At Analysis Scope |
248 | ||
249 | TH3F *f3dHistPrimAnalysisPtVsYVsMultLambda; //! Lambda | |
250 | TH3F *f3dHistPrimAnalysisPtVsYVsMultAntiLambda; //! AntiLambda | |
251 | TH3F *f3dHistPrimAnalysisPtVsYVsMultK0Short; //! K0Short | |
252 | ||
388b34d3 | 253 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultLambda; //! Lambda |
254 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultAntiLambda; //! AntiLambda | |
255 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultK0Short; //! K0Short | |
256 | ||
bcb6ffdb | 257 | //---> Containers for monte carlo information for calculating efficiency! |
258 | ||
259 | TH3F *f3dHistPrimRawPtVsYVsMultLambda; //! Lambda | |
260 | TH3F *f3dHistPrimRawPtVsYVsMultAntiLambda; //! AntiLambda | |
261 | TH3F *f3dHistPrimRawPtVsYVsMultK0Short; //! K0Short | |
262 | ||
388b34d3 | 263 | TH3F *f3dHistPrimRawPtVsYCMSVsMultLambda; //! Lambda |
264 | TH3F *f3dHistPrimRawPtVsYCMSVsMultAntiLambda; //! AntiLambda | |
265 | TH3F *f3dHistPrimRawPtVsYCMSVsMultK0Short; //! K0Short | |
266 | ||
df661cb7 | 267 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjLambda; //! Non-injected Lambda |
268 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjAntiLambda; //! Non-injected AntiLambda | |
269 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjK0Short; //! Non-injected K0Short | |
270 | ||
a18e0a39 | 271 | TH3F *f3dHistPrimRawPtVsYVsMultMCLambda; //! Lambda |
272 | TH3F *f3dHistPrimRawPtVsYVsMultMCAntiLambda; //! AntiLambda | |
273 | TH3F *f3dHistPrimRawPtVsYVsMultMCK0Short; //! K0Short | |
274 | ||
275 | TH3F *f3dHistPrimRawPtVsYVsVertexZLambda; //! Lambda | |
276 | TH3F *f3dHistPrimRawPtVsYVsVertexZAntiLambda; //! AntiLambda | |
277 | TH3F *f3dHistPrimRawPtVsYVsVertexZK0Short; //! K0Short | |
278 | ||
7e81bb32 | 279 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultLambda; //! Lambda |
280 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultAntiLambda; //! AntiLambda | |
281 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultK0Short; //! K0Short | |
282 | ||
bcb6ffdb | 283 | //---> Filled vs Decay Length |
284 | ||
285 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthLambda; //! Lambda | |
286 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthAntiLambda; //! AntiLambda | |
287 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthK0Short; //! K0Short | |
288 | ||
289 | //---> Needed for FeedDown Corrections | |
290 | ||
291 | TH3F *f3dHistGenPtVsYVsMultXiMinus; //! Generated Xi- Distrib | |
292 | TH3F *f3dHistGenPtVsYVsMultXiPlus; //! Generated Xi+ Distrib | |
c228a5af | 293 | TH3F *f3dHistGenPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib |
294 | TH3F *f3dHistGenPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib | |
295 | ||
296 | TH3F *f3dHistGenSelectedPtVsYVsMultXiMinus; //! Generated Xi- Distrib, at event selection level | |
297 | TH3F *f3dHistGenSelectedPtVsYVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level | |
298 | TH3F *f3dHistGenSelectedPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level | |
299 | TH3F *f3dHistGenSelectedPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level | |
e039e6b3 | 300 | |
301 | TH3F *f3dHistGenPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib | |
302 | TH3F *f3dHistGenPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib | |
303 | TH3F *f3dHistGenPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib | |
304 | TH3F *f3dHistGenPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib | |
305 | ||
306 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib, at event selection level | |
307 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level | |
308 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level | |
309 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level | |
bcb6ffdb | 310 | |
311 | TH1F *fHistPVx; //! PVx distrib | |
312 | TH1F *fHistPVy; //! PVy distrib | |
313 | TH1F *fHistPVz; //! PVz distrib | |
314 | TH1F *fHistPVxAnalysis; //! PVx distrib | |
315 | TH1F *fHistPVyAnalysis; //! PVy distrib | |
316 | TH1F *fHistPVzAnalysis; //! PVz distrib | |
317 | TH1F *fHistPVxAnalysisHasHighPtLambda; //! PVx distrib | |
318 | TH1F *fHistPVyAnalysisHasHighPtLambda; //! PVy distrib | |
319 | TH1F *fHistPVzAnalysisHasHighPtLambda; //! PVz distrib | |
320 | ||
4f8ed367 | 321 | TH1F *fHistSwappedV0Counter; //! Swapped v0 counter |
322 | ||
bcb6ffdb | 323 | AliAnalysisTaskExtractPerformanceV0(const AliAnalysisTaskExtractPerformanceV0&); // not implemented |
324 | AliAnalysisTaskExtractPerformanceV0& operator=(const AliAnalysisTaskExtractPerformanceV0&); // not implemented | |
325 | ||
326 | ClassDef(AliAnalysisTaskExtractPerformanceV0, 11); | |
327 | }; | |
328 | ||
329 | #endif |