with plt.style.context(pubStyle):
    fig, ax = plt.subplots(1,1,figsize=(10,7))
    ax.scatter(color, mag, s=1, alpha=0.1, color='grey')
    for aID, eID in tqdm(ids):
        Apick = sAR[aID]
        Spick = sER[eID]
        AMags, EMags, optA, optE, funA, funE = correct_iso_set(Apick, Spick)
        plot_iso_set(AMags, EMags, alpha=0.1)
    ax.set_xlim(1, 6)
    ax.set_ylim(15, 22)
    ax.invert_yaxis()
    # ax.legend()
    ax.set_xlabel("F275W-F814W", fontsize=27)
    ax.set_ylabel("F814W", fontsize=27)
    # fig.savefig("Figures/BestFitResults.pdf")