From a3c91e09d88300b7c1455f3ea535dca35024a3d2 Mon Sep 17 00:00:00 2001
From: Francesca Ercolessi <francesca.ercolessi@cern.ch>
Date: Sat, 1 Nov 2025 10:49:36 +0100
Subject: [PATCH] add centrality in MC

---
 PWGLF/Tasks/Strangeness/strangenessInJets.cxx | 7 ++++---
 1 file changed, 4 insertions(+), 3 deletions(-)

diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
index 4ae80aa2d8f..057ed5a2638 100644
--- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
+++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
@@ -24,6 +24,7 @@
 #include "PWGJE/DataModel/Jet.h"
 #include "PWGJE/DataModel/JetReducedData.h"
 #include "PWGLF/DataModel/LFStrangenessTables.h"
+#include "PWGLF/DataModel/mcCentrality.h"
 
 #include "Common/Core/RecoDecay.h"
 #include "Common/Core/trackUtilities.h"
@@ -74,7 +75,7 @@ using std::array;
 
 // Define convenient aliases for joined AOD tables
 using SelCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms>;
-using SimCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::CentFT0Ms, aod::McCollisionLabels>;
+using SimCollisions = soa::Join<aod::Collisions, aod::EvSels, aod::McCentFT0Ms, aod::McCollisionLabels>;
 using DaughterTracks = soa::Join<aod::Tracks, aod::TracksIU, aod::TracksExtra, aod::TracksCovIU, aod::TracksDCA,
                                  aod::pidTPCFullPi, aod::pidTPCFullKa, aod::pidTPCFullPr,
                                  aod::pidTOFFullPi, aod::pidTOFFullKa, aod::pidTOFFullPr>;
@@ -1157,7 +1158,7 @@ struct StrangenessInJets {
   Preslice<DaughterTracksMC> perCollisionTrk = o2::aod::track::collisionId;
 
   // Generated MC events
-  void processMCgenerated(aod::McCollisions const& collisions, aod::McParticles const& mcParticles)
+  void processMCgenerated(soa::Join<aod::McCollisions, aod::McCentFT0Ms> const& collisions, aod::McParticles const& mcParticles)
   {
     // Define per-event particle containers
     std::vector<fastjet::PseudoJet> fjParticles;
@@ -1190,7 +1191,7 @@ struct StrangenessInJets {
       registryMC.fill(HIST("number_of_events_mc_gen"), 2.5);
 
       // Multiplicity of generated event
-      double genMultiplicity = 0.0;
+      double genMultiplicity = collision.centFT0M();
 
       // MC particles per collision
       auto mcParticlesPerColl = mcParticles.sliceBy(perMCCollision, collision.globalIndex());