ISAJET/code/logphw.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION LOGPHW(IERR,DELPH)
   3 C
   4 C       SET AND CHECK LIMITS FOR W(Z0) PHI
   5 C
   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"
  14       DIMENSION DELPH(3)
  15       DATA UNDEF/-.9E9/
  16 C
  17       LOGPHW=.TRUE.
  18       FIXPHW=.FALSE.
  19 C
  20       IF(FIXPHI(3)) THEN
  21           FIXPHW=.TRUE.
  22           PHIW=AMOD(PHI(3)+PI,2.*PI)
  23       ELSEIF(PHWMIN.LT.UNDEF.AND.PHWMAX.LT.UNDEF) THEN
  24           PHWMIN=0.
  25           PHWMAX=2.*PI
  26       ELSEIF(PHWMAX.LT.UNDEF) THEN
  27           FIXPHW=.TRUE.
  28           PHW=PHWMIN
  29           FIXPHI(3)=.TRUE.
  30           PHWMAX=PHWMIN
  31           PHI(3)=PHIW+PI
  32           PHIMIN(3)=PHIW
  33           PHIMAX(3)=PHIW
  34       ENDIF
  35 C
  36       DELPHW=PHWMAX-PHWMIN
  37 C
  38       IF(DELPHW.LT.0.OR.DELPHW.GT.2.*PI) THEN
  39         CALL LOGERR(110,1,IERR)
  40         LOGPHW=.FALSE.
  41       ENDIF
  42 C
  43       IF(DELPHW.LE.DELPH(3)) THEN
  44         PHIMIN(3)=PHWMIN+PI
  45         PHIMAX(3)=PHIMIN(3)+DELPHW
  46       ELSE
  47         PHWMIN=PHIMIN(3)+PI
  48         PHWMAX=PHWMIN+DELPH(3)
  49       ENDIF
  50 C
  51       RETURN
  52       END