Skip to content

[PWGHF] Add possibility to run more loops on MC particles in dataCreatorCharmHadPiReduced #18067

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 1 commit into from
Sep 26, 2026
+18 −0
Merged

[PWGHF] Add possibility to run more loops on MC particles in dataCreatorCharmHadPiReduced #18067

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
18 changes: 18 additions & 0 deletions PWGHF/D2H/TableProducer/dataCreatorCharmHadPiReduced.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -177,6 +177,7 @@ struct HfDataCreatorCharmHadPiReduced {
Configurable<double> invMassWindowCharmHadPi{"invMassWindowCharmHadPi", 0.3, "invariant-mass window for CharmHad-Pi pair preselections (GeV/c2)"};
// MC extra
Configurable<bool> checkDecayTypeMc{"checkDecayTypeMc", false, "flag to enable MC checks on decay type"};
Configurable<int> nMcLoops{"nMcLoops", 0, "number of extra loops over generated particles per MC collision"};
} configs;
// vertexing
struct : o2::framework::ConfigurableGroup {
Expand Down Expand Up @@ -394,6 +395,9 @@ struct HfDataCreatorCharmHadPiReduced {
hCandidatesLc = registry.add<TH1>("hCandidatesLc", "Lc candidate counter", {HistType::kTH1D, {axisCands}});
hCandidatesD0FromDstar = registry.add<TH1>("hCandidatesD0FromDstar", "D0 from D* candidate counter", {HistType::kTH1D, {axisCands}});
hCandidatesBHadron = registry.add<TH1>("hCandidatesBHadron", "B hadron candidate counter", {HistType::kTH1D, {axisCands}});
if (configs.nMcLoops > 0) {
registry.add("hNMatchedMcLoops", "matches found in extra MC loops;matches per MC collision;entries", {HistType::kTH1D, {{100, -0.5, 99.5}}});
}

setLabelHistoCands(hCandidatesD0);
setLabelHistoCands(hCandidatesDPlus);
Expand Down Expand Up @@ -1793,6 +1797,20 @@ struct HfDataCreatorCharmHadPiReduced {
ptProngs[1], yProngs[1], etaProngs[1], hfRejMap, centFT0C, centFT0M);
}
} // gen

// extra loops over generated particles (increases CPU time)
if (configs.nMcLoops > 0) {
int nMatched{0};
for (int iLoop = 0; iLoop < configs.nMcLoops; ++iLoop) {
for (const auto& particle : mcParticlesPerMcColl) {
nMatched += RecoDecay::isMatchedMCGen<true>(particlesMc, particle, Pdg::kB0, std::array{-static_cast<int>(Pdg::kDPlus), +kPiPlus}, true);
nMatched += RecoDecay::isMatchedMCGen<true>(particlesMc, particle, Pdg::kBPlus, std::array{-static_cast<int>(Pdg::kDPlus), +kPiPlus}, true);
nMatched += RecoDecay::isMatchedMCGen<true>(particlesMc, particle, Pdg::kBS, std::array{-static_cast<int>(Pdg::kDPlus), +kPiPlus}, true);
nMatched += RecoDecay::isMatchedMCGen<true>(particlesMc, particle, Pdg::kLambdaB0, std::array{-static_cast<int>(Pdg::kDPlus), +kPiPlus}, true);
}
}
registry.fill(HIST("hNMatchedMcLoops"), nMatched);
}
}

////////////////////////////////////////////////////////////////////////////////////////////////////
Expand Down
Loading