| e74335a4 |
1 | * $Id$ |
| 2 | |
| 3 | |
| 4 | SUBROUTINE AR3JET(S,X1,X3,JL) |
| 5 | C |
| bc676b8e |
6 | #define BLANKET_SAVE |
| e74335a4 |
7 | #include "hiparnt.inc" |
| 8 | #include "lujets_hijing.inc" |
| 9 | SAVE |
| 10 | C |
| 11 | C=1./3. |
| 12 | IF(K(JL,2).NE.21 .AND. K(JL+1,2).NE.21) C=8./27. |
| 13 | EXP1=3 |
| 14 | EXP3=3 |
| 15 | IF(K(JL,2).NE.21) EXP1=2 |
| 16 | IF(K(JL+1,2).NE.21) EXP3=2 |
| 17 | A=0.24**2/S |
| 18 | YMA=ALOG(.5/SQRT(A)+SQRT(.25/A-1)) |
| 19 | D=4.*C*YMA |
| 20 | SM1=P(JL,5)**2/S |
| 21 | SM3=P(JL+1,5)**2/S |
| 22 | XT2M=(1.-2.*SQRT(SM1)+SM1-SM3)*(1.-2.*SQRT(SM3)-SM1+SM3) |
| 23 | XT2M=MIN(.25,XT2M) |
| 24 | NTRY=0 |
| 25 | 1 IF(NTRY.EQ.5000) THEN |
| 26 | X1=.5*(2.*SQRT(SM1)+1.+SM1-SM3) |
| 27 | X3=.5*(2.*SQRT(SM3)+1.-SM1+SM3) |
| 28 | RETURN |
| 29 | ENDIF |
| 30 | NTRY=NTRY+1 |
| 31 | |
| 32 | XT2=A*(XT2M/A)**(RLU_HIJING(0)**(1./D)) |
| 33 | |
| 34 | YMAX=ALOG(.5/SQRT(XT2)+SQRT(.25/XT2-1.)) |
| 35 | Y=(2.*RLU_HIJING(0)-1.)*YMAX |
| 36 | X1=1.-SQRT(XT2)*EXP(Y) |
| 37 | X3=1.-SQRT(XT2)*EXP(-Y) |
| 38 | X2=2.-X1-X3 |
| 39 | NEG=0 |
| 40 | IF(K(JL,2).NE.21 .OR. K(JL+1,2).NE.21) THEN |
| 41 | IF((1.-X1)*(1.-X2)*(1.-X3)-X2*SM1*(1.-X1)-X2*SM3*(1.-X3). |
| 42 | & LE.0..OR.X1.LE.2.*SQRT(SM1)-SM1+SM3.OR.X3.LE.2.*SQRT(SM3) |
| 43 | & -SM3+SM1) NEG=1 |
| 44 | X1=X1+SM1-SM3 |
| 45 | X3=X3-SM1+SM3 |
| 46 | ENDIF |
| 47 | IF(NEG.EQ.1) GOTO 1 |
| 48 | |
| 49 | FG=2.*YMAX*C*(X1**EXP1+X3**EXP3)/D |
| 50 | XT2M=XT2 |
| 51 | IF(FG.LT.RLU_HIJING(0)) GOTO 1 |
| 52 | |
| 53 | RETURN |
| 54 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame