Skip to content

[PWGHF] Some fixes in SigmaC MC matching #14093

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
mfaggin merged 1 commit into from
Dec 4, 2025
Merged

[PWGHF] Some fixes in SigmaC MC matching #14093

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
43 changes: 40 additions & 3 deletions PWGHF/D2H/Tasks/taskSigmac.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -585,6 +585,24 @@ struct HfTaskSigmac {
} /// end THn for candidate Λc+ cut variation w/o Σc0,++ mass-window cut
}; /// end fillHistosData

/// @brief function to remap the value of the resonant decay channel to fit the binning of the thnAxisChannel axis
/// @param channel the value obtained from candidateLc.flagMcDecayChanGen() or particleLc.flagMcDecayChanGen()
int remapResoChannelLc(int channel)
{
switch (channel) {
case 0:
// direct channel
return 0;
case o2::hf_decay::hf_cand_3prong::DecayChannelResonant::LcToPKstar0:
return 1;
case o2::hf_decay::hf_cand_3prong::DecayChannelResonant::LcToDeltaplusplusK:
return 2;
case o2::hf_decay::hf_cand_3prong::DecayChannelResonant::LcToL1520Pi:
return 3;
}
return -1;
}

/// @brief function to fill the histograms needed in analysis (MC)
/// @param candidatesSc are the reconstructed candidate Σc0,++ with MC info
/// @param mcParticles are the generated particles with flags wheter they are Σc0,++ or not
Expand Down Expand Up @@ -679,6 +697,7 @@ struct HfTaskSigmac {
etaGenSoftPi = daugSoftPi.eta();
phiGenSoftPi = daugSoftPi.phi();
}
channel = remapResoChannelLc(channel);

