with plt.style.context(pubStyle):
    fig, ax = plt.subplots(1,1,figsize=(10,7))
    ax.scatter(HUGScolor, HUGSmag,s=1,c=HUGSdensity, cmap="Reds")
    ax.scatter(ACScolor, ACSmag, s=1, c=ACSdensity, cmap='Blues', alpha=0.25)
    ax.set_xlim(0,2)
    ax.set_ylim(16, 25)
    ax.invert_yaxis()
    
    ax.set_xlabel("F606W - F814W", fontsize=27)
    ax.set_ylabel("F606W", fontsize=27)
    
    
    legend_elements = [
        Line2D([0], [0], marker='o', color='w', label='HUGS', markerfacecolor='r', markersize=15),
        Line2D([0], [0], marker='o', color='w', label='ACS', markerfacecolor='b', markersize=15),
    ]
    ax.legend(handles = legend_elements, frameon=False, fontsize=27)
    fig.savefig("Figures/photometricOffset.pdf")