*-- 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 *