{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# How to create epochs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "So, in your experiment, participants undergo a number of trials (events) and these events are possibly of different conditions. And you are wondering how can you locate these events on your signals and perhaps make them into epochs for future analysis?\n", "\n", "This example can be referenced by [citing the package](https://github.com/neuropsychology/NeuroKit#citation).\n", "\n", "This example shows how to use Neurokit to extract epochs from data based on events localisation. In case you have multiple data files for each subject, this example also shows you how to create a loop through the subject folders and put the files together in an epoch format for further analysis." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Load NeuroKit and other useful packages\n", "import neurokit2 as nk\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "plt.rcParams['figure.figsize'] = [15, 5] # Bigger images" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In this example, we will use a short segment of data which has ECG, EDA and respiration (RSP) signals." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## One signal with multiple event markings" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Retrieve ECG data from data folder (sampling rate= 1000 Hz)\n", "data = nk.data(\"bio_eventrelated_100hz\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Besides the signal channels, this data also has a fourth channel which consists of a string of 0 and 5. This is a binary marking of the Digital Input channel in BIOPAC. \n", "\n", "Let's visualize the event-marking channel below." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2AAAAEvCAYAAAAjA6I0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAagElEQVR4nO3da6xsZ3kf8P/jc45vYGOIjwFjnAMNpXWSNiQnDglVLxCFmwX90EpEJaVtgqtUrUholYKokFA/pa3SqGrV1AWakJiLyyVBiCRYTShFDXaOAYMvODGXmBMbfCxisEnx9e2HPTvebGb7zJ69Zs1aa34/aevsPXtm9J7zvs9Z81/zzrOqtRYAAABW74x1DwAAAGBTCGAAAAA9EcAAAAB6IoABAAD0RAADAADoiQAGAADQk8OreNILL7ywHTt2bBVPDQAAMHg33HDDPa21o7tvX0kAO3bsWE6cOLGKpwYAABi8qvqTebfbgggAANATAQwAAKAnAhgAAEBPBDAAAICeCGAAAAA9EcAAAAB6IoABAAD0ZKHrgFXVl5Lcl+SRJA+31o6vclAAAABTtJ8LMf+d1to9KxsJAADAxO0ngI3Wfd96KB/6zF3rHkanzj3zUF7+/U/P4UP730X6mZP35uY7v7GCUcEwPO/SC/JXnnb+vh/3rYceyYc/e1ceePjRFYwKpuEgx58bv3xvbrnL8YdxuvQp5+YF33PhUo+99pav5p77H+h4RGz7W3/5aC6+4Jx1D2NhiwawluQjVdWS/LfW2lW771BVVya5MkkuvfTS7kbYgXv//KG88f2fXfcwOvfU88/O85/9Xft+3OuvuTG3333/CkYEw/BD3/3kvO9nf2zfj/vff3Qqr7/mxhWMCKbl4gvOyQ8fe8q+H/fz7/l0vnDPN1cwIli9I4cqt/3bl+aMM2pfj7v7G9/Ka99xYkWjIkl+9R//8CQD2Ataa3dW1UVJrq2qz7XWPrbzDrNQdlWSHD9+vHU8zgN5+pPOzife+KJ1D6MzN568N//012/IQ48sd5b+wYcfzYu/96l5yyu+r+ORwfq9/ppP5/4HHl7qsQ/O3vl652t/JM++8IldDgsm4VN3/Fl+9upP5qEl3yV+4OFH87Lvf1refMX3djwyWK23ffwL+e//54tLPXZ7V8W/eflfzRV/7eIuh8XMBeceWfcQ9mWhANZau3P2591V9YEklyf52OM/ajgOHzojT3vS2eseRmdO/tmZB36Oc888PKl/E9h21uEzctBdHhedd5b6gDme8oSDH3/OOeL4w/g88ayDv8B/0jlHrH2SLNCGvqqeUFXnbX+f5CeS3LTqgQEAAEzNIu+APTXJB6pq+/7vbK39zkpHBQAAMEGnDWCttS8k+es9jAWgE21Qn0IFAHjM/nvIMhjLvshs8eqU6Zq9W78UlQGLUStsKmufLghgAAAAPRHARugAJ/gfe46DPwVMmAqBeQ7yDvNjz9HBQKBnnbz2sviZEcAAAAB6IoABk+NzjgDAUAlgwKTY4AEADJkAtoG06Ib5muKAhSgVNtUyxwn1wm4C2IipZwDGxEkOAAFspLRBhMdz0Nd4GlXBfLrwsqm6WLfWPtsEMAAAgJ4IYMCkePcKABgyAQwAAKAnAhgAAEBPBLARW7ablCZUTJ01Dqu17MXOlSZjt8waXrZemC4BDAAAoCcC2Ah10wZYpwKm6uBrW3XAfJ204lZgjFAnr72sfWYEMAAAgJ4IYAAAAD0RwIDJ8XFnAGCoBDBgUg6yx173RFiMWmFTLbP21Qu7CWAjpp4BGBMvRAEEsFHShQpWqxQIzKULL5uqi+OCQwvbBDAAAICeCGDA5DT7nACAgRLAgEmxwwMAGDIBbAN5dwDma1rbwEJUCptqmeOEemE3AWzMVDQAI+IkB4AANkqddOLpYBwwVeoD9qITHCxLB1C2CWAAAAA9EcCASXF2HQAYMgEMAACgJwIYAABATwSwDaQHFVO37JUWXKEBFuNyJmyqZZa+emE3AWzEtPMFYEy8DgUQwEapix4DGhUwVV20+VUfMF8XtaG+GCNrny4JYAAAAD0RwAAAAHoigAGT4/ORAMBQCWAbyIegmbKD7LFXG7AYpQKLUy/sJoCNmBeLAIyJwxbAPgJYVR2qqk9V1YdWOSBOr5NOPJ30UoRpUh8wXzeVob4YH8cFurSfd8Bel+TWVQ0EAABg6hYKYFV1SZKXJ3nraocDcHC25wIAQ7XoO2C/nOQXkjy6wrEAHJgLXQIAQ3baAFZVVyS5u7V2w2nud2VVnaiqE6dOnepsgHRPi26YT2XAghQLG2qZHRZ2ZbDbIu+AvSDJK6rqS0neneSFVfUbu+/UWruqtXa8tXb86NGjHQ8TAABg/E4bwFprb2ytXdJaO5bkVUl+r7X26pWPjD1td+I5yBkV27Rgb+oD5qsDFkdr6otx6qQDtcXPjOuAAZNjtwcAMFSH93Pn1tpHk3x0JSMB6IBrtQAAQ+YdMAAAgJ4IYAAAAD0RwDaQdqgwX1McsBCXM2FTLbf21QvfTgAboe0mOgcpZ414mKwSpGBVDn7oaD6lySh1sW6tfbYJYAAAAD0RwAAAAHoigAEAAPREAAMAAOiJALaBtCdgyirLr3G1AYvR54ZNtczaVy/sJoABAAD0RAAbsYO12tYMFfbiMg0w30FrozX1xTh1sW6tfbYJYAAAAD0RwIDpsd8eABgoAQyYlLLHAwAYMAFsA+nGA3tQG7AQxxE21TJLX7mwmwAGAADQEwFshLZ3WB2oB6JdWrAn2xhhvjpgB93WwXPAOnSxbq19tglgwOTY7gEADJUABkyK84sAwJAJYAAAAD0RwAAAAHoigG0kn5CBeZragIWoFDZVW+IaDC7bwG4C2Ah104kHpqlquQPktz1HR2OBqemiQagmo4yRtU+XBLARc0YFgDE56MkRgCkQwAAAAHoigAEAAPREAAMAAOiJALaBbMFnyirLd2hTG7AYn+ViUy2z8nXYZTcBDAAAoCcC2Ag91sZ0+TMqWqHC3tQHrEaLyzywuax9tglgAAAAPRHAAAAAeiKAAZNSVZppAACDJYBtIK9NYT61AYtRK2yqZU7wOSnIbgIYAABATwSwEeqiQ1vpxQN7Uh8wXyfHH21GGaEu1q2lzzYBbMS8pQ3AmDhuAQhgwMRUkuYTKgDAQAlgAAAAPRHAAAAAenLaAFZVZ1fV9VV1Y1XdXFVv6WNgrE6zCR/mUhqwGLXCxtKGng4cXuA+DyR5YWvt/qo6kuTjVfXbrbVPrHhs7KGLDm068TBZdfCDnfqA+XQIZVN1s/LVD1tOG8Da1tsl989+PDL7kuUHwCQAMCZ2YAAs+BmwqjpUVZ9OcneSa1tr1612WAAAANOzUABrrT3SWvuBJJckubyqvm/3farqyqo6UVUnTp061fU4AQAARm9fXRBba/cm+WiSl8z53VWtteOtteNHjx7taHgAAADTsUgXxKNVdcHs+3OS/HiSz616YKyOHfhMWaWWbsLhAs6wKLXCZlrmOOHYwm6LdEF8epJfq6pD2Qps17TWPrTaYQEAAEzPIl0QP5PkeT2MhQV10SJbI1TYm/qA+To5/igwRsjap0v7+gwYw6KbLwBj4rAFIIABAAD0RgADJsUWDwBgyASwDWTrIsynNmAxaoVNtczaVy/sJoABAAD0RAAboS52WJV9WrA35QFzddOFV4ExPp289urgOZgGAWzEXNgPgFFx2AIQwIBpqSTNhnsAYKAEMAAAgJ4IYAAAAD0RwDaQ7Vkwn8qAxagVNpW1TxcEsBHSwBD2pksbrE4XteEYxhh10T1aB2q2CWDA5DhDCQAMlQA2YnYSAjAmDlsAAhgAAEBvBDAAAICeCGAbyBYQpuxATQLs64WFKBU21TKdpNULuwlgwOQ42AEAQyWAjZI2wLBK6gPm6+YyDzA+1j5dEsBGzEl+AMZkme1bAFMjgAEAAPREAAMmpSpp3h8GAAZKANtEXpvCXEoDFuMkB5tqmZWvXthNAAMAAOiJADZC3XTi0YsH9qI6YL4uakOXUcbI2qdLAhgAAEBPBLAR084XvlOVCzHDUClNAAEMAACgNwIYALAv3mVmU1n7dEEA20D+74D5HFgB6JpjC7sJYCOkEw88noMv7lIgMFcnXXjVF2PUwbq19NkmgAGT42QjADBUAhgA0AtbsQAEMAAAgN4IYAAAAD0RwDaQCzgzZQf5kLPagMWoFDZVW2L1qxd2E8CAyZGjAIChEsBGqIsWvjqhwt7UB+zF8YfN1MklgKx+ZgSwEXOWH4AxWWb7FsDUCGAAAAA9EcCASbHBAwAYstMGsKp6ZlX9flXdWlU3V9Xr+hgYq2MDCNO33CpXG7AYHUPZWEssffXCbocXuM/DSf5la+2TVXVekhuq6trW2i0rHhsAAMCknPYdsNbaXa21T86+vy/JrUmeseqBsbdOOvHYpwV7Uh8wXye1ob4YIWufLu3rM2BVdSzJ85Jct4rBAByU8AQADNnCAayqnpjkfUl+rrX2jTm/v7KqTlTViVOnTnU5RvagnS8AY+KjMAALBrCqOpKt8HV1a+398+7TWruqtXa8tXb86NGjXY4RYF+8yAMAhmqRLoiV5G1Jbm2t/dLqh8SqeXEK86kNAB7PMocJhxZ2W+QdsBck+akkL6yqT8++XrbicQEAAEzOadvQt9Y+Hn1bBqWLJgOlUwETVR38d9XFc8AUddMITn0xPt0cW2DLvrogAoyB7R4AwFAJYCPm8yoAjInDFoAABgAA0BsBDAAAoCcC2AZyAWem7CD9ZVQGLMYWeDbVMmtfvbCbAAZMTnO0AwAGSgAbIa1QYcUUCMzVxSVMXAWFMXIJILokgAEAAPREABsxu6wAGBXHLQABDJgWGzwAgCETwDaQd86YumWXuOYdsBjddNlUy6199cK3E8AAAAB6IoCNUCdNdOzTgj1pVAXzOfywqax9uiSAAZOizS8AMGQC2IjZUQzAmPjsGIAABkyQXhoAwFAJYBvIa1MADsJJDjbVMmtfvbCbAAYAANATAWxDlV48sCfVAfN10eNGnxzGyNqnSwIYAABATwQwYHKaDfcAwEAJYCPmRSYAY+KwBSCAAQAA9EYA20TOQDJhB/mQs7PzsBi1wqZaZukrF3YTwEZIJx5YrVIgMFcXHXR14WWMrH26JIABk+NsIwAwVAIYAABATwSwEXOWH4AxcdwCEMCAibHHHgAYMgFsAzXnIJm6JZe42oDFqBQ21TLXYNU1lN0EMAAAgJ4IYCPURYtsm7Rgb+oD5nMZFDaWtU+HBDBgUhzgAIAhE8AAAAB6IoCNmQ91wlxKA4ZpmQYGAFMjgG0gxz+YT23AYgQpNtUyS1+9sJsABgAA0BMBbIS66DGgUQFTpT5g2JQXY9TJsaWD52AaBDAAAICeCGDA5NhvDwAM1WkDWFW9varurqqb+hgQi2t6vQEwIo5aAIu9A/arSV6y4nEAAABM3mkDWGvtY0m+1sNY6IkzkEzZQRpoqA1YjFqBxakXdvMZsBHqokNb6cUDe1IfMF8nHUK1GWWEqpsXX5CkwwBWVVdW1YmqOnHq1KmunhZg35xtBACGqrMA1lq7qrV2vLV2/OjRo109LQAAwGTYgggAANCTRdrQvyvJHyR5blWdrKqfXv2wWIRLHcF36mSfPrASjlsAyeHT3aG19pN9DIT+uEgtzKc0YEFqhQ21zHHCsYXdbEEEJsfBDgAYKgFshLpokW2XFuxNfcB8XWzxVV6MUSdXYLD6mRHAgElxeAMAhkwAAwAA6IkANmI+5gIAAOMigG0gwY2pa0uu8mUfB5tGrbCplln76oXdBDAAAICeCGAj1EWHNo0KmCyLG1amk05wapQR6uS1l7XPjAAGAADQEwEMmBwXYgYAhkoAAwAA6IkANmLO8gMwFs1BCyCJALaRHAOZsjpAmwC1AYtRK2yqpda+emEXAWyEOmmioxUP7El5wHzddOFVYIyPDtR0SQADJsfJRgBgqAQwAACAnghgAAAAPRHARqzZaAXfwee3YJg07gDYIoABAPsiS7GpNEGkCwIYMD2OdgDAQAlgY6QVKqyUNtkwXxe1YZswY9TN2rf42SKAAZPi8AYADJkABgAA0BMBDAAAoCcC2Igt09K36QPMBlj2Eg3qAxaz1PGn+2FA75Y5Tji0sJsABgAA0BMBbIR0oYK9dbG21QfM10l9HfwpoHeOLXRJAAMAAOiJAAYAANATAQyYHB94BgCGSgAbMa8xARgLXUYBtghgG8bxj6k7SJMa9QGLWfZSDzB2y6x89cJuAtgIddOFSise2IvqgPm6qA2d4NhUlj7bBDAAAICeCGDA5NjsAQAMlQAGTIrtTQDAkAlgAAAAPRHAxmyJlm22ZsHe1AcsZpmOoeqLKVhq7Vv87CKAAZPjekMAwFAJYCOkDTCsVikQmK+Ly6CoL0aoi3Vr6bNNAAMmxfENABiyhQJYVb2kqm6rqtur6g2rHhQAAMAUnTaAVdWhJP8lyUuTXJbkJ6vqslUPDAAAYGoOL3Cfy5Pc3lr7QpJU1buTvDLJLascGKd3y1335dpbvrqvxzzyqOYETN+jLfuujST5/Kn7VzAamJ5b7/rGEsefR1c0GujP//38PfniPd/c12NuvvPrKxoNY7VIAHtGki/v+Plkkh/ZfaequjLJlUly6aWXdjI45jv7yKGceeiMvOv6O/Ku6+9Y6jnOP3uRqYfxOf+cI0mS177jxFKPP++swz5HBns498zDOXKocvV1d+Tq6xx/2BxPmh1b3vxbNy/9HOeffaSr4TBydbp2zVX195O8uLX2M7OffyrJ5a21f7HXY44fP95OnFjuxQ+LufPe/5evffPBpR57RlWe+7TzcugMLzOZnkcebbntK/fl0SVb0V903lm56PyzOx4VTIfjD5uotZbPn7o/33pouXdyn3jW4Ry78Akdj4qhq6obWmvHd9++yGmok0meuePnS5Lc2dXAWM7FF5yTiy84Z93DgME5dEblsovPX/cwYLIcf9hEVZXvuei8dQ+DiVikC+IfJnlOVT2rqs5M8qokH1ztsAAAAKbntO+AtdYerqp/nuR3kxxK8vbW2vIbYAEAADbUQp+Eba19OMmHVzwWAACASVvoQswAAAAcnAAGAADQEwEMAACgJwIYAABATwQwAACAnghgAAAAPRHAAAAAelKtte6ftOpUkj/p/IkP7sIk96x7EHwbczJM5mV4zMnwmJNhMi/DY06Gybys3ne31o7uvnElAWyoqupEa+34usfBY8zJMJmX4TEnw2NOhsm8DI85GSbzsj62IAIAAPREAAMAAOjJpgWwq9Y9AL6DORkm8zI85mR4zMkwmZfhMSfDZF7WZKM+AwYAALBOm/YOGAAAwNpsRACrqpdU1W1VdXtVvWHd45myqnpmVf1+Vd1aVTdX1etmtz+lqq6tqj+e/fnkHY9542xubquqF++4/Yeq6rOz3/2nqqp1/J2moqoOVdWnqupDs5/NyZpV1QVV9d6q+tysZn7UvKxXVf387P+um6rqXVV1tjnpX1W9varurqqbdtzW2TxU1VlV9Z7Z7ddV1bE+/35jtce8/PvZ/2GfqaoPVNUFO35nXlZs3pzs+N2/qqpWVRfuuM2cDEFrbdJfSQ4l+XySZyc5M8mNSS5b97im+pXk6Ul+cPb9eUn+KMllSf5dkjfMbn9Dkl+cfX/ZbE7OSvKs2Vwdmv3u+iQ/mqSS/HaSl6777zfmrySvT/LOJB+a/WxO1j8nv5bkZ2bfn5nkAvOy1vl4RpIvJjln9vM1Sf6ROVnLXPzNJD+Y5KYdt3U2D0n+WZJfmX3/qiTvWfffeQxfe8zLTyQ5PPv+F83L+udkdvszk/xutq7Le6E5GdbXJrwDdnmS21trX2itPZjk3UleueYxTVZr7a7W2idn39+X5NZsvah5ZbZebGb259+dff/KJO9urT3QWvtiktuTXF5VT09yfmvtD9pW1b9jx2PYp6q6JMnLk7x1x83mZI2q6vxsHTjfliSttQdba/fGvKzb4STnVNXhJOcmuTPmpHettY8l+dqum7uch53P9d4kL/Iu5enNm5fW2kdaaw/PfvxEkktm35uXHuxRK0nyH5P8QpKdzR7MyUBsQgB7RpIv7/j55Ow2Vmz2NvXzklyX5KmttbuSrZCW5KLZ3faan2fMvt99O8v55Wz9R/zojtvMyXo9O8mpJP+jtraGvrWqnhDzsjattT9N8h+S3JHkriRfb619JOZkKLqch794zCw8fD3Jd61s5Jvjn2Tr3ZPEvKxNVb0iyZ+21m7c9StzMhCbEMDmpXStH1esqp6Y5H1Jfq619o3Hu+uc29rj3M4+VdUVSe5urd2w6EPm3GZOunc4W9tG/mtr7XlJvpmtbVV7MS8rNvtM0SuztTXn4iRPqKpXP95D5txmTvq3zDyYo45V1ZuSPJzk6u2b5tzNvKxYVZ2b5E1J3jzv13NuMydrsAkB7GS29sFuuyRbW0pYkao6kq3wdXVr7f2zm786e4s7sz/vnt2+1/yczGPbGHbezv69IMkrqupL2dqC+8Kq+o2Yk3U7meRka+262c/vzVYgMy/r8+NJvthaO9VaeyjJ+5P8WMzJUHQ5D3/xmNl20ydl/jYuFlBVr0lyRZJ/MNvClpiXdflL2TqJdOPsuH9Jkk9W1dNiTgZjEwLYHyZ5TlU9q6rOzNYHCD+45jFN1mxf8NuS3Npa+6Udv/pgktfMvn9Nkt/acfurZl12npXkOUmun20vua+qnj97zn+44zHsQ2vtja21S1prx7K1/n+vtfbqmJO1aq19JcmXq+q5s5telOSWmJd1uiPJ86vq3Nm/5Yuy9TlWczIMXc7Dzuf6e9n6f9FZ/SVU1UuS/Oskr2it/fmOX5mXNWitfba1dlFr7djsuH8yW83RvhJzMhzr6v7R51eSl2WrG9/nk7xp3eOZ8leSv5Gtt6Y/k+TTs6+XZWu/8P9K8sezP5+y4zFvms3NbdnRKSzJ8SQ3zX73nzO7cLivA83P385jXRDNyfrn4weSnJjVy28mebJ5WfucvCXJ52b/nr+erW5h5qT/eXhXtj6H91C2XkD+dJfzkOTsJP8zW00Irk/y7HX/ncfwtce83J6tzwhtH/N/xbysd052/f5LmXVBNCfD+dr+xwUAAGDFNmELIgAAwCAIYAAAAD0RwAAAAHoigAEAAPREAAMAAOiJAAYAANATAQwAAKAnAhgAAEBP/j/BHikuWN32fwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Visualize the event-marking channel\n", "plt.plot(data['Photosensor'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Depends on how you set up your experiment, the onset of the event can either be marked by signal going from 0 to 5 or vice versa. Specific to this data, the onsets of the events are marked where the signal in the event-marking channel goes from 5 to 0 and the offsets of the events are marked where the signal goes from 0 to 5.\n", "\n", "As shown in the above figure, there are four times the signal going from 5 to 0, corresponding to the 4 events (4 trials) in this data. \n", "\n", "There were 2 types (the condition) of images that were shown to the participant: “Negative” vs. “Neutral” in terms of emotion. Each condition had 2 trials. The following list is the condition order." ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "condition_list = [\"Negative\", \"Neutral\", \"Neutral\", \"Negative\"]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Before we can epoch the data, we have to locate the events and extract their related information. This can be done using Neurokit function [events_find()](https://neurokit2.readthedocs.io/en/latest/functions.html#neurokit2.events_find>)." ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "{'onset': array([ 1024, 4957, 9224, 12984]),\n", " 'duration': array([300, 300, 300, 300]),\n", " 'label': array(['1', '2', '3', '4'], dtype=') gives you a `dictionary` that contains the information of event onsets, event duration, event label and event condition. \n", "\n", "As stated, as the event onsets of this data are marked by event channel going from 5 to **0**, the `threshold_keep` is set to `below`. Depends on your data, you can customize the `arguments` in [events_find()](https://neurokit2.readthedocs.io/en/latest/functions.html#neurokit2.events_find>) to correctly locate the events. \n", "\n", "You can use the [events_plot()](https://neurokit2.readthedocs.io/en/latest/functions.html#neurokit2.events_plot>) function to plot the events that have been found together with your event channel to confirm that it is correct." ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA2AAAAEvCAYAAAAjA6I0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deXxV1b338e8iZCBAgJAwhlEEZBBEBkULyOBAEbXFOrWV68DTWrW19tapderV0tantr0+7S3VVm9tsdbh0jqPiBaEAg4ggyCDBCKZIISEEBLW80eOlqsxnLP32WedvfN5v17nRXKSs/qV3cXev73WXstYawUAAAAACF4b1wEAAAAAoLWgAAMAAACAFKEAAwAAAIAUoQADAAAAgBShAAMAAACAFKEAAwAAAIAUaRtEowUFBbZ///5BNJ0eNm5s+nPIELc5gIiriHW1rnQ1IDTot0AScK0ZCatWrSq31hZ++v1ACrD+/ftr5cqVQTSdHm66qenPH//YbQ4g4l6KdbXpdDUgNOi3QBJwrRkJxpjtzb4fxEbMY8eOtZEuwAAAAACgBcaYVdbasZ9+n2fAAAAAACBFKMC8+PKXm14AAvXol5teAMKDfgskAdeakRbIM2CRV1HhOgHQKtTS1YDQod8CScC1ZqQxAgYAAAAAKUIBBgAAAAApQgEGAAAAACkS1zNgxphtkqolNUpqaG45xVZl2jTXCYBWYQBdDQgd+i2QBFxrRlpc+4DFCrCx1tryeBplHzAAAAAArdnn7QMWyCqIlTX1WrjiwyCaTgtt2xidNbKnOmR7++trPGz19JoS1RxsSHIyIHqG9uioE/p28fz5reU1enMLq0kBqdQvP1cTBxV4/nzxnlq9vimue75AZB1f1EnDe3Xy/Pn3dlXp3eKqJCZCssRbQVhJLxhjrKTfWmsXfPoXjDHzJM2TpKweg3TTE2uSlzLNPPjobdrXo4M6LF3s6fNrdlbp2oVvJTUTEEUzHh2npZmH9Z8bvbcx/9n1ev693ckLBaBFMx4dJ2MO6+TNVsYYT2384qVNemxVcZKTAeHx4KO36UBOW2ndm57buOHxd7V2574kpkKyxFuAnWKt3WWM6SbpRWPMBmvtkiN/IVaULZCkUSecaJ+/KZpzVytr6rXvzwfVps774GF9w2FJ0q8uOkHj++cnKxoQOff97ZCq6/yNFNc3HNbQHh314L+NT1IqAC357dONKqmq89VGfcNhFXVpp8e+MTFJqYBw2fs/t+pg/UFfbdQ3HNaUIYWa/6Xjk5QKier5k+bfj6uKsNbuiv1Zaox5UtJ4SUs+7/czM4x6dMpJPGUItM0wSta9hK7tsyL79wQkQxuTnGm6WW3b0NeAFMloU5uUdjIz6LdovfZ5Gzz+jHaZGfSjNHTUZeiNMe2NMR0//lrS6ZLWBh0MAAAAAKImnhGw7pKejM3jbivpz9ba5wJNBQAAAAARdNQCzFq7RdKoFGQJjZePGa8zR3RXD9dBgIhrPHGfSstrJHlfBRFAamWeVKPi98sl8Ywz4NWa0V9QZU09F+ARFcgy9FH3uwlfUp9zhutE10GAiGuYXa6t75dJOtZzG0ff6RBAMmWfv0drX9oia4fK4yKIQKv34syvaUv5fl3pOggCcdRnwPC/cS4BUidZF2/0WyB1DD0O8I2bF9HGCJgHj/z5RvV/ob202tveDNZyTx6IR9atAzX+QB/pZtdJAMSr+vreOnNPvnS36yRAeF131zwdONQoXfe26ygIACNgAAAgrXCbEvCP+/3piwLMIUaXAQBoHudIwD+mMqYnCjAAkcYdQAAAkE4owAAAQNJx7wMAmsciHAkyxuipoV/QrON7sg8YELDDp1Sp5KNq+d4HjDkYQMpkTd6vbRtKxT5ggHerJ8xQafVB9gGLKAowDx4e80UdO3u4TnIdBIi4xjMr9OGGUklDXUcBEKecc6q0IXe7pOGuowChtWTG+dq0e7+udR0EgWAKogc5h+qUUXfAdQwg+g4atannnykgTGydUcYh+i3gR+bBOmXX17mOgYAwAubBg3+9Xf1f9b4PGID4ZN01QONqi6TbXScBEK/9N/fSjD1dpJ+5TgKE17d+dm3TPmA3nuE6CgLALSqvfDxdzIPJQAJ8dhj6G+CG9bEEqZ/PAmhCL0pfFGAAAAAAkCIUYAlK6lpqLMwGpARdDQghOi5ateR0AENHSksUYAAAAACQIizC4cFjI6fr3NG92AcMCFjjlD0qLtkn3/uAAUiZ7DP2afN7u0W/BbxbPmmWSqrq2AcsohgB8+CxkdO1eeYc1zGAyGuculfFY0pcxwCQgOwzqrV5ZLHrGECovTlptl4Yf5brGAgIBZgHXWqrlLW30nUMIPr2ZahtTaavJlhNDUitw1VtlF2byQpsgA/tq/eo4/69rmMgIExB9OA3//Nj9Xu9vfRF9gEDgpR1T1+Nrekp3e06CYB47b+jp07b01m613USILyu+OUNqq1vkG49x3UUBIARsASZJCwmww15ILWS0W8BxIn+BvhHP4o0CjAAAJBWuE8J+McU/PRFAeYQezMAANA8zpAAoooCDEBa4/4dAAAecScjLbEIhwcPnzBTXx5TpJ6ugwAR13hGpbYXV4n9hIDwyDm7ShvfLZG19FvAq9enzdHOvbUa7ToIAsEImAdPHTdJW6fNch0DiLzDp1Zp18jdvtvhBiCQOjlTa7T1OPbvA/x46+TTtXj0VNcxEBAKsAQZGfXcV6b2u3e5jgJEninPVE5VtusYABLQWNpW7ffluI4BhFrnio9UuLfUdQwEhCmIHtz71P9Vv2XtpS+zDxgQpLa/LNIJNd2ln7lOAiBe1T/uri9UdpR+5ToJEF6X/uZW1RxskH70FddREABGwAAAAAAgRSjAHLCs6wYkwF9/YRsUAACQTijAAABA0vm62ciNE8A3ulH6ogBzyLA0G5AShs4GhA79FkBUsQhHooz0u/Hn6YJxfdkHDAjY4XPKtWX7Hkn5rqMAiFO78/fovbd2ShrjOgoQWq/M/Ko+rPS/Dxi3MdITBZgHLw+aoImnDHMdA4g8O65au9uVu44BIAHZE2u1o5rlswE/1p44We/trHIdAwGhAPNgYEWxOn2YI2mA6yhAtO3MUvuyXNcpACSg4cNM5VW0dx0DCLVuu7ZpX2m16xgICAWYB3c/f5/6rWwvXTzNdRQg0jJ+01vH1xRKv/TeBquOAqm1/95umljZwXUMINQufOCupn3Afvp111EQABbhAAAASccWEADQPAowFzgpASnFQ8hA6tDfAKBlFGAJYlVcIHUMl3IAACBiKMAc4tISAIDPsrKcIwG/mHGVtliEw4P/nHihLpnAPmBA0A6fX6r3t1aKfcCA8Mj9WqXe+ecOSeNcRwFC6/lzr9DW8hrf+4AhPcVdgBljMiStlLTTWjsruEjp7x/9R+u0sce5jgFEnh1Vo7I2lf7a4A4gkFJZJx5QSVmF6xhAqL0/coLeLd7rux3DszNpKZEpiN+WtD6oIGEybPcW5W9a5zoGEH1bc5RX4n85a84/QOo0bM5W/u481zGAUCvavlGDdm12HQMBiWsEzBhTJOmLku6S9N1AE4XArS8vUL+320uXteqBQCBwGQ/01Mgaph8CYbL/vgKNr2QDdcCP8/77nqZ9wO690nUUBCDeEbBfSPq+pMMBZgkFbqQDAAAA8OqoBZgxZpakUmvtqqP83jxjzEpjzMqysrKkBQQAAACAqIhnBOwUSbONMdskPSJpqjHm4U//krV2gbV2rLV2bGFhYZJjRgtrAgCpwyIcAIAw4vQVXUctwKy1N1lri6y1/SVdKOkVa+1XA08GAAAAABHDPmAe/HTSpbp0Yj+d47MdlgYFWnb4ax9p3aZySRN8tWN4ehNImQ5XVmj1sg9lrfd+ay2rl6J1e/rCq/X+7mr92kcbjKClr4QKMGvtYkmLA0kSEsYYrS46TmeNYB8wIGhm6AFVHqxyHQNAArJGHFTpjj2uYwChtm3wKL3XLgn7gCUhC5IvkX3AEDOmeL26rW1xTRIAybChnfI/7OQ6BYAE1K/NVrfiLq5jAKHW//13NHzbWtcxEBCmIHrw/SUPqe977aWrzncdBYg088ceGlbDhRwQJvt/11VjKtu5jgGE2qxH7lP1wQZJ17iOggAwAgYgrfmdw26ZBQ8AANIIBRgAAAAApAgFmGfe76qzLxGQYjyFDKQco88A0DwKsARxHQekEB0OAABEDItweHDntHm6/NT++rLPdtjjBDiKyz/Smo2lkk52nQRAnPKuLteKf2yTn35rLfv3oXV78tJ/1/qSffq9jzYsU67SFgWYB+u6D1TFsUNdxwCib2CdqvZVu04BIAGZx9arcvM+1zGAUNvVf4g2G/bTiyqmIHpwyra31XvlP1zHAKLv7fYq2JzvqwluAAKpdXBlO/Xc1tV1DCDUBq95U2M2rfTdDrOt0hMjYB5cs/QR9d2QK333666jANH210INrc5znQJAAvb/sYtGVWZz8wPwYcYT9+ukgw2S/t11FASAETAAkccNQCB1uOMOAC2jAEsQJxYAAAAAXlGAAQAAAECKUIA5wOaUAAAAQOvEIhwe3HzG1Zo3aaAudB0EiLpv7tLqdR9J+oLnJqx4BgxIpU7Xl2npa1skneo6ChBaf73yB3p3Z5UWug6CQFCAebCla5H29h3oux0uCoGjKKrX/tJa383w7CaQOpl9G7Sva42vuR5Wln6LVq2sd3/tqK/01QbzrdIXUxA9mLZ5ufotfdl1DCD6VnRQz42FrlMASEDdP3LVZ3M31zGAUBu26jWdsm6p6xgICCNgCTIyunLFk+q7OVe6aZ7rOECkmUUFOra6o+sYABKw/9FOGl6R6ToGEGpTnvqjTjzQIOkHvtphIDk9MQIGAAAAAClCAQYg2pgEDwAIIU5f0UUBBgAAAAApQgHmka/VnbilAaSUYRY8kHKWkx0ANItFODy4btb1+uaUY/R110GAqLuuWCvWfCRpsuskAOLU5ZZSPfzKZkmTXEcBQmvht+7S6g/36EnXQRAICrAEGSOV5BWqpluvpLQF4POZwkYd6FTnOgaABLTt3qiaPPot4Mfegh4q3cdqolFFAebBrPVLdEzWB9KUa11HASLNvp6noi09XMcAkIDal9trwPqevtpg9iJau9FLn1f70mpJUz23QT9KXxRgHnz1rWfUZ3uudCcFGBCoZ7vomOpcX01YWbETCpA6NYvyNKSCywvAj5NffFTDDzRI+pHrKAgAi3AAAAAAEWR43iUtUYABiDzOP0DqMfsJ8I7TVrRRgAEAAABAilCAAQAAAECK8JSsB9889yZdM22QLvP4eaZlAPExNxZr2ds75WcVKACplX/nbj344vuSTnMdBQith757j1Zu26NnXAdBICjAPNiT20l1nfJdxwAiz3Rq1MHcQ77aYBleILUyOh/23W+B1q42r4v25ja6joGAUIB5MGfNSxqitdKUG3y2xCOWQEvsS53UbysnICBMap7toEFriny3w+ptaM3GLl6kzrv3S5rmOgoCQAHmwZw1L6nPzlxpvt8CDEBL7Mud1b+6ne92uI4DUqf22Y4aVJHha/SZgWu0duMW/02DDxyS9FPPbVh6UtpiEQ4AAAAASBEKsARxJx0AAABhwGVreqIAAwAAAIAUoQADEGnMgAcAAOmEAsyDueffrqd+fL/nz1vWxQbi0ua2D/X6RSt9t2OYhAGkTOHPduvF81e4jgGE2u9vvk/fvfRu1zEQEAowD+oyc9SQ439lNgAtMzlWjZmHXccAkIA2H/db7jUCnh3KbqeDWTmuYyAgR12G3hiTI2mJpOzY7z9mrb0t6GDpysjoq6uf1oi6ldJp/v4aWNADaJl9pouO2cZVHBAm1U921NA1/Xy3wykSrdnJz/9FhaXsAxZV8ewDdlDSVGvtfmNMpqQ3jDHPWmvfDDhb2pq14XX1Kc2V1GrrUCAlDr+epz7V2a5jAEjAgVfaq39FT9cxgFA7fukL6nfgkOsYCMhRCzDb9MDS/ti3mbEXt6QBhALPXALhQ7cF/KMfpa+4ngEzxmQYY96WVCrpRWvt8mBjAQAAAPCFubxpKa4CzFrbaK0dLalI0nhjzIhP/44xZp4xZqUxZmVZWVmycwKAZzxvCQAIF05cUZbQKojW2r2SFks6s5mfLbDWjrXWji0sLExSPAAAEEaWpxUAoFlHLcCMMYXGmM6xr9tJmi5pQ9DB0pUx0oUXz9eTP3/YdRQg8jJ/8qFe/Tr7CQFh0v2+j/Tcxa12nS4gKRbc8YCuuuLnrmMgIPGsgthT0kPGmAw1FWyPWmufCjZWtHFPEAAAAGid4lkF8V1JJ6QgS2hcufwJnbC3UJo633UUINIaH8/X4O3+5sFzwwNIrX1/ztOIdQNdxwBCbdLfHlLvMvYBi6p4RsDwKdM+WKE+e3J9t8PjlUDL7D87qHdVpusYABJwYGmuiiq6+W6HxXPQmg1dtUQ9DtS7joGAJLQIB/6FvRUAAPh8/s6TnGTRuiXj/gPXqumLAgwAAAAAUoQCLEHMiAAAAEAYGK5c0xLPgHlQ1zZbjdnZrmMA0Zdl1ZDZ6DoFgASYbKvGtvRbwI9DWdk6eIjiKaoowDyY+5U7dP2MwbrGdRAg4jLv3KHXlm2TdJbnNpgDD6RW9/+7W7/5+zpJM1xHAULrDz/4tV7ZUCp2wowmpiC6wAUhkFKG5dSAlKG/Af7RjaKNETAPrvnHQo3dXSBN+4XrKECkNS4s0PAdGa5jAEjA3j900qgNg7jXCPgw9a+/Vb+yGumW6a6jIAAUYB6csv0dFe1Lwj5g3N4AWnT4nVx1r2KgHgiTA6vaqWd5ge92OEWiNTtmzQrl17APWFRxZZMgiiYAAAAAXlGAAYg0pkEB4cPiOQCijAIMAAAAAFKEZ8A82NMuTwV5HVzHAKKvY6MO1h/y3QwTh4HUyejUqIM8uwL4Utuhk6rsQd/t8ORMemIEzINvnneznr7tPtcxgMjL/MFOvTHnLdcxACSg+91levW81bLMIwQ8+9P379UNF9/mOgYCQgEGAAAAAClCAZYgI+n7rz2oiQ/c47kNy7IAQFwa/9BNo14d7DoGgARU/qazTnxtiOsYQKid8fAv9a0XHnAdAwHhGTAPxuzcoJ61/vcBA9CywxvaqaDKZyNMgwJS6uDaHBWWd3EdAwi1vu+/o448SxlZjIA5xHORAAA0z3CWBBBRFGAAIo9VoAAAQLqgAAMAAEnnZ/IvE4cB/1iJNH3xDJgHJR0L1Kmwo+sYQOSZgkOqtXWuYwBIQNvCBtUerJPEfpmAV1Vdu2t3Bue/qKIAS5Ax0nVnf0/fmX6sjnMdBoi4rO+XaOkbWyX1dh0FQJy63V6uJYvekzTddRQgtP76nR/rhfd2a6bPdpiBn56YggggvfmcQcEEDABAOHEGiyoKMA9ufWmBJv/6Ls+fZ0ouEJ9D/9VNJ77IWDMQJuW/yNf4l4a5jgGE2hcf+Im++8xvXMdAQJiC6MGw0i0qPMQ+YEDQDm/JUZcq/3csmIIBpE79pizll+e5jgGEWs9tG5XDPmCRxQiYR8kYxGJpbABAVPmd7cE5EkBUUYABAAAAQIpQgCXIcEsOSBl6GwAAiBoKMA+25PfW3t79XccAIs/0rte+/BpfbbDoDZBamX0PJaHf0nHRulX06qftXf1twUIvSl8swuHBzWdeo29PO1YjXQcBIi7r2x9p+etbJPX11Q4j10DqdLuxUkv/Z62k7q6jAKG16Krb9Oyaj/Ql10EQCEbAAAAAgAji/mN6ogDz4O7n/lPT7v2h6xhA5NX/socmPDPCdQwACSidn6+Jz42UZQIU4Nk5v75Dtyy613UMBIQpiB4MrNypLnaP588ztR2Ij92Zpbyq9q5jAEjAoQ8zlVdJvwX8KNi1XYZ9wCKLETAAac3v/QruwgMAgHRCAeaQYZFtAACaxRkSQFRRgAGIPC7kAABAuqAA82Bdt4EqO2ao6xhA5LU5pk57uu1zHQNAArIH16uSfgv4UjJgqN7vcYzrGAgIi3B4cOf0ebp26iCNcR0EiLisb5Rq5WtbJA1wHQVAnAqv26MVT66TbE/PbfDkJlq7Zy+/QU+vKdHFPtpg0bf0xQiYB+ypAKQGz0kCAOAd59H0xAiYB/f+/R4NWdFROv1vrqMAkVb3k546pbSDdLf3NrgDCKTWR7d11aQdo6WbXScBwmvOvTfqpL0HpB/OcB0FAaAA86DHvnJ1yKjx/HmuB4H42PK2yq3OcR0DQAIaSum3gF95FbvVnX3AIuuoUxCNMX2MMa8aY9YbY94zxnw7FcEAIFmYNgykDt0NAFoWzwhYg6TrrbWrjTEdJa0yxrxorV0XcLbI46IQABBVvmd7cJIEEFFHHQGz1pZYa1fHvq6WtF5S76CDAQAAAEDUJLQKojGmv6QTJC0PIkxYvNV7qHYdN9p1DCDyMoYdUFnvva5jAEhAzsiDKuu9x3UMINR2DB2ld4uOcx0DAYm7ADPGdJD0uKTvWGs/s8OiMWaeMWalMWZlWVlZMjOmnZ9Nmas3LrvedQwg8rIvK9fqyRt8tcEqiEBqFVxVpVWTN/pqg36L1u6lr31Hv5pxma82LMu+pa24CjBjTKaaiq8/WWufaO53rLULrLVjrbVjCwsLk5kRAAAAACLhqItwGGOMpAckrbfW/jz4SOnv10/crYGvt5dOf8F1FCDSDtzZS1PK8qQf+22Jh/mBVCm5sUCn7RrDPmCADxfOv06Tqg5It57uqx3WsklP8YyAnSLpa5KmGmPejr1mBpwrrXU+sE/t9vFcChA0W52h7ANZrmMASEDjvjbKPpDFNELAh3bVe9X5QLXrGAjIUUfArLVviNvHSWU5KwEAAACtUkKrICKG8VwgJehpAAAgaijAAEQa481AOHEDBq0dE6aiiwLMg2X9Rmn76JNdxwAiL2N0rXb1L/fdDoPWQOrkjq1TST///RZozbaOOkkrBrLnbFRRgHlw36kX6c1LrnIdA4i87Esq9O4pm13HAJCAgsv36R36LeDLkgu+oQWTL3EdAwGhAPOIze0AAPh8fs6TnGEB/5jCmL6OugoiPuv3f7lVfV7Mlc5Y4joKEGm1t/TW9MouSdgHDECqFH+nUDM+Gifd5DoJEF6X3PENTd9XJ93mbx8wpCcKMA9yGurV9iCDh0Dg6tso41CGrybY9gFIrcMHjTIa/PVboLXLrK9T9qGDvtvhGej0RBXhAJeDAAAAQOtEAQYg8rgBCAAA0gUFmEMMCwMA0DzOkQCiigLMg1cGjdMH46e4jgFEXtsJNSo+ptR1DAAJ6HDqARUfU8oKbIAP74+drCWDJ7iOgYBQgHnwwIQv659zLnMdA4i8rPP3aO2ELa5jAEhA169W028Bn5ae9296aOIc1zEQEAowD5gWAaQGXQ0A0FoxiBxdLEPvwcMP36ief8+RzlzhOgoQaTXfK9KZe7pK810nARCv7d/opjPLTpJu9N4G20egtZt7y1x9sfqgdPsZntugF6UvRsAARB6j1kDq0N0AoGUUYAAAAEAkcUskHVGAOcDMCgBA1HGqA4DmUYA5ZLgrAQBAszhDAogqCjAPnj3uC9ow6UzXMYDIy5xcrW1DS3w9kM+IM5BaHacf0LahJa5jAKH23qln6vlhk1zHQEAowDz484mz9Nasi13HACIve3aVNozZ7rsdRpuB1Mk/f39S+i3Qmq2aeZH+Mu5s1zEQEAowD3IO1alt3QHXMYDIs3VGGYf4ZwoIk8P0W8C3tgcPKOdQnesYCAj7gHnwu7/cpp7P5Ujnsg8YEKSaW3prxp586R7XSQDEa8e3CzSjbLz0Q9dJgPC6+I5vakZ1nXT3Oa6jIADcogIAAEnHZsqAW3TB9EUBBiDt+TmJWBbDBgAAaYQCzAkuCAEAABAswxpUaYkCDEDkcQICAADpgkU4PHji+OmadGyh+vpsh4tCoGVZp+/T5nW7JeW7jgIgTp1n1WjzymJJx/pqx3CSRCv2zrRztGJrpca5DoJAMALmwZOjZmjNjPNcxwAiL/uMam0eWew6BoAEdD67ln4L+PTu9PO0aPTprmMgIBRgHnSprVJOVaXrGEDkHa5qo+zaTNcxACSgYW9Tv2UFNsC7dlV71Lm2ynUMBIQpiB786vG71eOlHOlC9gEDglR7Z0+dtqez7C+8t8FFIJBaxTd01WllJ0q3u04ChNec+ddpUnWd9LOvuI6CADACBgAAAAApQgEGIPJ4lh9IJf8djpFrtHbJOW3RkdIVBRgAAAAQQdx/TE8UYA5wZw8AAABonViEw4OFJ87UlMHd1M9nO0yLAlqWfXaVNq4pEfuAAeGRP2e/nnlzu6ShvtrhFInWbNXMC/SPzeUa7zoIAsEImAfPDZ+s9ZPOch0DiLzs0/Zr63Elsj6GjRlwBlKr0+kHtPW4EtcxgFBb94Wz9Ozwya5jICAUYB70qCpTxzJOLkDQGkvbqv2+HNcxACTg0EcZ9FvAp7yyEvWoKnMdAwFhCqIHP110j7ovzpEuZR8wIEg187vrC3vypP/0145hMhOQMsW35usLZaOl/3CdBAivc35+k06qrpN++TXXURAARsAAAAAAIEUowAAAQNL5WfHX8vQm4BurbqevoxZgxpjfG2NKjTFrUxEIAAAAAKIqnhGwByWdGXCOVoUbEkCcYo9u+ekzflZQBADAmSScvtjyKD0dtQCz1i6RVJmCLKHxh5O+pOXnzXUdA4i8dnP26r3xW/w3xAkISJnCr1Unp98CrdjyL83VH076kusYCAirIHqwePBJ6jPB7zbMrMwGHE3WxBrt2F/qOgaABORNqtOOj0olDffVDnfu0ZptnnCaFh/e7joGApK0RTiMMfOMMSuNMSvLyqK9b0H/imLlF291HQOIvMYdmcqraO86BoAEHNzWln4L+JRfvFX9K4pdx0BAkjYCZq1dIGmBJI0dOzbSD13c8fSv1P31HGneTNdRgEjbf283Tazs6DoGgATsvLuLJpaOlJ0f6UsBIFBn3Xe7xlTVSbrCdRQEgGXoAaQ9f8tZAwAApI94lqFfKGmZpCHGmGJjzOXBx0pzzEsHUoTOBgAAouWoUxCttRelIggABIUyDkfOeRoAABKjSURBVAgXdo8A/KMbpS+mIAIAAABAirAMvQe/PfUiTT2um7wuRM+dPSA+uV+t1Dsrd0ga7zoKgDh1u2KfHl/8gaSRrqMAobX0gv+jF9bt9n32Y8uj9EQB5sGygSdowKi+vtthjxOgZVknHlBJeYXrGAAS0HHCQZVsq/B9s5ELR7Rm206YqGV121zHQECYgujB0I8+UPct613HACKvYXOW8nfnyfqZyc6IM5BSBzZmKn93nusYQKh1+2C9hpR84DoGAsIImAc3vvBbdV+WI137ZddRgEjb//8KNb6SDV2BMNl1T2eNLx3mOgYQatN/N18jquokXe06SosOHTqk4uJi1dXVuY7iVE5OjoqKipSZmRnX71OAAYg8w3xfAACSrri4WB07dlT//v1b7bnWWquKigoVFxdrwIABcX2GKYgAAAAAElZXV6euXbu22uJLarrJ27Vr14RGASnAAAAAAHjSmouvjyX6d0AB5hHP9QMA8Pn8nCfZrgXwz7aSjnTXXXdp+PDhOv744zV69GgtX75cV1xxhdatW5f0/60OHTokpR2eAfPgl6fN1fRh3fVvHj/va0U3oBVpf3mFVi//UNZO8NwGvQ1IrZ5X79OjL22SNMp1FCC0Xrv0O3p6TYn+4jpImlu2bJmeeuoprV69WtnZ2SovL1d9fb3uv/9+19FaxAiYB+/0Ha7i405wHQOIvKyRdSot2uM6BoAEtB9VT78FfNo57AS9VeR/NdGozw4sKSlRQUGBsrOzJUkFBQXq1auXpkyZopUrV0qSHnjgAQ0ePFhTpkzRlVdeqauvblpZcu7cubr22ms1ceJEDRw4UI899pgkaf/+/Zo2bZrGjBmjkSNHatGiRUnPzQiYB6N3rFNRh3Jplr+OEfE+Afh2aG2OuhV38d0OfQ1Indp3spLSb+m4aM2K1r2tE4p3STrLdZS43fH397Ru176ktjmsV55uO3v45/789NNP15133qnBgwdr+vTpuuCCCzR58uRPfr5r1y796Ec/0urVq9WxY0dNnTpVo0b9a3S+pKREb7zxhjZs2KDZs2drzpw5ysnJ0ZNPPqm8vDyVl5frpJNO0uzZs5P6rBsFmAfffuUPKvxnjvTvl7iOAkRazf1dNaaynesYABKw6748jSkd4joGEGqTH7pXg6vqJF3vOkpa69Chg1atWqXXX39dr776qi644ALNnz//k5+vWLFCkydPVn5+viTp/PPP1/vvv//Jz88991y1adNGw4YN0+7duyU1PTt38803a8mSJWrTpo127typ3bt3q0ePHknLTQEGAAAAwJeWRqqClJGRoSlTpmjKlCkaOXKkHnrooU9+drSFSD6eunjk7/7pT39SWVmZVq1apczMTPXv3z/pG03zDBgAAEi61rICGwB3Nm7cqE2bNn3y/dtvv61+/fp98v348eP12muvac+ePWpoaNDjjz9+1DarqqrUrVs3ZWZm6tVXX9X27duTnpsRMACRxkUgAADRtH//fl1zzTXau3ev2rZtq0GDBmnBggWaM2eOJKl37966+eabNWHCBPXq1UvDhg1Tp06dWmzzkksu0dlnn62xY8dq9OjRGjp0aNJzU4ABiLyorwIFAEBrdOKJJ2rp0qWfeX/x4sWffH3xxRdr3rx5amho0HnnnafTTz9dkvTggw/+r8/s379fUtNKisuWLWv2f+/j3/GLKYge/OSMb+j5K25wHQOIvA5Xl2vFtORvpAggOL3/vcp3v2W/TLR2L/+fm3T3jHm+2qAXNbn99ts1evRojRgxQgMGDNC5557rOhIjYF5s7HGMjh/Yx/PnmREFxCdzUL0qNyd3SVsAwcodckiVb9FvAT9KjzlO63fluI4RCffcc4/rCJ/BCJgHJ21ZrQFvNz80mQimRQEtq1/VTj23dXUdA0ACqt/MTkq/5RSJ1qzfW0s1cetbvtuhH6UnRsA8mPf6QhW+nS398HLXUYBIq/ljF42qzPY1asyAM5BaH93fUaNKj6XvAT5MXPhfGlBVJ+lm11EQAEbAAAAAACBFKMAARB5TMIDUob8BQMsowAAAAACE0nPPPachQ4Zo0KBBmj9/vus4caEAAwAAABA6jY2N+ta3vqVnn31W69at08KFC7VuXfpvX0MB5sGds67VU1fd5joGEHl515dp6RlrXMcAkIA+P9xLvwV8euHaO/TDs652HSPtrVixQoMGDdLAgQOVlZWlCy+8UIsWLXId66hYBdGDbQV9VN67yPPnWRkKiE/bvoe0r2uNr01Z2XcPSK2c/o1N/dbP6qX0W7RylUUDtGWLv47goh89OOWz7w3/ijTuKulQrfSnmZ/9+ei5Ta/acunROf/7Z3MXt/y/t3PnTvXp86+9eYuKirR8+fIEU6ceI2AeTNn4pgavWOw6BhB5B5fmqs/mbq5jAEhA1WvZ9FvAp2PefFVTN/kvJEzEN521zVSZYfhvZgTMg68ve1yFa7Ol//iWz5bS//8ggEs1f+ms4ZVZvtsJwz/GQFSU/ncHDS8d6Lsdui1as3FP/EFFVXWS7nAdJSEtjVhl5rb889yCo494fVpRUZF27NjxyffFxcXq1atXYo04wAgYAAAAgNAZN26cNm3apK1bt6q+vl6PPPKIZs+e7TrWUTECBgAAACB02rZtq/vuu09nnHGGGhsbddlll2n48OGuYx0VBRgAAAgAK2kACN7MmTM1c2Yzq3ukMaYgAkh7vlZT4yIQABBCnL2iiwLMg1u+9H098Z0fu44BRF6XH+zW67Pe9t0Oz/IDqTPg7j1J6bdAa/bs93+q7519vesYCAhTED3Y3amb9hX2dB0DiLyMbo2qyatzHQNAArJ6HPbdb624cYLWrbqwp0ryan210dwS7UgPjIB5cMbaxRr++rOeP0+HAOJz4JX2GrCemx1AmFQ+l0O/BXwa8toz+uK6Ja5jICCMgHlw/j+fVuGGLOln3/PVDnucAC2rXdRJQyoyXccAkICyv7bXkN39fLdjGANDKzb6qYUqrKqTNN91FASAETCPGMQCUsdPd6OvAm7Q9wCf6ESRRQEGAAAAIJQuu+wydevWTSNGjHAdJW4UYB4wdRAIGfosACBMOG/Fbe7cuXruuedcx0gIBRgAAACAUJo0aZLy8/Ndx0hIXItwGGPOlPRLSRmS7rfWtuonAr/3lR/orBE9dIPrIEDEdbljtx586X1JU11HARCnY+7ZowWPv6vrdbLrKEBo/f0Hv9IDb2zRKtdBEjVlymff+8pXpKuukmprpZkzP/vzuXObXuXl0pw5//tnixcnP2MaOOoImDEmQ9L/k3SWpGGSLjLGDAs6WDrb276TavO6uI4BRF5G58M6mHvIdQwACcjsYum3gE8HOnXRntxOrmMgIPGMgI2XtNlau0WSjDGPSDpH0rogg6Wz2W+9oD7bc/XisYWePr92Z1WSEwHRVPtsBw1aU6RXNpSqXWaGpzYO1DcmORWAlpQvaqdBa4q09IMKbavwtpHsnpp6FXTITnIyIDyGv/CEvrS+VC+uG+u5jUONDlZRbGnEKje35Z8XFER2xOvT4inAekvaccT3xZImfPqXjDHzJM2TpL59+yYlXLo6792XdPBQoy4s/MxfQ0I6ZrMNG9CS+hfyNKgsQ9cufNNXO3k57CUGpMrep9prUEmRbvubv347sLB9khIB4TPq5UXqXFGjC/97uq928tpx/ktH8VQAza3D8pmS2lq7QNICSRo7dmykNy4Y0StP9Q2H9dQ1p3puo1O7THXLy0liKiB6Cjtmq312W13to69J0rHdOyQpEYCj6dSurUb16eS73w4ooABD69WjU47y2mX6utaUpCE9OiYpUfq66KKLtHjxYpWXl6uoqEh33HGHLr/8ctexWhRPAVYsqc8R3xdJ2hVMnHDIzGijzIw2GtGbublAsIxys9rS14BQMWqXSb8F/DCS2mdl0I/isHDhQtcREhbPMvT/lHSsMWaAMSZL0oWS/hZsLAAAAACInqOOgFlrG4wxV0t6Xk3L0P/eWvte4MkAAAAAIGLiWgXCWvuMpGcCzhIez/BXAaTCJXQ1IHTot0AScK0ZaSzD50VurusEQKuQSVcDQod+CyRBiK41rbUyprk1+1oPaxNbfzCeZ8Dwab/+ddMLQKD++eumF4DwoN8CSRCSa82cnBxVVFQkXIBEibVWFRUVysmJf3VzRsC8ePTRpj+vusptDiDi3ot1tXF0NSA06LdAEoTkWrOoqEjFxcUqKytzHcWpnJwcFRUVxf37FGAAAAAAEpaZmakBAwa4jhE6TEEEAAAAgBShAAMAAACAFKEAAwAAAIAUMUGsWmKMqZa0MekNIxUKJJW7DgHPOH7hxbELN45fuHH8wotjF25RP379rLWFn34zqEU4NlprxwbUNgJkjFnJsQsvjl94cezCjeMXbhy/8OLYhVtrPX5MQQQAAACAFKEAAwAAAIAUCaoAWxBQuwgexy7cOH7hxbELN45fuHH8wotjF26t8vgFsggHAAAAAOCzmIIIAAAAACmS1ALMGHOmMWajMWazMebGZLYNb4wxfYwxrxpj1htj3jPGfDv2fr4x5kVjzKbYn12O+MxNsWO40RhzxhHvn2iMWRP72a+MMcbFf1NrY4zJMMa8ZYx5KvY9xy4kjDGdjTGPGWM2xPrgyRy/8DDGXBf7d3OtMWahMSaH45e+jDG/N8aUGmPWHvFe0o6XMSbbGPOX2PvLjTH9U/nfF2Wfc+x+Fvu3811jzJPGmM5H/Ixjl0aaO35H/Ox7xhhrjCk44j2On7U2KS9JGZI+kDRQUpakdyQNS1b7vDwfl56SxsS+7ijpfUnDJP1U0o2x92+U9JPY18Nixy5b0oDYMc2I/WyFpJMlGUnPSjrL9X9fa3hJ+q6kP0t6KvY9xy4kL0kPSboi9nWWpM4cv3C8JPWWtFVSu9j3j0qay/FL35ekSZLGSFp7xHtJO16SrpL0X7GvL5T0F9f/zVF5fc6xO11S29jXP+HYpe+rueMXe7+PpOclbZdUwPH71yuZI2DjJW221m6x1tZLekTSOUlsHx5Ya0ustatjX1dLWq+mC4tz1HRxqNif58a+PkfSI9bag9barZI2SxpvjOkpKc9au8w29YD/PuIzCIgxpkjSFyXdf8TbHLsQMMbkqemk9IAkWWvrrbV7xfELk7aS2hlj2krKlbRLHL+0Za1dIqnyU28n83gd2dZjkqYxmpkczR07a+0L1tqG2LdvSiqKfc2xSzOf0/ck6V5J35d05IITHD8ldwpib0k7jvi+OPYe0kRsyPYEScsldbfWlkhNRZqkbrFf+7zj2Dv29affR7B+oaZ/vA4f8R7HLhwGSiqT9AfTNIX0fmNMe3H8QsFau1PSPZI+lFQiqcpa+4I4fmGTzOP1yWdihUGVpK6BJceRLlPTiIjEsQsFY8xsSTutte986kccPyW3AGuuEmWJxTRhjOkg6XFJ37HW7mvpV5t5z7bwPgJijJklqdRauyrejzTzHsfOnbZqmpLxG2vtCZJq1DQF6vNw/NJI7Fmhc9Q0RaaXpPbGmK+29JFm3uP4pS8vx4tj6YAx5hZJDZL+9PFbzfwaxy6NGGNyJd0i6dbmftzMe63u+CWzACtW01zPjxWpaboGHDPGZKqp+PqTtfaJ2Nu7Y8O9iv1ZGnv/845jsf41/H/k+wjOKZJmG2O2qWlK71RjzMPi2IVFsaRia+3y2PePqakg4/iFw3RJW621ZdbaQ5KekDRRHL+wSebx+uQzsWmpndT8tCskiTHmUkmzJF0Sm5YmcezC4Bg13bx6J3YNUyRptTGmhzh+kpJbgP1T0rHGmAHGmCw1PST3tyS2Dw9ic2QfkLTeWvvzI370N0mXxr6+VNKiI96/MLbizABJx0paEZu6UW2MOSnW5teP+AwCYK29yVpbZK3tr6b+9Iq19qvi2IWCtfYjSTuMMUNib02TtE4cv7D4UNJJxpjc2N/7NDU9Q8vxC5dkHq8j25qjpn+TQ30XPp0ZY86UdIOk2dba2iN+xLFLc9baNdbabtba/rFrmGI1LQj3kTh+TZK5ooekmWpaZe8DSbcks21eno/JqWoapn1X0tux10w1zZ19WdKm2J/5R3zmltgx3KgjVuuSNFbS2tjP7lNsI29eKTmOU/SvVRA5diF5SRotaWWs//2PpC4cv/C8JN0haUPs7/6Palq1i+OXpi9JC9X0vN4hNV3wXZ7M4yUpR9Jf1bRowApJA13/N0fl9TnHbrOanvv5+Nrlvzh26flq7vh96ufbFFsFkePX9Pr4PwwAAAAAELCkbsQMAAAAAPh8FGAAAAAAkCIUYAAAAACQIhRgAAAAAJAiFGAAAAAAkCIUYAAAAACQIhRgAAAAAJAiFGAAAAAAkCL/H9fdLph+FNOiAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot = nk.events_plot(events, data['Photosensor'])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Or you can visualize the events together with the all other signals." ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3AAAAExCAYAAAA9cUA3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzdeXgURfoH8G/lPkgCJOFMMAn3HSSgKEgEREAUT5AVlQVFUVbZVVfRdVf9uYo36+Kq6CrqIuCFBwIqKIKKYlBuCGeAEAgkIRcJOev3xztDAuSamZ7pOb6f5+mnSU9PTwGp6Xq7qt5SWmsQERERERGR+/MzuwBERERERETUNAzgiIiIiIiIPAQDOCIiIiIiIg/BAI6IiIiIiMhDMIAjIiIiIiLyEAzgiIiIiIiIPERAYycopd4CMBbAMa11L8uxxQC6Wk5pDiBfa51cx3szABQBqAJQqbVOMajcREREREREPkc1tg6cUuoSAMUA3rUGcGe9/gKAAq31E3W8lgEgRWudY0xxiYiIiIiIfFejPXBa6zVKqYS6XlNKKQDjAQwzslAxMTE6IaHOj/QO6emy79q14fOIyCG5lqoWzapG5DFYb4kMwLamV9iwYUOO1jr27OONBnCNGAIgW2u9u57XNYCvlVIawOta63lNuWhCQgLS0tIcLJobmzVL9k8/bW45iLzcSktVG8GqRuQxWG+JDMC2pldQSh2o83hjQygtb04AsPTsIZRKqVcB7NFav1DP+9pprbOUUq0AfAPgT1rrNfWcOw3ANADo0KFD/wMH6iwvERERERGR11NKbagrh4jdWSiVUgEArgWwuL5ztNZZlv0xAEsADGzg3Hla6xStdUps7Dk9hURERERERD7PkWUERgDYqbXOrOtFpVS4UirC+mcAIwFsdeDzvMd118lGRE71wXWyEZHnYL0lMgDbml6t0QBOKbUQwDoAXZVSmUqpqZaXbgSw8Kxz2ymllll+bA3gB6XUJgDrAXyptV5hXNE9WG6ubETkVCW5shGR52C9JTIA25perSlZKCfWc3xyHceyAIyx/HkfgL4Olo+IiIiIiIgsHBlCSURERERERC7EAI6IiIiIiMhDOLoOHNlj+HCzS0DkExJZ1Yg8DustkQHY1vRqTVoHztVSUlK0Vy/kTURERERE1ADD14EjIiJymbJiYPOHwJHNZpeEiIjIVBxCaYbRo2W/fLm55SDyRCV5QHUV0Cy20VMXWKraTaxqnq3iFDB/DHBkE6D8geveBHpda3apyElYb4kMwLamV2MPnBlKS2UjoqbTGlj5GPBsIvB8Z2D1bDnWgIpS2cjD/fKqBG+jnwXiBwIfTwV2LDW7VOQkrLdEBmBb06sxgPNW1dXA8oeA5zoD88cChUfMLhGRY7Z9AvzwEtDjaqD7lcDqp2Uj71ZVAfwyD0hKBS64A7jpI6BdP+CTacDRLWaXjoiIyOUYwHmrDW/JU+uWicDBn4EFNwCVZWaXish+a18EWvUArn8buOEdoOsYYM3znBPl7Q78BBRlAQNuk5+DmwE3vg+ENgfevxEoPmZu+YiIiFyMAZw30hr4+VUg/gJgylfADW8D2VuAn142u2RE9snZA2RvBVKmAH5+so17BQhtAXzzqNmlI2fav0bmvSUOrTkW0UaCuJJcYMmdjQ6lJSIi8iYM4MwwdqxsznJsO5C7B+gzAVBKhpt1HQP8NBcoK3Le5xI5y+6vZN/l8ppjYS2BwX8G9q0GMn6o821dxspGHmz/GqD9+UBI5JnH2yUDI/4B7F0FbP/MnLKRU7DeEhnA2W1NMhXXgfNGP8wBVv4DuG8XENFajh3eALwxDBj+D2DIX8wtH5Gt3rlKhsrd/fOZxytKgZfPB6LaA1O/kQcW5D0qTgFPtQMuvgcY8di5r1dVAm+kAidzgT+lAUHhLi4gERGR83AdOF9yeAPQIqEmeAOA9v2BTiOAn/4t6ykReYrKcpnH2XHYua8FhgKpDwKZvwI7vnB92ci5ju8EdBXQtm/dr/sHAKOfkzlyP3KIOBE5WekJzrslt8AAzgypqbI5y+HfJGA729AHgdI84Pf3nPfZREbL3gJUlUn6+LokTwJiu8sSA1WVZ7w0P1U28lDZ22Tfunf955w3SDKT/vgvoOCwa8pFTsV6S27pp38DzyTIMjbvXSMjBNyZs9uaZCoGcN6m+BhQmAm0O//c1+IHAh0GAev+c05Dl8htZW6Qfdw5IwiEfwAw7G9A3l7gt/kuKxa5QPY2ICBUsuk25LLHpadu1ROuKRcR+ZbDvwFfPwq0TwFSpgJ7vwW+mmV2qciHMYDzNod/k31dPXAAcPG9QMFBYMsHrisTkSMyfwUi2gKR7es/p9sVkqVw5RMc3uJNsrcCrboDfv4Nn9ciARh0N7B5kQwhJyIyitbAillAWDRw8yfA2BcliEt7C9j1ldmlI3e35SNg6V+AQ78aelkGcN7GurBtm3qGHHUZBbTuJWtqsReOPEHmr9L71lCCEqWAK14EKkuBrx5xXdnIuY6nSwDXFIP/AoTHAise5rICRGScHZ8Dh34Ghv8dCImSY6Nmy9D9pX9hdm+q348vAx9PBdL+C/x3BLDhHcMuzQDO2+Tulp6K4GZ1v64UMPSvct62T1xbNiJbncwFTuyXYSuNiekEXDxTepf3fe/8spFzlZ8Eio8CLZOadn5IpAylPfQzsJkjDIjIAFoDP7wERHcG+k2qOR4QBFz1MlB4GPj2n+aVj9xXZpqsU9vjauD+PcB5FwNLZwIH1hlyeQZwZhg/XjZnyN0DRHdq+JxuVwKtegKrZ7MXjtxblmVIcH3z38425C8ynO7L+4DKMvQcD/R0UlUjJ8vbJ/umBnAA0O9moG0y8PUjEvyTR2K9JbeRvRXI+h244I5zh3LHDwQGTAV+ea1mrrY7cWZbkxq39kUZdjtuLtAsFvjDYiAqHvh8BlBZ5vDlGcCZ4a67ZDOa1kDOHiCmc8Pn+fkBlz4sSR82LTS+HERGyUwDlJ80ypsiMBQY84L0MK98HAPuAgacXdWyt0sPnTVAIPeUu1f20R2b/h4/f2DcK0BpPvDJ7TVZ4o5uBdY8D3xxL/Dd0zI0k9xWnfWWyAzbP5N7UM9r6n59+D9kjvbSmUB1lWvL1hhH2pqVZcDOZcCeVUB1tbHl8gUnMoD0ZUD/yUBwhBwLjgCueEE6Wr5/xuGPCHD4CmS7khLZh4UZe92Tx4GyAunqb0y3K2RY2qongB7jZPgRkbs5vEHmGdQ3JLgunUcAA6cBP7+CiqB44MLpCAwDUHwc+PxPwK7lNee2T5Ev1HZNDBDJdezpgQOANr2A0bOlF/afrQH/IKCqXF4LbSnrOH3/jCQ9ueyJxhOkkMtVWG6RgQbfIolstuMLGfoWHlP36yGRwOVPAh9NATbMlx45d2FvW7O6Cnh/ArDvO/k5YQhw00dAYIix5fNmv74pgX/KWb8PnS8DuowG1r4A9P2DTP2wE3vgzDBmjGxGy9kt+6b8QigFjHlOgr5vnzS+LO6q+Djw01xg3StA0VGzS0MN0VoCuLh6Mqo2ZNRsIOlSLJjSCwtStsvv+L/7S+rnlCnA+PeAC6YDx3YA868Adiw1vvzkmLy9QHirmqeXthhwG3DDfKD5eRIApj4M3LcLeHA/cF860HcisG4u8GJ34FSB4UUnxywYIxuRqY7vAo7vBLpf1fB5Pa+VIGfVE+41dNvetuaPcyR4G3AbcOnfgIy1MuyPmqa6Gti0SDpKourInn25Zc7ksvscSrjFAM6b5FmGHLVs4pCj9ufL06L182SMt9kKMoFPpgHvXAns/c746+fuBV4ZKPNjvnoY+FeyNOjJPZ3YLwvP17ckRkP8/IGJCyWhT0kusOY5GVo8bTUw9iWgx1XSSzNttRxffJOh2aHIALn7bBs+ebae1wAzNwN3/wKkPghEtJbjEa2Ba14FLn0EKM4G/nu5JEwhIqptx2ey7z624fOUAkY/K9kov/XwtSgrTslD7s6XA2OeB4Y+IBl+t3wI7F5pduk8w9HN0jnS7Yq6X4/uKL8v+1bLEF07MYDzJoVZAFTD62WdbfjfJfX2l/ebO865+Bjw35HA5sXA/jXA/66VNVaMojXw5V+Aqgrg9u+A6T/Jk5EP/8i5UO6qsTUNGxMYKr0v8QOBP/0G3LYSaN3jzHNiuwCTlwEtEmXI3c5ljpWZjJO3r+kPo+wx9K/ANfOA4zvkwZG7zV8hInNt/xyIGwhEtmv83NY9gAvulAeB1nuXJ9rxuTw4vfDOmqV7Uh+Se+mKB4HKcnPL5wn2rpJ9x2H1nzPgNqBVD2DlY0B5iV0fwwDOmxQeBpq1kvS2TRUSJfNADqcB3892Xtkas/yv0lPyxxXArEyZm7T0z8alA9+3WrZhf5Oex9Y9gZs+lC+oBeOBU4XGfA4ZJ+t3ICAEiO3m2HX8g+WJV33ryAWFAdO+k7UTF08CNi127PPIcWXFliUEEp37OX0nyEOsnUuBj/5Yk/SEiHxb3n7pSenRyPDJ2lIfkrlyX9wjD4s90Y4vgIh2QGJqzbGAYJmWkLtHHoRTw/Z8K+2JZq3qP8fPHxj9jIw0WvW4XR/DAM6bFGY17UnR2freCLTpI8PM9q8xvlyNydwAbFsCDP4zcN4gmfMy6WMgOFIyyR3b4fhn/P4eENIcSPljzbGWSTIXKm8f8PFtzLTkbg7/Jr+X/oHO/6zQFsCtnwMdBgFLpkm2QvbImMfaK+7IEMqmGnIfMOJxGcoyL1UyVhKRb9tpmRfd/cqmvyckErjyX8DRLZL11tNUVUqG5k7DJFt5bV0uBzpcJG2pI5vNKZ8nKCuStUg7Dm/83MRLgIF3yDIUdrS9GcCZYfJk2YxWmGXb8EkrpYDJX8r6cYsmuX4+3G/vAIHhwKBak2RDIoHbLN3QCydKpbBXab4kqeh9gzxJqi1xCDDiMWD3V+b2QNKZqquAI5ukt9QByZNlaxLrg4OEIZJN7O3RQOERhz6f7GRvBkp7DZ4JXPO6DKd8cziQ8aNrPpfqZFO9JXKGrZ/IA8QWCba9r9sVQJ8bgTXPSjBkJlvbmlm/SSbz+oKPGxcAYTEy5NxTexidbf9aoLoS6NSEAA6Q9mfLjsCS6TaPBGs0gFNKvaWUOqaU2lrr2GNKqcNKqY2Wrc40N0qpUUqpdKXUHqXUQzaVzJs5LYA7bF8PHCAB08RFMm9owQ3G9Ho1RUWp9L71uOrcVPGxXYBbv5Au5q8ftf8ztn0CVJUB/W6q+/WL/iRfuN8/A+z+xv7PIePk7AIqTgLt+jl0GZsbgoEh8jt3xYtA1kbg1YuAQ786VAayQ64lo64z58Cdre+NwC2fSc//e1cDaW87lCGM7McAjkx1dIsEM8l/sO/9V7wgyzl9NMWSm8AktrY196wCoICk1LpfD2sJXDlHHnStm+tw8bzS3lXSIRF/YdPODwoDrnlN2u8LJ9p0z2lKD9x8AKPqOP6S1jrZsp0z818p5Q/gFQCjAfQAMFEp1ePs83xSTo5sRiorlnTY9gZwgAxXunmJ/Pn1ocDn9wBbPpKMkPmHnPPEJX0ZUFYojae6JF4iAdaGt4GfX7XvM35fIJNF61sMWin5UmrdC1h8M3DwZ/s+h4xj7QVu51gPXEmObDZRSrKzTlkuwzf/exmw7AGgJM+hspANcvbIPAxb1v8zQlIqMGO9DKVdOlN64/hQx+XsqrdERtkwX+ZO95lg3/uDmwET/gdUngI+uNW8xB+2tjX3fScPTcNa1n9Ot7GyrXyMy+/UZc8qGdllSy6K+IEyhejAD8BXjzQ5iGs0gNNarwFgT8tlIIA9Wut9WutyAIsAjLPjOt7n+utlM5L1KY89Qyhra90DuGMN0HU0sPF94OOp8jR6Ti9g9nnAm5cBq58Bjqc7XmZAhilEtAUSLqn/nOGPyTj0FQ8B69+w7frH0yVBS/If6k9iAUjP48RFgF+AJDXZv9a2zyFjZf0OBDWTYb0O+OB62ezSvj8wI03mTa6fB7ycDHzzd8mYSs6Vu9uhBU4dEtoCmPSJpHnO3QssuB5492pZcoTzZF3CoXpL5IjMNMmA3eu6hgOZxsR2AcbNBTLXy1x+M+ZU29LWrCyXUSfnXdTweUoBY+fIItVLZ5rbw+hu8vbJiLGmzH8727BHgcShwM+vWDKmVzb6FkfmwM1QSm22DLFsUcfr7QEcqvVzpuUYOUPhYdk70gNnFdkOGP+OZIOc/pNlSNkL0kt28hiw+ilZT+3dq4H0FfZ/zqlCebrd4+pzJ8zW5h8AXPdf6Y1bdr88+WmqTQsB5d+0J2nN44E710rmoHfHAT/MYYPNLFm/S49pQ78XrhASKevGTV0p8yF+/BfwYg8ZGrNvNYfYOYPW0gMX3dm8MvgHABfcAczcAgy5X34f37sGmNNbvn9y95pXNiIy3sFfgEU3AW+OkIfKl9mXGfAMPa8Bhj4IbP9U7hnunII/e4tMNWnKsj3NYoGbPwVKT0hbqTTf+eXzBPtWy77jpba/189PRsBddI88QFhwfaNz8O1tHb0KoCOAZABHALxQxzl1dXfU29pRSk1TSqUppdKOHz9uZ7F8mFE9cLUFhki6/cRLZM2KsS8C926yNGruk6GGCycAH062b3hZ+nL5wuh5dePnBgQDExfLEKcfXgK+/lvjjWetga0fS2VqKJ1rbS0TgdtXyXtW/gP43zXyJUWuU1UhcxDa1TPk1QzxA4DJS4Gp38iirls/lhvXKxfIoqeOJNmhM53MkYn0MSYGcFYhkcDwR4H7dgLXvyUjFH58Gfj3+cBrgyVzJR/yEHm2nN3AO2Ml82Tv6+V7vqlthsakzgIu/ZsEcQtvlHn/7ihzg+zjUpp2ftJQGbWUs0vuhe7693Klgz8DzVrbP3LIzx8Y+X+WRb6/k/bF1k/qP92ez9BaZ2utq7TW1QDegAyXPFsmgPhaP8cBqLevVWs9T2udorVOiY2NtadYvs0awEW0df5nNe8gayfdnw5cfC+w7VNgTh/psbJlHaXtn8k8l7i6fn3qEBQG3PQx0Ot64Kd/yxOKhp6EZ/4K5B+U820REgXc9BEw8knJRjd3oIyJZ9Yl1zi+U+YOOJjAxCniBwI3zAce2Cu/HxWlwNePAM8kAPPHAutekRshlyCwX45leLaZPXBnCwyVIVU3fQj8ZbvMyz2ZA3xwC/BcR+CzuyWozz9odkmJyBbVVVJ/q8qB274FrnsTiDLwQbhSwNAHgNHPSYKLl/vJUEV3czgNCG8FRMU3fq5V58vkwdaRjRLEnSpwXvk8wYF1QIcLG56u0xQX3AFM+14eInz0x3pPsyuAU0rVjhKuAVDXwjm/AuislEpUSgUBuBHA5/Z8HjVB4WFJ7xoY4rrPtC4CPvUboG0f6bH6Vx/gl3mNBzsnc4E930jvmy3D5PwD5At21GwZ8vCfC2VeUlH2ueeuf0OyAXWrM0lqw5SSRtqUr2RI6Rf3As8kAisfl4YbOc/pBCZuGMBZhcfI78fMzcCEBUDfiUD+AeCrh4E3h8lQu8/vkUXBCzLNLq1nObpF9m16mVuO+kS0keD93s3AuP/ICIXtX8gQqTm9gRe6A5/eBfz2LpehIHJ36cuAQ79IXY5rwvBBe10wTaaCFB2R9SbXPA9Uljnv82yVmQbEDbA9+Oh1HXDly8DhDcBrQ2TviwoygYKDkgDLCO2SgWmrgWvm1XuK0o0MQ1NKLQSQCiAGQDaAf1h+ToYMicwAcIfW+ohSqh2AN7XWYyzvHQNgDgB/AG9prf/ZlHKnpKTotLS0ppzqmRYvlv0EOzMc1WXBePliuNOk5Btay/CDNc/L05jIOODCO4F+kyQpwNlWzwZWPw3c9QvQqpt9n1mUDXzzKLDZ8u/ZZ4IkO+k8Up6GfzoduHim42PZtQZ2fA6segLI3SPHki4FUqYAXcdIUEnGWfpnYMvHwEMHHH6StdXyq9HLwKrWoILD0iDYs1KexpVZnkjGdJVhmG36yhdz22TbslT5kiXT5d/vgd1ml6Tpqiql4ZKxFtj1lSQusGp+HhB/gSQHaN9fhtcEhZlXVg/g8npLvuudK4G8/TI9xM/f+Z+Xf1B6/PavkdFMo2ZLO8LRXpu6NLWtWZIHPJsoI6uG3GffZ2X8AHx8O3DyOHDJAzI6y5UdCmbb8pEk/Zv2veHTP5RSG7TW54xtbTSAM4PXB3DO8Opg6fb/w2Jzy6G1rLv2w0vyJD0gBOgyShbR7jpaviCPpwNvDJOMOxPfd/wzs7fLopnblpx5PLQF8OftxjaW9n0P/PqmBHQAENoS6DRChhJ0G8uGmRHmpcqi2rd+YXZJHFNdBWRvlRvbrq+AY9vl5gYAQRGySHnzDtK4734lENrc3PK6i1cvll6uSR+bXRL7VVUAB9cBB36SJ9sHfpJ1Da3CYoBW3WWeX3RnGZrbvr9zGnFEVLcTGcC/+joWuNhDa5kTt2KWPHiP7gyMni1tCTPsXgksuA645XOZ22av4mPyAHbnUlkA/aJ7gP5/ND8ZmSt8eb8kzXvwgOEP9RnAuZNDluSc8TaMNW7MMwlAz2sl0Yi72L8W2LwI2LkMKM2TtLMdLgKObZPX7/zR2LHmZUXyFKTwMBAYJsPaIp00J7CsGNjyoUw03b1SGmf+wZL8JClVFiBOSmUvi60qTgFPxwGD7pLhuQ4qsFQ1W4b1O43WMlf1wI/yIOD4DnnIUVUudUNXS69u615AWLRMJo9o5xs3P6uyYmB2B1kTZ/ijZpfGOFWVsjTCsR3A0c3yEOt4OpBXaw5vaAtZ/LVVN6BtX1kDsXkHnwzq3Krekvf67V3g8z8Bd68HYru6/vMry4Ff3wC+ewooLwa6XiFrkHYcZky9b2pb87ungDXPSfAREun45+79TpZ9Or5TpkKkPix/J28erfTqxUB4LHDLp4ZfmgGcO0lNlf3q1cZcr7wEeKqt658iNVV1FfD7exJclRdLr9WIx2TenDeorpYnTls/lsb5yVpZVM+7WBpjbXpLspaWia4ZpuGpDq2XhbMnLJBsjw6anyr7yasdvpRzVJySyeNbP5EMVvkHpI6cpoCYLtJb07avZJltmSSN/OAI04rtNNYnwTcvkRu+t6sslyBu1wpZZy57O1CSi9MJmwNCZT2puAHSEIrpKv//4dGmFtvZ3L7eknf45A5JLHL/bnMflJwqBFY9LnOmyy0ZjTuPlDZD+/OBVj3sexjd1Lbm/LGSgMTIKThaA7+8Dnz3T6CsUBKkdB8LDJwGxHbzrgdTpfnSiXLpw8DQvxp++foCOC8Oh33I6TXg4swtR338/IH+k2XzRn5+QI+rZNNaGuEH1kkwl70VSHsbqLSk2A0IkflPnYbLl3Kb3kBUHIM6q8xfZd/UVMaeLjAESBgsGyAPA4qOyNCe7K1A0VHZH/xZhtzU1iJB5lO16iFBXnQn6dGOjPPcXrt93wF+gdIT5QsCgiQ4b9Vdeh0ByWx64Ccge5vMlzn4M/D7/2TotlWz1pJcKaaLNIZa9ZAhuAEhEuAZ8RSdyJtpLffo8y4yP5gIiZS1di97QnoFD/4sc6l3f11zTmhLoNsVkuWw43DjRheVl8gw7/63GnM9K6UkD0LKH+UB1U9zZX2ztLfk+6v3DbIGcPv+nnu/ssr8FYCW/xsXYgDnDaxZ7owcjkj2UUoa1i0SgOSJcqyqEjiySTJdHdkkc2O+q5XPxz9Ielda95SGWHQn+YJrcZ40yJSf7wR4h9YDUR1kDpQv8vOTehzVHki4+MzXyoqAEwckS2dOugzHy94uCT/OoGQ5kYjWkkCjRYL03EW0sfTe9WjaDbP4uMz1PJEhv4OABAlt+8pQaKMnqFeWS0KiLpf79lzSwFB5wNNpeM2xqgpJoJS9DcjbJ0kXsrfKEi5VdWSyi2gn/+8tkyQVdcsk+V5p3ZPBHREgD1oLDsk8LXcRFA5cOF226ip5gHNsh8yfPrxBvo9/fw+Akiy93ccBfSfIUGt77V0lD5i7jjbsr3GGgGCgxzjZCo9IjoRNi4B1c2ULbQEkDJG5d1EdZBqKf6BzyuIsB9cBfgFNWwTdQAzgvMHpHjgGcG7JP0DSE9dOUXyqUL6YM9ZKo+zIJllrri7KT+bUxQ2QTIbRnaXnzhuTXmSmSUIHOldwhNy0z06vX1YkN/r8QxJsFWfLXLsTGZab/heArrUunX+wNO7jB0qDvnVPmXMX2V7mjh7dDKyfV5Pd1S8AgJJr6FqLVncdA5x/q1wnrKXjf7/Ni2T4sbf21DvCP7Cmp+5sBYelp/bUCXmYV5In3y056TIss+is5Vcj2gHRHaXRFBUnf27ZEWjG9VfJh2T8KHvr6Ad34+cvUy5aJtYshVRdLUPuty2RedTfPSlb3EDp6eoyyvbv4u2fSe/eeS74d4hsCwy6W7aio8Dub2TURcaPNYnh/INkqPj5N9cEpsoPgKWXNCAYaB4v96uAYOeXuSkOrJMHm0HhLv1YBnDeoMAawLUztxzUdCGRQIcLZLOqqpT5MLl7ZSHrE/vlC7swUxpkmxcDm2pl7QxqJsOnWveQnpaYLtK70jLJM4ck5O6Vv2sHN3oi6gmCI2oCsbpYh/UWWH6PcvcCObtkTurv79V/3VY9ZUhPx2E1v095+6UneesnQPpyGeYDSEDY7nxJnxzbHYjpBDRPaPrv4bEdwNePSkPErExsnsraY1uf0hOWJ/k75f89b6/0dGecNd/F39IwCouRYF1Zev2VksA9INjSo9dG5mC2TAJaJLJHjzzTgR8lcIm1cxkjM/j5yQMz60PO3L3ygO7n/8iySYAEPR2Hyfz7yrKGg5zKMiB9hazH6+oEIxFtJEg7/2a5R53IkIzNB36UfALLmzCXrHkHuWckXCwPtuMGuH7pgtJ8WTZm0N2u/VwwgDPHfQYnGinMlAmi7vI0guzjHyCZsOrLhlVVIU+tjm6W7E45u6VBtuMLaaRZ+QVII6tlogylC2spw6ianydBfuueEvy527BMazDQZZRhlxzkhjl9XK72sN7aT5urq+SmmX9QEmcc2y6/Y216y1j+uoblWJ8I971R5k4cTgMO/yY9fZm/yvCY05/rJz17MV2A4Ei5nnVoTECwDBOurpJAYvc3QEgUcPWr5s9H8TahLQ0zFz4AACAASURBVGRr27fmmNbSc5uz2zIkc68kMcg/KMetQZv1XKWkobL5A6Ci5Mzrh8VIMBfbBQiOkt+zZq1kKGhAiOxbJEgZlH+TgnrWW3K6jLUy/80TH3ZaRXcEBs8EBs0A9q+WAOjgz8Bv78mInvgK+Z5dcqecG9NFgjtrAqxfXpOkKT2vNvNvId8v1nvL+TcDV82VrL0VltwBWuN0YqfSfPnOytklw8jTlwNbP5LXAkIkIG9/vgwZ7zBIlmpxZsKv9GVAdaVkEHUxZqH0Bu9dK2n6p602uyRkluLjMmTqeLr0tBRkyhdgSZ5s1sxWtQU1k0Zz4iXSwPILkC/5Zm1kHx7r2sb0G8PlieD0H1z3mWSs0718eyQoPLpVhkUWZkkmMkACg8qymmGdwVGyDl7qQ9IDRO7tZI48SMrbV7Md3SzBYFW5bPUJDJfgLiRKAvrI9rKPai8PmyIb6U0kMkL+IWBOL1lE+8LpZpfGeBWlMld63/dA1m8yRaM4u+b1yPbyQCV7K9D5clk/2FMfnFkfRGb9DuxbLdMwju8485ywGPmeOe8iGTYeEAyczJV70cnjNZnDdZX820R3ksAvYUjDD7qrKoF5Q+XB592/OO3fkFko3Ul6uuy7GrTuSOFh+YUj39UsVrb6xvMXHbXMkzoojevDv8lToxMZsmZL8dFz3xPaUnrs/IOkUWX9ErQO1wsKB0KaG/OldXiD9OaMmu34tWrJsVS1GBOW+PFJUXGydb6s4fOqqyXAq66UoTSe2njwReExsp09FxOoGQpVUSJLZFSWyuiA3L3SyCnOtvT2WtbD27n0zHmVAACFnNJBQEQbxPQIlGGasV1kSG94rIwo4O8LOeKAm89/c1RgqAQrp6KBdtdKW7OsCNj5pfRcZW2UoOWC6bL8lCfXJz9/eeAc3RHofb0c01pGKR3dIlvGD/Lds27uue8PCKkZbVJWDBR9UfNadCdg4B01GTq1PnOI5rdPSBB8/dum/BuyB84MRq8D91Qc0O8mYPQzxlyPfE91lfTUFRySuXcFmdLoKj4mPScnDgAlOTI3r7bAcMkYqPzkCbp1qIJSklHKmto8OMLS8PKTYRyhLeWLMDBMrv/BrfLZ9/wuT+cNwvWkiNxYVaUkWbF+15zIAEpyMf/+yUDlKUy+7XrLuni1BEdKL15UnJzf81r5fglqZsnk213+DMgyDURn+2yGTD34637PHkLZGKPbmp6uqkIC2ZI8WUsz2NLWqP07UFVpWQrqR+CXeUD2lprXQpoDt62S6S7fPgls+RDocyNw7etOLTZ74LzVqQIZHscMlOQIP/+aXrz259d9jjXIO/iTPL2rOCW9vxWlshVkSo8KII2uw7/LkITKU3U8ZT+bAm6Yb2jwRkRuzj9Ann6fPd/yWcv+r/vkO+fEfiB3nwRsBYekgVV4RIZB/Tinnu8XJUFeUDMZIh4UZtk3kyfrIVEyjCyijQSE5DsyfvD8+W9kO/9AeZDcUKZO/4CaHr3zbwH2rJI52oEhwM+vAv+7Fig/KQ+0z78FGP1s/ddyMgZwnu5EhuwdWQeEqCmsQV6Pcba/tyRPvvRKT0jjyzq8qqJUgr6k1PqzKBKR77I2uOpbY0lrmVNZfFSGhhcdlSGclWUyZDxvb00WVus8zLP5BdQMGQ+OkMDPP1CWV/DzB6Bq1kJUfhIMRsXJeyLayNBOd0sKRXU7cUAeCFxwp9klIU9Qe03Otn2BT6ZJe3vSR0C7fqYWjQGcp8vdK3vOgSN3Zm2ENY8H2vYxuzRE5C2Ukqfj1kyr9dFaeuoKD0smu7JCGSmQs1seLhVmyciC/APyYLToKFBd0bQyBIZLBr3wWOnVC4uW77noThIchjQHWpwn84k9eb6RN9i3WvZJqSYWgjxSz2sk26R/oFvUYwZwns4awLVMMrccRERE7kopWcagriGbdbHOl/GzNJN0dc1QzfwDMnqgMMuSGGqvJEkozavp9atLeKw8xbcurxAVL2neI1pLEGgN8iLaun49K1+xb7VkWq5vuR6ihrjRvFoGcGb429+Mu1buHpn/FhRm3DWJvMQlBlY1InINt6i31vkydWloDg0gwzdz90hQV1UpS7qUnpAhnoVZck7efqCqrP75wda5e8GRsnB6ULgkfQoMBSLaSS/fqXzpKVTKslaWhVLSLgiPlWA1ok3NGoy+rLoa2P890Okyt+hBcToj25rkdhjAmWHECOOulbtbJlsS0TmSDKxqROQaHl9vA4JrllsBAIys/9yCTOBUoQR2p/Il+CvKksCs4LAM+dz+mSVZVEn912mIX4AEfdEdZRh7WLT0/kW2lwQuMZ2l5y8g2L7re4qMtTJstrFlTryFkW1NcjsM4MywcaPsk5Mdu05VJZC9DUiZ6niZiLzQUUtVa+NgVSMi1/GpehsVB0QBaN2j8XO1lmGdubslK7CfP9D8PAnQlAKgAGhZzyr/oAzpPJEhW2GWJPA4slF6A8+hJKgLCpeev9AWZyZ1CYmSpC2hzWXKhif2YP32rvRodrvC7JK4hlFtTXJLDODMMHOm7B1dmyNnl6Rob8fKSVSXFZaqxnXgiDwH6209lAJCIuvPyGkVEgVENbC0UFmRJGzJ2yc9fSdzZHmGU4Uy7DP/ILDvu3PX/azNmoFTWbJvRraVoM8qPFYCwubx0rPXwpLkJSSy6X9fI+XtB7Z/CgycJsNQfYFRbU1ySwzgPFnmetmbnMqUiIiIPERwhGwNJT+zzqkrzJKAr6xQloPJ2ys9eIVZsg6trpbXc3bVep8Gio/LUjFnC2omnx3TRebnBYVLL2KzVhJ4hsdKRs+AEOOGdJYVAZ/eBfgFAhffa8w1iUzGAM6T7f4GiIzjEgJERERkHOsQyYZ68hqitWVZhkOyRt/J4zKEs6JEgr3SfGDXClkTtLyorgLUZONs3qFmnb6gZpKt0y9QhnoGBElPX3CkDCmtPbRTayDzV2DZA8DRLcB1b0ivIZEXcM8ALnsb8O8U6eaO7iiTa2M6S/d8QKiM927VXTJBBYWf2W3vKaorgZ/mAhk/SEpiwDLB+DygdS9ZLLmhzJJlxZIOt/f1njkWnYiIiLyTUtKr1qxVw+dpLYlFSvIkicuxHdJjln9Q1ufL2yc9fScOAOnLJMkLdN3XCo6SoZxB4ZKoJXubJIEJjgImvOc7c9/IJ7hnABcUDrTpDRQdkacmhUeAipP1n+8fLF3tQeHydCYoXFLmBlgW91RKgr4WiXI8Kl6ewrTu6fqx0NVVknUq/yDw9RZZ4LNtXyn/kU3Aji/kvM9nyPCG+IHAeRcDXUadmbp42QNAeTGQPMm15SciIiIyglJAeIxsgLR5GlKaL8lZyk/KvLbKMiA/Q5YIyNsnPXz5B6X92KoHkPwHYOAdQLNYp/9ViFxJaV3PkwwTpaSk6LS0tJoD1dUSwFnT6Fon3lZVyITbskLJyFheLFtZsfRwnciQSu7nL5N0qyvO/bDIOAniItvJGisBQdLjd8YWaFmrRVkWA/UDoKVLPypeuvNDomRMd0O9YfmHgC//Any7DGh3PnD78/JlVfs9leXA/jXA7q+APatkYnFVubzWMkmGEpQVA4fTgKRU4JbPHP73JvJWh36SffxF5paDiJqO9ZbIAD9ZKtJFrEieTCm1QWudcs5xjwjgjFBZXhME5u0Dio8BB9dJgFeSJ+Ozq8okKKwql/OrrFuFJWiDBHK6uv4FOFsmSY9gcIQEhQHBMoTg2E7J6qT8gMueAC6c3vRy7/8e2PWV9NyV5ErAl3QpMPjPMj6ciIiIiIi8CgM4o2ktAVVhFnDyGJC1UYLBk8csXfoHJQVvUbZM0I1oK4tHXngXsDtXrsGnIkROxSf5RJ6H9ZbIAOyB8wr1BXDuOQfOEygl65s0j5efu19Z93nWALn2MMkbUmXPtTmInGrVw7LnelJEnoP1lsgAD1sqEtuaXokBnLMxQyQRERERERnEr7ETlFJvKaWOKaW21jr2nFJqp1Jqs1JqiVKqeT3vzVBKbVFKbVRKufmYSCIiIiIiIvfWaAAHYD6AUWcd+wZAL611HwC7AMxq4P2Xaq2T6xq/SURERERERE3XaACntV4DIO+sY19rrSstP/4MIM4JZSMiIiIiIqJajJgDNwXA4npe0wC+VkppAK9rrecZ8Hmeb84cs0tA5BNGsaoReRzWWyIDsK3p1RwK4JRSjwCoBLCgnlMu1lpnKaVaAfhGKbXT0qNX17WmAZgGAB06dHCkWO4vOdnsEhD5hDasakQeh/WWyABsa3q1psyBq5NS6lYAYwHcpOtZTE5rnWXZHwOwBMDA+q6ntZ6ntU7RWqfExsbaWyzPsHKlbETkVPtWykZEnoP1lsgAbGt6Nbt64JRSowA8CGCo1rqknnPCAfhprYssfx4J4Am7S+pNnnxS9iNGmFsOIi+3xlLVkljViDwG6y2RAdjW9GpNWUZgIYB1ALoqpTKVUlMBzAUQARkWuVEp9Zrl3HZKqWWWt7YG8INSahOA9QC+1FqvcMrfgoiIiIiIyAc02gOntZ5Yx+H/1nNuFoAxlj/vA9DXodIRERERERHRaXbPgSMiIiIiIiLXYgBHRERERETkIYxYB45s9frrZpeAyCeMZVUj8jist0QGYFvTqzGAM0PXrmaXgMgnxLCqEXkc1lsiA7Ct6dU4hNIMX3whGxE5VfoXshGR52C9JTIA25pejT1wZnjhBdlfeaW55SDycussVa0rqxqRx2C9JTIA25pejT1wREREREREHoIBHBERERERkYdgAEdEREREROQhGMARERERERF5CKW1NrsM50hJSdFpaWlnHKuoqEBmZiZOnTplUqkMVFkp+wDPySETEhKCuLg4BAYGml0UoiYrOCT7qHhzy0FETcd6S2SAQ5aKFM+K5MmUUhu01ilnH/eYCCIzMxMRERFISEiAUsrs4vgUrTVyc3ORmZmJxMREs4tD1GRsABJ5HtZbIgMwcPNqHjOE8tSpU4iOjvaO4C0vTzYPoZRCdHS0d/R+kk/Zulg2IvIcrLdEBli8WDbySh7TAwfAO4I3ADh+XPYtW5pbDht4zb89+ZS0V2Xfa4K55SCipmO9JTLAq5aKNIEVyRt5TA8cERERERGRr2MAZwN/f38kJyef3mbPng1AEqw89NBD6Ny5M3r16oWBAwdi+fLlAIDi4mJMnz4dHTt2RL9+/dC/f3+88cEHZv41iIiIiIjIQ3nUEEqzhYaGYuPGjeccf/TRR3HkyBFs3boVwcHByM7Oxvfffw8AuO2225CUlITdu3fDz88Px48fx1vPPOPqohMRERERkRdgAOegkpISvPHGG9i/fz+Cg4MBAK1bt8b48eOxd+9erF+/Hu+//z78/KSzMzY2Fg/efruZRSYiIiIiIg/lkQHc419sw/asQkOv2aNdJP5xZc8GzyktLUVycvLpn2fNmoXu3bujQ4cOiIyMPOf8bdu2oW/fvqeDt9OSkgwpMxE1bPxHZpeAiGzFektkgI9YkbyZRwZwZqlrCOXmzZub/P5//vOf+PDDD3Hs2DFkZWUZXTwiOktYjNklICJbsd4SGSCGFcmbeWQA11hPmSt16tQJBw8eRFFRESIiIs54rUePHti0aROqq6vh5+eHRx55BI888giahYcDOTmsXEROtnG+7JMnm1kKIrIF6y2RAebPl/3kyWaWgpyEWSgdFBYWhqlTp+Kee+5BeXk5AODIkSP43//+h06dOiElJQV/+9vfUFVVBUAWJNdaA7m5ZhabyCdsnF/TGCQiz8B6S2SA+fNrgjjyOgzgbGCdA2fdHnroIQDAk08+idjYWPTo0QO9evXC1VdfjdjYWADAm2++idzcXHTq1An9+/fHiBEj8Mx995n51yAiIiIiIg/lkUMozWLtRTtbUFAQnn32WTz77LPnvBYZGYnXX3/9zIPp6c4oHhEREREReTn2wBEREREREXkIBnBEREREREQegkMozdCpk9klIPIJNy0zuwREZCvWWyIDLGNF8maN9sAppd5SSh1TSm2tdaylUuobpdRuy75FPe8dpZRKV0rtUUo9ZGTBPZq/v2xE5FSBYbIRkedgvSUyQFiYbOSVmjKEcj6AUWcdewjAKq11ZwCrLD+fQSnlD+AVAKMB9AAwUSnVw6HSeotjx2QjIqf69T+yEZHnYL0lMsB//iMbeaVGAzit9RoAeWcdHgfgHcuf3wFwdR1vHQhgj9Z6n9a6HMAiy/voxAnZiMiptn0gGxF5DtZbIgN88IFs5JXsTWLSWmt9BAAs+1Z1nNMewKFaP2dajnksf3//M9aBmz17NgAgNTUVXbt2RZ8+fdCtWzfMmDED+fn5Z7z3pZdeQkhICAoKCswoOhEREREReQFnZqFUdRzT9Z6s1DSlVJpSKu348eNOLJb9QkNDsXHjxtObdSFvAFiwYAE2b96MzZs3Izg4GOPGndnZuHDhQgwYMABLlixxdbGJiIiIiMhL2BvAZSul2gKAZV/XhK5MAPG1fo4DkFXfBbXW87TWKVrrlNjYWDuLZT7rot4HDx7Epk2bAAB79+5FcXExnnzySSxcuNDkEhIRERERkaeydxmBzwHcCmC2Zf9ZHef8CqCzUioRwGEANwL4g52fd6blDwFHtxhyqdPa9AZGz27wlNLSUiQnJ5/+edasWZgwYcI55/n7+6Nv377YuXMn+vbti4ULF2LixIkYMmQI0tPTcSw3F62io40tPxEREREReb1GAzil1EIAqQBilFKZAP4BCdw+UEpNBXAQwA2Wc9sBeFNrPUZrXamUmgHgKwD+AN7SWm9zzl/DNaxDKJtC65rRoosWLcKSJUvg5+eHa6+9Fh/+/jvuvvtuZxWTiCwmrza7BERkK9ZbIgOsXm12CciJGg3gtNYT63lpeB3nZgEYU+vnZQCMX0mwkZ4ys1VVVWHLli3o3r07Nm/ejN27d+Oyyy4DAJSXlyMpKYkBHBERERER2cyZSUx8UkVFBWbNmoX4+Hj06dMHCxcuxGOPPYaMjAxkZGQgKysLhw8exIFffzW7qERe76fnZSMiz8F6S2SA55+XjbwSAzgbWOfAWbfaWShvuukm9OnTB7169cLJkyfx2WcyLXDRokW45pprzrjONcOGYdHixS4tO5Ev2rVUNiLyHKy3RAZYulQ28kr2JjHxSVVVVXUeX93AOOP9+/efc+zFWbOMKhIREREREfkQ9sARERERERF5CAZwREREREREHoJDKM2glNklIPIJgaFml4CIbMV6S2SAUFYkb8YAzgxduphdAiKfcNNys0tARLZivSUywHJWJG/GIZREREREREQeggGcGbKyZCMip/r+/2QjIs/BektkgP/7P9nIKzGAs4G/vz+Sk5PRq1cv3HDDDSgpKUFGRgZ69epl03Xmv/MOsvbscVIpichq/yrZiMhzsN4SGWDVKtnIKzGAs0FoaCg2btyIrVu3IigoCK+99ppd15m/ZAmyjh0zuHTGq2/dOyIiIiIiMgcDODsNGTIEeyy9aFVVVbj99tvRs2dPjBw5EqWlpQCAjRs34sILL0SfPn1wzTXX4MSJE/joo4+Qtm0bbnrgASQnJ6O0tBSrVq1Cv3790Lt3b0yZMgVlZWUAgIceegg9evRAnz59cP/99wMAjh8/juuuuw4DBgzAgAED8OOPPwIAHnvsMUyZMgWpqalISkrCyy+/DAA4efIkrrjiCvTt2xe9evXC4sWLAaDez0xISMATTzyBwYMH48MPP3TdPygRERERETXKI7NQPrP+GezM22noNbu17IYHBz7YpHMrKyuxfPlyjBo1CgCwe/duLFy4EG+88QbGjx+Pjz/+GJMmTcItt9yCf//73xg6dCj+/ve/4/HHH8ecOXMw99ln8fxf/4qU66/HqVOnMHnyZKxatQpdunTBLbfcgldffRW33HILlixZgp07d0Iphfz8fADAvffeiz//+c8YPHgwDh48iMsvvxw7duwAAOzcuRPfffcdioqK0LVrV0yfPh0rVqxAu3bt8OWXXwIACgoK6v3MmTNnAgBCQkLwww8/GPrvS0REREREjmMPnA1KS0uRnJyMlJQUdOjQAVOnTgUAJCYmIjk5GQDQv39/ZGRkoKCgAPn5+Rg6dCgA4NZbb8WaNWvkQkoB/v4AgPT0dCQmJqKLZWkB63mRkZEICQnBbbfdhk8++QRhYWEAgJUrV2LGjBlITk7GVVddhcLCQhQVFQEArrjiCgQHByMmJgatWrVCdnY2evfujZUrV+LBBx/E2rVrERUVVe9nWk2YMMHJ/5JErhEWLRsReQ7WWyIDREfLRl7JI3vgmtpTZjTrHLizBQcHn/6zv7//6SGUDVwIiI8HAGit6zwlICAA69evx6pVq7Bo0SLMnTsX3377Laqrq7Fu3TqE1rFA49nlqKysRJcuXbBhwwYsW7YMs2bNwsiRI3HVVVc1WLzw8PCGy0/kIcZ/bHYJiMhWrLdEBviYFcmbsQfOSaKiotCiRQusXbsWAPDee++d7o2LiIg43WvWrVs3ZGRknJ5PZz2vuLgYBQUFGDNmDObMmXM6cBw5ciTmzp17+nPqCihry8rKQlhYGCZNmoT7778fv/32W72fSURERERE7s0je+A8xTvvvIM777wTJSUlSEpKwttvvw0AmHzllbjzttsQGhGBdevW4e2338YNN9yAyspKDBgwAHfeeSfy8vIwbtw4nDp1ClprvPTSSwCAl19+GXfffTf69OmDyspKXHLJJQ1mw9yyZQseeOAB+Pn5ITAwEK+++ipCQkLq/Ewib7NyluxHPG1uOYio6VhviQwwy1KRnmZF8kaqviF8ZkpJSdFpaWlnHNuxYwe6d+9uUokMlp4u+65dzS2Hjbzq/4B8wvxU2U9ebWYpiMgWrLdEBkhNlf3q1WaWghyklNqgtU45+ziHUBIREREREXkIBnBEREREREQeggEcERERERGRh2ASEzMEBppdAiKfEBlndgmIyFast0QGiGNF8mYM4MyQlGR2CYh8wrX/M7sERGQr1lsiA/yPFcmbcQglERERERGRh2APnA38/f3Ru3dvVFZWIjExEe+99x6aN2+O6upqzJw5E99++y2UUggJCcEHH3yAxMREJCQkICIiAn5+fmjdujXeffddtCkvlwt26GDuX4jIy62YKftRc8wtBxE1HestkQFmWirSHFYkb8QeOBuEhoZi48aN2Lp1K1q2bIlXXnkFALB48WJkZWVh8+bN2LJlC5YsWYLmzZufft93332HTZs2ISUlBU899RRQWiobETnV0Y2yEZHnYL0lMsDGjbKRV2IAZ6dBgwbh8OHDAIAjR46gbdu28POTf864uDi0aNHinPdccskl2LNnj0vLSURERERE3sPuIZRKqa4AFtc6lATg71rrObXOSQXwGYD9lkOfaK2fsPczrY4+9RTKdux09DJnCO7eDW0efrhJ51ZVVWHVqlWYOnUqAGD8+PEYPHgw1q5di+HDh2PSpEno16/fOe9bunQpevfubWi5iYiIiIjId9jdA6e1TtdaJ2utkwH0B1ACYEkdp661nmdE8Gam0tJSJCcnIzo6Gnl5ebjssssASI9beno6nn76afj5+WH48OFYtWrV6fddeumlSE5ORmFhIWbNmmVW8YmIiIiIyMMZlcRkOIC9WusDBl2vQU3tKTOadQ5cQUEBxo4di1deeQX33HMPACA4OBijR4/G6NGj0bp1a3z66acYPnw4AJkDFxMTU3Oh/Hwzik/kc6K7mF0CIrIV6y2RAbqwInkzowK4GwEsrOe1QUqpTQCyANyvtd5m0GeaJioqCi+//DLGjRuH6dOnY8uWLWjTpg3atWuH6upqbN68GX369Kn/AgkJLisrkS+7cp7ZJSAiW7HeEhlgHiuSN3M4gFNKBQG4CkBdYwN/A3Ce1rpYKTUGwKcAOtdznWkApgFABw9Ir9+vXz/07dsXixYtQmxsLG6//XaUlZUBAAYOHIgZM2aYXEIiIiIiIvI2Smvt2AWUGgfgbq31yCacmwEgRWud09B5KSkpOi0t7YxjO3bsQPfu3R0pqvvIyJC9h/XEedX/AfmEL6bJnk/0iTwH6y2RAaZZKhJ74jyaUmqD1jrl7ONGDKGciHqGTyql2gDI1lprpdRASNKUXAM+07NZeuqIyLlyd5ldAiKyFestkQF2sSJ5M4cCOKVUGIDLANxR69idAKC1fg3A9QCmK6UqAZQCuFE72uVHRERERETkoxwK4LTWJQCizzr2Wq0/zwUw15HPICIiIiIiImH3OnBERERERETkWkYtI0C2CA01uwREPqFNstklICJbsd4SGSCZFcmbMYAzgwcsk0DkDUbNMbsERGQr1lsiA8xhRfJmHEJpgxUrVqBr167o1KkTZs+ebXZxiIiIiIjIxzCAa6KqqircfffdWL58ObZv346FCxdi+/bt9l1s3z7ZiMipPpkkGxF5DtZbIgNMmiQbeSUOoWyi9evXo1OnTkhKSgIA3Hjjjfjss8/Qo0cP2y9WUWFw6YioLoWZZpeAiGzFektkgExWJG/msQHc/NRzj/UcDwy4C6goARaMOff15MmyleQAH1x/5muTVzf8eYcPH0Z8fPzpn+Pi4vDLL7/YWGoiIiIiIiL7cQhlE9W1/rhSyoSSEBERERGRr/LYHriGeswCwxp+PSym8R63s8XFxeHQoUOnf87MzES7du1suwgREREREZEDPDaAc7UBAwZg9+7d2L9/P9q3b49Fixbh/ffft+9i4eHGFo6I6hQ3yOwSEJGtWG+JDDCIFcmbMYBrooCAAMydOxeXX345qqqqMGXKFPTs2dO+i8XFGVs4IqrTiKfNLgER2Yr1lsgAT7MieTMGcDYYM2YMxoypIzsKERERERGRCzCJiRn27JGNiJzqg+tkIyLPwXpLZIDrrpONvBJ74MxQVWV2CYh8Qkmu2SUgIlux3hIZIJcVyZuxB46IiIiIiMhDMIAjIiIi27kySwAAIABJREFUIiLyEAzgiIiIiIiIPATnwJkhIsLsEhD5hMThZpeAiGzFektkgOGsSN6MAZwNpkyZgqVLl6JVq1bYunWr/Rdq1864QhFRvYY+anYJiMhWrLdEBniUFcmbcQilDSZPnowVK1aYXQwiIiIiIvJRDOBscMkll6Bly5aOX2jXLtmIyKkWjJaNiDwH6y2RAUaPlo28kucOoUxNPffY+PHAXXcBJSXAmDHnvj55smw5OcD115/52urVxpexPlq77rOIfFhFqdklICJbsd4SGaCUFcmbsQeOiIiIiIjIQ3huD1xDPWZhYQ2/HhPj2h43IiIiIiIiA7AHjoiIiIiIyEN4bg+cCSZOnIjVq1cjJycHcXFxePzxxzF16lTbLxQVZXzhiOgcXcaaXQIishXrLZEBxrIieTMGcDZYuHChMRdq08aY6xBRgy663+wSEJGtWG+JDHA/K5I3c2gIpVIqQym1RSm1USmVVsfrSin1slJqj1Jqs1LqfEc+j4iIiIiIyJcZ0QN3qdY6p57XRgPobNkuAPCqZe/b0tNl37WrueUg8nLzU2U/ebWZpSAiW7DeEhnAutwWk/Z5JWcPoRwH4F2ttQbws1KquVKqrdb6SENvOlV5Cjtyd5xxrKK6AiUVJVBKObG4rhGkq+Gn/ODI36S8qhxVusqwMjVGa42K6opz/l+I3FlJRQcE+QfDka+60spSZBRkGFYmImpYSUUHBPoHAQi0+xrlVeXYm7/XuEIReZgOFSUIDgh2qKFfVF6EzKJMw8pExnE0gNMAvlZKaQCva63nnfV6ewCHav2caTnWYAC3t2Avxi8df8axvyT+BZWHKhEUEeTxQVxCZRmC/e2vVKcqT7n0xqS1RnlROTblbsKLaS+67HOJHDU09y2EBoQC6GX3NR5f9zi+3PelcYUiogYNzX0L/n7+0Lqf3ff7Z399FovTFxtcMiLP8VbuPjQLikB3B64xY9UM/HbsN8PKRMZxNIC7WGudpZRqBeAbpdROrfWaWq/X9c2r67qQUmoagGkA0Oa8NvjXpf868/VqhaqSKlScrHCwyOaq1tXYfTwHASoA/nZOQSyvKkdOaQ4igiIQ6Gf/E0pbVPlXIaFDAv6V8K/GTyZyEytfjERp5SmHrlFYVoj4iHjcn8IJ4USusPqllsg7lefQNQrLCxETGoNHL3zUoFIReZZmL0xEVXWlQ9coLC9En9g+mNrLjozrZIjhGF7ncYcCOK11lmV/TCm1BMBAALUDuEwA8bV+jgOQVc+15gGYBwApKSl6WIdhjhTNbeWW5mLf+EvQNvI8tF6/3a5rbMjegJkrZuKNkW+gX9t+BpeQyHus8d+P0spSh66hoREZFAlv/U4icjc/BxwGIHVPOTDZIDwwnPWWfNYevyBUVTt2/wOA1mGtWY/ckN1ZKJVS4UqpCOufAYwEsPWs0z4HcIslG+WFAAoam//mC74aGIWDo+zP5SJTComoMXrYNhy/8HvHruFgI5KIbBM8ci8ODfzKsXsdb5Pk4zamdsGPF7d36Bpsb7ovR5YRaA3gB6XUJgDrAXyptV6hlLpTKXWn5ZxlAPYB2APgDQB3OVRaL6CUwuLh0dg98TLHr8VGJVGD1LVpOHL552YXg4hsEDZhG/YNd3z+Gu+R5Mt+urovlo9ONLsY5CR2D6HUWu8D0LeO46/V+rMGcLe9n+GNFBRCyqrhV2r/vBzNR4tETXMqEKosyLFraHh84iQij1IaCP+yEIfudbxPkq8LOlWJwDLH8kawHrkvZy8jQHX4z4sZ6BD5HLB+itlFIfJu99+E3mVXADebXRAiaqoTM8ZgcHFf4DazS0LkuW57aInMAZ9gdknIGRwZQkl2MHJIB4eHEDXMiBrCOXBE5uDTfyL7GXXX4v3PPTGAIyKvZsgcbN6/iFzGiBHLDP6IAEez+bAeuS8GcC5m5Fwazsshcj5m4SIyiYNVj/dIIsexHrknBnAmceS+xAYlkS0cfwLJISRErudQEhPeJ8nnOX7fYj1yX0xiYoLPBrfA6MQhaGN2QYi8nBqzCUdytwPggvdEniLsql3IOLASdSS6JqImShvVC0dOZqG32QUhp2APnIsppfDZkBbYe/Vgs4tC5PX8xm7B0dQVDl2DPXBErhV+9R4cGPKZw/NvWG/Jl20Y0wsrh8U5fB3WI/fEAM4EzYsqEXyiyO73c1IpUdPo/FAEFEQ6fB3OASBynaoTwQgqau7QNXifJF8Xll+CiIIys4tBTsIhlC6moPDi3IOIj/g38Ktja5zzqQhRw6ofvh69Tw0DpjlwEbYDiVwq774RGFTUDfpO9sAR2evmf3yOkxUngclml4ScgT1wLmbEDYVPFolch0MoiczBex2R/bgOqndjAEdEREREROQhGMC5GNeBI3IdBcdHQLIXgMi1rHc2R1OY8x5JPs+IZQBYjdwSAzgPxHU5iGzheH1hQ5DIs/A+SeQ41iP3xQDOxRQUFg+Lxs4JlxpyLSKqn9+1v+PQiE8dugZvYESu1WxCOvYOW2x2MYg82i/jzsfSy+PNLgY5CQM4E3x1QRT2jRpg9/s5pIuoaQJG7sTRQSsdugYncRO5VrNRGci84CuH7nW8T5Kv2zqiJ9Zc3Naha/D+574YwJmgdW45wo/mOXwdDusiapjOjkBIbiuzi0FENqg8GobQ3NYOX4f3SPJlUdmFiMkpNbsY5CRcB87FlFJ4el4m4iLeAEbcb3ZxiLxaxWNXonfpIGCG/dfQWsPPj8+6iFwlZ9YlGFiUCD2DvWhE9rr+yc9xeXkRcIdj1+GDEPfEVokn4j2NyKU4hITIs3DuKpHjWI/cFwM4FzOyIchGJREReStH57HxHkm+jL/93o0BnAfi5GyipjNiHTg2BIlc5/Q6cExiQuQgx+oB73/uiwGci3EsMZHrsLYR+S42PInIWzGAczEFhXdGxWDzzZeZXRQirxdwUxr2j17g0DW01owEiVwoavIO7Br1DuffEDngpxsH4cOxHRy+Dh+EuCdmoTTB9/0i0bdfX7vfz6EhRE0TMGQvsnetdfg6vIERuU6zSw/jSLPvHboG75Pk69IHd8G6rONmF4OchD1wLqagkHCkDJEZRxy/FodjEjVIH2iJ8CzHnkCyIUjkWuX7I9HsSILD1+E9knxZzMFcxGWdNLsY5CTsgTPB3+cfRvuIBcBVfze7KERerezpkehdkgw84Nh12ANH5DrHH78A/QvbAn82uyREnuvK577EJWUFwL2OXYcPQtwTe+Bc7XR6LQeya1ney0YlkfOxB47IHA7NgdO8R5JvM+K3n/NQ3RcDOBfjDYX+n73rDo+q+NrvTaOFXm2IWNJBitTQlaIgokhViiBdpCpVpBcFpUnvCvizUKSFLp1AqMmm91422c1utu893x83d7N39252NwmCfvs+Tx7I3Lnnnpk5c+aUmYkL/xwqZLaRKwLpggtPA67giQsuuOCCOMrswDEM8xLDMJcYholkGCaCYRirJC3DMF0ZhpEzDPOg+Me1Z7AC4FrUXHDBcVTEbHEFXlxw4Z9DhWQOXOukCy6gIv4OnAvPJspzBs4AYCYR3WMYpjqAMIZhzhGRxKLeVSLqW47v/KfAG4Iuo9IFF/4dcC1gLrjwdOCaey644IIL4ihzBo6IMonoXvH/FQAiAbxQUYz9l7H9/QYI+7zP02bDBRf+8/Aacxsx7+8sPyFXrMQFF/4x1JkQgcj3tz9tNlxw4V+NK6M648CAl8tNx5UseDZRIWfgGIZpAqAFgNsij9szDPOQYZjTDMMElEJjHMMwdxmGuZub+9/9uxUMw+BWgDfS2/pWCC0XXHDBNjzbpiI3UEwtOQ7XIW4XXPhnUa1DNnICbpVr7hHItUa68P8aSW+9hrtBtcpFw5UFf3ZRbgeOYRhvAH8AmEZEhRaP7wF4mYiaA9gI4KgtOkS0nYhaE1Hr+vXrl5etZxYMGPgkq1EnOuVps+KCC/95sNENUCP5jXLRIJArAumCC/8gtJG1UTPZx2U8uuBCOdAoJguvJSnLTccVCHk2US4HjmEYT3DO2y9E9KflcyIqJCJl8f9PAfBkGKZeeb75X8DXBzMR/N0fZX7flRFwwQXHoF7bBUG/zCo3HZcD54IL/xxyVrbEmwe/LhcN1zrpwv939Fx/GlP2xZeLhmsePbsozy2UDIBdACKJaJ2NOo2K64FhmDbF35OW9Zv/BZREMso/KVxGpQsu2EEFTBFXFsAFF/6dcK2RLvx/hkv6/9sozy2UHQF8CuAxwzAPisvmAWgMAES0FcBAABMZhjEAUAMYQi533gUXXPiXwbWFxAUX/jm4ZpsLLjw7cAVCnk2U2YEjomuwo2eJaBOATWX9hgviMGUEXHPKBRdKBYPy57pdMScXXHg6KPclJq5F0oX/93D9Hbj/KirkFkoXnAXjmhIuuPAvgssQdMGFfx4u49EFF54+XDtQnk24HLingA0DG+LmlLL/bXNXRsAFFxxDlck3EDFwQ7npuBw4F1z459BgRjjCB64vFw2X8+fC/3dcmtAT24aU8+/AuabRM4vynIFzoYx49IY33gp8pdx0XEalCy6UDs83s5Hv/rBcNFyG4JOHJjISinPnARAMublwq1EDDWbNckV+/5+iaksppNqH5Q5WuuTnvw3143DIjx9HtbZtwGq1gNEI0htglBWg0muvgfR6AIBRJoNHw0YAw4AMeoAlGAvy4dWkCVeHZUFEYBUKuNeqDY+6deD12mv/evnJaP4ywmtEPG02XHhCcDlwTwHNY4vQyJgItHzanLjgwn8b+oeNUDuuWfkJ/bvX8Wca+QcPInvFSsBgABgGKDbaq7Zqherduz9l7lx4GlDfq4e6sc2fNhsuPMPQ5+Qg6eOPAQAFBw5UOP16U79A/UmTKpzuP4kXHiUjIENebjquZMGzCZcD9xQw9fcsNKp2Ahj9Q5ne5zMCrknlggulQ72pAwKVTYBlZadB5LoM4UlB/egRspcug3eXLnh+1Uq416qFvC1bkLt+A/J+2gLv4GAwXl5Pm00X/mFk/xCIQPmXoLnlyMCRa438L6Pw1CnT/2sNGYxaH34I95o1kbt+PQpPnUbjPbvhXrMmwDBIHPAhAODFn36CR726MBYUIHX8BFRt0wb1Jk0E4+kJAEge/gkAoEqrVsjbsBGGrGy416wBMhjh0bABqrVpg8r+/v98Y8uIrlvOIUgjBRaVnYatHShEBNJq4Va5ctmJAzBIpdCnpqJys2Zg3Bw71UVGI5RXrsCjfgNUCQwo1/dF6RNBEx4Oj4YN4dmgQbnpGZVFAGuEe40aFcBdCVwO3FODa1uWCy78G+DaQvnkkP/zz3CrVg0vrP0ebtWqAQDqTZwIryZNkD59BtJnzcZzy5fBvXr1p8ypCy48HRBRubfyGeVyaOPjUaV5czDu7hXE2dMDEaHw1GlU8vFB02NHBc+e//57NFqyFO7e1Uxlr9+8ATcvL5OOAYDXb1yHe+3agr5949ZNMF5eYDw8kDF/AWR//AG4uYFxcwNptQCAF35Yhxp9+jzhFj7b0KWlIXXCBOgSElF37FjUnz6tTDKqkUiQPGIkWKUS3m/3wIs//gjGw75bkrVkKWS//gowDJ5bvhy1PhxQlmbYRO66HyDdsQNMlSpovGM7qrZuXWZaqrAwpI4bDzIY8PyqlRUqO65LTP6FMGXg/uX7s11w4d8CVyS/4mHIz4fi9BnU7N9fYFgBQI0+fdBw7hwozp5FTJu2iGnXHnG9eiFt6pfQZ+eI0tJnZ/9TrJcKYlkU3Q6FJjr6abPyRMCq1SCD4Yl+w5nZps/IgEYisTovR1wKzikYlcoKuSRMcfESUidMRP7Bg+Wil7//AKJbtETiwI+hz8wsEw31w4eI790HycOGI2nIUBjlzm+pY7Va5Kxfj4z586HPyioTHwCgiY5GymefIXX8BGjj48tEg4xGZK9cCc2jR6g1eJDVc8bNTeC8AYBH7dpWOsajTh0rG8q9Vi24Va0KxssLL3y3Br6PHnI/Dx/g9atX4NW0KaR795aJb6t2EEFx/jwKz4SYzupVNBgGABHytm5D6pQpUF677jQNywAmsSwyvp4DQ3YOvLt3g3T7dhSeOOk0XVarRfr0GXCrXh11Px8L5fkLkO7ZY/c95bXrkP36K2oNHYKqbdsi69tvyyxLYtAmJEK6eze8e/SAZ8OGSJ81G8bCwjLRYrVaZHz1Ndxr10alV19F5sJvYCgoqDBeXQ7c04IrqO+CC08cFeF2uTJwQqju3kXe1q3Qp6eXi07BwUMgvR61hwwWfV5n5Eg0+f131Js8GTXe7YPK/v5QXr2K1M8/h1GhAMA5Eznr1yOuS1fEdemK7NVrQEZjmfhRh0dAcelSuZwTIkLGrNlIGTkSif0/QMGhQ2WmBQC6tHTkbtgAdXj5LiLQZ2cj7YupSOjXD7Lffy8TDdLrkbN2LaLbtEVM23aQ/XmkzPzI/vgTScOGI2vpMrAaje1v2nF+5CdOIr53HyR++BEy58wFsWyZ+NFnZCBpyFDEtH4LyZ9+Wi4jqzDkLNImT4YqNBTZS5Yif7d9o1SUzrlzyF6xApUDA6BLSkLq+AlgVSqnaHCZkolw8/ZGg1kzoYmKQvrs2U71ExEhc8FCSLdshfzYcSR/8in0OdZBFHvQxsYiZeQoaCSRUD94gOQRI6FLSXGaTubCb1Cw/wBqj/gUtYcOdfp9Z8C4u5ucPI/69VHro4+gefgI2sTEctPO370HaVO+QPq0aUidPPmJBUXqylnk/vgjVHfuInXCBKgflu9SL9lvv0MdFoaGc+fixQ0bUDkwENlrVsOoVDpFR7prF3TJyXhu6VI0mDkT3m/3QN6mzdClptp8hwwG5KxeBc/GjdFw7ly88P13YCpXRvby5RV2O7t0xw4wXl54bsliPP/ddzDk5nJntMuA/L37oE9Px3PLluL5NavBqlTI21hxfxr7X+fA6bNzRCOwZYEuLQ15W7ZA/ehRhdADAH16OtJnfwXp3r3/6ev+WY0GisuXYZTJnjYrzxSMMtl/etz/rXjat9mxWi2kO3dC/tdfZeZFm5CI7JWroLp/v0zvGwoKkD5rNrJWrADpdGWiUfC//yH5k0+R++N6JA4a7JQhZ26gGAsLkb9vH7zf7oFKr79u850qgQGoP2UyGn3zDV784Qe8tHkTtAkJSB03HtK9e5HQtx+kW7aieq9eqDlgAPL37EHmwm+cNuTz9+1D0sCBSJs4CckjR3FnFsqAggM/o/DUKdQZ8xmqde6ErOUrymwwGQsLkTpmDPJ+2oLkTz6BKiysTHT0GRlIHjYcyuvXwXh6IXPBQhQcPuwUDSJC1tJlkO7YiRq9e6Oynx8y589H4blzTtPJXrUamfPnw5ifj4KDB5E+fUaZjFfp7j3ImDULlZsFoc6oUZAfO4b8ffudpqPPzkbyiJHQxsWhzsiR0Dx6jLQvvijTHFE/DkfGV1+hSvPmeP36NVTv1Qs569Y5LQO6tHRkzp2HyoGBaLx7N1748Udo4+KQuWChw/qDd7xIr0fjHdtRd+xYNJo3F0VXriJvk+NGZP7efSj86y/U/3Iqmhw6CEN+PlLHjTcFURwBW1SE1ImTwHh6osmvh/HyoUOAwYDUceOdcpYLz5yB/M8/UXfCeDSaN+8f34VUo29fwM0N8uPHy0VHFRaGnHXr4N2jBxrOnYOiK1eRu758fzZDDJ4qHerJWNR47z28du4sPBs0QPrsr8AWOaff+H425OYi5/vvUbVtW9Qc8AEYNzc0WvQNjHlSpxwTXVoapNu2o3rv3vAO7ggAaLRgARh3d2QtWWpTxmW//QZtbBwazJoJNy8veNSrh/pTp6Loxk3BmciywpCXh8ITJ1BrwAB41K2LKkGBqDd+HORHj6LwTIhztAoKIN2xA949eqBau3ao9NprqD10KAoOHy63E20CET1zP61atSIxFN27R5FBzSgyMIgKL1wQreModOnpFN22HUl8fEkSGERFoaHloscj4eNBHE0fX8r/3/9E6wxe6kMHfpnpFF1DYSFJ9x8g1f37dD7pPAXuDaQoaVRFsOw0WJalpJGjSOLjS7Fvv0MGmax89PR6ylq1mjIWLCRDYWGZ+Mndtp0yFy8uFy/GoiJSPXhArF5fpvez1qwhiY8vpUyeTKzBUGY+WJalnE2bKGn4J6S6f7/MdP4/gWVZkp8JoYLffyejVmsqX3/4V+q4+CMyssYy0czZuImutQuiQ5N7l2tMU7+YWqIXDh1y+n2DTEYxwZ1I4uNLkc2ak1oicZpGyvgJJh4yFi1y+n3V/fskCQik5DFjSfXwIUUGNaPUKV/YfY81GCj1y2kk8fGl5M/GkFGlouzvvyeJjy+pIyOd5kN+8iRFtWpNEh9finvvPVLeum16lrN+A0l8fClzyVJiWdYxeqfPkMTHl1KnfEH5hw6b2sganZMZdUQERQYGUcr4CcSyLBlkMort1p1iunUjfX6+U7RYo5FSxo0nSUAgyc+EUFyv3hTVpi1pk5OdppP0yacU1ao1qR6HE6vTcXLg50+F5845TEe6dy9JfHwp+/u1RERkVKkoYdAgimzWnFQPHjhOZ/+BkvExGEj6yy/c70uXCer9cvQCdVz8EaUr0kXpyE+e5Mbsy2lk1GqJZVlKmTCRIoOakSY2loiIPg/5nIafHF4qP6xeT4nDh1Nki5akevSIiIhkf53g5siChQ7LEBG3fsT17MWNt1RKREQGuZxiu3Wn2Hd6kkGhdIgOy7KUPPozimrRkrSpaaby3G3bSeLjS3m79zhEp+D334v1zWEB7fS580ji4+vQ+Kujo0kSEEipU6aY+kJx9RpJAgIp6dMRAl1bGjIWLSKJrx8V3b1rKiu6e5ciA4MocdhwMmo0dmnoc3Mpum07Shj4cZnX6IpA8mdjKKZzFzKq1WV6X5eZSdHBwRTbs6fJ3sn4ZlGpNmNZoLx5k6J836DbzZuaZK8oNJQkvn6UPGasw7qky+Eu9O2Nb4nV60vmWEKCoE7Gwm84nXLxol16rF5PSaNGUWSLlqTLzBQ84/VM1qrVVv2ry86mqNZvUdKoUYJ5yer1lDDgQ5L4+lHazFmkvH3bqXlrjpwNG0ni40ua+JL2sTodJXw8iNOhDx86TCtz8WKS+PqRJibGVMbrg+j2Hago7J5DdLSpaQTgLon4Sk/dWRP7EXPgWIOB4t59j2I6daa4Pu9SVJu2ZJDLHeoAK1osSynjxlPkmy1Idf8+p3S7diOj0jEFawuqx+Ek8fEl6b793KLZ+i3S5+RY1WuxvwWtu7vOYbpGrZYTUB9fkvj60dU/NlaIA2csKiqToCuuXCGJjy+lzZ5NEv8AyliwoFx88Iu4xMfXanI6At4Ak/j4UsKgQU4bXkScQRL/fv8SHnQ6p95XR0WbDEqJjy/l7djhNA88ZMeOmcY6uk1b0mVllZlW0Z07lLdrt8mgcBasXk+yY8dIeeNGmZUiazCQ4spVp41Ph+nrdJQ+Z65ABnjd8NP9nyhwb6DTDhzLspS5eAlJfHzpfNfm3KKyek2Z+Cu8cJEkPr6Uu2UrJY8ZS5HN3yRdurhhaguZS5aSxM+fFFeuUnRwMCUOHebUeKgePeJ42L6dslatJomPLylv37b/YjFYvZ7i3+9PMV26moyO3K3bOEPw/PlS31X8/TcX2Bg/gSS+fhQ/YABJAoMo/es5Dn/fEka1mvS5uVZ9wLKsqX3Z339vt49U9+9TZLPmlDhkqMmIzD90iOurrdsc50eppLjefSgmuJPAWVM9ekyRgUGcQ+hEACBn/XpuLfn5ZyIi0iYnU3SbthTft69T6x6vWwt++62E16Iik/NVFBZml4bi8mWS+PlzRryZbtXn5VHs2+9QdPsOpE1JsUtHFR5OkoBASpkwUUAna+Uqrq37D5jK/oz5kwL3Boo6cHwgN3HoMIHhr8/J4Qz8QYOINRgccuB4h1927JigPHvtOkH/OwLeSbGcV0V375LEz5/SZs12aM4W/PEn9+1ffhGUsyzLBYL8A0h581apNPT5+RTdpi0lDh9utR4aNRpK+GggRbVsRZr4eJs0WIOBEgcNpuh27a0CELLjx00OtD25Vly7xunPlausnvGOeNLIUWRQKGzzwrKUMmky5zzExZX6vScN5a3bJn0mO3qU8g8dprwdOyhv9x7SJiaW+q5RpaKEDz/i+t7MsGe1WkoePZokPr6UOHQYSfcfIE1sLKklElJcuUoFR45QwZ9HSB0VRawdp5llWcrbs4ckfv50q3s76ruxteB5/q+/UmRQM5L4+FL8gAGUuWQp5e3aTfKQEFLeuk1FYWGkuHKVCs+fp8KLF2nsqva0efdESh77OSeXB6znhFGp5GxUP39K/+prUly9RrrMTDLIZGQoKCB9Tg4ZZDJiDQbKWLCQ00m//27Nu9FocmZjOnehnI2bqPDiRSoKC6PEIUMpsllz0T42KJSUtWo1RbV+iyQ+vhT1VhtKHj2a5KfPkCY+ntRRUaSOiCi17/S5uRTVshWlTJ5s9UyXkUGxb79DksAgyly8mOSnz5AuO9smrbxdu7nA1PLlVs80CQkcLf8Aylm/vlRbU3n7NkkCg/79DpzsBBcVk586ReqICKcXWHPIQ0IEkayisDCbCsYZZHyziCKbNSeDTEaa+ASKDAyitJmzrOpN+Oo1+m3zeIfp5m7ZWizwf1B83770KLg9tdoeUC4HLnfbdpL4+nFRNJXKqXeTRo6imE6didVqKWs1l3VyNJpgCdZgoNiePSlh0CCS7ttPEh9fUly56vj7RiPF93uf4vq8a1r45CdPOs1H3s6dJPHxNTkClounPaR+OY2iWrYiQ0EBpUyeTJFvtnDaSCfinJHY7j0o4aOBpImLo8jmb1La9OlO0yEqnvx+/pxj2bMXGQoKnKaRPn++yTEqiwPDGo2UMmkylzkKDKL8gwedpmEPWStWkMTHl3LWb+Bd95dJAAAgAElEQVSi54FBlPTpCGJ1Ovpx6zHqPnssGYzOZc/yDx40RQI/ONKffhvXs1g2rzhFx6hSUWz3HhT33nvE6nSkS0ujyMAgp4Ie+rw8igxqRhkLv+F4O3SY4+XyZYdppE2fTlGtWpNBoSCjWk0xXbtR/IABDgc78vbs4eZWSIipjNXpKL5vP86pKyW7kDZjJkW1aUtGrdZEJzo4WDS4VRFgWZYzpovns1oiEY3Ya+LiKLp9B4rt2VNgoLIsS2nTZ5DE14/yDx12yOhO/+przngXMazzD//KyefGTQ7xz69P6fPmCb6tvHGDJIFBFN+3H6mj7Ot+XXo6RbVsRcmjP7Nqgz4/35TVMzeIWZYlbWoa6dLTyVhURLLjxymqRUuKHzCAjEVFVt/QxCdQVJu2FNezF6mjom3yYlSpuABs5y5Weog1GChl4iQugn/pEhER7d97jbrPHktpijRBXV16OkV36Eix7/QUzWrKjv9lClSMDRlLn5z8xCZPylucfhQLJLBGI6VMmEgS/wBSXL1mkwaPwosXS9WROZs3czKwYWOp8qTLyqaot9qIOl5EnKEa9957FN2ufanrS8aiRSTxDyB1tPiY6DIyKLp9B4rr2Yt0WeKGKJ8ttXRuefBGasaiRTb1iKGwkGK6dqO4Pu/azFjJjh4lSUAgxff/wCYvsqNHOZtt127R5/808nbvMTlBlj/R7TtQTJeuFNezFyV8NJCSR4+mrNVrSPbXCUocMpQkvn6imSpWp6O8Xbspvm8/UbqmH18/iu4YTAmDBlHKpMmUtWo1yY4eJfmpUyQ98DMljRhZvKNgCv1v7SiaPMfH6lu6rGzK3bqNSzS0aFn694p/It9sUeruEYNCSVkrVlBkKfQiA4O44Nra0hMYRaGhlPTpCJL4+pW8HxhE8tOnS33PqFKR7OhRyli0iGK6dbP6fuw7PUl+6pTVHDQoFJQ4bDhFNmtuM6ihl0op/auvBeMe270Hpc2cRfKQEFJLJFR44SKlTZ/O9f8XU20GNwxyObdm+HC72KT79pHqcTgpr18n2bFjlLd7D2UuXkKRLVpSXO8+/24HjjUYKO699zgjqFhRJI/9nKLbtis1aiMGo1JJMV26Uvz7/QWLesY3i0ji5+/UdhBLulEtW1H6V1+ZyvjonvLGDUHdu77elNLyVYfo6qVSimrRklImcVGBorB7JPHxpdWf+pXZgeMVc1yfd7kowTLrKIEtqMLDBRkmW/3pKArPnTM55qxWyzkvAz502LDkjR3Z8eOcnPTqTQkffuTc1heViqI7dKTk0aOJZVlKHDqMy8g6uD1EHc1l37J//JGIiDPSm79JqV9MdZgHHvy2F96Qydm4yelsCRFnoMV06kxxPXtR4cWLXDBh+gynaBReusQZJStWmBw51aPHTtGQnz5dnA1Zy20J8/GlgiNHnKJRGvjgi/kWLH6xz1y6jJa3SqfJPqFOOXDa1FSKDGpGyZ9/TqzRSP2P9KdZZ6dS3LvvUWy37qKGrC2Yxs9sm1/mkqUkCQgkbWqqYzSKjT9+Wwer01Fst+5cpsEBOdemppHEz5+y1pQYlwV/HuHk7IL9LS8GmYyiWrSk5HHjrL5XdO8eSXz9KHPJUvF3FUqKbP6macsmy7Kkun+f9Lm5dr9bHrBGI2Wv+4Ek/gHcYuvnTzHBnSjho4GUsWgRZS5ZSlFt2lJ0x2DRiK6xqMgUEY/t2ZPS58ylnA0bqfDSJdKmpgkWZt6QzdmwUZwXljU5eBkLFlLBH3+S7MQJUt68ZbV+FYWGUuSbLShx0GDRbWXK69cpun0HkgQGUd6u3TYNBKNWS4lDhnLb8Gxkx7SpqRQdHEzRHTpS1qrVlPrlNIoODrYyehKHDLVpWBNxa1J0x2Buq920aaS4dk3QLqNazTlovn6kuCbuDBmVSoofMIAiA4Moc8lSWhuURpN9Qim1sGSOqB6HU2y37hTVqrXNLAzLspQ69UuS+PrRiu8H2HTg9Dk5nH7s1dvmzhuDQkHx/d6nyKBmgsCFJbQpKRTdvgPFv9/f5lZA1mAwbV3MXLxYdK1kDQZua1nzN0vN5GjiEyiqVWuK7dmTCi9csJIBVXg4NydFov/mKAoLo6gWLSn27XesAg/qiAiKbNacksd+XqqO4bdCp82YadV21miktBkzHbKrFFeuUmSLlhQdHEyy43+RXio1fZdvb+LQYeXaxl7RMCiUpE1MJF1mJhmVStJlZlLezp2UsWgRpc+dR2nTZ1Dy559TwsCPSVLsuES1bGXXCSHibArZiRMkDwmhorB7pE1KIk18PMmO/0U56zdQxoIFlDRyFMX37Wdyisydirw9e4g1Gim5ZVO661u91G+xLEsGuZxU4eGkvHmLFFeuUlFYGKnCw0n18CGNW96Otm4b73Dm36BQkPLGDco/dIik+/aRdN8+yj940LT7Qx4S4rB9ppdKSfXwIRVevCjYTuwIWIOBlDdukOwvrh9lR49SfL/3uexe126UMPBjin+/PxcQaduOJP4BJD91yj5drZZUDx5Q3p49lDr1S4ru0FHQ/1EtWnKZNQdktfDSJc6pF3N2mzWn1ClTSJeV/e924Piomnlmhd8SlLN+vd1OMoetjJFBoeAiRT17lWlrZsHvf3B079wxlRk1Gop9pyfFdu8hMFicceAyFy8hiZ+/YMG6M30MPfb1pcib9oXNEsrr17n0/+TJxOp0lLlsuZWBWRrSps/gMk1mZ9XkZ0JKFiYnti+yLEsJHw+i2O49TAtawZEjJofOkffj+39Acb1Kzifx0W5bhoIY+Mwffw6S3/Ih3bffofdTp5Vk33jwWVOnsol6PcW+/Y7AATVlcPq867BDaToPEhhE6ogIIirZ220un6XBqFZTbI+3Ka7Pu8RqtWRQKCi6fQdKGjXK4fYQESUOGkyxPXsSazAQq9VS0qhRgjOnrE5HBoWCNHFxpMvKdsrxZo1GShjwIXcewcKpylzOyfVqn3ia7BNKeqPjwYWMBQsoMjDItD+//5H+NP3SdG4rlI8vZa/7wSE6RqWSotq0pdQpUwTluqwsigxqRunz59unodVSdMdgSh43TlBuknMHsgNZK1aQJCBQcN7AlOl1wAnks2a2zt3xOkR+xtrI5R3FsmboywtdRkax0bOe0ufPp+TRoymqRUuKatWakkePLnVbL2s0UsGfRyh57OcU06mzIBoc9VYbSpk0mYsS+/hS6tQvS12wjWo1t0PDImof3bYdSfftI7VEQnm795AkMIji+rxbanZSn59PKZMnm95PGjWKMhYspOzvvqO8PXtI9tcJUxTe3m4EdXS0yVGN6daN0mbNpvyDByn/118pd+s2kp044ZAhopdKuS1Mbdqa2hbTrRslj/6M4nr15vSpnV0Npgh3YBBtrnKb1tW4TTHTp1Dm4sWcU+bnTzGdOpPqcXipdIxFRZTw0UB67OdHuyZ04Yzgu3dNOkKfm0sJAz+myDdb2D1Lqs/Pp4SPBnKO7NBhVoE0TXw8xXTtRtFt2pa6HZGIkyfe/ojr8y7lbtlCRaGhZCgoIJZluYCDj/jWMksU3blDsT25XQGx3bpT7k8/kepxOOkyMynu3fcoumOwQ3aM6v59iu3WneOpV29KnzuP24USGMRtB7YTaGFZlnK3c2fzYjp3ofSvvibpvn0kO3HCtOXO0Z1S6qhoiu/bt2SOFW8ZlgQEUnTbdqRLc86Af5ZgUChIHRXlVPDPUbB6PWni4kgTE2O1hjriwNlD58OdacmNJeVl85kAazBQwZ9HOOd6zFhKmTiJUqdMobRZs8u+i0yvJ+XNWyQ/E0JF9+6VaYzVUVGcngoNJU18PBnkcsE4/msdOGNxVia+/wdWzkHajJkkCQh0+KIHdVQ0SfwDbBpORXfucIdzR45yyttnWZYSBg3ijF0LY0j16BFFvtmCYjp1JvnJk8SybLED19QuXd7Zsjzgfe7xUbrRwpckA/o5FZEyKpUU2607F3Us3s5gVKko9u13KK5nL7sHibXJySTxD6CsVasF5ebnTtKmz7C7R5sH7ygJDlkbDBTfty/nlNnJ6PHZO/OMjlGj4aIrHw10yBkwajQU07kLJQ0XRmqTR4+m6Hbt7V6qoomJIYmvH2X/IDTqjVotxfXqTbE9ezp0QJvIPCsivKCHP0OUs8nBLVhnz1ptNzGqVNy2uX7vO3S+j3f4zCOz/AFjy4yyLRTd47LF5udIDDIZxfXuQxL/AIo2M/b4n+TRnzm81ZPPVsqO/2X1zKjRUNx779HG6nfpi9dDSWd07EyjNjWVJAGBgozS+0fepxmXuOxl2uzZFBkYZPesA1FJpltMP2UuX04SP3+7l3jwAQ3LgIRRq+Xk/ONBpeoAg0xGkS1aUtrs2VbP+G2ipZ2nYY1Gin2nJyUOs32WyKjVUsKgQRTVoqXVlq2kUaMo9p2eZT4/+STAGo1l4seoVpPy5i3K//VXSv/qa4rt2ZPiPxhAOZs3O7z7wKjVkjYlhTRxcaS4csXkaJnk//PPHbqIib+0J33+fEr4aCBFBwcLnMOoFi2p4E/HM90VldkwqtVUePEi5W7dRmkzZ1HChx9R4pChpPj7b4dpGGQy+uG1eFpT+zZFdutC0W3aUmy37pS5fLnjWQC5nH4f3onCfc30i58/xXTpSpFvtqDIZs3tnt00tUmjIen+AxTTpStJfHwpYcCHlLVmDWUuXkyRzZpTdPsOpkCZI5CfPUuJgwYLdV9xpthy22xpYHU6kp8+bbpQzBS5D2pGyps3HebHqFaTdP8BSpkwkaLbtaeo1m9R+rx5pZ7xsYTyxg1KmTiJy8SayaD055+dDsopb94i6d69lLFoEaVMmEhZa9Y4dMbSBWukuBy4/wT+tQ5caeeiDHI5xXbvwd2EaOeGJ9ZgoISPB4keyDVHwR9/UmSz5iTx9aP0r+c4dD6MP9Rq68Cz6tEj7uB+sYPjiAOnz82l6OBgLvNiwcPZpLM0frZ/8da07+3yxyNz6TLuJiiLg+u8I2Uvm5k+dx5FNmsuqthZlqW8HTtI4sMdSHbkhqbE4cMppktXq8xS4fnznGNXyo1MLMtS/AcDuOyOhfHEZ0PFMgKW4OXL0ilRR0SQxNfP7rmv1GnTKKpFS1GZ4h1wR86+sDodxb7DGYRiC17a9Omc45CUVCodo0pFsd26c46aRb/wDm/WylWlZkq1KSkUGdTMasulyTke8KFDBmvq1C+5s08WESltcjIljxnLbUvbvJnydu6igj+PUO6WLRQZGERxvXrbjWYb5HKK7tCREgcPsWkgqMLDaVPV27Si7m3SGhwLKqTPn0+RQc0EF8eYO3Cmm7CGf1Kq0cvq9RTb421KHDJUnH+ZjKLbd6DEIUNt0mFZlhIGfMhtHRdpoyPnQvjouFimwajRUHRwMCWPHm3zfcXlyw5lcnRZWRQT3InbbVB8YY4uK4skvn42txa6wI2xOjKS5KfPkOrhw3I5uvzNl+ro6DKdd32W8H2bXJrsE0ophWU33MeEjKFx/xvCnU25eJFyNmyk9K/nUObiJYILJByFUaUi6d69JVvi/AMobfqMMl8ypc/NJcXly5S3azdlr11HBUeOlOkCLiLuTJPsxAnK3b69TG2rKLAsS/rcXFJHRTt9vMWFiofLgftv4F/pwBkKC03btmwtbKYbnmbOKnXxk/78s81ovSV0GRncVodiJ84ekkePpujg4FKdFtZgMGU1Ep6vZteBS536JUU2f1P0EHJIYggF7g2kyDnTHG4Tv/3L1nm3tNmzSRIYZPPQsyYhgST+AXbPyxX8eYQkvn6UOm1aqYtRUWgo5/Sa3TzGg2VZShw0mGI6dbbpQPNOnliU2XRmsnefUh0No0pF0R2DuQsvRGQnfd48kgQGWV2Zy0MdHS2afTNH2vQZFBnUzK7jxd98x599s4QuO5vb9vXZmFLlnL+9TuzPYrAsa7rhKWXceJtGXsqkyaJX/BKV3BqWu217qe3RJiWRxM/fdO24oyi6e5ei27WnyDdbkHT/AZsylLlsOUl8/UgVXvp2qrUvR9LmKrcpZfp0u8EYbVISJ+MWGe9+R/rRzMslf/aDvywnff58m/zJT53iHJ+zZ21+j8+u2Qqc8OddbR0cZ1mW20rn5y86D4xaLcV06lyqg5a3cxcndzauf07+/HOKCe7kUNZW9fAhd6Pj8OHEarWUu2ULSXx8n9jtoy78d2Fy4OTlcODOjKERp0ZUIFclYFn2mcoqu+CCGDgHzrtcNDod6kRLb4qfcXbhn4EtB+6Z/UPeRITs1athlMnQcPZsm3+wsWqrVqg/9QsUnjiB3A0bRP/4pi4pCbnrN6Bqu3ao0fc9u9/2fO45NPxqNupNnAj50aOQHTlqs6764UMU3biJuqNGw61yZZv1GHd31JsyGdXfeRu6Kg1wbUT3UmkqQkJQ9/OxqPzGGzbruU3/HFVat0LGnDnIWbsW2thYziu3AKvRIHP+Ani+8AIaTPtSlFbDuXPh7u2N9C+mQpeWLnhGLIvMhQvhVrUq6o373CY/AFBrwAdoMHMGFKfPIHfjRpv18rZsgXu9eqj18UCrZwzDoMHsWTDk5CB//wGr58SyyN20GZ6NG6Nmv77W77u7o8G0adAlJoq+zyN//wEY8/JQ/8upovLVYNo0uFWujMwFC8GKyFXeli1wq1oVdUaOtPmNBl9/DcbTE5nfLALp9aJ1WLUaeZt/QpWWLeHdpYtoHc8GDVB/+jQUXb8O6bZtonW0CYmQ7tyFGn37oupbb1k9ZxgGjb5dhIYLF0B54waShgyFLjVVUEf5999QXriA+pMmwrNRIysa1fv0QfVevZC3cSM0EonNduf99BMYDw/UGfGpzTpiqNqqFV45egRVW7dG9vLlSP50BBQXLkCXkgJtQiLUEREoOHwYBT//jNpDh6BKQECp9Kqsf4isfouhPH0GScOGQ5uQaLNu7oaNYDw8UNdCxi3nVK0PB6DepImQ//4Hspcts3pORMjbsQNer7yC6t1tz/Oa/fuj5kcfIu+nLVCcP2/1XLp7F9yqV0fN998XfZ9hGLzw3Xeo+tZbyJw7F+kzZwn+sLbsf7/BkJODumPH2uSh9qefoJKvLzLnzoM2IUHwTJecjKIrV1Fr8GAwnp42afCo0qwZnlu+HOq7YUgcMgTSnbtQrWNHeDVubPddF1wwR+NljxA2agkIZfuj908aDMP8439E2gUXnMX5uR9j5WdNnjYbLjwheDxtBsRgyJMic+48yI8eRd3x41HZ37/U+nXHjYMuMQnSLVsh/+NP1Bk5EjU/HAD3WrWgjY1F2sRJYNzd8dzib51SuvUmT4Lq7l1kLVmCKs2boVLTplZ18rZug3vNmqg9ZLBdegzDoNGSJbhz8wJ8D4ZB0yfGykEjIuSsXQf3unVRd9So0ul5euKlrduQvXQJpDt3QbpjJzwbN0aNPn1Qs19fVHrtNQBA7g8/QpeUhMZ7dsOtWjVRWh61a+PFzZuROm4cEj/4APUmTkDtIUPgVq0aCg4cgPpuGJ5bvgwe9evbbWedMWOgTUyEdMtWsIUKNJg5A25Vq5qeq8LCUHTjJhp89ZVNp7dq69bw7t4deVu3wuvlxqjeq5dp7ApPnoQ2MhLPf7cGjIe4CHv36AHvt3sgZ+1aGBWFqN7jbVT28wXj7s7xcO8+8n76CdXfeRtVW7US75P69dFowXxkfD0Hie/3R63Bg1H97R7weuklqO7eheJMCOqOGweP2rVt9oVnwwZoOH8+MufNQ/qMGWi0aBE86tXj098grRZZy5fDkJuLF35YV6p81h42DOp795H743qoQkNRyccXns8/j0qvvw7Gwx1ZS5aCqVIFDb/+yiYNhmFQZ/hwVPb1ReqkyUjo2w/Ve/VElWbNQXo9pDt2wOvVV1FnxAib7zda9A0S7oUhaegw1Hj3XVRt3RrudWqD9HqQTo+iWzchP3ac6xsH5MW6zxripe3bID92DDmrViNt8hSrOpUDAtBg9my7tDxeKcSRd1MwfcJG5Hw1Fwl9+8K7a1dUDvCHW6VKIIMBRoUC+owMKE6fQb1JE+HZoIF1uyEcl3pffAFWo0X+7t1QR0TAu0sXVGneHJ4NG0Jx/gK0kkg8t3KlSd7EwDAMGn3zDbTRMUj7chpqDxsG7+COYKpUgSY8AsrzF1Bv6heCuWMJtypV0HjnDuRt3Qbp9u1QXr6MWoMHw61yJeTv3Yeqbduiavv2tt/38sKLP/6ApE8+RdLgIag9fBiqtmyJSr6+yF2/HoynJ2oN+thuP/OoWRwgy1m7Fp7PNUKjxd86/K4LLvCo3LQIyoykp82GCy78q1HQpAGSNJWeNhsuPCEwYhmbp43AylXot6ZNUXvYMDScNxeMm/1EIRGh6PoNSHfthOrmLa7QzQ1gWbjVqIHGu3aiSlCQ07zos3OQ+MEH8GjQAE1+PSxwODTR0Ujs/wHqTf0C9SdNcpjm8rGvYchVAuPfHK/8+YfAaFdcuoS0iZPQcOEC1Bk+XPT9kKQQzPp7Fo68fwSv1eacNH1mJpRXrkJx9iyKbt4EWBZVW7dG5aAg5O/Zg9rDhqHRNwvt8qZLTUX2suVQ/v033KpVQ9X27aC8/De8O3fGi5s3OewAk06H7O+/R8GBn1HptVfReP9+k6OTPGo0tHFxeO3cWbhVqWKThiE3FynjxkMbGYnK/v6oM+YzMB6eyPr2W3g+9xya/P5bqbJhVCqROW8+FGfPAgDcqlaF16uvgowGaCOj4PFcIzQ5fFjUYDeH4tIl5G3ZCs2jRwAAzxdegD47G57PPYemR4/YdIrNkb9/P7JXrwGMRsDDAzAYBM/rjhuHBjOm26VDRiOkO3dBfvQo9BkZIK3W9MytalW8sHEDvDt2tEsHAHQpKZDu2AHFufMwymQAAK9XX8WLGzeIBivMoc/OQe66dVBcvgxWLhc+9PBAnREj0GDal2C8vBzixRZIp4Pq/gPoMzPAuLvDrUoVuNepg8qBgXBzgPbGjWfxV/xx/LV2LdwLFMjfuxeFp89An16SZWYqV4Z79eqo2r4dnl+2zIrnfkf6wa+OH9Z0WSPkjQiy//2G/AP7oYuLFzyrFhyMl7Zvc0h3GQsLkb16NeRHj3HyUYyqb72Fl3ZsLzWzbw5dUhKyVq5E0bXrAMuickAAXty0UTSTavVuSgqyli9H0ZWrgNmaUG/KFNSfMtmh77vgQkXh5x2h2B+xHzuWzMbLNV4uE40xIWNgYA3Y12dfBXPnwr8NmXI1Knu4o3a18q1H/zYc/WEMLqdexo/r4u1XtoFOhzuhV5NeWNBuQQVy5oIzYBgmjIhaW5U/iw5c65YtKfTmTejdPVDJw3YE2xbU4RFQhYbCWCiHR/36qP722/Bs2NBm/egsBR6kFmDwW+JbfZRXriB13HjU7P8+Ks2aixp1asDT3Q1pX0xF0Y0bePncOTyUsWjbtK7pHSNLeJwux5sv1bKiF+ZXAy/qvZHvURMvbFiPGj17AgBIr0dC/w8AIjQ9fgyMpyduJUih0hnQ3beEfzEHzhz6nBwUnjiJ/J8PwJCRiWod2qPGug345X4WJnZ5FW5u9p2wqAs34HHqKNhbN1C1VSs8t2I5mKrVEJFRiKAXa9p939R3168jbdJkeDVtipe2bYUmIgJpEyehwZyv7WYYAYAMBsiPHUfe1q3QF2/3c2/YEC/v2YNKTV9xiAdDbi6KbtyA+tFj6BITwXh5odJrr6LuhAlw9/Y21TsTngUiQp+g50Tp6JKTofz7b6jC7qGwak1gyCfwb/66QzwA3BZHxdkQsCo1WDd3uLm7wc3LC5X9/VEtuKNN5zhdpsZzNSoLxk1nYGFkWXjKC6COiABYQpXmzSD1rIaGNYQG/8WobLRvWg9VvMTnErEsDHl5YNzd4V67tkNOh+ldoxG6lBSwSiUMbh7QubujZsP6Jme9SGuAkQg1KtvfgldReJwmR+ALNcAwDFa+lYl0RRrWRjRDJfeSSCTpdCCjEXB3t3IECzV6pEhVCHyBk3NbDpw5jIWF0EREwJCfD/caNVGtXVuHth0KaMhk0CUlgVWpkM1UQYM3A+BdpYS3giIdorMVaGemZ8RAOh238czDEywRPNwdH0+jUgmNRAKNRAKvxi/j/gsBaFizCnwaVXeqLbaQq9BCozfipTpcVlGtM2L9hVhMe/t1VPbk5JPb3w+H9NSzDrXOCHc3Bl4eT/a0gs7AwsONeab7TKHRo7qDemBtOykSZQmYHdpA4MAptQYk5hZZrUHnJdlo9mJNNDDTfWNCxqBIp8PhfsJt9ESEB6kytGhse+eEOR6kyhDwPLfmlxUKjR4ebm42dbAlkqVFqFHZU+BwJOUVQWMwwrdRDYdoqHVGrD0bjZk9fQTffZgqwwu1q6Cet2OZGbXOiNuJUnT1EQY69UYWRpZM89YelFoD3BnG4T6wBBFh+5UE9Gv+PJ6vVcVUtv5CLIa2aWxa9w7cTMLCYxGIW97HpPuazDmJSh5uiF7WBwBnny06Ho4xwU3xSr2SAOynu24j6IWa+Kq3r8N8Zck1qF7ZA9UqCXcDsSxBZ2St+odlCQduJWNQ65fK3Be2UKjRw2Ak1CmWm9RWryFblY1K59JQvbIHmtQTBpvD0+Wo4uWOV+t7i5ED4LgDJ1fp4V3ZA+4O6qDEvCLU8/YS6AS5Wo/sQg3eaChcb5RaAyp5uAnmYEKuEm4MI2iTkSXkKDR4rqbtxEBFQak1oKqnu0M6N0+pdXi+icGWA/dsnoFzc8P6K8nwWXAGSq0wU/HnvTQk5hUJynIVWhy5n2b6vUpgAOp+Nhrvq4LQL/U5gfO270YSmsw5CY2+JNLd68cr+PqPxwKaqfkqrDwdCZYleHfujHqTJkF+7DjSuwQjLLg7Evr1g+LcOdQdOwbf3cjE4O23EJ5ekonYdDEOH2y+jvspBQK6chV3DkpbozK8XnsVOWvXglWrkaPQIGXXPugSElx8riwAACAASURBVNBg9iyT8Tdk+y18tveugIbByAIAQhPzrbru5KNMXMhhUfez0Xjt7Fm8cvwYXtqxA9+eicV3IdH4OzbXqp2XonOs6PQ+V4C3PbvgjevX8OKG9XD39sb2qwnot+ka7iSVfFeh0WPd2WjoDKzg/eMPM7D4rwh4d+yIFzdtgi4hAXHdeyBt4iRk1H4etYcONdXdcz3Rakx0BhZN5pzErpspqPXRh3j11Em8fPAgToxdjL5tplk5b03mnMS2v4VRJoORRWq+Ch7166Nm//5otHABGu/eBe2SNQguCECGTij+E34Ow8Rf7gnKkqVF6LDyAjLlani9/DLqjBiBF9f/iHepDd49FCOou+DoY7RfeUFQFpqYj0Fbb0JvZFGp6SuoN2ECGsyYjnZZr2JWtTaoN34cvDsFg2EYaPRGrDsXA62hpB8ScpXouOoitli0rc/6K/D7JoQLUHTtiurdu+FqrhFtV1zAxahsU73IzEJ8tvcuFh4LF7x/O0GKS1HcuDNubvBs0AAedeviVmIBjKwwqDNk+02cfJQpKIvJViA6SwHG3R2VXnkFVYKCMOW2Aq32xgq2lL61/DyafXtW8O79lALBWCflFeGX28mm3+UqPYIWhQjkO6dQA6lSK6CjM7DYdS0RemOJ7J16nIl+m67h2IMMqzasO1cyXoyXF4oYD4zYfx9pBSpB3dF77qDvxmtgi/uBQFBqDXh13ikr3ROamI8bcXlwr1ED1dq3R8333oN3p2CckORCrhaeeQxefRHLTwrPDRqMLIqKdZx7rVqo8uabqNahA7ofy8CIPXcEdYfvvI0h22+Z+OJx8HaKoP8YLy+4eXlh+v8e4LX5pwV1FRo9un53CQ9SZYLy+Fwlfr6VDHdvb1Rr0wZ1R41C9e7dMGLPHfT68Yqg7rfHI/DjeaHsJ+UV4aEFzbQCFWb8+kCgG95afh6d1lwy/b79SgK2/h2P/TeTTGWf7gpF03mnYA/HHqSb9CmPh6kyFGqEZdfj8rDlsnUEesvleMTlKEy/ExGWnpAgKqtQUI8/MG4PiXlF+DtGqF/9vjmDvhuvCsoGbb2JA7eSBWV/hKUhIVcpKDsUmgJJhpCXkIgs3EqQWn37jQWnMf1/DwRlxx6ko8mck8gp1AjK7yblW8kmwBnkYu3cfCkOVyzaFZejxIpTkVb1L0fnoMmck1Z9eC+lAEHfnsXZiCxB+bwjj9F97WWrb/J0LelP+uUe+m26BrWuRH+wLGHs/rsYvP2WoK60SIuHqXKcl2QLyg+FpmLATzesyj/eegODtt0UlMVmK/DB5utYeSrKisctl+ORZKEPNHojvv79EXIVQl0V9O1ZdPveup2d1lzE9F8fWJV3+e6yVb90/f4yev941apuTLZCYHvw2H09ETuvJWLnVeHZ1v6br6P/putW9c9LstF88VmBbgaABUfDMWrPHcRmKwTlH/50A74Lz1jRWXQsHE3mnLQqD1wUIpj7AGfYnnqcaVVXptIhU64WlKXmq7HydBTGHSixhx6lyfHj+VhMO1zSh6tOc2OltmiH1kwPSTIK8fOtFHxxSLjeX43Nw08iumLglhvoLjJ+ANBu5QX032zdn4v/ioDvwjNW62lIRBYWHY/A92ejBeVyFee8WOKvhxlWeu7I/TQrmwkAWi89j5ZLz1nR6LfpGrqK8N934zX0WPu3aLvM1xnzIwSSjEKrsVFo9Gi+5CxWnY60ovN3TK6VPgCAbt9fxqBtwjk7ZPst9PzhilXdwEUhmGRhm3Vf+7dVm74LiUb7lRetdJ7lGKh0Buy4kmBVLskoRIiFjsqUq630ulytR+CiEIFNAXBjbin7j9PkaL3sPP4ISxOUTzt8X3SeDNtxCx1XXbQqF8Oz6cAB+L24sTKV8PKIGf97iL4bhEps7L47mP7rQ+QXCetmFWqQLhMKGr+QW9a1xNTD97Ht7wRIMjnBqz/1C7x48CCOvNoZ96o9D8/nX0Ddz8ei7pgxJuEsMOM1IoNTqDkWirz72stgiQAGaDRvHvQpqUgaNAjn3h0Mxfof4N21K7y7dSuVN573ny7HWT2bfPAeJvwcBgBgPDxQ+Y03wLi7Q6HhjESjUSiw0399gNEWhqItRBX3hbnBu+VyPDZcjMNvYcLLMKYeuo8915MAAN6dgtH02FHU+fRTHH6jB75sPwFulUqiEZsvce3geeT+ry9uIzdejKcnqrZsgc151aBz9xQs7LyiWXlauMiuPB2FTmsuWSnFX++kQmdkcTrceuGwxM+3kpEh1+C4hUMgXjcFmXLht2b//hChSflIL1Bb1bc0in66HI8NF2Jx8HaKqSxDxtG7EZ8nqBufKzQcAG4xA4D7KSWGNN+nyVJh/cHbb2H0XuG4P06TY+iOW1hrtrAQEW4l5GPyQaHy7PnDFSvD3lLJAYBKJ1xg5Co9Bvx0A18evm8q+2TXbcw/Em4yyh6myaDQGrDxYqypTpsVF9BqmfCij/03k7D0hAT7b5YYw3E5nBEcb2EM30+RYcOFWEHZ6cdZuBaXhx/PC8vDkrmgC2tmRKYVqGFkCSceCuVg0LabGLbztqAsIVeJLw7dx8z/PRSUpxWoseOq8BKVyQfvIWBRCMRwL0XoEPG6iLUwaucdeYz5R4QOOgArJxbgsglJUhW+CxHOlQ82XceCo+EOOSp7byRZ9VnX7y9bGTHzjoTjz/vpVrJrDt7I0pvppWtx1vUtAzSp+Sp8efgBppgZYEaW0H/zdYzaHSp4d/jO21h9RthercGI1WeiMHBricGeq9Ri17VEfLJT+P4rc09h2A7hGK8/H2tl0HX7/jJGWnwbAGKyhbIYmpSPhUeF4zXzt4dWxvncPx/jXYu1bvyBMAyxcFR4WI734VBOJ8fmlHxfazBi4NabGLtPOPfVOiNen38aP1iMK8AZRSMs2jVu/11sv5KA1HyhXjsTzhk/95KFsss799ctxvbg7RQkiOiytGK6lpeY3Cuemwa2xBg3FOv/JAsdxwd/IzOFxmNMsSOSki8M3NxJKrAKiuYWB40kmUIHKb9Ih9VnovDpbqFcnHyUiV/vpmKliCGbJWKcp+arceR+ulU5ABSoxC+9skTPH66g78ZrVuXa4rllFJnTlnYRAKw4HQm5Wo80i7UqMY+TH0un/7GI0wgA+24mi5YDnMNmjsm/3MOkX+4hRyHsmw6rLqL9SqERqy8e8yJtyZrCy4HGIFxnADh0/Y2jG9DuJhcgIc9aTnnE5Sityn4pXsfNZRUokUuZxfh2XH0RbVcIA8CZcjW+OHQfE38JE5TzTmqBhX2sMwq/VdZ8fEhEFprOO4W4HKXVHHx3w1WrseHtjL8eWttUI3eHigYeAOu5afm7Oc5ZBFzEcLk4GSE1s+9DE/Px6rxTpnUd4HTa8lORVjbguxuuYvwBYV/33XDNSq8XFNP/65FQ5/I2rzmii/WNpe47asOmvBEvFZ2fYnhmHTh+N5n5BOMN9SILozC92Mi1NGzEwKd3LT1vS/CPzaP7nkHNsCfgPXz31nC8tG0rGsycCcbLy8Sjm9kWOP59ywlkLljVOnTACxvWg6lSFd56Fc42fou7mMPOOTO2eEI5cyEL3zeWr9xNLhCpLQ6+feY6gjfA1DprBWoOryZN0PDrr7DPvw+UXsJLGVhT/1mX2WqhwWz89BYKkgfvIFkqSlZkvGyBT4+LLYKOgO93R77FZ2IMRutvOfJ5fneBuWyLzSNb4KNq5gannWniNPgFxlyZmrLSxYsw31f25jPvHBaYzSmTnFvVtqZlTxeYjzn/P0e2S/AGhmWUUgwhEdaLkmWGrTS+nAW/BUVvENJQFMteRY43OSD7BHG9ZA6xAA0/zuZGO28oWTq+YuDnmHkE292k36z1yU2LrNcP52NKNejKAkvjq7zgd0KbzyO+3eHpQkOJz1oeCk2BI+ApWureEn1nwUtx3zoqXlobfWHKzJmV8eNuJWe8nrdghmysheLf4/61vMSI/6ZaJ+TTaNI/T38rK99HjvJSYscIyx21mcqCjGIdaWk/WAb+zPljRMrEdIxlsNoczqyL5YUN88QKlrvNgBJ9bLnzo7R2W8P5Rp4tXpfuFe8gs2dr8jLiiA3+JGGSETN2+YCgecDcZHPo7Q+OVCTZY9vOsIbJLnsCffNM3kIJiBtxtgx1cqIz7RkKJmVvUvIlL/DOnOWkEXOOHDFeAKDGO++gxjvvoFtxKnVGdftnTUiEN3twxmmxiSe0JolNenttNBgJ/NZyMYcHKJkwlkcXWCf6j1/8yjr3eJF1pNv5BdLc4ODfc0QxiimvssD8fcvoYXlRWoaHb7+pzXY+zbNpHiE0Of7FRErrCjED1xz89wkEECPgrTSUyF3ZBsJgL7hUEUNigzUjS2Xm2xIOBS8c0Eti/SE2L53pF15niH33CdippcKew17W9/k+Mn9c0m5hXX7uuTuoPGzVsrfOOKtHLaP//G/mxjkvH5ZftBUcKHFsHPi+yQm0UW79UfHypwAxHkuTNWeCjRUNx5xMaw9OLPjLwxGD+Uk2lafNraElZ93I4nlZYC/AXcJE2b8BAKDS12xLVJTuJCLTGu7M901z3k7DTbJexhRWyRja0nNkZoOUz4YsDc9sBk5sIPR2DHVnhEesM8Ui7ubDY7SxUJRMJvMMnG0BmTu2Kf76dpDjzFrAueiLkB9brzg0SZ6QYcMvKub9bzQtJjbeEYkqW9bl6Vr2k83FVwSOZAhKfd+JaC/fJndBIID7t6yLKiviFNpDWY1iR1Daoso/K5cjamGE1Fl5BaHj5pZW1WbfCnl13Lgx2pA7R2Ev0v0kInk8KjKCystOaaLnSCRTNMAjooyc6Rc+GGjusJSMm8NknIaYni3veNoKbPIwbw5f11IfONt2W05QRQWRrg88LjpvefrmTn1pzjjHo7j+dySAVyKftvpLPJj7NJwgS5TYICW8lBYcsmUjONNfwvccDzo68vf+xDJwYsFYvoml6bLy6mhnYKXPnZjuthxusUy0GEKWfIIF4152/INlhGXw1R7syYZ5s8tr15dG35nxN+dZTBbNYT7m5bUhS8Mz68CJrUl6g/hCxQu5M1kKsUXT/H0xZaa3sVCUZNvMaaH4fetRy6rjhcIGjt/kaIvPsgyeLYF1ZJKIZbTKqtzNwdMwH3N7i7L5+OltbKExiih3wLlFttwOFL9AOdC/YguuM9kcU12Rxcw5o9TMOCqDB1eachbdHlr8r9HCkS9LJqhk3nH/ejZSQV03G2JLnb0oJs8PEYGKaznCElsG/klEnm3B2a1MYgtPRdEuDY5kuh0x+MWMTjE2neGdl0Nhttk+v46gtCyHWBBSbE44k5Wz1W4xPWdLr4rpHsD2XLY1trZ2TpgyfA7OCVVNBTdvLT5vfi6VB6+jxIeNsQ7s2QkOitW1pC0mP1x97t+yRvYrEs7OEVu6wV4w1Znv24Ijc1fM6BYbS76sNDXqjAyUhtLWOp5fy7aVRf5sldvrt6JGdZBV1/kboE1dzPD/lM4sPx8czeDbSsbwMG+X3omt5Y72rTO7sMR4sheYNx8Wnn1H+8YZlEvNMAzTm2GYaIZh4hiGmSPynGEYZkPx80cMw7R0lDbfWQUqnWkxMx9I80N+vIznKbWiC5/5pRt8x6cVqKwEw3wh5R8VaY0mXviFgkCCm9X47yu1BhP90pzJ3qH5ePXMfdFnlrcKiZXztNV6VvQ2MaDkoLblO/lFOtMtluawvCyGh/ltPnw/KDUGU/t5ujoDK9r3d5LyRcvNx89oRoMHPzY5Cg1SpMLD5gBwI05qOjPFj5uBJdPlJ0DJNpssuUYw+Xie02Vq0X3n5t/j39MbWFFlKXa7kqDPir9VqNGXyIYZHfPv88NCBCs5UuuMoucMBXJRTFeu1pvKeToaPSsqW4LLYExyZTT1v7mi1YocFFdorOVP7CwPPwfNjWTLNuYX6UBEpt/zlDqry4bMb30zd44tgzhGlsCyhKIzr+DF271M75jLPiuikM3B9ze3hZIrkyp1ov1ozicvJ0Vag+lskTl9lU5E5vJVpjrm28OisxRWdSMzC0Vl0fJGPB73U2VWt/rlKrRWFwoA3FkLsb6wvGGSr2sJ8wtLeDL3Ugqs9BR/oyjfjDVnogW3pwLczYwsS4LAHc8HL59Ksz421+di/Wauc/i6BpZM48yXFar1yJJbXzghJutifaUUGV/+tkvzgAh/uZK5w87LkfkcyhA50H4zvuRMnvkcFdNdJx5lILX4wg7zfjPXUyZ9W6gVXHhk7jzfMDuEz+vcq7G5gls/+fmXLC0S6Da+PEOmtrqhESg5a8Pj5cd+ePF2L9xLyRcNQFyNzTPpI76dWgOLR2kys7pc+aN0uaCtPLU7SQWit/6ZX4DE132UJjddTAaUjJlUqRPYF7y+j8goFF23rsZaX/QECMfTHGEiZ9SPP8wQXVNPPsq06Cvu/7cSpCZezGXt1zvC8478q7cSpAJZ5z91I14qOnZnwrNE+TkTniWqK8V0xMlHmaJr8eXoHKt1QlqkM81P/v2CIr1p7pRkaa3XId4u4rnNKtSInlU2v53WHJZyat5syws4eLmMyCi0yMhwSMlXWV3eAkBwyyw/LwtUQjnjyV2PyxPtY57/18/eQ8/bJTxbXu7F426S8PIevq/TzNYlS1yLLRlHXq6yCjWiN6ICQhk3H5vHadb1z0myTXrZXAeJ8W9+OQjP6t3kApMM8/yHJReYdBtfLymvSFTuxOy66/FSk/3Az/PsQq1Jt5rjUZrMyiYp1OhN9xyYI0OmFu1jc31jC2V24BiGcQewGUAfAP4AhjIM429RrQ+A14t/xgHY4ih9voP6b76O3dcTAQivgTW/ZpOv+96Ga/j9nvCqTgAIXn3JZKzwdT/dFWp1tfSbS86abqHiD7IP3XELmy5ytyTyi5beSIIrgfkBGrPvrulGRX48Ru+5g73F/PMYfDkXLx+4ZnUtLwD4LjwjekOm78IzJkHhBSOrUIu31/1tVRfgbqYyX5x4fqb9+gBLT0is6rdadl7w5wF4tFlxwaQg+HYuPBaBRccjBGWrz0Rh7bloq/c/3noTJ0SuCe646qLp2myeRufvLpluGjL8X3vnHV5VsfXh36T3HtIr6SEkISFACCkkgYTQQbrSEQERBRTs/arYC3YvlmtBrn7eq15BqSLSBJUOoZcASYD0nvn+2OXs2bMjSJEcnPd5zpOcOfvsM3vWzJq1pqxRK7+Urmf6R7/gQznaldZwG7p4g/q/IuvRb29kQikrZfHWukOY/QnvSGcuWq0qS0XxPvf9fjzy313ctQUv/shFDEp7cqWqlKimbn66RYoKp+1EUx773rSEVE5/6D+78PJKqR4pn209et4wwl3MA9+pxqryvO//fBTzPv+NSfv1+AWMfpuPXhe28FvVQFTK/KficvW3tA5P9qI13PcTHl7BOf/R93/HGfwZT6/GqQt1qpNeWtWAaXK0J6WMi15ej//+XoJ6WXHvKanEwFfXM8ot+9k1qhJUnm3xmoN4So4yqOj6l1YewKNf70b1ZzHouGqkOpKYtWgNN9P3xfaTuEsXMRKQwkOv3CPVRyUHb647hMe/4dtPl8e+V3WMIrO9p6swVo5eqO3AUx//gTN4shatwRp5k7W2fugjfQJSqOWlW49z6dnPrjE0SIcu3qAGS1Ge+XBZjWEo8f6vrFej4mnLfdBrP3FHheRoZKEw5u1Nauh7pV0vWr6Pazspj/+A4+dq1eUlza2Ui9o59/PfsOFgOePMKJEulbxcqG1Cn+elMtLqAaNy6/nUKtV4VK6tbWzBkMU/MWnNrRQDXpWi+mnllPjICm7wq+sTKzmjovPDKzjDO+/5dWhsbmWMESXinNZBVUJoa5/FKMz+6Lc3qsF7tHVLqyuV8v90y3Hc++UO5r4trRS9NeHDlb61saUVWZp2rh3UHPPOJs4oefCrXXh+hSmUtnL1K6uK8dBXJpkr16/YfYYL1Q9IdVQbATJqUxd0XDUSdy/7nYluq9xn7ue/4c21h7g8DtTUaSX1m99LMFUTWU6p1//57ZQasVlL7nNr1b5WubairglFL5siPSq/2djSytgCqtN4ogIDX+MjQ9787mbDgYDRb280HOwZ9voGzkCf/cl2/G5gJM/8eBsTwEe5y48HyjD+n5I+1w4O6Y9OUsr2nn/vwDPfmfpyRebPf78fC7/4nfvd6R/9Yhg5dubH2/C5gZ4a8/Ymw/5VW48UJvxzC7YckSOPtpgG+vJlu0cpqn1nqtTo5MpzZDy9mjuioM8L61DX2KKW8ZnKBmQ9s4b73bzn1xk6RkMXb2Dau7adFr70o+F3bnlvMz7SHBui5G/joXMY9voG7vpRmiOpGuVgU/VNrch4mm/b85f9zkUEVvK/p6QSCct+wk2rTI5T7nNrDZ2I4W/8zNidisxfXlWM2qYWw5mqce9uUr+jbYPaiKjaujz67Y1q+WgHnQa8up7rD2d+vE09Gkerc3OfW8s5O2Pf2WSyjeXPFn6xA09+K0WCVerY+uIyVfco10l1mm0HgGTX6QdVxr+3Gd/I9Ul53uqGZvR7iY+wOfyNn/Fv2RdRnu2HPWe5Y6oAKeKqPkgWINlD+mNk9FxJEJM0AMWU0kMAQAj5FMAgAFrrZhCAD6hU4hsJIW6EED9K6R/Gby+pqIe1ZqRn46FyJAS4qp2QQmV9E1btOctELNp+7AIC3e2Z0QFA8vKrGppxptJ0350nK1QDDZAq1bZj59FKKROJb31xKfLiOjAK/OSFOpy6UIflu04zSnPT4XNIP3aeUWifbjmOzChvzJWNalCAkCZ8uf9/sComKD5bAytn09kTH+2QZoz2llTBytmUv0921qC5leJfv62FhXw+eGlVA/adls6B0Tt+ZysbsKekEjtPVmCDpjL+dLAcB0ur8dZa9nyYbUfPw9HGCi+tZBXpsXO12FNSiW93mPK4eu9Z+ewo0yjeRxuPYViXQO7MsNKqBqzbX8qdfXS4rAY/HSxHvSYS0E/FZfBzteOcjfM1jfhWF/L1TGW9dM6axiA4cLYaFbVN+GZHCSPr305cQElFHX48UKaG9wWkELO7TlWoYXkVjp+rw4Gz1Viy4Yia9vGmY5iZE8F1NEfLarBsKztwUFxahdrGZpRVm2Sy6VA5EgPd8M56U7k3NLficHkN9pZUYanmHsu2HUd+nA8mv28682bzkXM4XVGPj3WR4naerMCOkxV4bbVpQOLr30swO7cKN79rcvq2H7uA0xX1eHMdO3Cx7dh57DtdxZTBhoPlOFJWg0maowZKKupxvqZRVUwK+89Uc6Nu64vLcKiMNWx3napkzvv6fvcZbD92nhmY2XL4HFNPTpyv4xSsUh+1xwJ88PMRDOjsjzc0oeaXbDiCyS3SmXRWTrsAao0zLcDS3bU4daEOr689CCs5ZtB/i3cgffcp/HriPKycTUbHF/tKUdVYg/KzplmJb3ecxtRe4XhHdyTA/tNV2Fp/njEKd5ysQElFHV783pTX2sYWnLxQh1V72fMXj5TVYKN1OWbqlPzJC3X4QFMPlWt/OXpeHTDSls3hshoupPyxczX45eg53K4ZsDh5oQ4VdU1Yogt9fKhUurd+oOfE+Vq8s5595uKz1ZyzcvJCHZpaWpkIt5V1TXh7HatvdpysYMIuh3s7cQZQSUUdZ2QfLWfr5enKepyprMe9unpSXt2gDv6ZnqEOGw6WYfYnpnOj9p6uQmV9E/KeNzl9pVUNqKxvwiuaOtZKgTNVDfjqV1PY97LqBmw4WI6DZ/XHVpxnZqyk367FW7oyqKxvYgyesuoGlFY1oOsTpiMz6ptaUVnfhIU6g/tURR32nq7Ckg3sM56tko49UQxfQNYzZTUYpznyorqhGSUVdYa6ubK+Cct3nsZLuqM3jpTXYNPhc8xxKRV1Tdh5sgJbjpxj6t3PB8twprIea/eV4knNOWqHy2pw4EwV/r2NDZ9/pKwGVfVNsn7xAgBYOh7A2hMrsbfKTnLwHEtVo2Vb+Wm8sWUnnl2xT23HAPDVAUv8cuQ8ztSUA5AOKD5xrhYHS6uxZl8pPtlsat9KH/6R7vkPlkr1Wu2zZarqm7B6Xykz8NfUQlFe3YAvtp3EE9+ajg+4UNuE3acq8epqtgyPlkv654Cuzpw4V4ej5WVYtJwdCD1YWs0Zk6VVDdh65Bxz5Axg6hP3n6li6lp5dQN2n6rEhxuPMNdX1jdh18lKbD9+npHp+dpGHC6rwdKtx5njAmobW3CotJrr42samvHzwXJuJqe0qgE/7D7DhYA/V9OIDQfLmaMVztc24mh5DVcvTlfW47udJXhW0+9WNTTj5IU65ky4UxX1ePzr3cwgyfJdp7F2HzvreeBsFSZojk5qbGnFhdpGPKrTd8Vnq7HjZIUaml6hvLoRra1V2Hemiqs3pyvqsfd0JbYeYXVi8dlq7DhRgd9PXmCOfDl+rg5nK+vxue6MsGPnalFe04jHdXmqrG/CdztPM0cuVTc0YeuRc/ipuJy7R2RzK0AprJxNuuPTXc04U1mPX4+fh5WzyXb9v/3NIJBmOL8+aOobW2gzquqasHTrcU6+py7UYe3+s1j0HVtnS6sasGL3aW5W8nRFPY6U1zBH/wDA4XL+HNHqhmb8cvQcpn/E9oenKurxjO5YmKPltdh5sgJHNLPeZdWN2H2qkjnX70xlAw6X1eB/O0227O5TFdh5soLrS4+fq+WOaSitktqRMoEBSHWxrrEFL+jORj1aXos1+85igaZfWre/FGcr67mzBovPVoOAcEcNnLxQpx4fZAT5MxFemC8SMhxAAaV0ivz+ZgDdKKWzNNd8DeApSul6+f1KAPdQSrca3G8apFk62PhGpPiNf/Gy8mUOLP1yLAhpwqSF4Zd9D0oJag7cC9py8aiVNwIW5K+PDtceIeTKohlZWpBL3it0pb/VHihcGg9Lm1KsXTjpiu7TeK4nGs4MvPnMCwAAIABJREFUuKzvWlmQi0aWvF7YWVswAyjmio2lBbd098/UdWtLctF9GdeCP5PHvxJHG0vuuJ7L4XL1duFnibC0LbnidttUGY/6kzdf0T20tNWWba0smIGoG5lrqc/+TJ9zNfunq1XfrzV/9pk/+/w2WFiXX5GtCQANpbloLMu/onvcyDjZWhkuxbxaHH26/y+U0lR9+pXsgTPakaevWpdyjZRI6VuU0lRKaWqML++UjOkWjDfG8Vvo/jmhK5c2OMkf/7ujF5f+2hj+++N7hHBpY7sFc5sTByT6Y9n0Hty1H0/phlg/FyYtK8qb26w9OMkf38zOAAC01gegpT4QNYfmYITfi3i/z2eoOTSHeU0OfQ2v9PqIS58Q8grGBb6MmuIFqvPmZGuFr2/PwKq5WVz+HhkYjzXzshHZwYlJz4/zwUqD6yf2DMW/b0vn0t+flIbbsjsyaUlBbuga6s49u9F95+ZH4fs7M7n0O/Oi4GhjyaQNSQ7A+5PSmLRWCnx9ewb3/UFJ/vjXlG5c+scGaXF+Lvh2di8kBrIBZLqHe2D7A7xyenJIAvol+DJpmVHeXN4A4NmbErm0eX2iuLQRqYH48e4cLv2tm1PQM8KTSesT54O7C6KZNH9Xe8NyuKVHCGb3jmDSMiK88OigeCatg7MtvpvDt42FhTF4ZlhnJs3Z1grPGTzXx1P5sjW6blLPMC4tP86HeW9rdXEVZHTNu+M5XYZ4fxf0junApLlYuyPCLY5rR1PDXsPE0FeZtMDaB/BCzw+5a7tZP4EQjGTu2yPc01AfPD8iEbNyWDn4u9njv7N4mb0xLoVLm5HdEf+cyOs0o/K9KSUQa+Zlc+nz+kQZ6oIRqYF4ckgCk+bpaIuvZvbkrh3VNQirDe5tJPu3b+FlMbRLANeujfT6paCv1wDwju43G1tasXgsq99bWik+m9ad++7Cwhi8eTNb9k0tFI8P7sSkWVoQvDeBfzajfujWTN5Ayon25tJu1z1LSyvl2h0AQx1j9LsvjEzEFzN4ff3v23ogzMuRSesb74PhKYHctcvnZGLRcDYPTa0UX85IN7y3UZkqei0rin3mzoFu+G5OL9yaxZaPm4O1oR57ckgC1s3PAW2xR0tdkNoGC92fxaLuH3BtM8fpafyj2xIu/cm0JSjyeBb1p0zttluYB76ZnYFhXdgy6BHuiU335nJ5eWFkIif/5laKz6f34GTWSim+mtkTM3R9pLOtlaG+f3RQvOHzvzI6GR8YyP5TgzKf3zcay+fwferCwhismZeNogQ/Jt3P1a6N+pxi+Px9432w4s5MOOjacVqYh+H19xTEYJtBPzq1V5ihnrorn+8f+yX4qXaSlnv7xeBLg7popHv0dAvz4NIGJfkz7xtbWvHiyCTuuo+ndMOOh/tw6XfmReHr2zPwYH/9riEpT3seLeDShyQHYP09OXh+BKvLvZxsDdvZ/L7RXD9HKQyvHdU1CD8t6M2lBzr7IdKg/5sT8yaWFi3j0scGvoRHUt7j0hOdh2PFnZmGdfa9CamGfcXXt2dgs0E9eXl0MtbO56836uOm9goz7J+M+v+HBsRxfcygJH+MTgti0hKD3GCjsylc7a2x/h6+nd5fFGuYJ73f0djSis8N7IGhyQGYq6vnXk62+HY2b4NNyQgzbD9GeVC4EgfuBABtyQQC0B8tfinXcFhbWnCGxqSeoSjoxCqk1BB35MR0wAO6RnRXfjTnVAFAUWc/vDw6mUmb1TsSI1NZAS/sF4vFOmfvjtwIpIayiqBrqDvSI7zw8AD29+f3jeY6uTvyohDv7wpfFzuAWgCtNmht8MX9+b3RxS8OrQ2+zGtOViaywxPhiEAm/a6sLNyTmwPabHJChqcEolOAK8K9WScNAManhyLUyxFPDmXLc2qvcHT0dsLQLgFM+pzcKKSEsE4ZIDlm+g7hjtxIPD+CVXq3Zkn3zYjwYn8vMxyRPrwBd3vvCCzox1bQ+4piOUOgZ4QnOgW4Ymov1il4eEA8eup+q1ekF9IjvDhF+eKoJMT5u+AFnaK+pUco3B1tuLyN6RaM+X1jmLRJPUO5vAGSDPQdxfSsjviHrtzvzI9CkAd7kLm9tSX6xPvi8cHstWO7h+C2LNYgKOrsJ8uaNc4eGhCPkWnBbF4zQnFLj1AmrU+cD2J8XRDgZs+kT8sMx02prGEzPCUIw1J4Yye9I1veADAsJRADE9lOcW6fKHg5seWq73CHpQTiPp38f7k/j3m/sJCVAQDkxvpgy33sdbf3jsBjOiM8wN0BbnaOXPuanZmJKd3SmbRbUrojLyIJU3Xpzw0pxDu3sI7LtMxwTh8AUic9ry/rdOfF+iBBN2jgbGeFgk6+GJLMtr9ZvSOQE806oSGeDhiWEsgNCDwwIA6hOiMdAGZkRyDc24kzVB/oH4cx3dg6UtDJF4lBbpicwbarhYWxnAMAwFD2+XE+CPVk6/RzNyXiLZ2sXxiZhCAPtt6tm892msum90BnXVnNyWOfu0uwG/J0AwGZUd7op9NPWVHe6BbuiVydUz8pIwx94325a8d1ZwfzRqcFoXcM+zt21hYo6OTH6cI786O48n7VYMBwbh+2bnRwtsWIrmz/4+9qx+mYGF9nrv8DgCHJgegSzOprFzsrpIR4cLpnQWEsN9CUG9MB0b7OuEnXB2ZHeSM52J27NwB0C/fkjN1pmeEI8nDAAl1bvTUzHDG+LpjWi3Xg5uZHoVOAKzeYNrJrEII9HeBsZwVKrdU2+ECf3iiITsa8nCymbd7dOxv9Y1KYNH+HMAyITcGM9F4ANemf0WnBiPd35QYiR3YNgo+LHeyt9QOJgVxb7N/ZD11DPTiZpYZ4IDHIDdN1957ZO4LT9wBwc/cQdArgI1EPSPRHpkH/0j3cE566PmpGdkdE+zpzOnZKr3CEejlyspjbJ5qrz862kh7ycbHDMzon/vbekYjyceYGO27vHQkfFztuUHhyRhg8HG0wWOccTZXzo2d2biReGsXWoztyIxHv78rpngnpYUjW1cWCeF9uQHBc92A8pLPHXjcYKHtpFGsHjuwahME6XQwA6RFecLaz5gz+2bkR6BTgikkZxoOU9jaWmJ0bqftOJALdHbg8T84IQ5dgd4TodOitmeHIjWWvTe/oiS7B7twzzsmLQoCbPdwc2IiT/q72hv3f5G7piPWK5tIX5PbG0E5d8c30m9jre3ZElI+zYZ3tHePDySs3pgM6Bbiig4sdd/3ARH+EePL1YVhKIOdszciOQGKQG5Pm5WTDlQsATEgPxQzd4OmsnAiu/7gtqyM3GPDy6GQEujtwNtzEnmGcXp+TF4UcXZ+SGekl6QWdDTW3bzQm9Axl0sakBSHO34UbsJ+e3REejjZwtWdleLPBJJPCleyB2wIgkhASBuAkgFEAxuiu+Q+AWfL+uG4AKi62/01BiTBlaUGQE+2NjjrnZFCSP+bJnaGyb8LdwRojUqUOQMvQ5AAMkI1L5VpnWyvcltMR3s626m9F+zhjcHIAnGyt1OUQtlYWmJAeiogOrPMR6+eCB/tLsxv1mqUT47oHI97fhdlzNi0zXK3gFBS3DV6IrChv3J/Ykdscmh/no+ZVKgfp3kUJfsiO9uauH5Doz436AZJBMyTZdB9tAILEQFd0CZYahXISfVqoB/rE+8BVpwC6h3tgcobU+WrD9+bF+iAj0ouJSpUR4YXUEA8531KZJgS4Ij/OB3a6znFolwCkd/SChQVBnRy1jRCpcXs52XLP8oRsmCvl4Wxrhaxob87xyoryxiMDZblolob1S/BVOxzt7GhGhBc3azMhPRTpHT3lZzalD07yR4+ObKMb1iUQubEd5N+TnrkowQ994n1gZWmBBjnN1d4aI1ID4efKGrD5cT6YIncC2jCzBfG+6BHuyci7KMFPdegUuSUHuyEv1geWFmy47P6d/TjHdmRqEO7Kl9qMstzsppRAZER6GdarO3Qd0MjUIMzSzSAMTQ5QnTzl+fNiOyAzyhuOtlbq/tSRqUEo0M1m3poZjhk5Eeo+gtQQd0zNDIeHRqaS0xGCh/8r7QUY3yNE7WitNA88IT0UubE+KNfsOXx0UDyGzbUCIcDDz0pp/RJ81fqsLe8J6aEYnBTAPEdKiDsGJwfA1d6aCTwwrnswZ2QNSvLH4KQArhy7h3tgejZrvObH+agzNspvhXk5YnhKIBxsWJWcHOyGxwZJdb+6waS7pmaGw8WObasxvs4YkOivhoKva5LaVbi3IzIjveGsuz4jwgvTdPnwdbHD8JRATg8kBblx9SEtzEM1zJW2NjDRH3lxPiCEMFFLixL8EOPrrNZbRxtLDOkSwIRbf2lUElJC3FVdVdjJFznRHZjQ9klBbnh0EOuk+7naYUEBa6z2ivTCg7KRo9TBtFAP9O3kC2tL1hgrSvDjnO7JGWGYqTMGpvYKwwjZ0VHKq2uoO3JjJf1WLz9vTrQ3ijr7w9HWJMtRXYMwUmf090vw5YyL5GA3PDSAnTUv6uyHO3XXDUryV/OixdPRBk/Ig59anT8rJ4Izsjo42+JOg1kQFzsrzslReEp2Cus0exX7xPlw/bNCRqSkg7R9h6ejDfLjJF2gyCbQ3R4ZEV6qbj4wcTf2lFTB3toSqaHuat1t0Oj0pCA3BLqzfb2rvbU6+KsP5a8YXfoQ44qzrMgv1NNBddC07ZkQGBrsADA1k9fhtlYWnKMPADNz+H4/yseJG2wDJJvmXnmAS7tEc2iySddo5ZwTbVr9o18F1CuS7Q+CPRyYGSTt/QPc7BHvLw2C65dYp8qrbrTXJ2lmNZR0BxtLpHf04oz4jt6OnD0GSHVA6aOVdEKkwQS9A+XjYotbdYObrvbWmNQzTA3YBEi6T9ufDOsSiHHd2UGsaB9nTNUNMKSFeTArSJQy9nO1Q0EnX8OgHp6ONpijaU+K/nNzsEbPCC+1/WnbQmQHJ3U2UNvv9IzwgpVOT3VwtlWdckUmznZWKErwg6+rHXOP1BB3SQ/PWyZ9+Vlp72usnwsm6hwKJR+3aBwFfXn3iuQHFawsCKZn8XrCggAzcvh0G0sLQ30DAE8Pk9qs0r7dHayREcnbdx6ONvjHUHagIS3MA8O7BEp9jqY+DUj0R0QHJ5RrYkPkxfogK8qbCUQT7++i2nuKXksKckNOdAdYWhBVLwS4STpKq9cBqazvkfsfpe4rfW6Amz0TNCY3poOqQ5RnDfdyRGaUt2r3KjLMjvZGr0hv2FqxtrOWy3bgKKXNhJBZAJZDOmb+PUrpLkLIdPnzNwB8C6AfgGIAtQAmXur9lYKYkhGGhf34KUTtCIoitPHpoVyHCADPjUhUG5xSOAOT/DEjO4L5rdm5kSjq7CenSdcNTgow/H3tEs1aWUB94nzUWZQGTWd+r+b7jc2tOO/gigcmZMJT56gAwOtjuzANVzG0nxqWwBlggLTkwoj3J3ZllIzy3DnR3vjnxDQufWpmODcyBACfTjNNCyt56Rzoind0U9gdnG3xkWbZoqLwHh0Uz42cAcDTwzqrhpTSic/KieBGqAEwS0qUSn9vUSxG62acAHbpkfJs43uE4BGN0aes3w/1dGDyrHBfUSxn5AW42ePFUXxZLxreWTUwleeYnRuJaHkqX3HuR3YNYuqBgnZGSonIF+Rhjzdu5kcNX9MsEVPq1+tjU1TlrRgP7g7WhqP/T2tGWBX5LCiMMayHRvXq6eH8Uq9nb0pUn1951nHdQ5Atj1wrCvHeolhuZGmhzjDJiPTiZkb0M0NaOSq/a2NlgYdlp10x2twdrA0NosVjTeWq7UiV7wOmznF4SqBax5RrOzjbcjOlgFQO+joDsO1HQStzpXwe6B/LjY4DwJczTMtHlI7grj5RmGiwPPV/d/Ri2nyN7PDdX2R8b23dV575rvwobnYBAP7PYBnL0ltNz6Z0cg/0j4O3s9IRSfcsiPdV664i69Xzs9HB2U4N4e3paINBsgOtXDO/bzS3qsAoHz8v5JfpfDjZ9Gy1chnfUxhjuLrgtbF8W9Gv6gCA+4pMacrzzs6NVI0bRTeNTgtGH109fspgmaS2Lipo5a3mz6At62cQFH7RLMGpUwdUfDgHFQA262awFbY9kM8ZjwBw+B/91PqlhPiekB7KtB1FryQEuGLprT1gLy+/0zrhm+/LU50LRWd+PKU7M/BabduIBocmLJ+ZiSgfUx1Q6sbdBdGYnskbib8+mM8Z1/6udvhpQW81XYlAF+LpgK9m9oSbg2QkKnuLvpuTyQ04AsCeRwsM0/c+ZkrXOk37Hi/krj3wRCEz8KSw4k5+yTMgyULvqG17IB/umgEWpV/+7cE+cLYzmXSKzvJyssHKudmc/l2nW9qp2FGTeobh7oJok6yVgcnOfnh5VLL6jEr/+uPdOfBzNTlpioxeHpXMzZQDwMq52dy147oHM3pVuffme/O4mUcA2Lgwl5Pzbw/1kfMrLfKakxfJ2YPPjeCX6S032Nqh1W1a1szPbtOg/kW3/E2R1+29I5l+TFtHvr/LJHdFX751c4phn6xtr4pMJvYMY5aiKuX53sSu3AAfAMOtRfp8AKbomlE+Tlg+J9PQYd3/eCF3bqSSbqQ/9jxWYHgOpFavKPl/bHAn9O/sz11rtMTw06ndORtkdFqwugJBaQcudlaq3arYT12C3fCFRucq+uiZ4Z0RJa8YU/T6gwPiOPsEAJZpthwp187Nj1Z9Ce0zv6tZeqk86wsjk5hZRiX9tTFdOGdRz5XMwIFS+i0kJ02b9obmfwpg5uXcW9nYbWugMPUoQjNSrgA7iqbc10lTMIqy0I46WMqWoIv9xYtIEbq2sBVB6hs7BTB8xw9w+OQ0MHUydy+jig+AG5m/GPoGV99GfhrUsrv4alqlUWuXmjS0UfbK7+t/LzHIDb8dv8AYu4oxYG9zcVkriutS8qt0bPq81av3MP49bd4UJdzWXi2tAlPqAVM+crnbXcJeL+X7DtYXl7VRmSn7yvUdtRFK8KJLaV9/hPb5lXLVylwxivT7obQoxyfoR/0uFW9NZ6cEolDy8OuStr+nOHsdnNnOUqlj2r0f9QZpWoyct0tBqV92fzDKpqA4cG0pdX2bV+vIJdQnxSGxvYR2Zfh9g/JRBkq0RqVynaOsz5RwzNq2qNajK6ybCsoMf1uyuxyUMtX2I0q7u5T291fi4fjn8tNWH6StX4oT1LZutWB0k/JVLycbxqBRr7dhf9P9Z29E1LrBztqC+V1FlnZWlobGo/ZadRWNtaXuHlK6k62V6rwpeaS0bV3fVn+hTdcfkK7nz+oJbb6V/szBhn0eWytL1De1wt6GLRNT32V5SXVSawton0npKxxtLBnZqath7KyYOnOx/pV5PvmvPn+KveJgc3E56/kzv/9nsfkT8qtvw1ZpK+smHXpxfd2WXaL0tw7Ksy9ZIn/Cz6D9EU3Niu62brOsjeQCtK0/jJw3gJWlat9dQn9olI8GgzJvNrAvlPrVlv7S9hWKs3cp9cnIPm1L3nUXqaf6Jd1GXJEDdy0Z2y0Yh0prMEW358nP1Y7bvzO0SwDeWncI/XR7BCb1DOMOqxzVNRhHymqYpWAjUoOwel8psx6/X4Ivth87z02tx/g6q565Ql6sDxJ1S4w6yx71aN2ek39N6QanZffD7l/bGQfus2ndmXNwFObkReLFHw5wld/LyRZFCfxowPSsjvi/7Se59GR5yeQtuvW08/pG43RlPTdLlhXlzYWATQhwRb8EX3XpKiAte+gR7om5+v05/eOweHUxt1frpZFJ3LlJQ7oE4M11h9A/gR9x0a9/npHTEbtOVXB7EwrifRHgztaLsd2CsbekklsOFOLpCFsrC8zXjUqPSA1kQvkDQJCHA8K9HLllW2O7BTPHEQDAvL5RmLv0N/i4mhyCnBhvvPDDfm5EfnCSPxMKGJCWMQxJDjBcvqSfbXxhZBJeX3sQLhrj2MfFFtMywzFGd627gzXGdmPl/sm07vj69xLOsVpQGIPu4ewy0Tl5kdyI4ILCGHyhO05gUs9QbD58jtlIfFd+FJ7/fj+j1FND3Jk9qtN6dcTZygZuxky7NLgowU89g0XB1d4a9xfFok+cqWzdZENACZqgOHBvvpbCHaBta2WJp4YmqMu8FBb2i0Erpcx9/d3s4e9qh4cGssvbHuwfh8+28Gcdje8Rwi23fevmFOZYCQD4x9AEvLKqGF11a+/H9wjhjJrJvcLwY3EZt+T3g0lphrrj8SGd8Mx3e9ElhG1D9xfFIkG3l+HuvtGoqm9Gnm5vwYLCGO7w3sJOvkwYZkCaUX9uxX6m0ymI98XUXmGYlWPSi9Myw/HKqmK1g/R1tUOvSC/MyTNdM7FnKJ78di8z8v7qmGS4O7Aj8R9P7cY474A0E6rv+Gb1jsTsT7Zze0zuL4rlluCN7RbMLTV8Y1wKF0766WEJ+HK7B5I0+unhgXFIDXVHmkaWs3IiuOM0JqSHcrMwjw3uxLRlQJrd1+v9Dyen4ZjBwbGzciK4AbDsaG/cmhWOGVnsUtDbsjtyexYBaf9RtEGgmRdGJjLnPAFA/87+WLH7DMbq+rdOAa6I9nFmZisBadnvHbmRXPCIBwfE4eH/7IKHTrZZJ6JQcqGOk89t2RE4XF7L7UsJ93ZErC+7793P1Q6WFgR36/R8jJ8zhiQHcEtkv53dC1uOnOMMVjtrC24PHyAtSdTP6FpZEAR52GNOLtsfFnbyNWyjjw6KNzzwe9HwzrhQy54XN79vNBYt38cZfJ9P74HVe89yA2DKTLg+EFaUj5PhIFBmpDcWLd/HzTLkxvoY7u29r18s7v9qJzOIAUh7n0s+/x1Jwaze+WBSmnr2q8KIrkE4W9WgroZS+GhKGj7bcpwbdEkIcDUMDqRdFtknzgdvrTvE2AjJwW7YfuwC973CTuyzBrrb48R5/mDvDyal4Yc9ZwydmYJ4X9Q08hEIZ+dGory6kdsbbm1hATcHa8aOAqSAfB9uPMo5fPcXxSJYt4+yT7wP3lh7kOsLFhTGsP2t7MBNmruYs+cASYa/nuDLJd7fBdMywzEhPZT77NUxydhbwp9hvKAwxlB/3JbdESUX+DJ97qZEWFmy5TkiNRCr9p7l9kE/MaQTp9dvSgnkjl/Ij/PB49/sYewlbydb3NIjBKO6mtI6BbjC3cGaW9L59i2p+GTzMfhrtrosKIxBdcNOLmDfqK5B6KJr/wsKY1DT0MzYUDaWFhiU5M/prEcGxuPxb/Yg1IuV7bM3JeK99YfbdJC1XPYxAteS1NRUunUrd9LAjUN2tvR3zZrrmYt2T31TCywtyGXPbtwoKG30j0YdBcYsyZb+TlhzPXMhEAj+DKLdCq4mlFJQys7WNLe0woIQbhVJcyvlnFKzRdiaNwSEEMNjBG6QWiq4EbkWSyDMEeG4CQQCgUBweRBCuKVsRkv9hM0hMCf+3lMbAoFAIBAIBAKBQGBGCAdOIBAIBAKBQCAQCMwEsQfuelArb1p24DeSCwSCq0eT3NSsRVMTCMwG0W4FgquAsDVvCMQeuPaEaEwCwV+CMAAFAvNDtFuB4CogbM0bGrGE8nqweLH0EggE15Qti6WXQCAwH0S7FQiuAsLWvKERDtz1YOlS6SUQCK4pu5ZKL4FAYD6IdisQXAWErXlDIxw4gUAgEAgEAoFAIDAThAMnEAgEAoFAIBAIBGaCcOAEAoFAIBAIBAKBwEwQDpxAIBAIBAKBQCAQmAnt8hw4QkgVgH3XOx+Cy8ILQNn1zoTgshHyM1+E7MwbIT/zRsjPfBGyM29udPmFUEq99Ynt9Ry4fUaH1gnaP4SQrUJ25ouQn/kiZGfeCPmZN0J+5ouQnXnzd5WfWEIpEAgEAoFAIBAIBGaCcOAEAoFAIBAIBAKBwExorw7cW9c7A4LLRsjOvBHyM1+E7MwbIT/zRsjPfBGyM2/+lvJrl0FMBAKBQCAQCAQCgUDA015n4AQCgUAgEAgEAoFAoKNdOXCEkAJCyD5CSDEhZMH1zo8AIIQEEUJWE0L2EEJ2EULukNM9CCHfE0IOyH/dNd9ZKMtwHyGkryY9hRCyQ/7sZUIIuR7P9HeDEGJJCNlOCPlafi9kZyYQQtwIIcsIIXvlNthDyM98IITcKevNnYSQTwghdkJ+7RdCyHuEkLOEkJ2atKsmL0KILSHkMzl9EyEk9K98vhuZNmS3SNadvxNCviSEuGk+E7JrRxjJT/PZPEIIJYR4adL+9vJrNw4cIcQSwGsACgHEARhNCIm7vrkSAGgGMJdSGgugO4CZslwWAFhJKY0EsFJ+D/mzUQDiARQAWCzLFgBeBzANQKT8KvgrH+RvzB0A9mjeC9mZDy8B+I5SGgMgEZIchfzMAEJIAIDZAFIppZ0AWEKSj5Bf+2UJ+LK9mvKaDOA8pTQCwAsAnr5mT/L3Ywl42X0PoBOltDOA/QAWAkJ27ZQlMNBrhJAgAPkAjmnShPzQjhw4AGkAiimlhyiljQA+BTDoOufpbw+ltIRSuk3+vwqSARkASTbvy5e9D2Cw/P8gAJ9SShsopYcBFANII4T4AXChlP5MpY2XH2i+I7hGEEICARQBeEeTLGRnBhBCXABkAngXACiljZTSCxDyMyesANgTQqwAOAA4BSG/dguldB2Ac7rkqykv7b2WAcgVs6lXByPZUUpXUEqb5bcbAQTK/wvZtTPaaHuA5GzdDUAbsEPID+3LgQsAcFzz/oScJmgnyFPOyQA2AfChlJYAkpMHoIN8WVtyDJD/16cLri0vQlJ+rZo0ITvzIBxAKYB/EmkJ7DuEEEcI+ZkFlNKTAJ6FNHJcAqCCUroCQn7mxtWUl/od2bGoAOB5zXIu0DIJwP/k/4XszABCyEAAJymlv+k+EvJD+3LgjDxhESKznUAIcQLwbwBzKKWVf3SpQRr9g3TBNYIQ0h/AWUrpL5f6FYM0IbvrhxWALgBep5QmA6iBvHyrDYT82hHyXqlBAMIA+ANwJISM+6OvGKQJ+bV1noq/AAACl0lEQVRfLkdeQpbXAULIfZC2g/xLSTK4TMiuHUEIcQBwH4AHjT42SPvbya89OXAnAARp3gdCWm4iuM4QQqwhOW//opR+ISefkaerIf89K6e3JccTMC1f0KYLrh09AQwkhByBtCS5NyHkIwjZmQsnAJyglG6S3y+D5NAJ+ZkHeQAOU0pLKaVNAL4AkA4hP3PjaspL/Y68rNYVxsvGBFcJQsh4AP0BjKWmc7OE7No/HSENfv0m2zCBALYRQnwh5AegfTlwWwBEEkLCCCE2kDYo/uc65+lvj7xG+F0Aeyilz2s++g+A8fL/4wF8pUkfJUf8CYO0iXSzvPSkihDSXb7nLZrvCK4BlNKFlNJASmkopPa0ilI6DkJ2ZgGl9DSA44SQaDkpF8BuCPmZC8cAdCeEOMjlngtpD7GQn3lxNeWlvddwSDrZrGcB2jOEkAIA9wAYSCmt1XwkZNfOoZTuoJR2oJSGyjbMCQBd5H5RyA8AKKXt5gWgH6RIQQcB3He98yNeFAAyIE0z/w7gV/nVD9La4ZUADsh/PTTfuU+W4T4AhZr0VAA75c9ehXyQvHj9JXLMBvC1/L+QnZm8ACQB2Cq3v/8D4C7kZz4vAI8A2CuX/YcAbIX82u8LwCeQ9is2QTIYJ19NeQGwA/A5pKALmwGEX+9nvlFebciuGNK+J8V2eUPIrn2+jOSn+/wIAC8hP9NLeTCBQCAQCAQCgUAgELRz2tMSSoFAIBAIBAKBQCAQ/AHCgRMIBAKBQCAQCAQCM0E4cAKBQCAQCAQCgUBgJggHTiAQCAQCgUAgEAjMBOHACQQCgUAgEAgEAoGZIBw4gUAgEAgEAoFAIDAThAMnEAgEAoFAIBAIBGaCcOAEAoFAIBAIBAKBwEz4f30c4hWFH0k1AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plot = nk.events_plot(events, data)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "After you have located the events, you can now create epochs using the NeuroKit [epochs_create()](https://neurokit2.readthedocs.io/en/latest/functions.html#neurokit2.epochs_create>) function. However, we recommend to process your signal first before cutting them to smaller epochs. You can read more about processing of physiological signals using NeuroKit in [Custom your Processing Pipeline](https://neurokit2.readthedocs.io/en/latest/examples/custom.html>) Example." ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "# Process the signal\n", "df, info = nk.bio_process(ecg=data[\"ECG\"], rsp=data[\"RSP\"], eda=data[\"EDA\"], sampling_rate=100)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now, let's think about how we want our epochs to be like. For this example, we want:\n", " \n", " 1. Epochs to start *1 second before the event onset*\n", "\n", " 2. Epochs to end *6 seconds* afterwards\n", "\n", "These are passed into the `epochs_start` and `epochs_end` arguments, respectively. \n", "\n", "Our epochs will then cover the region from **-1 s** to **+6 s** relative to the onsets of events (i.e., 700 data points since the signal is sampled at 100Hz)." ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "# Build and plot epochs\n", "epochs = nk.epochs_create(df, events, sampling_rate=100, epochs_start=-1, epochs_end=6)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And as easy as that, you have created a dictionary of four dataframes, each correspond to an epoch of the event. \n", "\n", "Here, in the above example, all your epochs have the same starting time and ending time, specified by `epochs_start` and `epochs_end`. Nevertheless, you can also pass a list of different timings to these two arguments to customize the duration of the epochs for each of your events. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## One subject with multiple data files" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In some experimental designs, instead of having one signal file with multiple events, each subject can have multiples files where each file is the record of one event. \n", "\n", "In the following example, we will show you how to create a loop through the subject folders and put the files together in an epoch format for further analysis. " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Firstly, let's say your data is arranged as the following where each subject has a folder and in each folder there are multiple data files corresponding to different events:" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "```\n", "[Experiment folder]\n", "|\n", "└── Data\n", "| |\n", "| └── Subject_001/\n", "| | │ event_1.[csv]\n", "| | │ event_2.[csv]\n", "| | |__ ......\n", "| └── Subject_002/\n", "| │ event_1.[csv]\n", "| │ event_2.[csv]\n", "| |__ ......\n", "└── analysis_script.py\n", "```" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following will illustrate how your analysis script might look like. Try to re-create such data structure and the analysis script in your computer! " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now, in our analysis scripts, let's load the necessary packages:" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [], "source": [ "# Load packages\n", "import pandas as pd\n", "import os" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Assuming that your working directory is now at your analysis script, and you want to read all the data files of `Subject_001`. \n", "\n", "Your analysis script should look something like below:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "\n", "# Your working directory should be at Experiment folder\n", "participant = 'Subject_001'\n", "\n", "sampling_rate=100\n", "\n", "# List all data files in Subject_001 folder\n", "all_files = os.listdir('Data/' + participant)\n", "\n", "# Create an empty directory to store your files (events)\n", "epochs = {}\n", "\n", "# Loop through each file in the subject folder\n", "for i, file in enumerate(all_files):\n", " # Read the file\n", " data = pd.read_csv('Data/' + participant + '/' + file)\n", " # Add a Label column (e.g Label 1 for epoch 1)\n", " data['Label'] = np.full(len(data), str(i+1))\n", " # Set index of data to time in seconds\n", " index = data.index/sampling_rate\n", " data = data.set_index(pd.Series(index))\n", " # Append the file into the dictionary\n", " epochs[str(i + 1)] = data\n", " \n", "epochs" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And tah dah! You now should have a `dictionary` called epochs that resembles the output of NeuroKit [epochs_create()](https://neurokit2.readthedocs.io/en/latest/functions.html#neurokit2.epochs_create>). Each `DataFrame` in the epochs corresponds to an event (a trial) that `Subject_001` underwent. \n", "\n", "The epochs is now ready for further analysis!" ] } ], "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.7.1" } }, "nbformat": 4, "nbformat_minor": 2 }