37 #include <TParticle.h>
38 #include <TParticlePDG.h>
42 #include <THnSparse.h>
43 #include <TLorentzVector.h>
46 #include <TDatabasePDG.h>
47 #include <AliAnalysisDataSlot.h>
48 #include <AliAnalysisDataContainer.h>
50 #include "AliMCEvent.h"
51 #include "AliAnalysisManager.h"
52 #include "AliAODMCHeader.h"
53 #include "AliAODHandler.h"
55 #include "AliExternalTrackParam.h"
56 #include "AliAODVertex.h"
57 #include "AliESDVertex.h"
58 #include "AliAODRecoDecay.h"
61 #include "AliESDtrack.h"
62 #include "AliAODTrack.h"
64 #include "AliAODcascade.h"
65 #include "AliAODMCParticle.h"
66 #include "AliAnalysisTaskSE.h"
68 #include "AliPIDResponse.h"
69 #include "AliPIDCombined.h"
70 #include "AliTOFPIDResponse.h"
72 #include "AliInputEventHandler.h"
73 #include "AliESDtrackCuts.h"
74 #include "AliNeutralTrackParam.h"
75 #include "AliKFParticle.h"
76 #include "AliKFVertex.h"
77 #include "AliExternalTrackParam.h"
78 #include "AliESDtrack.h"
79 #include "AliCentrality.h"
80 #include "AliVertexerTracks.h"
81 #include "AliEventPoolManager.h"
102 fIsEventSelected(kFALSE),
103 fWriteVariableTree(kFALSE),
104 fWriteEachVariableTree(kFALSE),
105 fWriteMCVariableTree(kFALSE),
107 fEleVariablesTree(0),
110 fMCEleVariablesTree(0),
111 fMCV0VariablesTree(0),
112 fMCGenPairVariablesTree(0),
113 fCorrelationVariablesTree(0),
114 fReconstructPrimVert(kFALSE),
120 fCandidateVariables(),
121 fCandidateEleVariables(),
122 fCandidateV0Variables(),
123 fCandidateMCVariables(),
124 fCandidateMCEleVariables(),
125 fCandidateMCV0Variables(),
126 fCandidateMCGenPairVariables(),
127 fCorrelationVariables(),
135 fUseCentralityV0M(kFALSE),
138 fMCDoPairAnalysis(kFALSE),
139 fHistoEleLambdaMass(0),
140 fHistoEleLambdaMassRS(0),
141 fHistoEleLambdaMassRS1(0),
142 fHistoEleLambdaMassRS2(0),
143 fHistoEleLambdaMassWS(0),
144 fHistoEleLambdaMassWS1(0),
145 fHistoEleLambdaMassWS2(0),
146 fHistoEleLambdaMassRSMix(0),
147 fHistoEleLambdaMassRSMix1(0),
148 fHistoEleLambdaMassRSMix2(0),
149 fHistoEleLambdaMassWSMix(0),
150 fHistoEleLambdaMassWSMix1(0),
151 fHistoEleLambdaMassWSMix2(0),
152 fHistoEleLambdaMassRSSide(0),
153 fHistoEleLambdaMassRSSide1(0),
154 fHistoEleLambdaMassRSSide2(0),
155 fHistoEleLambdaMassWSSide(0),
156 fHistoEleLambdaMassWSSide1(0),
157 fHistoEleLambdaMassWSSide2(0),
158 fHistoEleLambdaMassRSAway(0),
159 fHistoEleLambdaMassRS1Away(0),
160 fHistoEleLambdaMassRS2Away(0),
161 fHistoEleLambdaMassWSAway(0),
162 fHistoEleLambdaMassWS1Away(0),
163 fHistoEleLambdaMassWS2Away(0),
164 fHistoEleLambdaMassRSMixAway(0),
165 fHistoEleLambdaMassRSMix1Away(0),
166 fHistoEleLambdaMassRSMix2Away(0),
167 fHistoEleLambdaMassWSMixAway(0),
168 fHistoEleLambdaMassWSMix1Away(0),
169 fHistoEleLambdaMassWSMix2Away(0),
170 fHistoEleLambdaMassRSSideAway(0),
171 fHistoEleLambdaMassRSSide1Away(0),
172 fHistoEleLambdaMassRSSide2Away(0),
173 fHistoEleLambdaMassWSSideAway(0),
174 fHistoEleLambdaMassWSSide1Away(0),
175 fHistoEleLambdaMassWSSide2Away(0),
176 fHistoEleLambdaMassvsElePtRS(0),
177 fHistoEleLambdaMassvsElePtWS(0),
178 fHistoEleLambdaMassvsElePtRSMix(0),
179 fHistoEleLambdaMassvsElePtWSMix(0),
180 fHistoEleLambdaMassvsElePtRSSide(0),
181 fHistoEleLambdaMassvsElePtWSSide(0),
182 fHistoEleLambdaMassvsElePtRS1(0),
183 fHistoEleLambdaMassvsElePtWS1(0),
184 fHistoEleLambdaMassvsElePtRSMix1(0),
185 fHistoEleLambdaMassvsElePtWSMix1(0),
186 fHistoEleLambdaMassvsElePtRSSide1(0),
187 fHistoEleLambdaMassvsElePtWSSide1(0),
188 fHistoEleLambdaMassvsElePtRS2(0),
189 fHistoEleLambdaMassvsElePtWS2(0),
190 fHistoEleLambdaMassvsElePtRSMix2(0),
191 fHistoEleLambdaMassvsElePtWSMix2(0),
192 fHistoEleLambdaMassvsElePtRSSide2(0),
193 fHistoEleLambdaMassvsElePtWSSide2(0),
198 fHistoEleLambdaMassMCS(0),
199 fHistoEleLambdaMassLambdabMCS(0),
200 fHistoEleLambdaMassXibMCS(0),
201 fHistoEleLambdaMassMCGen(0),
202 fHistoEleLambdaMassvsElePtMCS(0),
203 fHistoEleLambdaMassvsElePtMCGen(0),
204 fHistoEleLambdaMassvsElePtMCS1(0),
205 fHistoEleLambdaMassvsElePtMCGen1(0),
206 fHistoEleLambdaMassvsElePtMCS2(0),
207 fHistoEleLambdaMassvsElePtMCGen2(0),
210 fHistoElePtvsEtaRS(0),
211 fHistoElePtvsEtaWS(0),
212 fHistoElePtvsEtaRSMix(0),
213 fHistoElePtvsEtaWSMix(0),
214 fHistoElePtvsEtaMCS(0),
215 fHistoElePtvsEtaMCGen(0),
216 fHistoElePtvsLambdaPtRS(0),
217 fHistoElePtvsLambdaPtWS(0),
218 fHistoElePtvsLambdaPtRSMix(0),
219 fHistoElePtvsLambdaPtWSMix(0),
220 fHistoElePtvsLambdaPtMCS(0),
221 fHistoElePtvsLambdaPtvsLcPtMCS(0),
222 fHistoElePtvsLambdaPtMCGen(0),
223 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
224 fHistoElePtvsLambdaPtMCLcGen(0),
225 fHistoElePtvsd0RS(0),
226 fHistoElePtvsd0WS(0),
227 fHistoElePtvsd0RSMix(0),
228 fHistoElePtvsd0WSMix(0),
229 fHistoElePtvsd0MCS(0),
230 fHistoElePtvsd0PromptMCS(0),
231 fHistoElePtvsd0BFeeddownMCS(0),
232 fHistoEleLambdaMassFeeddownXic0MCS(0),
233 fHistoEleLambdaMassFeeddownXic0MCGen(0),
234 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
235 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
236 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
237 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
238 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
239 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
240 fHistoElePtFeeddownXic0MCS(0),
241 fHistoElePtFeeddownXic0MCGen(0),
242 fHistoElePtvsEtaFeeddownXic0MCS(0),
243 fHistoElePtvsEtaFeeddownXic0MCGen(0),
244 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
245 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
246 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
247 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
248 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
249 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
250 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
251 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
252 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
253 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
254 fHistoElePtFeeddownXicPlusMCS(0),
255 fHistoElePtFeeddownXicPlusMCGen(0),
256 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
257 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
258 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
259 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
262 fHistoBachPtMCGen(0),
264 fHistoLambdaMassvsPt(0),
265 fHistoLambdaMassvsPtMCS(0),
266 fHistoLambdaMassvsPtMCGen(0),
267 fHistoLambdaPtvsEta(0),
268 fHistoLambdaPtvsEtaMCS(0),
269 fHistoLambdaPtvsEtaMCGen(0),
270 fHistoLambdaArmenteros(0),
271 fHistoLambdaArmenterosMCS(0),
272 fHistoLambdaPtvsDl(0),
273 fHistoLambdaPtvsDlSide(0),
274 fHistoLambdaPtvsDlMCS(0),
275 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
276 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
277 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
278 fHistoK0sMassvsPt(0),
279 fHistoElectronTPCPID(0),
280 fHistoElectronTOFPID(0),
281 fHistoElectronTPCSelPID(0),
282 fHistoElectronTOFSelPID(0),
283 fHistoElectronTPCPIDSelTOF(0),
284 fHistoElectronTPCPIDSelTOFSmallEta(0),
285 fHistoElectronTPCPIDSelTOFLargeEta(0),
286 fHistoMassConversionsMin(0),
287 fHistoMassConversionsSameSignMin(0),
288 fHistoElectronQovPtvsPhi(0),
289 fHistoLambdaQovPtvsPhi(0),
296 fHistoFeedDownXic0MCGen(0),
297 fHistoFeedDownXic0MCGen1(0),
298 fHistoFeedDownXic0MCGen2(0),
299 fHistoFeedDownXic0MCS(0),
300 fHistoFeedDownXic0MCS1(0),
301 fHistoFeedDownXic0MCS2(0),
302 fHistoFeedDownXicPlusMCGen(0),
303 fHistoFeedDownXicPlusMCGen1(0),
304 fHistoFeedDownXicPlusMCGen2(0),
305 fHistoFeedDownXicPlusMCS(0),
306 fHistoFeedDownXicPlusMCS1(0),
307 fHistoFeedDownXicPlusMCS2(0),
308 fHistoFeedDownSigmaMCGen(0),
309 fHistoFeedDownSigmaMCGen1(0),
310 fHistoFeedDownSigmaMCGen2(0),
311 fHistoFeedDownSigmaMCS(0),
312 fHistoFeedDownSigmaMCS1(0),
313 fHistoFeedDownSigmaMCS2(0),
314 fHistoLcElectronMCGen(0),
315 fHistoLcElectronMCGen1(0),
316 fHistoLcElectronMCGen2(0),
317 fHistoLcElectronMCS(0),
318 fHistoLcElectronMCS1(0),
319 fHistoLcElectronMCS2(0),
320 fHistoElectronFeedDownXic0MCGen(0),
321 fHistoElectronFeedDownXic0MCGen1(0),
322 fHistoElectronFeedDownXic0MCGen2(0),
323 fHistoElectronFeedDownXic0MCS(0),
324 fHistoElectronFeedDownXic0MCS1(0),
325 fHistoElectronFeedDownXic0MCS2(0),
326 fHistoElectronFeedDownXicPlusMCGen(0),
327 fHistoElectronFeedDownXicPlusMCGen1(0),
328 fHistoElectronFeedDownXicPlusMCGen2(0),
329 fHistoElectronFeedDownXicPlusMCS(0),
330 fHistoElectronFeedDownXicPlusMCS1(0),
331 fHistoElectronFeedDownXicPlusMCS2(0),
332 fHistoElectronMCGen(0),
333 fHistoBottomElectronMCGen(0),
334 fHistoCharmElectronMCGen(0),
335 fHistoLambdaMCGen(0),
336 fHistoElePtvsV0dlRS(0),
337 fHistoElePtvsV0dlRS1(0),
338 fHistoElePtvsV0dlRS2(0),
339 fHistoElePtvsV0dlRSSide(0),
340 fHistoElePtvsV0dlRSSide1(0),
341 fHistoElePtvsV0dlRSSide2(0),
342 fHistoElePtvsV0dlRSMix(0),
343 fHistoElePtvsV0dlRSMix1(0),
344 fHistoElePtvsV0dlRSMix2(0),
345 fHistoElePtvsV0dlWS(0),
346 fHistoElePtvsV0dlWS1(0),
347 fHistoElePtvsV0dlWS2(0),
348 fHistoElePtvsV0dlWSSide(0),
349 fHistoElePtvsV0dlWSSide1(0),
350 fHistoElePtvsV0dlWSSide2(0),
351 fHistoElePtvsV0dlWSMix(0),
352 fHistoElePtvsV0dlWSMix1(0),
353 fHistoElePtvsV0dlWSMix2(0),
354 fHistoElePtvsV0dlMCS(0),
355 fHistoElePtvsV0dlMCS1(0),
356 fHistoElePtvsV0dlMCS2(0),
357 fHistoElePtvsV0dlFeedDownXic0MCS(0),
358 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
359 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
360 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
361 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
362 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
363 fHistoElePtvsV0dcaRS(0),
364 fHistoElePtvsV0dcaRS1(0),
365 fHistoElePtvsV0dcaRS2(0),
366 fHistoElePtvsV0dcaRSSide(0),
367 fHistoElePtvsV0dcaRSSide1(0),
368 fHistoElePtvsV0dcaRSSide2(0),
369 fHistoElePtvsV0dcaRSMix(0),
370 fHistoElePtvsV0dcaRSMix1(0),
371 fHistoElePtvsV0dcaRSMix2(0),
372 fHistoElePtvsV0dcaWS(0),
373 fHistoElePtvsV0dcaWS1(0),
374 fHistoElePtvsV0dcaWS2(0),
375 fHistoElePtvsV0dcaWSSide(0),
376 fHistoElePtvsV0dcaWSSide1(0),
377 fHistoElePtvsV0dcaWSSide2(0),
378 fHistoElePtvsV0dcaWSMix(0),
379 fHistoElePtvsV0dcaWSMix1(0),
380 fHistoElePtvsV0dcaWSMix2(0),
381 fHistoElePtvsV0dcaMCS(0),
382 fHistoElePtvsV0dcaMCS1(0),
383 fHistoElePtvsV0dcaMCS2(0),
384 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
385 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
386 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
387 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
388 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
389 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
390 fHistoEleLambdaPtvsV0dlRS(0),
391 fHistoEleLambdaPtvsV0dlRS1(0),
392 fHistoEleLambdaPtvsV0dlRS2(0),
393 fHistoEleLambdaPtvsV0dlRSSide(0),
394 fHistoEleLambdaPtvsV0dlRSSide1(0),
395 fHistoEleLambdaPtvsV0dlRSSide2(0),
396 fHistoEleLambdaPtvsV0dlRSMix(0),
397 fHistoEleLambdaPtvsV0dlRSMix1(0),
398 fHistoEleLambdaPtvsV0dlRSMix2(0),
399 fHistoEleLambdaPtvsV0dlWS(0),
400 fHistoEleLambdaPtvsV0dlWS1(0),
401 fHistoEleLambdaPtvsV0dlWS2(0),
402 fHistoEleLambdaPtvsV0dlWSSide(0),
403 fHistoEleLambdaPtvsV0dlWSSide1(0),
404 fHistoEleLambdaPtvsV0dlWSSide2(0),
405 fHistoEleLambdaPtvsV0dlWSMix(0),
406 fHistoEleLambdaPtvsV0dlWSMix1(0),
407 fHistoEleLambdaPtvsV0dlWSMix2(0),
408 fHistoEleLambdaPtvsV0dlMCS(0),
409 fHistoEleLambdaPtvsV0dlMCS1(0),
410 fHistoEleLambdaPtvsV0dlMCS2(0),
411 fHistoEleLambdaPtvsV0dlPromptMCS(0),
412 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
413 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
414 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
415 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
416 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
417 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
418 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
419 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
420 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
421 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
422 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
423 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
424 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
425 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
426 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
427 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
428 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
429 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
430 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
431 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
432 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
433 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
434 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
435 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
436 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
437 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
438 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
439 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
440 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
441 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
442 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
443 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
444 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
445 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
446 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
447 fHistoEleLambdaPtvsV0dlXibMCS(0),
448 fHistoEleLambdaPtvsV0dlXibMCS1(0),
449 fHistoEleLambdaPtvsV0dlXibMCS2(0),
450 fHistoEleLambdaPtvsV0dlRSAway(0),
451 fHistoEleLambdaPtvsV0dlRS1Away(0),
452 fHistoEleLambdaPtvsV0dlRS2Away(0),
453 fHistoEleLambdaPtvsV0dlRSSideAway(0),
454 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
455 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
456 fHistoEleLambdaPtvsV0dlRSMixAway(0),
457 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
458 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
459 fHistoEleLambdaPtvsV0dlWSAway(0),
460 fHistoEleLambdaPtvsV0dlWS1Away(0),
461 fHistoEleLambdaPtvsV0dlWS2Away(0),
462 fHistoEleLambdaPtvsV0dlWSSideAway(0),
463 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
464 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
465 fHistoEleLambdaPtvsV0dlWSMixAway(0),
466 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
467 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
468 fHistoResponseElePt(0),
469 fHistoResponseElePt1(0),
470 fHistoResponseElePt2(0),
471 fHistoResponseEleLambdaPt(0),
472 fHistoResponseEleLambdaPt1(0),
473 fHistoResponseEleLambdaPt2(0),
474 fHistoResponseEleLambdaPtFeeddownXic0(0),
475 fHistoResponseEleLambdaPtFeeddownXic01(0),
476 fHistoResponseEleLambdaPtFeeddownXic02(0),
477 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
478 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
479 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
480 fHistoResponseEleLambdaPtFeeddownSigma(0),
481 fHistoResponseEleLambdaPtFeeddownSigma1(0),
482 fHistoResponseEleLambdaPtFeeddownSigma2(0),
483 fHistoResponseLambdaPt(0),
484 fHistoResponseLambdaPtFeeddownXic0(0),
485 fHistoResponseLambdaPtFeeddownXicPlus(0),
486 fHistoResponseLambdaPtFeeddownSigma(0),
487 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
488 fHistoEleLambdaPtvsRapidityRS(0),
489 fHistoEleLambdaPtvsRapidityWS(0),
490 fHistoEleLambdaPtvsRapidityMCS(0),
491 fHistoElectronPi0Total(0),
492 fHistoElectronPi0Tag(0),
493 fHistoElectronEtaTotal(0),
494 fHistoElectronEtaTag(0),
496 fHistonEvtvsRunNumber(0),
497 fHistonElevsRunNumber(0),
498 fHistonLambdavsRunNumber(0),
499 fHistoMCEventType(0),
500 fHistoMCDeltaPhiccbar(0),
501 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
502 fHistodPhiSdEtaSElectronProtonR125RS(0),
503 fHistodPhiSdEtaSElectronProtonR125WS(0),
504 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
505 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
506 fHistodPhiSdEtaSElectronPionR125RS(0),
507 fHistodPhiSdEtaSElectronPionR125WS(0),
508 fHistodPhiSdEtaSElectronPionR125RSMix(0),
509 fHistodPhiSdEtaSElectronPionR125WSMix(0),
511 fMixWithoutConversionFlag(kTRUE),
512 fNumberOfEventsForMixing (5),
518 fElectronTracks(0x0),
525 fElectronCutVarsArray(0x0),
526 fV0CutVarsArray1(0x0),
527 fV0CutVarsArray2(0x0)
532 for(Int_t i=0;i<17;i++){
537 for(Int_t i=0;i<8;i++){
545 Bool_t writeVariableTree) :
546 AliAnalysisTaskSE(name),
555 fIsEventSelected(kFALSE),
556 fWriteVariableTree(writeVariableTree),
557 fWriteEachVariableTree(kFALSE),
558 fWriteMCVariableTree(kFALSE),
560 fEleVariablesTree(0),
563 fMCEleVariablesTree(0),
564 fMCV0VariablesTree(0),
565 fMCGenPairVariablesTree(0),
566 fCorrelationVariablesTree(0),
567 fReconstructPrimVert(kFALSE),
573 fCandidateVariables(),
574 fCandidateEleVariables(),
575 fCandidateV0Variables(),
576 fCandidateMCVariables(),
577 fCandidateMCEleVariables(),
578 fCandidateMCV0Variables(),
579 fCandidateMCGenPairVariables(),
580 fCorrelationVariables(),
588 fUseCentralityV0M(kFALSE),
591 fMCDoPairAnalysis(kFALSE),
592 fHistoEleLambdaMass(0),
593 fHistoEleLambdaMassRS(0),
594 fHistoEleLambdaMassRS1(0),
595 fHistoEleLambdaMassRS2(0),
596 fHistoEleLambdaMassWS(0),
597 fHistoEleLambdaMassWS1(0),
598 fHistoEleLambdaMassWS2(0),
599 fHistoEleLambdaMassRSMix(0),
600 fHistoEleLambdaMassRSMix1(0),
601 fHistoEleLambdaMassRSMix2(0),
602 fHistoEleLambdaMassWSMix(0),
603 fHistoEleLambdaMassWSMix1(0),
604 fHistoEleLambdaMassWSMix2(0),
605 fHistoEleLambdaMassRSSide(0),
606 fHistoEleLambdaMassRSSide1(0),
607 fHistoEleLambdaMassRSSide2(0),
608 fHistoEleLambdaMassWSSide(0),
609 fHistoEleLambdaMassWSSide1(0),
610 fHistoEleLambdaMassWSSide2(0),
611 fHistoEleLambdaMassRSAway(0),
612 fHistoEleLambdaMassRS1Away(0),
613 fHistoEleLambdaMassRS2Away(0),
614 fHistoEleLambdaMassWSAway(0),
615 fHistoEleLambdaMassWS1Away(0),
616 fHistoEleLambdaMassWS2Away(0),
617 fHistoEleLambdaMassRSMixAway(0),
618 fHistoEleLambdaMassRSMix1Away(0),
619 fHistoEleLambdaMassRSMix2Away(0),
620 fHistoEleLambdaMassWSMixAway(0),
621 fHistoEleLambdaMassWSMix1Away(0),
622 fHistoEleLambdaMassWSMix2Away(0),
623 fHistoEleLambdaMassRSSideAway(0),
624 fHistoEleLambdaMassRSSide1Away(0),
625 fHistoEleLambdaMassRSSide2Away(0),
626 fHistoEleLambdaMassWSSideAway(0),
627 fHistoEleLambdaMassWSSide1Away(0),
628 fHistoEleLambdaMassWSSide2Away(0),
629 fHistoEleLambdaMassvsElePtRS(0),
630 fHistoEleLambdaMassvsElePtWS(0),
631 fHistoEleLambdaMassvsElePtRSMix(0),
632 fHistoEleLambdaMassvsElePtWSMix(0),
633 fHistoEleLambdaMassvsElePtRSSide(0),
634 fHistoEleLambdaMassvsElePtWSSide(0),
635 fHistoEleLambdaMassvsElePtRS1(0),
636 fHistoEleLambdaMassvsElePtWS1(0),
637 fHistoEleLambdaMassvsElePtRSMix1(0),
638 fHistoEleLambdaMassvsElePtWSMix1(0),
639 fHistoEleLambdaMassvsElePtRSSide1(0),
640 fHistoEleLambdaMassvsElePtWSSide1(0),
641 fHistoEleLambdaMassvsElePtRS2(0),
642 fHistoEleLambdaMassvsElePtWS2(0),
643 fHistoEleLambdaMassvsElePtRSMix2(0),
644 fHistoEleLambdaMassvsElePtWSMix2(0),
645 fHistoEleLambdaMassvsElePtRSSide2(0),
646 fHistoEleLambdaMassvsElePtWSSide2(0),
651 fHistoEleLambdaMassMCS(0),
652 fHistoEleLambdaMassLambdabMCS(0),
653 fHistoEleLambdaMassXibMCS(0),
654 fHistoEleLambdaMassMCGen(0),
655 fHistoEleLambdaMassvsElePtMCS(0),
656 fHistoEleLambdaMassvsElePtMCGen(0),
657 fHistoEleLambdaMassvsElePtMCS1(0),
658 fHistoEleLambdaMassvsElePtMCGen1(0),
659 fHistoEleLambdaMassvsElePtMCS2(0),
660 fHistoEleLambdaMassvsElePtMCGen2(0),
663 fHistoElePtvsEtaRS(0),
664 fHistoElePtvsEtaWS(0),
665 fHistoElePtvsEtaRSMix(0),
666 fHistoElePtvsEtaWSMix(0),
667 fHistoElePtvsEtaMCS(0),
668 fHistoElePtvsEtaMCGen(0),
669 fHistoElePtvsLambdaPtRS(0),
670 fHistoElePtvsLambdaPtWS(0),
671 fHistoElePtvsLambdaPtRSMix(0),
672 fHistoElePtvsLambdaPtWSMix(0),
673 fHistoElePtvsLambdaPtMCS(0),
674 fHistoElePtvsLambdaPtvsLcPtMCS(0),
675 fHistoElePtvsLambdaPtMCGen(0),
676 fHistoElePtvsLambdaPtvsLcPtMCGen(0),
677 fHistoElePtvsLambdaPtMCLcGen(0),
678 fHistoElePtvsd0RS(0),
679 fHistoElePtvsd0WS(0),
680 fHistoElePtvsd0RSMix(0),
681 fHistoElePtvsd0WSMix(0),
682 fHistoElePtvsd0MCS(0),
683 fHistoElePtvsd0PromptMCS(0),
684 fHistoElePtvsd0BFeeddownMCS(0),
685 fHistoEleLambdaMassFeeddownXic0MCS(0),
686 fHistoEleLambdaMassFeeddownXic0MCGen(0),
687 fHistoEleLambdaMassvsElePtFeeddownXic0MCS(0),
688 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen(0),
689 fHistoEleLambdaMassvsElePtFeeddownXic0MCS1(0),
690 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen1(0),
691 fHistoEleLambdaMassvsElePtFeeddownXic0MCS2(0),
692 fHistoEleLambdaMassvsElePtFeeddownXic0MCGen2(0),
693 fHistoElePtFeeddownXic0MCS(0),
694 fHistoElePtFeeddownXic0MCGen(0),
695 fHistoElePtvsEtaFeeddownXic0MCS(0),
696 fHistoElePtvsEtaFeeddownXic0MCGen(0),
697 fHistoElePtvsLambdaPtFeeddownXic0MCS(0),
698 fHistoElePtvsLambdaPtFeeddownXic0MCGen(0),
699 fHistoEleLambdaMassFeeddownXicPlusMCS(0),
700 fHistoEleLambdaMassFeeddownXicPlusMCGen(0),
701 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS(0),
702 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen(0),
703 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS1(0),
704 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen1(0),
705 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCS2(0),
706 fHistoEleLambdaMassvsElePtFeeddownXicPlusMCGen2(0),
707 fHistoElePtFeeddownXicPlusMCS(0),
708 fHistoElePtFeeddownXicPlusMCGen(0),
709 fHistoElePtvsEtaFeeddownXicPlusMCS(0),
710 fHistoElePtvsEtaFeeddownXicPlusMCGen(0),
711 fHistoElePtvsLambdaPtFeeddownXicPlusMCS(0),
712 fHistoElePtvsLambdaPtFeeddownXicPlusMCGen(0),
715 fHistoBachPtMCGen(0),
717 fHistoLambdaMassvsPt(0),
718 fHistoLambdaMassvsPtMCS(0),
719 fHistoLambdaMassvsPtMCGen(0),
720 fHistoLambdaPtvsEta(0),
721 fHistoLambdaPtvsEtaMCS(0),
722 fHistoLambdaPtvsEtaMCGen(0),
723 fHistoLambdaArmenteros(0),
724 fHistoLambdaArmenterosMCS(0),
725 fHistoLambdaPtvsDl(0),
726 fHistoLambdaPtvsDlSide(0),
727 fHistoLambdaPtvsDlMCS(0),
728 fHistoLambdaPtvsDlFeeddownXi0MCS(0),
729 fHistoLambdaPtvsDlFeeddownXiMinusMCS(0),
730 fHistoLambdaPtvsDlFeeddownOmegaMCS(0),
731 fHistoK0sMassvsPt(0),
732 fHistoElectronTPCPID(0),
733 fHistoElectronTOFPID(0),
734 fHistoElectronTPCSelPID(0),
735 fHistoElectronTOFSelPID(0),
736 fHistoElectronTPCPIDSelTOF(0),
737 fHistoElectronTPCPIDSelTOFSmallEta(0),
738 fHistoElectronTPCPIDSelTOFLargeEta(0),
739 fHistoMassConversionsMin(0),
740 fHistoMassConversionsSameSignMin(0),
741 fHistoElectronQovPtvsPhi(0),
742 fHistoLambdaQovPtvsPhi(0),
749 fHistoFeedDownXic0MCGen(0),
750 fHistoFeedDownXic0MCGen1(0),
751 fHistoFeedDownXic0MCGen2(0),
752 fHistoFeedDownXic0MCS(0),
753 fHistoFeedDownXic0MCS1(0),
754 fHistoFeedDownXic0MCS2(0),
755 fHistoFeedDownXicPlusMCGen(0),
756 fHistoFeedDownXicPlusMCGen1(0),
757 fHistoFeedDownXicPlusMCGen2(0),
758 fHistoFeedDownXicPlusMCS(0),
759 fHistoFeedDownXicPlusMCS1(0),
760 fHistoFeedDownXicPlusMCS2(0),
761 fHistoFeedDownSigmaMCGen(0),
762 fHistoFeedDownSigmaMCGen1(0),
763 fHistoFeedDownSigmaMCGen2(0),
764 fHistoFeedDownSigmaMCS(0),
765 fHistoFeedDownSigmaMCS1(0),
766 fHistoFeedDownSigmaMCS2(0),
767 fHistoLcElectronMCGen(0),
768 fHistoLcElectronMCGen1(0),
769 fHistoLcElectronMCGen2(0),
770 fHistoLcElectronMCS(0),
771 fHistoLcElectronMCS1(0),
772 fHistoLcElectronMCS2(0),
773 fHistoElectronFeedDownXic0MCGen(0),
774 fHistoElectronFeedDownXic0MCGen1(0),
775 fHistoElectronFeedDownXic0MCGen2(0),
776 fHistoElectronFeedDownXic0MCS(0),
777 fHistoElectronFeedDownXic0MCS1(0),
778 fHistoElectronFeedDownXic0MCS2(0),
779 fHistoElectronFeedDownXicPlusMCGen(0),
780 fHistoElectronFeedDownXicPlusMCGen1(0),
781 fHistoElectronFeedDownXicPlusMCGen2(0),
782 fHistoElectronFeedDownXicPlusMCS(0),
783 fHistoElectronFeedDownXicPlusMCS1(0),
784 fHistoElectronFeedDownXicPlusMCS2(0),
785 fHistoElectronMCGen(0),
786 fHistoBottomElectronMCGen(0),
787 fHistoCharmElectronMCGen(0),
788 fHistoLambdaMCGen(0),
789 fHistoElePtvsV0dlRS(0),
790 fHistoElePtvsV0dlRS1(0),
791 fHistoElePtvsV0dlRS2(0),
792 fHistoElePtvsV0dlRSSide(0),
793 fHistoElePtvsV0dlRSSide1(0),
794 fHistoElePtvsV0dlRSSide2(0),
795 fHistoElePtvsV0dlRSMix(0),
796 fHistoElePtvsV0dlRSMix1(0),
797 fHistoElePtvsV0dlRSMix2(0),
798 fHistoElePtvsV0dlWS(0),
799 fHistoElePtvsV0dlWS1(0),
800 fHistoElePtvsV0dlWS2(0),
801 fHistoElePtvsV0dlWSSide(0),
802 fHistoElePtvsV0dlWSSide1(0),
803 fHistoElePtvsV0dlWSSide2(0),
804 fHistoElePtvsV0dlWSMix(0),
805 fHistoElePtvsV0dlWSMix1(0),
806 fHistoElePtvsV0dlWSMix2(0),
807 fHistoElePtvsV0dlMCS(0),
808 fHistoElePtvsV0dlMCS1(0),
809 fHistoElePtvsV0dlMCS2(0),
810 fHistoElePtvsV0dlFeedDownXic0MCS(0),
811 fHistoElePtvsV0dlFeedDownXic0MCS1(0),
812 fHistoElePtvsV0dlFeedDownXic0MCS2(0),
813 fHistoElePtvsV0dlFeedDownXicPlusMCS(0),
814 fHistoElePtvsV0dlFeedDownXicPlusMCS1(0),
815 fHistoElePtvsV0dlFeedDownXicPlusMCS2(0),
816 fHistoElePtvsV0dcaRS(0),
817 fHistoElePtvsV0dcaRS1(0),
818 fHistoElePtvsV0dcaRS2(0),
819 fHistoElePtvsV0dcaRSSide(0),
820 fHistoElePtvsV0dcaRSSide1(0),
821 fHistoElePtvsV0dcaRSSide2(0),
822 fHistoElePtvsV0dcaRSMix(0),
823 fHistoElePtvsV0dcaRSMix1(0),
824 fHistoElePtvsV0dcaRSMix2(0),
825 fHistoElePtvsV0dcaWS(0),
826 fHistoElePtvsV0dcaWS1(0),
827 fHistoElePtvsV0dcaWS2(0),
828 fHistoElePtvsV0dcaWSSide(0),
829 fHistoElePtvsV0dcaWSSide1(0),
830 fHistoElePtvsV0dcaWSSide2(0),
831 fHistoElePtvsV0dcaWSMix(0),
832 fHistoElePtvsV0dcaWSMix1(0),
833 fHistoElePtvsV0dcaWSMix2(0),
834 fHistoElePtvsV0dcaMCS(0),
835 fHistoElePtvsV0dcaMCS1(0),
836 fHistoElePtvsV0dcaMCS2(0),
837 fHistoElePtvsV0dcaFeedDownXic0MCS(0),
838 fHistoElePtvsV0dcaFeedDownXic0MCS1(0),
839 fHistoElePtvsV0dcaFeedDownXic0MCS2(0),
840 fHistoElePtvsV0dcaFeedDownXicPlusMCS(0),
841 fHistoElePtvsV0dcaFeedDownXicPlusMCS1(0),
842 fHistoElePtvsV0dcaFeedDownXicPlusMCS2(0),
843 fHistoEleLambdaPtvsV0dlRS(0),
844 fHistoEleLambdaPtvsV0dlRS1(0),
845 fHistoEleLambdaPtvsV0dlRS2(0),
846 fHistoEleLambdaPtvsV0dlRSSide(0),
847 fHistoEleLambdaPtvsV0dlRSSide1(0),
848 fHistoEleLambdaPtvsV0dlRSSide2(0),
849 fHistoEleLambdaPtvsV0dlRSMix(0),
850 fHistoEleLambdaPtvsV0dlRSMix1(0),
851 fHistoEleLambdaPtvsV0dlRSMix2(0),
852 fHistoEleLambdaPtvsV0dlWS(0),
853 fHistoEleLambdaPtvsV0dlWS1(0),
854 fHistoEleLambdaPtvsV0dlWS2(0),
855 fHistoEleLambdaPtvsV0dlWSSide(0),
856 fHistoEleLambdaPtvsV0dlWSSide1(0),
857 fHistoEleLambdaPtvsV0dlWSSide2(0),
858 fHistoEleLambdaPtvsV0dlWSMix(0),
859 fHistoEleLambdaPtvsV0dlWSMix1(0),
860 fHistoEleLambdaPtvsV0dlWSMix2(0),
861 fHistoEleLambdaPtvsV0dlMCS(0),
862 fHistoEleLambdaPtvsV0dlMCS1(0),
863 fHistoEleLambdaPtvsV0dlMCS2(0),
864 fHistoEleLambdaPtvsV0dlPromptMCS(0),
865 fHistoEleLambdaPtvsV0dlPromptMCS1(0),
866 fHistoEleLambdaPtvsV0dlPromptMCS2(0),
867 fHistoEleLambdaPtvsV0dlBFeeddownMCS(0),
868 fHistoEleLambdaPtvsV0dlBFeeddownMCS1(0),
869 fHistoEleLambdaPtvsV0dlBFeeddownMCS2(0),
870 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS(0),
871 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS1(0),
872 fHistoEleLambdaPtvsV0dlFeedDownXic0MCS2(0),
873 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS(0),
874 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS1(0),
875 fHistoEleLambdaPtvsV0dlFeedDownXic0PromptMCS2(0),
876 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS(0),
877 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS1(0),
878 fHistoEleLambdaPtvsV0dlFeedDownXic0BFeeddownMCS2(0),
879 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS(0),
880 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS1(0),
881 fHistoEleLambdaPtvsV0dlFeedDownXicPlusMCS2(0),
882 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS(0),
883 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS1(0),
884 fHistoEleLambdaPtvsV0dlFeedDownXicPlusPromptMCS2(0),
885 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS(0),
886 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS1(0),
887 fHistoEleLambdaPtvsV0dlFeedDownXicPlusBFeeddownMCS2(0),
888 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS(0),
889 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS1(0),
890 fHistoEleLambdaPtvsV0dlFeedDownSigmaMCS2(0),
891 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS(0),
892 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS1(0),
893 fHistoEleLambdaPtvsV0dlFeedDownSigmaPromptMCS2(0),
894 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS(0),
895 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS1(0),
896 fHistoEleLambdaPtvsV0dlFeedDownSigmaBFeeddownMCS2(0),
897 fHistoEleLambdaPtvsV0dlLambdabMCS(0),
898 fHistoEleLambdaPtvsV0dlLambdabMCS1(0),
899 fHistoEleLambdaPtvsV0dlLambdabMCS2(0),
900 fHistoEleLambdaPtvsV0dlXibMCS(0),
901 fHistoEleLambdaPtvsV0dlXibMCS1(0),
902 fHistoEleLambdaPtvsV0dlXibMCS2(0),
903 fHistoEleLambdaPtvsV0dlRSAway(0),
904 fHistoEleLambdaPtvsV0dlRS1Away(0),
905 fHistoEleLambdaPtvsV0dlRS2Away(0),
906 fHistoEleLambdaPtvsV0dlRSSideAway(0),
907 fHistoEleLambdaPtvsV0dlRSSide1Away(0),
908 fHistoEleLambdaPtvsV0dlRSSide2Away(0),
909 fHistoEleLambdaPtvsV0dlRSMixAway(0),
910 fHistoEleLambdaPtvsV0dlRSMix1Away(0),
911 fHistoEleLambdaPtvsV0dlRSMix2Away(0),
912 fHistoEleLambdaPtvsV0dlWSAway(0),
913 fHistoEleLambdaPtvsV0dlWS1Away(0),
914 fHistoEleLambdaPtvsV0dlWS2Away(0),
915 fHistoEleLambdaPtvsV0dlWSSideAway(0),
916 fHistoEleLambdaPtvsV0dlWSSide1Away(0),
917 fHistoEleLambdaPtvsV0dlWSSide2Away(0),
918 fHistoEleLambdaPtvsV0dlWSMixAway(0),
919 fHistoEleLambdaPtvsV0dlWSMix1Away(0),
920 fHistoEleLambdaPtvsV0dlWSMix2Away(0),
921 fHistoResponseElePt(0),
922 fHistoResponseElePt1(0),
923 fHistoResponseElePt2(0),
924 fHistoResponseEleLambdaPt(0),
925 fHistoResponseEleLambdaPt1(0),
926 fHistoResponseEleLambdaPt2(0),
927 fHistoResponseEleLambdaPtFeeddownXic0(0),
928 fHistoResponseEleLambdaPtFeeddownXic01(0),
929 fHistoResponseEleLambdaPtFeeddownXic02(0),
930 fHistoResponseEleLambdaPtFeeddownXicPlus(0),
931 fHistoResponseEleLambdaPtFeeddownXicPlus1(0),
932 fHistoResponseEleLambdaPtFeeddownXicPlus2(0),
933 fHistoResponseEleLambdaPtFeeddownSigma(0),
934 fHistoResponseEleLambdaPtFeeddownSigma1(0),
935 fHistoResponseEleLambdaPtFeeddownSigma2(0),
936 fHistoResponseLambdaPt(0),
937 fHistoResponseLambdaPtFeeddownXic0(0),
938 fHistoResponseLambdaPtFeeddownXicPlus(0),
939 fHistoResponseLambdaPtFeeddownSigma(0),
940 fHistoLcPtvseleLambdaPtvsElePtvsLambdaPt(0),
941 fHistoEleLambdaPtvsRapidityRS(0),
942 fHistoEleLambdaPtvsRapidityWS(0),
943 fHistoEleLambdaPtvsRapidityMCS(0),
944 fHistoElectronPi0Total(0),
945 fHistoElectronPi0Tag(0),
946 fHistoElectronEtaTotal(0),
947 fHistoElectronEtaTag(0),
949 fHistonEvtvsRunNumber(0),
950 fHistonElevsRunNumber(0),
951 fHistonLambdavsRunNumber(0),
952 fHistoMCEventType(0),
953 fHistoMCDeltaPhiccbar(0),
954 fGTI(0),fGTIndex(0), fTrackBuffSize(19000),
955 fHistodPhiSdEtaSElectronProtonR125RS(0),
956 fHistodPhiSdEtaSElectronProtonR125WS(0),
957 fHistodPhiSdEtaSElectronProtonR125RSMix(0),
958 fHistodPhiSdEtaSElectronProtonR125WSMix(0),
959 fHistodPhiSdEtaSElectronPionR125RS(0),
960 fHistodPhiSdEtaSElectronPionR125WS(0),
961 fHistodPhiSdEtaSElectronPionR125RSMix(0),
962 fHistodPhiSdEtaSElectronPionR125WSMix(0),
964 fMixWithoutConversionFlag(kTRUE),
965 fNumberOfEventsForMixing (5),
971 fElectronTracks(0x0),
978 fElectronCutVarsArray(0x0),
979 fV0CutVarsArray1(0x0),
980 fV0CutVarsArray2(0x0)
985 Info(
"AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Constructor");
987 for(Int_t i=0;i<17;i++){
992 for(Int_t i=0;i<8;i++){
996 DefineOutput(1,TList::Class());
997 DefineOutput(2,TList::Class());
998 DefineOutput(3,TList::Class());
999 DefineOutput(4,TTree::Class());
1000 DefineOutput(5,TTree::Class());
1001 DefineOutput(6,TTree::Class());
1002 DefineOutput(7,TTree::Class());
1003 DefineOutput(8,AliNormalizationCounter::Class());
1004 DefineOutput(9,TTree::Class());
1005 DefineOutput(10,TTree::Class());
1006 DefineOutput(11,TTree::Class());
1007 DefineOutput(12,TTree::Class());
1015 Info(
"~AliAnalysisTaskSELc2eleLambdafromAODtracks",
"Calling Destructor");
1110 if (fDebug > 1) AliInfo(
"Init");
1129 AliError(
"NO EVENT FOUND!");
1132 AliAODEvent* aodEvent =
dynamic_cast<AliAODEvent*
>(fInputEvent);
1140 fBzkG = (Double_t)aodEvent->GetMagneticField();
1141 AliKFParticle::SetField(
fBzkG);
1142 if (TMath::Abs(
fBzkG)<0.001) {
1153 TClonesArray *mcArray = 0;
1154 AliAODMCHeader *mcHeader=0;
1157 mcArray =
dynamic_cast<TClonesArray*
>(aodEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1159 AliError(
"Could not find Monte-Carlo in AOD");
1165 mcHeader = (AliAODMCHeader*)aodEvent->GetList()->FindObject(AliAODMCHeader::StdBranchName());
1167 AliError(
"AliAnalysisTaskSELc2eleLambdafromAODtracks::UserExec: MC header branch not found!\n");
1172 Double_t zMCVertex = mcHeader->GetVtxZ();
1174 AliDebug(2,Form(
"Event rejected: abs(zVtxMC)=%f > fAnalCuts->GetMaxVtxZ()=%f",zMCVertex,
fAnalCuts->
GetMaxVtxZ()));
1188 fVtx1 = (AliAODVertex*)aodEvent->GetPrimaryVertex();
1191 Double_t pos[3],cov[6];
1193 fVtx1->GetCovarianceMatrix(cov);
1194 fV1 =
new AliESDVertex(pos,cov,100.,100,
fVtx1->GetName());
1198 if(!fIsTriggerNotOK)
fCEvents->Fill(3);
1205 fIsMB=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kMB)==(AliVEvent::kMB);
1206 fIsSemi=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kSemiCentral)==(AliVEvent::kSemiCentral);
1207 fIsCent=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kCentral)==(AliVEvent::kCentral);
1208 fIsINT7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kINT7)==(AliVEvent::kINT7);
1209 fIsEMC7=(((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected()&AliVEvent::kEMC7)==(AliVEvent::kEMC7);
1223 AliCentrality *cent = aodEvent->GetCentrality();
1224 fCentrality = cent->GetCentralityPercentile(
"V0M");
1228 if(fCentrality<0.||fCentrality>100.-0.0000001) {
1235 Int_t runnumber_offset = 0;
1236 Int_t runnumber = aodEvent->GetRunNumber();
1237 if(runnumber<=131000&&runnumber>=114000){
1238 runnumber_offset = 114000;
1239 }
else if(runnumber<=196000&&runnumber>=195000){
1240 runnumber_offset = 195000;
1241 }
else if(runnumber<=170593&&runnumber>=167902){
1242 runnumber_offset = 167902;
1287 AliAnalysisTaskSE::Terminate();
1289 fOutput =
dynamic_cast<TList*
> (GetOutputData(1));
1291 AliError(
"fOutput not available");
1295 fOutputAll =
dynamic_cast<TList*
> (GetOutputData(3));
1297 AliError(
"fOutputAll not available");
1352 TString normName=
"NormalizationCounter";
1353 AliAnalysisDataContainer *cont = GetOutputSlot(8)->GetContainer();
1354 if(cont)normName=(TString)cont->GetName();
1376 fEventBuffer[i]=
new TTree(Form(
"EventBuffer_%d",i),
"Temporary buffer for event mixing");
1400 AliAODEvent *aodEvent, TClonesArray *mcArray
1407 if(fElectronTracks) fElectronTracks->Delete();
1408 if(fV0Tracks1) fV0Tracks1->Delete();
1409 if(fV0Tracks2) fV0Tracks2->Delete();
1410 if(fElectronCutVarsArray) fElectronCutVarsArray->Delete();
1411 if(fV0CutVarsArray1) fV0CutVarsArray1->Delete();
1412 if(fV0CutVarsArray2) fV0CutVarsArray2->Delete();
1413 fV0dlArray1.clear();
1414 fV0dlArray2.clear();
1415 fV0dcaArray1.clear();
1416 fV0dcaArray2.clear();
1419 ResetGlobalTrackReference();
1421 for (Int_t iTrack=0;iTrack<aodEvent->GetNumberOfTracks();iTrack++){
1423 AliAODTrack *track =
dynamic_cast<AliAODTrack *
>(aodEvent->GetTrack(iTrack));
1424 if (!track)
continue;
1427 StoreGlobalTrackReference(track,iTrack);
1434 Int_t nV0s= aodEvent->GetNumberOfV0s();
1435 Int_t nTracks= aodEvent->GetNumberOfTracks();
1437 Bool_t seleTrkFlags[nTracks];
1439 SelectTrack(aodEvent,nTracks,nSeleTrks,seleTrkFlags,mcArray);
1441 Bool_t seleV0Flags[nV0s];
1443 SelectV0(aodEvent,nV0s,nSeleV0,seleV0Flags,mcArray);
1445 Int_t runnumber_offset = 0;
1446 Int_t runnumber = aodEvent->GetRunNumber();
1447 if(runnumber<=131000&&runnumber>=114000){
1448 runnumber_offset = 114000;
1449 }
else if(runnumber<=196000&&runnumber>=195000){
1450 runnumber_offset = 195000;
1451 }
else if(runnumber<=170593&&runnumber>=167902){
1452 runnumber_offset = 167902;
1454 fHistonElevsRunNumber->Fill(runnumber-runnumber_offset,nSeleTrks);
1455 fHistonLambdavsRunNumber->Fill(runnumber-runnumber_offset,nSeleV0);
1460 for (Int_t iv0 = 0; iv0<nV0s; iv0++) {
1461 if(!seleV0Flags[iv0])
continue;
1462 AliAODv0 *v0 = aodEvent->GetV0(iv0);
1465 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1466 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1471 for (Int_t itrk = 0; itrk<nTracks; itrk++) {
1472 if(!seleTrkFlags[itrk])
continue;
1473 AliAODTrack *trk = (AliAODTrack*)aodEvent->GetTrack(itrk);
1481 AliAODTrack *trkpid = 0;
1482 if(fAnalCuts->GetProdAODFilterBit()==7){
1483 trkpid = fGTI[-trk->GetID()-1];
1488 Int_t cpid = cptrack->GetID();
1489 Int_t cnid = cntrack->GetID();
1490 Int_t lpid = trkpid->GetID();
1491 if((cpid==lpid)||(cnid==lpid))
continue;
1493 AliAODVertex *secVert = ReconstructSecondaryVertex(v0,trk,aodEvent);
1494 if(!secVert)
continue;
1501 FillROOTObjects(elobj,v0,trk,trkpid,aodEvent,mcArray);
1503 elobj->GetSecondaryVtx()->RemoveDaughters();
1505 delete elobj;elobj=NULL;
1511 fEventInfo->SetString(Form(
"Ev%d_esd%d_E%d_V%d",AliAnalysisManager::GetAnalysisManager()->GetNcalls(),((AliAODHeader*)aodEvent->GetHeader())->GetEventNumberESDFile(),fElectronTracks->GetEntries(),fV0Tracks1->GetEntries()+fV0Tracks2->GetEntries()));
1512 Int_t ind=GetPoolIndex(fVtxZ,fCentrality);
1513 if(ind>=0 && ind<fNOfPools){
1514 if(fEventBuffer[ind]->GetEntries() >= fNumberOfEventsForMixing){
1515 DoEventMixingWithPools(ind);
1516 if(fEventBuffer[ind]->GetEntries() >= 20*fNumberOfEventsForMixing){
1520 fEventBuffer[ind]->Fill();
1531 const char* nameoutput = GetOutputSlot(4)->GetContainer()->GetName();
1532 fVariablesTree =
new TTree(nameoutput,
"Candidates variables tree");
1535 TString * fCandidateVariableNames =
new TString[nVar];
1537 fCandidateVariableNames[ 0]=
"Centrality";
1538 fCandidateVariableNames[ 1]=
"InvMassEleLambda";
1539 fCandidateVariableNames[ 2]=
"EleLambdaPt";
1540 fCandidateVariableNames[ 3]=
"EleLambdaPx";
1541 fCandidateVariableNames[ 4]=
"EleLambdaPy";
1542 fCandidateVariableNames[ 5]=
"EleLambdaPz";
1543 fCandidateVariableNames[ 6]=
"ElePx";
1544 fCandidateVariableNames[ 7]=
"ElePy";
1545 fCandidateVariableNames[ 8]=
"ElePz";
1546 fCandidateVariableNames[ 9]=
"V0Px";
1547 fCandidateVariableNames[10]=
"V0Py";
1548 fCandidateVariableNames[11]=
"V0Pz";
1549 fCandidateVariableNames[12]=
"AntiLambdaFlag";
1550 fCandidateVariableNames[13]=
"MassLambda";
1551 fCandidateVariableNames[14]=
"MassAntiLambda";
1552 fCandidateVariableNames[15]=
"Eled0";
1553 fCandidateVariableNames[16]=
"V0d0";
1554 fCandidateVariableNames[17]=
"nSigmaTPCele";
1555 fCandidateVariableNames[18]=
"nSigmaTOFele";
1556 fCandidateVariableNames[19]=
"nSigmaTPCv0pr";
1557 fCandidateVariableNames[20]=
"nSigmaTOFv0pr";
1558 fCandidateVariableNames[21]=
"EleCharge";
1559 fCandidateVariableNames[22]=
"ProtonPx";
1560 fCandidateVariableNames[23]=
"ProtonPy";
1561 fCandidateVariableNames[24]=
"ProtonPz";
1562 fCandidateVariableNames[25]=
"PiPx";
1563 fCandidateVariableNames[26]=
"PiPy";
1564 fCandidateVariableNames[27]=
"PiPz";
1565 fCandidateVariableNames[28]=
"mcpdglc";
1566 fCandidateVariableNames[29]=
"mclablc";
1567 fCandidateVariableNames[30]=
"mcpdgmomele";
1568 fCandidateVariableNames[31]=
"mcpdgmomv0";
1569 fCandidateVariableNames[32]=
"Mixing";
1570 fCandidateVariableNames[33]=
"mcpdgele";
1571 fCandidateVariableNames[34]=
"nSigmaTPCpr_etrk";
1572 fCandidateVariableNames[35]=
"nSigmaTOFpr_etrk";
1573 fCandidateVariableNames[36]=
"nSigmaTPCka_etrk";
1574 fCandidateVariableNames[37]=
"nSigmaTOFka_etrk";
1575 fCandidateVariableNames[38]=
"MassK0Short";
1576 fCandidateVariableNames[39]=
"mcpdggrmomele";
1577 fCandidateVariableNames[40]=
"mcpdggrmomv0";
1578 fCandidateVariableNames[41]=
"mcngenele";
1579 fCandidateVariableNames[42]=
"mcngenv0";
1580 fCandidateVariableNames[43]=
"mclcpx";
1581 fCandidateVariableNames[44]=
"mclcpy";
1582 fCandidateVariableNames[45]=
"mclcpz";
1583 fCandidateVariableNames[46]=
"mcelepx";
1584 fCandidateVariableNames[47]=
"mcelepy";
1585 fCandidateVariableNames[48]=
"mcelepz";
1586 fCandidateVariableNames[49]=
"mcv0px";
1587 fCandidateVariableNames[50]=
"mcv0py";
1588 fCandidateVariableNames[51]=
"mcv0pz";
1589 fCandidateVariableNames[52]=
"nSigmaTPCpi_etrk";
1590 fCandidateVariableNames[53]=
"nSigmaTOFpi_etrk";
1591 fCandidateVariableNames[54]=
"PrimVertx";
1592 fCandidateVariableNames[55]=
"PrimVerty";
1593 fCandidateVariableNames[56]=
"PrimVertz";
1594 fCandidateVariableNames[57]=
"V0Vertx";
1595 fCandidateVariableNames[58]=
"V0Verty";
1596 fCandidateVariableNames[59]=
"V0Vertz";
1598 fCandidateVariableNames[60]=
"DcaV0PrToPrimVertex";
1599 fCandidateVariableNames[61]=
"DcaV0PiToPrimVertex";
1600 fCandidateVariableNames[62]=
"DcaV0daughters";
1601 fCandidateVariableNames[63]=
"V0CosPointingAngle";
1602 fCandidateVariableNames[64]=
"V0ProperDecayLength";
1603 fCandidateVariableNames[65]=
"MassK0Short2";
1605 fCandidateVariableNames[66]=
"nSigmaTPCv0pi";
1606 fCandidateVariableNames[67]=
"nSigmaTOFv0pi";
1608 fCandidateVariableNames[68]=
"EleITSMatch";
1609 fCandidateVariableNames[69]=
"V0PosITSMatch";
1610 fCandidateVariableNames[70]=
"V0NegITSMatch";
1611 fCandidateVariableNames[71]=
"IsV0PeakRegion";
1612 fCandidateVariableNames[72]=
"mcpdgv0";
1614 fCandidateVariableNames[73]=
"mcv0pospx";
1615 fCandidateVariableNames[74]=
"mcv0pospy";
1616 fCandidateVariableNames[75]=
"mcv0pospz";
1617 fCandidateVariableNames[76]=
"mcpdgv0pos";
1618 fCandidateVariableNames[77]=
"mcv0negpx";
1619 fCandidateVariableNames[78]=
"mcv0negpy";
1620 fCandidateVariableNames[79]=
"mcv0negpz";
1621 fCandidateVariableNames[80]=
"mcpdgv0neg";
1623 fCandidateVariableNames[81]=
"ProtonPxPV";
1624 fCandidateVariableNames[82]=
"ProtonPyPV";
1625 fCandidateVariableNames[83]=
"ProtonPzPV";
1626 fCandidateVariableNames[84]=
"PiPxPV";
1627 fCandidateVariableNames[85]=
"PiPyPV";
1628 fCandidateVariableNames[86]=
"PiPzPV";
1630 fCandidateVariableNames[87]=
"mcv0vertx";
1631 fCandidateVariableNames[88]=
"mcv0verty";
1632 fCandidateVariableNames[89]=
"mcv0vertz";
1634 fCandidateVariableNames[90]=
"OnFlyStatus";
1635 fCandidateVariableNames[91]=
"MagneticField";
1636 fCandidateVariableNames[92]=
"EvNumber";
1637 fCandidateVariableNames[93]=
"RunNumber";
1640 for (Int_t ivar=0; ivar<nVar; ivar++) {
1657 for(Int_t i=0;i<94;i++)
1662 Bool_t anti_lambda_flag = kFALSE;
1665 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
1666 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
1672 Double_t v0px = elobj->PxProng(1);
1673 Double_t v0py = elobj->PyProng(1);
1674 Double_t v0pz = elobj->PzProng(1);
1675 Double_t momv0 = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz);
1676 Double_t Ev0 = sqrt(momv0*momv0+1.115683*1.115683);
1677 Double_t epx = elobj->PxProng(0);
1678 Double_t epy = elobj->PyProng(0);
1679 Double_t epz = elobj->PzProng(0);
1680 Double_t mome = sqrt(epx*epx+epy*epy+epz*epz);
1681 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
1682 Double_t cosoa = (v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz);
1683 Double_t Esum = Ee + Ev0;
1685 Double_t uxe = epx/mome;
1686 Double_t uye = epy/mome;
1687 Double_t uze = epz/mome;
1688 Double_t lf = -2.*(v0px*uxe+v0py*uye+v0pz*uze);
1689 Double_t pxv_flip = v0px + lf * uxe;
1690 Double_t pyv_flip = v0py + lf * uye;
1691 Double_t pzv_flip = v0pz + lf * uze;
1692 Double_t pxsum_flip = epx + pxv_flip;
1693 Double_t pysum_flip = epy + pyv_flip;
1694 Double_t pzsum_flip = epz + pzv_flip;
1695 Double_t melam_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
1696 Double_t ptelam_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
1698 Double_t minmass_ee = 9999.;
1700 Double_t minmasslike_ee = 9999.;
1704 UInt_t pdgdg[2]={11,3122};
1722 Double_t nSigmaTPCele=-9999.;
1723 Double_t nSigmaTOFele=-9999.;
1724 Double_t nSigmaTPCv0pr=-9999.;
1725 Double_t nSigmaTOFv0pr=-9999.;
1726 Double_t nSigmaTPCv0pi=-9999.;
1727 Double_t nSigmaTOFv0pi=-9999.;
1738 if(anti_lambda_flag){
1756 if(anti_lambda_flag){
1772 AliAODMCParticle *mclc = 0;
1773 AliAODMCParticle *mcele = 0;
1774 AliAODMCParticle *mcv0 = 0;
1776 Int_t mcpdgele_array[100];
1777 Int_t mcpdgv0_array[100];
1778 Int_t mclabelele_array[100];
1779 Int_t mclabelv0_array[100];
1780 Int_t mcngen_ele=-9999;
1781 Int_t mcngen_v0=-9999;
1784 mclablc =
MatchToMC(elobj,mcArray,mcpdgele_array, mcpdgv0_array,mclabelele_array,mclabelv0_array,mcngen_ele,mcngen_v0);
1785 if(mclabelele_array[0]>=0){
1786 mcele = (AliAODMCParticle*) mcArray->At(mclabelele_array[0]);
1793 if(mclabelv0_array[0]>=0){
1794 mcv0 = (AliAODMCParticle*) mcArray->At(mclabelv0_array[0]);
1802 Int_t pdgdgv0[2]={211,211};
1803 Int_t labV0 = v0->MatchToMC(310,mcArray,2,pdgdgv0);
1805 mcv0 = (AliAODMCParticle*) mcArray->At(labV0);
1815 Int_t labPos = cptrack->GetLabel();
1817 AliAODMCParticle *mcv0pos = (AliAODMCParticle*)mcArray->At(labPos);
1826 Int_t labNeg = cntrack->GetLabel();
1828 AliAODMCParticle *mcv0neg = (AliAODMCParticle*)mcArray->At(labNeg);
1836 mclc = (AliAODMCParticle*) mcArray->At(mclablc);
1879 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
1880 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
1881 if(!anti_lambda_flag){
1889 Double_t posVtx[3] = {0.,0.,0.};
1890 fVtx1->GetXYZ(posVtx);
1892 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
1902 if(!anti_lambda_flag){
1927 Double_t dphis_ele_pr, detas_ele_pr,dphis_ele_pi, detas_ele_pi;
1928 dphis_ele_pr = 9999.;detas_ele_pr = 9999.;dphis_ele_pi = 9999.;detas_ele_pi = 9999.;
1932 cont[0] = elobj->InvMass(2,pdgdg);
1933 cont[1] = elobj->Pt();
1937 Double_t cont_flip[3];
1938 cont_flip[0] = melam_flip;
1939 cont_flip[1] = ptelam_flip;
1943 cont2[0] = elobj->InvMass(2,pdgdg);
1944 cont2[1] = trk->Pt();
1946 Double_t cont_eleptvseta[3];
1947 cont_eleptvseta[0] = trk->Pt();
1948 cont_eleptvseta[1] = trk->Eta();
1951 Double_t cont_eleptvslambdapt[3];
1952 cont_eleptvslambdapt[0] = trk->Pt();
1953 cont_eleptvslambdapt[1] = v0->Pt();
1956 Double_t cont_eleptvsd0[3];
1957 cont_eleptvsd0[0] = trk->Pt();
1958 cont_eleptvsd0[1] = elobj->Getd0Prong(0)*trk->Charge();
1961 Double_t cont_eleptvsv0dl[3];
1962 cont_eleptvsv0dl[0] = trk->Pt();
1963 cont_eleptvsv0dl[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
1966 Double_t cont_elelamptvsv0dl[4];
1967 cont_elelamptvsv0dl[0] = elobj->Pt();
1968 cont_elelamptvsv0dl[1] = v0->
DecayLengthV0(posVtx)*1.115683/ptotlam;
1969 cont_elelamptvsv0dl[2] = elobj->Getd0Prong(0)*trk->Charge();
1972 Double_t cont_elelamptvsv0dl_flip[4];
1973 cont_elelamptvsv0dl_flip[0] = ptelam_flip;
1974 cont_elelamptvsv0dl_flip[1] = v0->DecayLengthV0(posVtx)*1.115683/ptotlam;
1975 cont_elelamptvsv0dl_flip[2] = 0.0;
1979 Double_t cont_eleptvsv0dca[3];
1980 cont_eleptvsv0dca[0] = trk->Pt();
1981 cont_eleptvsv0dca[1] = v0->DcaV0ToPrimVertex();
1984 Double_t elobj_mass = elobj->InvMass(2,pdgdg);
1985 Double_t elobj_px = elobj->Px();
1986 Double_t elobj_py = elobj->Py();
1987 Double_t elobj_pz = elobj->Pz();
1988 Double_t elobj_E = sqrt(elobj_mass*elobj_mass+elobj_px*elobj_px+elobj_py*elobj_py+elobj_pz*elobj_pz);
1989 Double_t elobj_rap = 0.5*log((elobj_E+elobj_pz)/(elobj_E-elobj_pz));
1996 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2023 for(Int_t ih=0;ih<17;ih++){
2024 Double_t cont_eleptvscutvars[3];
2025 cont_eleptvscutvars[0] = elobj->Pt();
2029 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2031 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2033 cont_eleptvscutvars[1] = nSigmaTPCele;
2035 cont_eleptvscutvars[1] = nSigmaTOFele;
2037 cont_eleptvscutvars[1] = trk->Eta();
2039 cont_eleptvscutvars[1] = trk->GetITSNcls();
2041 if(!anti_lambda_flag)
2042 cont_eleptvscutvars[1] = v0->MassLambda();
2044 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2047 lPosV0[0] = v0->DecayVertexV0X();
2048 lPosV0[1] = v0->DecayVertexV0Y();
2049 lPosV0[2] = v0->DecayVertexV0Z();
2050 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2052 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2054 if(!anti_lambda_flag)
2055 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2057 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2059 if(!anti_lambda_flag)
2060 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2062 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2064 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2066 cont_eleptvscutvars[1] = v0->MassK0Short();
2068 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2070 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2072 cont_eleptvscutvars[1] = v0->Eta();
2074 Double_t v0px = elobj->PxProng(1);
2075 Double_t v0py = elobj->PyProng(1);
2076 Double_t v0pz = elobj->PzProng(1);
2077 Double_t epx = elobj->PxProng(0);
2078 Double_t epy = elobj->PyProng(0);
2079 Double_t epz = elobj->PzProng(0);
2080 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2082 cont_eleptvscutvars[1] = -9999.;
2090 Int_t pdgcode = mclc->GetPdgCode();
2091 Double_t cont_mclc[3];
2092 cont_mclc[0] = mclc->Pt();
2093 cont_mclc[1] = mclc->Y();
2095 Double_t cont_mcele[3];
2096 cont_mcele[0] = mcele->Pt();
2097 cont_mcele[1] = mcele->Eta();
2100 Int_t labmotherlc = mclc->GetMother();
2101 Bool_t isbottomfd = kFALSE;
2103 AliAODMCParticle *motherlc = (AliAODMCParticle*)mcArray->At(labmotherlc);
2104 Int_t pdgmotherlc = motherlc->GetPdgCode();
2105 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){
2109 cont2[1] = mcele->Pt();
2111 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2162 Double_t cont_eleptvslambdaptvslcpt[4];
2163 cont_eleptvslambdaptvslcpt[0] = cont_eleptvslambdapt[0];
2164 cont_eleptvslambdaptvslcpt[1] = cont_eleptvslambdapt[1];
2165 cont_eleptvslambdaptvslcpt[2] = mclc->Pt();
2166 cont_eleptvslambdaptvslcpt[3] = cont_eleptvslambdapt[2];
2169 Double_t cont_allpt[4];
2170 cont_allpt[0] = mclc->Pt();
2171 cont_allpt[1] = elobj->Pt();
2172 cont_allpt[2] = trk->Pt();
2173 cont_allpt[3] = v0->Pt();
2182 for(Int_t ih=0;ih<17;ih++){
2183 Double_t cont_eleptvscutvars[3];
2184 cont_eleptvscutvars[0] = elobj->Pt();
2188 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2190 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2192 cont_eleptvscutvars[1] = nSigmaTPCele;
2194 cont_eleptvscutvars[1] = nSigmaTOFele;
2196 cont_eleptvscutvars[1] = trk->Eta();
2198 cont_eleptvscutvars[1] = trk->GetITSNcls();
2200 if(!anti_lambda_flag)
2201 cont_eleptvscutvars[1] = v0->MassLambda();
2203 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2206 lPosV0[0] = v0->DecayVertexV0X();
2207 lPosV0[1] = v0->DecayVertexV0Y();
2208 lPosV0[2] = v0->DecayVertexV0Z();
2209 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2211 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2213 if(!anti_lambda_flag)
2214 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2216 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2218 if(!anti_lambda_flag)
2219 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2221 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2223 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2225 cont_eleptvscutvars[1] = v0->MassK0Short();
2227 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2229 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2231 cont_eleptvscutvars[1] = v0->Eta();
2233 Double_t v0px = elobj->PxProng(1);
2234 Double_t v0py = elobj->PyProng(1);
2235 Double_t v0pz = elobj->PzProng(1);
2236 Double_t epx = elobj->PxProng(0);
2237 Double_t epy = elobj->PyProng(0);
2238 Double_t epz = elobj->PzProng(0);
2239 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2241 cont_eleptvscutvars[1] = -9999.;
2248 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2295 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2342 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2372 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2399 for(Int_t ih=0;ih<17;ih++){
2400 Double_t cont_eleptvscutvars[3];
2401 cont_eleptvscutvars[0] = elobj->Pt();
2405 cont_eleptvscutvars[1] = trk->GetTPCNcls();
2407 cont_eleptvscutvars[1] = trk->GetTPCsignalN();
2409 cont_eleptvscutvars[1] = nSigmaTPCele;
2411 cont_eleptvscutvars[1] = nSigmaTOFele;
2413 cont_eleptvscutvars[1] = trk->Eta();
2415 cont_eleptvscutvars[1] = trk->GetITSNcls();
2417 if(!anti_lambda_flag)
2418 cont_eleptvscutvars[1] = v0->MassLambda();
2420 cont_eleptvscutvars[1] = v0->MassAntiLambda();
2423 lPosV0[0] = v0->DecayVertexV0X();
2424 lPosV0[1] = v0->DecayVertexV0Y();
2425 lPosV0[2] = v0->DecayVertexV0Z();
2426 cont_eleptvscutvars[1] = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
2428 cont_eleptvscutvars[1] = v0->DcaV0Daughters();
2430 if(!anti_lambda_flag)
2431 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2433 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2435 if(!anti_lambda_flag)
2436 cont_eleptvscutvars[1] = v0->DcaNegToPrimVertex();
2438 cont_eleptvscutvars[1] = v0->DcaPosToPrimVertex();
2440 cont_eleptvscutvars[1] = v0->CosPointingAngle(posVtx);
2442 cont_eleptvscutvars[1] = v0->MassK0Short();
2444 cont_eleptvscutvars[1] = nSigmaTPCv0pr;
2446 cont_eleptvscutvars[1] = nSigmaTPCv0pi;
2448 cont_eleptvscutvars[1] = v0->Eta();
2450 Double_t v0px = elobj->PxProng(1);
2451 Double_t v0py = elobj->PyProng(1);
2452 Double_t v0pz = elobj->PzProng(1);
2453 Double_t epx = elobj->PxProng(0);
2454 Double_t epy = elobj->PyProng(0);
2455 Double_t epz = elobj->PzProng(0);
2456 cont_eleptvscutvars[1] = acos((v0px*epx+v0py*epy+v0pz*epz)/sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz)/sqrt(epx*epx+epy*epy+epz*epz));
2458 cont_eleptvscutvars[1] = -9999.;
2466 Int_t pdgcode = mclc->GetPdgCode();
2467 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2474 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2480 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2486 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2492 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2506 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsPeakRegion(v0))
2508 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2512 if(cont_flip[0]<2.3){
2521 if(cont_flip[0]<2.3){
2531 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2552 }
else if((trk->Charge()<0 && !anti_lambda_flag) || (trk->Charge()>0 && anti_lambda_flag)){
2576 if( melam_flip < 10. && cosoa < 0. && fAnalCuts->IsSideBand(v0))
2578 if((trk->Charge()>0 && !anti_lambda_flag) || (trk->Charge()<0 && anti_lambda_flag)){
2582 if(cont_flip[0]<2.3){
2591 if(cont_flip[0]<2.3){
2602 for(Int_t iv=0;iv<11;iv++){
2614 if(trk->Charge()>0){
2617 }
else if(trk->Charge()<0){
2626 Int_t pdgcode = mclc->GetPdgCode();
2627 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && abs(mcpdgv0_array[1])==4122){
2630 if(abs(pdgcode)==4122 && abs(mcpdgele_array[1])==4122 && ( abs(mcpdgv0_array[1])==3212 || abs(mcpdgv0_array[1])==3214 || abs(mcpdgv0_array[1])==3224) ){
2633 if(abs(pdgcode)==4132 && abs(mcpdgele_array[1])==4132 && abs(mcpdgv0_array[1])==3312){
2636 if(abs(pdgcode)==4232 && abs(mcpdgele_array[1])==4232 && abs(mcpdgv0_array[1])==3322){
2639 if(abs(pdgcode)==5122 && abs(mcpdgele_array[1])==5122 && abs(mcpdgv0_array[1])==4122 && abs(mcpdgv0_array[2])==5122){
2642 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5132){
2645 if(abs(pdgcode)==5132 && abs(mcpdgele_array[1])==5132 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5132){
2648 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4132 && abs(mcpdgv0_array[2])==5232){
2651 if(abs(pdgcode)==5232 && abs(mcpdgele_array[1])==5232 && abs(mcpdgv0_array[1])==4232 && abs(mcpdgv0_array[2])==5232){
2687 for(Int_t i=0;i<94;i++)
2693 Double_t pxe = trke->Px();
2694 Double_t pye = trke->Py();
2695 Double_t pze = trke->Pz();
2696 Double_t mome = sqrt(pxe*pxe+pye*pye+pze*pze);
2697 Double_t Ee = sqrt(mome*mome+0.000510998928*0.000510998928);
2699 Double_t pxv = v0->Px();
2700 Double_t pyv = v0->Py();
2701 Double_t pzv = v0->Pz();
2702 Double_t momv = sqrt(pxv*pxv+pyv*pyv+pzv*pzv);
2703 Double_t Ev = sqrt(momv*momv+1.115683*1.115683);
2705 Double_t cosoa = (pxe*pxv+pye*pyv+pze*pzv)/mome/momv;
2707 Double_t pxsum = pxe + pxv;
2708 Double_t pysum = pye + pyv;
2709 Double_t pzsum = pze + pzv;
2710 Double_t Esum = Ee + Ev;
2711 Double_t mel = sqrt(Esum*Esum-pxsum*pxsum-pysum*pysum-pzsum*pzsum);
2713 Double_t uxe = pxe/mome;
2714 Double_t uye = pye/mome;
2715 Double_t uze = pze/mome;
2716 Double_t lf = -2.*(pxv*uxe+pyv*uye+pzv*uze);
2717 Double_t pxv_flip = pxv + lf * uxe;
2718 Double_t pyv_flip = pyv + lf * uye;
2719 Double_t pzv_flip = pzv + lf * uze;
2720 Double_t pxsum_flip = pxe + pxv_flip;
2721 Double_t pysum_flip = pye + pyv_flip;
2722 Double_t pzsum_flip = pze + pzv_flip;
2723 Double_t mel_flip = sqrt(Esum*Esum-pxsum_flip*pxsum_flip-pysum_flip*pysum_flip-pzsum_flip*pzsum_flip);
2724 Double_t ptel_flip = sqrt(pxsum_flip*pxsum_flip+pysum_flip*pysum_flip);
2727 UInt_t pdgdg[2]={11,3122};
2780 cont[1] = sqrt(pxsum*pxsum+pysum*pysum);
2784 Double_t cont_flip[3];
2785 cont_flip[0] = mel_flip;
2786 cont_flip[1] = ptel_flip;
2791 cont2[1] = sqrt(pxe*pxe+pye*pye);
2793 Double_t cont_eleptvseta[3];
2794 cont_eleptvseta[0] = trke->Pt();
2795 cont_eleptvseta[1] = trke->Eta();
2798 Double_t cont_eleptvslambdapt[3];
2799 cont_eleptvslambdapt[0] = trke->Pt();
2800 cont_eleptvslambdapt[1] = v0->Pt();
2803 Double_t cont_eleptvsd0[3];
2804 cont_eleptvsd0[0] = trke->Pt();
2805 cont_eleptvsd0[1] = 0.;
2808 Double_t cont_eleptvsv0dl[3];
2809 cont_eleptvsv0dl[0] = trke->Pt();
2810 cont_eleptvsv0dl[1] = v0info[0];
2813 Double_t cont_eleptvsv0dca[3];
2814 cont_eleptvsv0dca[0] = trke->Pt();
2815 cont_eleptvsv0dca[1] = v0info[1];
2818 Double_t cont_elelamptvsv0dl[4];
2819 cont_elelamptvsv0dl[0] = sqrt(pxsum*pxsum+pysum*pysum);
2820 cont_elelamptvsv0dl[1] = v0info[0];
2821 cont_elelamptvsv0dl[2] = 0.0;
2824 Double_t cont_elelamptvsv0dl_flip[4];
2825 cont_elelamptvsv0dl_flip[0] = ptel_flip;
2826 cont_elelamptvsv0dl_flip[1] = v0info[0];
2827 cont_elelamptvsv0dl_flip[2] = 0.0;
2830 Double_t xyzR125_ele[3], xyzR125_pr[3], xyzR125_pi[3];
2831 xyzR125_ele[0] = (*elevars)[0];
2832 xyzR125_ele[1] = (*elevars)[1];
2833 xyzR125_ele[2] = (*elevars)[2];
2834 xyzR125_pr[0] = (*v0vars)[0];
2835 xyzR125_pr[1] = (*v0vars)[1];
2836 xyzR125_pr[2] = (*v0vars)[2];
2837 xyzR125_pi[0] = (*v0vars)[3];
2838 xyzR125_pi[1] = (*v0vars)[4];
2839 xyzR125_pi[2] = (*v0vars)[5];
2841 Double_t rdhfcutvars[9];
2842 rdhfcutvars[0] = xyzR125_ele[0];
2843 rdhfcutvars[1] = xyzR125_ele[1];
2844 rdhfcutvars[2] = xyzR125_ele[2];
2845 rdhfcutvars[3] = xyzR125_pr[0];
2846 rdhfcutvars[4] = xyzR125_pr[1];
2847 rdhfcutvars[5] = xyzR125_pr[2];
2848 rdhfcutvars[6] = xyzR125_pi[0];
2849 rdhfcutvars[7] = xyzR125_pi[1];
2850 rdhfcutvars[8] = xyzR125_pi[2];
2861 if(((
int)trke->T())*chargepr>0){
2921 if(mel_flip < 10. && cosoa<0. && fAnalCuts->IsPeakRegion(v0))
2923 if(((
int)trke->T())*chargepr>0){
2927 if(cont_flip[0]<2.3){
2936 if(cont_flip[0]<2.3){
2947 for(Int_t iv=0;iv<11;iv++){
2955 fCorrelationVariables[4] = v0info[0];
2956 fCorrelationVariables[5] = (*elevars)[5];
2957 fCorrelationVariables[6] = 2;
2959 if(chargepr>0) fCorrelationVariables[7] = 0;
2960 else fCorrelationVariables[7] = 2;
2961 }
else if(trke->T()<0){
2962 if(chargepr>0) fCorrelationVariables[7] = 1;
2963 else fCorrelationVariables[7] = 3;
2981 const char* nameoutput = GetOutputSlot(5)->GetContainer()->GetName();
2985 TString * fCandidateVariableNames =
new TString[nVar];
2987 fCandidateVariableNames[ 0]=
"ElePx";
2988 fCandidateVariableNames[ 1]=
"ElePy";
2989 fCandidateVariableNames[ 2]=
"ElePz";
2990 fCandidateVariableNames[ 3]=
"TPCChi2overNDF";
2991 fCandidateVariableNames[ 4]=
"ITSNcls";
2992 fCandidateVariableNames[ 5]=
"TPCNcls";
2993 fCandidateVariableNames[ 6]=
"TPCNclsPID";
2994 fCandidateVariableNames[ 7]=
"TPCNclsRatio";
2995 fCandidateVariableNames[ 8]=
"d0R";
2996 fCandidateVariableNames[ 9]=
"d0Z";
2997 fCandidateVariableNames[10]=
"ITSClusterMap";
2998 fCandidateVariableNames[11]=
"nSigmaTPCele";
2999 fCandidateVariableNames[12]=
"nSigmaTOFele";
3000 fCandidateVariableNames[13]=
"nSigmaTPCpi";
3001 fCandidateVariableNames[14]=
"nSigmaTPCka";
3002 fCandidateVariableNames[15]=
"nSigmaTPCpr";
3003 fCandidateVariableNames[16]=
"EvNumber";
3004 fCandidateVariableNames[17]=
"EleCharge";
3005 fCandidateVariableNames[18]=
"ElePdgCode";
3006 fCandidateVariableNames[19]=
"EleMotherPdgCode";
3007 fCandidateVariableNames[20]=
"mcelepx";
3008 fCandidateVariableNames[21]=
"mcelepy";
3009 fCandidateVariableNames[22]=
"mcelepz";
3010 fCandidateVariableNames[23]=
"Centrality";
3011 fCandidateVariableNames[24]=
"PrimVertZ";
3012 fCandidateVariableNames[25]=
"RunNumber";
3014 for (Int_t ivar=0; ivar<nVar; ivar++) {
3032 Double_t d0z0[2],covd0z0[3];
3033 trk->PropagateToDCA(
fVtx1,
fBzkG,kVeryBig,d0z0,covd0z0);
3037 pv[0] =
fVtx1->GetX();
3038 pv[1] =
fVtx1->GetY();
3039 pv[2] =
fVtx1->GetZ();
3040 Double_t xyzR125[3]={9999.,9999.,9999.};
3042 TVector *varvec =
new TVector(6);
3043 (*varvec)[0] = xyzR125[0];
3044 (*varvec)[1] = xyzR125[1];
3045 (*varvec)[2] = xyzR125[2];
3048 (*varvec)[5] = d0z0[0];
3050 fElectronTracks->AddLast(
new TLorentzVector(trk->Px(),trk->Py(),trk->Pz(),trk->Charge()));
3054 Int_t pdgEle = -9999;
3055 Int_t pdgEleMother = -9999;
3056 Float_t mcelepx = -9999;
3057 Float_t mcelepy = -9999;
3058 Float_t mcelepz = -9999;
3061 Int_t labEle = trk->GetLabel();
3062 if(labEle<0)
return;
3063 AliAODMCParticle *mcetrk = (AliAODMCParticle*)mcArray->At(labEle);
3065 pdgEle = mcetrk->GetPdgCode();
3066 mcelepx = mcetrk->Px();
3067 mcelepy = mcetrk->Py();
3068 mcelepz = mcetrk->Pz();
3085 for(Int_t i=0;i<26;i++){
3096 if(trk->GetTPCNclsF()>0)
3101 Int_t itsmap = trk->GetITSClusterMap();
3104 Bool_t spdfirst = (itsmap & bit1) == bit1;
3105 Bool_t spdsecond = (itsmap & bit2) == bit2;
3143 const char* nameoutput = GetOutputSlot(6)->GetContainer()->GetName();
3147 TString * fCandidateVariableNames =
new TString[nVar];
3149 fCandidateVariableNames[ 0]=
"V0Px";
3150 fCandidateVariableNames[ 1]=
"V0Py";
3151 fCandidateVariableNames[ 2]=
"V0Pz";
3152 fCandidateVariableNames[ 3]=
"MassLambda";
3153 fCandidateVariableNames[ 4]=
"MassAntiLambda";
3154 fCandidateVariableNames[ 5]=
"ProtonPx";
3155 fCandidateVariableNames[ 6]=
"ProtonPy";
3156 fCandidateVariableNames[ 7]=
"ProtonPz";
3157 fCandidateVariableNames[ 8]=
"PionPx";
3158 fCandidateVariableNames[ 9]=
"PionPy";
3159 fCandidateVariableNames[10]=
"PionPz";
3160 fCandidateVariableNames[11]=
"RfidV0";
3161 fCandidateVariableNames[12]=
"DcaV0PrToPrimVertex";
3162 fCandidateVariableNames[13]=
"DcaV0PiToPrimVertex";
3163 fCandidateVariableNames[14]=
"DcaV0daughters";
3164 fCandidateVariableNames[15]=
"V0CosPointingAngle";
3165 fCandidateVariableNames[16]=
"V0ProperDecayLength";
3166 fCandidateVariableNames[17]=
"MassK0Short";
3167 fCandidateVariableNames[18]=
"nSigmaTPCpr";
3168 fCandidateVariableNames[19]=
"nSigmaTPCpi";
3169 fCandidateVariableNames[20]=
"TPCNCrossV0Pr";
3170 fCandidateVariableNames[21]=
"TPCNCrossV0Pi";
3171 fCandidateVariableNames[22]=
"TPCNCrossRatioV0Pr";
3172 fCandidateVariableNames[23]=
"TPCNCrossRatioV0Pi";
3173 fCandidateVariableNames[24]=
"V0PdgCode";
3174 fCandidateVariableNames[25]=
"V0MotherPdgCode";
3175 fCandidateVariableNames[26]=
"mcv0px";
3176 fCandidateVariableNames[27]=
"mcv0py";
3177 fCandidateVariableNames[28]=
"mcv0pz";
3178 fCandidateVariableNames[29]=
"mcv0vertx";
3179 fCandidateVariableNames[30]=
"mcv0verty";
3180 fCandidateVariableNames[31]=
"mcv0vertz";
3181 fCandidateVariableNames[32]=
"EvNumber";
3182 fCandidateVariableNames[33]=
"Centrality";
3183 fCandidateVariableNames[34]=
"PrimVertZ";
3184 fCandidateVariableNames[35]=
"RunNumber";
3186 for (Int_t ivar=0; ivar<nVar; ivar++) {
3200 AliAODTrack *cptrack = (AliAODTrack*)(v0->GetDaughter(0));
3201 AliAODTrack *cntrack = (AliAODTrack*)(v0->GetDaughter(1));
3202 if(!cptrack)
return;
3203 if(!cntrack)
return;
3222 Double_t mlamPDG = TDatabasePDG::Instance()->GetParticle(3122)->Mass();
3223 Double_t posVtx[3] = {0.,0.,0.};
3224 fVtx1->GetXYZ(posVtx);
3225 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3226 Double_t v0propdl = v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam;
3227 Double_t cont_arm[3];
3228 cont_arm[0] = v0->Pt();
3229 cont_arm[1] = v0->PtArmV0();
3230 cont_arm[2] = v0->AlphaV0();
3249 Double_t momv0x = v0->MomV0X();
3250 Double_t momv0y = v0->MomV0Y();
3251 Double_t phi_alice = atan2(momv0y,momv0x);
3252 if(phi_alice<0.) phi_alice += 2 * M_PI;
3255 Int_t v0pdgcode = -9999;
3256 Int_t v0motherpdgcode = -9999;
3257 Float_t mcv0px = -9999;
3258 Float_t mcv0py = -9999;
3259 Float_t mcv0pz = -9999;
3260 Float_t mcv0vertx = -9999;
3261 Float_t mcv0verty = -9999;
3262 Float_t mcv0vertz = -9999;
3265 Int_t pdgdgv0[2]={2212,211};
3266 Int_t labV0 = v0->MatchToMC(3122,mcArray,2,pdgdgv0);
3280 AliAODMCParticle *mcv0trk = (AliAODMCParticle*)mcArray->At(labV0);
3281 if(!mcv0trk)
return;
3283 Bool_t hfv0_flag = kFALSE;
3284 Bool_t promptv0_flag = kFALSE;
3285 v0pdgcode = mcv0trk->GetPdgCode();
3286 Int_t labv0mother = mcv0trk->GetMother();
3288 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3290 v0motherpdgcode = motherv0->GetPdgCode();
3291 if(abs(v0motherpdgcode)>4000&&abs(v0motherpdgcode)<4400){
3294 if(abs(v0motherpdgcode)==3322){
3296 }
else if(abs(v0motherpdgcode)==3312){
3298 }
else if(abs(v0motherpdgcode)==3334){
3302 promptv0_flag = kTRUE;
3306 if(!hfv0_flag)
return;
3309 mcv0px = mcv0trk->Px();
3310 mcv0py = mcv0trk->Py();
3311 mcv0pz = mcv0trk->Pz();
3313 Int_t labptrk = cptrack->GetLabel();
3315 AliAODMCParticle *mcv0ptrk = (AliAODMCParticle*)mcArray->At(labptrk);
3317 mcv0vertx = mcv0ptrk->Xv();
3318 mcv0verty = mcv0ptrk->Yv();
3319 mcv0vertz = mcv0ptrk->Zv();
3326 Double_t posVtx[3] = {0.,0.,0.};
3327 fVtx1->GetXYZ(posVtx);
3328 TLorentzVector *lv =
new TLorentzVector();
3329 Double_t ptotlam = TMath::Sqrt(pow(v0->Px(),2)+pow(v0->Py(),2)+pow(v0->Pz(),2));
3330 Double_t xyzR125pr[3] = {9999.,9999.,9999.};
3331 Double_t xyzR125pi[3] = {9999.,9999.,9999.};
3333 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassLambda());
3335 fV0dlArray1.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3339 TVector *varvec =
new TVector(15);
3340 (*varvec)[0] = xyzR125pr[0];
3341 (*varvec)[1] = xyzR125pr[1];
3342 (*varvec)[2] = xyzR125pr[2];
3343 (*varvec)[3] = xyzR125pi[0];
3344 (*varvec)[4] = xyzR125pi[1];
3345 (*varvec)[5] = xyzR125pi[2];
3346 (*varvec)[6] = v0->DecayVertexV0X();
3347 (*varvec)[7] = v0->DecayVertexV0Y();
3348 (*varvec)[8] = v0->DecayVertexV0Z();
3349 (*varvec)[9] = v0->MomPosX();
3350 (*varvec)[10] = v0->MomPosY();
3351 (*varvec)[11] = v0->MomPosZ();
3352 (*varvec)[12] = v0->MomNegX();
3353 (*varvec)[13] = v0->MomNegY();
3354 (*varvec)[14] = v0->MomNegZ();
3357 lv->SetXYZM(v0->Px(),v0->Py(),v0->Pz(),v0->MassAntiLambda());
3359 fV0dlArray2.push_back(v0->DecayLengthV0(posVtx)*mlamPDG/ptotlam);
3363 TVector *varvec =
new TVector(15);
3364 (*varvec)[0] = xyzR125pr[0];
3365 (*varvec)[1] = xyzR125pr[1];
3366 (*varvec)[2] = xyzR125pr[2];
3367 (*varvec)[3] = xyzR125pi[0];
3368 (*varvec)[4] = xyzR125pi[1];
3369 (*varvec)[5] = xyzR125pi[2];
3370 (*varvec)[6] = v0->DecayVertexV0X();
3371 (*varvec)[7] = v0->DecayVertexV0Y();
3372 (*varvec)[8] = v0->DecayVertexV0Z();
3373 (*varvec)[9] = v0->MomNegX();
3374 (*varvec)[10] = v0->MomNegY();
3375 (*varvec)[11] = v0->MomNegZ();
3376 (*varvec)[12] = v0->MomPosX();
3377 (*varvec)[13] = v0->MomPosY();
3378 (*varvec)[14] = v0->MomPosZ();
3385 for(Int_t i=0;i<36;i++){
3395 Bool_t isparticle = kTRUE;
3415 lPosV0[0] = v0->DecayVertexV0X();
3416 lPosV0[1] = v0->DecayVertexV0Y();
3417 lPosV0[2] = v0->DecayVertexV0Z();
3418 Double_t decayvertV0 = TMath::Sqrt(lPosV0[0]*lPosV0[0]+lPosV0[1]*lPosV0[1]);
3421 Double_t lDcaPosToPrimVertex = v0->DcaPosToPrimVertex();
3422 Double_t lDcaNegToPrimVertex = v0->DcaNegToPrimVertex();
3452 if(cptrack->GetTPCNclsF()>0)
3453 fCandidateV0Variables[22] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3454 if(cntrack->GetTPCNclsF()>0)
3455 fCandidateV0Variables[23] =(Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3459 if(cntrack->GetTPCNclsF()>0)
3460 fCandidateV0Variables[22] = (Float_t) cntrack->GetTPCClusterInfo(2,1)/(Float_t)cntrack->GetTPCNclsF();
3461 if(cptrack->GetTPCNclsF()>0)
3462 fCandidateV0Variables[23] = (Float_t) cptrack->GetTPCClusterInfo(2,1)/(Float_t)cptrack->GetTPCNclsF();
3487 const char* nameoutput = GetOutputSlot(7)->GetContainer()->GetName();
3491 TString * fCandidateVariableNames =
new TString[nVar];
3493 fCandidateVariableNames[ 0]=
"Centrality";
3494 fCandidateVariableNames[ 1]=
"DecayType";
3495 fCandidateVariableNames[ 2]=
"LcPx";
3496 fCandidateVariableNames[ 3]=
"LcPy";
3497 fCandidateVariableNames[ 4]=
"LcPz";
3498 fCandidateVariableNames[ 5]=
"ElePx";
3499 fCandidateVariableNames[ 6]=
"ElePy";
3500 fCandidateVariableNames[ 7]=
"ElePz";
3501 fCandidateVariableNames[ 8]=
"V0Px";
3502 fCandidateVariableNames[ 9]=
"V0Py";
3503 fCandidateVariableNames[10]=
"V0Pz";
3505 for (Int_t ivar=0; ivar<nVar; ivar++) {
3517 if(!mcepart)
return;
3518 if(!mcv0part)
return;
3520 for(Int_t i=0;i<11;i++){
3536 Double_t epx = mcepart->Px();
3537 Double_t epy = mcepart->Py();
3538 Double_t epz = mcepart->Pz();
3539 Double_t eE = sqrt(epx*epx+epy*epy+epz*epz+0.000511*0.000511);
3540 Double_t v0px = mcv0part->Px();
3541 Double_t v0py = mcv0part->Py();
3542 Double_t v0pz = mcv0part->Pz();
3543 Double_t v0E = sqrt(v0px*v0px+v0py*v0py+v0pz*v0pz+1.1157*1.1157);
3545 Double_t InvMassEleLambda = sqrt(pow(eE+v0E,2)-pow(epx+v0px,2)-pow(epy+v0py,2)-pow(epz+v0pz,2));
3548 cont[0] = InvMassEleLambda;
3549 cont[1] = mcpart->Pt();
3552 cont2[0] = InvMassEleLambda;
3553 cont2[1] = mcepart->Pt();
3555 Double_t cont_eleptvseta[3];
3556 cont_eleptvseta[0] = mcepart->Pt();
3557 cont_eleptvseta[1] = mcepart->Eta();
3559 Double_t cont_eleptvslambdapt[3];
3560 cont_eleptvslambdapt[0] = mcepart->Pt();
3561 cont_eleptvslambdapt[1] = mcv0part->Pt();
3563 Double_t cont_eleptvslambdaptvslcpt[4];
3564 cont_eleptvslambdaptvslcpt[0] = mcepart->Pt();
3565 cont_eleptvslambdaptvslcpt[1] = mcv0part->Pt();
3566 cont_eleptvslambdaptvslcpt[2] = mcpart->Pt();
3569 contmc[0] = mcpart->Pt();
3570 contmc[1] = mcpart->Y();
3572 Double_t contmcele[3];
3573 contmcele[0] = mcepart->Pt();
3574 contmcele[1] = mcepart->Eta();
3579 esdcuts->GetEtaRange(etamin,etamax);
3589 if(fabs(mcepart->Eta())<etamax){
3593 if(InvMassEleLambda<2.3){
3599 if(fabs(mcpart->Y())<0.7){
3600 if(InvMassEleLambda<2.3){
3605 }
else if(decaytype==1){
3613 if(fabs(mcepart->Eta())<etamax){
3617 if(InvMassEleLambda<2.3){
3623 }
else if(decaytype==2){
3631 if(fabs(mcepart->Eta())<etamax){
3635 if(InvMassEleLambda<2.3){
3641 }
else if(decaytype==3){
3657 const char* nameoutput = GetOutputSlot(9)->GetContainer()->GetName();
3661 TString * fCandidateVariableNames =
new TString[nVar];
3663 fCandidateVariableNames[ 0]=
"Centrality";
3664 fCandidateVariableNames[ 1]=
"ElePx";
3665 fCandidateVariableNames[ 2]=
"ElePy";
3666 fCandidateVariableNames[ 3]=
"ElePz";
3667 fCandidateVariableNames[ 4]=
"ElePdgCode";
3668 fCandidateVariableNames[ 5]=
"EleMotherPdgCode";
3669 fCandidateVariableNames[ 6]=
"RunNumber";
3670 fCandidateVariableNames[ 7]=
"EvNumber";
3672 for (Int_t ivar=0; ivar<nVar; ivar++) {
3683 if(!mcepart)
return;
3685 Int_t pdgarray_ele[100], labelarray_ele[100], ngen_ele;
3691 contmc[0] = mcepart->Pt();
3692 contmc[1] = mcepart->Eta();
3695 if(ele_from_bottom){
3697 }
else if(ele_from_charm){
3701 Bool_t hfe_flag = kFALSE;
3702 Int_t labemother = mcepart->GetMother();
3703 Int_t pdgmotherele = -9999;
3705 AliAODMCParticle *motherele = (AliAODMCParticle*)mcArray->At(labemother);
3706 pdgmotherele = motherele->GetPdgCode();
3707 if(abs(pdgmotherele)>4000&&abs(pdgmotherele)<4400){
3712 if(!hfe_flag)
return;
3716 for(Int_t i=0;i<8;i++){
3740 const char* nameoutput = GetOutputSlot(10)->GetContainer()->GetName();
3744 TString * fCandidateVariableNames =
new TString[nVar];
3746 fCandidateVariableNames[ 0]=
"Centrality";
3747 fCandidateVariableNames[ 1]=
"V0Px";
3748 fCandidateVariableNames[ 2]=
"V0Py";
3749 fCandidateVariableNames[ 3]=
"V0Pz";
3750 fCandidateVariableNames[ 4]=
"V0PdgCode";
3751 fCandidateVariableNames[ 5]=
"V0MotherPdgCode";
3752 fCandidateVariableNames[ 6]=
"RunNumber";
3753 fCandidateVariableNames[ 7]=
"EvNumber";
3755 for (Int_t ivar=0; ivar<nVar; ivar++) {
3766 if(!mcv0part)
return;
3768 for(Int_t i=0;i<8;i++){
3772 Bool_t hfv0_flag = kFALSE;
3773 Int_t labv0mother = mcv0part->GetMother();
3774 Int_t pdgmotherv0 = -9999;
3776 AliAODMCParticle *motherv0 = (AliAODMCParticle*)mcArray->At(labv0mother);
3778 pdgmotherv0 = motherv0->GetPdgCode();
3779 if(abs(pdgmotherv0)>4000&&abs(pdgmotherv0)<4400){
3784 if(!hfv0_flag)
return;
3787 contmc[0] = mcv0part->Pt();
3788 contmc[1] = mcv0part->Eta();
3812 fCEvents =
new TH1F(
"fCEvents",
"conter",18,-0.5,17.5);
3814 fCEvents->GetXaxis()->SetBinLabel(1,
"X1");
3815 fCEvents->GetXaxis()->SetBinLabel(2,
"Analyzed events");
3816 fCEvents->GetXaxis()->SetBinLabel(3,
"AliAODVertex exists");
3817 fCEvents->GetXaxis()->SetBinLabel(4,
"TriggerOK");
3818 fCEvents->GetXaxis()->SetBinLabel(5,
"IsEventSelected");
3819 fCEvents->GetXaxis()->SetBinLabel(6,
"CascadesHF exists");
3820 fCEvents->GetXaxis()->SetBinLabel(7,
"MCarray exists");
3821 fCEvents->GetXaxis()->SetBinLabel(8,
"MCheader exists");
3822 fCEvents->GetXaxis()->SetBinLabel(9,
"triggerClass!=CINT1");
3823 fCEvents->GetXaxis()->SetBinLabel(10,
"triggerMask!=kAnyINT");
3824 fCEvents->GetXaxis()->SetBinLabel(11,
"triggerMask!=kAny");
3825 fCEvents->GetXaxis()->SetBinLabel(12,
"vtxTitle.Contains(Z)");
3826 fCEvents->GetXaxis()->SetBinLabel(13,
"vtxTitle.Contains(3D)");
3827 fCEvents->GetXaxis()->SetBinLabel(14,
"vtxTitle.Doesn'tContain(Z-3D)");
3829 fCEvents->GetXaxis()->SetBinLabel(16,
"!IsEventSelected");
3830 fCEvents->GetXaxis()->SetBinLabel(17,
"triggerMask!=kAnyINT || triggerClass!=CINT1");
3833 fCEvents->GetYaxis()->SetTitle(
"counts");
3835 fHTrigger =
new TH1F(
"fHTrigger",
"counter",18,-0.5,17.5);
3837 fHTrigger->GetXaxis()->SetBinLabel(1,
"X1");
3838 fHTrigger->GetXaxis()->SetBinLabel(2,
"kMB");
3839 fHTrigger->GetXaxis()->SetBinLabel(3,
"kSemiCentral");
3840 fHTrigger->GetXaxis()->SetBinLabel(4,
"kCentral");
3841 fHTrigger->GetXaxis()->SetBinLabel(5,
"kINT7");
3842 fHTrigger->GetXaxis()->SetBinLabel(6,
"kEMC7");
3844 fHTrigger->GetXaxis()->SetBinLabel(8,
"kMB|kSemiCentral|kCentral");
3845 fHTrigger->GetXaxis()->SetBinLabel(9,
"kINT7|kEMC7");
3846 fHTrigger->GetXaxis()->SetBinLabel(11,
"kMB&kSemiCentral");
3847 fHTrigger->GetXaxis()->SetBinLabel(12,
"kMB&kCentral");
3848 fHTrigger->GetXaxis()->SetBinLabel(13,
"kINT7&kEMC7");
3850 fHCentrality =
new TH1F(
"fHCentrality",
"conter",100,0.,100.);
3869 Int_t bins_base[3]= {23,100 ,10};
3870 Double_t xmin_base[3]={1.1,0 ,0.00};
3871 Double_t xmax_base[3]={5.7,20. ,100};
3872 fHistoEleLambdaMass =
new THnSparseF(
"fHistoEleLambdaMass",
"",3,bins_base,xmin_base,xmax_base);
3948 Int_t bins_base_elept[3]= {10,100 ,10};
3949 Double_t xmin_base_elept[3]={1.1,0 ,0.00};
3950 Double_t xmax_base_elept[3]={3.1,10. ,100};
4009 fHistoElePtRS =
new TH2F(
"fHistoElePtRS",
"",100,0,10,10,0,100);
4011 fHistoElePtWS =
new TH2F(
"fHistoElePtWS",
"",100,0.,10.,10,0,100);
4017 fHistoElePtMCS =
new TH2F(
"fHistoElePtMCS",
"",100,0,10,10,0,100);
4022 Int_t bins_eleptvseta[3]= {10,10 ,10};
4023 Double_t xmin_eleptvseta[3]={0.,-1. ,0.0};
4024 Double_t xmax_eleptvseta[3]={5.,1. ,100};
4026 fHistoElePtvsEtaRS =
new THnSparseF(
"fHistoElePtvsEtaRS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4028 fHistoElePtvsEtaWS =
new THnSparseF(
"fHistoElePtvsEtaWS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4030 fHistoElePtvsEtaRSMix =
new THnSparseF(
"fHistoElePtvsEtaRSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4032 fHistoElePtvsEtaWSMix =
new THnSparseF(
"fHistoElePtvsEtaWSMix",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4034 fHistoElePtvsEtaMCS =
new THnSparseF(
"fHistoElePtvsEtaMCS",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4036 fHistoElePtvsEtaMCGen =
new THnSparseF(
"fHistoElePtvsEtaMCGen",
"",3,bins_eleptvseta,xmin_eleptvseta,xmax_eleptvseta);
4039 Int_t bins_eleptvslambdapt[3]= {10,10 ,10};
4040 Double_t xmin_eleptvslambdapt[3]={0.,0. ,0.0};
4041 Double_t xmax_eleptvslambdapt[3]={5.,5. ,100};
4043 fHistoElePtvsLambdaPtRS =
new THnSparseF(
"fHistoElePtvsLambdaPtRS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4045 fHistoElePtvsLambdaPtWS =
new THnSparseF(
"fHistoElePtvsLambdaPtWS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4047 fHistoElePtvsLambdaPtRSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtRSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4049 fHistoElePtvsLambdaPtWSMix =
new THnSparseF(
"fHistoElePtvsLambdaPtWSMix",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4051 fHistoElePtvsLambdaPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtMCS",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4053 fHistoElePtvsLambdaPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4055 fHistoElePtvsLambdaPtMCLcGen =
new THnSparseF(
"fHistoElePtvsLambdaPtMCLcGen",
"",3,bins_eleptvslambdapt,xmin_eleptvslambdapt,xmax_eleptvslambdapt);
4058 Int_t bins_eleptvslambdaptvslcpt[4]= {10,10,10,10};
4059 Double_t xmin_eleptvslambdaptvslcpt[4]={0.,0.,0.,0.0};
4060 Double_t xmax_eleptvslambdaptvslcpt[4]={5.,5.,10.,100};
4061 fHistoElePtvsLambdaPtvsLcPtMCS =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCS",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4063 fHistoElePtvsLambdaPtvsLcPtMCGen =
new THnSparseF(
"fHistoElePtvsLambdaPtvsLcPtMCGen",
"",4,bins_eleptvslambdaptvslcpt,xmin_eleptvslambdaptvslcpt,xmax_eleptvslambdaptvslcpt);
4066 Int_t bins_allpt[4]= {10,10,20,20};
4067 Double_t xmin_allpt[4]={0.,0.,0.,0.0};
4068 Double_t xmax_allpt[4]={20.,20.,10.,10};
4079 Int_t bins_eleptvsd0[3]= {10 ,10 ,10};
4080 Double_t xmin_eleptvsd0[3]={0.,-0.2 ,0.0};
4081 Double_t xmax_eleptvsd0[3]={5.,0.2 ,100};
4083 fHistoElePtvsd0RS =
new THnSparseF(
"fHistoElePtvsd0RS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4085 fHistoElePtvsd0WS =
new THnSparseF(
"fHistoElePtvsd0WS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4087 fHistoElePtvsd0RSMix =
new THnSparseF(
"fHistoElePtvsd0RSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4089 fHistoElePtvsd0WSMix =
new THnSparseF(
"fHistoElePtvsd0WSMix",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4091 fHistoElePtvsd0MCS =
new THnSparseF(
"fHistoElePtvsd0MCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4093 fHistoElePtvsd0PromptMCS =
new THnSparseF(
"fHistoElePtvsd0PromptMCS",
"",3,bins_eleptvsd0,xmin_eleptvsd0,xmax_eleptvsd0);
4159 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4165 fHistoBachPt =
new TH1F(
"fHistoBachPt",
"Bachelor p_{T}",100,0.0,5.0);
4167 fHistoBachPtMCS =
new TH1F(
"fHistoBachPtMCS",
"Bachelor p_{T}",100,0.0,5.0);
4171 fHistod0Bach =
new TH1F(
"fHistod0Bach",
"Bachelor d_{0}",100,-0.5,0.5);
4173 fHistoLambdaMassvsPt=
new TH2F(
"fHistoLambdaMassvsPt",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4175 fHistoLambdaMassvsPtMCS=
new TH2F(
"fHistoLambdaMassvsPtMCS",
"Lambda mass",100,1.116-0.05,1.116+0.05,20,0.,10.);
4185 fHistoK0sMassvsPt=
new TH2F(
"fHistoK0sMassvsPt",
"K0s mass",100,0.497-0.05,0.497+0.05,20,0.,10.);
4187 fHistoLambdaPtvsDl=
new TH2F(
"fHistoLambdaPtvsDl",
"Lambda pt vs dl",20,0.,10.,20,0.,40.);
4200 Int_t bins_armenteros[3]= {20 ,100 ,200};
4201 Double_t xmin_armenteros[3]={0.,0. ,-1.};
4202 Double_t xmax_armenteros[3]={20.,1. ,1.};
4203 fHistoLambdaArmenteros =
new THnSparseF(
"fHistoLambdaArmenteros",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4205 fHistoLambdaArmenterosMCS =
new THnSparseF(
"fHistoLambdaArmenterosMCS",
"",3,bins_armenteros,xmin_armenteros,xmax_armenteros);
4227 for(Int_t i=0;i<8;i++){
4236 Int_t bins_lcmcgen[3]= {100 ,20 ,10};
4237 Double_t xmin_lcmcgen[3]={0.,-1.0 ,0.0};
4238 Double_t xmax_lcmcgen[3]={20.,1.0 ,100};
4239 fHistoLcMCGen =
new THnSparseF(
"fHistoLcMCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4241 fHistoLcMCGen1 =
new THnSparseF(
"fHistoLcMCGen1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4243 fHistoLcMCGen2 =
new THnSparseF(
"fHistoLcMCGen2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4245 fHistoFeedDownXic0MCGen =
new THnSparseF(
"fHistoFeedDownXic0MCGen",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4264 fHistoLcMCS =
new THnSparseF(
"fHistoLcMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4266 fHistoLcMCS1 =
new THnSparseF(
"fHistoLcMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4268 fHistoLcMCS2 =
new THnSparseF(
"fHistoLcMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4270 fHistoFeedDownXic0MCS =
new THnSparseF(
"fHistoFeedDownXic0MCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4272 fHistoFeedDownXic0MCS1 =
new THnSparseF(
"fHistoFeedDownXic0MCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4274 fHistoFeedDownXic0MCS2 =
new THnSparseF(
"fHistoFeedDownXic0MCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4282 fHistoFeedDownSigmaMCS =
new THnSparseF(
"fHistoFeedDownSigmaMCS",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4284 fHistoFeedDownSigmaMCS1 =
new THnSparseF(
"fHistoFeedDownSigmaMCS1",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4286 fHistoFeedDownSigmaMCS2 =
new THnSparseF(
"fHistoFeedDownSigmaMCS2",
"",3,bins_lcmcgen,xmin_lcmcgen,xmax_lcmcgen);
4289 Int_t bins_elemcgen[3]= {100 ,20 ,10};
4290 Double_t xmin_elemcgen[3]={0.,-1.0 ,0.0};
4291 Double_t xmax_elemcgen[3]={10.,1.0 ,100};
4292 fHistoElectronMCGen =
new THnSparseF(
"fHistoElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4296 fHistoCharmElectronMCGen =
new THnSparseF(
"fHistoCharmElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4298 fHistoLcElectronMCGen =
new THnSparseF(
"fHistoLcElectronMCGen",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4300 fHistoLcElectronMCGen1 =
new THnSparseF(
"fHistoLcElectronMCGen1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4302 fHistoLcElectronMCGen2 =
new THnSparseF(
"fHistoLcElectronMCGen2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4316 fHistoLcElectronMCS =
new THnSparseF(
"fHistoLcElectronMCS",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4318 fHistoLcElectronMCS1 =
new THnSparseF(
"fHistoLcElectronMCS1",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4320 fHistoLcElectronMCS2 =
new THnSparseF(
"fHistoLcElectronMCS2",
"",3,bins_elemcgen,xmin_elemcgen,xmax_elemcgen);
4335 Int_t bins_lambdamcgen[3]= {50 ,20 ,10};
4336 Double_t xmin_lambdamcgen[3]={0.,-1.0 ,0.0};
4337 Double_t xmax_lambdamcgen[3]={10.,1.0 ,100};
4338 fHistoLambdaMCGen =
new THnSparseF(
"fHistoLambdaMCGen",
"",3,bins_lambdamcgen,xmin_lambdamcgen,xmax_lambdamcgen);
4341 Int_t bins_eleptvsv0dl[3]= {100 ,20 ,10};
4342 Double_t xmin_eleptvsv0dl[3]={0.,0. ,0.0};
4343 Double_t xmax_eleptvsv0dl[3]={10.,40. ,100};
4344 fHistoElePtvsV0dlRS =
new THnSparseF(
"fHistoElePtvsV0dlRS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4346 fHistoElePtvsV0dlRS1 =
new THnSparseF(
"fHistoElePtvsV0dlRS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4348 fHistoElePtvsV0dlRS2 =
new THnSparseF(
"fHistoElePtvsV0dlRS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4350 fHistoElePtvsV0dlRSSide =
new THnSparseF(
"fHistoElePtvsV0dlRSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4352 fHistoElePtvsV0dlRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4354 fHistoElePtvsV0dlRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlRSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4356 fHistoElePtvsV0dlRSMix =
new THnSparseF(
"fHistoElePtvsV0dlRSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4358 fHistoElePtvsV0dlRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4360 fHistoElePtvsV0dlRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlRSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4362 fHistoElePtvsV0dlWS =
new THnSparseF(
"fHistoElePtvsV0dlWS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4364 fHistoElePtvsV0dlWS1 =
new THnSparseF(
"fHistoElePtvsV0dlWS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4366 fHistoElePtvsV0dlWS2 =
new THnSparseF(
"fHistoElePtvsV0dlWS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4368 fHistoElePtvsV0dlWSSide =
new THnSparseF(
"fHistoElePtvsV0dlWSSide",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4370 fHistoElePtvsV0dlWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4372 fHistoElePtvsV0dlWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dlWSSide2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4374 fHistoElePtvsV0dlWSMix =
new THnSparseF(
"fHistoElePtvsV0dlWSMix",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4376 fHistoElePtvsV0dlWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4378 fHistoElePtvsV0dlWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dlWSMix2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4380 fHistoElePtvsV0dlMCS =
new THnSparseF(
"fHistoElePtvsV0dlMCS",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4382 fHistoElePtvsV0dlMCS1 =
new THnSparseF(
"fHistoElePtvsV0dlMCS1",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4384 fHistoElePtvsV0dlMCS2 =
new THnSparseF(
"fHistoElePtvsV0dlMCS2",
"",3,bins_eleptvsv0dl,xmin_eleptvsv0dl,xmax_eleptvsv0dl);
4399 Int_t bins_eleptvsv0dca[3]= {100 ,20 ,10};
4400 Double_t xmin_eleptvsv0dca[3]={0.,0. ,0.0};
4401 Double_t xmax_eleptvsv0dca[3]={10.,1. ,100};
4402 fHistoElePtvsV0dcaRS =
new THnSparseF(
"fHistoElePtvsV0dcaRS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4404 fHistoElePtvsV0dcaRS1 =
new THnSparseF(
"fHistoElePtvsV0dcaRS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4406 fHistoElePtvsV0dcaRS2 =
new THnSparseF(
"fHistoElePtvsV0dcaRS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4408 fHistoElePtvsV0dcaRSSide =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4410 fHistoElePtvsV0dcaRSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4412 fHistoElePtvsV0dcaRSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4414 fHistoElePtvsV0dcaRSMix =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4416 fHistoElePtvsV0dcaRSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4418 fHistoElePtvsV0dcaRSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaRSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4420 fHistoElePtvsV0dcaWS =
new THnSparseF(
"fHistoElePtvsV0dcaWS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4422 fHistoElePtvsV0dcaWS1 =
new THnSparseF(
"fHistoElePtvsV0dcaWS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4424 fHistoElePtvsV0dcaWS2 =
new THnSparseF(
"fHistoElePtvsV0dcaWS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4426 fHistoElePtvsV0dcaWSSide =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4428 fHistoElePtvsV0dcaWSSide1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4430 fHistoElePtvsV0dcaWSSide2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSSide2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4432 fHistoElePtvsV0dcaWSMix =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4434 fHistoElePtvsV0dcaWSMix1 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4436 fHistoElePtvsV0dcaWSMix2 =
new THnSparseF(
"fHistoElePtvsV0dcaWSMix2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4438 fHistoElePtvsV0dcaMCS =
new THnSparseF(
"fHistoElePtvsV0dcaMCS",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4440 fHistoElePtvsV0dcaMCS1 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS1",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4442 fHistoElePtvsV0dcaMCS2 =
new THnSparseF(
"fHistoElePtvsV0dcaMCS2",
"",3,bins_eleptvsv0dca,xmin_eleptvsv0dca,xmax_eleptvsv0dca);
4457 Int_t bins_elelamptvsv0dl[4]= {100 ,20 ,100, 10};
4458 Double_t xmin_elelamptvsv0dl[4]={0.,0. ,-0.5, 0.0};
4459 Double_t xmax_elelamptvsv0dl[4]={20.,40.,0.5 ,100};
4460 fHistoEleLambdaPtvsV0dlRS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4462 fHistoEleLambdaPtvsV0dlRS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4464 fHistoEleLambdaPtvsV0dlRS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4466 fHistoEleLambdaPtvsV0dlRSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4468 fHistoEleLambdaPtvsV0dlRSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4470 fHistoEleLambdaPtvsV0dlRSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4472 fHistoEleLambdaPtvsV0dlRSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4474 fHistoEleLambdaPtvsV0dlRSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4476 fHistoEleLambdaPtvsV0dlRSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4478 fHistoEleLambdaPtvsV0dlWS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4480 fHistoEleLambdaPtvsV0dlWS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4482 fHistoEleLambdaPtvsV0dlWS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4484 fHistoEleLambdaPtvsV0dlWSSide =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4486 fHistoEleLambdaPtvsV0dlWSSide1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4488 fHistoEleLambdaPtvsV0dlWSSide2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSSide2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4490 fHistoEleLambdaPtvsV0dlWSMix =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4492 fHistoEleLambdaPtvsV0dlWSMix1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4494 fHistoEleLambdaPtvsV0dlWSMix2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSMix2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4496 fHistoEleLambdaPtvsV0dlMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4498 fHistoEleLambdaPtvsV0dlMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4500 fHistoEleLambdaPtvsV0dlMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4576 fHistoEleLambdaPtvsV0dlXibMCS =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4578 fHistoEleLambdaPtvsV0dlXibMCS1 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS1",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4580 fHistoEleLambdaPtvsV0dlXibMCS2 =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlXibMCS2",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4583 fHistoEleLambdaPtvsV0dlRSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4585 fHistoEleLambdaPtvsV0dlRS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4587 fHistoEleLambdaPtvsV0dlRS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlRS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4601 fHistoEleLambdaPtvsV0dlWSAway =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWSAway",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4603 fHistoEleLambdaPtvsV0dlWS1Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS1Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4605 fHistoEleLambdaPtvsV0dlWS2Away =
new THnSparseF(
"fHistoEleLambdaPtvsV0dlWS2Away",
"",4,bins_elelamptvsv0dl,xmin_elelamptvsv0dl,xmax_elelamptvsv0dl);
4687 for(Int_t ih=0;ih<17;ih++){
4688 Int_t bins_eleptvscutvars[3];
4689 Double_t xmin_eleptvscutvars[3];
4690 Double_t xmax_eleptvscutvars[3];
4692 bins_eleptvscutvars[0] = 20;
4693 xmin_eleptvscutvars[0] = 0.;
4694 xmax_eleptvscutvars[0] = 20.;
4695 bins_eleptvscutvars[2] = 10;
4696 xmin_eleptvscutvars[2] = 0.;
4697 xmax_eleptvscutvars[2] = 100.;
4701 bins_eleptvscutvars[1] = 40;
4702 xmin_eleptvscutvars[1] = 0.;
4703 xmax_eleptvscutvars[1] = 160.;
4704 }
else if(ih==2 || ih==3){
4706 bins_eleptvscutvars[1] = 20;
4707 xmin_eleptvscutvars[1] = -5.;
4708 xmax_eleptvscutvars[1] = 5.;
4711 bins_eleptvscutvars[1] = 30;
4712 xmin_eleptvscutvars[1] = -1.5;
4713 xmax_eleptvscutvars[1] = 1.5;
4716 bins_eleptvscutvars[1] = 7;
4717 xmin_eleptvscutvars[1] = -0.5;
4718 xmax_eleptvscutvars[1] = 6.5;
4721 bins_eleptvscutvars[1] = 50;
4722 xmin_eleptvscutvars[1] = 1.1156-0.03;
4723 xmax_eleptvscutvars[1] = 1.1156+0.03;
4726 bins_eleptvscutvars[1] = 20;
4727 xmin_eleptvscutvars[1] = 0.;
4728 xmax_eleptvscutvars[1] = 5.;
4731 bins_eleptvscutvars[1] = 20;
4732 xmin_eleptvscutvars[1] = 0.;
4733 xmax_eleptvscutvars[1] = 2.;
4734 }
else if(ih==9 || ih==10 ){
4736 bins_eleptvscutvars[1] = 20;
4737 xmin_eleptvscutvars[1] = 0.;
4738 xmax_eleptvscutvars[1] = 0.5;
4741 bins_eleptvscutvars[1] = 20;
4742 xmin_eleptvscutvars[1] = 0.95;
4743 xmax_eleptvscutvars[1] = 1.0;
4746 bins_eleptvscutvars[1] = 50;
4747 xmin_eleptvscutvars[1] = 0.497-0.03;
4748 xmax_eleptvscutvars[1] = 0.497+0.03;
4749 }
else if(ih==13 || ih==14){
4751 bins_eleptvscutvars[1] = 20;
4752 xmin_eleptvscutvars[1] = -5;
4753 xmax_eleptvscutvars[1] = 5;
4756 bins_eleptvscutvars[1] = 30;
4757 xmin_eleptvscutvars[1] = -1.5;
4758 xmax_eleptvscutvars[1] = 1.5;
4761 bins_eleptvscutvars[1] = 20;
4762 xmin_eleptvscutvars[1] = 0.;
4763 xmax_eleptvscutvars[1] = 3.141592/2;
4766 fHistoElePtvsCutVarsRS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsRS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4768 fHistoElePtvsCutVarsWS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsWS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4770 fHistoElePtvsCutVarsMCS[ih] =
new THnSparseF(Form(
"fHistoElePtvsCutVarsMCS[%d]",ih),
"",3,bins_eleptvscutvars,xmin_eleptvscutvars,xmax_eleptvscutvars);
4774 Int_t bins_eletag[3]= {20 ,40 ,10};
4775 Double_t xmin_eletag[3]={0.,0. ,0.0};
4776 Double_t xmax_eletag[3]={10.,20 ,100};
4777 fHistoElectronPi0Total =
new THnSparseF(
"fHistoElectronPi0Total",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4779 fHistoElectronPi0Tag =
new THnSparseF(
"fHistoElectronPi0Tag",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4781 fHistoElectronEtaTotal =
new THnSparseF(
"fHistoElectronEtaTotal",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4783 fHistoElectronEtaTag =
new THnSparseF(
"fHistoElectronEtaTag",
"",3,bins_eletag,xmin_eletag,xmax_eletag);
4797 if(!part)
return 0x0;
4798 if(!aod)
return 0x0;
4803 AliAODVertex *primVertexAOD;
4804 Bool_t unsetvtx = kFALSE;
4808 primVertexAOD =
fVtx1;
4813 primVertexAOD =
fVtx1;
4815 if(!primVertexAOD)
return 0x0;
4816 Double_t posprim[3]; primVertexAOD->GetXYZ(posprim);
4821 AliESDtrack *esdtrack =
new AliESDtrack((AliVTrack*)partpid);
4823 AliNeutralTrackParam *trackV0=NULL;
4824 const AliVTrack *trackVV0 =
dynamic_cast<const AliVTrack*
>(v0);
4825 if(trackVV0) trackV0 =
new AliNeutralTrackParam(trackVV0);
4827 Double_t xdummy, ydummy;
4828 Double_t dca = esdtrack->GetDCA(trackV0,
fBzkG,xdummy,ydummy);
4835 Double_t d0z0bach[2],covd0z0bach[3];
4836 if(sqrt(pow(secVert->GetX(),2)+pow(secVert->GetY(),2))<1.){
4837 part->PropagateToDCA(secVert,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4838 trackV0->PropagateToDCA(secVert,
fBzkG,kVeryBig);
4840 part->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4841 trackV0->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig);
4843 Double_t momv0_new[3]={-9999,-9999,-9999.};
4844 trackV0->GetPxPyPz(momv0_new);
4846 Double_t px[2],py[2],pz[2];
4847 px[0] = part->Px(); py[0] = part->Py(); pz[0] = part->Pz();
4848 px[1] = momv0_new[0]; py[1] = momv0_new[1]; pz[1] = momv0_new[2];
4853 Double_t d0[3],d0err[3];
4855 part->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0bach,covd0z0bach);
4857 d0err[0] = TMath::Sqrt(covd0z0bach[0]);
4859 Double_t d0z0v0[2],covd0z0v0[3];
4860 trackV0->PropagateToDCA(primVertexAOD,
fBzkG,kVeryBig,d0z0v0,covd0z0v0);
4862 d0err[1] = TMath::Sqrt(covd0z0v0[0]);
4867 Short_t
charge = part->Charge();
4871 if(unsetvtx)
delete primVertexAOD; primVertexAOD=NULL;
4872 if(esdtrack)
delete esdtrack;
4873 if(trackV0)
delete trackV0;
4877 UShort_t
id[2]={(UShort_t)part->GetID(),(UShort_t)trackV0->GetID()};
4880 theCascade->GetSecondaryVtx()->AddDaughter(part);
4881 theCascade->GetSecondaryVtx()->AddDaughter(v0);
4883 if(unsetvtx)
delete primVertexAOD; primVertexAOD=NULL;
4884 if(esdtrack)
delete esdtrack;
4885 if(trackV0)
delete trackV0;