Skip to content

[PWGLF] added event loss information. #12992

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
maciacco merged 8 commits into from
Sep 16, 2025
Merged

[PWGLF] added event loss information. #12992

Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
e3854c6
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 5, 2025
5ac67c1
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 5, 2025
a60d351
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 5, 2025
ee9f264
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 5, 2025
dc3b42c
Merge branch 'AliceO2Group:master' into master
jtanwar2212 Sep 10, 2025
4619be8
Merge branch 'AliceO2Group:master' into master
jtanwar2212 Sep 16, 2025
869f88f
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 16, 2025
4d32f32
Update nucleitpcpbpb.cxx
jtanwar2212 Sep 16, 2025
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
97 changes: 75 additions & 22 deletions PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.

Check failure on line 1 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/workflow-file] 

Name of a workflow file must match the name of the main struct in it (without the PWG prefix). (Class implementation files should be in "Core" directories.)
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
Expand Down Expand Up @@ -124,7 +124,7 @@
Configurable<float> cfgCutEta{"cfgCutEta", 0.9f, "Eta range for tracks"};
Configurable<bool> cfgetaRequire{"cfgetaRequire", true, "eta cut require"};
Configurable<bool> cfgetaRequireMC{"cfgetaRequireMC", true, "eta cut require for generated particles"};
Configurable<bool> cfgrapidityRequireMC{"cfgrapidityRequireMC", true, "rapidity cut require for generated particles"};
Configurable<bool> cfgRapidityRequireMC{"cfgRapidityRequireMC", true, "rapidity cut require for generated particles"};
Configurable<bool> cfgUsePVcontributors{"cfgUsePVcontributors", true, "use tracks that are PV contibutors"};
Configurable<bool> cfgITSrequire{"cfgITSrequire", true, "Additional cut on ITS require"};
Configurable<bool> cfgTPCrequire{"cfgTPCrequire", true, "Additional cut on TPC require"};
Expand All @@ -146,6 +146,7 @@
Configurable<bool> cfgDCAwithptRequire{"cfgDCAwithptRequire", true, "Require DCA cuts with pt dependance"};
Configurable<bool> cfgRequirebetaplot{"cfgRequirebetaplot", true, "Require beta plot"};
Configurable<bool> cfgIncludeMaterialInEfficiency{"cfgIncludeMaterialInEfficiency", true, "Require from material in efficiency"};
Configurable<bool> cfgMasscut{"cfgMasscut", true, "Require mass cut on He4 particles"};

Configurable<LabeledArray<double>> cfgBetheBlochParams{"cfgBetheBlochParams", {kBetheBlochDefault[0], nParticles, nBetheParams, particleNames, betheBlochParNames}, "TPC Bethe-Bloch parameterisation for light nuclei"};
Configurable<LabeledArray<double>> cfgTrackPIDsettings{"cfgTrackPIDsettings", {kTrackPIDSettings[0], nParticles, nTrkSettings, particleNames, trackPIDsettingsNames}, "track selection and PID criteria"};
Expand All @@ -155,7 +156,7 @@
Configurable<bool> cfgFillmassnsigma{"cfgFillmassnsigma", true, "Fill mass vs nsigma histograms"};
Configurable<float> centcut{"centcut", 80.0f, "centrality cut"};
Configurable<float> cfgCutRapidity{"cfgCutRapidity", 0.5f, "Rapidity range"};
Configurable<float> cfgtpcNClsFindable{"cfgtpcNClsFindable", 0.8f, "tpcNClsFindable over crossedRows"}; /////////////
Configurable<float> cfgtpcNClsFindable{"cfgtpcNClsFindable", 0.8f, "tpcNClsFindable over crossedRows"};
Configurable<float> cfgZvertex{"cfgZvertex", 10, "Min Z Vertex"};
Configurable<bool> cfgZvertexRequire{"cfgZvertexRequire", true, "Pos Z cut require"};
Configurable<bool> cfgZvertexRequireMC{"cfgZvertexRequireMC", true, "Pos Z cut require for generated particles"};
Expand All @@ -179,7 +180,7 @@
ConfigurableAxis speciesTrackingAxis{"speciesTrackingAxis", {11, -0.5, 10.5}, "particle type 0: pion, 1: proton, 2: deuteron, 3: triton, 4:He3, 5:He4"};
ConfigurableAxis axisDCA{"axisDCA", {400, -10., 10.}, "DCA axis"};
ConfigurableAxis particleAntiAxis{"particleAntiAxis", {2, 0, 2}, "Particle/Anti-particle"}; // 0 = particle, 1 = anti-particle
ConfigurableAxis decayTypeAxis{"decayTypeAxis", {3, 0, 3}, "Decay type"}; // 0 = primary, 1 = from decay
ConfigurableAxis decayTypeAxis{"decayTypeAxis", {3, -0.5, 2.5}, "Decay type"}; // 0 = primary, 1 = from decay, 2 = material

