Casuality cuts added (Yu.Belikov)

[u/mrichter/AliRoot.git] / TOF / AliTOFpidESD.cxx

diff --git a/TOF/AliTOFpidESD.cxx b/TOF/AliTOFpidESD.cxx
index 3987e44d95c4dd476a40525295a992c4a1a7e841..38e7b8a6c3f3d3eea8369082b7e73d244396b311 100644 (file)
--- a/TOF/AliTOFpidESD.cxx
+++ b/TOF/AliTOFpidESD.cxx
@@ -170,6 +170,8 @@ Int_t AliTOFpidESD::MakePID(AliESD *event)
     if ((t->GetStatus()&AliESDtrack::kTRDout)==0) continue;
     if ((t->GetStatus()&AliESDtrack::kTRDStop)!=0) continue;
 
+    Double_t time[10]; t->GetIntegratedTimes(time);
+
     Double_t x,par[5]; t->GetExternalParameters(x,par);
     Double_t cov[15]; t->GetExternalCovariance(cov);
 
@@ -189,8 +191,12 @@ Double_t dz=5*TMath::Sqrt(cov[2]) + 0.5*fDz + 2.5*TMath::Abs(par[3]);
       if (c->GetZ() > z+dz) break;
       if (c->IsUsed()) continue;
 
+      Double_t tof=50*c->GetTDC()+32;
+      if (t->GetIntegratedLength()/tof > 0.031) continue;
+      if (tof>35000) continue;
+
       Double_t dph=TMath::Abs(c->GetPhi()-phi);
-      if (dph>TMath::Pi()) dph-=2*TMath::Pi();
+      if (dph>TMath::Pi()) dph=2*TMath::Pi()-dph; //Thanks to B.Zagreev
       if (dph>dphi) continue;
 
       Double_t d2=dph*dph*fR*fR + (c->GetZ()-z)*(c->GetZ()-z);
@@ -216,8 +222,6 @@ Double_t dz=5*TMath::Sqrt(cov[2]) + 0.5*fDz + 2.5*TMath::Abs(par[3]);
 
     if ((t->GetStatus()&AliESDtrack::kTIME)==0) continue;
 
-    Double_t time[10]; t->GetIntegratedTimes(time);
-
     //track length correction
     Double_t rc=TMath::Sqrt(c->GetR()*c->GetR() + c->GetZ()*c->GetZ());
     Double_t rt=TMath::Sqrt(x*x + par[0]*par[0] + par[1]*par[1]);