{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Final Processing of HDF-N Blind source catalogue\n", "\n", "![HELP LOGO](https://avatars1.githubusercontent.com/u/7880370?s=100&v=4>)\n", "\n", "\n", "The final processing stage requires:\n", "1. Quick validation of blind catalogues and Bayesian Pvalue maps\n", "2. Skewness level\n", "3. Adding flag to catalogue\n", "4. Merging MF catalogue with XID+ flux densities" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import seaborn as sns\n", "from astropy.table import Table,hstack\n", "%matplotlib inline\n", "import numpy as np\n", "import pylab as plt\n", "\n", "from astropy import units as u\n", "from astropy.table import Column\n", "\n", "import herschelhelp_internal\n", "from herschelhelp_internal.utils import gen_help_id\n", "import numpy.core.defchararray as np_f\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read tables" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "cat=Table.read('./data/dmu22_XID+SPIRE_HDF-N_BLIND.fits')\n", "cat['RA'].unit=u.deg\n", "cat['Dec'].unit=u.deg" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "Table length=10\n", "
HELP_ID | RA | Dec | F_SPIRE_250 | FErr_SPIRE_250_u | FErr_SPIRE_250_l | F_SPIRE_350 | FErr_SPIRE_350_u | FErr_SPIRE_350_l | F_SPIRE_500 | FErr_SPIRE_500_u | FErr_SPIRE_500_l | Bkg_SPIRE_250 | Bkg_SPIRE_350 | Bkg_SPIRE_500 | Sig_conf_SPIRE_250 | Sig_conf_SPIRE_350 | Sig_conf_SPIRE_500 | Rhat_SPIRE_250 | Rhat_SPIRE_350 | Rhat_SPIRE_500 | n_eff_SPIRE_250 | n_eff_SPIRE_500 | n_eff_SPIRE_350 | Pval_res_250 | Pval_res_350 | Pval_res_500 |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
deg | deg | mJy | mJy | mJy | mJy | mJy | mJy | mJy | mJy | mJy | mJy/Beam | mJy/Beam | mJy/Beam | mJy/Beam | mJy/Beam | mJy/Beam | ||||||||||
bytes27 | float64 | float64 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 |
75 | 189.47207051330057 | 61.83578062744778 | 34.711334 | 36.23167 | 33.190384 | 19.386705 | 21.187244 | 17.51052 | 3.3328526 | 5.1653214 | 1.470307 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.99905246 | 0.9989022 | 1.0012667 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
112 | 189.5238629258716 | 61.836233003789474 | 25.403498 | 26.987457 | 23.738743 | 21.146824 | 23.24672 | 19.17715 | 9.826296 | 11.910365 | 7.7862062 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.9989988 | 0.99837524 | 0.99857944 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
162 | 189.48615133968767 | 61.83075473094658 | 43.63474 | 45.232418 | 41.948498 | 22.383675 | 24.397493 | 20.364653 | 8.66891 | 11.0693035 | 6.2495065 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.9990405 | 0.9983873 | 0.9987331 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 0.999 | 0.998 |
239 | 189.56380822099334 | 61.882527454707855 | 27.448551 | 28.943016 | 25.994795 | 12.093296 | 13.931255 | 10.147441 | 0.49339843 | 1.2023259 | 0.13614236 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.999261 | 0.9984511 | 0.9984796 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
270 | 189.48726585833506 | 61.82353070859996 | 9.781449 | 11.297208 | 8.303107 | 0.79481995 | 1.827133 | 0.21942252 | 0.44249716 | 1.1876771 | 0.11953026 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.99903667 | 0.99965024 | 0.9985619 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
311 | 189.50077124630025 | 61.88572402559173 | 26.289492 | 27.816673 | 24.73249 | 22.274868 | 24.099386 | 20.348732 | 17.361462 | 19.541756 | 15.092925 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.99901617 | 0.9986254 | 0.9992568 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
335 | 189.5137444621643 | 61.88625288078024 | 38.04361 | 39.532288 | 36.50855 | 28.59526 | 30.520775 | 26.655268 | 12.135117 | 14.264226 | 9.908193 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 1.0000267 | 0.9999896 | 0.9993363 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |
348 | 189.4806409579133 | 61.87520780875087 | 24.721268 | 26.571983 | 22.983389 | 23.526894 | 26.054724 | 21.080145 | 10.578096 | 13.497189 | 7.3615203 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.99937123 | 1.0007389 | 0.9985527 | 2000.0 | 2000.0 | 2000.0 | 0.999 | 1.0 | 1.0 |
356 | 189.49149797544604 | 61.836577514074236 | 14.518369 | 16.236134 | 12.866267 | 11.295403 | 13.394885 | 9.354682 | 8.656886 | 10.960348 | 6.47291 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.9982637 | 0.99849766 | 0.9993845 | 2000.0 | 2000.0 | 2000.0 | 0.998 | 1.0 | 0.971 |
868 | 189.46331500943853 | 61.8449624106938 | 7.858822 | 9.55664 | 6.318145 | 0.74476194 | 1.6726309 | 0.22208466 | 0.24079666 | 0.5806592 | 0.07355579 | -0.97889066 | -1.285639 | -1.6167406 | 3.091677 | 3.7658827 | 3.5072732 | 0.99962616 | 0.9991022 | 0.99939114 | 2000.0 | 2000.0 | 2000.0 | 1.0 | 1.0 | 1.0 |