Skip to content

[PWGCF] FemtoUniverse: Check process that generated MC particle when processing fractions #15503

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
prchakra merged 1 commit into from
Mar 24, 2026
+16 −2
Merged

[PWGCF] FemtoUniverse: Check process that generated MC particle when processing fractions #15503

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
9 changes: 7 additions & 2 deletions PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -884,6 +884,9 @@ struct FemtoUniverseProducerTask {
if (particle.isPhysicalPrimary()) {
return 0;
} else if (particle.has_mothers()) {
if (particle.getProcess() == 20 || particle.getProcess() == 23) { // treat particles from hadronic scattering (20, 23) as primary
return 0;
}
auto motherparticlesMC = particle.template mothers_as<aod::McParticles>();
auto motherparticleMC = motherparticlesMC.front();
return motherparticleMC.pdgCode();
Expand Down Expand Up @@ -2142,6 +2145,7 @@ struct FemtoUniverseProducerTask {
// instead of the bitmask, the PDG of the particle is stored as uint32_t

int32_t variablePDG = confStoreMCmothers ? getMotherPDG(particle) : particle.pdgCode();
int32_t variableCut = confStoreMCmothers ? particle.getProcess() : 0;

// now the table is filled
if constexpr (resolveDaughs) {
Expand All @@ -2156,7 +2160,7 @@ struct FemtoUniverseProducerTask {
particle.eta(),
particle.phi(),
aod::femtouniverseparticle::ParticleType::kMCTruthTrack,
0,
variableCut,
pdgCode,
variablePDG,
childIDs,
Expand Down Expand Up @@ -2195,13 +2199,14 @@ struct FemtoUniverseProducerTask {
}

int32_t variablePDG = confStoreMCmothers ? getMotherPDG(particle) : particle.pdgCode();
int32_t variableCut = confStoreMCmothers ? particle.getProcess() : 0;

outputParts(outputCollision.lastIndex(),
particle.pt(),
particle.eta(),
particle.phi(),
aod::femtouniverseparticle::ParticleType::kMCTruthTrack,
0,
variableCut,
static_cast<uint32_t>(particle.pdgCode()),
variablePDG,
childIDs,
Expand Down
6 changes: 6 additions & 0 deletions PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackTrackMultKtExtended.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -366,6 +366,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
}
if (doprocessFractionsMCTruth) {
mixedMultRegistryPM.add("MCtruth/motherParticle", "pair fractions; part1 mother PDG; part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
mixedMultRegistryPM.add("MCtruth/mcProcess", "pair ; part1 VMC physics code; part2 VMC physics code", {HistType::kTH2F, {{50, 0, 50}, {50, 0, 50}}});
}
}

Expand All @@ -376,6 +377,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
}
if (doprocessFractionsMCTruth) {
mixedMultRegistryPP.add("MCtruth/motherParticle", "pair fractions; part1 mother PDG; part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
mixedMultRegistryPP.add("MCtruth/mcProcess", "pair ; part1 VMC physics code; part2 VMC physics code", {HistType::kTH2F, {{50, 0, 50}, {50, 0, 50}}});
}
sameEventContPP.init(&resultRegistryPP, confkstarBins, confMultBins, confkTBins, confmTBins, confmultBins3D, confmTBins3D, twotracksconfigs.confIsMC, twotracksconfigs.confUse3D);
mixedEventContPP.init(&resultRegistryPP, confkstarBins, confMultBins, confkTBins, confmTBins, confmultBins3D, confmTBins3D, twotracksconfigs.confIsMC, twotracksconfigs.confUse3D);
Expand All @@ -395,6 +397,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
}
if (doprocessFractionsMCTruth) {
mixedMultRegistryMM.add("MCtruth/motherParticle", "pair fractions; part1 mother PDG; part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
mixedMultRegistryMM.add("MCtruth/mcProcess", "pair ; part1 VMC physics code; part2 VMC physics code", {HistType::kTH2F, {{50, 0, 50}, {50, 0, 50}}});
}
sameEventContMM.init(&resultRegistryMM, confkstarBins, confMultBins, confkTBins, confmTBins, confmultBins3D, confmTBins3D, twotracksconfigs.confIsMC, twotracksconfigs.confUse3D);
mixedEventContMM.init(&resultRegistryMM, confkstarBins, confMultBins, confkTBins, confmTBins, confmultBins3D, confmTBins3D, twotracksconfigs.confIsMC, twotracksconfigs.confUse3D);
Expand Down Expand Up @@ -955,6 +958,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
const auto& pdgParticleTwo = pdg->GetParticle(pdgCodePartTwo);
if (pdgParticleOne && pdgParticleTwo && (pdgCodePartOne == trackonefilter.confPDGCodePartOne) && (pdgCodePartTwo == tracktwofilter.confPDGCodePartTwo)) {
mixedMultRegistryPM.fill(HIST("MCtruth/motherParticle"), p1.tempFitVar(), p2.tempFitVar());
mixedMultRegistryPM.fill(HIST("MCtruth/mcProcess"), p1.cut(), p2.cut());
}
}
break;
Expand All @@ -967,6 +971,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
const auto& pdgParticleTwo = pdg->GetParticle(pdgCodePartTwo);
if (pdgParticleOne && pdgParticleTwo && (pdgCodePartOne == trackonefilter.confPDGCodePartOne) && (pdgCodePartTwo == trackonefilter.confPDGCodePartOne)) {
mixedMultRegistryPP.fill(HIST("MCtruth/motherParticle"), p1.tempFitVar(), p2.tempFitVar());
mixedMultRegistryPP.fill(HIST("MCtruth/mcProcess"), p1.cut(), p2.cut());
}
}
break;
Expand All @@ -979,6 +984,7 @@ struct FemtoUniversePairTaskTrackTrackMultKtExtended {
const auto& pdgParticleTwo = pdg->GetParticle(pdgCodePartTwo);
if (pdgParticleOne && pdgParticleTwo && (pdgCodePartOne == tracktwofilter.confPDGCodePartTwo) && (pdgCodePartTwo == tracktwofilter.confPDGCodePartTwo)) {
mixedMultRegistryMM.fill(HIST("MCtruth/motherParticle"), p1.tempFitVar(), p2.tempFitVar());
mixedMultRegistryMM.fill(HIST("MCtruth/mcProcess"), p1.cut(), p2.cut());
}
}
break;
Expand Down
3 changes: 3 additions & 0 deletions PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Extended.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -336,6 +336,7 @@ struct FemtoUniversePairTaskTrackV0Extended {

if (doprocessPairFractionsMCTruthV0 || doprocessPairFractionsMCTruth) {
registryMCtruth.add("mothersTruth/motherParticle", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {confMotherPDGBins, confMotherPDGBins}});
registryMCtruth.add("mothersTruth/mcProcess", "pair ; part1 VMC physics code; part2 VMC physics code", {HistType::kTH2F, {{50, 0, 50}, {50, 0, 50}}});
}

// MC reco
Expand Down Expand Up @@ -1547,6 +1548,7 @@ struct FemtoUniversePairTaskTrackV0Extended {
continue;

registryMCtruth.fill(HIST("mothersTruth/motherParticle"), p1.tempFitVar(), p2.tempFitVar());
registryMCtruth.fill(HIST("mothersTruth/mcProcess"), p1.cut(), p2.cut());
}
};

Expand Down Expand Up @@ -1580,6 +1582,7 @@ struct FemtoUniversePairTaskTrackV0Extended {
continue;

registryMCtruth.fill(HIST("mothersTruth/motherParticle"), p1.tempFitVar(), p2.tempFitVar());
registryMCtruth.fill(HIST("mothersTruth/mcProcess"), p1.cut(), p2.cut());
}
};

Expand Down
Loading