/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/mini.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/mini.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.10 by pam-fi, Thu Nov 30 17:01:52 2006 UTC revision 1.13 by pam-fi, Mon Feb 19 16:28:39 2007 UTC
# Line 121  c$$$      DELETA2 = 0.016/0.3/0.4/0.4451 Line 121  c$$$      DELETA2 = 0.016/0.3/0.4/0.4451
121  *      *    
122        ISTEP=0                   !num. steps to minimize chi^2        ISTEP=0                   !num. steps to minimize chi^2
123        JFAIL=0                   !error flag        JFAIL=0                   !error flag
124          CHI2=0
125    
126        if(TRKDEBUG) print*,'guess: ',al        if(TRKDEBUG) print*,'guess: ',al
127        if(TRKDEBUG) print*,'mini2: step ',istep,chi2,1./AL(5)        if(TRKDEBUG) print*,'mini2: step ',istep,chi2,1./AL(5)
# Line 272  c$$$         if(TRKVERBOSE)print*,i,' -- Line 273  c$$$         if(TRKVERBOSE)print*,i,' --
273  c$$$      ENDDO  c$$$      ENDDO
274  c$$$      print*,'chi2 -- ',DCHI2  c$$$      print*,'chi2 -- ',DCHI2
275    
276          IF(ISTEP.LT.3) GOTO 10 ! ***PP***
277        DO I=1,5        DO I=1,5
278           IF(ABS(DAL(I)).GT.ALTOL(I))GOTO 10 !>>>> new step!           IF(ABS(DAL(I)).GT.ALTOL(I))GOTO 10 !>>>> new step!
279        ENDDO        ENDDO
# Line 696  c$$$            if(.TRUE.)print*,'step', Line 698  c$$$            if(.TRUE.)print*,'step',
698           AYV(I)=DATAN(VOUT(5)/VOUT(6))*180./ACOS(-1.)           AYV(I)=DATAN(VOUT(5)/VOUT(6))*180./ACOS(-1.)
699  *     -----------------------------------------------  *     -----------------------------------------------
700    
701             IF(TRACKMODE.EQ.1) THEN
702    *     -----------------------------------------------
703    *        change of energy by bremsstrahlung for electrons
704                VOUT(7) = VOUT(7) * 0.997 !0.9968
705    *     -----------------------------------------------
706             ENDIF
707    
708        ENDDO        ENDDO
709    
710  c$$$      print*,'POSXY (dopo) ',vout  c$$$      print*,'POSXY (dopo) ',vout
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23