{ "cells": [ { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/mc741/anaconda3/lib/python3.6/site-packages/mpl_toolkits/axes_grid/__init__.py:12: MatplotlibDeprecationWarning: \n", "The mpl_toolkits.axes_grid module was deprecated in Matplotlib 2.1 and will be removed two minor releases later. Use mpl_toolkits.axes_grid1 and mpl_toolkits.axisartist provies the same functionality instead.\n", " obj_type='module')\n" ] } ], "source": [ "import pylab\n", "import pymoc\n", "import xidplus\n", "import numpy as np\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook uses all the raw data from the masterlist, 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. Here we use the SERVS depth." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Sel_func=pymoc.MOC()\n", "Sel_func.read('../../dmu4/dmu4_sm_ELAIS-N2/data/holes_ELAIS-N2_irac_i1_O16_20180921_MOC.fits')\n", "SWIRE_MOC=pymoc.MOC()\n", "SWIRE_MOC.read('../../dmu0/dmu0_DataFusion-Spitzer/data/DF-SWIRE_ELAIS-N2_MOC.fits')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Final=Sel_func.intersection(SWIRE_MOC)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Final.write('./data/testMoc.fits', overwrite=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in Masterlist\n", "Next step is to read in Masterlist and select only sources that are detected in mid-infrared and at least one other wavelength domain (i.e. optical or nir). This will remove most of the objects in the catalogue that are artefacts. We can do this by using the `flag_optnir_det` flag and selecting sources that have a binary value of $>= 5$" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from astropy.io import fits\n", "masterfile='master_catalogue_elais-n2_20180218.fits'\n", "masterlist=fits.open('../../dmu1/dmu1_ml_ELAIS-N2/data/'+masterfile)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "good=masterlist[1].data['flag_optnir_det']>=5" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create uninformative (i.e. conservative) upper and lower limits based on IRAC fluxes\n", "As the default flux prior for XID+ is a uniform distribution, it makes sense to set reasonable upper and lower 24 micron flux limits based on the longest wavelength IRAC flux available. For a lower limit I take IRAC/500.0 and for upper limit I take IRACx500." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": true }, "outputs": [], "source": [ "MIPS_lower=np.full(good.sum(),0.0)\n", "MIPS_upper=np.full(good.sum(),1E5)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "#masterlist[1].header" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": true }, "outputs": [], "source": [ "#for i in range(0,good.sum()):\n", "# if masterlist[1].data['flag_irac4'][good][i]>0:\n", "# MIPS_lower[i]=masterlist[1].data['f_irac4'][good][i]/500.0\n", "# MIPS_upper[i]=masterlist[1].data['f_irac4'][good][i]*500.0\n", "# elif masterlist[1].data['flag_irac3'][good][i]>0:\n", "# MIPS_lower[i]=masterlist[1].data['f_irac3'][good][i]/500.0\n", "# MIPS_upper[i]=masterlist[1].data['f_irac3'][good][i]*500.0\n", "# elif masterlist[1].data['flag_irac2'][good][i]>0:\n", "# MIPS_lower[i]=masterlist[1].data['f_irac2'][good][i]/500.0\n", "# MIPS_upper[i]=masterlist[1].data['f_irac2'][good][i]*500.0\n", "# elif masterlist[1].data['flag_irac1'][good][i]>0:\n", "# MIPS_lower[i]=masterlist[1].data['f_irac1'][good][i]/500.0\n", "# MIPS_upper[i]=masterlist[1].data['f_irac1'][good][i]*500.0" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": true }, "outputs": [], "source": [ "for i in range(0,good.sum()):\n", " if ~np.isnan(masterlist[1].data['f_irac_i4'][good][i]):\n", " MIPS_lower[i]=masterlist[1].data['f_irac_i4'][good][i]/500.0\n", " MIPS_upper[i]=masterlist[1].data['f_irac_i4'][good][i]*500.0\n", " elif ~np.isnan(masterlist[1].data['f_irac_i3'][good][i]):\n", " MIPS_lower[i]=masterlist[1].data['f_irac_i3'][good][i]/500.0\n", " MIPS_upper[i]=masterlist[1].data['f_irac_i3'][good][i]*500.0\n", " elif ~np.isnan(masterlist[1].data['f_irac_i2'][good][i]):\n", " MIPS_lower[i]=masterlist[1].data['f_irac_i2'][good][i]/500.0\n", " MIPS_upper[i]=masterlist[1].data['f_irac_i2'][good][i]*500.0\n", " elif ~np.isnan(masterlist[1].data['f_irac_i1'][good][i]):\n", " MIPS_lower[i]=masterlist[1].data['f_irac_i1'][good][i]/500.0\n", " MIPS_upper[i]=masterlist[1].data['f_irac_i1'][good][i]*500.0" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in Map\n", "We are now ready to read in the MIPS map\n" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [ { "ename": "NameError", "evalue": "name 'fits' is not defined", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mMIPS_Map\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mfits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'../../dmu17/dmu17_HELP_Legacy_Maps/ELAIS-N2/data/input_data/swire_EN2_M1_v4_help.fits'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'fits' is not defined" ] } ], "source": [ "MIPS_Map=fits.open('../../dmu17/dmu17_HELP_Legacy_Maps/ELAIS-N2/data/input_data/swire_EN2_M1_v4_help.fits')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAEJFJREFUeJzt3W+MnWlZx/Hvj+JgAroIuxrSP7ba2lB9AeykEFSyGv5MgVokRFtJBNO0WUL9E31BCSZIfLO+UCNSWYuULsS0aXCVjjtaDZEUk6qdJQgtTWGsazq22c5askg0lsLlizmVwzgzPWfOmZ7du99PMpl5rjnP/VznyZkr91zPc+6TqkKS1K7njDoBSdLqstBLUuMs9JLUOAu9JDXOQi9JjbPQS1LjLPSS1DgLvSQ17rnDHjDJBuBDwFPAl6vqoWEfQ5LUu54KfZIjwJuBa1X1Y13xCeAPgDXAn3SK+o8Aj1XVHyf5eC/j33vvvbVx48Z+c5eku9rjjz/+VFXdd7vHpZclEJK8Bvg68PFbhT7JGuDLwOuAWeAssAd4EvgkUMAnqupjtxt/fHy8pqenb5uHJOnbkjxeVeO3e1xPPfqqOg1cXxDeDsxU1aWqugEcB3YBvwS8v6p+GnhTf2lLkoZtkIuxa4HLXduzndhfA7+S5GHgiaV2TrI/yXSS6bm5uQHSkCQtZ5CLsVkkVlV1Dnjb7XauqsNJrgI7x8bG7h8gD0nSMgaZ0c8C67u21wFXBktHkjRsgxT6s8CWJJuSjAG7gZP9DFBVk1W1/5577hkgDUnScnoq9EmOAWeArUlmk+ytqpvAAeAUcAE4UVXn+zl4kp1JDj/99NP95i1J6lFPt1euNm+vlKT+DfX2ytXijF6SVt/Ql0DoR1VNApPj4+P7VjrGxoOPDTGj/jzxkG8TkPTM54xekho30kLvXTeStPpcpliSGmfrRpIaZ+tGkhpn60aSGmfrRpIaZ+tGkhpn60aSGmehl6TGWeglqXFejJWkxnkxVpIaZ+tGkhpnoZekxg19PfokPwm8vTP2tqp69bCPIUnqXa+fGXskybUk5xbEJ5JcTDKT5CBAVX22qh4E/hJ4ZPgpS5L60Wvr5igw0R1IsgY4BOwAtgF7kmzresgvAMeGkKMkaQA9FfqqOg1cXxDeDsxU1aWqugEcB3YBJNkAPF1VXxtmspKk/g1yMXYtcLlre7YTA9gLfGy5nZPsTzKdZHpubm6ANCRJyxnkYmwWiRVAVb3/djtX1eEkV4GdY2Nj9w+QhyRpGYPM6GeB9V3b64Ar/QzgG6YkafUNUujPAluSbEoyBuwGTvYzgEsgSNLq6/X2ymPAGWBrktkke6vqJnAAOAVcAE5U1fnVS1WStBI99eiras8S8SlgaqUHr6pJYHJ8fHzfSseQJC3P1SslqXGuXilJjXNGL0mNc0YvSY1zmWJJapytG0lqnK0bSWqcrRtJapytG0lqnK0bSWqcrRtJapyFXpIaZ6GXpMZ5MVaSGufFWElqnK0bSWqchV6SGmehl6TG9fRRgv1I8hzgt4HvBaar6pFhH0OS1LtePxz8SJJrSc4tiE8kuZhkJsnBTngXsBb4BjA73HQlSf3qtXVzFJjoDiRZAxwCdgDbgD1JtgFbgTNV9evAu4aXqiRpJXoq9FV1Gri+ILwdmKmqS1V1AzjO/Gx+Fvhq5zHfXGrMJPuTTCeZnpub6z9zSVJPBrkYuxa43LU924k9CrwhyR8Cp5fauaoOAx8APjc2NjZAGpKk5QxyMTaLxKqq/gvY28sAVTUJTI6Pj+8bIA9J0jIGmdHPAuu7ttcBV/oZwCUQJGn1DVLozwJbkmxKMgbsBk4OJy1J0rD0envlMeAMsDXJbJK9VXUTOACcAi4AJ6rqfD8Hd60bSVp9PfXoq2rPEvEpYGqlB0+yE9i5efPmlQ4hSboNV6+UpMa5Hr0kNc4ZvSQ1ztUrJalxtm4kqXG2biSpcbZuJKlxtm4kqXG2biSpcbZuJKlxFnpJapyFXpIa58VYSWqcF2MlqXG2biSpcRZ6SWqchV6SGjf0Qp/kgSSfTfJwkgeGPb4kqT+9fmbskSTXkpxbEJ9IcjHJTJKDnXABXwe+G5gdbrqSpH71OqM/Ckx0B5KsAQ4BO4BtwJ4k24DPVtUO4D3AB4aXqiRpJXoq9FV1Gri+ILwdmKmqS1V1AzgO7Kqqb3V+/1XgeUPLVJK0Is8dYN+1wOWu7VnglUneCrwBeCHwoaV2TrIf2A+wYcOGAdKQJC1nkEKfRWJVVY8Cj95u56o6nOQqsHNsbOz+AfKQJC1jkLtuZoH1XdvrgCuDpSNJGrZBCv1ZYEuSTUnGgN3AyX4GcAkESVp9vd5eeQw4A2xNMptkb1XdBA4Ap4ALwImqOt/PwV3UTJJWX089+qras0R8Cpha6cGrahKYHB8f37fSMSRJy3OZYklqnMsUS1LjnNFLUuOc0UtS45zRS1LjnNFLUuP84BFJapyFXpIaZ49ekhpnj16SGmfrRpIaZ6GXpMZZ6CWpcV6MlaTGeTFWkhpn60aSGmehl6TGrUqhT/L8JI8nefNqjC9J6l2vnxl7JMm1JOcWxCeSXEwyk+Rg16/eA5wYZqKSpJXpdUZ/FJjoDiRZAxwCdgDbgD1JtiV5LfAl4Mkh5ilJWqFePxz8dJKNC8LbgZmqugSQ5DiwC3gB8Hzmi/9/J5mqqm8NLWNJUl96KvRLWAtc7tqeBV5ZVQcAkrwTeGqpIp9kP7AfYMOGDQOkIUlaziCFPovE6v9+qDq63M5VdTjJVWDn2NjY/QPkIUlaxiB33cwC67u21wFXBktHkjRsgxT6s8CWJJuSjAG7gZP9DOA7YyVp9fV6e+Ux4AywNclskr1VdRM4AJwCLgAnqup8Pwd3rRtJWn293nWzZ4n4FDC10oNX1SQwOT4+vm+lY0iSlufqlZLUOFevlKTGOaOXpMY5o5ekxrlMsSQ1ztaNJDXO1o0kNc7WjSQ1ztaNJDXO1o0kNc7WjSQ1zkIvSY2z0EtS47wYK0mN82KsJDXO1o0kNc5CL0mNG3qhT/LSJA8n+WSSdw17fElSf3r9zNgjSa4lObcgPpHkYpKZJAcBqupCVT0I/BwwPvyUJUn96HVGfxSY6A4kWQMcAnYA24A9SbZ1fvczwN8Dnx5appKkFemp0FfVaeD6gvB2YKaqLlXVDeA4sKvz+JNV9Wrg7cNMVpLUv+cOsO9a4HLX9izwyiQPAG8FngdMLbVzkv3AfoANGzYMkIYkaTmDFPosEquq+gzwmdvtXFWHk1wFdo6Njd0/QB6SpGUMctfNLLC+a3sdcKWfAXzDlCStvkEK/VlgS5JNScaA3cDJfgZwCQRJWn293l55DDgDbE0ym2RvVd0EDgCngAvAiao6v3qpSpJWoqcefVXtWSI+xTIXXHsYdxKYHB8f37fSMSRJy3P1SklqnKtXSlLjnNFLUuOc0UtS41ymWJIaZ+tGkhpn60aSGmfrRpIaZ+tGkhpn60aSGmfrRpIaZ6GXpMZZ6CWpcV6MlaTGeTFWkhpn60aSGmehl6TGWeglqXGrUuiTvCXJR5J8KsnrV+MYkqTe9FzokxxJci3JuQXxiSQXk8wkOQhQVX9RVfuAdwI/P9SMJUl96enDwTuOAh8CPn4rkGQNcAh4HTALnE1ysqq+1HnIb3Z+36SNBx8byXGfeOhNIzmupGennmf0VXUauL4gvB2YqapLVXUDOA7syrzfAf6qqj632HhJ9ieZTjI9Nze30vwlSbcxaI9+LXC5a3u2E/tl4LXA25I8uNiOVXW4qsaravy+++4bMA1J0lL6ad0sJovEqqo+CHzwtjsnO4GdmzdvHjANSdJSBp3RzwLru7bXAVcGHFOSNESDFvqzwJYkm5KMAbuBk73u7BIIkrT6+rm98hhwBtiaZDbJ3qq6CRwATgEXgBNVdb6PMV3UTJJWWc89+qras0R8CpgaWkaSpKFy9UpJapxr3UhS4/zgEUlqnK0bSWqcM3pJapwzeklqnBdjJalxFnpJapw9eklqnD16SWqcrRtJapyFXpIaZ49ekhpnj16SGmfrRpIaZ6GXpMZZ6CWpcT1/wlSvkvwQ8D7gnqp627DHF2w8+NjIjv3EQ28a2bElrUxPM/okR5JcS3JuQXwiycUkM0kOAlTVparauxrJSpL612vr5igw0R1IsgY4BOwAtgF7kmwbanaSpIH1VOir6jRwfUF4OzDTmcHfAI4Du4acnyRpQINcjF0LXO7angXWJnlxkoeBlyd571I7J9mfZDrJ9Nzc3ABpSJKWM8jF2CwSq6r6D+DB2+1cVYeTXAV2jo2N3T9AHpKkZQxS6GeB9V3b64Arg6UjLW5Udxp5l5FaMEjr5iywJcmmJGPAbuBkPwO4BIIkrb6eZvRJjgEPAPcmmQXeX1UfTXIAOAWsAY5U1fl+Dp5kJ7Bz8+bN/WWtkRnlPfySVqanQl9Ve5aITwFTKz14VU0Ck+Pj4/tWOoYkaXkuUyxJjXOZYklq3NDXuumHPXo907mukFrgjF6SGucyxZLUOFs30jOUbxK7c1pv0dm6kaTG2bqRpMbZupH0HXz3c3ts3UhS42zdSFLjLPSS1DgLvSQ1zkXNJKlxXoyVpMbZupGkxlnoJalxFnpJalyqatQ5kGQO+LcBh7kXeGoI6TzbeR48B7d4Hua1fB5+sKruu92DnhGFfhiSTFfV+KjzGDXPg+fgFs/DPM+DrRtJap6FXpIa11KhPzzqBJ4hPA+eg1s8D/Pu+vPQTI9ekrS4lmb0kqRFNFHok0wkuZhkJsnBUedzpyR5IskXk3w+yXQn9qIkf5vkK53v3zfqPIctyZEk15Kc64ot+rwz74Od18YXkrxidJkP1xLn4beS/HvnNfH5JG/s+t17O+fhYpI3jCbr4UqyPsnfJbmQ5HySX+3E77rXw3Ke9YU+yRrgELAD2AbsSbJttFndUT9VVS/run3sIPDpqtoCfLqz3ZqjwMSC2FLPewewpfO1H/jwHcrxTjjK/z8PAL/feU28rKqmADp/E7uBH+3s80edv51nu5vAb1TVS4FXAe/uPNe78fWwpGd9oQe2AzNVdamqbgDHgV0jzmmUdgGPdH5+BHjLCHNZFVV1Gri+ILzU894FfLzm/QPwwiQvuTOZrq4lzsNSdgHHq+p/qupfgRnm/3ae1arqalV9rvPzfwIXgLXcha+H5bRQ6NcCl7u2Zzuxu0EBf5Pk8ST7O7EfqKqrMP9HAHz/yLK7s5Z63nfj6+NApy1xpKt11/x5SLIReDnwj/h6+A4tFPosErtbbiX68ap6BfP/jr47yWtGndAz0N32+vgw8MPAy4CrwO924k2fhyQvAP4M+LWq+tpyD10k1sx5WEoLhX4WWN+1vQ64MqJc7qiqutL5fg34c+b/FX/y1r+ine/XRpfhHbXU876rXh9V9WRVfbOqvgV8hG+3Z5o9D0m+i/ki/6dV9Wgn7OuhSwuF/iywJcmmJGPMX3A6OeKcVl2S5yf5nls/A68HzjH/3N/Redg7gE+NJsM7bqnnfRL4xc7dFq8Cnr71L32LFvSbf5b51wTMn4fdSZ6XZBPzFyP/6U7nN2xJAnwUuFBVv9f1K18P3arqWf8FvBH4MvAvwPtGnc8des4/BPxz5+v8recNvJj5uwy+0vn+olHnugrP/RjzbYlvMD9D27vU82b+X/VDndfGF4HxUee/yufhE53n+QXmi9pLuh7/vs55uAjsGHX+QzoHP8F86+ULwOc7X2+8G18Py335zlhJalwLrRtJ0jIs9JLUOAu9JDXOQi9JjbPQS1LjLPSS1DgLvSQ1zkIvSY37X8LJAORrehJ1AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from matplotlib import pyplot as plt\n", "data=MIPS_Map[1].data\n", "plt.hist(data.flatten());\n", "plt.yscale('log')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "XTENSION= 'IMAGE ' / Image extension \n", "BITPIX = -32 / array data type \n", "NAXIS = 2 / number of array dimensions \n", "NAXIS1 = 6038 \n", "NAXIS2 = 9283 \n", "PCOUNT = 0 / number of parameters \n", "GCOUNT = 1 / number of groups \n", "TELESCOP= 'Spitzer ' / Spitzer Space Telescope \n", "INSTRUME= 'MIPS ' / Spitzer Space Telescope instrument ID \n", "CHNLNUM = 1 / This image: 1=24um,2=70um,3=160um \n", "EXPTYPE = 'scn ' / Exposure Type \n", "FOVID = 113 / Field of View ID for Commanded Pointing \n", "FOVNAME = 'MIPS_70um_scan' / Field of View Name for Commanded Pointing \n", "OBSRVR = 'Carol Lonsdale' / Observer Name (Last, First) \n", "OBSRVRID= 88 / Observer ID of Principal Investigator \n", "PROCYCL = 2 / Proposal Cycle \n", "PROGID = 183 / Program ID \n", "DATE_OBS= '2004-07-09T14:29:09.876' / Date & time at DCE start \n", "CRVAL1 = 249.1511577735449 / [deg] RA at CRPIX1,CRPIX2 \n", "CRVAL2 = 40.9865025 / [deg] DEC at CRPIX1,CRPIX2 \n", "RADESYS = 'ICRS ' / International Celestial Reference System \n", "EQUINOX = 2000. / Equinox for ICRS celestial coord. system \n", "CTYPE1 = 'RA---TAN' / RA projection type \n", "CTYPE2 = 'DEC--TAN' / DEC projection type \n", "CRPIX1 = 3019.628767 / Reference pixel along axis 1 \n", "CRPIX2 = 4642.180189 / Reference pixel along axis 2 \n", "CDELT1 = -0.000333333 / [deg/pix] Scale for axis 1 \n", "CDELT2 = 0.000333333 / [deg/pix] Scale for axis 2 \n", "CROTA2 = 124.96831 / [deg] Orientation of axis 2 (W of N, +=CW) \n", "BUNIT = 'MJy/sr ' / Units of image data \n", "AORKEY = 10055424 / AOR or IER key. Astro. Obs Req/Instr Eng Req \n", "DPID = 18543388 / Data Product Instance ID \n", "COMMENT MOSAIC_INT Module Version 3.1 image created \n", "COMMENT Mon Nov 22 17:51:58 2004 \n", "COMMENT MOSAIC_REINT Module Version 2.1 image created \n", "COMMENT Thu Nov 25 00:28:39 2004 \n", "COMMENT MOSAIC_COADD Module Version 3.3 image created \n", "COMMENT Thu Nov 25 07:06:17 2004 \n", "COMMENT MOSAIC_COMBINE Module Version 1.7 image created \n", "COMMENT Thu Nov 25 07:09:51 2004 \n", "TOTALBCD= 33695 / Total number of BCD's used in this mosaic \n", "EXTNAME = 'IMAGE ' \n", "CHECKSUM= '9DQjECPg9CPgCCPg' / HDU checksum updated 2018-10-16T15:29:52 \n", "DATASUM = '405582022' / data unit checksum updated 2018-10-16T15:29:52 " ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "MIPS_Map[1].header" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAERZJREFUeJzt3W+MnWlZx/Hvz66zJguCsKshbdcW2jTOK9mdLASV8IJAu1qKaLSFRDBNmyXWaIwJJZiI8YVioi/QCimhFg1p0+AqrTtkNURSTCruLOFPa1MZVsgObejgksV/cSlcvpizMhnnTM+/mdO5z/eTTGae+5znfq57n9mr91zPc+4nVYUkqV3fN+4AJEnry0QvSY0z0UtS40z0ktQ4E70kNc5EL0mNM9FLUuNM9JLUOBO9JDXurnEHAHDvvffWjh07xh2GJG0qTz755Deq6r7bvW+siT7JfmD/rl27mJubG2cokrTpJPlqL+8ba+mmqi5U1dEXvehF4wxDkppmjV6SGjfWRJ9kf5KTzz777DjDkKSmWbqRpMY5o5ekxjmjl6TGeTFWkhpn6UaSGjdRpZsdxx/bkONI0p3E0o0kNW7iEr2zekmTZuISvSRNGi/GSlLjJupirCRNIks3ktQ4E70kNW7kDx5J8lPA2zp9T1fVa0Z9DElS73qa0Sc5leRmkssr2vcmuZZkPslxgKr6dFU9AvwN8JHRhyxJ6kevpZvTwN7lDUm2ACeAfcA0cCjJ9LK3vBU4M4IYJUlD6CnRV9VF4JkVzQ8B81X1VFU9B5wFDgAkuR94tqq+NcpgJUn9G+Zi7Fbg6WXbC502gMPAn621c5KjSeaSzC0uLg4RhiRpLcNcjM0qbQVQVb99u52r6mSSG8D+qampB4eIQ5K0hmFm9AvA9mXb24Dr/XTgB6Ykaf0Nk+ifAHYn2ZlkCjgInO+nA5dAkKT11+vtlWeAS8CeJAtJDlfVLeAY8DhwFThXVVfWL1RJ0iB6qtFX1aEu7bPA7KAHr6oLwIWZmZkjg/YhSVqbq1dKUuNcvVKSGueMXpIa54xekhrnMsWS1DhLN5LUOEs3ktQ4SzeS1DhLN5LUOEs3ktQ4SzeS1DgTvSQ1zkQvSY3zYqwkNc6LsZLUOEs3ktQ4E70kNa6nJ0z1I8n3Ab8L/CAwV1UfGfUxJEm96/WZsaeS3ExyeUX73iTXkswnOd5pPgBsBb4NLIw2XElSv3ot3ZwG9i5vSLIFOAHsA6aBQ0mmgT3Apar6DeCdowtVkjSInhJ9VV0EnlnR/BAwX1VPVdVzwFmWZvMLwDc77/nOqAKVJA1mmIuxW4Gnl20vdNoeBd6Y5I+Bi912TnI0yVySucXFxSHCkCStZZiLsVmlrarqv4DDt9u5qk4muQHsn5qaenCIOCRJaxhmRr8AbF+2vQ24Plw4kqRRGybRPwHsTrIzyRRwEDjfTwd+MlaS1l+vt1eeAS4Be5IsJDlcVbeAY8DjwFXgXFVd6efgrnUjSeuvpxp9VR3q0j4LzA568Kq6AFyYmZk5MmgfkqS1uXqlJDXO1SslqXHO6CWpcc7oJalxzuglqXHO6CWpcT54RJIaZ+lGkhpn6UaSGmfpRpIaZ6KXpMaZ6CWpcV6MlaTGeTFWkhpn6UaSGmeil6TGjTzRJ3ldkk8n+WCS1426f0lSf3p9lOCpJDeTXF7RvjfJtSTzSY53mgv4D+AHWHqAuCRpjHqd0Z8G9i5vSLIFOAHsA6aBQ0mmgU9X1T7gXcDvjC5USdIgekr0VXUReGZF80PAfFU9VVXPAWeBA1X13c7r3wTuHlmkkqSB9PRw8C62Ak8v214AXpXkLcAbgRcDf9Jt5yRHgaMA999//xBhSJLWMkyizyptVVWPAo/ebueqOpnkBrB/amrqwSHikCStYZi7bhaA7cu2twHX++nAD0xJ0vobJtE/AexOsjPJFHAQON9PBy6BIEnrr9fbK88Al4A9SRaSHK6qW8Ax4HHgKnCuqq6sX6iSpEH0VKOvqkNd2meB2UEPXlUXgAszMzNHBu1DkrQ2V6+UpMa5eqUkNc4ZvSQ1zhm9JDXOZYolqXGWbiSpcZZuJKlxlm4kqXGWbiSpcZZuJKlxlm4kqXEmeklqnIlekhrnxVhJapwXYyWpcZZuJKlxE5nodxx/bNwhSNKGWZdEn+SeJE8m+Zn16F+S1Ltenxl7KsnNJJdXtO9Nci3JfJLjy156F3BulIFKkgbT64z+NLB3eUOSLcAJYB8wDRxKMp3k9cA/A18fYZySpAH1+nDwi0l2rGh+CJivqqcAkpwFDgAvAO5hKfn/d5LZqvruyCKWJPWlp0TfxVbg6WXbC8CrquoYQJJ3AN/oluSTHAWOAtx///1DhCFJWsswiT6rtNX//VB1eq2dq+pkkhvA/qmpqQeHiEOStIZh7rpZALYv294GXB8uHEnSqA2T6J8AdifZmWQKOAic76cDPxkrSeuv19srzwCXgD1JFpIcrqpbwDHgceAqcK6qrvRzcNe6kaT11+tdN4e6tM8Cs4MevKouABdmZmaODNqHJGltrl4pSY1z9UpJapwzeklqnDN6SWqcM3pJatzEzuh3HH/MdeklTYSJfPCIJE0SE70kNc4avSQ1bmJr9JI0KSzdSFLjTPSS1Dhr9JLUOGv0ktQ4SzeS1DgTvSQ1buITvcsgSGrdyBN9kh9L8sEkH0vyzlH3L0nqT6/PjD2V5GaSyyva9ya5lmQ+yXGAqrpaVY8AvwDMjD5kSVI/ep3Rnwb2Lm9IsgU4AewDpoFDSaY7r70J+AfgkyOLVJI0kJ4SfVVdBJ5Z0fwQMF9VT1XVc8BZ4EDn/eer6jXA27r1meRokrkkc4uLi4NFL0m6rbuG2Hcr8PSy7QXgVUleB7wFuBuY7bZzVZ1McgPYPzU19eAQcUiS1jBMos8qbVVVnwI+1UsHVXUBuDAzM3NkiDgkSWsY5q6bBWD7su1twPV+OnAJBElaf8Mk+ieA3Ul2JpkCDgLnRxPWxvJeekkt6/X2yjPAJWBPkoUkh6vqFnAMeBy4Cpyrqiv9HNy1biRp/fVUo6+qQ13aZ1njguvtJNkP7N+1a9egXUiSbsPVKyWpca5H32GdXlKrnNFLUuMmfvVKSWqdpRtJapylG0lqnKUbSWqcpZtlvPNGUoss3UhS4yzdSFLjTPQrWL6R1BoTvSQ1zouxktQ4L8ZKUuMs3UhS40z0q9hx/DEvykpqxrok+iRvTvKhJB9P8ob1OMZGMNlLakHPiT7JqSQ3k1xe0b43ybUk80mOA1TVX1fVEeAdwC+ONGJJUl/6mdGfBvYub0iyBTgB7AOmgUNJppe95bc6r0uSxqTnRF9VF4FnVjQ/BMxX1VNV9RxwFjiQJe8DPlFVnx1duJKkfg1bo98KPL1se6HT9qvA64GfT/LIajsmOZpkLsnc4uLikGHc3jD1di/OStrM7hpy/6zSVlX1fuD9a+1YVSeT3AD2T01NPThkHOvGBC9psxt2Rr8AbF+2vQ24PmSfdyyTvqTNaNhE/wSwO8nOJFPAQeB8rztvxk/GWsaRtNn0c3vlGeASsCfJQpLDVXULOAY8DlwFzlXVlT76dK0bSVpnPdfoq+pQl/ZZYHaQg1fVBeDCzMzMkUH2lyTdnqtXSlLjXL1SkhrnjH5AXpCVtFk4o5ekxjmjH4KzekmbgTN6SWqcDx6RpMaZ6CWpcdboJalx1uiH5AVZSXc6SzeS1DgT/Qg4q5d0J7NGPyIuXyzpTmWNXpIaZ+lmHTizl3QnMdFLUuNM9OvEmr2kO8XIE32Slyf5cJKPjbpvSVL/ekr0SU4luZnk8or2vUmuJZlPchygqp6qqsPrEexm4Cxe0p2m1xn9aWDv8oYkW4ATwD5gGjiUZHqk0TXIfwgkbbSeEn1VXQSeWdH8EDDfmcE/B5wFDvR64CRHk8wlmVtcXOw54BaY7CVtpGFq9FuBp5dtLwBbk7w0yQeBVyZ5d7edq+pkVc1U1cx99903RBiSpLUMk+izSltV1b9V1SNV9Yqq+r01O2jok7HdOHuXNG53DbHvArB92fY24Ppw4bTJZC9pnIZJ9E8Au5PsBL4GHATeOpKoJsjyfwS+8vs/PcZIJLWq19srzwCXgD1JFpIcrqpbwDHgceAqcK6qrvRzcNe6kaT119OMvqoOdWmfBWYHPXiS/cD+Xbt2DdrFprXj+GPO4CVtCFevlKTGuR69JDXOGb0kNc7VKyWpcZZuxmiQ++u9J19SvyzdSFLjLN1IUuMs3UhS4yzdSFLjLN1IUuNM9JLUOBO9JDXOi7GbQLd753ccf+z/vbZa21p9rLXPIDHdqTZbvNIoeTFWkhpn6UaSGmeil6TGDfMowVUluQf4U+A54FNV9dFRH0OS1LteHyV4KsnNJJdXtO9Nci3JfJLjnea3AB+rqiPAm0YcrySpT72Wbk4De5c3JNkCnAD2AdPAoSTTwDbg6c7bvjOaMCVJg+op0VfVReCZFc0PAfNV9VRVPQecBQ4ACywl+577lyStn2ES8Va+N3OHpQS/FXgU+LkkHwAudNs5ydEkc0nmFhcXhwijHc/fz77y6/nXetlvtfbV3t9LLL0e53bHG+RY3Y7byxjXiq/XOHrt83b9jer+/V7+2/fb37Cv9/t7dLsxjEIv/78M0udm/xzGMBdjs0pbVdV/Ar98u52r6mSSG8D+qampB4eIQ5K0hmFm9AvA9mXb24Drw4UjSRq1YRL9E8DuJDuTTAEHgfP9dOAnYyVp/fV6e+UZ4BKwJ8lCksNVdQs4BjwOXAXOVdWVfg7uWjeStP56qtFX1aEu7bPA7KAHr6oLwIWZmZkjg/YhSVqbq1dKUuNcvVKSGueMXpIa54xekhqXqhp3DCRZBL66QYe7F/jGBh1r3BxreyZlnOBYe/GjVXXf7d50RyT6jZRkrqpmxh3HRnCs7ZmUcYJjHSUXHZOkxpnoJalxk5joT447gA3kWNszKeMExzoyE1ejl6RJM4kzekmaKBOV6Ls847YJSb6S5ItJPpdkrtP2kiR/l+RLne8/NO44B7HaM4u7jS1L3t85x19I8sD4Iu9fl7G+N8nXOuf2c0keXvbauztjvZbkjeOJun9Jtif5+yRXk1xJ8mud9ubO6xpj3bjzWlUT8QVsAb4MvByYAj4PTI87rhGO7yvAvSva/gA43vn5OPC+ccc54NheCzwAXL7d2ICHgU+w9GCcVwOfGXf8Ixjre4HfXOW9053f47uBnZ3f7y3jHkOP43wZ8EDn5xcC/9IZT3PndY2xbth5naQZfbdn3LbsAPCRzs8fAd48xlgGVqs/s7jb2A4Af15L/hF4cZKXbUykw+sy1m4OAGer6n+q6l+BeZZ+z+94VXWjqj7b+fnfWVrqfCsNntc1xtrNyM/rJCX6bs+4bUUBf5vkySRHO20/UlU3YOmXDfjhsUU3et3G1up5PtYpWZxaVoJrYqxJdgCvBD5D4+d1xVhhg87rJCX6VZ9xu+FRrJ+fqKoHgH3AryR57bgDGpMWz/MHgFcAPw7cAP6w077px5rkBcBfAr9eVd9a662rtG32sW7YeZ2kRN/0M26r6nrn+03gr1j6U+/rz/952/l+c3wRjly3sTV3nqvq61X1nar6LvAhvvdn/KYea5LvZynxfbSqHu00N3leVxvrRp7XSUr0Qz/j9k6V5J4kL3z+Z+ANwGWWxvf2ztveDnx8PBGui25jOw/8UucujVcDzz5fCtisVtSif5alcwtLYz2Y5O4kO4HdwD9tdHyDSBLgw8DVqvqjZS81d167jXVDz+u4r0hv8NXvh1m64v1l4D3jjmeE43o5S1fpPw9ceX5swEuBTwJf6nx/ybhjHXB8Z1j60/bbLM12DncbG0t/9p7onOMvAjPjjn8EY/2Lzli+0EkCL1v2/vd0xnoN2Dfu+PsY50+yVI74AvC5ztfDLZ7XNca6YefVT8ZKUuMmqXQjSRPJRC9JjTPRS1LjTPSS1DgTvSQ1zkQvSY0z0UtS40z0ktS4/wXLmDzxphsguQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "data=MIPS_Map[1].data\n", "plt.hist(data.flatten(),bins=np.arange(-10.0,250.0,1.0));\n", "plt.yscale('log')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in PSF" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": true }, "outputs": [], "source": [ "MIPS_psf=fits.open('../../dmu17/dmu17_HELP_Legacy_Maps/ELAIS-N2/data/dmu17_MIPS_ELAIS-N2_20181029.fits')" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAD5xJREFUeJzt3W+spGdZx/Hvz9aCllr+G9Ltum12qewLKWTSghiDiGZLWWoI0W5IBLPZDYYaTEx0G42JL4zwRrCxohup+4a01qK4S1cLKTSNSQPd5e8u68qh1vSkyG4trAkx1sLli3m2HY/nz8yZmZ2Ze7+fZHLOc88z91xzztOr917Pfe47VYUkqV0/NOsAJEnTZaKXpMaZ6CWpcSZ6SWqciV6SGmeil6TGmeglqXEmeklqnIlekhp36SzfPMluYPcVV1yx79WvfvUsQ5GkhXP8+PGnquoVG52XeVgCodfr1bFjx2YdhiQtlCTHq6q30XmWbiSpcTNN9El2Jzl47ty5WYYhSU2baaKvqiNVtf/KK6+cZRiS1DRH9JLUOEf0ktQ4b8ZKUuMs3UhS4yzdSFLjmizdbDtwP9sO3D/rMCRpLli6kaTGzXStm6o6Ahzp9Xr7JtGfo3hJ+v+aLN1Ikp5nopekxlmjl6TGOb1Skhpn6UaSGmeil6TGmeglqXHejJWkxnkzVpIaZ+lGkhpnopekxpnoJalxJnpJapyJXpIaZ6KXpMY5j16SGuc8eklqnKUbSWqciV6SGmeil6TGmeglqXEmeklqnIlekhpnopekxk0l0Se5PMnxJG+fRv+SpOENleiT3JXkTJITK9p3JTmdZCnJgYGnfge4d5KBSpI2Z9gR/SFg12BDkkuAO4GbgJ3AniQ7k7wV+Drw7QnGKUnapEuHOamqHk6ybUXzDcBSVT0GkOQe4BbgRcDl9JP/fyU5WlU/WNlnkv3AfoCtW7duNn5J0gaGSvRruAp4YuB4Gbixqm4DSPJe4KnVkjxAVR0EDgL0er0aIw5J0jrGSfRZpe25hF1VhzbsINkN7N6+ffsYYUiS1jPOrJtl4OqB4y3Ak6N04OqVkjR94yT6R4EdSa5JchlwK3B4lA5cj16Spm/Y6ZV3A48A1yVZTrK3qp4FbgMeAE4B91bVyVHefNoj+m0H7n/uIUkXq2Fn3exZo/0ocHSiEUmSJsqtBCWpcW4lKEmNc1EzSWqcpRtJapylG0lqnKUbSWqcpRtJapylG0lqnKUbSWqciV6SGmeNXpIaZ41ekhpn6UaSGmeil6TGmeglqXHj7Bk7tknsGeumIpK0Pm/GSlLjLN1IUuNM9JLUOBO9JDXORC9JjTPRS1LjTPSS1DgXNZOkxjmPXpIaZ+lGkhpnopekxpnoJalxJnpJapyJXpIaZ6KXpMaZ6CWpcRNP9Elek+TPk9yX5Ncn3b8kaTRDJfokdyU5k+TEivZdSU4nWUpyAKCqTlXV+4BfBnqTD1mSNIphR/SHgF2DDUkuAe4EbgJ2AnuS7OyeewfwT8CDE4tUkrQpQyX6qnoYeHpF8w3AUlU9VlXPAPcAt3TnH66qnwbevVafSfYnOZbk2NmzZzcXvSRpQ+NsDn4V8MTA8TJwY5I3A+8EXgAcXevFVXUQOAjQ6/VqjDgkSesYJ9FnlbaqqoeAh4bqINkN7N6+ffsYYUiS1jPOrJtl4OqB4y3Ak6N04OqVkjR94yT6R4EdSa5JchlwK3B4lA5cj16Spm/Y6ZV3A48A1yVZTrK3qp4FbgMeAE4B91bVyVHe3BG9JE3fUDX6qtqzRvtR1rnhKkmaPbcSlKTGuZWgJDXORc0kqXGWbiSpcZZuJKlxlm4kqXGWbiSpcZZuJKlxlm4kqXEmeklqnDV6SWqcNXpJapylG0lqnIlekhp30ST6bQfuZ9uB+0d+TpIWnTdjJalx3oyVpMZdNKUbSbpYmeglqXEmeklqnIlekhpnopekxpnoJalxzqMf4B9OSWqR8+glqXGWbiSpcSZ6SWqciV6SGnfprAOYFW+6SrpYOKKXpMaZ6CWpcVMp3ST5JeBm4JXAnVX16Wm8z2ZYspF0sRl6RJ/kriRnkpxY0b4ryekkS0kOAFTVJ6tqH/Be4FcmGrEkaSSjlG4OAbsGG5JcAtwJ3ATsBPYk2Tlwyu91z0uSZmToRF9VDwNPr2i+AViqqseq6hngHuCW9H0I+Ieq+uLkwpUkjWrcm7FXAU8MHC93bb8BvBV4V5L3rfbCJPuTHEty7OzZs2OGIUlay7g3Y7NKW1XVHcAd672wqg4CBwF6vV6NGYckaQ3jjuiXgasHjrcATw774nlbvVKSWjRuon8U2JHkmiSXAbcCh4d9satXStL0jTK98m7gEeC6JMtJ9lbVs8BtwAPAKeDeqjo5Qp+O6CVpyoau0VfVnjXajwJHN/PmVXUEONLr9fZt5vWSpI25BIIkNc6tBCWpcW4lKEmNc0S/CjcJl9QSR/SS1DhvxkpS40z0ktQ4a/SS1Dhr9JLUOEs3ktQ4E70kNc4avSQ1zhq9JDXO0o0kNc5EL0mNM9FLUuO8GStJjfNmrCQ1ztKNJDXORC9JjRt6c/CL0eDmI49/8OYZRiJJm+eIXpIaZ6KXpMY5vVKSGuf0SklqnKUbSWqciV6SGmeil6TGmeglqXEm+gnYduD+//PHVZI0T0z0ktQ4E/2Qhhm1O7KXNI8mnuiTXJvkY0num3TfkqTRDZXok9yV5EySEyvadyU5nWQpyQGAqnqsqvZOI1hJ0uiGHdEfAnYNNiS5BLgTuAnYCexJsnOi0UmSxjZUoq+qh4GnVzTfACx1I/hngHuAWyYcnyRpTOOsR38V8MTA8TJwY5KXAX8IvC7J7VX1R6u9OMl+YD/A1q1bxwhjNrzpKmlRjJPos0pbVdV/AO/b6MVVdRA4CNDr9WqMOCRJ6xhn1s0ycPXA8RbgyVE6cJliSZq+cRL9o8COJNckuQy4FTg8SgcuUyxJ0zfs9Mq7gUeA65IsJ9lbVc8CtwEPAKeAe6vq5Chv7oh+bf7xlaRJGapGX1V71mg/Chzd7JtX1RHgSK/X27fZPiRJ63MJBElq3DizbsaWZDewe/v27bMM44JYrQzz+AdvnkEkki427hkrSY1zRD+iadwgPd/nLEf4Kz/XuLEM9jdKX5t9naS1OaKXpMZ5M1aSGmeil6TGWaOfsvVq+iufW+3c9WbrjDqTZx7uBUi68KzRS1LjLN1IUuNM9JLUOGv0UzBPi5ENW/cfpa9JzrE/b1bz9sd5n5VWe1/vi2geWKOXpMZZupGkxpnoJalxJnpJapw3YxfQrG/2LvLCYy3fHF3v9zKJ39msf3azuu4W+Xo/z5uxktQ4SzeS1DgTvSQ1zkQvSY0z0UtS40z0ktQ4p1dqVZOcUjaptXVG7WeY6YAXeo2cYd5j2JiG+Xmst9bRND7vZvte7XWbuW5G/X2uF++405inPd11FE6vlKTGWbqRpMaZ6CWpcSZ6SWqciV6SGmeil6TGmeglqXEmeklq3MT/YCrJ5cCfAc8AD1XVxyf9HpKk4Q01ok9yV5IzSU6saN+V5HSSpSQHuuZ3AvdV1T7gHROOV5I0omFLN4eAXYMNSS4B7gRuAnYCe5LsBLYAT3SnfX8yYUqSNmuoRF9VDwNPr2i+AViqqseq6hngHuAWYJl+sh+6f0nS9IxTo7+K50fu0E/wNwJ3AH+a5GbgyFovTrIf2A+wdevWMcLQoAu9r+fKhZ+muVDWJM9fb7GvzRp3EbVhFiAbNsbN/gxWxrbZRdHGXchu1PM3u4jaqO+73u9j1AXrLqRxEn1Waauq+h7waxu9uKoOAgcBer1ejRGHJGkd45RWloGrB463AE+O0kGS3UkOnjt3bowwJEnrGSfRPwrsSHJNksuAW4HDo3TgMsWSNH3DTq+8G3gEuC7JcpK9VfUscBvwAHAKuLeqTo7y5o7oJWn6hqrRV9WeNdqPAkc3++ZVdQQ40uv19m22D0nS+mY6/dERvSRNn1sJSlLj/IMmSWqcpRtJalyqZv+3SknOAv82wS5fDjw1wf4ulEWNGxY39kWNGxY39kWNG+Yv9p+oqldsdNJcJPpJS3KsqnqzjmNUixo3LG7sixo3LG7sixo3LG7s1uglqXEmeklqXKuJ/uCsA9ikRY0bFjf2RY0bFjf2RY0bFjT2Jmv0kqTntTqilyR1mkr0a+xhOzdW23s3yUuTfCbJN7qvL+nak+SO7rN8NcnrZxj31Uk+l+RUkpNJPrBAsb8wyReSfKWL/Q+69muSfL6L/a+7FVhJ8oLueKl7ftusYu/iuSTJl5J8asHifjzJ15J8Ocmxrm0RrpcXJ7kvyT931/sbFyHujTST6LP2Hrbz5BAr9t4FDgAPVtUO4MHuGPqfY0f32A989ALFuJpngd+qqtcAbwDe3/1sFyH2/wbeUlWvBa4HdiV5A/Ah4MNd7N8B9nbn7wW+U1XbgQ93583SB+ivDnveosQN8HNVdf3AdMRFuF7+BPjHqvpJ4LX0f/aLEPf6qqqJB/BG4IGB49uB22cd1ypxbgNODByfBl7Vff8q4HT3/V8Ae1Y7b9YP4O+BX1i02IEfBb5If8vLp4BLV1479JfdfmP3/aXdeZlRvFvoJ5a3AJ+iv6vb3MfdxfA48PIVbXN9vQA/Bvzryp/bvMc9zKOZET2r72F71YxiGcWPV9W3ALqvr+za5/LzdCWB1wGfZ0Fi78ofXwbOAJ8Bvgl8t/p7KqyM77nYu+fPAS+7sBE/5yPAbwM/6I5fxmLEDVDAp5McT39/aJj/6+Va4CzwV1257C+TXM78x72hlhL9qnvYXvAoJmfuPk+SFwGfAH6zqv5zvVNXaZtZ7FX1/aq6nv4I+QbgNaud1n2di9iTvB04U1XHB5tXOXWu4h7wpqp6Pf3yxvuT/Ow6585L7JcCrwc+WlWvA77H82Wa1cxL3BtqKdGPvYftjHw7yasAuq9nuva5+jxJfph+kv94Vf1t17wQsZ9XVd8FHqJ/n+HFSc5vvDMY33Oxd89fCTx9YSMF4E3AO5I8DtxDv3zzEeY/bgCq6snu6xng7+j/D3ber5dlYLmqPt8d30c/8c973BtqKdGPvYftjBwG3tN9/x769e/z7b/a3dl/A3Du/D8fL7QkAT4GnKqqPx54ahFif0WSF3ff/wjwVvo32D4HvKs7bWXs5z/Tu4DPVleAvZCq6vaq2lJV2+hfy5+tqncz53EDJLk8yRXnvwd+ETjBnF8vVfXvwBNJruuafh74OnMe91BmfZNgkg/gbcC/0K/B/u6s41klvruBbwH/Q380sJd+HfVB4Bvd15d254b+LKJvAl8DejOM+2fo/5P0q8CXu8fbFiT2nwK+1MV+Avj9rv1a4AvAEvA3wAu69hd2x0vd89fOwXXzZuBTixJ3F+NXusfJ8/8tLsj1cj1wrLtePgm8ZBHi3ujhX8ZKUuNaKt1IklZhopekxpnoJalxJnpJapyJXpIaZ6KXpMaZ6CWpcSZ6SWrc/wIcEwsYdZZ4+gAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "dat=MIPS_psf[1].data\n", "plt.hist(dat.flatten(),bins=np.arange(-10.0,650.0,5.0));\n", "plt.yscale('log')" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "ename": "NameError", "evalue": "name 'np' is not defined", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmean\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m~\u001b[0m\u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0misnan\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdat\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mNameError\u001b[0m: name 'np' is not defined" ] } ], "source": [ "np.mean(~np.isnan(dat))" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": true }, "outputs": [], "source": [ "centre=np.long((MIPS_psf[1].header['NAXIS1']-1)/2)\n", "radius=20" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAASwAAAD8CAYAAADNNJnuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJztnX2wJFd53p93Pu69uysJIVZIa63kFfKWLXDsxWwUElIVEMJWFIIgMS6TD6tiJTIpqMJ8xJbYxB9JSCCA5C+CszYEuQpbYDAlFSaGRYai/IcFK1gJybItoQiQtNFqkVa6+3HvnZl+88f0FffOec7e09MzPd13n99W1945c0736Z6eM6ef97zva+4OIYRoAq1Zd0AIIVLRgCWEaAwasIQQjUEDlhCiMWjAEkI0Bg1YQojGoAFLCDFVzOxqM/sbM3vIzG4stS+twxJCTAszawP4WwCvAfAogK8BeJO7/9U4+9MMSwgxTa4A8JC7P+zuKwBuA3DtuDvrlOmJmV0N4DcBtAH8vru/93T1t2/f7rt27SpzyIlz/2NP0HJvkzIjFWlZZNZasr2xurReensDqZt4nMKQQ0V6ypsnPgw4/aAi7Vnd1HrRDpCyjFdt9cOyl+y8IP1YI9x9991H3f38sXcA4Kdetc2/99Qg7Xj3Lt8PYGlN0X5337/m9UUAvrvm9aMA/t64fRt7wMqneh/Cmqmemd1xuqnerl27cPDgwXEPORUu/4+30PL+tvCuG8yHZdlcWOadyDdrLrxrrUvK2rx9ux3WbZGybpffbJ1WWLdD2rdJPYAPhLHBgZGRuoMsLMsyPvHvk/KsSPte+Cs06IV1fSmsZ/3Iw0ji4NQ+xdsvPBn2/+D/eDs/VgJm9u2xG+d876kBvvr5S5Lqtnc8uOTue0/XJVI2tg5V5pFwolM9IUQ9cABZ4r8EHgVw8ZrXOwE8Pm7fyjwSJk31zOwGADcAwCWXpI3aQojZ4XD0PO2RMIGvAdhtZpcCeAzAzwL4F+PurMwMK2mq5+773X2vu+89//xSj9ZCiIqY1AzL3fsA3grg8wAeAPBJd79/3H6VmWFNdKpXBZfvC/Wq5e38og+2kHKmQRG9iulKANDupGlQTKuKlbeJrtTt8F/HbjssZ+1jGlaL1GW6VAymdzFdahDZ54DVJRpWf0AsJuAa3FCKHWnPdPheRHYhx8eA6Gp93j7rTsvCMT4Ox2CCy53c/XMAPjeJfZWZYT031TOzOQynendMolNCiNmSwZO2qhl7huXufTNbneq1AXy0zFRPCFEPHMBgBoNRCqXWYU1yqieEqA+zmD2lUGrAEkJsPhxAr6YuexqwhBDrcPjmfCSsMy/7tzcHZcsvIavXz+YWtdaW0Gci1crXavEPu0Osd2z1eTvSnlnvmOWrG7Hysfb0+JZuJSxrOWSWv75zWxCty8oiVlbWV+rGRMoGLW55dGIR9GVSN3JJWj1ePlMcGNRzvNq8A5YQYjyGK93riQYsIcQIhsHUPODLoQFLCLGOoeiuAUsI0QCG67A0YFXKs5eGZYOzQtG7tZUEJAIwNx+WM9GcurvERHdSlwnBzF1muN80ZYEJ6QDQIpYfJjrH2neIGN+KCPSp9ImY3Y38unM3nrCsV8g1h5WFbXuR7++AuPYMiBtOxI6ArMvLZ00Rw0mVbNoBSwgxHpphCSEag8MwqGn0dA1YQogAPRIKIRqBw7DCkhrUgMYPWD/x78IV7QAw+BFS2CUCdyx2VDcU3efZSnUipBdZaZ66ehzgojkjloSC7rPI8anont6e7pP4rMV+3VtE4M6c9CmW2iKxr6niPACskDKfC8uyef6I1Tu7fjOZ4cJRPRIKIRqCRHchRCNwN7pcpA5owBJCBGSaYQkhmsBQdK/n0FA28/MjABYBDAD0N0ioKIRoAJtddH+Vux+dwH42ZO/PhxbBxct43cE24jKzQCx/xAUHABaIlXCOWQRJWcy1JRZnapQiVrZZr5cpZBFMrBu1UiZaFGPtLTHIE7uirchlTrXIrkQyRw+2huUv/ffhff6ND78j6TiTIpa5aNbUc94nhJgZdV7pXrZXDuALZnZ3nuFZCLEJyLyVtFVN2RnWK9z9cTN7IYADZvbX7v6VtRWUql6IZjF0ft6EMyx3fzz//wiAzwC4gtRRqnohGoTD0PN20lY1Y8+wzGwbgJa7L+Z//ySA/3y6Ng888gRedv16QfH4D0ZiH20hiQB+IhStnaWUB08isbAQRvxn4joALHSYaw5JTEGE9JjoXpZUgZ2lhI/us8BxmOWoT3YQO/8iAnkqtH3s8pPvF2vP+h8V8kk5u/7ZgM8N+tvC8sVLw6/lnrdwF7RDH5q8GO/O44zVgTKPhBcA+IwNo511APyhu//ZRHolhJghtvkWjrr7wwB+fIJ9EULUAMfmnGEJITYpdRXdNWAJIdbhsJkvSI5R6YA1mAcWd62/ECvP5/GomJhuc2SlOSkblocC+Za5UHTf0uWpd5nA3m2lrWqPxWNiugC7MaKiNymnAm8B/YHly2DJHgCg1SaVyaMDE+KBCax0T2wfTcJBVsq3iBJf1hAwyMLrP5jj15R9fj2STfrZLm+/+7+vF+PnL9r5spQ+no5hmq96zmXq2SshxAxRIlUhRENwYCar2FPQgCWECNAMSwjRCNyttjOsevZKCDEzhqL79F1zzOyNZna/mWVmlhRLr9IZlneAlRest7T5Nu4a05lnFkHiLhNxremSDDdbiUVwvs3bMyshs/6xTDJFfp2KWAlTF/PFs86kZYMpayUrQlnzebGsP2l1WxEraSosllQW2SezEvo2ss82/6qOZuOZjHtfZTHd7wPwzwD8r9QGeiQUQqxjKLpPX8Ny9wcAIHfvS0IDlhAioMBK9+1mdnDN6/3uvn8KXQKgAUsIMULBle5HT5fLwcy+COBC8tY+d7+9aN80YAkhAiaVhMLdr5rIjnKqHbDM4VvWi+FdkhgCAOZIcgjmWsPE8Vg5E9jniLsNAHQi5SlkEc06VWCPxbNioj8Tzefb3DWlQ1PNp8fzShXji7gWFTE6lHXtSa3Lzj/zSIwvcp/MtcMve48YgYCYQE/cdSLXJOj9BAwm7kCvpOFhWtSzV0KImTF8JJx+THcze4OZPQrg7wP4UzP7/EZt9EgohAioYqW7u38Gw9DqyWjAEkKso6plDeOgAUsIMUJ9XXM2HLDM7KMAXgvgiLv/aF52HoBPANgF4BEAP+PuT294NANsJKZSKyIQs4zKc4mJIQAusC+0Q9G+iMBcRCBOha60jvSJieZzBQwJ862wLjMutAsIt0w07md8uXWPfAlWsvAWjMXjSr3+ZeNp8baRz4SkhGbXv0fuZyASz4ycfzSJxWhG6QlNjOoa0z1lGP0YgKtHym4EcKe77wZwZ/5aCLEJGFoJ20lb1Ww4YOWJUZ8aKb4WwK3537cCeP2E+yWEmBGrC0dTtqoZV8O6wN0PA4C7H84zP1PWZn5uv+DcMQ8nhKiSuj4STl10z/2K9gPA/KU7qwsDIIQYi81oJXzCzHbks6sdAI4ktTKHddaLl+3YqmxS3iVidCw8DBPY54jwGUsYwUhOIlHg14mJ/qzvALCtsxKWtcOyLaQM4KJ715jonp65msVEWsq6tO6pQVh+oj8f1gNvvzQIb1dmzYp5GjDhvHRiDJLYos3u04hxiBktWq3wmjKPBgBhlusJTQnqaiUct1d3ALgu//s6AIWdGIUQ9cTd0PdW0lY1Kcsa/gjAKzEMI/EogF8F8F4AnzSz6wF8B8Abp9lJIUS1NPaR0N3fFHnr1RPuixCiBmxGDUsIsYnRgCWEaARKVb+G1khu9LiVkCShYO460XhWaSnkY9afaXxgzLVmoRNaBJ/XPUXbn9NZCut2wrpntcN6ALBg4bGKWAlZYgJqJXRu5Ts+WAjKFolF9FhvC20PhO2Z5bAfycSQkXJm5SsSC43HEyOp6rOICxq5T9l93uvw9r0RV7fN7pqjGZYQYh3ucX/OWaMBSwgRoEdCIUQjkIYlhGgUsbwCs6bSAcsMsBGRMRZ7iZUzMTgap6iEGwZQwA2HCNEdImQDwAJxIzq7GwrkL+ieoO2fT8rPa4dl21rLtP1WUj4X6WsqK0TIPpmF7jYAcKIVnut8K0xz3A78TeJkCAX6bDRGVA5dmU2+l6VjnBHjDnPXAbjA3m0T40Ck/Wh8OYnuQogzCndpWEKIxmAYyEoohGgK0rCEEI1AvoRrGF3pHk24QMpTy4BIcgcixJeN+8P2yeJuATye1blkVfv27nHafnvn2bB9+2R4nIjovs3C47OV7kVgK91PxER/IsZ3LTREtCNBnQYkGhKPhxXJPE3icfF4ZuE+Y8YdBo+7xe8zdp8ygT7mEdIOVsBPICCWD3WsOqIZlhAiQFZCIUQjcInuQogmoUdCIURjqKuVcMN5n5l91MyOmNl9a8p+zcweM7ND+XbNdLsphKgK9+GAlbJVTcoM62MAfgfAH4yU3+LuHyjbARI6CMBpsoQkwi1FBSxKiR8Gs1JuiWS9OYfErjqvQ9xtOtxKyCyCrOwc41a6eWIRnGMWLdo6TNACACvESrfgPENMauwtZqUDuEWShUFZjmQkXmFZd5gLFhOcI9bkVOthNOsOc+MhdZkLDwD059Zfa2tN5lmuscsa3P0rZrZr+l0RQtSFumpYZUwBbzWze/NHxudPrEdCiJniMGRZK2mrmnGP+GEAlwHYA+AwgA/GKprZDWZ20MwODp7lUQiEEPXCE7eqGWvAcvcn3H3g7hmA3wNwxWnq7nf3ve6+t31OGEpECFEzGi66B6ymqc9fvgHAfaer//12HrgYxMR17lozeSE+lr02WXQnomtMdD+rE4rhz0sU0mPl55IYU9uIuwsAzJNTmjOSKp225qL7wMPSk6QM4AkfGL02vy2XstC15mRnLig7PuDxuJbaYft+P2xPP/vI7VDVQ1HMODXq6jaxIaSmGta4mZ9faWZ7MDytRwD8whT7KISomCpmT2b2fgD/FMAKgG8B+Dfufux0bcbN/PyRsXoohKg9DiDLKnncOwDgJnfvm9n7ANwE4JdP16CeDkNCiNnhANzStjKHcf+C+3OL9v4SwM6N2sg1RwgRUGAd1nYzO7jm9X533z/GIX8ewCc2qlT5gDUquseC85cW2BPlxyIr3VmfWJbguRYXvbe2wnhULDHE1shKdRbPignsWyOnvkBiMnVJWSsy8c6I7D5gn1MWi7EV9rVHzmkpcv5nk4zWi1lYtqXNM0cvJiYmYavn4xnCq3lIiRqnIt+f0qR//Y66+97Ym2b2RQAXkrf2ufvteZ19GN4cH9/oYJphCSFGmNySBXe/6rRHMrsOwGsBvNp943mdBiwhREgFyxrM7GoMRfZ/5O58Lc8IGrCEEOtxwKuxEv4OgHkAB2y4HvAv3f3Np2ugAUsIQZj+gOXuP1S0jQYsIURIU1e6T5rRcbtInCBWVmXcHprhhJTNR6yE863QZWchsQwAFkg8KeZuw6yBw7rhx901khY98uuakevf87BPC5HPdECsjOycouefkWtFrKSx688suuyeYsTus1Rr9rTu06nd/RqwhBCNYHXhaA3RgCWECKhrAD8NWEKIkGqshIXRgCWECCjpaDI1Kh2wzMq73EyauGsOE65ZwgYixNPIUTwJwxwrA3dt6RLXEpZugbnbDMvD2qwsRmpN5sIDAHPkWrFzip0/S1jBkkDErj8zkJSlbIbksu2n8m2aVTjRBDTDEkKMUD4Sw7TQgCWECNEMSwjRGKYUBKIsGrCEEOtp8josM7sYw6zPF2I47u539980s/MwDLi1C8O47j/j7k9Pr6ucKkX81NXKsT4xMbhV8qesTZJIbFbYtaLXtID3RCqzXukeO860Pv2a2caeIyX6WB/AO939cgAvB/AWM3sxgBsB3OnuuwHcmb8WQmwGapqYcMMBy90Pu/vX878XATwA4CIA1wK4Na92K4DXT6uTQggBFNSwzGwXgJcCuAvABau5Cd39sJm9MNLmBgA3AED3/HPK9FUIURFNfiQEAJjZWQA+DeAX3f3Z1HZrMz93nqfMz0LUHsfQNSdlq5ikAcvMuhgOVh939z/Ji58wsx35+zsAHJlOF4UQlVNTDSvFSmgYJk59wN1vXvPWHQCuA/De/P/bx+lA1DWG2D9oWUnrzTSI9WlAfh8yUjaI2H4GZL8sVTzNZINI1huyz3bEtYcdi+4zcicPaAiA9F9pfq3INZ3CPRXPmpNet0nU9RRSNKxXAPjXAL5pZofysndjOFB90syuB/AdAG+cTheFEJXT1AHL3f8C8Z/BV0+2O0KIWtDUAUsIcWZh3uxHQiHEmYYC+A3Drqa6KFSZXCIV1icm2jIhGAB6HkaUWiFlPecfS4/st0dE7x4RxwGgzeJMkcucReLjMoGdJaGItWepJXok7thKJPIWv1ZpZQDQy8Lyadxndbx3i6IZlhCiOWjAEkI0AmlYQohGoQFLCNEUSKj8WlD5gFXTgXtsmJC7nPHLupx1g7IlVtYKywBgiYjJS0T0jsd9CuuyVent2Ep7mvk5vLOXIqL7EhGj2TmdzOZ5e3Kt2DWNXf8+TSySRjTGVkXxsOLJUpov8BdBMywhREhNZxYasIQQ65HoLoRoFBqwhBCNQQOWEKIJGGQlfI4sW2+p8QLWj9QyIN16E49zxMpCKxOzPK1ErFQns7mgbDHbEpRty1Zo+wULnVva7KewxZxggIyliidWxiKwIzFrIACcJNdlMVsIyk5ErITsWi0OwvanBuF1BoCVQZprDrsnimTiKZJ+PvWeHkR8+7Jp+PxJwxJCNIqaDljjL0wRQmxeKgiRbGb/xczuNbNDZvYFM/uBjdpowBJCBKzGxNpoK8n73f3H3H0PgM8C+JWNGmw4YJnZxWb2JTN7wMzuN7O35eW/ZmaP5aPjITO7pnT3hRD1oIIZ1kj2rW0pe0zRsFYzP3/dzM4GcLeZHcjfu8XdP5DcQYQi4SDjYyYT47lAyduzmExFkgO0iJmEiu6k/6cG3LXmeD8Uk7e2QtH4GBHXAaBN+sQSVsSSWCxYKLB3SYyrIrAYXczdBuAuN8eyrWHZICwDgGf6YfnxQbjP2PVnn1+ycabkt7OIaw1LDDJqrHqu7iDNiFUIL2Ql3G5mB9e83u/u+1Mbm9l7APwcgGcAvGqj+ikx3Q8DWE2Yumhmq5mfhRCblfTx+ai77429aWZfBHAheWufu9/u7vsA7DOzmwC8FcCvnu5gZTI/vwLAW83s5wAcxHAW9nSR/Qkh6smkljW4+1WJVf8QwJ9igwGrTObnDwO4DMAeDGdgH4y0u8HMDprZwcEzJ1MPJ4SYJdVYCXevefk6AH+9UZukGRbL/OzuT6x5//cwVPkD8ufZ/QCw8EM/UNPVHUKI56guq/N7zeyHAWQAvg3gzRs1GDvzs5ntyPUtAHgDgPs27J4b+v31k7rYCl4mZjOBPmuVE42LCPFs9XufxMM60ecrrTtEyewU6D9LbrHSTk9iwVbKd62ffHwGO9aSc9H7WbKqfXEQrl4/2j+btn+aiO6LvXCfSwN+/qkr0MtmbmZCekwMZ/d0n6zI75EyAOj315dHQpEVwlDNSnd3/+dF25TJ/PwmM9uD4Vj8CIBfKHpwIUQ9aaxrzmkyP39u8t0RQtSCpg5YQogzEA1YQohGoGgNQohGoQFraCnpr6w/5Mocj8fUaYfWsy6xiLUzbmVrkXTnnRbJMDMFd52ViJVq0Xicp3CfEYsSORbLGnOyfYq251bC8JowF6DY8Vla+JiVkMW5Yu42zBoIAE+tbAvKFom7UyweWXLsqwLfVmZ5ZMdZIdZkAOhRKyEp60dc2AbTyZqjAH5CiMagR0IhRDOobuFoYTRgCSFCNGAJIZpAVSvdx6HaASsD/NR68XGly7vQJi4r3XYoEMdcW6iYStKqF0n1zQT2IgkHlvpcjE7tExO4aWILkpgBABZIcgomujPjwrBfxF2KnP9JEqMK4H19tk/cdUgZACz20gT2IolJUgX22OfMjtUjAnvMtWaln1Y3I0I8APho+YRS1xvzQ6sBmmEJIdYjDUsI0ST0SCiEaA4asIQQTUEzLAAYGDrPrD9krxvJqEtWuneI6N4ukJG3iOgaE55TiCYcYAI1iZ3FYoEBwDIRmE+0QyF6S5tnjp5vhbGv2DVpF0hMwQwBrJ8Az8jMEkacjBgnYivYRymSpTnVeyFqCEmMXRX7TGldsqo9loTCp5H5GdAMSwjREIplzakUDVhCiHVoHZYQollMItbyFNCAJYQIqOsMKyVV/YKZfdXM7slT1f96Xn6pmd1lZg+a2SfMjGdeEEI0i9QUXzMY1FJmWMsArnT343m6r78ws/8D4B0Ypqq/zcx+F8D1GOYqjNJaAbY+vt6qsbgl4rLQCS1FnU6oBLZb/Kox1x4bhHWjrj3EoketTOxTixhumKWJWY+ySDwplqGHZYg50Ypk7SHxwLpEXY1Z2Vj/meWTuaYAPCYUix3WJ1a62PHLxrOix0mMcQXwvrIYV8wFBwiz3gBARix/xKssrzxaMVKvIHUV3TecYfmQ4/nLbr45gCsBfCovvxXA66fSQyFE5ViWtlVNUuZnM2vnKb6OADgA4FsAjrn76sKeRwFcFGn7/czPp05Mos9CiGniGIruKVvFJA1Y7j5w9z0AdgK4AsDlrFqk7X533+vue9tbwhC3Qoj6YZ62VU3SgLWKux8D8GUALwdwrpmtChA7ATw+2a4JIWZGU0V3MzsfQM/dj5nZFgBXAXgfgC8B+GkAtwG4DsDtG+3rJRdfgIM3v31d2e7/djOty+JkrcyFriXtFj8FFjuLCbSxFOJMtKcCb6EkFomxl2ICL/l9YQJtNB4UEb2nkZY9LlATgblA/Kay1z+1T8wQEjMkLPfD+2+pR+5dUg8A+r1wvwPmmkPqDTs2Gg+LVytC0xeO7gBwq5m1MZyRfdLdP2tmfwXgNjP7rwC+AeAjU+ynEKIq3JsbwM/d7wXwUlL+MIZ6lhBis1HP8Uor3YUQIU1+JBRCnEk4gKY+Ek6bB9/9Dlq++z2hGL8yH64Ab0VWuqeKsd3ISvcu0jIis5XyReIxsVXxhdoXoIhAXu446cbnsqJ52SzNTGAfkDImrsfKl3vhfbq0xL0XBktkv72wnzYqrue0lkbqTio+Vj3Hq9kPWEKI+qFHQiFEY2islVAIcYZR4zRfhVa6CyE2P8OFo560TeR4Zu8yMzez7RvVre0M68F9oRh/2QdCIX45IhqzFexMYJ7vhqvnAWCB1GWr5xlzkXos4QH7xYiK7onCQhEhezqie3p4Gl4v/Xe0UCiYxBXsLDHEKSKkA8CplTSBfbDI27ePE++D5bD/nRP8nOaeHdkfzz9SnIoiMZjZxQBeA+A7KfU1wxJCBFQ4w7oFwC8h8SG0tjMsIcSMqEjDMrPXAXjM3e8xS5t9a8ASQoxQyJdwu5kdXPN6v7vvX31hZl8EcCFptw/AuwH8ZJGeacASQoSkP+4ddfe98d34VazczP4OgEsBrM6udgL4upld4e7/L7Y/DVhCiPVUkEjV3b8J4IWrr83sEQB73f3o6do1asD61rtCy+GLPhiJp0XK2I9GzKKUamVkVzBq5SNuRCwxRDwxxuTjQTErW1myiB6R7hrEzz81nhaz8gHAMkl4wRJGsHhWS8QaCADLy2kWwbmj/Kv24E3cNW1c7LfecfdEdqS8hEKIxlDxeOXuu1LqacASQgRYVs88XxqwhBDrcVS2cLQoZTI/f8zM/q+ZHcq3PdPvrhBi2hjSFo1OyjWnCGUyPwPAf3D3T52m7dR5+J1ctLzs/aEY32OeSiVdU5jo3YmYWLKImF7mWMzdp9g+SzXnbjSx+5gciwvx6ZmfmcDOMkwDwDKpy8R0VrayzL8q2amwvPNMWHbeffUUsaM0VXR3dwfAMj8LITYrNR2wxsr87O535W+9x8zuNbNbzGx+ar0UQlTHqoaVslXMWJmfzexHAdwE4EcA/F0A5wH4ZdZ2bar6J598ckLdFkJME8uypK1qxs38fLW7H/YhywD+NyIpv9amqj///PNLd1gIMW18+EiYslXM2JmfzWyHux+2oSPQ6wHcN+W+FmLr46FAu7iNZPnt8l8JI6vSjYjeqRmiY+2LrF7nq9LD35yyQvy0SF3p3o+I7kxgZ6vXlyIJI04SMf3UyVDJ6BMh3U5wIb97Muzrtu+G53TXxye7on2qOGqrYZXJ/Pzn+WBmAA4BePMU+ymEqJJ6/uaVyvx85VR6JISYObNYY5WCVroLIUI0YAkhGoE7MKjnM6EGLCFEiGZY1fLNW94elLHYWYN5bv3pEyshi2d1qkCfmJWQEYtRNdcKM/xwd51IWvNEB4WolZJZ+UpmrVnJSIyqAvGsmLsNy2QD8Gw2/ZPhPtvHwrL57/Fr2j0Rlt3z2+G91zg0YAkhGoEDUOZnIUQzcMClYQkhmoBDorsQokFIw5o9nZOhGDzYFkmY0AlF1n4rFHhThfQYAyJEDzwU14fHZ2446fG4mMtO6SQW1LUmYsgg50oTQ0RE9xXmmkMSRqyQMgDok5hWdiIs2/p42M/7PrAJhPQiaMASQjSD2Tg2p6ABSwixHgegJBRCiMagGZYQohnINacWMC3a+pHMz71QeM064Q76vVAIZlmjAR4Iv9cOj8OyEQPAfDvMEs3iccVEd1aXUUSIZ0L6oFA8K5JEoh8R3Umcq2UipPfJinYAwDNh+fMeDPt6z2+dYQL7KA641mEJIRqDVroLIRqDNCwhRCNwl5VQCNEgajrDSs6ak+cm/IaZfTZ/famZ3WVmD5rZJ8xsbnrdFEJUh8MHg6StaorMsN4G4AEA5+Sv3wfgFne/zcx+F8D1AD484f6NzWtabwzK7D/9g7BsELHokfKMWAQZWcTKNyD77HfCffYiN8IKsRJ22+HUPWYNbNPYWWnZgWJw1yJ+TamVj1gEezHXGnL9+8vkMznFP6e5xbBfOw4coXXPaGocXiY18/NOAP8EwO/nrw3AlQA+lVe5FcNUX0KIzYBnaVvFpM6wfgPALwE4O3/9AgDH3J/z0n0UwEWsoZndAOAGALjkkkvG76kQohIcgDd1hmVmrwVwxN3gdlpIAAADZElEQVTvXltMqtIzVOZnIRqGe6NnWK8A8DozuwbAAoYa1m8AONfMOvksayeAx6fXTSFElcxCUE/BvID50sxeCeBd7v5aM/tjAJ9eI7rf6+7/c4P2TwL4dv5yO4Cj43W7tuicmsFmPqcfdPdSjzJm9mf5/lI46u5XlzleEcoMWC8CcBuA8wB8A8C/cvflAvs66O57C/a31uicmoHOqbkUWjjq7l8G8OX874cBXDH5LgkhBCd54agQQsyaWQ5Y+2d47Gmhc2oGOqeGUkjDEkKIWaJHQiFEY9CAJYRoDJUPWGZ2tZn9jZk9ZGY3Vn38SWBmHzWzI2Z235qy88zsQB694oCZPX+WfSyKmV1sZl8yswfM7H4ze1te3tjzMrMFM/uqmd2Tn9Ov5+WNjzRypkZPqXTAMrM2gA8B+McAXgzgTWb24ir7MCE+BmB0sdyNAO50990A7sxfN4k+gHe6++UAXg7gLfln0+TzWgZwpbv/OIA9AK42s5fj+5FGdgN4GsNII01jNXrKKpvhnDak6hnWFQAecveH3X0Fw4Wn11bch9K4+1cAPDVSfC2GUSuABkavcPfD7v71/O9FDL8MF6HB5+VDjucvu/nmaHikkTM5ekrVA9ZFAL675nU0ykMDucDdDwPDLz+AF864P2NjZrsAvBTAXWj4eeWPTocAHAFwAMC3kBhppMasRk9Z9T5Ojp7SdKoesJKjPIjZYGZnAfg0gF9092dn3Z+yuPvA3fdg6KB/BYDLWbVqezU+ZaOnNJ2qY7o/CuDiNa83U5SHJ8xsh7sfNrMdGP6iNwoz62I4WH3c3f8kL278eQGAux8zsy9jqM81OdLIGR09peoZ1tcA7M4tGnMAfhbAHRX3YVrcAeC6/O/rANw+w74UJtdBPgLgAXe/ec1bjT0vMzvfzM7N/94C4CoMtbkvAfjpvFqjzsndb3L3ne6+C8Pvz5+7+79Eg8+pEO5e6QbgGgB/i6GWsK/q40/oHP4IwGEAPQxnjddjqCPcCeDB/P/zZt3Pguf0DzF8jLgXwKF8u6bJ5wXgxzCMJHIvgPsA/Epe/iIAXwXwEIA/BjA/676OeX6vBPDZzXROG21yzRFCNAatdBdCNAYNWEKIxqABSwjRGDRgCSEagwYsIURj0IAlhGgMGrCEEI3h/wPAg+VanjwH1QAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import pylab as plt\n", "plt.imshow(np.log10(MIPS_psf[1].data[centre-radius:centre+radius+1,centre-radius:centre+radius+1]/np.max(MIPS_psf[1].data[centre-radius:centre+radius+1,centre-radius:centre+radius+1])))\n", "plt.colorbar()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Set XID+ prior class" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "prior_MIPS=xidplus.prior(MIPS_Map[1].data,MIPS_Map[2].data,MIPS_Map[0].header,MIPS_Map[1].header,moc=Final)\n", "prior_MIPS.prior_cat(masterlist[1].data['ra'][good],masterlist[1].data['dec'][good],masterfile,flux_lower=MIPS_lower,\n", " flux_upper=MIPS_upper,ID=masterlist[1].data['help_id'][good])\n" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": true }, "outputs": [], "source": [ "prior_MIPS.set_prf(MIPS_psf[1].data[centre-radius:centre+radius+1,centre-radius:centre+radius+1]/1.0E6,np.arange(0,41/2.0,0.5),np.arange(0,41/2.0,0.5))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculate tiles\n", "As fitting the whole map would be too computationally expensive, I split based on HEALPix pixels. For MIPS, the optimum order is 11. So that I don't have to read the master prior based on the whole map into memory each time (which requires a lot more memory) I also create another layer of HEALPix pixels based at the lower order of 6." ] }, { "cell_type": "code", "execution_count": 16, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "----- There are 5442 tiles required for input catalogue and 13 large tiles\n" ] }, { "ename": "SystemExit", "evalue": "", "output_type": "error", "traceback": [ "An exception has occurred, use %tb to see the full traceback.\n", "\u001b[0;31mSystemExit\u001b[0m\n" ] } ], "source": [ "import pickle\n", "#from moc, get healpix pixels at a given order\n", "from xidplus import moc_routines\n", "order=11\n", "tiles=moc_routines.get_HEALPix_pixels(order,prior_MIPS.sra,prior_MIPS.sdec,unique=True)\n", "order_large=6\n", "tiles_large=moc_routines.get_HEALPix_pixels(order_large,prior_MIPS.sra,prior_MIPS.sdec,unique=True)\n", "print('----- There are '+str(len(tiles))+' tiles required for input catalogue and '+str(len(tiles_large))+' large tiles')\n", "output_folder='./data/'\n", "outfile=output_folder+'Master_prior.pkl'\n", "with open(outfile, 'wb') as f:\n", " pickle.dump({'priors':[prior_MIPS],'tiles':tiles,'order':order,'version':xidplus.io.git_version()},f)\n", "outfile=output_folder+'Tiles.pkl'\n", "with open(outfile, 'wb') as f:\n", " pickle.dump({'tiles':tiles,'order':order,'tiles_large':tiles_large,'order_large':order_large,'version':xidplus.io.git_version()},f)\n", "raise SystemExit()" ] }, { "cell_type": "code", "execution_count": 17, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "128252" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "prior_MIPS.nsrc" ] }, { "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.6" } }, "nbformat": 4, "nbformat_minor": 2 }