Skip to content

[PWGLF] NucleiTask - Add mothers histo for MC #12863

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 5 commits into from
Sep 4, 2025
Merged

[PWGLF] NucleiTask - Add mothers histo for MC #12863

Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
cccfd9b
NucleiTask - Add mothers histo for MC
giovannimalfattore Sep 4, 2025
2efa62c
NucleiTask - Fix MEGAlinter error
giovannimalfattore Sep 4, 2025
f98003d
NucleiTask - Fix MEGAlinter error
giovannimalfattore Sep 4, 2025
bcf6b1a
NucleiTask - Fix format issues
giovannimalfattore Sep 4, 2025
9eaba2d
NucleiTask - Partial fit of O2Linter
giovannimalfattore Sep 4, 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
192 changes: 139 additions & 53 deletions PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,8 +42,9 @@
#include "ReconstructionDataFormats/Track.h"

#include <TF1.h>
// #include <TLorentzVector.h>

#include <gsl/span>
#include <memory>
#include <string>

using namespace o2;
Expand Down Expand Up @@ -221,12 +222,30 @@
static constexpr int PDGTriton = o2::constants::physics::Pdg::kTriton;
static constexpr int PDGHelium = o2::constants::physics::Pdg::kHelium3;
static constexpr int PDGAlpha = o2::constants::physics::Pdg::kAlpha;
static constexpr int PDGHyperTriton = o2::constants::physics::Pdg::kHyperTriton;
static constexpr float MassProtonVal = o2::constants::physics::MassProton;
static constexpr float MassDeuteronVal = o2::constants::physics::MassDeuteron;
static constexpr float MassTritonVal = o2::constants::physics::MassTriton;
static constexpr float MassHeliumVal = o2::constants::physics::MassHelium3;
static constexpr float MassAlphaVal = o2::constants::physics::MassAlpha;

// PDG of Mothers
static constexpr int kPdgMotherlist[] = {
PDGProton, // proton
PDGPion, // pi+
PDGKaon, // K+
311, // K0
PDGDeuteron, // deuteron
PDGTriton, // triton
PDGHelium, // He-3
PDGAlpha, // Alpha
1000130270, // Aluminium
1000140280, // Silicon
1000260560 // Iron
};
static constexpr int kNumMotherlist = sizeof(kPdgMotherlist) / sizeof(kPdgMotherlist[0]);
static constexpr int kMaxNumMom = 4; // X: 0..4, overflow=5

template <typename TrackType>
float averageClusterSizeTrk(const TrackType& track)
{
Expand Down Expand Up @@ -825,24 +844,12 @@
histos.add<TH1>("tracks/triton/h1antiTritonSpectraTrueTransport", "#it{p}_{T} (#bar{t})", HistType::kTH1F, {ptAxis});
}
if (enableHe) {
// histos.add<TH1>("tracks/helium/h1HeliumSpectraTrue", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueWPID", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1HeliumSpectraTruePrim", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueSec", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueTransport", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});

histos.add<TH1>("tracks/helium/h1HeliumSpectraTrue_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueWPID_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1HeliumSpectraTruePrim_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueSec_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1HeliumSpectraTrueTransport_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});

// histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrue", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrueWPID", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTruePrim", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrueSec", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});
// histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrueTransport", "#it{p}_{T}/z (He)", HistType::kTH1F, {ptZHeAxis});

histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrue_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTrueWPID_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
histos.add<TH1>("tracks/helium/h1antiHeliumSpectraTruePrim_Z2", "#it{p}_{T} (He)", HistType::kTH1F, {ptHeAxis});
Expand Down Expand Up @@ -1138,6 +1145,25 @@
histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtHeliumTruePrim", "DCAxy vs Pt (He); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});
histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueSec", "DCAxy vs Pt (He); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});
histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueMaterial", "DCAxy vs Pt (He); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});

histos.add<TH2>("tracks/helium/dca/before/hMomTrueMaterial", "MC mothers;mother index;mother PDG", HistType::kTH2I, {{kMaxNumMom + 2, -0.5, static_cast<double>(kMaxNumMom) + 1.5}, {kNumMotherlist + 2, -1.5, static_cast<double>(kNumMotherlist) + 0.5}});

// Fix for getting TH2 pointer
std::shared_ptr<TH2> hTemp = histos.get<TH2>(HIST("tracks/helium/dca/before/hMomTrueMaterial"));
TH2* hPDG = hTemp.get();

