Lugovtsov/cosmic-rays/Obrabotka.ipynb

484 lines
93 KiB
Plaintext

{
"cells": [
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"id": "kaayLGB4AyPP"
},
"outputs": [],
"source": [
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"df = pd.read_csv('/content/Untitled_angle_09_12.dat',delimiter = '\\t')\n",
"df.columns=['Count','Time']\n",
"ns = 0.064\n",
"df['Time'] = df.index * ns\n",
"df['Speed'] = 0.46/df['Time']*10**9\n",
"correct_df = df[(df['Time']>1.6) & (df['Time']<15)]"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 206
},
"id": "QTTbLpFSBWIG",
"outputId": "09b8c4f2-764e-40e7-f1cc-fce0589a6d15"
},
"outputs": [
{
"data": {
"text/html": [
"\n",
" <div id=\"df-8adb9e0c-ac81-43d8-84bb-060e6a6df807\" class=\"colab-df-container\">\n",
" <div>\n",
"<style scoped>\n",
" .dataframe tbody tr th:only-of-type {\n",
" vertical-align: middle;\n",
" }\n",
"\n",
" .dataframe tbody tr th {\n",
" vertical-align: top;\n",
" }\n",
"\n",
" .dataframe thead th {\n",
" text-align: right;\n",
" }\n",
"</style>\n",
"<table border=\"1\" class=\"dataframe\">\n",
" <thead>\n",
" <tr style=\"text-align: right;\">\n",
" <th></th>\n",
" <th>Count</th>\n",
" <th>Time</th>\n",
" <th>Speed</th>\n",
" </tr>\n",
" </thead>\n",
" <tbody>\n",
" <tr>\n",
" <th>26</th>\n",
" <td>31</td>\n",
" <td>1.664</td>\n",
" <td>2.764423e+08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>27</th>\n",
" <td>31</td>\n",
" <td>1.728</td>\n",
" <td>2.662037e+08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>28</th>\n",
" <td>23</td>\n",
" <td>1.792</td>\n",
" <td>2.566964e+08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>29</th>\n",
" <td>26</td>\n",
" <td>1.856</td>\n",
" <td>2.478448e+08</td>\n",
" </tr>\n",
" <tr>\n",
" <th>30</th>\n",
" <td>22</td>\n",
" <td>1.920</td>\n",
" <td>2.395833e+08</td>\n",
" </tr>\n",
" </tbody>\n",
"</table>\n",
"</div>\n",
" <div class=\"colab-df-buttons\">\n",
"\n",
" <div class=\"colab-df-container\">\n",
" <button class=\"colab-df-convert\" onclick=\"convertToInteractive('df-8adb9e0c-ac81-43d8-84bb-060e6a6df807')\"\n",
" title=\"Convert this dataframe to an interactive table.\"\n",
" style=\"display:none;\">\n",
"\n",
" <svg xmlns=\"http://www.w3.org/2000/svg\" height=\"24px\" viewBox=\"0 -960 960 960\">\n",
" <path d=\"M120-120v-720h720v720H120Zm60-500h600v-160H180v160Zm220 220h160v-160H400v160Zm0 220h160v-160H400v160ZM180-400h160v-160H180v160Zm440 0h160v-160H620v160ZM180-180h160v-160H180v160Zm440 0h160v-160H620v160Z\"/>\n",
" </svg>\n",
" </button>\n",
"\n",
" <style>\n",
" .colab-df-container {\n",
" display:flex;\n",
" gap: 12px;\n",
" }\n",
"\n",
" .colab-df-convert {\n",
" background-color: #E8F0FE;\n",
" border: none;\n",
" border-radius: 50%;\n",
" cursor: pointer;\n",
" display: none;\n",
" fill: #1967D2;\n",
" height: 32px;\n",
" padding: 0 0 0 0;\n",
" width: 32px;\n",
" }\n",
"\n",
" .colab-df-convert:hover {\n",
" background-color: #E2EBFA;\n",
" box-shadow: 0px 1px 2px rgba(60, 64, 67, 0.3), 0px 1px 3px 1px rgba(60, 64, 67, 0.15);\n",
" fill: #174EA6;\n",
" }\n",
"\n",
" .colab-df-buttons div {\n",
" margin-bottom: 4px;\n",
" }\n",
"\n",
" [theme=dark] .colab-df-convert {\n",
" background-color: #3B4455;\n",
" fill: #D2E3FC;\n",
" }\n",
"\n",
" [theme=dark] .colab-df-convert:hover {\n",
" background-color: #434B5C;\n",
" box-shadow: 0px 1px 3px 1px rgba(0, 0, 0, 0.15);\n",
" filter: drop-shadow(0px 1px 2px rgba(0, 0, 0, 0.3));\n",
" fill: #FFFFFF;\n",
" }\n",
" </style>\n",
"\n",
" <script>\n",
" const buttonEl =\n",
" document.querySelector('#df-8adb9e0c-ac81-43d8-84bb-060e6a6df807 button.colab-df-convert');\n",
" buttonEl.style.display =\n",
" google.colab.kernel.accessAllowed ? 'block' : 'none';\n",
"\n",
" async function convertToInteractive(key) {\n",
" const element = document.querySelector('#df-8adb9e0c-ac81-43d8-84bb-060e6a6df807');\n",
" const dataTable =\n",
" await google.colab.kernel.invokeFunction('convertToInteractive',\n",
" [key], {});\n",
" if (!dataTable) return;\n",
"\n",
" const docLinkHtml = 'Like what you see? Visit the ' +\n",
" '<a target=\"_blank\" href=https://colab.research.google.com/notebooks/data_table.ipynb>data table notebook</a>'\n",
" + ' to learn more about interactive tables.';\n",
" element.innerHTML = '';\n",
" dataTable['output_type'] = 'display_data';\n",
" await google.colab.output.renderOutput(dataTable, element);\n",
" const docLink = document.createElement('div');\n",
" docLink.innerHTML = docLinkHtml;\n",
" element.appendChild(docLink);\n",
" }\n",
" </script>\n",
" </div>\n",
"\n",
"\n",
"<div id=\"df-0afaef0d-7249-44f0-8914-0d49d8966d79\">\n",
" <button class=\"colab-df-quickchart\" onclick=\"quickchart('df-0afaef0d-7249-44f0-8914-0d49d8966d79')\"\n",
" title=\"Suggest charts\"\n",
" style=\"display:none;\">\n",
"\n",
"<svg xmlns=\"http://www.w3.org/2000/svg\" height=\"24px\"viewBox=\"0 0 24 24\"\n",
" width=\"24px\">\n",
" <g>\n",
" <path d=\"M19 3H5c-1.1 0-2 .9-2 2v14c0 1.1.9 2 2 2h14c1.1 0 2-.9 2-2V5c0-1.1-.9-2-2-2zM9 17H7v-7h2v7zm4 0h-2V7h2v10zm4 0h-2v-4h2v4z\"/>\n",
" </g>\n",
"</svg>\n",
" </button>\n",
"\n",
"<style>\n",
" .colab-df-quickchart {\n",
" --bg-color: #E8F0FE;\n",
" --fill-color: #1967D2;\n",
" --hover-bg-color: #E2EBFA;\n",
" --hover-fill-color: #174EA6;\n",
" --disabled-fill-color: #AAA;\n",
" --disabled-bg-color: #DDD;\n",
" }\n",
"\n",
" [theme=dark] .colab-df-quickchart {\n",
" --bg-color: #3B4455;\n",
" --fill-color: #D2E3FC;\n",
" --hover-bg-color: #434B5C;\n",
" --hover-fill-color: #FFFFFF;\n",
" --disabled-bg-color: #3B4455;\n",
" --disabled-fill-color: #666;\n",
" }\n",
"\n",
" .colab-df-quickchart {\n",
" background-color: var(--bg-color);\n",
" border: none;\n",
" border-radius: 50%;\n",
" cursor: pointer;\n",
" display: none;\n",
" fill: var(--fill-color);\n",
" height: 32px;\n",
" padding: 0;\n",
" width: 32px;\n",
" }\n",
"\n",
" .colab-df-quickchart:hover {\n",
" background-color: var(--hover-bg-color);\n",
" box-shadow: 0 1px 2px rgba(60, 64, 67, 0.3), 0 1px 3px 1px rgba(60, 64, 67, 0.15);\n",
" fill: var(--button-hover-fill-color);\n",
" }\n",
"\n",
" .colab-df-quickchart-complete:disabled,\n",
" .colab-df-quickchart-complete:disabled:hover {\n",
" background-color: var(--disabled-bg-color);\n",
" fill: var(--disabled-fill-color);\n",
" box-shadow: none;\n",
" }\n",
"\n",
" .colab-df-spinner {\n",
" border: 2px solid var(--fill-color);\n",
" border-color: transparent;\n",
" border-bottom-color: var(--fill-color);\n",
" animation:\n",
" spin 1s steps(1) infinite;\n",
" }\n",
"\n",
" @keyframes spin {\n",
" 0% {\n",
" border-color: transparent;\n",
" border-bottom-color: var(--fill-color);\n",
" border-left-color: var(--fill-color);\n",
" }\n",
" 20% {\n",
" border-color: transparent;\n",
" border-left-color: var(--fill-color);\n",
" border-top-color: var(--fill-color);\n",
" }\n",
" 30% {\n",
" border-color: transparent;\n",
" border-left-color: var(--fill-color);\n",
" border-top-color: var(--fill-color);\n",
" border-right-color: var(--fill-color);\n",
" }\n",
" 40% {\n",
" border-color: transparent;\n",
" border-right-color: var(--fill-color);\n",
" border-top-color: var(--fill-color);\n",
" }\n",
" 60% {\n",
" border-color: transparent;\n",
" border-right-color: var(--fill-color);\n",
" }\n",
" 80% {\n",
" border-color: transparent;\n",
" border-right-color: var(--fill-color);\n",
" border-bottom-color: var(--fill-color);\n",
" }\n",
" 90% {\n",
" border-color: transparent;\n",
" border-bottom-color: var(--fill-color);\n",
" }\n",
" }\n",
"</style>\n",
"\n",
" <script>\n",
" async function quickchart(key) {\n",
" const quickchartButtonEl =\n",
" document.querySelector('#' + key + ' button');\n",
" quickchartButtonEl.disabled = true; // To prevent multiple clicks.\n",
" quickchartButtonEl.classList.add('colab-df-spinner');\n",
" try {\n",
" const charts = await google.colab.kernel.invokeFunction(\n",
" 'suggestCharts', [key], {});\n",
" } catch (error) {\n",
" console.error('Error during call to suggestCharts:', error);\n",
" }\n",
" quickchartButtonEl.classList.remove('colab-df-spinner');\n",
" quickchartButtonEl.classList.add('colab-df-quickchart-complete');\n",
" }\n",
" (() => {\n",
" let quickchartButtonEl =\n",
" document.querySelector('#df-0afaef0d-7249-44f0-8914-0d49d8966d79 button');\n",
" quickchartButtonEl.style.display =\n",
" google.colab.kernel.accessAllowed ? 'block' : 'none';\n",
" })();\n",
" </script>\n",
"</div>\n",
" </div>\n",
" </div>\n"
],
"text/plain": [
" Count Time Speed\n",
"26 31 1.664 2.764423e+08\n",
"27 31 1.728 2.662037e+08\n",
"28 23 1.792 2.566964e+08\n",
"29 26 1.856 2.478448e+08\n",
"30 22 1.920 2.395833e+08"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"correct_df.head()"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 466
},
"id": "5yKqwV9SxLHP",
"outputId": "735be114-40dd-449f-8f4a-b48be1d2a5b5"
},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 0, 'Time, ns')"
]
},
"execution_count": 12,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjIAAAGwCAYAAACzXI8XAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABFmUlEQVR4nO3de3xU9Z3/8fcAIQmXDHLLhBouIl7CRYVWi1hX5V4Xa611tYJIu27LA9CodS1bEVltKbq9aNWw2q3Yn6W03ZUKVuBBEfDSIApGpSgVCoKaQAWSQDCBJvP7AyZmkrmd+zmT1/PxmIcyc873+znfOWE+zMx5JxSNRqMCAAAIoA5eFwAAAGAWjQwAAAgsGhkAABBYNDIAACCwaGQAAEBg0cgAAIDAopEBAACB1cnrApzW1NSkjz/+WN27d1coFPK6HAAAkIFoNKojR46oX79+6tAh+fsuWd/IfPzxxyouLva6DAAAYMK+fft0+umnJ3086xuZ7t27Szq5EAUFBR5XAwAAMlFbW6vi4uLm1/Fksr6RiX2cVFBQQCMDAEDApPtaCF/2BQAAgUUjAwAAAotGBgAABBaNDAAACCwaGQAAEFg0MgAAILBoZAAAQGDRyAAAgMCikQEAAIGV9cm+QdbYFNXm3Yd04Ei9+nbP04WDeqpjB37xJQAAMTQyPrV6W6UWrNyuypr65vuKwnmaP6VEk4YVeVgZAAD+wUdLPrR6W6VmPrM1romRpKqaes18ZqtWb6v0qDIAAPyFRsZnGpuiWrByu6IJHovdt2DldjU2JdoCAID2hUbGZzbvPtTmnZiWopIqa+q1efch94oCAMCnaGR85sCR5E2Mme0AAMhmNDI+07d7nq3bAQCQzWhkfObCQT1VFM5TsousQzp59dKFg3q6WRYAAL5EI+MzHTuENH9KiSS1aWZif54/pYQ8GQAARCPjS5OGFals6khFwvEfH0XCeSqbOpIcGQAATiEQz6cmDSvS+JIIyb4AAKTg6TsyZWVlGjFihAoKClRQUKDRo0dr1apVzY9fdtllCoVCcbfvfOc7Hlbsro4dQho9uJe+cv7nNHpwL5oYAABa8fQdmdNPP10/+tGPNGTIEEWjUT399NP6yle+ojfffFNDhw6VJN1yyy36z//8z+Z9unTp4lW5AADAZzxtZKZMmRL35x/84AcqKyvTpk2bmhuZLl26KBKJZDxmQ0ODGhoamv9cW1trT7EAAMB3fPNl38bGRi1btkx1dXUaPXp08/2//vWv1bt3bw0bNkxz587VsWPHUo6zcOFChcPh5ltxcbHTpQMAAI+EotGop7+055133tHo0aNVX1+vbt26aenSpfryl78sSXriiSc0YMAA9evXT2+//bbuvvtuXXjhhXr22WeTjpfoHZni4mLV1NSooKDA8eMBAADW1dbWKhwOp3399ryROX78uPbu3auamhr97//+r37xi19o48aNKikpabPtiy++qLFjx2rnzp0aPHhwRuNnuhAAAMA/Mn399vyjpc6dO+vMM8/UqFGjtHDhQp133nl6+OGHE2570UUXSZJ27tzpZokAAMCnPG9kWmtqaor7aKiliooKSVJREYFwAADA46uW5s6dq8mTJ6t///46cuSIli5dqg0bNmjNmjXatWtX8/dlevXqpbffflu33367Lr30Uo0YMcLLsgEAgE942sgcOHBAN910kyorKxUOhzVixAitWbNG48eP1759+/SnP/1JP/vZz1RXV6fi4mJ97Wtf0z333ONlyQAAwEc8/7Kv0/iyLwAAwROYL/sCAACYRSMDAAACi0YGAAAEFo0MAAAILBoZAAAQWDQyAAAgsGhkAABAYNHIAACAwKKRAQAAgUUjAwAAAotGBgAABBaNDAAACCwaGQAAEFg0MgAAILBoZAAAQGDRyAAAgMCikQEAAIFFIwMAAAKLRgYAAAQWjQwAAAgsGhkAABBYNDIAACCwaGQAAEBg0cgAAIDAopEBAACB1cnrAoKosSmqzbsP6cCRevXumiuFpE+ONqhv9zxdOKinOnYIeV0iAADtAo2MQau3VWrByu2qrKlP+HhROE/zp5Ro0rAilysDAKD94aMlA1Zvq9TMZ7YmbWIkqaqmXjOf2arV2ypdrAwAgPaJRiZDjU1RLVi5XdE028UeX7Byuxqb0m0NAACsoJHJ0Obdh1K+E9NSVFJlTb027z7kbFEAALRzNDIZOnAksybG6j4AACBzNDIZ6ts9z5V9AABA5mhkMnThoJ4qCucpkwurQzp59dKFg3o6XRYAAO0ajUyGOnYIaf6UEklK2czEHps/pYQ8GQAAHEYjY8CkYUUqmzpSkXDyj4wi4TyVTR3pqxyZxqaoyncd1HMVH6l810GupgIAZA1PA/HKyspUVlamPXv2SJKGDh2qe++9V5MnT5Yk1dfX684779SyZcvU0NCgiRMn6vHHH1dhYaFnNU8aVqTxJZHAJPsmCvAjtA8AkC1C0WjUs3+er1y5Uh07dtSQIUMUjUb19NNP66GHHtKbb76poUOHaubMmfrjH/+oJUuWKBwOa/bs2erQoYNeffXVjOeora1VOBxWTU2NCgoKHDwa/4kF+LV+gmNtlt/eOQIAICbT129PG5lEevbsqYceekjXXnut+vTpo6VLl+raa6+VJL333ns699xzVV5eri9+8YsZjddeG5nGpqguWfRi0uybkE5+DPbK3Vf46h0kAACkzF+/ffMdmcbGRi1btkx1dXUaPXq0tmzZohMnTmjcuHHN25xzzjnq37+/ysvLk47T0NCg2trauFt7lC7Aj9A+AEA28LyReeedd9StWzfl5ubqO9/5jpYvX66SkhJVVVWpc+fO6tGjR9z2hYWFqqqqSjrewoULFQ6Hm2/FxcUOH4E/ZRrGR2gfACDIPG9kzj77bFVUVOi1117TzJkzNX36dG3fvt30eHPnzlVNTU3zbd++fTZWGxyZhvER2gcACDJPr1qSpM6dO+vMM8+UJI0aNUqvv/66Hn74Yf3Lv/yLjh8/rurq6rh3Zfbv369IJJJ0vNzcXOXm5jpdtu/FAvyqauoT/qLL2HdkCO0DAASZ5+/ItNbU1KSGhgaNGjVKOTk5WrduXfNjO3bs0N69ezV69GgPKwyGVAF+hPYBALKFp+/IzJ07V5MnT1b//v115MgRLV26VBs2bNCaNWsUDof1rW99S3fccYd69uypgoICzZkzR6NHj874iiW/amyKNufQOJk9Ewvwa50jEyFHBgCQJTxtZA4cOKCbbrpJlZWVCofDGjFihNasWaPx48dLkn7605+qQ4cO+trXvhYXiBdkbgfUtQ7w82NoHwAAZvkuR8ZufsqRIaAOAIDMBC5HJts1NkW1YOX2hF+8jd23YOV2fg8SAAAG0Mi4hIA6AADsRyPjEgLqAACwH42MSwioAwDAfjQyLokF1CW7Viikk1cvEVAHAEDmaGRcQkAdAAD2o5FxUSygLhKO//goEs7j0msAAEzw/HctZbvWKb7jSyKWAurcSgUGACAIaGQcZHeKr9upwAAA+B0fLTkkluLbOjumqqZeM5/ZqtXbKj0dDwCAbEAj4wC7U3xJBQYAIDEaGQfYneJLKjAAAInRyDjA7hRfUoEBAEiMRsYBdqf4kgoMAEBiNDIOsDvFl1RgAAASo5FxgN0pvqQCAwCQWCgajWb1pS61tbUKh8OqqalRQUGBq3ObzX1JFnqXarxYyF5Vzac6VHdcPbvlKlKQOjAvNo+RfQAAcEOmr980Mg4zmsSbrvlJNN7a7VVt9km0b7p50u0DAIBbaGRO8bqRMSIWetf6CYm1PYl+H1OyfVrv33JfM/sAAOCmTF+/+Y6MT5gJvUu1T2uxfc3sAwCAX9HI+ISZ0Lt0+yTa18w+AAD4Fb800ifMhN4ZDcAzE5hHyB4AwM9oZHzCTOid0QA8M4F5hOwBAPyMj5Z8wkzoXWyfdFrua2YfAAD8ikbGJ8yE3sX2ySTxJbavmX0AAPArGhmbNTZFVb7roJ6r+Ejluw42X/WT7P6WxpdEVDruLIXzc+Luj4Tz9Ng3LlA4v3Ob/ScNK1LZ1JFJ32UpCue1uYzayD6Z1A0AgFfIkbFRsjC7q84r0oq3KlMm/Cbat0d+jmaMGaghfbvr/j+mTgg2k9Kbbh+zycQAAFhFIN4pbjUymYTMtdQy5E5S0iC8ZOOlCsmzg5lwPgAA7EIgnouMhMzFxLa9b8VfdN+K1EF4qfZ3IrTOTDgfAABeoJGxQaYhc61FJVXVNqiq1lxWi1OhdWbC+QAA8AKNjA28Do2ze34z4XwAAHiBRsYGXofG2T2/mXA+AAC8QCNjg3RhdsmEJEUKTl4pZCatxanQOjPhfAAAeIFGxgapwuySiW1331VDdd9VqYPwUj3mRGidmXA+AAC8QCNjk1jIXKRVyFxROE/fvnRQm/C5SIvQuWT7RsJ5Wjx1pBYneczJS6BT1cSl1wAAvyBHxmaxkLkDR+rVt/tn4XLJ7s9k33SPuXE8RoL2/MCr9QIA2CPT129Pf/v1woUL9eyzz+q9995Tfn6+Lr74Yi1atEhnn3128zaXXXaZNm7cGLfft7/9bS1evNjtcjPSsUNIowf3yvj+TLfJZH8ndOwQUs2nx/Xgmh2BSfglkRgA2g9PP1rauHGjZs2apU2bNmnt2rU6ceKEJkyYoLq6urjtbrnlFlVWVjbfHnzwQY8qbn9iCb+tc2Wqauo185mtWr2t0qPKEgtavQAAazx9R2b16tVxf16yZIn69u2rLVu26NJLL22+v0uXLopEIm6X1+6lS/gN6WTC7/iSiC8+tglavQAA63z1Zd+amhpJUs+e8Zf1/vrXv1bv3r01bNgwzZ07V8eOHUs6RkNDg2pra+NuMCdoCb9BqxcAYJ2n78i01NTUpNLSUo0ZM0bDhg1rvv8b3/iGBgwYoH79+untt9/W3XffrR07dujZZ59NOM7ChQu1YMECt8rOakFL+A1avQAA63zTyMyaNUvbtm3TK6+8Enf/v/3bvzX///Dhw1VUVKSxY8dq165dGjx4cJtx5s6dqzvuuKP5z7W1tSouLnau8CwWtITfoNULALDOFx8tzZ49W88//7zWr1+v008/PeW2F110kSRp586dCR/Pzc1VQUFB3A3mBC3hN2j1AgCs87SRiUajmj17tpYvX64XX3xRgwYNSrtPRUWFJKmoiMtonRa0hN+g1QsAsM7TRmbWrFl65plntHTpUnXv3l1VVVWqqqrSp59+KknatWuX7r//fm3ZskV79uzRihUrdNNNN+nSSy/ViBEjvCy93Qhawm/Q6gUAWONpsm8olPhfxk899ZRuvvlm7du3T1OnTtW2bdtUV1en4uJiffWrX9U999yT8UdGbif7ZqugJeUGrV4AQLxMX7/5FQUAAMB3Mn399sWXfQEAAMygkQEAAIFFIwMAAAKLRgYAAAQWjQwAAAgsGhkAABBYvvldS0FHbglrEBQ8TwCyCY2MDVZvq9SCldtVWfPZb1UuCudp/pSSdpMkyxoEA88TgGzDR0sWrd5WqZnPbI17YZCkqpp6zXxmq1Zvq/SoMvewBsHA8wQgG9HIWNDYFNWClduVKBo5dt+CldvV2JS94cmsQTDwPAHIVjQyFmzefajNv25bikqqrKnX5t2H3CvKZaxBMPA8AchWNDIWHDiS/IXBzHZBxBoEA88TgGxFI2NB3+55tm4XRKxBMPA8AchWNDIWXDiop4rCeUp24WpIJ68IuXBQTzfLchVrEAw8TwCyFY2MBR07hDR/SokktXmBiP15/pSSrM7oYA2CgecJQLaikbFo0rAilU0dqUg4/i35SDhPZVNHtotsDtYgGHieAGSjUDQazerrLWtraxUOh1VTU6OCggLH5mmdljpqwGna8sFhT9NTrSS4mtmXxNhg4HkCEASZvn6T7GuTjh1CGj24l6STwWP/9NB6T9NTrSS4mt235RrAv3ieAGQTPlqymR/SU63U4If6AQDIFI2MjfyQnmqlBj/UDwCAETQyNvJDeqqVGvxQPwAARtDI2MgP6alWavBD/QAAGEEjYyM/pKdaqcEP9QMAYASNjI38kJ5qpQY/1A8AgBE0MjbyQ3qqlRr8UD8AAEYQiOeATHNY7AgmSzaGUzky40si2rz7kKpqPtWhuuPq2S1XfbvlSiHpk6MNBKwBAGyR6es3jYxD0jUpVhqNTMewO9l37faqNvMl4nb4HwAg+9DInOJVI5NKLHSu9cLHWoxMfu+NHWMYkWy+RJyqAQDQfmT6+s13ZFxmR+ic28F1qeZLhPA8AIBbaGRcZkfonNvBdenmc6MGAAASoZFxmR2hc24H11kZh/A8AICTaGRcZkfonNvBdVbGITwPAOAkGhmX2RE653ZwXWw+IwjPAwC4gUbGZXaEzrkdXBebL9PRCM8DALiFRsYDk4YVqWzqSEVavcsRCedlfMmylTEam6Iq33VQz1V8pPJdBzO6sig2XybvzBg5DgAArCBHxkNOJvsmYzWILzYfyb4AACcFIhBv4cKFevbZZ/Xee+8pPz9fF198sRYtWqSzzz67eZv6+nrdeeedWrZsmRoaGjRx4kQ9/vjjKiwszGgOPzcybnM7RA8AALMCEYi3ceNGzZo1S5s2bdLatWt14sQJTZgwQXV1dc3b3H777Vq5cqV+//vfa+PGjfr44491zTXXeFh1MLkdogcAgBt89dHS3//+d/Xt21cbN27UpZdeqpqaGvXp00dLly7VtddeK0l67733dO6556q8vFxf/OIX24zR0NCghoaG5j/X1taquLi43b8jU77roG54clPa7X5zyxc1enAvFyoCACC5QLwj01pNTY0kqWfPk5fsbtmyRSdOnNC4ceOatznnnHPUv39/lZeXJxxj4cKFCofDzbfi4mLnCw8At0P0AABwg6lG5owzztDBgwfb3F9dXa0zzjjDVCFNTU0qLS3VmDFjNGzYMElSVVWVOnfurB49esRtW1hYqKqqqoTjzJ07VzU1Nc23ffv2maon27gdogcAgBs6mdlpz549amxsbHN/Q0ODPvroI1OFzJo1S9u2bdMrr7xiav+Y3Nxc5ebmWhojG8VC7apq6hN+Tyakk5dNE2AHAAgSQ43MihUrmv9/zZo1CofDzX9ubGzUunXrNHDgQMNFzJ49W88//7xeeuklnX766c33RyIRHT9+XNXV1XHvyuzfv1+RSMTwPO1ZLNRu5jNbFZLimhkC7AAAQWWokbn66qslSaFQSNOnT497LCcnRwMHDtSPf/zjjMeLRqOaM2eOli9frg0bNmjQoEFxj48aNUo5OTlat26dvva1r0mSduzYob1792r06NFGSoc+C7VrnSMTMZAjAwCAnxhqZJqamiRJgwYN0uuvv67evXtbmnzWrFlaunSpnnvuOXXv3r35ey/hcFj5+fkKh8P61re+pTvuuEM9e/ZUQUGB5syZo9GjRye8Yimo7AjGy9SkYUUaXxJxbb6WWh5n766E6AEArPP08utQKPEL11NPPaWbb75Z0meBeL/5zW/iAvEy/WjJ74F4VpN2gyLRcbaUjccMADDP8WTfdevWad26dTpw4EDzOzUxv/zlL80M6Qg/NzLtJWk32XG2lG3HDACwxtEcmQULFmjChAlat26dPvnkEx0+fDjuhvTaS9JuquNsKZuOGQDgHlOXXy9evFhLlizRtGnT7K6n3di8+1DSj1mkky/slTX12rz7UKCTdtMdZ0vZcswAAPeYekfm+PHjuvjii+2upV1pL0m7ZuoP+jEDANxjqpH513/9Vy1dutTuWtqV9pK0a6b+oB8zAMA9pj5aqq+v1xNPPKE//elPGjFihHJycuIe/8lPfmJLcdmsvSTtpjvOlrLlmAEA7jHVyLz99ts6//zzJUnbtm2LeyzZJdWI116SdlMdZ0vZdMwAAPd4miPjBj9ffi2Zz5FxKkTPqXHT5cj0yM/RjDEDNfuKIZ41Mm4GEwIAUnM8RyYo/N7ISMZfQJ0K0XM6nK91su/rew5pyZ/3qPrTE47MZ0R7CSYEgKBwtJG5/PLLU36E9OKLLxod0jFBaGSMcCpEz+1wPj+FAfqpFgDASY4G4p1//vk677zzmm8lJSU6fvy4tm7dquHDh5suGqk5FaLndjifn8IA/VQLAMA4U1/2/elPf5rw/vvuu09Hjx61VBCScypEz+1wPj+FAfqpFgCAcabekUlm6tSpvvo9S9nGqRA9t8P5/BQG6KdaAADG2drIlJeXKy+PMDOnOBWi53Y4n5/CAP1UCwDAOFMfLV1zzTVxf45Go6qsrNQbb7yhefPm2VIY2nIqRM/tcD4/hQH6qRYAgHGm3pEJh8Nxt549e+qyyy7TCy+8oPnz59tdI06JhctJn11RE2MlUM6pcf0yX1BqAQAYR45MAAU1RyaT+VoH41kNqct0/0yOPdFYklwJ0SOsD0B740og3pYtW/Tuu+9KkoYOHaoLLrjA7FCOycZGRgpesm+q+R59caeeenV3m2C8q84r0oq3Kk03VkYbs1THnrDp6nLyd4xVH3M20I+wPgDtkaONzIEDB3T99ddrw4YN6tGjhySpurpal19+uZYtW6Y+ffqYLtxu2drIZItkYXTJZBpSZ2fInZEa7Q7RI6wPQHvlaCDenDlzdOTIEf3lL3/RoUOHdOjQIW3btk21tbW69dZbTReN9iVVGF0ymYTU2RlyZ7RGO0P0COsDgPRMNTKrV6/W448/rnPPPbf5vpKSEj322GNatWqVbcUhu6ULo0umZUidmXHT7W+1RiPjW5nbrnkAIMhMNTJNTU3Kyclpc39OTo6amposF4X2wWrIXLL97Qy5s1KjU8dn9zwAEGSmGpkrrrhCt912mz7++OPm+z766CPdfvvtGjt2rG3FIbtZDZlLtr+dIXdWanTq+OyeBwCCzFQj8+ijj6q2tlYDBw7U4MGDNXjwYA0aNEi1tbX6+c9/bneNyFKxMDqj10SFdPKqnWQhdenGTbe/1RqNjG9lbrvmAYAgM9XIFBcXa+vWrfrjH/+o0tJSlZaW6oUXXtDWrVt1+umn210jslSqMLpkMgmpszPkzmiNdoboEdYHAOkZamRefPFFlZSUqLa2VqFQSOPHj9ecOXM0Z84cfeELX9DQoUP18ssvO1UrstCkYUUqmzpSkXD8xyNF4Tx9+9JBKmp1fyScl9Elx8nGzXT/TMbq0SWnOUvGyvhm5rZ7HgAIKkM5MldddZUuv/xy3X777Qkff+SRR7R+/XotX77ctgKtIkfGfWZC9ZLtYzVN186Av0xqGTXgNG354LCqaj7Vobrj6tktV5EC68GCdhwH6cAAgsSRQLwBAwZo9erVcZddt/Tee+9pwoQJ2rt3r/GKHUIj4y6nU2j9nHKbqLYYr2v087oBQCKOBOLt378/4WXXMZ06ddLf//53I0Mii8RSaFu/kFfV1GvmM1u1elulr8e3IlltMZUe1ujndQMAqww1Mp/73Oe0bdu2pI+//fbbKiriX3ftkdMptH5OuTWS/ut2jX5eNwCwg6FG5stf/rLmzZun+vq2/+r89NNPNX/+fP3zP/+zbcUhOJxOofVzym2m6b9e1OjndQMAO3QysvE999yjZ599VmeddZZmz56ts88+W9LJ78Y89thjamxs1Pe//31HCoW/OZ1C6+eUW6Nzulmjn9cNAOxgqJEpLCzUn//8Z82cOVNz585V7HvCoVBIEydO1GOPPabCwkJHCoW/OZ1C6+eUW6Nzulmjn9cNAOxgqJGRTl659MILL+jw4cPauXOnotGohgwZotNOO82J+hAQsRTaqpr6hN/HCOlk9onZFFqnx7ciVlu6j5e8qNHP6wYAdjCV7CtJp512mr7whS/owgsvpImB4ym0fk65jdWWycxu1+jndQMAO5huZIDWkqXQFhbkqnTcEDX8o0nluw6avkIm2fjh/ByVjhuiK84pVPmug3qu4qM28zQ2RZM+ZofxJRGVjjtLPfITxxMUGUjiNVNrqn1itYXznU0hBgAvGArEs9tLL72khx56SFu2bFFlZaWWL1+uq6++uvnxm2++WU8//XTcPhMnTtTq1asznoNAPPe1TJDd88kx/WbzXlXV2hfE1tgU1aMv7tRTr+5W9acnmu/vEJJavubH5pHkekhfOK+TxpcUasyQPoaSfc0E16XaR2p77D3yczRjzEDNvmII78QA8C1Hkn3ttmrVKr366qsaNWqUrrnmmoSNzP79+/XUU08135ebm2vooywaGe/Egthan2Cxl06z7wYkG7e1kJR0G6s1pKvFzPhmxkq1j9PHDgBOciTZ126TJ0/WAw88oK9+9atJt8nNzVUkEmm+8X2cYHAqiM1I+FyqbfwW0mdmrEz2SYQgPADZxPffkdmwYYP69u2rs88+WzNnztTBgwdTbt/Q0KDa2tq4G9znVBBbpuFzmfBTSJ+ZsaysBUF4ALKFrxuZSZMm6Ve/+pXWrVunRYsWaePGjZo8ebIaGxuT7rNw4UKFw+HmW3FxsYsVI8apIDYngtv8ENJnZiw71oIgPABBZzhHxk3XX3998/8PHz5cI0aM0ODBg7VhwwaNHTs24T5z587VHXfc0fzn2tpamhkPOBXE5kRwmx9C+syMZcdaEIQHIOh8/Y5Ma2eccYZ69+6tnTt3Jt0mNzdXBQUFcTe4LxbEluyamJBOXlljNIgt3bhGmK0h01qMjG9mLCtrYfXYAcAvAtXIfPjhhzp48CC/YTsAnApiSzVua6Ek/2+1hkxqMTq+mbEy2ceO2gDAzzxtZI4ePaqKigpVVFRIknbv3q2Kigrt3btXR48e1V133aVNmzZpz549Wrdunb7yla/ozDPP1MSJE70sO7CcDoVrLVmAndUgtmTjtn5NjoTztHjqSC12oIZ0tbQeP5O1N7NeqUICbx83RI9/w7ljBwA/8DRHZsOGDbr88svb3D99+nSVlZXp6quv1ptvvqnq6mr169dPEyZM0P3332/oF1OSI3OSmaA1u7QMyOvbPfNwOKPjjhpwmrZ8cDjhPE7VkKyWluMbXXsztSYLCSwK52nelefqtK65jh07ADghEIF4bqCRcS6YDum5tfY8xwCyTSAC8eA8p4LpkJ5ba89zDKA9o5HJck4F0yE9t9ae5xhAe0Yjk+WcCqZDem6tPc8xgPaMRibLORVMh/TcWnueYwDtGY1MlnMqmA7pubX2PMcA2jMamSznVDAd0nNr7XmOAbRnNDLtgFPBdEjPrbXnOQbQXpEj0444HQrn9XxmuBWUV1XzqQ7VHVfPbrmKFHz2MY/dc7c8nt5dc6WQ9MnRhoTjB+H5AdB+EYh3Co2MN7xMEs6U0zWmGl+SZ3NPGlYUiOcHQPtGI3MKjYz7gpAy63SNqcZP9gPnxtyS9G+XDtITL+329fMDACT7whNBSJl1usZMxk/Ejbmjkp58uW0TY9f8AOA2GhnYKggps07XmG78VNyYO1WP4ofnBwCMoJGBrYKQMut0jXYcm5dz2zkOADiNRga2CkLKrNM12nFsXs5t5zgA4DQaGdgqCCmzTteYbvxU3Ji7Q6htcJ5d8wOA22hkYKsgpMw6XWMm43s1d0jSLV8a5Nj8AOA2Lr+GI/yeU9LYFNWjL+7UU6/uVvWnJ5rvLwrnad6V5+q0rrnNQXGjBpymLR8cTvrnZEFyRnNkwnmdNL6kUGOG9GkOzWs9bqYhdm7myKSrKVUooJGGiQA/GMH5EnzkyJxCI+Mdv/5FkuhFvEd+jmaMGaghfbvr/j/GP9YhFH+lT+s/p2oAUq1BsmYq2bhGm49MGwwrz4+ZhimT2o3OA7TE+ZIdaGROoZFBS2aC6tIxGySXrJbWY5dNHSlJvgsZNBu813rbdLUHIWAR/sH5kj0IxANaMRtUl46ZILlUtbR234q/6L4V/goZtBK811qq2oMQsAj/4Hxpn2hk0G5YCapLx2iQXKa1RCVV1TaoqtZfIYNWg/di0tUehIBF+AfnS/tEI4N2w42QNy8DAd0MsbN7rmTjBSFgEf7B+dI+0cig3XAj5M3LQEA3Q+zsnivZeEEIWIR/cL60TzQyaDesBNWlYzRILlZLJuNGCk5eruynELtMg/fSSVd7EAIW4R+cL+0TjQzaDbNBdemYCZKL1ZLJ1vddNVT3XeWvkMFMg/cyqShV7UEIWIR/cL60T1x+jcCwK5fGaFCdlRwZM7UkG9eNbAyja2w1RyYWPpguLM+pY/dr1hGsIUcmO5AjcwqNTHaw+y+mdEF1LR8zm+xrtJZMkm+dfOE1u8Zmk30P1x1vEz6Yal67j50Xu+xGkxp8NDKn0MgEHwFXznN7jY2EATrx3HJOAf5HIB6yAgFXznN7jY2EATrx3HJOAdmFRga+RsCV89xeYyNhgE48t5xTQHahkYGvEXDlPLfX2Og4boXvOT0vAGfQyMDXCLhynttrbHQct8L3nJ4XgDNoZOBrBFw5z+01NhIG6MRzyzkFZBcaGfgaAVfOc3uNjYQBOvHcck4B2YVGBr43aViRyqaOVKTVv+Ij4Twuk7WJ22scmy/ZOzNFDj+3nFNA9iBHBr6SKhxOEgFXDksWBphJWF+6sVIF/MXG79Gls6qPffZfp+ZLFW4ocZ4BfpDp63cnF2tq46WXXtJDDz2kLVu2qLKyUsuXL9fVV1/d/Hg0GtX8+fP15JNPqrq6WmPGjFFZWZmGDBniXdFwjJG4fjijY4eQRg/uJenk8/FPD6039Xxkmprbej6zz38m86Xa5ivnf85Q3QD8w9OPlurq6nTeeefpscceS/j4gw8+qEceeUSLFy/Wa6+9pq5du2rixImqr+eyyGwTS1pNlu9RWVOvmc9s1eptlS5X1j5ZeT6S7VtlYh+75rNrGwD+45uPlkKhUNw7MtFoVP369dOdd96p7373u5KkmpoaFRYWasmSJbr++uszGpePlvyvsSmqSxa9mDYkLaST32F45e4reKvfQVaej3T7mtnH6nyFBbmSQqqqtbYN5x7grsD/ioLdu3erqqpK48aNa74vHA7roosuUnl5edL9GhoaVFtbG3eDv3md9Ip4Vp4PM6m5Ts9XVduQtEExsg3nHuBPvm1kqqqqJEmFhYVx9xcWFjY/lsjChQsVDoebb8XFxY7WCeu8TnpFPCvPh5nUXDfmswvnHuA/vm1kzJo7d65qamqab/v27fO6JKThddIr4ll5Psyk5roxn1049wD/8W0jE4lEJEn79++Pu3///v3NjyWSm5urgoKCuBv8zeukV8Sz8nyYSc11er5IwclLuK1uw7kH+JNvG5lBgwYpEolo3bp1zffV1tbqtdde0+jRoz2sDHbzOukV8aw8H2ZSc52e776rhuq+q6xvw7kH+JOnjczRo0dVUVGhiooKSSe/4FtRUaG9e/cqFAqptLRUDzzwgFasWKF33nlHN910k/r16xeXNYPs4FbSa2NTVOW7Duq5io9UvuugGpt8cdFeQqlqjT22fOuH+p+X/6blb5o/nkTzjC+JqHTcWeqRn5Nwn0TPR2ycT4836msjP6eCVvu2TM1tPef4kojp5z9ZSm84P0e3jj1T3fNy1PCPJpWOO0uFBW23KR03RONLIpbSfpM9V0E634Cg8vTy6w0bNujyyy9vc//06dO1ZMmS5kC8J554QtXV1brkkkv0+OOP66yzzsp4Di6/DpZUyb5W/zUcpLCzVLVKsi04MNE8PbqcbECqj51ovi+c10njSwo1ZkifhM9HqjC7Hvk5mjFmoGZfMUQdO4RSHtv4kojp57+xKapHX9ypp17drepPTyTcJlKQq88P7KlX3v8kbpuW65ZJQnC6NSwK5+mq84q04q3KQJxvgB9l+vrtmxwZp9DIQPosNK31yR57efLT79dJVWsmP6whZXY8yeZJNqaSjJvJOLGaJDn2PBg5nkT1mZnf6Jx+PN8Avwp8jgxgl8amqBas3J7wxSZ234KV233xtn8mtWYi3fGkmieRZOtkZJz7VvxF961w5nkwejx2zG9mTr+db0A2oJFB1jMT0uaVTMPhUsnkeMzMYzXMzsnQObfWzY45/XS+AdmARgZZz0xIm1fsrCHVWFbmcTqQzsyYbq2bnXP64XwDsgGNDLKemZA2r9hZQ6qxrMzjdCCdmTHdWjc75/TD+QZkAxoZZD0zIW1eSVdrJjI5HjPzWA2zczJ0zq11s2NOP51vQDagkUHWMxPS5pVMas1EuuNJNU8idoTZORk6Z/R4WjMzv5k5/Xa+AdmAy68ReJnmfhjNkTGaJ2Inv+XItB639docrjuu+/+YWU125fkken7Wbq9KujaxeezOd7GaI5PsPEt1/nl5bmYT1tHfyJE5hUYmuznVnPghPC+TFzI7ggMTzSMp6dzJ1mbelefqtK65GdVk9QUkk1C9A0fq1btrrhSSPjnakHGTYIaZZiTVcaRqgqS2jSxBe8b54WccqdHInEIjk72cCrkLUnie2/ywNn6owQ5mwvSSbRu0Y/datpxD2Y5APGQ1p0LughSe5zY/rI0farCDlTC9VI8F4di9li3nED5DI4NAcirkLkjheW7zw9r4oQY72BHg11pQjt1r2XIO4TM0Mggkp0LLghSe5zY/rI0farCDH9aovcqWcwifoZFBIDkVWhak8Dy3+WFt/FCDHfywRu1VtpxD+AyNDALJqZC7IIXnuc0Pa+OHGuxgR4Bfa0E5dq9lyzmEz9DIIJCcCrkLUnie2/ywNn6owQ5WwvQS7ROkY/datpxD+AyNDAJr0rAilU0dqUiriPxIOM/S5ZNOjZsN/LA2fqjBDsmOoyicp29fOqjNr36IhPO0eOpILc6CY/datpxDOIkcGQSeU+mcpH4m59bapArrsyMM0GwNyeYxs23r4xg14DRt+eBwyuOzsv52hilmMk+sxthx+enniZ9xfyMQ7xQaGSCYzP76BKdrSDafHdva/esTMjket+bpEJJaRrOQoot0aGROoZEBgsdI6q1TaaxG0l/t2DYZu44vk3lDLs0Tm0s2zIfsRbIvgEAymnrrRBqrkfRXu7ZNxo7jMzKvW/OQogu70MgA8BUzqbd2p7EaSX+1c9tMxjAj03ndmseu+QBJ6uR1AQDQkpVEVbvSWJ1If7WjNrNjOJVw7fV+gEQjA8BnrCSq2pXG6kT6qx21mR3DqYRrr/cDJD5aAuAzZlJv7U5jNZL+aue2yVg9vti8bs1jJOSPFF1YRSMDwFeMpt46kcZqJP3Vrm2TseP4YvNmsrcd80jpj48UXdiFRgZwWWNTVOW7Duq5io9UvusgV2wkkCx5tUeXnOYsmZjCglyVjhuihn802bqe40siKh13lsL5ief79Hij/uflv2n5mx8pnN9Zj32jbb3h/ByVjhui8SWRtMeWKtE3domylXMnNm+yd2YiKdbRyLzJjq91rxIJ5+mxb1ygcH5nfhZgCTkygIuMhKYhdbLvgSP12vPJMf1m815V1dq7ngnD+PJzdMmQ3npjz+G4+VrOO+/Kc/X+gTo99epuVX+aOrQvWapssvvtOncSJfvuPZh8HSWZmjddsu/huuO6/4/8LCA5AvFOoZGBXxgJTUN6Tq2n0cC6TDhVkx3nTqqxk62Bn48H2YNAPMBHjISmIT2n1tNMYF0mnKrJ6rmTydiJ+PV40D7RyAAuMBKahvScWk+zgXWZcKomK+eOleP14/GgfaKRAVzgRMBae+bUerqx/k7VZKZ2L0L6+FmA3WhkABc4EbDWnjm1nm6sv1M1mandi5A+fhZgNxoZwAVGQtOQnlPraTawLhNO1WTl3LFyvH48HrRPNDKAC4yEpiE9p9bTTGBdIm7VZPXcyWRsu+flZwF24/JrwEVBz5FJlnHiFbuzVdJlnFx1XpFWvFWZ8MuqqXJXIgW5uuHC/hrYu6vhdYvVtnZ7lf5Q8bEO1R03NW6q5y7ROvbsmqOvnv85FeR3TpoxM74kYvp8MPLctay9d9dcKSQdqK1vzsGJFLTNqWlZi9fnrdfzO8Xp4yJH5hQaGfhNUP9S82sTZnU9kx3XvCvP1Wldc5MG1rUMlIsUJH/htBLaZ7bBaD1uJs+d0YZp7fYqy+dDJs9dotoT6RCSWl6xbTXQzy5+/bmxyo3jyopG5r777tOCBQvi7jv77LP13nvvZTwGjQxgXbYGmDl9XFbGtyuozkgNmW7r1vlgJZzQyUC/TPFzY03WBOINHTpUlZWVzbdXXnnF65KAdiVbA8ycPi4r49sVVGekhky3Pf6PJlfOB6vhhE4F+mWKnxv3jsv3jUynTp0UiUSab7179065fUNDg2pra+NuAMzL1gAzp4/Lyvh2BdUZqSHTbf9f+R5Xzgcnwwkl589bfm7cOy7fNzLvv/+++vXrpzPOOEM33nij9u7dm3L7hQsXKhwON9+Ki4tdqhTITtkaYOb0cVkZ366gOiM1ZLrtB4eOZTymFW6dT07Nw8+Ne8fl60bmoosu0pIlS7R69WqVlZVp9+7d+tKXvqQjR44k3Wfu3Lmqqalpvu3bt8/FioHsk60BZk4fl5Xx7QqqM1JDptsO6Nkl4zGtcOt8cmoefm7cOy5fNzKTJ0/W17/+dY0YMUITJ07UCy+8oOrqav3ud79Luk9ubq4KCgribgDMy9YAM6ePy8r4dgXVGakh022njR7oyvngZDih5Px5y8+Ne8fl60amtR49euiss87Szp07vS4FaDeyNcDM6eOyMr5dQXVGash0286dOrhyPlgNJ3Qq0C9T/Ny4d1yBamSOHj2qXbt2qagoeJerAUE2aViRyqaOVCQc/3ZxJJwX2EtIJeePy8r4qfZdPHWkFmc4rpEaMt3WrfMh2TyJtH7dNLpOTuDnxh2+zpH57ne/qylTpmjAgAH6+OOPNX/+fFVUVGj79u3q06dPRmOQIwPYJ6hhfunqzuS4rBy7lfFT7ZssnK9lym2yJNxU40vK6FjNrEmq+ZIdR8v7+3Yzn+wbG6dHl86qPpZ632Q1mXneM3l+Pjna4Mh5ZeW5snN+MzJ9/e5k24wO+PDDD3XDDTfo4MGD6tOnjy655BJt2rQp4yYGgL06dghp9OBeXpdhSCYJpOmOy2qKqdXxk+3bsUNINZ8e14NrdsTt2zrltqXYuB07hBw/rtYSzdejS44kqfrYiTbbJ0vrTVab0XVKNk+qmqykMsd+zcUdv6tIegmz1fHNJjub4Ze/D3z9jowdeEcGaL/sSCANYvpvKrFx/+3SQXripd2upc5aSemNMVubHXO3rsPu56Xl2DI5vtVkZ7/JmmRfADDDjgTSoKb/phI9dXvy5bZNTCbzmmE1pTfGTG12zd2a3c9LjB2pz2aSnYOMRgZAVrIjgTRb03+l5B89pZvXDDtTeo3W5kRCsJPPi9XxzSY7B5mvvyMDAGbZkUCa7em/Zub1chwzYzq5Tk4/L1bGN1JH0NKFW6ORAZCV7Eggzfb0XzPzejmOmTGdXCennxcr4xupI2jpwq3x0RKArGRHAmm2pv9KJ6/UcSud1c6UXqO1xea2k5PPi9XxzSY7BxmNDICsZEcCaVDTf1MJnbrd8qVBpuY1w2pKb4yZ2mJz2512ZPfzEmNH6rOZZOcg4/JrAFnNjgwNu3I4kgWIWRk/0b4hKelVMy3HTTevkZA+KX2AnNEcmdbHkWhNMg1lSzR3TOscmS6dO6pDKKSjDf9os204r5PGlxRq9ODecaF6rYME93xyTL/ZvFdVtW1zZFa8VZn0S7g98nM0Y8xAzbzszKQBf8mOJ1KQqxsu7K/+PbvEBQi+vueQlvx5j6o/PZFy20TH0bd728DAVOGDdsr09ZtGBkDWsyOB1OoYZpuGTGt79MWdeurV3XEvVj3yczT94gG6cFCvpMmxRpurRC/ERgLk0jVBr+78RGvfPaCaVscxY8xAzb5iSNoX81QNYKqk3bINu9qsX8+uOfrKef10pP4fbWpqLdEaxJqFgb27Jm0Ee3dN3GhkEgbYcpxEjVPCOvNzdMmQ3npjz+GE2yY6jta1GA0qNItG5hQaGQBe83OonpHxzEoXIJduXqeD3lKN5XSAX6ZrnWosu58vo5wK1yMQDwB8wM+hekbHs8LO8EE7jzmTscyyM9Qw2VhOPV9GeB2uRyMDAA7yc6iemfHMsDt80M5jduJ4M63F6NyJxnK6/kx5Ga5HjgwAOMjPoXp21mF1bK+C3twKg7MzPK/lfn4Ls/OiHhoZAHCQn0P17KzD6theBb25FQZnZ3hey/38FmbnRT18tAQADvJzqJ6Z8cywO3zQzmN24ngzrcXo3InGciLwzwwvw/VoZADAQX4O1TM6nhV2hg/aecyZjGWWnaGGycZyKvDPCK/D9WhkAMBhk4YVqWzqSEVa/cs5Es6z5ZJVu8dPOl5Brv55RJF65OfE3d+jS05z/khrkYJclY4booZ/NKl818GUV7UYOY5k2xYamC+TeRdPHanFU0emfdcj0RqE83NUOm6IxpdEDM/duh8I5+fo1rFnqntejpZv/VD/8/LftPzNj1S+66DGl0RUlkGN0sl3Tb596aCk23bp3FHdcuO/ddK6NWldm13nsVnkyACAS+wI5nNz/HSBay2D6qS2yb57DyZOt00XnmbkONLVaCSsLdW8rcP0enTpnDDZN1EwoZljThbSl0hs/PElkTbPQd9uuVJIbQIRWx5PohDCbrkdFY1Kdccbm+/LNHXYLgTinUIjAwDWmAmfczoE0I4a7eZGSF8ymYYOWpnHzbWUCMQDANjATPic0yGAdtRoN7dC+lIxeoxG5/E6+C4ZGhkAQFJmwuecDgG0o0a7eR3SZ+YY3ZrHaTQyAICkzATVuR2+5/ewP6PbWanT7sBAJ/a1G40MACApM0F1bofv+T3sz+h2Vuq0OzDQiX3tRiMDAEjKTPic0yGAdtRoN69D+swco1vzOI1GBgCQlJnwOadDAO2o0W5uhfSlYvQYjc7jdfBdMjQyAICUzATuOR0C6PV8TteQbKxEiiwcY7J5EgX8eR18lww5MgCAjJgJ3HM6BNCP80lKWIPV9evd9WS43YHa+ubAu0jByfA8qwF1mRyHHfMYQSDeKTQyAAAnrN5WqQUrt8ddwpwsxdfItk7VYIVb87REI3MKjQwAwG5GUnydSh12K83Yq9Rkkn0BAHCAkRRfp1KH3Uoz9kNqcjo0MgAAGGAkxdep1GG30oz9kJqcTqf0mwAAgBgnkoSNJuW6lWbsh9TkdGhkAAAwwIkkYaNJuW6lGfshNTkdPloCAMAAIym+TqUOu5Vm7IfU5HRoZAAAMMBIiq9TqcNupRn7ITU5HRoZAAAMMpLia2Tbxqaoyncd1HMVH6l818GUVwOZSRI2Mr6VedwUiByZxx57TA899JCqqqp03nnn6ec//7kuvPDCjPYlRwYA4BQjab3ptjUbOpdpDVZD7dxOTc6aQLzf/va3uummm7R48WJddNFF+tnPfqbf//732rFjh/r27Zt2fxoZAIDfOR0651WonRVZE4j3k5/8RLfccotmzJihkpISLV68WF26dNEvf/lLr0sDAMAyp0PnghBqZ4WvG5njx49ry5YtGjduXPN9HTp00Lhx41ReXp5wn4aGBtXW1sbdAADwK6dD54IQameFrxuZTz75RI2NjSosLIy7v7CwUFVVVQn3WbhwocLhcPOtuLjYjVIBADDF6dC5IITaWeHrRsaMuXPnqqampvm2b98+r0sCACApp0PnghBqZ4Wvk3179+6tjh07av/+/XH379+/X5FIJOE+ubm5ys3NdaM8AAAsi4XOVdXUJ/weS0gnL3U2Gzrn9Phe8/U7Mp07d9aoUaO0bt265vuampq0bt06jR492sPKAACwh9Ohc0EItbPC142MJN1xxx168skn9fTTT+vdd9/VzJkzVVdXpxkzZnhdGgAAtnA6dM7voXZW+D5HRpIeffTR5kC8888/X4888oguuuiijPYlRwYAEBROh865HWpnRdYE4llFIwMAQPBkTSAeAABAMjQyAAAgsGhkAABAYNHIAACAwKKRAQAAgUUjAwAAAotGBgAABBaNDAAACCwaGQAAEFi+/u3XdogFF9fW1npcCQAAyFTsdTvdLyDI+kbmyJEjkqTi4mKPKwEAAEYdOXJE4XA46eNZ/7uWmpqa9PHHH6t79+4Khbz7xVi1tbUqLi7Wvn37+J1PYj1aYi3isR7xWI94rMdnsn0totGojhw5on79+qlDh+TfhMn6d2Q6dOig008/3esymhUUFGTlCWcW6/EZ1iIe6xGP9YjHenwmm9ci1TsxMXzZFwAABBaNDAAACCwaGZfk5uZq/vz5ys3N9boUX2A9PsNaxGM94rEe8ViPz7AWJ2X9l30BAED24h0ZAAAQWDQyAAAgsGhkAABAYNHIAACAwKKRcdDChQv1hS98Qd27d1ffvn119dVXa8eOHV6X5Rs/+tGPFAqFVFpa6nUpnvnoo480depU9erVS/n5+Ro+fLjeeOMNr8vyRGNjo+bNm6dBgwYpPz9fgwcP1v3335/296xki5deeklTpkxRv379FAqF9Ic//CHu8Wg0qnvvvVdFRUXKz8/XuHHj9P7773tTrMNSrcWJEyd09913a/jw4eratav69eunm266SR9//LF3BTss3bnR0ne+8x2FQiH97Gc/c60+r9HIOGjjxo2aNWuWNm3apLVr1+rEiROaMGGC6urqvC7Nc6+//rr++7//WyNGjPC6FM8cPnxYY8aMUU5OjlatWqXt27frxz/+sU477TSvS/PEokWLVFZWpkcffVTvvvuuFi1apAcffFA///nPvS7NFXV1dTrvvPP02GOPJXz8wQcf1COPPKLFixfrtddeU9euXTVx4kTV19e7XKnzUq3FsWPHtHXrVs2bN09bt27Vs88+qx07duiqq67yoFJ3pDs3YpYvX65NmzapX79+LlXmE1G45sCBA1FJ0Y0bN3pdiqeOHDkSHTJkSHTt2rXRf/qnf4redtttXpfkibvvvjt6ySWXeF2Gb1x55ZXRb37zm3H3XXPNNdEbb7zRo4q8Iym6fPny5j83NTVFI5FI9KGHHmq+r7q6Opqbmxv9zW9+40GF7mm9Fols3rw5Kin6wQcfuFOUh5Ktx4cffhj93Oc+F922bVt0wIAB0Z/+9Keu1+YV3pFxUU1NjSSpZ8+eHlfirVmzZunKK6/UuHHjvC7FUytWrNDnP/95ff3rX1ffvn11wQUX6Mknn/S6LM9cfPHFWrdunf76179Kkt566y298sormjx5sseVeW/37t2qqqqK+5kJh8O66KKLVF5e7mFl/lBTU6NQKKQePXp4XYonmpqaNG3aNN11110aOnSo1+W4Lut/aaRfNDU1qbS0VGPGjNGwYcO8Lsczy5Yt09atW/X66697XYrn/va3v6msrEx33HGH/uM//kOvv/66br31VnXu3FnTp0/3ujzXfe9731Ntba3OOeccdezYUY2NjfrBD36gG2+80evSPFdVVSVJKiwsjLu/sLCw+bH2qr6+XnfffbduuOGGrP3FieksWrRInTp10q233up1KZ6gkXHJrFmztG3bNr3yyitel+KZffv26bbbbtPatWuVl5fndTmea2pq0uc//3n98Ic/lCRdcMEF2rZtmxYvXtwuG5nf/e53+vWvf62lS5dq6NChqqioUGlpqfr169cu1wPpnThxQtddd52i0ajKysq8LscTW7Zs0cMPP6ytW7cqFAp5XY4n+GjJBbNnz9bzzz+v9evX6/TTT/e6HM9s2bJFBw4c0MiRI9WpUyd16tRJGzdu1COPPKJOnTqpsbHR6xJdVVRUpJKSkrj7zj33XO3du9ejirx111136Xvf+56uv/56DR8+XNOmTdPtt9+uhQsXel2a5yKRiCRp//79cffv37+/+bH2JtbEfPDBB1q7dm27fTfm5Zdf1oEDB9S/f//mv1c/+OAD3XnnnRo4cKDX5bmCd2QcFI1GNWfOHC1fvlwbNmzQoEGDvC7JU2PHjtU777wTd9+MGTN0zjnn6O6771bHjh09qswbY8aMaXM5/l//+lcNGDDAo4q8dezYMXXoEP9vq44dO6qpqcmjivxj0KBBikQiWrdunc4//3xJUm1trV577TXNnDnT2+I8EGti3n//fa1fv169evXyuiTPTJs2rc33DSdOnKhp06ZpxowZHlXlLhoZB82aNUtLly7Vc889p+7duzd/lh0Oh5Wfn+9xde7r3r17m+8Hde3aVb169WqX3xu6/fbbdfHFF+uHP/yhrrvuOm3evFlPPPGEnnjiCa9L88SUKVP0gx/8QP3799fQoUP15ptv6ic/+Ym++c1vel2aK44ePaqdO3c2/3n37t2qqKhQz5491b9/f5WWluqBBx7QkCFDNGjQIM2bN0/9+vXT1Vdf7V3RDkm1FkVFRbr22mu1detWPf/882psbGz+u7Vnz57q3LmzV2U7Jt250bqRy8nJUSQS0dlnn+12qd7w+rKpbCYp4e2pp57yujTfaM+XX0ej0ejKlSujw4YNi+bm5kbPOeec6BNPPOF1SZ6pra2N3nbbbdH+/ftH8/LyomeccUb0+9//frShocHr0lyxfv36hH9fTJ8+PRqNnrwEe968edHCwsJobm5udOzYsdEdO3Z4W7RDUq3F7t27k/7dun79eq9Ld0S6c6O19nb5dSgabSexmQAAIOvwZV8AABBYNDIAACCwaGQAAEBg0cgAAIDAopEBAACBRSMDAAACi0YGAAAEFo0MAAAILBoZAK66+eabszJWH4A3+F1LAGwTCoVSPj5//nw9/PDDIlAcgF1oZADYprKysvn/f/vb3+ree++N+w3f3bp1U7du3bwoDUCW4qMlALaJRCLNt3A4rFAoFHdft27d2ny0dNlll2nOnDkqLS3VaaedpsLCQj355JOqq6vTjBkz1L17d5155platWpV3Fzbtm3T5MmT1a1bNxUWFmratGn65JNPDNUbCoX0i1/8Ql/96lfVpUsXDRkyRCtWrGh+/PDhw7rxxhvVp08f5efna8iQIXrqqacsrREAe9HIAPDc008/rd69e2vz5s2aM2eOZs6cqa9//eu6+OKLtXXrVk2YMEHTpk3TsWPHJEnV1dW64oordMEFF+iNN97Q6tWrtX//fl133XWG516wYIGuu+46vf322/ryl7+sG2+8UYcOHZIkzZs3T9u3b9eqVav07rvvqqysTL1797b12AFYw2+/BuCIJUuWqLS0VNXV1XH333zzzaqurtYf/vAHSSffkWlsbNTLL78sSWpsbFQ4HNY111yjX/3qV5KkqqoqFRUVqby8XF/84hf1wAMP6OWXX9aaNWuax/3www9VXFysHTt26KyzzsqoxlAopHvuuUf333+/JKmurk7dunXTqlWrNGnSJF111VXq3bu3fvnLX1pcDQBO4R0ZAJ4bMWJE8/937NhRvXr10vDhw5vvKywslCQdOHBAkvTWW29p/fr1zd+56datm8455xxJ0q5du0zP3bVrVxUUFDTPM3PmTC1btkznn3++/v3f/11//vOfzR0gAMfwZV8AnsvJyYn7cygUirsvdjVUU1OTJOno0aOaMmWKFi1a1GasoqIiy3PH5pk8ebI++OADvfDCC1q7dq3Gjh2rWbNm6b/+678MzQHAOTQyAAJn5MiR+r//+z8NHDhQnTo5+9dYnz59NH36dE2fPl1f+tKXdNddd9HIAD7CR0sAAmfWrFk6dOiQbrjhBr3++uvatWuX1qxZoxkzZqixsdG2ee69914999xz2rlzp/7yl7/o+eef17nnnmvb+ACso5EBEDj9+vXTq6++qsbGRk2YMEHDhw9XaWmpevTooQ4dTv61tmTJkrQBfel07txZc+fO1YgRI3TppZeqY8eOWrZsmR2HAMAmXLUEICvNnz9fGzdu1IYNG7wuBYCD+I4MgKy0atUqPfroo16XAcBhvCMDAAACi+/IAACAwKKRAQAAgUUjAwAAAotGBgAABBaNDAAACCwaGQAAEFg0MgAAILBoZAAAQGDRyAAAgMD6/58/49JsUCClAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.scatter(correct_df['Time'], correct_df['Count'])\n",
"plt.ylabel('Count')\n",
"plt.xlabel('Time, ns')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"id": "SCrC88HHpDFt"
},
"outputs": [],
"source": [
"n = 15\n",
"df_avg = pd.Series(correct_df['Count']).rolling(window=n).mean().iloc[n-1:].values"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 447
},
"id": "11Dml6iApz8n",
"outputId": "394a6c70-bd3c-4c44-ae3f-2fa0fef081ff"
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7ae94ba16f20>]"
]
},
"execution_count": 95,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABFNElEQVR4nO3dd3xV9eH/8dcdyc1OyE5IAgkj7L1BREGUKkjFPWu1juK2ftW2rrZKtb/WVmtBra1VRK1VHNRR9lD2RjYEAtkh42aQec/vj4RoJEACyT03N+/n43EfNfeee+87t8h9+zmf8/lYDMMwEBEREXETq9kBREREpGNR+RARERG3UvkQERERt1L5EBEREbdS+RARERG3UvkQERERt1L5EBEREbdS+RARERG3spsd4IdcLheZmZkEBwdjsVjMjiMiIiLNYBgGJSUlxMfHY7WefmzD48pHZmYmiYmJZscQERGRs3DkyBESEhJOe4zHlY/g4GCgLnxISIjJaURERKQ5nE4niYmJDd/jp+Nx5ePEqZaQkBCVDxERkXamOVMmNOFURERE3ErlQ0RERNxK5UNERETcSuVDRERE3ErlQ0RERNxK5UNERETcSuVDRERE3ErlQ0RERNxK5UNERETcSuVDRERE3ErlQ0RERNxK5UNERETcSuUDMAyD/2w8ylffZpsdRURExOt53K627lZeVcMj/9nGf7dlYbda2PzkRQT7+ZgdS0RExGt1uPJRUlHNx1sy2ZxeCAZszyhmX24pADUug28znYxKiTA5pYiIiPfqMOXjSEE5s5cf4JPNGZRV1TZ6LCLQl+gQP3ZlOdl+tFjlQ0REpA11mPJRUlHDvLXpAHSLCmTawM4E+NrwsVmY0j+O/2w8yq4sJ9syik1OKiIi4t1aVD5mzZrFRx99xO7du/H392fMmDE8//zzpKamNhwzYcIEli9f3uh5d955J3PmzGmdxGepT3wI917YnTHdIhmVEo7FYmn0+ICEUAC2Hy0yIZ2IiEjH0aKrXZYvX87MmTNZs2YNCxcupLq6msmTJ1NWVtbouJ/97GdkZWU13F544YVWDX22Hp6cyuhuEScVD4D+nevKx6Fj5RSXV7s7moiISIfRopGPL7/8stHPb775JtHR0WzcuJHx48c33B8QEEBsbGzrJHSTsABfksIDSC8oZ3tGMeN6RJ50zIZDBcxedoCnp/UlMTzAhJQiIiLt3zmt81FcXDc/Ijw8vNH977zzDpGRkfTr14/HH3+c8vLyU75GZWUlTqez0c0s/etPvWzLKGry8Re+2sPi3bm8uuKAG1OJiIh4l7MuHy6XiwceeICxY8fSr1+/hvuvv/565s6dy9KlS3n88cd5++23ufHGG0/5OrNmzSI0NLThlpiYeLaRztmAzifmfZw86bSwrIoNhwoAWLQzF8Mw3JpNRETEW5z11S4zZ85kx44drFq1qtH9d9xxR8M/9+/fn7i4OCZOnMiBAwfo1q3bSa/z+OOP89BDDzX87HQ6TSsgDSMfTZSPJbtzcdX3jWxnBd9mOulXX1ZERESk+c5q5OOee+5hwYIFLF26lISEhNMeO3LkSAD279/f5OMOh4OQkJBGN7OcmHSaUXScY6WVjR5btCsHAGv9XNX/7cxxazYRERFv0aLyYRgG99xzD/Pnz2fJkiUkJyef8TlbtmwBIC4u7qwCulOwnw8pUYEAjdb7qKiuZfnePABuGtUFgEUqHyIiImelReVj5syZzJ07l3nz5hEcHEx2djbZ2dkcP34cgAMHDvDb3/6WjRs3cujQIT799FNuvvlmxo8fz4ABA9rkF2htgxM7AbDm4LGG+9YcPEZ5VS0xIQ7undgDqwV2ZjnJKDpuVkwREZF2q0XlY/bs2RQXFzNhwgTi4uIabu+//z4Avr6+LFq0iMmTJ9OrVy8efvhhZsyYwWeffdYm4dvCuB51S6uv2pffcN/C+lGOib1jiAxyMLRLXUFZvEujHyIiIi3VogmnZ7rCIzEx8aTVTdubsd3r1vf4NtNJfmklnQJ8G+Z7XNQnpuF/1x8q5Ivt2dw8uqtZUUVERNqlc1rnwxtFB/vRO65u0uvX+/NZsTePHGclof4+jK7fcG5KvzhsVgurDx5jm5ZjFxERaRGVjyacV7+66cp9+cxdcxiAK4cm4OdjAyAxPIDLB8UD8PKSpq/iERERkaapfDThRPlYtCuHpXtyAbh+ZFKjY2Ze0B2LpW4+yM5M81ZlFRERaW9UPpowvGs4vnYrReXVuAwY0y2CblFBjY7pFhXEZQPqRj/+unSfGTFFRETaJZWPJvj52BiZ/N1+NTeM7NLkcfdc0B2AL3Zkc6Tg1PvXiIiIyHdUPk7hxKmXqGAHk/vGNHlMamwww7t2wjDgmwP5TR4jIiIijal8nMLVwxKZ0i+W303vh4/t1B/TyOS6K2DWHixwVzQREZF2TeXjFMICfJl941Au7ht72uNGptSdnlmbpvIhIiLSHCof52hIUidsVgsZRcc5Wqh5HyIiImei8nGOAh32ht1w12n0Q0RE5IxUPlrBiStjNO9DRETkzFQ+WsGJeR/rDql8iIiInInKRysY2iUciwXS8svIdVaYHUdERMSjqXy0glB/H/rUb0b38ZYMXK7T7/4rIiLSkal8tJITO94+9/luLvjjMhbtzDE5kYiIiGdS+Wgl91zYnZ+M6Uqww87hY+Xc995myqtqzI4lIiLicVQ+WklYgC9PT+vL2l9NJCk8gPKqWhbtyjU7loiIiMdR+WhlAb52pg2s2+320y2ZJqcRERHxPHazA3ijaYPi+evS/Szfm0tReRV5JZU8/+UeCsurABjTLYKHJ6eanFJERMQcKh9toGdMML1ig9mdXcI7a9OZtzadjKLjDY9vPFzIgIQwLurT9G65IiIi3kynXdrItEF1p17+8NUeMoqO0yUigDk3DuWaYYkAPP3pt5qQKiIiHZLKRxuZOiC+4Z+DHHb+fvMwLukXy1PT+tA5zJ+MouP8dcl+ExOKiIiYQ+WjjSSGBzAhNQq71cKL1wyiR0wwUDch9elpfQF4feVB0vLLzIwpIiLidprz0Ybm3DgU5/FqokP8Gt1/UZ8YxnSL4JsDx1i4M5s7xnczKaGIiIj7aeSjDfn52E4qHiec1yMKgE2Hi9yYSERExHwqHyYZkhQGwKb0QgxDe8GIiEjHofJhkgEJYditFnJLKsks1k64IiLScah8mMTf10bv+p1wNx0uNDmNiIiI+6h8mOj7p15EREQ6CpUPEw3p0gmATelF5gYRERFxI5UPEw1JqisfOzOLqaiuNTmNiIiIe6h8mCihkz+RQQ6qaw12ZBSbHUdERMQtVD5MZLFYNO9DREQ6HJUPkw2uP/Wy4ZDKh4iIdAwqHyYblRIOwJqDx6ipdZmcRkREpO2pfJhsQEIYIX52nBU1bNO8DxER6QBUPkxms1oY2z0SgFX78k1OIyIi0vZUPjzAiU3mVu7LMzmJiIhI21P58ADn9agb+dicXkRJRbXJaURERNqWyocHSAwPoGtEADUugzUHC8yOIyIi0qZUPjzEiVMvq3TqRUREvJzKh4cYV3/qZaUmnYqIiJdT+fAQo7tFYLXAwfwy8koqzY4jIiLSZlQ+PESInw+J4QEA7MstMTmNiIhI21H58CA9ooMAOJBbanISERGRtqPy4UG6RwcDsE/lQ0REvJjKhwfpXj/ysS9H5UNERLyXyocHOXHaZX+eyoeIiHgvlQ8P0q2+fOSVVFJUXmVyGhERkbah8uFBghx2Oof5A7Bf8z5ERMRLqXx4mBOjHyofIiLirVQ+PMyJeR/7ckvJKDrOlL+s5OZ/rOOrb7OpqXWZnE5EROTc2c0OII19v3z8ZdFedmU52ZUFK/bm0a9zCB/cOQZ/X5vJKUVERM6eRj48zInLbbekF/LRpgwArhmWSJDDzo4MJ6v2a+8XERFp31Q+PMyJ8uGsqKHGZTCueyTPXzmAHw/uDMBK7XorIiLtnMqHhwkL8CUq2NHw870Xdge+2/V2lXa9FRGRdk7lwwOdmPcxIjmckSkRQN2utzarhYP5ZRwtLDcznoiIyDlR+fBAVw9LJDkykF9f2rvhvhA/HwYlhgEa/RARkfZN5cMDTR/cmaW/mMCAhLBG959Xf+plpcqHiIi0Yyof7ciJ8vH1gXxqXYbJaURERM6Oykc7MjAhjGCHnaLyanZkFJsdR0RE5Ky0qHzMmjWL4cOHExwcTHR0NNOnT2fPnj2NjqmoqGDmzJlEREQQFBTEjBkzyMnJadXQHZXdZmV0t7oJqFrvQ0RE2qsWlY/ly5czc+ZM1qxZw8KFC6murmby5MmUlZU1HPPggw/y2Wef8cEHH7B8+XIyMzO54oorWj14RzW2e92plzUHj5mcRERE5OxYDMM468kDeXl5REdHs3z5csaPH09xcTFRUVHMmzePK6+8EoDdu3fTu3dvVq9ezahRo874mk6nk9DQUIqLiwkJCTnbaF5rV5aTKX9ZSYCvja1PTcbHpjNnIiJivpZ8f5/TN1dxcd28g/DwcAA2btxIdXU1kyZNajimV69eJCUlsXr16nN5K6mXGhNMqL8P5VW1fJvpNDuOiIhIi511+XC5XDzwwAOMHTuWfv36AZCdnY2vry9hYWGNjo2JiSE7O7vJ16msrMTpdDa6yalZrRaGd60re2t16kVERNqhsy4fM2fOZMeOHbz33nvnFGDWrFmEhoY23BITE8/p9TqCUSn15SOtwOQkIiIiLXdW5eOee+5hwYIFLF26lISEhIb7Y2NjqaqqoqioqNHxOTk5xMbGNvlajz/+OMXFxQ23I0eOnE2kDmVEcl35WH+oQOt9iIhIu9Oi8mEYBvfccw/z589nyZIlJCcnN3p86NCh+Pj4sHjx4ob79uzZQ3p6OqNHj27yNR0OByEhIY1ucnp94kIIctgpqahhV5ZOU4mISPtib8nBM2fOZN68eXzyyScEBwc3zOMIDQ3F39+f0NBQbrvtNh566CHCw8MJCQnh3nvvZfTo0c260kWax26zMqxrJ5btyWNtWgH9OoeaHUlERKTZWjTyMXv2bIqLi5kwYQJxcXENt/fff7/hmBdffJHLLruMGTNmMH78eGJjY/noo49aPXhHd+LUy7o0TToVEZH2pUUjH81ZEsTPz49XXnmFV1555axDyZmNSqlb6XT1gWNU17q03oeIiLQb+sZqpwYmhNEpwAdnRQ3rD+mqFxERaT9UPtopm9XChb1iAFi0M9fkNCIiIs2n8tGOXdSnrnws3JXdrFNiIiIinkDlox07r0ckvnYrRwqOsy+31Ow4IiIizaLy0Y4FOuyM7VY38XThzhyT04iIiDSPykc7N+nEqReVDxERaSdUPtq5Sb3ryseWI0XkllSYnEZEROTMVD7auZgQP/rE1S1Jv+FQoclpREREzkzlwwsMTAwDYNvRYnODiIiINIPKhxcYmFC3t8v2jCJzg4iIiDSDyocX6F9fPrYdLdZ6HyIi4vFUPrxAz5hgfO1WSipqOHys3Ow4IiIip6Xy4QV8bNaGSadbjxaZG0ZEROQMVD68xIAT8z406VRERDycyoeXGJAQBsC2DJUPERHxbCofXuLEyMe3GcXUujTpVEREPJfKh5foFhWEv4+Nsqpa0vK1yZyIiHgulQ8vYbNa6Ne5ftLpEZ16ERERz6Xy4UX6dw4D6vZ5ERER8VQqH15kRHI4AF8fyDc5iYiIyKmpfHiR0d0isFrgYF4ZGUXHzY4jIiLSJJUPLxLq78Og+k3mVu3LMzeMiIjIKah8eJlxPaIAWLlPp15ERMQzqXx4mfE9IgH4en8+Lq33ISIiHkjlw8sMTAwjyGGnsLyabzOdZscRERE5icqHl/GxWRndLQKAFZr3ISIiHkjlwwudV3/qZZXmfYiIiAdS+fBCY7vXlY+N6YXU1LpMTiMiItKYyocXSo4IJMDXRlWNi0PHysyOIyIi0ojKhxeyWi30iAkGYE+2NpkTERHPovLhpXo1lA9d8SIiIp5F5cNLpcbWlY/d2SUmJxEREWlM5cNL9aovH3tyVD5ERMSzqHx4qZ715SO9oJzyqhqT04iIiHxH5cNLRQY5iAzyxTBgX44mnYqIiOdQ+fBiJ+Z97NG8DxER8SAqH16sZ4zmfYiIiOdR+fBivTTyISIiHkjlw4ulxoYAutxWREQ8i8qHF+sRHQRAfmklx0orTU4jIiJSR+XDiwU67CSFBwCwPaPY5DQiIiJ1VD683Ikdbt/85pC5QUREROqpfHi5u85PwWa1sGxPHluPFJkdR0REROXD23WJCOTygfEAvLxkv8lpREREVD46hJkXdsdigUW7ctiZqV1uRUTEXCofHUC3qCAuG1A3+jF7+QGT04iISEen8tFB3Dk+BYAvd2SRV6LLbkVExDwqHx1Ev86hDEoMo7rW4IONR8yOIyIiHZjKRwdyw8gkAOatTafWZZicRkREOiqVjw5k6sB4QvzsHC08zop9eWbHERGRDkrlowPx87Fx5dBEAN5Zc9jkNCIi0lGpfHQw19efelmyO5eSimqT04iISEek8tHBdI8OIi7UD5eh3W5FRMQcKh8dUJ+4EAAtOCYiIqZQ+eiA+sSrfIiIiHlUPjqghpGPLJUPERFxP5WPDqhvfCgAe3JKqK51mZxGREQ6GpWPDiihkz/BDjtVNS4O5pWZHUdERDoYu9kBxP2sVgu940JYd6iAnVnFpMYGn3TM2oPHeHHRXqprDSzANcMTuWpYovvDioiI11H56KD6xNeXj0wnPx7c+LGyyhoefH8LmcUVDfcdyCtlxpAErFaLm5OKiIi30WmXDup0k05fWryPzOIKEjr5M+fGofj72Cgsr2Z/Xqm7Y4qIiBdS+eigTlxu+22mE8P4bpO5PdklvLEqDYBnpvXlkn6xDO3SCag7FSMiInKuWlw+VqxYwdSpU4mPj8disfDxxx83evwnP/kJFoul0e2SSy5prbzSSrpHB2G3Wigqrybre6dXnv70W2pcBpP7xDCxdwwAI5LDAViTVmBKVhER8S4tLh9lZWUMHDiQV1555ZTHXHLJJWRlZTXc3n333XMKKa3Pz8dG9+ggoG70AyCr+DirDx7DaoEnp/ZpOHZkfflYl1bQaJRERETkbLR4wumUKVOYMmXKaY9xOBzExsaedShxj0GJYezOLuHr/flc1CeGlfvyAeifEEZCp4CG4wYmhuFrt5JXUklafhkpUUFmRRYRES/QJnM+li1bRnR0NKmpqdx9990cO3bquQKVlZU4nc5GN3GPC3tFA7BoVw6GYbCqvnyM7xHZ6Dg/HxuDEsMAWKtTLyIico5avXxccsklvPXWWyxevJjnn3+e5cuXM2XKFGpra5s8ftasWYSGhjbcEhO1loS7jOsRicNu5WjhcXZllbBqf135GNc98qRjR33v1IuIiMi5aPXyce211zJt2jT69+/P9OnTWbBgAevXr2fZsmVNHv/4449TXFzccDty5EhrR5JTCPC1c179KMfLS/ZRUFZFoK+NwUmdTjp2RHIEUHfFi+Z9iIjIuWjzS21TUlKIjIxk//79TT7ucDgICQlpdBP3mVR/RcsXO7IBGJUSga/95D8WQ7qEYbdayCyuIL2g3K0ZRUTEu7R5+Th69CjHjh0jLi6urd9KzsKFvaMb/Xxej5NPuUDdKMnwrnWnXj7fnt3muURExHu1uHyUlpayZcsWtmzZAkBaWhpbtmwhPT2d0tJSHnnkEdasWcOhQ4dYvHgxl19+Od27d+fiiy9u7ezSCqKD/RomkwKM6xF1ymOnDowH4NOtmW0dS0REvFiLy8eGDRsYPHgwgwfXbQjy0EMPMXjwYJ588klsNhvbtm1j2rRp9OzZk9tuu42hQ4eycuVKHA5Hq4eX1nFRn7pTL/GhfnSLCjzlcVP6xeJjs7Ary8m+nJI2y5NfWkmOs+LMB4qISLvU4nU+JkyYcNoJh1999dU5BRL3u2Z4ImsOHmP6oM5YLKfeOK5ToC/je0SxeHcun27N5OHJqa2aY11aAf9YlcbCXTn42a0sevh84kL9W/U9RETEfNrbRYgMcvD2bSOZMTThjMdOG/TdqZfWvOplc3oh1762mi+/zabWZVBWVcvS3Xmt9voiIuI5VD6kRSb1jsHfx8bhY+VsPVrcaq/77w1HcRkwOiWCa4bVrfWycp/Kh4iIN1L5kBYJdNiZVD9H5P31rbMmS1WNi8+3ZwFw74XduWZEXfn4en8+tS6tKSIi4m1UPqTFbhrVBYAPNx4lu7jxxNCK6loWbMskv7Sy2a+3cl8exceriQ52MDIlggGdQwnxs+OsqGHb0aLWjC4iIh5A5UNabERyOCOSw6mqdTFn+YGG+7OKj3PVnNXcM28zP/rLSrYcKWr0PJfL4JMtGWQVH290/ydb6i7dvWxAPDarBbvNyphudeuNnNhvRkREvIfKh5yV+y7sAcC769LJdVaweFcOU1/+mu0ZdfNAcksqufrV1XyyJaPhOe+uT+f+97Zw0xvrqKpxAVBeVcPCnTnAd5NZAc7rWVc+Vqp8iIh4HZUPOStju0cwOCmMyhoXE/7fMm771wbySyvpFRvMlw+cx0V9YqiqcfHwv7dytLAcwzB4e/VhAPbnlvL3VQcBWLgzh+PVtXSJCGBgQmjD64+vX+xsU3ohpZU17v8FRUSkzah8yFmxWCzcN7Fu9KO8qpZgPzu3jUvmw7vH0Cs2hFdvHMqYbhHUuAxmLzvApvQidmd/tzDZS4v3MW9tOr+evwOAaQPjG60xkhgeQJeIAGpcBmsOHHPvLyciIm2qxYuMiZwwoWcUf7xqIAZwaf84/H1tDY9ZrRbun9iDbw4c44MNRxs2o7tyaALpBeWsSyvgl/O3AzC8ayduPy/lpNcf3yOKt48d5l+rDzGxd/RpF0ATEZH2QyMfctYsFgszhiZw5dCERsXjhJEpEQ0TU0/M3bhxVBd+N70fdmtdkbhuRBLv3D6KUH+fk55/27hkfO1WVu7L12Z2IiJeROVD2tSJiakAfeNDGJgQSs+YYN6/czRv3jqc537cD197038Mu0YG8vMJ3QD4zYJvKamodktmERFpWyof0qbGdo9gaJdOANwyumvDqZOhXToxIfXMp1LuOr8bXSICyHFW8tiH2zlSf/pGRETaL4vRmht0tAKn00loaCjFxcWEhISYHUdawbHSSjalFzHpLOdtLN+bxy3/WAeAxQKT+8Twh6sGEuJ38qmaM3FWVPPsgl1c2Duai/vGtvj5p/KXRftYd6huYmx0sB//d0mqNsUTkQ6lJd/fKh/SLizdk8sbK9NYtb9u7siFvaJ5/eZh2KwtKzMvfLmbvy07QLCfna8fu/CsCswP5TorGPHc4kb39YkL4cO7xzQ5F0ZExBu15Ptbp12kXbggNZq5t4/kw7vH4LBbWbI7lxe+2s3CnTn8/J2NPPj+FlYfOHbanXaLyqt4q36tkZKKGt765lCrZEvLLwMgJsTBn64eSESgLzuznDzyn62tuvOviIi3UPmQdmVol068cOUAAF5dfpCfvbWBz7dnM39zBte9voZLX1pFXknT+8r88+tDlFbWEFg/GvHGqjTKWmEBs0PH6spHamwIVwxJYPaNQ/GxWViwLYsH3t/CjozW2/1XRMQbqHxIu3P5oM7cXX8VTKcAH+4Yn8L1I5MI8LWxM8vJ31cePOk5JRXV/PPrNABmzRhA14gACsurmbvm8DnnScuvmwTbNSIAqNv75jeX9wPq9q257OVV3P6v9VTXus75vUREvIEWGZN26f8uTuWKwZ1JDA/Az6duJOOC1Gh+9tYG/r3hCA9N7onDbsMwDDYeLuRvyw7grKihe3QQl/WPo7K6lkf+s43XVhzkmuGJhAX4nnWWQ/WnXbpGBDbcd92IJLpHB/HW6sN8uSOLRbtyWbQzhyn9487tFxcR8QIa+ZB2yWKx0CMmuKF4AFyQGkVcqB+F5dV8sT2b41W1XPf6Gq6cs5olu3MB+MXkVKxWC9MHdyY5MpBjZVXcM28zNecwKnHitEtyZGCj+4d3Defl6wZzx/i61VvfWZt+1u8hIuJNVD7Ea9htVq4dngTA3DWHeeQ/W1lzsAA/HytXD0vg03vGckm/ustrfWxWXrl+CP4+Nlbtz2fWF7vP6j1dLqOhfHT9Qfk44drhSVgssGp/fsPkVBGRjkzlQ7zKNcMTsVktbDhcyIJtWfjYLLx920heuHIgAxLCGh3bJz6EP109EKibfPq/b1u+hHtOSQUV1S5sVgsJnZpe1yMxPIALUqMBmLf23OeYiIi0dyof4lViQ/24qHdMw8/PTOvH8K7hpzx+Sv84bh3bFYAPNh5t8fsdqp9smtjJHx/bqf91umFkUsN7PPqfbfR/+it+Vb+xnohIR6PyIV7n7gndCHLYubP+KpgzuXJoAgAr9+VxvKq2Re91plMuJ0xIjSY+1I+i8mre33CEkooa3lt/hGOlTV8WLCLizVQ+xOsMTAxj+9OTefxHvZt1fJ+4EDqH+VNR7eLr+hVUm6upK12aYrNaeHhyKuGBvkwbGE/36CBqXQafb89q0fuJiHgDlQ/xSi3ZQ8ZisTCpd92cjEW7clr0PmkN5SPgjMfOGJrApicu4qXrBnPt8EQAPt2a2aL3ExHxBiofIsCkPnXzRBbtysXlav6S6M097fJDlw2Ix2KB9YcKySg63qLnioi0dyofIsDI5AiCHXbySyvZcrSoWc9xuQwOH6ubcPrDNT7OJDbUj5HJdRNhP9Poh4h0MCofIoCv3cr5qVEAzF19mKW7c9mbU3La52Q7K6iscWG3Wugc1vRltqczbWBnoG4JdhGRjkTlQ6TeRfWnXj7anMGtb65n8osr+OP/9pzyNMyJyaZJ4QHYT3OZ7alM6ReL3WphV5aT/bmnLzoiIt5E5UOk3sV9Y7m0fxz9O4fSOy4EgJeX7OfOuRsp/cHut86KauasqNvArqXzPU7oFOjL+T3rRls+1eiHiHQgKh8i9fx8bLxywxA+u3ccX9x/Hn+8aiC+NisLd+Zw99yNDfu/HMwr5cevfM2KvXk47FZuH5d81u85bVA8AJ9szcQwmj/RVUSkPVP5EDmFGUMTePeOkfj72Fi5L5/nv9zNir15TH/law7klREb4scHd41mTPfIs36PSb1j8PexcfhYOduOFrdiehERz6XyIXIaQ7uE8/+uqtv/5fWVadzyz3U4K2oYkhTGp/eOPWm/mJYKdNgbLvP94cTTg3mlLNimERER8T4qHyJncOmAOO69sDsAhgFXDU3g3TtGER3s1yqvf/nAulMvC7ZlUls/ufVIQTlXzP6Ge+Zt5n87W7bwmYiIp7ObHUCkPXhwUk/CAnyJCPTl8kHxLVpB9UzG94wi1N+H3JJKvtiRxQWp0fzsrQ0UlVcD8Pbqw1zcN7bV3k9ExGwqHyLNYLVauO0cJpaejq/dyo/6x/HuunTumbeZTgE+FJZXExHoS0F5Fav255OWX9bihcxERDyVTruIeIDHpvTiisGd8bVbKSyvxtdm5bWbhzVcivvuunSTE4qItB6L4WGz2ZxOJ6GhoRQXFxMSEmJ2HBG3Kiir4vPtWfSIDmJkSgSLduZw+1sb6BTgw+rHJ+LnYzM7oohIk1ry/a2RDxEPEh7oy42jujAyJQKAC3pFEx/qR2F5NV/syDI5nYhI61D5EPFgNquF60YkAfD6ijRddisiXkHlQ8TD3TiqC4G+NnZmOVm8K9fsOCIi50zlQ8TDdQr05abRXQF4eck+jX6ISLun8iHSDtx+XjJ+Pla2Hi1m+d48s+OIiJwTlQ+RdiAyyMENI7sA8OdF+xpWQhURaY9UPkTaiTvHp+DvY2PLkSJe+HK32XFERM6ayodIOxEd4scLVw4A4NUVB5m/+ajJiUREzo7Kh0g7MnVgPD+f0A2ARz/czrajReYGEhE5CyofIu3MLyanMrFXNFU1Lu54ayO5zgqzI4mItIjKh0g7Y7Va+PO1g+geHUS2s4K75m6ksqbW7FgiIs2m8iHSDgX7+fD6zcMI8bOzKb2IZ/+7y+xIIiLNpvIh0k4lRwby0nWDAXh7zWG2HCkyN5CISDOpfIi0YxNSo7liSGcMA341f7vW/xCRdkHlQ6Sd++WPehPiZ+fbTCdvrz5kdhwRkTNS+RBp5yKDHPzfJb0A+OP/9lJRrcmnIuLZVD5EvMB1I5IID/SlpLKGvTklZscRETktlQ8RL2CzWugVGwzAnmyVDxHxbCofIl4iVeVDRNoJlQ8RL5EaU18+dNpFRDycyoeIlzgx8rFbIx8i4uFUPkS8RI/6kY+8kkoKyqpMTiMicmoqHyJeIshhJzHcH9C8DxHxbC0uHytWrGDq1KnEx8djsVj4+OOPGz1uGAZPPvkkcXFx+Pv7M2nSJPbt29daeUXkNFJjQgDYk+00OYmIyKm1uHyUlZUxcOBAXnnllSYff+GFF3jppZeYM2cOa9euJTAwkIsvvpiKCm37LdLWUmODAE06FRHPZm/pE6ZMmcKUKVOafMwwDP785z/z61//mssvvxyAt956i5iYGD7++GOuvfbac0srIqeVGls38qFJpyLiyVp1zkdaWhrZ2dlMmjSp4b7Q0FBGjhzJ6tWrm3xOZWUlTqez0U1Ezs6Jhcb2ZpdgGNpkTkQ8U4tHPk4nOzsbgJiYmEb3x8TENDz2Q7NmzeKZZ55pzRgiHVZyZCA+NgtlVbUcLTxOYnhAi1+jqsbF31cdZH9OKQARQb48PDkVPx9ba8cVkQ6qVcvH2Xj88cd56KGHGn52Op0kJiaamEik/fKxWekWFcTu7BK+zSwmMTwAwzCYu+YwafnlzBjamb7xoWw9UsTHWzIYlBjG5YM6Nzw/v7SSn8/dxLpDBY1eNybEj9vPS3H3ryMiXqpVy0dsbCwAOTk5xMXFNdyfk5PDoEGDmnyOw+HA4XC0ZgyRDm1QYhi7s0t49MPt2K1WPtmayWdbMwH4x9dpxIb4ke2smwBusdRdojuxdwzfZhZzx1sbySg6TrDDzh3jUzhcUM5/Nh5l3tp0bhuXjMViMfNXExEv0apzPpKTk4mNjWXx4sUN9zmdTtauXcvo0aNb861E5BQenpzKoMQwio9Xc/tbG/hsayZ2q4UJqVH42CxkOyvwtVvpExeCYcD9723h9RUHuXL2ajKKjpMcGcj8mWO5d2IPnp7WlyCHnYP5Zaw+eMzsX01EvESLRz5KS0vZv39/w89paWls2bKF8PBwkpKSeOCBB/jd735Hjx49SE5O5oknniA+Pp7p06e3Zm4ROYWoYAfv3TGKX87fzkebMggP9OVvNwxhVEoEeSWVbD1SxJAunQhy2Lnx72tZd6iAZz/fBcB5PSL563VDCA3wAepGRaYPjmfumnTeWZPOmG6RZv5qIuIlLEYLp8QvW7aMCy644KT7b7nlFt58800Mw+Cpp57itddeo6ioiHHjxvG3v/2Nnj17Nuv1nU4noaGhFBcXExIS0pJoIvI9hmGw4XAhyZGBRAY1fWozv7SSy//6NRlFx7l9XDKPTemF3dZ4QHRnppMfvbQSu9XCi9cM4osdWRw+Vg5AgK+NS/vHccXQBEL8fNr8dxIRz9WS7+8Wl4+2pvIh4l7Fx6vJLDpO77hT//t2xd++ZlN60Skf9/exMeuK/kwf3PmUx4iId2vJ97fpV7uIiLlC/X0I9T/9qMUd41O4a+4mgv3szBiSwPk9o7BaLaTllTJvXTp7c0p57vNdTBsYj9WqSakicnoa+RCRZjmUX0Z0iIMA38b/zVJZU8uQ3yykrKqWT2aOZWBimDkBRcRULfn+1q62ItIsXSMDTyoeAA67jfNTowBYtCvH3bFEpB1S+RCRc3ZRn7pVjRfuVPkQkTNT+RCRc3ZBajQ2q4Xd2SUcKSg3O46IeDiVDxE5Z2EBvgzr0gnQqRcROTOVDxFpFSdOvah8iMiZqHyISKs4UT7WHixg+d48k9OIiCfTOh8i0iq6RARyYa9oluzO5dZ/ruPRS3pxXo+oRsdEBTuICtZGkiIdndb5EJFWU1lTyxMf7+DfG442+bjDbuWrB8bTNTLQzclEpK1pnQ8RMYXDbuP5GQN4ZlpfksIDiA52NNwCfG1U1riYu+aw2TFFxGQa+RARt1iyO4efvrmBUH8f1v5yIn4+NrMjiUgr0siHiHic83tG0znMn+Lj1fx3W5bZcUTERCofIuIWNquF60cmAfDOWp16EenIVD5ExG2uGpaA3WphU3oRry4/wHvr0tmfW2J2LBFxM11qKyJuEx3sx8X9YvnvtixmfbEbgBA/O0t+MYHIIF2CK9JRaORDRNzq/y5OZerAeCb1jqFzmD/Oihqe+3yX2bFExI10tYuImGZzeiFXzP4Gw4D37hjFqJQIsyOJyFnS1S4i0i4MTurEdSPqJqE+9uE2fv/Fbv70vz0cLdTOuCLeTHM+RMRUj17ci692ZHPoWDlzlh8AYFd2Ca/fPMzkZCLSVlQ+RMRUoQE+vPGT4fx3WyallbW8uy6dFXvzOF5Vi7+vFiIT8UYqHyJiukGJYQxKDMMwDFbszSOj6Dir9uc37JQrIt5Fcz5ExGNYLJaGwrFoZ47JaUSkrah8iIhHOVE+Fu/OodblURfjiUgrUfkQEY8yIjmcYD87+aVVbDlSZHYcEWkDKh8i4lF8bFYmpEYDsGiXTr2IeCOVDxHxOJN615WPTzZnsGpfPq7vnX5xuQxe+HI3v/54O5U1tWZFFJFzoKtdRMTjTEiNJjzQl8ziCm58Yy0pUYH84cqBDO3SiRcX7eVvy+rWA6muMfj9jP5YLBaTE4tIS2h5dRHxSIfyy/jn12l8tCmDksoafG1WrhqWwDtr0wGwWMAw4OmpffjJ2GST04pIS76/VT5ExKOVVFTziw+28tW3383/uG1cMrEhfjz7+S5sVgtv/XQEY7tHmphSRLS3i4h4jWA/H2bfMJT7JvbAYoELe0Xz+JRe3H5eMlcM7kyty2DmvE2kH9N+MCLthUY+RKTdKCqvItTfp2GOR0V1Lde8upqtR4tJjQnmg7tHE+BjY8PhQuauOcyGQ4U8f+UAzu8ZZXJyEe+n0y4i0mFkF1cw7a+ryC2pbPLxEcnh/PvO0W5OJdLx6LSLiHQYsaF+zLlpKJFBjob7AnxtzBiSAMD6QwVkFR83K56INEGX2opIuzckqRPfPHYh5VU1AAT42vG1WzlSUM66QwUs2JrFz8anmJxSRE7QyIeIeAVfu5WwAF/CAnzxtdf91TZ1UDwAn27NNDOaiPyAyoeIeK1L+8dht1rYnlHMwbxSs+OISD2VDxHxWuGBvozrUbf+h0Y/RDyHyoeIeLVpA+tOvby9+rB2yRXxECofIuLVpvSLo1dsMMfKqrj61dV8tOmo2ZFEOjyVDxHxav6+Nv5z9xgu6hNDVY2Lh/69lfmbVUBEzKTyISJeL8hh59Ubh3Lr2K4APPrhdrbqFIyIaVQ+RKRDsFotPHFpHyb2iqaqxsUdb29o8X4wLpfBkt057MspaaOUIh2DllcXkQ6lpKKaH//tG/bnlmK1wIW9Yrh7QjeGdul02ueVVdbw0L+3NOyuO6JrOHdNSOHCXjHuiC3i8bS8uojIKQT7+fCPW4YzplsELgMW7crhutfWsCm9EKjbrO7DjUc5UvDdqMiRgnJmzP6Gr77NwcdmwWa1sO5QAT99c4NO34icBY18iEiHtT+3lN8s2MmKvXlEBTv4+83DePKTHWw9WkyIn52/Xj8EX7uVn7+ziYKyKqKCHcy5cSidw/z55fztLNmdy6Te0fz9luFm/yoiptOutiIizVRaWcOMv33DnibmcVgtYLVYqHEZ9O8cyms3DyUu1B+Ag3mlTPrTclwG/Pe+cfSND3V3dBGPotMuIiLNFOSw8/rNwwgL8AEgNSaYRQ+dz1VDE3AZUOMyuHxQPB/cNbqheACkRAVx2YC6Bcz+umS/KdlF2iuNfIiIAHtzSlixN49rRyQR5LBjGAYLtmVRXevix4M7Y7FYmnzO5BdXAPDVA+NJjQ12d2wRj6GRDxGRFuoZE8zt56UQ5LADYLFYmDowniuGJDRZPE48Z0q/WACe/XwXHvbfciIeS+VDROQcPHJxKr42Kyv25vH59myz44i0CyofIiLnICUqiLsmdAPgNwu+pbSyxuREIp5P5UNE5Bz9fEI3ukQEkOOs5Nn/7sLl0ukXkdNR+RAROUd+PjZ+c3k/AN5dl87MeZsor9IIiMipqHyIiLSC83tG8cKMAfjYLHyxI5srZ6/WKRiRU7CbHUBExFtcPTyRlKhA7nx7IzuznPzrm0PMvKB7o2PeXn2It1YfxmUY2K1W7p3YvWG9EJGOQiMfIiKtaFjXcH59WW8A/r7yIGXfG/3Yl1PCM5/tZF9uKQfyytiTU8LjH24nt6TCrLgiplD5EBFpZVMHxNM1IoDC8mreWXsYAMMweOKTHdS4DM7vGcX7d4xiQEIoJZU1PFc/SfXPi/ZyyZ9XsHRPrsm/gUjbUvkQEWlldpuVn9efbnltRRrHq2r5eEsGaw4W4Odj5XfT+zEyJYJnp/fHYoGPt2Ry5Zxv+POifezOLuGnb67n1eUHtGiZeC2VDxGRNvDjwZ1J6ORPfmkl/Z7+igff3wrAvRf2IDE8AID+CaHcNKoLAJvSi/C1W7mwVzSGAbO+2M0LX+0xLb9IW1L5EBFpAz42K7+YnApAbf26H33jQ/jZeSmNjnt4cirdogLpHObPv+8czRu3DOOpqX2AujkjmUXH3RtcxA20sZyISBsqKq+iotoFQGSQL3bbyf/NV13rwm61NNpD5ppXV7M2rYBbRnfhmfo1REQ8makbyz399NNYLJZGt169erX224iItAthAb7EhvoRG+rXZPGAulGSH25ed9/EHgC8u/4Iuc4KDMPQwmXiNdpknY++ffuyaNGi797EruVERERaYky3CIZ26cTGw4U8/MFW8koq2ZNTwotXD2L64M5mxxM5J20y58NutxMbG9twi4yMbIu3ERHxWhaLhXsvrLtiZuW+fHZnl2AY8PRn31JQVmVyOpFz0yblY9++fcTHx5OSksINN9xAenr6KY+trKzE6XQ2uomISN2S7VcM6Uy/ziE8eVkfesUGU1Reze+/2GV2tHNmGAY5Ti2u1lG1+oTTL774gtLSUlJTU8nKyuKZZ54hIyODHTt2EBwcfNLxTz/9NM8888xJ92vCqYhIYxsPFzBj9moA/nPXaIZ1DTc50dmpdRncNXcjC3fm8OtLe3P7D64Aqql1sS2jmIrqWgBSY4KJCHIAdZNztx0torKmbhJv37hQQgN83PsLSJNaMuG0za92KSoqokuXLvzpT3/itttuO+nxyspKKisrG352Op0kJiaqfIiINOGxD7fx3voj9I0PYcG9406aqNoePP/lbmYvOwCA1QL/+MlwJqRGA5BfWsndczey/lBhw/E+NgtT+sWRGO7PBxuOklvy3XdGqL8Pr1w/hHE9dHrfbKZe7fJDYWFh9OzZk/379zf5uMPhICQkpNFNRESa9n+X9CLQ18a3mU4W72p/y7B/siWjoXgMTgrDZcC9725m/uajzN98lGkvr2L9oUICfG30jAkiKTyA6lqDT7dm8srSA+SWVNIpwIeeMUHEhDgoPl7NLf9cxz+/TtOKsO1Im1+GUlpayoEDB7jpppva+q1ERLxeeKAvN43uypzlB3h5yT4m9o4+7eiHy2WwNq2A/gmhBDnq/sovKq9ib04pw7t2cuvISX5pJY9/tB2AO89P4aGLenL962vZeLiwYQVYgJTIQF67eRjdo4MA2JFRzLx16eSXVDJtUDyT+8Tia7dSUV3LL+dv56NNGTzz2U52Z5Xwm+l9cdhtbvud5Oy0+mmXX/ziF0ydOpUuXbqQmZnJU089xZYtW9i5cydRUVFnfL4WGRMROb380krGPb+EimoX//rpCM7v2fTfrYZh8Mv5O3h3XTpdIwL4+y3DKK2s5c63N5DjrOSeC7rzi4tT3ZZ71he7eHX5QQYkhDL/52OxWS3klVTy3Oe7OFpYDkD36CAeu6R3s+dxGIbBG6vSeO7zXbgMGNqlE3NuHEpUsKMtfxVpgqlzPq699lpWrFjBsWPHiIqKYty4cTz77LN069atWc9X+RARObPfLtjJG6vS6BMXwtXDEvCxWzm/ZxQJnQIajnl79SGe+OTbhp+DHHaqal1U1U/WBHjl+iFcOiCuzfMWllUx9vkllFfV8vebhzGpT0yrvv7yvXncM28TJRU1xIX68frNw+jXObRV30NOz6MmnLaUyoeIyJnlOisY98LSRkXCaoELUqMZ1jWc49W1/G3pfmpcBjMv6MaGQ4WsTSsAYFLvaOLD/Hlr9WH8fWy88ZNhjE6JOKdTMCUV1Xy5I5vzU6OIDvYD4GBeKesPFXBRn1j++XUaLy/ZT5+4EP57X9tMlD2YV8rtb23gYF4Zfj5W/nDlQKYOjG/195GmqXyIiHQAC7Zl8uWObABySypZV18uvm/6oHhevGYQNS6D11YcJMDXxi2ju+IyDG59cz0r9+UD0K9zCBekRmOzWgjx8+HKYQmE+DXv1MfBvFJ+9tYGDuSVER3s4NWbhpJbUsmD72+hvKoWX5sViwUqa1zMuXEIl/Rru5EWZ0U19727mWV78gCYeUE3Hr4oFau1/V0V1N6ofIiIdEAH80r5z8bvLkXtEh7Az8an4OfT9ATM4uPV/HbBTj7dmtloBAUgJSqQv988jJSooNO+59f787l77kacFd/tO+Nrs1JVW/d6EYG+HKtfkTU1Jpgv7j+vzYtArcvgha928+ryg0DdSM+L1wwiuJllSs6OyoeIiDRbQVkV8zdncDCvFIAlu3PJKq4g2M/OX68fcsoJrbuynFzxt284Xl3L0C6d+ONVA3n2810s3JkDwE/GdOVXl/ZmV1bdZcFTB8bRPfrkxSbbyvzNR3n0w+1U1bjoER3E6zcPo2tkoNvev6NR+RARkbOWV1LJXXM3svFwIVYL/PJHvbl1bDIr9+Wx4VAhI1PC6R0XwvRXvuZo4XHGdY/kjZ8Mw2G34XIZfLDxCGEBvlzcN9bsX4WtR4q4o/7qnlB/H37UPxaw0L9zKNePTDI7nldR+RARkXNSWVPLkx9/y/sbjgDQKcCHwvLqhscddiuVNS66RATwycyxhAX4mhX1jHKcFdzx9ka2HilqdP//HhxPzxj3jcR4O49a4VRERNofh93G72f05+mpfbBZLRSWVxPiZ2dKv1iCHXYqa1wE+tp4/eZhHl08AGJC/Hj/jlE89+P+PHxRT3rH1X0xnjg9JO7X5iuciohI+2SxWPjJ2GQGJ3UivaCcSb1j8Pe1UVZZw+LdufSIDmo3Iwd+PraG0yzhQb78av4OFu3KYeYF3U1O1jFp5ENERE5rYGIYUwfG4+9bd9VMoMPOtIHxDSMI7c3EXnULnG05UkRuSYXJaTomlQ8REelQYkP9GJAQimHA0t3tb3M+b6DyISIiHc5FvetGPzTvwxwqHyIi0uGc2Ftm5b58jlfVmpym49GEUxER6XB6xQbTOcyfjKLjnP+HpTh8rFw7PEkTUN1EIx8iItLhWCwWrhjSGajbF+dIwXH+8NUe5q1NNzlZx6CRDxER6ZAenNSTKf3iqKp1sXBnNq8sPcBTn+6gR0wQw7uGmx3Pq6l8iIhIh2S1WugTX3e58MCEUA4dK+e/27K4+tXV2K2WhoXWLhsQb3JS76PTLiIi0uFZLBb+cOUARnQNxzCgutagtLKGJz7eQWH9rrzSelQ+REREgABfO+/fOYq1v5zIN49dSGpMMIXl1bzw1W6zo3kdlQ8REZF6FouFmBA/4sP8+d2P+wHw7rojbEovNDmZd1H5EBERacLwruFcOTQBgJ++uZ5ZX+ziSEG5yam8g8qHiIjIKTw+pRcpUYEUlVfz6vKDTPzjclbtyzc7Vrun8iEiInIKEUEO/vfAeF67aSjDu3aiqtbFzHmbOHyszOxo7ZrKh4iIyGnYbVYm943l7dtGMigxjOLj1fzsrQ04K6rNjtZuqXyIiIg0g5+PjVdvGkp0sIO9OaWc9/xSfrtgJ+nHNA+kpVQ+REREmikmxI83bhlOYrg/xcereWNVGhf/eQVfbM8yO1q7ovIhIiLSAv0TQln2iwv450+GM6JrOMera7n7nU388X97qKhu2Q65BWVVLN+bR3F5xzqFYzEMwzA7xPc5nU5CQ0MpLi4mJCTE7DgiIiKnVFPrYtYXu3ljVRoAnQJ8uGpYItePSKJrZCAARwvL+ebAMVyu775uDWBdWgH/3Z5FVY0LPx8r0wbGc8PILgxMDDPhNzl3Lfn+VvkQERE5R/M3H+X/fbWXjKLjDfed1yMSH5uVpXtyOd03bUSgL8e+t4R7/86h3DgqiakD4wnwbd0t2Mqravh0Sya+ditXDElo1ddW+RAREXGzWpfBsj25zF1zmGV78xoVjmFdOhEW4Nvo+JgQB1cPS2RAQigbDxfyztp0/rsti6paFwDBfnZmDEnghpFJ9IgJPuX7fnMgH+fxGi7sFY2v3UpBWRXL9uTSJSKQIUlhWCwW9uWU8M7adD7cdJSSihoSOvmz/JELsFktrfb7q3yIiIiY6EhBOR9uOkqty+CKIQkk15+COZOCsio+2HCEeevSOfy9q2hGJodzw6guXNI3Fl/7d9M1l+zO4bZ/bcAwIDLIl0GJnVixN6+hwPSKDSbE34d1aQUNz+kSEcANI5O4ZUxXHHZbK/3GKh8iIiLtmstlsGp/PnPXHGbRrhxOTBeJDPLl6mGJXDciicqaWn78yjeUVNbg72Pj+Pcmu/aIDiK9oJzKmroSYrNamNQ7mhtGdmFc90isrTjicYLKh4iIiJfIKj7Ou+uO8N66dHJLKgGwWCDYYcdZUcOI5HD+desIlu/NZVdWCRf2imZgYhjF5dV8sjWD41W1TBsUT1yof5vmVPkQERHxMtW1LhbvymHumnRW7a/bX6ZzmD+f3DOWyCCHyela9v3dutNoRUREpE342Kxc0i+OS/rFkZZfxqKdOVzcN9YjikdLqXyIiIi0M8mRgfxsfIrZMc6aVjgVERERt1L5EBEREbdS+RARERG3UvkQERERt1L5EBEREbdS+RARERG3UvkQERERt1L5EBEREbdS+RARERG3UvkQERERt1L5EBEREbdS+RARERG3UvkQERERt/K4XW0NwwDA6XSanERERESa68T39onv8dPxuPJRUlICQGJioslJREREpKVKSkoIDQ097TEWozkVxY1cLheZmZkEBwdjsVjOeLzT6SQxMZEjR44QEhLihoTtjz6j5tHn1Dz6nJpHn1Pz6HM6s/byGRmGQUlJCfHx8Vitp5/V4XEjH1arlYSEhBY/LyQkxKP/T/EE+oyaR59T8+hzah59Ts2jz+nM2sNndKYRjxM04VRERETcSuVDRERE3Krdlw+Hw8FTTz2Fw+EwO4rH0mfUPPqcmkefU/Poc2oefU5n5o2fkcdNOBURERHv1u5HPkRERKR9UfkQERERt1L5EBEREbdS+RARERG3apflY9asWQwfPpzg4GCio6OZPn06e/bsMTuWx/v973+PxWLhgQceMDuKx8nIyODGG28kIiICf39/+vfvz4YNG8yO5VFqa2t54oknSE5Oxt/fn27duvHb3/62Wfs4eLMVK1YwdepU4uPjsVgsfPzxx40eNwyDJ598kri4OPz9/Zk0aRL79u0zJ6xJTvcZVVdX8+ijj9K/f38CAwOJj4/n5ptvJjMz07zAJjnTn6Xvu+uuu7BYLPz5z392W77W1C7Lx/Lly5k5cyZr1qxh4cKFVFdXM3nyZMrKysyO5rHWr1/Pq6++yoABA8yO4nEKCwsZO3YsPj4+fPHFF+zcuZM//vGPdOrUyexoHuX5559n9uzZ/PWvf2XXrl08//zzvPDCC7z88stmRzNVWVkZAwcO5JVXXmny8RdeeIGXXnqJOXPmsHbtWgIDA7n44oupqKhwc1LznO4zKi8vZ9OmTTzxxBNs2rSJjz76iD179jBt2jQTkprrTH+WTpg/fz5r1qwhPj7eTcnagOEFcnNzDcBYvny52VE8UklJidGjRw9j4cKFxvnnn2/cf//9ZkfyKI8++qgxbtw4s2N4vEsvvdT46U9/2ui+K664wrjhhhtMSuR5AGP+/PkNP7tcLiM2Ntb4wx/+0HBfUVGR4XA4jHfffdeEhOb74WfUlHXr1hmAcfjwYfeE8kCn+pyOHj1qdO7c2dixY4fRpUsX48UXX3R7ttbQLkc+fqi4uBiA8PBwk5N4ppkzZ3LppZcyadIks6N4pE8//ZRhw4Zx1VVXER0dzeDBg3n99dfNjuVxxowZw+LFi9m7dy8AW7duZdWqVUyZMsXkZJ4rLS2N7OzsRv/uhYaGMnLkSFavXm1iMs9WXFyMxWIhLCzM7CgexeVycdNNN/HII4/Qt29fs+OcE4/bWK6lXC4XDzzwAGPHjqVfv35mx/E47733Hps2bWL9+vVmR/FYBw8eZPbs2Tz00EP88pe/ZP369dx33334+vpyyy23mB3PYzz22GM4nU569eqFzWajtraWZ599lhtuuMHsaB4rOzsbgJiYmEb3x8TENDwmjVVUVPDoo49y3XXXefwmau72/PPPY7fbue+++8yOcs7affmYOXMmO3bsYNWqVWZH8ThHjhzh/vvvZ+HChfj5+Zkdx2O5XC6GDRvGc889B8DgwYPZsWMHc+bMUfn4nn//+9+88847zJs3j759+7JlyxYeeOAB4uPj9TlJq6iurubqq6/GMAxmz55tdhyPsnHjRv7yl7+wadMmLBaL2XHOWbs+7XLPPfewYMECli5dSkJCgtlxPM7GjRvJzc1lyJAh2O127HY7y5cv56WXXsJut1NbW2t2RI8QFxdHnz59Gt3Xu3dv0tPTTUrkmR555BEee+wxrr32Wvr3789NN93Egw8+yKxZs8yO5rFiY2MByMnJaXR/Tk5Ow2NS50TxOHz4MAsXLtSoxw+sXLmS3NxckpKSGv4+P3z4MA8//DBdu3Y1O16LtcuRD8MwuPfee5k/fz7Lli0jOTnZ7EgeaeLEiWzfvr3Rfbfeeiu9evXi0UcfxWazmZTMs4wdO/akS7X37t1Lly5dTErkmcrLy7FaG//3is1mw+VymZTI8yUnJxMbG8vixYsZNGgQAE6nk7Vr13L33XebG86DnCge+/btY+nSpURERJgdyePcdNNNJ83bu/jii7npppu49dZbTUp19tpl+Zg5cybz5s3jk08+ITg4uOHcaWhoKP7+/ian8xzBwcEnzYMJDAwkIiJC82O+58EHH2TMmDE899xzXH311axbt47XXnuN1157zexoHmXq1Kk8++yzJCUl0bdvXzZv3syf/vQnfvrTn5odzVSlpaXs37+/4ee0tDS2bNlCeHg4SUlJPPDAA/zud7+jR48eJCcn88QTTxAfH8/06dPNC+1mp/uM4uLiuPLKK9m0aRMLFiygtra24e/08PBwfH19zYrtdmf6s/TDUubj40NsbCypqanujnruzL7c5mwATd7++c9/mh3N4+lS26Z99tlnRr9+/QyHw2H06tXLeO2118yO5HGcTqdx//33G0lJSYafn5+RkpJi/OpXvzIqKyvNjmaqpUuXNvn30S233GIYRt3ltk888YQRExNjOBwOY+LEicaePXvMDe1mp/uM0tLSTvl3+tKlS82O7lZn+rP0Q+35UluLYXTw5QlFRETErdr1hFMRERFpf1Q+RERExK1UPkRERMStVD5ERETErVQ+RERExK1UPkRERMStVD5ERETErVQ+RERExK1UPkRERMStVD5ERETErVQ+RERExK1UPkRERMSt/j9l9R9mXokgWwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(np.array(correct_df['Time'])[7:-7],df_avg)"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"id": "y6b3kv0kBtHV"
},
"outputs": [],
"source": [
"L=0.46\n",
"me = 9.11 * 10**-31\n",
"c = 299792458\n",
"m_mu = 1.883531627 * 10**-28\n",
"def gamma(u):\n",
" return 1/np.sqrt(1-(u/c)**2)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 449
},
"id": "Eg95ViryBwzG",
"outputId": "62f8ee4c-7e67-4ec0-8626-90bf84ad290d"
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjIAAAGwCAYAAACzXI8XAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABdx0lEQVR4nO3deXxTZfY/8M/N3r20hS5Q9n3HslhhZJFFdBgXZlQUKIvrgI7ijMpPHEUdGec7bqOIy0ABFXBwG3UUBlBAR9YCssi+L13Y2nRN0uT+/kjubdKmbZImuUn6eb9eeUmTNDkXKjmc5zznEURRFEFEREQUhlRKB0BERETkKyYyREREFLaYyBAREVHYYiJDREREYYuJDBEREYUtJjJEREQUtpjIEBERUdjSKB1AoNlsNly4cAFxcXEQBEHpcIiIiMgDoiiitLQUGRkZUKnqr7tEfCJz4cIFZGZmKh0GERER+eDs2bNo06ZNvY9HfCITFxcHwP4bER8fr3A0RERE5Amj0YjMzEz5c7w+EZ/ISMtJ8fHxTGSIiIjCTGNtIWz2JSIiorDFRIaIiIjCFhMZIiIiCltMZIiIiChsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkSEiIqKwxUSGiIiIwhYTGSIiIgpbTGSIiIgobDGRIUVVmq1Kh0BERGGMiQwpZvmWU+j93Fp8f6hI6VCIiChMMZEhxew4dRVWm4gdp64oHQoREYUpJjKkmJJKCwDgSrlZ4UiIiChcMZEhxRgdicxlJjJEROQjJjKkGCmRucpEhoiIfMREhhRjrOLSEhERNQ0TGVKEKIpyjwyXloiIyFeKJjKLFi1C3759ER8fj/j4eGRnZ+Pbb7+VHx8xYgQEQXC5PfjggwpGTP5SZbHBYhUB2Jt+LVabwhEREVE40ij55m3atMFf//pXdOnSBaIoYtmyZbjllluwe/du9OrVCwBw33334fnnn5e/Jzo6WqlwyY+kaozkaoUZreIMCkVDREThStFEZsKECS5f/+Uvf8GiRYuwdetWOZGJjo5GWlqax69pMplgMpnkr41Go3+CJb+S+mMkV8qZyBARkfdCpkfGarVi1apVKC8vR3Z2tnz/Rx99hJSUFPTu3Rtz585FRUVFg6+zYMECJCQkyLfMzMxAh04+qF2RuVLGPhkiIvKeohUZANi3bx+ys7NRVVWF2NhYfP755+jZsycA4O6770a7du2QkZGBvXv34sknn8Thw4fx2Wef1ft6c+fOxZw5c+SvjUYjk5kQZKyVyLDhl4iIfKF4ItOtWzfs2bMHJSUl+OSTT5CTk4NNmzahZ8+euP/+++Xn9enTB+np6bjhhhtw/PhxdOrUye3r6fV66PX6YIVPPnLXI0NEROQtxZeWdDodOnfujKysLCxYsAD9+vXDG2+84fa5Q4YMAQAcO3YsmCFSANSpyHBpiYiIfKB4IlObzWZzadZ1tmfPHgBAenp6ECOiQCiprHb5mkPxiIjIF4ouLc2dOxfjx49H27ZtUVpaihUrVmDjxo1Yu3Ytjh8/jhUrVuCmm25CcnIy9u7di8ceewzXX389+vbtq2TY5AfSriWDVoUqi42JDBER+UTRRKaoqAhTp05Ffn4+EhIS0LdvX6xduxZjxozB2bNnsX79erz++usoLy9HZmYmJk6ciHnz5ikZMvmJ1CPTPjkGhwpKcbncfRWOiIioIYomMosXL673sczMTGzatCmI0VAwST0yHVLsiQwrMkRE5IuQ65Gh5qHEKZEB2CNDRES+YSJDijBW2Zt92zsSmasVFthsopIhERFRGGIiQ4owOvXIAIDVJtaZLUNERNQYJjKkCCmRSYnVIU5vb9W6wqF4RETkJSYyFHRWm4hSk31pKSFKi6RYHQD2yRARkfeYyFDQlTqdfB0fpUVSjD2R4XRfIiLyFhMZCjqjY6pvtE4NrVqF5BhWZIiIyDdMZCjopKbeeIMWAOSKzBUOxSMiIi8xkaGgk44nSIiyJzItpKUlVmSIiMhLTGQo6OSKTJR9txKXloiIyFdMZCjopK3XUkUmKUYPgIkMERF5j4kMBV3tHhlWZIiIyFdMZCjopB6Z+Kjazb5MZIiIyDtMZCjoanpkXBOZy+VmiCLPWyIiIs8xkaGgk+bISD0yyY7JvuZqG8rNVsXiIiKi8MNEhoKupkfGvmspWqeBQWv/UbzC6b5EROQFJjIUdLXnyABAsmPn0mUOxSMiIi8wkaGgq90jA7Dhl4iIfMNEhoKudo8MwOm+RETkGyYyFFSiKMoD8eJdlpbsicxVJjJEROQFJjIUVKZqG8xWGwDXigyXloiIyBdMZCiopP4YtUpAjE4t35/EpSUiIvIBExkKKqPT1mtBEOT7eUwBERH5gokMBZW7HUsAKzJEROQbJjIUVPI5SwbXREaa7nuFc2SIiMgLTGQoqKSKTEKdiox9IB4n+xIRkTeYyFBQSTNk4qM0LvdLS0vlZiuqLDxviYiIPMNEhoKqvopMvEEDjcre/MuGXyIi8pSm8acQ+U/NriXXREYQBLSI0eFiqQlXys3ISIxSIjwiCqAzlyvw318KYBNFpUMhPxvWuSV6ZsQr8t5MZCio6tu1BNi3YEuJDBFFlsMFpfjdOz/BWFWtdCgUAC/dpmUiQ82DvGvJTSLD6b5Ekenc1QpMXbINxqpqdEuNQ6/WynzgUeC0T4lW7L2ZyFBQ1dcjA3CWDFEkulxmwtTF21FoNKFraiw+fuBaJEbrlA6LIgibfSmo5F1Lhro5dM10X86SIYoEZaZqTF+6AyculaN1YhSWzxjCJIb8jokMBZW0tOS+IuOYJcOKDFHYM1fb8OAHedh7rgQtorVYNmMw0hIMSodFEYiJDAVVQ82+SY7pvpc5FI8orNlsIh5f/TN+PHYJ0To1cqcPRudWsUqHRRFK0URm0aJF6Nu3L+Lj4xEfH4/s7Gx8++238uNVVVWYNWsWkpOTERsbi4kTJ6KwsFDBiKkpbDYRZSb70pK7igwPjiQKf6IoYv5XB/DVzxegVQt4Z3IW+mcmKh0WRTBFE5k2bdrgr3/9K/Ly8rBz506MGjUKt9xyCw4cOAAAeOyxx/DVV19h9erV2LRpEy5cuIDbb79dyZCpCUqrqiGNj6g9RwZw2rVUwUSGKFy99d0xLNtyGoIA/P13/XB915ZKh0QRTtFdSxMmTHD5+i9/+QsWLVqErVu3ok2bNli8eDFWrFiBUaNGAQByc3PRo0cPbN26Fddee60SIUcMc7UNKgHQqIOXy0r9MVFaNXSauu/L7ddE4W3FtjN4Zd0RAMCzv+6JW/q3Vjgiag5CpkfGarVi1apVKC8vR3Z2NvLy8mCxWDB69Gj5Od27d0fbtm2xZcuWel/HZDLBaDS63MhVtdWGca9vxm/e+h/EIE7YrOmPcZ8/S4lMcYUF1VZb0OIioqZbsz8f877YBwCYPbIzpg3toHBE1Fwonsjs27cPsbGx0Ov1ePDBB/H555+jZ8+eKCgogE6nQ2JiosvzU1NTUVBQUO/rLViwAAkJCfItMzMzwFcQfkoqLTh5qRy/5Bthqg5ewmBsYIYMALSI1kGwH7eEqxWWYIVFRE205fhlPLJyD2wiMGlwJh4f21XpkKgZUTyR6datG/bs2YNt27bhoYceQk5ODn755RefX2/u3LkoKSmRb2fPnvVjtJGh2lZThQnmSdMl9ZyzJFGrBCQ6khwuLxGFh/3nS3Df8p0wW20Y1ysVL97aB4L0LxKiIFB8sq9Op0Pnzp0BAFlZWdixYwfeeOMN3HnnnTCbzSguLnapyhQWFiItLa3e19Pr9dDr9YEOO6w5JzKVFisSg/S+Dc2QkSTF6HC1woLL5SYAcUGKjIh8cfpyOabl7kCZqRpDOiThjbsGQK1iEkPBpXhFpjabzQaTyYSsrCxotVps2LBBfuzw4cM4c+YMsrOzFYww/Dn3n1SaFajINJDIJHMoHlFYKCqtwpTF23GpzIQe6fF4P2cgDFq10mFRM6RoRWbu3LkYP3482rZti9LSUqxYsQIbN27E2rVrkZCQgJkzZ2LOnDlISkpCfHw8Hn74YWRnZ3PHUhPVrsgES0PHE0i4c4ko9BmrLMhZsgNnrlSgbVI0ls0YVO+SMVGgKZrIFBUVYerUqcjPz0dCQgL69u2LtWvXYsyYMQCA1157DSqVChMnToTJZMK4cePw9ttvKxlyRKi2OvfIBK/Zt6EDIyXSdF8mMkShqcpixX3LduJgvhEpsXp8MHMwWsXx6AFSjqKJzOLFixt83GAwYOHChVi4cGGQImoeqm01yUswm32lHpmGl5aYyBCFKqtNxB9W7ca2k1cQq9dg6fRBaJcco3RY1MyFXI8MBZ5zRSbUemSkpaXLTGSIQoooipj3xT6sPVAInVqF96cORO/WCUqHRcREpjlSrkem4e3XgFOPTIgfHCmKIl74+he8ueGo0qEQBcWr645g5fazUAnAPyb1R3anZKVDIgIQAtuvKfisCs+RaWz7NRD6S0tHi8qw+MeTAIBr2rXA0M4pCkdEFDhL/3cSb353DADw4q19cGPvdIUjIqrBikwz5Lz9Org9Mo5dS/UcUQCEz9LS4YJS+dcv/uegS3JIFEn+vec8nvvKPqT08TFdcfeQtgpHROSKiUwzpNTSkicVGWmOzNUKM2whnBwcKaxJZA7mG/H57vMKRkMUGJuPXMQfV/8MAMjJbofZozorHBFRXUxkmiHnXUuV5uBsv66yWGF2nOvUULNvixj7Y1abKO9yCkVSRaZjS/uOjb+vPRzUxmmiQNtzthgPfpgHi1XEhH4ZeHZCLx49QCGJiUwz5DJHpjo4H75So69KAGJ19S8t6TVqxOntj4fy8pJUkZl3cw+0aRGFAmMV/vnDCYWjIvKPY0VlmJ67HRVmK37VJQWv/K4fVDx6gEIUE5lmyGVpKUhVBKm6EmfQNvoXojQU72qIJjKVZitOX6kAAPRpnYgnbuwOAFi06TiKSquUDI2oyfJLKpGzZDuuVljQr00CFk3Ogk7DjwoKXfzpbIaUOP3ak/4YSag3/B4rKoMo2uNMidVhQt909M9MRIXZitfWcTs2ha/iCjNylmzH+eJKdEyJwZJpgxCr5+ZWCm1MZJohl0Mjg5TIyOcsNbBjSRLq030PO5aVuqbGQhAECIKAp2/uAQD4eMcZlx1NROGi0mzFzGU7caSwDKnxeiyfORjJsXqlwyJqFBOZZiiQFZn950uw8XBRnfu9qci0iA7tREbqj+mWGiffN6h9Em7slQabCCz49mBQ4xFFEd/uy8exIiZQ5BuL1YZZK3Yh7/RVxBs0WD5jCNq0iFY6LCKPMJFphlyOKPDzoZEPfJCH6Ut34FhRmcv98jlLHpyQK/XIXA7R6b5SxaVbWrzL/U+N7w6NSsDGwxfxw9GLQYnFZhPx9Bf78dBHuzBl8XZ5ZxiRp2w2EU9+uhffHSqCXqPCkmmD0C0trvFvJAoRTGSaIavzoZF+bPYVRRGFxiqIIvDdoUKXx0oqPK/I1CwtmfwWmz/JFZm0WJf726fEYEp2OwDAX4IwJM9qE/HEp3uxYtsZAEB+SRW+3Z8f0PekyPPXNYfw2a7zUKsEvH3PNRjYPknpkIi8wkSmGQrUQDxTtU1+7Y2HXSsSnpx8LUlyDMULxWbfkkoL8kvsO5O6pNb9V+sjo7og3qDBoYJSfLrrXMDisFhteOzjPfgk7xzUKgFDO9vPvVn840mIYugOEqTQ8u6m43hvs31swMsT++KGHqkKR0TkPSYyzZDLHBk/JjLlpmr51ztOXUGZ09fe9MiEcrPvUUc1JiPB4HaZrEWMDg+P6gLAPiSvwlxd5zlNZa62YfaKXfjy5wvQqgW8NWkA3rhrAHQaFfaeK0He6at+f0+KPJ/kncOCbw8BAOaO747fZrVROCIi3zCRaYYCVZEpN9W8lsUqYsvxy/LX8q4lQ+O7lqTt16E4R0besdRAD8HU69ohMykKRaUmvL/5pF/fv8pixYMf5mHtgULo1Cq8MzkL4/ukIyVWj9v6twYALPmff9+TIs+Gg4V48tO9AID7r++IB4Z3UjgiIt8xkWmGAnVopHMFBoDL7iXvlpZq5siE2jLJkYK6O5Zq02vUeNIxJO/dzcdRZPTPkLxKsxX3LtuJ7w4VwaBV4Z85A12WAqYPaw8AWLO/AGcdA/uIatt56gp+/9EuWG0ibr+mNZ5y/KwShSsmMs1QoCb7lptrJzIX5UREWlryJJFJduxaMlXbUBFi5xfVzJBpeFfHzX1qhuS9uu5Ik9+3zFSNnNzt+PHYJUTr1Fg6fTCu79rS5Tnd0+IxrHMKbCKwfMupJr8nRZ7DBaWYsXQHTNU2jOreCi9P7MujByjsMZFphpwPjayqtvmt6iFVZDq2jIFOrcL54kocv2jfhu3N9utonQYGrf1HM5T6ZERRdNp63XAiIwgCnvm1fUjev3aexaECo8/vW1JpwZTF27D95BXE6TX4YOZgXNsx2e1zZziqMqu2n61TIaPm7dzVCkxdsg3GqmpktWuBhXdfA62aHwEU/vhTHIHKTdU4d7UC54sr3SYpzhUZq02EpVbz77mrFfJNqqR4+r4AkBKrx5CO9i2c0u4lb7ZfA0ByCO5culRmxtUKCwQB6NwqttHnZ7VLwk197EPyXvrmkE/vebXcjHv+uRW7zxQjIUqLj+4bgqx29W+PHdG1FTqmxKDUVI1Pdp716T0p8lwuM2Hq4u0oNJrQNTUWi3MGIkqnVjosIr9gIhNhCkqqMPgv6zHs5e8x9K/f4bGP99R5jvOuJaCm4bfKYsXIv2/EsJe/l29ZL6xD3ukrHr23lMjE6jUY7lj22HTkImw2EaUmz48oAIAWMfaE53ATKhn+Js2PaZ8cA4PWsw+BJ2/sDq1awOYjF7HpiHdD8i6VmTDp/a3Yf96I5BgdVt53Lfq2SWzwe1QqAdOHtgcA5P50KuCzbCj0lZmqMX3pDpy4VI7WiVFYPmMIEh3Ts4kiAROZCHOksBTlTn0lu84U13lO7Q83kyOROXe1Qp6RoteooBLs1Zvdbl7DnTLHrqUYvQYjurUCAGw7cQUXy0yQCkOeLC0B9soCAMz/6hf8ciE0kplDBTVnLHmqXXIMpma3BwAs+MbzIXmFxirc+e4WHCooRas4PVbdfy16ZsQ3/o0Abr+mDeINGpy+XIHvDtU9LoKaD3O1DQ9+kIe950rQIlqLZTMGIy3BoHRYRH7FRCbCWBw7krRqewOfu11Jzj0yQE1FRlpGykyKwuEXx+OeIe1c7m9MTUVGjU4tY9CmRRTMVhvW7C8AYE+OPK1k/GF0FwztnIwKsxX3Ltvht50/TeHJjiV3Hh7VWR6S90le48s954srcce7W3D8YjkyEgz4+IFst8P36hOj12DSkLYAgCU/cit2c2WziXh89c9yg3ju9MEeLYkShRsmMhFG6neJ1duXcExuzt6pb2mp9tC6xGj7f4srvEtkYnQaCIIgLy99+fMFl9f1hFatwtt3Z6FjyxhcKKnCfct3+v2AS295MkPGncRoHR65wTEk779HXAYH1nb6cjnueGcLTl+uQGZSFD5+IBsdUmK8jjUnuz3UKgFbTlzGgQslXn8/hTdRFDH/qwP4yjE08Z3JWeifmah0WEQBwUQmwkgVmTjHEo77ikytRMbsPpGR/utpRUbaJRPjSKKk5SVp0qwnW6+dJURrkTttEFpEa/HzuRI8/q+fYVOo58NmE+Wpvt5WZABgSnY7tE2KxsVSkzwSvrbjF8tw57tbcb64Eh1SYvCvB7KRmeTbCcQZiVEY3zsNAJD7v1M+vQaFr7e+O4ZlW05DEIC//65fna36RJGEiUyEqUlkaioytXcuOQ/EA4AqxwnYtXcWSf8t9nppyf7e13VKhs5pe6c3FRlJu+QYvDM5C1q1gP/sy8dr65s+k8UX54srUW62QqsW0N6HConzkLz3Np9AYa2lssMFpbjz3a0oMFahS6tYfHz/tUhPiGpSzDOGdQAAfLnnAopKlV+ao+BYse0MXnHMLnr21z1xi2PiM1GkYiITYaRERkomgLrLS7UrMlXy0pI9EUmI0jn+621FpqbZV/rvoA4t5Mc9OZ7AnSEdk/HSbX0AAG9+dwyf7w7cYYz1kXYsdWoZ6/PsjZv6pOGatomotFjxyn8Py/fvP1+Cu97bgktlJvRMj8eq+69Fq/imN2Re07YFBrRNhNlqw0dbzzT59Sj0rdmfj3lf7AMAzB7ZGdOGdlA4IqLAYyITYaQemTinpMFkqZXIeNgjI/3X6GVFJkZf09Ar7T5yfj1f/G5gJh4aYT8P5slP9mHnKc+2hPuLpxN9GyIIAp6+uScAYHXeORzMN2L3mau4+/2tuFphQb82CVhx3xAkx+r9EjMAzHB8kH249bTiPUYUWD8dv4RHVu6BTQQmDc7E42O7Kh0SUVAwkYkwUkXGoFVDmjxuqnb9AKuvR6a40j58rqbZ116ZKa7wbCiddESBczVoRLeatXlve2Rq+9PYbhjXKxVmqw33f5CHM5eDd57QEQ8n+jYmq10L3NwnHaIIPP6vnzFl8XZ50uoH9/p/vsf43mnISDDgcrlZbrqmyLP/fAnuX54Hs9WGcb1S8eKtfSAIPHqAmgcmMmHs1KVy3LtsJ3aduSrfJyUyOnXNVueq2hWZWtuvqxyJjrG+ikxVtUfHGNRu9gXsE3AzHHMrmlKRAezD3l67sz96t47HlXIzZizbUSdJC5TDhfajFppSkZFIQ/J+yTeizFSN7I7JWD5jsMczdryhUasw9br2AOxbsUPtEE5PVVmsePKTvXjuywMheSq6kk5fLse03B0oM1VjSIckvHHXAKh5fhI1I0xkwtiMZTuw/mAhchZvl++Tlpa0ahX0Gvsfb+0P+9pD2erbtSRtv7baRI/O7and7AvYl1NuGWBvNuye5tlAt4ZE6zT459RBSIzW4lhRmbwjKpCqrTYcL7InMr7sWKqtbXI07v1VRwDA9V1bInf6IJfkz98mDWqLKK0ahwpKseX45YC9T6BYbSIeXbUHH+88i6U/ncINr27CJ3nnwjYp86ei0ipMWbwdl8pM6JEej/dzBno8q4koUjCRCWMnLpYDgDz+H3AaiKcR6q3I1G4IrqqnR8agVUPnSIY8mSVTXqvZV/L4mK745pFf4aY+aZ5eWoPSEgzo5xjVH4zlpVOXK2C22hClVaNNi6btJJI8Ma4b1jz6K+ROGxTwD56EaC1+m9UGALDkf+E3IO+lbw5izYEC6NQqdG4ViyvlZvxx9c+4672tOFZUqnR4ijFWWZCzZAfOXKlA26RoLJsxKCBVPaJQx0QmAkiVF6AmSdGonJaW6qnISE259TX7Ov+6sZ1LoijKPTLOzb6AfXmjZ0a8X9fs2yfb56ucCkIic6Sw5mgClZ9K9oIgoHtafNCWAKY5zl/acKgIJy+VB+U9/WHp/05isWM68f/9ri++/cOv8NT47jBoVdh28grGv/ED/r72cLNrZK6yWHHfsp04mG9ESqweH8wcjFZxPHqAmicmMmHKebmodWJNlUDakaTTOC0t1emRcZ3+K8+RcZPIJHqYyFSYrfJ5SrEBXCaRtEu2z3I5fTnwH8qHC5q+Y0lpnVrGYlT3VhBFe3IQDv57oADzv/4FAPCncd1wS//W0KpVeHB4J6x7bDhu6N4KFquIt74/hrGvbfb6UM5wZbWJ+MOq3dh28gpi9RosnT5I/v+BqDliIhOmTl2qqUTEOm21NjudtaSXl5Zq7VqSt2jbk5RKixWmaquc0PhSkZGqMYIARAVhjb59SvArMk3dsaQ0aSv26rxzHs8GUsqes8V4ZNVuiI6txL93bL2XZCZF4585A/HO5CykxRtw5koFcpZsx+wVu0LiXK5AEUUR877Yh7UHCqFTq/D+1IHo3TpB6bCIFMVEJkwddeoNqHA67brm0EjnZl/3FRlp1kyV2Sp/sAmC6wwajxMZqT/Gcc5SoDlXZALd9OmPGTKhYGjnZHRLjUOF2YqPd4TugLyzVypw77IdqLLYMLxrS7xwS2+3P1OCIODG3mlY//hwzBzWASoB+HpvPm54ZROWbznl8Unj4eTVdUewcvtZqATgH5P6I7tTstIhESlO0URmwYIFGDRoEOLi4tCqVSvceuutOHz4sMtzRowYAUEQXG4PPvigQhGHjmOOXTRAza4jALBU1+xaMtRbkXFt9q20WOWt13F6jUsfSIKHB0e6G4YXSG1aREEl2JO4i2WmgL1PlcWKU46eknCvyAiCgBnD2gMAlv10us5RFaGguMKMnNztuFRmRs/0eCy85xpoGpmkHKvX4Jlf98SXs4ehX2YiSk3V+PO/D+D2t/+H/ecj58DMpf87iTe/OwYAePHWPrixd7rCERGFBkUTmU2bNmHWrFnYunUr1q1bB4vFgrFjx6K83LXv4b777kN+fr58+9vf/qZQxKHDOZGpMDvtWrLVLC0Z6qnIWOv0yNRUZGoPZPO0IuNuhkwg6TVqZDh6g04HcHnp+MUy2ET770OrOP9N3FXKLf1bIylGh/PFlVh7oFDpcFyYqq24/4M8nLhYjvQEA3KnD/Kq36p36wR89tB1eOHW3ogzaPDzuRL85q0fMf+rAx6NDwhl/95zHs99Ze8XenxMV9w9pK3CERGFDkUTmTVr1mDatGno1asX+vXrh6VLl+LMmTPIy8tzeV50dDTS0tLkW3x80+eRKOnUpXKs2V/QpCUR10TGeWnJaY5MPRUZKdlx7pFx1+gLAImOc5caX1qqO0Mm0NpJO5cCuAvniNOJ15EwKdWgVWOy40MwlLZi22winvhkL7afvII4vQa50wch1YfzptQqAVOubYcNjw/Hb/plwCbaT/8e/comfLsvPyxnz2w+chF/XP0zACAnux1mj+qscEREoSWkemRKSuxl4KSkJJf7P/roI6SkpKB3796YO3cuKirq/xe4yWSC0Wh0uYWaEX/fiAc/zMP3h4t8fo1TTrt1TNU2ucpicVRfNGqVXJGps/3akexITcKVFlu9iUxClP05JZUNT1OVKzK6YCYyUp9M4CoyhwscE33TYgP2HsE2ObsdtGoBeaevYs/ZYqXDAQC8su4w/r3nAjQqAYsmZzV5eGKrOAP+MWkAls8YjHbJ0SgwVuGhj3ZhxtIdOHsleEdbNNWes8V48MM8WKwiJvTLwLMTekVEQk3kTyGTyNhsNjz66KMYOnQoevfuLd9/991348MPP8T333+PuXPn4oMPPsDkyZPrfZ0FCxYgISFBvmVmZgYjfJ/sPlPs0/dZbWKdIXfS8pJ0/IBOLUCvdb/92iI1++qdmn0r6klkor1s9g1iRUaaJXM6gB9MzhWZSNEqzoAJ/TIA2I8tUNrK7Wew8PvjAICXbu+DYV1S/Pba13dtibWPXo9HRnWGVi3g+8MXMea1TXh74zG5MT5UHSsqw/Tc7agwW/GrLil45Xf9/DbHiCiShEwiM2vWLOzfvx+rVq1yuf/+++/HuHHj0KdPH9xzzz1Yvnw5Pv/8cxw/ftzt68ydOxclJSXy7ezZs8EIP6jcDf+SGn7NTktLBk3DA/GkikxVtRUllfZEqPbBjt4vLQVvPHowZslEwgwZd6St2N/sy0d+SaVicWw8XIR5X+wHADwyqjPuGOj/f3gYtGrMGdsN3/7helzbMQlVFhv+tuYwbv7HD9gR5FPUPZVfUomcJdvlU9EXTc6Sp2wTkauQ+D9j9uzZ+Prrr/H999+jTZs2DT53yJAhAIBjx465fVyv1yM+Pt7lFmkqnRKZaJ09cZD6ZJyXluqtyNTetWSuv0dGSmwa27UU7GZfAGjvSGROXgrMFuzSKgvOF9s/5CMtkendOgFDOiSh2iZi+ZbTisRw4EIJZn20C1abiNsHtMZjY7oG9P06t4rFyvuuxat39ENyjA5HCsvwu3e24MlP9obUQZTFFWbkLNmO88WV6JgSgyXTvGt6JmpuFE1kRFHE7Nmz8fnnn+O7775Dhw4dGv2ePXv2AADS08N/66GvRWKp+mLQqhDt6EmREhnnpSWpIlPfoZFxhprt18WOHpi6PTKeLi0Fv9m3bZJ9aam0qtqjs6C8ddTRUN0qTo8WMbpGnh1+Zgyz//+2YtsZl51vwXChuBIzlu5AudmK7I7J+OvEvkHp/RAEAbdf0wYbHh+OSYPt1Z+Pd57FqFc2YvXOs4o3A1earZi5bCeOFJYhNV6P5TMHIzk2/HfLEQWSoonMrFmz8OGHH2LFihWIi4tDQUEBCgoKUFlp/1fw8ePH8cILLyAvLw+nTp3Cl19+ialTp+L6669H3759lQxdUVJiYtCq5YpMpcX+QeSytFTPoZF1BuI5zZGps2vJ0SNTWlXd4ICxmnOWgpfIROnUSHPsbDkVgOWlIwWRMdG3PqN7pKJtUjRKKi34bNf5oL2vscqCGUt3oNBoQpdWsXhnSvCXTRKjdVhwe198+lA2uqXG4WqFBX/6ZK+iB1FarDbMWrELeaevIt6gwfIZQ9CmRbQisRCFE0UTmUWLFqGkpAQjRoxAenq6fPv4448BADqdDuvXr8fYsWPRvXt3PP7445g4cSK++uorJcNWXKXZnphEOSUytZeWtM5LS9X1DcSzJykWq4gr5Q1XZADIyY47ZQo0+wI1W7ADsXMpUib61ketEjDtuvYA7FuxbUGYhGux2jDro104VFCKlnF65E4fVOdnLpiy2iXh60eGYe747ojSqhU7iNJmE/Hkp3vx3aEi6DUqLJk2KGITaCJ/U3ThtbEybmZmJjZt2hSkaBTgYyldat6N0qoRVTuRkU6/dlpacq7I2GwipM8r5zOaCo326bi1P1S0ahVidGqUO/po6ltiUaLZF7D3yWw7eSUwFZkI3LFU2x2DMvHauiM4cbEcm45exMhurQL2XqIo4unP9+GHo5cQpVVjSc6gkKg4aNUqPDC8E27um45n/30AGw4V4a3vj+HLny/ghVt7Y3jXlgGP4a9rDuGzXeehVgl4+55rMLB9UuPfREQAQqTZl7wj9cjonZeW5B4Zx+nX9VRkqp3+1R3rNPOlqNR+0J67fx1L9xU3WJEJ/tISALQNZEVGniETuYlMrF6DOwbZe0UCvRX7re+O4V87z0ElAG/dPQB92oTWYYdtWtQcRJmeUHMQ5awVu1AYwIMo3910HO9tPgEAeHliX9zQIzVg70UUiZjIhCFp11KUVoUorWuzr9l5aclNRca5z0WjFmBwJDvSRGCpJ8ZZvAcNv+UKJTLSziV/V2Qul5lwyXGGU5dWkTMMz51p17WHSgB+OHpJrkL52+e7z+GVdUcAAPNv6R2yH9bSQZTr5gzHvcM6QK0S8J+9+Rj9yiYs+8n/B1F+kncOC749BACYO747fpvV8K5NIqqLiYyCfN2jIa3dR+mce2TsiYTz0pJUkXFe65eOJ5CeE6V1XQqqPUcGqEluPElkgr1NVOqROePnisyRQns1JjMpKujJWbBlJkVjbM80AEBuAI4t+On4JTzxyV4AwAPXd8SUa9v5/T38LVavwbxf98SXs4fKB1E+++UB3ObHgyg3HCzEk5/af1/uv74jHhjeyS+vS9TcMJEJQ1JiYtA0vLRUs/3aqSJjdarIqFQuiYwg1Ez7dSZvwa6of9aG3OwbxCMKgJpE5nK5GcYq/23BrumPibw5RO7M/JV9K/Znu87Ljd/+cLSwFA98YB+xf3OfdDx5Y3e/vXYw9MpwPYhyr+Mgyue+PIDSJvy87Tx1Bb+XZuhc0xpPhdnvC1EoYSIThqSlIoPOqdnXUnfXkqGBiowg2HetGJwSmXiD1u0IdE9myShVkYkzaJESa29A9mdVRtqx1C2CzlhqyMB2LdCndQJM1Tas2OafAXlFpVWYlrsDpVXVyGrXAq/cEZ4j9t0dRLn0p1MY/eomfOPDQZSHC0oxY+kOmKptGNW9FV6e2Dcsf1+IQgUTmTBU2UBFxuy8tOSuIuOo2Ggcf3E6JzL1bYNNjLYnCvUNnbPaRDmmmCDvWgJqjirwZ5/MkQg9mqA+giBgxrD2AIDlW07LvVa+KjdVY+bSnThfXIkOKTF4f+pAl5+1cCQdRPnBzMFonxyNQqMJv/fyIMpzVyswdck2GB3J3cK7r4FWzb+GiZqC/wcpyNdBplLSEqVznuwrHRrptLTkpiJT7VhaUjsSGamiA9SfyDRWkSl3mgqrRD+Jv2fJiKLoVJFpHokMANzcJwOt4vQoKjXhP/su+Pw61VYbHlm5G/vOlyApRofcaYOQFEGTkX/VpSXWPHo9HrmhC3RqlctBlA0lgJfLTJi6eDsKjSZ0TY3F4pyBLv//EZFvmMiEIZc5MtqaOTJWmyhXXOwD8epWZKRER6tSya8h8TmRcSwrqVUC9AocbCfvXLrkn4pMgbEKpVXV0KgEdExpHktLAKDTqDA1296Iu/jHkz6N6xdFEfO/+gUbHIPd3p86EO1TYvwdquIMWjXmjOmKb/7wK2R3TJYPovz1mz9g+8m6B1GWmaoxfekOnLhUjtaJUVg2Y7Bc6SSipmEiE4aq5LOWXJeWpB1LAKDVqGBwJBXmaps8tdXq6JFRq6WlpZofgcYSmfrmyMhbr3XqoJyXU5u/KzLSidcdUmKa3YnDdw9pB71Ghf3njdhx6qrX3//PH07ig62nIQjA63f2R1a7FgGIMnR0bhWLFfcNwWt31hxEece7W/DEJz/LB1Gaq2148IM87D1XghbRWiybMRjpCVEKR04UOZrX39IRQu6RqTXZ1zmR0dRq5JWqMtK8GI2jIuPS7NtIIlPfEQXljh1LSp3Q6+8eGWnHUiQPwqtPUowOt1/TGoD3A/K+2ZePv3xzEADw9E09ML5P+B/s6glBEHDbANeDKP+185x8EOXjq3/Gj8cuIVqnRu70wegc4XOJiIKNiYyfbTtxGQ98sBP5JZV+e80F3x7Eq45hYoDTriWtuqZHxmKV+18AaSBezR+vNN23drOvJ0tLjc2RUWoYnqS9oyJTVGryyynO0kTfSD6aoCEzhtq3Yv/3lwKPm1jzTl/Box/vAQDkZLfDzGGNn2QfaZwPouyeVnMQ5Vc/X4BWLeCdyVnon5modJhEEYeJjJ/d+d5WrD1QiD+u/rnR5woejMS7XGbCu5tO4B8bjsoJQ81kX+elpWq5IqNWCVCrBGjUKjlhkZIf54F5gHfNvvXtWlLqeAJJYrROjvGMhx+8DTkS4YdFNqZLahx+1SVF3mbcmJOXynHvsp0wV9swukcq/jyhlyJLjKEiq10Svnq45iBKQQD+/rt+uD4IZzYRNUdMZALk/FX/VGSc2y2lhKFmsq/KZWlJ3nrtNJNCqsrUV5FxXlpydzwBUJPIVFqsdU7SBmp2LSm1tATUVGVOXWpaImO1iTha1Px2LNUmVVQ+3nG2wcFvV8rNmJ67HVcrLOjbJgH/mNRf3hHXnEkHUW5+YiTWPTYct/RvrXRIRBGLiUyIszmd7SJ9oNQ32VdaWtI5zaWQEpWaiowjkVHX7ZGpryITZ9DKW8XdLS/JU30VmCEjkfpkTjexT+bslQpUWWzQa1Rom6T8ycxKub5LS3RqGYMyUzVW7zzn9jlVFivuXbYDpy5XoHViFP6ZM1Be6iS7lnF69sQQBRgTGQV5Un13Pq26pNJ1acmgUyPa6dBIadlI69QbI1VkpOTHlx4ZtUqQjy5w1/CrdI8M4FSRaeLOJWl+TJfU2GZdWVCpBEx39Mrk/nSyzmGJNpuIxz7eg11nihFv0GDZjEFoFWdQIlQiauaYyIQ45wZeo1yRcTT7amp2LdmXfeouLRlqzZKRjihQqzzffg3UTPd1V5FR6ngCZ/6qyDS3ib4NmXhNGyREaXH2SiXWHyx0eWzBtwfx7f4CaNUC3ps6EJ1b8feLiJTBRCbEVTudVl1a5ajIyJN9a5aWnB93Hnmul5eWHBWZWktLnlRknB9z1/CrdLMvALRP8c8sGXmiLxMZROnUuHtIWwCuW7GXbzmF93+wf/1/v+2HazsmKxIfERHARCbkOZf0pWUdk5vJvkBNtUTnZmlJqshU115a0jU+RwZoeLpvKFRk2ibZKzIXSipdjmTwVnOeIePO1Ox20KgEbDt5BfvPl2D9L4V47ssDAIA/ju2KWwewiZWIlMVEJsQ598hIS0uV8mRfFVQqQV4ekhIdrdp5acm1R0aq8NTetSQIkPtg3EloYJaMNBAvRsFzY1JidYjRqSGK9oP5fGGutuHERfvSFCsydukJUbjJMdju+a9+wcMrd8MmAncNysSskZ0Vjo6IiImMojxpJXXpkamshiiKLnNkAMg7RaQkQ5raC6DOCdhys6/aNZGJN2ihaqC5NdSXlgRBcOqT8S2ROXGpDNU2EXF6DdIT2LgqmeHYir391BVUWqy4vmtLvHBr72Y9K4aIQgcTmSDz9iA+1x4ZCyxWEVKRRi8nMvb/SomM866l2hWZ2kcUSAlKSmzDB9iF+tISUNMn4+vOJemMpa5pcfyQdtI/M1E+M6l7WhwW3j3ApQ+LiEhJHPoQIPWlK94eKOzSI1NVLVdjgJqKTEytioxOXXfXUs32a9elpb6tEzBnTFdc07bhw/0SGzhvKRQqMkDTdy4194m+DXl5Yh+s3nkOM4d1QJyh/l4qIqJgYyITZDanTMbbOTLGSgtMjoRErRLkXphoxyC6YrdLS67NvjUD8ezfq1IJeOSGLo3G0dAJ2NJkX6UTmabOkqk5Y4kDzGrr3CoOc2/qoXQYRER1sD4cIPXlKF4WZOrMkZGH4WlU8vJH7YqM69KSo0emzkA87/7oGzo4UunTryV+q8hwxxIRUdhgIhNktib0yBgraxIZ523T0q9LKswAXJeW6lZkXAfieSpebvY113msZmlJuV1LANDekcicu1opX6enKszV8oGT3LFERBQ+mMgEmXMe40lDqdXlrKVqeaqvtBsJqNn27G5pqW6PjOvSkqdqmn2rXe63WG0wO5IkpSsyreL0MGhVsNpErw/tPFpoX1ZKidUhOVYfiPCIiCgAmMgEmbfNvrXnyDhP9ZVEOxIIaWu0+7OW3A/E85R0RIGx0uKy80rasQRA8QMDVSpBPujxlJfLS4fZ6EtEFJaYyASZ10tLTj0yVRabPBTPeaKvVJGRHtO62bUkTQOurnVEgaekiozZanPZOSUtK+nUKpeJwkrxdZYMz1giIgpPyn/yRKh6t197+TrOPTIAcKzIvgTifNijVAmRciSd27OWpIF4dQ+W9ESMTi1/j3PDrzzVV+H+GIm0c8nbREY+Y4mNvkREYYWJTJB5W5Fx7pEBgP9bexhATaUFqJtEaNw2+zoG4vm4a0kQBLfTfUNlhozEl51L5mobfrlgBMCKDBFRuGEiE2Sid5tpXHpk9BoV9BoVYvUa+fwboG5vivPUVUOdioxvzb4AkBpvH9t/zqmRNlSm+kqknUve9Mgs33IKl8vNSInVoVdGfKBCIyKiAAiNT59mRPRycUnqaRnTMxXvTx3o9jm1KzIuS0tSs69UkbH6trQEAB1bxuCXfCNOXioDkArAvm3ZHkNo/Ci1cywtnb1SCatNbHSb+aUyE95YfxQA8MS47i6VLiIiCn2syASZzesjChpPPGpXZDTumn1rV2R8SWRS7NUO6YRoACiTe2RCI5HJSIyCVi3AbLUhv6TxLdiv/PcISk3V6N06Hr/NahOECImIyJ+YyASZ94dGNr7LKFrnWkXQNliR8W3XEgB0bGkf3X/iUk0iU7O0FBqVDLVKQGaSZw2/By6UYNWOMwCAZyf0avD0byIiCk1MZAKkvo9E7ysyjVdQPOmRMdXateTtZF/AvrQE1K7IOJaWFJ4h48yTPhlRFDH/q18gisCEfhkY1D4pWOEREZEfMZEJkPpPvxbd/ro+UkWmocSjdo+M6xwZ111L8hwZHxKZDo6lpUtlJnlmTXmI7VoCIA/Fa6gi8+3+Amw/eQUGrQpPje8erNCIiMjPFE1kFixYgEGDBiEuLg6tWrXCrbfeisOHD7s8p6qqCrNmzUJycjJiY2MxceJEFBYWKhRx0zmnLp5UZ6o9aM6tXQ1xXVpyrch4slRVnziDFi3j7OP7TzqqMqG2awlwOgX7kvuKTJXFir/85yAA4IHrO6F1YlTQYiMiIv9SNJHZtGkTZs2aha1bt2LdunWwWCwYO3YsystrPoAee+wxfPXVV1i9ejU2bdqECxcu4Pbbb1cw6qaxuVRkGn9+tQfbpRvqkZEqMlKPTLWPA/EkcsPvJftgvlBr9gWAdo4YpUMga3t/8wmcL65EeoIBDw7vFMzQiIjIzxT99FmzZo3L10uXLkWrVq2Ql5eH66+/HiUlJVi8eDFWrFiBUaNGAQByc3PRo0cPbN26Fddee22d1zSZTDCZTPLXRqMxsBfhJecqjCdbsa0eDLCrnUTo3FRkLFYRVpvodESBj4lMyxhsO3nFTUUmNJp9AdceGVEUXQ7nLCipwtsbjwMAnhrf3eXMKiIiCj8h1SNTUlICAEhKsjde5uXlwWKxYPTo0fJzunfvjrZt22LLli1uX2PBggVISEiQb5mZmYEP3AvOfTEeLS150COj16jg/LBzkuJ8JlOVxdqk7dcA0DHFvnPpuGPZpjzE5sgAQOvEKKhVAqosNhSVmlwee3nNIVRarBjYrgV+0y9DoQiJiMhffEpkOnbsiMuXL9e5v7i4GB07dvQpEJvNhkcffRRDhw5F7969AQAFBQXQ6XRITEx0eW5qaioKCgrcvs7cuXNRUlIi386ePetTPIHispzkSbOvBz0ygiC49Mm4234NABVmq89HFEhq71wKtSMKAECnUcl9L859MrvOXMXnu89DEOzbrZ0rNUREFJ58+vQ5deoUrFZrnftNJhPOnz/vUyCzZs3C/v378eOPP/r0/RK9Xg+9Xt+k1wgk0WVpqXGe9MgAQLRejVJHUuGcyKhUAqK0alRarKg0W2sG7Pm4tCTtXDp1qRw2mxiSzb6AfcLvmSsVOH25AkM6JsNms2+3BoDfZbVBnzYJCkdIRET+4NWnz5dffin/eu3atUhIqPkwsFqt2LBhA9q3b+91ELNnz8bXX3+NzZs3o02bmumqaWlpMJvNKC4udqnKFBYWIi0tzev3CQU2l6Ulz3tk1I1UUOyzZOzLKNpaSUq0zp7IVFiqawbi+ViRyUyKhkYloNJiRYGxyun069BKZNonx+CHo5fkWTKf7z6Pn88WI1avwR/HdVM4OiIi8hevPn1uvfVWAPaljJycHJfHtFot2rdvj1deecXj1xNFEQ8//DA+//xzbNy4ER06dHB5PCsrC1qtFhs2bMDEiRMBAIcPH8aZM2eQnZ3tTehBV/9APB93LTXS0+K8c0lba2t1lE4NlNuXlqwe9Nw0RKtWoW1yNE5cLMeJi+Xy0lIoNfsCNWcunb5cgXJTNV5ecwgAMHtUZ7SKMygZGhER+ZFXiYzNsSzRoUMH7NixAykpKU1681mzZmHFihX497//jbi4OLnvJSEhAVFRUUhISMDMmTMxZ84cJCUlIT4+Hg8//DCys7Pd7lgKJfUOxPPgOc6kHpnGEo/6emSAmiSnymz1qOemMR1TYuyJzKWykByIBwDtnHYuvb3xGIpKTWiXHI3pQ9srGxgREfmVT58+J0+e9MubL1q0CAAwYsQIl/tzc3Mxbdo0AMBrr70GlUqFiRMnwmQyYdy4cXj77bf98v5KEL1cWvK4IqN3rsi4PjfKkeRUmK0e99w0pGPLWOBgEQ4VlMqvV/uYBKVJQ/GOFZXhaJF95s3TN/WQt6MTEVFk8PnTZ8OGDdiwYQOKiorkSo1kyZIlHr2GJyP6DQYDFi5ciIULF/oUZ6gRvSzJWD2cxNtgRcaxBbvCYnU6osD3nffSULx950qc3j+0EoTMpGgIAmCqtv9sDuucgjE9UxWOioiI/M2nT7P58+dj7Nix2LBhAy5duoSrV6+63Kh+tkZ2LR0tLMX7m09g+8krAPzYIwOg0lxdM9m3CRUZaefS4YJSAPbpwb4ceRBIBq0a6fH2Xhi1SsAzv+7J7dZERBHIp4rMO++8g6VLl2LKlCn+jifiuexacjMRb/aK3ThcaE8QDr1wo+c9Mk49KjpN7aUlR0XG3PSBeIBjaQmA2RFbqG29lnRqFYsLJVW4Z0hbdEuLUzocIiIKAJ8+gcxmM6677jp/x9IsNDZHxvl8oHJTtceJh3NFpvaykby0ZLbWbL9uQgUlJVaHOIMGpVWh2egreWp8dwxoW4gHrvdtSCMREYU+nz7N7r33XqxYscLfsTQLDW2/rjRbUWmpGTRotto8Pq3aOZnQatzvWqr0U0VGEAS5Twaoe/p2qOiVkYA5Y7qGbKJFRERN59Pf8FVVVXjvvfewfv169O3bF1qt1uXxV1991S/BRbrau5Yul7ueC2SpFn2qyDS8a6npPTKAfXnpZ0ezb6guLRERUeTz6RNo79696N+/PwBg//79Lo+xobJhDW25vlxmdvnabLXB4mGPjEsiU3tpSarIWKweNw83xqUiE2LD8IiIqPnwKZH5/vvv/R1Hs+Gya6mxiozV5kVFxpOlpWp5+3VjRx40pkNL50SGFRkiIlJGaO2ZbQZcB+K5PlanIlPtTY9MQ0tLNc2+8tJSkysysfKvubRERERK8ekTaOTIkQ0uIX333Xc+BxTpXOfI1K7IuCYyPldkGlpaknctNS2R6ZDCigwRESnPp08gqT9GYrFYsGfPHuzfv7/OYZJUo9JshfOm69rtMlfK3fXIeHbIo7RzSK0SoKr13Citc0Wm6ZN9AXuVJyPBgAslVUxkiIhIMT59Ar322mtu73/uuedQVlbWpIAi1b5zJZjw1o/okR4v31d7aelSWe0eGRFWD5eCpLOW3D1P2rVU5pj74snreaJjS/vAuVA7+ZqIiJoPv/bITJ482eNzlpqb19cfAQAczDc63VtraalWj4zFix6ZdknRGNS+BW7t37rOY9LSkrHKIt/X1KUlAJiY1Rptk6IxrHPLJr8WERGRL/y6JrBlyxYYDAZ/vmRE82RpSeqRaWxpSaNWYfWD7qctS0tLrhWZpuewtw1og9sGtGny6xAREfnKp0Tm9ttvd/laFEXk5+dj586deOaZZ/wSWHNQZyCeY2kpVq9BmakaFqvN6bRq3ysoUkWm1OSUyPihIkNERKQ0nxKZhIQEl69VKhW6deuG559/HmPHjvVLYM2BWGumjLRrKS3BgGNFZY7t154NxGtItJsjBPzRI0NERKQ0nxKZ3Nxcf8cR8dztVneux5SbrTBV25OWdEciY2/2tT9L24RDHqN0rs24apXACcxERBQRmtQjk5eXh4MHDwIAevXqhQEDBvglqObCeWlJWlaK0qoRH2U/u8ridGhk0yoydRMZIiKiSOBTIlNUVIS77roLGzduRGJiIgCguLgYI0eOxKpVq9CyJXexeMSpJCMtKyXH6qBzVF/M1TZY/dAjo1WroFUL8kwaLisREVGk8Gm94uGHH0ZpaSkOHDiAK1eu4MqVK9i/fz+MRiMeeeQRf8cYsZyXlqSt18kxOvmIAbPVBosfemQAwKCtqcowkSEiokjhU0VmzZo1WL9+PXr06CHf17NnTyxcuJDNvg6nL1cg938n8dusNshe8B3KnHYMSWyiiGNFZfjbmkNIjbdvW0+O1UPnOPTR5YiCJu4yitapUerYft3YTBoiIqJw4VMiY7PZoNVq69yv1Wphc1QQCJj/1S/IL6lym8QA9l1LOUu243xxpXxfvEEjN/Y698g0de6LfeeSyfFarMgQEVFk8OnTcdSoUfjDH/6ACxcuyPedP38ejz32GG644Qa/BRcJCkqq6n3MJoouSQwAqFUquUemymKTt2g3NfmI4tISERFFIJ8SmbfeegtGoxHt27dHp06d0KlTJ3To0AFGoxFvvvmmv2MMa7WH3jlz94haVbPVusJsrbnfD0tLEi4tERFRpPBpaSkzMxO7du3C+vXrcejQIQBAjx49MHr0aL8GFwkayGPcZjJqlSAnMlWWmkSmyRUZHSsyREQUebz6p/l3332Hnj17wmg0QhAEjBkzBg8//DAefvhhDBo0CL169cIPP/wQqFjDkrX2EddO3FVrVIIgN/tWmp0Tmab2yDhXZJjIEBFRZPDq0/H111/Hfffdh/j4+DqPJSQk4IEHHsCrr77qt+AiQYNLS/VWZOyJRqUfKzLOxxSo/XBgJBERUSjw6hPt559/xo033ljv42PHjkVeXl6Tg4okDffIeFaREQRAxaUlIiKiOrxKZAoLC91uu5ZoNBpcvHixyUFFkgZWltw+phJqemQqLI65L35IPFx2LXFpiYiIIoRXiUzr1q2xf//+eh/fu3cv0tPTmxxUJNl0pP7Ezt3WbOddS1JFpqn9MUCtHhlWZIiIKEJ49Ql500034ZlnnkFVVd0P4MrKSjz77LP49a9/7bfgIkFDzb77zpfUuU+lqru05JeKjEsiwx4ZIiKKDF5tv543bx4+++wzdO3aFbNnz0a3bt0AAIcOHcLChQthtVrx9NNPByTQ5kItCNA5ln6kIwX0Wj9UZLi0REREEcirRCY1NRU//fQTHnroIcydOxeio5FVEASMGzcOCxcuRGpqakACbS6c58iUOo42cD7w0VfOu5a4tERERJHC64F47dq1wzfffIOrV6/i2LFjEEURXbp0QYsWLQIRX7Pj3Owr8Uci47y0xO3XREQUKXya7AsALVq0wKBBg/wZC6G+RMa/zb5aLi0REVGE4D/NQ4xaBbnZV2LQ+Lsiw0SGiIgig6KJzObNmzFhwgRkZGRAEAR88cUXLo9PmzYNgiC43BoayBcJVCpBPv1a4pyE+Io9MkREFIkUTWTKy8vRr18/LFy4sN7n3HjjjcjPz5dvK1euDGKEwacWBGg1romG3g8VGZ5+TUREkcjnHhl/GD9+PMaPH9/gc/R6PdLS0oIUkfKcdy1J/NEj4zLZlxUZIiKKECH/T/ONGzeiVatW6NatGx566CFcvny5weebTCYYjUaXWzhRCW6Wlvy8a4lzZIiIKFKEdCJz4403Yvny5diwYQNefvllbNq0CePHj4fVaq33exYsWICEhAT5lpmZGcSIm04l1N1u7Z85MpzsS0REkUfRpaXG3HXXXfKv+/Tpg759+6JTp07YuHEjbrjhBrffM3fuXMyZM0f+2mg0hlUyo1YJiNW7/rH4Y2nJeecTl5aIiChShNU/zTt27IiUlBQcO3as3ufo9XrEx8e73MKJSiXAoFXBOdfwR0VGpRLkJSo1l5aIiChChFUic+7cOVy+fDmiT9hWO7aZxzhVZfyRyAA1y0taLi0REVGEUHRpqayszKW6cvLkSezZswdJSUlISkrC/PnzMXHiRKSlpeH48eN44okn0LlzZ4wbN07BqANL5SjFxOk18qGR/kpkonRqoJwD8YiIKHIomsjs3LkTI0eOlL+WeltycnKwaNEi7N27F8uWLUNxcTEyMjIwduxYvPDCC9Dr9UqFHHBqwZ5kuFZk/FNBkSoy7JEhIqJIoWgiM2LECPkEbXfWrl0bxGhCg7Tq45zI+GP7NQBEOab7ciAeERFFCn6ihRiVoyITG4AemShHZYcVGSIiihRMZEKM1L8So69JXvy3tCRVZJjIEBFRZGAiE2Lc9sj44awlAGifHAMAaJ0Y5ZfXIyIiUlpID8RrjqRdSy5LS344/RoAnhzfDROzWqNnenjN1iEiIqoPE5kQE8iKjF6jRq+MBL+8FhERUSjg0lKIkXYtxQZg+zUREVGk4SdkiHG3aynKT0tLREREkYaJTIiRdi05z47x19ISERFRpGEiE2KkHhnnYwT8NUeGiIgo0jCRCTHSriXnWS96Df+YiIiI3OEnZIiRKjHOFRkVJ/ESERG5xUQmxEg5SwaH1hERETWKc2T8rG1SNM5cqfD5+6VdS9e0bYGnb+qB9ikx/gqNiIgo4jCRCTHOS0r3Xd9RwUiIiIhCH5eW/EyE2KTvlyoyRERE1DgmMiFGzcZeIiIijzGR8TOxaQUZJjJEREReYCITYpjHEBEReY6JjJ81tSLDHhkiIiLPMZFpAlEUUWWxAgCsNhEllRb5a19xaYmIiMhz3H7dBPf8cxt+On4ZO+eNxsAX1/vlNVmRISIi8hwrMk3w0/HLAIBv9xf47TVZkSEiIvIcE5kguL5rS4+fy0SGiIjIc0xk/KCx1MOb3IQrS0RERJ5jIhME3vS9qJnJEBEReYyJjB80lnt4U5Hh0hIREZHnmMj4gdDI4pLgRZVFxUSGiIjIY0xkgsCb5SIuLREREXmOiYwfNJZ7eJObcGmJiIjIc0xk/ODHo5cafNybZl8WZIiIiDzHRMYP/rMvv8HHvarIMJMhIiLyGBOZIGioInNz33SXr7m0RERE5DkmMkHQUG7SNina5WtvdjgRERE1d0xkgoDJCRERUWAwkQkC5jFERESBwUQmCLzZtURERESeUzSR2bx5MyZMmICMjAwIgoAvvvjC5XFRFPHnP/8Z6enpiIqKwujRo3H06FFlgm0C9u8SEREFhqKJTHl5Ofr164eFCxe6ffxvf/sb/vGPf+Cdd97Btm3bEBMTg3HjxqGqqirIkTYNdyIREREFhkbJNx8/fjzGjx/v9jFRFPH6669j3rx5uOWWWwAAy5cvR2pqKr744gvcddddwQy1SdjsS0REFBgh2yNz8uRJFBQUYPTo0fJ9CQkJGDJkCLZs2VLv95lMJhiNRpeb0hpLY1Ji9QBYuSEiIvJWyCYyBQUFAIDU1FSX+1NTU+XH3FmwYAESEhLkW2ZmZkDj9ERjzb6v39kfv81qg/emZAUpIiIiosig6NJSIMydOxdz5syRvzYajYonM40VWoZ1ScGwLinBCYaIiCiChGxFJi0tDQBQWFjocn9hYaH8mDt6vR7x8fEuN6WxR4aIiCgwQjaR6dChA9LS0rBhwwb5PqPRiG3btiE7O1vByLzHOTJERESBoejSUllZGY4dOyZ/ffLkSezZswdJSUlo27YtHn30Ubz44ovo0qULOnTogGeeeQYZGRm49dZblQvaBw3lMUxxiIiIfKdoIrNz506MHDlS/lrqbcnJycHSpUvxxBNPoLy8HPfffz+Ki4sxbNgwrFmzBgaDQamQfcLNSERERIGhaCIzYsQIiKJY7+OCIOD555/H888/H8So/I9LS0RERIERsj0ykUTFkgwREVFAMJEJAuYxREREgcFEJgicl5ZaJ0Zh2nXtlQuGiIgogjCRCQCVAJz6683y184FmX/PHornftMr+EERERFFICYyQcbGXyIiIv9hIhNk7JchIiLyHyYyQcbjCoiIiPyHiUyQ1a7IMK8hIiLyHROZAKg94s/5a1ZkiIiI/IeJDBEREYUtJjJEREQUtpjIBEDt46NSYvXyrw0a+295cowOADC6R2rQ4iIiIoo0ih4a2VwYtCr89NQoqAQBGrU9kdn4pxEoKKlCl9Q4haMjIiIKX0xkgkAUgYzEKJf74gxaxBm0CkVEREQUGbi0RERERGGLiQwRERGFLSYyREREFLaYyBAREVHYYiJDREREYYuJTBDUPrKAiIiI/IOJDBEREYUtJjJBwGMiiYiIAoOJDBEREYUtJjJEREQUtpjI+EisfTJkQ88NYBxERETNGRMZH3mRxxAREVGAMJEhIiKisMVExkcsyBARESmPiYyPvOmRISIiosBgIhMEaoGTZIiIiAKBiYyPPKnH3PerDujdOh6/6Z8R8HiIiIiaI43SAYQrT1aWnr65Z+ADISIiasZYkSEiIqKwxUTGRyL3LRERESmOiYyPuGmJiIhIeUxkiIiIKGyFdCLz3HPPQRAEl1v37t2VDouIiIhCRMjvWurVqxfWr18vf63RhEbIXFoiIiJSXmhkBQ3QaDRIS0vz+Pkmkwkmk0n+2mg0BiIsvPnd0YC8LhEREXkupJeWAODo0aPIyMhAx44dcc899+DMmTMNPn/BggVISEiQb5mZmQGJ651Nx+t97KERnQLynkRERORKEEP40KBvv/0WZWVl6NatG/Lz8zF//nycP38e+/fvR1xcnNvvcVeRyczMRElJCeLj4/0WW/un/uP2/q8fHoae6fFQqXgsARERka+MRiMSEhIa/fwO6aWl8ePHy7/u27cvhgwZgnbt2uFf//oXZs6c6fZ79Ho99Hp9sEKso3frBMXem4iIqLkJ+aUlZ4mJiejatSuOHTumdChEREQUAsIqkSkrK8Px48eRnp6udChEREQUAkI6kfnjH/+ITZs24dSpU/jpp59w2223Qa1WY9KkSUqHRkRERCEgpHtkzp07h0mTJuHy5cto2bIlhg0bhq1bt6Jly5ZKh0ZEREQhIKQTmVWrVikdAhEREYWwkF5aIiIiImoIExkiIiIKW0xk/Gjzn0YqHQIREVGzwkTGj9omRysdAhERUbPCRIaIiIjCFhMZIiIiCltMZIiIiChsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkfGTV+/op3QIREREzQ4TGT+5/Zo2SodARETU7DCRISIiorDFRIaIiIjCFhMZIiIiCltMZIiIiChsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkfGDfc+NVToEIiKiZomJjB/EGbRKh0BERNQsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkSEiIqKwxUSGiIiIwhYTGSIiIgpbTGSIiIgobDGRISIiorDFRIaIiIjCFhMZIiIiCltMZIiIiChshUUis3DhQrRv3x4GgwFDhgzB9u3blQ6JiIiIQkDIJzIff/wx5syZg2effRa7du1Cv379MG7cOBQVFSkdGhERESks5BOZV199Fffddx+mT5+Onj174p133kF0dDSWLFmidGhERESksJBOZMxmM/Ly8jB69Gj5PpVKhdGjR2PLli1uv8dkMsFoNLrciIiIKDKFdCJz6dIlWK1WpKamutyfmpqKgoICt9+zYMECJCQkyLfMzMyAxJY7bRAAIKtdi4C8PhERETUupBMZX8ydOxclJSXy7ezZswF5n5HdW+HUX2/Gpw9dF5DXJyIiosZplA6gISkpKVCr1SgsLHS5v7CwEGlpaW6/R6/XQ6/XByM8IiIiUlhIV2R0Oh2ysrKwYcMG+T6bzYYNGzYgOztbwciIiIgoFIR0RQYA5syZg5ycHAwcOBCDBw/G66+/jvLyckyfPl3p0IiIiEhhIZ/I3Hnnnbh48SL+/Oc/o6CgAP3798eaNWvqNAATERFR8yOIoigqHUQgGY1GJCQkoKSkBPHx8UqHQ0RERB7w9PM7pHtkiIiIiBrCRIaIiIjCFhMZIiIiCltMZIiIiChsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkSEiIqKwFfJHFDSVNLjYaDQqHAkRERF5SvrcbuwAgohPZEpLSwEAmZmZCkdCRERE3iotLUVCQkK9j0f8WUs2mw0XLlxAXFwcBEHw+PuMRiMyMzNx9uzZiDyjidcX3nh94SuSrw3g9YW7ULo+URRRWlqKjIwMqFT1d8JEfEVGpVKhTZs2Pn9/fHy84n+YgcTrC2+8vvAVydcG8PrCXahcX0OVGAmbfYmIiChsMZEhIiKisMVEph56vR7PPvss9Hq90qEEBK8vvPH6wlckXxvA6wt34Xh9Ed/sS0RERJGLFRkiIiIKW0xkiIiIKGwxkSEiIqKwxUSGiIiIwlazTmQWLlyI9u3bw2AwYMiQIdi+fXu9z126dCkEQXC5GQyGIEbrPW+uDwCKi4sxa9YspKenQ6/Xo2vXrvjmm2+CFK33vLm+ESNG1PnzEwQBN998cxAj9o63f36vv/46unXrhqioKGRmZuKxxx5DVVVVkKL1jjfXZrFY8Pzzz6NTp04wGAzo168f1qxZE8RovbN582ZMmDABGRkZEAQBX3zxRaPfs3HjRlxzzTXQ6/Xo3Lkzli5dGvA4feXt9eXn5+Puu+9G165doVKp8OijjwYlTl95e32fffYZxowZg5YtWyI+Ph7Z2dlYu3ZtcIL1krfX9uOPP2Lo0KFITk5GVFQUunfvjtdeey04wXqh2SYyH3/8MebMmYNnn30Wu3btQr9+/TBu3DgUFRXV+z3x8fHIz8+Xb6dPnw5ixN7x9vrMZjPGjBmDU6dO4ZNPPsHhw4fx/vvvo3Xr1kGO3DPeXt9nn33m8me3f/9+qNVq/O53vwty5J7x9vpWrFiBp556Cs8++ywOHjyIxYsX4+OPP8b/+3//L8iRN87ba5s3bx7effddvPnmm/jll1/w4IMP4rbbbsPu3buDHLlnysvL0a9fPyxcuNCj5588eRI333wzRo4ciT179uDRRx/FvffeG7Ifht5en8lkQsuWLTFv3jz069cvwNE1nbfXt3nzZowZMwbffPMN8vLyMHLkSEyYMCEkfz69vbaYmBjMnj0bmzdvxsGDBzFv3jzMmzcP7733XoAj9ZLYTA0ePFicNWuW/LXVahUzMjLEBQsWuH1+bm6umJCQEKToms7b61u0aJHYsWNH0Ww2ByvEJvH2+mp77bXXxLi4OLGsrCxQITaJt9c3a9YscdSoUS73zZkzRxw6dGhA4/SFt9eWnp4uvvXWWy733X777eI999wT0Dj9AYD4+eefN/icJ554QuzVq5fLfXfeeac4bty4AEbmH55cn7Phw4eLf/jDHwIWj795e32Snj17ivPnz/d/QH7k67Xddttt4uTJk/0fUBM0y4qM2WxGXl4eRo8eLd+nUqkwevRobNmypd7vKysrQ7t27ZCZmYlbbrkFBw4cCEa4XvPl+r788ktkZ2dj1qxZSE1NRe/evfHSSy/BarUGK2yP+frn52zx4sW46667EBMTE6gwfebL9V133XXIy8uTl2hOnDiBb775BjfddFNQYvaUL9dmMpnqLONGRUXhxx9/DGiswbJlyxaX3w8AGDdunMc/yxRabDYbSktLkZSUpHQofrd792789NNPGD58uNKhuGiWicylS5dgtVqRmprqcn9qaioKCgrcfk+3bt2wZMkS/Pvf/8aHH34Im82G6667DufOnQtGyF7x5fpOnDiBTz75BFarFd988w2eeeYZvPLKK3jxxReDEbJXfLk+Z9u3b8f+/ftx7733BirEJvHl+u6++248//zzGDZsGLRaLTp16oQRI0aE3NKSL9c2btw4vPrqqzh69ChsNhvWrVsnLxVGgoKCAre/H0ajEZWVlQpFRb76+9//jrKyMtxxxx1Kh+I3bdq0gV6vx8CBAzFr1qyQ+7uzWSYyvsjOzsbUqVPRv39/DB8+HJ999hlatmyJd999V+nQ/MJms6FVq1Z47733kJWVhTvvvBNPP/003nnnHaVD87vFixejT58+GDx4sNKh+M3GjRvx0ksv4e2338auXbvw2Wef4T//+Q9eeOEFpUNrsjfeeANdunRB9+7dodPpMHv2bEyfPh0qFf/6otCyYsUKzJ8/H//617/QqlUrpcPxmx9++AE7d+7EO++8g9dffx0rV65UOiQXGqUDUEJKSgrUajUKCwtd7i8sLERaWppHr6HVajFgwAAcO3YsECE2iS/Xl56eDq1WC7VaLd/Xo0cPFBQUwGw2Q6fTBTRmbzTlz6+8vByrVq3C888/H8gQm8SX63vmmWcwZcoU+V9Kffr0QXl5Oe6//348/fTTIfOh78u1tWzZEl988QWqqqpw+fJlZGRk4KmnnkLHjh2DEXLApaWluf39iI+PR1RUlEJRkbdWrVqFe++9F6tXr66zVBjuOnToAMD+90phYSGee+45TJo0SeGoaoTG325BptPpkJWVhQ0bNsj32Ww2bNiwAdnZ2R69htVqxb59+5Cenh6oMH3my/UNHToUx44dg81mk+87cuQI0tPTQyqJAZr257d69WqYTCZMnjw50GH6zJfrq6ioqJOsSEmpGELHqTXlz85gMKB169aorq7Gp59+iltuuSXQ4QZFdna2y+8HAKxbt87jv4tIeStXrsT06dOxcuXKkB7p4A82mw0mk0npMFwp3W2slFWrVol6vV5cunSp+Msvv4j333+/mJiYKBYUFIiiKIpTpkwRn3rqKfn58+fPF9euXSseP35czMvLE++66y7RYDCIBw4cUOoSGuTt9Z05c0aMi4sTZ8+eLR4+fFj8+uuvxVatWokvvviiUpfQIG+vTzJs2DDxzjvvDHa4XvP2+p599lkxLi5OXLlypXjixAnxv//9r9ipUyfxjjvuUOoS6uXttW3dulX89NNPxePHj4ubN28WR40aJXbo0EG8evWqQlfQsNLSUnH37t3i7t27RQDiq6++Ku7evVs8ffq0KIqi+NRTT4lTpkyRn3/ixAkxOjpa/NOf/iQePHhQXLhwoahWq8U1a9YodQkN8vb6RFGUn5+VlSXefffd4u7du0P2705vr++jjz4SNRqNuHDhQjE/P1++FRcXK3UJ9fL22t566y3xyy+/FI8cOSIeOXJE/Oc//ynGxcWJTz/9tFKX4FazTWREURTffPNNsW3btqJOpxMHDx4sbt26VX5s+PDhYk5Ojvz1o48+Kj83NTVVvOmmm8Rdu3YpELXnvLk+URTFn376SRwyZIio1+vFjh07in/5y1/E6urqIEftOW+v79ChQyIA8b///W+QI/WNN9dnsVjE5557TuzUqZNoMBjEzMxM8fe//33Ifth7c20bN24Ue/ToIer1ejE5OVmcMmWKeP78eQWi9sz3338vAqhzk64pJydHHD58eJ3v6d+/v6jT6cSOHTuKubm5QY/bU75cn7vnt2vXLuixe8Lb6xs+fHiDzw8l3l7bP/7xD7FXr15idHS0GB8fLw4YMEB8++23RavVqswF1EMQxRCqOxMRERF5oVn2yBAREVFkYCJDREREYYuJDBEREYUtJjJEREQUtpjIEBERUdhiIkNERERhi4kMERERhS0mMkRERBS2mMgQERFR2GIiQ0QBN23aNAiCUOd24403evwap06dgiAIUKvVOH/+vMtj+fn50Gg0EAQBp06d8uj1JkyYUO/7//DDDxAEAXv37vU4PiJSBhMZIgqKG2+8Efn5+S63lStXev06rVu3xvLly13uW7ZsGVq3bu3V68ycORPr1q3DuXPn6jyWm5uLgQMHom/fvl7HR0TBxUSGiIJCr9cjLS3N5daiRQuvXycnJwe5ubku9+Xm5iInJ6fOc/fv34/x48cjNjYWqampmDJlCi5dugQA+PWvf42WLVti6dKlLt9TVlaG1atXY+bMmV7HRkTBx0SGiMLKb37zG1y9ehU//vgjAODHH3/E1atXMWHCBJfnFRcXY9SoURgwYAB27tyJNWvWoLCwEHfccQcAQKPRYOrUqVi6dCmcz85dvXo1rFYrJk2aFLyLIiKfMZEhoqD4+uuvERsb63J76aWXvH4drVaLyZMnY8mSJQCAJUuWYPLkydBqtS7Pe+uttzBgwAC89NJL6N69OwYMGIAlS5bg+++/x5EjRwAAM2bMwPHjx7Fp0yb5+3JzczFx4kQkJCQ04WqJKFg0SgdARM3DyJEjsWjRIpf7kpKSfHqtGTNm4LrrrsNLL72E1atXY8uWLaiurnZ5zs8//4zvv/8esbGxdb7/+PHj6Nq1K7p3747rrrsOS5YswYgRI3Ds2DH88MMPeP75532Ki4iCj4kMEQVFTEwMOnfu7JfX6tOnD7p3745JkyahR48e6N27N/bs2ePynLKyMkyYMAEvv/xyne9PT0+Xfz1z5kw8/PDDWLhwIXJzc9GpUycMHz7cL3ESUeBxaYmIwtKMGTOwceNGzJgxw+3j11xzDQ4cOID27dujc+fOLreYmBj5eXfccQdUKhVWrFiB5cuXY8aMGRAEIViXQURNxESGiILCZDKhoKDA5SbtIPLFfffdh4sXL+Lee+91+/isWbNw5coVTJo0CTt27MDx48exdu1aTJ8+HVarVX5ebGws7rzzTsydOxf5+fmYNm2azzERUfAxkSGioFizZg3S09NdbsOGDZMfnzZtGkaMGOHx62k0GqSkpECjcb9CnpGRgf/973+wWq0YO3Ys+vTpg0cffRSJiYlQqVz/6ps5cyauXr2KcePGISMjw6frIyJlCKLzvkMiIoUMHz4cI0eOxHPPPad0KEQURpjIEJHiSkpK0KtXLxw6dMjtLiMiovowkSEiIqKwxR4ZIiIiCltMZIiIiChsMZEhIiKisMVEhoiIiMIWExkiIiIKW0xkiIiIKGwxkSEiIqKwxUSGiIiIwhYTGSIiIgpb/x8PJsVLP7vyOgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot((gamma(correct_df['Speed']))*me*c**2/1.6*10**13, correct_df['Count'])\n",
"plt.ylabel('Count')\n",
"plt.xlabel('E, MeV')\n",
"plt.show()"
]
}
],
"metadata": {
"colab": {
"provenance": []
},
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"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.10.12"
}
},
"nbformat": 4,
"nbformat_minor": 4
}