nA = 0 
nE = 4
with plt.style.context(pubStyle):
    A = fitResults[nA][2][0]
    E = fitResults[nE][2][1]
    print(fitResults[nA][0])
    print(fitResults[nE][0])

    Acolor = A['WFC3_UVIS_F275W'] - A["WFC3_UVIS_F814W"]
    Amag = A["WFC3_UVIS_F814W"]

    Ecolor = E["WFC3_UVIS_F275W"] - E["WFC3_UVIS_F814W"]
    Emag = E["WFC3_UVIS_F814W"]
    fig, ax = plt.subplots(1,1,figsize=(10,7))
    
    ASotedIDX = np.argsort(Amag)
    ax.plot(fiducial['Acolor'], fiducial['mag'], color=colorscheme.BonfireRed, linestyle='dashed', alpha=0.5, label="Fiducial E")
    ax.plot(savgol_filter(Acolor[ASotedIDX], 10, 2), Amag[ASotedIDX], color=colorscheme.BonfireRed, label="Best Fit E")
    ax.plot(fiducial['Bcolor'], fiducial['mag'], color=colorscheme.Violet, linestyle='dashed', alpha=0.5, label="Fiducial A")
    # ax.set_xlim(1.5, 5)
    # ax.set_ylim(18, 21)
    ax.plot(Ecolor, Emag, color=colorscheme.Violet, label="Best Fit A")
    ax.invert_yaxis()
    
    ax.legend(frameon=False)
    
#     print(fitResults[nA][1][0][3].x)
#     print(fitResults[nE][1][1][3].x)
    
#     print(fitResults[nA][1][0][1])
#     print(fitResults[nE][1][1][1])
    
#     print(fitResults[nA][1][0][2])
#     print(fitResults[nE][1][1][2])
    
    ax.set_xlabel("F275W - F814W", fontsize=23)
    ax.set_ylabel("F814W", fontsize=23)
    
    fig.savefig("Figures/BestFit.png", dpi=300)