ApJdataFrames 010: Luhman2007Title: The Stellar Population of the Chamaeleon I Star-forming Region
Authors: Kevin L Luhman
Data is from this paper:
http://iopscience.iop.org/0067-0049/173/1/104/fulltext/71585.tables.html
In [1]:
%pylab inline
import seaborn as sns
import warnings
warnings.filterwarnings("ignore")
In [2]:
import pandas as pd
In [3]:
names = ["Name_2MASS","Name_alt","Spectral_Type","T_eff","AJ","L_bol","IMF"]
tbl6 = pd.read_csv("http://iopscience.iop.org/0067-0049/173/1/104/fulltext/71585.tb6.txt",
sep='\t', names=names, na_values='\ldots')
tbl6.head()
Out[3]:
In [4]:
upper_limits = (tbl6.T_eff == "\\leq 2400")
measurements = ~upper_limits
print "There are {} upper limits and {} measurements".format(upper_limits.sum(), measurements.sum())
In [5]:
sns.set_context("notebook", font_scale=1.5)
In [6]:
plt.plot(tbl6.T_eff[measurements], tbl6.L_bol[measurements], '.')
#plt.errorbar([2400.0]*upper_limits.sum(), tbl6.L_bol[upper_limits], xerr=0.0001, yerr=None, xlolims=True)
plt.ylabel(r"$L/L_{sun}$")
plt.xlabel(r"$T_{eff} (K)$")
plt.yscale("log")
plt.title("Luhman 2007 Chamaeleon I Members")
plt.xlim(5000,2000)
Out[6]:
In [7]:
from astroquery.simbad import Simbad
import astropy.coordinates as coord
import astropy.units as u
In [8]:
customSimbad = Simbad()
customSimbad.add_votable_fields('otype', 'id', 'sptype')
In [9]:
query_list = '2MASS'+tbl6.Name_2MASS.dropna().values
result = customSimbad.query_objects(query_list, verbose=True)
In [11]:
print "There were {} sources queried, and {} sources found.".format(len(query_list), len(result))
if len(query_list) == len(result):
print "Hooray! Everything matched"
else:
print "Which ones were not found?"
In [12]:
result[0:2]
Out[12]:
!mkdir ../data/Luhman2007
In [14]:
tbl6.to_csv("../data/Luhman2007/tbl6.csv", sep="\t")
The end.