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