Skip to content

[PWGCF] Add pi-d systematic configurable in HadNucleiFemto.cxx #18414

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
blacwovie wants to merge 2 commits into
base: master
Choose a base branch
from
+25 −15
Open

[PWGCF] Add pi-d systematic configurable in HadNucleiFemto.cxx #18414

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
b4711b8
Update HadNucleiFemto.cxx
blacwovie Oct 22, 2026
90a9a48
Update HadronNucleiTables.h
blacwovie Oct 22, 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
4 changes: 3 additions & 1 deletion PWGCF/Femto/FemtoNuclei/DataModel/HadronNucleiTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -183,6 +183,7 @@ DECLARE_SOA_COLUMN(ZVertex, zVertex, float);
DECLARE_SOA_COLUMN(Multiplicity, multiplicity, uint16_t);
DECLARE_SOA_COLUMN(CentFT0C, centFT0C, float);
DECLARE_SOA_COLUMN(MultiplicityFT0C, multiplicityFT0C, float);
DECLARE_SOA_COLUMN(MultNTracksPV, multNTracksPV, int32_t);

} // namespace hadron_nuclei_tables

Expand Down Expand Up @@ -268,7 +269,8 @@ DECLARE_SOA_TABLE(HadronNucleiMult, "AOD", "HADNUCLEIMULT",
hadron_nuclei_tables::ZVertex,
hadron_nuclei_tables::Multiplicity,
hadron_nuclei_tables::CentFT0C,
hadron_nuclei_tables::MultiplicityFT0C)
hadron_nuclei_tables::MultiplicityFT0C,
hadron_nuclei_tables::MultNTracksPV)

