Be sure to load mapping when needed
[u/mrichter/AliRoot.git] / TRD / AliTRDpidESD.cxx
index 1087b0e..705103a 100644 (file)
@@ -159,7 +159,7 @@ Int_t AliTRDpidESD::MakePID(AliESDEvent * const event)
 
             
     // Skip tracks which have no TRD signal at all
-    if (t->GetTRDsignal() == 0.) continue;
+    if (t->GetTRDsignal()<1.e-5) continue;
   
     // Loop over detector layers
     mom          = 0.;
@@ -281,13 +281,13 @@ Bool_t AliTRDpidESD::RecalculateTrackSegmentKine(AliESDtrack * const esd
     mom    = op->GetP();
     s      = op->GetSnp();
     t      = op->GetTgl();
-          if (s < 1.) length /= TMath::Sqrt((1. - s*s) / (1. + t*t));
+    if (s < 1.) length /= TMath::Sqrt((1.-s)*(1.+s) / (1. + t*t));
     return kFALSE;
   }
   mom        = param->GetP();
   s = param->GetSnp();
   t = param->GetTgl();
-  if (s < 1.) length    /= TMath::Sqrt((1. - s*s) / (1. + t*t));
+  if (s < 1.) length    /= TMath::Sqrt((1.-s)*(1.+s) / (1. + t*t));
 
   // check if track is crossing tracking sector by propagating to chamber exit- maybe is too much :)
   Double_t alpha = param->GetAlpha();