R_h_left, t_edges_left = np.histogram(th_left, bins=bins, range=(
                0, np.pi), weights=RR_left,    density=False)