{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Final Processing\n",
    "The final processing stage requires:\n",
    "1. Quick validation of catalogues and Bayesian Pvalue maps\n",
    "2. Skewness level\n",
    "3. Adding flag to catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import seaborn as sns\n",
    "from astropy.table import Table\n",
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import pylab as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read tables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: UnitsWarning: 'degrees' did not parse as fits unit: At col 0, Unit 'degrees' not supported by the FITS standard.  [astropy.units.core]\n",
      "WARNING: UnitsWarning: 'mJy/Beam' did not parse as fits unit: At col 4, Unit 'Beam' not supported by the FITS standard. Did you mean beam? [astropy.units.core]\n"
     ]
    }
   ],
   "source": [
    "cat=Table.read('../data/CDFS-SWIRE/PACS/dmu26_XID+PACS_CDFS-SWIRE_20171019.fits')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<Table length=10>\n",
       "<table id=\"table4480558696\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>RA</th><th>Dec</th><th>F_PACS_100</th><th>FErr_PACS_100_u</th><th>FErr_PACS_100_l</th><th>F_PACS_160</th><th>FErr_PACS_160_u</th><th>FErr_PACS_160_l</th><th>Bkg_PACS_100</th><th>Bkg_PACS_160</th><th>Sig_conf_PACS_100</th><th>Sig_conf_PACS_160</th><th>Rhat_PACS_100</th><th>Rhat_PACS_160</th><th>n_eff_PACS_100</th><th>n_eff_PACS_160</th><th>Pval_res_100</th><th>Pval_res_160</th></tr></thead>\n",
       "<thead><tr><th></th><th>degrees</th><th>degrees</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy/Beam</th><th>mJy/Beam</th><th>mJy/Beam</th><th>mJy/Beam</th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>str27</th><th>float64</th><th>float64</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th></tr></thead>\n",
       "<tr><td>HELP_J033650.824-295620.783</td><td>54.21176639</td><td>-29.939106486</td><td>6.87275</td><td>16.2245</td><td>1.78597</td><td>13.6253</td><td>27.3919</td><td>4.32662</td><td>5.85037</td><td>4.14517</td><td>3.95679</td><td>0.289239</td><td>1.00055</td><td>1.00093</td><td>2000.0</td><td>2000.0</td><td>0.003</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033655.442-295449.596</td><td>54.2310102146</td><td>-29.9137765914</td><td>20.6963</td><td>36.8386</td><td>7.17254</td><td>52.5302</td><td>71.8413</td><td>33.7153</td><td>5.07407</td><td>4.70184</td><td>2.50966</td><td>0.273807</td><td>0.999138</td><td>1.00028</td><td>2000.0</td><td>1069.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033616.314-295548.181</td><td>54.0679762</td><td>-29.930050306</td><td>20.4869</td><td>36.9846</td><td>7.26307</td><td>10.0954</td><td>21.9491</td><td>3.04819</td><td>4.86817</td><td>3.09368</td><td>0.162121</td><td>0.150051</td><td>0.998687</td><td>0.998595</td><td>2000.0</td><td>2000.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033601.094-295546.425</td><td>54.0045569</td><td>-29.929562436</td><td>26.0394</td><td>45.0952</td><td>10.3994</td><td>35.1703</td><td>55.4868</td><td>17.5868</td><td>3.95185</td><td>4.14887</td><td>0.168648</td><td>0.160877</td><td>0.999148</td><td>0.99955</td><td>1100.0</td><td>1228.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033606.888-295518.496</td><td>54.0287013</td><td>-29.921804496</td><td>19.1804</td><td>37.2356</td><td>6.11049</td><td>12.3666</td><td>26.795</td><td>3.84751</td><td>4.58486</td><td>3.36094</td><td>0.137129</td><td>0.13689</td><td>1.00177</td><td>1.0003</td><td>2000.0</td><td>2000.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033607.350-295516.790</td><td>54.03062359</td><td>-29.921330436</td><td>37.0479</td><td>57.8534</td><td>18.3395</td><td>16.7424</td><td>31.4419</td><td>5.55697</td><td>4.58486</td><td>3.36094</td><td>0.137129</td><td>0.13689</td><td>0.998624</td><td>1.00164</td><td>2000.0</td><td>1629.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033610.856-295506.940</td><td>54.04523478</td><td>-29.918594406</td><td>38.9425</td><td>60.466</td><td>19.7534</td><td>39.9349</td><td>58.8555</td><td>20.5389</td><td>4.58486</td><td>3.36094</td><td>0.137129</td><td>0.13689</td><td>0.99846</td><td>1.00126</td><td>1430.0</td><td>891.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033555.160-295410.647</td><td>53.97983175</td><td>-29.902957406</td><td>140.717</td><td>159.889</td><td>121.065</td><td>98.4569</td><td>115.299</td><td>81.2531</td><td>4.326</td><td>3.08357</td><td>0.150315</td><td>0.130711</td><td>0.999733</td><td>0.999802</td><td>2000.0</td><td>2000.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033559.324-295522.307</td><td>53.99718296</td><td>-29.922863036</td><td>21.6159</td><td>40.2529</td><td>7.36929</td><td>29.3267</td><td>47.2794</td><td>13.6595</td><td>4.326</td><td>3.08357</td><td>0.150315</td><td>0.130711</td><td>1.00198</td><td>1.00001</td><td>1416.0</td><td>2000.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "<tr><td>HELP_J033600.747-295242.908</td><td>54.00311348</td><td>-29.878585636</td><td>24.9269</td><td>43.857</td><td>9.58364</td><td>63.3859</td><td>83.0581</td><td>44.418</td><td>3.7181</td><td>4.17204</td><td>0.132656</td><td>0.128657</td><td>1.00019</td><td>0.999564</td><td>2000.0</td><td>1347.0</td><td>0.0</td><td>0.0</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=10>\n",
       "          help_id                 RA      ... Pval_res_100 Pval_res_160\n",
       "                               degrees    ...                          \n",
       "           str27               float64    ...   float32      float32   \n",
       "--------------------------- ------------- ... ------------ ------------\n",
       "HELP_J033650.824-295620.783   54.21176639 ...        0.003          0.0\n",
       "HELP_J033655.442-295449.596 54.2310102146 ...          0.0          0.0\n",
       "HELP_J033616.314-295548.181    54.0679762 ...          0.0          0.0\n",
       "HELP_J033601.094-295546.425    54.0045569 ...          0.0          0.0\n",
       "HELP_J033606.888-295518.496    54.0287013 ...          0.0          0.0\n",
       "HELP_J033607.350-295516.790   54.03062359 ...          0.0          0.0\n",
       "HELP_J033610.856-295506.940   54.04523478 ...          0.0          0.0\n",
       "HELP_J033555.160-295410.647   53.97983175 ...          0.0          0.0\n",
       "HELP_J033559.324-295522.307   53.99718296 ...          0.0          0.0\n",
       "HELP_J033600.747-295242.908   54.00311348 ...          0.0          0.0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cat[0:10]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy import units as u\n",
    "cat['RA'].unit=u.degree\n",
    "cat['Dec'].unit=u.degree"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Look at Symmetry of PDFs to determine depth level of catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAa4AAAGoCAYAAAAerAGHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xt81NWd//HXXJIJIQkECMpFShVycL0gghXEFqSgrT/d\narGyUrVqrfa6ttvWWnuxunZru7a129qLrVZrS9Xq2rW7WhVBVJAqEBDaegLeKjcJSciVTJLJ/P6Y\nGZiEmW8myUxmvjPv5+PBQ2e+c/l8E5jPnPP9nM/xhMNhRERE3MKb7QBEREQGQolLRERcRYlLRERc\nRYlLRERcRYlLRERcxZ/tAIaqrq4lXFlZSmNje7ZDGRaFcq6Fcp6gc81XQz3XqqpyTxrDySt5MeLy\n+33ZDmHYFMq5Fsp5gs41XxXSuQ63vEhcIiJSOFw/VSjp8+zmXSk/duEpkzIYiYhIchpxiYiIq2jE\nJYOSaHSmUZiIDAclrgI0kClBEZFco6lCERFxFSUuERFxFSUuERFxFV3jkrRJdu1MRRsikk4acYmI\niKtoxJXHVD0oIvlIIy4REXEVJS4REXEVJS4REXEVXeOSjFO1oYikk0ZcIiLiKkpcIiLiKkpcIiLi\nKkpcIiLiKirOyBNabCwihUIjLhERcRUlLhERcRVNFUrWJJre1NouEemPRlwiIuIqSlwiIuIqSlwi\nIuIqSlwiIuIqSlwiIuIqqip0mT+/+CYtrR3ZDkNEJGuUuCSnxErky8tKeiVolcmLSIymCkVExFWU\nuERExFWUuERExFWUuERExFWUuERExFVUVSiuoIa8IhKjxJXDEn1Yl5eVZCESEZHcoalCERFxFSUu\nERFxFSUuERFxFSUuERFxFRVniGslKl4BVRuK5DuNuERExFU04soByUYOIiJyJI24RETEVTTikryj\na18i+U0jLhERcRUlLhERcRVNFQ4jFWFklxr1iuQHjbhERMRVNOLKEI2u3EGFHCLuo8QlkoCmFUVy\nlxLXAGgUVdgy9fsvLyuhpbUj5ccrgUqh84TD4WzHICIikjIVZ4iIiKsocYmIiKsocYmIiKsocYmI\niKsocYmIiKsocYmIiKsocYmIiKsocYmIiKsocYmIiKsocYmIiKu4vldhXV1LuLKylMbG9myHMiwK\n5VwL5TxB55qvhnquVVXlnlQfW1fXkne9+5zOPy9GXH6/L9shDJtCOddCOU/QuearQjrX4ZYXiUtE\nRApHxqYKjTFFwD3AVCAA3GqtfSzu+BeAq4G66F3XAtuBnwIzgSBwtbV2R6ZiFBER98nkiOtSoN5a\n+17gA8BP+hyfDVxurV0Y/WOBC4ASa+084Abg+xmMT0REXCiTiesPwDei/+8Buvscnw181RjzgjHm\nq9H7zgT+DGCtXQ/MyWB8IiLiQhmbKrTWtgIYY8qBh4Gv93nIA8CdQDPwqDHmPKACaIp7TMgY47fW\n9k16h1RWlgJQVVWevuBzXKGca6GcJ+hc89VwnWtlZWlBFYNktBzeGHMM8CjwU2vtirj7PcAd1tqm\n6O3/A2YRSWLxv2mvU9ICaGxsp6qqnLq6lrTHn4sK5VwL5TxB55qvhnquA0l6+bjEwOn8M1mccRTw\nFPBZa+0zfQ5XANuMMccDbcAiIoUcI4DzgYeMMXOBrZmKT0RE3CmTI64bgUrgG8aY2LWuXwIjrbV3\nGWNuBFYTqR58xlr7uDHGCywxxqwjcl3sygzGJyJ5Ztu2rfzoR7fj9/s47bS5XHXVNb2ONzc3ccst\n36CtrY1Ro0bxla98ncrKMRmLJxjs4JZbvkFjYyOlpaV87Ws3U1lZmbH3KxSecNjdC67r6lrCmn7I\nP4VynqBzTacrrljOt7/9PSZOnMSXv3wd11zzaaqrZxw6/pOf3EFFRQWXX34VL7/8F5555iluuOEb\nDq84eFVV5fz4xz+jra2Nj3/8WlaufJJt27by+c9/KdXnq3NGEq5v+SQiA/P443/i+eefpb29nQMH\nDnDllVezcOH7qanZyF13/RSfz8fEiZO4/vqvEQx2cNttt9La2sL+/XV8+MMXc+GFF/HZz15DZeUY\nmpub+eIXr+c737kFn89PT08PN910K0cddTQ//vEPeeWVzQAsWfIBLr74Em644QZCIdi7dw/19fu5\n8cZvYcwMli49j3e9aypTp76bf/3XLx6K9frrP097++HrN1OnHsuXvnRDwvNqa2ulq6uTSZMmA/Ce\n98xjw4aXeiWuN998nWuu+TQAJ588kx/+8HsA3H//vUyfXs3cuWf0+3OK2bnzbW677d97xbBkyQf4\n0Ic+fOj2K69sYfnyywGYO3c+99579wB+U5KMElcKgl0hmlqDjCoLEChKvXKnv+cN9nVFhurgwYP8\n8Id3cuBAI5/4xMc488wFfPe73+ZnP/sVlZVj+OUvf8bjj/8JY45n8eKzWbBgEfv31/HZz17DhRde\nBMDixeewYMFZPPLIQxx//Al8+tPXsWVLDW1traxd+zx79uzmrrvuJRQK8alPfZzZs08D4OijJ3D9\n9V/jscce5bHH/psvf/lG9u17h3vu+S2jRo3uFef3vndHyufU1tZGaenIQ7dLS0vZvXtXr8dMn254\n4YXnqK6ewQsvPEdHRwcAl112Rco/J78/8rE5efIx/OQnd/UbU1lZ2aF42tpaUz4fSU6Jy0Gop4cH\nV+2gpraOhuYgYyoCzKquYtmiafi8yZfA9fe8wb6uSLqccsqpeL1exowZS3l5Bfv311Ffv59vfCMy\nmgkGg5x22unMmzefhx5awZo1qyktHUl39+Ei3ylT3gXAeed9iN/97j6++MXPMXJkGdde+xneeusN\nZs48BY/Hg9/v54QTTuLNN18HIskDYPz4o9i6dQsAo0aNPiJpQf8jrkceeZDVqyO1X1//+s0cPHj4\nse3t7ZSV9a5Mu+yyK7jjjtv5zGc+wbx58znqqKMG9HM6cOAA48aNA1IbcY0cOZL29ra4eMoc309S\no8Tl4MFVO1i5Yeeh2/XNwUO3ly+uHvTzBvu6Iuli7asANDTU09bWRlXVeMaPH89tt/2AsrIyXnhh\nDSNGlPLAA7/lxBNP5sILL2LTpg28+OILh17DG/2S9cILa5g5cxZXXXUNTz/9Z373u/tYsGARjz/+\nGMuWfZTu7m62bXuFD37wPDZvfhmP58hLF94kX9j6G3EtXbqMpUuXHbrt9xexa9dOJk6cxEsvvciV\nV/Yuzti8uYbzz7+Ak06aybPPPsNJJ80c0M8pvrAilRHXSSfN5MUX1/JP/3Qi69evZebMWY6Pl9Qo\ncSUR7ApRU1uX8FhN7X6WLjiu1/RebNpvRMDv+Lzzz5ia9PgmW8f7Zk6kavQITR1KRjU01HPddZ+i\ntbWVL37xK/h8Pq677kt8+cvXEQ6HKS0dyTe+cTMej4cf/vB7PPPMU5SVleHz+ejs7Oz1WjNm/BO3\n3noT9913Nz09PXzuc/+GMTOoqdnItddeSVdXF4sWLcaYGfzv/2b2vL70pa9y881fp6enh9NOO50T\nTjgRgC984TN873t3MGXKu7j11psAGDeuiq9+NVKYkegaV7Kf00BceOFF3HrrTXzqUx+nqKiIm266\nNQ1nKaoqTGJfYztf/cV6Ev10vB74j2vmMr6y9Ihpv1FlxRxo7UzwrMjzvrjsFG5/YHPC14XIGgCn\nqcNCqUArlPOE4T/Xxx//E2+99Saf+tTnhu09Y3L19xobYcauw8HQf05pWICsqsIkNOJKYlRZgDEV\nAeqbg0ccqywvYVRZADhyWjBZ0oo9b/L4sqSvCxBGU4ciw23aNMPRRx+d7TAkRUpcSQSKfMycPo5V\nG3cdcWzm9LEEinyO04mJzKoeR3lpMbOqq3olu2QSTUmKDNW5556f7RByTqKkpZ9T7lLicpBsnBq7\nv6k1SEOSkRNAZVmAprYgleUlzKoex7JF0wAO/bemdj8NzR1Jpw0bWzpoag0yPtpIWERElLiSCnaF\n2Lx9f8Jjm7fXc9HCkON04tiKEm746Cz2NR5k8vgyykuLDx3rDoVZPHsy558xlabWIHf8YQsNLUdO\nMcZPSYqISIQSVxJOo6n4kVCyab/SEj+3/W5Tr3Va58+fykPP7ODVfzQeur+0pIi2jsQN8GdVR9aL\n7Gts1yJlEZEoJa4kUi3OiJ/2a2zpoLK8hNISP2/vO7xCPlZssXrTTkI99Lo/0euXFPs446SjCYfD\nfP2X63slv89erHUgIlLYlLiSCBT5ko6mZlWPOzT68Xm9LF9czdIFx1F34CCdXd387I9/Tfia8UnL\nycgSP+GeMKtqdh+6L5b8SkcUc8H8qQM+HxGRfKHE5SDRaCq+yCKmPdjFiqe38+pbDQmvVQ1UQ0uQ\nmiTX19Zv28MH33OMpg1FpGApcTmIH00laoYb6ulhxcrtrHtlD8HuFIdTKagYmXwR8/4DB1VpKCIF\nTYkrBYEi3xGJItTTwy33buh1LStdmhwWMY8bPUKVhiJS0NSKfJBWPF075KTl86bc0eWQuSdO0DSh\niBQ0Ja4UBbtC7GtsJ9gVinTMSHINKhVjygPMPeEoKkpTH/COrQiweM5krjr/hEG/r4hIPtBUYT8S\n7Z01ffIox56ETkaXFXPTlafRerCL9X99J6XneIDrLjqZyePL8fn0XUNECpsSVz8eeGY7z8T1K6xv\nDlL/t32Dfr2m1k4OBrt58uV/pPycyvIAVSrGEBEBNFV4hL5Tgmu37k3r63s88F+PvMILW/ak/Jz2\nYDePrHmNUE/6KhdFRNxKI66oZFOCHZ2htL5PTxh272/v/4FxOjpDhxZCX3fJ7LTGIyLiNkpcUX33\n1RrqlGAm1NTup6MzcV9DEZFCoalCGPC+WtlS39zB/gMHsx2GiEhWKXHR/75aueS7v3lZ17pEpKAp\ncXG4E3wiY8oDnHXqJIr9A18snAlv7mlhxcrt2Q5DRCRrlLgAv89DaUlRwmMzp43lrFMm4vHkRuIC\n2Fy7n2BXeotGRETcQsUZRAozErVvKhvh55XX6lkdt71ILmhsDarRrogUrIwlLmNMEXAPMBUIALda\nax+LO34J8HmgG9gKfNpa22OM2QQ0Rx/2hrX2ykzFCM6FGa0Hu2k9mHtVfJHBXzjbYYiIZEUmR1yX\nAvXW2suMMWOAzcBjAMaYEcCtwEnW2nZjzO+B84wxTwEea+3CDMbVi5sKM2LCYbjhF3/hmPFlfO3y\nUyn2a+AsIoUjk594fwAejv6/h8jIKiYInGGtja3E9QMdwEygNJrA/MCN1tr1Tm9SGZ0uq6oqH1SQ\n5aNGUFU5gn2N7iszf3tfK99bsZkfffGsbIeSEYP9nbqRzjU/Dde5VlaW4vcXzq4RGUtc1tpWAGNM\nOZEE9vW4Yz3AO9HjnwPKgKeBE4HbgV8B04EnjDHGWpt0vq6xsZ2qqnLq6loGHevJx43ttfjYTd7c\n08zrb9VTXlqc7VDSaqi/UzfRueanoZ7rQJJeY+PAuvG4gdP5Z7Sq0BhzDLAauN9au6LPMa8x5nZg\nCbDUWhsGaoHfWmvD1tpaoB6YkMkYAZYtmsbiOZMZGy2Jz6ECwn71hGFnBjazFBHJVRlLXMaYo4Cn\ngK9Ya+9J8JBfACXABXFThlcB348+fyJQAaTejXaQfF4vyxZNoyQQGYCGXVT34PXA+MoR7KxrZee+\nFpXJi0jey+Q1rhuBSuAbxphvRO/7JTAS2AB8HHgeWGWMAfgRcDdwrzHmBSJlc1c5TROm04qV29lV\n1zYcb5VWpSV+vvbL9XR2R7JtSbGXM06awCXvn47Pq2V6IpJ/MnmN6zrgOoeHJPtUXZ6BcBwFu0Js\nrnXe0bjY56EzlFtDsRHF3iPK9Ts6e1i1cRdej4fli6uzFJmISOboKzmRkvgDrc4l8bmWtACCncl7\nFm6ydZo2FJG8pMSFc6/CXObUarehJdJdQ0Qk3yhxAYEiH7Oqq7IdRlqVlfgZVea+ZCwi0h+1XIha\ntmgaABtfraMxD0YqM6aOIVBUOAsSRaRwaMQV5fN6Wb64mm9ddRqjRibuFO8myxdPY19ju65ziUje\n0Yirj/LSYipGBmhq68p2KIM2aVwp3/7NRhqag4ypCDBjSiWXLKmmNKBft4i4nz7J+gh2hWjvcG/S\nGj+6hF37D7d/qW8OsnbbXjbW7uPMkyeybNE0re8SEVfTJ1gfTa1B6l3WLT7evgMdCe/v6Oxh5Yad\nPLhqxzBHJCKSXkpcfYwqC1BSnL8/lhrtniySd57dvCvbIQyr/P2EHoBgV+hQIUOoJ0xXt9MKKXdr\nbOnQ+i4RcbWCvsYV6unhwVU7qKmtO1TIUOz3EcrfvEVleUDru0TE1Qo6cT24akevfbjcfG0rVR2d\n3fh9Ltq3RUSkj4KdKgx2haiprct2GMOurSPEipXbsx2GiMigFWziamoN0lAAI6xENqtAQ0RcrGAT\nl1sb66bDgTY14BXJN4VUWViwiStQ5GPGlMpsh5EVZSV+3t7XSkt7Z6+KShERNyjo4oxLllSzsXYf\nHQ77WuWjloPd3PnoNgB8Xgj1wNiKALOqq9RZQ0RyXkF/QpUG/Jx58sRsh5FVsdL/+uagOmuIiCsU\ndOKCyHYmi+dMZmxFCV4PBIoK+0eizhoikusKeqoQDm9nsnTBcTS1BhlVFuCh1TtYvalwLnTGi3XW\nGF9Zmu1QREQSKuzhRZxAkY/xlaX4fR48HvAW6BrdyvISddYQcalCqSxU4urjwVU7WLVxFz3hbEeS\nHbOqx6W8c7IqEkUkGwp+qjBeS3snG17dl+0wssLrgUlVZVy08Nh+H5uox6MqEkVkuChxcfiDeOOr\ndRxo7cx2OFnRE4a397Xy8LOvs3xxteNjE/V4jN3u77kiIkOlr8cc/iBuVDeJfqsKnXo8qiJRRIZD\nwSeuQm22m0x/+3U59XjUXl8i2ffs5l15X6RR8ImrkJvtJtJfVaFTj0dVJIrIcMjYNS5jTBFwDzAV\nCAC3Wmsfizt+PvBNoBu4x1r7S2OMF/gpMBMIAldbazPayiH2QVwIe3Glor+qwkCRj1nVVb2ucaX6\nXBGRdMjkiOtSoN5a+17gA8BPYgeiSe2HwNnAAuAaY8xRwAVAibV2HnAD8P0Mxgcc/iAudBWlxZx1\n6iSWLZrW72P7dhsZW1HC4jmTU3quiMhQZbKq8A/Aw9H/9xAZWcUcD+yw1jYCGGNeAN4HzAP+DGCt\nXW+MmZPB+A6JfeDW1O6nvrljON4y5zS3d7K5dh8+r6ffsvZE3UY00hKR4ZKxxGWtbQUwxpQTSWBf\njztcATTF3W4BRiW4P2SM8Vtr45NeL5XR1kRVVeVDive6S2bT0dmNfauRr/983ZBey60aW7tYuWEn\npSOK+cQFJ6X0nMkZjGeov1M30bnmp+E615GlxXj7fNnM559zSonLGDMfOAn4NXC6tfa5FJ93DPAo\n8FNr7Yq4Q81A/E+1HDiQ4H6vU9ICaGxsp6qqnLq6llRC6tfoET4CRV6CXYW11Um8J9e/ydlzJlMa\nGJ5lfsGu0BEjt3T+TnOdzjU/DfVcB5J42tqPXH/6h6dfBWDhKZMGHUM2OZ1/v59MxpjriFx7mkRk\n+u8Xxpi7rbW39/O8o4CngM9aa5/pc/jvwHRjzBiglcg04e1AGDgfeMgYMxfY2l986fbH598o6KQF\nEOzq4dv3beCWq9+T0U4YTh04RESSSeVT6QrgHKDNWlsPnAZclcLzbgQqgW8YY56N/vmoMeYaa20X\n8G/Ak8CLRKoKdxEZnXUYY9YRKd74woDPaAi0puuwPQ3trHi6NqPvEVv4Xd8cJIz2BBOR1KQyFxSy\n1nYaY2K3O4B+2yNYa68DrnM4/ifgT33u6wE+mUJMGaE1Xb1tqq3j4kXTjyi8SDS1N1D9deDo6HSc\nIRaRApZK4lpjjLkdGGmMuQC4Bug79ZcXtKart6a2Ln77pOWKc2fg83rT2ly3vw4cjc1BNdIUkYRS\n+bT5MrAd2AJcDjwOfCmTQWWL1nQdae22vYem7tI5tddfB47KJMdERJJ+qTXGTIm7+UT0T8xE4B+Z\nCiqbli2aRqgnzAtbdqF+sREbX63jnNOOcZzaW7rguAFNG/bXgaOk2E9h1J6JZFZ830K3Vhj25TQb\ns4ZIlV/8XsCx22Gg/42bXCbU08Pvn9nOi9v2KGnFaWwN8u/3baQ5QcktHG6uOz66pi5V8Qu/G1s6\nqCwvYVb1OFUVioijpInLWvvu4QwkF8R2P5YjJUtaMPjmuurAISKD4TRV+C1r7beMMb8mMsLqxVqb\nSkl8zupbGRfsCrHJFubux0M11Oa6gSLfgEdrIlK4nKYKN0b/++wwxDFsklXGnTVrEg0thbn78UCN\nLiumua1TU3sikhVOU4WxNVYTrbXfiT9mjPmPjEaVQcm2nX/1rUa8nsgW9pKc1wMzp43j7NOOYUxF\niab2RFwkXwo1nKYKbwPGA/9sjJne5zlziXTGcBWnRa8769qGORp36gnDms27KfJHrk+JiAw3p3Vc\njxCpLGyL/jf250ng/2U+tPRTZ4z02WTrCKr0UkSywGmq8GXgZWPMH621Tcke5ybqjJE+DS3BQZXA\ni4gMVSpddS4wxnyfSMNciK7jsta67uKG06JX5+d56eru0fWvOF4PjBimbU9EROKl8slzE7DQWrst\n08EMh0SLXktL/Ly9rzXpczq7eo5cD1DgesJwMNhNeWlxtkMRkQKTSuLalS9JCxIvevX7PKxYuZ01\nNbsSjqoqywNAWOXyccZWBBhVFkhLp3gRGX7xFYYxbqk0TCVxbTTGPExkU8iO2J3W2t9kLKph0HfR\n62VnGwiHWV2z+4jHjhxRRN2Bg8MZXs4bEfDz8LM72Lx9/5A7xYuIDEQqiWsU0ALMi7svDLg6cSWy\nfEk1Pp93QNOIhWpnXVuvJQSx9XCAyuRFJKP6TVzW2isBjDGV1trGzIeUPX2nEUcE/Nxy78vZDstV\namrrBtwpXkRkIPpNXMaYmcCDQKkxZi7wHHCxtXZTpoPLltg04r7GdpXOD1B9s8rkRSSzUrkY8WPg\nQqDeWrsb+BTw84xGlSNGlQUo8nn6f6AcojJ5Ecm0VBJXqbX277Eb1tqngYLZntajvDUgsTJ5EXGf\nZzfv6vUnV6WSuBqi04VhAGPMR4GGjEaVI5pag3R2awXXQP34v7fS2a3kJSKZkUri+hRwJ3CCMeYA\n8HngkxmNKkeMKgswtqJgBpdps6uujW//Jm8vgYpIlvWbuKy1rwFLgTHAFOASa63NdGC5INYiSgZu\nV10rLQ67JouIDFa/icsY86/AE9baNiL9Cv9kjLkm45HliGWLpvH+2ZPweXWxayB6wrBT699EJANS\nKf+6BjgdwFr7ljFmNvAX4K5MBpYrfF4vFy2cFtkxWS2fUub1wOTxZUmPq1WUSO7L1bZQqSSuIiB+\nMVMnFFbP2abWII1KWgMyqaosYQPeUE8PD67aEfkioFZRIjIIqSSuPwKrjDEPRW9/GPifzIWUe7SP\n18BMqhrJ1y4/NeGxB1ft6LWtjFpFichApZK4biSyAHkB0AX8l7X2j6m8uDHmdOC71tqFcfcdDTwQ\n97BTgBustT83xmwCmqP3vxFrN5Vtg93Hq1Bdee4MDrR0MqrM02saMNgVoqa2LuFzamr3q1WUiKQk\nlcT1srX2VODhgbywMeZ64DKgLf5+a+1eYGH0MfOAbwO/NMaUAJ74JJdLzp8/lVWbdtLTk+1Ict9/\n/GYjPWGoLCvi+KljWb5kOqWBIppagzQkGbU2tnSoVZSIpCSViwrvGGPea4wZ6IKm14hMKyZkjPEQ\naSf1KWttCJhJpB/iU8aYVdG+iDkh1NPDbfdvUtJKUWxPs8bWLtZt28uX7lzHipW1lJUWMybJurjK\n8hJGlWnNnIj0L5UR1xxgDYAxJgx4gLC11nFOx1r7iDFmqsNDzgf+GrcmrB24HfgVMB14whhjrLWO\nLRgqo9/Qq6rK+z+TQfrZI1vY09CesdfPdx2dIVZu2ElxsZ/5Myfx2POvH/GY+TMnMnni6F73ZfJ3\nmmt0rvlpuM51ZGkx3mEqbtq4oz7h/R+YN3VY3h9S29YkUytwLwV+FHe7FthhrQ0DtcaYemAC8LbT\nizQ2tlNVVU5dXUtGggx2hVj3ypGbS8rA/fnFN3nfKRNYNHsSW7bXH9rzbFb1OM6fN6XX7zCTv9Nc\no3PNT0M914EkvbYcWOyf7t+r0/mnsq1JMfAlwACfI9Ly6TZr7VB/UnOAdXG3rwJOAj5tjJkIVAB7\nhvgeQ9bUGuRAa/b/UuSDnjA8W7OHxXMmc+snTtc6LhEZlFTGlncCZcBsoBuYBtw90DcyxiyPddww\nxlQBzdHRVczdwGhjzAtE9v+6qr9pwuGgfoXpt8lGKgvHV5YqaYnIgKVyjWu2tfZUY8wHrbXtxpiP\nAVtTeXFr7ZvA3Oj/r4i7v45IGXz8YzuB5akGPlxUCp9+DS1Bfvuk5YpzZ6S06FhdNkQkXiqJKxyd\nLoyNjsZRYJ0zli2aBsALr+yhozOU8DGBIi+zq6tY99d3hjM011q7bS8jSvyOi47VZUNEEknlX/8d\nwEpggjHmDmAD8MOMRpVjfF4vyxdXc/tn5nPGiUcTKD7yxxbs6mHdX99Bn6epq6ndT7Ar8RcBONxl\no745SJjDXTYeXLXjiMcGu0Lsa2x3fD0RyZzh3IQylarC+40xG4GziCS68621r2Q0qhxVGvBz9Xn/\nxLL2adx0z0sJizbCWuuVMqdFxx2d3Sl12dCoTKTwpLKtSRFwNvABIsnr9Oji4YJ1MNhNU5JKw4Ka\nQx0iv8/DiEDi706Nzf132YCBjcpEJD+k8pX0V8AZRLYxuY9IAiuoqcK+Yk13ZWg6u8Pc/OuXWLGy\nllCftiSVFcl/xrEuG/31PtS0oUh+SiVxnW6tvdha+ydr7f8AHyEyAitY2hk5fRpaOhOOkEqK/Ul/\nxrOqxxEo8qXU+1BE8k8qVYVvG2OmWWtjnyxHAZm98uYCsUrDTbaOhhZ9QA5Vou7wsZ9xTe3+Xl02\nYvc7bTej3oci2ZXJTShT3UhyizHmOSILkM8E9hhjVgFYaxelJRKXiVUavm/mRG66+yVd2xqiRIUa\nsZ/x0gW+p2qhAAAgAElEQVTHJVzH5bTGLjYqE5H8k0riuqnP7dszEYhbVY0eoU0m08BphBQo8iXd\n7qS/UZmI5J9UyuHXDEcgbqXOGukx2BFSf6MyEck/qYy4pB99v/VXjCxWY94UeT0wqaqMixYeO6TX\ncRqViUh+UeJKg77f+kcE/Nxy78uaPkxBTxje3tfKw8++7tj+SUQkJpVtTd4NnEdkc8ceYAfwJ2vt\nWxmOzXXiv/XPmFLJ2m17sxyReySqKhSR/BJfaTiUCsOk67iMMROMMQ8CvwemEklYFpgCPGSMedAY\nM3nQ75znLllSTUmCnoaSWH1zBw3NHWl7PfUuFMlfTiOu24CbrbV/S3TQGDMT+A5wWSYCc7vSgJ8z\nT56ooo0BeGrD23zsnBlDeg31LhTJf0kTl7X2Y05PtNZuQUnLUaJS7aIiD3vrD2Y5stz0l7++w78s\nmj6k6cJY78KYWO9CQNfQRPJEKte4jgWuJbIP16HmutbaqzIYV15IVKq9fecBfvDglmyHlpM6OkPU\nNbYzeXz5oJ7fX+9CXUMTyQ+pVBU+QmQ/rudR8/NBiS/aeNdR5Xg9kWo6OVIoPPgfTF1je9JKTqct\nVETEXVJJXB5r7ZczHkmBKC8tprTET+vB7myHkpOe/Ms/+NgHjx/Qc+KvayWj3oUiuSVWYTiY6sJU\nEtc6Y8yFwP9Ya7VN4hAFu0IU+Qp6OzNH6/+2j+07m5g/cxLnz5uSUkFF3+taiah3oUj+SJq4jDE9\nRKYGPcAngbAxhujtsLVWnwKD0NQa5EBrV7bDyGn1zUEee/512g929ltQ4XRdC2BMeYBTTZV6F4rk\nEaeqQqc1XppzGSSnrTikt1QKKpz25PJ44DMXnsjIEUV0h8L4VA0vkhf6/adsjHmxz20vsCFjEeU5\nbUKZuoYUNoN02o06UOTjzke38tVfrOfrv1yfcKdlEXEfp1HVquh04enGmJ7YH6CDSAcNGaRli6Zx\n1qnp2VAtn3mAJ1/6x6Fkk6gbhtMXgY7OEA0tnYQ5vJ6r707LIpJdz27elXDTSSdOU4WLAIwxP7LW\nXjfE2CSOz+vlsrMNze1BNr66P9vh5KyeMKyu2Y3H68Hr8STthnHkQu8AbR1ddHQeObrSei4R93Os\nKjTGnEN0I0ljzOXAe4AN1tp7Mx9a/rvyg8ezufZ5Qpq9crRu6146Og+Psvp2w+i70LuzK8RN97yc\n8LW0nkvE/ZymCu8AbgRKjDH/DnwU+CvwYWPMfw1TfHmtNFDEglmaMuxPfNKKV1O7/4hpw/GVpVRV\nlia97qX1XCLu51SccTawyFq7l8i2Jv9srf0ZcCGwZDiCKwSXvH86Z506Ca3sGrjGJMUbTte9tJ5L\nxP2cpgrbgfHAHuAdYCQQjP43pbYPxpjTge9aaxf2uf8LwNVAbAHOtcB24KfAzOj7XG2tzfsr6bHr\nXYTDrK7Zne1wclLA7yXYfeR8qtPoKVGD41nV47SeSyQPOCWum4GXjTEPAK8Ca4wxK4FzgO/198LG\nmOuJdI9vS3B4NnC5tXZj3OM/DJRYa+cZY+YC3wc+lPKZuNzyJdW0B0P85W/vZDuUnFLs93BK9Tj+\n8rd9RxxzGj0lanCskZZI7hpICyhP2KGpaXT34wuBaUSS3F7gf621L/X3wsaYpcArwP3W2rl9jv2d\nyPWyo4H/s9Z+xxjzA+Ala+0D0cfsstb2ewbd3aGw358fH0gdnd1cetMTBBNUwxUqnxdCPTAi4AM8\nBDu7GTd6BCceN45rLjiRkSOKsx2iSKakfAXh8bWvh715st/cB+ZNjf1v0vN3rCq01r4B/GAwb26t\nfcQYMzXJ4QeAO4Fm4FFjzHlABdAU95iQMcZvrXWclmxsbKeqqpy6upbBhJlz5p5wNGs0ZXhIrOLy\nYDBShDFhTCnBrm5Wb3ibLbX78mKTyHz6+9sfnevAnp+qtvbOQb9Pron9zJzOf9j/tRtjPMAd1tr9\n1tpO4P+AWUSSWHyk3v6SVj66dEk1x4wvy3YYOWtPQ7sWFYsUOKcmu7uBoxIcGmqT3QpgmzHmeCLX\nvxYB9wAjgPOBh6LXuLYO8vVdzef18s0r5rDi6VrWbt1DZ7c27uqPFhWLFBanqcLTgFXAhdbavw31\njYwxy4Eya+1dxpgbgdVEqgefsdY+Hu2BuMQYs45IcrxyqO/pVj6vl4sXTWfLjv00tOTPFECmaFGx\nSGFxavm0K1q2fjPwkcG8uLX2TWBu9P9XxN1/P3B/n8f2ENk+RYh2PVfSSokWFYsUlv6KMx4HHh+m\nWCTOqLIAY8qLlbxSoEXFIoVlQMUZxphrMhWI9BYo8nGqGZ/tMHLa2IoAi+dM1qJikQLjOOJK4JPA\nXZkIRI60bNE0esLhI5rMSsR1F53M5PGplwyLSH4YaDm8WuoNI5/Xy6VLDD/83JncsPyUbIeTc56p\n2aWNIUUK0EAT188zEoU4ChT5qJ4yhtKAruPEW1OzW2u4RPJIKu2eIIWpQmNMEbAYGAccjO7LhbX2\nN0MJUAampb1T04UJbLJ1WsMlUmBSucb1B2AC8Hcgtho2DChxDaOd+1rp0VrkIzS0BHl9VxPHThql\n5CVSIFJJXDOstTMyHok4mjy+DI8HHHoiF6z/fGAzY8qLOdWM79W3MNgVUmd4kTyUSuJ6zRgzxVr7\nj4xHI0mVlxYzadxIdtYl2iVGGlo6WblhJz3hMJe8fzoPrtpBTW0dDc1BxlQE8qIZr4hEOPUqXE1k\nSnA8sNUYs4W4DSSttYsyH57E+/rHZnP9z9bT3KZFycms27qXcBhWb9p16L5YM16A5YursxWaiKSJ\n04jrW8MVhKSm2O/nNzedww9WbGD9tnfoTLArcKHr6AxRU1uX8Jia8YrkplSrCWOcehWuATDG/Nha\n+7n4Y8aY+4A1gwlQhsbn83LFB47n4rOmc/OvX6LuQEe2Q8o5Ta2JR6RqxiuSH5ymCn8FHAvMMcac\n0Oc5ozMdmDjzeT30qMwwoSK/N+FoVM14RfKD01ThrcBU4EdEOsTHdBMpjZcsamoN0tAczHYYOSnZ\nFGppiR+/T81fRNwuaYlVdEuS54DPElnHdRQQAjZYaxuGJTpJalRZgDEVGj0MxNv7WtVpQyQPJE1c\nxpgzAAt8E/h/RHYn/hZQa4x5/7BEJ0kFinzMqq7KdhiuU1O7n2CXOpCIuJnTVOFdwLnW2u3xdxpj\npgGPAidlMjDpX2w7j5ra/TS2dFDk9xLsUqWhExVoiOSOgVYTxjglLl/fpBX1OuoSnxN8Xi/LF1ez\ndMFxNLUG6ewKcdM9L6OSjeRUoCHifk6J63+NMX8CHgD2RO87Gvgo2hU5pwSKfIyvLCXYFWJMRYB6\nFW0kNSLgU4GGiMs5FWd8GbgXeD9wPfAVIl3if22tvX5YopMB0XWv/u2sa+OWezdoHy8RF3PsVWit\nfQR4JLq1yfFAp7X21WGJTAYldt3ruZpddIY0aZjI2/taWfF0LZedo97RIm7kVFX4RPS/JxJZt3Uf\n8KAxZmv0PslBseteX1o+K9uh5LSa7aouFMmGhadMOvRnsJxaZR8V/e8dwOestbOstTOJrOu6e9Dv\nKMPimPHllBSrJ18yTa2d1DW2s6+xXQlMxGVS2dZklLX2idgNa+0aY4xqiXNcoMjH/JOO5pmNu/p/\ncAEKFPv40cOvaNsTERdy+lc6zRjzMyBojPkEgDGm0hjzJQ5XGUoO+5f3T+esUyfhVRHdETo6Q9Q3\nBwlzeNsTddUQcQenxHUC8DTwPJGWTxAphT8duCKzYUk6+LxeLjvbsOCUidkOJWeMKQ8knUJVVw0R\nd3CaKqyz1v438N+xO6y1PwF+ErttjCmx1mpfjRy3fEk1Pp+Xja/uozHJlh+F4F/OOpaxo0u589Ft\nCY+rq4aIOzglrt8ZY/4MPGCtbYk/YIwpBy4nsq7rwmQvYIw5HfiutXZhn/svAT5PpNP8VuDT1toe\nY8wmoDn6sDestVcO8HwkgVil4fyTjubmX2/IdjhZ88Dq1x2Pq6uGSOYMpYqwL6fE9RHgU8DLxpgD\nwE4iiWYqMJbIdicfSfZkY8z1wGVAW5/7RxDZMuUka227Meb3wHnGmKcAT98kJ+nz3Obd2Q4hp82Y\nom3mRNzAaQfkHuBO4E5jzExgOtADvGat3ZLCa78GfBi4v8/9QeAMa217XAwdwEygNJrA/MCN1tr1\nAzkZSS7YFeKV1+qzHUbO8nlh7ba9vPqPRlUYiuQ4TzicuLuCMWYtsB14AnjKWts40Bc3xkwlMtU4\nN8nxzwHnRv+cCMwFfkUkST4BGGttt9N7dHeHwn6/1iv1Z8/+Nq69bSVJft3Sxz+/91g+cYE2QJCs\nSrke+PG1r4e9Of5F6wPzpg70KUnP32nENT+6hckHgd9Ep/ieBZ6w1m4caATxjDFe4HtANbDUWhs2\nxtQCO6y1YSJ7ftUTqWZ82+m1Ghvbqaoqp66uxelheWOw5xrqCjGmXA14U7V2y24++J5jCBQd/lIU\n7ArR1BpkVFmg1/1Dpb+/+Wmo51pVVZ7yY9vac7/oaqA/C6fz769X4Q7gx8CPjTElwELgcmPMj6y1\nZw4oit5+QWTK8ILolCTAVUT2+Pq0MWYiUIHWi6VNrAHvyg07sx2KK8RXGIZ6enhw1Q5qauu0YFkk\nB/TbOcMY827gPA5f49oO/GCgb2SMWQ6UARuAjxNZH7bKGAORQo+7gXuNMS8AYeCq/qYJZWD6bjxZ\nPqKIpvauLEeVm+IrDB9ctaNXwo8tWAZYvrg6K/GJ5Lp0VhH2lTRxGWMmEOlT+C5gLbAD6ALeDTxk\njHkT+KK1NulXeGvtm0SuW2GtXRF3KNnX1OWphy4D1XfjyREBP7fc+7KmDxM4edpYAkU+Wto72fhq\nXcLH1NTuZ+mC49I6bSgi/XMacd0G3Gyt/Vuig9FKw+8QKXkXF4ltPAlo+jCJzbX72LGziabWDprb\nEw/8tWBZJDucijM+5vTEaEm8kpbLLVs0jfaObtZt25vtUHJKY2sXja3O06hasCySHalc4zoWuBYY\nR1x5orX2qgzGJcPE5/Vy2TmGv73ZwIECbgc1GLOqx2maUCQLUtnW5BFgJZFiCq0CykOxisPVm7QF\nSiqKfR7OOHkCZ82aRLArNCwl8yJukMmCjHipJC6PtfbLGY9Esmr54uns2NnE2/tasx1KzusMhXl+\ny26erdnN2Ghp/EULj+XhZ18f1pL5XEqSuRSL5L9UEtc6Y8yFwP/ErbmSPOPzevnmFXNY8XQtm2rr\naGpTmbyTUPRfQqw03v7jQK+kP5SS+f6SQC6tK8ulWKRwOJXD9xCZGvQAnwTC0TVXHiBsrdXXqjwT\nud41A4DVNWrIOxA7k4xUUymZ7+jsZl9jO2WlRfzx+Tf6TQLJ1pUd7Ojm0nPMsI54tMZNssGpqjDp\n1yVjjEqp8pSa8Q5Osou/TiXzsdHKK6/VU9d4kECxj47OwxtZJkoCwa4QNbWJ15Wt3baXv7/VwKlm\n/BHJLtgVoq6xHTweqkaPSEtyc4pFa9wkk1KpKnzRWjsv7raXSPcLdSDNQ02tQRq0IDlt4kvm+yaP\nR9a81mu0Ep+04m2ydYeSQH+/n4aWTlZu2Eko1MNl58wg1NPD75/Zzrqte+jojMxvlhT7mH/S0fzL\n+6cPaTrPKZaGZq1xk8xxmipcRaQ3YWzaMKYbeCyzYUm2jCoLMKZCzXjTZVb1OPw+D/c9+Xde3LaX\nzq7I2CxQ5MWTYu/vhpbgoSSQ6u9nzebd4PHg9cCqjb2rRTs6QzyzcRcej2dI03lOsXg88OTLb7N8\n8dCSo+Su4aogTMRpqnARQLSh7nXDF5Jkk5rxps/cE47i/Pnv5pu/eok9De29jgW7Uq9z8npgRCDy\nTzXV309PGFZv2kWgKHnSqKmNjOQAx2KQZMUiTrHE3t/ndU6OqkaUwXCcKjTGzAD2G2N+TqTB7m7S\nsK2J5LZli6YRCvXwbM1uLdwbgvV/fYe//PWdIf8Me8KRUdfBYDejygJxzZLr+h15OSXIhuYg9/35\nVba/fSBhMUgqFYOxvytrNu+mJ8GJJrvWpWpEGQqnqcJPA9cAfwBejt49AfiVMea31trvD0N8kgWx\n6sJQOMxzm7WzzFCkK/Hfdv9GOrt7en3AL11wHPc/aQfdritMJLnG9C0GSVYxGAr1cM57phwaJZ3z\nnik8m6QKNVlxSjaqETW6yx9OI67PA6dYa3vNcRhjfgBsApS48txlZxve2N2iRck5INgdGTnFPuBb\n27v42AdncOW5M/D74Lkt6es1+cIrezh37pSkFYNrNkcWX8eS6Pnzpx5RERlTXOSjrLQocg5dIfbs\nb+Nge+ewViNqdJd/nBJXF1CU4P4R0WOS52KLkn/7VC3PbdlNWPOGOWP9396hZnsd80+ewIJTJqc1\ncXV0hvjtk7VJKwZjU4KxJPrX1xuSVkR2dIZ49LnX8Xg8kcTREmTUyOKkfTEz0XFfa83yj1Pi+jZQ\nY4x5hsM7EU8AFgFfy3Rgkht8Xi8f+8AMvB4tSs41wa4eVm3cxbMZ6DH5xp6mlKtL+xae9LV2695e\nic2pmXO6O+5rrVl6ZLOCMJGk4+Toxo9nEmmu2w4cjP7/e621DwxPeJIrli+p5qxTc+svr0QkKooY\nqsbWLop86ZlGSzYaSyTdHfed1prFRnfiPk7dMSZZa3dba38DPAq0AaWAvp4UIJ/Xy2VnGxbMmpDt\nUGSY7G08yDHjyxhTnvlGOSXFPhbNnnSoYjJVwa4Q+xrbCXYlTo6xtWaJaD8193KaKvwTcKox5iLg\nB0S2N/EC1xtjrrfWPjwcAUpuuXSJwYsnafmz5JfW9s5Da8gyqaMzRDgM3aEw9U3tR1T+9e064vd5\nUiq4cFprpv3U3CuVv5E3EJkefAvAGPMfwGpAiasA+bxex/JnyS+NrZ00DtMGo6s37WLd1t0Eu8K9\ntot5aPVrR7SsGje6hJ372g4916ng4vC6t/00tnRQWV7CrOpxAx7dSe5IJXG1AW/H3d6HNpQsaGoL\nJYNR5PfQ1e380RGMtsRKtl0MREZn8UkrXqKCC5/Xy/LF1SxdcFzCdVyDXd+Vr+vCcq0QIxGnxDXK\nGPM3ItODPwA+b4yZDdxEpEhDCpTaQslg9Je0Ekm2XUwyTuX0gSJfr/sHu75L68Kyz6lX4XHGmDHA\naUBJ9O6jgSeBnw9DbJLDli2aRld3pNWPSKYMNNUNpOBisOu7tC4s+xy/HlhrG6y1T1pr/yd61yRr\n7Z3W2tTrWyUv+bxePnj6lGyHIdJLrOAiWbVh7P6Wfrp3JKtS7G9dWLLnSXoNtFzok8BdmQhE3GdU\nWYCxutYlWXLM+DLaO7p7FVxctPBY7n/Ksrl2PwdaD0/jXbTwWB5+9vVD03ujypy7d9Q1tlNc5Ds0\neotdy0plXZj2IMu8gSauFHcQkkKga12SLe87eQLnnD6FYr+XfY0HmTy+jNISP7fcu6FXMUeyIg+n\n7h3FRT5+9PArNDQHCRT7gDAdnT2MrQhw8nFjkxYmaV3Y8Blo4tK1LeklvtS4oaUDwio5lczbtL2O\n517Zg9cT6RwytiLAiICfnXWJqw131aVe5NHRGTrU7SO+60d9c5DVNbs5ZnxZwsTlxnVhbqggTKTf\nxGWMKQIWA+OAg8aYywGiHTX6e+7pwHettQv73H8+8E0iuynfY639pTHGC/wUmAkEgauttTsGdjoy\n3PqWGj/50j/U01AyrvVgN9C74W/kYyOxdC6WbzvYxVmnTuKVHfVaF5YlqYy4/kCkue7fOfxlOgw4\nJi5jzPXAZUTWgcXfXwT8kEi1Yhuw1hjzGDAfKLHWzjPGzCWybcqHUj8VyaZYqfHyJdV4vB5e2LKH\nzu7Ud/kVcYsDrUHOOe0YLj5rWl6u43KDVBLXDGvtjEG89mvAh4H7+9x/PLDDWtsIYIx5AXgfMA/4\nM4C1dr0xZs4g3lOyzOf1cukSw4fmv5tv3v0Xmtq0A47kl9i1rL7rwmT4pJK4XjPGTLHW/mMgL2yt\nfcQYMzXBoQqgKe52CzAqwf0hY4zfWtvt9D6V0b84VVXlAwnP1dxwrlXAglOP4bHnX892KCJpNX/m\nRCZPHJ3SY4fr3+rI0mK8g1j87IbPkkSSJi5jzGoiU4Ljga3GmC1ErkkBYK1dNMj3bAbif1rlwIEE\n93v7S1oAjY3tVFWVU1fXMshw3MVN53r+vCns3NvMpu37sx2KyJCVFPs48+QJnD9vSkr/Bof6b3Ug\nSaWtfXD9JHP5s8Tp/J1GXN9KeyQRfwemR7tytBKZJrydSJI8H3goeo1ra4beX4aJz+vlYx+cQc32\nF1RpKDnH64H3nTIRr9fDlu2RQovRZQHag90J9xAbWeJn6YLjXNfWya2Vg06cWj6tATDG/Nha+7n4\nY8aY+4A1A3kjY8xyoMxae5cx5t+ItI7yEqkq3GWMeRRYYoxZR2S92JUDOxXJReWlxUweX3ZEo1SR\noRo9spiTp49l5nHj+PEjA/+eu+CUiVx2TuTy/UcWRhrmdnb3cNPdLyV8fGNLUAuMc4TTVOGvgGOB\nOcaYE/o8J6UJXmvtm8Dc6P+viLv/T0T2+4p/bA+RzhySZ752+an8+30b2ZVkjY3IQE0YU8rXPjab\n0kARwa4QJcW+lHZa9nhgTILy9VihRbArpAXGLuA0VXgrMBX4EXBz3P3dRKb7RFJS7Pfz7x8/nZb2\nTuzbjXjwcFTlCL55z8vZDk1cak9DO398/g2WL64mUORj/klH88zGXY7PGVMe4PMXz6Rq9Iik5eup\nbjyZr1uauIXTVOGbwJtEFgQDYIy5xFr7+8yHJfmovLSYOeYoAPY1tmc5GnG7+L23LnzfcayN22wy\nkVNNFZOryvp9XaeNJ7WlSW5wmiq8PMHdt0QXEKfUOUMkGTXolaGKb2rb2t5J0CFpnXHi0Sl3tnDa\neHLFytqc39IkH4sx+nL6ivBJIt0rFgJnRf+UR/+7MNOBSX6LTcmIDFb8NafYrtyJjK0IcNk5ZsAj\noth1r/jpQW1pkhucfpPvBX5GpBDjemvtlcA/rLVXWmuvGpboJG8Fu0KcNWsSZ506ibEVJdp2QAbs\n5OPGHEoqTl+EZlVXpeU6VCpbmsjwcLrGFQK+aYyZDzxmjPk2avwtQ5ToGsHJx43lfadM5L/+sIXG\nVrWIcru5J4zn7X1t7Nnfltbmtn0tnnNMr9tO16bSITaqU8Vh9vXb8slau9YYcw7wEyJdNEQGLdG2\n56trduPzeTl+6ljWbdubxehkKLxeKPJ6+ctf9zGmIsB7Z05kTnUVR48t5WBniDse2kxDy+A6PPQ1\ntqKEMRUlve7re21qRMDPwWA33aEwvjTUTaRacSiZl9J+XNbaZmPM14H3GWOOs9a+luG4JA/1d43g\na5fPZlNtXUrrcSR3FPs9+H0+2oPdBHsiBRL1zUHWbN5Nkd/LCceOBeCU6ipW9VOyDpHWSiNL/DS2\nBCkuSrw+yylR+H0eVm7cmZHKv0yP6iQ1TlWF7wfuBdqBrwD/BawHbjbGfNVa+9CwRCh5o79rBJ1d\noZTW40j6eD2R+f/RIyNb2Q92Zq89mLitaHzJeqrXMTu7Qtx46akUF/koKy3ij8+/MaBEkWhUn67K\nP6eKw2wqhErCeE4jrv8EPgiUAauBk621240xVcDTgBKXDEgq1wh0EXV4LThlIue8ZwqjygJ4PGFu\nuXcDu/cPbI1dZ3fy31pDtGhhVFmAzSk2W64sL6EqrppvIImiv1F9LIkOlbY0yS6ncbPfWrsNeAlo\nstZuB7DW1pHiFKNIPOfKr3EAbFEn+WExtqKExXMms3xJ9aGS72K/n1uvnssPPnsGp04fx+iRRUN+\nn9EjA4wqCziOtvtKNA3YtzQ9GVX+FQanBPSKMWYFMJLInlz/CfwauBDYMRzBSf5xukZQ39SR8odb\nPkm1z95gxaYDx5SXcPK0sSyePZkxFSVJk8DoshI+u/Rkgl0hGpo7WLlxJ1u276ehJfHvxin+U6JJ\nyGm07fVAOAxjKoZ+vUiVf4XBKXFdCVwG+Ihc6/oWkenBzcA1mQ5M8pPTNQKnD51iv8dxSsqtKkr9\nzJ85ic7ObjbZOhpb+6+68wATxpWmPKUXPx04kGmyQJGPCWNHctnZhovPmhZJYhve5pXXGnp96QiH\nwwmvSx4zvozli6cfeq1kFXkLZk3inNOOScv1IlX+FQanxOWz1t4Td/tr0T+HGGNKrLUdGYlM8lqi\nawROHzqnzTiKddv25t01sOb2bp548S0Wz5nMf1w7j98+aVmbwpKAz1x4EqtrdvUauZ4yfSzdoR62\n7KinqbWz1whmqNV0h5LYOTOOaDAb6unB4/FQU7ufhuYORpUVM2v6OJYvqe71vssWTaN0RDFrt+w+\nYrSdzj5/qvzLf55wOPFHgTHmEeDPwAPW2pY+x8qBy4HF1toLMx6lg7q6lrCbdgUeqnw/19gC5Vde\nq2f/gYOHPnQueO+x3HT3X/K2t+HYihK+ecUcWg92sXLjTl7ZsT/puY6tKOHWT5xOoMiXsEt5tjqX\np/K+VVXl7Nx9YFjiy3YH9zTsgJxyQ5m6upZ8+07neP5OI66PAJ8CXjbGHAB2EtnSZCowlsh2Jx9J\nX5gih6cSr106gtferO/1oZNsNBYz/8SjKS7ysrpm93CFmzb1zR18656XOdB6uJvIwa4Q67e9c8Rj\n46e8ko1cs1Hxlur7Dld8qvzLX04tn3qAO4E7jTEzgelAD/CatXbLMMUnBaqk2H/Eh05sqmeTraOh\nJYjXAz3hSBPV2AJTAJ/PS03tfuqbh38W+5jxZbQd7EpayOCkMVrxFusm8v7Zk1g8Z7KmvET6SDpV\n6BaaKsxPTucZmwKKtfRJNBUUe8wTf3mLNZv3ZCzOkmIfnV2hXkmlOxRO+VqVk9iUIJBTi12HolD+\n/oKmCodqsFOFIjkpfgqovLTY8TGXnm3o6g4PuQfi2IoApSVFtB3s4kBrMO7a27tpbe/qlVR8Xrji\n3ADXWKsAAAx/SURBVBmMKPEfGi2NLgvQHuweUNl7/H5TmvISOUyJS/Kaz+vlsnMMf3+rkcYUp+/i\nR1F91z0luuBfGjhyoW6isv9H1ryW8BrdiICPg8EjE5rWHYkkpsQleS9Q5GO2SV7Y0Xe6L9EoKv61\nBjL6iX98sjLtkpIi/veFN454rtYdiSSmxCUFYdmiaYTDYdZu3Xtouq6k2McZJx3Nh993HK3tnf2O\nooYq2eLrMWNG0tHRpSIMkRQpcUlB8Hm9fHSJ4aKF06g7cBDC4V6NXEsDw/dPoe+ozefLzY7jIrlK\niUsKSqDIx+SqsmyHkZDWHYmkJn19VkRERIaBEpeIiLiKEpeIiLiKEpeIiLhKxoozjDFe4KfATCAI\nXG2t3RE9djTwQNzDTwFusNb+3BizCWiO3v+GtfbKTMUoIiLuk8mqwguAEmvtPGPMXOD7wIcArLV7\ngYUAxph5wLeBXxpjSgCPtXZhBuMSEREXy1iTXWPMD4CXrLUPRG/vstZO6vMYD/Ay8FFrrTXGnA78\nBniLSFK90Vq73ul9urtDYb9fa15EJO+k3GQ3Tz8Hs9JktwJoirsdMsb4rbXdcfedD/zVWmujt9uB\n24FfEdlG5QljjOnznF4aG9vVcToPFcp5gs41X6WhO3zKj21sbB/0++Qqp/PPZOJqBuLf2ZsgAV1K\nZEPKmFpgh7U2DNQaY+qBCcDbGYxTRERcJJNVhWuBcwGi17i2JnjMHGBd3O2riFwLwxgzkcioLXOb\nKYmIiOtkcsT1KLDEGLOOyFzllcaY5UCZtfYuY0wV0BwdXcXcDdxrjHkBCANXOU0TiohI4clY4rLW\n9gCf7HP3q3HH64iUwcc/pxNYnqmYRETE/bQAWUREXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJ\nS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0RE\nXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJS0REXEWJ\nS0QkzYJdIfbsbyPYFcp2KHnJn+0ARETyRainhwdX7aCmto6GliBjygPMqq5i2aJp+LwaJ6SLEpeI\nSJo8uGoHKzfsPHS7vjl46PbyxdXZCivvZCxxGWO8wE+BmUAQuNpauyPu+BeAq4G66F3XAtudniMi\nkquCXSFqausSHqup3c/SBccRKPINc1T5KZNj1wuAEmvtPOAG4Pt9js8GLrfWLoz+sSk8R0QkJzW1\nBmloDiY81tjSQVNr4mMycJlMXGcCfwaw1q4H5vQ5Phv4qjHmBWPMV1N8johIThpVFmBMRSDhscry\nEkaVJT4mA5fJa1wVQFPc7ZAxxm+t7Y7efgC4E2gGHjXGnJfCc45QWVkKQFVVeTpjz2mFcq6Fcp6g\nc80X82dO4rHnX09w/0QmTxydsfetrCzF7y+cachMJq5mIP5vqDeWgIwxHuAOa21T9Pb/AbOcnpNM\nY2M7VVXl1NW1pDX4XFUo51oo5wk613yyePZEHl/3Bt2h8KH7/D4Pi2dPHPB5DyTBNza2D+i13cDp\n/DM5VbgWOBfAGDMX2Bp3rALYZowpiyaxRcDGfp4jIpLTbvzF+l5JC6A7FObGX6zPUkT5KZMjrkeB\nJcaYdYAHuNIYsxwos9beZYy5EVhNpHrwGWvt49FKxF7PyWB8IiJpU990kNaDiSeIWg92U990kLGj\nRgxzVPkpY4nLWtsDfLLP3a/GHb8fuD+F54iI5Dz7jwP9Hj/jJCWudNBSbhGRNDBTnIsv+jsuqVPi\nEhFJg7GjRuD3eRIe8/s8miZMIyUuEZE0CHaF8BBOeMzjCavhbhopcYmIpEHdgYMky01d3ZHjkh5K\nXCIiadDZ5bjktN/jkjolLhGRNCguci7S7u+4pE6JS0QkDapGjyBQlPgjNVDkpWq0ijPSRYlLRCQN\nAkU+zjx5QsJjZ548QVuapJHGriIiafIv75+Ox+NJuAOypI8Sl4hImvi8XpYvrmbpguPwFRcR6uzS\nSCsDNFUoIpJmgSIfE8aNVNLKECUuERFxFSUuERFxFSUuERFxFSUuEZH/396dh+pVnHEc/0YNJmRR\n6x82ophCzYNSjWLjQrRYsdQtLgWhpCbmSmJCqmCqCG4JisVu1mJdwEBUXGOrlStuLVGiUWrbgGJI\n+jOKS01EtLWt1q1Jrn/MXDi+vVvSm/vemff3gQs5M+fkzNyH5Hln5rxnrChOXGZmVhQnLjMzK4oT\nl5mZFcWJy8zMiuLEZWZmRXHiMjOzoozp6el7q2kzM7PRyCMuMzMrihOXmZkVxYnLzMyK4sRlZmZF\nceIyM7OiOHGZmVlRnLjMzKwou7W7AUMVEbsAtwDTgc+A+ZJebdQvAeYD7+WihZI04g0dRhFxFPBT\nSce3lM8ClgJbgBWSlrehecNqgL5WE9eIGAusAKYCuwPXSupu1FcT1yH0taa47gosBwLoARZJWteo\nryauo0UxiQs4Exgn6ZiIOBq4HjijUX8EMFfS2ra0bphFxKXAHOA/LeVjgRuAGbnuuYjolvTuyLdy\nePTX16ymuJ4D/F3SnIj4CvAi0A1VxrXfvmY1xXUWgKSZEXE88GPy/00VxnVUKGmq8FjgCQBJfwS+\n2VJ/BHBZRKyJiMtGunE7wWvA9/ooPwh4VdIHkj4H1gDfGtGWDb/++gp1xfU3wFX5z2NIn8B71RbX\ngfoKFcVV0sPA+fnwAOCfjera4joqlJS4JgP/ahxvjYjmiPF+YBFwAnBsRJw2ko0bbpIeBP7bR1Xr\n7+FDYI8RadROMkBfoaK4SvpI0ocRMQn4LXBlo7qquA7SV6gorgCStkTEncCvgXsaVVXFdbQoKXH9\nG5jUON5F0haAiBgD/ErS+/lTzaPA4W1o40ho/T1M4suf8KpRY1wjYn/gaeAuSfc2qqqLa399rTGu\nAJLOBaYByyNiQi6uLq6jQUlrXM+R5pIfyGtcLzfqJgPrIuIg0jzyCaSF4RptAA7M6wYfkaYdftHe\nJu00VcU1IvYBfg9cIGlVS3VVcR2kr7XFdQ6wn6TrgI+BbfkHKovraFFS4vod8J2IeJ40Z94VEbOB\niZJui4jLSZ/uPgNWSXqsjW0ddi19/RHwJGnEvELSpva2bnhVHNfLgb2AqyKid/1nOTChwrgO1tea\n4voQcHtEPAOMBS4CzoqIjvj32g7e1sTMzIpS0hqXmZmZE5eZmZXFicvMzIrixGVmZkVx4jIzs6KU\n9Di8dbiImAq8AqxvqZol6W99nD8P+CXwVi4aD6wGFje+vL43sAm4QtL1LdefSnqseyKwK+krGcsk\nbYuIPYGbgUPz6ZuACyVtHEI/JgPPA6dJeiOXnZjbOh5YKenKXH4Y6THyPYBnSC9wbX19kllH8YjL\nSrNZ0mEtP/+TtBq6e88DDia9O+6HjfrZwCPA+fmNDgBExEnATUCXpOmkl6ROB67Op1wHrJN0iKRD\ngDuBlYM1Pr8Ffw3pDQu9ZeNJX8A9I7dvRkScnKvvJiXEaaTvLy4Y7B5mtXPiso4haStppDOtUdxF\nGjl9Dny7UX4FcLWkV/K1nwCLSSM2gK8C4/J2O5CS1rIhNGMBKXFubpQdCWyU9HoeTd0NnB0RBwDj\n80ulAe4Azh7CPcyq5qlCK82+EfFi4/geST8fyoV5WvBk4Cf5eDowBXiWlHgWAU/l0w8HXmheL+lt\n4O18eC3wMLA4Ip4C/kBKOAOSND/f+0t9At5pHL8D7DdAuVlHc+Ky0mzO035DdXpOdGNIMwwPAffl\nui7gAUlbI2IlsDQi9sl7JW3L1/RJ0tqI+BowEzgRuBhYGBHH7MAaVF/36e/+2/ooM+soTlxWu25J\n81oL8wZ/PwC2RETvhqQ9wHmk9au/kPZ8W9+4Zhppe45zSbtxL5G0GlgdEdcAG0kjtT9vZxs3kaYe\ne00hTSX2V27W0bzGZZ1qFvCepCmSpkqaCiwEFuSHNH4GLIuIAwEiYiL5CUVJPaQHPS5prHHtS/og\n+NoOtOWFdIv4et4GfjbwuKQ3gU8jYmY+by7w+I501qwmTlzWqbqAW1vK7gXGAd+V9ATpAY2VEfES\n8CdgLbA0n/t94BvA6xGxnrQx4mxJ/9jehkj6FJgHPEga4f2VtPkipFHhDRGxAZgA3Li9f79Zbfx2\neDMzK4rXuKxoEbGEtObUarOkU9rQnuNI27f35RRJXqMy+z95xGVmZkXxGpeZmRXFicvMzIrixGVm\nZkVx4jIzs6I4cZmZWVG+AJAMxnvBWRqsAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x115f3a978>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "skew=(cat['FErr_PACS_100_u']-cat['F_PACS_100'])/(cat['F_PACS_100']-cat['FErr_PACS_100_l'])\n",
    "skew.name='(84th-50th)/(50th-16th) percentile'\n",
    "g=sns.jointplot(x=np.log10(cat['F_PACS_100']),y=skew, kind='scatter')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For 100 $\\mathrm{\\mu m}$ depth is ~ 30mJy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAGoCAYAAAATsnHAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xt81NWd//HXzCSZJCSBBAKCqCiQg/UCKFZbtSjFut3V\nrpZWHmWrVbdqL7q22+r+6nrt2rXb1dZua7urq11rS6XV6mrX1oogBe9AEKl6Algv3CQkIVcySWbm\n98fMhEmYmcwkmcx3vvN+Ph4+ZL7fuZzvDMxnzjmf8zmecDiMiIiI03hz3QAREZFEFKBERMSRFKBE\nRMSRFKBERMSRFKBERMSRinLdgHQ1NrZnnG5YXV1OS0tXNprjGG6/RrdfH7j/Gt1+fTCya6ytrfSM\ncnNcw9U9qKIiX66bkHVuv0a3Xx+4/xrdfn1QGNeYC64OUCIikr/yZohPsu+5TTsPOXbWvMNz0BIR\nEQWogpQoEImIOI0ClKSULJipZyUi2aY5KBERcSQFKBERcSQN8blYNuealFAhItmmHpSIiDiSApSI\niDiSApSIiDiS5qBk1CglXURGk3pQIiLiSOpBuYSqQ4iI26gHJSIijqQAJSIijqQhPsk6JU+IyHCo\nByUiIo6kACUiIo6kACUiIo6kACUiIo6kJIk8MzjhoLKilPaO7hy1RkQkexSgJGe0ZYeIpKIhPhER\ncSQFKBERcSQFKBERcSQFKBERcSQlSYijJMtSVPKESOFRD0pERBxJAUpERBxJQ3wOpQ0IRaTQqQcl\nIiKOpB6U5AVVnRApPOpBiYiIIylAiYiIIylAiYiIIylAiYiIIylAiYiIIymLzwG05ml4kr1vyu4T\ncQf1oERExJEUoERExJEUoERExJE0ByWuo7kpEXdQD0pERBxJPagxpGw9EZH0qQclIiKOpB6UFAxV\nRBfJL+pBiYiIIylAiYiII2mIL0uUEJEflJIu4lzqQYmIiCOpByWSgBIqRHJPAWqENJRXODQcKDK2\nFKASUNCRTIz070tlRSntHd2HHFfgk0LnCYfDuW6DiIjIIZQkISIijqQAJSIijqQAJSIijqQAJSIi\njqQAJSIijqQAJSIijqQAJSIijqQAJSIijqQAJSIijqQAJSIijpQ3tfgaG9szrslUXV1OS0tXNprj\nGG6/RrdfH7j/Gt1+fTCya6ytrfSke9/hfA/mg2Tvgat7UEVFvlw3Ievcfo1uvz5w/zW6/fqgMK4x\nF1wdoEREJH8pQImIiCMpQImIiCMpQImIiCMpQImIiCMpQImIiCMpQImIiCPlzUJdESlMW7a8zg9/\neCdFRT5OOeU0Lr/8ygHn29pa+fa3b6Kzs5Px48fzT/90I9XVNaxf/wr33fdTioqKqK6u5sYbv01p\naWnO2imZUw9KRBztzjvv4NZbv8NPfnI/b7yxhYaGtwac//nPf8aJJ87jpz+9nyVLlvJf/3UPAHfd\n9V3uuONO7rnnPqZPP5Inn3w8p+2UzBVMDyrQG6S1I8D4Cj/+4sxWfY/ksSJj4amnnmTt2ufo6upi\n//79XHbZFznrrI9TX7+Be+/9CT6fj2nTDuf66/+ZQKCb7373djo62tm3r5FPf/oiLrzwM1x99ZVU\nV9fQ1tbGN75xPXfc8W18viJCoRC33HI7U6Ycxo9+9AM2b94EwDnn/BUXXfQ5vvOdWykuLmbPnt00\nNe3jhhtuxZg5LFlyHkcdNYNjjzVcccU1/W29/vqv0dV1sCzQjBnH8M1v/r+E19XZ2UFvbw+HHz4d\ngA9/+COsX/8KdXVz+u/zzjtvc+WVXwHgxBPn8oMffA+AH/3oXmpqJgIQDAYpKSkB4Oqrr+THP753\nwOtcffWVHHXUDN599x0AbrvtX5k4cVL/+UcfXcHq1c8OeMyNN36bww47DICOjqHbKZlzfYAKhkKs\nWLWN+oZGmtsC1FT5mV9Xy9JFs/B5U3cgR/JYkbF24MABfvCDe9i/v4UrrvgCZ5yxkH/7t+/w05/+\nN9XVNdx330956qknMeZYFi/+BAsXLmLfvkauvvpKLrzwMwAsXnwuCxeezaOP/ppjjz2Or3zlWl57\nrZ7Ozg6ef34tu3fv4t57/4dgMMiXv/z3nHzyKQAcdthUrr/+n3niicd44onfct11N7B37wc88MAv\nmDXrCBob2/vb+b3v3Z32NXV2dlJePq7/dnl5Obt27Rxwn9mzDevW/Ym6ujmsW/cnuru7AZg0KRJg\n1qxZxcaN6/niF78EcEhwijn++BO57rob+O1vf8NDD/2Mr33tuv5zS5YsZcmSpUnb2dHRMWQ7JXOu\nD1ArVm1j5fod/beb2gKsXL+Dru4+Lj7XpOwRJXsswLLFddlrtMgwzJt3El6vl5qaiVRWVrFvXyNN\nTfu46aZI7yQQCHDKKafykY+czq9/vZw1a1ZTXj6Ovr6+/uc48sijADjvvL/ll798kG984xrGjavg\nqqu+yrvv/oW5c+fh8XgoKiriuONO4J133gYiQQJg8uQpvP76awCMHz+B8eMnHNLOoXpQ8b2VG2+8\njQMHDt63q6uLiorKAc938cWXcvfdd/LVr17BRz5yOlOmTOk/t2LFL3nuuWe5664f4ff7U75/sWB7\nwgknsm7dmgHnhupBVVRUDNlOyVzWApQxphh4AJgB+IHbrbVPxJ3/OvBFoDF66CprrR3NNnT39FHf\n0Jjw3Atb9mDfa0naIwr0BpM+tr5hH0sWztRwnziKtZE5j+bmJjo7O6mtnczkyZP57ne/T0VFBevW\nraGsrJyHH/4Fxx9/Ihde+Bk2blzPiy+u638Ob/Tfwbp1a5g7dz6XX34lzzzzB375ywdZuHARTz31\nBEuX/h19fX1s2bKZT37yPOAFPJ5Di1F7k4wyDNWDGtxbKSoqZufOHUybdjivvPIil102MPlg06Z6\nzj//Ak44YS7PPfcsJ5wwF4AHH7wfa9/i7rt/gt8/dHKEtW8yefIUNm9+jaOPPiZlmwarqKgYsp2S\nuWz2oD4PNFlrLzbG1ACbgCfizp8MXGKt3ZCtBrS0BWhuCyQ9n6pH1NqR/LEt7d20dgSYXF0+eo0V\nGaHm5iauvfbLdHR08I1v/BM+n49rr/0m1113LeFwmPLycdx00214PB5+8IPv8eyzf6SiogKfz0dP\nT8+A55oz50PcfvstPPjg/YRCIa655h8xZg719Ru46qrL6O3tZdGixRiT/TmWb37zW9x2242EQiFO\nOeVUjjvueAC+/vWv8r3v3c2RRx7F7bffAsCkSbV861s30dzcxM9+dh91dXP4xjf+AYCPf/wT/XNt\niYb5nnrqd6xYsZzS0lJuuunbo9ZOGT5POJyd7UWMMRWAx1rbboyZCLxqrT0m7vybwJ+Bw4D/s9be\nker5hrMPSuX4Mr50x0qaUgQpgIlVpdx+xakDekSB3iA33vdSwscOvn8ukyhqaysHjO+7jduvD0bn\nGp966kneffcdvvzla4a+8xhz2mf4wx/exbXXfmPAsauvvpLrrruBo46aMaznHMk1aj+o5O9B1npQ\n1toOAGNMJfAIcOOguzwM3AO0AY8ZY86z1v4u2fNVV5cPa8+V0+cezhNr3055n5b2bnwlxdROGjfg\neLLHnj53GtOnTSAYDPHAk3/mpS27adx/gNoJZZx2/FQuP/84fL6xS6KorXX3WLfbrw9Gfo2VlaWU\nl5c49r1yUru++tWrDmlPSUkR1dXlI2rnWFzjcL8H81XWelAAxpgjgMeAn1hrH4g77gGqrLWt0dtf\nASZaa/8l2XMN55dDbW0lez5oZcWqbWy0jTS3J+5JJepBQXwW3z5a2ruprvQz58hqPndOHeX+Ipav\nbBiQRBGzeMH0MUuicNqv09Hm9usD91+j268P1IMaqTHvQRljpgB/BK621j476HQVsMUYcyzQCSwi\nklAx6nxeL8sW17Fk4Ux+8bTl+S17DrnP/LpJCYfmYo+94MyjWf7MVt56t5kXtuzhrfdaOHHWJF7b\nqiQKEZFsyWaSxA1ANXCTMeam6LH7gHHW2nuNMTcAq4EA8Ky19qkstgV/sY9L/3oOZaVFcT2iUubX\nTWLpolkpH/v42r/wQlxga2oLsHpj8jUOSqIQERm5bM5BXQtcm+L8Q8BD2Xr9ROJ7U+kmNaRKN/d6\nIJSgw11dWcr4ioNrLlSJQkQkc65fqJuIv9iXdu8mVbp5ouAEYI6MLE5MVonigjOPpqOrVwFLRCSF\nggxQiSTr5ZT5i5hQ4aelI1G6uZ8TZ05k8/Zmmtu68ZdEHvdidBFweWkx7+/t6L9/bN3Vus27CfQE\nVTpJRCSFgg9QyXo5nznrGB557m3qGxoTBieAOUdWc9Gi2Vy0iEMSMJraAknXX3X3BPvvo9JJIiKJ\nFXyASlZvz763f0DvJ56/2Es4DM9HM/rmHDWBzW83D7sNyvoTETlUQY8rpUqA2NmYODh5PRDoDdHT\nFwIiAe351z+gvat32O2IZf2JiMhBBR2gmtu6kw7DJUuASHZ8JAZn/YmISIEHqJXr3096zpv22u7U\ninxDP1GyhcIiIoWsYANUoDfI5u1NSc+PVk+pL5j6iXxeD6FQiGAoNDovKCLiEgUboFKtb0qkprIk\nK+0IhsKs2riLFau2ZeX5RUTyVcEGqPEVfmqq0pv3mVBRwlc+fUJW21Pf0EigN5jV1xARyScFG6D8\nxT7m19Wmdd/Wzh5KfN6s9aIAmtsDyuQTEYlTsAEKYOmiWSxeMJ3qITLoJozzU1tdzrw0A9pw+Iu9\nVJRnLwCKiOSbgg5QPq+XpYtmUVaaOoNuXjTLbpQS+xLq7gnx+BAbK4qIFJKCDVCB3iA79rZz3xNv\nsGtfV9L7Ta8dx7LFswn0Btm0dV9W21TfsE/zUCIiUQVX6igYCvGrZ7fywuu76e4ZOrX7yxccj8/r\npam1K6Osv+Fo1j5SIiL9Cq4HtWLVNlZt2JlWcBo/roSaqtLInyv8jK/I7hyRB3j6lfe0JkpEhAIL\nUIHeIBvt3rTvX1FePKDCw/Taimw0q18oDKvrD66JCvQG2dvSpWE/ESlIBTXE19oRoLm9J+37dwf6\n6Ar08fjayLYbyer2jbaNdi/BUJjXtjbS3N5DTWUJJ5nJ2jdKRApKQQWo8RV+aipL0g5SLe0BfvG0\n5aU3PshyywZqbu9h9cadA26vXL+DUDjM588x/ccDvUF27+sk2BtULT8RcZ2CClAAc46q4YW4jQVT\nKS7y8sqbYxucUnnh9T189qxZFPk8BzdZbA9QU6mdeUXEfQoiQA3eNddf7KUvGCIYzUXweT0EE1SH\nDfQ6K1mhuydIY0sXf9q8O+Emi6CdeUXEPQri53Zs19ymtgBhIoEnGILTPjSFb//9h/n3r3yEqTX5\nkdrd2xdKusmi1lGJiJu4vgeVatfchvf3U1ZaxOZt+8YsAWKknokG2kRatI5KRFzE9QEq1bYaze2B\nAckI+eClNz6gtMSbcB2XduYVETdx/RBfqm01RmvX3LGXuOHamVdE3MT1ASrVthqjtWvuWOvpDfLR\n4w9jYlUpXg9MrCpl8YLpLF00K9dNExEZNa4f4gP6v7jrG/bR0t5NdWUpJ86soX7rPvZ3pL9w1ymq\nK0u5+NzIeihfSTHBnl71nETEdQoiQPm8XpYtrmPJwpm0dgQYX+GnyOdh647WvAxQ8UN5lVV+tr/T\nwfgKv4KUiMs9t2knZ807PNfNGDMFEaBi/MW+/gy35Ssb2NHYmeMWZcbrhY/NncrZ8w/vL8G0eXsT\njS0HqKnSYl0RcZeCClAQSTtvbOlKmnruZKEQvPD6B6yp341/UCafFuuKiNsUTICKryaRL2ueEunp\niwSlZNuF1DfsY8nCmRruE5G8VzBjQfHVJNwstlhXRCTfFUSASlVNwm2qK/1arCsirlAQASpVNQm3\nKS8tTji8p80PRSTfFMQcVKyahNuH9wDaO3to7+qhsjyyPf3gSu7K9hORfFEQ31Cpqkm4zf7OHm55\n4BWWr2zoD07xldxj2X6xbeVFRJyqIHpQcGg1iQkVfjq7ex2359No2N8R2YE3GAyxeXtTwvso209E\nnK5gAtTgahI9fSFuuf+VXDcrq+q37qM1SaUMbc0hIk5XEEN88fzFvkiWWzictMq5W7R29DAhSUaf\ntuYQEacrmB4UQFegj18908Bb77VEtn4vcXd8rq70M3f2pIR7XmlrDpH8VEj1+AoiQMWSBdZt3jWg\nAkOyagxuMeeoapYtno3P6xlQyX1+3SRtzSEijlcQAWr5Mw2srt+V62aMuc+cdUzCSu7qOYlIPnD1\nGFcwGOKhP1rWbCq84ARwIHBwUW6skruCk4jkC1cHqAee/DOrN+4ccufc/N36PbWV69/PdRNERIbN\ntQEq0BvkpS2707pvvm79PpTXtjWptJGICz236dDEJzdybYBq7QjQuP9ArpuRU83tAVo7AqrDJyJ5\nybVJEuMr/NROKGNvS+EGKZ8Hfv/ye2x5u0l1+EQk77j2W8pf7OO046fmuhk5FQzDmk27VIdPRPKS\nawMUwOXnH8fiBdOZWFWK1xNZuCqwbvNuugK9uW6GiEhKrg5QPl9kDdBtf38Kpx13GGG3ZkNkqLsn\nyPJntua6GSIiKbl2Dire42v/wgtb9uS6GY7y1rstBHqDWhclkqdimXxuLnvk6h4UFNZ275nY3xHJ\n8BMRcSrXB6ihtnsvKXb9W5DQhAq/qpmLiKO5/ts5tt17ItUVfr571WmcdvyUMW5V7tUdNUHDeyLi\naK4PUKm2ez95Ti0TKkr5wrlzCi7Dr6TIpfWdRMQ1XB+gILLde3y6+cSqUhYvmN6/5YS/2Me82ZNy\n3Mqx9fKf96qyhIg4WlpZfMaY04ETgJ8Bp1pr/5TVVo2yVFtOxPaKem1rYSVSBHpDNLZ0MX1yZa6b\nIiIj4OZsviF7UMaYa4HbgX8EKoD/MsZ8M9sNy4ZEW06sWLWNlet30Nzek8OW5caKVVu1YFdEHCud\nIb5LgXOBTmttE3AKcHk2GzVWAr1BNtq9uW5Gzvz5nf1c+8N1/PIZSzDk7t2FRST/pBOggtba+O5F\nN+CKyYvWjkBB9pziBUNhnt2wk1/8sSHXTRERGSCdALXGGHMnMM4YcwHwBPBsdps1Nsr8Ra7drDBT\nazbt4qb7X6anry/XTRERAdJLkrgOuAJ4DbgEeAr4z6EeZIwpBh4AZgB+4HZr7RNx588Hbgb6gAes\ntfdl2viRiCVHqDzfQTsbO/nOzzdy2+UfznVTRCRDgzcxdEPSRNIAZYw5Mu7m76P/xUwD3hviuT8P\nNFlrLzbG1ACbiPS+YsHrB0TmszqB540xT1hrP8j8EoZnxaptqs+XwM7GDppaDxAMhQdkO4qIjLVU\nPag1QBiIHwSL3Q4Dxwzx3L8BHon+2UOkpxRzLLDNWtsCYIxZB3ws+piEqqvLKSrK/MuytvbQNOru\nnj42b2/K+LkKQSgM33loI62dAWonlHHa8VO5/Pzj8Plyt2Qu0WfoNm6/RrdfH4zNNY4rL8Gb5maj\nbnjPkwYoa+3RI3lia20HgDGmkkigujHudBXQGne7HRif6vlaWroybkNtbSWNje2HHN/b0lXQO+0O\nZX+0iOzelgM8sfZtug70sGxxXU7akuwzdBO3X6Pbrw9Gdo2ZBJLOrvSTuvLpPU/2HqQa4rvVWnur\nMeZnRHpMA1hrh0w1N8YcATwG/MRauzzuVBsQ36JKYP9Qzzdaxlf4mVBRwv6Ows7gS1d9wz6WLJyp\n4T4RGVOphvg2RP//3HCe2BgzBfgjcLW1dnDW35vA7OjcVAeR4b07h/M6w+Ev9jF/9iRW1+8aq5fM\nay3t3bR2BJhcXZ7rpohIAUk1xPdk9I/TrLV3xJ8zxvxrGs99A1AN3GSMuSl67D5gnLX2XmPMPwJP\nE0l1f8BauzPJ82TFsnPq2Lazjff3dozly+al6spSbc0hkmcGZ/VB/mX2pRri+y4wGfiUMWb2oMec\nRiQAJWWtvRa4NsX5J4Enk53PNp/Xy82XLmD5yq28+uYHdBzQ+p9k5tdN0vCeiIy5VOkgjxLJ5OuM\n/j/239PA32S/adnn83q5+BOGWy49Ba3XPZTHA0dMruAzZx1DoDfI3pYuVUAXkTGTaojvVeBVY8zj\n1trWZPdzg4njyxhXVqRe1CDhMLy/t4PbH9zAgUAfzW0Baqr8zK+rZemiWfjSTHcVERmOdCpJXGCM\nuYvIfBJE10FZa/N+zCfQG6S1I0CZv4hin/pQyexo7Oz/c1NbgJXrdwDkLPVcRApDOgHqFuAsa+2W\nbDdmrMTKHNU3NNLcFqBqXDGtndp2IhNKPReRbEsnQO10U3CCg3tAxSg4Za65TannIvkmUWYfODe7\nL50AtcEY8wiRNU3dsYPW2p9nrVVZFOgNUt9QWLvnZoPHA0+/+j7LFs/WXJSIZEU6AWo8kVJEH4k7\nFgbyMkC1dgRobgvkuhl5LxSG1Rt34vN6NBclIlkxZICy1l4GYIypjhV3zWfjK/zUVPlpUpAaFZqL\nEpFsGXJsxhgz1xjzFvCaMWaaMWabMeakMWhbVviLfcyvq811M1wjVgZJRGS0pTN58CPgQiJ7O+0C\nvkwaGxY62dJFs1i8YDoTq0rxeqCm0s/hteO0u+4wqAySiGRLOnNQ5dbaN40xAFhrn4luAZ+3fF4v\nyxbXsWThTFo7Av0b8937xBZeemNvrpuXV06cNVHDeyJ5zqm78abTg2o2xswluuWGMebvgOastmqM\n+It9TK4u7/+C/fy5cygtUUZaJk6YUc2Ove0qgSQioy6dHtSXgQeB44wx+4GtRLZzd51yfxEfPWEq\nqzaMaWH1vPYfv40skSst8fLRE6byuY8r7VxERseQ3yTW2u3AEqAGOBL4nLXWZrthuaJpqOHp7gmx\nasNOVqzaluumiIhLpJPF9w/A7621nUTq8T1pjLky6y3LgUBvkE1b9+W6GXlto92r4T4RGRXpDPFd\nCZwKYK191xhzMvAycG82G5YLWsQ7cs3tPSqBJJLnnFISKZ3JgmIg/lu7h2jChNvEFvHKyPiUry8i\noyCdAPU4sMoYc7Ux5moiNfn+N7vNyg0t4h0d//rQRpavbCAYCuW6KSKSx9IZ4ruByELdhUAv8B/W\n2sez2qocWrpoFhAp4dPS3o3H4yEYcmWHMWtaOrRnlIiMXDoB6lVr7UnAI9lujBPEL+J96GnLC1v2\n5LpJeUt1+kRkJNIZ4vvAGHOmMabgJmfse3lfGzenVKdPREYinR7UAmANgDEmjIu2fE9FGX3pKfZC\naYmP9u5DU8tVp0/EXRJl92Uzsy+d7TYKMmtA23KkJwRMqCqjvbvjkHPz6yYNGN4L9AYH1D4UEUll\nyABljCkBvgkY4Brga8B3rbU9WW5bTsUy+uK3hpdDBUPw/t4OptaU09MXpKU9QHVlKfPrJnHBmcew\nt6WLivISHl/7NvUNjTS3Baip8jO/rpali2apLJKIJJXOEN89QCNwMtAHzALuBy7OYrscYXBG3/hx\nJbR0uDouD9vu5i6qK4o55dgpfOKUw1m3eQ+33P8yzW0B/CVeunsOppw3tSnLT0SGlk6AOtlae5Ix\n5pPW2i5jzBeA17PdMCcYvC1HT1+Im+9/JdfNcqyWjl5efuMDXn7jgwHH44NTPGX5iUgq6YyvhKPD\nfLHFQJNwaSWJZGLbctROKKOmsiTXzXENZfmJSCrp9KDuBlYCU40xdxNZtHtbVlvlIIMn9k8ykzUv\nNUqU5SeS/7KZ2ZdOFt9DxpgNwNlEelznW2s3j8qrO1gwFGLFqm2HTOx/6vSjaOsK8MobjbluYt4b\nnOUnIhIvnSy+YuATwMeJlDrqNsa8bq119TDfilXbBvSUYhP76zbvprtH20kMl8cDNdEsv1gSiohI\nIukM8f03UEZkew0vcAlwHJF0c1cK9Aapb0jcQ1JwGr7qihK+vnQetRPK1HMSkSGlE6BOtdbOid0w\nxjwJbMlek3JPVSSy44SZNUyvrch1M0QkT6QToN43xsyy1sb28p4CJN7NyiVURSI7FsyZkusmiMgY\nGJw4MdykiXQCVDHwmjHmT0QW6p4B7DbGrAKw1i4a1is7mKpIZMfUGu2yKyLpSydA3TLo9p3ZaIjT\nDK4iUVzkJdCrDfhGIhgKqx6fiKQtnTTzNWPREKcZXEXi96+8x5r6XbluVt7yAk+9/C5/frtZ9fhE\nJC3p9KAKmr/Yx/gKP1u2N+W6KXktBPxp0+7+26rHJyJD0U/XNCirL3vqG/YR6FXqvogcKp2FukcD\n5wGzifwQ3gY8aa19N8ttc4xMsvomjCtif2ffGLTKHWL1+CZXK4FCxK3is/oyyehL2oMyxkw1xqwA\nfgXMIBKYLHAk8GtjzApjzPRhtjevxLL60tHapeCUCdXjE5FkUvWgvgvcZq19I9FJY8xc4A4KYF8o\nGJjV19TWnfR+YVcXgBp9qscnIskkDVDW2i+keqC19jUKJDjBwKy+5rZuVm7YwWtbG2lu1waGw1Fd\n4efkObWqxyciSaUzB3UMcBWRfaA8sePW2suz2C7H8hf7mDpxHBd/wjD3mBrufqQg9m4cVR4P3HTp\nyUyoKM11U0TEwdLJ4ns0er+1wJq4/wpaMBRi47Z9uW5GXgqH4fG1b+e6GSLicOmsg/JYa6/Lekvy\nzIpV2was65HMrH1tDz6fj2WLZ2uhrkgBiWX0pZPNl843wwvGmAuNMfoWiUq1HYekJwys3riTFau2\nEegNsrel65D1UIHeILv3dWqdlEiBStqDMsaEiHyPeIAvAWFjDNHbYWttwaZeaeHu6Fm3efchuxZ/\n5qxjeOS5tyPH2wPUVKoskkghSpXFl2qNVEEvXNF2HKOnuyfYvwlkrPyRfW8/7+/t6L+PyiKJFKYh\nf44aY14cdNsLrM9ai/JAJgt3JXM7GzsSHldZJJHCkmqIbxVwVvTP8ftM9AFPZLdZzrd00SyCoTBr\n6ncS0uLcUZXs/VRZJBF3SLfcUaohvkUAxpgfWmuvHaV2uYbP6+XiTxgIh1mtbThGldeTOEipLJJI\nYUmZZm6MOZfohoXGmEuADwPrrbX/k/2m5Ydl59Th83n7NzasriylvLRowByKZObw2oqE7198WSRt\nfCjifqkcgjNWAAAgAElEQVSG+O4G5gNLjTH/QiQ4PQ582hhzkrX2H8aojY42eGPD8RV+inweVqza\nxvq39rK/Q6WQhuL1QjgENVWlzK+bFJfFdzDoz6+bFB1WDbFi1bZDMv+U4SfiPql6UJ8ATrDWBo0x\n5wGnWWsDxph7gS1j07z84S/2DZgbWba4jvM/OoNbH3iVlg5l+6USCsFpH5rCFz45p783FAv6vpJi\ngj29/ceXr2zoz+gDZfiJuFmqn5xdwOTonz8AxkX/PI5IooQMobK8hJPnKNsvHRvtXoKh0IBj/mIf\nUyeNGzCsl2yBtDL8RNwnVYC6DXjVGHMn8BawxhjzA+Al4K6xaJwbLF00i8ULplOtyf2UeoJhbn9w\nwyFBKl6qBdKxDD8Rcb74DQxTSRqgrLVPAmcCu4AS4EWgHbhUSRLpi81R3Xr5KUyoKMl1cxxtT3MX\nt/3s1aRBqsxfxIQkgV4ZfiLukzKLz1r7F+D7Y9QWV6ssL+EkU8uqDen9cihUOxo7Wf5MAxctmk1r\nR4DK8WUDEiOSzedp40MR90mnmrmMEs/QdxHg+df3sHl7E81tAWqry/AX+5Km7U+sOpjhJyLukirN\nfBcwJcGpgi8WOxyB3iCbtmr/qHT09IX66xzubTmQ9H4TKkq4+dIFVJZr6FTEjVL1oE4BVgEXWmvf\nGM6TG2NOBf7NWnvWoONfB74IxFKyrrLW2uG8Rr4YqgL6sUdOYE/LAVraNdGfrrbOHg4E+hSgRFwq\nVamjndFAchvw2Uyf2BhzPXAx0Jng9MnAJdbaDZk+b74aqgL6m+/tx6MxwIwoMUIkP424Fh+AtfYp\n4KlhtmE78GngoQTnTga+ZYw5DPg/a+0dQz1ZdXU5RUWZjyrW1lZm/JhsmVs3mVXr3096Pqyisxk5\nfe40pk+bkOtmjAon/T3NBrdfH4zNNY4rL8Hrgoop6b5XGSVJGGOutNbem859rbWPGmNmJDn9MHAP\n0AY8Zow5z1r7u1TP19LSlUlTgcib0NjYnvHjsuXTZx7NC5t30t2TfK2PJHbE5Aq6uvsGlD46/yNH\nOurzHS6n/T0dbW6/PhjZNWYS2Dq73FE6bfB7lew9yDSL70tAWgEqGWOMB7jbWtsavf1/RGr+pQxQ\nblDuL+KME6cNKNUjqU2uLuPEmRNZumgWfcGwCsSKFJBMA9RozJJUAVuMMccSmZ9aBDwwCs+bF2Lp\n0PUNjdqRNw13fOUMPMFICaO+oEoZiRSSTAPUfw73hYwxy4AKa+29xpgbgNVAAHg2OtdVEOKrnz/0\ntOWFLXty3SRH293UyZSqElUwF3GJdBMkADzhIWbmjTHFwGJgEnE9KGvtz4fZvmFpbGzPOIXA6WPf\nsQoJ6zbvprtHvYNEzj31SPr6gjyboALH4gXTXVHB3Ol/T0fK7dcHI56DSntk6jfPvJX3qVSJAlSy\n9yCdHtRvgKnAm0DszQkDYxqg3CjWm7rgzGP41TMNvPVeCy3tAaor/ZQU+djdnHliiNs8/fJ7lJYk\nnm/aaBv52Nxp1E4oSzgnpU0NRfJbOgFqjrV2TtZbUsDK/UX8/XkfGvCFCnDjfS9pngqS9i6b2wPc\ncv8rhwz5aVNDEXdI51/rdmPMkVlvifRveugv9uEv9jG/TntJDSXMwU0LV6zaBsCKVdtYuX4HTW2B\nhOdFJD+kqsW3msi//8nA68aY14jbqNBauyj7zStsBzP+9tHc1s34ihLKSovYva+whv5KS3xpzdHV\nN+zj/I/OSLmp4ZKFMzXcJ5InUg3x3TpWjZDE4jP+YkN/rR0BvvVfL5H3M6UZOO34KRR5vf2BOtm1\nt7R3s2Nvx5CbGk6uLs9eY0UkqUwy+CB1Lb41AMaYH1lrr4k/Z4x5EFgznAZK5mJDfzB0TT83am7t\n5lOnz+D8j86gtSPADx/ZnPD6qytLmT65Iun7o9p9Ivkl1RDffwPHAAuMMccNeow7CqDlodjcVCFV\no9i8vZnN25vxeeHMedOYN3tSwrTz+XWTqCwvSfr+aFNDkfySaojvdmAG8EMiFc1j+oiknEuOLF00\niwPdfTxfYIt8gyF4buMuFp18OIsXTKe+Yd+A2nyxObv4ubtE50UkP6Qa4nvHGPMecDVwOBACdgHr\nrbWFM77kQD6vl8+fa3j1rT309A19f7fZ8OYe/uEz8zjtQ5MpKS46ZB1Uork79ZxE8k+qIb6PAg8C\n7wCxn+pTgdnGmMuttc9mv3mSitfrAwqvAkVrV5B/+XlkK7GSIi+nfqiWvzp1BjVVpQMCUfzcnYjk\nn1RDfPcCf22t3Rp/0BgzC3gMOCGbDZPUWjsCBFQeiZ6+EGs3f8DazR8wUQtyRRwr0ww+SL1Q1zc4\nOEW9zehUNZcRiGXzyUFakCviLql6UL8zxjxJZHPB3dFjhwF/x/B32ZVRUojZfOnSglwRd0iVJHGd\nMWYJ8DfANCK9pp3Az6y1j4xR+ySFpYtm8fauNt7e1ZbrpjiKFuSKuEPKYrHW2keBR6NbbhwL9Fhr\n3xqTlsmQ+oJhOg705roZjqMFuSLukHQOyhjz++j/jyey7ulBYIUx5vXoMcmx1o4AjfsP5LoZjnPi\nzJpIEkmvkkhEcu2seYcPK0ECUvegpkT/fzdwjbU2FrAWAvcDpw7rFWXUjK/wUzuhjL0tClIQGYM+\nvHYcm7c38Vz9Lm2zIZLn0vlXOz4WnKC/Rp8G9x3AX+zjtOOn5roZjhEGdjR2apsNEZdIFaBmGWN+\nCgSMMVcAGGOqjTHf5GBWn+TY5ecfx+IF05lYVYrHA14tADhEfcM+DfeJ5KFUQ3zHAacA+4lUkIBI\nivmpwKXZbZaky+cbWNbn6VfeY3X9rlw3y1GU1SeSn1IFqEZr7W+B38YOWGt/DPw4dtsYU2qt7c5i\n+yRNsbI+y86pw+fzstE20tyukokAEyr89PSFCPQGtTZKJI+kClC/NMb8AXjYWtsef8IYUwlcAiwG\nLsxi+yRD8YVSH/rDW7zw5w9y3aSc6wr0ccv9ryhpQmSMDTd7LyZVgPos8GXgVWPMfmAHka02ZgAT\niWzD8dkRvbpkjb/Yx2V/cyzlZcW8+uYHtHYW7nqp2HbxsaQJgGWL63LZJBFJQ6pKEiHgHuAeY8xc\nYDaRLTe2W2tfG6P2yQjEelPnnnIE1//0xYLaJj4VlUISyQ+pttt4HtgK/B74o4JS/gqGwgpOcZQ0\nIZIfUvWgTo9urfFJ4OfGmDLgOeD31toNY9Q+GQWxyufNbUqaAJVCEskXKWeKrbXbrLU/staeD5wH\nrAcuMcasG5PWyajwF/s49sjqXDfDMU6YWa3hPZE8kLJYLIAx5mgiwSk2B7UV+H6W2yWj7HPn1LGh\nYS/dPaFcNyXnXtzyAUU+n7L5RLJkpNl7MamKxU41xqwAfkUkc28bYIEjgV8bY1YYY6aPSisk68r9\nRZxx4rRcN8MRAr0hVq7fwfKVifbjFBGnSNWD+i5wm7X2jUQno5l9dwAXZ6NhMvqWLppFOBxm7ebd\n9PSqJ7WmficHuvv4/LmGcv+QgwkiMsZSJUl8IdUDo1l9Ck55xOf18nfnGD51+tHc8sAr7O/oyXWT\ncioUhpfe+IBN2xo548RpGvITcZh05qCOAa4CJhHZ0QAAa+3lWWyXZFFleQkL5kzWdvFR3T0hLeAV\ncaB0xjUeBVYCa0HLadxi6aJZHOju4/kte3LdFMfYaBtZsnAmENkMcnyFX9l+IhkYreSImHQClMda\ne92ovqrknM/r5fPnGt54p4mWjsItgxSvuT3At3/2Kj19QZrbAqrdV8ACvUH9SHGAdALUC8aYC4H/\njZY/EpfwF/s4dsZEXlAvqt/u5q7+Pzupdp++MMdGMBRixapt1Dc06keKA6QqdRQiMqTnAb4EhI0x\nRG+HrbX6V+ICy86ZzQa7l4Cy+pLKZe0+fWGOrRWrtg2Ym3XSj5RClCqLL9UaKdWJcYlyfzEnm8nq\nRaWQy9p9+sIcO4HeIPUNjQnPqcBwbgz5E8wY8+Kg214iJY/EJZadM5vSEv3DSybd2n2B3iB7W7pG\nbXv5ob4wtY396GrtCCStVxn7kSJjK9UQ3yrgrOif48d/+oAnstssGUvl/mLOOHGq0s6TmF83acAv\n58HzQekOw2U6j5TOF6Yqso+eWFHlpgTvuQoMD220M/gg9RDfIgBjzA+ttdeO+iuLo1xw5jGs27xL\ntfoGKS3xEgqHCYYi70uiQBQKh1m1YWf/Y2LDcMFgiHM/fCQV5cU8vvYvSQNYd08fe1u6Dglc2fjC\nzEWyRbLrcxp/sY/5dbUJf6gN/pEiYyNlFp8xZg6wzxjzn0QKxe5C2224UkdXDwEFp0N094RYtWEn\n4TD09AQHrBuLBaLSksQj5Ws27eK5+l34S3z9u/rGP66zu5fSYh9/fqeFxpYDhwSu0fzCzEWyRew1\nN29vSnh9TrR00SwgMoTa0t5NdWUp8+sm9R+XsZVqiO8rwJXAb4BXo4enAv9tjPmFtfauMWifjJFU\nv9YFVm/cmfRcsl5nKBw7n3iu6MUtHwy4nSgBItEX5omzJnL2/MMJ9AbTDlLpJFuMdu8qHxM8YrtQ\nL1k4U2n9DpCqB/U1YJ61tiv+oDHm+8BGQAHKRVL9WpextW7zbi4482jK/cUDvjD3NHfx9MvvsXnb\nPp7buDPtHslQyRYXnHl0yiHI4cj3jDh/sU/zew6QKkD1AsUJjpdFz4nLHPprXT2qXOjuCbL8ma18\n8bwPAZGhskfXbD9kjjDdHslQyRbLn9k6YJnBaPR0lOAhoyFVgPoOUG+MeRbYHT02FVgE/HO2GyZj\nL9Hwxu0PvsrOfV1DP1hG1Z//0kx7Vw+V5SWHDJUNNlSPpKK8+JB5sJjqSj9vvduc8fMONRyojDgZ\nDamy+JYbY54DFgPTiFSQWAvcYq3dNTbNk1yIDW8EeoNJ508ku1o7e7j1gVeZVzeJ17YmHiqLaWrr\nprmtm6kTxyU8//javyT9HOccWZ10kXaink66yRbKiJPRkKpaxOHW2l3W2p8DjwGdQDmgv1kFItUw\njWRfS0eA1Rt30tw+9L5dKzck7mGlmgvyF3s5a940qitLEp4f3NNp7+rhp49vYeX6HTS1BQhzcDhw\nxapthzx+6aJZLF4wncnVZXg9MLGqlMULpnPBmceM6oJmca9UQ3xPAicZYz4DfJ/Ithte4HpjzPXW\n2kfGooGSO8rsyx+vbd3H2fMPZ/y4Eg4E+vqH3lL9yAj0hvjXX2ykpDjx79S5syfiL/bR09fHd36+\nkR17O5Lut5NoODA2ZHzVkjK2v9PUvx7slvtfVl1BSUs61cz/H3CmtfZdAGPMvwKrAQUol1NmX/5o\nbg9w8/2v4PVE0ttrKks49qgalpw1M+WPjDAkLRQc2530Oz/fyPt7O1K/fls3b+9s5ZjDxx8yfFda\nUsTk6nKWr2zIu7Rzya10AlQn8H7c7b1o48KCEZ/Z19zWrQ/e4WJrr5rbe3h+yx42NOyldkL5sHrB\nm7Y28VcfPsDOxtTBCcDjgTsf3pSyzFM+p51LatkocwSpi8WON8a8AUwhMsSHMeZk4H+JJEtIAYgN\n09x+xancdvkpTKxS9lU+6e4J8f7eDqZOKsOX4ShaS3s39r39/UEvlVCYlHNSrR2BpEGyuW10C7GO\ndtFeyZ1UWXwzjTE1wClAafTwYcDTwH+OQdvEQfzFPqZPrtSQX57ave9Axo8pLvImTaAYyuBe0fgK\nP6Ul3oRVN/wlvlFJO9feWe6TcojPWttMJCDFHG6tvSe7TRInW7poFva9/UPOSUj+C/SG+PeHXxvW\nY2O9ooGLcT0J7xsOh2ls6aK2ujzpMF86ZZjysbSSpJbOHFS8LwH3ZqMhkh/6gmG6ulVIRFLzeODp\nV99n2eLZQGSIL5BkLVagN8TND7zKxAQ9nky2MtEcl/tk2u9N/BNICkaquQSRmFA4UmD3/t+9wTu7\nWynzF1EzxPxlovmrWK9oqHVX2mzQnTLtQWnuqcCNr/AzoaKE/R1DLx4VeemNvbz0xl5KS3xMmlAK\nafy4ifV4In9Or1ek0kpjJ1sZe4mks+V7sTHmk8aYi4EDxphLjDGXjEHbxIH8xT7mz56U62ZInunu\nCbJjbydHTK5gYlUpnhRjMbEeT6peUfOgXlFszV4iiUorKdMvP6TTg/oNkSKxb3Jw/VMY+Hm2GiXO\ntuycOrbtbFOihGSsq7uXmy89hdbOHu7+9aaEZZziezzJekUe4OlX3mPZOXX4vF6CoRChcHhApmBp\niY/TTzhswGaDyvTLL+kEqDnW2jlZb4nkDZ/Xy82XLmD5yq3UNzRquE/S1twe4ECgj+m1FZxkJg9Z\nTDbZsoZQGFbX78Lni6zTW7FqG6s2DNxUsrsniMfjGRB4lOmXX9L5ybDdGHNk1lsiecXn9XLxJwx3\nXPURPnr8YblujuSJmkp/f+8oVkx2YlXpgGKy8T2epYtmcfZJh+NNMiRY37CP9q6elHNVsWG8oTL9\nNNznPKm2fF9NZChvMvC6MeY1oC923lq7aKgnN8acCvybtfasQcfPB26OPt8D1tr7htV6yTl/sY/L\n/noO5aVFvPzGB7R3KQVdkovNE+1t6WJ8hZ8lC2fysbnTIBymNrpmqqm1u3+9k8/r5dxTjmD1xp0J\nn6+lvZsdezvS2hwx3U0U01lzJWMj1RDfrSN5YmPM9cDFRGr5xR8vBn5ApEJFJ/C8MeYJa+0HI3k9\nyZ1YOaRzTzmC6//zRcIq2CcJnPahKYTCYW687yWa2gKUlngBD4GeIDVVfspLi+k80ENLew81VX5O\nnDWJxSdPp8xflLQKRXVlKdMnV6SVwTdUpl9FeQnLVzYU7PzUWGbnpStVqaM1AMaYH1lrr4k/Z4x5\nEFgzxHNvBz4NPDTo+LHANmttS/S51gEfI5KMIXls4vgyptdWKHlCDuEv9lJWWjRgnmjw9vXxgaOp\nLbIX1uqNO5MGJ4jMV1WWl6S1OeJQmyg+vvZtzU85TKohvv8GjgEWGGOOG/SYCUM9sbX2UWPMjASn\nqoDWuNvtwPihnq+6upyiosy727W1lRk/Jt846Rrv/vpCrvvRWt7Z00Yo8XeKFCCPB7b8JfHW8kNJ\nFpzK/D6uuOAEystKuPqi+ZSXlfDSlt3s23+ASRPKOO34qVx+/nH0BkO0tAWorvInvd/fnWu45q7n\nEr7O5u1NXLWkjNKS1DllY/HvcFx5Cd4s9eac9D0Sk+odvx2YAfwQuC3ueB+RlPPhagPi34lKYP9Q\nD2pp6cr4hWprK2lsbM/4cfnEidd44yULaO/q4b097axvaGTL2800tXXnulmSQ909IQI9mResTeVA\nIMgPH67nsr+eg8/r5YLTZ/DJDx/RP39U5PPw41/XJxyyi7+fv9jHX95vobElcfv27T/A9neaBtUV\nHGgk/w4zCQydXdnLmM3l90iy9yDVEN87wDvA3NgxY8znrLW/GmFb3gRmRyuldxAZ3rtzhM8pDlNZ\nXsJxx0zkuGMmEugN0tzWzU8f38KOxs6hHyyuU1Ppx+Nh1MtkvbBlD+WlRf1DcP5iX38gGWqDxPiA\no0oUzpRqiC9RtYhvR5McsNZmtFDXGLMMqLDW3muM+UciVdK9RLL4EqfoiCtEvjTKMEdOoHH/gaQ7\nuIp7lfl9eDweYPRr4iUqBptp8dih5qfyMZvPiUkPmUo1xPclYDbwJAeLxFYCZ5NmJYloL+y06J+X\nxx1/Mvq8UiBWrNrGsxv0O6RQ7dyXeoi+tMTLpAlldB3oo7k9syAWnyIek25KebwLzjyGA919vPVe\nCy3tAaorS5lfN2nAuiwZW6kC1JnALcDxwFXW2kZjTL219rKxaZq4RapfsyIQmaMyR0zgU6cfzc0P\nvEJrBtVJEg3BZTJkN7j8UXVlCacddxjLzplNub94+BclI5Y0HcRaG7TW3gzcBTxhjDmPg7X4RNKW\n6tesSMzzr++JFonNLBEg0RBcJsVjB2/p0dzewwtb9vD42r9kfA0yuobMV7TWPg+cC1xEpKqESEZi\nv2bFfXyjmPHc3ROkty/EhIrU28yPH1eStDRSvHRKKan8kbOltR+UtbbNGHMj8DFjzExr7fYst0tc\nJNUEtOS3M+dNo8jrZcNbjbSMwqaAxUVe5s+exOr6XQnPT6wq5eZLF3Ag0DdkKaJYhZMlC2cmLV00\nnLkqGTupsvg+DvwP0AX8E/AfwEvAbcaYb1lrfz0mLRRXiP1qrW9opKktgNcTqUjtL/Yqqy8P+Iu9\n1FaX0dnVy/6OHqor/ZxkImuK+oJhTq6bxPeWbxrRHEBpiY/a6vKU27nEKkdUlqfuZQ1suy9pkMnn\n9HI3ZOkNJVUP6t+BTwIVwGrgRGvtVmNMLfAMoAAlaRv8a7bMX8SBQB8V5SU8vvZt6hv20dzWTbHP\nQ09QU51OUVHq4+sXzWNabQX+Yt+AQqpFPg/Ln2mgfuu+Udly5aMnHNbfw4lt57KpYR/7OwPUZCmj\nzo3p5W6SKkAVWWu3GGO8QKu1ditANJsv063iRYCBv2Zjv4IHB65bH3iZlg5VRc8Wjwfmz5rEO3va\naWkPpOz1nHLsFI6edrASWezzC4ZCfPt/1o9K3cWauN5YTGw7l4vOnpX1yuIHe/f7aGnvVnq5g6QK\nNJuNMcuBcUT2hPp34GfAhcC2sWicFI74wFVaUgwoQGWLB9i4dR8Tq/yc+qEp1Dc0Eug7dJjV54Ul\nZyX+kl6+cuuIg1NNpZ+vXTSX2gllSYNPquG50ZLOXJXkRqocnMuAlcDvgLOAHiLDescCV2a9ZVKQ\nAr1BunsUnLIpFO0yNbUFeOmNDxIGJ4BQCDoS1H4L9AbZ1LBvxO04ydQyPTp06ASxYOiU9kjqHpTP\nWvtA3O1/jv7XzxhTaq1VFVAZNa0dAfZreG9Yaqr8HH90NX96bU9Gj4slrBz6fImTBCKfUfKMvfHj\nSvjaRXP502u72LytiZb2bkqiX/qRvZ80hCbpSRWgfmmM+QPwsLV2QJlbY0wlcAmwmMiQn8ioSJVV\nJal1dfdSXOTjrJOm8dKWD+juiazh8Rd7+dDRNdQn6fUkCk6QPElgqM/opLpJHDWlkos/YQicHcRX\nUkww2ivWENrIFUL2XkyqIb7PAiXAq8aYl4wxjxhjHjbGvARsJDKU/dmxaKQUjlQVACRi/uyJnDV/\nGqUlA7/ku3tCkXqH4TDzZk2kqjzy+3NcaRETxpUwMcli6YlVfs6ePy3lgtZ4qT6jIyZXsOycugH3\nnTppHP5in4bQJGOpttsIAfcA9xhj5hIpHBsCtltrXxuj9kkBSpRVVeb3aasOoKq8mCs/dTwAm7ft\n6+8lxXuufveA283tPayu38URkysS9nrm19WybHHdgBTyoYJI/GfU3N7NhHF+5tVNYtni2QWxPbqM\njXQrSbwGKCjJmIjPqooNDxX5PNGCnvsKevPDY2fU4C/2sbeli5b2zNYedR7o5eyTDu+fFxqcTp1J\nxpwy32QsaD2TOJa/2EftpHH9O30uW1zH+R+dwS0PvDIqC0Pzjc8LF58bGT4bzlzd/o4A555yxCFr\niwK9QZpau4YVZMYiDVwKlwKU5JUDgb6Mq127xcL5h/dv/zCc+oax0j3xi22Xr2xIuCW6humcqZAS\nJCCNauYiTpIPldGnTx7HGScelvT84OSGodRU+lm8YDqf+/jsAcfjq3V7iKSLpzLUNhOxLdFXrNI6\nfHEG9aAkr2Tacygp8tKTZCGqv9hDoDdxjvW40iI6u/uSPq/XA5MnlHHklAq27WqjpS3A+IqSaMJB\nJJC8u6cjYbWFj55wGF6Ppz/BoLrSz+zp4yn1F7Fle3P//NCJsyay+OTp1FSVJhx6GzwP9PSr77N6\n46G7FpeW+DjjxKkZbTMxeEt0kVxQgJK8s3TRLMLhMM+/vmfAWp/aCWUcCPQN2K77gjOPobUjwMr1\n77M57ss//twfX32fzdua2N8Z6M9GW7poJo8893b/8FdJkXdAxYVQGPa0HOD4mRP5zhWnJUwUSFXw\n1Of19geWmTMm0t56ACCjTLqY2JBdJIPOE5f96GfOkdV87pw6yv0D/6lrmwnJB55wOD8qRzc2tmfc\n0Nrayv4Jdrdy+zWmur5Ab5DG/QcgHKY2ur4m1Rd8Judi24BvtHtpbu/BQ+LtpCdWlXL7FaemDCZD\nBZ3R/gzTCXKB3iA33vdSwiSLdK4pE27/Owoju8ba2sohBmcPGs73YD5I9h5oDkrylr/Yx/TaCqZP\nruz/Mk21GDSTc7H5meZoKneyb4VYb2Oodo7lAtV0Xi+TLdFFckVDfCKDpJqfGWyoTe2GM2Q3VrTN\nhDidApTIIKnmZwZL1tuIDRE6OYVbi23F6RSgRAZJtQjW64kM9w21w2tsiDAmlsINkQXHTqLFtuJU\nClAig6RKZV84bxrnfvjIIRMQlMItMnIKUCIJpJqfGWqITincIqNDAUokgZHMz6QaIhwqqUJEDnLG\nbK2IQw0nRVwp3CKjQz0okSxQCrfIyClAiWSBUrhFRk4BSiSLlMItMnyagxIREUdSgBIREUdSgBIR\nEUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdS\ngBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIR\nEUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUcqytYTG2O8wE+AuUAA\n+KK1dlvc+a8DXwQao4eustbabLVHRETyS9YCFHABUGqt/Ygx5jTgLuBv486fDFxird2QxTaIiEie\nyuYQ3xnAHwCstS8BCwadPxn4ljFmnTHmW1lsh4iI5KFs9qCqgNa420FjTJG1ti96+2HgHqANeMwY\nc5619nfJnqy6upyiIl/Gjaitrcz4MfnG7dfo9usD91+j268PxuYah/s9mK+yGaDagPhPzBsLTsYY\nDy8jnJAAAAhhSURBVHC3tbY1evv/gPlA0gDV0tKVcQNqaytpbGzP+HH5xO3X6PbrA/dfo9uvL9Ab\nxFdSTLCnF39xdn9ED+d7MB8kew+yGaCeB84Hfh2dg3o97lwVsMUYcyzQCSwCHshiW0RERlUwFGLF\nqm3UNzTS3B6gptLP/Lpali6ahc+rBOnRkM0A9RhwjjHmBcADXGaMWQZUWGvvNcbcAKwmkuH3rLX2\nqSy2RURkVK1YtY2V63f0325qC/TfXra4LlfNcpWsBShrbQj40qDDb8Wdfwh4KFuvLyKSLYHeIPUN\njQnP1Tc0smThzGEN98lA6oeKiGSotSNAU1sg4bmmtgCtHYnPSWYUoEREMlTmL8LrSXzO64mcl5FT\ngBIRydCBQB+hcOJzoXDkvIycApSISIbK/EUk6UDhQT2o0aIAJSKSoQOBPpJ0oAijHtRoUYASEcnQ\n+Ao/E6v8Cc9NrPIzviLxOcmMApSISIb8xT7m19UmPDe/rlYp5qNEA6UiIsOwdNEsAOob9tHS3k11\nZSnz6yb1H5eRU4ASERkGn9fLssV1LFk4c0S1+CQ5DfGJiIyAv9jH1EnjFJyyQAFKREQcSQFKREQc\nSQFKREQcSQFKREQcSQFKREQcSQFKREQcSQFKREQcSQFKREQcSQFKREQcSQFKREQcyRMOJ9vVRERE\nJHfUgxIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUdSgBIREUcqynUDRsoY4wV+AswF\nAsAXrbXb4s5/Hfgi0Bg9dJW11o55Q0eBMeZU4N+stWcNOn4+cDPQBzxgrb0vB80bsRTXl/efoTGm\nGHgAmAH4gduttU/Enc/7zzCNa8zrz9EY4wPuAwwQBr5krd0Sdz7vP0OnyfsABVwAlFprP2KMOQ24\nC/jbuPMnA5dYazfkpHWjxBhzPXAx0DnoeDHwA+CU6LnnjTFPWGs/GPtWDl+y64tyw2f4eaDJWnux\nMaYG2AQ8Ae75DElxjVH5/jmeD2CtPd0YcxbwHaLfNS76DB3FDUN8ZwB/ALDWvgQsGHT+ZOBbxph1\nxphvjXXjRtF24NMJjh8LbLPWtlhre4B1wMfGtGWjI9n1gTs+w98AN0X/7CHyKzvGLZ9hqmuEPP8c\nrbWPA1dGbx4F7I877ZbP0FHcEKCqgNa420FjTHzP8GHgS8Ai4AxjzHlj2bjRYq19FOhNcGrw9bcD\n48ekUaMoxfWBCz5Da22HtbbdGFMJPALcGHfaLZ9hqmsEd3yOfcaYB4EfAb+MO+WKz9Bp3BCg2oDK\nuNtea20fgDHGA9xtrd0X/VXzf8D8HLQxmwZffyUDf9nlNTd9hsaYI4DVwEPW2uVxp1zzGSa7Rjd9\njtbaLwB1wH3GmHHRw675DJ3EDXNQzxMZG/51dA7q9bhzVcAWY8yxRMaFFxGZxHWTN4HZ0TH/DiLD\nCnfmtkmjyhWfoTFmCvBH4Gpr7bODTrviMxziGvP+czTGXAxMt9beAXQBoeh/4JLP0GncEKAeA84x\nxrxAZNz7MmPMMqDCWnuvMeYGIr/oAsCz1tqnctjWUTPoGv8ReJpIj/gBa+3O3LZu5Fz4Gd4AVAM3\nGWNi8zT3AeNc9BkOdY35/jn+FviZMeZPQDHwNeBCY4xr/x3mmrbbEBERR3LDHJSIiLiQApSIiDiS\nApSIiDiSApSIiDiSApSIiDiSG9LMxWWMMTOABuCNQafOt9a+n+D+lwLfB96LHioD1gBfiVu0PRHY\nCfyztfauQY//GyIp0hWAj8jShVustSFjzATgHuDE6N13AtdYa7emcR1VwAvAedbad6LHPkKkZlsl\nsBn4grW2xxgzj0hK9njgT0QKkQ4uFSRSUNSDEqfaZa2dN+i/Q4JTnCdi9wM+RKQ22lfjzi8DngSu\njFY1AMAY81fAj4HLrLVziRT7nAvcFr3LHcAWa+0J1toTgAeBFUM1PlqZfR2RigOxY1VE1tJcaa09\nLnr476P//wWRwFdHZD3fFUO9hojbKUCJ61hrg0R6LnVxhy8j0hPqAc6OO/7PwG3W2oboYw8AXyHS\nAwM4DCiNbusCkeB0SxrNuIJIgNwVd+wc4EVr7ebo7WuAx4wxRwFl0WLHAP8DfDaN1xBxNQ3xiVNN\nM8Zsirv9S2vtv6fzwOhw3ieB70ZvzwWmAmuJBJgvAauid58PvBz/eGvtDmBH9ObtwOPAV4wxq4Bn\niPR2UrLWfjH62vGHZwEdxpjHgJnR9nwj2obdcffbDUxP51pF3EwBSpxqV3S4Ll2figY0D5GRgd8C\nv4qeuwz4tbU2aIxZAdxsjJkS3asn9P/bu1/QquIwjOPfYBA0rc20BX2KzeIYNmEwtIuD+QfkFsvA\nYFIwWgYrJrEJGjUotjWZ3mCZexAZFhEEs2E4w/sbzMst9xo8nj2f+OOec3/tPb/zvJy3XTOW7aGk\neWARuEgVlIGkhSkyomPAEnCeysseA3ep79eN+jVmLeJISYGKvnhh+/roYhsktwLsSToYZLkP3KTy\npffUDLHtQ9ecoUZFXKOmNa/Z3gQ2JT0APlGnnncT7vEb8Nb2bvuf58Bt4An1KvHALH++Gow4kpJB\nRd9dBr7bnrU9Z3sOGAC3WrPEQ+C+pNMAkk7SOgJt71MNF3cOZVCnqAe7z1Ps5Q1wro2kALgEDG1/\nAX5KWmzrq8CrKe4f0SspUNF3N4BHI2tPgePAku3XVKPEM0kfgC1gCNxrv70CnAV2JW1TQ/eu2v4x\n6UZaF+IAeClpB5ihTnFQp7x1SR+BE8DGpPeP6Jt8zTwiIjopGVT8FyStUZnQqK+2l//Bfi5QY7/H\nWbadDCniL+UEFRERnZQMKiIiOikFKiIiOikFKiIiOikFKiIiOikFKiIiOuk3UlfbVV0hvT4AAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x116622908>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "skew=(cat['FErr_PACS_160_u']-cat['F_PACS_160'])/(cat['F_PACS_160']-cat['FErr_PACS_160_l'])\n",
    "skew.name='(84th-50th)/(50th-16th) percentile'\n",
    "g=sns.jointplot(x=np.log10(cat['F_PACS_160']),y=skew, kind='scatter')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For 160 $\\mathrm{\\mu m}$ depth is ~ 30mJy"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Add flag to catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from astropy.table import Column\n",
    "cat.add_column(Column(np.zeros(len(cat), dtype=bool),name='flag_PACS_100'))\n",
    "cat.add_column(Column(np.zeros(len(cat), dtype=bool),name='flag_PACS_160'))\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "ind_100=(cat['Pval_res_100']>0.5) | (cat['F_PACS_100'] < 30)\n",
    "ind_160=(cat['Pval_res_160']>0.5) | (cat['F_PACS_160'] < 30) \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "205845 203736 240489\n"
     ]
    }
   ],
   "source": [
    "print(ind_100.sum(),ind_160.sum(),len(cat))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cat['flag_PACS_100'][ind_100]=True\n",
    "cat['flag_PACS_160'][ind_160]=True\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "cat.write('../data/CDFS-SWIRE/PACS/dmu26_XID+PACS_CDFS-SWIRE_cat_20171019.fits', format='fits',overwrite=True)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "## Check maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from astropy.io import fits"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "pacs100_pval=fits.open('../data/CDFS-SWIRE/PACS/dmu26_XID+PACS_PACS_100_CDFS-SWIRE_Bayes_Pval_20171005.fits')\n",
    "pacs160_pval=fits.open('../data/CDFS-SWIRE/PACS/dmu26_XID+PACS_PACS_160_CDFS-SWIRE_Bayes_Pval_20171005.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD3CAYAAADyvkg2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEVFJREFUeJzt3X2MXNV5x/HvrhcMRGvLUtahVSEJVHpCpZJAopgGGyNE\n4gBJ6H9FVloStwiQK0KSVuXFiFJRYUUJESi8pLQuIKiihAglpQIjxSq1XYLbKI1AJU8EaovSKOrg\nGtsRsV3M9I97OVlbuzuzs8POndnvR4riOXvmznl28P3NOWfu9Vi73UaSJIDxQQ9AktQchoIkqTAU\nJEmFoSBJKgwFSVIxMegBLFSrdbDnr0+tWnUK+/a93s/hDIy1NJO1NJO1wNTU5NhM7Ut6pjAxsWzQ\nQ+gba2kma2kma5ndkg4FSdKxDAVJUmEoSJIKQ0GSVBgKkqTCUJAkFYaCJKkwFCRJhaEgSSoMBelt\n9MkvfmfQQ5DmxVCQ+mzT1h1s2rpj1sdSkxkKkqTCUJAkFYaCtEhcQtIwMBSkPjp+L0EaNoaCJKkw\nFCRJhaEgSSoMBUlSYShIkgpDQZJUGAqSpMJQkBaR1y6o6QwFqU884WsUGAqSpMJQkCQVhoIkqTAU\npEXm3oOazFCQJBWGgtQHfvrXqDAUJEmFoSBJKgwFSVJhKEiSCkNBklQYCpKkwlCQJBUTnTpExAnA\nQ8B7gKPAVcAbwINAG3gB2JyZb0bEVcDV9c9vz8wnIuJk4BFgNXAQuDIzWxFxHnBX3ffpzLytfr1b\ngcvq9uszc0//ypWa4a3rGrbdcNGARyIdq5uZwqXARGZ+BPgL4C+BO4EtmbkOGAMuj4hTgeuA84EN\nwB0RsRy4Fni+7vswsKU+7v3ARmAtsCYizomIc4H1wBrgCuCe/pQpSepGN6HwE2AiIsaBFcD/AR8E\nnql//iRwMfBhYHdmHs7M/cBLwNlUJ/2npveNiBXA8sx8OTPbwPb6GGupZg3tzHylft2pfhQqSeqs\n4/IR8AuqpaMfA+8EPgFcUJ/MoVoSWkkVGPunPW+m9ultB47rewZwCNg7wzFasw1u1apTmJhY1kUZ\nM5uamuz5uU1jLYvvk1/8zoKePyx1vmXYxjsXa5lZN6HweWB7Zt4YEacBO4ATp/18EniN6iQ/2aG9\nU98js7TPat++17soYWZTU5O0Wgd7fn6TWMtwGqY6R+l9sZbZg6Sb5aN9/OqT/v8CJwA/jIgL67ZL\ngJ3AHmBdRJwUESuBs6g2oXdT7UuUvpl5ADgSEWdGxBjVHsTOuu+GiBiPiNOB8cx8dV6VSpJ61s1M\n4avAtojYSTVDuAn4V+CBiDgReBF4LDOPRsTdVCf3ceDmzDwUEfcBD0XELqqZwMb6uNcAjwLLqPYR\nngOoX+fZ+hib+1SnJKkLY+12u3OvBmu1DvZcgFPIZhqmWvpxy+xh+VrqML0vnVgLTE1Njs3U7sVr\nkqTCUJAkFYaC1CP/tTWNIkNBklQYCpKkwlCQJBWGgiSpMBQkSYWhIEkqDAVJUmEoSJIKQ0EaMC+C\nU5MYCpKkwlCQJBWGgtQDl3w0qgwFSVJhKEiSCkNBklQYCpKkwlCQGsCNazWFoSBJKgwFSVJhKEjz\n5FKPRpmhIEkqDAVJUmEoSJIKQ0GSVBgKkqTCUJAkFYaCJKkwFCRJhaEgSSoMBUlSYShIkgpDQZJU\nGApSQ2zausOb7WngJrrpFBE3Ap8CTgTuBZ4BHgTawAvA5sx8MyKuAq4G3gBuz8wnIuJk4BFgNXAQ\nuDIzWxFxHnBX3ffpzLytfq1bgcvq9uszc0+/ipUkza3jTCEiLgQ+ApwPrAdOA+4EtmTmOmAMuDwi\nTgWuq/ttAO6IiOXAtcDzdd+HgS31oe8HNgJrgTURcU5EnFu/xhrgCuCePtUpSepCNzOFDcDzwOPA\nCuBPgauoZgsATwIfA44CuzPzMHA4Il4CzqY66X9pWt9bImIFsDwzXwaIiO3AxcBhqllDG3glIiYi\nYiozWwsvVVoYl3a0FHQTCu8E3g18Angv8F1gvD5xQ7UktJIqMPZPe95M7dPbDhzX9wzgELB3hmPM\nGgqrVp3CxMSyLsqY2dTUZM/PbRprGQ1Nrr3JY5sva5lZN6GwF/hxZh4BMiIOUS0hvWUSeI3qJD/Z\nob1T3yOztM9q377XuyhhZlNTk7RaB3t+fpNYy+hoau2j9L5Yy+xB0s23j3YBH4+IsYj4deAdwPfq\nvQaAS4CdwB5gXUScFBErgbOoNqF3A5dO75uZB4AjEXFmRIxRLVHtrPtuiIjxiDidakby6ryrlST1\npONMof4G0QVUJ/1xYDPwH8ADEXEi8CLwWGYejYi7qU7u48DNmXkoIu4DHoqIXVQzgY31oa8BHgWW\nUe0jPAcQETuBZ6e9liRpkYy12+3OvRqs1TrYcwFOIZupqbUs5kbzthsuWrTX6lZT35deWAtMTU2O\nzdTuxWuSpMJQkCQVhoIkqTAUpC544ZqWCkNBklQYCpKkwlCQJBWGgiSpMBQkSYWhIEkqDAVJUmEo\nSA3kdREaFENBklQYCpKkwlCQOnApR0uJoSBJKgwFSVJhKEiSCkNBklQYCpKkwlCQJBWGgiSpMBSk\nhtq0dYfXSGjRGQqSpMJQkCQVhoIkqZgY9ACkpnI9X0uRMwVJUmEoSJIKQ0GSVBgKkqTCUJAkFYaC\nJKkwFKSG86uxWkyGgiSpMBQkSUVXVzRHxGrgB8BHgTeAB4E28AKwOTPfjIirgKvrn9+emU9ExMnA\nI8Bq4CBwZWa2IuI84K6679OZeVv9OrcCl9Xt12fmnr5VKs2DSzZaqjrOFCLiBODrwC/rpjuBLZm5\nDhgDLo+IU4HrgPOBDcAdEbEcuBZ4vu77MLClPsb9wEZgLbAmIs6JiHOB9cAa4Argnv6UKEnqVjfL\nR1+mOon/rH78QeCZ+s9PAhcDHwZ2Z+bhzNwPvAScTXXSf2p634hYASzPzJczsw1sr4+xlmrW0M7M\nV4CJiJhacIWSpK7NuXwUEZ8BWpm5PSJurJvH6pM5VEtCK4EVwP5pT52pfXrbgeP6ngEcAvbOcIzW\nXGNcteoUJiaWzdVlTlNTkz0/t2msZXQ15ffRlHH0g7XMrNOewiagHREXAx+gWgJaPe3nk8BrVCf5\nyQ7tnfoemaV9Tvv2vd6py6ympiZptQ72/PwmsZbR1oTfxyi9L9Yye5DMuXyUmRdk5vrMvBD4N+AP\ngCcj4sK6yyXATmAPsC4iToqIlcBZVJvQu4FLp/fNzAPAkYg4MyLGqPYgdtZ9N0TEeEScDoxn5qvz\nrlSS1LNe/j2FLwIPRMSJwIvAY5l5NCLupjq5jwM3Z+ahiLgPeCgidlHNBDbWx7gGeBRYRrWP8BxA\nROwEnq2PsXkBdUmSejDWbrc792qwVutgzwU4hWymJtTStK+kbrvhokEPoRHvS79YC0xNTY7N1O7F\na5KkwlCQhkDTZi4aXYaCJKkwFCRJhaEgSSoMBek4rt9rKTMUJEmFoSBJKgwFaUi4rKXFYChIkgpD\nQZJUGAqSpMJQkKZx3V5LnaEgSSoMBUlSYShIkgpDQZJUGAqSpMJQkCQVhoI0RDZt3eHXZvW2MhQk\nSYWhIEkqDAVJUjEx6AFITeA6vVRxpiBJKgwFSVJhKEiSCkNBGkLugejtYihIkgpDQZJUGApa8lyK\nkX7FUJAkFYaCJKkwFCRJhaEgDSn3QvR2MBQkScWcN8SLiBOAbcB7gOXA7cC/Aw8CbeAFYHNmvhkR\nVwFXA28At2fmExFxMvAIsBo4CFyZma2IOA+4q+77dGbeVr/ercBldfv1mbmnv+VKkubSaabwaWBv\nZq4DPg58DbgT2FK3jQGXR8SpwHXA+cAG4I6IWA5cCzxf930Y2FIf935gI7AWWBMR50TEucB6YA1w\nBXBP/8qUZuYSjHSsTqHwLeCW+s9jVJ/gPwg8U7c9CVwMfBjYnZmHM3M/8BJwNtVJ/6npfSNiBbA8\nM1/OzDawvT7GWqpZQzszXwEmImKqH0VKkroz5/JRZv4CICImgceoPul/uT6ZQ7UktBJYAeyf9tSZ\n2qe3HTiu7xnAIWDvDMdozTXGVatOYWJi2Vxd5jQ1Ndnzc5vGWpaexf49jdL7Yi0z6/iP7ETEacDj\nwL2Z+XcR8aVpP54EXqM6yU92aO/U98gs7XPat+/1Tl1mNTU1Sat1sOfnN4m1LE2L+XsapffFWmYP\nkjmXjyLiXcDTwJ9l5ra6+YcRcWH950uAncAeYF1EnBQRK4GzqDahdwOXTu+bmQeAIxFxZkSMUe1B\n7Kz7boiI8Yg4HRjPzFfnXam0hLgnon7rNFO4CVgF3BIRb+0tfA64OyJOBF4EHsvMoxFxN9XJfRy4\nOTMPRcR9wEMRsYtqJrCxPsY1wKPAMqp9hOcAImIn8Gx9jM39KlKS1J2xdrvduVeDtVoHey7AKWQz\nLVYto/Ipe9sNFy3K6/jfWDMtYPlobKZ2L16TJBWGgiSpMBQkSYWhIA25TVt3jMz+iAbPUJAkFR0v\nXpNGjZ+qpdk5U5AkFYaCJKkwFCRJhaEgjQj3StQPhoIkqTAUtKT4aVqam6EgSSoMBUlSYShII8Tl\nMS2UoSBJKgwFLRl+ipY6MxSkEeNdU7UQhoIkqfAuqRp5fmqWuudMQZJUGArSiHKGpF4YCpKkwlDQ\nSPPTsjQ/hoI0wgxFzZehoJHlCVGaP0NBklQYCtKIc8ak+TAUNJI8ER7LW1+oW4aCJKnwNhcaKX4a\nlhbGmYK0hBia6sRQ0MjwhNcdf0+ai6GgkeCJTuoP9xQ09AyE+Zv+O9t2w0UDHImaxpmChpZfs+wP\nf4earnEzhYgYB+4F3g8cBv4oM18a7KjUNJ7I+mvT1h3OGAQ0MBSA3wVOyszfiYjzgK8Alw94TBow\nQ+Dt55KSoJmhsBZ4CiAzvx8RHxrweLRIPPE3RzfvhcExmsba7fagx3CMiPhr4NuZ+WT9+BXgjMx8\nY7Ajk6TR18SN5gPA5LTH4waCJC2OJobCbuBSgHpP4fnBDkeSlo4m7ik8Dnw0Iv4ZGAM+O+DxSNKS\n0bg9BUnS4DRx+UiSNCCGgiSpMBQkSUUTN5oXTUQsA+4EPgQsB/48M58Y7KgWJiLeBzwHvCszDw16\nPL2IiJXAI8AK4ETgC5n57GBHNT+jcruWiDgB2Aa8h+rvyO2Z+d2BDmqBImI18APgo5n540GPp1cR\ncSPwKaq/I/dm5t/047hLfabw+8AJmXk+1a00fnPA41mQiFhBdVuQw4MeywJ9AfheZq4HPgPcM9jh\n9KTcrgW4gep9GUafBvZm5jrg48DXBjyeBalD7uvALwc9loWIiAuBjwDnA+uB0/p17KUeChuA/46I\nfwAeAP5+wOPpWUSMAX8F3AS8PuDhLNRXqf7iQjWbHcYZzzG3a6GajQ6jbwG31H8eA4b9QtIvA/cD\nPxv0QBZoA9U1XI9Tnbf6tsKxZJaPIuIPgc8f19yiOuF8ArgA+Nv6/xttllr+C/hGZv4oIgYwqt7M\nUstnM/NfIuJUqmWk6xd/ZAu2Atg/7fHRiJgYtqvzM/MXABExCTwGbBnsiHoXEZ8BWpm5vV56GWbv\nBN5Nde56L/DdiHhfZi74GoMlfZ1CRHwD+FZmfrt+/PPMPHXAw+pJRLwE/LR+eB6wJzMbH3CziYjf\nBr4B/Mlb98EaJhFxJ/D9zPxm/finmfkbAx5WTyLiNKpPpPdm5rZBj6dXEfFPQLv+3weAnwCfysyf\nD3RgPYiIrVQB95X68Y+o9kj+Z6HHXjIzhVnsorqlxrcj4v3AKwMeT88ys+yHRMR/Ah8b2GAWKCJ+\ni2rZ4vcy80eDHk+PdgOfBL45zLdriYh3AU8Df5yZ3xv0eBZi+oekiPhH4JphDITaLuBz9YePXwPe\nAeztx4GXeig8ANwXEd+nWi+9ZsDjUeUO4CTgrnopbH9mDtu/qTEqt2u5CVgF3BIRb+0tXJKZQ71R\nO+wy84mIuADYQ7U3vDkzj/bj2Et6+UiSdKyl/u0jSdI0hoIkqTAUJEmFoSBJKgwFSVJhKEiSCkNB\nklT8Pz6lLV+FlmcHAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1163a8f60>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(pacs100_pval[1].data[np.isfinite(pacs100_pval[1].data)],bins=np.arange(-6,6,0.05));"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAD3CAYAAADyvkg2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEw9JREFUeJzt3X2MXNV5x/HvrhdsqNaupY5Dq0JSSPuEViUB0pgWHCwE\nOEAT+iYVoaQktBSQK6BJ1RIwIkhUUESoQOElJXGAQpUm0LSpI4KroBDbJbhN0whU+kSmLyiNog6W\nsR05toWZ/nGvTwdrd2c9u3juzH4/UpSZM+feOc8Ovr8998y9O9bpdJAkCWB80AOQJDWHoSBJKgwF\nSVJhKEiSCkNBklRMDHoAc9Vu7+7761PLlx/Ljh175nM4A2MtzWQtzWQt0GpNjk3VvqBnChMTiwY9\nhHljLc1kLc1kLdNb0KEgSXojQ0GSVBgKkqTCUJAkFYaCJKkwFCRJhaEgSSoMBUlSYShIkgpDQXoT\nvf9jfzfoIUiHxVCQ5tnltz/N5bc/PehhSH2Z1Q3xIuJfgF310/8E/hR4COgALwBrM/P1iLgCuBJ4\nDbg1MzdExDHAo8AKYDdwWWa2I+IM4O6678bMvKV+r5uBi+r26zJz67xUKknqqWcoRMQSYCwzV3e1\nfRlYl5lfj4gHgIsj4lngGuDdwBJgc0T8A3A18HxmfiIiLgHWAdcCDwC/CfwH8JWIOBUYA84GVgLH\nA08AvzRfxUqDcHDWsP76cwY8Eqm32cwU3gkcGxEb6/43AKcDz9SvPwmcDxwAtmTmPmBfRGwDTgHO\nAu7o6ntTRCwFFmfmSwAR8RRwLrCPatbQAV6OiImIaGVme7rBLV9+7JzuEthqTfa9bdNYS7McWsMo\n1jTMrGVqswmFPcCdwGeAn6U6sI/VB26oTgktA5YCO7u2m6q9u23XIX1PBPYC26fYx7ShMJd7orda\nk7Tbu/vevkmspRm61xIOXWQe1poOGubP5VDWMn2QzCYUvgtsq0PguxGxnWqmcNAk8CrVQX6yR3uv\nvvunaZckHQGz+fbR5cAnASLip6h+y98YEavr1y8ANgFbgVURsSQilgEnUy1CbwEu7O6bmbuA/RFx\nUkSMAWvqfWwB1kTEeEScAIxn5ivzUKckaRZmEwqfBX48IjYDf00VEtcCt9SLy0cDj2fmD4B7qA7u\nTwM3ZuZe4H7gF+rtfx+4pd7vVcBjVGHy7cx8LjO/VW//LNUi89r5KVMaPL+mqmEw1un0/SeOG2Eu\nf6PZ84rNNMy19DrwD/M3kIb5czmUtfg3miVJs2AoSJIKQ0GSVBgKkqTCUJAkFYaCJKkwFKR54nUI\nGgWGgiSpMBQkSYWhIEkqDAXpCHLdQU1nKEiSCkNBklQYCpKkwlCQ5oFrBRoVhoIkqTAUJEmFoSBJ\nKgwFSVJhKEhHmIvSajJDQZJUGAqSpMJQkCQVhoIkqTAUJEmFoSBJKgwFSVJhKEiSiolBD0AaZl6I\nplHjTEGSVBgKkqTCUJAkFbNaU4iIFcC3gPOA14CHgA7wArA2M1+PiCuAK+vXb83MDRFxDPAosALY\nDVyWme2IOAO4u+67MTNvqd/nZuCiuv26zNw6b5VKDXJwLWL99ecMeCTSG/WcKUTEUcCngR/VTXcB\n6zJzFTAGXBwRxwHXAGcCa4DbImIxcDXwfN33EWBdvY8HgEuBs4CVEXFqRJwGnA2sBC4B7p2fEiVJ\nszWb00d3Uh3Ev18/Px14pn78JHAu8B5gS2buy8ydwDbgFKqD/le7+0bEUmBxZr6UmR3gqXofZ1HN\nGjqZ+TIwERGtOVcoSZq1GU8fRcSHgXZmPhURH6+bx+qDOVSnhJYBS4GdXZtO1d7dtuuQvicCe4Ht\nU+yjPdMYly8/lomJRTN1mVGrNdn3tk1jLcNn2OoctvHOxFqm1mtN4XKgExHnAu+iOgW0ouv1SeBV\nqoP8ZI/2Xn33T9M+ox079vTqMq1Wa5J2e3ff2zeJtQynYapzlD4Xa5k+SGY8fZSZ783MszNzNfCv\nwO8AT0bE6rrLBcAmYCuwKiKWRMQy4GSqRegtwIXdfTNzF7A/Ik6KiDGqNYhNdd81ETEeEScA45n5\nymFXKknqWz9XNH8MeDAijgZeBB7PzAMRcQ/VwX0cuDEz90bE/cDDEbGZaiZwab2Pq4DHgEVU6wjP\nAUTEJuDZeh9r51CX9KbzamaNorFOp9O7V4O127v7LsApZDMNSy3zEQrD9JXUYflcZsNaoNWaHJuq\n3YvXJEmFoSBJKgwFSVJhKEgD5GK1msZQkCQVhoIkqTAUJEmFoSBJKgwFSVJhKEiSCkNBklQYCpKk\nwlCQJBWGgiSpMBSkPnh7Co0qQ0GSVBgKkqTCUJAkFYaCJKkwFKQBc9FaTWIoSJIKQ0GSVBgKkqTC\nUJAkFYaCJKkwFCRJhaEgHSa/QqpRZihIkgpDQZJUGAqSpMJQkBrAdQo1haEgSSoMBUlSMdGrQ0Qs\nAh4EAugAVwF7gYfq5y8AazPz9Yi4ArgSeA24NTM3RMQxwKPACmA3cFlmtiPiDODuuu/GzLylfr+b\ngYvq9usyc+s81itJmsFsZgrvB8jMM4F1wJ8CdwHrMnMVMAZcHBHHAdcAZwJrgNsiYjFwNfB83feR\neh8ADwCXAmcBKyPi1Ig4DTgbWAlcAtw7L1VKkmal50whM/82IjbUT98KvAqcCzxTtz0JnA8cALZk\n5j5gX0RsA06hOujf0dX3pohYCizOzJcAIuKpep/7qGYNHeDliJiIiFZmtqcb3/LlxzIxseiwiu7W\nak32vW3TWMtwG4aah2GMs2UtU+sZCgCZ+VpEPAz8OvBbwHn1gRuqU0LLgKXAzq7Npmrvbtt1SN8T\nqU5LbZ9iH9OGwo4de2ZTwpRarUna7d19b98k1jL8ml7zKH0u1jJ9kMx6oTkzLwN+jmp94Ziulyap\nZg+76scztR9O3+52SdIR0DMUIuJDEfHx+uke4HXgnyNidd12AbAJ2AqsioglEbEMOJlqEXoLcGF3\n38zcBeyPiJMiYoxqDWJT3XdNRIxHxAnAeGa+Mh+FSpJ6m83po78BPhcR3wCOAq4DXgQejIij68eP\nZ+aBiLiH6uA+DtyYmXsj4n7g4YjYDOynWlyG6ltMjwGLqNYRngOIiE3As/U+1s5TnZKkWRjrdDq9\nezVYu7277wI8r9hMTa/lzbr6eP3157wp+50vTf9cDoe1QKs1OTZVuxevSZIKQ0GSVBgKkqRiVtcp\nSPJOploYnClIkgpDQZJUGAqSpMJQkBri8tufdt1CA2coSJIKQ0GSVBgKkqTCUJAkFYaCJKkwFCRJ\nhaEgSSoMBWkWvH5AC4WhIEkqDAWpYZyVaJAMBUlSYShIkgpDQZJUGAqSpMJQkCQVhoIkqTAUJEmF\noSBJKgwFSVJhKEiSCkNBklQYClIP3otIC4mhIEkqDAWpgZydaFAmZnoxIo4C1gNvAxYDtwL/BjwE\ndIAXgLWZ+XpEXAFcCbwG3JqZGyLiGOBRYAWwG7gsM9sRcQZwd913Y2beUr/fzcBFdft1mbl1fsuV\nJM2k10zhg8D2zFwFvA/4FHAXsK5uGwMujojjgGuAM4E1wG0RsRi4Gni+7vsIsK7e7wPApcBZwMqI\nODUiTgPOBlYClwD3zl+ZkqTZ6BUKXwRuqh+PUf0GfzrwTN32JHAu8B5gS2buy8ydwDbgFKqD/le7\n+0bEUmBxZr6UmR3gqXofZ1HNGjqZ+TIwERGt+ShSkjQ7M54+yswfAkTEJPA41W/6d9YHc6hOCS0D\nlgI7uzadqr27bdchfU8E9gLbp9hHe6YxLl9+LBMTi2bqMqNWa7LvbZvGWkZLE38GTRxTv6xlajOG\nAkBEHA98CbgvM/8qIu7oenkSeJXqID/Zo71X3/3TtM9ox449vbpMq9WapN3e3ff2TWIto6dpP4NR\n+lysZfogmfH0UUS8BdgI/Elmrq+bvx0Rq+vHFwCbgK3AqohYEhHLgJOpFqG3ABd2983MXcD+iDgp\nIsao1iA21X3XRMR4RJwAjGfmK4ddqSSpb71mCjcAy4GbIuLg2sK1wD0RcTTwIvB4Zh6IiHuoDu7j\nwI2ZuTci7gcejojNVDOBS+t9XAU8BiyiWkd4DiAiNgHP1vtYO19FSpJmZ6zT6fTu1WDt9u6+C3AK\n2UxNq2VQ1wysv/6cgbzvdJr2ucyFtUCrNTk2VbsXr0mSip4LzdJC5VXFWoicKUiSCkNBaqjLb3/a\n2YqOOENBklQYCpKkwlCQJBWGgiSpMBQkSYWhIEkqDAVJUmEoSJIKQ0GSVBgK0hS8klgLlaEgNZwB\npSPJUJAkFYaCJKkwFCRJhaEgSSoMBUlSYShIkgpDQZJUGAqSpMJQkCQVhoIkqTAUpEM08bYSTRyT\nRpOhIEkqDAVJUmEoSJIKQ0GSVBgKkqTCUJAkFYaC1KXJX/1s8tg0OiZm0ykiVgJ/lpmrI+LtwENA\nB3gBWJuZr0fEFcCVwGvArZm5ISKOAR4FVgC7gcsysx0RZwB31303ZuYt9fvcDFxUt1+XmVvnsVZJ\nUg89ZwoR8cfAZ4AlddNdwLrMXAWMARdHxHHANcCZwBrgtohYDFwNPF/3fQRYV+/jAeBS4CxgZUSc\nGhGnAWcDK4FLgHvnp0RJ0mzNZqbwEvAbwF/Wz08HnqkfPwmcDxwAtmTmPmBfRGwDTqE66N/R1fem\niFgKLM7MlwAi4ingXGAf1ayhA7wcERMR0crM9kyDW778WCYmFs2u2im0WpN9b9s01jL6Bv1zGfT7\nzydrmVrPUMjMJyLibV1NY/WBG6pTQsuApcDOrj5TtXe37Tqk74nAXmD7FPuYMRR27NjTq4RptVqT\ntNu7+96+SaxlYRjkz2WUPhdrmT5IZrWmcIjXux5PAq9SHeQne7T36rt/mnZJ0hHSz7ePvh0Rq+vH\nFwCbgK3AqohYEhHLgJOpFqG3ABd2983MXcD+iDgpIsao1iA21X3XRMR4RJwAjGfmK/0WJkk6fP3M\nFD4GPBgRRwMvAo9n5oGIuIfq4D4O3JiZeyPifuDhiNhMNRO4tN7HVcBjwCKqdYTnACJiE/BsvY+1\nc6hLktSHsU6n07tXg7Xbu/suwPOKzTTIWpp+LcD6688Z2Hv731gzzWFNYWyq9n5mCtLIaXoYHHRw\nnIMMB402r2iWJBWGgiSpMBQkSYWhIEkqDAVJUmEoaMEblm8edRvGMWs4GAqSpMJQkCQVhoIkqTAU\nJEmFoSANKReb9WYwFLSgeWCV3shQkCQVhoIkqTAUJEmFoaAFaxTWE0ahBjWLoSBJKgwFSVJhKEiS\nCv9GsxacUTsP799t1nxypiBJKgwFSVJhKEgjYtROi2kwDAUtKB44pZkZCtIIMfQ0V4aCFgwPmFJv\nhoI0Yi6//WkDUH3zOgWNPA+Q0uw5U5BGlGGofhgKGmkL/cC40OvX4TMUNLI8IEqHr3FrChExDtwH\nvBPYB/xeZm4b7Kg0bAyE/9f9s/D+SOqlcaEA/BqwJDN/OSLOAD4JXDzgMWkIGAS9efM89dLEUDgL\n+CpAZn4zIt494PGogQyAuZnp52dgLGxjnU5n0GN4g4j4DPBEZj5ZP38ZODEzXxvsyCRp9DVxoXkX\nMNn1fNxAkKQjo4mhsAW4EKBeU3h+sMORpIWjiWsKXwLOi4h/BMaAjwx4PJK0YDRuTUGSNDhNPH0k\nSRoQQ0GSVBgKkqSiiQvNR0xELALuAt4NLAY+kZkbBjuquYmIdwDPAW/JzL2DHk8/ImIZ8CiwFDga\n+GhmPjvYUR2eUbldS0QcBawH3kb1b+TWzPzyQAc1RxGxAvgWcF5m/vugx9OviPg48AGqfyP3ZeZn\n52O/C32m8CHgqMw8k+pWGm8f8HjmJCKWUt0WZN+gxzJHHwW+lplnAx8G7h3scPpSbtcCXE/1uQyj\nDwLbM3MV8D7gUwMez5zUIfdp4EeDHstcRMRq4FeAM4GzgePna98LPRTWAP8TEV8BHgT+fsDj6VtE\njAF/AdwA7BnwcObqz6n+4UI1mx3GGc8bbtdCNRsdRl8EbqofjwHDfiHpncADwPcHPZA5WkN1DdeX\nqI5b83aGY8GcPoqI3wX+8JDmNtUB51eB9wKfq/+/0aap5b+Bz2fmdyJiAKPqzzS1fCQz/ykijqM6\njXTdkR/ZnC0FdnY9PxARE8N2dX5m/hAgIiaBx4F1gx1R/yLiw0A7M5+qT70Ms58A3kp17PoZ4MsR\n8Y7MnPM1Bgv6OoWI+Dzwxcx8on7+g8w8bsDD6ktEbAO+Vz89A9iamY0PuOlExC8Cnwf+6OB9sIZJ\nRNwFfDMzv1A//15m/vSAh9WXiDie6jfS+zJz/aDH06+I+AbQqf/3LuC7wAcy8wcDHVgfIuJ2qoD7\nZP38O1RrJP87130vmJnCNDZT3VLjiYh4J/DygMfTt8ws6yER8V/A+QMbzBxFxM9Tnbb47cz8zqDH\n06ctwPuBLwzz7Voi4i3ARuAPMvNrgx7PXHT/khQRXweuGsZAqG0Grq1/+fhJ4MeA7fOx44UeCg8C\n90fEN6nOl1414PGochuwBLi7PhW2MzOH7W9qjMrtWm4AlgM3RcTBtYULMnOoF2qHXWZuiIj3Alup\n1obXZuaB+dj3gj59JEl6o4X+7SNJUhdDQZJUGAqSpMJQkCQVhoIkqTAUJEmFoSBJKv4Prfj0VI6V\nII0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1171ebc50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(pacs160_pval[1].data[np.isfinite(pacs160_pval[1].data)],bins=np.arange(-6,6,0.05));"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Check cat"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: UnitsWarning: 'Degrees' did not parse as fits unit: At col 0, Unit 'Degrees' not supported by the FITS standard.  [astropy.units.core]\n"
     ]
    }
   ],
   "source": [
    "old_cat=Table.read('../data/CDFS-SWIRE/PACS/WP4-CDFS-SWIRE-PACSxID24-v1.fits.gz')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "old_cat['RA'].unit=u.degree\n",
    "old_cat['Dec'].unit=u.degree"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;Table length=87953&gt;\n",
       "<table id=\"table4682590584\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>XID</th><th>RA</th><th>Dec</th><th>F_PACS_100__A4</th><th>F_PACS_100__A5</th><th>F_PACS_100__A6</th><th>F_PACS_100__A7</th><th>F_PACS_100__A8</th><th>F_PACS_100__A10</th><th>F_PACS_100</th><th>Ferr_PACS_100__A4</th><th>Ferr_PACS_100__A5</th><th>Ferr_PACS_100__A6</th><th>Ferr_PACS_100__A7</th><th>Ferr_PACS_100__A8</th><th>Ferr_PACS_100__A10</th><th>Ferr_PACS_100</th><th>F_PACS_100__SKY</th><th>F_PACS_160__A4</th><th>F_PACS_160__A5</th><th>F_PACS_160__A6</th><th>F_PACS_160__A7</th><th>F_PACS_160__A8</th><th>F_PACS_160__A10</th><th>F_PACS_160</th><th>Ferr_PACS_160__A4</th><th>Ferr_PACS_160__A5</th><th>Ferr_PACS_160__A6</th><th>Ferr_PACS_160__A7</th><th>Ferr_PACS_160__A8</th><th>Ferr_PACS_160__A10</th><th>Ferr_PACS_160</th><th>F_PACS_160__SKY</th><th>HELP_ID</th></tr></thead>\n",
       "<thead><tr><th></th><th>deg</th><th>deg</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 / pix</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 / pix</th><th></th></tr></thead>\n",
       "<thead><tr><th>int32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>float32</th><th>str21</th></tr></thead>\n",
       "<tr><td>20269</td><td>51.3391</td><td>-29.6326</td><td>5.80525</td><td>9.10904</td><td>9.52371</td><td>8.90982</td><td>9.87109</td><td>16.1481</td><td>0.960175</td><td>10.9545</td><td>10.2805</td><td>10.4431</td><td>10.9296</td><td>11.7329</td><td>13.4598</td><td>8.54818</td><td>0.0</td><td>17.5657</td><td>16.0669</td><td>17.1881</td><td>18.7219</td><td>23.2457</td><td>30.4306</td><td>21.9574</td><td>18.1397</td><td>15.2603</td><td>14.225</td><td>13.5566</td><td>13.3353</td><td>14.2112</td><td>10.8895</td><td>0.0</td><td>CDS-PACSxID24-1-00001</td></tr>\n",
       "<tr><td>9465</td><td>51.3404</td><td>-29.5366</td><td>37.1872</td><td>36.8413</td><td>32.099</td><td>29.5196</td><td>28.2173</td><td>30.7685</td><td>30.7021</td><td>10.9781</td><td>10.2845</td><td>10.4078</td><td>10.9001</td><td>11.7103</td><td>13.4265</td><td>8.53753</td><td>0.0</td><td>63.0557</td><td>59.6926</td><td>59.6764</td><td>60.0185</td><td>60.3393</td><td>61.5365</td><td>53.9875</td><td>19.0148</td><td>15.9282</td><td>14.8685</td><td>14.2821</td><td>14.1632</td><td>15.0268</td><td>10.9488</td><td>0.0</td><td>CDS-PACSxID24-1-00002</td></tr>\n",
       "<tr><td>32785</td><td>51.3417</td><td>-29.6128</td><td>40.1389</td><td>41.569</td><td>41.6262</td><td>41.6549</td><td>48.1906</td><td>63.5627</td><td>25.1915</td><td>11.258</td><td>10.5272</td><td>10.7061</td><td>11.1753</td><td>11.9642</td><td>13.7515</td><td>8.81871</td><td>0.0</td><td>2.46996</td><td>-3.51323</td><td>-4.20904</td><td>-2.81796</td><td>-0.26346</td><td>1.12468</td><td>16.2734</td><td>18.6486</td><td>15.8986</td><td>14.9312</td><td>14.2132</td><td>13.8468</td><td>14.5681</td><td>11.1448</td><td>0.0</td><td>CDS-PACSxID24-1-00003</td></tr>\n",
       "<tr><td>60323</td><td>51.3436</td><td>-29.6118</td><td>22.0892</td><td>21.0967</td><td>22.4059</td><td>20.2713</td><td>15.4865</td><td>23.8448</td><td>3.80425</td><td>11.086</td><td>10.509</td><td>10.735</td><td>11.1713</td><td>11.9713</td><td>13.6328</td><td>8.65524</td><td>0.0</td><td>6.87234</td><td>1.42008</td><td>1.1628</td><td>2.5351</td><td>0.156342</td><td>-4.5212</td><td>8.33736</td><td>18.4275</td><td>15.682</td><td>14.7175</td><td>14.044</td><td>13.7558</td><td>14.517</td><td>11.1362</td><td>0.0</td><td>CDS-PACSxID24-1-00004</td></tr>\n",
       "<tr><td>22332</td><td>51.3443</td><td>-29.5597</td><td>21.7853</td><td>19.5587</td><td>19.0401</td><td>23.1344</td><td>31.0047</td><td>41.9829</td><td>25.3307</td><td>11.3134</td><td>10.5788</td><td>10.8004</td><td>11.2801</td><td>12.0131</td><td>13.7747</td><td>8.78858</td><td>0.0</td><td>-3.77195</td><td>-2.73229</td><td>-3.90126</td><td>-4.01867</td><td>-1.35882</td><td>-1.09863</td><td>10.3223</td><td>18.2565</td><td>15.6621</td><td>14.7207</td><td>14.0263</td><td>13.6915</td><td>14.5277</td><td>11.1758</td><td>0.0</td><td>CDS-PACSxID24-1-00005</td></tr>\n",
       "<tr><td>18577</td><td>51.3449</td><td>-29.5822</td><td>-9.65322</td><td>-6.59912</td><td>-5.24355</td><td>-3.83109</td><td>-5.384</td><td>-16.1549</td><td>-0.769118</td><td>11.1482</td><td>10.3816</td><td>10.4948</td><td>10.9838</td><td>11.7589</td><td>13.499</td><td>8.6653</td><td>0.0</td><td>7.4133</td><td>8.82751</td><td>12.8263</td><td>15.9227</td><td>21.4362</td><td>25.5777</td><td>42.2368</td><td>18.1161</td><td>15.3303</td><td>14.3075</td><td>13.6264</td><td>13.38</td><td>14.2146</td><td>10.9573</td><td>0.0</td><td>CDS-PACSxID24-1-00006</td></tr>\n",
       "<tr><td>51911</td><td>51.3451</td><td>-29.6181</td><td>14.2468</td><td>11.2671</td><td>10.6499</td><td>6.22995</td><td>2.58583</td><td>-1.99709</td><td>15.2052</td><td>11.2722</td><td>10.4884</td><td>10.7116</td><td>11.2265</td><td>11.9926</td><td>13.7069</td><td>8.66916</td><td>0.0</td><td>5.21845</td><td>7.09556</td><td>10.6759</td><td>10.4271</td><td>9.68967</td><td>13.277</td><td>14.2788</td><td>18.702</td><td>15.8801</td><td>14.8524</td><td>14.1734</td><td>13.862</td><td>14.6394</td><td>11.2637</td><td>0.0</td><td>CDS-PACSxID24-1-00007</td></tr>\n",
       "<tr><td>19642</td><td>51.3451</td><td>-29.6764</td><td>-10.7898</td><td>-12.8643</td><td>-13.062</td><td>-11.6732</td><td>-10.9978</td><td>2.1182</td><td>-8.33675</td><td>11.1092</td><td>10.368</td><td>10.5964</td><td>11.1072</td><td>11.8695</td><td>13.6347</td><td>8.68843</td><td>0.0</td><td>50.9411</td><td>51.5546</td><td>49.9911</td><td>47.5871</td><td>47.1265</td><td>51.2473</td><td>26.5581</td><td>18.2558</td><td>15.642</td><td>14.7067</td><td>14.0133</td><td>13.692</td><td>14.3928</td><td>11.1234</td><td>0.0</td><td>CDS-PACSxID24-1-00008</td></tr>\n",
       "<tr><td>17590</td><td>51.3457</td><td>-29.5624</td><td>11.417</td><td>12.115</td><td>7.56328</td><td>2.81013</td><td>6.40629</td><td>10.1592</td><td>18.007</td><td>11.1196</td><td>10.4087</td><td>10.5542</td><td>10.9821</td><td>11.8081</td><td>13.6093</td><td>8.65609</td><td>0.0</td><td>-13.3776</td><td>-3.94243</td><td>0.20211</td><td>2.54558</td><td>1.6211</td><td>4.05115</td><td>13.4845</td><td>18.7138</td><td>15.9023</td><td>14.9397</td><td>14.243</td><td>13.9395</td><td>14.6481</td><td>11.4432</td><td>0.0</td><td>CDS-PACSxID24-1-00009</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><td>...</td></tr>\n",
       "<tr><td>5205</td><td>54.6421</td><td>-26.8742</td><td>33.5527</td><td>36.4854</td><td>39.1846</td><td>37.6279</td><td>35.524</td><td>26.889</td><td>nan</td><td>9.55403</td><td>8.83414</td><td>8.85859</td><td>9.20788</td><td>9.86909</td><td>11.3105</td><td>0.0</td><td>0.0</td><td>54.1851</td><td>46.2404</td><td>41.182</td><td>37.1592</td><td>36.8146</td><td>39.0444</td><td>38.5617</td><td>15.5697</td><td>13.5964</td><td>12.837</td><td>12.2169</td><td>11.7886</td><td>12.4189</td><td>9.0548</td><td>0.0</td><td>CDS-PACSxID24-1-87944</td></tr>\n",
       "<tr><td>43597</td><td>54.6425</td><td>-27.0098</td><td>-5.24445</td><td>-5.88225</td><td>-6.76164</td><td>-6.75699</td><td>-3.02663</td><td>8.63295</td><td>-1.76605</td><td>11.1267</td><td>10.3275</td><td>10.4275</td><td>10.9124</td><td>11.7341</td><td>13.5267</td><td>8.67068</td><td>0.0</td><td>58.0044</td><td>61.9884</td><td>67.8042</td><td>69.6219</td><td>71.4516</td><td>76.7022</td><td>59.5624</td><td>18.3986</td><td>15.5561</td><td>14.5839</td><td>13.954</td><td>13.7231</td><td>14.5855</td><td>11.013</td><td>0.0</td><td>CDS-PACSxID24-1-87945</td></tr>\n",
       "<tr><td>63883</td><td>54.6431</td><td>-26.8946</td><td>-1.46345</td><td>0.955591</td><td>3.84553</td><td>7.02278</td><td>9.80761</td><td>15.918</td><td>-8.29514</td><td>9.09804</td><td>8.5168</td><td>8.72154</td><td>9.15914</td><td>9.85406</td><td>11.3178</td><td>7.04574</td><td>0.0</td><td>19.6981</td><td>24.2293</td><td>25.6439</td><td>26.4106</td><td>28.554</td><td>39.0982</td><td>39.1694</td><td>15.359</td><td>12.8145</td><td>11.9521</td><td>11.4844</td><td>11.3624</td><td>12.1181</td><td>9.04582</td><td>0.0</td><td>CDS-PACSxID24-1-87946</td></tr>\n",
       "<tr><td>9212</td><td>54.6439</td><td>-26.9477</td><td>22.2942</td><td>18.002</td><td>14.5635</td><td>16.4403</td><td>17.8126</td><td>21.0479</td><td>23.6774</td><td>11.3235</td><td>10.5241</td><td>10.6373</td><td>11.0775</td><td>11.8638</td><td>13.6721</td><td>8.78468</td><td>0.0</td><td>77.1115</td><td>67.7446</td><td>68.0229</td><td>68.5946</td><td>67.9523</td><td>65.9078</td><td>53.4683</td><td>18.1003</td><td>15.4794</td><td>14.5165</td><td>13.8027</td><td>13.4961</td><td>14.1797</td><td>10.883</td><td>0.0</td><td>CDS-PACSxID24-1-87947</td></tr>\n",
       "<tr><td>38007</td><td>54.6443</td><td>-26.9441</td><td>4.17411</td><td>4.24658</td><td>1.30328</td><td>1.41589</td><td>4.67424</td><td>11.284</td><td>8.3487</td><td>11.094</td><td>10.3741</td><td>10.6084</td><td>11.1131</td><td>11.8719</td><td>13.566</td><td>8.61118</td><td>0.0</td><td>28.1981</td><td>23.8154</td><td>19.4331</td><td>16.7743</td><td>13.9798</td><td>11.2945</td><td>20.1674</td><td>18.005</td><td>15.547</td><td>14.6348</td><td>13.9185</td><td>13.5599</td><td>14.242</td><td>10.8467</td><td>0.0</td><td>CDS-PACSxID24-1-87948</td></tr>\n",
       "<tr><td>12521</td><td>54.6451</td><td>-26.8673</td><td>-34.4655</td><td>-31.9906</td><td>-28.1668</td><td>-30.8689</td><td>-37.3356</td><td>-47.0382</td><td>nan</td><td>9.99409</td><td>9.34578</td><td>9.54043</td><td>9.8956</td><td>10.5509</td><td>12.185</td><td>0.0</td><td>0.0</td><td>31.3043</td><td>26.5653</td><td>26.3379</td><td>28.5822</td><td>33.2361</td><td>45.1108</td><td>nan</td><td>16.3754</td><td>13.6032</td><td>12.6096</td><td>12.0886</td><td>11.9589</td><td>12.8544</td><td>0.0</td><td>0.0</td><td>CDS-PACSxID24-1-87949</td></tr>\n",
       "<tr><td>57428</td><td>54.6479</td><td>-26.9004</td><td>12.8088</td><td>16.3453</td><td>23.4027</td><td>31.2341</td><td>36.1239</td><td>44.6128</td><td>12.6636</td><td>9.29365</td><td>8.66666</td><td>8.82462</td><td>9.22949</td><td>9.84607</td><td>11.3037</td><td>7.11582</td><td>0.0</td><td>7.25308</td><td>6.87022</td><td>6.84851</td><td>3.47506</td><td>-1.08814</td><td>-12.8405</td><td>7.46674</td><td>15.8314</td><td>13.5912</td><td>12.7602</td><td>12.1233</td><td>11.7925</td><td>12.3422</td><td>9.14821</td><td>0.0</td><td>CDS-PACSxID24-1-87950</td></tr>\n",
       "<tr><td>53063</td><td>54.6502</td><td>-26.8884</td><td>11.9815</td><td>17.8446</td><td>24.5976</td><td>29.9165</td><td>33.2239</td><td>34.3721</td><td>3.96745</td><td>9.08998</td><td>8.5749</td><td>8.67761</td><td>9.01221</td><td>9.66609</td><td>11.0995</td><td>7.08585</td><td>0.0</td><td>26.0078</td><td>26.0833</td><td>28.5762</td><td>28.6108</td><td>29.9842</td><td>33.9096</td><td>26.4185</td><td>16.3557</td><td>14.035</td><td>13.0336</td><td>12.2927</td><td>11.9087</td><td>12.4729</td><td>9.29037</td><td>0.0</td><td>CDS-PACSxID24-1-87951</td></tr>\n",
       "<tr><td>11880</td><td>54.6507</td><td>-26.8708</td><td>3.71808</td><td>-4.14208</td><td>-3.31731</td><td>-2.81346</td><td>-8.62721</td><td>-9.37119</td><td>nan</td><td>10.7531</td><td>10.1343</td><td>10.2507</td><td>10.6596</td><td>11.3782</td><td>13.0308</td><td>0.0</td><td>0.0</td><td>36.5679</td><td>41.095</td><td>46.3715</td><td>50.8598</td><td>51.7177</td><td>50.0848</td><td>nan</td><td>17.8796</td><td>15.0778</td><td>14.0428</td><td>13.6003</td><td>13.4318</td><td>14.2455</td><td>0.0</td><td>0.0</td><td>CDS-PACSxID24-1-87952</td></tr>\n",
       "<tr><td>64057</td><td>54.6527</td><td>-26.8789</td><td>-21.9937</td><td>-23.0053</td><td>-14.8471</td><td>-1.97126</td><td>5.45709</td><td>5.22542</td><td>-6.5803</td><td>9.62018</td><td>9.06912</td><td>9.27284</td><td>9.7431</td><td>10.4056</td><td>11.8825</td><td>7.33741</td><td>0.0</td><td>4.0536</td><td>8.68244</td><td>12.4786</td><td>14.3206</td><td>15.9504</td><td>11.8151</td><td>-0.897587</td><td>17.2565</td><td>14.6522</td><td>13.6642</td><td>12.961</td><td>12.605</td><td>13.2316</td><td>9.69027</td><td>0.0</td><td>CDS-PACSxID24-1-87953</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=87953>\n",
       " XID     RA     Dec    ... Ferr_PACS_160 F_PACS_160__SKY        HELP_ID       \n",
       "        deg     deg    ...      mJy         mJy / pix                         \n",
       "int32 float32 float32  ...    float32        float32             str21        \n",
       "----- ------- -------- ... ------------- --------------- ---------------------\n",
       "20269 51.3391 -29.6326 ...       10.8895             0.0 CDS-PACSxID24-1-00001\n",
       " 9465 51.3404 -29.5366 ...       10.9488             0.0 CDS-PACSxID24-1-00002\n",
       "32785 51.3417 -29.6128 ...       11.1448             0.0 CDS-PACSxID24-1-00003\n",
       "60323 51.3436 -29.6118 ...       11.1362             0.0 CDS-PACSxID24-1-00004\n",
       "22332 51.3443 -29.5597 ...       11.1758             0.0 CDS-PACSxID24-1-00005\n",
       "18577 51.3449 -29.5822 ...       10.9573             0.0 CDS-PACSxID24-1-00006\n",
       "51911 51.3451 -29.6181 ...       11.2637             0.0 CDS-PACSxID24-1-00007\n",
       "19642 51.3451 -29.6764 ...       11.1234             0.0 CDS-PACSxID24-1-00008\n",
       "17590 51.3457 -29.5624 ...       11.4432             0.0 CDS-PACSxID24-1-00009\n",
       "  ...     ...      ... ...           ...             ...                   ...\n",
       " 5205 54.6421 -26.8742 ...        9.0548             0.0 CDS-PACSxID24-1-87944\n",
       "43597 54.6425 -27.0098 ...        11.013             0.0 CDS-PACSxID24-1-87945\n",
       "63883 54.6431 -26.8946 ...       9.04582             0.0 CDS-PACSxID24-1-87946\n",
       " 9212 54.6439 -26.9477 ...        10.883             0.0 CDS-PACSxID24-1-87947\n",
       "38007 54.6443 -26.9441 ...       10.8467             0.0 CDS-PACSxID24-1-87948\n",
       "12521 54.6451 -26.8673 ...           0.0             0.0 CDS-PACSxID24-1-87949\n",
       "57428 54.6479 -26.9004 ...       9.14821             0.0 CDS-PACSxID24-1-87950\n",
       "53063 54.6502 -26.8884 ...       9.29037             0.0 CDS-PACSxID24-1-87951\n",
       "11880 54.6507 -26.8708 ...           0.0             0.0 CDS-PACSxID24-1-87952\n",
       "64057 54.6527 -26.8789 ...       9.69027             0.0 CDS-PACSxID24-1-87953"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "old_cat"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.coordinates import SkyCoord\n",
    "from astropy import units as u\n",
    "c = SkyCoord(ra=cat['RA'], dec=cat['Dec'])  \n",
    "catalog = SkyCoord(ra=old_cat['RA'], dec=old_cat['Dec'])  \n",
    "idx, d2d, d3d = c.match_to_catalog_sky(catalog)  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x118b41550>]"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD3CAYAAAAT+Z8iAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1w1Nd97/H3PuhZK6EnJCQtAhs4BmIMBhsDEkhg7Ni4\ncSczvm2T3KZt4roZt256e5ve1G5mmrGbmyZNb907eSiJx0mdZuY6aWYSM46TMQIkHgQGHGMDB7AN\n7EpCCD0/S7v7u3/sCi9iV9pd7Wp3f/t9zXis/e1v93cOu/vZ355zfudYDMNACCGEeVmTXQAhhBCJ\nJUEvhBAmJ0EvhBAmJ0EvhBAmJ0EvhBAmZ092AULp7h6KeShQSUk+fX2j8SxOypM6Z4ZMq3Om1Rfm\nX+eKCocl1HbTndHb7bZkF2HBSZ0zQ6bVOdPqC4mrs+mCXgghxK0k6IUQwuQk6IUQwuQk6IUQwuQk\n6IUQwuRScnhlPPz4Nxc4cLodr88/UrMg185nHlJsXlOZ5JIJIcSt2s52se/oZTp6Rqkuy2fPlmVx\nzSpTBv2Pf3OBN0+6b9k2Mu7he794D0DCXgiRMtrOdt3MJgB390jcs8qUTTeH3m4Pe9++o1cWsCRC\nCDG7fUcvh9kev6yK6IxeKbUZ+LrWujFo26eAv9BabwncfhJ4CvAAz2utX1NK5QGvAIuBIeCzWuvu\nuJU+jClv+AtrO3tGEn14IYSIWMeN0FfCxjOr5jyjV0p9Cfg+kBu0bQPwOcASuF0FPANsAx4GvqaU\nygG+AJzRWjcAPwKei1vJZ2GzhrwKGIDiwuyFKIIQQkRkUZhMimdWRdJ08z7wyekbSqky4B+BLwbt\ncz9wWGs9obUeAC4B64B64FeBfV4HHoxHoeeSmz3LZcSyoJYQIh3EMavmbLrRWv9MKbUMQCllA34A\n/A9gLGi3ImAg6PYQUDxj+/S2OZWU5M9rzoexSW/Y+wZGJqmocMT83KnKjHWai9TZ/DKhvv0jkyG3\nxzOroh11sxFYCXwHf1POGqXU/wH2A8ElcgD9wGDQ9ultc5rn7G1Ul+Xj7g7dvrWkrIDu7qGYnz8V\nVVQ4TFenuUidzS9T6hsur2LJqnBfDFGNutFaH9darw10yv4+cFZr/UXgONCglMpVShUDq4F3gcPA\no4GHPwK0RFXqGO3ZsmyW++oWoghCCBGRcHkVz6yKyzh6rfU1pdSL+IPcCjyrtR5XSn0H+KFSqhWY\nBD4Vj+PNZXrs6asHLtE7OAFAqSOHJ5pWyBh6IURKmc6kfUev0NkzwpKyAvZsqYtrVlkMI/V6J+ez\n8Eim/NwLJnXODJlW50yrL8y/zhmz8IgQQohbSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJ\nSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdALIYTJSdAL\nIYTJSdALIYTJSdALIYTJSdALIYTJRbQ4uFJqM/B1rXWjUmo98G+AF5gA/lBr3aWUehJ4CvAAz2ut\nX1NK5QGvAIuBIeCzWuvuRFRECCFEaHOe0SulvgR8H8gNbPpX4C+01o3AfwF/q5SqAp4BtgEPA19T\nSuUAXwDOaK0bgB8Bz8W9BkIIIWYVSdPN+8Ang27/vtb67cDfdmAcuB84rLWe0FoPAJeAdUA98KvA\nvq8DD8al1EIIISI2Z9ON1vpnSqllQbc7AZRSW4E/B7bjP4sfCHrYEFAMFAVtn942p5KSfOx2WyS7\nhlRR4Yj5selK6pwZMq3OmVZfSEydI2qjn0kp9XvAs8AerXW3UmoQCC6dA+gHgrdPb5tTX99oLMUC\n/P9I3d1DMT8+HUmdM0Om1TmT6tt2tot9Ry/T0TNKdVk+e7YsY/OayqifJ9yXRNRBr5T6DP5O10at\ndW9g83HgBaVULpADrAbeBQ4DjwbufwRoibrkQghhYm1nu/jeL967edvdPXLzdixhH0pUwyuVUjbg\nRfxn5/+llDqglPoHrfW1wPYWYD/wrNZ6HPgOsFYp1Qr8KfAPcSm1EEKYxL6jl8NsvxK3Y0R0Rq+1\nvgw8ELhZGmafvcDeGdtGgSfmUT4hhDC1jhuhm6o7e0bidgy5YEoIIZKoujw/5PYlZQVxO4YEvRBC\nJNGeLcvCbK+L2zFiGnUjhBAiPqY7XPcdvUJnzwhLygrYs6Uubh2xIEEvhBBJt3lNJZvXVCZsSKk0\n3QghhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0Ash\nhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlJ0AshhMlFNB+9Umoz8HWtdaNSagXwMmAA7wJPa619\nSqkngacAD/C81vo1pVQe8AqwGBgCPqu17k5APYQQQoQx5xm9UupLwPeB3MCmbwHPaa0bAAvwuFKq\nCngG2AY8DHxNKZUDfAE4E9j3R8Bz8a+CEEKkN8MwuNT/IRd7PkzI80dyRv8+8EngPwK3NwIHA3+/\nDjwEeIHDWusJYEIpdQlYB9QD/xS0799HUqiSknzsdltEFQilosIR82PTldQ5M2Ranc1e30nvFEeu\nvsW+C/u50u+msrCCf9vz1bgfZ86g11r/TCm1LGiTRWttBP4eAoqBImAgaJ9Q26e3zamvbzSS3UJK\n1FJcqUzqnBkyrc5mru/g5BAt7cdocR9laGoYq8XKvYvX8fvrf2dedQ73xRjLmrG+oL8dQD8wGPh7\ntu3T24QQIiO5hzpodrfy1rXTeAwvefZcHly6gx21WynNLaGiJDFfbrEE/WmlVKPW+gDwCNAMHAde\nUErlAjnAavwdtYeBRwP3PwK0xKPQQgiRLnyGj/d6zrP/agsX+t8HYHFeOY3OejZXbSTXnpPwMsQS\n9H8N7FVKZQPngJ9qrb1KqRfxB7kVeFZrPa6U+g7wQ6VUKzAJfCpeBRdCiFQ27png2LW3OOBqpXus\nB4BVJSvY6axnbdldWC0LN7rdYhjG3HstsO7uoZgLZeZ2vXCkzpkh0+qcrvXtHe/jgPswRzqOM+YZ\nx26xsalqAzudDdQULpn1sfOtc0WFwxJqeyxn9EIIIWb4YOAK+10t/Lb7XXyGD0dWIY8u301DzQMU\nZSd39JAEvRBCxMjr83K6+wzNrlYuD14FoKZwCU3OBjZVrifLmhoRmxqlEEKINDIyNcrhjjYOuo/Q\nPzGABQt3l69mp7OBlYvuxGIJ2YKSNBL0QggRoa6R6zS7D9PW+RaTvimybdnsqN1GY+1WFudXJLt4\nYUnQCyHELAzD4HzfRZpdrbzXcx6AkpxF7HFuY+uS+8nPyktyCecmQS+EECFMeqc40XWKZlcrnSNd\nANxRXEeTs4F7ytdis8Y+TctCk6AXQoggAxODHGo/Smv7MYanRrBarGyqXE+Ts55lRUuTXbyYSNAL\nIQTgGmqn2dXKW11v4zW85NvzeKiuie01WyjJXZTs4s2LBL0QImP5DB9nbpyl2dXKxf4PAKjMr6DJ\nWc/9VRvJsWUnuYTxIUEvhMg4455xjnb6pye4Md4LwF0lK9m5tIHVpasWdHqChSBBL4TIGDfGejno\nPsyRjhOMe8exW+1sXXI/Tc56qgurkl28hJGgF0KYmmEYvD9wmWZXC7/tfg8Dg6JsBw8u3UF9zWYc\n2YXJLmLCSdALIUzJ4/Nw6vo7NLtauTrkBsBZWE2Ts4GNlfdgT5HpCRZC5tRUCJERhidHaO1o45D7\nCAOTg1iwcE/Fx2iqrWfFouUpNz3BQpCgF0KYQudIF82uVo5fO8mUz0OuLYcmZz2NtdsozytLdvGS\nSoJeCJG2DMPgbO8Fml0tnOu9AEBZbimNzm1sWXIfefbcJJcwNUjQCyHSzqR3kuPX/NMTXBu9DsCd\nxcvZubSBdeVrTDc8cr4k6IUQaaN/YoBD7qO0dhxjZGoUq8XKfZX3stNZz9Ki2mQXL2XFFPRKqSzg\nh8AywAs8CXiAlwED/8LgT2utfUqpJ4GnAvc/r7V+bf7FFkJkkiuDLppdrZy8/lt8ho+CrHw+XreT\nhtotLMopTnbxUl6sZ/SPAnat9Val1G7gBSALeE5rfUAp9V3gcaXUUeAZYBOQC7QqpX6jtZ6IR+GF\nEObl8/k4ff0Mza4W3h+4DEBVQSU7a+u5r+pesm1ZyS1gGok16C8AdqWUFSgCpoAHgIOB+18HHsJ/\ntn84EOwTSqlLwDrgxLxKLYQwrTHPGEc6TtDSdpTukR4A1pQpdtY2cFfpyowcHjlfsQb9MP5mm/NA\nOfAYsF1rbQTuHwKK8X8JDAQ9bnr7rEpK8rHbY5/ruaIiuQvxJoPUOTOYuc7Xhrt5/UIzzR8eYdwz\nQbYtiwfvbODRVU3UFi1JdvEWTCJe41iD/q+AN7TWX1ZKOYH9QPA0bw6gHxgM/D1z+6z6+kZjLJb/\nH6m7eyjmx6cjqXNmMGOdDcPgUv8H7He1cubGWQwMirOLeOiOJh5ft4vxQQMmMF29w5nvaxzuSyLW\noO/D31wD0Iu/ff60UqpRa30AeARoBo4DLyilcoEcYDX+jlohRAab8nk41fVbml0tuIY7AFjqqGWn\ns4F7F6/DZrXhyClknMwI+ESLNej/BXhJKdWC/0z+74C3gL1KqWzgHPBTrbVXKfUi0AJYgWe11uNx\nKLcQIg0NTQ7T2n6MQ+1HGZwcwoKFDRV30+Rs4I7iOml/T5CYgl5rPQz8txB37Qix715gbyzHiVXb\n2S5ebb5E75B/cI/NasHrM7BZwBvoRSh15PBE0wo2r6mM2zH3Hb1Mx41Rqsvz2bNlWdyeW4h01zF8\njWZXC8e7TuPxeci15bLLuZ0dtVspyytNdvGS7mZ+9IxSXRb//DDdBVOHTrv53i/eu2Wb1+dP9+mQ\nB+gdmri533z/QdvOdt1yTHf3SNyeW4h05TN8nO3RNLtaOd93EYDyvDKaaut5YMlGcmV6AmBh8sN0\nQf/qmxej2n/f0Svz/sfcd/Rywp5biHQz4Z2krfMkB9ytdI12A7By0R3sdDbwsfLVMj3BDK82Xwq9\n/cAlCfpwrnZF13nT2TMy72N23Ag9Sigezy1Euugb7+eg+wiHO9oY9Yxht9jYXLWRJmcDTkd1souX\nsqabmG/bPhi/60pNF/RLKx1c7hyMeP8lZQXzPmZ1eT7u7ttDPR7PLUSq+3DgKs2uFk53n8Fn+CjM\nKuCRZQ/SULOF4hzzjvtPJ6YL+id2reQbr5yMeP89W+rmfcw9W5bd1i8Qr+cWIhV5fV7e7n6XZlcr\nHw5eAaC6oIomZwP3Va4nS6YniFipIyfkWX2pIyduxzBd0G/fUMvg4Pgto27CKXXkxKUNbPo59h29\nQmfPCEvKCtizpU7a54XpjE6NcbijjYPuI/RN+K99/FjZapqc9aiSFTI8MgZPNK0IeaL4RNOKuB3D\ndEEP/uCdDtmZPdrBBkYmE3JMIczm+mg3za7DHLv2FpPeSbKtWWyv2UqjcxuV+RXJLl5aW4gTRVMG\nfbDNayrZd/SytKELESXDMLjQ9z77XS2813MeA4NFOcU8uuxBtlXfT35WfrKLaBrTJ4qJmubC9EEP\n0oYuRDSmvFO81fU2ze5W2oc7AVhetJQmZz3rK+7GZo19wkGRHKYP+ukrziyA3WbF6/NRXV4obehC\nzDA4OUSL+ygt7ccYmhrGarGycfE9NDnrWV4sJ0XxFupq+sd2JGaUkqmDfmb7/JTXByAhL0QQ91AH\nza5W3uo6jcfwkmfPY/fSRrbXbqE0tyTZxTOlcFfDFhXlsro2/itmmTro5YpVIULzGT7evXGOZlcr\nF/rfB2BxXjlNznrur9pIrj1+Q/vE7cJl06tvXuQrn90U9+OZOujlilUhbjXumeBY51sccLfSPeZf\nvUmVrGCns4E1ZUqmJ1gg4bLpyrXIL/aMhqmDXq5YFcKvZ6yPg+2HOdJxnDHPOHarnS1L7qPJWU9N\nYeas3pQqwmWTYfibdeLd4mDqoJfRNiKTGYbBh4NX2H+1hbe738XAwJFdyJ7lu2mo2YIjuzDZRcxY\n4bIJEtO0bOqglytWRSby+rycvv4O+92tXBl0AVBTuISdzgY2Vq4ny2rqj31a2Lymkn//5XsYxu33\nJaJp2fSvuFyxKjLFyNQoh9vbONh+hP6JASxYWFe+liZnPSsX3SHTE6SYmvKCBWtaNn3QC2F210au\n0+xupa3zJFO+KXJs2TTWbmNH7TYW55cnu3gijIVsWpagFyINGYbB+b6L7He1cLZHA1CaW8KO2q1s\nXXI/+Vl5SS6hmEuopuU/eFil1jh6pdSXgU/gXxz828BB4GXAAN4FntZa+5RSTwJPAR7gea31a/Mt\ntBCZatI7xYmuUzS7Wukc6QLgjuJlNDnruad8rUxPkGZmNi2n1Fw3SqlGYCuwDcgH/ifwLeA5rfUB\npdR3gceVUkeBZ4BNQC7QqpT6jdY6fkunCJEBBiYGefNMM7++eIjhqRGsFiubKtez09lAXZEz2cUT\nKS7WM/qHgTPAz4Ei4G+AJ/Gf1QO8DjwEeIHDgWCfUEpdAtYBJ+ZTaCEyxdUhN82uVk52/Rav4aXA\nns9DdU1sr9lCSe6iZBdPpIlYg74cqAMeA5YDvwCsWuvpwUJDQDH+L4GBoMdNb59VSUk+dnvsP0Er\nKjJv+TKps3n4fD7e6niHfRf2c67bv9h9jaOKR1ftZPuyzeTYs5NcwoVj1td4Nomoc6xB3wOc11pP\nAlopNQ4E/350AP3AYODvmdtn1dcX+vLgSCSqjSuVSZ3NYcwz7p+ewNXKjfFeAFaXrqLJ2cDq0pVU\nLi4O1DkzWj7N+BrPZb51DvclEWvQtwJ/qZT6FrAEKADeVEo1aq0PAI8AzcBx4AWlVC6QA6zG31Er\nhAi4MdbLQbd/eoJx7wRZVjvbqu+nsbae6sKqZBdPmEBMQa+1fk0ptR1/kFuBp4EPgb1KqWzgHPBT\nrbVXKfUi0BLY71mt9Xh8ii5E+jIMg/cHLtPsauG33e9hYFCc7WB3XSP11Q9QmC3zMUUr1PzucrGk\nX8zDK7XWXwqxeUeI/fYCe2M9jhBm4vF5OHX9HZpdLVwdagfA6ahhp7OBexevwy7TE8Qk3PzugIQ9\ncsGUEAtieHKE1o5jHHIfYWByCAsW1ld8jCZnA3cWL5PpCeZJ1p6YnQS9EAnUOdJFs6uF49dOMeXz\nkGvLYaezgR212yjPK0128UxjrrUnMr1ZR4JeiDjzGT7O9V6k2dXCud4LAJTlltLkrOeBJZvIs+cm\nuYTmM9vaE9Ksk6FBn+nf7iIxJr2TtF3zT0/QNXodgBWLlrPT2cDd5Wtk9aaARHz+ZpsgTJp1MjDo\n5dtdxFv/xAAH3Uc43N7GiGcUm8XG/VX30uSsZ6mjNtnFW1BzhXiiPn+zrT2x95dnQz4mk5YUzbig\nl293ES9XBl3sd7Vw6vo7+AwfhVkFfHzZLhpqHmBRTvxnIEx1kYR4Ij9/4daekCVFMzDoM33BcGm2\nmh+vz8s7N86y39XCBwOXAVhSUEmTs577Ku8l25aV3ALOItGvfSQhnozPnywpmoFBn8nf7tJsFbsx\nzxhHOk5wwH2Y3vE+ANaW3UWTs567Slam/PDIhXjtIwnxZHz+ZEnRDAz6TP52l2ar6HWP9tDsbuVY\n5wkmvJNkWbOor3mAptp6qgoWJ7t4EVuI1z6SEE/W5y/TlxTNuKDP5G/3TG+2ipRhGFzs/4BmVytn\nbpzFwGBRTjEfr9vFtprNFGTlJ7uIUVuI1z6SEM/kz18yZVzQQ+Z+u2dys1UkpnweTna9TbOrFfdw\nBwB1Dic7nfVsWLwurVdvWojXPtIQz9TPXzKZPujbznbxavMleoc+mtq11JHDE00r0v7NdrNzrWeU\n6rK5O9cyudlqNkOTw7S0H+VQ+1GGJoexYGHD4nXsdNazvKguqvb3VO3sTsRrH66uqVBfcStTB/3M\nDqhpvUMTad8JGUvnmvxsvlX7cCfNrlZOdJ3G4/OQZ89l19Lt7KjZRlleSdTPl8qd3fF+7VO5ruJ2\npg76cB1QH92fvp2QsXauZfoZl8/wcbZHs9/Vgu67BEBFXhmNznoeqNpI7jymJ0j1zu54vvapXldx\nK1MHfbgOqGnp3AkpHavRmfBO0tb5Fs3uVq6P3gBg1aI72bm0gbVld8VleoJMek0yqa5mYOqgD9cB\nNS2dOyGlYzUyfeP9HHQfobWjjTHPGHaLjQeqNtHkrKfWUR3XYy3Ea5IqfQDy/ksvpg76cB1Q00bH\np/j815tTqtMsUtKxOrsPB67Q7GrldPeZm9MTPLrsQeprtlCck5gFp+P9mswM9fVqMa+1fnjz/mS2\ni8v7L72YOuin3/wzR90U5mUxPDZ1c1s6diRJx+rtvD4vb3efodnVyoeDVwGoLqhip7OBTZXryUrw\n9ATxfE1CdXa6uz8MuW8y2sXl/ZdeTBv0M8+GNqyqQF/to+PGKBOTnpCPSbeOpOnOtfmuHJ/uRqdG\nOdxxnIPuI/RN9GPBwt3lq2mqbWBVyZ0LOj1BvDo8X22+FPG+sbSLx6MJKNM79tPJvIJeKbUYOAns\nBjzAy4ABvAs8rbX2KaWeBJ4K3P+81vq1eZU4AqHPhj76MPi8oR8nHUnppWu0mwOuVo51vsWkb4ps\naxbba7bS5NzG4vyKZBcvZm1nu275BTqXaNvFZWhk5ok56JVSWcD3gLHApm8Bz2mtDyilvgs8rpQ6\nCjwDbAJygVal1G+01pG/i2Mw17DKcKQjKfUZhoHuu8QPzh/jVMcZAEpyFvFo7Va2Vd9PfhpOTzBT\ntO/faNvFZWhk5pnPGf03ge8CXw7c3ggcDPz9OvAQ4AUOB4J9Qil1CVgHnJjHcec017DKcKQjKXVN\neac40fU2za4WOkauAbC8qI4mZz3rKz6W1tMTTJtuTpltpNiujbXoq/3zaheXoZGZJ6agV0r9EdCt\ntX5DKTUd9BattRH4ewgoBoqAgaCHTm+fVUlJPnZ77B/c0uJcbvSPzblflt2Kz2dQUuS/SGbva2d5\n44SLJ3atZPuG9FoZqKIiMSNJkq1/fJBfXzrIry8dYnBiGKvFytalm9izaicry5Ynu3hhHTrt5tU3\nL3K1a4illY4531OHTrtnHSEGUL4ojy9+auO8y7a0ysHlzsHbtjsrHSn3Pkq18iyERNQ51jP6PwEM\npdSDwHrgR0DwnK0OoB8YDPw9c/us+vpiOyMHOOceiCjkARx5WWxYVcGbJ903t13uHOQbr5xkcHA8\nbX7GmrEz1jXUQbOrhZNdb+MxvOTb89i9tJEdtVspyV1ERVlq1Tm4c3NRYfYtbezT76lT57r49O5V\nIR//kzfOz3mMP35sTVzq/PB9zpBfKg/f50ypf1Mzvq/nMt86h/uSiCnotdbbp/9WSh0A/gz4hlKq\nUWt9AHgEaAaOAy8opXKBHGA1/o7ahHn1zYsR79s7NHFLyAeT9sqF5zN8vHvjHPtdLVzs/wCAxfnl\nNNU2sHnJRnJs2Uku4UdmC/ZwHalvnnSzoqY45PtqtubG2opC9mypY/uG2rgEnwyNzDzxHF7518Be\npVQ2cA74qdbaq5R6EWgBrMCzWuvxOB7zNle74nMGIO2VC2fcM86xzpM0u1u5MdYDwF0lK2ly1rOm\nTMVleoJ4mjlqJZoRMuFOIMJdaVpbUchXP3d/bAWdhQyNzCzzDnqtdWPQzR0h7t8L7J3vcSK1tDJ0\n+2O0ZARO4vWM9XLQfYQjnccZ84xjt9rZuuQ+mpwNVBdWJbt4YcU6qgvCn0DIlaYikUx3wdQTu1by\njVdOzvt55AOWGIZh8MHAFZpdLbzd/S4GBo7sQh5b/hD1NQ/gyC4M+9jg5pKlVQ4evs8JsOBzv8Q6\nqgvCn0BIc4pIJNMF/fYNtXzzxycxjLn3DWfXxtqIP2CpMslUqvP6vJy6/g7NrlauDLkAqC2sZqez\ngXsr7yHLOvtbcWZzyeXOwdvOgBfqwp+5JsubzWwnENKcIhLFdEEPUFKYE1W76Uz66pwDgwC5wjAS\nI1OjHG5v42D7EfonBrBg4Z7ytTQ561mx6I6IpyeIprkk0R3pc02WZ7NY8IY40ygtypH3hUgKUwb9\nfEXaEStXGIZ3beQ6ze5W2jpPMuWbIseWTWPtNhpr66nIL4v6+aJpLkl0R/r0a/vSvnNMeX2372DB\nPxHIDE80rkhouYQIx5RB3z88Oa/HFxdENoxPrjC8lWEYnO+9yH53C2d7NABluSU01m5jS/V95Nnz\nYn7uaJpLFqojPWTIA16fP+VLHTkMjExKe7tIOlMG/cxxzVGLcLJDWXzBb9I7xYlrp9jvbuXaSBcA\ndxYvY6ezgbvL18RleoK5mktu3Te6jvRo+1nCrUU8U35uFt98eltUZREiEUwX9IdOu+cX8sBAhL8I\nMn1I3MDEIIfcR2jpOMbI1ChWi5X7KjfQ5KynrsgZ12PNHJXirAwedRPZSJVQgQ5E3c8SaX9Bpv6y\nE6nHdEEfzZWxNgt4Q7SlRnpGnqlD4q4Oudl/tZVT13+L1/BSYM/n4bqdbK/dwqKcOacyilnwqJSK\nCgevHbx0y9Wpo+NT/Psv3uOlfWfx+Axqygtunp2H6zgvdeSEPNZs/SyR9hdk2i87kbpMF/RXrkV+\nsVSokIfozsgzZUicz/Dxzo2z7L/awvsD/pWOqvIX0+Ss5/6qe8mOYXqCWIamTj+m/cbILUNog3/F\nTQVe2OCz83Bn4eF+/XXcGA5bhkj7CzLll51IfaYLervNypQndCfZXEqLcniicUVGBHekxjzjHO08\nwQHXYXrGewFYXbqKnc4GVpeuinn1pliGpkbaNj7Tqwcu0TsYXXOezRp+2oVwTXbS+SpSlemC3hNm\nJEQk8nOy5MMZcGOshwPuwxztOMG4d4Isq51t1ZtpctazpCBxS+WFajJpO9t127q/0Yg25AE8vvDv\no0xtshPpy3RBX1acF/E0xTNleueZYRhc6v+QZncr73S/h4FBcbaD3XVN1NdspjArPm3Osy2VN/M1\niPUsPlKljtAX19WUh5+KATKnyU6Yg+mCfj4ytfPM4/Nw6vo77He14BpqB2Cpo4YmZwP3Ll6HPcz0\nBLFO/zChz3ZFAAAQEklEQVTbqJWZr8F8JhCby1OfWAuQ0SOnRGYwXdDHejYPmffhHp4coaX9GIfa\njzA4OYQFC+sr7qbJWc+dxctmbX+fz/QP7bN0ZM58DdpvJOZXVm1F4S3llGYYYWamC/r5mD57NPuH\nvGP4GgfcrRy/doopn4dcWy47nQ3sqN1GeV5pRM8R6/QPbWe7Qs0OAPjniAH4yg/abg6ZnGtyul0b\na8MuHjOb4C8UaYYRZidBH8TMk5L5DB/nei+w/2oL5/v81xqU55bS6KzngSWbyLPnRvV8kU7/MNtK\nTDN5DSOqBT2sFm4uzTdb2Fss/jZ3OWMXmUqCPgQzTUo26Z2k7dpJml2H6Rq9DsDKRXfQ5Gzg7vLV\nMa/eFG4seXFhNl/5QRvt3SNYZ1yQNt8rlm8vg7/DVF/tm3W/mvLErNIkRLowVdC3ne2Ky/PMdlY6\nW6djKs1N3z8xwEH3EVrbjzHqGcNmsbG5aiNNznqcjpqoyhxqn3BjyXsHJ+jFH+jhLkiLl+nml7mu\nVM20vhchZjJV0MdrhEbwyI9IOx2j2S+RXwZXBl3sd7Vw6vo7+AwfhVkFPLJsFw01WynOuXWF+EjK\nPNvUAbs21qKv9tPZM0JxYTZ9gxNh298jkWWzhp0RMlhpUQ6f+8THWF3rn24h3K+LLJuVP9mz2jS/\nzoSIVUxBr5TKAl4ClgE5wPPAWeBl/DNxvws8rbX2KaWeBJ4CPMDzWuvX5l/s0OazxFuw4DPASDsd\nI9kvUQuVfLR6UwsfDFwBoLqgiiZnPZsqN5Bty7pl/7kuQAouc7gLm3qHJnjzpJunPrGWS+0DMXWI\nzrR9fXXI5wl1xWlFhYPubv9C8OF+XUjIC+EX6xn9Z4AerfV/V0qVAm8H/ntOa31AKfVd4HGl1FHg\nGWATkAu0KqV+o7WOb2NtwHynJ7ZZLDTeW3NLOIQb3ufuHuYrP2i7eUYeSefkfEaqhPoVMDo1xpHO\n47QeO0r3qH96grVld7HT2YAqWRFyeGQkFyC13xi+OfLFN8ewl3hdzFRalMOnd69iRU1x1EMd5UpV\nIWYXa9C/Cvw08LcF/9n6RuBgYNvrwEOAFzgcCPYJpdQlYB1wYrYnLynJx26Pfg5zqy22jsVpXsPg\nzZNu7l1dyfYNtcDsc+dMn5EXFeWytMrB5c7bJ1RzVjqoqPA3mXT0hP8ymN5npkOn3bf9Cvj3N47z\n9vgYF4bPMO6ZIMeWzUMrtvPoyiaqi6rCPs+rb14MWcaZDIOY10SN1cDwJBUVDh7b4eCxHZGtxBT8\nbxbN49JZuPeJWWVafSExdY4p6LXWwwBKKQf+wH8O+KbWevr0bwgoBoqAgaCHTm+fVV9fbE0wvQPj\nMT1upp+8oW+2/3oimCDtB794N+TScQAP3+e82cRQXRZ+oZLpfW4vy/nAXwZWRy/2qitYF13nnX5Y\nlFPMw3U7efzuXYwN+mCCkM/z499ciEvTSiLN9m8QSnDTTabItDpnWn1h/nUO9yURc2esUsoJ/Bz4\nttb6P5VS/xR0twPoBwYDf8/cnhDRLDc3m+DmlpqKgjmfM9SkWaFmwoxloZKOniFs5R3YK69gLfC/\nAXzDxUxdW8Zj25rYUldNYU4BY9z65mg728Urv9aMjHtmLXuqkJExQiROrJ2xlcCvgT/XWr8Z2Hxa\nKdWotT4APAI0A8eBF5RSufg7bVfj76hNiGiWm5uN12fwuf+9n5qKAtTSkpi+PELNhBlNW/LQ5DCH\n2o+Su+EQhm0Cw7Dg6anC21WHb7gEgJ8d/JAta6tvedx8Z3pMhl0ba6U9XYgEshhzXWMeglLqX4Hf\nA84Hbf5L4EUgGzgHPKm19gZG3fwpYAX+UWv9s7mev7t7KOZRep/7+v45L5uP1vQwwvYbwxE/t9UC\n1eUFUQ+jbB/upNnVyomu03h8HvDambruxNu1FGPy9sW1Sx059I9MUl2WT3FhDu992Btt9RaMzWoh\nN9t281fGfOb/l5/15pdp9YW4NN2EnKAqpqBPtFiDPlFT2tZWfHRlpX8EzEdn5KPjUxGfPT/1ibUh\nQ81n+Hiv5zz7Xa1c6PMPZ3TYSuh5vwrvjRrwpc/lDrUVhaili0L2CYSrfywkBMwv0+oLiQv69EmQ\nCLz8+rmEPG970LJyMyfAiqajc+YwynHPBG3XTnLA1cr1sRsAFBvV3Li0hOt95fgHNKWP4C/EWIZJ\nCiESw1RBPzEV++pSszEM/5l8qKkBohnNMt3J2zvex0H3EQ53HGfMM4bdYuOBJZsYczs5djL2aZaT\nTWaEFCI1mSroE+l7v3iPfUcv39LWHu2UC+XV4/zg3Vc4ff0MBgbGVDb5w6t5bNUOGlffwVOvNce9\n3AtBphoQIrVJ0Edh5pQFEU25YPFhK+nCXnWZwcIBTl0H36gDz7U6vD3VjBtWfnTxMj/adzmxhY+D\ncHO/S8gLkdok6GPw0r5z7P3lWWxW8HnD7GSbwl7hwlZ5FWvOOIYB3r4KPNeW4RsqJZXb321WC4Zh\nYLNa8fh81JQX3mxjl7Z3IdKPBH0MpmdYDBXyltxh7JVXsZW3Y7F5Mbw2PNeW4umqw5hIjzVpP//Y\nmrDhLW3vQqQfUwW9zWrB60vGcFEDa1EP9qor2BZ1A+CbyMXTvgJPdy14s+Z4fHJMTwucZbPi9fmo\nDjpzF0KYh6mCfsFD3uLFVtaJveoy1nz/EEzv0CJ/80zfYvzXiKUmm9XC9/6mMdnFEEIsAFMF/YLJ\nmsC++Cr2xS4sWZMYPgueniV4rtVhjCxKduluEe5XTvDiKkIIc5Ogj4IlfxB75WVsZZ1YrAaGJ4up\njuV4uupgKrrFtRPBbrNQVJDNwPBHi3RA6DnjZRIxITKHBP2cDKyLrvvb34v888j4xgqYCgyPTJXp\nCeaaXkBGygiRuVIjpVKR1YOtwo298irWXP94ee9Amb/9fSB1pieYuX5qKMkeKZNKi6YLkYkk6Gew\nZI9iq7yKvcKNxe7B8FnxXK/1D48cS43Vbgry7Hxmt7oZlqk8+VOi1skVQkROgh4AA2thv3/0TEkX\nFgsYkzlMuZfjue4ET3ayC8iujbV8eveqZBcjarGukyuEiJ/MDnqLD1vpNeyVl7EW+tdS9Y0U+dvf\ne5eAkZzhkQV5dkbGAnO2O3J4oim2OdtTQSSLpgshEss0Qd92tivyne2T2Ctc2CuvYsme8E9P0FuJ\np6sO31AJC93+brH4jzjzgqXptu29vzx724Rq6SLc8o4yvFOIhWOaoI9kLnpL7rB/9ExZOxabLzA9\nQV1geoL8BSjlrQpy7XzmIRVy+uOZywHG0radCp2gsayTOx+pUGchUo1pgj78XPQG1uIb2CuvYFvk\nX9zDN57HlLsOb3ftgg6PjGTpvLlWyQo1XXIkz5OsTtBo1smdr1SpsxCpJuEpp5SyAt8G7gEmgM9r\nrS8l+rhYvNjKO7BXXvloeoLBEjxd09MTLFzzjM1qYe+XmiLaN5I57iMJsFTqBF2o4Z2pVGchUslC\nnM7+LpCrtd6ilHoA+Gfg8YQdLWs8aHqCKf/0BDeq/dMTjIYfa55I0bRHRzTHfcBsAZaJnaCZWGch\nIrEQQV8P/ApAa31MKbUpUQey12rsVZf90xNMZTHVfgee60uTPj1BNO3R4TovQ5ktwCLpBDVbe7Z0\n/AoR2kIEfREwEHTbq5Sya6094R5QUpKP3W6L+kDWgkGM8QKmuqanJ4j+OWK1bIkDsHDl2iB2mxWv\n18fSqiKe2LWS7RtqI36eP3j4Lr7xysmI9nVWOqio8F/ENf3/uZ7nDx5WVFQ4OHTaHbI9u6goN6ry\nJlO0dTYDs9QjUplWX0hMnRci6AeB4JJbZwt5gL6+yJsvgk3q+2J6XDx4PAZf/Vzo40dz1erq2mKe\n+sTaWzov1dJFIZfwe/g+J93dQyGvjA31PHu21LG6tpju7iF+8sb5kMf/yRt61ukUUkUsdU53qXwF\ndCJkWn1h/nUO9yWxEEF/GPgd4P8F2ujPLMAxF1w824FDdV7GsoTfbJ2gZm3PTva8PkKkooUI+p8D\nu5VSR/APdfnjBTjmgkt0O3C8A0zas4XIHAkPeq21D/izRB8n2dJtfveFvpBJCJE8prlgKlkiuQgq\nFS3khUxCiOSSoJ+H2opCvvq5+5NdjJhJe7YQmSF1V69OA9LMIYRIB3JGPwuLBWoCM0qCNHMIIdKT\naYI+y2ZhymvE5bkK87L49O5VtwW5BLsQIh2ZJui3r68JeVHRXNYuL2VgeFLO1IUQpmWaoJ9eZi/S\nsK+tKJRQF0JkBNMEPfjD/ouf2phxl00LIcRsZNSNEEKYnAS9EEKYnAS9EEKYnAS9EEKYnAS9EEKY\nnMUw4nORkRBCiNQkZ/RCCGFyEvRCCGFyEvRCCGFyEvRCCGFyEvRCCGFyEvRCCGFyEvRCCGFyppi9\nUillBb4N3ANMAJ/XWl9KbqniRymVBbwELANygOeBs8DLgAG8CzyttfYppZ4EngI8wPNa69eSUeZ4\nUEotBk4Cu/HX52VMXF8ApdSXgU8A2fjf0wcxab0D7+sf4n9fe4EnMfHrrJTaDHxda92olFpBhPVU\nSuUBrwCLgSHgs1rr7miObZYz+t8FcrXWW4D/BfxzkssTb58BerTWDcDHgf8LfAt4LrDNAjyulKoC\nngG2AQ8DX1NK5SSpzPMSCIHvAWOBTaauL4BSqhHYir8+OwAn5q73o4Bda70V+CrwAiatr1LqS8D3\ngdzApmjq+QXgTGDfHwHPRXt8swR9PfArAK31MWBTcosTd68Cfx/424L/234j/rM9gNeBB4H7gcNa\n6wmt9QBwCVi3wGWNl28C3wU6ArfNXl/wf7jPAD8Hfgm8hrnrfQGwB36RFwFTmLe+7wOfDLodTT1v\n5lvQvlExS9AXAQNBt71KKVM0SwForYe11kNKKQfwU/zf6Bat9fT8FUNAMbf/O0xvTytKqT8CurXW\nbwRtNm19g5TjP0l5Avgz4MeA1cT1HsbfbHMe2Au8iElfZ631z/B/kU2Lpp7B22Oqu1mCfhBwBN22\naq09ySpMIiilnEAz8B9a6/8EfEF3O4B+bv93mN6ebv4E2K2UOgCsx/9zdXHQ/War77Qe4A2t9aTW\nWgPj3PqhNlu9/wp/fVfh71/7If6+iWlmq2+waD6/wdtjqrtZgv4w/vY+lFIP4P/5axpKqUrg18Df\naq1fCmw+HWjTBXgEaAGOAw1KqVylVDGwGn9HT1rRWm/XWu/QWjcCbwN/CLxu1voGaQU+rpSyKKWq\ngQLgTRPXu4+PzlR7gSxM/L6eIZp63sy3oH2jYpbmjZ/jPwM8gr8N+4+TXJ54+zugBPh7pdR0W/1f\nAi8qpbKBc8BPtdZepdSL+N8IVuBZrfV4Ukocf38N7DVzfQMjLLbj/8BbgaeBDzFvvf8FeEkp1YL/\nTP7vgLcwb32DRfx+Vkp9B/ihUqoVmAQ+Fe3BZJpiIYQwObM03QghhAhDgl4IIUxOgl4IIUxOgl4I\nIUxOgl4IIUxOgl4IIUxOgl4IIUzu/wPpxGyjNgDelAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x118b41320>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(cat['F_PACS_100'],old_cat['F_PACS_100'][idx],'o' )\n",
    "plt.plot(np.arange(0,1000),np.arange(0,1000))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x118d2f860>]"
      ]
     },
     "execution_count": 33,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD3CAYAAAAT+Z8iAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAH51JREFUeJzt3X+wlNWd5/F333uRK9AgYIviVZkq9DtoVSKLIyFRxCWu\nozg6m1pq83PMJCpxsjGZmRpnN5pkkjIz6/pjKrobx6CsGpNsljjWGihjHEJEGIIJwSoV81UyokEI\ndfl9AUF+9P7R3dfm3ufpH0//vKc/ryqK208/3c8593Z/+vR5znNOKpvNIiIi4epqdQFERKSxFPQi\nIoFT0IuIBE5BLyISOAW9iEjgelpdgCj9/QOJhgJNnDiG3bsP1rs4bU117gyqc2eotc6ZTDoVtT2o\nFn1PT3eri9B0qnNnUJ07Q6PqHFTQi4jIcAp6EZHAKehFRAJX0clYM5sN3Onu88zs/wCn5++aBvzC\n3T9qZt8CLgEG8vddB7wLPA6clt9+vbv317H8IiJSRtmgN7NbgU8BBwDc/aP57ROBlcBf5nedBVzp\n7juKHvtXwEvu/ndm9lHgduCLda2BiIiUVEmL/rfAR4DvDtn+deB+d99mZl3AucB3zGwK8LC7LyHX\nwv8f+f2fBr5SSaEmThyT+OxzJpNm1YYtLF3xOm9tH2DS+F4Adu07xNlT0iycfy7A4P2FbXNn9iU6\nXjvIZNKtLkLTqc6doVPqXJxZjcikskHv7k+Y2bTibWZ2GjCf91rzY4H7gXuBbmClmf0KGA/sze8z\nAEyopFBJx5FmMmmWPbeJB596ZXDbjj3vDP68eds+7np8/QmPKWzbt+8Qs8+fkui4rZTJpOnvHyi/\nY0BU587QKXVet3H7CZlVSybFfTAmPRn7n4Dvu/ux/O2DwLfc/aC7DwA/A94P7AMKR04DexIer2LL\n125O+Lg361oOEZFKxGVWPTMpadB/mFxXTMF5wBoz6zazUeS6bH4NrAGuzu9zFfB80oJWauuOZN8G\ntu08UOeSiIiUF5dZ9cykpEFvwL8Vbrj7q+T68H8BPAc85u6vAA8AF5jZauAmcv36DTX11DGJHnfG\n5LF1LomISHlxmVXPTKpoeKW7bwY+UHT7goh97gLuGrLtILCwtiJWZ8GcaSf0d1X+uHMaUBoRkdLi\nMquemdSWk5rVonDyYvnaN9m28wATxp4EKdi7/13OmDx28JdXuL+wbSSeiBWRkW9oZjUik4ILesj9\n4sr9khTsItIuCpnVqJFGmgJBRCRwCnoRkcAp6EVEAqegFxEJnIJeRCRwCnoRkcAp6EVEAqegFxEJ\nnIJeRCRwCnoRkcAp6EVEAqegFxEJnIJeRCRwCnoRkcAp6EVEAhfkfPTrNm5n6cpN7Bo4DMCk9GgW\nXj5dc9CLSFtat3E7y9duZuvOg0ydPIYFc6Zp4ZFS1m3cPmxZrl0Dhwe3KexFpJ0Mzawt/QfqnlcV\nBb2ZzQbudPd5ZjYTWAa8nr/7AXf/oZndCCwCjgJ3uPsyMzsZeBw4DRgArnf3/rqUPMbytZtL3Pem\ngl5E2kpcZtUzr8oGvZndCnwKOJDfNAu4193vKdrndOAW4CKgF1htZs8CNwMvufvfmdlHgduBL9al\n5DG27jgYe9+2nQdi7xMRaYW4zKpnXlVyMva3wEeKbs8CFpjZKjN72MzSwMXAGnc/7O57gU3A+4BL\ngJ/kH/c08OG6lTzG1FPHxN53xuSxjT68iEhV4jKrnnlVtkXv7k+Y2bSiTS8AD7n7ejO7Dfga8CKw\nt2ifAWACML5oe2FbWRMnjqGnp7uSXYf52JV/yF2Pr4+5z8hk0omet52FWKdyVOfO0Al1jsuseuZV\nkpOxT7r7nsLPwP3AKqC4RGlgD7CvaHthW1m7d8d3v5SSyaSZ0TeBRddewNKfb2LXvhNH3czom9CQ\nFdZbqVGrxrcz1bkzdEqdC5m1fO2bbNt5gDMmj2XBnHMS5VXcB0OSoH/GzL7g7i8A84H15Fr53zSz\nXmA0MAN4GVgDXJ2//yrg+QTHq9rs86fopKuIjBiFzGrUh1uSoL8ZuN/MjgC/B25y931mdh+5IO8C\nbnP3Q2b2APComa0G3gU+Xq+Ci4hIZVLZbLbVZRimv38gUaE65ateMdW5M6jOnaHWOmcy6VTUdk2B\nICISOAW9iEjgFPQiIoFT0IuIBE5BLyISOAW9iEjgFPQiIoFT0IuIBE5BLyISOAW9iEjgFPQiIoFT\n0IuIBE5BLyISOAW9iEjgFPQiIoFT0IuIBE5BLyISOAW9iEjgFPQiIoGraHFwM5sN3Onu88zsQuB+\n4BhwGPgzd99uZt8CLgEKCx5eR25B8MeB0/Lbr3f3/jrXQURESijbojezW4GHgN78pm8BX3D3ecA/\nA3+b3z4LuNLd5+X/7QVuBl5y90uBx4Db61x+EREpo5Kum98CHym6/VF3fzH/cw9wyMy6gHOB75jZ\nGjP7TP7+S4Cf5H9+GvhwHcosIiJVKNt14+5PmNm0otvbAMzsg8B/AeYCY8l159wLdAMrzexXwHhg\nb/6hA8CESgo1ceIYenq6K69FkUwmnehxI5nq3BlU587QiDpX1Ec/lJn9Z+A2YIG795tZN/Atdz+Y\nv/9nwPuBfUCh1GlgTyXPv3v3wSTFIpNJ098/UH7HgKjOnUF17gy11jnuQ6LqoDezTwKLgHnuviu/\n+Tzgh2Y2k1x30CXAo+ROwl4NvABcBTxfdclFRKQmVQV9vuV+H/AW8M9mBvCcu3/NzL4L/AI4Ajzm\n7q+Y2RvAo2a2mtwInI/XtfQiIlJWKpvNtroMw/T3DyQqlL7qdQbVuTN0Up3XbdzO8rWb2brzIFMn\nj2HBnGnMPn9K1c+TyaRTUdsT9dGLiEh9rNu4nQefemXw9pb+A4O3k4R9FF0ZKyLSQsvXbo7Z/mbd\njqGgFxFpoa07okcZbtt5oG7HUNCLiLTQ1FPHRG4/Y/LYuh1DQS8i0kIL5kyL2X5O3Y6hk7EiIi1U\nOOG6fO2bbNt5gDMmj2XBnHPqdiIWFPQiIi03+/wpzD5/SsOGlKrrRkQkcAp6EZHAKehFRAKnoBcR\nCZyCXkQkcAp6EZHAKehFRAKnoBcRCZyCXkQkcAp6EZHAKehFRAKnoBcRCVxFk5qZ2WzgTnefZ2bT\ngUeALPAy8Hl3P25mNwKLgKPAHe6+zMxOBh4HTgMGgOvdvb8B9RARkRhlW/RmdivwENCb33QvcLu7\nXwqkgOvM7HTgFuBDwJXAP5jZaOBm4KX8vo8Bt9e/CiIiUkolXTe/BT5SdHsW8Fz+56eBDwMXA2vc\n/bC77wU2Ae8DLgF+MmRfERFporJdN+7+hJlNK9qUcvds/ucBYAIwHthbtE/U9sK2siZOHENPT3cl\nuw6TyaQTPW4kU507g+rcGRpR5yQLjxwv+jkN7AH25X8utb2wrazdu6MXyy2nUZP2tzPVuTOozp2h\n1jrHfUgkGXWzwczm5X++CngeeAG41Mx6zWwCMIPcido1wNVD9hURkSZKEvR/DXzdzNYCJwE/cvff\nA/eRC/KfAbe5+yHgAeACM1sN3AR8vT7FFhGRSqWy2Wz5vZqsv38gUaH0Va8zqM6dQXVO9PhU1HZd\nMCUiEjgFvYhI4BT0IiKBU9CLiAROQS8iEjgFvYhI4BT0IiKBU9CLiAQuyVw3I9K6jdtZvnYzW3cc\nZOqpY1gwZxqzz5/S6mKJiDRcRwT9uo3befCpVwZvb+k/MHhbYS8ioeuIrpvlazfHbH+zqeUQEWmF\njgj6rTuipz3etvNAk0siItJ8HRH0U08dE7n9jMljm1wSEZHm64igXzBnWsz2c5pbEBGRFuiIk7GF\nE67L177Jtp0HOGPyWBbMOUcnYkWkI3RE0EMu7BXsItKJOqLrRkSkkynoRUQCp6AXEQmcgl5EJHCJ\nTsaa2aeBT+dv9gIXAnOAZcDr+e0PuPsPzexGYBFwFLjD3ZfVUmAREalOoqB390eARwDM7H8BS4BZ\nwL3ufk9hPzM7HbgFuIjcB8JqM3vW3Q/XVmwREalUKpvNJn6wmV0E3O3u88zsAcDIfXi8DnwJuBy4\n2t0/l9//SeDv3f2XpZ736NFj2Z6e7sTlEhHpUKmojbWOo/8y8PX8zy8AD7n7ejO7Dfga8CKwt2j/\nAWBCuSfdvTt6bppyMpk0/f0DiR47UqnOnUF17gy11jmTSUduTxz0ZnYKYO6+Mr/pSXffU/gZuB9Y\nBRQfOQ3soQk0/7zIyLZqwxZ+8Mxv9B6ug1pa9HOBFUW3nzGzL7j7C8B8YD25Vv43zawXGA3MAF6u\n4ZgViZt/funKTcw8L4O/tfuEFw9Q9YeCPkhE6qv4PXXKuJPYNfDeqTytIVGbWoLegH8run0zcL+Z\nHQF+D9zk7vvM7D7geXJDOW9z90M1HLMij//UI7fvGjjMivVbBm8Xv3iitsW9oLSQiUh9DX1PFYd8\nseVr39R7LIHEQe/udw25/WvgQxH7LQYWJz1OEgcOHa35OUq9oEotZKIXoUj14t5TQ2kNiWR0wVSM\nUi8oLWQiUl9x76mhtIZEMkEGfXdX5AijqpR6QWkhE5H6intPDaU1JJIJMujnzTyz5uco9YLSQiYi\n9RX3npqUHk13V4q+zDgWXXuBukYTCnI++k9ccR4Aq17cypFjx4HcVQRnZsZhZ5+Cv7XnhAVIoLpF\nSbSQiUh9Rb2nPnalMaOv7GU3UoGaroxtlP7+gUSF0gUWnUF17gyqc6LHR/ZbB9l1IyIi71HQi4gE\nTkEvIhI4Bb2ISOCCG3WjOWhEqqf3TdiCCvpVG7YEOQeN3oSNod9rjuZuCl9QQb90xeuR20fyHDR6\nEzbGSP+91vNDSnM3hS+ooH9re/T405E8B43ehI3Rit9rvcK53h9SmrspfEGdjD17SvTqKiN5Dhq9\nCRuj2b/XQjhv6T/A8Wx2MJzXbdxe9XOV+pBKQnM3hS+ooF84/9zI7SN5Dhq9CRuj2b/XuHBesvzV\nqsO+3h9SmrspfEEF/dyZfSy69gL6MuOCmQhJb8LGaPbvNS6cjxw7XnXLvt4fUrPPnxLc+0ZOFFQf\nPeRetCG9QDWBWmM0+/c69dQxbOmPb3FXc25gwZxpw1ZGy21P/iEV2vtGThRc0FdjpAyv05uwvob+\n3W+45vyG/37jwrmgmm4XffhLtTo26Ef68DpJplV/98JzL1n+6uDU2cWq7XbRh79UI3HQm9mvgX35\nm28A3wQeAbLAy8Dn3f24md0ILAKOAne4+7KaSlwnGrbYmVr5dy88f727XUTKSRT0ZtYLpNx9XtG2\np4Db3f3nZvZPwHVmtha4BbgI6AVWm9mz7h69xHsdFb6el+oXjbJ1x36++vC6tu/Okeqt27g99vWw\nbeeBpnTlqdtFWiHRwiNmNht4DHiT3IfFl4EngD53z5rZdcB/AJ4Brnb3z+Uf9yTw9+7+y1LPf/To\nsWxPT3fV5SpYtWELdz2+PvHjh/qbT85i7sy+uj2fNF+510TmlJPp3/POsO3628sIE7nwSNKum4PA\n3cBDwLnA0+Ra+IVPjQFgAjAe2Fv0uML2knbvrmxF+KEKq7P84JnfJHp8nB884227pJlW4alMudfE\nsYh+89zj2uNvr79zZ6jDClOR25MG/WvApnywv2ZmO4FZRfengT3k+vDTEdsbKm7McpzurhRnTB7L\n2/37ifp+005XoQ7tXvjYlX/YFkHU7kq9JhZdewGLf7wx8r52+tuLJJX0gqnPAPcAmNlUci33n5rZ\nvPz9VwHPAy8Al5pZr5lNAGaQO1HbUHEXlETpy4xj8a2X843PXsyZmeiRD12pFDfcuZKvPrwu0SXr\n9RJ1Gf1dj69vaZna2bqN2/nqw+u44c6VdMe80vsy45h9/hRdgSxBS9qifxh4xMxWkxtl8xlgB7DY\nzE4CXgV+5O7HzOw+cqHfBdzm7ofqUO6Syo1ZPnHfc4p+jn5cYThcq4dgaqRQvKHfdOzsiaxYv2Xw\n/uPHoh9X+Pvb2RMjT9RqNIyEIFHQu/u7wMcj7rosYt/FwOIkx2mkSeNHs3De9BMCcuiIiK5UKnLM\nc6uCVROcRYsaGx83umZUdxfHs9kTRrus27j9hA+Fgvmz+jr+A1TCEOQFU3EtX8h9Vf/GZy+Ovb/4\nQpQb7lwZuU+rgjXuMvpO714o9fce6ng2y+JbL6/o8f5Ww08niTRFkEH/domx89WEdLsFayPmOGkH\ntY5fr+bke9TfTt+UJHRBzV4J8L1nX4scOVNQTUi328yRUbMM/s0nZ43o7oV6zNNezcn3qL+dTsRK\n6IJq0a/asCWyr7VYNSHdjlcxDp3jZKSMNY5rtVdzgnnwOXYeZOrk956j0pPvcVPvhvpNSaQgqKCP\nWzO22INPvcLSlZtYePl0gGEjNfyt3cPCqDgcCkP2NEVC5eImEtv09t6SUxJU8hyQ+/BbunITuwbi\nZ9YoDKOM0o4f6CL1FFTQx60ZO9SugcPDWnBDR2pEDaXUjJfJLF25KXJ7qW9fQ7tNyrX89+x/t2QZ\nyrXONRukhCyoPvq4NWNrUbwOZ73X6uwUpVracQ4eOnJCP325E6Zx/eyjuru0WpJ0vKCCPm7N2FoU\ndyFodEbzFL51FcK+3AnTuBPnn1kwQyEvHS+ooJ87s4/uyLnbkivuQtDojBMVTzFQanqISenRiY9R\n+LZUbgSU1j0ViRdUHz3AsepnXS6pkikSOnF0RjXnKxZePr3iKSmGKnxbquSEqfrZRaIFF/SjulMc\nKZP2hekP4MTgsLNPwd/aUzJIhj6mU0dnxJ2vWLL8VYBhv7dyo2LiFH9bKgT5SBlSKtIuggv6uRee\nGTtvySeuOG/Y9mpDWq3GnLd3RJ+XOHLseGTLPmmrvhO/LYnUW3BB/4krzmPDa/3DWo8r1g+/mKov\nM3aw77fRS8hFqdfSdc1YAq/4WEtXbqLcwmQPLdvI0pWb2LP/3cFrFKLMn9U3+C1qwriTIAt7D7zb\n0d+WROot0VKCjdbfP5CoUJlMmpv/+79UvU5slEafyBvax530uK9u2Ru5RF6t5f/es6+x6sW3OXIs\ny6juFHMvPJPpZ05I3Ncep9wkc1E6setmpNa5lkbISK1zLeqwwlRdlxJsS6s2bKlLyEPjpyKu19zy\ncVcDl3ueqDdgoVxDf4dHjmVZsX4LP/t16eklktDQ1HDpAsP2EVTQVzIFQqVqDaByLZl6jcmPuxo4\n6nkKZXp7x4ETul6K34ClNOLLX6cOTe0EWiinfQQ1jn7ztn11e65aAqiSGRnrNSY/7mrgoc9TXKZ2\n6q3TydZw6QLD9hFU0NeTnX1K4sdWMlVCvaZAjrsaeOjzVLM4R6PMn9WnC5o6iC4wbB9Bdd3UUy2r\nC8W1ZLbu2H/CzJcX/MEkXntrD0eOHWdUdxdzL5xaNviGdgldaKdF7vedH7/C8rWbB7uMSi3GUk/d\nqRTHyTIxPVojaDqcLjBsH4mC3sxGAUuAacBo4A7gd8AyoNBR/oC7/9DMbgQWAUeBO9x9Wa2FboZa\nvl7GrUx1PMvg9qGzZR45dpwV67cw/cwJsYEYvTbqG5H7ZrPv9b2vfmlbycVYatGXGacQl0i6wLB9\nJG3RfxLY6e6fMrNJwIvAN4B73f2ewk5mdjpwC3AR0AusNrNn3b36SySbrJavl5UuhBEl7kTVuo3b\nWbJ8Y6LnfOWNXYkeVwm9caUUXWDYHpIG/VLgR/mfU+Ra67MAM7PryLXqvwRcDKzJB/thM9sEvA/4\nZU2lboLCNLnVDlEsfmEXt2Te7t9fUas66pvE9559rezKWa2iERQi7S9R0Lv7fgAzS5ML/NvJdeE8\n5O7rzew24GvkWvp7ix46AEwo9/wTJ46hp6e7qjKt2lDfICxMk/v2roMs+o/vizxe1Bjh8eN7mTuz\nj2suS3PNZdMH7//C3SsrGhV01pQ0mUx68Bj/e9lGdux5pw41aoxtOw8MlreZWnHMVlOdO0Mj6pz4\nZKyZnQU8CXzb3b9vZqe4e+EM5pPA/cAqoLjUaaDsWc7du6NPZpby8P97uerHVGLZ6jc4c9KYYa3W\nHzzzm8j9f/CMM6Nv+GfZlX90VkXdOVf+0Vkse25T4knAGmFUdxfpMaMiy3PG5LFNv3pRV0x2BtU5\n2eOjJBpeaWZTgJ8Cf+vuS/KbnzGzwrXs84H1wAvApWbWa2YTgBlAQxK5kaEYtYJUtWOEZ58/peS8\n7KO6u5g/q4/VL23jwadeaZuQh9ziHYU1dofSCAqR9pe0Rf9lYCLwFTP7Sn7bXwH/aGZHgN8DN7n7\nPjO7D3ie3IfKbe5+qNZCN9uW/v3DtsWNrCl1ErfUuqaFUTftYFJ6dOywSI2gEBl5kvbRfxH4YsRd\nH4rYdzGwOMlx2lmpMcJxJ2njPhzayaTxo7n7L4b9GQGNoBAZqXTBVIU+e+fP6OlKcfRYljG9PRw8\nfPSE+8ee3MPonm6+8+NXYueRqWXYZbMUFmQRkXBoCoQKZbO5WRyzwIFDR4fNF3PgnaPsGjgcO4/M\n0p9vAiBV5zVtazX6pG5NSSASOLXom2TXvsMtb80XL/KhPnaRzqGg7wCT0qNZePl0hbpIh1LQB6q7\nK6VWu4gACvogxS2ELiKdKYigL17Qo5OlUvDv/51CXkROFETQt8OiGq3U3ZVi3swzFfAiEimIoI+b\njiBUxXPAd+J8ICJSnSCC/pRxJ7XV3DCNpLHuIlItXTA1gsyf1aeQF5GqBdGi75TWfC3r2FYqbp4e\nERm5ggj6TlHLOraViFqTtnBbYS8ycqnrpk10p1Jc8AeT6MuMi92nlnVsKxE3eilqPn4RGTnUom+x\nvsw4vvHZiwdvD21VF2v0Ih/VLqYiIiODWvQtNjRE41rVk8aPbnj3ydRTx0Rub/Q3CRFpLLXoW2xo\niMa1qveWWJ2qoNYTqaUWUxGRkUtB32JDQzTJEoVQnxOphf20XKBIWBT0dTYpPbri4Z5RFz8lbVWX\nOpFaTVBruUCR8DQ86M2sC/g28H7gMHCDu29q9HFbIZWCMb09FQV9X2ZcZKAmbVXrRKqIxGlGi/5P\ngV53n2NmHwDuAa5rwnGb7sxTx7F1R2XBWqqFnqRVnbTLR0TC14xRN5cAPwFw918AFzXhmC2xYM45\nsSNXRnV3NXRt1gVzpsWWSUQ6WzNa9OOBvUW3j5lZj7sfjXvAxIlj6OnpbnzJ6qgrBddcNp3x43u5\n6/H1w+7/0sdmMndmX03HWLVhC0tXvM5b2wc4e0qahfPPJZNJk8mkueayNOPH97J0xev8bvsAZ+Xv\nr/WY7SqTSbe6CE2nOneGRtS5GUG/DygueVepkAfYvXvkTTs89dRx9PcPMKNvAouuvWBYH/uMvgk1\nTSc8dFTN5m37Bj9QZvRNGPz/q9ef+IUpxCmMO3FqZtW5M9Ra57gPiWYE/RrgT4D/m++jf6kJx2y6\n4i6SRoxciRtVs3TF68PCXUSkWDOC/kngCjP7VyAF/HkTjtk0k8aPZuG86Q0fkhg3quZ32zurxSMi\n1Wt40Lv7ceBzjT5Oq9z9Fx9qynHiRtWcNaXz+jBFpDqa66YGk9Kjm3asuFE1C+ef27QyiMjIpCtj\na7Dw8ulNO1bchVRzZ/Z13AkrEamOgr5KXancCJtWzAGj6QlEJAkFfRW0MLeIjEQK+hK6UymyZFvW\nghcRqQcFPbnJyLpSKU3LKyJBCiLox508iv3vHEn8+Jv+RF0yIhKuIIZXfuKK8yrab2xvD6lUboKx\nrhQNm2BMRKSdBNGiLwT19//lNQYODm/Zj+ruYu6FUyv+QBARCUkQQQ+5sL/msukaUy4iMkQQXTci\nIhJPQS8iEjgFvYhI4BT0IiKBU9CLiAQulc1mW10GERFpILXoRUQCp6AXEQmcgl5EJHAKehGRwCno\nRUQCp6AXEQmcgl5EJHAjfvZKM+sCvg28HzgM3ODum1pbqvoxs1HAEmAaMBq4A9gIPAJkgZeBz7v7\ncTO7EVgEHAXucPdlrShzvZjZacB64ApydXqEgOtsZv8NuBY4idxr+jkCrnP+tf0oudf2MeBGAv47\nm9ls4E53n2dm06mwnmZ2MvA4cBowAFzv7v3VHDuEFv2fAr3uPgf4r8A9LS5PvX0S2OnulwJ/DPxP\n4F7g9vy2FHCdmZ0O3AJ8CLgS+AczG92iMtcsHwIPAu/kNwVdZzObB3yQXF0uA84i8DoDVwM97v5B\n4BvANwm0zmZ2K/AQ0JvfVE09bwZeyu/7GHB7tccPIegvAX4C4O6/AC5qbXHqbinwlfzPKXKf9LPI\ntfYAngY+DFwMrHH3w+6+F9gEvK/JZa2nu4F/Arbmb4de5yuBl4AngR8Dywi/zq8BPflv5eOBI4Rb\n598CHym6XU09BzOuaN+qhBD044G9RbePmdmI75IqcPf97j5gZmngR+Q+zVPuXpi7YgCYwPDfQ2H7\niGNmnwb63f2Zos1B1xk4lVwjZSHwOeB7QFfgdd5PrtvmN8Bi4D4C/Tu7+xPkPsgKqqln8fZEdQ8h\n6PcB6aLbXe5+tFWFaQQzOwtYCXzX3b8PHC+6Ow3sYfjvobB9JPoMcIWZ/Ry4kNzX1dOK7g+xzjuB\nZ9z9XXd34BAnvqFDrPNfkqvzeeTOsT1K7vxEQYh1LqjmPVy8PVHdQwj6NeT6+jCzD5D7+hsMM5sC\n/BT4W3dfkt+8Id+nC3AV8DzwAnCpmfWa2QRgBrmTPCOOu89198vcfR7wIvBnwNMh1xlYDfyxmaXM\nbCowFlgReJ13815LdRcwisBf20WqqedgxhXtW5UQujieJNf6+1dyfdh/3uLy1NuXgYnAV8ys0Ff/\nReA+MzsJeBX4kbsfM7P7yL0IuoDb3P1QS0rcGH8NLA61zvnRFXPJvdm7gM8DbxBwnYF/BJaY2fPk\nWvJfBn5F2HUuqPj1bGYPAI+a2WrgXeDj1R5M0xSLiAQuhK4bEREpQUEvIhI4Bb2ISOAU9CIigVPQ\ni4gETkEvIhI4Bb2ISOD+P+4fjHENwSnjAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x117e07940>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(cat['F_PACS_160'],old_cat['F_PACS_160'][idx],'o' )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [conda env:new]",
   "language": "python",
   "name": "conda-env-new-py"
  },
  "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.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}