// Reduced hadron-hypertriton pair output. Reconstructable quantities such as
// k*, invariant mass and pointing variables are intentionally derived offline.
Expand Down
36 changes: 22 additions & 14 deletions PWGCF/Femto/FemtoNuclei/TableProducer/HadNucleiFemto.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -89,12 +89,11 @@ using std::array;
using CollBracket = o2::math_utils::Bracket<int>;
using HyperCandidates = aod::DataHypCandsWColl;
using HyperCandidatesMC = aod::MCHypCandsWColl;
using CollisionsFull = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::FT0Mults>;
// The standard Femto producer stores multNTracksPV as its collision
// multiplicity. Subscribe to PVMults only in the mixed-event process so the
// same vertex--multiplicity pool definition can optionally be reproduced.
using CollisionsFullWithPVMult = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::FT0Mults, aod::PVMults>;
using CollisionsFullMC = soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::CentFT0Cs, aod::FT0Mults>;
// Save the PV-track multiplicity alongside numContrib in the event-info output
// for both same-event and mixed-event pairs.
using CollisionsFull = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Cs, aod::FT0Mults, aod::PVMults>;
using CollisionsFullWithPVMult = CollisionsFull;
using CollisionsFullMC = soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::CentFT0Cs, aod::FT0Mults, aod::PVMults>;
using HadHyperCollisionsFull = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0As, aod::CentFT0Cs, aod::CentFT0Ms, aod::FT0Mults>;
using HadHyperCollisionsFullMC = soa::Join<aod::Collisions, aod::McCollisionLabels, aod::EvSels, aod::CentFT0As, aod::CentFT0Cs, aod::CentFT0Ms, aod::FT0Mults>;
using TrackCandidates = soa::Join<aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::TracksDCA, aod::TrackSelection, aod::pidTPCFullDe, aod::pidTPCFullTr, aod::pidTOFFullDe, aod::pidTOFFullTr, aod::pidTOFFullHe, aod::pidTPCFullPr, aod::pidTOFFullPr, aod::pidTPCFullPi, aod::pidTOFFullPi, aod::pidTPCFullKa, aod::pidTOFFullKa, aod::TOFSignal, aod::TOFEvTime>;
Expand Down Expand Up @@ -317,6 +316,7 @@ struct BufferedCollision {
uint16_t numContrib{0u};
float centFT0C{0.f};
float multFT0C{0.f};
int32_t multNTracksPV{0};
std::vector<BufferedTrack> nuclei;
std::vector<BufferedTrack> hadrons;
};
Expand Down Expand Up @@ -403,6 +403,12 @@ struct HadNucleiFemto {
// cppcheck-suppress unusedStructMember
std::string prefix{"deuteronPid"};
// Deuteron purity and PID cuts
Configurable<int> settingTPCNClsFoundMinDe{"settingTPCNClsFoundMinDe", 110, "Minimum found TPC clusters for deuterons"};
Configurable<int> settingTPCCrossedRowsMinDe{"settingTPCCrossedRowsMinDe", 100, "Minimum crossed TPC rows for deuterons"};
Configurable<float> settingDeDCAxyOffset{"settingDeDCAxyOffset", 0.004f, "Deuteron DCAxy cut offset in cm"};
Configurable<float> settingDeDCAxyPtCoeff{"settingDeDCAxyPtCoeff", 0.013f, "Deuteron DCAxy inverse-pT cut coefficient"};
Configurable<float> settingDeDCAzOffset{"settingDeDCAzOffset", 0.004f, "Deuteron DCAz cut offset in cm"};
Configurable<float> settingDeDCAzPtCoeff{"settingDeDCAzPtCoeff", 0.013f, "Deuteron DCAz inverse-pT cut coefficient"};
Configurable<float> settingCutPinMinDe{"settingCutPinMinDe", 0.0f, "Minimum Pin for De"};
Configurable<float> settingCutDeptMin{"settingCutDeptMin", 0.6f, "Minimum PT cut on De"};
Configurable<float> settingCutDeptMax{"settingCutDeptMax", 1.6f, "Maximum PT cut on De"};
Expand Down Expand Up @@ -1230,16 +1236,14 @@ struct HadNucleiFemto {
return false;
}

constexpr int minTPCNClsFound = 110;
constexpr int minTPCNClsCrossedRows = 100;
constexpr float minTPCCrossedRowsOverFound = 0.f;
constexpr int maxTPCNClsShared = 160;
constexpr float maxSharedTPCFraction = 1.f;
constexpr int minITSNClsInnerBarrel = 1;
const float tpcCrossedRowsOverFound = candidate.tpcNClsFound() > 0 ? static_cast<float>(candidate.tpcNClsCrossedRows()) / candidate.tpcNClsFound() : 0.f;

return !(candidate.tpcNClsFound() < minTPCNClsFound ||
candidate.tpcNClsCrossedRows() < minTPCNClsCrossedRows ||
return !(candidate.tpcNClsFound() < deuteronPid.settingTPCNClsFoundMinDe ||
candidate.tpcNClsCrossedRows() < deuteronPid.settingTPCCrossedRowsMinDe ||
tpcCrossedRowsOverFound < minTPCCrossedRowsOverFound ||
candidate.tpcNClsShared() > maxTPCNClsShared ||
candidate.tpcFractionSharedCls() > maxSharedTPCFraction ||
Expand Down Expand Up @@ -1840,8 +1844,8 @@ struct HadNucleiFemto {
if (absPt <= 0.f) {
return false;
}
const float deDCAxyMax = 0.004f + 0.013f / absPt;
const float deDCAzMax = 0.004f + 0.013f / absPt;
const float deDCAxyMax = deuteronPid.settingDeDCAxyOffset + deuteronPid.settingDeDCAxyPtCoeff / absPt;
const float deDCAzMax = deuteronPid.settingDeDCAzOffset + deuteronPid.settingDeDCAzPtCoeff / absPt;
if (std::abs(candidate.dcaXY()) > deDCAxyMax || std::abs(candidate.dcaZ()) > deDCAzMax) {
return false;
}
Expand Down Expand Up @@ -2645,7 +2649,8 @@ struct HadNucleiFemto {
collision.posZ(),
collision.numContrib(),
collision.centFT0C(),
collision.multFT0C());
collision.multFT0C(),
collision.multNTracksPV());
}
}

Expand All @@ -2658,7 +2663,8 @@ struct HadNucleiFemto {
collision.posZ,
collision.numContrib,
collision.centFT0C,
collision.multFT0C);
collision.multFT0C,
collision.multNTracksPV);
}
}

Expand Down Expand Up @@ -3371,6 +3377,7 @@ struct HadNucleiFemto {
currentCollision.numContrib = collision.numContrib();
currentCollision.centFT0C = collision.centFT0C();
currentCollision.multFT0C = collision.multFT0C();
currentCollision.multNTracksPV = collision.multNTracksPV();

const uint64_t collIdx = collision.globalIndex();
auto tracksThisCollision = tracks.sliceBy(mPerCol, collIdx);
Expand Down Expand Up @@ -3548,6 +3555,7 @@ struct HadNucleiFemto {
current.numContrib = collision.numContrib();
current.centFT0C = collision.centFT0C();
current.multFT0C = collision.multFT0C();
current.multNTracksPV = collision.multNTracksPV();
row.eventId = current.eventId;
auto tracksThisCollision = tracks.sliceBy(mPerCol, collision.globalIndex());
tracksThisCollision.bindExternalIndices(&tracks);
Expand Down
Loading