X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=HIJING%2Fhijing1_36%2Fg.F;h=67c05ca6766553b2d3bad88ab506afeb8ef8a3b1;hb=b5738316e7fb48c59aafe665b1c3e657790fec19;hp=fd4898578fd58937a93b4a456b33d2394350b767;hpb=e74335a4fa134f7a8870dae77bf84fbc68d77184;p=u%2Fmrichter%2FAliRoot.git

diff --git a/HIJING/hijing1_36/g.F b/HIJING/hijing1_36/g.F
index fd4898578fd..67c05ca6766 100644
--- a/HIJING/hijing1_36/g.F
+++ b/HIJING/hijing1_36/g.F
@@ -5,13 +5,16 @@ C
 C
 	FUNCTION G(Y1,Y2,PT2)
 	IMPLICIT REAL*8 (A-H,O-Z)
+#define BLANKET_SAVE
 #include "hiparnt.inc"
 	DIMENSION F(2,7)
 	SAVE
 	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