Perform check on cos(crossing_angle_with_padrow) in GetCrossingAngle().
authorvestbo <vestbo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 2 Oct 2002 15:24:15 +0000 (15:24 +0000)
committervestbo <vestbo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 2 Oct 2002 15:24:15 +0000 (15:24 +0000)
Every once in a while this number is equal to one, which results in
return value being NaN.

HLT/src/AliL3Track.cxx

index 6a5a485..03e9bab 100644 (file)
@@ -189,8 +189,9 @@ Double_t AliL3Track::GetCrossingAngle(Int_t padrow)
   tangent[1] = (xyz[0] - GetCenterX())/GetRadius();
   
   Double_t perp_padrow[2] = {1,0}; //locally in slice
-
+  
   Double_t cos_beta = fabs(tangent[0]*perp_padrow[0] + tangent[1]*perp_padrow[1]);
+  if(cos_beta > 1) cos_beta=1;
   return acos(cos_beta);
 }