Skip to content

[PWGLF] Update DCA cuts for systematic studies and add V0 mass plots #16120

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

Open
PaolaVT wants to merge 3 commits into
base: master
Choose a base branch
from
+127 −31
Open

[PWGLF] Update DCA cuts for systematic studies and add V0 mass plots #16120

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
1399f53
New hist was added
Apr 24, 2026
83e925b
Merge branch 'AliceO2Group:master' into master
PaolaVT May 5, 2026
7e25b0d
DCA cuts were changed, mass plots were added
May 5, 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
158 changes: 127 additions & 31 deletions PWGLF/Tasks/Nuspex/dedxPidAnalysis.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -83,7 +83,6 @@ struct DedxPidAnalysis {
static constexpr int EtaIntervals = 8;
static constexpr int ParticlesType = 4;
static constexpr int CentralityClasses = 10;
float tpcCut = 0.6;
float pionMin = 0.35;
float pionMax = 0.45;
float elTofCut = 0.1;
Expand Down Expand Up @@ -197,6 +196,11 @@ struct DedxPidAnalysis {
Configurable<float> etaMin{"etaMin", -0.8f, "etaMin"};
Configurable<float> etaMax{"etaMax", +0.8f, "etaMax"};
Configurable<float> minNCrossedRowsOverFindableClustersTPC{"minNCrossedRowsOverFindableClustersTPC", 0.8f, "Additional cut on the minimum value of the ratio between crossed rows and findable clusters in the TPC"};
Configurable<float> nSigmaDCAxy{"nSigmaDCAxy", 3.0, "nSigma DCAxy selection"};
Configurable<float> dcaXYp0{"dcaXYp0", 0.0105f, "DCAxy formula: p0 + p1/pt^p2"};
Configurable<float> dcaXYp1{"dcaXYp1", 0.0350f, "DCAxy p1 parameter"};
Configurable<float> dcaXYp2{"dcaXYp2", 1.1f, "DCA_xy p2 parameter"};
Configurable<float> nSigmaDCAz{"nSigmaDCAz", 3.0, "nSigma DCAz selection"};
Configurable<float> maxDCAz{"maxDCAz", 0.1f, "maxDCAz"};
// v0 cuts
Configurable<float> v0cospaMin{"v0cospaMin", 0.999f, "Minimum V0 CosPA"};
Expand All @@ -215,7 +219,8 @@ struct DedxPidAnalysis {
Configurable<float> v0rapidityCut{"v0rapidityCut", 0.5f, "V0 rapidity cut"};
Configurable<float> v0ProperLifetimeCutK0s{"v0ProperLifetimeCutK0s", 20.f, "V0 proper lifetime cut for K0s"};
Configurable<float> v0ProperLifetimeCutLambda{"v0ProperLifetimeCutLambda", 30.f, "V0 proper lifetime cut for Lambda"};
Configurable<float> nsigmaTOFmax{"nsigmaTOFmax", 3.0f, "Maximum nsigma TOF"};
Configurable<float> nsigmaMax{"nsigmaMax", 3.0f, "Maximum nsigma cut"};
Configurable<float> tpcMomentumCut{"tpcMomentumCut", 0.6f, "Momentum threshold above which TOF is required"};
Configurable<float> invMassCutK0s{"invMassCutK0s", 0.015f, "invariant Mass Cut for K0s"};
Configurable<float> invMassCutLambda{"invMassCutLambda", 0.015f, "invariant Mass Cut for Lambda"};
Configurable<float> invMassCutGamma{"invMassCutGamma", 0.015f, "invariant Mass Cut for Gamma"};
Expand Down Expand Up @@ -276,8 +281,8 @@ struct DedxPidAnalysis {
selectedTracks.SetMaxChi2PerClusterTPC(maxChi2TPC);
selectedTracks.SetRequireHitsInITSLayers(1, {0, 1, 2});
selectedTracks.SetMaxChi2PerClusterITS(maxChi2ITS);
selectedTracks.SetMaxDcaXYPtDep([](float pt) { return 0.0105f + 0.0350f / std::pow(pt, 1.1f); });
selectedTracks.SetMaxDcaZ(maxDCAz);
// selectedTracks.SetMaxDcaXYPtDep([](float pt) { return 0.0105f + 0.0350f / std::pow(pt, 1.1f); });
// selectedTracks.SetMaxDcaZ(maxDCAz);
selectedTracks.SetRequireGoldenChi2(true);

return selectedTracks;
Expand All @@ -299,7 +304,8 @@ struct DedxPidAnalysis {
LOGF(info, "GoodZvtxFT0vsPV cut disabled");
}
if (nINELSelectionMode == NoSelINEL) {
LOGF(info, "INEL cut disabled");
LOGF(info, "Applying just INEL");
label = "INEL";
} else if (nINELSelectionMode == SelINELgt0) {
LOGF(info, "Applying INEL > 0 cut");
label = "INEL > 0";
Expand Down Expand Up @@ -483,6 +489,30 @@ struct DedxPidAnalysis {
for (int i = 0; i < CentralityClasses + 1; ++i) {
hDedxVspTMomentumVsCent[i] = registryDeDx.add<TH3>(DedxvspTMomentumvsCent[i].data(), "dE/dx", HistType::kTH3F, {{ptAxis}, {dedxAxis}, {etaAxis}});
}

// Invariant Mass
registryDeDx.add("hMassK0s", "Invariant mass K0s;m_{#pi#pi} (GeV/c^{2});Counts",
HistType::kTH1F, {{200, 0.4, 0.6}});
registryDeDx.add("hMassLambda", "Invariant mass #Lambda;m_{p#pi} (GeV/c^{2});Counts",
HistType::kTH1F, {{200, 1.08, 1.16}});
registryDeDx.add("hMassAntiLambda", "Invariant mass #bar{#Lambda};m_{#bar{p}#pi} (GeV/c^{2});Counts",
HistType::kTH1F, {{200, 1.08, 1.16}});
registryDeDx.add("hMassGamma", "Invariant mass #gamma;m_{ee} (GeV/c^{2});Counts",
HistType::kTH1F, {{200, 0.0, 0.1}});

// Armenteros-Podolanski plot
registryDeDx.add("hArmenterosAll", "Armenteros-Podolanski (all V0s);#alpha;q_{T} (GeV/c)",
HistType::kTH2F, {{200, -1, 1}, {200, 0, 0.3}});

// Armenteros-Podolanski plot by particle
registryDeDx.add("hArmenterosK0s", "Armenteros-Podolanski K0s;#alpha;q_{T} (GeV/c)",
HistType::kTH2F, {{200, -1, 1}, {200, 0, 0.3}});
registryDeDx.add("hArmenterosLambda", "Armenteros-Podolanski #Lambda;#alpha;q_{T} (GeV/c)",
HistType::kTH2F, {{200, -1, 1}, {200, 0, 0.3}});
registryDeDx.add("hArmenterosAntiLambda", "Armenteros-Podolanski #bar{#Lambda};#alpha;q_{T} (GeV/c)",
HistType::kTH2F, {{200, -1, 1}, {200, 0, 0.3}});
registryDeDx.add("hArmenterosGamma", "Armenteros-Podolanski #gamma;#alpha;q_{T} (GeV/c)",
HistType::kTH2F, {{200, -1, 1}, {200, 0, 0.3}});
}

registryDeDx.add(
Expand Down Expand Up @@ -663,6 +693,23 @@ struct DedxPidAnalysis {
registryDeDx.add("Tracks_vs_pT_all", "pT All", HistType::kTH1F, {{ptAxis}});
registryDeDx.add("Tracks_vs_pT_all_cuts", "pT All + cuts", HistType::kTH1F, {{ptAxis}});

// NCluster distributions
registryDeDx.add("hTPCClustersBefore", "N clusters TPC found Before", HistType::kTH1F, {{200, 0, 200, "N_{cl,found, Before}"}});
registryDeDx.add("hTPCClustersAfter", "N clusters TPC found After", HistType::kTH1F, {{200, 0, 200, "N_{cl,found, After}"}});

registryDeDx.add("hTPCPIDBefore", "N clusters TPC PID Before", HistType::kTH1F, {{200, 0, 200, "N_{cl,PID, Before}"}});
registryDeDx.add("hTPCPIDAfter", "N clusters TPC PID After", HistType::kTH1F, {{200, 0, 200, "N_{cl,PID, After}"}});

// DCA cut
registryDeDx.add("hDCAxyVsPt_before", "DCAxy vs pT before cut;#it{p}_{T} (GeV/c);DCA_{xy} (cm)",
HistType::kTH2F, {{ptAxis}, {200, -0.5, 0.5}});
registryDeDx.add("hDCAzVsPt_before", "DCAz vs pT before cut;#it{p}_{T} (GeV/c);DCA_{z} (cm)",
HistType::kTH2F, {{ptAxis}, {200, -0.5, 0.5}});
registryDeDx.add("hDCAxyVsPt_after", "DCAxy vs pT after cut;#it{p}_{T} (GeV/c);DCA_{xy} (cm)",
HistType::kTH2F, {{ptAxis}, {200, -0.5, 0.5}});
registryDeDx.add("hDCAzVsPt_after", "DCAz vs pT after cut;#it{p}_{T} (GeV/c);DCA_{z} (cm)",
HistType::kTH2F, {{ptAxis}, {200, -0.5, 0.5}});

// Event Counter
registryDeDx.add("evsel", "events selected", HistType::kTH1F, {{6, 0.5, 6.5, ""}});
auto hstat = registryDeDx.get<TH1>(HIST("evsel"));
Expand Down Expand Up @@ -721,6 +768,7 @@ struct DedxPidAnalysis {
xSec->SetBinLabel(V0TypeGamma, "V0TypeGamma");
xSec->SetBinLabel(V0RapidityGamma, "V0RapidityGamma");
xSec->SetBinLabel(MassCutGamma, "MassCutGamma");

mySelectionPrim = myTrackSelection();
}

Expand All @@ -740,7 +788,20 @@ struct DedxPidAnalysis {

return true;
}

// DCA xy cut
template <typename T1>
bool passesDCAxyCut(const T1& track) const
{
const float maxDcaXY = nSigmaDCAxy.value * (dcaXYp0.value + dcaXYp1.value / std::pow(track.pt(), dcaXYp2.value)) / 3.0;
return std::abs(track.dcaXY()) < maxDcaXY;
}
// DCA z cut
template <typename T1>
bool passesDCAzCut(const T1& track) const
{
const float maxiDcaZ = nSigmaDCAz.value * (maxDCAz.value) / 3.0;
return std::abs(track.dcaZ()) < maxiDcaZ;
}
// Momentum
template <typename T1>
float getMomentum(const T1& track)
Expand Down Expand Up @@ -780,7 +841,7 @@ struct DedxPidAnalysis {
return false;
if (!passedSingleTrackSelection(ntrack, collision))
return false;
/*double sigmap = 0.0;
double sigmap = 0.0;
double sigman = 0.0;

if (v0SelectionMode == V0TPC) {
Expand All @@ -794,19 +855,20 @@ struct DedxPidAnalysis {
sigman = std::hypot(ntrack.tpcNSigmaPi(), ntrack.tofNSigmaPi());
}

if (ptrack.tpcInnerParam() > tpcCut) {
if (ptrack.tpcInnerParam() > tpcMomentumCut) {
if (!ptrack.hasTOF())
return false;
if (std::abs(sigmap) > nsigmaTOFmax)
if (std::abs(sigmap) > nsigmaMax)
return false;
}

if (ntrack.tpcInnerParam() > tpcCut) {
if (ntrack.tpcInnerParam() > tpcMomentumCut) {
if (!ntrack.hasTOF())
return false;
if (std::abs(sigman) > nsigmaTOFmax)
if (std::abs(sigman) > nsigmaMax)
return false;
}*/
}

if (fillHist)
registryDeDx.fill(HIST("trackselSec"), TrkSecCutLabel::SingleTrackSelectionK0s);

Expand Down Expand Up @@ -860,7 +922,7 @@ struct DedxPidAnalysis {
if (!passedSingleTrackSelection(ntrack, collision))
return false;

/*double sigmap = 0.0;
double sigmap = 0.0;
double sigman = 0.0;

if (v0SelectionMode == V0TPC) {
Expand All @@ -874,19 +936,19 @@ struct DedxPidAnalysis {
sigman = std::hypot(ntrack.tpcNSigmaPi(), ntrack.tofNSigmaPi());
}

if (ptrack.tpcInnerParam() > tpcCut) {
if (ptrack.tpcInnerParam() > tpcMomentumCut) {
if (!ptrack.hasTOF())
return false;
if (std::abs(sigmap) > nsigmaTOFmax)
if (std::abs(sigmap) > nsigmaMax)
return false;
}

if (ntrack.tpcInnerParam() > tpcCut) {
if (ntrack.tpcInnerParam() > tpcMomentumCut) {
if (!ntrack.hasTOF())
return false;
if (std::abs(sigman) > nsigmaTOFmax)
if (std::abs(sigman) > nsigmaMax)
return false;
}*/
}

if (fillHist)
registryDeDx.fill(HIST("trackselSec"), TrkSecCutLabel::SingleTrackSelectionLambda);
Expand Down Expand Up @@ -935,7 +997,7 @@ struct DedxPidAnalysis {
return false;
if (!passedSingleTrackSelection(ntrack, collision))
return false;
/*double sigmap = 0.0;
double sigmap = 0.0;
double sigman = 0.0;

if (v0SelectionMode == V0TPC) {
Expand All @@ -948,19 +1010,19 @@ struct DedxPidAnalysis {
sigmap = std::hypot(ptrack.tpcNSigmaPi(), ptrack.tofNSigmaPi());
sigman = std::hypot(ntrack.tpcNSigmaPr(), ntrack.tofNSigmaPr());
}
if (ptrack.tpcInnerParam() > tpcCut) {
if (ptrack.tpcInnerParam() > tpcMomentumCut) {
if (!ptrack.hasTOF())
return false;
if (std::abs(sigmap) > nsigmaTOFmax)
if (std::abs(sigmap) > nsigmaMax)
return false;
}

if (ntrack.tpcInnerParam() > tpcCut) {
if (ntrack.tpcInnerParam() > tpcMomentumCut) {
if (!ntrack.hasTOF())
return false;
if (std::abs(sigman) > nsigmaTOFmax)
if (std::abs(sigman) > nsigmaMax)
return false;
}*/
}

if (fillHist)
registryDeDx.fill(HIST("trackselSec"), TrkSecCutLabel::SingleTrackSelectionAntiLambda);
Expand Down Expand Up @@ -1011,7 +1073,7 @@ struct DedxPidAnalysis {
if (!passedSingleTrackSelection(ntrack, collision))
return false;

/*double sigmap = 0.0;
double sigmap = 0.0;
double sigman = 0.0;

if (v0SelectionMode == V0TPC) {
Expand All @@ -1025,19 +1087,19 @@ struct DedxPidAnalysis {
sigman = std::hypot(ntrack.tpcNSigmaEl(), ntrack.tofNSigmaEl());
}

if (ptrack.tpcInnerParam() > tpcCut) {
if (ptrack.tpcInnerParam() > tpcMomentumCut) {
if (!ptrack.hasTOF())
return false;
if (std::abs(sigmap) > nsigmaTOFmax)
if (std::abs(sigmap) > nsigmaMax)
return false;
}

if (ntrack.tpcInnerParam() > tpcCut) {
if (ntrack.tpcInnerParam() > tpcMomentumCut) {
if (!ntrack.hasTOF())
return false;
if (std::abs(sigman) > nsigmaTOFmax)
if (std::abs(sigman) > nsigmaMax)
return false;
}*/
}
const float gammaMass = 2 * MassElectron; // GeV/c^2

if (fillHist)
Expand Down Expand Up @@ -1462,6 +1524,7 @@ struct DedxPidAnalysis {
}

if (nINELSelectionMode == NoSelINEL) {
registryDeDx.fill(HIST("evsel"), EvCutLabel::INELgt);
} else if (nINELSelectionMode == SelINELgt0) {
if (!collision.isInelGt0())
return;
Expand Down Expand Up @@ -1490,10 +1553,28 @@ struct DedxPidAnalysis {
for (const auto& trk : tracks) {
registryDeDx.fill(HIST("trackselAll"), TrkPriCutLabel::AllPri);

// track Selection
// Ncl distribution before cuts
registryDeDx.fill(HIST("hTPCClustersBefore"), trk.tpcNClsFound());
registryDeDx.fill(HIST("hTPCPIDBefore"), trk.tpcNClsPID());

// track Selection wo DCA
if (!mySelectionPrim.IsSelected(trk))
continue;

// Before DCA cuts
registryDeDx.fill(HIST("hDCAxyVsPt_before"), trk.pt(), trk.dcaXY());
registryDeDx.fill(HIST("hDCAzVsPt_before"), trk.pt(), trk.dcaZ());

// DCA cuts
if (!passesDCAxyCut(trk))
continue;
if (!passesDCAzCut(trk))
continue;

// After DCA cuts
registryDeDx.fill(HIST("hDCAxyVsPt_after"), trk.pt(), trk.dcaXY());
registryDeDx.fill(HIST("hDCAzVsPt_after"), trk.pt(), trk.dcaZ());

registryDeDx.fill(HIST("trackselAll"), TrkPriCutLabel::SelectionPrim);
// For pt
registryDeDx.fill(HIST("Tracks_vs_pT_all"), trk.pt());
Expand Down Expand Up @@ -1527,6 +1608,10 @@ struct DedxPidAnalysis {
// For pt + cuts
registryDeDx.fill(HIST("Tracks_vs_pT_all_cuts"), trk.pt());

// Ncl distribution After all cuts
registryDeDx.fill(HIST("hTPCClustersAfter"), trk.tpcNClsFound());
registryDeDx.fill(HIST("hTPCPIDAfter"), trk.tpcNClsPID());

float signedP = trk.sign() * getMomentum(trk);
float signedpT = trk.sign() * trk.pt();

Expand Down Expand Up @@ -1718,12 +1803,17 @@ struct DedxPidAnalysis {

fillHist = true;

// Armenteros for all V0
registryDeDx.fill(HIST("hArmenterosAll"), v0.alpha(), v0.qtarm());

// K0s Selection
if (passedK0Selection(v0, negTrack, posTrack, collision)) {
registryDeDx.fill(HIST("heta_vs_NclPID_TPC_After_PionsK0s"), posTrack.eta(), posTrack.tpcNClsPID());
registryDeDx.fill(HIST("heta_vs_NclPID_TPC_After_PionsK0s"), negTrack.eta(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_PionsK0s"), posTrack.p(), posTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_PionsK0s"), negTrack.p(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hMassK0s"), v0.mK0Short());
registryDeDx.fill(HIST("hArmenterosK0s"), v0.alpha(), v0.qtarm());
for (int i = 0; i < EtaIntervals; ++i) {
if (negTrack.eta() > EtaCut[i] && negTrack.eta() < EtaCut[i + 1]) {
registryDeDx.fill(HIST(DedxvsMomentumNeg[1]), std::abs(signedPneg), negTrack.tpcSignal() * 50 / calibrationFactorNeg->at(i), negTrack.eta());
Expand All @@ -1740,6 +1830,8 @@ struct DedxPidAnalysis {
registryDeDx.fill(HIST("heta_vs_NclPID_TPC_After_PionsLambda"), negTrack.eta(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_ProtonsLambda"), posTrack.p(), posTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_PionsLambda"), negTrack.p(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hMassLambda"), v0.mLambda());
registryDeDx.fill(HIST("hArmenterosLambda"), v0.alpha(), v0.qtarm());
for (int i = 0; i < EtaIntervals; ++i) {
if (negTrack.eta() > EtaCut[i] && negTrack.eta() < EtaCut[i + 1]) {
registryDeDx.fill(HIST(DedxvsMomentumNeg[1]), std::abs(signedPneg), negTrack.tpcSignal() * 50 / calibrationFactorNeg->at(i), negTrack.eta());
Expand All @@ -1756,6 +1848,8 @@ struct DedxPidAnalysis {
registryDeDx.fill(HIST("heta_vs_NclPID_TPC_After_ProtonsLambda"), negTrack.eta(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_PionsLambda"), posTrack.p(), posTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_ProtonsLambda"), negTrack.p(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hMassAntiLambda"), v0.mAntiLambda());
registryDeDx.fill(HIST("hArmenterosAntiLambda"), v0.alpha(), v0.qtarm());
for (int i = 0; i < EtaIntervals; ++i) {
if (negTrack.eta() > EtaCut[i] && negTrack.eta() < EtaCut[i + 1]) {
registryDeDx.fill(HIST(DedxvsMomentumNeg[2]), std::abs(signedPneg), negTrack.tpcSignal() * 50 / calibrationFactorNeg->at(i), negTrack.eta());
Expand All @@ -1772,6 +1866,8 @@ struct DedxPidAnalysis {
registryDeDx.fill(HIST("heta_vs_NclPID_TPC_After_ElectronsGamma"), negTrack.eta(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_ElectronsGamma"), posTrack.p(), posTrack.tpcNClsPID());
registryDeDx.fill(HIST("hp_vs_NclPID_TPC_After_ElectronsGamma"), negTrack.p(), negTrack.tpcNClsPID());
registryDeDx.fill(HIST("hMassGamma"), v0.mGamma());
registryDeDx.fill(HIST("hArmenterosGamma"), v0.alpha(), v0.qtarm());
for (int i = 0; i < EtaIntervals; ++i) {
if (negTrack.eta() > EtaCut[i] && negTrack.eta() < EtaCut[i + 1]) {
registryDeDx.fill(HIST(DedxvsMomentumNeg[3]), std::abs(signedPneg), negTrack.tpcSignal() * 50 / calibrationFactorNeg->at(i), negTrack.eta());
Expand Down
Loading