docs: def _g_a (self, x, a, b, s):
          """Asymmetric width term
          x: frequency coordinate
          a: peak position
          b: half width
          s: asymmetry parameter
          """
          return 2*b/(1.0+np.exp(s*(x-a)))