{ "cells": [ { "cell_type": "markdown", "metadata": { "collapsed": true }, "source": [ "# XMM-13hr master catalogue\n", "## Checks and diagnostics\n", "There is very little data on XMM-13hr so there is very little data to compare to." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "This notebook was run with herschelhelp_internal version: \n", "0246c5d (Thu Jan 25 17:01:47 2018 +0000) [with local modifications]\n", "This notebook was executed on: \n", "2018-02-21 14:21:47.320364\n" ] } ], "source": [ "from herschelhelp_internal import git_version\n", "print(\"This notebook was run with herschelhelp_internal version: \\n{}\".format(git_version()))\n", "import datetime\n", "print(\"This notebook was executed on: \\n{}\".format(datetime.datetime.now()))" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "#%config InlineBackend.figure_format = 'svg'\n", "\n", "import matplotlib.pyplot as plt\n", "plt.rc('figure', figsize=(10, 6))\n", "plt.style.use('ggplot')\n", "\n", "import locale\n", "locale.setlocale(locale.LC_ALL, 'en_GB')\n", "\n", "import os\n", "import time\n", "import itertools\n", "\n", "from astropy.coordinates import SkyCoord\n", "from astropy.table import Table\n", "from astropy import units as u\n", "from astropy import visualization as vis\n", "import numpy as np\n", "from matplotlib_venn import venn3, venn2\n", "\n", "from herschelhelp_internal.masterlist import (nb_compare_mags, nb_ccplots, nb_histograms, find_last_ml_suffix,\n", " quick_checks)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Diagnostics done using: master_catalogue_xmm-13hr_20180221.fits\n" ] } ], "source": [ "OUT_DIR = os.environ.get('OUT_DIR', \"./data\")\n", "SUFFIX = find_last_ml_suffix()\n", "#SUFFIX = \"20171016\"\n", "\n", "master_catalogue_filename = \"master_catalogue_xmm-13hr_{}.fits\".format(SUFFIX)\n", "master_catalogue = Table.read(\"{}/{}\".format(OUT_DIR, master_catalogue_filename))\n", "\n", "print(\"Diagnostics done using: {}\".format(master_catalogue_filename))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 0 - Quick checks" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Table shows only problematic columns.\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/opt/anaconda3/envs/herschelhelp_internal/lib/python3.6/site-packages/numpy/core/numeric.py:301: FutureWarning: in the future, full(32, False) will return an array of dtype('bool')\n", " format(shape, fill_value, array(fill_value).dtype), FutureWarning)\n", "/opt/anaconda3/envs/herschelhelp_internal/lib/python3.6/site-packages/numpy/core/numeric.py:301: FutureWarning: in the future, full(32, 0) will return an array of dtype('int64')\n", " format(shape, fill_value, array(fill_value).dtype), FutureWarning)\n" ] }, { "data": { "text/html": [ "<Table length=9>\n", "
idx | Column | All nan | #Measurements | #Zeros | #Negative | Minimum value |
---|---|---|---|---|---|---|
0 | merr_ap_90prime_g | False | 37185 | 0 | 2055 | -2633.04052734 |
1 | f_90prime_g | False | 37693 | 463 | 1631 | -20246.65625 |
2 | merr_90prime_g | False | 37230 | 0 | 1631 | -17914.5019531 |
3 | merr_ap_90prime_r | False | 37674 | 0 | 662 | -2159.73193359 |
4 | f_90prime_r | False | 37693 | 0 | 502 | -319.387512207 |
5 | merr_90prime_r | False | 37693 | 0 | 502 | -1078.2043457 |
6 | merr_ap_mosaic_z | False | 37323 | 0 | 821 | -874.305114746 |
7 | f_mosaic_z | False | 37693 | 337 | 283 | -2101.95874023 |
8 | merr_mosaic_z | False | 37356 | 0 | 283 | -288.257385254 |