ISAJET/code/rejjet.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION REJJET()
   3 C-----------------------------------------------------------------------
   4 C-
   5 C-         This is called after EVOLVE for TWOJET and DRELLYAN events
   6 C-         to test the partons (jets). REJJET=.FALSE. keeps the event.
   7 C-
   8 C-----------------------------------------------------------------------
   9 #if defined(CERNLIB_IMPNONE)
  10       IMPLICIT NONE
  11 #endif
  12 #include "isajet/jetset.inc"
  13       INTEGER I,IFLAV
  14       REJJET=.FALSE.
  15 C***************************************
  16 C           Sample REJJET function which keeps the event if one of the
  17 C           outgoing partons is a heavy quark.
  18 C      REJJET=.TRUE.
  19 C      DO 1 I=1,NJSET
  20 C        IF(JDCAY(I).NE.0) GO TO 1
  21 C        IFLAV=IABS(JTYPE(I))
  22 C        IF(IFLAV.GE.4.AND.IFLAV.LT.9) THEN
  23 C          REJJET=.FALSE.
  24 C          RETURN
  25 C        ENDIF
  26 C  1   CONTINUE
  27 C***************************************
  28       RETURN
  29       END