37 #include <TParticle.h>
38 #include <TParticlePDG.h>
42 #include <THnSparse.h>
43 #include <TLorentzVector.h>
46 #include <TDatabasePDG.h>
47 #include <AliAnalysisDataSlot.h>
48 #include <AliAnalysisDataContainer.h>
50 #include "AliMCEvent.h"
51 #include "AliAnalysisManager.h"
52 #include "AliAODMCHeader.h"
53 #include "AliAODHandler.h"
55 #include "AliExternalTrackParam.h"
56 #include "AliAODVertex.h"
57 #include "AliESDVertex.h"
58 #include "AliAODRecoDecay.h"
61 #include "AliESDtrack.h"
62 #include "AliAODTrack.h"
64 #include "AliAODcascade.h"
65 #include "AliAODMCParticle.h"
66 #include "AliAnalysisTaskSE.h"
68 #include "AliPIDResponse.h"
69 #include "AliPIDCombined.h"
70 #include "AliTOFPIDResponse.h"
72 #include "AliInputEventHandler.h"
73 #include "AliESDtrackCuts.h"
74 #include "AliNeutralTrackParam.h"
75 #include "AliKFParticle.h"
76 #include "AliKFVertex.h"
77 #include "AliExternalTrackParam.h"
78 #include "AliESDtrack.h"
79 #include "AliCentrality.h"
80 #include "AliVertexerTracks.h"
81 #include "AliEventPoolManager.h"
102 fIsEventSelected(kFALSE),
103 fWriteVariableTree(kFALSE),
104 fWriteEachVariableTree(kFALSE),
105 fWriteMCVariableTree(kFALSE),
107 fEleVariablesTree(0),
110 fMCEleVariablesTree(0),
111 fMCV0VariablesTree(0),
112 fMCGenPairVariablesTree(0),
113 fCorrelationVariablesTree(0),
114 fReconstructPrimVert(kFALSE),
120 fCandidateVariables(),
121 fCandidateEleVariables(),
122 fCandidateV0Variables(),
123 fCandidateMCVariables(),
124 fCandidateMCEleVariables(),
125 fCandidateMCV0Variables(),
126 fCandidateMCGenPairVariables(),
127 fCorrelationVariables(),
135 fUseCentralityV0M(kFALSE),
138 fMCDoPairAnalysis(kFALSE),
139 fHistoEleLambdaMass(0),
140 fHistoEleLambdaMassRS(0),
141 fHistoEleLambdaMassRS1(0),
142 fHistoEleLambdaMassRS2(0),
143 fHistoEleLambdaMassWS(0),
144 fHistoEleLambdaMassWS1(0),
145 fHistoEleLambdaMassWS2(0),
146 fHistoEleLambdaMassRSMix(0),
147 fHistoEleLambdaMassRSMix1(0),
148 fHistoEleLambdaMassRSMix2(0),
149 fHistoEleLambdaMassWSMix(0),
150 fHistoEleLambdaMassWSMix1(0),
151 fHistoEleLambdaMassWSMix2(0),
152 fHistoEleLambdaMassRSSide(0),
153 fHistoEleLambdaMassRSSide1(0),
154 fHistoEleLambdaMassRSSide2(0),
155 fHistoEleLambdaMassWSSide(0),
156 fHistoEleLambdaMassWSSide1(0),
157 fHistoEleLambdaMassWSSide2(0),
158 fHistoEleLambdaMassRSAway(0),
159 fHistoEleLambdaMassRS1Away(0),
160 fHistoEleLambdaMassRS2Away(0),
161 fHistoEleLambdaMassWSAway(0),
162 fHistoEleLambdaMassWS1Away(0),
163 fHistoEleLambdaMassWS2Away(0),
164 fHistoEleLambdaMassRSMixAway(0),
165 fHistoEleLambdaMassRSMix1Away(0),
166 fHistoEleLambdaMassRSMix2Away(0),
167 fHistoEleLambdaMassWSMixAway(0),
168 fHistoEleLambdaMassWSMix1Away(0),
169 fHistoEleLambdaMassWSMix2Away(0),
170 fHistoEleLambdaMassRSSideAway(0),
171 fHistoEleLambdaMassRSSide1Away(0),
172 fHistoEleLambdaMassRSSide2Away(0),
173 fHistoEleLambdaMassWSSideAway(0),
174 fHistoEleLambdaMassWSSide1Away(0),
175 fHistoEleLambdaMassWSSide2Away(0),
176 fHistoEleLambdaMassvsElePtRS(0),
177 fHistoEleLambdaMassvsElePtWS(0),
178 fHistoEleLambdaMassvsElePtRSMix(0),
179 fHistoEleLambdaMassvsElePtWSMix(0),
180 fHistoEleLambdaMassvsElePtRSSide(0),
181 fHistoEleLambdaMassvsElePtWSSide(0),
182 fHistoEleLambdaMassvsElePtRS1(0),
183 fHistoEleLambdaMassvsElePtWS1(0),
184 fHistoEleLambdaMassvsElePtRSMix1(0),
185 fHistoEleLambdaMassvsElePtWSMix1(0),
186 fHistoEleLambdaMassvsElePtRSSide1(0),
187 fHistoEleLambdaMassvsElePtWSSide1(0),
188 fHistoEleLambdaMassvsElePtRS2(0),
189 fHistoEleLambdaMassvsElePtWS2(0),
190 fHistoEleLambdaMassvsElePtRSMix2(0),
191 fHistoEleLambdaMassvsElePtWSMix2(0),
192 fHistoEleLambdaMassvsElePtRSSide2(0),
193 fHistoEleLambdaMassvsElePtWSSide2(0),
198 fHistoEleLambdaMassMCS(0),
199 fHistoEleLambdaMassLambdabMCS(0),
200 fHistoEleLambdaMassXibMCS(0),
201 fHistoEleLambdaMassMCGen(0),
202 fHistoEleLambdaMassvsElePtMCS(0),
203 fHistoEleLambdaMassvsElePtMCGen(0),
204 fHistoEleLambdaMassvsElePtMCS1(0),
205 fHistoEleLambdaMassvsElePtMCGen1(0),
206 fHistoEleLambdaMassvsElePtMCS2(0),
207 fHistoEleLambdaMassvsElePtMCGen2(0),
210 fHistoElePtvsEtaRS(0),
211 fHistoElePtvsEtaWS(0),
212 fHistoElePtvsEtaRSMix(0),
213 fHistoElePtvsEtaWSMix(0),
214 fHistoElePtvsEtaMCS(0),
215 fHistoElePtvsEtaMCGen(0),
216 fHistoElePtvsLambdaPtRS(0),
217 fHistoElePtvsLambdaPtWS(0),
218 fHistoElePtvsLambdaPtRSMix(0),
219 fHistoElePtvsLambdaPtWSMix(0),
220 fHistoElePtvsLambdaPtMCS(0),
221 fHistoElePtvsLambdaPtvsLcPtMCS(0),
222 fHistoElePtvsLambdaPtMCGen(0),
223 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
224 fHistoElePtvsLambdaPtMCLcGen(0),
225 fHistoElePtvsd0RS(0),
226 fHistoElePtvsd0WS(0),
227 fHistoElePtvsd0RSMix(0),
228 fHistoElePtvsd0WSMix(0),
229 fHistoElePtvsd0MCS(0),
230 fHistoElePtvsd0PromptMCS(0),
231 fHistoElePtvsd0BFeeddownMCS(0),
232 fHistoEleLambdaMassFeeddownXic0MCS(0),
233 fHistoEleLambdaMassFeeddownXic0MCGen(0),
234 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
235 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
236 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
237 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
238 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
239 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
240 fHistoElePtFeeddownXic0MCS(0),
241 fHistoElePtFeeddownXic0MCGen(0),
242 fHistoElePtvsEtaFeeddownXic0MCS(0),
243 fHistoElePtvsEtaFeeddownXic0MCGen(0),
244 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
245 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
246 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
247 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
248 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
249 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
250 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
251 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
252 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
253 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
254 fHistoElePtFeeddownXicPlusMCS(0),
255 fHistoElePtFeeddownXicPlusMCGen(0),
256 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
257 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
258 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
259 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
262 fHistoBachPtMCGen(0),
264 fHistoLambdaMassvsPt(0),
265 fHistoLambdaMassvsPtMCS(0),
266 fHistoLambdaMassvsPtMCGen(0),
267 fHistoLambdaPtvsEta(0),
268 fHistoLambdaPtvsEtaMCS(0),
269 fHistoLambdaPtvsEtaMCGen(0),
270 fHistoLambdaArmenteros(0),
271 fHistoLambdaArmenterosMCS(0),
272 fHistoLambdaPtvsDl(0),
273 fHistoLambdaPtvsDlSide(0),
274 fHistoLambdaPtvsDlMCS(0),
275 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
276 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
277 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
278 fHistoK0sMassvsPt(0),
279 fHistoElectronTPCPID(0),
280 fHistoElectronTOFPID(0),
281 fHistoElectronTPCSelPID(0),
282 fHistoElectronTOFSelPID(0),
283 fHistoElectronTPCPIDSelTOF(0),
284 fHistoElectronTPCPIDSelTOFSmallEta(0),
285 fHistoElectronTPCPIDSelTOFLargeEta(0),
286 fHistoMassConversionsMin(0),
287 fHistoMassConversionsSameSignMin(0),
288 fHistoElectronQovPtvsPhi(0),
289 fHistoLambdaQovPtvsPhi(0),
296 fHistoFeedDownXic0MCGen(0),
297 fHistoFeedDownXic0MCGen1(0),
298 fHistoFeedDownXic0MCGen2(0),
299 fHistoFeedDownXic0MCS(0),
300 fHistoFeedDownXic0MCS1(0),
301 fHistoFeedDownXic0MCS2(0),
302 fHistoFeedDownXicPlusMCGen(0),
303 fHistoFeedDownXicPlusMCGen1(0),
304 fHistoFeedDownXicPlusMCGen2(0),
305 fHistoFeedDownXicPlusMCS(0),
306 fHistoFeedDownXicPlusMCS1(0),
307 fHistoFeedDownXicPlusMCS2(0),
308 fHistoFeedDownSigmaMCGen(0),
309 fHistoFeedDownSigmaMCGen1(0),
310 fHistoFeedDownSigmaMCGen2(0),
311 fHistoFeedDownSigmaMCS(0),
312 fHistoFeedDownSigmaMCS1(0),
313 fHistoFeedDownSigmaMCS2(0),
314 fHistoLcElectronMCGen(0),
315 fHistoLcElectronMCGen1(0),
316 fHistoLcElectronMCGen2(0),
317 fHistoLcElectronMCS(0),
318 fHistoLcElectronMCS1(0),
319 fHistoLcElectronMCS2(0),
320 fHistoElectronFeedDownXic0MCGen(0),
321 fHistoElectronFeedDownXic0MCGen1(0),
322 fHistoElectronFeedDownXic0MCGen2(0),
323 fHistoElectronFeedDownXic0MCS(0),
324 fHistoElectronFeedDownXic0MCS1(0),
325 fHistoElectronFeedDownXic0MCS2(0),
326 fHistoElectronFeedDownXicPlusMCGen(0),
327 fHistoElectronFeedDownXicPlusMCGen1(0),
328 fHistoElectronFeedDownXicPlusMCGen2(0),
329 fHistoElectronFeedDownXicPlusMCS(0),
330 fHistoElectronFeedDownXicPlusMCS1(0),
331 fHistoElectronFeedDownXicPlusMCS2(0),
332 fHistoElectronMCGen(0),
333 fHistoLambdaMCGen(0),
334 fHistoElePtvsV0dlRS(0),
335 fHistoElePtvsV0dlRS1(0),
336 fHistoElePtvsV0dlRS2(0),
337 fHistoElePtvsV0dlRSSide(0),
338 fHistoElePtvsV0dlRSSide1(0),
339 fHistoElePtvsV0dlRSSide2(0),
340 fHistoElePtvsV0dlRSMix(0),
341 fHistoElePtvsV0dlRSMix1(0),
342 fHistoElePtvsV0dlRSMix2(0),
343 fHistoElePtvsV0dlWS(0),
344 fHistoElePtvsV0dlWS1(0),
345 fHistoElePtvsV0dlWS2(0),
346 fHistoElePtvsV0dlWSSide(0),
347 fHistoElePtvsV0dlWSSide1(0),
348 fHistoElePtvsV0dlWSSide2(0),
349 fHistoElePtvsV0dlWSMix(0),
350 fHistoElePtvsV0dlWSMix1(0),
351 fHistoElePtvsV0dlWSMix2(0),
352 fHistoElePtvsV0dlMCS(0),
353 fHistoElePtvsV0dlMCS1(0),
354 fHistoElePtvsV0dlMCS2(0),
355 fHistoElePtvsV0dlFeedDownXic0MCS(0),
356 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
357 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
358 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
359 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
360 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
361 fHistoElePtvsV0dcaRS(0),
362 fHistoElePtvsV0dcaRS1(0),
363 fHistoElePtvsV0dcaRS2(0),
364 fHistoElePtvsV0dcaRSSide(0),
365 fHistoElePtvsV0dcaRSSide1(0),
366 fHistoElePtvsV0dcaRSSide2(0),
367 fHistoElePtvsV0dcaRSMix(0),
368 fHistoElePtvsV0dcaRSMix1(0),
369 fHistoElePtvsV0dcaRSMix2(0),
370 fHistoElePtvsV0dcaWS(0),
371 fHistoElePtvsV0dcaWS1(0),
372 fHistoElePtvsV0dcaWS2(0),
373 fHistoElePtvsV0dcaWSSide(0),
374 fHistoElePtvsV0dcaWSSide1(0),
375 fHistoElePtvsV0dcaWSSide2(0),
376 fHistoElePtvsV0dcaWSMix(0),
377 fHistoElePtvsV0dcaWSMix1(0),
378 fHistoElePtvsV0dcaWSMix2(0),
379 fHistoElePtvsV0dcaMCS(0),
380 fHistoElePtvsV0dcaMCS1(0),
381 fHistoElePtvsV0dcaMCS2(0),
382 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
383 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
384 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
385 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
386 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
387 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
388 fHistoEleLambdaPtvsV0dlRS(0),
389 fHistoEleLambdaPtvsV0dlRS1(0),
390 fHistoEleLambdaPtvsV0dlRS2(0),
391 fHistoEleLambdaPtvsV0dlRSSide(0),
392 fHistoEleLambdaPtvsV0dlRSSide1(0),
393 fHistoEleLambdaPtvsV0dlRSSide2(0),
394 fHistoEleLambdaPtvsV0dlRSMix(0),
395 fHistoEleLambdaPtvsV0dlRSMix1(0),
396 fHistoEleLambdaPtvsV0dlRSMix2(0),
397 fHistoEleLambdaPtvsV0dlWS(0),
398 fHistoEleLambdaPtvsV0dlWS1(0),
399 fHistoEleLambdaPtvsV0dlWS2(0),
400 fHistoEleLambdaPtvsV0dlWSSide(0),
401 fHistoEleLambdaPtvsV0dlWSSide1(0),
402 fHistoEleLambdaPtvsV0dlWSSide2(0),
403 fHistoEleLambdaPtvsV0dlWSMix(0),
404 fHistoEleLambdaPtvsV0dlWSMix1(0),
405 fHistoEleLambdaPtvsV0dlWSMix2(0),
406 fHistoEleLambdaPtvsV0dlMCS(0),
407 fHistoEleLambdaPtvsV0dlMCS1(0),
408 fHistoEleLambdaPtvsV0dlMCS2(0),
409 fHistoEleLambdaPtvsV0dlPromptMCS(0),
410 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
411 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
412 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
413 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
414 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
415 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
416 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
417 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
418 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
419 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
420 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
421 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
422 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
423 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
424 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
425 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
426 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
427 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
428 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
429 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
430 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
431 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
432 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
433 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
434 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
435 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
436 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
437 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
438 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
439 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
440 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
441 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
442 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
443 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
444 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
445 fHistoEleLambdaPtvsV0dlXibMCS(0),
446 fHistoEleLambdaPtvsV0dlXibMCS1(0),
447 fHistoEleLambdaPtvsV0dlXibMCS2(0),
448 fHistoEleLambdaPtvsV0dlRSAway(0),
449 fHistoEleLambdaPtvsV0dlRS1Away(0),
450 fHistoEleLambdaPtvsV0dlRS2Away(0),
451 fHistoEleLambdaPtvsV0dlRSSideAway(0),
452 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
453 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
454 fHistoEleLambdaPtvsV0dlRSMixAway(0),
455 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
456 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
457 fHistoEleLambdaPtvsV0dlWSAway(0),
458 fHistoEleLambdaPtvsV0dlWS1Away(0),
459 fHistoEleLambdaPtvsV0dlWS2Away(0),
460 fHistoEleLambdaPtvsV0dlWSSideAway(0),
461 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
462 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
463 fHistoEleLambdaPtvsV0dlWSMixAway(0),
464 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
465 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
466 fHistoResponseElePt(0),
467 fHistoResponseElePt1(0),
468 fHistoResponseElePt2(0),
469 fHistoResponseEleLambdaPt(0),
470 fHistoResponseEleLambdaPt1(0),
471 fHistoResponseEleLambdaPt2(0),
472 fHistoResponseEleLambdaPtFeeddownXic0(0),
473 fHistoResponseEleLambdaPtFeeddownXic01(0),
474 fHistoResponseEleLambdaPtFeeddownXic02(0),
475 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
476 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
477 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
478 fHistoResponseEleLambdaPtFeeddownSigma(0),
479 fHistoResponseEleLambdaPtFeeddownSigma1(0),
480 fHistoResponseEleLambdaPtFeeddownSigma2(0),
481 fHistoResponseLambdaPt(0),
482 fHistoResponseLambdaPtFeeddownXic0(0),
483 fHistoResponseLambdaPtFeeddownXicPlus(0),
484 fHistoResponseLambdaPtFeeddownSigma(0),
485 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
486 fHistoEleLambdaPtvsRapidityRS(0),
487 fHistoEleLambdaPtvsRapidityWS(0),
488 fHistoEleLambdaPtvsRapidityMCS(0),
489 fHistoElectronPi0Total(0),
490 fHistoElectronPi0Tag(0),
491 fHistoElectronEtaTotal(0),
492 fHistoElectronEtaTag(0),
494 fHistonEvtvsRunNumber(0),
495 fHistonElevsRunNumber(0),
496 fHistonLambdavsRunNumber(0),
497 fHistoMCEventType(0),
498 fHistoMCDeltaPhiccbar(0),
499 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
500 fHistodPhiSdEtaSElectronProtonR125RS(0),
501 fHistodPhiSdEtaSElectronProtonR125WS(0),
502 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
503 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
504 fHistodPhiSdEtaSElectronPionR125RS(0),
505 fHistodPhiSdEtaSElectronPionR125WS(0),
506 fHistodPhiSdEtaSElectronPionR125RSMix(0),
507 fHistodPhiSdEtaSElectronPionR125WSMix(0),
509 fMixWithoutConversionFlag(kTRUE),
510 fNumberOfEventsForMixing (5),
516 fElectronTracks(0x0),
523 fElectronCutVarsArray(0x0),
524 fV0CutVarsArray1(0x0),
525 fV0CutVarsArray2(0x0)
530 for(Int_t i=0;i<17;i++){
535 for(Int_t i=0;i<8;i++){
543 Bool_t writeVariableTree) :
544 AliAnalysisTaskSE(name),
553 fIsEventSelected(kFALSE),
554 fWriteVariableTree(writeVariableTree),
555 fWriteEachVariableTree(kFALSE),
556 fWriteMCVariableTree(kFALSE),
558 fEleVariablesTree(0),
561 fMCEleVariablesTree(0),
562 fMCV0VariablesTree(0),
563 fMCGenPairVariablesTree(0),
564 fCorrelationVariablesTree(0),
565 fReconstructPrimVert(kFALSE),
571 fCandidateVariables(),
572 fCandidateEleVariables(),
573 fCandidateV0Variables(),
574 fCandidateMCVariables(),
575 fCandidateMCEleVariables(),
576 fCandidateMCV0Variables(),
577 fCandidateMCGenPairVariables(),
578 fCorrelationVariables(),
586 fUseCentralityV0M(kFALSE),
589 fMCDoPairAnalysis(kFALSE),
590 fHistoEleLambdaMass(0),
591 fHistoEleLambdaMassRS(0),
592 fHistoEleLambdaMassRS1(0),
593 fHistoEleLambdaMassRS2(0),
594 fHistoEleLambdaMassWS(0),
595 fHistoEleLambdaMassWS1(0),
596 fHistoEleLambdaMassWS2(0),
597 fHistoEleLambdaMassRSMix(0),
598 fHistoEleLambdaMassRSMix1(0),
599 fHistoEleLambdaMassRSMix2(0),
600 fHistoEleLambdaMassWSMix(0),
601 fHistoEleLambdaMassWSMix1(0),
602 fHistoEleLambdaMassWSMix2(0),
603 fHistoEleLambdaMassRSSide(0),
604 fHistoEleLambdaMassRSSide1(0),
605 fHistoEleLambdaMassRSSide2(0),
606 fHistoEleLambdaMassWSSide(0),
607 fHistoEleLambdaMassWSSide1(0),
608 fHistoEleLambdaMassWSSide2(0),
609 fHistoEleLambdaMassRSAway(0),
610 fHistoEleLambdaMassRS1Away(0),
611 fHistoEleLambdaMassRS2Away(0),
612 fHistoEleLambdaMassWSAway(0),
613 fHistoEleLambdaMassWS1Away(0),
614 fHistoEleLambdaMassWS2Away(0),
615 fHistoEleLambdaMassRSMixAway(0),
616 fHistoEleLambdaMassRSMix1Away(0),
617 fHistoEleLambdaMassRSMix2Away(0),
618 fHistoEleLambdaMassWSMixAway(0),
619 fHistoEleLambdaMassWSMix1Away(0),
620 fHistoEleLambdaMassWSMix2Away(0),
621 fHistoEleLambdaMassRSSideAway(0),
622 fHistoEleLambdaMassRSSide1Away(0),
623 fHistoEleLambdaMassRSSide2Away(0),
624 fHistoEleLambdaMassWSSideAway(0),
625 fHistoEleLambdaMassWSSide1Away(0),
626 fHistoEleLambdaMassWSSide2Away(0),
627 fHistoEleLambdaMassvsElePtRS(0),
628 fHistoEleLambdaMassvsElePtWS(0),
629 fHistoEleLambdaMassvsElePtRSMix(0),
630 fHistoEleLambdaMassvsElePtWSMix(0),
631 fHistoEleLambdaMassvsElePtRSSide(0),
632 fHistoEleLambdaMassvsElePtWSSide(0),
633 fHistoEleLambdaMassvsElePtRS1(0),
634 fHistoEleLambdaMassvsElePtWS1(0),
635 fHistoEleLambdaMassvsElePtRSMix1(0),
636 fHistoEleLambdaMassvsElePtWSMix1(0),
637 fHistoEleLambdaMassvsElePtRSSide1(0),
638 fHistoEleLambdaMassvsElePtWSSide1(0),
639 fHistoEleLambdaMassvsElePtRS2(0),
640 fHistoEleLambdaMassvsElePtWS2(0),
641 fHistoEleLambdaMassvsElePtRSMix2(0),
642 fHistoEleLambdaMassvsElePtWSMix2(0),
643 fHistoEleLambdaMassvsElePtRSSide2(0),
644 fHistoEleLambdaMassvsElePtWSSide2(0),
649 fHistoEleLambdaMassMCS(0),
650 fHistoEleLambdaMassLambdabMCS(0),
651 fHistoEleLambdaMassXibMCS(0),
652 fHistoEleLambdaMassMCGen(0),
653 fHistoEleLambdaMassvsElePtMCS(0),
654 fHistoEleLambdaMassvsElePtMCGen(0),
655 fHistoEleLambdaMassvsElePtMCS1(0),
656 fHistoEleLambdaMassvsElePtMCGen1(0),
657 fHistoEleLambdaMassvsElePtMCS2(0),
658 fHistoEleLambdaMassvsElePtMCGen2(0),
661 fHistoElePtvsEtaRS(0),
662 fHistoElePtvsEtaWS(0),
663 fHistoElePtvsEtaRSMix(0),
664 fHistoElePtvsEtaWSMix(0),
665 fHistoElePtvsEtaMCS(0),
666 fHistoElePtvsEtaMCGen(0),
667 fHistoElePtvsLambdaPtRS(0),
668 fHistoElePtvsLambdaPtWS(0),
669 fHistoElePtvsLambdaPtRSMix(0),
670 fHistoElePtvsLambdaPtWSMix(0),
671 fHistoElePtvsLambdaPtMCS(0),
672 fHistoElePtvsLambdaPtvsLcPtMCS(0),
673 fHistoElePtvsLambdaPtMCGen(0),
674 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
675 fHistoElePtvsLambdaPtMCLcGen(0),
676 fHistoElePtvsd0RS(0),
677 fHistoElePtvsd0WS(0),
678 fHistoElePtvsd0RSMix(0),
679 fHistoElePtvsd0WSMix(0),
680 fHistoElePtvsd0MCS(0),
681 fHistoElePtvsd0PromptMCS(0),
682 fHistoElePtvsd0BFeeddownMCS(0),
683 fHistoEleLambdaMassFeeddownXic0MCS(0),
684 fHistoEleLambdaMassFeeddownXic0MCGen(0),
685 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
686 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
687 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
688 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
689 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
690 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
691 fHistoElePtFeeddownXic0MCS(0),
692 fHistoElePtFeeddownXic0MCGen(0),
693 fHistoElePtvsEtaFeeddownXic0MCS(0),
694 fHistoElePtvsEtaFeeddownXic0MCGen(0),
695 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
696 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
697 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
698 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
699 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
700 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
701 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
702 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
703 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
704 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
705 fHistoElePtFeeddownXicPlusMCS(0),
706 fHistoElePtFeeddownXicPlusMCGen(0),
707 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
708 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
709 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
710 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
713 fHistoBachPtMCGen(0),
715 fHistoLambdaMassvsPt(0),
716 fHistoLambdaMassvsPtMCS(0),
717 fHistoLambdaMassvsPtMCGen(0),
718 fHistoLambdaPtvsEta(0),
719 fHistoLambdaPtvsEtaMCS(0),
720 fHistoLambdaPtvsEtaMCGen(0),
721 fHistoLambdaArmenteros(0),
722 fHistoLambdaArmenterosMCS(0),
723 fHistoLambdaPtvsDl(0),
724 fHistoLambdaPtvsDlSide(0),
725 fHistoLambdaPtvsDlMCS(0),
726 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
727 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
728 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
729 fHistoK0sMassvsPt(0),
730 fHistoElectronTPCPID(0),
731 fHistoElectronTOFPID(0),
732 fHistoElectronTPCSelPID(0),
733 fHistoElectronTOFSelPID(0),
734 fHistoElectronTPCPIDSelTOF(0),
735 fHistoElectronTPCPIDSelTOFSmallEta(0),
736 fHistoElectronTPCPIDSelTOFLargeEta(0),
737 fHistoMassConversionsMin(0),
738 fHistoMassConversionsSameSignMin(0),
739 fHistoElectronQovPtvsPhi(0),
740 fHistoLambdaQovPtvsPhi(0),
747 fHistoFeedDownXic0MCGen(0),
748 fHistoFeedDownXic0MCGen1(0),
749 fHistoFeedDownXic0MCGen2(0),
750 fHistoFeedDownXic0MCS(0),
751 fHistoFeedDownXic0MCS1(0),
752 fHistoFeedDownXic0MCS2(0),
753 fHistoFeedDownXicPlusMCGen(0),
754 fHistoFeedDownXicPlusMCGen1(0),
755 fHistoFeedDownXicPlusMCGen2(0),
756 fHistoFeedDownXicPlusMCS(0),
757 fHistoFeedDownXicPlusMCS1(0),
758 fHistoFeedDownXicPlusMCS2(0),
759 fHistoFeedDownSigmaMCGen(0),
760 fHistoFeedDownSigmaMCGen1(0),
761 fHistoFeedDownSigmaMCGen2(0),
762 fHistoFeedDownSigmaMCS(0),
763 fHistoFeedDownSigmaMCS1(0),
764 fHistoFeedDownSigmaMCS2(0),
765 fHistoLcElectronMCGen(0),
766 fHistoLcElectronMCGen1(0),
767 fHistoLcElectronMCGen2(0),
768 fHistoLcElectronMCS(0),
769 fHistoLcElectronMCS1(0),
770 fHistoLcElectronMCS2(0),
771 fHistoElectronFeedDownXic0MCGen(0),
772 fHistoElectronFeedDownXic0MCGen1(0),
773 fHistoElectronFeedDownXic0MCGen2(0),
774 fHistoElectronFeedDownXic0MCS(0),
775 fHistoElectronFeedDownXic0MCS1(0),
776 fHistoElectronFeedDownXic0MCS2(0),
777 fHistoElectronFeedDownXicPlusMCGen(0),
778 fHistoElectronFeedDownXicPlusMCGen1(0),
779 fHistoElectronFeedDownXicPlusMCGen2(0),
780 fHistoElectronFeedDownXicPlusMCS(0),
781 fHistoElectronFeedDownXicPlusMCS1(0),
782 fHistoElectronFeedDownXicPlusMCS2(0),
783 fHistoElectronMCGen(0),
784 fHistoLambdaMCGen(0),
785 fHistoElePtvsV0dlRS(0),
786 fHistoElePtvsV0dlRS1(0),
787 fHistoElePtvsV0dlRS2(0),
788 fHistoElePtvsV0dlRSSide(0),
789 fHistoElePtvsV0dlRSSide1(0),
790 fHistoElePtvsV0dlRSSide2(0),
791 fHistoElePtvsV0dlRSMix(0),
792 fHistoElePtvsV0dlRSMix1(0),
793 fHistoElePtvsV0dlRSMix2(0),
794 fHistoElePtvsV0dlWS(0),
795 fHistoElePtvsV0dlWS1(0),
796 fHistoElePtvsV0dlWS2(0),
797 fHistoElePtvsV0dlWSSide(0),
798 fHistoElePtvsV0dlWSSide1(0),
799 fHistoElePtvsV0dlWSSide2(0),
800 fHistoElePtvsV0dlWSMix(0),
801 fHistoElePtvsV0dlWSMix1(0),
802 fHistoElePtvsV0dlWSMix2(0),
803 fHistoElePtvsV0dlMCS(0),
804 fHistoElePtvsV0dlMCS1(0),
805 fHistoElePtvsV0dlMCS2(0),
806 fHistoElePtvsV0dlFeedDownXic0MCS(0),
807 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
808 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
809 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
810 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
811 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
812 fHistoElePtvsV0dcaRS(0),
813 fHistoElePtvsV0dcaRS1(0),
814 fHistoElePtvsV0dcaRS2(0),
815 fHistoElePtvsV0dcaRSSide(0),
816 fHistoElePtvsV0dcaRSSide1(0),
817 fHistoElePtvsV0dcaRSSide2(0),
818 fHistoElePtvsV0dcaRSMix(0),
819 fHistoElePtvsV0dcaRSMix1(0),
820 fHistoElePtvsV0dcaRSMix2(0),
821 fHistoElePtvsV0dcaWS(0),
822 fHistoElePtvsV0dcaWS1(0),
823 fHistoElePtvsV0dcaWS2(0),
824 fHistoElePtvsV0dcaWSSide(0),
825 fHistoElePtvsV0dcaWSSide1(0),
826 fHistoElePtvsV0dcaWSSide2(0),
827 fHistoElePtvsV0dcaWSMix(0),
828 fHistoElePtvsV0dcaWSMix1(0),
829 fHistoElePtvsV0dcaWSMix2(0),
830 fHistoElePtvsV0dcaMCS(0),
831 fHistoElePtvsV0dcaMCS1(0),
832 fHistoElePtvsV0dcaMCS2(0),
833 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
834 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
835 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
836 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
837 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
838 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
839 fHistoEleLambdaPtvsV0dlRS(0),
840 fHistoEleLambdaPtvsV0dlRS1(0),
841 fHistoEleLambdaPtvsV0dlRS2(0),
842 fHistoEleLambdaPtvsV0dlRSSide(0),
843 fHistoEleLambdaPtvsV0dlRSSide1(0),
844 fHistoEleLambdaPtvsV0dlRSSide2(0),
845 fHistoEleLambdaPtvsV0dlRSMix(0),
846 fHistoEleLambdaPtvsV0dlRSMix1(0),
847 fHistoEleLambdaPtvsV0dlRSMix2(0),
848 fHistoEleLambdaPtvsV0dlWS(0),
849 fHistoEleLambdaPtvsV0dlWS1(0),
850 fHistoEleLambdaPtvsV0dlWS2(0),
851 fHistoEleLambdaPtvsV0dlWSSide(0),
852 fHistoEleLambdaPtvsV0dlWSSide1(0),
853 fHistoEleLambdaPtvsV0dlWSSide2(0),
854 fHistoEleLambdaPtvsV0dlWSMix(0),
855 fHistoEleLambdaPtvsV0dlWSMix1(0),
856 fHistoEleLambdaPtvsV0dlWSMix2(0),
857 fHistoEleLambdaPtvsV0dlMCS(0),
858 fHistoEleLambdaPtvsV0dlMCS1(0),
859 fHistoEleLambdaPtvsV0dlMCS2(0),
860 fHistoEleLambdaPtvsV0dlPromptMCS(0),
861 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
862 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
863 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
864 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
865 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
866 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
867 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
868 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
869 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
870 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
871 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
872 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
873 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
874 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
875 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
876 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
877 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
878 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
879 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
880 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
881 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
882 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
883 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
884 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
885 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
886 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
887 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
888 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
889 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
890 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
891 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
892 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
893 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
894 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
895 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
896 fHistoEleLambdaPtvsV0dlXibMCS(0),
897 fHistoEleLambdaPtvsV0dlXibMCS1(0),
898 fHistoEleLambdaPtvsV0dlXibMCS2(0),
899 fHistoEleLambdaPtvsV0dlRSAway(0),
900 fHistoEleLambdaPtvsV0dlRS1Away(0),
901 fHistoEleLambdaPtvsV0dlRS2Away(0),
902 fHistoEleLambdaPtvsV0dlRSSideAway(0),
903 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
904 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
905 fHistoEleLambdaPtvsV0dlRSMixAway(0),
906 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
907 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
908 fHistoEleLambdaPtvsV0dlWSAway(0),
909 fHistoEleLambdaPtvsV0dlWS1Away(0),
910 fHistoEleLambdaPtvsV0dlWS2Away(0),
911 fHistoEleLambdaPtvsV0dlWSSideAway(0),
912 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
913 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
914 fHistoEleLambdaPtvsV0dlWSMixAway(0),
915 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
916 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
917 fHistoResponseElePt(0),
918 fHistoResponseElePt1(0),
919 fHistoResponseElePt2(0),
920 fHistoResponseEleLambdaPt(0),
921 fHistoResponseEleLambdaPt1(0),
922 fHistoResponseEleLambdaPt2(0),
923 fHistoResponseEleLambdaPtFeeddownXic0(0),
924 fHistoResponseEleLambdaPtFeeddownXic01(0),
925 fHistoResponseEleLambdaPtFeeddownXic02(0),
926 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
927 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
928 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
929 fHistoResponseEleLambdaPtFeeddownSigma(0),
930 fHistoResponseEleLambdaPtFeeddownSigma1(0),
931 fHistoResponseEleLambdaPtFeeddownSigma2(0),
932 fHistoResponseLambdaPt(0),
933 fHistoResponseLambdaPtFeeddownXic0(0),
934 fHistoResponseLambdaPtFeeddownXicPlus(0),
935 fHistoResponseLambdaPtFeeddownSigma(0),
936 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
937 fHistoEleLambdaPtvsRapidityRS(0),
938 fHistoEleLambdaPtvsRapidityWS(0),
939 fHistoEleLambdaPtvsRapidityMCS(0),
940 fHistoElectronPi0Total(0),
941 fHistoElectronPi0Tag(0),
942 fHistoElectronEtaTotal(0),
943 fHistoElectronEtaTag(0),
945 fHistonEvtvsRunNumber(0),
946 fHistonElevsRunNumber(0),
947 fHistonLambdavsRunNumber(0),
948 fHistoMCEventType(0),
949 fHistoMCDeltaPhiccbar(0),
950 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
951 fHistodPhiSdEtaSElectronProtonR125RS(0),
952 fHistodPhiSdEtaSElectronProtonR125WS(0),
953 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
954 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
955 fHistodPhiSdEtaSElectronPionR125RS(0),
956 fHistodPhiSdEtaSElectronPionR125WS(0),
957 fHistodPhiSdEtaSElectronPionR125RSMix(0),
958 fHistodPhiSdEtaSElectronPionR125WSMix(0),
960 fMixWithoutConversionFlag(kTRUE),
961 fNumberOfEventsForMixing (5),
967 fElectronTracks(0x0),
974 fElectronCutVarsArray(0x0),
975 fV0CutVarsArray1(0x0),
976 fV0CutVarsArray2(0x0)
981 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
983 for(Int_t i=0;i<17;i++){
988 for(Int_t i=0;i<8;i++){
992 DefineOutput(1,TList::Class());
993 DefineOutput(2,TList::Class());
994 DefineOutput(3,TList::Class());
995 DefineOutput(4,TTree::Class());
996 DefineOutput(5,TTree::Class());
997 DefineOutput(6,TTree::Class());
998 DefineOutput(7,TTree::Class());
999 DefineOutput(8,AliNormalizationCounter::Class());
1000 DefineOutput(9,TTree::Class());
1001 DefineOutput(10,TTree::Class());
1002 DefineOutput(11,TTree::Class());
1003 DefineOutput(12,TTree::Class());
1011 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1106 if (fDebug > 1) AliInfo(
"Init");
1125 AliError(
"NO EVENT FOUND!");
1128 AliAODEvent* aodEvent =
dynamic_cast<AliAODEvent*
>(fInputEvent);
1136 fBzkG = (Double_t)aodEvent->GetMagneticField();
1137 AliKFParticle::SetField(
fBzkG);
1138 if (TMath::Abs(
fBzkG)<0.001) {
1149 TClonesArray *mcArray = 0;
1150 AliAODMCHeader *mcHeader=0;
1153 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1155 AliError(
"Could not find Monte-Carlo in AOD");
1161 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1163 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1168 Double_t zMCVertex = mcHeader->GetVtxZ();
1170 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1184 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1187 Double_t pos[3],cov[6];
1189 fVtx1->GetCovarianceMatrix(cov);
1190 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1194 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1201 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1202 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1203 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1204 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1205 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1219 AliCentrality *cent = aodEvent->GetCentrality();
1220 fCentrality = cent->GetCentralityPercentile(
"V0M");
1224 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1231 Int_t runnumber_offset = 0;
1232 Int_t runnumber = aodEvent->GetRunNumber();
1233 if(runnumber<=131000&&runnumber>=114000){
1234 runnumber_offset = 114000;
1235 }
else if(runnumber<=196000&&runnumber>=195000){
1236 runnumber_offset = 195000;
1237 }
else if(runnumber<=170593&&runnumber>=167902){
1238 runnumber_offset = 167902;
1283 AliAnalysisTaskSE::Terminate();
1285 fOutput =
dynamic_cast<TList*
> (GetOutputData(1));
1287 AliError(
"fOutput not available");
1291 fOutputAll =
dynamic_cast<TList*
> (GetOutputData(3));
1293 AliError(
"fOutputAll not available");
1348 TString normName=
"NormalizationCounter";
1349 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1350 if(cont)normName=(TString)cont->GetName();
1372 fEventBuffer[i]=
new TTree(Form(
"EventBuffer_%d",i),
"Temporary buffer for event mixing");
1396 AliAODEvent *aodEvent, TClonesArray *mcArray
1403 if(fElectronTracks) fElectronTracks->Delete();
1404 if(fV0Tracks1) fV0Tracks1->Delete();
1405 if(fV0Tracks2) fV0Tracks2->Delete();
1406 if(fElectronCutVarsArray) fElectronCutVarsArray->Delete();
1407 if(fV0CutVarsArray1) fV0CutVarsArray1->Delete();
1408 if(fV0CutVarsArray2) fV0CutVarsArray2->Delete();
1409 fV0dlArray1.clear();
1410 fV0dlArray2.clear();
1411 fV0dcaArray1.clear();
1412 fV0dcaArray2.clear();
1415 ResetGlobalTrackReference();
1417 for (Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1419 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1420 if (!track)
continue;
1423 StoreGlobalTrackReference(track,iTrack);
1430 Int_t nV0s= aodEvent->GetNumberOfV0s();
1431 Int_t nTracks= aodEvent->GetNumberOfTracks();
1433 Bool_t seleTrkFlags[nTracks];
1435 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1437 Bool_t seleV0Flags[nV0s];
1439 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1441 Int_t runnumber_offset = 0;
1442 Int_t runnumber = aodEvent->GetRunNumber();
1443 if(runnumber<=131000&&runnumber>=114000){
1444 runnumber_offset = 114000;
1445 }
else if(runnumber<=196000&&runnumber>=195000){
1446 runnumber_offset = 195000;
1447 }
else if(runnumber<=170593&&runnumber>=167902){
1448 runnumber_offset = 167902;
1450 fHistonElevsRunNumber->Fill(runnumber-runnumber_offset,nSeleTrks);
1451 fHistonLambdavsRunNumber->Fill(runnumber-runnumber_offset,nSeleV0);
1456 for (Int_t iv0 = 0; iv0<nV0s; iv0++) {
1457 if(!seleV0Flags[iv0])
continue;
1458 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1461 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1462 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1467 for (Int_t itrk = 0; itrk<nTracks; itrk++) {
1468 if(!seleTrkFlags[itrk])
continue;
1469 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1477 AliAODTrack *trkpid = 0;
1478 if(fAnalCuts->GetProdAODFilterBit()==7){
1479 trkpid = fGTI[-trk->GetID()-1];
1484 Int_t cpid = cptrack->GetID();
1485 Int_t cnid = cntrack->GetID();
1486 Int_t lpid = trkpid->GetID();
1487 if((cpid==lpid)||(cnid==lpid))
continue;
1489 AliAODVertex *secVert = ReconstructSecondaryVertex(v0,trk,aodEvent);
1490 if(!secVert)
continue;
1497 FillROOTObjects(elobj,v0,trk,trkpid,aodEvent,mcArray);
1499 elobj->GetSecondaryVtx()->RemoveDaughters();
1501 delete elobj;elobj=NULL;
1507 fEventInfo->SetString(Form(
"Ev%d_esd%d_E%d_V%d",AliAnalysisManager::GetAnalysisManager()->GetNcalls(),((AliAODHeader*)aodEvent->GetHeader())->GetEventNumberESDFile(),fElectronTracks->GetEntries(),fV0Tracks1->GetEntries()+fV0Tracks2->GetEntries()));
1508 Int_t ind=GetPoolIndex(fVtxZ,fCentrality);
1509 if(ind>=0 && ind<fNOfPools){
1510 if(fEventBuffer[ind]->GetEntries() >= fNumberOfEventsForMixing){
1511 DoEventMixingWithPools(ind);
1512 if(fEventBuffer[ind]->GetEntries() >= 20*fNumberOfEventsForMixing){
1516 fEventBuffer[ind]->Fill();
1527 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1528 fVariablesTree =
new TTree(nameoutput,
"Candidates variables tree");
1531 TString * fCandidateVariableNames =
new TString[nVar];
1533 fCandidateVariableNames[ 0]=
"Centrality";
1534 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1535 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1536 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1537 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1538 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1539 fCandidateVariableNames[ 6]=
"ElePx";
1540 fCandidateVariableNames[ 7]=
"ElePy";
1541 fCandidateVariableNames[ 8]=
"ElePz";
1542 fCandidateVariableNames[ 9]=
"V0Px";
1543 fCandidateVariableNames[10]=
"V0Py";
1544 fCandidateVariableNames[11]=
"V0Pz";
1545 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1546 fCandidateVariableNames[13]=
"MassLambda";
1547 fCandidateVariableNames[14]=
"MassAntiLambda";
1548 fCandidateVariableNames[15]=
"Eled0";
1549 fCandidateVariableNames[16]=
"V0d0";
1550 fCandidateVariableNames[17]=
"nSigmaTPCele";
1551 fCandidateVariableNames[18]=
"nSigmaTOFele";
1552 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1553 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1554 fCandidateVariableNames[21]=
"EleCharge";
1555 fCandidateVariableNames[22]=
"ProtonPx";
1556 fCandidateVariableNames[23]=
"ProtonPy";
1557 fCandidateVariableNames[24]=
"ProtonPz";
1558 fCandidateVariableNames[25]=
"PiPx";
1559 fCandidateVariableNames[26]=
"PiPy";
1560 fCandidateVariableNames[27]=
"PiPz";
1561 fCandidateVariableNames[28]=
"mcpdglc";
1562 fCandidateVariableNames[29]=
"mclablc";
1563 fCandidateVariableNames[30]=
"mcpdgmomele";
1564 fCandidateVariableNames[31]=
"mcpdgmomv0";
1565 fCandidateVariableNames[32]=
"Mixing";
1566 fCandidateVariableNames[33]=
"mcpdgele";
1567 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1568 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1569 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1570 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1571 fCandidateVariableNames[38]=
"MassK0Short";
1572 fCandidateVariableNames[39]=
"mcpdggrmomele";
1573 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1574 fCandidateVariableNames[41]=
"mcngenele";
1575 fCandidateVariableNames[42]=
"mcngenv0";
1576 fCandidateVariableNames[43]=
"mclcpx";
1577 fCandidateVariableNames[44]=
"mclcpy";
1578 fCandidateVariableNames[45]=
"mclcpz";
1579 fCandidateVariableNames[46]=
"mcelepx";
1580 fCandidateVariableNames[47]=
"mcelepy";
1581 fCandidateVariableNames[48]=
"mcelepz";
1582 fCandidateVariableNames[49]=
"mcv0px";
1583 fCandidateVariableNames[50]=
"mcv0py";
1584 fCandidateVariableNames[51]=
"mcv0pz";
1585 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1586 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1587 fCandidateVariableNames[54]=
"PrimVertx";
1588 fCandidateVariableNames[55]=
"PrimVerty";
1589 fCandidateVariableNames[56]=
"PrimVertz";
1590 fCandidateVariableNames[57]=
"V0Vertx";
1591 fCandidateVariableNames[58]=
"V0Verty";
1592 fCandidateVariableNames[59]=
"V0Vertz";
1594 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1595 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1596 fCandidateVariableNames[62]=
"DcaV0daughters";
1597 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1598 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1599 fCandidateVariableNames[65]=
"MassK0Short2";
1601 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1602 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1604 fCandidateVariableNames[68]=
"EleITSMatch";
1605 fCandidateVariableNames[69]=
"V0PosITSMatch";
1606 fCandidateVariableNames[70]=
"V0NegITSMatch";
1607 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1608 fCandidateVariableNames[72]=
"mcpdgv0";
1610 fCandidateVariableNames[73]=
"mcv0pospx";
1611 fCandidateVariableNames[74]=
"mcv0pospy";
1612 fCandidateVariableNames[75]=
"mcv0pospz";
1613 fCandidateVariableNames[76]=
"mcpdgv0pos";
1614 fCandidateVariableNames[77]=
"mcv0negpx";
1615 fCandidateVariableNames[78]=
"mcv0negpy";
1616 fCandidateVariableNames[79]=
"mcv0negpz";
1617 fCandidateVariableNames[80]=
"mcpdgv0neg";
1619 fCandidateVariableNames[81]=
"ProtonPxPV";
1620 fCandidateVariableNames[82]=
"ProtonPyPV";
1621 fCandidateVariableNames[83]=
"ProtonPzPV";
1622 fCandidateVariableNames[84]=
"PiPxPV";
1623 fCandidateVariableNames[85]=
"PiPyPV";
1624 fCandidateVariableNames[86]=
"PiPzPV";
1626 fCandidateVariableNames[87]=
"mcv0vertx";
1627 fCandidateVariableNames[88]=
"mcv0verty";
1628 fCandidateVariableNames[89]=
"mcv0vertz";
1630 fCandidateVariableNames[90]=
"OnFlyStatus";
1631 fCandidateVariableNames[91]=
"MagneticField";
1632 fCandidateVariableNames[92]=
"EvNumber";
1633 fCandidateVariableNames[93]=
"RunNumber";
1636 for (Int_t ivar=0; ivar<nVar; ivar++) {
1653 for(Int_t i=0;i<94;i++)
1658 Bool_t anti_lambda_flag = kFALSE;
1661 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1662 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1668 Double_t v0px = elobj->PxProng(1);
1669 Double_t v0py = elobj->PyProng(1);
1670 Double_t v0pz = elobj->PzProng(1);
1671 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1672 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1673 Double_t epx = elobj->PxProng(0);
1674 Double_t epy = elobj->PyProng(0);
1675 Double_t epz = elobj->PzProng(0);
1676 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1677 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1678 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
1679 Double_t Esum = Ee + Ev0;
1681 Double_t uxe = epx/mome;
1682 Double_t uye = epy/mome;
1683 Double_t uze = epz/mome;
1684 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
1685 Double_t pxv_flip = v0px + lf * uxe;
1686 Double_t pyv_flip = v0py + lf * uye;
1687 Double_t pzv_flip = v0pz + lf * uze;
1688 Double_t pxsum_flip = epx + pxv_flip;
1689 Double_t pysum_flip = epy + pyv_flip;
1690 Double_t pzsum_flip = epz + pzv_flip;
1691 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
1692 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
1694 Double_t minmass_ee = 9999.;
1696 Double_t minmasslike_ee = 9999.;
1700 UInt_t pdgdg[2]={11,3122};
1718 Double_t nSigmaTPCele=-9999.;
1719 Double_t nSigmaTOFele=-9999.;
1720 Double_t nSigmaTPCv0pr=-9999.;
1721 Double_t nSigmaTOFv0pr=-9999.;
1722 Double_t nSigmaTPCv0pi=-9999.;
1723 Double_t nSigmaTOFv0pi=-9999.;
1734 if(anti_lambda_flag){
1752 if(anti_lambda_flag){
1768 AliAODMCParticle *mclc = 0;
1769 AliAODMCParticle *mcele = 0;
1770 AliAODMCParticle *mcv0 = 0;
1772 Int_t mcpdgele_array[100];
1773 Int_t mcpdgv0_array[100];
1774 Int_t mclabelele_array[100];
1775 Int_t mclabelv0_array[100];
1776 Int_t mcngen_ele=-9999;
1777 Int_t mcngen_v0=-9999;
1780 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
1781 if(mclabelele_array[0]>=0){
1782 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
1789 if(mclabelv0_array[0]>=0){
1790 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
1798 Int_t pdgdgv0[2]={211,211};
1799 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
1801 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
1811 Int_t labPos = cptrack->GetLabel();
1813 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
1822 Int_t labNeg = cntrack->GetLabel();
1824 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
1832 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
1875 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
1876 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
1877 if(!anti_lambda_flag){
1885 Double_t posVtx[3] = {0.,0.,0.};
1886 fVtx1->GetXYZ(posVtx);
1888 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
1898 if(!anti_lambda_flag){
1923 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
1924 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
1928 cont[0] = elobj->InvMass(2,pdgdg);
1929 cont[1] = elobj->Pt();
1933 Double_t cont_flip[3];
1934 cont_flip[0] = melam_flip;
1935 cont_flip[1] = ptelam_flip;
1939 cont2[0] = elobj->InvMass(2,pdgdg);
1940 cont2[1] = trk->Pt();
1942 Double_t cont_eleptvseta[3];
1943 cont_eleptvseta[0] = trk->Pt();
1944 cont_eleptvseta[1] = trk->Eta();
1947 Double_t cont_eleptvslambdapt[3];
1948 cont_eleptvslambdapt[0] = trk->Pt();
1949 cont_eleptvslambdapt[1] = v0->Pt();
1952 Double_t cont_eleptvsd0[3];
1953 cont_eleptvsd0[0] = trk->Pt();
1954 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
1957 Double_t cont_eleptvsv0dl[3];
1958 cont_eleptvsv0dl[0] = trk->Pt();
1959 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
1962 Double_t cont_elelamptvsv0dl[4];
1963 cont_elelamptvsv0dl[0] = elobj->Pt();
1964 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
1965 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
1968 Double_t cont_elelamptvsv0dl_flip[4];
1969 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
1970 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
1971 cont_elelamptvsv0dl_flip[2] = 0.0;
1975 Double_t cont_eleptvsv0dca[3];
1976 cont_eleptvsv0dca[0] = trk->Pt();
1977 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
1980 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
1981 Double_t elobj_px = elobj->Px();
1982 Double_t elobj_py = elobj->Py();
1983 Double_t elobj_pz = elobj->Pz();
1984 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
1985 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
1992 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2019 for(Int_t ih=0;ih<17;ih++){
2020 Double_t cont_eleptvscutvars[3];
2021 cont_eleptvscutvars[0] = elobj->Pt();
2025 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2027 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2029 cont_eleptvscutvars[1] = nSigmaTPCele;
2031 cont_eleptvscutvars[1] = nSigmaTOFele;
2033 cont_eleptvscutvars[1] = trk->Eta();
2035 cont_eleptvscutvars[1] = trk->GetITSNcls();
2037 if(!anti_lambda_flag)
2038 cont_eleptvscutvars[1] = v0->MassLambda();
2040 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2043 lPosV0[0] = v0->DecayVertexV0X();
2044 lPosV0[1] = v0->DecayVertexV0Y();
2045 lPosV0[2] = v0->DecayVertexV0Z();
2046 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2048 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2050 if(!anti_lambda_flag)
2051 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2053 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2055 if(!anti_lambda_flag)
2056 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2058 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2060 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2062 cont_eleptvscutvars[1] = v0->MassK0Short();
2064 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2066 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2068 cont_eleptvscutvars[1] = v0->Eta();
2070 Double_t v0px = elobj->PxProng(1);
2071 Double_t v0py = elobj->PyProng(1);
2072 Double_t v0pz = elobj->PzProng(1);
2073 Double_t epx = elobj->PxProng(0);
2074 Double_t epy = elobj->PyProng(0);
2075 Double_t epz = elobj->PzProng(0);
2076 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2078 cont_eleptvscutvars[1] = -9999.;
2086 Int_t pdgcode = mclc->GetPdgCode();
2087 Double_t cont_mclc[3];
2088 cont_mclc[0] = mclc->Pt();
2089 cont_mclc[1] = mclc->Y();
2091 Double_t cont_mcele[3];
2092 cont_mcele[0] = mcele->Pt();
2093 cont_mcele[1] = mcele->Eta();
2096 Int_t labmotherlc = mclc->GetMother();
2097 Bool_t isbottomfd = kFALSE;
2099 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2100 Int_t pdgmotherlc = motherlc->GetPdgCode();
2101 if(TMath::Abs(pdgmotherlc)==511||TMath::Abs(pdgmotherlc)==521||TMath::Abs(pdgmotherlc)==5122||TMath::Abs(pdgmotherlc)==5132||TMath::Abs(pdgmotherlc)==5232||TMath::Abs(pdgmotherlc)==5332){
2105 cont2[1] = mcele->Pt();
2107 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2158 Double_t cont_eleptvslambdaptvslcpt[4];
2159 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2160 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2161 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2162 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2165 Double_t cont_allpt[4];
2166 cont_allpt[0] = mclc->Pt();
2167 cont_allpt[1] = elobj->Pt();
2168 cont_allpt[2] = trk->Pt();
2169 cont_allpt[3] = v0->Pt();
2178 for(Int_t ih=0;ih<17;ih++){
2179 Double_t cont_eleptvscutvars[3];
2180 cont_eleptvscutvars[0] = elobj->Pt();
2184 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2186 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2188 cont_eleptvscutvars[1] = nSigmaTPCele;
2190 cont_eleptvscutvars[1] = nSigmaTOFele;
2192 cont_eleptvscutvars[1] = trk->Eta();
2194 cont_eleptvscutvars[1] = trk->GetITSNcls();
2196 if(!anti_lambda_flag)
2197 cont_eleptvscutvars[1] = v0->MassLambda();
2199 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2202 lPosV0[0] = v0->DecayVertexV0X();
2203 lPosV0[1] = v0->DecayVertexV0Y();
2204 lPosV0[2] = v0->DecayVertexV0Z();
2205 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2207 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2209 if(!anti_lambda_flag)
2210 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2212 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2214 if(!anti_lambda_flag)
2215 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2217 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2219 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2221 cont_eleptvscutvars[1] = v0->MassK0Short();
2223 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2225 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2227 cont_eleptvscutvars[1] = v0->Eta();
2229 Double_t v0px = elobj->PxProng(1);
2230 Double_t v0py = elobj->PyProng(1);
2231 Double_t v0pz = elobj->PzProng(1);
2232 Double_t epx = elobj->PxProng(0);
2233 Double_t epy = elobj->PyProng(0);
2234 Double_t epz = elobj->PzProng(0);
2235 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2237 cont_eleptvscutvars[1] = -9999.;
2244 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2291 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2338 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2368 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2395 for(Int_t ih=0;ih<17;ih++){
2396 Double_t cont_eleptvscutvars[3];
2397 cont_eleptvscutvars[0] = elobj->Pt();
2401 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2403 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2405 cont_eleptvscutvars[1] = nSigmaTPCele;
2407 cont_eleptvscutvars[1] = nSigmaTOFele;
2409 cont_eleptvscutvars[1] = trk->Eta();
2411 cont_eleptvscutvars[1] = trk->GetITSNcls();
2413 if(!anti_lambda_flag)
2414 cont_eleptvscutvars[1] = v0->MassLambda();
2416 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2419 lPosV0[0] = v0->DecayVertexV0X();
2420 lPosV0[1] = v0->DecayVertexV0Y();
2421 lPosV0[2] = v0->DecayVertexV0Z();
2422 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2424 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2426 if(!anti_lambda_flag)
2427 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2429 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2431 if(!anti_lambda_flag)
2432 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2434 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2436 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2438 cont_eleptvscutvars[1] = v0->MassK0Short();
2440 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2442 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2444 cont_eleptvscutvars[1] = v0->Eta();
2446 Double_t v0px = elobj->PxProng(1);
2447 Double_t v0py = elobj->PyProng(1);
2448 Double_t v0pz = elobj->PzProng(1);
2449 Double_t epx = elobj->PxProng(0);
2450 Double_t epy = elobj->PyProng(0);
2451 Double_t epz = elobj->PzProng(0);
2452 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2454 cont_eleptvscutvars[1] = -9999.;
2462 Int_t pdgcode = mclc->GetPdgCode();
2463 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2470 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2476 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2482 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2488 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2502 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2504 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2508 if(cont_flip[0]<2.3){
2517 if(cont_flip[0]<2.3){
2527 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2548 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2572 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2574 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2578 if(cont_flip[0]<2.3){
2587 if(cont_flip[0]<2.3){
2598 for(Int_t iv=0;iv<11;iv++){
2610 if(trk->Charge()>0){
2613 }
else if(trk->Charge()<0){
2622 Int_t pdgcode = mclc->GetPdgCode();
2623 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2626 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2629 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2632 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2635 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2638 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2641 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2644 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2647 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2683 for(Int_t i=0;i<94;i++)
2689 Double_t pxe = trke->Px();
2690 Double_t pye = trke->Py();
2691 Double_t pze = trke->Pz();
2692 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
2693 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
2695 Double_t pxv = v0->Px();
2696 Double_t pyv = v0->Py();
2697 Double_t pzv = v0->Pz();
2698 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
2699 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
2701 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
2703 Double_t pxsum = pxe + pxv;
2704 Double_t pysum = pye + pyv;
2705 Double_t pzsum = pze + pzv;
2706 Double_t Esum = Ee + Ev;
2707 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
2709 Double_t uxe = pxe/mome;
2710 Double_t uye = pye/mome;
2711 Double_t uze = pze/mome;
2712 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
2713 Double_t pxv_flip = pxv + lf * uxe;
2714 Double_t pyv_flip = pyv + lf * uye;
2715 Double_t pzv_flip = pzv + lf * uze;
2716 Double_t pxsum_flip = pxe + pxv_flip;
2717 Double_t pysum_flip = pye + pyv_flip;
2718 Double_t pzsum_flip = pze + pzv_flip;
2719 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
2720 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
2723 UInt_t pdgdg[2]={11,3122};
2776 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
2780 Double_t cont_flip[3];
2781 cont_flip[0] = mel_flip;
2782 cont_flip[1] = ptel_flip;
2787 cont2[1] = sqrt(pxe*pxe+pye*pye);
2789 Double_t cont_eleptvseta[3];
2790 cont_eleptvseta[0] = trke->Pt();
2791 cont_eleptvseta[1] = trke->Eta();
2794 Double_t cont_eleptvslambdapt[3];
2795 cont_eleptvslambdapt[0] = trke->Pt();
2796 cont_eleptvslambdapt[1] = v0->Pt();
2799 Double_t cont_eleptvsd0[3];
2800 cont_eleptvsd0[0] = trke->Pt();
2801 cont_eleptvsd0[1] = 0.;
2804 Double_t cont_eleptvsv0dl[3];
2805 cont_eleptvsv0dl[0] = trke->Pt();
2806 cont_eleptvsv0dl[1] = v0info[0];
2809 Double_t cont_eleptvsv0dca[3];
2810 cont_eleptvsv0dca[0] = trke->Pt();
2811 cont_eleptvsv0dca[1] = v0info[1];
2814 Double_t cont_elelamptvsv0dl[4];
2815 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
2816 cont_elelamptvsv0dl[1] = v0info[0];
2817 cont_elelamptvsv0dl[2] = 0.0;
2820 Double_t cont_elelamptvsv0dl_flip[4];
2821 cont_elelamptvsv0dl_flip[0] = ptel_flip;
2822 cont_elelamptvsv0dl_flip[1] = v0info[0];
2823 cont_elelamptvsv0dl_flip[2] = 0.0;
2826 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
2827 xyzR125_ele[0] = (*elevars)[0];
2828 xyzR125_ele[1] = (*elevars)[1];
2829 xyzR125_ele[2] = (*elevars)[2];
2830 xyzR125_pr[0] = (*v0vars)[0];
2831 xyzR125_pr[1] = (*v0vars)[1];
2832 xyzR125_pr[2] = (*v0vars)[2];
2833 xyzR125_pi[0] = (*v0vars)[3];
2834 xyzR125_pi[1] = (*v0vars)[4];
2835 xyzR125_pi[2] = (*v0vars)[5];
2837 Double_t rdhfcutvars[9];
2838 rdhfcutvars[0] = xyzR125_ele[0];
2839 rdhfcutvars[1] = xyzR125_ele[1];
2840 rdhfcutvars[2] = xyzR125_ele[2];
2841 rdhfcutvars[3] = xyzR125_pr[0];
2842 rdhfcutvars[4] = xyzR125_pr[1];
2843 rdhfcutvars[5] = xyzR125_pr[2];
2844 rdhfcutvars[6] = xyzR125_pi[0];
2845 rdhfcutvars[7] = xyzR125_pi[1];
2846 rdhfcutvars[8] = xyzR125_pi[2];
2857 if(((
int)trke->T())*chargepr>0){
2917 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
2919 if(((
int)trke->T())*chargepr>0){
2923 if(cont_flip[0]<2.3){
2932 if(cont_flip[0]<2.3){
2943 for(Int_t iv=0;iv<11;iv++){
2951 fCorrelationVariables[4] = v0info[0];
2952 fCorrelationVariables[5] = (*elevars)[5];
2953 fCorrelationVariables[6] = 2;
2955 if(chargepr>0) fCorrelationVariables[7] = 0;
2956 else fCorrelationVariables[7] = 2;
2957 }
else if(trke->T()<0){
2958 if(chargepr>0) fCorrelationVariables[7] = 1;
2959 else fCorrelationVariables[7] = 3;
2977 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
2981 TString * fCandidateVariableNames =
new TString[nVar];
2983 fCandidateVariableNames[ 0]=
"ElePx";
2984 fCandidateVariableNames[ 1]=
"ElePy";
2985 fCandidateVariableNames[ 2]=
"ElePz";
2986 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
2987 fCandidateVariableNames[ 4]=
"ITSNcls";
2988 fCandidateVariableNames[ 5]=
"TPCNcls";
2989 fCandidateVariableNames[ 6]=
"TPCNclsPID";
2990 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
2991 fCandidateVariableNames[ 8]=
"d0R";
2992 fCandidateVariableNames[ 9]=
"d0Z";
2993 fCandidateVariableNames[10]=
"ITSClusterMap";
2994 fCandidateVariableNames[11]=
"nSigmaTPCele";
2995 fCandidateVariableNames[12]=
"nSigmaTOFele";
2996 fCandidateVariableNames[13]=
"nSigmaTPCpi";
2997 fCandidateVariableNames[14]=
"nSigmaTPCka";
2998 fCandidateVariableNames[15]=
"nSigmaTPCpr";
2999 fCandidateVariableNames[16]=
"EvNumber";
3000 fCandidateVariableNames[17]=
"EleCharge";
3001 fCandidateVariableNames[18]=
"ElePdgCode";
3002 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3003 fCandidateVariableNames[20]=
"mcelepx";
3004 fCandidateVariableNames[21]=
"mcelepy";
3005 fCandidateVariableNames[22]=
"mcelepz";
3006 fCandidateVariableNames[23]=
"Centrality";
3007 fCandidateVariableNames[24]=
"PrimVertZ";
3008 fCandidateVariableNames[25]=
"RunNumber";
3010 for (Int_t ivar=0; ivar<nVar; ivar++) {
3028 Double_t d0z0[2],covd0z0[3];
3029 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3033 pv[0] =
fVtx1->GetX();
3034 pv[1] =
fVtx1->GetY();
3035 pv[2] =
fVtx1->GetZ();
3036 Double_t xyzR125[3]={9999.,9999.,9999.};
3038 TVector *varvec =
new TVector(6);
3039 (*varvec)[0] = xyzR125[0];
3040 (*varvec)[1] = xyzR125[1];
3041 (*varvec)[2] = xyzR125[2];
3044 (*varvec)[5] = d0z0[0];
3046 fElectronTracks->AddLast(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3050 Int_t pdgEle = -9999;
3051 Int_t pdgEleMother = -9999;
3052 Float_t mcelepx = -9999;
3053 Float_t mcelepy = -9999;
3054 Float_t mcelepz = -9999;
3057 Int_t labEle = trk->GetLabel();
3058 if(labEle<0)
return;
3059 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3061 pdgEle = mcetrk->GetPdgCode();
3062 mcelepx = mcetrk->Px();
3063 mcelepy = mcetrk->Py();
3064 mcelepz = mcetrk->Pz();
3081 for(Int_t i=0;i<26;i++){
3092 if(trk->GetTPCNclsF()>0)
3097 Int_t itsmap = trk->GetITSClusterMap();
3100 Bool_t spdfirst = (itsmap & bit1) == bit1;
3101 Bool_t spdsecond = (itsmap & bit2) == bit2;
3139 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3143 TString * fCandidateVariableNames =
new TString[nVar];
3145 fCandidateVariableNames[ 0]=
"V0Px";
3146 fCandidateVariableNames[ 1]=
"V0Py";
3147 fCandidateVariableNames[ 2]=
"V0Pz";
3148 fCandidateVariableNames[ 3]=
"MassLambda";
3149 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3150 fCandidateVariableNames[ 5]=
"ProtonPx";
3151 fCandidateVariableNames[ 6]=
"ProtonPy";
3152 fCandidateVariableNames[ 7]=
"ProtonPz";
3153 fCandidateVariableNames[ 8]=
"PionPx";
3154 fCandidateVariableNames[ 9]=
"PionPy";
3155 fCandidateVariableNames[10]=
"PionPz";
3156 fCandidateVariableNames[11]=
"RfidV0";
3157 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3158 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3159 fCandidateVariableNames[14]=
"DcaV0daughters";
3160 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3161 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3162 fCandidateVariableNames[17]=
"MassK0Short";
3163 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3164 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3165 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3166 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3167 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3168 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3169 fCandidateVariableNames[24]=
"V0PdgCode";
3170 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3171 fCandidateVariableNames[26]=
"mcv0px";
3172 fCandidateVariableNames[27]=
"mcv0py";
3173 fCandidateVariableNames[28]=
"mcv0pz";
3174 fCandidateVariableNames[29]=
"mcv0vertx";
3175 fCandidateVariableNames[30]=
"mcv0verty";
3176 fCandidateVariableNames[31]=
"mcv0vertz";
3177 fCandidateVariableNames[32]=
"EvNumber";
3178 fCandidateVariableNames[33]=
"Centrality";
3179 fCandidateVariableNames[34]=
"PrimVertZ";
3180 fCandidateVariableNames[35]=
"RunNumber";
3182 for (Int_t ivar=0; ivar<nVar; ivar++) {
3196 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3197 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3198 if(!cptrack)
return;
3199 if(!cntrack)
return;
3218 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3219 Double_t posVtx[3] = {0.,0.,0.};
3220 fVtx1->GetXYZ(posVtx);
3221 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3222 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3223 Double_t cont_arm[3];
3224 cont_arm[0] = v0->Pt();
3225 cont_arm[1] = v0->PtArmV0();
3226 cont_arm[2] = v0->AlphaV0();
3245 Double_t momv0x = v0->MomV0X();
3246 Double_t momv0y = v0->MomV0Y();
3247 Double_t phi_alice = atan2(momv0y,momv0x);
3248 if(phi_alice<0.) phi_alice += 2 * M_PI;
3251 Int_t v0pdgcode = -9999;
3252 Int_t v0motherpdgcode = -9999;
3253 Float_t mcv0px = -9999;
3254 Float_t mcv0py = -9999;
3255 Float_t mcv0pz = -9999;
3256 Float_t mcv0vertx = -9999;
3257 Float_t mcv0verty = -9999;
3258 Float_t mcv0vertz = -9999;
3261 Int_t pdgdgv0[2]={2212,211};
3262 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
3276 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
3277 if(!mcv0trk)
return;
3279 Bool_t hfv0_flag = kFALSE;
3280 Bool_t promptv0_flag = kFALSE;
3281 v0pdgcode = mcv0trk->GetPdgCode();
3282 Int_t labv0mother = mcv0trk->GetMother();
3284 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3286 v0motherpdgcode = motherv0->GetPdgCode();
3287 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
3290 if(abs(v0motherpdgcode)==3322){
3292 }
else if(abs(v0motherpdgcode)==3312){
3294 }
else if(abs(v0motherpdgcode)==3334){
3298 promptv0_flag = kTRUE;
3302 if(!hfv0_flag)
return;
3305 mcv0px = mcv0trk->Px();
3306 mcv0py = mcv0trk->Py();
3307 mcv0pz = mcv0trk->Pz();
3309 Int_t labptrk = cptrack->GetLabel();
3311 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
3313 mcv0vertx = mcv0ptrk->Xv();
3314 mcv0verty = mcv0ptrk->Yv();
3315 mcv0vertz = mcv0ptrk->Zv();
3322 Double_t posVtx[3] = {0.,0.,0.};
3323 fVtx1->GetXYZ(posVtx);
3324 TLorentzVector *lv =
new TLorentzVector();
3325 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3326 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
3327 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
3329 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
3331 fV0dlArray1.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3335 TVector *varvec =
new TVector(15);
3336 (*varvec)[0] = xyzR125pr[0];
3337 (*varvec)[1] = xyzR125pr[1];
3338 (*varvec)[2] = xyzR125pr[2];
3339 (*varvec)[3] = xyzR125pi[0];
3340 (*varvec)[4] = xyzR125pi[1];
3341 (*varvec)[5] = xyzR125pi[2];
3342 (*varvec)[6] = v0->DecayVertexV0X();
3343 (*varvec)[7] = v0->DecayVertexV0Y();
3344 (*varvec)[8] = v0->DecayVertexV0Z();
3345 (*varvec)[9] = v0->MomPosX();
3346 (*varvec)[10] = v0->MomPosY();
3347 (*varvec)[11] = v0->MomPosZ();
3348 (*varvec)[12] = v0->MomNegX();
3349 (*varvec)[13] = v0->MomNegY();
3350 (*varvec)[14] = v0->MomNegZ();
3353 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
3355 fV0dlArray2.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3359 TVector *varvec =
new TVector(15);
3360 (*varvec)[0] = xyzR125pr[0];
3361 (*varvec)[1] = xyzR125pr[1];
3362 (*varvec)[2] = xyzR125pr[2];
3363 (*varvec)[3] = xyzR125pi[0];
3364 (*varvec)[4] = xyzR125pi[1];
3365 (*varvec)[5] = xyzR125pi[2];
3366 (*varvec)[6] = v0->DecayVertexV0X();
3367 (*varvec)[7] = v0->DecayVertexV0Y();
3368 (*varvec)[8] = v0->DecayVertexV0Z();
3369 (*varvec)[9] = v0->MomNegX();
3370 (*varvec)[10] = v0->MomNegY();
3371 (*varvec)[11] = v0->MomNegZ();
3372 (*varvec)[12] = v0->MomPosX();
3373 (*varvec)[13] = v0->MomPosY();
3374 (*varvec)[14] = v0->MomPosZ();
3381 for(Int_t i=0;i<36;i++){
3391 Bool_t isparticle = kTRUE;
3411 lPosV0[0] = v0->DecayVertexV0X();
3412 lPosV0[1] = v0->DecayVertexV0Y();
3413 lPosV0[2] = v0->DecayVertexV0Z();
3414 Double_t decayvertV0 = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
3417 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
3418 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
3448 if(cptrack->GetTPCNclsF()>0)
3449 fCandidateV0Variables[22] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3450 if(cntrack->GetTPCNclsF()>0)
3451 fCandidateV0Variables[23] =(Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3455 if(cntrack->GetTPCNclsF()>0)
3456 fCandidateV0Variables[22] = (Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3457 if(cptrack->GetTPCNclsF()>0)
3458 fCandidateV0Variables[23] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3483 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
3487 TString * fCandidateVariableNames =
new TString[nVar];
3489 fCandidateVariableNames[ 0]=
"Centrality";
3490 fCandidateVariableNames[ 1]=
"DecayType";
3491 fCandidateVariableNames[ 2]=
"LcPx";
3492 fCandidateVariableNames[ 3]=
"LcPy";
3493 fCandidateVariableNames[ 4]=
"LcPz";
3494 fCandidateVariableNames[ 5]=
"ElePx";
3495 fCandidateVariableNames[ 6]=
"ElePy";
3496 fCandidateVariableNames[ 7]=
"ElePz";
3497 fCandidateVariableNames[ 8]=
"V0Px";
3498 fCandidateVariableNames[ 9]=
"V0Py";
3499 fCandidateVariableNames[10]=
"V0Pz";
3501 for (Int_t ivar=0; ivar<nVar; ivar++) {
3513 if(!mcepart)
return;
3514 if(!mcv0part)
return;
3516 for(Int_t i=0;i<11;i++){
3532 Double_t epx = mcepart->Px();
3533 Double_t epy = mcepart->Py();
3534 Double_t epz = mcepart->Pz();
3535 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
3536 Double_t v0px = mcv0part->Px();
3537 Double_t v0py = mcv0part->Py();
3538 Double_t v0pz = mcv0part->Pz();
3539 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
3541 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
3544 cont[0] = InvMassEleLambda;
3545 cont[1] = mcpart->Pt();
3548 cont2[0] = InvMassEleLambda;
3549 cont2[1] = mcepart->Pt();
3551 Double_t cont_eleptvseta[3];
3552 cont_eleptvseta[0] = mcepart->Pt();
3553 cont_eleptvseta[1] = mcepart->Eta();
3555 Double_t cont_eleptvslambdapt[3];
3556 cont_eleptvslambdapt[0] = mcepart->Pt();
3557 cont_eleptvslambdapt[1] = mcv0part->Pt();
3559 Double_t cont_eleptvslambdaptvslcpt[4];
3560 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
3561 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
3562 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
3565 contmc[0] = mcpart->Pt();
3566 contmc[1] = mcpart->Y();
3568 Double_t contmcele[3];
3569 contmcele[0] = mcepart->Pt();
3570 contmcele[1] = mcepart->Eta();
3575 esdcuts->GetEtaRange(etamin,etamax);
3585 if(fabs(mcepart->Eta())<etamax){
3589 if(InvMassEleLambda<2.3){
3595 if(fabs(mcpart->Y())<0.7){
3596 if(InvMassEleLambda<2.3){
3601 }
else if(decaytype==1){
3609 if(fabs(mcepart->Eta())<etamax){
3613 if(InvMassEleLambda<2.3){
3619 }
else if(decaytype==2){
3627 if(fabs(mcepart->Eta())<etamax){
3631 if(InvMassEleLambda<2.3){
3637 }
else if(decaytype==3){
3653 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
3657 TString * fCandidateVariableNames =
new TString[nVar];
3659 fCandidateVariableNames[ 0]=
"Centrality";
3660 fCandidateVariableNames[ 1]=
"ElePx";
3661 fCandidateVariableNames[ 2]=
"ElePy";
3662 fCandidateVariableNames[ 3]=
"ElePz";
3663 fCandidateVariableNames[ 4]=
"ElePdgCode";
3664 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
3665 fCandidateVariableNames[ 6]=
"RunNumber";
3666 fCandidateVariableNames[ 7]=
"EvNumber";
3668 for (Int_t ivar=0; ivar<nVar; ivar++) {
3679 if(!mcepart)
return;
3681 Bool_t hfe_flag = kFALSE;
3682 Int_t labemother = mcepart->GetMother();
3683 Int_t pdgmotherele = -9999;
3685 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
3686 pdgmotherele = motherele->GetPdgCode();
3687 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
3691 if(!hfe_flag)
return;
3694 contmc[0] = mcepart->Pt();
3695 contmc[1] = mcepart->Eta();
3699 for(Int_t i=0;i<8;i++){
3723 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
3727 TString * fCandidateVariableNames =
new TString[nVar];
3729 fCandidateVariableNames[ 0]=
"Centrality";
3730 fCandidateVariableNames[ 1]=
"V0Px";
3731 fCandidateVariableNames[ 2]=
"V0Py";
3732 fCandidateVariableNames[ 3]=
"V0Pz";
3733 fCandidateVariableNames[ 4]=
"V0PdgCode";
3734 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
3735 fCandidateVariableNames[ 6]=
"RunNumber";
3736 fCandidateVariableNames[ 7]=
"EvNumber";
3738 for (Int_t ivar=0; ivar<nVar; ivar++) {
3749 if(!mcv0part)
return;
3751 for(Int_t i=0;i<8;i++){
3755 Bool_t hfv0_flag = kFALSE;
3756 Int_t labv0mother = mcv0part->GetMother();
3757 Int_t pdgmotherv0 = -9999;
3759 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3761 pdgmotherv0 = motherv0->GetPdgCode();
3762 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
3767 if(!hfv0_flag)
return;
3770 contmc[0] = mcv0part->Pt();
3771 contmc[1] = mcv0part->Eta();
3795 fCEvents =
new TH1F(
"fCEvents",
"conter",18,-0.5,17.5);
3797 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
3798 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
3799 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
3800 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
3801 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
3802 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
3803 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
3804 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
3805 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
3806 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
3807 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
3808 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
3809 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
3810 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
3812 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
3813 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
3816 fCEvents->GetYaxis()->SetTitle(
"counts");
3818 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
3820 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
3821 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
3822 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
3823 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
3824 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
3825 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
3827 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
3828 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
3829 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
3830 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
3831 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
3833 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
3852 Int_t bins_base[3]= {23,100 ,10};
3853 Double_t xmin_base[3]={1.1,0 ,0.00};
3854 Double_t xmax_base[3]={5.7,20. ,100};
3855 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
3931 Int_t bins_base_elept[3]= {10,100 ,10};
3932 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
3933 Double_t xmax_base_elept[3]={3.1,10. ,100};
3992 fHistoElePtRS =
new TH2F(
"fHistoElePtRS",
"",100,0,10,10,0,100);
3994 fHistoElePtWS =
new TH2F(
"fHistoElePtWS",
"",100,0.,10.,10,0,100);
4000 fHistoElePtMCS =
new TH2F(
"fHistoElePtMCS",
"",100,0,10,10,0,100);
4005 Int_t bins_eleptvseta[3]= {10,10 ,10};
4006 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4007 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4009 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4011 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4013 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4015 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4017 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4019 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4022 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4023 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4024 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4026 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4028 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4030 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4032 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4034 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4036 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4038 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4041 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4042 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4043 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4044 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4046 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4049 Int_t bins_allpt[4]= {10,10,20,20};
4050 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4051 Double_t xmax_allpt[4]={20.,20.,10.,10};
4062 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4063 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4064 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4066 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4068 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4070 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4072 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4074 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4076 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4142 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4148 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4150 fHistoBachPtMCS =
new TH1F(
"fHistoBachPtMCS",
"Bachelor p_{T}",100,0.0,5.0);
4154 fHistod0Bach =
new TH1F(
"fHistod0Bach",
"Bachelor d_{0}",100,-0.5,0.5);
4156 fHistoLambdaMassvsPt=
new TH2F(
"fHistoLambdaMassvsPt",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4158 fHistoLambdaMassvsPtMCS=
new TH2F(
"fHistoLambdaMassvsPtMCS",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4168 fHistoK0sMassvsPt=
new TH2F(
"fHistoK0sMassvsPt",
"K0s mass",100,0.497-0.05,0.497+0.05,20,0.,10.);
4170 fHistoLambdaPtvsDl=
new TH2F(
"fHistoLambdaPtvsDl",
"Lambda pt vs dl",20,0.,10.,20,0.,40.);
4183 Int_t bins_armenteros[3]= {20 ,100 ,200};
4184 Double_t xmin_armenteros[3]={0.,0. ,-1.};
4185 Double_t xmax_armenteros[3]={20.,1. ,1.};
4186 fHistoLambdaArmenteros =
new THnSparseF(
"fHistoLambdaArmenteros",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4188 fHistoLambdaArmenterosMCS =
new THnSparseF(
"fHistoLambdaArmenterosMCS",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4210 for(Int_t i=0;i<8;i++){
4219 Int_t bins_lcmcgen[3]= {100 ,20 ,10};
4220 Double_t xmin_lcmcgen[3]={0.,-1.0 ,0.0};
4221 Double_t xmax_lcmcgen[3]={20.,1.0 ,100};
4222 fHistoLcMCGen =
new THnSparseF(
"fHistoLcMCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4224 fHistoLcMCGen1 =
new THnSparseF(
"fHistoLcMCGen1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4226 fHistoLcMCGen2 =
new THnSparseF(
"fHistoLcMCGen2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4228 fHistoFeedDownXic0MCGen =
new THnSparseF(
"fHistoFeedDownXic0MCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4247 fHistoLcMCS =
new THnSparseF(
"fHistoLcMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4249 fHistoLcMCS1 =
new THnSparseF(
"fHistoLcMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4251 fHistoLcMCS2 =
new THnSparseF(
"fHistoLcMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4253 fHistoFeedDownXic0MCS =
new THnSparseF(
"fHistoFeedDownXic0MCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4255 fHistoFeedDownXic0MCS1 =
new THnSparseF(
"fHistoFeedDownXic0MCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4257 fHistoFeedDownXic0MCS2 =
new THnSparseF(
"fHistoFeedDownXic0MCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4265 fHistoFeedDownSigmaMCS =
new THnSparseF(
"fHistoFeedDownSigmaMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4267 fHistoFeedDownSigmaMCS1 =
new THnSparseF(
"fHistoFeedDownSigmaMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4269 fHistoFeedDownSigmaMCS2 =
new THnSparseF(
"fHistoFeedDownSigmaMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4272 Int_t bins_elemcgen[3]= {100 ,20 ,10};
4273 Double_t xmin_elemcgen[3]={0.,-1.0 ,0.0};
4274 Double_t xmax_elemcgen[3]={10.,1.0 ,100};
4275 fHistoElectronMCGen =
new THnSparseF(
"fHistoElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4277 fHistoLcElectronMCGen =
new THnSparseF(
"fHistoLcElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4279 fHistoLcElectronMCGen1 =
new THnSparseF(
"fHistoLcElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4281 fHistoLcElectronMCGen2 =
new THnSparseF(
"fHistoLcElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4295 fHistoLcElectronMCS =
new THnSparseF(
"fHistoLcElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4297 fHistoLcElectronMCS1 =
new THnSparseF(
"fHistoLcElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4299 fHistoLcElectronMCS2 =
new THnSparseF(
"fHistoLcElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4314 Int_t bins_lambdamcgen[3]= {50 ,20 ,10};
4315 Double_t xmin_lambdamcgen[3]={0.,-1.0 ,0.0};
4316 Double_t xmax_lambdamcgen[3]={10.,1.0 ,100};
4317 fHistoLambdaMCGen =
new THnSparseF(
"fHistoLambdaMCGen",
"",3,bins_lambdamcgen,xmin_lambdamcgen,xmax_lambdamcgen);
4320 Int_t bins_eleptvsv0dl[3]= {100 ,20 ,10};
4321 Double_t xmin_eleptvsv0dl[3]={0.,0. ,0.0};
4322 Double_t xmax_eleptvsv0dl[3]={10.,40. ,100};
4323 fHistoElePtvsV0dlRS =
new THnSparseF(
"fHistoElePtvsV0dlRS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4325 fHistoElePtvsV0dlRS1 =
new THnSparseF(
"fHistoElePtvsV0dlRS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4327 fHistoElePtvsV0dlRS2 =
new THnSparseF(
"fHistoElePtvsV0dlRS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4329 fHistoElePtvsV0dlRSSide =
new THnSparseF(
"fHistoElePtvsV0dlRSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4331 fHistoElePtvsV0dlRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4333 fHistoElePtvsV0dlRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4335 fHistoElePtvsV0dlRSMix =
new THnSparseF(
"fHistoElePtvsV0dlRSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4337 fHistoElePtvsV0dlRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4339 fHistoElePtvsV0dlRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4341 fHistoElePtvsV0dlWS =
new THnSparseF(
"fHistoElePtvsV0dlWS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4343 fHistoElePtvsV0dlWS1 =
new THnSparseF(
"fHistoElePtvsV0dlWS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4345 fHistoElePtvsV0dlWS2 =
new THnSparseF(
"fHistoElePtvsV0dlWS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4347 fHistoElePtvsV0dlWSSide =
new THnSparseF(
"fHistoElePtvsV0dlWSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4349 fHistoElePtvsV0dlWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4351 fHistoElePtvsV0dlWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4353 fHistoElePtvsV0dlWSMix =
new THnSparseF(
"fHistoElePtvsV0dlWSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4355 fHistoElePtvsV0dlWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4357 fHistoElePtvsV0dlWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4359 fHistoElePtvsV0dlMCS =
new THnSparseF(
"fHistoElePtvsV0dlMCS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4361 fHistoElePtvsV0dlMCS1 =
new THnSparseF(
"fHistoElePtvsV0dlMCS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4363 fHistoElePtvsV0dlMCS2 =
new THnSparseF(
"fHistoElePtvsV0dlMCS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4378 Int_t bins_eleptvsv0dca[3]= {100 ,20 ,10};
4379 Double_t xmin_eleptvsv0dca[3]={0.,0. ,0.0};
4380 Double_t xmax_eleptvsv0dca[3]={10.,1. ,100};
4381 fHistoElePtvsV0dcaRS =
new THnSparseF(
"fHistoElePtvsV0dcaRS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4383 fHistoElePtvsV0dcaRS1 =
new THnSparseF(
"fHistoElePtvsV0dcaRS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4385 fHistoElePtvsV0dcaRS2 =
new THnSparseF(
"fHistoElePtvsV0dcaRS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4387 fHistoElePtvsV0dcaRSSide =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4389 fHistoElePtvsV0dcaRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4391 fHistoElePtvsV0dcaRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4393 fHistoElePtvsV0dcaRSMix =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4395 fHistoElePtvsV0dcaRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4397 fHistoElePtvsV0dcaRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4399 fHistoElePtvsV0dcaWS =
new THnSparseF(
"fHistoElePtvsV0dcaWS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4401 fHistoElePtvsV0dcaWS1 =
new THnSparseF(
"fHistoElePtvsV0dcaWS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4403 fHistoElePtvsV0dcaWS2 =
new THnSparseF(
"fHistoElePtvsV0dcaWS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4405 fHistoElePtvsV0dcaWSSide =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4407 fHistoElePtvsV0dcaWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4409 fHistoElePtvsV0dcaWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4411 fHistoElePtvsV0dcaWSMix =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4413 fHistoElePtvsV0dcaWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4415 fHistoElePtvsV0dcaWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4417 fHistoElePtvsV0dcaMCS =
new THnSparseF(
"fHistoElePtvsV0dcaMCS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4419 fHistoElePtvsV0dcaMCS1 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4421 fHistoElePtvsV0dcaMCS2 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4436 Int_t bins_elelamptvsv0dl[4]= {100 ,20 ,100, 10};
4437 Double_t xmin_elelamptvsv0dl[4]={0.,0. ,-0.5, 0.0};
4438 Double_t xmax_elelamptvsv0dl[4]={20.,40.,0.5 ,100};
4439 fHistoEleLambdaPtvsV0dlRS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4441 fHistoEleLambdaPtvsV0dlRS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4443 fHistoEleLambdaPtvsV0dlRS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4445 fHistoEleLambdaPtvsV0dlRSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4447 fHistoEleLambdaPtvsV0dlRSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4449 fHistoEleLambdaPtvsV0dlRSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4451 fHistoEleLambdaPtvsV0dlRSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4453 fHistoEleLambdaPtvsV0dlRSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4455 fHistoEleLambdaPtvsV0dlRSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4457 fHistoEleLambdaPtvsV0dlWS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4459 fHistoEleLambdaPtvsV0dlWS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4461 fHistoEleLambdaPtvsV0dlWS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4463 fHistoEleLambdaPtvsV0dlWSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4465 fHistoEleLambdaPtvsV0dlWSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4467 fHistoEleLambdaPtvsV0dlWSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4469 fHistoEleLambdaPtvsV0dlWSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4471 fHistoEleLambdaPtvsV0dlWSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4473 fHistoEleLambdaPtvsV0dlWSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4475 fHistoEleLambdaPtvsV0dlMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4477 fHistoEleLambdaPtvsV0dlMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4479 fHistoEleLambdaPtvsV0dlMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4555 fHistoEleLambdaPtvsV0dlXibMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4557 fHistoEleLambdaPtvsV0dlXibMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4559 fHistoEleLambdaPtvsV0dlXibMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4562 fHistoEleLambdaPtvsV0dlRSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4564 fHistoEleLambdaPtvsV0dlRS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4566 fHistoEleLambdaPtvsV0dlRS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4580 fHistoEleLambdaPtvsV0dlWSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4582 fHistoEleLambdaPtvsV0dlWS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4584 fHistoEleLambdaPtvsV0dlWS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4666 for(Int_t ih=0;ih<17;ih++){
4667 Int_t bins_eleptvscutvars[3];
4668 Double_t xmin_eleptvscutvars[3];
4669 Double_t xmax_eleptvscutvars[3];
4671 bins_eleptvscutvars[0] = 20;
4672 xmin_eleptvscutvars[0] = 0.;
4673 xmax_eleptvscutvars[0] = 20.;
4674 bins_eleptvscutvars[2] = 10;
4675 xmin_eleptvscutvars[2] = 0.;
4676 xmax_eleptvscutvars[2] = 100.;
4680 bins_eleptvscutvars[1] = 40;
4681 xmin_eleptvscutvars[1] = 0.;
4682 xmax_eleptvscutvars[1] = 160.;
4683 }
else if(ih==2 || ih==3){
4685 bins_eleptvscutvars[1] = 20;
4686 xmin_eleptvscutvars[1] = -5.;
4687 xmax_eleptvscutvars[1] = 5.;
4690 bins_eleptvscutvars[1] = 30;
4691 xmin_eleptvscutvars[1] = -1.5;
4692 xmax_eleptvscutvars[1] = 1.5;
4695 bins_eleptvscutvars[1] = 7;
4696 xmin_eleptvscutvars[1] = -0.5;
4697 xmax_eleptvscutvars[1] = 6.5;
4700 bins_eleptvscutvars[1] = 50;
4701 xmin_eleptvscutvars[1] = 1.1156-0.03;
4702 xmax_eleptvscutvars[1] = 1.1156+0.03;
4705 bins_eleptvscutvars[1] = 20;
4706 xmin_eleptvscutvars[1] = 0.;
4707 xmax_eleptvscutvars[1] = 5.;
4710 bins_eleptvscutvars[1] = 20;
4711 xmin_eleptvscutvars[1] = 0.;
4712 xmax_eleptvscutvars[1] = 2.;
4713 }
else if(ih==9 || ih==10 ){
4715 bins_eleptvscutvars[1] = 20;
4716 xmin_eleptvscutvars[1] = 0.;
4717 xmax_eleptvscutvars[1] = 0.5;
4720 bins_eleptvscutvars[1] = 20;
4721 xmin_eleptvscutvars[1] = 0.95;
4722 xmax_eleptvscutvars[1] = 1.0;
4725 bins_eleptvscutvars[1] = 50;
4726 xmin_eleptvscutvars[1] = 0.497-0.03;
4727 xmax_eleptvscutvars[1] = 0.497+0.03;
4728 }
else if(ih==13 || ih==14){
4730 bins_eleptvscutvars[1] = 20;
4731 xmin_eleptvscutvars[1] = -5;
4732 xmax_eleptvscutvars[1] = 5;
4735 bins_eleptvscutvars[1] = 30;
4736 xmin_eleptvscutvars[1] = -1.5;
4737 xmax_eleptvscutvars[1] = 1.5;
4740 bins_eleptvscutvars[1] = 20;
4741 xmin_eleptvscutvars[1] = 0.;
4742 xmax_eleptvscutvars[1] = 3.141592/2;
4745 fHistoElePtvsCutVarsRS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsRS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4747 fHistoElePtvsCutVarsWS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsWS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4749 fHistoElePtvsCutVarsMCS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsMCS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4753 Int_t bins_eletag[3]= {20 ,40 ,10};
4754 Double_t xmin_eletag[3]={0.,0. ,0.0};
4755 Double_t xmax_eletag[3]={10.,20 ,100};
4756 fHistoElectronPi0Total =
new THnSparseF(
"fHistoElectronPi0Total",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4758 fHistoElectronPi0Tag =
new THnSparseF(
"fHistoElectronPi0Tag",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4760 fHistoElectronEtaTotal =
new THnSparseF(
"fHistoElectronEtaTotal",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4762 fHistoElectronEtaTag =
new THnSparseF(
"fHistoElectronEtaTag",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4776 if(!part)
return 0x0;
4777 if(!aod)
return 0x0;
4782 AliAODVertex *primVertexAOD;
4783 Bool_t unsetvtx = kFALSE;
4787 primVertexAOD =
fVtx1;
4792 primVertexAOD =
fVtx1;
4794 if(!primVertexAOD)
return 0x0;
4795 Double_t posprim[3]; primVertexAOD->GetXYZ(posprim);
4800 AliESDtrack *esdtrack =
new AliESDtrack((AliVTrack*)partpid);
4802 AliNeutralTrackParam *trackV0=NULL;
4803 const AliVTrack *trackVV0 =
dynamic_cast<const AliVTrack*
>(v0);
4804 if(trackVV0) trackV0 =
new AliNeutralTrackParam(trackVV0);
4806 Double_t xdummy, ydummy;
4807 Double_t dca = esdtrack->GetDCA(trackV0,
fBzkG,xdummy,ydummy);
4814 Double_t d0z0bach[2],covd0z0bach[3];
4815 if(sqrt(pow(secVert->GetX(),2)+pow(secVert->GetY(),2))<1.){
4816 part->PropagateToDCA(secVert,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4817 trackV0->PropagateToDCA(secVert,
fBzkG,kVeryBig);
4819 part->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4820 trackV0->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig);
4822 Double_t momv0_new[3]={-9999,-9999,-9999.};
4823 trackV0->GetPxPyPz(momv0_new);
4825 Double_t px[2],py[2],pz[2];
4826 px[0] = part->Px(); py[0] = part->Py(); pz[0] = part->Pz();
4827 px[1] = momv0_new[0]; py[1] = momv0_new[1]; pz[1] = momv0_new[2];
4832 Double_t d0[3],d0err[3];
4834 part->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4836 d0err[0] = TMath::Sqrt(covd0z0bach[0]);
4838 Double_t d0z0v0[2],covd0z0v0[3];
4839 trackV0->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0v0,covd0z0v0);
4841 d0err[1] = TMath::Sqrt(covd0z0v0[0]);
4846 Short_t
charge = part->Charge();
4850 if(unsetvtx)
delete primVertexAOD; primVertexAOD=NULL;
4851 if(esdtrack)
delete esdtrack;
4852 if(trackV0)
delete trackV0;
4856 UShort_t
id[2]={(UShort_t)part->GetID(),(UShort_t)trackV0->GetID()};
4859 theCascade->GetSecondaryVtx()->AddDaughter(part);
4860 theCascade->GetSecondaryVtx()->AddDaughter(v0);
4862 if(unsetvtx)
delete primVertexAOD; primVertexAOD=NULL;
4863 if(esdtrack)
delete esdtrack;
4864 if(trackV0)
delete trackV0;
4877 TObjArray *TrackArray =
new TObjArray(3);
4879 AliESDtrack *cptrk1 =
new AliESDtrack((AliVTrack*)trk);
4880 TrackArray->AddAt(cptrk1,0);
4882 AliESDtrack *cascptrack =
new AliESDtrack((AliVTrack*)v0->GetDaughter(0));
4883 TrackArray->AddAt(cascptrack,1);
4884 AliESDtrack *cascntrack =
new AliESDtrack((AliVTrack*)v0->GetDaughter(1));
4885 TrackArray->AddAt(cascntrack,2);
4889 for(Int_t i=0;i<3;i++)
4891 AliESDtrack *tesd = (AliESDtrack*)TrackArray->UncheckedAt(i);
4894 TrackArray->Clear();
<