Skip to content

[PWGCF] FemtoUniverse -- Add hInvMassK0Short to helicity analysis. #16102

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 2 commits into from
May 4, 2026
+30 −19
Merged

[PWGCF] FemtoUniverse -- Add hInvMassK0Short to helicity analysis. #16102

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
9cf21c8
Add hInvMassK0Short to helicity analysis
a-m-andrushko May 4, 2026
5ed8a8b
Merge branch 'AliceO2Group:master' into hInvMassK0Short
a-m-andrushko May 4, 2026
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
49 changes: 30 additions & 19 deletions PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,7 @@

#include <TFile.h>
#include <TH1.h>
#include <TPDGCode.h>

#include <cmath>
#include <cstddef>
Expand Down Expand Up @@ -188,6 +189,7 @@
Configurable<bool> cfgProcessHel4{"cfgProcessHel4", false, "Process particle pairs from the helicity range 4"}; // -0.5 > cosineTheta >= -1.0
ConfigurableAxis confInvMassMotherpTBinsHel{"confInvMassMotherpTBinsHel", {5, 0, 5}, "pT binning in the pT vs. InvMassMother plot for helicity"};
ConfigurableAxis confInvMassMotherBinsHel{"confInvMassMotherBinsHel", {1000, 0.8, 1.4}, "InvMassMother binning in the pT vs. InvMassMother plot for helicity"};
ConfigurableAxis confInvMassK0Short{"confInvMassK0s", {1000, 0.2, 0.8}, "Invariant mass binning for K0 Short"};

/// Efficiency
Configurable<std::string> confLocalEfficiency{"confLocalEfficiency", "", "Local path to efficiency .root file"};
Expand Down Expand Up @@ -321,6 +323,7 @@
thetaRegistry.add("Theta/Mother/hInvMassMotherHel2", " ; p_{T} (GeV/#it{c}); M_{#Lambda};", kTH2F, {confInvMassMotherpTBinsHel, confInvMassMotherBinsHel});
thetaRegistry.add("Theta/Mother/hInvMassMotherHel3", " ; p_{T} (GeV/#it{c}); M_{#Lambda};", kTH2F, {confInvMassMotherpTBinsHel, confInvMassMotherBinsHel});
thetaRegistry.add("Theta/Mother/hInvMassMotherHel4", " ; p_{T} (GeV/#it{c}); M_{#Lambda};", kTH2F, {confInvMassMotherpTBinsHel, confInvMassMotherBinsHel});
thetaRegistry.add("Theta/Mother/hInvMassK0Short", " ; M_{K^{0}_{S}}; ;", kTH1F, {confInvMassK0Short});

/// MC Truth
registryMCtruth.add("plus/MCtruthLambda", "MC truth Lambdas;#it{p}_{T} (GeV/c); #eta", {HistType::kTH2F, {{500, 0, 5}, {400, -1.0, 1.0}}});
Expand Down Expand Up @@ -450,7 +453,7 @@
}

/// This function processes the same event for Track-V0
template <typename PartType, typename PartitionType, typename MCParticles = std::nullptr_t>
template <bool confIsMC, typename PartType, typename PartitionType, typename MCParticles = std::nullptr_t>
void doSameEvent(FilteredFDCollision const& col, PartType const& parts, PartitionType& groupPartsOne, PartitionType& groupPartsTwo, int helRange, [[maybe_unused]] MCParticles mcParts = nullptr)
{
const auto& magFieldTesla = col.magField();
Expand Down Expand Up @@ -488,14 +491,22 @@
thetaRegistry.fill(HIST("Theta/NegativeChild/hThetaEta"), negChild.eta(), cosineTheta);
thetaRegistry.fill(HIST("Theta/NegativeChild/hThetaPhi"), negChild.phi(), cosineTheta);

if (cosineTheta <= 1.0 && cosineTheta >= 0.1)

Check failure on line 494 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
thetaRegistry.fill(HIST("Theta/Mother/hInvMassMotherHel1"), part.pt(), part.mLambda());
else if (cosineTheta < 0.1 && cosineTheta >= -0.1)

Check failure on line 496 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
thetaRegistry.fill(HIST("Theta/Mother/hInvMassMotherHel2"), part.pt(), part.mLambda());
else if (cosineTheta < -0.1 && cosineTheta >= -0.5)

Check failure on line 498 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
thetaRegistry.fill(HIST("Theta/Mother/hInvMassMotherHel3"), part.pt(), part.mLambda());
else if (cosineTheta < -0.5 && cosineTheta >= -1)

Check failure on line 500 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
thetaRegistry.fill(HIST("Theta/Mother/hInvMassMotherHel4"), part.pt(), part.mLambda());

if constexpr (confIsMC) {
if (part.has_fdMCParticle()) {
if ((part.fdMCParticle()).pdgMCTruth() == kK0Short) {
thetaRegistry.fill(HIST("Theta/Mother/hInvMassK0Short"), part.mKaon());
}
}
}
}

