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:
N | |
---|---|
COSMOS | 2980 |
GEMS | 246 |
GOODS-S | 4455 |
Total | 7681 |
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 [ ]: