docs: def _function(self, x, a, b, c, d):
          """Gaussian PDF function
          x: coordinate
          a: peak position
          b: deviation (sigma)
          c: area proportional parameter
          d: base line
          """
          return c/b*np.sqrt(2*np.pi)*np.exp(-(x-a)**2/(2*b**2))+d