Update NonId function. Change proton PID cuts

[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoESDTrackCut.cxx

diff --git a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
index 4e94a17ebe479c944b8d8dbabe1fe2cd5c7f404a..dcab41e73ff90d4d7e6ab5a16104ce3475f17942 100644 (file)
--- a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
+++ b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoESDTrackCut.cxx
@@ -354,7 +354,8 @@ bool AliFemtoESDTrackCut::Pass(const AliFemtoTrack* track)
   
          }
          else if (fMostProbable == 4) { // proton nsigma-PID required contour adjusting
-           if (IsProtonNSigma(track->P().Mag(), track->NSigmaTPCP(), track->NSigmaTOFP()) && IsProtonTPCdEdx(track->P().Mag(), track->TPCsignal()))
+           if (IsProtonNSigma(track->P().Mag(), track->NSigmaTPCP(), track->NSigmaTOFP()) // && IsProtonTPCdEdx(track->P().Mag(), track->TPCsignal())
+            )
              imost = 4;
          }
 
@@ -956,20 +957,14 @@ bool AliFemtoESDTrackCut::IsPionNSigma(float mom, float nsigmaTPCPi, float nsigm
 
 bool AliFemtoESDTrackCut::IsProtonNSigma(float mom, float nsigmaTPCP, float nsigmaTOFP)
 {
-  if(mom<0.8)
-    {
-      if(nsigmaTOFP<-999.)
-       {
-         if(TMath::Abs(nsigmaTPCP)<3.0) return true;
-         else return false;
+    if (mom > 0.8) {
+        if (TMath::Hypot( nsigmaTOFP, nsigmaTPCP )/TMath::Sqrt(2) < 3.0)
+            return true;
        }
-      else if(TMath::Abs(nsigmaTOFP)<3.0 && TMath::Abs(nsigmaTPCP)<3.0) return true;
+    else {
+        if (TMath::Abs(nsigmaTPCP) < 3.0)
+            return true;
     }
-  else if(nsigmaTOFP<-999.)
-    {
-      return false;
-    } 
-  else if(TMath::Abs(nsigmaTOFP)<3.0 && TMath::Abs(nsigmaTPCP)<3.0) return true;
 
   return false;
 }