for (const auto& part : groupPartsOne) {
Expand Down Expand Up @@ -560,28 +571,28 @@
}

case 1: {
if (cosineTheta <= 1.0 && cosineTheta >= 0.1)

Check failure on line 574 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
sameEventContHel1.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
}

case 2: {
if (cosineTheta < 0.1 && cosineTheta >= -0.1)

Check failure on line 581 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
sameEventContHel2.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
}

case 3: {
if (cosineTheta < -0.1 && cosineTheta >= -0.5)

Check failure on line 588 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
sameEventContHel3.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
}

case 4: {
if (cosineTheta < -0.5 && cosineTheta >= -1.0)

Check failure on line 595 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
sameEventContHel4.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
Expand All @@ -599,27 +610,27 @@
auto groupPartsTwo = partsTwo->sliceByCached(aod::femtouniverseparticle::fdCollisionId, col.globalIndex(), cache);

if (cfgProcessHel)
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 0);
doSameEvent<false>(col, parts, groupPartsOne, groupPartsTwo, 0);

if (cfgProcessHel1)
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 1);
doSameEvent<false>(col, parts, groupPartsOne, groupPartsTwo, 1);

if (cfgProcessHel2)
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 2);
doSameEvent<false>(col, parts, groupPartsOne, groupPartsTwo, 2);

if (cfgProcessHel3)
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 3);
doSameEvent<false>(col, parts, groupPartsOne, groupPartsTwo, 3);

if (cfgProcessHel4)
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 4);
doSameEvent<false>(col, parts, groupPartsOne, groupPartsTwo, 4);
}
PROCESS_SWITCH(FemtoUniversePairTaskTrackV0Helicity, processSameEvent, "Enable processing same event for track - V0", false);

void processSameEventMCReco(FilteredFDCollision const& col, FemtoRecoParticles const& parts, aod::FdMCParticles const& mcparts)
{
auto groupPartsOne = partsOneMCReco->sliceByCached(aod::femtouniverseparticle::fdCollisionId, col.globalIndex(), cache);
auto groupPartsTwo = partsTwoMCReco->sliceByCached(aod::femtouniverseparticle::fdCollisionId, col.globalIndex(), cache);
doSameEvent(col, parts, groupPartsOne, groupPartsTwo, 0, mcparts);
doSameEvent<true>(col, parts, groupPartsOne, groupPartsTwo, 0, mcparts);
}
PROCESS_SWITCH(FemtoUniversePairTaskTrackV0Helicity, processSameEventMCReco, "Enable processing same event for track - V0 MC Reco", false);

Expand Down Expand Up @@ -858,14 +869,14 @@
}

case 1: {
if (cosineTheta <= 1.0 && cosineTheta >= 0.1)

Check failure on line 872 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
mixedEventContHel1.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
}

