Fix MC flagging for c-deuteron

fgrosa · fgrosa · commit 5795b0ecac29 · 2026-07-07T16:55:42.000+02:00
diff --git a/PWGHF/TableProducer/candidateCreator3Prong.cxx b/PWGHF/TableProducer/candidateCreator3Prong.cxx
@@ -1358,6 +1358,9 @@ struct HfCandidateCreator3ProngExpressions {
         auto bHadMother = mcParticles.rawIteratorAt(idxBhadMothers[0]);
         rowMcMatchRec(flagChannelMain, origin, swapping, flagChannelResonant, bHadMother.pt(), bHadMother.pdgCode(), nKinkedTracks, nInteractionsWithMaterial);
       } else {
+        if (std::abs(flagChannelMain) == DecayChannelMain::CDeuteronToDeKPi) {
+          origin == RecoDecay::OriginType::Prompt;
+        }
         rowMcMatchRec(flagChannelMain, origin, swapping, flagChannelResonant, -1.f, 0, nKinkedTracks, nInteractionsWithMaterial);
       }
     }
diff --git a/PWGHF/Utils/utilsMcMatching.h b/PWGHF/Utils/utilsMcMatching.h
@@ -338,7 +338,7 @@ inline void flipPdgSign(const int pdgMother, const int pdgToFlip, std::array<int
 template <typename Part>
 inline int getResonantDecayCDeuteron(Part const& particle)
 {
-  auto statusCode = std::abs(particle.getGenStatusCode());
+  auto statusCode = particle.getGenStatusCode();
   if (statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[0] || statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[1]) {
     return o2::hf_decay::hf_cand_3prong::DecayChannelResonant::CDeuteronToDeKstar0;
   }

Original file line number	Diff line number	Diff line change
`@@ -1358,6 +1358,9 @@ struct HfCandidateCreator3ProngExpressions {`
`1358`	`1358`	`auto bHadMother = mcParticles.rawIteratorAt(idxBhadMothers[0]);`
`1359`	`1359`	`rowMcMatchRec(flagChannelMain, origin, swapping, flagChannelResonant, bHadMother.pt(), bHadMother.pdgCode(), nKinkedTracks, nInteractionsWithMaterial);`
`1360`	`1360`	`} else {`
	`1361`	`+ if (std::abs(flagChannelMain) == DecayChannelMain::CDeuteronToDeKPi) {`
	`1362`	`+ origin == RecoDecay::OriginType::Prompt;`
	`1363`	`+ }`
`1361`	`1364`	`rowMcMatchRec(flagChannelMain, origin, swapping, flagChannelResonant, -1.f, 0, nKinkedTracks, nInteractionsWithMaterial);`
`1362`	`1365`	`}`
`1363`	`1366`	`}`
Original file line number	Diff line number	Diff line change
`@@ -338,7 +338,7 @@ inline void flipPdgSign(const int pdgMother, const int pdgToFlip, std::array<int`
`338`	`338`	`template <typename Part>`
`339`	`339`	`inline int getResonantDecayCDeuteron(Part const& particle)`
`340`	`340`	`{`
`341`		`- auto statusCode = std::abs(particle.getGenStatusCode());`
	`341`	`+ auto statusCode = particle.getGenStatusCode();`
`342`	`342`	`if (statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[0] \|\| statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[1]) {`
`343`	`343`	`return o2::hf_decay::hf_cand_3prong::DecayChannelResonant::CDeuteronToDeKstar0;`
`344`	`344`	`}`