Some small corrections. (Piotr Skowronski)

diff --git a/MEVSIM/common_facfac.inc b/MEVSIM/common_facfac.inc
index 870b350ac43f91930eb63612bffc4961cda07b49..8ef92cf41653d11782c6e6af87531a7d5acc0a98 100644 (file)
--- a/MEVSIM/common_facfac.inc
+++ b/MEVSIM/common_facfac.inc
@@ -1,7 +1,7 @@
 CCC   Common FACFAC for Factorials
 
       integer factorial_max
-      parameter (factorial_max = 10000) ! max # multiplicity per event;
+      parameter (factorial_max = 80000) ! max # multiplicity per event;
 CCC                                     ! for any specific particle ID;
 CCC                                     ! also used for log(n!).
       Common/FACFAC/ FACLOG(factorial_max)

diff --git a/MEVSIM/multiplicity_gen.F b/MEVSIM/multiplicity_gen.F
index 4c7e10bcb6d7d633d69b31731b8dc23c15ec7a0a..17dcbaceb542fab351bc1e1a8a349ce6f3533ce3 100644 (file)
--- a/MEVSIM/multiplicity_gen.F
+++ b/MEVSIM/multiplicity_gen.F
@@ -1692,9 +1692,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')