diff --git a/PWGLF/Tasks/Resonances/k892analysis_PbPb.cxx b/PWGLF/Tasks/Resonances/k892analysis_PbPb.cxx index aeefe798b82..d1866ffb594 100644 --- a/PWGLF/Tasks/Resonances/k892analysis_PbPb.cxx +++ b/PWGLF/Tasks/Resonances/k892analysis_PbPb.cxx @@ -279,19 +279,19 @@ struct k892analysis_PbPb { if (TofandTpcPID) { - if (candidate.hasTOF() && std::abs(candidate.tofNSigmaKa()) < cMaxTOFnSigmaKaon && candidate.hasTPC() && std::abs(candidate.tpcNSigmaKa()) < cMaxTPCnSigmaKaon) { // tof and tpc cut + if (candidate.hasTOF() && std::abs(candidate.tofNSigmaKa()) <= cMaxTOFnSigmaKaon && candidate.hasTPC() && std::abs(candidate.tpcNSigmaKa()) <= cMaxTPCnSigmaKaon) { // tof and tpc cut return true; } } else { - if (candidate.hasTPC() && std::abs(candidate.tpcNSigmaKa()) < cMaxTPCnSigmaKaon) { // tpc cut, tof when available + if (candidate.hasTPC() && std::abs(candidate.tpcNSigmaKa()) <= cMaxTPCnSigmaKaon) { // tpc cut, tof when available if (cByPassTOF) // skip tof selection return true; if (candidate.hasTOF()) { - if (std::abs(candidate.tofNSigmaKa()) < cMaxTOFnSigmaKaon) { + if (std::abs(candidate.tofNSigmaKa()) <= cMaxTOFnSigmaKaon) { return true; } } else { @@ -309,19 +309,19 @@ struct k892analysis_PbPb { if (TofandTpcPID) { - if (candidate.hasTOF() && std::abs(candidate.tofNSigmaPi()) < cMaxTOFnSigmaPion && candidate.hasTPC() && std::abs(candidate.tpcNSigmaPi()) < cMaxTPCnSigmaPion) { // tof and tpc cut + if (candidate.hasTOF() && std::abs(candidate.tofNSigmaPi()) <= cMaxTOFnSigmaPion && candidate.hasTPC() && std::abs(candidate.tpcNSigmaPi()) <= cMaxTPCnSigmaPion) { // tof and tpc cut return true; } } else { - if (candidate.hasTPC() && std::abs(candidate.tpcNSigmaPi()) < cMaxTPCnSigmaPion) { // tpc cut, tof when available + if (candidate.hasTPC() && std::abs(candidate.tpcNSigmaPi()) <= cMaxTPCnSigmaPion) { // tpc cut, tof when available if (cByPassTOF) // skip tof selection return true; if (candidate.hasTOF()) { - if (std::abs(candidate.tofNSigmaPi()) < cMaxTOFnSigmaPion) { + if (std::abs(candidate.tofNSigmaPi()) <= cMaxTOFnSigmaPion) { return true; } } else { @@ -382,10 +382,10 @@ struct k892analysis_PbPb { if constexpr (IsMC) { if (tpclowpt) { - if (trk1ptPi > cMaxPtTPC || trk2ptKa > cMaxPtTPC) + if (trk1ptPi >= cMaxPtTPC || trk2ptKa >= cMaxPtTPC) continue; } else if (tofhighpt) { - if (trk1ptPi < cMinPtTOF || trk2ptKa < cMinPtTOF) + if (trk1ptPi <= cMinPtTOF || trk2ptKa <= cMinPtTOF) continue; } } @@ -552,8 +552,8 @@ struct k892analysis_PbPb { Filter collisionFilter = nabs(aod::collision::posZ) <= cfgCutVertex; Filter centralityFilter = nabs(aod::cent::centFT0C) <= cfgCutCentrality; - Filter acceptanceFilter = (nabs(aod::track::eta) < cfgCutEta && nabs(aod::track::pt) > cfgCutPT); - Filter DCAcutFilter = (nabs(aod::track::dcaXY) < cfgCutDCAxy) && (nabs(aod::track::dcaZ) < cfgCutDCAz); + Filter acceptanceFilter = (nabs(aod::track::eta) < cfgCutEta && nabs(aod::track::pt) >= cfgCutPT); + Filter DCAcutFilter = (nabs(aod::track::dcaXY) <= cfgCutDCAxy) && (nabs(aod::track::dcaZ) <= cfgCutDCAz); using EventCandidates = soa::Filtered>; using TrackCandidates = soa::Filtered cfgCutVertex) - return; + continue; histos.fill(HIST("hMCrecCollSels"), 5); auto centrality = RecCollision.centRun2V0M();