In [1]:

    

import absorberspec


    

In [3]:

    

masterwave, allflux, allivar = absorberspec.rest_allspec_readin()


    

    




Reading /Users/Benjamin/AstroData/AllInOne/AIO_MgIIAbsorber_SDSS_SDSSRestFrame_Wave00900_01800A.fits.
Reading /Users/Benjamin/AstroData/AllInOne/AIO_MgIIAbsorber_SDSS_SDSSRestFrame_Wave01800_03600A.fits.
Reading /Users/Benjamin/AstroData/AllInOne/AIO_MgIIAbsorber_SDSS_SDSSRestFrame_Wave03600_07200A.fits.

In [4]:

    

absorbers = absorberspec.absorber_readin()


    

In [5]:

    

absorbers.shape


    

    
Out[5]:





(35752,)

In [13]:

    

wave_pos = np.array([2200., 2900.])
rest_loc = np.searchsorted(masterwave, wave_pos)
outwave = masterwave[rest_loc[0]:rest_loc[1]]
index = (where(absorbers['REW_MGII_2796']>0.8))[0]
tmpflux = allflux[rest_loc[0]:rest_loc[1], index]
tmpivar = allivar[rest_loc[0]:rest_loc[1], index]
fluxmean = np.zeros((tmpflux.shape)[0])
fluxmedian = np.zeros((tmpflux.shape)[0])
fluxflag = np.ones(fluxmedian.size)
for i in np.arange((tmpflux.shape)[0]):
    iuse = (np.where(tmpivar[i,:]>0))[0]
    fluxmedian[i] = np.median(tmpflux[i,iuse])
    fluxmean[i] = np.mean(tmpflux[i,iuse])


    

In [19]:

    

allflux[rest_loc[0]:rest_loc[1], 0]


    

    
Out[19]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.])

In [1]:

    

import fitsio


    

In [2]:

    

data = (fitsio.read('/Users/Benjamin/AstroData/AllInOne/AIO_MgIIAbsorber_SDSS_SDSSRestFrame_Wave01800_03600A.fits'))[0]


    

In [3]:

    

data.shape


    

    
Out[3]:





()

In [ ]: