Skip to content

Commit b4a58bc

Browse files
authored
[PWGLF] Adjusted generated MC event counter for consistency with reconstructed MC (#12012)
1 parent 109c012 commit b4a58bc

File tree

1 file changed

+5
-3
lines changed

1 file changed

+5
-3
lines changed

PWGLF/Tasks/Strangeness/strangenessInJets.cxx

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -982,12 +982,14 @@ struct StrangenessInJets {
982982
registryMC.fill(HIST("number_of_events_mc_gen"), 0.5);
983983

984984
// Need to apply event selection to simulated events
985+
registryMC.fill(HIST("number_of_events_mc_gen"), 1.5);
986+
985987
// Require vertex position within the allowed z range
986988
if (std::fabs(collision.posZ()) > zVtx)
987989
continue;
988990

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

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

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

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

0 commit comments

Comments
 (0)