7272#include < type_traits>
7373#include < unordered_map>
7474#include < utility>
75+ #include < variant>
7576#include < vector>
7677
7778using namespace o2 ;
@@ -511,7 +512,7 @@ struct PhiStrangenessCorrelation {
511512 const bool applyK0sMassCut = (analysisMode == kDeltaYvsDeltaPhi ) && k0sConfigs.selectK0sInSigRegion ;
512513 const auto & [minMassK0s, maxMassK0s] = k0sConfigs.rangeMK0sSignal .value ;
513514 auto isK0sValid = [&](const auto & k0s) {
514- return (!applyEfficiency || k0s.pt () <= binspTK0S->back ()) && (!applyK0sMassCut || k0s.inMassRegion (minMassK0s, maxMassK0s));
515+ return (!applyEfficiency || k0s.pt () < binspTK0S->back ()) && (!applyK0sMassCut || k0s.inMassRegion (minMassK0s, maxMassK0s));
515516 };
516517
517518 const bool applyPionNSigmaCut = (analysisMode == kDeltaYvsDeltaPhi ) && pionConfigs.selectPionInSigRegion ;
@@ -520,11 +521,11 @@ struct PhiStrangenessCorrelation {
520521 // const float& tofPIDThreshold = pionConfigs.tofPIDThreshold;
521522
522523 auto isPionValid = [&](const auto & pion) {
523- return (!applyEfficiency || pion.pt () <= binspTPi->back ()) && (!applyPionNSigmaCut || pion.inNSigmaRegion (pidTPCMax, pidTOFMax));
524+ return (!applyEfficiency || pion.pt () < binspTPi->back ()) && (!applyPionNSigmaCut || pion.inNSigmaRegion (pidTPCMax, pidTOFMax));
524525 };
525526
526527 for (const auto & phiCand : phiCandidates) {
527- if (applyEfficiency && phiCand.pt () > binspTPhi->back ())
528+ if (applyEfficiency && phiCand.pt () >= binspTPhi->back ())
528529 continue ;
529530
530531 float weightPhi = computeWeight (BoundEfficiencyMap (effMapPhi, multiplicity, phiCand.pt (), phiCand.y ()));
@@ -651,7 +652,7 @@ struct PhiStrangenessCorrelation {
651652 const auto & [minMassK0s, maxMassK0s] = k0sConfigs.rangeMK0sSignal .value ;
652653
653654 auto isK0sValid = [&](const auto & k0s) {
654- return (!applyEfficiency || k0s.pt () <= binspTK0S->back ()) && (!applyK0sMassCut || k0s.inMassRegion (minMassK0s, maxMassK0s));
655+ return (!applyEfficiency || k0s.pt () < binspTK0S->back ()) && (!applyK0sMassCut || k0s.inMassRegion (minMassK0s, maxMassK0s));
655656 };
656657
657658 auto tuplePhiK0S = std::make_tuple (phiCandidates, k0sReduced);
@@ -662,7 +663,7 @@ struct PhiStrangenessCorrelation {
662663 float multiplicity = c1.centFT0M ();
663664
664665 for (const auto & [phiCand, k0s] : o2::soa::combinations (o2::soa::CombinationsFullIndexPolicy (phiCands, k0sRed))) {
665- if (applyEfficiency && phiCand.pt () > binspTPhi->back ())
666+ if (applyEfficiency && phiCand.pt () >= binspTPhi->back ())
666667 continue ;
667668 if (!isK0sValid (k0s))
668669 continue ;
@@ -746,7 +747,7 @@ struct PhiStrangenessCorrelation {
746747 // const float& tofPIDThreshold = pionConfigs.tofPIDThreshold;
747748
748749 auto isPionValid = [&](const auto & pion) {
749- return (!applyEfficiency || pion.pt () <= binspTPi->back ()) && (!applyPionNSigmaCut || pion.inNSigmaRegion (pidTPCMax, pidTOFMax));
750+ return (!applyEfficiency || pion.pt () < binspTPi->back ()) && (!applyPionNSigmaCut || pion.inNSigmaRegion (pidTPCMax, pidTOFMax));
750751 };
751752
752753 auto tuplePhiPion = std::make_tuple (phiCandidates, pionTracks);
@@ -757,7 +758,7 @@ struct PhiStrangenessCorrelation {
757758 float multiplicity = c1.centFT0M ();
758759
759760 for (const auto & [phiCand, piTrack] : o2::soa::combinations (o2::soa::CombinationsFullIndexPolicy (phiCands, piTracks))) {
760- if (applyEfficiency && phiCand.pt () > binspTPhi->back ())
761+ if (applyEfficiency && phiCand.pt () >= binspTPhi->back ())
761762 continue ;
762763 if (!isPionValid (piTrack))
763764 continue ;
0 commit comments