extract_reference_data.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import sys\n",
    "sys.path.append(\"/Users/papr/work/pupil/pupil_src/shared_modules\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import seaborn as sns\n",
    "\n",
    "import file_methods as fm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "datapath = \"/Volumes/cluster/recordings/reference_recordings/5-vs-10-point-2d-calib_v2.2\"\n",
    "notification_data = fm.load_pldata_file(datapath, \"notify\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "def from_ref_data(ref_data):\n",
    "    col_index = pd.MultiIndex.from_product(\n",
    "        [[\"normalized\", \"pixels\"], [\"x\", \"y\"]], names=[\"unit\", \"coordinate\"]\n",
    "    )\n",
    "    return pd.DataFrame(\n",
    "        [[*ref[\"norm_pos\"], *ref[\"screen_pos\"]] for ref in ref_data], columns=col_index\n",
    "    ).stack(level=[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [],
   "source": [
    "ref_data = [\n",
    "    from_ref_data(calib_data[\"ref_list\"])\n",
    "    for calib_data, ts, topic in zip(*notification_data)\n",
    "    if topic == \"notify.calibration.calibration_data\"\n",
    "]\n",
    "ref_data = pd.concat(\n",
    "    ref_data, keys=notification_data.timestamps, names=[\"timestamp\"]\n",
    ")\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": "<div>\n<style scoped>\n    .dataframe tbody tr th:only-of-type {\n        vertical-align: middle;\n    }\n\n    .dataframe tbody tr th {\n        vertical-align: top;\n    }\n\n    .dataframe thead th {\n        text-align: right;\n    }\n</style>\n<table border=\"1\" class=\"dataframe\">\n  <thead>\n    <tr style=\"text-align: right;\">\n      <th></th>\n      <th></th>\n      <th>coordinate</th>\n      <th>x</th>\n      <th>y</th>\n    </tr>\n    <tr>\n      <th>timestamp</th>\n      <th></th>\n      <th>unit</th>\n      <th></th>\n      <th></th>\n    </tr>\n  </thead>\n  <tbody>\n    <tr>\n      <th rowspan=\"5\" valign=\"top\">36.722786</th>\n      <th rowspan=\"2\" valign=\"top\">0</th>\n      <th>normalized</th>\n      <td>0.424593</td>\n      <td>0.255945</td>\n    </tr>\n    <tr>\n      <th>pixels</th>\n      <td>543.479258</td>\n      <td>535.719555</td>\n    </tr>\n    <tr>\n      <th rowspan=\"2\" valign=\"top\">1</th>\n      <th>normalized</th>\n      <td>0.426099</td>\n      <td>0.256685</td>\n    </tr>\n    <tr>\n      <th>pixels</th>\n      <td>545.407200</td>\n      <td>535.186604</td>\n    </tr>\n    <tr>\n      <th>2</th>\n      <th>normalized</th>\n      <td>0.426993</td>\n      <td>0.257668</td>\n    </tr>\n    <tr>\n      <th>...</th>\n      <th>...</th>\n      <th>...</th>\n      <td>...</td>\n      <td>...</td>\n    </tr>\n    <tr>\n      <th rowspan=\"5\" valign=\"top\">269735.971200</th>\n      <th>387</th>\n      <th>pixels</th>\n      <td>746.284958</td>\n      <td>275.000000</td>\n    </tr>\n    <tr>\n      <th rowspan=\"2\" valign=\"top\">388</th>\n      <th>normalized</th>\n      <td>0.583721</td>\n      <td>0.618056</td>\n    </tr>\n    <tr>\n      <th>pixels</th>\n      <td>747.162724</td>\n      <td>275.000000</td>\n    </tr>\n    <tr>\n      <th rowspan=\"2\" valign=\"top\">389</th>\n      <th>normalized</th>\n      <td>0.583836</td>\n      <td>0.617982</td>\n    </tr>\n    <tr>\n      <th>pixels</th>\n      <td>747.310480</td>\n      <td>275.053028</td>\n    </tr>\n  </tbody>\n</table>\n<p>1170 rows × 2 columns</p>\n</div>",
      "text/plain": "coordinate                             x           y\ntimestamp         unit                              \n36.722786     0   normalized    0.424593    0.255945\n                  pixels      543.479258  535.719555\n              1   normalized    0.426099    0.256685\n                  pixels      545.407200  535.186604\n              2   normalized    0.426993    0.257668\n...                                  ...         ...\n269735.971200 387 pixels      746.284958  275.000000\n              388 normalized    0.583721    0.618056\n                  pixels      747.162724  275.000000\n              389 normalized    0.583836    0.617982\n                  pixels      747.310480  275.053028\n\n[1170 rows x 2 columns]"
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ref_data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": "<seaborn.axisgrid.FacetGrid at 0x12f2c6430>"
     },
     "execution_count": 25,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAALWCAYAAABiENTGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdebxdZX3o/893nyEnOZmTkxASQhIIKoMgjXEsUqgWqRe0KmJrFUtLbVHppVqx10strffWetVi5dribG9VKP2p0UuxluFiVYYoYQhTYgghgczzSc64v78/9krYOZxMsHf2OSef9+u1X2etZz1rre+e1vnuZz3rWZGZSJIkSXphSo0OQJIkSRoJTKwlSZKkGjCxliRJkmrAxFqSJEmqARNrSZIkqQZMrCVJkqQaMLGWJEmSasDEeoCImBgRf1w1f2xE3FTH/Z0REefXa/u1EBFfjoj7I+KBiLgpIsZWLbsoIh6OiKUR8c1B1h0XEUuqHhsj4u+KZVcW6z4QEbdGxPFF+RkR8bNimw9ExDuqtvfjqm09HRHfLconRMT3iziXRsR7q9bpr1pnUVX5ORHxi4h4KCK+HhHNRfmFxX6XRMTiiHjtgOc0PiJWR8Tnq8reUayzNCI+WVU+OyJuj4j7iuXnVy37aEQsj4jHIuI3BuyjqVjnBweLd5DXfOWgb+TzEBE/Lf7OiYjfrtV2GykiPh4RHyqmr4mIX6/BNne+8MhUCx7Dn8tjuMfwI3EMj4gLIuKq57nuHRGxoNYxNURm+qh6AHOAh47g/i4BPt/o532QGMdXTX8GuKqYng/cB0wq5qcdwrZ+DpxVTP8aMKaY/iPghmL6JGB+MX0s8AwwcZBt/Svw7mL6z4FPFtMdwGagtZjfOci6JeAp4KRi/hrg0mJ6LBDF9EuBRwesey3wzT3vGzAFWAV0FPNfB84tpq8H/qiYPhlYWTV9PzAKmAv8Emiq2seVxT5+cLB4B3luK+vwGTh7TywN/iw212AbHwc+VOO4nvMZ89GYh8fwQWP0GL7vuh7Dh9gDuANY0Og4avGwxfq5/gY4ofil+6niV95DABFxSUR8NyJ+FBErI+L9xS/2+yLiroiYXNQ7ISJuiYifF7/OX1yUv734pXp/RNwZEa1UvlzvKPb3johYWPzSvy8ifhoRLzrMfd8REdcW23soIha+0BckM7cX2w5gNLDndp1/AFyXmVuKeusPtJ2IOAmYBvy4qH97Zu4qFt8FzCrKH8/MZcX008B6Kgfa6m2NB84BvrsnTGBcEeNYKgflvgOEMwXoyczHi/kfAW8t9rkzi2860F71fImIXwGmA/9eta15wLLM3FDM/8eebRXrji+mJwBPF9MXAt/OzO7MfAJYDiws9jEL+E3gS4cS76EqPjdTi+kFEXFHMf3xiPhK8dlZEREfrFpnT0vs3wC/Wnyu/uvh7HeQOO6IiE9GxD0R8XhE/GpR3hYRX42IB4vP9a8V5ZdExKKIuA249TC+C38QEfcW37d/jYgxg8TytYh4W/F67GkRezAisli+v+/y3OJ7+mBE/PULeT1Ucx7DB/AY7jGc2h/Dn/MZLT7jny+mvxcR7y6m/zAi/rmYfkPx/fhFRPxLVJ09KZY3Fcflh4rj6wuKtSEandkPtQcDWjuq56m0TCwHxlE5SGwD3lcs+yzwJ8X0rTz7a/0VwG3F9IPAzGJ6YtU2P1+1v/EUrXLArwP/epj7vgP4YjF9FoO03AAvApbs5/GcVoVina8C64DbebaF4rvA3wI/oXJQPe8gr+3VwP/az7LPAx8bpHwh8AhQGlD+buCmqvlxRWzPADuB36xa1gcsLmJ8c1EWwJMUv5CptGA8WLXOW4BHqRzcX1WUlYrXd1b1+wZMAlYXn5VmKq0w3y+WzSje99XAFuBXqp7vu6r292XgbcX0TcCvUNXCcLB4B7w2K/dXDkwtphcAdxTTHwd+SqXlZSqwCWgplu0s/u6NZZDtjjvA5+nkQerfAXy6mD4f+I9i+k+BrxTTL6bSgtRWvNargcmH+V2YUrXPvwY+UPV8P1RMf23P615V91PApw7yXV7Esy1tl2OL9ZB54DHcY3h6DKf+x/DnfEYHvKbTqXzefxV4HJhcxHYn0F7U+QhwddU2FxSv24+q9jXo53koPwbt36MDuj0zdwA7ImIb8P2i/EHgpcWvr1cD/1L54Q1UPuxQOXh9LSJuBP6//Wx/AvD1iJhP5Zdyy6Huu6retwAy886o9CWbmJlb9yzMzMeAMw7nSWfmeyOiCfh74B1UDtLNVE4lnk3lQHVnRJxWva8BLgZ+d2BhRLyLyhfqdQPKZwD/BLwnM8sDVnsn+7YG/AaVg8A5wAnAjyLix1lpqTk+M9dExDzgtoh4MDN/GREXA5+NiFFUWi/6q57vd4DvRMRZwF9R+Qf5x8DNmbm66r0lM7dExB8BNwBlKge4E6ri/FpmfjoiXgX8U0Scup/Xh4h4E7A+M38eEWdX7SMPFG8N/N/M7Aa6I2I9lYPi6kNZsfhMHtbniWc//z+n8s8M4LVUPl9k5qMR8SSVU8pQOdBurlr/UL4LpxatyROptID98GBBRaUv6JnAGw7yXX4Nz7Y2/RPwSTRceAz3GO4xvMrzPIY/5zM6YJvrIuJqKj+W3pKZm4vX5mTgJ8Xr3wr8bMB2VwDzIuLvgf/LvmcWhgUT68PXXTVdrpovU3k9S8DWzHzOhzQz3xcRr6ByiujnxSmpgf6KysH3LRExh8qvuEPd995dDdx19UxUTk3eMMi+Ac7e30E1M/sj4tvAn1E5KK8G7s7MXuCJiHicykH63oHrRsTpVFpxfj6g/NeB/wa8rjgo7CkfT+VL9d8y864B60yl0grylqri9wJ/k5WfuMsj4gkqrZ73ZOaaIv4VxamzlwG/zMyfUfk1TUS8gWeTuOrnfGdEzCv2+Soqp9L+mEqi1hoROzPzqsz8PsU/yYi4jGcPmJcC5xXb+llEtFH51b4GOK5qV7OKsguAC6JygUwbMD4i/k9mvutQ4j2IPp69YLltwLLqz1Y/h3FsiIhxFKeGB/HbmfnwIOV79neo++rcz/qw/+/C16i0bt0fEZdQSR72q/hn+XEq/Uf7I2K/3+XCwO+ZhgeP4R7DPYZXeZ7H8AN+RgunUWk9P3bPrqg0krxzf7EUP3JOp/JD633ARcDvHSD8Icc+1s+1g8ppkeel+HX9RES8HSp92ooPCRFxQmbenZlXAxuofCkH7m8ClS8nVE6rPB/vKPb3WmBbZm4bEONjmXnGfh77HJCL+E/cM03loPFosfi7FMlKcdA6icqvzcG8k+IXbtW2Xwb8I3BBVvXti0q/xe8A38jMwa7mfxuVU1pdVWWrgHOL9adTOVW6IiImFa0De2J8DfBwMT+t+DuKyimpfyjmTyyeKxFxJpXWqk2Z+TuZOTsz5wAfKuK7asC2JlFpFdnTElMd10uoHAw3UOlKcHFEjIqIuVT+md2TmR/NzFnFPi6mcgr6XQeK9zCspHKaDQ6zbx8H+F5k5o4DfJ4GOyDvz4+B34G9fTlnA48dZpzVxgHPRETLnu3uT9Ha8i0q3Ts2wIG/y1RaLi8upg+4bR1xHsOreAz3GF6o9TH8gJ/RqPS7fiOVH0EfKl6ju4DXVH0e24tjffV6U6l0G/pX4GNUziAOKybWA2TmJiqnKR6KiE89z838DnBpRNwPLKVykQPAp6LSGf8hKqea7qdymuTkKC58odLf7X9GxH08/zMKXcX6/0Dl1/YLEVROaz5I5XTlDCoX60Dl1PqmiHiYyvP4cPH6ERFLBmznIgYclKn0ZR1L5ZRr9TBKF1Hpt3VJPHtBWXXr0cWDbOuvgFcXcd4KfCQzNwIvARYX78XtVFpE9hwoPhwRjwAPUOlPd1tR/lbgoeI5XAe8o2hFOZBri9fhJ8U+9lyg8qfAHxT7/xZwSVYsBW6k8g/iFuDyzDzYacH9xXuo/rKIczGHfwryAaA/Khdt1etikv8NlIr38AYqr1X3QdY5kP8O3E3lPXn0IHUvBI4HvrjnM1eU7++7fAVweRHrzBcQo2rMY/hzeAz3GA61P4bv9zNa/HD4IvB7Wbl49U+BrwAbqfzY/FZEPEClG8iLB2x3JnBH8d79H+CjNYj1iIqDf9Y0nETlNNmHMnNxo2NRY0TEyqK1RNIw4zFcQ/0Y7mf0wGyxliRJkmrAixdHmMw8u9ExqOH+rtEBSHp+PIaLIX4M9zN6YHYFkSRJkmpg2LVYn3feeXnLLbc0OgxJOtLi4FWqKj93SLZ5VG7w8Y2ifA6VUQYuKoa4Cio3rTgf2EXlIq1fHGgfHo8lHcUGPSYPuz7WGzdubHQIkjTkVQ/JRmWIrl1UhkC7Crg1M+dTGX3hqmKVN1IZMmw+cBnwhYPtw+OxJO1r2CXWkqTDdi6VG2o8SWXouK8X5V8H3lxMX0hlXN/Myg09JkblznmSpENkYi1JI1/1uMHTM/OZYnotlVsfQ2X82Keq1lnNIGN0R8RlEbE4IhZv2LChXvFK0rBkYi1JI1hxF7wLgH8ZuKy4acZhXcGemddn5oLMXNDR0VGjKCVpZDCxlqSR7Y3ALzJzXTG/bk8Xj+LvnltRr6Fyi+49ZvHsrbklSYfAxFqSRrZ3su/toxcB7ymm3wN8r6r83VHxSmBbVZcRSdIhGHbD7UmSDk1EtAOvB/6wqvhvgBsj4lLgSeCiovxmKkPtLacygsh7j2CokjQimFhL0giVmZ3AlAFlm6iMEjKwbgKXH6HQJGlEsiuIJEmSVAMm1pIkSVINmFhLkiRJNWAfa0mqs56efh54ehtrt3cxY3wbpx07gdbWpkaHJUkj2tbdXTy+tpN127uZPn4UJx3TzsTRbXXdpy3WklRHO3d38d37n+Zv/u1hygnLN+xk8arN7N7d2+jQJGnE2rq7i2e27CpugVW5D9YzW3axdXdXXfdb18Q6Is6LiMciYnlEXLWfOhdFxMMRsTQivlnPeCTpSNq9u5eHn+nkxsVP8ufnv5hjxo1iTGsTzU0lHt+4jb6+cqNDlKQRaWtnFw+u2cm7v3oPH/jWEt79lXt4cM1OtnYO08Q6IpqA66jc9etk4J0RcfKAOvOBjwKvycxTgD+pVzySdKQ9uHY763Z08z9/61R6+oO127uYOnYUqzbt5MlNXazeur3RIUrSiLR+ez9XL1pKV2+lAaOrt8zVi5ayfnt/XfdbzxbrhcDyzFyRmT3At4ELB9T5A+C6zNwCkJnrkaQRYt32buZOGcPmzh4CiIBSBLMmjwaSzTttsZakeli3o3tvUr1HV2+ZdTu667rfel68OBN4qmp+NfCKAXVOAoiInwBNwMcz85aBG4qIy4DLAGbPnl2XYCWp1qaPH0VLM6zb2MPyDZ2VPtbrd3JCRzsnTBvLxp32s5akepg+fhRtLaV9kuu2lhLTx4+q634bffFiMzAfOBt4J/DFiJg4sFJmXp+ZCzJzQUdHxxEOUZKen9OOGU9nV5mtu/ZNoLfu6qW7J+nutcVakurhpGPaueaCU2lrqaS6bS0lrrngVE46pr2u+61ni/Ua4Liq+VlFWbXVwN2Z2Qs8ERGPU0m0761jXJJ0RIwe3UJXf6U/39jWJmZMHMPunj4mt7fSR5k5U0Y3OMKhaffu3kr/9GKIrNOOGc/o0S2NDkvSMDJxdBtvOLWDOVMXHtHh9uqZWN8LzI+IuVQS6ouB3x5Q57tUWqq/GhFTqXQNWVHHmCTpiCpFANCf8OGb7qert0xbS4n/8ZbTOPOlkxsc3dCzc3cXDxfjzk4rRlH55eYdnDB5nMm1pMMycXQbC+fWN5EeqG5dQTKzD3g/8EPgEeDGzFwaEddExAVFtR8CmyLiYeB24MOZuanWsZTLyYoNO/nZLzeyYsNOyuWs9S4kaVD95WRjZw/X3rpsn6vT//w7D7Jqy64GRze07N7dy5ptXXuGnCUCunt76ezuY/2uzsYGJ0mHoK53XszMm4GbB5RdXTWdwJXFoy7K5eTWR9fxxIadzJg4hqc27+LprbtZePxk73wmqe66e8uUk0GvTl+/o4t5HWMbFNnQ0tdX5uH1W/nl+t17h8iq9Ik8hZfNHs/WzuT4KY2OUpIOrNEXL9bdExt38szW3Uwc00IAmbB1dy93Lt/gzRkk1d3xU9ppCvZeQLNHW0uJaeOO7CnKoWzpM9vo7y8NOu7sls5+OnvrO/asJNXCiE+sn9nWRSlg4bwJlf56o5qZNm4ULzm2jZWbvTmDpPqaO7Wd02ZN4Ipz5+9zdfpnLjqDOVPqe3X6cPLMti7W7eja77izWzp7GhSZJB26unYFGQp2dPVx9ksm8dNlO7h60UNVpxdP5RXzxtHXV6a5ecT/vpDUIKVScM6LpnNix1jOnD2JXT19zJ7cztyp7ZRK0ejwhowZE0bT3de/33FnW5p8rSQNfSM+o5w6tpWnt/TvTaphz+nFh1i3vZ+lz2xrcISSRrpSKZgzdSyvmDeFX3vxdE6YNtakeoBTZoynqVTmmgtOec64s1Pamzjt2Ofc4kCShpwR32I9fnQTj6/bNfjpxe3dtDYFpx+3n5UlSUdEc3OJ02dOYfKY7XzjvQtZt6My7uzMSU1MHzveM4uShoURn1jPnTyObbv3c3px3ChGtXiwlqShoLm5xNyOicz1BruShqkRn1W2tjYxcUzTIKcXT2HCmCZOmTGhwRFKkiRpJBjxLdYAM8ePZsOkHr56ycvZuLObqWNHUc4yx05o8/SiJEmSauKoyCrbR4/ijJkTaCoFEDSVgjNmTmBcne8XL0mSpKPHUdFiDZXkeuHcUY0OQ5IkSSPUUdFiLUmSJNWbibUkSZJUAybWkiRJUg2YWEvSCBUREyPipoh4NCIeiYhXRcTHI2JNRCwpHudX1f9oRCyPiMci4jcaGbskDUdHzcWLknQUuha4JTPfFhGtwBjgN4DPZub/qq4YEScDFwOnAMcC/xERJ2Vm/5EOWpKGK1usJWkEiogJwFnAlwEysycztx5glQuBb2dmd2Y+ASwHFtY/UkkaOUysJWlkmgtsAL4aEfdFxJcior1Y9v6IeCAivhIRk4qymcBTVeuvLsr2ERGXRcTiiFi8YcOGuj4BSRpuTKwlaWRqBs4EvpCZLwM6gauALwAnAGcAzwCfPpyNZub1mbkgMxd0dHTUOGRJGt5MrCVpZFoNrM7Mu4v5m4AzM3NdZvZnZhn4Is9291gDHFe1/qyiTJJ0iEysJWkEysy1wFMR8aKi6Fzg4YiYUVXtLcBDxfQi4OKIGBURc4H5wD1HLGBJGgEcFUSSRq4PAP9cjAiyAngv8LmIOANIYCXwhwCZuTQibgQeBvqAyx0RRJIOj4m1JI1QmbkEWDCg+HcPUP8TwCfqGpQkjWB2BZEkSZJqwMRakiRJqgETa0mSJKkGTKwlSZKkGjCxliRJkmrAxFqSJEmqARNrSZIkqQZMrCVJkqQaMLGWJEmSasDEWpIkSaoBE2tJkiSpBkysJUmSpBowsZYkSZJqwMRakiRJqgETa0mSJKkGTKwlSZKkGjCxliRJkmrAxFqSJEmqARNrSZIkqQZMrCVJkqQaqGtiHRHnRcRjEbE8Iq4aZPklEbEhIpYUj9+vZzySJElSvTTXa8MR0QRcB7weWA3cGxGLMvPhAVVvyMz31ysOSZIk6UioZ4v1QmB5Zq7IzB7g28CFddyfJEmS1DD1TKxnAk9Vza8uygZ6a0Q8EBE3RcRxg20oIi6LiMURsXjDhg31iFWSJEl6QRp98eL3gTmZ+VLgR8DXB6uUmddn5oLMXNDR0XFEA5QkSZIORT0T6zVAdQv0rKJsr8zclJndxeyXgF+pYzySJElS3dQzsb4XmB8RcyOiFbgYWFRdISJmVM1eADxSx3gk6agSEROLbnaPRsQjEfGqiJgcET+KiGXF30lF3YiIzxWjOD0QEWc2On5JGm7qllhnZh/wfuCHVBLmGzNzaURcExEXFNU+GBFLI+J+4IPAJfWKR5KOQtcCt2Tmi4HTqRyLrwJuzcz5wK3FPMAbgfnF4zLgC0c+XEka3iIzGx3DYVmwYEEuXry40WFI0pEWh1U5YgKwBJiXVQf6iHgMODsznynOGt6RmS+KiH8spr81sN7+9uHxWNJRbNBjcqMvXpQk1cdcYAPw1Yi4LyK+FBHtwPSqZHktML2YPqSRnBylSZL2z8RakkamZuBM4AuZ+TKgk2e7fQBQtGQf1mlLR2mSpP0zsZakkWk1sDoz7y7mb6KSaK/bc+F48Xd9sfygIzlJkg7MxFqSRqDMXAs8FREvKorOBR6mMjrTe4qy9wDfK6YXAe8uRgd5JbDtQP2rJUnP1dzoACRJdfMB4J+LIU9XAO+l0qByY0RcCjwJXFTUvRk4H1gO7CrqSpIOg4m1JI1QmbkEWDDIonMHqZvA5XUPSpJGMLuCSJIkSTVgYi1JkiTVgIm1JEmSVAMm1pIkSVINmFhLkiRJNWBiLUmSJNWAibUkSZJUAybWkiRJUg2YWEuSJEk1YGItSZIk1YCJtSRJklQDJtaSJElSDZhYS5IkSTVgYi1JkiTVgIm1JEmSVAMm1pIkSVINmFhLkiRJNWBiLUmSJNWAibUkSZJUAybWkiRJUg2YWEuSJEk1YGItSZIk1YCJtSRJklQDJtaSJElSDZhYS5IkSTVgYi1JkiTVgIm1JEmSVAMm1pI0QkXEyoh4MCKWRMTiouzjEbGmKFsSEedX1f9oRCyPiMci4jcaF7kkDU/NjQ5AklRXv5aZGweUfTYz/1d1QUScDFwMnAIcC/xHRJyUmf1HKE5JGvZssZYkAVwIfDszuzPzCWA5sLDBMUnSsGJiLUkjVwL/HhE/j4jLqsrfHxEPRMRXImJSUTYTeKqqzuqibB8RcVlELI6IxRs2bKhf5JI0DJlYS9LI9drMPBN4I3B5RJwFfAE4ATgDeAb49OFsMDOvz8wFmbmgo6Oj5gFL0nBmYi1JI1Rmrin+rge+AyzMzHWZ2Z+ZZeCLPNvdYw1wXNXqs4oySdIhMrGWpBEoItojYtyeaeANwEMRMaOq2luAh4rpRcDFETEqIuYC84F7jmTMkjTcOSqIJI1M04HvRARUjvXfzMxbIuKfIuIMKv2vVwJ/CJCZSyPiRuBhoA+43BFBJOnwmFhL0giUmSuA0wcp/90DrPMJ4BP1jEuSRrK6dgWJiPOKGw0sj4irDlDvrRGREbGgnvFIkiRJ9VK3FuuIaAKuA15PZdimeyNiUWY+PKDeOOAK4O56xSJJkoa+cjlZuamTddu7mD6+jTlT2imVotFhSYesni3WC4HlmbkiM3uAb1O5AcFAfwV8EuiqYyySJGkIK5eTW5au5c9uWsKGnT3cu3Izdz+xiZ4eu/pr+KhnYn3Qmw1ExJnAcZn5f+sYhyRJGuJWburkK//5S9565mw+/e+Psmrzbn62YhM/WbHR5FrDRsOG24uIEvAZ4E8Poa53+pIkaQRbt72Ld796Hv9w53LesWA2P3hgDf1luO+prfz0iU309ZUbHaJ0UPVMrA92s4FxwKnAHRGxEnglsGiwCxi905ckSSPb9PFt7O7p400vnckNi1ftk1z/YtUWfvLLjSbXGvLqmVjfC8yPiLkR0QpcTOUGBABk5rbMnJqZczJzDnAXcEFmLq5jTJIkaQiaM6Wd2ZPH0FRib3L9e6+eS1MJyllJrv/j0XWUy9noUKX9qltinZl9wPuBHwKPADcWNyC4JiIuqNd+JUnS8FMqBb9y3CRedtxEmkpw8ctnEwMGBFm/vYtVmzsbE6B0COp6g5jMvBm4eUDZ1fupe3Y9Y5EkSUNba2sTr5k3lYigpQkyg/Gje5jc3sKqTZ2UEzbt7GHO1EZHKg2uYRcvSpIkDdTa2sQrZk+iuy+5Z+Vmlm/YyX//3kM0NTXR1lKiL+1nraHLxFqSJA0Z5XLykyc2sau7n5kTRvPiY8Zx+etO4B/+33LGtLbQ128faw1dde0KIkmSdDie2tIJAeNHN9NXTsa3NdPWUuLS18yls6eP9lFNjQ5R2i9brCVJ0pCxq7uPEsHiJ7fs7QayaWcPp80aT3trM7Mntzc6RGm/TKwlSdKQUC4nKzbuYsnqrZQTmgJ+79Vz+d93LKenD8aNbmLuVBNrDV12BZEkSUPCqs2drNq8i+vvXEFXb5m2lhJXnDufi18+mw07ujn/tBmUSnHwDUkNYmItSZKGhHXbu7n21mVMGtPKb505iwjo6u3nlGMnMLm91aRaQ56JtSRJGhI6e/qYNKaV333l8XzutmV7W63/+s2nctYJDl6toc8+1pIkaUg4fnI7b18wa29SDdDVW+Zj332Ip3d0NTg66eBMrCVJ0pAwd2o7J00btzep3qOrt8x6E2sNAybWkiRpSCiVgpfMGE9by77pSVtLiWnj2hoUlXTo7GMtSZKGjLlT2/nMRWfwyVse4U0vnUlTCV5+/GRmTxrT6NCkgzKxliRJQ0apFLzhJdPp7S/zkX99YO8FjJ+56AzOO+UYRwbRkGZXEEmSNKSs2rJrb1INlT7WV964hJWbOhscmXRgJtaSJGlIWbe9ywsYNSzZFUSSJA0p08e3cfyU0bzppTOJoufH9+9f4wWMGvJMrCVJ0pAye9IYPnDOfD723Yf2uUmMFzBqqDOxlqQRKiJWAjuAfqAvMxdExGTgBmAOsBK4KDO3REQA1wLnA7uASzLzF42IW1q1ZRd/f9syLn3tvL0t1n9/2zLOnD2JeR1jGxucdAD2sZakke3XMvOMzFxQzF8F3JqZ84Fbi3mANwLzi8dlwBeOeKQS0NdXZuPObq58/Yt4yYxx3PnYer704xW8Y8FsNnd2Nzo86YBssZako8uFwNnF9NeBO4CPFOXfyMwE7oqIiRExIzOfaUiUOir19ZX5t6XPsGz9TsoJTQG/f9Y8oMyn/30ZrzlxSqNDlA7IxFqSRq4E/j0iEvjHzLwemF6VLK8FphfTM4GnqtZdXZSZWOuIeXTddlZv2c31d67Y27f6inPnc+bxE7n8dSfQ018++EakBjKxlqSR67WZuSYipgE/iohHqxdmZhZJ9yGLiMuodBVh9uzZtYtUArbs6uXaW5cxaUwrv3XmLCKgq7efLCfHTBzN1l29rNiwkweM+OoAACAASURBVDlT2r1RjIYk+1hL0giVmWuKv+uB7wALgXURMQOg+Lu+qL4GOK5q9VlF2cBtXp+ZCzJzQUdHRz3D11Goq7efSWNa+d1XHs8PHljDtLGtnP2iqTy9rZvL/unn/ME3fs75n/sxtyxdS7l8WL8JpSPCxFqSRqCIaI+IcXumgTcADwGLgPcU1d4DfK+YXgS8OypeCWyzf7WOtHlT2nn7glncsHgVl7/uBKaOG8WWzj7+/DsPPucujE9s9C6MGnrsCiJJI9N04DuVUfRoBr6ZmbdExL3AjRFxKfAkcFFR/2YqQ+0tpzLc3nuPfMg62s3tGMv8aeN400tnMmZUC+UyLFm9ddC7MD61uZMTpjn0noYWE2tJGoEycwVw+iDlm4BzBylP4PIjEJq0X6VScPKM8Sxbv4PO7j4AygltLaV9kuu2lhJjWk1hNPTYFUSSJA0Zc6e28/LjJ9Pe1kx7WzPfv38NHzxnPm0tlZSlraXEf3/TyWzv7m1wpNJzmVhLkqQho1QKXjVvCmNaS5QCPnjuSdyweBWXvnYeHzz3RD5z0Rl85xdP0dE+qtGhSs9hYi1JkoaU5uYSZ580nVOOHc/Jx4zlL/7LKTSVoL8Mn7zlES5eeDynHDuh0WFKz2EHJUmSNOSUSsHxUyoXJ86fVmZKeytrt3XxhpOnccqMCTQ32zaoocfEWpIkDWnNzSVOP24Spx938LpSIx30515EfCAiJh2JYCRJkqTh6lDOo0wH7o2IGyPivCgGRZUkSZL0rIMm1pn5MWA+8GXgEmBZRPyPiDihzrFJkiRJw8Yh9fwvbhywtnj0AZOAmyLib+sYmyRJkjRsHPTixYi4Ang3sBH4EvDhzOyNiBKwDPiz+oYoHVnlcrJ6aydrt3azfkc3Mya0cdqxE2htbWp0aJIkaQg7lFFBJgO/lZlPVhdmZjki3lSfsKTG6Onp52dPbGTr7l5+uaGTcsKy9TtYs3U35518jMm1JEnar0PpY/0XA5PqqmWP1D4kqTF6evr5zxUbGdXcxLZd+94qd+uuHh5eu71BkUmSpOHAcawloK+vzOJVW1jy1FZeOXcynT39XH/nCrp6y7S1lLji3Pns7u1rdJiSJGkI87ZFOuqVy8lPV2ziqS27KCf0J1x76zK6essAdPWWufbWZfSXGxyoJEka0kysddRbsWEni5/czJjWZr5//xo6u/v2JtV7dPWW6ertb1CEkiRpODCx1lHviU2VixS//tMVvO+sE9m0s5u2ln2/Gm0tJeZMaW9QhJIkaTgwsdZRr62lie/fv4ZzXnwM//qLVUwd18Zfv/nUvcl1W0uJT7/9DOZ1jG1wpJIkaSir68WLEXEecC3QBHwpM/9mwPL3AZcD/cBO4LLMfLieMUkDjWoJLn75bL597yre9NKZLH16G6fPmsA3fm8h9zyxmYVzJrNgzmRKpWh0qJIkaQirW4t1RDQB1wFvBE4G3hkRJw+o9s3MPC0zzwD+FvhMveKR9qejvY321iYuPGMmUeTOa7bsZsX6HTSVgo5xo0yqJUnSQdWzK8hCYHlmrsjMHuDbwIXVFTKzemDgdiDrGI80qOOntHPspNH7lI0b3cL2rj7aW5v2JtuSJEkHUs/EeibwVNX86qJsHxFxeUT8kkqL9QcH21BEXBYRiyNi8YYNG+oSrI5epVJwzoum85oTptLaVOLEaeNYt72LbV39fOH/rWDt9q5GhyhJkoaBht8gJjOvA66LiN8GPga8Z5A61wPXAyxYsMBWbdVcqRRMHTuK/33H8n2G2mtrKTFtXFsDI5MkScNFPVus1wDHVc3PKsr259vAm+sYj7Rf5XLyxKadXHHu/H1GA/nMRWc4zJ4kSTok9WyxvheYHxFzqSTUFwO/XV0hIuZn5rJi9jeBZUgNsHJTJ+//5n1MGtPKpa+dRwSUAk6eMc4LFyVJ0iGpW2KdmX0R8X7gh1SG2/tKZi6NiGuAxZm5CHh/RPw60AtsYZBuINKRsG57F129ZZ7Z1sV1ty/fW/7qE6YwZ6rjV0uSpIOrax/rzLwZuHlA2dVV01fUc//SoZo+vo22lpL9qyVJ0vPmnRclYPakMc+52+Jfv/lUZk8a0+DIJEnScNHwUUGkoWDVll38/W3L9vavzoS/v20ZZ86e5K3MNawVN+taDKzJzDdFxNeA1wHbiiqXZOaSiAgqd8o9H9hVlP+iETFL0nBlYi1R6WP95Kbd+/SvBli/o8vEWsPdFcAjwPiqsg9n5k0D6r0RmF88XgF8ofgrSTpEdgXRUa9cTvr6c283kD3sY63hLiJmURlx6UuHUP1C4BtZcRcwMSJm1DVASRphTKx11HtiYycf+96DfPCcfcew/sRbTmPWhNEHWVsa0v4O+DOgPKD8ExHxQER8NiJGFWWHerdc74QrSfthVxAd9Z7c3MmTm3bzT3c9uU8f65YSPLx2O2fMntToEKXDFhFvAtZn5s8j4uyqRR8F1gKtVO5o+xHgmkPdrnfClaT9M7HWUa+9tZm2ltI+Y1i3tZT41NtOZ83W3SbWGq5eA1wQEecDbcD4iPg/mfmuYnl3RHwV+FAxf7h3y5UkDWBXEB31po8f9ZxbmV9x7nye2bqLye2tlMs2ymn4ycyPZuaszJxD5c63t2Xmu/b0my5GAXkz8FCxyiLg3VHxSmBbZj7TiNglabiyxVpHvdmT2zl+8hguO2se5azcyry9tYlxo1toaylx22PrOOdF0721uUaKf46IDiCAJcD7ivKbqQy1t5zKcHvvbUx4kjR8mVjrqFcqBcdObGNsWzPbdvfR3tpEXznpK5fp7u3ngdXbmDtlLCdMc9g9DU+ZeQdwRzF9zn7qJHD5kYtKkkYeu4JIQHNTibXbdhPAhp3d9PSV2dnVy7aufsoJT27qbHSIkiRpiDOxloAXTx9PW2szj6/fwVNbdvP4+h20NDXxxTuXUwpobfarIkmSDsyuIBLQ3Fzi9SdNY0r7KNZs2UVbazPf+OkKzn3JMbS3NtHcZP9qSZJ0YDbD6ahXLierNu/krlVb6O3vZ9KYVrp6+njPq+exJ58eO8rfoJIk6cDMFnRUK5eTHy9fT18/3L96K3c8up73vmYuq7fu3jtCyDETR/Pi6eMbHaokSRribLHWUW3lpk527O7n/tVbKSc8vn4nX/3JE5w4bRzHTRrNSdPGcfqsCTTbx1qSJB2E2YKOauu2d9HZ3Uc54fv3r+GD58zn8fU7+eC37uMvFi2lP5PjJrU3OkxJkjQM2BVER7Xp49vY1NnD9+9fwzsWzOaGxau49LXzaCrBmbMn8doTpnpjGEmSdEhssdZRbc6Udsa1NfG+153IDYtX8aaXzqSpBC87biKvnjvFLiCSJOmQ2WKto1qpFPzqidN4aksnJ0w9nY2d3cycMJpTjrVftSRJOjwm1jrqlUrB8VPGcvwUb1kuSZKeP5vkJEmSpBowsZYkSZJqwMRakiRJqgETa0mSJKkGTKwlSZKkGjCxliRJkmrAxFqSJEmqARNrSZIkqQZMrCVJkqQaMLGWJEmSasDEWpIkSaoBE2tJkiSpBkysJUmSpBowsZYkSZJqwMRakiRJqgETa0mSJKkGTKwlaQSLiKaIuC8iflDMz42IuyNieUTcEBGtRfmoYn55sXxOI+OWpOHIxFqSRrYrgEeq5j8JfDYzTwS2AJcW5ZcCW4ryzxb1JEmHwcRakkaoiJgF/CbwpWI+gHOAm4oqXwfeXExfWMxTLD+3qC9JOkQm1pI0cv0d8GdAuZifAmzNzL5ifjUws5ieCTwFUCzfVtTfR0RcFhGLI2Lxhg0b6hm7JA07JtaSNAJFxJuA9Zn581puNzOvz8wFmbmgo6OjlpuWpGGvrol1RJwXEY8VF8NcNcjyKyPi4Yh4ICJujYjj6xmPJB1FXgNcEBErgW9T6QJyLTAxIpqLOrOANcX0GuA4gGL5BGDTkQxYkoa7uiXWEdEEXAe8ETgZeGdEnDyg2n3Agsx8KZU+fX9br3gk6WiSmR/NzFmZOQe4GLgtM38HuB14W1HtPcD3iulFxTzF8tsyM49gyJI07NWzxXohsDwzV2RmD5UWkwurK2Tm7Zm5q5i9i0rriSSpfj4CXBkRy6n0of5yUf5lYEpRfiXwnLOMkqQDaz54ledt74UwhdXAKw5Q/1Lg3wZbEBGXAZcBzJ49u1bxSdJRITPvAO4opldQafgYWKcLePsRDUySRpghcfFiRLwLWAB8arDlXiwjSZKkoa6eLdZ7L4QpVF8ks1dE/Drw34DXZWZ3HeORJEmS6qaeLdb3AvOL2+e2Url4ZlF1hYh4GfCPwAWZub6OsUiSJEl1VbfEurjBwPuBH1K5ne6Nmbk0Iq6JiAuKap8CxgL/EhFLImLRfjYnSZIkDWn17ApCZt4M3Dyg7Oqq6V+v5/4lSZKkI2VIXLwoSZIkDXcm1pIkSVINmFhLkiRJNWBiLUmSJNWAibUkSZJUAybWkiRJUg2YWEuSJEk1YGItSZIk1YCJtSRJklQDJtaSJElSDZhYS5IkSTVgYi1JkiTVgIm1JEmSVAMm1pIkSVINmFhLkiRJNWBiLUmSJNWAibUkSZJUAybWkiRJUg2YWEuSJEk1YGItSZIk1YCJtSRJklQDJtaSNAJFRFtE3BMR90fE0oj4y6L8axHxREQsKR5nFOUREZ+LiOUR8UBEnNnYZyBJw09zowOQJNVFN3BOZu6MiBbgPyPi34plH87MmwbUfyMwv3i8AvhC8VeSdIhssZakESgrdhazLcUjD7DKhcA3ivXuAiZGxIx6xylJI4mJtSSNUBHRFBFLgPXAjzLz7mLRJ4ruHp+NiFFF2UzgqarVVxdlA7d5WUQsjojFGzZsqGv8kjTcmFhL0giVmf2ZeQYwC1gYEacCHwVeDLwcmAx85DC3eX1mLsjMBR0dHTWPWZKGMxNrSRrhMnMrcDtwXmY+U3T36Aa+Ciwsqq0BjqtabVZRJkk6RCbWkjQCRURHREwspkcDrwce3dNvOiICeDPwULHKIuDdxeggrwS2ZeYzDQhdkoYtRwWRpJFpBvD1iGii0ohyY2b+ICJui4gOIIAlwPuK+jcD5wPLgV3AexsQsyQNaybWkjQCZeYDwMsGKT9nP/UTuLzecUnSSGZXEEmSJKkGTKwlSZKkGjCxliRJkmrAxFqSJEmqARNrSZIkqQZMrCVJkqQaMLGWJEmSasDEWpIkSaoBE2tJkiSpBkysJUmSpBowsZYkSZJqwMRakiRJqoG6JtYRcV5EPBYRyyPiqkGWnxURv4iIvoh4Wz1jkSRJkuqpbol1RDQB1wFvBE4G3hkRJw+otgq4BPhmveKQJEmSjoTmOm57IbA8M1cARMS3gQuBh/dUyMyVxbJyHeOQJEmS6q6eXUFmAk9Vza8uyg5bRFwWEYsjYvGGDRtqEpwkSZJUS8Pi4sXMvD4zF2Tmgo6OjkaHI0mSJD1HPRPrNcBxVfOzijJJkiRpxKlnYn0vMD8i5kZEK3AxsKiO+5MkSZIapm6JdWb2Ae8Hfgg8AtyYmUsj4pqIuAAgIl4eEauBtwP/GBFL6xWPJEmSVE/1HBWEzLwZuHlA2dVV0/dS6SIiSZIkDWvD4uJFSZIkaagzsZYkSZJqwMRakiRJqgETa0mSJKkG6nrxoiRJA23d3cWytZ2s29HNlPZWxo5q4sQpYxk9uqXRoUnSC2JiLUk6Yrbu7uLOxzexfP1OyglNAXOmtvP01i7OOmGqybWkYc2uIJI0AkVEW0TcExH3R8TSiPjLonxuRNwdEcsj4obiBl5ExKhifnmxfE6tYyqXkxXrd7Fmy26uv3MFn79tOf945wrWbusiInhw7fZa71KSjigTa2D37l7ueWIT37//ae55YhO7d/c2OiRJeqG6gXMy83TgDOC8iHgl8Engs5l5IrAFuLSofymwpSj/bFGvplZt7mRHdx/X3rqMrt4yAF29Za69dRnNpWDd9u5a71KSjqijPrHevruLO1dsZP32bnZ197F+Rzd3rtjI1t1djQ5Nkp63rNhZzLYUjwTOAW4qyr8OvLmYvrCYp1h+bkRELWPa3dNHV2//3qR6j67eMp09/UwfN6qWu5OkI+6o7mO9e3cvy9Z3srmzl7/8/lK6esu0tZT4i/9yCk9t2k3LlKB9tAd6ScNTRDQBPwdOBK4Dfglszcy+ospqYGYxPRN4CiAz+yJiGzAF2Dhgm5cBlwHMnj37sOLZ0dXPzImjaGsp7ZNct7WUmDimhWMne7yVVH/lcrJyUyfrtncxfXwbc6a0UyrVph3hqG6xfnDtdrp6y3uTamDv/M7ufpau3XmQLUjS0JWZ/Zl5BjALWAi8uAbbvD4zF2Tmgo6OjsNad92Obrp6y3zyrS+lraXy76etpcT/fMtpjG8rcez49hcaniQdULmc3LJ0Led/7se884t3c/7nfswtS9dSLmdNtn9Ut1jv6c832GnJTZ09ZG1eY0lqqMzcGhG3A68CJkZEc9FqPQtYU1RbAxwHrI6IZmACsKmWcRwzfhSbd/USlPni7y5g6+5eJre3ECRN0VSzFiNJ2p8nNnZy5Y1L9mlQvfLGJbzoA7/KCdPGvuDtH9Ut1tPHj6KtubS35WSPtpYS08eNsr+fpGErIjoiYmIxPRp4PfAIcDvwtqLae4DvFdOLinmK5bdl1rZ54bQZE8hy0p/BvU9u5rF1O7j7ic1s7OzlpGnjarkrSRrUk5s6B21QXbW5sybbP6oT69OOGU85k4//l1P2OS35VxeeyugWmDK2qcERStLzNgO4PSIeAO4FfpSZPwA+AlwZEcup9KH+clH/y8CUovxK4KpaB9TW1szrTuzghKntvHLeFF5yzDhed1IH573kGFpbPd5Kqr9RLYM3qLY21yYlPqq7gowe3cJZJ0zl8Y07+colL2dLZw/Tx49iSnsT23f30d56VL88koaxzHwAeNkg5Suo9LceWN4FvL3ecbW1NfPS4ybVezeSNKjmUnDFufP3DvvZ1lLiinPn01yjrmhHfeY4enQLp82cyKot22iO4JltXZTLbUyb0Mw0L6SRJEkaEXp6+iln0t7axGVnzaOcUApob22iv0Y93476xBqgVApmT5pAudxJmWTquFHMnlS7oVckHR3K5WT99p2s2tzDuh3dTB8/ipOOaWfi6LZGhyZJR7Wenn6+e//TnDht8EbTtqbadEczsS6USsG8jrHM63jhV4RKOvpUkupt3LlsO1cvenZc/GsuOJU3nNphci1JDfTA09t4fO0WTj9uHBPHtLCxs2fvsoljWmhvs4+1JA0ZKzd1snFH/96kGipXml+96CHmTF3Iwrkm1pLUKGu3d/Gm04+jtw+mjW9l+vjRrN/RxbRxbST99PQdfBuH4qgeFUSSamXd9i7W7egadBinPWPmS5IaY8b4Nnb19kPArp4y3X39lBO6+/rZ1VOmVp1/TawlqQamj29j+vi2/Y6LL0lqnNOOncDmzh6am3hOEh3AqJba7MfEWpJqYM6UdmZPauKaC/YdF/+aC05l9mTHaJakRmptbWLWxNEsfmIL7W1NtDU3UQpoa26iva2Jlqba9I62j7Uk1UCpFExsG8tZ8+Eb711YGRVk3ChmT25i8hjvKihJjXbqsRNYtWUXKzfsYkxrC7t6+lm/o5txbU3Mnl2bIZZNrCWpRtramplcGsfqrdsoRSXZnjxmnHcVlKQhoLm5xPmnzODRddvZsquXUS0l5k5pZ17H2JoNsWxiLUk11NraxII5kxsdhiRpEM3NJU6dObFu27ePtSRJklQDJtaSJElSDZhYS5IkSTVgYi1JkiTVgIm1JEmSVAMm1pIkSVINmFhLkiRJNWBiLUmSJNWAibUkSZJUA5GZjY7hsETEBuDJRsdxEFOBjY0O4gXyOQwNPoeho9HPY2NmntfA/T9HjY/HjX59D2SoxjZU44KhG5txHb6hGluj4xr0mDzsEuvhICIWZ+aCRsfxQvgchgafw9AxUp7HUDWUX9+hGttQjQuGbmzGdfiGamxDNS67gkiSJEk1YGItSZIk1YCJdX1c3+gAasDnMDT4HIaOkfI8hqqh/PoO1diGalwwdGMzrsM3VGMbknHZx1qSJEmqAVusJUmSpBowsZYkSZJqwMT6eYqI8yLisYhYHhFXDbL8yoh4OCIeiIhbI+L4RsR5MIfwPN4XEQ9GxJKI+M+IOLkRcR7IwZ5DVb23RkRGxJAbnucQ3odLImJD8T4siYjfb0ScB3Io70NEXFR8L5ZGxDePdIwHcwjvw2er3oPHI2JrI+IcjiLiuIi4ver9v6IonxwRP4qIZcXfSUV5RMTnivfigYg4s05xtUXEPRFxfxHXXxblcyPi7mL/N0REa1E+qphfXiyfU4+4quJrioj7IuIHQyyulVX/GxYXZQ19L4t9TYyImyLi0Yh4JCJeNUTielHVsWNJRGyPiD8ZIrH91+Kz/1BEfKv4TjT8cxYRVxQxLY2IPynKGv56HVRm+jjMB9AE/BKYB7QC9wMnD6jza8CYYvqPgBsaHffzfB7jq6YvAG5pdNyH+xyKeuOAO4G7gAWNjvt5vA+XAJ9vdKwv8DnMB+4DJhXz0xod9/P5LFXV/wDwlUbHPVwewAzgzGJ6HPA4cDLwt8BVRflVwCeL6fOBfwMCeCVwd53iCmBsMd0C3F3s70bg4qL8H4A/Kqb/GPiHYvrieh/bgSuBbwI/KOaHSlwrgakDyhr6Xhb7+jrw+8V0KzBxKMQ1IMYmYC1wfKNjA2YCTwCjqz5flzT6cwacCjwEjAGagf8ATmz063UoD1usn5+FwPLMXJGZPcC3gQurK2Tm7Zm5q5i9C5h1hGM8FIfyPLZXzbYDQ+1q14M+h8JfAZ8Euo5kcIfoUJ/DUHYoz+EPgOsycwtAZq4/wjEezOG+D+8EvnVEIhsBMvOZzPxFMb0DeITKP/ULqSRDFH/fXExfCHwjK+4CJkbEjDrElZm5s5j9/9m78zg7rvLA+7+n19urWkurJUuyJRsBtmVsY0U4wxIPDoyBIEEgiZlJsBkSv5lAQoZkgplkBDaZhCxvGPLGk8QBgsmQGI9ZLAeHJcYeJwzGlsGLZGEsvGDJUqutpbvV6r3P+0dVy1ft1mbd27f76vf9fO6nq07VrXrudvqpU6dO1eePBLweuPUocU3GeytweUREqeMCiIjlwFuAT+XzMRviOoaKfpYRMQ94HfBpgJTSSErpQKXjmsblwI9SSk/PktjqgKaIqCNLZHdR+e/ZuWTJ8aGU0hjwf4CfZXa8X8dkYv3iLAOeKZrfkZcdzXvJjqRmmxN6HRHxvoj4EdmR4m/MUGwn6rivIT8ltCKl9NWZDOwknOj36R35Ka5bI2LFzIR2wk7kNbwUeGlEfDsi7o2IWXV7bk7idx1Z165VwLdmIK6qk58+vpisdbgrpbQrX7Qb6MqnT7aePZV4aiPiQWAP8E2yMxcH8n/oU/d9OK58eS+wsBxxAf8D+B1gIp9fOEviguzg4xsR8UBEXJOXVfqzXAX0AH8bWfeZT0VEyyyIa6oref6gvKKxpZR2An8K/Jgsoe4FHqDy37MtwGsjYmFENJO1SK9g9n2WL2BiXWYR8YvAWuBPKh3Li5VSuiGldA7wIeD3Kh3PyYiIGuDPgN+qdCyn6HZgZUrpFWT/+G86zvqzUR1Zd5DLyFp7/yYiOioa0Yt3JXBrSmm80oHMNRHRCnwR+M0pZ8RI2TndGT8rllIaTyldRHZmcR3w8pmOYaqI+BlgT0rpgUrHchSvSSm9EngT8L6IeF3xwgp9lnXAK4G/TCldDAyQdReodFyH5X2V1wP/e+qySsSW91HeQHZQcgbZmemKN3qklLaRnWX+BvA14EFgfMo6Ff0sj8bE+sXZSXbkNGl5XnaEiPhp4HeB9Sml4RmK7WSc0OsocjPPn3aZLY73GtrI+mrdHRFPkfW92hSz6wLG434OKaW9Rd+hTwGXzFBsJ+pEvks7gE0ppdGU0pNkfWxXz1B8J+Jkfg/FLU46QRFRT5ZUfz6l9KW8uHvylG3+d7KL0MnWT6cs7zZwF/CTZKeS66bZ9+G48uXzgL1lCOfVwPq83rqZ7NT8J2dBXMDhls7JLl1fJjsgqfRnuQPYkVL6bj5/K1miXem4ir0J+F5KqTufr3RsPw08mVLqSSmNAl8i++5V/HuWUvp0SumSlNLrgP1k/zMq/X4dl4n1i3M/sDq/araB7J/spuIVIuJi4K/JkurZ1pd00om8juLE5y3A4zMY34k45mtIKfWmlBallFamlFaS9Xdfn1LaXJlwp3Uin0NxX7H1ZP1TZ5PjvgbgK2St1UTEIrKuIU/MZJDHcSKvgYh4OTAf+M4Mxzen5f0wPw1sSyn9WdGiTcBV+fRVwG1F5e/Or/a/FOgtOgVcyrg6J8+cREQT8Aay39ddwDuPEtdkvO8EvpW3nJVUSunDKaXleb11Zb6f/1DpuAAioiUi2iangTeSnbqv6GeZUtoNPBMRL8uLLgcerXRcU0y9NqPSsf0YuDQimvPf6OR7Nhu+Z4vzv2eS9a/+eyr/fh1fqtBVk3P9Qdbf54dkffF+Ny+7nixpg+wK1m6y0xcPkrXUVTzuF/E6PglszV/DXcD5lY75ZF/DlHXvZpaNCnKCn8Mf5p/DQ/nn8PJKx/wiXkOQdct5FHiE/Irz2fQ4ke8S8FHg45WOda49gNeQnbZ9uKhefDNZ/8w7yQ7a/xlYUPR9uSH/LB4p1+8WeAXZaDUPkyWHG/Pys4H7gO1kp+0b8/JCPr89X372DLx3l/H8qCAVjyuP4aH8sbXot1LRzzLf10XA5vzz/ArZQXDF48r310LWujuvqKzisQHXAT/Iv/9/BzTOku/Zv5D9v3gIuHy2vF/He3hLc0mSJKkE7AoiSZIklYCJtSRJklQCJtaSJElSCZhYS5IkSSVgYi1JkiSVgIm1JEmSVAIm1pIkSVIJmFhLkiRJJWBim++puwAAIABJREFULUmSJJWAibUkSZJUAibWkiRJUgmYWEuSJEklYGINRERHRPxa0fwZEXFrGfd3UUS8uVzbL4WI+HxEPBYRWyLiMxFRX7Tssoh4MCK2RsT/KSr/QL7+1oj4zaLyL+TrPxgRT0XEg3n5uqLyhyLi7UXPeSoiHsmXbT5KjPMj4ssR8XBE3BcRa/LylxVt98GI6JuMJyJ+Lo9vIiLWFm3rDRHxQL7PByLi9UXLLsnLt0fEn0dE5OULIuKbEfF4/nd+Xh75etvz2F5ZtK0/zve/bcq2vpa/B1sj4q8iovY4+9iQb/vBiNgcEa85ynt0d0SsPP4nfnwR8amIOC+f/q+l2Gal5d/lf8yn10fEtSXY5t3F3y2Vn3X4C1mHW4dPs60Zq8NP5TcYER+NiN8udUwzJqV02j+AlcCWGdzf1cBfVPp1HyfGNwORP/4B+E95eQfwKHBmPr84/7sG2AI0A3XAPwMvmWa7/y+wMZ9uBury6aXAnqL5p4BFx4nxT4CP5NMvB+6cZp1aYDdwVj5/LvAy4G5gbdF6FwNnFL2WnUXL7gMuzd+LfwLelJf/MXBtPn0t8EdF790/5etfCnw3L/83wLfzmGqB7wCX5cva878BfBG48jj7aAUin34F8IOjvEd3AyvL8P04OAu+owHUnOI2LgP+scRxHfHd8jEj3wXr8BfGaB3+/DLr8Bdut+J1+DFi+yjw25WO48U+bLHOfBw4Jz9y/JOIWBkRWwAi4uqI+Ep+pPlURLw/Ij4YEd+PiHsjYkG+3jn5EesDEfEvEfHyvPzn8haAhyLinohoAK4HfiHf3y/kR/3fybf5fyPiZSe577sj4pP59rZExLpTfUNSSnekHFmltDxf9O+BL6WUfpyvtycvP5es8jmUUhoD/g/ws8XbzI/sf56skqdoXYACkE4yzPOAb+Xb+gGwMiK6pqxzOfCjlNLT+XrbUkqPTfN6v59Sejaf3Qo0RURjRCwlqzDvzd+LzwFvy9fbANyUT980pfxz+dt3L9CRbyflr7MBaATqge58/335c+vy5aloWy/YR0rpYB4PQAsn+d5FxGcj4p1F8wfzv5fl36dbI+IHkbV6TbbI3B0RayPi4/n782BEfP5k9jtNHB+NrDXt7oh4IiJ+o2jZB/Pv85Z4vrVqZWStcJ8jSwJem8f52Yj4YR7vT0fEt/MWonX586b9jU2J5eqI+It8uri1bDAifioiWvJY78u3syFftykibo6sBevLQNOpvCd6UazDp7AOtw6fwTr87/Lv/+MR8St5efFv8D9HxGfy6Qvy73jz0X5zU7b/GxHxaGSt+zefSqwzptKZ/Wx4MKW1o3ierGViO9AGdAK9wK/myz4B/GY+fSewOp9+FfCtfPoRYFk+3VG0zb8o2l87zx/l/zTwxZPc993A3+TTr2OalhuyI/wHj/LoOMZ7Uw98D3htPv8/gBvyfT4AvDsvPxf4IbCQrBXjO8D/N2VbrwM2Tyl7FVkleBB4e1H5k/l+HwCuOUpsfwB8Ip9eB4wBl0xZ5zPA+6d57t0cpVUReCfwz/n02snpfP615C2cwIGi8picB/4ReE3Rsjsn9wX8KXAg/yz/+5T9fh3YD/w9UHusfeTzbwd+AOwDfvIor+VupmntAD4LvLNo/mD+97I8tuVkXcW+M/lait8zjtHaAXzhKN+zd0+z7keB/0v2T2oRsDf/zl1C9ttpIWvZ2UrWIrUSmAAuLfqtjgEX5PE+kH/mQfYP7SvH+Y1dVvR5Xs2UVkjgrcC/5DH9AfCLk79lsu97C/BB4DN5+SvyeGyxnsEH1uHW4Ucusw6f2Tr8IbIGhUXAM8AZHPkbrAHuyV/vZuDVx/nNfZS8xRp4Fmgs/v3N9kcdOhF3pZT6gf6I6AVuz8sfAV4REa1kp4j+d35gCFmiANlpo89GxC3Al46y/XnATRGxmuyotb5o2TH3XbTeZAvCPRHRHhEdKaUDkwtTdoR/0cm+cOB/AveklP4ln68jS3ouJ/shfSci7k0pbYuIPwK+AQyQ/QjHp2zrXZNxFsX1XeD8iDg3fw/+KaU0RFYR7IyIxcA3I+IHKaV7pmzv48AnI+vv9wjw/eJ95i1L64EPn+iLjYjzgT8C3niiz8lfR4qIY7Y4RMRLyP55TbYcfTMiXjv53qaU/l1EFIDPA68HvnmsfaSUvgx8OSJeB3yM7B96KdyXUtqRx/wgWQX5ryf65JTSL5zk/r6aUhoGhiNiD9AFvAb4ckppII/jS2T/EDcBT6esFWnSkymlR/L1tpKdTk4R8UgeOxz7NzatfN0/Af5tSmk0It4IrI/n+/4VgDPJko0/z1/7wxHx8Em+fpWfdbh1+DFZhx8R58nW4bellAaBwYi4i+wg6cGi7U1ExNXAw8Bfp5S+fZzfXLGHgc9HxFeAr5xkXBVhYn1ihoumJ4rmJ8jewxqyo9AXVHoppV+NiFcBbwEeiIhLptn+x8gq37dHdpHC3Sex78O7mrrr4pn81OQXptk3ZH3EDkwtjIiPkLWy/D9FxTuAvXnCMxAR9wAXAj9MKX0a+HT+3D/I153cVh3ZacXpXj95pX6QrG/c5pTSzrx8T2Sn19eRHfEWP6cPeE++/SBrIXmiaJU3Ad9LKXUf5XVPfb3LgS+THZX/KC/eyfOVKPn0zny6OyKWppR25acJ9xQ9Z8U0z/lF4N6U0uQpu38CfpKsRXTyNQ1FxG1kra3fPMY+it+HeyLi7IhYlFJ67kReK1nLUE0eRw3ZqctJxd+5cU6ynoiIL5C1rk31Zymlz01TfrL7GzjG84/2GznWb+wF8kr/FuBXUkq7JouBd6Qpp6GL/ilo9rIOf551uHX4Mb2IOvyY393carKzGmfk80f9zU3xFrLGi7cCvxsRF6Tnux/NSvaxzvSTnap7UfLK4cmI+Dk4fEXxhfn0OSml76aUNgI9ZD/Wqfubx/M/9KtfZBi/kO/vNUBvSql3SoyPpZQuOspjugr5l4F/B7wrpTRRtOg24DURURcRzWSnb7blz1mc/z2TrAL++6Ln/TTZxRnFFfWqvLImIs4iu3jlqcj6srbl5S1kLQ9bpomxI2/RAPhlslaZvqJVXtC6cjQR0QF8lewik29PludJVV9EXJpX/O/O3wPIWk+vyqevmlL+7vx7cCnZ57EL+DHwU/l7Vw/8FLAtIlrzCnfyn9dbyE4PHnUfEfGSPB4iu2K9kawbxYl6iuf/Qa7nBFpwpxiNolEGiqWUfuEo37PpKuSj+RfgbZH1w2shO4X4L8d5zrGc7G/sM8DfFrXyQXaa99eL3veL8/J7yPqtEtmoBsWtkJoZ1uFTWIdnrMOPqpR1+IaIKETEQrKuKPcXL4yIeWRn9V4HLIyIdx7rN1f0vBpgRUrpLuBDZL+z1pN8nTPOxBpIKe0Fvh1Zh/o/eZGb+Q/AeyPiIbL+Zhvy8j+JbJifLWR9SR8C7gLOi/zCF7Krhv8wIr7Piz+LMJQ//6+A977IbRT7K7JT8t/J49wIWasE8DWy0zP3AZ9KKU1WmF+MiEfJTnW+b0plfyUvrCBfAzwU2amqLwO/lh+tdwH/mr+X95F1FfgaQET8akT8av78c4EtEfEYWcvGByY3nFfmb2DKqduIeHtE7CBrZfhqRHw9X/R+4CXAxnj+orXF+bJfAz5F1lfyR2RXi0N2GvMNEfE42T+dj+fld5C1umwH/iZ/PsCt+fMfIfsePJRSup2sn+6myLoQPEjWovFXx9nHO/LX/iBZf8lfSCkd8zTmFH9D9g/iofy9mNoKfDw3Ag/HKV74cjQppe+R9SG8D/gu2ffs+6ewyRP+jeUJwjuB/1j0XVhL1ipZT/a6t+bzAH8JtEbENrKL2h44hTj1IliHT8s63Dr8WEpZhz9M9pu4F/hYev4i0kmfAG5IKf2Q7Lv98fyzOdpvblIt8L8i69b3feDPpzuInG3i5D5HzUYRcTdZR/9pxwrV6Sv/blydUnqqwqFIOgrrcB3NbK/DI+KjZBdC/mmlY5ktbLGWJEmSSsAWa6mKRXYl9lfmwukzSdKRrMPnHhNrSZIkqQTsCiJJkiSVwJwbx/qKK65IX/va1yodhiTNtFk3YLb1saTT2LR18pxrsX7uuRMdO12SVE7Wx5J0pDmXWEuSJEmzkYm1JEmSVAIm1pIkSVIJmFhLkiRJJWBiLUmSJJWAibUkSZJUAibWkiRJUgmYWEuSJEklYGItSZIklUBZb2keEVcAnwRqgU+llD4+ZflZwGeATmAf8IsppR3ljEk6USMj4zy5r58Dh8bZ0z9MV3sj5y5poa2pUOnQpKo2MjLO9uf66Rsap6d/mCXtjZy/pI3mpoZKhyZJx1S2xDoiaoEbgDcAO4D7I2JTSunRotX+FPhcSummiHg98IfAL5UrJulEjYyMc9+Pn+PZ/cNs3LSVodEJCvU1XL9+DW9a02lyLZXJyMg4P3qul7HxIKVERFb+xN5+zl5oci1pditni/U6YHtK6QmAiLgZ2AAUJ9bnAR/Mp+8CvlLGeKQT9sizvdRG7eGkGmBodIKNm7awctE61q0ysZbK4ZkD/fx4/zD7Bka57vaig9oNa2hprGWVifVpb2IisX1PH0/uPURTfS2tjXW8fHGrB12aFcrZx3oZ8EzR/I68rNhDwM/m028H2iJi4dQNRcQ1EbE5Ijb39PSUJVip2K6+IXr6hw8n1ZOGRifo7huuUFRS9dt7cJzRsXQ4qYb8oPa2LfT0j1c4OlXaxETi3id62N03zOBI9v3oHxrlkV0HODQ4UuHopMpfvPjbwE9FxPeBnwJ2Ai+oOVNKN6aU1qaU1nZ2ds50jDoNLWxpoLOtkUL9kT+RQn0Ni9saKxSVVP26+4cZGB6b/qC234Pa092OAwOMp0R9bQ01Qfa3Bmpranj6wKFKhyeVNbHeCawoml+elx2WUno2pfSzKaWLgd/Nyw6UMSbphDQ31DIyPs71688/nFxnfazPp72ptsLRSdWrq72RlkLdtAe1Xe0e1J7uRsfGaKyr47mDI3S1F2gr1DG/qY6mhjr6Bj2jocorZx/r+4HVEbGKLKG+Evj3xStExCJgX0ppAvgw2QghUsWdvaiJ//vEEGctauKm96xjT/9w1lId48xvqvSJHql6nbekhf2HRvjIW89/QR/rM+d7UHs6GxkZZ+uuAbbvOchEgu17+jlncSttjbWsXFjLEz2e0VDllS2xTimNRcT7ga+TDbf3mZTS1oi4HticUtoEXAb8YUQk4B7gfeWKRzoZbU0FfvLs+fxw9wA9AyN0NNdTqK+hvVDL4va2SocnVa3W/Lf31HOD/O3VP8HegRG62hpZ0lHLotb2SoenCtrW3cfO/YPceM8Thw+4PnD5al61agE9/eMs8YyGZoGyNr2llO5IKb00pXROSum/52Ub86SalNKtKaXV+Tq/nFLycFOzRntTgVeetYCXLWmjvjZoK9Rz1sJ51NREpUOTqlp7U4Hzls6juaGW+pqgoa6GpW3t1NV5tuh01js4yifvfPyIi1o/eefjHBweo7t/mAuWzqtwhFKZbxAjzXU1NcHZna2c3dla6VCk00pdXQ0XrpjPhSuOv65OD4dGxqe9qPXQyDhd7Y0UCqY0qjwP/yVJ0qzX2Tr9SE2LWhs5b0lLhaKSjmRiLUmSZr3VXc38/tvWHDFS0++/bQ0v6Wqm1bvhapbwvIkklUnv4BBP9BxiT98wTQ3ZHeJesriZdpMA6aS1NxV443mdnLlgHd19w3S1N/LyJS3+njSrmFhLUhn0Dg6xu3eI0fHEyHiio66WhrrggacPcMlZHSYD0ovQ3lRg3Sp/O5q9TKwlqQwOHBrikZ39/N5XthweGuz337aGha31PPncIBeuMDmQpGpjH2tJKoOe/onDSTVkoxf83le20Fhby6ER7xAnSdXIFmtJKoPnDg5POzRY39AY4xOpQlFJksrJxFqSyqC5oZZCfc0RyXWhvoamhloavdGJJFUla3dJKoOO5gY+cPnqI4YG+8Dlq2luqGVJu/2rJaka2WItSWVw3pJ2nto7wDWvO5uJBDUBy+Y3kZjgzAXezEKSqpGJtSSVQV1dDW8+fylb5/eys3eQRS2NdM1rZMX8FmpqotLhSZLKwMRaksqkrq6GC8+cz4XMr8j+I+IpoB8YB8ZSSmsjYgHwBWAl8BTw8yml/RERwCeBNwOHgKtTSt+rRNySNFfZx1qSqtu/TSldlFJam89fC9yZUloN3JnPA7wJWJ0/rgH+csYjlaQ5zsRakk4vG4Cb8umbgLcVlX8uZe4FOiJiaSUClKS5ysRakqpXAr4REQ9ExDV5WVdKaVc+vRvoyqeXAc8UPXdHXnaEiLgmIjZHxOaenp5yxS1Jc5J9rCWper0mpbQzIhYD34yIHxQvTCmliDipu9WklG4EbgRYu3atd7qRNKuNjIzz8LO97O4bYml7gQvOmEdDQ23Z9mdiLUlVKqW0M/+7JyK+DKwDuiNiaUppV97VY0+++k5gRdHTl+dlkjQnjYyMs617P13tdUxMNPJs7xDjKfHSJS10NJXnfgJ2BZGkKhQRLRHRNjkNvBHYAmwCrspXuwq4LZ/eBLw7MpcCvUVdRiRpztnZ209zQx2PPjvInv5hDg2Psad/mO/+aD8HBofKsk9brCWpOnUBX85G0aMO+PuU0tci4n7gloh4L/A08PP5+neQDbW3nWy4vffMfMiSVDo9/eMU6mvYd2iU627fytDoBIX6Gj7y1vP58XODdKwofau1ibUkVaGU0hPAhdOU7wUun6Y8Ae+bgdAkaUZ09w+zsLXhcFINMDQ6wXW3b+Vvr/6JsuzTriCSJEmqOl3tjewbGDmcVE8aGp1g38BIWfZpi7UkSZKqznlLWqgJKNTXHJFcF+prWNzWWJZ92mItSZKkqtPaVOAlXS1cv2ENhfos5S3U13D9hjWcu6SlLPu0xVqSJElVqaOpwJvP72TlwnV09w3T1d7IeUtaaC3TcHsm1pIkSaparU0F1q0qTyI9lV1BJEmSpBI4LVqs+waH6OkfYu/Bcbr7y38aQJIkSaefqk+s+waH2NU7yFPPDTE6nvK77sD+gRFefc58k2tJkiSVRNUn1jsPDDEyyrR33fnRc4NcWIa77kiSJOn0U/V9rPsGxxkYHZ/2rjsDI+MVjk6SJEnVoupbrPf0DxPBtHfd2V+mu+5IkiTp9FP1iXVXeyPB9Hfd6Wovz113JEmSdPqp+q4gL13SwvDY2AvuuvOxDWtYs2RehaOTJElStaj6FuuOpgJrlsPC1iFues869vQPs6S9kQuWzqNQqPqXL0mSpBlS1swyIq4APgnUAp9KKX18yvIzgZuAjnyda1NKd5Q6jo6mAh0OqydJkqQyKltXkIioBW4A3gScB7wrIs6bstrvAbeklC4GrgT+Z7nikSRJksqpnH2s1wHbU0pPpJRGgJuBDVPWSUB7Pj0PeLaM8UiSJEllU86uIMuAZ4rmdwCvmrLOR4FvRMSvAy3AT5cxHkmSJKlsKj0qyLuAz6aUlgNvBv4uIl4QU0RcExGbI2JzT0/PjAcpSZIkHU85E+udwIqi+eV5WbH3ArcApJS+AxSARVM3lFK6MaW0NqW0trOzs0zhSpIkSS9eORPr+4HVEbEqIhrILk7cNGWdHwOXA0TEuWSJtU3SkiRJmnPKllinlMaA9wNfB7aRjf6xNSKuj4j1+Wq/BfxKRDwE/ANwdUoplSsmSZIkqVzKOo51Pib1HVPKNhZNPwq8upwxSJIkSTOh0hcvSpIkSVXBxFqSJEkqARNrSZIkqQRMrCVJkqQSMLGWJEmSSsDEWpIkSSoBE2tJkiSpBEysJUmSpBIwsZYkSZJKwMRakiRJKgETa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQRMrCVJkqQSMLGWJEmSSsDEWpIkSSoBE2tJkiSpBEysJUmSpBIwsZYkSZJKwMRakiRJKgETa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQRMrCVJkqQSMLGWJEmSSsDEWpKqWETURsT3I+If8/lVEfHdiNgeEV+IiIa8vDGf354vX1nJuCVpLjKxlqTq9gFgW9H8HwGfSCm9BNgPvDcvfy+wPy//RL6eJOkkmFhLUpWKiOXAW4BP5fMBvB64NV/lJuBt+fSGfJ58+eX5+pKkE2RiLUnV638AvwNM5PMLgQMppbF8fgewLJ9eBjwDkC/vzdc/QkRcExGbI2JzT09POWOXpDnHxFqSqlBE/AywJ6X0QCm3m1K6MaW0NqW0trOzs5SblqQ5r67SAUiSyuLVwPqIeDNQANqBTwIdEVGXt0ovB3bm6+8EVgA7IqIOmAfsnfmwJWnussVakqpQSunDKaXlKaWVwJXAt1JK/wG4C3hnvtpVwG359KZ8nnz5t1JKaQZDlqQ5r6yJdURcERGP5cM3XTvN8k9ExIP544cRcaCc8UiS+BDwwYjYTtaH+tN5+aeBhXn5B4EX1NmSpGMrW1eQiKgFbgDeQHaBzP0RsSml9OjkOiml/1y0/q8DF5crHkk6XaWU7gbuzqefANZNs84Q8HMzGpgkVZlytlivA7anlJ5IKY0AN5MN53Q07wL+oYzxSJIkSWVTzsT68NBNueJhnY4QEWcBq4BvHWW5wztJkiRpVpstFy9eCdyaUhqfbqHDO0mSJGm2K2diPTl006TiYZ2muhK7gUiSJGkOK2difT+wOiJWRUQDWfK8aepKEfFyYD7wnTLGIkmSJJVV2RLr/OYD7we+DmwDbkkpbY2I6yNifdGqVwI3O16qJEmS5rKy3nkxpXQHcMeUso1T5j9azhgkSZKkmTBbLl6UJEmS5jQTa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQRMrCVJkqQSMLGWJEmSSsDEWpIkSSoBE2tJkiSpBEysJUmSpBIwsZYkSZJKwMRakiRJKgETa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQTqKh2AJOn0MTY2wQ+6+9h/aJSh0XHOXtjCqs5Wamqi0qFJ0imzxVqSNCPGxibY1d/H+ESitiYYHpugd2iEXX29TEykSocnSafMFmtJ0ox47mAfj+4cZN+hUa67fStDoxMU6mv42IY1vPYldSzpaK10iJJ0SmyxliTNiB/vG2d0PB1OqgGGRif4b7dt4cf7hyscnSSdOlusJUkzort/mEPDY4eT6klDoxN095lYS5r7bLGWJM2IrvZGWgp1FOqP/NdTqK+hq62xQlFJUumYWEuSZsRLl7RQXxN85K3nH06uJ/tYr1hQW+HoJOnU2RVEkjQjOpoKvOqc+Tz93CCfufon2D8wwpL2Rha21tLZ2l7p8CTplJlYS5JmTEdTgdalDWzd1Uv/YFBXW8Pyjnbq6k6fE6j9g0Ns2z1Ad98wXe2NnLukhbamQqXDklQCJtaSSmpgcJituw8eThqWz29kSXuLNwDRYXV1NVy4Yj4Xrqh0JDPvwOAQPyxKqhe21LJj/yGWg8m1VAVMrCWVzMDgMF/dsoeNm7YcHqP4+vVruGDFKC9b3GFyrdPawOAw39jSwy2bn+bd/+ZsnnxugLHxJs7oqGfvwJCJtVQFTp9zb5LKbuvug4eTasiGUdu4aQv9hyZ4eu9AhaOTKmvr7oPcsvlp3vHKM/mdWx/iQ198hP9402bue6qf8YnjP1/S7GdiLalkuvuGpx+juH+YZ3sHKxSVNDt09w3z7n9zNtf945E3yNl42xb2DYxXODpJpWBiLalkutobjzpGcf/QWIWikmaHrvZGBke8QY5UzUysJZXMwpZarl9/5BjF168/n4WttXQ011c4Oqmyls9vZMX8pmkPPhd5gxypKnjxoqSSqa+Dl3U187n3rKO7f4iutgKNdYnnBkZoqjex1ultSXsLvYOjXL9hDRtvK77A93xaG71BjlQNTKwllcyCpgLDo0MMj48D2QgghYY6xoZGWbWoqbLBSRVWUxOc2dFM/9A4n75qLfsGRuloqmdsYpyzFzRXOjxJJVDWriARcUVEPBYR2yPi2qOs8/MR8WhEbI2Ivy9nPJLKq7WpwNL2I4cMm0iwuquNdocSk2hpamTNknbqa7N/v431NbzqrAU0NzVUODJJpVC2FuuIqAVuAN4A7ADuj4hNKaVHi9ZZDXwYeHVKaX9ELC5XPJJmRmtTgXWrTKKlo2luamDdqoWVDkNSGZSzxXodsD2l9ERKaQS4GdgwZZ1fAW5IKe0HSCntKWM8kiRJUtmUs4/1MuCZovkdwKumrPNSgIj4NlALfDSl9LWpG4qIa4BrAM4888yyBCtJ5TI0NMYju3rZ3TfM4rZGOpprOWdRO3V1DswkSdWk0hcv1gGrgcuA5cA9EXFBSulA8UoppRuBGwHWrl2bZjpISXqxhobG2LRl15GjQGxYw7O9Q7z2nMUm15JURcpZo+8EVhTNL8/Liu0ANqWURlNKTwI/JEu0JakqbNndeziphufvtNdYV8sjzx44zrMlSXNJORPr+4HVEbEqIhqAK4FNU9b5CllrNRGxiKxryBNljEmSZtTuo9zmff/AKM8eGKpQVJKkcihbV5CU0lhEvB/4Oln/6c+klLZGxPXA5pTSpnzZGyPiUWAc+C8ppb3likmSZtrClgYK9TVHJNeF+hrmtzi8miRVm7L2sU4p3QHcMaVsY9F0Aj6YPySp6jQ31PKRt57PdbdvPdzH+iNvPZ+W+trJe+hIkqpEpS9elKSqtrgt2N1bz5++80IGRsZoaaijvjZobwo6Wx3vW5KqiYm1JJXRwuY2Xr4UuvvG6e6DxW2NLOmoZV5TgVbvRilJVeW4Fy9GxK9HxPyZCEaSqk1DQy1L29qoiaAmoKYmWNLaxrwyJ9URUYiI+yLioYjYGhHX5eWrIuK7EbE9Ir6QX1xORDTm89vz5SvLGqAkVaGssK4dAAAeSUlEQVQTGRWki+x25LdExBURYa9ASToJDQ21rF25gLe84gzWrlxAQ0PtTOx2GHh9SulC4CLgioi4FPgj4BMppZcA+4H35uu/F9ifl38iX0+SdBKOm1inlH6PbGzpTwNXA49HxB9ExDlljk2S9CKlzMF8tj5/JOD1wK15+U3A2/LpDfk8+fLLbUiRpJNzQuNY56N37M4fY8B84NaI+OMyxiZJOgURURsRDwJ7gG8CPwIOpJTG8lV2AMvy6WXAM5ANlwr0Agun2eY1EbE5Ijb39PSU+yVIUsn1Dw5x35N7uf2hZ7nvyb30D5bungIn0sf6AxHxAPDHwLeBC1JK/wm4BHhHySKRJJVUSmk8pXQR2Z1v1wEvL8E2b0wprU0pre3s7DzlGCVpJh0YHGLb7gG6+4bpam9kXlMtP9zTX7Lk+kRGBVkA/GxK6eniwpTSRET8TEmikCSVTUrpQETcBfwk0BERdXmr9HJgZ77aTmAFsCMi6oB5gDfsklQ1+geH+O6P9jM6njg0PMaeftg/ECxqq2NP/xBtJbio/LiJdUrpI8dYtu2UI5AklVxEdAKjeVLdBLyB7ILEu4B3AjcDVwG35U/ZlM9/J1/+rbwboCRVhSeeG2TfodEX3LBrSXuBg0Olqe5OqI+1JGnOWQrcFREPA/cD30wp/SPwIeCDEbGdrA/1p/P1Pw0szMs/CFxbgZglqWwGRsYPJ9UAQ6MTXHf7VgZGxxkYHS/JPrxBjCRVoZTSw8DF05Q/Qdbfemr5EPBzMxCaJFXEvoGRw0n1pKHRCfYNjFCqIZBMrCVJklT1lnc0UaivOSK5LtTX0NxQy/zm+pLsw64gkiRJqnprzpjHxzasoVCfpb+F+hquW38+o+MTXHBGR0n2YYu1JEmSql5dXQ1vv2gZqxe3svPAIAtaGuhormN1Zzt1daVpazaxliRJ0mmhrq6Gi86cz0Vnzi/P9suyVUmSJFXc0NAYj+zqZXffMEvaG7lg6TwKBdO/cvGdPQ1MTCR29g6wa/8w3f3ZnYbWLGmjuamh0qFJkqQy6R0c4rGiuwy2FmrZ/lwvL1lkcl0uXrxY5SYmEo/39PLs/mG6+4dY3NZId98h/vmxHg4NjlQ6PEmSVAYDg8Ps7h2ksa6Gha0N7OkfZmB4nKaGWnb3H6x0eFXLw5Uqt7N3gId+3M/GTVsO32Xo+vVrmNcEP9hzkFeetaDSIUqSpBJ77tAgT+0dYt/AkXca/NiGNaw9q63S4VUtW6yr3O4Dw4eTasgGQt+4aQsdzQX6h8cqHJ0kSSqHnv4JRsfSC+40+N9u20JPf2nuMqgXssW6yu3uG572LkM9/cM01HpcJUlSNXru4DADw2PT5gDd/cMViqr6mVhXucVtjdPeZWhRWyOFulLdwFOSJM0mzQ21jE2kaXOArvbGCkZW3WyyrHItjbVct/78F9xlqLWhlqYGP35JkqrR/OYGDg2P8pG3HpkDfGzDGs5b0lLh6KqXLdZV7uwFzfT0D3HjL13C/kOjzG+uZ3xigsZ6aPDTlySpKp27pJ0n9w7Qd2iEv/6lSzhwaJTOtkbWnNFKa1Oh0uFVLZssq1xzUwMXn9VBob6WmggKdbUsn19gIiU6m5srHZ4kSSqDuroa3nz+Ui5Y3sHQyDgrFzaz7qwFtJtUl9Vp0WY5NjbBo7t72XVgiKaGWuYV6jl3STsNDbWVDm1GdDQVWLMkO4bq7huG8AYxkiRVu7q6Gi5cMZ8LV1Q6ktNH1SfWY2MT/NOWXfyXLz58eAzHP37HKxgaH+XiZQtPm+S6uamBdasWVjoMSZKkqlX1XUEe6z5wOKmGbJiZ3/niwxTq6nhyX3+Fo5MkSVK1qPrEesf+6cdx7u4dofeQA6RLkiSpNKo+sW5qqD08zMykQn0NjQ01DpAuSZKkkqn6xLq1sY4PXL76iDEcP3D5atoa6ljiAOmSJEkqkaq/ePGMeTUs62jimtedzUSCmoBlHQXmt9SwtK210uFJkiSpSlR9Yr2guY11K2HJvALdfcN0tTWybEEtHY1NFApV//IlSZI0Q6q+K0hDQy0LmtuoiaAmoKYm6Gxuo6XJbiCSJEkqndOiybahoZa1KxdUOgxJkiRVsapvsZYkSZJmQlkT64i4IiIei4jtEXHtNMuvjoieiHgwf/xyOeORJEmSyqVsXUEioha4AXgDsAO4PyI2pZQenbLqF1JK7y9XHJIkSdJMKGeL9Tpge0rpiZTSCHAzsKGM+5MkSZIqppyJ9TLgmaL5HXnZVO+IiIcj4taIWDHdhiLimojYHBGbe3p6yhGrJEmSdEoqffHi7cDKlNIrgG8CN023UkrpxpTS2pTS2s7OzhkNUJIkSToR5UysdwLFLdDL87LDUkp7U0rD+eyngEvKGI8kSZJUNuVMrO8HVkfEqohoAK4ENhWvEBFLi2bXA9vKGI8kSZJUNmUbFSSlNBYR7we+DtQCn0kpbY2I64HNKaVNwG9ExHpgDNgHXF2ueCRJkqRyKuudF1NKdwB3TCnbWDT9YeDD5YxBkiRJmgmVvnhRkiRJqgom1pIkSVIJmFhLkiRJJWBiLUmSJJWAibUkSZJUAibWkiRJUgmYWEuSJEklYGItSZIklYCJtSRJklQCJtaSJElSCZhYS5IkSSVgYi1JkiSVgIm1JEmSVAIm1pIkSVIJmFhLkiRJJWBiLUmSJJWAibUkSZJUAibWkiRJUgmYWEtSFYqIFRFxV0Q8GhFbI+IDefmCiPhmRDye/52fl0dE/HlEbI+IhyPilZV9BZI095hYS1J1GgN+K6V0HnAp8L6IOA+4FrgzpbQauDOfB3gTsDp/XAP85cyHLElzm4m1JFWhlNKulNL38ul+YBuwDNgA3JSvdhPwtnx6A/C5lLkX6IiIpTMctiTNaSbWklTlImIlcDHwXaArpbQrX7Qb6MqnlwHPFD1tR14mSTpBJtaSVMUiohX4IvCbKaW+4mUppQSkk9zeNRGxOSI29/T0lDBSSZr7TKwlqUpFRD1ZUv35lNKX8uLuyS4e+d89eflOYEXR05fnZUdIKd2YUlqbUlrb2dlZvuAlaQ4ysZakKhQRAXwa2JZS+rOiRZuAq/Lpq4DbisrfnY8OcinQW9RlRJJ0AuoqHYAkqSxeDfwS8EhEPJiX/Vfg48AtEfFe4Gng5/NldwBvBrYDh4D3zGy4kjT3mVhLUhVKKf0rEEdZfPk06yfgfWUNSpKqnF1BJEmSpBIwsZYkSZJKwMRakiRJKgETa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQTKmlhHxBUR8VhEbI+Ia4+x3jsiIkXE2nLGI0mSJJVL2RLriKgFbgDeBJwHvCsizptmvTbgA8B3yxWLJEmSVG7lbLFeB2xPKT2RUhoBbgY2TLPex4A/AobKGIskSZJUVuVMrJcBzxTN78jLDouIVwIrUkpfPdaGIuKaiNgcEZt7enpKH6kkSZJ0iip28WJE1AB/BvzW8dZNKd2YUlqbUlrb2dlZ/uAkSZKkk1TOxHonsKJofnleNqkNWAPcHRFPAZcCm7yAUZIkSXNRORPr+4HVEbEqIhqAK4FNkwtTSr0ppUUppZUppZXAvcD6lNLmMsYkSZIklUXZEuuU0hjwfuDrwDbglpTS1oi4PiLWl2u/kiRJUiXUlXPjKaU7gDumlG08yrqXlTMWSZIkqZy886IkSZJUAibWkiRJUgmYWEuSJEklYGItSZIklYCJtSRJklQCJtaSJElSCZhYS5IkSSVgYi1JkiSVgIm1JEmSVAIm1pIkSVIJmFhLkiRJJWBiLUmSJJWAibUkSZJUAibWkiRJUgmYWEuSJEklYGItSZIklUBdpQOQJJ0+JiYST+0doLtviK72AisXtlBTE5UOS5JKwhZrSdKMmJhI7OztpYYxaiPY+mwfm5/ex4HBoUqHJkklYYu1JGlG7Ok7SO+hcR599hAbN21haHSCQn0N169fwxvXdNLRVKh0iJJ0SmyxliTNiB/vH+HQCIeTaoCh0Qk2btrCD3YdZGxsosIRStKpMbGWJM2I7r5huvuGDyfVk4ZGJ+juG2brrt4KRSZJpWFiLUmaEV3tjSxub6RQf+S/nkJ9DR3N9ezuta+1pLnNxFqSNCOWzqulpSHrUz2ZXBfqa7hu/fl8/t4nWTLPPtaS5jYvXpQkzYjO5lbGxg+yqrPAZ9/zE3T3DdPRXM/n732SN5x3BucvnVfpECXplJhYS5JmRKFQx1Ja6R/qZ7wmsbitkX0DI/zaZS9lzRnzqKvzJKqkuc3EWpI0YwqFOl6xYn6lw5CksrB5QJIkSSoBW6wlqcQGBofZ1n2Q3X3DLGxpoLmhlpcuaqWpqb7SoUmSysjEWpJKaGBwmK9u3cPG256/s+BH3no+u3uHeN05i0yuJamK2RVEkkpo6+6Dh5NqyG5+ct3tWxkdTzyyu6/C0UmSyskWa0kqoaPdWXBgZIxu82pJqmq2WEtSCXUd5c6CLQ11dLU1VigqSdJMMLGWpBJa2FLL9RuOvLPgR956PvW1wbzm2gpHJ0kqJ7uCSFIJtTTU8rLFTfyv966ju2+Y+S0NtNTXUlszwTJv2S1JVc0Wa0kqocXtrSxoq2MiQUpQE0FtbXDWwhZam0ysJamalTWxjogrIuKxiNgeEddOs/xXI+KRiHgwIv41Is4rZzySVG41NcGyefNY1NrIorYGFrc1ct7SebSZVEtS1StbYh0RtcANwJuA84B3TZM4/31K6YKU0kXAHwN/Vq54JGmm1NQEZ3e2cunZizi7s5WampjxGCLiMxGxJyK2FJUtiIhvRsTj+d/5eXlExJ/njSAPR8QrZzxgSaoC5WyxXgdsTyk9kVIaAW4GNhSvkFIqHnyqBUhljEeSTiefBa6YUnYtcGdKaTVwZz4PWQPI6vxxDfCXMxSjJFWVcibWy4BniuZ35GVHiIj3RcSPyFqsf6OM8UjSaSOldA+wb0rxBuCmfPom4G1F5Z9LmXuBjohYOjORSlL1qPjFiymlG1JK5wAfAn5vunUi4pqI2BwRm3t6emY2QEmqHl0ppV359G6gK58+oYYQsD6WpGMpZ2K9E1hRNL88Lzuam3m+9eQIKaUbU0prU0prOzs7SxiiJJ2eUkqJF9H9zvpY0lw0NDTG/U/u5faHnuW7T+zlsd0HGBubOP4TT1I5x7G+H1gdEavIEuorgX9fvEJErE4pPZ7PvgV4HElSuXRHxNKU0q68q8eevPxkG0Ikac4YGhrjGz/oZnvPQSYSbN/TzzmdLfQcHObSlZ3U1ZWunblsLdYppTHg/cDXgW3ALSmlrRFxfUSsz1d7f0RsjYgHgQ8CV5UrHkkSm3i+nr0KuK2o/N356CCXAr1FXUYkaU7b1t3HzgOD3HjPE/zFt7bz1/c8wc4DQzTW1bFtd9/xN3ASynrnxZTSHcAdU8o2Fk1/oJz7l6TTVUT8A3AZsCgidgAfAT4O3BIR7wWeBn4+X/0O4M3AduAQ8J4ZD1iSyqRveIxP3vk4Q6NZ14+h0Qk+eefj/PUvXcLgyHhJ9+UtzSWpCqWU3nWURZdPs24C3lfeiCSpMgZHxg8n1ZOGRidKnlTDLBgVRJIkSSqXxW2NFOqPTHkL9TU01tWwoKWhpPsysZYkSVLVWtBcw++/bc3h5LpQX8PHNqxhZHyCxhJeuAh2BZGOa2Ii8dTeAbr7huhqL7ByYUtFblEtSZJOXldrC5ecCX/3H9fR3T/MwpYG6mph38FRzlrUVNJ9mVhLxzAxkXis+wD9gxM8d3CY2ppgeGyUl3V1mFxLkjQHNDXVs4QWRscHSKmR7v5hutobedU58+loKpR0XybW0jHs6jvIIzv62bhpK0OjExTqa7h+/fm0NdWwvGNepcOTJEknoKmpnpc2dZR9P/axlo5h576Rw0k1ZFcRb9y0lWf3jdE/OFTh6CRJ0mxiYi0dQ3f/8LRD9HT3D7Nt90CFopIkSbORXUGKDA2N8ciuXnb3DbOkvZELls6jUPAtOp11tWdD9BQn14X6GrraG+nuG+bQ4AjNTaUdqkeS9EIDg8Ns3X2Q7r6sf+z5S1ppaWqsdFjSEWyxzh0YHOLhPKnuam9kUWstuw8eZGhorNKhqYLOmF/L9euPHKLn+vVrWNhSS1dbI1t291c4Qkmqfn2DQ+zuH4T0fFn/8CADg8OVC0qahs2xZEn1N7b0sHHTliMuUDtvaQv76g5yRqH8nd01Oy1pbefiM+Fz71lHd/8QXW0F5rXUsn3PIV66pJktOw9VOkRJqmojI+P8YHc/+w+NMTqWODQ8xp5+2H8ouHAZtlprVrHFGvjh7oHDSTU8f4HaoRHYsb/0t7vU3FFXV8PYBERATQRjKTE0nFi2oEDfoURXmxW6JJXTw8/20lRfx76BUX771of40Jce4bf/90PsGxil52A2LKo0W5hYA919R7tAbYjuPk8zne5WdrSwb2CElOCZfYd4Zv8gO/cNMTQ2ysLW2kqHJ0lVbXffEAeHx7nu9iNHaLru9q0cHB7jqb1eSK7Zw8Sa5y9QK1aor6GrrWCLpGhqqufV58xncXsjzQ11LG5rZOWiZupqgxUdbZUOT5Kq2tL2AnsPjkzbALZ3YIQ9/Q59qtnDxBp46ZKWaS5QO5/mBlhgi6SA1qYCFy3rYOm8Aj0Hh+kfmuCCpQtoaPD7IUnldMEZ81h81AawRha3lfbOedKp8OJFoKOpwBvXdLJy0bpsGJ+2RjrbatlzcISXL5lf6fA0SzQ01LJ25YJKhyFJp5WGhlrOzRvAjhxkYA0N9bByYUulQ5QOM7HOdTQVeOWKBrbu6mV37xCN9QUuWbGIujob9SVJqqS2pgI/s6aLlYuaDzeALVvQyNL2FmpqotLhSYeZWBepq6vhwhXzuXBFpSORJEnFmpsaWLdqYaXDkI7J5lhJkiSpBEysJUmSpBIwsZYkSZJKwMRakiRJKgETa0mSJKkETKwlSZKkEjCxliRJkkrAxFqSJEkqARNrSZIkqQQipVTpGE5KRPQAT5/g6ouA58oYTrnM1bhh7sY+V+OGuRu7cZ+c51JKV1Rgv0eV18cDzK7PcTZ+r2ZbTMZzbLMtHph9MRnPUerkOZdYn4yI2JxSWlvpOE7WXI0b5m7sczVumLuxG3d1mG3vx2yLB2ZfTMZzbLMtHph9MRnP0dkVRJIkSSoBE2tJkiSpBKo9sb6x0gG8SHM1bpi7sc/VuGHuxm7c1WG2vR+zLR6YfTEZz7HNtnhg9sVkPEdR1X2sJUmSpJlS7S3WkiRJ0owwsZYkSZJKoCoS64i4IiIei4jtEXHtNMtfFxHfi4ixiHhnJWKczgnE/cGIeDQiHo6IOyPirErEOZ0TiP1XI+KRiHgwIv41Is6rRJxTHS/uovXeEREpImbF8D0n8H5fHRE9+fv9YET8ciXinM6JvOfx/7d3/7FXV3Ucx5+v+PIbAoR0JBRQKJlzgj+GZczE2tQGrVyj1dTNajMrqbay2VpWm1nNsuWyCTZtARaZkS2nBGarQEF+COIQhBTHD2eCpTOJ3v1xzoXPvn6/fD9fd7/3XL68Htvd/XzOvbuf1z6fc889n885917pY7mub5a0qNUZu1Jjn/+wsr+3StpfImdfk3SHpH2SNlXKTpD0oKSn8v2YXC5JP877bKOkGX2QZ6KklZX6cm3JTJKGSHpE0oac54ZcPlnS6rzduyUNyuWD8/q2/PikZuap5BogaZ2k+0rnkbSz8nmwJpcVq0N5O6MlLZX0pKQtks4rWIdOrbQl6yW9JGl+4ffZF3N93iRpca7nJevQtTnLZknzc1nROtStiDimb8AAYDswBRgEbABO6/ScScAZwF3AZaUz9yL3+4Fheflq4O7SuXuR/c2V5TnA/cdC7vy8kcDDwCrg7GMhN3Al8JPSWd9g9qnAOmBMXj/xWMjd6fmfB+4onbuP9sUsYAawqVL2PeC6vHwdcFNevgT4IyBgJrC6D/KMB2bk5ZHAVuC0Upny647IywOB1Xk7vwLm5fLbgKvz8meB2/LyvL5q14EvAYuA+/J6sTzATmBcp7JidShv507gU3l5EDC6dKa8rQHAHuDtBev0ycAOYGil7lxZqg4BpwObgGFAB7AceGc7HK+ubv3hivW5wLaIeDoiXgOWAHOrT4iInRGxEfhfiYDdqJN7ZUS8kldXARNanLE7dbK/VFkdDrTDt2R7zJ19G7gJeLWV4Y6ibu52VCf7p4FbI+JFgIjY1+KMXentPv84sLglyVosIh4G/tmpeC6pY0K+/3Cl/K5IVgGjJY1vcp7dEfFYXv4XsIXUESiSKb/uv/PqwHwL4EJgaTd5GjmXArMlqVl5ACRNAC4FFuR1lczTjWJ1SNIo0gnjQoCIeC0i9pfMVDEb2B4R/yicpwMYKqmD1KHdTbk69C5S5/iViPgv8GfgI7TH8Xqd/tCxPhl4trK+K5e1u97mvop0BtYOamWXdI2k7aSzyi+0KNvR9Jg7DxlNjIg/tDJYD+rWlY/mYa+lkia2JlqP6mQ/BThF0l8lrZLUDn/bXfv9qTRFazKwogW52sVJEbE7L+8BTsrLLW2P85DzdNJV4mKZlKZdrAf2AQ+SRjv2505A520ezpMfPwCMbWYe4EfAVzhyMWls4TwBPCBpraTP5LKSdWgy8Dzwc6XpMgskDS+cqWEeR07Si+SJiOeAHwDPkDrUB4C1lKtDm4D3SRoraRjpivRE2uN4vU5/6Fj3e5I+CZwNfL90lt6IiFsj4h3AV4Gvl87TE0lvAm4Gvlw6yxvwe2BSRJxB+mC/s4fnt5MO0nSQC0hXfm+XNLpoot6ZByyNiEOlg5QQaey15SNSkkYAvwHmdxoha3mmiDgUEWeSRhXPBaa1atudSfoQsC8i1pbK0IXzI2IGcDFwjaRZ1QcL1KEO0vSmn0bEdOBl0lSCkpnIc5bnAL/u/Fgr8+S5ynNJJyBvJY06F7vgERFbSKPIDwD3A+uBQ52eU6Qd6kp/6Fg/RzpzaZiQy9pdrdySLgKuB+ZExH9alK0nvd3nSzgyRFNST7lHkuZyPSRpJ2lu1jKV/wJjj/s7Il6o1I8FwFktytaTOnVlF7AsIg5GxA7SnNmpLcrXnd7U8eoVpuPF3sbQar5vTN9pSXssaSCpU/3LiLinHTIB5OkEK4HzSMPPHV1s83Ce/Pgo4IUmxngvMCe3YUtIw/e3FMzTuALamOb1W9LJR8njtQvYFRGr8/pSUke7dB26GHgsIvbm9VJ5LgJ2RMTzEXEQuIdUr0rWoYURcVZEzAJeJH1OlD5eXeoPHetHgan526qDSB9yywpnqqPH3JKmAz8jdarbYd5pQ53s1Y7RpcBTLczXnaPmjogDETEuIiZFxCTSvPY5EbGmTNzD6uzv6vyxOaR5p+2gzvvzXtLVaiSNI00NebqVIbtQq12RNA0YA/y9xflKWwZckZevAH5XKb88fyt/JnCgMlTbFHnu5kJgS0TcXDqTpLc0RlgkDQU+QHr/rQQav0LVOU8j52XAiny1rSki4msRMSG3YfPy63+iVB5JwyWNbCwDHyQN7RerQxGxB3hW0qm5aDbwRMlMWefvapTK8wwwU9Kw/H5r7J8idQhA0on5/m2k+dWLKH+8uhYt/KZkX91I8222kua1XZ/LvkXqFAGcQzpDfZl0FrW5dOaauZcDe0nDHutJV/WK566Z/RZgc869Enh36cx1cnd67kO0wa+C1NzfN+b9vSHv72mlM/ciu0hTcJ4AHid/67z0rU5dAb4JfLd01j7eD4tJ8ywP5nb0KtL8yT+RTpiXAydUjuWteZ893hfvH+B80pDvxkrbeEmpTKRfnFqX82wCvpHLpwCPANtIQ/uDc/mQvL4tPz6lD4/dBRz5VZAiefJ2N+Tb5sp7qVgdyts5E1iTj9u9pBPkkvV6OKl/MqpSVjLPDcCTuU7/Ahhcsk4DfyF9RmwAZpfeP0e7+S/NzczMzMyaoD9MBTEzMzMzK84dazMzMzOzJnDH2szMzMysCdyxNjMzMzNrAneszczMzMyawB1rMzMzM7MmcMfazMzMzKwJ3LE264GkcyRtlDQk/4vYZkmnl85lZna8cXts7c5/EGNWg6TvkP5daiiwKyJuLBzJzOy45PbY2pk71mY1SBoEPAq8CrwnIg4VjmRmdlxye2ztzFNBzOoZC4wARpKulJiZWRluj61t+Yq1WQ2SlgFLgMnA+Ij4XOFIZmbHJbfH1s46Sgcwa3eSLgcORsQiSQOAv0m6MCJWlM5mZnY8cXts7c5XrM3MzMzMmsBzrM3MzMzMmsAdazMzMzOzJnDH2szMzMysCdyxNjMzMzNrAneszczMzMyawB1rMzMzM7MmcMfazMzMzKwJ/g8uktprD1nA0gAAAABJRU5ErkJggg==\n",
      "text/plain": "<Figure size 720x720 with 4 Axes>"
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Plots one recorded calibration per row, identified by `timestamp`\n",
    "# Left column shows normalized coordinates, righ column shows pixel coordinates\n",
    "sns.relplot(\n",
    "    data=ref_data,\n",
    "    x=\"x\",\n",
    "    y=\"y\",\n",
    "    col=\"unit\",\n",
    "    row=\"timestamp\",\n",
    "    facet_kws=dict(sharex=False, sharey=False),\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.8.2 64-bit ('pupil': pyenv)",
   "name": "python382jvsc74a57bd005b700fd9ecfc0c5b40ebbb44ab050033c7129b8b5c85a3f55e3f873f86603ea"
  },
  "language_info": {
   "name": "python",
   "version": ""
  },
  "orig_nbformat": 2
 },
 "nbformat": 4,
 "nbformat_minor": 2
}