[PWGLF] Corrected centrality calibration in RecoMC process (#15148)

nkaratze · web-flow · commit 13d990ae7db8 · 2026-02-24T21:10:27.000+01:00
diff --git a/PWGLF/Tasks/Strangeness/v0ptinvmassplots.cxx b/PWGLF/Tasks/Strangeness/v0ptinvmassplots.cxx
@@ -378,6 +378,7 @@ struct V0PtInvMassPlots {
     rPtAnalysis.fill(HIST("hNEvents"), 5.5, collision.centFT0M());
     rPtAnalysis.get<TH2>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(6, "cutZVertex");
     if (doisInelGt0 && !(collision.multNTracksPVeta1() > 0)) {
+      // if (doisInelGt0 && !(collision.multMCNParticlesEta10() > 0)) { //CHANGE TO THIS
       return false;
     }
     rPtAnalysis.fill(HIST("hNEvents"), 6.5, collision.centFT0M());
@@ -694,17 +695,18 @@ struct V0PtInvMassPlots {
   using DaughterTracksDerived = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;
   o2::framework::Service<o2::framework::O2DatabasePDG> pdgDB;
 
-  void genMCProcess(
-    soa::Join<aod::McCollisions, aod::McCentFT0Ms>::iterator const& mcCollision,
-    soa::SmallGroups<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>> const& collisions,
-    aod::McParticles const& mcParticles)
+  // This is the process for Generated Particles
+  void genMCProcess(soa::Join<aod::McCollisions, aod::McCentFT0Ms, aod::MultMCExtras>::iterator const& mcCollision,
+                    soa::SmallGroups<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>> const& collisions,
+                    aod::McParticles const& mcParticles)
   {
     // Event Efficiency, Event Split and V0 Signal Loss Corrections
     rMCCorrections.fill(HIST("hNEvents_Corrections"), 0.5, mcCollision.centFT0M()); // All Events
     if (std::abs(mcCollision.posZ()) > cutZVertex) {
       return;
     }
-    if (!pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) {
+    if (!(mcCollision.multMCNParticlesEta10() > 0)) { // TRY TO CHANGE TO THIS
+      // if (!pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) {
       return;
     }
     rMCCorrections.fill(HIST("hNEvents_Corrections"), 1.5, mcCollision.centFT0M()); // Event Efficiency Denominator
@@ -787,7 +789,8 @@ struct V0PtInvMassPlots {
     // End of Signal Loss Numenator Loop
   }
   // This is the Process for the MC reconstructed Data
-  void recMCProcess(soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>::iterator const& collision,
+  void recMCProcess(soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::MultsExtraMC, aod::CentFT0Ms>::iterator const& collision,
+                    soa::Join<aod::McCollisions, aod::McCentFT0Ms> const& /*mcCollisions*/,
                     soa::Join<aod::V0Datas, aod::McV0Labels> const& V0s,
                     DaughterTracks const&, // no need to define a variable for tracks, if we don't access them directly
                     aod::McParticles const& /*mcParticles*/)
@@ -808,6 +811,9 @@ struct V0PtInvMassPlots {
     std::vector<double> lambdaptedgevalues(nLambdaHistograms + 1);
     std::vector<double> antilambdaptedgevalues(nAntilambdaHistograms + 1);
 
+    // For centrality estimation
+    const auto& mcCollision = collision.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();
+
     for (int i = 0; i < nKaonHistograms + 1; i++) {
       kaonptedgevalues[i] = std::stod(pthistos::kaonPtBins[i]);
     }
@@ -820,21 +826,21 @@ struct V0PtInvMassPlots {
     if (!acceptEvent(collision)) { // Event Selection
       return;
     }
-    rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Event Split Numenator
+    rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, mcCollision.centFT0M()); // Event Split Numenator
     for (const auto& v0 : V0s) {
       // Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
       const auto& posDaughterTrack = v0.template posTrack_as<DaughterTracks>();
       const auto& negDaughterTrack = v0.template negTrack_as<DaughterTracks>();
-      if (!acceptV0(v0, posDaughterTrack, negDaughterTrack, collision)) { // V0 Selections
+      if (!acceptV0(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // V0 Selections
         continue;
       }
       // kzero analysis
       if (kzeroAnalysis == true) {
-        if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, collision)) { // K0sh Selection
+        if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // K0sh Selection
           // K0sh Signal Split Numerator Start
           for (int i = 0; i < nKaonHistograms; i++) {
             if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator
-              pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M());       // filling the k0s namespace histograms for K0sh Splitting Numerator
+              pthistos::kaonSplit[i]->Fill(v0.mK0Short(), mcCollision.centFT0M());     // filling the k0s namespace histograms for K0sh Splitting Numerator
             }
           }
           // K0sh Signla Split Numerator End
@@ -844,17 +850,17 @@ struct V0PtInvMassPlots {
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nKaonHistograms; i++) {
                   if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges
-                    pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M());          // filling the k0s namespace histograms
+                    pthistos::kaonPt[i]->Fill(v0.mK0Short(), mcCollision.centFT0M());        // filling the k0s namespace histograms
                   }
                 }
               }
               if (!v0mcParticle.isPhysicalPrimary()) {
                 auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers
                 if (!v0mothers.empty()) {
                   auto& v0mcParticleMother = v0mothers.front(); // First mother
-                  rFeeddownMatrices.fill(HIST("hK0shFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                  rFeeddownMatrices.fill(HIST("hK0shFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   if (v0mcParticleMother.pdgCode() == kPhi) { // Phi Mother Matched
-                    rFeeddownMatrices.fill(HIST("hK0shPhiFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hK0shPhiFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                 }
               }
@@ -864,11 +870,11 @@ struct V0PtInvMassPlots {
       }
       // lambda analysis
       if (lambdaAnalysis == true) {
-        if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Lambda Selections
+        if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // Lambda Selections
           // Lambda Signal Split Numerator Start
           for (int i = 0; i < nLambdaHistograms; i++) {
             if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {
-              pthistos::lambdaSplit[i]->Fill(v0.mLambda(), collision.centFT0M());
+              pthistos::lambdaSplit[i]->Fill(v0.mLambda(), mcCollision.centFT0M());
             }
           }
           // Lambda Signal Split Numerator End
@@ -878,23 +884,23 @@ struct V0PtInvMassPlots {
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nLambdaHistograms; i++) {
                   if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {
-                    pthistos::lambdaPt[i]->Fill(v0.mLambda(), collision.centFT0M());
+                    pthistos::lambdaPt[i]->Fill(v0.mLambda(), mcCollision.centFT0M());
                   }
                 }
               }
               if (!v0mcParticle.isPhysicalPrimary()) {
                 auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers
                 if (!v0mothers.empty()) {
                   auto& v0mcParticleMother = v0mothers.front(); // First mother
-                  rFeeddownMatrices.fill(HIST("hLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                  rFeeddownMatrices.fill(HIST("hLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   if (v0mcParticleMother.pdgCode() == kXiMinus) { // Xi Minus Mother Matched
-                    rFeeddownMatrices.fill(HIST("hLambdaXiMinusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hLambdaXiMinusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                   if (v0mcParticleMother.pdgCode() == kXi0) { // Xi Zero Mother Matched
-                    rFeeddownMatrices.fill(HIST("hLambdaXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hLambdaXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                   if (v0mcParticleMother.pdgCode() == kOmegaMinus) { // Omega Mother Matched
-                    rFeeddownMatrices.fill(HIST("hLambdaOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hLambdaOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                 }
               }
@@ -904,11 +910,11 @@ struct V0PtInvMassPlots {
       }
       // antilambda analysis
       if (antiLambdaAnalysis == true) {
-        if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Antilambda Selections
+        if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // Antilambda Selections
           // Antilambda Signal Split Numerator End
           for (int i = 0; i < nAntilambdaHistograms; i++) {
             if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {
-              pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), collision.centFT0M());
+              pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), mcCollision.centFT0M());
             }
           }
           // Antilambda Signal Split Numerator End
@@ -918,23 +924,23 @@ struct V0PtInvMassPlots {
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nAntilambdaHistograms; i++) {
                   if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {
-                    pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), collision.centFT0M());
+                    pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), mcCollision.centFT0M());
                   }
                 }
               }
               if (!v0mcParticle.isPhysicalPrimary()) {
                 auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers
                 if (!v0mothers.empty()) {
                   auto& v0mcParticleMother = v0mothers.front(); // First mother
-                  rFeeddownMatrices.fill(HIST("hAntiLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                  rFeeddownMatrices.fill(HIST("hAntiLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   if (v0mcParticleMother.pdgCode() == kXiPlusBar) { // Xi Plus Mother Matched
-                    rFeeddownMatrices.fill(HIST("hAntiLambdaXiPlusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hAntiLambdaXiPlusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                   if (v0mcParticleMother.pdgCode() == -kXi0) { // Anti-Xi Zero Mother Matched
-                    rFeeddownMatrices.fill(HIST("hAntiLambdaAntiXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hAntiLambdaAntiXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                   if (v0mcParticleMother.pdgCode() == kOmegaPlusBar) { // Anti-Omega (minus) Mother Matched
-                    rFeeddownMatrices.fill(HIST("hAntiLambdaAntiOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());
+                    rFeeddownMatrices.fill(HIST("hAntiLambdaAntiOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());
                   }
                 }
               }
@@ -945,7 +951,7 @@ struct V0PtInvMassPlots {
     }
   }
   // This is the process for Real Data
-  void dataProcess(soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::CentFT0Ms /*,aod::CentNGlobals*/>::iterator const& collision,
+  void dataProcess(soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::CentFT0Ms, aod::MultsExtra /*,aod::CentNGlobals*/>::iterator const& collision,
                    aod::V0Datas const& V0s,
                    DaughterTracks const&)
   {
@@ -975,9 +981,9 @@ struct V0PtInvMassPlots {
       antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);
     }
 
-    if (!acceptEvent(collision)) { // Event Selection
-      return;
-    }
+    // if (!acceptEvent(collision)) { // Event Selection
+    //   return;
+    // }
     rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Number of recorded events
     for (const auto& v0 : V0s) {
       // Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
@@ -1049,9 +1055,9 @@ struct V0PtInvMassPlots {
       antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);
     }
 
-    if (!acceptEvent(collision)) { // Event Selection
-      return;
-    }
+    // if (!acceptEvent(collision)) { // Event Selection
+    //   return;
+    // }
     rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Number of recorded events
     for (const auto& v0 : V0s) {
       // Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
@@ -1092,10 +1098,9 @@ struct V0PtInvMassPlots {
       }
     }
   }
-  // This is the Process for the MC reconstructed Data
   void recMCProcessDerived(soa::Join<aod::StraCollisions, aod::StraEvSels, aod::StraCollLabels, aod::StraCents>::iterator const& collision,
-                           // <aod::McCollisions> const& /*mcCollisions*/,
-                           soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels, aod::V0MCMothers> const& V0s,
+                           // To add McCentFT0Ms
+                           soa::Join<aod::V0CollRefs, aod::V0MCCores, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels, aod::V0MCMothers> const& V0s,
                            DaughterTracksDerived const&)
   {
     // tokenise strings into individual values
@@ -1123,9 +1128,9 @@ struct V0PtInvMassPlots {
     for (int i = 0; i < nAntilambdaHistograms + 1; i++) {
       antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);
     }
-    if (!acceptEvent(collision)) { // Event Selection
-      return;
-    }
+    // if (!acceptEvent(collision)) { // Event Selection
+    //   return;
+    // }
     rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Event Split Numenator
     for (const auto& v0 : V0s) {
       // Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)
@@ -1139,18 +1144,18 @@ struct V0PtInvMassPlots {
         if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, collision)) { // K0sh Selection
           // K0sh Signal Split Numerator Start
           for (int i = 0; i < nKaonHistograms; i++) {
-            if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator
-              pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M());       // filling the k0s namespace histograms for K0sh Splitting Numerator
+            if (kaonptedgevalues[i] <= v0.ptMC() && v0.ptMC() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator
+              pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M());           // filling the k0s namespace histograms for K0sh Splitting Numerator
             }
           }
-          // K0sh Signla Split Numerator End
+          // K0sh SignaL Split Numerator End
           if (v0.has_v0MCCore()) {
             auto v0mcParticle = v0.v0MCCore_as<aod::V0MCCores>();
             if (dotruthk0sh && (v0mcParticle.pdgCode() == kK0Short)) { // kzero matched
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nKaonHistograms; i++) {
-                  if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges
-                    pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M());          // filling the k0s namespace histograms
+                  if (kaonptedgevalues[i] <= v0.ptMC() && v0.ptMC() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges
+                    pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M());              // filling the k0s namespace histograms
                   }
                 }
               }
@@ -1170,7 +1175,7 @@ struct V0PtInvMassPlots {
         if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Lambda Selections
           // Lambda Signal Split Numerator Start
           for (int i = 0; i < nLambdaHistograms; i++) {
-            if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {
+            if (lambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < lambdaptedgevalues[i + 1]) {
               pthistos::lambdaSplit[i]->Fill(v0.mLambda(), collision.centFT0M());
             }
           }
@@ -1180,7 +1185,7 @@ struct V0PtInvMassPlots {
             if (dotruthLambda && (v0mcParticle.pdgCode() == kLambda0)) { // lambda matched
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nLambdaHistograms; i++) {
-                  if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {
+                  if (lambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < lambdaptedgevalues[i + 1]) {
                     pthistos::lambdaPt[i]->Fill(v0.mLambda(), collision.centFT0M());
                   }
                 }
@@ -1207,7 +1212,7 @@ struct V0PtInvMassPlots {
         if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Antilambda Selections
           // Antilambda Signal Split Numerator End
           for (int i = 0; i < nAntilambdaHistograms; i++) {
-            if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {
+            if (antilambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < antilambdaptedgevalues[i + 1]) {
               pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), collision.centFT0M());
             }
           }
@@ -1217,7 +1222,7 @@ struct V0PtInvMassPlots {
             if (dotruthAntilambda && (v0mcParticle.pdgCode() == kLambda0Bar)) { // antilambda matched
               if (v0mcParticle.isPhysicalPrimary()) {
                 for (int i = 0; i < nAntilambdaHistograms; i++) {
-                  if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {
+                  if (antilambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < antilambdaptedgevalues[i + 1]) {
                     pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), collision.centFT0M());
                   }
                 }

Original file line number	Diff line number	Diff line change
`@@ -378,6 +378,7 @@ struct V0PtInvMassPlots {`
`378`	`378`	`rPtAnalysis.fill(HIST("hNEvents"), 5.5, collision.centFT0M());`
`379`	`379`	`rPtAnalysis.get<TH2>(HIST("hNEvents"))->GetXaxis()->SetBinLabel(6, "cutZVertex");`
`380`	`380`	`if (doisInelGt0 && !(collision.multNTracksPVeta1() > 0)) {`
	`381`	`+ // if (doisInelGt0 && !(collision.multMCNParticlesEta10() > 0)) { //CHANGE TO THIS`
`381`	`382`	`return false;`
`382`	`383`	`}`
`383`	`384`	`rPtAnalysis.fill(HIST("hNEvents"), 6.5, collision.centFT0M());`
`@@ -694,17 +695,18 @@ struct V0PtInvMassPlots {`
`694`	`695`	`using DaughterTracksDerived = soa::Join<aod::DauTrackExtras, aod::DauTrackTPCPIDs>;`
`695`	`696`	`o2::framework::Service<o2::framework::O2DatabasePDG> pdgDB;`
`696`	`697`
`697`		`- void genMCProcess(`
`698`		`- soa::Join<aod::McCollisions, aod::McCentFT0Ms>::iterator const& mcCollision,`
`699`		`- soa::SmallGroups<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>> const& collisions,`
`700`		`- aod::McParticles const& mcParticles)`
	`698`	`+ // This is the process for Generated Particles`
	`699`	`+ void genMCProcess(soa::Join<aod::McCollisions, aod::McCentFT0Ms, aod::MultMCExtras>::iterator const& mcCollision,`
	`700`	`+ soa::SmallGroups<soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>> const& collisions,`
	`701`	`+ aod::McParticles const& mcParticles)`
`701`	`702`	`{`
`702`	`703`	`// Event Efficiency, Event Split and V0 Signal Loss Corrections`
`703`	`704`	`rMCCorrections.fill(HIST("hNEvents_Corrections"), 0.5, mcCollision.centFT0M()); // All Events`
`704`	`705`	`if (std::abs(mcCollision.posZ()) > cutZVertex) {`
`705`	`706`	`return;`
`706`	`707`	`}`
`707`		`- if (!pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) {`
	`708`	`+ if (!(mcCollision.multMCNParticlesEta10() > 0)) { // TRY TO CHANGE TO THIS`
	`709`	`+ // if (!pwglf::isINELgtNmc(mcParticles, 0, pdgDB)) {`
`708`	`710`	`return;`
`709`	`711`	`}`
`710`	`712`	`rMCCorrections.fill(HIST("hNEvents_Corrections"), 1.5, mcCollision.centFT0M()); // Event Efficiency Denominator`
`@@ -787,7 +789,8 @@ struct V0PtInvMassPlots {`
`787`	`789`	`// End of Signal Loss Numenator Loop`
`788`	`790`	`}`
`789`	`791`	`// This is the Process for the MC reconstructed Data`
`790`		`- void recMCProcess(soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::McCentFT0Ms>::iterator const& collision,`
	`792`	`+ void recMCProcess(soa::Join<aod::Collisions, aod::EvSels, aod::McCollisionLabels, aod::PVMults, aod::MultsExtraMC, aod::CentFT0Ms>::iterator const& collision,`
	`793`	`+ soa::Join<aod::McCollisions, aod::McCentFT0Ms> const& /mcCollisions/,`
`791`	`794`	`soa::Join<aod::V0Datas, aod::McV0Labels> const& V0s,`
`792`	`795`	`DaughterTracks const&, // no need to define a variable for tracks, if we don't access them directly`
`793`	`796`	`aod::McParticles const& /mcParticles/)`
`@@ -808,6 +811,9 @@ struct V0PtInvMassPlots {`
`808`	`811`	`std::vector<double> lambdaptedgevalues(nLambdaHistograms + 1);`
`809`	`812`	`std::vector<double> antilambdaptedgevalues(nAntilambdaHistograms + 1);`
`810`	`813`
	`814`	`+ // For centrality estimation`
	`815`	`+ const auto& mcCollision = collision.mcCollision_as<soa::Join<aod::McCollisions, aod::McCentFT0Ms>>();`
	`816`	`+`
`811`	`817`	`for (int i = 0; i < nKaonHistograms + 1; i++) {`
`812`	`818`	`kaonptedgevalues[i] = std::stod(pthistos::kaonPtBins[i]);`
`813`	`819`	`}`
`@@ -820,21 +826,21 @@ struct V0PtInvMassPlots {`
`820`	`826`	`if (!acceptEvent(collision)) { // Event Selection`
`821`	`827`	`return;`
`822`	`828`	`}`
`823`		`- rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Event Split Numenator`
	`829`	`+ rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, mcCollision.centFT0M()); // Event Split Numenator`
`824`	`830`	`for (const auto& v0 : V0s) {`
`825`	`831`	`// Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)`
`826`	`832`	`const auto& posDaughterTrack = v0.template posTrack_as<DaughterTracks>();`
`827`	`833`	`const auto& negDaughterTrack = v0.template negTrack_as<DaughterTracks>();`
`828`		`- if (!acceptV0(v0, posDaughterTrack, negDaughterTrack, collision)) { // V0 Selections`
	`834`	`+ if (!acceptV0(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // V0 Selections`
`829`	`835`	`continue;`
`830`	`836`	`}`
`831`	`837`	`// kzero analysis`
`832`	`838`	`if (kzeroAnalysis == true) {`
`833`		`- if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, collision)) { // K0sh Selection`
	`839`	`+ if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // K0sh Selection`
`834`	`840`	`// K0sh Signal Split Numerator Start`
`835`	`841`	`for (int i = 0; i < nKaonHistograms; i++) {`
`836`	`842`	`if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator`
`837`		`- pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms for K0sh Splitting Numerator`
	`843`	`+ pthistos::kaonSplit[i]->Fill(v0.mK0Short(), mcCollision.centFT0M()); // filling the k0s namespace histograms for K0sh Splitting Numerator`
`838`	`844`	`}`
`839`	`845`	`}`
`840`	`846`	`// K0sh Signla Split Numerator End`
`@@ -844,17 +850,17 @@ struct V0PtInvMassPlots {`
`844`	`850`	`if (v0mcParticle.isPhysicalPrimary()) {`
`845`	`851`	`for (int i = 0; i < nKaonHistograms; i++) {`
`846`	`852`	`if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges`
`847`		`- pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms`
	`853`	`+ pthistos::kaonPt[i]->Fill(v0.mK0Short(), mcCollision.centFT0M()); // filling the k0s namespace histograms`
`848`	`854`	`}`
`849`	`855`	`}`
`850`	`856`	`}`
`851`	`857`	`if (!v0mcParticle.isPhysicalPrimary()) {`
`852`	`858`	`auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers`
`853`	`859`	`if (!v0mothers.empty()) {`
`854`	`860`	`auto& v0mcParticleMother = v0mothers.front(); // First mother`
`855`		`- rFeeddownMatrices.fill(HIST("hK0shFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`861`	`+ rFeeddownMatrices.fill(HIST("hK0shFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`856`	`862`	`if (v0mcParticleMother.pdgCode() == kPhi) { // Phi Mother Matched`
`857`		`- rFeeddownMatrices.fill(HIST("hK0shPhiFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`863`	`+ rFeeddownMatrices.fill(HIST("hK0shPhiFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`858`	`864`	`}`
`859`	`865`	`}`
`860`	`866`	`}`
`@@ -864,11 +870,11 @@ struct V0PtInvMassPlots {`
`864`	`870`	`}`
`865`	`871`	`// lambda analysis`
`866`	`872`	`if (lambdaAnalysis == true) {`
`867`		`- if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Lambda Selections`
	`873`	`+ if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // Lambda Selections`
`868`	`874`	`// Lambda Signal Split Numerator Start`
`869`	`875`	`for (int i = 0; i < nLambdaHistograms; i++) {`
`870`	`876`	`if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {`
`871`		`- pthistos::lambdaSplit[i]->Fill(v0.mLambda(), collision.centFT0M());`
	`877`	`+ pthistos::lambdaSplit[i]->Fill(v0.mLambda(), mcCollision.centFT0M());`
`872`	`878`	`}`
`873`	`879`	`}`
`874`	`880`	`// Lambda Signal Split Numerator End`
`@@ -878,23 +884,23 @@ struct V0PtInvMassPlots {`
`878`	`884`	`if (v0mcParticle.isPhysicalPrimary()) {`
`879`	`885`	`for (int i = 0; i < nLambdaHistograms; i++) {`
`880`	`886`	`if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {`
`881`		`- pthistos::lambdaPt[i]->Fill(v0.mLambda(), collision.centFT0M());`
	`887`	`+ pthistos::lambdaPt[i]->Fill(v0.mLambda(), mcCollision.centFT0M());`
`882`	`888`	`}`
`883`	`889`	`}`
`884`	`890`	`}`
`885`	`891`	`if (!v0mcParticle.isPhysicalPrimary()) {`
`886`	`892`	`auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers`
`887`	`893`	`if (!v0mothers.empty()) {`
`888`	`894`	`auto& v0mcParticleMother = v0mothers.front(); // First mother`
`889`		`- rFeeddownMatrices.fill(HIST("hLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`895`	`+ rFeeddownMatrices.fill(HIST("hLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`890`	`896`	`if (v0mcParticleMother.pdgCode() == kXiMinus) { // Xi Minus Mother Matched`
`891`		`- rFeeddownMatrices.fill(HIST("hLambdaXiMinusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`897`	`+ rFeeddownMatrices.fill(HIST("hLambdaXiMinusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`892`	`898`	`}`
`893`	`899`	`if (v0mcParticleMother.pdgCode() == kXi0) { // Xi Zero Mother Matched`
`894`		`- rFeeddownMatrices.fill(HIST("hLambdaXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`900`	`+ rFeeddownMatrices.fill(HIST("hLambdaXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`895`	`901`	`}`
`896`	`902`	`if (v0mcParticleMother.pdgCode() == kOmegaMinus) { // Omega Mother Matched`
`897`		`- rFeeddownMatrices.fill(HIST("hLambdaOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`903`	`+ rFeeddownMatrices.fill(HIST("hLambdaOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`898`	`904`	`}`
`899`	`905`	`}`
`900`	`906`	`}`
`@@ -904,11 +910,11 @@ struct V0PtInvMassPlots {`
`904`	`910`	`}`
`905`	`911`	`// antilambda analysis`
`906`	`912`	`if (antiLambdaAnalysis == true) {`
`907`		`- if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Antilambda Selections`
	`913`	`+ if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, mcCollision)) { // Antilambda Selections`
`908`	`914`	`// Antilambda Signal Split Numerator End`
`909`	`915`	`for (int i = 0; i < nAntilambdaHistograms; i++) {`
`910`	`916`	`if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {`
`911`		`- pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), collision.centFT0M());`
	`917`	`+ pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), mcCollision.centFT0M());`
`912`	`918`	`}`
`913`	`919`	`}`
`914`	`920`	`// Antilambda Signal Split Numerator End`
`@@ -918,23 +924,23 @@ struct V0PtInvMassPlots {`
`918`	`924`	`if (v0mcParticle.isPhysicalPrimary()) {`
`919`	`925`	`for (int i = 0; i < nAntilambdaHistograms; i++) {`
`920`	`926`	`if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {`
`921`		`- pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), collision.centFT0M());`
	`927`	`+ pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), mcCollision.centFT0M());`
`922`	`928`	`}`
`923`	`929`	`}`
`924`	`930`	`}`
`925`	`931`	`if (!v0mcParticle.isPhysicalPrimary()) {`
`926`	`932`	`auto v0mothers = v0mcParticle.mothers_as<aod::McParticles>(); // Get mothers`
`927`	`933`	`if (!v0mothers.empty()) {`
`928`	`934`	`auto& v0mcParticleMother = v0mothers.front(); // First mother`
`929`		`- rFeeddownMatrices.fill(HIST("hAntiLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`935`	`+ rFeeddownMatrices.fill(HIST("hAntiLambdaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`930`	`936`	`if (v0mcParticleMother.pdgCode() == kXiPlusBar) { // Xi Plus Mother Matched`
`931`		`- rFeeddownMatrices.fill(HIST("hAntiLambdaXiPlusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`937`	`+ rFeeddownMatrices.fill(HIST("hAntiLambdaXiPlusFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`932`	`938`	`}`
`933`	`939`	`if (v0mcParticleMother.pdgCode() == -kXi0) { // Anti-Xi Zero Mother Matched`
`934`		`- rFeeddownMatrices.fill(HIST("hAntiLambdaAntiXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`940`	`+ rFeeddownMatrices.fill(HIST("hAntiLambdaAntiXiZeroFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`935`	`941`	`}`
`936`	`942`	`if (v0mcParticleMother.pdgCode() == kOmegaPlusBar) { // Anti-Omega (minus) Mother Matched`
`937`		`- rFeeddownMatrices.fill(HIST("hAntiLambdaAntiOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), collision.centFT0M());`
	`943`	`+ rFeeddownMatrices.fill(HIST("hAntiLambdaAntiOmegaFeeddownMatrix"), v0mcParticle.pt(), v0mcParticleMother.pt(), mcCollision.centFT0M());`
`938`	`944`	`}`
`939`	`945`	`}`
`940`	`946`	`}`
`@@ -945,7 +951,7 @@ struct V0PtInvMassPlots {`
`945`	`951`	`}`
`946`	`952`	`}`
`947`	`953`	`// This is the process for Real Data`
`948`		`- void dataProcess(soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::CentFT0Ms /,aod::CentNGlobals/>::iterator const& collision,`
	`954`	`+ void dataProcess(soa::Join<aod::Collisions, aod::EvSels, aod::PVMults, aod::CentFT0Ms, aod::MultsExtra /,aod::CentNGlobals/>::iterator const& collision,`
`949`	`955`	`aod::V0Datas const& V0s,`
`950`	`956`	`DaughterTracks const&)`
`951`	`957`	`{`
`@@ -975,9 +981,9 @@ struct V0PtInvMassPlots {`
`975`	`981`	`antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);`
`976`	`982`	`}`
`977`	`983`
`978`		`- if (!acceptEvent(collision)) { // Event Selection`
`979`		`- return;`
`980`		`- }`
	`984`	`+ // if (!acceptEvent(collision)) { // Event Selection`
	`985`	`+ // return;`
	`986`	`+ // }`
`981`	`987`	`rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Number of recorded events`
`982`	`988`	`for (const auto& v0 : V0s) {`
`983`	`989`	`// Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)`
`@@ -1049,9 +1055,9 @@ struct V0PtInvMassPlots {`
`1049`	`1055`	`antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);`
`1050`	`1056`	`}`
`1051`	`1057`
`1052`		`- if (!acceptEvent(collision)) { // Event Selection`
`1053`		`- return;`
`1054`		`- }`
	`1058`	`+ // if (!acceptEvent(collision)) { // Event Selection`
	`1059`	`+ // return;`
	`1060`	`+ // }`
`1055`	`1061`	`rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Number of recorded events`
`1056`	`1062`	`for (const auto& v0 : V0s) {`
`1057`	`1063`	`// Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)`
`@@ -1092,10 +1098,9 @@ struct V0PtInvMassPlots {`
`1092`	`1098`	`}`
`1093`	`1099`	`}`
`1094`	`1100`	`}`
`1095`		`- // This is the Process for the MC reconstructed Data`
`1096`	`1101`	`void recMCProcessDerived(soa::Join<aod::StraCollisions, aod::StraEvSels, aod::StraCollLabels, aod::StraCents>::iterator const& collision,`
`1097`		`- // <aod::McCollisions> const& /mcCollisions/,`
`1098`		`- soa::Join<aod::V0CollRefs, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels, aod::V0MCMothers> const& V0s,`
	`1102`	`+ // To add McCentFT0Ms`
	`1103`	`+ soa::Join<aod::V0CollRefs, aod::V0MCCores, aod::V0Cores, aod::V0Extras, aod::V0CoreMCLabels, aod::V0MCMothers> const& V0s,`
`1099`	`1104`	`DaughterTracksDerived const&)`
`1100`	`1105`	`{`
`1101`	`1106`	`// tokenise strings into individual values`
`@@ -1123,9 +1128,9 @@ struct V0PtInvMassPlots {`
`1123`	`1128`	`for (int i = 0; i < nAntilambdaHistograms + 1; i++) {`
`1124`	`1129`	`antilambdaptedgevalues[i] = std::stod(pthistos::antilambdaPtBins[i]);`
`1125`	`1130`	`}`
`1126`		`- if (!acceptEvent(collision)) { // Event Selection`
`1127`		`- return;`
`1128`		`- }`
	`1131`	`+ // if (!acceptEvent(collision)) { // Event Selection`
	`1132`	`+ // return;`
	`1133`	`+ // }`
`1129`	`1134`	`rPtAnalysis.fill(HIST("hNRecEvents"), 0.5, collision.centFT0M()); // Event Split Numenator`
`1130`	`1135`	`for (const auto& v0 : V0s) {`
`1131`	`1136`	`// Checking that the V0 is a true K0s/Lambdas/Antilambdas and then filling the parameter histograms and the invariant mass plots for different cuts (which are taken from namespace)`
`@@ -1139,18 +1144,18 @@ struct V0PtInvMassPlots {`
`1139`	`1144`	`if (acceptK0sh(v0, posDaughterTrack, negDaughterTrack, collision)) { // K0sh Selection`
`1140`	`1145`	`// K0sh Signal Split Numerator Start`
`1141`	`1146`	`for (int i = 0; i < nKaonHistograms; i++) {`
`1142`		`- if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator`
`1143`		`- pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms for K0sh Splitting Numerator`
	`1147`	`+ if (kaonptedgevalues[i] <= v0.ptMC() && v0.ptMC() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges for K0sh Splitting Numerator`
	`1148`	`+ pthistos::kaonSplit[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms for K0sh Splitting Numerator`
`1144`	`1149`	`}`
`1145`	`1150`	`}`
`1146`		`- // K0sh Signla Split Numerator End`
	`1151`	`+ // K0sh SignaL Split Numerator End`
`1147`	`1152`	`if (v0.has_v0MCCore()) {`
`1148`	`1153`	`auto v0mcParticle = v0.v0MCCore_as<aod::V0MCCores>();`
`1149`	`1154`	`if (dotruthk0sh && (v0mcParticle.pdgCode() == kK0Short)) { // kzero matched`
`1150`	`1155`	`if (v0mcParticle.isPhysicalPrimary()) {`
`1151`	`1156`	`for (int i = 0; i < nKaonHistograms; i++) {`
`1152`		`- if (kaonptedgevalues[i] <= v0.pt() && v0.pt() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges`
`1153`		`- pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms`
	`1157`	`+ if (kaonptedgevalues[i] <= v0.ptMC() && v0.ptMC() < kaonptedgevalues[i + 1]) { // finding v0s with pt within the range of our bin edges`
	`1158`	`+ pthistos::kaonPt[i]->Fill(v0.mK0Short(), collision.centFT0M()); // filling the k0s namespace histograms`
`1154`	`1159`	`}`
`1155`	`1160`	`}`
`1156`	`1161`	`}`
`@@ -1170,7 +1175,7 @@ struct V0PtInvMassPlots {`
`1170`	`1175`	`if (acceptLambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Lambda Selections`
`1171`	`1176`	`// Lambda Signal Split Numerator Start`
`1172`	`1177`	`for (int i = 0; i < nLambdaHistograms; i++) {`
`1173`		`- if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {`
	`1178`	`+ if (lambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < lambdaptedgevalues[i + 1]) {`
`1174`	`1179`	`pthistos::lambdaSplit[i]->Fill(v0.mLambda(), collision.centFT0M());`
`1175`	`1180`	`}`
`1176`	`1181`	`}`
`@@ -1180,7 +1185,7 @@ struct V0PtInvMassPlots {`
`1180`	`1185`	`if (dotruthLambda && (v0mcParticle.pdgCode() == kLambda0)) { // lambda matched`
`1181`	`1186`	`if (v0mcParticle.isPhysicalPrimary()) {`
`1182`	`1187`	`for (int i = 0; i < nLambdaHistograms; i++) {`
`1183`		`- if (lambdaptedgevalues[i] <= v0.pt() && v0.pt() < lambdaptedgevalues[i + 1]) {`
	`1188`	`+ if (lambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < lambdaptedgevalues[i + 1]) {`
`1184`	`1189`	`pthistos::lambdaPt[i]->Fill(v0.mLambda(), collision.centFT0M());`
`1185`	`1190`	`}`
`1186`	`1191`	`}`
`@@ -1207,7 +1212,7 @@ struct V0PtInvMassPlots {`
`1207`	`1212`	`if (acceptAntilambda(v0, posDaughterTrack, negDaughterTrack, collision)) { // Antilambda Selections`
`1208`	`1213`	`// Antilambda Signal Split Numerator End`
`1209`	`1214`	`for (int i = 0; i < nAntilambdaHistograms; i++) {`
`1210`		`- if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {`
	`1215`	`+ if (antilambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < antilambdaptedgevalues[i + 1]) {`
`1211`	`1216`	`pthistos::antilambdaSplit[i]->Fill(v0.mAntiLambda(), collision.centFT0M());`
`1212`	`1217`	`}`
`1213`	`1218`	`}`
`@@ -1217,7 +1222,7 @@ struct V0PtInvMassPlots {`
`1217`	`1222`	`if (dotruthAntilambda && (v0mcParticle.pdgCode() == kLambda0Bar)) { // antilambda matched`
`1218`	`1223`	`if (v0mcParticle.isPhysicalPrimary()) {`
`1219`	`1224`	`for (int i = 0; i < nAntilambdaHistograms; i++) {`
`1220`		`- if (antilambdaptedgevalues[i] <= v0.pt() && v0.pt() < antilambdaptedgevalues[i + 1]) {`
	`1225`	`+ if (antilambdaptedgevalues[i] <= v0.ptMC() && v0.ptMC() < antilambdaptedgevalues[i + 1]) {`
`1221`	`1226`	`pthistos::antilambdaPt[i]->Fill(v0.mAntiLambda(), collision.centFT0M());`
`1222`	`1227`	`}`
`1223`	`1228`	`}`