Removing useless const to avoid warnings on alphacxx6

[u/mrichter/AliRoot.git] / MEVSIM / multiplicity_gen.F

diff --git a/MEVSIM/multiplicity_gen.F b/MEVSIM/multiplicity_gen.F
index 4c7e10bcb6d7d633d69b31731b8dc23c15ec7a0a..8d98273b608f460b9423b94d5bbfcadd4af3e689 100644 (file)
--- a/MEVSIM/multiplicity_gen.F
+++ b/MEVSIM/multiplicity_gen.F
@@ -1,3 +1,10 @@
+#ifdef __APPLE__
+#define STOP CALL EXIT !
+#define stop CALL EXIT !
+#endif
+      subroutine ah
+      STOP
+      end
 C      Program Mult_Gen
        SUBROUTINE multgen
       implicit none
@@ -423,7 +430,7 @@ CCC   Variables associated with trigger fluctuations:
       real*4 MultFac_xfunc_save(nmax_integ)
 CCC  Open I/O Files:
 
-      open(unit=4,type='old',access='sequential',name='mult_gen.in')
+      open(unit=4,status='old',access='sequential',file='mult_gen.in')
 C-->  Commented for AliRoot direct COMMON block access
 C      open(unit=7,type='new',access='sequential',name='mult_gen.out')
 C      open(unit=8,type='new',access='sequential',name='mult_gen.log')
@@ -1521,6 +1528,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 +1701,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')