def plot_iso_set(AMags, EMags, **kwargs): ax = plt.gca() Acolor = AMags["WFC3_UVIS_F275W"] - AMags["WFC3_UVIS_F814W"] Amag = AMags["WFC3_UVIS_F814W"] Ecolor = EMags["WFC3_UVIS_F275W"] - EMags["WFC3_UVIS_F814W"] Emag = EMags["WFC3_UVIS_F814W"] ax.plot(Acolor, Amag, label="Population A", color='black', **kwargs) ax.plot(Ecolor, Emag, '--', label="Population E", color='black', **kwargs)