(file) Return to h_track_fit.f CVS log (file) (dir) Up to [HallC] / Analyzer / HTRACKING

Diff for /Analyzer/HTRACKING/h_track_fit.f between version 1.2 and 1.3

version 1.2, 1994/02/22 05:25:00 version 1.3, 1994/08/18 02:45:53
Line 9 
Line 9 
 *                              remove minuit. Make fit linear *                              remove minuit. Make fit linear
 *                              still does not do errors properly *                              still does not do errors properly
 * $Log$ * $Log$
   * Revision 1.3  1994/08/18 02:45:53  cdaq
   * (DM) Add calculation of residuals
   *
 * Revision 1.2  1994/02/22 05:25:00  cdaq * Revision 1.2  1994/02/22 05:25:00  cdaq
 * (SAW) Remove dfloat calls with floating args * (SAW) Remove dfloat calls with floating args
 * *
Line 20 
Line 23 
       include "gen_data_structures.cmn"       include "gen_data_structures.cmn"
       include "hms_tracking.cmn"       include "hms_tracking.cmn"
       include "hms_geometry.cmn"       include "hms_geometry.cmn"
         external H_DPSIFUN
 * *
 * *
 *     local variables *     local variables
Line 32 
Line 36 
       integer*4 ihit,ierr       integer*4 ihit,ierr
       integer*4 hit,plane       integer*4 hit,plane
       integer*4 i,j                             ! loop index       integer*4 i,j                             ! loop index
         real*8   H_DPSIFUN
         real*8   pos
       real*4   ray(hnum_fpray_param)       real*4   ray(hnum_fpray_param)
         real*8   ray1(4)
         real*8   ray2(4)
       real*8 TT(hnum_fpray_param)       real*8 TT(hnum_fpray_param)
       real*8 AA(hnum_fpray_param,hnum_fpray_param)       real*8 AA(hnum_fpray_param,hnum_fpray_param)
       real*8 dray(hnum_fpray_param)       real*8 dray(hnum_fpray_param)
Line 46 
Line 54 
 * *
       ABORT= .FALSE.       ABORT= .FALSE.
       ierr=0       ierr=0
   *  initailize residuals
   
         do itrack=1,HNTRACKS_MAX
           do plane=1,HMAX_NUM_DC_PLANES
             hdc_residual(itrack,plane)=1000
             hdc_sing_res(itrack,plane)=1000
           enddo
         enddo
   
   
 *     test for no tracks *     test for no tracks
       if(HNTRACKS_FP.ge.1) then       if(HNTRACKS_FP.ge.1) then
         do itrack=1,HNTRACKS_FP         do itrack=1,HNTRACKS_FP
           chi2= dummychi2           chi2= dummychi2
           htrack_fit_num=itrack           htrack_fit_num=itrack
   
 *     are there enough degrees of freedom *     are there enough degrees of freedom
           HNFREE_FP(itrack)=HNTRACK_HITS(itrack,1)-hnum_fpray_param           HNFREE_FP(itrack)=HNTRACK_HITS(itrack,1)-hnum_fpray_param
           if(HNFREE_FP(itrack).gt.0) then           if(HNFREE_FP(itrack).gt.0) then
   
 *     initialize parameters *     initialize parameters
            do i=1,hnum_fpray_param            do i=1,hnum_fpray_param
             TT(i)=0.             TT(i)=0.
Line 89 
Line 109 
       else       else
 *      calculate chi2 *      calculate chi2
         chi2=0.         chi2=0.
   
         ray(1)=dray(1)         ray(1)=dray(1)
         ray(2)=dray(2)         ray(2)=dray(2)
         ray(3)=dray(3)         ray(3)=dray(3)
Line 96 
Line 117 
          do ihit=2,HNTRACK_HITS(itrack,1)+1          do ihit=2,HNTRACK_HITS(itrack,1)+1
               hit=HNTRACK_HITS(itrack,ihit)               hit=HNTRACK_HITS(itrack,ihit)
               plane=HDC_PLANE_NUM(hit)               plane=HDC_PLANE_NUM(hit)
   
 * note chi2 is single precision * note chi2 is single precision
            chi2=chi2+             hdc_sing_res(itrack,plane)=
      &         (HDC_WIRE_COORD(hit)      &         (HDC_WIRE_COORD(hit)
      &         -hplane_coeff(remap(1),plane)*ray(1)      &         -hplane_coeff(remap(1),plane)*ray(1)
      &         -hplane_coeff(remap(2),plane)*ray(2)      &         -hplane_coeff(remap(2),plane)*ray(2)
      &         -hplane_coeff(remap(3),plane)*ray(3)      &         -hplane_coeff(remap(3),plane)*ray(3)
      &         -hplane_coeff(remap(4),plane)*ray(4))**2       &         -hplane_coeff(remap(4),plane)*ray(4))
      &         /(hdc_sigma(plane)*hdc_sigma(plane))             chi2=chi2+((hdc_sing_res(itrack,plane))**2
        &         )/(hdc_sigma(plane)*hdc_sigma(plane))
   
         enddo         enddo
       endif       endif
  
Line 116 
Line 140 
           HCHI2_FP(itrack)=chi2           HCHI2_FP(itrack)=chi2
         enddo                                 ! end loop over tracks         enddo                                 ! end loop over tracks
       endif       endif
   
   * calculate residuals for each chamber if in single stub mode
   * and there were 2 tracks found one in first chanber and one in the second
   
         if (hsingle_stub.ne.0) then
           if (HNTRACKS_FP.eq.2) then
             itrack=1
             ihit=2
             hit=HNTRACK_HITS(itrack,ihit)
             plane=HDC_PLANE_NUM(hit)
             if (plane.le.6) then
               itrack=2
               hit=HNTRACK_HITS(itrack,ihit)
               plane=HDC_PLANE_NUM(hit)
               if (plane.ge.7) then
   
   * condition of above met calculating residuals
   * assigning rays to tracks in each chamber
   * ray1 is ray from first chamber fit
   * ray2 is ray from second chamber fit
   
               ray1(1)=dble(HX_FP(1))
               ray1(2)=dble(HY_FP(1))
               ray1(3)=dble(HXP_FP(1))
               ray1(4)=dble(HYP_FP(1))
               ray2(1)=dble(HX_FP(2))
               ray2(2)=dble(HY_FP(2))
               ray2(3)=dble(HXP_FP(2))
               ray2(4)=dble(HYP_FP(2))
   
               itrack=1
   * loop over hits in second chamber
               do ihit=1,HNTRACK_HITS(itrack+1,1)
   
   * calculate residual in second chamber from first chamber track
                 hit=HNTRACK_HITS(itrack+1,ihit+1)
                 plane=HDC_PLANE_NUM(hit)
                 pos=H_DPSIFUN(ray1,plane)
                 hdc_residual(itrack,plane)=HDC_WIRE_COORD(hit)-pos
   
               enddo
   
               itrack=2
   * loop over hits in first chamber
               do ihit=1,HNTRACK_HITS(itrack-1,1)
   
   * calculate residual in first chamber form second chamber track
                 hit=HNTRACK_HITS(itrack-1,ihit+1)
                 plane=HDC_PLANE_NUM(hit)
                 pos=H_DPSIFUN(ray2,plane)
                 hdc_residual(itrack,plane)=HDC_WIRE_COORD(hit)-pos
   
               enddo
               endif
             endif
           endif
         endif
   
 *     test if we want to dump out trackfit results *     test if we want to dump out trackfit results
       if(hdebugtrackprint.ne.0) then       if(hdebugtrackprint.ne.0) then
          call h_print_tracks          call h_print_tracks
Line 124 
Line 206 
       end       end
  
 * *
   

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
Analyzer/Replay: Mark Jones, Documents: Stephen Wood
 Powered by
ViewCVS 0.9.2-cvsgraph-1.4.0