with plt.style.context(pubStyle):
    for n in tqdm(range(len(fitResults))):
        fig, ax = plt.subplots(1,1,figsize=(10,7))
        ax.plot(fiducial['Acolor'], fiducial['mag'], linewidth = 5, color=colorscheme.BonfireRed, label="E")
        ax.plot(fiducial['Bcolor'], fiducial['mag'], linewidth = 5, color=colorscheme.Violet, label="A", linestyle='dashed')
        ax.set_xlim(1.5, 5)
        ax.set_ylim(15.5, 22)
        # print(fitResults[n][0])
        ax.set_title(f"n: {n}, chi2: {fitResults[n][0]:0.3f}")
        A = fitResults[n][2][0]
        E = fitResults[n][2][1]

        Acolor = A['F275W'] - A["F814W"]
        Amag = A["F814W"]

        Ecolor = E["F275W"] - E["F814W"]
        Emag = E["F814W"]
        ax.plot(Acolor, Amag)
        ax.plot(Ecolor, Emag)
        ax.invert_yaxis()

        plt.savefig(f"Figures/bfr-{n}.png")
        plt.close(fig)