Ingesting Prochaska+11 CGM dataset [v1.0]¶
# imports
from astropy.table import Table
from astropy.coordinates import SkyCoord
import pyigm
from pyigm.cgm.cgmsurvey import CGMAbsSurvey
Summary¶
ovi_file = pyigm.__path__[0]+'/data/CGM/P11/lowovidat.fits'
ovidat = Table.read(ovi_file)
ovidat[0:3]
QSO | QSO_RA | QSO_DEC | QSO_ZEM | QSO_VMAG | QSO_UV | FLG_GAL | GAL_FIL | R_LIMIT | N_GAL [4] | COMPLETE [2,4] | FLG_FUSE | FUSE_EXP | FUSE_SNR | FLG_STIS | STIS_COMM [4] | FLG_GHRS | GHRS_COMM [4] |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
str15 | str15 | str15 | float32 | float32 | float32 | int16 | str58 | float32 | int32 | int32 | int16 | float32 | float32 | int16 | str60 | int16 | str60 |
Q0026+1259 | 00:29:13.8 | +13:16:04 | 0.142 | 14.78 | 20.0 | 1 | /u/xavier/LCO/OVI/FUSE/data/Q0026+1259/Q0026+1259_gal.fits | 20.0 | 131 .. 47 | 22 .. 85 | 1 | 20000.0 | 7.0 | 0 | .. | 1 | G270M 5222s .. |
TONS180 | 00:57:20.0 | -22:22:56 | 0.06198 | 16.6 | 30.0 | 1 | /u/xavier/LCO/OVI/FUSE/data/TONS180/TONS180_gal.fits | 19.7 | 110 .. 4 | 15 .. 92 | 1 | 132453.0 | 15.0 | 2 | G140M 7000s 15. .. | 0 | .. |
TONS210 | 01:21:51.5 | -28:20:57 | 0.116 | 14.7 | 70.0 | 1 | /u/xavier/LCO/OVI/FUSE/data/TONS210/TONS210_gal.fits | 20.0 | 71 .. 5 | 6 .. 87 | 1 | 56500.0 | 20.0 | 2 | E140M 22000s .. | 0 | .. |
qso_radec = SkyCoord(ra=ovidat['QSO_RA'], dec=ovidat['QSO_DEC'], unit=(u.hourangle, u.deg))
Dwarfs¶
cgm_dwarf_file = pyigm.__path__[0]+'/data/CGM/P11/dwarf_galabs_strct.fits'
cgm_dwarfs = Table.read(cgm_dwarf_file)
cgm_dwarfs[0:3]
FIELD | ID | OBJ_ID | FLG_ANLY | FLG_SURVEY | OBJ_TYPE | MAG [10] | MAGERR [10] | FILTER [10] | IMG_FIL [10] | XYPIX [2] | RA | DEC | AREA | STARGAL | GAL_TYPE | GAL_COEFF [10] | Z | VCIRC | FSPEC_FIL [10] | DRA | DDEC |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
str11 | int32 | str1 | int16 | int16 | int16 | float32 | float32 | str1 | str27 | float32 | float64 | float64 | float32 | float32 | str5 | float32 | float64 | float32 | str29 | float64 | float64 |
Q0026+1259 | 1303 | a | 7 | 1 | 0 | 19.733 .. 0.0 | 0.125 .. 9.99 | B .. | Images/Q0026+1259XB.fits .. | 1014.76 .. 1034.29 | 7.28886499155 | 13.2746021629 | 10.43 | 0.18 | Late | 0.331349 .. 0.0 | 0.0329451337457 | 0.0 | .. | 45.6990947951 | 0.0191390593786 |
TONS180 | 2295 | a | 7 | 1 | 0 | 18.923 .. 1.0 | 0.088 .. 0.05 | B .. | Images/TONS180XB.fits .. | 1318.89 .. 607.18 | 14.2667432785 | -22.44755991 | 10.92 | 0.19 | Late | -0.0115093 .. 0.0 | 0.0233643911779 | 0.0 | .. | 154.07390626 | 0.0207054292147 |
PKS0405-123 | 90033 | a | 7 | 1 | 0 | 0.0 .. 0.44 | 0.0 .. 0.06 | .. | .. | 61.9512 .. -12.1839 | 61.9512481689 | -12.1838884354 | 0.0 | 0.0 | Late | 0.0 .. 0.0 | 0.167 | 1.0 | .. | 97.7760719383 | 0.0877822129587 |
Funny columns – Renaming¶
cgm_dwarfs.rename_column('DRA', 'rho(kpc)')
cgm_dwarfs.rename_column('DDEC', 'L(L*)')
MAG has HI¶
2 == Ref
3 == z
4 == Lya EW
5 == sigEW Lya
8 == NHI
9 == sigNHI
cgm_dwarfs[1]['MAG']
array([ 1.89230003e+01, 1.80701466e+01, 5.00000000e+00,
2.34109219e-02, 2.22000000e+02, 2.90000000e+01,
5.60000000e+01, 4.00000000e+00, 1.38000002e+01,
1.00000000e+00], dtype=float32)
MAGERR has OVI¶
2 == Ref
3 == z
4 == EW 1031
5 == sigEW 1031
8 == NOVI
9 == sigNOVI
cgm_dwarfs[1]['MAGERR']
array([ 8.79999995e-02, 6.49999976e-02, 4.00000000e+00,
2.33999994e-02, 4.30000000e+01, 1.50000000e+01,
3.00000000e+01, 0.00000000e+00, 1.34799995e+01,
5.00000007e-02], dtype=float32)
Refs¶
refdict = {1: 'tripp08', 2: 'tc08a', 3: 'ds08', 4: 'dsr+06', 5: 'pss04', 6: 'cm09', 9: 'p+11'}
Ingest¶
python ingest_lit.py
Read¶
p11_tarfile = pyigm.__path__[0]+'/data/CGM/P11/P11_sys.tar'
p11 = CGMAbsSurvey.from_tarball(p11_tarfile, chk_lowz=False)
WARNING: UnitsWarning: The unit 'Angstrom' has been deprecated in the FITS standard. Suggested: 10**-1 nm. [astropy.units.format.utils]
/Users/xavier/local/Python/linetools/linetools/lists/linelist.py:374: RuntimeWarning: divide by zero encountered in log10
self._data['log(w*f)'] = np.log10(qm_strength)
/Users/xavier/anaconda/lib/python2.7/site-packages/numpy/ma/core.py:824: RuntimeWarning: invalid value encountered in less_equal
return umath.less_equal(x, self.critical_value)
Loading abundances from Asplund2009
Abundances are relative by number on a logarithmic scale with H=12
Skipping a likely folder: CGM_JSON
/Users/xavier/local/Python/linetools/linetools/isgm/abssystem.py:288: UserWarning: Input AbsComponent with Zion=(8, 6) does not match AbsSystem rules. Not appending
warnings.warn('Input AbsComponent with Zion={} does not match AbsSystem rules. Not appending'.format(abscomp.Zion))
/Users/xavier/local/Python/linetools/linetools/isgm/abssystem.py:294: UserWarning: Failed velocity overlap
warnings.warn('Failed velocity overlap')
/Users/xavier/local/Python/linetools/linetools/isgm/abssystem.py:288: UserWarning: Input AbsComponent with Zion=(1, 1) does not match AbsSystem rules. Not appending
warnings.warn('Input AbsComponent with Zion={} does not match AbsSystem rules. Not appending'.format(abscomp.Zion))
p11
<CGM_Survey: nsys=54, ref=>
<IGMSystem: IGMSystem 00:29:13.8 13:16:04, 0.0329451, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 00:29:13.8 13:16:04, 0.039311, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 00:57:20 -22:22:56, 0.0233644, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 00:57:20 -22:22:56, 0.045619, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 03:11:55.2 -76:51:51, 0.202643, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 03:11:55.2 -76:51:51, 0.0593531, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.0964516, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.297609, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.352, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.153212, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.203022, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.36124, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.167, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 04:07:48.4 -12:11:37, 0.167043, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 10:07:26.1 12:48:56, 0.0296659, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 10:07:26.1 12:48:56, 0.0092207, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 10:31:54.3 -14:16:51, 0.0508329, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 11:19:08.7 21:19:18, 0.165951, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 11:19:08.7 21:19:18, 0.0600208, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 11:19:08.7 21:19:18, 0.0593876, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 11:19:08.7 21:19:18, 0.13829, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:14:17.7 14:03:13, 0.06438, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:14:17.7 14:03:13, 0.0646084, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:14:17.7 14:03:13, 0.0519877, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:14:17.7 14:03:13, 0.0511341, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:19:20.9 06:38:38, 0.124102, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:19:20.9 06:38:38, 0.0131789, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:19:20.9 06:38:38, 0.00666757, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:19:20.9 06:38:38, 0.0080957, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:29:06.7 02:03:09, 0.00620912, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:33:25.8 09:31:23, 0.0118122, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:33:25.8 09:31:23, 0.125431, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:33:25.8 09:31:23, 0.0905207, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 12:33:25.8 09:31:23, 0.206801, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:05:33 -10:33:19, 0.191709, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:05:33 -10:33:19, 0.0935802, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:05:33 -10:33:19, 0.0425621, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:05:33 -10:33:19, 0.145303, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:09:47 08:19:49, 0.127621, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 13:09:47 08:19:49, 0.0337313, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 14:29:06.4 01:17:06, 0.0299413, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 14:29:06.4 01:17:06, 0.0281113, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 15:55:43 11:11:24, 0.0150682, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 15:55:43 11:11:24, 0.0395311, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 15:55:43 11:11:24, 0.0420751, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 15:55:43 11:11:24, 0.0416811, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.13262, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.0503651, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.155532, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.0807993, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.0788036, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.0809754, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:55:01.5 -09:22:25, 0.0516533, NHI=0, Z/H=0>
<IGMSystem: IGMSystem 21:58:51.8 -30:13:30, 0.0169383, NHI=0, Z/H=0>
p11.rho
\[[46.32007,~231.08717,~153.91919,~281.05024,~34.785816,~241.08236,~271.56657,~263.85542,~170.30404,~193.11905,~277.85574,~228.92416,~101.21597,~122.3351,~181.21741,~78.670472,~306.79467,~159.32049,~133.01256,~221.58367,~137.97841,~151.94839,~72.602911,~182.7867,~138.10197,~94.401827,~106.07147,~92.907301,~37.14687,~85.547486,~41.258645,~260.27253,~120.93128,~243.45857,~225.98147,~69.865838,~227.78876,~88.463791,~99.493909,~274.85217,~162.60742,~309.02131,~176.70156,~282.59707,~293.98086,~184.38542,~228.32783,~305.40178,~266.93147,~292.71387,~236.98361,~35.061178,~271.77462,~108.99879] \; \mathrm{kpc}\]