Examples for LLSSurvey (v1.5) ============================= .. code:: python # imports import seaborn as sns; sns.set_style("white") import imp import h5py from bokeh.io import output_notebook, show, hplot, output_file from bokeh.plotting import figure from bokeh.models import Range1d output_notebook() from linetools import utils as ltu from pyigm.surveys.llssurvey import LLSSurvey from pyigm.surveys import llssurvey as llss from pyigm.abssys.igmsys import AbsSubSystem from pyigm.abssys.lls import LLSSystem from pyigm.abssys import lls as pylls from pyigm.surveys import lls_literature as llit .. parsed-literal:: /Users/xavier/anaconda/lib/python2.7/site-packages/matplotlib/__init__.py:872: UserWarning: axes.color_cycle is deprecated and replaced with axes.prop_cycle; please use the latter. warnings.warn(self.msg_depr % (key, alt_key)) .. raw:: html
BokehJS successfully loaded.
LLS Tree (JXP) -- Likely to Deprecate ------------------------------------- .. code:: python #reload(llsu) lls_survey = LLSSurvey.from_flist('Lists/lls_metals.lst', tree=os.getenv('LLSTREE')) .. parsed-literal:: Read 165 files from Lists/lls_metals.lst in the tree /u/xavier/LLS/ HD-LLS (Prochaska+15, ApJS, 221, 22) ------------------------------------ Simple Init ~~~~~~~~~~~ .. code:: python # Includes metallicity PDFs from Fumagalli+16 reload(llss) hdlls = llss.LLSSurvey.load_HDLLS(load_sys=False) hdlls .. parsed-literal:: HD-LLS: Loading summary file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1.fits HD-LLS: Loading ions file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_ions.json HD-LLS: Loading metallicity file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1_dustnhi.hdf5 .. parsed-literal:: With Systems (takes ~2 min) ~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code:: python hdlls = LLSSurvey.load_HDLLS() hdlls .. parsed-literal:: HD-LLS: Loading summary file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1.fits HD-LLS: Loading ions file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_ions.json HD-LLS: Loading metallicity file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1_dustnhi.hdf5 Loading systems from /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_sys.tar.gz Skipping a likely folder: SYS .. parsed-literal:: WARNING: UnitsWarning: The unit 'Angstrom' has been deprecated in the FITS standard. Suggested: 10**-1 nm. [astropy.units.format.utils] WARNING:astropy:UnitsWarning: The unit 'Angstrom' has been deprecated in the FITS standard. Suggested: 10**-1 nm. .. parsed-literal:: linetools.lists.parse: Reading linelist --- /Users/xavier/local/Python/linetools/linetools/data/lines/morton03_table2.fits.gz linetools.lists.parse: Reading linelist --- /Users/xavier/local/Python/linetools/linetools/data/lines/morton00_table2.fits.gz linetools.lists.parse: Reading linelist --- /Users/xavier/local/Python/linetools/linetools/data/lines/verner96_tab1.fits.gz linetools.lists.parse: Reading linelist --- /Users/xavier/local/Python/linetools/linetools/data/lines/verner94_tab6.fits .. parsed-literal:: WARNING: UnitsWarning: '0.1nm' did not parse as fits unit: Numeric factor not supported by FITS [astropy.units.core] WARNING:astropy:UnitsWarning: '0.1nm' did not parse as fits unit: Numeric factor not supported by FITS .. parsed-literal:: linetools.lists.parse: Reading linelist --- /Users/xavier/local/Python/linetools/linetools/data/lines/EUV_lines.ascii read_sets: Using set file -- /Users/xavier/local/Python/linetools/linetools/lists/sets/llist_v1.0.ascii .. parsed-literal:: Spectra ======= .. code:: python hdlls = LLSSurvey.load_HDLLS(grab_spectra=True) .. parsed-literal:: HD-LLS: Loading summary file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1.fits HD-LLS: Loading ions file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_ions.json HD-LLS: Loading metallicity file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_DR1_dustnhi.hdf5 Loading systems from /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/HD-LLS_sys.tar.gz Skipping a likely folder: SYS HD-LLS: Using files in /Users/xavier/local/Python/pyigm/pyigm/data/LLS/HD-LLS/Spectra/ Simple Attributes ~~~~~~~~~~~~~~~~~ .. code:: python hdlls.nsys .. parsed-literal:: 157 .. code:: python hdlls.NHI[0:10] .. parsed-literal:: array([ 19.65, 20.05, 17.55, 19.1 , 20. , 19.05, 19.1 , 19.05, 19.25, 20.2 ]) .. code:: python hdlls.name[0:5] .. parsed-literal:: array([u'J000345-232346.5_z2.187', u'J003454.8+163920_z3.754', u'J004049.5-402514_z2.816', u'J010355.3-300946_z2.908', u'J010516.8-184642_z2.927'], dtype=' ZionAEjzvminvmaxflag_NlogNsig_logN km / skm / s int64int64int64float64float64float64float64int64float64float64 6100.03.72296-150.460978217225.539021783313.0610.0 6200.03.72296-150.460978217225.539021783214.680.01 62063.423.72296-74.4609782171225.53902178300.00.0 6400.03.72296-150.460978217225.539021783113.9060.021 13200.03.72296-150.460978217225.539021783112.7420.027 13300.03.72296-150.460978217225.539021783312.4710.0 14200.03.72296-150.460978217225.539021783114.0310.02 14400.03.72296-150.460978217225.539021783113.5680.018 24200.03.72296-150.460978217225.539021783313.420.0 28200.03.72296-150.460978217225.539021783313.4720.0 30200.03.72296-150.460978217225.539021783312.9710.0 .. code:: python CII_clms = hdlls.ions((6,2)) CII_clms[70:80] .. raw:: html <Table length=10>
nameZionAEjzvminvmaxflag_NlogNsig_logN
km / skm / s
unicode32int64int64int64float64float64float64float64int64float64float64
J111008.61+024458.1_z3.4760000.00.00.00.000.00.0
J111113.64-080402.47_z3.4810000.00.00.00.000.00.0
J111113.64-080402.47_z3.8116200.03.8112-21.021.0113.4910.047
J113130.41+604420.7_z2.3620000.00.00.00.000.00.0
J113418.96+574204.7_z3.4106200.03.4144.9584916128.798769771312.79967719080.0
J113621+005021_z3.2486200.03.24829-89.294322049685.7056779504113.9340.03
J115659.59+551308.1_z2.6160000.00.00.00.000.00.0
J115906.48+133737.7_z3.7236200.03.72296-150.460978217225.539021783214.680.01
J115940.7-003203.5_z1.9046200.01.9044-380.090.0215.3790.909
J120331.29+152254.7_z2.7086200.02.708-220.0303.0214.6270.009
.. code:: python gdCII = CII_clms['flag_N']>0 np.sum(gdCII) .. parsed-literal:: 103 Simple Plot ~~~~~~~~~~~ .. code:: python # NHI p_NHI = figure(plot_width=400, plot_height=400, title='HD-LLS NHI')#,background_fill="#E8DDCB") hist, edges = np.histogram(hdlls.NHI, range=(17.,20.4), density=True, bins=20) p_NHI.quad(top=hist, bottom=0, left=edges[:-1], right=edges[1:],fill_color='blue') p_NHI.xaxis.axis_label = 'N_HI' # z p_z = figure(plot_width=400, plot_height=400, title='HD-LLS z')#,background_fill="#E8DDCB") hist, edges = np.histogram(hdlls.zabs, range=(2.4, 4.5), density=True, bins=20) p_z.quad(top=hist, bottom=0, left=edges[:-1], right=edges[1:],fill_color='red') p_z.xaxis.axis_label = 'z_LLS' # Show show(hplot(p_NHI,p_z)) .. raw:: html
-------------- Metallicity PDFs (Fumagalli+16, MNRAS, 455, 4100) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. code:: python print('Mean metallicity of {:s} is {:g}'.format(hdlls._abs_sys[0], hdlls._abs_sys[0].metallicity.meanZH)) .. parsed-literal:: Mean metallicity of is -1.47315 SDSS LLS (Prochaska+10, ApJ, 718, 391) -------------------------------------- All ~~~ .. code:: python sdss_dr7_all = LLSSurvey.load_SDSS_DR7(sample='all') sdss_dr7_all .. parsed-literal:: SDSS-DR7: Loading LLS file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/SDSS/lls_dr7_stat_LLS.fits.gz SDSS-DR7: Loading QSOs file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/SDSS/lls_dr7_qsos_sn2050.fits.gz .. parsed-literal:: .. code:: python sdss_dr7_all.sightlines[0:5] .. raw:: html <QTable length=5>
DRPLATEFIBERMJDRADECIMAGUUGGRSNRFLG_LLSFLG_EXTRALLS_FLUXZT2ZT0ZLLSZEMFLG_QSO
str3int32int32int32float64float64float32float32float32float32float32int16int16float32float32float32float32float32int16
DR7118732452708129.12736.345786619.849923.61632.880790.7638171.46418300.03.454694.813230.03.400010
DR716974453142199.1668510.98245720.135925.7154.935950.5628781.52782200.03.454694.793190.03.400060
DR7129118252735226.3885341.32954619.157823.75573.84410.6317772.45875200.04521413.449563.901373.304473.400180
DR7174229253053145.0249910.98884519.52224.93174.6010.6132771.33537300.03.454694.759930.03.400440
DR7213759954206220.7820328.057319.851525.75785.37150.4490011.64206400.03.454694.761260.03.40060
Stat ~~~~ .. code:: python sdss_dr7_stat = LLSSurvey.load_SDSS_DR7() .. parsed-literal:: SDSS-DR7: Loading LLS file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/SDSS/lls_dr7_stat_LLS.fits.gz SDSS-DR7: Loading QSOs file /Users/xavier/local/Python/pyigm/pyigm/data/LLS/SDSS/lls_dr7_qsos_sn2050.fits.gz SDSS-DR7: Performing stats (~60s) SDSS-DR7: Loaded .. code:: python len(sdss_dr7_stat.NHI) .. parsed-literal:: 254 Sky Map ~~~~~~~ .. code:: python all_coord = sdss_dr7_all.coord stat_coord = sdss_dr7_stat.coord .. code:: python %matplotlib inline .. code:: python plt.figure(figsize=(12, 7)) plt.clf() # Setup ax = plt.axes(projection='mollweide') ax = plt.axes() ax.set_xlabel('RA') ax.set_ylabel('DEC') ax.set_xticklabels(np.arange(30,331,30)) ax.grid(True) # All plt.scatter((all_coord.ra.value-180.)*np.pi/180., all_coord.dec.value*np.pi/180., marker='o', s=3., lw=0.5, edgecolors='gray', facecolors='none', label='all') # Stat plt.scatter((stat_coord.ra.value-180.)*np.pi/180., stat_coord.dec.value*np.pi/180., marker='o', s=3., lw=0.5, edgecolors='blue', facecolors='blue', label='stat') # Legend legend = plt.legend(loc='upper right', scatterpoints=1, borderpad=0.2, handletextpad=0.1, fontsize='large') plt.show() .. image:: LLSSurvey_examples_files/LLSSurvey_examples_33_0.png -------------- :math:`z \sim 2.5` HST (O'Meara et al. 2013, ApJ, 765, 137) ----------------------------------------------------------- ACS ~~~ .. code:: python acs = LLSSurvey.load_HST_ACS() acs .. parsed-literal:: HST-ACS: Loaded .. parsed-literal:: WFC3 ~~~~ .. code:: python wfc3 = LLSSurvey.load_HST_WFC3() wfc3 .. parsed-literal:: HST-WFC3: Loaded .. parsed-literal:: Combined ~~~~~~~~ .. code:: python HST_LLS = wfc3 + acs HST_LLS .. parsed-literal:: -------------- :math:`z \sim 3` MagE (Fumagalli et al. 2013, ApJ, 775, 78) ----------------------------------------------------------- Load ~~~~ .. code:: python z3mage = LLSSurvey.load_mage_z3() z3mage .. parsed-literal:: .. code:: python z3mage_NC = LLSSurvey.load_mage_z3(sample='non-color') z3mage_NC .. parsed-literal:: g(z) plot ~~~~~~~~~ .. code:: python zeval, gz = z3mage.calculate_gz() .. code:: python zeval[4000], gz[4000] .. parsed-literal:: (2.8705998897560931, 67) .. code:: python plt.clf() plt.plot(zeval, gz) plt.xlabel('z') plt.ylabel('g(z)') plt.show() .. image:: LLSSurvey_examples_files/LLSSurvey_examples_50_0.png Literature ---------- .. code:: python from pyigm.surveys.lls_literature import log_sum .. code:: python reload(llit) zonak04 = llit.zonak2004() print(zonak04) zonak04._ionN .. parsed-literal:: .. raw:: html <Table length=3>
logNsig_logNflag_NZion
float64float64int64int64int64
12.44735191210.1291787930711144
13.43294394280.1709971841741143
12.4152948020.051122
.. code:: python reload(llit) jenkins05 = llit.jenkins2005() jenkins05._ionN .. parsed-literal:: No error for N I No error for O I No error for O I WARNING: Using 1393.7550 Angstrom for your input 1393.7600 Angstrom WARNING: Using 1402.7700 Angstrom for your input 1402.7730 Angstrom WARNING: Using 1012.4950 Angstrom for your input 1012.5010 Angstrom WARNING: Using 1190.2030 Angstrom for your input 1190.1910 Angstrom .. raw:: html <Table length=13>
ionZsig_logNflag_NlogN
int64int64float64int64float64
3160.06114.19
110.05117.98
460.4113.9
2140.05113.95
270.4113.92
2160.2113.7
180.05114.47
360.0213.5
170.0312.44
4140.25113.53
1180.0312.6
260.0214.7
2260.12113.59
.. code:: python reload(llit) tripp05 = llit.tripp2005() tripp05._ionN .. raw:: html <Table length=7>
logNionsig_logNZflag_N
float64int64float64int64int64
13.3640.063
13.506683163920.0521636459228141
14.405923951410.07559112568281
13.2810.073
12.8840.0143
14.267382219320.26305521569361
13.4320.11261
.. code:: python reload(llit) peroux06a = llit.peroux06a() print(peroux06a) peroux06a._ionN .. parsed-literal:: .. raw:: html <Table length=5>
ionZsig_logNflag_NlogN
int64int64float64int64float64
2240.0312.6380948725
2300.0312.0878283736
1120.0333114495397112.2347702952
2260.0165813621851114.4774468535
2120.0213.6989700043
.. code:: python reload(llit) peroux06b = llit.peroux06b() print(peroux06b) peroux06b._ionN .. parsed-literal:: .. raw:: html <Table length=7>
ionZsig_logNflag_NlogN
int64int64float64int64float64
2240.161460798648113.0820669343
2250.0187835949381113.3100982719
2120.0215.1537904822
1120.0129070876912113.2627950666
2220.0656476491959112.4092739152
2300.044243954439113.4082908562
2260.0260564242229115.1549410447
.. code:: python reload(llit) meiring06 = llit.meiring06() print(meiring06) meiring06._ionN .. parsed-literal:: .. raw:: html <Table length=3>
ionZsig_logNflag_NlogN
int64int64float64int64float64
2220.0313.01
2240.0312.76
2300.0312.08
.. code:: python reload(llit) meiring07 = llit.meiring07() for imeiring07 in meiring07: print(imeiring07) print(imeiring07._ionN) .. parsed-literal:: ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 24 1 13.25 2 0.03 25 1 12.82 1 0.01 12 1 12.7 2 0.01 26 1 15.03 2 0.0 12 2 14.39 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 20 1 11.42 1 0.02 12 1 12.06 2 0.12 26 1 13.43 2 0.0 12 2 13.71 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.03 26 1 14.37 2 0.04 25 1 12.28 1 0.05 12 1 11.8 3 0.03 13 1 12.74 2 0.0 12 2 13.87 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.06 26 1 14.33 2 0.04 20 1 12.1 1 0.02 12 1 12.43 3 0.02 13 1 13.0 2 0.0 12 2 14.32 ion Z sig_logN flag_N logN --- --- -------- ------ ----- 2 13 0.0 2 13.72 3 13 0.02 1 13.29 2 14 0.15 1 14.86 1 12 0.02 1 12.49 2 12 0.0 2 14.26 2 26 0.02 1 14.5 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 26 1 14.36 3 0.14 13 1 12.74 1 0.08 12 1 12.21 2 0.0 13 2 13.56 2 0.0 12 2 14.25 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.05 26 1 14.44 3 0.14 13 1 12.86 1 0.03 12 1 12.47 2 0.0 13 2 13.33 2 0.0 12 2 14.39 .. code:: python reload(llit) meiring08 = llit.meiring08() for imeiring08 in meiring08: print(imeiring08) print(imeiring08._ionN) .. parsed-literal:: ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.11 25 1 12.57 2 0.0 13 2 13.6 3 0.03 13 1 13.15 2 0.05 14 1 15.05 2 0.0 12 2 14.44 2 0.09 26 1 14.96 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.01 26 1 14.35 2 0.1 25 1 12.31 1 0.02 12 1 12.54 3 0.04 13 1 13.65 2 0.0 12 2 14.3 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 1 0.07 12 1 12.04 2 0.08 26 1 13.49 2 0.0 12 2 13.59 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.11 26 1 14.35 2 0.0 12 2 14.25 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.01 26 1 14.35 1 0.05 12 1 12.1 3 0.12 13 1 12.68 2 0.0 12 2 14.37 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 1 0.05 12 1 12.32 2 0.02 26 1 14.55 2 0.0 12 2 14.36 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.07 26 1 13.76 1 0.09 12 1 12.37 3 0.15 13 1 12.52 2 0.0 12 2 14.12 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.01 20 1 12.79 2 0.03 25 1 13.05 1 0.0 12 2 12.9 2 0.02 26 1 14.92 2 0.0 12 2 14.26 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.17 24 1 12.75 2 0.1 25 1 12.2 2 0.0 13 2 13.71 3 0.01 13 1 13.67 2 0.12 14 1 14.84 1 0.02 12 1 12.51 2 0.0 12 2 14.45 2 0.02 26 1 14.61 .. code:: python reload(llit) nestor08 = llit.nestor08() print(nestor08) nestor08._ionN .. parsed-literal:: .. raw:: html <Table length=2>
ionZsig_logNflag_NlogN
int64int64float64int64float64
2240.0312.59
2300.0312.13
.. code:: python reload(llit) meiring09 = llit.meiring09() for imeiring09 in meiring09: print(imeiring09) print(imeiring09._ionN) .. parsed-literal:: ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.02 26 1 13.75 1 0.04 12 1 12.24 3 0.06 13 1 13.13 2 0.0 12 2 14.13 ion Z sig_logN flag_N logN --- --- -------- ------ ----- 2 13 0.0 2 13.07 3 13 0.04 1 12.83 2 14 0.04 1 14.45 1 12 0.05 1 11.75 2 12 0.0 2 13.81 2 26 0.01 1 14.25 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 26 1 14.68 2 0.0 14 2 14.86 1 0.06 12 1 12.26 2 0.0 13 2 13.68 2 0.0 12 2 14.5 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 14 1 13.56 2 0.02 26 1 13.33 2 0.03 13 1 12.21 2 0.0 12 2 13.25 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.1 20 1 12.22 1 0.05 12 1 12.44 2 0.03 26 1 14.17 2 0.0 12 2 14.04 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.1 24 1 13.08 2 0.04 25 1 13.24 2 0.0 13 2 13.77 3 0.02 13 1 13.52 2 0.12 14 1 15.31 1 0.02 12 1 12.67 2 0.04 20 1 12.59 2 0.0 12 2 14.44 2 0.02 26 1 15.0 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 1 0.14 12 1 11.84 2 0.02 26 1 13.47 2 0.0 12 2 14.05 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.26 24 1 12.94 3 0.05 13 1 13.41 2 0.06 14 1 15.49 1 0.03 12 1 12.6 2 0.0 12 2 14.58 2 0.01 26 1 14.88 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.09 26 1 14.19 3 0.03 13 1 13.43 1 0.07 12 1 12.32 2 0.02 13 1 13.42 2 0.0 12 2 14.39 .. code:: python reload(llit) dessauges09 = llit.dessauges09() for ills in dessauges09: print(ills) print(ills._ionN) .. parsed-literal:: ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 26 1 14.32 2 0.0 30 3 10.87 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.05 26 1 13.17 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.14 26 1 14.62 2 0.0 30 3 11.6 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.08 26 1 14.57 2 0.11 30 1 12.1 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.1 26 1 14.48 2 0.16 30 1 12.02 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.15 26 1 15.1 2 0.15 30 1 12.86 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 26 1 13.45 2 0.0 30 3 10.83 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.06 26 1 14.48 2 0.0 30 3 11.42 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.09 26 1 14.37 2 0.15 30 1 12.35 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.08 26 1 15.06 2 0.0 30 3 12.38 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.07 26 1 15.1 2 0.0 30 3 12.49 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.1 26 1 14.7 2 0.1 30 1 12.18 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.02 26 1 13.71 2 0.0 30 3 11.25 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.04 26 1 13.9 2 0.0 30 3 11.57 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.1 26 1 14.19 ion Z sig_logN flag_N logN --- --- -------- ------ ----- 2 26 0.0 2 15.22 2 30 0.09 1 12.96 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.11 26 1 13.5 2 0.0 30 3 12.03 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.14 26 1 14.96 2 0.0 30 3 11.94 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.24 26 1 14.28 2 0.0 30 3 11.1 .. code:: python reload(llit) tumlinson11 = llit.tumlinson11() print(tumlinson11) tumlinson11._ionN .. parsed-literal:: .. raw:: html <Table length=22>
logNionsig_logNZflag_N
float64int64float64int64int64
14.77023661430.062
14.952411889220.18879226753261
13.990661670320.0196862138948121
14.820.0163
14.040.1161
13.710.073
14.15390189120.142347356383141
14.208965947830.0142
11.310.0263
11.710.0143
12.080628722310.0336401521748121
13.650.073
14.410.063
14.624396880130.0623370581385161
13.927541022620.0242723547038261
11.75390189120.0835967450643201
14.399693455520.064051284841371
14.780628722310.13838773968581
14.55390189130.0711736781914261
15.089936511630.051306059610271
11.520.0223
14.94139690560.051785093050181
.. code:: python reload(llit) kacprzak12 = llit.kacprzak12() print(kacprzak12) kacprzak12._ionN .. parsed-literal:: .. raw:: html <Table length=8>
ionZsig_logNflag_NlogN
int64int64float64int64float64
2140.11113.16
1120.06111.54
460.05114.41
2120.07113.11
4140.0312.4
260.0213.39
1140.0311.8
680.05114.49
.. code:: python reload(llit) battisti12 = llit.battisti12() for ills in battisti12: print(ills) print(ills._ionN) .. parsed-literal:: ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.0 15 3 13.5 2 0.06 14 1 14.62 2 0.0 7 2 14.95 3 0.0 14 2 13.74 2 0.0 16 3 14.72 1 0.09 8 1 15.95 3 0.0 6 2 14.17 1 0.04 7 1 14.75 3 0.0 26 3 14.23 4 0.11 14 1 13.54 2 0.0 28 3 14.41 2 0.0 6 2 15.18 2 0.09 26 1 14.22 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 4 0.0 6 2 14.1 2 0.0 15 3 14.17 2 0.0 14 2 14.39 2 0.0 7 2 14.85 3 0.0 14 2 13.77 2 0.0 16 3 14.65 1 0.0 8 2 15.08 1 0.06 12 1 12.7 2 0.0 22 3 11.94 2 0.06 20 1 12.81 1 0.09 7 1 14.1 3 0.11 26 1 14.59 2 0.0 12 2 13.99 4 0.12 14 1 13.86 2 0.0 28 3 13.67 2 0.0 6 2 14.91 2 0.08 26 1 14.9 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.07 15 1 12.81 2 0.03 14 1 14.73 2 0.0 7 2 14.95 3 0.0 14 2 14.0 2 0.0 16 3 14.53 3 0.09 26 1 14.14 1 0.02 8 1 15.64 3 0.0 6 2 14.52 1 0.17 7 1 13.65 2 0.04 26 1 14.3 2 0.0 28 3 14.23 2 0.0 6 2 15.06 6 0.05 8 1 14.34 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 2 0.0 15 3 12.79 2 0.0 14 2 14.11 3 0.0 14 2 13.41 2 0.12 16 1 14.6 1 0.0 8 2 15.58 3 0.0 6 2 14.33 1 0.14 7 1 14.6 3 0.0 26 3 13.69 4 0.13 14 1 13.2 2 0.0 28 3 13.99 2 0.0 6 2 14.52 2 0.08 26 1 14.2 ion sig_logN Z flag_N logN --- -------- --- ------ ----- 4 0.08 6 1 13.99 2 0.05 14 1 14.22 2 0.07 7 1 14.16 3 0.0 14 2 13.3 2 0.0 16 3 14.24 1 0.0 8 2 14.56 1 0.0 7 3 13.74 4 0.08 14 1 13.59 2 0.0 28 3 13.92 2 0.0 6 2 14.35 2 0.07 26 1 14.01 Load them All ~~~~~~~~~~~~~ .. code:: python lls_lit = llit.load_lls_lit() lls_lit .. parsed-literal:: No error for N I No error for O I No error for O I WARNING: Using 1393.7550 Angstrom for your input 1393.7600 Angstrom WARNING: Using 1402.7700 Angstrom for your input 1402.7730 Angstrom WARNING: Using 1012.4950 Angstrom for your input 1012.5010 Angstrom WARNING: Using 1190.2030 Angstrom for your input 1190.1910 Angstrom .. parsed-literal:: .. code:: python lls_lit.ref .. parsed-literal:: u'Zon04,Jen05,Tri05,Prx06a,Prx06b,Mei06,Mei07,Mei08,Nes08,Mei09,DZ09,Tum11,Kcz12,Bat12' Plot ~~~~ .. code:: python # NHI p_NHI = figure(plot_width=400, plot_height=400, title='Lit-LLS NHI')#,background_fill="#E8DDCB") hist, edges = np.histogram(lls_lit.NHI, range=(17.,20.4), density=True, bins=20) p_NHI.quad(top=hist, bottom=0, left=edges[:-1], right=edges[1:],fill_color='blue') p_NHI.xaxis.axis_label = 'N_HI' # z p_z = figure(plot_width=400, plot_height=400, title='Lit-LLS z')#,background_fill="#E8DDCB") hist, edges = np.histogram(lls_lit.zabs, range=(0., 4.5), bins=20) p_z.quad(top=hist, bottom=0, left=edges[:-1], right=edges[1:],fill_color='red') p_z.xaxis.axis_label = 'z_LLS' # Show show(hplot(p_NHI,p_z)) .. raw:: html