Skip to content

Commit 38cea88

Browse files
author
nkaratze
committed
Bug fix, removed signal splitting
1 parent f7f3381 commit 38cea88

File tree

1 file changed

+39
-71
lines changed

1 file changed

+39
-71
lines changed

‎PWGLF/Tasks/Strangeness/v0ptinvmassplots.cxx‎

Lines changed: 39 additions & 71 deletions
Original file line numberDiff line numberDiff line change
@@ -339,42 +339,42 @@ struct V0PtInvMassPlots {
339339
}
340340
rPtAnalysis.fill(HIST("hNEvents"), 1.5);
341341
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(2, "sel 8");
342-
if (!(collision.selection_bit(aod::evsel::kNoTimeFrameBorder) && doNoTimeFrameBorder)) {
342+
if (!(doNoTimeFrameBorder && collision.selection_bit(aod::evsel::kNoTimeFrameBorder))) {
343343
return false;
344344
}
345345
rPtAnalysis.fill(HIST("hNEvents"), 2.5);
346346
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(3, "NoTimeFrameBorder");
347-
if (!(collision.selection_bit(aod::evsel::kNoITSROFrameBorder) && doNoITSROFrameBorder)) {
347+
if (!(doNoITSROFrameBorder && collision.selection_bit(aod::evsel::kNoITSROFrameBorder))) {
348348
return false;
349349
}
350350
rPtAnalysis.fill(HIST("hNEvents"), 3.5);
351351
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(4, "NoITSROFrameBorder");
352-
if (!(collision.selection_bit(aod::evsel::kIsTriggerTVX) && doIsTriggerTVX)) {
352+
if (!(doIsTriggerTVX && collision.selection_bit(aod::evsel::kIsTriggerTVX))) {
353353
return false;
354354
}
355355
rPtAnalysis.fill(HIST("hNEvents"), 4.5);
356356
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(5, "IsTriggerTVX");
357-
if (!(std::abs(collision.posZ()) < cutZVertex && docutZVertex)) {
357+
if (!(docutZVertex && std::abs(collision.posZ()) < cutZVertex)) {
358358
return false;
359359
}
360360
rPtAnalysis.fill(HIST("hNEvents"), 5.5);
361361
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(6, "cutZVertex");
362-
if (!(collision.selection_bit(aod::evsel::kIsVertexTOFmatched) && doIsVertexTOFmatched)) {
362+
if (!(doIsVertexTOFmatched && collision.selection_bit(aod::evsel::kIsVertexTOFmatched))) {
363363
return false;
364364
}
365365
rPtAnalysis.fill(HIST("hNEvents"), 6.5);
366366
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(7, "IsVertexTOFmatched");
367-
if (!(collision.selection_bit(aod::evsel::kNoSameBunchPileup) && doNoSameBunchPileup)) {
367+
if (!(doNoSameBunchPileup && collision.selection_bit(aod::evsel::kNoSameBunchPileup))) {
368368
return false;
369369
}
370370
rPtAnalysis.fill(HIST("hNEvents"), 7.5);
371371
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(8, "NoSameBunchPileup");
372-
if (!(collision.selection_bit(aod::evsel::kIsVertexITSTPC) && doIsVertexITSTPC)) {
372+
if (!(doIsVertexITSTPC && collision.selection_bit(aod::evsel::kIsVertexITSTPC))) {
373373
return false;
374374
}
375375
rPtAnalysis.fill(HIST("hNEvents"), 8.5);
376376
rPtAnalysis.get<TH1>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(9, "IsVertexITSTPC");
377-
if (!(collision.isInelGt0() && doisInelGt0)) {
377+
if (!(doisInelGt0 && collision.isInelGt0())) {
378378
return false;
379379
}
380380
rPtAnalysis.fill(HIST("hNEvents"), 9.5);
@@ -394,22 +394,22 @@ struct V0PtInvMassPlots {
394394

395395
rPtAnalysis.fill(HIST("hNV0s"), 0.5);
396396
rPtAnalysis.get<TH1>(HIST("hNV0s"))->GetXaxis()->SetBinLabel(1, "All V0s");
397-
if (std::abs(v0.y()) > rapidityCut && doRapidityCut) { // V0 Rapidity Cut
397+
if (doRapidityCut && std::abs(v0.y()) > rapidityCut) { // V0 Rapidity Cut
398398
return false;
399399
}
400400
rPtAnalysis.fill(HIST("hNV0s"), 1.5);
401401
rPtAnalysis.get<TH1>(HIST("hNV0s"))->GetXaxis()->SetBinLabel(2, "Rapidity");
402-
if ((std::abs(posDaughterTrack.eta()) > etadau && std::abs(negDaughterTrack.eta()) > etadau) && doDaughterPseudorapidityCut) { // Daughters Pseudorapidity Cut
402+
if (doDaughterPseudorapidityCut && (std::abs(posDaughterTrack.eta()) > etadau && std::abs(negDaughterTrack.eta()) > etadau)) { // Daughters Pseudorapidity Cut
403403
return false;
404404
}
405405
rPtAnalysis.fill(HIST("hNV0s"), 2.5);
406406
rPtAnalysis.get<TH1>(HIST("hNV0s"))->GetXaxis()->SetBinLabel(3, "Dau Pseudorapidity");
407-
if ((posDaughterTrack.isITSAfterburner() || negDaughterTrack.isITSAfterburner()) && !doisITSAfterburner) { // ITS After Burner on daughter tracks
407+
if (!doisITSAfterburner && (posDaughterTrack.isITSAfterburner() || negDaughterTrack.isITSAfterburner())) { // ITS After Burner on daughter tracks
408408
return false;
409409
}
410410
rPtAnalysis.fill(HIST("hNV0s"), 3.5);
411411
rPtAnalysis.get<TH1>(HIST("hNV0s"))->GetXaxis()->SetBinLabel(4, "ITS Afterburner");
412-
if (posDaughterTrack.itsNCls() <= itsMinHits && negDaughterTrack.itsNCls() <= itsMinHits && doitsMinHits) { // Minimum hits in the ITS
412+
if (doitsMinHits && posDaughterTrack.itsNCls() <= itsMinHits && negDaughterTrack.itsNCls() <= itsMinHits) { // Minimum hits in the ITS
413413
return false;
414414
rPtAnalysis.fill(HIST("hNV0s"), 4.5);
415415
rPtAnalysis.get<TH1>(HIST("hNV0s"))->GetXaxis()->SetBinLabel(5, "ITS Min Hits");
@@ -431,55 +431,55 @@ struct V0PtInvMassPlots {
431431
rPtAnalysis.fill(HIST("hNK0sh"), 0.5);
432432
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(1, "All");
433433
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 0.5, v0.mK0Short());
434-
if ((std::abs(posDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion && std::abs(negDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion) && doK0shTPCPID) { // TPC PID for two pions
434+
if (doK0shTPCPID && (std::abs(posDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion && std::abs(negDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion)) { // TPC PID for two pions
435435
return false;
436436
}
437437
rPtAnalysis.fill(HIST("hNK0sh"), 1.5);
438438
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(2, "TPC_PID");
439439
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 1.5, v0.mK0Short());
440-
if (std::abs(v0.mLambda() - o2::constants::physics::MassLambda0) < compv0masscut && std::abs(v0.mAntiLambda() - o2::constants::physics::MassLambda0) < compv0masscut && doK0shcomptmasscut) { // Kzero competitive v0 mass cut (cut out Lambdas and Anti-Lambdas)
440+
if (doK0shcomptmasscut && std::abs(v0.mLambda() - o2::constants::physics::MassLambda0) < compv0masscut && std::abs(v0.mAntiLambda() - o2::constants::physics::MassLambda0) < compv0masscut) { // Kzero competitive v0 mass cut (cut out Lambdas and Anti-Lambdas)
441441
return false;
442442
}
443443
rPtAnalysis.fill(HIST("hNK0sh"), 2.5);
444444
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(3, "Compt_Mass");
445445
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 2.5, v0.mK0Short());
446-
if (v0.v0radius() > kaonshmaxct && doK0shMaxct) { // K0sh max ct
446+
if (doK0shMaxct && v0.v0radius() > kaonshmaxct) { // K0sh max ct
447447
return false;
448448
}
449449
rPtAnalysis.fill(HIST("hNK0sh"), 3.5);
450450
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(4, "Max_ct");
451451
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 3.5, v0.mK0Short());
452-
if (v0.qtarm() < (k0shparamArmenterosCut * std::abs(v0.alpha())) && doK0shArmenterosCut) { // K0sh Armenteros Cut
452+
if (doK0shArmenterosCut && v0.qtarm() < (k0shparamArmenterosCut * std::abs(v0.alpha()))) { // K0sh Armenteros Cut
453453
return false;
454454
}
455455
rPtAnalysis.fill(HIST("hNK0sh"), 4.5);
456456
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(5, "Armenteros");
457457
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 4.5, v0.mK0Short());
458-
if (v0.v0cosPA() < kaonshSettingcosPA && doK0shcosPACut) { // K0sh cosPA Topological Cut
458+
if (doK0shcosPACut && v0.v0cosPA() < kaonshSettingcosPA) { // K0sh cosPA Topological Cut
459459
return false;
460460
}
461461
rPtAnalysis.fill(HIST("hNK0sh"), 5.5);
462462
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(6, "cosPA");
463463
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 5.5, v0.mK0Short());
464-
if (v0.dcaV0daughters() > kaonshSettingdcav0dau && doK0shDCAdauCut) { // K0sh DCAdaughters Topological Cut
464+
if (doK0shDCAdauCut && v0.dcaV0daughters() > kaonshSettingdcav0dau) { // K0sh DCAdaughters Topological Cut
465465
return false;
466466
}
467467
rPtAnalysis.fill(HIST("hNK0sh"), 6.5);
468468
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(7, "DCAdau");
469469
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 6.5, v0.mK0Short());
470-
if (v0.v0radius() < kaonshSettingradius && doK0shv0radiusCut) { // K0sh v0radius Topological Cut
470+
if (doK0shv0radiusCut && v0.v0radius() < kaonshSettingradius) { // K0sh v0radius Topological Cut
471471
return false;
472472
}
473473
rPtAnalysis.fill(HIST("hNK0sh"), 7.5);
474474
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(8, "v0radius");
475475
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 7.5, v0.mK0Short());
476-
if (std::abs(v0.dcapostopv()) < kaonshSettingdcapostopv && doK0shdcaposdautopv) { // K0sh DCAPosDaughterToPV Topological Cut
476+
if (doK0shdcaposdautopv && std::abs(v0.dcapostopv()) < kaonshSettingdcapostopv) { // K0sh DCAPosDaughterToPV Topological Cut
477477
return false;
478478
}
479479
rPtAnalysis.fill(HIST("hNK0sh"), 8.5);
480480
rPtAnalysis.get<TH1>(HIST("hNK0sh"))->GetXaxis()->SetBinLabel(9, "DCAPosDautoPV");
481481
rPtAnalysis.fill(HIST("hMassK0ShortvsCuts"), 8.5, v0.mK0Short());
482-
if (std::abs(v0.dcanegtopv()) < kaonshSettingdcanegtopv && doK0shdcanegdautopv) { // K0sh DCANegDaughterToPV Topological Cut
482+
if (doK0shdcanegdautopv && std::abs(v0.dcanegtopv()) < kaonshSettingdcanegtopv) { // K0sh DCANegDaughterToPV Topological Cut
483483
return false;
484484
}
485485
rPtAnalysis.fill(HIST("hNK0sh"), 9.5);
@@ -508,55 +508,55 @@ struct V0PtInvMassPlots {
508508
rPtAnalysis.fill(HIST("hNLambda"), 0.5);
509509
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(1, "All");
510510
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 0.5, v0.mLambda());
511-
if (std::abs(posDaughterTrack.tpcNSigmaPr()) > nSigmaTPCProton && std::abs(negDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion && doLambdaTPCPID) { // TPC PID on daughter pion and proton for Lambda
511+
if (doLambdaTPCPID && std::abs(posDaughterTrack.tpcNSigmaPr()) > nSigmaTPCProton && std::abs(negDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion) { // TPC PID on daughter pion and proton for Lambda
512512
return false;
513513
}
514514
rPtAnalysis.fill(HIST("hNLambda"), 1.5);
515515
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(2, "TPC_PID");
516516
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 1.5, v0.mLambda());
517-
if (std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < compv0masscut && doLambdacomptmasscut) { // Lambda competitive v0 mass cut (cut out Kaons)
517+
if (doLambdacomptmasscut && std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < compv0masscut) { // Lambda competitive v0 mass cut (cut out Kaons)
518518
return false;
519519
}
520520
rPtAnalysis.fill(HIST("hNLambda"), 2.5);
521521
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(3, "Compt_Mass");
522522
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 2.5, v0.mLambda());
523-
if (v0.v0radius() > lambdamaxct && doLambdaMaxct) { // Lambda max ct
523+
if (doLambdaMaxct && v0.v0radius() > lambdamaxct) { // Lambda max ct
524524
return false;
525525
}
526526
rPtAnalysis.fill(HIST("hNLambda"), 3.5);
527527
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(4, "Max_ct");
528528
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 3.5, v0.mLambda());
529-
if (v0.qtarm() < (lambdaparamArmenterosCut * std::abs(v0.alpha())) && doLambdaArmenterosCut) { // Lambda Armenteros Cut
529+
if (doLambdaArmenterosCut && v0.qtarm() < (lambdaparamArmenterosCut * std::abs(v0.alpha()))) { // Lambda Armenteros Cut
530530
return false;
531531
}
532532
rPtAnalysis.fill(HIST("hNLambda"), 4.5);
533533
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(5, "Armenteros");
534534
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 4.5, v0.mLambda());
535-
if (v0.v0cosPA() < lambdaSettingcosPA && doLambdacosPACut) { // Lambda cosPA Topological Cut
535+
if (doLambdacosPACut && v0.v0cosPA() < lambdaSettingcosPA) { // Lambda cosPA Topological Cut
536536
return false;
537537
}
538538
rPtAnalysis.fill(HIST("hNLambda"), 5.5);
539539
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(6, "cosPA");
540540
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 5.5, v0.mLambda());
541-
if (v0.dcaV0daughters() > lambdaSettingdcav0dau && doLambdaDCAdauCut) { // Lambda DCAdaughters Topological Cut
541+
if (doLambdaDCAdauCut && v0.dcaV0daughters() > lambdaSettingdcav0dau) { // Lambda DCAdaughters Topological Cut
542542
return false;
543543
}
544544
rPtAnalysis.fill(HIST("hNLambda"), 6.5);
545545
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(7, "DCAdau");
546546
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 6.5, v0.mLambda());
547-
if (v0.v0radius() < lambdaSettingradius && doLambdav0radiusCut) { // Lambda v0radius Topological Cut
547+
if (doLambdav0radiusCut && v0.v0radius() < lambdaSettingradius) { // Lambda v0radius Topological Cut
548548
return false;
549549
}
550550
rPtAnalysis.fill(HIST("hNLambda"), 7.5);
551551
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(8, "v0radius");
552552
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 7.5, v0.mLambda());
553-
if (std::abs(v0.dcapostopv()) < lambdaSettingdcapostopv && doLambdadcaposdautopv) { // Lambda DCAPosDaughterToPV Topological Cut
553+
if (doLambdadcaposdautopv && std::abs(v0.dcapostopv()) < lambdaSettingdcapostopv) { // Lambda DCAPosDaughterToPV Topological Cut
554554
return false;
555555
}
556556
rPtAnalysis.fill(HIST("hNLambda"), 8.5);
557557
rPtAnalysis.get<TH1>(HIST("hNLambda"))->GetXaxis()->SetBinLabel(9, "DCAPosDautoPV");
558558
rPtAnalysis.fill(HIST("hMassLambdavsCuts"), 8.5, v0.mLambda());
559-
if (std::abs(v0.dcanegtopv()) < lambdaSettingdcanegtopv && doLambdadcanegdautopv) { // Lambda DCANegDaughterToPV Topological Cut
559+
if (doLambdadcanegdautopv && std::abs(v0.dcanegtopv()) < lambdaSettingdcanegtopv) { // Lambda DCANegDaughterToPV Topological Cut
560560
return false;
561561
}
562562
rPtAnalysis.fill(HIST("hNLambda"), 9.5);
@@ -585,55 +585,55 @@ struct V0PtInvMassPlots {
585585
rPtAnalysis.fill(HIST("hNAntilambda"), 0.5);
586586
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(1, "All");
587587
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 0.5, v0.mAntiLambda());
588-
if (std::abs(negDaughterTrack.tpcNSigmaPr()) > nSigmaTPCProton && std::abs(posDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion) { // TPC PID on daughter pion and proton for AntiLambda
588+
if (doAntilambdaTPCPID && std::abs(negDaughterTrack.tpcNSigmaPr()) > nSigmaTPCProton && std::abs(posDaughterTrack.tpcNSigmaPi()) > nSigmaTPCPion) { // TPC PID on daughter pion and proton for AntiLambda
589589
return false;
590590
}
591591
rPtAnalysis.fill(HIST("hNAntilambda"), 1.5);
592592
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(2, "TPC_PID");
593593
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 1.5, v0.mAntiLambda());
594-
if (std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < compv0masscut && doAntilambdacomptmasscut) { // Antilambda competitive v0 mass cut (cut out Kaons)
594+
if (doAntilambdacomptmasscut && std::abs(v0.mK0Short() - o2::constants::physics::MassK0Short) < compv0masscut) { // Antilambda competitive v0 mass cut (cut out Kaons)
595595
return false;
596596
}
597597
rPtAnalysis.fill(HIST("hNAntilambda"), 2.5);
598598
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(3, "Compt_Mass");
599599
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 2.5, v0.mAntiLambda());
600-
if (v0.v0radius() > antilambdamaxct && doAntilambdaMaxct) { // Antilambda max ct
600+
if (doAntilambdaMaxct && v0.v0radius() > antilambdamaxct) { // Antilambda max ct
601601
return false;
602602
}
603603
rPtAnalysis.fill(HIST("hNAntilambda"), 3.5);
604604
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(4, "Max_ct");
605605
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 3.5, v0.mAntiLambda());
606-
if (v0.qtarm() < (antilambdaparamArmenterosCut * std::abs(v0.alpha())) && doAntilambdaArmenterosCut) { // Antilambda Armenteros Cut
606+
if (doAntilambdaArmenterosCut && v0.qtarm() < (antilambdaparamArmenterosCut * std::abs(v0.alpha()))) { // Antilambda Armenteros Cut
607607
return false;
608608
}
609609
rPtAnalysis.fill(HIST("hNAntilambda"), 4.5);
610610
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(5, "Armenteros");
611611
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 4.5, v0.mAntiLambda());
612-
if (v0.v0cosPA() < antilambdaSettingcosPA && doAntilambdacosPACut) { // Antilambda cosPA Topological Cut
612+
if (doAntilambdacosPACut && v0.v0cosPA() < antilambdaSettingcosPA) { // Antilambda cosPA Topological Cut
613613
return false;
614614
}
615615
rPtAnalysis.fill(HIST("hNAntilambda"), 5.5);
616616
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(6, "cosPA");
617617
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 5.5, v0.mAntiLambda());
618-
if (v0.dcaV0daughters() > antilambdaSettingdcav0dau && doAntilambdaDCAdauCut) { // Antilambda DCAdaughters Topological Cut
618+
if (doAntilambdaDCAdauCut && v0.dcaV0daughters() > antilambdaSettingdcav0dau) { // Antilambda DCAdaughters Topological Cut
619619
return false;
620620
}
621621
rPtAnalysis.fill(HIST("hNAntilambda"), 6.5);
622622
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(7, "DCAdau");
623623
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 6.5, v0.mAntiLambda());
624-
if (v0.v0radius() < antilambdaSettingradius && doAntilambdav0radiusCut) { // Antilambda v0radius Topological Cut
624+
if (doAntilambdav0radiusCut && v0.v0radius() < antilambdaSettingradius) { // Antilambda v0radius Topological Cut
625625
return false;
626626
}
627627
rPtAnalysis.fill(HIST("hNAntilambda"), 7.5);
628628
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(8, "v0radius");
629629
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 7.5, v0.mAntiLambda());
630-
if (std::abs(v0.dcapostopv()) < antilambdaSettingdcapostopv && doAntilambdadcaposdautopv) { // Antilambda DCAPosDaughterToPV Topological Cut
630+
if (doAntilambdadcaposdautopv && std::abs(v0.dcapostopv()) < antilambdaSettingdcapostopv) { // Antilambda DCAPosDaughterToPV Topological Cut
631631
return false;
632632
}
633633
rPtAnalysis.fill(HIST("hNAntilambda"), 8.5);
634634
rPtAnalysis.get<TH1>(HIST("hNAntilambda"))->GetXaxis()->SetBinLabel(9, "DCAPosDautoPV");
635635
rPtAnalysis.fill(HIST("hMassAntilambdavsCuts"), 8.5, v0.mAntiLambda());
636-
if (std::abs(v0.dcanegtopv()) < antilambdaSettingdcanegtopv && doAntilambdadcanegdautopv) { // Antilambda DCANegDaughterToPV Topological Cut
636+
if (doAntilambdadcanegdautopv && std::abs(v0.dcanegtopv()) < antilambdaSettingdcanegtopv) { // Antilambda DCANegDaughterToPV Topological Cut
637637
return false;
638638
}
639639
rPtAnalysis.fill(HIST("hNAntilambda"), 9.5);
@@ -767,43 +767,11 @@ struct V0PtInvMassPlots {
767767
}
768768
rMCCorrections.fill(HIST("hNRecEvents_MC"), 0.5); // Event Split Numenator
769769

770-
// v0 Signal Splitting Numenator Start
771-
for (const auto& mcParticle : mcParticles) {
772-
if (mcParticle.isPhysicalPrimary()) {
773-
if (std::abs(mcParticle.y()) < rapidityCut) {
774-
if (mcParticle.pdgCode() == kK0Short) { // kzero matched
775-
rMCCorrections.fill(HIST("hK0ShSplitNumenatorPtSpectrum"), mcParticle.pt());
776-
}
777-
if (mcParticle.pdgCode() == kLambda0) { // lambda matched
778-
rMCCorrections.fill(HIST("hLambdaSplitNumenatorPtSpectrum"), mcParticle.pt());
779-
}
780-
if (mcParticle.pdgCode() == kLambda0Bar) { // antilambda matched
781-
rMCCorrections.fill(HIST("hAntilambdaSplitNumenatorPtSpectrum"), mcParticle.pt());
782-
}
783-
}
784-
}
785-
}
786-
// V0 Signal Splitting Numenator End
787-
788770
for (const auto& v0 : V0s) {
789771
// Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
790772
if (v0.has_mcParticle()) {
791773
auto v0mcParticle = v0.mcParticle();
792774

793-
// signal splitting demoninator
794-
if (v0mcParticle.isPhysicalPrimary()) {
795-
if (v0mcParticle.pdgCode() == kK0Short) { // kzero matched
796-
rMCCorrections.fill(HIST("hK0ShSplitDenominatorPtSpectrum"), v0mcParticle.pt());
797-
}
798-
if (v0mcParticle.pdgCode() == kLambda0) { // lambda matched
799-
rMCCorrections.fill(HIST("hLambdaSplitDenominatorPtSpectrum"), v0mcParticle.pt());
800-
}
801-
if (v0mcParticle.pdgCode() == kLambda0Bar) { // antilambda matched
802-
rMCCorrections.fill(HIST("hAntilambdaSplitDenominatorPtSpectrum"), v0mcParticle.pt());
803-
}
804-
}
805-
// signal splitting demoninator end
806-
807775
if (!acceptV0(v0)) { // V0 Selections
808776
continue;
809777
}

0 commit comments

Comments
 (0)