case 2: {
if (cosineTheta < 0.1 && cosineTheta >= -0.1)

Check failure on line 879 in PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackV0Helicity.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[magic-number] 

Avoid magic numbers in expressions. Assign the value to a clearly named variable or constant.
mixedEventContHel2.setPair<false>(p1, p2, multCol, confUse3D, weight);

break;
Expand Down Expand Up @@ -1099,7 +1110,7 @@
if (!pdgParticle)
continue;

if (pdgCode == 3122) {
if (pdgCode == kLambda0) {
registryMCtruth.fill(HIST("plus/MCtruthLambda"), part.pt(), part.eta());

// Helicity angle
Expand All @@ -1125,31 +1136,31 @@
registryMCtruth.fill(HIST("ThetaMCTruth/NegativeChild/hThetaPhi"), negChild.phi(), cosineTheta);

continue;
} else if (pdgCode == -3122) {
} else if (pdgCode == -kLambda0) {
registryMCtruth.fill(HIST("minus/MCtruthLambda"), part.pt(), part.eta());
continue;
}

if (pdgParticle->Charge() > 0.0) {
registryMCtruth.fill(HIST("plus/MCtruthAllPt"), part.pt());
}
if (pdgCode == 211) {
if (pdgCode == kPiPlus) {
registryMCtruth.fill(HIST("plus/MCtruthPi"), part.pt(), part.eta());
registryMCtruth.fill(HIST("plus/MCtruthPiPt"), part.pt());
}
if (pdgCode == 2212) {
if (pdgCode == kProton) {
registryMCtruth.fill(HIST("plus/MCtruthPr"), part.pt(), part.eta());
registryMCtruth.fill(HIST("plus/MCtruthPrPt"), part.pt());
}

if (pdgParticle->Charge() < 0.0) {
registryMCtruth.fill(HIST("minus/MCtruthAllPt"), part.pt());
}
if (pdgCode == -211) {
if (pdgCode == kPiMinus) {
registryMCtruth.fill(HIST("minus/MCtruthPi"), part.pt(), part.eta());
registryMCtruth.fill(HIST("minus/MCtruthPiPt"), part.pt());
}
if (pdgCode == -2212) {
if (pdgCode == -kProton) {
registryMCtruth.fill(HIST("minus/MCtruthPr"), part.pt(), part.eta());
registryMCtruth.fill(HIST("minus/MCtruthPrPt"), part.pt());
}
Expand All @@ -1166,7 +1177,7 @@
const auto& mcpart = mcparts.iteratorAt(mcPartId);
//
if (part.partType() == aod::femtouniverseparticle::ParticleType::kV0) {
if (mcpart.pdgMCTruth() == 3122) {
if (mcpart.pdgMCTruth() == kLambda0) {
const auto& posChild = parts.iteratorAt(part.globalIndex() - 2);
const auto& negChild = parts.iteratorAt(part.globalIndex() - 1);
/// Daughters that do not pass this condition are not selected
Expand All @@ -1181,7 +1192,7 @@
registryMCreco.fill(HIST("plus/MCrecoLambdaChildPi"), mcpartChild.pt(), mcpartChild.eta()); // lambda pion child
}
}
} else if (mcpart.pdgMCTruth() == -3122) {
} else if (mcpart.pdgMCTruth() == -kLambda0) {
const auto& posChild = parts.iteratorAt(part.globalIndex() - 2);
const auto& negChild = parts.iteratorAt(part.globalIndex() - 1);
/// Daughters that do not pass this condition are not selected
Expand All @@ -1200,21 +1211,21 @@
} else if (part.partType() == aod::femtouniverseparticle::ParticleType::kTrack) {
if (part.sign() > 0) {
registryMCreco.fill(HIST("plus/MCrecoAllPt"), mcpart.pt());
if (mcpart.pdgMCTruth() == 211 && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePi()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePi()))) {
if (mcpart.pdgMCTruth() == kPiPlus && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePi()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePi()))) {
registryMCreco.fill(HIST("plus/MCrecoPi"), mcpart.pt(), mcpart.eta());
registryMCreco.fill(HIST("plus/MCrecoPiPt"), mcpart.pt());
} else if (mcpart.pdgMCTruth() == 2212 && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePr()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePr()))) {
} else if (mcpart.pdgMCTruth() == kProton && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePr()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePr()))) {
registryMCreco.fill(HIST("plus/MCrecoPr"), mcpart.pt(), mcpart.eta());
registryMCreco.fill(HIST("plus/MCrecoPrPt"), mcpart.pt());
}
}

if (part.sign() < 0) {
registryMCreco.fill(HIST("minus/MCrecoAllPt"), mcpart.pt());
if (mcpart.pdgMCTruth() == -211 && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePi()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePi()))) {
if (mcpart.pdgMCTruth() == kPiMinus && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePi()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePi()))) {
registryMCreco.fill(HIST("minus/MCrecoPi"), mcpart.pt(), mcpart.eta());
registryMCreco.fill(HIST("minus/MCrecoPiPt"), mcpart.pt());
} else if (mcpart.pdgMCTruth() == -2212 && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePr()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePr()))) {
} else if (mcpart.pdgMCTruth() == -kProton && isNSigmaCombined(part.p(), aod::pidtpc_tiny::binning::unPackInTable(part.tpcNSigmaStorePr()), aod::pidtof_tiny::binning::unPackInTable(part.tofNSigmaStorePr()))) {
registryMCreco.fill(HIST("minus/MCrecoPr"), mcpart.pt(), mcpart.eta());
registryMCreco.fill(HIST("minus/MCrecoPrPt"), mcpart.pt());
}
Expand Down
Loading