| 141 |
JFAIL=0 !error flag |
JFAIL=0 !error flag |
| 142 |
CHI2=0 |
CHI2=0 |
| 143 |
|
|
| 144 |
if(TRKDEBUG) print*,'guess: ',al |
if(TRKDEBUG) print*,'mini : guess ',al |
| 145 |
if(TRKDEBUG) print*,'mini2: step ',istep,chi2,AL(5) |
if(TRKDEBUG) print*,'mini : step ',istep,chi2,AL(5) |
| 146 |
|
|
| 147 |
* |
* |
| 148 |
* ----------------------- |
* ----------------------- |
| 252 |
ENDDO |
ENDDO |
| 253 |
ENDIF |
ENDIF |
| 254 |
|
|
| 255 |
if(TRKDEBUG) print*,'mini2: step ',istep,chi2,AL(5) |
if(TRKDEBUG) print*,'mini : step ',istep,chi2,AL(5) |
| 256 |
|
|
| 257 |
c$$$ PRINT*,'DAL ',(DAL(K),K=1,5) |
c$$$ PRINT*,'DAL ',(DAL(K),K=1,5) |
| 258 |
c$$$ PRINT*,'CHI2DOLD ',(CHI2DOLD(K),K=1,5) |
c$$$ PRINT*,'CHI2DOLD ',(CHI2DOLD(K),K=1,5) |
| 563 |
|
|
| 564 |
11 CONTINUE |
11 CONTINUE |
| 565 |
|
|
| 566 |
if(TRKDEBUG) print*,'mini2: -ok- ',istep,chi2,AL(5) |
if(TRKDEBUG) print*,'mini : -ok- ',istep,chi2,AL(5) |
| 567 |
|
|
| 568 |
NSTEP=ISTEP ! ***PP*** |
NSTEP=ISTEP ! ***PP*** |
| 569 |
|
|
| 1142 |
|
|
| 1143 |
REAL*4 XP(NPLANES),ZP(NPLANES),AP(NPLANES),RP(NPLANES) |
REAL*4 XP(NPLANES),ZP(NPLANES),AP(NPLANES),RP(NPLANES) |
| 1144 |
REAL*4 CHI,XC,ZC,RADIUS |
REAL*4 CHI,XC,ZC,RADIUS |
| 1145 |
|
|
| 1146 |
|
c$$$ DO I=1,nplanes |
| 1147 |
|
c$$$ print *,i,' - ',XGOOD(I),YGOOD(I) |
| 1148 |
|
c$$$ print *,i,' - ',xm(i),ym(i),zm(i) |
| 1149 |
|
c$$$ print *,i,' A ',xm_a(i),ym_a(i),zm_a(i) |
| 1150 |
|
c$$$ print *,i,' B ',xm_b(i),ym_b(i),zm_b(i) |
| 1151 |
|
c$$$ enddo |
| 1152 |
|
|
| 1153 |
|
|
| 1154 |
* ---------------------------------------- |
* ---------------------------------------- |
| 1155 |
* Y view |
* Y view |
| 1156 |
* ---------------------------------------- |
* ---------------------------------------- |