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> 49 #include "AliMCEvent.h" 50 #include "AliAnalysisManager.h" 51 #include "AliAODMCHeader.h" 52 #include "AliAODHandler.h" 54 #include "AliExternalTrackParam.h" 55 #include "AliAODVertex.h" 56 #include "AliESDVertex.h" 57 #include "AliAODRecoDecay.h" 60 #include "AliESDtrack.h" 61 #include "AliAODTrack.h" 63 #include "AliAODcascade.h" 64 #include "AliAODMCParticle.h" 65 #include "AliAnalysisTaskSE.h" 67 #include "AliPIDResponse.h" 68 #include "AliPIDCombined.h" 69 #include "AliTOFPIDResponse.h" 71 #include "AliInputEventHandler.h" 72 #include "AliESDtrackCuts.h" 73 #include "AliNeutralTrackParam.h" 74 #include "AliKFParticle.h" 75 #include "AliKFVertex.h" 76 #include "AliExternalTrackParam.h" 77 #include "AliESDtrack.h" 78 #include "AliCentrality.h" 79 #include "AliVertexerTracks.h" 80 #include "AliEventPoolManager.h" 83 #include "AliQnCorrectionsManager.h" 84 #include "AliQnCorrectionsQnVector.h" 85 #include "AliAnalysisTaskFlowVectorCorrections.h" 106 fHNTrackletCorrvsZ(0),
108 fIsEventSelected(kFALSE),
109 fWriteVariableTree(kFALSE),
110 fWriteEachVariableTree(kFALSE),
111 fWriteMCVariableTree(kFALSE),
113 fEleVariablesTree(0),
115 fSingleVariablesTree(0),
117 fMCEleVariablesTree(0),
118 fMCV0VariablesTree(0),
119 fMCGenPairVariablesTree(0),
120 fCorrelationVariablesTree(0),
121 fReconstructPrimVert(kFALSE),
127 fCandidateVariables(),
128 fCandidateEleVariables(),
129 fCandidateV0Variables(),
130 fCandidateSingleVariables(),
131 fCandidateMCVariables(),
132 fCandidateMCEleVariables(),
133 fCandidateMCV0Variables(),
134 fCandidateMCGenPairVariables(),
135 fCorrelationVariables(),
147 fUseCentralityV0M(kFALSE),
148 fUseCentralitySPDTracklet(kFALSE),
152 fMCDoPairAnalysis(kFALSE),
153 fHistoEleLambdaMass(0),
154 fHistoEleLambdaMassRS(0),
155 fHistoEleLambdaMassRS1(0),
156 fHistoEleLambdaMassRS2(0),
157 fHistoEleLambdaMassWS(0),
158 fHistoEleLambdaMassWS1(0),
159 fHistoEleLambdaMassWS2(0),
160 fHistoEleLambdaMassRSMix(0),
161 fHistoEleLambdaMassRSMix1(0),
162 fHistoEleLambdaMassRSMix2(0),
163 fHistoEleLambdaMassWSMix(0),
164 fHistoEleLambdaMassWSMix1(0),
165 fHistoEleLambdaMassWSMix2(0),
166 fHistoEleLambdaMassRSSide(0),
167 fHistoEleLambdaMassRSSide1(0),
168 fHistoEleLambdaMassRSSide2(0),
169 fHistoEleLambdaMassWSSide(0),
170 fHistoEleLambdaMassWSSide1(0),
171 fHistoEleLambdaMassWSSide2(0),
172 fHistoEleLambdaMassRSAway(0),
173 fHistoEleLambdaMassRS1Away(0),
174 fHistoEleLambdaMassRS2Away(0),
175 fHistoEleLambdaMassWSAway(0),
176 fHistoEleLambdaMassWS1Away(0),
177 fHistoEleLambdaMassWS2Away(0),
178 fHistoEleLambdaMassRSMixAway(0),
179 fHistoEleLambdaMassRSMix1Away(0),
180 fHistoEleLambdaMassRSMix2Away(0),
181 fHistoEleLambdaMassWSMixAway(0),
182 fHistoEleLambdaMassWSMix1Away(0),
183 fHistoEleLambdaMassWSMix2Away(0),
184 fHistoEleLambdaMassRSSideAway(0),
185 fHistoEleLambdaMassRSSide1Away(0),
186 fHistoEleLambdaMassRSSide2Away(0),
187 fHistoEleLambdaMassWSSideAway(0),
188 fHistoEleLambdaMassWSSide1Away(0),
189 fHistoEleLambdaMassWSSide2Away(0),
190 fHistoEleLambdaMassvsElePtRS(0),
191 fHistoEleLambdaMassvsElePtWS(0),
192 fHistoEleLambdaMassvsElePtRSMix(0),
193 fHistoEleLambdaMassvsElePtWSMix(0),
194 fHistoEleLambdaMassvsElePtRSSide(0),
195 fHistoEleLambdaMassvsElePtWSSide(0),
196 fHistoEleLambdaMassvsElePtRS1(0),
197 fHistoEleLambdaMassvsElePtWS1(0),
198 fHistoEleLambdaMassvsElePtRSMix1(0),
199 fHistoEleLambdaMassvsElePtWSMix1(0),
200 fHistoEleLambdaMassvsElePtRSSide1(0),
201 fHistoEleLambdaMassvsElePtWSSide1(0),
202 fHistoEleLambdaMassvsElePtRS2(0),
203 fHistoEleLambdaMassvsElePtWS2(0),
204 fHistoEleLambdaMassvsElePtRSMix2(0),
205 fHistoEleLambdaMassvsElePtWSMix2(0),
206 fHistoEleLambdaMassvsElePtRSSide2(0),
207 fHistoEleLambdaMassvsElePtWSSide2(0),
212 fHistoEleLambdaMassMCS(0),
213 fHistoEleLambdaMassLambdabMCS(0),
214 fHistoEleLambdaMassXibMCS(0),
215 fHistoEleLambdaMassMCGen(0),
216 fHistoEleLambdaMassvsElePtMCS(0),
217 fHistoEleLambdaMassvsElePtMCGen(0),
218 fHistoEleLambdaMassvsElePtMCS1(0),
219 fHistoEleLambdaMassvsElePtMCGen1(0),
220 fHistoEleLambdaMassvsElePtMCS2(0),
221 fHistoEleLambdaMassvsElePtMCGen2(0),
224 fHistoElePtvsEtaRS(0),
225 fHistoElePtvsEtaWS(0),
226 fHistoElePtvsEtaRSMix(0),
227 fHistoElePtvsEtaWSMix(0),
228 fHistoElePtvsEtaMCS(0),
229 fHistoElePtvsEtaMCGen(0),
230 fHistoElePtvsLambdaPtRS(0),
231 fHistoElePtvsLambdaPtWS(0),
232 fHistoElePtvsLambdaPtRSMix(0),
233 fHistoElePtvsLambdaPtWSMix(0),
234 fHistoElePtvsLambdaPtMCS(0),
235 fHistoElePtvsLambdaPtvsLcPtMCS(0),
236 fHistoElePtvsLambdaPtMCGen(0),
237 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
238 fHistoElePtvsLambdaPtMCLcGen(0),
239 fHistoElePtvsd0RS(0),
240 fHistoElePtvsd0WS(0),
241 fHistoElePtvsd0RSMix(0),
242 fHistoElePtvsd0WSMix(0),
243 fHistoElePtvsd0MCS(0),
244 fHistoElePtvsd0PromptMCS(0),
245 fHistoElePtvsd0BFeeddownMCS(0),
246 fHistoEleLambdaMassFeeddownXic0MCS(0),
247 fHistoEleLambdaMassFeeddownXic0MCGen(0),
248 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
249 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
250 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
251 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
252 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
253 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
254 fHistoElePtFeeddownXic0MCS(0),
255 fHistoElePtFeeddownXic0MCGen(0),
256 fHistoElePtvsEtaFeeddownXic0MCS(0),
257 fHistoElePtvsEtaFeeddownXic0MCGen(0),
258 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
259 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
260 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
261 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
262 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
263 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
264 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
265 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
266 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
267 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
268 fHistoElePtFeeddownXicPlusMCS(0),
269 fHistoElePtFeeddownXicPlusMCGen(0),
270 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
271 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
272 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
273 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
276 fHistoBachPtMCGen(0),
277 fHistoBachFracSharedITS(0),
278 fHistoBachChi2NclsITS(0),
280 fHistoLambdaMassvsPt(0),
281 fHistoLambdaMassvsPtMCS(0),
282 fHistoLambdaMassvsPtMCGen(0),
283 fHistoLambdaPtvsEta(0),
284 fHistoLambdaPtvsEtaMCS(0),
285 fHistoLambdaPtvsEtaMCGen(0),
286 fHistoLambdaArmenteros(0),
287 fHistoLambdaArmenterosMCS(0),
288 fHistoLambdaPtvsDl(0),
289 fHistoLambdaPtvsDlSide(0),
290 fHistoLambdaPtvsDlMCS(0),
291 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
292 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
293 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
294 fHistoK0sMassvsPt(0),
295 fHistoElectronTPCPID(0),
296 fHistoElectronTOFPID(0),
297 fHistoElectronITSPID(0),
298 fHistoElectronTPCSelPID(0),
299 fHistoElectronTOFSelPID(0),
300 fHistoElectronITSSelPID(0),
301 fHistoElectronTPCPIDSelTOF(0),
302 fHistoElectronTPCPIDSelTOFSmallEta(0),
303 fHistoElectronTPCPIDSelTOFLargeEta(0),
304 fHistoElectronITSPIDSelTPCTOF(0),
305 fHistoElectronITSPIDSelTPCTOFSmallEta(0),
306 fHistoElectronITSPIDSelTPCTOFLargeEta(0),
307 fHistoMassConversionsMin(0),
308 fHistoMassConversionsSameSignMin(0),
309 fHistoElectronQovPtvsPhi(0),
310 fHistoLambdaQovPtvsPhi(0),
317 fHistoLcMCGenWithEleA(0),
318 fHistoLcMCGenWithEleB(0),
319 fHistoLbMCGenWithEleA(0),
320 fHistoLbMCGenWithEleB(0),
321 fHistoLbMCGenWithLambdac(0),
328 fHistoFeedDownXic0MCGen(0),
329 fHistoFeedDownXic0MCGen1(0),
330 fHistoFeedDownXic0MCGen2(0),
331 fHistoFeedDownXic0MCS(0),
332 fHistoFeedDownXic0MCS1(0),
333 fHistoFeedDownXic0MCS2(0),
334 fHistoFeedDownXicPlusMCGen(0),
335 fHistoFeedDownXicPlusMCGen1(0),
336 fHistoFeedDownXicPlusMCGen2(0),
337 fHistoFeedDownXicPlusMCS(0),
338 fHistoFeedDownXicPlusMCS1(0),
339 fHistoFeedDownXicPlusMCS2(0),
340 fHistoFeedDownSigmaMCGen(0),
341 fHistoFeedDownSigmaMCGen1(0),
342 fHistoFeedDownSigmaMCGen2(0),
343 fHistoFeedDownSigmaMCS(0),
344 fHistoFeedDownSigmaMCS1(0),
345 fHistoFeedDownSigmaMCS2(0),
346 fHistoLcElectronMCGen(0),
347 fHistoLcElectronMCGen1(0),
348 fHistoLcElectronMCGen2(0),
349 fHistoLcElectronMCS(0),
350 fHistoLcElectronMCS1(0),
351 fHistoLcElectronMCS2(0),
352 fHistoLbElectronMCGen(0),
353 fHistoLbElectronMCGen1(0),
354 fHistoLbElectronMCGen2(0),
355 fHistoLbElectronMCS(0),
356 fHistoLbElectronMCS1(0),
357 fHistoLbElectronMCS2(0),
358 fHistoElectronFeedDownXic0MCGen(0),
359 fHistoElectronFeedDownXic0MCGen1(0),
360 fHistoElectronFeedDownXic0MCGen2(0),
361 fHistoElectronFeedDownXic0MCS(0),
362 fHistoElectronFeedDownXic0MCS1(0),
363 fHistoElectronFeedDownXic0MCS2(0),
364 fHistoElectronFeedDownXicPlusMCGen(0),
365 fHistoElectronFeedDownXicPlusMCGen1(0),
366 fHistoElectronFeedDownXicPlusMCGen2(0),
367 fHistoElectronFeedDownXicPlusMCS(0),
368 fHistoElectronFeedDownXicPlusMCS1(0),
369 fHistoElectronFeedDownXicPlusMCS2(0),
370 fHistoElectronFeedDownXibMCGen(0),
371 fHistoElectronFeedDownXibMCGen1(0),
372 fHistoElectronFeedDownXibMCGen2(0),
373 fHistoElectronFeedDownXibMCS(0),
374 fHistoElectronFeedDownXibMCS1(0),
375 fHistoElectronFeedDownXibMCS2(0),
376 fHistoElectronMCGen(0),
377 fHistoBottomElectronMCGen(0),
378 fHistoCharmElectronMCGen(0),
379 fHistoLambdaMCGen(0),
380 fHistoElePtvsV0dlRS(0),
381 fHistoElePtvsV0dlRS1(0),
382 fHistoElePtvsV0dlRS2(0),
383 fHistoElePtvsV0dlRSSide(0),
384 fHistoElePtvsV0dlRSSide1(0),
385 fHistoElePtvsV0dlRSSide2(0),
386 fHistoElePtvsV0dlRSMix(0),
387 fHistoElePtvsV0dlRSMix1(0),
388 fHistoElePtvsV0dlRSMix2(0),
389 fHistoElePtvsV0dlWS(0),
390 fHistoElePtvsV0dlWS1(0),
391 fHistoElePtvsV0dlWS2(0),
392 fHistoElePtvsV0dlWSSide(0),
393 fHistoElePtvsV0dlWSSide1(0),
394 fHistoElePtvsV0dlWSSide2(0),
395 fHistoElePtvsV0dlWSMix(0),
396 fHistoElePtvsV0dlWSMix1(0),
397 fHistoElePtvsV0dlWSMix2(0),
398 fHistoElePtvsV0dlMCS(0),
399 fHistoElePtvsV0dlMCS1(0),
400 fHistoElePtvsV0dlMCS2(0),
401 fHistoElePtvsV0dlFeedDownXic0MCS(0),
402 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
403 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
404 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
405 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
406 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
407 fHistoElePtvsV0dcaRS(0),
408 fHistoElePtvsV0dcaRS1(0),
409 fHistoElePtvsV0dcaRS2(0),
410 fHistoElePtvsV0dcaRSSide(0),
411 fHistoElePtvsV0dcaRSSide1(0),
412 fHistoElePtvsV0dcaRSSide2(0),
413 fHistoElePtvsV0dcaRSMix(0),
414 fHistoElePtvsV0dcaRSMix1(0),
415 fHistoElePtvsV0dcaRSMix2(0),
416 fHistoElePtvsV0dcaWS(0),
417 fHistoElePtvsV0dcaWS1(0),
418 fHistoElePtvsV0dcaWS2(0),
419 fHistoElePtvsV0dcaWSSide(0),
420 fHistoElePtvsV0dcaWSSide1(0),
421 fHistoElePtvsV0dcaWSSide2(0),
422 fHistoElePtvsV0dcaWSMix(0),
423 fHistoElePtvsV0dcaWSMix1(0),
424 fHistoElePtvsV0dcaWSMix2(0),
425 fHistoElePtvsV0dcaMCS(0),
426 fHistoElePtvsV0dcaMCS1(0),
427 fHistoElePtvsV0dcaMCS2(0),
428 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
429 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
430 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
431 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
432 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
433 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
434 fHistoEleLambdaPtvsV0dlRS(0),
435 fHistoEleLambdaPtvsV0dlRS1(0),
436 fHistoEleLambdaPtvsV0dlRS2(0),
437 fHistoEleLambdaPtvsV0dlRSSide(0),
438 fHistoEleLambdaPtvsV0dlRSSide1(0),
439 fHistoEleLambdaPtvsV0dlRSSide2(0),
440 fHistoEleLambdaPtvsV0dlRSMix(0),
441 fHistoEleLambdaPtvsV0dlRSMix1(0),
442 fHistoEleLambdaPtvsV0dlRSMix2(0),
443 fHistoEleLambdaPtvsV0dlWS(0),
444 fHistoEleLambdaPtvsV0dlWS1(0),
445 fHistoEleLambdaPtvsV0dlWS2(0),
446 fHistoEleLambdaPtvsV0dlWSSide(0),
447 fHistoEleLambdaPtvsV0dlWSSide1(0),
448 fHistoEleLambdaPtvsV0dlWSSide2(0),
449 fHistoEleLambdaPtvsV0dlWSMix(0),
450 fHistoEleLambdaPtvsV0dlWSMix1(0),
451 fHistoEleLambdaPtvsV0dlWSMix2(0),
452 fHistoEleLambdaPtvsV0dlMCS(0),
453 fHistoEleLambdaPtvsV0dlMCS1(0),
454 fHistoEleLambdaPtvsV0dlMCS2(0),
455 fHistoEleLambdaPtvsV0dlPromptMCS(0),
456 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
457 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
458 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
459 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
460 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
461 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
462 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
463 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
464 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
465 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
466 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
467 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
468 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
469 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
470 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
471 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
472 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
473 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
474 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
475 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
476 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
477 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
478 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
479 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
480 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
481 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
482 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
483 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
484 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
485 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
486 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
487 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
488 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
489 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
490 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
491 fHistoEleLambdaPtvsV0dlXibMCS(0),
492 fHistoEleLambdaPtvsV0dlXibMCS1(0),
493 fHistoEleLambdaPtvsV0dlXibMCS2(0),
494 fHistoEleLambdaPtvsV0dlRSAway(0),
495 fHistoEleLambdaPtvsV0dlRS1Away(0),
496 fHistoEleLambdaPtvsV0dlRS2Away(0),
497 fHistoEleLambdaPtvsV0dlRSSideAway(0),
498 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
499 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
500 fHistoEleLambdaPtvsV0dlRSMixAway(0),
501 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
502 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
503 fHistoEleLambdaPtvsV0dlWSAway(0),
504 fHistoEleLambdaPtvsV0dlWS1Away(0),
505 fHistoEleLambdaPtvsV0dlWS2Away(0),
506 fHistoEleLambdaPtvsV0dlWSSideAway(0),
507 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
508 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
509 fHistoEleLambdaPtvsV0dlWSMixAway(0),
510 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
511 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
512 fHistoCorrelationVariablesvsEleLambdaPt(0),
513 fHistoCorrelationVariablesvsEleLambdaPtMix(0),
514 fHistoCorrelationVariablesvsEleLambdaPtMC(0),
515 fHistoCorrelationVariablesvsElePt(0),
516 fHistoCorrelationVariablesvsElePtMix(0),
517 fHistoCorrelationVariablesvsElePtMC(0),
518 fHistoCorrelationVariablesvsLambdaPt(0),
519 fHistoCorrelationVariablesvsLambdaPtMix(0),
520 fHistoCorrelationVariablesvsLambdaPtMC(0),
521 fHistoMassVariablesvsEleLambdaPt(0),
522 fHistoMassVariablesvsEleLambdaPtMix(0),
523 fHistoMassVariablesvsEleLambdaPtMC(0),
524 fHistoMassVariablesvsElePt(0),
525 fHistoMassVariablesvsElePtMix(0),
526 fHistoMassVariablesvsElePtMC(0),
527 fHistoMassVariablesvsLambdaPt(0),
528 fHistoMassVariablesvsLambdaPtMix(0),
529 fHistoMassVariablesvsLambdaPtMC(0),
530 fHistoResponseElePt(0),
531 fHistoResponseElePt1(0),
532 fHistoResponseElePt2(0),
533 fHistoResponseEleLambdaPt(0),
534 fHistoResponseEleLambdaPt1(0),
535 fHistoResponseEleLambdaPt2(0),
536 fHistoResponseEleLambdaPtFeeddownXic0(0),
537 fHistoResponseEleLambdaPtFeeddownXic01(0),
538 fHistoResponseEleLambdaPtFeeddownXic02(0),
539 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
540 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
541 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
542 fHistoResponseEleLambdaPtFeeddownSigma(0),
543 fHistoResponseEleLambdaPtFeeddownSigma1(0),
544 fHistoResponseEleLambdaPtFeeddownSigma2(0),
545 fHistoResponseEleLambdaPtLambdab(0),
546 fHistoResponseElePtLambdab(0),
547 fHistoResponseLambdaPtLambdab(0),
548 fHistoResponseLambdaPt(0),
549 fHistoResponseLambdaPtFeeddownXic0(0),
550 fHistoResponseLambdaPtFeeddownXicPlus(0),
551 fHistoResponseLambdaPtFeeddownSigma(0),
552 fHistoResponseLambdaPtvsEleLambdaPt(0),
553 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
554 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
555 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
556 fHistoEleLambdaPtvsRapidityRS(0),
557 fHistoEleLambdaPtvsRapidityWS(0),
558 fHistoEleLambdaPtvsRapidityMCS(0),
560 fHistoElectronPi0Total(0),
561 fHistoElectronPi0Tag(0),
563 fHistoElectronEtaTotal(0),
564 fHistoElectronEtaTag(0),
568 fHistonEvtvsRunNumber(0),
569 fHistonElevsRunNumber(0),
570 fHistonLambdavsRunNumber(0),
571 fHistoMCEventType(0),
572 fHistoMCDeltaPhiccbar(0),
574 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
575 fHistodPhiSdEtaSElectronProtonR125RS(0),
576 fHistodPhiSdEtaSElectronProtonR125WS(0),
577 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
578 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
579 fHistodPhiSdEtaSElectronPionR125RS(0),
580 fHistodPhiSdEtaSElectronPionR125WS(0),
581 fHistodPhiSdEtaSElectronPionR125RSMix(0),
582 fHistodPhiSdEtaSElectronPionR125WSMix(0),
584 fMixWithoutConversionFlag(kFALSE),
585 fNumberOfEventsForMixing (5),
591 fHistoPoolMonitor(0),
592 fHistoPoolIDCounter(0),
605 for(
Int_t i=0;i<17;i++){
610 for(
Int_t i=0;i<8;i++){
620 Bool_t writeVariableTree) :
1130 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
1132 for(
Int_t i=0;i<17;i++){
1137 for(
Int_t i=0;i<8;i++){
1143 DefineOutput(1,TList::Class());
1144 DefineOutput(2,TList::Class());
1145 DefineOutput(3,TList::Class());
1146 DefineOutput(4,TTree::Class());
1147 DefineOutput(5,TTree::Class());
1148 DefineOutput(6,TTree::Class());
1149 DefineOutput(7,TTree::Class());
1150 DefineOutput(8,AliNormalizationCounter::Class());
1151 DefineOutput(9,TTree::Class());
1152 DefineOutput(10,TTree::Class());
1153 DefineOutput(11,TTree::Class());
1154 DefineOutput(12,TTree::Class());
1162 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1272 if (fDebug > 1) AliInfo(
"Init");
1291 AliError(
"NO EVENT FOUND!");
1303 AliKFParticle::SetField(
fBzkG);
1304 if (TMath::Abs(
fBzkG)<0.001) {
1313 AliAODTracklets* tracklets=aodEvent->GetTracklets();
1314 Int_t nTr=tracklets->GetNumberOfTracklets();
1315 for(
Int_t iTr=0; iTr<nTr; iTr++){
1316 Double_t theta=tracklets->GetTheta(iTr);
1317 Double_t eta=-TMath::Log(TMath::Tan(theta/2.));
1318 if(eta>-1.0 && eta<1.0) countTr++;
1320 AliAODVertex *vtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1323 if(vtx1->GetNContributors()>0){
1347 TClonesArray *mcArray = 0;
1348 AliAODMCHeader *mcHeader=0;
1351 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1353 AliError(
"Could not find Monte-Carlo in AOD");
1359 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1361 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1366 Double_t zMCVertex = mcHeader->GetVtxZ();
1368 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1382 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1387 fVtx1->GetCovarianceMatrix(cov);
1388 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1392 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1402 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1403 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1404 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1405 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1406 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1422 if(countCorr>=0 && countCorr<=0)
fCentrality = 5.;
1423 else if(countCorr>=1 && countCorr<=8)
fCentrality = 15.;
1424 else if(countCorr>= 9 && countCorr<=13)
fCentrality = 25.;
1425 else if(countCorr>=14 && countCorr<=19)
fCentrality = 35.;
1426 else if(countCorr>=20 && countCorr<=30)
fCentrality = 45.;
1427 else if(countCorr>=31 && countCorr<=49)
fCentrality = 55.;
1432 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1442 AliEventplane *pl=aodEvent->GetEventplane();
1444 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec:no eventplane! v2 analysis without eventplane not possible!\n");
1463 fQ = pl->GetQVector();
1465 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: no Q vectors");
1472 AliQnCorrectionsManager *flowQnVectorMgr;
1473 AliAnalysisTaskFlowVectorCorrections *flowQnVectorTask =
1474 dynamic_cast<AliAnalysisTaskFlowVectorCorrections *
>(AliAnalysisManager::GetAnalysisManager()->GetTask(
"FlowQnVectorCorrections"));
1475 if (flowQnVectorTask != NULL) {
1476 flowQnVectorMgr = flowQnVectorTask->GetAliQnCorrectionsManager();
1479 AliFatal(
"This task needs the Flow Qn vector corrections framework and it is not present. Aborting!!!");
1483 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: No such event plane method");
1486 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: No such event plane method");
1489 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: No such event plane method");
1492 const AliQnCorrectionsQnVector *myQnVector;
1493 myQnVector = flowQnVectorMgr->GetDetectorQnVector(
"TPC");
1495 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: no Q vectors");
1499 Double_t ep_tpc = myQnVector->EventPlane(2);
1508 Int_t runnumber_offset = 0;
1509 Int_t runnumber = aodEvent->GetRunNumber();
1510 if(runnumber<=131000&&runnumber>=114000){
1511 runnumber_offset = 114000;
1512 }
else if(runnumber<=196000&&runnumber>=195000){
1513 runnumber_offset = 195000;
1514 }
else if(runnumber<=170593&&runnumber>=167902){
1515 runnumber_offset = 167902;
1516 }
else if(runnumber<=246994&&runnumber>=244824){
1517 runnumber_offset = 244824;
1562 AliAnalysisTaskSE::Terminate();
1567 AliError(
"fOutput not available");
1573 AliError(
"fOutputAll not available");
1630 TString normName=
"NormalizationCounter";
1631 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1632 if(cont)normName=(
TString)cont->GetName();
1714 for (
Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1716 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1717 if (!track)
continue;
1728 Int_t nV0s= aodEvent->GetNumberOfV0s();
1729 Int_t nTracks= aodEvent->GetNumberOfTracks();
1731 Bool_t seleTrkFlags[nTracks];
1733 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1735 Bool_t seleV0Flags[nV0s];
1737 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1739 Int_t runnumber_offset = 0;
1740 Int_t runnumber = aodEvent->GetRunNumber();
1741 if(runnumber<=131000&&runnumber>=114000){
1742 runnumber_offset = 114000;
1743 }
else if(runnumber<=196000&&runnumber>=195000){
1744 runnumber_offset = 195000;
1745 }
else if(runnumber<=170593&&runnumber>=167902){
1746 runnumber_offset = 167902;
1747 }
else if(runnumber<=246994&&runnumber>=244824){
1748 runnumber_offset = 244824;
1758 for (
Int_t itrk = 0; itrk<nTracks; itrk++) {
1759 if(!seleTrkFlags[itrk])
continue;
1760 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1765 for (
Int_t iv0 = 0; iv0<nV0s; iv0++) {
1766 if(!seleV0Flags[iv0])
continue;
1767 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1779 for (
Int_t iv0 = 0; iv0<nV0s; iv0++) {
1780 if(!seleV0Flags[iv0])
continue;
1781 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1784 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1785 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1790 for (
Int_t itrk = 0; itrk<nTracks; itrk++) {
1791 if(!seleTrkFlags[itrk])
continue;
1792 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1800 AliAODTrack *trkpid = 0;
1802 trkpid =
fGTI[-trk->GetID()-1];
1807 Int_t cpid = cptrack->GetID();
1808 Int_t cnid = cntrack->GetID();
1809 Int_t lpid = trkpid->GetID();
1810 if((cpid==lpid)||(cnid==lpid))
continue;
1813 if(!secVert)
continue;
1822 elobj->GetSecondaryVtx()->RemoveDaughters();
1824 delete elobj;elobj=NULL;
1848 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1854 fCandidateVariableNames[ 0]=
"Centrality";
1855 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1856 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1857 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1858 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1859 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1860 fCandidateVariableNames[ 6]=
"ElePx";
1861 fCandidateVariableNames[ 7]=
"ElePy";
1862 fCandidateVariableNames[ 8]=
"ElePz";
1863 fCandidateVariableNames[ 9]=
"V0Px";
1864 fCandidateVariableNames[10]=
"V0Py";
1865 fCandidateVariableNames[11]=
"V0Pz";
1866 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1867 fCandidateVariableNames[13]=
"MassLambda";
1868 fCandidateVariableNames[14]=
"MassAntiLambda";
1869 fCandidateVariableNames[15]=
"Eled0";
1870 fCandidateVariableNames[16]=
"V0d0";
1871 fCandidateVariableNames[17]=
"nSigmaTPCele";
1872 fCandidateVariableNames[18]=
"nSigmaTOFele";
1873 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1874 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1875 fCandidateVariableNames[21]=
"EleCharge";
1876 fCandidateVariableNames[22]=
"ProtonPx";
1877 fCandidateVariableNames[23]=
"ProtonPy";
1878 fCandidateVariableNames[24]=
"ProtonPz";
1879 fCandidateVariableNames[25]=
"PiPx";
1880 fCandidateVariableNames[26]=
"PiPy";
1881 fCandidateVariableNames[27]=
"PiPz";
1882 fCandidateVariableNames[28]=
"mcpdglc";
1883 fCandidateVariableNames[29]=
"mclablc";
1884 fCandidateVariableNames[30]=
"mcpdgmomele";
1885 fCandidateVariableNames[31]=
"mcpdgmomv0";
1886 fCandidateVariableNames[32]=
"Mixing";
1887 fCandidateVariableNames[33]=
"mcpdgele";
1888 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1889 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1890 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1891 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1892 fCandidateVariableNames[38]=
"MassK0Short";
1893 fCandidateVariableNames[39]=
"mcpdggrmomele";
1894 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1895 fCandidateVariableNames[41]=
"mcngenele";
1896 fCandidateVariableNames[42]=
"mcngenv0";
1897 fCandidateVariableNames[43]=
"mclcpx";
1898 fCandidateVariableNames[44]=
"mclcpy";
1899 fCandidateVariableNames[45]=
"mclcpz";
1900 fCandidateVariableNames[46]=
"mcelepx";
1901 fCandidateVariableNames[47]=
"mcelepy";
1902 fCandidateVariableNames[48]=
"mcelepz";
1903 fCandidateVariableNames[49]=
"mcv0px";
1904 fCandidateVariableNames[50]=
"mcv0py";
1905 fCandidateVariableNames[51]=
"mcv0pz";
1906 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1907 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1908 fCandidateVariableNames[54]=
"PrimVertx";
1909 fCandidateVariableNames[55]=
"PrimVerty";
1910 fCandidateVariableNames[56]=
"PrimVertz";
1911 fCandidateVariableNames[57]=
"V0Vertx";
1912 fCandidateVariableNames[58]=
"V0Verty";
1913 fCandidateVariableNames[59]=
"V0Vertz";
1915 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1916 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1917 fCandidateVariableNames[62]=
"DcaV0daughters";
1918 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1919 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1920 fCandidateVariableNames[65]=
"MassK0Short2";
1922 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1923 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1925 fCandidateVariableNames[68]=
"EleITSMatch";
1926 fCandidateVariableNames[69]=
"V0PosITSMatch";
1927 fCandidateVariableNames[70]=
"V0NegITSMatch";
1928 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1929 fCandidateVariableNames[72]=
"mcpdgv0";
1931 fCandidateVariableNames[73]=
"mcv0pospx";
1932 fCandidateVariableNames[74]=
"mcv0pospy";
1933 fCandidateVariableNames[75]=
"mcv0pospz";
1934 fCandidateVariableNames[76]=
"mcpdgv0pos";
1935 fCandidateVariableNames[77]=
"mcv0negpx";
1936 fCandidateVariableNames[78]=
"mcv0negpy";
1937 fCandidateVariableNames[79]=
"mcv0negpz";
1938 fCandidateVariableNames[80]=
"mcpdgv0neg";
1940 fCandidateVariableNames[81]=
"ProtonPxPV";
1941 fCandidateVariableNames[82]=
"ProtonPyPV";
1942 fCandidateVariableNames[83]=
"ProtonPzPV";
1943 fCandidateVariableNames[84]=
"PiPxPV";
1944 fCandidateVariableNames[85]=
"PiPyPV";
1945 fCandidateVariableNames[86]=
"PiPzPV";
1947 fCandidateVariableNames[87]=
"mcv0vertx";
1948 fCandidateVariableNames[88]=
"mcv0verty";
1949 fCandidateVariableNames[89]=
"mcv0vertz";
1951 fCandidateVariableNames[90]=
"OnFlyStatus";
1952 fCandidateVariableNames[91]=
"MagneticField";
1953 fCandidateVariableNames[92]=
"EvNumber";
1954 fCandidateVariableNames[93]=
"RunNumber";
1974 for(
Int_t i=0;i<94;i++)
1979 Bool_t anti_lambda_flag = kFALSE;
1982 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1983 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1992 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1993 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1997 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1998 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1999 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
2005 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
2006 Double_t pxv_flip = v0px + lf * uxe;
2007 Double_t pyv_flip = v0py + lf * uye;
2008 Double_t pzv_flip = v0pz + lf * uze;
2009 Double_t pxsum_flip = epx + pxv_flip;
2010 Double_t pysum_flip = epy + pyv_flip;
2011 Double_t pzsum_flip = epz + pzv_flip;
2012 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
2013 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
2021 UInt_t pdgdg[2]={11,3122};
2055 if(anti_lambda_flag){
2073 if(anti_lambda_flag){
2089 AliAODMCParticle *mclc = 0;
2090 AliAODMCParticle *mcele = 0;
2091 AliAODMCParticle *mcv0 = 0;
2093 Int_t mcpdgele_array[100];
2094 Int_t mcpdgv0_array[100];
2095 Int_t mclabelele_array[100];
2096 Int_t mclabelv0_array[100];
2097 Int_t mcngen_ele=-9999;
2098 Int_t mcngen_v0=-9999;
2101 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
2102 if(mclabelele_array[0]>=0){
2103 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
2110 if(mclabelv0_array[0]>=0){
2111 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
2119 Int_t pdgdgv0[2]={211,211};
2120 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
2122 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
2132 Int_t labPos = cptrack->GetLabel();
2134 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
2143 Int_t labNeg = cntrack->GetLabel();
2145 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
2153 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
2196 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
2197 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
2198 if(!anti_lambda_flag){
2207 fVtx1->GetXYZ(posVtx);
2209 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
2219 if(!anti_lambda_flag){
2244 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
2245 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
2249 cont[0] = elobj->InvMass(2,pdgdg);
2250 cont[1] = elobj->Pt();
2255 cont_flip[0] = melam_flip;
2256 cont_flip[1] = ptelam_flip;
2260 cont2[0] = elobj->InvMass(2,pdgdg);
2261 cont2[1] = trk->Pt();
2264 cont_eleptvseta[0] = trk->Pt();
2265 cont_eleptvseta[1] = trk->Eta();
2269 cont_eleptvslambdapt[0] = trk->Pt();
2270 cont_eleptvslambdapt[1] = v0->Pt();
2274 cont_eleptvsd0[0] = trk->Pt();
2275 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
2279 cont_eleptvsv0dl[0] = trk->Pt();
2280 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2284 cont_elelamptvsv0dl[0] = elobj->Pt();
2285 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
2286 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
2289 Double_t cont_elelamptvsv0dl_flip[4];
2290 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
2291 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2292 cont_elelamptvsv0dl_flip[2] = 0.0;
2297 cont_eleptvsv0dca[0] = trk->Pt();
2298 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
2301 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
2305 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
2306 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
2313 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2340 for(
Int_t ih=0;ih<17;ih++){
2342 cont_eleptvscutvars[0] = elobj->Pt();
2346 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2348 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2350 cont_eleptvscutvars[1] = nSigmaTPCele;
2352 cont_eleptvscutvars[1] = nSigmaTOFele;
2354 cont_eleptvscutvars[1] = trk->Eta();
2356 cont_eleptvscutvars[1] = trk->GetITSNcls();
2358 if(!anti_lambda_flag)
2359 cont_eleptvscutvars[1] = v0->MassLambda();
2361 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2364 lPosV0[0] = v0->DecayVertexV0X();
2365 lPosV0[1] = v0->DecayVertexV0Y();
2366 lPosV0[2] = v0->DecayVertexV0Z();
2367 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2369 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2371 if(!anti_lambda_flag)
2372 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2374 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2376 if(!anti_lambda_flag)
2377 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2379 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2381 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2383 cont_eleptvscutvars[1] = v0->MassK0Short();
2385 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2387 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2389 cont_eleptvscutvars[1] = v0->Eta();
2397 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2399 cont_eleptvscutvars[1] = -9999.;
2407 Int_t pdgcode = mclc->GetPdgCode();
2409 cont_mclc[0] = mclc->Pt();
2410 cont_mclc[1] = mclc->Y();
2413 cont_mcele[0] = mcele->Pt();
2414 cont_mcele[1] = mcele->Eta();
2417 Int_t labmotherlc = mclc->GetMother();
2418 Bool_t isbottomfd = kFALSE;
2420 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2421 Int_t pdgmotherlc = motherlc->GetPdgCode();
2422 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){
2426 cont2[1] = mcele->Pt();
2428 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2480 Double_t cont_eleptvslambdaptvslcpt[4];
2481 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2482 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2483 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2484 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2488 cont_allpt[0] = mclc->Pt();
2489 cont_allpt[1] = elobj->Pt();
2490 cont_allpt[2] = trk->Pt();
2491 cont_allpt[3] = v0->Pt();
2500 for(
Int_t ih=0;ih<17;ih++){
2502 cont_eleptvscutvars[0] = elobj->Pt();
2506 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2508 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2510 cont_eleptvscutvars[1] = nSigmaTPCele;
2512 cont_eleptvscutvars[1] = nSigmaTOFele;
2514 cont_eleptvscutvars[1] = trk->Eta();
2516 cont_eleptvscutvars[1] = trk->GetITSNcls();
2518 if(!anti_lambda_flag)
2519 cont_eleptvscutvars[1] = v0->MassLambda();
2521 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2524 lPosV0[0] = v0->DecayVertexV0X();
2525 lPosV0[1] = v0->DecayVertexV0Y();
2526 lPosV0[2] = v0->DecayVertexV0Z();
2527 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2529 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2531 if(!anti_lambda_flag)
2532 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2534 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2536 if(!anti_lambda_flag)
2537 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2539 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2541 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2543 cont_eleptvscutvars[1] = v0->MassK0Short();
2545 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2547 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2549 cont_eleptvscutvars[1] = v0->Eta();
2557 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2559 cont_eleptvscutvars[1] = -9999.;
2566 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2613 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2660 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2690 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2717 for(
Int_t ih=0;ih<17;ih++){
2719 cont_eleptvscutvars[0] = elobj->Pt();
2723 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2725 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2727 cont_eleptvscutvars[1] = nSigmaTPCele;
2729 cont_eleptvscutvars[1] = nSigmaTOFele;
2731 cont_eleptvscutvars[1] = trk->Eta();
2733 cont_eleptvscutvars[1] = trk->GetITSNcls();
2735 if(!anti_lambda_flag)
2736 cont_eleptvscutvars[1] = v0->MassLambda();
2738 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2741 lPosV0[0] = v0->DecayVertexV0X();
2742 lPosV0[1] = v0->DecayVertexV0Y();
2743 lPosV0[2] = v0->DecayVertexV0Z();
2744 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2746 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2748 if(!anti_lambda_flag)
2749 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2751 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2753 if(!anti_lambda_flag)
2754 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2756 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2758 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2760 cont_eleptvscutvars[1] = v0->MassK0Short();
2762 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2764 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2766 cont_eleptvscutvars[1] = v0->Eta();
2774 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2776 cont_eleptvscutvars[1] = -9999.;
2784 Int_t pdgcode = mclc->GetPdgCode();
2786 cont_mclb[0] = mclc->Pt();
2787 cont_mclb[1] = mclc->Y();
2790 cont_mcele[0] = mcele->Pt();
2791 cont_mcele[1] = mcele->Eta();
2794 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2815 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2827 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2839 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2851 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2871 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2873 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2896 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2917 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2941 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2943 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2967 for(
Int_t iv=0;iv<15;iv++){
2972 for(
Int_t iv=0;iv<7;iv++){
2973 cont_cor_nd[iv] = -9999.;
2976 for(
Int_t iv=0;iv<8;iv++){
2977 cont_mass_nd[iv] = -9999.;
2988 if(trk->Charge()>0){
2991 }
else if(trk->Charge()<0){
3001 cont_cor_nd[0] = elobj->Pt();
3003 cont_cor_nd[2] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;;
3004 if(trk->Charge()>0){
3005 if(!anti_lambda_flag) cont_cor_nd[3] = 0;
3006 else cont_cor_nd[3] = 2;
3007 }
else if(trk->Charge()<0){
3008 if(!anti_lambda_flag) cont_cor_nd[3] = 3;
3009 else cont_cor_nd[3] = 1;
3017 Int_t pdgcode = mclc->GetPdgCode();
3018 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
3022 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
3026 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
3030 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
3034 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
3038 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
3042 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
3046 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
3050 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
3099 cont_cor_nd[0] = trk->Pt();
3101 cont_cor_nd[0] = v0->Pt();
3106 cont_cor_nd[0] = trk->Pt();
3108 cont_cor_nd[0] = v0->Pt();
3113 cont_mass_nd[0] = elobj->InvMass(2,pdgdg);
3114 cont_mass_nd[1] = cont_cor_nd[0];
3115 cont_mass_nd[4] = cont_cor_nd[3];
3116 cont_mass_nd[5] = cont_cor_nd[4];
3117 cont_mass_nd[6] = cont_cor_nd[5];
3118 cont_mass_nd[7] = cont_cor_nd[6];
3123 if(melam_flip < 10.&& cosoa < 0.) cont_mass_nd[2]=0;
3126 cont_mass_nd[1] = trk->Pt();
3128 cont_mass_nd[1] = v0->Pt();
3132 cont_mass_nd[1] = trk->Pt();
3134 cont_mass_nd[1] = v0->Pt();
3150 for(
Int_t i=0;i<94;i++)
3159 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
3160 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
3165 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
3166 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
3168 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
3174 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
3179 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
3180 Double_t pxv_flip = pxv + lf * uxe;
3181 Double_t pyv_flip = pyv + lf * uye;
3182 Double_t pzv_flip = pzv + lf * uze;
3183 Double_t pxsum_flip = pxe + pxv_flip;
3184 Double_t pysum_flip = pye + pyv_flip;
3185 Double_t pzsum_flip = pze + pzv_flip;
3186 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
3187 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
3190 UInt_t pdgdg[2]={11,3122};
3243 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
3248 cont_flip[0] = mel_flip;
3249 cont_flip[1] = ptel_flip;
3254 cont2[1] = sqrt(pxe*pxe+pye*pye);
3257 cont_eleptvseta[0] = trke->Pt();
3258 cont_eleptvseta[1] = trke->Eta();
3262 cont_eleptvslambdapt[0] = trke->Pt();
3263 cont_eleptvslambdapt[1] = v0->Pt();
3267 cont_eleptvsd0[0] = trke->Pt();
3268 cont_eleptvsd0[1] = 0.;
3272 cont_eleptvsv0dl[0] = trke->Pt();
3273 cont_eleptvsv0dl[1] = (*v0vars)[15];
3277 cont_eleptvsv0dca[0] = trke->Pt();
3278 cont_eleptvsv0dca[1] = (*v0vars)[16];
3282 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
3283 cont_elelamptvsv0dl[1] = (*v0vars)[15];
3284 cont_elelamptvsv0dl[2] = 0.0;
3287 Double_t cont_elelamptvsv0dl_flip[4];
3288 cont_elelamptvsv0dl_flip[0] = ptel_flip;
3289 cont_elelamptvsv0dl_flip[1] = (*v0vars)[15];
3290 cont_elelamptvsv0dl_flip[2] = 0.0;
3293 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
3294 xyzR125_ele[0] = (*elevars)[0];
3295 xyzR125_ele[1] = (*elevars)[1];
3296 xyzR125_ele[2] = (*elevars)[2];
3297 xyzR125_pr[0] = (*v0vars)[0];
3298 xyzR125_pr[1] = (*v0vars)[1];
3299 xyzR125_pr[2] = (*v0vars)[2];
3300 xyzR125_pi[0] = (*v0vars)[3];
3301 xyzR125_pi[1] = (*v0vars)[4];
3302 xyzR125_pi[2] = (*v0vars)[5];
3305 rdhfcutvars[0] = xyzR125_ele[0];
3306 rdhfcutvars[1] = xyzR125_ele[1];
3307 rdhfcutvars[2] = xyzR125_ele[2];
3308 rdhfcutvars[3] = xyzR125_pr[0];
3309 rdhfcutvars[4] = xyzR125_pr[1];
3310 rdhfcutvars[5] = xyzR125_pr[2];
3311 rdhfcutvars[6] = xyzR125_pi[0];
3312 rdhfcutvars[7] = xyzR125_pi[1];
3313 rdhfcutvars[8] = xyzR125_pi[2];
3324 if(((
int)trke->T())*chargepr>0){
3384 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
3386 if(((
int)trke->T())*chargepr>0){
3410 for(
Int_t iv=0;iv<15;iv++){
3415 for(
Int_t iv=0;iv<7;iv++){
3416 cont_cor_nd[iv] = -9999.;
3419 for(
Int_t iv=0;iv<8;iv++){
3420 cont_mass_nd[iv] = -9999.;
3427 fCorrelationVariables[4] = (*v0vars)[15];
3428 fCorrelationVariables[5] = (*elevars)[5];
3429 fCorrelationVariables[6] = 2;
3431 if(chargepr>0) fCorrelationVariables[7] = 0;
3432 else fCorrelationVariables[7] = 2;
3433 }
else if(trke->T()<0){
3434 if(chargepr>0) fCorrelationVariables[7] = 1;
3435 else fCorrelationVariables[7] = 3;
3437 fCorrelationVariables[8] = (*elevars)[6];
3438 fCorrelationVariables[9] = (*elevars)[7];
3440 fCorrelationVariables[11] = sqrt(pxsum*pxsum+pysum*pysum);
3441 fCorrelationVariables[12] = mel;
3442 fCorrelationVariables[13] = cosoa;
3444 cont_cor_nd[0] = sqrt(pxsum*pxsum+pysum*pysum);
3445 cont_cor_nd[1] = TVector2::Phi_mpi_pi(v0->Phi()-trke->Phi());
3446 cont_cor_nd[2] = (*v0vars)[15];
3448 if(chargepr>0) cont_cor_nd[3] = 0;
3449 else cont_cor_nd[3] = 2;
3450 }
else if(trke->T()<0){
3451 if(chargepr>0) cont_cor_nd[3] = 3;
3452 else cont_cor_nd[3] = 1;
3454 cont_cor_nd[4] = fCorrelationVariables[8];
3456 if(fabs(fCorrelationVariables[9]-1013)<0.001) cont_cor_nd[5] = 7;
3457 if(fabs(fCorrelationVariables[9]-1016)<0.001) cont_cor_nd[5] = 8;
3458 if(fabs(fCorrelationVariables[9]-1019)<0.001) cont_cor_nd[5] = 9;
3466 if(mel<fAnalCuts->GetEleLambdaMassMax())
3468 cont_cor_nd[0] = trke->Pt();
3470 cont_cor_nd[0] = v0->Pt();
3474 cont_mass_nd[0] = mel;
3475 cont_mass_nd[1] = cont_cor_nd[0];
3476 cont_mass_nd[4] = cont_cor_nd[3];
3477 cont_mass_nd[5] = cont_cor_nd[4];
3478 cont_mass_nd[6] = cont_cor_nd[5];
3479 cont_mass_nd[7] = cont_cor_nd[6];
3484 if(mel_flip < 10.&& cosoa < 0.) cont_mass_nd[2]=0;
3486 cont_mass_nd[0] = trke->Pt();
3488 cont_mass_nd[0] = v0->Pt();
3500 const char* nameoutput = GetOutputSlot(11)->GetContainer()->GetName();
3506 fCandidateVariableNames[ 0]=
"Px";
3507 fCandidateVariableNames[ 1]=
"Py";
3508 fCandidateVariableNames[ 2]=
"Pz";
3509 fCandidateVariableNames[ 3]=
"Charge";
3510 fCandidateVariableNames[ 4]=
"LambdaMass";
3511 fCandidateVariableNames[ 5]=
"Bz";
3512 fCandidateVariableNames[ 6]=
"Centrality";
3513 fCandidateVariableNames[ 7]=
"PrimVertZ";
3514 fCandidateVariableNames[ 8]=
"EvNumber";
3515 fCandidateVariableNames[ 9]=
"RunNumber";
3516 fCandidateVariableNames[10]=
"EventPlane";
3517 fCandidateVariableNames[11]=
"Vars0";
3518 fCandidateVariableNames[12]=
"Vars1";
3519 fCandidateVariableNames[13]=
"Vars2";
3520 fCandidateVariableNames[14]=
"Vars3";
3521 fCandidateVariableNames[15]=
"Vars4";
3522 fCandidateVariableNames[16]=
"MCVars0";
3523 fCandidateVariableNames[17]=
"MCVars1";
3524 fCandidateVariableNames[18]=
"MCVars2";
3525 fCandidateVariableNames[19]=
"MCVars3";
3526 fCandidateVariableNames[20]=
"MCVars4";
3527 fCandidateVariableNames[21]=
"Vars5";
3528 fCandidateVariableNames[22]=
"Vars6";
3529 fCandidateVariableNames[23]=
"Vars7";
3530 fCandidateVariableNames[24]=
"Vars8";
3531 fCandidateVariableNames[25]=
"Vars9";
3532 fCandidateVariableNames[26]=
"Vars10";
3533 fCandidateVariableNames[27]=
"Vars11";
3534 fCandidateVariableNames[28]=
"Vars12";
3535 fCandidateVariableNames[29]=
"Vars13";
3554 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
3560 fCandidateVariableNames[ 0]=
"ElePx";
3561 fCandidateVariableNames[ 1]=
"ElePy";
3562 fCandidateVariableNames[ 2]=
"ElePz";
3563 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
3564 fCandidateVariableNames[ 4]=
"ITSNcls";
3565 fCandidateVariableNames[ 5]=
"TPCNcls";
3566 fCandidateVariableNames[ 6]=
"TPCNclsPID";
3567 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
3568 fCandidateVariableNames[ 8]=
"d0R";
3569 fCandidateVariableNames[ 9]=
"d0Z";
3570 fCandidateVariableNames[10]=
"ITSClusterMap";
3571 fCandidateVariableNames[11]=
"nSigmaTPCele";
3572 fCandidateVariableNames[12]=
"nSigmaTOFele";
3573 fCandidateVariableNames[13]=
"nSigmaTPCpi";
3574 fCandidateVariableNames[14]=
"nSigmaTPCka";
3575 fCandidateVariableNames[15]=
"nSigmaTPCpr";
3576 fCandidateVariableNames[16]=
"EvNumber";
3577 fCandidateVariableNames[17]=
"EleCharge";
3578 fCandidateVariableNames[18]=
"ElePdgCode";
3579 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3580 fCandidateVariableNames[20]=
"mcelepx";
3581 fCandidateVariableNames[21]=
"mcelepy";
3582 fCandidateVariableNames[22]=
"mcelepz";
3583 fCandidateVariableNames[23]=
"Centrality";
3584 fCandidateVariableNames[24]=
"PrimVertZ";
3585 fCandidateVariableNames[25]=
"RunNumber";
3602 AliAODTrack *trkpid = 0;
3604 trkpid =
fGTI[-trk->GetID()-1];
3614 for(
Int_t i=0;i<6;i++){
3615 if(trk->HasPointOnITSLayer(i)) ncls_its++;
3617 Int_t ncls_its_shared = 0;
3618 for(
Int_t i=0;i<6;i++){
3619 if(trk->HasPointOnITSLayer(i) && trk->HasSharedPointOnITSLayer(i)) ncls_its_shared++;
3631 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3640 Int_t pdgEle = -9999;
3641 Int_t pdgEleMother = -9999;
3645 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
3646 for(
Int_t i=0;i<100;i++){
3647 pdgarray_ele[i]=-9999;
3648 labelarray_ele[i]=-9999;
3651 Int_t mcetype = -9999;
3654 Int_t labEle = trk->GetLabel();
3655 if(labEle<0)
return;
3656 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3658 if(abs(mcetrk->GetPdgCode())!=11)
return;
3660 pdgEle = mcetrk->GetPdgCode();
3661 mcelepx = mcetrk->Px();
3662 mcelepy = mcetrk->Py();
3663 mcelepz = mcetrk->Pz();
3665 Int_t labEleMother = mcetrk->GetMother();
3666 if(labEleMother>-1){
3667 AliAODMCParticle *mcemothertrk = (AliAODMCParticle*)mcArray->At(labEleMother);
3669 pdgEleMother = mcemothertrk->GetPdgCode();
3676 Bool_t gamma_flag = kFALSE;
3677 Bool_t pi0_flag = kFALSE;
3678 Bool_t eta_flag = kFALSE;
3681 if(abs(pdgarray_ele[0])>400&&abs(pdgarray_ele[0])<440){
3684 if(abs(pdgarray_ele[0])>4000&&abs(pdgarray_ele[0])<4400){
3687 if(abs(pdgarray_ele[0])>500&&abs(pdgarray_ele[0])<540){
3690 if(abs(pdgarray_ele[0])>5000&&abs(pdgarray_ele[0])<5400){
3694 if(abs(pdgarray_ele[0])==22){
3701 if((abs(pdgarray_ele[0])==22) && (abs(pdgarray_ele[1])==111)){
3703 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[1]);
3704 pt_pi0 = mctrkm->Pt();
3706 if(abs(pdgarray_ele[0])==111){
3708 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[0]);
3709 pt_pi0 = mctrkm->Pt();
3711 if((abs(pdgarray_ele[0])==22) && (abs(pdgarray_ele[1])==221)){
3713 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[1]);
3714 pt_eta = mctrkm->Pt();
3716 if(abs(pdgarray_ele[0])==221){
3718 AliAODMCParticle *mctrkm = (AliAODMCParticle*)mcArray->At(labelarray_ele[0]);
3719 pt_eta = mctrkm->Pt();
3724 cont_pi0[0] = trk->Pt();
3725 cont_pi0[1] = pt_pi0;
3733 cont_eta[0] = trk->Pt();
3734 cont_eta[1] = pt_eta;
3740 if(hfe_flag==0)
return;
3755 pv[0] =
fVtx1->GetX();
3756 pv[1] =
fVtx1->GetY();
3757 pv[2] =
fVtx1->GetZ();
3758 Double_t xyzR125[3]={9999.,9999.,9999.};
3760 TVector *varvec =
new TVector(8);
3761 (*varvec)[0] = xyzR125[0];
3762 (*varvec)[1] = xyzR125[1];
3763 (*varvec)[2] = xyzR125[2];
3766 (*varvec)[5] = d0z0[0];
3767 (*varvec)[6] = (
Int_t)isconv + 2*(
Int_t)isconv_like;
3768 (*varvec)[7] = mcetype;
3771 m_ReservoirE[
fPoolIndex][nextRes].push_back(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3790 for(
Int_t i=0;i<26;i++){
3793 for(
Int_t i=0;i<30;i++){
3900 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3906 fCandidateVariableNames[ 0]=
"V0Px";
3907 fCandidateVariableNames[ 1]=
"V0Py";
3908 fCandidateVariableNames[ 2]=
"V0Pz";
3909 fCandidateVariableNames[ 3]=
"MassLambda";
3910 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3911 fCandidateVariableNames[ 5]=
"ProtonPx";
3912 fCandidateVariableNames[ 6]=
"ProtonPy";
3913 fCandidateVariableNames[ 7]=
"ProtonPz";
3914 fCandidateVariableNames[ 8]=
"PionPx";
3915 fCandidateVariableNames[ 9]=
"PionPy";
3916 fCandidateVariableNames[10]=
"PionPz";
3917 fCandidateVariableNames[11]=
"RfidV0";
3918 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3919 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3920 fCandidateVariableNames[14]=
"DcaV0daughters";
3921 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3922 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3923 fCandidateVariableNames[17]=
"MassK0Short";
3924 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3925 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3926 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3927 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3928 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3929 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3930 fCandidateVariableNames[24]=
"V0PdgCode";
3931 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3932 fCandidateVariableNames[26]=
"mcv0px";
3933 fCandidateVariableNames[27]=
"mcv0py";
3934 fCandidateVariableNames[28]=
"mcv0pz";
3935 fCandidateVariableNames[29]=
"mcv0vertx";
3936 fCandidateVariableNames[30]=
"mcv0verty";
3937 fCandidateVariableNames[31]=
"mcv0vertz";
3938 fCandidateVariableNames[32]=
"EvNumber";
3939 fCandidateVariableNames[33]=
"Centrality";
3940 fCandidateVariableNames[34]=
"PrimVertZ";
3941 fCandidateVariableNames[35]=
"RunNumber";
3957 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3958 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3959 if(!cptrack)
return;
3960 if(!cntrack)
return;
3979 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3981 fVtx1->GetXYZ(posVtx);
3982 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3983 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3985 cont_arm[0] = v0->Pt();
3986 cont_arm[1] = v0->PtArmV0();
3987 cont_arm[2] = v0->AlphaV0();
4008 Double_t phi_alice = atan2(momv0y,momv0x);
4009 if(phi_alice<0.) phi_alice += 2 * M_PI;
4015 fVtx1->GetXYZ(posVtx);
4016 TLorentzVector *lv =
new TLorentzVector();
4017 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
4018 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
4019 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
4021 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
4025 TVector *varvec =
new TVector(17);
4026 (*varvec)[0] = xyzR125pr[0];
4027 (*varvec)[1] = xyzR125pr[1];
4028 (*varvec)[2] = xyzR125pr[2];
4029 (*varvec)[3] = xyzR125pi[0];
4030 (*varvec)[4] = xyzR125pi[1];
4031 (*varvec)[5] = xyzR125pi[2];
4032 (*varvec)[6] = v0->DecayVertexV0X();
4033 (*varvec)[7] = v0->DecayVertexV0Y();
4034 (*varvec)[8] = v0->DecayVertexV0Z();
4035 (*varvec)[9] = v0->MomPosX();
4036 (*varvec)[10] = v0->MomPosY();
4037 (*varvec)[11] = v0->MomPosZ();
4038 (*varvec)[12] = v0->MomNegX();
4039 (*varvec)[13] = v0->MomNegY();
4040 (*varvec)[14] = v0->MomNegZ();
4041 (*varvec)[15] = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
4042 (*varvec)[16] = v0->DcaV0ToPrimVertex();
4045 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
4049 TVector *varvec =
new TVector(17);
4050 (*varvec)[0] = xyzR125pr[0];
4051 (*varvec)[1] = xyzR125pr[1];
4052 (*varvec)[2] = xyzR125pr[2];
4053 (*varvec)[3] = xyzR125pi[0];
4054 (*varvec)[4] = xyzR125pi[1];
4055 (*varvec)[5] = xyzR125pi[2];
4056 (*varvec)[6] = v0->DecayVertexV0X();
4057 (*varvec)[7] = v0->DecayVertexV0Y();
4058 (*varvec)[8] = v0->DecayVertexV0Z();
4059 (*varvec)[9] = v0->MomNegX();
4060 (*varvec)[10] = v0->MomNegY();
4061 (*varvec)[11] = v0->MomNegZ();
4062 (*varvec)[12] = v0->MomPosX();
4063 (*varvec)[13] = v0->MomPosY();
4064 (*varvec)[14] = v0->MomPosZ();
4065 (*varvec)[15] = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
4066 (*varvec)[16] = v0->DcaV0ToPrimVertex();
4073 Int_t v0pdgcode = -9999;
4074 Int_t v0motherpdgcode = -9999;
4081 Int_t pdgarray_v0[100], labelarray_v0[100], ngen_v0;
4084 Int_t pdgdgv0[2]={2212,211};
4085 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
4098 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
4100 mcv0px = mcv0trk->Px();
4101 mcv0py = mcv0trk->Py();
4102 mcv0pz = mcv0trk->Pz();
4104 Int_t labptrk = cptrack->GetLabel();
4106 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
4108 mcv0vertx = mcv0ptrk->Xv();
4109 mcv0verty = mcv0ptrk->Yv();
4110 mcv0vertz = mcv0ptrk->Zv();
4116 Bool_t hfv0_flag = kFALSE;
4117 Bool_t promptv0_flag = kFALSE;
4118 v0pdgcode = mcv0trk->GetPdgCode();
4119 Int_t labv0mother = mcv0trk->GetMother();
4121 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
4123 v0motherpdgcode = motherv0->GetPdgCode();
4124 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
4127 if(abs(v0motherpdgcode)==3322){
4129 }
else if(abs(v0motherpdgcode)==3312){
4131 }
else if(abs(v0motherpdgcode)==3334){
4135 promptv0_flag = kTRUE;
4143 for(
Int_t i=0;i<100;i++){
4144 pdgarray_v0[i] = -9999;
4145 labelarray_v0[i] = -9999;
4152 for(
Int_t i=0;i<36;i++){
4156 for(
Int_t i=0;i<30;i++){
4160 Bool_t isparticle = kFALSE;
4334 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
4340 fCandidateVariableNames[ 0]=
"Centrality";
4341 fCandidateVariableNames[ 1]=
"DecayType";
4342 fCandidateVariableNames[ 2]=
"LcPx";
4343 fCandidateVariableNames[ 3]=
"LcPy";
4344 fCandidateVariableNames[ 4]=
"LcPz";
4345 fCandidateVariableNames[ 5]=
"ElePx";
4346 fCandidateVariableNames[ 6]=
"ElePy";
4347 fCandidateVariableNames[ 7]=
"ElePz";
4348 fCandidateVariableNames[ 8]=
"V0Px";
4349 fCandidateVariableNames[ 9]=
"V0Py";
4350 fCandidateVariableNames[10]=
"V0Pz";
4364 if(!mcepart)
return;
4365 if(!mcv0part)
return;
4367 for(
Int_t i=0;i<11;i++){
4386 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
4390 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
4392 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
4395 cont[0] = InvMassEleLambda;
4396 cont[1] = mcpart->Pt();
4399 cont2[0] = InvMassEleLambda;
4400 cont2[1] = mcepart->Pt();
4403 cont_eleptvseta[0] = mcepart->Pt();
4404 cont_eleptvseta[1] = mcepart->Eta();
4407 cont_eleptvslambdapt[0] = mcepart->Pt();
4408 cont_eleptvslambdapt[1] = mcv0part->Pt();
4410 Double_t cont_eleptvslambdaptvslcpt[4];
4411 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
4412 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
4413 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
4416 contmc[0] = mcpart->Pt();
4417 contmc[1] = mcpart->Y();
4420 contmcele[0] = mcepart->Pt();
4421 contmcele[1] = mcepart->Eta();
4424 contmc_withele[0] = mcpart->Pt();
4425 contmc_withele[1] = mcpart->Y();
4426 contmc_withele[2] = mcepart->Y();
4430 esdcuts->GetEtaRange(etamin,etamax);
4437 if(fabs(mcepart->Eta())<0.8){
4444 if(fabs(mcepart->Eta())<etamax){
4448 if(InvMassEleLambda<fAnalCuts->GetEleLambdaMassMax()){
4454 if(fabs(mcpart->Y())<0.7){
4455 if(InvMassEleLambda<fAnalCuts->GetEleLambdaMassMax()){
4460 }
else if(decaytype==1){
4468 if(fabs(mcepart->Eta())<etamax){
4472 if(InvMassEleLambda<fAnalCuts->GetEleLambdaMassMax()){
4478 }
else if(decaytype==2){
4486 if(fabs(mcepart->Eta())<etamax){
4490 if(InvMassEleLambda<fAnalCuts->GetEleLambdaMassMax()){
4496 }
else if(decaytype==3){
4500 }
else if(decaytype==10){
4509 if(fabs(mcepart->Eta())<0.8){
4512 }
else if(decaytype==11){
4528 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
4534 fCandidateVariableNames[ 0]=
"Centrality";
4535 fCandidateVariableNames[ 1]=
"ElePx";
4536 fCandidateVariableNames[ 2]=
"ElePy";
4537 fCandidateVariableNames[ 3]=
"ElePz";
4538 fCandidateVariableNames[ 4]=
"ElePdgCode";
4539 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
4540 fCandidateVariableNames[ 6]=
"RunNumber";
4541 fCandidateVariableNames[ 7]=
"EvNumber";
4554 if(!mcepart)
return;
4556 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
4561 Int_t semi_flag = 0;
4562 if(pdgarray_ele[0]>400&&pdgarray_ele[0]<450) semi_flag = 1;
4563 if(pdgarray_ele[0]>4000&&pdgarray_ele[0]<4500) semi_flag = 1;
4564 if(pdgarray_ele[0]>500&&pdgarray_ele[0]<550) semi_flag = 2;
4565 if(pdgarray_ele[0]>5000&&pdgarray_ele[0]<5500) semi_flag = 2;
4568 contmc[0] = mcepart->Pt();
4569 contmc[1] = mcepart->Eta();
4572 if(semi_flag==1 && !ele_from_bottom){
4575 if(semi_flag==1 && ele_from_bottom){
4582 Bool_t hfe_flag = kFALSE;
4583 Int_t labemother = mcepart->GetMother();
4584 Int_t pdgmotherele = -9999;
4586 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
4587 pdgmotherele = motherele->GetPdgCode();
4588 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
4593 if(!hfe_flag)
return;
4597 for(
Int_t i=0;i<8;i++){
4621 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
4627 fCandidateVariableNames[ 0]=
"Centrality";
4628 fCandidateVariableNames[ 1]=
"V0Px";
4629 fCandidateVariableNames[ 2]=
"V0Py";
4630 fCandidateVariableNames[ 3]=
"V0Pz";
4631 fCandidateVariableNames[ 4]=
"V0PdgCode";
4632 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
4633 fCandidateVariableNames[ 6]=
"RunNumber";
4634 fCandidateVariableNames[ 7]=
"EvNumber";
4647 if(!mcv0part)
return;
4649 for(
Int_t i=0;i<8;i++){
4653 Bool_t hfv0_flag = kFALSE;
4654 Int_t labv0mother = mcv0part->GetMother();
4655 Int_t pdgmotherv0 = -9999;
4657 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
4659 pdgmotherv0 = motherv0->GetPdgCode();
4660 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
4665 if(!hfv0_flag)
return;
4668 contmc[0] = mcv0part->Pt();
4669 contmc[1] = mcv0part->Eta();
4693 fCEvents =
new TH1F(
"fCEvents",
"conter",20,-0.5,19.5);
4695 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
4696 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
4697 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
4698 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
4699 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
4700 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
4701 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
4702 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
4703 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
4704 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
4705 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
4706 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
4707 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
4708 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
4710 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
4711 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
4713 fCEvents->GetXaxis()->SetBinLabel(19,
"No RP");
4714 fCEvents->GetXaxis()->SetBinLabel(20,
"No Qvector (only for TPC event plane)");
4716 fCEvents->GetYaxis()->SetTitle(
"counts");
4718 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
4720 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
4721 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
4722 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
4723 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
4724 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
4725 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
4727 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
4728 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
4729 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
4730 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
4731 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
4733 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
4734 fHEventPlane =
new TH1F(
"fHEventPlane",
"conter",100,-3.14,3.14);
4735 fHNTrackletvsZ =
new TH2F(
"fHNTrackletvsZ",
"N_{tracklet} vs z",30,-15.,15.,120,-0.5,119.5);
4761 Int_t bins_base[3]= {23,100 ,10};
4762 Double_t xmin_base[3]={1.1,0 ,0.00};
4763 Double_t xmax_base[3]={5.7,20. ,100};
4764 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
4840 Int_t bins_base_elept[3]= {10,100 ,10};
4841 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
4842 Double_t xmax_base_elept[3]={3.1,10. ,100};
4914 Int_t bins_eleptvseta[3]= {10,10 ,10};
4915 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4916 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4918 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4920 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4922 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4924 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4926 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4928 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4931 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4932 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4933 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4935 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4937 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4939 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4941 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4943 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4945 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4947 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4950 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4951 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4952 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4953 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4955 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4958 Int_t bins_allpt[4]= {10,10,20,20};
4959 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4960 Double_t xmax_allpt[4]={20.,20.,10.,10};
4971 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4972 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4973 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4975 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4977 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4979 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4981 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4983 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4985 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);