Skip to content

[PWGHF] Fix MC flagging for c-deuteron #16953

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fgrosa merged 3 commits into from
Jul 7, 2026
+10 −7
Merged

[PWGHF] Fix MC flagging for c-deuteron #16953

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
5795b0e
Fix MC flagging for c-deuteron
fgrosa Jul 7, 2026
d6b4b86
Fix variable assigment and camel case
fgrosa Jul 7, 2026
fcc840e
Use constexpr for const arrays
fgrosa Jul 7, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 3 additions & 0 deletions PWGHF/TableProducer/candidateCreator3Prong.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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);
}
}
Expand Down
14 changes: 7 additions & 7 deletions PWGHF/Utils/utilsMcMatching.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -160,9 +160,9 @@ static const std::unordered_map<DecayChannelMain, const std::vector<int>> daught
{DecayChannelMain::CDeuteronToDeKPi, {+o2::constants::physics::Pdg::kDeuteron, +PDG_t::kKMinus, +PDG_t::kPiPlus}}};

/// resonances in c-deuteron decay are not stored in the particle stack for c-deuteron, but tagged with specific status codes
static const std::array<int, 2> StatusCodeCDeuteronToDeKstar0{-85, -95};
static const std::array<int, 2> StatusCodeCDeuteronToNeDeltaplusK{-86, -96};
static const std::array<int, 2> StatusCodeCDeuteronToNeL1520Pi{-87, -97};
static constexpr std::array<int, 2> statusCodeCDeuteronToDeKstar0{-85, -95};
static constexpr std::array<int, 2> statusCodeCDeuteronToNeDeltaplusK{-86, -96};
static constexpr std::array<int, 2> statusCodeCDeuteronToNeL1520Pi{-87, -97};

/// Returns a map of the possible final states for a specific 3-prong particle specie
/// \param pdgMother PDG code of the mother particle
Expand Down Expand Up @@ -338,14 +338,14 @@ 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());
if (statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[0] || statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToDeKstar0[1]) {
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;
}
if (statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToNeDeltaplusK[0] || statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToNeDeltaplusK[1]) {
if (statusCode == o2::hf_decay::hf_cand_3prong::statusCodeCDeuteronToNeDeltaplusK[0] || statusCode == o2::hf_decay::hf_cand_3prong::statusCodeCDeuteronToNeDeltaplusK[1]) {
return o2::hf_decay::hf_cand_3prong::DecayChannelResonant::CDeuteronToNeDeltaplusK;
}
if (statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToNeL1520Pi[0] || statusCode == o2::hf_decay::hf_cand_3prong::StatusCodeCDeuteronToNeL1520Pi[1]) {
if (statusCode == o2::hf_decay::hf_cand_3prong::statusCodeCDeuteronToNeL1520Pi[0] || statusCode == o2::hf_decay::hf_cand_3prong::statusCodeCDeuteronToNeL1520Pi[1]) {
return o2::hf_decay::hf_cand_3prong::DecayChannelResonant::CDeuteronToNeL1520Pi;
}
return 0;
Expand Down
Loading