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