@@ -2307,7 +2307,7 @@ struct HfTrackIndexSkimCreator {
23072307 if constexpr (doPvRefit) {
23082308 // fill table row with coordinates of PV refit
23092309 rowProng2PVrefit (pvRefitCoord2Prong[0 ], pvRefitCoord2Prong[1 ], pvRefitCoord2Prong[2 ],
2310- pvRefitCovMatrix2Prong[0 ], pvRefitCovMatrix2Prong[1 ], pvRefitCovMatrix2Prong[2 ], pvRefitCovMatrix2Prong[3 ], pvRefitCovMatrix2Prong[4 ], pvRefitCovMatrix2Prong[5 ]);
2310+ pvRefitCovMatrix2Prong[0 ], pvRefitCovMatrix2Prong[1 ], pvRefitCovMatrix2Prong[2 ], pvRefitCovMatrix2Prong[3 ], pvRefitCovMatrix2Prong[4 ], pvRefitCovMatrix2Prong[5 ]);
23112311 }
23122312
23132313 if (debug) {
@@ -3017,10 +3017,9 @@ struct HfTrackIndexSkimCreatorCascades {
30173017 double massK0s{0 .};
30183018 double massPi{0 .};
30193019 double massLc{0 .};
3020- double mass2K0sP{0 .}; // WHY HERE?
30213020
30223021 Filter filterSelectCollisions = (aod::hf_sel_collision::whyRejectColl == 0 );
3023- Filter filterSelectTrackIds = (aod::hf_sel_track::isSelProng & static_cast <uint32_t >(BIT(CandidateType::CandV0bachelor))) != 0u && (applyProtonPid == false || (aod::hf_sel_track::isProton & static_cast <uint8_t >(BIT(ChannelsProtonPid::LcToPK0S))) != 0u );
3022+ Filter filterSelectTrackIds = (aod::hf_sel_track::isSelProng & static_cast <uint32_t >(BIT(CandidateType::CandV0bachelor))) != 0u && (applyProtonPid == false || (aod::hf_sel_track::isProton & static_cast <uint32_t >(BIT(ChannelsProtonPid::LcToPK0S))) != 0u );
30243023
30253024 using SelectedCollisions = soa::Filtered<soa::Join<aod::Collisions, aod::HfSelCollision>>;
30263025 using FilteredTrackAssocSel = soa::Filtered<soa::Join<aod::TrackAssoc, aod::HfSelTrack>>;
@@ -3165,7 +3164,7 @@ struct HfTrackIndexSkimCreatorCascades {
31653164
31663165 // invariant-mass cut: we do it here, before updating the momenta of bach and V0 during the fitting to save CPU
31673166 // TODO: but one should better check that the value here and after the fitter do not change significantly!!!
3168- mass2K0sP = RecoDecay::m (std::array{std::array{bach.px (), bach.py (), bach.pz ()}, momentumV0}, std::array{massP, massK0s});
3167+ double mass2K0sP = RecoDecay::m (std::array{std::array{bach.px (), bach.py (), bach.pz ()}, momentumV0}, std::array{massP, massK0s});
31693168 if ((cutInvMassCascLc >= 0 .) && (std::abs (mass2K0sP - massLc) > cutInvMassCascLc)) {
31703169 continue ;
31713170 }
0 commit comments