67 |
|
|
68 |
c LOGICAL TRKDEBUG,TRKVERBOSE |
c LOGICAL TRKDEBUG,TRKVERBOSE |
69 |
c COMMON/TRKD/TRKDEBUG,TRKVERBOSE |
c COMMON/TRKD/TRKDEBUG,TRKVERBOSE |
70 |
LOGICAL TRKDEBUG,TRKVERBOSE,STUDENT |
LOGICAL TRKDEBUG,TRKVERBOSE,STUDENT,FIRSTSTEPS,FIRSTSTUDENT |
71 |
COMMON/TRKD/TRKDEBUG,TRKVERBOSE |
COMMON/TRKD/TRKDEBUG,TRKVERBOSE |
72 |
|
|
73 |
DIMENSION AL0(5) |
DIMENSION AL0(5) |
76 |
* define kind of minimization (0x=chi2+gaussian or 1x=likelihood+student) |
* define kind of minimization (0x=chi2+gaussian or 1x=likelihood+student) |
77 |
* |
* |
78 |
STUDENT = .false. |
STUDENT = .false. |
79 |
|
FIRSTSTEPS = .true. |
80 |
|
FIRSTSTUDENT = .true. |
81 |
IF(MOD(INT(TRACKMODE/10),10).EQ.1) STUDENT = .true. |
IF(MOD(INT(TRACKMODE/10),10).EQ.1) STUDENT = .true. |
82 |
|
|
83 |
IF(IPRINT.EQ.1) THEN |
IF(IPRINT.EQ.1) THEN |
146 |
* **** Chi2+gaussian minimization |
* **** Chi2+gaussian minimization |
147 |
* ------------------------------- |
* ------------------------------- |
148 |
|
|
149 |
IF(.NOT.STUDENT) THEN |
IF(.NOT.STUDENT.OR.FIRSTSTEPS) THEN |
150 |
|
|
151 |
|
IF(ISTEP.GE.3) FIRSTSTEPS = .false. |
152 |
|
|
153 |
CALL CHISQ(IFLAG,JFAIL) !chi^2 and its derivatives |
CALL CHISQ(IFLAG,JFAIL) !chi^2 and its derivatives |
154 |
IF(JFAIL.NE.0) THEN |
IF(JFAIL.NE.0) THEN |
256 |
* **** Likelihood+Student minimization |
* **** Likelihood+Student minimization |
257 |
* ------------------------------- |
* ------------------------------- |
258 |
|
|
259 |
IF(STUDENT) THEN |
IF(STUDENT.AND.(.NOT.FIRSTSTEPS)) THEN |
260 |
|
|
261 |
|
IF(FIRSTSTUDENT) THEN |
262 |
|
FIRSTSTUDENT = .false. |
263 |
|
ISTEP = 1 |
264 |
|
ENDIF |
265 |
|
|
266 |
CALL CHISQSTT(1,JFAIL) |
CALL CHISQSTT(1,JFAIL) |
267 |
DO I=1,5 |
DO I=1,5 |
268 |
DAL(I)=0. |
DAL(I)=0. |