TAxis* axPDG = hPDG->GetXaxis();
for (int i = 0; i <= kMaxNumMom; ++i) {
axPDG->SetBinLabel(i + 1, Form("%d", i));
}
axPDG->SetBinLabel(kMaxNumMom + 2, ">=5");
TAxis* ayPDG = hPDG->GetYaxis();
ayPDG->SetBinLabel(1, "-1"); // undefined
ayPDG->SetBinLabel(2, "0"); // other
for (int i = 0; i < kNumMotherlist; ++i) {
ayPDG->SetBinLabel(i + 3, Form("%d", kPdgMotherlist[i]));
}

histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueTransport", "DCAxy vs Pt (He); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});

histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtantiHeliumTrue", "DCAxy vs Pt (#bar{He}); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});
Expand Down Expand Up @@ -2099,7 +2125,7 @@
template <bool IsMC, bool IsFilteredData, typename CollisionType, typename TracksType, typename ParticleType>
void fillHistograms(const CollisionType& event,
const TracksType& tracks,
const ParticleType& /*particles*/)
const ParticleType& particles)
{
histos.fill(HIST("event/eventSkimming"), 0.5);
// Apply skimming
Expand Down Expand Up @@ -2206,7 +2232,7 @@

tracks.copyIndexBindings(tracksWithITS);

for (auto& track : tracksWithITS) {
for (auto const& track : tracksWithITS) {
if constexpr (!IsFilteredData) {
if (!track.isGlobalTrackWoDCA() && filterOptions.enableIsGlobalTrack) {
continue;
Expand All @@ -2219,14 +2245,13 @@
histos.fill(HIST("tracks/avgClusterSizePerCoslInvVsITSlayers"), track.p(), averageClusterSizePerCoslInv(track), track.itsNCls());
}

if (track.itsNCls() < trkqcOptions.cfgCutITSClusters)
continue;
if (track.tpcNClsCrossedRows() < trkqcOptions.cfgCutTPCXRows)
continue;
if (track.tpcNClsFound() < trkqcOptions.cfgCutTPCClusters)
continue;
if (track.tpcCrossedRowsOverFindableCls() < trkqcOptions.cfgCutTPCCROFnd)
if (track.itsNCls() < trkqcOptions.cfgCutITSClusters ||
track.tpcNClsCrossedRows() < trkqcOptions.cfgCutTPCXRows ||
track.tpcNClsFound() < trkqcOptions.cfgCutTPCClusters ||
track.tpcCrossedRowsOverFindableCls() < trkqcOptions.cfgCutTPCCROFnd) {
continue;
}

auto tpcChi2NclRange = (std::vector<float>)trkqcOptions.tpcChi2NclCuts;
if ((track.tpcChi2NCl() < tpcChi2NclRange[0]) || (track.tpcChi2NCl() > tpcChi2NclRange[1]))
continue;
Expand All @@ -2252,21 +2277,21 @@
float shiftPtNeg = 0.f;

if (enablePtShiftHe && !fShiftPtHe) {
fShiftPtHe = new TF1("fShiftPtHe", "[0] * TMath::Exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto par = (std::vector<float>)parShiftPtHe;
fShiftPtHe->SetParameters(par[0], par[1], par[2], par[3], par[4]);
fShiftPtHe = new TF1("fShiftPtHe", "[0] * exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto parHe = (std::vector<float>)parShiftPtHe; // NOLINT
fShiftPtHe->SetParameters(parHe[0], parHe[1], parHe[2], parHe[3], parHe[4]);
}

if (enablePtShiftHe && !fShiftPtantiHe) {
fShiftPtantiHe = new TF1("fShiftPtantiHe", "[0] * TMath::Exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto par = (std::vector<float>)parShiftPtAntiHe;
fShiftPtantiHe->SetParameters(par[0], par[1], par[2], par[3], par[4]);
fShiftPtantiHe = new TF1("fShiftPtantiHe", "[0] * exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto parAntiHe = (std::vector<float>)parShiftPtAntiHe; // NOLINT
fShiftPtantiHe->SetParameters(parAntiHe[0], parAntiHe[1], parAntiHe[2], parAntiHe[3], parAntiHe[4]);
}

if (enablePtShiftAntiD && !fShiftAntiD) {
fShiftAntiD = new TF1("fShiftAntiD", "[0] * TMath::Exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto par = (std::vector<float>)parShiftPtAntiD;
fShiftAntiD->SetParameters(par[0], par[1], par[2], par[3], par[4]);
fShiftAntiD = new TF1("fShiftAntiD", "[0] * exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto parAntiD = (std::vector<float>)parShiftPtAntiD; // NOLINT
fShiftAntiD->SetParameters(parAntiD[0], parAntiD[1], parAntiD[2], parAntiD[3], parAntiD[4]);
}

switch (unableAntiDPtShift) {
Expand All @@ -2282,9 +2307,9 @@
}

if (enablePtShiftD && !fShiftD) {
fShiftD = new TF1("fShiftD", "[0] * TMath::Exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto par = (std::vector<float>)parShiftPtD;
fShiftD->SetParameters(par[0], par[1], par[2], par[3], par[4]);
fShiftD = new TF1("fShiftD", "[0] * exp([1] + [2] * x) + [3] + [4] * x", 0.f, 8.f);
auto parD = (std::vector<float>)parShiftPtD; // NOLINT
fShiftD->SetParameters(parD[0], parD[1], parD[2], parD[3], parD[4]);
}

switch (unableDPtShift) {
Expand Down Expand Up @@ -2475,11 +2500,16 @@
}

// Rapidity cuts
prRapCut = track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton)) > kinemOptions.cfgRapidityCutLow && track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton)) < kinemOptions.cfgRapidityCutHigh;
deRapCut = track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron)) > kinemOptions.cfgRapidityCutLow && track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron)) < kinemOptions.cfgRapidityCutHigh;
trRapCut = track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Triton)) > kinemOptions.cfgRapidityCutLow && track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Triton)) < kinemOptions.cfgRapidityCutHigh;
heRapCut = track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)) > kinemOptions.cfgRapidityCutLow && track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)) < kinemOptions.cfgRapidityCutHigh;
alRapCut = track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Alpha)) > kinemOptions.cfgRapidityCutLow && track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Alpha)) < kinemOptions.cfgRapidityCutHigh;
auto rapCheck = [&](float m2z) {
const float rap = track.rapidity(m2z);
return (rap > kinemOptions.cfgRapidityCutLow) && (rap < kinemOptions.cfgRapidityCutHigh);
};

