def _function(self, x, a, b, c, d):
        """Damped harmonic oscillator PS function
        x: frequency coordinate
        a: peak position
        b: half width
        c: area proportional parameter
        d: base line
        """
        return c/((a**2-x**2)**2 + (b*x)**2)+d