From 90fd4917c57b007245c3dc923d8a7399fe0ce866 Mon Sep 17 00:00:00 2001 From: Deependra Sharma Date: Tue, 11 Nov 2025 22:36:05 +0000 Subject: [PATCH 1/2] Correcting delta mass definition in McRec proccess function --- PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx b/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx index 7bc69aaaf0f..0b4e388c69d 100644 --- a/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx +++ b/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx @@ -569,7 +569,12 @@ struct HfTaskDstarToD0Pi { } 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)) { // Aplly all selection to study D*->D0(piKpi0)pi channel same as signal channel // MC Matched but to D*->D0(piKpi0)pi channel - auto deltaMDstar = std::abs(candDstarMcRec.invMassDstar() - candDstarMcRec.invMassD0()); + auto deltaMDstar = -999.; + if (candDstarMcRec.signSoftPi() < 0) { + deltaMDstar = candDstarMcRec.invMassAntiDstar() - candDstarMcRec.invMassD0Bar(); + } else { + deltaMDstar = candDstarMcRec.invMassDstar() - candDstarMcRec.invMassD0(); + } if constexpr (ApplyMl) { auto bdtScore = candDstarMcRec.mlProbDstarToD0Pi(); // inclusive study From 24bd851f5fc6b59ce73299ae469de468343b74c0 Mon Sep 17 00:00:00 2001 From: Fabrizio Date: Wed, 12 Nov 2025 13:36:54 +0100 Subject: [PATCH 2/2] Update PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx --- PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx b/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx index 0b4e388c69d..c6d11b6d5a7 100644 --- a/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx +++ b/PWGHF/D2H/Tasks/taskDstarToD0Pi.cxx @@ -569,7 +569,7 @@ struct HfTaskDstarToD0Pi { } 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)) { // Aplly all selection to study D*->D0(piKpi0)pi channel same as signal channel // MC Matched but to D*->D0(piKpi0)pi channel - auto deltaMDstar = -999.; + double deltaMDstar = -999.; if (candDstarMcRec.signSoftPi() < 0) { deltaMDstar = candDstarMcRec.invMassAntiDstar() - candDstarMcRec.invMassD0Bar(); } else {