Examples for reading SDSS Quasar Database



In [41]:

    
# imports
from xastropy.sdss import quasars as sdssq
from xastropy.sdss import qso as sdssqso

Generate the Class



In [37]:

    
reload(sdssq)
sdss_dr7 = sdssq.SdssQuasars()









    



SDSS_QUASAR: Using summary file /u/xavier/SDSS//DR7_QSO/dr7_qso.fits.gz



In [38]:

    
sdss_dr7









    Out[38]:





[SdssQuasars: Version=DR7]

Grab by Plate/Fiber or Name



In [35]:

    
row = sdss_dr7.get_qso((287,264))
row









    Out[35]:




<QTable masked=False length=1>

Z Z_ERR Z_CONF Z_WARN PLATE MJD FIBERID FLG_TARG RAOBJ DECOBJ PSF_U PSF_G PSF_R PSF_I PSF_Z PSF_SU PSF_SG PSF_SR PSF_SI PSF_SZ
float32 float32 float32 int16 int32 int32 int32 int16 float64 float64 float32 float32 float32 float32 float32 float32 float32 float32 float32 float32
0.331188 0.00179699 0.0 0 287 52023 264 0 182.68148 -0.65201896 20.0461 20.1288 19.46 19.2803 18.6533 0.049576 0.03244 0.034281 0.038232 0.042613



In [36]:

    
## By name
row = sdss_dr7.get_qso('J000009.42-102751.9')
row









    Out[36]:




<QTable masked=False length=1>

Z Z_ERR Z_CONF Z_WARN PLATE MJD FIBERID FLG_TARG RAOBJ DECOBJ PSF_U PSF_G PSF_R PSF_I PSF_Z PSF_SU PSF_SG PSF_SR PSF_SI PSF_SZ
float32 float32 float32 int16 int32 int32 int32 int16 float64 float64 float32 float32 float32 float32 float32 float32 float32 float32 float32 float32
1.84493 0.00199866 0.0 0 650 52143 199 0 0.03926445 -10.46441 19.2416 19.0187 18.9659 18.7754 18.7046 0.035869 0.027196 0.021356 0.018248 0.046806

Generate a QSO Object



In [80]:

    
## Simple
reload(sdssqso)
qso = sdssqso.SdssQso()
qso









    Out[80]:





[SdssQso: 00:00:00 +00:00:00, z=0]



In [82]:

    
## From the Database Class
reload(sdssq)
reload(sdssqso)
sdss_dr7 = sdssq.SdssQuasars()
qso = sdss_dr7[(287,264)]
qso









    



SDSS_QUASAR: Using summary file /u/xavier/SDSS//DR7_QSO/dr7_qso.fits.gz






    Out[82]:





[SdssQso: 12:10:43.5552 -00:39:07.2683, z=0.331188]

Load Spectrum



In [73]:

    
qso.load_spec()
#qso.spec.plot()









    



SdssQso: Loading spectrum from /u/xavier/SDSS//DR7_QSO/spectro/1d_26/0287/1d/spSpec-52023-0287-264.fit



In [70]:

Sandbox



In [3]:

    
path = os.getenv('SDSSPATH')+'/DR7_QSO/' 
summf = path+'dr7_qso.fits.gz'
datdir = path+'spectro/1d_26/'
sdss = QTable.read(summf)
sdss[0:3]









    Out[3]:




<QTable masked=False length=3>

Z Z_ERR Z_CONF Z_WARN PLATE MJD FIBERID FLG_TARG RAOBJ DECOBJ PSF_U PSF_G PSF_R PSF_I PSF_Z PSF_SU PSF_SG PSF_SR PSF_SI PSF_SZ
float32 float32 float32 int16 int32 int32 int32 int16 float64 float64 float32 float32 float32 float32 float32 float32 float32 float32 float32 float32
0.331188 0.00179699 0.0 0 287 52023 264 0 182.68148 -0.65201896 20.0461 20.1288 19.46 19.2803 18.6533 0.049576 0.03244 0.034281 0.038232 0.042613
0.332475 0.000957323 0.0 0 277 51908 70 0 166.25829 -0.74849605 19.8404 19.7542 19.2126 18.9581 18.4499 0.043795 0.0422 0.034018 0.026109 0.036945
0.37957 0.00156732 0.0 0 303 51615 54 0 213.86771 -1.0875114 20.0745 19.9341 19.4695 18.9985 18.2755 0.04626 0.025193 0.019835 0.029198 0.030035



In [5]:

    
tmp = 'DR7'
tmp.lower()









    Out[5]:





'dr7'



In [ ]:

Z	Z_ERR	Z_CONF	Z_WARN	PLATE	MJD	FIBERID	FLG_TARG	RAOBJ	DECOBJ	PSF_U	PSF_G	PSF_R	PSF_I	PSF_Z	PSF_SU	PSF_SG	PSF_SR	PSF_SI	PSF_SZ
float32	float32	float32	int16	int32	int32	int32	int16	float64	float64	float32	float32	float32	float32	float32	float32	float32	float32	float32	float32
0.331188	0.00179699	0.0	0	287	52023	264	0	182.68148	-0.65201896	20.0461	20.1288	19.46	19.2803	18.6533	0.049576	0.03244	0.034281	0.038232	0.042613