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