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