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

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

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

revision 1.6 by pam-fi, Mon Oct 2 13:12:48 2006 UTC revision 1.8 by pam-fi, Wed Oct 25 16:18:41 2006 UTC
# Line 242  c         iflag=0 Line 242  c         iflag=0
242           IDCAND = icand         !fitted track-candidate           IDCAND = icand         !fitted track-candidate
243           ifail=0                !error flag in chi2 computation           ifail=0                !error flag in chi2 computation
244           jstep=0                !# minimization steps           jstep=0                !# minimization steps
245             iprint=0
246           call mini_2(jstep,ifail)           if(DEBUG)iprint=1
247             call mini2(jstep,ifail,iprint)
248           if(ifail.ne.0) then           if(ifail.ne.0) then
249              if(DEBUG)then              if(DEBUG)then
250                 print *,                 print *,
251       $              '*** MINIMIZATION FAILURE *** (mini_2) '       $              '*** MINIMIZATION FAILURE *** (mini2) '
252       $              ,iev       $              ,iev
253              endif              endif
254              chi2=-chi2              chi2=-chi2
# Line 667  c      double precision xi_B,yi_B,zi_B Line 668  c      double precision xi_B,yi_B,zi_B
668              resxPAM = resxPAM*fbad_cog(2,icx)              resxPAM = resxPAM*fbad_cog(2,icx)
669           elseif(PFAx.eq.'ETA2')then           elseif(PFAx.eq.'ETA2')then
670  c            cog2 = cog(2,icx)  c            cog2 = cog(2,icx)
671  c            etacorr = pfa_eta2(cog2,viewx,nldx,angx)              c            etacorr = pfaeta2(cog2,viewx,nldx,angx)            
672  c            stripx = stripx + etacorr  c            stripx = stripx + etacorr
673              stripx = stripx + pfa_eta2(icx,angx)            !(3)              stripx = stripx + pfaeta2(icx,angx)            !(3)
674              resxPAM = risx_eta2(angx)                       !   (4)              resxPAM = risx_eta2(angx)                       !   (4)
675              if(DEBUG.and.fbad_cog(2,icx).ne.1)              if(DEBUG.and.fbad_cog(2,icx).ne.1)
676       $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)       $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)
677              resxPAM = resxPAM*fbad_cog(2,icx)              resxPAM = resxPAM*fbad_cog(2,icx)
678           elseif(PFAx.eq.'ETA3')then                         !(3)           elseif(PFAx.eq.'ETA3')then                         !(3)
679              stripx = stripx + pfa_eta3(icx,angx)            !(3)              stripx = stripx + pfaeta3(icx,angx)            !(3)
680              resxPAM = risx_eta3(angx)                       !   (4)              resxPAM = risx_eta3(angx)                       !   (4)
681              if(DEBUG.and.fbad_cog(3,icx).ne.1)              !(3)              if(DEBUG.and.fbad_cog(3,icx).ne.1)              !(3)
682       $           print*,'BAD icx >>> ',viewx,fbad_cog(3,icx)!(3)       $           print*,'BAD icx >>> ',viewx,fbad_cog(3,icx)!(3)
683              resxPAM = resxPAM*fbad_cog(3,icx)               !(3)              resxPAM = resxPAM*fbad_cog(3,icx)               !(3)
684           elseif(PFAx.eq.'ETA4')then                         !(3)           elseif(PFAx.eq.'ETA4')then                         !(3)
685              stripx = stripx + pfa_eta4(icx,angx)            !(3)              stripx = stripx + pfaeta4(icx,angx)            !(3)
686              resxPAM = risx_eta4(angx)                       !   (4)              resxPAM = risx_eta4(angx)                       !   (4)
687              if(DEBUG.and.fbad_cog(4,icx).ne.1)              !(3)              if(DEBUG.and.fbad_cog(4,icx).ne.1)              !(3)
688       $           print*,'BAD icx >>> ',viewx,fbad_cog(4,icx)!(3)       $           print*,'BAD icx >>> ',viewx,fbad_cog(4,icx)!(3)
689              resxPAM = resxPAM*fbad_cog(4,icx)               !(3)              resxPAM = resxPAM*fbad_cog(4,icx)               !(3)
690           elseif(PFAx.eq.'ETA')then                          !(3)           elseif(PFAx.eq.'ETA')then                          !(3)
691              stripx = stripx + pfa_eta(icx,angx)             !(3)              stripx = stripx + pfaeta(icx,angx)             !(3)
692              resxPAM = ris_eta(icx,angx)                     !   (4)              resxPAM = ris_eta(icx,angx)                     !   (4)
693              if(DEBUG.and.fbad_cog(2,icx).ne.1)              !(3)              if(DEBUG.and.fbad_cog(2,icx).ne.1)              !(3)
694       $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)!(3)       $           print*,'BAD icx >>> ',viewx,fbad_cog(2,icx)!(3)
# Line 731  c     $     print*,PFAx,icx,angx,stripx, Line 732  c     $     print*,PFAx,icx,angx,stripx,
732              resyPAM = resyPAM*fbad_cog(2,icy)              resyPAM = resyPAM*fbad_cog(2,icy)
733           elseif(PFAy.eq.'ETA2')then           elseif(PFAy.eq.'ETA2')then
734  c            cog2 = cog(2,icy)  c            cog2 = cog(2,icy)
735  c            etacorr = pfa_eta2(cog2,viewy,nldy,angy)  c            etacorr = pfaeta2(cog2,viewy,nldy,angy)
736  c            stripy = stripy + etacorr  c            stripy = stripy + etacorr
737              stripy = stripy + pfa_eta2(icy,angy)            !(3)              stripy = stripy + pfaeta2(icy,angy)            !(3)
738              resyPAM = risy_eta2(angy)                       !   (4)              resyPAM = risy_eta2(angy)                       !   (4)
739              resyPAM = resyPAM*fbad_cog(2,icy)              resyPAM = resyPAM*fbad_cog(2,icy)
740              if(DEBUG.and.fbad_cog(2,icy).ne.1)              if(DEBUG.and.fbad_cog(2,icy).ne.1)
741       $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)       $           print*,'BAD icy >>> ',viewy,fbad_cog(2,icy)
742           elseif(PFAy.eq.'ETA3')then                         !(3)           elseif(PFAy.eq.'ETA3')then                         !(3)
743              stripy = stripy + pfa_eta3(icy,angy)            !(3)              stripy = stripy + pfaeta3(icy,angy)            !(3)
744              resyPAM = resyPAM*fbad_cog(3,icy)               !(3)              resyPAM = resyPAM*fbad_cog(3,icy)               !(3)
745              if(DEBUG.and.fbad_cog(3,icy).ne.1)              !(3)              if(DEBUG.and.fbad_cog(3,icy).ne.1)              !(3)
746       $           print*,'BAD icy >>> ',viewy,fbad_cog(3,icy)!(3)       $           print*,'BAD icy >>> ',viewy,fbad_cog(3,icy)!(3)
747           elseif(PFAy.eq.'ETA4')then                         !(3)           elseif(PFAy.eq.'ETA4')then                         !(3)
748              stripy = stripy + pfa_eta4(icy,angy)            !(3)              stripy = stripy + pfaeta4(icy,angy)            !(3)
749              resyPAM = resyPAM*fbad_cog(4,icy)               !(3)              resyPAM = resyPAM*fbad_cog(4,icy)               !(3)
750              if(DEBUG.and.fbad_cog(4,icy).ne.1)              !(3)              if(DEBUG.and.fbad_cog(4,icy).ne.1)              !(3)
751       $           print*,'BAD icy >>> ',viewy,fbad_cog(4,icy)!(3)       $           print*,'BAD icy >>> ',viewy,fbad_cog(4,icy)!(3)
752           elseif(PFAy.eq.'ETA')then                          !(3)           elseif(PFAy.eq.'ETA')then                          !(3)
753              stripy = stripy + pfa_eta(icy,angy)             !(3)              stripy = stripy + pfaeta(icy,angy)             !(3)
754              resyPAM = ris_eta(icy,angy)                     !   (4)              resyPAM = ris_eta(icy,angy)                     !   (4)
755  c            resyPAM = resyPAM*fbad_cog(2,icy)              !(3)TEMPORANEO  c            resyPAM = resyPAM*fbad_cog(2,icy)              !(3)TEMPORANEO
756              resyPAM = resyPAM*fbad_eta(icy,angy)            !   (4)              resyPAM = resyPAM*fbad_eta(icy,angy)            !   (4)
# Line 1708  c      common/dbg/DEBUG Line 1709  c      common/dbg/DEBUG
1709              goto 10              goto 10
1710           endif           endif
1711  *     ----------------------------------------------------  *     ----------------------------------------------------
1712    *     cut on multiplicity (X VIEW)
1713    *     ----------------------------------------------------
1714             if(mult(icx).ge.mult_x_max)then
1715                cl_single(icx)=0
1716                goto 10
1717             endif
1718    *     ----------------------------------------------------
1719  *     cut BAD (X VIEW)              *     cut BAD (X VIEW)            
1720  *     ----------------------------------------------------  *     ----------------------------------------------------
1721           badseed=BAD(VIEW(icx),nvk(MAXS(icx)),nst(MAXS(icx)))           badseed=BAD(VIEW(icx),nvk(MAXS(icx)),nst(MAXS(icx)))
# Line 1753  c     endif Line 1761  c     endif
1761                 goto 20                 goto 20
1762              endif              endif
1763  *     ----------------------------------------------------  *     ----------------------------------------------------
1764    *     cut on multiplicity (X VIEW)
1765    *     ----------------------------------------------------
1766                if(mult(icy).ge.mult_y_max)then
1767                   cl_single(icy)=0
1768                   goto 20
1769                endif
1770    *     ----------------------------------------------------
1771  *     cut BAD (Y VIEW)              *     cut BAD (Y VIEW)            
1772  *     ----------------------------------------------------  *     ----------------------------------------------------
1773              badseed=BAD(VIEW(icy),nvk(MAXS(icy)),nst(MAXS(icy)))              badseed=BAD(VIEW(icy),nvk(MAXS(icy)),nst(MAXS(icy)))
# Line 2998  c     $                                 Line 3013  c     $                                
3013                                enddo                                enddo
3014                                ifail=0 !error flag in chi^2 computation                                ifail=0 !error flag in chi^2 computation
3015                                jstep=0 !number of  minimization steps                                jstep=0 !number of  minimization steps
3016                                call mini_2(jstep,ifail)                                iprint=0
3017                                  if(DEBUG)iprint=1
3018                                  call mini2(jstep,ifail,iprint)
3019                                if(ifail.ne.0) then                                if(ifail.ne.0) then
3020                                   if(DEBUG)then                                   if(DEBUG)then
3021                                      print *,                                      print *,
3022       $                              '*** MINIMIZATION FAILURE *** '       $                              '*** MINIMIZATION FAILURE *** '
3023       $                              //'(mini_2 in clouds_to_ctrack)'       $                              //'(mini2 in clouds_to_ctrack)'
3024                                   endif                                   endif
3025                                   chi2=-chi2                                   chi2=-chi2
3026                                endif                                endif
# Line 3752  c*************************************** Line 3769  c***************************************
3769    
3770  c      good2=1!.true.  c      good2=1!.true.
3771        chi2_nt(ntr)        = sngl(chi2)        chi2_nt(ntr)        = sngl(chi2)
3772        nstep_nt(ntr)       = 0!nstep        nstep_nt(ntr)       = nstep
3773    
3774        phi   = al(4)             !(4)        phi   = al(4)             !(4)
3775        sinth = al(3)             !(4)        sinth = al(3)             !(4)
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.23