prRapCut = rapCheck(MassProtonVal);
deRapCut = rapCheck(MassDeuteronVal);
trRapCut = rapCheck(MassTritonVal);
heRapCut = rapCheck(MassHeliumVal / 2.0);
alRapCut = rapCheck(MassAlphaVal / 2.0);

isDeuteron = enableDe && deRapCut;
isHelium = enableHe && heRapCut;
Expand Down Expand Up @@ -2615,6 +2645,7 @@
}

if constexpr (IsMC) {
// auto const& mcParticles = particles;
bool isPhysPrim = false;
bool isProdByGen = false;
bool isWeakDecay = false;
Expand All @@ -2624,15 +2655,15 @@
if constexpr (IsFilteredData) {
isPhysPrim = track.isPhysicalPrimary();
isProdByGen = track.producedByGenerator();
isWeakDecay = track.getProcess() == 4;
isWeakDecay = track.getProcess() == 4; // NOLINT

Check failure on line 2658 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
pdgCode = track.pdgCode();
} else {
if (!track.has_mcParticle()) {
continue;
}
isPhysPrim = track.mcParticle().isPhysicalPrimary();
isProdByGen = track.mcParticle().producedByGenerator();
isWeakDecay = track.mcParticle().getProcess() == 4;
isWeakDecay = track.mcParticle().getProcess() == 4; // NOLINT

Check failure on line 2666 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
pdgCode = track.mcParticle().pdgCode();
}

Expand Down Expand Up @@ -3037,6 +3068,8 @@
break;
}
}
} else {
(void)particles;
}
// Tracks DCA histos fill
if (outFlagOptions.makeDCABeforeCutPlots) {
Expand Down Expand Up @@ -3114,19 +3147,27 @@
}

