In [ ]:

    

from astropy.coordinates import angles


    

In [ ]:

    

import astropy.units as u
import astropy.coordinates as coord
from astropy.table import Table
import matplotlib.pyplot as pl
import numpy as np
%matplotlib inline


    

In [ ]:

    

c = coord.SkyCoord(ra=[167]*u.deg, dec=[-14]*u.deg)


    

In [ ]:

    

d = np.genfromtxt('/Users/adrian/projects/triand-rrlyrae/data/Catalina_all_RRLyr.txt', 
                  dtype=None, usecols=range(11), names=True)
d_c = coord.SkyCoord(ra=d['ra']*u.degree, dec=d['dec']*u.degree)


    

In [ ]:

    

sep = d_c.separation(c)


    

In [ ]:

    

ix = (sep < 4*u.degree) & (d['helio_dist'] > 10.) & (d['helio_dist'] < 25.)
d[ix]['helio_dist']


    

In [ ]:

    

pl.plot(d['ra'], d['dec'], ls='none', marker=',', alpha=1.)
pl.plot(d['ra'][ix], d['dec'][ix], ls='none', marker='o', ms=4.)
pl.xlim(175, 160)
pl.ylim(-25, 0)


    

In [ ]:

    




    

In [ ]: