def METHOD_NAME(self, f, d, pos=True):
    """Get the distribution of target positions
    Args:
        f (array):
            Distribution function evaluated
        d (float):
            Star distance
        pos (boolean):
            North or south
    Returns:
        tuple:
            array:
                Right ascension values
            array:
                Declination values
            array:
                Distances of the star
    """
    n = int(len(f))
    #        flips = np.arange(1,n,2) if f[0] == 0 else np.arange(0,n,2)
    flips = np.arange(0, n, 2)
    # angular separations from reference star
    psi = np.pi * f
    cosPsi = np.cos(psi)
    # calculating phi angle (i.e. DEC)
    sinPhi = np.abs(cosPsi) + (1 - np.abs(cosPsi)) * np.random.rand(n)
    phi = np.arcsin(sinPhi)  # only returns angles from 0 to pi/2
    # calculating phi angle (i.e. RA)
    cosTheta = cosPsi / sinPhi
    theta = np.arccos(cosTheta)
    # moving stars to southern hemisphere
    phi[flips] = np.pi - phi[flips]
    if pos:
        theta = 2 * np.pi - theta
    # final transforms
    dec = np.pi / 2.0 - phi
    dists = d * np.ones(n)
    return theta, dec, dists