37 #include <TParticle.h>
38 #include <TParticlePDG.h>
42 #include <THnSparse.h>
43 #include <TLorentzVector.h>
46 #include <TDatabasePDG.h>
47 #include <AliAnalysisDataSlot.h>
48 #include <AliAnalysisDataContainer.h>
50 #include "AliMCEvent.h"
51 #include "AliAnalysisManager.h"
52 #include "AliAODMCHeader.h"
53 #include "AliAODHandler.h"
55 #include "AliExternalTrackParam.h"
56 #include "AliAODVertex.h"
57 #include "AliESDVertex.h"
58 #include "AliAODRecoDecay.h"
61 #include "AliESDtrack.h"
62 #include "AliAODTrack.h"
64 #include "AliAODcascade.h"
65 #include "AliAODMCParticle.h"
66 #include "AliAnalysisTaskSE.h"
68 #include "AliPIDResponse.h"
69 #include "AliPIDCombined.h"
70 #include "AliTOFPIDResponse.h"
72 #include "AliInputEventHandler.h"
73 #include "AliESDtrackCuts.h"
74 #include "AliNeutralTrackParam.h"
75 #include "AliKFParticle.h"
76 #include "AliKFVertex.h"
77 #include "AliExternalTrackParam.h"
78 #include "AliESDtrack.h"
79 #include "AliCentrality.h"
80 #include "AliVertexerTracks.h"
81 #include "AliEventPoolManager.h"
102 fIsEventSelected(kFALSE),
103 fWriteVariableTree(kFALSE),
104 fWriteEachVariableTree(kFALSE),
105 fWriteMCVariableTree(kFALSE),
107 fEleVariablesTree(0),
110 fMCEleVariablesTree(0),
111 fMCV0VariablesTree(0),
112 fMCGenPairVariablesTree(0),
113 fCorrelationVariablesTree(0),
114 fReconstructPrimVert(kFALSE),
120 fCandidateVariables(),
121 fCandidateEleVariables(),
122 fCandidateV0Variables(),
123 fCandidateMCVariables(),
124 fCandidateMCEleVariables(),
125 fCandidateMCV0Variables(),
126 fCandidateMCGenPairVariables(),
127 fCorrelationVariables(),
135 fUseCentralityV0M(kFALSE),
138 fMCDoPairAnalysis(kFALSE),
139 fHistoEleLambdaMass(0),
140 fHistoEleLambdaMassRS(0),
141 fHistoEleLambdaMassRS1(0),
142 fHistoEleLambdaMassRS2(0),
143 fHistoEleLambdaMassWS(0),
144 fHistoEleLambdaMassWS1(0),
145 fHistoEleLambdaMassWS2(0),
146 fHistoEleLambdaMassRSMix(0),
147 fHistoEleLambdaMassRSMix1(0),
148 fHistoEleLambdaMassRSMix2(0),
149 fHistoEleLambdaMassWSMix(0),
150 fHistoEleLambdaMassWSMix1(0),
151 fHistoEleLambdaMassWSMix2(0),
152 fHistoEleLambdaMassRSSide(0),
153 fHistoEleLambdaMassRSSide1(0),
154 fHistoEleLambdaMassRSSide2(0),
155 fHistoEleLambdaMassWSSide(0),
156 fHistoEleLambdaMassWSSide1(0),
157 fHistoEleLambdaMassWSSide2(0),
158 fHistoEleLambdaMassRSAway(0),
159 fHistoEleLambdaMassRS1Away(0),
160 fHistoEleLambdaMassRS2Away(0),
161 fHistoEleLambdaMassWSAway(0),
162 fHistoEleLambdaMassWS1Away(0),
163 fHistoEleLambdaMassWS2Away(0),
164 fHistoEleLambdaMassRSMixAway(0),
165 fHistoEleLambdaMassRSMix1Away(0),
166 fHistoEleLambdaMassRSMix2Away(0),
167 fHistoEleLambdaMassWSMixAway(0),
168 fHistoEleLambdaMassWSMix1Away(0),
169 fHistoEleLambdaMassWSMix2Away(0),
170 fHistoEleLambdaMassRSSideAway(0),
171 fHistoEleLambdaMassRSSide1Away(0),
172 fHistoEleLambdaMassRSSide2Away(0),
173 fHistoEleLambdaMassWSSideAway(0),
174 fHistoEleLambdaMassWSSide1Away(0),
175 fHistoEleLambdaMassWSSide2Away(0),
176 fHistoEleLambdaMassvsElePtRS(0),
177 fHistoEleLambdaMassvsElePtWS(0),
178 fHistoEleLambdaMassvsElePtRSMix(0),
179 fHistoEleLambdaMassvsElePtWSMix(0),
180 fHistoEleLambdaMassvsElePtRSSide(0),
181 fHistoEleLambdaMassvsElePtWSSide(0),
182 fHistoEleLambdaMassvsElePtRS1(0),
183 fHistoEleLambdaMassvsElePtWS1(0),
184 fHistoEleLambdaMassvsElePtRSMix1(0),
185 fHistoEleLambdaMassvsElePtWSMix1(0),
186 fHistoEleLambdaMassvsElePtRSSide1(0),
187 fHistoEleLambdaMassvsElePtWSSide1(0),
188 fHistoEleLambdaMassvsElePtRS2(0),
189 fHistoEleLambdaMassvsElePtWS2(0),
190 fHistoEleLambdaMassvsElePtRSMix2(0),
191 fHistoEleLambdaMassvsElePtWSMix2(0),
192 fHistoEleLambdaMassvsElePtRSSide2(0),
193 fHistoEleLambdaMassvsElePtWSSide2(0),
198 fHistoEleLambdaMassMCS(0),
199 fHistoEleLambdaMassLambdabMCS(0),
200 fHistoEleLambdaMassXibMCS(0),
201 fHistoEleLambdaMassMCGen(0),
202 fHistoEleLambdaMassvsElePtMCS(0),
203 fHistoEleLambdaMassvsElePtMCGen(0),
204 fHistoEleLambdaMassvsElePtMCS1(0),
205 fHistoEleLambdaMassvsElePtMCGen1(0),
206 fHistoEleLambdaMassvsElePtMCS2(0),
207 fHistoEleLambdaMassvsElePtMCGen2(0),
210 fHistoElePtvsEtaRS(0),
211 fHistoElePtvsEtaWS(0),
212 fHistoElePtvsEtaRSMix(0),
213 fHistoElePtvsEtaWSMix(0),
214 fHistoElePtvsEtaMCS(0),
215 fHistoElePtvsEtaMCGen(0),
216 fHistoElePtvsLambdaPtRS(0),
217 fHistoElePtvsLambdaPtWS(0),
218 fHistoElePtvsLambdaPtRSMix(0),
219 fHistoElePtvsLambdaPtWSMix(0),
220 fHistoElePtvsLambdaPtMCS(0),
221 fHistoElePtvsLambdaPtvsLcPtMCS(0),
222 fHistoElePtvsLambdaPtMCGen(0),
223 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
224 fHistoElePtvsLambdaPtMCLcGen(0),
225 fHistoElePtvsd0RS(0),
226 fHistoElePtvsd0WS(0),
227 fHistoElePtvsd0RSMix(0),
228 fHistoElePtvsd0WSMix(0),
229 fHistoElePtvsd0MCS(0),
230 fHistoElePtvsd0PromptMCS(0),
231 fHistoElePtvsd0BFeeddownMCS(0),
232 fHistoEleLambdaMassFeeddownXic0MCS(0),
233 fHistoEleLambdaMassFeeddownXic0MCGen(0),
234 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
235 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
236 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
237 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
238 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
239 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
240 fHistoElePtFeeddownXic0MCS(0),
241 fHistoElePtFeeddownXic0MCGen(0),
242 fHistoElePtvsEtaFeeddownXic0MCS(0),
243 fHistoElePtvsEtaFeeddownXic0MCGen(0),
244 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
245 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
246 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
247 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
248 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
249 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
250 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
251 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
252 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
253 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
254 fHistoElePtFeeddownXicPlusMCS(0),
255 fHistoElePtFeeddownXicPlusMCGen(0),
256 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
257 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
258 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
259 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
262 fHistoBachPtMCGen(0),
264 fHistoLambdaMassvsPt(0),
265 fHistoLambdaMassvsPtMCS(0),
266 fHistoLambdaMassvsPtMCGen(0),
267 fHistoLambdaPtvsEta(0),
268 fHistoLambdaPtvsEtaMCS(0),
269 fHistoLambdaPtvsEtaMCGen(0),
270 fHistoLambdaArmenteros(0),
271 fHistoLambdaArmenterosMCS(0),
272 fHistoLambdaPtvsDl(0),
273 fHistoLambdaPtvsDlSide(0),
274 fHistoLambdaPtvsDlMCS(0),
275 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
276 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
277 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
278 fHistoK0sMassvsPt(0),
279 fHistoElectronTPCPID(0),
280 fHistoElectronTOFPID(0),
281 fHistoElectronTPCSelPID(0),
282 fHistoElectronTOFSelPID(0),
283 fHistoElectronTPCPIDSelTOF(0),
284 fHistoElectronTPCPIDSelTOFSmallEta(0),
285 fHistoElectronTPCPIDSelTOFLargeEta(0),
286 fHistoMassConversionsMin(0),
287 fHistoMassConversionsSameSignMin(0),
288 fHistoElectronQovPtvsPhi(0),
289 fHistoLambdaQovPtvsPhi(0),
296 fHistoLcMCGenWithEleA(0),
297 fHistoLcMCGenWithEleB(0),
298 fHistoLbMCGenWithEleA(0),
299 fHistoLbMCGenWithEleB(0),
300 fHistoLbMCGenWithLambdac(0),
307 fHistoFeedDownXic0MCGen(0),
308 fHistoFeedDownXic0MCGen1(0),
309 fHistoFeedDownXic0MCGen2(0),
310 fHistoFeedDownXic0MCS(0),
311 fHistoFeedDownXic0MCS1(0),
312 fHistoFeedDownXic0MCS2(0),
313 fHistoFeedDownXicPlusMCGen(0),
314 fHistoFeedDownXicPlusMCGen1(0),
315 fHistoFeedDownXicPlusMCGen2(0),
316 fHistoFeedDownXicPlusMCS(0),
317 fHistoFeedDownXicPlusMCS1(0),
318 fHistoFeedDownXicPlusMCS2(0),
319 fHistoFeedDownSigmaMCGen(0),
320 fHistoFeedDownSigmaMCGen1(0),
321 fHistoFeedDownSigmaMCGen2(0),
322 fHistoFeedDownSigmaMCS(0),
323 fHistoFeedDownSigmaMCS1(0),
324 fHistoFeedDownSigmaMCS2(0),
325 fHistoLcElectronMCGen(0),
326 fHistoLcElectronMCGen1(0),
327 fHistoLcElectronMCGen2(0),
328 fHistoLcElectronMCS(0),
329 fHistoLcElectronMCS1(0),
330 fHistoLcElectronMCS2(0),
331 fHistoLbElectronMCGen(0),
332 fHistoLbElectronMCGen1(0),
333 fHistoLbElectronMCGen2(0),
334 fHistoLbElectronMCS(0),
335 fHistoLbElectronMCS1(0),
336 fHistoLbElectronMCS2(0),
337 fHistoElectronFeedDownXic0MCGen(0),
338 fHistoElectronFeedDownXic0MCGen1(0),
339 fHistoElectronFeedDownXic0MCGen2(0),
340 fHistoElectronFeedDownXic0MCS(0),
341 fHistoElectronFeedDownXic0MCS1(0),
342 fHistoElectronFeedDownXic0MCS2(0),
343 fHistoElectronFeedDownXicPlusMCGen(0),
344 fHistoElectronFeedDownXicPlusMCGen1(0),
345 fHistoElectronFeedDownXicPlusMCGen2(0),
346 fHistoElectronFeedDownXicPlusMCS(0),
347 fHistoElectronFeedDownXicPlusMCS1(0),
348 fHistoElectronFeedDownXicPlusMCS2(0),
349 fHistoElectronFeedDownXibMCGen(0),
350 fHistoElectronFeedDownXibMCGen1(0),
351 fHistoElectronFeedDownXibMCGen2(0),
352 fHistoElectronFeedDownXibMCS(0),
353 fHistoElectronFeedDownXibMCS1(0),
354 fHistoElectronFeedDownXibMCS2(0),
355 fHistoElectronMCGen(0),
356 fHistoBottomElectronMCGen(0),
357 fHistoCharmElectronMCGen(0),
358 fHistoLambdaMCGen(0),
359 fHistoElePtvsV0dlRS(0),
360 fHistoElePtvsV0dlRS1(0),
361 fHistoElePtvsV0dlRS2(0),
362 fHistoElePtvsV0dlRSSide(0),
363 fHistoElePtvsV0dlRSSide1(0),
364 fHistoElePtvsV0dlRSSide2(0),
365 fHistoElePtvsV0dlRSMix(0),
366 fHistoElePtvsV0dlRSMix1(0),
367 fHistoElePtvsV0dlRSMix2(0),
368 fHistoElePtvsV0dlWS(0),
369 fHistoElePtvsV0dlWS1(0),
370 fHistoElePtvsV0dlWS2(0),
371 fHistoElePtvsV0dlWSSide(0),
372 fHistoElePtvsV0dlWSSide1(0),
373 fHistoElePtvsV0dlWSSide2(0),
374 fHistoElePtvsV0dlWSMix(0),
375 fHistoElePtvsV0dlWSMix1(0),
376 fHistoElePtvsV0dlWSMix2(0),
377 fHistoElePtvsV0dlMCS(0),
378 fHistoElePtvsV0dlMCS1(0),
379 fHistoElePtvsV0dlMCS2(0),
380 fHistoElePtvsV0dlFeedDownXic0MCS(0),
381 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
382 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
383 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
384 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
385 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
386 fHistoElePtvsV0dcaRS(0),
387 fHistoElePtvsV0dcaRS1(0),
388 fHistoElePtvsV0dcaRS2(0),
389 fHistoElePtvsV0dcaRSSide(0),
390 fHistoElePtvsV0dcaRSSide1(0),
391 fHistoElePtvsV0dcaRSSide2(0),
392 fHistoElePtvsV0dcaRSMix(0),
393 fHistoElePtvsV0dcaRSMix1(0),
394 fHistoElePtvsV0dcaRSMix2(0),
395 fHistoElePtvsV0dcaWS(0),
396 fHistoElePtvsV0dcaWS1(0),
397 fHistoElePtvsV0dcaWS2(0),
398 fHistoElePtvsV0dcaWSSide(0),
399 fHistoElePtvsV0dcaWSSide1(0),
400 fHistoElePtvsV0dcaWSSide2(0),
401 fHistoElePtvsV0dcaWSMix(0),
402 fHistoElePtvsV0dcaWSMix1(0),
403 fHistoElePtvsV0dcaWSMix2(0),
404 fHistoElePtvsV0dcaMCS(0),
405 fHistoElePtvsV0dcaMCS1(0),
406 fHistoElePtvsV0dcaMCS2(0),
407 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
408 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
409 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
410 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
411 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
412 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
413 fHistoEleLambdaPtvsV0dlRS(0),
414 fHistoEleLambdaPtvsV0dlRS1(0),
415 fHistoEleLambdaPtvsV0dlRS2(0),
416 fHistoEleLambdaPtvsV0dlRSSide(0),
417 fHistoEleLambdaPtvsV0dlRSSide1(0),
418 fHistoEleLambdaPtvsV0dlRSSide2(0),
419 fHistoEleLambdaPtvsV0dlRSMix(0),
420 fHistoEleLambdaPtvsV0dlRSMix1(0),
421 fHistoEleLambdaPtvsV0dlRSMix2(0),
422 fHistoEleLambdaPtvsV0dlWS(0),
423 fHistoEleLambdaPtvsV0dlWS1(0),
424 fHistoEleLambdaPtvsV0dlWS2(0),
425 fHistoEleLambdaPtvsV0dlWSSide(0),
426 fHistoEleLambdaPtvsV0dlWSSide1(0),
427 fHistoEleLambdaPtvsV0dlWSSide2(0),
428 fHistoEleLambdaPtvsV0dlWSMix(0),
429 fHistoEleLambdaPtvsV0dlWSMix1(0),
430 fHistoEleLambdaPtvsV0dlWSMix2(0),
431 fHistoEleLambdaPtvsV0dlMCS(0),
432 fHistoEleLambdaPtvsV0dlMCS1(0),
433 fHistoEleLambdaPtvsV0dlMCS2(0),
434 fHistoEleLambdaPtvsV0dlPromptMCS(0),
435 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
436 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
437 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
438 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
439 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
440 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
441 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
442 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
443 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
444 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
445 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
446 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
447 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
448 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
449 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
450 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
451 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
452 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
453 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
454 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
455 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
456 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
457 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
458 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
459 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
460 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
461 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
462 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
463 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
464 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
465 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
466 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
467 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
468 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
469 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
470 fHistoEleLambdaPtvsV0dlXibMCS(0),
471 fHistoEleLambdaPtvsV0dlXibMCS1(0),
472 fHistoEleLambdaPtvsV0dlXibMCS2(0),
473 fHistoEleLambdaPtvsV0dlRSAway(0),
474 fHistoEleLambdaPtvsV0dlRS1Away(0),
475 fHistoEleLambdaPtvsV0dlRS2Away(0),
476 fHistoEleLambdaPtvsV0dlRSSideAway(0),
477 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
478 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
479 fHistoEleLambdaPtvsV0dlRSMixAway(0),
480 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
481 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
482 fHistoEleLambdaPtvsV0dlWSAway(0),
483 fHistoEleLambdaPtvsV0dlWS1Away(0),
484 fHistoEleLambdaPtvsV0dlWS2Away(0),
485 fHistoEleLambdaPtvsV0dlWSSideAway(0),
486 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
487 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
488 fHistoEleLambdaPtvsV0dlWSMixAway(0),
489 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
490 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
491 fHistoResponseElePt(0),
492 fHistoResponseElePt1(0),
493 fHistoResponseElePt2(0),
494 fHistoResponseEleLambdaPt(0),
495 fHistoResponseEleLambdaPt1(0),
496 fHistoResponseEleLambdaPt2(0),
497 fHistoResponseEleLambdaPtFeeddownXic0(0),
498 fHistoResponseEleLambdaPtFeeddownXic01(0),
499 fHistoResponseEleLambdaPtFeeddownXic02(0),
500 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
501 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
502 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
503 fHistoResponseEleLambdaPtFeeddownSigma(0),
504 fHistoResponseEleLambdaPtFeeddownSigma1(0),
505 fHistoResponseEleLambdaPtFeeddownSigma2(0),
506 fHistoResponseEleLambdaPtLambdab(0),
507 fHistoResponseElePtLambdab(0),
508 fHistoResponseLambdaPtLambdab(0),
509 fHistoResponseLambdaPt(0),
510 fHistoResponseLambdaPtFeeddownXic0(0),
511 fHistoResponseLambdaPtFeeddownXicPlus(0),
512 fHistoResponseLambdaPtFeeddownSigma(0),
513 fHistoResponseLambdaPtvsEleLambdaPt(0),
514 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
515 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
516 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
517 fHistoEleLambdaPtvsRapidityRS(0),
518 fHistoEleLambdaPtvsRapidityWS(0),
519 fHistoEleLambdaPtvsRapidityMCS(0),
521 fHistoElectronPi0Total(0),
522 fHistoElectronPi0Tag(0),
524 fHistoElectronEtaTotal(0),
525 fHistoElectronEtaTag(0),
527 fHistonEvtvsRunNumber(0),
528 fHistonElevsRunNumber(0),
529 fHistonLambdavsRunNumber(0),
530 fHistoMCEventType(0),
531 fHistoMCDeltaPhiccbar(0),
532 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
533 fHistodPhiSdEtaSElectronProtonR125RS(0),
534 fHistodPhiSdEtaSElectronProtonR125WS(0),
535 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
536 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
537 fHistodPhiSdEtaSElectronPionR125RS(0),
538 fHistodPhiSdEtaSElectronPionR125WS(0),
539 fHistodPhiSdEtaSElectronPionR125RSMix(0),
540 fHistodPhiSdEtaSElectronPionR125WSMix(0),
542 fMixWithoutConversionFlag(kFALSE),
543 fNumberOfEventsForMixing (5),
549 fElectronTracks(0x0),
556 fElectronCutVarsArray(0x0),
557 fV0CutVarsArray1(0x0),
558 fV0CutVarsArray2(0x0)
563 for(Int_t i=0;i<17;i++){
568 for(Int_t i=0;i<8;i++){
576 Bool_t writeVariableTree) :
577 AliAnalysisTaskSE(name),
586 fIsEventSelected(kFALSE),
587 fWriteVariableTree(writeVariableTree),
588 fWriteEachVariableTree(kFALSE),
589 fWriteMCVariableTree(kFALSE),
591 fEleVariablesTree(0),
594 fMCEleVariablesTree(0),
595 fMCV0VariablesTree(0),
596 fMCGenPairVariablesTree(0),
597 fCorrelationVariablesTree(0),
598 fReconstructPrimVert(kFALSE),
604 fCandidateVariables(),
605 fCandidateEleVariables(),
606 fCandidateV0Variables(),
607 fCandidateMCVariables(),
608 fCandidateMCEleVariables(),
609 fCandidateMCV0Variables(),
610 fCandidateMCGenPairVariables(),
611 fCorrelationVariables(),
619 fUseCentralityV0M(kFALSE),
622 fMCDoPairAnalysis(kFALSE),
623 fHistoEleLambdaMass(0),
624 fHistoEleLambdaMassRS(0),
625 fHistoEleLambdaMassRS1(0),
626 fHistoEleLambdaMassRS2(0),
627 fHistoEleLambdaMassWS(0),
628 fHistoEleLambdaMassWS1(0),
629 fHistoEleLambdaMassWS2(0),
630 fHistoEleLambdaMassRSMix(0),
631 fHistoEleLambdaMassRSMix1(0),
632 fHistoEleLambdaMassRSMix2(0),
633 fHistoEleLambdaMassWSMix(0),
634 fHistoEleLambdaMassWSMix1(0),
635 fHistoEleLambdaMassWSMix2(0),
636 fHistoEleLambdaMassRSSide(0),
637 fHistoEleLambdaMassRSSide1(0),
638 fHistoEleLambdaMassRSSide2(0),
639 fHistoEleLambdaMassWSSide(0),
640 fHistoEleLambdaMassWSSide1(0),
641 fHistoEleLambdaMassWSSide2(0),
642 fHistoEleLambdaMassRSAway(0),
643 fHistoEleLambdaMassRS1Away(0),
644 fHistoEleLambdaMassRS2Away(0),
645 fHistoEleLambdaMassWSAway(0),
646 fHistoEleLambdaMassWS1Away(0),
647 fHistoEleLambdaMassWS2Away(0),
648 fHistoEleLambdaMassRSMixAway(0),
649 fHistoEleLambdaMassRSMix1Away(0),
650 fHistoEleLambdaMassRSMix2Away(0),
651 fHistoEleLambdaMassWSMixAway(0),
652 fHistoEleLambdaMassWSMix1Away(0),
653 fHistoEleLambdaMassWSMix2Away(0),
654 fHistoEleLambdaMassRSSideAway(0),
655 fHistoEleLambdaMassRSSide1Away(0),
656 fHistoEleLambdaMassRSSide2Away(0),
657 fHistoEleLambdaMassWSSideAway(0),
658 fHistoEleLambdaMassWSSide1Away(0),
659 fHistoEleLambdaMassWSSide2Away(0),
660 fHistoEleLambdaMassvsElePtRS(0),
661 fHistoEleLambdaMassvsElePtWS(0),
662 fHistoEleLambdaMassvsElePtRSMix(0),
663 fHistoEleLambdaMassvsElePtWSMix(0),
664 fHistoEleLambdaMassvsElePtRSSide(0),
665 fHistoEleLambdaMassvsElePtWSSide(0),
666 fHistoEleLambdaMassvsElePtRS1(0),
667 fHistoEleLambdaMassvsElePtWS1(0),
668 fHistoEleLambdaMassvsElePtRSMix1(0),
669 fHistoEleLambdaMassvsElePtWSMix1(0),
670 fHistoEleLambdaMassvsElePtRSSide1(0),
671 fHistoEleLambdaMassvsElePtWSSide1(0),
672 fHistoEleLambdaMassvsElePtRS2(0),
673 fHistoEleLambdaMassvsElePtWS2(0),
674 fHistoEleLambdaMassvsElePtRSMix2(0),
675 fHistoEleLambdaMassvsElePtWSMix2(0),
676 fHistoEleLambdaMassvsElePtRSSide2(0),
677 fHistoEleLambdaMassvsElePtWSSide2(0),
682 fHistoEleLambdaMassMCS(0),
683 fHistoEleLambdaMassLambdabMCS(0),
684 fHistoEleLambdaMassXibMCS(0),
685 fHistoEleLambdaMassMCGen(0),
686 fHistoEleLambdaMassvsElePtMCS(0),
687 fHistoEleLambdaMassvsElePtMCGen(0),
688 fHistoEleLambdaMassvsElePtMCS1(0),
689 fHistoEleLambdaMassvsElePtMCGen1(0),
690 fHistoEleLambdaMassvsElePtMCS2(0),
691 fHistoEleLambdaMassvsElePtMCGen2(0),
694 fHistoElePtvsEtaRS(0),
695 fHistoElePtvsEtaWS(0),
696 fHistoElePtvsEtaRSMix(0),
697 fHistoElePtvsEtaWSMix(0),
698 fHistoElePtvsEtaMCS(0),
699 fHistoElePtvsEtaMCGen(0),
700 fHistoElePtvsLambdaPtRS(0),
701 fHistoElePtvsLambdaPtWS(0),
702 fHistoElePtvsLambdaPtRSMix(0),
703 fHistoElePtvsLambdaPtWSMix(0),
704 fHistoElePtvsLambdaPtMCS(0),
705 fHistoElePtvsLambdaPtvsLcPtMCS(0),
706 fHistoElePtvsLambdaPtMCGen(0),
707 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
708 fHistoElePtvsLambdaPtMCLcGen(0),
709 fHistoElePtvsd0RS(0),
710 fHistoElePtvsd0WS(0),
711 fHistoElePtvsd0RSMix(0),
712 fHistoElePtvsd0WSMix(0),
713 fHistoElePtvsd0MCS(0),
714 fHistoElePtvsd0PromptMCS(0),
715 fHistoElePtvsd0BFeeddownMCS(0),
716 fHistoEleLambdaMassFeeddownXic0MCS(0),
717 fHistoEleLambdaMassFeeddownXic0MCGen(0),
718 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
719 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
720 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
721 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
722 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
723 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
724 fHistoElePtFeeddownXic0MCS(0),
725 fHistoElePtFeeddownXic0MCGen(0),
726 fHistoElePtvsEtaFeeddownXic0MCS(0),
727 fHistoElePtvsEtaFeeddownXic0MCGen(0),
728 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
729 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
730 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
731 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
732 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
733 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
734 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
735 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
736 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
737 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
738 fHistoElePtFeeddownXicPlusMCS(0),
739 fHistoElePtFeeddownXicPlusMCGen(0),
740 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
741 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
742 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
743 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
746 fHistoBachPtMCGen(0),
748 fHistoLambdaMassvsPt(0),
749 fHistoLambdaMassvsPtMCS(0),
750 fHistoLambdaMassvsPtMCGen(0),
751 fHistoLambdaPtvsEta(0),
752 fHistoLambdaPtvsEtaMCS(0),
753 fHistoLambdaPtvsEtaMCGen(0),
754 fHistoLambdaArmenteros(0),
755 fHistoLambdaArmenterosMCS(0),
756 fHistoLambdaPtvsDl(0),
757 fHistoLambdaPtvsDlSide(0),
758 fHistoLambdaPtvsDlMCS(0),
759 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
760 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
761 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
762 fHistoK0sMassvsPt(0),
763 fHistoElectronTPCPID(0),
764 fHistoElectronTOFPID(0),
765 fHistoElectronTPCSelPID(0),
766 fHistoElectronTOFSelPID(0),
767 fHistoElectronTPCPIDSelTOF(0),
768 fHistoElectronTPCPIDSelTOFSmallEta(0),
769 fHistoElectronTPCPIDSelTOFLargeEta(0),
770 fHistoMassConversionsMin(0),
771 fHistoMassConversionsSameSignMin(0),
772 fHistoElectronQovPtvsPhi(0),
773 fHistoLambdaQovPtvsPhi(0),
780 fHistoLcMCGenWithEleA(0),
781 fHistoLcMCGenWithEleB(0),
782 fHistoLbMCGenWithEleA(0),
783 fHistoLbMCGenWithEleB(0),
784 fHistoLbMCGenWithLambdac(0),
791 fHistoFeedDownXic0MCGen(0),
792 fHistoFeedDownXic0MCGen1(0),
793 fHistoFeedDownXic0MCGen2(0),
794 fHistoFeedDownXic0MCS(0),
795 fHistoFeedDownXic0MCS1(0),
796 fHistoFeedDownXic0MCS2(0),
797 fHistoFeedDownXicPlusMCGen(0),
798 fHistoFeedDownXicPlusMCGen1(0),
799 fHistoFeedDownXicPlusMCGen2(0),
800 fHistoFeedDownXicPlusMCS(0),
801 fHistoFeedDownXicPlusMCS1(0),
802 fHistoFeedDownXicPlusMCS2(0),
803 fHistoFeedDownSigmaMCGen(0),
804 fHistoFeedDownSigmaMCGen1(0),
805 fHistoFeedDownSigmaMCGen2(0),
806 fHistoFeedDownSigmaMCS(0),
807 fHistoFeedDownSigmaMCS1(0),
808 fHistoFeedDownSigmaMCS2(0),
809 fHistoLcElectronMCGen(0),
810 fHistoLcElectronMCGen1(0),
811 fHistoLcElectronMCGen2(0),
812 fHistoLcElectronMCS(0),
813 fHistoLcElectronMCS1(0),
814 fHistoLcElectronMCS2(0),
815 fHistoLbElectronMCGen(0),
816 fHistoLbElectronMCGen1(0),
817 fHistoLbElectronMCGen2(0),
818 fHistoLbElectronMCS(0),
819 fHistoLbElectronMCS1(0),
820 fHistoLbElectronMCS2(0),
821 fHistoElectronFeedDownXic0MCGen(0),
822 fHistoElectronFeedDownXic0MCGen1(0),
823 fHistoElectronFeedDownXic0MCGen2(0),
824 fHistoElectronFeedDownXic0MCS(0),
825 fHistoElectronFeedDownXic0MCS1(0),
826 fHistoElectronFeedDownXic0MCS2(0),
827 fHistoElectronFeedDownXicPlusMCGen(0),
828 fHistoElectronFeedDownXicPlusMCGen1(0),
829 fHistoElectronFeedDownXicPlusMCGen2(0),
830 fHistoElectronFeedDownXicPlusMCS(0),
831 fHistoElectronFeedDownXicPlusMCS1(0),
832 fHistoElectronFeedDownXicPlusMCS2(0),
833 fHistoElectronFeedDownXibMCGen(0),
834 fHistoElectronFeedDownXibMCGen1(0),
835 fHistoElectronFeedDownXibMCGen2(0),
836 fHistoElectronFeedDownXibMCS(0),
837 fHistoElectronFeedDownXibMCS1(0),
838 fHistoElectronFeedDownXibMCS2(0),
839 fHistoElectronMCGen(0),
840 fHistoBottomElectronMCGen(0),
841 fHistoCharmElectronMCGen(0),
842 fHistoLambdaMCGen(0),
843 fHistoElePtvsV0dlRS(0),
844 fHistoElePtvsV0dlRS1(0),
845 fHistoElePtvsV0dlRS2(0),
846 fHistoElePtvsV0dlRSSide(0),
847 fHistoElePtvsV0dlRSSide1(0),
848 fHistoElePtvsV0dlRSSide2(0),
849 fHistoElePtvsV0dlRSMix(0),
850 fHistoElePtvsV0dlRSMix1(0),
851 fHistoElePtvsV0dlRSMix2(0),
852 fHistoElePtvsV0dlWS(0),
853 fHistoElePtvsV0dlWS1(0),
854 fHistoElePtvsV0dlWS2(0),
855 fHistoElePtvsV0dlWSSide(0),
856 fHistoElePtvsV0dlWSSide1(0),
857 fHistoElePtvsV0dlWSSide2(0),
858 fHistoElePtvsV0dlWSMix(0),
859 fHistoElePtvsV0dlWSMix1(0),
860 fHistoElePtvsV0dlWSMix2(0),
861 fHistoElePtvsV0dlMCS(0),
862 fHistoElePtvsV0dlMCS1(0),
863 fHistoElePtvsV0dlMCS2(0),
864 fHistoElePtvsV0dlFeedDownXic0MCS(0),
865 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
866 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
867 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
868 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
869 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
870 fHistoElePtvsV0dcaRS(0),
871 fHistoElePtvsV0dcaRS1(0),
872 fHistoElePtvsV0dcaRS2(0),
873 fHistoElePtvsV0dcaRSSide(0),
874 fHistoElePtvsV0dcaRSSide1(0),
875 fHistoElePtvsV0dcaRSSide2(0),
876 fHistoElePtvsV0dcaRSMix(0),
877 fHistoElePtvsV0dcaRSMix1(0),
878 fHistoElePtvsV0dcaRSMix2(0),
879 fHistoElePtvsV0dcaWS(0),
880 fHistoElePtvsV0dcaWS1(0),
881 fHistoElePtvsV0dcaWS2(0),
882 fHistoElePtvsV0dcaWSSide(0),
883 fHistoElePtvsV0dcaWSSide1(0),
884 fHistoElePtvsV0dcaWSSide2(0),
885 fHistoElePtvsV0dcaWSMix(0),
886 fHistoElePtvsV0dcaWSMix1(0),
887 fHistoElePtvsV0dcaWSMix2(0),
888 fHistoElePtvsV0dcaMCS(0),
889 fHistoElePtvsV0dcaMCS1(0),
890 fHistoElePtvsV0dcaMCS2(0),
891 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
892 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
893 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
894 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
895 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
896 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
897 fHistoEleLambdaPtvsV0dlRS(0),
898 fHistoEleLambdaPtvsV0dlRS1(0),
899 fHistoEleLambdaPtvsV0dlRS2(0),
900 fHistoEleLambdaPtvsV0dlRSSide(0),
901 fHistoEleLambdaPtvsV0dlRSSide1(0),
902 fHistoEleLambdaPtvsV0dlRSSide2(0),
903 fHistoEleLambdaPtvsV0dlRSMix(0),
904 fHistoEleLambdaPtvsV0dlRSMix1(0),
905 fHistoEleLambdaPtvsV0dlRSMix2(0),
906 fHistoEleLambdaPtvsV0dlWS(0),
907 fHistoEleLambdaPtvsV0dlWS1(0),
908 fHistoEleLambdaPtvsV0dlWS2(0),
909 fHistoEleLambdaPtvsV0dlWSSide(0),
910 fHistoEleLambdaPtvsV0dlWSSide1(0),
911 fHistoEleLambdaPtvsV0dlWSSide2(0),
912 fHistoEleLambdaPtvsV0dlWSMix(0),
913 fHistoEleLambdaPtvsV0dlWSMix1(0),
914 fHistoEleLambdaPtvsV0dlWSMix2(0),
915 fHistoEleLambdaPtvsV0dlMCS(0),
916 fHistoEleLambdaPtvsV0dlMCS1(0),
917 fHistoEleLambdaPtvsV0dlMCS2(0),
918 fHistoEleLambdaPtvsV0dlPromptMCS(0),
919 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
920 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
921 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
922 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
923 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
924 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
925 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
926 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
927 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
928 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
929 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
930 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
931 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
932 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
933 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
934 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
935 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
936 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
937 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
938 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
939 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
940 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
941 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
942 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
943 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
944 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
945 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
946 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
947 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
948 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
949 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
950 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
951 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
952 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
953 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
954 fHistoEleLambdaPtvsV0dlXibMCS(0),
955 fHistoEleLambdaPtvsV0dlXibMCS1(0),
956 fHistoEleLambdaPtvsV0dlXibMCS2(0),
957 fHistoEleLambdaPtvsV0dlRSAway(0),
958 fHistoEleLambdaPtvsV0dlRS1Away(0),
959 fHistoEleLambdaPtvsV0dlRS2Away(0),
960 fHistoEleLambdaPtvsV0dlRSSideAway(0),
961 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
962 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
963 fHistoEleLambdaPtvsV0dlRSMixAway(0),
964 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
965 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
966 fHistoEleLambdaPtvsV0dlWSAway(0),
967 fHistoEleLambdaPtvsV0dlWS1Away(0),
968 fHistoEleLambdaPtvsV0dlWS2Away(0),
969 fHistoEleLambdaPtvsV0dlWSSideAway(0),
970 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
971 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
972 fHistoEleLambdaPtvsV0dlWSMixAway(0),
973 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
974 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
975 fHistoResponseElePt(0),
976 fHistoResponseElePt1(0),
977 fHistoResponseElePt2(0),
978 fHistoResponseEleLambdaPt(0),
979 fHistoResponseEleLambdaPt1(0),
980 fHistoResponseEleLambdaPt2(0),
981 fHistoResponseEleLambdaPtFeeddownXic0(0),
982 fHistoResponseEleLambdaPtFeeddownXic01(0),
983 fHistoResponseEleLambdaPtFeeddownXic02(0),
984 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
985 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
986 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
987 fHistoResponseEleLambdaPtFeeddownSigma(0),
988 fHistoResponseEleLambdaPtFeeddownSigma1(0),
989 fHistoResponseEleLambdaPtFeeddownSigma2(0),
990 fHistoResponseEleLambdaPtLambdab(0),
991 fHistoResponseElePtLambdab(0),
992 fHistoResponseLambdaPtLambdab(0),
993 fHistoResponseLambdaPt(0),
994 fHistoResponseLambdaPtFeeddownXic0(0),
995 fHistoResponseLambdaPtFeeddownXicPlus(0),
996 fHistoResponseLambdaPtFeeddownSigma(0),
997 fHistoResponseLambdaPtvsEleLambdaPt(0),
998 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
999 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
1000 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
1001 fHistoEleLambdaPtvsRapidityRS(0),
1002 fHistoEleLambdaPtvsRapidityWS(0),
1003 fHistoEleLambdaPtvsRapidityMCS(0),
1005 fHistoElectronPi0Total(0),
1006 fHistoElectronPi0Tag(0),
1008 fHistoElectronEtaTotal(0),
1009 fHistoElectronEtaTag(0),
1011 fHistonEvtvsRunNumber(0),
1012 fHistonElevsRunNumber(0),
1013 fHistonLambdavsRunNumber(0),
1014 fHistoMCEventType(0),
1015 fHistoMCDeltaPhiccbar(0),
1016 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
1017 fHistodPhiSdEtaSElectronProtonR125RS(0),
1018 fHistodPhiSdEtaSElectronProtonR125WS(0),
1019 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
1020 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
1021 fHistodPhiSdEtaSElectronPionR125RS(0),
1022 fHistodPhiSdEtaSElectronPionR125WS(0),
1023 fHistodPhiSdEtaSElectronPionR125RSMix(0),
1024 fHistodPhiSdEtaSElectronPionR125WSMix(0),
1026 fMixWithoutConversionFlag(kFALSE),
1027 fNumberOfEventsForMixing (5),
1033 fElectronTracks(0x0),
1040 fElectronCutVarsArray(0x0),
1041 fV0CutVarsArray1(0x0),
1042 fV0CutVarsArray2(0x0)
1047 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
1049 for(Int_t i=0;i<17;i++){
1054 for(Int_t i=0;i<8;i++){
1058 DefineOutput(1,TList::Class());
1059 DefineOutput(2,TList::Class());
1060 DefineOutput(3,TList::Class());
1061 DefineOutput(4,TTree::Class());
1062 DefineOutput(5,TTree::Class());
1063 DefineOutput(6,TTree::Class());
1064 DefineOutput(7,TTree::Class());
1065 DefineOutput(8,AliNormalizationCounter::Class());
1066 DefineOutput(9,TTree::Class());
1067 DefineOutput(10,TTree::Class());
1068 DefineOutput(11,TTree::Class());
1069 DefineOutput(12,TTree::Class());
1077 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1172 if (fDebug > 1) AliInfo(
"Init");
1191 AliError(
"NO EVENT FOUND!");
1194 AliAODEvent* aodEvent =
dynamic_cast<AliAODEvent*
>(fInputEvent);
1202 fBzkG = (Double_t)aodEvent->GetMagneticField();
1203 AliKFParticle::SetField(
fBzkG);
1204 if (TMath::Abs(
fBzkG)<0.001) {
1215 TClonesArray *mcArray = 0;
1216 AliAODMCHeader *mcHeader=0;
1219 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1221 AliError(
"Could not find Monte-Carlo in AOD");
1227 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1229 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1234 Double_t zMCVertex = mcHeader->GetVtxZ();
1236 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1250 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1253 Double_t pos[3],cov[6];
1255 fVtx1->GetCovarianceMatrix(cov);
1256 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1260 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1267 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1268 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1269 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1270 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1271 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1285 AliCentrality *cent = aodEvent->GetCentrality();
1286 fCentrality = cent->GetCentralityPercentile(
"V0M");
1290 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1297 Int_t runnumber_offset = 0;
1298 Int_t runnumber = aodEvent->GetRunNumber();
1299 if(runnumber<=131000&&runnumber>=114000){
1300 runnumber_offset = 114000;
1301 }
else if(runnumber<=196000&&runnumber>=195000){
1302 runnumber_offset = 195000;
1303 }
else if(runnumber<=170593&&runnumber>=167902){
1304 runnumber_offset = 167902;
1349 AliAnalysisTaskSE::Terminate();
1351 fOutput =
dynamic_cast<TList*
> (GetOutputData(1));
1353 AliError(
"fOutput not available");
1357 fOutputAll =
dynamic_cast<TList*
> (GetOutputData(3));
1359 AliError(
"fOutputAll not available");
1414 TString normName=
"NormalizationCounter";
1415 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1416 if(cont)normName=(TString)cont->GetName();
1438 fEventBuffer[i]=
new TTree(Form(
"EventBuffer_%d",i),
"Temporary buffer for event mixing");
1462 AliAODEvent *aodEvent, TClonesArray *mcArray
1469 if(fElectronTracks) fElectronTracks->Delete();
1470 if(fV0Tracks1) fV0Tracks1->Delete();
1471 if(fV0Tracks2) fV0Tracks2->Delete();
1472 if(fElectronCutVarsArray) fElectronCutVarsArray->Delete();
1473 if(fV0CutVarsArray1) fV0CutVarsArray1->Delete();
1474 if(fV0CutVarsArray2) fV0CutVarsArray2->Delete();
1475 fV0dlArray1.clear();
1476 fV0dlArray2.clear();
1477 fV0dcaArray1.clear();
1478 fV0dcaArray2.clear();
1481 ResetGlobalTrackReference();
1483 for (Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1485 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1486 if (!track)
continue;
1489 StoreGlobalTrackReference(track,iTrack);
1496 Int_t nV0s= aodEvent->GetNumberOfV0s();
1497 Int_t nTracks= aodEvent->GetNumberOfTracks();
1499 Bool_t seleTrkFlags[nTracks];
1501 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1503 Bool_t seleV0Flags[nV0s];
1505 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1507 Int_t runnumber_offset = 0;
1508 Int_t runnumber = aodEvent->GetRunNumber();
1509 if(runnumber<=131000&&runnumber>=114000){
1510 runnumber_offset = 114000;
1511 }
else if(runnumber<=196000&&runnumber>=195000){
1512 runnumber_offset = 195000;
1513 }
else if(runnumber<=170593&&runnumber>=167902){
1514 runnumber_offset = 167902;
1516 fHistonElevsRunNumber->Fill(runnumber-runnumber_offset,nSeleTrks);
1517 fHistonLambdavsRunNumber->Fill(runnumber-runnumber_offset,nSeleV0);
1522 for (Int_t iv0 = 0; iv0<nV0s; iv0++) {
1523 if(!seleV0Flags[iv0])
continue;
1524 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1527 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1528 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1533 for (Int_t itrk = 0; itrk<nTracks; itrk++) {
1534 if(!seleTrkFlags[itrk])
continue;
1535 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1543 AliAODTrack *trkpid = 0;
1544 if(fAnalCuts->GetProdAODFilterBit()==7){
1545 trkpid = fGTI[-trk->GetID()-1];
1550 Int_t cpid = cptrack->GetID();
1551 Int_t cnid = cntrack->GetID();
1552 Int_t lpid = trkpid->GetID();
1553 if((cpid==lpid)||(cnid==lpid))
continue;
1555 AliAODVertex *secVert = ReconstructSecondaryVertex(v0,trk,aodEvent);
1556 if(!secVert)
continue;
1563 FillROOTObjects(elobj,v0,trk,trkpid,aodEvent,mcArray);
1565 elobj->GetSecondaryVtx()->RemoveDaughters();
1567 delete elobj;elobj=NULL;
1573 fEventInfo->SetString(Form(
"Ev%d_esd%d_E%d_V%d",AliAnalysisManager::GetAnalysisManager()->GetNcalls(),((AliAODHeader*)aodEvent->GetHeader())->GetEventNumberESDFile(),fElectronTracks->GetEntries(),fV0Tracks1->GetEntries()+fV0Tracks2->GetEntries()));
1574 Int_t ind=GetPoolIndex(fVtxZ,fCentrality);
1575 if(ind>=0 && ind<fNOfPools){
1576 if(fEventBuffer[ind]->GetEntries() >= fNumberOfEventsForMixing){
1577 DoEventMixingWithPools(ind);
1578 if(fEventBuffer[ind]->GetEntries() >= 20*fNumberOfEventsForMixing){
1582 fEventBuffer[ind]->Fill();
1593 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1594 fVariablesTree =
new TTree(nameoutput,
"Candidates variables tree");
1597 TString * fCandidateVariableNames =
new TString[nVar];
1599 fCandidateVariableNames[ 0]=
"Centrality";
1600 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1601 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1602 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1603 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1604 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1605 fCandidateVariableNames[ 6]=
"ElePx";
1606 fCandidateVariableNames[ 7]=
"ElePy";
1607 fCandidateVariableNames[ 8]=
"ElePz";
1608 fCandidateVariableNames[ 9]=
"V0Px";
1609 fCandidateVariableNames[10]=
"V0Py";
1610 fCandidateVariableNames[11]=
"V0Pz";
1611 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1612 fCandidateVariableNames[13]=
"MassLambda";
1613 fCandidateVariableNames[14]=
"MassAntiLambda";
1614 fCandidateVariableNames[15]=
"Eled0";
1615 fCandidateVariableNames[16]=
"V0d0";
1616 fCandidateVariableNames[17]=
"nSigmaTPCele";
1617 fCandidateVariableNames[18]=
"nSigmaTOFele";
1618 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1619 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1620 fCandidateVariableNames[21]=
"EleCharge";
1621 fCandidateVariableNames[22]=
"ProtonPx";
1622 fCandidateVariableNames[23]=
"ProtonPy";
1623 fCandidateVariableNames[24]=
"ProtonPz";
1624 fCandidateVariableNames[25]=
"PiPx";
1625 fCandidateVariableNames[26]=
"PiPy";
1626 fCandidateVariableNames[27]=
"PiPz";
1627 fCandidateVariableNames[28]=
"mcpdglc";
1628 fCandidateVariableNames[29]=
"mclablc";
1629 fCandidateVariableNames[30]=
"mcpdgmomele";
1630 fCandidateVariableNames[31]=
"mcpdgmomv0";
1631 fCandidateVariableNames[32]=
"Mixing";
1632 fCandidateVariableNames[33]=
"mcpdgele";
1633 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1634 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1635 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1636 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1637 fCandidateVariableNames[38]=
"MassK0Short";
1638 fCandidateVariableNames[39]=
"mcpdggrmomele";
1639 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1640 fCandidateVariableNames[41]=
"mcngenele";
1641 fCandidateVariableNames[42]=
"mcngenv0";
1642 fCandidateVariableNames[43]=
"mclcpx";
1643 fCandidateVariableNames[44]=
"mclcpy";
1644 fCandidateVariableNames[45]=
"mclcpz";
1645 fCandidateVariableNames[46]=
"mcelepx";
1646 fCandidateVariableNames[47]=
"mcelepy";
1647 fCandidateVariableNames[48]=
"mcelepz";
1648 fCandidateVariableNames[49]=
"mcv0px";
1649 fCandidateVariableNames[50]=
"mcv0py";
1650 fCandidateVariableNames[51]=
"mcv0pz";
1651 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1652 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1653 fCandidateVariableNames[54]=
"PrimVertx";
1654 fCandidateVariableNames[55]=
"PrimVerty";
1655 fCandidateVariableNames[56]=
"PrimVertz";
1656 fCandidateVariableNames[57]=
"V0Vertx";
1657 fCandidateVariableNames[58]=
"V0Verty";
1658 fCandidateVariableNames[59]=
"V0Vertz";
1660 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1661 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1662 fCandidateVariableNames[62]=
"DcaV0daughters";
1663 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1664 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1665 fCandidateVariableNames[65]=
"MassK0Short2";
1667 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1668 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1670 fCandidateVariableNames[68]=
"EleITSMatch";
1671 fCandidateVariableNames[69]=
"V0PosITSMatch";
1672 fCandidateVariableNames[70]=
"V0NegITSMatch";
1673 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1674 fCandidateVariableNames[72]=
"mcpdgv0";
1676 fCandidateVariableNames[73]=
"mcv0pospx";
1677 fCandidateVariableNames[74]=
"mcv0pospy";
1678 fCandidateVariableNames[75]=
"mcv0pospz";
1679 fCandidateVariableNames[76]=
"mcpdgv0pos";
1680 fCandidateVariableNames[77]=
"mcv0negpx";
1681 fCandidateVariableNames[78]=
"mcv0negpy";
1682 fCandidateVariableNames[79]=
"mcv0negpz";
1683 fCandidateVariableNames[80]=
"mcpdgv0neg";
1685 fCandidateVariableNames[81]=
"ProtonPxPV";
1686 fCandidateVariableNames[82]=
"ProtonPyPV";
1687 fCandidateVariableNames[83]=
"ProtonPzPV";
1688 fCandidateVariableNames[84]=
"PiPxPV";
1689 fCandidateVariableNames[85]=
"PiPyPV";
1690 fCandidateVariableNames[86]=
"PiPzPV";
1692 fCandidateVariableNames[87]=
"mcv0vertx";
1693 fCandidateVariableNames[88]=
"mcv0verty";
1694 fCandidateVariableNames[89]=
"mcv0vertz";
1696 fCandidateVariableNames[90]=
"OnFlyStatus";
1697 fCandidateVariableNames[91]=
"MagneticField";
1698 fCandidateVariableNames[92]=
"EvNumber";
1699 fCandidateVariableNames[93]=
"RunNumber";
1702 for (Int_t ivar=0; ivar<nVar; ivar++) {
1719 for(Int_t i=0;i<94;i++)
1724 Bool_t anti_lambda_flag = kFALSE;
1727 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1728 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1734 Double_t v0px = elobj->PxProng(1);
1735 Double_t v0py = elobj->PyProng(1);
1736 Double_t v0pz = elobj->PzProng(1);
1737 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1738 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1739 Double_t epx = elobj->PxProng(0);
1740 Double_t epy = elobj->PyProng(0);
1741 Double_t epz = elobj->PzProng(0);
1742 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1743 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1744 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
1745 Double_t Esum = Ee + Ev0;
1747 Double_t uxe = epx/mome;
1748 Double_t uye = epy/mome;
1749 Double_t uze = epz/mome;
1750 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
1751 Double_t pxv_flip = v0px + lf * uxe;
1752 Double_t pyv_flip = v0py + lf * uye;
1753 Double_t pzv_flip = v0pz + lf * uze;
1754 Double_t pxsum_flip = epx + pxv_flip;
1755 Double_t pysum_flip = epy + pyv_flip;
1756 Double_t pzsum_flip = epz + pzv_flip;
1757 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
1758 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
1760 Double_t minmass_ee = 9999.;
1762 Double_t minmasslike_ee = 9999.;
1766 UInt_t pdgdg[2]={11,3122};
1784 Double_t nSigmaTPCele=-9999.;
1785 Double_t nSigmaTOFele=-9999.;
1786 Double_t nSigmaTPCv0pr=-9999.;
1787 Double_t nSigmaTOFv0pr=-9999.;
1788 Double_t nSigmaTPCv0pi=-9999.;
1789 Double_t nSigmaTOFv0pi=-9999.;
1800 if(anti_lambda_flag){
1818 if(anti_lambda_flag){
1834 AliAODMCParticle *mclc = 0;
1835 AliAODMCParticle *mcele = 0;
1836 AliAODMCParticle *mcv0 = 0;
1838 Int_t mcpdgele_array[100];
1839 Int_t mcpdgv0_array[100];
1840 Int_t mclabelele_array[100];
1841 Int_t mclabelv0_array[100];
1842 Int_t mcngen_ele=-9999;
1843 Int_t mcngen_v0=-9999;
1846 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
1847 if(mclabelele_array[0]>=0){
1848 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
1855 if(mclabelv0_array[0]>=0){
1856 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
1864 Int_t pdgdgv0[2]={211,211};
1865 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
1867 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
1877 Int_t labPos = cptrack->GetLabel();
1879 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
1888 Int_t labNeg = cntrack->GetLabel();
1890 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
1898 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
1941 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
1942 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
1943 if(!anti_lambda_flag){
1951 Double_t posVtx[3] = {0.,0.,0.};
1952 fVtx1->GetXYZ(posVtx);
1954 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
1964 if(!anti_lambda_flag){
1989 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
1990 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
1994 cont[0] = elobj->InvMass(2,pdgdg);
1995 cont[1] = elobj->Pt();
1999 Double_t cont_flip[3];
2000 cont_flip[0] = melam_flip;
2001 cont_flip[1] = ptelam_flip;
2005 cont2[0] = elobj->InvMass(2,pdgdg);
2006 cont2[1] = trk->Pt();
2008 Double_t cont_eleptvseta[3];
2009 cont_eleptvseta[0] = trk->Pt();
2010 cont_eleptvseta[1] = trk->Eta();
2013 Double_t cont_eleptvslambdapt[3];
2014 cont_eleptvslambdapt[0] = trk->Pt();
2015 cont_eleptvslambdapt[1] = v0->Pt();
2018 Double_t cont_eleptvsd0[3];
2019 cont_eleptvsd0[0] = trk->Pt();
2020 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
2023 Double_t cont_eleptvsv0dl[3];
2024 cont_eleptvsv0dl[0] = trk->Pt();
2025 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2028 Double_t cont_elelamptvsv0dl[4];
2029 cont_elelamptvsv0dl[0] = elobj->Pt();
2030 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
2031 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
2034 Double_t cont_elelamptvsv0dl_flip[4];
2035 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
2036 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2037 cont_elelamptvsv0dl_flip[2] = 0.0;
2041 Double_t cont_eleptvsv0dca[3];
2042 cont_eleptvsv0dca[0] = trk->Pt();
2043 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
2046 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
2047 Double_t elobj_px = elobj->Px();
2048 Double_t elobj_py = elobj->Py();
2049 Double_t elobj_pz = elobj->Pz();
2050 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
2051 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
2058 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2085 for(Int_t ih=0;ih<17;ih++){
2086 Double_t cont_eleptvscutvars[3];
2087 cont_eleptvscutvars[0] = elobj->Pt();
2091 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2093 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2095 cont_eleptvscutvars[1] = nSigmaTPCele;
2097 cont_eleptvscutvars[1] = nSigmaTOFele;
2099 cont_eleptvscutvars[1] = trk->Eta();
2101 cont_eleptvscutvars[1] = trk->GetITSNcls();
2103 if(!anti_lambda_flag)
2104 cont_eleptvscutvars[1] = v0->MassLambda();
2106 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2109 lPosV0[0] = v0->DecayVertexV0X();
2110 lPosV0[1] = v0->DecayVertexV0Y();
2111 lPosV0[2] = v0->DecayVertexV0Z();
2112 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2114 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2116 if(!anti_lambda_flag)
2117 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2119 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2121 if(!anti_lambda_flag)
2122 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2124 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2126 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2128 cont_eleptvscutvars[1] = v0->MassK0Short();
2130 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2132 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2134 cont_eleptvscutvars[1] = v0->Eta();
2136 Double_t v0px = elobj->PxProng(1);
2137 Double_t v0py = elobj->PyProng(1);
2138 Double_t v0pz = elobj->PzProng(1);
2139 Double_t epx = elobj->PxProng(0);
2140 Double_t epy = elobj->PyProng(0);
2141 Double_t epz = elobj->PzProng(0);
2142 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2144 cont_eleptvscutvars[1] = -9999.;
2152 Int_t pdgcode = mclc->GetPdgCode();
2153 Double_t cont_mclc[3];
2154 cont_mclc[0] = mclc->Pt();
2155 cont_mclc[1] = mclc->Y();
2157 Double_t cont_mcele[3];
2158 cont_mcele[0] = mcele->Pt();
2159 cont_mcele[1] = mcele->Eta();
2162 Int_t labmotherlc = mclc->GetMother();
2163 Bool_t isbottomfd = kFALSE;
2165 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2166 Int_t pdgmotherlc = motherlc->GetPdgCode();
2167 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){
2171 cont2[1] = mcele->Pt();
2173 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2225 Double_t cont_eleptvslambdaptvslcpt[4];
2226 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2227 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2228 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2229 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2232 Double_t cont_allpt[4];
2233 cont_allpt[0] = mclc->Pt();
2234 cont_allpt[1] = elobj->Pt();
2235 cont_allpt[2] = trk->Pt();
2236 cont_allpt[3] = v0->Pt();
2245 for(Int_t ih=0;ih<17;ih++){
2246 Double_t cont_eleptvscutvars[3];
2247 cont_eleptvscutvars[0] = elobj->Pt();
2251 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2253 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2255 cont_eleptvscutvars[1] = nSigmaTPCele;
2257 cont_eleptvscutvars[1] = nSigmaTOFele;
2259 cont_eleptvscutvars[1] = trk->Eta();
2261 cont_eleptvscutvars[1] = trk->GetITSNcls();
2263 if(!anti_lambda_flag)
2264 cont_eleptvscutvars[1] = v0->MassLambda();
2266 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2269 lPosV0[0] = v0->DecayVertexV0X();
2270 lPosV0[1] = v0->DecayVertexV0Y();
2271 lPosV0[2] = v0->DecayVertexV0Z();
2272 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2274 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2276 if(!anti_lambda_flag)
2277 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2279 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2281 if(!anti_lambda_flag)
2282 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2284 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2286 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2288 cont_eleptvscutvars[1] = v0->MassK0Short();
2290 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2292 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2294 cont_eleptvscutvars[1] = v0->Eta();
2296 Double_t v0px = elobj->PxProng(1);
2297 Double_t v0py = elobj->PyProng(1);
2298 Double_t v0pz = elobj->PzProng(1);
2299 Double_t epx = elobj->PxProng(0);
2300 Double_t epy = elobj->PyProng(0);
2301 Double_t epz = elobj->PzProng(0);
2302 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2304 cont_eleptvscutvars[1] = -9999.;
2311 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2358 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2405 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2435 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2462 for(Int_t ih=0;ih<17;ih++){
2463 Double_t cont_eleptvscutvars[3];
2464 cont_eleptvscutvars[0] = elobj->Pt();
2468 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2470 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2472 cont_eleptvscutvars[1] = nSigmaTPCele;
2474 cont_eleptvscutvars[1] = nSigmaTOFele;
2476 cont_eleptvscutvars[1] = trk->Eta();
2478 cont_eleptvscutvars[1] = trk->GetITSNcls();
2480 if(!anti_lambda_flag)
2481 cont_eleptvscutvars[1] = v0->MassLambda();
2483 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2486 lPosV0[0] = v0->DecayVertexV0X();
2487 lPosV0[1] = v0->DecayVertexV0Y();
2488 lPosV0[2] = v0->DecayVertexV0Z();
2489 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2491 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2493 if(!anti_lambda_flag)
2494 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2496 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2498 if(!anti_lambda_flag)
2499 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2501 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2503 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2505 cont_eleptvscutvars[1] = v0->MassK0Short();
2507 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2509 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2511 cont_eleptvscutvars[1] = v0->Eta();
2513 Double_t v0px = elobj->PxProng(1);
2514 Double_t v0py = elobj->PyProng(1);
2515 Double_t v0pz = elobj->PzProng(1);
2516 Double_t epx = elobj->PxProng(0);
2517 Double_t epy = elobj->PyProng(0);
2518 Double_t epz = elobj->PzProng(0);
2519 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2521 cont_eleptvscutvars[1] = -9999.;
2529 Int_t pdgcode = mclc->GetPdgCode();
2530 Double_t cont_mclb[3];
2531 cont_mclb[0] = mclc->Pt();
2532 cont_mclb[1] = mclc->Y();
2534 Double_t cont_mcele[3];
2535 cont_mcele[0] = mcele->Pt();
2536 cont_mcele[1] = mcele->Eta();
2539 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2558 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2568 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2578 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2588 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2606 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2608 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2612 if(cont_flip[0]<2.3){
2621 if(cont_flip[0]<2.3){
2631 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2652 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2676 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2678 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2682 if(cont_flip[0]<2.3){
2691 if(cont_flip[0]<2.3){
2702 for(Int_t iv=0;iv<13;iv++){
2714 if(trk->Charge()>0){
2717 }
else if(trk->Charge()<0){
2728 Int_t pdgcode = mclc->GetPdgCode();
2729 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2732 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2735 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2738 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2741 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2744 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2747 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2750 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2753 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2795 for(Int_t i=0;i<94;i++)
2801 Double_t pxe = trke->Px();
2802 Double_t pye = trke->Py();
2803 Double_t pze = trke->Pz();
2804 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
2805 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
2807 Double_t pxv = v0->Px();
2808 Double_t pyv = v0->Py();
2809 Double_t pzv = v0->Pz();
2810 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
2811 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
2813 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
2815 Double_t pxsum = pxe + pxv;
2816 Double_t pysum = pye + pyv;
2817 Double_t pzsum = pze + pzv;
2818 Double_t Esum = Ee + Ev;
2819 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
2821 Double_t uxe = pxe/mome;
2822 Double_t uye = pye/mome;
2823 Double_t uze = pze/mome;
2824 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
2825 Double_t pxv_flip = pxv + lf * uxe;
2826 Double_t pyv_flip = pyv + lf * uye;
2827 Double_t pzv_flip = pzv + lf * uze;
2828 Double_t pxsum_flip = pxe + pxv_flip;
2829 Double_t pysum_flip = pye + pyv_flip;
2830 Double_t pzsum_flip = pze + pzv_flip;
2831 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
2832 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
2835 UInt_t pdgdg[2]={11,3122};
2888 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
2892 Double_t cont_flip[3];
2893 cont_flip[0] = mel_flip;
2894 cont_flip[1] = ptel_flip;
2899 cont2[1] = sqrt(pxe*pxe+pye*pye);
2901 Double_t cont_eleptvseta[3];
2902 cont_eleptvseta[0] = trke->Pt();
2903 cont_eleptvseta[1] = trke->Eta();
2906 Double_t cont_eleptvslambdapt[3];
2907 cont_eleptvslambdapt[0] = trke->Pt();
2908 cont_eleptvslambdapt[1] = v0->Pt();
2911 Double_t cont_eleptvsd0[3];
2912 cont_eleptvsd0[0] = trke->Pt();
2913 cont_eleptvsd0[1] = 0.;
2916 Double_t cont_eleptvsv0dl[3];
2917 cont_eleptvsv0dl[0] = trke->Pt();
2918 cont_eleptvsv0dl[1] = v0info[0];
2921 Double_t cont_eleptvsv0dca[3];
2922 cont_eleptvsv0dca[0] = trke->Pt();
2923 cont_eleptvsv0dca[1] = v0info[1];
2926 Double_t cont_elelamptvsv0dl[4];
2927 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
2928 cont_elelamptvsv0dl[1] = v0info[0];
2929 cont_elelamptvsv0dl[2] = 0.0;
2932 Double_t cont_elelamptvsv0dl_flip[4];
2933 cont_elelamptvsv0dl_flip[0] = ptel_flip;
2934 cont_elelamptvsv0dl_flip[1] = v0info[0];
2935 cont_elelamptvsv0dl_flip[2] = 0.0;
2938 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
2939 xyzR125_ele[0] = (*elevars)[0];
2940 xyzR125_ele[1] = (*elevars)[1];
2941 xyzR125_ele[2] = (*elevars)[2];
2942 xyzR125_pr[0] = (*v0vars)[0];
2943 xyzR125_pr[1] = (*v0vars)[1];
2944 xyzR125_pr[2] = (*v0vars)[2];
2945 xyzR125_pi[0] = (*v0vars)[3];
2946 xyzR125_pi[1] = (*v0vars)[4];
2947 xyzR125_pi[2] = (*v0vars)[5];
2949 Double_t rdhfcutvars[9];
2950 rdhfcutvars[0] = xyzR125_ele[0];
2951 rdhfcutvars[1] = xyzR125_ele[1];
2952 rdhfcutvars[2] = xyzR125_ele[2];
2953 rdhfcutvars[3] = xyzR125_pr[0];
2954 rdhfcutvars[4] = xyzR125_pr[1];
2955 rdhfcutvars[5] = xyzR125_pr[2];
2956 rdhfcutvars[6] = xyzR125_pi[0];
2957 rdhfcutvars[7] = xyzR125_pi[1];
2958 rdhfcutvars[8] = xyzR125_pi[2];
2969 if(((
int)trke->T())*chargepr>0){
3029 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
3031 if(((
int)trke->T())*chargepr>0){
3035 if(cont_flip[0]<2.3){
3044 if(cont_flip[0]<2.3){
3055 for(Int_t iv=0;iv<13;iv++){
3063 fCorrelationVariables[4] = v0info[0];
3064 fCorrelationVariables[5] = (*elevars)[5];
3065 fCorrelationVariables[6] = 2;
3067 if(chargepr>0) fCorrelationVariables[7] = 0;
3068 else fCorrelationVariables[7] = 2;
3069 }
else if(trke->T()<0){
3070 if(chargepr>0) fCorrelationVariables[7] = 1;
3071 else fCorrelationVariables[7] = 3;
3073 fCorrelationVariables[8] = (*elevars)[6];
3074 fCorrelationVariables[9] = (*elevars)[7];
3076 fCorrelationVariables[11] = sqrt(pxsum*pxsum+pysum*pysum);
3077 fCorrelationVariables[12] = mel;
3093 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
3097 TString * fCandidateVariableNames =
new TString[nVar];
3099 fCandidateVariableNames[ 0]=
"ElePx";
3100 fCandidateVariableNames[ 1]=
"ElePy";
3101 fCandidateVariableNames[ 2]=
"ElePz";
3102 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
3103 fCandidateVariableNames[ 4]=
"ITSNcls";
3104 fCandidateVariableNames[ 5]=
"TPCNcls";
3105 fCandidateVariableNames[ 6]=
"TPCNclsPID";
3106 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
3107 fCandidateVariableNames[ 8]=
"d0R";
3108 fCandidateVariableNames[ 9]=
"d0Z";
3109 fCandidateVariableNames[10]=
"ITSClusterMap";
3110 fCandidateVariableNames[11]=
"nSigmaTPCele";
3111 fCandidateVariableNames[12]=
"nSigmaTOFele";
3112 fCandidateVariableNames[13]=
"nSigmaTPCpi";
3113 fCandidateVariableNames[14]=
"nSigmaTPCka";
3114 fCandidateVariableNames[15]=
"nSigmaTPCpr";
3115 fCandidateVariableNames[16]=
"EvNumber";
3116 fCandidateVariableNames[17]=
"EleCharge";
3117 fCandidateVariableNames[18]=
"ElePdgCode";
3118 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3119 fCandidateVariableNames[20]=
"mcelepx";
3120 fCandidateVariableNames[21]=
"mcelepy";
3121 fCandidateVariableNames[22]=
"mcelepz";
3122 fCandidateVariableNames[23]=
"Centrality";
3123 fCandidateVariableNames[24]=
"PrimVertZ";
3124 fCandidateVariableNames[25]=
"RunNumber";
3126 for (Int_t ivar=0; ivar<nVar; ivar++) {
3144 Double_t d0z0[2],covd0z0[3];
3145 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3149 pv[0] =
fVtx1->GetX();
3150 pv[1] =
fVtx1->GetY();
3151 pv[2] =
fVtx1->GetZ();
3152 Double_t xyzR125[3]={9999.,9999.,9999.};
3154 TVector *varvec =
new TVector(8);
3155 (*varvec)[0] = xyzR125[0];
3156 (*varvec)[1] = xyzR125[1];
3157 (*varvec)[2] = xyzR125[2];
3160 (*varvec)[5] = d0z0[0];
3161 (*varvec)[6] = convtype;
3162 (*varvec)[7] = mcetype;
3164 fElectronTracks->AddLast(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3168 Int_t pdgEle = -9999;
3169 Int_t pdgEleMother = -9999;
3170 Float_t mcelepx = -9999;
3171 Float_t mcelepy = -9999;
3172 Float_t mcelepz = -9999;
3175 Int_t labEle = trk->GetLabel();
3176 if(labEle<0)
return;
3177 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3179 pdgEle = mcetrk->GetPdgCode();
3180 mcelepx = mcetrk->Px();
3181 mcelepy = mcetrk->Py();
3182 mcelepz = mcetrk->Pz();
3199 for(Int_t i=0;i<26;i++){
3210 if(trk->GetTPCNclsF()>0)
3215 Int_t itsmap = trk->GetITSClusterMap();
3218 Bool_t spdfirst = (itsmap & bit1) == bit1;
3219 Bool_t spdsecond = (itsmap & bit2) == bit2;
3257 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3261 TString * fCandidateVariableNames =
new TString[nVar];
3263 fCandidateVariableNames[ 0]=
"V0Px";
3264 fCandidateVariableNames[ 1]=
"V0Py";
3265 fCandidateVariableNames[ 2]=
"V0Pz";
3266 fCandidateVariableNames[ 3]=
"MassLambda";
3267 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3268 fCandidateVariableNames[ 5]=
"ProtonPx";
3269 fCandidateVariableNames[ 6]=
"ProtonPy";
3270 fCandidateVariableNames[ 7]=
"ProtonPz";
3271 fCandidateVariableNames[ 8]=
"PionPx";
3272 fCandidateVariableNames[ 9]=
"PionPy";
3273 fCandidateVariableNames[10]=
"PionPz";
3274 fCandidateVariableNames[11]=
"RfidV0";
3275 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3276 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3277 fCandidateVariableNames[14]=
"DcaV0daughters";
3278 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3279 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3280 fCandidateVariableNames[17]=
"MassK0Short";
3281 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3282 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3283 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3284 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3285 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3286 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3287 fCandidateVariableNames[24]=
"V0PdgCode";
3288 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3289 fCandidateVariableNames[26]=
"mcv0px";
3290 fCandidateVariableNames[27]=
"mcv0py";
3291 fCandidateVariableNames[28]=
"mcv0pz";
3292 fCandidateVariableNames[29]=
"mcv0vertx";
3293 fCandidateVariableNames[30]=
"mcv0verty";
3294 fCandidateVariableNames[31]=
"mcv0vertz";
3295 fCandidateVariableNames[32]=
"EvNumber";
3296 fCandidateVariableNames[33]=
"Centrality";
3297 fCandidateVariableNames[34]=
"PrimVertZ";
3298 fCandidateVariableNames[35]=
"RunNumber";
3300 for (Int_t ivar=0; ivar<nVar; ivar++) {
3314 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3315 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3316 if(!cptrack)
return;
3317 if(!cntrack)
return;
3336 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3337 Double_t posVtx[3] = {0.,0.,0.};
3338 fVtx1->GetXYZ(posVtx);
3339 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3340 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3341 Double_t cont_arm[3];
3342 cont_arm[0] = v0->Pt();
3343 cont_arm[1] = v0->PtArmV0();
3344 cont_arm[2] = v0->AlphaV0();
3363 Double_t momv0x = v0->MomV0X();
3364 Double_t momv0y = v0->MomV0Y();
3365 Double_t phi_alice = atan2(momv0y,momv0x);
3366 if(phi_alice<0.) phi_alice += 2 * M_PI;
3369 Int_t v0pdgcode = -9999;
3370 Int_t v0motherpdgcode = -9999;
3371 Float_t mcv0px = -9999;
3372 Float_t mcv0py = -9999;
3373 Float_t mcv0pz = -9999;
3374 Float_t mcv0vertx = -9999;
3375 Float_t mcv0verty = -9999;
3376 Float_t mcv0vertz = -9999;
3379 Int_t pdgdgv0[2]={2212,211};
3380 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
3394 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
3395 if(!mcv0trk)
return;
3397 Bool_t hfv0_flag = kFALSE;
3398 Bool_t promptv0_flag = kFALSE;
3399 v0pdgcode = mcv0trk->GetPdgCode();
3400 Int_t labv0mother = mcv0trk->GetMother();
3402 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3404 v0motherpdgcode = motherv0->GetPdgCode();
3405 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
3408 if(abs(v0motherpdgcode)==3322){
3410 }
else if(abs(v0motherpdgcode)==3312){
3412 }
else if(abs(v0motherpdgcode)==3334){
3416 promptv0_flag = kTRUE;
3420 if(!hfv0_flag)
return;
3423 mcv0px = mcv0trk->Px();
3424 mcv0py = mcv0trk->Py();
3425 mcv0pz = mcv0trk->Pz();
3427 Int_t labptrk = cptrack->GetLabel();
3429 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
3431 mcv0vertx = mcv0ptrk->Xv();
3432 mcv0verty = mcv0ptrk->Yv();
3433 mcv0vertz = mcv0ptrk->Zv();
3440 Double_t posVtx[3] = {0.,0.,0.};
3441 fVtx1->GetXYZ(posVtx);
3442 TLorentzVector *lv =
new TLorentzVector();
3443 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3444 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
3445 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
3447 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
3449 fV0dlArray1.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3453 TVector *varvec =
new TVector(15);
3454 (*varvec)[0] = xyzR125pr[0];
3455 (*varvec)[1] = xyzR125pr[1];
3456 (*varvec)[2] = xyzR125pr[2];
3457 (*varvec)[3] = xyzR125pi[0];
3458 (*varvec)[4] = xyzR125pi[1];
3459 (*varvec)[5] = xyzR125pi[2];
3460 (*varvec)[6] = v0->DecayVertexV0X();
3461 (*varvec)[7] = v0->DecayVertexV0Y();
3462 (*varvec)[8] = v0->DecayVertexV0Z();
3463 (*varvec)[9] = v0->MomPosX();
3464 (*varvec)[10] = v0->MomPosY();
3465 (*varvec)[11] = v0->MomPosZ();
3466 (*varvec)[12] = v0->MomNegX();
3467 (*varvec)[13] = v0->MomNegY();
3468 (*varvec)[14] = v0->MomNegZ();
3471 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
3473 fV0dlArray2.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3477 TVector *varvec =
new TVector(15);
3478 (*varvec)[0] = xyzR125pr[0];
3479 (*varvec)[1] = xyzR125pr[1];
3480 (*varvec)[2] = xyzR125pr[2];
3481 (*varvec)[3] = xyzR125pi[0];
3482 (*varvec)[4] = xyzR125pi[1];
3483 (*varvec)[5] = xyzR125pi[2];
3484 (*varvec)[6] = v0->DecayVertexV0X();
3485 (*varvec)[7] = v0->DecayVertexV0Y();
3486 (*varvec)[8] = v0->DecayVertexV0Z();
3487 (*varvec)[9] = v0->MomNegX();
3488 (*varvec)[10] = v0->MomNegY();
3489 (*varvec)[11] = v0->MomNegZ();
3490 (*varvec)[12] = v0->MomPosX();
3491 (*varvec)[13] = v0->MomPosY();
3492 (*varvec)[14] = v0->MomPosZ();
3499 for(Int_t i=0;i<36;i++){
3509 Bool_t isparticle = kTRUE;
3529 lPosV0[0] = v0->DecayVertexV0X();
3530 lPosV0[1] = v0->DecayVertexV0Y();
3531 lPosV0[2] = v0->DecayVertexV0Z();
3532 Double_t decayvertV0 = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
3535 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
3536 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
3566 if(cptrack->GetTPCNclsF()>0)
3567 fCandidateV0Variables[22] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3568 if(cntrack->GetTPCNclsF()>0)
3569 fCandidateV0Variables[23] =(Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3573 if(cntrack->GetTPCNclsF()>0)
3574 fCandidateV0Variables[22] = (Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3575 if(cptrack->GetTPCNclsF()>0)
3576 fCandidateV0Variables[23] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3601 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
3605 TString * fCandidateVariableNames =
new TString[nVar];
3607 fCandidateVariableNames[ 0]=
"Centrality";
3608 fCandidateVariableNames[ 1]=
"DecayType";
3609 fCandidateVariableNames[ 2]=
"LcPx";
3610 fCandidateVariableNames[ 3]=
"LcPy";
3611 fCandidateVariableNames[ 4]=
"LcPz";
3612 fCandidateVariableNames[ 5]=
"ElePx";
3613 fCandidateVariableNames[ 6]=
"ElePy";
3614 fCandidateVariableNames[ 7]=
"ElePz";
3615 fCandidateVariableNames[ 8]=
"V0Px";
3616 fCandidateVariableNames[ 9]=
"V0Py";
3617 fCandidateVariableNames[10]=
"V0Pz";
3619 for (Int_t ivar=0; ivar<nVar; ivar++) {
3631 if(!mcepart)
return;
3632 if(!mcv0part)
return;
3634 for(Int_t i=0;i<11;i++){
3650 Double_t epx = mcepart->Px();
3651 Double_t epy = mcepart->Py();
3652 Double_t epz = mcepart->Pz();
3653 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
3654 Double_t v0px = mcv0part->Px();
3655 Double_t v0py = mcv0part->Py();
3656 Double_t v0pz = mcv0part->Pz();
3657 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
3659 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
3662 cont[0] = InvMassEleLambda;
3663 cont[1] = mcpart->Pt();
3666 cont2[0] = InvMassEleLambda;
3667 cont2[1] = mcepart->Pt();
3669 Double_t cont_eleptvseta[3];
3670 cont_eleptvseta[0] = mcepart->Pt();
3671 cont_eleptvseta[1] = mcepart->Eta();
3673 Double_t cont_eleptvslambdapt[3];
3674 cont_eleptvslambdapt[0] = mcepart->Pt();
3675 cont_eleptvslambdapt[1] = mcv0part->Pt();
3677 Double_t cont_eleptvslambdaptvslcpt[4];
3678 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
3679 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
3680 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
3683 contmc[0] = mcpart->Pt();
3684 contmc[1] = mcpart->Y();
3686 Double_t contmcele[3];
3687 contmcele[0] = mcepart->Pt();
3688 contmcele[1] = mcepart->Eta();
3690 Double_t contmc_withele[3];
3691 contmc_withele[0] = mcpart->Pt();
3692 contmc_withele[1] = mcpart->Y();
3693 contmc_withele[2] = mcepart->Y();
3697 esdcuts->GetEtaRange(etamin,etamax);
3704 if(fabs(mcepart->Eta())<0.8){
3711 if(fabs(mcepart->Eta())<etamax){
3715 if(InvMassEleLambda<2.3){
3721 if(fabs(mcpart->Y())<0.7){
3722 if(InvMassEleLambda<2.3){
3727 }
else if(decaytype==1){
3735 if(fabs(mcepart->Eta())<etamax){
3739 if(InvMassEleLambda<2.3){
3745 }
else if(decaytype==2){
3753 if(fabs(mcepart->Eta())<etamax){
3757 if(InvMassEleLambda<2.3){
3763 }
else if(decaytype==3){
3767 }
else if(decaytype==10){
3776 if(fabs(mcepart->Eta())<0.8){
3779 }
else if(decaytype==11){
3795 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
3799 TString * fCandidateVariableNames =
new TString[nVar];
3801 fCandidateVariableNames[ 0]=
"Centrality";
3802 fCandidateVariableNames[ 1]=
"ElePx";
3803 fCandidateVariableNames[ 2]=
"ElePy";
3804 fCandidateVariableNames[ 3]=
"ElePz";
3805 fCandidateVariableNames[ 4]=
"ElePdgCode";
3806 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
3807 fCandidateVariableNames[ 6]=
"RunNumber";
3808 fCandidateVariableNames[ 7]=
"EvNumber";
3810 for (Int_t ivar=0; ivar<nVar; ivar++) {
3821 if(!mcepart)
return;
3823 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
3828 Int_t semi_flag = 0;
3829 if(pdgarray_ele[0]>400&&pdgarray_ele[0]<450) semi_flag = 1;
3830 if(pdgarray_ele[0]>4000&&pdgarray_ele[0]<4500) semi_flag = 1;
3831 if(pdgarray_ele[0]>500&&pdgarray_ele[0]<550) semi_flag = 2;
3832 if(pdgarray_ele[0]>5000&&pdgarray_ele[0]<5500) semi_flag = 2;
3835 contmc[0] = mcepart->Pt();
3836 contmc[1] = mcepart->Eta();
3839 if(semi_flag==1 && !ele_from_bottom){
3842 if(semi_flag==1 && ele_from_bottom){
3849 Bool_t hfe_flag = kFALSE;
3850 Int_t labemother = mcepart->GetMother();
3851 Int_t pdgmotherele = -9999;
3853 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
3854 pdgmotherele = motherele->GetPdgCode();
3855 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
3860 if(!hfe_flag)
return;
3864 for(Int_t i=0;i<8;i++){
3888 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
3892 TString * fCandidateVariableNames =
new TString[nVar];
3894 fCandidateVariableNames[ 0]=
"Centrality";
3895 fCandidateVariableNames[ 1]=
"V0Px";
3896 fCandidateVariableNames[ 2]=
"V0Py";
3897 fCandidateVariableNames[ 3]=
"V0Pz";
3898 fCandidateVariableNames[ 4]=
"V0PdgCode";
3899 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
3900 fCandidateVariableNames[ 6]=
"RunNumber";
3901 fCandidateVariableNames[ 7]=
"EvNumber";
3903 for (Int_t ivar=0; ivar<nVar; ivar++) {
3914 if(!mcv0part)
return;
3916 for(Int_t i=0;i<8;i++){
3920 Bool_t hfv0_flag = kFALSE;
3921 Int_t labv0mother = mcv0part->GetMother();
3922 Int_t pdgmotherv0 = -9999;
3924 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3926 pdgmotherv0 = motherv0->GetPdgCode();
3927 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
3932 if(!hfv0_flag)
return;
3935 contmc[0] = mcv0part->Pt();
3936 contmc[1] = mcv0part->Eta();
3960 fCEvents =
new TH1F(
"fCEvents",
"conter",18,-0.5,17.5);
3962 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
3963 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
3964 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
3965 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
3966 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
3967 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
3968 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
3969 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
3970 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
3971 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
3972 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
3973 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
3974 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
3975 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
3977 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
3978 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
3981 fCEvents->GetYaxis()->SetTitle(
"counts");
3983 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
3985 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
3986 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
3987 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
3988 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
3989 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
3990 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
3992 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
3993 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
3994 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
3995 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
3996 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
3998 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
4017 Int_t bins_base[3]= {23,100 ,10};
4018 Double_t xmin_base[3]={1.1,0 ,0.00};
4019 Double_t xmax_base[3]={5.7,20. ,100};
4020 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
4096 Int_t bins_base_elept[3]= {10,100 ,10};
4097 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
4098 Double_t xmax_base_elept[3]={3.1,10. ,100};
4157 fHistoElePtRS =
new TH2F(
"fHistoElePtRS",
"",100,0,10,10,0,100);
4159 fHistoElePtWS =
new TH2F(
"fHistoElePtWS",
"",100,0.,10.,10,0,100);
4165 fHistoElePtMCS =
new TH2F(
"fHistoElePtMCS",
"",100,0,10,10,0,100);
4170 Int_t bins_eleptvseta[3]= {10,10 ,10};
4171 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4172 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4174 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4176 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4178 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4180 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4182 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4184 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4187 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4188 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4189 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4191 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4193 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4195 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4197 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4199 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4201 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4203 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4206 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4207 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4208 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4209 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4211 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4214 Int_t bins_allpt[4]= {10,10,20,20};
4215 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4216 Double_t xmax_allpt[4]={20.,20.,10.,10};
4227 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4228 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4229 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4231 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4233 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4235 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4237 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4239 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4241 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4307 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4313 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4315 fHistoBachPtMCS =
new TH1F(
"fHistoBachPtMCS",
"Bachelor p_{T}",100,0.0,5.0);
4319 fHistod0Bach =
new TH1F(
"fHistod0Bach",
"Bachelor d_{0}",100,-0.5,0.5);
4321 fHistoLambdaMassvsPt=
new TH2F(
"fHistoLambdaMassvsPt",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4323 fHistoLambdaMassvsPtMCS=
new TH2F(
"fHistoLambdaMassvsPtMCS",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4333 fHistoK0sMassvsPt=
new TH2F(
"fHistoK0sMassvsPt",
"K0s mass",100,0.497-0.05,0.497+0.05,20,0.,10.);
4335 fHistoLambdaPtvsDl=
new TH2F(
"fHistoLambdaPtvsDl",
"Lambda pt vs dl",20,0.,10.,20,0.,40.);
4348 Int_t bins_armenteros[3]= {20 ,100 ,200};
4349 Double_t xmin_armenteros[3]={0.,0. ,-1.};
4350 Double_t xmax_armenteros[3]={20.,1. ,1.};
4351 fHistoLambdaArmenteros =
new THnSparseF(
"fHistoLambdaArmenteros",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4353 fHistoLambdaArmenterosMCS =
new THnSparseF(
"fHistoLambdaArmenterosMCS",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4375 for(Int_t i=0;i<8;i++){
4384 Int_t bins_lcmcgen_withelea[3]= {20 ,100 ,100};
4385 Double_t xmin_lcmcgen_withelea[3]={0.,-5.,-5.};
4386 Double_t xmax_lcmcgen_withelea[3]={20.,5.,5.};
4387 fHistoLcMCGenWithEleA =
new THnSparseF(
"fHistoLcMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4392 fHistoLbMCGenWithEleA =
new THnSparseF(
"fHistoLbMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4397 Int_t bins_lbmcgen_withlc[3]= {50 ,100 ,100};
4398 Double_t xmin_lbmcgen_withlc[3]={0.,-5.,-5.};
4399 Double_t xmax_lbmcgen_withlc[3]={50.,5.,5.};
4400 fHistoLbMCGenWithLambdac =
new THnSparseF(
"fHistoLbMCGenWithLambdac",
"",3,bins_lbmcgen_withlc,xmin_lbmcgen_withlc,xmax_lbmcgen_withlc);
4403 Int_t bins_lcmcgen[3]= {100 ,20 ,10};
4404 Double_t xmin_lcmcgen[3]={0.,-1.0 ,0.0};
4405 Double_t xmax_lcmcgen[3]={20.,1.0 ,100};
4406 Int_t bins_lbmcgen[3]= {100 ,20 ,10};
4407 Double_t xmin_lbmcgen[3]={0.,-1.0 ,0.0};
4408 Double_t xmax_lbmcgen[3]={50.,1.0 ,100};
4409 fHistoLcMCGen =
new THnSparseF(
"fHistoLcMCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4411 fHistoLcMCGen1 =
new THnSparseF(
"fHistoLcMCGen1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4413 fHistoLcMCGen2 =
new THnSparseF(
"fHistoLcMCGen2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4415 fHistoLbMCGen =
new THnSparseF(
"fHistoLbMCGen",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4417 fHistoLbMCGen1 =
new THnSparseF(
"fHistoLbMCGen1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4419 fHistoLbMCGen2 =
new THnSparseF(
"fHistoLbMCGen2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4421 fHistoFeedDownXic0MCGen =
new THnSparseF(
"fHistoFeedDownXic0MCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4440 fHistoLcMCS =
new THnSparseF(
"fHistoLcMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4442 fHistoLcMCS1 =
new THnSparseF(
"fHistoLcMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4444 fHistoLcMCS2 =
new THnSparseF(
"fHistoLcMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4446 fHistoLbMCS =
new THnSparseF(
"fHistoLbMCS",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4448 fHistoLbMCS1 =
new THnSparseF(
"fHistoLbMCS1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4450 fHistoLbMCS2 =
new THnSparseF(
"fHistoLbMCS2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4452 fHistoFeedDownXic0MCS =
new THnSparseF(
"fHistoFeedDownXic0MCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4454 fHistoFeedDownXic0MCS1 =
new THnSparseF(
"fHistoFeedDownXic0MCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4456 fHistoFeedDownXic0MCS2 =
new THnSparseF(
"fHistoFeedDownXic0MCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4464 fHistoFeedDownSigmaMCS =
new THnSparseF(
"fHistoFeedDownSigmaMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4466 fHistoFeedDownSigmaMCS1 =
new THnSparseF(
"fHistoFeedDownSigmaMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4468 fHistoFeedDownSigmaMCS2 =
new THnSparseF(
"fHistoFeedDownSigmaMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4471 Int_t bins_elemcgen[3]= {100 ,20 ,10};
4472 Double_t xmin_elemcgen[3]={0.,-1.0 ,0.0};
4473 Double_t xmax_elemcgen[3]={10.,1.0 ,100};
4474 fHistoElectronMCGen =
new THnSparseF(
"fHistoElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4478 fHistoCharmElectronMCGen =
new THnSparseF(
"fHistoCharmElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4480 fHistoLcElectronMCGen =
new THnSparseF(
"fHistoLcElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4482 fHistoLcElectronMCGen1 =
new THnSparseF(
"fHistoLcElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4484 fHistoLcElectronMCGen2 =
new THnSparseF(
"fHistoLcElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4486 fHistoLbElectronMCGen =
new THnSparseF(
"fHistoLbElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4488 fHistoLbElectronMCGen1 =
new THnSparseF(
"fHistoLbElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4490 fHistoLbElectronMCGen2 =
new THnSparseF(
"fHistoLbElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4510 fHistoLcElectronMCS =
new THnSparseF(
"fHistoLcElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4512 fHistoLcElectronMCS1 =
new THnSparseF(
"fHistoLcElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4514 fHistoLcElectronMCS2 =
new THnSparseF(
"fHistoLcElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4516 fHistoLbElectronMCS =
new THnSparseF(
"fHistoLbElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4518 fHistoLbElectronMCS1 =
new THnSparseF(
"fHistoLbElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4520 fHistoLbElectronMCS2 =
new THnSparseF(
"fHistoLbElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4541 Int_t bins_lambdamcgen[3]= {50 ,20 ,10};
4542 Double_t xmin_lambdamcgen[3]={0.,-1.0 ,0.0};
4543 Double_t xmax_lambdamcgen[3]={10.,1.0 ,100};
4544 fHistoLambdaMCGen =
new THnSparseF(
"fHistoLambdaMCGen",
"",3,bins_lambdamcgen,xmin_lambdamcgen,xmax_lambdamcgen);
4547 Int_t bins_eleptvsv0dl[3]= {100 ,20 ,10};
4548 Double_t xmin_eleptvsv0dl[3]={0.,0. ,0.0};
4549 Double_t xmax_eleptvsv0dl[3]={10.,40. ,100};
4550 fHistoElePtvsV0dlRS =
new THnSparseF(
"fHistoElePtvsV0dlRS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4552 fHistoElePtvsV0dlRS1 =
new THnSparseF(
"fHistoElePtvsV0dlRS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4554 fHistoElePtvsV0dlRS2 =
new THnSparseF(
"fHistoElePtvsV0dlRS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4556 fHistoElePtvsV0dlRSSide =
new THnSparseF(
"fHistoElePtvsV0dlRSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4558 fHistoElePtvsV0dlRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4560 fHistoElePtvsV0dlRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4562 fHistoElePtvsV0dlRSMix =
new THnSparseF(
"fHistoElePtvsV0dlRSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4564 fHistoElePtvsV0dlRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4566 fHistoElePtvsV0dlRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4568 fHistoElePtvsV0dlWS =
new THnSparseF(
"fHistoElePtvsV0dlWS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4570 fHistoElePtvsV0dlWS1 =
new THnSparseF(
"fHistoElePtvsV0dlWS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4572 fHistoElePtvsV0dlWS2 =
new THnSparseF(
"fHistoElePtvsV0dlWS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4574 fHistoElePtvsV0dlWSSide =
new THnSparseF(
"fHistoElePtvsV0dlWSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4576 fHistoElePtvsV0dlWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4578 fHistoElePtvsV0dlWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4580 fHistoElePtvsV0dlWSMix =
new THnSparseF(
"fHistoElePtvsV0dlWSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4582 fHistoElePtvsV0dlWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4584 fHistoElePtvsV0dlWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4586 fHistoElePtvsV0dlMCS =
new THnSparseF(
"fHistoElePtvsV0dlMCS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4588 fHistoElePtvsV0dlMCS1 =
new THnSparseF(
"fHistoElePtvsV0dlMCS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4590 fHistoElePtvsV0dlMCS2 =
new THnSparseF(
"fHistoElePtvsV0dlMCS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4605 Int_t bins_eleptvsv0dca[3]= {100 ,20 ,10};
4606 Double_t xmin_eleptvsv0dca[3]={0.,0. ,0.0};
4607 Double_t xmax_eleptvsv0dca[3]={10.,1. ,100};
4608 fHistoElePtvsV0dcaRS =
new THnSparseF(
"fHistoElePtvsV0dcaRS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4610 fHistoElePtvsV0dcaRS1 =
new THnSparseF(
"fHistoElePtvsV0dcaRS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4612 fHistoElePtvsV0dcaRS2 =
new THnSparseF(
"fHistoElePtvsV0dcaRS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4614 fHistoElePtvsV0dcaRSSide =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4616 fHistoElePtvsV0dcaRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4618 fHistoElePtvsV0dcaRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4620 fHistoElePtvsV0dcaRSMix =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4622 fHistoElePtvsV0dcaRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4624 fHistoElePtvsV0dcaRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4626 fHistoElePtvsV0dcaWS =
new THnSparseF(
"fHistoElePtvsV0dcaWS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4628 fHistoElePtvsV0dcaWS1 =
new THnSparseF(
"fHistoElePtvsV0dcaWS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4630 fHistoElePtvsV0dcaWS2 =
new THnSparseF(
"fHistoElePtvsV0dcaWS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4632 fHistoElePtvsV0dcaWSSide =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4634 fHistoElePtvsV0dcaWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4636 fHistoElePtvsV0dcaWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4638 fHistoElePtvsV0dcaWSMix =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4640 fHistoElePtvsV0dcaWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4642 fHistoElePtvsV0dcaWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4644 fHistoElePtvsV0dcaMCS =
new THnSparseF(
"fHistoElePtvsV0dcaMCS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4646 fHistoElePtvsV0dcaMCS1 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4648 fHistoElePtvsV0dcaMCS2 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4663 Int_t bins_elelamptvsv0dl[4]= {100 ,20 ,100, 10};
4664 Double_t xmin_elelamptvsv0dl[4]={0.,0. ,-0.5, 0.0};
4665 Double_t xmax_elelamptvsv0dl[4]={20.,40.,0.5 ,100};
4666 fHistoEleLambdaPtvsV0dlRS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4668 fHistoEleLambdaPtvsV0dlRS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4670 fHistoEleLambdaPtvsV0dlRS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4672 fHistoEleLambdaPtvsV0dlRSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4674 fHistoEleLambdaPtvsV0dlRSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4676 fHistoEleLambdaPtvsV0dlRSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4678 fHistoEleLambdaPtvsV0dlRSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4680 fHistoEleLambdaPtvsV0dlRSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4682 fHistoEleLambdaPtvsV0dlRSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4684 fHistoEleLambdaPtvsV0dlWS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4686 fHistoEleLambdaPtvsV0dlWS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4688 fHistoEleLambdaPtvsV0dlWS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4690 fHistoEleLambdaPtvsV0dlWSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4692 fHistoEleLambdaPtvsV0dlWSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4694 fHistoEleLambdaPtvsV0dlWSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4696 fHistoEleLambdaPtvsV0dlWSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4698 fHistoEleLambdaPtvsV0dlWSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4700 fHistoEleLambdaPtvsV0dlWSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4702 fHistoEleLambdaPtvsV0dlMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4704 fHistoEleLambdaPtvsV0dlMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4706 fHistoEleLambdaPtvsV0dlMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4782 fHistoEleLambdaPtvsV0dlXibMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4784 fHistoEleLambdaPtvsV0dlXibMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4786 fHistoEleLambdaPtvsV0dlXibMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4789 fHistoEleLambdaPtvsV0dlRSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4791 fHistoEleLambdaPtvsV0dlRS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4793 fHistoEleLambdaPtvsV0dlRS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4807 fHistoEleLambdaPtvsV0dlWSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4809 fHistoEleLambdaPtvsV0dlWS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4811 fHistoEleLambdaPtvsV0dlWS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);