ISAJET/code/logp.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION LOGP(IERR)
   3 C
   4 C       SET AND CHECK LIMITS FOR JET MOMENTA
   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 PLOW/1.0/
  15       DATA UNDEF/-.9E9/
  16 C
  17       LOGP=.TRUE.
  18       DO 10 I=1,NJET
  19       FIXP(I)=.FALSE.
  20       IF(PMIN(I).LT.UNDEF.AND.PMAX(I).LT.UNDEF)  PMAX(I)=HALFE
  21       IF(PMAX(I).GT.HALFE) PMAX(I)=HALFE
  22       IF(PMAX(I).LT.UNDEF) FIXP(I)=.TRUE.
  23       IF(PMIN(I).LT.UNDEF) PMIN(I)=PLOW
  24       IF(FIXP(I)) THEN
  25         PMAX(I)=PMIN(I)
  26         P(I)=PMIN(I)
  27       ENDIF
  28       IF(KEYS(3).AND.I.EQ.3.AND.QTMIN.GT.0) PMIN(I)=QTMIN
  29    10 CONTINUE
  30 C
  31       RETURN
  32       END