n_h, t_edges = np.histogram(theta, bins=bins, range=(
            0, np.pi), weights=None, density=False)