Skip to content

Commit 1be094f

Browse files
committed
FU: adding colCheck, mother PDG check, fixing a bug in PurityQA process
1 parent fbbe848 commit 1be094f

3 files changed

Lines changed: 113 additions & 149 deletions

File tree

PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -2886,6 +2886,8 @@ struct FemtoUniverseProducerTask {
28862886
aod::BCsWithTimestamps const&,
28872887
aod::McParticles const& mcParts)
28882888
{
2889+
std::set<int> mcColIds;
2890+
mcColIds.clear();
28892891
// MC Reco
28902892
for (const auto& col : collisions) {
28912893
auto groupedTracks = tracks.sliceBy(perCollisionTracks, col.globalIndex());
@@ -2895,12 +2897,16 @@ struct FemtoUniverseProducerTask {
28952897
// fill the tables
28962898
const auto colcheck = fillCollisions<true>(col, tracks);
28972899
if (colcheck) {
2900+
mcColIds.insert(col.mcCollisionId());
28982901
fillTracks<true>(groupedTracks);
28992902
fillD0D0barMcMl<true>(col, groupedTracks, groupedD0s, mcParts);
29002903
}
29012904
}
29022905
// MC Truth
29032906
for (const auto& mccol : mccols) {
2907+
if (confCollMCTruthOnlyReco && !mcColIds.contains(mccol.globalIndex())) {
2908+
continue;
2909+
}
29042910
auto groupedMCParticles = hfMcGenCands.sliceBy(mcPartPerMcColl, mccol.globalIndex());
29052911
auto groupedCollisions = collisions.sliceBy(recoCollsPerMCColl, mccol.globalIndex());
29062912
fillMCTruthCollisions(groupedCollisions, groupedMCParticles); // fills the reco collisions for mc collision

0 commit comments

Comments
 (0)