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"
104 fHNTrackletCorrvsZ(0),
106 fIsEventSelected(kFALSE),
107 fWriteVariableTree(kFALSE),
108 fWriteEachVariableTree(kFALSE),
109 fWriteMCVariableTree(kFALSE),
111 fEleVariablesTree(0),
113 fSingleVariablesTree(0),
115 fMCEleVariablesTree(0),
116 fMCV0VariablesTree(0),
117 fMCGenPairVariablesTree(0),
118 fCorrelationVariablesTree(0),
119 fReconstructPrimVert(kFALSE),
125 fCandidateVariables(),
126 fCandidateEleVariables(),
127 fCandidateV0Variables(),
128 fCandidateSingleVariables(),
129 fCandidateMCVariables(),
130 fCandidateMCEleVariables(),
131 fCandidateMCV0Variables(),
132 fCandidateMCGenPairVariables(),
133 fCorrelationVariables(),
145 fUseCentralityV0M(kFALSE),
146 fUseCentralitySPDTracklet(kFALSE),
150 fMCDoPairAnalysis(kFALSE),
151 fHistoEleLambdaMass(0),
152 fHistoEleLambdaMassRS(0),
153 fHistoEleLambdaMassRS1(0),
154 fHistoEleLambdaMassRS2(0),
155 fHistoEleLambdaMassWS(0),
156 fHistoEleLambdaMassWS1(0),
157 fHistoEleLambdaMassWS2(0),
158 fHistoEleLambdaMassRSMix(0),
159 fHistoEleLambdaMassRSMix1(0),
160 fHistoEleLambdaMassRSMix2(0),
161 fHistoEleLambdaMassWSMix(0),
162 fHistoEleLambdaMassWSMix1(0),
163 fHistoEleLambdaMassWSMix2(0),
164 fHistoEleLambdaMassRSSide(0),
165 fHistoEleLambdaMassRSSide1(0),
166 fHistoEleLambdaMassRSSide2(0),
167 fHistoEleLambdaMassWSSide(0),
168 fHistoEleLambdaMassWSSide1(0),
169 fHistoEleLambdaMassWSSide2(0),
170 fHistoEleLambdaMassRSAway(0),
171 fHistoEleLambdaMassRS1Away(0),
172 fHistoEleLambdaMassRS2Away(0),
173 fHistoEleLambdaMassWSAway(0),
174 fHistoEleLambdaMassWS1Away(0),
175 fHistoEleLambdaMassWS2Away(0),
176 fHistoEleLambdaMassRSMixAway(0),
177 fHistoEleLambdaMassRSMix1Away(0),
178 fHistoEleLambdaMassRSMix2Away(0),
179 fHistoEleLambdaMassWSMixAway(0),
180 fHistoEleLambdaMassWSMix1Away(0),
181 fHistoEleLambdaMassWSMix2Away(0),
182 fHistoEleLambdaMassRSSideAway(0),
183 fHistoEleLambdaMassRSSide1Away(0),
184 fHistoEleLambdaMassRSSide2Away(0),
185 fHistoEleLambdaMassWSSideAway(0),
186 fHistoEleLambdaMassWSSide1Away(0),
187 fHistoEleLambdaMassWSSide2Away(0),
188 fHistoEleLambdaMassvsElePtRS(0),
189 fHistoEleLambdaMassvsElePtWS(0),
190 fHistoEleLambdaMassvsElePtRSMix(0),
191 fHistoEleLambdaMassvsElePtWSMix(0),
192 fHistoEleLambdaMassvsElePtRSSide(0),
193 fHistoEleLambdaMassvsElePtWSSide(0),
194 fHistoEleLambdaMassvsElePtRS1(0),
195 fHistoEleLambdaMassvsElePtWS1(0),
196 fHistoEleLambdaMassvsElePtRSMix1(0),
197 fHistoEleLambdaMassvsElePtWSMix1(0),
198 fHistoEleLambdaMassvsElePtRSSide1(0),
199 fHistoEleLambdaMassvsElePtWSSide1(0),
200 fHistoEleLambdaMassvsElePtRS2(0),
201 fHistoEleLambdaMassvsElePtWS2(0),
202 fHistoEleLambdaMassvsElePtRSMix2(0),
203 fHistoEleLambdaMassvsElePtWSMix2(0),
204 fHistoEleLambdaMassvsElePtRSSide2(0),
205 fHistoEleLambdaMassvsElePtWSSide2(0),
210 fHistoEleLambdaMassMCS(0),
211 fHistoEleLambdaMassLambdabMCS(0),
212 fHistoEleLambdaMassXibMCS(0),
213 fHistoEleLambdaMassMCGen(0),
214 fHistoEleLambdaMassvsElePtMCS(0),
215 fHistoEleLambdaMassvsElePtMCGen(0),
216 fHistoEleLambdaMassvsElePtMCS1(0),
217 fHistoEleLambdaMassvsElePtMCGen1(0),
218 fHistoEleLambdaMassvsElePtMCS2(0),
219 fHistoEleLambdaMassvsElePtMCGen2(0),
222 fHistoElePtvsEtaRS(0),
223 fHistoElePtvsEtaWS(0),
224 fHistoElePtvsEtaRSMix(0),
225 fHistoElePtvsEtaWSMix(0),
226 fHistoElePtvsEtaMCS(0),
227 fHistoElePtvsEtaMCGen(0),
228 fHistoElePtvsLambdaPtRS(0),
229 fHistoElePtvsLambdaPtWS(0),
230 fHistoElePtvsLambdaPtRSMix(0),
231 fHistoElePtvsLambdaPtWSMix(0),
232 fHistoElePtvsLambdaPtMCS(0),
233 fHistoElePtvsLambdaPtvsLcPtMCS(0),
234 fHistoElePtvsLambdaPtMCGen(0),
235 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
236 fHistoElePtvsLambdaPtMCLcGen(0),
237 fHistoElePtvsd0RS(0),
238 fHistoElePtvsd0WS(0),
239 fHistoElePtvsd0RSMix(0),
240 fHistoElePtvsd0WSMix(0),
241 fHistoElePtvsd0MCS(0),
242 fHistoElePtvsd0PromptMCS(0),
243 fHistoElePtvsd0BFeeddownMCS(0),
244 fHistoEleLambdaMassFeeddownXic0MCS(0),
245 fHistoEleLambdaMassFeeddownXic0MCGen(0),
246 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
247 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
248 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
249 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
250 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
251 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
252 fHistoElePtFeeddownXic0MCS(0),
253 fHistoElePtFeeddownXic0MCGen(0),
254 fHistoElePtvsEtaFeeddownXic0MCS(0),
255 fHistoElePtvsEtaFeeddownXic0MCGen(0),
256 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
257 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
258 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
259 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
260 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
261 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
262 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
263 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
264 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
265 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
266 fHistoElePtFeeddownXicPlusMCS(0),
267 fHistoElePtFeeddownXicPlusMCGen(0),
268 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
269 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
270 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
271 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
274 fHistoBachPtMCGen(0),
276 fHistoLambdaMassvsPt(0),
277 fHistoLambdaMassvsPtMCS(0),
278 fHistoLambdaMassvsPtMCGen(0),
279 fHistoLambdaPtvsEta(0),
280 fHistoLambdaPtvsEtaMCS(0),
281 fHistoLambdaPtvsEtaMCGen(0),
282 fHistoLambdaArmenteros(0),
283 fHistoLambdaArmenterosMCS(0),
284 fHistoLambdaPtvsDl(0),
285 fHistoLambdaPtvsDlSide(0),
286 fHistoLambdaPtvsDlMCS(0),
287 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
288 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
289 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
290 fHistoK0sMassvsPt(0),
291 fHistoElectronTPCPID(0),
292 fHistoElectronTOFPID(0),
293 fHistoElectronTPCSelPID(0),
294 fHistoElectronTOFSelPID(0),
295 fHistoElectronTPCPIDSelTOF(0),
296 fHistoElectronTPCPIDSelTOFSmallEta(0),
297 fHistoElectronTPCPIDSelTOFLargeEta(0),
298 fHistoMassConversionsMin(0),
299 fHistoMassConversionsSameSignMin(0),
300 fHistoElectronQovPtvsPhi(0),
301 fHistoLambdaQovPtvsPhi(0),
308 fHistoLcMCGenWithEleA(0),
309 fHistoLcMCGenWithEleB(0),
310 fHistoLbMCGenWithEleA(0),
311 fHistoLbMCGenWithEleB(0),
312 fHistoLbMCGenWithLambdac(0),
319 fHistoFeedDownXic0MCGen(0),
320 fHistoFeedDownXic0MCGen1(0),
321 fHistoFeedDownXic0MCGen2(0),
322 fHistoFeedDownXic0MCS(0),
323 fHistoFeedDownXic0MCS1(0),
324 fHistoFeedDownXic0MCS2(0),
325 fHistoFeedDownXicPlusMCGen(0),
326 fHistoFeedDownXicPlusMCGen1(0),
327 fHistoFeedDownXicPlusMCGen2(0),
328 fHistoFeedDownXicPlusMCS(0),
329 fHistoFeedDownXicPlusMCS1(0),
330 fHistoFeedDownXicPlusMCS2(0),
331 fHistoFeedDownSigmaMCGen(0),
332 fHistoFeedDownSigmaMCGen1(0),
333 fHistoFeedDownSigmaMCGen2(0),
334 fHistoFeedDownSigmaMCS(0),
335 fHistoFeedDownSigmaMCS1(0),
336 fHistoFeedDownSigmaMCS2(0),
337 fHistoLcElectronMCGen(0),
338 fHistoLcElectronMCGen1(0),
339 fHistoLcElectronMCGen2(0),
340 fHistoLcElectronMCS(0),
341 fHistoLcElectronMCS1(0),
342 fHistoLcElectronMCS2(0),
343 fHistoLbElectronMCGen(0),
344 fHistoLbElectronMCGen1(0),
345 fHistoLbElectronMCGen2(0),
346 fHistoLbElectronMCS(0),
347 fHistoLbElectronMCS1(0),
348 fHistoLbElectronMCS2(0),
349 fHistoElectronFeedDownXic0MCGen(0),
350 fHistoElectronFeedDownXic0MCGen1(0),
351 fHistoElectronFeedDownXic0MCGen2(0),
352 fHistoElectronFeedDownXic0MCS(0),
353 fHistoElectronFeedDownXic0MCS1(0),
354 fHistoElectronFeedDownXic0MCS2(0),
355 fHistoElectronFeedDownXicPlusMCGen(0),
356 fHistoElectronFeedDownXicPlusMCGen1(0),
357 fHistoElectronFeedDownXicPlusMCGen2(0),
358 fHistoElectronFeedDownXicPlusMCS(0),
359 fHistoElectronFeedDownXicPlusMCS1(0),
360 fHistoElectronFeedDownXicPlusMCS2(0),
361 fHistoElectronFeedDownXibMCGen(0),
362 fHistoElectronFeedDownXibMCGen1(0),
363 fHistoElectronFeedDownXibMCGen2(0),
364 fHistoElectronFeedDownXibMCS(0),
365 fHistoElectronFeedDownXibMCS1(0),
366 fHistoElectronFeedDownXibMCS2(0),
367 fHistoElectronMCGen(0),
368 fHistoBottomElectronMCGen(0),
369 fHistoCharmElectronMCGen(0),
370 fHistoLambdaMCGen(0),
371 fHistoElePtvsV0dlRS(0),
372 fHistoElePtvsV0dlRS1(0),
373 fHistoElePtvsV0dlRS2(0),
374 fHistoElePtvsV0dlRSSide(0),
375 fHistoElePtvsV0dlRSSide1(0),
376 fHistoElePtvsV0dlRSSide2(0),
377 fHistoElePtvsV0dlRSMix(0),
378 fHistoElePtvsV0dlRSMix1(0),
379 fHistoElePtvsV0dlRSMix2(0),
380 fHistoElePtvsV0dlWS(0),
381 fHistoElePtvsV0dlWS1(0),
382 fHistoElePtvsV0dlWS2(0),
383 fHistoElePtvsV0dlWSSide(0),
384 fHistoElePtvsV0dlWSSide1(0),
385 fHistoElePtvsV0dlWSSide2(0),
386 fHistoElePtvsV0dlWSMix(0),
387 fHistoElePtvsV0dlWSMix1(0),
388 fHistoElePtvsV0dlWSMix2(0),
389 fHistoElePtvsV0dlMCS(0),
390 fHistoElePtvsV0dlMCS1(0),
391 fHistoElePtvsV0dlMCS2(0),
392 fHistoElePtvsV0dlFeedDownXic0MCS(0),
393 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
394 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
395 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
396 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
397 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
398 fHistoElePtvsV0dcaRS(0),
399 fHistoElePtvsV0dcaRS1(0),
400 fHistoElePtvsV0dcaRS2(0),
401 fHistoElePtvsV0dcaRSSide(0),
402 fHistoElePtvsV0dcaRSSide1(0),
403 fHistoElePtvsV0dcaRSSide2(0),
404 fHistoElePtvsV0dcaRSMix(0),
405 fHistoElePtvsV0dcaRSMix1(0),
406 fHistoElePtvsV0dcaRSMix2(0),
407 fHistoElePtvsV0dcaWS(0),
408 fHistoElePtvsV0dcaWS1(0),
409 fHistoElePtvsV0dcaWS2(0),
410 fHistoElePtvsV0dcaWSSide(0),
411 fHistoElePtvsV0dcaWSSide1(0),
412 fHistoElePtvsV0dcaWSSide2(0),
413 fHistoElePtvsV0dcaWSMix(0),
414 fHistoElePtvsV0dcaWSMix1(0),
415 fHistoElePtvsV0dcaWSMix2(0),
416 fHistoElePtvsV0dcaMCS(0),
417 fHistoElePtvsV0dcaMCS1(0),
418 fHistoElePtvsV0dcaMCS2(0),
419 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
420 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
421 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
422 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
423 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
424 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
425 fHistoEleLambdaPtvsV0dlRS(0),
426 fHistoEleLambdaPtvsV0dlRS1(0),
427 fHistoEleLambdaPtvsV0dlRS2(0),
428 fHistoEleLambdaPtvsV0dlRSSide(0),
429 fHistoEleLambdaPtvsV0dlRSSide1(0),
430 fHistoEleLambdaPtvsV0dlRSSide2(0),
431 fHistoEleLambdaPtvsV0dlRSMix(0),
432 fHistoEleLambdaPtvsV0dlRSMix1(0),
433 fHistoEleLambdaPtvsV0dlRSMix2(0),
434 fHistoEleLambdaPtvsV0dlWS(0),
435 fHistoEleLambdaPtvsV0dlWS1(0),
436 fHistoEleLambdaPtvsV0dlWS2(0),
437 fHistoEleLambdaPtvsV0dlWSSide(0),
438 fHistoEleLambdaPtvsV0dlWSSide1(0),
439 fHistoEleLambdaPtvsV0dlWSSide2(0),
440 fHistoEleLambdaPtvsV0dlWSMix(0),
441 fHistoEleLambdaPtvsV0dlWSMix1(0),
442 fHistoEleLambdaPtvsV0dlWSMix2(0),
443 fHistoEleLambdaPtvsV0dlMCS(0),
444 fHistoEleLambdaPtvsV0dlMCS1(0),
445 fHistoEleLambdaPtvsV0dlMCS2(0),
446 fHistoEleLambdaPtvsV0dlPromptMCS(0),
447 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
448 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
449 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
450 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
451 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
452 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
453 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
454 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
455 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
456 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
457 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
458 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
459 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
460 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
461 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
462 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
463 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
464 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
465 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
466 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
467 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
468 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
469 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
470 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
471 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
472 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
473 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
474 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
475 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
476 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
477 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
478 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
479 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
480 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
481 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
482 fHistoEleLambdaPtvsV0dlXibMCS(0),
483 fHistoEleLambdaPtvsV0dlXibMCS1(0),
484 fHistoEleLambdaPtvsV0dlXibMCS2(0),
485 fHistoEleLambdaPtvsV0dlRSAway(0),
486 fHistoEleLambdaPtvsV0dlRS1Away(0),
487 fHistoEleLambdaPtvsV0dlRS2Away(0),
488 fHistoEleLambdaPtvsV0dlRSSideAway(0),
489 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
490 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
491 fHistoEleLambdaPtvsV0dlRSMixAway(0),
492 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
493 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
494 fHistoEleLambdaPtvsV0dlWSAway(0),
495 fHistoEleLambdaPtvsV0dlWS1Away(0),
496 fHistoEleLambdaPtvsV0dlWS2Away(0),
497 fHistoEleLambdaPtvsV0dlWSSideAway(0),
498 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
499 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
500 fHistoEleLambdaPtvsV0dlWSMixAway(0),
501 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
502 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
503 fHistoCorrelationVariablesvsEleLambdaPt(0),
504 fHistoCorrelationVariablesvsEleLambdaPtMix(0),
505 fHistoCorrelationVariablesvsEleLambdaPtMC(0),
506 fHistoCorrelationVariablesvsElePt(0),
507 fHistoCorrelationVariablesvsElePtMix(0),
508 fHistoCorrelationVariablesvsElePtMC(0),
509 fHistoCorrelationVariablesvsLambdaPt(0),
510 fHistoCorrelationVariablesvsLambdaPtMix(0),
511 fHistoCorrelationVariablesvsLambdaPtMC(0),
512 fHistoMassVariablesvsEleLambdaPt(0),
513 fHistoMassVariablesvsEleLambdaPtMix(0),
514 fHistoMassVariablesvsEleLambdaPtMC(0),
515 fHistoMassVariablesvsElePt(0),
516 fHistoMassVariablesvsElePtMix(0),
517 fHistoMassVariablesvsElePtMC(0),
518 fHistoMassVariablesvsLambdaPt(0),
519 fHistoMassVariablesvsLambdaPtMix(0),
520 fHistoMassVariablesvsLambdaPtMC(0),
521 fHistoResponseElePt(0),
522 fHistoResponseElePt1(0),
523 fHistoResponseElePt2(0),
524 fHistoResponseEleLambdaPt(0),
525 fHistoResponseEleLambdaPt1(0),
526 fHistoResponseEleLambdaPt2(0),
527 fHistoResponseEleLambdaPtFeeddownXic0(0),
528 fHistoResponseEleLambdaPtFeeddownXic01(0),
529 fHistoResponseEleLambdaPtFeeddownXic02(0),
530 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
531 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
532 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
533 fHistoResponseEleLambdaPtFeeddownSigma(0),
534 fHistoResponseEleLambdaPtFeeddownSigma1(0),
535 fHistoResponseEleLambdaPtFeeddownSigma2(0),
536 fHistoResponseEleLambdaPtLambdab(0),
537 fHistoResponseElePtLambdab(0),
538 fHistoResponseLambdaPtLambdab(0),
539 fHistoResponseLambdaPt(0),
540 fHistoResponseLambdaPtFeeddownXic0(0),
541 fHistoResponseLambdaPtFeeddownXicPlus(0),
542 fHistoResponseLambdaPtFeeddownSigma(0),
543 fHistoResponseLambdaPtvsEleLambdaPt(0),
544 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
545 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
546 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
547 fHistoEleLambdaPtvsRapidityRS(0),
548 fHistoEleLambdaPtvsRapidityWS(0),
549 fHistoEleLambdaPtvsRapidityMCS(0),
551 fHistoElectronPi0Total(0),
552 fHistoElectronPi0Tag(0),
554 fHistoElectronEtaTotal(0),
555 fHistoElectronEtaTag(0),
557 fHistonEvtvsRunNumber(0),
558 fHistonElevsRunNumber(0),
559 fHistonLambdavsRunNumber(0),
560 fHistoMCEventType(0),
561 fHistoMCDeltaPhiccbar(0),
563 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
564 fHistodPhiSdEtaSElectronProtonR125RS(0),
565 fHistodPhiSdEtaSElectronProtonR125WS(0),
566 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
567 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
568 fHistodPhiSdEtaSElectronPionR125RS(0),
569 fHistodPhiSdEtaSElectronPionR125WS(0),
570 fHistodPhiSdEtaSElectronPionR125RSMix(0),
571 fHistodPhiSdEtaSElectronPionR125WSMix(0),
573 fMixWithoutConversionFlag(kFALSE),
574 fNumberOfEventsForMixing (5),
592 for(
Int_t i=0;i<17;i++){
597 for(
Int_t i=0;i<8;i++){
606 Bool_t writeVariableTree) :
617 fHNTrackletCorrvsZ(0),
619 fIsEventSelected(kFALSE),
620 fWriteVariableTree(writeVariableTree),
621 fWriteEachVariableTree(kFALSE),
622 fWriteMCVariableTree(kFALSE),
624 fEleVariablesTree(0),
626 fSingleVariablesTree(0),
628 fMCEleVariablesTree(0),
629 fMCV0VariablesTree(0),
630 fMCGenPairVariablesTree(0),
631 fCorrelationVariablesTree(0),
632 fReconstructPrimVert(kFALSE),
638 fCandidateVariables(),
639 fCandidateEleVariables(),
640 fCandidateV0Variables(),
641 fCandidateSingleVariables(),
642 fCandidateMCVariables(),
643 fCandidateMCEleVariables(),
644 fCandidateMCV0Variables(),
645 fCandidateMCGenPairVariables(),
646 fCorrelationVariables(),
658 fUseCentralityV0M(kFALSE),
659 fUseCentralitySPDTracklet(kFALSE),
663 fMCDoPairAnalysis(kFALSE),
664 fHistoEleLambdaMass(0),
665 fHistoEleLambdaMassRS(0),
666 fHistoEleLambdaMassRS1(0),
667 fHistoEleLambdaMassRS2(0),
668 fHistoEleLambdaMassWS(0),
669 fHistoEleLambdaMassWS1(0),
670 fHistoEleLambdaMassWS2(0),
671 fHistoEleLambdaMassRSMix(0),
672 fHistoEleLambdaMassRSMix1(0),
673 fHistoEleLambdaMassRSMix2(0),
674 fHistoEleLambdaMassWSMix(0),
675 fHistoEleLambdaMassWSMix1(0),
676 fHistoEleLambdaMassWSMix2(0),
677 fHistoEleLambdaMassRSSide(0),
678 fHistoEleLambdaMassRSSide1(0),
679 fHistoEleLambdaMassRSSide2(0),
680 fHistoEleLambdaMassWSSide(0),
681 fHistoEleLambdaMassWSSide1(0),
682 fHistoEleLambdaMassWSSide2(0),
683 fHistoEleLambdaMassRSAway(0),
684 fHistoEleLambdaMassRS1Away(0),
685 fHistoEleLambdaMassRS2Away(0),
686 fHistoEleLambdaMassWSAway(0),
687 fHistoEleLambdaMassWS1Away(0),
688 fHistoEleLambdaMassWS2Away(0),
689 fHistoEleLambdaMassRSMixAway(0),
690 fHistoEleLambdaMassRSMix1Away(0),
691 fHistoEleLambdaMassRSMix2Away(0),
692 fHistoEleLambdaMassWSMixAway(0),
693 fHistoEleLambdaMassWSMix1Away(0),
694 fHistoEleLambdaMassWSMix2Away(0),
695 fHistoEleLambdaMassRSSideAway(0),
696 fHistoEleLambdaMassRSSide1Away(0),
697 fHistoEleLambdaMassRSSide2Away(0),
698 fHistoEleLambdaMassWSSideAway(0),
699 fHistoEleLambdaMassWSSide1Away(0),
700 fHistoEleLambdaMassWSSide2Away(0),
701 fHistoEleLambdaMassvsElePtRS(0),
702 fHistoEleLambdaMassvsElePtWS(0),
703 fHistoEleLambdaMassvsElePtRSMix(0),
704 fHistoEleLambdaMassvsElePtWSMix(0),
705 fHistoEleLambdaMassvsElePtRSSide(0),
706 fHistoEleLambdaMassvsElePtWSSide(0),
707 fHistoEleLambdaMassvsElePtRS1(0),
708 fHistoEleLambdaMassvsElePtWS1(0),
709 fHistoEleLambdaMassvsElePtRSMix1(0),
710 fHistoEleLambdaMassvsElePtWSMix1(0),
711 fHistoEleLambdaMassvsElePtRSSide1(0),
712 fHistoEleLambdaMassvsElePtWSSide1(0),
713 fHistoEleLambdaMassvsElePtRS2(0),
714 fHistoEleLambdaMassvsElePtWS2(0),
715 fHistoEleLambdaMassvsElePtRSMix2(0),
716 fHistoEleLambdaMassvsElePtWSMix2(0),
717 fHistoEleLambdaMassvsElePtRSSide2(0),
718 fHistoEleLambdaMassvsElePtWSSide2(0),
723 fHistoEleLambdaMassMCS(0),
724 fHistoEleLambdaMassLambdabMCS(0),
725 fHistoEleLambdaMassXibMCS(0),
726 fHistoEleLambdaMassMCGen(0),
727 fHistoEleLambdaMassvsElePtMCS(0),
728 fHistoEleLambdaMassvsElePtMCGen(0),
729 fHistoEleLambdaMassvsElePtMCS1(0),
730 fHistoEleLambdaMassvsElePtMCGen1(0),
731 fHistoEleLambdaMassvsElePtMCS2(0),
732 fHistoEleLambdaMassvsElePtMCGen2(0),
735 fHistoElePtvsEtaRS(0),
736 fHistoElePtvsEtaWS(0),
737 fHistoElePtvsEtaRSMix(0),
738 fHistoElePtvsEtaWSMix(0),
739 fHistoElePtvsEtaMCS(0),
740 fHistoElePtvsEtaMCGen(0),
741 fHistoElePtvsLambdaPtRS(0),
742 fHistoElePtvsLambdaPtWS(0),
743 fHistoElePtvsLambdaPtRSMix(0),
744 fHistoElePtvsLambdaPtWSMix(0),
745 fHistoElePtvsLambdaPtMCS(0),
746 fHistoElePtvsLambdaPtvsLcPtMCS(0),
747 fHistoElePtvsLambdaPtMCGen(0),
748 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
749 fHistoElePtvsLambdaPtMCLcGen(0),
750 fHistoElePtvsd0RS(0),
751 fHistoElePtvsd0WS(0),
752 fHistoElePtvsd0RSMix(0),
753 fHistoElePtvsd0WSMix(0),
754 fHistoElePtvsd0MCS(0),
755 fHistoElePtvsd0PromptMCS(0),
756 fHistoElePtvsd0BFeeddownMCS(0),
757 fHistoEleLambdaMassFeeddownXic0MCS(0),
758 fHistoEleLambdaMassFeeddownXic0MCGen(0),
759 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
760 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
761 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
762 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
763 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
764 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
765 fHistoElePtFeeddownXic0MCS(0),
766 fHistoElePtFeeddownXic0MCGen(0),
767 fHistoElePtvsEtaFeeddownXic0MCS(0),
768 fHistoElePtvsEtaFeeddownXic0MCGen(0),
769 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
770 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
771 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
772 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
773 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
774 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
775 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
776 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
777 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
778 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
779 fHistoElePtFeeddownXicPlusMCS(0),
780 fHistoElePtFeeddownXicPlusMCGen(0),
781 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
782 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
783 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
784 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
787 fHistoBachPtMCGen(0),
789 fHistoLambdaMassvsPt(0),
790 fHistoLambdaMassvsPtMCS(0),
791 fHistoLambdaMassvsPtMCGen(0),
792 fHistoLambdaPtvsEta(0),
793 fHistoLambdaPtvsEtaMCS(0),
794 fHistoLambdaPtvsEtaMCGen(0),
795 fHistoLambdaArmenteros(0),
796 fHistoLambdaArmenterosMCS(0),
797 fHistoLambdaPtvsDl(0),
798 fHistoLambdaPtvsDlSide(0),
799 fHistoLambdaPtvsDlMCS(0),
800 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
801 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
802 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
803 fHistoK0sMassvsPt(0),
804 fHistoElectronTPCPID(0),
805 fHistoElectronTOFPID(0),
806 fHistoElectronTPCSelPID(0),
807 fHistoElectronTOFSelPID(0),
808 fHistoElectronTPCPIDSelTOF(0),
809 fHistoElectronTPCPIDSelTOFSmallEta(0),
810 fHistoElectronTPCPIDSelTOFLargeEta(0),
811 fHistoMassConversionsMin(0),
812 fHistoMassConversionsSameSignMin(0),
813 fHistoElectronQovPtvsPhi(0),
814 fHistoLambdaQovPtvsPhi(0),
821 fHistoLcMCGenWithEleA(0),
822 fHistoLcMCGenWithEleB(0),
823 fHistoLbMCGenWithEleA(0),
824 fHistoLbMCGenWithEleB(0),
825 fHistoLbMCGenWithLambdac(0),
832 fHistoFeedDownXic0MCGen(0),
833 fHistoFeedDownXic0MCGen1(0),
834 fHistoFeedDownXic0MCGen2(0),
835 fHistoFeedDownXic0MCS(0),
836 fHistoFeedDownXic0MCS1(0),
837 fHistoFeedDownXic0MCS2(0),
838 fHistoFeedDownXicPlusMCGen(0),
839 fHistoFeedDownXicPlusMCGen1(0),
840 fHistoFeedDownXicPlusMCGen2(0),
841 fHistoFeedDownXicPlusMCS(0),
842 fHistoFeedDownXicPlusMCS1(0),
843 fHistoFeedDownXicPlusMCS2(0),
844 fHistoFeedDownSigmaMCGen(0),
845 fHistoFeedDownSigmaMCGen1(0),
846 fHistoFeedDownSigmaMCGen2(0),
847 fHistoFeedDownSigmaMCS(0),
848 fHistoFeedDownSigmaMCS1(0),
849 fHistoFeedDownSigmaMCS2(0),
850 fHistoLcElectronMCGen(0),
851 fHistoLcElectronMCGen1(0),
852 fHistoLcElectronMCGen2(0),
853 fHistoLcElectronMCS(0),
854 fHistoLcElectronMCS1(0),
855 fHistoLcElectronMCS2(0),
856 fHistoLbElectronMCGen(0),
857 fHistoLbElectronMCGen1(0),
858 fHistoLbElectronMCGen2(0),
859 fHistoLbElectronMCS(0),
860 fHistoLbElectronMCS1(0),
861 fHistoLbElectronMCS2(0),
862 fHistoElectronFeedDownXic0MCGen(0),
863 fHistoElectronFeedDownXic0MCGen1(0),
864 fHistoElectronFeedDownXic0MCGen2(0),
865 fHistoElectronFeedDownXic0MCS(0),
866 fHistoElectronFeedDownXic0MCS1(0),
867 fHistoElectronFeedDownXic0MCS2(0),
868 fHistoElectronFeedDownXicPlusMCGen(0),
869 fHistoElectronFeedDownXicPlusMCGen1(0),
870 fHistoElectronFeedDownXicPlusMCGen2(0),
871 fHistoElectronFeedDownXicPlusMCS(0),
872 fHistoElectronFeedDownXicPlusMCS1(0),
873 fHistoElectronFeedDownXicPlusMCS2(0),
874 fHistoElectronFeedDownXibMCGen(0),
875 fHistoElectronFeedDownXibMCGen1(0),
876 fHistoElectronFeedDownXibMCGen2(0),
877 fHistoElectronFeedDownXibMCS(0),
878 fHistoElectronFeedDownXibMCS1(0),
879 fHistoElectronFeedDownXibMCS2(0),
880 fHistoElectronMCGen(0),
881 fHistoBottomElectronMCGen(0),
882 fHistoCharmElectronMCGen(0),
883 fHistoLambdaMCGen(0),
884 fHistoElePtvsV0dlRS(0),
885 fHistoElePtvsV0dlRS1(0),
886 fHistoElePtvsV0dlRS2(0),
887 fHistoElePtvsV0dlRSSide(0),
888 fHistoElePtvsV0dlRSSide1(0),
889 fHistoElePtvsV0dlRSSide2(0),
890 fHistoElePtvsV0dlRSMix(0),
891 fHistoElePtvsV0dlRSMix1(0),
892 fHistoElePtvsV0dlRSMix2(0),
893 fHistoElePtvsV0dlWS(0),
894 fHistoElePtvsV0dlWS1(0),
895 fHistoElePtvsV0dlWS2(0),
896 fHistoElePtvsV0dlWSSide(0),
897 fHistoElePtvsV0dlWSSide1(0),
898 fHistoElePtvsV0dlWSSide2(0),
899 fHistoElePtvsV0dlWSMix(0),
900 fHistoElePtvsV0dlWSMix1(0),
901 fHistoElePtvsV0dlWSMix2(0),
902 fHistoElePtvsV0dlMCS(0),
903 fHistoElePtvsV0dlMCS1(0),
904 fHistoElePtvsV0dlMCS2(0),
905 fHistoElePtvsV0dlFeedDownXic0MCS(0),
906 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
907 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
908 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
909 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
910 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
911 fHistoElePtvsV0dcaRS(0),
912 fHistoElePtvsV0dcaRS1(0),
913 fHistoElePtvsV0dcaRS2(0),
914 fHistoElePtvsV0dcaRSSide(0),
915 fHistoElePtvsV0dcaRSSide1(0),
916 fHistoElePtvsV0dcaRSSide2(0),
917 fHistoElePtvsV0dcaRSMix(0),
918 fHistoElePtvsV0dcaRSMix1(0),
919 fHistoElePtvsV0dcaRSMix2(0),
920 fHistoElePtvsV0dcaWS(0),
921 fHistoElePtvsV0dcaWS1(0),
922 fHistoElePtvsV0dcaWS2(0),
923 fHistoElePtvsV0dcaWSSide(0),
924 fHistoElePtvsV0dcaWSSide1(0),
925 fHistoElePtvsV0dcaWSSide2(0),
926 fHistoElePtvsV0dcaWSMix(0),
927 fHistoElePtvsV0dcaWSMix1(0),
928 fHistoElePtvsV0dcaWSMix2(0),
929 fHistoElePtvsV0dcaMCS(0),
930 fHistoElePtvsV0dcaMCS1(0),
931 fHistoElePtvsV0dcaMCS2(0),
932 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
933 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
934 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
935 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
936 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
937 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
938 fHistoEleLambdaPtvsV0dlRS(0),
939 fHistoEleLambdaPtvsV0dlRS1(0),
940 fHistoEleLambdaPtvsV0dlRS2(0),
941 fHistoEleLambdaPtvsV0dlRSSide(0),
942 fHistoEleLambdaPtvsV0dlRSSide1(0),
943 fHistoEleLambdaPtvsV0dlRSSide2(0),
944 fHistoEleLambdaPtvsV0dlRSMix(0),
945 fHistoEleLambdaPtvsV0dlRSMix1(0),
946 fHistoEleLambdaPtvsV0dlRSMix2(0),
947 fHistoEleLambdaPtvsV0dlWS(0),
948 fHistoEleLambdaPtvsV0dlWS1(0),
949 fHistoEleLambdaPtvsV0dlWS2(0),
950 fHistoEleLambdaPtvsV0dlWSSide(0),
951 fHistoEleLambdaPtvsV0dlWSSide1(0),
952 fHistoEleLambdaPtvsV0dlWSSide2(0),
953 fHistoEleLambdaPtvsV0dlWSMix(0),
954 fHistoEleLambdaPtvsV0dlWSMix1(0),
955 fHistoEleLambdaPtvsV0dlWSMix2(0),
956 fHistoEleLambdaPtvsV0dlMCS(0),
957 fHistoEleLambdaPtvsV0dlMCS1(0),
958 fHistoEleLambdaPtvsV0dlMCS2(0),
959 fHistoEleLambdaPtvsV0dlPromptMCS(0),
960 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
961 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
962 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
963 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
964 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
965 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
966 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
967 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
968 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
969 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
970 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
971 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
972 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
973 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
974 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
975 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
976 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
977 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
978 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
979 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
980 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
981 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
982 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
983 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
984 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
985 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
986 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
987 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
988 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
989 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
990 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
991 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
992 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
993 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
994 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
995 fHistoEleLambdaPtvsV0dlXibMCS(0),
996 fHistoEleLambdaPtvsV0dlXibMCS1(0),
997 fHistoEleLambdaPtvsV0dlXibMCS2(0),
998 fHistoEleLambdaPtvsV0dlRSAway(0),
999 fHistoEleLambdaPtvsV0dlRS1Away(0),
1000 fHistoEleLambdaPtvsV0dlRS2Away(0),
1001 fHistoEleLambdaPtvsV0dlRSSideAway(0),
1002 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
1003 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
1004 fHistoEleLambdaPtvsV0dlRSMixAway(0),
1005 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
1006 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
1007 fHistoEleLambdaPtvsV0dlWSAway(0),
1008 fHistoEleLambdaPtvsV0dlWS1Away(0),
1009 fHistoEleLambdaPtvsV0dlWS2Away(0),
1010 fHistoEleLambdaPtvsV0dlWSSideAway(0),
1011 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
1012 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
1013 fHistoEleLambdaPtvsV0dlWSMixAway(0),
1014 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
1015 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
1016 fHistoCorrelationVariablesvsEleLambdaPt(0),
1017 fHistoCorrelationVariablesvsEleLambdaPtMix(0),
1018 fHistoCorrelationVariablesvsEleLambdaPtMC(0),
1019 fHistoCorrelationVariablesvsElePt(0),
1020 fHistoCorrelationVariablesvsElePtMix(0),
1021 fHistoCorrelationVariablesvsElePtMC(0),
1022 fHistoCorrelationVariablesvsLambdaPt(0),
1023 fHistoCorrelationVariablesvsLambdaPtMix(0),
1024 fHistoCorrelationVariablesvsLambdaPtMC(0),
1025 fHistoMassVariablesvsEleLambdaPt(0),
1026 fHistoMassVariablesvsEleLambdaPtMix(0),
1027 fHistoMassVariablesvsEleLambdaPtMC(0),
1028 fHistoMassVariablesvsElePt(0),
1029 fHistoMassVariablesvsElePtMix(0),
1030 fHistoMassVariablesvsElePtMC(0),
1031 fHistoMassVariablesvsLambdaPt(0),
1032 fHistoMassVariablesvsLambdaPtMix(0),
1033 fHistoMassVariablesvsLambdaPtMC(0),
1034 fHistoResponseElePt(0),
1035 fHistoResponseElePt1(0),
1036 fHistoResponseElePt2(0),
1037 fHistoResponseEleLambdaPt(0),
1038 fHistoResponseEleLambdaPt1(0),
1039 fHistoResponseEleLambdaPt2(0),
1040 fHistoResponseEleLambdaPtFeeddownXic0(0),
1041 fHistoResponseEleLambdaPtFeeddownXic01(0),
1042 fHistoResponseEleLambdaPtFeeddownXic02(0),
1043 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
1044 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
1045 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
1046 fHistoResponseEleLambdaPtFeeddownSigma(0),
1047 fHistoResponseEleLambdaPtFeeddownSigma1(0),
1048 fHistoResponseEleLambdaPtFeeddownSigma2(0),
1049 fHistoResponseEleLambdaPtLambdab(0),
1050 fHistoResponseElePtLambdab(0),
1051 fHistoResponseLambdaPtLambdab(0),
1052 fHistoResponseLambdaPt(0),
1053 fHistoResponseLambdaPtFeeddownXic0(0),
1054 fHistoResponseLambdaPtFeeddownXicPlus(0),
1055 fHistoResponseLambdaPtFeeddownSigma(0),
1056 fHistoResponseLambdaPtvsEleLambdaPt(0),
1057 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
1058 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
1059 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
1060 fHistoEleLambdaPtvsRapidityRS(0),
1061 fHistoEleLambdaPtvsRapidityWS(0),
1062 fHistoEleLambdaPtvsRapidityMCS(0),
1064 fHistoElectronPi0Total(0),
1065 fHistoElectronPi0Tag(0),
1067 fHistoElectronEtaTotal(0),
1068 fHistoElectronEtaTag(0),
1070 fHistonEvtvsRunNumber(0),
1071 fHistonElevsRunNumber(0),
1072 fHistonLambdavsRunNumber(0),
1073 fHistoMCEventType(0),
1074 fHistoMCDeltaPhiccbar(0),
1076 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
1077 fHistodPhiSdEtaSElectronProtonR125RS(0),
1078 fHistodPhiSdEtaSElectronProtonR125WS(0),
1079 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
1080 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
1081 fHistodPhiSdEtaSElectronPionR125RS(0),
1082 fHistodPhiSdEtaSElectronPionR125WS(0),
1083 fHistodPhiSdEtaSElectronPionR125RSMix(0),
1084 fHistodPhiSdEtaSElectronPionR125WSMix(0),
1086 fMixWithoutConversionFlag(kFALSE),
1087 fNumberOfEventsForMixing (5),
1099 m_ReservoirVarsL1(),
1105 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
1107 for(
Int_t i=0;i<17;i++){
1112 for(
Int_t i=0;i<8;i++){
1117 DefineOutput(1,TList::Class());
1118 DefineOutput(2,TList::Class());
1119 DefineOutput(3,TList::Class());
1120 DefineOutput(4,TTree::Class());
1121 DefineOutput(5,TTree::Class());
1122 DefineOutput(6,TTree::Class());
1123 DefineOutput(7,TTree::Class());
1124 DefineOutput(8,AliNormalizationCounter::Class());
1125 DefineOutput(9,TTree::Class());
1126 DefineOutput(10,TTree::Class());
1127 DefineOutput(11,TTree::Class());
1128 DefineOutput(12,TTree::Class());
1136 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1246 if (fDebug > 1) AliInfo(
"Init");
1265 AliError(
"NO EVENT FOUND!");
1277 AliKFParticle::SetField(
fBzkG);
1278 if (TMath::Abs(
fBzkG)<0.001) {
1287 AliAODTracklets* tracklets=aodEvent->GetTracklets();
1288 Int_t nTr=tracklets->GetNumberOfTracklets();
1289 for(
Int_t iTr=0; iTr<nTr; iTr++){
1290 Double_t theta=tracklets->GetTheta(iTr);
1291 Double_t eta=-TMath::Log(TMath::Tan(theta/2.));
1292 if(eta>-1.0 && eta<1.0) countTr++;
1294 AliAODVertex *vtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1297 if(vtx1->GetNContributors()>0){
1321 TClonesArray *mcArray = 0;
1322 AliAODMCHeader *mcHeader=0;
1325 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1327 AliError(
"Could not find Monte-Carlo in AOD");
1333 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1335 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1340 Double_t zMCVertex = mcHeader->GetVtxZ();
1342 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1356 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1361 fVtx1->GetCovarianceMatrix(cov);
1362 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1366 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1376 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1377 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1378 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1379 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1380 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1394 AliCentrality *cent = aodEvent->GetCentrality();
1395 fCentrality = cent->GetCentralityPercentile(
"V0M");
1397 if(countCorr>=0 && countCorr<=0)
fCentrality = 5.;
1398 else if(countCorr>=1 && countCorr<=8)
fCentrality = 15.;
1399 else if(countCorr>= 9 && countCorr<=13)
fCentrality = 25.;
1400 else if(countCorr>=14 && countCorr<=19)
fCentrality = 35.;
1401 else if(countCorr>=20 && countCorr<=30)
fCentrality = 45.;
1402 else if(countCorr>=31 && countCorr<=49)
fCentrality = 55.;
1407 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1414 AliEventplane *pl=aodEvent->GetEventplane();
1416 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec:no eventplane! v2 analysis without eventplane not possible!\n");
1434 fQ = pl->GetQVector();
1436 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: no Q vectors");
1446 Int_t runnumber_offset = 0;
1447 Int_t runnumber = aodEvent->GetRunNumber();
1448 if(runnumber<=131000&&runnumber>=114000){
1449 runnumber_offset = 114000;
1450 }
else if(runnumber<=196000&&runnumber>=195000){
1451 runnumber_offset = 195000;
1452 }
else if(runnumber<=170593&&runnumber>=167902){
1453 runnumber_offset = 167902;
1499 AliAnalysisTaskSE::Terminate();
1504 AliError(
"fOutput not available");
1510 AliError(
"fOutputAll not available");
1567 TString normName=
"NormalizationCounter";
1568 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1569 if(cont)normName=(
TString)cont->GetName();
1616 fPoolIndex=GetPoolIndex(fVtxZ,fCentrality,fEventPlane);
1617 Int_t nextRes( nextResVec[fPoolIndex] );
1618 while(!m_ReservoirE[fPoolIndex][nextRes].empty()){
1619 delete m_ReservoirE[fPoolIndex][nextRes].back();
1620 m_ReservoirE[fPoolIndex][nextRes].pop_back();
1622 while(!m_ReservoirL1[fPoolIndex][nextRes].empty()){
1623 delete m_ReservoirL1[fPoolIndex][nextRes].back();
1624 m_ReservoirL1[fPoolIndex][nextRes].pop_back();
1626 while(!m_ReservoirL2[fPoolIndex][nextRes].empty()){
1627 delete m_ReservoirL2[fPoolIndex][nextRes].back();
1628 m_ReservoirL2[fPoolIndex][nextRes].pop_back();
1630 while(!m_ReservoirVarsE[fPoolIndex][nextRes].empty()){
1631 delete m_ReservoirVarsE[fPoolIndex][nextRes].back();
1632 m_ReservoirVarsE[fPoolIndex][nextRes].pop_back();
1634 while(!m_ReservoirVarsL1[fPoolIndex][nextRes].empty()){
1635 delete m_ReservoirVarsL1[fPoolIndex][nextRes].back();
1636 m_ReservoirVarsL1[fPoolIndex][nextRes].pop_back();
1638 while(!m_ReservoirVarsL2[fPoolIndex][nextRes].empty()){
1639 delete m_ReservoirVarsL2[fPoolIndex][nextRes].back();
1640 m_ReservoirVarsL2[fPoolIndex][nextRes].pop_back();
1645 if(fAnalCuts->GetProdAODFilterBit()==7){
1646 ResetGlobalTrackReference();
1648 for (
Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1650 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1651 if (!track)
continue;
1654 StoreGlobalTrackReference(track,iTrack);
1662 Int_t nV0s= aodEvent->GetNumberOfV0s();
1663 Int_t nTracks= aodEvent->GetNumberOfTracks();
1665 Bool_t seleTrkFlags[nTracks];
1667 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1669 Bool_t seleV0Flags[nV0s];
1671 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1673 Int_t runnumber_offset = 0;
1674 Int_t runnumber = aodEvent->GetRunNumber();
1675 if(runnumber<=131000&&runnumber>=114000){
1676 runnumber_offset = 114000;
1677 }
else if(runnumber<=196000&&runnumber>=195000){
1678 runnumber_offset = 195000;
1679 }
else if(runnumber<=170593&&runnumber>=167902){
1680 runnumber_offset = 167902;
1682 fHistonElevsRunNumber->Fill(runnumber-runnumber_offset,nSeleTrks);
1683 fHistonLambdavsRunNumber->Fill(runnumber-runnumber_offset,nSeleV0);
1685 if(nSeleTrks==0 || nSeleV0 ==0)
return;
1690 for (
Int_t itrk = 0; itrk<nTracks; itrk++) {
1691 if(!seleTrkFlags[itrk])
continue;
1692 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1694 FillElectronROOTObjects(trk,aodEvent,mcArray);
1697 for (
Int_t iv0 = 0; iv0<nV0s; iv0++) {
1698 if(!seleV0Flags[iv0])
continue;
1699 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1702 FillV0ROOTObjects(v0, aodEvent, mcArray);
1705 if(fWriteEachVariableTree)
1711 for (
Int_t iv0 = 0; iv0<nV0s; iv0++) {
1712 if(!seleV0Flags[iv0])
continue;
1713 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1716 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1717 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1722 for (
Int_t itrk = 0; itrk<nTracks; itrk++) {
1723 if(!seleTrkFlags[itrk])
continue;
1724 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1732 AliAODTrack *trkpid = 0;
1733 if(fAnalCuts->GetProdAODFilterBit()==7){
1734 trkpid = fGTI[-trk->GetID()-1];
1739 Int_t cpid = cptrack->GetID();
1740 Int_t cnid = cntrack->GetID();
1741 Int_t lpid = trkpid->GetID();
1742 if((cpid==lpid)||(cnid==lpid))
continue;
1744 AliAODVertex *secVert = ReconstructSecondaryVertex(v0,trk,aodEvent);
1745 if(!secVert)
continue;
1752 FillROOTObjects(elobj,v0,trk,trkpid,aodEvent,mcArray);
1754 elobj->GetSecondaryVtx()->RemoveDaughters();
1756 delete elobj;elobj=NULL;
1762 DoEventMixingWithPools(fPoolIndex);
1764 Int_t nextRes( nextResVec[fPoolIndex] );
1766 if( nextRes>=fNumberOfEventsForMixing ){
1768 reservoirsReady[fPoolIndex] = kTRUE;
1770 nextResVec[fPoolIndex] = nextRes;
1780 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1786 fCandidateVariableNames[ 0]=
"Centrality";
1787 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1788 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1789 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1790 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1791 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1792 fCandidateVariableNames[ 6]=
"ElePx";
1793 fCandidateVariableNames[ 7]=
"ElePy";
1794 fCandidateVariableNames[ 8]=
"ElePz";
1795 fCandidateVariableNames[ 9]=
"V0Px";
1796 fCandidateVariableNames[10]=
"V0Py";
1797 fCandidateVariableNames[11]=
"V0Pz";
1798 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1799 fCandidateVariableNames[13]=
"MassLambda";
1800 fCandidateVariableNames[14]=
"MassAntiLambda";
1801 fCandidateVariableNames[15]=
"Eled0";
1802 fCandidateVariableNames[16]=
"V0d0";
1803 fCandidateVariableNames[17]=
"nSigmaTPCele";
1804 fCandidateVariableNames[18]=
"nSigmaTOFele";
1805 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1806 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1807 fCandidateVariableNames[21]=
"EleCharge";
1808 fCandidateVariableNames[22]=
"ProtonPx";
1809 fCandidateVariableNames[23]=
"ProtonPy";
1810 fCandidateVariableNames[24]=
"ProtonPz";
1811 fCandidateVariableNames[25]=
"PiPx";
1812 fCandidateVariableNames[26]=
"PiPy";
1813 fCandidateVariableNames[27]=
"PiPz";
1814 fCandidateVariableNames[28]=
"mcpdglc";
1815 fCandidateVariableNames[29]=
"mclablc";
1816 fCandidateVariableNames[30]=
"mcpdgmomele";
1817 fCandidateVariableNames[31]=
"mcpdgmomv0";
1818 fCandidateVariableNames[32]=
"Mixing";
1819 fCandidateVariableNames[33]=
"mcpdgele";
1820 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1821 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1822 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1823 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1824 fCandidateVariableNames[38]=
"MassK0Short";
1825 fCandidateVariableNames[39]=
"mcpdggrmomele";
1826 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1827 fCandidateVariableNames[41]=
"mcngenele";
1828 fCandidateVariableNames[42]=
"mcngenv0";
1829 fCandidateVariableNames[43]=
"mclcpx";
1830 fCandidateVariableNames[44]=
"mclcpy";
1831 fCandidateVariableNames[45]=
"mclcpz";
1832 fCandidateVariableNames[46]=
"mcelepx";
1833 fCandidateVariableNames[47]=
"mcelepy";
1834 fCandidateVariableNames[48]=
"mcelepz";
1835 fCandidateVariableNames[49]=
"mcv0px";
1836 fCandidateVariableNames[50]=
"mcv0py";
1837 fCandidateVariableNames[51]=
"mcv0pz";
1838 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1839 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1840 fCandidateVariableNames[54]=
"PrimVertx";
1841 fCandidateVariableNames[55]=
"PrimVerty";
1842 fCandidateVariableNames[56]=
"PrimVertz";
1843 fCandidateVariableNames[57]=
"V0Vertx";
1844 fCandidateVariableNames[58]=
"V0Verty";
1845 fCandidateVariableNames[59]=
"V0Vertz";
1847 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1848 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1849 fCandidateVariableNames[62]=
"DcaV0daughters";
1850 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1851 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1852 fCandidateVariableNames[65]=
"MassK0Short2";
1854 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1855 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1857 fCandidateVariableNames[68]=
"EleITSMatch";
1858 fCandidateVariableNames[69]=
"V0PosITSMatch";
1859 fCandidateVariableNames[70]=
"V0NegITSMatch";
1860 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1861 fCandidateVariableNames[72]=
"mcpdgv0";
1863 fCandidateVariableNames[73]=
"mcv0pospx";
1864 fCandidateVariableNames[74]=
"mcv0pospy";
1865 fCandidateVariableNames[75]=
"mcv0pospz";
1866 fCandidateVariableNames[76]=
"mcpdgv0pos";
1867 fCandidateVariableNames[77]=
"mcv0negpx";
1868 fCandidateVariableNames[78]=
"mcv0negpy";
1869 fCandidateVariableNames[79]=
"mcv0negpz";
1870 fCandidateVariableNames[80]=
"mcpdgv0neg";
1872 fCandidateVariableNames[81]=
"ProtonPxPV";
1873 fCandidateVariableNames[82]=
"ProtonPyPV";
1874 fCandidateVariableNames[83]=
"ProtonPzPV";
1875 fCandidateVariableNames[84]=
"PiPxPV";
1876 fCandidateVariableNames[85]=
"PiPyPV";
1877 fCandidateVariableNames[86]=
"PiPzPV";
1879 fCandidateVariableNames[87]=
"mcv0vertx";
1880 fCandidateVariableNames[88]=
"mcv0verty";
1881 fCandidateVariableNames[89]=
"mcv0vertz";
1883 fCandidateVariableNames[90]=
"OnFlyStatus";
1884 fCandidateVariableNames[91]=
"MagneticField";
1885 fCandidateVariableNames[92]=
"EvNumber";
1886 fCandidateVariableNames[93]=
"RunNumber";
1906 for(
Int_t i=0;i<94;i++)
1911 Bool_t anti_lambda_flag = kFALSE;
1914 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1915 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1924 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1925 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1929 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1930 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1931 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
1937 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
1938 Double_t pxv_flip = v0px + lf * uxe;
1939 Double_t pyv_flip = v0py + lf * uye;
1940 Double_t pzv_flip = v0pz + lf * uze;
1941 Double_t pxsum_flip = epx + pxv_flip;
1942 Double_t pysum_flip = epy + pyv_flip;
1943 Double_t pzsum_flip = epz + pzv_flip;
1944 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
1945 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
1953 UInt_t pdgdg[2]={11,3122};
1987 if(anti_lambda_flag){
2005 if(anti_lambda_flag){
2021 AliAODMCParticle *mclc = 0;
2022 AliAODMCParticle *mcele = 0;
2023 AliAODMCParticle *mcv0 = 0;
2025 Int_t mcpdgele_array[100];
2026 Int_t mcpdgv0_array[100];
2027 Int_t mclabelele_array[100];
2028 Int_t mclabelv0_array[100];
2029 Int_t mcngen_ele=-9999;
2030 Int_t mcngen_v0=-9999;
2033 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
2034 if(mclabelele_array[0]>=0){
2035 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
2042 if(mclabelv0_array[0]>=0){
2043 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
2051 Int_t pdgdgv0[2]={211,211};
2052 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
2054 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
2064 Int_t labPos = cptrack->GetLabel();
2066 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
2075 Int_t labNeg = cntrack->GetLabel();
2077 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
2085 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
2128 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
2129 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
2130 if(!anti_lambda_flag){
2139 fVtx1->GetXYZ(posVtx);
2141 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
2151 if(!anti_lambda_flag){
2176 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
2177 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
2181 cont[0] = elobj->InvMass(2,pdgdg);
2182 cont[1] = elobj->Pt();
2187 cont_flip[0] = melam_flip;
2188 cont_flip[1] = ptelam_flip;
2192 cont2[0] = elobj->InvMass(2,pdgdg);
2193 cont2[1] = trk->Pt();
2196 cont_eleptvseta[0] = trk->Pt();
2197 cont_eleptvseta[1] = trk->Eta();
2201 cont_eleptvslambdapt[0] = trk->Pt();
2202 cont_eleptvslambdapt[1] = v0->Pt();
2206 cont_eleptvsd0[0] = trk->Pt();
2207 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
2211 cont_eleptvsv0dl[0] = trk->Pt();
2212 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2216 cont_elelamptvsv0dl[0] = elobj->Pt();
2217 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
2218 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
2221 Double_t cont_elelamptvsv0dl_flip[4];
2222 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
2223 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2224 cont_elelamptvsv0dl_flip[2] = 0.0;
2229 cont_eleptvsv0dca[0] = trk->Pt();
2230 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
2233 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
2237 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
2238 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
2245 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2272 for(
Int_t ih=0;ih<17;ih++){
2274 cont_eleptvscutvars[0] = elobj->Pt();
2278 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2280 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2282 cont_eleptvscutvars[1] = nSigmaTPCele;
2284 cont_eleptvscutvars[1] = nSigmaTOFele;
2286 cont_eleptvscutvars[1] = trk->Eta();
2288 cont_eleptvscutvars[1] = trk->GetITSNcls();
2290 if(!anti_lambda_flag)
2291 cont_eleptvscutvars[1] = v0->MassLambda();
2293 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2296 lPosV0[0] = v0->DecayVertexV0X();
2297 lPosV0[1] = v0->DecayVertexV0Y();
2298 lPosV0[2] = v0->DecayVertexV0Z();
2299 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2301 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2303 if(!anti_lambda_flag)
2304 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2306 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2308 if(!anti_lambda_flag)
2309 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2311 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2313 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2315 cont_eleptvscutvars[1] = v0->MassK0Short();
2317 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2319 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2321 cont_eleptvscutvars[1] = v0->Eta();
2329 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2331 cont_eleptvscutvars[1] = -9999.;
2339 Int_t pdgcode = mclc->GetPdgCode();
2341 cont_mclc[0] = mclc->Pt();
2342 cont_mclc[1] = mclc->Y();
2345 cont_mcele[0] = mcele->Pt();
2346 cont_mcele[1] = mcele->Eta();
2349 Int_t labmotherlc = mclc->GetMother();
2350 Bool_t isbottomfd = kFALSE;
2352 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2353 Int_t pdgmotherlc = motherlc->GetPdgCode();
2354 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){
2358 cont2[1] = mcele->Pt();
2360 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2412 Double_t cont_eleptvslambdaptvslcpt[4];
2413 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2414 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2415 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2416 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2420 cont_allpt[0] = mclc->Pt();
2421 cont_allpt[1] = elobj->Pt();
2422 cont_allpt[2] = trk->Pt();
2423 cont_allpt[3] = v0->Pt();
2432 for(
Int_t ih=0;ih<17;ih++){
2434 cont_eleptvscutvars[0] = elobj->Pt();
2438 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2440 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2442 cont_eleptvscutvars[1] = nSigmaTPCele;
2444 cont_eleptvscutvars[1] = nSigmaTOFele;
2446 cont_eleptvscutvars[1] = trk->Eta();
2448 cont_eleptvscutvars[1] = trk->GetITSNcls();
2450 if(!anti_lambda_flag)
2451 cont_eleptvscutvars[1] = v0->MassLambda();
2453 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2456 lPosV0[0] = v0->DecayVertexV0X();
2457 lPosV0[1] = v0->DecayVertexV0Y();
2458 lPosV0[2] = v0->DecayVertexV0Z();
2459 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2461 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2463 if(!anti_lambda_flag)
2464 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2466 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2468 if(!anti_lambda_flag)
2469 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2471 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2473 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2475 cont_eleptvscutvars[1] = v0->MassK0Short();
2477 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2479 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2481 cont_eleptvscutvars[1] = v0->Eta();
2489 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2491 cont_eleptvscutvars[1] = -9999.;
2498 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2545 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2592 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2622 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2649 for(
Int_t ih=0;ih<17;ih++){
2651 cont_eleptvscutvars[0] = elobj->Pt();
2655 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2657 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2659 cont_eleptvscutvars[1] = nSigmaTPCele;
2661 cont_eleptvscutvars[1] = nSigmaTOFele;
2663 cont_eleptvscutvars[1] = trk->Eta();
2665 cont_eleptvscutvars[1] = trk->GetITSNcls();
2667 if(!anti_lambda_flag)
2668 cont_eleptvscutvars[1] = v0->MassLambda();
2670 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2673 lPosV0[0] = v0->DecayVertexV0X();
2674 lPosV0[1] = v0->DecayVertexV0Y();
2675 lPosV0[2] = v0->DecayVertexV0Z();
2676 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2678 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2680 if(!anti_lambda_flag)
2681 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2683 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2685 if(!anti_lambda_flag)
2686 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2688 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2690 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2692 cont_eleptvscutvars[1] = v0->MassK0Short();
2694 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2696 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2698 cont_eleptvscutvars[1] = v0->Eta();
2706 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2708 cont_eleptvscutvars[1] = -9999.;
2716 Int_t pdgcode = mclc->GetPdgCode();
2718 cont_mclb[0] = mclc->Pt();
2719 cont_mclb[1] = mclc->Y();
2722 cont_mcele[0] = mcele->Pt();
2723 cont_mcele[1] = mcele->Eta();
2726 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2747 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2759 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2771 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2783 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2803 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2805 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2809 if(cont_flip[0]<2.3){
2818 if(cont_flip[0]<2.3){
2828 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2849 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2873 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2875 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2879 if(cont_flip[0]<2.3){
2888 if(cont_flip[0]<2.3){
2899 for(
Int_t iv=0;iv<14;iv++){
2904 for(
Int_t iv=0;iv<7;iv++){
2905 cont_cor_nd[iv] = -9999.;
2908 for(
Int_t iv=0;iv<8;iv++){
2909 cont_mass_nd[iv] = -9999.;
2920 if(trk->Charge()>0){
2923 }
else if(trk->Charge()<0){
2933 cont_cor_nd[0] = elobj->Pt();
2935 cont_cor_nd[2] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;;
2936 if(trk->Charge()>0){
2937 if(!anti_lambda_flag) cont_cor_nd[3] = 0;
2938 else cont_cor_nd[3] = 2;
2939 }
else if(trk->Charge()<0){
2940 if(!anti_lambda_flag) cont_cor_nd[3] = 3;
2941 else cont_cor_nd[3] = 1;
2949 Int_t pdgcode = mclc->GetPdgCode();
2950 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2954 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2958 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2962 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2966 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2970 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2974 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2978 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2982 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
3020 if(elobj->InvMass(2,pdgdg)<2.3)
3022 cont_cor_nd[0] = trk->Pt();
3024 cont_cor_nd[0] = v0->Pt();
3027 if(elobj->InvMass(2,pdgdg)<2.3)
3029 cont_cor_nd[0] = trk->Pt();
3031 cont_cor_nd[0] = v0->Pt();
3036 cont_mass_nd[0] = elobj->InvMass(2,pdgdg);
3037 cont_mass_nd[1] = cont_cor_nd[0];
3038 cont_mass_nd[4] = cont_cor_nd[3];
3039 cont_mass_nd[5] = cont_cor_nd[4];
3040 cont_mass_nd[6] = cont_cor_nd[5];
3041 cont_mass_nd[7] = cont_cor_nd[6];
3046 if(melam_flip < 10.&& cosoa < 0.) cont_mass_nd[2]=0;
3049 cont_mass_nd[1] = trk->Pt();
3051 cont_mass_nd[1] = v0->Pt();
3055 cont_mass_nd[1] = trk->Pt();
3057 cont_mass_nd[1] = v0->Pt();
3073 for(
Int_t i=0;i<94;i++)
3082 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
3083 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
3088 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
3089 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
3091 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
3097 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
3102 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
3103 Double_t pxv_flip = pxv + lf * uxe;
3104 Double_t pyv_flip = pyv + lf * uye;
3105 Double_t pzv_flip = pzv + lf * uze;
3106 Double_t pxsum_flip = pxe + pxv_flip;
3107 Double_t pysum_flip = pye + pyv_flip;
3108 Double_t pzsum_flip = pze + pzv_flip;
3109 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
3110 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
3113 UInt_t pdgdg[2]={11,3122};
3166 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
3171 cont_flip[0] = mel_flip;
3172 cont_flip[1] = ptel_flip;
3177 cont2[1] = sqrt(pxe*pxe+pye*pye);
3180 cont_eleptvseta[0] = trke->Pt();
3181 cont_eleptvseta[1] = trke->Eta();
3185 cont_eleptvslambdapt[0] = trke->Pt();
3186 cont_eleptvslambdapt[1] = v0->Pt();
3190 cont_eleptvsd0[0] = trke->Pt();
3191 cont_eleptvsd0[1] = 0.;
3195 cont_eleptvsv0dl[0] = trke->Pt();
3196 cont_eleptvsv0dl[1] = (*v0vars)[15];
3200 cont_eleptvsv0dca[0] = trke->Pt();
3201 cont_eleptvsv0dca[1] = (*v0vars)[16];
3205 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
3206 cont_elelamptvsv0dl[1] = (*v0vars)[15];
3207 cont_elelamptvsv0dl[2] = 0.0;
3210 Double_t cont_elelamptvsv0dl_flip[4];
3211 cont_elelamptvsv0dl_flip[0] = ptel_flip;
3212 cont_elelamptvsv0dl_flip[1] = (*v0vars)[15];
3213 cont_elelamptvsv0dl_flip[2] = 0.0;
3216 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
3217 xyzR125_ele[0] = (*elevars)[0];
3218 xyzR125_ele[1] = (*elevars)[1];
3219 xyzR125_ele[2] = (*elevars)[2];
3220 xyzR125_pr[0] = (*v0vars)[0];
3221 xyzR125_pr[1] = (*v0vars)[1];
3222 xyzR125_pr[2] = (*v0vars)[2];
3223 xyzR125_pi[0] = (*v0vars)[3];
3224 xyzR125_pi[1] = (*v0vars)[4];
3225 xyzR125_pi[2] = (*v0vars)[5];
3228 rdhfcutvars[0] = xyzR125_ele[0];
3229 rdhfcutvars[1] = xyzR125_ele[1];
3230 rdhfcutvars[2] = xyzR125_ele[2];
3231 rdhfcutvars[3] = xyzR125_pr[0];
3232 rdhfcutvars[4] = xyzR125_pr[1];
3233 rdhfcutvars[5] = xyzR125_pr[2];
3234 rdhfcutvars[6] = xyzR125_pi[0];
3235 rdhfcutvars[7] = xyzR125_pi[1];
3236 rdhfcutvars[8] = xyzR125_pi[2];
3247 if(((
int)trke->T())*chargepr>0){
3307 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
3309 if(((
int)trke->T())*chargepr>0){
3313 if(cont_flip[0]<2.3){
3322 if(cont_flip[0]<2.3){
3333 for(
Int_t iv=0;iv<14;iv++){
3338 for(
Int_t iv=0;iv<7;iv++){
3339 cont_cor_nd[iv] = -9999.;
3342 for(
Int_t iv=0;iv<8;iv++){
3343 cont_mass_nd[iv] = -9999.;
3350 fCorrelationVariables[4] = (*v0vars)[15];
3351 fCorrelationVariables[5] = (*elevars)[5];
3352 fCorrelationVariables[6] = 2;
3354 if(chargepr>0) fCorrelationVariables[7] = 0;
3355 else fCorrelationVariables[7] = 2;
3356 }
else if(trke->T()<0){
3357 if(chargepr>0) fCorrelationVariables[7] = 1;
3358 else fCorrelationVariables[7] = 3;
3360 fCorrelationVariables[8] = (*elevars)[6];
3361 fCorrelationVariables[9] = (*elevars)[7];
3363 fCorrelationVariables[11] = sqrt(pxsum*pxsum+pysum*pysum);
3364 fCorrelationVariables[12] = mel;
3365 fCorrelationVariables[13] = cosoa;
3367 cont_cor_nd[0] = sqrt(pxsum*pxsum+pysum*pysum);
3368 cont_cor_nd[1] = TVector2::Phi_mpi_pi(v0->Phi()-trke->Phi());
3369 cont_cor_nd[2] = (*v0vars)[15];
3371 if(chargepr>0) cont_cor_nd[3] = 0;
3372 else cont_cor_nd[3] = 2;
3373 }
else if(trke->T()<0){
3374 if(chargepr>0) cont_cor_nd[3] = 3;
3375 else cont_cor_nd[3] = 1;
3377 cont_cor_nd[4] = fCorrelationVariables[8];
3379 if(fabs(fCorrelationVariables[9]-1013)<0.001) cont_cor_nd[5] = 7;
3380 if(fabs(fCorrelationVariables[9]-1016)<0.001) cont_cor_nd[5] = 8;
3381 if(fabs(fCorrelationVariables[9]-1019)<0.001) cont_cor_nd[5] = 9;
3391 cont_cor_nd[0] = trke->Pt();
3393 cont_cor_nd[0] = v0->Pt();
3397 cont_mass_nd[0] = mel;
3398 cont_mass_nd[1] = cont_cor_nd[0];
3399 cont_mass_nd[4] = cont_cor_nd[3];
3400 cont_mass_nd[5] = cont_cor_nd[4];
3401 cont_mass_nd[6] = cont_cor_nd[5];
3402 cont_mass_nd[7] = cont_cor_nd[6];
3407 if(mel_flip < 10.&& cosoa < 0.) cont_mass_nd[2]=0;
3409 cont_mass_nd[0] = trke->Pt();
3411 cont_mass_nd[0] = v0->Pt();
3423 const char* nameoutput = GetOutputSlot(11)->GetContainer()->GetName();
3429 fCandidateVariableNames[ 0]=
"Px";
3430 fCandidateVariableNames[ 1]=
"Py";
3431 fCandidateVariableNames[ 2]=
"Pz";
3432 fCandidateVariableNames[ 3]=
"Charge";
3433 fCandidateVariableNames[ 4]=
"LambdaMass";
3434 fCandidateVariableNames[ 5]=
"Bz";
3435 fCandidateVariableNames[ 6]=
"Centrality";
3436 fCandidateVariableNames[ 7]=
"PrimVertZ";
3437 fCandidateVariableNames[ 8]=
"EvNumber";
3438 fCandidateVariableNames[ 9]=
"RunNumber";
3439 fCandidateVariableNames[10]=
"EventPlane";
3440 fCandidateVariableNames[11]=
"Vars0";
3441 fCandidateVariableNames[12]=
"Vars1";
3442 fCandidateVariableNames[13]=
"Vars2";
3443 fCandidateVariableNames[14]=
"Vars3";
3444 fCandidateVariableNames[15]=
"Vars4";
3463 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
3469 fCandidateVariableNames[ 0]=
"ElePx";
3470 fCandidateVariableNames[ 1]=
"ElePy";
3471 fCandidateVariableNames[ 2]=
"ElePz";
3472 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
3473 fCandidateVariableNames[ 4]=
"ITSNcls";
3474 fCandidateVariableNames[ 5]=
"TPCNcls";
3475 fCandidateVariableNames[ 6]=
"TPCNclsPID";
3476 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
3477 fCandidateVariableNames[ 8]=
"d0R";
3478 fCandidateVariableNames[ 9]=
"d0Z";
3479 fCandidateVariableNames[10]=
"ITSClusterMap";
3480 fCandidateVariableNames[11]=
"nSigmaTPCele";
3481 fCandidateVariableNames[12]=
"nSigmaTOFele";
3482 fCandidateVariableNames[13]=
"nSigmaTPCpi";
3483 fCandidateVariableNames[14]=
"nSigmaTPCka";
3484 fCandidateVariableNames[15]=
"nSigmaTPCpr";
3485 fCandidateVariableNames[16]=
"EvNumber";
3486 fCandidateVariableNames[17]=
"EleCharge";
3487 fCandidateVariableNames[18]=
"ElePdgCode";
3488 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3489 fCandidateVariableNames[20]=
"mcelepx";
3490 fCandidateVariableNames[21]=
"mcelepy";
3491 fCandidateVariableNames[22]=
"mcelepz";
3492 fCandidateVariableNames[23]=
"Centrality";
3493 fCandidateVariableNames[24]=
"PrimVertZ";
3494 fCandidateVariableNames[25]=
"RunNumber";
3511 AliAODTrack *trkpid = 0;
3513 trkpid =
fGTI[-trk->GetID()-1];
3522 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3531 Int_t pdgEle = -9999;
3532 Int_t pdgEleMother = -9999;
3536 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
3537 for(
Int_t i=0;i<100;i++){
3538 pdgarray_ele[i]=-9999;
3539 labelarray_ele[i]=-9999;
3542 Int_t mcetype = -9999;
3545 Int_t labEle = trk->GetLabel();
3546 if(labEle<0)
return;
3547 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3549 if(abs(mcetrk->GetPdgCode())!=11)
return;
3551 pdgEle = mcetrk->GetPdgCode();
3552 mcelepx = mcetrk->Px();
3553 mcelepy = mcetrk->Py();
3554 mcelepz = mcetrk->Pz();
3556 Int_t labEleMother = mcetrk->GetMother();
3557 if(labEleMother>-1){
3558 AliAODMCParticle *mcemothertrk = (AliAODMCParticle*)mcArray->At(labEleMother);
3560 pdgEleMother = mcemothertrk->GetPdgCode();
3567 Bool_t gamma_flag = kFALSE;
3568 Bool_t pi0_flag = kFALSE;
3569 Bool_t eta_flag = kFALSE;
3572 if(abs(pdgarray_ele[0])>400&&abs(pdgarray_ele[0])<440){
3575 if(abs(pdgarray_ele[0])>4000&&abs(pdgarray_ele[0])<4400){
3578 if(abs(pdgarray_ele[0])>500&&abs(pdgarray_ele[0])<540){
3581 if(abs(pdgarray_ele[0])>5000&&abs(pdgarray_ele[0])<5400){
3585 if(abs(pdgarray_ele[0])==22){
3592 if((abs(pdgarray_ele[0])==22) && (abs(pdgarray_ele[1])==111)){
3594 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[1]);
3595 pt_pi0 = mctrkm->Pt();
3597 if(abs(pdgarray_ele[0])==111){
3599 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[0]);
3600 pt_pi0 = mctrkm->Pt();
3602 if((abs(pdgarray_ele[0])==22) && (abs(pdgarray_ele[1])==221)){
3604 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[1]);
3605 pt_eta = mctrkm->Pt();
3607 if(abs(pdgarray_ele[0])==221){
3609 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[0]);
3610 pt_eta = mctrkm->Pt();
3615 cont_pi0[0] = trk->Pt();
3616 cont_pi0[1] = pt_pi0;
3624 cont_eta[0] = trk->Pt();
3625 cont_eta[1] = pt_eta;
3631 if(hfe_flag==0)
return;
3646 pv[0] =
fVtx1->GetX();
3647 pv[1] =
fVtx1->GetY();
3648 pv[2] =
fVtx1->GetZ();
3649 Double_t xyzR125[3]={9999.,9999.,9999.};
3651 TVector *varvec =
new TVector(8);
3652 (*varvec)[0] = xyzR125[0];
3653 (*varvec)[1] = xyzR125[1];
3654 (*varvec)[2] = xyzR125[2];
3657 (*varvec)[5] = d0z0[0];
3658 (*varvec)[6] = (
Int_t)isconv + 2*(
Int_t)isconv_like;
3659 (*varvec)[7] = mcetype;
3662 m_ReservoirE[
fPoolIndex][nextRes].push_back(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3681 for(
Int_t i=0;i<26;i++){
3684 for(
Int_t i=0;i<16;i++){
3779 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3785 fCandidateVariableNames[ 0]=
"V0Px";
3786 fCandidateVariableNames[ 1]=
"V0Py";
3787 fCandidateVariableNames[ 2]=
"V0Pz";
3788 fCandidateVariableNames[ 3]=
"MassLambda";
3789 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3790 fCandidateVariableNames[ 5]=
"ProtonPx";
3791 fCandidateVariableNames[ 6]=
"ProtonPy";
3792 fCandidateVariableNames[ 7]=
"ProtonPz";
3793 fCandidateVariableNames[ 8]=
"PionPx";
3794 fCandidateVariableNames[ 9]=
"PionPy";
3795 fCandidateVariableNames[10]=
"PionPz";
3796 fCandidateVariableNames[11]=
"RfidV0";
3797 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3798 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3799 fCandidateVariableNames[14]=
"DcaV0daughters";
3800 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3801 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3802 fCandidateVariableNames[17]=
"MassK0Short";
3803 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3804 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3805 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3806 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3807 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3808 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3809 fCandidateVariableNames[24]=
"V0PdgCode";
3810 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3811 fCandidateVariableNames[26]=
"mcv0px";
3812 fCandidateVariableNames[27]=
"mcv0py";
3813 fCandidateVariableNames[28]=
"mcv0pz";
3814 fCandidateVariableNames[29]=
"mcv0vertx";
3815 fCandidateVariableNames[30]=
"mcv0verty";
3816 fCandidateVariableNames[31]=
"mcv0vertz";
3817 fCandidateVariableNames[32]=
"EvNumber";
3818 fCandidateVariableNames[33]=
"Centrality";
3819 fCandidateVariableNames[34]=
"PrimVertZ";
3820 fCandidateVariableNames[35]=
"RunNumber";
3836 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3837 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3838 if(!cptrack)
return;
3839 if(!cntrack)
return;
3858 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3860 fVtx1->GetXYZ(posVtx);
3861 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3862 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3864 cont_arm[0] = v0->Pt();
3865 cont_arm[1] = v0->PtArmV0();
3866 cont_arm[2] = v0->AlphaV0();
3887 Double_t phi_alice = atan2(momv0y,momv0x);
3888 if(phi_alice<0.) phi_alice += 2 * M_PI;
3894 fVtx1->GetXYZ(posVtx);
3895 TLorentzVector *lv =
new TLorentzVector();
3896 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3897 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
3898 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
3900 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
3904 TVector *varvec =
new TVector(17);
3905 (*varvec)[0] = xyzR125pr[0];
3906 (*varvec)[1] = xyzR125pr[1];
3907 (*varvec)[2] = xyzR125pr[2];
3908 (*varvec)[3] = xyzR125pi[0];
3909 (*varvec)[4] = xyzR125pi[1];
3910 (*varvec)[5] = xyzR125pi[2];
3911 (*varvec)[6] = v0->DecayVertexV0X();
3912 (*varvec)[7] = v0->DecayVertexV0Y();
3913 (*varvec)[8] = v0->DecayVertexV0Z();
3914 (*varvec)[9] = v0->MomPosX();
3915 (*varvec)[10] = v0->MomPosY();
3916 (*varvec)[11] = v0->MomPosZ();
3917 (*varvec)[12] = v0->MomNegX();
3918 (*varvec)[13] = v0->MomNegY();
3919 (*varvec)[14] = v0->MomNegZ();
3920 (*varvec)[15] = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3921 (*varvec)[16] = v0->DcaV0ToPrimVertex();
3924 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
3928 TVector *varvec =
new TVector(17);
3929 (*varvec)[0] = xyzR125pr[0];
3930 (*varvec)[1] = xyzR125pr[1];
3931 (*varvec)[2] = xyzR125pr[2];
3932 (*varvec)[3] = xyzR125pi[0];
3933 (*varvec)[4] = xyzR125pi[1];
3934 (*varvec)[5] = xyzR125pi[2];
3935 (*varvec)[6] = v0->DecayVertexV0X();
3936 (*varvec)[7] = v0->DecayVertexV0Y();
3937 (*varvec)[8] = v0->DecayVertexV0Z();
3938 (*varvec)[9] = v0->MomNegX();
3939 (*varvec)[10] = v0->MomNegY();
3940 (*varvec)[11] = v0->MomNegZ();
3941 (*varvec)[12] = v0->MomPosX();
3942 (*varvec)[13] = v0->MomPosY();
3943 (*varvec)[14] = v0->MomPosZ();
3944 (*varvec)[15] = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3945 (*varvec)[16] = v0->DcaV0ToPrimVertex();
3952 Int_t v0pdgcode = -9999;
3953 Int_t v0motherpdgcode = -9999;
3960 Int_t pdgarray_v0[100], labelarray_v0[100], ngen_v0;
3963 Int_t pdgdgv0[2]={2212,211};
3964 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
3978 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
3979 if(!mcv0trk)
return;
3984 for(
Int_t i=0;i<100;i++){
3985 pdgarray_v0[i] = -9999;
3986 labelarray_v0[i] = -9999;
3991 Bool_t hfv0_flag = kFALSE;
3992 Bool_t promptv0_flag = kFALSE;
3993 v0pdgcode = mcv0trk->GetPdgCode();
3994 Int_t labv0mother = mcv0trk->GetMother();
3996 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3998 v0motherpdgcode = motherv0->GetPdgCode();
3999 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
4002 if(abs(v0motherpdgcode)==3322){
4004 }
else if(abs(v0motherpdgcode)==3312){
4006 }
else if(abs(v0motherpdgcode)==3334){
4010 promptv0_flag = kTRUE;
4017 mcv0px = mcv0trk->Px();
4018 mcv0py = mcv0trk->Py();
4019 mcv0pz = mcv0trk->Pz();
4021 Int_t labptrk = cptrack->GetLabel();
4023 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
4025 mcv0vertx = mcv0ptrk->Xv();
4026 mcv0verty = mcv0ptrk->Yv();
4027 mcv0vertz = mcv0ptrk->Zv();
4032 for(
Int_t i=0;i<36;i++){
4036 for(
Int_t i=0;i<16;i++){
4040 Bool_t isparticle = kFALSE;
4180 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
4186 fCandidateVariableNames[ 0]=
"Centrality";
4187 fCandidateVariableNames[ 1]=
"DecayType";
4188 fCandidateVariableNames[ 2]=
"LcPx";
4189 fCandidateVariableNames[ 3]=
"LcPy";
4190 fCandidateVariableNames[ 4]=
"LcPz";
4191 fCandidateVariableNames[ 5]=
"ElePx";
4192 fCandidateVariableNames[ 6]=
"ElePy";
4193 fCandidateVariableNames[ 7]=
"ElePz";
4194 fCandidateVariableNames[ 8]=
"V0Px";
4195 fCandidateVariableNames[ 9]=
"V0Py";
4196 fCandidateVariableNames[10]=
"V0Pz";
4210 if(!mcepart)
return;
4211 if(!mcv0part)
return;
4213 for(
Int_t i=0;i<11;i++){
4232 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
4236 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
4238 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
4241 cont[0] = InvMassEleLambda;
4242 cont[1] = mcpart->Pt();
4245 cont2[0] = InvMassEleLambda;
4246 cont2[1] = mcepart->Pt();
4249 cont_eleptvseta[0] = mcepart->Pt();
4250 cont_eleptvseta[1] = mcepart->Eta();
4253 cont_eleptvslambdapt[0] = mcepart->Pt();
4254 cont_eleptvslambdapt[1] = mcv0part->Pt();
4256 Double_t cont_eleptvslambdaptvslcpt[4];
4257 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
4258 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
4259 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
4262 contmc[0] = mcpart->Pt();
4263 contmc[1] = mcpart->Y();
4266 contmcele[0] = mcepart->Pt();
4267 contmcele[1] = mcepart->Eta();
4270 contmc_withele[0] = mcpart->Pt();
4271 contmc_withele[1] = mcpart->Y();
4272 contmc_withele[2] = mcepart->Y();
4276 esdcuts->GetEtaRange(etamin,etamax);
4283 if(fabs(mcepart->Eta())<0.8){
4290 if(fabs(mcepart->Eta())<etamax){
4294 if(InvMassEleLambda<2.3){
4300 if(fabs(mcpart->Y())<0.7){
4301 if(InvMassEleLambda<2.3){
4306 }
else if(decaytype==1){
4314 if(fabs(mcepart->Eta())<etamax){
4318 if(InvMassEleLambda<2.3){
4324 }
else if(decaytype==2){
4332 if(fabs(mcepart->Eta())<etamax){
4336 if(InvMassEleLambda<2.3){
4342 }
else if(decaytype==3){
4346 }
else if(decaytype==10){
4355 if(fabs(mcepart->Eta())<0.8){
4358 }
else if(decaytype==11){
4374 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
4380 fCandidateVariableNames[ 0]=
"Centrality";
4381 fCandidateVariableNames[ 1]=
"ElePx";
4382 fCandidateVariableNames[ 2]=
"ElePy";
4383 fCandidateVariableNames[ 3]=
"ElePz";
4384 fCandidateVariableNames[ 4]=
"ElePdgCode";
4385 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
4386 fCandidateVariableNames[ 6]=
"RunNumber";
4387 fCandidateVariableNames[ 7]=
"EvNumber";
4400 if(!mcepart)
return;
4402 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
4407 Int_t semi_flag = 0;
4408 if(pdgarray_ele[0]>400&&pdgarray_ele[0]<450) semi_flag = 1;
4409 if(pdgarray_ele[0]>4000&&pdgarray_ele[0]<4500) semi_flag = 1;
4410 if(pdgarray_ele[0]>500&&pdgarray_ele[0]<550) semi_flag = 2;
4411 if(pdgarray_ele[0]>5000&&pdgarray_ele[0]<5500) semi_flag = 2;
4414 contmc[0] = mcepart->Pt();
4415 contmc[1] = mcepart->Eta();
4418 if(semi_flag==1 && !ele_from_bottom){
4421 if(semi_flag==1 && ele_from_bottom){
4428 Bool_t hfe_flag = kFALSE;
4429 Int_t labemother = mcepart->GetMother();
4430 Int_t pdgmotherele = -9999;
4432 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
4433 pdgmotherele = motherele->GetPdgCode();
4434 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
4439 if(!hfe_flag)
return;
4443 for(
Int_t i=0;i<8;i++){
4467 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
4473 fCandidateVariableNames[ 0]=
"Centrality";
4474 fCandidateVariableNames[ 1]=
"V0Px";
4475 fCandidateVariableNames[ 2]=
"V0Py";
4476 fCandidateVariableNames[ 3]=
"V0Pz";
4477 fCandidateVariableNames[ 4]=
"V0PdgCode";
4478 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
4479 fCandidateVariableNames[ 6]=
"RunNumber";
4480 fCandidateVariableNames[ 7]=
"EvNumber";
4493 if(!mcv0part)
return;
4495 for(
Int_t i=0;i<8;i++){
4499 Bool_t hfv0_flag = kFALSE;
4500 Int_t labv0mother = mcv0part->GetMother();
4501 Int_t pdgmotherv0 = -9999;
4503 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
4505 pdgmotherv0 = motherv0->GetPdgCode();
4506 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
4511 if(!hfv0_flag)
return;
4514 contmc[0] = mcv0part->Pt();
4515 contmc[1] = mcv0part->Eta();
4539 fCEvents =
new TH1F(
"fCEvents",
"conter",20,-0.5,19.5);
4541 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
4542 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
4543 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
4544 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
4545 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
4546 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
4547 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
4548 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
4549 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
4550 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
4551 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
4552 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
4553 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
4554 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
4556 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
4557 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
4559 fCEvents->GetXaxis()->SetBinLabel(19,
"No RP");
4560 fCEvents->GetXaxis()->SetBinLabel(20,
"No Qvector (only for TPC event plane)");
4562 fCEvents->GetYaxis()->SetTitle(
"counts");
4564 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
4566 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
4567 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
4568 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
4569 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
4570 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
4571 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
4573 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
4574 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
4575 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
4576 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
4577 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
4579 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
4580 fHEventPlane =
new TH1F(
"fHEventPlane",
"conter",100,-3.14,3.14);
4581 fHNTrackletvsZ =
new TH2F(
"fHNTrackletvsZ",
"N_{tracklet} vs z",30,-15.,15.,120,-0.5,119.5);
4603 Int_t bins_base[3]= {23,100 ,10};
4604 Double_t xmin_base[3]={1.1,0 ,0.00};
4605 Double_t xmax_base[3]={5.7,20. ,100};
4606 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
4682 Int_t bins_base_elept[3]= {10,100 ,10};
4683 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
4684 Double_t xmax_base_elept[3]={3.1,10. ,100};
4756 Int_t bins_eleptvseta[3]= {10,10 ,10};
4757 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4758 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4760 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4762 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4764 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4766 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4768 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4770 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4773 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4774 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4775 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4777 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4779 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4781 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4783 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4785 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4787 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4789 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4792 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4793 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4794 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4795 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4797 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4800 Int_t bins_allpt[4]= {10,10,20,20};
4801 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4802 Double_t xmax_allpt[4]={20.,20.,10.,10};
4813 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4814 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4815 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4817 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4819 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4821 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4823 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4825 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4827 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4893 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4899 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4901 fHistoBachPtMCS =
new TH1F(
"fHistoBachPtMCS",
"Bachelor p_{T}",100,0.0,5.0);
4905 fHistod0Bach =
new TH1F(
"fHistod0Bach",
"Bachelor d_{0}",100,-0.5,0.5);
4934 Int_t bins_armenteros[3]= {20 ,100 ,200};
4935 Double_t xmin_armenteros[3]={0.,0. ,-1.};
4936 Double_t xmax_armenteros[3]={20.,1. ,1.};
4937 fHistoLambdaArmenteros =
new THnSparseF(
"fHistoLambdaArmenteros",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4939 fHistoLambdaArmenterosMCS =
new THnSparseF(
"fHistoLambdaArmenterosMCS",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4961 for(
Int_t i=0;i<8;i++){
4970 Int_t bins_lcmcgen_withelea[3]= {20 ,100 ,100};
4971 Double_t xmin_lcmcgen_withelea[3]={0.,-5.,-5.};
4972 Double_t xmax_lcmcgen_withelea[3]={20.,5.,5.};
4973 fHistoLcMCGenWithEleA =
new THnSparseF(
"fHistoLcMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4978 fHistoLbMCGenWithEleA =
new THnSparseF(
"fHistoLbMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4983 Int_t bins_lbmcgen_withlc[3]= {50 ,100 ,100};
4984 Double_t xmin_lbmcgen_withlc[3]={0.,-5.,-5.};
4985 Double_t xmax_lbmcgen_withlc[3]={50.,5.,5.};
4986 fHistoLbMCGenWithLambdac =
new THnSparseF(
"fHistoLbMCGenWithLambdac",
"",3,bins_lbmcgen_withlc,xmin_lbmcgen_withlc,xmax_lbmcgen_withlc);
4989 Int_t bins_lcmcgen[3]= {100 ,20 ,10};
4990 Double_t xmin_lcmcgen[3]={0.,-1.0 ,0.0};
4991 Double_t xmax_lcmcgen[3]={20.,1.0 ,100};
4992 Int_t bins_lbmcgen[3]= {100 ,20 ,10};
4993 Double_t xmin_lbmcgen[3]={0.,-1.0 ,0.0};
4994 Double_t xmax_lbmcgen[3]={50.,1.0 ,100};
4995 fHistoLcMCGen =
new THnSparseF(
"fHistoLcMCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4997 fHistoLcMCGen1 =
new THnSparseF(
"fHistoLcMCGen1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4999 fHistoLcMCGen2 =
new THnSparseF(
"fHistoLcMCGen2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5001 fHistoLbMCGen =
new THnSparseF(
"fHistoLbMCGen",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5003 fHistoLbMCGen1 =
new THnSparseF(
"fHistoLbMCGen1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5005 fHistoLbMCGen2 =
new THnSparseF(
"fHistoLbMCGen2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5007 fHistoFeedDownXic0MCGen =
new THnSparseF(
"fHistoFeedDownXic0MCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5026 fHistoLcMCS =
new THnSparseF(
"fHistoLcMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5028 fHistoLcMCS1 =
new THnSparseF(
"fHistoLcMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5030 fHistoLcMCS2 =
new THnSparseF(
"fHistoLcMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5032 fHistoLbMCS =
new THnSparseF(
"fHistoLbMCS",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5034 fHistoLbMCS1 =
new THnSparseF(
"fHistoLbMCS1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5036 fHistoLbMCS2 =
new THnSparseF(
"fHistoLbMCS2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
5038 fHistoFeedDownXic0MCS =
new THnSparseF(
"fHistoFeedDownXic0MCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5040 fHistoFeedDownXic0MCS1 =
new THnSparseF(
"fHistoFeedDownXic0MCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5042 fHistoFeedDownXic0MCS2 =
new THnSparseF(
"fHistoFeedDownXic0MCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5050 fHistoFeedDownSigmaMCS =
new THnSparseF(
"fHistoFeedDownSigmaMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5052 fHistoFeedDownSigmaMCS1 =
new THnSparseF(
"fHistoFeedDownSigmaMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5054 fHistoFeedDownSigmaMCS2 =
new THnSparseF(
"fHistoFeedDownSigmaMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
5057 Int_t bins_elemcgen[3]= {100 ,20 ,10};
5058 Double_t xmin_elemcgen[3]={0.,-1.0 ,0.0};
5059 Double_t xmax_elemcgen[3]={10.,1.0 ,100};
5060 fHistoElectronMCGen =
new THnSparseF(
"fHistoElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5064 fHistoCharmElectronMCGen =
new THnSparseF(
"fHistoCharmElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5066 fHistoLcElectronMCGen =
new THnSparseF(
"fHistoLcElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5068 fHistoLcElectronMCGen1 =
new THnSparseF(
"fHistoLcElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5070 fHistoLcElectronMCGen2 =
new THnSparseF(
"fHistoLcElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5072 fHistoLbElectronMCGen =
new THnSparseF(
"fHistoLbElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5074 fHistoLbElectronMCGen1 =
new THnSparseF(
"fHistoLbElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5076 fHistoLbElectronMCGen2 =
new THnSparseF(
"fHistoLbElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5096 fHistoLcElectronMCS =
new THnSparseF(
"fHistoLcElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5098 fHistoLcElectronMCS1 =
new THnSparseF(
"fHistoLcElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5100 fHistoLcElectronMCS2 =
new THnSparseF(
"fHistoLcElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5102 fHistoLbElectronMCS =
new THnSparseF(
"fHistoLbElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5104 fHistoLbElectronMCS1 =
new THnSparseF(
"fHistoLbElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
5106 fHistoLbElectronMCS2 =
new THnSparseF(
"fHistoLbElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);