ax = par_set.plot_geometry()
par_set.plot_mass_centers(ax=ax)