| 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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame