From 0c701bc335f6c53de8567f7131b9eb91baa39d68 Mon Sep 17 00:00:00 2001
From: alcaliva <32872606+alcaliva@users.noreply.github.com>
Date: Mon, 18 May 2026 16:32:44 +0200
Subject: [PATCH 1/2] Replace cone radius calculation with rJet

---
 PWGLF/Tasks/Strangeness/strangenessInJets.cxx | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
index 1be5e3c83d6..b6fadbaa314 100644
--- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
+++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
@@ -1455,7 +1455,7 @@ struct StrangenessInJets {
 
         // Set up two perpendicular cone axes for underlying event estimation
         const TVector3 jetAxis(jet.px(), jet.py(), jet.pz());
-        const double coneRadius = std::sqrt(jet.area() / PI);
+        const double coneRadius = rJet;
         TVector3 ueAxis1(0, 0, 0), ueAxis2(0, 0, 0);
         getPerpendicularDirections(jetAxis, ueAxis1, ueAxis2);
         if (ueAxis1.Mag() == 0 || ueAxis2.Mag() == 0) {

From 273943ad4f91f770399f02102bc91deae7943e58 Mon Sep 17 00:00:00 2001
From: alcaliva <32872606+alcaliva@users.noreply.github.com>
Date: Mon, 18 May 2026 16:44:24 +0200
Subject: [PATCH 2/2] Refactor pseudorapidity threshold to a constant

---
 PWGLF/Tasks/Strangeness/strangenessInJets.cxx | 24 ++++++++++---------
 1 file changed, 13 insertions(+), 11 deletions(-)

diff --git a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
index b6fadbaa314..245ce30e294 100644
--- a/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
+++ b/PWGLF/Tasks/Strangeness/strangenessInJets.cxx
@@ -1664,7 +1664,8 @@ struct StrangenessInJets {
           if (!motherPos.isPhysicalPrimary())
             continue;
 
-          if (std::abs(motherPos.eta()) > 0.8)
+          const double maxPseudorap = 0.8;
+          if (std::abs(motherPos.eta()) > maxPseudorap)
             continue;
 
           // K0s
@@ -1744,8 +1745,8 @@ struct StrangenessInJets {
           if (!mcParticle.isPhysicalPrimary()) {
             continue;
           }
-
-          if (std::abs(mcParticle.eta()) > 0.8) {
+          const double maxPseudorap = 0.8;
+          if (std::abs(mcParticle.eta()) > maxPseudorap) {
             continue;
           }
 
@@ -1786,9 +1787,10 @@ struct StrangenessInJets {
         }
       }
 
-      for (auto& particle : mcParticlesPerColl) {
+      const double maxPseudorap = 0.8;
+      for (auto const& particle : mcParticlesPerColl) {
 
-        if (particle.isPhysicalPrimary() && std::abs(particle.eta()) <= 0.8) {
+        if (particle.isPhysicalPrimary() && std::abs(particle.eta()) <= maxPseudorap) {
           switch (particle.pdgCode()) {
             case kK0Short:
               if (enabledSignals.value[ParticleOfInterest::kV0Particles]) {
@@ -1889,8 +1891,8 @@ struct StrangenessInJets {
             continue;
           if (!motherPos.isPhysicalPrimary())
             continue;
-
-          if (std::abs(motherPos.eta()) > 0.8)
+          double maxPseudorap = 0.8;
+          if (std::abs(motherPos.eta()) > maxPseudorap)
             continue;
 
           // K0s
@@ -1908,9 +1910,9 @@ struct StrangenessInJets {
         }
       }
 
-      for (auto& particle : mcParticlesPerColl) {
-
-        if (particle.isPhysicalPrimary() && std::abs(particle.eta()) <= 0.8) {
+      for (auto const& particle : mcParticlesPerColl) {
+        double maxPseudorap = 0.8;
+        if (particle.isPhysicalPrimary() && std::abs(particle.eta()) <= maxPseudorap) {
           switch (particle.pdgCode()) {
             case kK0Short:
               if (enabledSignals.value[ParticleOfInterest::kV0Particles]) {
@@ -1936,7 +1938,7 @@ struct StrangenessInJets {
       // Loop over selected jets
       for (int i = 0; i < static_cast<int>(selectedJet.size()); i++) {
 
-        for (auto& particle : mcParticlesPerColl) {
+        for (auto const& particle : mcParticlesPerColl) {
 
           const ParticlePositionWithRespectToJet positionMC{particle.px(), particle.py(), particle.pz(), selectedJet[i], ue1[i], ue2[i]};