| 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. |