R_h_right, t_edges_right = np.histogram(th_right, bins=bins, range=(
                0, np.pi), weights=RR_right, density=False)