ApJdataFramesTitle: THE SDSS-III APOGEE SPECTRAL LINE LIST FOR H-BAND SPECTROSCOPY
Authors: M Shetrone, D Bizyaev, J E Lawler, C Allende Prieto, J A Johnson, V V Smith, K Cunha, J. Holtzman, A E García Pérez, Sz Mészáros, J Sobeck, O Zamora, D A Garcia Hernandez, D Souto, D Chojnowski, L Koesterke, S Majewski, and G Zasowski
Data is from this paper:
http://iopscience.iop.org/0067-0049/221/2/24/
In [1]:
import pandas as pd
In [2]:
from astropy.io import ascii, votable, misc
In [16]:
#! mkdir ../data/Shetrone2015
#! wget http://iopscience.iop.org/0067-0049/221/2/24/suppdata/apjs521087t7_mrt.txt
#! mv apjs521087t7_mrt.txt ../data/Shetrone2015/
#! du -hs ../data/Shetrone2015/apjs521087t7_mrt.txt
The file is about 24 MB.
In [19]:
dat = ascii.read('../data/Shetrone2015/apjs521087t7_mrt.txt')
In [20]:
! head ../data/Shetrone2015/apjs521087t7_mrt.txt
In [29]:
dat.info
Out[29]:
In [22]:
df = dat.to_pandas()
In [23]:
df.head()
Out[23]:
In [24]:
df.columns
Out[24]:
In [69]:
sns.distplot(df.Wave, norm_hist=False, kde=False)
Out[69]:
In [27]:
df.count()
Out[27]:
In [34]:
sns.lmplot('orggf', 'newgf', df, fit_reg=False)
Out[34]:
In [36]:
from astropy import units as u
In [37]:
u.cm
Out[37]:
In [49]:
EP1 = df.EP1.values*1.0/u.cm
EP2 = df.EP2.values*1.0/u.cm
In [50]:
EP1_eV = EP1.to(u.eV, equivalencies=u.equivalencies.spectral())
EP2_eV = EP2.to(u.eV, equivalencies=u.equivalencies.spectral())
In [51]:
deV = EP1_eV - EP2_eV
In [52]:
sns.distplot(deV)
Out[52]:
In [67]:
plt.plot(df.Wave, deV, '.', alpha=0.05)
plt.xlabel('$\lambda (\AA)$')
plt.ylabel('$\Delta E \;(\mathrm{eV})$')
Out[67]:
There are finite differences between the difference in the energy levels and the emitted wavelength based on other properties of the transition.