/[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.7 by pam-fi, Wed Nov 15 14:54:06 2006 UTC revision 1.8 by pam-fi, Tue Nov 21 14:00:40 2006 UTC
# Line 41  c      DATA XGOOD,YGOOD/nplanes*1.,nplan Line 41  c      DATA XGOOD,YGOOD/nplanes*1.,nplan
41  *                               !the tracking procedure  *                               !the tracking procedure
42        DATA STEPMAX/100./        !maximum number of steps in the trackin gprocess        DATA STEPMAX/100./        !maximum number of steps in the trackin gprocess
43    
44        DATA ALMAX/dinf,dinf,1.,dinf,dinf/ !limits on alpha vector components  c      DATA ALMAX/dinf,dinf,1.,dinf,dinf/ !limits on alpha vector components
45        DATA ALMIN/-dinf,-dinf,-1.,-dinf,-dinf/ !"  c      DATA ALMIN/-dinf,-dinf,-1.,-dinf,-dinf/ !"
46          DATA ALMAX/dinf,dinf,dinf,dinf,dinf/ !limits on alpha vector components
47          DATA ALMIN/-dinf,-dinf,-dinf,-dinf,-dinf/ !"
48    
49        DIMENSION DAL(5)                    !increment of vector alfa        DIMENSION DAL(5)                    !increment of vector alfa
50        DIMENSION CHI2DD_R(4,4),CHI2D_R(4) !hessiano e gradiente di chi2        DIMENSION CHI2DD_R(4,4),CHI2D_R(4) !hessiano e gradiente di chi2
# Line 798  c      IMPLICIT DOUBLE PRECISION (A-H,O- Line 800  c      IMPLICIT DOUBLE PRECISION (A-H,O-
800           ENDIF           ENDIF
801        ENDDO        ENDDO
802        CALL TRICIRCLE(NP,XP,ZP,AP,RP,CHI,XC,ZC,RADIUS,IFLAG)        CALL TRICIRCLE(NP,XP,ZP,AP,RP,CHI,XC,ZC,RADIUS,IFLAG)
803  c      print*,' circle: ',XC,ZC,RADIUS,' --- ',CHI        print*,' circle: ',XC,ZC,RADIUS,' --- ',CHI,IFLAG
804        IF(IFLAG.NE.0)GOTO 10 !straigth fit        IF(IFLAG.NE.0)GOTO 10 !straigth fit
805          if(CHI.gt.100)GOTO 10 !straigth fit
806        ARG = RADIUS**2-(ZINI-ZC)**2        ARG = RADIUS**2-(ZINI-ZC)**2
807        IF(ARG.LT.0)GOTO 10       !straigth fit        IF(ARG.LT.0)GOTO 10       !straigth fit
808        DC = SQRT(ARG)              DC = SQRT(ARG)      
# Line 808  c      print*,' circle: ',XC,ZC,RADIUS,' Line 811  c      print*,' circle: ',XC,ZC,RADIUS,'
811        AX = -(ZINI-ZC)/DC        AX = -(ZINI-ZC)/DC
812        DEF=100./(RADIUS*0.3*0.43)        DEF=100./(RADIUS*0.3*0.43)
813        IF(XC.GT.0)DEF=-DEF        IF(XC.GT.0)DEF=-DEF
814          
815          IF(ABS(X0).GT.30)THEN
816             PRINT*,'STRANGE GUESS: XC,ZC,R ',XC,ZC,RADIUS
817         $     ,' - CHI ',CHI,' - X0,AX,DEF ',X0,AX,DEF
818             GOTO 10       !straigth fit
819          ENDIF
820        GOTO 20                   !guess is ok        GOTO 20                   !guess is ok
821    
822  *     ----------------------------------------  *     ----------------------------------------
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.23