diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx index 057ed5a2638..3bf011a7220 100644 --- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx +++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx @@ -75,7 +75,7 @@ using std::array; // Define convenient aliases for joined AOD tables using SelCollisions = soa::Join; -using SimCollisions = soa::Join; +using SimCollisions = soa::Join; using DaughterTracks = soa::Join; @@ -1429,9 +1429,9 @@ struct StrangenessInJets { PROCESS_SWITCH(StrangenessInJets, processMCgenerated, "process generated events", false); // Reconstructed MC events - void processMCreconstructed(SimCollisions const& collisions, DaughterTracksMC const& mcTracks, - aod::V0Datas const& fullV0s, aod::CascDataExt const& Cascades, - const aod::McParticles&) + void processMCreconstructed(SimCollisions const& collisions, soa::Join const&, + DaughterTracksMC const& mcTracks, aod::V0Datas const& fullV0s, + aod::CascDataExt const& Cascades, const aod::McParticles&) { // Define per-event containers std::vector fjParticles; @@ -1446,6 +1446,12 @@ struct StrangenessInJets { // Loop over reconstructed collisions for (const auto& collision : collisions) { + if (!collision.has_mcCollision()) { + continue; + } + + const auto& mcCollision = collision.mcCollision_as>(); + // Clear containers at the start of the event loop fjParticles.clear(); selectedJet.clear(); @@ -1466,7 +1472,7 @@ struct StrangenessInJets { registryMC.fill(HIST("number_of_events_mc_rec"), 2.5); // Event multiplicity - const float multiplicity = collision.centFT0M(); + const float multiplicity = mcCollision.centFT0M(); // Number of V0 and cascades per collision auto v0sPerColl = fullV0s.sliceBy(perCollisionV0, collision.globalIndex());