In [1]:

    

# Import any necessary packages

from matplotlib import pyplot as plt
from astropy.io import fits
#from astropy.table import Table
#import requests
import numpy as np


    

In [2]:

    

# I did a positional match with a radius of 1.0 arcsec on 
# Table 6 in Willett+(in prep) and Table 2 in Simmons+(submitted). 

# Split by survey, these are the matches:


    

In [3]:

    

# Get the matched dataset

data = fits.getdata('/Users/willettk/Dropbox/gzhubble/gzh_gzcandels.fits',1)


    

In [ ]:

    

gzc = data['t00_smooth_or_featured_a1_features_weighted_frac']
gzh = data['t01_smooth_or_features_a02_features_or_disk_best_fraction']


    

In [ ]:

    

# Make a Simmons-style plot of the averaged vote fractions, akin to Figure 7 in the GZC paper.

fig,ax = plt.subplots(1,1)

ax.hexbin(gzc,gzh,cmap=plt.cm.gray)

plt.show()


    

In [ ]: