[PWGLF] NucleiTask - ITS cut fix

giovannimalfattore · Jan 10, 2025 · b5e4525 · b5e4525
1 parent a6ac2a5
commit b5e4525
Showing 1 changed file with 4 additions and 5 deletions.
diff --git a/PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx b/PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx
@@ -2147,6 +2147,8 @@ struct LFNucleiBATask {
       }
 
       float nITSTr, nITSHe;
+      nITSTr = track.itsNSigmaTr();
+      nITSHe = track.itsNSigmaHe();
 
       heP = track.p();
       antiheP = track.p();
@@ -2313,8 +2315,8 @@ struct LFNucleiBATask {
         isHe = isHelium && track.sign() > 0;
         isAntiHe = isHelium && track.sign() < 0;
       } else {
-        isHe = isHelium && track.sign() > 0 && std::abs(track.itsNSigmaTr()) > nsigmaITSvar.nsigmaITSTr;
-        isAntiHe = isHelium && track.sign() < 0 && std::abs(track.itsNSigmaTr()) > nsigmaITSvar.nsigmaITSTr;
+        isHe = isHelium && track.sign() > 0 && std::abs(nITSTr) > nsigmaITSvar.nsigmaITSTr;
+        isAntiHe = isHelium && track.sign() < 0 && std::abs(nITSTr) > nsigmaITSvar.nsigmaITSTr;
       }
 
       isDeWoDCAxy = isDe && passDCAzCutDe;
@@ -3469,9 +3471,6 @@ struct LFNucleiBATask {
         continue;
       }
 
-      nITSTr = track.itsNSigmaTr();
-      nITSHe = track.itsNSigmaHe();
-
       if (outFlagOptions.makeDCAAfterCutPlots) {
 
         if (isHeWTPCpid) {