xmax2 = width + max_dist * np.cos(rng.uniform(0, 2 * np.pi))