fix for ITS SA nuclie PID (F.Prino,B.Hess,F.Noferini,S.Trogolo,J.Wiechula)

diff --git a/STEER/STEERBase/AliPIDResponse.cxx b/STEER/STEERBase/AliPIDResponse.cxx
index cdf975769aae00383a33ce2fa58473cf60cd2f00..16a851d3caf505bfcc620f042cf3739e4a93c4bd 100644 (file)
--- a/STEER/STEERBase/AliPIDResponse.cxx
+++ b/STEER/STEERBase/AliPIDResponse.cxx
@@ -2344,10 +2344,11 @@ AliPIDResponse::EDetPidStatus AliPIDResponse::GetComputeITSProbability  (const A
     //      for electrons also for ITS+TPC tracks
     Double_t bethe=fITSResponse.Bethe(momITS,(AliPID::EParticleType)j,isSA || (j==(Int_t)AliPID::kElectron))*chargeFactor;
     Double_t sigma=fITSResponse.GetResolution(bethe,nPointsForPid,isSA || (j==(Int_t)AliPID::kElectron));
-    if (TMath::Abs(dedx-bethe) > fRange*sigma) {
+    Double_t nSigma=fITSResponse.GetNumberOfSigmas(track, (AliPID::EParticleType)j);
+    if (TMath::Abs(nSigma) > fRange) {
       p[j]=TMath::Exp(-0.5*fRange*fRange)/sigma;
     } else {
-      p[j]=TMath::Exp(-0.5*(dedx-bethe)*(dedx-bethe)/(sigma*sigma))/sigma;
+      p[j]=TMath::Exp(-0.5*nSigma*nSigma)/sigma;
       mismatch=kFALSE;
     }
   }