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