4 SUBROUTINE AR3JET(S,X1,X3,JL)
7 #include "lujets_hijing.inc"
11 IF(K(JL,2).NE.21 .AND. K(JL+1,2).NE.21) C=8./27.
14 IF(K(JL,2).NE.21) EXP1=2
15 IF(K(JL+1,2).NE.21) EXP3=2
17 YMA=ALOG(.5/SQRT(A)+SQRT(.25/A-1))
21 XT2M=(1.-2.*SQRT(SM1)+SM1-SM3)*(1.-2.*SQRT(SM3)-SM1+SM3)
24 1 IF(NTRY.EQ.5000) THEN
25 X1=.5*(2.*SQRT(SM1)+1.+SM1-SM3)
26 X3=.5*(2.*SQRT(SM3)+1.-SM1+SM3)
31 XT2=A*(XT2M/A)**(RLU_HIJING(0)**(1./D))
33 YMAX=ALOG(.5/SQRT(XT2)+SQRT(.25/XT2-1.))
34 Y=(2.*RLU_HIJING(0)-1.)*YMAX
35 X1=1.-SQRT(XT2)*EXP(Y)
36 X3=1.-SQRT(XT2)*EXP(-Y)
39 IF(K(JL,2).NE.21 .OR. K(JL+1,2).NE.21) THEN
40 IF((1.-X1)*(1.-X2)*(1.-X3)-X2*SM1*(1.-X1)-X2*SM3*(1.-X3).
41 & LE.0..OR.X1.LE.2.*SQRT(SM1)-SM1+SM3.OR.X3.LE.2.*SQRT(SM3)
48 FG=2.*YMAX*C*(X1**EXP1+X3**EXP3)/D
50 IF(FG.LT.RLU_HIJING(0)) GOTO 1