]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/STRANGENESS/LambdaK0/AliAnalysisTaskExtractPerformanceV0.h
Small modifications for getting track status flags
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / LambdaK0 / AliAnalysisTaskExtractPerformanceV0.h
CommitLineData
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
28class TList;
29class TH1F;
30class TH2F;
31class TH3F;
32class TVector3;
33class THnSparse;
34
35class AliESDpid;
36class AliESDtrackCuts;
37class AliESDEvent;
38class AliPhysicsSelection;
39class AliCFContainer;
40
41//#include "TString.h"
42//#include "AliESDtrackCuts.h"
43#include "AliAnalysisTaskSE.h"
44
45class 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; }
bcb6ffdb 63
64 private:
65 // Note : In ROOT, "//!" means "do not stream the data from Master node to Worker node" ...
66 // your data member object is created on the worker nodes and streaming is not needed.
67 // http://root.cern.ch/download/doc/11InputOutput.pdf, page 14
68 TList *fListHistV0; //! List of Cascade histograms
c228a5af 69 TTree *fTree; //! Output Tree, V0
bcb6ffdb 70
548bd092 71 AliPIDResponse *fPIDResponse; // PID response object
76029adc 72 AliESDtrackCuts *fESDtrackCuts; // ESD track cuts used for primary track definition
9083ce79 73
548bd092 74
75 //Objects Controlling Task Behaviour
76
77 Bool_t fkIsNuclear; //if true, replace multiplicity est. by centrality (default FALSE)
048f4f8f 78 Bool_t fkSwitchINT7 ; //if true, skip FASTOnly (default FALSE)
548bd092 79 Bool_t fkUseOnTheFly; //if true, will use On-the-fly V0s instead of Offline V0s (default FALSE)
7f4db1cd 80 Bool_t fkTakeAllTracks; // if true, no TPC crossed rows and ratio cut
9083ce79 81 Double_t fpArapidityShift; //pA rapidity shift (should be 0.465, usually)
e039e6b3 82 TString fCentralityEstimator; //Centrality Estimator String value (default V0M)
bcb6ffdb 83
c228a5af 84//===========================================================================================
85// Variables for Tree, V0s
86//===========================================================================================
bcb6ffdb 87 Int_t fTreeVariablePrimaryStatus; //!
88 Int_t fTreeVariablePrimaryStatusMother; //!
89 Float_t fTreeVariableChi2V0; //!
90 Float_t fTreeVariableDcaV0Daughters; //!
91 Float_t fTreeVariableDcaV0ToPrimVertex; //!
92 Float_t fTreeVariableDcaPosToPrimVertex; //!
93 Float_t fTreeVariableDcaNegToPrimVertex; //!
94 Float_t fTreeVariableV0CosineOfPointingAngle; //!
95 Float_t fTreeVariableV0Radius; //!
96 Float_t fTreeVariablePt; //!
97 Float_t fTreeVariablePtMC; //!
98 Float_t fTreeVariableRapK0Short; //!
99 Float_t fTreeVariableRapLambda; //!
100 Float_t fTreeVariableRapMC; //!
101 Float_t fTreeVariableInvMassK0s; //!
102 Float_t fTreeVariableInvMassLambda; //!
103 Float_t fTreeVariableInvMassAntiLambda; //!
104 Float_t fTreeVariableAlphaV0; //!
105 Float_t fTreeVariablePtArmV0;//!
106 Float_t fTreeVariableNegTotMomentum; //!
107 Float_t fTreeVariablePosTotMomentum; //!
108 Float_t fTreeVariableNegTransvMomentum; //!
109 Float_t fTreeVariablePosTransvMomentum; //!
110 Float_t fTreeVariableNegTransvMomentumMC; //!
111 Float_t fTreeVariablePosTransvMomentumMC; //!
112
113 Float_t fTreeVariableNSigmasPosProton; //!
114 Float_t fTreeVariableNSigmasPosPion; //!
115 Float_t fTreeVariableNSigmasNegProton; //!
116 Float_t fTreeVariableNSigmasNegPion; //!
117
118 Float_t fTreeVariablePtMother; //!
119 Float_t fTreeVariableV0CreationRadius; //!
120 Int_t fTreeVariablePID; //!
121 Int_t fTreeVariablePIDPositive; //!
122 Int_t fTreeVariablePIDNegative; //!
123 Int_t fTreeVariablePIDMother; //!
124 Int_t fTreeVariableIndexStatus; //!
125 Int_t fTreeVariableIndexStatusMother; //!
126
0eebce93 127 Int_t fTreeVariableRunNumber; //!
128 ULong64_t fTreeVariableEventNumber; //!
129
bcb6ffdb 130 //Note: TDistOverTotMom needs a mass hypothesis to be converted to proper decaylength.
131 Float_t fTreeVariableDistOverTotMom;//!
132
133 Float_t fTreeVariablePosEta; //!
134 Float_t fTreeVariableNegEta; //!
135
a18e0a39 136 Float_t fTreeVariableVertexZ; //!
137
bcb6ffdb 138 Int_t fTreeVariableLeastNbrCrossedRows;//!
139 Float_t fTreeVariableLeastRatioCrossedRowsOverFindable;//!
140 Int_t fTreeVariableMultiplicity;//!
a18e0a39 141 Int_t fTreeVariableMultiplicityMC;//!
bcb6ffdb 142
11960ef4 143 Float_t fTreeVariableV0x; //!
144 Float_t fTreeVariableV0y; //!
145 Float_t fTreeVariableV0z; //!
146
147 Float_t fTreeVariableV0Px; //!
148 Float_t fTreeVariableV0Py; //!
149 Float_t fTreeVariableV0Pz; //!
150
151 Float_t fTreeVariableMCV0x; //!
152 Float_t fTreeVariableMCV0y; //!
153 Float_t fTreeVariableMCV0z; //!
154
155 Float_t fTreeVariableMCV0Px; //!
156 Float_t fTreeVariableMCV0Py; //!
157 Float_t fTreeVariableMCV0Pz; //!
158
653265df 159 Float_t fTreeVariablePVx; //!
160 Float_t fTreeVariablePVy; //!
161 Float_t fTreeVariablePVz; //!
162
163 Float_t fTreeVariableMCPVx; //!
164 Float_t fTreeVariableMCPVy; //!
165 Float_t fTreeVariableMCPVz; //!
166
df661cb7 167 Bool_t fTreeVariableIsNonInjected; //!
f068724a 168 //Decay Length issue debugging: ULong_t with track status
169 ULong64_t fTreeVariableNegTrackStatus;
170 ULong64_t fTreeVariablePosTrackStatus;
df661cb7 171
c228a5af 172//===========================================================================================
173// Histograms
174//===========================================================================================
175
bcb6ffdb 176 TH1F *fHistV0MultiplicityBeforeTrigSel; //! V0 multiplicity distribution
177 TH1F *fHistV0MultiplicityForTrigEvt; //! V0 multiplicity distribution
178 TH1F *fHistV0MultiplicityForSelEvt; //! V0 multiplicity distribution
179 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnly; //! V0 multiplicity distribution
180 TH1F *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution
181
182 TH1F *fHistMultiplicityBeforeTrigSel; //! multiplicity distribution
183 TH1F *fHistMultiplicityForTrigEvt; //! multiplicity distribution
184 TH1F *fHistMultiplicity; //! multiplicity distribution
185 TH1F *fHistMultiplicityNoTPCOnly; //! multiplicity distribution
186 TH1F *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution
187
1d94e33b 188 //Raw Data for J/Psi paper Technique
189 TH2F *f2dHistMultiplicityVsTrueBeforeTrigSel; //! multiplicity distribution
190 TH2F *f2dHistMultiplicityVsTrueForTrigEvt; //! multiplicity distribution
191 TH2F *f2dHistMultiplicityVsTrue; //! multiplicity distribution
192 TH2F *f2dHistMultiplicityVsTrueNoTPCOnly; //! multiplicity distribution
193 TH2F *f2dHistMultiplicityVsTrueNoTPCOnlyNoPileup; //! multiplicity distribution
194
195 //Raw Data for Vertex Z position estimator change
196 TH2F *f2dHistMultiplicityVsVertexZBeforeTrigSel; //! multiplicity distribution
197 TH2F *f2dHistMultiplicityVsVertexZForTrigEvt; //! multiplicity distribution
198 TH2F *f2dHistMultiplicityVsVertexZ; //! multiplicity distribution
199 TH2F *f2dHistMultiplicityVsVertexZNoTPCOnly; //! multiplicity distribution
200 TH2F *f2dHistMultiplicityVsVertexZNoTPCOnlyNoPileup; //! multiplicity distribution
201
202 TH1F *fHistGenVertexZBeforeTrigSel; //! multiplicity distribution
203 TH1F *fHistGenVertexZForTrigEvt; //! multiplicity distribution
204 TH1F *fHistGenVertexZ; //! multiplicity distribution
205 TH1F *fHistGenVertexZNoTPCOnly; //! multiplicity distribution
206 TH1F *fHistGenVertexZNoTPCOnlyNoPileup; //! multiplicity distribution
207
bcb6ffdb 208//---> Filled At Analysis Scope
209
210 TH3F *f3dHistPrimAnalysisPtVsYVsMultLambda; //! Lambda
211 TH3F *f3dHistPrimAnalysisPtVsYVsMultAntiLambda; //! AntiLambda
212 TH3F *f3dHistPrimAnalysisPtVsYVsMultK0Short; //! K0Short
213
388b34d3 214 TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultLambda; //! Lambda
215 TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultAntiLambda; //! AntiLambda
216 TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultK0Short; //! K0Short
217
bcb6ffdb 218//---> Containers for monte carlo information for calculating efficiency!
219
220 TH3F *f3dHistPrimRawPtVsYVsMultLambda; //! Lambda
221 TH3F *f3dHistPrimRawPtVsYVsMultAntiLambda; //! AntiLambda
222 TH3F *f3dHistPrimRawPtVsYVsMultK0Short; //! K0Short
223
388b34d3 224 TH3F *f3dHistPrimRawPtVsYCMSVsMultLambda; //! Lambda
225 TH3F *f3dHistPrimRawPtVsYCMSVsMultAntiLambda; //! AntiLambda
226 TH3F *f3dHistPrimRawPtVsYCMSVsMultK0Short; //! K0Short
227
df661cb7 228 TH3F *f3dHistPrimRawPtVsYVsMultNonInjLambda; //! Non-injected Lambda
229 TH3F *f3dHistPrimRawPtVsYVsMultNonInjAntiLambda; //! Non-injected AntiLambda
230 TH3F *f3dHistPrimRawPtVsYVsMultNonInjK0Short; //! Non-injected K0Short
231
a18e0a39 232 TH3F *f3dHistPrimRawPtVsYVsMultMCLambda; //! Lambda
233 TH3F *f3dHistPrimRawPtVsYVsMultMCAntiLambda; //! AntiLambda
234 TH3F *f3dHistPrimRawPtVsYVsMultMCK0Short; //! K0Short
235
236 TH3F *f3dHistPrimRawPtVsYVsVertexZLambda; //! Lambda
237 TH3F *f3dHistPrimRawPtVsYVsVertexZAntiLambda; //! AntiLambda
238 TH3F *f3dHistPrimRawPtVsYVsVertexZK0Short; //! K0Short
239
7e81bb32 240 TH3F *f3dHistPrimCloseToPVPtVsYVsMultLambda; //! Lambda
241 TH3F *f3dHistPrimCloseToPVPtVsYVsMultAntiLambda; //! AntiLambda
242 TH3F *f3dHistPrimCloseToPVPtVsYVsMultK0Short; //! K0Short
243
bcb6ffdb 244//---> Filled vs Decay Length
245
246 TH3F *f3dHistPrimRawPtVsYVsDecayLengthLambda; //! Lambda
247 TH3F *f3dHistPrimRawPtVsYVsDecayLengthAntiLambda; //! AntiLambda
248 TH3F *f3dHistPrimRawPtVsYVsDecayLengthK0Short; //! K0Short
249
250//---> Needed for FeedDown Corrections
251
252 TH3F *f3dHistGenPtVsYVsMultXiMinus; //! Generated Xi- Distrib
253 TH3F *f3dHistGenPtVsYVsMultXiPlus; //! Generated Xi+ Distrib
c228a5af 254 TH3F *f3dHistGenPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib
255 TH3F *f3dHistGenPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib
256
257 TH3F *f3dHistGenSelectedPtVsYVsMultXiMinus; //! Generated Xi- Distrib, at event selection level
258 TH3F *f3dHistGenSelectedPtVsYVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level
259 TH3F *f3dHistGenSelectedPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level
260 TH3F *f3dHistGenSelectedPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level
e039e6b3 261
262 TH3F *f3dHistGenPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib
263 TH3F *f3dHistGenPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib
264 TH3F *f3dHistGenPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib
265 TH3F *f3dHistGenPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib
266
267 TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib, at event selection level
268 TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level
269 TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level
270 TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level
bcb6ffdb 271
272 TH1F *fHistPVx; //! PVx distrib
273 TH1F *fHistPVy; //! PVy distrib
274 TH1F *fHistPVz; //! PVz distrib
275 TH1F *fHistPVxAnalysis; //! PVx distrib
276 TH1F *fHistPVyAnalysis; //! PVy distrib
277 TH1F *fHistPVzAnalysis; //! PVz distrib
278 TH1F *fHistPVxAnalysisHasHighPtLambda; //! PVx distrib
279 TH1F *fHistPVyAnalysisHasHighPtLambda; //! PVy distrib
280 TH1F *fHistPVzAnalysisHasHighPtLambda; //! PVz distrib
281
4f8ed367 282 TH1F *fHistSwappedV0Counter; //! Swapped v0 counter
283
bcb6ffdb 284 AliAnalysisTaskExtractPerformanceV0(const AliAnalysisTaskExtractPerformanceV0&); // not implemented
285 AliAnalysisTaskExtractPerformanceV0& operator=(const AliAnalysisTaskExtractPerformanceV0&); // not implemented
286
287 ClassDef(AliAnalysisTaskExtractPerformanceV0, 11);
288};
289
290#endif