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 fHistoCorrelationVariablesvsEleLambdaPt(0),
492 fHistoCorrelationVariablesvsEleLambdaPtMix(0),
493 fHistoCorrelationVariablesvsEleLambdaPtMC(0),
494 fHistoCorrelationVariablesvsElePt(0),
495 fHistoCorrelationVariablesvsElePtMix(0),
496 fHistoCorrelationVariablesvsElePtMC(0),
497 fHistoCorrelationVariablesvsLambdaPt(0),
498 fHistoCorrelationVariablesvsLambdaPtMix(0),
499 fHistoCorrelationVariablesvsLambdaPtMC(0),
500 fHistoResponseElePt(0),
501 fHistoResponseElePt1(0),
502 fHistoResponseElePt2(0),
503 fHistoResponseEleLambdaPt(0),
504 fHistoResponseEleLambdaPt1(0),
505 fHistoResponseEleLambdaPt2(0),
506 fHistoResponseEleLambdaPtFeeddownXic0(0),
507 fHistoResponseEleLambdaPtFeeddownXic01(0),
508 fHistoResponseEleLambdaPtFeeddownXic02(0),
509 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
510 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
511 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
512 fHistoResponseEleLambdaPtFeeddownSigma(0),
513 fHistoResponseEleLambdaPtFeeddownSigma1(0),
514 fHistoResponseEleLambdaPtFeeddownSigma2(0),
515 fHistoResponseEleLambdaPtLambdab(0),
516 fHistoResponseElePtLambdab(0),
517 fHistoResponseLambdaPtLambdab(0),
518 fHistoResponseLambdaPt(0),
519 fHistoResponseLambdaPtFeeddownXic0(0),
520 fHistoResponseLambdaPtFeeddownXicPlus(0),
521 fHistoResponseLambdaPtFeeddownSigma(0),
522 fHistoResponseLambdaPtvsEleLambdaPt(0),
523 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
524 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
525 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
526 fHistoEleLambdaPtvsRapidityRS(0),
527 fHistoEleLambdaPtvsRapidityWS(0),
528 fHistoEleLambdaPtvsRapidityMCS(0),
530 fHistoElectronPi0Total(0),
531 fHistoElectronPi0Tag(0),
533 fHistoElectronEtaTotal(0),
534 fHistoElectronEtaTag(0),
536 fHistonEvtvsRunNumber(0),
537 fHistonElevsRunNumber(0),
538 fHistonLambdavsRunNumber(0),
539 fHistoMCEventType(0),
540 fHistoMCDeltaPhiccbar(0),
541 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
542 fHistodPhiSdEtaSElectronProtonR125RS(0),
543 fHistodPhiSdEtaSElectronProtonR125WS(0),
544 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
545 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
546 fHistodPhiSdEtaSElectronPionR125RS(0),
547 fHistodPhiSdEtaSElectronPionR125WS(0),
548 fHistodPhiSdEtaSElectronPionR125RSMix(0),
549 fHistodPhiSdEtaSElectronPionR125WSMix(0),
551 fMixWithoutConversionFlag(kFALSE),
552 fNumberOfEventsForMixing (5),
557 fEventInfo(new TObjString(
"")),
558 fElectronTracks(0x0),
565 fElectronCutVarsArray(0x0),
566 fV0CutVarsArray1(0x0),
567 fV0CutVarsArray2(0x0)
572 for(Int_t i=0;i<17;i++){
577 for(Int_t i=0;i<8;i++){
585 Bool_t writeVariableTree) :
586 AliAnalysisTaskSE(name),
595 fIsEventSelected(kFALSE),
596 fWriteVariableTree(writeVariableTree),
597 fWriteEachVariableTree(kFALSE),
598 fWriteMCVariableTree(kFALSE),
600 fEleVariablesTree(0),
603 fMCEleVariablesTree(0),
604 fMCV0VariablesTree(0),
605 fMCGenPairVariablesTree(0),
606 fCorrelationVariablesTree(0),
607 fReconstructPrimVert(kFALSE),
613 fCandidateVariables(),
614 fCandidateEleVariables(),
615 fCandidateV0Variables(),
616 fCandidateMCVariables(),
617 fCandidateMCEleVariables(),
618 fCandidateMCV0Variables(),
619 fCandidateMCGenPairVariables(),
620 fCorrelationVariables(),
628 fUseCentralityV0M(kFALSE),
631 fMCDoPairAnalysis(kFALSE),
632 fHistoEleLambdaMass(0),
633 fHistoEleLambdaMassRS(0),
634 fHistoEleLambdaMassRS1(0),
635 fHistoEleLambdaMassRS2(0),
636 fHistoEleLambdaMassWS(0),
637 fHistoEleLambdaMassWS1(0),
638 fHistoEleLambdaMassWS2(0),
639 fHistoEleLambdaMassRSMix(0),
640 fHistoEleLambdaMassRSMix1(0),
641 fHistoEleLambdaMassRSMix2(0),
642 fHistoEleLambdaMassWSMix(0),
643 fHistoEleLambdaMassWSMix1(0),
644 fHistoEleLambdaMassWSMix2(0),
645 fHistoEleLambdaMassRSSide(0),
646 fHistoEleLambdaMassRSSide1(0),
647 fHistoEleLambdaMassRSSide2(0),
648 fHistoEleLambdaMassWSSide(0),
649 fHistoEleLambdaMassWSSide1(0),
650 fHistoEleLambdaMassWSSide2(0),
651 fHistoEleLambdaMassRSAway(0),
652 fHistoEleLambdaMassRS1Away(0),
653 fHistoEleLambdaMassRS2Away(0),
654 fHistoEleLambdaMassWSAway(0),
655 fHistoEleLambdaMassWS1Away(0),
656 fHistoEleLambdaMassWS2Away(0),
657 fHistoEleLambdaMassRSMixAway(0),
658 fHistoEleLambdaMassRSMix1Away(0),
659 fHistoEleLambdaMassRSMix2Away(0),
660 fHistoEleLambdaMassWSMixAway(0),
661 fHistoEleLambdaMassWSMix1Away(0),
662 fHistoEleLambdaMassWSMix2Away(0),
663 fHistoEleLambdaMassRSSideAway(0),
664 fHistoEleLambdaMassRSSide1Away(0),
665 fHistoEleLambdaMassRSSide2Away(0),
666 fHistoEleLambdaMassWSSideAway(0),
667 fHistoEleLambdaMassWSSide1Away(0),
668 fHistoEleLambdaMassWSSide2Away(0),
669 fHistoEleLambdaMassvsElePtRS(0),
670 fHistoEleLambdaMassvsElePtWS(0),
671 fHistoEleLambdaMassvsElePtRSMix(0),
672 fHistoEleLambdaMassvsElePtWSMix(0),
673 fHistoEleLambdaMassvsElePtRSSide(0),
674 fHistoEleLambdaMassvsElePtWSSide(0),
675 fHistoEleLambdaMassvsElePtRS1(0),
676 fHistoEleLambdaMassvsElePtWS1(0),
677 fHistoEleLambdaMassvsElePtRSMix1(0),
678 fHistoEleLambdaMassvsElePtWSMix1(0),
679 fHistoEleLambdaMassvsElePtRSSide1(0),
680 fHistoEleLambdaMassvsElePtWSSide1(0),
681 fHistoEleLambdaMassvsElePtRS2(0),
682 fHistoEleLambdaMassvsElePtWS2(0),
683 fHistoEleLambdaMassvsElePtRSMix2(0),
684 fHistoEleLambdaMassvsElePtWSMix2(0),
685 fHistoEleLambdaMassvsElePtRSSide2(0),
686 fHistoEleLambdaMassvsElePtWSSide2(0),
691 fHistoEleLambdaMassMCS(0),
692 fHistoEleLambdaMassLambdabMCS(0),
693 fHistoEleLambdaMassXibMCS(0),
694 fHistoEleLambdaMassMCGen(0),
695 fHistoEleLambdaMassvsElePtMCS(0),
696 fHistoEleLambdaMassvsElePtMCGen(0),
697 fHistoEleLambdaMassvsElePtMCS1(0),
698 fHistoEleLambdaMassvsElePtMCGen1(0),
699 fHistoEleLambdaMassvsElePtMCS2(0),
700 fHistoEleLambdaMassvsElePtMCGen2(0),
703 fHistoElePtvsEtaRS(0),
704 fHistoElePtvsEtaWS(0),
705 fHistoElePtvsEtaRSMix(0),
706 fHistoElePtvsEtaWSMix(0),
707 fHistoElePtvsEtaMCS(0),
708 fHistoElePtvsEtaMCGen(0),
709 fHistoElePtvsLambdaPtRS(0),
710 fHistoElePtvsLambdaPtWS(0),
711 fHistoElePtvsLambdaPtRSMix(0),
712 fHistoElePtvsLambdaPtWSMix(0),
713 fHistoElePtvsLambdaPtMCS(0),
714 fHistoElePtvsLambdaPtvsLcPtMCS(0),
715 fHistoElePtvsLambdaPtMCGen(0),
716 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
717 fHistoElePtvsLambdaPtMCLcGen(0),
718 fHistoElePtvsd0RS(0),
719 fHistoElePtvsd0WS(0),
720 fHistoElePtvsd0RSMix(0),
721 fHistoElePtvsd0WSMix(0),
722 fHistoElePtvsd0MCS(0),
723 fHistoElePtvsd0PromptMCS(0),
724 fHistoElePtvsd0BFeeddownMCS(0),
725 fHistoEleLambdaMassFeeddownXic0MCS(0),
726 fHistoEleLambdaMassFeeddownXic0MCGen(0),
727 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
728 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
729 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
730 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
731 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
732 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
733 fHistoElePtFeeddownXic0MCS(0),
734 fHistoElePtFeeddownXic0MCGen(0),
735 fHistoElePtvsEtaFeeddownXic0MCS(0),
736 fHistoElePtvsEtaFeeddownXic0MCGen(0),
737 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
738 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
739 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
740 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
741 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
742 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
743 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
744 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
745 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
746 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
747 fHistoElePtFeeddownXicPlusMCS(0),
748 fHistoElePtFeeddownXicPlusMCGen(0),
749 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
750 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
751 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
752 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
755 fHistoBachPtMCGen(0),
757 fHistoLambdaMassvsPt(0),
758 fHistoLambdaMassvsPtMCS(0),
759 fHistoLambdaMassvsPtMCGen(0),
760 fHistoLambdaPtvsEta(0),
761 fHistoLambdaPtvsEtaMCS(0),
762 fHistoLambdaPtvsEtaMCGen(0),
763 fHistoLambdaArmenteros(0),
764 fHistoLambdaArmenterosMCS(0),
765 fHistoLambdaPtvsDl(0),
766 fHistoLambdaPtvsDlSide(0),
767 fHistoLambdaPtvsDlMCS(0),
768 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
769 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
770 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
771 fHistoK0sMassvsPt(0),
772 fHistoElectronTPCPID(0),
773 fHistoElectronTOFPID(0),
774 fHistoElectronTPCSelPID(0),
775 fHistoElectronTOFSelPID(0),
776 fHistoElectronTPCPIDSelTOF(0),
777 fHistoElectronTPCPIDSelTOFSmallEta(0),
778 fHistoElectronTPCPIDSelTOFLargeEta(0),
779 fHistoMassConversionsMin(0),
780 fHistoMassConversionsSameSignMin(0),
781 fHistoElectronQovPtvsPhi(0),
782 fHistoLambdaQovPtvsPhi(0),
789 fHistoLcMCGenWithEleA(0),
790 fHistoLcMCGenWithEleB(0),
791 fHistoLbMCGenWithEleA(0),
792 fHistoLbMCGenWithEleB(0),
793 fHistoLbMCGenWithLambdac(0),
800 fHistoFeedDownXic0MCGen(0),
801 fHistoFeedDownXic0MCGen1(0),
802 fHistoFeedDownXic0MCGen2(0),
803 fHistoFeedDownXic0MCS(0),
804 fHistoFeedDownXic0MCS1(0),
805 fHistoFeedDownXic0MCS2(0),
806 fHistoFeedDownXicPlusMCGen(0),
807 fHistoFeedDownXicPlusMCGen1(0),
808 fHistoFeedDownXicPlusMCGen2(0),
809 fHistoFeedDownXicPlusMCS(0),
810 fHistoFeedDownXicPlusMCS1(0),
811 fHistoFeedDownXicPlusMCS2(0),
812 fHistoFeedDownSigmaMCGen(0),
813 fHistoFeedDownSigmaMCGen1(0),
814 fHistoFeedDownSigmaMCGen2(0),
815 fHistoFeedDownSigmaMCS(0),
816 fHistoFeedDownSigmaMCS1(0),
817 fHistoFeedDownSigmaMCS2(0),
818 fHistoLcElectronMCGen(0),
819 fHistoLcElectronMCGen1(0),
820 fHistoLcElectronMCGen2(0),
821 fHistoLcElectronMCS(0),
822 fHistoLcElectronMCS1(0),
823 fHistoLcElectronMCS2(0),
824 fHistoLbElectronMCGen(0),
825 fHistoLbElectronMCGen1(0),
826 fHistoLbElectronMCGen2(0),
827 fHistoLbElectronMCS(0),
828 fHistoLbElectronMCS1(0),
829 fHistoLbElectronMCS2(0),
830 fHistoElectronFeedDownXic0MCGen(0),
831 fHistoElectronFeedDownXic0MCGen1(0),
832 fHistoElectronFeedDownXic0MCGen2(0),
833 fHistoElectronFeedDownXic0MCS(0),
834 fHistoElectronFeedDownXic0MCS1(0),
835 fHistoElectronFeedDownXic0MCS2(0),
836 fHistoElectronFeedDownXicPlusMCGen(0),
837 fHistoElectronFeedDownXicPlusMCGen1(0),
838 fHistoElectronFeedDownXicPlusMCGen2(0),
839 fHistoElectronFeedDownXicPlusMCS(0),
840 fHistoElectronFeedDownXicPlusMCS1(0),
841 fHistoElectronFeedDownXicPlusMCS2(0),
842 fHistoElectronFeedDownXibMCGen(0),
843 fHistoElectronFeedDownXibMCGen1(0),
844 fHistoElectronFeedDownXibMCGen2(0),
845 fHistoElectronFeedDownXibMCS(0),
846 fHistoElectronFeedDownXibMCS1(0),
847 fHistoElectronFeedDownXibMCS2(0),
848 fHistoElectronMCGen(0),
849 fHistoBottomElectronMCGen(0),
850 fHistoCharmElectronMCGen(0),
851 fHistoLambdaMCGen(0),
852 fHistoElePtvsV0dlRS(0),
853 fHistoElePtvsV0dlRS1(0),
854 fHistoElePtvsV0dlRS2(0),
855 fHistoElePtvsV0dlRSSide(0),
856 fHistoElePtvsV0dlRSSide1(0),
857 fHistoElePtvsV0dlRSSide2(0),
858 fHistoElePtvsV0dlRSMix(0),
859 fHistoElePtvsV0dlRSMix1(0),
860 fHistoElePtvsV0dlRSMix2(0),
861 fHistoElePtvsV0dlWS(0),
862 fHistoElePtvsV0dlWS1(0),
863 fHistoElePtvsV0dlWS2(0),
864 fHistoElePtvsV0dlWSSide(0),
865 fHistoElePtvsV0dlWSSide1(0),
866 fHistoElePtvsV0dlWSSide2(0),
867 fHistoElePtvsV0dlWSMix(0),
868 fHistoElePtvsV0dlWSMix1(0),
869 fHistoElePtvsV0dlWSMix2(0),
870 fHistoElePtvsV0dlMCS(0),
871 fHistoElePtvsV0dlMCS1(0),
872 fHistoElePtvsV0dlMCS2(0),
873 fHistoElePtvsV0dlFeedDownXic0MCS(0),
874 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
875 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
876 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
877 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
878 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
879 fHistoElePtvsV0dcaRS(0),
880 fHistoElePtvsV0dcaRS1(0),
881 fHistoElePtvsV0dcaRS2(0),
882 fHistoElePtvsV0dcaRSSide(0),
883 fHistoElePtvsV0dcaRSSide1(0),
884 fHistoElePtvsV0dcaRSSide2(0),
885 fHistoElePtvsV0dcaRSMix(0),
886 fHistoElePtvsV0dcaRSMix1(0),
887 fHistoElePtvsV0dcaRSMix2(0),
888 fHistoElePtvsV0dcaWS(0),
889 fHistoElePtvsV0dcaWS1(0),
890 fHistoElePtvsV0dcaWS2(0),
891 fHistoElePtvsV0dcaWSSide(0),
892 fHistoElePtvsV0dcaWSSide1(0),
893 fHistoElePtvsV0dcaWSSide2(0),
894 fHistoElePtvsV0dcaWSMix(0),
895 fHistoElePtvsV0dcaWSMix1(0),
896 fHistoElePtvsV0dcaWSMix2(0),
897 fHistoElePtvsV0dcaMCS(0),
898 fHistoElePtvsV0dcaMCS1(0),
899 fHistoElePtvsV0dcaMCS2(0),
900 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
901 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
902 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
903 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
904 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
905 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
906 fHistoEleLambdaPtvsV0dlRS(0),
907 fHistoEleLambdaPtvsV0dlRS1(0),
908 fHistoEleLambdaPtvsV0dlRS2(0),
909 fHistoEleLambdaPtvsV0dlRSSide(0),
910 fHistoEleLambdaPtvsV0dlRSSide1(0),
911 fHistoEleLambdaPtvsV0dlRSSide2(0),
912 fHistoEleLambdaPtvsV0dlRSMix(0),
913 fHistoEleLambdaPtvsV0dlRSMix1(0),
914 fHistoEleLambdaPtvsV0dlRSMix2(0),
915 fHistoEleLambdaPtvsV0dlWS(0),
916 fHistoEleLambdaPtvsV0dlWS1(0),
917 fHistoEleLambdaPtvsV0dlWS2(0),
918 fHistoEleLambdaPtvsV0dlWSSide(0),
919 fHistoEleLambdaPtvsV0dlWSSide1(0),
920 fHistoEleLambdaPtvsV0dlWSSide2(0),
921 fHistoEleLambdaPtvsV0dlWSMix(0),
922 fHistoEleLambdaPtvsV0dlWSMix1(0),
923 fHistoEleLambdaPtvsV0dlWSMix2(0),
924 fHistoEleLambdaPtvsV0dlMCS(0),
925 fHistoEleLambdaPtvsV0dlMCS1(0),
926 fHistoEleLambdaPtvsV0dlMCS2(0),
927 fHistoEleLambdaPtvsV0dlPromptMCS(0),
928 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
929 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
930 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
931 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
932 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
933 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
934 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
935 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
936 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
937 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
938 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
939 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
940 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
941 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
942 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
943 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
944 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
945 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
946 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
947 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
948 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
949 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
950 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
951 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
952 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
953 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
954 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
955 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
956 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
957 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
958 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
959 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
960 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
961 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
962 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
963 fHistoEleLambdaPtvsV0dlXibMCS(0),
964 fHistoEleLambdaPtvsV0dlXibMCS1(0),
965 fHistoEleLambdaPtvsV0dlXibMCS2(0),
966 fHistoEleLambdaPtvsV0dlRSAway(0),
967 fHistoEleLambdaPtvsV0dlRS1Away(0),
968 fHistoEleLambdaPtvsV0dlRS2Away(0),
969 fHistoEleLambdaPtvsV0dlRSSideAway(0),
970 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
971 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
972 fHistoEleLambdaPtvsV0dlRSMixAway(0),
973 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
974 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
975 fHistoEleLambdaPtvsV0dlWSAway(0),
976 fHistoEleLambdaPtvsV0dlWS1Away(0),
977 fHistoEleLambdaPtvsV0dlWS2Away(0),
978 fHistoEleLambdaPtvsV0dlWSSideAway(0),
979 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
980 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
981 fHistoEleLambdaPtvsV0dlWSMixAway(0),
982 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
983 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
984 fHistoCorrelationVariablesvsEleLambdaPt(0),
985 fHistoCorrelationVariablesvsEleLambdaPtMix(0),
986 fHistoCorrelationVariablesvsEleLambdaPtMC(0),
987 fHistoCorrelationVariablesvsElePt(0),
988 fHistoCorrelationVariablesvsElePtMix(0),
989 fHistoCorrelationVariablesvsElePtMC(0),
990 fHistoCorrelationVariablesvsLambdaPt(0),
991 fHistoCorrelationVariablesvsLambdaPtMix(0),
992 fHistoCorrelationVariablesvsLambdaPtMC(0),
993 fHistoResponseElePt(0),
994 fHistoResponseElePt1(0),
995 fHistoResponseElePt2(0),
996 fHistoResponseEleLambdaPt(0),
997 fHistoResponseEleLambdaPt1(0),
998 fHistoResponseEleLambdaPt2(0),
999 fHistoResponseEleLambdaPtFeeddownXic0(0),
1000 fHistoResponseEleLambdaPtFeeddownXic01(0),
1001 fHistoResponseEleLambdaPtFeeddownXic02(0),
1002 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
1003 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
1004 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
1005 fHistoResponseEleLambdaPtFeeddownSigma(0),
1006 fHistoResponseEleLambdaPtFeeddownSigma1(0),
1007 fHistoResponseEleLambdaPtFeeddownSigma2(0),
1008 fHistoResponseEleLambdaPtLambdab(0),
1009 fHistoResponseElePtLambdab(0),
1010 fHistoResponseLambdaPtLambdab(0),
1011 fHistoResponseLambdaPt(0),
1012 fHistoResponseLambdaPtFeeddownXic0(0),
1013 fHistoResponseLambdaPtFeeddownXicPlus(0),
1014 fHistoResponseLambdaPtFeeddownSigma(0),
1015 fHistoResponseLambdaPtvsEleLambdaPt(0),
1016 fHistoResponseLambdaPtvsEleLambdaPtLambdab(0),
1017 fHistoResponseMcGenLambdabPtvsLambdacPt(0),
1018 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
1019 fHistoEleLambdaPtvsRapidityRS(0),
1020 fHistoEleLambdaPtvsRapidityWS(0),
1021 fHistoEleLambdaPtvsRapidityMCS(0),
1023 fHistoElectronPi0Total(0),
1024 fHistoElectronPi0Tag(0),
1026 fHistoElectronEtaTotal(0),
1027 fHistoElectronEtaTag(0),
1029 fHistonEvtvsRunNumber(0),
1030 fHistonElevsRunNumber(0),
1031 fHistonLambdavsRunNumber(0),
1032 fHistoMCEventType(0),
1033 fHistoMCDeltaPhiccbar(0),
1034 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
1035 fHistodPhiSdEtaSElectronProtonR125RS(0),
1036 fHistodPhiSdEtaSElectronProtonR125WS(0),
1037 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
1038 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
1039 fHistodPhiSdEtaSElectronPionR125RS(0),
1040 fHistodPhiSdEtaSElectronPionR125WS(0),
1041 fHistodPhiSdEtaSElectronPionR125RSMix(0),
1042 fHistodPhiSdEtaSElectronPionR125WSMix(0),
1044 fMixWithoutConversionFlag(kFALSE),
1045 fNumberOfEventsForMixing (5),
1050 fEventInfo(new TObjString(
"")),
1051 fElectronTracks(0x0),
1058 fElectronCutVarsArray(0x0),
1059 fV0CutVarsArray1(0x0),
1060 fV0CutVarsArray2(0x0)
1065 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
1067 for(Int_t i=0;i<17;i++){
1072 for(Int_t i=0;i<8;i++){
1076 DefineOutput(1,TList::Class());
1077 DefineOutput(2,TList::Class());
1078 DefineOutput(3,TList::Class());
1079 DefineOutput(4,TTree::Class());
1080 DefineOutput(5,TTree::Class());
1081 DefineOutput(6,TTree::Class());
1082 DefineOutput(7,TTree::Class());
1083 DefineOutput(8,AliNormalizationCounter::Class());
1084 DefineOutput(9,TTree::Class());
1085 DefineOutput(10,TTree::Class());
1086 DefineOutput(11,TTree::Class());
1087 DefineOutput(12,TTree::Class());
1095 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1190 if (fDebug > 1) AliInfo(
"Init");
1209 AliError(
"NO EVENT FOUND!");
1212 AliAODEvent* aodEvent =
dynamic_cast<AliAODEvent*
>(fInputEvent);
1220 fBzkG = (Double_t)aodEvent->GetMagneticField();
1221 AliKFParticle::SetField(
fBzkG);
1222 if (TMath::Abs(
fBzkG)<0.001) {
1233 TClonesArray *mcArray = 0;
1234 AliAODMCHeader *mcHeader=0;
1237 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1239 AliError(
"Could not find Monte-Carlo in AOD");
1245 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1247 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1252 Double_t zMCVertex = mcHeader->GetVtxZ();
1254 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1268 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1271 Double_t pos[3],cov[6];
1273 fVtx1->GetCovarianceMatrix(cov);
1274 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1278 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1285 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1286 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1287 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1288 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1289 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1303 AliCentrality *cent = aodEvent->GetCentrality();
1304 fCentrality = cent->GetCentralityPercentile(
"V0M");
1308 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1315 Int_t runnumber_offset = 0;
1316 Int_t runnumber = aodEvent->GetRunNumber();
1317 if(runnumber<=131000&&runnumber>=114000){
1318 runnumber_offset = 114000;
1319 }
else if(runnumber<=196000&&runnumber>=195000){
1320 runnumber_offset = 195000;
1321 }
else if(runnumber<=170593&&runnumber>=167902){
1322 runnumber_offset = 167902;
1367 AliAnalysisTaskSE::Terminate();
1369 fOutput =
dynamic_cast<TList*
> (GetOutputData(1));
1371 AliError(
"fOutput not available");
1375 fOutputAll =
dynamic_cast<TList*
> (GetOutputData(3));
1377 AliError(
"fOutputAll not available");
1432 TString normName=
"NormalizationCounter";
1433 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1434 if(cont)normName=(TString)cont->GetName();
1456 fEventBuffer[i]=
new TTree(Form(
"EventBuffer_%d",i),
"Temporary buffer for event mixing");
1480 AliAODEvent *aodEvent, TClonesArray *mcArray
1487 if(fElectronTracks) fElectronTracks->Delete();
1488 if(fV0Tracks1) fV0Tracks1->Delete();
1489 if(fV0Tracks2) fV0Tracks2->Delete();
1490 if(fElectronCutVarsArray) fElectronCutVarsArray->Delete();
1491 if(fV0CutVarsArray1) fV0CutVarsArray1->Delete();
1492 if(fV0CutVarsArray2) fV0CutVarsArray2->Delete();
1493 fV0dlArray1.clear();
1494 fV0dlArray2.clear();
1495 fV0dcaArray1.clear();
1496 fV0dcaArray2.clear();
1499 ResetGlobalTrackReference();
1501 for (Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1503 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1504 if (!track)
continue;
1507 StoreGlobalTrackReference(track,iTrack);
1514 Int_t nV0s= aodEvent->GetNumberOfV0s();
1515 Int_t nTracks= aodEvent->GetNumberOfTracks();
1517 Bool_t seleTrkFlags[nTracks];
1519 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1521 Bool_t seleV0Flags[nV0s];
1523 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1525 Int_t runnumber_offset = 0;
1526 Int_t runnumber = aodEvent->GetRunNumber();
1527 if(runnumber<=131000&&runnumber>=114000){
1528 runnumber_offset = 114000;
1529 }
else if(runnumber<=196000&&runnumber>=195000){
1530 runnumber_offset = 195000;
1531 }
else if(runnumber<=170593&&runnumber>=167902){
1532 runnumber_offset = 167902;
1534 fHistonElevsRunNumber->Fill(runnumber-runnumber_offset,nSeleTrks);
1535 fHistonLambdavsRunNumber->Fill(runnumber-runnumber_offset,nSeleV0);
1540 for (Int_t iv0 = 0; iv0<nV0s; iv0++) {
1541 if(!seleV0Flags[iv0])
continue;
1542 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1545 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1546 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1551 for (Int_t itrk = 0; itrk<nTracks; itrk++) {
1552 if(!seleTrkFlags[itrk])
continue;
1553 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1561 AliAODTrack *trkpid = 0;
1562 if(fAnalCuts->GetProdAODFilterBit()==7){
1563 trkpid = fGTI[-trk->GetID()-1];
1568 Int_t cpid = cptrack->GetID();
1569 Int_t cnid = cntrack->GetID();
1570 Int_t lpid = trkpid->GetID();
1571 if((cpid==lpid)||(cnid==lpid))
continue;
1573 AliAODVertex *secVert = ReconstructSecondaryVertex(v0,trk,aodEvent);
1574 if(!secVert)
continue;
1581 FillROOTObjects(elobj,v0,trk,trkpid,aodEvent,mcArray);
1583 elobj->GetSecondaryVtx()->RemoveDaughters();
1585 delete elobj;elobj=NULL;
1591 fEventInfo->SetString(Form(
"Ev%d_esd%d_E%d_V%d",AliAnalysisManager::GetAnalysisManager()->GetNcalls(),((AliAODHeader*)aodEvent->GetHeader())->GetEventNumberESDFile(),fElectronTracks->GetEntries(),fV0Tracks1->GetEntries()+fV0Tracks2->GetEntries()));
1592 Int_t ind=GetPoolIndex(fVtxZ,fCentrality);
1593 if(ind>=0 && ind<fNOfPools){
1594 if(fEventBuffer[ind]->GetEntries() >= fNumberOfEventsForMixing){
1595 DoEventMixingWithPools(ind);
1596 if(fEventBuffer[ind]->GetEntries() >= 20*fNumberOfEventsForMixing){
1600 fEventBuffer[ind]->Fill();
1611 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1612 fVariablesTree =
new TTree(nameoutput,
"Candidates variables tree");
1615 TString * fCandidateVariableNames =
new TString[
nVar];
1617 fCandidateVariableNames[ 0]=
"Centrality";
1618 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1619 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1620 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1621 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1622 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1623 fCandidateVariableNames[ 6]=
"ElePx";
1624 fCandidateVariableNames[ 7]=
"ElePy";
1625 fCandidateVariableNames[ 8]=
"ElePz";
1626 fCandidateVariableNames[ 9]=
"V0Px";
1627 fCandidateVariableNames[10]=
"V0Py";
1628 fCandidateVariableNames[11]=
"V0Pz";
1629 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1630 fCandidateVariableNames[13]=
"MassLambda";
1631 fCandidateVariableNames[14]=
"MassAntiLambda";
1632 fCandidateVariableNames[15]=
"Eled0";
1633 fCandidateVariableNames[16]=
"V0d0";
1634 fCandidateVariableNames[17]=
"nSigmaTPCele";
1635 fCandidateVariableNames[18]=
"nSigmaTOFele";
1636 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1637 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1638 fCandidateVariableNames[21]=
"EleCharge";
1639 fCandidateVariableNames[22]=
"ProtonPx";
1640 fCandidateVariableNames[23]=
"ProtonPy";
1641 fCandidateVariableNames[24]=
"ProtonPz";
1642 fCandidateVariableNames[25]=
"PiPx";
1643 fCandidateVariableNames[26]=
"PiPy";
1644 fCandidateVariableNames[27]=
"PiPz";
1645 fCandidateVariableNames[28]=
"mcpdglc";
1646 fCandidateVariableNames[29]=
"mclablc";
1647 fCandidateVariableNames[30]=
"mcpdgmomele";
1648 fCandidateVariableNames[31]=
"mcpdgmomv0";
1649 fCandidateVariableNames[32]=
"Mixing";
1650 fCandidateVariableNames[33]=
"mcpdgele";
1651 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1652 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1653 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1654 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1655 fCandidateVariableNames[38]=
"MassK0Short";
1656 fCandidateVariableNames[39]=
"mcpdggrmomele";
1657 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1658 fCandidateVariableNames[41]=
"mcngenele";
1659 fCandidateVariableNames[42]=
"mcngenv0";
1660 fCandidateVariableNames[43]=
"mclcpx";
1661 fCandidateVariableNames[44]=
"mclcpy";
1662 fCandidateVariableNames[45]=
"mclcpz";
1663 fCandidateVariableNames[46]=
"mcelepx";
1664 fCandidateVariableNames[47]=
"mcelepy";
1665 fCandidateVariableNames[48]=
"mcelepz";
1666 fCandidateVariableNames[49]=
"mcv0px";
1667 fCandidateVariableNames[50]=
"mcv0py";
1668 fCandidateVariableNames[51]=
"mcv0pz";
1669 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1670 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1671 fCandidateVariableNames[54]=
"PrimVertx";
1672 fCandidateVariableNames[55]=
"PrimVerty";
1673 fCandidateVariableNames[56]=
"PrimVertz";
1674 fCandidateVariableNames[57]=
"V0Vertx";
1675 fCandidateVariableNames[58]=
"V0Verty";
1676 fCandidateVariableNames[59]=
"V0Vertz";
1678 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1679 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1680 fCandidateVariableNames[62]=
"DcaV0daughters";
1681 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1682 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1683 fCandidateVariableNames[65]=
"MassK0Short2";
1685 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1686 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1688 fCandidateVariableNames[68]=
"EleITSMatch";
1689 fCandidateVariableNames[69]=
"V0PosITSMatch";
1690 fCandidateVariableNames[70]=
"V0NegITSMatch";
1691 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1692 fCandidateVariableNames[72]=
"mcpdgv0";
1694 fCandidateVariableNames[73]=
"mcv0pospx";
1695 fCandidateVariableNames[74]=
"mcv0pospy";
1696 fCandidateVariableNames[75]=
"mcv0pospz";
1697 fCandidateVariableNames[76]=
"mcpdgv0pos";
1698 fCandidateVariableNames[77]=
"mcv0negpx";
1699 fCandidateVariableNames[78]=
"mcv0negpy";
1700 fCandidateVariableNames[79]=
"mcv0negpz";
1701 fCandidateVariableNames[80]=
"mcpdgv0neg";
1703 fCandidateVariableNames[81]=
"ProtonPxPV";
1704 fCandidateVariableNames[82]=
"ProtonPyPV";
1705 fCandidateVariableNames[83]=
"ProtonPzPV";
1706 fCandidateVariableNames[84]=
"PiPxPV";
1707 fCandidateVariableNames[85]=
"PiPyPV";
1708 fCandidateVariableNames[86]=
"PiPzPV";
1710 fCandidateVariableNames[87]=
"mcv0vertx";
1711 fCandidateVariableNames[88]=
"mcv0verty";
1712 fCandidateVariableNames[89]=
"mcv0vertz";
1714 fCandidateVariableNames[90]=
"OnFlyStatus";
1715 fCandidateVariableNames[91]=
"MagneticField";
1716 fCandidateVariableNames[92]=
"EvNumber";
1717 fCandidateVariableNames[93]=
"RunNumber";
1720 for (Int_t ivar=0; ivar<
nVar; ivar++) {
1737 for(Int_t i=0;i<94;i++)
1742 Bool_t anti_lambda_flag = kFALSE;
1745 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1746 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1752 Double_t v0px = elobj->PxProng(1);
1753 Double_t v0py = elobj->PyProng(1);
1754 Double_t v0pz = elobj->PzProng(1);
1755 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1756 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1757 Double_t epx = elobj->PxProng(0);
1758 Double_t epy = elobj->PyProng(0);
1759 Double_t epz = elobj->PzProng(0);
1760 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1761 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1762 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
1763 Double_t Esum = Ee + Ev0;
1765 Double_t uxe = epx/mome;
1766 Double_t uye = epy/mome;
1767 Double_t uze = epz/mome;
1768 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
1769 Double_t pxv_flip = v0px + lf * uxe;
1770 Double_t pyv_flip = v0py + lf * uye;
1771 Double_t pzv_flip = v0pz + lf * uze;
1772 Double_t pxsum_flip = epx + pxv_flip;
1773 Double_t pysum_flip = epy + pyv_flip;
1774 Double_t pzsum_flip = epz + pzv_flip;
1775 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
1776 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
1778 Double_t minmass_ee = 9999.;
1780 Double_t minmasslike_ee = 9999.;
1784 UInt_t pdgdg[2]={11,3122};
1802 Double_t nSigmaTPCele=-9999.;
1803 Double_t nSigmaTOFele=-9999.;
1804 Double_t nSigmaTPCv0pr=-9999.;
1805 Double_t nSigmaTOFv0pr=-9999.;
1806 Double_t nSigmaTPCv0pi=-9999.;
1807 Double_t nSigmaTOFv0pi=-9999.;
1818 if(anti_lambda_flag){
1836 if(anti_lambda_flag){
1852 AliAODMCParticle *mclc = 0;
1853 AliAODMCParticle *mcele = 0;
1854 AliAODMCParticle *mcv0 = 0;
1856 Int_t mcpdgele_array[100];
1857 Int_t mcpdgv0_array[100];
1858 Int_t mclabelele_array[100];
1859 Int_t mclabelv0_array[100];
1860 Int_t mcngen_ele=-9999;
1861 Int_t mcngen_v0=-9999;
1864 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
1865 if(mclabelele_array[0]>=0){
1866 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
1873 if(mclabelv0_array[0]>=0){
1874 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
1882 Int_t pdgdgv0[2]={211,211};
1883 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
1885 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
1895 Int_t labPos = cptrack->GetLabel();
1897 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
1906 Int_t labNeg = cntrack->GetLabel();
1908 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
1916 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
1959 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
1960 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
1961 if(!anti_lambda_flag){
1969 Double_t posVtx[3] = {0.,0.,0.};
1970 fVtx1->GetXYZ(posVtx);
1972 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
1982 if(!anti_lambda_flag){
2007 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
2008 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
2012 cont[0] = elobj->InvMass(2,pdgdg);
2013 cont[1] = elobj->Pt();
2017 Double_t cont_flip[3];
2018 cont_flip[0] = melam_flip;
2019 cont_flip[1] = ptelam_flip;
2023 cont2[0] = elobj->InvMass(2,pdgdg);
2024 cont2[1] = trk->Pt();
2026 Double_t cont_eleptvseta[3];
2027 cont_eleptvseta[0] = trk->Pt();
2028 cont_eleptvseta[1] = trk->Eta();
2031 Double_t cont_eleptvslambdapt[3];
2032 cont_eleptvslambdapt[0] = trk->Pt();
2033 cont_eleptvslambdapt[1] = v0->Pt();
2036 Double_t cont_eleptvsd0[3];
2037 cont_eleptvsd0[0] = trk->Pt();
2038 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
2041 Double_t cont_eleptvsv0dl[3];
2042 cont_eleptvsv0dl[0] = trk->Pt();
2043 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2046 Double_t cont_elelamptvsv0dl[4];
2047 cont_elelamptvsv0dl[0] = elobj->Pt();
2048 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
2049 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
2052 Double_t cont_elelamptvsv0dl_flip[4];
2053 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
2054 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
2055 cont_elelamptvsv0dl_flip[2] = 0.0;
2059 Double_t cont_eleptvsv0dca[3];
2060 cont_eleptvsv0dca[0] = trk->Pt();
2061 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
2064 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
2065 Double_t elobj_px = elobj->Px();
2066 Double_t elobj_py = elobj->Py();
2067 Double_t elobj_pz = elobj->Pz();
2068 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
2069 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
2076 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2103 for(Int_t ih=0;ih<17;ih++){
2104 Double_t cont_eleptvscutvars[3];
2105 cont_eleptvscutvars[0] = elobj->Pt();
2109 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2111 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2113 cont_eleptvscutvars[1] = nSigmaTPCele;
2115 cont_eleptvscutvars[1] = nSigmaTOFele;
2117 cont_eleptvscutvars[1] = trk->Eta();
2119 cont_eleptvscutvars[1] = trk->GetITSNcls();
2121 if(!anti_lambda_flag)
2122 cont_eleptvscutvars[1] = v0->MassLambda();
2124 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2127 lPosV0[0] = v0->DecayVertexV0X();
2128 lPosV0[1] = v0->DecayVertexV0Y();
2129 lPosV0[2] = v0->DecayVertexV0Z();
2130 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2132 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2134 if(!anti_lambda_flag)
2135 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2137 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2139 if(!anti_lambda_flag)
2140 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2142 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2144 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2146 cont_eleptvscutvars[1] = v0->MassK0Short();
2148 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2150 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2152 cont_eleptvscutvars[1] = v0->Eta();
2154 Double_t v0px = elobj->PxProng(1);
2155 Double_t v0py = elobj->PyProng(1);
2156 Double_t v0pz = elobj->PzProng(1);
2157 Double_t epx = elobj->PxProng(0);
2158 Double_t epy = elobj->PyProng(0);
2159 Double_t epz = elobj->PzProng(0);
2160 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2162 cont_eleptvscutvars[1] = -9999.;
2170 Int_t pdgcode = mclc->GetPdgCode();
2171 Double_t cont_mclc[3];
2172 cont_mclc[0] = mclc->Pt();
2173 cont_mclc[1] = mclc->Y();
2175 Double_t cont_mcele[3];
2176 cont_mcele[0] = mcele->Pt();
2177 cont_mcele[1] = mcele->Eta();
2180 Int_t labmotherlc = mclc->GetMother();
2181 Bool_t isbottomfd = kFALSE;
2183 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2184 Int_t pdgmotherlc = motherlc->GetPdgCode();
2185 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){
2189 cont2[1] = mcele->Pt();
2191 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2243 Double_t cont_eleptvslambdaptvslcpt[4];
2244 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2245 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2246 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2247 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2250 Double_t cont_allpt[4];
2251 cont_allpt[0] = mclc->Pt();
2252 cont_allpt[1] = elobj->Pt();
2253 cont_allpt[2] = trk->Pt();
2254 cont_allpt[3] = v0->Pt();
2263 for(Int_t ih=0;ih<17;ih++){
2264 Double_t cont_eleptvscutvars[3];
2265 cont_eleptvscutvars[0] = elobj->Pt();
2269 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2271 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2273 cont_eleptvscutvars[1] = nSigmaTPCele;
2275 cont_eleptvscutvars[1] = nSigmaTOFele;
2277 cont_eleptvscutvars[1] = trk->Eta();
2279 cont_eleptvscutvars[1] = trk->GetITSNcls();
2281 if(!anti_lambda_flag)
2282 cont_eleptvscutvars[1] = v0->MassLambda();
2284 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2287 lPosV0[0] = v0->DecayVertexV0X();
2288 lPosV0[1] = v0->DecayVertexV0Y();
2289 lPosV0[2] = v0->DecayVertexV0Z();
2290 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2292 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2294 if(!anti_lambda_flag)
2295 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2297 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2299 if(!anti_lambda_flag)
2300 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2302 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2304 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2306 cont_eleptvscutvars[1] = v0->MassK0Short();
2308 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2310 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2312 cont_eleptvscutvars[1] = v0->Eta();
2314 Double_t v0px = elobj->PxProng(1);
2315 Double_t v0py = elobj->PyProng(1);
2316 Double_t v0pz = elobj->PzProng(1);
2317 Double_t epx = elobj->PxProng(0);
2318 Double_t epy = elobj->PyProng(0);
2319 Double_t epz = elobj->PzProng(0);
2320 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2322 cont_eleptvscutvars[1] = -9999.;
2329 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2376 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2423 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2453 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2480 for(Int_t ih=0;ih<17;ih++){
2481 Double_t cont_eleptvscutvars[3];
2482 cont_eleptvscutvars[0] = elobj->Pt();
2486 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2488 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2490 cont_eleptvscutvars[1] = nSigmaTPCele;
2492 cont_eleptvscutvars[1] = nSigmaTOFele;
2494 cont_eleptvscutvars[1] = trk->Eta();
2496 cont_eleptvscutvars[1] = trk->GetITSNcls();
2498 if(!anti_lambda_flag)
2499 cont_eleptvscutvars[1] = v0->MassLambda();
2501 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2504 lPosV0[0] = v0->DecayVertexV0X();
2505 lPosV0[1] = v0->DecayVertexV0Y();
2506 lPosV0[2] = v0->DecayVertexV0Z();
2507 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2509 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2511 if(!anti_lambda_flag)
2512 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2514 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2516 if(!anti_lambda_flag)
2517 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2519 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2521 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2523 cont_eleptvscutvars[1] = v0->MassK0Short();
2525 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2527 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2529 cont_eleptvscutvars[1] = v0->Eta();
2531 Double_t v0px = elobj->PxProng(1);
2532 Double_t v0py = elobj->PyProng(1);
2533 Double_t v0pz = elobj->PzProng(1);
2534 Double_t epx = elobj->PxProng(0);
2535 Double_t epy = elobj->PyProng(0);
2536 Double_t epz = elobj->PzProng(0);
2537 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2539 cont_eleptvscutvars[1] = -9999.;
2547 Int_t pdgcode = mclc->GetPdgCode();
2548 Double_t cont_mclb[3];
2549 cont_mclb[0] = mclc->Pt();
2550 cont_mclb[1] = mclc->Y();
2552 Double_t cont_mcele[3];
2553 cont_mcele[0] = mcele->Pt();
2554 cont_mcele[1] = mcele->Eta();
2557 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2578 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2590 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2602 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2614 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2634 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2636 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2640 if(cont_flip[0]<2.3){
2649 if(cont_flip[0]<2.3){
2659 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2680 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2704 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2706 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2710 if(cont_flip[0]<2.3){
2719 if(cont_flip[0]<2.3){
2730 for(Int_t iv=0;iv<13;iv++){
2734 Double_t cont_cor_nd[7];
2735 for(Int_t iv=0;iv<7;iv++){
2736 cont_cor_nd[iv] = -9999.;
2747 if(trk->Charge()>0){
2750 }
else if(trk->Charge()<0){
2759 cont_cor_nd[0] = elobj->Pt();
2761 cont_cor_nd[2] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;;
2762 if(trk->Charge()>0){
2763 if(!anti_lambda_flag) cont_cor_nd[3] = 0;
2764 else cont_cor_nd[3] = 2;
2765 }
else if(trk->Charge()<0){
2766 if(!anti_lambda_flag) cont_cor_nd[3] = 3;
2767 else cont_cor_nd[3] = 1;
2775 Int_t pdgcode = mclc->GetPdgCode();
2776 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2780 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2784 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2788 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2792 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2796 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2800 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2804 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2808 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2846 if(elobj->InvMass(2,pdgdg)<2.3)
2848 cont_cor_nd[0] = trk->Pt();
2850 cont_cor_nd[0] = v0->Pt();
2853 if(elobj->InvMass(2,pdgdg)<2.3)
2855 cont_cor_nd[0] = trk->Pt();
2857 cont_cor_nd[0] = v0->Pt();
2874 for(Int_t i=0;i<94;i++)
2880 Double_t pxe = trke->Px();
2881 Double_t pye = trke->Py();
2882 Double_t pze = trke->Pz();
2883 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
2884 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
2886 Double_t pxv = v0->Px();
2887 Double_t pyv = v0->Py();
2888 Double_t pzv = v0->Pz();
2889 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
2890 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
2892 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
2894 Double_t pxsum = pxe + pxv;
2895 Double_t pysum = pye + pyv;
2896 Double_t pzsum = pze + pzv;
2897 Double_t Esum = Ee + Ev;
2898 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
2900 Double_t uxe = pxe/mome;
2901 Double_t uye = pye/mome;
2902 Double_t uze = pze/mome;
2903 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
2904 Double_t pxv_flip = pxv + lf * uxe;
2905 Double_t pyv_flip = pyv + lf * uye;
2906 Double_t pzv_flip = pzv + lf * uze;
2907 Double_t pxsum_flip = pxe + pxv_flip;
2908 Double_t pysum_flip = pye + pyv_flip;
2909 Double_t pzsum_flip = pze + pzv_flip;
2910 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
2911 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
2914 UInt_t pdgdg[2]={11,3122};
2967 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
2971 Double_t cont_flip[3];
2972 cont_flip[0] = mel_flip;
2973 cont_flip[1] = ptel_flip;
2978 cont2[1] = sqrt(pxe*pxe+pye*pye);
2980 Double_t cont_eleptvseta[3];
2981 cont_eleptvseta[0] = trke->Pt();
2982 cont_eleptvseta[1] = trke->Eta();
2985 Double_t cont_eleptvslambdapt[3];
2986 cont_eleptvslambdapt[0] = trke->Pt();
2987 cont_eleptvslambdapt[1] = v0->Pt();
2990 Double_t cont_eleptvsd0[3];
2991 cont_eleptvsd0[0] = trke->Pt();
2992 cont_eleptvsd0[1] = 0.;
2995 Double_t cont_eleptvsv0dl[3];
2996 cont_eleptvsv0dl[0] = trke->Pt();
2997 cont_eleptvsv0dl[1] = v0info[0];
3000 Double_t cont_eleptvsv0dca[3];
3001 cont_eleptvsv0dca[0] = trke->Pt();
3002 cont_eleptvsv0dca[1] = v0info[1];
3005 Double_t cont_elelamptvsv0dl[4];
3006 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
3007 cont_elelamptvsv0dl[1] = v0info[0];
3008 cont_elelamptvsv0dl[2] = 0.0;
3011 Double_t cont_elelamptvsv0dl_flip[4];
3012 cont_elelamptvsv0dl_flip[0] = ptel_flip;
3013 cont_elelamptvsv0dl_flip[1] = v0info[0];
3014 cont_elelamptvsv0dl_flip[2] = 0.0;
3017 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
3018 xyzR125_ele[0] = (*elevars)[0];
3019 xyzR125_ele[1] = (*elevars)[1];
3020 xyzR125_ele[2] = (*elevars)[2];
3021 xyzR125_pr[0] = (*v0vars)[0];
3022 xyzR125_pr[1] = (*v0vars)[1];
3023 xyzR125_pr[2] = (*v0vars)[2];
3024 xyzR125_pi[0] = (*v0vars)[3];
3025 xyzR125_pi[1] = (*v0vars)[4];
3026 xyzR125_pi[2] = (*v0vars)[5];
3028 Double_t rdhfcutvars[9];
3029 rdhfcutvars[0] = xyzR125_ele[0];
3030 rdhfcutvars[1] = xyzR125_ele[1];
3031 rdhfcutvars[2] = xyzR125_ele[2];
3032 rdhfcutvars[3] = xyzR125_pr[0];
3033 rdhfcutvars[4] = xyzR125_pr[1];
3034 rdhfcutvars[5] = xyzR125_pr[2];
3035 rdhfcutvars[6] = xyzR125_pi[0];
3036 rdhfcutvars[7] = xyzR125_pi[1];
3037 rdhfcutvars[8] = xyzR125_pi[2];
3048 if(((
int)trke->T())*chargepr>0){
3108 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
3110 if(((
int)trke->T())*chargepr>0){
3114 if(cont_flip[0]<2.3){
3123 if(cont_flip[0]<2.3){
3134 for(Int_t iv=0;iv<13;iv++){
3138 Double_t cont_cor_nd[7];
3139 for(Int_t iv=0;iv<7;iv++){
3140 cont_cor_nd[iv] = -9999.;
3147 fCorrelationVariables[4] = v0info[0];
3148 fCorrelationVariables[5] = (*elevars)[5];
3149 fCorrelationVariables[6] = 2;
3151 if(chargepr>0) fCorrelationVariables[7] = 0;
3152 else fCorrelationVariables[7] = 2;
3153 }
else if(trke->T()<0){
3154 if(chargepr>0) fCorrelationVariables[7] = 1;
3155 else fCorrelationVariables[7] = 3;
3157 fCorrelationVariables[8] = (*elevars)[6];
3158 fCorrelationVariables[9] = (*elevars)[7];
3160 fCorrelationVariables[11] = sqrt(pxsum*pxsum+pysum*pysum);
3161 fCorrelationVariables[12] = mel;
3163 cont_cor_nd[0] = sqrt(pxsum*pxsum+pysum*pysum);
3164 cont_cor_nd[1] = TVector2::Phi_mpi_pi(v0->Phi()-trke->Phi());
3165 cont_cor_nd[2] = v0info[0];
3167 if(chargepr>0) cont_cor_nd[3] = 0;
3168 else cont_cor_nd[3] = 2;
3169 }
else if(trke->T()<0){
3170 if(chargepr>0) cont_cor_nd[3] = 3;
3171 else cont_cor_nd[3] = 1;
3173 cont_cor_nd[4] = fCorrelationVariables[8];
3175 if(fabs(fCorrelationVariables[9]-1013)<0.001) cont_cor_nd[5] = 7;
3176 if(fabs(fCorrelationVariables[9]-1016)<0.001) cont_cor_nd[5] = 8;
3177 if(fabs(fCorrelationVariables[9]-1019)<0.001) cont_cor_nd[5] = 9;
3187 cont_cor_nd[0] = trke->Pt();
3189 cont_cor_nd[0] = v0->Pt();
3202 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
3206 TString * fCandidateVariableNames =
new TString[
nVar];
3208 fCandidateVariableNames[ 0]=
"ElePx";
3209 fCandidateVariableNames[ 1]=
"ElePy";
3210 fCandidateVariableNames[ 2]=
"ElePz";
3211 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
3212 fCandidateVariableNames[ 4]=
"ITSNcls";
3213 fCandidateVariableNames[ 5]=
"TPCNcls";
3214 fCandidateVariableNames[ 6]=
"TPCNclsPID";
3215 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
3216 fCandidateVariableNames[ 8]=
"d0R";
3217 fCandidateVariableNames[ 9]=
"d0Z";
3218 fCandidateVariableNames[10]=
"ITSClusterMap";
3219 fCandidateVariableNames[11]=
"nSigmaTPCele";
3220 fCandidateVariableNames[12]=
"nSigmaTOFele";
3221 fCandidateVariableNames[13]=
"nSigmaTPCpi";
3222 fCandidateVariableNames[14]=
"nSigmaTPCka";
3223 fCandidateVariableNames[15]=
"nSigmaTPCpr";
3224 fCandidateVariableNames[16]=
"EvNumber";
3225 fCandidateVariableNames[17]=
"EleCharge";
3226 fCandidateVariableNames[18]=
"ElePdgCode";
3227 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3228 fCandidateVariableNames[20]=
"mcelepx";
3229 fCandidateVariableNames[21]=
"mcelepy";
3230 fCandidateVariableNames[22]=
"mcelepz";
3231 fCandidateVariableNames[23]=
"Centrality";
3232 fCandidateVariableNames[24]=
"PrimVertZ";
3233 fCandidateVariableNames[25]=
"RunNumber";
3235 for (Int_t ivar=0; ivar<
nVar; ivar++) {
3253 Double_t d0z0[2],covd0z0[3];
3254 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3258 pv[0] =
fVtx1->GetX();
3259 pv[1] =
fVtx1->GetY();
3260 pv[2] =
fVtx1->GetZ();
3261 Double_t xyzR125[3]={9999.,9999.,9999.};
3263 TVector *varvec =
new TVector(8);
3264 (*varvec)[0] = xyzR125[0];
3265 (*varvec)[1] = xyzR125[1];
3266 (*varvec)[2] = xyzR125[2];
3269 (*varvec)[5] = d0z0[0];
3270 (*varvec)[6] = convtype;
3271 (*varvec)[7] = mcetype;
3273 fElectronTracks->AddLast(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3277 Int_t pdgEle = -9999;
3278 Int_t pdgEleMother = -9999;
3279 Float_t mcelepx = -9999;
3280 Float_t mcelepy = -9999;
3281 Float_t mcelepz = -9999;
3284 Int_t labEle = trk->GetLabel();
3285 if(labEle<0)
return;
3286 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3288 pdgEle = mcetrk->GetPdgCode();
3289 mcelepx = mcetrk->Px();
3290 mcelepy = mcetrk->Py();
3291 mcelepz = mcetrk->Pz();
3308 for(Int_t i=0;i<26;i++){
3319 if(trk->GetTPCNclsF()>0)
3324 Int_t itsmap = trk->GetITSClusterMap();
3327 Bool_t spdfirst = (itsmap & bit1) == bit1;
3328 Bool_t spdsecond = (itsmap & bit2) == bit2;
3366 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3370 TString * fCandidateVariableNames =
new TString[
nVar];
3372 fCandidateVariableNames[ 0]=
"V0Px";
3373 fCandidateVariableNames[ 1]=
"V0Py";
3374 fCandidateVariableNames[ 2]=
"V0Pz";
3375 fCandidateVariableNames[ 3]=
"MassLambda";
3376 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3377 fCandidateVariableNames[ 5]=
"ProtonPx";
3378 fCandidateVariableNames[ 6]=
"ProtonPy";
3379 fCandidateVariableNames[ 7]=
"ProtonPz";
3380 fCandidateVariableNames[ 8]=
"PionPx";
3381 fCandidateVariableNames[ 9]=
"PionPy";
3382 fCandidateVariableNames[10]=
"PionPz";
3383 fCandidateVariableNames[11]=
"RfidV0";
3384 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3385 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3386 fCandidateVariableNames[14]=
"DcaV0daughters";
3387 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3388 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3389 fCandidateVariableNames[17]=
"MassK0Short";
3390 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3391 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3392 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3393 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3394 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3395 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3396 fCandidateVariableNames[24]=
"V0PdgCode";
3397 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3398 fCandidateVariableNames[26]=
"mcv0px";
3399 fCandidateVariableNames[27]=
"mcv0py";
3400 fCandidateVariableNames[28]=
"mcv0pz";
3401 fCandidateVariableNames[29]=
"mcv0vertx";
3402 fCandidateVariableNames[30]=
"mcv0verty";
3403 fCandidateVariableNames[31]=
"mcv0vertz";
3404 fCandidateVariableNames[32]=
"EvNumber";
3405 fCandidateVariableNames[33]=
"Centrality";
3406 fCandidateVariableNames[34]=
"PrimVertZ";
3407 fCandidateVariableNames[35]=
"RunNumber";
3409 for (Int_t ivar=0; ivar<
nVar; ivar++) {
3423 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3424 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3425 if(!cptrack)
return;
3426 if(!cntrack)
return;
3445 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3446 Double_t posVtx[3] = {0.,0.,0.};
3447 fVtx1->GetXYZ(posVtx);
3448 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3449 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3450 Double_t cont_arm[3];
3451 cont_arm[0] = v0->Pt();
3452 cont_arm[1] = v0->PtArmV0();
3453 cont_arm[2] = v0->AlphaV0();
3472 Double_t momv0x = v0->MomV0X();
3473 Double_t momv0y = v0->MomV0Y();
3474 Double_t phi_alice = atan2(momv0y,momv0x);
3475 if(phi_alice<0.) phi_alice += 2 * M_PI;
3478 Int_t v0pdgcode = -9999;
3479 Int_t v0motherpdgcode = -9999;
3480 Float_t mcv0px = -9999;
3481 Float_t mcv0py = -9999;
3482 Float_t mcv0pz = -9999;
3483 Float_t mcv0vertx = -9999;
3484 Float_t mcv0verty = -9999;
3485 Float_t mcv0vertz = -9999;
3488 Int_t pdgdgv0[2]={2212,211};
3489 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
3503 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
3504 if(!mcv0trk)
return;
3506 Bool_t hfv0_flag = kFALSE;
3507 Bool_t promptv0_flag = kFALSE;
3508 v0pdgcode = mcv0trk->GetPdgCode();
3509 Int_t labv0mother = mcv0trk->GetMother();
3511 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3513 v0motherpdgcode = motherv0->GetPdgCode();
3514 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
3517 if(abs(v0motherpdgcode)==3322){
3519 }
else if(abs(v0motherpdgcode)==3312){
3521 }
else if(abs(v0motherpdgcode)==3334){
3525 promptv0_flag = kTRUE;
3529 if(!hfv0_flag)
return;
3532 mcv0px = mcv0trk->Px();
3533 mcv0py = mcv0trk->Py();
3534 mcv0pz = mcv0trk->Pz();
3536 Int_t labptrk = cptrack->GetLabel();
3538 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
3540 mcv0vertx = mcv0ptrk->Xv();
3541 mcv0verty = mcv0ptrk->Yv();
3542 mcv0vertz = mcv0ptrk->Zv();
3549 Double_t posVtx[3] = {0.,0.,0.};
3550 fVtx1->GetXYZ(posVtx);
3551 TLorentzVector *lv =
new TLorentzVector();
3552 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3553 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
3554 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
3556 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
3558 fV0dlArray1.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3562 TVector *varvec =
new TVector(15);
3563 (*varvec)[0] = xyzR125pr[0];
3564 (*varvec)[1] = xyzR125pr[1];
3565 (*varvec)[2] = xyzR125pr[2];
3566 (*varvec)[3] = xyzR125pi[0];
3567 (*varvec)[4] = xyzR125pi[1];
3568 (*varvec)[5] = xyzR125pi[2];
3569 (*varvec)[6] = v0->DecayVertexV0X();
3570 (*varvec)[7] = v0->DecayVertexV0Y();
3571 (*varvec)[8] = v0->DecayVertexV0Z();
3572 (*varvec)[9] = v0->MomPosX();
3573 (*varvec)[10] = v0->MomPosY();
3574 (*varvec)[11] = v0->MomPosZ();
3575 (*varvec)[12] = v0->MomNegX();
3576 (*varvec)[13] = v0->MomNegY();
3577 (*varvec)[14] = v0->MomNegZ();
3580 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
3582 fV0dlArray2.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3586 TVector *varvec =
new TVector(15);
3587 (*varvec)[0] = xyzR125pr[0];
3588 (*varvec)[1] = xyzR125pr[1];
3589 (*varvec)[2] = xyzR125pr[2];
3590 (*varvec)[3] = xyzR125pi[0];
3591 (*varvec)[4] = xyzR125pi[1];
3592 (*varvec)[5] = xyzR125pi[2];
3593 (*varvec)[6] = v0->DecayVertexV0X();
3594 (*varvec)[7] = v0->DecayVertexV0Y();
3595 (*varvec)[8] = v0->DecayVertexV0Z();
3596 (*varvec)[9] = v0->MomNegX();
3597 (*varvec)[10] = v0->MomNegY();
3598 (*varvec)[11] = v0->MomNegZ();
3599 (*varvec)[12] = v0->MomPosX();
3600 (*varvec)[13] = v0->MomPosY();
3601 (*varvec)[14] = v0->MomPosZ();
3608 for(Int_t i=0;i<36;i++){
3618 Bool_t isparticle = kTRUE;
3638 lPosV0[0] = v0->DecayVertexV0X();
3639 lPosV0[1] = v0->DecayVertexV0Y();
3640 lPosV0[2] = v0->DecayVertexV0Z();
3641 Double_t decayvertV0 = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
3644 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
3645 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
3675 if(cptrack->GetTPCNclsF()>0)
3676 fCandidateV0Variables[22] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3677 if(cntrack->GetTPCNclsF()>0)
3678 fCandidateV0Variables[23] =(Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3682 if(cntrack->GetTPCNclsF()>0)
3683 fCandidateV0Variables[22] = (Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3684 if(cptrack->GetTPCNclsF()>0)
3685 fCandidateV0Variables[23] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3710 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
3714 TString * fCandidateVariableNames =
new TString[
nVar];
3716 fCandidateVariableNames[ 0]=
"Centrality";
3717 fCandidateVariableNames[ 1]=
"DecayType";
3718 fCandidateVariableNames[ 2]=
"LcPx";
3719 fCandidateVariableNames[ 3]=
"LcPy";
3720 fCandidateVariableNames[ 4]=
"LcPz";
3721 fCandidateVariableNames[ 5]=
"ElePx";
3722 fCandidateVariableNames[ 6]=
"ElePy";
3723 fCandidateVariableNames[ 7]=
"ElePz";
3724 fCandidateVariableNames[ 8]=
"V0Px";
3725 fCandidateVariableNames[ 9]=
"V0Py";
3726 fCandidateVariableNames[10]=
"V0Pz";
3728 for (Int_t ivar=0; ivar<
nVar; ivar++) {
3740 if(!mcepart)
return;
3741 if(!mcv0part)
return;
3743 for(Int_t i=0;i<11;i++){
3759 Double_t epx = mcepart->Px();
3760 Double_t epy = mcepart->Py();
3761 Double_t epz = mcepart->Pz();
3762 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
3763 Double_t v0px = mcv0part->Px();
3764 Double_t v0py = mcv0part->Py();
3765 Double_t v0pz = mcv0part->Pz();
3766 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
3768 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
3771 cont[0] = InvMassEleLambda;
3772 cont[1] = mcpart->Pt();
3775 cont2[0] = InvMassEleLambda;
3776 cont2[1] = mcepart->Pt();
3778 Double_t cont_eleptvseta[3];
3779 cont_eleptvseta[0] = mcepart->Pt();
3780 cont_eleptvseta[1] = mcepart->Eta();
3782 Double_t cont_eleptvslambdapt[3];
3783 cont_eleptvslambdapt[0] = mcepart->Pt();
3784 cont_eleptvslambdapt[1] = mcv0part->Pt();
3786 Double_t cont_eleptvslambdaptvslcpt[4];
3787 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
3788 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
3789 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
3792 contmc[0] = mcpart->Pt();
3793 contmc[1] = mcpart->Y();
3795 Double_t contmcele[3];
3796 contmcele[0] = mcepart->Pt();
3797 contmcele[1] = mcepart->Eta();
3799 Double_t contmc_withele[3];
3800 contmc_withele[0] = mcpart->Pt();
3801 contmc_withele[1] = mcpart->Y();
3802 contmc_withele[2] = mcepart->Y();
3806 esdcuts->GetEtaRange(etamin,etamax);
3813 if(fabs(mcepart->Eta())<0.8){
3820 if(fabs(mcepart->Eta())<etamax){
3824 if(InvMassEleLambda<2.3){
3830 if(fabs(mcpart->Y())<0.7){
3831 if(InvMassEleLambda<2.3){
3836 }
else if(decaytype==1){
3844 if(fabs(mcepart->Eta())<etamax){
3848 if(InvMassEleLambda<2.3){
3854 }
else if(decaytype==2){
3862 if(fabs(mcepart->Eta())<etamax){
3866 if(InvMassEleLambda<2.3){
3872 }
else if(decaytype==3){
3876 }
else if(decaytype==10){
3885 if(fabs(mcepart->Eta())<0.8){
3888 }
else if(decaytype==11){
3904 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
3908 TString * fCandidateVariableNames =
new TString[
nVar];
3910 fCandidateVariableNames[ 0]=
"Centrality";
3911 fCandidateVariableNames[ 1]=
"ElePx";
3912 fCandidateVariableNames[ 2]=
"ElePy";
3913 fCandidateVariableNames[ 3]=
"ElePz";
3914 fCandidateVariableNames[ 4]=
"ElePdgCode";
3915 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
3916 fCandidateVariableNames[ 6]=
"RunNumber";
3917 fCandidateVariableNames[ 7]=
"EvNumber";
3919 for (Int_t ivar=0; ivar<
nVar; ivar++) {
3930 if(!mcepart)
return;
3932 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
3937 Int_t semi_flag = 0;
3938 if(pdgarray_ele[0]>400&&pdgarray_ele[0]<450) semi_flag = 1;
3939 if(pdgarray_ele[0]>4000&&pdgarray_ele[0]<4500) semi_flag = 1;
3940 if(pdgarray_ele[0]>500&&pdgarray_ele[0]<550) semi_flag = 2;
3941 if(pdgarray_ele[0]>5000&&pdgarray_ele[0]<5500) semi_flag = 2;
3944 contmc[0] = mcepart->Pt();
3945 contmc[1] = mcepart->Eta();
3948 if(semi_flag==1 && !ele_from_bottom){
3951 if(semi_flag==1 && ele_from_bottom){
3958 Bool_t hfe_flag = kFALSE;
3959 Int_t labemother = mcepart->GetMother();
3960 Int_t pdgmotherele = -9999;
3962 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
3963 pdgmotherele = motherele->GetPdgCode();
3964 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
3969 if(!hfe_flag)
return;
3973 for(Int_t i=0;i<8;i++){
3997 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
4001 TString * fCandidateVariableNames =
new TString[
nVar];
4003 fCandidateVariableNames[ 0]=
"Centrality";
4004 fCandidateVariableNames[ 1]=
"V0Px";
4005 fCandidateVariableNames[ 2]=
"V0Py";
4006 fCandidateVariableNames[ 3]=
"V0Pz";
4007 fCandidateVariableNames[ 4]=
"V0PdgCode";
4008 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
4009 fCandidateVariableNames[ 6]=
"RunNumber";
4010 fCandidateVariableNames[ 7]=
"EvNumber";
4012 for (Int_t ivar=0; ivar<
nVar; ivar++) {
4023 if(!mcv0part)
return;
4025 for(Int_t i=0;i<8;i++){
4029 Bool_t hfv0_flag = kFALSE;
4030 Int_t labv0mother = mcv0part->GetMother();
4031 Int_t pdgmotherv0 = -9999;
4033 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
4035 pdgmotherv0 = motherv0->GetPdgCode();
4036 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
4041 if(!hfv0_flag)
return;
4044 contmc[0] = mcv0part->Pt();
4045 contmc[1] = mcv0part->Eta();
4069 fCEvents =
new TH1F(
"fCEvents",
"conter",18,-0.5,17.5);
4071 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
4072 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
4073 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
4074 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
4075 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
4076 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
4077 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
4078 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
4079 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
4080 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
4081 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
4082 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
4083 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
4084 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
4086 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
4087 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
4090 fCEvents->GetYaxis()->SetTitle(
"counts");
4092 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
4094 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
4095 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
4096 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
4097 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
4098 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
4099 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
4101 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
4102 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
4103 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
4104 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
4105 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
4107 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
4126 Int_t bins_base[3]= {23,100 ,10};
4127 Double_t xmin_base[3]={1.1,0 ,0.00};
4128 Double_t xmax_base[3]={5.7,20. ,100};
4129 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
4205 Int_t bins_base_elept[3]= {10,100 ,10};
4206 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
4207 Double_t xmax_base_elept[3]={3.1,10. ,100};
4266 fHistoElePtRS =
new TH2F(
"fHistoElePtRS",
"",100,0,10,10,0,100);
4268 fHistoElePtWS =
new TH2F(
"fHistoElePtWS",
"",100,0.,10.,10,0,100);
4274 fHistoElePtMCS =
new TH2F(
"fHistoElePtMCS",
"",100,0,10,10,0,100);
4279 Int_t bins_eleptvseta[3]= {10,10 ,10};
4280 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4281 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4283 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4285 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4287 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4289 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4291 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4293 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4296 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4297 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4298 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4300 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4302 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4304 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4306 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4308 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4310 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4312 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4315 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4316 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4317 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4318 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4320 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4323 Int_t bins_allpt[4]= {10,10,20,20};
4324 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4325 Double_t xmax_allpt[4]={20.,20.,10.,10};
4336 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4337 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4338 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4340 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4342 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4344 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4346 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4348 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4350 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4416 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4422 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4424 fHistoBachPtMCS =
new TH1F(
"fHistoBachPtMCS",
"Bachelor p_{T}",100,0.0,5.0);
4428 fHistod0Bach =
new TH1F(
"fHistod0Bach",
"Bachelor d_{0}",100,-0.5,0.5);
4430 fHistoLambdaMassvsPt=
new TH2F(
"fHistoLambdaMassvsPt",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4432 fHistoLambdaMassvsPtMCS=
new TH2F(
"fHistoLambdaMassvsPtMCS",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4442 fHistoK0sMassvsPt=
new TH2F(
"fHistoK0sMassvsPt",
"K0s mass",100,0.497-0.05,0.497+0.05,20,0.,10.);
4444 fHistoLambdaPtvsDl=
new TH2F(
"fHistoLambdaPtvsDl",
"Lambda pt vs dl",20,0.,10.,20,0.,40.);
4457 Int_t bins_armenteros[3]= {20 ,100 ,200};
4458 Double_t xmin_armenteros[3]={0.,0. ,-1.};
4459 Double_t xmax_armenteros[3]={20.,1. ,1.};
4460 fHistoLambdaArmenteros =
new THnSparseF(
"fHistoLambdaArmenteros",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4462 fHistoLambdaArmenterosMCS =
new THnSparseF(
"fHistoLambdaArmenterosMCS",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4484 for(Int_t i=0;i<8;i++){
4493 Int_t bins_lcmcgen_withelea[3]= {20 ,100 ,100};
4494 Double_t xmin_lcmcgen_withelea[3]={0.,-5.,-5.};
4495 Double_t xmax_lcmcgen_withelea[3]={20.,5.,5.};
4496 fHistoLcMCGenWithEleA =
new THnSparseF(
"fHistoLcMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4501 fHistoLbMCGenWithEleA =
new THnSparseF(
"fHistoLbMCGenWithEleA",
"",3,bins_lcmcgen_withelea,xmin_lcmcgen_withelea,xmax_lcmcgen_withelea);
4506 Int_t bins_lbmcgen_withlc[3]= {50 ,100 ,100};
4507 Double_t xmin_lbmcgen_withlc[3]={0.,-5.,-5.};
4508 Double_t xmax_lbmcgen_withlc[3]={50.,5.,5.};
4509 fHistoLbMCGenWithLambdac =
new THnSparseF(
"fHistoLbMCGenWithLambdac",
"",3,bins_lbmcgen_withlc,xmin_lbmcgen_withlc,xmax_lbmcgen_withlc);
4512 Int_t bins_lcmcgen[3]= {100 ,20 ,10};
4513 Double_t xmin_lcmcgen[3]={0.,-1.0 ,0.0};
4514 Double_t xmax_lcmcgen[3]={20.,1.0 ,100};
4515 Int_t bins_lbmcgen[3]= {100 ,20 ,10};
4516 Double_t xmin_lbmcgen[3]={0.,-1.0 ,0.0};
4517 Double_t xmax_lbmcgen[3]={50.,1.0 ,100};
4518 fHistoLcMCGen =
new THnSparseF(
"fHistoLcMCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4520 fHistoLcMCGen1 =
new THnSparseF(
"fHistoLcMCGen1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4522 fHistoLcMCGen2 =
new THnSparseF(
"fHistoLcMCGen2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4524 fHistoLbMCGen =
new THnSparseF(
"fHistoLbMCGen",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4526 fHistoLbMCGen1 =
new THnSparseF(
"fHistoLbMCGen1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4528 fHistoLbMCGen2 =
new THnSparseF(
"fHistoLbMCGen2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4530 fHistoFeedDownXic0MCGen =
new THnSparseF(
"fHistoFeedDownXic0MCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4549 fHistoLcMCS =
new THnSparseF(
"fHistoLcMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4551 fHistoLcMCS1 =
new THnSparseF(
"fHistoLcMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4553 fHistoLcMCS2 =
new THnSparseF(
"fHistoLcMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4555 fHistoLbMCS =
new THnSparseF(
"fHistoLbMCS",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4557 fHistoLbMCS1 =
new THnSparseF(
"fHistoLbMCS1",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4559 fHistoLbMCS2 =
new THnSparseF(
"fHistoLbMCS2",
"",3,bins_lbmcgen,xmin_lbmcgen,xmax_lbmcgen);
4561 fHistoFeedDownXic0MCS =
new THnSparseF(
"fHistoFeedDownXic0MCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4563 fHistoFeedDownXic0MCS1 =
new THnSparseF(
"fHistoFeedDownXic0MCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4565 fHistoFeedDownXic0MCS2 =
new THnSparseF(
"fHistoFeedDownXic0MCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4573 fHistoFeedDownSigmaMCS =
new THnSparseF(
"fHistoFeedDownSigmaMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4575 fHistoFeedDownSigmaMCS1 =
new THnSparseF(
"fHistoFeedDownSigmaMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4577 fHistoFeedDownSigmaMCS2 =
new THnSparseF(
"fHistoFeedDownSigmaMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4580 Int_t bins_elemcgen[3]= {100 ,20 ,10};
4581 Double_t xmin_elemcgen[3]={0.,-1.0 ,0.0};
4582 Double_t xmax_elemcgen[3]={10.,1.0 ,100};
4583 fHistoElectronMCGen =
new THnSparseF(
"fHistoElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4587 fHistoCharmElectronMCGen =
new THnSparseF(
"fHistoCharmElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4589 fHistoLcElectronMCGen =
new THnSparseF(
"fHistoLcElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4591 fHistoLcElectronMCGen1 =
new THnSparseF(
"fHistoLcElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4593 fHistoLcElectronMCGen2 =
new THnSparseF(
"fHistoLcElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4595 fHistoLbElectronMCGen =
new THnSparseF(
"fHistoLbElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4597 fHistoLbElectronMCGen1 =
new THnSparseF(
"fHistoLbElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4599 fHistoLbElectronMCGen2 =
new THnSparseF(
"fHistoLbElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4619 fHistoLcElectronMCS =
new THnSparseF(
"fHistoLcElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4621 fHistoLcElectronMCS1 =
new THnSparseF(
"fHistoLcElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4623 fHistoLcElectronMCS2 =
new THnSparseF(
"fHistoLcElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4625 fHistoLbElectronMCS =
new THnSparseF(
"fHistoLbElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4627 fHistoLbElectronMCS1 =
new THnSparseF(
"fHistoLbElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4629 fHistoLbElectronMCS2 =
new THnSparseF(
"fHistoLbElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4650 Int_t bins_lambdamcgen[3]= {50 ,20 ,10};
4651 Double_t xmin_lambdamcgen[3]={0.,-1.0 ,0.0};
4652 Double_t xmax_lambdamcgen[3]={10.,1.0 ,100};
4653 fHistoLambdaMCGen =
new THnSparseF(
"fHistoLambdaMCGen",
"",3,bins_lambdamcgen,xmin_lambdamcgen,xmax_lambdamcgen);
4656 Int_t bins_eleptvsv0dl[3]= {100 ,20 ,10};
4657 Double_t xmin_eleptvsv0dl[3]={0.,0. ,0.0};
4658 Double_t xmax_eleptvsv0dl[3]={10.,40. ,100};
4659 fHistoElePtvsV0dlRS =
new THnSparseF(
"fHistoElePtvsV0dlRS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4661 fHistoElePtvsV0dlRS1 =
new THnSparseF(
"fHistoElePtvsV0dlRS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4663 fHistoElePtvsV0dlRS2 =
new THnSparseF(
"fHistoElePtvsV0dlRS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4665 fHistoElePtvsV0dlRSSide =
new THnSparseF(
"fHistoElePtvsV0dlRSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4667 fHistoElePtvsV0dlRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4669 fHistoElePtvsV0dlRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4671 fHistoElePtvsV0dlRSMix =
new THnSparseF(
"fHistoElePtvsV0dlRSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4673 fHistoElePtvsV0dlRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4675 fHistoElePtvsV0dlRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4677 fHistoElePtvsV0dlWS =
new THnSparseF(
"fHistoElePtvsV0dlWS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4679 fHistoElePtvsV0dlWS1 =
new THnSparseF(
"fHistoElePtvsV0dlWS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4681 fHistoElePtvsV0dlWS2 =
new THnSparseF(
"fHistoElePtvsV0dlWS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4683 fHistoElePtvsV0dlWSSide =
new THnSparseF(
"fHistoElePtvsV0dlWSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4685 fHistoElePtvsV0dlWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4687 fHistoElePtvsV0dlWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4689 fHistoElePtvsV0dlWSMix =
new THnSparseF(
"fHistoElePtvsV0dlWSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4691 fHistoElePtvsV0dlWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4693 fHistoElePtvsV0dlWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4695 fHistoElePtvsV0dlMCS =
new THnSparseF(
"fHistoElePtvsV0dlMCS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4697 fHistoElePtvsV0dlMCS1 =
new THnSparseF(
"fHistoElePtvsV0dlMCS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4699 fHistoElePtvsV0dlMCS2 =
new THnSparseF(
"fHistoElePtvsV0dlMCS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4714 Int_t bins_eleptvsv0dca[3]= {100 ,20 ,10};
4715 Double_t xmin_eleptvsv0dca[3]={0.,0. ,0.0};
4716 Double_t xmax_eleptvsv0dca[3]={10.,1. ,100};
4717 fHistoElePtvsV0dcaRS =
new THnSparseF(
"fHistoElePtvsV0dcaRS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4719 fHistoElePtvsV0dcaRS1 =
new THnSparseF(
"fHistoElePtvsV0dcaRS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4721 fHistoElePtvsV0dcaRS2 =
new THnSparseF(
"fHistoElePtvsV0dcaRS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4723 fHistoElePtvsV0dcaRSSide =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4725 fHistoElePtvsV0dcaRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4727 fHistoElePtvsV0dcaRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4729 fHistoElePtvsV0dcaRSMix =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4731 fHistoElePtvsV0dcaRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4733 fHistoElePtvsV0dcaRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4735 fHistoElePtvsV0dcaWS =
new THnSparseF(
"fHistoElePtvsV0dcaWS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4737 fHistoElePtvsV0dcaWS1 =
new THnSparseF(
"fHistoElePtvsV0dcaWS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4739 fHistoElePtvsV0dcaWS2 =
new THnSparseF(
"fHistoElePtvsV0dcaWS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4741 fHistoElePtvsV0dcaWSSide =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4743 fHistoElePtvsV0dcaWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4745 fHistoElePtvsV0dcaWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4747 fHistoElePtvsV0dcaWSMix =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4749 fHistoElePtvsV0dcaWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4751 fHistoElePtvsV0dcaWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4753 fHistoElePtvsV0dcaMCS =
new THnSparseF(
"fHistoElePtvsV0dcaMCS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4755 fHistoElePtvsV0dcaMCS1 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4757 fHistoElePtvsV0dcaMCS2 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4772 Int_t bins_elelamptvsv0dl[4]= {100 ,20 ,100, 10};
4773 Double_t xmin_elelamptvsv0dl[4]={0.,0. ,-0.5, 0.0};
4774 Double_t xmax_elelamptvsv0dl[4]={20.,40.,0.5 ,100};
4775 fHistoEleLambdaPtvsV0dlRS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4777 fHistoEleLambdaPtvsV0dlRS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4779 fHistoEleLambdaPtvsV0dlRS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4781 fHistoEleLambdaPtvsV0dlRSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4783 fHistoEleLambdaPtvsV0dlRSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4785 fHistoEleLambdaPtvsV0dlRSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4787 fHistoEleLambdaPtvsV0dlRSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4789 fHistoEleLambdaPtvsV0dlRSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4791 fHistoEleLambdaPtvsV0dlRSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4793 fHistoEleLambdaPtvsV0dlWS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4795 fHistoEleLambdaPtvsV0dlWS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4797 fHistoEleLambdaPtvsV0dlWS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4799 fHistoEleLambdaPtvsV0dlWSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4801 fHistoEleLambdaPtvsV0dlWSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4803 fHistoEleLambdaPtvsV0dlWSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4805 fHistoEleLambdaPtvsV0dlWSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4807 fHistoEleLambdaPtvsV0dlWSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4809 fHistoEleLambdaPtvsV0dlWSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4811 fHistoEleLambdaPtvsV0dlMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4813 fHistoEleLambdaPtvsV0dlMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4815 fHistoEleLambdaPtvsV0dlMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4891 fHistoEleLambdaPtvsV0dlXibMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4893 fHistoEleLambdaPtvsV0dlXibMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4895 fHistoEleLambdaPtvsV0dlXibMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4898 fHistoEleLambdaPtvsV0dlRSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4900 fHistoEleLambdaPtvsV0dlRS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4902 fHistoEleLambdaPtvsV0dlRS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4916 fHistoEleLambdaPtvsV0dlWSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4918 fHistoEleLambdaPtvsV0dlWS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4920 fHistoEleLambdaPtvsV0dlWS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);