ISAJET/code/logpt.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION LOGPT(IERR)
   3 C
   4 C       SET AND CHECK LIMITS FOR JET PT
   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       DATA PTLOW/1.0/
  15       DATA UNDEF/-.9E9/
  16 C
  17       LOGPT=.TRUE.
  18       DO 20 I=1,NJET
  19       FIXPT(I)=.FALSE.
  20 C
  21       IF(PTMIN(I).LT.UNDEF.AND.PTMAX(I).LT.UNDEF) THEN
  22         PTMAX(I)=PMAX(I)
  23         PTMIN(I)=PTLOW
  24         IF(KEYS(3).AND.I.EQ.3.AND.QTMIN.GT.0.) PTMIN(I)=QTMIN
  25         IF(PMIN(I).LT.PTMIN(I)) PMIN(I)=PTMIN(I)
  26       ELSEIF(PTMAX(I).LT.UNDEF) THEN
  27         FIXPT(I)=.TRUE.
  28         PTMAX(I)=PTMIN(I)
  29       ELSEIF(PTMIN(I).LT.UNDEF) THEN
  30         PTMIN(I)=PTLOW
  31         IF(KEYS(3).AND.I.EQ.3.AND.QTMIN.GT.0.) PTMIN(I)=QTMIN
  32       ENDIF
  33 C
  34       IF(FIXPT(I)) PTMAX(I)=PTMIN(I)
  35       IF(FIXPT(I)) PT(I)=PTMIN(I)
  36       IF(PTMAX(I).GT.PMAX(I)) PTMAX(I)=PMAX(I)
  37       IF(PMIN(I).LT.PTMIN(I)) PMIN(I)=PTMIN(I)
  38 C
  39    20 CONTINUE
  40 C
  41       RETURN
  42       END