Skip to content

Commit 68c3dbf

Browse files
[PWGHF/D2H] Correcting delta mass definition in McRec proccess function (#13792)
Co-authored-by: Fabrizio <fabrizio.grosa@cern.ch>
1 parent 39a9c37 commit 68c3dbf

File tree

1 file changed

+6
-1
lines changed

1 file changed

+6
-1
lines changed

PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx

Lines changed: 6 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -569,7 +569,12 @@ struct HfTaskDstarToD0Pi {
569569
} else if (studyD0ToPiKPi0 && candDstarMcRec.isSelDstarToD0Pi() && (std::abs(candDstarMcRec.flagMcMatchRec()) == hf_decay::hf_cand_dstar::DecayChannelMain::DstarToPiKPiPi0) && (std::abs(candDstarMcRec.flagMcMatchRecD0()) == hf_decay::hf_cand_2prong::DecayChannelMain::D0ToPiKPi0)) {
570570
// Aplly all selection to study D*->D0(piKpi0)pi channel same as signal channel
571571
// MC Matched but to D*->D0(piKpi0)pi channel
572-
auto deltaMDstar = std::abs(candDstarMcRec.invMassDstar() - candDstarMcRec.invMassD0());
572+
double deltaMDstar = -999.;
573+
if (candDstarMcRec.signSoftPi() < 0) {
574+
deltaMDstar = candDstarMcRec.invMassAntiDstar() - candDstarMcRec.invMassD0Bar();
575+
} else {
576+
deltaMDstar = candDstarMcRec.invMassDstar() - candDstarMcRec.invMassD0();
577+
}
573578
if constexpr (ApplyMl) {
574579
auto bdtScore = candDstarMcRec.mlProbDstarToD0Pi();
575580
// inclusive study

0 commit comments

Comments
 (0)