Additional protection
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Feb 2004 16:33:56 +0000 (16:33 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Feb 2004 16:33:56 +0000 (16:33 +0000)
HIJING/hijing1_36/g.F

index fd48985..ce26aca 100644 (file)
@@ -11,7 +11,9 @@ C
        XT=2.0*DSQRT(PT2)/HINT1(1)
        X1=0.50*XT*(DEXP(Y1)+DEXP(Y2))
        X2=0.50*XT*(DEXP(-Y1)+DEXP(-Y2))
-       Z=DSQRT(1.D0-XT**2/X1/X2)
+       DUMMY = 1.D0-XT**2/X1/X2
+       if (DUMMY.LT.0.D0) DUMMY = 0.D0
+       Z=DSQRT(DUMMY)
        SS=X1*X2*HINT1(1)**2
        T=-(1.0-Z)/2.0
        U=-(1.0+Z)/2.0