Skip to content

Commit 7f0818e

Browse files
scannitoscannito
committed
Try again
1 parent 12c308f commit 7f0818e

1 file changed

Lines changed: 98 additions & 28 deletions

File tree

‎PWGLF/Tasks/Strangeness/phiStrangeCorrelation.cxx‎

Lines changed: 98 additions & 28 deletions
Original file line numberDiff line numberDiff line change
@@ -648,6 +648,76 @@ struct PhiStrangeCorrelation {
648648
const std::array<std::pair<float, float>, kPhiMassRegions> phiMassRegions = {phiConfigs.rangeMPhiSignal, phiConfigs.rangeMPhiSideband};
649649

650650
if (analysisMode == kMassvsMass) {
651+
if constexpr (PartType == kK0S) {
652+
if constexpr (IsME) {
653+
customFillTHn(HIST("phiK0S/h6PhiK0SDataME"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.m());
654+
} else {
655+
customFillTHn(HIST("phiK0S/h6PhiK0SData"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.m());
656+
}
657+
} else if constexpr (PartType == kXi) {
658+
if constexpr (IsME) {
659+
customFillTHn(HIST("phiXi/h6PhiXiDataME"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.m());
660+
} else {
661+
customFillTHn(HIST("phiXi/h6PhiXiData"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.m());
662+
}
663+
} else if constexpr (PartType == kPion) {
664+
if constexpr (IsME) {
665+
customFillTHn(HIST("phiPi/h6PhiPiTPCDataME"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.nSigmaTPC());
666+
customFillTHn(HIST("phiPi/h6PhiPiTOFDataME"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.nSigmaTOF());
667+
} else {
668+
customFillTHn(HIST("phiPi/h6PhiPiTPCData"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.nSigmaTPC());
669+
customFillTHn(HIST("phiPi/h6PhiPiTOFData"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assoc.nSigmaTOF());
670+
}
671+
}
672+
} else if (analysisMode == kDeltaYvsDeltaPhi) {
673+
static_for<0, kPhiMassRegions - 1>([&](auto i_idx) {
674+
constexpr unsigned int Idx = i_idx.value;
675+
const auto& [minMassPhi, maxMassPhi] = phiMassRegions[Idx];
676+
if (!phiCand.inMassRegion(minMassPhi, maxMassPhi)) {
677+
return;
678+
}
679+
680+
if constexpr (PartType == kK0S) {
681+
if constexpr (IsME) {
682+
if constexpr (Idx == kSignalRegion)
683+
customFillTHn(HIST("phiK0S/h5PhiK0SDataMESignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
684+
else
685+
customFillTHn(HIST("phiK0S/h5PhiK0SDataMESideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
686+
} else {
687+
if constexpr (Idx == kSignalRegion)
688+
customFillTHn(HIST("phiK0S/h5PhiK0SDataSignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
689+
else
690+
customFillTHn(HIST("phiK0S/h5PhiK0SDataSideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
691+
}
692+
} else if constexpr (PartType == kXi) {
693+
if constexpr (IsME) {
694+
if constexpr (Idx == kSignalRegion)
695+
customFillTHn(HIST("phiXi/h5PhiXiDataMESignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
696+
else
697+
customFillTHn(HIST("phiXi/h5PhiXiDataMESideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
698+
} else {
699+
if constexpr (Idx == kSignalRegion)
700+
customFillTHn(HIST("phiXi/h5PhiXiDataSignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
701+
else
702+
customFillTHn(HIST("phiXi/h5PhiXiDataSideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
703+
}
704+
} else if constexpr (PartType == kPion) {
705+
if constexpr (IsME) {
706+
if constexpr (Idx == kSignalRegion)
707+
customFillTHn(HIST("phiPi/h5PhiPiDataMESignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
708+
else
709+
customFillTHn(HIST("phiPi/h5PhiPiDataMESideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
710+
} else {
711+
if constexpr (Idx == kSignalRegion)
712+
customFillTHn(HIST("phiPi/h5PhiPiDataSignal"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
713+
else
714+
customFillTHn(HIST("phiPi/h5PhiPiDataSideband"), weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), getDeltaPhi(phiCand.phi(), assoc.phi()));
715+
}
716+
}
717+
});
718+
}
719+
720+
/*if (analysisMode == kMassvsMass) {
651721
auto fillMass = [&](auto histID, auto assocVal) {
652722
customFillTHn(histID, weight, multiplicity, phiCand.pt(), assoc.pt(), phiCand.y() - assoc.y(), phiCand.m(), assocVal);
653723
};
@@ -721,35 +791,35 @@ struct PhiStrangeCorrelation {
721791
else
722792
fillDelta(HIST("phiPi/h5PhiPiDataSideband"));
723793
}
724-
}
725-
726-
/*if constexpr (PartType == kK0S) {
727-
if constexpr (IsME) {
728-
auto t = std::make_tuple(HIST("phiK0S/h5PhiK0SDataMESignal"), HIST("phiK0S/h5PhiK0SDataMESideband"));
729-
fillDelta(std::get<Idx>(t));
730-
} else {
731-
auto t = std::make_tuple(HIST("phiK0S/h5PhiK0SDataSignal"), HIST("phiK0S/h5PhiK0SDataSideband"));
732-
fillDelta(std::get<Idx>(t));
733-
}
734-
} else if constexpr (PartType == kXi) {
735-
if constexpr (IsME) {
736-
auto t = std::make_tuple(HIST("phiXi/h5PhiXiDataMESignal"), HIST("phiXi/h5PhiXiDataMESideband"));
737-
fillDelta(std::get<Idx>(t));
738-
} else {
739-
auto t = std::make_tuple(HIST("phiXi/h5PhiXiDataSignal"), HIST("phiXi/h5PhiXiDataSideband"));
740-
fillDelta(std::get<Idx>(t));
741-
}
742-
} else if constexpr (PartType == kPion) {
743-
if constexpr (IsME) {
744-
auto t = std::make_tuple(HIST("phiPi/h5PhiPiDataMESignal"), HIST("phiPi/h5PhiPiDataMESideband"));
745-
fillDelta(std::get<Idx>(t));
746-
} else {
747-
auto t = std::make_tuple(HIST("phiPi/h5PhiPiDataSignal"), HIST("phiPi/h5PhiPiDataSideband"));
748-
fillDelta(std::get<Idx>(t));
749-
}
750794
}*/
751-
});
752-
}
795+
796+
/*if constexpr (PartType == kK0S) {
797+
if constexpr (IsME) {
798+
auto t = std::make_tuple(HIST("phiK0S/h5PhiK0SDataMESignal"), HIST("phiK0S/h5PhiK0SDataMESideband"));
799+
fillDelta(std::get<Idx>(t));
800+
} else {
801+
auto t = std::make_tuple(HIST("phiK0S/h5PhiK0SDataSignal"), HIST("phiK0S/h5PhiK0SDataSideband"));
802+
fillDelta(std::get<Idx>(t));
803+
}
804+
} else if constexpr (PartType == kXi) {
805+
if constexpr (IsME) {
806+
auto t = std::make_tuple(HIST("phiXi/h5PhiXiDataMESignal"), HIST("phiXi/h5PhiXiDataMESideband"));
807+
fillDelta(std::get<Idx>(t));
808+
} else {
809+
auto t = std::make_tuple(HIST("phiXi/h5PhiXiDataSignal"), HIST("phiXi/h5PhiXiDataSideband"));
810+
fillDelta(std::get<Idx>(t));
811+
}
812+
} else if constexpr (PartType == kPion) {
813+
if constexpr (IsME) {
814+
auto t = std::make_tuple(HIST("phiPi/h5PhiPiDataMESignal"), HIST("phiPi/h5PhiPiDataMESideband"));
815+
fillDelta(std::get<Idx>(t));
816+
} else {
817+
auto t = std::make_tuple(HIST("phiPi/h5PhiPiDataSignal"), HIST("phiPi/h5PhiPiDataSideband"));
818+
fillDelta(std::get<Idx>(t));
819+
}
820+
}*/
821+
//});
822+
//}
753823
}
754824

755825
template <AssociatedParticleType PartType, bool IsME, typename TPhiCand, typename TAssoc>

0 commit comments

Comments
 (0)