dmu7/1_Make_sxds_CIGALE_input.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "5ecae634-7d42-46b1-9ca9-e6494bc7de93",
   "metadata": {},
   "source": [
    "# Make CIGALE test input from VIDEO SXDS data\n",
    "\n",
    "This notebook uses the HELP python environment to make a CIGALE input catalogue:\n",
    "\n",
    "https://github.com/H-E-L-P/herschelhelp_internal\n",
    "\n",
    "to make the CIGALE inputs we must correct for galactic extinction, convert flux measurements from nano-Janksy to milli-Jansky, merge in the spectroscopic redshifts.\n",
    "\n",
    "This is an early test on the catalogues as part of our validation process."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "cb0728b0-8431-426f-bb42-9860f1c7421e",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/_collections_abc.py:841: MatplotlibDeprecationWarning: \n",
      "The examples.directory rcparam was deprecated in Matplotlib 3.0 and will be removed in 3.2. In the future, examples will be found relative to the 'datapath' directory.\n",
      "  self[key] = other[key]\n",
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/_collections_abc.py:841: MatplotlibDeprecationWarning: \n",
      "The savefig.frameon rcparam was deprecated in Matplotlib 3.1 and will be removed in 3.3.\n",
      "  self[key] = other[key]\n",
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/_collections_abc.py:841: MatplotlibDeprecationWarning: \n",
      "The text.latex.unicode rcparam was deprecated in Matplotlib 3.0 and will be removed in 3.2.\n",
      "  self[key] = other[key]\n",
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/_collections_abc.py:841: MatplotlibDeprecationWarning: \n",
      "The verbose.fileo rcparam was deprecated in Matplotlib 3.1 and will be removed in 3.3.\n",
      "  self[key] = other[key]\n",
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/_collections_abc.py:841: MatplotlibDeprecationWarning: \n",
      "The verbose.level rcparam was deprecated in Matplotlib 3.1 and will be removed in 3.3.\n",
      "  self[key] = other[key]\n",
      "/Users/raphaelshirley/miniconda3/envs/herschelhelp_internal/lib/python3.6/site-packages/seaborn/apionly.py:9: UserWarning: As seaborn no longer sets a default style on import, the seaborn.apionly module is deprecated. It will be removed in a future version.\n",
      "  warnings.warn(msg, UserWarning)\n"
     ]
    }
   ],
   "source": [
    "\n",
    "%matplotlib inline\n",
    "#%config InlineBackend.figure_format = 'svg'\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "plt.rc('figure', figsize=(10, 6))\n",
    "from herschelhelp.filters import correct_galactic_extinction\n",
    "from herschelhelp.external import convert_table_for_cigale\n",
    "from herschelhelp_internal.masterlist import merge_catalogues, nb_merge_dist_plot, specz_merge\n",
    "from herschelhelp_internal.utils import coords_to_hpidx, ebv, gen_help_id, inMoc\n",
    "\n",
    "import numpy as np\n",
    "import os\n",
    "import time\n",
    "\n",
    "from astropy.table import Table\n",
    "import astropy.units as u\n",
    "from astropy.coordinates import SkyCoord"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "8d558663-1054-4a71-9960-e60efa5ac8c8",
   "metadata": {},
   "outputs": [],
   "source": [
    "SUFFIX = os.environ.get('SUFFIX', time.strftime(\"%Y%m%d\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ef277a05-970c-479d-a2fa-a8f10210851c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'20210906'"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "SUFFIX"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "0e8fe7f2-d70a-4c6b-a04d-477dab414fd2",
   "metadata": {},
   "outputs": [],
   "source": [
    "hsc_bands = 'GRIZY'\n",
    "vista_bands = 'ZYJHK'\n",
    "bands = ['HSC_{}'.format(b) for b in hsc_bands]\n",
    "bands += ['VISTA_{}'.format(b) for b in vista_bands]\n",
    "cols = ['id', 'VISTA_Ks_m_coord_ra','VISTA_Ks_m_coord_dec','VISTA_Ks_m_detect_isPrimary']\n",
    "cols+=['{}_m_slot_ModelFlux_flux'.format(b.replace('K','Ks')) for b in bands]\n",
    "cols+=['{}_m_slot_ModelFlux_fluxErr'.format(b.replace('K','Ks'))  for b in bands]\n",
    "cat = Table.read('../dmu5/dmu5_VIDEO/data/full_reduced_cat_SXDS.fits')[cols]\n",
    "cat['VISTA_Ks_m_coord_ra'].name='ra'\n",
    "cat['VISTA_Ks_m_coord_dec'].name='dec'\n",
    "#Some HELP functions require a help_id col\n",
    "cat['id'].name='help_id'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "e856cd37-caa4-4520-b708-66447b3e4ada",
   "metadata": {},
   "outputs": [],
   "source": [
    "#Rename to CIGALE/HELP format\n",
    "for c in cat.colnames:\n",
    "    if c.endswith('m_slot_ModelFlux_flux'):\n",
    "        mask = cat[c]>1e10\n",
    "        cat[c][mask]=np.nan\n",
    "        f_name='f_{}_{}'.format(\n",
    "            c.split('_')[0].lower().replace('hsc','suprime'),\n",
    "            c.split('_')[1].lower()\n",
    "        )\n",
    "        cat[c].name=f_name\n",
    "        cat[f_name]=cat[f_name].to(u.uJy)\n",
    "    if c.endswith('m_slot_ModelFlux_fluxErr'):\n",
    "        mask = cat[c]>1e10\n",
    "        cat[c][mask]=np.nan\n",
    "        ferr_name='ferr_{}_{}'.format(\n",
    "            c.split('_')[0].lower().replace('hsc','suprime'),\n",
    "            c.split('_')[1].lower()\n",
    "        )\n",
    "        cat[c].name=ferr_name\n",
    "        cat[ferr_name]=cat[ferr_name].to(u.uJy)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "702240bb-d351-4d7d-af27-d97bd7b6b108",
   "metadata": {},
   "outputs": [],
   "source": [
    "mask = (cat['ra']<1e10) | (cat['dec']<1e10)\n",
    "mask&=cat['VISTA_Ks_m_detect_isPrimary']==1\n",
    "mask&=(cat['f_suprime_i']/cat['ferr_suprime_i'] > 5)\n",
    "mask&=(cat['f_vista_ks']/cat['ferr_vista_ks'] > 5)\n",
    "cat = cat[mask]\n",
    "    \n",
    "cat['ra']=cat['ra'].to(u.deg)\n",
    "cat['dec']=cat['dec'].to(u.deg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "6753b192-a84a-4aee-ab17-e31ff31e6cb6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=5</i>\n",
       "<table id=\"table140196264302632\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>ra</th><th>dec</th><th>VISTA_Ks_m_detect_isPrimary</th><th>f_suprime_g</th><th>f_suprime_r</th><th>f_suprime_i</th><th>f_suprime_z</th><th>f_suprime_y</th><th>f_vista_z</th><th>f_vista_y</th><th>f_vista_j</th><th>f_vista_h</th><th>f_vista_ks</th><th>ferr_suprime_g</th><th>ferr_suprime_r</th><th>ferr_suprime_i</th><th>ferr_suprime_z</th><th>ferr_suprime_y</th><th>ferr_vista_z</th><th>ferr_vista_y</th><th>ferr_vista_j</th><th>ferr_vista_h</th><th>ferr_vista_ks</th></tr></thead>\n",
       "<thead><tr><th></th><th>deg</th><th>deg</th><th></th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>float64</th><th>float64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>149940984795367161</td><td>34.03198792407938</td><td>-5.636746033451176</td><td>1</td><td>0.1918380120633462</td><td>0.24114609138728316</td><td>0.23910295624418496</td><td>0.3943670650640425</td><td>0.539275444256114</td><td>3.3093421221815893</td><td>0.5509247061854611</td><td>0.41870752165663716</td><td>0.3170050549873974</td><td>1.4210650305132755</td><td>0.01933043244302869</td><td>0.030078245621829574</td><td>0.027566898414982985</td><td>0.0778258370449352</td><td>0.1583567806404285</td><td>0.06724782055241774</td><td>0.024002575829053612</td><td>0.023418802960395863</td><td>0.035446290741617044</td><td>0.06803846082798298</td></tr>\n",
       "<tr><td>149940984795367228</td><td>34.03239043843923</td><td>-5.63513409597789</td><td>1</td><td>0.5193620245085407</td><td>0.8073257364329538</td><td>0.6775046773308321</td><td>1.1696392994439777</td><td>1.376104949778619</td><td>0.7212169913401356</td><td>1.1436669705612827</td><td>1.5161937043715057</td><td>0.7363723324967524</td><td>1.9766101985524038</td><td>0.022960893798619434</td><td>0.05553958305818498</td><td>0.03954325169680315</td><td>0.0898612977774623</td><td>0.19791099115015662</td><td>0.00928557764376431</td><td>0.018521803949740616</td><td>0.029229104310411657</td><td>0.03511634243528409</td><td>0.054673791964521476</td></tr>\n",
       "<tr><td>149940984795367270</td><td>34.092669339372655</td><td>-5.634361990777769</td><td>1</td><td>0.36384861296157717</td><td>0.7892301676025678</td><td>1.2315274181560765</td><td>1.4359311773782693</td><td>1.8894775961718346</td><td>1.0766499332264072</td><td>1.0725330148649175</td><td>1.2880861811268447</td><td>2.062232090325329</td><td>1.6832976225622505</td><td>0.017392853109826446</td><td>0.02630454007725901</td><td>0.02704338733473185</td><td>0.05569826571177505</td><td>0.19013509621684996</td><td>0.0107364834608036</td><td>0.018037655838292282</td><td>0.029016609691475705</td><td>0.048036069653102864</td><td>0.06734791599744612</td></tr>\n",
       "<tr><td>149940984795367288</td><td>34.107736780652246</td><td>-5.634215412521896</td><td>1</td><td>0.1386149795474109</td><td>0.05213683841882216</td><td>0.17740765789056012</td><td>0.3595270423291327</td><td>0.8112687409022378</td><td>0.2476538292474139</td><td>0.31713807754942025</td><td>0.8667135068230658</td><td>1.2445507718205489</td><td>0.9546880507128258</td><td>0.022092521907655677</td><td>0.026488851509346475</td><td>0.029234692678244106</td><td>0.057606192085304514</td><td>0.15846595810001673</td><td>0.010854198644272166</td><td>0.018258807749051895</td><td>0.03785618111103003</td><td>0.11295631819774835</td><td>0.08459614884109998</td></tr>\n",
       "<tr><td>149940984795367292</td><td>33.96461795640932</td><td>-5.633879365508137</td><td>1</td><td>0.5398559354684481</td><td>0.713206483879246</td><td>0.7658640563342345</td><td>1.1315693305752763</td><td>1.535094461866041</td><td>0.5064604919459583</td><td>0.6504969996625733</td><td>1.0177887608324134</td><td>nan</td><td>1.034369760515549</td><td>0.019511755613238497</td><td>0.0456005272285144</td><td>0.030271430955299255</td><td>0.06703133567861314</td><td>0.20062454166051713</td><td>0.010458156541137244</td><td>0.01701051807225166</td><td>0.026377929025613298</td><td>nan</td><td>0.06331285964531365</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=5>\n",
       "     help_id               ra         ...    ferr_vista_ks    \n",
       "                          deg         ...         uJy         \n",
       "      int64             float64       ...       float64       \n",
       "------------------ ------------------ ... --------------------\n",
       "149940984795367161  34.03198792407938 ...  0.06803846082798298\n",
       "149940984795367228  34.03239043843923 ... 0.054673791964521476\n",
       "149940984795367270 34.092669339372655 ...  0.06734791599744612\n",
       "149940984795367288 34.107736780652246 ...  0.08459614884109998\n",
       "149940984795367292  33.96461795640932 ...  0.06331285964531365"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cat[0:5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "82f1ed4d-940a-4e94-844f-98622dfe3a39",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mkdir: data: File exists\n",
      "--2021-09-06 13:38:25--  http://hedam.lam.fr/HELP/dataproducts/dmu23/dmu23_XMM-LSS/data/XMM-LSS-specz-v2.91.fits.gz\n",
      "Resolving hedam.lam.fr (hedam.lam.fr)... 195.221.212.141\n",
      "Connecting to hedam.lam.fr (hedam.lam.fr)|195.221.212.141|:80... connected.\n",
      "HTTP request sent, awaiting response... 200 OK\n",
      "Length: 3391102 (3.2M) [application/x-gzip]\n",
      "Saving to: ‘./data/XMM-LSS-specz-v2.91.fits.gz.9’\n",
      "\n",
      "XMM-LSS-specz-v2.91 100%[===================>]   3.23M  5.61MB/s    in 0.6s    \n",
      "\n",
      "2021-09-06 13:38:26 (5.61 MB/s) - ‘./data/XMM-LSS-specz-v2.91.fits.gz.9’ saved [3391102/3391102]\n",
      "\n"
     ]
    }
   ],
   "source": [
    "!mkdir data\n",
    "!wget http://hedam.lam.fr/HELP/dataproducts/dmu23/dmu23_XMM-LSS/data/XMM-LSS-specz-v2.91.fits.gz -P ./data/"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "08c6e594-4a87-4371-bbd1-63c1973dc41a",
   "metadata": {},
   "source": [
    "## Merge Specz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "19b84aad-3093-42f2-93bc-bc64dc209f6c",
   "metadata": {},
   "outputs": [],
   "source": [
    "specz=Table.read('./data/XMM-LSS-specz-v2.91.fits.gz')\n",
    "for c in specz.colnames:\n",
    "    specz[c].name=c.lower()\n",
    "specz['objid'].name='specz_id'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "02467758-646e-41ae-a556-e2c499332892",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "127682"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(specz)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "3cee166e-600e-47b3-9bd3-1f8b503ad642",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAFzCAYAAAAaBYxDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXicZ3nv8d892kYz2ndZkiXZlrc4cRbH2TdIICyHABdLEggFkqahgQJd6XV6yulyDrQ9h1MKgTQNEGgJAZq0pJAQWpIQyOLYcWwnXiVbsiTL1r7v0jznjxkJoXgZe0Z6Z6Tv57rmkmbmnXdujxPp5+e53+cx55wAAABwbnxeFwAAAJDMCFMAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMUj16o2LiopcTU2NV28PAAAQtVdeeaXLOVd8suc8C1M1NTXasWOHV28PAAAQNTM7eqrnmOYDAACIAWEKAAAgBoQpAACAGBCmAAAAYkCYAgAAiAFhCgAAIAaEKQAAgBicMUyZ2TfNrMPMXj/F8x8ysz2R2wtmtjn+ZQIAACSmaEamHpJ082meb5R0nXPuAkl/JemBONQFAACQFM64Arpz7jkzqznN8y/MufuSpMrYywIAAEgO8e6ZulPSk6d60szuNrMdZrajs7Mzzm8NAACw+OIWpszsBoXD1J+c6hjn3APOuS3OuS3FxSfdKxAAACCpxGWjYzO7QNKDkt7mnOuOxzkBAACSQcxhysxWSnpM0h3OuUOxl4RYPbytOarjbr9s5QJXAgDA0nfGMGVm35N0vaQiM2uV9HlJaZLknLtf0p9LKpT0NTOTpCnn3JaFKhgAACCRRHM1321neP4uSXfFrSIAAIAkwgroAAAAMSBMAQAAxIAwBQAAEAPCFAAAQAwIUwAAADEgTAEAAMSAMAUAABADwhQAAEAMCFMAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMSBMAQAAxIAwBQAAEAPCFAAAQAwIUwAAADEgTAEAAMSAMAUAABADwhQAAEAMCFMAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMSBMAQAAxIAwBQAAEAPCFAAAQAwIUwAAADEgTAEAAMSAMAUAABADwhQAAEAMCFMAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMSBMAQAAxIAwBQAAEAPCFAAAQAwIUwAAADEgTAEAAMSAMAUAABCDM4YpM/ummXWY2euneN7M7B/MrMHM9pjZxfEvEwAAIDFFMzL1kKSbT/P82yTVRW53S/p67GUBAAAkhzOGKefcc5J6TnPILZK+48JekpRnZuXxKhAAACCRxaNnqkJSy5z7rZHH3sDM7jazHWa2o7OzMw5vDQAA4K14hCk7yWPuZAc65x5wzm1xzm0pLi6Ow1sDAAB4Kx5hqlVS1Zz7lZLa4nBeAACAhBePMPW4pI9Eruq7XFK/c+54HM4LAACQ8FLPdICZfU/S9ZKKzKxV0uclpUmSc+5+SU9IerukBkkjkj62UMUCAAAkmjOGKefcbWd43km6N24VAQAAJBFWQAcAAIgBYQoAACAGhCkAAIAYEKYAAABiQJgCAACIAWEKAAAgBoQpAACAGBCmAAAAYkCYAgAAiAFhCgAAIAaEKQAAgBgQpgAAAGJAmAIAAIgBYQoAACAGhCkAAIAYEKYAAABiQJgCAACIAWEKAAAgBoQpAACAGBCmAAAAYkCYAgAAiAFhCgAAIAaEKQAAgBgQpgAAAGJAmAIAAIgBYQoAACAGhCkAAIAYEKYAAABiQJgCAACIAWEKAAAgBoQpAACAGBCmAAAAYkCYAgAAiAFhCgAAIAaEKQAAgBgQpgAAAGKQ6nUBiN7D25q9LgEAAMzDyBQAAEAMCFMAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMYgqTJnZzWZ20MwazOxzJ3k+18z+w8x2m9leM/tY/EsFAABIPGcMU2aWIuk+SW+TtFHSbWa2cd5h90ra55zbLOl6Sf/XzNLjXCsAAEDCiWZkaqukBufcEefchKRHJN0y7xgnKdvMTFKWpB5JU3GtFAAAIAFFE6YqJLXMud8aeWyur0raIKlN0muSPu2cC8WlQgAAgAQWTZiykzzm5t1/q6RdklZIulDSV80s5w0nMrvbzHaY2Y7Ozs6zLhYAACDRRBOmWiVVzblfqfAI1Fwfk/SYC2uQ1Chp/fwTOececM5tcc5tKS4uPteaAQAAEkY0YWq7pDozq400ld8q6fF5xzRLerMkmVmppHWSjsSzUAAAgESUeqYDnHNTZvZJSU9JSpH0TefcXjO7J/L8/ZL+StJDZvaawtOCf+Kc61rAugEAABLCGcOUJDnnnpD0xLzH7p/zfZukt8S3NAAAgMTHCugAAAAxIEwBAADEgDAFAAAQA8IUAABADAhTAAAAMSBMAQAAxIAwBQAAEAPCFAAAQAwIUwAAADEgTAEAAMSAMAUAABADwhQAAEAMCFMAAAAxIEwBAADEINXrAhA/nYPjOtY3qvaBMXUMjGlofEq3bV2pvEC616UBALBkEaaWiPqOQT30fJOcJJ9JRVkZ6hoa1y8buvTfLljhdXkAACxZhKkl4qUjPQpkpOrOq2pVlJ2uVJ9PP9zRoleaevXm9SUKpPNXDQDAQqBnagkYGJvUwRMDumRlnspy/Ur1hf9ar64r0sR0SC839nhcIQAASxdhagl4tblPISdtqS74jcfLczNVV5KlFw93a2o65FF1AAAsbYSpJOec046mHtUUBlSUnfGG56+uK9Lg+JR2tfR5UB0AAEsfYSrJNXYPq3t4QltqCk76/JriLJXn+vXLhi6FnFvk6gAAWPoIU0luR1OvMlJ92rQi96TPm5muXlOkzsFx1bcPLnJ1AAAsfYSpJDY6Ma3Xj/Xrwqo8paee+q/ygso85Wam6Zf1XYtYHQAAywNhKontbu3TVMi9ofF8vhSf6crVhTrSNaz2gbFFqg4AgOWBMJXEdjT1qDzXrxV5/jMee35FeBrwcOfQQpcFAMCyQphKUsf7R9XWP6YtNQUyszMenxdIV15mmpq6hhehOgAAlg/CVJJqjISijeU5Ub+mpiiopu4ROa7qAwAgbghTSaq5Z0Q5/lTlZqZF/ZrqwoCGxqfUPTyxgJUBALC8EKaSVGvvqKoKAmf1mprCoCQx1QcAQBwRppLQ0PiUeoYnVJV/dmGqJDtDgfQUNXWPLFBlAAAsP4SpJNTaEw5DZzsyZWaqKQyqqZuRKQAA4oUwlYSae0fkM6kiL/OsX1tTGFDP8IQGxiYXoDIAAJYfwlQSau0ZVVmO/7Srnp9KTRF9UwAAxBNhKsmEnFNL74gqz3KKb0Z5bqbSU3xM9QEAECeEqSTTOTiu8amQVp5l8/mMFJ9pZUFATV00oQMAEA+EqSTT2hsOQZUFZ98vNaO6KKD2gTH1j9I3BQBArAhTSaa5Z1T+NJ+KsjLO+Rw1hUE5Sa8c7YlfYQAALFOEqSTT2juiqvyAfFHsx3cqVfkBpZjp5cbeOFYGAMDyRJhKIuNT0zrRP6bKc+yXmpGe6lNFfqa2NzEyBQBArAhTSeRY36icpJUx9EvNqCkMaE9rn8Ymp2MvDACAZYwwlURaekYlKeaRKUmqLgxqctrptWP9MZ8LAIDljDCVRFp6RlQYTFcwIzXmc62IrJ5+4PhAzOcCAGA5I0wlCRdZrPNs9+M7lRx/qnIz07Tv+GBczgcAwHIVVZgys5vN7KCZNZjZ505xzPVmtsvM9prZL+JbJo73j2lwbEpV+bH3S0nhTY/Xl2XrwAlGpgAAiMUZw5SZpUi6T9LbJG2UdJuZbZx3TJ6kr0l6l3PuPEnvX4Bal7XXI71N8eiXmrGhPEcHTwwqFHJxOycAAMtNNCNTWyU1OOeOOOcmJD0i6ZZ5x9wu6THnXLMkOec64lsm6juGJEkl2ee+WOd8G8qzNTIxreYetpYBAOBcRROmKiS1zLnfGnlsrrWS8s3sWTN7xcw+Eq8CEVbfPqi8zDRlpKXE7Zzry3Ikiak+AABiEE2YOtlS2/PnhVIlXSLpHZLeKul/mNnaN5zI7G4z22FmOzo7O8+62OXsUPuQSnLiNyolSWtLs+Uz0YQOAEAMoglTrZKq5tyvlNR2kmN+6pwbds51SXpO0ub5J3LOPeCc2+Kc21JcXHyuNS870yGnw51DKsn2x/W8mekpqikKsjwCAAAxiCZMbZdUZ2a1ZpYu6VZJj8875keSrjGzVDMLSLpM0v74lrp8tfSMaHwqFNd+qRkbynN04AQjUwAAnKszhinn3JSkT0p6SuGA9APn3F4zu8fM7okcs1/STyXtkfSypAedc68vXNnLy0zzeWlOfEemJGlDWbaae0Y0ODYZ93MDALAcRLWUtnPuCUlPzHvs/nn3/07S38WvNMw41B4eOSpegJGpmSb0Q+2DuqS6IO7nBwBgqWMF9CTQ0DGkFbl++eN4Jd+MDSvCYYomdAAAzg1hKgnUdwxqTWn2gpx7Ra5f2f5UmtABADhHhKkEFwo5NXQMaW1J1oKc38y0oYwmdAAAzhVhKsG19o5qbDKkutKFCVNSeCX0A8cH2FYGAIBzQJhKcDPN52tKFmaaT5LWl+doeGJarb2jC/YeAAAsVYSpBDezLMJCjkytLwsHtf1sKwMAwFkjTCW4+vZBleX4leNPW7D3WFeWLTNpP03oAACcNcJUgqvvGFrQUSlJCqSnqqYwqAMsjwAAwFkjTCWwmSv56hawX2rG+rJspvkAADgHhKkEdqxvVKOT01q7wCNTUniPvqPdIxoen1rw9wIAYCkhTCWw+o7wtNtCT/NJv25CP9jOVB8AAGeDMJXADrWHr+RbyGURZqyNrLDeEHlPAAAQHcJUAqtvH1JpToZyMxfuSr4ZVQUBZaT6Zte1AgAA0SFMJbD6jsFFaT6XpBSfaXVx1uy6VgAAIDqEqQQ1eyXfIvRLzagrzVI9I1MAAJwVwlSCOtY3qpGJ6UUbmZLCfVNt/WMaHJtctPcEACDZpXpdAKSHtzW/4bGDJ8IjRE1dwyd9fiGsKQmPgh3uHNaFVXmL8p4AACQ7RqYSVNfQuCSpKDtj0d6zLhKmaEIHACB6hKkE1TU0rsy0FAXTUxbtPVcWBJSe6lMDTegAAESNMJWgOofGVZSVLjNbtPdMTfFpVVGQkSkAAM4CYSpBdQ2Oqyhr8ab4ZqwtzVY9C3cCABA1wlQCGp+a1sDYlIoXsV9qRl1Jlo71jbJHHwAAUSJMJaDuoQlJUqEHI1N1M9vK0DcFAEBUCFMJaOZKvmJPwlT4ij5WQgcAIDqEqQTUOTQuk1SYlb7o711dEFB6io+V0AEAiBJhKgF1D00oN5CmtJTF/+tJTfFpVXGQkSkAAKJEmEpAnR5dyTdjTUkWyyMAABAlwlSCcc6pa8jbMLW2NFutvaMameCKPgAAzoQwlWCGxqc0PhVSsQf9UjNmtpU53DHsWQ0AACQLwlSC6ZzZk8/DkamZK/qY6gMA4MwIUwmmezC8xtRibnA8X3VhUGkpRhM6AABRIEwlmM6hcaX6TLmZaZ7VkJbiU21RkOURAACIAmEqwcw0n/sWcYPjk6krzWZkCgCAKBCmEkzX0Lgni3XOV1eSpZbeEY1OTHtdCgAACY0wlUCmQ049wxOebCMz39rSbDknHe5kdAoAgNMhTCWQ3uEJhZy3zeczZpZHqO+gbwoAgNMhTCWQRFgWYUZ1YVCpPtOhdkamAAA4HcJUAumaDVPe90ylp85c0UeYAgDgdAhTCaRraFyB9BQF0lO9LkVSePFOpvkAADg9wlQC6RpKjObzGXUl2WruGdHYJFf0AQBwKoSpBNI16O0Gx/PVlWbJOamB9aYAADglwlSCGJuc1uD4VEJcyTdjbWm2JMIUAACnQ5hKEDPN58UJ0Hw+o6YwqBSf0TcFAMBpEKYSxEyYKkygab70VJ9qCgMsjwAAwGlEFabM7GYzO2hmDWb2udMcd6mZTZvZ++JX4vLQNTQhk1QYTJyRKSk81cc0HwAAp3bGa/DNLEXSfZJuktQqabuZPe6c23eS4/5G0lMLUehS1zk4rvxgulJTFm+w8OFtzWc8ZmwypKPdwxqbnJY/LWURqgIAILlE85t7q6QG59wR59yEpEck3XKS4z4l6VFJHXGsb9noHhpPiMU65yvNyVDISUc6h70uBQCAhBRNmKqQ1DLnfmvksVlmViHpPZLuP92JzOxuM9thZjs6OzvPttYlyzmnrqGJhFoWYUZJtl8Se/QBAHAq0YQpO8ljbt79v5f0J865067u6Jx7wDm3xTm3pbi4ONoal7yBsSlNTIcSMkwVZaWHr+ijCR0AgJOKZt+SVklVc+5XSmqbd8wWSY+YmSQVSXq7mU055/49LlUucV0JtMHxfKkpPlUXBhiZAgDgFKIJU9sl1ZlZraRjkm6VdPvcA5xztTPfm9lDkn5MkIpe52BkjakEWrBzrrqSLEamAAA4hTNO8znnpiR9UuGr9PZL+oFzbq+Z3WNm9yx0gctB99C40lJMOf7E2OB4vrWl2WrqHtb4FHv0AQAwX1S/vZ1zT0h6Yt5jJ202d859NPaylpfOofCefJFp0oSzpiRr9oq+DeU5XpcDAEBCYQX0BJCoV/LNmNmjr57FOwEAeAPClMfGp6bVO5zYYaq2KCifSfXtNKEDADAfYcpjzd0jcpKKsxNvwc4Z/rQUVRcGaUIHAOAkCFMeO9IVXlk8kUempMgVfSyPAADAGxCmPDazTUvCh6nSLDV1j3BFHwAA8xCmPNbYNaTsjNSE30R4bWm2pkNOjV3s0QcAwFyEKY8d6RxWYYKPSknS+rLwkgj7jw94XAkAAImFMOWxxq7hhG4+n7GqOKj0VJ/2tRGmAACYizDlof6RSXUn+LIIM9JSfFpXmq39x2lCBwBgLsKUhw53hZcaSIYwJUkby3O07/iAnHNelwIAQMIgTHmoMUmu5JuxcUWOeoYn1D4w7nUpAAAkDMKUh450DSnVZyoIJn7PlBQOU5K073i/x5UAAJA4CFMeauwa1sqCgFJ8ibnB8Xzry8J79NGEDgDArxGmPHSkc1i1RUGvy4hatj9N1YUB7WN5BAAAZhGmPBKKLIC5qjh5wpQUaUJnZAoAgFmEKY+09Y9qfCqk2qIsr0s5KxvLc9TUPaKh8SmvSwEAICEQpjxyOHIlX9KNTEWa0A8w1QcAgCTClGfq28OLX64tzfa4krPz6yv6CFMAAEiEKc8cah9UUVZG0iyLMKMsx6/8QBp9UwAARBCmPHKwfUhrS5OrX0qSzEwbV+QwMgUAQARhygOhkFND+2DSTfHN2Fieo4MnBjU1HfK6FAAAPEeY8sCxvlENT0wnbZjaUJ6j8amQGruGvS4FAADPEaY8UN8x03yefNN8Ek3oAADMRZjywMETQ5KkuiQdmVpdnKX0FB9N6AAAiDDlifr2QZXl+JWbmeZ1KeckLcWntWVZjEwBACDClCcOdQyqLkmn+GbMbCvjnPO6FAAAPEWYWmTTIaf69iGtS9Ipvhkby3PUPTyh9oFxr0sBAMBThKlF1tIzovGpUNJeyTfj/MpcSdLu1j6PKwEAwFuEqUV2cGYbmbLkDlPnrchVeopPO4/2el0KAACeIkwtspk9+epKkrtnyp+Wok0VOXqFMAUAWOZSvS5guTnYPqSKvEwFM5L/o7+kOl/ffvGoJqZCSk8llwOA1x7e1hz3c95+2cq4n3Op4TfgIqtvH9S6JJ/im3HxynxNTIW0t63f61IAAPAMYWoRTU6HdKRzOOmXRZhxcXW+JDHVBwBY1ghTi+ho97AmpkNaW7I0RqZKc/yqzM/UzmbCFABg+Ur+xp0kcqg9vI3MUpnmk8J9Uy8d6ZZzTmbmdTkAkPCccxqZmFb30IS6hsfVPTSh7qFxdQ9PqGsofH9kYkrjUyFNToc0Oe3UOzKhQHqqAukpCqSnqCgrQ7WFQeUF0vjZmwAIU4vo4IlBmYX3tlsqLqnO1492tamtf0wVeZlelwMAnnLOaWBsSif6x9TWP6oT/WM63jeq4/1jkVv4+5GJ6ZO+PiPVp2BGqvypPqX4bPY2Pe3U3j+mkYkpjU5OKxTZfCI3M001hQGdtyJXG1fkyEew8gRhahHVdwyquiCgzPQUr0uJm4tX/rpvijAFYKl7eFuzxian1T08oZ7hCfVERpT6RifVPzKp/rFJTUyFfuM1PpNKsv0qy/VrXVm2rltbopKcDBVlZei11j4FM1KVlZGqYEaq0lLO3H0Tck4dA+Nq7B5WU9ewjnQOa3drv0qyM3T9uhJdUJlLqFpkhKlFdPDEoOqSfOXz+daXZSszLUU7j/bqXZtXeF0OAMTF+NS0GruGVd8+pIaOITV2Detoz4jq2wffMKoUzEhVfiBNJTkZWluapZzMNOXOuWX705Tie2O4mZgKaV1ZzlnX5jNTWW44nF2xqlAh5/Raa7+eOdihH+xo0c/3t+ttm8q1ccXZnxvnhjC1SManptXUPaKbN5V5XUpcpab4dGFVHlf0AUhKw+NTOtw5pO+8eFSdg+PqGBhTx+C4eoYnNLONu0nKC6SpIJiuTStyVRBMV0EwXYVZ6SoIpCsjzdvZBp+ZNlfl6fzKXO1rG9DTBzr0L9uO6vp1xbpxQymjVIuAMLVIGruGNR1ySb8n38lcUp2vr//isEYmphRI5z8pAImnb2RCDR1Dqu8Ymv16uGNIx/pGZ49JMVNhVrrKc/26oDJPJTkZKskOT8dFM/3mNZ+ZNlXkan1Ztn60u03PHuzUif4xfWBLlfweB76ljt98i+T1YwOSpI3lS2/Y9ZLqfE2HnPa09uvyVYVelwNgmZqaDul4/5iOdo/ocOdMaBpUQ8ewuobGZ4/zp/m0ujhLl9bk67aSKq0pydaBEwMqDGacdDou2aSm+PTeiyq0Ii9TP9nTpq89e1gfubxaRdkZXpe2ZBGmFsme1j4F01O0agldyTfjopV5ksJN6IQpAAtpfGpaLT2jau4ZVlPXiI52h3uZjnaPqLV3RJPTbvbY7IxUrSnN0pvWF2tNSZbqSrK1piRLFXmZ8s0LTT3DE4v9R1lQZqYrVhWqNCdD39vWrG8836hPXLdaOZlpXpe2JBGmFsnu1n5tqshdEv/qmS8vkK7VxUHtpG8KQIwe3tasyemQeoYnwusvDY+rK/K1Z2hC/aOTcnOOz8pIVXVhQBvKs3XzpjLVFAa0siCo2qJgOEi83DJ77MzyBMvJqqIsfeyqWj3w3BF958Um/fa1q5SRypRfvEUVpszsZklflpQi6UHn3BfnPf8hSX8SuTsk6RPOud3xLDSZTUyFtL9tQB+9qsbrUhbMJdX5+tm+dhbvBBCViamQWnpH1NQ1rMauYTV1h0eaXm/rV//IbwamQHqKCoPpqikKhhu/I7eCrAwF01N+42fOdCjco9rYNbz4f6gEtSIvU7durdI/v3hU39/eog9fXk1TepydMUyZWYqk+yTdJKlV0nYze9w5t2/OYY2SrnPO9ZrZ2yQ9IOmyhSg4GR08MaiJ6ZAuqMz1upQFc0l1vn6wo1WHO4e0ZolslwMgdtMhp+aeER08MaADJwZ1MHJr6h6eXXhSknL8qaotCqq6IKDClRkqykpXYTDc/L2U1ubzyvqyHL3zgnL9x57jevK143rHBSxlE0/RjExtldTgnDsiSWb2iKRbJM2GKefcC3OOf0lSZTyLTHa7W/skSZsr8zyuZOFcubpIkvSLQ12EKWCZ6h4a14ETg5HQNBAOTu2DGpsML2JpJtUUBrWuNFvvuKBcNYVB1RSFp+TyI9uiPLyt2eM/xdJ1xeoidQ1P6PnD3SrKztBltfS4xks0YapCUsuc+606/ajTnZKejKWopWZPa5/yA2mqzF+6K4RXFQS0piRLzx7s0J1X13pdDoAFNDY5rfr2IT30QpNO9I+qfWBcJwbGNDQ+NXtMMD1FZbl+XbIyX6U54QUmS7L9Sk/99RID41Oh2ZEqLI53nF+ursFxPfHaca0pzlJhFlf4xUM0YepkE6vuJI/JzG5QOExdfYrn75Z0tyStXLkyyhKT357Wfp1fmbfke4luWFesb79wVMPjUwpmcG0DkOxCkSm6mem5A5HRprlTdKk+U2mOX2tLs8Orcuf4VZqToWw/V40lIp+Z3ntxpf7+vw7psVeP6c6ra+mfioNofuO1Sqqac79SUtv8g8zsAkkPSnqbc677ZCdyzj2gcD+VtmzZctJAttSMTkyrvmNIN20s9bqUBXfDuhL90y8b9cLh7mXx5wWWku6h8UhgigSn9kEdOjGo0cnw1ilmUnVBQOvKsvXOzSu0oSxb9e1DKshK55dxksnNTNPbzy/Xv716TC839rCkTRxEE6a2S6ozs1pJxyTdKun2uQeY2UpJj0m6wzl3KO5VJrG9bf2aDjldsIT7pWZsqSlQMD1FzxzsIEwBCWpmiu7AnJ6mAycG1Tn460UtA5EpuotW5kVGmsK3uVN0vSOTLAKZxLZU5+u1Y/366d4TWleWrfxAutclJbUzhinn3JSZfVLSUwovjfBN59xeM7sn8vz9kv5cUqGkr0Wmsqacc1sWruzksbu1X5K0eQlfyTcjPdWnq+uK9OyBDpZIADwWCjm19Ian6A4cH9TB9gG93Nir7qHx2T6NVJ+pJCdDVfmZurQ6X6WRabqsjFT+/13izEzvuahCX/55vf7t1WP62JU1/J3HIKrGFufcE5KemPfY/XO+v0vSXfEtbWnY09qnshy/SnL8XpeyKG5YV6Kn9rbrUPuQ1pVxVR+wGMYmp3XwxKD2HR/Q3rZ+7W0LjzqNTPx6im5lQUAl2Rk6vyJ3trepkCm6ZS0/kK6bzyvT47vb9MrRXm2pKfC6pKRFl/AC29Pav6TXl5rv+nUlkqRnDnYQpoAF0D8yqb3H+7WvbUD72ga0t21A9R2Dsw3hGak+ledmanNl3mxoKsnJYNVrnNTW2gLtae3TU3tPaFNFLhsinyPC1ALqH51UY9ew3nfJ8ll2qyzXrw3lOXrmQIfuuW611+UASa1/ZFJ7jvVpd0ufXjsWHnFq7R2dfb4kO0PnrchReZ5fK3IzVZ7rV36Q0SZEz2emd5y/Qvc924Wji2YAABl+SURBVKBfHOrUW88r87qkpESYWkCvHwv3Sy2nkSkpvETCPz53RANjk8rh8mggKmOT09rbNqA9reHwtLu1/ze2RCkMpmtFXqbOK89ReV44OLH8AOKhIj9TF1Xl6fmGLm2tLaAZ/RwQphbQzMrnF1Qs/Sv55rphfYm+9uxh/aq+S28/v9zrcoCEMx1yaugY0u6WPu1q7dOe1j4dOD6oqchcXWlOhjZX5qmuJEuV+QFV5GWypQoW1E0bS/XasX49tfeEbr10+awDGS+EqQW0p6VfNYUB5QaW178eL6rKU44/Vc8c6CBMYdlzzulY36h2t/Rrd2ufdrX0aVdznyamw1us+NN8qsjL1FVrilSVn6mK/IByM5fXzwx4Ly+QHr4a+2Cnrlo9oqqCgNclJRXC1ALa09qnS5bA1RHR7pV1+2Xhf82kpvh07dpiPXuoU6GQk89H/waWj6HxKe2JhKZXm8NfZ9ZwSk/xaeOKHF1cna+q/ExV5ge4og4J47q6Ym1v6tUTrx3X3deuYqmEs0CYWiCdg+Nq6x/Tx5dZv9SMN60v0Y/3HNerLX26pDrf63KABREKOTV0DunV5l7taunTMwc61T4wNruOU2EwXVUFAV2+qlBV+Zkqy/Ur1ec77TkBr2SkpejGDSX60a427W0b0KaK5fn761wQphbIK0d7JEkXVi2vfqkZN20sVWZain64o4UwhSWja2hcu5r79GpLODztbumf3dw3x5+q0hy/Nq4oUVV+QFX5mQqwRyWSzJbqAr14uFs/23dCG8pzlMLMQlT4P32BPFffpayMVG1epmEq25+md15Qrv/Y3ab/8c6NbHyMpDM+Na19bQN68JeNaukdUUvPiHpHJiVJPpPKcvw6b0VOODgVMF2HpSHFZ3rLxlL9y7ZmvdrMQp7R4jfcAnDO6blDnbpidaHSUpbvkP6tW6v0w1da9ZM9x/WBS6vO/ALAQ91D49pxtFevHO3VjqYevX5sYLZJPDczTZX5mbp8VeHs1XVz96kDlpIN5TmqzM/Uzw90LNvZlbNFmFoAR7tH1No7qt+5dpXXpXjq4pX5Wl0c1Pd3tBCmkFCcczrSNaxXmnq142iPdjT16khkTaf0FJ/Or8zVR6+q0cUr89TYNcLVdVhWzExv2Vimbz7fqJebevSRK2u8LinhEaYWwC/rOyVJ19QVe1yJt8xMt166Uv/rif2qbx9UXSnby8Ab41PTev3YgF452qPtTb3aebRX3cMTkqTMtBRVFwb01vPKVFMY0Iq8zNkR5Z7hSYIUlqXVxUGtKgrqmYOdGh6folXjDPh0FsBz9V2qKshUdSHrdLzn4gr9zU8P6PvbW/Rn79zodTlYJvpHJvVKczg4vdLUq12tfZqYCk/Z1RQGdP26EoWcU3VBQEXZGfQ6AfOER6dKdf9zR/TQC02694Y1XpeU0AhTcTY5HdKLh7v1rgtXsEaHpKKsDN20sVSPvXpMf3zzevpMEHfOObX0jGp7U492HO3Vz/e3qyOyrpPPpBV5mdpaU6CVBQFVFwbYggWI0srCoNaXZev+XxzWhy+rXnYLUJ8NwlSc7Wrp09D4lK6tK/K6lITxwUur9OTrJ/Rf+9tZER0xm5wOaV/bgHZEGsV3HO2dXRQz25+q8ly/NlflqbowoMq8AAEeiMFNG0v1lacb9MAvD+uP3rre63ISFmEqzn55qFM+k65YTZiacU1dsVbk+vXI9hbCFM7awNikXm3u00PPN+po94haekc0OR1eFjM/kKbqwqCuXF2o6oKgSnKYsgPiqTw3U+/avELf/FWTPnplrYqzM7wuKSERpuLsufouXViVR9PqHCk+0/u2VOkrT9erqWtYNUVBr0tCAmsfGNO2xh5tbwyPOh04MSDnJFN4ym5LTYGqCwKqLgzy/xmwCD5701r95LXjuu+ZBv3Pd53ndTkJiTAVR30jE9rT2qdPvanO61ISzocvW6l/eu6I/u9/HtJXbrvI63KQIJxzau0d1bbGHr3c2K2XG3vU1D0iSQqmp+ji6nx9+s11urSmQIfaB5WRmuJxxcDyU1sU1Ae2VOrhbc2665paVeZzcdV8hKk4euFwt0JOuoZ+qTcoyfHrrmtq9ZWnG3TX1bXLdmX45c45p8OdQ5HwFL4d7x+TFF6ioLYoqLefX67awqDKcv2zW1kc7R4hSAEe+tSb6vTozmP6h5/X62/ft9nrchIOYSqOflnfqexlvIXMmdx97So9vK1ZX3hyv77325dzteMyMB1yOnBiQC839mjbkR5tb+qZXd+pJDtDW2sL5DNTbVFQxSxRACSsFXmZuuPyan3r+Ub9znWrtbo4y+uSEgphKk7CW8h0LfstZE4n25+m33tznT7/+F49e7BTN6wv8bokxFko5HTgxKCeb+jSv77SqqM9wxqbDK/vlB9IU21RUNetLVZtUVAFwXQCNZBEfvf61Xrk5WZ96T8P6b7bL/a6nIRCmIqTI13DOtY3qnuuX+11KQnttq0r9a3nG/XFJw/o2rXF7Ei+BLT0jOj5hi79qqFLLxzuVk9k5KkoK0PnV+SptiigmsKg8gLpHlcKIBaFWRm68+pa/cPTDfrEdf3aVJHrdUkJgzAVJ/+285h8Jt24gdGW00lP9emP3rpe9z68U4/ubNUHtrBnX7L5p+eO6HDnkA53Dutw59BseMrxp2p1cZbetL5Eq4uzuNIOWILuunaVvv3iUX3xyQP65zu3MrocQZiKg6npkH74SouuX1ei8txMr8tJeG8/v0ybq/L0pZ8d0jvOL2fPpwQ3OjGtl5t69EJk9Glv24AkKSPVp1VF4TWe1hRnqTg7gx+swBKX40/TZ26s01/8xz79bF+73npemdclJQR+i8XBLw51qn1gXH95C6Ms0TAz/fk7N+j997+oP3l0j75y20X8Ek4gU9Mh7TnWr+fru/T84S7tPNqniemQ0lJMF6/M140bSrWmJEsVeZlM0wLL0Icvr9b3Xm7WX/9kn65bWyx/GlfaEqbi4JHtLSrKytCbaKiO2iXVBfrDt67T3/70oC6sytNd16zyuqRlyzmn+o4hffm/6nW4c0iNXcMaj2wKXJ7r12W1BVpdkqWawiBbswBQWopPn/9v5+lDD27TN37VyCbIIkzFrGNgTE8f6NBvX7Nq2V/F9/C25qiOu/2ylZKkT1y3Wrtb+vSFJw9oU0WuLl9VuJDlYY7W3hG90NCt5w+Hm8Zn9rYrCKbrgso8rS4OalVxlrKYggVwEletKdLN55Xpq0836L0XVyz7Fhd+UsboX3e2ajrk9MFLmeI7W2am//P+zXr3fc/rkw/v1H986upl/z/kQukaGte2Iz3h8NTQNbvKeFFWuq5cXaSr1hSqa3BC+UGuuAMQnf/+jg16+mCHvvjkAX351uW9swVhKgahkNP3t7fostoC1bLf3DnJ9qfpH++4RLd89Xn97nd36rt3XaZAOv9ZxsI5p6PdI9reFF4kc0dTr450DUuSsjJSdVltge64okZXrSnUutLs2X61aEcWAUCSqgoC+p1rV+krTzfow5dX69KaAq9L8gy/tWLwUmO3jnaP6LM3rvW6lKS2piRb/+f9m3Xvwzt16wMv6cHf2qKSbL/XZSWNqemQDpwY1MuNPdpxtEfbm3pnp+0y01JUXRjQzeeVqaYwoIr8wGzT+M6jfdp5tM/L0gEkuU9cv1qPvtKqzz26Rz/+1DXKTF+ezeiEqRh8f3uLcvypunkTl4bG6m3nl+uBO7boU997Ve+57wV986OXal1ZttdlJaSBsUn9/X/Wq7lnWEe7R9TcMzLbMJ4XSFNNYXi5gppCtmgBsLAC6an62/dt1oe/sU1ffHK//uKWTV6X5AnC1DnqHZ7Qk6+f0G2XVnFZaJzcuLFUP7znCn38oe1639df0H0fuljXri32uixPTU6HdOD4oHa19GpXS792tfTqcGd4ys4kleb4dWFVnmoKg6ouDLDKOIBFd3Vdke68ulbf+FWjrl9fohvWLb8r2wlT5+j//dchTU2H9KHLq70uZUnZVJGrf7/3Kn38oe366Lde1ocuq9Znb1qrgmXQGO2cU0vPqHa19mlXc592tfRqb9vA7KhTYTBdF1bl6d0XVqhneEKV+YFlO6QOILH80VvX6Vf1Xfrjf92jn376GhVmZXhd0qIiTJ2D11r79S8vHdUdl1drbSlTUfG2Ii9TP7znCv3dUwf13W3N+vddx/TpN9fpI1fULJl1jqZDTk3dw9rbNqB9bQPad3xAe4/1qzuyNUuqz1SRl6lLawpUmZ+pqvyA8gJps83iy+0HFYDE5k9L0d/feqFu+erz+tPHXtM/3nHJslqMmTB1lkIhpz/70esqCGbo99+yzutylqxsf5r+8pZNuuPyav3VT/brr3+yX//80lHdtnWl3nNRhUpzkqNB3TmnnuEJHekaVn37kPYd79fetgEdOD6o0clpSVJaiqmuJFs3rC/R5HRIVfkBleb4WV0cQFLZUJ6jP3rrOv2vJ/bru9ua9eFlNHNDmDpLj2xv0e6WPv2/D25mI9dzdDaX4N9+2Up9+2OX6tmDnfrK0/X64pMH9Lc/PaCr1hTp3RdW6PLVhVqR6/f8X0AjE1M61juqw53DOtI1pMMd4a9HOofVPzo5e5w/zafy3ExdtDJP5bmZWpHnV3F2hlJ9S2PEDcDydufVtfpVQ5c+//heleb4ddPGUq9LWhSEqbPQMzyhv33qgC6rLdC7L6zwupxlw8x0w/oS3bC+RI1dw/q3na16dOcx/cEPd0uSynL8uqQ6XxdW5am6MKCqgvAt1tW7p0NOg2OTGhidUv/opHpHJnSif0zH+8d0YmA0/LV/TG19oxoYm/qN15bmZGhVUZbeeUG5VhVnqaVnRMVZGb8xVQcAS43PZ/rahy7W7Q9u070P79R3Pr51WexuQZg6C3/z5AENjU3pr969iV+IHqktCur337JOn7lxrfYdH9ArR3tnbz957fhvHJvjT1VuIE05/jRl+1OVlZGqFJ8pxWfyWfg2OR3S5HRIE9NOY5PTGhid1ODYVPjr+NQpqgivHF4WGRHbuCJHuf405QbSVJSVoaKsjDdc4UlvHYDlIpiRqm999FK9//4X9Nvf3qHv3X25NlXkel3WgiJMRenHe9r0/R0tuvvaVfxiTAA+n2lTRa42VeTqt66skRQeOWzpGVFL74hae0d1PDJiNBOQjvWNKRRymnZOoZBTyDmlpfiUnuqb/VqZH1BuZpqO94/Kn5aizMjNn5aizPQU5WamKcefqtRlvg8jAJxOQTBd/3znZXrf11/QR7/1sn7wO1doVXGW12UtGMJUFH6065g++/1d2lpToM/cWOd1OcvKuW5xkuNPU05ZLD1t+TG8FgCwIi9T37nzMn3gH1/Ue772gr70gc1684al2UPFP6/P4LGdreEgVVughz5+KfvGAQAQpTUlWfq3371SlfmZuvPbO/Q3Pz2gqemQ12XFHWHqNH64o0V/8MPdumJ1ob710a0EKQAAzlJ1YVCPfuJK3bZ1pb7+7GF96MFtOt4/6nVZcUWYOomWnhF9+pFX9Uf/ukdXrynSN37rUlaaBgDgHPnTUvSF956vL31gs/a09uu6v3tWf/6j19XWtzRCVVRDLWZ2s6QvS0qR9KBz7ovznrfI82+XNCLpo865nXGudcH1DE/oK0/X619eOqoUn+l3r1+t33tzHXvvAQAQB++9uFKX1hToa8826OFtzfrey816/5Yq3XF5tdaXZSftlfJnDFNmliLpPkk3SWqVtN3MHnfO7Ztz2Nsk1UVul0n6euRrwmvpGdFz9Z167lCnflXfpdHJaX1gS5U+c+NaleUmxyrbAAAki6qCgL7w3gt07w1rdP8vDusH21v18LZmrcj16/r1JXrTuhJtqclPqo3boxmZ2iqpwTl3RJLM7BFJt0iaG6ZukfQd55yT9JKZ5ZlZuXPu+BtPtzgmpkJqHxjT0PiUhsenNDwxrcGxSR3rHZ29dP5I57Cae0YkSSty/XrXhSv08atqVcfSBwAALKjK/ID++t3n6/feXKen93fo6QMd+vdXj81exV0YTNfqkiytLs5SRZ5f+cF0FQTSlR9MVzA9VempPmWkhpe1yfanKtvv3a4k0YSpCkktc+636o2jTic7pkKSZ2HqUPug3vmVX530udzMNFXmZ+q8FTn62FU1uqauWKuLg0k7vAgAQLIqyfbr1q0rdevWlRqfmtaOpl7taxtQQ8eQDncO6cnXj6tvZPK057j9spX63+85f5EqfqNowtTJEoY7h2NkZndLujtyd8jMDkbx/gtiT+Tr/V4V8JuKJHV5XcQSw2caf3ym8cdnGn98pnH2Ia8LiMIXIrcFdsqdm6MJU62Squbcr5TUdg7HyDn3gKQHonjPZcXMdjjntnhdx1LCZxp/fKbxx2caf3ym8EI0SyNsl1RnZrVmli7pVkmPzzvmcUkfsbDLJfV72S8FAACwWM44MuWcmzKzT0p6SuGlEb7pnNtrZvdEnr9f0hMKL4vQoPDSCB9buJIBAAASR1TrTDnnnlA4MM197P453ztJ98a3tGWFqc/44zONPz7T+OMzjT8+Uyw6C+cgAAAAnAu2kwEAAIgBYcpDZnazmR00swYz+5zX9SwFZvZNM+sws9e9rmUpMLMqM3vGzPab2V4z+7TXNSU7M/Ob2ctmtjvymf6F1zUtFWaWYmavmtmPva4FywthyiNztul5m6SNkm4zs43eVrUkPCTpZq+LWEKmJP2Bc26DpMsl3ct/pzEbl/Qm59xmSRdKujlyFTRi92lJ+70uAssPYco7s9v0OOcmJM1s04MYOOeek9TjdR1LhXPu+Mym5c65QYV/UVV4W1Vyc2FDkbtpkRvNqzEys0pJ75D0oNe1YPkhTHnnVFvwAAnJzGokXSRpm7eVJL/IdNQuSR2S/tM5x2cau7+X9MeSQl4XguWHMOWdqLbgARKBmWVJelTSZ5xzA17Xk+ycc9POuQsV3i1iq5lt8rqmZGZm75TU4Zx7xetasDwRprwT1RY8gNfMLE3hIPVd59xjXtezlDjn+iQ9K/r8YnWVpHeZWZPCLRNvMrN/8bYkLCeEKe9Es00P4CkzM0nfkLTfOfclr+tZCsys2MzyIt9nSrpR0gFvq0puzrk/dc5VOudqFP5Z+rRz7sMel4VlhDDlEefclKSZbXr2S/qBc26vt1UlPzP7nqQXJa0zs1Yzu9PrmpLcVZLuUPhf+rsit7d7XVSSK5f0jJntUfgfVf/pnONSfiCJsQI6AABADBiZAgAAiAFhCgAAIAaEKQAAgBgQpgAAAGJAmAIAAIgBYQoAACAGhCkAZ2Rm05E1pvaa2W4z+30z80We22Jm/3Ca19aY2e2LV+0b3ns0sg9eQjCzD5pZg5mxthSwRBCmAERj1Dl3oXPuPEk3SXq7pM9LknNuh3Pu907z2hpJnoSpiMORffCiZmYpC1WMc+77ku5aqPMDWHyEKQBnxTnXIeluSZ+0sOtnRlnM7Lo5K6W/ambZkr4o6ZrIY5+NjBb90sx2Rm5XRl57vZk9a2b/amYHzOy7ke1sZGaXmtkLkVGxl80s28xSzOzvzGy7me0xs9+Jpn4z+3czeyUyynb3nMeHzOwvzWybpCtO8Z7nRb7fFXnPushrPzzn8X+cCWNmdnPkz7jbzH4ex78GAAkk1esCACQf59yRyDRfybyn/lDSvc65580sS9KYpM9J+kPn3DslycwCkm5yzo1Fwsj3JG2JvP4iSecpvOn385KuMrOXJX1f0gedc9vNLEfSqKQ7JfU75y41swxJz5vZz5xzjWco/+POuZ7IvnjbzexR51y3pKCk151zfx7ZL/PASd7zHklfds59N3JMipltkPRBSVc55ybN7GuSPmRmT0r6J0nXOucazazgrD9oAEmBMAXgXNlJHnte0pfM7LuSHnPOtUYGl+ZKk/RVM7tQ0rSktXOee9k51ypJkT6nGkn9ko4757ZLknNuIPL8WyRdYGbvi7w2V1KdpDOFqd8zs/dEvq+KvKY7UsujkcfXneI9X5T0382sMvLnqzezN0u6ROFgJkmZkjokXS7puZlw55zrOUNdAJIUYQrAWTOzVQqHjw5JG2Yed8590cx+onBP1UtmduNJXv5ZSe2SNivcajA257nxOd9PK/wzyiSdbBNRk/Qp59xTZ1H39ZJulHSFc27EzJ6V5I88Peacm55z7je8p3Pu4cg04DskPWVmd0WO/bZz7k/nvde7TlE3gCWGnikAZ8XMiiXdL+mrbt5O6Wa22jn3mnPubyTtkLRe0qCk7DmH5So86hOSdIekMzV7H5C0wswujbxHtpmlSnpK0ifMLC3y+FozC57hXLmSeiNBar3Co0dRv2ckRB5xzv2DpMclXSDp55LeZ2YlkWMLzKxa0ouSrjOz2pnHz1AbgCTFyBSAaGRGpt3SJE1J+mdJXzrJcZ8xsxsUHlXaJ+lJSSFJU2a2W9JDkr4m6VEze7+kZyQNn+6NnXMTZvZBSV+J9DmNKjy69KDC04A7I43qnZLefYY/x08l3WNmeyQdlPTSWb7nByV92MwmJZ2Q9JeR/qs/k/SzSB/ZpMJ9Yy9FGtwfizzeofCVkACWGJv3D0sAWDLMrEbSj51zmzwu5TdEphtnm/IBJDem+QAsZdOSci3BFu1UeHSu1+taAMQHI1MAAAAxYGQKAAAgBoQpAACAGBCmAAAAYkCYAgAAiAFhCgAAIAb/H+mH43uwwFQTAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x432 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "nb_merge_dist_plot(\n",
    "    SkyCoord(cat['ra'], cat['dec']),\n",
    "    SkyCoord(specz['ra'] * u.deg, specz['dec'] * u.deg)\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "a26f0b9d-0c4e-4457-a88b-dd0c7f278988",
   "metadata": {},
   "outputs": [],
   "source": [
    "cat = specz_merge(cat, specz, radius=1. * u.arcsec)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "720260c2-8966-47bf-8841-16673c0a8609",
   "metadata": {},
   "outputs": [],
   "source": [
    "cat=cat[cat['specz_id']!='']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "f726057d-4ed9-43ff-9376-eb14db40098d",
   "metadata": {},
   "outputs": [],
   "source": [
    "cat['redshift']=cat['zspec']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "bbf5fd12-671d-4c68-9cbc-59ef76418b4a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "83687"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(cat)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "a06a49c2-9b04-4c4c-8dfe-43993f5d8c41",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=5</i>\n",
       "<table id=\"table140195925109392\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>ra</th><th>dec</th><th>VISTA_Ks_m_detect_isPrimary</th><th>f_suprime_g</th><th>f_suprime_r</th><th>f_suprime_i</th><th>f_suprime_z</th><th>f_suprime_y</th><th>f_vista_z</th><th>f_vista_y</th><th>f_vista_j</th><th>f_vista_h</th><th>f_vista_ks</th><th>ferr_suprime_g</th><th>ferr_suprime_r</th><th>ferr_suprime_i</th><th>ferr_suprime_z</th><th>ferr_suprime_y</th><th>ferr_vista_z</th><th>ferr_vista_y</th><th>ferr_vista_j</th><th>ferr_vista_h</th><th>ferr_vista_ks</th><th>specz_id</th><th>zspec</th><th>zspec_qual</th><th>zspec_association_flag</th><th>redshift</th></tr></thead>\n",
       "<thead><tr><th></th><th>deg</th><th>deg</th><th></th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th>uJy</th><th></th><th></th><th></th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>float64</th><th>float64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>str33</th><th>float64</th><th>int64</th><th>bool</th><th>float64</th></tr></thead>\n",
       "<tr><td>149940984795368488</td><td>34.10407025372452</td><td>-5.605978251120504</td><td>1</td><td>0.34723895518347586</td><td>1.3840965971098</td><td>3.9475580319571364</td><td>6.025872568788343</td><td>6.975954853809888</td><td>5.983723849711964</td><td>6.675733786120917</td><td>10.721057893784982</td><td>14.98851788474913</td><td>18.405063522879285</td><td>0.017079866914470943</td><td>0.029173640006174003</td><td>0.03066456777209846</td><td>0.06388906821369413</td><td>0.16999562415197259</td><td>0.011375732332700046</td><td>0.01812977345981327</td><td>0.029668155255377117</td><td>0.04940991001460621</td><td>0.06957982257254518</td><td>VIPERS-105217258</td><td>0.6381</td><td>4</td><td>False</td><td>0.6381</td></tr>\n",
       "<tr><td>149940984795369896</td><td>34.03893298344629</td><td>-5.571325213602774</td><td>1</td><td>2.2021062818504467</td><td>2.8041652666715726</td><td>3.5174863865610306</td><td>5.526770264870579</td><td>6.770216848940957</td><td>3.804265082589713</td><td>4.590691597626179</td><td>5.793308461435549</td><td>8.10252777223114</td><td>8.973638516449446</td><td>0.031022190625259245</td><td>0.055287615142117455</td><td>0.06814548605888027</td><td>0.15140730599479107</td><td>0.317848329753599</td><td>0.014083530133696762</td><td>0.019421278083787123</td><td>0.03644891667642288</td><td>0.06227872095955973</td><td>0.08572107358215024</td><td>VIPERS-114013662</td><td>1.3838</td><td>2</td><td>False</td><td>1.3838</td></tr>\n",
       "<tr><td>149940984795384280</td><td>34.07550074924483</td><td>-5.63786147015687</td><td>1</td><td>7.908911855784077</td><td>37.724888903945974</td><td>71.32971427834585</td><td>99.74206019180589</td><td>117.40106197113136</td><td>61.06424558897672</td><td>87.27580447200684</td><td>132.73155447859673</td><td>196.13805315489753</td><td>247.99742952619835</td><td>0.04336680535575246</td><td>0.07633462747352954</td><td>0.06154238812315937</td><td>0.11829435786263728</td><td>0.2562789379547039</td><td>0.01631740354432525</td><td>0.02660168944989076</td><td>0.037652646404786064</td><td>0.06915811945471749</td><td>0.1048592132537531</td><td>OzDES_001924</td><td>0.41443</td><td>4</td><td>False</td><td>0.41443</td></tr>\n",
       "<tr><td>149940984795384292</td><td>34.11936529228525</td><td>-5.640340119641305</td><td>1</td><td>4.055106681849048</td><td>10.153298519473163</td><td>27.712999267158203</td><td>40.73576718298628</td><td>55.56889384470375</td><td>nan</td><td>nan</td><td>38.85576229544384</td><td>nan</td><td>75.835338283525</td><td>0.05970076247366647</td><td>0.057817778655414764</td><td>0.047464833136954</td><td>0.09978751602416282</td><td>0.24092955651370862</td><td>nan</td><td>nan</td><td>0.06695707503320905</td><td>nan</td><td>0.1851097861005107</td><td>VIPERS-105209717</td><td>0.7046</td><td>5</td><td>False</td><td>0.7046</td></tr>\n",
       "<tr><td>149940984795384296</td><td>34.05322981757711</td><td>-5.627701179292539</td><td>1</td><td>3.0493001722782167</td><td>6.750550582274663</td><td>9.827000670781645</td><td>12.393846109540867</td><td>21.746535343817182</td><td>9.486476727247677</td><td>12.851044609859064</td><td>18.581404448347776</td><td>25.700209515855086</td><td>28.660933057947968</td><td>0.018789770737188287</td><td>0.034602460221046696</td><td>0.03179320131099759</td><td>0.0710747402489922</td><td>0.18186859704618394</td><td>0.01058423994539642</td><td>0.017349735800718302</td><td>0.02926319464407901</td><td>0.04940178362927275</td><td>0.06864640159051165</td><td>SDSS 1237679321787138765</td><td>0.4584428</td><td>4</td><td>False</td><td>0.4584428</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=5>\n",
       "     help_id               ra        ... zspec_association_flag  redshift\n",
       "                          deg        ...                                 \n",
       "      int64             float64      ...          bool           float64 \n",
       "------------------ ----------------- ... ---------------------- ---------\n",
       "149940984795368488 34.10407025372452 ...                  False    0.6381\n",
       "149940984795369896 34.03893298344629 ...                  False    1.3838\n",
       "149940984795384280 34.07550074924483 ...                  False   0.41443\n",
       "149940984795384292 34.11936529228525 ...                  False    0.7046\n",
       "149940984795384296 34.05322981757711 ...                  False 0.4584428"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cat[0:5]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "458ded97-b58e-4f73-9b54-bd39ea419232",
   "metadata": {},
   "source": [
    "## Add EBV column\n",
    "\n",
    "This is used to correct for galactic extinction."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "0869d235-f98d-4e82-9abb-28a59538524f",
   "metadata": {},
   "outputs": [],
   "source": [
    "cat.add_column(\n",
    "    ebv(cat['ra'], cat['dec'])\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "96e77d85-4633-44aa-9d0e-ca151c3d4670",
   "metadata": {},
   "source": [
    "## Convert to CIGALE and write catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "55aff124-9166-4004-adb0-2c3c860c2463",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Correction for galactic extinction\n",
    "best_catalogue = correct_galactic_extinction(cat, inplace=True)\n",
    "# Convertion to CIGALE format\n",
    "best_catalogue = convert_table_for_cigale(best_catalogue, inplace=True, remove_zerofluxes=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "9f3e6972-b746-4cb1-bf4f-2a053b06e67f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=5</i>\n",
       "<table id=\"table140195925110288\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>id</th><th>ra</th><th>dec</th><th>ebv</th><th>redshift</th><th>suprime_g</th><th>suprime_g_err</th><th>suprime_r</th><th>suprime_r_err</th><th>suprime_i</th><th>suprime_i_err</th><th>vista_z</th><th>vista_z_err</th><th>suprime_z</th><th>suprime_z_err</th><th>suprime_y</th><th>suprime_y_err</th><th>vista_y</th><th>vista_y_err</th><th>vista_j</th><th>vista_j_err</th><th>vista_h</th><th>vista_h_err</th><th>vista_ks</th><th>vista_ks_err</th></tr></thead>\n",
       "<thead><tr><th></th><th>deg</th><th>deg</th><th></th><th></th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>149940984795368488</td><td>34.10407025372452</td><td>-5.605978251120504</td><td>0.019408485287841733</td><td>0.6381</td><td>0.0003707103401241331</td><td>1.8234369095463047e-05</td><td>0.0014492702502430339</td><td>3.054735387727702e-05</td><td>0.0040795585295098115</td><td>3.1689945529787287e-05</td><td>0.00614249646057396</td><td>1.1677576931196508e-05</td><td>0.006181673600107129</td><td>6.554094229567316e-05</td><td>0.007128869396925194</td><td>0.00017372196753342423</td><td>0.006811587734607005</td><td>1.849872186137994e-05</td><td>0.010874506057748394</td><td>3.0092789092561986e-05</td><td>0.015127951678390143</td><td>4.986955594155991e-05</td><td>0.018521238977526833</td><td>7.001902059604628e-05</td></tr>\n",
       "<tr><td>149940984795369896</td><td>34.03893298344629</td><td>-5.571325213602774</td><td>0.01935036037226215</td><td>1.3838</td><td>0.0023504957120446704</td><td>3.3112627961640014e-05</td><td>0.002935801829130571</td><td>5.788299412721663e-05</td><td>0.0036347478684042043</td><td>7.041723349384336e-05</td><td>0.003904901490083189</td><td>1.4456089838846577e-05</td><td>0.005669233414295567</td><td>0.00015531012095257475</td><td>0.006918172300662536</td><td>0.00032479454643423596</td><td>0.004683831470995584</td><td>1.9815313566878886e-05</td><td>0.005875976733322391</td><td>3.696902862520745e-05</td><td>0.008177676436429351</td><td>6.285633856478707e-05</td><td>0.009030111268119673</td><td>8.626053200722899e-05</td></tr>\n",
       "<tr><td>149940984795384280</td><td>34.07550074924483</td><td>-5.63786147015687</td><td>0.018948701557299608</td><td>0.41443</td><td>0.008430436372969434</td><td>4.6226472100993716e-05</td><td>0.03945822569863658</td><td>7.984195704700455e-05</td><td>0.07365745832959542</td><td>6.355073666769944e-05</td><td>0.06264565221890002</td><td>1.673998225465305e-05</td><td>0.10225906859798124</td><td>0.00012127953675879002</td><td>0.1199129080703559</td><td>0.0002617621358049348</td><td>0.08900941266288202</td><td>2.713009370809775e-05</td><td>0.13458599441743818</td><td>3.8178704971420046e-05</td><td>0.19791924755087797</td><td>6.978616716309361e-05</td><td>0.24952562760764904</td><td>0.00010550537175960117</td></tr>\n",
       "<tr><td>149940984795384292</td><td>34.11936529228525</td><td>-5.640340119641305</td><td>0.019494952102083473</td><td>0.7046</td><td>0.004330470615423</td><td>6.37547709330979e-05</td><td>0.010633571554853827</td><td>6.055268495217597e-05</td><td>0.028643877981386446</td><td>4.90591752872079e-05</td><td>nan</td><td>nan</td><td>0.04179375730946512</td><td>0.00010237919930399822</td><td>0.056792462842597966</td><td>0.00024623457368483365</td><td>nan</td><td>nan</td><td>0.039414391709676706</td><td>6.79197171072505e-05</td><td>nan</td><td>nan</td><td>0.07631616136994657</td><td>0.0001862834481516632</td></tr>\n",
       "<tr><td>149940984795384296</td><td>34.05322981757711</td><td>-5.627701179292539</td><td>0.018787766424658486</td><td>0.4584428</td><td>0.0032486128155618348</td><td>2.001793413886617e-05</td><td>0.007058023176420897</td><td>3.6178525473552676e-05</td><td>0.01014492362369567</td><td>3.2821774400794524e-05</td><td>0.009730038904305172</td><td>1.0855986832857301e-05</td><td>0.012703917698155636</td><td>7.285290155780674e-05</td><td>0.022207818496026925</td><td>0.00018572635730118835</td><td>0.013104122860695391</td><td>1.7691407697612162e-05</td><td>0.018838791785484567</td><td>2.9668544829878598e-05</td><td>0.025931610457753225</td><td>4.984659009111149e-05</td><td>0.028836041526832543</td><td>6.906580755516208e-05</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=5>\n",
       "        id                 ra        ...      vista_ks_err     \n",
       "                          deg        ...          mJy          \n",
       "      int64             float64      ...        float64        \n",
       "------------------ ----------------- ... ----------------------\n",
       "149940984795368488 34.10407025372452 ...  7.001902059604628e-05\n",
       "149940984795369896 34.03893298344629 ...  8.626053200722899e-05\n",
       "149940984795384280 34.07550074924483 ... 0.00010550537175960117\n",
       "149940984795384292 34.11936529228525 ...  0.0001862834481516632\n",
       "149940984795384296 34.05322981757711 ...  6.906580755516208e-05"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "best_catalogue[0:5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "a50ca80e-195d-4bcf-a985-0b966e78413d",
   "metadata": {},
   "outputs": [],
   "source": [
    "best_catalogue.write('./data/lsst_ir_fusion_sxds_specz_cigale_input_{}.fits'.format(SUFFIX),overwrite=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d2727198-e3f3-4c1d-a493-cf143a5eb529",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "a8ce751f-8c63-4a13-af3c-4980997a2c9c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "asdf\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    cat['ebv'].to(u.mJy)\n",
    "except UnitConversionError:\n",
    "    print('asdf')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "2f9cad21-8ea8-490c-b3a8-da78ce8b6e2e",
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.units import UnitConversionError"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6afd0294-ecf4-4a0f-b9ec-0c8ec64e42f8",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python (herschelhelp_internal)",
   "language": "python",
   "name": "helpint"
  },
  "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.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}