Additional protection
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 13 Aug 2007 17:59:47 +0000 (17:59 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 13 Aug 2007 17:59:47 +0000 (17:59 +0000)
HIJING/hijing1_36/fjet.F

index 3cc16cf..34bb1a0 100644 (file)
@@ -10,6 +10,7 @@ C
        SAVE
        PT2=(HINT1(1)**2/4.0-HIPR1(8)**2)*X(1)+HIPR1(8)**2
        XT=2.0*DSQRT(PT2)/HINT1(1)
+       IF (XT.GT.1.0) XT=1.0 
        YMX1=DLOG(1.0/XT+DSQRT(1.0/XT**2-1.0))
        Y1=2.0*YMX1*X(2)-YMX1
        YMX2=DLOG(2.0/XT-DEXP(Y1))