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

File: [HallC] / Analyzer / HTRACKING / h_fcnchisq.f (download)
Revision: 1.2, Wed Apr 13 19:38:06 1994 UTC (30 years, 5 months ago) by cdaq
Branch: MAIN
CVS Tags: oct1194, nov2894, may1894, mar1495, jun1794, jun1594, jun1394, jun0794, jul0794, jan2795, dec0694, aug1794, aug0394, apr1594, apr0695
Changes since 1.1: +6 -3 lines
(SAW) Change name of dummy arg to dumarg

      subroutine H_FCNCHISQ(npar,grad,fval,ray,iflag,dumarg)
*     This subroutine calculates chi**2 for MINUIT for HMS. The
*     arguments are determined by MINUIT
*
*     d.f. geesaman             17 January 1994
* $Log: h_fcnchisq.f,v $
* Revision 1.2  1994/04/13 20:38:06  cdaq
* (SAW) Change name of dummy arg to dumarg
*
* Revision 1.1  1994/02/19  06:14:15  cdaq
* Initial revision
*
*
      implicit none
      external H_DPSIFUN
      real*8 H_DPSIFUN
      include "gen_data_structures.cmn"
      include "hms_tracking.cmn"
      include "hms_geometry.cmn"
*
*     input
      real*8 ray(*),grad(*),dumarg
      integer*4 npar,iflag
*     output
      real*8 fval                              ! value of chi2
*
*     local variables
      real*8 diff
      integer*4 ihit
      integer*4 hitnum,planenum

      fval=0.0d0
      do ihit=1,HNTRACK_HITS(htrack_fit_num,1)
         hitnum=HNTRACK_HITS(htrack_fit_num,ihit+1)
         planenum=HDC_PLANE_NUM(hitnum)
         diff=(dble(HDC_WIRE_COORD(hitnum))-H_DPSIFUN(ray,planenum))
     &        /dble(hdc_sigma(planenum))
         fval=fval+diff*diff
      enddo
      return
      end

Analyzer/Replay: Mark Jones, Documents: Stephen Wood
Powered by
ViewCVS 0.9.2-cvsgraph-1.4.0