Returns half-Cauchy random variates.

def rhalf_cauchy(alpha, beta, size=None):
    """
    Returns half-Cauchy random variates.
    """

    return abs(alpha + beta * np.tan(pi * random_number(size) - pi / 2.0))