From: shahoian Date: Mon, 2 Feb 2015 09:08:42 +0000 (+0100) Subject: fix for ITS SA nuclie PID (F.Prino,B.Hess,F.Noferini,S.Trogolo,J.Wiechula) X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=04f7af3e2053dbd51e781bd764355e4f95e18ca1 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 cdf975769aa..16a851d3caf 100644 --- 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; } }