/// Fill histograms
int sigmacSpecies = -1;
Expand Down Expand Up @@ -770,7 +789,8 @@ struct HfTaskSigmac {
}
double ptGenLc(particle.pt()), ptGenLcBMother(-1.);
int const origin = particle.originMcGen();
int const channel = particle.flagMcDecayChanGen();
int channel = particle.flagMcDecayChanGen();
channel = remapResoChannelLc(channel);
if (origin == RecoDecay::OriginType::Prompt) {
registry.fill(HIST("MC/generated/hnLambdaCGen"), ptGenLc, ptGenLcBMother, origin, channel);
} else {
Expand Down Expand Up @@ -802,7 +822,21 @@ struct HfTaskSigmac {
bool const isTrueScStar0Reco = std::abs(candSc.flagMcMatchRec()) == o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi;
bool const isTrueScPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi;
bool const isTrueScStarPlusPlusReco = std::abs(candSc.flagMcMatchRec()) == o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi;
if (!isTrueSc0Reco && !isTrueScStar0Reco && !isTrueScPlusPlusReco && !isTrueScStarPlusPlusReco) {
continue;
}
int sigmacSpecies = -1;

/// debug
if ((isTrueSc0Reco || isTrueScStar0Reco) && chargeSc != o2::aod::hf_cand_sigmac::ChargeNull) {
/// this should never happen
LOG(fatal) << "isTrueSc0Reco=" << isTrueSc0Reco << ", isTrueScStar0Reco=" << isTrueScStar0Reco << ", but chargeSc = " << static_cast<int>(chargeSc) << "! Not possible, abort...";
}
if ((isTrueScPlusPlusReco || isTrueScStarPlusPlusReco) && std::abs(chargeSc) != o2::aod::hf_cand_sigmac::ChargePlusPlus) {
/// this should never happen
LOG(fatal) << "isTrueScPlusPlusReco=" << isTrueScPlusPlusReco << ", isTrueScStarPlusPlusReco=" << isTrueScStarPlusPlusReco << ", but chargeSc = " << static_cast<int>(chargeSc) << "! Not possible, abort...";
}

if ((isTrueSc0Reco || isTrueScStar0Reco) && (chargeSc == o2::aod::hf_cand_sigmac::ChargeNull)) {
/// Reconstructed Σc0 signal
// Get the corresponding MC particle for Sc, found as the mother of the soft pion
Expand All @@ -820,7 +854,7 @@ struct HfTaskSigmac {
// Get the corresponding MC particle for Lc
auto arrayDaughtersLc = std::array{candidateLc.template prong0_as<aod::TracksWMc>(), candidateLc.template prong1_as<aod::TracksWMc>(), candidateLc.template prong2_as<aod::TracksWMc>()};
int8_t sign = 0;
int const indexMcLcRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersLc, o2::constants::physics::Pdg::kLambdaCPlus, std::array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
int const indexMcLcRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughtersLc, o2::constants::physics::Pdg::kLambdaCPlus, std::array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
auto particleLc = mcParticles.rawIteratorAt(indexMcLcRec);
// Get the corresponding MC particle for soft pion
auto particleSoftPi = candSc.prong1_as<aod::TracksWMc>().mcParticle();
Expand All @@ -837,6 +871,7 @@ struct HfTaskSigmac {
double cpaLc(candidateLc.cpa()), cpaXYLc(candidateLc.cpaXY());
int const origin = candSc.originMcRec();
auto channel = candidateLc.flagMcDecayChanRec(); /// 0: direct; 1: Λc± → p± K*; 2: Λc± → Δ(1232)±± K∓; 3: Λc± → Λ(1520) π±
channel = remapResoChannelLc(channel);

/// candidate Λc+ → pK-π+ (and charge conjugate) within the range of M(pK-π+) chosen in the Σc0,++ builder
if ((TESTBIT(isCandPKPiPiKP, o2::aod::hf_cand_sigmac::Decays::PKPi)) && std::abs(candidateLc.template prong0_as<aod::TracksWMc>().mcParticle().pdgCode()) == kProton) {
Expand Down Expand Up @@ -1026,7 +1061,7 @@ struct HfTaskSigmac {
// Get the corresponding MC particle for Lc
auto arrayDaughtersLc = std::array{candidateLc.template prong0_as<aod::TracksWMc>(), candidateLc.template prong1_as<aod::TracksWMc>(), candidateLc.template prong2_as<aod::TracksWMc>()};
int8_t sign = 0;
int const indexMcLcRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughtersLc, o2::constants::physics::Pdg::kLambdaCPlus, std::array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
int const indexMcLcRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughtersLc, o2::constants::physics::Pdg::kLambdaCPlus, std::array{+kProton, -kKPlus, +kPiPlus}, true, &sign, 2);
auto particleLc = mcParticles.rawIteratorAt(indexMcLcRec);
// Get the corresponding MC particle for soft pion
auto particleSoftPi = candSc.prong1_as<aod::TracksWMc>().mcParticle();
Expand All @@ -1043,6 +1078,7 @@ struct HfTaskSigmac {
double cpaLc(candidateLc.cpa()), cpaXYLc(candidateLc.cpaXY());
int const origin = candSc.originMcRec();
auto channel = candidateLc.flagMcDecayChanRec(); /// 0: direct; 1: Λc± → p± K*; 2: Λc± → Δ(1232)±± K∓; 3: Λc± → Λ(1520) π±; FIXME: DecayChannelResonant
channel = remapResoChannelLc(channel);

/// candidate Λc+ → pK-π+ (and charge conjugate) within the range of M(pK-π+) chosen in the Σc0,++ builder
if ((TESTBIT(isCandPKPiPiKP, o2::aod::hf_cand_sigmac::Decays::PKPi)) && std::abs(candidateLc.template prong0_as<aod::TracksWMc>().mcParticle().pdgCode()) == kProton) {
Expand Down Expand Up @@ -1229,6 +1265,7 @@ struct HfTaskSigmac {
double cpaLc(candidateLc.cpa()), cpaXYLc(candidateLc.cpaXY());
int const origin = candidateLc.originMcRec();
auto channel = candidateLc.flagMcDecayChanRec(); /// 0: direct; 1: Λc± → p± K*; 2: Λc± → Δ(1232)±± K∓; 3: Λc± → Λ(1520) π±
channel = remapResoChannelLc(channel);
int pdgAbs = -1;
if (candidateLc.template prong0_as<aod::TracksWMc>().has_mcParticle()) {
pdgAbs = std::abs(candidateLc.template prong0_as<aod::TracksWMc>().mcParticle().pdgCode());
Expand Down
21 changes: 9 additions & 12 deletions PWGHF/TableProducer/candidateCreatorSigmac0plusplus.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -174,7 +174,7 @@
LOG(info) << "### ITS hitmap for soft pion";
LOG(info) << " >>> setSoftPiItsHitMap.size(): " << setSoftPiItsHitMap.size();
LOG(info) << " >>> Custom ITS hitmap dfchecked: ";
for (const auto it : setSoftPiItsHitMap) {

Check failure on line 177 in PWGHF/TableProducer/candidateCreatorSigmac0plusplus.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop] 

Use constant references for non-modified iterators in range-based for loops.
LOG(info) << " Layer " << static_cast<int>(it) << " ";
}
LOG(info) << "############";
Expand Down Expand Up @@ -516,20 +516,18 @@
/// 3. in case of (ii): resonant channel to pK-π+

/// look for Σc0(2455)
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaC0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
indexRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughters, Pdg::kSigmaC0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
flag = sign * o2::hf_decay::hf_cand_sigmac::DecayChannelMain::Sc0ToPKPiPi;
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaC0);
}
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaC0);

/// look for Σc0(2520)
if (flag == 0) {
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCStar0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
indexRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughters, Pdg::kSigmaCStar0, std::array{+kProton, -kKPlus, +kPiPlus, -kPiPlus}, true, &sign, 3);
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
flag = sign * o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStar0ToPKPiPi;
}
if (particleAntiparticle < 0) {
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCStar0);
}
Expand All @@ -543,20 +541,18 @@
/// 3. in case of (ii): resonant channel to pK-π+

/// look for Σc++(2455)
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
indexRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughters, Pdg::kSigmaCPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
flag = sign * o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScplusplusToPKPiPi;
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCPlusPlus);
}
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCPlusPlus);

/// look for Σc++(2520)
if (flag == 0) {
indexRec = RecoDecay::getMatchedMCRec(mcParticles, arrayDaughters, Pdg::kSigmaCStarPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
indexRec = RecoDecay::getMatchedMCRec<false, false, false, true, true>(mcParticles, arrayDaughters, Pdg::kSigmaCStarPlusPlus, std::array{+kProton, -kKPlus, +kPiPlus, +kPiPlus}, true, &sign, 3);
if (indexRec > -1) { /// due to (*) no need to check anything for LambdaC
flag = sign * o2::hf_decay::hf_cand_sigmac::DecayChannelMain::ScStarPlusPlusToPKPiPi;
}
if (particleAntiparticle < 0) {
auto particle = mcParticles.rawIteratorAt(indexRec);
particleAntiparticle = isParticleAntiparticle(particle, Pdg::kSigmaCStarPlusPlus);
}
Expand All @@ -581,6 +577,7 @@
for (const auto& particle : mcParticles) {
flag = 0;
origin = 0;
sign = 0;
std::vector<int> idxBhadMothers{};
int8_t particleAntiparticle = -1;

Expand Down
Loading