Skip to content

[PWGCF] FemtoUniverse: Correct getting PDG of particle's mother #15411

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 16, 2026
+14 −10
Merged

[PWGCF] FemtoUniverse: Correct getting PDG of particle's mother #15411

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
16 changes: 9 additions & 7 deletions PWGCF/FemtoUniverse/TableProducer/femtoUniverseProducerTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -881,12 +881,14 @@ struct FemtoUniverseProducerTask {
template <typename ParticleType>
int32_t getMotherPDG(ParticleType particle)
{
auto motherparticlesMC = particle.template mothers_as<aod::McParticles>();
if (!motherparticlesMC.empty()) {
if (particle.isPhysicalPrimary()) {
return 0;
} else if (particle.has_mothers()) {
auto motherparticlesMC = particle.template mothers_as<aod::McParticles>();
auto motherparticleMC = motherparticlesMC.front();
return particle.isPhysicalPrimary() ? 0 : motherparticleMC.pdgCode();
return motherparticleMC.pdgCode();
} else {
return 9999;
return 999;
}
}

Expand Down Expand Up @@ -926,7 +928,7 @@ struct FemtoUniverseProducerTask {
outputPartsMCLabels(outputPartsMC.lastIndex());
} else {
outputPartsMCLabels(-1);
outputDebugPartsMC(9999);
outputDebugPartsMC(-999);
}
}

Expand Down Expand Up @@ -2169,7 +2171,7 @@ struct FemtoUniverseProducerTask {
// aligned, so that they can be joined in the task.
if constexpr (transientLabels) {
outputPartsMCLabels(-1);
outputDebugPartsMC(9999);
outputDebugPartsMC(-999);
}
}
if constexpr (resolveDaughs) {
Expand Down Expand Up @@ -2214,7 +2216,7 @@ struct FemtoUniverseProducerTask {
// aligned, so that they can be joined in the task.
if constexpr (transientLabels) {
outputPartsMCLabels(-1);
outputDebugPartsMC(9999);
outputDebugPartsMC(-999);
}
}
}
Expand Down
8 changes: 5 additions & 3 deletions PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Extended.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -161,6 +161,7 @@ struct FemtoUniversePairTaskTrackV0Extended {
Configurable<int> confEtaBins{"confEtaBins", 29, "Number of eta bins in deta dphi"};
ConfigurableAxis confmTBins3D{"confmTBins3D", {VARIABLE_WIDTH, 1.02f, 1.14f, 1.20f, 1.26f, 1.38f, 1.56f, 1.86f, 4.50f}, "mT Binning for the 3Dimensional plot: k* vs multiplicity vs mT (set <<confUse3D>> to true in order to use)"};
ConfigurableAxis confMultBins3D{"confMultBins3D", {VARIABLE_WIDTH, 0.0f, 20.0f, 30.0f, 40.0f, 99999.0f}, "multiplicity Binning for the 3Dimensional plot: k* vs multiplicity vs mT (set <<confUse3D>> to true in order to use)"};
ConfigurableAxis confMotherPDGBins{"confMotherPDGBins", {8001, -4000, 4000}, "Binning for the mothers' PDG code in pair fractions histogram"};

struct : o2::framework::ConfigurableGroup {
Configurable<bool> confIsCPR{"confIsCPR", true, "Close Pair Rejection"};
Expand Down Expand Up @@ -328,7 +329,7 @@ struct FemtoUniversePairTaskTrackV0Extended {
registryMCtruth.add("minus/MCtruthPrPt", "MC truth protons;#it{p}_{T} (GeV/c)", {HistType::kTH1F, {{500, 0, 5}}});

if (doprocessPairFractionsMCTruthV0 || doprocessPairFractionsMCTruth) {
registryMCtruth.add("mothersTruth/motherParticle", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
registryMCtruth.add("mothersTruth/motherParticle", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {confMotherPDGBins, confMotherPDGBins}});
}

// MC reco
Expand All @@ -354,8 +355,8 @@ struct FemtoUniversePairTaskTrackV0Extended {
registryMCreco.add("minus/MCrecoPrPt", "MC reco protons;#it{p}_{T} (GeV/c)", {HistType::kTH1F, {{500, 0, 5}}});

if (doprocessPairFractions || doprocessPairFractionsV0) {
registryMCreco.add("mothersReco/motherParticle", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
registryMCreco.add("mothersReco/motherParticlePDGCheck", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {{8001, -4000, 4000}, {8001, -4000, 4000}}});
registryMCreco.add("mothersReco/motherParticle", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {confMotherPDGBins, confMotherPDGBins}});
registryMCreco.add("mothersReco/motherParticlePDGCheck", "pair fractions;part1 mother PDG;part2 mother PDG", {HistType::kTH2F, {confMotherPDGBins, confMotherPDGBins}});
}
sameEventCont.init(&resultRegistry, confkstarBins, confMultBins, confkTBins, confmTBins, confMultBins3D, confmTBins3D, confEtaBins, confPhiBins, confIsMC, confUse3D);
sameEventCont.setPDGCodes(ConfTrkSelection.confTrkPDGCodePartOne, ConfV0Selection.confV0PDGCodePartTwo);
Expand Down Expand Up @@ -1450,6 +1451,7 @@ struct FemtoUniversePairTaskTrackV0Extended {
continue;
if ((ConfV0Selection.confV0Type2 == 0 && mcParticle2.pdgMCTruth() != kLambda0) || (ConfV0Selection.confV0Type2 == 1 && mcParticle2.pdgMCTruth() != kLambda0Bar))
continue;

registryMCreco.fill(HIST("mothersReco/motherParticlePDGCheck"), p1.motherPDG(), p2.motherPDG());
}
};
Expand Down
Loading