ApJdataFramesTitle: IDENTIFICATION, CLASSIFICATIONS, AND ABSOLUTE PROPERTIES OF 773 ECLIPSING BINARIES FOUND IN THE TRANS-ATLANTIC EXOPLANET SURVEY
Authors: Jonathan Devor, David Charbonneau, Francis T O'Donovan, Georgi Mandushev, and Guillermo Torres
Data is from this paper:
http://iopscience.iop.org/article/10.1088/0004-6256/135/3/850/
In [1]:
import pandas as pd
In [2]:
from astropy.io import ascii, votable, misc
In [3]:
#! mkdir ../data/Devor2008
In [4]:
#! curl http://iopscience.iop.org/1538-3881/135/3/850/suppdata/aj259648_mrt7.txt >> ../data/Devor2008/aj259648_mrt7.txt
In [5]:
! du -hs ../data/Devor2008/aj259648_mrt7.txt
Not too big at all.
In [6]:
dat = ascii.read('../data/Devor2008/aj259648_mrt7.txt')
In [7]:
! head ../data/Devor2008/aj259648_mrt7.txt
In [8]:
dat.info
Out[8]:
In [9]:
df = dat.to_pandas()
In [10]:
df.head()
Out[10]:
In [11]:
df.columns
Out[11]:
In [12]:
sns.distplot(df.Per, norm_hist=False, kde=False)
Out[12]:
In [13]:
gi = (df.RAh == 4) & (df.RAm == 16) & (df.DEd == 28) & (df.DEm == 7)
In [14]:
gi.sum()
Out[14]:
In [15]:
df[gi].T
Out[15]:
The source is named T-Tau0-01262
The light curve files have the following 3-column format:
Column 1 - the Heliocentric Julian date (HJD), minus 2400000
Column 2 - normalized r-band magnitude
Column 3 - magnitude uncertainty
In [16]:
! head ../data/Devor2008/T-Tau0-01262.lc
In [17]:
cols = ['HJD-2400000', 'r_band', 'r_unc']
lc_raw = pd.read_csv('../data/Devor2008/T-Tau0-01262.lc', names=cols, delim_whitespace=True)
In [18]:
lc_raw.head()
Out[18]:
In [19]:
lc_raw.count()
Out[19]:
In [20]:
sns.set_context('talk')
In [21]:
plt.plot(lc_raw['HJD-2400000'], lc_raw.r_band, '.')
plt.ylim(0.6, -0.6)
Out[21]:
In [22]:
plt.plot(np.mod(lc_raw['HJD-2400000'], 3.375)/3.375, lc_raw.r_band, '.', alpha=0.5)
plt.xlabel('phase')
plt.ylabel('$\Delta \;\; r$')
plt.ylim(0.6, -0.6)
Out[22]:
In [23]:
plt.plot(np.mod(lc_raw['HJD-2400000'], 6.74215), lc_raw.r_band, '.')
plt.ylim(0.6, -0.6)
Out[23]:
The Devor et al. period is just twice the photometric period of 3.375 days.
Are those large vertical drops flares?
In [30]:
! ls /Users/gully/Downloads/catalog/T-Tau0-* | head -n 10
In [34]:
lc2 = pd.read_csv('/Users/gully/Downloads/catalog/T-Tau0-00397.lc', names=cols, delim_whitespace=True)
plt.plot(lc2['HJD-2400000'], lc2.r_band, '.')
plt.ylim(0.6, -0.6)
Out[34]:
In [46]:
this_p = df.Per[df.Name == 'T-Tau0-00397']
plt.plot(np.mod(lc2['HJD-2400000'], this_p), lc2.r_band, '.', alpha=0.5)
plt.xlabel('phase')
plt.ylabel('$\Delta \;\; r$')
plt.ylim(0.6, -0.6)
Out[46]:
In [ ]: