Skip to content

[PWGLF] Adjusted generated MC event counter for consistency with reconstructed MC #12012

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
alcaliva merged 1 commit into from
Jul 11, 2025
Merged

[PWGLF] Adjusted generated MC event counter for consistency with reconstructed MC #12012

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
8 changes: 5 additions & 3 deletions PWGLF/Tasks/Strangeness/strangenessInJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -982,12 +982,14 @@ struct StrangenessInJets {
registryMC.fill(HIST("number_of_events_mc_gen"), 0.5);

// Need to apply event selection to simulated events
registryMC.fill(HIST("number_of_events_mc_gen"), 1.5);

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

// Fill event counter after selection on z-vertex
registryMC.fill(HIST("number_of_events_mc_gen"), 1.5);
registryMC.fill(HIST("number_of_events_mc_gen"), 2.5);

// Multiplicity of generated event
double genMultiplicity = 0.0;
Expand All @@ -1011,7 +1013,7 @@ struct StrangenessInJets {
// Skip events with no particles
if (fjParticles.size() < 1)
continue;
registryMC.fill(HIST("number_of_events_mc_gen"), 2.5);
registryMC.fill(HIST("number_of_events_mc_gen"), 3.5);

// Cluster MC particles into jets using anti-kt algorithm
fastjet::JetDefinition jetDef(fastjet::antikt_algorithm, rJet);
Expand All @@ -1036,7 +1038,7 @@ struct StrangenessInJets {
// Apply jet pT threshold
if (jetMinusBkg.pt() < minJetPt)
continue;
registryMC.fill(HIST("number_of_events_mc_gen"), 3.5);
registryMC.fill(HIST("number_of_events_mc_gen"), 4.5);

// Set up two perpendicular cone axes for underlying event estimation
TVector3 jetAxis(jet.px(), jet.py(), jet.pz());
Expand Down
Loading