]>
Commit | Line | Data |
---|---|---|
0795afa3 | 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 |