ISAJET/code/rejfrg.F

   1 #include "isajet/pilot.h"
   2       LOGICAL FUNCTION REJFRG()
   3 C-----------------------------------------------------------------------
   4 C-
   5 C-         This is called after FRGMNT for TWOJET and DRELLYAN events
   6 C-         to test the fragmentation. REJFRG=.FALSE. keeps the event.
   7 C-
   8 C-----------------------------------------------------------------------
   9 #if defined(CERNLIB_IMPNONE)
  10       IMPLICIT NONE
  11 #endif
  12 #include "isajet/partcl.inc"
  13       INTEGER I,ID
  14       REAL PTL,CUTLEP,CUTNU
  15       DATA CUTLEP/50./,CUTNU/100./
  16       REJFRG=.FALSE.
  17 C***************************************
  18 C           Sample REJFRG function which keeps the event if it contains
  19 C           any lepton satisfying
  20 C             PT > CUTLEP (charged lepton)
  21 C             PT > CUTNU  (neutrino)
  22 C           Appropriate values of the cuts must be set by the user.
  23 C      REJFRG=.TRUE.
  24 C      DO 1 I=1,NPTCL
  25 C        IF(IDCAY(I).NE.0) GO TO 1
  26 C        ID=IABS(IDENT(I))
  27 C        IF(ID.LE.10.OR.ID.GE.20) GO TO 1
  28 C        PTL=SQRT(PPTCL(1,I)**2+PPTCL(2,I)**2)
  29 C        IF((ID.EQ.11.OR.ID.EQ.13.OR.ID.EQ.15).AND.PTL.GT.CUTNU) THEN
  30 C          REJFRG=.FALSE.
  31 C          RETURN
  32 C        ELSEIF((ID.EQ.12..OR.ID.EQ.14).AND.PTL.GT.CUTLEP) THEN
  33 C          REJFRG=.FALSE.
  34 C          RETURN
  35 C        ENDIF
  36 C   1  CONTINUE
  37 C***************************************
  38       RETURN
  39       END