*-- Author : John V. Morris
FUNCTION FM2SOV(T)
************************************************************************
* *
* Symmetric double Gaussian plus quadratic background. *
* *
* JVM 31/7/91 *
************************************************************************
*KEEP,FCOMF.
PARAMETER( NP = 8 ) ! number of fit parameters
COMMON/FCOMF/PAR(NP),STAG
*KEND.
* PAR(1 to 3) ; Background = PAR(1) + PAR(2)*T + PAR(3)*T**2
* PAR(4) ; Lower peak height
* PAR(5) ; Upper peak height
* PAR(6) ; Drift Velocity in microns/nsec
* PAR(7) ; Timing Resolution in nsec
* PAR(8) ; Asymmetry in nsec
RMS = 0.5*SQRT(5.)*PAR(7)
SEP = 2.0*STAG/PAR(6)
FM2SOV = 0.0
B = PAR(1) + PAR(2)*T + PAR(3)*T*T
IF( B.GT.0. ) THEN
FM2SOV = B
ENDIF
F1 = ( (T+SEP+PAR(8))/RMS )**2
IF( F1.LT.50. )THEN
FM2SOV = FM2SOV + PAR(4)*EXP( -0.5*F1 )
ENDIF
F2 = ( (T-SEP+PAR(8))/RMS )**2
IF( F2.LT.50. )THEN
FM2SOV = FM2SOV + PAR(5)*EXP( -0.5*F2 )
ENDIF
RETURN
END
*