def METHOD_NAME(mean, sigma_low, sigma_high, size=1000):
    '''
    Samples from a skewed Gaussian distribution (hacky).
    Assumes different sigma left and right from the mean, samples them separately
    and combines into one sample, stitched together at the mean.
    Parameters
    ----------
    mean: float
        Mean of the Gaussian distribution
    sigma_low: float
        Standard deviation of the left Gaussian
    sigma_high: float
        Standard deviation of the right Gaussian
    size: int
        Number of samples to return
    Returns
    -------
    An array of the samples drawn from the distribution.
    '''
    samples_low = distl.gaussian(mean, sigma_low).sample(size)
    samples_high = distl.gaussian(mean, sigma_high).sample(size)
    samples_low_cut = samples_low[samples_low < mean]
    samples_high_cut = samples_high[samples_high >= mean]
    return np.hstack((samples_low_cut, samples_high_cut))