]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | LOGICAL FUNCTION LOGMGY(IERR) | |
3 | C | |
4 | C Set and check limits for dijet masses. | |
5 | C | |
6 | C Note we use the convention that not setting an upper limit | |
7 | C gives a fixed variable, even though that currently is not | |
8 | C implemented in N-jet phase space. | |
9 | C | |
10 | #if defined(CERNLIB_IMPNONE) | |
11 | IMPLICIT NONE | |
12 | #endif | |
13 | #include "isajet/itapes.inc" | |
14 | #include "isajet/jetlim.inc" | |
15 | #include "isajet/primar.inc" | |
16 | #include "isajet/jetpar.inc" | |
17 | #include "isajet/const.inc" | |
18 | #include "isajet/dylim.inc" | |
19 | #include "isajet/keys.inc" | |
20 | #include "isajet/q1q2.inc" | |
21 | #include "isajet/mglims.inc" | |
22 | C | |
23 | REAL UNDEF | |
24 | INTEGER IERR | |
25 | DATA UNDEF/-.9E9/ | |
26 | C | |
27 | LOGMGY=.TRUE. | |
28 | C | |
29 | C Attempt to fix YHMG | |
30 | C | |
31 | IF(YHMGMN.GT.UNDEF.AND.YHMGMX.LT.UNDEF) THEN | |
32 | LOGMGY=.FALSE. | |
33 | RETURN | |
34 | ENDIF | |
35 | C | |
36 | C No limits | |
37 | C | |
38 | IF(EHMGMN.LT.0.OR.EHMGMX.LT.0) THEN | |
39 | LOGMGY=.FALSE. | |
40 | RETURN | |
41 | ENDIF | |
42 | IF(YHMGMN.LT.UNDEF) THEN | |
43 | YHMGMX=LOG(ECM/EHMGMN) | |
44 | YHMGMN=-YHMGMX | |
45 | ENDIF | |
46 | C | |
47 | RETURN | |
48 | END |