diff --git a/PWGCF/FemtoDream/Core/femtoDreamContainer.h b/PWGCF/FemtoDream/Core/femtoDreamContainer.h index 75582365f37..2b58f0e9c5c 100644 --- a/PWGCF/FemtoDream/Core/femtoDreamContainer.h +++ b/PWGCF/FemtoDream/Core/femtoDreamContainer.h @@ -270,7 +270,7 @@ class FemtoDreamContainer } if (extendedplots) { mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("/relPairkstarmTPtPart1PtPart2MultPercentile"), femtoObs, mT, part1.pt(), part2.pt(), multPercentile); - if constexpr (std::is_same_v && std::is_same_v) { + if constexpr (requires { part1.mLambda(); part2.mLambda(); }) { mHistogramRegistry->fill(HIST(mFolderSuffix[mEventType]) + HIST(o2::aod::femtodreamMCparticle::MCTypeName[mc]) + HIST("/invMassPart1invMassPart2kstar"), part1.mLambda(), part2.mLambda(), femtoObs); } } diff --git a/PWGCF/FemtoDream/Core/femtoDreamResoSelection.h b/PWGCF/FemtoDream/Core/femtoDreamResoSelection.h index b9bff7ac239..ad18bbf9ff8 100644 --- a/PWGCF/FemtoDream/Core/femtoDreamResoSelection.h +++ b/PWGCF/FemtoDream/Core/femtoDreamResoSelection.h @@ -65,7 +65,7 @@ class FemtoDreamResoSelection virtual ~FemtoDreamResoSelection() = default; template - uint32_t getType(V const& track1, V const& track2, bool resoIsNotAnti); + int getType(V const& track1, V const& track2, bool resoIsNotAnti); /// assigns value from configurbale to private class member template @@ -218,7 +218,7 @@ class FemtoDreamResoSelection }; // namespace femtoDream template -uint32_t FemtoDreamResoSelection::getType(V const& track1, V const& track2, bool resoIsNotAnti) +int FemtoDreamResoSelection::getType(V const& track1, V const& track2, bool resoIsNotAnti) { float posThresh = 0.; float negThresh = 0.; @@ -403,13 +403,9 @@ float FemtoDreamResoSelection::getNSigTotal(T const& track, V const& pid, float return std::abs(nSigTPC); } - float nSigTOF = 0.; - if (!track.hasTOF()) { - nSigTOF = 999.f; - } else { - nSigTOF = o2::aod::pidutils::tofNSigma(pid, track); - } - return std::sqrt(nSigTPC * nSigTPC + nSigTOF * nSigTOF); + float nSigTOF = track.hasTOF() ? o2::aod::pidutils::tofNSigma(pid, track) : 999.f; + + return std::hypot(nSigTPC, nSigTOF); } //// new getCutContainer