--- DarthVader/TrackerLevel2/src/F77/mini.f 2008/12/05 08:26:47 1.23 +++ DarthVader/TrackerLevel2/src/F77/mini.f 2008/12/23 11:28:36 1.24 @@ -141,7 +141,7 @@ CHI2=0 if(TRKDEBUG) print*,'guess: ',al - if(TRKDEBUG) print*,'mini2: step ',istep,chi2,1./AL(5) + if(TRKDEBUG) print*,'mini2: step ',istep,chi2,AL(5) * * ----------------------- @@ -251,7 +251,7 @@ ENDDO ENDIF - if(TRKDEBUG) print*,'mini2: step ',istep,chi2,1./AL(5) + if(TRKDEBUG) print*,'mini2: step ',istep,chi2,AL(5) c$$$ PRINT*,'DAL ',(DAL(K),K=1,5) c$$$ PRINT*,'CHI2DOLD ',(CHI2DOLD(K),K=1,5) @@ -418,6 +418,10 @@ * --------------------------------------------- *------------------------------------------------------------* c$$$ ALTOL(5) = DSQRT(DELETA1**2+DELETA2**2*AL(5)**2)/FACT + IF(FACT.EQ.0)THEN + IFAIL=1 + RETURN + ENDIF ALTOL(5) = DSQRT((DELETA1*AVRESX)**2+DELETA2**2*AL(5)**2)/FACT ALTOL(1) = ALTOL(5)/DELETA1 ALTOL(2) = ALTOL(1) @@ -554,12 +558,11 @@ * ------------------------------------ * Reduced chi^2 CHI2 = CHI2/dble(ndof) - c print*,'mini2: chi2 ',chi2 11 CONTINUE - if(TRKDEBUG) print*,'mini2: -ok- ',istep,chi2,1./AL(5) + if(TRKDEBUG) print*,'mini2: -ok- ',istep,chi2,AL(5) NSTEP=ISTEP ! ***PP***