changes in the PID of protons + new config for the proton femtoscopy train
[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / AliFemtoUser / AliFemtoESDTrackCut.cxx
index f2b64445c4f0fe65e934c0d3c803468cfc6f112b..f0a9c78c86453714f0a080be94c81b2e866af35c 100644 (file)
@@ -377,12 +377,16 @@ bool AliFemtoESDTrackCut::Pass(const AliFemtoTrack* track)
            }
 
          }
-      else if (fMostProbable == 4) { // proton nsigma-PID required contour adjusting (in LHC10h)
-        if ( IsProtonNSigma(track->P().Mag(), track->NSigmaTPCP(), track->NSigmaTOFP()) && (TMath::Abs(track->NSigmaTPCP()) < TMath::Abs(track->NSigmaTPCPi())) && (TMath::Abs(track->NSigmaTPCP()) < TMath::Abs(track->NSigmaTPCK())) && (TMath::Abs(track->NSigmaTOFP()) < TMath::Abs(track->NSigmaTOFPi())) && (TMath::Abs(track->NSigmaTOFP()) < TMath::Abs(track->NSigmaTOFK()))
-             // && IsProtonTPCdEdx(track->P().Mag(), track->TPCsignal())
-            )
+    else if (fMostProbable == 4) { // proton nsigma-PID required contour adjusting (in LHC10h)
+      if ( IsProtonNSigma(track->P().Mag(), track->NSigmaTPCP(), track->NSigmaTOFP()) // && (TMath::Abs(track->NSigmaTPCP()) < TMath::Abs(track->NSigmaTPCPi())) && (TMath::Abs(track->NSigmaTPCP()) < TMath::Abs(track->NSigmaTPCK())) && (TMath::Abs(track->NSigmaTOFP()) < TMath::Abs(track->NSigmaTOFPi())) && (TMath::Abs(track->NSigmaTOFP()) < TMath::Abs(track->NSigmaTOFK()))
+           // && IsProtonTPCdEdx(track->P().Mag(), track->TPCsignal())
+        )
              imost = 4;
          }
+    else if (fMostProbable == 5) { // no-protons
+      if ( !IsProtonNSigma(track->P().Mag(), track->NSigmaTPCP(), track->NSigmaTOFP()) )
+             imost = 5;
+         }
 
        }