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