Skip to content

[PWGHF] Σc0,++: Remove duplicated consuming of McParticles. #11855

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 3 commits into from
Jul 3, 2025
Merged

[PWGHF] Σc0,++: Remove duplicated consuming of McParticles. #11855

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
6e3fe8e
Remove duplicated consuming of McParticles.
Jul 1, 2025
d8a1da7
Remove unneccessary comments.
Jul 1, 2025
557164f
Add comment about new data type.
Jul 2, 2025
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
18 changes: 9 additions & 9 deletions PWGHF/TableProducer/candidateCreatorSigmac0plusplus.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -417,7 +417,7 @@ struct HfCandidateSigmac0plusplusMc {

using BCsInfo = soa::Join<aod::BCs, aod::Timestamps, aod::BcSels>;
using LambdacMc = soa::Join<aod::HfCand3Prong, aod::HfSelLc, aod::HfCand3ProngMcRec>;
using McParticlesLcGenMatch = soa::Join<aod::McParticles, aod::HfCand3ProngMcGen>;
using McParticlesLcGenMatch = soa::Join<aod::McParticles, aod::HfCand3ProngMcGen>; // including response of particle matching to MC from candidate-creator-3-prong
using McCollisionsNoCents = soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels>;

PresliceUnsorted<McCollisionsNoCents> colPerMcCollision = aod::mccollisionlabel::mcCollisionId;
Expand Down Expand Up @@ -466,9 +466,9 @@ struct HfCandidateSigmac0plusplusMc {
/// @brief process function for MC matching of Σc0,++ → Λc+(→pK-π+) π- reconstructed candidates and counting of generated ones
/// @param candidatesSigmac reconstructed Σc0,++ candidates
/// @param mcParticles table of generated particles
void processMc(aod::McParticles const& mcParticles,
void processMc(McParticlesLcGenMatch const& mcParticles,
aod::TracksWMc const& tracks,
LambdacMc const& candsLc, McParticlesLcGenMatch const& mcParticlesLcGenMatch,
LambdacMc const& candsLc,
McCollisionsNoCents const& collInfos,
aod::McCollisions const&,
BCsInfo const&)
Expand Down Expand Up @@ -576,7 +576,7 @@ struct HfCandidateSigmac0plusplusMc {
} /// end loop over reconstructed Σc0,++ candidates

/// Match generated Σc0,++ candidates
for (const auto& particle : mcParticlesLcGenMatch) {
for (const auto& particle : mcParticles) {
flag = 0;
origin = 0;
std::vector<int> idxBhadMothers{};
Expand All @@ -602,7 +602,7 @@ struct HfCandidateSigmac0plusplusMc {
/// → here we check level 1. first, and then levels 2. and 3. are inherited by the Λc+ → pK-π+ MC matching in candidateCreator3Prong.cxx

/// look for Σc0,++(2455)
if (RecoDecay::isMatchedMCGen(mcParticlesLcGenMatch, particle, Pdg::kSigmaC0, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiMinus)}, true, &sign, 1)) {
if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kSigmaC0, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiMinus)}, true, &sign, 1)) {
// generated Σc0(2455)
for (const auto& daughter : particle.daughters_as<McParticlesLcGenMatch>()) {
// look for Λc+ daughter decaying in pK-π+
Expand All @@ -615,7 +615,7 @@ struct HfCandidateSigmac0plusplusMc {
break;
}
}
} else if (RecoDecay::isMatchedMCGen(mcParticlesLcGenMatch, particle, Pdg::kSigmaCPlusPlus, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiPlus)}, true, &sign, 1)) {
} else if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kSigmaCPlusPlus, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiPlus)}, true, &sign, 1)) {
// generated Σc++(2455)
for (const auto& daughter : particle.daughters_as<McParticlesLcGenMatch>()) {
// look for Λc+ daughter decaying in pK-π+
Expand All @@ -632,7 +632,7 @@ struct HfCandidateSigmac0plusplusMc {

/// look for Σc0,++(2520)
if (flag == 0) {
if (RecoDecay::isMatchedMCGen(mcParticlesLcGenMatch, particle, Pdg::kSigmaCStar0, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiMinus)}, true, &sign, 1)) {
if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kSigmaCStar0, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiMinus)}, true, &sign, 1)) {
// generated Σc0(2520)
for (const auto& daughter : particle.daughters_as<McParticlesLcGenMatch>()) {
// look for Λc+ daughter decaying in pK-π+
Expand All @@ -645,7 +645,7 @@ struct HfCandidateSigmac0plusplusMc {
break;
}
}
} else if (RecoDecay::isMatchedMCGen(mcParticlesLcGenMatch, particle, Pdg::kSigmaCStarPlusPlus, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiPlus)}, true, &sign, 1)) {
} else if (RecoDecay::isMatchedMCGen(mcParticles, particle, Pdg::kSigmaCStarPlusPlus, std::array{static_cast<int>(Pdg::kLambdaCPlus), static_cast<int>(kPiPlus)}, true, &sign, 1)) {
// generated Σc++(2520)
for (const auto& daughter : particle.daughters_as<McParticlesLcGenMatch>()) {
// look for Λc+ daughter decaying in pK-π+
Expand All @@ -663,7 +663,7 @@ struct HfCandidateSigmac0plusplusMc {

/// check the origin (prompt vs. non-prompt)
if (flag != 0) {
origin = RecoDecay::getCharmHadronOrigin(mcParticlesLcGenMatch, particle, false, &idxBhadMothers);
origin = RecoDecay::getCharmHadronOrigin(mcParticles, particle, false, &idxBhadMothers);
}
/// fill the table with results of generation level MC matching
if (origin == RecoDecay::OriginType::NonPrompt) {
Expand Down
Loading