14-jun-2001 NvE Memberfunction GetRapidity() introduced for AliTrack and AliJet
[u/mrichter/AliRoot.git] / RALICE / Ali3Vector.cxx
index cdf6d31e141a7cecbb7e0c5ae82d7aa38c16dd75..a4cbdcb0160ffd8f4cd57b886b027c118f73c35d 100644 (file)
@@ -443,11 +443,13 @@ Double_t Ali3Vector::GetPseudoRapidity()
 // Provide the pseudo-rapidity w.r.t. the z-axis.
 // In other words : eta=-log(tan(theta/2))
 // The error on the scalar result (pseudo-rap.) is updated accordingly
+ Double_t pi=acos(-1.);
  Double_t v[3];
  GetVector(v,"sph");
  Double_t thetahalf=v[1]/2.;
- Double_t arg=tan(thetahalf);
- Double_t eta=0;
+ Double_t arg=0;
+ if (v[1]<pi) arg=tan(thetahalf);
+ Double_t eta=9999;
  if (arg>0) eta=-log(arg);
  Double_t e[3];
  GetErrors(e,"sph");