]> git.uio.no Git - u/mrichter/AliRoot.git/blob - ISAJET/code/logqt.F
Bug in V0A fixed (Guillermo)
[u/mrichter/AliRoot.git] / ISAJET / code / logqt.F
1 #include "isajet/pilot.h"
2       LOGICAL FUNCTION LOGQT(IERR)
3 C
4 C      SET AND CHECK W(Z0) PT RANGE
5 C
6 #include "isajet/itapes.inc"
7 #include "isajet/jetlim.inc"
8 #include "isajet/primar.inc"
9 #include "isajet/jetpar.inc"
10 #include "isajet/const.inc"
11 #include "isajet/dylim.inc"
12 #include "isajet/keys.inc"
13 #include "isajet/q1q2.inc"
14       DATA UNDEF/-.9E9/
15 C
16       LOGQT=.TRUE.
17       FIXQT=.FALSE.
18       IF(QTMIN.LT.UNDEF.AND.QTMAX.LT.UNDEF) THEN
19         QTMAX=PTMAX(3)
20         QTMIN=PTMIN(3)
21 C          SET DEFAULT QTW LIMITS IF NONE WERE SET
22         IF(QTMAX.GT.0.99*HALFE) THEN
23           NJET=2
24           QTMIN=0.
25           QTMAX=0.
26           QTW=0.
27           STDDY=.TRUE.
28           FIXQT=.TRUE.
29           PTMIN(3)=0.
30           PTMAX(3)=0.
31           FIXPT(3)=.TRUE.
32           CALL LOGERR(0,1,IERR)
33         ENDIF
34       ELSEIF(FIXPT(3)) THEN
35         QTW=PT(3)
36         QTMIN=PTMIN(3)
37         QTMAX=QTMIN
38         FIXQT=.TRUE.
39       ELSEIF(QTMAX.LT.UNDEF) THEN
40         FIXQT=.TRUE.
41         QTW=QTMIN
42         QTMAX=QTMIN
43         FIXPT(3)=.TRUE.
44         PT(3)=QTW
45         PTMIN(3)=QTMIN
46         PTMAX(3)=QTMAX
47       ELSE
48         IF(QTMAX.LT.PTMAX(3)) PTMAX(3)=QTMAX
49         IF(QTMIN.GT.PTMIN(3)) PTMIN(3)=QTMIN
50         IF(QTMAX.GT.PTMAX(3)) QTMAX=PTMAX(3)
51         IF(QTMIN.LT.PTMIN(3)) QTMIN=PTMIN(3)
52       ENDIF
53 C
54       RETURN
55       END