if constexpr (IsMC) {
// auto const& mcParticles = particles;
bool isPhysPrim = false;
bool isProdByGen = false;
bool isWeakDecay = false;
bool hasFakeHit = false;

// PID
int pdgCode = 0;
int pdgMom = 0;
// gen Pt
float genPt = 0;
// Mothers variables
[[maybe_unused]] int firstMotherId = -1;
[[maybe_unused]] int firstMotherPdg = -1;
[[maybe_unused]] int pdgList[8];
[[maybe_unused]] int nSaved = 0;

if constexpr (IsFilteredData) {
isPhysPrim = track.isPhysicalPrimary();
isProdByGen = track.producedByGenerator();
isWeakDecay = track.getProcess() == 4;
isWeakDecay = track.getProcess() == 4; // NOLINT

Check failure on line 3170 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
pdgCode = track.pdgCode();
genPt = std::sqrt(std::pow(track.px(), 2) + std::pow(track.py(), 2));

Expand All @@ -3136,11 +3177,37 @@
}
isPhysPrim = track.mcParticle().isPhysicalPrimary();
isProdByGen = track.mcParticle().producedByGenerator();
isWeakDecay = track.mcParticle().getProcess() == 4;
isWeakDecay = track.mcParticle().getProcess() == 4; // NOLINT

Check failure on line 3180 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
pdgCode = track.mcParticle().pdgCode();
genPt = track.mcParticle().pt();

// Access to MC particles mother
o2::aod::McParticles::iterator mc = particles.iteratorAt(track.mcParticleId());
gsl::span<const int> motherIds = mc.mothersIds();
const int nMothers = static_cast<int>(motherIds.size());
firstMotherId = -1;
firstMotherPdg = -1;

nSaved = 0;

for (int iMom = 0; iMom < nMothers; ++iMom) {
int motherId = motherIds[iMom];
if (motherId < 0 || motherId >= particles.size()) {
continue; // added check on mother
}
o2::aod::McParticles::iterator mother = particles.iteratorAt(motherId);
pdgMom = mother.pdgCode();

if (iMom == 0) {
firstMotherId = motherId;
firstMotherPdg = pdgMom;
}
if (nSaved < 8) {

Check failure on line 3204 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
pdgList[nSaved++] = pdgMom;
}
}

genPt = track.mcParticle().pt();
for (int i = 0; i < 10; i++) { // From ITS to TPC

Check failure on line 3210 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
if (track.mcMask() & 1 << i) {
hasFakeHit = true;
break;
Expand Down Expand Up @@ -3400,10 +3467,30 @@
}
if (!isPhysPrim && !isProdByGen && outFlagOptions.makeDCABeforeCutPlots) {
histos.fill(HIST("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueTransport"), hePt, track.dcaXY());
if (isWeakDecay)
if (isWeakDecay) {
histos.fill(HIST("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueSec"), hePt, track.dcaXY());
else
} else {
histos.fill(HIST("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueMaterial"), hePt, track.dcaXY());
if (!IsFilteredData) {
if (nSaved > 0) {
for (int i = 0; i < nSaved; ++i) {
int idxComp = (i <= kMaxNumMom) ? i : (kMaxNumMom + 1);
int pdgMom = pdgList[i];
int yVal = -1;
if (pdgMom != -1) {
yVal = 0;
for (int j = 0; j < kNumMotherlist; ++j) {
if (kPdgMotherlist[j] == pdgMom) {
yVal = j + 1;
break;
}
}
}
histos.fill(HIST("tracks/helium/dca/before/hMomTrueMaterial"), idxComp, yVal);
}
}
}
}
if (track.hasTOF() && outFlagOptions.doTOFplots) {
histos.fill(HIST("tracks/helium/dca/before/TOF/hDCAxyVsPtHeliumTrueTransport"), hePt, track.dcaXY());
if (isWeakDecay)
Expand Down Expand Up @@ -3682,6 +3769,8 @@
}
break;
}
} else {
(void)particles;
}

// DCA Cut
Expand Down Expand Up @@ -3813,10 +3902,10 @@
debugHistos.fill(HIST("debug/qa/h2TPCncrVsPtPos"), track.tpcInnerParam(), track.tpcNClsCrossedRows());
debugHistos.fill(HIST("debug/qa/h2TPCncrVsTPCsignalPos"), track.tpcSignal(), track.tpcNClsCrossedRows());

if (track.tpcInnerParam() < 0.5f) {

Check failure on line 3905 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
debugHistos.fill(HIST("debug/qa/h1TPCncrLowPPos"), track.tpcNClsCrossedRows());
}
if ((track.tpcInnerParam() >= 0.5f) && (track.tpcInnerParam() < 1.f)) {

Check failure on line 3908 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
debugHistos.fill(HIST("debug/qa/h1TPCncrMidPPos"), track.tpcNClsCrossedRows());
}
if (track.tpcInnerParam() >= 1.f) {
Expand All @@ -3826,7 +3915,7 @@
debugHistos.fill(HIST("debug/qa/h2TPCncrVsPtNeg"), track.tpcInnerParam(), track.tpcNClsCrossedRows());
debugHistos.fill(HIST("debug/qa/h2TPCncrVsTPCsignalNeg"), track.tpcSignal(), track.tpcNClsCrossedRows());

if (track.tpcInnerParam() < 0.5f) {

Check failure on line 3918 in PWGLF/Tasks/Nuspex/LFNucleiBATask.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.
debugHistos.fill(HIST("debug/qa/h1TPCncrLowPNeg"), track.tpcNClsCrossedRows());
}
if ((track.tpcInnerParam() >= 0.5f) && (track.tpcInnerParam() < 1.f)) {
Expand Down Expand Up @@ -4553,11 +4642,8 @@
histos.fill(HIST("tracks/eff/helium/hPtHe"), 2 * hePt);
histos.fill(HIST("tracks/eff/helium/h2pVsTPCmomentumHe"), heTPCmomentum, heP);
}
// histos.fill(HIST("tracks/helium/h1HeliumSpectra"), hePt);
histos.fill(HIST("tracks/helium/h1HeliumSpectra_Z2"), 2 * hePt);
// histos.fill(HIST("tracks/helium/h2HeliumYvsPt"), track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)), hePt);
histos.fill(HIST("tracks/helium/h2HeliumYvsPt_Z2"), track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)), 2 * hePt);
// histos.fill(HIST("tracks/helium/h2HeliumEtavsPt"), track.eta(), hePt);
histos.fill(HIST("tracks/helium/h2HeliumEtavsPt_Z2"), track.eta(), 2 * hePt);
if (outFlagOptions.enablePIDplot)
histos.fill(HIST("tracks/helium/h2TPCsignVsTPCmomentumHelium"), heTPCmomentum, track.tpcSignal());
Expand All @@ -4567,11 +4653,8 @@
histos.fill(HIST("tracks/eff/helium/hPtantiHe"), 2 * antihePt);
histos.fill(HIST("tracks/eff/helium/h2pVsTPCmomentumantiHe"), antiheTPCmomentum, antiheP);
}
// histos.fill(HIST("tracks/helium/h1antiHeliumSpectra"), antihePt);
histos.fill(HIST("tracks/helium/h1antiHeliumSpectra_Z2"), 2 * antihePt);
// histos.fill(HIST("tracks/helium/h2antiHeliumYvsPt"), track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)), antihePt);
histos.fill(HIST("tracks/helium/h2antiHeliumYvsPt_Z2"), track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3)), 2 * antihePt);
// histos.fill(HIST("tracks/helium/h2antiHeliumEtavsPt"), track.eta(), antihePt);
histos.fill(HIST("tracks/helium/h2antiHeliumEtavsPt_Z2"), track.eta(), 2 * antihePt);
if (outFlagOptions.enablePIDplot)
histos.fill(HIST("tracks/helium/h2TPCsignVsTPCmomentumantiHelium"), antiheTPCmomentum, track.tpcSignal());
Expand Down Expand Up @@ -4811,6 +4894,7 @@
}

if constexpr (IsMC) {
// auto const& mcParticles = particles;
bool isPhysPrim = false;
bool isProdByGen = false;
bool isWeakDecay = false;
Expand Down Expand Up @@ -5733,6 +5817,8 @@
default:
break;
}
} else {
(void)particles;
}
}
}
Expand Down Expand Up @@ -5973,7 +6059,7 @@
spectraGen.fill(HIST("histGenVetxZ"), mcCollision.posZ());
if (mcCollision.centFT0M() < cfgMultCutLow || mcCollision.centFT0M() > cfgMultCutHigh)
return;
for (auto& mcParticleGen : mcParticles) {
for (auto& mcParticleGen : mcParticles) { // NOLINT

Check failure on line 6062 in PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[const-ref-in-for-loop] 

Use constant references for non-modified iterators in range-based for loops.
if (mcParticleGen.y() > kinemOptions.cfgRapidityCutHigh || mcParticleGen.y() < kinemOptions.cfgRapidityCutLow) {
continue;
}
Expand Down
Loading