ISAJET/code/ygenj.F

   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