Skip to content

[PWGLF] Add process function for data in minimum bias events #16310

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
romainschotter merged 2 commits into from
May 19, 2026
+190 −6
Merged

[PWGLF] Add process function for data in minimum bias events #16310

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
2b7af55
Use long unsigned int in for loop in AddV0sForJetReconstruction
lorber98 May 19, 2026
ec6ab32
Add process function for data in minimum bias events
lorber98 May 19, 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
196 changes: 190 additions & 6 deletions PWGLF/Tasks/Strangeness/strangenessInJetsIons.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,6 +106,7 @@

// Define histogram registries
HistogramRegistry registryData{"registryData", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryDataMB{"registryDataMB", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryMC{"registryMC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};
HistogramRegistry registryQC{"registryQC", {}, OutputObjHandlingPolicy::AnalysisObject, true, true};

Expand Down Expand Up @@ -549,6 +550,45 @@
registryMC.add("h3dAntiLambdaFeeddownFromXi0", "h3dAntiLambdaFeeddownFromXi0", kTH3D, {multAxis, ptAxis, ptAxis});
}
}

// Histograms for real data in MB events
if (doprocessDataMB) {

// Event counters
registryDataMB.add("number_of_events_data", "number of events in data", HistType::kTH1D, {{20, 0, 20, "Event Cuts"}});
registryDataMB.add("number_of_events_vsmultiplicity", "number of events in data vs multiplicity", HistType::kTH1D, {{101, 0, 101, "Multiplicity percentile"}});

registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(1, "All collisions");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(2, "Zorro selection");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(3, "sel8");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(4, "posZ cut");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(5, "kNoSameBunchPileup");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(6, "kIsGoodZvtxFT0vsPV");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(7, "No empty events");
registryDataMB.get<TH1>(HIST("number_of_events_data"))->GetXaxis()->SetBinLabel(8, "At least one jet");

// Histograms for analysis of strange hadrons
if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) {
registryDataMB.add("Lambda_in_MB", "Lambda_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassLambdaAxis});
registryDataMB.add("AntiLambda_in_MB", "AntiLambda_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassLambdaAxis});
registryDataMB.add("K0s_in_MB", "K0s_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassK0sAxis});
}
if (particleOfInterestDict[ParticleOfInterest::kCascades]) {
registryDataMB.add("XiPos_in_MB", "XiPos_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassXiAxis});
registryDataMB.add("XiNeg_in_MB", "XiNeg_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassXiAxis});
registryDataMB.add("OmegaPos_in_MB", "OmegaPos_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassOmegaAxis});
registryDataMB.add("OmegaNeg_in_MB", "OmegaNeg_in_MB", HistType::kTH3F, {multAxis, ptAxis, invMassOmegaAxis});
}
if (particleOfInterestDict[ParticleOfInterest::kPions]) {
registryDataMB.add("Pion_in_MB", "Pion_in_MB", HistType::kTHnSparseF, {multAxis, ptAxisLongLived, nsigmaTPCAxis, nsigmaTOFAxis, dcaAxis});
}
if (particleOfInterestDict[ParticleOfInterest::kKaons]) {
registryDataMB.add("Kaon_in_MB", "Kaon_in_MB", HistType::kTHnSparseF, {multAxis, ptAxisLongLived, nsigmaTPCAxis, nsigmaTOFAxis, dcaAxis});
}
if (particleOfInterestDict[ParticleOfInterest::kProtons]) {
registryDataMB.add("Proton_in_MB", "Proton_in_MB", HistType::kTHnSparseF, {multAxis, ptAxisLongLived, nsigmaTPCAxis, nsigmaTOFAxis, dcaAxis});
}
}
}

// Delta phi calculation
Expand Down Expand Up @@ -1318,7 +1358,7 @@
if (motherPos != motherNeg)
continue;

if (std::abs(motherPos.eta()) > 0.8)

Check failure on line 1361 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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.
continue;

// Vertex position vector
Expand Down Expand Up @@ -1417,7 +1457,7 @@
if (!motherBach.isPhysicalPrimary())
continue;

if (std::abs(motherPos.eta()) > 0.8)

Check failure on line 1460 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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.
continue;

// Xi+
Expand Down Expand Up @@ -1455,7 +1495,7 @@
continue;
}

if (std::abs(mcParticle.eta()) > 0.8) {

Check failure on line 1498 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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.
continue;
}

Expand Down Expand Up @@ -1512,7 +1552,7 @@
if (std::abs(v0mother.pdgCode()) == o2::constants::physics::Pdg::kXi0)
rapidityXi = RecoDecay::y(std::array{v0mother.px(), v0mother.py(), v0mother.pz()}, o2::constants::physics::MassXi0);

if (std::fabs(rapidityXi) > 0.5f)

Check failure on line 1555 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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; // not a valid mother rapidity (PDG selection is later)

// __________________________________________
Expand Down Expand Up @@ -1635,7 +1675,7 @@
bool isV0 = isK0S || isLambda || isAntiLambda;
if (!isV0)
continue;
for (int i = 0; i < int(fjTracks.size()); ++i) {
for (long unsigned int i = 0; i < fjTracks.size(); ++i) {
if (isV0DaughterTrack(fjTracks[i], v0)) {
// LOG(info) << "[AddV0sForJetReconstructionData] V0 daughter track found in fjTracks.";
isTrackReplaced[i] = true;
Expand All @@ -1645,7 +1685,7 @@

std::vector<fastjet::PseudoJet> cleanFjInput;
cleanFjInput.reserve(fjInput.size());
for (int i = 0; i < int(fjInput.size()); ++i) {
for (long unsigned int i = 0; i < fjInput.size(); ++i) {
if (!isTrackReplaced[i])
cleanFjInput.push_back(fjInput[i]);
}
Expand Down Expand Up @@ -1733,7 +1773,7 @@
bool isV0 = isK0S || isLambda || isAntiLambda;
if (!isV0)
continue;
for (int i = 0; i < int(fjTracks.size()); ++i) {
for (long unsigned int i = 0; i < fjTracks.size(); ++i) {
if (isV0DaughterTrack(fjTracks[i], v0)) {
// LOG(info) << "[AddV0sForJetReconstructionMCD] V0 daughter track found in fjTracks.";
isTrackReplaced[i] = true;
Expand All @@ -1743,7 +1783,7 @@

std::vector<fastjet::PseudoJet> cleanFjInput;
cleanFjInput.reserve(fjInput.size());
for (int i = 0; i < int(fjInput.size()); ++i) {
for (long unsigned int i = 0; i < fjInput.size(); ++i) {
if (!isTrackReplaced[i])
cleanFjInput.push_back(fjInput[i]);
}
Expand Down Expand Up @@ -1789,7 +1829,7 @@
// LOG(info) << "[AddV0sForJetReconstructionMCP] Add V0 as input for jet finder.";

// Remove V0 daughter particles if already in the input list for the jet finder
for (int i = 0; i < int(fjParticleObj.size()); ++i) {
for (long unsigned int i = 0; i < fjParticleObj.size(); ++i) {
const auto& mcPart = fjParticleObj[i];
if (!mcPart.has_mothers())
continue;
Expand All @@ -1805,7 +1845,7 @@

std::vector<fastjet::PseudoJet> cleanFjInput;
cleanFjInput.reserve(fjInput.size());
for (int i = 0; i < int(fjInput.size()); ++i) {
for (long unsigned int i = 0; i < fjInput.size(); ++i) {
if (!isTrackReplaced[i])
cleanFjInput.push_back(fjInput[i]);
}
Expand Down Expand Up @@ -2185,7 +2225,7 @@
}
}

if (multiplicity == -999)

Check failure on line 2228 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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.
continue; */

// Clear containers at the start of the event loop
Expand Down Expand Up @@ -2623,7 +2663,7 @@
// ------------------------------------------------
// --- Generated hadrons in reconstructed jets ----
for (const auto& particle : mcParticlesPerColl) {
if (!particle.isPhysicalPrimary() || std::abs(particle.eta()) > 0.8)

Check failure on line 2666 in PWGLF/Tasks/Strangeness/strangenessInJetsIons.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.
continue;

int absPdg = std::abs(particle.pdgCode());
Expand Down Expand Up @@ -3181,6 +3221,150 @@
}
}
PROCESS_SWITCH(StrangenessInJetsIons, processMCmodels, "Process MC generated events (with models)", false);

// --- Process Minimum Bias events ---
// Process data MB
void processDataMB(SelCollisions::iterator const& collision, aod::V0Datas const& fullV0s,
aod::CascDataExt const& Cascades, DaughterTracks const& tracks,
aod::BCsWithTimestamps const&)
{
// Vertex position vector
TVector3 vtxPos(collision.posX(), collision.posY(), collision.posZ());

// Fill event counter before event selection
registryDataMB.fill(HIST("number_of_events_data"), 0.5);

// Get the bunch crossing (BC) information associated with the collision
auto bc = collision.template bc_as<aod::BCsWithTimestamps>();

// Initialize CCDB objects using the BC info
initCCDB(bc);

// If skimmed processing is enabled, skip this event unless it passes Zorro selection
if (cfgSkimmedProcessing && !zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC())) {
return;
}

// Fill event counter after zorro selection
registryDataMB.fill(HIST("number_of_events_data"), 1.5);

// Event selection
if (!collision.sel8())
return;

// Fill event counter after sel8 selection
registryDataMB.fill(HIST("number_of_events_data"), 2.5);

// Require vertex position within the allowed z range
if (std::fabs(collision.posZ()) > zVtx)
return;

// Fill event counter after z vertex selection
registryDataMB.fill(HIST("number_of_events_data"), 3.5);

// Reject collisions associated to the same found BC
if (requireNoSameBunchPileup && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup))
return;

// Fill event counter after selection kNoSameBunchPileup
registryDataMB.fill(HIST("number_of_events_data"), 4.5);

// Compatible z_vtx from FT0 and from PV
if (requireGoodZvtxFT0vsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
return;

// Fill event counter after selection kIsGoodZvtxFT0vsPV
registryDataMB.fill(HIST("number_of_events_data"), 5.5);

// Event multiplicity
float centrality;
if (centrEstimator == 0) {
centrality = collision.centFT0C();
} else {
centrality = collision.centFT0M();
}
registryDataMB.fill(HIST("number_of_events_vsmultiplicity"), centrality);

if (particleOfInterestDict[ParticleOfInterest::kV0Particles]) { // V0s
for (const auto& v0 : fullV0s) {

// Get V0 daughters
const auto& pos = v0.posTrack_as<DaughterTracks>();
const auto& neg = v0.negTrack_as<DaughterTracks>();
TVector3 v0dir(v0.px(), v0.py(), v0.pz());

// K0s
if (passedK0ShortSelection(v0, pos, neg, vtxPos)) {
registryDataMB.fill(HIST("K0s_in_MB"), centrality, v0.pt(), v0.mK0Short());
}
// Lambda
if (passedLambdaSelection(v0, pos, neg, vtxPos)) {
registryDataMB.fill(HIST("Lambda_in_MB"), centrality, v0.pt(), v0.mLambda());
}
// AntiLambda
if (passedAntiLambdaSelection(v0, pos, neg, vtxPos)) {
registryDataMB.fill(HIST("AntiLambda_in_MB"), centrality, v0.pt(), v0.mAntiLambda());
}
}
}

if (particleOfInterestDict[ParticleOfInterest::kCascades]) { // Cascades
for (const auto& casc : Cascades) {
// Get cascade daughters
const auto& bach = casc.bachelor_as<DaughterTracks>();
const auto& pos = casc.posTrack_as<DaughterTracks>();
const auto& neg = casc.negTrack_as<DaughterTracks>();
TVector3 cascadeDir(casc.px(), casc.py(), casc.pz());

// Xi+
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
registryDataMB.fill(HIST("XiPos_in_MB"), centrality, casc.pt(), casc.mXi());
}
// Xi-
if (passedXiSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
registryDataMB.fill(HIST("XiNeg_in_MB"), centrality, casc.pt(), casc.mXi());
}
// Omega+
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() > 0) {
registryDataMB.fill(HIST("OmegaPos_in_MB"), centrality, casc.pt(), casc.mOmega());
}
// Omega-
if (passedOmegaSelection(casc, pos, neg, bach, collision) && bach.sign() < 0) {
registryDataMB.fill(HIST("OmegaNeg_in_MB"), centrality, casc.pt(), casc.mOmega());
}
}
}
if (particleOfInterestDict[ParticleOfInterest::kPions] ||
particleOfInterestDict[ParticleOfInterest::kKaons] ||
particleOfInterestDict[ParticleOfInterest::kProtons]) {
for (const auto& trk : tracks) {

if (!passedSingleTrackSelection(trk)) {
continue;
}

float nsigmaTPC = 0.f;
float nsigmaTOF = 0.f;

if (particleOfInterestDict[ParticleOfInterest::kPions]) {
nsigmaTPC = trk.tpcNSigmaPi();
nsigmaTOF = trk.tofNSigmaPi();
registryDataMB.fill(HIST("Pion_in_MB"), centrality, trk.pt() * trk.sign(), nsigmaTPC, nsigmaTOF, trk.dcaXY());
}
if (particleOfInterestDict[ParticleOfInterest::kKaons]) {
nsigmaTPC = trk.tpcNSigmaKa();
nsigmaTOF = trk.tofNSigmaKa();
registryDataMB.fill(HIST("Kaon_in_MB"), centrality, trk.pt() * trk.sign(), nsigmaTPC, nsigmaTOF, trk.dcaXY());
}
if (particleOfInterestDict[ParticleOfInterest::kProtons]) {
nsigmaTPC = trk.tpcNSigmaPr();
nsigmaTOF = trk.tofNSigmaPr();
registryDataMB.fill(HIST("Proton_in_MB"), centrality, trk.pt() * trk.sign(), nsigmaTPC, nsigmaTOF, trk.dcaXY());
}
}
}
}
PROCESS_SWITCH(StrangenessInJetsIons, processDataMB, "Process data in minimum bias events", false);
};

WorkflowSpec defineDataProcessing(ConfigContext const& cfgc)
Expand Down
Loading