r1_calib_example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "from ctapipe.io.lsteventsource import LSTEventSource\n",
    "from ctapipe.io import EventSeeker        \n",
    "\n",
    "from calib_drs4 import LSTR1Calibrator"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING:ctapipe.io.eventseeker.EventSeeker:Seeking to event by looping through events... (potentially long process)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6101\n",
      "event id =  6102\n",
      "event id =  6103\n",
      "event id =  6104\n",
      "event id =  6105\n",
      "event id =  6106\n",
      "event id =  6107\n",
      "event id =  6108\n",
      "event id =  6109\n",
      "event id =  6110\n",
      "event id =  6111\n",
      "event id =  6112\n",
      "event id =  6113\n",
      "event id =  6114\n",
      "event id =  6115\n"
     ]
    }
   ],
   "source": [
    "calib = LSTR1Calibrator(pedestal_path=\"/media/pawel1/ADATA HD330/20181006/pedestal_00080000.dat\")\n",
    "path_to_data_dir = '/media/pawel1/ADATA HD330/20181006/'\n",
    "file_path_list = [path_to_data_dir + \"LST-1.Run0008.0001.fits.fz\"]\n",
    "\n",
    "nr_of_events = 15\n",
    "reader = LSTEventSource(input_url=file_path_list[0], max_events=100)\n",
    "seeker = EventSeeker(reader)\n",
    "event = seeker[0:nr_of_events]\n",
    "\n",
    "for i, ev in enumerate(event):\n",
    "    print(\"event id = \", ev.r0.event_id)  \n",
    "    calib.calibrate(ev) \n",
    "                  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from plot_drs4_functions import plot_waveform #This is function to plot waveform"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6101\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAF3CAYAAABJzllyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X+8HWV94PHPV34FTRASaBoI20SNVooYuEDAnzfIryYE1NWuVlsQMNhl98UWS8VuK6FdV9xVUNQqsCBRqtGqVCAowXijpUpoQlMCBpACQkgqGIIkRX5/948zF26S+2POzZl7zpl83q/XvO6ZOc/M+T53zjnf88w880xkJpIkqV5e0u4AJElS65ngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBqqPMFHxE4R8S8RcV0xPz0ilkfEzyPiGxGxa7F8t2L+nuL5aVXHJklSXY1FC/4sYM2A+U8CF2XmDGAjcFqx/DRgY2a+CrioKCdJkkah0gQfEVOBucD/K+YDOAr4VlFkIfD24vFJxTzF828rykuSpCZV3YL/DPDnwPPF/CTgscx8tphfC+xXPN4PeBCgeP7XRXlJktSknavacEScADycmSsjord/8SBFs8RzA7c7H5gPsPvuu/fsv//+LYi24fnnn+clL6lHv0Pr0pnqUpe61AOsS6eqS12qqMfdd9/9q8zcZ8SCmVnJBHyCRgv9fuDfgSeAvwN+BexclDkSuKF4fANwZPF456JcDPcaPT092Up9fX0t3V47WZfOVJe61KUemdalU9WlLlXUA1iRJfJwZT+PMvOjmTk1M6cB7wF+mJnvA/qAdxXFTga+Wzy+ppineP6HRUUkSVKT2nH84yPA2RFxD41z7JcXyy8HJhXLzwbObUNskiTVQmXn4AfKzGXAsuLxvcDhg5R5Enj3WMQjSVLdjUmClySpWc888wxr167lySefbHcoo/byl7+cNWvWjFxwEOPGjWPq1Knssssuo1rfBC9J6khr165lwoQJTJs2jW4dFmXTpk1MmDCh6fUykw0bNrB27VqmT58+qtfu/msQJEm19OSTTzJp0qSuTe7bIyKYNGnSdh29MMFLkjrWjpjc+21v3U3wkiQNYaeddmLmzJkceOCBzJs3j8cee+yF5xYuXMiMGTOYMWMGCxcuHGYr7WGClyRpCLvvvjurVq3i9ttvZ+LEiXzhC18A4NFHH+X8889n+fLl3HLLLZx//vls3LixzdFuyQQvSVIJRx55JA899BAAN9xwA8cccwwTJ05kr7324phjjuH73/9+myPckglektTxIqqZynruuedYunQpJ554IgAPPfQQA++FMnXq1BeSf6cwwUuSNITf/OY3zJw5k0mTJvHoo49yzDHHADDYSOqd1iHQBC9J6niZ1Uwj6T8H/4tf/IKnn376hXPwU6dO5cEHH3yh3Nq1a9l3332rqv6omOAlSRrBy1/+ci6++GI+9alP8cwzz3DcccexZMkSNm7cyMaNG1myZAnHHXdcu8PcgglekqQSDj74YF7/+tezaNEiJk6cyF/91V9x2GGHcdhhh/Gxj32MiRMntjvELThUrSRJQ9i8efMW89dee+0Lj0899VROPfXUsQ6pNFvwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJQxjudrHHH388e+65JyeccEIbIxyaCV6SpCEMdbtYgHPOOYevfvWrbYxueCZ4SVJnq+pWck3eUm7g7WIB3va2tzFhwoQqatwSJnhJkkaw9e1iu4EJXpKkIQx1u9huYIKXJGkIQ90uthuY4CVJna2qm8GXvSk8294uthuY4CVJKmHg7WIB3vzmN/Pud7+bpUuXMnXqVG644YY2R7glbxcrSdIQhrtd7D/+4z+OdThNsQUvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5LUpIsvvpjXvva1vO997+Opp57i6KOPZubMmXzjG99oajvLli3jJz/5SSUxeh28JElN+tu//Vu+973vMX36dG6++WaeeeYZVq1a1fR2li1bxvjx43nDG97Q8hhN8JIkDePCCy/kiiuuAOD000/nzjvv5N577+XEE0/k/e9/P5dddhmPPPIIM2fO5Nvf/jaXXXYZ11xzDTvvvDO9vb1cfPHFPPLII3zoQx/igQceAOAzn/kM++23H1/60pfYaaeduOqqq/jc5z7Hm9/85pbFbYKXJHWFBcsWcP6Pzi9V9oOHfJBL5126xbL5187nslsve2H+vLeex4LeBcNuZ+XKlXz5y19m+fLlZCazZs3iqquu4vvf/z59fX3svffezJo1i0996lNcd911PProo1x99dXceeedRAQPPvggAGeddRZ/+qd/ypve9CYeeOABjjvuONasWcOHPvQhxo8fz5/92Z81988owQQvSdIQbrrpJt7xjnfwspe9DIB3vvOdww5Ru8ceezBu3DhOP/105s6dy1vf+lYAfvCDH/Czn/3shXKPP/44mzZtqjR2E7wkSUPIkneb67fzzjtzyy23sHTpUhYtWsRnP/tZfvSjH/H888/z05/+lN13372iSAeJZcxeSZKk7bCgd8GIh9SHc+m8S7c5bD+St7zlLZxyyimce+65ZCZXX301X/3qV7nwwgsHLb9582aeeOIJ5syZwxFHHMErX/lKAI499lg+//nPc8455wCwatUqZs6cyYQJE3j88cdHXafheJmcJElDOOSQQzjllFM4/PDDmTVrFqeffjoHH3zwkOU3bdrECSecwEEHHcRb3/pWPvGJTwCNy+pWrFjBQQcdxAEHHMCXvvQlAObNm8fVV1/NzJkzW353OlvwkiQN4+yzz+bss8/eYtn999//wuPe3l56e3sBmDJlCrfccssLz/WfZ997770HvUb+1a9+Nbfddlvrg8YWvCRJtWSClySphkzwkiTVUGUJPiLGRcQtEfGvEXFHRJxfLL8yIu6LiFXFNLNYHhFxcUTcExG3RcQhVcUmSeoOzV6mVifbW/cqO9k9BRyVmZsjYhfgpoj4XvHcOZn5ra3K/z4wo5hmAV8s/kqSdkDjxo1jw4YNTJo0iYhodzhjKjPZsGED48aNG/U2Kkvw2fjpsbmY3aWYhvs5chLwlWK9myNiz4iYkpnrq4pRktS5pk6dytq1a3nkkUfaHcqoPfnkk6NO0uPGjWPq1Kmjfu2o8vBHROwErAReBXwhMz8SEVcCR9Jo4S8Fzs3MpyLiOuCCzLypWHcp8JHMXLHVNucD8wEmT57cs2jRopbFu3nzZsaPH9+y7bWTdelMdalLXeoB1qVT1aUuVdRj9uzZKzPz0BELZmblE7An0AccCEwBAtgNWAh8rCizGHjTgHWWAj3Dbbenpydbqa+vr6Xbayfr0pnqUpe61CPTunSqutSlinoAK7JE7h2TXvSZ+RiwDDg+M9cXMT4FfBk4vCi2Fth/wGpTgXVjEZ8kSXVTZS/6fSJiz+Lx7sDRwJ0RMaVYFsDbgduLVa4B/rjoTX8E8Ov0/LskSaNSZS/6KcDC4jz8S4BvZuZ1EfHDiNiHxmH6VcCHivLXA3OAe4AngA9UGJskSbVWZS/624BtRuTPzKOGKJ/AmVXFI0nSjsSR7CRJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNbRzuwOorYgt5zPbE4ckaYdkC16SpBoywUuSVEMm+LESMfz0ilfAjTe2O0pJUk2Y4DvFfffBGWe0OwpJUk2Y4KsyfXrz69x3ny19SVJLmOCrcsklo0vyI7nvPjj22BF/CMz6wz8c+x8CN97Y+AEy0o8Uf6xIUuVM8FU55hi4997G5XHDTVX8CAB2X79+7A/5n3FG4wdIWZ6WkKTKmODbraqWPjSXbNv1emMdo9Qsj0ypS5ng261sS3/JktH9EPDLSBpckbh7Z88e/jNy7LHNH5kqcRqtqcnPqUbBBN8tyv4QaFaVh8lbFaPUjLIt7mYTdzt5OmtbzR5Z6Z92IJUl+IgYFxG3RMS/RsQdEXF+sXx6RCyPiJ9HxDciYtdi+W7F/D3F89Oqiq3WRtt7X6qLZvuCdIvh6lT2aESdjgjUdT+3UJUt+KeAozLz9cBM4PiIOAL4JHBRZs4ANgKnFeVPAzZm5quAi4pyalaV5/Sr0olXBLTajvgF3GplW2xVfelPn944Vdbq02itMJoOrq0+jTDW72uT+4gix+CQaUS8FLgJ+BNgMfDbmflsRBwJLMjM4yLihuLxTyNiZ+DfgX1ymAAPPfTQXLFiRcviPOXKU1j4i4Wlyn7wkA9y6bxLt1g2/9r5XHbrZaXWP++t57Ggd8EWy+Z9fR7X3X1dqfUvOeES5vfM32JZz6U93Lr+1lLrX/M1mHf3gAWZ7PvpfVm/eX2p9Vd8cAU9+/ZssSzOL3/466FPw76bXpxfNwH2+3Dp1cnztnxbrFy3kkMvO7TUulPGT2Hdh9dtsezau67lxEUnllr/kCmHsHL+ysbMjTfCGWdw6cT7OGNeqdU54S649utbLlvQC+f3llv/gyvh0muLmenT4ZJLmP/k34/de+8nk5i/ZMMWy3rmw637llqda95zDfNeM+CfdeON7PuD32f9S58rtf6KS6Bnq7dpLCj32rDte+/fXvlbvOqPHi69flvfe+tg5ZZfO1zaQ3vee4X58+CynsHLb+28ZbDgF9MbpxsLo/7eKw63N/Xe+xrMu2vL/Vf1997Jv3MyV55yZbkAS4qIlZk54puu0rvJRcROwErgVcAXgH8DHsvMZ4sia4H9isf7AQ8CFMn/18Ak4FdbbXM+MB9g8uTJLFu2rGXxPv3006XLrlu/bpvXXrd+3eCFB3H//fdvs/6GX20YvPAg7rrrLpZt2nL9TZs2DV64hGXLljVV/xUrV7Dp7tG/3vba+n9316a7Sq/79NNPb7P+6l+tLr3+pk2bXlh/1sknNy5JnFh69da67z5+c/LJrFtwSOlVtve9x4Ymyg5i9erVTFg/4YX5WSefDH9YLrm3wk/+/u/Ze7e9X5j/xaO/gNWnlF6/ne+9Wrjvvi3+B6P93usdxUuv/vjHmbDV/7/q773B9vlYqTTBZ+ZzwMyI2BO4GnjtYMWKv4P9DNqm9Z6ZlwKXQqMF39vb25pggSvvv7J02X2n7MvWr/21TV+Dcj8EmTZt2jbrT1o/CR4tt/5rXvMaenu2XH/C3RNgc7n1t9bb28uuK3eFku/1Q3sO3eaXLD9q4gX33x9+9mATK2xp6//dhHUToNzBC3bddddt1t901ya4o9z6EyZMeHH99SV3eIV2X7+efa9bDCVbUdv73tter/uf/5Peu0cuV5U3vOEN7DvhxSbft274VlPrt/O9t92OmAU3P9z2w9sD/wfb+73XjNe97nX0vmbL9av+3htsn4+ZzByTCTgPOIdGi3znYtmRwA3F4xuAI4vHOxflYrht9vT0ZCv19fW1dHvtNGJdRu6Pv+00fXrmkiXlt7m9ym5vyZJGbKOp0/bUd7AYR5iemDJl+G1WVZf++myvKuKqYr80oaM/902+H0Z8f7XLKPbziHVp9fdNRap4fwErMkfOu5W14CNiH+CZzHwsInYHjqbRca4PeBewCDgZ+G6xyjXF/E+L539YVESdov9SnQHnz8bUWF/i0t8RqRkjvGWXL1s2/K/5/sshR/KKVzTfCuu/10ErDVff0cQ4fXr73l+dqOz7oTDi+6uL7L5+ffOfP22hyl70U4C+iLgN+Gfgxsy8DvgIcHZE3EPjHPvlRfnLgUnF8rOBcyuMTaPt7WvP1c7QDVdLNBtj0WFQNdTp79WaqqwFn5m3AQcPsvxe4PBBlj8JvLuqeLSVSy7p/OtIp0/v7Pi2NpZfYlW29Msaqb5Ntj5VY1V/3/gDYlCOZLejKjsyXjvPkoymBTjcdcplp9Fcz9yprc+KWvq/mTKlM+urztTM902zn79O/ex1gEp70UvbpfhSWDbW5xXr1PKsqC51OterDtOuz30N2YKXJKmGbMGreTvYDRskqRvZgpckqYZM8BrZaDpp2atVktrKBK+ReT2zJHUdz8FrZHXqVS5JOwhb8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINDdmLPiLOLrH+f2Sm10NJktRhhmvBnwOMByYMM3246gAlSVLzhrsO/quZ+dfDrRwRL2txPJIkqQWGbMFn5p+PtHKZMpIkaeyN2MkuIs6KiD2i4fKIuDUijh2L4CRJ0uiU6UV/amY+DhwL7AN8ALig0qgkSdJ2KZPg+2/+PQf4cmb+64BlkiSpA5VJ8CsjYgmNBH9DREwAnq82LEmStD3K3E3uNGAmcG9mPhERk2gcppckSR2qTAv+xsy8NTMfA8jMDcBF1YYlSZK2x3Aj2Y0DXgrsHRF78eJ59z2AfccgNkmSNErDteDPAFYCv1v87Z++C3yh+tA0WnPnwuzZvUQw4jR3brujlSRVYcgWfGZ+FvhsRPz3zPzcGMak7XT99dWUlSR1jxE72WXm5yLiDcC0geUz8ysVxqUWyBz++fBiR0mqrRETfER8FXglsAp4rlicgAlekqQOVeYyuUOBAzJHag9KkqROUeYyuduB3646EEmS1DplWvB7Az+LiFuAp/oXZuaJlUUlSZK2S5kEv6DqICRJUmuV6UX/o7EIRJK0rblzy1/OOmvW67j55mrjUfco04t+E41e8wC7ArsA/5GZe1QZmCSpubEqli+fVF0g6jplWvATBs5HxNuBwyuLqEb6rzP3+gNJ28txLdSsMr3ot5CZ/wAcVUEskiSpRcocon/ngNmX0Lgu3japJEkdrEwv+nkDHj8L3A+cVEk0kiSpJcqcg//AWAQiSZJaZ8Rz8BExNSKujoiHI+KXEfHtiJg6FsFJkqTRKdPJ7svANcC+wH7AtcUySRrW3LmN3t1lprlz2x2tVC9lEvw+mfnlzHy2mK4E9qk4Lkk10Mw13M2UlTSyMgn+VxHx/ojYqZjeD2yoOjBJ9ZE5/CSp9cok+FOBPwD+HVgPvKtYJkmSOtSICT4zH8jMEzNzn8z8rcx8e2b+YiyCk1rNc8KSdhRletEvjIg9B8zvFRFXlFhv/4joi4g1EXFHRJxVLF8QEQ9FxKpimjNgnY9GxD0RcVdEHDfaSklD8ZywpB1FmYFuDsrMx/pnMnNjRBxcYr1ngQ9n5q0RMQFYGRE3Fs9dlJmfGlg4Ig4A3gP8Ho0e+z+IiFdn5nOlaiI1wXG9JdVdmXPwL4mIvfpnImIi5QbIWZ+ZtxaPNwFraFxmN5STgEWZ+VRm3gfcgze1qa1mDpWfe+7r2h2uJHWdMgn+08BPIuJvIuKvgZ8A/6eZF4mIacDBwPJi0X+LiNsi4ooBPx72Ax4csNpahv9BoC7mLTC3n/0JJA2nTEv8KxGxgsYd5AJ4Z2b+rOwLRMR44NvA/8jMxyPii8Df0Lhhzd/Q+AFxarHtbV5+kO3NB+YDTJ48mWXLlpUNZUSbN29u6fagF6DF22zl65YtV4XGa/f1Df/as2c3yrUuxrLba/XrNrTyPXb99b1NlG1tXcrXoxfo7Pdi6z/3rdYLdPb/sAqdv1/KaWs9MrOyCdgFuAE4e4jnpwG3F48/Cnx0wHM3AEcOt/2enp5spb6+vpZur/8q37FW9nXbFV8zrz3yFdQvTnPmtP51W62V77F21qVsPdr9/y6j1Z/7VmvX/3DOnNZ+9prV6fulrCrqAazIEjl4yEP0EXHrSD8OhisTEQFcDqzJzAsHLJ8yoNg7gNuLx9cA74mI3SJiOjADuGWkGFRvc+aMXKafvd6l1vGKk+433CH610bEbcM8H8DLh3n+jcAfAasjYlWx7C+A90bETBqH3+8HzgDIzDsi4pvAz2j0wD8za9KDvkyP7DlzYPHi6mPpNosXNw459vb2Dluu/39s73eptbzipHsNl+B/t8T6QybgzLyJwc+rD/lbLzM/Dny8xOt2hTlzyv+y9Rfw9mnmf93MUQGpbubO7fzvm0aMvaXK2jga2pAJPh2tbruVfdP5C3j7+QGXymkmubfrx7CnB1qjzEA3qrmyPzD8pSzVRzfc5Kfs6YGy32HdUOdWKnMdvGqq2V/n/lJW1bymX83wdNvwTPA7sMWLy1+EJlXJqyU0Gs18h+2I32NDHqKPiE0MMtAMjY5zmZl7VBaVpI7VTAeosuyvIrXekC34zJyQmXsMMk0wuVfDw5PqBs22oDv5MOrcuY3REv3sqY5Kd7KLiN8CxvXPZ+YDlUS0A/JyOnWjOhzytLe26qzM/eBPjIifA/cBP6IxOM33Ko5rh1L2PJLGXp1u5lKXelTBz57qqEwnu78BjgDuzszpwNuAf6o0KqnN6nSFgR3YpB1TmQT/TGZuoHFf+JdkZh8ws+K4pLaq0xUGHiGSdkxlzsE/Vtzy9cfA30XEwzTGipdUQ/ZUl+qhTAv+JOA3wJ8C3wf+DZhXZVCSxl4zh/JnzdpQXSCSWmLEBJ+Z/5GZz2Xms5m5MDMvLg7ZS6rA3LnlOsS1WjOnJS64YHXrA2hCp3YYLLvvdvROjRobZXrRvzMifh4Rv46IxyNiU0Q8PhbBSTuibrgZSLt0eofBZl/TTo2qUplz8P8HmJeZa6oORtKL7Pi2rW4Z8a7Mvmv2RilSs8qcg/+lyV2SWquZoxE72pEatUaZFvyKiPgG8A/AU/0LM/M7lUUlSTXnrZdVtTIJfg/gCeDYAcsSMMFL2qF4OF3dZMQEn5kfGItAJKlTNXO/CA+nq1OMmOAj4uJBFv8aWJGZ3219SJLUGq1qcXs4Xd2oTCe7cTSGpv15MR0ETAROi4jPVBibulC7ruGWBnLQHqncOfhXAUdl5rMAEfFFYAlwDNDe0S7UcbyGW52gmRb3smWrgd6qQpHapkyC3w94GY3D8hSP983M5yLiqaFX047Ma7glqb3KDnSzKiKWAQG8BfjfEfEy4AcVxiZJkkapTC/6yyPieuBwGgn+LzJzXfH0OVUGJ0mSRmfITnYR8bvF30OAKcCDwAPAbxfLJEkdxg6u6jdcC/5sYD7w6UGeS+CoSiKSJDXNa/W1tSETfGbOL/7OHrtwJEmjsXgxLFu2jN7e3naHog5R5nax746ICcXjv4yI70TEwdWHJkmSRqvMQDd/lZmbIuJNwHHAQuBL1YYlSaqTsoNg2U+gdcok+OeKv3OBLxbD0+5aXUiSpLppZhAscITBVihzHfxDEXEJcDTwyYjYjXI/DCRJ2kLZQbAcYXD7lUnUfwDcAByfmY/RGIfe698lSepgZQa6eYIB937PzPXA+iqDkiRJ26fMIXpJJZTpHDRr1uu4+ebqY5Ekz6VL26mZQUOWL59UXSCSNIAt+C7kZSSdpeytSd1vksaSLfgu0kxL0UtMJI0lr23vPCb4LrJ4ceMSkzLTBResbne4knYAzTQ8HAN/bHmIXk3xl7ikgcqeotLYswWvUvyVLkndxRa8SvFXuiR1F1vwkiTVkAlekqQaMsFLklRDJnhJkmqosgQfEftHRF9ErImIOyLirGL5xIi4MSJ+Xvzdq1geEXFxRNwTEbdFxCFVxSZJUt1V2YJ/FvhwZr4WOAI4MyIOAM4FlmbmDGBpMQ/w+8CMYpoPfLHC2CRJqrXKEnxmrs/MW4vHm4A1wH7AScDCothC4O3F45OAr2TDzcCeETGlqvgkSaqzyMzqXyRiGvBj4EDggczcc8BzGzNzr4i4DrggM28qli8FPpKZK7ba1nwaLXwmT57cs2jRopbFuXnzZsaPH9+y7bWTdek8s2f3AtDXt6wl5dqpLvsErEunqktdqqjH7NmzV2bmoSMWzMxKJ2A8sBJ4ZzH/2FbPbyz+LgbeNGD5UqBnuG339PRkK/X19bV0e+1kXTpP/50CWlWuneqyTzKtS6eqS12qqAewIkvk30p70UfELsC3gb/LzO8Ui3/Zf+i9+PtwsXwtsP+A1acC66qMT5KkuqqyF30AlwNrMvPCAU9dA5xcPD4Z+O6A5X9c9KY/Avh1Zq6vKj5JkuqsyrHo3wj8EbA6IlYVy/4CuAD4ZkScBjwAvLt47npgDnAP8ATwgQpjkySp1ipL8NnoLDfUzUXfNkj5BM6sKh5JknYkjmQnSVINmeAlSaohE7wkSTVkgpckqYaq7EUvaRAxVNdTSWohW/DSGJkzp5qykjQYE7w0RhYvbowv/+JgtENPixe3O1pJ3c4EL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBqqLMFHxBUR8XBE3D5g2YKIeCgiVhXTnAHPfTQi7omIuyLiuKrikiRpR1BlC/5K4PhBll+UmTOL6XqAiDgAeA/we8U6fxsRO1UYmyRJtVZZgs/MHwOPlix+ErAoM5/KzPuAe4DDq4pNkqS6i8ysbuMR04DrMvPAYn4BcArwOLAC+HBmboyIzwM3Z+ZVRbnLge9l5rcG2eZ8YD7A5MmTexYtWtSyeDdv3sz48eNbtr12si6dqS51qUs9wLp0qrrUpYp6zJ49e2VmHjpiwcysbAKmAbcPmJ8M7ETjyMHHgSuK5V8A3j+g3OXAfx5p+z09PdlKfX19Ld1eO1mXzlSXutSlHpnWpVPVpS5V1ANYkSVy8Jj2os/MX2bmc5n5PHAZLx6GXwvsP6DoVGDdWMYmSVKdjGmCj4gpA2bfAfT3sL8GeE9E7BbWES7dAAAIiElEQVQR04EZwC1jGZskSXWyc1UbjoivA73A3hGxFjgP6I2ImUAC9wNnAGTmHRHxTeBnwLPAmZn5XFWxSZJUd5Ul+Mx87yCLLx+m/MdpnJeXJEnbyZHsJEmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBqqLMFHxBUR8XBE3D5g2cSIuDEifl783atYHhFxcUTcExG3RcQhVcUlSdKOoMoW/JXA8VstOxdYmpkzgKXFPMDvAzOKaT7wxQrjkiSp9ipL8Jn5Y+DRrRafBCwsHi8E3j5g+Vey4WZgz4iYUlVskiTV3Vifg5+cmesBir+/VSzfD3hwQLm1xTJJkjQKO7c7gEIMsiwHLRgxn8ZhfIDNEXFXC+PYG/hVC7fXTtalM9WlLnWpB1iXTlWXulRRj98pU2isE/wvI2JKZq4vDsE/XCxfC+w/oNxUYN1gG8jMS4FLqwguIlZk5qFVbHusWZfOVJe61KUeYF06VV3q0s56jPUh+muAk4vHJwPfHbD8j4ve9EcAv+4/lC9JkppXWQs+Ir4O9AJ7R8Ra4DzgAuCbEXEa8ADw7qL49cAc4B7gCeADVcUlSdKOoLIEn5nvHeKptw1SNoEzq4qlCZUc+m8T69KZ6lKXutQDrEunqktd2laPaORWSZJUJw5VK0lSDZngCxFxfETcVQyXe+7Ia3SuiLg/IlZHxKqIWNHueJrRzBDHnWyIeiyIiIeK/bIqIua0M8ayImL/iOiLiDURcUdEnFUs78b9MlRdumrfRMS4iLglIv61qMf5xfLpEbG82CffiIhd2x3rSIapy5URcd+AfTKz3bGWERE7RcS/RMR1xXzb9okJnsYOAb5AY8jcA4D3RsQB7Y1qu83OzJldeJnJlZQf4riTXcm29QC4qNgvMzPz+jGOabSeBT6cma8FjgDOLD4f3bhfhqoLdNe+eQo4KjNfD8wEji+uQPokjXrMADYCp7UxxrKGqgvAOQP2yar2hdiUs4A1A+bbtk9M8A2HA/dk5r2Z+TSwiMbwuRpjTQ5x3LGGqEdXysz1mXlr8XgTjS+v/ejO/TJUXbpKMaz35mJ2l2JK4CjgW8XybtknQ9Wl60TEVGAu8P+K+aCN+8QE31C3oXITWBIRK4uR/7rdUEMcd6P/Vtwx8YpuOKS9tYiYBhwMLKfL98tWdYEu2zfFoeBVNAYMuxH4N+CxzHy2KNI132Nb1yUz+/fJx4t9clFE7NbGEMv6DPDnwPPF/CTauE9M8A2lh8rtEm/MzENonHI4MyLe0u6ABDTukvhKGoch1wOfbm84zYmI8cC3gf+RmY+3O57tMUhdum7fZOZzmTmTxsifhwOvHazY2EY1OlvXJSIOBD4K/C5wGDAR+EgbQxxRRJwAPJyZKwcuHqTomO0TE3xD6aFyu0Fmriv+PgxcTePD381+2X93wa2GOO4qmfnL4ovseeAyumi/RMQuNBLi32Xmd4rFXblfBqtLN++bzHwMWEajT8GeEdE/vknXfY8NqMvxxemUzMyngC/T+fvkjcCJEXE/jdO8R9Fo0bdtn5jgG/4ZmFH0dtwVeA+N4XO7TkS8LCIm9D8GjgVuH36tjjfUEMddJba8BfI76JL9UpxHvBxYk5kXDniq6/bLUHXptn0TEftExJ7F492Bo2n0J+gD3lUU65Z9Mlhd7hzw4zFonLfu6H2SmR/NzKmZOY1GDvlhZr6PNu4TB7opFJfFfAbYCbgiMz/e5pBGJSJeQaPVDo2RCr/WTXWJAUMcA7+kMcTxPwDfBP4TxRDHmdnRHdiGqEcvjUPACdwPnNEN91yIiDcB/wis5sVzi39B49x1t+2XoeryXrpo30TEQTQ6bO1Eo6H2zcz86+Lzv4jGIe1/Ad5ftIA71jB1+SGwD43D3KuADw3ojNfRIqIX+LPMPKGd+8QEL0lSDXmIXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7ykykTEsojothseSbVggpckqYZM8NIOphjtcHFx/+3bI+K/RMTHIuKfi/lLi9HD+lvgF0XEj4t7qB8WEd8p7m39v4oy0yLizohYWNwY5FsR8dJBXvfYiPhpRNwaEX9fjAcvqSImeGnHczywLjNfn5kHAt8HPp+ZhxXzuwMnDCj/dGa+BfgSjWE2zwQOBE6JiElFmdcAl2bmQcDjwH8d+IIRsTfwl8DRxY2QVgBnV1ZDSSZ4aQe0Gjg6Ij4ZEW/OzF8DsyNieUSspnGTjN8bUP6aAevdUdwE5CngXl68SdODmflPxeOrgDdt9ZpHAAcA/1TcFvRk4HdaXjNJL9h55CKS6iQz746IHmAO8ImIWEKjVX5oZj4YEQuAcQNW6R83+/kBj/vn+79Dth7zeuv5oHGf7/e2oAqSSrAFL+1gImJf4InMvAr4FHBI8dSvivPi7xpy5aH9p4g4snj8XuCmrZ6/GXhjRLyqiOGlEfHqUbyOpJJswUs7ntcB/zcingeeAf6Exu04V9O4k9o/j2Kba4CTI+IS4OfAFwc+mZmPRMQpwNcjYrdi8V8Cd4+mApJG5t3kJG2XiJgGXFd00JPUITxEL0lSDdmClySphmzBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYb+P8hxaAd19wu9AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6102\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAF3CAYAAABJzllyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xu4XGV96PHvTy6Gmigm0DQQ2kSNVqoY2FzrbQfl0kBAPdqj1RYEjPbQ89CiVOzTaqjHo+1RUJQqoSCptEarUoGggHFvW6qGJjTlYrgVIoSkAiFIUuT+O3/M2riT7MuavWftmVn7+3meeTJrzbvW/N6s2fOb913veldkJpIkqV6e1+4AJElS65ngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBqqPMFHxC4R8e8RcVWxPDciVkXEnRHxtYjYvVj//GL5ruL1OVXHJklSXU1EC/4MYN2g5b8CzsvMecAW4NRi/anAlsx8GXBeUU6SJI1BpQk+ImYDxwF/WywHcCTwjaLIMuAtxfMTi2WK199UlJckSU2qugX/WeBPgWeL5RnAI5n5dLG8Adi3eL4vcB9A8frPi/KSJKlJu1a144g4HnggM9dERO/A6iGKZonXBu93MbAYYI899ujZb7/9WhBtw7PPPsvznlePcYfWpTPVpS51qQdYl05Vl7pUUY877rjjoczce9SCmVnJA/gkjRb6euC/gMeAvwceAnYtyhwBXFM8vwY4oni+a1EuRnqPnp6ebKW+vr6W7q+drEtnqktd6lKPTOvSqepSlyrqAazOEnm4sp9HmfmRzJydmXOAdwLfz8x3A33A24tiJwHfLp5fUSxTvP79oiKSJKlJ7ej/+DBwZkTcReMc+8XF+ouBGcX6M4Gz2xCbJEm1UNk5+MEysx/oL57fDRw6RJnHgXdMRDySJNXdhCR4SZKa9dRTT7FhwwYef/zxdocyZi960YtYt27d6AWHMGXKFGbPns1uu+02pu1N8JKkjrRhwwamTZvGnDlz6NZpUbZu3cq0adOa3i4z2bx5Mxs2bGDu3Lljeu/uvwZBklRLjz/+ODNmzOja5D4eEcGMGTPG1XthgpckdazJmNwHjLfuJnhJkoaxyy67MH/+fF71qlexaNEiHnnkkedeW7ZsGfPmzWPevHksW7ZshL20hwlekqRh7LHHHqxdu5ZbbrmF6dOnc8EFFwDw8MMPc84557Bq1SpuuOEGzjnnHLZs2dLmaLdngpckqYQjjjiC+++/H4BrrrmGo446iunTp/PiF7+Yo446iu9+97ttjnB7JnhJUseLqOZR1jPPPMPKlSs54YQTALj//vsZfC+U2bNnP5f8O4UJXpKkYfziF79g/vz5zJgxg4cffpijjjoKgKFmUu+0AYEmeElSx8us5jGagXPwP/3pT3nyySefOwc/e/Zs7rvvvufKbdiwgX322aeq6o+JCV6SpFG86EUv4vzzz+fTn/40Tz31FMcccwzXXnstW7ZsYcuWLVx77bUcc8wx7Q5zOyZ4SZJKOPDAA3nNa17D8uXLmT59On/xF3/BIYccwiGHHMJHP/pRpk+f3u4Qt+NUtZIkDWPbtm3bLV955ZXPPT/llFM45ZRTJjqk0mzBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIkDWOk28Uee+yx7Lnnnhx//PFtjHB4JnhJkoYx3O1iAc466yy+8pWvtDG6kZngJUmdrapbyTV5S7nBt4sFeNOb3sS0adOqqHFLmOAlSRrFjreL7QYmeEmShjHc7WK7gQlekqRhDHe72G5ggpckdbaqbgZf9qbw7Hy72G5ggpckqYTBt4sFeP3rX8873vEOVq5cyezZs7nmmmvaHOH2vF2sJEnDGOl2sf/yL/8y0eE0xRa8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSk84//3xe+cpX8u53v5snnniCN7/5zcyfP5+vfe1rTe2nv7+fH/7wh5XE6HXwkiQ16W/+5m/4zne+w9y5c/nxj3/MU089xdq1a5veT39/P1OnTuW3f/u3Wx6jCV6SpBGce+65XHLJJQCcdtpp3Hbbbdx9992ccMIJvOc97+Giiy7iwQcfZP78+Xzzm9/koosu4oorrmDXXXelt7eX888/nwcffJAPfOAD3HvvvQB89rOfZd999+VLX/oSu+yyC5dddhmf//znef3rX9+yuE3wkqSusKR/Cef84JxSZd930PtYumjpdusWX7mYi2686Lnlj73xYyzpXTLiftasWcOXv/xlVq1aRWZy2GGHcdlll/Hd736Xvr4+9tprLw477DA+/elPc9VVV/Hwww9z+eWXc9tttxER3HfffQCcccYZ/Mmf/Amve93ruPfeeznmmGNYt24dH/jAB5g6dSof+tCHmvvPKMEEL0nSMK6//nre+ta38oIXvACAt73tbSNOUfvCF76QKVOmcNppp3Hcccfxxje+EYDvfe97/OQnP3mu3KOPPsrWrVsrjd0EL0nSMLLk3eYG7Lrrrtxwww2sXLmS5cuX87nPfY4f/OAHPPvss/zoRz9ijz32qCjSIWKZsHeSJGkclvQuGbVLfSRLFy3dqdt+NG94wxs4+eSTOfvss8lMLr/8cr7yla9w7rnnDll+27ZtPPbYYyxcuJDDDz+cl770pQAcffTRfOELX+Css84CYO3atcyfP59p06bx6KOPjrlOI/EyOUmShnHQQQdx8sknc+ihh3LYYYdx2mmnceCBBw5bfuvWrRx//PEccMABvPGNb+STn/wk0LisbvXq1RxwwAHsv//+fOlLXwJg0aJFXH755cyfP7/ld6ezBS9J0gjOPPNMzjzzzO3WrV+//rnnvb299Pb2AjBr1ixuuOGG514bOM++1157DXmN/Mtf/nJuuumm1geNLXhJkmrJBC9JUg2Z4CVJqqHKEnxETImIGyLiPyLi1og4p1h/aUTcExFri8f8Yn1ExPkRcVdE3BQRB1UVmySpOzR7mVqdjLfuVQ6yewI4MjO3RcRuwPUR8Z3itbMy8xs7lP8dYF7xOAz4YvGvJGkSmjJlCps3b2bGjBlERLvDmVCZyebNm5kyZcqY91FZgs/GT49txeJuxWOknyMnAn9XbPfjiNgzImZl5qaqYpQkda7Zs2ezYcMGHnzwwXaHMmaPP/74mJP0lClTmD179pjfO6rs/oiIXYA1wMuACzLzwxFxKXAEjRb+SuDszHwiIq4CPpWZ1xfbrgQ+nJmrd9jnYmAxwMyZM3uWL1/esni3bdvG1KlTW7a/drIunakudalLPcC6dKq61KWKeixYsGBNZh48asHMrPwB7An0Aa8CZgEBPB9YBny0KLMCeN2gbVYCPSPtt6enJ1upr6+vpftrJ+vSmepSl7rUI9O6dKq61KWKegCrs0TunZBR9Jn5CNAPHJuZm4oYnwC+DBxaFNsA7Ddos9nAxomIT5KkuqlyFP3eEbFn8XwP4M3AbRExq1gXwFuAW4pNrgD+oBhNfzjw8/T8uyRJY1LlKPpZwLLiPPzzgK9n5lUR8f2I2JtGN/1a4ANF+auBhcBdwGPAeyuMTZKkWqtyFP1NwE4z8mfmkcOUT+D0quKRJGkycSY7SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeDb7brr4CUvgYhyj5e8pLGNJEkjMMG32/vfD/fcU778Pfc0tpEkaQQm+HZrJrmPZxtJ0qRigtfk4ekQSZOICb4qZZPJjjKHfuxolP0e9nu/Z3LakadDJE0iJviqNJtMWmyPTZtMTjvydIikScQEX5WxJIa5c8f22lhiaLa72i5rVcVTJ1IlTPCdYu5cuPDC4V+/8MKxJfnhviSPPrr5HyF167IuezpE1fLUiVSJXdsdwKQx3sRx1FFw992jlxvqvH4rdWKX9XXXtf2UiMbBUyf1N4a/0d6qYplEP+JtwdfNWFr53a7q5N7N3cat7v4u9te7YEFzp3eaGWiq+vEHeFuY4Oum2a78uXPh2muH767uhl+7rR7vMNYYOrHbuNXd3xPxRd2tn0MNz+TeFpUl+IiYEhE3RMR/RMStEXFOsX5uRKyKiDsj4msRsXux/vnF8l3F63Oqiq3Wiq78/r6+kZP2wOPuuxvbNKNVLbdWtRabNdp4h7Em/078Emt193fVdWzm/75Fn8NaXFLabM9Kp/Y4qaWqbME/ARyZma8B5gPHRsThwF8B52XmPGALcGpR/lRgS2a+DDivKNe9bH2M7p57GoP9WjkQEMb/o2asAxo1PqP98KrIqJeUdsMo/24bqFii8VG6kdLsYxKJnIAKR8SvANcDfwisAH4tM5+OiCOAJZl5TERcUzz/UUTsCvwXsHeOEODBBx+cq1evblmcJ196Mst+uqxU2fcd9D6WLlq63brFVy7mohsvKrX9x974MZb0Ltlu3aKvLuKqO64qtf2Fx1/I4p7F263rWdrDjZtuLLX9Fe+8gkWvWLTdun0+sw+btm0qtf3qC6Fnh6KxpNSmANz/Gdhn6y+XN06DfT9Yfvvc4b3WzIKDS35fzZo6i40f3Ljduitvv5ITlp9QavuDZh3EmsVrGgtFD8LSHnj/ohE2GuT4lx/Ple+6crt1S/qXcM4Pzim1/fvWwNKBzYukuPjxfxzfZ+/3gqteUWpzLrwSBqo/8IXpZ6/cts999gYNOrvy5XDC75XbfrvPXmHpmqW8/6pyARx/O1z51UErMpv77I33e68flvSzXaLtqu+9962mZ5+e7dbFOSP3Ip70Gydx6cmXltp/WRGxJjMPHq1cpaPoI2IXYA3wMuAC4D+BRzLz6aLIBmDf4vm+wH0ARfL/OTADeGiHfS4GFgPMnDmT/v7+lsX75JNPli67cdPGnd5746aNQxcewvr163fafvNDm0tvf/vtt9O/dfvtt27dOnThIdx8881M2zRtu3XN1L+bPfnkkzv939/80M2lt9+6detz2/eO4f03P7R5p/dfv379GPYE3HMPvzjpJDYuOaj0JkN99sZqYD9+9soZ+OwddtJJjZ6DJg3+7A24fePtY46nv7+/qc/eeL/3Br/vgG763lu9ZjVb7yj/fgP7b2WeakalCT4znwHmR8SewOXAK4cqVvw71M+gnVrvmbkUWAqNFnxvb29rggUuXX9p6bL7zNqHHd/7H7b+A5T8m50zZ85O28/YNAMeLrf9K17xCnp7tt9+2h3TYFu57V/96lfT+4rtt999ze5Q9rO+ejXs8EuWUX7Jbme//eAn95Uvv6MdO3Y2roGLRv1BC8Duu+++0//91tu3wq3l3nratGk7bd+MGfnsTtv39/fDT8e2vz02bWKfq1ZAz+hlYejP3meWDl12NAP76arP3v33w7R92jKC/7nP3hiSO8C0O+6k90MLtlt3Rw9QsvdoR729vU199sb7vTf4fQd00/fewT0H79SC5wcjbzPU982EycwJeQAfA86i0SLftVh3BHBN8fwa4Iji+a5FuRhpnz09PdlKfX19Ld1fO1mXCTKWs4Bz57Z2f808Wv2+3Wys/zfjLTdU2U49znPnZl577ch1KWMMn5uO/rtvQhX1AFZnibxb5Sj6vYuWOxGxB/BmYB3QB7y9KHYS8O3i+RXFMsXr3y8qInWuVk8hvKPhvnrH8r6tngq5bgMR23md/kQd52a1ezCexqXKUfSzgL6IuAn4N+C6zLwK+DBwZkTcReMc+8VF+YuBGcX6M4GzK4xNao12jbgfy3wHLZwK+RezZrVlxPuk0+rj0q57WqgtKjsHn5k3AQcOsf5u4NAh1j8OvKOqeKRKtGsK4bLvW9H+VvX3t++8YqvMndv8ZZgT/WOu1cflwgtbO1mRM9R1NGeykzQ5tboXpBsM/GAY7ex7WZ0wi6SG5c1mJE1OA7M+VtUbYdf0zurwI6mLmOCldvDLX3XkuOiOYhe9JLXCWLue7bJWRUzw0kRo9+VOqt5Yrqiwy1oVsotemgjNjl72i7/7tPrKBmmcTPDSRKh6QJck7cAuekmSasgEL0lSDdlFL0kamZd1diVb8JIk1ZAJXpK0PS/rrAUTvCRpe5Nxnv4aGvYcfEScWWL7/85Mj6ok1YnX9NfCSC34s4CpwLQRHh+sOkBJktS8kUbRfyUz/3KkjSPiBS2OR5IktcCwLfjM/NPRNi5TRpIkTbxRB9lFxBkR8cJouDgiboyIoyciOEmSNDZlRtGfkpmPAkcDewPvBT5VaVSSJGlcyiT4gSmMFgJfzsz/GLROkiR1oDIJfk1EXEsjwV8TEdOAZ6sNS5IkjUeZuehPBeYDd2fmYxExg0Y3vSRJ6lBlWvDXZeaNmfkIQGZuBs6rNixJkjQeI81kNwX4FWCviHgxvzzv/kJgnwmITZIkjdFIXfTvB/6YRjJfwy8T/KPABRXHJUmSxmHYBJ+ZnwM+FxH/OzM/P4ExSZKkcRp1kF1mfj4ifhuYM7h8Zv5dhXFJkqRxGDXBR8RXgJcCa4FnitUJmOAlSepQZS6TOxjYPzOz6mAkSVJrlLlM7hbg16oORJIktU6ZFvxewE8i4gbgiYGVmXlCZVFJkqRxKZPgl1QdhCRJaq0yo+h/MBGBSJKk1ikzin4rjVHzALsDuwH/nZkvrDIwSZI0dmVa8NMGL0fEW4BDK4tIkiSNW5lR9NvJzH8CjqwgFkmS1CJluujfNmjxeTSui/eaeEmSOliZUfSLBj1/GlgPnFhJNJIkqSXKnIN/70QEIkmSWmfUc/ARMTsiLo+IByLiZxHxzYiYPRHBSZKksSkzyO7LwBU07gu/L3BlsU6SJHWoMgl+78z8cmY+XTwuBfauOC5JkjQOZRL8QxHxnojYpXi8B9hcdWCSJGnsyiT4U4DfBf4L2AS8vVgnSZI6VJlR9PcC3jlOkqQuUmYU/bKI2HPQ8osj4pIS2+0XEX0RsS4ibo2IM4r1SyLi/ohYWzwWDtrmIxFxV0TcHhHHjLVSkiRNdmUmujkgMx8ZWMjMLRFxYIntngY+mJk3RsQ0YE1EXFe8dl5mfnpw4YjYH3gn8Fs0Rux/LyJenpnPlKqJJEl6Tplz8M+LiBcPLETEdMp17W/KzBuL51uBdTQusxvOicDyzHwiM+8B7sKb2kiSNCZlEvxngB9GxMcj4i+BHwJ/3cybRMQc4EBgVbHqjyLipoi4ZNCPh32B+wZttoGRfxBIkqRhRObo940pus+PBAJYmZk/Kf0GEVOBHwCfyMxvRcRM4CEaN6z5ODArM0+JiAuAH2XmZcV2FwNXZ+Y3d9jfYmAxwMyZM3uWL19eNpRRbdu2jalTp7Zsf+1kXTpTXepSl3qAdelUdalLFfVYsGDBmsw8eNSCmVnZA9gNuAY4c5jX5wC3FM8/Anxk0GvXAEeMtP+enp5spb6+vpbur52sS2eqS13qUo9M69Kp6lKXKuoBrM4SOXjYLvqIuHG0HwcjlYmIAC4G1mXmuYPWzxpU7K3ALcXzK4B3RsTzI2IuMA+4YbQYJEnSzkYaLPfKiLhphNcDeNEIr78W+H3g5ohYW6z7M+BdETGfRhf9euD9AJl5a0R8HfgJjRH4p6cj6CVJGpOREvxvlth+2AScmdfT+BGwo6tH2OYTwCdKvK8kSRrBsAk+M386kYFIkqTWKXOZnCRJ6jImeEmSasgEL0lSDQ17Dj4ittIY6b7TS0Bm5gsri0qSJI3LSIPspk1kIJIkqXXK3E0OgIj4VWDKwHI27hMvSZI6UJn7wZ8QEXcC99CYU3498J2K45IkSeNQZpDdx4HDgTsycy7wJuBfK41KkiSNS5kE/1RmbqZxX/jnZWYfML/iuCRJ0jiUOQf/SHHL138G/j4iHqAxV7wkSepQZVrwJwK/AP4E+C7wn8CiKoOSNPGOOw4iyj3OPvvV7Q5X0ihGTfCZ+d+Z+UxmPp2ZyzLz/KLLXlKNXD3sbaB2tmrVjOoC6WLN/EiKaJSXqlJmFP3bIuLOiPh5RDwaEVsj4tGJCE7SxMsc+dGMZhJeHZJdMz+SxlJeakaZc/B/DSzKzHVVByOpXppJYHVKdmV+CMVQN9OWWqhMgv+ZyV3SeIyW8Ex2UuuVSfCrI+JrwD8BTwyszMxvVRaVJEkalzIJ/oXAY8DRg9YlYIKXJKlDjZrgM/O9ExGIJElqnTKj6M8f4vHxiDhxIgKUpAGTbVR+N/CYdK4yE91MoTE17Z3F4wBgOnBqRHy2wtgmDf9ApHIm66j8TuYx6VxlzsG/DDgyM58GiIgvAtcCRwE3VxjbpOEfiNQcR+V3Ho9J5ynTgt8XeMGg5RcA+2TmMwwaVa/xa+UEI5Kkya3sRDdrI6IfCOANwP+NiBcA36swNkmSNEZlRtFfHBFXA4fSSPB/lpkbi5fPqjI4SZI0NsN20UfEbxb/HgTMAu4D7gV+rVgnqQnHHQcLFvQ6mLKDeExUZyO14M8EFgOfGeK1BI6sJCKpphxM2Xk8JqqzYVvwmbm4+HfBEA+TuzpKN11q6GDKzjOZjkk3/a1ofMpMdPOOiJhWPP/ziPhWRBxYfWhSebbEpHL8W5k8ylwm9xeZuTUiXgccAywDvlRtWNLYTKaWmDQe/q3UX5kE/0zx73HAFzPz28Du1YVUHwPdXJIkTbQyCf7+iLgQ+F3g6oh4fsntJElSm5RJ1L8LXAMcm5mP0JiH3uvfJUnqYGUmunmMQfd+z8xNwKYqg5IkSeNjV7skSTVkgpckqYZM8JIk1ZAJXupSzkgmaSQmeKlLOSOZpJGUuR+8pA422qxjTrYkTU624KUhNNP9bRe4pE5kgpeG0GyXtl3gkjqNXfTSCMrcdMMucEmdyBa8SnHEtiR1F1vwKqVuI7ZtdUuqO1vwakq330N64cJqykpSp6kswUfEfhHRFxHrIuLWiDijWD89Iq6LiDuLf19crI+IOD8i7oqImyLioKpi0+S1YsXoP1IGHitWtDtaSRq7KlvwTwMfzMxXAocDp0fE/sDZwMrMnAesLJYBfgeYVzwWA1+sMDZJkmqtsgSfmZsy88bi+VZgHbAvcCKwrCi2DHhL8fxE4O+y4cfAnhExq6r4JEmqswkZZBcRc4ADgVXAzOKe8mTmpoj41aLYvsB9gzbbUKzb7t7zEbGYRgufmTNn0t/f37I4t23b1tL9QS9AiX22ulw31KW8cnVp/fuWV/a9J1u5dr53p5drtmw5rf1baVe5htZ/h7VHW+uRmZU+gKnAGuBtxfIjO7y+pfh3BfC6QetXAj0j7bunpydbqa+vr6X7GzibO9HlMju/Ls0oU5cq3resdh3nTi/Xzvfu9HLNli2rlX8r7fy/yWz9d1i7VFEPYHWWyL+VjqKPiN2AbwJ/n5nfKlb/bKDrvfj3gWL9BmC/QZvPBjZWGZ8kdbpm5qA4++xXtztcdZAqR9EHcDGwLjPPHfTSFcBJxfOTgG8PWv8HxWj6w4GfZ9GVL0mTVTPzSqxaNaO6QNR1qjwH/1rg94GbI2Jtse7PgE8BX4+IU4F7gXcUr10NLATuAh4D3lthbJLUVbxroJpVWYLPzOuB4T5ybxqifAKnVxWP1G38wpY0Hs5kJ3UYZ9uT1ArORS91GGfQk9QKtuDVFo4MlqRqmeDVFo4MlqRq2UWvtnJksCRVwxa8JEk1ZAteUm3ZA6TJzBa8pNpp5vLBww7bXF0gGpPjjoMFC3pLDcI97rh2R9u5TPA15B+HJrsVKwZubTL641Ofurnd4WoHzQzCbabsZGMXfQ35xyGpDhyEOz4m+Brzj0OSJi+76CVJqiETvCRJNWQXvaQx8RSP1NlswUtqSrN3sPOOd1J7mOAlNWXFCujr6y99GZp3x5PawwQvSVINmeAlSaohE7wkSTXkKHqpRRxVLqmT2IKXxskbm6jORrufhTqXCV4aJ29sojpq5oerl0J2JrvoJUk78fLG7mcLXpKkGrIFL0klVXHO2fPYqooteEkdo1MHdFVxPtpz3KqaLXhJbbdwIVx9dfmyE62K89Ge41bVTPCS2s5kJ7WeXfSSJNWQLXhVwoFDktRetuDVUg4ckqTOYAteLVXVuVR7BCSpObbg1dHsEZCksTHBq6OtWAF9ff2l5nl3JLY0OY02f0K751FoFxO8JKkr2Ws3Ms/BS5K6kr12I7MFL0lSDdmC70KT8VySJKk5tuC7SDPnmw47bHN1gUiSOp4t+C7SzPmm/v6bgd6qQpEkdTgTvDRJeGpHmlzsopdqzsmC1Em8Xn3i2IKXas5LidQJFi6Eq68uX74xjmhGZfFMBib4CeAvUkmTXbM/NB1HNH520VfIrlFJUrtUluAj4pKIeCAibhm0bklE3B8Ra4vHwkGvfSQi7oqI2yPimKrimkgrVow+f7rzqEuSqlBlC/5S4Ngh1p+XmfOLx9UAEbE/8E7gt4pt/iYidqkwNkmqJQewaUBlCT4z/xl4uGTxE4HlmflEZt4D3AUcWlVsklQ3nhLUjiIzq9t5xBzgqsx8VbG8BDgZeBRYDXwwM7dExBeAH2fmZUW5i4HvZOY3htjnYmAxwMyZM3uWL1/esni3bdvG1KlTW7a/dipTlwULeoHG7VhH00zZVptsx6Ub1KUeYF06VV3qUkU9FixYsCYzDx61YGZW9gDmALcMWp4J7EKj5+ATwCXF+guA9wwqdzHwP0bbf09PT7ZSX19fS/fXTmXqMjACoIxmyrbaZDsu3aAu9ci0Lp2qLnWpoh7A6iyRgyd0FH1m/iwzn8nMZ4GL+GU3/AZgv0FFZwMbJzI2SZLqZEITfETMGrT4VmBghP0VwDsj4vkRMReYB9wwkbFJklQnlU10ExFfpTFLwV4RsQH4GNAbEfOBBNYD7wfIzFsj4uvAT4CngdMz85mqYtP2HFkrSfVTWYLPzHcNsfriEcp/gsZ5eU2QZqeOdOStJHUPp6qdxJxcR5Lqy6kSpvS0AAAIMUlEQVRqJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNVQZQk+Ii6JiAci4pZB66ZHxHURcWfx74uL9RER50fEXRFxU0QcVFVckiRNBlW24C8Fjt1h3dnAysycB6wslgF+B5hXPBYDX6wwLkmSaq+yBJ+Z/ww8vMPqE4FlxfNlwFsGrf+7bPgxsGdEzKoqNkmS6m6iz8HPzMxNAMW/v1qs3xe4b1C5DcU6SZI0Bru2O4BCDLEuhywYsZhGNz7Atoi4vYVx7AU81ML9tZN16Ux1qUtd6gHWpVPVpS5V1OM3yhSa6AT/s4iYlZmbii74B4r1G4D9BpWbDWwcageZuRRYWkVwEbE6Mw+uYt8Tzbp0prrUpS71AOvSqepSl3bWY6K76K8ATiqenwR8e9D6PyhG0x8O/HygK1+SJDWvshZ8RHwV6AX2iogNwMeATwFfj4hTgXuBdxTFrwYWAncBjwHvrSouSZImg8oSfGa+a5iX3jRE2QROryqWJlTS9d8m1qUz1aUudakHWJdOVZe6tK0e0citkiSpTpyqVpKkGjLBFyLi2Ii4vZgu9+zRt+hcEbE+Im6OiLURsbrd8TSjmSmOO9kw9VgSEfcXx2VtRCxsZ4xlRcR+EdEXEesi4taIOKNY343HZbi6dNWxiYgpEXFDRPxHUY9zivVzI2JVcUy+FhG7tzvW0YxQl0sj4p5Bx2R+u2MtIyJ2iYh/j4iriuW2HRMTPI0DAlxAY8rc/YF3RcT+7Y1q3BZk5vwuvMzkUspPcdzJLmXnegCcVxyX+Zl59QTHNFZPAx/MzFcChwOnF38f3XhchqsLdNexeQI4MjNfA8wHji2uQPorGvWYB2wBTm1jjGUNVxeAswYdk7XtC7EpZwDrBi237ZiY4BsOBe7KzLsz80lgOY3pczXBmpziuGMNU4+ulJmbMvPG4vlWGl9e+9Kdx2W4unSVYlrvbcXibsUjgSOBbxTru+WYDFeXrhMRs4HjgL8tloM2HhMTfEPdpspN4NqIWFPM/NfthpviuBv9UXHHxEu6oUt7RxExBzgQWEWXH5cd6gJddmyKruC1NCYMuw74T+CRzHy6KNI132M71iUzB47JJ4pjcl5EPL+NIZb1WeBPgWeL5Rm08ZiY4BtKT5XbJV6bmQfROOVwekS8od0BCWjcJfGlNLohNwGfaW84zYmIqcA3gT/OzEfbHc94DFGXrjs2mflMZs6nMfPnocArhyo2sVGNzY51iYhXAR8BfhM4BJgOfLiNIY4qIo4HHsjMNYNXD1F0wo6JCb6h9FS53SAzNxb/PgBcTuOPv5v9bODugjtMcdxVMvNnxRfZs8BFdNFxiYjdaCTEv8/MbxWru/K4DFWXbj42mfkI0E9jTMGeETEwv0nXfY8NqsuxxemUzMwngC/T+cfktcAJEbGexmneI2m06Nt2TEzwDf8GzCtGO+4OvJPG9LldJyJeEBHTBp4DRwO3jLxVxxtuiuOuEtvfAvmtdMlxKc4jXgysy8xzB73UdcdluLp027GJiL0jYs/i+R7Am2mMJ+gD3l4U65ZjMlRdbhv04zFonLfu6GOSmR/JzNmZOYdGDvl+Zr6bNh4TJ7opFJfFfBbYBbgkMz/R5pDGJCJeQqPVDo2ZCv+hm+oSg6Y4Bn5GY4rjfwK+Dvw6xRTHmdnRA9iGqUcvjS7gBNYD7++Gey5ExOuAfwFu5pfnFv+Mxrnrbjsuw9XlXXTRsYmIA2gM2NqFRkPt65n5l8Xf/3IaXdr/DrynaAF3rBHq8n1gbxrd3GuBDwwajNfRIqIX+FBmHt/OY2KClySphuyilySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC+pMhHRHxHddsMjqRZM8JIk1ZAJXppkitkOVxT3374lIv5nRHw0Iv6tWF5azB420AI/LyL+ubiH+iER8a3i3tb/pygzJyJui4hlxY1BvhERvzLE+x4dET+KiBsj4h+L+eAlVcQEL00+xwIbM/M1mfkq4LvAFzLzkGJ5D+D4QeWfzMw3AF+iMc3m6cCrgJMjYkZR5hXA0sw8AHgU+F+D3zAi9gL+HHhzcSOk1cCZldVQkglemoRuBt4cEX8VEa/PzJ8DCyJiVUTcTOMmGb81qPwVg7a7tbgJyBPA3fzyJk33Zea/Fs8vA163w3seDuwP/GtxW9CTgN9oec0kPWfX0YtIqpPMvCMieoCFwCcj4loarfKDM/O+iFgCTBm0ycC82c8Oej6wPPAdsuOc1zsuB437fL+rBVWQVIIteGmSiYh9gMcy8zLg08BBxUsPFefF3z7sxsP79Yg4onj+LuD6HV7/MfDaiHhZEcOvRMTLx/A+kkqyBS9NPq8G/l9EPAs8Bfwhjdtx3kzjTmr/NoZ9rgNOiogLgTuBLw5+MTMfjIiTga9GxPOL1X8O3DGWCkganXeTkzQuETEHuKoYoCepQ9hFL0lSDdmClySphmzBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYb+P9E0Psvf+5hYAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6103\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAF3CAYAAABJzllyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xu4HXV56PHvK7dQE4QEmiaE00SNVooYSLh53UFuDQTUoz1aPQW5RFp6HlqUij2thPb4SHsQFLVCKJco1ehRqdwUNO6NpUgwoSkXw60QISSVW5CkyP09f6wJ3Qn7MmtnzV5rzf5+nmc9e82s38y8vz2z5l2/md/MRGYiSZLq5VXtDkCSJLWeCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqocoTfERsExH/GhFXF8MzImJZRNwbEd+MiO2L8TsUw/cVn0+vOjZJkupqNFrwpwKr+g3/LXBeZs4E1gMnFONPANZn5uuB84pykiRpBCpN8BExDTgS+IdiOICDgW8XRRYD7yneH1MMU3z+7qK8JElqUtUt+M8Dfw68VAxPAp7MzBeK4TXA7sX73YGHAIrPf1WUlyRJTdq2qhlHxFHAI5m5IiJ6No0eoGiW+Kz/fBcACwB23HHH2XvssUcLom146aWXeNWr6tHv0Lp0prrUpS71AOvSqepSlyrqcc899zyWmbsNWzAzK3kBn6XRQl8N/AfwNPCPwGPAtkWZg4DrivfXAQcV77ctysVQy5g9e3a2Um9vb0vn107WpTPVpS51qUemdelUdalLFfUAlmeJPFzZz6PM/FRmTsvM6cAHgR9n5oeBXuD9RbFjge8V768shik+/3FREUmS1KR2HP/4JHBaRNxH4xz7xcX4i4FJxfjTgDPaEJskSbVQ2Tn4/jKzD+gr3t8P7D9AmWeAD4xGPJIk1d2oJHhJkpr1/PPPs2bNGp555pl2hzJir3nNa1i1atXwBQcwbtw4pk2bxnbbbTei6U3wkqSOtGbNGiZMmMD06dPp1tuibNiwgQkTJjQ9XWby+OOPs2bNGmbMmDGiZXf/NQiSpFp65plnmDRpUtcm960REUyaNGmrjl6Y4CVJHWssJvdNtrbuJnhJkgaxzTbbMGvWLPbaay/mz5/Pk08++fJnixcvZubMmcycOZPFixcPMZf2MMFLkjSIHXfckZUrV3LHHXcwceJEvvzlLwPwxBNPcNZZZ7Fs2TJuueUWzjrrLNavX9/maDdngpckqYSDDjqIhx9+GIDrrruOQw89lIkTJ7LLLrtw6KGH8oMf/KDNEW7OBC9J6ngR1bzKevHFF1m6dClHH300AA8//DD9n4Uybdq0l5N/pzDBS5I0iF//+tfMmjWLSZMm8cQTT3DooYcCMNCd1DutQ6AJXpLU8TKreQ1n0zn4X/ziFzz33HMvn4OfNm0aDz300Mvl1qxZw9SpU6uq/oiY4CVJGsZrXvMazj//fM455xyef/55Dj/8cK6//nrWr1/P+vXruf766zn88MPbHeZmTPCSJJWwzz778Ja3vIUlS5YwceJE/uqv/or99tuP/fbbj09/+tNMnDix3SFuxlvVSpI0iI0bN242fNVVV738/vjjj+f4448f7ZBKswUvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5I0iKEeF3vEEUew8847c9RRR7UxwsGZ4CVJGsRgj4sFOP300/na177WxuiGZoKXJHW2qh4l1+Qj5fo/Lhbg3e9+NxMmTKiixi1hgpckaRhbPi62G5jgJUkaxGCPi+0GJnhJkgYx2ONiu4EJXpLU2ap6GHzZh8LzysfFdgMTvCRJJfR/XCzAO97xDj7wgQ+wdOlSpk2bxnXXXdfmCDfn42IlSRrEUI+L/ed//ufRDqcptuAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZKadP755/OmN72JD3/4wzz77LMccsghzJo1i29+85tNzaevr4+bbrqpkhi9Dl6SpCb9/d//Pd///veZMWMGN998M88//zwrV65sej59fX2MHz+et771rS2P0QQvSdIQzj33XC655BIATjzxRO666y7uv/9+jj76aD7ykY9w0UUX8eijjzJr1iy+853vcNFFF3HllVey7bbb0tPTw/nnn8+jjz7KySefzIMPPgjA5z//eXbffXcuuOACttlmGy6//HK++MUv8o53vKNlcZvgJUldYWHfQs664axSZU/a9yQWzV+02bgFVy3golsvenn4zHedycKehUPOZ8WKFVx66aUsW7aMzOSAAw7g8ssv5wc/+AG9vb3suuuuHHDAAZxzzjlcffXVPPHEE1xxxRXcddddRAQPPfQQAKeeeip/9md/xtvf/nYefPBBDj/8cFatWsXJJ5/M+PHj+cQnPtHcP6MEE7wkSYO48cYbee9738urX/1qAN73vvcNeYvanXbaiXHjxnHiiSdy5JFH8q53vQuAH/3oR/z85z9/udxTTz3Fhg0bKo3dBC9J0iCy5NPmNtl222255ZZbWLp0KUuWLOELX/gCN9xwAy+99BI//elP2XHHHSuKdIBYRm1JkiRthYU9C4c9pD6URfMXveKw/XDe+c53ctxxx3HGGWeQmVxxxRV87Wtf49xzzx2w/MaNG3n66aeZN28eBx54IK973esAOOyww/jSl77E6aefDsDKlSuZNWsWEyZM4KmnnhpxnYbiZXKSJA1i33335bjjjmP//ffngAMO4MQTT2SfffYZtPyGDRs46qij2HvvvXnXu97FZz/7WaBxWd3y5cvZe++92XPPPbngggsAmD9/PldccQWzZs1q+dPpbMFLkjSE0047jdNOO22zcatXr375fU9PDz09PQBMmTKFW2655eXPNp1n33XXXQe8Rv4Nb3gDt912W+uDxha8JEm1ZIKXJKmGTPCSJNVQZQk+IsZFxC0R8W8RcWdEnFWMvywiHoiIlcVrVjE+IuL8iLgvIm6LiH2rik2S1B2avUytTra27lV2snsWODgzN0bEdsCNEfH94rPTM/PbW5T/PWBm8ToA+ErxV5I0Bo0bN47HH3+cSZMmERHtDmdUZSaPP/4448aNG/E8Kkvw2fjpsbEY3K54DfVz5Bjgq8V0N0fEzhExJTPXVRWjJKlzTZs2jTVr1vDoo4+2O5QRe+aZZ0acpMeNG8e0adNGvOyo8vBHRGwDrABeD3w5Mz8ZEZcBB9Fo4S8FzsjMZyPiauDszLyxmHYp8MnMXL7FPBcACwAmT548e8mSJS2Ld+PGjYwfP75l82sn69KZ6lKXutQDrEunqktdqqjH3LlzV2TmnGELZmblL2BnoBfYC5gCBLADsBj4dFHmGuDt/aZZCswear6zZ8/OVurt7W3p/NrJunSmutSlLvXItC6dqi51qaIewPIskXtHpRd9Zj4J9AFHZOa6IsZngUuB/Ytia4A9+k02DVg7GvFJklQ3Vfai3y0idi7e7wgcAtwVEVOKcQG8B7ijmORK4A+L3vQHAr9Kz79LkjQiVfainwIsLs7Dvwr4VmZeHRE/jojdaBymXwmcXJS/FpgH3Ac8DXy0wtgkSaq1KnvR3wa84o78mXnwIOUTOKWqeCRJGku8k50kSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyATfbj/8Ibz2tRBR7vXa1zamkTqd27Y2aXZbiOCAP/gDt4etZIJvt499DB54oHz5Bx5oTCN1Orft+iubuA87rLltAdhx3Tq3h61kgm+3Jjf6EU9TZ2OxpdgNdXbb7l4VJu6muD1sFRO8hteuZFIst2fu3NbuZOrQUrR1rCo1u32pI5ngq1I2KW4pc+BXO40kmRx22Nb/EKhyJ9PtO69ubB2X3bY78WhEXZTdL1W1rcyYAddfP/i24PbQUib4qnT6L+BmWuVVJtmhWpWd/P9rVtmjEe60hjfcdtPkEacx1Zmrqv1SmcSdCfffD4ce2tpll21QDNaoqjETfFVG8iWaMaP1cQymU36APPBA676Iw+1ktjSaSbbqQ+qjdeSnmR8qZY1kux9qu2nylM2wnblG0AO89Gu0f8iNdL80mol7NPeDNVdZgo+IcRFxS0T8W0TcGRFnFeNnRMSyiLg3Ir4ZEdsX43cohu8rPp9eVWwdacYMuPDC0VveSJP7UF/266+v9ss52q2DZloGw+2oR3pIvdWtkK1NSFV0qrrwwvbv1IeqU9Wnijq5b8Ro75egM7aHmqiyBf8scHBmvgWYBRwREQcCfwucl5kzgfXACUX5E4D1mfl64LyiXH20Ojm147DjcF/2Qw9t1GO4uo7ky9uKL3yVO41O31G321D/+yq3m2aM9imqTdp5JK0dh9SHU2wPfb29Q8dWdYOiBiJHoQNXRPwGcCPwR8A1wG9l5gsRcRCwMDMPj4jrivc/jYhtgf8AdsshApwzZ04uX768ZXEed9lxLP7F4lJlT9r3JBbNX7TZuAVXLeCiWy8qNf2Z7zqThT0LNxs3/xvzufqeq0tNf+FVsGDF5uNmL4Bbp5aanCu/DvPv6Tcik6mfm8q6jetKTb/8pOXMnjp7s3FxVvmW5cPf2oOpP3/o5eG1E2D3j5eenDxz881ixdoVzLloTqlpp2yAtZ/bfNxVb4Cj/6DcsvddCys2X/Usmg0fm19u+qPuhqu+vnn8C/sWctYNZ5Wa/qQVsOiqfiNmzGDBFw4pv+31wcK+zcfN/xBc/cZSk2/9tvfBK5n/xs3/WU1texfC7C2KxsJyywZ4+LSHmTph6stHQpre9rZY1oopMKfkb7uO2Pa+0W9EZnPbXrv3e0ddyILZCzYbN3vRbG5dd2up6bd62xvBfu/Y3z6Wy467rNT8y4qIFZk57A5v25Yu9ZVBbAOsAF4PfBn4d+DJzHyhKLIG2L14vzvwEECR/H8FTAIe22KeC4AFAJMnT6avr69l8T733HOly65dt/YVy167bm3p6VevXv2K6R9/7PHS07daX19fU/VfvmI5G+7ZMOLl3XTuuey6w64vDz/27GNw8wdKT7/l/+7uDXeXnva5SZPo6/32ZuNuf+wmuPN/l57H1toy/tWrV49oPr+eMoV7/viPWbuub9iyrfSfkydz3yc+wfo5jX3MhhULYOO9paa9/fbbmbBuwmbjmtr2LryADRO2+DVyw9zS0990003susOuHDBlSuP8e5P6ens3G757w91w68mlpt207fXMLR9vlfr6+pra9tq937v77rvp27D59Bs2lN8PbfW2N4L93nPPPdfSPNWMShN8Zr4IzIqInYErgDcNVKz4O9DPoFe03jNzEbAIGi34np6e1gQLXLb6stJlp06ZypbL/vqGr0PJ/cX06dNfMf2kdZPgidIhtFRPTw/br9geSm7rc2bPecUvWW4ov7y3vvWtjVZUYe2GtXBz+em3/N9NWDsByv2IZ/vtt3/F9Bvu3gB3llz47H0hiybsCM+Hb7n8vr4++EXJiU86Ca5stKJ2BN4CTL3qntLbHmeeCb0LNx/3jflQshXFhRfysw1v2KwOE+6ZABvLTf7mN7+Znjf2bDauLdve4sWN0yqPNXeIvCXb3owZHdHJtaenp6ltr937vTe+8Y30zN58+k7f9gba34yazByVF3AmcDqNFvm2xbiDgOuK99cBBxXvty3KxVDznD17drZSb29vS+fXFsOfVRv41cE6er2U/T920f+7jI5eJ01qS12uvz5zxoyRf19b9epgddnGqqgHsDxL5N0qe9HvVrTciYgdgUOAVUAv8P6i2LHA94r3VxbDFJ//uKiIJNVL2Y6FZV/t6riqjlZlL/opQG9E3Ab8DPhhZl4NfBI4LSLuo3GO/eKi/MXApGL8acAZFcZWX37RpbGn2UvL2nH5m0ZdZefgM/M2YJ8Bxt8P7D/A+GeA8r2sNLALL2zuul2/6FL323RpWV9f+873quNU2slObeAXXZKEt6qVJKmWTPCSJNWQh+ilVhpjT6uS1LlswUuSVEMmeGlreFmipA5lgpe2RpPXH/96yhQvS5Q0KjwHL22NTXckK2mZly9KGiW24CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaGrQXfUScVmL6/8xMr/mRJKnDDNWCPx0YD0wY4vXxqgOUJEnNG+o6+K9l5l8PNXFEvLrF8UiSpBYYtAWfmX8+3MRlykiSpNE3bCe7iDg1InaKhosj4taIOGw0gpMkSSNTphf98Zn5FHAYsBvwUeDsSqOSJElbpUyC3/SA63nApZn5b/3GSZKkDlQmwa+IiOtpJPjrImIC8FK1YUmSpK1R5mlyJwCzgPsz8+mImETjML0kSepQZVrwP8zMWzPzSYDMfBw4r9qwJEnS1hjqTnbjgN8Ado2IXfiv8+47AVNHITZJkjRCQ7XgPwasAH6n+Lvp9T3gy9WHJqkTHXkkzJ3bQwSlXkce2e6IpbFp0BZ8Zn4B+EJE/K/M/OIoxiSpg117bbXlJbXGsJ3sMvOLEfFWYHr/8pn51QrjktThMocvE15QK7XNsAk+Ir4GvA5YCbxYjE7ABC9JUocqc5ncHGDPzDK/1yVJUicoc5ncHcBvVR2IJElqnTIt+F2Bn0fELcCzm0Zm5tGVRSVJkrZKmQS/sOog6mpTByNPbkiSRluZXvQ3jEYgkiSpdcr0ot9Ao9c8wPbAdsB/ZuZOVQYmSZJGrkwLfkL/4Yh4D7B/ZRFJkqStVqYX/WYy85+AgyuIRZIktUiZQ/Tv6zf4KhrXxdttTJKkDlamF/38fu9fAFYDx1QSjSSNoiOPhGuv7SlVdt48uOaaauORWmnYQ/SZ+dF+r5My8zOZ+choBCc18+Syujy17Mgjyz2lrU51bpdmHoTjQ3NGz1j83ldh2AQfEdMi4oqIeCQifhkR34mIaaMRnNTqHXA3JE+TzujLHPql0eV3oDXKdLK7FLgSmArsDlxVjJNGTat2wN204zDpvFI3/EBT6/gd2DplEvxumXlpZr5QvC4Ddqs4LmlEhtvpb+KOozt10w80qd3KJPjHIuIjEbFN8foI8HjVgUnNmDevmrLqTP5Ak4ZXJsEfD/w+8B/AOuD9xTipY1xzzfA7/U0ve0JrpDw10JnKnrYZa8r0on8wM4/OzN0y8zcz8z2Z+YvRCE6SOkEzR308NTB6PBo3tDK96BdHxM79hneJiEtKTLdHRPRGxKqIuDMiTi3GL4yIhyNiZfGa12+aT0XEfRFxd0QcPtJKSVIrlT1CpNHVzJG7sbh+ytzoZu/MfHLTQGauj4h9Skz3AvDxzLw1IiYAKyLih8Vn52XmOf0LR8SewAeB36XRY/9HEfGGzHyxVE0kSdLLypyDf1VE7LJpICImUu4hNesy89bi/QZgFY3L7AZzDLAkM5/NzAeA+xgjD7Xx0h9JUquVSfCfA26KiL+JiL8GbgL+rpmFRMR0YB9gWTHqTyLitoi4pN+Ph92Bh/pNtoahfxDUhpf+SJJarUxL/KsRsZzGE+QCeF9m/rzsAiJiPPAd4E8z86mI+ArwNzQeWPM3NH5AHF/M+xWLH2B+C4AFAJMnT6avr69sKMPauHFjS+cHPUD53pu9vUMve+7cxvzKxNj6urRLD1Cuzt2w3HLrpeyyy5U744w3s2zZpOGDAw444HHOPvv2liy3+bJllJ1fq8uV1er5lddYzz0tnWe57aEKPUB7/o+t1s59cWSFPQ8iYjvgauC6zDx3gM+nA1dn5l4R8SmAzPxs8dl1wMLM/Olg858zZ04uX768ZfH29fXR09PTsvk1HmRRrmyZB1ls+qEw3Cpr9XLbqWydu2W5ZbaxsstutlxZrVpus2XLqOp/M9rxVaGqy8DaWZc6dIxrdV4BiIgVmTlnuHKDHqKPiFtLLGTQMhERwMXAqv7JPSKm9Cv2XuCO4v2VwAcjYoeImAHMBG4ZLoZO1q5rsz3kr4GU7WHstcTdrZle5fY4r7ehDtG/KSJuG+LzAF4zxOdvA/4ncHtErCzG/QXwoYiYRePw+2rgYwCZeWdEfAv4OY0e+KfYg37rtLplp3qbN6/8D74DDngcKHfYX1J7DJXgf6fE9IMm4My8kYHPqw+6C8nMzwCfKbFcSS3WzFGkvr7b2XSeVFJnGjTBe7c6SZK6V5nL5CRJUpcxwUuSVEMmeEmSamjQc/ARsYEBbjRDo+NcZuZOlUWlIdn7XZI0nEFb8Jk5ITN3GuA1weTeHs08GrFxGZMkaawq8zQ5ACLiN4Fxm4Yz88FKItKgvIxJklRWmefBHx0R9wIPADfQuDnN9yuOSzVX9gl6kqSRKdPJ7m+AA4F7MnMG8G7gXyqNSrXXzC1yPd0gSc0rk+Cfz8zHaTwX/lWZ2QvMqjgujRFl7ondnqdZSSqj7NG4iEZZjZ4y5+CfLB75+hPgHyPiERr3ipck1VwrT5X5cKvRVaYFfwzwa+DPgB8A/w7MrzIoSVJ7NXPVzrx5Pp2uEw3bgs/M/+w3uLjCWCRJHaKVj7BWe5TpRf++iLg3In4VEU9FxIaIeGo0gpMkSSNT5hz83wHzM3NV1cFIkqTWKHMO/pcmd0lV8r4IUuuVacEvj4hvAv8EPLtpZGZ+t7KoJI0pzfSubqbzlzSWlUnwOwFPA4f1G5eACV5SS9nbWmqdMr3oPzoagUidqOxh4Xnz7HUsqbMMm+Aj4vwBRv8KWJ6Z32t9SFL7zZvX3GFjb+AhqdOU6WQ3jsatae8tXnsDE4ETIuLzFcYmtc0115S7ja6HlCV1qjLn4F8PHJyZLwBExFeA64FDAW8SLqlj2fNeY1mZFvzuwKv7Db8amJqZL9KvV726lw+KUN0009PepxWqrsok+L8DVkbEpRFxGfCvwDkR8WrgR1UGp2o1e7mR55nVLZo5xdKupxX6FDZVrUwv+osj4lpgfyCAv8jMtcXHp1cZnKrVTK9vD3VKrWUnTlVt0AQfEb+TmXdFxL7FqIeKv78VEb+VmbdWH56kOvAH4uCG66jp/04jNVQL/jRgAfC5AT5L4OBKIpJUG81cbugd6qTWGjTBZ+aC4u/c0QtHUrM6uYXnzX+k9inzuNgPRMSE4v1fRsR3I2Kf6kOTNJRmWry2jqWxp8x18H+Vmf8vIt4OHA6cA1wAHFBpZJKGZOtY0lDKXCb3YvH3SOArxe1pt68uJKk7lbnc6Ywz3tzuMCWNEWUS/MMRcSHw+8C1EbFDyemkMaGZw9/Llk2qLhCpC5T5IazWKJOofx+4DjgiM5+kcR96r3+XCmVvqiKNZc32A/EOg1uvzI1unqbfs98zcx2wrsqgpDqzhaKxqNk+I319twM9VYQyZnioXS1V9vabY5G93iWNJhO8WqqZW2qOtSR2zTXQ29tX6nC+PeQlba0yl8lJTfOcs8aysXqUSp3FFrwktUizR6XG2lEsjS5b8GqKLRNpcJ5aUSexBa9S7CAmSd3FFrxKsWUiSd3FFrwkSTVkgpckqYZM8JIk1ZDn4CWpC3gFi5plC16SOlgzV6X4gBb1V1mCj4g9IqI3IlZFxJ0RcWoxfmJE/DAi7i3+7lKMj4g4PyLui4jbImLfqmKTpG5R9mmFmXD22be3O1x1kCpb8C8AH8/MNwEHAqdExJ7AGcDSzJwJLC2GAX4PmFm8FgBfqTA2SZJqrbIEn5nrMvPW4v0GYBWwO3AMsLgothh4T/H+GOCr2XAzsHNETKkqPkmS6ixyFJ4KEhHTgZ8AewEPZubO/T5bn5m7RMTVwNmZeWMxfinwycxcvsW8FtBo4TN58uTZS5YsaVmcGzduZPz48S2bXztZl85Ul7rUpR5gXTpVXepSRT3mzp27IjPnDFeu8l70ETEe+A7wp5n5VAzeFXSgD17x6yMzFwGLAObMmZM9PT0tihT6+vpo5fzaybp0prrUpS71AOvSqepSl3bWo9Je9BGxHY3k/o+Z+d1i9C83HXov/j5SjF8D7NFv8mnA2irjkySprqrsRR/AxcCqzDy330dXAscW748Fvtdv/B8WvekPBH6Vmeuqik+SpDqr8hD924D/CdweESuLcX8BnA18KyJOAB4EPlB8di0wD7gPeBr4aIWxSZJUa5Ul+KKz3GAn3N89QPkETqkqHkmSxhLvZCdJUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVUGUJPiIuiYhHIuKOfuMWRsTDEbGyeM3r99mnIuK+iLg7Ig6vKi5JksaCKlvwlwFHDDD+vMycVbyuBYiIPYEPAr9bTPP3EbFNhbFJklRrlSX4zPwJ8ETJ4scASzLz2cx8ALgP2L+q2CRJqrvIzOpmHjEduDoz9yqGFwLHAU8By4GPZ+b6iPgScHNmXl6Uuxj4fmZ+e4B5LgAWAEyePHn2kiVLWhbvxo0bGT9+fMvm107WpTPVpS51qQdYl05Vl7pUUY+5c+euyMw5wxbMzMpewHTgjn7Dk4FtaBw5+AxwSTH+y8BH+pW7GPjvw81/9uzZ2Uq9vb0tnV87WZfOVJe61KUemdalU9WlLlXUA1ieJXLwqPaiz8xfZuaLmfkScBH/dRh+DbBHv6LTgLWjGZskSXUyqgk+Iqb0G3wvsKmH/ZXAByNih4iYAcwEbhnN2CRJqpOqMYLUAAAIfklEQVRtq5pxRHwD6AF2jYg1wJlAT0TMAhJYDXwMIDPvjIhvAT8HXgBOycwXq4pNkqS6qyzBZ+aHBhh98RDlP0PjvLwkSdpK3slOkqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqqHKEnxEXBIRj0TEHf3GTYyIH0bEvcXfXYrxERHnR8R9EXFbROxbVVySJI0FVbbgLwOO2GLcGcDSzJwJLC2GAX4PmFm8FgBfqTAuSZJqr7IEn5k/AZ7YYvQxwOLi/WLgPf3GfzUbbgZ2jogpVcUmSVLdjfY5+MmZuQ6g+PubxfjdgYf6lVtTjJMkSSOwbbsDKMQA43LAghELaBzGB9gYEXe3MI5dgcdaOL92si6dqS51qUs9wLp0qrrUpYp6/HaZQqOd4H8ZEVMyc11xCP6RYvwaYI9+5aYBaweaQWYuAhZVEVxELM/MOVXMe7RZl85Ul7rUpR5gXTpVXerSznqM9iH6K4Fji/fHAt/rN/4Pi970BwK/2nQoX5IkNa+yFnxEfAPoAXaNiDXAmcDZwLci4gTgQeADRfFrgXnAfcDTwEerikuSpLGgsgSfmR8a5KN3D1A2gVOqiqUJlRz6bxPr0pnqUpe61AOsS6eqS13aVo9o5FZJklQn3qpWkqQaMsEXIuKIiLi7uF3uGcNP0bkiYnVE3B4RKyNiebvjaUYztzjuZIPUY2FEPFysl5URMa+dMZYVEXtERG9ErIqIOyPi1GJ8N66XwerSVesmIsZFxC0R8W9FPc4qxs+IiGXFOvlmRGzf7liHM0RdLouIB/qtk1ntjrWMiNgmIv41Iq4uhtu2TkzwNFYI8GUat8zdE/hQROzZ3qi22tzMnNWFl5lcRvlbHHeyy3hlPQDOK9bLrMy8dpRjGqkXgI9n5puAA4FTiu9HN66XweoC3bVungUOzsy3ALOAI4orkP6WRj1mAuuBE9oYY1mD1QXg9H7rZGX7QmzKqcCqfsNtWycm+Ib9gfsy8/7MfA5YQuP2uRplTd7iuGMNUo+ulJnrMvPW4v0GGjuv3enO9TJYXbpKcVvvjcXgdsUrgYOBbxfju2WdDFaXrhMR04AjgX8ohoM2rhMTfEPdbpWbwPURsaK481+3G+wWx93oT4onJl7SDYe0txQR04F9gGV0+XrZoi7QZeumOBS8ksYNw34I/DvwZGa+UBTpmv3YlnXJzE3r5DPFOjkvInZoY4hlfR74c+ClYngSbVwnJviG0rfK7RJvy8x9aZxyOCUi3tnugAQ0npL4OhqHIdcBn2tvOM2JiPHAd4A/zcyn2h3P1higLl23bjLzxcycRePOn/sDbxqo2OhGNTJb1iUi9gI+BfwOsB8wEfhkG0McVkQcBTySmSv6jx6g6KitExN8Q+lb5XaDzFxb/H0EuILGl7+b/XLT0wW3uMVxV8nMXxY7speAi+ii9RIR29FIiP+Ymd8tRnflehmoLt28bjLzSaCPRp+CnSNi0/1Num4/1q8uRxSnUzIznwUupfPXyduAoyNiNY3TvAfTaNG3bZ2Y4Bt+BswsejtuD3yQxu1zu05EvDoiJmx6DxwG3DH0VB1vsFscd5XY/BHI76VL1ktxHvFiYFVmntvvo65bL4PVpdvWTUTsFhE7F+93BA6h0Z+gF3h/Uaxb1slAdbmr34/HoHHeuqPXSWZ+KjOnZeZ0Gjnkx5n5Ydq4TrzRTaG4LObzwDbAJZn5mTaHNCIR8VoarXZo3Knw691Ul+h3i2PglzRucfxPwLeA/0Zxi+PM7OgObIPUo4fGIeAEVgMf64ZnLkTE24F/Bm7nv84t/gWNc9fdtl4Gq8uH6KJ1ExF70+iwtQ2Nhtq3MvOvi+//EhqHtP8V+EjRAu5YQ9Tlx8BuNA5zrwRO7tcZr6NFRA/wicw8qp3rxAQvSVINeYhekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvKTKRERfRHTbA4+kWjDBS5JUQyZ4aYwp7nZ4TfH87Tsi4n9ExKcj4mfF8KLi7mGbWuDnRcRPimeo7xcR3y2ebf1/ijLTI+KuiFhcPBjk2xHxGwMs97CI+GlE3BoR/6+4H7ykipjgpbHnCGBtZr4lM/cCfgB8KTP3K4Z3BI7qV/65zHwncAGN22yeAuwFHBcRk4oybwQWZebewFPAH/dfYETsCvwlcEjxIKTlwGmV1VCSCV4ag24HDomIv42Id2Tmr4C5EbEsIm6n8ZCM3+1X/sp+091ZPATkWeB+/ushTQ9l5r8U7y8H3r7FMg8E9gT+pXgs6LHAb7e8ZpJetu3wRSTVSWbeExGzgXnAZyPiehqt8jmZ+VBELATG9Ztk032zX+r3ftPwpn3Ilve83nI4aDzn+0MtqIKkEmzBS2NMREwFns7My4FzgH2Ljx4rzou/f9CJB/ffIuKg4v2HgBu3+Pxm4G0R8foiht+IiDeMYDmSSrIFL409bwb+b0S8BDwP/BGNx3HeTuNJaj8bwTxXAcdGxIXAvcBX+n+YmY9GxHHANyJih2L0XwL3jKQCkobn0+QkbZWImA5cXXTQk9QhPEQvSVIN2YKXJKmGbMFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphv4/uKsUEnHhVXQAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6104\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAF3CAYAAABJzllyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XuYHWWd4PHvT25BEwQCm0kIO4kaHRnEQHP12kG5TEJAXZnV0RkQMLrL7sMOyog+MxJm1tXZRVCUUcNyyciM0VUZuSlg7NZhFJiEyXAx3AYihGTkFiQRuf/2j1PB7iTdXaf7VJ9zqr+f5zlPn3rPW3V+b9c59Ttv1VtVkZlIkqR6eVm7A5AkSa1ngpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqqPIEHxHbRcS/RMRVxfTsiLgpIu6JiG9GxI5F+U7F9L3F67Oqjk2SpLoajx78acDqAdN/DZyXmXOADcDJRfnJwIbMfA1wXlFPkiSNQqUJPiJmAguA/1tMB3A48O2iylLgXcXz44ppitffUdSXJElNqroH/wXgz4AXi+mpwBOZ+XwxvRbYq3i+F/AgQPH6r4r6kiSpSdtXteCIOAZ4ODNXRkTv5uJtVM0Srw1c7iJgEcDOO+/cs/fee7cg2oYXX3yRl72sHuMObUtnqktb6tIOsC2dqi5tqaIdd99996OZueeIFTOzkgfwWRo99DXAvwNPAX8HPApsX9Q5DLi2eH4tcFjxfPuiXgz3Hj09PdlKfX19LV1eO9mWzlSXttSlHZm2pVPVpS1VtANYkSXycGU/jzLzk5k5MzNnAe8DfpSZHwD6gPcW1U4Avlc8v6KYpnj9R0VDJElSk9qx/+MTwOkRcS+NY+wXFeUXAVOL8tOBM9sQmyRJtVDZMfiBMrMf6C+e3wccvI06TwPHj0c8kiTV3bgkeEmSmvXcc8+xdu1ann766XaHMmqvfOUrWb169cgVt2HSpEnMnDmTHXbYYVTzm+AlSR1p7dq1TJkyhVmzZtGtl0XZuHEjU6ZMaXq+zOSxxx5j7dq1zJ49e1Tv3f3nIEiSaunpp59m6tSpXZvcxyIimDp16pj2XpjgJUkdayIm983G2nYTvCRJQ9huu+2YO3cu++67LwsXLuSJJ5546bWlS5cyZ84c5syZw9KlS4dZSnuY4CVJGsLOO+/MqlWruP3229l999254IILAHj88cc5++yzuemmm7j55ps5++yz2bBhQ5ujHcwEL0lSCYcddhgPPfQQANdeey1HHHEEu+++O7vtthtHHHEEP/jBD9oc4WAmeElSx4uo5lHWCy+8wPLlyzn22GMBeOihhxh4L5SZM2e+lPw7hQlekqQh/OY3v2Hu3LlMnTqVxx9/nCOOOAKAbV1JvdMGBJrgJUkdL7Oax0g2H4P/xS9+wbPPPvvSMfiZM2fy4IMPvlRv7dq1zJgxo6rmj4oJXpKkEbzyla/k/PPP55xzzuG5557jqKOO4rrrrmPDhg1s2LCB6667jqOOOqrdYQ5igpckqYT999+fN77xjSxbtozdd9+dv/iLv+Cggw7ioIMO4tOf/jS77757u0McxEvVSpI0hE2bNg2avvLKK196ftJJJ3HSSSeNd0il2YOXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUkawnC3iz366KPZddddOeaYY9oY4dBM8JIkDWGo28UCnHHGGXz9619vY3TDM8FLkjpbVbeSa/KWcgNvFwvwjne8gylTplTR4pYwwUuSNIItbxfbDUzwkiQNYajbxXYDE7wkSUMY6nax3cAEL0nqbFXdDL7sTeHZ+nax3cAEL0lSCQNvFwvw1re+leOPP57ly5czc+ZMrr322jZHOJi3i5UkaQjD3S72H//xH8c7nKbYg5ckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSWrS+eefz+tf/3o+8IEP8Mwzz/DOd76TuXPn8s1vfrOp5fT39/PTn/60khg9D16SpCb9zd/8Dd///veZPXs2N954I8899xyrVq1qejn9/f1MnjyZN73pTS2P0QQvSdIwzj33XC6++GIATjnlFO68807uu+8+jj32WD74wQ9y4YUX8sgjjzB37ly+853vcOGFF3LFFVew/fbb09vby/nnn88jjzzCRz/6UR544AEAvvCFL7DXXnvx1a9+le22247LLruML33pS7z1rW9tWdwmeElSV1jcv5izf3x2qbofPuDDLFm4ZFDZoisXceEtF740fdbbz2Jx7+Jhl7Ny5UouueQSbrrpJjKTQw45hMsuu4wf/OAH9PX1sccee3DIIYdwzjnncNVVV/H4449z+eWXc+eddxIRPPjggwCcdtpp/Omf/ilvectbeOCBBzjqqKNYvXo1H/3oR5k8eTIf//jHm/tnlGCClyRpCDfccAPvfve7ecUrXgHAe97znmEvUbvLLrswadIkTjnlFBYsWMDb3/52AH74wx/y85///KV6Tz75JBs3bqw0dhO8JElDyJJ3m9ts++235+abb2b58uUsW7aML37xi/z4xz/mxRdf5Gc/+xk777xzRZFuI5ZxeydJksZgce/iEXepD2fJwiVb7bYfydve9jZOPPFEzjzzTDKTyy+/nK9//euce+6526y/adMmnnrqKebPn8+hhx7Kq1/9agCOPPJIvvzlL3PGGWcAsGrVKubOncuUKVN48sknR92m4XianCRJQzjggAM48cQTOfjggznkkEM45ZRT2H///Yesv3HjRo455hj2228/3v72t/PZz34WaJxWt2LFCvbbbz/22WcfvvrVrwKwcOFCLr/8cubOndvyu9PZg5ckaRinn346p59++qCyNWvWvPS8t7eX3t5eAKZPn87NN9/80mubj7Pvscce2zxH/rWvfS233npr64PGHrwkSbVkgpckqYZM8JIk1VBlCT4iJkXEzRHxrxFxR0ScXZRfGhH3R8Sq4jG3KI+IOD8i7o2IWyPigKpikyR1h2ZPU6uTsba9ykF2zwCHZ+amiNgBuCEivl+8dkZmfnuL+n8AzCkehwBfKf5KkiagSZMm8dhjjzF16lQiot3hjKvM5LHHHmPSpEmjXkZlCT4bPz02FZM7FI/hfo4cB/xtMd+NEbFrREzPzPVVxShJ6lwzZ85k7dq1PPLII+0OZdSefvrpUSfpSZMmMXPmzFG/d1S5+yMitgNWAq8BLsjMT0TEpcBhNHr4y4EzM/OZiLgK+Fxm3lDMuxz4RGau2GKZi4BFANOmTetZtmxZy+LdtGkTkydPbtny2sm2dKa6tKUu7QDb0qnq0pYq2jFv3ryVmXngiBUzs/IHsCvQB+wLTAcC2AlYCny6qHM18JYB8ywHeoZbbk9PT7ZSX19fS5fXTralM9WlLXVpR6Zt6VR1aUsV7QBWZIncOy6j6DPzCaAfODoz1xcxPgNcAhxcVFsL7D1gtpnAuvGIT5KkuqlyFP2eEbFr8Xxn4J3AnRExvSgL4F3A7cUsVwB/UoymPxT4VXr8XZKkUalyFP10YGlxHP5lwLcy86qI+FFE7EljN/0q4KNF/WuA+cC9wFPAhyqMTZKkWqtyFP2twFZX5M/Mw4eon8CpVcUjSdJE4pXsJEmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPAT1fXXw6teBRHlH696VWM+SVLHM8FPVB/5CNx/f3Pz3H9/Yz5JUsczwU9UzSb3sc4nSRpXJnhJkmpo+3YHoA6ROfRrEeMXhySpJezBS5JUQyZ4SZJqyASv5nk6nSR1PBO8Ws/T6SSp7UzwGtns2c3P4+l0ktRWJniN7GtfG12SlyS1jafJaWRHHAH33TdyPU+nk6SOYQ++boprzPfOmzf8QDip2zR7/4TxHuw5Ee/vMJo2l3yMuA2ry/+wQib4uhnNNear0mlfzAo3Rm5kxkGzn+3xHuw5Ee/v0Anbm27/H1bIBF+VdvU2RvNla9fx9W7YAJd1//1w5JEjrudD/uiP/CEwWqNZd61Y32W/yxPx/g6dEnunxNFhTPBV6fTexmazZzcG0bVqWc0azy9mB2wEdl6/vjN7G03+IC29+7QOezc6oZcqjYIJviqj7W20emOZOfzjvvsag+hawdH25XRisuiAvRst+2xvaaw/Vka7V2y4712TMXbNj6SRtjdNPPr7+pr7H2orJvhu08nHmzaPth/py9spWrEhuu660f2oGa8NetW7l1tppB8C3aKVe8U2G+9DQGU/N+polSX4iJgUETdHxL9GxB0RcXZRPjsiboqIeyLimxGxY1G+UzF9b/H6rKpia4uhEkSn79bW8Kr6UdOqH3ITcfdylXuRRuqZl90rVlGMIx4CKpu4jzxy4n1uaqjKHvwzwOGZ+UZgLnB0RBwK/DVwXmbOATYAJxf1TwY2ZOZrgPOKevXnbu2JYbQ/5Mba069i93KZ3adj2btRNsahVPWdamXPvMrv/XCfmyoTt9uxjhM5DrtMI+LlwA3AfwGuBn4nM5+PiMOAxZl5VERcWzz/WURsD/w7sGcOE+CBBx6YK1asaFmcJ156Ikt/sbRU3Q8f8GGWLFwyqGzRlYu48JYLS81/1tvPYnHv4kFlC7+xkKvuvqrU/F875mss6lk0qKxnSQ+3rL+l1PxXvO8KFr5u4aCyGZ+fwfpN60vNv+LDK+iZ0TOoLM4uv8vuoc/DjI2/nV43Bfb6WOnZybMGfCyuv56VnzyRAxeuKzXv9I2w7vMM6llfedeVHLvs2FLzHzD9AFYuWjmobMnKJXzkqnI97mPugiu/MbhscS+c3Vtqdj68EpZcObhs0UK4sGfb9bd0Vj8s7h9ctvD9cNXrys3/tWO+xms3vpbe3t6Xysb02bv+emb88A9Y//IXSs2/4soZ9Hz20kG95KY+e6c/xIwpM16a/va13+b4G48vPf+gzx6wct1KDrzwwFLzTp88nXUfG/w5beqztw5WDt7ssKQHPrJw2/W31PbPXiu3e8Uhgp5FcMuMYWYaoB3bvRN+9wQuPfHScgGWFBErM3PED12lV7KLiO2AlcBrgAuAfwOeyMzniyprgb2K53sBDwIUyf9XwFTg0S2WuQhYBDBt2jT6+/tbFu+zzz5buu669eu2eu9168slGIA1a9ZsNf9jjz5Wev677rqL/o2D59+4ceO2K2/DbbfdxpT1UwaVNdP+FStXsPHu8u/XagP/d4eccAJQ7gs61DJue/S20vNt3Lhxq3V317q7Ss//2GGH0b/ofw0qW7PmUij543Ks1pxwAv1nnTg4pts+BY//rNT8d911FzN2mTHofzCmz94OO/DsrrvCs+U+/yvO+jQbd9gBRvnd/+lPf8oeO+3x0vSvf/3rpubfat1vLL/un3322a3mb+qz99o59PcNzvB3rbsS7jm39DLGYt0xC+g//eODy+4+B9ZfXWr+Vm73ekvP9Vvt2O5ta52Pl0oTfGa+AMyNiF2By4HXb6ta8XdbP4O26r1n5hJgCTR68AN7EWN16ZpLS9edMX0GW77332/8+9J5ZtasWVvNP3X9VHi83Pyve93r6O0ZPP+Uu6fApnLzv+ENb6D3dYPn33HljlDys35gz4Fb/ZLlx+XmbYXeefMGF0xvcgGzZw/6/2+8ayPcUW7WKVOmbLXu7l55N9xTbv6pe0zdav7+/n74Rbn5x6oVn73JGycPWkY3ffbe9KY3bdWDb8aW/7sp66ZAuZ0X7LjjjlvNP56fvZdsPtxwxBHQvxh+fHap2Tpxu9eMdnz2trXOx01mjssDOAs4g0aPfPui7DDg2uL5tcBhxfPti3ox3DJ7enqylfr6+lq6vJba8gjnWOu103XXZc6e3Yox7KN7zJ7diKENWvYZa/Z/2OI2d/R3pUm2pTON2JZu2NZlNesEWJEl8m5lPfiI2BN4LjOfiIidgXfSGDjXB7wXWAacAHyvmOWKYvpnxes/Khqiuil785pXvar5AUGzZ5dbdrcr+z+UNGFVOYp+OtAXEbcC/wxcn5lXAZ8ATo+Ie2kcY7+oqH8RMLUoPx04s8LY1A2aHWlcxfnHktSlKuvBZ+atwP7bKL8POHgb5U8D5Yeyqv6KXmp/f3/7jmFJUpfyfvCSpO5X9sp6E+jIr5eqlSSphkzw3cprQ0uayLxy3ohM8JKk7uNlvkdkgu8Wo/kg++GXVFdlb/S05WMCMcF3iyZ/rf5m+nRPGZOkCcxR9N2iyQub3OSpZZI0odmDlySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg0NeZpcRJxeYv5fZ6YnW0uS1GGG68GfAUwGpgzz+FjVAUqSpOYNd6Gbr2fmXw43c0S8osXxSJKkFhiyB5+ZfzbSzGXqSJKk8TfiILuIOC0idomGiyLilog4cjyCkyRJo1NmFP1JmfkkcCSwJ/Ah4HOVRiVJksakTIKP4u984JLM/NcBZZIkqQOVSfArI+I6Ggn+2oiYArxYbViSJGksytwu9mRgLnBfZj4VEVNp7KaXJEkdqkwP/vrMvCUznwDIzMeA86oNS5IkjcVwV7KbBLwc2CMiduO3x913AWaMQ2ySJGmUhttF/xHgf9BI5iv5bYJ/Erig4rgkSdIYDJngM/OLwBcj4r9n5pfGMSZJkjRGIw6yy8wvRcSbgFkD62fm31YYlyRJGoMRE3xEfB14NbAKeKEoTsAEL0lShypzmtyBwD6ZmVUHI0mSWqPMaXK3A79TdSCSJKl1yvTg9wB+HhE3A89sLszMYyuLSpIkjUmZBL+46iAkSVJrlRlF/+PxCESSJLVOmVH0G2mMmgfYEdgB+HVm7lJlYJIkafTK9OCnDJyOiHcBB1cWkSRJGrMyo+gHycx/AA6vIBZJktQiZXbRv2fA5MtonBfvOfGSJHWwMqPoFw54/jywBjiukmgkSVJLlDkG/6HxCESSJLXOiMfgI2JmRFweEQ9HxC8j4jsRMXM8gpMkSaNTZpDdJcAVNO4LvxdwZVEmSZI6VJkEv2dmXpKZzxePS4E9K45LkiSNQZkE/2hEfDAitiseHwQeqzowSZI0emUS/EnAHwL/DqwH3luUSZKkDlVmFP0DgHeOkySpi5QZRb80InYdML1bRFxcYr69I6IvIlZHxB0RcVpRvjgiHoqIVcVj/oB5PhkR90bEXRFx1GgbJUnSRFfmQjf7ZeYTmycyc0NE7F9ivueBj2XmLRExBVgZEdcXr52XmecMrBwR+wDvA36fxoj9H0bEazPzhVItkSRJLylzDP5lEbHb5omI2J1yu/bXZ+YtxfONwGoap9kN5ThgWWY+k5n3A/fiTW0kSRqVMgn+88BPI+KvIuIvgZ8C/7uZN4mIWcD+wE1F0X+LiFsj4uIBPx72Ah4cMNtahv9BIKkmFiyAiHKPBQvaHa3UHSJz5PvGFLvPDwcCWJ6ZPy/9BhGTgR8Dn8nM70bENOBRGjes+StgemaeFBEXAD/LzMuK+S4CrsnM72yxvEXAIoBp06b1LFu2rGwoI9q0aROTJ09u2fLaybZ0prq0pdXtmDevt6n6fX39LXvvuqwTsC2dqIp2zJs3b2VmHjhixcys7AHsAFwLnD7E67OA24vnnwQ+OeC1a4HDhlt+T09PtlJfX19Ll9dOtqUz1aUtrW4HNB6tqteMuqyTTNvSiapoB7AiS+TgIXfRR8QtI/04GK5ORARwEbA6M88dUD59QLV3A7cXz68A3hcRO0XEbGAOcPNIMUiSpK0NN1ju9RFx6zCvB/DKYV5/M/DHwG0Rsaoo+xTw/oiYS2MX/RrgIwCZeUdEfAv4OY0R+KemI+glSRqV4RL875WYf8gEnJk30PgRsKVrhpnnM8BnSryvJEkaxpAJPjN/MZ6BSJKk1ilzmpwkSeoyJnhJkmrIBC9JUg0NeQw+IjbSGOm+1Us0TkTdpbKoJEnSmAw3yG7KeAYiSZJap8zd5ACIiP8ATNo8nY37xEuSpA5U5n7wx0bEPcD9NK4pvwb4fsVxSZKkMSgzyO6vgEOBuzNzNvAO4J8qjUqSJI1JmQT/XGY+RuO+8C/LzD5gbsVxSZKkMShzDP6J4pavPwH+LiIepnGteEmS1KHK9OCPA34D/CnwA+DfgIVVBiVpZAsWQES5x4IF7Y62u/m/VjcaMcFn5q8z84XMfD4zl2bm+cUue0ltdM2Qt20aW11trYr/tT8aVLUyo+jfExH3RMSvIuLJiNgYEU+OR3BSWRN5Y5k5/EOt08r/tT/QVLUyu+j/N3BsZr4yM3fJzClexU6dpp0by2Z+XJx55hta++bqep36A20i/2iuizKD7H6Zmasrj0RqgZE2iBGtf89mfjDcdNPU1gcgVcA9DN2vTIJfERHfBP4BeGZzYWZ+t7KopC7Ujh8XUtX8XHevMgl+F+Ap4MgBZQmY4KVRKLNBnD8frr66+lgk1deICT4zPzQegUh1N39++V2Z11xjz0gT14IFcM01vaXq+mN4aCMm+Ig4fxvFvwJWZOb3Wh+SVE9XXw39/f309vYOW6+xcSu3zPnzxx7XRNZMItH48fh/a5QZRT+JxqVp7yke+wG7AydHxBcqjE1qm2ZGELfa1VePPLJ688Oey9g0kxz8MTX+OvUMg25R5hj8a4DDM/N5gIj4CnAdcARwW4WxSW3TbK/AjX93M1l0t7I/tCfaei7Tg98LeMWA6VcAMzLzBQaMqpfqyJ60RqMde34WLIB583on1Hnr/rAeXtkL3ayKiEsi4lLgX4BzIuIVwA+rDE6jMxG/6FInaDbhtDJBNXvcutU/QNqxvWnmcNZE671DuVH0F0XENcDBQACfysx1xctnVBmcRscBKlJ7dMKenJESWasHcTZ7dojGz5AJPiJ+LzPvjIgDiqIHi7+/ExG/k5m3VB+eBmrmiwn1uECFo5yl1mr1j5Cyy9u8vemG7U5dDLeL/vTi7+e38Tin4ri0Dc0k90MOqccN/xzl3HmaOQRU1WEgDz91n2a/n3XZhrXTkD34zFxU/J03fuGojDLHkvr7bwN6qw5l3EzE42edqtndrK3cLevu4O7V7J6Dum3D2qHM7WKPj4gpxfM/j4jvRsT+1YcmqZO1Y1BT2UFVksqNov+LzNwYEW8BjgKWAl+tNixJkjQWZRL8C8XfBcBXisvT7lhdSBpv3vNZkuqnTIJ/KCK+BvwhcE1E7FRyPnW4Zge9eExTkrpHmUT9h8C1wNGZ+QSN69B7/nsNNHORCElSdxkxwWfmU5n53cy8p5hen5nXVR+aOlGrduWXvZmLpAa/K2qWu9pVSjO788vsyp+I5/RLo9HMd8/vigYqczc5qemrVZU1Ec/pb5cy62b+/M643Kp+q5n14XdFA9mD14Q0kXZ3tnrvy2hMpP+31ClM8JpQmkl2dbn0bTsvDjMR/99Sp3AXvSYUdz+PrJU9av/fUvvYg5cEOJhLqht78JIAB3NJdWMPXpKkGjLBS5JUQyZ4SZJqyGPwkmrLc+w1kdmDl1Q7nhEgVZjgI2LviOiLiNURcUdEnFaU7x4R10fEPcXf3YryiIjzI+LeiLg1Ig6oKjZJ9dbMnRI/97nb2h2uVIkqe/DPAx/LzNcDhwKnRsQ+wJnA8sycAywvpgH+AJhTPBYBX6kwNkmSaq2yBF/cVvaW4vlGYDWwF3AcsLSothR4V/H8OOBvs+FGYNeImF5VfJIk1VlkVRehHvgmEbOAnwD7Ag9k5q4DXtuQmbtFxFXA5zLzhqJ8OfCJzFyxxbIW0ejhM23atJ5ly5a1LM5NmzYxefLkli2v1ebN6wWgr69/xLrtakvZGLuhLVWoS1vq0g6wLZ2qLm2poh3z5s1bmZkHjlgxMyt9AJOBlcB7iukntnh9Q/H3auAtA8qXAz3DLbunpydbqa+vr6XLa7XNRw3LaFdbysbYDW2pQl3aUpd2ZNqWTlWXtlTRDmBFlsi/lY6ij4gdgO8Af5eZ3y2Kf7l513vx9+GifC2w94DZZwLrqoxPkqS6qnIUfQAXAasz89wBL10BnFA8PwH43oDyPylG0x8K/Coz11cVnyRJdVblhW7eDPwxcFtErCrKPgV8DvhWRJwMPAAcX7x2DTAfuBd4CvhQhbFJklRrlSX4bAyWG+o6Uu/YRv0ETq0qHkmSJhKvZFehCC+VKUlqDxO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywY+DzTedGe6xYEG7o5Qk1YkJvkLz55eve8011cUhSZp4KrsfvODqq8vV23xLWW8tK0lqFXvwHaCZnn4zdSVJE5c9+A5QtqcvSVJZ9uAlSaohE7wkSTVkgpckqYZM8JIk1ZCD7FQJT/mTpPayB6+W8pQ/SeoM9uDVUp7yJ0mdwR68JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBqqLMFHxMUR8XBE3D6gbHFEPBQRq4rH/AGvfTIi7o2IuyLiqKrikiRpIqiyB38pcPQ2ys/LzLnF4xqAiNgHeB/w+8U8fxMR21UYmyRJtVZZgs/MnwCPl6x+HLAsM5/JzPuBe4GDq4pNkqS6i8ysbuERs4CrMnPfYnoxcCLwJLAC+FhmboiILwM3ZuZlRb2LgO9n5re3scxFwCKAadOm9Sxbtqxl8W7atInJkye3bHntZFs6U13aUpd2gG3pVHVpSxXtmDdv3srMPHDEiplZ2QOYBdw+YHoasB2NPQefAS4uyi8APjig3kXAfxpp+T09PdlKfX19LV1eO9mWzlSXttSlHZm2pVPVpS1VtANYkSVy8LiOos/MX2bmC5n5InAhv90NvxbYe0DVmcC68YxNkqQ6GdcEHxHTB0zVj+/AAAAImElEQVS+G9g8wv4K4H0RsVNEzAbmADePZ2ySJNXJ9lUtOCK+AfQCe0TEWuAsoDci5gIJrAE+ApCZd0TEt4CfA88Dp2bmC1XFJklS3VWW4DPz/dsovmiY+p+hcVxekiSNkVeykySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQaMsFLklRDJnhJkmqosgQfERdHxMMRcfuAst0j4vqIuKf4u1tRHhFxfkTcGxG3RsQBVcUlSdJEUGUP/lLg6C3KzgSWZ+YcYHkxDfAHwJzisQj4SoVxSZJUe5Ul+Mz8CfD4FsXHAUuL50uBdw0o/9tsuBHYNSKmVxWbJEl1N97H4Kdl5nqA4u9/KMr3Ah4cUG9tUSZJkkZh+3YHUIhtlOU2K0YsorEbH2BTRNzVwjj2AB5t4fLaybZ0prq0pS7tANvSqerSlira8btlKo13gv9lREzPzPXFLviHi/K1wN4D6s0E1m1rAZm5BFhSRXARsSIzD6xi2ePNtnSmurSlLu0A29Kp6tKWdrZjvHfRXwGcUDw/AfjegPI/KUbTHwr8avOufEmS1LzKevAR8Q2gF9gjItYCZwGfA74VEScDDwDHF9WvAeYD9wJPAR+qKi5JkiaCyhJ8Zr5/iJfesY26CZxaVSxNqGTXf5vYls5Ul7bUpR1gWzpVXdrStnZEI7dKkqQ68VK1kiTVkAm+EBFHR8RdxeVyzxx5js4VEWsi4raIWBURK9odTzOaucRxJxuiHYsj4qFivayKiPntjLGsiNg7IvoiYnVE3BERpxXl3bhehmpLV62biJgUETdHxL8W7Ti7KJ8dETcV6+SbEbFju2MdyTBtuTQi7h+wTua2O9YyImK7iPiXiLiqmG7bOjHB01ghwAU0Lpm7D/D+iNinvVGN2bzMnNuFp5lcSvlLHHeyS9m6HQDnFetlbmZeM84xjdbzwMcy8/XAocCpxfejG9fLUG2B7lo3zwCHZ+YbgbnA0cUZSH9Nox1zgA3AyW2Msayh2gJwxoB1sqp9ITblNGD1gOm2rRMTfMPBwL2ZeV9mPgsso3H5XI2zJi9x3LGGaEdXysz1mXlL8XwjjY3XXnTnehmqLV2luKz3pmJyh+KRwOHAt4vyblknQ7Wl60TETGAB8H+L6aCN68QE31C3S+UmcF1ErCyu/NfthrrEcTf6b8UdEy/uhl3aW4qIWcD+wE10+XrZoi3QZeum2BW8isYFw64H/g14IjOfL6p0zXZsy7Zk5uZ18plinZwXETu1McSyvgD8GfBiMT2VNq4TE3xD6Uvldok3Z+YBNA45nBoRb2t3QAIad0l8NY3dkOuBz7c3nOZExGTgO8D/yMwn2x3PWGyjLV23bjLzhcycS+PKnwcDr99WtfGNanS2bEtE7At8Evg94CBgd+ATbQxxRBFxDPBwZq4cWLyNquO2TkzwDaUvldsNMnNd8fdh4HIaX/5u9svNdxfc4hLHXSUzf1lsyF4ELqSL1ktE7EAjIf5dZn63KO7K9bKttnTzusnMJ4B+GmMKdo2Izdc36brt2IC2HF0cTsnMfAa4hM5fJ28Gjo2INTQO8x5Oo0fftnVigm/4Z2BOMdpxR+B9NC6f23Ui4hURMWXzc+BI4Pbh5+p4Q13iuKvE4Fsgv5suWS/FccSLgNWZee6Al7puvQzVlm5bNxGxZ0TsWjzfGXgnjfEEfcB7i2rdsk621ZY7B/x4DBrHrTt6nWTmJzNzZmbOopFDfpSZH6CN68QL3RSK02K+AGwHXJyZn2lzSKMSEa+i0WuHxpUK/76b2hIDLnEM/JLGJY7/AfgW8B8pLnGcmR09gG2IdvTS2AWcwBrgI91wz4WIeAvwj8Bt/PbY4qdoHLvutvUyVFveTxetm4jYj8aAre1odNS+lZl/WXz/l9HYpf0vwAeLHnDHGqYtPwL2pLGbexXw0QGD8TpaRPQCH8/MY9q5TkzwkiTVkLvoJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuqTET0R0S33fBIqgUTvCRJNWSClyaY4mqHVxf33749Iv5zRHw6Iv65mF5SXD1scw/8vIj4SXEP9YMi4rvFva3/Z1FnVkTcGRFLixuDfDsiXr6N9z0yIn4WEbdExP8rrgcvqSImeGniORpYl5lvzMx9gR8AX87Mg4rpnYFjBtR/NjPfBnyVxmU2TwX2BU6MiKlFndcBSzJzP+BJ4L8OfMOI2AP4c+CdxY2QVgCnV9ZCSSZ4aQK6DXhnRPx1RLw1M38FzIuImyLiNho3yfj9AfWvGDDfHcVNQJ4B7uO3N2l6MDP/qXh+GfCWLd7zUGAf4J+K24KeAPxuy1sm6SXbj1xFUp1k5t0R0QPMBz4bEdfR6JUfmJkPRsRiYNKAWTZfN/vFAc83T2/ehmx5zestp4PGfb7f34ImSCrBHrw0wUTEDOCpzLwMOAc4oHjp0eK4+HuHnHlo/zEiDiuevx+4YYvXbwTeHBGvKWJ4eUS8dhTvI6kke/DSxPMG4P9ExIvAc8B/oXE7ztto3Entn0exzNXACRHxNeAe4CsDX8zMRyLiROAbEbFTUfznwN2jaYCkkXk3OUljEhGzgKuKAXqSOoS76CVJqiF78JIk1ZA9eEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNXQ/wdqCgpquTUgpAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "event id =  6105\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfgAAAF3CAYAAABJzllyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xu4HWV96PHvT25JJQoJNA2E00SNVooYCBDA2wbl0oSA+kiPVFsQMNBDz6FFqdjTSmiPj7YHQbFUDeWSSil6VCo3BYzZWAqEJjTlYriVRAhJlUuQRCTcfuePNYl7J9l7z9p7zV5rTb6f51nPXjPrnVm/d8+s9VvvzDvzRmYiSZLq5TXtDkCSJLWeCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqocoTfERsFxH/HhHXF9NTI2JxRDwcEd+IiB2L+TsV048Ur0+pOjZJkupqNFrwZwLL+0z/NXBhZk4D1gKnFPNPAdZm5puAC4tykiRpGCpN8BExGZgN/H0xHcDhwLeKIguA9xfPjyumKV5/b1FekiQ1qeoW/BeBPwVeLaYnAM9m5svF9Cpgz+L5nsDjAMXrPy/KS5KkJm1f1Yoj4hjgZ5m5NCJ6Ns7eStEs8Vrf9c4F5gKMHTt2xl577dWCaBteffVVXvOaevQ7tC6dqS51qUs9wLp0qrrUpYp6PPTQQ09l5u5DFszMSh7A52i00FcC/wU8D/wj8BSwfVHmEOCm4vlNwCHF8+2LcjHYe8yYMSNbadGiRS1dXztZl85Ul7rUpR6Z1qVT1aUuVdQDWJIl8nBlP48y89OZOTkzpwAfBn6YmR8BFgEfKoqdCHy3eH5tMU3x+g+LikiSpCa14/jHp4CzIuIRGufYLy3mXwpMKOafBZzThtgkSaqFys7B95WZvUBv8fxR4KCtlHkBOH404pEkqe5GJcFLktSsl156iVWrVvHCCy+0O5Rhe/3rX8/y5cuHLrgVY8aMYfLkyeywww7DWt4EL0nqSKtWrWLcuHFMmTKFbr0tyrp16xg3blzTy2UmTz/9NKtWrWLq1KnDeu/uvwZBklRLL7zwAhMmTOja5D4SEcGECRNGdPTCBC9J6ljbYnLfaKR1N8FLkjSA7bbbjunTp7PPPvswZ84cnn322U2vLViwgGnTpjFt2jQWLFgwyFrawwQvSdIAxo4dy7Jly7jvvvsYP348F198MQDPPPMM5513HosXL+auu+7ivPPOY+3atW2Otj8TvCRJJRxyyCE88cQTANx0000cccQRjB8/nl133ZUjjjiC73//+22OsD8TvCSp40VU8yjrlVdeYeHChRx77LEAPPHEE/QdC2Xy5Mmbkn+nMMFLkjSAX/7yl0yfPp0JEybwzDPPcMQRRwCwtTupd1qHQBO8JKnjZVbzGMrGc/A/+clPePHFFzedg588eTKPP/74pnKrVq1ijz32qKr6w2KClyRpCK9//eu56KKLOP/883nppZc46qijuPnmm1m7di1r167l5ptv5qijjmp3mP2Y4CVJKmG//fbj7W9/O1dffTXjx4/nL/7iLzjwwAM58MAD+cxnPsP48ePbHWI/3qpWkqQBrF+/vt/0ddddt+n5ySefzMknnzzaIZVmC16SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJGkAgw0Xe/TRR7PLLrtwzDHHtDHCgZngJUkawEDDxQKcffbZfP3rX29jdIMzwUuSOltVQ8k1OaRc3+FiAd773vcybty4KmrcEiZ4SZKGsPlwsd3ABC9J0gAGGi62G5jgJUkawEDDxXYDE7wkqbNVNRh82UHh2XK42G5ggpckqYS+w8UCvOtd7+L4449n4cKFTJ48mZtuuqnNEfbncLGSJA1gsOFi/+Vf/mW0w2mKLXhJkmrIBC9JUg2Z4CVJqiETvCRJNWSClySphkzwkiTVkAlekqQmXXTRRbz1rW/lIx/5CBs2bOB973sf06dP5xvf+EZT6+nt7eX222+vJEavg5ckqUl/93d/x/e+9z2mTp3KnXfeyUsvvcSyZcuaXk9vby8777wzhx56aMtjNMFLkjSICy64gMsuuwyAU089lQceeIBHH32UY489lo9+9KNccsklPPnkk0yfPp1vf/vbXHLJJVx77bVsv/329PT0cNFFF/Hkk09y+umn89hjjwHwxS9+kT333JOvfvWrbLfddlx55ZV8+ctf5l3velfL4jbBS5K6wrzeeZx363mlyn58/48zf878fvPmXjeXS+6+ZNP0ue85l3k98wZdz9KlS7n88stZvHgxmcnMmTO58sor+f73v8+iRYvYbbfdmDlzJueffz7XX389zzzzDNdccw0PPPAAEcHjjz8OwJlnnsmf/Mmf8M53vpPHHnuMo446iuXLl3P66aez884788lPfrK5f0YJJnhJkgZw22238YEPfIDXvva1AHzwgx8c9Ba1r3vd6xgzZgynnnoqs2fP5j3veQ8AP/jBD/jxj3+8qdxzzz3HunXrKo3dBC9J0gCy5GhzG22//fbcddddLFy4kKuvvpovfelL3Hrrrbz66qvccccdjB07tqJItxLLqL2TJEkjMK9n3pCH1Aczf878LQ7bD+Xd7343J510Eueccw6ZyTXXXMPXv/51Lrjggq2WX79+Pc8//zyzZs3i4IMP5o1vfCMARx55JH/7t3/L2WefDcCyZcuYPn0648aN47nnnht2nQbjZXKSJA1g//3356STTuKggw5i5syZnHrqqey3334Dll+3bh3HHHMM++67L+95z3v43Oc+BzQuq1uyZAn77rsve++9N1/96lcBmDNnDtdccw3Tp09v+eh0tuAlSRrEWWedxVlnndVv3sqVKzc97+npoaenB4BJkyZx1113bXpt43n23XbbbavXyL/5zW/mnnvuaX3Q2IKXJKmWTPCSJNWQCV6SpBqqLMFHxJiIuCsi/iMi7o+I84r5V0TEiohYVjymF/MjIi6KiEci4p6I2L+q2CRJ3aHZy9TqZKR1r7KT3Qbg8MxcHxE7ALdFxPeK187OzG9tVv53gGnFYybwleKvJGkbNGbMGJ5++mkmTJhARLQ7nFGVmTz99NOMGTNm2OuoLMFn46fH+mJyh+Ix2M+R44B/KJa7MyJ2iYhJmbmmqhglSZ1r8uTJrFq1iieffLLdoQzbCy+8MOwkPWbMGCZPnjzs944qD39ExHbAUuBNwMWZ+amIuAI4hEYLfyFwTmZuiIjrgc9n5m3FsguBT2Xmks3WOReYCzBx4sQZV199dcviXb9+PTvvvHPL1tdO1qUz1aUudakHWJdOVZe6VFGPww47bGlmHjBkwcys/AHsAiwC9gEmAQHsBCwAPlOUuQF4Z59lFgIzBlvvjBkzspUWLVrU0vW1k3XpTHWpS13qkWldOlVd6lJFPYAlWSL3jkov+sx8FugFjs7MNUWMG4DLgYOKYquAvfosNhlYPRrxSZJUN1X2ot89InYpno8F3gc8EBGTinkBvB+4r1jkWuAPit70BwM/T8+/S5I0LFX2op8ELCjOw78G+GZmXh8RP4yI3Wkcpl8GnF6UvxGYBTwCPA98rMLYJEmqtSp70d8DbHFH/sw8fIDyCZxRVTySJG1LvJOdJEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEMmeEmSasgEL0lSDZngJUmqIRO8hnbLLfCGN0BEuccb3tBYRpLUNib4btFkkp35e7/XuiR72mmwYkX58itWNJaRJLWNCb4qrW71Nplkx65Z07ok20xyH8kykqSWMcFXpdWtXpOsJKkJJviqdEpCruJ8eebWH9K2rNmjds087NeiYagswUfEmIi4KyL+IyLuj4jzivlTI2JxRDwcEd+IiB2L+TsV048Ur0+pKrZaaFWS9Xy51BrNHrVrxooVcOSRo9v3Rl2vyhb8BuDwzHw7MB04OiIOBv4auDAzpwFrgVOK8qcAazPzTcCFRbn6GK1W79SpzS+zYsXgXxydzl7+A/N/M3Jl/4cdcEqspX1v1PUiR+HQakT8GnAb8IfADcBvZObLEXEIMC8zj4qIm4rnd0TE9sB/AbvnIAEecMABuWTJkpbFedIVJ7HgJwtKlf34/h9n/pz5/ebNvW4ul9x9Sanlz33Puczrmdcvgc45Aa5/S7lYv3bM15g7Y26/eTPmz+DuNXeXWv7aq2DOQ/3n7fEJWDOu3Psv+fgSZuwxo9+8OK/8j4EnznqCPcbtsWl69brV7HnBnqWXz3P77BZveANLX1jBASW/1yatg9Vf6D/vujfDsb9Xbvn9J+3P0rlL+82bv3Q+p11fLoBj3nwM151wXb9583rncd6t55Va/uP7f5z5Y47v12KcOwcumTHEgoVze2Feb/95Te1718Hvr57E2AUL4IgjgCb3vQ9fy5y3zOk3b48v7MGa9WtKLd/qfe9bN32L4+88vvTyOa//9NJJjN6+txqW9v/aYf4MOG3O1stv7pgH4bp/otEQ+NrX4Igjmt/3RvK91zvCfW+k33tt2PdO/M0TueKkK0qtv6yIWJqZBwxVbvuWvuuWQWwHLAXeBFwM/CfwbGa+XBRZBWz8Vt8TeBygSP4/ByYAT222zrnAXICJEyfS29vbsnhffPHF0mVXr1m9xXuvXrO69PIrV66kt7eXntJL9Pfggw/Su67/+69bt6708hvG70rjAMrwLFm6hHUPlX+/zd1+++3sttNum6af2vDUIKW31Pd/37NiBUwadihNW7du3Rbb/sHVD5Ze/umnnt5i+ZUrV5ZefvWa1fxy3omN1lqbjF2zhl+eeCKLr7oKaG7fu/feexm32S/JZj57rd73fvGLXwx7Xc16ccIEehd9q9+8e5+6He7/36MWAwArVmzafs3ueyP53hupkX7vtWPfe/HFF1uap5pRaYLPzFeA6RGxC3AN8NatFSv+bu1n0Bat98ycD8yHRgu+p6enNcECV6y8onTZPSbtwebvfdW6q6Dkd+6UKVMay0+dOqxDe295y1vomdH//cc9NA7Wl1t+pz/+Y/j0FcM+rHjAjAO2+CXLreWXP/TQQ7dowXNn+eV7DjusfOEWG/fQw/R8sv/7P3TkBDi03PIT7riDnv/Vf/neHij7a2+P629gbPty+yZj16zZ9BloZt9729veRs9bevrN23HpjlDye/aA005nxub1n1duWYBDjz+ePfp8R//nG38dfr/88iOx4447bvG9se7BdXB/yRXM2B+y/9Ejls6HkkeP+tq4/Xp7e+En5ZbZ4nvvllu46tYfwZubfvthGen33oj3vWF8721tm4+azByVB3AucDaNFvn2xbxDgJuK5zcBhxTPty/KxWDrnDFjRrbSokWLWreyzc+2D+TmmzOnTh3oDP2Wj6lTG8uMZl3KKluHsnVp5v9Sxf+6To+h/tdV/m9K7rNDqnLbTZ068PuW3b/aqartV9VntJlHl6viuxhYkjl03q2yF/3uRcudiBgLvA9YDiwCPlQUOxH4bvH82mKa4vUfFhWptyOOgEcfLb+7P/ropvOeXW+onsHDObowWCfDsv/r4XRUbKepU0e+35T93wzHUFdqdEIntsE6mnaDYvv1LlrU2u1X1We0mRjtGDpsVfainwQsioh7gH8DbsnM64FPAWdFxCM0zrFfWpS/FJhQzD8LOKfC2FSFdifGjR2HRuprX2t/XcpqVZ2beb/hGCwRVHl5mfpr935dxf7qpb4DqizBZ+Y9mblfZu6bmftk5l8W8x/NzIMy802ZeXxmbijmv1BMv6l4/dGqYhsVrWj5dJsqE2MrWqllNXNU5eabm6rzLydNaizTqgOYo31Ep4ptXEUrsMxjOPVod4IcqW74jA73Ul9twTvZqXXKJsbhJMXRbKU2o8lTLIuvuqq7T7GUPRRc5Y/adh2pGe2jJVWo6DPa0v/NcH+ElL3Xwzak0l700lZt/JIpaXFvb/t6oWr0jdYRr40/Vty/ttTkZ7SK9x5yu2xjyXo4bMFLkrpPt58uGQUmeElS9+mmzrBt4iF6SVL3aedphC5hC16SpBqyBS9pdNgpShpVtuAlSaohE7ykamyLN5KROogJXlI1tsUbyUgdxHPwkqphL2eprWzBS5JUQyZ4SZJqaMBD9BFxVonlf5GZnjSTJKnDDNaCPxvYGRg3yOMTVQcoSZKaN1gnu69vHMN9IBHx2hbHI0mSWmDAFnxm/ulQC5cpI0mSRt+Qnewi4syIeF00XBoRd0fEkaMRnCRJGp4yvehPzszngCOB3YGPAZ+vNCpJkjQiZRL8xhEiZgGXZ+Z/9JknSZI6UJkEvzQibqaR4G+KiHHAq9WGJUmSRqLMrWpPAaYDj2bm8xExgcZhekmS1KHKtOBvycy7M/NZgMx8Griw2rAkSdJIDHYnuzHArwG7RcSu/Oq8++uAPUYhNkmSNEyDHaI/DfhjGsl8Kb9K8M8BF1cclyRJGoEBE3xmfgn4UkT8z8z88ijGJEmSRmjITnaZ+eWIOBSY0rd8Zv5DhXFJkqQRGDLBR8TXgTcCy4BXitkJmOAlSepQZS6TOwDYOzOz6mAkSVJrlLlM7j7gN6oORJIktU6ZFvxuwI8j4i5gw8aZmXlsZVFJkqQRKZPg51UdhCRJaq0yvehvHY1AJElS65TpRb+ORq95gB2BHYBfZObrqgxMkiQNX5kW/Li+0xHxfuCgyiKSJEkjVqYXfT+Z+c/A4RXEIkmSWqTMIfoP9pl8DY3r4r0mXpKkDlamF/2cPs9fBlYCx1USjSRJaoky5+A/NhqBSJKk1hnyHHxETI6IayLiZxHx04j4dkRMHo3gJEnS8JTpZHc5cC2NceH3BK4r5kmSpA5VJsHvnpmXZ+bLxeMKYPeK45IkSSNQJsE/FREfjYjtisdHgaerDkySJA1fmQR/MvC7wH8Ba4APFfMkSVKHKtOL/jHAkeMkSeoiZXrRL4iIXfpM7xoRl5VYbq+IWBQRyyPi/og4s5g/LyKeiIhlxWNWn2U+HRGPRMSDEXHUcCslSdK2rsyNbvbNzGc3TmTm2ojYr8RyLwOfyMy7I2IcsDQibileuzAzz+9bOCL2Bj4M/DaNHvs/iIg3Z+YrpWoiSZI2KXMO/jURsevGiYgYT7lD+2sy8+7i+TpgOY3L7AZyHHB1Zm7IzBXAIziojSRJw1ImwX8BuD0i/ioi/hK4HfibZt4kIqYA+wGLi1l/FBH3RMRlfX487Ak83mexVQz+g0CSJA0gMoceN6Y4fH44EMDCzPxx6TeI2Bm4FfhsZn4nIiYCT9EYsOavgEmZeXJEXAzckZlXFstdCtyYmd/ebH1zgbkAEydOnHH11VeXDWVI69evZ+edd27Z+trJunSmutSlLvUA69Kp6lKXKupx2GGHLc3MA4YsmJmVPYAdgJuAswZ4fQpwX/H808Cn+7x2E3DIYOufMWNGttKiRYtaur52si6dqS51qUs9Mq1Lp6pLXaqoB7AkS+TgAQ/RR8TdQ/04GKxMRARwKbA8My/oM39Sn2IfAO4rnl8LfDgidoqIqcA04K6hYpAkSVsarLPcWyPinkFeD+D1g7z+DuD3gXsjYlkx78+AEyJiOo1D9CuB0wAy8/6I+CbwYxo98M9Ie9BLkjQsgyX43yqx/IAJODNvo/EjYHM3DrLMZ4HPlnhfSZI0iAETfGb+ZDQDkSRJrVPmMjlJktRlTPCSJNWQCV6SpBoa8Bx8RKyj0dN9i5eAzMzXVRaVJEkakcE62Y0bzUAkSVLrlBlNDoCI+HVgzMbpbIwTL0mSOlCZ8eCPjYiHgRU07im/EvhexXFJkqQRKNPJ7q+Ag4GHMnMq8F7gXyuNSpIkjUiZBP9SZj5NY1z412TmImB6xXFJkqQRKJPgny2GfP0R8I8R8SUa94qX1AVmz4aIco/Zs9sdraRWKZPgjwN+CfwJ8H3gP4E5VQYlqXVuHHD0h5GVldTZhuxFn5m/6DO5oMJYJFUot3ZXiz5ia0NDSepaZXrRfzAiHo6In0fEcxGxLiKeG43gJEnS8JS5Dv5vgDmZubzqYCRJUmuUOQf/U5O7JEndpUwLfklEfAP4Z2DDxpmZ+Z3KopIkSSNSJsG/DngeOLLPvARM8JIkdagyveg/NhqB1NHs2eUvO5o1C264odp4JEnbjiETfERctJXZPweWZOZ3Wx9SfXj9sSSpXcp0shtD49a0DxePfYHxwCkR8cUKY6uNzMEfkiS1Wplz8G8CDs/MlwEi4ivAzcARwL0VxiZJkoapTAt+T+C1faZfC+yRma/Qp1e9JKl9Zs+Gww7rccwBbVImwf8NsCwiLo+IK4B/B86PiNcCP6gyOEnbBgfEGTn7/GhzZXrRXxoRNwIHAQH8WWauLl4+u8rgpG7QzNUSM2e+jTvvrDaebmRyah3HHNBGA7bgI+K3ir/7A5OAx4HHgN8o5kmiuYSzePGE6gKpATukSq0zWAv+LGAu8IWtvJbA4ZVEJHUpW06SOsmALfjMnFv8PWwrD5O7upLnejUcVew37ouqWpnhYo+PiHHF8z+PiO9ExH7Vhya1nud6NRxV7Dfui6pamV70f5GZ6yLincBRwALgq9WGJVXLc70ajir2G/dFVaVMgn+l+Dsb+Epxe9odqwtJkiSNVJkE/0REfA34XeDGiNip5HKSJKlNyiTq3wVuAo7OzGdp3Ife698lSepgZW508zx9xn7PzDXAmiqDkiRJI+OhdknSFryMr/uZ4CVJW/Ayvu5XZrhYSdI2yjs0di9b8JKkrtTMaYRt8YeICV6S1JU8NTA4D9FLkrqad/zbOlvwkiTVkAlekqQaMsFLklRDJnhJkmrIBC9JUg2Z4CVJqqHKEnxE7BURiyJieUTcHxFnFvPHR8QtEfFw8XfXYn5ExEUR8UhE3BMR+1cVmyRJdVdlC/5l4BOZ+VbgYOCMiNgbOAdYmJnTgIXFNMDvANOKx1zgKxXGJklSrVWW4DNzTWbeXTxfBywH9gSOAxYUxRYA7y+eHwf8QzbcCewSEZOqik+SNHqava3sOee8rd0hd71RuZNdREwB9gMWAxOLMeXJzDUR8etFsT2Bx/sstqqY12/s+YiYS6OFz8SJE+nt7W1ZnOvXr2/p+qAHoMQ6y5Yrr0xdzjnnbSxePKH0OmfOfJrPf/7eEUbWvNZulx6g1dukinW2Utn3LVtu2/usdMd+09nlbryxZ4j19Ld48YQO/kyV1/rPShMys9IHsDOwFPhgMf3sZq+vLf7eALyzz/yFwIzB1j1jxoxspUWLFrV0fY0bKLau3KxZvyo71GPmzKdKv28zj3Zo5XZp9Tapap2tVEV8nf5ZaUaZunTDflOXclWts11a/VnJzASWZIn8W2kv+ojYAfg28I+Z+Z1i9k83Hnov/v6smL8K2KvP4pOB1VXG122aGVihmZZ5mdSuzlP2kKekbVOVvegDuBRYnpkX9HnpWuDE4vmJwHf7zP+Dojf9wcDPsziUr/46ORk3c55t9uz2xtrtmvnBN2tWdXFI6kxVnoN/B/D7wL0RsayY92fA54FvRsQpwGPA8cVrNwKzgEeA54GPVRibKtJM0nGox9Zo9486SZ2psgSfmbcBAx0gfO9WyidwRlXxaHQNlXQ8dCxJ1fJOdpIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvqSmzZ8Nhh/V4OaTU4UzwkprS7OWNXg4ptceo3IteUv2Uuf7eyyGl9rEFX2PewrQ7eQtaSa1ggq+hZm5L6i1MB9eOROstaCW1ggm+hm64ARYt6i01iMwNN7Q72s7UbOKsItG6/SSNhOfgpa0wcUrqdrbgJUmqIRO8Oloz11x7vbU2ambY4nPOeVu7w5UqYYJXR3P4WQ1HM/vC4sUTqgtEaiPPwasrOPys+iq7vd1vRo//y85jC15S1/AS0M7jNulctuAldQ2vbug8bpPOZQtekqQaMsFLklRDHqKXpC5gJzY1yxa8JHWwZjqmzZz5dHWBqOvYgpekDtZMJ7be3nuBnqpCUZexBS9JUg2Z4CVJqiEP0UuqnB3EpNFnC15SZbzLmdQ+tuAl9dPK1rZ3OZPaxxa8JMDLsaS6sQWvWvFc7/B5OVbruB+qE9iCVy14rledoNl9y31RVbIFr1rwXK86gfuhOokteEmSasgEL0lSDZngJUmqIRO8JEk1ZCe7DuKlNVJ7+NlTHdmC7wBe4iW1h5891Zkt+A7QTZfW2NJRndxwA/T29tLT09PuUKSWswWvUmzpSFJ3sQVfocx2R9A6VR1l8IiAJFXDFrzawoFNJKlatuDVFg5sIknVsgUvSVINmeAlSaohE7wkSTVUWYKPiMsi4mcRcV+fefMi4omIWFY8ZvV57dMR8UhEPBgRR1UVlyRJ24IqW/BXAEdvZf6FmTm9eNwIEBF7Ax8GfrtY5u8iYrsKY5MkqdYqS/CZ+SPgmZLFjwOuzswNmbkCeAQ4qKrYJEmqu3ZcJvdHEfEHwBLgE5m5FtgTuLNPmVXFvC1ExFxgLsDEiRPp7e1tWWDr169v6fraybp0oh6AEnUpW6596rNNtsW69ACdvh+Wfe+y5dqnnftXZIW3W4uIKcD1mblPMT0ReApI4K+ASZl5ckRcDNyRmVcW5S4FbszMbw+2/gMOOCCXLFnSsnjrdE9q69J5Nt61b6iPXNly7VSXbQLbXl26YT/shhjLqmL/ioilmXnAUOVGtRd9Zv40M1/JzFeBS/jVYfhVwF59ik4GVo9mbJIk1cmoJviImNRn8gPAxh721wIfjoidImIqMA24azRjk6RtScTgD3W/ys7BR8Q/0ThBsltErALOBXoiYjqNQ/QrgdMAMvP+iPgm8GPgZeCMzHylqtgkaVs1axbceGP5supelSX4zDxhK7MvHaT8Z4HPVhWP1ClsHamdqhoZUp3HO9lJo6TMVSHjAAAIWElEQVSZ1pAtJ0kjZYKXRskNN8CiRb1kMuTDVpakkTLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQ5Xdi16SpJFw3IaRsQUvSeoojtvQGiZ4SVJHcdyG1jDBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1ZAJXpKkGjLBS5JUQyZ4SZJqyAQvSVINmeAlSaohE7wkSTVkgpckqYZM8JIk1VBlCT4iLouIn0XEfX3mjY+IWyLi4eLvrsX8iIiLIuKRiLgnIvavKi5JkrYFVbbgrwCO3mzeOcDCzJwGLCymAX4HmFY85gJfqTAuSZJqr7IEn5k/Ap7ZbPZxwILi+QLg/X3m/0M23AnsEhGTqopNkqS6G+1z8BMzcw1A8ffXi/l7Ao/3KbeqmCdJkoZh+3YHUIitzMutFoyYS+MwPsD6iHiwhXHsBjzVwvW1k3XpTHWpS13qAdalU9WlLlXU4zfLFBrtBP/TiJiUmWuKQ/A/K+avAvbqU24ysHprK8jM+cD8KoKLiCWZeUAV6x5t1qUz1aUudakHWJdOVZe6tLMeo32I/lrgxOL5icB3+8z/g6I3/cHAzzceypckSc2rrAUfEf8E9AC7RcQq4Fzg88A3I+IU4DHg+KL4jcAs4BHgeeBjVcUlSdK2oLIEn5knDPDSe7dSNoEzqoqlCZUc+m8T69KZ6lKXutQDrEunqktd2laPaORWSZJUJ96qVpKkGjLBFyLi6Ih4sLhd7jlDL9G5ImJlRNwbEcsiYkm742lGM7c47mQD1GNeRDxRbJdlETGrnTGWFRF7RcSiiFgeEfdHxJnF/G7cLgPVpau2TUSMiYi7IuI/inqcV8yfGhGLi23yjYjYsd2xDmWQulwRESv6bJPp7Y61jIjYLiL+PSKuL6bbtk1M8DQ2CHAxjVvm7g2cEBF7tzeqETssM6d34WUmV1D+Fsed7Aq2rAfAhcV2mZ6ZN45yTMP1MvCJzHwrcDBwRvH56MbtMlBdoLu2zQbg8Mx8OzAdOLq4AumvadRjGrAWOKWNMZY1UF0Azu6zTZa1L8SmnAks7zPdtm1igm84CHgkMx/NzBeBq2ncPlejrMlbHHesAerRlTJzTWbeXTxfR+PLa0+6c7sMVJeuUtzWe30xuUPxSOBw4FvF/G7ZJgPVpetExGRgNvD3xXTQxm1igm+o261yE7g5IpYWd/7rdgPd4rgb/VExYuJl3XBIe3MRMQXYD1hMl2+XzeoCXbZtikPBy2jcMOwW4D+BZzPz5aJI13yPbV6XzNy4TT5bbJMLI2KnNoZY1heBPwVeLaYn0MZtYoJvKH2r3C7xjszcn8YphzMi4t3tDkhAY5TEN9I4DLkG+EJ7w2lOROwMfBv448x8rt3xjMRW6tJ12yYzX8nM6TTu/HkQ8NatFRvdqIZn87pExD7Ap4HfAg4ExgOfamOIQ4qIY4CfZebSvrO3UnTUtokJvqH0rXK7QWauLv7+DLiGxoe/m/104+iCm93iuKtk5k+LL7JXgUvoou0SETvQSIj/mJnfKWZ35XbZWl26edtk5rNAL40+BbtExMb7m3Td91ifuhxdnE7JzNwAXE7nb5N3AMdGxEoap3kPp9Gib9s2McE3/BswrejtuCPwYRq3z+06EfHaiBi38TlwJHDf4Et1vIFucdxVov8QyB+gS7ZLcR7xUmB5Zl7Q56Wu2y4D1aXbtk1E7B4RuxTPxwLvo9GfYBHwoaJYt2yTrdXlgT4/HoPGeeuO3iaZ+enMnJyZU2jkkB9m5kdo4zbxRjeF4rKYLwLbAZdl5mfbHNKwRMQbaLTaoXGnwqu6qS7R5xbHwE9p3OL4n4FvAv+N4hbHmdnRHdgGqEcPjUPACawETuuGMRci4p3AvwD38qtzi39G49x1t22XgepyAl20bSJiXxodtraj0VD7Zmb+ZfH5v5rGIe1/Bz5atIA71iB1+SGwO43D3MuA0/t0xutoEdEDfDIzj2nnNjHBS5JUQx6ilySphkzwkiTVkAlekqQaMsFLklRDJnhJkmrIBC+pMhHRGxHdNuCRVAsmeEmSasgEL21jirsd3lCMv31fRPz3iPhMRPxbMT2/uHvYxhb4hRHxo2IM9QMj4jvF2Nb/pygzJSIeiIgFxcAg34qIX9vK+x4ZEXdExN0R8f+K+8FLqogJXtr2HA2szsy3Z+Y+wPeBv83MA4vpscAxfcq/mJnvBr5K4zabZwD7ACdFxISizFuA+Zm5L/Ac8D/6vmFE7Ab8OfC+YiCkJcBZldVQkgle2gbdC7wvIv46It6VmT8HDouIxRFxL41BMn67T/lr+yx3fzEIyAbgUX41SNPjmfmvxfMrgXdu9p4HA3sD/1oMC3oi8Jstr5mkTbYfuoikOsnMhyJiBjAL+FxE3EyjVX5AZj4eEfOAMX0W2Xjf7Ff7PN84vfE7ZPN7Xm8+HTTG+T6hBVWQVIIteGkbExF7AM9n5pXA+cD+xUtPFefFPzTgwgP7bxFxSPH8BOC2zV6/E3hHRLypiOHXIuLNw3gfSSXZgpe2PW8D/m9EvAq8BPwhjeE476Uxktq/DWOdy4ETI+JrwMPAV/q+mJlPRsRJwD9FxE7F7D8HHhpOBSQNzdHkJI1IREwBri866EnqEB6ilySphmzBS5JUQ7bgJUmqIRO8JEk1ZIKXJKmGTPCSJNWQCV6SpBoywUuSVEP/H33nXFkBrHd5AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x432 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "for i, ev in enumerate(event[:5]):\n",
    "    print(\"event id = \", ev.r0.event_id)  \n",
    "    plot_waveform(ev, gain=1, nr_clus=200, pixel=0, ylim=[100, 400]) "
   ]
  },
  {
   "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.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}