Skip to content

Instantly share code, notes, and snippets.

@pllim

pllim/synphot_issue_159_postmortem2.ipynb

Created October 19, 2018 21:34
Show Gist options
  • Save pllim/2d053b327d384c587b0297950644e7d1 to your computer and use it in GitHub Desktop.
Save pllim/2d053b327d384c587b0297950644e7d1 to your computer and use it in GitHub Desktop.
Download ZIP
synphot issue 159 (4) without patch 166
Display the source blob
Display the rendered blob
Raw
synphot_issue_159_postmortem2.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import os\n",
"\n",
"import numpy as np\n",
"from astropy import constants as const\n",
"from astropy import units as u\n",
"from astropy.table import Table\n",
"\n",
"import pysynphot as S\n",
"import synphot\n",
"import stsynphot"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Parameters to use for calculations.\n",
"mjd = 57754\n",
"filters = ['f435w', 'f475w', 'f502n', 'f550m', 'f555w', \n",
" 'f606w', 'f625w', 'f658n', 'f660n', 'f814w', \n",
" 'f850lp', 'f892n']\n",
"spectra = {'Vega': 'alpha_lyr_stis_008.fits',\n",
" 'g191b2b': 'g191b2b_stisnic_006.fits'}"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Cache bandpass objects so we don't have to rebuild in loops.\n",
"bp_pysyn_dict = {}\n",
"bp_stsyn_dict = {}\n",
"\n",
"for filt in filters:\n",
" obsmode = 'acs,wfc1,{},mjd#{}'.format(filt, mjd)\n",
" \n",
" bp_pysyn = S.ObsBandpass(obsmode)\n",
" bp_pysyn_dict[filt] = bp_pysyn\n",
" \n",
" bp_stsyn = stsynphot.band(obsmode)\n",
" bp_stsyn_dict[filt] = bp_stsyn"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"# Table to store results.\n",
"t = Table(names=('source', 'obsmode',\n",
" 'plam_pysyn', 'plam_st', 'd_plam',\n",
" 'flam_pysyn', 'flam_st', 'd_flam',\n",
" 'st_pysyn', 'st_st', 'd_st',\n",
" 'pnu_pysyn', 'pnu_st', 'd_pnu',\n",
" 'fnu_pysyn', 'fnu_st', 'd_fnu',\n",
" 'jy_pysyn', 'jy_st', 'd_jy',\n",
" 'njy_pysyn', 'njy_st', 'd_njy',\n",
" 'ab_pysyn', 'ab_st', 'd_ab',\n",
" 'vmag_pysyn', 'vmag_st', 'd_vmag',\n",
" 'ct_pysyn', 'ct_st', 'd_ct',\n",
" 'cr_pysyn', 'cr_st', 'd_cr',\n",
" 'obmag_pysyn', 'obmag_st', 'd_ob'),\n",
" dtype=('S25', 'S25',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4',\n",
" 'f4', 'f4', 'f4'))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"for name, spec in spectra.items():\n",
" specfile = os.path.join(os.environ['PYSYN_CDBS'], 'calspec', spec)\n",
" sp_pysyn = S.FileSpectrum(specfile)\n",
" sp_stsyn = synphot.SourceSpectrum.from_file(specfile)\n",
" \n",
" for filt in filters:\n",
" # PySynphot\n",
" bp_pysyn = bp_pysyn_dict[filt]\n",
" obs_pysyn = S.Observation(sp_pysyn, bp_pysyn)\n",
" plam_pysyn = obs_pysyn.effstim()\n",
" flam_pysyn = obs_pysyn.effstim('flam')\n",
" st_pysyn = obs_pysyn.effstim('stmag')\n",
" pnu_pysyn = obs_pysyn.effstim('photnu')\n",
" fnu_pysyn = obs_pysyn.effstim('fnu')\n",
" jy_pysyn = obs_pysyn.effstim('jy')\n",
" njy_pysyn = obs_pysyn.effstim('njy')\n",
" ab_pysyn = obs_pysyn.effstim('abmag')\n",
" vmag_pysyn = obs_pysyn.effstim('vegamag')\n",
" ct_pysyn = obs_pysyn.effstim('count')\n",
" cr_pysyn = obs_pysyn.countrate()\n",
" obmag_pysyn = obs_pysyn.effstim('obmag')\n",
" \n",
" # synphot\n",
" bp_stsyn = bp_stsyn_dict[filt]\n",
" obs_stsyn = synphot.Observation(sp_stsyn, bp_stsyn, binset=bp_stsyn.binset)\n",
" plam_st = obs_stsyn.effstim()\n",
" flam_st = obs_stsyn.effstim('flam')\n",
" st_st = obs_stsyn.effstim('stmag')\n",
" pnu_st = obs_stsyn.effstim('photnu')\n",
" fnu_st = obs_stsyn.effstim('fnu')\n",
" jy_st = obs_stsyn.effstim('Jy')\n",
" njy_st = obs_stsyn.effstim('nJy')\n",
" ab_st = obs_stsyn.effstim('abmag')\n",
" vmag_st = obs_stsyn.effstim('vegamag', vegaspec=stsynphot.Vega)\n",
" ct_st = obs_stsyn.effstim('count', area=stsynphot.conf.area)\n",
" cr_st = obs_stsyn.countrate(stsynphot.conf.area)\n",
" obmag_st = obs_stsyn.effstim('obmag', area=stsynphot.conf.area)\n",
" \n",
" # Differences: Percentage for flux, delta for mag\n",
" d_plam = abs(plam_st.value - plam_pysyn) * 100 / plam_pysyn\n",
" d_flam = abs(flam_st.value - flam_pysyn) * 100 / flam_pysyn\n",
" d_st = abs(st_st.value - st_pysyn)\n",
" d_pnu = abs(pnu_st.value - pnu_pysyn) * 100 / pnu_pysyn\n",
" d_fnu = abs(fnu_st.value - fnu_pysyn) * 100 / fnu_pysyn\n",
" d_jy = abs(jy_st.value - jy_pysyn) * 100 / jy_pysyn\n",
" d_njy = abs(njy_st.value - njy_pysyn) * 100 / njy_pysyn\n",
" d_ab = abs(ab_st.value - ab_pysyn)\n",
" d_vmag = abs(vmag_st.value - vmag_pysyn)\n",
" d_ct = abs(ct_st.value - ct_pysyn) * 100 / ct_pysyn\n",
" d_cr = abs(cr_st.value - cr_pysyn) * 100 / cr_pysyn\n",
" d_ob = abs(obmag_st.value - obmag_pysyn)\n",
" \n",
" t.add_row((name, filt,\n",
" plam_pysyn, plam_st, d_plam,\n",
" flam_pysyn, flam_st, d_flam,\n",
" st_pysyn, st_st, d_st,\n",
" pnu_pysyn, pnu_st, d_pnu,\n",
" fnu_pysyn, fnu_st, d_fnu,\n",
" jy_pysyn, jy_st, d_jy,\n",
" njy_pysyn, njy_st, d_njy,\n",
" ab_pysyn, ab_st, d_ab,\n",
" vmag_pysyn, vmag_st, d_vmag,\n",
" ct_pysyn, ct_st, d_ct,\n",
" cr_pysyn, cr_st, d_cr,\n",
" obmag_pysyn, obmag_st, d_ob))"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"#t['plam_pysyn'].format = '.4f'\n",
"#t['plam_st'].format = '.4f'\n",
"t['d_plam'].format = '.4f'\n",
"\n",
"#t['flam_pysyn'].format = '.4f'\n",
"#t['flam_st'].format = '.4f'\n",
"t['d_flam'].format = '.4f'\n",
"\n",
"t['st_pysyn'].format = '.4f'\n",
"t['st_st'].format = '.4f'\n",
"t['d_st'].format = '.4f'\n",
"\n",
"#t['pnu_pysyn'].format = '.4f'\n",
"#t['pnu_st'].format = '.4f'\n",
"t['d_pnu'].format = '.4f'\n",
"\n",
"#t['fnu_pysyn'].format = '.4f'\n",
"#t['fnu_st'].format = '.4f'\n",
"t['d_fnu'].format = '.4f'\n",
"\n",
"#t['jy_pysyn'].format = '.4f'\n",
"#t['jy_st'].format = '.4f'\n",
"t['d_jy'].format = '.4f'\n",
"\n",
"#t['njy_pysyn'].format = '.4f'\n",
"#t['njy_st'].format = '.4f'\n",
"t['d_njy'].format = '.4f'\n",
"\n",
"t['ab_pysyn'].format = '.4f'\n",
"t['ab_st'].format = '.4f'\n",
"t['d_ab'].format = '.4f'\n",
"\n",
"t['vmag_pysyn'].format = '.4f'\n",
"t['vmag_st'].format = '.4f'\n",
"t['d_vmag'].format = '.4f'\n",
"\n",
"#t['ct_pysyn'].format = '.4f'\n",
"#t['ct_st'].format = '.4f'\n",
"t['d_ct'].format = '.4f'\n",
"\n",
"#t['cr_pysyn'].format = '.4f'\n",
"#t['cr_st'].format = '.4f'\n",
"t['d_cr'].format = '.4f'\n",
"\n",
"t['obmag_pysyn'].format = '.4f'\n",
"t['obmag_st'].format = '.4f'\n",
"t['d_ob'].format = '.4f'"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode plam_pysyn plam_st d_plam\n",
"------- ------- ------------ ------------ ------\n",
" Vega f435w 1388.9543 1389.7065 0.0542\n",
" Vega f475w 1262.1 1247.3917 1.1654\n",
" Vega f502n 1174.1124 1174.0823 0.0026\n",
" Vega f550m 954.75726 953.2209 0.1609\n",
" Vega f555w 1024.2319 1016.5271 0.7523\n",
" Vega f606w 854.60974 833.35205 2.4874\n",
" Vega f625w 744.1318 737.2676 0.9224\n",
" Vega f658n 585.37463 585.5664 0.0328\n",
" Vega f660n 638.09894 638.15344 0.0085\n",
" Vega f814w 460.78983 455.1136 1.2318\n",
" Vega f850lp 374.4472 372.78082 0.4450\n",
" Vega f892n 389.4863 389.56317 0.0197\n",
"g191b2b f435w 0.037219618 0.036767814 1.2139\n",
"g191b2b f475w 0.029373676 0.028755099 2.1059\n",
"g191b2b f502n 0.02471785 0.02471605 0.0073\n",
"g191b2b f550m 0.01838102 0.018336529 0.2420\n",
"g191b2b f555w 0.02076922 0.020515472 1.2218\n",
"g191b2b f606w 0.01620097 0.015612465 3.6325\n",
"g191b2b f625w 0.01307267 0.012908479 1.2560\n",
"g191b2b f658n 0.010810098 0.010810378 0.0026\n",
"g191b2b f660n 0.011000068 0.010999745 0.0029\n",
"g191b2b f814w 0.0066147675 0.006489574 1.8926\n",
"g191b2b f850lp 0.0046830717 0.0046371873 0.9798\n",
"g191b2b f892n 0.004802293 0.004801422 0.0181\n"
]
}
],
"source": [
"# PHOTLAM\n",
"t['source', 'obsmode', 'plam_pysyn', 'plam_st', 'd_plam'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode flam_pysyn flam_st d_flam\n",
"------- ------- -------------- -------------- ------\n",
" Vega f435w 6.3583143e-09 6.3583143e-09 0.0000\n",
" Vega f475w 5.2617617e-09 5.2617617e-09 0.0000\n",
" Vega f502n 4.643267e-09 4.643267e-09 0.0000\n",
" Vega f550m 3.3965732e-09 3.3965732e-09 0.0000\n",
" Vega f555w 3.786644e-09 3.786644e-09 0.0000\n",
" Vega f606w 2.8481941e-09 2.8481941e-09 0.0000\n",
" Vega f625w 2.3368611e-09 2.3368611e-09 0.0000\n",
" Vega f658n 1.7661446e-09 1.7661446e-09 0.0000\n",
" Vega f660n 1.9207116e-09 1.9207116e-09 0.0000\n",
" Vega f814w 1.1339126e-09 1.1339126e-09 0.0000\n",
" Vega f850lp 8.21991e-10 8.21991e-10 0.0000\n",
" Vega f892n 8.6786933e-10 8.6786933e-10 0.0000\n",
"g191b2b f435w 1.7038287e-13 1.7038287e-13 0.0000\n",
"g191b2b f475w 1.224604e-13 1.224604e-13 0.0000\n",
"g191b2b f502n 9.775178e-14 9.775178e-14 0.0000\n",
"g191b2b f550m 6.5390944e-14 6.5390944e-14 0.0000\n",
"g191b2b f555w 7.678498e-14 7.678498e-14 0.0000\n",
"g191b2b f606w 5.3993658e-14 5.3993658e-14 0.0000\n",
"g191b2b f625w 4.1053234e-14 4.1053234e-14 0.0000\n",
"g191b2b f658n 3.2615346e-14 3.2615346e-14 0.0000\n",
"g191b2b f660n 3.311079e-14 3.311079e-14 0.0000\n",
"g191b2b f814w 1.6277635e-14 1.6277635e-14 0.0000\n",
"g191b2b f850lp 1.02803355e-14 1.02803355e-14 0.0000\n",
"g191b2b f892n 1.0700667e-14 1.0700667e-14 0.0000\n"
]
}
],
"source": [
"# FLAM\n",
"t['source', 'obsmode', 'flam_pysyn', 'flam_st', 'd_flam'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode st_pysyn st_st d_st \n",
"------- ------- -------- ------- ------\n",
" Vega f435w -0.6084 -0.6084 0.0000\n",
" Vega f475w -0.4028 -0.4028 0.0000\n",
" Vega f502n -0.2671 -0.2671 0.0000\n",
" Vega f550m 0.0724 0.0724 0.0000\n",
" Vega f555w -0.0456 -0.0456 0.0000\n",
" Vega f606w 0.2636 0.2636 0.0000\n",
" Vega f625w 0.4784 0.4784 0.0000\n",
" Vega f658n 0.7824 0.7824 0.0000\n",
" Vega f660n 0.6913 0.6913 0.0000\n",
" Vega f814w 1.2636 1.2636 0.0000\n",
" Vega f850lp 1.6128 1.6128 0.0000\n",
" Vega f892n 1.5539 1.5539 0.0000\n",
"g191b2b f435w 10.8214 10.8214 0.0000\n",
"g191b2b f475w 11.1800 11.1800 0.0000\n",
"g191b2b f502n 11.4247 11.4247 0.0000\n",
"g191b2b f550m 11.8612 11.8612 0.0000\n",
"g191b2b f555w 11.6868 11.6868 0.0000\n",
"g191b2b f606w 12.0691 12.0691 0.0000\n",
"g191b2b f625w 12.3666 12.3666 0.0000\n",
"g191b2b f658n 12.6164 12.6164 0.0000\n",
"g191b2b f660n 12.6001 12.6001 0.0000\n",
"g191b2b f814w 13.3710 13.3710 0.0000\n",
"g191b2b f850lp 13.8700 13.8700 0.0000\n",
"g191b2b f892n 13.8265 13.8265 0.0000\n"
]
}
],
"source": [
"# STMAG\n",
"t['source', 'obsmode', 'st_pysyn', 'st_st', 'd_st'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode pnu_pysyn pnu_st d_pnu \n",
"------- ------- -------------- -------------- --------\n",
" Vega f435w 7.730035e+28 8.847673e-09 100.0000\n",
" Vega f475w 8.389408e+28 9.444223e-09 100.0000\n",
" Vega f502n 8.880312e+28 9.881408e-09 100.0000\n",
" Vega f550m 8.90155e+28 9.907053e-09 100.0000\n",
" Vega f555w 8.744736e+28 9.7731e-09 100.0000\n",
" Vega f606w 8.752616e+28 9.762139e-09 100.0000\n",
" Vega f625w 8.810981e+28 9.786893e-09 100.0000\n",
" Vega f658n 7.606809e+28 8.4731875e-09 100.0000\n",
" Vega f660n 8.331032e+28 9.272819e-09 100.0000\n",
" Vega f814w 8.820841e+28 9.851898e-09 100.0000\n",
" Vega f850lp 9.096622e+28 1.01971365e-08 100.0000\n",
" Vega f892n 9.278707e+28 1.0332995e-08 100.0000\n",
"g191b2b f435w 2.0714068e+24 2.286499e-13 100.0000\n",
"g191b2b f475w 1.9525216e+24 2.1372029e-13 100.0000\n",
"g191b2b f502n 1.8695163e+24 2.0799981e-13 100.0000\n",
"g191b2b f550m 1.7137296e+24 1.902663e-13 100.0000\n",
"g191b2b f555w 1.7732441e+24 1.953939e-13 100.0000\n",
"g191b2b f606w 1.659247e+24 1.7873878e-13 100.0000\n",
"g191b2b f625w 1.5478852e+24 1.7020971e-13 100.0000\n",
"g191b2b f658n 1.4047474e+24 1.5633577e-13 100.0000\n",
"g191b2b f660n 1.436171e+24 1.5980342e-13 100.0000\n",
"g191b2b f814w 1.2662566e+24 1.384726e-13 100.0000\n",
"g191b2b f850lp 1.1376806e+24 1.2543943e-13 100.0000\n",
"g191b2b f892n 1.14404725e+24 1.2727157e-13 100.0000\n"
]
}
],
"source": [
"# PHOTNU\n",
"t['source', 'obsmode', 'pnu_pysyn', 'pnu_st', 'd_pnu'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode fnu_pysyn fnu_st d_fnu \n",
"------- ------- ------------- ------------- ------\n",
" Vega f435w 3.9749725e-20 4.0147043e-20 0.9995\n",
" Vega f475w 3.9536202e-20 3.9227186e-20 0.7816\n",
" Vega f502n 3.9077815e-20 3.9076726e-20 0.0028\n",
" Vega f550m 3.5295506e-20 3.5241002e-20 0.1544\n",
" Vega f555w 3.630073e-20 3.6080113e-20 0.6077\n",
" Vega f606w 3.3317976e-20 3.251441e-20 2.4118\n",
" Vega f625w 3.1055585e-20 3.0750833e-20 0.9813\n",
" Vega f658n 2.5537552e-20 2.5554763e-20 0.0674\n",
" Vega f660n 2.7903187e-20 2.7908117e-20 0.0177\n",
" Vega f814w 2.4479938e-20 2.4208402e-20 1.1092\n",
" Vega f850lp 2.2373305e-20 2.2328678e-20 0.1995\n",
" Vega f892n 2.3006927e-20 2.3017541e-20 0.0461\n",
"g191b2b f435w 1.0651679e-24 1.0494221e-24 1.4782\n",
"g191b2b f475w 9.201518e-25 8.958352e-25 2.6427\n",
"g191b2b f502n 8.226807e-25 8.225842e-25 0.0117\n",
"g191b2b f550m 6.7951034e-25 6.773608e-25 0.3163\n",
"g191b2b f555w 7.3610066e-25 7.2473583e-25 1.5439\n",
"g191b2b f606w 6.3161406e-25 6.0209187e-25 4.6741\n",
"g191b2b f625w 5.4557467e-25 5.3659477e-25 1.6460\n",
"g191b2b f658n 4.716013e-25 4.7163844e-25 0.0079\n",
"g191b2b f660n 4.8101785e-25 4.8099975e-25 0.0038\n",
"g191b2b f814w 3.514164e-25 3.4276916e-25 2.4607\n",
"g191b2b f850lp 2.7981458e-25 2.7622793e-25 1.2818\n",
"g191b2b f892n 2.836711e-25 2.8359818e-25 0.0257\n"
]
}
],
"source": [
"# FNU\n",
"t['source', 'obsmode', 'fnu_pysyn', 'fnu_st', 'd_fnu'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode jy_pysyn jy_st d_jy \n",
"------- ------- ----------- ----------- ------\n",
" Vega f435w 3974.9727 4014.7043 0.9995\n",
" Vega f475w 3953.6204 3922.7188 0.7816\n",
" Vega f502n 3907.7817 3907.6726 0.0028\n",
" Vega f550m 3529.5508 3524.1003 0.1544\n",
" Vega f555w 3630.0732 3608.0115 0.6077\n",
" Vega f606w 3331.7976 3251.441 2.4118\n",
" Vega f625w 3105.5586 3075.0835 0.9813\n",
" Vega f658n 2553.7551 2555.4763 0.0674\n",
" Vega f660n 2790.3186 2790.8118 0.0177\n",
" Vega f814w 2447.994 2420.84 1.1092\n",
" Vega f850lp 2237.3303 2232.8677 0.1995\n",
" Vega f892n 2300.6926 2301.7542 0.0461\n",
"g191b2b f435w 0.10651679 0.10494221 1.4782\n",
"g191b2b f475w 0.092015184 0.089583516 2.6427\n",
"g191b2b f502n 0.08226807 0.082258426 0.0117\n",
"g191b2b f550m 0.06795103 0.06773608 0.3163\n",
"g191b2b f555w 0.07361007 0.072473586 1.5439\n",
"g191b2b f606w 0.0631614 0.060209185 4.6741\n",
"g191b2b f625w 0.054557465 0.053659476 1.6460\n",
"g191b2b f658n 0.04716013 0.04716384 0.0079\n",
"g191b2b f660n 0.048101783 0.048099972 0.0038\n",
"g191b2b f814w 0.03514164 0.034276918 2.4607\n",
"g191b2b f850lp 0.027981458 0.027622793 1.2818\n",
"g191b2b f892n 0.028367111 0.02835982 0.0257\n"
]
}
],
"source": [
"# Jy\n",
"t['source', 'obsmode', 'jy_pysyn', 'jy_st', 'd_jy'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode njy_pysyn njy_st d_njy \n",
"------- ------- --------------- --------------- ------\n",
" Vega f435w 3974972600000.0 4014704200000.0 0.9995\n",
" Vega f475w 3953620400000.0 3922718600000.0 0.7816\n",
" Vega f502n 3907781700000.0 3907672600000.0 0.0028\n",
" Vega f550m 3529550600000.0 3524100400000.0 0.1544\n",
" Vega f555w 3630073000000.0 3608011300000.0 0.6077\n",
" Vega f606w 3331797500000.0 3251441000000.0 2.4118\n",
" Vega f625w 3105558600000.0 3075083300000.0 0.9813\n",
" Vega f658n 2553755200000.0 2555476400000.0 0.0674\n",
" Vega f660n 2790318700000.0 2790811800000.0 0.0177\n",
" Vega f814w 2447993700000.0 2420840000000.0 1.1092\n",
" Vega f850lp 2237330400000.0 2232868000000.0 0.1995\n",
" Vega f892n 2300692700000.0 2301754000000.0 0.0461\n",
"g191b2b f435w 106516790.0 104942210.0 1.4782\n",
"g191b2b f475w 92015180.0 89583510.0 2.6427\n",
"g191b2b f502n 82268070.0 82258424.0 0.0117\n",
"g191b2b f550m 67951030.0 67736080.0 0.3163\n",
"g191b2b f555w 73610060.0 72473580.0 1.5439\n",
"g191b2b f606w 63161404.0 60209184.0 4.6741\n",
"g191b2b f625w 54557468.0 53659476.0 1.6460\n",
"g191b2b f658n 47160132.0 47163840.0 0.0079\n",
"g191b2b f660n 48101784.0 48099972.0 0.0038\n",
"g191b2b f814w 35141640.0 34276916.0 2.4607\n",
"g191b2b f850lp 27981458.0 27622794.0 1.2818\n",
"g191b2b f892n 28367112.0 28359820.0 0.0257\n"
]
}
],
"source": [
"# nJy\n",
"t['source', 'obsmode', 'njy_pysyn', 'njy_st', 'd_njy'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode ab_pysyn ab_st d_ab \n",
"------- ------- -------- ------- ------\n",
" Vega f435w -0.0983 -0.1091 0.0108\n",
" Vega f475w -0.0925 -0.0840 0.0085\n",
" Vega f502n -0.0798 -0.0798 0.0000\n",
" Vega f550m 0.0307 0.0324 0.0017\n",
" Vega f555w 0.0002 0.0068 0.0066\n",
" Vega f606w 0.0933 0.1198 0.0265\n",
" Vega f625w 0.1697 0.1804 0.0107\n",
" Vega f658n 0.3821 0.3813 0.0007\n",
" Vega f660n 0.2859 0.2857 0.0002\n",
" Vega f814w 0.4280 0.4401 0.0121\n",
" Vega f850lp 0.5257 0.5278 0.0022\n",
" Vega f892n 0.4954 0.4949 0.0005\n",
"g191b2b f435w 11.3315 11.3476 0.0162\n",
"g191b2b f475w 11.4904 11.5194 0.0291\n",
"g191b2b f502n 11.6119 11.6120 0.0001\n",
"g191b2b f550m 11.8195 11.8229 0.0034\n",
"g191b2b f555w 11.7327 11.7496 0.0169\n",
"g191b2b f606w 11.8989 11.9508 0.0520\n",
"g191b2b f625w 12.0579 12.0759 0.0180\n",
"g191b2b f658n 12.2161 12.2160 0.0001\n",
"g191b2b f660n 12.1946 12.1946 0.0000\n",
"g191b2b f814w 12.5354 12.5625 0.0271\n",
"g191b2b f850lp 12.7828 12.7968 0.0140\n",
"g191b2b f892n 12.7680 12.7682 0.0003\n"
]
}
],
"source": [
"# ABMAG\n",
"t['source', 'obsmode', 'ab_pysyn', 'ab_st', 'd_ab'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode vmag_pysyn vmag_st d_vmag\n",
"------- ------- ---------- ------- ------\n",
" Vega f435w 0.0000 0.0000 0.0000\n",
" Vega f475w 0.0000 0.0000 0.0000\n",
" Vega f502n 0.0000 0.0000 0.0000\n",
" Vega f550m 0.0000 0.0000 0.0000\n",
" Vega f555w 0.0000 -0.0000 0.0000\n",
" Vega f606w 0.0000 0.0000 0.0000\n",
" Vega f625w 0.0000 0.0000 0.0000\n",
" Vega f658n 0.0000 0.0000 0.0000\n",
" Vega f660n 0.0000 -0.0000 0.0000\n",
" Vega f814w 0.0000 0.0000 0.0000\n",
" Vega f850lp 0.0000 -0.0000 0.0000\n",
" Vega f892n 0.0000 0.0000 0.0000\n",
"g191b2b f435w 11.4298 11.3855 0.0443\n",
"g191b2b f475w 11.5828 11.5924 0.0096\n",
"g191b2b f502n 11.6917 11.6918 0.0000\n",
"g191b2b f550m 11.7888 11.7913 0.0025\n",
"g191b2b f555w 11.7324 11.7419 0.0095\n",
"g191b2b f606w 11.8056 11.8391 0.0335\n",
"g191b2b f625w 11.8882 11.8975 0.0093\n",
"g191b2b f658n 11.8340 11.8167 0.0174\n",
"g191b2b f660n 11.9087 11.9058 0.0029\n",
"g191b2b f814w 12.1075 12.1256 0.0182\n",
"g191b2b f850lp 12.2571 12.2638 0.0066\n",
"g191b2b f892n 12.2726 12.2619 0.0107\n"
]
}
],
"source": [
"# VEGAMAG\n",
"t['source', 'obsmode', 'vmag_pysyn', 'vmag_st', 'd_vmag'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode ct_pysyn ct_st d_ct \n",
"------- ------- ------------- ------------- ------\n",
" Vega f435w 20189444000.0 20189840000.0 0.0020\n",
" Vega f475w 28787247000.0 28787812000.0 0.0020\n",
" Vega f502n 882509500.0 882526850.0 0.0020\n",
" Vega f550m 8508281300.0 8508448300.0 0.0020\n",
" Vega f555w 19277093000.0 19277472000.0 0.0020\n",
" Vega f606w 36447630000.0 36448340000.0 0.0020\n",
" Vega f625w 19664169000.0 19664554000.0 0.0020\n",
" Vega f658n 895943400.0 895961000.0 0.0020\n",
" Vega f660n 372359600.0 372366900.0 0.0020\n",
" Vega f814w 16085079000.0 16085394000.0 0.0020\n",
" Vega f850lp 5404658000.0 5404764000.0 0.0020\n",
" Vega f892n 577606200.0 577617540.0 0.0020\n",
"g191b2b f435w 541013.75 541024.3 0.0020\n",
"g191b2b f475w 669984.4 669997.5 0.0020\n",
"g191b2b f502n 18578.918 18579.281 0.0020\n",
"g191b2b f550m 163801.73 163804.94 0.0020\n",
"g191b2b f555w 390897.94 390905.62 0.0020\n",
"g191b2b f606w 690943.4 690956.94 0.0020\n",
"g191b2b f625w 345453.84 345460.62 0.0020\n",
"g191b2b f658n 16545.363 16545.69 0.0020\n",
"g191b2b f660n 6419.038 6419.164 0.0020\n",
"g191b2b f814w 230905.83 230910.36 0.0020\n",
"g191b2b f850lp 67594.05 67595.375 0.0020\n",
"g191b2b f892n 7121.777 7121.9165 0.0020\n"
]
}
],
"source": [
"# Count rate (from effstim)\n",
"t['source', 'obsmode', 'ct_pysyn', 'ct_st', 'd_ct'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode cr_pysyn cr_st d_cr \n",
"------- ------- ------------- ------------- ------\n",
" Vega f435w 20189434000.0 20189830000.0 0.0020\n",
" Vega f475w 28787243000.0 28787806000.0 0.0020\n",
" Vega f502n 882509100.0 882526460.0 0.0020\n",
" Vega f550m 8508282000.0 8508449000.0 0.0020\n",
" Vega f555w 19277095000.0 19277474000.0 0.0020\n",
" Vega f606w 36447633000.0 36448346000.0 0.0020\n",
" Vega f625w 19664169000.0 19664554000.0 0.0020\n",
" Vega f658n 895961150.0 895978700.0 0.0020\n",
" Vega f660n 372347200.0 372354500.0 0.0020\n",
" Vega f814w 16085079000.0 16085394000.0 0.0020\n",
" Vega f850lp 5404658700.0 5404764700.0 0.0020\n",
" Vega f892n 577605800.0 577617150.0 0.0020\n",
"g191b2b f435w 541013.6 541024.2 0.0020\n",
"g191b2b f475w 669984.3 669997.5 0.0020\n",
"g191b2b f502n 18578.914 18579.28 0.0020\n",
"g191b2b f550m 163801.73 163804.95 0.0020\n",
"g191b2b f555w 390897.97 390905.66 0.0020\n",
"g191b2b f606w 690943.44 690956.94 0.0020\n",
"g191b2b f625w 345453.88 345460.66 0.0020\n",
"g191b2b f658n 16545.441 16545.766 0.0020\n",
"g191b2b f660n 6418.988 6419.114 0.0020\n",
"g191b2b f814w 230905.83 230910.36 0.0020\n",
"g191b2b f850lp 67594.055 67595.375 0.0020\n",
"g191b2b f892n 7121.7783 7121.918 0.0020\n"
]
}
],
"source": [
"# Count rate (from countrate)\n",
"t['source', 'obsmode', 'cr_pysyn', 'cr_st', 'd_cr'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
" source obsmode obmag_pysyn obmag_st d_ob \n",
"------- ------- ----------- -------- ------\n",
" Vega f435w -25.7628 -25.7628 0.0000\n",
" Vega f475w -26.1480 -26.1480 0.0000\n",
" Vega f502n -22.3643 -22.3643 0.0000\n",
" Vega f550m -24.8246 -24.8246 0.0000\n",
" Vega f555w -25.7126 -25.7126 0.0000\n",
" Vega f606w -26.4042 -26.4042 0.0000\n",
" Vega f625w -25.7342 -25.7342 0.0000\n",
" Vega f658n -22.3807 -22.3807 0.0000\n",
" Vega f660n -21.4274 -21.4274 0.0000\n",
" Vega f814w -25.5161 -25.5161 0.0000\n",
" Vega f850lp -24.3319 -24.3319 0.0000\n",
" Vega f892n -21.9041 -21.9041 0.0000\n",
"g191b2b f435w -14.3330 -14.3330 0.0000\n",
"g191b2b f475w -14.5652 -14.5652 0.0000\n",
"g191b2b f502n -10.6726 -10.6726 0.0000\n",
"g191b2b f550m -13.0358 -13.0358 0.0000\n",
"g191b2b f555w -13.9802 -13.9802 0.0000\n",
"g191b2b f606w -14.5986 -14.5986 0.0000\n",
"g191b2b f625w -13.8460 -13.8460 0.0000\n",
"g191b2b f658n -10.5467 -10.5467 0.0000\n",
"g191b2b f660n -9.5187 -9.5187 0.0000\n",
"g191b2b f814w -13.4086 -13.4086 0.0000\n",
"g191b2b f850lp -12.0748 -12.0748 0.0000\n",
"g191b2b f892n -9.6315 -9.6315 0.0000\n"
]
}
],
"source": [
"# OBMAG\n",
"t['source', 'obsmode', 'obmag_pysyn', 'obmag_st', 'd_ob'].pprint(max_lines=-1, max_width=-1)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.4"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment