1 #include "isajet/pilot.h"
2 LOGICAL FUNCTION LOGX(IERR)
4 C SET AND CHECK LIMITS FOR JET FEYNMAN X
6 #include "isajet/itapes.inc"
7 #include "isajet/jetlim.inc"
8 #include "isajet/primar.inc"
9 #include "isajet/jetpar.inc"
10 #include "isajet/const.inc"
11 #include "isajet/dylim.inc"
12 #include "isajet/keys.inc"
13 #include "isajet/q1q2.inc"
21 IF(FIXYJ(I).AND.(FIXP(I).OR.FIXPT(I)))FIXXJ(I)=.TRUE.
24 IF(XJMIN(I).LT.UNDEF.AND.XJMAX(I).LT.UNDEF) THEN
29 IF(XJMAX(I).LT.UNDEF) FIXXJ(I)=.TRUE.
30 IF(FIXXJ(I)) XJMAX(I)=XJMIN(I)
32 IF(.NOT.FIXXJ(I)) THEN
33 IF(THMIN(I).LT.HALFPI) X1=PMAX(I)*COS(THMIN(I))/HALFE
34 IF(THMIN(I).GE.HALFPI) X1=PMIN(I)*COS(THMIN(I))/HALFE
35 IF(THMAX(I).GT.HALFPI) X2=PMAX(I)*COS(THMAX(I))/HALFE
36 IF(THMAX(I).LT.HALFPI) X2=PMIN(I)*COS(THMAX(I))/HALFE
37 IF(X1.LT.XJMAX(I)) XJMAX(I)=X1
38 IF(X2.GT.XJMIN(I)) XJMIN(I)=X2
44 CTH(I)=XJ(I)*HALFE/P(I)
45 IF(ABS(CTH(I)).LE.1.0) THEN
46 STH(I)=SQRT(1.-CTH(I)**2)
47 TH(I)=ATAN2(STH(I),CTH(I))
48 YJ(I)=-ALOG(TAN(TH(I)/2.))
63 TH(I)=ATAN(PT(I)/XJ(I)/HALFE)
65 YJ(I)=-ALOG(TAN(TH(I)/2.))
78 P(I)=XJ(I)*HALFE/CTH(I)