impoved num precision

[u/mrichter/AliRoot.git] / TUHKMgen / UHKM / InitialStateHydjet.cxx

diff --git a/TUHKMgen/UHKM/InitialStateHydjet.cxx b/TUHKMgen/UHKM/InitialStateHydjet.cxx
index c2c50d31bb0500cb0fe317844bad6614b085ad20..64e9106b6a943416c2d49801276a5cdb7fc59fa0 100644 (file)
--- a/TUHKMgen/UHKM/InitialStateHydjet.cxx
+++ b/TUHKMgen/UHKM/InitialStateHydjet.cxx
@@ -239,9 +239,9 @@ void InitialStateHydjet::Initialize(List_t &source, ParticleAllocator & allocato
                
            Double_t php0 = TMath::TwoPi() * gRandom->Rndm();
            Double_t ctp0 = 2. * gRandom->Rndm() - 1.;
-           Double_t stp0 = TMath::Sqrt(1. - ctp0 * ctp0); 
+           Double_t stp0 = TMath::Sqrt((1.-ctp0)*(1.+ctp0)); 
            e = mass + (eMax - mass) * arrayFunctDistE(); 
-           Double_t pp0 = TMath::Sqrt(e * e - mass * mass);
+           Double_t pp0 = TMath::Sqrt((e-mass)*(e+mass));
            px0 = pp0 * stp0 * TMath::Sin(php0); 
            py0 = pp0 * stp0 * TMath::Cos(php0);
            pz0 = pp0 * ctp0;