bestFitPhotF1 = optimizationResults['phot']["F275W"]
bestFitPhotF2 = optimizationResults['phot']["F814W"]

comp = order_best_fit_result(optimizationResults['r'])
AShiftParams = comp['A'][0]
EShiftParams = comp['E'][0]

bestFitPhotColor = bestFitPhotF1-bestFitPhotF2
bestFitPhotMag = bestFitPhotF2
testA = optimizationResults['r']['A'][AShiftParams[3]][AShiftParams[4]]
testB = optimizationResults['r']['E'][0.24][1.901]
fig, ax = plt.subplots(1,1,figsize=(10,7))
ax.scatter(bestFitPhotColor, bestFitPhotMag, s=1, c='yellow', alpha=0.1)
ax.plot(testA['fiducial'][:,0], testA['fiducial'][:,1], 'o', color='green')
ax.plot(testB['fiducial'][:,0], testB['fiducial'][:,1], 'o', color='red')

iso = testA['iso']
isoAAtAge, _ = interp_isochrone_age(iso, AShiftParams[1][0])

isoAColor = isoAAtAge['WFC3_UVIS_F275W_MAG'] - isoAAtAge["WFC3_UVIS_F814W_MAG"]
isoAMag = isoAAtAge["WFC3_UVIS_F814W_MAG"]
isoAAptMag, isoAAptColor = shift_isochrone(isoAColor, isoAMag, AShiftParams[1][1], AShiftParams[1][2])

ax.plot(isoAAptColor, isoAAptMag, color='tab:cyan')
ax.invert_yaxis()