*-- Author : John V. Morris
SUBROUTINE FVFITP(NC,IDF)
************************************************************************
* *
* Fit the Class I plot IDF to determine Drift Velocity and resolution *
* JVM 24/9/91 *
* *
* Scale errors on fitted parameters by SQRT(CHISQ/NDF) for bad Chisq!? *
* JVM 29/7/92 not very sanitary, but takes care of some systematics *
************************************************************************
*KEEP,FCOMF.
*KEEP,FVOUTP.
*KEEP,FHLUN.
*KEND.
* is there anything to fit ??
* WRITE(LUNH,1002)IDF,MXC
* WRITE(LUNH,1000)IDF,HSUM(IDF)
CALL UCOPY( P,PAR,NP )
CALL UCOPY( ST0,ST,NP )
CALL UCOPY( PMI0,PMI,NP )
CALL UCOPY( PMA0,PMA,NP )
CALL HFIT(IDF,FM2SOV,NP,PAR,CHISQ,IC,SIG,COV,ST,PMI,PMA)
* compute scaling factor for errors .. ?? ..
CALL HGIVE(IDF,TITL,NBIN,XMIN,XMAX,NY,YMI,YMA,NWT,LOC) ! get bin s
* compute resolution in microns and print results
* WRITE(LUNH,9600)PAR(6),SIG(6),PAR(7),SIG(7),RMIC(NC),EMIC(NC)
* & ,ASYM(NC),EASYM(NC)
*9600 FORMAT(/,1X,'Drift Velocity =',F10.2,' +/-',F6.2,' micr/nsec',
* & /,1X,'Resolution =',F10.2,' +/-',F6.2,' nsec ',
* & /,1X,' =',F10.2,' +/-',F6.2,' microns ',
* & /,1X,'Asymmetry =',F10.2,' +/-',F6.2,' nsec ')
*
* recompute RMIC in nsecs for compatibility with FSUMP
* compute number of fitted 4-hit track segments and error
* BIN = 3.0
* note than this constant includes the conversion factor of
* sqrt(5)/2 which converts PAR(7) to the actual RMS of the peaks.
* WRITE(LUNH,9601)SEGS(NC),ESEG(NC)
* WRITE(LUNH,9602)SEGF(NC),ESEGF(NC)
* WRITE(LUNH,9603)SEGN(NC),ESEGN(NC)