First commit.
[u/mrichter/AliRoot.git] / ISAJET / code / logmgy.F
CommitLineData
0795afa3 1#include "isajet/pilot.h"
2 LOGICAL FUNCTION LOGMGY(IERR)
3C
4C Set and check limits for dijet masses.
5C
6C Note we use the convention that not setting an upper limit
7C gives a fixed variable, even though that currently is not
8C implemented in N-jet phase space.
9C
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"
22C
23 REAL UNDEF
24 INTEGER IERR
25 DATA UNDEF/-.9E9/
26C
27 LOGMGY=.TRUE.
28C
29C Attempt to fix YHMG
30C
31 IF(YHMGMN.GT.UNDEF.AND.YHMGMX.LT.UNDEF) THEN
32 LOGMGY=.FALSE.
33 RETURN
34 ENDIF
35C
36C No limits
37C
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
46C
47 RETURN
48 END