]>
Commit | Line | Data |
---|---|---|
0795afa3 | 1 | #include "isajet/pilot.h" |
2 | LOGICAL FUNCTION LOGMIJ(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 AMLOW,UNDEF | |
24 | INTEGER I,J,IERR | |
25 | DATA AMLOW/1.0/ | |
26 | DATA UNDEF/-.9E9/ | |
27 | C | |
28 | LOGMIJ=.TRUE. | |
29 | C | |
30 | DO 100 I=1,MXLIM | |
31 | DO 101 J=I+1,MXLIM | |
32 | FIXMIJ(I,J)=.FALSE. | |
33 | FIXMIJ(J,I)=.FALSE. | |
34 | IF(AMIJMN(I,J).LT.UNDEF.AND.AMIJMX(I,J).LT.UNDEF) THEN | |
35 | AMIJMX(I,J)=ECM | |
36 | AMIJMX(J,I)=ECM | |
37 | ENDIF | |
38 | IF(AMIJMX(I,J).GT.ECM) THEN | |
39 | AMIJMX(I,J)=ECM | |
40 | AMIJMX(J,I)=ECM | |
41 | ENDIF | |
42 | IF(AMIJMX(I,J).LT.UNDEF) THEN | |
43 | AMIJMX(I,J)=AMIJMN(I,J) | |
44 | FIXMIJ(I,J)=.TRUE. | |
45 | AMIJMX(J,I)=AMIJMN(I,J) | |
46 | FIXMIJ(J,I)=.TRUE. | |
47 | ENDIF | |
48 | IF(AMIJMN(I,J).LT.UNDEF) THEN | |
49 | AMIJMN(I,J)=AMLOW | |
50 | AMIJMN(J,I)=AMLOW | |
51 | ENDIF | |
52 | 101 CONTINUE | |
53 | 100 CONTINUE | |
54 | C | |
55 | RETURN | |
56 | END |