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