Skip to content

[PWGLF] Change Generated Event selection in NucleiSpectra.cxx #15158

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
Cristian-Moscatelli wants to merge 2 commits into
base: master
Choose a base branch
from
+58 −4
Open

[PWGLF] Change Generated Event selection in NucleiSpectra.cxx #15158

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
9ea81c0
Change Gen Event selection
Feb 25, 2026
45987a4
Clang-format correction
Feb 25, 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
62 changes: 58 additions & 4 deletions PWGLF/TableProducer/Nuspex/nucleiSpectra.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,8 +52,10 @@
#include "ReconstructionDataFormats/Track.h"

#include <Math/Vector4D.h>
#include <TDatabasePDG.h>
#include <TMCProcess.h>
#include <TPDGCode.h> // for PDG codes
#include <TParticlePDG.h>
#include <TRandom3.h>

#include <algorithm>
Expand Down Expand Up @@ -334,6 +336,7 @@ struct nucleiSpectra {
ConfigurableAxis cfgNTPCClusBins{"cfgNTPCClusBins", {3, 89.5, 159.5}, "N TPC clusters binning"};

Configurable<bool> cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Skimmed dataset processing"};
Configurable<std::string> cfgTriggerList{"cfgTriggerList", "fHe", "Trigger List"};

// running variables for track tuner
o2::dataformats::DCA mDcaInfoCov;
Expand Down Expand Up @@ -453,7 +456,7 @@ struct nucleiSpectra {
return;
}
if (cfgSkimmedProcessing) {
zorro.initCCDB(ccdb.service, bc.runNumber(), bc.timestamp(), "fHe");
zorro.initCCDB(ccdb.service, bc.runNumber(), bc.timestamp(), cfgTriggerList);
zorro.populateHistRegistry(spectra, bc.runNumber());
}
auto timestamp = bc.timestamp();
Expand Down Expand Up @@ -911,21 +914,72 @@ struct nucleiSpectra {
PROCESS_SWITCH(nucleiSpectra, processDataFlowAlternative, "Data analysis with flow - alternative framework", false);

Preslice<TrackCandidates> tracksPerCollisions = aod::track::collisionId;
Preslice<aod::McParticles> particlesPerMcCollision = aod::mcparticle::mcCollisionId;
void processMC(soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels> const& collisions, aod::McCollisions const& mcCollisions, soa::Join<TrackCandidates, aod::McTrackLabels> const& tracks, aod::McParticles const& particlesMC, aod::BCsWithTimestamps const&)
{
nuclei::candidates.clear();

bool selectINELgt0 = cfgEventSelections->get(nuclei::evSel::kINELgt0);
std::vector<bool> goodCollisions(mcCollisions.size(), false);

auto* pdgDB = TDatabasePDG::Instance(); // Useful for evaluating the particle charge

for (const auto& c : mcCollisions) {
spectra.fill(HIST("hGenVtxZ"), c.posZ());

// Apply the |z| < 10 cm condition
if (std::abs(c.posZ()) > 10.f)
Copy link
Collaborator

@mpuccio mpuccio Feb 25, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You should keep track of the events you discard, also note that for the TVX definition has a different Z vertex window

continue;

const auto& slicedParticles = particlesMC.sliceBy(particlesPerMcCollision, c.globalIndex());

bool hasHitFT0A(false);
bool hasHitFT0C(false);
bool acceptEvent = !selectINELgt0;
Copy link
Collaborator

@mpuccio mpuccio Feb 25, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This, in combination with line 950, makes the configurable not useful as the statement in line 950 will be always true.


for (const auto& p : slicedParticles) {
if (!p.isPhysicalPrimary())
continue;

// Apply the TVX trigger condition
if (p.eta() > 3.5f && p.eta() < 4.9f)
Copy link
Collaborator

@mpuccio mpuccio Feb 25, 2026

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

You should request that the particle in acceptance is charged

hasHitFT0A = true;
else if (p.eta() > -3.3f && p.eta() < -2.1f)
hasHitFT0C = true;

// Apply the INEL>0 selection (only in case of active configurable)
if (selectINELgt0 && !acceptEvent) {
if (std::abs(p.eta()) < 1.0f) {
auto* pdg = pdgDB->GetParticle(p.pdgCode());
if (pdg && pdg->Charge() != 0)
acceptEvent = true;
}
}

if (hasHitFT0A && hasHitFT0C && acceptEvent)
break;
}

if (hasHitFT0A && hasHitFT0C && acceptEvent) {
goodCollisions[c.globalIndex()] = true;
spectra.fill(HIST("hGenVtxZ"), c.posZ());
}
}
std::vector<bool> goodCollisions(mcCollisions.size(), false);

for (const auto& collision : collisions) {
if (!eventSelectionWithHisto(collision)) {
continue;
}
goodCollisions[collision.mcCollisionId()] = true;

int mcId = collision.mcCollisionId();
if (mcId < 0)
continue;
if (!goodCollisions[mcId])
continue;

const auto& slicedTracks = tracks.sliceBy(tracksPerCollisions, collision.globalIndex());
fillDataInfo(collision, slicedTracks);
}

std::vector<bool> isReconstructed(particlesMC.size(), false);
for (auto& c : nuclei::candidates) {
auto label = tracks.iteratorAt(c.globalIndex);
Expand Down
Loading