Analyzer/HTRACKING/h_fcnchisq.f - view

Return to h_fcnchisq.f CVS log

Up to [HallC] / Analyzer / HTRACKING

File: [HallC] / Analyzer / HTRACKING / h_fcnchisq.f (download)
Revision: 1.3, Mon May 22 18:39:10 1995 UTC (29 years, 4 months ago) by cdaq
Branch: MAIN
CVS Tags: spring03, sep0596, sep-26-2002, sep-25-2002, sep-24-2002, sep-09-2002, sane, pionct, online07, online04, online03, oct1199, nov2696, mduality, may2495, mar-24-2003, jul2895, jan2496, jan1896, jan1796, gep_online, gep3, fpi2, emc, e01004, dec0198, bigcal, baryon, aug-12-2003, apr3096, apr-02-2003, Initial-CVS-Release, HEAD, Extra_Shower_Tubes_on_HMS_not_SOS
Changes since 1.2: +4 -1 lines

(SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts"

      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.3  1995/05/22 19:39:10  cdaq
* (SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts"
*
* 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 "hms_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