Protection against division by 0 added.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Jan 2005 10:58:03 +0000 (10:58 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Jan 2005 10:58:03 +0000 (10:58 +0000)
HIJING/hipyset1_35/pyremn_hijing.F

index d9bff03..278ee55 100644 (file)
@@ -261,6 +261,7 @@ C...Relative distribution of energy for particle into two jets.
           ENDIF 
 C...Relative distribution of energy for particle into jet plus particle.    
         ELSE    
+  111     CONTINUE 
           IF(MSTP(92).LE.1) THEN    
             IF(IMB.EQ.1) CHI(JT)=RLU_HIJING(0) 
             IF(IMB.EQ.2) CHI(JT)=1.-SQRT(RLU_HIJING(0))    
@@ -268,6 +269,7 @@ C...Relative distribution of energy for particle into jet plus particle.
             CHI(JT)=1.-RLU_HIJING(0)**(1./(1.+PARP(93+2*IMB))) 
           ENDIF 
           IF(MOD(KFLCH(JT)/1000,10).NE.0) CHI(JT)=1.-CHI(JT)    
+          IF (CHI(JT) .EQ. 1. .OR. CHI(JT) .EQ. 0.) GOTO 111
         ENDIF   
         PMS(JT)=PMS(JT+4)/CHI(JT)+PMS(JT+2)/(1.-CHI(JT))    
         KFLS=KCHG(LUCOMP_HIJING(KFLCH(JT)),2)*ISIGN(1,KFLCH(JT))