More volume overlaps corrected
[u/mrichter/AliRoot.git] / ISAJET / code / rejfrg.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 LOGICAL FUNCTION REJFRG()
3C-----------------------------------------------------------------------
4C-
5C- This is called after FRGMNT for TWOJET and DRELLYAN events
6C- to test the fragmentation. REJFRG=.FALSE. keeps the event.
7C-
8C-----------------------------------------------------------------------
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.
17C***************************************
18C Sample REJFRG function which keeps the event if it contains
19C any lepton satisfying
20C PT > CUTLEP (charged lepton)
21C PT > CUTNU (neutrino)
22C Appropriate values of the cuts must be set by the user.
23C REJFRG=.TRUE.
24C DO 1 I=1,NPTCL
25C IF(IDCAY(I).NE.0) GO TO 1
26C ID=IABS(IDENT(I))
27C IF(ID.LE.10.OR.ID.GE.20) GO TO 1
28C PTL=SQRT(PPTCL(1,I)**2+PPTCL(2,I)**2)
29C IF((ID.EQ.11.OR.ID.EQ.13.OR.ID.EQ.15).AND.PTL.GT.CUTNU) THEN
30C REJFRG=.FALSE.
31C RETURN
32C ELSEIF((ID.EQ.12..OR.ID.EQ.14).AND.PTL.GT.CUTLEP) THEN
33C REJFRG=.FALSE.
34C RETURN
35C ENDIF
36C 1 CONTINUE
37C***************************************
38 RETURN
39 END