// CCDB
Service<o2::ccdb::BasicCCDBManager> ccdb;
Expand Down Expand Up @@ -281,6 +282,8 @@
histomc.add("histPIDtrackanti", " delta pt vs pt rec", HistType::kTH2F, {{1000, 0, 10, "p_{T}(reco)"}, {9, -0.5, 8.5, "p_{T}(reco) - p_{T}(gen)"}});
histomc.add("histPIDtrackhe4", " delta pt vs pt rec", HistType::kTH2F, {{1000, 0, 10, "p_{T}(reco)"}, {9, -0.5, 8.5, "p_{T}(reco) - p_{T}(gen)"}});
histomc.add("histPIDtrackantihe4", " delta pt vs pt rec", HistType::kTH2F, {{1000, 0, 10, "p_{T}(reco)"}, {9, -0.5, 8.5, "p_{T}(reco) - p_{T}(gen)"}});
histomc.add("hEventLossDenom", "Event loss denominator", kTH1F, {axisCent});
histomc.add("hEventLossNumer", "Event loss numerator", kTH1F, {axisCent});

histomc.add("histWeakDecayPtHe3", "Pt distribution of He3 from weak decays", kTH2F, {ptAxis, axisCent});
histomc.add("histWeakDecayPtAntiHe3", "Pt distribution of Anti-He3 from weak decays", kTH2F, {ptAxis, axisCent});
Expand All @@ -296,8 +299,6 @@
histomc.add("histProcess23Details", "Process 23 details", kTH2F, {{4, 0.5, 4.5, "particle type"}, {100, 0, 10, "p_{T}"}});
histomc.add("histAllMaterialSecondariesGen", "All material secondaries (gen)", kTH3F, {{100, 0, 10, "p_{T}"}, {20, -1, 1, "y"}, {5, -0.5, 4.5, "type"}});
histomc.add("histAllMaterialSecondariesReco", "All material secondaries (reco)", kTH3F, {{100, 0, 10, "p_{T}"}, {20, -1, 1, "y"}, {5, -0.5, 4.5, "type"}});

// Add axis labels for type: 0=unknown, 1=He3, 2=anti-He3, 3=He4, 4=anti-He4
}
}
//----------------------------------------------------------------------------------------------------------------
Expand Down Expand Up @@ -435,20 +436,50 @@

mcCollInfos.clear();
mcCollInfos.resize(mcCollisions.size());

// Store centrality regardless of cuts FIRST
for (auto const& collision : collisions) {
int mcCollIdx = collision.mcCollisionId();
if (mcCollIdx < 0 || mcCollIdx >= static_cast<int>(mcCollisions.size())) {
continue;
}
if (std::abs(collision.posZ()) > cfgZvertex)
// STORE CENTRALITY WITHOUT ANY CUTS
mcCollInfos[mcCollIdx].centrality = collision.centFT0C();
}

// FILL DENOMINATOR: ONCE per MC collision
for (size_t i = 0; i < mcCollInfos.size(); i++) {
if (mcCollInfos[i].centrality >= 0) { // Only if we found a matching collision
histomc.fill(HIST("hEventLossDenom"), mcCollInfos[i].centrality);
}
}

