ISAJET/code/logmij.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION LOGMIJ(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 AMLOW,UNDEF
  24       INTEGER I,J,IERR
  25       DATA AMLOW/1.0/
  26       DATA UNDEF/-.9E9/
  27 C
  28       LOGMIJ=.TRUE.
  29 C
  30       DO 100 I=1,MXLIM
  31         DO 101 J=I+1,MXLIM
  32           FIXMIJ(I,J)=.FALSE.
  33           FIXMIJ(J,I)=.FALSE.
  34           IF(AMIJMN(I,J).LT.UNDEF.AND.AMIJMX(I,J).LT.UNDEF) THEN
  35             AMIJMX(I,J)=ECM
  36             AMIJMX(J,I)=ECM
  37           ENDIF
  38           IF(AMIJMX(I,J).GT.ECM) THEN
  39             AMIJMX(I,J)=ECM
  40             AMIJMX(J,I)=ECM
  41           ENDIF
  42           IF(AMIJMX(I,J).LT.UNDEF) THEN
  43             AMIJMX(I,J)=AMIJMN(I,J)
  44             FIXMIJ(I,J)=.TRUE.
  45             AMIJMX(J,I)=AMIJMN(I,J)
  46             FIXMIJ(J,I)=.TRUE.
  47           ENDIF
  48           IF(AMIJMN(I,J).LT.UNDEF) THEN
  49             AMIJMN(I,J)=AMLOW
  50             AMIJMN(J,I)=AMLOW
  51           ENDIF
  52 101     CONTINUE
  53 100   CONTINUE
  54 C
  55       RETURN
  56       END