From be96d8ce5af82a331bbc122f6f305d27fc6ef825 Mon Sep 17 00:00:00 2001
From: Marcello Di Costanzo <mdicosta@aliceml.cern.ch>
Date: Wed, 11 Nov 2026 16:08:55 +0100
Subject: [PATCH 1/2] Use generated kinematics in correlations efficiencies

---
 .../TableProducer/correlatorFlowCharmHadronsReduced.cxx  | 1 +
 PWGHF/HFC/Tasks/taskCorrelationsEfficiencies.cxx         | 9 +++++----
 2 files changed, 6 insertions(+), 4 deletions(-)

diff --git a/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx b/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
index 07db185bd39..84c560a3de5 100644
--- a/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
+++ b/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
@@ -46,6 +46,7 @@
 
 #include <chrono>
 #include <cmath>
+#include <chrono>
 #include <cstddef>
 #include <cstdint>
 #include <cstdlib>
diff --git a/PWGHF/HFC/Tasks/taskCorrelationsEfficiencies.cxx b/PWGHF/HFC/Tasks/taskCorrelationsEfficiencies.cxx
index 351d8897a29..a5bb2e8cd66 100644
--- a/PWGHF/HFC/Tasks/taskCorrelationsEfficiencies.cxx
+++ b/PWGHF/HFC/Tasks/taskCorrelationsEfficiencies.cxx
@@ -416,7 +416,7 @@ struct HfTaskCorrelationsEfficiencies {
         }
         int assTrkOrigin = mcParticle.isPhysicalPrimary() ? 0 : 1; // 0 for primary, 1 for secondary
 
-        std::vector<double> values = {assTrk.pt(), assTrk.eta(), assTrk.phi(), mcCollision.posZ(), static_cast<double>(matchedSpecie), static_cast<double>(assTrkOrigin)};
+        std::vector<double> values = {mcParticle.pt(), mcParticle.eta(), mcParticle.phi(), mcCollision.posZ(), static_cast<double>(matchedSpecie), static_cast<double>(assTrkOrigin)};
 
         // Fill histograms for associated tracks
         registry.get<THnSparse>(HIST("hAssTrackReco"))->Fill(values.data());
@@ -424,12 +424,13 @@ struct HfTaskCorrelationsEfficiencies {
         if (fillRecoTrees) {
           // Downsample tracks based on pt
           if (downSampleTracks < 1.) {
-            float const pseudoRndm = assTrk.pt() * 1000. - static_cast<int64_t>(assTrk.pt() * 1000);
-            if (assTrk.pt() < ptMaxDownSampleTracks && pseudoRndm >= downSampleTracks) {
+            // Generate a pseudo-random number between 0 and 1
+            float const pseudoRndm = mcParticle.pt() * 1000. - static_cast<int64_t>(mcParticle.pt() * 1000);
+            if (mcParticle.pt() < ptMaxDownSampleTracks && pseudoRndm >= downSampleTracks) {
               continue;
             }
           }
-          rowAssocMc(assTrk.pt(), assTrk.eta(), assTrk.phi(), mcCollision.posZ(), matchedSpecie, assTrkOrigin);
+          rowAssocMc(mcParticle.pt(), mcParticle.eta(), mcParticle.phi(), mcCollision.posZ(), matchedSpecie, assTrkOrigin);
           rowAssocTrkSels(assTrk.tpcNClsCrossedRows(), assTrk.itsClusterMap(), assTrk.itsNCls(), assTrk.dcaXY(), assTrk.dcaZ());
         }
       }

From 20f702be4b94cb73399f1313a71257787d1a0dcf Mon Sep 17 00:00:00 2001
From: Marcello Di Costanzo <mdicosta@aliceml.cern.ch>
Date: Sat, 14 Nov 2026 21:57:28 +0100
Subject: [PATCH 2/2] Remove header

---
 PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx | 1 -
 1 file changed, 1 deletion(-)

diff --git a/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx b/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
index 84c560a3de5..07db185bd39 100644
--- a/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
+++ b/PWGHF/HFC/TableProducer/correlatorFlowCharmHadronsReduced.cxx
@@ -46,7 +46,6 @@
 
 #include <chrono>
 #include <cmath>
-#include <chrono>
 #include <cstddef>
 #include <cstdint>
 #include <cstdlib>