for (auto const& collision : collisions) {
int mcCollIdx = collision.mcCollisionId();
if (mcCollIdx < 0 || mcCollIdx >= static_cast<int>(mcCollisions.size())) {
continue;
if (!collision.sel8())
}
if (std::abs(collision.posZ()) > cfgZvertex && cfgZvertexRequireMC)
continue;
if (!collision.sel8() && cfgsel8Require)
continue;
if (collision.centFT0C() > centcut)
continue;

// Additional cuts
if (removeITSROFrameBorder && !collision.selection_bit(aod::evsel::kNoITSROFrameBorder))
continue;
if (removeNoSameBunchPileup && !collision.selection_bit(aod::evsel::kNoSameBunchPileup))
continue;
if (requireIsGoodZvtxFT0vsPV && !collision.selection_bit(aod::evsel::kIsGoodZvtxFT0vsPV))
continue;
if (requireIsVertexITSTPC && !collision.selection_bit(aod::evsel::kIsVertexITSTPC))
continue;
if (removeNoTimeFrameBorder && !collision.selection_bit(aod::evsel::kNoTimeFrameBorder))
continue;

mcCollInfos[mcCollIdx].passedEvSel = true;
mcCollInfos[mcCollIdx].centrality = collision.centFT0C();
histomc.fill(HIST("hEventLossNumer"), mcCollInfos[mcCollIdx].centrality);
}
for (auto const& mcCollision : mcCollisions) {
size_t idx = mcCollision.globalIndex();
Expand All @@ -464,10 +495,20 @@
int pdgCode = mcParticle.pdgCode();
bool isHe3 = (std::abs(pdgCode) == particlePdgCodes.at(4));
bool isHe4 = (std::abs(pdgCode) == particlePdgCodes.at(5));

if (std::abs(mcParticle.eta()) > cfgCutEta)
if (mcParticle.isPhysicalPrimary()) {
if (pdgCode == particlePdgCodes.at(4)) {
histomc.fill(HIST("histPtgenHe3"), mcParticle.pt());
} else if (pdgCode == -particlePdgCodes.at(4)) {
histomc.fill(HIST("histPtgenAntiHe3"), mcParticle.pt());
} else if (pdgCode == particlePdgCodes.at(5)) {
histomc.fill(HIST("histPtgenHe4"), mcParticle.pt());
} else if (pdgCode == -particlePdgCodes.at(5)) {
histomc.fill(HIST("histPtgenAntiHe4"), mcParticle.pt());
}
}
if (std::abs(mcParticle.eta()) > cfgCutEta && cfgetaRequireMC)
continue;
if (std::abs(mcParticle.y()) > cfgCutRapidity)
if (std::abs(mcParticle.y()) > cfgCutRapidity && cfgRapidityRequireMC)
continue;
bool isMaterialSecondary = false;
if (!mcParticle.isPhysicalPrimary() && (isHe3 || isHe4)) {
Expand Down Expand Up @@ -543,7 +584,7 @@
particleType = he4;

if (particleType >= 0) {
float centrality = mcCollInfos[idx].passedEvSel ? mcCollInfos[idx].centrality : -1.0f;
float centrality = mcCollInfos[idx].centrality; // Always use actual centrality
histomc.fill(HIST("hDenomSignalLoss"), particleType, mcParticle.pt(), mcParticle.y(), centrality, particleAnti, decayType);

if (mcCollInfos[idx].passedEvSel) {
Expand All @@ -563,15 +604,6 @@
histomc.fill(HIST("histWeakDecayPtAntiHe4"), mcParticle.pt(), centrality);
}
}
if (pdgCode == particlePdgCodes.at(4)) {
histomc.fill(HIST("histPtgenHe3"), mcParticle.pt());
} else if (pdgCode == -particlePdgCodes.at(4)) {
histomc.fill(HIST("histPtgenAntiHe3"), mcParticle.pt());
} else if (pdgCode == particlePdgCodes.at(5)) {
histomc.fill(HIST("histPtgenHe4"), mcParticle.pt());
} else if (pdgCode == -particlePdgCodes.at(5)) {
histomc.fill(HIST("histPtgenAntiHe4"), mcParticle.pt());
}
}
}
for (auto const& mcCollision : mcCollisions) {
Expand Down Expand Up @@ -933,6 +965,7 @@
{
if (!track.hasTOF() || !cfgFillmass)
return;

float beta{o2::pid::tof::Beta::GetBeta(track)};
const float eps = 1e-6f;
if (beta < eps || beta > 1.0f - eps)
Expand All @@ -942,7 +975,16 @@
float massTOF = p * charge * std::sqrt(1.f / (beta * beta) - 1.f);
// get PDG mass
float pdgMass = particleMasses[species];
float massDiff = massTOF - pdgMass;
float massDiff = 0.0;
if (species != he4) {
massDiff = massTOF - pdgMass;
}
if (species == he4) {
if (cfgMasscut && (massTOF * massTOF > 6.5 && massTOF * massTOF < 9.138))

Check failure on line 983 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.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.
return;
massDiff = massTOF - pdgMass;
}

float ptMomn;
setTrackParCov(track, mTrackParCov);
mTrackParCov.setPID(track.pidForTracking());
Expand All @@ -966,8 +1008,19 @@
float charge = (species == he3 || species == he4) ? 2.f : 1.f;
float p = getRigidity(track);
float massTOF = p * charge * std::sqrt(1.f / (beta * beta) - 1.f);

// get PDG mass
float masssquare = massTOF * massTOF;

if (species != he4) {
masssquare = massTOF * massTOF;
}
if (species == he4) {
if (cfgMasscut && (massTOF * massTOF > 6.5 && massTOF * massTOF < 9.138))

Check failure on line 1019 in PWGLF/Tasks/Nuspex/nucleitpcpbpb.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.
return;
masssquare = massTOF * massTOF;
}

if (track.sign() > 0) {
hmassnsigma[2 * species]->Fill(sigma, masssquare);
} else if (track.sign() < 0) {
Expand Down
Loading