{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "WARNING: AstropyDeprecationWarning: block_reduce was moved to the astropy.nddata.blocks module. Please update your import statement. [astropy.nddata.utils]\n" ] } ], "source": [ "import pylab\n", "import pymoc\n", "import xidplus\n", "import numpy as np\n", "%matplotlib inline\n", "from astropy.table import Table" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import seaborn as sns" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook uses all the raw data from the XID+MIPS catalogue, maps, PSF and relevant MOCs to create XID+ prior object and relevant tiling scheme" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in MOCs\n", "The selection functions required are the main MOC associated with the masterlist. As the prior for XID+ is based on IRAC detected sources coming from two different surveys at different depths (SERVS and SWIRE) I will split the XID+ run into two different runs. Here we use the SERVS depth." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "Sel_func=pymoc.MOC()\n", "# Sel_func.read('../data/CDFS-SWIRE/holes_CDFS-SWIRE_irac1_O16_MOC.fits')\n", "Sel_func.read('../../dmu4/dmu4_sm_CDFS-SWIRE/data/holes_CDFS-SWIRE_irac1_O16_MOC.fits') # Update masterlist" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in XID+MIPS catalogue" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# XID_MIPS=Table.read('../data/CDFS-SWIRE/MIPS/dmu26_XID+MIPS_CDFS-SWIRE_cat_20170901.fits')\n", "XID_MIPS=Table.read('../dmu26_XID+MIPS_CDFS-SWIRE/data/dmu26_XID+MIPS_CDFS-SWIRE_cat_20201205.fits') # Update masterlist" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/html": [ "Table length=10\n", "
help_id | RA | Dec | F_MIPS_24 | FErr_MIPS_24_u | FErr_MIPS_24_l | Bkg_MIPS_24 | Sig_conf_MIPS_24 | Rhat_MIPS_24 | n_eff_MIPS_24 | Pval_res_24 | flag_mips_24 |
---|---|---|---|---|---|---|---|---|---|---|---|
degrees | degrees | muJy | muJy | muJy | MJy / sr | MJy / sr | |||||
bytes27 | float64 | float64 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | float32 | bool |
HELP_J033650.824-295620.783 | 54.211766390022845 | -29.9391064859918 | 272.50247 | 540.70447 | 86.236374 | -0.003483794 | 5.029496e-06 | 1.0027044 | 1124.0 | 0.0 | False |
HELP_J033655.442-295449.596 | 54.23101021455793 | -29.91377659143576 | 236.55391 | 506.73215 | 68.16577 | -0.0006361296 | 5.0581343e-06 | 0.9994119 | 1476.0 | 0.0 | False |
HELP_J033616.314-295548.181 | 54.067976200022855 | -29.9300503059918 | 301.9198 | 318.07288 | 286.00903 | -0.004160782 | 5.0500876e-06 | nan | 1700.0 | 0.0 | False |
HELP_J033601.094-295546.425 | 54.00455690002285 | -29.929562435991798 | 81.58568 | 95.428375 | 67.528366 | 0.0016651336 | 4.938521e-06 | nan | 1122.0 | 0.0 | False |
HELP_J033608.786-295409.821 | 54.03660934002285 | -29.9027279559918 | 6.9124126 | 15.632443 | 2.0185032 | -0.0052455342 | 5.063621e-06 | nan | 2000.0 | 0.0 | True |
HELP_J033606.888-295518.496 | 54.028701300022846 | -29.9218044959918 | 36.028355 | 50.140163 | 21.682087 | -0.0052455342 | 5.063621e-06 | nan | 719.0 | 0.0 | False |
HELP_J033610.856-295506.940 | 54.04523478002285 | -29.918594405991797 | 250.45511 | 263.87708 | 237.46614 | -0.0052455342 | 5.063621e-06 | nan | 1863.0 | 0.0 | False |
HELP_J033607.350-295516.790 | 54.030623590022856 | -29.9213304359918 | 184.89633 | 199.12378 | 170.45247 | -0.0052455342 | 5.063621e-06 | nan | 1419.0 | 0.0 | False |
HELP_J033555.160-295410.647 | 53.97983175002285 | -29.9029574059918 | 1659.3827 | 1672.8367 | 1645.4224 | -0.0065137288 | 5.1315847e-06 | nan | 2000.0 | 0.993 | True |
HELP_J033557.604-295401.282 | 53.99001862002285 | -29.900356185991797 | 15.390374 | 26.97467 | 5.390874 | -0.0065137288 | 5.1315847e-06 | 0.9992996 | 1482.0 | 0.0 | True |