In [1]:

    

%matplotlib inline
%config InlineBackend.figure_format='retina'
# %config InlineBackend.figure_format='svg'

from collections import namedtuple
import matplotlib as mpl
import matplotlib.pyplot as plt
from matplotlib.gridspec import GridSpec

from androcmd.plot import contour_hess
from androcmd.phat_sixcolor import Catalog
from androcmd.phat_sixcolor import make_f475w_f814w
from androcmd.phat_sixcolor import make_f475w_f110w
from androcmd.phat_sixcolor import make_f475w_f160w
from androcmd.phat_sixcolor import make_f814w_f110w
from androcmd.phat_sixcolor import make_f814w_f160w
from androcmd.phat_sixcolor import make_f110w_f160w

# for bootstrapping planes
from starfisher import ColorPlane
Lim = namedtuple('Lim', 'x y')
wfc3_bands = ['F475W', 'F814W', 'F110W', 'F160W']
dpix = 0.05
mag_lim = 30.

catalog = Catalog(23)


    

    




/Users/jsick/phat_data

In [2]:

    

plane = make_f475w_f814w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f475w_vega'] - catalog.data['f814w_vega'],
             catalog.data['f814w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

    




/Users/jsick/.virtualenvs/starfish/lib/python2.7/site-packages/matplotlib/contour.py:931: UserWarning: linewidths is ignored by contourf
  warnings.warn('linewidths is ignored by contourf')
/Users/jsick/.virtualenvs/starfish/lib/python2.7/site-packages/matplotlib/figure.py:387: UserWarning: matplotlib is currently using a non-GUI backend, so cannot show the figure
  "matplotlib is currently using a non-GUI backend, "






    

In [3]:

    

plane = make_f110w_f160w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f110w_vega'] - catalog.data['f160w_vega'],
             catalog.data['f160w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

In [4]:

    

plane = make_f475w_f160w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f475w_vega'] - catalog.data['f160w_vega'],
             catalog.data['f160w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

In [5]:

    

plane = make_f475w_f110w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f475w_vega'] - catalog.data['f110w_vega'],
             catalog.data['f110w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

In [8]:

    

plane = make_f814w_f110w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f814w_vega'] - catalog.data['f110w_vega'],
             catalog.data['f110w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

In [7]:

    

plane = make_f814w_f160w()
fig = plt.figure(figsize=(6, 6))
gs = GridSpec(1, 1, wspace=0.4, bottom=0.2)
ax = fig.add_subplot(gs[0, 0])
contour_hess(ax,
             catalog.data['f814w_vega'] - catalog.data['f160w_vega'],
             catalog.data['f160w_vega'],
             plane.x_span, plane.y_span,
             plot_args={'ms': 3})
plane.plot_mask(ax)
ax.set_xlabel(plane.x_label)
ax.set_ylabel(plane.y_label)
ax.set_xlim(plane.xlim)
ax.set_ylim(plane.ylim)
fig.show()


    

In [ ]: