[u/mrichter/AliRoot.git] / ISAJET / code / rejjet.F

#include "isajet/pilot.h"
      LOGICAL FUNCTION REJJET() 
C-----------------------------------------------------------------------
C-
C-         This is called after EVOLVE for TWOJET and DRELLYAN events 
C-         to test the partons (jets). REJJET=.FALSE. keeps the event.
C-
C-----------------------------------------------------------------------
#if defined(CERNLIB_IMPNONE)
      IMPLICIT NONE
#endif
#include "isajet/jetset.inc"
      INTEGER I,IFLAV  
      REJJET=.FALSE.
C***************************************
C           Sample REJJET function which keeps the event if one of the 
C           outgoing partons is a heavy quark. 
C      REJJET=.TRUE. 
C      DO 1 I=1,NJSET    
C        IF(JDCAY(I).NE.0) GO TO 1   
C        IFLAV=IABS(JTYPE(I))    
C        IF(IFLAV.GE.4.AND.IFLAV.LT.9) THEN
C          REJJET=.FALSE.  
C          RETURN  
C        ENDIF
C  1   CONTINUE  
C***************************************
      RETURN    
      END
Commit	Line	Data
0795afa3	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