From dd0bf2c1371e99d52f95d837a59d757f3846e876 Mon Sep 17 00:00:00 2001
From: Yakiv <yakiv.paroviak@cern.ch>
Date: Fri, 21 Nov 2025 12:50:04 +0100
Subject: [PATCH] Fix Purity interpolation

---
 PWGLF/Tasks/Strangeness/strangeCascTrack.cxx | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/PWGLF/Tasks/Strangeness/strangeCascTrack.cxx b/PWGLF/Tasks/Strangeness/strangeCascTrack.cxx
index 45bea61a7c6..2f03174760e 100644
--- a/PWGLF/Tasks/Strangeness/strangeCascTrack.cxx
+++ b/PWGLF/Tasks/Strangeness/strangeCascTrack.cxx
@@ -490,11 +490,11 @@ struct StrangeCascTrack {
 
       if (doApplyPurity1D) {
         if constexpr (requires { cascade.topologyChi2(); }) {
-          purityOmega = hPurityOmegaTra1D->Interpolate(cascade.pt(), mult);
-          purityXi = hPurityXiTra1D->Interpolate(cascade.pt(), mult);
+          purityOmega = hPurityOmegaTra1D->Interpolate(cascade.pt());
+          purityXi = hPurityXiTra1D->Interpolate(cascade.pt());
           if (doPropagatePurity1D) {
-            purityOmegaErr = hPurityErrOmegaTra1D->Interpolate(cascade.pt(), mult);
-            purityXiErr = hPurityErrXiTra1D->Interpolate(cascade.pt(), mult);
+            purityOmegaErr = hPurityErrOmegaTra1D->Interpolate(cascade.pt());
+            purityXiErr = hPurityErrXiTra1D->Interpolate(cascade.pt());
           }
           if (purityOmega == 0) { // check for zero purity, do not apply if the case
             purityOmega = 1.;
@@ -505,11 +505,11 @@ struct StrangeCascTrack {
             purityXiErr = 0.;
           }
         } else {
-          purityOmega = hPurityOmegaStd1D->Interpolate(cascade.pt(), mult);
-          purityXi = hPurityXiStd1D->Interpolate(cascade.pt(), mult);
+          purityOmega = hPurityOmegaStd1D->Interpolate(cascade.pt());
+          purityXi = hPurityXiStd1D->Interpolate(cascade.pt());
           if (doPropagatePurity1D) {
-            purityOmegaErr = hPurityErrOmegaStd1D->Interpolate(cascade.pt(), mult);
-            purityXiErr = hPurityErrXiStd1D->Interpolate(cascade.pt(), mult);
+            purityOmegaErr = hPurityErrOmegaStd1D->Interpolate(cascade.pt());
+            purityXiErr = hPurityErrXiStd1D->Interpolate(cascade.pt());
           }
           if (purityOmega == 0) { // check for zero purity, do not apply if the case
             purityOmega = 1.;