ran declared external
[u/mrichter/AliRoot.git] / MEVSIM / multiplicity_gen.F
index 4c7e10bcb6d7d633d69b31731b8dc23c15ec7a0a..23aa4bdd404908adfaafe0d6296237213b1e80c8 100644 (file)
@@ -1521,6 +1521,8 @@ C-OUT      Close(Unit=10)
       integer reac_plane_cntrl,iphi
       integer n_integ_pts
 
+      external ran
+
       do i = 1,factorial_max
       do j = 1,4
          pout(pid,j,i) = 0.0
@@ -1692,9 +1694,11 @@ CCC   This Function computes the pseudorapidty for a given mass, pt, y:
          pz = (y/abs(y))*pzmag
       end if
       pmag = sqrt(pt*pt + pz*pz)
-      if(pt.ne.0.0) then
+      if( (pt.ne.0.0) .and. (pmag .ne.pz) .and. (pmag.ne.(-pz)) ) then
+      
          pseudorapidity = 0.5*log((pmag+pz)/(pmag-pz))
-      else if (pt.eq.0.0) then
+      else 
+CCC      if (pt.eq.0.0) then
          pseudorapidity = 86.0
 C-->         write(8,10)
 10       Format(10x,'Function pseudorapidity called with pt=0')