]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | LOGICAL FUNCTION YGENJ(I) | |
3 | C | |
4 | C GENERATE Y FOR TWOJET | |
5 | C | |
6 | #include "isajet/itapes.inc" | |
7 | #include "isajet/jetlim.inc" | |
8 | #include "isajet/jetpar.inc" | |
9 | #include "isajet/primar.inc" | |
10 | #include "isajet/ptpar.inc" | |
11 | #include "isajet/totals.inc" | |
12 | ACOSH(X)=ALOG(X+SQRT(X**2-1.0)) | |
13 | YGENJ=.TRUE. | |
14 | YMAX=ACOSH(HALFE/PT(I)) | |
15 | YMIN=-YMAX | |
16 | IF(YMAX.LT.YJMIN(I).OR.YMIN.GT.YJMAX(I)) GOTO 10 | |
17 | YJ(I)=YJMIN(I)+(YJMAX(I)-YJMIN(I))*RANF() | |
18 | IF(YJ(I).LT.YMIN.OR.YJ(I).GT.YMAX) GOTO 10 | |
19 | TH(I)=2.*ATAN(EXP(-YJ(I))) | |
20 | CTH(I)=COS(TH(I)) | |
21 | STH(I)=SIN(TH(I)) | |
22 | WT=WT*(YJMAX(I)-YJMIN(I)) | |
23 | RETURN | |
24 | 10 YGENJ=.FALSE. | |
25 | RETURN | |
26 | END |