transform_to_catalog_faintest_13.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import numpy as np\n",
    "from astropy import units as u\n",
    "from astropy.table import Table \n",
    "from astropy.nddata import CCDData\n",
    "from astropy.coordinates import SkyCoord, EarthLocation\n",
    "from astropy.time import Time, TimeDeltaSec\n",
    "\n",
    "import glowing_waffles.differential_photometry as gwdf"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## This should be any one of the images in your table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: FITSFixedWarning: RADECSYS= 'FK5 ' / Equatorial coordinate system \n",
      "the RADECSYS keyword is deprecated, use RADESYSa. [astropy.wcs.wcs]\n"
     ]
    }
   ],
   "source": [
    "ccd = CCDData.read('kelt-16-S001-R001-C100-r.fit')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "C:\\Anaconda3\\envs\\ast366\\lib\\site-packages\\astropy\\table\\column.py:965: RuntimeWarning: invalid value encountered in less\n",
      "  return getattr(self.data, op)(other)\n"
     ]
    }
   ],
   "source": [
    "apass, apass_x, apass_y, apass_for_transform, in_apass_x, in_apass_y = gwdf.find_apass_stars(ccd)\n",
    "vsx, vsx_x, vsx_y, vsx_names = gwdf.find_known_variables(ccd)\n",
    "\n",
    "apass_all_coords = SkyCoord(apass['RAJ2000'], apass['DEJ2000'], unit='deg')\n",
    "apass_for_transform_coords = SkyCoord(apass_for_transform['RAJ2000'], apass_for_transform['DEJ2000'], unit='deg')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Change the file name below to match the one your magnitudes are in"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "input_mags = Table.read('calibrated_magnitudes_kelt_16_faint_13.fits')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get our \"star_id\" for Kelt-16"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "802"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "kelt16_coords = SkyCoord.from_name('kelt-16')\n",
    "\n",
    "input_coords = SkyCoord(input_mags['RA'], input_mags['Dec'])\n",
    "idx, d2d, d3d = kelt16_coords.match_to_catalog_sky(input_coords)\n",
    "\n",
    "kelt_id = input_mags['star_id'][idx]\n",
    "kelt_id"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "input_mags['mag_cal'] = np.zeros(len(input_mags))\n",
    "input_mags['zero_point'] = np.zeros(len(input_mags))\n",
    "input_mags['color_term'] = np.zeros(len(input_mags))\n",
    "input_mags['t_start_tdb'] = np.zeros(len(input_mags))\n",
    "input_mags['t_mid_tdb'] = np.zeros(len(input_mags))\n",
    "\n",
    "mags_grouped = input_mags.group_by('file')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "feder = EarthLocation(\"263.54672d\", \"46.86678d\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0\n",
      "1\n",
      "2\n",
      "3\n",
      "4\n",
      "5\n",
      "6\n",
      "7\n",
      "8\n",
      "9\n",
      "10\n",
      "11\n",
      "12\n",
      "13\n",
      "14\n",
      "15\n",
      "16\n",
      "17\n",
      "18\n",
      "19\n",
      "20\n",
      "21\n",
      "22\n",
      "23\n",
      "24\n",
      "25\n",
      "26\n",
      "27\n",
      "28\n",
      "29\n",
      "30\n",
      "31\n",
      "32\n",
      "33\n",
      "34\n",
      "35\n",
      "36\n",
      "37\n",
      "38\n",
      "39\n",
      "40\n",
      "41\n",
      "42\n",
      "43\n",
      "44\n",
      "45\n",
      "46\n",
      "47\n",
      "48\n",
      "49\n",
      "50\n",
      "51\n",
      "52\n",
      "53\n",
      "54\n",
      "55\n",
      "56\n",
      "57\n",
      "58\n",
      "59\n",
      "60\n",
      "61\n",
      "62\n",
      "63\n",
      "64\n",
      "65\n",
      "66\n",
      "67\n",
      "68\n",
      "69\n",
      "70\n",
      "71\n",
      "72\n",
      "73\n",
      "74\n",
      "75\n",
      "76\n",
      "77\n",
      "78\n",
      "79\n",
      "80\n",
      "81\n",
      "82\n",
      "83\n",
      "84\n",
      "85\n",
      "86\n",
      "87\n",
      "88\n",
      "89\n",
      "90\n",
      "91\n",
      "92\n",
      "93\n",
      "94\n",
      "95\n",
      "96\n",
      "97\n",
      "98\n",
      "99\n",
      "100\n",
      "101\n",
      "102\n",
      "103\n",
      "104\n",
      "105\n",
      "106\n",
      "107\n",
      "108\n",
      "109\n",
      "110\n",
      "111\n",
      "112\n",
      "113\n",
      "114\n",
      "115\n",
      "116\n",
      "117\n",
      "118\n",
      "119\n",
      "120\n",
      "121\n",
      "122\n",
      "123\n",
      "124\n",
      "125\n",
      "126\n",
      "127\n",
      "128\n",
      "129\n",
      "130\n",
      "131\n",
      "132\n",
      "133\n",
      "134\n",
      "135\n",
      "136\n",
      "137\n",
      "138\n",
      "139\n",
      "140\n",
      "141\n",
      "142\n",
      "143\n",
      "144\n",
      "145\n",
      "146\n",
      "147\n",
      "148\n",
      "149\n",
      "150\n",
      "151\n",
      "152\n",
      "153\n",
      "154\n",
      "155\n",
      "156\n",
      "157\n",
      "158\n",
      "159\n",
      "160\n",
      "161\n",
      "162\n",
      "163\n",
      "164\n",
      "165\n"
     ]
    }
   ],
   "source": [
    "for idx, our_mags in enumerate(mags_grouped.groups):\n",
    "    print(idx)\n",
    "    new_mags, matches, transform = gwdf.transform_magnitudes(our_mags, apass, apass_for_transform, \n",
    "                                                faintest_mag_for_transform=13,\n",
    "                                                input_mag_colum='mag_inst_r',\n",
    "                                                catalog_mag_column='r_mag',\n",
    "                                                catalog_color_column='B-V'\n",
    "                                               )\n",
    "    our_mags['mag_cal'][matches] = new_mags.data.filled(np.nan)\n",
    "    our_mags['mag_cal'][~matches] = np.nan\n",
    "    our_mags['zero_point'] = transform.parameters[0]\n",
    "    our_mags['color_term'] = transform.parameters[1]\n",
    "    t_utc = Time(our_mags['date-obs'], scale='utc', location=feder)\n",
    "    our_coords = SkyCoord(our_mags['RA'], our_mags['Dec'])\n",
    "    delta_t_tdb = t_utc.light_travel_time(our_coords)\n",
    "    t_tdb = t_utc.tdb + delta_t_tdb\n",
    "    t_tdb_mid = t_tdb + our_mags['exposure'] / 2\n",
    "    our_mags['t_start_tdb'][:] = t_tdb.jd\n",
    "    our_mags['t_mid_tdb'][:] = t_tdb_mid.jd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=166</i>\n",
       "<table id=\"table1992735490400\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>id</th><th>xcenter</th><th>ycenter</th><th>aperture_sum</th><th>annulus_sum</th><th>RA</th><th>Dec</th><th>sky_per_pix_avg</th><th>sky_per_pix_med</th><th>sky_per_pix_std</th><th>aperture</th><th>aperture_area</th><th>annulus_inner</th><th>annulus_outer</th><th>annulus_area</th><th>exposure</th><th>date-obs</th><th>night</th><th>aperture_net_flux</th><th>mag_inst_r</th><th>airmass</th><th>filter</th><th>file</th><th>star_id</th><th>mag_cal</th><th>zero_point</th><th>color_term</th><th>t_start_tdb</th><th>t_mid_tdb</th><th>vsx_name</th><th>mag_error</th></tr></thead>\n",
       "<thead><tr><th></th><th>pix</th><th>pix</th><th>adu</th><th></th><th>deg</th><th>deg</th><th>adu</th><th>adu</th><th>adu</th><th>pix</th><th></th><th>pix</th><th>pix</th><th></th><th>s</th><th></th><th></th><th>adu</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>adu</th></tr></thead>\n",
       "<thead><tr><th>int32</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>bytes19</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>bytes1</th><th>bytes28</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>bytes50</th><th>float64</th></tr></thead>\n",
       "<tr><td>781</td><td>1787.948845056098</td><td>2058.150237692041</td><td>20808.57221698761</td><td>121391.86424392462</td><td>314.2147445679225</td><td>31.656110247210144</td><td>32.26209999996634</td><td>31.9692325592041</td><td>9.08347604498193</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>6213.542022692223</td><td>-5.456263310539181</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>813</td><td>15.05631881193547</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264892617</td><td>2458323.726836484</td><td></td><td>0.030919110943438567</td></tr>\n",
       "<tr><td>784</td><td>1785.938622456044</td><td>2059.316144539606</td><td>20858.7404088974</td><td>119886.584471941</td><td>314.2148507716824</td><td>31.656105195711376</td><td>31.59096561546277</td><td>31.26988983154297</td><td>9.053706194473842</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:21:26</td><td>58322</td><td>6567.3242572796335</td><td>-5.516386360458807</td><td>1.119</td><td>r</td><td>kelt-16-S001-R001-C002-r.fit</td><td>813</td><td>15.005234437800425</td><td>20.53078419934884</td><td>-0.017932292930569062</td><td>2458323.727565684</td><td>2458323.727912906</td><td></td><td>0.029255632919843334</td></tr>\n",
       "<tr><td>784</td><td>1784.9537896367412</td><td>2060.239942443822</td><td>21061.263086795807</td><td>120885.11024236679</td><td>314.21478100810253</td><td>31.6560875312381</td><td>32.19647718973571</td><td>31.614980697631836</td><td>9.213407786951601</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:22:57</td><td>58322</td><td>6495.919952448519</td><td>-5.504516872328452</td><td>1.116</td><td>r</td><td>kelt-16-S001-R001-C003-r.fit</td><td>813</td><td>15.01128793813637</td><td>20.526286396213315</td><td>-0.020511909950228484</td><td>2458323.728618961</td><td>2458323.7289661835</td><td></td><td>0.02964272349927733</td></tr>\n",
       "<tr><td>784</td><td>1783.1993768753148</td><td>2061.4069123946338</td><td>20873.835624694824</td><td>121417.08053278923</td><td>314.2147811389474</td><td>31.656040493470496</td><td>32.61946269674147</td><td>32.11604309082031</td><td>9.489787568159514</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:24:30</td><td>58322</td><td>6117.1383551082035</td><td>-5.439285969201987</td><td>1.114</td><td>r</td><td>kelt-16-S001-R001-C004-r.fit</td><td>813</td><td>15.080823197245817</td><td>20.530155649465307</td><td>-0.019660436888239836</td><td>2458323.729695388</td><td>2458323.7300426103</td><td></td><td>0.03143177253830521</td></tr>\n",
       "<tr><td>785</td><td>1781.9027919620642</td><td>2062.327054731321</td><td>21005.655146598816</td><td>120592.66210007668</td><td>314.214782618286</td><td>31.656011348308038</td><td>32.36319065343766</td><td>31.94162654876709</td><td>9.438249448743965</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:26:02</td><td>58322</td><td>6364.892618085367</td><td>-5.482392913616714</td><td>1.112</td><td>r</td><td>kelt-16-S001-R001-C005-r.fit</td><td>813</td><td>15.033770061252127</td><td>20.525782352422304</td><td>-0.01882461414139766</td><td>2458323.7307602395</td><td>2458323.7311074617</td><td></td><td>0.030240423898910528</td></tr>\n",
       "<tr><td>785</td><td>1781.8996289049244</td><td>2063.213519686476</td><td>20566.57011795044</td><td>118674.8752913475</td><td>314.2147977308057</td><td>31.656118462713398</td><td>31.46723964718453</td><td>31.064083099365234</td><td>9.305156870437177</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:27:34</td><td>58322</td><td>6331.126275724846</td><td>-5.4766176501815975</td><td>1.109</td><td>r</td><td>kelt-16-S001-R001-C006-r.fit</td><td>813</td><td>15.040195374230713</td><td>20.527967277448095</td><td>-0.021828284615877658</td><td>2458323.73182508</td><td>2458323.732172302</td><td></td><td>0.030263093600696123</td></tr>\n",
       "<tr><td>786</td><td>1781.9960678450905</td><td>2064.1817343334387</td><td>20088.65239715576</td><td>118884.719876647</td><td>314.21486049807714</td><td>31.656116236978438</td><td>31.41591850157064</td><td>31.083864212036133</td><td>9.020285637033622</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:29:05</td><td>58322</td><td>5876.425694231124</td><td>-5.395698333386536</td><td>1.107</td><td>r</td><td>kelt-16-S001-R001-C007-r.fit</td><td>813</td><td>15.120952998078382</td><td>20.526803112606675</td><td>-0.019866499297297162</td><td>2458323.7328783534</td><td>2458323.7332255756</td><td></td><td>0.032459084644107856</td></tr>\n",
       "<tr><td>786</td><td>1780.9012803195455</td><td>2065.175090603007</td><td>20302.57527732849</td><td>118600.19006443024</td><td>314.21478413365315</td><td>31.656117568027526</td><td>31.640581857604882</td><td>31.033759117126465</td><td>9.35914577622116</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:30:37</td><td>58322</td><td>5988.713266569741</td><td>-5.416249010933759</td><td>1.105</td><td>r</td><td>kelt-16-S001-R001-C008-r.fit</td><td>813</td><td>15.097258774530607</td><td>20.523206242548962</td><td>-0.018979368275593916</td><td>2458323.7339432035</td><td>2458323.7342904257</td><td></td><td>0.03191831759392079</td></tr>\n",
       "<tr><td>786</td><td>1780.9761818454338</td><td>2065.3296235526727</td><td>20549.849619865417</td><td>117465.8494348675</td><td>314.21488175894893</td><td>31.656035546312157</td><td>31.194802083855446</td><td>30.957958221435547</td><td>9.005083850011944</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:32:09</td><td>58322</td><td>6437.653627682208</td><td>-5.494734226270779</td><td>1.103</td><td>r</td><td>kelt-16-S001-R001-C009-r.fit</td><td>813</td><td>15.018026080723029</td><td>20.521371378848638</td><td>-0.016851414497126738</td><td>2458323.7350080577</td><td>2458323.73535528</td><td></td><td>0.02975297284829818</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></tr>\n",
       "<tr><td>770</td><td>1611.0025764166523</td><td>2028.997086936642</td><td>978484.3872070312</td><td>7195702.549194336</td><td>314.2147698855968</td><td>31.65603475049374</td><td>2223.641084423466</td><td>2223.4371337890625</td><td>40.33916579421037</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:19:47</td><td>58322</td><td>-27467.140079477802</td><td>nan</td><td>1.134</td><td>r</td><td>kelt-16-S001-R001-C157-r.fit</td><td>813</td><td>nan</td><td>19.22386755074025</td><td>-1.3122305157762144</td><td>2458323.89309175</td><td>2458323.893438972</td><td></td><td>-0.03219959136840294</td></tr>\n",
       "<tr><td>771</td><td>1609.0167501979035</td><td>2028.970385588987</td><td>17660.42583656311</td><td>120162.84250068665</td><td>314.21471767273204</td><td>31.656066034512143</td><td>37.073671779148945</td><td>36.899871826171875</td><td>9.163134357192101</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:21:20</td><td>58322</td><td>888.6918505349167</td><td>-3.3447932049801885</td><td>1.137</td><td>r</td><td>kelt-16-S001-R001-C158-r.fit</td><td>813</td><td>15.732807020970863</td><td>19.83712445213971</td><td>-1.4863488176044457</td><td>2458323.89416817</td><td>2458323.894515393</td><td></td><td>0.21043499593900344</td></tr>\n",
       "<tr><td>769</td><td>1607.0763680291975</td><td>2028.9298179938169</td><td>16006.943586826324</td><td>107325.86852121353</td><td>314.2147728498083</td><td>31.65605442725308</td><td>32.97938626398661</td><td>32.56431198120117</td><td>8.827952164727288</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:22:53</td><td>58322</td><td>1087.420731441296</td><td>-3.563909232099036</td><td>1.139</td><td>r</td><td>kelt-16-S001-R001-C159-r.fit</td><td>813</td><td>15.513852815011731</td><td>19.84223292237398</td><td>-1.4960291346697931</td><td>2458323.895244593</td><td>2458323.895591815</td><td></td><td>0.16875060029672403</td></tr>\n",
       "<tr><td>767</td><td>1606.0239253348616</td><td>2028.0788539658713</td><td>15373.285944461823</td><td>102784.02175226063</td><td>314.21489624226655</td><td>31.656021188103455</td><td>31.680692012101385</td><td>31.486581802368164</td><td>8.884703995935839</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:24:26</td><td>58322</td><td>1041.278527298191</td><td>-3.516832492635058</td><td>1.142</td><td>r</td><td>kelt-16-S001-R001-C160-r.fit</td><td>813</td><td>15.568678756963507</td><td>19.844465595778527</td><td>-1.4852336786894473</td><td>2458323.896321016</td><td>2458323.896668238</td><td></td><td>0.17495033851078678</td></tr>\n",
       "<tr><td>767</td><td>1603.9189257886737</td><td>2027.0813115007163</td><td>15596.40848493576</td><td>103324.87474060059</td><td>314.21476145345406</td><td>31.656033574566514</td><td>31.761498997575146</td><td>31.498794555664062</td><td>8.630915948370786</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:25:58</td><td>58322</td><td>1227.8448487752012</td><td>-3.6957739423016793</td><td>1.145</td><td>r</td><td>kelt-16-S001-R001-C161-r.fit</td><td>813</td><td>15.352808750648256</td><td>19.806226273386642</td><td>-1.4826684845010059</td><td>2458323.8973858645</td><td>2458323.8977330867</td><td></td><td>0.14871392384966306</td></tr>\n",
       "<tr><td>767</td><td>1600.9945391005856</td><td>2027.075366027383</td><td>16011.204504013062</td><td>107841.29452800751</td><td>314.21480618925074</td><td>31.655996195719585</td><td>33.23752366072162</td><td>32.966156005859375</td><td>8.889355503121656</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:27:31</td><td>58322</td><td>974.9030415433062</td><td>-3.445318774217707</td><td>1.148</td><td>r</td><td>kelt-16-S001-R001-C162-r.fit</td><td>813</td><td>15.65602812274351</td><td>19.857938574626928</td><td>-1.4806099596626967</td><td>2458323.8984622895</td><td>2458323.8988095117</td><td></td><td>0.18826581026294584</td></tr>\n",
       "<tr><td>766</td><td>1600.9229128917393</td><td>2026.0006649220559</td><td>17137.788482666016</td><td>116840.57440662384</td><td>314.21488244433203</td><td>31.655999820292646</td><td>35.98125847035921</td><td>35.62561798095703</td><td>8.834313830952802</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:29:03</td><td>58322</td><td>860.250634720991</td><td>-3.3094777152356487</td><td>1.15</td><td>r</td><td>kelt-16-S001-R001-C163-r.fit</td><td>813</td><td>15.7748457272073</td><td>19.843868828280502</td><td>-1.4863902762975267</td><td>2458323.899527136</td><td>2458323.899874358</td><td></td><td>0.21613832227821309</td></tr>\n",
       "<tr><td>766</td><td>1596.9283741093948</td><td>2025.0790093711425</td><td>19530.575635910034</td><td>129572.53356838226</td><td>314.2147110877354</td><td>31.65599317733549</td><td>39.82762808309098</td><td>39.52415657043457</td><td>9.338279734710198</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:30:36</td><td>58322</td><td>1512.9811693227311</td><td>-3.9224990178990393</td><td>1.153</td><td>r</td><td>kelt-16-S001-R001-C164-r.fit</td><td>813</td><td>15.171330684749162</td><td>19.86638863018421</td><td>-1.5118570518404364</td><td>2458323.900603561</td><td>2458323.9009507834</td><td></td><td>0.12603747820164996</td></tr>\n",
       "<tr><td>765</td><td>1595.8771719236654</td><td>2024.0943777667505</td><td>21233.034595489502</td><td>144025.15901470184</td><td>314.2148337494568</td><td>31.6560119995558</td><td>44.258371079705995</td><td>43.61780548095703</td><td>9.259050718196388</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:32:08</td><td>58322</td><td>1211.0192195743366</td><td>-3.6807928001662433</td><td>1.156</td><td>r</td><td>kelt-16-S001-R001-C165-r.fit</td><td>813</td><td>15.40728163906529</td><td>19.851143115226407</td><td>-1.493285186133748</td><td>2458323.9016684056</td><td>2458323.902015628</td><td></td><td>0.16035605423862012</td></tr>\n",
       "<tr><td>762</td><td>1592.9315468740322</td><td>2025.065667632063</td><td>24200.955696105957</td><td>165725.3926334381</td><td>314.21482693892267</td><td>31.656108436633787</td><td>50.959213893596946</td><td>50.72254943847656</td><td>9.637534843931036</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1147.5504479857045</td><td>-3.6223446780254838</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>813</td><td>15.479855871630765</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027448203</td><td>2458323.9030920425</td><td></td><td>0.1743219979053475</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=166>\n",
       "  id       xcenter            ycenter       ... vsx_name      mag_error      \n",
       "             pix                pix         ...                  adu         \n",
       "int32      float64            float64       ... bytes50        float64       \n",
       "----- ------------------ ------------------ ... -------- --------------------\n",
       "  781  1787.948845056098  2058.150237692041 ...          0.030919110943438567\n",
       "  784  1785.938622456044  2059.316144539606 ...          0.029255632919843334\n",
       "  784 1784.9537896367412  2060.239942443822 ...           0.02964272349927733\n",
       "  784 1783.1993768753148 2061.4069123946338 ...           0.03143177253830521\n",
       "  785 1781.9027919620642  2062.327054731321 ...          0.030240423898910528\n",
       "  785 1781.8996289049244  2063.213519686476 ...          0.030263093600696123\n",
       "  786 1781.9960678450905 2064.1817343334387 ...          0.032459084644107856\n",
       "  786 1780.9012803195455  2065.175090603007 ...           0.03191831759392079\n",
       "  786 1780.9761818454338 2065.3296235526727 ...           0.02975297284829818\n",
       "  ...                ...                ... ...      ...                  ...\n",
       "  770 1611.0025764166523  2028.997086936642 ...          -0.03219959136840294\n",
       "  771 1609.0167501979035  2028.970385588987 ...           0.21043499593900344\n",
       "  769 1607.0763680291975 2028.9298179938169 ...           0.16875060029672403\n",
       "  767 1606.0239253348616 2028.0788539658713 ...           0.17495033851078678\n",
       "  767 1603.9189257886737 2027.0813115007163 ...           0.14871392384966306\n",
       "  767 1600.9945391005856  2027.075366027383 ...           0.18826581026294584\n",
       "  766 1600.9229128917393 2026.0006649220559 ...           0.21613832227821309\n",
       "  766 1596.9283741093948 2025.0790093711425 ...           0.12603747820164996\n",
       "  765 1595.8771719236654 2024.0943777667505 ...           0.16035605423862012\n",
       "  762 1592.9315468740322  2025.065667632063 ...            0.1743219979053475"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mags_grouped[mags_grouped['star_id'] == 813]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "mags_grouped.sort('t_start_tdb')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=228416</i>\n",
       "<table id=\"table1992579081832\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>id</th><th>xcenter</th><th>ycenter</th><th>aperture_sum</th><th>annulus_sum</th><th>RA</th><th>Dec</th><th>sky_per_pix_avg</th><th>sky_per_pix_med</th><th>sky_per_pix_std</th><th>aperture</th><th>aperture_area</th><th>annulus_inner</th><th>annulus_outer</th><th>annulus_area</th><th>exposure</th><th>date-obs</th><th>night</th><th>aperture_net_flux</th><th>mag_inst_r</th><th>airmass</th><th>filter</th><th>file</th><th>star_id</th><th>mag_cal</th><th>zero_point</th><th>color_term</th><th>t_start_tdb</th><th>t_mid_tdb</th><th>vsx_name</th><th>mag_error</th></tr></thead>\n",
       "<thead><tr><th></th><th>pix</th><th>pix</th><th>adu</th><th></th><th>deg</th><th>deg</th><th>adu</th><th>adu</th><th>adu</th><th>pix</th><th></th><th>pix</th><th>pix</th><th></th><th>s</th><th></th><th></th><th>adu</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>adu</th></tr></thead>\n",
       "<thead><tr><th>int32</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>bytes19</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>bytes1</th><th>bytes28</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>bytes50</th><th>float64</th></tr></thead>\n",
       "<tr><td>1422</td><td>4011.091526838012</td><td>3969.142916177813</td><td>19740.155430793762</td><td>103860.89329612255</td><td>314.5432670134417</td><td>32.01973136240193</td><td>31.808280954476942</td><td>31.461752891540527</td><td>10.803624996914255</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>5350.428135927457</td><td>-5.293886548969706</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1479</td><td>15.214845041740617</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.726455484</td><td>2458323.7268027063</td><td></td><td>0.03554203721814884</td></tr>\n",
       "<tr><td>1431</td><td>3863.942898239571</td><td>3998.9143066290862</td><td>33952.74383163452</td><td>102687.83718004823</td><td>314.51558576760186</td><td>32.01961637265539</td><td>31.504446114354263</td><td>31.2497501373291</td><td>10.520896497204099</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>19700.46818020352</td><td>-6.709106579064407</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1489</td><td>13.799244435989024</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264559236</td><td>2458323.726803146</td><td></td><td>0.010865317328612463</td></tr>\n",
       "<tr><td>1411</td><td>3963.0904975394483</td><td>3951.90486856516</td><td>15249.548798561096</td><td>104428.82190561295</td><td>314.5352479090502</td><td>32.015567615887285</td><td>31.808294335373855</td><td>31.52717399597168</td><td>10.731908677068274</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>859.8154503196401</td><td>-3.3089283230064974</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1468</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264559385</td><td>2458323.7268031607</td><td></td><td>0.21162200178674756</td></tr>\n",
       "<tr><td>1421</td><td>3923.008965618607</td><td>3969.8875728993826</td><td>20822.206637382507</td><td>101420.80074810982</td><td>314.52733781030184</td><td>32.017048316638714</td><td>31.151097199642955</td><td>30.903825759887695</td><td>10.454439681919503</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>6729.782269015484</td><td>-5.542917744838749</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1478</td><td>14.962746853484397</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.726455944</td><td>2458323.7268031663</td><td></td><td>0.028532498737706097</td></tr>\n",
       "<tr><td>1405</td><td>4015.1792324262733</td><td>3923.173770284605</td><td>30234.403473317623</td><td>104178.87234064937</td><td>314.54572817533784</td><td>32.01282162890551</td><td>31.792475467154404</td><td>31.450231552124023</td><td>10.459912367244256</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>15851.826412462999</td><td>-6.473113480632851</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1461</td><td>14.035685270542295</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.726455993</td><td>2458323.726803215</td><td></td><td>0.013127131925856725</td></tr>\n",
       "<tr><td>1417</td><td>3848.017315598605</td><td>3961.9101672544393</td><td>16185.751720428467</td><td>115460.04435215145</td><td>314.5140998526454</td><td>32.013443881563674</td><td>31.558219960992364</td><td>30.922298431396484</td><td>10.712969735233127</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>1909.1493538937557</td><td>-4.175014973009543</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1474</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264564354</td><td>2458323.7268036576</td><td></td><td>0.09620434270592455</td></tr>\n",
       "<tr><td>1419</td><td>3723.004118947109</td><td>3968.894287448843</td><td>16342.364827394485</td><td>102438.52531147748</td><td>314.49127182053485</td><td>32.01053792907989</td><td>31.50461414878146</td><td>31.270769119262695</td><td>10.508239826067275</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>2090.013158979509</td><td>-4.273287762138963</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1476</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264570193</td><td>2458323.7268042415</td><td></td><td>0.08801476817405915</td></tr>\n",
       "<tr><td>1438</td><td>3604.9472263665893</td><td>4021.055869255223</td><td>16588.801261663437</td><td>99055.38167959452</td><td>314.4680014835236</td><td>32.014766868890604</td><td>30.492322640378536</td><td>30.28482151031494</td><td>9.544357245216066</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>2794.399482765415</td><td>-4.5886364415189576</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1497</td><td>15.915617787974224</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264570445</td><td>2458323.7268042667</td><td></td><td>0.06594961678474937</td></tr>\n",
       "<tr><td>1382</td><td>3932.032866076106</td><td>3855.864058718527</td><td>16039.92415869236</td><td>101568.43381509185</td><td>314.5332432307907</td><td>31.999874377180166</td><td>31.223004554285843</td><td>31.05585479736328</td><td>10.306347741275532</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>1914.969669465072</td><td>-4.178319960213055</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1436</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.726457212</td><td>2458323.7268044343</td><td></td><td>0.09575037644968758</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></tr>\n",
       "<tr><td>107</td><td>182.0508911360561</td><td>303.05677630682095</td><td>22842.986793518066</td><td>158728.4168729782</td><td>314.0267663163187</td><td>31.347071117831405</td><td>48.82892485615355</td><td>48.830955505371094</td><td>9.676146609681743</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>753.3016015657413</td><td>-3.1653374378374917</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>134</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027719204</td><td>2458323.9031191426</td><td></td><td>0.26216596171042117</td></tr>\n",
       "<tr><td>39</td><td>541.0001708205461</td><td>135.01320125036852</td><td>22581.07301712036</td><td>158686.21473121643</td><td>314.09743360865417</td><td>31.33297236070604</td><td>48.77210934622453</td><td>48.38126182556152</td><td>10.314663448435162</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>517.0905563268607</td><td>-2.7568317266208298</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>64</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027719847</td><td>2458323.903119207</td><td></td><td>0.3810401372674711</td></tr>\n",
       "<tr><td>97</td><td>193.00622531012922</td><td>280.0069870918759</td><td>23552.252098083496</td><td>158362.16340827942</td><td>314.02960264845143</td><td>31.34390001098602</td><td>48.58026260878469</td><td>48.42496299743652</td><td>9.545562662829555</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1575.0590566276842</td><td>-3.966157316521702</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>124</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027721277</td><td>2458323.90311935</td><td>ASASSN-V J205607.11+312038.1</td><td>0.126141972104651</td></tr>\n",
       "<tr><td>41</td><td>503.9609772205891</td><td>138.00226718367384</td><td>24672.248279571533</td><td>158899.16868400574</td><td>314.0906849908972</td><td>31.33223335878</td><td>48.33081605587856</td><td>48.28989028930664</td><td>9.58830815061887</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>2807.9022000782643</td><td>-4.593870153756231</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>66</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027721426</td><td>2458323.9031193648</td><td></td><td>0.07142793901418203</td></tr>\n",
       "<tr><td>14</td><td>602.9782312415216</td><td>75.94704453974062</td><td>21924.57398033142</td><td>155915.59290409088</td><td>314.1107678715254</td><td>31.325936902368422</td><td>48.065853973192866</td><td>48.05634117126465</td><td>9.302963672902422</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>180.09392311026386</td><td>-1.6116628576991223</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>38</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027723395</td><td>2458323.9031195617</td><td></td><td>1.0873130222307894</td></tr>\n",
       "<tr><td>80</td><td>264.0154239709242</td><td>209.05195133865692</td><td>22760.32354927063</td><td>157921.01183128357</td><td>314.04501177829536</td><td>31.33534466061794</td><td>48.68451682124661</td><td>48.531917572021484</td><td>9.813231068559666</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>735.9670132262545</td><td>-3.1400610840687304</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>105</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027725724</td><td>2458323.9031197946</td><td></td><td>0.26812608072554045</td></tr>\n",
       "<tr><td>98</td><td>83.04174328655962</td><td>284.8475961628155</td><td>23363.15751838684</td><td>160498.44156837463</td><td>314.0097193604797</td><td>31.34107838620185</td><td>49.18918774371364</td><td>49.234052658081055</td><td>9.755796233160744</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1110.493235742062</td><td>-3.5867050046852746</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>125</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027726413</td><td>2458323.9031198635</td><td></td><td>0.17867986111223463</td></tr>\n",
       "<tr><td>18</td><td>306.91673956570503</td><td>80.1025961804836</td><td>23178.463884353638</td><td>156871.81446647644</td><td>314.0575777132189</td><td>31.317007021227244</td><td>48.20240802911838</td><td>48.1168098449707</td><td>9.73735439939401</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1372.2082276089386</td><td>-3.816465258361166</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>42</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.902773827</td><td>2458323.903121049</td><td></td><td>0.1442924276151822</td></tr>\n",
       "<tr><td>37</td><td>102.38009225045143</td><td>130.41017764070943</td><td>80841.77911376953</td><td>162045.0770702362</td><td>314.0190387942966</td><td>31.318082075902257</td><td>49.66731724692211</td><td>49.48946762084961</td><td>9.88967002386121</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>58372.81413972158</td><td>-7.8884417885370715</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>61</td><td>9.262596422862341</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027743097</td><td>2458323.903121532</td><td></td><td>0.004799602972304384</td></tr>\n",
       "<tr><td>10</td><td>113.92856473017127</td><td>60.970711391059524</td><td>23356.94402694702</td><td>163692.00853157043</td><td>314.02373840609386</td><td>31.30783464202709</td><td>50.31406455831275</td><td>50.1636962890625</td><td>9.865457764126704</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>595.3974621831621</td><td>-2.909932658865633</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>34</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.902775045</td><td>2458323.903122267</td><td></td><td>0.33344694042018297</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=228416>\n",
       "  id       xcenter       ...           vsx_name                mag_error      \n",
       "             pix         ...                                      adu         \n",
       "int32      float64       ...           bytes50                  float64       \n",
       "----- ------------------ ... ---------------------------- --------------------\n",
       " 1422  4011.091526838012 ...                               0.03554203721814884\n",
       " 1431  3863.942898239571 ...                              0.010865317328612463\n",
       " 1411 3963.0904975394483 ...                               0.21162200178674756\n",
       " 1421  3923.008965618607 ...                              0.028532498737706097\n",
       " 1405 4015.1792324262733 ...                              0.013127131925856725\n",
       " 1417  3848.017315598605 ...                               0.09620434270592455\n",
       " 1419  3723.004118947109 ...                               0.08801476817405915\n",
       " 1438 3604.9472263665893 ...                               0.06594961678474937\n",
       " 1382  3932.032866076106 ...                               0.09575037644968758\n",
       "  ...                ... ...                          ...                  ...\n",
       "  107  182.0508911360561 ...                               0.26216596171042117\n",
       "   39  541.0001708205461 ...                                0.3810401372674711\n",
       "   97 193.00622531012922 ... ASASSN-V J205607.11+312038.1    0.126141972104651\n",
       "   41  503.9609772205891 ...                               0.07142793901418203\n",
       "   14  602.9782312415216 ...                                1.0873130222307894\n",
       "   80  264.0154239709242 ...                               0.26812608072554045\n",
       "   98  83.04174328655962 ...                               0.17867986111223463\n",
       "   18 306.91673956570503 ...                                0.1442924276151822\n",
       "   37 102.38009225045143 ...                              0.004799602972304384\n",
       "   10 113.92856473017127 ...                               0.33344694042018297"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mags_grouped"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Identify known variables in the field "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### First, the list of variables..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table masked=True length=4</i>\n",
       "<table id=\"table1992664347088\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>_q</th><th>OID</th><th>n_OID</th><th>Name</th><th>V</th><th>Type</th><th>l_max</th><th>max</th><th>u_max</th><th>n_max</th><th>f_min</th><th>l_min</th><th>min</th><th>u_min</th><th>n_min</th><th>l_Period</th><th>Period</th><th>u_Period</th><th>RAJ2000</th><th>DEJ2000</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>mag</th><th></th><th></th><th></th><th></th><th>mag</th><th></th><th></th><th></th><th>d</th><th></th><th>deg</th><th>deg</th></tr></thead>\n",
       "<thead><tr><th>int32</th><th>int32</th><th>bytes1</th><th>bytes30</th><th>uint8</th><th>bytes30</th><th>bytes1</th><th>float32</th><th>bytes1</th><th>bytes6</th><th>bytes1</th><th>bytes1</th><th>float32</th><th>bytes1</th><th>bytes6</th><th>bytes1</th><th>float64</th><th>bytes3</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>1</td><td>605228</td><td>B</td><td>ASASSN-V J205607.11+312038.1</td><td>0</td><td>L:</td><td></td><td>13.670</td><td></td><td>V</td><td></td><td></td><td>14.240</td><td></td><td>V</td><td></td><td>--</td><td></td><td>314.02964</td><td>31.34391</td></tr>\n",
       "<tr><td>1</td><td>605239</td><td>B</td><td>ASASSN-V J205648.40+314708.0</td><td>0</td><td>L</td><td></td><td>14.600</td><td></td><td>V</td><td></td><td></td><td>15.250</td><td></td><td>V</td><td></td><td>--</td><td></td><td>314.20167</td><td>31.78556</td></tr>\n",
       "<tr><td>1</td><td>605260</td><td>B</td><td>ASASSN-V J205811.97+315125.4</td><td>0</td><td>L</td><td></td><td>11.150</td><td></td><td>V</td><td></td><td></td><td>11.650</td><td></td><td>V</td><td></td><td>--</td><td></td><td>314.54987</td><td>31.85706</td></tr>\n",
       "<tr><td>1</td><td>605262</td><td>B</td><td>ASASSN-V J205816.59+315327.0</td><td>0</td><td>L</td><td></td><td>12.980</td><td></td><td>V</td><td></td><td></td><td>13.500</td><td></td><td>V</td><td></td><td>--</td><td></td><td>314.56913</td><td>31.89084</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table masked=True length=4>\n",
       "  _q   OID   n_OID  ... u_Period  RAJ2000   DEJ2000 \n",
       "                    ...             deg       deg   \n",
       "int32 int32  bytes1 ...  bytes3   float64   float64 \n",
       "----- ------ ------ ... -------- --------- ---------\n",
       "    1 605228      B ...          314.02964  31.34391\n",
       "    1 605239      B ...          314.20167  31.78556\n",
       "    1 605260      B ...          314.54987  31.85706\n",
       "    1 605262      B ...          314.56913  31.89084"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "vsx"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Do the matching..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "mag_coords = SkyCoord(mags_grouped['RA'], mags_grouped['Dec'])\n",
    "vsx_coords = SkyCoord(vsx['RAJ2000'], vsx['DEJ2000'])\n",
    "vsx_idx, d2d, d3d = mag_coords.match_to_catalog_sky(vsx_coords)\n",
    "matches = d2d < 3 * u.arcsec"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Add VSX name column to output table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "vsx_names = np.array([''] * len(mags_grouped), dtype='S50')\n",
    "vsx_names[matches] = vsx['Name'][vsx_idx[matches]]\n",
    "mags_grouped['vsx_name'] = vsx_names"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=664</i>\n",
       "<table id=\"table1992804742648\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>id</th><th>xcenter</th><th>ycenter</th><th>aperture_sum</th><th>annulus_sum</th><th>RA</th><th>Dec</th><th>sky_per_pix_avg</th><th>sky_per_pix_med</th><th>sky_per_pix_std</th><th>aperture</th><th>aperture_area</th><th>annulus_inner</th><th>annulus_outer</th><th>annulus_area</th><th>exposure</th><th>date-obs</th><th>night</th><th>aperture_net_flux</th><th>mag_inst_r</th><th>airmass</th><th>filter</th><th>file</th><th>star_id</th><th>mag_cal</th><th>zero_point</th><th>color_term</th><th>t_start_tdb</th><th>t_mid_tdb</th><th>vsx_name</th><th>mag_error</th></tr></thead>\n",
       "<thead><tr><th></th><th>pix</th><th>pix</th><th>adu</th><th></th><th>deg</th><th>deg</th><th>adu</th><th>adu</th><th>adu</th><th>pix</th><th></th><th>pix</th><th>pix</th><th></th><th>s</th><th></th><th></th><th>adu</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>adu</th></tr></thead>\n",
       "<thead><tr><th>int32</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>bytes19</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>bytes1</th><th>bytes28</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>bytes50</th><th>float64</th></tr></thead>\n",
       "<tr><td>1127</td><td>3981.024900930935</td><td>3133.711912110588</td><td>104532.0606918335</td><td>101442.7775220871</td><td>314.56911871176715</td><td>31.890840658692333</td><td>31.127654286931467</td><td>30.994300842285156</td><td>9.364067750930671</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>90450.24164732533</td><td>-8.363939539563152</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1175</td><td>12.126367705966778</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.726465278</td><td>2458323.7268125</td><td>ASASSN-V J205816.59+315327.0</td><td>0.003386311193159725</td></tr>\n",
       "<tr><td>1071</td><td>3835.1134500924204</td><td>2943.469287848386</td><td>492788.5991668701</td><td>110800.89044713974</td><td>314.5499296601638</td><td>31.857078164023573</td><td>32.90284882853618</td><td>32.45624923706055</td><td>10.056446641754372</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>477903.70103155583</td><td>-10.171266195517138</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1116</td><td>10.317496358390535</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264682455</td><td>2458323.7268154677</td><td>ASASSN-V J205811.97+315125.4</td><td>0.001252010506684385</td></tr>\n",
       "<tr><td>1044</td><td>1888.2009196467923</td><td>2882.2815871039766</td><td>28041.29564189911</td><td>100018.43565130234</td><td>314.2015583002371</td><td>31.78546252570611</td><td>30.74093960672709</td><td>30.49066734313965</td><td>8.263990686974921</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>14134.422197155558</td><td>-6.348610359933802</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1087</td><td>14.13995071537337</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7264792533</td><td>2458323.7268264755</td><td>ASASSN-V J205648.40+314708.0</td><td>0.014458509677046733</td></tr>\n",
       "<tr><td>111</td><td>387.3365579918771</td><td>313.29412291824985</td><td>48053.009115219116</td><td>103461.68193542957</td><td>314.0295452895356</td><td>31.34387651792485</td><td>31.42676253360653</td><td>31.090402603149414</td><td>9.102557474960875</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>33835.87668777587</td><td>-7.296358794756332</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>124</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287169795</td><td>2458323.7265166207</td><td>2458323.726863843</td><td>ASASSN-V J205607.11+312038.1</td><td>0.006956718274117619</td></tr>\n",
       "<tr><td>1132</td><td>3978.345365708148</td><td>3134.935049437815</td><td>104407.17543792725</td><td>100156.27068307996</td><td>314.5691108334105</td><td>31.89082572060043</td><td>30.750724102461827</td><td>30.511476516723633</td><td>9.39499511773456</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:21:26</td><td>58322</td><td>90495.87559159531</td><td>-8.364487176918031</td><td>1.119</td><td>r</td><td>kelt-16-S001-R001-C002-r.fit</td><td>1175</td><td>12.139290990485462</td><td>20.53078419934884</td><td>-0.017932292930569062</td><td>2458323.7275417033</td><td>2458323.7278889255</td><td>ASASSN-V J205816.59+315327.0</td><td>0.0033827794818791354</td></tr>\n",
       "<tr><td>1075</td><td>3832.579448965366</td><td>2944.8305248558004</td><td>496636.5559768677</td><td>110109.05312919617</td><td>314.5499411144663</td><td>31.857088957236957</td><td>32.50487637735971</td><td>31.927480697631836</td><td>9.991011148412639</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:21:26</td><td>58322</td><td>481931.69633692177</td><td>-10.180378937182802</td><td>1.119</td><td>r</td><td>kelt-16-S001-R001-C002-r.fit</td><td>1116</td><td>10.322161901125089</td><td>20.53078419934884</td><td>-0.017932292930569062</td><td>2458323.727544668</td><td>2458323.7278918903</td><td>ASASSN-V J205811.97+315125.4</td><td>0.0012460473535955033</td></tr>\n",
       "<tr><td>1048</td><td>1885.9828031816776</td><td>2883.371024340974</td><td>27823.397190093994</td><td>97537.35446536541</td><td>314.20162991032146</td><td>31.785436576935204</td><td>29.960256093650624</td><td>29.562904357910156</td><td>8.601142130370645</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:21:26</td><td>58322</td><td>14269.696646232638</td><td>-6.358952062712039</td><td>1.119</td><td>r</td><td>kelt-16-S001-R001-C002-r.fit</td><td>1087</td><td>14.143427386706897</td><td>20.53078419934884</td><td>-0.017932292930569062</td><td>2458323.7275556778</td><td>2458323.7279029</td><td>ASASSN-V J205648.40+314708.0</td><td>0.014290265142871666</td></tr>\n",
       "<tr><td>114</td><td>384.99479692707763</td><td>314.95309246255306</td><td>48350.27584075928</td><td>100976.08187294006</td><td>314.02958310247396</td><td>31.34393711594204</td><td>30.85836694644707</td><td>30.577152252197266</td><td>9.057798615473851</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:21:26</td><td>58322</td><td>34390.27951905326</td><td>-7.314004475203867</td><td>1.119</td><td>r</td><td>kelt-16-S001-R001-C002-r.fit</td><td>124</td><td>nan</td><td>20.53078419934884</td><td>-0.017932292930569062</td><td>2458323.727593038</td><td>2458323.7279402604</td><td>ASASSN-V J205607.11+312038.1</td><td>0.006855548127122464</td></tr>\n",
       "<tr><td>1132</td><td>3977.865619318615</td><td>3136.006973719115</td><td>104312.06343460083</td><td>102337.80106922984</td><td>314.5691239348477</td><td>31.890829979352105</td><td>31.406738637146347</td><td>31.27564239501953</td><td>9.422474400267523</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:22:57</td><td>58322</td><td>90103.98960450382</td><td>-8.359775263410631</td><td>1.116</td><td>r</td><td>kelt-16-S001-R001-C003-r.fit</td><td>1175</td><td>12.135620194908102</td><td>20.526286396213315</td><td>-0.020511909950228484</td><td>2458323.728594978</td><td>2458323.7289422</td><td>ASASSN-V J205816.59+315327.0</td><td>0.003396828328599868</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></tr>\n",
       "<tr><td>1030</td><td>1697.0920941276536</td><td>2849.9346524283856</td><td>18315.172921180725</td><td>125235.62050628662</td><td>314.2015442805204</td><td>31.78546977637469</td><td>38.51718330794241</td><td>38.44398498535156</td><td>8.615355215557955</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:30:36</td><td>58322</td><td>890.4097043034526</td><td>-3.346889922399776</td><td>1.153</td><td>r</td><td>kelt-16-S001-R001-C164-r.fit</td><td>1087</td><td>14.12471711827638</td><td>19.86638863018421</td><td>-1.5118570518404364</td><td>2458323.900593541</td><td>2458323.9009407633</td><td>ASASSN-V J205648.40+314708.0</td><td>0.21154653010203212</td></tr>\n",
       "<tr><td>98</td><td>197.99324494658933</td><td>281.1145194721711</td><td>18105.081740379333</td><td>123628.33135795593</td><td>314.02960153786603</td><td>31.343892791570525</td><td>37.704612881725375</td><td>37.582550048828125</td><td>9.127951745208897</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:30:36</td><td>58322</td><td>1047.9167240420575</td><td>-3.523732139558261</td><td>1.153</td><td>r</td><td>kelt-16-S001-R001-C164-r.fit</td><td>124</td><td>nan</td><td>19.86638863018421</td><td>-1.5118570518404364</td><td>2458323.9006308587</td><td>2458323.900978081</td><td>ASASSN-V J205607.11+312038.1</td><td>0.17929932984750513</td></tr>\n",
       "<tr><td>1112</td><td>3787.9527167352517</td><td>3100.0854107295104</td><td>23679.531801223755</td><td>140086.71602773666</td><td>314.5691594073565</td><td>31.8907674690204</td><td>43.13014656026375</td><td>43.129438400268555</td><td>9.559063559797423</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:32:08</td><td>58322</td><td>4167.913173419256</td><td>-5.022711868659153</td><td>1.156</td><td>r</td><td>kelt-16-S001-R001-C165-r.fit</td><td>1175</td><td>12.579543707657342</td><td>19.851143115226407</td><td>-1.493285186133748</td><td>2458323.9016445093</td><td>2458323.9019917315</td><td>ASASSN-V J205816.59+315327.0</td><td>0.04758048794015593</td></tr>\n",
       "<tr><td>1056</td><td>3641.9638714304388</td><td>2910.221984666655</td><td>45752.761713027954</td><td>142691.65343379974</td><td>314.54992476462706</td><td>31.857057201653955</td><td>43.59621496835564</td><td>43.469970703125</td><td>9.674811507447986</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:32:08</td><td>58322</td><td>26030.29870470528</td><td>-7.011613090371983</td><td>1.156</td><td>r</td><td>kelt-16-S001-R001-C165-r.fit</td><td>1116</td><td>10.487605843869925</td><td>19.851143115226407</td><td>-1.493285186133748</td><td>2458323.901647468</td><td>2458323.9019946903</td><td>ASASSN-V J205811.97+315125.4</td><td>0.008957217967426507</td></tr>\n",
       "<tr><td>1029</td><td>1696.0686240183486</td><td>2849.0870063170582</td><td>20262.845050811768</td><td>140621.7524547577</td><td>314.2016651276276</td><td>31.785510218634233</td><td>43.25923712676382</td><td>42.967451095581055</td><td>8.94330766355854</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:32:08</td><td>58322</td><td>692.8272265548003</td><td>-3.0744775762686105</td><td>1.156</td><td>r</td><td>kelt-16-S001-R001-C165-r.fit</td><td>1087</td><td>14.411301941087924</td><td>19.851143115226407</td><td>-1.493285186133748</td><td>2458323.9016583837</td><td>2458323.9020056063</td><td>ASASSN-V J205648.40+314708.0</td><td>0.27766849844528935</td></tr>\n",
       "<tr><td>98</td><td>195.96100066832054</td><td>279.99589279214837</td><td>20216.607152938843</td><td>137892.50992965698</td><td>314.0295663687923</td><td>31.34391448851532</td><td>42.27441364349705</td><td>42.1519660949707</td><td>9.391972952517587</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:32:08</td><td>58322</td><td>1092.1129763782337</td><td>-3.568584129359999</td><td>1.156</td><td>r</td><td>kelt-16-S001-R001-C165-r.fit</td><td>124</td><td>nan</td><td>19.851143115226407</td><td>-1.493285186133748</td><td>2458323.901695705</td><td>2458323.902042927</td><td>ASASSN-V J205607.11+312038.1</td><td>0.17597524853284996</td></tr>\n",
       "<tr><td>1109</td><td>3785.021643076381</td><td>3100.0952943883944</td><td>26207.028230667114</td><td>162080.23638153076</td><td>314.5691933997279</td><td>31.890761997285306</td><td>49.89945617459544</td><td>49.56210136413574</td><td>10.179819997869881</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>3633.0460798492895</td><td>-4.873592476623279</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>1175</td><td>12.750004471577281</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.902720932</td><td>2458323.903068154</td><td>ASASSN-V J205816.59+315327.0</td><td>0.05596242871036663</td></tr>\n",
       "<tr><td>1053</td><td>3638.9205418558713</td><td>2910.265401649082</td><td>48697.16188621521</td><td>163903.4020500183</td><td>314.54993639670187</td><td>31.85704421116024</td><td>50.18069327020513</td><td>49.87424850463867</td><td>10.212886151543175</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>25995.951070735642</td><td>-7.0101794880913</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>1116</td><td>10.510881020213631</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027238917</td><td>2458323.903071114</td><td>ASASSN-V J205811.97+315125.4</td><td>0.00915560066834782</td></tr>\n",
       "<tr><td>1026</td><td>1692.9911226337838</td><td>2849.009035103361</td><td>22954.716468811035</td><td>161000.73784828186</td><td>314.2016670050808</td><td>31.78545328377892</td><td>49.70097566556886</td><td>49.64897537231445</td><td>9.364158337596784</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>470.52478489478017</td><td>-2.6543714729844363</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>1087</td><td>14.853314945125792</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027348105</td><td>2458323.9030820327</td><td>ASASSN-V J205648.40+314708.0</td><td>0.42033236398949014</td></tr>\n",
       "<tr><td>97</td><td>193.00622531012922</td><td>280.0069870918759</td><td>23552.252098083496</td><td>158362.16340827942</td><td>314.02960264845143</td><td>31.34390001098602</td><td>48.58026260878469</td><td>48.42496299743652</td><td>9.545562662829555</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1575.0590566276842</td><td>-3.966157316521702</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>124</td><td>nan</td><td>19.861563724094726</td><td>-1.4860336530214768</td><td>2458323.9027721277</td><td>2458323.90311935</td><td>ASASSN-V J205607.11+312038.1</td><td>0.126141972104651</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=664>\n",
       "  id       xcenter       ...           vsx_name                 mag_error      \n",
       "             pix         ...                                       adu         \n",
       "int32      float64       ...           bytes50                   float64       \n",
       "----- ------------------ ... ---------------------------- ---------------------\n",
       " 1127  3981.024900930935 ... ASASSN-V J205816.59+315327.0  0.003386311193159725\n",
       " 1071 3835.1134500924204 ... ASASSN-V J205811.97+315125.4  0.001252010506684385\n",
       " 1044 1888.2009196467923 ... ASASSN-V J205648.40+314708.0  0.014458509677046733\n",
       "  111  387.3365579918771 ... ASASSN-V J205607.11+312038.1  0.006956718274117619\n",
       " 1132  3978.345365708148 ... ASASSN-V J205816.59+315327.0 0.0033827794818791354\n",
       " 1075  3832.579448965366 ... ASASSN-V J205811.97+315125.4 0.0012460473535955033\n",
       " 1048 1885.9828031816776 ... ASASSN-V J205648.40+314708.0  0.014290265142871666\n",
       "  114 384.99479692707763 ... ASASSN-V J205607.11+312038.1  0.006855548127122464\n",
       " 1132  3977.865619318615 ... ASASSN-V J205816.59+315327.0  0.003396828328599868\n",
       "  ...                ... ...                          ...                   ...\n",
       " 1030 1697.0920941276536 ... ASASSN-V J205648.40+314708.0   0.21154653010203212\n",
       "   98 197.99324494658933 ... ASASSN-V J205607.11+312038.1   0.17929932984750513\n",
       " 1112 3787.9527167352517 ... ASASSN-V J205816.59+315327.0   0.04758048794015593\n",
       " 1056 3641.9638714304388 ... ASASSN-V J205811.97+315125.4  0.008957217967426507\n",
       " 1029 1696.0686240183486 ... ASASSN-V J205648.40+314708.0   0.27766849844528935\n",
       "   98 195.96100066832054 ... ASASSN-V J205607.11+312038.1   0.17597524853284996\n",
       " 1109  3785.021643076381 ... ASASSN-V J205816.59+315327.0   0.05596242871036663\n",
       " 1053 3638.9205418558713 ... ASASSN-V J205811.97+315125.4   0.00915560066834782\n",
       " 1026 1692.9911226337838 ... ASASSN-V J205648.40+314708.0   0.42033236398949014\n",
       "   97 193.00622531012922 ... ASASSN-V J205607.11+312038.1     0.126141972104651"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mags_grouped[mags_grouped['vsx_name'] != '']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The file name below is probably too generic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [],
   "source": [
    "#mags_grouped.write('calibrated_magnitudes_kelt_16_faint_13.fits')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The light curve we want to make is a scatter plot\n",
    "\n",
    "with the column 't_mid_tdb' on the horizontal axis and 'mag_cal' on the vertical axis."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: Cannot aggregate column 'date-obs' with type '|S19' [astropy.table.groups]\n",
      "WARNING: Cannot aggregate column 'filter' with type '|S1' [astropy.table.groups]\n",
      "WARNING: Cannot aggregate column 'vsx_name' with type '|S50' [astropy.table.groups]\n"
     ]
    }
   ],
   "source": [
    "regrouped = mags_grouped.group_by('file')\n",
    "agg_regrouped = regrouped.groups.aggregate(np.mean)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3wAAAEyCAYAAACh2dIXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X2wbWddH/Dvj1yCI6ES3q7IpSSQqBPtjCZXoGO1J+UtONZYASfR0UxLJnUGOrWOU0OtQBE7QLG0jtROCplJHSSgljGlUOTttjOdBpOLTCVC5BITuZKBQm7VK0Lenv5x9sXDZu979tl77bP3fs7nM3Pnnr3OOns/az1rPWt91/Ostau1FgAAAPrziFUXAAAAgOUQ+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdOrTqAszjCU94QrvgggtWXYyv85d/+Zd59KMfvepiMAB12Q912Rf12Q912Q912Q91uVmOHz/+hdbaE3ebbyMD3wUXXJDbb7991cX4OseOHcvW1taqi8EA1GU/1GVf1Gc/1GU/1GU/1OVmqap7ZpnPkE4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANCpQQJfVV1RVXdW1Ymqun7C7x9VVe8Y/f4jVXXB2O//ZlWdrqqfHaI8AAAADBD4quqcJG9O8sIklyS5uqouGZvtpUlOtdYuSvKmJK8f+/2bkrx30bIAAADw14bo4XtmkhOttbtaa/cnuTnJlWPzXJnkptHPv5XkOVVVSVJVP5zkriR3DFAWAAAARqq1ttgbVL04yRWttWtHr38iybNaay/fMc/HR/OcHL3+dJJnJfmrJB9I8rwkP5vkdGvtjVM+57ok1yXJ4cOHL7v55psXKvcynD59Ouedd96qi8EA1GU/1GVf1Gc/1GU/1GU/1OVmufzyy4+31o7uNt+hAT6rJkwbT5HT5vlXSd7UWjs96vCbqrV2Q5IbkuTo0aNta2tr7yVdsmPHjmUdy8Xeqct+qMu+qM9+qMt+qMt+qMs+DRH4TiZ56o7XR5J8dso8J6vqUJJvSnJftnv5XlxVb0jy2CQPV9WXW2u/OkC5AAAADrQhAt9tSS6uqguT/GmSq5L82Ng8tyS5Jsn/TvLiJB9q22NJv+/MDFX16mwP6RT2AAAABrBw4GutPVhVL0/yviTnJLmxtXZHVb0mye2ttVuSvDXJr1fViWz37F216OcCAABwdkP08KW19p4k7xmb9sodP385yUt2eY9XD1EWAAAAtg3yxesAAACsH4EPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdGiTwVdUVVXVnVZ2oqusn/P5RVfWO0e8/UlUXjKY/r6qOV9UfjP7/e0OUBwAAgAECX1Wdk+TNSV6Y5JIkV1fVJWOzvTTJqdbaRUnelOT1o+lfSPL3W2t/K8k1SX590fIAAACwbYgevmcmOdFau6u1dn+Sm5NcOTbPlUluGv38W0meU1XVWvv91tpnR9PvSPINVfWoAcoEAABw4A0R+J6S5DM7Xp8cTZs4T2vtwSR/luTxY/O8KMnvt9a+MkCZAAAADrxDA7xHTZjW9jJPVX1Htod5Pn/qh1Rdl+S6JDl8+HCOHTu254Iu2+nTp9eyXOyduuyHuuyL+uyHuuyHuuyHuuzTEIHvZJKn7nh9JMlnp8xzsqoOJfmmJPclSVUdSfKuJD/ZWvv0tA9prd2Q5IYkOXr0aNva2hqg6MM6duxY1rFc7J267Ie67Iv67Ie67Ie67Ie67NMQQzpvS3JxVV1YVecmuSrJLWPz3JLth7IkyYuTfKi11qrqsUn+W5JXtNb+1wBlAQAAYGThwDe6J+/lSd6X5BNJ3tlau6OqXlNVPzSa7a1JHl9VJ5L8TJIzX93w8iQXJfmFqvrY6N+TFi0TAAAAwwzpTGvtPUneMzbtlTt+/nKSl0z4u9cmee0QZQAAAOBrDfLF6wAAAKwfgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADo1SOCrqiuq6s6qOlFV10/4/aOq6h2j33+kqi7Y8btXjKbfWVUvGKI8AAAADBD4quqcJG9O8sIklyS5uqouGZvtpUlOtdYuSvKmJK8f/e0lSa5K8h1JrkjyH0bvBwAAwIKG6OF7ZpITrbW7Wmv3J7k5yZVj81yZ5KbRz7+V5DlVVaPpN7fWvtJa++MkJ0bvBwAAwIIODfAeT0nymR2vTyZ51rR5WmsPVtWfJXn8aPqtY3/7lEkfUlXXJbkuSQ4fPpxjx44NUPRhnT59ei3Lxd6py36oy76oz36oy36oy36oyz4NEfhqwrQ24zyz/O32xNZuSHJDkhw9erRtbW3toYj749ixY1nHcrF36rIf6rIv6rMf6rIf6rIf6rJPQwzpPJnkqTteH0ny2WnzVNWhJN+U5L4Z/xYAAIA5DBH4bktycVVdWFXnZvshLLeMzXNLkmtGP784yYdaa200/arRUzwvTHJxkt8boEwAAAAH3sJDOkf35L08yfuSnJPkxtbaHVX1miS3t9ZuSfLWJL9eVSey3bN31ehv76iqdyb5wyQPJnlZa+2hRcsEAADAMPfwpbX2niTvGZv2yh0/fznJS6b87S8l+aUhygEAAMBfG+SL1wEAAFg/Ah8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvhgBsfvOZU3f/hEjt9zatVFAQCWwLGeXh1adQFg3R2/51R+/C235v4HH865hx6Rt1377Fz2tPNXXSwA2FjH7zmVW+/6Yp799MevxTHVsZ6e6eGDXdx61xdz/4MP5+GWPPDgw7n1ri+uukisAVeCAeZzJlz98u/emR9/y61r0Y461tMzPXywi2c//fE599Aj8sCDD+eRhx6RZz/98asuEiu2yVeC1+2qOnDwTApX87ZHQ7VpjvX0TOCDXVz2tPPztmuf7SSZrxryZGU/7UdQFSiB3QwVroZs0xzr6ZnABzO47Gnna/z5qk29ErzsoLrJPZ/A/hkqXA3dpjnW0yuBD2CPNvVK8LKD6qb2fAL7b4hwtakX32C/CXwAc9jEK8HLDqpOvuZnKCzs3aZefIP9JvDRLSdQ8+tl3R2/51Te/en785gLT230cgxpmUHVydd8DIXdHL20jT3ZxItvsN8EPjbSbgfdaSdQDta76+Xk88xyfOWBh/Puu2/d2OU4m1m25/3e5sdPvnra55a1LIbCboZe2sZN0lP7Aask8HWs14ZyloPutO/TOWgH63m2gV5OPs8sR8tmL8c0s+wHs56gLqutWPYJ8n62cctcFkNhN0MvbeM627lPJwfvmA3LIvB1qucrkbMcdCedQB20g/W820AvJ59nluP+BzZ7OaaZZXueZZ5l9oYvc5/b7zZu2rIM0ctqKOxm6KVtXFfj+/SLLj1yoI7ZsEwCX6fW8UtNhzLLQXfaCdRBOljPuw30cvJ5Zjne/oHbcvVzv2djl2OaSfvB+L46y76yzN7woU+Qd96Tud8XcKat76F6Wee9D2nd2udJNqGMs+ilbVxX4/t0y/zH7F62ORiKwLckq25s1vFLTYcy60F3/ATqoB2sF9kGerkJ/rKnnZ+/eMa5uy7LqvfXeYxvz8nkkLbbNr/M3vAh97nxezJf+YPfsdQLOOPbxKRlefOHTwzSy7pIGdetfR63CWXci1W3jZvYVs1qvC160aVH8qJLj+x5eXvb5mAIAt8SrENjs65fajqUeQ+6qz5YL9MsJ6hDvXdPVn2f2yJ2bs/Twsdu2/yye8OH6rkavyfz1Jfun1juIepp2jYxviyzXFQZspdz2joZ751dp+10kWPIOu5zq7QO5xbLNK0t2usyrut5C6ySwLcEqzjATfq7IcLNtJOVXg7EQ67vVZZp1hPUecvT00nGLCfN8w7LW6Uhe3RX3Rs+aX1PuidzvNxD1dOsbfgs62modXm2dXKmzs//xnPXbjudd7vchH1uvx2EILPM85ZN0Mu5FetH4FuC/T7ALXJgnOdhAkMeiFfZuK1ifS/rvZd5ItDTScYsJ82z3ue2370U+/3gj/3sDZ8lhL/s8ot2vSdzqG11L234LOtpiF7Os62Ts82zin11vD7n2S5nXZZ1OEHerzIIMrNZ9QWrebnIwTItFPiq6nFJ3pHkgiR3J/nR1tqpCfNdk+Rfjl6+trV2U1V9Y5LfTPKMJA8l+a+ttesXKc+62Etjs9sBfZkn+/M+TGCok4pVD6Vb9vrez69EmPVEYJ4yraKXd1nfLzfLSfOk91p1L8WyH/yxSnsJ4Zc97ez3ZM7yIJtZLHLCOGTA37lept2zOF7nqw4EZxttsBez7HN72b829atHdlqXILPXdbmKIDNkW7hfx7p1uWBDnxbt4bs+yQdba6+rqutHr39u5wyjUPiqJEeTtCTHq+qWJF9J8sbW2oer6twkH6yqF7bW3rtgmdbCLI3NrAf0SX+3s/GZ92R03vs/hjqpWvVQurMFmTNPApz35H/eE5F563KWE4F51+Wye3nnKefQXzmx2/667F6K/XqfdTRvCJ9kvJ6S+Z82Os8J45D7xfh6mXbP4niZl3Vf47zlnndbnWWf28vFt2W1V/u9bw7RWzzkhYhZ1uUmt1+LjASa5aLlPOeAMI9FA9+VSbZGP9+U5FjGAl+SFyR5f2vtviSpqvcnuaK19vYkH06S1tr9VfXRJEcWLM9GmeeAPq3xmefkaN77P4Y6qZql52iZB4qzBZkzTwKctCzLPhGZ9wrubicCi6zLIXt5dzsQzvLe837+IlfIl9VLMclQF3XW1SwXOeY9sd35d7M8RXNIQ/b+T1ovs6yT8Xn2u3dlyG110vLOc4FsmceRZe+by3wI0TzmWZeb3H7Ns7yzfrfpPOeAMK9qrc3/x1X/r7X22B2vT7XWzh+b52eTfENr7bWj17+Q5K9aa2/cMc9jk3w0yXNba3dN+azrklyXJIcPH77s5ptvnrvcy3L69Omcd955U39/4tRD+eR9D+XbH3dOLjr/nJw49VDecNuX8+DDyaFHJP/8e74hF51/zlk/492fvj+//akH0pI8IsmPXPzI/OAzzp27zDvL9Mn7HprrvaeVaXx5d/v8nevkgYeTRz4i+bFvPze/8cn7Z1pHs3zevMsyi52fn2Ri3Y6Xcej63K18e93e9vJeSWaq7531O6kMs5RzL8uy2365bJO28bOtp2nraIjte2jzlGnS8iW7bztn7KU+h9zmh/q8WfaBnfOusk2b16zlHq/LefaNZPZ2Z6/bwazLsax9cy/bytkMuQ1MW5e71eU6tl+zGOo87dsfd87EbXc/26dZrfqYyd5cfvnlx1trR3ebb9cevqr6QJJvnvCrn5+xLDVh2ldTZlUdSvL2JL8yLewlSWvthiQ3JMnRo0fb1tbWjB+/f44dO5Zp5Tp+z6m88YNnrvg8lLdd++xcu3V+vvvSvV29e8yFp/Luu2/96pWyvXyh9KQrhVtjv5/nvSeVKcnXLe+k99oae33Hh0/kwXZnWpKHWvLEp16Yt//dx880NGKWz5t1We5/4OGc+8jZ18Gkz3/7pV87PHbSPFdfmLnrc6+2kj1vb7O+V7Kzl/ehqVc0x+v3K499Wra2LtpzOfeyLGfbL/fDzk+eZTudto62sl7m3ecmLd/LLr9o1787Yy/1uZXhtvmhPm+WfWDn+y1qkWPGvLZmnG9nXc67b8yy7Wxltu1gZ3uVzHYM28vy7tVetpWzGXIb2MrkdblbXV67tZm9VVvZexsyaX3fetcXJ267+9k+zWrVx0yWY9fA11p77rTfVdXnqurJrbV7q+rJST4/YbaT+dr28Ei2h36ecUOST7XW/t1MJd5Q04YF7HXY0rxD0mYZ0jHve0/6u2lDqXYbnjLvMKah7xs525MAZ/38l11+0a7DIOe9V2le8w6T2+29JtV38vVDfWcd2jPP0LVk9U/sG2K46qYMf1r2A4aGst/byW7b7iqWfxOGiS1735g01PVs9yW/6NIjK7/vbKhtZehtYJm3D6yjoc7Thhy6Dnu16D18tyS5JsnrRv//zoR53pfkX1fVmS36+UlekSRV9dok35Tk2gXLsfYWabgnfaH2XhuIWRvgIe6bSaY/2GVZoXPo+0bO9iTAeT9/6HuV1smkZZs14K7yYQJDmuXzZ9lONuUEfZkPGFqmWe+vWZZVLP+8bcykULSscu/nvjFpGxhvr1pW/8TTIbeV/TzOrOtFq/3+aoidn7Hqdg8WDXyvS/LOqnppkj9J8pIkqaqjSX6qtXZta+2+qvrFJLeN/uY1o2lHsj0s9JNJPlpVSfKrrbW3LFimtbTMnrlZrMNV5VkfoDDPgWnVjeksn79IGVf9FQi7mfeK5qofJjCkWT5/1m1gEy4CLLI9r3L5JtVTMv/TPOexCfU76QmCr3n3HUtbR/u5b0zaBsaPkS+69EhedOmRlZ+gr3pbmfVpkzufbL3q4/G0Mq76O+5WXZccbAsFvtbaF5M8Z8L027Oj1661dmOSG8fmOZnJ9/dtpN0e5Z/Mt7MPPVRxlVeVlx06V92YzvL585Rx1V+BMKt5rmgOGdJWfVV5yOGqm2ITl2XW3uhNW6692uvw4/d+/N6lr6P92p6m3Towqb3qeTvYbRvYy1fkjD/Zet3ahoO4j8NOi/bwkekN3hCW/Yjr/bSOV/02wbQD1RA9c8s+CO62zQ29fa97L+8mW/X9kUPZS290r+YZfvzC73xybrv7vi7W0dnC3SZv23sxyzYwy/HhzDwt6x2kVn1BEFZN4BvAMhu83k4iD9IBdSjz3g8573vvp6G371VvX6v+/GVZh+FQQ5qnN7on8w4//rZvfszXraNNvRDQ6746q6EeknNmnvsfWO8gddD2cRgn8A1g2Q3eQT8wHXSTDlRDfaH0OhwEbd/rf9J8EIZDHaTtcN7hx+OvF7kQMMttEJtq0v68bvv4UA/JOTPPXp9svQoHaR+HcQLfADapwWMzLfN+SAfB1dqE3rNV9wQzrKEu9Mx7IWCZt0Gs2qT9OdnfhwLNYtZtYNZ70/f6ZOverFugh3EC30A0eOyndeiZYxib0Htme+vPEBd65r0QsOz7vlZ58j3tKbDL3MfnXV4X+4axCRftQOCDDeVg3YdN6T2zvTFu3gsBy7wNYtUn39P250n3Yffw/aOrsG69aZtw0Q4EPoAV0nvGJpvnQsAyb4NY9cn3tP1557RkuCGeq17e/bYOAXc8cG7KRTsONoEPYMX0nnHQLOs2iHU4+Z60P++cNtRDt5L1WN79tOqAOy1wumjHuhP4AMas25AhYDabcPK9H98/2msbtuqAOy1wumjHuhP4AHZYhyFDwPzW/eR76JA25NdlrLtVB/pVB06Yl8AHsMOqhwzBKvTaI7SulhnSem/DVhnoVx04YV4CH8AOruBy0PTcI7Qphgxp2rDlWvceZJhE4APW3n72PriCy0HTe4/QJtiP+/qAg0vgA9baKnofXMHlINEjtHpDhzRtGLCTwAesNb0PsFx6hNbDqkPaQbuP86AtLwebwAesNb0PsHyrDhus1kG7j/OgLS8IfMBa0/sAsFwHbSTFQVteEPiAtaf3AWB5DtpIioO2vCDwAQAcYAdtJMVBW14Q+AAADriDNpLioC0vB9sjVl0AAAAAlkPgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADolMAHAADQKYEPAACgUwIfAABApwQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdGqhwFdVj6uq91fVp0b/nz9lvmtG83yqqq6Z8Ptbqurji5QFAACAr7VoD9/1ST7YWrs4yQdHr79GVT0uyauSPCvJM5O8amcwrKofSXJ6wXIAAAAwZtHAd2WSm0Y/35TkhyfM84Ik72+t3ddaO5Xk/UmuSJKqOi/JzyR57YLlAAAAYMyige9wa+3eJBn9/6QJ8zwlyWd2vD45mpYkv5jkl5N8acFyAAAAMObQbjNU1QeSfPOEX/38jJ9RE6a1qvquJBe11v5ZVV0wQzmuS3Jdkhw+fDjHjh2b8eP3z+nTp9eyXOyduuyHuuyL+uyHuuyHuuyHuuzTroGvtfbcab+rqs9V1ZNba/dW1ZOTfH7CbCeTbO14fSTJsSR/O8llVXX3qBxPqqpjrbWtTNBauyHJDUly9OjRtrU1cbaVOnbsWNaxXOyduuyHuuyL+uyHuuyHuuyHuuzTokM6b0ly5qmb1yT5nQnzvC/J86vq/NHDWp6f5H2ttV9rrX1La+2CJH8nyR9NC3sAAADs3aKB73VJnldVn0ryvNHrVNXRqnpLkrTW7sv2vXq3jf69ZjQNAACAJdp1SOfZtNa+mOQ5E6bfnuTaHa9vTHLjWd7n7iTfuUhZAAAA+FqL9vABAACwpgQ+AACATgl8AAAAnRL4AAAAOiXwAQAAdErgAwAA6JTABwAA0CmBDwAAoFMCHwAAQKcEPgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBTAh8AAECnBD4AAIBOCXwAAACdEvgAAAA6JfABAAB0SuADAADoVLXWVl2GPauq/5vknlWXY4InJPnCqgvBINRlP9RlX9RnP9RlP9RlP9TlZnlaa+2Ju820kYFvXVXV7a21o6suB4tTl/1Ql31Rn/1Ql/1Ql/1Ql30ypBMAAKBTAh8AAECnBL5h3bDqAjAYddkPddkX9dkPddkPddkPddkh9/ABAAB0Sg8fAABApwQ+AACATgl8A6iqK6rqzqo6UVXXr7o8zK6qnlpVH66qT1TVHVX1T0fTX11Vf1pVHxv9+4FVl5XZVNXdVfUHo3q7fTTtcVX1/qr61Oj/81ddTs6uqr5tx/73sar686r6afvm5qiqG6vq81X18R3TJu6Lte1XRsfR/1NVl66u5IybUpf/pqo+Oaqvd1XVY0fTL6iqv9qxj/7H1ZWccVPqcmq7WlWvGO2Xd1bVC1ZTahblHr4FVdU5Sf4oyfOSnExyW5KrW2t/uNKCMZOqenKSJ7fWPlpVj0lyPMkPJ/nRJKdba29caQHZs6q6O8nR1toXdkx7Q5L7WmuvG12UOb+19nOrKiN7M2pn/zTJs5L8w9g3N0JVfX+S00n+c2vtO0fTJu6LoxPMf5LkB7Jdz/++tfasVZWdrzWlLp+f5EOttQer6vVJMqrLC5K8+8x8rJcpdfnqTGhXq+qSJG9P8swk35LkA0m+tbX20L4WmoXp4VvcM5OcaK3d1Vq7P8nNSa5ccZmYUWvt3tbaR0c//0WSTyR5ympLxRJcmeSm0c83ZTvUszmek+TTrbV7Vl0QZtda+59J7hubPG1fvDLbJ6CttXZrkseOLsixBibVZWvtd1trD45e3prkyL4XjD2bsl9Oc2WSm1trX2mt/XGSE9k+72XDCHyLe0qSz+x4fTICw0YaXZX87iQfGU16+Wioyo2GAG6UluR3q+p4VV03mna4tXZvsh3ykzxpZaVjHldl+yrzGfbNzTVtX3Qs3Wz/KMl7d7y+sKp+v6r+R1V936oKxZ5Malftl50Q+BZXE6YZJ7thquq8JL+d5Kdba3+e5NeSPCPJdyW5N8kvr7B47M33ttYuTfLCJC8bDV9hQ1XVuUl+KMlvjibZN/vkWLqhqurnkzyY5G2jSfcm+Zutte9O8jNJfqOq/saqysdMprWr9stOCHyLO5nkqTteH0ny2RWVhTlU1SOzHfbe1lr7L0nSWvtca+2h1trDSf5TDGHYGK21z47+/3ySd2W77j53ZnjY6P/Pr66E7NELk3y0tfa5xL7ZgWn7omPpBqqqa5L8YJIfb6OHQoyG/31x9PPxJJ9O8q2rKyW7OUu7ar/shMC3uNuSXFxVF46uRF+V5JYVl4kZVVUleWuST7TW/u2O6TvvHfkHST4+/resn6p69OjhO6mqRyd5frbr7pYk14xmuybJ76ymhMzh6uwYzmnf3HjT9sVbkvzk6Gmdz07yZ2eGfrKequqKJD+X5Idaa1/aMf2JowctpaqenuTiJHetppTM4izt6i1JrqqqR1XVhdmuy9/b7/KxuENrQHIqAAAA+klEQVSrLsCmGz2d6uVJ3pfknCQ3ttbuWHGxmN33JvmJJH9QVR8bTfsXSa6uqu/K9tCFu5P849UUjz06nORd2zk+h5L8Rmvtv1fVbUneWVUvTfInSV6ywjIyo6r6xmw/AXnn/vcG++ZmqKq3J9lK8oSqOpnkVUlel8n74nuy/YTOE0m+lO2nsbImptTlK5I8Ksn7R23ura21n0ry/UleU1UPJnkoyU+11mZ9SAhLNqUutya1q621O6rqnUn+MNvDdl/mCZ2bydcyAAAAdMqQTgAAgE4JfAAAAJ0S+AAAADol8AEAAHRK4AMAAOiUwAcAANApgQ8AAKBT/x+qsexbrF/3YQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15, 5))\n",
    "plt.plot(agg_regrouped['color_term'], '.')\n",
    "plt.ylim(-0.05, 0.05)\n",
    "plt.grid()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x1cf81d8bda0>]"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3QAAAE/CAYAAAAOkIE9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X+cZXV93/HXZxkWqwthXHWDLtllI6KIFdnRbMLDdhcCEqRKS5pIicEopUlsovVHg1Xzs2m1aU30oY+Hj+2qYNyyrYCFEK0S3AkPG9bCIAq4KLiyuEJDAktlJWEZ9tM/7hm8e/fcvWfO3PnxnXk9H495zL3nnu+533M/98w97/s950xkJpIkSZKk8iyb7w5IkiRJktox0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhRqZ7w7Uec5znpNr166d8+f94Q9/yLOe9aw5f161Y73KYr3KYa3KYr3KYr3KYr3KstjqNTEx8beZ+dxB8w0MdBFxPPBp4MeBA8DmzPxwRDwb+O/AWuA+4Bcyc29N+58AtgDHAwmcm5n3He45165dy6233jqoa0M3Pj7Oxo0b5/x51Y71Kov1Koe1Kov1Kov1Kov1Kstiq1dE7G4yX5NDLieBd2bmS4ANwFsj4mTgMuDGzDwRuLG6X+fTwB9V7V8FPNSkY5IkSZKkwxsY6DLzwcy8rbr9GLATeAHweuCKarYrgPN721bBbyQzb6ja78vMx4fUd0mSJEla0iIzm88csRa4CTgFuD8zj+16bG9mjvbMfz5wCbAfOAH4C+CyzHyqZtmXApcCrFq1av22bdumuy4ztm/fPlasWDHnz6t2rFdZrFc5rFVZrFdZrFdZrFdZFlu9Nm3aNJGZY4Pma3xRlIhYAVwNvD0zfxARTZqNAK8GXgHcT+ecuzcBn+idMTM3A5sBxsbGcj6Of11sx90udtarLNarHNaqLNarLNarLNarLEu1Xo3+bUFEHEknzG3NzGuqyX8dEcdVjx9H/blxe4CvZeauzJwE/idw2sy7LUmSJEkaGOiiMxT3CWBnZn6o66HrgIur2xcD19Y0vwUYjYipy22eAXyzfXclSZIkSVOajNCdDrwROCMibq9+zgU+AJwVEfcAZ1X3iYixiNgCUJ0r9y7gxoi4Awjgv87CekiSJEnSkjPwHLrM/AqdIFbnzJr5b6VzIZSp+zcA/7BtByVJkiRJ9RqdQydJWpgmdu/lY9vvZWL33vnuiiRJmgeNr3IpSZp/E7v3smPXw2xYtxKAi7bsYP/kAZaPLGPrJRtYv2Z0YLv1a0YPud89z1GPPsXGuVqhJaLu9ZYkaRgMdJI0ZMPaea8LYd0B7oLTVrN/8gAHEp6cPMCOXQ/XhrXedr993kv5/evvOigIwo/C4UjAK07bW0zwmM2w1CQIN2nTNHhPtz8A9+59iru231sbzufiNZEkzS8DnaQlqe1O6aB2w9p5r1vOjl0PHxTgElg+sownJw9w5MgyNqxb2ajdF+588JAgCDw9bTJ5OhwudHXrC/QdfZxuMBsUhHuXBYeOmva+/tN5bQctG+A/3fL3TOa3asP5dN6DbV6Tfss29EnS3DHQaVqHYvnhfDC/BW9mYvderv/Ofo4+YW5GfWYrdDUJD3U7773zNOl33XI2rFt5UIC74LTVXHDa6oOW/bHt9w5s93OnHMct9z1yUBCEH4XDI4Knw+FcvgfbhK7e1+nq2/ZwzW17+o4+Ng1mdbXsF4QHjZr2vv5Tr3eT9R+0bIAnD0BSH86bvgf7bRODXpO6cDrMEUlJ0mAGupaahJ6mO0Nt2w1rPaZzKNbhPuiXmpnstMzVKM9CMLUuTzx5gOvv2zHr69LktWsaugbtzNaFh96d99FnLm+0/dRti70hYP2a0adHfLr72b1+deGhrt1JP370IcuZmueoR3cD7UZ56l63NnVrGrp61zeoDzPTDWZ1tawLwk1GTfvVbboje3XLBjhyGTyV1Ibzw70Hm2wTg16TuuA/rC81/JyRpGYMdA01Oeylzc5I252YpqNog9pN91Csw33Qw+Ad4iZ9aqruvJHD1W3YOwP9vqketM4wOCDP5BCttoY1MtvvPdY9gjCb69JkZ7JJ6ILBO7N14eGtm1540M570x3l3vn2Pr6/b3g73Ot3uNDX3a5uOVPTxsf3tH4Ptj0kr81o2NSyu9cX4Orb9vQdfWwazOpq2S8IDxo1rXu9m4y+NRmRBfi3r3wGTxy7pjacN/17UhfUmrwmdTUZ5pcadZ+HTXikydIyzNr6PlGJDHQNND3spc3OSJudmKajaE3aTfdQrH4f9E0Oc5pO8G1Sk97zRpoe9tN25LR3Wr9vqgetc7/3z6DRmab9bLMj06TfTUJ7XZup12n/kwevy2x9E990Z7JJ6GqyM1sXHnp33pvsKPcbWWvzWrRtd7jXsV/derU9JK/NaNjUsnvXty7QTjeY9atl3f1Bo6Z1ZjKy17vsF44ewcaNLzxo2qD3YK9+zzXoNak7xHdYX2rUfR42+VvR9O+ZO+tza7aObILh1bbuvdPky2NpvhnoGmh62EubnZE2OzHTGUUb1K7pN9CDPuibHObUtE9NRiR37Hr4kPNGmoxyQbM//E13CHpfl7qdm951rnv/NB2dGdTP6Rwye7j3eF2/m4T2urA69R678i9u4cKffeWsfxPfJKzVhYBB23S/HfxBdWq6o9xvvvlS158mo2/9gmC3YY2G9fui43Cjj4e73/b1bxOgm46+DSOc91u3um1nuq9Jv3oP40uNus/DmZzH2mTEuc1VTDVYk8/Utkc29bvab10fBm0DdV9WX1Xz5bHvAS00BroGmn7wttkZqfvAHNSu6Sha03aDdnT6zTPdw5ya9qnJzv6GdSsPOW+kyTlIhxs9aBIG63ZCu1+Xfjs3Td4/bUZnZhL2u5+/Sb+bhPZ+X3asXzPKYz+5/On1afpNfK9BH+xTH7bDGp1osoPfpE7TGdVZSDsHvf3pt/301mVYIWCYoavN+s6mw70HZuv5upfd9NDYQa9Jk5rM5EuN3s/D3i/NmpzHerjPo25NwuJ8nN/c9gqtw9pO2i5rup+pbY9s6ve5M91TZfodVl930aEmR69Ic8lA10DTD962OyNt2jUdRWvSbiavS/cymu4QT/cDu9/IYu95I73t+o1y9fvD3+SQx2EGgzY7qoMuCNEk7Dc5BLGu3zA4tPcLq72a9nu6F42oC4JNX9smX2QM00ILb03U7SgfLlR3m27o66fE162f+VyXJttOU7P5pUbv/bod7iajvXDo3+E2Fz0a5uvWRJujGYZ5wa4mIbdJuyafqW2PbKr73Glzqkzdewfgs7fcf9CXx03eJ4vlb5TKYaBraFgfvG2XM+hDrum0ud5BbdOnpjv7veeN1O1sNt2J6P0D3S8MzmYwGDRPvw/p6Yb9fjskTfrdJqz2W9dB/W5z0Yi6kN309ddgdXWrO8y46Q6mNZk/Tbed2dTmPdD7HoRm57H2TmtyHmddWOz3uk1n1OyoR59iY8P1bXo0wzAuqlX3mjQJL9D+Ik/DOrJpqr79Xrcmp8r0e+/UXXRo0Puk7guDpu8TqQ0DnRacpiGlSbt+8/U+1iQM9ms7V5oGsSb9brsjN8zXZFC/m3wgN625hmfQFzB176e5HtXQYCVvO02+aBqk7j3Z5KJHda/bdA/VHAl4xWn1/5OzzVEYTUfD2lzQqEnIpWZa28/Upl/CDvrcaTKKB83eO71fHjf5UsEL82iuGei0ILUdWZxJuFjoOzfD+ka9hHWFub1ohNpr8n5aCKNBOtRi2XbarEe/92STsDjoy6dBh2pOJk8HoSbn5033NIW60bAmV33uF8LanC8/zM+ZYYzk9jt6ZFhfSLa5QJpfbGmYDHRSZaHv3Mz3B+Rca/qBrPk36P1UypcIWjpmckRHr+keqnlEUHsuc9ujMJqMhtUtGwZfubnu+aZz6P18butz+fxNj1zwiy3NFgOdVJD5/oCca0ttfRcza6mFZljvySajWN3zHPXobvY+vr/1OcGDnr/pKHm/w04HvSZNDp1c6vrVxC+2NFsMdJIkSTMwnUM1x8f3cPQJgw9vnM4Of9tRckeMZo/BV3PJQCdJkjREbQPWbO7w1x266YiRtDgY6CRJkubYQhitWQh9kDRzy+a7A5IkSZKkdgx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhRoY6CLi+IjYHhE7I+KuiHhbNf3ZEXFDRNxT/e77nykj4piI+H5EfHSYnZckSZKkpazJCN0k8M7MfAmwAXhrRJwMXAbcmJknAjdW9/v5A+AvZ9pZSZIkSdKPDAx0mflgZt5W3X4M2Am8AHg9cEU12xXA+XXtI2I9sAr40jA6LEmSJEnqiMxsPnPEWuAm4BTg/sw8tuuxvZk52jP/MuDLwBuBM4GxzPzXfZZ9KXApwKpVq9Zv27ZtWisyDPv27WPFihVz/rxqx3qVxXqVw1qVxXqVxXqVxXqVZbHVa9OmTROZOTZovpGmC4yIFcDVwNsz8wcR0aTZrwOfz8zvDZo/MzcDmwHGxsZy48aNTbs2NOPj48zH86od61UW61UOa1UW61UW61UW61WWpVqvRoEuIo6kE+a2ZuY11eS/jojjMvPBiDgOeKim6U8Dr46IXwdWAMsjYl9mHu58O0mSJElSAwMDXXSG1j4B7MzMD3U9dB1wMfCB6ve1vW0z86Ku5byJziGXhjlJkiRJGoImV7k8nc45cGdExO3Vz7l0gtxZEXEPcFZ1n4gYi4gts9ZjSZIkSRLQYIQuM78C9DsB7sya+W8FLqmZfjlw+fS6J0mSJEnqp8kInSRJkiRpATLQSZIkSVKhDHSSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhTLQSZIkSVKhDHSSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhTLQSZIkSVKhDHSSJEmSVKiBgS4ijo+I7RGxMyLuioi3VdOfHRE3RMQ91e/RmranRsTNVbtvRMQvzsZKSJIkSdJS1GSEbhJ4Z2a+BNgAvDUiTgYuA27MzBOBG6v7vR4HfjkzXwqcA/xJRBw7nK5LkiRJ0tI2MNBl5oOZeVt1+zFgJ/AC4PXAFdVsVwDn17T9dmbeU91+AHgIeO5wui5JkiRJS1tkZvOZI9YCNwGnAPdn5rFdj+3NzEMOu+x6/FV0gt9LM/NAzeOXApcCrFq1av22bdsa92tY9u3bx4oVK+b8edWO9SqL9SqHtSqL9SqL9SqL9SrLYqvXpk2bJjJzbNB8I00XGBErgKuBt2fmDyKicWci4jjgT4GL68IcQGZuBjYDjI2N5caNGxsvf1jGx8eZj+dVO9arLNarHNaqLNarLNarLNarLEu1Xo2uchkRR9IJc1sz85pq8l9XQW0qsD3Up+0xwJ8D78vMHTPvsiRJkiQJml3lMoBPADsz80NdD10HXFzdvhi4tqbtcuBzwKcz87Mz764kSZIkaUqTEbrTgTcCZ0TE7dXPucAHgLMi4h7grOo+ETEWEVuqtr8A/CPgTV1tTx3+akiSJEnS0jPwHLrM/ArQ74S5M2vmvxW4pLr9GeAzM+mgJEmSJKleo3PoJEmSJEkLj4FOkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEIZ6CRJkiSpUAY6SZIkSSqUgU6SJEmaYxO79/Kx7fcysXvvfHdFhRuZ7w5IkiRJS8nE7r1ctGUH+ycPsHxkGVsv2cD6NaPz3S0VyhE6SZIkaQ7t2PUw+ycPcCDhyckD7Nj18Hx3SQUz0EmSJElzaMO6lSwfWcYRAUeOLGPDupXz3SUVzEMuJUmSpDm0fs0oWy/ZwI5dD7Nh3UoPt9SMGOgkSZKkObZ+zahBTkPhIZeSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhRoY6CLi+IjYHhE7I+KuiHhbNf3ZEXFDRNxT/R7t0/7iap57IuLiYa+AJEmSJC1VTUboJoF3ZuZLgA3AWyPiZOAy4MbMPBG4sbp/kIh4NvA7wE8BrwJ+p1/wkyRJkiRNz8BAl5kPZuZt1e3HgJ3AC4DXA1dUs10BnF/T/DXADZn5SGbuBW4AzhlGxyVJkiRpqZvWOXQRsRZ4BfBVYFVmPgid0Ac8r6bJC4Dvdd3fU02TJEmSJM1QZGazGSNWAH8J/GFmXhMRj2bmsV2P783M0Z427waOysx/X91/P/B4Zv6XmuVfClwKsGrVqvXbtm1ru06t7du3jxUrVsz586od61UW61UOa1UW61UW61UW61WWxVavTZs2TWTm2KD5RposLCKOBK4GtmbmNdXkv46I4zLzwYg4DniopukeYGPX/dXAeN1zZOZmYDPA2NhYbty4sW62WTU+Ps58PK/asV5lsV7lsFZlsV5lsV5lsV5lWar1anKVywA+AezMzA91PXQdMHXVyouBa2uafxE4OyJGq4uhnF1NkyRJkiTNUJNz6E4H3gicERG3Vz/nAh8AzoqIe4CzqvtExFhEbAHIzEeAPwBuqX5+v5omSZIkSZqhgYdcZuZXgOjz8Jk1898KXNJ1/5PAJ9t2UJIkSZJUb1pXuZQkSZIkLRwGOkmSJEkqlIFOkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEIZ6CRJkiSpUAY6SZIkSSqUgU6SJEmSCmWgkyRJkqRCGegkSZIkqVAGOkmSJEkqlIFOkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRpCZrYvZePbb+Xid1757srQ9d03RbDazAy3x2QJEmSNLcmdu/loi072D95gOUjy9h6yQYAdux6mA3rVrJ+zeg897C/id17D9vPunWbyXwLnYFOkiRJWmJ27HqY/ZMHOJDw5OQBrr5tD9fctmfBh5smIax33Xbsepj1a0YPCYL95iuNgU6SJElaYjasW8nykWU8OXmAI0eWEVBEuGkSwnrXbcO6lbVBsG6+EhnoJEmSpIL1jjwNOiQRYP2aUbZesuHp+QCuvm3PnIabun4O6nu/EDaxey/Xf2c/R5+w95B1W79mlI9tv/eQIPjWTS88ZL4SGegkSZKkBWq654v99nkv5fevv6vRoZPr14we9FjTcNMkMA5qB9SewzfocMq6sDb1Gjzx5AGuv2/H0+262/YLgr3zlchAJ0mSJC1Abc4X+8KdD7Y+dLIu3NSN/rW5kEhvuwtOW31IP6HZYZ+9/Zx6DZL+7eqC4GJhoJMkSZIWoDbni/3cKcdxy32PDOXQybrw1vZCIr3tEmpHzNqc0zb1Gux/8vDtFsNoXB0DnSRJkrQANbloR93I00k/fvRQRqLqwlvbC4n0trvgtNVccNrqQ/rZZhRt6jW48i9u4cKffeWiDG2HY6CTJEmSFqCmhwn2jjwNaySqLry1PXSxX7u6QyPb9H39mlEe+8nlSy7MQYNAFxGfBM4DHsrMU6ppLwc+DqwA7gMuyswf1LT9N8AlQAJ3AL+SmX8/tN5LkiRJi0TdxUbm8zDBw4Ww7j41vUhKm3VpewGWpaTJCN3lwEeBT3dN2wK8KzP/MiLeDLwbeH93o4h4AfCbwMmZ+XcR8T+AN1TLkyRJklRpe7GR2TYohM1mvxfqa7LQLBs0Q2beBDzSM/kk4Kbq9g3ABX2ajwD/ICJGgGcCD7TspyRJkrRo1Z2vVoLZ7Hepr8lca3sO3Z3A64BrgX8OHN87Q2Z+PyL+M3A/8HfAlzLzS/0WGBGXApcCrFq1ivHx8ZZda2/fvn3z8rxqx3qVxXqVw1qVxXqVxXqVZS7rddSjTzESMJlwRMBRj+5mfHzPnDz3TMxmv6e77KW6fUVmDp4pYi1wfdc5dC8GPgKsBK4DfjMzV/a0GQWuBn4ReBT4LHBVZn5m0PONjY3lrbfeOq0VGYbx8XE2btw458+rdqxXWaxXOaxVWaxXWaxXWea6XqWeLzab/Z7Oshfb9hURE5k5Nmi+ViN0mXk3cHb1RC8CXlsz288C383Mv6nmuwb4GWBgoJMkSZKWmlL/T9ps9rvU12QuDTyHrk5EPK/6vQx4H50rXva6H9gQEc+MiADOBHa27agkSZIk6WADA11EXAncDJwUEXsi4i3AhRHxbeBuOhc6+VQ17/Mj4vMAmflV4CrgNjr/smAZsHlW1kKSJEmSlqCBh1xm5oV9HvpwzbwPAOd23f8d4Hda906SJEmS1FerQy4lSZIkSfPPQCdJkiRJhTLQSZIkSVKhDHSSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhTLQSZIkSVKhDHSSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhTLQSZIkSVKhBga6iPhkRDwUEXd2TXt5RNwcEXdExJ9FxDF92h4bEVdFxN0RsTMifnqYnZckSZKkpazJCN3lwDk907YAl2Xmy4DPAe/u0/bDwP/KzBcDLwd2tuynJEmSJKnHwECXmTcBj/RMPgm4qbp9A3BBb7tq1O4fAZ+olrM/Mx+dUW8lSZIkSU+LzBw8U8Ra4PrMPKW6/1fABzPz2oh4B/B7mXl0T5tTgc3AN+mMzk0Ab8vMH/Z5jkuBSwFWrVq1ftu2bW3XqbV9+/axYsWKOX9etWO9ymK9ymGtymK9ymK9ymK9yrLY6rVp06aJzBwbNF/bQPdi4CPASuA64Dczc2VPmzFgB3B6Zn41Ij4M/CAz3z/o+cbGxvLWW28d2K9hGx8fZ+PGjXP+vGrHepXFepXDWpXFepXFepXFepVlsdUrIhoFulZXuczMuzPz7MxcD1wJfKdmtj3Ansz8anX/KuC0Ns8nSZIkSTpUq0AXEc+rfi8D3gd8vHeezPy/wPci4qRq0pl0Dr+UJEmSJA1Bk39bcCVwM3BSROyJiLcAF0bEt4G7gQeAT1XzPj8iPt/V/DeArRHxDeBU4D8MewUkSZIkaakaGTRDZl7Y56EP18z7AHBu1/3bgYHHfUqSJEmSpq/VIZeSJEmSpPlnoJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEIZ6CRJkiSpUAY6SZIkSSqUgU6SJEmSCmWgkyRJkqRCGegkSZIkqVAGOkmSJEkqlIFOkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEIZ6CRJkiSpUAY6SZIkSSqUgU6SJEmSCmWgkyRJkqRCGegkSZIkqVAGOkmSJEkqlIFOkiRJkgploJMkSZKkQhnoJEmSJKlQAwNdRHwyIh6KiDu7pr08Im6OiDsi4s8i4pjDtD8iIr4WEdcPq9OSJEmSpGYjdJcD5/RM2wJclpkvAz4HvPsw7d8G7GzVO0mSJElSXwMDXWbeBDzSM/kk4Kbq9g3ABXVtI2I18Fo6AVCSJEmSNESRmYNnilgLXJ+Zp1T3/wr4YGZeGxHvAH4vM4+uaXcV8B+Bo4F3ZeZ5h3mOS4FLAVatWrV+27Zt01+bGdq3bx8rVqyY8+dVO9arLNarHNaqLNarLNarLNarLIutXps2bZrIzLFB8420XP6bgY9ExG8D1wH7e2eIiPOAhzJzIiI2DlpgZm4GNgOMjY3lxo0Dmwzd+Pg48/G8asd6lcV6lcNalcV6lcV6lcV6lWWp1qtVoMvMu4GzASLiRXQOq+x1OvC6iDgXeAZwTER8JjN/qW1nJUmSJEk/0urfFkTE86rfy4D3AR/vnScz35OZqzNzLfAG4MuGOUmSJEkanib/tuBK4GbgpIjYExFvAS6MiG8DdwMPAJ+q5n1+RHx+NjssSZIkSeoYeMhlZl7Y56EP18z7AHBuzfRxYHyafZMkSZIkHUarQy4lSZIkSfPPQCdJkiRJhTLQSZIkSVKhDHSSJEmSVCgDnSRJkiQVykAnSZIkSYUy0EmSJElSoQx0kiRJklQoA50kSZIkFcpAJ0mSJBViYvdePrb9XiZ2753vrmiBGJnvDkiSJEkabGL3Xi7asoP9kwdYPrKMrZdsYP2a0fnuluaZI3SSJElSAXbsepj9kwc4kPDk5AF27Hp4vrukBcBAJ0mSJBVgw7qVLB9ZxhEBR44sY8O6lfPdJS0AHnIpSZIkFWD9mlG2XrKBHbseZsO6lR5uKcBAJ0mSJBVj/ZpRg5wO4iGXkiRJklQoA50kSZIkFcpAJ0mSJEmFMtBJkiRJUqEMdJIkSZJUKAOdJEmSJBXKQCdJkiRJhTLQSZIkSVKhDHSSJEmSVCgDnSRJkqShmti9l49tv5eJ3XvnuyuL3sh8d0CSJEnS4jGxey8XbdnB/skDLB9ZxtZLNrB+zeh8d2vRcoROkiRJ0tDs2PUw+ycPcCDhyckD7Nj18Hx3aVEbGOgi4pMR8VBE3Nk17eURcXNE3BERfxYRx9S0Oz4itkfEzoi4KyLeNuzOS5IkSVpYNqxbyfKRZRwRcOTIMjasWznfXVrUmozQXQ6c0zNtC3BZZr4M+Bzw7pp2k8A7M/MlwAbgrRFx8gz6KkmSJGmBW79mlK2XbOAdZ5/k4ZZzYOA5dJl5U0Ss7Zl8EnBTdfsG4IvA+3vaPQg8WN1+LCJ2Ai8AvjmzLkuSJElayNavGTXIzZHIzMEzdQLd9Zl5SnX/r4APZua1EfEO4Pcy8+gB7W8CTsnMH/SZ51Ip5sJZAAAMlElEQVTgUoBVq1at37Zt27RWZBj27dvHihUr5vx51Y71Kov1Koe1Kov1Kov1Kov1Kstiq9emTZsmMnNs0Hxtr3L5ZuAjEfHbwHXA/n4zRsQK4Grg7f3CHEBmbgY2A4yNjeXGjRtbdq298fFx5uN51Y71Kov1Koe1Kov1Kov1Kov1KstSrVerQJeZdwNnA0TEi4DX1s0XEUfSCXNbM/Oatp2UJEmS1N7E7r3s2PUwG9at9FDIRaZVoIuI52XmQxGxDHgf8PGaeQL4BLAzMz80s25KkiRJasP/C7e4Nfm3BVcCNwMnRcSeiHgLcGFEfBu4G3gA+FQ17/Mj4vNV09OBNwJnRMTt1c+5s7IWkiRJkmr5f+EWtyZXubywz0Mfrpn3AeDc6vZXgJhR7yRJkiTNyNT/hXty8oD/F24RantRFEmSJEkFmPq/cJ5DtzgZ6CRJkqRFzv8Lt3gNPIdOkiRJkrQwGegkSZIkqVAGOkmSJEkqlIFOkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEJFZs53Hw4REX8D7J6Hp34O8Lfz8Lxqx3qVxXqVw1qVxXqVxXqVxXqVZbHVa01mPnfQTAsy0M2XiLg1M8fmux9qxnqVxXqVw1qVxXqVxXqVxXqVZanWy0MuJUmSJKlQBjpJkiRJKpSB7mCb57sDmhbrVRbrVQ5rVRbrVRbrVRbrVZYlWS/PoZMkSZKkQjlCJ0mSJEmFMtBJkiRJUqEWbaCLiHMi4lsRcW9EXFbz+B9HxO3Vz7cj4tFq+qkRcXNE3BUR34iIX+xqc3lEfLer3alzuU6LWdt6VY891fXYdV3TT4iIr0bEPRHx3yNi+Vytz2I3g+1rU9f02yPi7yPi/Ooxt69Z0qBePxER2yPia9XfvXO7HntP1e5bEfGapstUe23rFRFnRcRERNxR/T6jq814tcyp7et5c7lOi9kM6rU2Iv6uqyYf72qzvqrjvRHxkYiIuVynxWwG9bqo5/PrwNTnlNvX7GlQrzURcWNVq/GIWN312MXVPuA9EXFx1/TFt31l5qL7AY4AvgOsA5YDXwdOPsz8vwF8srr9IuDE6vbzgQeBY6v7lwM/P9/rt9h+ZlKv6v6+PvP9D+AN1e2PA7823+u6GH5mWq+u6c8GHgGeWd13+5qnetE5ifzXqtsnA/d13f46cBRwQrWcI6b7HvBnzur1CuD51e1TgO93tRkHxuZ7/RbbzwzrtRa4s89y/w/w00AAXwB+br7XdTH8zKRePfO8DNjVdd/ta/7q9Vng4ur2GcCfVrefDeyqfo9Wt0erxxbd9rVYR+heBdybmbsycz+wDXj9Yea/ELgSIDO/nZn3VLcfAB4CBv6Hds1I63r1U33bcgZwVTXpCuD8IfRVw6vXzwNfyMzHZ6GP+pEm9UrgmOr2jwEPVLdfD2zLzCcy87vAvdXypvseUHOt65WZX6s+twDuAp4REUfNQZ+XsplsX7Ui4jjgmMy8OTt7n5/Gz69hGVa9Bu6HaCia1Otk4Mbq9vaux18D3JCZj2TmXuAG4JzFun0t1kD3AuB7Xff3VNMOERFr6Hzz/OWax15F5xuB73RN/sNqWPeP/aAcmpnW6xkRcWtE7Jg6fA9YCTyamZODlqlpG8r2BbyBQz8Q3b6Gr0m9fhf4pYjYA3yezqjq4do2fg9o2mZSr24XAF/LzCe6pn2qOhzs/YviEKOFYab1OqE6tO8vI+LVXcvcM2CZamdY29cvcujnl9vX8DWp19fp/L0D+KfA0RGx8jBtF+X2tVgDXd2G1O//M7wBuCoznzpoAZ0E/6fAr2TmgWrye4AXA6+kM4T7W8Pp7pI303r9RGaOAf8C+JOI+MlpLlPTM6zt62XAF7smu33Njib1uhC4PDNXA+cCfxoRyw7T1u1r9sykXp0FRLwU+CDwr7raXJSZLwNeXf28cai9XrpmUq8H6Xx+vQJ4B/DfIuKYhstUO8PYvn4KeDwz7+xq4/Y1O5rU613AP46IrwH/GPg+MHmYtoty+1qsgW4PcHzX/dX0P8ThkFGC6g/qnwPvy8wdU9Mz88HseAL4FJ2hYM3cjOo1dYhRZu6icxz7K4C/BY6NiJEGy9T0zKhelV8APpeZT05NcPuaNU3q9RY655ySmTcDzwCec5i203kPaHpmUi+qCwJ8DvjlzHz66JLM/H71+zHgv+H2NSyt61UdyvxwNX2CztFAL6qWubqrvdvX8Mxo+6rU7Ye4fc2OgfXKzAcy859VX4y8t5r2/w7TdlFuX4s10N0CnBidqxwup7PxXdc7U0ScROdEyZu7pi2n82H46cz8bM/8x1W/g87xtt3fzqi9mdRrdOrQvIh4DnA68M3quOjtdM7TArgYuHZW12LpaF2vLoecf+D2NWua1Ot+4EyAiHgJnR2Yv6nme0NEHBURJwAn0jmZvNF7QK20rldEHEvny8j3ZOb/npo5Ikaqv49ExJHAebh9DctM6vXciDiimr6Ozva1KzMfBB6LiA3V38Nfxs+vYZnJ30Oqkbp/TudcLqppbl+zZ2C9IuI5XSOo7wE+Wd3+InB2tZ84CpwNfHHRbl/zfVWW2fqhM0z+bTrfeL23mvb7wOu65vld4AM97X4JeBK4vevn1OqxLwN30NlQPwOsmO/1XCw/M6jXz1Q1+Xr1+y1dj62js/N5L52rIB013+u5WH7a1quavpbOIRHLeqa7fc1TveicVP6/q+3oduDsrrbvrdp9i64rgdUt05/5rRfwPuCHPZ9fzwOeBUwA36BzsZQPA0fM93oulp8Z1OuCqh5fB24D/knXMseqv4XfAT4KxHyv52L5meHfw43Ajp7luX3Nb71+HrinmmcLXft6wJvp7APeS+cUqqnpi277imrFJEmSJEmFWayHXEqSJEnSomegkyRJkqRCGegkSZIkqVAGOkmSJEkLWkT8UUTcHRHfiIjPVVf27TfvERHxtYi4vmva5RHx3eofwN8eEadW038sIv4sIr4eEXdFxK9U09dExEQ1710R8avV9GdGxJ9XfbkrIj7Q9Ry/GhF3VG2+EhEnN1ivn4iIL0XEzoj4ZkSsnfZr40VRJEmSJC0UEbEReFNmvqlr2tnAlzNzMiI+CJCZv9Wn/TvoXM3ymMw8r5p2OXB9Zl7VM++/A34sM38rIp5L56rOPz71cGY+EREr6FwZ82eAR4Gfyszt1b9TuBH4D5n5hYg4JjN/UC33dcCvZ+Y5A9Z1HPjDzLyhep4Dmfl4s1eqwxE6SZIkSQtaZn4pMyeruzs4+B+EPy0iVgOvpfNvDBotGji6+r90K4BHgMnM3J+ZT1TzHEWVmzLz8czcXt3eT+ffjqyu7v+ga7nPqpY9NWL4RxFxSzXC+K+q6ScDI5l5Q9V+33TDHBjoJEmSJJXlzcAX+jz2J8C/BQ7UPPaHVaD644g4qpr2UeAlwAN0/h/u2zLzAEBEHB8R3wC+B3wwMx/oXlh12Oc/oTNKNzXtrRHxHeA/Ab9ZTX4L8P8y85XAK4F/GREnAC8CHo2Ia6pDRP8oIo6Y3kthoJMkSZK0AETEVyPidjqja6/rOt/tNV3zvBeYBLbWtD8PeCgzJ2oW/x7gxXQC1bOBqcM1X0Pnn8g/HzgV+GhEHAOQmd/LzH8IvBC4OCJWdT3XCHAl8JHM3DU1PTM/lpk/WS3/fdXks4Ffrtbtq8BK4ERgBHg18K6qX+uANzV8uZ5moJMkSZI07zLzpzLzVOAS4LrMPLX6+SJARFwMnAdclPUXAjmdThC8D9gGnBERn6mW/WB2PAF8CnhV1eZXgGuqx+4Fvksn+HX36wHgLjrha8pm4J7M/JM+q7MNOL+6HcBvdK3PCZn5JWAP8LXM3FUdTvo/gdMavVhdDHSSJEmSFrSIOIfOqNfr+p1nlpnvyczVmbkWeAOdi6j8UtX+uOp30Alad1bN7gfOrB5bBZwE7IqI1RHxD6rpo3TC4req+/8e+DHg7T19PLHr7muBe6rbXwR+LSKOrOZ7UUQ8C7gFGK0uxgJwBvDNab40jEy3gSRJkiTNsY/SuTjJDZ1Mxo7M/NWIeD6wJTPPHdB+axWcgs4hlr9aTf8D4PKIuKN67Lcy828j4izgv0REVtP/c2beUV105b3A3cBtVV8+mplbgH8dET8LPAnsBS6unmMLsLaaP4C/Ac7PzKci4l3AjdX0CeC/TveF8d8WSJIkSVKhPORSkiRJkgploJMkSZKkQhnoJEmSJKlQBjpJkiRJKpSBTpIkSZIKZaCTJEmSpEIZ6CRJkiSpUAY6SZIkSSrU/weGL5JUGyBlMAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1080x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15, 5))\n",
    "plt.grid()\n",
    "plt.plot(agg_regrouped['t_mid_tdb'], agg_regrouped['zero_point'], '.')\n",
    "#plt.ylim(20, 21)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3QAAAEyCAYAAACyDpLjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X10XHd95/H3V5YNBSUb4YAJebBxCyEhLQGpIEgLdtKwaeAUtpwthJSGguvtKS1hedhNe6CctvuQLoXTcMppME5wWNy4bRKakA0PaRrh5RBlsQLbPDiQ4I2DG4PBEU2ctMiyvvuHRkZWZjQPGs3M1X2/zvGR7p177/zufO9Yv8/93TsTmYkkSZIkqXj6ut0ASZIkSVJrDHSSJEmSVFAGOkmSJEkqKAOdJEmSJBWUgU6SJEmSCspAJ0mSJEkFZaCTJEmSpIIy0EmSJElSQRnoJEmSJKmg+rvdgGpOPPHEXLduXbeb0XFPPPEEz3zmM7vdDHWI9S4X610u1rs8rHW5WO9y6Xa9x8fHf5iZz663XN1AFxGnAp8BngtMA1sy84qIeBbw18A64CHg1zJzosr6pwFbgVOBBC7MzIcWes5169axa9euek1bdkZHR9mwYUO3m6EOsd7lYr3LxXqXh7UuF+tdLt2ud0TsbWS5Ri65nALel5lnACPAuyLiTOAy4LbMfAFwW2W6ms8AH6ms/3LgQCMNkyRJkiQtrG6gy8z9mXlX5ffHgd3AycAbgGsqi10DvHH+upXg15+Zt1bWP5SZT7ap7ZIkSZJUapGZjS8csQ7YCZwFPJyZJ8x5bCIzB+ct/0ZgEzAJPB/4e+CyzDxSZdubgc0Aa9asGdqxY0ez+1J4hw4dYmBgoNvNUIdY73Kx3uVivcvDWpeL9S6Xbtd748aN45k5XG+5hj8UJSIGgOuB92TmYxHRyGr9wC8CLwUeZuaeu7cDV81fMDO3AFsAhoeHs4zXJ3f7Ol11lvUuF+tdLta7PKx1uVjvcilKvRv62oKIWMlMmNuemTdUZn8/Ik6qPH4S1e+N2wd8IzP3ZOYU8HfAyxbfbEmSJElS3UAXM0NxVwG7M/Njcx66Cbik8vslwI1VVv86MBgRsx+3eS5wX+vNlSRJkiTNamSE7hzgbcC5EfHNyr8LgcuB8yPiAeD8yjQRMRwRWwEq98q9H7gtIu4GAvjUEuyHJEmSJJVO3XvoMvOrzASxas6rsvwuZj4IZXb6VuDnWm2gJEmSJKm6hu6hU/PG907widsfZHzvU75rXZIkSZLaouFPuVR1l9+ymy/e+z0uePFzuezCM4CZMPfmLXcwdSTpXxH89eZXMrR2sM6WJEmSJKk5jtA1qNqI2+W37ObKnXt46OCTXLlzD5ffshuAT37lO0wdmfl+v6kjySe/8p0FtzO7rf+888mj22h3WyVJkiQtP47QNWB87wS/duXXOJKwIuBvfvtVDK0dZPude49Zbvude7nswjO475F/Pmb+7PT43gku+tQYh6emWdnfx7W/NcLQ2sGjwRA4+nN2tG92vbE9BxlZv/qYkb5q82s9R739q7Z9SZIkSb3NEboGXP6F3VQG3DiSM9Ozv881O/30lSuOmT87fcNd+5icmiaByalpbrhrHwBfvPd7xyw/d3p87wRv/uTX+MiXvsWbP/m1o6Nu43snuHjrGB/98re4eOvY0fm1nmN2nfkjd7W2026OGh5rfO8EN39nshCvh7VTuzV7TLVyDHrcSpLKwhG6Btz3yGNVp9etfgb37X/86Px1q58BwDt+YT1/8Lm7j85/xy+sB+CB7z/OXLPTZ596Ag8dfPLo/LNPPeHo71d+5TtMTc/8PjU9M/2p3xhmbM9Bfnx4Jrj9+PA0Y3sOMrR2kAOP//iY55idHt87wUVb7uDwkWTliuDayn19c7czOWc79TQ7anjx1jEmp6ZZ1d/H9k31Rw2Xs9nX48eHp7n5obFFvx5LOcJq7VRPs8dfs8dUK8egx23jWqlft67oaOW529Ver2SR1MsMdA2YOjJddfrs0waPCXRnnzbzn/zpzz2OlSviaHg6/bnHAfDoE5PHbGd2+gVrjjtm/tzpA4/96zGPzU4//i+HmR0gzMo0wHOOe9oxy89OX3/XPiYrQ4iTR5Lr79rH0NpBBp+x6uh2poHBZ6w6uu5Coa1aZ6nW/LE9B5mcmmY64fDUsaHxr+58mC/cs59fPusk3vqK0+o+90KaCZndNPt6JE99PdrZOW7H67FQ7XrVcujAdaLjWuu918z2WwlOzR5TrZx0auW47USwWern6MVw3a72djPYt/rcvfR3R9LyZqBrwPE/tZIfHJo8ZhrgTS87het2ffdocHvTy04BZjoTR6ZnYtL0dB7tTDz/2QM8+IMnjm7n+c8eAGBk/WqevrKPycPTrFrZx8j61UeXefPPn8b/3Xf3MdMA9+4/dtRwdvpXX3YKfzu+7+g9dL9aadMP543czU5PPDlJX8B0Ql/MTMPCf8DmdpYm53SWao0ajqxfzar+vqNtmt2/v7rz4aMjmf/7gR8C8NZXnFb3PsBmOpbt/kPcjpA0+3pMHj729Whn57hdr0et2rXz9WjVUo4GtxKU6z22lPtQawS+llrvvWZfj1aCU71jar6FTjot9Bz9fTMn1Vb0Rd3n6MVRw060qZVw3a4TPM22t5Xnbld7m92OI8SSOs1A14Dnn/jMYwLd8098JgBDawe5dvMrn9Ihq9Vh+e3X/DS33/99pqahv29menY72zeNcO3ff52Lfunnj/mP/62vOI2HDz5x9KsRZjtdv3zWSUc7YrPTR9v0WyNPadOJ80buZqdH1q+mf8VMW/tX/KStC/0BG3zGKip5len8SQer1qjh0NpB/vD1Lz46GjC7nS/cs/+YNn3hnv289RWnHb0PEH5yH2C9zm6ts/gL7Ueznel2haRa9W5n57jWtpp9jtm2tmOktp3aGTCqHQfNBuXZ7TTzgUS1lm/2mIXaI/C11Hrv1VKrTa0E/nrH1Pz5tU461RUBZOXnwpodBexEuOhE2Go2XDe7/EKabW8rz92u9ja7nSJe2SCp2Ax0DfiZNcfxfx6aOGZ61tDawaf8R12rwzK0dpC//g+vqtqRGVo7yOM/varqGf9tdzzE5NQ02+54iPNf/FyG1g4e7XxVu2SqWptqjSYCkHnsT44NeitWHPsHbOLJSSpdJfr4SQer1qjh+N4J/vjme5mcmubrDz3K6c89jqG1gzVD6bzPmjlmutYfylpn8WudqW+lM92ukATV671Qp6FWR75WWK61rVafo9r+tPP1aHakr9WAUe15q41sNRuUgQVPRFRTa/laz71QmJwfWepFmFrvvVpqvY9aCfyz61X7v67a8q10ysf2HGTqyExAO3Kk/jHY7ChgsyOAs+ssZXhq5XVaqH7tWL6d7W3ludvV3ma3087gK0mNMNA14Kzn/ZsFp6up1QmuNb+WhTqQb33FaQ3d+zL7vNVGE8f2HGRqOmc6PnMuDwWYzpn503lsxBpZv5qnrXzqH6tancRa+1ArlC4UPmt1pBY8i1/lTH2znenZ525XSKpVo2qdhnqjQtXC8kInFZp9jlra9Xq0MtLXasCYP7/WyFat7Sy0bwudiKim1vK1nnuh/wtqXWpdy0InhGqqMeLVbOCvpdbyrXTKm+1QtzQK2MQIICx9eGo1vDT7N6nZ5RfaTisBrdnnbmd7m3lN2xV8JakRBroG3DPve+XmTy+ldp7pq/YHqdb2b7hr3zFfjj53tKHWH6tancSF9qFaKK0VPo+q0pFaaESl2pn6hTrT1Ua8FtrvWuu0EpKq1aiVUcNa26ql1VHGdoTGlkf6mggYtZ57oZGtattZqKO24Ch4FQst38x7dXb5apdaL6SZE0LNjnjVa2+zy7cSOpZyRKWV16MT2hVeOqVo7W3Gct43Sb3HQNeAZi9naqelPtNXa/v1Rhtq/bGqFdDadSa2Vkeq2RGVWp3pWiNeC6m1TrsuPWxl1HCh7Td7WdtCo4y16tRMMG1lpK/ZDnWt5252ZKvePi94IqINyy/0PlrKDmTRLueb3d5Sjai08np0+xMlJUnLl4GuAa10+tppqc/0Vdt+s6MNrTxHK5o9i7/QKFKtS1AX+kCKZkaY6oWkm78zyXHPn6i7/YU6m812RJu9rK3VTmUzwbTZkb56x0E1Cz13syNbC2llJGkpl2+XVsNWUfav2edu5fVo1yWokiTNZ6BrQLs7fUXQ7OhBp7RztK/Zy9qaDW71QtL8LxZv9fLJZjqizQbiVkcZWwmmjdZodn67Rn+6GSKKxNfpWM2+Hu28BFWSpLkMdA0qY2emV/d5KdvVyqhQs2FhNiTN/2LxTnTgOnFpWbvu62tkX4oy+iN1+xJUSdLyZaCT5mllVKiVUbL5XyzeqQ7cUl9a5siCVJ0nISRJS8FAJzWhHR2s2ZBU7Yvke7ED10on1JEFSZKkzjDQSV0wtLb6F8kvF70YTCVJkpajvm43QJIkSZLUGgOdJEmSJBWUgU6SJEmSCspAJ0mSJEkFZaCTJEmSpIIy0EmSJElSQRnoJEmSJKmgDHSSJEmSVFB1A11EnBoRt0fE7oi4NyIurcx/VkTcGhEPVH7W/BbhiDg+Iv4pIv6inY2XJEmSpDJrZIRuCnhfZp4BjADviogzgcuA2zLzBcBtlela/gT4ymIbK0mSJEn6ibqBLjP3Z+Zdld8fB3YDJwNvAK6pLHYN8MZq60fEELAG+HI7GixJkiRJmhGZ2fjCEeuAncBZwMOZecKcxyYyc3De8n3APwBvA84DhjPzd2tsezOwGWDNmjVDO3bsaGpHloNDhw4xMDDQ7WaoQ6x3uVjvcrHe5WGty8V6l0u3671x48bxzByut1x/oxuMiAHgeuA9mflYRDSy2u8At2Tmd+stn5lbgC0Aw8PDuWHDhkabtmyMjo5Sxv0uK+tdLta7XKx3eVjrcrHe5VKUejcU6CJiJTNhbntm3lCZ/f2IOCkz90fEScCBKqu+EvjFiPgdYABYFRGHMnOh++0kSZIkSQ2oG+hiZmjtKmB3Zn5szkM3AZcAl1d+3jh/3cy8eM523s7MJZeGOUmSJElqg0Y+5fIcZu6BOzcivln5dyEzQe78iHgAOL8yTUQMR8TWJWuxJEmSJAloYIQuM78K1LoB7rwqy+8CNlWZvw3Y1lzzJEmSJEm1NDJCJ0mSJEnqQQY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFVTdQBcRp0bE7RGxOyLujYhLK/OfFRG3RsQDlZ+DVdY9OyLuqKz3jxHx5qXYCUmSJEkqo0ZG6KaA92XmGcAI8K6IOBO4DLgtM18A3FaZnu9J4Dcy88XABcCfR8QJ7Wm6JEmSJJVb3UCXmfsz867K748Du4GTgTcA11QWuwZ4Y5V1v52ZD1R+fwQ4ADy7PU2XJEmSpHKLzGx84Yh1wE7gLODhzDxhzmMTmfmUyy7nPP5yZoLfizNzusrjm4HNAGvWrBnasWNHw+1aLg4dOsTAwEC3m6EOsd7lYr3LxXqXh7UuF+tdLt2u98aNG8czc7jecv2NbjAiBoDrgfdk5mMR0XBjIuIk4H8Cl1QLcwCZuQXYAjA8PJwbNmxoePvLxejoKGXc77Ky3uVivcvFepeHtS4X610uRal3Q59yGRErmQlz2zPzhsrs71eC2mxgO1Bj3eOB/wV8MDPHFt9kSZIkSRI09imXAVwF7M7Mj8156CbgksrvlwA3Vll3FfA54DOZ+beLb64kSZIkaVYjI3TnAG8Dzo2Ib1b+XQhcDpwfEQ8A51emiYjhiNhaWffXgFcDb5+z7tnt3w1JkiRJKp+699Bl5leBWjfMnVdl+V3ApsrvnwU+u5gGSpIkSZKqa+geOkmSJElS7zHQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZK0zI3vneATtz/I+N6JbjdFktRm/d1ugCRJWjrjeye4eOsYk1PTrOrvY/umEYbWDna7WZKkNnGETpKkZWxsz0Emp6aZTjg8Nc3YnoPdbpIkqY0MdJIkLWMj61ezqr+PFQEr+/sYWb+6202SJLWRl1xKkrSMDa0dZPumEcb2HGRk/Wovt5SkZcZAJ0nSMje0dtAgJ0nLlJdcSpIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUHUDXUScGhG3R8TuiLg3Ii6tzH9WRNwaEQ9Ufg7WWP+SyjIPRMQl7d4BSZIkSSqrRkbopoD3ZeYZwAjwrog4E7gMuC0zXwDcVpk+RkQ8C/gw8Arg5cCHawU/SZIkSVJz6ga6zNyfmXdVfn8c2A2cDLwBuKay2DXAG6us/m+BWzPz0cycAG4FLmhHwyVJkiSp7Jq6hy4i1gEvBe4E1mTmfpgJfcBzqqxyMvDdOdP7KvMkSZIkSYsUmdnYghEDwFeA/5qZN0TEjzLzhDmPT2Tm4Lx1PgA8LTP/S2X6Q8CTmfnRKtvfDGwGWLNmzdCOHTta3afCOnToEAMDA91uhjrEepeL9S4X610e1rpcrHe5dLveGzduHM/M4XrL9TeysYhYCVwPbM/MGyqzvx8RJ2Xm/og4CThQZdV9wIY506cAo9WeIzO3AFsAhoeHc8OGDdUWW9ZGR0cp436XlfUuF+tdLta7PKx1uVjvcilKvRv5lMsArgJ2Z+bH5jx0EzD7qZWXADdWWf1LwGsjYrDyYSivrcyTJEmSJC1SI/fQnQO8DTg3Ir5Z+XchcDlwfkQ8AJxfmSYihiNiK0BmPgr8CfD1yr8/rsyTJEmSJC1S3UsuM/OrQNR4+Lwqy+8CNs2Zvhq4utUGSpIkSZKqa+pTLiVJkiRJvcNAJ0mSJEkFZaCTJEmSpIIy0EmSJElSQRnoJEmSJKmgDHSSJEmSVFAGOkmSJEkqKAOdJEmSJBWUgU6SJEmSCspAJ0mSJEkFZaCTJEmSpIIy0EmSJElSQRnoJEmSJKmgDHSSJEmSVFAGOkmSJEkqKAOdJEmSJBWUgU6SJEmSCspAJ0mS1Ebjeyf4xO0PMr53ottNkXpGr74verVdzejvdgMkSZKWi/G9E1y8dYzJqWlW9fexfdMIQ2sHF7W9sT0HGVm/elHbkaB7x1O73xfLvV3NMtBJkiS1ydieg0xOTTOdcHhqmrE9B1vuIC6XzqZ6QzePp3a+L6B9wbTd7eoWL7mUJElqk5H1q1nV38eKgJX9fYysX93ytqp1NqVWdfN4auf7YjaYfvTL3+LirWOLulSyne3qJkfoJElSz2vnpWJLednZ0NpBtm8aacv2Zzubh6emC93ZLJp2HR+9drlsq8dTO/ajne+Ldo6qtbNd3WSgkyRJi7LUHdd2XirWicvOhtYOtmWb3e5s9logaUWz+9Cu46MXL5dt5Xhq5360633R7hMd7WpXNxnoJElSyzrRcW3nGfmi3TPTbGeznaNLvRZImtXKPrTr+OjV46zZ46kX96PbJzp6kYFOkiS1rBMdvnaekV/OlzG2M4T1Yke+Wa3sQ7uOj+VynPXqfiyHUbV2MtBJkqSWdaLD184z8sv57H47Q1ivduSb0co+tOv4WC7H2XLZj+XOQCdJklrWqQ5fO8/IL9ez++0MYcuhI9/qPrTzHsgivm7zLZf9WM7qBrqIuBp4PXAgM8+qzHsJcCUwADwEXJyZj1VZ9z8Cm4AE7gZ+MzP/tW2tlyRJbdfsfVh2+HpDu0NYt+7fa6flcmz24mur3tHICN024C+Az8yZtxV4f2Z+JSLeAXwA+NDclSLiZODdwJmZ+S8R8TfAWyrbkyRJPWg5fBhGmXUrwHjcLB1fW9VT94vFM3Mn8Oi82acDOyu/3wq8qcbq/cBPRUQ/8AzgkRbbKUmSOsAvs1YrPG6Wjq+t6mn1Hrp7gF8BbgT+PXDq/AUy858i4s+Ah4F/Ab6cmV+utcGI2AxsBlizZg2jo6MtNq24Dh06VMr9LivrXS7Wu1yKXO+n/egI/QFTCSsCnvajvYyO7ut2s3pWkWvdTmU5brpR77K8tr2oKO/vyMz6C0WsA26ecw/di4CPA6uBm4B3Z+bqeesMAtcDbwZ+BPwtcF1mfrbe8w0PD+euXbua2pHlYHR0lA0bNnS7GeoQ610u1rtcil5v79dpXNFr3U5lOG66Ve8yvLa9qNvv74gYz8zhesu1NEKXmfcDr6080QuB11VZ7JeA/5eZP6gsdwPwKqBuoJMkSd2zXD5IQp3lcbN0fG21kLr30FUTEc+p/OwDPsjMJ17O9zAwEhHPiIgAzgN2t9pQSZIkSdKx6ga6iLgWuAM4PSL2RcQ7gYsi4tvA/cx80MmnK8s+LyJuAcjMO4HrgLuY+cqCPmDLkuyFJEmSJJVQ3UsuM/OiGg9dUWXZR4AL50x/GPhwy62TJEmSJNXU0iWXkiRJkqTuM9BJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoOoGuoi4OiIORMQ9c+a9JCLuiIi7I+LzEXF8jXVPiIjrIuL+iNgdEa9sZ+MlSZIkqcwaGaHbBlwwb95W4LLM/Fngc8AHaqx7BfDFzHwR8BJgd4vtlCRJkiTNUzfQZeZO4NF5s08HdlZ+vxV40/z1KqN2rwauqmxnMjN/tKjWSpIkSZKOisysv1DEOuDmzDyrMv014E8z88aIeC/wR5l53Lx1zga2APcxMzo3DlyamU/UeI7NwGaANWvWDO3YsaPVfSqsQ4cOMTAw0O1mqEOsd7lY73Kx3uVhrcvFepdLt+u9cePG8cwcrrdcq4HuRcDHgdXATcC7M3P1vHWGgTHgnMy8MyKuAB7LzA/Ve77h4eHctWtX3XYtN6Ojo2zYsKHbzVCHWO9ysd7lYr3Lw1qXi/Uul27XOyIaCnQtfcplZt6fma/NzCHgWuA7VRbbB+zLzDsr09cBL2vl+SRJkiRJT9VSoIuI51R+9gEfBK6cv0xmfg/4bkScXpl1HjOXX0qSJEmS2qCRry24FrgDOD0i9kXEO4GLIuLbwP3AI8CnK8s+LyJumbP67wHbI+IfgbOB/9buHZAkSZKksuqvt0BmXlTjoSuqLPsIcOGc6W8Cda/7lCRJkiQ1r6VLLiVJkiRJ3WegkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkEZ6CRJkiSpoAx0kiRJklRQBjpJkiRJKigDnSRJkiQVlIFOkiRJkgrKQCdJkiRJBWWgkyRJkqSCMtBJkiRJUkHVDXQRcXVEHIiIe+bMe0lE3BERd0fE5yPi+AXWXxER34iIm9vVaEmSJElSYyN024AL5s3bClyWmT8LfA74wALrXwrsbql1kiRJkqSa6ga6zNwJPDpv9unAzsrvtwJvqrZuRJwCvI6ZAChJkiRJaqPIzPoLRawDbs7MsyrTXwP+NDNvjIj3An+UmcdVWe864L8DxwHvz8zXL/Acm4HNAGvWrBnasWNH83tTcIcOHWJgYKDbzVCHWO9ysd7lYr3Lw1qXi/Uul27Xe+PGjeOZOVxvuf4Wt/8O4OMR8YfATcDk/AUi4vXAgcwcj4gN9TaYmVuALQDDw8O5YUPdVZad0dFRyrjfZWW9y8V6l4v1Lg9rXS7Wu1yKUu+WAl1m3g+8FiAiXsjMZZXznQP8SkRcCDwdOD4iPpuZv95qYyVJkiRJP9HS1xZExHMqP/uADwJXzl8mM38/M0/JzHXAW4B/MMxJkiRJUvs08rUF1wJ3AKdHxL6IeCdwUUR8G7gfeAT4dGXZ50XELUvZYEmSJEnSjLqXXGbmRTUeuqLKso8AF1aZPwqMNtk2SZIkSdICWrrkUpIkSZLUfQY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZLUceN7J/jE7Q8yvnei202RpELr73YDJElSuYzvneDirWNMTk2zqr+P7ZtGGFo72O1mSVIhOUInSZI6amzPQSanpplOODw1zdieg91ukiQVloFOkiR11Mj61azq72NFwMr+PkbWr+52kySpsLzkUpIkddTQ2kG2bxphbM9BRtav9nJLSVoEA50kSeq4obWDBjlJagMvuZQkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIKykAnSZIkSQVloJMkSZKkgjLQSZIkqaeM753gE7c/yPjeiW43Rep5/d1ugCRJkjRrfO8EF28dY3JqmlX9fWzfNMLQ2sFuN0vqWY7QSZIkqWeM7TnI5NQ00wmHp6YZ23Ow202SelrdQBcRV0fEgYi4Z868l0TEHRFxd0R8PiKOr7LeqRFxe0Tsjoh7I+LSdjdekiRJy8vI+tWs6u9jRcDK/j5G1q/udpOkntbICN0LooYGAAAHQElEQVQ24IJ587YCl2XmzwKfAz5QZb0p4H2ZeQYwArwrIs5cRFslSZK0zA2tHWT7phHe+9rTvdxSakDde+gyc2dErJs3+3RgZ+X3W4EvAR+at95+YH/l98cjYjdwMnDf4posSZKk5Wxo7aBBTmpQZGb9hWYC3c2ZeVZl+mvAn2bmjRHxXuCPMvO4OuvvBM7KzMdqLLMZ2AywZs2aoR07djS1I8vBoUOHGBgY6HYz1CHWu1ysd7lY7/Kw1uVivcul2/XeuHHjeGYO11uu1U+5fAfw8Yj4Q+AmYLLWghExAFwPvKdWmAPIzC3AFoDh4eHcsGFDi00rrtHRUcq432VlvcvFepeL9S4Pa10u1rtcilLvlgJdZt4PvBYgIl4IvK7achGxkpkwtz0zb2i1kZIkSd0yvneCsT0HGVm/2ssAJfWclgJdRDwnMw9ERB/wQeDKKssEcBWwOzM/trhmSpIkdZ7fiSap1zXytQXXAncAp0fEvoh4J3BRRHwbuB94BPh0ZdnnRcQtlVXPAd4GnBsR36z8u3BJ9kKSJGkJ+J1oknpdI59yeVGNh66osuwjwIWV378KxKJaJ0mS1EWz34l2eGra70ST1JNa/VAUSZKkZW/2O9G8h05SrzLQSZIkLcDvRJPUy+reQydJkiRJ6k0GOkmSJEkqKAOdJEmSJBWUgU6SJEmSCspAJ0mSJEkFZaCTJEmSpIIy0EmSJElSQRnoJEmSJKmgDHSSJEmSVFCRmd1uw1NExA+Avd1uRxecCPyw241Qx1jvcrHe5WK9y8Nal4v1Lpdu13ttZj673kI9GejKKiJ2ZeZwt9uhzrDe5WK9y8V6l4e1LhfrXS5FqbeXXEqSJElSQRnoJEmSJKmgDHS9ZUu3G6COst7lYr3LxXqXh7UuF+tdLoWot/fQSZIkSVJBOUInSZIkSQVloJMkSZKkgjLQdUBEXB0RByLinhqPR0R8PCIejIh/jIiXzXv8+Ij4p4j4i860WIuxmHpHxGkR8eWI2B0R90XEuk61W61ZZL3/R0TcW6n3xyMiOtdytaKBer8oIu6IiB9HxPvnPXZBRHyrcixc1pkWq1Wt1joiTo2I2yvv63sj4tLOtVqtWsx7u/L4ioj4RkTcvPSt1WIt8v/yEyLiuoi4v/I+f2VnWl2bga4ztgEXLPD4LwMvqPzbDPzlvMf/BPjKkrRMS2Ebrdf7M8BHMvMM4OXAgSVqo9pnGy3UOyJeBZwD/BxwFvDzwGuWsqFqi20sXO9HgXcDfzZ3ZkSsAD7BzPFwJnBRRJy5RG1Ue2yjhVoDU8D7Kv+PjwDvstaFsI3W6j3rUmB3m9ukpbON1ut9BfDFzHwR8BJ6oO4Gug7IzJ3MHBi1vAH4TM4YA06IiJMAImIIWAN8eelbqnZotd6VP/j9mXlrZTuHMvPJDjRZi7CI93cCTwdWAU8DVgLfX+r2anHq1TszD2Tm14HD8x56OfBgZu7JzElgBzPHhnpUq7XOzP2ZeVfl98eZ6eydvJRt1eIt4r1NRJwCvA7YunQtVDu1Wu+IOB54NXBVZbnJzPzRUra1EQa63nAy8N050/uAkyOiD/go8IGutEpLpWq9gRcCP4qIGyqXbXykclZfxVa13pl5B3A7sL/y70uZ2fWzfFoytd73WsYql82/FLizuy3REvtz4D8B091uiJbceuAHwKcrfbWtEfHMbjfKQNcbqt03k8DvALdk5nerPK7iqlXvfuAXgfczc/ndeuDtnWuWlkjVekfEzwBnAKcw07E/NyJe3dGWqZNqve+1TEXEAHA98J7MfKzb7dHSiIjXAwcyc7zbbVFH9AMvA/4yM18KPAF0/Z5oA11v2AecOmf6FOAR4JXA70bEQ8xcw/sbEXF555unNqtV733ANyqXZE0Bf8fMfxoqtlr1/nfAWOXS2kPAF5i530bLU63jQMtQRKxkJsxtz8wbut0eLalzgF+p9NV2MHNy7rPdbZKW0D5gX2bOjrpfRw/01Qx0veEmZsJaRMQI8M+Va/AvzszTMnMdM6M2n8nMrp8F0KJVrTfwdWAwIp5dWe5c4L5uNVJtU6veDwOviYj+SufvNfTAjdVaMl8HXhARz4+IVcBbmDk2tMxUPq32KmB3Zn6s2+3R0srM38/MUyp9tbcA/5CZv97lZmmJZOb3gO9GxOmVWefRA321/m43oAwi4lpgA3BiROwDPszMByCQmVcCtwAXAg8CTwK/2Z2Wqh1arXdmHql8NO5tlQ7BOPCpju+AmrKI9/d1zIT2u5m59O6Lmfn5jjZeTatX74h4LrALOB6Yjoj3AGdm5mMR8bvAl4AVwNWZeW839kGNabXWzHxy7duAuyPim5XN/UFm3tLhXVATFvPe7lKTtQiLrPfvAdsrJ+f20AP99sj0En5JkiRJKiIvuZQkSZKkgjLQSZIkSVJBGegkSZIkqaAMdJIkSZJUUAY6SZIkSSooA50kSZIkFZSBTpIkSZIK6v8Dxk1X9ZCxfEYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1080x360 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(15, 5))\n",
    "\n",
    "plt.plot(agg_regrouped['airmass'], agg_regrouped['zero_point'], '.')\n",
    "#plt.ylim(20, 21)\n",
    "plt.grid()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "gain = 1.47\n",
    "read_noise = 9.42\n",
    "dark_current = 0.01"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "mags = Table.read('calibrated_magnitudes_kelt_16_faint_13.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "snr = (gain * mags['aperture_net_flux'] / \n",
    "       np.sqrt(gain * mags['aperture_net_flux'] + \n",
    "               mags['aperture_area'] * (1 + mags['aperture_area']/mags['annulus_area']) *\n",
    "               (gain * mags['sky_per_pix_avg'] + \n",
    "                gain * dark_current * mags['exposure'] +\n",
    "                read_noise**2\n",
    "               )\n",
    "              ))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1612.506328384899"
      ]
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "snr.max()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;Column name=&apos;aperture_net_flux&apos; dtype=&apos;float64&apos; unit=&apos;adu&apos; length=228416&gt;\n",
       "<table>\n",
       "<tr><td>0.03554203721814884</td></tr>\n",
       "<tr><td>0.010865317328612463</td></tr>\n",
       "<tr><td>0.21162200178674756</td></tr>\n",
       "<tr><td>0.028532498737706097</td></tr>\n",
       "<tr><td>0.013127131925856725</td></tr>\n",
       "<tr><td>0.09620434270592455</td></tr>\n",
       "<tr><td>0.08801476817405915</td></tr>\n",
       "<tr><td>0.06594961678474937</td></tr>\n",
       "<tr><td>0.09575037644968758</td></tr>\n",
       "<tr><td>0.005352178717360018</td></tr>\n",
       "<tr><td>0.03523151162172711</td></tr>\n",
       "<tr><td>0.08240132602802241</td></tr>\n",
       "<tr><td>...</td></tr>\n",
       "<tr><td>0.2080079958139937</td></tr>\n",
       "<tr><td>0.41598596723909026</td></tr>\n",
       "<tr><td>0.26216596171042117</td></tr>\n",
       "<tr><td>0.3810401372674711</td></tr>\n",
       "<tr><td>0.126141972104651</td></tr>\n",
       "<tr><td>0.07142793901418203</td></tr>\n",
       "<tr><td>1.0873130222307894</td></tr>\n",
       "<tr><td>0.26812608072554045</td></tr>\n",
       "<tr><td>0.17867986111223463</td></tr>\n",
       "<tr><td>0.1442924276151822</td></tr>\n",
       "<tr><td>0.004799602972304384</td></tr>\n",
       "<tr><td>0.33344694042018297</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Column name='aperture_net_flux' dtype='float64' unit='adu' length=228416>\n",
       " 0.03554203721814884\n",
       "0.010865317328612463\n",
       " 0.21162200178674756\n",
       "0.028532498737706097\n",
       "0.013127131925856725\n",
       " 0.09620434270592455\n",
       " 0.08801476817405915\n",
       " 0.06594961678474937\n",
       " 0.09575037644968758\n",
       "0.005352178717360018\n",
       " 0.03523151162172711\n",
       " 0.08240132602802241\n",
       "                 ...\n",
       "  0.2080079958139937\n",
       " 0.41598596723909026\n",
       " 0.26216596171042117\n",
       "  0.3810401372674711\n",
       "   0.126141972104651\n",
       " 0.07142793901418203\n",
       "  1.0873130222307894\n",
       " 0.26812608072554045\n",
       " 0.17867986111223463\n",
       "  0.1442924276151822\n",
       "0.004799602972304384\n",
       " 0.33344694042018297"
      ]
     },
     "execution_count": 27,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mag_error = 1 / snr\n",
    "mag_error"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [],
   "source": [
    "mags['mag_error'] = mag_error"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: VerifyWarning: Keyword name 'aperture_photometry_args' is greater than 8 characters or contains characters not allowed by the FITS standard; a HIERARCH card will be created. [astropy.io.fits.card]\n"
     ]
    }
   ],
   "source": [
    "mags.write('calibrated_magnitudes_kelt_16_faint_13.fits', overwrite=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<i>Table length=228416</i>\n",
       "<table id=\"table1992664016992\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>id</th><th>xcenter</th><th>ycenter</th><th>aperture_sum</th><th>annulus_sum</th><th>RA</th><th>Dec</th><th>sky_per_pix_avg</th><th>sky_per_pix_med</th><th>sky_per_pix_std</th><th>aperture</th><th>aperture_area</th><th>annulus_inner</th><th>annulus_outer</th><th>annulus_area</th><th>exposure</th><th>date-obs</th><th>night</th><th>aperture_net_flux</th><th>mag_inst_r</th><th>airmass</th><th>filter</th><th>file</th><th>star_id</th><th>mag_cal</th><th>zero_point</th><th>color_term</th><th>t_start_tdb</th><th>t_mid_tdb</th><th>vsx_name</th><th>mag_error</th></tr></thead>\n",
       "<thead><tr><th></th><th>pix</th><th>pix</th><th>adu</th><th></th><th>deg</th><th>deg</th><th>adu</th><th>adu</th><th>adu</th><th>pix</th><th></th><th>pix</th><th>pix</th><th></th><th>s</th><th></th><th></th><th>adu</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>adu</th></tr></thead>\n",
       "<thead><tr><th>int32</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>bytes19</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>bytes1</th><th>bytes28</th><th>int32</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>bytes50</th><th>float64</th></tr></thead>\n",
       "<tr><td>1422</td><td>4011.091526838012</td><td>3969.142916177813</td><td>19740.155430793762</td><td>103860.89329612255</td><td>314.5432670134417</td><td>32.01973136240193</td><td>31.808280954476942</td><td>31.461752891540527</td><td>10.803624996914255</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>5350.428135927457</td><td>-5.293886548969706</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1479</td><td>15.214845041740617</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.726455484</td><td>2458323.7268027063</td><td></td><td>0.03554203721814884</td></tr>\n",
       "<tr><td>1431</td><td>3863.942898239571</td><td>3998.9143066290862</td><td>33952.74383163452</td><td>102687.83718004823</td><td>314.51558576760186</td><td>32.01961637265539</td><td>31.504446114354263</td><td>31.2497501373291</td><td>10.520896497204099</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>19700.46818020352</td><td>-6.709106579064407</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1489</td><td>13.799244435989024</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.7264559236</td><td>2458323.726803146</td><td></td><td>0.010865317328612463</td></tr>\n",
       "<tr><td>1411</td><td>3963.0904975394483</td><td>3951.90486856516</td><td>15249.548798561096</td><td>104428.82190561295</td><td>314.5352479090502</td><td>32.015567615887285</td><td>31.808294335373855</td><td>31.52717399597168</td><td>10.731908677068274</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>859.8154503196401</td><td>-3.3089283230064974</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1468</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.7264559385</td><td>2458323.7268031607</td><td></td><td>0.21162200178674756</td></tr>\n",
       "<tr><td>1421</td><td>3923.008965618607</td><td>3969.8875728993826</td><td>20822.206637382507</td><td>101420.80074810982</td><td>314.52733781030184</td><td>32.017048316638714</td><td>31.151097199642955</td><td>30.903825759887695</td><td>10.454439681919503</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>6729.782269015484</td><td>-5.542917744838749</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1478</td><td>14.962746853484397</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.726455944</td><td>2458323.7268031663</td><td></td><td>0.028532498737706097</td></tr>\n",
       "<tr><td>1405</td><td>4015.1792324262733</td><td>3923.173770284605</td><td>30234.403473317623</td><td>104178.87234064937</td><td>314.54572817533784</td><td>32.01282162890551</td><td>31.792475467154404</td><td>31.450231552124023</td><td>10.459912367244256</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>15851.826412462999</td><td>-6.473113480632851</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1461</td><td>14.035685270542295</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.726455993</td><td>2458323.726803215</td><td></td><td>0.013127131925856725</td></tr>\n",
       "<tr><td>1417</td><td>3848.017315598605</td><td>3961.9101672544393</td><td>16185.751720428467</td><td>115460.04435215145</td><td>314.5140998526454</td><td>32.013443881563674</td><td>31.558219960992364</td><td>30.922298431396484</td><td>10.712969735233127</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>1909.1493538937557</td><td>-4.175014973009543</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1474</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.7264564354</td><td>2458323.7268036576</td><td></td><td>0.09620434270592455</td></tr>\n",
       "<tr><td>1419</td><td>3723.004118947109</td><td>3968.894287448843</td><td>16342.364827394485</td><td>102438.52531147748</td><td>314.49127182053485</td><td>32.01053792907989</td><td>31.50461414878146</td><td>31.270769119262695</td><td>10.508239826067275</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>2090.013158979509</td><td>-4.273287762138963</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1476</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.7264570193</td><td>2458323.7268042415</td><td></td><td>0.08801476817405915</td></tr>\n",
       "<tr><td>1438</td><td>3604.9472263665893</td><td>4021.055869255223</td><td>16588.801261663437</td><td>99055.38167959452</td><td>314.4680014835236</td><td>32.014766868890604</td><td>30.492322640378536</td><td>30.28482151031494</td><td>9.544357245216066</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>2794.399482765415</td><td>-4.5886364415189576</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1497</td><td>15.915617787974224</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.7264570445</td><td>2458323.7268042667</td><td></td><td>0.06594961678474937</td></tr>\n",
       "<tr><td>1382</td><td>3932.032866076106</td><td>3855.864058718527</td><td>16039.92415869236</td><td>101568.43381509185</td><td>314.5332432307907</td><td>31.999874377180166</td><td>31.223004554285843</td><td>31.05585479736328</td><td>10.306347741275532</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T05:19:53</td><td>58322</td><td>1914.969669465072</td><td>-4.178319960213055</td><td>1.121</td><td>r</td><td>kelt-16-S001-R001-C001-r.fit</td><td>1436</td><td>nan</td><td>20.524021782107393</td><td>-0.022386809287167866</td><td>2458323.726457212</td><td>2458323.7268044343</td><td></td><td>0.09575037644968758</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></tr>\n",
       "<tr><td>107</td><td>182.0508911360561</td><td>303.05677630682095</td><td>22842.986793518066</td><td>158728.4168729782</td><td>314.0267663163187</td><td>31.347071117831405</td><td>48.82892485615355</td><td>48.830955505371094</td><td>9.676146609681743</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>753.3016015657413</td><td>-3.1653374378374917</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>134</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027719204</td><td>2458323.9031191426</td><td></td><td>0.26216596171042117</td></tr>\n",
       "<tr><td>39</td><td>541.0001708205461</td><td>135.01320125036852</td><td>22581.07301712036</td><td>158686.21473121643</td><td>314.09743360865417</td><td>31.33297236070604</td><td>48.77210934622453</td><td>48.38126182556152</td><td>10.314663448435162</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>517.0905563268607</td><td>-2.7568317266208298</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>64</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027719847</td><td>2458323.903119207</td><td></td><td>0.3810401372674711</td></tr>\n",
       "<tr><td>97</td><td>193.00622531012922</td><td>280.0069870918759</td><td>23552.252098083496</td><td>158362.16340827942</td><td>314.02960264845143</td><td>31.34390001098602</td><td>48.58026260878469</td><td>48.42496299743652</td><td>9.545562662829555</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1575.0590566276842</td><td>-3.966157316521702</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>124</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027721277</td><td>2458323.90311935</td><td>ASASSN-V J205607.11+312038.1</td><td>0.126141972104651</td></tr>\n",
       "<tr><td>41</td><td>503.9609772205891</td><td>138.00226718367384</td><td>24672.248279571533</td><td>158899.16868400574</td><td>314.0906849908972</td><td>31.33223335878</td><td>48.33081605587856</td><td>48.28989028930664</td><td>9.58830815061887</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>2807.9022000782643</td><td>-4.593870153756231</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>66</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027721426</td><td>2458323.9031193648</td><td></td><td>0.07142793901418203</td></tr>\n",
       "<tr><td>14</td><td>602.9782312415216</td><td>75.94704453974062</td><td>21924.57398033142</td><td>155915.59290409088</td><td>314.1107678715254</td><td>31.325936902368422</td><td>48.065853973192866</td><td>48.05634117126465</td><td>9.302963672902422</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>180.09392311026386</td><td>-1.6116628576991223</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>38</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027723395</td><td>2458323.9031195617</td><td></td><td>1.0873130222307894</td></tr>\n",
       "<tr><td>80</td><td>264.0154239709242</td><td>209.05195133865692</td><td>22760.32354927063</td><td>157921.01183128357</td><td>314.04501177829536</td><td>31.33534466061794</td><td>48.68451682124661</td><td>48.531917572021484</td><td>9.813231068559666</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>735.9670132262545</td><td>-3.1400610840687304</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>105</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027725724</td><td>2458323.9031197946</td><td></td><td>0.26812608072554045</td></tr>\n",
       "<tr><td>98</td><td>83.04174328655962</td><td>284.8475961628155</td><td>23363.15751838684</td><td>160498.44156837463</td><td>314.0097193604797</td><td>31.34107838620185</td><td>49.18918774371364</td><td>49.234052658081055</td><td>9.755796233160744</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1110.493235742062</td><td>-3.5867050046852746</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>125</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027726413</td><td>2458323.9031198635</td><td></td><td>0.17867986111223463</td></tr>\n",
       "<tr><td>18</td><td>306.91673956570503</td><td>80.1025961804836</td><td>23178.463884353638</td><td>156871.81446647644</td><td>314.0575777132189</td><td>31.317007021227244</td><td>48.20240802911838</td><td>48.1168098449707</td><td>9.73735439939401</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>1372.2082276089386</td><td>-3.816465258361166</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>42</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.902773827</td><td>2458323.903121049</td><td></td><td>0.1442924276151822</td></tr>\n",
       "<tr><td>37</td><td>102.38009225045143</td><td>130.41017764070943</td><td>80841.77911376953</td><td>162045.0770702362</td><td>314.0190387942966</td><td>31.318082075902257</td><td>49.66731724692211</td><td>49.48946762084961</td><td>9.88967002386121</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>58372.81413972158</td><td>-7.8884417885370715</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>61</td><td>9.262596422862337</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.9027743097</td><td>2458323.903121532</td><td></td><td>0.004799602972304384</td></tr>\n",
       "<tr><td>10</td><td>113.92856473017127</td><td>60.970711391059524</td><td>23356.94402694702</td><td>163692.00853157043</td><td>314.02373840609386</td><td>31.30783464202709</td><td>50.31406455831275</td><td>50.1636962890625</td><td>9.865457764126704</td><td>12.0</td><td>452.3893421169302</td><td>27.0</td><td>42.0</td><td>3251.548396465436</td><td>60.0</td><td>2018-07-24T09:33:41</td><td>58322</td><td>595.3974621831621</td><td>-2.909932658865633</td><td>1.159</td><td>r</td><td>kelt-16-S001-R001-C166-r.fit</td><td>34</td><td>nan</td><td>19.861563724094722</td><td>-1.4860336530214842</td><td>2458323.902775045</td><td>2458323.903122267</td><td></td><td>0.33344694042018297</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=228416>\n",
       "  id       xcenter       ...           vsx_name                mag_error      \n",
       "             pix         ...                                      adu         \n",
       "int32      float64       ...           bytes50                  float64       \n",
       "----- ------------------ ... ---------------------------- --------------------\n",
       " 1422  4011.091526838012 ...                               0.03554203721814884\n",
       " 1431  3863.942898239571 ...                              0.010865317328612463\n",
       " 1411 3963.0904975394483 ...                               0.21162200178674756\n",
       " 1421  3923.008965618607 ...                              0.028532498737706097\n",
       " 1405 4015.1792324262733 ...                              0.013127131925856725\n",
       " 1417  3848.017315598605 ...                               0.09620434270592455\n",
       " 1419  3723.004118947109 ...                               0.08801476817405915\n",
       " 1438 3604.9472263665893 ...                               0.06594961678474937\n",
       " 1382  3932.032866076106 ...                               0.09575037644968758\n",
       "  ...                ... ...                          ...                  ...\n",
       "  107  182.0508911360561 ...                               0.26216596171042117\n",
       "   39  541.0001708205461 ...                                0.3810401372674711\n",
       "   97 193.00622531012922 ... ASASSN-V J205607.11+312038.1    0.126141972104651\n",
       "   41  503.9609772205891 ...                               0.07142793901418203\n",
       "   14  602.9782312415216 ...                                1.0873130222307894\n",
       "   80  264.0154239709242 ...                               0.26812608072554045\n",
       "   98  83.04174328655962 ...                               0.17867986111223463\n",
       "   18 306.91673956570503 ...                                0.1442924276151822\n",
       "   37 102.38009225045143 ...                              0.004799602972304384\n",
       "   10 113.92856473017127 ...                               0.33344694042018297"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mags"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5,0,'t_mid_tdb')"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAELCAYAAADz6wBxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xu4JHV95/H3Zy7CGYzMKJPEOQQHbxAukQknihB3uUUIeEEIATZ4YXEJSVyUKBFiVtAky7jgg6smcYlBjBKCIs6yQh7Gx8GQ8ABxhhmcQUARBOeMWcaFMVFGnMt3/6hq6Omp6q7uququPufzep7zTJ/qqurfqen6fet3V0RgZmY2qDmjToCZmY03BxIzMyvFgcTMzEpxIDEzs1IcSMzMrBQHEjMzK6W2QCLpGklPSNrQtu10SfdL2ilpqsuxF6b7bZB0vaQ90+37S7pH0nck3SDpeXWl38zMiqmzRHItcGLHtg3AqcAdeQdJmgQuAKYi4hBgLnBm+vZHgKsi4hXAU8C5FafZzMz6VFsgiYg7gCc7tj0QEQ8VOHweMCFpHrAA2CRJwLHAjek+nwVOqTDJZmY2gHmjTkCniJiWdCXwOLAVWBkRKyXtA2yJiO3prhuBybzzSDoPOA9gr732OvzAAw+sOeVmZjPLmjVrfhgRi3vt17hAImkR8GZgf2AL8EVJZwO3ZeyeO79LRFwNXA0wNTUVq1evriG1ZmYzl6THiuzXxF5bxwOPRsTmiNgG3AQcCfwQWJhWdwHsC2waURrNzCzVxEDyOHCEpAVpu8hxwAORzC55O/Bb6X5vB/73iNJoZmapOrv/Xg/cBRwgaaOkcyW9RdJG4LXALZJuS/ddIulWgIi4h6RB/V5gfZrGq9PTvh/4Q0kPAy8C/qau9JuZWTGaDdPIu43EzKx/ktZERO6Yv5YmVm2ZmdkYcSAxM7NSHEjMzKwUBxIzMyvFgcTMzEpxIDEzs1IaN0WKmdlss2LtNFfc9hCbtmxlycIJLjrhAE5ZljuVYOM4kJiZjdCKtdNcctN6tm7bAcD0lq1cctN6gLEJJq7aMjMboStue+jZINKyddsOrrityIobzeBAYmY2Qpu2bO1rexM5kJiZjdCShRN9bW8iBxIzsxG66IQDmJg/d5dtE/PnctEJB4woRf1zY7uZ2Qi1GtTda8vMzAZ2yrLJsQocnVy1ZWZmpTiQmJlZKQ4kZmZWigOJmZmV4kBiZmalOJCYmVkpDiRmZlaKA4mZmZXiQGJmZqU4kJiZWSkOJGZmVooDiZmZleJJGxtq3NdwNpspfC/25kDSQDNhDWezphkkIPheLMZVWw00E9ZwNmuSVkCY3rKV4LmAsGLtdNfjfC8W40DSQDNhDWezJhk0IPheLMZVWw20ZOEE0xlf1G5rOLse1yzfoAFhkHtxNnKJpGYr1k5z1PJV7H/xLRy1fFXPojT0v4bzoMV2s9kiL+Nv3551r86E9dSHwYGkRoNm8Kcsm+TyUw9lcuEEAiYXTnD5qYfmljBcj2vWXa+AkHevAn3di7OVq7Zq1C2D7/VF7GcNZ9fjmnXXupfyqn+73at3XnxsbYFjplRJO5DUaFgZvOtxzXrr9nA27IexFWunuezm+9mydduz28a5a3FtVVuSrpH0hKQNbdtOl3S/pJ2Sproce2G63wZJ10vaM93+LkkPSwpJ+9SV9qoUqZetgutxzcqp+15tb3857EMruejG+3YJIi3jWiVdZxvJtcCJHds2AKcCd+QdJGkSuACYiohDgLnAmenbdwLHA49Vndg6DCuD77dNxcx2Vee92tn+smXrNrbtiNz9x7FKuraqrYi4Q9LSjm0PAEjqdfg8YELSNmABsCk9fm3B4xshq172mAMXc8VtD3HhDesqrRPtp03FzHbVqw2ljKz2l27GsUq6cW0kETEt6UrgcWArsDIiVvZ7HknnAecB7LffftUmsg/tGbynWzBrrn4fxoo2lPdTwhjXKunGdf+VtAh4M7A/sATYS9LZ/Z4nIq6OiKmImFq8eHHVyRyIu+mazQz9dO0vWsJYtGD+2FZJNy6QkLSBPBoRmyNiG3ATcOSI01SJvCeT6S1bWdrHgEUz622QwcBF9fNQmNX+Mn+OWLRg/rNtmh874zDWfvD1YxlEoIFVWyRVWkdIWkBStXUcsHq0SepfVrE3r5tuSz9VXTOl/7lZHequRu72ULhi7fQun1Fn+0tTKCK/90CpE0vXA0cD+wD/F7gUeBL4BLAY2AKsi4gTJC0BPh0RJ6XHfgg4A9gOrAXeGRHPSLoA+CPgF4EngFsj4p290jI1NRWrVw8vFnV+iSGp+zzt8Em+tGa6Z8Pb5MIJ7rz42L7PP67FYrOqHbV8VeZDW697q+z5YWbdi5LWRETuUI2W2qq2IuKsiHhxRMyPiH0j4m8i4svp6z0i4hci4oR0302tIJL+fmlEHBgRh0TEWyPimXT7x9Pj50XEkiJBZBTyir23P7iZy089tOfxvRrn3NZi1l3dAwyzqqtaZuO92MQ2krHX7Ut8yrJJJns0vvVqnPOUKGbd5d1Dc6Td2koGaUtpjd3KM9vuRQeSGvQaJdvtaaZI979hjZg3G1d599iOiF16V5WZObvbQ2HAQA38dXYQqJMDSQ16jZJtH4kOMDcdYFl0RLqnRDHrrnWPzc0YvNxe9VS2mrjbQ2G/yzmM83IQTey1NfaK9NIoMxK9bC8Q9/iy2eCUZZNceMO6zPdaVU9lq4nb78WsxvfO2b673XtlZgsfNQeSmtQ9Zcmg5/foeptNes2MXcXM2a17cf+LbyGrD2wrKPW698a57dNVWw1Rdd1o3vnc48tmk17VwFVWE/dqu+x1741z26cDSQNUXTfa7Xzj/NRj1q9eM2NXOXN2r6DU694b57ZPV201QNV1o93O50WwbLbpVQ1cVTV0r7bLXvdeZ3vLXGmXEkuTq54dSBqg6lJCt/NddcZhmaPiLzrhADfCm5XULShddMIBufde+/HA2LVjumqrAaquG+12vryiPDC2XQ/NxkHRarRxbMd0iaQBijypVHm+rKemo5avGtuuh2bjokg12ji2YzqQNEDVs4MOcr5eX15Xe5kNxzi2YzqQNETV4076PV+3L6/HnpgNZpAHsKprKIbBbSQNN6y5d7p1PRzHOluzURu0W3+VXZKHxSWSBhtmSaBbdVivaSbMbHdluvXXPTNG1RxIGizvi/jeL9zHhTesG6itoltRO+/LO451tmZ5htHet2LtdO7CVzPxAcxVWw2W94XbETFQF91Bi9rjPOLWrN0wZthtfUaemfgA5kDSYEW+cP20VQza1jGOdbY2+xRpTxxGe1/WZ7TM1AcwV201WFbvjSxFi8pl+qePW52tzS5F2xOHMUaj27lm6gOYSyQN1lkSyFqkB5LlQ4v06hrn2UXNuila0qjiHuhV8sk712Q6s8RM5EDScKcsm+TOi4/l0eUn89HfflXu8qFF6nvd1mEzVdGSRtl7oEgby2y8zxxIhqSK8SBFSijd6nvd1mEzVdGSRpl7YMXaad77hft6lnxm432miKw1vWaWqampWL169cg+v7P+FpInlLJfrrwV2QQ8uvzkgc9rNg7au/HuPTGfn/xsO9t2PHdHVHGPtX9Wt/bKmXrPSVoTEVO99nOJZAjq6iniNg+brTqrmLZs3QYBixbMr6UU0K0nFviec6+tIairp0hVc/J4QkYbN1kZ+7adwYLnzWPtB19f+ed1u1dnevtHEQ4kQ1DXyPAqZg32hIw2jqp4OMt7gMranncPz5VmfPtHEQ4kQ1DnbJ5lx3dUvcyvWT8GLQ2XfTjLe4Ba/diTfGnN9G7bTzt8cpftUG0bzLhzG8kQNLkXxzguomMzQ5npSsp2sc17gLr+nu9nbr/9wc2NvYebwCWSIWnqyHBPyGijUnZ23NY5BqnW7TaPXd7+Tb2Hm8AlklluNg6esmYoWxpuDda96ozDALjwhnWFx2jlPSjlzR7hB6vuHEhmuSZXu9nM0z4wd04FmXbVM1qf9Zpf8oPVAFy1ZS6y21B0NnBnVSP1m2kPWj3WrWps6iUvdHf4PjmQzHAeI2JNkTeob67EzoiBvp91zGjtB6v+OZCMiUECgseIWJPkZe47IwaeXsSdRZqhtjYSSddIekLShrZtp0u6X9JOSbnzt0i6MN1vg6TrJe2Zbr9O0kPp9mskza8r/U3STz1wex10kQnmzIaljil93FmkGepsbL8WOLFj2wbgVOCOvIMkTQIXAFMRcQgwFzgzffs64EDgUGACeGe1SW6monN1dQacbl0ZzYatjkzfnUWaobaqrYi4Q9LSjm0PACint0abecCEpG3AAmBTevytrR0k/Quwb3Upbq6i9cC9JpZrmSOxYu20bzYbqiqm9Mk7r7/Lo9W4NpKImJZ0JfA4sBVYGREr2/dJq7TeCrw77zySzgPOA9hvv/3qS/AQFK0HLlrS2BHhthIbCWf6M1PjxpFIWgS8GdgfWALsJensjt3+ErgjIv4p7zwRcXVETEXE1OLFi+tL8BAUrRLop67ZbSU2rqpYJM6q1bgSCXA88GhEbAaQdBNwJPD59PdLgcXA744shTXq1jurV5VA3uSQedVd01u2ctTyVe4abGPDPRGbqYmB5HHgCEkLSKq2jgNWA0h6J3ACcFxE7BxdEuvR6yYZdP6hK257KLNqTOlndH5W1jl8k1oTeLbqZqotkEi6Hjga2EfSRuBS4EngEyQlilskrYuIEyQtAT4dESdFxD2SbgTuBbYDa4Gr09N+CngMuCttsL8pIj5c198wbFXcJHkBp7OkIthtmd6t23Zw2c3388z2nX7is0bybNXNVGevrbNy3vpyxr6bgJPafr+UJPB07tfEElRl6rpJskoqWSUUSJcs7VDXE59H3Vu/hjkA0d/P4mZ0xjxu6rxJOksqRy1flRtMslT9xFdHXbdv/Jkj7/+yzkXiOj/fbTHFNa7X1mw2zFG6eZ+1aEH2ZAFVP/EVHWRZVNbo/4u+eB/LPrzSvXsaqFvPq24zOQxrAGLV38+ZziWSBqlrwFY/nwW7t6fUEcyqrsbLuvG37QyeejqpqpvNT5RNK6n1etrv1VY4jLEobovpjwNJwwxzwFa3z6o746m6Gq/IDV5n756mZdbt6WpaFU2vQNGETNyTQfbHgcR2M4xgVnVdd7cOBO3yMqMygaCJmXVLE7vL9goUTcjEh9UWM1O4jcRGouq67qw2nyxZmdGgq+y1NLk+vQlP9516zQLchBl9+/1+zvbR9i6R2MhUWfLpbPPZe2I+P/nZdrbteG60TF5mVPapvYmZdUsTnu479XraH2ZbYTdFv59NLpEOiwOJzRidN37R6qqygaCJmXVLE6toigSKcZrcsYnVh8NWKJBI+tWMzT8CHouI7dUmycZBUxuX2xXNjLoFgiJ/ZxMz65amPN1npWvUaahKk0ukw1K0RPKXwK8C3ySZXeOQ9PWLJJ3fOc27zWwzrSifFwiOOXBxob+zqZl1y0zKtIsY9kNOk0ukw1I0kHwPODci7geQdBBwEfCnwE2AA8kMlHdDzrSifLfJLov+nU3IrMehlFhGkb9vFA85TS6RDkvRQHJgK4gARMS3JC2LiEcKrHZoY6jbDTkTi/JZgeDCG9Zl7tvEv3OmlRI7Ff37RvGQ0/QS6TAUDSQPSfor4O/T388Avi1pD2D3Wf5s7HW7IWdLUX6c/s6ZVkrsVPTvG9VDThNKpKNUdBzJO4CHgfcAFwKPpNu2AcfUkTAbrW43ZBP6+Q/DOP2dM7GU2K7o39drjIrVo1CJJCK2Ah9Nfzr9uNIUWSN0exqfLUX5cfo7+yk9NbEtpVeaiv59bq8YDUV0Lm+UsZP0CuBy4CBgz9b2iHhpfUmrztTUVKxevXrUyRgrnXXSkNyQdcy0auVl/X8BLFown0vfePCz/2dN/H8tkqZ+0t3EQDmuJK2JiKle+xVtI/kMyUJTV5FUZZ1D0g3YZqhxehpvsmFlaq1zXnbz/bssTvbU09v6mll3FIqkqZ/v42xvrxiFooFkIiK+JkkR8RhwmaR/ImMVQ5s5fEOWM+yeVK1A0bnKZdNm1i362Zu2bHXpYkwUbWz/qaQ5wHckvUvSW4CfrzFdZmNvFJM5FplZN0v79mFPQJiXpr0n5peaTNOGp2iJ5D3AAuACkkGIxwBvqytRNj7qeGJsnXN6y1bmSuyIYLLic9f1hNt+/rzWxzqf/ns1SvdqjG7SgD6JxlXDWbaiJZIAPgfcDEwBrwT+uq5E2XgoO/16r3MC7Eg7g3Q7d9En6DrS2+38eersitqry3Kv6dFHUYrKS9OWp7OHqM2ULs0zSdESyXUkU6KsB3bWlxwbJ3U03Gads9u5+3mCrruhuVvaW+ruilp2Zt0mDehrlUo7eUxI8xQNJJsj4uZaU2Jjp45Mp9exne/3ExzqziS7nUdQedVfHVOwlx3NX2XVoceEjI+igeRSSZ8GvgY809oYETfVkiobC3mZTgBHLV81UCbSa8nczgytn+BQ95QneeefXDjBnRcfW8ln1N2GkZV5C1j6ogmOWr5qqBMmugv6+CjaRnIOcBhwIvDG9OcNdSXKxkO35W0HbX/ods6sp9F+psSoe8qTYUypUncbxinLJjnt8MldBokFcOd3n+zZtlRH2k5ZNsmdFx/Lo8tP5s6Lj3UQaaiiJZJXRcShtabExk77E2PWk/gg7Q+d5+zVa6tI9Ud7dcveE/PZc/4ctjy9rfIn3GE8QQ+jDeP2Bzd37SwAzZow0UavaCC5W9JBEfGtWlNjY6dVH7/04lsy3x8kE+mnjr9X5t1Z3bJl6zYm5s/lqjMOq21Q4LgvolT0/6x9vxVrp5mTBv0602bNVDSQ/DrwdkmPkrSRCIiI+JXaUmZjY8Xa6eQLkfFeFZlImcblqnpqNWWE9TAaoHu1U7XMkZ6t3rrkpvWZQcSN47ND0UByYq2psLF2xW0PZQYRQelMpGwDbhXVLYOkoa7AU0f1WWdajzlwMV9aM92zK/OOCC65aT17zp+Tue9cyZN8zhJFp5F/rO6E2PjKy5SD8j2JypYoqqgK6jcNdfes6iyBtQZkDhJYstL6pTXTnHb4JLc/uHmX4HL9Pd/frdSxdduO3ICzM9130LTZ+ChaIjHL1a3ba1llSxRVVAX1m4ZhzrBbNmjlpfX2Bzfv1mX5ursf7yttrbmyZuryv/acot1/zXLV2e217Ip3vaYEqSMNw+y91E+X26ypZPodh5Nl4cT8zP//bnNl2cziQGKlVZFZ56kiSPUzFiErs+03DcNc7rVoIMibZ2zhgvmZx/czDueyNx3c11xZ0+n08DZzuGrLKlFXt9dhjm7Oqya6/NRDufzUQwunYZhTexRtA8oruewxbw4T8+cWSmuv/4uic2UBuVVcTekdZ/0ptNTuuPNSu1bEUctXVTbFybAyxKJL0O5/8S25PeuuOuOwWqbtz1v+t6XzujZxGeDZruqldgdJwDUk06g8ERGHpNtOBy4Dfhl4dURk5u6SLgTeSdLxZz1wTkT8VNLfkExjL+DbwDsi4sd1/Q02u1TZtjGs1SWLlti6lVxa+1bdMN467j03rMt8v8wEnNYsdbaRXMvu4082AKcCd+QdJGmSZAGtqTQAzQXOTN++MCJelQ6EfBx4V9WJttlrkLaNYa8mmKVIG1Cvdp665vA6Zdlkbu+9MhNwWrPUViKJiDskLe3Y9gCApKxD2s0DJiRtI1mZcVN6/L+lxwuYIHswtTVU1SsfVl191G/bxihWExxUr5JLnZl40es6jOlfrB6Na2yPiGlJV5KUOLYCKyNiZet9SZ8BTgK+Bbw37zySzgPOA9hvv/1qTbPtKiuDh12rTjpXPoT+Mt86MvF+GvZXrJ3mvV+4L3OAXlOrYrpVtw2aiRcJ5kWvq9cfGV+1NranJZKvtNpI2rZ/HXhfVhuJpEXAl4AzgC3AF4EbI+LzbfvMBT4BfCMiPtMrHW5sH568BtM958/hqZzuoMCzJZOiJYxuDeP9nGcQvRqRBTy6/ORKPmdYPZgGaeiuo3HcvbaaZeSN7SUcDzwaEZsBJN0EHAk8G0giYoekG0iW/+0ZSGx48urae83b1CpRFC1h5FW59HueQfRaUreqiSqHWW02SDfrOhrH80pNDjDN1sRA8jhwhKQFJFVbxwGr03aRl0XEw+nrNwIPjjCdlmHQOvW5Ul+ZUl5VTL/n6deKtdNdZ8atoiqmW7XZe79wH1BfMOnnvMNqHB+ntqjZqrZeW5KuB+4CDpC0UdK5kt4iaSPwWuAWSbel+y6RdCtARNwD3AjcS9L1dw5wNUmNwWclrU+3vxj4cF3pt8H0M41Gy8T8uZlTkEN+ppTXC6nf8+TJ6o3VytDyVDHbbesz8v6O1oy7TRgZPqwR/HWvCmnl1dlr66yct76cse8mkgb01u+XApdmHHtUNamzuuQ1mF72poOB/JUP80ZB52VKeVUx/Z4nS94TcN506a2/sUgQ6VVF06vaDIbboN8tvcNqHHe34OZrYtWWjbF+p9GAJLP6yTPbd9uelyl1Zm6dqx2WzdwGaecpGkR6VdEMsjphXXqlt9v/dZVtGu4W3HwOJFa5fura83pALVown0vfePBu625cdvP9bNn6XO+vfjK3ovrNpCfbRofnKdpduOjqhMNYebJIY3rW/3XVbRruFtx8DiQ2UnlVOQueN2+3IJLX5bZo5lb0yTkvM184MZ9ntu/sO0Pr1e7RHriyMs0sT/9sOyvWTvecybjoWvb9lI56Bdqqe3MNc+JOG4wDidWmSPVG0cyqV9tBt8wtL9Nc/diTuywp29p+2uGTuy0129nO00+G1k934axM85gDF/OV+36wS0nsqae3dX3K7xUoimT2g1Yp1dGmMay5y2wwDiTWl6J130WrN4pmVr0yoW6ZW16mmbd07O0Pbu46bXy/GVq3tGeVaLIyzdsf3LxLIGmlNe8pv1egKJLZD1ql5DaN2ccLW1lheYsjZXVFLdpls+iiUd0yoV6ZW16m2a2qqZ/FsHrJS3s/3YX7fcrvtb1I191BFyyrc8VMayYHEiusn/78RTO4oplVVuYESaN8r8ytW0bez/6DystYP/rbryocoPods9Fre9HMfpCAOmgAsvHlqi0rrN/1vYtWbxSp/y7T4JpXRZPXFlL1k3MVjcX9VjP12r9smnpVcbpNY3ZxILHC+gkOdXTZHDRz6pZpTr3khUPpDVQ2Y+034y+y/6Bp8pQl1slL7Vph/c726on2ZqYqlyS2Zhvn2X+toQZ5KnbgqM+oArWnLLFODiTWFweH+vQTGEZZveTuvdbJvbbMGqCfrtUw2hlx3b3XOrlEYjYCnaWPp3+2va9pRUZZveQpS6yTA4nZkGVVS+XJCwyjrl5yFae1c9WW2ZAVWXOkJS8wuHrJmsQlErMhK1r91C0wuHrJmsSBxGrn8SS76jZN/V57zCt8nVy9ZE3hQGK18ijo3XVbjni2XhMbb24jsVqNsptqU3lSQ5tpXCKxWnkUdDZXS9lM4hKJ1arfac7NbPw4kFit3E3VbOZz1ZbVyt1UzWY+BxKrndsDunP3aBt3DiRmI+Tu0TYTOJDYrNSUUkC37tEOJDYuHEhs1mlKKWDF2uncCRtne/doGy/utWWzThMGSbaCWR53j7Zx4kBis04TBkl2mwHY3aNt3DiQ2KzThEGS3YKWp0uxceNAYrNOEwZJ5gWtyYUTDiI2dhxIbNZpwqSJTQhmZlVxry2blUY9SNIj/m0mqS2QSLoGeAPwREQckm47HbgM+GXg1RGxOufYC4F3AgGsB86JiJ+2vf+JdNvz60q/Wd1GHczMqlJn1da1wIkd2zYApwJ35B0kaRK4AJhKA9Bc4My296eAhVUn1szMBlNbIImIO4AnO7Y9EBFFOuvPAyYkzQMWAJsAJM0FrgD+qOLkmpnZgBrX2B4R08CVwOPAD4AfRcTK9O13ATdHxA96nUfSeZJWS1q9efPm+hJsZjbLNS6QSFoEvBnYH1gC7CXpbElLgNOBTxQ5T0RcHRFTETG1ePHi+hJsZjbLNbHX1vHAoxGxGUDSTcCRwFPAy4GHJQEskPRwRLx8ZCk1M7NGBpLHgSMkLQC2AscBqyPiFuAXWztJ+rGDiJnZ6NVWtSXpeuAu4ABJGyWdK+ktkjYCrwVukXRbuu8SSbcCRMQ9wI3AvSRdf+cAV9eVTjMzK0cRMeo01G5qaipWr84csmJmZjkkrYmIqV77Na6x3czMxosDiZmZleJAYmZmpTiQmJlZKQ4kZmZWigOJmZmV4kBiZmalOJCYmVkpDiRmZlaKA4mZmZXiQGJmZqU4kJiZWSkOJGZmVooDiZmZleJAYmZmpTiQmJlZKQ4kZmZWigOJmZmV4kBiZmalOJCYmVkpDiRmZlaKA4mZmZXiQGJmZqU4kJiZWSkOJGZmVooDiZmZleJAYmZmpTiQmJlZKQ4kZmZWigOJmZmV4kBiZmalOJCYmVkpDiRmZlaKA4mZmZXiQGJmZqXUFkgkXSPpCUkb2radLul+STslTXU59sJ0vw2Srpe0Z7r9WkmPSlqX/hxWV/rNzKyYOksk1wIndmzbAJwK3JF3kKRJ4AJgKiIOAeYCZ7btclFEHJb+rKs2yWZm1q95dZ04Iu6QtLRj2wMAknodPg+YkLQNWABsqiGJZmZWgdoCyaAiYlrSlcDjwFZgZUSsbNvlzyV9EPgacHFEPJN1HknnAeelv/5Y0kN1prvDPsAPh/h5VXG6h2cc0wxO97CNOt0vKbJT4wKJpEXAm4H9gS3AFyWdHRGfBy4B/hV4HnA18H7gw1nniYir032GTtLqiMhtA2oqp3t4xjHN4HQP27iku4m9to4HHo2IzRGxDbgJOBIgIn4QiWeAzwCvHmE6zcyMZgaSx4EjJC1Q0phyHNBqW3lx+q+AU0ga783MbITq7P57PXAXcICkjZLOlfQWSRuB1wK3SLot3XeJpFsBIuIe4EbgXmB9msZWFdV1ktan2/cB/qyu9Jc0kiq1CjjdwzOOaQane9jGIt2KiFGnwczMxlgTq7bMzGyMOJCYmVkpDiQ9SDpR0kOSHpZ0ccb7V7VN2fJtSVvS7YdJuiud6uWbks5oO6b2qV69ixvxAAAKDUlEQVQGTXf63o62925u276/pHskfUfSDZKe15R0Szqmbfs6ST+VdEr6XhOu936Sbpe0Nv0+nNT23iXpcQ9JOqHoOUeVZkm/IWmNpPXpv8e2HfP19Jyta/3zDUr3Uklb29L2qbZjDk//noclfTzt0NOUdP9Ox3d7Z+s7PIzrXUhE+Cfnh2R6lu8CLyUZu3IfcFCX/f8rcE36+pXAK9LXS4AfAAvT368FfquJ6U5//3HOfl8Azkxffwr4vSalu237C4EngQVNud4kjaa/l74+CPhe2+v7gD1Ixk59Nz1fX9diyGleBixJXx8CTLcd83WS6Y2aeK2XAhtyzvsvJJ2ABPwD8JtNSXfHPocCjwzrehf9cYmku1cDD0fEIxHxM+DvSQZL5jkLuB4gIr4dEd9JX28CngAW15zeloHTnSd9QjuWpEcdwGdJumBXqap0/xbwDxHxdMXpy1Mk3QG8IH29N89N+/Nm4O8j4pmIeBR4OD1fv9diaGmOiLXpdxrgfmBPSXtUmLZuylzrTEqGFbwgIu6KJHf+W0bz3S6S7p736ig4kHQ3CXy/7feN6bbdSHoJyRPlqoz3Xk3yFPLdts1/nhZfr6rhJiyb7j0lrZZ0d6t6CHgRsCUitvc6ZwmVXG+SST47b7ZRX+/LgLOVdH+/laQ01e3YwtdiQGXS3O40YG3sOlXRZ9Jqlv9WQxVR2XTvn1Yd/aOk17Wdc2OPc5ZV1fU+g92/23Ve70IcSLrL+k/J6y99JnBjROzY5QTJ087ngHMiYme6+RLgQODXSKph3l9Ncp/72Ixt/aR7v0imZfhPwMckvazPcw6qqut9KHBb2+YmXO+zgGsjYl/gJOBzkuZ0Obbu610mzckJpIOBjwC/23bM70TEocDr0p+3VphmKJfuH5B8t5cBfwj8naQXFDxnWVVc79cAT0dE+0Dsuq93IQ4k3W0Efqnt933JLybv9hScfklvAf4kIu5ubY/6p3ople5WtUVEPEJSB7uMZOK4hZJa87N1O+egSqU79dvAlyOZXgdozPU+l6SNiYi4C9iTZFBt3rH9XIthpxlJ+wJfBt4WEc+WtCNiOv3334G/o0HXOq0+/H/p9jUkNQSvTM+5b49zjizdbe9n3at1X+9iRt1I0+QfkkktHyGpQmk1kB2csd8BwPdIB3im255HMkPxezL2f3H6r4CPAcsblO5FwB7p632A75A2CgJfZNfG9t9vSrrb3rsbOKZp15ukAfcd6etfJslEBBzMro3tj5A0zBa6FiNK88J0/9MyzrlP+no+SXva+Q261ouBuen2lwLTwAvT378BHMFzje0nNSXd6e9zSILRS4d5vQv/faP40HH6ISlifpvk6eUD6bYPA29q2+eyzswJOBvYBqxr+zksfW8VyTQvG4DPA89vULqPTNN2X/rvuW3vvZSkd8vDJEFlj6akO92+NM0c5nRsH/n1JumFc2d6XdcBr2879gPpcQ/R1lso65xNSDPwJ8BPOr7bPw/sBawBvknSCP8/STPuhqT7tDRd95FMwfTGtnNOpd+P7wKfJOMhZcTfkaOBuzvON5TrXeTHU6SYmVkpbiMxM7NSHEjMzKwUBxIzMyvFgcTMzEpxIDEzK0DSFZIeTGdI+LKkhV32nZuOoP9K27bMyUMl7S3p/0i6T8kkr+ek21+SToq5Lt1+frp9gaRb0rTcL2l522ecn04+uU7SP0s6qMDftZ+klZIekPQtSUv7vjbutWVmtitJR5OM6XhH27bXA6siYrukjwBEROYsCZL+kKRL8Qsi4g3ptmuBr0TEjR37/jGwd0S8X9Jikm7gv9h6OyKekfR8ku7JRwJbgNdExO1KZuD+GvDfI+IfJL0gIv4tPe+bSMZ6ndjjb/068OcR8dX0c3ZGn/PUuURi1oWkhZJ+v4LznC/pbRnbl0rakHVM+v5haptyPuP970nap9d5rLyIWBnPzTV3N7uOhn9WOur/ZODTRU8N/Fw6T9bzSWau3h4RP4vn5jDbgzS/joinI+L29PXPSMbE7Jv+/m9t590rPXerhHSFpG+kJarfTbcfBMyLiK+mx/+43yACDiRmvSwESgeSiPhURPztAIceRjKQzZrlP5OMRM/yMeCPgJ0Z72VNHvpJnhvJvh54d6Tz8kn6JUnfJJnw8SPx3KzLpO8vBN5IUippbfsDSd8F/gdwQbr5XOBHEfFrJHPO/RdJ+5NMEbNF0k1pVdwVkub2dykcSMx6WQ68LK1zvqLzTUlHpzPJfkHJQlvLlSxE9C9pXfXL0v0uk/S+9PXhaX34XcAf5H1wWm3xYeCM9PPPkPSitD57raT/xa6TAc6T9Nk0o7pR0oIqL8RsoGThtnUkpYk3tbVntC849gFgO3BdxvFvAJ6IZC6vTnmTh55AMpJ9CcmDwyfTefqIiO9HxK8ALwfeLukX2j5rHsncWx+PZF480mP+IiJelp7/T9LNrwfelv5t95DM5v0KkmlWXge8L03XS4F3FLxczxnFcHr/+GdcfuiyGFL6/tEkddYvJql+mAY+lL73buBj6evLgPelr78J/Mf09RU9zv8O4JNtv38c+GD6+mSSqot90nQGcFT63jWtz/PPQP/vR5PMxNu5/e3AXaSLpmW8fznJnFjfA/4VeBr4fM75v5K+vgV4Xdt7q4BXZxzzGdoWaEv/jz/e5W+YQ1IKAfgScELGPkcAX2/7/a3AX/R7vVwiMSvvG5HMMPwMyTxKK9Pt60ky+GdJ2ptkpcx/TDd9rs/P+g8k84UREbcAT7W99/2IuDN9/Xng1/s8t3Uh6USSp/w3RU47QkRcEhH7RsRSktl6V0XE2enxL07/FcnCWa02rceB49L3foFkUtJHJO0raSLdvgg4iqQhHkl/RrL41Xs60viKtl9PJpl0FZJlFX5P0vx0v1dK2otksspFaSM/JIvXfavPS8O83ruYWQ/tizrtbPt9J7vfY6L8Whd5x3dud5fMan2SpNT51SQWcHdEnC9pCfDpiOjVlnVdmmGLpCrr/HT7nwLXSlqfvvf+iPihpN8APiqptT7NlRGxPm3M/wDwIHBvmpZPRsSngXdJOp5kwtinSEpQkFTVLU33F7AZOCUidqRVrl9Lt68B/rrfC+NAYtbdvwM/V9XJImKLpB9J+vWI+Gfgd/r8/DvSY/5M0m+STPvfsp+k10aylsVZwD9Xle7ZJiK+TrIWT/u2l+fsu4mMDhGd54iIY7sc//qM7V8FfiVj+0ayF8oiIt6ds30n8MfpT6HP6Yertsy6iGQhpDslbchqbB/QOcBfpI3tW3vseztwUKuxHfgQ8B8k3UuS+Tzetu8DJA2y3yRpzP2ritJr1pUHJJqZWSkukZiZWSluIzErQNKh7N7D6pmIeE1F5z8B+EjH5kcj4i1VnN+sTq7aMjOzUly1ZWZmpTiQmJlZKQ4kZmZWigOJmZmV8v8B2F+SC6Dt9GAAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "this_star = mags['star_id'] == 802\n",
    "plt.scatter(mags['t_mid_tdb'][this_star][0:140], mags['mag_cal'][this_star][0:140])\n",
    "plt.ylim(11.8, 11.85)\n",
    "plt.ylim(reversed(plt.ylim()))\n",
    "plt.ylabel('mag')\n",
    "plt.xlabel('t_mid_tdb')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}