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);
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);
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]);