]>
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; } |
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; //! |
168 | ||
c228a5af | 169 | //=========================================================================================== |
170 | // Histograms | |
171 | //=========================================================================================== | |
172 | ||
bcb6ffdb | 173 | TH1F *fHistV0MultiplicityBeforeTrigSel; //! V0 multiplicity distribution |
174 | TH1F *fHistV0MultiplicityForTrigEvt; //! V0 multiplicity distribution | |
175 | TH1F *fHistV0MultiplicityForSelEvt; //! V0 multiplicity distribution | |
176 | TH1F *fHistV0MultiplicityForSelEvtNoTPCOnly; //! V0 multiplicity distribution | |
177 | TH1F *fHistV0MultiplicityForSelEvtNoTPCOnlyNoPileup; //! V0 multiplicity distribution | |
178 | ||
179 | TH1F *fHistMultiplicityBeforeTrigSel; //! multiplicity distribution | |
180 | TH1F *fHistMultiplicityForTrigEvt; //! multiplicity distribution | |
181 | TH1F *fHistMultiplicity; //! multiplicity distribution | |
182 | TH1F *fHistMultiplicityNoTPCOnly; //! multiplicity distribution | |
183 | TH1F *fHistMultiplicityNoTPCOnlyNoPileup; //! multiplicity distribution | |
184 | ||
1d94e33b | 185 | //Raw Data for J/Psi paper Technique |
186 | TH2F *f2dHistMultiplicityVsTrueBeforeTrigSel; //! multiplicity distribution | |
187 | TH2F *f2dHistMultiplicityVsTrueForTrigEvt; //! multiplicity distribution | |
188 | TH2F *f2dHistMultiplicityVsTrue; //! multiplicity distribution | |
189 | TH2F *f2dHistMultiplicityVsTrueNoTPCOnly; //! multiplicity distribution | |
190 | TH2F *f2dHistMultiplicityVsTrueNoTPCOnlyNoPileup; //! multiplicity distribution | |
191 | ||
192 | //Raw Data for Vertex Z position estimator change | |
193 | TH2F *f2dHistMultiplicityVsVertexZBeforeTrigSel; //! multiplicity distribution | |
194 | TH2F *f2dHistMultiplicityVsVertexZForTrigEvt; //! multiplicity distribution | |
195 | TH2F *f2dHistMultiplicityVsVertexZ; //! multiplicity distribution | |
196 | TH2F *f2dHistMultiplicityVsVertexZNoTPCOnly; //! multiplicity distribution | |
197 | TH2F *f2dHistMultiplicityVsVertexZNoTPCOnlyNoPileup; //! multiplicity distribution | |
198 | ||
199 | TH1F *fHistGenVertexZBeforeTrigSel; //! multiplicity distribution | |
200 | TH1F *fHistGenVertexZForTrigEvt; //! multiplicity distribution | |
201 | TH1F *fHistGenVertexZ; //! multiplicity distribution | |
202 | TH1F *fHistGenVertexZNoTPCOnly; //! multiplicity distribution | |
203 | TH1F *fHistGenVertexZNoTPCOnlyNoPileup; //! multiplicity distribution | |
204 | ||
bcb6ffdb | 205 | //---> Filled At Analysis Scope |
206 | ||
207 | TH3F *f3dHistPrimAnalysisPtVsYVsMultLambda; //! Lambda | |
208 | TH3F *f3dHistPrimAnalysisPtVsYVsMultAntiLambda; //! AntiLambda | |
209 | TH3F *f3dHistPrimAnalysisPtVsYVsMultK0Short; //! K0Short | |
210 | ||
388b34d3 | 211 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultLambda; //! Lambda |
212 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultAntiLambda; //! AntiLambda | |
213 | TH3F *f3dHistPrimAnalysisPtVsYCMSVsMultK0Short; //! K0Short | |
214 | ||
bcb6ffdb | 215 | //---> Containers for monte carlo information for calculating efficiency! |
216 | ||
217 | TH3F *f3dHistPrimRawPtVsYVsMultLambda; //! Lambda | |
218 | TH3F *f3dHistPrimRawPtVsYVsMultAntiLambda; //! AntiLambda | |
219 | TH3F *f3dHistPrimRawPtVsYVsMultK0Short; //! K0Short | |
220 | ||
388b34d3 | 221 | TH3F *f3dHistPrimRawPtVsYCMSVsMultLambda; //! Lambda |
222 | TH3F *f3dHistPrimRawPtVsYCMSVsMultAntiLambda; //! AntiLambda | |
223 | TH3F *f3dHistPrimRawPtVsYCMSVsMultK0Short; //! K0Short | |
224 | ||
df661cb7 | 225 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjLambda; //! Non-injected Lambda |
226 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjAntiLambda; //! Non-injected AntiLambda | |
227 | TH3F *f3dHistPrimRawPtVsYVsMultNonInjK0Short; //! Non-injected K0Short | |
228 | ||
a18e0a39 | 229 | TH3F *f3dHistPrimRawPtVsYVsMultMCLambda; //! Lambda |
230 | TH3F *f3dHistPrimRawPtVsYVsMultMCAntiLambda; //! AntiLambda | |
231 | TH3F *f3dHistPrimRawPtVsYVsMultMCK0Short; //! K0Short | |
232 | ||
233 | TH3F *f3dHistPrimRawPtVsYVsVertexZLambda; //! Lambda | |
234 | TH3F *f3dHistPrimRawPtVsYVsVertexZAntiLambda; //! AntiLambda | |
235 | TH3F *f3dHistPrimRawPtVsYVsVertexZK0Short; //! K0Short | |
236 | ||
7e81bb32 | 237 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultLambda; //! Lambda |
238 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultAntiLambda; //! AntiLambda | |
239 | TH3F *f3dHistPrimCloseToPVPtVsYVsMultK0Short; //! K0Short | |
240 | ||
bcb6ffdb | 241 | //---> Filled vs Decay Length |
242 | ||
243 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthLambda; //! Lambda | |
244 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthAntiLambda; //! AntiLambda | |
245 | TH3F *f3dHistPrimRawPtVsYVsDecayLengthK0Short; //! K0Short | |
246 | ||
247 | //---> Needed for FeedDown Corrections | |
248 | ||
249 | TH3F *f3dHistGenPtVsYVsMultXiMinus; //! Generated Xi- Distrib | |
250 | TH3F *f3dHistGenPtVsYVsMultXiPlus; //! Generated Xi+ Distrib | |
c228a5af | 251 | TH3F *f3dHistGenPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib |
252 | TH3F *f3dHistGenPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib | |
253 | ||
254 | TH3F *f3dHistGenSelectedPtVsYVsMultXiMinus; //! Generated Xi- Distrib, at event selection level | |
255 | TH3F *f3dHistGenSelectedPtVsYVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level | |
256 | TH3F *f3dHistGenSelectedPtVsYVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level | |
257 | TH3F *f3dHistGenSelectedPtVsYVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level | |
e039e6b3 | 258 | |
259 | TH3F *f3dHistGenPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib | |
260 | TH3F *f3dHistGenPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib | |
261 | TH3F *f3dHistGenPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib | |
262 | TH3F *f3dHistGenPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib | |
263 | ||
264 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiMinus; //! Generated Xi- Distrib, at event selection level | |
265 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultXiPlus; //! Generated Xi+ Distrib, at event selection level | |
266 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaMinus; //! Generated Omega- Distrib, at event selection level | |
267 | TH3F *f3dHistGenSelectedPtVsYCMSVsMultOmegaPlus; //! Generated Omega+ Distrib, at event selection level | |
bcb6ffdb | 268 | |
269 | TH1F *fHistPVx; //! PVx distrib | |
270 | TH1F *fHistPVy; //! PVy distrib | |
271 | TH1F *fHistPVz; //! PVz distrib | |
272 | TH1F *fHistPVxAnalysis; //! PVx distrib | |
273 | TH1F *fHistPVyAnalysis; //! PVy distrib | |
274 | TH1F *fHistPVzAnalysis; //! PVz distrib | |
275 | TH1F *fHistPVxAnalysisHasHighPtLambda; //! PVx distrib | |
276 | TH1F *fHistPVyAnalysisHasHighPtLambda; //! PVy distrib | |
277 | TH1F *fHistPVzAnalysisHasHighPtLambda; //! PVz distrib | |
278 | ||
4f8ed367 | 279 | TH1F *fHistSwappedV0Counter; //! Swapped v0 counter |
280 | ||
bcb6ffdb | 281 | AliAnalysisTaskExtractPerformanceV0(const AliAnalysisTaskExtractPerformanceV0&); // not implemented |
282 | AliAnalysisTaskExtractPerformanceV0& operator=(const AliAnalysisTaskExtractPerformanceV0&); // not implemented | |
283 | ||
284 | ClassDef(AliAnalysisTaskExtractPerformanceV0, 11); | |
285 | }; | |
286 | ||
287 | #endif |