{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# XID+PACS GAMA-09 Prior"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "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, which provide 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\n",
    "from astropy.table import Table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This notebook uses all the raw data from the XID+PACS 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."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "Sel_func=pymoc.MOC()\n",
    "Sel_func.read('../../dmu4/dmu4_sm_GAMA-09/data/holes_GAMA-09_ukidss_k_O16_20180417.fits')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in CIGALE predictions catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cigale=Table.read('../../dmu28/dmu28_GAMA-09/data/GAMA9_Ldust_prediction_results.fits')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "ename": "KeyError",
     "evalue": "'id'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyError\u001b[0m                                  Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-5-cc26fb1fff49>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mcigale\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'id'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m'help_id'\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m~/anaconda3/lib/python3.6/site-packages/astropy/table/table.py\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, item)\u001b[0m\n\u001b[1;32m   1223\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mitem\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1224\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1225\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcolumns\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1226\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mint\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteger\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1227\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mRow\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mitem\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.6/site-packages/astropy/table/table.py\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, item)\u001b[0m\n\u001b[1;32m    105\u001b[0m         \"\"\"\n\u001b[1;32m    106\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 107\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mOrderedDict\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mitem\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    108\u001b[0m         \u001b[0;32melif\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mint\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minteger\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    109\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvalues\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mitem\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyError\u001b[0m: 'id'"
     ]
    }
   ],
   "source": [
    "cigale['id'].name = 'help_id'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=1815905</i>\n",
       "<table id=\"table140684224522112\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>bayes.dust.luminosity</th><th>bayes.dust.luminosity_err</th><th>best.omegacam_u</th><th>best.megacam_u</th><th>best.omegacam_g</th><th>best.suprime_g</th><th>best.decam_g</th><th>best.megacam_g</th><th>best.gpc1_g</th><th>best.suprime_r</th><th>best.gpc1_r</th><th>best.megacam_r</th><th>best.omegacam_r</th><th>best.decam_r</th><th>best.gpc1_i</th><th>best.megacam_i</th><th>best.omegacam_i</th><th>best.suprime_i</th><th>best.gpc1_z</th><th>best.vista_z</th><th>best.megacam_z</th><th>best.suprime_z</th><th>best.decam_z</th><th>best.gpc1_y</th><th>best.suprime_y</th><th>best.vista_y</th><th>best.ukidss_y</th><th>best.ukidss_j</th><th>best.vista_j</th><th>best.ukidss_h</th><th>best.vista_h</th><th>best.vista_ks</th><th>best.ukidss_k</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th></tr></thead>\n",
       "<thead><tr><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>1.7756811583627856e+39</td><td>1.8560994213383705e+39</td><td>0.0006849655648856183</td><td>0.0008213525215570016</td><td>0.0014418320784829225</td><td>0.0014851635254180297</td><td>0.0014918097398899804</td><td>0.0015248420785701489</td><td>0.0015223639970083685</td><td>0.0029236244063870925</td><td>0.0029323451806340392</td><td>0.0029942044550754767</td><td>0.003071945806420398</td><td>0.003209946563396643</td><td>0.00570727344151689</td><td>0.005874273794688017</td><td>0.006124423478050316</td><td>0.006517263751715815</td><td>0.012043668064701779</td><td>0.012387236998768944</td><td>0.012730939515909855</td><td>0.012885496989621615</td><td>0.013894760724993507</td><td>0.015677503791890205</td><td>0.015995150517214757</td><td>0.01681551604212792</td><td>0.017040868125072547</td><td>0.024214296725186508</td><td>0.02429380670766754</td><td>0.03624442746030166</td><td>0.03653520866519151</td><td>0.053322189508958964</td><td>0.05449461857253117</td></tr>\n",
       "<tr><td>2.7021265703371776e+37</td><td>2.4317979538966422e+36</td><td>0.05731104329672883</td><td>0.056223578961906474</td><td>0.06033796980008266</td><td>0.06082318037150582</td><td>0.06095714255843882</td><td>0.060962559918838935</td><td>0.06123168837954967</td><td>0.05832850068459655</td><td>0.05801617455383524</td><td>0.0574112860668829</td><td>0.05741336122906761</td><td>0.056988475834606624</td><td>0.05282675100407725</td><td>0.053128606958725344</td><td>0.05306034143787001</td><td>0.05267269867326263</td><td>0.051136801552879335</td><td>0.051096163964210975</td><td>0.051204577645507635</td><td>0.05120398359763581</td><td>0.05122226385822606</td><td>0.05100542972748624</td><td>0.05089417995814178</td><td>0.051669357982216245</td><td>0.052035689812129096</td><td>0.053585319826126916</td><td>0.05354605081471752</td><td>0.05234953875464276</td><td>0.052339845325697625</td><td>0.052655239140812916</td><td>0.0514389634370146</td></tr>\n",
       "<tr><td>4.570194861817776e+37</td><td>2.8856447896428874e+37</td><td>0.001562775024576193</td><td>0.0019241849445430818</td><td>0.004073465839980223</td><td>0.004389458937042463</td><td>0.004445338813423105</td><td>0.0046145338160607985</td><td>0.004617749979019916</td><td>0.013338348306415565</td><td>0.013418082449615516</td><td>0.013756858071279767</td><td>0.014021819449388108</td><td>0.01455585883157852</td><td>0.021016747199825565</td><td>0.02094968585405312</td><td>0.02126166115867816</td><td>0.021866551844849095</td><td>0.027435222606925637</td><td>0.02793834579131693</td><td>0.02841994636277121</td><td>0.02873496116158678</td><td>0.030069647897269628</td><td>0.032687682671318524</td><td>0.03361467379577717</td><td>0.03654262073705172</td><td>0.03715744099421227</td><td>0.05233128843384468</td><td>0.0525172784307798</td><td>0.06935969566164578</td><td>0.06970494617425149</td><td>0.09128148854208384</td><td>0.09383798173773936</td></tr>\n",
       "<tr><td>1.5495572578977246e+39</td><td>2.0426412306543896e+39</td><td>0.0007113458085920682</td><td>0.0008620998241307844</td><td>0.0014636804186210113</td><td>0.0015183653964489588</td><td>0.001528907940829464</td><td>0.0015650614718851665</td><td>0.001565504843687772</td><td>0.003486666574064081</td><td>0.0034925073419556116</td><td>0.003569124910119692</td><td>0.003707632994832172</td><td>0.003974755257606487</td><td>0.009555373494452808</td><td>0.009514978147316</td><td>0.009834455087319656</td><td>0.01053562651452195</td><td>0.017090714828527045</td><td>0.01772392654238113</td><td>0.017995607058220953</td><td>0.018406453207120137</td><td>0.01938260611424127</td><td>0.021299714542906932</td><td>0.022002392212932025</td><td>0.024543109914095355</td><td>0.025282734427615784</td><td>0.03732825099814391</td><td>0.037421056282661354</td><td>0.061449745967043015</td><td>0.06211407956664785</td><td>0.08968285075940517</td><td>0.09202998934088379</td></tr>\n",
       "<tr><td>1.9836831174912827e+39</td><td>1.1247053696095027e+38</td><td>6.298024535939316</td><td>6.138427217228513</td><td>6.478874044154541</td><td>6.499129387376567</td><td>6.5054687822726445</td><td>6.517273282028814</td><td>6.524689489260027</td><td>5.877294815665693</td><td>5.8475216229471565</td><td>5.792011289163387</td><td>5.782899379492132</td><td>5.7419449044740025</td><td>5.257985530146233</td><td>5.269362159585479</td><td>5.26062085048225</td><td>5.222712494633159</td><td>5.05455001321386</td><td>5.061177790248698</td><td>5.070187866004547</td><td>5.080311092127211</td><td>5.070331782010034</td><td>5.02012905273381</td><td>5.034770027035303</td><td>5.152091010483606</td><td>5.171181713596566</td><td>5.150449480333911</td><td>5.154095789248365</td><td>5.051034687088109</td><td>5.0549338616517465</td><td>4.9150757859246115</td><td>4.793075837926555</td></tr>\n",
       "<tr><td>9.008484458912574e+35</td><td>6.91995525495998e+35</td><td>0.07220251585052244</td><td>0.08880290109439372</td><td>0.20647972445855084</td><td>0.2105730067769274</td><td>0.21132906029022494</td><td>0.21734403031818966</td><td>0.21624840155682482</td><td>0.3434780306298811</td><td>0.344869563689463</td><td>0.3492399029340926</td><td>0.3541050196162975</td><td>0.3623248632048691</td><td>0.44410266529498543</td><td>0.4446860018573015</td><td>0.4487933614356649</td><td>0.4556164993451814</td><td>0.5359516185086711</td><td>0.546621011028694</td><td>0.5575971562099775</td><td>0.5552240278410377</td><td>0.5872793521796774</td><td>0.6383340563539209</td><td>0.65371752934474</td><td>0.687683157958246</td><td>0.6968285431579705</td><td>0.7690133643775512</td><td>0.7701625983507598</td><td>0.9232093024469973</td><td>0.9285849692123795</td><td>0.8145639697519402</td><td>0.7928518585568228</td></tr>\n",
       "<tr><td>3.601881487628594e+37</td><td>4.012016413741961e+37</td><td>0.008997100841526664</td><td>0.009501296483101336</td><td>0.012026686590607742</td><td>0.012569315929390514</td><td>0.012666843647848283</td><td>0.012809705273908052</td><td>0.012888039869490579</td><td>0.022415620052108315</td><td>0.022408800900497558</td><td>0.022587861833334038</td><td>0.022688343131727302</td><td>0.02287622471601989</td><td>0.024927102045823362</td><td>0.024932846246666984</td><td>0.025059341539773586</td><td>0.025231373467572496</td><td>0.026671696021761917</td><td>0.026791201171422318</td><td>0.027022632545495754</td><td>0.027012946429377915</td><td>0.027468761062184953</td><td>0.02825663518425728</td><td>0.028473775075190704</td><td>0.0294764504635975</td><td>0.02964108445294334</td><td>0.03419187268274122</td><td>0.03426414598595954</td><td>0.03691896256737512</td><td>0.036993528658823896</td><td>0.04139509485770477</td><td>0.04205066456324092</td></tr>\n",
       "<tr><td>6.326384780716312e+36</td><td>3.877157504172003e+36</td><td>0.029649954914199736</td><td>0.030813125127983475</td><td>0.05216607238625812</td><td>0.053551778604082216</td><td>0.05386316172944426</td><td>0.05454625387995087</td><td>0.055069418037298166</td><td>0.07015553214799677</td><td>0.06993405942528258</td><td>0.06994755202859998</td><td>0.0705539112371537</td><td>0.07142931833612215</td><td>0.07735553050171329</td><td>0.07731336357104028</td><td>0.07754269658589161</td><td>0.07786567842405112</td><td>0.08303982768621315</td><td>0.08378339577544028</td><td>0.08431096835942949</td><td>0.0850653096013099</td><td>0.086231003665275</td><td>0.08801115072636899</td><td>0.08904786259734587</td><td>0.09289337581208525</td><td>0.09380086781264434</td><td>0.10296717810120207</td><td>0.10312088885743537</td><td>0.11041042678028722</td><td>0.11066694911260547</td><td>0.1147225158896037</td><td>0.1119708009075359</td></tr>\n",
       "<tr><td>2.049659370130672e+38</td><td>1.7284443198752115e+38</td><td>0.06889481876409997</td><td>0.10840735518362506</td><td>0.45640410374469503</td><td>0.4788100798894797</td><td>0.48325877729157957</td><td>0.5036470484027452</td><td>0.5028229625625584</td><td>1.2877909267241812</td><td>1.2945004629626646</td><td>1.327779431278502</td><td>1.3701111628827567</td><td>1.4443441319317145</td><td>2.2948907637329037</td><td>2.3072969361333215</td><td>2.3491644358825132</td><td>2.4320992623691766</td><td>3.4235506956418336</td><td>3.54382172986511</td><td>3.6634625994652437</td><td>3.6647486508378697</td><td>4.032076107166555</td><td>4.686821407852258</td><td>4.900439463354127</td><td>5.54922765089344</td><td>5.7593400330451185</td><td>8.107505873555839</td><td>8.139940600050963</td><td>12.47782923107716</td><td>12.600650806937768</td><td>13.807142344135947</td><td>13.750692964949224</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>1.100216810929157e+37</td><td>5.896175988586009e+36</td><td>0.018592579601624082</td><td>0.020919925442333502</td><td>0.05001991326381216</td><td>0.052150489514316285</td><td>0.05260887209600786</td><td>0.05355893539008891</td><td>0.05420767587979453</td><td>0.0916988524545678</td><td>0.09158483163979131</td><td>0.09263653124022615</td><td>0.09437655045604403</td><td>0.09680966358839119</td><td>0.12091590915083761</td><td>0.12138962051225902</td><td>0.12226870965585715</td><td>0.12381957353135514</td><td>0.1459635259385529</td><td>0.14886209575768367</td><td>0.15086903584476083</td><td>0.15163037418287462</td><td>0.15714981016777657</td><td>0.1673550028740296</td><td>0.17058026370137358</td><td>0.18365713436150233</td><td>0.18691925893670955</td><td>0.22420426512402597</td><td>0.22470514203246972</td><td>0.2628864665632054</td><td>0.2640164030172578</td><td>0.2833373605386113</td><td>0.27644285679340813</td></tr>\n",
       "<tr><td>1.4206447502514093e+38</td><td>1.5785427761039392e+38</td><td>0.0010961636735176343</td><td>0.0012785054320938445</td><td>0.0025510925605829135</td><td>0.002678409217417819</td><td>0.0026980935737772375</td><td>0.0028340358800738145</td><td>0.0027905933139637617</td><td>0.008951242957455124</td><td>0.009070721793538732</td><td>0.00939015587071178</td><td>0.009654242875984131</td><td>0.010176472359692068</td><td>0.014786159688862512</td><td>0.014900246743430467</td><td>0.015126310254990867</td><td>0.015511031664951639</td><td>0.02020876996193327</td><td>0.02059542910670686</td><td>0.02095925872160439</td><td>0.021068221780402678</td><td>0.02216930093097443</td><td>0.02396711096546658</td><td>0.02458632549652058</td><td>0.026340242989286777</td><td>0.027029454515670386</td><td>0.03782988953354445</td><td>0.03799088458818508</td><td>0.053407596012437365</td><td>0.05361325984174868</td><td>0.06795484121987896</td><td>0.07058803064886479</td></tr>\n",
       "<tr><td>5.695247881763603e+37</td><td>1.1818187020180186e+37</td><td>0.48377129277079023</td><td>0.5895290888085368</td><td>1.5996456578518246</td><td>1.6487132364207824</td><td>1.6583874164923507</td><td>1.6942647945307476</td><td>1.7048917135447643</td><td>2.985469298961774</td><td>2.9958774821464162</td><td>3.0629960548625617</td><td>3.111944308697739</td><td>3.2109182601390183</td><td>4.133843632775286</td><td>4.124373713292741</td><td>4.1638873651925845</td><td>4.242162829051041</td><td>5.083748301049679</td><td>5.17485401848794</td><td>5.301499351005254</td><td>5.322851498098887</td><td>5.649763338515066</td><td>6.1720861140939824</td><td>6.382154860708415</td><td>6.813723990009727</td><td>6.910248457317847</td><td>8.028070665860378</td><td>8.046543792597635</td><td>9.858762054365634</td><td>9.919408009897797</td><td>9.514860119291388</td><td>9.317719358284283</td></tr>\n",
       "<tr><td>1.3552394175194585e+37</td><td>3.328740922492263e+36</td><td>0.007603298327210055</td><td>0.00931318609612951</td><td>0.029686863716538225</td><td>0.031365911604301945</td><td>0.031726746129917466</td><td>0.03243265915722482</td><td>0.03286298008690914</td><td>0.07074757234791672</td><td>0.07095506386509463</td><td>0.07227177403644384</td><td>0.07415352555725889</td><td>0.07713097429627166</td><td>0.10808333450625882</td><td>0.10896631669979505</td><td>0.1102114205678424</td><td>0.11272866727972175</td><td>0.14484780175892786</td><td>0.14873884753221064</td><td>0.15135491876416204</td><td>0.15164920343420854</td><td>0.16042368348909616</td><td>0.17659527708167252</td><td>0.18160096937317266</td><td>0.2015675282393538</td><td>0.20628974455350854</td><td>0.26916066593180815</td><td>0.2698963099795379</td><td>0.3448327226525966</td><td>0.3470370919408714</td><td>0.3922769597954111</td><td>0.3849792090702955</td></tr>\n",
       "<tr><td>1.328833590853375e+38</td><td>1.1914099343470326e+38</td><td>0.0006022683666690246</td><td>0.0007242911562616252</td><td>0.001582807793070364</td><td>0.0016451138443243443</td><td>0.0016579900950068388</td><td>0.001705854708666505</td><td>0.0017070424271012062</td><td>0.005980817140453547</td><td>0.006023205326675074</td><td>0.006317180214080813</td><td>0.006627180317476183</td><td>0.0071909695913864756</td><td>0.01208658706142171</td><td>0.012174052162286491</td><td>0.012333110000705265</td><td>0.012684823246904996</td><td>0.01744860880730458</td><td>0.017961729475200243</td><td>0.018233019958200668</td><td>0.018622847671863365</td><td>0.019668656125166325</td><td>0.021527603057021054</td><td>0.02226721239711716</td><td>0.02427949412047337</td><td>0.024648939932823292</td><td>0.03595833781293178</td><td>0.03622647541592231</td><td>0.05745830046414026</td><td>0.057764606777718296</td><td>0.07550539804206177</td><td>0.07825585808166652</td></tr>\n",
       "<tr><td>2.3710471636398194e+37</td><td>1.0020287258827178e+37</td><td>0.002592997456266232</td><td>0.003112399018737056</td><td>0.0061023112445072894</td><td>0.006503126451388559</td><td>0.006571282293492896</td><td>0.006758427208956666</td><td>0.00679344155648139</td><td>0.017205281656161787</td><td>0.017338333704537877</td><td>0.017703077481458575</td><td>0.017989121147900026</td><td>0.018643610360086025</td><td>0.025634580574466704</td><td>0.025477599946148805</td><td>0.025786242395717848</td><td>0.02648958083486864</td><td>0.03193053221901246</td><td>0.0325488395352815</td><td>0.03298686193119809</td><td>0.033292811743212365</td><td>0.03471680866000769</td><td>0.03720493346955144</td><td>0.038189529182115003</td><td>0.041053536242836114</td><td>0.04173014398620597</td><td>0.05648998465600138</td><td>0.056675709894982954</td><td>0.07076205068910989</td><td>0.0710360873721939</td><td>0.09006982246535605</td><td>0.09165088082958911</td></tr>\n",
       "<tr><td>1.5472336547762588e+38</td><td>1.5288694427606837e+38</td><td>0.00369256077195507</td><td>0.004642622314752422</td><td>0.011904422010889202</td><td>0.012770184106767188</td><td>0.012949281781316388</td><td>0.013135596133908856</td><td>0.013426204662944141</td><td>0.03466129412074146</td><td>0.034773130634512144</td><td>0.03570198564295854</td><td>0.03653541436520804</td><td>0.03809949105152831</td><td>0.054732980778248605</td><td>0.054723443558180596</td><td>0.055378889911757684</td><td>0.05681130514990006</td><td>0.07099196854384845</td><td>0.07228140848763041</td><td>0.07409012059666893</td><td>0.07422049580136385</td><td>0.07918996361429657</td><td>0.08674601452280511</td><td>0.08994579167580108</td><td>0.09605466456436049</td><td>0.09751370092086376</td><td>0.14024188830228226</td><td>0.14075950776964433</td><td>0.18575023545761693</td><td>0.1867613064242339</td><td>0.2518923327391029</td><td>0.2537738344197868</td></tr>\n",
       "<tr><td>8.403366521447583e+37</td><td>6.959359922574835e+37</td><td>0.0030236094392045326</td><td>0.0038458442562886394</td><td>0.010759143167665674</td><td>0.011432802980803155</td><td>0.01158362627669089</td><td>0.011921752767770497</td><td>0.01205361321269079</td><td>0.031881735652508375</td><td>0.03203586715595845</td><td>0.032627801545245506</td><td>0.0335940276088581</td><td>0.03509817005949716</td><td>0.051943536930588556</td><td>0.05194372536754738</td><td>0.05257599852768525</td><td>0.05402226436659846</td><td>0.06813881840738426</td><td>0.06976206444064686</td><td>0.07163647237235655</td><td>0.07175467390295569</td><td>0.07656731915947898</td><td>0.08560022622515245</td><td>0.0877416765505191</td><td>0.09402372330222404</td><td>0.09574400438394556</td><td>0.14091797871203673</td><td>0.1412710275967499</td><td>0.1868194510991229</td><td>0.18776910143114114</td><td>0.2532260109821104</td><td>0.25316678423230976</td></tr>\n",
       "<tr><td>1.617886580946589e+38</td><td>1.6853216553185113e+38</td><td>0.0020014530995138484</td><td>0.0023247639869776143</td><td>0.003766347028174967</td><td>0.003846383383959872</td><td>0.0038621873451095546</td><td>0.003945623940047657</td><td>0.003923687642800045</td><td>0.008099955099837623</td><td>0.008178518818648683</td><td>0.00843678999758594</td><td>0.008598102208437536</td><td>0.008930433514208254</td><td>0.011199635714184381</td><td>0.01124408068727464</td><td>0.011319347378660515</td><td>0.011463854512217181</td><td>0.013277923708261245</td><td>0.013451052810052594</td><td>0.013595504147970151</td><td>0.013670617796268986</td><td>0.01409491492982668</td><td>0.014756253205193153</td><td>0.015043323503812208</td><td>0.015756623953759392</td><td>0.015915399705521464</td><td>0.020170506197416213</td><td>0.02027972413428115</td><td>0.03567201601300925</td><td>0.03591339587577191</td><td>0.04546733192306832</td><td>0.04635943708894945</td></tr>\n",
       "<tr><td>1.2800960529174906e+38</td><td>1.4743806275086965e+38</td><td>0.0023776253388504633</td><td>0.00249984743414064</td><td>0.003201621302028554</td><td>0.0032462057265160876</td><td>0.003254030974748186</td><td>0.003276313896207675</td><td>0.0032812111105748254</td><td>0.006168121866117798</td><td>0.006208274967756541</td><td>0.006450901975866853</td><td>0.006617406116026633</td><td>0.006977490311363135</td><td>0.009475158662437172</td><td>0.009531143754364514</td><td>0.009580341132715331</td><td>0.009712410655140884</td><td>0.01149972460620663</td><td>0.01170134362975242</td><td>0.011804325853170795</td><td>0.011967755272609328</td><td>0.012345217290329023</td><td>0.012997286555072212</td><td>0.013170925173820623</td><td>0.013868981688387348</td><td>0.01403587348827635</td><td>0.017668205425703545</td><td>0.01775597691467449</td><td>0.024615802059289694</td><td>0.024687226493178557</td><td>0.0285686824681277</td><td>0.029216990406834058</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=1815905>\n",
       "bayes.dust.luminosity  bayes.dust.luminosity_err ...    best.ukidss_k    \n",
       "                                                 ...         mJy         \n",
       "       float64                  float64          ...       float64       \n",
       "---------------------- ------------------------- ... --------------------\n",
       "1.7756811583627856e+39    1.8560994213383705e+39 ...  0.05449461857253117\n",
       "2.7021265703371776e+37    2.4317979538966422e+36 ...   0.0514389634370146\n",
       " 4.570194861817776e+37    2.8856447896428874e+37 ...  0.09383798173773936\n",
       "1.5495572578977246e+39    2.0426412306543896e+39 ...  0.09202998934088379\n",
       "1.9836831174912827e+39    1.1247053696095027e+38 ...    4.793075837926555\n",
       " 9.008484458912574e+35      6.91995525495998e+35 ...   0.7928518585568228\n",
       " 3.601881487628594e+37     4.012016413741961e+37 ...  0.04205066456324092\n",
       " 6.326384780716312e+36     3.877157504172003e+36 ...   0.1119708009075359\n",
       " 2.049659370130672e+38    1.7284443198752115e+38 ...   13.750692964949224\n",
       "                   ...                       ... ...                  ...\n",
       " 1.100216810929157e+37     5.896175988586009e+36 ...  0.27644285679340813\n",
       "1.4206447502514093e+38    1.5785427761039392e+38 ...  0.07058803064886479\n",
       " 5.695247881763603e+37    1.1818187020180186e+37 ...    9.317719358284283\n",
       "1.3552394175194585e+37     3.328740922492263e+36 ...   0.3849792090702955\n",
       " 1.328833590853375e+38    1.1914099343470326e+38 ...  0.07825585808166652\n",
       "2.3710471636398194e+37    1.0020287258827178e+37 ...  0.09165088082958911\n",
       "1.5472336547762588e+38    1.5288694427606837e+38 ...   0.2537738344197868\n",
       " 8.403366521447583e+37     6.959359922574835e+37 ...  0.25316678423230976\n",
       " 1.617886580946589e+38    1.6853216553185113e+38 ...  0.04635943708894945\n",
       "1.2800960529174906e+38    1.4743806275086965e+38 ... 0.029216990406834058"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cigale"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in photoz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "photoz=Table.read('../../dmu24/dmu24_GAMA-09/data/master_catalogue_gama-09_20171206_photoz_20180213_r_optimised.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "photoz"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Join CIGALE and photoz tables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "prior=join(cigale,photoz,keys='help_id')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.cosmology import Planck15 as cosmo\n",
    "from astropy import units as u\n",
    "f_pred=prior['bayes.dust.luminosity']/(4*np.pi*cosmo.luminosity_distance(prior['z1_median']).to(u.cm))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "prior=prior[np.isfinite(f_pred.value)][np.log10(f_pred.value[np.isfinite(f_pred.value)])>8.5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "prior['DEC'].name='Dec'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in Maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "PACS_100='../../dmu18/dmu18_HELP-PACS-maps/data/GAMA-09_PACS100_v0.9.fits'\n",
    "PACS_160='../../dmu18/dmu18_HELP-PACS-maps/data/GAMA-09_PACS160_v0.9.fits'\n",
    "#output folder\n",
    "output_folder='./'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "from astropy import wcs\n",
    "\n",
    "#-----100-------------\n",
    "hdulist = fits.open(PACS_100)\n",
    "im100phdu=hdulist[0].header\n",
    "im100hdu=hdulist[1].header\n",
    "im100=hdulist[1].data\n",
    "w_100 = wcs.WCS(hdulist[1].header)\n",
    "pixsize100=3600.0*np.abs(hdulist[1].header['CDELT1']) #pixel size (in arcseconds)\n",
    "nim100=hdulist[2].data\n",
    "hdulist.close()\n",
    "\n",
    "#-----160-------------\n",
    "hdulist = fits.open(PACS_160)\n",
    "im160phdu=hdulist[0].header\n",
    "im160hdu=hdulist[1].header\n",
    "\n",
    "im160=hdulist[1].data #convert to mJy\n",
    "w_160 = wcs.WCS(hdulist[1].header)\n",
    "pixsize160=3600.0*np.abs(hdulist[1].header['CDELT1']) #pixel size (in arcseconds)\n",
    "nim160=hdulist[2].data\n",
    "hdulist.close()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in PSF"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "pacs100_psf=fits.open('../../dmu18/dmu18_GAMA-09/dmu18_PACS_100_PSF_GAMA-09_20190131.fits')\n",
    "pacs160_psf=fits.open('../../dmu18/dmu18_GAMA-09/dmu18_PACS_160_PSF_GAMA-09_20190131.fits')\n",
    "\n",
    "centre100=np.long((pacs100_psf[1].header['NAXIS1']-1)/2)\n",
    "radius100=15\n",
    "centre160=np.long((pacs160_psf[1].header['NAXIS1']-1)/2)\n",
    "radius160=25\n",
    "\n",
    "pind100=np.arange(0,radius100+1+radius100,1)*3600*np.abs(pacs100_psf[1].header['CDELT1'])/pixsize100 #get 100 scale in terms of pixel scale of map\n",
    "pind160=np.arange(0,radius160+1+radius160,1)*3600*np.abs(pacs160_psf[1].header['CDELT1'])/pixsize160 #get 160 scale in terms of pixel scale of map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[  0.           0.33333333   0.66666667   1.           1.33333333\n",
      "   1.66666667   2.           2.33333333   2.66666667   3.           3.33333333\n",
      "   3.66666667   4.           4.33333333   4.66666667   5.           5.33333333\n",
      "   5.66666667   6.           6.33333333   6.66666667   7.           7.33333333\n",
      "   7.66666667   8.           8.33333333   8.66666667   9.           9.33333333\n",
      "   9.66666667  10.        ]\n"
     ]
    }
   ],
   "source": [
    "print(pind100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x11fd211d0>"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABGMAAAItCAYAAACdL7ciAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3W2sZVd95/nfvec+nPtcD35iFIKDjJbISLRbJoEJ2C5o\nMDHdGU+TmcyoBd1pZmgYuQeYJmMGcCmaBiaCJs6MiRKiIrQT1FH3xMT9AsUDUWMbwyQNoUFqC/cC\nowYj4udy3bpV9/lhXtxT5vjevfe/9lpnr7PvWd+PVNI9e9Xae+21H84663Fsb29PAAAAAAAASGN8\n2AkAAAAAAADICZUxAAAAAAAACVEZAwAAAAAAkBCVMQAAAAAAAAlRGQMAAAAAAJDQxLATAABAW42N\njSVbcnBvb28s1bEAAADaKmX5SxpeGYyeMQAAAAAAAAlRGQMAAAAAAJAQw5QAACgxNsbIIQAAgJRy\nKX/RMwYAAAAAACAhesYAAFAil5YZAACAtsil/EXPGAAAAAAAgIToGQMAQIlcWmYAAADaIpfyFz1j\nAAAAAAAAEqJnDAAAJcbHabMAAABIKZfyVx5nCQAAAAAA0BL0jAEAoEQuY5YBAADaIpfyFz1jAAAA\nAAAAEqJnDAAAJXJpmQEAAGiLXMpf9IwBAAAAAABIiMoYAAAAAACAhBimBABAiVy6yQIAALRFLuUv\nesYAAAAAAAAkRM8YAABK5NIyAwAA0Ba5lL/oGQMAAAAAAJAQPWMAACiRS8sMAABAW+RS/qJnDAAA\nAAAAQEL0jAEAoMT4OG0WAAAAKeVS/srjLAEAAAAAAFqCnjEAAJTIZcwyAABAW+RS/qJnDAAAAAAA\nQEL0jAEAoEQuLTMAAABtkUv5i54xAAAAAAAACdEzBgCAErm0zAAAALRFLuUvesYAAAAAAAAkRM8Y\nAABK5NIyAwAA0Ba5lL/oGQMAAAAAAFDCOfca59yDB7b9A+fcX/Z9fpdz7q+dc3/lnPt71j6pjAEA\nAAAAACjgnLtD0mcldfu2/W1J/6Oksd7nayS9V9LrJL1F0m8556ar9ktlDAAAJcbGxpL9AwAAQNry\n12WWwX4g6W2XPjjnTkr6PyW9v+///KKkr3vvN7z3y5Iek/Sqqp1SGQMAAAAAAFDAe/8FSVuS5Jzr\nSPpDSf9M0krff1uUtNz3eUXSUtV+mcAXAIAS4+O0WQAAAKTU8vLXDZJeIen3tT9s6eedc/+XpK9I\nWuj7fwuSzlXtiMoYAAAAAAAAg/f+G5L+S0lyzl0r6V9779/fmzPm4865rqRpSa+U9EjVvqiMAQCg\nBHO5AAAApHUUy1/e+yedc3dLelj708F8xHu/XhVnbG9vL0niAAA4aq688spkX5LPPPPM0St5AAAA\nDFjK8pc0vDIYPWMAAChxFFtmAAAAjrJcyl+tnhkHAAAAAABg1NAzBgCAErm0zAAAALRFLuUvesYA\nAAAAAAAkRM8YAABK5NIyAwAA0Ba5lL/oGQMAAAAAAJAQPWMAACiRS8sMAABAW+RS/qJnDAAAAAAA\nQEL0jAEAoMT4OG0WAAAAKeVS/srjLAEAAAAAAFqCyhgAAAAAAICEGKYEAECJXCaQAwAAaItcyl/0\njAEAAAAAAEiInjEAAJTIpWUGAACgLXIpf9EzBgAAAAAAICF6xgAAUCKXlhkAAIC2yKX8Rc8YAAAA\nAACAhOgZAwBAiVxaZgAAANoil/IXPWMAAAAAAAASomcMAAAlcmmZAQAAaItcyl/0jAEAAAAAAEiI\nnjEAAJQYH6fNAgAAIKVcyl95nCUAAAAAAEBL0DMGAIASuYxZBgAAaItcyl/0jAEAAAAAAEiInjEA\nAJTIpWUGAACgLXIpf9EzBgAAAAAAICEqYwAAAAAAABJimBIAACVy6SYLAADQFrmUv+gZAwAAAAAA\nkFCSnjGrq6t7ZWHdblfr6+sDOU5IDdr4eP36qE6nUztOVdrGx8e1u7tbGLazs1P7WGX7GrSQ/K7K\nu06nU3i+e3ult0+pkOsaktchqs5ncnJSW1tbhWEh+R2Sd22IU/ZeCDlOiJDjDPr6zMzMaG1trfY+\ny4Tc36nyu+q9MDs7q9XV1UPbFxYWkjSZ5NIyg9G1u7v7ogd5bGys1rNt/V+rzGGFx75nrGc0Nvxg\n+iYmJrS9vX15iTP2X/fYgw63DPo7YGpqSpubmy98jn2/DvveaTp/D4ZPT09rY2Mj6pj9ms7/2Pyp\nKsvXPfbBcqVVJop971ms3ykTE9U/0UN+58Q4+Ox2u93GC0e5lL+G3jMm9c3URrncbBbygefhEvIh\nrNJ3FJEPwGDxXRuH76dw5F0c8i8ceReH/GsOc8YAAFCCAgjaxjk3Lun3JP0tSRuS/ifv/WPDTRUA\nAIOTS/krj7MEAAAYDf+NpK73/r+S9L9L+u0hpwcAAAQI6hlDqwwAIAcM6UALvV7S/ytJ3vu/cs69\nuuo/j42NHbqPB9nimONQwqmpqWEn4cjqdrvDTsKRNjMzM+wkHFmzs7PDTsKRlvrZzaX8FTpM6YVW\nGefca7XfKnPb4JIFAACAAouSlvs+7zjnJrz3hbPK7u3tvWgyyqpFA8riV8ltAt+DE1lamMD3pw5O\nosoEvvWu36An9s9pAt+DiwEwgW89B59dKlUHJ7QyplarDAAAR1EuLTM4Us5LWuj7PF5WEQMAwFGU\nS/krtDKmVqtMt9utrMGj21h5N9/cuv9aNcFNakteT09PDzsJrcB7QZqfnx92ElphYWHB/k9APr4u\n6Vck/T+93sn/ccjpAQAAAUJ/+dZqlenv1nTQwW5jMUJq0EK6eYX8aK9KW6fTKe0uZ3WjKxLbde5y\nheR3Vd5NTExoe/vwbRTSzTHkuobkdYiq85mentbGxkZhWEh+h+RdG+KUvRcG3WW6TMhxBn195ufn\ndeHChdr7LBNyf6fK76r3wsLCglZWVgq3p5DLbP44Uu6T9Gbn3P8naUzSP676z0XPcZ1n+6gPlbHi\np/puHYamW5ZD8rZ/W+wwIau82/Qwo1ghQ2EGOcSw6e+3pt8NbWbde1beNz380lJ0bVL3VMml/BVa\nGUOrDAAAQGLe+11J7xl2OgAAQJzQypharTIAABxFuYxZBgAAaItcyl9BlTG0ygAAAAAAAIQZ3myp\nAAC0XC5jlgEAANoil/JXksqYVBM4hXRnShWn7UJu+CbiFIWnmpA4RBP3T1ketXkC35DJYUMmH2vz\nZG8habPu7aLwlO+fkGe8zdcIAAAAaIs8qpwAAAAAAABagmFKAACUGMWekAAAAG2WS/mLyhgAAIAR\ndXC4Y6fTSTr81ypQxw5ttOI3cfxUwzFDhvP2i50moIlr078t9tq1XWz+Fj2ndZ7dpvOv6fsr5j1l\nDXWPfbYs1jDvTqcTFR465cHlKro2VpoQhsoYAABKHPUfAwAAAEdNLuUv5owBAAAAAABIiJ4xAACU\nyGVpRQAAgLbIpfyVx1kCAAAAAAC0BD1jAAAokcuYZQAAgLbIpfxFzxgAAAAAAICE6BkDAECJXMYs\nAwAAtEUu5S8qYwAAAEZUUVfvlN2/9/b2KsNj02Ltf9hi0hebN7HxY38MFR1/YuKnPz12d3ej9t/0\ntbfSZ4Vb6QvZf508s67fUX52rLQXnXt/nNj3kpW3nU4nKn4Tz16/nZ2dqPgYHCpjAAAoQYEEAAAg\nrVzKX0kqY6zMLApvc21tSNpS3lAhtakhcaxa3yIhNc0h+R3S2hJyjULyzTpOSL6WSZV3Iay0FYW3\n/dmrK+TdGHo+g7yvqrT53T0KnHNXSfqWpDdL2pZ0j6Q9SY9Iut17v+uce5ekd/fCP+a9/+KQkgsA\nAIAS9IwBAKBEm8YsO+cmJf2BpLXeprsk3em9f9A59xlJtznn/lLSeyW9WlJX0tecc3/hvd8YSqIB\nAABqalP5q0l5nCUAAEffpyR9RtLf9D7fIOmh3t/3S3qTpF+U9HXv/Yb3flnSY5JelTqhAAAAqEZl\nDAAAJcbGxpL9q+Kc+3VJz3jvv9SfPO/9pXFhK5KWJC1KWu77P5e2AwAAHAkpy1/DnNKAYUoAALTf\nOyXtOefeJOl6SX8s6aq+8AVJ5ySd7/19cDsAAABahMoYAABKtGUCaO/9TZf+ds49KOk9kv6Fc+6U\n9/5BSbdKekDSNyR93DnXlTQt6ZXan9wXAADgSGhL+atpVMYAAHA0fUDSGefclKRHJd3rvd9xzt0t\n6WHtD0X+iPd+fZiJRLtZBV5rhbSmVwi0VvhregVAa1W72PypErLq4CDFTqDZdN7Fhsfmb2z+s/pg\ntar8sfJuZ2fnsraVse59Kzz23rRw7wyHc+41kj7hvT/lnLte0qcl7UjakPQPvfdP1V3RksoYAACO\nEO/9qb6PNxeEn5F0JlmCAAAARphz7g5J75B0sbfp/5b0v3jvv+Oce7ekDzrnPqmaK1pSGQMAQIlc\nllYEAABoixaWv34g6W2SPt/7/D9475/o/T0haV19K1pK2nDOXVrR8ptlO23dWQIAAAAAALSB9/4L\nkrb6Pj8hSc65X5L0TyX9jgJWtKRnDAAAJXKZQA4AAKAtjkL5yzn330v6iKS/671/xjlXe0VLKmMA\nAAAAAAAug3Pu7dqfqPeU9/5sb3PtFS2pjAEAoEQLxywDAACMtDaXv5xzHUl3S3pc0p855yTpIe/9\nb9Zd0TJJZUyqZQlDlvkKOfaglxPrdDra3t4e2P5CunWFxBl0fnc6ncJl50LyJiRtqfLNMsj7q84y\nfpeEPBNNPEdNL1daJeQaNBGnKA9Sfjm15Zk4Cl1VAQAAMLq89z+U9NrexxMl/6fWipb0jAEAoAQV\nQTjqiu7h/m3WPR77DFiVzrH7t+I3XbHf5DvCSnts403TeV+Uvv5to/5+jb1+RfFTNlQ1/W6wVOVP\nqkazUE2/N3Mw6u+HS9rb/wcAAAAAAGAE0TMGAIASubTMAAAAtEUu5S96xgAAAAAAACREzxgAAEq0\neTZ/AACAUZRL+SuPswQAAAAAAGgJesYAAFAilzHLAAAAbZFL+YueMQAAAAAAAAnRMwYAgBK5jFlG\nvqzWx9jWyb29vaj41jO4u7sbtf+Q9NWJ0+Q7xEpHbN7v7OxExT94/ImJCW1vb5eGHxR7b8bmj7V/\n69rGxu/Pq6I4sednPTux5xf77hh0z4j+/cVem9jnOvbes+LHPrttkEv5K4+zBAAAAAAAaAkqYwAA\nAAAAABJimBIAACVymUAOAACgLXIpfyWpjLHGJBaFh4wTCxkbGzvWeFDKxvalGi+XKu+q4nQ6ncLx\nsSHjHkPOJ+ShH/RxyvJACsvvVHnXhEGlo+n5AJqMUxSe8vqEPBNNfHnm8oUMAACAfNAzBgCAElQE\nAQAApJVL+Ys5YwAAAAAAABKiZwwAACVyWVoRAACgLXIpf1EZAwAAMKKsuaeanodq2PuPDY+NU/V/\nrW74TYdbYq+dde81fX5NpL/O8a1way7Aovj926z4VvqtH7uxz05s/lSlL2QexaM07CX23rXmjOx0\nOpXh1r2HwaEyBgCAEhQ+AAAA0sql/JVH/x8AAAAAAICWoGcMAAAlchmzDAAA0Ba5lL/yOEsAAAAA\nAICWoGcMAAAlchmzDAAA0Ba5lL/oGQMAAAAAAJAQPWMAACiRy5hlAACAtsil/EVlDAAAwIja3d2t\n3BbbFdyKH7v/vb29qPix+x9mV3nr2E3/WLHyJvbaWOkf9r3T9LUPeXb6tzX97MZqcv+dTqfROMN+\nrzUdbjkYv9PpaGdn50WfMRhJKmNiXzaXK+RLqekv+Vip0tf/gA0rzvT0tDY3Nw9tLypINqEN98L0\n9LTW1tYKw7a2tmrvL+QaheTD5ORk7TjW81p0viHvhVQF6SbeWUXhKVsKUuX3sH9wVMllzDIAAEBb\n5FL+yqP/DwAAAAAAQEtQGQMAAAAAAJAQc8YAAFAil26yAAAAbZFL+YueMQAAAAAAAAkF94xxzv0H\nSed7H/+z9/4fDyZJAAC0Qy4tMwAAAG2RS/krqDLGOdeVNOa9PzXY5AAAAAAAAIy20J4xf0vSrHPu\ny719fNh7/1dl/3l2drZyPfKFhYXAZIyO2dnZYSehFbgXpGPHjg07Ca1w4sSJYSdh6BYXF4edhFaY\nmZkZ2rFzaZnB6Nrd3b2sbWWspeWtZyQ2fG9vrzLcOhcrfqyYd4SV9qbz1rq2sXlbFN5/zKrfBoNQ\n5z5vQuy9Hbv/psMtTeZ/09/NsfsfdtmhifSnPqdh52EqoZUxq5I+Jemzkl4h6X7nnPPebxf+59XV\n0h0tLCxoZWXl0HbrC6JISJymv6Qvx+zsbGkehdyIqeLs7OwMNE7ZvZDqy7QN98KxY8d07ty5wrCt\nra3a+wu5RiH5MDk5WTtO1fN64sQJnT179tD2VPd2iEGnbXFxUefPnz+0PeQ9F6oN75+ZmRmtra0V\nbgcAAACOqtDKmO9Jesx7vyfpe8655yS9RNKPB5YyAACGLJeWGQAAgLbIpfwV2sT6Tkm/LUnOuf9C\n0qKkJwaVKAAAAAAAgFEV2jPmDyXd45z7mqQ9Se8sG6IEAMBRlUvLDAAAQFvkUv4Kqozx3m9K+gcD\nTgsAAAAAAMDIC+0ZAwDAyMulZQYAAKAtcil/JamMsVb/KAoPuQBtXk0p5fJ6qY4VsrrP9nb1aLbN\nzc1D20LOJ+S6WmkrEpK2qhWOjh07puXl5drxBqnppSYvsZ7Xovsr1UpCExP1X41NxJmamjq0LeWX\nU6rVlKznNdU9CQAAAKRCzxgAAErk0jKD0VXUaDDIRpvYSvLYZ6zpZ7Ro/7m8F6zzTNWgGSo2/bHn\nZ8W3nsOi+P3brPMbdnhs/sbkv9XQb6U9pCNB2bGKWOfW9DvGuvfa8N7L5T1LZQwAAEBLOedeI+kT\n3vtTzrnrJN2j/cUTHpF0u/c+XddbAAAwMFTGAABQItXQOKCIc+4OSe+QdLG36S5Jd3rvH3TOfUbS\nbZLuG1b6AABoQi7lrzzOEgAA4Oj5gaS39X2+QdJDvb/vl/Sm5CkCAAADQc8YAACAFvLef8E5d23f\npjHv/aXJBlYkLVn76Ha7h1oYZ2dnB5bGplkTeA9jgu9ut5v8mKNieno62bFGcfL3xcXFYSfhyDpK\n7702ClmoAjZyFQCAErlMIIcjo39+mAVJ56wI6+vrL/o8Ozur1dXVFz7HTlQZO5GlFd+aaNJaadAK\ntybSPJi+brd7KE+rVJ1/7CSeseFW3g96Atrp6WltbGy88LnpH3fWtQ+ZQHeQ6t7bi4uLOn/+/GXv\nf5Qn8K17bx987zGBb/W9d/D8JiYmXrTqbIqKmVzKXwxTAgAAOBq+7Zw71fv7VkkPDzEtAAAgAj1j\nAAAokUvLDI6MD0g645ybkvSopHuHnB4AAAYul/IXlTEAAAAt5b3/oaTX9v7+nqSb68Qv6g5fZ/hF\n7FADy7CHiliK0te/LXaoT0zctv9Yse4969rHir23Yoc5xbLyr+nr3/SzHzMMLuS569/W9PDKJt8L\nUvy1DxmClsvqRqlRGQMAQIm2/9gBAAAYNbmUv6jiAgAAAAAASIieMQAAlMilZQYAAKAtcil/JamM\niV0a7XKlGlcckt7YsYVNs8bFFgkZK7u1tVU73IoTcpxhxulfGu6gl73sZXr++edr77NMyH01NTVV\nO07I/WMti9e//OXlxikyOTlZO06bvwA6nU5QvFTvmJC8q3omQvcJAAAAtBk9YwAAKEFFEAAAQFq5\nlL+YMwYAAAAAACAhesYAAFAil5YZAACAtsil/EVlDAAAwIgqmi+qzhxS1nx8IXOG9YstcFvps+a3\nC5lvsH+fqeZFbEJs3oXkbf+2kLkH62j63ohNv7X/ontnkPNjxu4r9t6PCQ+J278t9rls+txD7o06\nQuYRbPO77CijMgYAgBIUPgAAANLKpfzFnDEAAAAAAAAJ0TMGAIASubTMAAAAtEUu5S96xgAAAAAA\nACREZQwAAAAAAEBCDFMCAKBELt1kAQAA2qKN5S/n3GskfcJ7f8o5d52keyTtSXpE0u3e+13n3Lsk\nvVvStqSPee+/WLVPKmMAAABGlLW0dezytlb8QS7FGyJ2+WJreeGml6Ct0vSxm7g36uRdyP77NbE0\n9yDDLUXx27SsurU8cpP3Z9P3dtPLmlvpb3r/TccfRc65OyS9Q9LF3qa7JN3pvX/QOfcZSbc55/5S\n0nslvVpSV9LXnHN/4b3fKNsvlTEAAJSgQAIAAJBWC8tfP5D0Nkmf732+QdJDvb/vl3SLpB1JX+9V\nvmw45x6T9CpJ3yzbaZLKmJ2dndrhIbXJnU6ndhyrVrdISNqarmGNZV2jIltbW7XjbGyUVgyWhq+t\nrdU+zurqau04IcfZ3t6uHce6F86dO1e4fWKi/uM6OTlZO06IkPvXSlvRvRDyYp6amqodJ+S9EBIn\npZC8C4kT8n6Mbb0EAAAAmuK9/4Jz7tq+TWPe+0uF3hVJS5IWJS33/Z9L20vRMwYAgBItbJkBAAAY\naUeg/NXfUrgg6Zyk872/D24v1e6mXAAAAAAAgPb4tnPuVO/vWyU9LOkbkm50znWdc0uSXqn9yX1L\n0TMGAIASR6BlBgAAYKQcgfLXBySdcc5NSXpU0r3e+x3n3N3ar5gZl/QR7/161U6ojAEAAAAAACjh\nvf+hpNf2/v6epJsL/s8ZSWcud59UxgAAUOIItMwAAACMlFzKX1TGAAAAjKiiAm3/NmsFtJAVMfuF\nrLDWL7ZAHpu+opU6+1d4i1lJz1opLvbcY/dv5U1IeP82K32xx49dqa/plf5Cno06q59aq8xa17/p\nH8Mxxw+JO8jzib032378g3nV6XRedL+HrGCMYlTGAABQIpeWGQAAgLbIpfzFakoAAAAAAAAJ0TMG\nAIASubTMAAAAtEUu5S96xgAAAAAAACREzxgAAErk0jIDAADQFrmUv+gZAwAAAAAAkFCSnjFbW1u1\nw0OWi5ucnKwdJ2RprpCautDavTpL2F0SspzZxsZGkjirq6u1w5eXl2sf59lnn60d5/nnn68dx7q3\ni0xMVD92jz/+eOH2ubm52scKiRNyXUOevZmZmcrwixcvHtoW8rx2u93acUKEvLOsONvb26HJqX2s\nIiH5HfLOqjrP6enpwvDp6enaxwEAAADagmFKAACUyKWbLPJlVdRa4SENQP2sZyz2GYzdf1H4oN4L\nsWkbH4/r4B57fOvaF4X3b7PiW5X71r1pxbcaPGKPb12fkPzvT7PVYBJ7/zT9/WflX9XxrXun6bRb\nx4+9N2Lfq7GsZzeFXMpfDFMCAAAAAABIiJ4xAACUyKVlBgAAoC1yKX/RMwYAAAAAACAhesYAAFAi\nl5YZAACAtsil/EXPGAAAAAAAgIToGQMAQIlcWmYAAADaIpfyFz1jAAAAAAAAEqJnDAAAJdrSMuOc\n60g6I8lJ2pP0Hknrku7pfX5E0u3e+13n3LskvVvStqSPee+/OJREoxX29vYua1tTrGcoNtw6l/Hx\n6nZHK37R8VO9F6zjWOcWy8qb3d3dyvCdnZ3KbUXh/TY3NyvDt7e3hxpu5c/ERPXPrE6nUzu8P08m\nJyejjm+JffZCnq068WM0/Qw3nTdNS/kdUWbYeZAKPWMAAGi/X5Ek7/3rJN0p6eOS7pJ0p/f+Rklj\nkm5zzl0j6b2SXifpLZJ+yzk3PZwkAwAAoAw9YwAAKNGWlhnv/b91zl3q4fIySeckvUnSQ71t90u6\nRdKOpK977zckbTjnHpP0KknfTJxkAACAIG0pfzWNnjEAABwB3vtt59wfSfq0pH8lacx7f6kv8Yqk\nJUmLkpb7ol3aDgAAgBZJ0jPGGhNaFB4yVs06TpGQ8bYhNXWhYwdDzmlra6t2nIsXL9aOs7q6WjvO\n2bNnK8OfffbZQ9uefvrp2sf58Y9/XDvO8vKy/Z8OCMnr6enqEQM/+tGPCrfPzMzUPtbc3FztOPPz\n87XjzM7O1o5z5ZVXVoafP3/+0DZrfHURK7+LWOOwizTxLikKDx3Ha43tLxJyrJDjWHFC9jkobWuZ\n8d7/I+fcByX9e0n9L4UF7feWOd/7++B2AACAI6Ft5a+m0DMGAICWc869wzn3od7HVUm7kv7aOXeq\nt+1WSQ9L+oakG51zXefckqRXan9yXwAAALQIc8YAAFCiRS0zfybpXzrnvippUtL7JT0q6Yxzbqr3\n973e+x3n3N3ar5gZl/QR7/36sBINAABQV4vKX42iMgYAgJbz3l+U9GsFQTcX/N8z2l8GGwAAAC11\nWZUxzrnXSPqE9/6Uc+46SfdI2tN+1+fbvffDG9APAACAIFbrozUfVsjcV4MUu39rbqyi/ffniXX8\nqrn/hp23ofMZXmLN57W5uVm5zZp3z5qbcHt7uzJ8Y2MjKr4VbrHmrLPmpysKt86p38RE9c886/6L\nvT8sMfFjjz3sXhch751Bhofk37DzbFSZc8Y45+6Q9FlJ3d6muyTd6b2/UdKYpNuaSx4AAMMzNjaW\n7B8AAADSlr+GWQa7nAl8fyDpbX2fb5D0UO/v+yW9adCJAgAAAAAAGFXmMCXv/Recc9f2bRrz3l/q\n27Qiacnax/Hjxyu7yllL3OYgZAniQWrLNbjxxhuHnYShO3369LCT0Ao33XTTsJMwdMN+L7TFwsKC\n/Z8aErJkOQAAAMLlUv4KmcC3f4DogqRzVoTnn3++NOzKK6/UM888c2h7yFg2a2xkEWu8ZpFB3xxz\nc3O6ePFiYZg1nnZQcS5cuFA7jjWWt8jZs2dLw2688UY9/PDDh7Y//fTTtY/z4x//uHac5eXl2nFC\n8rpqDPHp06f10Y9+tDBsZmam9rFCfszPz8/XjjM7O1s7TlUF4E033aSvfvWrh7YfP3689nFOnDhR\nO05IvoW8f6q6RZa9F0LfPyFdMEPiWHMIFKkal7+wsKCVlZXC7QAAAMBRFVKq/7Zz7lTv71u1v3wm\nAAAjJ4fxygAAAG2Sy5wxIT1jPiDpjHNuStKjku4dbJIAAAAAAABG12VVxnjvfyjptb2/vyfp5gbT\nBABAK9BjBQAAIK1cyl8hPWMAAACAI19gDpmHq/+crTkOq+bEsuYas/I2NnxnZycq3Jo3b319vXKb\nNffg+fN+Ncp4AAAgAElEQVTna++/n5W+zc3NynBrDjTr3rHm06uaQ7AsvD/PrHsvdo5La15Na/9N\nvhtC7v2j9K4KmTsVR1OSypjYL4vLFXLjhkw22cTDXJb2kHMKmVR2Y2OjdhzrS7LIU089VTs8ZDLe\nkDhVkwuXCck368vt+9//fuF260u7SMhEtCEToy4tmYuqHWLdp3/zN39zaFtIwaLb7daO0+l0ascJ\nmWDZOs4gv4xTfbGHvFPb7CgV3gAAAEZBLuWvPNaMAgAAAAAAaAmGKQEAUCKXlhkAAIC2yKX8Rc8Y\nAAAAAACAhOgZAwBAiVxaZgAAANoil/IXPWMAAAAAAAASomcMAAAlcmmZAQAAaItcyl9UxgAAAIyo\nomXt+7dZBd7x8bhO1EXHTxk/dv9F4bu7uy/8vbOzUxm/Kjzk2IPUfx5Ftre3K8M3NjYqw9fX1yu3\nraysVMY/d+5cZfhzzz1XGb66uloZvrm5WRlumZycrAxfXFyMCp+bmzu07cKFCy/8bT27ExPVP/M6\nnU5UuPVuaPLHdNM/1K39N/1sWscfxvnnUjmSGsOUAAAAAAAAEqJnDAAAJWgJAgAASCuX8hc9YwAA\nAAAAABKiZwwAACVyaZkBAABoi1zKX/SMAQAAAAAASChJz5jp6ena4SGzVFuzwg+KNXN+Eat2ryzt\nIcfa2tqqHceacb7I888/XzvOE088UTv8ySefHPhxijz99NO14ywvL9eOY92njzzySOH2mZmZ2sda\nWlqqHefkyZO145w4caJ2HOuZ+MlPfnJoW0gezM/P147T7XZrx5mamqodx1qpYJBStTCEHGfYqwa0\n9dgAAAA5yqX8Rc8YAAAAAACAhJgzBgCAErm0zABlmn4GrP2H9JSuEz92/1b6m+wBGXtuVu/rzc3N\nyvCNjY3K8IsXL1ZuO3/+fGV8q2e0Fb6yslIZbqXfurZWL1qrd/IVV1xRGV7UU7k/zyYmqn/GWT12\nQ3r09hsfr27Tb7LXa8i++7fFvteafm9Zveib7lHchrJPG9KQAj1jAAAAAAAAEqJnDAAAJayWP6Ap\nzrlJSZ+TdK2kaUkfk/RdSfdI2pP0iKTbvfdpJswDACCRXMpfeZwlAADA0fJ2Sc9572+U9MuSflfS\nXZLu7G0bk3TbENMHAAAi0DMGAIASuYxZRiv9qaR7e3+PSdqWdIOkh3rb7pd0i6T7qnYyNTV1qIUx\nZMU4/FTIKn1H0cLCwsD3ef311w98nzl561vfOuwkHFnW6r6olrqnSi7lLypjAAAAWsZ7f0GSnHML\n2q+UuVPSp7z3l2aGXJFUPUOoDk/C2u12tb6+PtjEVrAmsmz7BL4Hw+fn53XhwoUXPlsTbcak35r8\n1wqPncB3bW2tMnx1dbUy/Ny5cy/6fP311+s73/nOC5/Pnj1bGf8nP/lJZXhuE/i+9a1v1Z//+Z+/\n8PnEiRNRx7cqFWdnZyvDrQmAm5xk1qoYOBg+PT39ouvddMVC7HvPknoC3/Hx8Re963IZQpQCOQkA\nQImxsbFk/4CDnHMvlfSApM977/9EUv8v/wVJ5wojAgBwhKUsfw2zDEZlDAAAQMs4566W9GVJH/Te\nf663+dvOuVO9v2+V9PAw0gYAAOIxTAkAgBL0WMEQfVjScUmnnXOne9veJ+lu59yUpEf10zllShXd\nw/3bmh4mZGn6GYsdpmRpsru+lTfWEClrGNLBIWwHWcOU+odrFTl//nzltqeffroy/lNPPVUZ/sQT\nT1SGP//885Xh1nA9K/9nZmYqw61hXNb1K7q3+s/JOr41jMqKb6XPYuVfzLNjxbXee01LPYyorpC8\nT10eyqX8RWUMAABAy3jv36f9ypeDbk6dFgAAMHgMUwIAAAAAAEgoSc+YkNngm+4WG3OckG57oV1N\nU+WD1ZW1SMhqDFaX1qJwq5tpkWeffbZ2HKu7a5GDKwVcDqtb8OOPP164fXJysvaxQpalDLmuIfep\nNcv/8vLyoW3W/VPEWi2hyPb2du04Ic9QSPf50G7Dqbqkhhyn6VUHYuTSTRYAAKAtcil/0TMGAAAA\nAAAgIeaMAQCgRC4tMwAAAG3RpvKXc25S0h9JulbSjqR3SdqWdI+kPUmPSLrde1+7+zo9YwAAAAAA\nAA57q6QJ7/0vSfrnkj4u6S5Jd3rvb5Q0Jum2kB3TMwYAgBJtapkBAADIQcvKX9+TNOGcG5e0KGlL\n0mslPdQLv1/SLZLuq7tjKmMAAABGVFGBtk4hN7ZAHBs/1UIGZaz8s86vKtw6N2vCdmvi+K2trcpw\na8L+tbW1yvCiifb7Pf3005XbrIUTrPCnnnqq9vH7hSwK0M9aiMDKf2vC++np6UPb+he1sBZpmJ2d\nrQy3Fiuw7q+QRQjqqHp2hv1D3Tp+03kTu/CB9W4JWYxhxF3Q/hCl/yTpCkl/T9JN3vtLF2JFUvUL\noQSVMQAAlBh2gQ8AACA3LSt//a+SvuS9/5Bz7qWSviJpqi98QVL9JXbFnDEAAAAAAABFnpd0qSvg\nWUmTkr7tnDvV23arpIdDdkzPGAAASrSsZQYAAGDktaz89TuSPuece1j7PWI+LOmvJZ1xzk1JelTS\nvSE7pjIGAAAAAADgAO/9BUm/VhB0c+y+qYwBAKBEy1pmAAAARl4u5S/mjAEAAAAAAEiInjEAAJRg\neUcAAIC0cil/JamM6XQ6tcND1l8P6c60s7NTO07s2vB1WOvAFwk5p5DjbG5uDjxOUXjIcdbW1mrH\nWV1drR3n/PnzAz/O2bNnC7dPTNR/XLe3t2vHWVxcrB3niiuuqB1na2urdnjI+YQ8r6meces4ReGh\naUt1Trl0KwWOiqJnMuVzGvKeS8nKCyv/YvLSKntZ33nW96hVflpfX68Mt8orFy5cqAwvKs/0b3vu\nuecq4z/zzDOV4U888URUeEi5r591/taPyYWFhcrwixcvVm6z0r+xsVEZHlIO62fdn1a51cqfJt9T\nse+l2LQNO751fgffTZ1O50XbrN/2uHz0jAEAoASVSwAAAGnlUv7Ko/8PAAAAAABAS1AZAwAAAAAA\nkBDDlAAAKJFLN1kAAIC2yKX8Rc8YAAAAAACAhOgZAwBAiVxaZgAAANoil/IXPWMAAAAAAAASomcM\nAAAlcmmZQb729vYajd90uKWJ/deJU/V/d3Z2KuNubW1FhVv7393djdr/+vp6ZfjKykrltuXl5cr4\nZ8+erQx/8sknK8OfeeaZynAr/RYrf+bm5irDr7rqqsrwtbW1ym0bGxuV8a3rb4W3+dm24lrf3cN+\nL8WWLVKXTTqdzovul06n0/gxcyl/0TMGAAAAAAAgIXrGAABQIpeWGQAAgLbIpfxFzxgAAAAAAICE\n6BkDAECJXFpmAAAA2iKX8leSypimJ0GK2U/IhR4fr9+hyJokrUybz2lycnLgxykKb3pyvUu2t7eT\nHCf0eQi5hwb1bFlCJvKamKh+/RSFh9xzIdr8BZDqmqZk5XebrwcAAAAQgp4xAACUoCIIAAAgrVzK\nX8wZAwAAAAAAkBA9YwAAKJFLywxGV9EQ1/5t1j1uDY20htBa4db+Y4cxxp5fUXj/tpj82dnZqYxr\nhcfmbWy4Nby7KLx/28bGRmX8lZWVyvCLFy9Whq+urlaGr62tVYZbpqamKsOt8wudwuCSYX8/xQ6b\njnl2rGkPivbdv7+mh3xb6Wt6eHYTU4SkHiY/7Ps7lcuqjHHOvUbSJ7z3p5xzf1vSFyV9vxf8+977\nf9NUAgEAAAAAAEaJWRnjnLtD0jskXap+vkHSXd77324yYQAADFsuLTMAAABtkUv563LmjPmBpLf1\nfb5B0t91zn3VOfeHzrmFZpIGAAAAAAAwesyeMd77Lzjnru3b9A1Jn/Xef8s59xFJvynpN6r2MTU1\nVTl2rtvtXl5qR9ixY8eGevyf+7mfqx3nF37hF2rH+dVf/dXK8NOnT9fe56gZxaWLQ3z0ox8ddhKG\nbnFxcdhJaIXZ2dmhHdsa9w0AAIDByqX8FTKB733e+3OX/pb0aSvC5uZmaVi329X6+npAMg6zJjor\nEvLDN2TCrao4x44d07lz5wrDQvImZEKyZ599tnacxx9/vHac7373u6Vhp0+fLvwB/thjj9U+Tkic\nH/3oR7XjPP/887XjVE0ot7e3V9otr9Pp1D5WSCXftddeWzvOK17xitpxrrvuutKwj370o4UVc698\n5StrH+flL3957ThXX3117Tjz8/O140xPT5eGLS4u6vz584e2h3bbDLl/Qo4VEqfq/Tg7O1v4zAyz\nggYAAACIFVLl9CXn3C/2/v47kr41wPQAAAAAAACMtJCeMf+zpE8757YkPSnpnww2SQAAtEMuE8gB\nAAC0RS7lr8uqjPHe/1DSa3t//wdJr2swTQAAABiAomGA/dusAq81nNsaum3Ft8JjC+Sx6bfixJxf\nyPD6ftacCk3nrZV329vbldus87fOzwq3huda4Vb+TE5OVoZPTFT/zAoJ7982NTVVGT/2/GPzN3bO\nj5hnp+ja1XnWY9Nu3TtNv9diwwcVB7aQnjEAAGQhl5YZAACAtsil/JXHNMUAAAAAAAAtkaRnjNUt\nrCg8VW1YquOE5IFU3MXTsrW1VTtOquOEaHO3uJBujFa30rJwK16Rubm52nFCVqkJWZ7e6t5bFB6S\n3yHdwEOehyZWcysKD+06m2plpJA4sV3Rm5RLywwAAEBb5FL+omcMAAAAAABAQswZAwBAiWH2ygEA\nAMhRLuWvPM4SAAAAAACgJegZAwBAiVzGLAMAALRFLuUvKmMAAABGVNFE4HUmxrcKxJ1OpzLcWsDA\nSkvTXdWt84tdRMA6/yqxeW+x8tY6vhV/YuLwz4z+bdPT05Xx5+fnK8OPHz9eGW5NrG8tPmBd+2PH\njlWGnzx5sjLcOr+iBRX6t1kLLlgLPxRdn37W/dX05PtV+R/yXPbfD7E/9Jt+bwy7IqLo+MNO06ii\nMgYAgBIUPgAAANLKpfzFnDEAAAAAAAAJ0TMGAIASubTMAAAAtEUu5S96xgAAAAAAACREzxgAAErk\n0jIDAADQFrmUv+gZAwAAAAAAkBCVMQAAAAAAAAklGaZkrbVeFB7SNSkkTsg68LFr29cJ397ern2s\nzc3N2nG2trZqxwkxPl5d/1cUPjFR/zadnp6uHafb7daOE8K6f6644orC7SHndOLEidpxjh8/XjvO\n3Nxc7TgzMzO1wycnJ2sfJ8Tu7m6SOCGsZ6hMLt09By00v4FcWO8W6xmyvhOt/YeUyQYp5vytvLG+\nV6xzt+Jbx7fKX1a5pKhs0L9taWmpMv7GxkbU8Y8dO1YZvrq6WhluscpYV111VVT8xcXFym1W2cu6\nfp1Op9H4TZY7Yp/72DJb7LnFli1iz99Kf1F46nJkLuWvPM4SAAAAAACgJZjAFwCAEvQoAgAASCuX\n8hc9YwAAAAAAABKiZwwAACVyaZkBAABoi1zKX/SMAQAAAAAASIieMQAAlGhLy4xzblLS5yRdK2la\n0sckfVfSPZL2JD0i6Xbv/a5z7l2S3i1pW9LHvPdfHEaaAQAAQrSl/NU0esYAANB+b5f0nPf+Rkm/\nLOl3Jd0l6c7etjFJtznnrpH0Xkmvk/QWSb/lnKte/xUAAADJ0TMGAIAS4+OtabP4U0n39v4e036v\nlxskPdTbdr+kWyTtSPq6935D0oZz7jFJr5L0zbTJRVsU3cN17murddLa197eXqvDrfMrCu8/Z2v/\nnU4n+NiWnZ2dyvCtra3KcOvaTU5OVoYvLCxUhp88ebJy2/b2dmX8qampyvDFxcXK8GPHjlWGb2xs\nVIZb13Zubq4y/IorrqgMv/rqqyvDr7zyyspt1vFnZmYqw63rW3XvplD1fFhpK7p2dZ5by7DfW5aQ\n95oVnrqnSovKX42iMgYAgJbz3l+QJOfcgvYrZe6U9Cnv/aUS24qkJUmLkpb7ol7aDgAAgBbJo8oJ\nAIAAY2Njyf5ZnHMvlfSApM977/9E0m5f8IKkc5LO9/4+uB0AAOBISFn+Gub8NFTGAADQcs65qyV9\nWdIHvfef623+tnPuVO/vWyU9LOkbkm50znWdc0uSXqn9yX0BAADQIgxTAgCgRItm8/+wpOOSTjvn\nTve2vU/S3c65KUmPSrrXe7/jnLtb+xUz45I+4r1fH0qKAQAAArSo/NWoJJUx1iRLReEhFyBksiNr\n8rNhCzkna0K0Iru7u/Z/OiDkGlmTsRWFz87O1j7O8ePHa8dZW1urHWd9vf5vHOt5eMlLXlK4fXq6\n/oIoIflgTXhXZH5+vnYca2K5onDr/iliTVBXJGTSsFy+NCwh7xJL7ER2o8B7/z7tV74cdHPB/z0j\n6UzjiQIAAEAwesYAAFCCSjYAAIC0cil/MWcMAAAAAABAQvSMAQAAyJTV+hjbOmkN+2xiaGO/2GGO\nReffvy0mf6y8sdJu5Z21fyvt1tBgaxj51VdfXbltYqL6Z4g1BHplZaUy/OLFi5XhW1tbleFW/s/N\nzVWGnzhxojL8mmuuqQwvGmp+8uTJF/628qfb7VaGW8O4rfsj9v6MefeExO1/Hqy0N/HeGGT8pnuN\nFJ0/Q8abQWUMAAAlcukmi/ZxznW0P/ePk7Qn6T2S1iXd0/v8iKTbvffN1mYAAJBYLuUvhikBAAC0\nz69Ikvf+dZLulPRxSXdJutN7f6OkMUm3DS95AAAgBj1jAAAoEbKyFjAI3vt/65z7Yu/jyySdk/Qm\nSQ/1tt0v6RZJ9w0heQAANKZt5S/n3Ick/deSpiT9nva/i+9RZE9VKmMAAABayHu/7Zz7I0l/X9J/\nK+nN3vtLA/dXJC1Z+5iamjpUqLXmkkA1a66QtlhcXBx2Eg554xvfOOwkHGmvf/3rh52EI2tmZmbY\nSTjSrDmkRplz7pSkX5L0Okmzkn5DP+2p+qBz7jPa76lau3GEyhgAAErkMmYZ7eW9/0fOuQ9K+veS\n+n9NLGi/t0ylzc3NF33udrtaX19/4XPsRJGx4dYkn9akkTs7O1H7txxM/9zcnDkx7OWKncB3e3u7\nMrz/OhdZW1uLCrcm0D0Y/sY3vlFf+cpXXvj83HPPVcZfXl6OOv6oTeD7+te/Xl/72tde+Hzs2LHK\n+Fblg5X+6enpynBrAuYm3x11v5tnZmZedD83PYFv7OTZw57A96CpqakXfZekqJhpWfnrLZL+o/Yr\nWxYl/W+S3qUB9FRtV/8fAAAAyDn3jl63aElalbQr6a97LXSSdKukh4eRNgAAMnKFpFdL+u+0P5n+\nv5I0XrenahF6xgAAUKJlLTPIy59J+pfOua9KmpT0fkmPSjrjnJvq/X3vENMHAEAjWlb+ek7Sf/Le\nb0ryzrl1SS/tC7+snqpFqIwBAABoGe/9RUm/VhB0c8p0tK27fF2x6SuK35ZzttJhDSOZnJysDI8d\nQlaUvv55bDqdTmV8axjN6upqZbg1zOrgEL6DrPydnZ2tDD958mRluDWnz/z8/KFt/UOXrPyxhpJY\n4bETqMYO9Rnmc9aWZzxUbN4XxY8d8nnEfU3S+5xzd0l6iaQ5Sf/OOXfKe/+g9nuqPhCyYypjAAAo\ncdQLZAAAAEdNm8pf3vsvOudukvQN7U/zcruk/6wB9FRtbWVMyAWIrQW8XCFps1ofysKteEWsCbcG\ndZyFhYXacazJzIrCNzY2ah8nRMjqEiFps1qqrrvuusLtIZNlhaz4EHJdr7jiitpxQu6FpaX6wzFD\n8iDkXgh57qxWwaLw0C+nkPdjqneqdZzMW2MAAAAwRN77Owo2R/dUbW1lDAAAwxbbTRsAAAD15FL+\nyuMsAQAAAAAAWoKeMQAAlGjTmGUAAIAc5FL+omcMAAAAAABAQvSMAQCgRC4tMwAAAG2RS/mLyhgA\nAAA0wlotbdjhliZ/EMSuFGedm5V2a0W/kFUc+xVNwNm/YqO1uqS1SqG1+uPOzk5luJV/Vv5YK5HO\nzs5GhRftv3+FSCt/rPRZ52eJvX+tCVpjnt2iuP3ptdIe+9zH5q11/Nj0pVotEzaGKQEAAAAAACRE\nzxgAAErk0k0WAACgLXIpf9EzBgAAAAAAICF6xgAAUCKXlhkAAIC2yKX8Rc8YAAAAAACAhOgZAwBA\nCWu1BwAAAAxWLuWvJJUxIcsOhiyXlmqZrpBuU6FLlIUsK9jtdmvHsZYHLBKybJsV52d+5mcObbOW\n5isSEqd/ucDLtbm5WTuOlbaf/dmfLdw+MzNT+1gh5xQS5/jx47XjvOQlL6kdvrS0VPs4i4uLteOE\n5EHIPWc9D0Xhod02Q+LFLls5KLl0VQUAAEA+6BkDAEAJKoJw1BU1VNVpvAppUGtTeBOVyv37HOY7\nwmo5ttJmxd/e3q4MtxoUJiYO/8yYnZ2tDO83PT1dGR57bWPzxzp/q3HUOv+i8P7GopAGnX7W+cU+\ne5bQhurLObb13hv2ucXmfawmOhYMWi7lrzz6/wAAAAAAALQEPWMAACiRS8sMAABAW+RS/qqsjHHO\nTUr6nKRrJU1L+pik70q6R9KepEck3e69b8fEAgAAAAAAAC1nDVN6u6TnvPc3SvplSb8r6S5Jd/a2\njUm6rdkkAgAwHGNjY8n+AQAAIG35a5hlMGuY0p9Kurf395ikbUk3SHqot+1+SbdIuq9qJ9PT05UT\nFYWsEjNqQlaJGUVvectbhp2EofvkJz857CS0whve8IZhJ2HoeDfu65/wEQAAABgFlZUx3vsLkuSc\nW9B+pcydkj7lvb80xfOKJLMWYWNjozRsZmZGa2trh7aH1FCFzDwdMhP7oOMsLS1peXm5MKwq78pc\nvHixdpyVlZXacVZXV2vHKTtPab8i5ktf+tKh7U899VTt4zzxxBO145w9e7Z2nEEvbf3JT35Sd9xx\nR2FYTktbv+ENb9ADDzxwaPuJEydqHyckTqqlratWUhjkuzE0Xqqlrave3bOzs4XvmlQVNNaKBwAA\nABisXMpf5lk6514q6QFJn/fe/4mk/tL5gqRzDaUNAAAAAABg5FgT+F4t6cuS/qn3/t/1Nn/bOXfK\ne/+gpFu1X1EDAMDIYS4XHHVFPc/q9CRuuoeclZaQXs+DZOVfzDsi9v1ixbeundXybIVbPUKL0jc1\nNfXC351OpzJ+Ve/Rsv33iz2/2P1b+RNy/P78s9I37O+v2ONXPftNvxeafi/Fxo+9ty1F1y51T5Vh\n37+pWHPGfFjScUmnnXOne9veJ+lu59yUpEf10zllAAAAAAAAYLDmjHmf9itfDrq5meQAAAAAAACM\ntjxmxgEAAAAAAGgJa5jSQOzs7NQOTzUuLWQ8Wkic0LF9ISu0hKy6E3JO3W63dhwrbddcc82hbSGr\npoTECVl1x7q3i1j59vM///OF20Oua8g1CokTkt/Wcu4nT548tG1hYaH2cfrHV18uaxx7kZB3Vsh7\nIeVqSiGaGMedy7hhAAAA5CNJZQwAAEcRFUEAAABp5VL+YpgSAAAAAABAQvSMAQCgRC4tMwAAAG2R\nS/mLyhgAAIARVTS3Wf82q8Db9Bx+sQVuK74V3sQ8V4MSe22avnZW3hUdv38et93d3cr4ExPVP1Ni\n86fpez82flH6BnlNretnhTed/1XhTczf2c86d+vetcJj87bp91ouFSFtQGUMAAAlKJAAAACklUv5\nizljAAAAAAAAEqJnDAAAJXJpmQEAAGiLXMpf9IwBAAAAAABIiJ4xAACUyKVlBgAAoC1yKX/RMwYA\nAAAAACAhesYAAFAil5YZAACAtsil/EVlDAAAwIja29u7rG1ldnd3K8NjC8xW/Nj9j48PvhN4nX3W\nyeu6rH3Hhlusa9PpdCq3WflYFL+O2HvLCrfyz4pvPVshxxzWvkJY+RPz7BadW//xrHsr5Nr029nZ\niYrftJB3Q2yeoFiSypiQCx7yggj5wg6JE/JyCH0oQ9IX8uXV7XZrx5mamqodZ3p6ujL82LFjteMU\nmZycrB3n5MmTteOE3KcTE9WP3ctf/vKgeEVC7oWQ/A65f+bm5irDjx8/fmhbyD0XEidVbXzTBeXY\nfaUqqDX9YyxGLi0zAAAAbZFL+Ys5YwAAAAAAABKiMgYAAAAAACAhKmMAAAAAAAASojIGAAAAAAAg\nIVZTAgCgRC4TyAEAALRFLuUvesYAAAAAAAAkRM8YAABK5NIyg9E1Pn643a1oW5mmn4G9vb2o+MNO\nnxW+u7sbfOw616nIsPO2KH7/OcWmz2Lln3V86/x3dnYqw61rb4UXHb8/jpW+YX9/NXl9Q86tThzr\n3mni2Rjk/i2x934KbUhDCvSMAQAAAAAASIieMQAAlMilZQYAAKAtcil/0TMGAAAAAAAgIXrGAABQ\nIpeWGQAAgLbIpfxFzxgAAAAAAICEkvSMsWakLgpPVRvWluMMMh0hs++nygcrbdPT07XjFOl0OrXj\nbG9v144Tkm/W+VxzzTUDO1ZInKmpqdpxJicna8cputb95ubmDm1rIr+HKeS90PTqE/1CVgFp4l2S\n8pwPyqVlBgAAoC1yKX8xTAkAAGBExS5tPepiGzpilohtetns2B8zTSwNPcilrWMbOmPz39J0Q0Ib\nly5PJeTa9997sfdG089e7L1thVuN1kXnn0vlSGpUxgAAUILCBwAAQFq5lL9oGgEAAAAAAEiInjEA\nAJTIpWUGAACgLXIpf9EzBgAAAAAAICEqYwAAAAAAABKiMgYAAAAAACAhKmMAAAAAAAASYgJfAABK\n5DKBHEbX+PjhdreibaH29vYGtq8Q1jNqhcemv8l3ROy5DfI6D0Ps+cdem93d3cpw696x4ocY5PPW\ndP4N8/4rSnv/NittsflsxbfCh33vp9pnm443LEf7LQ0AAAAAAHDE0DMGAIASubTMAAAAtEUby1/O\nuaskfUvSmyVtS7pH0p6kRyTd7r2v3R2OnjEAAAAAAAAFnHOTkv5A0lpv012S7vTe3yhpTNJtIftN\n0jNmcnKydnjIWL1U45aHPT66CSG1jyFjQa04U1NTh7Z1Op3ax5mYqH9rNzG2t4iV1/Pz8wM7Vsg1\nKoUVZl0AABixSURBVLoGlpBrZKXNem80KSTfQu456zhF4aHvn5D7O+RYbWzJiDFq5wMAANB2LSx/\nfUrSZyR9qPf5BkkP9f6+X9Itku6ru1N6xgAAAAAAABzgnPt1Sc9477/Ut3nMe3+p1XJF0lLIvpkz\nBgCAEi1smQEAABhpLSt/vVPSnnPuTZKul/THkq7qC1+QdC5kx/SMAQAAAAAAOMB7f5P3/mbv/SlJ\n35H0DyXd75w71fsvt0p6OGTf9IwBAKBEy1pmgNqK7uFB3tex+2p6Hj4rfVZ4Ufr65/OKSb8VNzbt\nTV8bK7xorrI685c1/f610hJ7/iHzz/UrOv9B3Xuhx68Tbp1/k+m33nux95Z1bjs7O5Xhw362j4Ij\ncI4fkHTGOTcl6VFJ94bshMoYAAAAAACACr3eMZfcHLs/KmMAAChxBFpmMOKcc1dJ+pakN0valnSP\npD1Jj0i63XufZilAAAASyaX8xZwxAAAALeScm5T0B5LWepvuknSn9/5GSWOSbhtW2gAAQBx6xgAA\nUCKXlhm01qckfUbSh3qfb5D0UO/v+yXdIum+qh1MTk4emt9genp6sKnMTLfbHXYSjqz5+flhJ+FI\nW1oKWj0XkiYm0v3sTXmsVCYnJ5MeL5fy1+jdKQAAAEecc+7XJT3jvf+Sc+5SZcyY9/7SrJcrksxf\nZltbWy/6PD09rY2NjUEmNcqwJyG1HExft9vV+vp6aXjMvg8a9iSfIRP0VoXPz8/rwoULl318a5LU\n2PO3Jlm1zs+KH+tg+peWlrS8vPzC59hnx8rf2PyPncC3Kn7dtE9MTGh7e3tgabPCY++Ntj37k5OT\nL/ouSV0xM8qojAEAAGifd0rac869SdL1kv5Y0lV94QuSzg0jYQAAIB6VMQAAlMilmyzax3t/06W/\nnXMPSnqPpH/hnDvlvX9Q0q2SHhhO6gAAaE4u5S8qYwAAAI6GD0g645ybkvSopHutCEUF2pSF3KaH\nIcUOpbAUDVUZVP7F7qfpoQp1hyEdVHTt+7dZ6YsZxnI58WPzL/b4FuvZHfb9M8zwkLh18q7p91as\nYQxBtJ53hBmpypiQByfVwxb6wrFe9IMS8oCFvAis8+l0Ope1zTI1NVU7TqqXjJVvs7OzAztWyP0T\nEifkXggpJIU8r03cp6niFIWH3qep8qGJH3m5tI4AZbz3p/o+3jysdAAAgMFhaWsAAAAAAICERqpn\nDAAAg0SvHAAAgLRyKX9RGQMAwBHhnHuNpE947085566TdI+kPUmPSLrde7/rnHuXpHdL2pb0Me/9\nF4eWYAAAABSqrIxxzk1K+pykayVNS/qYpB9L+qKk7/f+2+977/9Ng2kEAGAo2tQy45y7Q9I7JF3s\nbbpL0p3e+wedc5+RdJtz7i8lvVfSqyV1JX3NOfcX3vuNoSQaAACgpjaVv5pk9Yx5u6TnvPfvcM6d\nkPQdSf9c0l3e+99uPHUAAOCSH0h6m6TP9z7fIOmh3t/3S7pF0o6kr/cqXzacc49JepWkbyZOKwAA\nACpYlTF/qp8umzim/S7PN0hyzrnbtN875v3e+5XmkggAwHC0qWXGe/8F59y1fZvGvPeXlhhbkbQk\naVHSct//ubQdAADgSGhT+atJlZUx3vsLkuScW9B+pcyd2h+u9Fnv/beccx+R9JuSfqPyIBMTlUuk\nhixDPGrm5uaGnYRWWFriN8OxY8eGnYRWmJ+fH3YShm5QS70fdd1ud9hJaKv+dc4XJJ2TdL7398Ht\nyFRRgbZNhdy9vT37P1XY3d2tDLfemVZeFJVf+7dZx2+SlXYrPDbvY4Xk/SD3b7GOH7t/694p2n9s\nnlj7rxNupWWQaU1t2M+Otf+Qe6fO/pGOOYGvc+6lku6T9Hve+z9xzh3z3l8q2N0n6dPWPra3t0vD\npqamtLm5eWh7yE0S8oWY6masOs7c3JwuXrxYGLazs1P7WCH5EBIn5CVb9XJYWlrS8vJyaXjTUhWo\nqvLg2LFjOnducL+bQq7RoK9rmapnYn5+XhcuXKgVp0xI2kIqPCYm6s+HXpXXnU6n8PkPvU9TvRcG\n/SOv2+1qfX29cHsKbfrRWuDbzrlT3vsHJd0q6QFJ35D0cedcV/uNJ6/U/uS+AAAAR0LLy18DU/mr\nyzl3taQvS/qg9/5zvc1fcs79Yu/vvyPpWw2mDwAAFPuApP+jN2nvlKR7vfdPSrpb0sOSviLpI977\nw7VZAAAAGCqrKffDko5LOu2cO93b9s8k/Y5zbkvSk5L+SYPpAwBgaNrWMuO9/6Gk1/b+/p6kmwv+\nzxlJZ9KmDAAAYDDaVv5qijVnzPskva8g6HXNJAcAAAAAAGC01Z/kAACATPz/7d1fiCVZfQfwb3dP\nO+Ps7C6+KJGIPiSep7CRlWTFfwNxjRpEE0KeTIyLBBMflIREI7vkRQkBNaBiDKvrGqMQNK4EQd0H\n/ysS/BNQlLPZxIdAEogL7u7sZpedmc5D98TrzL1d3VV9T1V3fT4w0LdqqurcU6dunfur3zl3Lk9m\nAACmYi79rybBmK5JII9q4tRWk/72aRx9ZyzvM5Foq8bbZ6LXrvezbBLUVhPKtposucuqOmo1EW0f\nfa69VpPxttqmj646OMoJxlu1nz7HGfqrAAAAcNwc398cAwAAADiGDFMCADihlmWXHibjtCtzrSuD\nb2gGadf+uzLnjjLDsLV1l33d2ZnLzs3isr5Z41f0yZw+jKH1s47z1zJTdOj768o4X+f5W1a2ozwf\nQ8s+9HO1y3H+3JsbmTEAAAAADcmMAYAVzFcDANDWXPpfMmMAAAAAGpIZAwArzOXJDADAVMyl/yUz\nBgAAAKAhwRgAAACAhgRjAAAAABoyZwwArDCXMcucXMva8FG2652dnbVuP3T/ly9f3nd9n7ro2udB\n9z923XWtv3Tp0qDtNzf3f+Y7tB12HX/o/ru2H3r+xrbua6tr/fb2du9j9zk3i8uGvvd1t93DfMYs\n03XtHQdz6X8d/zMFAAAAcIzIjAGAFebyZAYAYCrm0v+SGQMAAADQUJPMmK4xp8vW9xnL12qbPpG6\nrm1WjQ1cx7GW6TO2cB1lW7a+T9lavZ91tLlTp5Zflq3aQqtrostxG+86lXpbpVV9rqPNjfl0ZC5P\nZgAApmIu/a/j9W0HAAAA4JgTjAEAAABoSDAGAAAAoCG/pgQAMFNdczYNXb9qTryDGjqn1NDyL9v/\n4jZD5jUYWvYuQ8/N0ON3mfr+T/qcFV3nf+i126Wrfre2tnpvu2z94rKhbb/vXKAH3X+XdXyuMQ7B\nGABYQYcFAKCtufS/DFMCAAAAaEhmDACsMJcnMwAAUzGX/pfMGAAAAICGBGMAAAAAGhKMAQAAAGjI\nnDEAsMJcxiwDAEzFXPpfgjEAACw1tEPctf3Ozs5atx+6ftnxL1269P9/b21tHXr7oyrb5cuXB20/\ndP99zs3isqHlH6qr/Ov+Mtjn/S1uM7T8m5vrHSAx9PhD6r+r7a372hpq3cfvc+7X/Z7navbBmFZR\nt743xD4Nv9V76nOcPjeOdRxnma4O1TJDb6SHKceUI8R9ytZ1ozl16tqPpz713XWcMfW5HvreDKd8\nEx36hWGdpnzdAQCcRHPpf5kzBgAAAKCh2WfGAMAqc3kyAwAwFXPpf8mMAQAAAGhIMAYAAACgIcEY\nAAAAgIbMGQMAK8xlzDIAwFTMpf8lGAMAcEIt+2n4xWVT/mn5ozB2h36/+uuq28uXL/fe90HWd+mq\nu3XX7djvb3Nz/wEEY7etobrKv7W1tdb9D6m/rmtj2b4X20vX9kOvvaHvvavtDT1+l2XbH/f2PlWC\nMQCwgs4HAEBbc+l/mTMGAAAAoCHBGAAAAICGDFMCgBXmkiYLADAVc+l/yYwBAAAAaKhJZkzXbNzL\n1veZIf2kbdP6WC30mb28T2S0Tx30Oc46tlm1fsoR4j5l6/O50NV+lmm1TR/71dvW1tbScvQtW6vP\nknVcr63OBwAAtCIzBgAAAKAhc8YAwApTzkiDo9CVmTY027ZvNuhR6Sp/n+Nvbv70WWbX9kPqb92Z\nzkPPTZ/1h6nvoXU79Nyvo+0M3X5K96TF62CZoWXdr/77nJvFLNuujNuxr72pb9/ClMpYStlOcleS\nZyU5neTtSX6Q5O4kO0m+n+SNtdZDp3LLjAEAAAC41muSPFBrfWGSlyV5X5J3J7l9b9lGklf12bHM\nGABYYUpPZgAA5mBi/a9PJPnk3t8bSS4muTnJl/eWfTbJS5Pcc9gdC8YAAAAAXKXWeiFJSinXZzco\nc3uSd9Zar4xnezjJjX32bZgSAAAAwBKllGck+WKSj9ZaP55kcX6Y65P8pM9+BWMAAAAArlJKeVqS\ne5O8pdZ6197i75ZSzu/9/fIkX+2zb8OUAGCFiY1ZBgA48SbW/3pbkqckuaOUcsfesjcleU8p5UlJ\nfpifzilzKIIxAAAAAFeptb4pu8GXq7146L4FYwBghYk9mYEjt7Oz0/2fBhh6DXVtv+7yL7NYpiHv\nr6vsY3/+bG4Om81gWfmP8j0NLV/X9mPX/7LjL5Z5aNsfu/2NWb9Tf29D2/Y6tD5fY19/rUzvTAMA\nAACcYDJjAGCFuTyZAQCYirn0v5oEY7a2tg69vk/q3eXLl7v/00jH6WpQq9LRLl26dOhj9XlPfRr8\nOrZZtn6MFOR16vvhMpW2elT2O87W1tbSsve5HvrUQR9HfZzt7e2l77fvcVp9Pq6j/fQ57wAAMGWG\nKQEAAAA0JBgDAAAA0JBgDAAAAEBDnXPGlFK2ktyZpCTZSfKGJI8luXvv9feTvLHW2mZiBgBoZC4T\nyAEATMVc+l8HmcD3lUlSa31+KeV8knck2Uhye631S6WUDyR5VZJ71lZKAAAObdlE3IvLuibqHjqJ\n/qofKLiiq8O97vKt237vb2jZh9ZNnx9VOEpjt62pv79l5Vtc1lX+oT8sMHT7rh9wGXJ+xv5cWMe5\nPcz+h7bdPuWf+mftcdUZjKm1frqU8pm9l89M8pMkL0ny5b1ln03y0uwTjDl16tS+jWZ7e/ug5T2x\nrrvuurGLMAlnz54duwijcz3sWlYPc6ubM2fOjF2ESTh37txox57LkxmmqZTynSQP7b38UXYfiN0d\nmckAnGBz6X8d6Keta60XSykfSfKbSX47ya211ivhsYeT3Ljf9hcvXly5bnt7O0888cQ1y0/aT1vv\n57rrrssjjzyydN2Uf9q6K+J92OOcPXs2jz766KG26XOcKW+z6nro67j+tPWqeuhzPbT6WeQ+9bbf\nU70zZ87kscceu2b53H7a+ty5c7lw4cLS5XCSlVLOJNmotZ5fWPZPkZkMACfCgSfwrbW+Nsmzszt/\nzJMXVl2f3WwZAACOxk1JzpZS7i2lfKGUckuSm/OzmckvGa10AMAgB5nA93eT/Hyt9S+TPJrkcpJv\nlVLO11q/lOTlSb641lICAMzLo0nemeSDSX4xu8GXjcNkJie7WaxXZ6ydOnVq6d8cjCGk/RmWP4y2\n15+6G+b06dNjF+FEOsgd+FNJPlxK+UqS7SRvTvLDJHeWUp609/cn11dEABjHXMYsM0n3Jbl/L/hy\nXynlgexmxlxxoMzkq4dqnjp16meGj3cN5Rw6aePQSTy7jt81BPOoJ528eghp1ySy++kq+9BzM7UJ\nfK8elr/uSVCP+wS+V7etVcOXV+lqX0Ovna6233XtD7l2Dtv2r667odfO0Gtr7MmnD9s2T58+nccf\nf/xnXq/bXPpfB5nA95Ekv7Nk1YuPvjgAACS5LckvJfmjUsrTk9yQ5F6ZyQBwMshNBYAV5vJkhkn6\nUJK7Sylfy+6vJ92W5MeRmQzACTeX/teG3wwHgOUefPDBZjfJG2+8cR49DwCAfbTsfyXj9cH6D9YD\nAAAA4NAEYwAAAAAaMmcMAKwwlzHLAABTMZf+l8wYAAAAgIYEYwAAAAAaEowBAAAAaEgwBgAAAKCh\n0SbwLaVsJnl/kpuSPJ7k9bXW+8cqz1hKKd9J8tDeyx/VWl83ZnlaKqX8apK/qrWeL6X8QpK7k+wk\n+X6SN9ZaL49ZvlauqofnJPlMkn/dW/03tdZ/GK9061VK2U5yV5JnJTmd5O1JfpCZtYUV9fAfmVFb\nSJJSylaSO5OU7J7/NyR5LCO2h7lMIMfJp9/Vj75KP+7vw0zxfnjclFKemuTbSW5NcjHq7sCu/n6a\n5B1pXH9z6X+NmRnz6iRnaq3PS/LWJO8asSyjKKWcSbJRaz2/929OgZg/S/LBJGf2Fr07ye211hcm\n2UjyqrHK1tKSerg5ybsX2sSJ/vKd5DVJHtg77y9L8r7Msy0sq4e5tYUkeWWS1Fqfn+T27N7859ge\nYB1m3+86LH2VQdzfh3E/HGAvGPi3Sf53b5G6O6AV30/V35qMGYx5QZLPJUmt9ZtJnjtiWcZyU5Kz\npZR7SylfKKXcMnaBGvq3JL+18PrmJF/e+/uzSV7SvETjWFYPv1FK+Uop5UOllOtHKlcrn0hyx97f\nG9l9cjHHtrCqHubUFlJr/XSSP9h7+cwkP8k82wOsg37X4emr9Of+PoD74WDvTPKBJP+591rdHdyy\n76fqb03GDMbckOTBhdeXSimjDZsayaPZ/bD49eymH35sLnVQa/3HJE8sLNqote7s/f1wkhvbl6q9\nJfXwz0n+tNb6oiT/nuQvRilYI7XWC7XWh/cCDZ/M7tOf2bWFFfUwq7ZwRa31YinlI0nem+RjmWF7\ngDXR7zokfZX+3N+Hcz/sp5Ty+0n+p9b6+YXF6u7grvl+GvW3NmMGYx5Ksvikd7PWenGswozkviR/\nX2vdqbXel+SBJD83cpnGsjju8PrsPgGYo3tqrd++8neS54xZmBZKKc9I8sUkH621fjwzbQtL6mF2\nbeGKWutrkzw7u+Pln7ywqnl72NjYaPYP1ky/a7hZ3p/6cn8fbkr3w2PktiS3llK+lOSXk/xdkqcu\nrFd3+1v2/fRpC+ub1F/L/teYfbAxgzFfT/KKJNlLf/reiGUZy23ZG7NdSnl6dp9a/deoJRrPd0sp\n5/f+fnmSr45YljF9vpTyK3t//1p2Jx47sUopT0tyb5K31Frv2ls8u7awoh5m1RaSpJTyu6WUP997\n+Wh2O+7fmlt7gDXR7xpudvenvtzfh3E/7K/W+qJa64trreeT/EuS30vyWXV3YMu+n96r/tZjzPTU\ne7IbtfxGdseSzmby2gUfSnJ3KeVr2Z2d+rYZP6X6kyR3llKelOSH2U1pnaM/TPLeUsoTSf47Px0v\nfFK9LclTktxRSrkytvxNSd4zs7awrB7+OMlfz6gtJMmnkny4lPKVJNtJ3pzdNuCzAYbT7xpOX+Xg\n3N+HcT88Wq7dg7vm+2mSH0f9rcXGzs5O9/8CgBm6cOFCs5vkuXPnjFUCAGavZf8rGa8PZuI2AFjB\nXC4AAG3Npf815pwxAAAAALMjGAMAAADQkGAMAAAAQEPmjAGAFeYyZhkAYCrm0v+SGQMAAADQkGAM\nAAAAQEOCMQAAAAANCcYAAAAANGQCXwBYYSoTyJVSNpO8P8lNSR5P8vpa6/3jlgoA4OhNpf+1bjJj\nAGD6Xp3kTK31eUnemuRdI5cHAIABBGMAYPpekORzSVJr/WaS545bHAAAhhCMAYDpuyHJgwuvL5VS\nDDUGADimdOQAYIUJjVl+KMn1C683a60XxyoMAMC6TKj/tVYyYwBg+r6e5BVJUkq5Jcn3xi0OAABD\nyIwBgOm7J8mtpZRvJNlI8rqRywMAwAAbOzs7Y5cBACbpsccea3aTPHPmzDxycgEA9tGy/5WM1wcz\nTAkAAACgIcEYAAAAgIbMGQMAK8xlNn8AgKmYS/9LZgwAAABAQ4IxAAAAAA0JxgAAAAA0JBgDAAAA\n0JAJfAFghblMIAcAMBVz6X/JjAEAAABoSDAGAAAAoCHBGAAAAICGzBkDACvMZcwyAMBUzKX/JTMG\nAAAAoCGZMQAAAABXKaVsJnl/kpuSPJ7k9bXW+49i3zJjAAAAAK716iRnaq3PS/LWJO86qh3LjAGA\nFeYyZhkAYCom1v96QZLPJUmt9ZullOce1Y5lxgAAAABc64YkDy68vlRKOZKkFsEYAAAAgGs9lOT6\nhdebtdaLR7FjwRgAAACAa309ySuSpJRyS5LvHdWOzRkDACtMbMwyAMCJN7H+1z1Jbi2lfCPJRpLX\nHdWON3Z2do5qXwBwoly6dKnZTXJra2tSPQ8AgDG07H8l4/XBDFMCAAAAaEgwBgAAAKAhwRgAAACA\nhkzgCwArTGwCOQCAE28u/S+ZMQAAAAANCcYAAAAANCQYAwAAANCQOWMAYIW5jFkGAJiKufS/ZMYA\nAAAANCQYAwAAANCQYAwAAABAQ+aMAYAV5jJmGQBgKubS/5IZAwAAANCQYAwAAABAQ4IxAAAAAA2Z\nMwYAVpvHoGUAgOmYRf9LZgwAAABAQ4IxAAAAAA0JxgAAAAA0JBgDAAAA0JBgDAAAAEBDgjEAAAAA\nDf0fZGKNM7PVbGYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x122ff1a20>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import pylab as plt\n",
    "plt.figure(figsize=(20,10))\n",
    "plt.subplot(1,2,1)\n",
    "plt.imshow(pacs100_psf[1].data[centre100-radius100:centre100+radius100+1,centre100-radius100:centre100+radius100+1])\n",
    "plt.colorbar()\n",
    "plt.subplot(1,2,2)\n",
    "plt.imshow(pacs160_psf[1].data[centre160-radius160:centre160+radius160+1,centre160-radius160:centre160+radius160+1])\n",
    "plt.colorbar()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Set XID+ prior class"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#---prior100--------\n",
    "prior100=xidplus.prior(im100,nim100,im100phdu,im100hdu, moc=Sel_func)#Initialise with map, uncertianty map, wcs info and primary header\n",
    "prior100.prior_cat(prior_cat(prior['RA'] ,prior['Dec'] ,'GAMA09_Ldust_prediction_results.fits',ID=prior['help_id'])#Set input catalogue\n",
    "prior100.prior_bkg(0.0,5)#Set prior on background (assumes Gaussian pdf with mu and sigma)\n",
    "\n",
    "#---prior160--------\n",
    "prior160=xidplus.prior(im160,nim160,im160phdu,im160hdu, moc=Sel_func)\n",
    "prior160.prior_cat(prior_cat(prior['RA'] ,prior['Dec'] ,'GAMA09_Ldust_prediction_results.fits',ID=prior['help_id'])\n",
    "prior160.prior_bkg(0.0,5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Divide by 1000 so that units are mJy\n",
    "prior100.set_prf(pacs100_psf[1].data[centre100-radius100:centre100+radius100+1,centre100-radius100:centre100+radius100+1]/1000.0,\n",
    "                pind100,pind100)\n",
    "prior160.set_prf(pacs160_psf[1].data[centre160-radius160:centre160+radius160+1,centre160-radius160:centre160+radius160+1]/1000.0,\n",
    "                pind160,pind160)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "----- There are 2165 tiles required for input catalogue and 18 large tiles\n",
      "writing total_bytes=781068219...\n",
      "writing bytes [0, 781068219)... done.\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,prior100.sra,prior100.sdec,unique=True)\n",
    "order_large=6\n",
    "tiles_large=moc_routines.get_HEALPix_pixels(order_large,prior100.sra,prior100.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='./'\n",
    "xidplus.io.pickle_dump({'priors':[prior100,prior160],'tiles':tiles,'order':order,'version':xidplus.io.git_version()},'Master_prior.pkl')\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": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "total 1525656\n",
      "drwxr-xr-x@ 3 pdh21  pdh21   102B 14 Feb 16:33 \u001b[34mdata\u001b[m\u001b[m/\n",
      "-rw-rw-r--  1 pdh21  pdh21    41K 15 Feb 11:16 XID+PACS_prior.ipynb\n",
      "-rw-rw-r--  1 pdh21  pdh21    18K 15 Feb 11:25 Tiles.pkl\n",
      "-rw-rw-r--  1 pdh21  pdh21   745M 15 Feb 11:25 Master_prior.pkl\n"
     ]
    }
   ],
   "source": [
    "ls -ltrh"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "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.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}