]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGLF/STRANGENESS/Cascades/AliAnalysisTaskStrangenessVsMultiplicityMC.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Introduce missing L4S4 chamber in sector 17
[u/mrichter/AliRoot.git] / PWGLF / STRANGENESS / Cascades / AliAnalysisTaskStrangenessVsMultiplicityMC.cxx
CommitLineData
af0d4791 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// This task is meant to explore the possibility of using a VZERO amplitude
19// based multiplicity estimator for proton-proton collisions. For this, two
20// main operation methods for this task are foreseen:
21//
22// 1) (under development) it should act as an auxiliary task and provide a
23// calibrated estimator
24//
25// 2) "Debug mode" which will also create a ROOT TTree object with event
26// by event info potentially used for exploration / calibration. This
27// includes the following info:
28//
29// --- All VZERO Amplitudes (saved as Float_t)
30// --- (optional) time for each channel
31// --- (optional) time width for each channel
32// --- GetReferenceMultiplicity Estimator, |eta|<0.5
33// --- GetReferenceMultiplicity Estimator, |eta|<0.8
34// --- (if MC) True Multiplicity, |eta|<0.5
35// --- (if MC) True Multiplicity, 2.8 < eta < 5.8 (VZEROA region)
36// --- (if MC) True Multiplicity, -3.7 < eta <-1.7 (VZEROC region)
37// --- Run Number
38//
39// +-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
40
41class TTree;
42class TParticle;
43class TVector3;
44
45//class AliMCEventHandler;
46//class AliMCEvent;
47//class AliStack;
48
49class AliESDVertex;
50class AliAODVertex;
51class AliESDv0;
52class AliAODv0;
53
54#include <Riostream.h>
55#include "TList.h"
56#include "TH1.h"
57#include "TH2.h"
58#include "TH3.h"
59#include "TFile.h"
60#include "THnSparse.h"
61#include "TVector3.h"
62#include "TCanvas.h"
63#include "TMath.h"
64#include "TLegend.h"
65//#include "AliLog.h"
66
67#include "AliESDEvent.h"
68#include "AliAODEvent.h"
69#include "AliV0vertexer.h"
70#include "AliCascadeVertexer.h"
71#include "AliESDpid.h"
72#include "AliESDtrack.h"
73#include "AliESDtrackCuts.h"
74#include "AliInputEventHandler.h"
75#include "AliAnalysisManager.h"
76#include "AliMCEventHandler.h"
77#include "AliMCEvent.h"
78#include "AliStack.h"
79#include "AliCentrality.h"
535c1213 80#include "AliPPVsMultUtils.h"
2ca41fef 81#include "AliPWG0Helper.h"
af0d4791 82#include "AliCFContainer.h"
83#include "AliMultiplicity.h"
84#include "AliAODMCParticle.h"
85#include "AliESDcascade.h"
86#include "AliAODcascade.h"
87#include "AliESDUtils.h"
88#include "AliGenEventHeader.h"
89#include "AliAnalysisTaskSE.h"
90#include "AliAnalysisUtils.h"
91#include "AliAnalysisTaskStrangenessVsMultiplicityMC.h"
92
535c1213 93
af0d4791 94using std::cout;
95using std::endl;
96
97ClassImp(AliAnalysisTaskStrangenessVsMultiplicityMC)
98
99AliAnalysisTaskStrangenessVsMultiplicityMC::AliAnalysisTaskStrangenessVsMultiplicityMC()
535c1213 100 : AliAnalysisTaskSE(), fListHist(0), fTreeEvent(0), fTreeV0(0), fTreeCascade(0), fPIDResponse(0), fESDtrackCuts(0), fPPVsMultUtils(0),
af0d4791 101 fkSaveV0Tree ( kFALSE ),
102 fkSaveCascadeTree ( kTRUE ),
dedb9266 103 fkRunVertexers ( kTRUE ),
104 fkSkipEventSelection( kFALSE ),
af0d4791 105 //---> Variables for fTreeEvent
106 fAmplitude_V0A (0),
107 fAmplitude_V0C (0),
108 fAmplitude_V0AEq (0),
109 fAmplitude_V0CEq (0),
110 fCentrality_V0A(0),
111 fCentrality_V0C(0),
112 fCentrality_V0M(0),
113 fCentrality_V0AEq(0),
114 fCentrality_V0CEq(0),
115 fCentrality_V0MEq(0),
535c1213 116 fCustomCentrality_V0M(0),
117 fCustomCentrality_V0MEq(0),
af0d4791 118 fRefMultEta5(0),
119 fRefMultEta8(0),
120 fTrueMultEta5(0),
121 fTrueMultEta8(0),
122 fTrueMultVZEROA(0),
123 fTrueMultVZEROC(0),
124 fRunNumber(0),
dedb9266 125 fEvSel_HasAtLeastSPDVertex(0),
126 fEvSel_VtxZCut(0),
127 fEvSel_IsNotPileup(0),
128 fEvSel_IsNotPileupMV(0),
129 fEvSel_IsNotPileupInMultBins(0),
130 fEvSel_HasVtxContributor(0),
131 fEvSel_Triggered(0),
2ca41fef 132 fEvSel_VtxZ(0),
133 fEvSel_MCType(0),
161c57ac 134 //---> Variables for fTreeV0
135 fTreeVariableChi2V0(0),
136 fTreeVariableDcaV0Daughters(0),
137 fTreeVariableDcaV0ToPrimVertex(0),
138 fTreeVariableDcaPosToPrimVertex(0),
139 fTreeVariableDcaNegToPrimVertex(0),
140 fTreeVariableV0CosineOfPointingAngle(0),
141 fTreeVariableV0Radius(0),
142 fTreeVariablePt(0),
143 fTreeVariablePtMC(0),
144 fTreeVariableRapK0Short(0),
145 fTreeVariableRapLambda(0),
146 fTreeVariableRapMC(0),
147 fTreeVariableInvMassK0s(0),
148 fTreeVariableInvMassLambda(0),
149 fTreeVariableInvMassAntiLambda(0),
150 fTreeVariableAlphaV0(0),
151 fTreeVariablePtArmV0(0),
152 fTreeVariableNegEta(0),
153 fTreeVariablePosEta(0),
154
155 fTreeVariableNSigmasPosProton(0),
156 fTreeVariableNSigmasPosPion(0),
157 fTreeVariableNSigmasNegProton(0),
158 fTreeVariableNSigmasNegPion(0),
159
160 fTreeVariableDistOverTotMom(0),
161 fTreeVariableLeastNbrCrossedRows(0),
162 fTreeVariableLeastRatioCrossedRowsOverFindable(0),
163
164 fTreeVariableCentV0A(0),
165 fTreeVariableCentV0C(0),
166 fTreeVariableCentV0M(0),
167 fTreeVariableCentV0AEq(0),
168 fTreeVariableCentV0CEq(0),
169 fTreeVariableCentV0MEq(0),
9a28e152 170 fTreeVariableCustomCentV0M(0),
161c57ac 171 fTreeVariableAmpV0A(0),
172 fTreeVariableAmpV0C(0),
173 fTreeVariableAmpV0AEq(0),
174 fTreeVariableAmpV0CEq(0),
175 fTreeVariableRefMultEta8(0),
176 fTreeVariableRefMultEta5(0),
177 fTreeVariableRunNumber(0),
178
179 fTreeVariablePtMother(0),
180 fTreeVariablePID(0),
181 fTreeVariablePIDPositive(0),
182 fTreeVariablePIDNegative(0),
183 fTreeVariablePIDMother(0),
184 fTreeVariablePrimaryStatus(0),
185 fTreeVariablePrimaryStatusMother(0),
af0d4791 186 //---> Variables for fTreeCascade
187 fTreeCascVarCharge(0),
188 fTreeCascVarMassAsXi(0),
189 fTreeCascVarMassAsOmega(0),
190 fTreeCascVarPt(0),
191 fTreeCascVarPtMC(0),
192 fTreeCascVarRapXi(0),
193 fTreeCascVarRapOmega(0),
194 fTreeCascVarRapMC(0),
195 fTreeCascVarNegEta(0),
196 fTreeCascVarPosEta(0),
197 fTreeCascVarBachEta(0),
198 fTreeCascVarDCACascDaughters(0),
199 fTreeCascVarDCABachToPrimVtx(0),
200 fTreeCascVarDCAV0Daughters(0),
201 fTreeCascVarDCAV0ToPrimVtx(0),
202 fTreeCascVarDCAPosToPrimVtx(0),
203 fTreeCascVarDCANegToPrimVtx(0),
204 fTreeCascVarCascCosPointingAngle(0),
205 fTreeCascVarCascRadius(0),
206 fTreeCascVarV0Mass(0),
207 fTreeCascVarV0CosPointingAngle(0),
208 fTreeCascVarV0CosPointingAngleSpecial(0),
209 fTreeCascVarV0Radius(0),
210 fTreeCascVarLeastNbrClusters(0),
211 fTreeCascVarDistOverTotMom(0),
212 fTreeCascVarNegNSigmaPion(0),
213 fTreeCascVarNegNSigmaProton(0),
214 fTreeCascVarPosNSigmaPion(0),
215 fTreeCascVarPosNSigmaProton(0),
216 fTreeCascVarBachNSigmaPion(0),
217 fTreeCascVarBachNSigmaKaon(0),
218 fTreeCascVarCentV0A(0),
219 fTreeCascVarCentV0C(0),
220 fTreeCascVarCentV0M(0),
221 fTreeCascVarCentV0AEq(0),
222 fTreeCascVarCentV0CEq(0),
223 fTreeCascVarCentV0MEq(0),
9a28e152 224 fTreeCascVarCustomCentV0M(0),
af0d4791 225 fTreeCascVarAmpV0A(0),
226 fTreeCascVarAmpV0C(0),
227 fTreeCascVarAmpV0AEq(0),
228 fTreeCascVarAmpV0CEq(0),
229 fTreeCascVarRefMultEta8(0),
230 fTreeCascVarRefMultEta5(0),
231 fTreeCascVarTrueMultEta5(0),
232 fTreeCascVarTrueMultEta8(0),
233 fTreeCascVarTrueMultVZEROA(0),
234 fTreeCascVarTrueMultVZEROC(0),
235 fTreeCascVarIsPhysicalPrimary(0),
236 fTreeCascVarPID(0),
26d39641 237 fTreeCascVarRunNumber(0),
af0d4791 238 //---> Histograms
239 fHistEventCounter(0),
240 //---> MC Generated Histo (analysis level)
161c57ac 241 fHistPt_GenK0Short(0),
242 fHistPt_GenLambda(0),
243 fHistPt_GenAntiLambda(0),
af0d4791 244 fHistPt_GenXiMinus(0),
245 fHistPt_GenXiPlus(0),
246 fHistPt_GenOmegaMinus(0),
247 fHistPt_GenOmegaPlus(0),
248
249 //VsRefMult
250 fHistPtVsRefMultEta5_GenXiMinus(0),
251 fHistPtVsRefMultEta5_GenXiPlus(0),
252 fHistPtVsRefMultEta5_GenOmegaMinus(0),
253 fHistPtVsRefMultEta5_GenOmegaPlus(0),
254 fHistPtVsRefMultEta8_GenXiMinus(0),
255 fHistPtVsRefMultEta8_GenXiPlus(0),
256 fHistPtVsRefMultEta8_GenOmegaMinus(0),
257 fHistPtVsRefMultEta8_GenOmegaPlus(0),
258
259 //VsCentralities
260 fHistPtVsCentV0A_GenXiMinus(0),
261 fHistPtVsCentV0A_GenXiPlus(0),
262 fHistPtVsCentV0A_GenOmegaMinus(0),
263 fHistPtVsCentV0A_GenOmegaPlus(0),
264 fHistPtVsCentV0C_GenXiMinus(0),
265 fHistPtVsCentV0C_GenXiPlus(0),
266 fHistPtVsCentV0C_GenOmegaMinus(0),
267 fHistPtVsCentV0C_GenOmegaPlus(0),
268 fHistPtVsCentV0M_GenXiMinus(0),
269 fHistPtVsCentV0M_GenXiPlus(0),
270 fHistPtVsCentV0M_GenOmegaMinus(0),
271 fHistPtVsCentV0M_GenOmegaPlus(0),
272
273 //Equalized
274 fHistPtVsCentV0AEq_GenXiMinus(0),
275 fHistPtVsCentV0AEq_GenXiPlus(0),
276 fHistPtVsCentV0AEq_GenOmegaMinus(0),
277 fHistPtVsCentV0AEq_GenOmegaPlus(0),
278 fHistPtVsCentV0CEq_GenXiMinus(0),
279 fHistPtVsCentV0CEq_GenXiPlus(0),
280 fHistPtVsCentV0CEq_GenOmegaMinus(0),
281 fHistPtVsCentV0CEq_GenOmegaPlus(0),
282 fHistPtVsCentV0MEq_GenXiMinus(0),
283 fHistPtVsCentV0MEq_GenXiPlus(0),
284 fHistPtVsCentV0MEq_GenOmegaMinus(0),
285 fHistPtVsCentV0MEq_GenOmegaPlus(0),
286
287 //VsAmp
288 fHistPtVsAmpV0A_GenXiMinus(0),
289 fHistPtVsAmpV0A_GenXiPlus(0),
290 fHistPtVsAmpV0A_GenOmegaMinus(0),
291 fHistPtVsAmpV0A_GenOmegaPlus(0),
292 fHistPtVsAmpV0C_GenXiMinus(0),
293 fHistPtVsAmpV0C_GenXiPlus(0),
294 fHistPtVsAmpV0C_GenOmegaMinus(0),
295 fHistPtVsAmpV0C_GenOmegaPlus(0),
296 fHistPtVsAmpV0M_GenXiMinus(0),
297 fHistPtVsAmpV0M_GenXiPlus(0),
298 fHistPtVsAmpV0M_GenOmegaMinus(0),
299 fHistPtVsAmpV0M_GenOmegaPlus(0),
300 //Equalized Amps
301 fHistPtVsAmpV0AEq_GenXiMinus(0),
302 fHistPtVsAmpV0AEq_GenXiPlus(0),
303 fHistPtVsAmpV0AEq_GenOmegaMinus(0),
304 fHistPtVsAmpV0AEq_GenOmegaPlus(0),
305 fHistPtVsAmpV0CEq_GenXiMinus(0),
306 fHistPtVsAmpV0CEq_GenXiPlus(0),
307 fHistPtVsAmpV0CEq_GenOmegaMinus(0),
308 fHistPtVsAmpV0CEq_GenOmegaPlus(0),
309 fHistPtVsAmpV0MEq_GenXiMinus(0),
310 fHistPtVsAmpV0MEq_GenXiPlus(0),
311 fHistPtVsAmpV0MEq_GenOmegaMinus(0),
1dac6d41 312 fHistPtVsAmpV0MEq_GenOmegaPlus(0),
313 fHistVZEROResponseStudy(0),
314 fHistVZEROResponseStudyTotal(0)
af0d4791 315
316//------------------------------------------------
317// Tree Variables
318{
319
320}
321
322AliAnalysisTaskStrangenessVsMultiplicityMC::AliAnalysisTaskStrangenessVsMultiplicityMC(const char *name)
535c1213 323 : AliAnalysisTaskSE(name), fListHist(0), fTreeEvent(0), fTreeV0(0), fTreeCascade(0), fPIDResponse(0), fESDtrackCuts(0), fPPVsMultUtils(0),
af0d4791 324 fkSaveV0Tree ( kFALSE ),
325 fkSaveCascadeTree ( kTRUE ),
326 fkRunVertexers ( kTRUE ),
dedb9266 327 fkSkipEventSelection( kFALSE ),
af0d4791 328 //---> Variables for fTreeEvent
329 fAmplitude_V0A (0),
330 fAmplitude_V0C (0),
331 fAmplitude_V0AEq (0),
332 fAmplitude_V0CEq (0),
333 fCentrality_V0A(0),
334 fCentrality_V0C(0),
335 fCentrality_V0M(0),
336 fCentrality_V0AEq(0),
337 fCentrality_V0CEq(0),
338 fCentrality_V0MEq(0),
535c1213 339 fCustomCentrality_V0M(0),
340 fCustomCentrality_V0MEq(0),
af0d4791 341 fRefMultEta5(0),
342 fRefMultEta8(0),
343 fTrueMultEta5(0),
344 fTrueMultEta8(0),
345 fTrueMultVZEROA(0),
346 fTrueMultVZEROC(0),
347 fRunNumber(0),
dedb9266 348 fEvSel_HasAtLeastSPDVertex(0),
349 fEvSel_VtxZCut(0),
350 fEvSel_IsNotPileup(0),
351 fEvSel_IsNotPileupMV(0),
352 fEvSel_IsNotPileupInMultBins(0),
353 fEvSel_HasVtxContributor(0),
354 fEvSel_Triggered(0),
2ca41fef 355 fEvSel_VtxZ(0),
356 fEvSel_MCType(0),
161c57ac 357 //---> Variables for fTreeV0
358 fTreeVariableChi2V0(0),
359 fTreeVariableDcaV0Daughters(0),
360 fTreeVariableDcaV0ToPrimVertex(0),
361 fTreeVariableDcaPosToPrimVertex(0),
362 fTreeVariableDcaNegToPrimVertex(0),
363 fTreeVariableV0CosineOfPointingAngle(0),
364 fTreeVariableV0Radius(0),
365 fTreeVariablePt(0),
366 fTreeVariablePtMC(0),
367 fTreeVariableRapK0Short(0),
368 fTreeVariableRapLambda(0),
369 fTreeVariableRapMC(0),
370 fTreeVariableInvMassK0s(0),
371 fTreeVariableInvMassLambda(0),
372 fTreeVariableInvMassAntiLambda(0),
373 fTreeVariableAlphaV0(0),
374 fTreeVariablePtArmV0(0),
375 fTreeVariableNegEta(0),
376 fTreeVariablePosEta(0),
377
378 fTreeVariableNSigmasPosProton(0),
379 fTreeVariableNSigmasPosPion(0),
380 fTreeVariableNSigmasNegProton(0),
381 fTreeVariableNSigmasNegPion(0),
382
383 fTreeVariableDistOverTotMom(0),
384 fTreeVariableLeastNbrCrossedRows(0),
385 fTreeVariableLeastRatioCrossedRowsOverFindable(0),
386
387 fTreeVariableCentV0A(0),
388 fTreeVariableCentV0C(0),
389 fTreeVariableCentV0M(0),
390 fTreeVariableCentV0AEq(0),
391 fTreeVariableCentV0CEq(0),
392 fTreeVariableCentV0MEq(0),
9a28e152 393 fTreeVariableCustomCentV0M(0),
161c57ac 394 fTreeVariableAmpV0A(0),
395 fTreeVariableAmpV0C(0),
396 fTreeVariableAmpV0AEq(0),
397 fTreeVariableAmpV0CEq(0),
398 fTreeVariableRefMultEta8(0),
399 fTreeVariableRefMultEta5(0),
400 fTreeVariableRunNumber(0),
401
402 fTreeVariablePtMother(0),
403 fTreeVariablePID(0),
404 fTreeVariablePIDPositive(0),
405 fTreeVariablePIDNegative(0),
406 fTreeVariablePIDMother(0),
407 fTreeVariablePrimaryStatus(0),
408 fTreeVariablePrimaryStatusMother(0),
af0d4791 409 //---> Variables for fTreeCascade
410 fTreeCascVarCharge(0),
411 fTreeCascVarMassAsXi(0),
412 fTreeCascVarMassAsOmega(0),
413 fTreeCascVarPt(0),
414 fTreeCascVarPtMC(0),
415 fTreeCascVarRapXi(0),
416 fTreeCascVarRapOmega(0),
417 fTreeCascVarRapMC(0),
418 fTreeCascVarNegEta(0),
419 fTreeCascVarPosEta(0),
420 fTreeCascVarBachEta(0),
421 fTreeCascVarDCACascDaughters(0),
422 fTreeCascVarDCABachToPrimVtx(0),
423 fTreeCascVarDCAV0Daughters(0),
424 fTreeCascVarDCAV0ToPrimVtx(0),
425 fTreeCascVarDCAPosToPrimVtx(0),
426 fTreeCascVarDCANegToPrimVtx(0),
427 fTreeCascVarCascCosPointingAngle(0),
428 fTreeCascVarCascRadius(0),
429 fTreeCascVarV0Mass(0),
430 fTreeCascVarV0CosPointingAngle(0),
431 fTreeCascVarV0CosPointingAngleSpecial(0),
432 fTreeCascVarV0Radius(0),
433 fTreeCascVarLeastNbrClusters(0),
434 fTreeCascVarDistOverTotMom(0),
435 fTreeCascVarNegNSigmaPion(0),
436 fTreeCascVarNegNSigmaProton(0),
437 fTreeCascVarPosNSigmaPion(0),
438 fTreeCascVarPosNSigmaProton(0),
439 fTreeCascVarBachNSigmaPion(0),
440 fTreeCascVarBachNSigmaKaon(0),
441 fTreeCascVarCentV0A(0),
442 fTreeCascVarCentV0C(0),
443 fTreeCascVarCentV0M(0),
444 fTreeCascVarCentV0AEq(0),
445 fTreeCascVarCentV0CEq(0),
446 fTreeCascVarCentV0MEq(0),
9a28e152 447 fTreeCascVarCustomCentV0M(0),
af0d4791 448 fTreeCascVarAmpV0A(0),
449 fTreeCascVarAmpV0C(0),
450 fTreeCascVarAmpV0AEq(0),
451 fTreeCascVarAmpV0CEq(0),
452 fTreeCascVarRefMultEta8(0),
453 fTreeCascVarRefMultEta5(0),
454 fTreeCascVarTrueMultEta5(0),
455 fTreeCascVarTrueMultEta8(0),
456 fTreeCascVarTrueMultVZEROA(0),
457 fTreeCascVarTrueMultVZEROC(0),
458 fTreeCascVarIsPhysicalPrimary(0),
459 fTreeCascVarPID(0),
26d39641 460 fTreeCascVarRunNumber(0),
af0d4791 461 //---> Histograms
462 fHistEventCounter(0),
463 //---> MC Generated Histo (analysis level)
161c57ac 464 fHistPt_GenK0Short(0),
465 fHistPt_GenLambda(0),
466 fHistPt_GenAntiLambda(0),
af0d4791 467 fHistPt_GenXiMinus(0),
468 fHistPt_GenXiPlus(0),
469 fHistPt_GenOmegaMinus(0),
470 fHistPt_GenOmegaPlus(0),
471
472 //VsRefMult
473 fHistPtVsRefMultEta5_GenXiMinus(0),
474 fHistPtVsRefMultEta5_GenXiPlus(0),
475 fHistPtVsRefMultEta5_GenOmegaMinus(0),
476 fHistPtVsRefMultEta5_GenOmegaPlus(0),
477 fHistPtVsRefMultEta8_GenXiMinus(0),
478 fHistPtVsRefMultEta8_GenXiPlus(0),
479 fHistPtVsRefMultEta8_GenOmegaMinus(0),
480 fHistPtVsRefMultEta8_GenOmegaPlus(0),
481
482 //VsCentralities
483 fHistPtVsCentV0A_GenXiMinus(0),
484 fHistPtVsCentV0A_GenXiPlus(0),
485 fHistPtVsCentV0A_GenOmegaMinus(0),
486 fHistPtVsCentV0A_GenOmegaPlus(0),
487 fHistPtVsCentV0C_GenXiMinus(0),
488 fHistPtVsCentV0C_GenXiPlus(0),
489 fHistPtVsCentV0C_GenOmegaMinus(0),
490 fHistPtVsCentV0C_GenOmegaPlus(0),
491 fHistPtVsCentV0M_GenXiMinus(0),
492 fHistPtVsCentV0M_GenXiPlus(0),
493 fHistPtVsCentV0M_GenOmegaMinus(0),
494 fHistPtVsCentV0M_GenOmegaPlus(0),
495
496 //Equalized
497 fHistPtVsCentV0AEq_GenXiMinus(0),
498 fHistPtVsCentV0AEq_GenXiPlus(0),
499 fHistPtVsCentV0AEq_GenOmegaMinus(0),
500 fHistPtVsCentV0AEq_GenOmegaPlus(0),
501 fHistPtVsCentV0CEq_GenXiMinus(0),
502 fHistPtVsCentV0CEq_GenXiPlus(0),
503 fHistPtVsCentV0CEq_GenOmegaMinus(0),
504 fHistPtVsCentV0CEq_GenOmegaPlus(0),
505 fHistPtVsCentV0MEq_GenXiMinus(0),
506 fHistPtVsCentV0MEq_GenXiPlus(0),
507 fHistPtVsCentV0MEq_GenOmegaMinus(0),
508 fHistPtVsCentV0MEq_GenOmegaPlus(0),
509
510 //VsAmp
511 fHistPtVsAmpV0A_GenXiMinus(0),
512 fHistPtVsAmpV0A_GenXiPlus(0),
513 fHistPtVsAmpV0A_GenOmegaMinus(0),
514 fHistPtVsAmpV0A_GenOmegaPlus(0),
515 fHistPtVsAmpV0C_GenXiMinus(0),
516 fHistPtVsAmpV0C_GenXiPlus(0),
517 fHistPtVsAmpV0C_GenOmegaMinus(0),
518 fHistPtVsAmpV0C_GenOmegaPlus(0),
519 fHistPtVsAmpV0M_GenXiMinus(0),
520 fHistPtVsAmpV0M_GenXiPlus(0),
521 fHistPtVsAmpV0M_GenOmegaMinus(0),
522 fHistPtVsAmpV0M_GenOmegaPlus(0),
523 //Equalized Amps
524 fHistPtVsAmpV0AEq_GenXiMinus(0),
525 fHistPtVsAmpV0AEq_GenXiPlus(0),
526 fHistPtVsAmpV0AEq_GenOmegaMinus(0),
527 fHistPtVsAmpV0AEq_GenOmegaPlus(0),
528 fHistPtVsAmpV0CEq_GenXiMinus(0),
529 fHistPtVsAmpV0CEq_GenXiPlus(0),
530 fHistPtVsAmpV0CEq_GenOmegaMinus(0),
531 fHistPtVsAmpV0CEq_GenOmegaPlus(0),
532 fHistPtVsAmpV0MEq_GenXiMinus(0),
533 fHistPtVsAmpV0MEq_GenXiPlus(0),
534 fHistPtVsAmpV0MEq_GenOmegaMinus(0),
1dac6d41 535 fHistPtVsAmpV0MEq_GenOmegaPlus(0),
536 fHistVZEROResponseStudy(0),
537 fHistVZEROResponseStudyTotal(0)
af0d4791 538{
539
540 //Re-vertex: Will only apply for cascade candidates
541
542 fV0VertexerSels[0] = 33. ; // max allowed chi2
543 fV0VertexerSels[1] = 0.02; // min allowed impact parameter for the 1st daughter (LHC09a4 : 0.05)
544 fV0VertexerSels[2] = 0.02; // min allowed impact parameter for the 2nd daughter (LHC09a4 : 0.05)
76ebc06c 545 fV0VertexerSels[3] = 2.0 ; // max allowed DCA between the daughter tracks (LHC09a4 : 0.5)
af0d4791 546 fV0VertexerSels[4] = 0.95; // min allowed cosine of V0's pointing angle (LHC09a4 : 0.99)
547 fV0VertexerSels[5] = 1.0 ; // min radius of the fiducial volume (LHC09a4 : 0.2)
548 fV0VertexerSels[6] = 200. ; // max radius of the fiducial volume (LHC09a4 : 100.0)
549
550 fCascadeVertexerSels[0] = 33. ; // max allowed chi2 (same as PDC07)
551 fCascadeVertexerSels[1] = 0.05 ; // min allowed V0 impact parameter (PDC07 : 0.05 / LHC09a4 : 0.025 )
552 fCascadeVertexerSels[2] = 0.010; // "window" around the Lambda mass (PDC07 : 0.008 / LHC09a4 : 0.010 )
553 fCascadeVertexerSels[3] = 0.03 ; // min allowed bachelor's impact parameter (PDC07 : 0.035 / LHC09a4 : 0.025 )
554 fCascadeVertexerSels[4] = 2.0 ; // max allowed DCA between the V0 and the bachelor (PDC07 : 0.1 / LHC09a4 : 0.2 )
555 fCascadeVertexerSels[5] = 0.95 ; // min allowed cosine of the cascade pointing angle (PDC07 : 0.9985 / LHC09a4 : 0.998 )
556 fCascadeVertexerSels[6] = 0.4 ; // min radius of the fiducial volume (PDC07 : 0.9 / LHC09a4 : 0.2 )
557 fCascadeVertexerSels[7] = 100. ; // max radius of the fiducial volume (PDC07 : 100 / LHC09a4 : 100 )
558
559
560 DefineOutput(1, TList::Class()); // Event Counter Histo
561 DefineOutput(2, TTree::Class()); // Event Tree
562 DefineOutput(3, TTree::Class()); // V0 Tree
563 DefineOutput(4, TTree::Class()); // Cascade Tree
564}
565
566
567AliAnalysisTaskStrangenessVsMultiplicityMC::~AliAnalysisTaskStrangenessVsMultiplicityMC()
568{
569//------------------------------------------------
570// DESTRUCTOR
571//------------------------------------------------
572
573 if (fListHist){
574 delete fListHist;
575 fListHist = 0x0;
576 }
577 if (fTreeEvent){
578 delete fTreeEvent;
579 fTreeEvent = 0x0;
580 }
581 if (fTreeV0){
582 delete fTreeV0;
583 fTreeV0 = 0x0;
584 }
585 if (fTreeCascade){
586 delete fTreeCascade;
587 fTreeCascade = 0x0;
588 }
535c1213 589 if (fPPVsMultUtils){
590 delete fPPVsMultUtils;
591 fPPVsMultUtils = 0x0;
592 }
af0d4791 593}
594
595//________________________________________________________________________
596void AliAnalysisTaskStrangenessVsMultiplicityMC::UserCreateOutputObjects()
597{
598
599 OpenFile(2);
600 // Called once
601
602//------------------------------------------------
603
604 fTreeEvent = new TTree("fTreeEvent","Event");
605
606//------------------------------------------------
607// fTree Branch definitions - Event by Event info
608//------------------------------------------------
609
610//-----------BASIC-INFO---------------------------
611
612 //--- VZERO Data (Integrated)
613 fTreeEvent->Branch("fAmplitude_V0A",&fAmplitude_V0A,"fAmplitude_V0A/F");
614 fTreeEvent->Branch("fAmplitude_V0C",&fAmplitude_V0C,"fAmplitude_V0C/F");
615 fTreeEvent->Branch("fAmplitude_V0AEq",&fAmplitude_V0AEq,"fAmplitude_V0AEq/F");
616 fTreeEvent->Branch("fAmplitude_V0CEq",&fAmplitude_V0CEq,"fAmplitude_V0CEq/F");
617
618 //Info from AliCentrality (not necessarily 'centrality' per se)
619 fTreeEvent->Branch("fCentrality_V0A",&fCentrality_V0A,"fCentrality_V0A/F");
620 fTreeEvent->Branch("fCentrality_V0C",&fCentrality_V0C,"fCentrality_V0C/F");
34442f2a 621 fTreeEvent->Branch("fCentrality_V0M",&fCentrality_V0M,"fCentrality_V0M/F");
af0d4791 622 fTreeEvent->Branch("fCentrality_V0AEq",&fCentrality_V0AEq,"fCentrality_V0AEq/F");
623 fTreeEvent->Branch("fCentrality_V0CEq",&fCentrality_V0CEq,"fCentrality_V0CEq/F");
34442f2a 624 fTreeEvent->Branch("fCentrality_V0MEq",&fCentrality_V0MEq,"fCentrality_V0MEq/F");
535c1213 625
626 fTreeEvent->Branch("fCustomCentrality_V0M",&fCustomCentrality_V0M,"fCustomCentrality_V0M/F");
c2026b92 627 fTreeEvent->Branch("fCustomCentrality_V0MEq",&fCustomCentrality_V0MEq,"fCustomCentrality_V0MEq/F");
af0d4791 628
629 //Official GetReferenceMultiplicity
630 fTreeEvent->Branch("fRefMultEta5",&fRefMultEta5,"fRefMultEta5/I");
631 fTreeEvent->Branch("fRefMultEta8",&fRefMultEta8,"fRefMultEta8/I");
632
633 fTreeEvent->Branch("fTrueMultEta5",&fTrueMultEta5,"fTrueMultEta5/I");
634 fTreeEvent->Branch("fTrueMultEta8",&fTrueMultEta8,"fTrueMultEta8/I");
635 fTreeEvent->Branch("fTrueMultVZEROA",&fTrueMultVZEROA,"fTrueMultVZEROA/I");
636 fTreeEvent->Branch("fTrueMultVZEROC",&fTrueMultVZEROC,"fTrueMultVZEROC/I");
637
638 //Run Number
639 fTreeEvent->Branch("fRunNumber", &fRunNumber, "fRunNumber/I");
dedb9266 640
641 //Booleans for Event Selection
642 fTreeEvent->Branch("fEvSel_HasAtLeastSPDVertex", &fEvSel_HasAtLeastSPDVertex, "fEvSel_HasAtLeastSPDVertex/O");
643 fTreeEvent->Branch("fEvSel_VtxZCut", &fEvSel_VtxZCut, "fEvSel_VtxZCut/O");
644 fTreeEvent->Branch("fEvSel_IsNotPileup", &fEvSel_IsNotPileup, "fEvSel_IsNotPileup/O");
645 fTreeEvent->Branch("fEvSel_IsNotPileupMV", &fEvSel_IsNotPileupMV, "fEvSel_IsNotPileupMV/O");
646 fTreeEvent->Branch("fEvSel_IsNotPileupInMultBins", &fEvSel_IsNotPileupInMultBins, "fEvSel_IsNotPileupInMultBins/O");
647 fTreeEvent->Branch("fEvSel_HasVtxContributor", &fEvSel_HasVtxContributor, "fEvSel_HasVtxContributor/O");
648 fTreeEvent->Branch("fEvSel_Triggered", &fEvSel_Triggered, "fEvSel_Triggered/O");
2ca41fef 649 fTreeEvent->Branch("fEvSel_VtxZ", &fEvSel_VtxZ, "fEvSel_VtxZ/F");
650 fTreeEvent->Branch("fEvSel_MCType", &fEvSel_MCType, "fEvSel_MCType/I");
dedb9266 651
af0d4791 652
653 //Create Basic V0 Output Tree
654 fTreeV0 = new TTree( "fTreeV0", "V0 Candidates");
655
161c57ac 656//------------------------------------------------
657// fTreeV0 Branch definitions
658//------------------------------------------------
659
660//-----------BASIC-INFO---------------------------
661 fTreeV0->Branch("fTreeVariableChi2V0",&fTreeVariableChi2V0,"fTreeVariableChi2V0/F");
662 fTreeV0->Branch("fTreeVariableDcaV0Daughters",&fTreeVariableDcaV0Daughters,"fTreeVariableDcaV0Daughters/F");
663 fTreeV0->Branch("fTreeVariableDcaV0ToPrimVertex",&fTreeVariableDcaV0ToPrimVertex,"fTreeVariableDcaV0ToPrimVertex/F");
664 fTreeV0->Branch("fTreeVariableDcaPosToPrimVertex",&fTreeVariableDcaPosToPrimVertex,"fTreeVariableDcaPosToPrimVertex/F");
665 fTreeV0->Branch("fTreeVariableDcaNegToPrimVertex",&fTreeVariableDcaNegToPrimVertex,"fTreeVariableDcaNegToPrimVertex/F");
666 fTreeV0->Branch("fTreeVariableV0Radius",&fTreeVariableV0Radius,"fTreeVariableV0Radius/F");
667 fTreeV0->Branch("fTreeVariablePt",&fTreeVariablePt,"fTreeVariablePt/F");
668 fTreeV0->Branch("fTreeVariablePtMC",&fTreeVariablePtMC,"fTreeVariablePtMC/F");
669 fTreeV0->Branch("fTreeVariableRapK0Short",&fTreeVariableRapK0Short,"fTreeVariableRapK0Short/F");
670 fTreeV0->Branch("fTreeVariableRapLambda",&fTreeVariableRapLambda,"fTreeVariableRapLambda/F");
671 fTreeV0->Branch("fTreeVariableRapMC",&fTreeVariableRapMC,"fTreeVariableRapMC/F");
672 fTreeV0->Branch("fTreeVariableInvMassK0s",&fTreeVariableInvMassK0s,"fTreeVariableInvMassK0s/F");
673 fTreeV0->Branch("fTreeVariableInvMassLambda",&fTreeVariableInvMassLambda,"fTreeVariableInvMassLambda/F");
674 fTreeV0->Branch("fTreeVariableInvMassAntiLambda",&fTreeVariableInvMassAntiLambda,"fTreeVariableInvMassAntiLambda/F");
675 fTreeV0->Branch("fTreeVariableV0CosineOfPointingAngle",&fTreeVariableV0CosineOfPointingAngle,"fTreeVariableV0CosineOfPointingAngle/F");
676 fTreeV0->Branch("fTreeVariableAlphaV0",&fTreeVariableAlphaV0,"fTreeVariableAlphaV0/F");
677 fTreeV0->Branch("fTreeVariablePtArmV0",&fTreeVariablePtArmV0,"fTreeVariablePtArmV0/F");
678 fTreeV0->Branch("fTreeVariableLeastNbrCrossedRows",&fTreeVariableLeastNbrCrossedRows,"fTreeVariableLeastNbrCrossedRows/I");
679 fTreeV0->Branch("fTreeVariableLeastRatioCrossedRowsOverFindable",&fTreeVariableLeastRatioCrossedRowsOverFindable,"fTreeVariableLeastRatioCrossedRowsOverFindable/F");
680 fTreeV0->Branch("fTreeVariableDistOverTotMom",&fTreeVariableDistOverTotMom,"fTreeVariableDistOverTotMom/F");
681 fTreeV0->Branch("fTreeVariableNSigmasPosProton",&fTreeVariableNSigmasPosProton,"fTreeVariableNSigmasPosProton/F");
682 fTreeV0->Branch("fTreeVariableNSigmasPosPion",&fTreeVariableNSigmasPosPion,"fTreeVariableNSigmasPosPion/F");
683 fTreeV0->Branch("fTreeVariableNSigmasNegProton",&fTreeVariableNSigmasNegProton,"fTreeVariableNSigmasNegProton/F");
684 fTreeV0->Branch("fTreeVariableNSigmasNegPion",&fTreeVariableNSigmasNegPion,"fTreeVariableNSigmasNegPion/F");
685 fTreeV0->Branch("fTreeVariableNegEta",&fTreeVariableNegEta,"fTreeVariableNegEta/F");
686 fTreeV0->Branch("fTreeVariablePosEta",&fTreeVariablePosEta,"fTreeVariablePosEta/F");
687//-----------MULTIPLICITY-INFO--------------------
688 fTreeV0->Branch("fTreeVariableCentV0A",&fTreeVariableCentV0A,"fTreeVariableCentV0A/F");
689 fTreeV0->Branch("fTreeVariableCentV0C",&fTreeVariableCentV0C,"fTreeVariableCentV0C/F");
690 fTreeV0->Branch("fTreeVariableCentV0M",&fTreeVariableCentV0M,"fTreeVariableCentV0M/F");
691 fTreeV0->Branch("fTreeVariableCentV0AEq",&fTreeVariableCentV0AEq,"fTreeVariableCentV0AEq/F");
692 fTreeV0->Branch("fTreeVariableCentV0CEq",&fTreeVariableCentV0CEq,"fTreeVariableCentV0CEq/F");
693 fTreeV0->Branch("fTreeVariableCentV0MEq",&fTreeVariableCentV0MEq,"fTreeVariableCentV0MEq/F");
9a28e152 694 fTreeV0->Branch("fTreeVariableCustomCentV0M",&fTreeVariableCustomCentV0M,"fTreeVariableCustomCentV0M/F");
161c57ac 695 fTreeV0->Branch("fTreeVariableAmpV0A",&fTreeVariableAmpV0A,"fTreeVariableAmpV0A/F");
696 fTreeV0->Branch("fTreeVariableAmpV0C",&fTreeVariableAmpV0C,"fTreeVariableAmpV0C/F");
697 fTreeV0->Branch("fTreeVariableAmpV0AEq",&fTreeVariableAmpV0AEq,"fTreeVariableAmpV0AEq/F");
698 fTreeV0->Branch("fTreeVariableAmpV0CEq",&fTreeVariableAmpV0CEq,"fTreeVariableAmpV0CEq/F");
699 fTreeV0->Branch("fTreeVariableRefMultEta8",&fTreeVariableRefMultEta8,"fTreeVariableRefMultEta8/I");
700 fTreeV0->Branch("fTreeVariableRefMultEta5",&fTreeVariableRefMultEta5,"fTreeVariableRefMultEta5/I");
701 fTreeV0->Branch("fTreeVariableRunNumber",&fTreeVariableRunNumber,"fTreeVariableRunNumber/I");
702//-----------MC Exclusive info--------------------
703 fTreeV0->Branch("fTreeVariablePtMother",&fTreeVariablePtMother,"fTreeVariablePtMother/F");
704 fTreeV0->Branch("fTreeVariablePID",&fTreeVariablePID,"fTreeVariablePID/I");
705 fTreeV0->Branch("fTreeVariablePIDPositive",&fTreeVariablePIDPositive,"fTreeVariablePIDPositive/I");
706 fTreeV0->Branch("fTreeVariablePIDNegative",&fTreeVariablePIDNegative,"fTreeVariablePIDNegative/I");
707 fTreeV0->Branch("fTreeVariablePIDMother",&fTreeVariablePIDMother,"fTreeVariablePIDMother/I");
708 fTreeV0->Branch("fTreeVariablePrimaryStatus",&fTreeVariablePrimaryStatus,"fTreeVariablePrimaryStatus/I");
709 fTreeV0->Branch("fTreeVariablePrimaryStatusMother",&fTreeVariablePrimaryStatusMother,"fTreeVariablePrimaryStatusMother/I");
710//------------------------------------------------
711
af0d4791 712 //Create Cascade output tree
713 fTreeCascade = new TTree("fTreeCascade","CascadeCandidates");
714
715//------------------------------------------------
716// fTreeCascade Branch definitions - Cascade Tree
717//------------------------------------------------
718
719//-----------BASIC-INFO---------------------------
720 fTreeCascade->Branch("fTreeCascVarCharge",&fTreeCascVarCharge,"fTreeCascVarCharge/I");
721 fTreeCascade->Branch("fTreeCascVarMassAsXi",&fTreeCascVarMassAsXi,"fTreeCascVarMassAsXi/F");
722 fTreeCascade->Branch("fTreeCascVarMassAsOmega",&fTreeCascVarMassAsOmega,"fTreeCascVarMassAsOmega/F");
723 fTreeCascade->Branch("fTreeCascVarPt",&fTreeCascVarPt,"fTreeCascVarPt/F");
724 fTreeCascade->Branch("fTreeCascVarPtMC",&fTreeCascVarPtMC,"fTreeCascVarPtMC/F");
725 fTreeCascade->Branch("fTreeCascVarRapXi",&fTreeCascVarRapXi,"fTreeCascVarRapXi/F");
726 fTreeCascade->Branch("fTreeCascVarRapOmega",&fTreeCascVarRapOmega,"fTreeCascVarRapOmega/F");
727 fTreeCascade->Branch("fTreeCascVarRapMC",&fTreeCascVarRapMC,"fTreeCascVarRapMC/F");
728 fTreeCascade->Branch("fTreeCascVarNegEta",&fTreeCascVarNegEta,"fTreeCascVarNegEta/F");
729 fTreeCascade->Branch("fTreeCascVarPosEta",&fTreeCascVarPosEta,"fTreeCascVarPosEta/F");
730 fTreeCascade->Branch("fTreeCascVarBachEta",&fTreeCascVarBachEta,"fTreeCascVarBachEta/F");
731//-----------INFO-FOR-CUTS------------------------
732 fTreeCascade->Branch("fTreeCascVarDCACascDaughters",&fTreeCascVarDCACascDaughters,"fTreeCascVarDCACascDaughters/F");
733 fTreeCascade->Branch("fTreeCascVarDCABachToPrimVtx",&fTreeCascVarDCABachToPrimVtx,"fTreeCascVarDCABachToPrimVtx/F");
734 fTreeCascade->Branch("fTreeCascVarDCAV0Daughters",&fTreeCascVarDCAV0Daughters,"fTreeCascVarDCAV0Daughters/F");
735 fTreeCascade->Branch("fTreeCascVarDCAV0ToPrimVtx",&fTreeCascVarDCAV0ToPrimVtx,"fTreeCascVarDCAV0ToPrimVtx/F");
736 fTreeCascade->Branch("fTreeCascVarDCAPosToPrimVtx",&fTreeCascVarDCAPosToPrimVtx,"fTreeCascVarDCAPosToPrimVtx/F");
737 fTreeCascade->Branch("fTreeCascVarDCANegToPrimVtx",&fTreeCascVarDCANegToPrimVtx,"fTreeCascVarDCANegToPrimVtx/F");
738 fTreeCascade->Branch("fTreeCascVarCascCosPointingAngle",&fTreeCascVarCascCosPointingAngle,"fTreeCascVarCascCosPointingAngle/F");
739 fTreeCascade->Branch("fTreeCascVarCascRadius",&fTreeCascVarCascRadius,"fTreeCascVarCascRadius/F");
740 fTreeCascade->Branch("fTreeCascVarV0Mass",&fTreeCascVarV0Mass,"fTreeCascVarV0Mass/F");
741 fTreeCascade->Branch("fTreeCascVarV0CosPointingAngle",&fTreeCascVarV0CosPointingAngle,"fTreeCascVarV0CosPointingAngle/F");
742 fTreeCascade->Branch("fTreeCascVarV0CosPointingAngleSpecial",&fTreeCascVarV0CosPointingAngleSpecial,"fTreeCascVarV0CosPointingAngleSpecial/F");
743 fTreeCascade->Branch("fTreeCascVarV0Radius",&fTreeCascVarV0Radius,"fTreeCascVarV0Radius/F");
744 fTreeCascade->Branch("fTreeCascVarLeastNbrClusters",&fTreeCascVarLeastNbrClusters,"fTreeCascVarLeastNbrClusters/I");
745//-----------MULTIPLICITY-INFO--------------------
746 fTreeCascade->Branch("fTreeCascVarCentV0A",&fTreeCascVarCentV0A,"fTreeCascVarCentV0A/F");
747 fTreeCascade->Branch("fTreeCascVarCentV0C",&fTreeCascVarCentV0C,"fTreeCascVarCentV0C/F");
748 fTreeCascade->Branch("fTreeCascVarCentV0M",&fTreeCascVarCentV0M,"fTreeCascVarCentV0M/F");
749 fTreeCascade->Branch("fTreeCascVarCentV0AEq",&fTreeCascVarCentV0AEq,"fTreeCascVarCentV0AEq/F");
750 fTreeCascade->Branch("fTreeCascVarCentV0CEq",&fTreeCascVarCentV0CEq,"fTreeCascVarCentV0CEq/F");
751 fTreeCascade->Branch("fTreeCascVarCentV0MEq",&fTreeCascVarCentV0MEq,"fTreeCascVarCentV0MEq/F");
9a28e152 752 fTreeCascade->Branch("fTreeCascVarCustomCentV0M",&fTreeCascVarCustomCentV0M,"fTreeCascVarCustomCentV0M/F");
af0d4791 753 fTreeCascade->Branch("fTreeCascVarAmpV0A",&fTreeCascVarAmpV0A,"fTreeCascVarAmpV0A/F");
754 fTreeCascade->Branch("fTreeCascVarAmpV0C",&fTreeCascVarAmpV0C,"fTreeCascVarAmpV0C/F");
755 fTreeCascade->Branch("fTreeCascVarAmpV0AEq",&fTreeCascVarAmpV0AEq,"fTreeCascVarAmpV0AEq/F");
756 fTreeCascade->Branch("fTreeCascVarAmpV0CEq",&fTreeCascVarAmpV0CEq,"fTreeCascVarAmpV0CEq/F");
757 fTreeCascade->Branch("fTreeCascVarRefMultEta8",&fTreeCascVarRefMultEta8,"fTreeCascVarRefMultEta8/I");
758 fTreeCascade->Branch("fTreeCascVarRefMultEta5",&fTreeCascVarRefMultEta5,"fTreeCascVarRefMultEta5/I");
759 fTreeCascade->Branch("fTreeCascVarTrueMultEta5",&fTreeCascVarTrueMultEta5,"fTreeCascVarTrueMultEta5/I");
760 fTreeCascade->Branch("fTreeCascVarTrueMultEta8",&fTreeCascVarTrueMultEta8,"fTreeCascVarTrueMultEta8/I");
761 fTreeCascade->Branch("fTreeCascVarTrueMultVZEROA",&fTreeCascVarTrueMultVZEROA,"fTreeCascVarTrueMultVZEROA/I");
762 fTreeCascade->Branch("fTreeCascVarTrueMultVZEROC",&fTreeCascVarTrueMultVZEROC,"fTreeCascVarTrueMultVZEROC/I");
763 fTreeCascade->Branch("fTreeCascVarIsPhysicalPrimary",&fTreeCascVarIsPhysicalPrimary,"fTreeCascVarIsPhysicalPrimary/I");
764 fTreeCascade->Branch("fTreeCascVarPID",&fTreeCascVarPID,"fTreeCascVarPID/I");
26d39641 765 fTreeCascade->Branch("fTreeCascVarRunNumber",&fTreeCascVarRunNumber,"fTreeCascVarRunNumber/I");
af0d4791 766//-----------DECAY-LENGTH-INFO--------------------
767 fTreeCascade->Branch("fTreeCascVarDistOverTotMom",&fTreeCascVarDistOverTotMom,"fTreeCascVarDistOverTotMom/F");
768//------------------------------------------------
769 fTreeCascade->Branch("fTreeCascVarNegNSigmaPion",&fTreeCascVarNegNSigmaPion,"fTreeCascVarNegNSigmaPion/F");
770 fTreeCascade->Branch("fTreeCascVarNegNSigmaProton",&fTreeCascVarNegNSigmaProton,"fTreeCascVarNegNSigmaProton/F");
771 fTreeCascade->Branch("fTreeCascVarPosNSigmaPion",&fTreeCascVarPosNSigmaPion,"fTreeCascVarPosNSigmaPion/F");
772 fTreeCascade->Branch("fTreeCascVarPosNSigmaProton",&fTreeCascVarPosNSigmaProton,"fTreeCascVarPosNSigmaProton/F");
773 fTreeCascade->Branch("fTreeCascVarBachNSigmaPion",&fTreeCascVarBachNSigmaPion,"fTreeCascVarBachNSigmaPion/F");
774 fTreeCascade->Branch("fTreeCascVarBachNSigmaKaon",&fTreeCascVarBachNSigmaKaon,"fTreeCascVarBachNSigmaKaon/F");
775
776//------------------------------------------------
777// Particle Identification Setup
778//------------------------------------------------
779
780 AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
781 AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
782 fPIDResponse = inputHandler->GetPIDResponse();
783
784 // Multiplicity
785 if(! fESDtrackCuts ){
786 fESDtrackCuts = new AliESDtrackCuts();
787 }
535c1213 788 if(! fPPVsMultUtils ){
789 fPPVsMultUtils = new AliPPVsMultUtils();
790 }
af0d4791 791
792//------------------------------------------------
793// V0 Multiplicity Histograms
794//------------------------------------------------
795
796 // Create histograms
797 OpenFile(1);
798 fListHist = new TList();
799 fListHist->SetOwner(); // See http://root.cern.ch/root/html/TCollection.html#TCollection:SetOwner
800
801 if(! fHistEventCounter ) {
802 //Histogram Output: Event-by-Event
803 fHistEventCounter = new TH1D( "fHistEventCounter", ";Evt. Sel. Step;Count",5,0,5);
804 fHistEventCounter->GetXaxis()->SetBinLabel(1, "Processed");
805 fHistEventCounter->GetXaxis()->SetBinLabel(2, "Phys-Sel");
806 fHistEventCounter->GetXaxis()->SetBinLabel(3, "Has Vtx");
807 fHistEventCounter->GetXaxis()->SetBinLabel(4, "Vtx |z|<10cm");
808 fHistEventCounter->GetXaxis()->SetBinLabel(5, "Isn't Pileup");
809 fListHist->Add(fHistEventCounter);
810 }
811
812 //Histograms for Efficiency corrections... a bunch of them
813 //1D Histograms - Fine if efficiency doesn't change vs mult (expected)
814 //---> Always filled for |y|<0.5
161c57ac 815 //V0s: basic Histos
816 if(! fHistPt_GenK0Short ) {
817 fHistPt_GenK0Short = new TH1D( "fHistPt_GenK0Short", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenK0Short); }
818 if(! fHistPt_GenLambda ) {
819 fHistPt_GenLambda = new TH1D( "fHistPt_GenLambda", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenLambda); }
820 if(! fHistPt_GenAntiLambda ) {
821 fHistPt_GenAntiLambda = new TH1D( "fHistPt_GenAntiLambda", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenAntiLambda); }
822 //Cascades: basic Histos
af0d4791 823 if(! fHistPt_GenXiMinus ) {
824 fHistPt_GenXiMinus = new TH1D( "fHistPt_GenXiMinus", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenXiMinus); }
825 if(! fHistPt_GenXiPlus ) {
826 fHistPt_GenXiPlus = new TH1D( "fHistPt_GenXiPlus", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenXiPlus); }
827 if(! fHistPt_GenOmegaMinus ) {
828 fHistPt_GenOmegaMinus = new TH1D( "fHistPt_GenOmegaMinus", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenOmegaMinus); }
829 if(! fHistPt_GenOmegaPlus ) {
830 fHistPt_GenOmegaPlus = new TH1D( "fHistPt_GenOmegaPlus", "Generated;p_{T} (GeV/c)",200,0,20); fListHist->Add(fHistPt_GenOmegaPlus); }
831 //2D Histos for vs Mult calculation
832 if(! fHistPtVsRefMultEta5_GenXiMinus ) {
833 fHistPtVsRefMultEta5_GenXiMinus = new TH2D( "fHistPtVsRefMultEta5_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta5_GenXiMinus); }
834 if(! fHistPtVsRefMultEta5_GenXiPlus ) {
835 fHistPtVsRefMultEta5_GenXiPlus = new TH2D( "fHistPtVsRefMultEta5_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta5_GenXiPlus); }
836 if(! fHistPtVsRefMultEta5_GenOmegaMinus ) {
837 fHistPtVsRefMultEta5_GenOmegaMinus = new TH2D( "fHistPtVsRefMultEta5_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta5_GenOmegaMinus); }
838 if(! fHistPtVsRefMultEta5_GenOmegaPlus ) {
839 fHistPtVsRefMultEta5_GenOmegaPlus = new TH2D( "fHistPtVsRefMultEta5_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta5_GenOmegaPlus); }
840 if(! fHistPtVsRefMultEta8_GenXiMinus ) {
841 fHistPtVsRefMultEta8_GenXiMinus = new TH2D( "fHistPtVsRefMultEta8_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta8_GenXiMinus); }
842 if(! fHistPtVsRefMultEta8_GenXiPlus ) {
843 fHistPtVsRefMultEta8_GenXiPlus = new TH2D( "fHistPtVsRefMultEta8_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta8_GenXiPlus); }
844 if(! fHistPtVsRefMultEta8_GenOmegaMinus ) {
b6b1e756 845 fHistPtVsRefMultEta8_GenOmegaMinus = new TH2D( "fHistPtVsRefMultEta8_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta8_GenOmegaMinus); }
af0d4791 846 if(! fHistPtVsRefMultEta8_GenOmegaPlus ) {
847 fHistPtVsRefMultEta8_GenOmegaPlus = new TH2D( "fHistPtVsRefMultEta8_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100); fListHist->Add(fHistPtVsRefMultEta8_GenOmegaPlus); }
848
849 //Centralities: V0A, V0C, V0M, +Eq
850 if(! fHistPtVsCentV0A_GenXiMinus ) {
851 fHistPtVsCentV0A_GenXiMinus = new TH2D(
852 "fHistPtVsCentV0A_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
853 fListHist->Add(fHistPtVsCentV0A_GenXiMinus); }
854 if(! fHistPtVsCentV0A_GenXiPlus ) {
855 fHistPtVsCentV0A_GenXiPlus = new TH2D(
856 "fHistPtVsCentV0A_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
857 fListHist->Add(fHistPtVsCentV0A_GenXiPlus); }
858 if(! fHistPtVsCentV0A_GenOmegaMinus ) {
859 fHistPtVsCentV0A_GenOmegaMinus = new TH2D(
860 "fHistPtVsCentV0A_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
861 fListHist->Add(fHistPtVsCentV0A_GenOmegaMinus); }
862 if(! fHistPtVsCentV0A_GenOmegaPlus ) {
863 fHistPtVsCentV0A_GenOmegaPlus = new TH2D(
864 "fHistPtVsCentV0A_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
865 fListHist->Add(fHistPtVsCentV0A_GenOmegaPlus); }
866
867 if(! fHistPtVsCentV0C_GenXiMinus ) {
868 fHistPtVsCentV0C_GenXiMinus = new TH2D(
869 "fHistPtVsCentV0C_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
870 fListHist->Add(fHistPtVsCentV0C_GenXiMinus); }
871 if(! fHistPtVsCentV0C_GenXiPlus ) {
872 fHistPtVsCentV0C_GenXiPlus = new TH2D(
873 "fHistPtVsCentV0C_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
874 fListHist->Add(fHistPtVsCentV0C_GenXiPlus); }
875 if(! fHistPtVsCentV0C_GenOmegaMinus ) {
876 fHistPtVsCentV0C_GenOmegaMinus = new TH2D(
877 "fHistPtVsCentV0C_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
878 fListHist->Add(fHistPtVsCentV0C_GenOmegaMinus); }
879 if(! fHistPtVsCentV0C_GenOmegaPlus ) {
880 fHistPtVsCentV0C_GenOmegaPlus = new TH2D(
881 "fHistPtVsCentV0C_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
882 fListHist->Add(fHistPtVsCentV0C_GenOmegaPlus); }
883
884 if(! fHistPtVsCentV0M_GenXiMinus ) {
885 fHistPtVsCentV0M_GenXiMinus = new TH2D(
886 "fHistPtVsCentV0M_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
887 fListHist->Add(fHistPtVsCentV0M_GenXiMinus); }
888 if(! fHistPtVsCentV0M_GenXiPlus ) {
889 fHistPtVsCentV0M_GenXiPlus = new TH2D(
890 "fHistPtVsCentV0M_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
891 fListHist->Add(fHistPtVsCentV0M_GenXiPlus); }
892 if(! fHistPtVsCentV0M_GenOmegaMinus ) {
893 fHistPtVsCentV0M_GenOmegaMinus = new TH2D(
894 "fHistPtVsCentV0M_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
895 fListHist->Add(fHistPtVsCentV0M_GenOmegaMinus); }
896 if(! fHistPtVsCentV0M_GenOmegaPlus ) {
897 fHistPtVsCentV0M_GenOmegaPlus = new TH2D(
898 "fHistPtVsCentV0M_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
899 fListHist->Add(fHistPtVsCentV0M_GenOmegaPlus); }
900
901 //Equalized
902 if(! fHistPtVsCentV0AEq_GenXiMinus ) {
903 fHistPtVsCentV0AEq_GenXiMinus = new TH2D(
904 "fHistPtVsCentV0AEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
905 fListHist->Add(fHistPtVsCentV0AEq_GenXiMinus); }
906 if(! fHistPtVsCentV0AEq_GenXiPlus ) {
907 fHistPtVsCentV0AEq_GenXiPlus = new TH2D(
908 "fHistPtVsCentV0AEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
909 fListHist->Add(fHistPtVsCentV0AEq_GenXiPlus); }
910 if(! fHistPtVsCentV0AEq_GenOmegaMinus ) {
911 fHistPtVsCentV0AEq_GenOmegaMinus = new TH2D(
912 "fHistPtVsCentV0AEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
913 fListHist->Add(fHistPtVsCentV0AEq_GenOmegaMinus); }
914 if(! fHistPtVsCentV0AEq_GenOmegaPlus ) {
915 fHistPtVsCentV0AEq_GenOmegaPlus = new TH2D(
916 "fHistPtVsCentV0AEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
917 fListHist->Add(fHistPtVsCentV0AEq_GenOmegaPlus); }
918
919 if(! fHistPtVsCentV0CEq_GenXiMinus ) {
920 fHistPtVsCentV0CEq_GenXiMinus = new TH2D(
921 "fHistPtVsCentV0CEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
922 fListHist->Add(fHistPtVsCentV0CEq_GenXiMinus); }
923 if(! fHistPtVsCentV0CEq_GenXiPlus ) {
924 fHistPtVsCentV0CEq_GenXiPlus = new TH2D(
925 "fHistPtVsCentV0CEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
926 fListHist->Add(fHistPtVsCentV0CEq_GenXiPlus); }
927 if(! fHistPtVsCentV0CEq_GenOmegaMinus ) {
928 fHistPtVsCentV0CEq_GenOmegaMinus = new TH2D(
929 "fHistPtVsCentV0CEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
930 fListHist->Add(fHistPtVsCentV0CEq_GenOmegaMinus); }
931 if(! fHistPtVsCentV0CEq_GenOmegaPlus ) {
932 fHistPtVsCentV0CEq_GenOmegaPlus = new TH2D(
933 "fHistPtVsCentV0CEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
934 fListHist->Add(fHistPtVsCentV0CEq_GenOmegaPlus); }
935
936 if(! fHistPtVsCentV0MEq_GenXiMinus ) {
937 fHistPtVsCentV0MEq_GenXiMinus = new TH2D(
938 "fHistPtVsCentV0MEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
939 fListHist->Add(fHistPtVsCentV0MEq_GenXiMinus); }
940 if(! fHistPtVsCentV0MEq_GenXiPlus ) {
941 fHistPtVsCentV0MEq_GenXiPlus = new TH2D(
942 "fHistPtVsCentV0MEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
943 fListHist->Add(fHistPtVsCentV0MEq_GenXiPlus); }
944 if(! fHistPtVsCentV0MEq_GenOmegaMinus ) {
945 fHistPtVsCentV0MEq_GenOmegaMinus = new TH2D(
946 "fHistPtVsCentV0MEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
947 fListHist->Add(fHistPtVsCentV0MEq_GenOmegaMinus); }
948 if(! fHistPtVsCentV0MEq_GenOmegaPlus ) {
949 fHistPtVsCentV0MEq_GenOmegaPlus = new TH2D(
950 "fHistPtVsCentV0MEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,100,0,100);
951 fListHist->Add(fHistPtVsCentV0MEq_GenOmegaPlus); }
952
953 //AMPLITUDES: V0A, V0C, V0M, +Eq
954 Double_t lMaxAmplitude = 2500;
955 Long_t lAmplitudeBins = 10000;
956 if(! fHistPtVsAmpV0A_GenXiMinus ) {
957 fHistPtVsAmpV0A_GenXiMinus = new TH2D(
958 "fHistPtVsAmpV0A_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
959 fListHist->Add(fHistPtVsAmpV0A_GenXiMinus); }
960 if(! fHistPtVsAmpV0A_GenXiPlus ) {
961 fHistPtVsAmpV0A_GenXiPlus = new TH2D(
962 "fHistPtVsAmpV0A_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
963 fListHist->Add(fHistPtVsAmpV0A_GenXiPlus); }
964 if(! fHistPtVsAmpV0A_GenOmegaMinus ) {
965 fHistPtVsAmpV0A_GenOmegaMinus = new TH2D(
966 "fHistPtVsAmpV0A_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
967 fListHist->Add(fHistPtVsAmpV0A_GenOmegaMinus); }
968 if(! fHistPtVsAmpV0A_GenOmegaPlus ) {
969 fHistPtVsAmpV0A_GenOmegaPlus = new TH2D(
970 "fHistPtVsAmpV0A_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
971 fListHist->Add(fHistPtVsAmpV0A_GenOmegaPlus); }
972
973 if(! fHistPtVsAmpV0C_GenXiMinus ) {
974 fHistPtVsAmpV0C_GenXiMinus = new TH2D(
975 "fHistPtVsAmpV0C_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
976 fListHist->Add(fHistPtVsAmpV0C_GenXiMinus); }
977 if(! fHistPtVsAmpV0C_GenXiPlus ) {
978 fHistPtVsAmpV0C_GenXiPlus = new TH2D(
979 "fHistPtVsAmpV0C_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
980 fListHist->Add(fHistPtVsAmpV0C_GenXiPlus); }
981 if(! fHistPtVsAmpV0C_GenOmegaMinus ) {
982 fHistPtVsAmpV0C_GenOmegaMinus = new TH2D(
983 "fHistPtVsAmpV0C_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
984 fListHist->Add(fHistPtVsAmpV0C_GenOmegaMinus); }
985 if(! fHistPtVsAmpV0C_GenOmegaPlus ) {
986 fHistPtVsAmpV0C_GenOmegaPlus = new TH2D(
987 "fHistPtVsAmpV0C_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
988 fListHist->Add(fHistPtVsAmpV0C_GenOmegaPlus); }
989
990 if(! fHistPtVsAmpV0M_GenXiMinus ) {
991 fHistPtVsAmpV0M_GenXiMinus = new TH2D(
992 "fHistPtVsAmpV0M_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
993 fListHist->Add(fHistPtVsAmpV0M_GenXiMinus); }
994 if(! fHistPtVsAmpV0M_GenXiPlus ) {
995 fHistPtVsAmpV0M_GenXiPlus = new TH2D(
996 "fHistPtVsAmpV0M_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
997 fListHist->Add(fHistPtVsAmpV0M_GenXiPlus); }
998 if(! fHistPtVsAmpV0M_GenOmegaMinus ) {
999 fHistPtVsAmpV0M_GenOmegaMinus = new TH2D(
1000 "fHistPtVsAmpV0M_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1001 fListHist->Add(fHistPtVsAmpV0M_GenOmegaMinus); }
1002 if(! fHistPtVsAmpV0M_GenOmegaPlus ) {
1003 fHistPtVsAmpV0M_GenOmegaPlus = new TH2D(
1004 "fHistPtVsAmpV0M_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1005 fListHist->Add(fHistPtVsAmpV0M_GenOmegaPlus); }
1006
1007 //Equalized
1008 if(! fHistPtVsAmpV0AEq_GenXiMinus ) {
1009 fHistPtVsAmpV0AEq_GenXiMinus = new TH2D(
1010 "fHistPtVsAmpV0AEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1011 fListHist->Add(fHistPtVsAmpV0AEq_GenXiMinus); }
1012 if(! fHistPtVsAmpV0AEq_GenXiPlus ) {
1013 fHistPtVsAmpV0AEq_GenXiPlus = new TH2D(
1014 "fHistPtVsAmpV0AEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1015 fListHist->Add(fHistPtVsAmpV0AEq_GenXiPlus); }
1016 if(! fHistPtVsAmpV0AEq_GenOmegaMinus ) {
1017 fHistPtVsAmpV0AEq_GenOmegaMinus = new TH2D(
1018 "fHistPtVsAmpV0AEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1019 fListHist->Add(fHistPtVsAmpV0AEq_GenOmegaMinus); }
1020 if(! fHistPtVsAmpV0AEq_GenOmegaPlus ) {
1021 fHistPtVsAmpV0AEq_GenOmegaPlus = new TH2D(
1022 "fHistPtVsAmpV0AEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1023 fListHist->Add(fHistPtVsAmpV0AEq_GenOmegaPlus); }
1024
1025 if(! fHistPtVsAmpV0CEq_GenXiMinus ) {
1026 fHistPtVsAmpV0CEq_GenXiMinus = new TH2D(
1027 "fHistPtVsAmpV0CEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1028 fListHist->Add(fHistPtVsAmpV0CEq_GenXiMinus); }
1029 if(! fHistPtVsAmpV0CEq_GenXiPlus ) {
1030 fHistPtVsAmpV0CEq_GenXiPlus = new TH2D(
1031 "fHistPtVsAmpV0CEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1032 fListHist->Add(fHistPtVsAmpV0CEq_GenXiPlus); }
1033 if(! fHistPtVsAmpV0CEq_GenOmegaMinus ) {
1034 fHistPtVsAmpV0CEq_GenOmegaMinus = new TH2D(
1035 "fHistPtVsAmpV0CEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1036 fListHist->Add(fHistPtVsAmpV0CEq_GenOmegaMinus); }
1037 if(! fHistPtVsAmpV0CEq_GenOmegaPlus ) {
1038 fHistPtVsAmpV0CEq_GenOmegaPlus = new TH2D(
1039 "fHistPtVsAmpV0CEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1040 fListHist->Add(fHistPtVsAmpV0CEq_GenOmegaPlus); }
1041
1042 if(! fHistPtVsAmpV0MEq_GenXiMinus ) {
1043 fHistPtVsAmpV0MEq_GenXiMinus = new TH2D(
1044 "fHistPtVsAmpV0MEq_GenXiMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1045 fListHist->Add(fHistPtVsAmpV0MEq_GenXiMinus); }
1046 if(! fHistPtVsAmpV0MEq_GenXiPlus ) {
1047 fHistPtVsAmpV0MEq_GenXiPlus = new TH2D(
1048 "fHistPtVsAmpV0MEq_GenXiPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1049 fListHist->Add(fHistPtVsAmpV0MEq_GenXiPlus); }
1050 if(! fHistPtVsAmpV0MEq_GenOmegaMinus ) {
1051 fHistPtVsAmpV0MEq_GenOmegaMinus = new TH2D(
1052 "fHistPtVsAmpV0MEq_GenOmegaMinus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1053 fListHist->Add(fHistPtVsAmpV0MEq_GenOmegaMinus); }
1054 if(! fHistPtVsAmpV0MEq_GenOmegaPlus ) {
1055 fHistPtVsAmpV0MEq_GenOmegaPlus = new TH2D(
1056 "fHistPtVsAmpV0MEq_GenOmegaPlus", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1057 fListHist->Add(fHistPtVsAmpV0MEq_GenOmegaPlus); }
1058
1dac6d41 1059 if(! fHistVZEROResponseStudy ) {
1060 fHistVZEROResponseStudy = new TH2D(
1061 "fHistVZEROResponseStudy", "Generated;p_{T} (GeV/c); Mult",200,0,20,lAmplitudeBins,0,lMaxAmplitude);
1062 fListHist->Add(fHistVZEROResponseStudy); }
1063
1064 if(! fHistVZEROResponseStudyTotal ) {
1065 fHistVZEROResponseStudyTotal = new TH2D(
a190e35b 1066 "fHistVZEROResponseStudyTotal", "Generated;p_{T} (GeV/c); Mult",5000,0,500,lAmplitudeBins,0,lMaxAmplitude);
1dac6d41 1067 fListHist->Add(fHistVZEROResponseStudyTotal); }
1068
af0d4791 1069 //List of Histograms: Normal
1070 PostData(1, fListHist);
1071
1072 //TTree Object: Saved to base directory. Should cache to disk while saving.
1073 //(Important to avoid excessive memory usage, particularly when merging)
1074 PostData(2, fTreeEvent);
1075 PostData(3, fTreeV0);
1076 PostData(4, fTreeCascade);
1077
1078}// end UserCreateOutputObjects
1079
1080
1081//________________________________________________________________________
1082void AliAnalysisTaskStrangenessVsMultiplicityMC::UserExec(Option_t *)
1083{
1084 // Main loop
1085 // Called for each event
1086
1087 AliESDEvent *lESDevent = 0x0;
1088 AliMCEvent *lMCevent = 0x0;
1089 AliStack *lMCstack = 0x0;
2ca41fef 1090
5f31c5d4 1091 //Zero all booleans, etc
1092 fEvSel_HasAtLeastSPDVertex = kFALSE;
1093 fEvSel_VtxZCut = kFALSE;
1094 fEvSel_IsNotPileup = kFALSE;
1095 fEvSel_IsNotPileupInMultBins = kFALSE;
1096 fEvSel_HasVtxContributor = kFALSE;
1097 fEvSel_Triggered = kFALSE;
2ca41fef 1098 fEvSel_VtxZ = -100;
1099 fEvSel_MCType = -100;
af0d4791 1100 // Connect to the InputEvent
1101 // After these lines, we should have an ESD/AOD event + the number of V0s in it.
1102
1103 // Appropriate for ESD analysis!
1104
1105 lESDevent = dynamic_cast<AliESDEvent*>( InputEvent() );
1106 if (!lESDevent) {
1107 AliWarning("ERROR: lESDevent not available \n");
1108 return;
1109 }
1110
1111 //Get VZERO Information for multiplicity later
1112 AliVVZERO* esdV0 = lESDevent->GetVZEROData();
1113 if (!esdV0) {
1114 AliError("AliVVZERO not available");
1115 return;
1116 }
1117
1118 lMCevent = MCEvent();
1119 if (!lMCevent) {
1120 Printf("ERROR: Could not retrieve MC event \n");
1121 cout << "Name of the file with pb :" << fInputHandler->GetTree()->GetCurrentFile()->GetName() << endl;
1122 return;
1123 }
1124
1125 lMCstack = lMCevent->Stack();
1126 if (!lMCstack) {
1127 Printf("ERROR: Could not retrieve MC stack \n");
1128 cout << "Name of the file with pb :" << fInputHandler->GetTree()->GetCurrentFile()->GetName() << endl;
1129 return;
1130 }
1131
1132 fRunNumber = lESDevent->GetRunNumber();
1133
1134 Double_t lMagneticField = -10;
1135 lMagneticField = lESDevent->GetMagneticField( );
1136
1137//------------------------------------------------
2ca41fef 1138// MC type (ND, SD, DD)
af0d4791 1139//------------------------------------------------
1140
2ca41fef 1141 AliGenEventHeader * header = lMCevent->GenEventHeader();
1142 Int_t processtype = AliPWG0Helper::GetPythiaEventProcessType(header);
1143 // non diffractive
1144 if (processtype !=92 && processtype !=93 && processtype != 94) fEvSel_MCType = 1;
1145 // single diffractive
1146 if ((processtype == 92 || processtype == 93)) fEvSel_MCType = 2;
1147 // double diffractive
1148 if (processtype == 94) fEvSel_MCType = 3;
1149
af0d4791 1150//------------------------------------------------
1151// Physics Selection
1152//------------------------------------------------
1153
1154 fHistEventCounter->Fill(0.5);
1155
1156 UInt_t maskIsSelected = ((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected();
1157 Bool_t isSelected = 0;
1158 isSelected = (maskIsSelected & AliVEvent::kMB) == AliVEvent::kMB;
dedb9266 1159 fEvSel_Triggered = isSelected;
af0d4791 1160
1161 //Standard Min-Bias Selection
dedb9266 1162 if ( (! isSelected) && (! fkSkipEventSelection ) ) {
af0d4791 1163 PostData(1, fListHist);
1164 PostData(2, fTreeEvent);
1165 PostData(3, fTreeV0);
1166 PostData(4, fTreeCascade);
1167 return;
1168 }
1169
1170 fHistEventCounter->Fill(1.5);
1171
1172 //------------------------------------------------
1173 // Primary Vertex Requirements Section:
1174 // ---> pp: has vertex, |z|<10cm
1175 //------------------------------------------------
1176
1177 //classical Proton-proton like selection
1178 const AliESDVertex *lPrimaryBestESDVtx = lESDevent->GetPrimaryVertex();
1179 const AliESDVertex *lPrimaryTrackingESDVtx = lESDevent->GetPrimaryVertexTracks();
1180 const AliESDVertex *lPrimarySPDVtx = lESDevent->GetPrimaryVertexSPD();
1181
1182 Double_t lBestPrimaryVtxPos[3] = {-100.0, -100.0, -100.0};
1183 lPrimaryBestESDVtx->GetXYZ( lBestPrimaryVtxPos );
2ca41fef 1184 fEvSel_VtxZ = lBestPrimaryVtxPos[2];
af0d4791 1185
1186 //Only accept if Tracking or SPD vertex is fine
dedb9266 1187 if (!lPrimarySPDVtx->GetStatus() && !lPrimaryTrackingESDVtx->GetStatus() && !fkSkipEventSelection ){
af0d4791 1188 AliWarning("Pb / No SPD prim. vertex nor prim. Tracking vertex ... return !");
1189 PostData(1, fListHist);
1190 PostData(2, fTreeEvent);
1191 PostData(3, fTreeV0);
1192 PostData(4, fTreeCascade);
1193 return;
1194 }
dedb9266 1195
1196 if(! (!lPrimarySPDVtx->GetStatus() && !lPrimaryTrackingESDVtx->GetStatus()) ){
1197 //Passed selection!
1198 fEvSel_HasAtLeastSPDVertex = kTRUE;
1199 }
1200
af0d4791 1201 //Has SPD or Tracking Vertex
1202 fHistEventCounter -> Fill(2.5);
1203
1204 //Always do Primary Vertex Selection
dedb9266 1205 if(TMath::Abs(lBestPrimaryVtxPos[2]) > 10.0 && !fkSkipEventSelection ) {
af0d4791 1206 AliWarning("Pb / | Z position of Best Prim Vtx | > 10.0 cm ... return !");
1207 PostData(1, fListHist);
1208 PostData(2, fTreeEvent);
1209 PostData(3, fTreeV0);
1210 PostData(4, fTreeCascade);
1211 return;
1212 }
dedb9266 1213
1214 if(TMath::Abs(lBestPrimaryVtxPos[2]) <= 10.0 ){
1215 //Passed selection!
1216 fEvSel_VtxZCut = kTRUE;
1217 }
af0d4791 1218
1219 //Fill Event selected counter
1220 fHistEventCounter -> Fill(3.5);
1221
1222 //------------------------------------------------
1223 // Check if this isn't pileup
1224 //------------------------------------------------
1225
535c1213 1226 if(lESDevent->IsPileupFromSPD() && !fkSkipEventSelection ){
af0d4791 1227 // minContributors=3, minZdist=0.8, nSigmaZdist=3., nSigmaDiamXY=2., nSigmaDiamZ=5.
1228 //-> see http://alisoft.cern.ch/viewvc/trunk/STEER/AliESDEvent.h?root=AliRoot&r1=41914&r2=42199&pathrev=42199
1229 AliWarning("Pb / Event tagged as pile-up by SPD... return !");
1230 PostData(1, fListHist);
1231 PostData(2, fTreeEvent);
1232 PostData(3, fTreeV0);
1233 PostData(4, fTreeCascade);
1234 return;
1235 }
dedb9266 1236
1237 if( !lESDevent->IsPileupFromSPD() ) fEvSel_IsNotPileup = kTRUE;
1238 if( !lESDevent->IsPileupFromSPDInMultBins() ) fEvSel_IsNotPileupInMultBins = kTRUE;
1239
1240 //First implementation of pileup from multi-vertexer (simple use of analysis utils)
1241 //if ( !fUtils->IsPileUpMV( lESDevent ) ) fEvSel_IsNotPileupMV = kTRUE;
1242 fEvSel_IsNotPileupMV = kFALSE ; //dummy
1243
af0d4791 1244 //Fill Event isn't pileup counter
1245 fHistEventCounter -> Fill(4.5);
1246
1247//------------------------------------------------
1248// Multiplicity Information Acquistion
1249//------------------------------------------------
1250
1251 //Monte Carlo Level information !
1252 //--------- GENERATED NUMBER OF CHARGED PARTICLES
1253 // ---> Variable Definition
1254
1255 Long_t lNchEta5 = 0;
1256 Long_t lNchEta8 = 0;
1257 Long_t lNchVZEROA = 0;
1258 Long_t lNchVZEROC = 0;
1259
1dac6d41 1260 Float_t lPtOfParticleInsideVZEROA = -1;
1261 Float_t lPOfParticleInsideVZEROA = -1;
1262
af0d4791 1263 //----- Loop on Stack ----------------------------------------------------------------
1264 for (Int_t iCurrentLabelStack = 0; iCurrentLabelStack < (lMCstack->GetNtrack()); iCurrentLabelStack++)
1265 {// This is the begining of the loop on tracks
1266 TParticle* particleOne = lMCstack->Particle(iCurrentLabelStack);
1267 if(!particleOne) continue;
1268 if(!particleOne->GetPDG()) continue;
1269 Double_t lThisCharge = particleOne->GetPDG()->Charge()/3.;
1270 if(TMath::Abs(lThisCharge)<0.001) continue;
1271 if(! (lMCstack->IsPhysicalPrimary(iCurrentLabelStack)) ) continue;
1272
1273 //Double_t gpt = particleOne -> Pt();
1274 Double_t geta = particleOne -> Eta();
1275
1276 if( TMath::Abs(geta) < 0.5 ) lNchEta5++;
1277 if( TMath::Abs(geta) < 0.8 ) lNchEta8++;
1278 if( 2.8 < geta && geta < 5.1 ) lNchVZEROA++;
1dac6d41 1279 if( 2.8 < geta && geta < 5.1 ) lPtOfParticleInsideVZEROA = particleOne->Pt();
1280 if( 2.8 < geta && geta < 5.1 ) lPOfParticleInsideVZEROA = particleOne->P();
af0d4791 1281 if(-3.7 < geta && geta <-1.7 ) lNchVZEROC++;
1282 }//End of loop on tracks
1dac6d41 1283
af0d4791 1284 //Attribution
1285 fTrueMultEta5 = lNchEta5;
1286 fTrueMultEta8 = lNchEta8;
1287 fTrueMultVZEROA = lNchVZEROA;
1288 fTrueMultVZEROC = lNchVZEROC;
1289 //----- End Loop on Stack ------------------------------------------------------------
1290
1291 //Standard GetReferenceMultiplicity Estimator (0.5 and 0.8)
1292 fRefMultEta5 = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,0.5);
1293 fRefMultEta8 = fESDtrackCuts->GetReferenceMultiplicity(lESDevent, AliESDtrackCuts::kTrackletsITSTPC,0.8);
1294
1295 // VZERO PART
1296 Float_t multV0A = 0; // multiplicity from V0 reco side A
1297 Float_t multV0C = 0; // multiplicity from V0 reco side C
1298 Float_t multV0AEq = 0; // multiplicity from V0 reco side A
1299 Float_t multV0CEq = 0; // multiplicity from V0 reco side C
1300 Float_t multV0ACorr = 0; // multiplicity from V0 reco side A
1301 Float_t multV0CCorr = 0; // multiplicity from V0 reco side C
1302
1303 //Non-Equalized Signal: copy of multV0ACorr and multV0CCorr from AliCentralitySelectionTask
1304 //Getters for uncorrected multiplicity
1305 multV0A=esdV0->GetMTotV0A();
1306 multV0C=esdV0->GetMTotV0C();
1307
1308 //Get Z vertex position of SPD vertex (why not Tracking if available?)
1309 Float_t zvtx = lPrimarySPDVtx->GetZ();
1310
1311 //Acquire Corrected multV0A
1312 multV0ACorr = AliESDUtils::GetCorrV0A(multV0A,zvtx);
1313 multV0CCorr = AliESDUtils::GetCorrV0C(multV0C,zvtx);
1314
1315 //Copy to Event Tree for extra information
1316 fAmplitude_V0A = multV0ACorr;
1317 fAmplitude_V0C = multV0CCorr;
1318
1dac6d41 1319 if( fTrueMultVZEROA == 1 ) fHistVZEROResponseStudy->Fill( lPtOfParticleInsideVZEROA , fAmplitude_V0A );
1320 if( fTrueMultVZEROA == 1 ) fHistVZEROResponseStudyTotal->Fill( lPOfParticleInsideVZEROA , fAmplitude_V0A );
1321
af0d4791 1322 // Equalized signals // From AliCentralitySelectionTask
1323 for(Int_t iCh = 4; iCh < 7; ++iCh) {
1324 Double_t mult = lESDevent->GetVZEROEqMultiplicity(iCh);
1325 multV0AEq += mult;
1326 }
1327 for(Int_t iCh = 0; iCh < 3; ++iCh) {
1328 Double_t mult = lESDevent->GetVZEROEqMultiplicity(iCh);
1329 multV0CEq += mult;
1330 }
1331 fAmplitude_V0AEq = multV0AEq;
1332 fAmplitude_V0CEq = multV0CEq;
1333
1334 fCentrality_V0A = -100;
1335 fCentrality_V0C = -100;
1336 fCentrality_V0M = -100;
1337 fCentrality_V0AEq = -100;
1338 fCentrality_V0CEq = -100;
1339 fCentrality_V0MEq = -100;
1340
1341 //AliCentrality... Check if working?
1342 AliCentrality* centrality;
1343 centrality = lESDevent->GetCentrality();
1344 if ( !(centrality->GetQuality()>1) ){
1345 fCentrality_V0A = centrality->GetCentralityPercentile( "V0A" );
1346 fCentrality_V0C = centrality->GetCentralityPercentile( "V0C" );
1347 fCentrality_V0M = centrality->GetCentralityPercentile( "V0M" );
1348 fCentrality_V0AEq = centrality->GetCentralityPercentile( "V0AEq" );
1349 fCentrality_V0CEq = centrality->GetCentralityPercentile( "V0CEq" );
1350 fCentrality_V0MEq = centrality->GetCentralityPercentile( "V0MEq" );
1351 }
1352
535c1213 1353 fCustomCentrality_V0M = fPPVsMultUtils -> GetMultiplicityPercentile(lESDevent, "V0M");
1354 fCustomCentrality_V0MEq = fPPVsMultUtils -> GetMultiplicityPercentile(lESDevent, "V0MEq");
1355
9a28e152 1356 fTreeVariableCustomCentV0M = fCustomCentrality_V0M;
1357 fTreeCascVarCustomCentV0M = fCustomCentrality_V0M;
1358
af0d4791 1359 //Event-level fill
535c1213 1360 fTreeEvent->Fill();
1361
1362 //STOP HERE if skipping event selections (no point in doing the rest...)
1363 if( fkSkipEventSelection ){
1364 PostData(1, fListHist);
1365 PostData(2, fTreeEvent);
1366 PostData(3, fTreeV0);
1367 PostData(4, fTreeCascade);
1368 return;
1369 }
1370
af0d4791 1371
1372//------------------------------------------------
1373
1374//------------------------------------------------
1375// Fill Efficiency Denominators, please
1376//------------------------------------------------
1377
1378 Int_t lThisPDG = 0;
1379 Double_t lThisRap = 0;
1380 Double_t lThisPt = 0;
1381
1382//----- Loop on Generated CASCADES ---------------
1383 for (Int_t ilab = 0; ilab < (lMCstack->GetNtrack()); ilab++)
1384 {// This is the begining of the loop on tracks
1385
1386 TParticle* lPart = 0x0;
1387 lPart = lMCstack->Particle( ilab );
1388 if(!lPart){
1389 Printf("Generated loop %d - MC TParticle pointer to current stack particle = 0x0 ! Skip ...\n", ilab );
1390 continue;
1391 }
1392
1393 lThisPDG = lPart->GetPdgCode();
1394
161c57ac 1395 if ( (TMath::Abs(lThisPDG) == 3312) || (TMath::Abs(lThisPDG) == 3334) || (TMath::Abs(lThisPDG) == 3122) || lThisPDG == 310 )
af0d4791 1396 {
1397 lThisRap = MyRapidity(lPart->Energy(),lPart->Pz());
1398 lThisPt = lPart->Pt();
1399
1400 //Use Physical Primaries only for filling These Histos
1401 if ( lMCstack->IsPhysicalPrimary(ilab)!=kTRUE ) continue;
1402
161c57ac 1403 if( lThisPDG == 310 && TMath::Abs(lThisRap) < 0.5 ) fHistPt_GenK0Short -> Fill ( lThisPt );
1404 if( lThisPDG == 3122 && TMath::Abs(lThisRap) < 0.5 ) fHistPt_GenLambda -> Fill ( lThisPt );
1405 if( lThisPDG == -3122 && TMath::Abs(lThisRap) < 0.5 ) fHistPt_GenAntiLambda -> Fill ( lThisPt );
1406
af0d4791 1407 if( lThisPDG == 3312 && TMath::Abs(lThisRap) < 0.5 ){
1408 fHistPt_GenXiMinus -> Fill (lThisPt);
1409 fHistPtVsRefMultEta5_GenXiMinus -> Fill (lThisPt, fRefMultEta5);
1410 fHistPtVsRefMultEta8_GenXiMinus -> Fill (lThisPt, fRefMultEta8);
1411 //Centralities
1412 fHistPtVsCentV0A_GenXiMinus -> Fill (lThisPt, fCentrality_V0A);
1413 fHistPtVsCentV0C_GenXiMinus -> Fill (lThisPt, fCentrality_V0C);
1414 fHistPtVsCentV0M_GenXiMinus -> Fill (lThisPt, fCentrality_V0M);
1415 fHistPtVsCentV0AEq_GenXiMinus -> Fill (lThisPt, fCentrality_V0AEq);
1416 fHistPtVsCentV0CEq_GenXiMinus -> Fill (lThisPt, fCentrality_V0CEq);
1417 fHistPtVsCentV0MEq_GenXiMinus -> Fill (lThisPt, fCentrality_V0MEq);
1418 //Amplitudes
1419 fHistPtVsAmpV0A_GenXiMinus -> Fill (lThisPt, fAmplitude_V0A);
1420 fHistPtVsAmpV0C_GenXiMinus -> Fill (lThisPt, fAmplitude_V0C);
1421 fHistPtVsAmpV0M_GenXiMinus -> Fill (lThisPt, fAmplitude_V0A + fAmplitude_V0C);
1422 fHistPtVsAmpV0AEq_GenXiMinus -> Fill (lThisPt, fAmplitude_V0AEq);
1423 fHistPtVsAmpV0CEq_GenXiMinus -> Fill (lThisPt, fAmplitude_V0CEq);
1424 fHistPtVsAmpV0MEq_GenXiMinus -> Fill (lThisPt, fAmplitude_V0AEq + fAmplitude_V0CEq);
1425 }
1426 if( lThisPDG == -3312 && TMath::Abs(lThisRap) < 0.5 ){
1427 fHistPt_GenXiPlus -> Fill (lThisPt);
1428 fHistPtVsRefMultEta5_GenXiPlus -> Fill (lThisPt, fRefMultEta5);
1429 fHistPtVsRefMultEta8_GenXiPlus -> Fill (lThisPt, fRefMultEta8);
1430 //Centralities
1431 fHistPtVsCentV0A_GenXiPlus -> Fill (lThisPt, fCentrality_V0A);
1432 fHistPtVsCentV0C_GenXiPlus -> Fill (lThisPt, fCentrality_V0C);
1433 fHistPtVsCentV0M_GenXiPlus -> Fill (lThisPt, fCentrality_V0M);
1434 fHistPtVsCentV0AEq_GenXiPlus -> Fill (lThisPt, fCentrality_V0AEq);
1435 fHistPtVsCentV0CEq_GenXiPlus -> Fill (lThisPt, fCentrality_V0CEq);
1436 fHistPtVsCentV0MEq_GenXiPlus -> Fill (lThisPt, fCentrality_V0MEq);
1437 //Amplitudes
1438 fHistPtVsAmpV0A_GenXiPlus -> Fill (lThisPt, fAmplitude_V0A);
1439 fHistPtVsAmpV0C_GenXiPlus -> Fill (lThisPt, fAmplitude_V0C);
1440 fHistPtVsAmpV0M_GenXiPlus -> Fill (lThisPt, fAmplitude_V0A + fAmplitude_V0C);
1441 fHistPtVsAmpV0AEq_GenXiPlus -> Fill (lThisPt, fAmplitude_V0AEq);
1442 fHistPtVsAmpV0CEq_GenXiPlus -> Fill (lThisPt, fAmplitude_V0CEq);
1443 fHistPtVsAmpV0MEq_GenXiPlus -> Fill (lThisPt, fAmplitude_V0AEq + fAmplitude_V0CEq);
1444 }
1445 if( lThisPDG == 3334 && TMath::Abs(lThisRap) < 0.5 ){
1446 fHistPt_GenOmegaMinus -> Fill (lThisPt);
1447 fHistPtVsRefMultEta5_GenOmegaMinus -> Fill (lThisPt, fRefMultEta5);
1448 fHistPtVsRefMultEta8_GenOmegaMinus -> Fill (lThisPt, fRefMultEta8);
1449 //Centralities
1450 fHistPtVsCentV0A_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0A);
1451 fHistPtVsCentV0C_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0C);
1452 fHistPtVsCentV0M_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0M);
1453 fHistPtVsCentV0AEq_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0AEq);
1454 fHistPtVsCentV0CEq_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0CEq);
1455 fHistPtVsCentV0MEq_GenOmegaMinus -> Fill (lThisPt, fCentrality_V0MEq);
1456 //Amplitudes
1457 fHistPtVsAmpV0A_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0A);
1458 fHistPtVsAmpV0C_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0C);
1459 fHistPtVsAmpV0M_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0A + fAmplitude_V0C);
1460 fHistPtVsAmpV0AEq_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0AEq);
1461 fHistPtVsAmpV0CEq_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0CEq);
1462 fHistPtVsAmpV0MEq_GenOmegaMinus -> Fill (lThisPt, fAmplitude_V0AEq + fAmplitude_V0CEq);
1463 }
1464 if( lThisPDG == -3334 && TMath::Abs(lThisRap) < 0.5 ){
1465 fHistPt_GenOmegaPlus -> Fill (lThisPt);
1466 fHistPtVsRefMultEta5_GenOmegaPlus -> Fill (lThisPt, fRefMultEta5);
1467 fHistPtVsRefMultEta8_GenOmegaPlus -> Fill (lThisPt, fRefMultEta8);
1468 //Centralities
1469 fHistPtVsCentV0A_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0A);
1470 fHistPtVsCentV0C_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0C);
1471 fHistPtVsCentV0M_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0M);
1472 fHistPtVsCentV0AEq_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0AEq);
1473 fHistPtVsCentV0CEq_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0CEq);
1474 fHistPtVsCentV0MEq_GenOmegaPlus -> Fill (lThisPt, fCentrality_V0MEq);
1475 //Amplitudes
1476 fHistPtVsAmpV0A_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0A);
1477 fHistPtVsAmpV0C_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0C);
1478 fHistPtVsAmpV0M_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0A + fAmplitude_V0C);
1479 fHistPtVsAmpV0AEq_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0AEq);
1480 fHistPtVsAmpV0CEq_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0CEq);
1481 fHistPtVsAmpV0MEq_GenOmegaPlus -> Fill (lThisPt, fAmplitude_V0AEq + fAmplitude_V0CEq);
1482 }
1483 }
1484 }//End of loop on tracks
1485//----- End Loop on Cascades ------------------------------------------------------------
1486
161c57ac 1487//------------------------------------------------
1488// Fill V0 Tree as needed
1489//------------------------------------------------
1490
1491//Variable definition
1492 Int_t lOnFlyStatus = 0;// nv0sOn = 0, nv0sOff = 0;
1493 Double_t lChi2V0 = 0;
1494 Double_t lDcaV0Daughters = 0, lDcaV0ToPrimVertex = 0;
1495 Double_t lDcaPosToPrimVertex = 0, lDcaNegToPrimVertex = 0;
1496 Double_t lV0CosineOfPointingAngle = 0;
1497 Double_t lV0Radius = 0, lPt = 0;
1498 Double_t lRapK0Short = 0, lRapLambda = 0;
1499 Double_t lInvMassK0s = 0, lInvMassLambda = 0, lInvMassAntiLambda = 0;
1500 Double_t lAlphaV0 = 0, lPtArmV0 = 0;
1501
1502 Double_t fMinV0Pt = 0;
1503 Double_t fMaxV0Pt = 100;
1504
1505 Int_t nv0s = 0;
1506 nv0s = lESDevent->GetNumberOfV0s();
1507
1508 for (Int_t iV0 = 0; iV0 < nv0s; iV0++) //extra-crazy test
1509 {// This is the begining of the V0 loop
1510 AliESDv0 *v0 = ((AliESDEvent*)lESDevent)->GetV0(iV0);
1511 if (!v0) continue;
1512
1513 Double_t tDecayVertexV0[3]; v0->GetXYZ(tDecayVertexV0[0],tDecayVertexV0[1],tDecayVertexV0[2]);
1514
1515 Double_t tV0mom[3];
1516 v0->GetPxPyPz( tV0mom[0],tV0mom[1],tV0mom[2] );
1517 Double_t lV0TotalMomentum = TMath::Sqrt(
1518 tV0mom[0]*tV0mom[0]+tV0mom[1]*tV0mom[1]+tV0mom[2]*tV0mom[2] );
1519
1520 lV0Radius = TMath::Sqrt(tDecayVertexV0[0]*tDecayVertexV0[0]+tDecayVertexV0[1]*tDecayVertexV0[1]);
1521
1522 lPt = v0->Pt();
1523 lRapK0Short = v0->RapK0Short();
1524 lRapLambda = v0->RapLambda();
1525 if ((lPt<fMinV0Pt)||(fMaxV0Pt<lPt)) continue;
1526
1527 UInt_t lKeyPos = (UInt_t)TMath::Abs(v0->GetPindex());
1528 UInt_t lKeyNeg = (UInt_t)TMath::Abs(v0->GetNindex());
1529
1530 Double_t lMomPos[3]; v0->GetPPxPyPz(lMomPos[0],lMomPos[1],lMomPos[2]);
1531 Double_t lMomNeg[3]; v0->GetNPxPyPz(lMomNeg[0],lMomNeg[1],lMomNeg[2]);
1532
1533 AliESDtrack *pTrack=((AliESDEvent*)lESDevent)->GetTrack(lKeyPos);
1534 AliESDtrack *nTrack=((AliESDEvent*)lESDevent)->GetTrack(lKeyNeg);
1535 if (!pTrack || !nTrack) {
1536 Printf("ERROR: Could not retreive one of the daughter track");
1537 continue;
1538 }
1539
1540 //Daughter Eta for Eta selection, afterwards
1541 fTreeVariableNegEta = nTrack->Eta();
1542 fTreeVariablePosEta = pTrack->Eta();
1543
1544 // Filter like-sign V0 (next: add counter and distribution)
1545 if ( pTrack->GetSign() == nTrack->GetSign()){
1546 continue;
1547 }
1548
1549 //________________________________________________________________________
1550 // Track quality cuts
1551 Float_t lPosTrackCrossedRows = pTrack->GetTPCClusterInfo(2,1);
1552 Float_t lNegTrackCrossedRows = nTrack->GetTPCClusterInfo(2,1);
1553 fTreeVariableLeastNbrCrossedRows = (Int_t) lPosTrackCrossedRows;
1554 if( lNegTrackCrossedRows < fTreeVariableLeastNbrCrossedRows )
1555 fTreeVariableLeastNbrCrossedRows = (Int_t) lNegTrackCrossedRows;
1556
1557 // TPC refit condition (done during reconstruction for Offline but not for On-the-fly)
1558 if( !(pTrack->GetStatus() & AliESDtrack::kTPCrefit)) continue;
1559 if( !(nTrack->GetStatus() & AliESDtrack::kTPCrefit)) continue;
1560
1561
1562 if ( ( ( pTrack->GetTPCClusterInfo(2,1) ) < 70 ) || ( ( nTrack->GetTPCClusterInfo(2,1) ) < 70 ) ) continue;
1563
1564 //GetKinkIndex condition
1565 if( pTrack->GetKinkIndex(0)>0 || nTrack->GetKinkIndex(0)>0 ) continue;
1566
1567 //Findable clusters > 0 condition
1568 if( pTrack->GetTPCNclsF()<=0 || nTrack->GetTPCNclsF()<=0 ) continue;
1569
1570 //Compute ratio Crossed Rows / Findable clusters
1571 //Note: above test avoids division by zero!
1572 Float_t lPosTrackCrossedRowsOverFindable = lPosTrackCrossedRows / ((double)(pTrack->GetTPCNclsF()));
1573 Float_t lNegTrackCrossedRowsOverFindable = lNegTrackCrossedRows / ((double)(nTrack->GetTPCNclsF()));
1574
1575 fTreeVariableLeastRatioCrossedRowsOverFindable = lPosTrackCrossedRowsOverFindable;
1576 if( lNegTrackCrossedRowsOverFindable < fTreeVariableLeastRatioCrossedRowsOverFindable )
1577 fTreeVariableLeastRatioCrossedRowsOverFindable = lNegTrackCrossedRowsOverFindable;
1578
1579 //Lowest Cut Level for Ratio Crossed Rows / Findable = 0.8, set here
1580 if ( fTreeVariableLeastRatioCrossedRowsOverFindable < 0.8 ) continue;
1581
1582 //End track Quality Cuts
1583 //________________________________________________________________________
1584
1585 lDcaPosToPrimVertex = TMath::Abs(pTrack->GetD(lBestPrimaryVtxPos[0],
1586 lBestPrimaryVtxPos[1],
1587 lMagneticField) );
1588
1589 lDcaNegToPrimVertex = TMath::Abs(nTrack->GetD(lBestPrimaryVtxPos[0],
1590 lBestPrimaryVtxPos[1],
1591 lMagneticField) );
1592
1593 lOnFlyStatus = v0->GetOnFlyStatus();
1594 lChi2V0 = v0->GetChi2V0();
1595 lDcaV0Daughters = v0->GetDcaV0Daughters();
1596 lDcaV0ToPrimVertex = v0->GetD(lBestPrimaryVtxPos[0],lBestPrimaryVtxPos[1],lBestPrimaryVtxPos[2]);
1597 lV0CosineOfPointingAngle = v0->GetV0CosineOfPointingAngle(lBestPrimaryVtxPos[0],lBestPrimaryVtxPos[1],lBestPrimaryVtxPos[2]);
1598 fTreeVariableV0CosineOfPointingAngle=lV0CosineOfPointingAngle;
1599
1600 // Getting invariant mass infos directly from ESD
1601 v0->ChangeMassHypothesis(310);
1602 lInvMassK0s = v0->GetEffMass();
1603 v0->ChangeMassHypothesis(3122);
1604 lInvMassLambda = v0->GetEffMass();
1605 v0->ChangeMassHypothesis(-3122);
1606 lInvMassAntiLambda = v0->GetEffMass();
1607 lAlphaV0 = v0->AlphaV0();
1608 lPtArmV0 = v0->PtArmV0();
1609
1610
1611//===============================================
1612// Monte Carlo Association starts here
1613//===============================================
1614
1615 //---> Set Everything to "I don't know" before starting
1616
1617 fTreeVariablePIDPositive = 0;
1618 fTreeVariablePIDNegative = 0;
1619
1620 fTreeVariablePtMother = -1;
1621 fTreeVariablePtMC = -1;
1622 fTreeVariableRapMC = -100;
1623
1624 fTreeVariablePID = -1;
1625 fTreeVariablePIDMother = -1;
1626
1627 fTreeVariablePrimaryStatus = 0;
1628 fTreeVariablePrimaryStatusMother = 0;
1629
1630 Int_t lblPosV0Dghter = (Int_t) TMath::Abs( pTrack->GetLabel() );
1631 Int_t lblNegV0Dghter = (Int_t) TMath::Abs( nTrack->GetLabel() );
1632
1633 TParticle* mcPosV0Dghter = lMCstack->Particle( lblPosV0Dghter );
1634 TParticle* mcNegV0Dghter = lMCstack->Particle( lblNegV0Dghter );
1635
1636 Int_t lPIDPositive = mcPosV0Dghter -> GetPdgCode();
1637 Int_t lPIDNegative = mcNegV0Dghter -> GetPdgCode();
1638
1639 fTreeVariablePIDPositive = lPIDPositive;
1640 fTreeVariablePIDNegative = lPIDNegative;
1641
1642 Int_t lblMotherPosV0Dghter = mcPosV0Dghter->GetFirstMother() ;
1643 Int_t lblMotherNegV0Dghter = mcNegV0Dghter->GetFirstMother();
1644
1645 if( lblMotherPosV0Dghter == lblMotherNegV0Dghter && lblMotherPosV0Dghter > -1 ){
1646 //either label is fine, they're equal at this stage
1647 TParticle* pThisV0 = lMCstack->Particle( lblMotherPosV0Dghter );
1648 //Set tree variables
1649 fTreeVariablePID = pThisV0->GetPdgCode(); //PDG Code
1650 fTreeVariablePtMC = pThisV0->Pt(); //Perfect Pt
1651
1652 //Only Interested if it's a Lambda, AntiLambda or K0s
1653 //Avoid the Junction Bug! PYTHIA has particles with Px=Py=Pz=E=0 occasionally,
1654 //having particle code 88 (unrecognized by PDG), for documentation purposes.
1655 //Even ROOT's TParticle::Y() is not prepared to deal with that exception!
1656 //Note that TParticle::Pt() is immune (that would just return 0)...
1657 //Though granted that that should be extremely rare in this precise condition...
1658 if( TMath::Abs(fTreeVariablePID) == 3122 || fTreeVariablePID==310 ){
1659 fTreeVariableRapMC = pThisV0->Y(); //Perfect Y
1660 }
1661 if( lMCstack->IsPhysicalPrimary (lblMotherPosV0Dghter) ) fTreeVariablePrimaryStatus = 1; //Is Primary!
1662 if( lMCstack->IsSecondaryFromWeakDecay(lblMotherPosV0Dghter) ) fTreeVariablePrimaryStatus = 2; //Weak Decay!
1663 if( lMCstack->IsSecondaryFromMaterial (lblMotherPosV0Dghter) ) fTreeVariablePrimaryStatus = 3; //Material Int!
1664
1665 //Now we try to acquire the V0 parent particle, if possible
1666 Int_t lblThisV0Parent = pThisV0->GetFirstMother();
1667 if ( lblThisV0Parent > -1 ){ //if it has a parent, get it and store specs
1668 TParticle* pThisV0Parent = lMCstack->Particle( lblThisV0Parent );
1669 fTreeVariablePIDMother = pThisV0Parent->GetPdgCode(); //V0 Mother PDG
1670 fTreeVariablePtMother = pThisV0Parent->Pt(); //V0 Mother Pt
1671 //Primary Status for the V0 Mother particle
1672 if( lMCstack->IsPhysicalPrimary (lblThisV0Parent) ) fTreeVariablePrimaryStatusMother = 1; //Is Primary!
1673 if( lMCstack->IsSecondaryFromWeakDecay(lblThisV0Parent) ) fTreeVariablePrimaryStatusMother = 2; //Weak Decay!
1674 if( lMCstack->IsSecondaryFromMaterial (lblThisV0Parent) ) fTreeVariablePrimaryStatusMother = 3; //Material Int!
1675 }
1676 }
1677
1678 fTreeVariablePt = v0->Pt();
1679 fTreeVariableChi2V0 = lChi2V0;
1680 fTreeVariableDcaV0ToPrimVertex = lDcaV0ToPrimVertex;
1681 fTreeVariableDcaV0Daughters = lDcaV0Daughters;
1682 fTreeVariableV0CosineOfPointingAngle = lV0CosineOfPointingAngle;
1683 fTreeVariableV0Radius = lV0Radius;
1684 fTreeVariableDcaPosToPrimVertex = lDcaPosToPrimVertex;
1685 fTreeVariableDcaNegToPrimVertex = lDcaNegToPrimVertex;
1686 fTreeVariableInvMassK0s = lInvMassK0s;
1687 fTreeVariableInvMassLambda = lInvMassLambda;
1688 fTreeVariableInvMassAntiLambda = lInvMassAntiLambda;
1689 fTreeVariableRapK0Short = lRapK0Short;
1690 fTreeVariableRapLambda = lRapLambda;
1691 fTreeVariableAlphaV0 = lAlphaV0;
1692 fTreeVariablePtArmV0 = lPtArmV0;
1693
1694 //Official means of acquiring N-sigmas
1695 fTreeVariableNSigmasPosProton = fPIDResponse->NumberOfSigmasTPC( pTrack, AliPID::kProton );
1696 fTreeVariableNSigmasPosPion = fPIDResponse->NumberOfSigmasTPC( pTrack, AliPID::kPion );
1697 fTreeVariableNSigmasNegProton = fPIDResponse->NumberOfSigmasTPC( nTrack, AliPID::kProton );
1698 fTreeVariableNSigmasNegPion = fPIDResponse->NumberOfSigmasTPC( nTrack, AliPID::kPion );
1699
1700//This requires an Invariant Mass Hypothesis afterwards
1701 fTreeVariableDistOverTotMom = TMath::Sqrt(
1702 TMath::Power( tDecayVertexV0[0] - lBestPrimaryVtxPos[0] , 2) +
1703 TMath::Power( tDecayVertexV0[1] - lBestPrimaryVtxPos[1] , 2) +
1704 TMath::Power( tDecayVertexV0[2] - lBestPrimaryVtxPos[2] , 2)
1705 );
1706 fTreeVariableDistOverTotMom /= (lV0TotalMomentum+1e-10); //avoid division by zero, to be sure
1707
1708 //Copy Multiplicity information
1709 fTreeVariableCentV0A = fCentrality_V0A;
1710 fTreeVariableCentV0C = fCentrality_V0C;
1711 fTreeVariableCentV0M = fCentrality_V0M;
1712 fTreeVariableCentV0AEq = fCentrality_V0AEq;
1713 fTreeVariableCentV0CEq = fCentrality_V0CEq;
1714 fTreeVariableCentV0MEq = fCentrality_V0MEq;
1715 fTreeVariableAmpV0A = fAmplitude_V0A;
1716 fTreeVariableAmpV0C = fAmplitude_V0C;
1717 fTreeVariableAmpV0AEq = fAmplitude_V0AEq;
1718 fTreeVariableAmpV0CEq = fAmplitude_V0CEq;
1719 fTreeVariableRefMultEta8 = fRefMultEta8;
1720 fTreeVariableRefMultEta5 = fRefMultEta5;
1721 fTreeVariableRunNumber = fRunNumber;
1722
1723//------------------------------------------------
1724// Fill Tree!
1725//------------------------------------------------
1726
1727 // The conditionals are meant to decrease excessive
1728 // memory usage!
1729
1730 //First Selection: Reject OnFly
1731 if( lOnFlyStatus == 0 ){
1732 //Second Selection: rough 20-sigma band, parametric.
1733 //K0Short: Enough to parametrize peak broadening with linear function.
1734 Double_t lUpperLimitK0Short = (5.63707e-01) + (1.14979e-02)*fTreeVariablePt;
1735 Double_t lLowerLimitK0Short = (4.30006e-01) - (1.10029e-02)*fTreeVariablePt;
1736 //Lambda: Linear (for higher pt) plus exponential (for low-pt broadening)
1737 //[0]+[1]*x+[2]*TMath::Exp(-[3]*x)
1738 Double_t lUpperLimitLambda = (1.13688e+00) + (5.27838e-03)*fTreeVariablePt + (8.42220e-02)*TMath::Exp(-(3.80595e+00)*fTreeVariablePt);
1739 Double_t lLowerLimitLambda = (1.09501e+00) - (5.23272e-03)*fTreeVariablePt - (7.52690e-02)*TMath::Exp(-(3.46339e+00)*fTreeVariablePt);
1740 //Do Selection
1741 if( (fTreeVariableInvMassLambda < lUpperLimitLambda && fTreeVariableInvMassLambda > lLowerLimitLambda ) ||
1742 (fTreeVariableInvMassAntiLambda < lUpperLimitLambda && fTreeVariableInvMassAntiLambda > lLowerLimitLambda ) ||
1743 (fTreeVariableInvMassK0s < lUpperLimitK0Short && fTreeVariableInvMassK0s > lLowerLimitK0Short ) ){
1744 //Pre-selection in case this is AA...
1745 if ( TMath::Abs(fTreeVariableNegEta)<0.8 && TMath::Abs(fTreeVariablePosEta)<0.8 && fkSaveV0Tree ) fTreeV0->Fill();
1746 }
1747 }
1748 }// This is the end of the V0 loop
1749
1750//------------------------------------------------
1751// Fill V0 tree over.
1752//------------------------------------------------
1753
af0d4791 1754//------------------------------------------------
1755// Rerun cascade vertexer!
1756//------------------------------------------------
1757
1758 if( fkRunVertexers ){
1759 lESDevent->ResetCascades();
1760 lESDevent->ResetV0s();
1761
1762 AliV0vertexer lV0vtxer;
1763 AliCascadeVertexer lCascVtxer;
1764
1765 lV0vtxer.SetDefaultCuts(fV0VertexerSels);
1766 lCascVtxer.SetDefaultCuts(fCascadeVertexerSels);
1767
1768 lV0vtxer.Tracks2V0vertices(lESDevent);
1769 lCascVtxer.V0sTracks2CascadeVertices(lESDevent);
1770 }
1771
1772//------------------------------------------------
1773// MAIN CASCADE LOOP STARTS HERE
1774//------------------------------------------------
1775// Code Credit: Antonin Maire (thanks^100)
1776// ---> This is an adaptation
1777
1778 Long_t ncascades = 0;
1779 ncascades = lESDevent->GetNumberOfCascades();
1780
1781 for (Int_t iXi = 0; iXi < ncascades; iXi++){
1782 //------------------------------------------------
1783 // Initializations
1784 //------------------------------------------------
1785 //Double_t lTrkgPrimaryVtxRadius3D = -500.0;
1786 //Double_t lBestPrimaryVtxRadius3D = -500.0;
1787
1788 // - 1st part of initialisation : variables needed to store AliESDCascade data members
1789 Double_t lEffMassXi = 0. ;
1790 //Double_t lChi2Xi = -1. ;
1791 Double_t lDcaXiDaughters = -1. ;
1792 Double_t lXiCosineOfPointingAngle = -1. ;
1793 Double_t lPosXi[3] = { -1000.0, -1000.0, -1000.0 };
1794 Double_t lXiRadius = -1000. ;
1795
1796 // - 2nd part of initialisation : Nbr of clusters within TPC for the 3 daughter cascade tracks
1797 Int_t lPosTPCClusters = -1; // For ESD only ...//FIXME : wait for availability in AOD
1798 Int_t lNegTPCClusters = -1; // For ESD only ...
1799 Int_t lBachTPCClusters = -1; // For ESD only ...
1800
1801 // - 3rd part of initialisation : about V0 part in cascades
1802 Double_t lInvMassLambdaAsCascDghter = 0.;
1803 //Double_t lV0Chi2Xi = -1. ;
1804 Double_t lDcaV0DaughtersXi = -1.;
1805
1806 Double_t lDcaBachToPrimVertexXi = -1., lDcaV0ToPrimVertexXi = -1.;
1807 Double_t lDcaPosToPrimVertexXi = -1.;
1808 Double_t lDcaNegToPrimVertexXi = -1.;
1809 Double_t lV0CosineOfPointingAngleXi = -1. ;
1810 Double_t lV0CosineOfPointingAngleXiSpecial = -1. ;
1811 Double_t lPosV0Xi[3] = { -1000. , -1000., -1000. }; // Position of VO coming from cascade
1812 Double_t lV0RadiusXi = -1000.0;
1813 Double_t lV0quality = 0.;
1814
1815 // - 4th part of initialisation : Effective masses
1816 Double_t lInvMassXiMinus = 0.;
1817 Double_t lInvMassXiPlus = 0.;
1818 Double_t lInvMassOmegaMinus = 0.;
1819 Double_t lInvMassOmegaPlus = 0.;
1820
1821 // - 6th part of initialisation : extra info for QA
1822 Double_t lXiMomX = 0. , lXiMomY = 0., lXiMomZ = 0.;
1823 Double_t lXiTransvMom = 0. ;
1824 //Double_t lXiTransvMomMC= 0. ;
1825 Double_t lXiTotMom = 0. ;
1826
1827 Double_t lBachMomX = 0., lBachMomY = 0., lBachMomZ = 0.;
1828 //Double_t lBachTransvMom = 0.;
1829 //Double_t lBachTotMom = 0.;
1830
1831 fTreeCascVarNegNSigmaPion = -100;
1832 fTreeCascVarNegNSigmaProton = -100;
1833 fTreeCascVarPosNSigmaPion = -100;
1834 fTreeCascVarPosNSigmaProton = -100;
1835 fTreeCascVarBachNSigmaPion = -100;
1836 fTreeCascVarBachNSigmaKaon = -100;
1837
1838 Short_t lChargeXi = -2;
1839 //Double_t lV0toXiCosineOfPointingAngle = 0. ;
1840
1841 Double_t lRapXi = -20.0, lRapOmega = -20.0, lRapMC = -20;// lEta = -20.0, lTheta = 360., lPhi = 720. ;
1842 //Double_t lAlphaXi = -200., lPtArmXi = -200.0;
1843
1844 // -------------------------------------
1845 // II.ESD - Calculation Part dedicated to Xi vertices (ESD)
1846
1847 AliESDcascade *xi = lESDevent->GetCascade(iXi);
1848 if (!xi) continue;
1849
1850 // - II.Step 2 : Assigning the necessary variables for specific AliESDcascade data members (ESD)
1851 //-------------
1852 lV0quality = 0.;
1853 xi->ChangeMassHypothesis(lV0quality , 3312); // default working hypothesis : cascade = Xi- decay
1854
1855 lEffMassXi = xi->GetEffMassXi();
1856 //lChi2Xi = xi->GetChi2Xi();
1857 lDcaXiDaughters = xi->GetDcaXiDaughters();
1858 lXiCosineOfPointingAngle = xi->GetCascadeCosineOfPointingAngle( lBestPrimaryVtxPos[0],
1859 lBestPrimaryVtxPos[1],
1860 lBestPrimaryVtxPos[2] );
1861 // Take care : the best available vertex should be used (like in AliCascadeVertexer)
1862
1863 xi->GetXYZcascade( lPosXi[0], lPosXi[1], lPosXi[2] );
1864 lXiRadius = TMath::Sqrt( lPosXi[0]*lPosXi[0] + lPosXi[1]*lPosXi[1] );
1865
1866 // - II.Step 3 : around the tracks : Bach + V0 (ESD)
1867 // ~ Necessary variables for ESDcascade data members coming from the ESDv0 part (inheritance)
1868 //-------------
1869
1870 UInt_t lIdxPosXi = (UInt_t) TMath::Abs( xi->GetPindex() );
1871 UInt_t lIdxNegXi = (UInt_t) TMath::Abs( xi->GetNindex() );
1872 UInt_t lBachIdx = (UInt_t) TMath::Abs( xi->GetBindex() );
1873 // Care track label can be negative in MC production (linked with the track quality)
1874 // However = normally, not the case for track index ...
1875
1876 // FIXME : rejection of a double use of a daughter track (nothing but just a crosscheck of what is done in the cascade vertexer)
1877 if(lBachIdx == lIdxNegXi) {
1878 AliWarning("Pb / Idx(Bach. track) = Idx(Neg. track) ... continue!"); continue;
1879 }
1880 if(lBachIdx == lIdxPosXi) {
1881 AliWarning("Pb / Idx(Bach. track) = Idx(Pos. track) ... continue!"); continue;
1882 }
1883
1884 AliESDtrack *pTrackXi = lESDevent->GetTrack( lIdxPosXi );
1885 AliESDtrack *nTrackXi = lESDevent->GetTrack( lIdxNegXi );
1886 AliESDtrack *bachTrackXi = lESDevent->GetTrack( lBachIdx );
1887
1888 if (!pTrackXi || !nTrackXi || !bachTrackXi ) {
1889 AliWarning("ERROR: Could not retrieve one of the 3 ESD daughter tracks of the cascade ...");
1890 continue;
1891 }
1892
1893 fTreeCascVarPosEta = pTrackXi->Eta();
1894 fTreeCascVarNegEta = nTrackXi->Eta();
1895 fTreeCascVarBachEta = bachTrackXi->Eta();
1896
1897 Double_t lBMom[3], lNMom[3], lPMom[3];
1898 xi->GetBPxPyPz( lBMom[0], lBMom[1], lBMom[2] );
1899 xi->GetPPxPyPz( lPMom[0], lPMom[1], lPMom[2] );
1900 xi->GetNPxPyPz( lNMom[0], lNMom[1], lNMom[2] );
1901
1902 //fTreeCascVarBachTransMom = TMath::Sqrt( lBMom[0]*lBMom[0] + lBMom[1]*lBMom[1] );
1903 //fTreeCascVarPosTransMom = TMath::Sqrt( lPMom[0]*lPMom[0] + lPMom[1]*lPMom[1] );
1904 //fTreeCascVarNegTransMom = TMath::Sqrt( lNMom[0]*lNMom[0] + lNMom[1]*lNMom[1] );
1905
1906 //------------------------------------------------
1907 // TPC dEdx information
1908 //------------------------------------------------
1909 fTreeCascVarNegNSigmaPion = fPIDResponse->NumberOfSigmasTPC( nTrackXi, AliPID::kPion );
1910 fTreeCascVarNegNSigmaProton = fPIDResponse->NumberOfSigmasTPC( nTrackXi, AliPID::kProton );
1911 fTreeCascVarPosNSigmaPion = fPIDResponse->NumberOfSigmasTPC( pTrackXi, AliPID::kPion );
1912 fTreeCascVarPosNSigmaProton = fPIDResponse->NumberOfSigmasTPC( pTrackXi, AliPID::kProton );
1913 fTreeCascVarBachNSigmaPion = fPIDResponse->NumberOfSigmasTPC( bachTrackXi, AliPID::kPion );
1914 fTreeCascVarBachNSigmaKaon = fPIDResponse->NumberOfSigmasTPC( bachTrackXi, AliPID::kKaon );
1915
1916 //------------------------------------------------
1917 // TPC Number of clusters info
1918 // --- modified to save the smallest number
1919 // --- of TPC clusters for the 3 tracks
1920 //------------------------------------------------
1921
1922 lPosTPCClusters = pTrackXi->GetTPCNcls();
1923 lNegTPCClusters = nTrackXi->GetTPCNcls();
1924 lBachTPCClusters = bachTrackXi->GetTPCNcls();
1925
1926 // 1 - Poor quality related to TPCrefit
1927 ULong_t pStatus = pTrackXi->GetStatus();
1928 ULong_t nStatus = nTrackXi->GetStatus();
1929 ULong_t bachStatus = bachTrackXi->GetStatus();
1930
1931 //fTreeCascVarkITSRefitBachelor = kTRUE;
1932 //fTreeCascVarkITSRefitNegative = kTRUE;
1933 //fTreeCascVarkITSRefitPositive = kTRUE;
1934
1935 if ((pStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Pos. track has no TPCrefit ... continue!"); continue; }
1936 if ((nStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / V0 Neg. track has no TPCrefit ... continue!"); continue; }
1937 if ((bachStatus&AliESDtrack::kTPCrefit) == 0) { AliWarning("Pb / Bach. track has no TPCrefit ... continue!"); continue; }
1938
1939 //Extra Debug Information: booleans for ITS refit
1940 //if ((pStatus&AliESDtrack::kITSrefit) == 0) { fTreeCascVarkITSRefitPositive = kFALSE; }
1941 //if ((nStatus&AliESDtrack::kITSrefit) == 0) { fTreeCascVarkITSRefitNegative = kFALSE; }
1942 //if ((bachStatus&AliESDtrack::kITSrefit) == 0) { fTreeCascVarkITSRefitBachelor = kFALSE; }
1943
1944 // 2 - Poor quality related to TPC clusters: lowest cut of 70 clusters
1945 if(lPosTPCClusters < 70) { AliWarning("Pb / V0 Pos. track has less than 70 TPC clusters ... continue!"); continue; }
1946 if(lNegTPCClusters < 70) { AliWarning("Pb / V0 Neg. track has less than 70 TPC clusters ... continue!"); continue; }
1947 if(lBachTPCClusters < 70) { AliWarning("Pb / Bach. track has less than 70 TPC clusters ... continue!"); continue; }
1948 Int_t leastnumberofclusters = 1000;
1949 if( lPosTPCClusters < leastnumberofclusters ) leastnumberofclusters = lPosTPCClusters;
1950 if( lNegTPCClusters < leastnumberofclusters ) leastnumberofclusters = lNegTPCClusters;
1951 if( lBachTPCClusters < leastnumberofclusters ) leastnumberofclusters = lBachTPCClusters;
1952
1953 lInvMassLambdaAsCascDghter = xi->GetEffMass();
1954 // This value shouldn't change, whatever the working hyp. is : Xi-, Xi+, Omega-, Omega+
1955 lDcaV0DaughtersXi = xi->GetDcaV0Daughters();
1956 //lV0Chi2Xi = xi->GetChi2V0();
1957
1958 lV0CosineOfPointingAngleXi = xi->GetV0CosineOfPointingAngle( lBestPrimaryVtxPos[0],
1959 lBestPrimaryVtxPos[1],
1960 lBestPrimaryVtxPos[2] );
1961 //Modification: V0 CosPA wrt to Cascade decay vertex
1962 lV0CosineOfPointingAngleXiSpecial = xi->GetV0CosineOfPointingAngle( lPosXi[0],
1963 lPosXi[1],
1964 lPosXi[2] );
1965
1966 lDcaV0ToPrimVertexXi = xi->GetD( lBestPrimaryVtxPos[0],
1967 lBestPrimaryVtxPos[1],
1968 lBestPrimaryVtxPos[2] );
1969
1970 lDcaBachToPrimVertexXi = TMath::Abs( bachTrackXi->GetD( lBestPrimaryVtxPos[0],
1971 lBestPrimaryVtxPos[1],
1972 lMagneticField ) );
1973 // Note : AliExternalTrackParam::GetD returns an algebraic value ...
1974
1975 xi->GetXYZ( lPosV0Xi[0], lPosV0Xi[1], lPosV0Xi[2] );
1976 lV0RadiusXi = TMath::Sqrt( lPosV0Xi[0]*lPosV0Xi[0] + lPosV0Xi[1]*lPosV0Xi[1] );
1977
1978 lDcaPosToPrimVertexXi = TMath::Abs( pTrackXi ->GetD( lBestPrimaryVtxPos[0],
1979 lBestPrimaryVtxPos[1],
1980 lMagneticField ) );
1981
1982 lDcaNegToPrimVertexXi = TMath::Abs( nTrackXi ->GetD( lBestPrimaryVtxPos[0],
1983 lBestPrimaryVtxPos[1],
1984 lMagneticField ) );
1985
1986 // - II.Step 4 : around effective masses (ESD)
1987 // ~ change mass hypotheses to cover all the possibilities : Xi-/+, Omega -/+
1988
1989 if( bachTrackXi->Charge() < 0 ) {
1990 lV0quality = 0.;
1991 xi->ChangeMassHypothesis(lV0quality , 3312);
1992 // Calculate the effective mass of the Xi- candidate.
1993 // pdg code 3312 = Xi-
1994 lInvMassXiMinus = xi->GetEffMassXi();
1995
1996 lV0quality = 0.;
1997 xi->ChangeMassHypothesis(lV0quality , 3334);
1998 // Calculate the effective mass of the Xi- candidate.
1999 // pdg code 3334 = Omega-
2000 lInvMassOmegaMinus = xi->GetEffMassXi();
2001
2002 lV0quality = 0.;
2003 xi->ChangeMassHypothesis(lV0quality , 3312); // Back to default hyp.
2004 }// end if negative bachelor
2005
2006
2007 if( bachTrackXi->Charge() > 0 ){
2008 lV0quality = 0.;
2009 xi->ChangeMassHypothesis(lV0quality , -3312);
2010 // Calculate the effective mass of the Xi+ candidate.
2011 // pdg code -3312 = Xi+
2012 lInvMassXiPlus = xi->GetEffMassXi();
2013
2014 lV0quality = 0.;
2015 xi->ChangeMassHypothesis(lV0quality , -3334);
2016 // Calculate the effective mass of the Xi+ candidate.
2017 // pdg code -3334 = Omega+
2018 lInvMassOmegaPlus = xi->GetEffMassXi();
2019
2020 lV0quality = 0.;
2021 xi->ChangeMassHypothesis(lV0quality , -3312); // Back to "default" hyp.
2022 }// end if positive bachelor
2023 // - II.Step 6 : extra info for QA (ESD)
2024 // miscellaneous pieces of info that may help regarding data quality assessment.
2025 //-------------
2026
2027 xi->GetPxPyPz( lXiMomX, lXiMomY, lXiMomZ );
2028 lXiTransvMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY );
2029 lXiTotMom = TMath::Sqrt( lXiMomX*lXiMomX + lXiMomY*lXiMomY + lXiMomZ*lXiMomZ );
2030
2031 xi->GetBPxPyPz( lBachMomX, lBachMomY, lBachMomZ );
2032 //lBachTransvMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY );
2033 //lBachTotMom = TMath::Sqrt( lBachMomX*lBachMomX + lBachMomY*lBachMomY + lBachMomZ*lBachMomZ );
2034
2035 lChargeXi = xi->Charge();
2036
2037 //lV0toXiCosineOfPointingAngle = xi->GetV0CosineOfPointingAngle( lPosXi[0], lPosXi[1], lPosXi[2] );
2038
2039 lRapXi = xi->RapXi();
2040 lRapOmega = xi->RapOmega();
2041 //lEta = xi->Eta();
2042 //lTheta = xi->Theta() *180.0/TMath::Pi();
2043 //lPhi = xi->Phi() *180.0/TMath::Pi();
2044 //lAlphaXi = xi->AlphaXi();
2045 //lPtArmXi = xi->PtArmXi();
2046
2047//------------------------------------------------
2048// Associate Cascade Candidates to Monte Carlo!
2049//------------------------------------------------
2050
2051//Warning: Not using Continues... Need to fill tree later!
2052
2053 Double_t lXiTransvMomMC= 0. ;
2054 Int_t lPDGCodeCascade = 0;
2055 Int_t lPID_BachMother = 0;
2056 Int_t lPID_NegMother = 0;
2057 Int_t lPID_PosMother = 0;
2058 fTreeCascVarIsPhysicalPrimary = 0; // 0: not defined, any candidate may have this
2059
2060 if(fDebug > 5)
2061 cout << "MC EventNumber : " << lMCevent->Header()->GetEvent()
2062 << " / MC event Number in Run : " << lMCevent->Header()->GetEventNrInRun() << endl;
2063
2064//----------------------------------------
2065// Regular MC ASSOCIATION STARTS HERE
2066//----------------------------------------
2067
2068 Int_t lblPosV0Dghter = (Int_t) TMath::Abs( pTrackXi->GetLabel() );
2069 // Abs value = needed ! question of quality track association ...
2070 Int_t lblNegV0Dghter = (Int_t) TMath::Abs( nTrackXi->GetLabel() );
2071 Int_t lblBach = (Int_t) TMath::Abs( bachTrackXi->GetLabel() );
2072
2073 TParticle* mcPosV0Dghter = lMCstack->Particle( lblPosV0Dghter );
2074 TParticle* mcNegV0Dghter = lMCstack->Particle( lblNegV0Dghter );
2075 TParticle* mcBach = lMCstack->Particle( lblBach );
2076
2077 //fTreeCascVarPosTransMomMC = mcPosV0Dghter->Pt();
2078 //fTreeCascVarNegTransMomMC = mcNegV0Dghter->Pt();
2079
2080 //fTreeCascVarPIDPositive = mcPosV0Dghter -> GetPdgCode();
2081 //fTreeCascVarPIDNegative = mcNegV0Dghter -> GetPdgCode();
2082 //fTreeCascVarPIDBachelor = mcBach->GetPdgCode();
2083
2084 // - Step 4.2 : level of the Xi daughters
2085
2086 Int_t lblMotherPosV0Dghter = mcPosV0Dghter->GetFirstMother() ;
2087 Int_t lblMotherNegV0Dghter = mcNegV0Dghter->GetFirstMother();
2088
2089 //Rather uncivilized: Open brackets for each 'continue'
2090 if(! (lblMotherPosV0Dghter != lblMotherNegV0Dghter) ) { // same mother
2091 if(! (lblMotherPosV0Dghter < 0) ) { // mother != primary (!= -1)
2092 if(! (lblMotherNegV0Dghter < 0) ) {
2093
2094 // mothers = Lambda candidate ... a priori
2095
2096 TParticle* mcMotherPosV0Dghter = lMCstack->Particle( lblMotherPosV0Dghter );
2097 TParticle* mcMotherNegV0Dghter = lMCstack->Particle( lblMotherNegV0Dghter );
2098
2099 // - Step 4.3 : level of Xi candidate
2100
2101 Int_t lblGdMotherPosV0Dghter = mcMotherPosV0Dghter->GetFirstMother() ;
2102 Int_t lblGdMotherNegV0Dghter = mcMotherNegV0Dghter->GetFirstMother() ;
2103
2104 if(! (lblGdMotherPosV0Dghter != lblGdMotherNegV0Dghter) ) {
2105 if(! (lblGdMotherPosV0Dghter < 0) ) { // primary lambda ...
2106 if(! (lblGdMotherNegV0Dghter < 0) ) { // primary lambda ...
2107
2108 // Gd mothers = Xi candidate ... a priori
2109
2110 TParticle* mcGdMotherPosV0Dghter = lMCstack->Particle( lblGdMotherPosV0Dghter );
2111 TParticle* mcGdMotherNegV0Dghter = lMCstack->Particle( lblGdMotherNegV0Dghter );
2112
2113 Int_t lblMotherBach = (Int_t) TMath::Abs( mcBach->GetFirstMother() );
2114
2115 // if( lblMotherBach != lblGdMotherPosV0Dghter ) continue; //same mother for bach and V0 daughters
2116 if(!(lblMotherBach != lblGdMotherPosV0Dghter)) { //same mother for bach and V0 daughters
2117
2118 TParticle* mcMotherBach = lMCstack->Particle( lblMotherBach );
2119
2120 lPID_BachMother = mcMotherBach->GetPdgCode();
2121 lPID_NegMother = mcGdMotherPosV0Dghter->GetPdgCode();
2122 lPID_PosMother = mcGdMotherNegV0Dghter->GetPdgCode();
2123
2124 if(lPID_BachMother==lPID_NegMother && lPID_BachMother==lPID_PosMother){
2125 lPDGCodeCascade = lPID_BachMother;
2126 lXiTransvMomMC = mcMotherBach->Pt();
2127 if( lMCstack->IsPhysicalPrimary (lblMotherBach) ) fTreeCascVarIsPhysicalPrimary = 1; //Is Primary!
2128 if( lMCstack->IsSecondaryFromWeakDecay(lblMotherBach) ) fTreeCascVarIsPhysicalPrimary = 2; //Weak Decay!
2129 if( lMCstack->IsSecondaryFromMaterial (lblMotherBach) ) fTreeCascVarIsPhysicalPrimary = 3; //From Material!
2130 if ( (mcMotherBach->Energy() + mcMotherBach->Pz()) / (mcMotherBach->Energy() - mcMotherBach->Pz() +1.e-13) !=0 ){
2131 lRapMC = 0.5*TMath::Log( (mcMotherBach->Energy() + mcMotherBach->Pz()) / (mcMotherBach->Energy() - mcMotherBach->Pz() +1.e-13) );
2132 }
2133 }
2134
2135 }}}}}}} //Ends all conditionals above...
2136
2137//----------------------------------------
2138// Regular MC ASSOCIATION ENDS HERE
2139//----------------------------------------
2140
2141
2142 //------------------------------------------------
2143 // Set Variables for adding to tree
2144 //------------------------------------------------
2145
2146 fTreeCascVarCharge = lChargeXi;
2147 fTreeCascVarPID = lPDGCodeCascade;
2148 if(lInvMassXiMinus!=0) fTreeCascVarMassAsXi = lInvMassXiMinus;
2149 if(lInvMassXiPlus!=0) fTreeCascVarMassAsXi = lInvMassXiPlus;
2150 if(lInvMassOmegaMinus!=0) fTreeCascVarMassAsOmega = lInvMassOmegaMinus;
2151 if(lInvMassOmegaPlus!=0) fTreeCascVarMassAsOmega = lInvMassOmegaPlus;
2152 fTreeCascVarPt = lXiTransvMom;
2153 fTreeCascVarPtMC = lXiTransvMomMC;
2154 fTreeCascVarRapXi = lRapXi ;
2155 fTreeCascVarRapMC = lRapMC ;
2156 fTreeCascVarRapOmega = lRapOmega ;
2157 fTreeCascVarDCACascDaughters = lDcaXiDaughters;
2158 fTreeCascVarDCABachToPrimVtx = lDcaBachToPrimVertexXi;
2159 fTreeCascVarDCAV0Daughters = lDcaV0DaughtersXi;
2160 fTreeCascVarDCAV0ToPrimVtx = lDcaV0ToPrimVertexXi;
2161 fTreeCascVarDCAPosToPrimVtx = lDcaPosToPrimVertexXi;
2162 fTreeCascVarDCANegToPrimVtx = lDcaNegToPrimVertexXi;
2163 fTreeCascVarCascCosPointingAngle = lXiCosineOfPointingAngle;
2164 fTreeCascVarCascRadius = lXiRadius;
2165 fTreeCascVarV0Mass = lInvMassLambdaAsCascDghter;
2166 fTreeCascVarV0CosPointingAngle = lV0CosineOfPointingAngleXi;
2167 fTreeCascVarV0CosPointingAngleSpecial = lV0CosineOfPointingAngleXiSpecial;
2168 fTreeCascVarV0Radius = lV0RadiusXi;
2169 fTreeCascVarLeastNbrClusters = leastnumberofclusters;
2170
2171 //Copy Multiplicity information
2172 fTreeCascVarCentV0A = fCentrality_V0A;
2173 fTreeCascVarCentV0C = fCentrality_V0C;
2174 fTreeCascVarCentV0M = fCentrality_V0M;
2175 fTreeCascVarCentV0AEq = fCentrality_V0AEq;
2176 fTreeCascVarCentV0CEq = fCentrality_V0CEq;
2177 fTreeCascVarCentV0MEq = fCentrality_V0MEq;
2178 fTreeCascVarAmpV0A = fAmplitude_V0A;
2179 fTreeCascVarAmpV0C = fAmplitude_V0C;
2180 fTreeCascVarAmpV0AEq = fAmplitude_V0AEq;
2181 fTreeCascVarAmpV0CEq = fAmplitude_V0CEq;
2182 fTreeCascVarRefMultEta8 = fRefMultEta8;
2183 fTreeCascVarRefMultEta5 = fRefMultEta5;
26d39641 2184 fTreeCascVarRunNumber = fRunNumber;
af0d4791 2185
2186 fTreeCascVarDistOverTotMom = TMath::Sqrt(
2187 TMath::Power( lPosXi[0] - lBestPrimaryVtxPos[0] , 2) +
2188 TMath::Power( lPosXi[1] - lBestPrimaryVtxPos[1] , 2) +
2189 TMath::Power( lPosXi[2] - lBestPrimaryVtxPos[2] , 2)
2190 );
2191 fTreeCascVarDistOverTotMom /= (lXiTotMom+1e-13);
2192
2193//All vars not specified here: specified elsewhere!
2194
2195//------------------------------------------------
2196// Fill Tree!
2197//------------------------------------------------
2198
2199// The conditional is meant to decrease excessive
2200// memory usage! Be careful when loosening the
2201// cut!
2202
2203 //Xi Mass window: 150MeV wide
2204 //Omega mass window: 150MeV wide
2205
2206 if( fkSaveCascadeTree && ( (fTreeCascVarMassAsXi<1.32+0.075&&fTreeCascVarMassAsXi>1.32-0.075) ||
2207 (fTreeCascVarMassAsOmega<1.68+0.075&&fTreeCascVarMassAsOmega>1.68-0.075) ) ){
2208 fTreeCascade->Fill();
2209 }
2210
2211//------------------------------------------------
2212// Fill tree over.
2213//------------------------------------------------
2214
2215 }// end of the Cascade loop (ESD or AOD)
2216
2217 // Post output data.
2218 PostData(1, fListHist);
2219 PostData(2, fTreeEvent);
2220 PostData(3, fTreeV0);
2221 PostData(4, fTreeCascade);
2222}
2223
2224//________________________________________________________________________
2225void AliAnalysisTaskStrangenessVsMultiplicityMC::Terminate(Option_t *)
2226{
2227 // Draw result to the screen
2228 // Called once at the end of the query
2229
2230 TList *cRetrievedList = 0x0;
2231 cRetrievedList = (TList*)GetOutputData(1);
2232 if(!cRetrievedList){
2233 Printf("ERROR - AliAnalysisTaskStrangenessVsMultiplicityMC : ouput data container list not available\n");
2234 return;
2235 }
2236
2237 fHistEventCounter = dynamic_cast<TH1D*> ( cRetrievedList->FindObject("fHistEventCounter") );
2238 if (!fHistEventCounter) {
2239 Printf("ERROR - AliAnalysisTaskStrangenessVsMultiplicityMC : fHistEventCounter not available");
2240 return;
2241 }
2242
2243 TCanvas *canCheck = new TCanvas("AliAnalysisTaskStrangenessVsMultiplicityMC","V0 Multiplicity",10,10,510,510);
2244 canCheck->cd(1)->SetLogy();
2245
2246 fHistEventCounter->SetMarkerStyle(22);
2247 fHistEventCounter->DrawCopy("E");
2248}
2249
2250//----------------------------------------------------------------------------
2251
2252Double_t AliAnalysisTaskStrangenessVsMultiplicityMC::MyRapidity(Double_t rE, Double_t rPz) const
2253{
2254 // Local calculation for rapidity
2255 Double_t ReturnValue = -100;
2256 if( (rE-rPz+1.e-13) != 0 && (rE+rPz) != 0 ){
2257 ReturnValue = 0.5*TMath::Log((rE+rPz)/(rE-rPz+1.e-13));
2258 }
2259 return ReturnValue;
2260}