Skip to content

[PWGLF] Possible change in MC trigger definition #16085

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
lhusova merged 2 commits into from
May 4, 2026
+64 −40
Merged

[PWGLF] Possible change in MC trigger definition #16085

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
a82720d
Possible change in trigger definition
scannito May 4, 2026
6250c24
Fix
scannito 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
104 changes: 64 additions & 40 deletions PWGLF/TableProducer/Strangeness/phiStrangeCorrelator.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/TableProducer/Strangeness/phiStrangeCorrelator.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 @@ -96,6 +96,8 @@
Configurable<float> cfgYAcceptance{"cfgYAcceptance", 0.5f, "Rapidity acceptance"};
} phiConfigs;

Configurable<bool> trueGenPhi{"trueGenPhi", false, "Phi trigger in MC Collisions: false - gen K+K- pair, true - gen phi meson"};

// Configurables on phi pT bins
Configurable<std::vector<double>> binspTPhi{"binspTPhi", {0.4, 0.8, 1.4, 2.0, 2.8, 4.0, 6.0, 10.0}, "pT bin limits for Phi"};

Expand Down Expand Up @@ -315,28 +317,41 @@

void processMCGen(aod::McCollisions::iterator const& mcCollision, aod::McParticles const& mcParticles)
{
for (const auto& mcParticle1 : mcParticles) {
if (!mcParticle1.isPhysicalPrimary() || std::abs(mcParticle1.eta()) > trackConfigs.etaMax)
continue;

for (const auto& mcParticle2 : mcParticles) {
if (!mcParticle2.isPhysicalPrimary() || std::abs(mcParticle2.eta()) > trackConfigs.etaMax)
if (trueGenPhi) {
for (const auto& mcParticle : mcParticles) {
if (std::abs(mcParticle.pdgCode()) != o2::constants::physics::Pdg::kPhi)
continue;
if (mcParticle.pt() < phiConfigs.minPhiPt)
continue;
if (std::abs(mcParticle.y()) > phiConfigs.cfgYAcceptance)
continue;

if (!(mcParticle1.pdgCode() == PDG_t::kKPlus && mcParticle2.pdgCode() == PDG_t::kKMinus) &&
!(mcParticle1.pdgCode() == PDG_t::kKMinus && mcParticle2.pdgCode() == PDG_t::kKPlus))
phimesonCandidatesMcGen(mcCollision.globalIndex(), 0, mcParticle.pt(), mcParticle.y(), mcParticle.phi());
}
} else {
for (const auto& mcParticle1 : mcParticles) {
if (!mcParticle1.isPhysicalPrimary() || std::abs(mcParticle1.eta()) > trackConfigs.etaMax)
continue;

ROOT::Math::PxPyPzMVector genKPair = recMother(mcParticle1, mcParticle2, massKa, massKa);
for (const auto& mcParticle2 : mcParticles) {
if (!mcParticle2.isPhysicalPrimary() || std::abs(mcParticle2.eta()) > trackConfigs.etaMax)
continue;

if (genKPair.Pt() < phiConfigs.minPhiPt)
continue;
if (genKPair.M() > phiConfigs.maxMPhi)
continue;
if (std::abs(genKPair.Rapidity()) > phiConfigs.cfgYAcceptance)
continue;
if (!(mcParticle1.pdgCode() == PDG_t::kKPlus && mcParticle2.pdgCode() == PDG_t::kKMinus) &&
!(mcParticle1.pdgCode() == PDG_t::kKMinus && mcParticle2.pdgCode() == PDG_t::kKPlus))
Comment on lines +340 to +341
continue;

ROOT::Math::PxPyPzMVector genKPair = recMother(mcParticle1, mcParticle2, massKa, massKa);

if (genKPair.Pt() < phiConfigs.minPhiPt)
continue;
if (genKPair.M() > phiConfigs.maxMPhi)
continue;
if (std::abs(genKPair.Rapidity()) > phiConfigs.cfgYAcceptance)
continue;

phimesonCandidatesMcGen(mcCollision.globalIndex(), genKPair.M(), genKPair.Pt(), genKPair.Rapidity(), genKPair.Phi());
phimesonCandidatesMcGen(mcCollision.globalIndex(), genKPair.M(), genKPair.Pt(), genKPair.Rapidity(), genKPair.Phi());
}
}
}
}
Expand Down Expand Up @@ -546,6 +561,7 @@
// Configurable<bool> forceTOF{"forceTOF", false, "force the TOF signal for the PID"};
Configurable<float> tofPIDThreshold{"tofPIDThreshold", 0.5, "minimum pT after which TOF PID is applicable"};
Configurable<std::vector<int>> trkPIDspecies{"trkPIDspecies", std::vector<int>{o2::track::PID::Kaon, o2::track::PID::Proton}, "Trk sel: Particles species for PID rejection, kaon, proton"};
Configurable<bool> applyElRejection{"applyElRejection", false, "Apply or not the electron rejection"};
Configurable<std::pair<float, float>> pidRangeTPCEl{"pidRangeTPCEl", {-3.0f, 5.0f}, "nSigma TPC range for electrons"};
Configurable<float> pidTPCMaxHadrons{"pidTPCMaxHadrons", 3.0f, "maximum nSigma TPC for hadrons"};
Configurable<float> pidTOFMaxHadrons{"pidTOFMaxHadrons", 3.0f, "maximum nSigma TOF for hadrons"};
Expand Down Expand Up @@ -632,23 +648,25 @@
template <typename T>
bool pidHypothesesRejection(const T& track)
{
// Electron rejection
auto nSigmaTPCEl = aod::pidutils::tpcNSigma(o2::track::PID::Electron, track);

if (nSigmaTPCEl > trackConfigs.pidRangeTPCEl->first && nSigmaTPCEl < trackConfigs.pidRangeTPCEl->second) {
auto nSigmaTPCPi = aod::pidutils::tpcNSigma(o2::track::PID::Pion, track);
auto nSigmaTPCKa = aod::pidutils::tpcNSigma(o2::track::PID::Kaon, track);
auto nSigmaTPCPr = aod::pidutils::tpcNSigma(o2::track::PID::Proton, track);

if (std::abs(nSigmaTPCPi) > trackConfigs.pidTPCMaxHadrons &&
std::abs(nSigmaTPCKa) > trackConfigs.pidTPCMaxHadrons &&
std::abs(nSigmaTPCPr) > trackConfigs.pidTPCMaxHadrons) {
return false;
// Electron rejection (if enabled)
if (trackConfigs.applyElRejection) {
auto nSigmaTPCEl = aod::pidutils::tpcNSigma(o2::track::PID::Electron, track);

if (nSigmaTPCEl > trackConfigs.pidRangeTPCEl->first && nSigmaTPCEl < trackConfigs.pidRangeTPCEl->second) {
auto nSigmaTPCPi = aod::pidutils::tpcNSigma(o2::track::PID::Pion, track);
auto nSigmaTPCKa = aod::pidutils::tpcNSigma(o2::track::PID::Kaon, track);
auto nSigmaTPCPr = aod::pidutils::tpcNSigma(o2::track::PID::Proton, track);

if (std::abs(nSigmaTPCPi) > trackConfigs.pidTPCMaxHadrons &&
std::abs(nSigmaTPCKa) > trackConfigs.pidTPCMaxHadrons &&
std::abs(nSigmaTPCPr) > trackConfigs.pidTPCMaxHadrons) {
return false;
}
}
}

// Other hadron species rejection
for (auto pid : {o2::track::PID::Kaon, o2::track::PID::Proton}) {

Check failure on line 669 in PWGLF/TableProducer/Strangeness/phiStrangeCorrelator.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.
auto nSigmaTPC = aod::pidutils::tpcNSigma(pid, track);

if (std::abs(nSigmaTPC) < trackConfigs.pidTPCMaxHadrons) {
Expand Down Expand Up @@ -806,6 +824,8 @@
Configurable<float> maxMPhiSignal{"maxMPhiSignal", 1.029f, "Upper limits on Phi mass for signal extraction"};
} phiConfigs;

Configurable<bool> withTrueGenPhi{"withTrueGenPhi", false, "Require the presence of a true generated phi meson in the event"};

// Defining the type of the collisions for data and MC
using SelCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::PVMults>;
using SimCollisions = soa::Join<SelCollisions, aod::McCollisionLabels>;
Expand Down Expand Up @@ -911,23 +931,27 @@
template <bool isMC, typename MC_T = void, typename T1, typename T2>
bool eventHasPhi(const T1& collision, const T2& phiCandidates)
{
uint16_t nPhi{0};
if (withTrueGenPhi) {
if (phiCandidates.size() < 1)
return false;
Comment on lines +934 to +936
} else {
uint16_t nPhi{0};

for (const auto& phiCand : phiCandidates) {
for (const auto& phiCand : phiCandidates) {
if (phiCand.inMassRegion(phiConfigs.minMPhiSignal, phiConfigs.maxMPhiSignal))
nPhi++;

if (phiCand.inMassRegion(phiConfigs.minMPhiSignal, phiConfigs.maxMPhiSignal))
nPhi++;
// histos.fill(HIST("hEta"), track1.eta());
// histos.fill(HIST("hNsigmaKaonTPC"), track1.tpcInnerParam(), track1.tpcNSigmaKa());
// histos.fill(HIST("hNsigmaKaonTOF"), track1.tpcInnerParam(), track1.tofNSigmaKa());
// histos.fill(HIST("h2DauTracksPhiDCAxy"), track1.pt(), track1.dcaXY());
// histos.fill(HIST("h2DauTracksPhiDCAz"), track1.pt(), track1.dcaZ());
}

// histos.fill(HIST("hEta"), track1.eta());
// histos.fill(HIST("hNsigmaKaonTPC"), track1.tpcInnerParam(), track1.tpcNSigmaKa());
// histos.fill(HIST("hNsigmaKaonTOF"), track1.tpcInnerParam(), track1.tofNSigmaKa());
// histos.fill(HIST("h2DauTracksPhiDCAxy"), track1.pt(), track1.dcaXY());
// histos.fill(HIST("h2DauTracksPhiDCAz"), track1.pt(), track1.dcaZ());
if (nPhi == 0)
return false;
}

if (nPhi == 0)
return false;

float multPercentile{0.0f};

if constexpr (!isMC) {
Expand Down
Loading