ApJdataFrames 013: Lada2006

Title: Spitzer Observations of IC 348: The Disk Population at 2-3 Million Years
Authors: Charles J Lada, August A Muench, Kevin L Luhman, Lori E Allen, Lee Hartmann, Tom Megeath, Philip Myers, Giovanni Fazio, Kenneth Wood, James Muzerolle, George Rieke, Nick Siegler, and Erick Young

Data is from this paper:
http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/204953.html


In [1]:
%pylab inline
import seaborn as sns


Populating the interactive namespace from numpy and matplotlib

In [2]:
import warnings
warnings.filterwarnings("ignore")

In [3]:
import pandas as pd

In [7]:
from astropy.io import ascii
from astropy.table import Table, join

Table 1 - Spitzer IRAC/MIPS IC348 catalog


In [5]:
tbl1 = ascii.read("http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/datafile1.txt")
tbl1[0:4]


Downloading http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/datafile1.txt [Done]
Out[5]:
<Table masked=True length=4>
SeqRAhRAmRAsDEdDEmDEs3.6IRACmage_3.6IRACmag4.5IRACmage_4.5IRACmag5.8IRACmage_5.8IRACmag8.0IRACmage_8.0IRACmag24MIPSmagerr24MIPSmagBLEND-IRACBLEND-MIPS
hminsdegarcminarcsecmagmagmagmagmagmagmagmagmagmag
int64int64int64float64int64int64float64float64float64float64float64float64float64float64float64float64float64int64int64
134434.21232946.696.740.016.540.026.580.026.50.030.89-9.030191--
234435.36432104.587.090.026.810.026.460.045.820.043.20.03----
334450.64932196.757.530.017.610.037.470.047.50.047.080.05----
434431.18832622.097.790.017.740.027.660.027.730.034.38-9.0----

Table 2 - SED Derived $\alpha_{IRAC}$ and $A_V$

But really... spectral types


In [6]:
tbl2 = ascii.read("http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/datafile2.txt")
tbl2[0:4]


Downloading http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/datafile2.txt [Done]
Out[6]:
<Table masked=True length=4>
SeqSpTypeAValphaIRAC-oerralphaIRAC-oalphaIRAC-derralphaIRAC-dNbandsSTARDISK
mag
int64string40float64float64float64float64float64int64string48
1B53.1-2.6380.102-2.7080.0914STAR
2A23.2-1.3960.127-1.4680.1434THICK
3A03.9-2.7940.11-2.8830.1154STAR
4F02.3-2.7860.091-2.8380.0824STAR

In [11]:
join_tbls = join(tbl1, tbl2, keys="Seq")
print "There are {} rows in tbl1, {} in tbl2, and {} in the joined table.".format(len(tbl1), len(tbl2), len(join_tbls))


There are 307 rows in tbl1, 307 in tbl2, and 307 in the joined table

In [12]:
join_tbls[0:4]


Out[12]:
<Table masked=True length=4>
SeqRAhRAmRAsDEdDEmDEs3.6IRACmage_3.6IRACmag4.5IRACmage_4.5IRACmag5.8IRACmage_5.8IRACmag8.0IRACmage_8.0IRACmag24MIPSmagerr24MIPSmagBLEND-IRACBLEND-MIPSSpTypeAValphaIRAC-oerralphaIRAC-oalphaIRAC-derralphaIRAC-dNbandsSTARDISK
hminsdegarcminarcsecmagmagmagmagmagmagmagmagmagmagmag
int64int64int64float64int64int64float64float64float64float64float64float64float64float64float64float64float64int64int64string40float64float64float64float64float64int64string48
134434.21232946.696.740.016.540.026.580.026.50.030.89-9.030191--B53.1-2.6380.102-2.7080.0914STAR
234435.36432104.587.090.026.810.026.460.045.820.043.20.03----A23.2-1.3960.127-1.4680.1434THICK
334450.64932196.757.530.017.610.037.470.047.50.047.080.05----A03.9-2.7940.11-2.8830.1154STAR
434431.18832622.097.790.017.740.027.660.027.730.034.38-9.0----F02.3-2.7860.091-2.8380.0824STAR

Table 3 - Convenient passbands table


In [14]:
names = ["PASSBAND","DATA SYSTEM","REFERENCES","center_wavelength","F_{nu} (Jy)","References"]
tbl3 = pd.read_csv("http://iopscience.iop.org/1538-3881/131/3/1574/fulltext/204953.tb3.txt", 
                   na_values="\ldots", names = names, sep='\t')
tbl3.head()


Out[14]:
PASSBAND DATA SYSTEM REFERENCES center_wavelength F_{nu} (Jy) References
0 V_{L} Landolt 1, 2, 3 0.5423 3723 4
1 R_{C} Cousins 1, 2, 3 0.6410 3064 5
2 I_{C} Cousins 1, 2, 3 0.7890 2416 5
3 J 2MASS 6, 7, 8, 9 1.2350 1594 10
4 H 2MASS 6, 7, 8, 9 1.6620 1024 10

The end.