Protection in calc. of ZDR1, ZDR2 as in Pythia 6.3
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Jun 2006 12:08:55 +0000 (12:08 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Jun 2006 12:08:55 +0000 (12:08 +0000)
HIJING/hipyset1_35/lushow_hijing.F

index 17eec7c..eaf9213 100644 (file)
@@ -537,8 +537,8 @@ C...Check if chosen multiplet m1,m2,z1,z2 is physical.
           ISL(3-ISLM)=0 
           ISLM=3-ISLM   
         ELSEIF(ISL(1).EQ.1.AND.ISL(2).EQ.1) THEN    
-          ZDR1=MAX(0.,V(N+1,3)/V(N+1,4)-1.) 
-          ZDR2=MAX(0.,V(N+2,3)/V(N+2,4)-1.) 
+          ZDR1=MAX(0D0,V(N+1,3)/MAX(1D-6,V(N+1,4))-1D0)
+          ZDR2=MAX(0D0,V(N+2,3)/MAX(1D-6,V(N+2,4))-1D0)
           IF(ZDR2.GT.RLU_HIJING(0)*(ZDR1+ZDR2)) ISL(1)=0   
           IF(ISL(1).EQ.1) ISL(2)=0  
           IF(ISL(1).EQ.0) ISLM=1