ISAJET/code/logmgy.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION LOGMGY(IERR)
   3 C
   4 C       Set and check limits for dijet masses.
   5 C
   6 C       Note we use the convention that not setting an upper limit
   7 C       gives a fixed variable, even though that currently is not
   8 C       implemented in N-jet phase space.
   9 C
  10 #if defined(CERNLIB_IMPNONE)
  11       IMPLICIT NONE
  12 #endif
  13 #include "isajet/itapes.inc"
  14 #include "isajet/jetlim.inc"
  15 #include "isajet/primar.inc"
  16 #include "isajet/jetpar.inc"
  17 #include "isajet/const.inc"
  18 #include "isajet/dylim.inc"
  19 #include "isajet/keys.inc"
  20 #include "isajet/q1q2.inc"
  21 #include "isajet/mglims.inc"
  22 C
  23       REAL UNDEF
  24       INTEGER IERR
  25       DATA UNDEF/-.9E9/
  26 C
  27       LOGMGY=.TRUE.
  28 C
  29 C          Attempt to fix YHMG
  30 C
  31       IF(YHMGMN.GT.UNDEF.AND.YHMGMX.LT.UNDEF) THEN
  32         LOGMGY=.FALSE.
  33         RETURN
  34       ENDIF
  35 C
  36 C          No limits
  37 C
  38       IF(EHMGMN.LT.0.OR.EHMGMX.LT.0) THEN
  39         LOGMGY=.FALSE.
  40         RETURN
  41       ENDIF
  42       IF(YHMGMN.LT.UNDEF) THEN
  43         YHMGMX=LOG(ECM/EHMGMN)
  44         YHMGMN=-YHMGMX
  45       ENDIF
  46 C
  47       RETURN
  48       END