|
24 | 24 | #include "PWGJE/DataModel/Jet.h" |
25 | 25 | #include "PWGJE/DataModel/JetReducedData.h" |
26 | 26 | #include "PWGLF/DataModel/LFStrangenessTables.h" |
| 27 | +#include "PWGLF/DataModel/mcCentrality.h" |
27 | 28 |
|
28 | 29 | #include "Common/Core/RecoDecay.h" |
29 | 30 | #include "Common/Core/trackUtilities.h" |
@@ -74,7 +75,7 @@ using std::array; |
74 | 75 |
|
75 | 76 | // Define convenient aliases for joined AOD tables |
76 | 77 | using SelCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>; |
77 | | -using SimCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::McCollisionLabels>; |
| 78 | +using SimCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::McCentFT0Ms, aod::McCollisionLabels>; |
78 | 79 | using DaughterTracks = soa::Join<aod::Tracks, aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::TracksDCA, |
79 | 80 | aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr, |
80 | 81 | aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>; |
@@ -1157,7 +1158,7 @@ struct StrangenessInJets { |
1157 | 1158 | Preslice<DaughterTracksMC> perCollisionTrk = o2::aod::track::collisionId; |
1158 | 1159 |
|
1159 | 1160 | // Generated MC events |
1160 | | - void processMCgenerated(aod::McCollisions const& collisions, aod::McParticles const& mcParticles) |
| 1161 | + void processMCgenerated(soa::Join<aod::McCollisions, aod::McCentFT0Ms> const& collisions, aod::McParticles const& mcParticles) |
1161 | 1162 | { |
1162 | 1163 | // Define per-event particle containers |
1163 | 1164 | std::vector<fastjet::PseudoJet> fjParticles; |
@@ -1190,7 +1191,7 @@ struct StrangenessInJets { |
1190 | 1191 | registryMC.fill(HIST("number_of_events_mc_gen"), 2.5); |
1191 | 1192 |
|
1192 | 1193 | // Multiplicity of generated event |
1193 | | - double genMultiplicity = 0.0; |
| 1194 | + double genMultiplicity = collision.centFT0M(); |
1194 | 1195 |
|
1195 | 1196 | // MC particles per collision |
1196 | 1197 | auto mcParticlesPerColl = mcParticles.sliceBy(perMCCollision, collision.globalIndex()); |
|
0 commit comments