Additional protection
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Dec 2003 16:02:00 +0000 (16:02 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 19 Dec 2003 16:02:00 +0000 (16:02 +0000)
HIJING/hipyset1_35/luindf_hijing.F

index 3943279..c577e15 100644 (file)
@@ -117,8 +117,9 @@ C...sharing energy according to Altarelli-Parisi splitting function.
         CALL LUPTDI_HIJING(0,PXO(1),PYO(1))    
         PXO(2)=-PXO(1)  
         PYO(2)=-PYO(1)  
-        WO(1)=WF*RLU_HIJING(0)**(1./3.)    
-        WO(2)=WF-WO(1)  
+ 111    WO(1)=WF*RLU_HIJING(0)**(1./3.)    
+        WO(2)=WF-WO(1)
+        IF (ABS(WO(2)).LE.0..OR.ABS(WO(1)).LE.0.) GOTO 111
       ENDIF 
     
 C...Initial values for rank, flavour, pT and W+.