{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Diamond detector" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import matplotlib.ticker as mticker\n", "from libra_toolbox.neutron_detection.diamond.process_data import DataProcessor" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Fabricate data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "total_time_s = 200 # s\n", "total_time_ps = total_time_s * 1e12 # s to ps\n", "\n", "# peak 1\n", "nb_counts_peak1 = int(2e4)\n", "size_peak1 = nb_counts_peak1\n", "time_values_out_peak1 = np.random.rand(size_peak1) * total_time_ps\n", "mean_energy_peak1 = 4e6\n", "std_energy_peak1 = 0.5e6\n", "energy_values_peak1 = np.random.normal(mean_energy_peak1, std_energy_peak1, size_peak1)\n", "\n", "# peak 2\n", "nb_counts_peak2 = int(7e4)\n", "size_peak2 = nb_counts_peak2\n", "time_values_out_peak2 = np.random.rand(size_peak2) * total_time_ps\n", "mean_energy_peak2 = 14e6\n", "std_energy_peak2 = 1e6\n", "\n", "energy_values_peak2 = np.random.normal(mean_energy_peak2, std_energy_peak2, size_peak2)\n", "\n", "time_values_out = np.concatenate((time_values_out_peak1, time_values_out_peak2))\n", "energy_values = np.concatenate((energy_values_peak1, energy_values_peak2))\n", "\n", "# remove values where t between 100 and 150 s\n", "mask = (time_values_out > 100e12) & (time_values_out < 150e12)\n", "time_values_out = time_values_out[~mask]\n", "energy_values = energy_values[~mask]\n", "\n", "# write data to files\n", "# make data directory if it doesn't exist\n", "import os\n", "\n", "if not os.path.exists(\"data\"):\n", " os.makedirs(\"data\")\n", "\n", "np.savetxt(\n", " \"data/data.csv\",\n", " np.column_stack((time_values_out, energy_values)),\n", " delimiter=\",\",\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read and process data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Added file: data/data.csv containing 67502 events\n" ] } ], "source": [ "data_proc = DataProcessor()\n", "data_proc.add_file(\"data/data.csv\", time_column=0, energy_column=1, delimiter=\",\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "res_avg = data_proc.get_avg_rate(0, 100)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot results" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAG0CAYAAADTmjjeAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA7UUlEQVR4nO3de1hU9d7//xcHGY+DmsFAkpmKiomVp2aX5gFBI+/c2Z2WKTs1s43thFLj1jS1suxgJ9NrZ4l7p5l6p7utpiEe2imdKPKMebq1rQOlwaApKKzvH/vH/BohgxGYgfV8XNe6Ltdan1nzfk8xvPjMWmv8DMMwBAAAYGL+3i4AAADA2whEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9LwaiBYsWKDo6GhZrVZZrVbZ7XZ9/PHHrv19+vSRn5+f2zJ+/Hi3Yxw7dkzx8fFq2LChQkJCNGnSJF28eNFtzNatW3XzzTfLYrGobdu2Sk1NrYn2AABALRHozSdv2bKlnn/+ebVr106GYWjJkiW666679O2336pTp06SpIceekizZs1yPaZhw4aufxcXFys+Pl42m007duzQyZMnNWrUKNWrV0/PPfecJOnIkSOKj4/X+PHjtXTpUqWnp2vs2LEKCwtTXFxcheosKSnRiRMn1KRJE/n5+VXhKwAAAKqLYRgqKChQeHi4/P1/Zw7I8DHNmjUzFi1aZBiGYdx+++3GY4899ptj169fb/j7+xsOh8O1bcGCBYbVajUKCwsNwzCMyZMnG506dXJ73LBhw4y4uLgK13T8+HFDEgsLCwsLC0stXI4fP/67v+u9OkP0a8XFxVq5cqXOnj0ru93u2r506VK99957stlsGjx4sJ566inXLFFGRoY6d+6s0NBQ1/i4uDg98sgj2rNnj2666SZlZGQoJibG7bni4uI0ceLE36ylsLBQhYWFrnXDMCRJx48fl9VqrYp2AQBANXM6nYqIiFCTJk1+d6zXA9GuXbtkt9t1/vx5NW7cWKtXr1ZUVJQk6f7771erVq0UHh6unTt3asqUKcrOztaHH34oSXI4HG5hSJJr3eFwXHaM0+nUuXPn1KBBgzI1zZkzRzNnziyzvfRcJwAAUHtU5HQXrwei9u3bKysrS/n5+Vq1apUSEhK0bds2RUVFady4ca5xnTt3VlhYmPr3769Dhw6pTZs21VZTSkqKkpOTXeulCRMAANRNXr/sPigoSG3btlXXrl01Z84cdenSRa+99lq5Y3v27ClJOnjwoCTJZrMpJyfHbUzpus1mu+wYq9Va7uyQJFksFtdsELNCAADUfV4PRJcqKSlxO3/n17KysiRJYWFhkiS73a5du3YpNzfXNSYtLU1Wq9X1sZvdbld6errbcdLS0tzOUwIAAObm1Y/MUlJSNGjQIF177bUqKCjQsmXLtHXrVm3cuFGHDh3SsmXLdMcdd+iqq67Szp07lZSUpN69eys6OlqSFBsbq6ioKI0cOVJz586Vw+HQtGnTlJiYKIvFIkkaP3683nzzTU2ePFmjR4/W5s2btWLFCq1bt86brQMAAB/i1UCUm5urUaNG6eTJkwoODlZ0dLQ2btyoAQMG6Pjx49q0aZNeffVVnT17VhERERo6dKimTZvmenxAQIDWrl2rRx55RHa7XY0aNVJCQoLbfYtat26tdevWKSkpSa+99ppatmypRYsWVfgeRAAAoO7zM0qvKcdvcjqdCg4OVn5+PucTAQBQS1Tm97fPnUMEAABQ0whEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9Lz61R0AAHhiXtqB3x2TNCCyBipBXcEMEQAAMD1miAAAddKls0jMGOFymCECAACmRyACAACmRyACAACmxzlEAACfVpEryoArxQwRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwvUBvFwAAwK/NSzvg7RJgQswQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0+MqMwCAKZR39VrSgEgvVAJfxAwRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPa8GogULFig6OlpWq1VWq1V2u10ff/yxa//58+eVmJioq666So0bN9bQoUOVk5Pjdoxjx44pPj5eDRs2VEhIiCZNmqSLFy+6jdm6datuvvlmWSwWtW3bVqmpqTXRHgAAqCW8Gohatmyp559/XpmZmfr666/Vr18/3XXXXdqzZ48kKSkpSf/85z+1cuVKbdu2TSdOnNDdd9/tenxxcbHi4+NVVFSkHTt2aMmSJUpNTdX06dNdY44cOaL4+Hj17dtXWVlZmjhxosaOHauNGzfWeL8AAMA3+RmGYXi7iF9r3ry5XnzxRd1zzz26+uqrtWzZMt1zzz2SpP3796tjx47KyMjQLbfcoo8//lh33nmnTpw4odDQUEnSwoULNWXKFP34448KCgrSlClTtG7dOu3evdv1HMOHD1deXp42bNhQbg2FhYUqLCx0rTudTkVERCg/P19Wq7UauwcA1OSXu3LZfd3mdDoVHBxcod/fPnMOUXFxsZYvX66zZ8/KbrcrMzNTFy5cUExMjGtMhw4ddO211yojI0OSlJGRoc6dO7vCkCTFxcXJ6XS6ZpkyMjLcjlE6pvQY5ZkzZ46Cg4NdS0RERFW2CgAAfIzXA9GuXbvUuHFjWSwWjR8/XqtXr1ZUVJQcDoeCgoLUtGlTt/GhoaFyOBySJIfD4RaGSveX7rvcGKfTqXPnzpVbU0pKivLz813L8ePHq6JVAADgo7x+p+r27dsrKytL+fn5WrVqlRISErRt2zav1mSxWGSxWLxaAwAAqDleD0RBQUFq27atJKlr16766quv9Nprr2nYsGEqKipSXl6e2yxRTk6ObDabJMlms+nLL790O17pVWi/HnPplWk5OTmyWq1q0KBBdbUFAABqEa9/ZHapkpISFRYWqmvXrqpXr57S09Nd+7Kzs3Xs2DHZ7XZJkt1u165du5Sbm+sak5aWJqvVqqioKNeYXx+jdEzpMQAAALw6Q5SSkqJBgwbp2muvVUFBgZYtW6atW7dq48aNCg4O1pgxY5ScnKzmzZvLarXq0Ucfld1u1y233CJJio2NVVRUlEaOHKm5c+fK4XBo2rRpSkxMdH3kNX78eL355puaPHmyRo8erc2bN2vFihVat26dN1sHAAA+xKuBKDc3V6NGjdLJkycVHBys6Ohobdy4UQMGDJAkzZs3T/7+/ho6dKgKCwsVFxent956y/X4gIAArV27Vo888ojsdrsaNWqkhIQEzZo1yzWmdevWWrdunZKSkvTaa6+pZcuWWrRokeLi4mq8XwAA4Jt87j5Evqgy9zEAAFwZ7kOEqlIr70MEAADgLQQiAABgegQiAABgegQiAABgegQiAABgegQiAABgegQiAABgegQiAABgel7/clcAgHnV5E0YgcthhggAAJgeM0QAANO6dIaKr/IwL2aIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6QV6uwAAAHzFvLQDZbYlDYj0QiWoacwQAQAA02OGCABQY8qbgQF8ATNEAADA9AhEAADA9AhEAADA9LwaiObMmaPu3burSZMmCgkJ0ZAhQ5Sdne02pk+fPvLz83Nbxo8f7zbm2LFjio+PV8OGDRUSEqJJkybp4sWLbmO2bt2qm2++WRaLRW3btlVqamp1twcAAGoJrwaibdu2KTExUZ9//rnS0tJ04cIFxcbG6uzZs27jHnroIZ08edK1zJ0717WvuLhY8fHxKioq0o4dO7RkyRKlpqZq+vTprjFHjhxRfHy8+vbtq6ysLE2cOFFjx47Vxo0ba6xXAADgu7x6ldmGDRvc1lNTUxUSEqLMzEz17t3btb1hw4ay2WzlHuOTTz7R3r17tWnTJoWGhurGG2/U7NmzNWXKFD399NMKCgrSwoUL1bp1a7388suSpI4dO+qzzz7TvHnzFBcXV30NAgCAWsGnziHKz8+XJDVv3txt+9KlS9WiRQvdcMMNSklJ0S+//OLal5GRoc6dOys0NNS1LS4uTk6nU3v27HGNiYmJcTtmXFycMjIyyq2jsLBQTqfTbQEAAHWXz9yHqKSkRBMnTtStt96qG264wbX9/vvvV6tWrRQeHq6dO3dqypQpys7O1ocffihJcjgcbmFIkmvd4XBcdozT6dS5c+fUoEEDt31z5szRzJkzq7xHAADgm3wmECUmJmr37t367LPP3LaPGzfO9e/OnTsrLCxM/fv316FDh9SmTZtqqSUlJUXJycmudafTqYiIiGp5LgAA4H0+8ZHZhAkTtHbtWm3ZskUtW7a87NiePXtKkg4ePChJstlsysnJcRtTul563tFvjbFarWVmhyTJYrHIarW6LQAAoO7yaiAyDEMTJkzQ6tWrtXnzZrVu3fp3H5OVlSVJCgsLkyTZ7Xbt2rVLubm5rjFpaWmyWq2KiopyjUlPT3c7Tlpamux2exV1AgAAajOvBqLExES99957WrZsmZo0aSKHwyGHw6Fz585Jkg4dOqTZs2crMzNTR48e1UcffaRRo0apd+/eio6OliTFxsYqKipKI0eO1HfffaeNGzdq2rRpSkxMlMVikSSNHz9ehw8f1uTJk7V//3699dZbWrFihZKSkrzWOwAA8B1eDUQLFixQfn6++vTpo7CwMNfywQcfSJKCgoK0adMmxcbGqkOHDnr88cc1dOhQ/fOf/3QdIyAgQGvXrlVAQIDsdrseeOABjRo1SrNmzXKNad26tdatW6e0tDR16dJFL7/8shYtWsQl9wAAQJLkZxiG4e0ifJ3T6VRwcLDy8/M5nwgArkBt/Lb7pAGR3i4BHqrM72+fOKkaAADAmwhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9AK9XQAAAL5sXtoBt/WkAZFeqgTViRkiAABgegQiAABgegQiAABgegQiAABgepxUDQCoFpeejAz4MmaIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6RGIAACA6Xk1EM2ZM0fdu3dXkyZNFBISoiFDhig7O9ttzPnz55WYmKirrrpKjRs31tChQ5WTk+M25tixY4qPj1fDhg0VEhKiSZMm6eLFi25jtm7dqptvvlkWi0Vt27ZVampqdbcHAABqCa8Gom3btikxMVGff/650tLSdOHCBcXGxurs2bOuMUlJSfrnP/+plStXatu2bTpx4oTuvvtu1/7i4mLFx8erqKhIO3bs0JIlS5Samqrp06e7xhw5ckTx8fHq27evsrKyNHHiRI0dO1YbN26s0X4BAIBv8jMMw/B2EaV+/PFHhYSEaNu2berdu7fy8/N19dVXa9myZbrnnnskSfv371fHjh2VkZGhW265RR9//LHuvPNOnThxQqGhoZKkhQsXasqUKfrxxx8VFBSkKVOmaN26ddq9e7fruYYPH668vDxt2LChTB2FhYUqLCx0rTudTkVERCg/P19Wq7WaXwUAqBvmpR3wdgnVImlApLdLQAU5nU4FBwdX6Pe3T51DlJ+fL0lq3ry5JCkzM1MXLlxQTEyMa0yHDh107bXXKiMjQ5KUkZGhzp07u8KQJMXFxcnpdGrPnj2uMb8+RumY0mNcas6cOQoODnYtERERVdckAADwOT4TiEpKSjRx4kTdeuutuuGGGyRJDodDQUFBatq0qdvY0NBQORwO15hfh6HS/aX7LjfG6XTq3LlzZWpJSUlRfn6+azl+/HiV9AgAAHxToLcLKJWYmKjdu3frs88+83Ypslgsslgs3i4DAADUEJ+YIZowYYLWrl2rLVu2qGXLlq7tNptNRUVFysvLcxufk5Mjm83mGnPpVWel6783xmq1qkGDBlXdDgAAqGW8GogMw9CECRO0evVqbd68Wa1bt3bb37VrV9WrV0/p6emubdnZ2Tp27JjsdrskyW63a9euXcrNzXWNSUtLk9VqVVRUlGvMr49ROqb0GAAAwNw8CkTffPONdu3a5Vr/xz/+oSFDhuh//ud/VFRUVOHjJCYm6r333tOyZcvUpEkTORwOORwO13k9wcHBGjNmjJKTk7VlyxZlZmbqwQcflN1u1y233CJJio2NVVRUlEaOHKnvvvtOGzdu1LRp05SYmOj62Gv8+PE6fPiwJk+erP379+utt97SihUrlJSU5En7AACgjvEoED388MM6cOA/l1MePnxYw4cPV8OGDbVy5UpNnjy5wsdZsGCB8vPz1adPH4WFhbmWDz74wDVm3rx5uvPOOzV06FD17t1bNptNH374oWt/QECA1q5dq4CAANntdj3wwAMaNWqUZs2a5RrTunVrrVu3TmlpaerSpYtefvllLVq0SHFxcZ60DwAA6hiP7kMUHBysb775Rm3atNELL7ygzZs3a+PGjdq+fbuGDx9e567Kqsx9DADArOrqfYcuxX2Iao/K/P726CozwzBUUlIiSdq0aZPuvPNOSVJERIR++uknTw4JAECtUF7wIyTVfh59ZNatWzc988wz+vvf/65t27YpPj5e0n++IuPS+/0AAAD4Oo9miObNm6cHHnhAa9as0dSpU9W2bVtJ0qpVq/SHP/yhSgtE7XHpX038xQQAqC08CkRdunRxu8qs1IsvvqjAQJ+51yMAAECFePSR2fXXX69Tp06V2X7+/HlFRjIrAAAAahePpnOOHj2q4uLiMtsLCwv1ww8/XHFRqBs48RAAUFtUKhB99NFHrn9v3LhRwcHBrvXi4mKlp6eXuds0AACAr6tUIBoyZIgkyc/PTwkJCW776tWrp+uuu04vv/xylRUHAABQEyoViErvPdS6dWt99dVXatGiRbUUBQAAUJM8OofoyJEjVV0HAACA13h8jXx6errS09OVm5vrmjkq9e67715xYQAAADXFo0A0c+ZMzZo1S926dVNYWJj8/Pyqui4AAIAa41EgWrhwoVJTUzVy5MiqrgcAAKDGeXRjxqKiIr6iAwAA1BkeBaKxY8dq2bJlVV0LAACAV3j0kdn58+f117/+VZs2bVJ0dLTq1avntv+VV16pkuIAAABqgkeBaOfOnbrxxhslSbt373bbxwnWAACgtvEoEG3ZsqWq6wAAAPAaj84hAgAAqEs8miHq27fvZT8a27x5s8cFAQAA1DSPAlHp+UOlLly4oKysLO3evbvMl74CAAD4Oo8C0bx588rd/vTTT+vMmTNXVBAAAEBNq9JziB544AG+xwwAANQ6VRqIMjIyVL9+/ao8JAAAQLXz6COzu+++223dMAydPHlSX3/9tZ566qkqKQwAAKCmeBSIgoOD3db9/f3Vvn17zZo1S7GxsVVSGAAAQE3xKBAtXry4qusAAADwGo8CUanMzEzt27dPktSpUyfddNNNVVIUAAC1yby0A27rSQMivVQJPOVRIMrNzdXw4cO1detWNW3aVJKUl5envn37avny5br66qurskYAAIBq5VEgevTRR1VQUKA9e/aoY8eOkqS9e/cqISFBf/nLX/T+++9XaZEAAN9y6YwIUNt5FIg2bNigTZs2ucKQJEVFRWn+/PmcVA0AAGodjwJRSUmJ6tWrV2Z7vXr1VFJScsVFwffx1yEAoC7xKBD169dPjz32mN5//32Fh4dLkv79738rKSlJ/fv3r9ICUbdw4iEAwBd5dKfqN998U06nU9ddd53atGmjNm3aqHXr1nI6nXrjjTequkYAAIBq5dEMUUREhL755htt2rRJ+/fvlyR17NhRMTExVVocAABATajUDNHmzZsVFRUlp9MpPz8/DRgwQI8++qgeffRRde/eXZ06ddK//vWv6qoVAACgWlQqEL366qt66KGHZLVay+wLDg7Www8/rFdeeaXKigMAAKgJlQpE3333nQYOHPib+2NjY5WZmXnFRQEAANSkSgWinJycci+3LxUYGKgff/zxiosCAACoSZUKRNdcc4127979m/t37typsLCwKy4KAACgJlUqEN1xxx166qmndP78+TL7zp07pxkzZujOO++ssuIAAABqQqUC0bRp03T69GlFRkZq7ty5+sc//qF//OMfeuGFF9S+fXudPn1aU6dOrfDxPv30Uw0ePFjh4eHy8/PTmjVr3Pb/6U9/kp+fn9ty6TlMp0+f1ogRI2S1WtW0aVONGTNGZ86ccRuzc+dO9erVS/Xr11dERITmzp1bmbYBAEAdV6n7EIWGhmrHjh165JFHlJKSIsMwJEl+fn6Ki4vT/PnzFRoaWuHjnT17Vl26dNHo0aN19913lztm4MCBWrx4sWvdYrG47R8xYoROnjyptLQ0XbhwQQ8++KDGjRunZcuWSZKcTqdiY2MVExOjhQsXateuXRo9erSaNm2qcePGVaZ9AABQR1X6xoytWrXS+vXr9fPPP+vgwYMyDEPt2rVTs2bNKv3kgwYN0qBBgy47xmKxyGazlbtv37592rBhg7766it169ZNkvTGG2/ojjvu0EsvvaTw8HAtXbpURUVFevfddxUUFKROnTopKytLr7zyCoEIAABI8vCrOySpWbNm6t69u3r06OFRGKqorVu3KiQkRO3bt9cjjzyiU6dOufZlZGSoadOmrjAkSTExMfL399cXX3zhGtO7d28FBQW5xsTFxSk7O1s///xzuc9ZWFgop9PptgAAgLrL40BUEwYOHKi//e1vSk9P1wsvvKBt27Zp0KBBKi4uliQ5HA6FhIS4PSYwMFDNmzeXw+Fwjbn0Y7zS9dIxl5ozZ46Cg4NdS0RERFW3BgAAfIhH32VWU4YPH+76d+fOnRUdHa02bdpo69at6t+/f7U9b0pKipKTk13rTqeTUAQAQB3m0zNEl7r++uvVokULHTx4UJJks9mUm5vrNubixYs6ffq067wjm82mnJwctzGl6791bpLFYpHVanVbAABA3VWrAtEPP/ygU6dOuW7+aLfblZeX5/Z1IZs3b1ZJSYl69uzpGvPpp5/qwoULrjFpaWlq3759tZ77BAAAag+vBqIzZ84oKytLWVlZkqQjR44oKytLx44d05kzZzRp0iR9/vnnOnr0qNLT03XXXXepbdu2iouLkyR17NhRAwcO1EMPPaQvv/xS27dv14QJEzR8+HCFh4dLku6//34FBQVpzJgx2rNnjz744AO99tprbh+JAQAAc/PqOURff/21+vbt61ovDSkJCQlasGCBdu7cqSVLligvL0/h4eGKjY3V7Nmz3e5FtHTpUk2YMEH9+/eXv7+/hg4dqtdff921Pzg4WJ988okSExPVtWtXtWjRQtOnT+eSex8xL+1AmW1JAyK9UAkAwMy8Goj69OnjurljeTZu3Pi7x2jevLnrJoy/JTo6Wv/6178qXR8AADCHWnUOEQAAQHUgEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANPz6pe7ovYo71vpAZgH7wGo65ghAgAApscMEQAAVay8GbWkAZFeqAQVxQwRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwvUBvFwBcal7aAbf1pAGRXqoEMKdLfwYBM2CGCAAAmB6BCAAAmB6BCAAAmB6BCAAAmB6BCAAAmB6BCAAAmJ5XA9Gnn36qwYMHKzw8XH5+flqzZo3bfsMwNH36dIWFhalBgwaKiYnR999/7zbm9OnTGjFihKxWq5o2baoxY8bozJkzbmN27typXr16qX79+oqIiNDcuXOruzUAAFCLePU+RGfPnlWXLl00evRo3X333WX2z507V6+//rqWLFmi1q1b66mnnlJcXJz27t2r+vXrS5JGjBihkydPKi0tTRcuXNCDDz6ocePGadmyZZIkp9Op2NhYxcTEaOHChdq1a5dGjx6tpk2baty4cTXab23BPUgAAGbj1UA0aNAgDRo0qNx9hmHo1Vdf1bRp03TXXXdJkv72t78pNDRUa9as0fDhw7Vv3z5t2LBBX331lbp16yZJeuONN3THHXfopZdeUnh4uJYuXaqioiK9++67CgoKUqdOnZSVlaVXXnmFQAQAACT58DlER44ckcPhUExMjGtbcHCwevbsqYyMDElSRkaGmjZt6gpDkhQTEyN/f3998cUXrjG9e/dWUFCQa0xcXJyys7P1888/l/vchYWFcjqdbgsAAFdiXtoBtwW+xWcDkcPhkCSFhoa6bQ8NDXXtczgcCgkJcdsfGBio5s2bu40p7xi/fo5LzZkzR8HBwa4lIiLiyhsCAAA+y2cDkTelpKQoPz/ftRw/ftzbJQEAgGrks4HIZrNJknJycty25+TkuPbZbDbl5ua67b948aJOnz7tNqa8Y/z6OS5lsVhktVrdFgAAUHf5bCBq3bq1bDab0tPTXducTqe++OIL2e12SZLdbldeXp4yMzNdYzZv3qySkhL17NnTNebTTz/VhQsXXGPS0tLUvn17NWvWrIa6AQAAvsyrgejMmTPKyspSVlaWpP+cSJ2VlaVjx47Jz89PEydO1DPPPKOPPvpIu3bt0qhRoxQeHq4hQ4ZIkjp27KiBAwfqoYce0pdffqnt27drwoQJGj58uMLDwyVJ999/v4KCgjRmzBjt2bNHH3zwgV577TUlJyd7qWsAAOBrvHrZ/ddff62+ffu61ktDSkJCglJTUzV58mSdPXtW48aNU15enm677TZt2LDBdQ8iSVq6dKkmTJig/v37y9/fX0OHDtXrr7/u2h8cHKxPPvlEiYmJ6tq1q1q0aKHp06dzyT0AAHDxMwzD8HYRvs7pdCo4OFj5+fmmOJ/I1y4HTRoQ6e0SAFPxtfeAuor3tupXmd/fPnsOEQAAQE0hEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANMjEAEAANPz6ld3AABgVuXdEZy7V3sPgQgATI6v6gD4yAwAAIBABAAAQCACAACmRyACAACmRyACAACmRyACAACmRyACAACmRyACAACmRyACAACmx52q4fO4vT0AoLoxQwQAAEyPQAQAAEyPQAQAAEyPQAQAAEyPQAQAAEyPQAQAAEyPy+5R7mXtAACYCTNEAADA9AhEAADA9AhEAADA9AhEAADA9AhEAADA9LjKDABMhKtKgfIxQwQAAEyPGSIAAHzEpTN4SQMivVSJ+TBDBAAATI9ABAAATI9ABAAATI9ABAAATI9ABAAATM+nA9HTTz8tPz8/t6VDhw6u/efPn1diYqKuuuoqNW7cWEOHDlVOTo7bMY4dO6b4+Hg1bNhQISEhmjRpki5evFjTrQAAAB/m85fdd+rUSZs2bXKtBwb+/yUnJSVp3bp1WrlypYKDgzVhwgTdfffd2r59uySpuLhY8fHxstls2rFjh06ePKlRo0apXr16eu6552q8FwAA4Jt8PhAFBgbKZrOV2Z6fn6933nlHy5YtU79+/SRJixcvVseOHfX555/rlltu0SeffKK9e/dq06ZNCg0N1Y033qjZs2drypQpevrppxUUFFTT7QAAAB/k0x+ZSdL333+v8PBwXX/99RoxYoSOHTsmScrMzNSFCxcUExPjGtuhQwdde+21ysjIkCRlZGSoc+fOCg0NdY2Ji4uT0+nUnj17fvM5CwsL5XQ63RYAAFB3+XQg6tmzp1JTU7VhwwYtWLBAR44cUa9evVRQUCCHw6GgoCA1bdrU7TGhoaFyOBySJIfD4RaGSveX7vstc+bMUXBwsGuJiIio2sYAAIBP8emPzAYNGuT6d3R0tHr27KlWrVppxYoVatCgQbU9b0pKipKTk13rTqeTUAQAQB3m0zNEl2ratKkiIyN18OBB2Ww2FRUVKS8vz21MTk6O65wjm81W5qqz0vXyzksqZbFYZLVa3RYAAFB31apAdObMGR06dEhhYWHq2rWr6tWrp/T0dNf+7OxsHTt2THa7XZJkt9u1a9cu5ebmusakpaXJarUqKiqqxusHAAC+yac/MnviiSc0ePBgtWrVSidOnNCMGTMUEBCg++67T8HBwRozZoySk5PVvHlzWa1WPfroo7Lb7brlllskSbGxsYqKitLIkSM1d+5cORwOTZs2TYmJibJYLF7uDgCq36Xfno7apbz/fkkDIr1QSd3n04Hohx9+0H333adTp07p6quv1m233abPP/9cV199tSRp3rx58vf319ChQ1VYWKi4uDi99dZbrscHBARo7dq1euSRR2S329WoUSMlJCRo1qxZ3moJVeTSNwneIAAAV8LPMAzD20X4OqfTqeDgYOXn59fJ84nqwl+QBCKgfHXh5xvueL+ruMr8/q5V5xABAABUBwIRAAAwPQIRAAAwPZ8+qRpVj/MJAAAoixkiAABgegQiAABgegQiAABgegQiAABgegQiAABgegQiAABgelx2DwB1BLfVADzHDBEAADA9ZogAAKhFLp0J5MteqwYzRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPQIRAAAwPS4ygwAgFqsvPtPceVZ5RGIUCfwhgAAuBJ8ZAYAAEyPQAQAAEyPQAQAAEyPQAQAAEyPk6rrML75Gqi7+PkGqhYzRAAAwPQIRAAAwPQIRAAAwPQ4hwgAagHOGQKqFzNEAADA9JghAgCgjrl0RpGvMvp9zBABAADTIxABAADTIxABAADT4xwi1Fl8hg4AqCgCEQAAdVx5t23gj0R3BKI6hPuUAHUDP8tAzeMcIgAAYHoEIgAAYHoEIgAAYHqcQwTT4KRC+CrOGYI3cCWuOwIRANQgwg/gm0wViObPn68XX3xRDodDXbp00RtvvKEePXp4uyyP8KYKAKhKFfm9UpdnkUwTiD744AMlJydr4cKF6tmzp1599VXFxcUpOztbISEh3i4PXsKUMaoSf6gAtZefYRiGt4uoCT179lT37t315ptvSpJKSkoUERGhRx99VE8++eRlH+t0OhUcHKz8/HxZrdaaKLcM3mhrBoEIV4KfU8C33kcr8/vbFDNERUVFyszMVEpKimubv7+/YmJilJGRUWZ8YWGhCgsLXev5+fmS/vPCVof5mw9Wy3FReXPWfOPR4xL7ta3iSuBr+DkFKqa891FvvUeW/t6uyNyPKQLRTz/9pOLiYoWGhrptDw0N1f79+8uMnzNnjmbOnFlme0RERLXViNrtf7xdAAD4MG+/RxYUFCg4OPiyY0wRiCorJSVFycnJrvWSkhKdPn1aV111lfz8/KrkOZxOpyIiInT8+HGvfQxX3ep6j3W9P4ke64q63mNd70+iR08ZhqGCggKFh4f/7lhTBKIWLVooICBAOTk5bttzcnJks9nKjLdYLLJYLG7bmjZtWi21Wa3WOvs/d6m63mNd70+ix7qirvdY1/uT6NETvzczVMoUd6oOCgpS165dlZ6e7tpWUlKi9PR02e12L1YGAAB8gSlmiCQpOTlZCQkJ6tatm3r06KFXX31VZ8+e1YMPPujt0gAAgJeZJhANGzZMP/74o6ZPny6Hw6Ebb7xRGzZsKHOidU2xWCyaMWNGmY/m6pK63mNd70+ix7qirvdY1/uT6LEmmOY+RAAAAL/FFOcQAQAAXA6BCAAAmB6BCAAAmB6BCAAAmB6BqBrNnz9f1113nerXr6+ePXvqyy+//M2xb7/9tnr16qVmzZqpWbNmiomJuex4X1GZHn9t+fLl8vPz05AhQ6q3wCtU2f7y8vKUmJiosLAwWSwWRUZGav369TVUrWcq2+Orr76q9u3bq0GDBoqIiFBSUpLOnz9fQ9VWzqeffqrBgwcrPDxcfn5+WrNmze8+ZuvWrbr55ptlsVjUtm1bpaamVnudV6KyPX744YcaMGCArr76almtVtntdm3cuLFmivWQJ/8dS23fvl2BgYG68cYbq62+quBJj4WFhZo6dapatWoli8Wi6667Tu+++271F+sBT/pbunSpunTpooYNGyosLEyjR4/WqVOnqq1GAlE1+eCDD5ScnKwZM2bom2++UZcuXRQXF6fc3Nxyx2/dulX33XeftmzZooyMDEVERCg2Nlb//ve/a7jyiqtsj6WOHj2qJ554Qr169aqhSj1T2f6Kioo0YMAAHT16VKtWrVJ2drbefvttXXPNNTVcecVVtsdly5bpySef1IwZM7Rv3z698847+uCDD/Q//+Ptbyoq39mzZ9WlSxfNnz+/QuOPHDmi+Ph49e3bV1lZWZo4caLGjh3r04Ghsj1++umnGjBggNavX6/MzEz17dtXgwcP1rffflvNlXqusj2WysvL06hRo9S/f/9qqqzqeNLjvffeq/T0dL3zzjvKzs7W+++/r/bt21djlZ6rbH/bt2/XqFGjNGbMGO3Zs0crV67Ul19+qYceeqj6ijRQLXr06GEkJia61ouLi43w8HBjzpw5FXr8xYsXjSZNmhhLliyprhKvmCc9Xrx40fjDH/5gLFq0yEhISDDuuuuuGqjUM5Xtb8GCBcb1119vFBUV1VSJV6yyPSYmJhr9+vVz25acnGzceuut1VpnVZBkrF69+rJjJk+ebHTq1Mlt27Bhw4y4uLhqrKzqVKTH8kRFRRkzZ86s+oKqQWV6HDZsmDFt2jRjxowZRpcuXaq1rqpUkR4//vhjIzg42Dh16lTNFFWFKtLfiy++aFx//fVu215//XXjmmuuqba6mCGqBkVFRcrMzFRMTIxrm7+/v2JiYpSRkVGhY/zyyy+6cOGCmjdvXl1lXhFPe5w1a5ZCQkI0ZsyYmijTY57099FHH8lutysxMVGhoaG64YYb9Nxzz6m4uLimyq4UT3r8wx/+oMzMTNfHaocPH9b69et1xx131EjN1S0jI8Pt9ZCkuLi4Cv/c1kYlJSUqKCjw2fcaTy1evFiHDx/WjBkzvF1Ktfjoo4/UrVs3zZ07V9dcc40iIyP1xBNP6Ny5c94urUrY7XYdP35c69evl2EYysnJ0apVq6r1vcY0d6quST/99JOKi4vL3AU7NDRU+/fvr9AxpkyZovDw8DJvzr7Ckx4/++wzvfPOO8rKyqqBCq+MJ/0dPnxYmzdv1ogRI7R+/XodPHhQf/7zn3XhwgWffFP2pMf7779fP/30k2677TYZhqGLFy9q/PjxPvuRWWU5HI5yXw+n06lz586pQYMGXqqs+rz00ks6c+aM7r33Xm+XUmW+//57Pfnkk/rXv/6lwMC6+Wvu8OHD+uyzz1S/fn2tXr1aP/30k/785z/r1KlTWrx4sbfLu2K33nqrli5dqmHDhun8+fO6ePGiBg8eXOmPTSuDGSIf9Pzzz2v58uVavXq16tev7+1yqkRBQYFGjhypt99+Wy1atPB2OdWipKREISEh+utf/6quXbtq2LBhmjp1qhYuXOjt0qrM1q1b9dxzz+mtt97SN998ow8//FDr1q3T7NmzvV0aPLBs2TLNnDlTK1asUEhIiLfLqRLFxcW6//77NXPmTEVGRnq7nGpTUlIiPz8/LV26VD169NAdd9yhV155RUuWLKkTs0R79+7VY489punTpyszM1MbNmzQ0aNHNX78+Gp7zroZnb2sRYsWCggIUE5Ojtv2nJwc2Wy2yz72pZde0vPPP69NmzYpOjq6Osu8IpXt8dChQzp69KgGDx7s2lZSUiJJCgwMVHZ2ttq0aVO9RVeCJ/8Nw8LCVK9ePQUEBLi2dezYUQ6HQ0VFRQoKCqrWmivLkx6feuopjRw5UmPHjpUkde7cWWfPntW4ceM0depU+fvX7r+xbDZbua+H1Wqtc7NDy5cv19ixY7Vy5UqfnYn2REFBgb7++mt9++23mjBhgqT/vNcYhqHAwEB98skn6tevn5ervHJhYWG65pprFBwc7NrWsWNHGYahH374Qe3atfNidVduzpw5uvXWWzVp0iRJUnR0tBo1aqRevXrpmWeeUVhYWJU/Z+1+9/JRQUFB6tq1q9LT013bSkpKlJ6eLrvd/puPmzt3rmbPnq0NGzaoW7duNVGqxyrbY4cOHbRr1y5lZWW5lv/6r/9yXc0TERFRk+X/Lk/+G9566606ePCgK+hJ0oEDBxQWFuZzYUjyrMdffvmlTOgpDYBGHfhaRLvd7vZ6SFJaWtplf25ro/fff18PPvig3n//fcXHx3u7nCpltVrLvNeMHz9e7du3V1ZWlnr27OntEqvErbfeqhMnTujMmTOubQcOHJC/v79atmzpxcqqhlfea6rtdG2TW758uWGxWIzU1FRj7969xrhx44ymTZsaDofDMAzDGDlypPHkk0+6xj///PNGUFCQsWrVKuPkyZOupaCgwFst/K7K9ngpX7/KrLL9HTt2zGjSpIkxYcIEIzs721i7dq0REhJiPPPMM95q4XdVtscZM2YYTZo0Md5//33j8OHDxieffGK0adPGuPfee73VwmUVFBQY3377rfHtt98akoxXXnnF+Pbbb43/+7//MwzDMJ588klj5MiRrvGHDx82GjZsaEyaNMnYt2+fMX/+fCMgIMDYsGGDt1r4XZXtcenSpUZgYKAxf/58t/eavLw8b7Xwuyrb46Vqw1Vmle2xoKDAaNmypXHPPfcYe/bsMbZt22a0a9fOGDt2rLdauKzK9rd48WIjMDDQeOutt4xDhw4Zn332mdGtWzejR48e1VYjgagavfHGG8a1115rBAUFGT169DA+//xz177bb7/dSEhIcK23atXKkFRmmTFjRs0XXgmV6fFSvh6IDKPy/e3YscPo2bOnYbFYjOuvv9549tlnjYsXL9Zw1ZVTmR4vXLhgPP3000abNm2M+vXrGxEREcaf//xn4+eff675witgy5Yt5f5clfaUkJBg3H777WUec+ONNxpBQUHG9ddfbyxevLjG666MyvZ4++23X3a8L/Lkv+Ov1YZA5EmP+/btM2JiYowGDRoYLVu2NJKTk41ffvml5ouvAE/6e/31142oqCijQYMGRlhYmDFixAjjhx9+qLYa/QyjDsxzAwAAXAHOIQIAAKZHIAIAAKZHIAIAAKZHIAIAAKZHIAIAAKZHIAIAAKZHIAIAAKZHIAIAAF7z6aefavDgwQoPD5efn5/WrFlTqcc//fTT8vPzK7M0atSoUschEAHAZRw9elR+fn7KysrydikV0qdPH02cONHbZQAVdvbsWXXp0kXz58/36PFPPPGETp486bZERUXpv//7vyt1HAIRgEr505/+VO5fYwMHDvR2aQBqoUGDBumZZ57RH//4x3L3FxYW6oknntA111yjRo0aqWfPntq6datrf+PGjWWz2VxLTk6O9u7dqzFjxlSqjsAraQKAOQ0cOFCLFy9222axWKr1OYuKihQUFFStzwHA90yYMEF79+7V8uXLFR4ertWrV2vgwIHatWuX2rVrV2b8okWLFBkZqV69elXqeZghAlBpFovF7S8ym82mZs2aufb7+flp0aJF+uMf/6iGDRuqXbt2+uijj9yOsXv3bg0aNEiNGzdWaGioRo4cqZ9++sm1v0+fPpowYYImTpyoFi1aKC4uTpL00UcfqV27dqpfv7769u2rJUuWyM/PT3l5eTp79qysVqtWrVrl9lxr1qxRo0aNVFBQUG4/JSUlmjt3rtq2bSuLxaJrr71Wzz77rNuYw4cPq2/fvmrYsKG6dOmijIwM175Tp07pvvvu0zXXXKOGDRuqc+fOev/9990e36dPH/3lL3/R5MmT1bx5c9lsNj399NNuY6ridQPqkmPHjmnx4sVauXKlevXqpTZt2uiJJ57QbbfdVuaPMkk6f/68li5dWunZIYlABKCazJw5U/fee6927typO+64QyNGjNDp06clSXl5eerXr59uuukmff3119qwYYNycnJ07733uh1jyZIlCgoK0vbt27Vw4UIdOXJE99xzj4YMGaLvvvtODz/8sKZOneoa36hRIw0fPrzMG+XixYt1zz33qEmTJuXWmpKSoueff15PPfWU9u7dq2XLlik0NNRtzNSpU/XEE08oKytLkZGRuu+++3Tx4kVJ/3kT7tq1q9atW6fdu3dr3LhxGjlypL788ssy/TRq1EhffPGF5s6dq1mzZiktLa3KXzegrti1a5eKi4sVGRmpxo0bu5Zt27bp0KFDZcavXr1aBQUFSkhIqPyTGQBQCQkJCUZAQIDRqFEjt+XZZ591jZFkTJs2zbV+5swZQ5Lx8ccfG4ZhGLNnzzZiY2Pdjnv8+HFDkpGdnW0YhmHcfvvtxk033eQ2ZsqUKcYNN9zgtm3q1KmGJOPnn382DMMwvvjiCyMgIMA4ceKEYRiGkZOTYwQGBhpbt24ttx+n02lYLBbj7bffLnf/kSNHDEnGokWLXNv27NljSDL27dv3m69TfHy88fjjj7vWb7/9duO2225zG9O9e3djypQprvWqet0ee+yx36wL8GWSjNWrV7vWly9fbgQEBBj79+83vv/+e7fl5MmTZR7fr18/Y8iQIR49N+cQAai0vn37asGCBW7bmjdv7rYeHR3t+nejRo1ktVqVm5srSfruu++0ZcsWNW7cuMyxDx06pMjISElS165d3fZlZ2ere/fubtt69OhRZr1Tp05asmSJnnzySb333ntq1aqVevfuXW4v+/btU2Fhofr373+5lt36CQsLkyTl5uaqQ4cOKi4u1nPPPacVK1bo3//+t4qKilRYWKiGDRv+5jFKj1P6mpQ3xtPXDagrbrrpJhUXFys3N/d3zwk6cuSItmzZUuZj5ooiEAGotEaNGqlt27aXHVOvXj23dT8/P5WUlEiSzpw5o8GDB+uFF14o87jSsFH6PJ4YO3as5s+fryeffFKLFy/Wgw8+KD8/v3LHNmjQoELH/HU/pccq7efFF1/Ua6+9pldffVWdO3dWo0aNNHHiRBUVFf3mMUqPU3qMioyp6OsG1CZnzpzRwYMHXetHjhxRVlaWmjdvrsjISI0YMUKjRo3Syy+/rJtuukk//vij0tPTFR0drfj4eNfj3n33XYWFhWnQoEEe1UEgAlDjbr75Zv3v//6vrrvuOgUGVvxtqH379lq/fr3btq+++qrMuAceeECTJ0/W66+/rr179172fIJ27dqpQYMGSk9P19ixYyvexK9s375dd911lx544AFJ/wlKBw4cUFRUlEfH+y2evm6AL/v666/Vt29f13pycrIkKSEhQampqVq8eLGeeeYZPf744/r3v/+tFi1a6JZbbtGdd97pekxJSYlSU1P1pz/9SQEBAR7VwUnVACqtsLBQDofDbanMlU6JiYk6ffq07rvvPn311Vc6dOiQNm7cqAcffFDFxcW/+biHH35Y+/fv15QpU3TgwAGtWLFCqampkuQ2A9SsWTPdfffdmjRpkmJjY9WyZcvfPGb9+vU1ZcoUTZ48WX/729906NAhff7553rnnXcq3E+7du2UlpamHTt2aN++fXr44YeVk5NT4cdXlKevG+DL+vTpI8MwyiylP9v16tXTzJkzdeTIERUVFenEiRP68MMP1blzZ9cx/P39dfz48TJXh1YGgQhApW3YsEFhYWFuy2233Vbhx4eHh2v79u0qLi5WbGysOnfurIkTJ6pp06by9//tt6XWrVtr1apV+vDDDxUdHa0FCxa4rjK79D5IY8aMUVFRkUaPHv279Tz11FN6/PHHNX36dHXs2FHDhg0rc27P5UybNk0333yz4uLi1KdPH9lsNg0ZMqTCj68oT183AL/P7/87qxsAaqVnn31WCxcu1PHjx922//3vf1dSUpJOnDjBDR0B/C4+hAZQq7z11lvq3r27rrrqKm3fvl0vvviiJkyY4Nr/yy+/6OTJk3r++ef18MMPE4YAVAhzrABqle+//1533XWXoqKiNHv2bD3++ONud3yeO3euOnToIJvNppSUFO8VCqBW4SMzAABgeswQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0yMQAQAA0/t/W61eJpmO6F4AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.hist(data_proc.energy_values, bins=100, alpha=0.5)\n", "plt.xlabel(\"Energy channel\")\n", "plt.ylabel(\"Counts\")\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Average count rate: 450.19 +/- 2.12 count/s \n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGyCAYAAAACgQXWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC6f0lEQVR4nOydd3hTZfvHvyere++WtpRVKFD2KHsPUVFw4GAo+irDheJERfBVXn4KDnAjqIh7l71l71UohUKhe9CR7jbj/P54ck6SNmmTNKVNuD/X1Quac3LyJD055/vc9/e+H47neR4EQRAEQRBOiqSlB0AQBEEQBNGckNghCIIgCMKpIbFDEARBEIRTQ2KHIAiCIAinhsQOQRAEQRBODYkdgiAIgiCcGhI7BEEQBEE4NSR2CIIgCIJwakjsEARBEATh1MhaegCtAa1Wi+zsbHh5eYHjuJYeDkEQBEEQFsDzPMrKyhAeHg6JpIH4Dd+CvPnmmzwAo5/Y2Fhxe1VVFT937lze39+f9/Dw4KdMmcLn5uYaHeP69ev8bbfdxru5ufFBQUH8Cy+8wKtUKqvGkZGRUW8c9EM/9EM/9EM/9OMYPxkZGQ3e51s8stO1a1fs2LFD/F0m0w/pueeew8aNG/HLL7/Ax8cH8+fPx5QpU3DgwAEAgEajwaRJkxAaGoqDBw8iJycHM2bMgFwuxzvvvGPxGLy8vAAAGRkZ8Pb2ttM7IwiCIAiiOSktLUVkZKR4HzdHi4sdmUyG0NDQeo8rlUqsWbMGGzZswKhRowAAa9euRZcuXXD48GEMHDgQ27Ztw4ULF7Bjxw6EhISgZ8+eWLp0KV566SUsXrwYCoXCojEIqStvb28SOwRBEAThYDRmQWlxg/Lly5cRHh6Odu3a4aGHHkJ6ejoA4MSJE1CpVBgzZoy4b+fOnREVFYVDhw4BAA4dOoTu3bsjJCRE3Gf8+PEoLS3F+fPnzb5mTU0NSktLjX4IgiAIgnBOWlTsDBgwAOvWrcOWLVvw6aefIi0tDUOHDkVZWRlyc3OhUCjg6+tr9JyQkBDk5uYCAHJzc42EjrBd2GaOd999Fz4+PuJPZGSkfd8YQRAEQRCthhZNY02cOFH8f3x8PAYMGIDo6Gj8/PPPcHNza7bXfeWVV7BgwQLxdyHnRxAEQRCE89HiaSxDfH190alTJ6SmpiI0NBS1tbUoKSkx2icvL0/0+ISGhiIvL6/edmGbOVxcXER/Dvl0CIIgCMK5aVVip7y8HFeuXEFYWBj69OkDuVyOnTt3ittTUlKQnp6OhIQEAEBCQgLOnTuH/Px8cZ/t27fD29sbcXFxN338BEEQBEG0Plo0jfXCCy/gjjvuQHR0NLKzs/Hmm29CKpXigQcegI+PD2bPno0FCxbA398f3t7eeOqpp5CQkICBAwcCAMaNG4e4uDhMnz4dy5cvR25uLhYtWoR58+bBxcWlJd8aQRAEQRCthBYVO5mZmXjggQdQWFiIoKAgDBkyBIcPH0ZQUBAAYOXKlZBIJJg6dSpqamowfvx4fPLJJ+LzpVIpEhMTMWfOHCQkJMDDwwMzZ87EkiVLWuotEQRBEATRyuB4nudbehAtTWlpKXx8fKBUKsm/QxAEQRAOgqX371bl2SEIgiAIgrA3JHYIgiAIgnBqSOwQBEEQBOHUkNghCIIgCMKpIbFDEARBEESzkZSlREllbYuOgcQOQRAEQRDNAs/zmL/hJPq+vQMHU2+02DhI7BAEQRAE0Syczy7FtcJKSCUc4iN9W2wcJHYI4hbnYm4pFv15Dvml1S09FIJoNZRWq5CUpWzpYTg8iWdzAAAjY4Ph6dJyfYxJ7BDELUxxRS0eWXsM6w+nY93Bay09HIJoNSz85Qxu/3h/i6ZeHB2e55F4NhsAcHuPsBYdC4kdgrhF4XkeC389ixwli+iczy5t4RERROugWqXB7pQCAMC2C3ktPBrH5WymEpnFVXCTSzGqc3CLjoXEDtEoPM/j5+MZOJtZ0tJDIezIt4euY0ey/kJ+IcfxxU5FjRpl1aqWHgbh4By/VoxatRYAsJ8iOzYjRHVGdwmGu6JFl+IksUM0zqErhXjx17OY8fXRFruR1Kg1qFZpWuS1nZHz2Ur8d2MyAGDh+FhwHFBQVoP8Msf17WSXVGHMir0Y+d4eVNaqW3o4hANz4Ipe4KTmlyOP/GxWw/M8Nur8OrfHt2wKCyCxQ1jAnkssnFtSqcK6A9du+utrtDwmfrAPo9/fK862CNupqFHjqQ2nUKvRYkyXYMwd0R4xgR4AgAsOmsqqqtXgP98dR46yGjfKax32fRCtgwO6aA7HGf9OWM7J9BJkK6vhoZBiRGzLprAAEjuEBey7rP+if7nvKpRVNze6k15Uias3KpBVUoX0osqb+trOyJt/n8fVGxUI9XbF/93TAxzHIS6MrRbsiKksnufxwq9nkJSlH3uyA74PonWgrFThnK4K6+5eEQCAA6mFLTkkh0RIYY2NC4GrXNrCoyGxc8vx7aFr+MaKqpuCshrxxhHl747SajW+3p/WTKMzTUpumfj/64UVzfIaPM83y3Gbk3UH0vDlv1etek5SlhK/nsiEhAM+nNYTfh4KAEBcuE7sOGBE5ONdqdh4NgdyKYehHQMBOKZoI1oHh67eAM8DHYI9MbV3GwAssuOI1whzZBRV4sVfzyCjmSaPWi2PTedYCmtSfHizvIa1kNi5hcgrrcYbf53Hm3+fR46yyqLnHNTlruPCvPHyxM4AgK/3p93U1t+X8vRiJ+2GfcWOVstj2heHMPHDfUgvdJyoUVKWEov/uYD/bkpGVollf0sA2KtLSY7uEoIB7QLEx7uG+wBwPLGzJSkHK7ZfAgC8fVc3TOsXBcDx3gfRehCiOEM6BKJPtB8UMglyS6txpaB5JlrWkF9WbZcJ33vbUvDz8Uzxu2Nvjl8vRl5pDbxcZRjWKbBZXsNaSOw4KLtT8vHupmT8e6kANWrLjLsnrxeL/z+dXmLRc4QU1pCOgZjQNRSdQ71QVqPGV/uaHt3heR5JWUqs2H4Jh6+aDxMbip3rDQiSldsvYcIH/1o1W7mQU4rDV4twMbcMUz87iPPZjtFE7Mt9+ojORSuiGPsuM7EzrKPxBUhIY6UVVqCixjHMvSm5ZXjupzMAgEcGt8X9/aLQJcwLAHAxtwxqDfm7COsR/DmD2gfAVS5Fv7Z+APQTv5aC53k88MVh3PbhviYVEqg0Wuy+mA8A2JOSD43W/hErIYU1Li4ULrKWT2EBJHYckqNpRfjPt8fx+b9XMePro+i9ZDvmfn8Cv53IRHkDN6qT6QZiJ6Ok0dfheR77BbHTIRASCYfnxnYCAKw9kIaiCtuiO4XlNVizPw0TP9yH2z/ej492XsaCn06b3d9Q7FxrYFbz/ZF0XMwtw2t/JlkcchYiHQBL2U37/DAOXWnd+fnM4kqxKynAbuyWUFmrxgmd4B3SMchoW5CXC4K8XMDzlh/v871XMOWTAy22wN+6g9dQpdJgcIcAvHZbFwBAdIAH3BVS1Ki1DZ4rBGGK7JIqXL1RAQkHDGzPIp+D2rOJwf7LLSt2bpTX4kpBBSpqNTh8tcjm4xxLK0JpNbtPFFeqLLoXWINGy2PTuVwAraMKS4DEjoORUVSJJ9efgErDo2u4N4K9XFBRq8Gmc7l4/pczeHTdMbPPPWkQzTllQWTnSkE5ckuroZBJ0D/GHwAwLi4EXcO9UVGrwRdW+kWqVRq8/mcSBryzE0sTL+BibhkUUgk4DshWVpss76xVa3HVIHxs7gZWXFGLG+U1AIB/LxXg7zPZFo1JEDsvTojFgBh/lNWoMfPro9iSlNPIM+1HrVrboEity9oD16DR8mKliKVm3CNpRVBpeET4uqFtgHu97V3DLTcp5yqr8d62FJxML8HulHyLx25PMotZBO+unhGQSdmlTCrh0DmURXcu5Fgm2ghCQIjq9Ij0hberHACb6AHAoauFzRIFsZTL+frz+cQ128VO3SaJQpTHXhxJK8SN8hr4uMkxuEPrSGEBJHYciooaNR7/9jiKKmrRLcIbvz45CIdfGY2/5g3GvJHtAQDHrhWZjLjUqDU4l6lP0ZzLUjYa5hdSWP3a+olueo7jsEAX3fnm4DVRYDRG2o0K3P3JQXx3+DrUWh7xbXywdHJXHH1tNDoFs5vTGRMzjLQbFVBrecgk7M6eVVxlsvw8taDc6PeliRcajTiUVqvE1N7t3cPxzaP9MaFrKGo1Wsz9/iQ2n7s5gufFX8+g79vbLcrFK6tU+PFoOgBg1qC2AIwN3A0hzEyHdgwEJyglA8SKLAtSeesOXoNKwy78l/LKG9m7ecjWeZXCfd2MHu8S5rhma6JlEcTO4Pb6m3S3CB94u8pQVq0Wq7Sag32XCxpMo1/J13/PjhtYEqyB53mxkeiErqEAgF12Fju/ncgSj6+QtR6J0XpGQjSIVsvjuZ9O42JuGQI9XfDF9L5wU0ghkXDoEemLheM7IzbECzxvuifE+exS1Gq08PdQwMtVhiqVBil5Dd8k9Sks45THqM7B6NHGB1UqDb61oLJr49kc3PHxfiTnlCLAQ4FvH+2Pv+cPwfSEtvB1VyC+DTPHmrqQCGOMb+MDd4UUWl4/ozfksu6Gm9AuAO2DPHCjvBb/23KxwXEdTC2EWsujXaAHogLc4SqXYvVDvXFvnzbQ8sBXN6HqrEatwaakXFSrtBaVt244ko6KWg1iQ7zwxDAmcK/eqLCo4eJ+A/+VKSytyCqvUeP7I9fF3y83ch41BzzPi8tchPm4Gm2LsyJCRRACPM/jgC6FbRiRkEo4JOhSWs3Vb+dKQTlmfH0Us9YeM5uCv2wgdpJzSq2KBgtczC1DZnEVXGQSLLq9CziOfU9ylfZpmlhYXoN/dH6d+/tH2uWY9oLETgtTUaPGnpR8aBsJj67ccQnbLuRBIZPgixl96s1mAf1NzFRuWYhg9I7yQ482vgAaTmWpNFrRNDy0zs2R4zg8NrQdAOCn4xlmI0QqjRZv/pWEeRtOorxGjf5t/bHpmaEY1slYPMVHsvGcyawvdi7pohaxod6IDmCN70ylslJ1F4K4cG+8OyUeAPDD0QwcTTMf7hVSWIbjkUo4PD26IxtPRkmzd+I9m6kUI1WXGhENNWoN1h5gAuzxYe0Q4u0CX3c5NFpefP/myC+tRkpeGTjOeNZqiBDZaczc++PRdJRVq+Gim7W1RGSntEqNylom8MJ8TEd2qNcOYQ2X88tRUFYDV7kEvaN9jbYJqazmEjv7LhWA55lvMLPYdHXlZYPvmZYHTqVbH93ZoUthDe0YiDZ+7uK9wF6p6B+PZaBWrUV8Gx/00l3XWwskdlqYJf9cwKy1x/Dz8Qyz+xxIvYGPd6UCAJZN6Y7eUX4m9xPFjomeEII5uXe0L3pF+QJo2KR8Kr0EFbUa+HsoxJugIeO7hiLAQ4G80hqzYdCv9qXhm0MsAjBnRHtseHwAQrxd6+3XQxfZOZtZUm/cQmQnNsRT9Jlcu2EisqPLZ3cM9kT/GH9M68dmFa/+cc5ktRrP8/hXJ3aGxxqLrzZ+bojwdYNay4uGXlOcz1Y2eXkFQzF2Mbfhm/Pfp7ORX1aDEG8X3NkjHByn96c0ZioW1vfpFu4j9tapS1sDc6+5En+VRiv2WRJEYUZxJapqb+5SHtm61gn+Hgq4KYyrPTqHejnF8heEfdiSlINub24Vb/Tm2C+m7f3rVRAN0omd49eLm2XZmoMGRRHmUllCqr6drtv58WvWi53tuhTW2LgQABAX59yZXP8arqxSYWniBWw7n2vRsdUaLdYfZtf7mQltTabKWxISOy2IWqPFFt2J1JCy3qjzjtzbpw2m6JpcmWJAjD8UUgmySqpwrU6J9snrJQBYZKenTnE3NDPYb1B+KZHUP2kVMgnu6cPG8oPOQ2JIZa1aLI9++65ueGlCZ9FEWpfYUC8opBKUVKqQUWQ8qxGiHZ1CvcTIjilvixDZ6BjiCQB4ZWIXBHoqkJpfjs/31jdSXykoR1ZJFRQyCQbGBBht4zgOA3SG7CNmqh6Sc0px56oDmPbF4SaZFg3FTkpumdkQNs/z4uf5yOAYMRfeOdRb99yGhVJjKSwAkEg4MSpibgX0jWdzkK2sRqCnC2YPiYG/hwI8zz7Pm4nQJ6puCgsA3BUycfmLZDIp3/J8d/g6ymvU+GjX5Qb3E0rLh5gw1bYL9ECYjytq1VqbREZDaLS8UesNU2nkkspaFJQxf+R9uonc8evWmZRzldU4m6kExwGjOhuLnQOpN+qJuHc3JWPN/jT857sTeH9bSqPZh+0X8pCjrEaAhwK392g9VVgCJHZakOPXi8WlF46mFZk9mYQvwnidocwc7gqZGH7df1lfUp1dUoXc0mpIJRx6tPEVxc6VggqzSz8Iz6+bwjJkWn/WwG3PpYJ6PprvD6ejqKIW0QHuYpTFHC4yKTrr+qOcMVhZvapWIy4P0SnECzGBLLKTVkfIlVWrRP9GhyB2HB93OV6/PQ4AsHp3ar3x7Ulh729AjH+9yAAADNQ13DPX/2fTuRxotDyuFlRg+wXLZj510dSJHBVXqlBgxvC9J6UAl/LK4ekiw4MDosTHLYns8DwvitehjVRHNLRsBM/zYgXerEHRcJVL0TGYicvGUnD2JrtE8OvUT+cClMq6VahRa3C1AaFdWavGsTT2HTubqTTrR1NrtGI5t6kKIo7jxMftvQr6+WylWAoOmP7uCZO5cB9XjNBFok+ll1jVS0owJveK9EWQlwsAiBW9VSoNjhhMvE6lF+PHY/psw8e7UjFvw8kG0/rrdP7NB/pHtZreOoaQ2GlBdibrw6rFlap6FUUA81pcLagAxwH9dNGGhhBmJYbrWQkprLgwb7gppAjwdEGUPxMOpiqgSqtVon+mbj8WQ2ICPTCofQB4HvjZ4ItRrdLgc91Ncd6IDmYjOobEG6SyBFLzy8HzQICHAoGeLmYjO0Jn0yAvF/i4y8XH7+wRjgEx/qhRa7Fss7FZWfDrDO9k+v0JYudMpmnfznaDkPjX+681+v5MIZgMvVxkiNal6C7lmr5w/3maVTjc1zdSLIkFgM7iTd282LiUV458nRehT1vTKVCBrg2YlA+kFuJCTinc5FI8NCAaABOhwmvcTPSVWPUjO4BhZRmJHWfmpV/PYtT7e816846kFaHWQBCYswuczihBeY0afu5yk2l7ABjcoXlMykIKSxAgpqKqgtjpEOKFTsFe8HKVobJWY3FPLEB/zRqjS2EBTMSN1C3SKZSga7Q83vjrPABgau82+L974iGXcticlIt7Pztksvt+ck4pjqQVQSrh8NDAqHrbWwMkdpqRapUG/xy/ioqKCpM/285mAgBc5ezP8O/5jHr77LmQCW1tNToHusDHTX+TM3fMPhHu0KpqcOhKoaj6D13Khra2Gt2CFeJ+XYNcoK2txtHLxuXVlZWV2H0uA6rqKkR7S+Ar14rPqaw0jo5UVVVhSvdAaGurseHgJShLy1BRUYF1ey8iv0iJCF833N07QtzX3JgrKioQrzPKnc1Uorq6GhUVFThzLQ/a2mq085WhoqICwa6AtrYamcVVUOneW3V1Nc7p9ovxkRods7KyEq/rKg4Sz+bgQEoOKioqUFhSisMp7DPpF+Eu7m+YQgr2kCDEjUdNVRUOJGcZHTclowDJOUpIJRxkEg5HruThyKVss+9Nq9VfbGtra8XH9+n+tj3CXNHBTw5tbTUuZBWb3PdYKhtvQpSH+JhGo0GnEE/mT1FW4HpekcnX33H2OnitBv1jAuAik0KlUpkda1tfKXitBhdySsHzPNRqtbht9Y7z0NZW4+7uAVCAHSMmgF2gL+eVGe1r6kel0kcRNRpNg/vW1tY2uG96XjG0tdXwV/BG+2q17HyN8ZVCW1uNpGt5Zo9LOD5C5Pe3E5kmt++7xISJ4HP541SWSc+NELEcERtsMm0PAAnt2ETyfLbSrh41QezMGMgmEDnK6nrtQ4RKrA5BnpBIONG3edxEvx2Nlq8XsS+vUYvNUscZiB0AGNWFiZ1dF/PB8zx+OJqOc1lKeLnK8PLEzri3byQ2PD4Q/h4KnM9m6ft9BpkDgK25CLByc3PR1haHJ3ilUskD4JVKpd2OWVGj4vu+vZ0HYPbHrV1fvv0rG/m3E8/z0S8l8jIXV7P7tu3Wz+j4gYGB5o8b3omPfimRP36tiOd5nnfzDzW7r294jNFx4+LizO4bHR1ttG/fvn3N7itx8+bXH74m7jt8+HCz+7q7u/MXc0r56JcS+bjXN/MTJ97W4OcW/VIif7WgnOd5nr/nnnsa3Le8vJx/+bczfPRLiXyb/hMa3Dc/P18c79y5cxvcN+LJNfz9nx/kn/nhJO/df0qD+yYlJYnHffPNNxvc9+F314v7Ll++vMF9d+/ezT7b5bt4/7FPNrhv0D1v8l/svcLzPM+vXbu2wX2D73qZj34pkc8uqeR//vnnBvd97X8f89EvJfJD/reTT0xMbHDfVatWie9t9+7dDe67fPlycd+jR482uO+bb74p7puUlNTgvi+88EKj313Ccbjz43189EuJfI+3tvK1ak297aPf38NHv5TI/3Mmi094Zwcf/VIi/+epTKN9Tl4v4qNfSuRjXk7kL+eVNvh6/d7erru2Ftpl/DUqDd950WY++qVE/mJOKT9s+S4++qVEfv/lAqP9pq85wke/lMhvOHKd53me/3jnJT76pUR+7vcn6h1zwU+n2bb1J/jU/DKe53l+49lsPvqlRH7E/+3mtVqt0f7l1Sq+46ub+OiXEvmjaYV8/OKtfPRLifza/VeN9ksvrODHrdjLR7+UyEe/lMgv+uMcX1Gj4osravjYRez5R67a53OxBkvv3xTZaSbcFTL0b9t42mlAO3+M7sKUtraB9Ku3m8zi1/Z0Yfvuv8xMZ6oG8rrVKo1RRMNeDUIlHCcamC2hQ7An3BVSVNRqLCr3tmYpgOfHxcLLRYbiStP+JFsZ0yUEjw6JsesxG1r7yxyCSbkxGjInGxKqq5izJAUU6s0iO5nFVahphioVgmgIte6CVVKpquevyy6pQmp+OSQcMLRDEO7py7yDhqksnuexfEsKAJay6aBrcGoOfbrdPs0FT2eUoEqlQYCHAp1CPM2mX4WGgoJHrk80u7ecuFZsdP0+nVGC306yKNfGczkYt/JfvPzbWfyqi3yN6RJcr0rKw0WGAe3Y8Z787gSUVSp0DvXCw7pIk0Ckvzv+mDcIMxPY498dvo6JH+7D2xuTUa3SokuYt7iOWGvE8jsoYTWPD2uHf577FXIph+3PDUeIQeXIzK+P4Nh1JUZ3DkHPSF8opBJEzF+PTU8PQUwQO6Hzy6oxfPkecByw9dWxRse+du2a2df9+Xgm3tqciv2pBRjUIQBhsz9BoKcC/y4cKZ7oNWoN+r+9E7VaLdKLKkU/zEPvrMcX/15FVIAb/p4/xMhoVvdL8u+//0Kr1eL6jQpM+HAfOA4I8FTgRlktXruts9FzN2/ebJTOqYtUwqFbuA+OXivCf5asxsa/wjDyvd3IVdZgw+P90SuKfRmf+uEkdqWW4vqNCiAW+O6775Dd4xFkFFVh3SP9jFbyBgB3d3d4cByeGt0B/62Yj05TnoGrXILMomp8NK0nxhqYvt3d9UsorFixAvNfehPjP9gHuZTDkVfHwE0hhbJShcH/2wWtVI6xcSGIDvDAmOlP49jgB/HE8Bg8Oya23ntzc9OHdV999VUsXLgQV/LLcfvH++Eik+Doa2OQUVSJ2z/ej3yFK7RaHhIJh2eeeQZz587Fiu0p+PLfNNzTJwJL7+ouHsvVlZ1PsaFe8OwxAQ9Pn4l3pnQ3eu0jVwsxa+0xBPl4imbmhx56CPfee6/Zv8WrfyXjr7N5uJBdikF9R6LrK3+itFqNHpE++GJGXyPPkEKhwNple1BYUYuYnoNRXm7eu6NQ6Evehw4disuZBXh4zRFkl1Sjjb8r8pQ1UGl4RPq7Ycw9/cR9e/fubXRcrZZHryXbUavRYseCYWgb7CNu69Kli7jv7HXHcPBKIZZM7op7dTc6uVw/dsLxMayE3HQuB0MNPIZCqqVHpC983OW4t08bfLzrMg6kFiK9sBJRAe7Yn3oDh64WQiGV4FldZ/iG6B7hix3J+Ubd6JuCUAGW0D4AHMchLswbm5NyjcrPy2vUyNJ51DroxE7PSF/IJBxyS6uRVVKFNn7s2rVc10h1jC41tSM538hoPDbOdJHLyNhg7Lt8A4W69NnSu7qZ9Fq6K2R4a3I3jI0LxYu/nsH1wkpxgjZrUHSrKzc3hCI7zUjPSF8MjA2HRuqCn87kw8PDAx4eHlBzCpzKqQYnU2BMlxC4yqXoGekLicIVSfk14n5JeTWQKFzRLToYoQHGs3dhH1M/o7uxC/up9BLsu1QAidwV/TqEwdPTU9zH38cb3doGQyJ3FZsLZhRV4tvjuZAoXPHGXb3h7+NtdFxDMQCwm7iHhwfiooMxNC4CnNwVRTUSBPv7YObwzib3NfcDAN11s6aLBdVQSxTIr+IgUbgiPiZU3K9jBItOiKX1UjmyK3jdfiH1jit8+WYNikFMiA+KaiTILgfkrm4YFR9lcl8AcHFxQefIIEQE+UIjdUFKYS08PDxwLKsCvMzFqMnhY8M7QaJwxa9nbkAid603BolE/zVTKBTw8PDA+QL2t+3dPhR+Pl6IiwqCq7s7qjWsb43Rvvls3wGdwo2OK5UyMdklzAucVIarJep6r30iuxIShSuGxupndHK5vMG/Rbc2TFhuSsrFrG9OoFwrQ/+OYdgwZzjCAnyN9pXL5eIF+GphVYPHNRQaxVVq/OfHJORWAh3CA/DXM2Pw29OjEBnsh6xyHvd/dVzs0iyVSo2OUw051FIFpC6uiAkLMBJREolE3C8+JgQShavR52K4L+H4GIqdrefzjKqT/tUVaQzTCaBIf3exgOOXExngeR7/t5VFdR4aGIUIE41a6yJGduy0bMTBVOOOzV0j6ldDClGdQE8X+Lqz89dNIUXXCDYWoRR+/+UbOHiFCbe3JnfDVzP74dcnE8RoS4SvG3rreqzVRShBB4ApvSPQr5GsxJCOgdjy3DBM1bVCCfJyweSeEZa/8RaAxE4z8x9dp+ENh9NRVs3SKHsu5UOj5dEx2BNRuiocIYxoWP4nhGUH1olWNEZUgDui/N2h1vJiUz9TjQjrNhd8Z1MyatVaDO4QIDadspQH+usd+E8MayeupWUNwoXkTKZSXIIg3MfVKJLQtk4X5SsFrGLL112OADPN8gDWF+i1SXHi732i/IyOawrjfjvsbyEsomf4+YzrGoo2fm4orlThj1NZFr1XoXpEOL5MKkEHXUTPcK0rtUaLMxnswmqumaSQxrqUV2Z08ddqebECo6GquroIyy0k55SirFqN/jH++ObR/vAy83mJFVn5llWGKKtUmLHmKK4WVCDcxxXfPTYAQV4u6Bnpi41PD8GYLsGo1Wjx2h9JJhtfCtUgwV4ukDdQ6RdH5edOj+H5XlRRK14/NVperJoa1kmfvr1PF+H75XgmNp7LwdlMJdwVUswb2cGi1+umExhXCsptWq7BkMpaNU5lMKEySLccRVyYcHz9EjCpdVJYAn2jdSbl60U64caiOobCrW9bf/z8RAL+mjcYPz+ZYLYytm2gB0bEBiHK3x2vTOxi0fi9XeV4/74e2PT0UPwxd5BN1/ybCYmdZmZU52C0D/JAWY0aP+nCiULHYcGrA0BcVfyoCbEzwIKS87oI/gzBld8nuv6N0rC54OGrhdiclAsJB7x+e5zV4chxcaHoFOKJdoEeNpceCq3Lk3NKxfLLjiHGOXShRFsInRpeCBob85guwWLfoNFdghvcV0Dfb6cItWot9ur684wx+NtJJZy4KOfXB9LMNgY0RPg7G7YTiA0Vyrj1oiElrwxVKg28XGRoH2R8sROI8neHm5x1Pjb0Mu1IzhN784ztYrl4NSy9HdjOH+se6QcPF/MZ7066Ro6pJsrPM4oq8eZfSZiz/gSmfXEIEz74F8OW78aFnFIEerrg+8cHGs2ofd0V+HJGX/HvZKrapLEeOwJdDMryG2uIRjgmgmdHOAeFBqznspQoqVTBy1UmXlcAYFzXEPi6y5FbWo2XfzsHAHhsSAwCPV0ser0gLxeE+7iC54HzTYzuHL9WDJWGR4Svm9gKJMTbBf4eCmi0vDjpuVynYaqAKHauFWPr+VycMSPcOI6tn9hY5GrdI/2xd+EIsQTeUuLCvcU0WmuGxE4zI5FweFwX3fl6fxqqVRqxoZ3hDbdPtB9kEg5ZJVXIKKpEQVkNruj66/S3QewYNo+TSzlxRmJIr0j2ZbmQU4rFf7O+Cg8OiLLY8GqIQibB5meGYdtzw+CusM0KFh3gDh83OWrVWvxzhi0mJwgAAaEzbkZRJdQarbheTGPGQoB96Vc92Bsr7uuBRwZbZiwWxM7pjBLsSclHeY0aQV4uRhdQgHU19VBIkZpfLobPzZFVUoWskipIDUpIAf17NeydIaQYe0b5mi2JlUg4dBKeq+u3w/O8uMTIjIRoo/5DjeHnocBTozpgWr9IrJ3Vv9G/Z8cGIjtLEi/gm0PXsTkpF4evFuFibhmUVSoEeCjw3ez+4t/TEI7Tfy6m+ogIkR1zPXYE2gV5QCGToLxGLaYGCedCiOzcER8OANialAuNlsc+XR+twe0DjaIZLjIp7u7F0i1CX53HhrWz6jW7N7BwsTUIJeeCXwdg537XOgvZpuq+Vx3qRHaEnlkpeWViHzFrhJspWrPnpqmQ2LkJ3NUrAoGeLshWVuOtfy5AWaWCn7vc6EbnrpCJguRIWhGOpLEvQpdQbzFPaw2D2gdCuDfGhfuYDDFG+rshwEMBlYbHxdwyeLvKsGBsfYOtpUglnEUNBM3BcZyYyjqu6yzcqU5kJ8TLFS4yCdRaHlm6agugfojXHD5uckzp3UZcbqExogPcEertilqNFu9vuwSARYjqCg9vVznu0l1EdyU3vAbPMV1Up1u4t1HEJDakfmRHaAjZy0wKS6CLKJTYBXJPSgHOZbGZnrBoqzU8Py4Wy6bGm+wuXRfhb5RRVGVUSXejvEZsVPbShM746IFe+G52fyQ+NQT7XhopRl5MIZipU0yIHaGhYGORHblUIn6mlMpyTgSxM6RjIHzd5SisqMWRtEL8K3SA71S/AvF+g47uc0d0aDSdXRehJ1jTxQ6bFAnNCgXqVmSJPXbqXOOCvVwRHeAOnmceRl8bhNutBImdm4CrXIpZg1i5nrCO1MjYYEjr3DAF387RtEJ9Cqud9VEdgC2X0F33pTRnSuM4TkxlAcAzYzrBvwHfy82ge50IVGwdsSORcAa+nUpxAdC6FwJ7wXGc+DcQFiU152dK0OXdjzeweCgAHNWlZupG7ITIztWCCnEl9NO6yE4vM3/Dus9NzmHrawnrAD08MLrZ/6b+HgrRL2W4+vqfp7Kg1vLoEemLOSPa484e4RjaMQjdInwajRYZpvTqrj2WrRTSWA1HdgBm3gaok7KzIqSxXGRSsVneL8czcVL3vRlmwqvWOdQbswa1xZguIZieEF1ve2MI16imVGQpK1VI0okloVmhQJxBZKdapUGGrnFiRxPRa0N7wtwR7a0WbrcSJHZuEg8NiIabQXRltAkPhWiGTSsS12mx1pxsyJzh7REb4mVkHq5Lb92XpV2QB2bY8MW3N/EG6SGOMy1iBN/O5bwysSqrbj7bnhj+DdzkUgxqb7pfjXDhSc4pRUUD5kXRr1On4iHMxxVerjKotTyu3ihHcUUtrupWH+9ZJ21WF3FB0LxSHEgtxKn0ErjIJHhsqH37AJlD+PyFZSN4nhd7e1jTb0kgOsADrnIJatTaesuD5IhLRTRePSN8Lpfz6/uJCMdHo2tnIZNyuK07W3zyj1NZ0Gh5xAR6INLftJdk8Z1d8dXMvjaZagWxc/VGBUqrbevddSStEFqeXXdD64j2rgYFAqn55dDyLCId6Fl/0iLcM0K9XTEjoa1NY7lVILFzk/DzUIjhU5mEM6oQEOjb1h8SjplvU/PLwXG2mZMFJnQLxdbnhtVLBRkyPSEaTwxrh88f7tNgZcvNokekPrIT7e9uMo3SVufz2HupABotD08XmdgIrzkwFDvDOgWavUCG+bgh3McVWt70mmMAUFheI0Y/6oodjuPESFZKbplYidQu0AN+jURnhLRPRlEVluuqMh7oH4Vgr+b7XAwRzjEh0nY+uxQXc8ugkElwp85PYQ1SCSces24qS1j01RKx4+XKIkimlgggHB8h6ieVcBjUPhDervqI4TALm2hai5+HApH+7NxLsjGVJfh1BpuYOMUEesJVLkFlrUZcvNNcAcZdvSIwb2R7fD69T6uvhmppWv7udgvx+LB2iPR3wwP9o0yW8Xq7ysUQJsBmpbb4dazB21WOV27rUq/qqaUI9XYVqwHMiTQhjSWk+tpbUInVFNoGuCNE1yl4TCNVTUKk7ISZVNaxa4IXydOkgOlk4FU5ZaFfB2AXYGGMZzOVUEgleHJ4+0afZy+E80cwjP+i61I7Li7EKnO0IYLwu1inFD+vVCd2LEhjCaliNVVjOSWC2JFJOChkEowzaBI61Ip2C9YSH+ELwLZUVnFFrdiiwlRXc6mEQ6wuIvm3rlDDXOTaRSbFwvGd0cPAjkCYhsTOTSTC1w37XhyFpXd1M7tP/7b6KEJTojqOCsdx6KEzKdetxBJoq0tjqTTsQmepObkpY3r7ru54ZHBb3Nmz4SiFUA56It202Nl7iRl2E8ykJw2Nuad0kZ3G/Dr65+qF8r1929QLjzcnwt/gUl4ZatQa/KW7SAudi20h1oRJOb+sBlqeVRhaUnUiiJ26vh/CORBErEQ32bmtOxM7cikneuiag+5NaC64csclKKtU6BLmbXbyJJiUrxawFK4l1aZEw9ByEa2MAe388fWBNABN8+s4Ms+M7gRXuRTTB5r2EEXXKVdubrEDMFOyJY0WhTVrTl4vFpd9EOB5HjuTmdgZZeYi18mgeqismvl+LBc7Xth7qQAyCYc5I25eVAfQjzuzuAr/nMlBSaUKod6uYsdaW9D7kPRiRyg7D/F2NVuKb4hM172aIjvOiRjZkbJzYVjHIMwa1Bbtgzwa7A3VVOJtNClfyivD90dYkcobt8fVK1IR6BpuXKnYXAUYtxIU2Wll9G/rD4VUAoVUcktGdgA2a1r1YG8Em/HhhHm7GpWON6c52Vo6h3nBTS5FabUaqQXGptikrFLkl9XAQyHFQDNVdkLqJltZjbIaNdwV0noVaeYQWr7PGtT2pjf58vdQiAbKldtZif6U3hFmL+aWIER2rhVWoKqWeW6EhoLhjZSdC1Bkx7nR8HrPDsA6kS++syumN7NZV1iqIb2oEiWVtRY9h+d5LE28AI2Wx4SuoQ1GnuLqiJ2bMaFzdkjstDL8PBRY90g/rHukX6Om1FsViYRDtEGVRYeg1hPilUslYjl/Xd+OYDYc2jHIaJFUQ/w8FAg26GAa38bH4t5FA9oF4Mwb4/DaJMvavdsboTRWWLTQliosQ4K8XBDgoQDP643PQmQnrJGGggIyEjtOi1bLQ2hWLpPc3FuZj5tcbIhp2G9Hq+Xx7qZkzPv+JNILjRtZ7kzOx77LN6CQSvDqbQ1/RzuHekGwIXoopBa1WSAahsROK2RQh0AMakL4/1ZAWITTVS5BhJ9ls/ybRR8zJuWdF5nYGdXIUhWGXiVLzMmG+LjLW6wLaieDCFufaD+0M7O8hTXU7Spt6VIRAlIpiR1nxTA1KW2Bc14oQT9rkMr6aNdlfP7vVWw8l4MJH/6L749cB8/zqFVr8d9NyQCA2UNjxDURzeGukKGdTkx1aOYCjFsFEjuEQxITyC4W7YM8m5QqaQ5MiZ1cZTWSskrBccYrDJvCMG3Vy4GqLDoYjLupUR2BuiZloXtyhJWRHfLsOB9agzXoBFF7MxG6vQu+nU3ncvDBDtbMs2OwJyprNXjtjyTMXHsMK7ZfQtqNCgR5uVi86GhcODs+mZPtA4kdwiERltbo2QrFgLAMSNqNChSW1wDQL/7aM9K30SqiTk2I7LQkQiWZq1yCSfFhdj2msAxGjtLKyA4nRHa0dhkP0XowFLCyFpjwiJ2Us5RIylJiwc+nAQCPDo7B1meHYdGkLlDIJPj3UgE+23sFALBwfCw8LTRO3983EqHerpjaO6JZxn+rQdVYhENye3w4grxcjDoutxZ83OXoGOyJy/nlOJlegrFxIdip8+s01qcHYGJJwrEKJ2tXIG5J+kb74elRHRAb6m23tvVCvxEhsmOtZ4f67DgvGo3+byppgTRP1wgfcBzzqD267hiqVVoM6xSEV2/rDImEw2ND22FEbBCe//kMzmQqEd/GB/f0tjziOaRjIA6/OroZ38GtBYkdwiEROqa2VvpE++FyfjlOXC/GkA6B2J/KFv0b3YhfB2A5+t/mDDJbjdZa4TgOC8bZvpCsKTqFeILjgBvltcguqcKNclb5Ymk1low8O06L2iBa1xKRHU8XGdoHeSI1vxz5ZTVoF+SBjx/oZVRQ0CHYC7/NGYR9qTfQO9LPonYJRPNAaSyCaAb0vp0iHEi9gRq1FhG+bhaXkfeK8kOEBcshODvuChmidJV3ey+xlaxd5RL4WtiVWSr02dGQ2HE2hLJzjkOLiQih3463qwxfzegLH7f656VMKsHI2GCbO4kT9oEiOwTRDAhi50ymEpuTcgGwqA5VVVhPbIgXrhdWYrfO9xTu42bx5yjM+A3NrIRzYLhUREsxY1Bb5Cir8cyYjnapPiSaD4rsEEQzEBPoAT93OWrVWvx1mq2DY2qle6JxBJOykAq01K8DkGfHmRGidS3h1xHoGemLH/4z8Jbtdu9IkNghiGaA4zgxuqPW8g12TSYaprNunaBKXRdlSyuxAOqg7MwI0bqWjOwQjgOJHYJoJoR1soCGuyYTDVN3QdhwK7xMYmRHQ6XnzoYQrWttfbaI1gmJHYJoJoTIDtB412TCPG0DPOBisBZauBWt82m5COdFQ2KHsAISOwTRTMS38YG3qwyuckmjXZMJ80glnNFir2G2RHZI7DgdgmdHepPXxSIcE6rGIohmwlUuxS9PDoJKo220azLRMLEh3kjKYl2UrYvssBshRXacD/LsENZAYocgmpG6fhPCNjobfI62RHY0VHrudJBnh7AGiv8RBNHqEUSjl6vM4rWFAP2sn+cBLUV3nAphvTMSO4QlUGSHIIhWT/8YfwxqH4C+BqZvSzDsrKvW8lDQjdFpEDw7lMYiLKHVRHaWLVsGjuPw7LPPio9VV1dj3rx5CAgIgKenJ6ZOnYq8vDyj56Wnp2PSpElwd3dHcHAwFi5cCLVafZNHTxBEc+Iql2LD4wOtXnvL8EZIvh3nQkhNUmSHsIRWIXaOHTuGzz//HPHx8UaPP/fcc/jnn3/wyy+/YO/evcjOzsaUKVPE7RqNBpMmTUJtbS0OHjyIb775BuvWrcMbb7xxs98CQRCtEKlRZId67TgTVHpOWEOLi53y8nI89NBD+PLLL+Hnpw9RK5VKrFmzBitWrMCoUaPQp08frF27FgcPHsThw4cBANu2bcOFCxewfv169OzZExMnTsTSpUuxevVq1NbWttRbIgiilUCRHeeFDMqENbS42Jk3bx4mTZqEMWPGGD1+4sQJqFQqo8c7d+6MqKgoHDp0CABw6NAhdO/eHSEh+jWHxo8fj9LSUpw/f97sa9bU1KC0tNTohyAI50Nax7NDOA/aVrAQKOE4tKhB+ccff8TJkydx7Nixettyc3OhUCjg6+tr9HhISAhyc3PFfQyFjrBd2GaOd999F2+99VYTR08QRGuH4zhIJRw0Wp6qsZwMiuwQ1tBikZ2MjAw888wz+P777+HqanmTMHvwyiuvQKlUij8ZGRk39fUJgrh5UBdl54Q8O4Q1tJjYOXHiBPLz89G7d2/IZDLIZDLs3bsXH330EWQyGUJCQlBbW4uSkhKj5+Xl5SE0NBQAEBoaWq86S/hd2McULi4u8Pb2NvohCMI5kXK0PpYzQpEdwhpaTOyMHj0a586dw+nTp8Wfvn374qGHHhL/L5fLsXPnTvE5KSkpSE9PR0JCAgAgISEB586dQ35+vrjP9u3b4e3tjbi4uJv+ngiCaH3IKLLjlOg9Oy1uPSUcgBbz7Hh5eaFbt25Gj3l4eCAgIEB8fPbs2ViwYAH8/f3h7e2Np556CgkJCRg4cCAAYNy4cYiLi8P06dOxfPly5ObmYtGiRZg3bx5cXGgtIoIgAKlUiOxQ6bkzQZEdwhpadQfllStXQiKRYOrUqaipqcH48ePxySefiNulUikSExMxZ84cJCQkwMPDAzNnzsSSJUtacNQEQbQmKLLjnNByEYQ1tCqxs2fPHqPfXV1dsXr1aqxevdrsc6Kjo7Fp06ZmHhlBEI6KaFDWkNhxJiiyQ1gDJTsJgnBqBE+HllY+dyqozw5hDSR2CIJwaqj03DmhyA5hDSR2CIJwaoSbIZWeOxfUZ4ewBhI7BEE4NeTZcU5I7BDWQGKHIAinRkaRHadETZ4dwgpI7BAE4dToPTvUZ8eZoMgOYQ0kdgiCcGoosuOckEGZsAYSOwRBODVkUHZOaLkIwhroLCEIwqkRboYkdpwLiuwQ1kBihyAIp0aY+FOfHeeClosgrIHEDkEQTg1FdpwTiuwQ1kBihyAIp4Y6KDsntFwEYQ0kdgiCcGr01VhUeu5MUGSHsAYSOwRBODUU2XFOxD47HIkdonFI7BAE4dTIpOxmqCWx41SIYkdKYodoHBI7BEE4NVKdQZkiO86Fhjw7hBWQ2CGIW5x9mftw7z/3IqUopaWH0iwIE3+qxnIuBPEqoTQWYQEkdgjiFkaj1eCdI+/gYtFF/H7595YeTrNAkR3nhCI7hDWQ2CGIW5id6TuRWZ4JALhYdLGFR9M80NpYzones0O3MaJx6CwhiFsUnuex7vw68feLRReh5Z2vPFswsKo1JHacCTVFdggrILHjoGy/vh3/O/o/1GpqW3oohINyMv8kzt04B4VEAYVEgUp1JTLLMlt6WHaH+uw4J+JyEeTZISyAxI4DotKqsPjgYqxPXo/NaZtbejiEg7IuaR0AYHKHyejk1wkAkFyU3IIjah7EVc95iuw4E9RUkLAGEjsOyKm8UyitLQWAVi92LhVfgkqjaulhEHW4qryKPZl7wIHDjLgZ6BzQGYBj+XbSS9Ox4vgKZJRlNLifjJoKOiVanXiVUZ8dwgJI7DgguzN2i/8/nHMYRdVFLTga82y9thVT/56KFSdW3JTXW5u0FjM2z4CyRnlTXs+R+fb8twCAkZEj0danLbr4dwHgWJGddefXYe35tXh408O4UHjB7H4SIbJDnh2nQvBgNVdk5/vk73HfP/chpzynWY5/K3Ew6yAuF19u0TGQ2HEweJ4XxY6L1AUaXoNt17a18KhMs+P6DgDAprRNN8X4uj55PU7ln8K/mf82+2s5MjeqbuDvK38DAB7p9ggAoLO/LrJT6DiRnazyLABAUXURHt36KI7kHDG5H0V2nJPmXC5CrVXjszOfIbkoGT+l/GT3499K1Gpq8cbBNzDl7yktem0mseNgXC65jKzyLCgkCszuPhtA60xl8TyPo7lHAbCb0fkb55v19SpVlcivzAcAnLtxrllfy9HZkLwBKq0KPYJ6oGdwTwBAR7+OkHASFFYXoqCyoGUHaCG5FbkAgDCPMFSoKjBnxxxsvba13n5Cnx0qPXcuBA9Wc0R2TuSdQElNCQBgy7Ut4MnvZTN/XfkLeZV5CHILwoCwAS02DhI7DsaejD0AgIHhA3F3h7vBgcPJ/JPihb85qFRV4rX9r2HliZUory236DmXSy4bpdf+zWqaom/sYmPo20i6kdSk13JmKlWV4kz1ka6PiI+7ydwQ4x0DwPJUFs/zLXYT4HlePOc/HvUxxkaPhUqrwsK9C/Fzys9G+1JkxzkRmwo2g2dHiEoDLIJ49sZZu79GS5JSlIK5O+Y2e9d0lVaFNefWAGBRZBepS7O+XkOQ2HEwdqezFNbIyJEI9QhF75DeAIAtaVua7TX/SP0Df1/5G18nfY07/7wTm65uavQmdzSHRXXkEjkANCl8qdFqMG3jNEz5ewpUWtNm52ul18T/Xyy6aJUpOqUopdkjT62Frde2orS2FJFekRgROcJom7Um5aWHlyLhh4RGDcLNQbmqHJXqSgBAlHcU/m/Y/+G+TveBB493jrwjRvkAg2osKj13KgTPjr2Xi9DyWuxM3wkACPcIB9A6o+dN4YOTH2Bf1j58eubTZn2djVc3Iqs8C/6u/rin0z3N+lqNQWLHgcivzEdSIYtaCDeq22JuA8B8Mc0Bz/PiTNlT7omCqgK8tO8lPLbtMVwtuWr2eUdymX/igc4PAAAuFF6wOT1yqfgSLhRewOXiy7iuvG5yn+ul+sdVWhUuFV+y6NgFlQV4eNPDeGjTQw5ViWQrwpIQUztOhVQiNdommJQt+RwqVZX4K/UvVKgqsCt9l/0H2ghCVMfHxQduMjdIJVIsGrgIXfy7QMNrcCz3mLivvs/OTR8m0Yzol4uw723sbMFZFFQVwFPuiRf7vQiATRI0Wo1dX6elyKvIw8HsgwCA/Vn7UamqbJbX0Wg1+OrcVwCAWV1nwU3m1iyvYykkdhwIIYUVHxiPQLdAAMDY6LGQcTIkFyXjmvKa3V/zeN5xXFVehZvMDRunbMT8nvPhInXB0dyjmPrPVKObioBaq8bx3OMAgNva3YauAV0BsC+WrWMQSFWmmtzHUOwAlvt2vrvwHao11dDwGrx18C2nuaCZ4mrJVZwuOA0pJ8XkDpPrbRdMysmFjaexjuQcQa2WNbQ8lX/KvgO1AEHshLiHiI9xHIf+of0BGJ8zFNlxTprLsyOksIZHDsewNsPg4+KDG1U3cCyv/rXOEUm8migWjNRoaprNNLzl2hZcL70OXxdf3B97f7O8hjWQ2HEghCqskVEjxcf8XP0wMHwggMZDrRqtBvN3zsf8nfMtro768eKPAIA72t0Bf1d/PNHjCfw5+U/0D+0PtVYt5mMNSS5MRrmqHF4KL3T264xhbYYBsD2VJQgnALhScsXkPoLYifFhvhNLfDvKGqXoX5FxMiQVJuHHlB9tGqMjIER1hrUZJoplQwSxk1meibLasgaPZejBOpV/6qZ7d/Iq8wAAoR6hRo/3De0LwPickZJnxylpDs8Oz/PYkc7EzpioMZBL5RgbPRaAc6SyeJ7Hn6l/AgAivSIBANuu27+aV8tr8cXZLwAA0+Omw13ubvfXsBYSO82IskaJf6780+A+G69uxN1/3d1oD4JKVaVYWjuizQijbYaprIZuOrszdmNv5l7szdxrUaqioLJATFHcF3uf+HgbrzZYPGgxAOBg9sF6SwwIKax+If0glUgxvM1wcV9TXpr8ynyzoVQtr8XJ/JPi7+bETnppOgBgUswkAJaJnZ9SfkKluhIdfDvg5f4vAwA+OvlRs5q9WwqVRoV/rrJzcUrHKSb38XHxQZhHGICGU1k8zxsJ16LqIiPP1M1A+BuFuhuLnd4hvcGBw7XSa6JvhxYCdU7UukidPT07F4suIqs8C65SVwwKHwRAf33dfn27wy/Pc6bgDK6VXoObzA1LBi0BAOzL3Gf3VNb269txVXkVXgov0crQ0pDYaSaUNUpM+WsKXt3/qtn+HznlOXjr0FtILUkVZ93mOJB9ACqtCpFekWjv295o28jIkXCRuuBa6bUGb1LfXfhO/P/hnMONvoffLv8GNa9Gr+BeiPWPNdoW6RWJhLAE8ODrjV0wJ/cPYymFLgFdEOAagEp1pVF6AWAz8Am/TcC8nfNMjuFKyRWxBFT4vS7KGiWKa4oBsLQZwDoEV6gqzL63KnUV1l9YDwCY3X027o29Fz2CeqBSXYl3j7xr9nmOyp7MPSiqLkKQWxCGRAwxu5/Yb6eB8yilOAX5lflwk7mhe2B3ADc/lSWmsTxCjB73VniL70GI7gil5xTZcS6EJpGNLQRaUl2C75O/bzRaCUCM6gyJGCJGI3oH90awWzDKastwIOtAE0fdsvx15S8AzP7QJ6QP2ni2QbWmGvuy9tntNQyjOg91eQheCi+7HbspkNhpJnxcfDAskqVvFh1YJC7vIMDzPP575L+oUlcBYH0dGkLw64yMHAmuzkzGU+EpporMhVrP3zhvFCE5nN2w2FFr1fjl0i8AYDbfem/svQBYtZZQJVWrqRVvfAPDWHpNwkkwtM1QAMaprLLaMry6/1WotCoczzuOjNL6VT3C5yKkp9JL0+tFh4SoTrBbMCK9IhHmEQYefINddX+//DuKa4oR4RmBCW0nQMJJ8GbCm5BxMuzK2CVWYzgLgiCd3GEyZBKZ2f0sMSkLf8MBoQPEv3Fj56+18Dzf4GzTXBoLAPqF9gOg9+1QZMc5sdSz89nZz7Ds6DK8efDNRo+58zr73o+OHi0+JpVIMT5mPADHTmVVqavEqt27OtwFjuMwti1L0dmzMe3O9J24VHwJ7jJ3PNzlYbsdt6mQ2GlGFvZdiEivSORW5NaLFmy/vh17M/dCxrEbT0pxitkeNmqtGnsz9wJAvXJhASHU+nvq77hRdaPe9u+SWVSnW0A3AGwmXqOpMTv2vZl7kV+ZD39XfzFnXZcRkSMQ4BqAG1U3sDeDje9MwRlUa6oR4BqAdj7txH0FMWY4g1h2dBlyKvSt2Lder98QTrhh3RZzGzzlnlDz6nopE+H3aJ9o9h4D2Xs0Z1JWaVX45vw3AFivGeHm39GvI2Z1mwUAeOfIOxb3FGrt5FbkijPSuzvc3eC+okm5gV47gtgZ2mYoegX3AmDfyE5JdQlmbpmJoT8OrWc8FzCXxgKAviHMtyOY58mz45xY6tkRzoPt17fjTMEZs/tdVV7FFeUVyCQyMfUuMLHtRAAsQtpc1UuWkFqcipN5Jxvf0QQ703eiXFWOCM8I9AnpAwAYH81E3L6sfeLEuylkl2djySGWHnuoy0PwcfFp8jHtBYmdZsRd7o53hrwDCSdB4tVEbLnGVHVZbRmWHV0GAHgs/jG08WwDLa/F6YLTJo9zOv80lDVK+Lj4iDeXugyPHI7O/p2hrFHinSPvGG3Lq8jD1jQmJBYNXIRAt0BUa6pxJt/8F/+ni8y4e3eHu6GQKkzuI5fIcXdHdvMUokBC1+T+Yf2NIlAJYQmQSWS4Xnod10uvY8f1Hfj7yt+QcBJMbs8qg7Zf3250fJ7nxYhB35C+aOfLxNMVpXEqS7ghRnlFAYCYWjHn29mcthk5FTkIcA3AXR3vMtr2RPwTaOPZBvmV+fjy3Jcmnw8wwbTs6DJx2YWWYG3SWnx25rNGzcF/pv4JHjz6hfZDlHdUg/t2CWCRnaslV02K4eLqYpwtYA3WhrUZhp7BPcGBQ0ZZhl06L+eU52DGlhk4lX8KtdpaI6OxAM/zYmSnbhoLMPbtFFQWiDdDLYkdp0JtwXIRZbVlRn7I94+/b/b7IkR1BoYNrJd66RbYDZFekahSV4lRdntzruCcGKU2hUarwWPbHsMjWx+xqRngX6kshTW5/WRIOHbrjwuIQ4RnBKrUVSarZa1Z5qdWU4vn9zyPkpoSxAXE4YkeT1g9xubEarGTnJyMN998E6NGjUL79u0RFhaG+Ph4zJw5Exs2bEBNjflowa1Iz+CeeKz7YwCApYeWIq8iDx+e/BAFVQVo690Wj3V/TFTZ5lIBQhXWsIhhZlMQcokcSwYtgYyTYfv17UZt839K+QlqXo3ewb3RNbCr2LLbnG/neul1HMo5BA6cmKoyh2B2FYzKgl9HSG8IeCo80SeYvc/fL/+Otw69BQB4tNujeK7Pc5BwElwovGDUoO566XXcqLoBhUSB7kHd0cG3A4D6vh3hAtHWuy0AfWTHlNjR8lqxguzhuIfrdfR0lbniuT7PAWjY8L0rfRe+T/4eiw8uNmpgd7M4nX8aK06swOrTq5GmTDO7n5bXitUXjUV1AFbK7eviCw2vQWpx/TL//Vn7wYNHJ79OCPUIhZfCC538OgGAUZrUFlKLU/Hw5oeN3o+p91ZaWyrOQg1LzwV8XHz0vp284waRHSo9dybEtbEaSGOdKTgDHjwC3QLhKnXFqfxTZvtCCZMtU5FsjuMwoe0EAM2TyjqccxgPbnoQj2591KzASC1JRWF1IbS8tl6XcEOO5hzFkkNLcCz3mHj9yinPEb2jd7S/Q9yX4zjx/W6/ZjzZ3HZtG4b9NAwv7H3BomjWsqPLkFSYBB8XH6wcsbJFuyWbwmKxc/LkSYwZMwa9evXC/v37MWDAADz77LNYunQpHn74YfA8j9deew3h4eH43//+R6LHgCd7PIku/l1QWluK+bvmiyfq6wNfh4vURRQ7psKThgt/Gpacm6JLQBdxvax3jryDouoiVKmr8PMl9nrT46YD0AsRc8ZpYXxD2wxFhGdEg68Z6RUpVi2sT14vzvqFfieGCL6dr5O+RklNCTr7d8bcHnMR4BaAfiHMZ2HYpl1IYXUP6g4XqQva+zBjdl2xI6SxhKhFXEAcOHDIqcipl9LbnbEbV5VX4Sn3NOtFGtpmKFylrsityEVKsekZlHDBVGlVRsZvaympLsGCPQsw+c/J9araGuKzs5+J/69r+jbkSM4RZJVnwUvuZTYdaQjHcQ2msvZlsjSkkJYEYJdU1un805i5ZSbyK/PRzqcdHu/+OAAgrbS+2BFSWL4uvnCVuZo8nlCCfiz3mDjzJ8+Oc2FJU0HhnBwUPki8/n1w8oN6ndjTlGlILkqGhJOYtQpMiGFix1xVqa1Uqiqx+OBiAMyLZq7i1PD7lXg10WSavVpdjRf/fRG/XPoFj259FA9sfABb0rbgj9Q/xOhuG682Rs8Rrgt7MvegWl0NgN0DXtj7ApQ1Smy9thUzNs9ocPX3v1L/wi+XfgEHDsuGLkO4Z7hVn8HNwGKxM3XqVEyZMgW5ubnYuXMn3n33XTz11FN47LHH8OKLL+Lbb79FWloaEhMTcerUKbz//vvNOW6HQi6R492h78JF6oKLRRfBg8fk9pPFaiVhyYdzN87VSx1cVV5FRlkG5BI5BocPbvS1noh/Ah18O6CougjLjizDP1f+gbJGiQjPCIyMZGJJEDtJhUn1KhQqVBX44/IfAMwbk+sitAH/4eIPUPNqRHhG1PtCAcY3SIVEgXeHvAu5lC0nMa7tOADGRjnDFBYAMbKTWqKPOPA8j/Qy48iOh9xD9AsZRndUGhU+OvkRAGBa52lmqwTcZG5ICE8AoF+ewxCVRiXe9AF2YVDWKE0eqyHOFJzBvYn3imWaH536yKLnJd1IMqoKMdXYUeCPVPa3vK3dbWaFQV3MmZTVWjX2Z7NQt+HfsiGxvuXaFkz4bUKD1X/nCs7h8W2Po7S2FD2CeuDbid+Kn7+pyE5D5mQB4ZwxjuyQ2HEm1Foecr/9eGLPfcguzza5z+n80wBYhP3Rbo/C39Uf10qv4bdLv4n7HMw+iFlbZgFgkzR/V3+Tx+ro2xHeCm/Uamst7tBuCStPrERWeZb4u7nJi6HYqVRXIvFqYr19frv8GwqrC+Gl8IKL1AXnC89j4b8LxWUh7upwV73ndA/sjlCPUFSpq3Ag6wA+Of0Jlh5eCh48JrSdgADXAKQUp2Daxmni52lISlEKlh5eCgCY02NOg9WeLYnFYufSpUuYO3cufH19G9wvISEBP/74IxYuXNjUsTkV7X3bi+kRPxc/vND3BXFblFcUAt0CodKqcK7A2FQrRHUGhA2wqDGTXCrH24PfhpSTYvO1zfjw5IcAmFlMWB4g1CMUbb3bQstr690of7v0G8pUZWjr3dbik1YwKgvhV1NRHYCJEaGq6pnez6CDXwdx26ioUZBwEiQVJiG7PBs8z4tfeuFmKnh2DCuyCqsLUaGqgISTGAksU6ms9cnrcVV5Ff6u/nikm34RTFMIwlD4/A05lnsMZaoy+Lv6o5NfJ1SqK8Xmi5bA8zy+u/AdZm2ehdyKXDF6tjlts0X9jz4/+zkAiILOMFxtSK2mVjSOC74oSzAX2TlTcAZltWXwcfFBfGC8+Liwcnpdk31ZbRneOfwOssqzsGj/IpOtAFRaFd44+AaqNdUYHD4YX477Ej4uPqJwzSrPqtfbpCFzskCfkD7gwCFNmYZKDWtLQJEd50Kj5SH3PofM8usmW3eotCox0twrqBc8FZ54sseTAIBPz3yK0tpSfHL6Ezy5/UkUVRch1i8WbyaYr9jiOA7dg5gf0F4Lgx7LPSY2Mm2sslEQO0Lk6aeUn4y+97WaWnyd9DUA4Nnez2LbPdswt8dc+Ln4AQC85F4YEzWm3nENU1mLDy0WhdGcHnOwfNhy/DDpB8T6xaKougiPbn0U6y+sx88pP2P5seWYu2MuZm+bjRpNDYZEDGl1Ph1DLBY7crm8we0lJSVW7X8r8mDnB7F82HJ8Nf4r+Lr6io9zHIfewSy6U/dEF1NYkQ2nsAzpGtgVs7rOAsD8DR5yj3p+DVO+HZVWhfXJrPfMzK4zRRNbYxgalQ2PXReO4/DBiA/wf8P+Dw/HGZckBroFiqJm+/XtyK7IRm5FLmScDD2CegBg/gxPuSc0vEZMXQlLZIR5hBkZqeualHMrcsUv8YI+C+Ct8G7wPQ1rMwwcOCQXJddrMrgrg6WwRkaOxOxuLG34ffL3FlUzVKoqsWDPAiw/thxqXo1x0ePw6x2/in4AIfJkjpSiFOzJ2AMOHJYPWw6FRIHC6kKTTf1O5J1ApboSAa4B6BrYtdGxCQgLgl4qumQ0kxOqsAaHDzZaVyvUIxQRnhHQ8lrx5gKwdKXQ/yivMg+rTq2q91rfnP8GqSWp8HPxw7Khy8T1cwLdAuEp94SW19YzbZrrsWOIj4uP2BvqajmbQAgLRxKOD8/zTLxybHkXU60iUopSUK2phrfCW5wo3dPpHrT1boui6iJM/nMyPj3zKXjwuKfTPVh/23qTEWlDBJFveJ7bSpW6SiyHn9pxKp6IZ0LheO7xepOX3Ipc5FTkQMpJsWjAIrjJ3JBakmrkk/sz9U/kV+Yj2D0Yd3W4C/6u/pjTcw623bMNy4ctx5fjvzQ7YR4XzSLrJTUl4MBh0YBFmNtzLjiOQ5hnGL6d+C1GRY6CSqvC/479D0sPL8V3F77Dvqx9UNYoEe0djWVDl1l8z2gJbBrZ//73P/z000/i7/fddx8CAgIQERGBM2fMV/jc6nAch4kxE0VDpyFiKsDg5L1RdUOM9JjLI5tjTs854sz/7g53w1PhabTdlG9n27VtyKnIgb+rv5GJzRKmdJwCDhwknMRsZAdg0ZkJMRNMfimE2cW269vEKpy4wDjxC8pxnNhQUchr101hCYiRncIk8DyP94+/jyp1FXoG9bTovQW4BYgRC8PojpbXin6d0VGjMa7tOLTxbIPimuJGG0MCwIoTK7AjfQdkEhleHfAq3hv+HjwVnpjfaz6knBT7svY1WFoqRHUmtJ2AWP9Y9AhmQtBUKsuwRNyai1C0VzSivaNRranG9M3T8dSup3C5+LJ4PMMUloAo1vOZWM+tyBW9TA92fhAAsOHiBqPV5TNKM/DZGeY9WthvYb0JgBAFrOvbsSSNBehTWVfK2DWJIjvOg/in1Imd1JLUeilPIRLSM7ineP7LJXI82/tZAOz66iZzwztD3sGbCW9alOaND2Jix9K19xrio5MfIaMsAyHuIXi+7/PoHtTd7ORFmHR08uuEEI8QsdWIUDWr0qrEwotHuz1qNPFzlbliYsxEcY1Cc++rnU87yCVyvDf8Pdzf2djC4C53x8qRKzGv5zzE+sViaMRQTI+bjtcHvo4149bg1zt+bVVl5qawSex89tlniIxk62ps374d27dvx+bNmzFx4kRKX9mIIHZO55+GWqsGwBoJ8uDRLaAbgt2DrTqei9QFq0avwrye8zC/1/x62/uF9gMHDleVV5FXkQee58XeMw92ftBqJ32kVyQ+GPkBVgxfgSD3IKueKzAmagw4cDhbcBYbr24EoL9hCdT17dQ1Jwt08usEuUQOZY0Sv17+FVuubYGEk+C1ga9ZfOMXBKZhqWnSjSQUVBXAQ+6BAWEDIJPIxJTYuvPrGjQuns4/LZq/V41ahQc6PyCW50d7R4vRsQ9PfmgyLXWl5Ipo4H48nhl4BWO3qRJtoaeRKXHSEFKJFGvGrcHUjlMh4STYk7EHU/+eitSSVEg4icn0Zq8QY5PyqlOrUKOpQe/g3ni5/8uYGDMRWl6LxYcWQ61Vg+d5LDm8BDWaGgwIG4Db291e75ii2KlzE8ur0JWdm6jEMkQwKV9W6sTOTV6/i2g+9JV1+sqlutEd4Vys265jVNQoTO04FQNCB+CHST9YNbETIsbXS6+jpLrE+oEbjO375O8BAIsHLRY9NoKYqhvhF96L4O8U/JTb07fjRtUNJF5JRHZFNgJcAzC141SrxyPhJPj+tu+x7Z5ton/S1D5P9ngSv975Kz4Z8wle7Pci7ou9D/3D+lvsB2xJbBI7ubm5othJTEzEfffdh3HjxuHFF1/EsWPOsTLszaaDbwd4KbxQqa4UeyiIXZMbqcIyR6RXJJ7s8SQ85B71tvm4+CAuIA4A641zJPcIkouS4SZzs3mF2lFRo4w6j1pLkHuQeGE6lHMIQH2xI0R2riqvAtCXnUd7RxvtJ5fKRe+J0Hfo/tj7xccsQUgdHs09Khq5hajO0Iih4uxpcofJCHANQG5FLjalbTJ5LJVWhbcOvSWa0wdH1DebPxn/JFykLjiZf9Jk+/Yvzn4BHjzGRI1BR7+OAAwWvswzDn0L/YxkEhkSwhIsfs8CIR4hWDxoMf6Y/AfGRo8FD3bsnkE9Tc7ghMjOuYJzOH/jvNh/6IW+L4DjOLzY70V4Kbxwsegivk/+HolXE3E45zBcpC54Y+Ab9bqCA/ponZCqFMit1Hl2LIjscOCQU3UdnLSMIjtOhKB1OF1kBzCu5OR53qzY4TgOiwctxlfjv6q39E5j+Lj4iNcac9Gd9NL0Bru312pq8ebBN8VrgeHkQZj01jUpG0apAFZ5Gx8UL3a6F3qCPdLtEZuFh6fC0+QCwc6CTWLHz88PGRmsH8qWLVswZgwzPfE8D41G09BTCTNIJVLxS3ki7wQqVZWin8baFJalCKmswzmHse78OgAs5WWYTrjZGM4qJJyk3oVKKD8XIjtCQ8G6YgfQp7LUWjX8Xf1NRrgaIsYnBm2920KtVeNANqt+EmaPo6JGifu5SF3Estavk7422SfD0JtiaE43JMQjRFw076OTH4nH4XkeZwrOiE0p/xP/H/E58UHxUEgUKKgqMOo2LKSc+oT0qZfCtIZ2Pu2wYsQK/DjpR0yPm45XB7xqdj9fF19Ua6rx/N7nwYPH+LbjRUNnoFsgnu/zPABg9enV+L9j/weAVQ+aa3RoKrLD87wY2WnIoAywG5OQMpa6p1GfHSdC/Fty+r/p+cLzYlVWZnkmblTdgEwiazB9YwuCb8eU2FFpVJi5ZSamJU7DoexDJp+/JmkN0pRpCHANwMJ+xpkQcfJi4NupUFWILTB6Bemvh9NipwEAvjjzBTLKMuDn4od7OzXcF+1WxiaxM2XKFDz44IMYO3YsCgsLMXEia6V96tQpdOjQoZFnE+YwNCkfyjmEGk0NIjwj0NG3Y7O8nmAk3pm+EweyDkDCScSbdkthWC3Q2b9zvRu1MBNLL01HtbrabGQH0IecActMyaYQq7LSWX+ea6XXIJPIMDRiqNF+98feDy+5F64qr2LliZVivwqgYW9KXWZ3mw1PuSdSilPw1qG38NSupzD8p+F4eNPD0PJaDG8zXOxyDMAo9H0sTx9VFf01EdalsMzRNbArXuz3Yr0FYQU4jhNnnVnlWZBJZHim9zNG+9zd8W70Du6NKnUVimuK0cG3g2ikN4WhZ0e48CtrlKjWsM822KPx1K6wTpbU46q4cCTh+IhROo6l/IVoozAZETwucQFxdk+xiBVZJkzKR3KP4EbVDfDg8er+V+v1+UpTpuHLsywK81L/l+pFSeMD4yHjZMirzBPL0c8WnIWW1yLcI9zIlD+u7Tj4uvhCzbPPYEbXGRZV7N6q2CR2Vq5cifnz5yMuLg7bt2+Hpye7IeXk5GDu3Ll2HeCthKFJWUiXmFr40170Cu4FhUQhlgSPjR7baDVCcxPiESJGc4TPw5Bg92B4yb2g4TU4knMEtdpayCQyhHmE1ds3ITwBPi4+GBox1GrDtYCQQtyXtU/sATQgbEA9Eeap8MSj3R8FwLw7d/11F/Zm7AXP81h6eGmD3hRDfF19MbPrTACs0/SejD0orimGQqLAgNABeLn/y/WeI9zQBZNyhapCDINb69dpCoJYB9isM9Ir0mi7sNiqXMIqNd9IeEPss2SKSK9ISDkpKlQVKKhiS1EIKSx/V3+LfGVCuk8iK6E+O06E+LfURXaESZKQyhIKPQwjIfbC0KRc11tn2CfsRtUNLNq/yChCu/TwUqi0KgyOGCxWYBriLncXKyeF73DdFJaAi9RFrLL1VniLkR7CNOaXPzbBG2+8gcmTJ6NPnz544YX6ofjnnnvObgO7Feka0BWuUleU1JSILcmtKTm3FleZK3oF98KRXFaR1dAs+2ayoM8CfHXuK5Mr5goVWacLTosl4JFekSaX0Qh0C8Se+/aIVWK2EB8YD39XfxRVF4mpvtFRpn1Js7vNRrR3NP539H/IKs/C/F3z0TWgK84XnodCosDrA1+3SLjOiJuBqyVXUautRa/gXugR1ANxAXFm1ygTfE0nck+A53kczj4MtVaNKK8otPVpa9P7tgWhSaaX3Esso61LO992WDthLWrUNWbXeRNQSBVo49UG10uv45ryGoLdgy02JwsIwgqchjw7ToSwzpng2RnXdhx+u/wbTuWfwo2qG2JkRzDO25NOfp3gInVBaW0prpdeF79jKq1KvCa9OuBVrDi+AgeyD+Cb89/gkW6P4M/UP3Es9xhcpa5YNGCR2WtB35C+OFNwBifyTuCuDnfpzckGkwmBmV1n4nrpddzW7rYmpatvBay6A2RmZmLixIlo06YN5syZgy1btqC2trbxJxIWIZfKxVmDSquCt8K7Wb6shgyKYEs99A3pK3pcWpqewT2xavQqsy3HhVSWYOA2lcISkElkRj1hrEUqkYrRkQpVBThwZgWo0Jzr77v+ZqupczKcL2Sl1k/0eKLBcRriLnfH8uHL8cHIDzCz60z0DO5pVugAbKYpl8iRX5WP9LJ0/JtlvkS8Oeka0BUrR6zEmvFrGkzV9QjqIQqjxhBMyoJvR2wo2Ig5WUAUwZyGIjtOhD6yw8ROpGckugV0Aw8ef6b+KXr6egb1tPtryyVysbjDsLngsZxjUNYo4e/qj/s63YeX+r8EgPnv9mbsxfsn2KoCc3vObTCCLpqUc49DrVWL6bK6kR2Atcj4cNSHGN92vF3emzNjldj5+uuvkZubix9++AFeXl545plnEBgYiKlTp+Lbb79FUVFRc43zlkEoLQRYfxRxZtpMPNTlIbzQ9wUsG7qsWV/Hnghip6ianW/RXpaJCFsxFDc9gno0WrHgLnfHgr4L8Msdv2BE5AiMjR6LR7o23LG5KbjKXEWRfCTniFF/nZvNmOgxRp6iplK3146QxrI0smModmjVc+dBH6VjKSKZRIYx0SyVJfSbifaORoBbQLO8vuAHNPTtbLvOUlhjosZAKpFiasepGN92PNS8GvN3zYeyRolYv9h6DVXr0iu4FyScBJnlmdiftR+V6kp4yj3FthuEbVgd25dIJBg6dCiWL1+OlJQUHDlyBAMGDMDnn3+O8PBwDBs2DO+99x6ysrIaPxhRD0OfSnOmsARcpC6Y2XVmg91oWxt1y0WjfZpX7AwMGyj6Q8ylsEzRwa8DPh71MVaMWNGgN8UeCL6dDckbxGZpdcv2HZG6FVliJZaFkR1hssBxWorsOBF1IzuGYqdcxZYsaY6ojkBdk7JaqxZ9lmPbsuaoHMfhzYQ3xeVgOHBGnjVzeCo8xRYZX537CgCbZDUlQk3YaFA2pEuXLnjxxRdx4MABZGRkYNasWdi3bx9++OEHe4zvliM+MB6+Lr7wcfGxaOHPW5G6M5zmjuy4y90xI24GOvh2wO3tGzYYtxRit2Al6yydEJbQYOrLURDEjtBrR4zsWCjODSM75NlxHtjfUguOE1Y+lyHaO1o0pAPGUXJ70yOQdS6/XHwZVeoqHM87juKaYvi5+BlNMrwUXnhv+HsI8wjDnB5zRJHUGMIxzhSwhpiN+duIxrHKoCzw6KOP4sMPP4SXl/Gq0e7u7ti/fz/++usvuwzuVsRd7o4fJjGhSIYz0wS5BcFL4SU2+rPUC9MUnu79NJ7u/XSzv46tCL4dlZZ1cL7Zfp3mQvDsZFdko0pdZXGPHQG9cV1DfXacCLYulv7vKfydx0SNweXiywBMe1zsRahHKALdAnGj6gaSC5PFKqxRUaPqFUt0C+yGbfdsM3UYs/QJ6YNvL3wr/k5ip+nYFNn55ptvUFVVf9HDqqoqfPvttyaeQVhDG682LV4C3prhOE5sLugmc7N6KQ1nxE3mZtRXqCX8Os2Bn6sffF18AbDojqXrYgnIOHbj4TgNtDzIt+MkMOGqb2ArCIzxbcdDwkkQ6hGKGO+YZnt9juPE5oKnC06L/X2EBTWbSp+QPuDAqrWknLTVFI84MlaJndLSUiiVSvA8j7KyMpSWloo/xcXF2LRpE4KD6cZDND+CbyfKK6rZ+hA5GoJvp4t/F6cSgEIq61T+KdRoagDA4venT2OxKACtj+UcaLUQ/TqA/u/c3rc9vpnwDb4Y+0WzXxeElNQPF39AUXURfFx80C+sn12O7ePiI6bkOvt3pmaBdsCqNJavry84jgPHcejUqf7K3RzH4a233rLb4AjCHHEBcfjt8m/o4EcVCgL3x96Pi0UXxSUnnIUYnxicyj8lLp8S4BpgsR9JNIbrbowaLQ85+TwdHrVWC84wjcXpb2XNmb4ypEcQ8+0I7RBGRY6ya/XswLCBuFR8Cf1DLWvTQDSMVWJn9+7d4Hkeo0aNwm+//QZ/f39xm0KhQHR0NMLDTfdGIQh7ckf7O6DhNRjRZkRLD6XVEOQehFWjV7X0MOyO4NsROkRbmsICDG6CBmKHcHyYZ0dXicXJWiS62zWgKyScROyQbG61cFuZ02MOIr0ibe7+ThhjldgZPnw4ACAtLQ2RkZGQSJpczEUQNuEmc3O6CAZhGiGNJZQUW9pjBzAsPWc3Rio/dw7UBmKnuds6mMNd7o72vu1xufgyvBReGBA6wK7H91R4YlpnWgLCXtikVqKjo1FaWopt27Zh/fr1+Pbbb41+LOXTTz9FfHw8vL294e3tjYSEBGzevFncXl1djXnz5iEgIACenp6YOnUq8vLyjI6Rnp6OSZMmwd3dHcHBwVi4cCHUarUtb4sgiFaIIHYErIrsGFRjARTZcRa0utJzwDiFdbMRevmMjBzZYqKLsAybzpJ//vkHDz30EMrLy+Ht7W0UQuQ4DjNmzLDoOG3atMGyZcvQsWNH8DyPb775BpMnT8apU6fQtWtXPPfcc9i4cSN++eUX+Pj4YP78+ZgyZQoOHDgAANBoNJg0aRJCQ0Nx8OBB5OTkYMaMGZDL5XjnnXdseWsEQbQyIjwjIJPIoNaySYxNYkfn76Dyc+dAreXFaJ2pdfFuFv+J/w/cZG54pFvzdUgn7INNkZ3nn38ejz76KMrLy1FSUoLi4mLxx5olI+644w7cdttt6NixIzp16oT//ve/8PT0xOHDh6FUKrFmzRqsWLECo0aNQp8+fbB27VocPHgQhw8zo+K2bdtw4cIFrF+/Hj179sTEiROxdOlSrF69mtbsIggnQSaRGTWOtCWNRZ4d58LIs9OCYifUIxQL+y1sdAkZouWxSexkZWXh6aefhru7/crhNBoNfvzxR1RUVCAhIQEnTpyASqXCmDFjxH06d+6MqKgoHDp0CABw6NAhdO/eHSEh+ovf+PHjUVpaivPnz9ttbARBtCyGq7fbEtlhnXa1UGtI7DgDaoOmgi0pdgjHwaazZPz48Th+/DjatWvX5AGcO3cOCQkJqK6uhqenJ/744w/ExcXh9OnTUCgU8PX1Ndo/JCQEubms1C83N9dI6AjbhW3mqKmpQU1Njfh7aWlpk98HQRDNh6FvxzbPDmjJCCeitUR2CMfBprNk0qRJWLhwIS5cuIDu3btDLjc2Zt15550WHys2NhanT5+GUqnEr7/+ipkzZ2Lv3r22DMti3n33XeoHRBAOhCB2OHAIcg+y+HnGYkdLTQWdBE0r8ewQjoNNZ8njjz8OAFiyZEm9bRzHQaPR1HvcHAqFAh06sMZwffr0wbFjx/Dhhx/i/vvvR21tLUpKSoyiO3l5eQgNZTO70NBQHD161Oh4QrWWsI8pXnnlFSxYsED8vbS0FJGRkRaPmSCIm4uw+GuYR5hVjduM96XIjrNguFwEiR3CEmzy7Gi1WrM/1ggdc8euqalBnz59IJfLsXPnTnFbSkoK0tPTkZCQAABISEjAuXPnkJ+fL+6zfft2eHt7Iy4uzuxruLi4iOXuwg9BEK2XLv5d8Er/V7B08FKrnifl9O2SOU5Dnh0nQcsbeHZasPSccBxa9Cx55ZVXMHHiRERFRaGsrAwbNmzAnj17sHXrVvj4+GD27NlYsGAB/P394e3tjaeeegoJCQkYOHAgAGDcuHGIi4vD9OnTsXz5cuTm5mLRokWYN28eXFxcWvKtEQRhRziOw4NdHrTpeWLZOnl2nAa1xqCpoB2XaCCcF5vEjqn0lSFvvPGGRcfJz8/HjBkzkJOTAx8fH8THx2Pr1q0YO3YsAGDlypWQSCSYOnUqampqMH78eHzyySfi86VSKRITEzFnzhwkJCTAw8MDM2fObHR8BEHcOsglclHsUJ8d54A8O4S12HSW/PHHH0a/q1QqpKWlQSaToX379haLnTVr1jS43dXVFatXr8bq1avN7hMdHY1NmzZZ9HoEQdx66NfH0lJkx0mg0nPCWmw6S06dOlXvsdLSUsyaNQt33313kwdFEARhL+RSOaDSeXZI7DgFWp4HGZQJa7DbSp7e3t5466238Prrr9vrkARBEE1Gb2Alz46zYOjZIbFDWIJdly1XKpVQKpX2PCRBEEST0K+PRWLHWWCeHarGIizHprPko48+Mvqd53nk5OTgu+++w8SJE+0yMIIgCHugXzKCPDvOgpo6KBNWYtNZsnLlSqPfJRIJgoKCMHPmTLzyyit2GRhBEIQ9MIzskGfHOWB9dkjsEJZj01mSlpZm73EQBEE0C4Yrn2uo9NwpUGt4cKBqLMJymuzZyczMRGZmpj3GQhAEYXf0kR0tRXacBI1WS00FCauwebmIJUuWwMfHB9HR0YiOjoavry+WLl0KLc2cCIJoReg9O2ry7DgJGkpjEVZi01ny2muvYc2aNVi2bBkGDx4MANi/fz8WL16M6upq/Pe//7XrIAmCIGxFfzPU0tpYToJhU0GK7BCWYJPY+eabb/DVV1/hzjvvFB+Lj49HREQE5s6dS2KHIIhWg5Fnhyex4wxoNLRcBGEdNqWxioqK0Llz53qPd+7cGUVFRU0eFEEQhL2gPjvOB5WeE9Zik9jp0aMHVq1aVe/xVatWoUePHk0eFEEQhL3Qe3ao9NxZoOUiCGux6SxZvnw5Jk2ahB07diAhIQEAcOjQIWRkZNCinARBtCr0aSwtNBoqoHAGjBYCpQ7KhAXYFNkZPnw4UlJScPfdd6OkpAQlJSWYMmUKUlJSMHToUHuPkSAIwmYM18aiyI5zQJ4dwlpsPksiIiLIiEwQRKtHLjVsKkhixxkgzw5hLTZFdtauXYtffvml3uO//PILvvnmmyYPiiAIwl4Yro1FkR3ngC0XQR2UCcuxSey8++67CAwMrPd4cHAw3nnnnSYPiiAIwl6IaSxODS2JHaeAiVbqoExYjk1iJz09HTExMfUej46ORnp6epMHRRAEYS9ouQjnQ6PVkmeHsAqbxE5wcDDOnj1b7/EzZ84gICCgyYMiCIKwF4al5+TZcQ40WkpjEdZhk9h54IEH8PTTT2P37t3QaDTQaDTYtWsXnnnmGUybNs3eYyQIgrAZfZqDqrGcBY2hQZlKzwkLsOksWbp0Ka5du4bRo0dDJmOH0Gq1mDFjBnl2CIJoVRimsTS0ULFTQNVYhLXYdJYoFAr89NNPePvtt3H69Gm4ubmhe/fuiI6Otvf4CIIgmoThchEU2XEONFoeHCiNRVhOk86Sjh07omPHjvYaC0EQhN0hz47zoaHIDmElFnt2li1bhqqqKov2PXLkCDZu3GjzoAiCIOyF0arnJHacAsM0FpWeE5Zgsdi5cOECoqKiMHfuXGzevBkFBQXiNrVajbNnz+KTTz7BoEGDcP/998PLy6tZBkwQBGENxp4dEjvOAFVjEdZi8Vny7bff4syZM1i1ahUefPBBlJaWQiqVwsXFBZWVlQCAXr164bHHHsOsWbPg6urabIMmCIKwFKrGcj40Wh6cjCI7hOVYJYl79OiBL7/8Ep9//jnOnj2L69evo6qqCoGBgejZs6fJrsoEQRAtCXl2nA/y7BDWYtNZIpFI0LNnT/Ts2dPOwyEIgrAvVI3lfBguF0Fih7AEm5oKEgRBOArGBmXqs+MMaLRavWeHmgoSFkBihyAIp8ZobSwNRXacAY2Wp7WxCKsgsUMQhFMjzPw5TgMtT2LHGSDPDmEtJHYIgnBqyLPjfKip9JywkiaJndTUVGzdulVsNsjTrIkgiFaG/mZI1VjOglqrpTQWYRU2iZ3CwkKMGTMGnTp1wm233YacnBwAwOzZs/H888/bdYAEQRBNQW9QJs+Os6DRasT/U58dwhJsEjvPPfccZDIZ0tPT4e7uLj5+//33Y8uWLXYbHEEQRFOhPjvOh4ZXi/+nyA5hCTadJdu2bcPWrVvRpk0bo8c7duyI69ev22VgBEEQ9sCw9FxNpedOgdogskNih7AEmyI7FRUVRhEdgaKiIri4uDR5UARBEPbC0KBMkR3ngDeM7FCfHcICbBI7Q4cOxbfffiv+znEctFotli9fjpEjR9ptcARBEE1FTGNBCw0VUTgFap3Y4cBBKpG28GgIR8AmSbx8+XKMHj0ax48fR21tLV588UWcP38eRUVFOHDggL3HSBAEYTNGaSwyKDsFGl4DGSiFRViOTZGdbt264dKlSxgyZAgmT56MiooKTJkyBadOnUL79u3tPUaCIAiboTSW86HhVQAAKaWwCAux6UxJT09HZGQkXnvtNZPboqKimjwwgiAIe0Bix7ngeR5aUENBwjpsiuzExMSgoKCg3uOFhYWIiYlp8qAIgiDshb70nIfKoIqHcEwM18WiHjuEpdgkdnieB8dx9R4vLy+Hq6trkwdFEARhLwxviBoSOw4PW/KDuicT1mHVmbJgwQIArPrq9ddfNyo/12g0OHLkCHr27GnXARIEQTQFwxuiWuf1IBwXWgSUsAWrzpRTp04BYJGdc+fOQaFQiNsUCgV69OiBF154wb4jJAiCaAKGN0TDzruEY6LheXC6RUDlJHYIC7HqTNm9ezcA4JFHHsGHH34Ib2/vZhkUQRCEvTBsOqfRkthxdDQaw8gOeXYIy7BJFq9du9be4yAIgmgWOI6DlJNCw2sosuMEqA3SWBTZISzF5jPl+PHj+Pnnn5Geno7a2lqjbb///nuTB0YQBGEvZBI5NBoNNDwZlB0dLc8DHJWeE9ZhUzXWjz/+iEGDBiE5ORl//PEHVCoVzp8/j127dsHHx8feYyQIgmgSQiqLIjuOj1rLg6NqLMJKbBI777zzDlauXIl//vkHCoUCH374IS5evIj77ruPGgoSBNHqEG6KJHYcH2PPDokdwjJsEjtXrlzBpEmTALAqrIqKCnAch+eeew5ffPGFXQdIEATRVISbopZXg6fFQB0atVZLaSzCamwSO35+figrKwMAREREICkpCQBQUlKCyspK+42OIAjCDuiXjNCCVoxwbJhnhyI7hHXYdKYMGzYM27dvR/fu3XHvvffimWeewa5du7B9+3aMHj3a3mMkCIJoEkIXZY7TQK3VQiqRtvCICFtRGy4XwVHpOWEZNomdVatWobq6GgDw2muvQS6X4+DBg5g6dSoWLVpk1wESBEE0FVoM1HlQa3iAFgIlrMTqM0WtViMxMRHjx48HAEgkErz88st2HxhBEIS9ENfH4jS6tZUIR4WWiyBswWrPjkwmw5NPPilGdgiCIFo7cgPPjkZDYseR0ZBnh7ABmwzK/fv3x+nTp+08FIIgiOZBLtV5dkCRHUdHo9WvjUVih7AUm86UuXPnYsGCBcjIyECfPn3g4eFhtD0+Pt4ugyMIgrAH5NlxHtTUZ4ewAZvOlGnTpgEAnn76afExjuPA8zw4joNGQy3ZCYJoPRiJHeqz49BoeR6AsDYWVWMRlmFTGistLa3ez9WrV8V/LeXdd99Fv3794OXlheDgYNx1111ISUkx2qe6uhrz5s1DQEAAPD09MXXqVOTl5Rntk56ejkmTJsHd3R3BwcFYuHAh1GrqlEoQBMNI7JBnx6FRk0GZsAGbzpTo6Gi7vPjevXsxb9489OvXD2q1Gq+++irGjRuHCxcuiKmx5557Dhs3bsQvv/wCHx8fzJ8/H1OmTMGBAwcAABqNBpMmTUJoaCgOHjyInJwczJgxA3K5HO+8845dxkkQhGMj9GPhOC3rwEs4LBqtVuyzQ2KHsJQWPVO2bNli9Pu6desQHByMEydOYNiwYVAqlVizZg02bNiAUaNGAQDWrl2LLl264PDhwxg4cCC2bduGCxcuYMeOHQgJCUHPnj2xdOlSvPTSS1i8eDEUCkVLvDWCIFoRgkGZPDuOD/Ps6AzKHIkdwjJsSmM1F0qlEgDg7+8PADhx4gRUKhXGjBkj7tO5c2dERUXh0KFDAIBDhw6he/fuCAkJEfcZP348SktLcf78+Zs4eoIgWiv6myJVYzk6hstFkGeHsJRWI4u1Wi2effZZDB48GN26dQMA5ObmQqFQwNfX12jfkJAQ5ObmivsYCh1hu7DNFDU1NaipqRF/Ly0ttdfbIAiiFULVWM4DeXYIW2g1kZ158+YhKSkJP/74Y7O/1rvvvgsfHx/xJzIystlfkyCIlkO4KTLPDokdR0aj5cHRchGEldgkdtq1a4fCwsJ6j5eUlKBdu3ZWH2/+/PlITEzE7t270aZNG/Hx0NBQ1NbWoqSkxGj/vLw8hIaGivvUrc4Sfhf2qcsrr7wCpVIp/mRkZFg9ZoIgHAeK7DgP1GeHsAWbxM61a9dM9tKpqalBVlaWxcfheR7z58/HH3/8gV27diEmJsZoe58+fSCXy7Fz507xsZSUFKSnpyMhIQEAkJCQgHPnziE/P1/cZ/v27fD29kZcXJzJ13VxcYG3t7fRD0EQzovh2lgkdhwbtlwERXYI67DqTPn777/F/2/duhU+Pj7i7xqNBjt37kTbtm0tPt68efOwYcMG/PXXX/Dy8hI9Nj4+PnBzc4OPjw9mz56NBQsWwN/fH97e3njqqaeQkJCAgQMHAgDGjRuHuLg4TJ8+HcuXL0dubi4WLVqEefPmwcXFxZq3RxCEk6K/KVLpuaNDC4EStmDVmXLXXXcBYN2SZ86cabRNLpejbdu2eP/99y0+3qeffgoAGDFihNHja9euxaxZswAAK1euhEQiwdSpU1FTU4Px48fjk08+EfeVSqVITEzEnDlzkJCQAA8PD8ycORNLliyx5q0RBOHE6D07aorsODhqLa/vs0Ol54SFWHWmaHUzopiYGBw7dgyBgYFNenHegrbtrq6uWL16NVavXm12n+joaGzatKlJYyEIwnmRGax6TgZlx0ar5QEyKBNWYtOZkpaWZu9xEARBNBtGnh1aLsKhMSw9pz47hKXYLIt37tyJnTt3Ij8/X4z4CHz99ddNHhhBEIS90KexqKmgo6PRasmzQ1iNTWfKW2+9hSVLlqBv374ICwsDx3H2HhdBEITd0Ed2tLpVswlHhXl2KI1FWIdNZ8pnn32GdevWYfr06fYeD0EQhN0x7LNDkR3HRkvVWIQN2NRnp7a2FoMGDbL3WAiCIJoFw7WxNFR67tAwsUpih7AOm8TOY489hg0bNth7LARBEM2CsOo5x2lYB17CYWF9dphgJYMyYSk2yeLq6mp88cUX2LFjB+Lj4yGXG59wK1assMvgCIIg7IFh6Tn12XFsjPrsUGSHsBCbzpSzZ8+iZ8+eAICkpCSjbWRWJgiitSGmsTg1eXYcHPLsELZg05mye/due4+DIAii2aDIjvNg2GeHOigTlmKTZ4cgCMKRMOyzQ2LHsdFoeXAgzw5hHTbJ4pEjRzaYrtq1a5fNAyIIgrA3+psiRXYcHVoIlLAFm84Uwa8joFKpcPr0aSQlJdVbIJQgCKKloT47zoOaxA5hAzadKStXrjT5+OLFi1FeXt6kAREEQdgbQ7FDfXYcG7ZcBHVQJqzDrp6dhx9+mNbFIgii1SGksWhtLMdHpdHSchGE1dhV7Bw6dAiurq72PCRBEESTMVr1nMSOQ6PRqsX/k9ghLMWmM2XKlClGv/M8j5ycHBw/fhyvv/66XQZGEARhLwxLzymy49ioDMUOlZ4TFmLTmeLj42P0u0QiQWxsLJYsWYJx48bZZWAEQRD2Qh8B0LCmdITDotaqxP9T6TlhKTaJnbVr19p7HARBEM2GYZ8diuw4NmpKYxE20KQz5cSJE0hOTgYAdO3aFb169bLLoAiCIOyJcTUWiR1HRs2rAQ7gIKHliQiLsUns5OfnY9q0adizZw98fX0BACUlJRg5ciR+/PFHBAUF2XOMBEEQTUJvUNZCTaXnDo1KqwakgJT8OoQV2FSN9dRTT6GsrAznz59HUVERioqKkJSUhNLSUjz99NP2HiNBEEST0KextFBrNC08GqIpCNVYUk7awiMhHAmbpPGWLVuwY8cOdOnSRXwsLi4Oq1evJoMyQRCtDkNvh0qjbmBPorWj5plYlVBkh7ACmyI7Wq0Wcnl9F7xcLoeWQsQEQbQyDEuUDUuXCcdDH9khsUNYjk1iZ9SoUXjmmWeQnZ0tPpaVlYXnnnsOo0ePttvgCIIg7IFhibKaJ7HjyKh5VnpOlViENdgkdlatWoXS0lK0bdsW7du3R/v27RETE4PS0lJ8/PHH9h4jQRBEkzBKY1Fkx6GhyA5hCzadLZGRkTh58iR27NiBixcvAgC6dOmCMWPG2HVwBEEQ9oDjOEgghRYaqMmz49BoQGKHsB6bzxaO4zB27FiMHTvWnuMhCIJoFiScDFpeAxWvanxnotUiNBWkNBZhDValsXbt2oW4uDiUlpbW26ZUKtG1a1fs27fPboMjCIKwF0KpskZLYseR0eqqsSiyQ1iDVWLngw8+wOOPPw5vb+9623x8fPDEE09gxYoVdhscQRCEvRBujuTZcWyENBZFdghrsErsnDlzBhMmTDC7fdy4cThx4kSTB0UQBGFvBLGjoWosh0aI7NCK54Q1WCV28vLyTPbXEZDJZCgoKGjyoAiCIOyNKHa01EHZkdGQZ4ewAavETkREBJKSksxuP3v2LMLCwpo8KIIgCHsj1d0c1WRQdmjEyA6JHcIKrBI7t912G15//XVUV1fX21ZVVYU333wTt99+u90GRxAEYS8ojeUckGeHsAWrzpZFixbh999/R6dOnTB//nzExsYCAC5evIjVq1dDo9Hgtddea5aBEgRBNAXh5khpLMeGp8gOYQNWnS0hISE4ePAg5syZg1deeQU8zwNgPXfGjx+P1atXIyQkpFkGShAE0RQEQystF+HYCJE5OYkdwgqsPluio6OxadMmFBcXIzU1FTzPo2PHjvDz82uO8REEQdgFwbOjJc+OQ6OFBhIAMon5YhmCqIvN0tjPzw/9+vWz51gIgiCaDZno2aE0liPDgyI7hPXYtBAoQRCEoyETIzuUxnJktGBiVS4lsUNYDokdgiBuCeS6tAdFdhwbofRcTmkswgpI7BAEcUsgRnZAkR1HRavlAU4QOxTZISyHxA5BELcEYuk5pbEcFg3Pg9OJHYVU0cKjIRwJEjsEQdwSyMXIDqWxHBWNlgfIs0PYAIkdgiBuCeRS5vHQ0qrnDotaywOcFgCgIM8OYQUkdgiCuCUQDK0U2XFcNIaeHYrsEFZAYocgiFsCwbPDk0HZYdFoDT07FNkhLIfEDkEQtwQKUexQZMdRUWu1+jQWiR3CCkjsEARxSyCTUhrL0TE0KNNCoIQ1kNghCOKWwEWMBGhYvxbC4TD07AjLfxCEJZDYIQjilkDsuMtpWFUP4XAwzw5LY1Fkh7AGEjsEQdwSKITqHU6rS4cQjobaMLJDYoewAhI7BEHcEgh9djhOzYyuhMOhNeizQ2KHsAYSOwRB3BKI1TucFqR1HBM1GZQJGyGxQxDELYG+466GIjsOimGfHRI7hDWQ2CEI4pZAH9nRkGfHQTFcLkJOy0UQVkBihyCIWwIhEsBRNZbDQqXnhK2Q2CEI4pZATHtQNZbDoqFqLMfk/B/AgY8AvuW+d3S2EIQjkHsO2P8B0OUOIG4ywHEtPSKHg/rsOD5qrZY8O47G9YPA7/8BNLWAbxTQ9a4WGQZFdgjCEdi5BEj6FfhlJvDlSODK7pYekcOhj+xooCGDskNCy0U4GAUpwA8PMKHT+XY2WWshSOwQN5ctrwIfdAfyLtyc18u/CFQU3pzXspWqYiDvvPnt6hrg2n72f5krkH0K+O4u4Js72f8JPSUZQOEVk5tEzw600JDWcUg0zmxQLrwClBe09CjsR1kesP4eoLoEaNMPmPIlIJG22HBaVOz8+++/uOOOOxAeHg6O4/Dnn38abed5Hm+88QbCwsLg5uaGMWPG4PLly0b7FBUV4aGHHoK3tzd8fX0xe/ZslJeX38R3QViMugY4/jVQkg78PB2oLm3e18u7AHw6CFg3CdC24sUfN0wDPh0MZBwzvT39MKCqBDxDgGfPAQPmAFIFkLYX+GoskLbv5o63tVJVAnwxHPgkgYncOhhGdqj03DFxWs9OxlFgdX/go17O8X2uKQc23Aso0wH/dsADPwIK9xYdUouKnYqKCvTo0QOrV682uX358uX46KOP8Nlnn+HIkSPw8PDA+PHjUV1dLe7z0EMP4fz589i+fTsSExPx77//4j//+c/Negs3B54HlJmwSye0C38B53613ihWcQNQ1zbttTOPA+oq9v/CVOCvuc1rWLvwF8BrgIJkIGWT5c/TaoGCSzdHIGWdBDIOA+CBcz+b3ufKLvZv+1GAZzAwcRnw1Amgw1hAqwJ+fAjIT27+sbZ2jn4BVBYCmhrgn2fqfV+M01jk2XFE1M7YZ0dVDfw1D9CqgdoyYP0Udu2qS8ElYMP9wNpJ7HrcWtGogV9mATlnAPcA4KFfAY/Alh5Vy4qdiRMn4u2338bdd99dbxvP8/jggw+waNEiTJ48GfHx8fj222+RnZ0tRoCSk5OxZcsWfPXVVxgwYACGDBmCjz/+GD/++COys7Nv8rtpJsryWM5zZVdgxxvm9+N54MQ3TFCYIzkR+HkG8NtsYNNCy2/mSb8D78cCX48HNCrrxm9I2l72b2g8IJEDyf8ABz+2/XiNcWmz/v+WvE7xdWDPMuCjnsDqfsD6qexCZAm1lUyQlmRYJ+BOrNP/PznRtKC9spP9236U/jHfKOD+9UDkQKBGycLFpTmWv6613LgMfD0R2PEWUFvRfK9jKzVlwOFPdL9wTECe/MZoFzIoOz6NLheh1bDvcWs8R82x513gxiXAIxiIncT8LT/PBI6tYdtrK4Adi1mU+tIW4Pp+JnpqK5v2umV51l2rNGogZQtwcBWL3JhCVQ38+giQuh2QuQEP/gwEtG/aOO1Eq/XspKWlITc3F2PGjBEf8/HxwYABA3Do0CEAwKFDh+Dr64u+ffuK+4wZMwYSiQRHjhwxe+yamhqUlpYa/bRKzv8BfDJQf9M++hVQWWR+33+eBtbdbjodUpLBZg8Cx74Efn2UpZYa4tI24PfH2awj+6R50aBRA/88C2x83vwXKO1f9m+/x4AJ77L/71is96NYS22l+dcqzWYzC3As5ZNxBEg3c06k7gC+uQP4MJ5deEqus8ev7mbCUKOu/5wbqUwMregKvB0KvBPGBOkH3YCk3ywbf00Zi7IBACcByrKBrBPG+5Tns0osAGg30nib3BV44AcgoCNQmgl8f2/zpQY3vwSkHwT2rwBW9Qcu/N2iZaT1OPYV8z4FdATGLWWPbX+TXdB1CGKHo9Jzh8VwuQi5RM58LjveYtHN1QOA/4ay7/FnQxzD/5J1Ajj4Efv/7SuB+78D+swCwAMbFwB/P8Xe1/6VLIrbYSzg5gdkHQd+e8y26HNtJfDLI8D7nYCfHmbp34YoSAG2vQ6sjAN+uB/Y9horkqjru6wuBb6/B0j+m11z710LtOlr+pgtQKsVO7m5uQCAkJAQo8dDQkLEbbm5uQgODjbaLpPJ4O/vL+5jinfffRc+Pj7iT2RkpJ1H30Qqi5gQ+WUWUFUEhHZnF3F1FXBqff39eV7/hVFXsRPS0KSpUTPBUl0ChPfWGcXkwIU/2clp7gZ5bT/z1mjVQEg39tje/wFFV+vvu2spcGItu+lknay/vbYCyNSJsJhhTPDE38/STL88Yl1UIv8iu7i9E8bGY4pLW9m/bfqy1wH0n5HRezzARIIgxGKGAXd/wWYkUhfgYqIuxGwQcUn6HfhiBBNJpZn61Bx05eAnv7XsfZz7BVBVAIGdgLi72GPJfxvvI1RdhcYDnkH1j+HuDzz8K5sV5p1jkbuGom/l+UwYW2MQv7afRZckMsAnir3nn6ezz83UuXCzqa1gs00AGPo88zSF9WARr62viLsZeXY0JHYcEU3dyM7mF5kAv5gIFFxkURGAnZcb7rM9wlOazcTy9/cB1w/ZafR1UNcAf84DeC3QbSrQ5XZm4L39A2D4S2yfk98Cygz2vZv2A/uuT/uBXZtSNrJJiDWTjtJsYN1twPnf2e8XE4HPh9UvdNCo2UTsy9HMS3TwI6A8D3APZN7BG5eAL0cBpzew/cvy2HGv7QMUXix1FTuxyR+RPWm1Yqc5eeWVV6BUKsWfjIyMlh4Sozwf2PU28HFvFh3gpMCwhcBju4BBT7F9jn1ZX81fP8BOVpkrE0aVhSzqIOR19/4PSD/ETsJ7vgbi72NfGoUnu8l/o4sGGXpyMk+wUKm6Gug0EXh8NxAznP1eN3qTnAgc+ED/u/BFMiT9EBNNvlGAfwzrE3P7B0BwV6Ain4U+TUVQDCnJYBeHTxPYlxQwH1K9tIX922mC/rO7uJFFZAQqi9jsiNeykshnzwEz/wF63A90Gg/c9w37G5z9Edi8kF2cNi1kY60tA6KHALN3AM+cAV7OAJ7Wibxr+y3LqQsprD6zWO8cgIkdw89W8Ot0GG3+OH5tgQd/AuTuLBr1fixLfe7/gF2oc5OAfe8DX40B3uvERPRng9mMuLHIHs+z/QCg90xg3hF2TkoVLFT9ySDT4vZmcmIdUHmDfQ7d7wWkMuCOj1i0LOk34PIOAOTZcQZqNWpwHPvbyTgZoMxiGwbOBR7+DXjmLDDvGODmzyLRvz7a+HXFkLwLwB9zgA/i2TXt8lZg7QTgz7n2jxT9+x7zE7oHAhP/T/84xwEjXwUmrWDXy6HPs+9d59vY9ugEYMrnADh2PzA1iTNF1kkmULJPsc/njg8B32gWyV4zDjj6pS4d/CkzSv82m0WQJDKWXrv/e2BBMjDnIEupq6uAP+ewHjpfj2MRaI8gYFYi0G64fT8rO9BqxU5oaCgAIC8vz+jxvLw8cVtoaCjy8/ONtqvVahQVFYn7mMLFxQXe3t5GPy1KwSXg76eBld2Af/+PheMDY4HZ24FRiwCZgl3EXX1ZJZMQtRA4oDvZez4IPPw7O4GL09jM5tJWdkwAuOMDJjQAoN0IdlK6B7J0z5oxwLJI5svY+hrw/VSgtpxFOu5dx8Zw+0o2o7iyS59+KbzCTniARY0A4Pyf9b0nV3V+nZhh+scU7ixsq/BiYshclIbn9SLw9HomTjrfzm5utWX1jb21lcDVPez/nSYAQbFMsIEHDq3SH/OveSx1FNABuOszdmExJHYiMOULsIvKV+wCePQLtm3o88CMv4DIfmwcrt6s6iCsB4tWCWLMHFkn2ecuVQA9HgA6jGFitfgakJekH6OhObkhInoD930LuHgzsZuyCdjxJrtQfzaY9enJPMY+A/927DPcv4Jd/IQ0mSkubQEyj7L8+/AX2d9s1CJgziEgKkF3wZvbuGgC2D6pO9lsNHEBkLK56b4DVRVw4EP2/6HPM6EDAOE9WYQHADY+B9RWGCwvoIGmNaXgCItRGwgXmUTGotUAi952GAP4RQNBnVj1j8yVnb+bFxpPILRaNpk7tZ6dO9vfBP6azywAnyYAZzawlFHUIPbdBIDT3wOr+rDrgD0KF3LOsu8fAEx6D/AIqL9Pv9lsAjb6jfqVTF3vBsb/l/1/+xvsO3VxY31BplGxCd7xtcDa24CyHCCoM/D4LjbJeuJfdi3V1AKbXgD+rwOw5WVWReUeCIx4lQmcBzawyJNMwczGD/0KjHyNTSjO/sSuW35tgUe3su9eK6TV2tljYmIQGhqKnTt3omfPngCA0tJSHDlyBHPmsItYQkICSkpKcOLECfTp0wcAsGvXLmi1WgwYMKClhm4dhz8Dtryk/z2iLzD4aXYCGvYkULgDvWcwFX/0C73Kz7/IZh/ggIT5rFrn4d+ANWNZPnjDfWy/XtOB7vcYv3Z4L2D2NnZTvHaApczSD7IfYSzTfmDeEIAZzYYvZMJjy8tA9CCW860pZTe+h34B3u/C0hyZx4Aog7+BmCaqo/gD2jMR9ttsJspihhoLIp4Htr6qN5+2HQqMWczSU4dWs23H1gB9HtF3FU7byyJQPpFASFf22KCnmPfpzA/sS3r+dyYIpArgnrWAi6fpv0/3e9hsJ/FZoDyX5cunfAl0HGt6/7jJTMRc+EuXezeDENWJm8xSUQDQfjQLTSf/wyJ0eUks6iV3ByItOJ87jgUWprILacYR/U9VCZtpxU5k4s87nI0v8Tn2Gl+MBEa8BAx+Ti8WAHZT2Knzvwx4AvAymEAEdmAzvU8GsNnpv//HRFBd1DUsbZayiQmdWoMo3PE1TES1G8HGFn8fIHcz/d40ahadkimA7vcBPhHs8ZPfsfC6TyQQP834OSNfZZGyknTgwIeQ93kIAFsbi5oKOiY1BilamUSm95u4+RrvGDUAmPoV8NN01u7CK4yl4lM2sQlghfEkWYSTsCjvoKf1fpO+jzL/TO45FtU+8xOLkPvaaH9QVQF/PMEi3V3uZMLFFhLmsWj3kU+BI5+xH4BNZvzbMwFSnMZeR6DDWDZ2V90E382XFToc/oSJJnU1m/wlzAd6TDP/fZRI2eQncgCb7HiHseuBV4jp/VsBLSp2ysvLkZqqTyukpaXh9OnT8Pf3R1RUFJ599lm8/fbb6NixI2JiYvD6668jPDwcd911FwCgS5cumDBhAh5//HF89tlnUKlUmD9/PqZNm4bw8PAWeldWcCOVnWAAuwkNfhaIGmh+KYB+j7HIxNXdLBoU1Ak4pDMMd56kd70HdmQzm2/uZGW4gbHARDNRk4D27GTneRalyTjCKlk0KmD8O/VFwKBngHO/sRvcp4PYzMojmEV/XLyYCDv7E7vBCWKnqlhnFgYTK3Xpfg+LxJz6joVEnzygn+nsWaYXOnd+zESb8Pn0fJDdjPOSWJ8K4fVSNus/U2Hf6EFARB8mADe9oC9FH/dfICze9Gcj0PcRQObC0lMjXmn4Ihd3F4uiXN3L0mSCkDGkulQfGevziMFz72Ri58Lf7EYtRHXaDmWvbwkyFxZtiuwHYL5+Rlv3nIqbzATqP8+y19z1Nqu0uPszdv4ArGNz/nnAxQcY8mz91/IIAG57j3V13reC3STCeui3VxYBP0xj55SAZwj7u0gVbNatzGAi9NJmlnKa/icgMRFwPvghsOcd9v+dS5hZu+eD+vTpkGeZEDLExRMYuoCJuoyjkPWbqfsstOTZcVBUWgOxo9Xo/XKuvvV37nIHMGEZm0zu/q/xNhdv1ujOI4h9R939WZl0uxFMLBgS2R94fA8T57veZpHOz4cyb1+ncda/ia2vAvkX2HVz0vvWP9+Q8e+w8V3dw66BBcnMr2TopZO7s+t85zuAYS/Ub+zHcUw4tRvBJg4xI0x/B03RbjiLPnFcq1/CpkXFzvHjxzFypL7CZMGCBQCAmTNnYt26dXjxxRdRUVGB//znPygpKcGQIUOwZcsWuLq6is/5/vvvMX/+fIwePRoSiQRTp07FRx9ZmMNsSXieRQs0NWxG/8CPjZ8sftEsHZOykUV3hr0AnNWlcAY/Y7xv1EAWejy1nkUyFB4NH5vj2Gw9sAPQ6yHz+8kULBLz9XgmdDgpEzrCrL/r3UzsXPiTfRElEl21Fc9El3eY6eNO/B+7Kd64xNJiD/7EIjd7l+m2/x+LbBni5seMfafXs/By1AD2uQppvtgJxu9v0FPMr3LhT932SUD/xxv+XAR6Psh+GiOgPRDSnZmFL24Eek+vv0/Sr3pjcvQg/eOdxrP8eEEyK/VONVFybi0NnVOewcC074EzP7IweNZxVsUyZjGbzQo3iCHPsM/aFF3vApLuZBGUP+ex8LhMwbwU66cw06irD9D/Pyx6E9ZLfyG97f9Y5+iUzSykn7aXzVIT5hm/Rv5FJnoB5vHKP88M00JJvlcY0PNh0+Pz1kWAqor1aSzy7DgstQaRHWl1me5/HBMvphj4JEvdHPiApaljb2NiO3pwfXHcEFIZi252Gs+uIdmnWNO8Ic8BIxcZR0Qb4vyfLNIEMN+NZ3CDuzeKRAJ0m8J+ADaxzDzBJhF+bdnExSvcMvES0lUfCbd2DA5Ai4qdESNGgG8gd85xHJYsWYIlS5aY3cff3x8bNmxojuE1L6c3MOe6zA24fYXlqnjAf5jYOfMD+11Ty0KJkf3r79thDPuxN1EDgQFPsrDpuLeBtoP129qPYpGAshwWIYoeZFzpZA6FB0snfTmKpeV+mKY3GY96nb1vU/SbzcTOhT9ZObsyg6WbFJ71o0hd7mQXgOJr7CY4eVXzzEa6TmZi58Kf9cUOz7P8OcDSXIav7+bH0nxXdjLBmK6rAmnInNxUOA7o+QBLH/41n0UNt7zMomkl6Wz2OeDJho8x6X0maPPOsRLZrncB301h6UyvcGD670BwF9OvHdqN/XgEsAjMjrdY1CYkju2j1TBvlaYW6DiOVckVXWXfnzM/AKVZLAomd61/fEAv0qqK9ctFcFqoaL0Ih0SlS8lwkIKrUbIHXX0avuGOfYtNBt38mv59F3wp2xaxCef+lazDeZc79dEhN3/mjawb1S2+zryZABNJTZnEmMPND+jYDNd8J8AxJJmjYq4KoOIG61UAACNeZl8gS4kZzgxmteXMiQ+w/PLNZsIyYMFFIGGu8eMyF2ZkA1iJNmCZ2AHYTW+CLlUhCJ3BzzDjqTkiejPvkaaWpcFSdM9rP7J+6kciBSb8jxmp7/vWdIrJHghl5Ff3sJmWIZnHgNyzemNyXYSF8g6uYu/JJ5Ll0JsbnzbA9D+YcJG7M6ED6EzJjUQFPYOBicvZ///9Pxb1K81k7RJmbzMtdOrS5xEmZjQ1LJUpVAYeWs0iTi7erHqP41j0bPTrLHz+fEr9iJ8hotgpgVyqX0uptinNMYkWQzAoc5Ca9+uYwt3ffhMbmQuLSt6zVl9csfUV5sP5/h7gq1Gs6vGPOfo17zQq5kusUbL02cjX7DMWwmJI7DQnm55nvpnLO4yrAba+xm6CId3qh+wbg+OMUy/+7VumnwHHmU9JCYa7C3+xvg4FFwFwQNshjR+372z98/s+Cox5q/GLVL/H2L/Hv2ZRL4CFqk0ROwH4z+7mbXYV2JGlW7Rq4KLBMhUVN9gFD2DpN1Niq/PtADi9F6H9qJuXC+c49lk+uZ99frGTWLm5JXS/h6UItCp2bkf0ZTNgS02cHAfcuYrNivPO6brKpupTaeP/qzclC0ikxqZpUwhip0YJGa+/3NVqzUxEiFZNrc6zI4FBJZYpv87NoNsU4Im9wMB57PvcbiTrh+Udwb4HZzYwX+N3U1hzwMxjLOo99SvAQHgTN4dWW43l8NSUMxNqbTnzIgTHMc+IRxDr2wKO9QKx5aSPnwbsWMJmCYPmt+hKsiZpN4LdZCry2U0LYCZgSyIpHAdM/ZrNfAI6WHaj7zqFCUghGgEO6Dje1tHbh7jJzFty4S/mgVLXsoZ/JemAXwzzM5nCM4il/q4fYL83R6i7MQLaM8+UNXAca01QlssMnnd+1HhEqC5eIaz3x8/TmcfiYiKrDmk/ihnTbcHgRihX6RvMqSmy45CoBbHDGUR2XH1abkAB7fXRaEMyj7PK2eR/9N4yALjzQ+si+YTdoMhOc+HiCcw9xFS/wpO57/+cw8KcAIvOtOlj+7GnfgUMe9G8MbMlkcr16ZiT37F/G0thGSKRsOiIpRENhTvQ08BU3aav6W7DNxOhSeCVXUC1kvX6uH6Ahb0f+LFh4Sd8dpykVTbnMotXKIua3bPGeqEjEHcn+1vyWmZWV3jqGgTaGN2SykTzqkw0tBqXMBOOg6FnR4zsWJLGutm06ctS5U+dYOZ8V19mN7C1zJxoMiR2mhPfKKb6n0sCRr8JeOpC7l7hzHTbFDqNA0a9Zl1Fwc1E/FLr0nd1++vYm76P6v9vLoV1MwnuzLxVWhVb1O/EOgAc63ER3Lnh53abytrDd7/PfBWUMzNhmb7B49gltvczEdDdDLnqEgiXPBWJHYdEpfPsSDmDHjstlcayBP92zN/z8nX9mm1Ei0BprJuBmx/r95Ewj830g+P0TZ2clbbDWAfOyhusnDoqoXlfL7ADEweXttRvnthSxE0G9l5kFU4AqwqxpC+HZzDwXANdjZ0dV2/gkc1AfrJ9qgnd/Fj6sKoYEsigRa3o/SAcCyGyY+TZMRPZ0Wg0UKno7+zoyOVySKVNt2qQ2LmZyFxa3eJozYZUxlISx79mZlVzHYrtyd2fA+Bbj4cp7i79Ehg9HmiZqjlHxacN+7EHwsy/qpilP6C/aRKOhYbXiR0jz46v0T48zyM3NxclJSU3dWxE8+Hr64vQ0FBwTSjWILFDNB+DnmZdmQffpJt8a2tuFdyF5eurlfqyaeLmY9BrRwIpNNCnQwjHQi1EdjjzkR1B6AQHB8Pd3b1JN0iiZeF5HpWVleIamGFhZiqALYDEDtF8+McAM/9u6VG0HBzH8vVEy2IodjgZwBsvO0A4DkI1ltRMZEej0YhCJyDAxOKahMPh5sbW58rPz0dwcLDNKa1WNhUmCIKwMwZiR0hjqUnsOCRqvmHPjuDRcXevs0o44dAIf8+meLBI7BAE4dwYiB2pbn0sairomAhprMaqsSh15VzY4+9JYocgCOfGYMkICUV2HBpR7EhaeZ+dZmDPnj3gOE40Xq9btw6+vr4tOiZHgsQOQRDOTV3PDvQ3TcKx0PAaAIAUEkBVyR5szX12rOTQoUOQSqWYNGlSSw/FbixevBg9e/Zs6WGQ2CEIwskxEju60nOqxnJIhNJzo0V2WnK5CDuzZs0aPPXUU/j333+RnZ3d0sNxKkjsEATh3Jjw7FBkxzER/m4KYWFlF5/W01eriZSXl+Onn37CnDlzMGnSJKxbt67Jx8zMzMQDDzwAf39/eHh4oG/fvjhy5Ii4/dNPP0X79u2hUCgQGxuL7777Ttx27do1cByH06dPi4+VlJSA4zjs2bMHgD61tnPnTvTt2xfu7u4YNGgQUlJSALBU21tvvYUzZ86A4zhwHId169aB53ksXrwYUVFRcHFxQXh4OJ5+unlblFDpOUEQzo3g6TAUOzyJHUdEH9nRiR23hqM6PM+jSqVp7mGZxE0utcpY+/PPP6Nz586IjY3Fww8/jGeffRavvPKKzebc8vJyDB8+HBEREfj7778RGhqKkydPQqvVAgD++OMPPPPMM/jggw8wZswYJCYm4pFHHkGbNm0wcuRIq17rtddew/vvv4+goCA8+eSTePTRR3HgwAHcf//9SEpKwpYtW7Bjxw4AgI+PD3777TesXLkSP/74I7p27Yrc3FycOXPGpvdpKSR2CIJwboTIDq+BDOzGoSGx45AInh2F7t/G/DpVKg3i3tjazKMyzYUl4+GusPwWu2bNGjz8MFvYecKECVAqldi7dy9GjBhh0+tv2LABBQUFOHbsGPz92cLDHTp0ELe/9957mDVrFubOnQsAWLBgAQ4fPoz33nvParHz3//+F8OHs/UPX375ZUyaNAnV1dVwc3ODp6cnZDIZQkNDxf3T09MRGhqKMWPGQC6XIyoqCv3797fpfVoKpbEIgnBu5G6AzJX9VxcRUNNCoA6JlhfSWCw64SyVWCkpKTh69CgeeOABAIBMJsP999+PNWvW2HzM06dPo1evXqLQqUtycjIGDx5s9NjgwYORnJxs9WvFx8eL/xe6HAtdj01x7733oqqqCu3atcPjjz+OP/74A2p1805AKLJDEITz4+YHlOWIXg9KYzkmGl4NcIBC+Ps1Etlxk0txYcn45h+Ymde2lDVr1kCtViM8PFx8jOd5uLi4YNWqVfDxsd6ELXQethWJbvkdXvBHwXxTP7lcbxkX0m5CuswUkZGRSPn/9u49Lqo6/x/4axhuAwgjch1EvISohYoYyPZNWWW9rJlubhmZtxSywEtq67p5WdsWfeTiNVd79BOyh4qaq/ZQNs0LWglKgiRokRrKrgqULqByn/n8/gBOTCAMyjDMmdfz8eDxgM85Z+b98QNz3r4/n3NObi6OHz+OY8eO4c0338SaNWtw+vRpvddqS6zsEJH81U1l2dZVdjiNZZ7qx822foF5C5UdhUIBB1trk3wZutampqYGn3zyCeLj45GVlSV9ffvtt9BoNEhKSnqkf6v+/fsjKysLd+/ebXJ73759cebMGb22M2fOoF+/fgAAd3d3AMDt27el7Q0XKxvK1tYWWm3jdVMqlQrjxo3Dxo0bcerUKaSlpSE7O7vVr28oVnaISP7qk5266Q9ejWWefpnGMqyyYw4OHz6M//3vf5g5c2ajCs7EiROxbds2zJ49u9WvGxkZibi4OEyYMAGrVq2Ct7c3Lly4AI1Gg7CwMLz99tt46aWXEBQUhIiICBw6dAj79++XFhKrVCoMGTIEq1evRo8ePVBUVISlS5e2Oo7u3bsjLy8PWVlZ6Nq1Kzp16oSkpCRotVqEhobCwcEBO3bsgEqlgp+fX6tf31Cs7BCR/EnJTu3/MFnZMU+6uvGzq78DtgzW7Gzbtg0RERFNTlVNnDgR58+fx8WLF1v9ura2tvjiiy/g4eGB3//+9wgMDMTq1aulB2lOmDABGzZswD/+8Q88+eST+PDDD5GYmKi3IDohIQE1NTUIDg7G/Pnz8d5777U6jokTJ2L06NH47W9/C3d3dyQlJUGtVuOjjz7CM888g/79++P48eM4dOiQUR/eqhANJ+QsVGlpKVxcXFBSUgJnZ2dTh0NEbe2zGODCDsT0/D98KfLhi4n497S/mjoqaqX/+39zUWKTgsgaV/zlP1nA2LXA0zOl7RUVFcjLy0OPHj1gb29vukCpTTU3roaev1nZISL5q5vuqJ/+qL+EmczLL5WdytoGGVR2qH0w2SEi+aubxrKrW6uj4zSWWdKhLtnRVtU2yGDNDrUPJjtEJH9SslN7kuSaHfMk6hco11TUNrCyQwZiskNE8lef7NRVBHScxjJL9ZUdKdlhZYcMxGSHiOSv/mosbe1aDy14B2VzVJ+k1o+j9CgQohYw2SEi+ZMqO7UnSVZ2zJOoq+xIN4izb/2dhckyMdkhIvmTkp3a6Q8mO+ZJoO6p50IAds6AleGPZCDLxmSHiORPmsaqnb7SgQuUzZFOquwIrtehVmGyQ0TyZ9cJUCil6Y/66RAyL9I0lgCg4hQWGY7JDhHJn0IBqNSwrrthPO+zY55+SXZY2XkU4eHhmD9/vqnDMAkmO0RkGVSda9d6gJUd89VggTLvsdMu5s6di+DgYNjZ2WHgwIGmDueRMdkhIsug6ixNY+mY7JgloWBlxxRee+01TJo0ydRhPBYmO0RkGVSdpWksVnbMk940lowqO+Hh4YiNjUVsbCxcXFzg5uaGZcuWoeFzuisrK7Fo0SL4+PjA0dERoaGhOHXqlLT9zp07iIyMhI+PDxwcHBAYGIikpKRm3zc5ORkuLi7YuXPnQ/fZuHEjYmJi0LNnz8fupylZt7wLEZEM6CU7XLNjnhpMYxlS2RECqC4zZkAPZ+NQu1bMQNu3b8fMmTORnp6O8+fPIzo6Gt26dUNUVBQAIDY2FpcvX8bu3buh0Whw4MABjB49GtnZ2fD390dFRQWCg4OxePFiODs7Izk5GVOmTEGvXr0QEhLS6P127dqF2bNnY9euXXjuuefarNsdFZMdIrIMDaaxWNkxP0IIoG4ay8bQyk51GRCnMW5gD/OXW4Cto8G7+/r6Yt26dVAoFAgICEB2djbWrVuHqKgo5OfnIzExEfn5+dBoavuzaNEiHDlyBImJiYiLi4OPjw8WLVokvd6cOXNw9OhR7N27t1Gys3nzZrzzzjs4dOgQhg0b1jb97eCY7BCRZeACZbOmE5CSHWsB2a3ZGTJkCBQNKkFhYWGIj4+HVqtFdnY2tFotevfurXdMZWUlunTpAgDQarWIi4vD3r17cfPmTVRVVaGyshIODg56x+zbtw9FRUU4c+YMnn76aeN3rINgskNElsFeXXuSBAAmO2ZHqxOAQgeg7qaChlR2bBxqKyymYOPQ8j4Gun//PpRKJTIyMqBU6t812snJCQCwZs0abNiwAevXr0dgYCAcHR0xf/58VFVV6e0fFBSEzMxMJCQkYPDgwXoJlpwx2SEiy6DqDBvUVXYUTHbMjVYnUJ+k2ggB2BvwEFCFolVTSaZ07tw5vZ/Pnj0Lf39/KJVKBAUFQavVoqioCM8++2yTx585cwbjx4/Hq6++CgDQ6XT44Ycf0K9fP739evXqhfj4eISHh0OpVOKDDz4wToc6GF6NRUSWocECZVZ2zE+NTgdFw2ksGV2NBQD5+flYsGABcnNzkZSUhE2bNmHevHkAgN69e2Py5MmYOnUq9u/fj7y8PKSnp2PVqlVITk4GAPj7++PYsWNITU3Fd999h9dffx2FhYVNvlfv3r2RkpKCf/3rXy3eZPDq1avIyspCQUEBysvLkZWVhaysrEYVo46OlR0isgyqzrCpy3VY2TE/Oh30p7FktmZn6tSpKC8vR0hICJRKJebNm4fo6Ghpe2JiIt577z0sXLgQN2/ehJubG4YMGSJdSbV06VL8+OOPGDVqFBwcHBAdHY0JEyagpKSkyfcLCAjAyZMnpQpPfHx8k/vNmjULp0+fln4OCgoCAOTl5aF79+5t1HvjY7JDRJZB1bn2JAkA0EKnE7Cysoz1CnLQqLJjL69nY9nY2GD9+vXYsmXLQ7evXLkSK1eubHK7q6srDh482Ox7NLwvDwD07dv3odWfhx1jrjiNRUSWQdVZWqCsUGihbXDDNur4arS6Xyo7No6Akv9XJ8Mx2SEiy2DvIq3ZUSi0dQteyVxUan+5EaS1zKo6ZHxMjYnIMiitobStvUwXTHbMTmVNtfS9jZ28kh25TBV1ZKzsEJHFsLF1rv1GoUUNkx2zUqX9JdmxltmVWGR8THaIyGJYq+oqAgpRuwaEzEZlw2RHZpUdMj4mO0RkMWwaXK7c8ORJHV99ZcdKCFg5GHBDQaIGmOwQkcWwUf1ykqyoMa+bolm6qro1O3J8LhYZH5MdIrIY1ipX6fuGC16p46uquxrL4OdiETXAZIeILEbDZKea01hmpUJbCQC1tw9gZYdaickOEVkMKwdXKOvutcM1O+alRldf2QGg4pqdRxEeHt7is7DkiskOEVmOBg8DZbJjXiqlNTus7LSXb7/9FpGRkfD19YVKpULfvn2xYcMGU4f1SHhTQSKyHPZqWAOoBFBZzWTHnFTXVXZsZPjE844qIyMDHh4e2LFjB3x9fZGamoro6GgolUrExsaaOrxWYWWHiCyHqjNs6io7VazsmJX68ZLjE8/Dw8MRGxuL2NhYuLi4wM3NDcuWLYNo8Py2yspKLFq0CD4+PnB0dERoaKjenZfv3LmDyMhI+Pj4wMHBAYGBgUhKSmr2fZOTk+Hi4oKdO3c2uf21117Dhg0bMGzYMPTs2ROvvvoqZsyYgf3797dJv9sTKztEZDkaPAy0Ssdkx5xU1y1QthGGX40lhEB5TbkRo3o4lbUKCoXC4P23b9+OmTNnIj09HefPn0d0dDS6deuGqKgoAEBsbCwuX76M3bt3Q6PR4MCBAxg9ejSys7Ph7++PiooKBAcHY/HixXB2dkZycjKmTJmCXr16ISQkpNH77dq1C7Nnz8auXbvw3HPPGRxnSUkJXF1dW96xg2GyQ0SWQ9W5tjIAoIr32TEruspSAPX32THsDsrlNeUI3RVqxKge7twr5+Bg42Dw/r6+vli3bh0UCgUCAgKQnZ2NdevWISoqCvn5+UhMTER+fj40Gg0AYNGiRThy5AgSExMRFxcHHx8fLFq0SHq9OXPm4OjRo9i7d2+jZGfz5s145513cOjQIQwbNszgGFNTU7Fnzx4kJycbfExHIZtprM2bN6N79+6wt7dHaGgo0tPTTR0SEXU0KrW0QFlbd/Ik86Crqh0vKygApY2Jo2l7Q4YM0asEhYWF4cqVK9BqtcjOzoZWq0Xv3r3h5OQkfZ0+fRrXrl0DAGi1Wvztb39DYGAgXF1d4eTkhKNHjyI/P1/vffbt24e33noLx44da1Wik5OTg/Hjx2PFihUYOXJk23S6HcmisrNnzx4sWLAAW7duRWhoKNavX49Ro0YhNzcXHh4epg6PiDoKGxWUovaEoq0sNm0s1Cqi6h4AQCEM/z+6ylqFc6+cM1ZILb53W7l//z6USiUyMjKgVCr1tjk5OQEA1qxZgw0bNmD9+vUIDAyEo6Mj5s+fj6oq/QpmUFAQMjMzkZCQgMGDBxs01Xb58mWMGDEC0dHRWLp0aZv1qz3JItlZu3YtoqKiMGPGDADA1q1bkZycjISEBPz5z382cXRE1JEo6wraRT9lIyvX18TRkKGKS68CABRQtrDnLxQKRaumkkzp3Dn9pOzs2bPw9/eHUqlEUFAQtFotioqK8OyzzzZ5/JkzZzB+/Hi8+uqrAACdTocffvgB/fr109uvV69eiI+PR3h4OJRKJT744INm47p06RKGDx+OadOm4e9///tj9NC0zD7ZqaqqQkZGBpYsWSK1WVlZISIiAmlpaU0eU1lZicrKSunn0lKWs4kshaIu2UkoPYSEs4dMHA21lsL8T1tNys/Px4IFC/D6668jMzMTmzZtQnx8PACgd+/emDx5MqZOnYr4+HgEBQXhp59+wokTJ9C/f3+MHTsW/v7+2LdvH1JTU9G5c2esXbsWhYWFjZKd+tdLSUlBeHg4rK2tsX79+iZjysnJwfDhwzFq1CgsWLAABQUFAAClUgl3d3ej/VsYg9n/1vz888/QarXw9PTUa/f09MT333/f5DGrVq3CypUr2yM8IupgBtr2wy1dDnQw/EoZ6hiUEBjYKdjUYRjF1KlTUV5ejpCQECiVSsybNw/R0dHS9sTERLz33ntYuHAhbt68CTc3NwwZMkS6kmrp0qX48ccfMWrUKDg4OCA6OhoTJkxASUlJk+8XEBCAkydPShWe+sSqoX379uGnn37Cjh07sGPHDqndz88P169fb9t/ACNTiIYX8puhW7duwcfHB6mpqQgLC5Pa//SnP+H06dONSoNA05UdX19flJSUwNnZuV3iJiKitlVRUYG8vDz06NED9vb2pg7HYOHh4Rg4cOBDKyyWrrlxLS0thYuLS4vnb7Ov7Li5uUGpVKKwsFCvvbCwEF5eXk0eY2dnBzs7u/YIj4iIiEzM7C89t7W1RXBwME6cOCG16XQ6nDhxQq/SQ0RERJbJ7Cs7ALBgwQJMmzYNgwcPRkhICNavX48HDx5IV2cRERF1VA0f+0DGIYtkZ9KkSfjpp5+wfPlyFBQUYODAgThy5EijRctERERkeWSR7ACQHqJGRERE1JDZr9khIiJqyMwvMqZfaYvxZLJDRESyYGNT+8yssrIyE0dCbal+POvH91HIZhqLiIgsm1KphFqtRlFREQDAwcHBoGc/UcckhEBZWRmKioqgVqsbPResNZjsEBGRbNTfX60+4SHzp1arH3rfPEMx2SEiItlQKBTw9vaGh4cHqqurTR0OPSYbG5vHqujUY7JDRESyo1Qq2+QkSfLABcpEREQka0x2iIiISNaY7BAREZGscc0OfrlhUWlpqYkjISIiIkPVn7dbuvEgkx0A9+7dAwD4+vqaOBIiIiJqrXv37sHFxeWh2xWC99WGTqfDrVu30KlTpza9AVVpaSl8fX3xn//8B87Ozm32uh2NJfTTEvoIWEY/LaGPgGX00xL6CFhGPx+1j0II3Lt3DxqNBlZWD1+Zw8oOACsrK3Tt2tVor+/s7CzbX9CGLKGfltBHwDL6aQl9BCyjn5bQR8Ay+vkofWyuolOPC5SJiIhI1pjsEBERkawx2TEiOzs7rFixAnZ2dqYOxagsoZ+W0EfAMvppCX0ELKOfltBHwDL6aew+coEyERERyRorO0RERCRrTHaIiIhI1pjsEBERkawx2SEiIiJZY7JjRJs3b0b37t1hb2+P0NBQpKenmzqkR7Zq1So8/fTT6NSpEzw8PDBhwgTk5ubq7RMeHg6FQqH3NXv2bBNF/Gj++te/NupDnz59pO0VFRWIiYlBly5d4OTkhIkTJ6KwsNCEEbde9+7dG/VRoVAgJiYGgPmO45dffolx48ZBo9FAoVDg4MGDetuFEFi+fDm8vb2hUqkQERGBK1eu6O1z9+5dTJ48Gc7OzlCr1Zg5cybu37/fjr1oXnN9rK6uxuLFixEYGAhHR0doNBpMnToVt27d0nuNpsZ/9erV7dyT5rU0ltOnT2/Uh9GjR+vtY85jCaDJv1GFQoE1a9ZI+3T0sTTkvGHIZ2p+fj7Gjh0LBwcHeHh44O2330ZNTU2rYmGyYyR79uzBggULsGLFCmRmZmLAgAEYNWoUioqKTB3aIzl9+jRiYmJw9uxZHDt2DNXV1Rg5ciQePHigt19UVBRu374tfb3//vsmivjRPfnkk3p9+Prrr6Vtb731Fg4dOoRPP/0Up0+fxq1bt/DCCy+YMNrW++abb/T6d+zYMQDAiy++KO1jjuP44MEDDBgwAJs3b25y+/vvv4+NGzdi69atOHfuHBwdHTFq1ChUVFRI+0yePBmXLl3CsWPHcPjwYXz55ZeIjo5ury60qLk+lpWVITMzE8uWLUNmZib279+P3NxcPP/88432fffdd/XGd86cOe0RvsFaGksAGD16tF4fkpKS9Lab81gC0Ovb7du3kZCQAIVCgYkTJ+rt15HH0pDzRkufqVqtFmPHjkVVVRVSU1Oxfft2fPzxx1i+fHnrghFkFCEhISImJkb6WavVCo1GI1atWmXCqNpOUVGRACBOnz4ttQ0bNkzMmzfPdEG1gRUrVogBAwY0ua24uFjY2NiITz/9VGr77rvvBACRlpbWThG2vXnz5olevXoJnU4nhJDHOAIQBw4ckH7W6XTCy8tLrFmzRmorLi4WdnZ2IikpSQghxOXLlwUA8c0330j7fP7550KhUIibN2+2W+yG+nUfm5Keni4AiBs3bkhtfn5+Yt26dcYNrg011c9p06aJ8ePHP/QYOY7l+PHjxfDhw/XazG0sf33eMOQz9d///rewsrISBQUF0j5btmwRzs7OorKy0uD3ZmXHCKqqqpCRkYGIiAipzcrKChEREUhLSzNhZG2npKQEAODq6qrXvnPnTri5ueGpp57CkiVLUFZWZorwHsuVK1eg0WjQs2dPTJ48Gfn5+QCAjIwMVFdX641rnz590K1bN7Md16qqKuzYsQOvvfaa3kNw5TCODeXl5aGgoEBv7FxcXBAaGiqNXVpaGtRqNQYPHiztExERASsrK5w7d67dY24LJSUlUCgUUKvVeu2rV69Gly5dEBQUhDVr1rR6SqAjOHXqFDw8PBAQEIA33ngDd+7ckbbJbSwLCwuRnJyMmTNnNtpmTmP56/OGIZ+paWlpCAwMhKenp7TPqFGjUFpaikuXLhn83nwQqBH8/PPP0Gq1eoMDAJ6envj+++9NFFXb0el0mD9/Pp555hk89dRTUvsrr7wCPz8/aDQaXLx4EYsXL0Zubi72799vwmhbJzQ0FB9//DECAgJw+/ZtrFy5Es8++yxycnJQUFAAW1vbRicOT09PFBQUmCbgx3Tw4EEUFxdj+vTpUpscxvHX6senqb/J+m0FBQXw8PDQ225tbQ1XV1ezHN+KigosXrwYkZGReg9WnDt3LgYNGgRXV1ekpqZiyZIluH37NtauXWvCaFtn9OjReOGFF9CjRw9cu3YNf/nLXzBmzBikpaVBqVTKbiy3b9+OTp06NZoyN6exbOq8YchnakFBQZN/t/XbDMVkh1otJiYGOTk5emtZAOjNhwcGBsLb2xsjRozAtWvX0KtXr/YO85GMGTNG+r5///4IDQ2Fn58f9u7dC5VKZcLIjGPbtm0YM2YMNBqN1CaHcbR01dXVeOmllyCEwJYtW/S2LViwQPq+f//+sLW1xeuvv45Vq1aZzeMIXn75Zen7wMBA9O/fH7169cKpU6cwYsQIE0ZmHAkJCZg8eTLs7e312s1pLB923mgvnMYyAjc3NyiVykYrygsLC+Hl5WWiqNpGbGwsDh8+jJSUFHTt2rXZfUNDQwEAV69ebY/QjEKtVqN37964evUqvLy8UFVVheLiYr19zHVcb9y4gePHj2PWrFnN7ieHcawfn+b+Jr28vBpdQFBTU4O7d++a1fjWJzo3btzAsWPH9Ko6TQkNDUVNTQ2uX7/ePgEaQc+ePeHm5ib9jsplLAHgq6++Qm5ubot/p0DHHcuHnTcM+Uz18vJq8u+2fpuhmOwYga2tLYKDg3HixAmpTafT4cSJEwgLCzNhZI9OCIHY2FgcOHAAJ0+eRI8ePVo8JisrCwDg7e1t5OiM5/79+7h27Rq8vb0RHBwMGxsbvXHNzc1Ffn6+WY5rYmIiPDw8MHbs2Gb3k8M49ujRA15eXnpjV1painPnzkljFxYWhuLiYmRkZEj7nDx5EjqdTkr4Orr6ROfKlSs4fvw4unTp0uIxWVlZsLKyajTtY07++9//4s6dO9LvqBzGst62bdsQHByMAQMGtLhvRxvLls4bhnymhoWFITs7Wy95rU/i+/Xr16pgyAh2794t7OzsxMcffywuX74soqOjhVqt1ltRbk7eeOMN4eLiIk6dOiVu374tfZWVlQkhhLh69ap49913xfnz50VeXp747LPPRM+ePcXQoUNNHHnrLFy4UJw6dUrk5eWJM2fOiIiICOHm5iaKioqEEELMnj1bdOvWTZw8eVKcP39ehIWFibCwMBNH3XparVZ069ZNLF68WK/dnMfx3r174sKFC+LChQsCgFi7dq24cOGCdCXS6tWrhVqtFp999pm4ePGiGD9+vOjRo4coLy+XXmP06NEiKChInDt3Tnz99dfC399fREZGmqpLjTTXx6qqKvH888+Lrl27iqysLL2/0/qrVlJTU8W6detEVlaWuHbtmtixY4dwd3cXU6dONXHP9DXXz3v37olFixaJtLQ0kZeXJ44fPy4GDRok/P39RUVFhfQa5jyW9UpKSoSDg4PYsmVLo+PNYSxbOm8I0fJnak1NjXjqqafEyJEjRVZWljhy5Ihwd3cXS5YsaVUsTHaMaNOmTaJbt27C1tZWhISEiLNnz5o6pEcGoMmvxMREIYQQ+fn5YujQocLV1VXY2dmJJ554Qrz99tuipKTEtIG30qRJk4S3t7ewtbUVPj4+YtKkSeLq1avS9vLycvHmm2+Kzp07CwcHB/GHP/xB3L5924QRP5qjR48KACI3N1ev3ZzHMSUlpcnf0WnTpgkhai8/X7ZsmfD09BR2dnZixIgRjfp/584dERkZKZycnISzs7OYMWOGuHfvngl607Tm+piXl/fQv9OUlBQhhBAZGRkiNDRUuLi4CHt7e9G3b18RFxenlyR0BM31s6ysTIwcOVK4u7sLGxsb4efnJ6Kiohr9R9Kcx7Lehx9+KFQqlSguLm50vDmMZUvnDSEM+0y9fv26GDNmjFCpVMLNzU0sXLhQVFdXtyoWRV1ARERERLLENTtEREQka0x2iIiISNaY7BAREZGsMdkhIiIiWWOyQ0RERLLGZIeIiIhkjckOERERyRqTHSIiIpI1JjtE1OFMnz4dEyZMMNn7T5kyBXFxcQbt+/LLLyM+Pt7IERHR4+AdlImoXSkUima3r1ixAm+99RaEEFCr1e0TVAPffvsthg8fjhs3bsDJyanF/XNycjB06FDk5eXBxcWlHSIkotZiskNE7aqgoED6fs+ePVi+fDlyc3OlNicnJ4OSDGOZNWsWrK2tsXXrVoOPefrppzF9+nTExMQYMTIielScxiKiduXl5SV9ubi4QKFQ6LU5OTk1msYKDw/HnDlzMH/+fHTu3Bmenp746KOP8ODBA8yYMQOdOnXCE088gc8//1zvvXJycjBmzBg4OTnB09MTU6ZMwc8///zQ2LRaLfbt24dx48bptf/zn/+Ev78/7O3t4enpiT/+8Y9628eNG4fdu3c//j8OERkFkx0iMgvbt2+Hm5sb0tPTMWfOHLzxxht48cUX8Zvf/AaZmZkYOXIkpkyZgrKyMgBAcXExhg8fjqCgIJw/fx5HjhxBYWEhXnrppYe+x8WLF1FSUoLBgwdLbefPn8fcuXPx7rvvIjc3F0eOHMHQoUP1jgsJCUF6ejoqKyuN03kieixMdojILAwYMABLly6Fv78/lixZAnt7e7i5uSEqKgr+/v5Yvnw57ty5g4sXLwIAPvjgAwQFBSEuLg59+vRBUFAQEhISkJKSgh9++KHJ97hx4waUSiU8PDyktvz8fDg6OuK5556Dn58fgoKCMHfuXL3jNBoNqqqq9KboiKjjYLJDRGahf//+0vdKpRJdunRBYGCg1Obp6QkAKCoqAlC70DglJUVaA+Tk5IQ+ffoAAK5du9bke5SXl8POzk5vEfXvfvc7+Pn5oWfPnpgyZQp27twpVY/qqVQqAGjUTkQdA5MdIjILNjY2ej8rFAq9tvoERafTAQDu37+PcePGISsrS+/rypUrjaah6rm5uaGsrAxVVVVSW6dOnZCZmYmkpCR4e3tj+fLlGDBgAIqLi6V97t69CwBwd3dvk74SUdtiskNEsjRo0CBcunQJ3bt3xxNPPKH35ejo2OQxAwcOBABcvnxZr93a2hoRERF4//33cfHiRVy/fh0nT56Utufk5KBr165wc3MzWn+I6NEx2SEiWYqJicHdu3cRGRmJb775BteuXcPRo0cxY8YMaLXaJo9xd3fHoEGD8PXXX0tthw8fxsaNG5GVlYUbN27gk08+gU6nQ0BAgLTPV199hZEjRxq9T0T0aJjsEJEsaTQanDlzBlqtFiNHjkRgYCDmz58PtVoNK6uHf/TNmjULO3fulH5Wq9XYv38/hg8fjr59+2Lr1q1ISkrCk08+CQCoqKjAwYMHERUVZfQ+EdGj4U0FiYgaKC8vR0BAAPbs2YOwsLAW99+yZQsOHDiAL774oh2iI6JHwcoOEVEDKpUKn3zySbM3H2zIxsYGmzZtMnJURPQ4WNkhIiIiWWNlh4iIiGSNyQ4RERHJGpMdIiIikjUmO0RERCRrTHaIiIhI1pjsEBERkawx2SEiIiJZY7JDREREssZkh4iIiGTt/wOABEQ5oWGujQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "rates, bins = data_proc.get_count_rate(bin_time=2)\n", "\n", "plt.plot(\n", " bins[:-1],\n", " rates,\n", " label=\"All counts\",\n", ")\n", "\n", "rates, bins = data_proc.get_count_rate(\n", " bin_time=2,\n", " energy_window=(\n", " mean_energy_peak1 - std_energy_peak1 * 2,\n", " mean_energy_peak1 + std_energy_peak1 * 2,\n", " ),\n", ")\n", "\n", "plt.plot(\n", " bins[:-1],\n", " rates,\n", " label=\"peak 1\",\n", ")\n", "\n", "rates, bins = data_proc.get_count_rate(\n", " bin_time=2,\n", " energy_window=(\n", " mean_energy_peak2 - std_energy_peak2 * 2,\n", " mean_energy_peak2 + std_energy_peak2 * 2,\n", " ),\n", ")\n", "\n", "plt.plot(\n", " bins[:-1],\n", " rates,\n", " label=\"peak 2\",\n", ")\n", "\n", "avg_count_rate, error = data_proc.get_avg_rate(0, 100)\n", "plt.hlines(avg_count_rate, 0, 100, color=\"black\", linestyle=\"--\")\n", "\n", "print(f\"Average count rate: {avg_count_rate:.2f} +/- {error:.2f} count/s \")\n", "\n", "plt.legend()\n", "plt.xlabel(\"Time (s)\")\n", "plt.ylabel(\"Count rate (count/s)\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "libra-toolbox", "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.12.5" } }, "nbformat": 4, "nbformat_minor": 2 }