{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Final Processing of SA13 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": 135, "metadata": { "collapsed": true }, "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": 136, "metadata": { "collapsed": true }, "outputs": [], "source": [ "cat=Table.read('./data/dmu22_XID+SPIRE_SA13_BLIND.fits')\n", "cat['RA'].unit=u.deg\n", "cat['Dec'].unit=u.deg" ] }, { "cell_type": "code", "execution_count": 137, "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 |
8 | 197.855149993 | 42.5671160212 | 67.8592 | 70.8763 | 64.6702 | 61.9552 | 63.8792 | 59.3428 | 45.0168 | 49.4403 | 40.5577 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999346 | 0.999348 | 0.998874 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.0 | 0.0 |
16 | 197.893045015 | 42.4977312992 | 44.0255 | 48.5561 | 39.9203 | 31.6237 | 35.3951 | 27.8582 | 19.0988 | 23.9203 | 14.2767 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999268 | 0.998481 | 0.999528 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.001 | 0.0 |
23 | 197.952630645 | 42.4605772106 | 47.8785 | 52.8323 | 42.9732 | 36.9955 | 41.9094 | 32.1833 | 44.4419 | 51.9702 | 37.1096 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999853 | 0.999552 | 0.998209 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.001 | 0.0 |
35 | 197.864409222 | 42.49879897 | 29.5685 | 34.166 | 25.0656 | 27.9472 | 31.6208 | 24.1863 | 7.34832 | 11.0765 | 3.43822 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999546 | 1.00016 | 1.00098 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.0 | 0.0 |
66 | 197.897155093 | 42.5113478736 | 47.7785 | 52.0938 | 43.4461 | 70.8185 | 74.0807 | 67.4937 | 58.0083 | 62.633 | 53.7125 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.99885 | 1.00044 | 1.00039 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.001 | 0.0 |
101 | 197.930695344 | 42.5099996023 | 25.2612 | 28.7756 | 21.9671 | 21.5853 | 25.5749 | 17.6083 | 15.1688 | 19.4054 | 10.742 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 1.00196 | 1.0018 | 0.999545 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.0 | 0.0 |
105 | 197.931033478 | 42.5297219563 | 21.7857 | 25.296 | 18.4972 | 30.4202 | 33.3535 | 27.3395 | 24.3927 | 28.5369 | 20.1483 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999197 | 1.00171 | 0.999674 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.007 | 0.0 |
123 | 198.002314994 | 42.4806095401 | 17.5904 | 22.2433 | 12.9469 | 2.89647 | 6.45764 | 0.853462 | 5.1742 | 9.7608 | 1.88636 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 1.0002 | 0.999584 | 1.00069 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.0 | 0.0 |
202 | 197.943605894 | 42.4539024374 | 19.699 | 24.5401 | 14.5753 | 10.6107 | 16.202 | 5.30745 | 2.29366 | 4.45961 | 0.690519 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 0.999728 | 1.00261 | 0.999986 | 2000.0 | 2000.0 | 2000.0 | 0.0 | 0.0 | 0.0 |
206 | 197.964198464 | 42.5336412393 | 8.40086 | 11.5089 | 5.13179 | 3.92352 | 6.83111 | 1.4877 | 1.53595 | 3.50255 | 0.389238 | -0.0708282 | -0.0872433 | -0.05502 | 0.005978 | 0.00867867 | 0.0115186 | 1.00035 | 0.999029 | 0.999893 | 2000.0 | 2000.0 | 2000.0 | 0.001 | 0.0 | 0.0 |