CP9_Geant4_detector_modifications&measurements.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "940f43dc",
   "metadata": {},
   "outputs": [],
   "source": [
    "import os\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5ad0d979",
   "metadata": {},
   "source": [
    "### Particle gun: Electrons with energy 300 MeV"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2c7a04eb",
   "metadata": {},
   "source": [
    "### Lead shield of 2cm thickness and water sphere of 50cm radius"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "9066ce98",
   "metadata": {},
   "outputs": [],
   "source": [
    "shield = pd.read_csv(\"output_nt_Shield.csv\")\n",
    "water = pd.read_csv(\"output_nt_Water.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8389b980",
   "metadata": {},
   "outputs": [],
   "source": [
    "shield_energies = shield[4:].values.astype(\"float\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "ff388024",
   "metadata": {},
   "outputs": [],
   "source": [
    "water_energies = water[4:].values.astype(\"float\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "bef33e40",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total events:  200\n"
     ]
    }
   ],
   "source": [
    "print(\"Total events: \",shield_energies.shape[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "f4e29e48",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHFCAYAAAAaD0bAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABC2ElEQVR4nO3deXQUVf7+8acDSYclCQTJHpKwBNkMCMiiIGEJRkFQEBxUwrAMLiyCDBqXIeg4QUXEBUUdSUQR0GFxAWUREmTAGRAQcJDlawIICQxrACFAcn9/+EuPbVZClg71fp3T51D33qr61KWgn1NV3W0zxhgBAABYiFtlFwAAAFDRCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEBABfnXv/6lu+66Sw0aNJDdbpe/v786deqkxx57zGlceHi4+vTpU+z2UlJSZLPZlJKScsW1XMm6w4YNU3h4+BXvI094eLiGDRtW6vUrWkJCgmw2m1Pbm2++qeTk5HLZn81mU0JCQrlsG0DhCEBABVi2bJk6d+6srKwsvfjii1q5cqVeffVV3XzzzVq4cGGptnnjjTdq48aNuvHGG8u4WmsbOXKkNm7c6NRWngEIQOWoXtkFAFbw4osvKiIiQitWrFD16v/7Z3fvvffqxRdfLNU2vb291bFjx7IqEf9fSEiIQkJCKrsMl3bp0iXZbDancxmoargCBFSA48eP67rrrivwDcPNreB/hl999ZVuvPFG1ahRQ9dff73mzJnj1F/YbazNmzfrzjvvlK+vrzw9PdWmTRt9/PHHJaozOTlZTZs2ld1uV7NmzTR37tySHaB+fVOcPHmyAgICVLNmTd1yyy3697//XeDYzMxMjR49WiEhIfLw8FBERISmTp2qy5cvO8akp6fLZrPpxRdf1PPPP68GDRrI09NT7dq109dff51vm+vXr1ePHj3k5eWlmjVrqnPnzlq2bJnTmF9++UWTJk1SRESEPD095evrq3bt2mn+/PmOMb+/BRYeHq4ffvhBqampstlsstlsTrcEs7KyHNv08PBQcHCwHn30UZ07d85p31lZWRo1apTq1aun2rVr67bbbtOePXtKPL8l3Y/NZtOYMWP0wQcfqFmzZqpZs6aioqL0xRdf5Nvm3r17NWTIEPn5+Tn+zmfNmuU0Ju88++CDD/TYY48pODhYdrtd+/btkyS9++67ioyMlN1uV/PmzfXRRx853TY1xqhJkybq3bt3vv2fPXtWPj4+euSRR0o8D0CZMQDK3ciRI40kM3bsWPPtt9+aixcvFjo2LCzMhISEmObNm5u5c+eaFStWmHvuucdIMqmpqY5xa9euNZLM2rVrHW1r1qwxHh4epkuXLmbhwoXmq6++MsOGDTOSTFJSUpHrJiUlGUmmX79+5vPPPzcffvihady4sQkNDTVhYWHFHmNcXJyx2Wzmz3/+s1m5cqWZMWOGCQ4ONt7e3iYuLs4xLiMjw7HNt99+26xevdo899xzxm63m2HDhjnGpaWlGUkmNDTU3HLLLWbRokXmk08+Me3btzfu7u5mw4YNjrEpKSnG3d3dtG3b1ixcuNAsXbrUxMTEGJvNZhYsWOAYN3r0aFOzZk0zY8YMs3btWvPFF1+YadOmmddff90xZsqUKea3/zVu2bLFNGzY0LRp08Zs3LjRbNy40WzZssUYY8y5c+dM69atzXXXXWdmzJhhVq9ebV599VXj4+NjunfvbnJzc40xxuTm5pro6Ghjt9vN888/b1auXGmmTJliGjZsaCSZKVOmFDm3Jd2PMcZIMuHh4eamm24yH3/8sVm+fLnp1q2bqV69uvm///s/x7gffvjB+Pj4mFatWpm5c+ealStXmscee8y4ubmZhIQEx7i8cyU4ONgMHDjQfPbZZ+aLL74wx48fN2+//baRZAYMGGC++OILM2/ePBMZGWnCwsKczplXX33V2Gw2s2fPHqfjmjVrlpFkfvjhhyKPHygPBCCgAhw7dszccsstRpKRZNzd3U3nzp1NYmKiOXPmjNPYsLAw4+npafbv3+9oO3/+vPH19TWjR492tBUUYq6//nrTpk0bc+nSJadt9unTxwQGBpqcnJwC183JyTFBQUHmxhtvdHozTU9PN+7u7sUGoF27dhlJZsKECU7t8+bNM5KcAtDo0aNN7dq1nY7PGGOmT5/u9GaYF4CCgoLM+fPnHeOysrKMr6+v6dmzp6OtY8eOxs/Pz2kuL1++bFq2bGlCQkIcx9SyZUvTv3//Io/l9wHIGGNatGhhbr311nxjExMTjZubm9m0aZNT+z/+8Q8jySxfvtwYY8yXX35pJJlXX33Vadzzzz9fogBU0v0Y82sA8vf3N1lZWY62zMxM4+bmZhITEx1tvXv3NiEhIeb06dNO2xwzZozx9PQ0J06cMMb871zp2rWr07icnBwTEBBgOnTo4NS+f//+fOdMVlaW8fLyMuPHj3ca27x5cxMdHV3ksQPlhVtgQAWoV6+evvnmG23atEnTpk1Tv379tGfPHsXHx6tVq1Y6duyY0/jWrVurQYMGjmVPT09FRkZq//79he5j3759+vHHH3XfffdJki5fvux43X777crIyNDu3bsLXHf37t06fPiwhgwZ4nT7JywsTJ07dy72+NauXStJjn3nGTRoUL7bfl988YWio6MVFBTkVGNsbKwkKTU11Wn83XffLU9PT8eyl5eX+vbtq3Xr1iknJ0fnzp3Tv/71Lw0cOFC1a9d2jKtWrZoeeOAB/fzzz47jvummm/Tll1/qiSeeUEpKis6fP1/ssRXliy++UMuWLdW6dWunY+ndu7fT7cnC5mfIkCFlup880dHR8vLyciz7+/vLz8/Pcf5cuHBBX3/9te666y7VrFkz37ly4cIFffvtt07bHDBggNPy7t27lZmZqUGDBjm1N2jQQDfffLNTm5eXl/74xz8qOTnZcctuzZo1+s9//qMxY8aUaA6AskYAAipQu3bt9Pjjj+uTTz7R4cOHNWHCBKWnp+d7ELpevXr51rXb7UW+YR85ckSSNGnSJLm7uzu9Hn74YUnKF7TyHD9+XJIUEBCQr6+gtpKuX7169XzHcuTIEX3++ef5amzRokWBNRZW08WLF3X27FmdPHlSxhgFBgbmGxcUFORU32uvvabHH39cS5cuVXR0tHx9fdW/f3/t3bu32GMsyJEjR7R9+/Z8x+Ll5SVjjONYjh8/XuBclGRur2Q/eYo7f44fP67Lly/r9ddfz7fN22+/XVL+v4ffz2/enPr7++fbV0FtY8eO1ZkzZzRv3jxJ0htvvKGQkBD169evRHMAlDUe4Qcqibu7u6ZMmaJXXnlFO3fuvOrtXXfddZKk+Ph43X333QWOadq0aYHteW+YmZmZ+foKaitq/eDgYEf75cuXHW+Uv63zhhtu0PPPP1/gtvJCS1H7z8zMlIeHh2rXrq3q1avLzc1NGRkZ+cYdPnzYsU9JqlWrlqZOnaqpU6fqyJEjjqtBffv21Y8//ljscf7eddddpxo1auR7QP23/dKv85M3F78NJyWZ2yvZT0nVrVvXcYWssAeQIyIinJZ//91IeceRF7x/q6Djaty4sWJjYzVr1izFxsbqs88+09SpU1WtWrUrqh0oKwQgoAJkZGQUeIVi165dkvK/6ZdG06ZN1aRJE33//ff629/+dsXrBgYGav78+Zo4caLjzW7//v3asGFDsfV169ZNkjRv3jy1bdvW0f7xxx87fbJLkvr06aPly5erUaNGqlu3brG1LV68WC+99JLjNtiZM2f0+eefq0uXLqpWrZpq1aqlDh06aPHixZo+fbpq1KghScrNzdWHH36okJAQRUZG5tuuv7+/hg0bpu+//14zZ87UL7/8opo1axZYQ2FX3/r06aO//e1vqlevXr7A8FvR0dF68cUXNW/ePI0bN87R/tFHHxV7/Feyn5KqWbOmoqOjtXXrVt1www3y8PC44m00bdpUAQEB+vjjjzVx4kRH+4EDBwo9Z8aPH6+YmBjFxcWpWrVqGjVq1FUdB3A1CEBABejdu7dCQkLUt29fXX/99crNzdW2bdv08ssvq3bt2ho/fnyZ7Oftt99WbGysevfurWHDhik4OFgnTpzQrl27tGXLFn3yyScFrufm5qbnnntOI0eO1F133aVRo0bp1KlTSkhIKNFtmmbNmun+++/XzJkz5e7urp49e2rnzp2aPn26vL29ncY+++yzWrVqlTp37qxx48apadOmunDhgtLT07V8+XLNnj3b6Xt4qlWrpl69emnixInKzc3VCy+8oKysLE2dOtUxJjExUb169VJ0dLQmTZokDw8Pvfnmm9q5c6fmz5/vCHQdOnRQnz59dMMNN6hu3bratWuXPvjgA3Xq1KnQ8CNJrVq10oIFC7Rw4UI1bNhQnp6eatWqlR599FEtWrRIXbt21YQJE3TDDTcoNzdXBw4c0MqVK/XYY4+pQ4cOiomJUdeuXTV58mSdO3dO7dq10z//+U998MEHxc6tpBLv50q8+uqruuWWW9SlSxc99NBDCg8P15kzZ7Rv3z59/vnnWrNmTZHru7m5aerUqRo9erQGDhyo4cOH69SpU5o6daoCAwML/HqHXr16qXnz5lq7dq3uv/9++fn5XVHNQJmq5IewAUtYuHChGTJkiGnSpImpXbu2cXd3Nw0aNDAPPPCA+c9//uM0NiwszNxxxx35tnHrrbc6fRKpoE+BGWPM999/bwYNGmT8/PyMu7u7CQgIMN27dzezZ88udt2///3vpkmTJsbDw8NERkaaOXPmmLi4uBJ9DD47O9s89thjxs/Pz3h6epqOHTuajRs3mrCwMKdPgRljzH//+18zbtw4ExERYdzd3Y2vr69p27ateeqpp8zZs2eNMf/7FNgLL7xgpk6dakJCQoyHh4dp06aNWbFiRb79f/PNN6Z79+6mVq1apkaNGqZjx47m888/dxrzxBNPmHbt2pm6desau91uGjZsaCZMmGCOHTvmGFPQp8DS09NNTEyM8fLyMpKc5uPs2bPm6aefNk2bNjUeHh6Oj5ZPmDDBZGZmOsadOnXKDB8+3NSpU8fUrFnT9OrVy/z4448l+hTYlexHknnkkUfyrV/Q30NaWpoZPny4CQ4ONu7u7qZ+/fqmc+fO5q9//atjTN658sknnxRY1zvvvGMaN27sdM7069fPtGnTpsDxCQkJRpL59ttviz1moDzZjDGm0tIXABQiPT1dEREReumllzRp0qTKLgcldOrUKUVGRqp///5655138vW3a9dONptNmzZtqoTqgP/hFhgAoFQyMzP1/PPPKzo6WvXq1dP+/fv1yiuv6MyZM063dbOysrRz50598cUX+u6777RkyZJKrBr4FQEIAFAqdrtd6enpevjhh3XixAnVrFlTHTt21OzZsx1fayBJW7ZscYSkKVOmqH///pVXNPD/cQsMAABYDl+ECAAALIcABAAALIcABAAALIeHoAuQm5urw4cPy8vLK9/XvwMAANdkjNGZM2cUFBRU4Jdx/hYBqACHDx9WaGhoZZcBAABK4eDBg07fKF8QAlABvLy8JP06gb//Gn8AAOCasrKyFBoa6ngfLwoBqAB5t728vb0JQAAAVDEleXyFh6ABAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlVK/sAgAA1hP+xLJix6RPu6MCKoFVcQUIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYDgEIAABYTqUGoHXr1qlv374KCgqSzWbT0qVLnfptNluBr5deeqnQbSYnJxe4zoULF8r5aAAAQFVRqQHo3LlzioqK0htvvFFgf0ZGhtNrzpw5stlsGjBgQJHb9fb2zreup6dneRwCAACogqpX5s5jY2MVGxtbaH9AQIDT8qeffqro6Gg1bNiwyO3abLZ86wIAAOSpMs8AHTlyRMuWLdOIESOKHXv27FmFhYUpJCREffr00datWyugQgAAUFVUmQD0/vvvy8vLS3fffXeR466//nolJyfrs88+0/z58+Xp6ambb75Ze/fuLXSd7OxsZWVlOb0AAMC1q8oEoDlz5ui+++4r9lmejh076v7771dUVJS6dOmijz/+WJGRkXr99dcLXScxMVE+Pj6OV2hoaFmXDwAAXEiVCEDffPONdu/erZEjR17xum5ubmrfvn2RV4Di4+N1+vRpx+vgwYNXUy4AAHBxlfoQdEm99957atu2raKioq54XWOMtm3bplatWhU6xm63y263X02JAACgCqnUAHT27Fnt27fPsZyWlqZt27bJ19dXDRo0kCRlZWXpk08+0csvv1zgNoYOHarg4GAlJiZKkqZOnaqOHTuqSZMmysrK0muvvaZt27Zp1qxZ5X9AAACgSqjUALR582ZFR0c7lidOnChJiouLU3JysiRpwYIFMsboD3/4Q4HbOHDggNzc/ncn79SpU/rTn/6kzMxM+fj4qE2bNlq3bp1uuumm8jsQAABQpdiMMaayi3A1WVlZ8vHx0enTp+Xt7V3Z5QDANSf8iWXFjkmfdkcFVIJryZW8f1eJh6ABAADKEgEIAABYDgEIAABYTpX4GDwAAK6AZ5euHVwBAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAllO9sgsAAKC8hT+xrNgx6dPuqIBK4Cq4AgQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACynUgPQunXr1LdvXwUFBclms2np0qVO/cOGDZPNZnN6dezYsdjtLlq0SM2bN5fdblfz5s21ZMmScjoCAABQFVVqADp37pyioqL0xhtvFDrmtttuU0ZGhuO1fPnyIre5ceNGDR48WA888IC+//57PfDAAxo0aJD+9a9/lXX5AACgiqpemTuPjY1VbGxskWPsdrsCAgJKvM2ZM2eqV69eio+PlyTFx8crNTVVM2fO1Pz586+qXgAAcG1w+WeAUlJS5Ofnp8jISI0aNUpHjx4tcvzGjRsVExPj1Na7d29t2LCh0HWys7OVlZXl9AIAANculw5AsbGxmjdvntasWaOXX35ZmzZtUvfu3ZWdnV3oOpmZmfL393dq8/f3V2ZmZqHrJCYmysfHx/EKDQ0ts2MAAACup1JvgRVn8ODBjj+3bNlS7dq1U1hYmJYtW6a777670PVsNpvTsjEmX9tvxcfHa+LEiY7lrKwsQhAAANcwlw5AvxcYGKiwsDDt3bu30DEBAQH5rvYcPXo031Wh37Lb7bLb7WVWJwAAcG0ufQvs944fP66DBw8qMDCw0DGdOnXSqlWrnNpWrlypzp07l3d5AACgiqjUK0Bnz57Vvn37HMtpaWnatm2bfH195evrq4SEBA0YMECBgYFKT0/Xk08+qeuuu0533XWXY52hQ4cqODhYiYmJkqTx48era9eueuGFF9SvXz99+umnWr16tdavX1/hxwcAAFxTpQagzZs3Kzo62rGc9xxOXFyc3nrrLe3YsUNz587VqVOnFBgYqOjoaC1cuFBeXl6OdQ4cOCA3t/9dyOrcubMWLFigp59+Ws8884waNWqkhQsXqkOHDhV3YAAAwKVVagDq1q2bjDGF9q9YsaLYbaSkpORrGzhwoAYOHHg1pQEAgGtYlXoGCAAAoCwQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOVU6m+BAQBwrQl/YlmF7i992h0Vur9rBVeAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5VSv7AIAACit8CeWVXYJqKK4AgQAACyHAAQAACyHAAQAACyHAAQAACynUgPQunXr1LdvXwUFBclms2np0qWOvkuXLunxxx9Xq1atVKtWLQUFBWno0KE6fPhwkdtMTk6WzWbL97pw4UI5Hw0AAKgqKjUAnTt3TlFRUXrjjTfy9f3yyy/asmWLnnnmGW3ZskWLFy/Wnj17dOeddxa7XW9vb2VkZDi9PD09y+MQAABAFVSpH4OPjY1VbGxsgX0+Pj5atWqVU9vrr7+um266SQcOHFCDBg0K3a7NZlNAQECZ1goAAK4dVeoZoNOnT8tms6lOnTpFjjt79qzCwsIUEhKiPn36aOvWrUWOz87OVlZWltMLAABcu6pMALpw4YKeeOIJDRkyRN7e3oWOu/7665WcnKzPPvtM8+fPl6enp26++Wbt3bu30HUSExPl4+PjeIWGhpbHIQAAABdRJQLQpUuXdO+99yo3N1dvvvlmkWM7duyo+++/X1FRUerSpYs+/vhjRUZG6vXXXy90nfj4eJ0+fdrxOnjwYFkfAgAAcCEu/1MYly5d0qBBg5SWlqY1a9YUefWnIG5ubmrfvn2RV4DsdrvsdvvVlgoAAKoIl74ClBd+9u7dq9WrV6tevXpXvA1jjLZt26bAwMByqBAAAFRFlXoF6OzZs9q3b59jOS0tTdu2bZOvr6+CgoI0cOBAbdmyRV988YVycnKUmZkpSfL19ZWHh4ckaejQoQoODlZiYqIkaerUqerYsaOaNGmirKwsvfbaa9q2bZtmzZpV8QcIAABcUqUGoM2bNys6OtqxPHHiRElSXFycEhIS9Nlnn0mSWrdu7bTe2rVr1a1bN0nSgQMH5Ob2vwtZp06d0p/+9CdlZmbKx8dHbdq00bp163TTTTeV78EAAIAqo1IDULdu3WSMKbS/qL48KSkpTsuvvPKKXnnllastDQAAXMNc+hkgAACA8kAAAgAAluPyH4MHALiG8CeWlWhc+rQ7KnR/QGlwBQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFhOqQLQli1btGPHDsfyp59+qv79++vJJ5/UxYsXy6w4AACA8lCqADR69Gjt2bNHkvTTTz/p3nvvVc2aNfXJJ59o8uTJZVogAABAWStVANqzZ49at24tSfrkk0/UtWtXffTRR0pOTtaiRYvKsj4AAIAyV700KxljlJubK0lavXq1+vTpI0kKDQ3VsWPHyq46AECFCH9iWWWXAFSoUl0Bateunf7617/qgw8+UGpqqu644w5JUlpamvz9/cu0QAAAgLJWqgD0yiuvaMuWLRozZoyeeuopNW7cWJL0j3/8Q507dy7TAgEAAMpaqW6BRUVFOX0KLM9LL72k6tVLtUkAAIAKU6orQA0bNtTx48fztV+4cEGRkZFXXRQAAEB5KlUASk9PV05OTr727Oxs/fzzz1ddFAAAQHm6ovtVn332mePPK1askI+Pj2M5JydHX3/9tSIiIsquOgAAgHJwRQGof//+kiSbzaa4uDinPnd3d4WHh+vll18us+IAAADKwxUFoLzv/omIiNCmTZt03XXXlUtRAAAA5alUH9lKS0sr6zoAAAAqTKk/s/7111/r66+/1tGjRx1XhvLMmTPnqgsDAAAoL6UKQFOnTtWzzz6rdu3aKTAwUDabrazrAgAAKDelCkCzZ89WcnKyHnjggbKuBwAAoNyV6nuALl68yE9eAACAKqtUAWjkyJH66KOPrnrn69atU9++fRUUFCSbzaalS5c69RtjlJCQoKCgINWoUUPdunXTDz/8UOx2Fy1apObNm8tut6t58+ZasmTJVdcKAACuHaW6BXbhwgW98847Wr16tW644Qa5u7s79c+YMaNE2zl37pyioqL0xz/+UQMGDMjX/+KLL2rGjBlKTk5WZGSk/vrXv6pXr17avXu3vLy8Ctzmxo0bNXjwYD333HO66667tGTJEg0aNEjr169Xhw4drvxgAQDANadUAWj79u1q3bq1JGnnzp1OfVfyQHRsbKxiY2ML7DPGaObMmXrqqad09913S5Lef/99+fv766OPPtLo0aMLXG/mzJnq1auX4uPjJUnx8fFKTU3VzJkzNX/+/BLXBgAArl2lCkBr164t6zrySUtLU2ZmpmJiYhxtdrtdt956qzZs2FBoANq4caMmTJjg1Na7d2/NnDmz0H1lZ2crOzvbsZyVlXV1xQMAAJdWqmeAKkJmZqYkyd/f36nd39/f0VfYele6TmJionx8fByv0NDQq6gcAAC4ulJdAYqOji7yVteaNWtKXdDv/X4/xphib7Nd6Trx8fGaOHGiYzkrK4sQBADANaxUASjv+Z88ly5d0rZt27Rz5858P5JaWgEBAZJ+vaITGBjoaD969Gi+Kzy/X+/3V3uKW8dut8tut19lxQAAoKooVQB65ZVXCmxPSEjQ2bNnr6qgPBEREQoICNCqVavUpk0bSb9+/1BqaqpeeOGFQtfr1KmTVq1a5fQc0MqVK/neIgAA4FDq3wIryP3336+bbrpJ06dPL9H4s2fPat++fY7ltLQ0bdu2Tb6+vmrQoIEeffRR/e1vf1OTJk3UpEkT/e1vf1PNmjU1ZMgQxzpDhw5VcHCwEhMTJUnjx49X165d9cILL6hfv3769NNPtXr1aq1fv74sDxUAAFRhZRqANm7cKE9PzxKP37x5s6Kjox3Lec/hxMXFKTk5WZMnT9b58+f18MMP6+TJk+rQoYNWrlzp9B1ABw4ckJvb/57l7ty5sxYsWKCnn35azzzzjBo1aqSFCxfyHUAAAMChVAEo73t58hhjlJGRoc2bN+uZZ54p8Xa6desmY0yh/TabTQkJCUpISCh0TEpKSr62gQMHauDAgSWuAwAAWEupApCPj4/Tspubm5o2bapnn33W6Xt7AAAAXFGpAlBSUlJZ1wEAAFBhruoZoO+++067du2SzWZT8+bNHZ/WAgAAcGWlCkBHjx7Vvffeq5SUFNWpU0fGGJ0+fVrR0dFasGCB6tevX9Z1AgAAlJlS/RTG2LFjlZWVpR9++EEnTpzQyZMntXPnTmVlZWncuHFlXSMAAECZKtUVoK+++kqrV69Ws2bNHG3NmzfXrFmzeAgaAAC4vFJdAcrNzZW7u3u+dnd3d+Xm5l51UQAAAOWpVAGoe/fuGj9+vA4fPuxoO3TokCZMmKAePXqUWXEAAADloVQB6I033tCZM2cUHh6uRo0aqXHjxoqIiNCZM2f0+uuvl3WNAAAAZapUzwCFhoZqy5YtWrVqlX788UcZY9S8eXP17NmzrOsDAAAoc1d0BWjNmjVq3ry5srKyJEm9evXS2LFjNW7cOLVv314tWrTQN998Uy6FAgAAlJUrCkAzZ87UqFGj5O3tna/Px8dHo0eP1owZM8qsOAAAgPJwRbfAvv/+e73wwguF9sfExGj69OlXXRQAABUt/IlllV0CKtAVXQE6cuRIgR9/z1O9enX997//veqiAAAAytMVBaDg4GDt2LGj0P7t27crMDDwqosCAAAoT1cUgG6//Xb95S9/0YULF/L1nT9/XlOmTFGfPn3KrDgAAIDycEXPAD399NNavHixIiMjNWbMGDVt2lQ2m027du3SrFmzlJOTo6eeeqq8agUAACgTVxSA/P39tWHDBj300EOKj4+XMUaSZLPZ1Lt3b7355pvy9/cvl0IBAADKyhV/EWJYWJiWL1+ukydPat++fTLGqEmTJqpbt2551AcAAFDmSvVN0JJUt25dtW/fvixrAQAAqBCl+i0wAACAqowABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALMflA1B4eLhsNlu+1yOPPFLg+JSUlALH//jjjxVcOQAAcFXVK7uA4mzatEk5OTmO5Z07d6pXr1665557ilxv9+7d8vb2dizXr1+/3GoEAABVi8sHoN8Hl2nTpqlRo0a69dZbi1zPz89PderUKcfKAABAVeXyt8B+6+LFi/rwww81fPhw2Wy2Ise2adNGgYGB6tGjh9auXVtBFQIAgKrA5a8A/dbSpUt16tQpDRs2rNAxgYGBeuedd9S2bVtlZ2frgw8+UI8ePZSSkqKuXbsWuE52drays7Mdy1lZWWVdOgAAcCFVKgC99957io2NVVBQUKFjmjZtqqZNmzqWO3XqpIMHD2r69OmFBqDExERNnTq1zOsFAACuqcrcAtu/f79Wr16tkSNHXvG6HTt21N69ewvtj4+P1+nTpx2vgwcPXk2pAADAxVWZK0BJSUny8/PTHXfcccXrbt26VYGBgYX22+122e32qykPAABUIVUiAOXm5iopKUlxcXGqXt255Pj4eB06dEhz586VJM2cOVPh4eFq0aKF46HpRYsWadGiRZVROgAAcEFVIgCtXr1aBw4c0PDhw/P1ZWRk6MCBA47lixcvatKkSTp06JBq1KihFi1aaNmyZbr99tsrsmQAAODCqkQAiomJkTGmwL7k5GSn5cmTJ2vy5MkVUBUAAKiqqsxD0AAAAGWFAAQAACyHAAQAACynSjwDBAAovfAnll3T+wNKgytAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAcghAAADAclw6ACUkJMhmszm9AgICilwnNTVVbdu2laenpxo2bKjZs2dXULUAAKCqqF7ZBRSnRYsWWr16tWO5WrVqhY5NS0vT7bffrlGjRunDDz/UP//5Tz388MOqX7++BgwYUBHlAgCAKsDlA1D16tWLveqTZ/bs2WrQoIFmzpwpSWrWrJk2b96s6dOnE4AAAICDS98Ck6S9e/cqKChIERERuvfee/XTTz8VOnbjxo2KiYlxauvdu7c2b96sS5cuFbpedna2srKynF4AAODa5dIBqEOHDpo7d65WrFihd999V5mZmercubOOHz9e4PjMzEz5+/s7tfn7++vy5cs6duxYoftJTEyUj4+P4xUaGlqmxwEAAFyLSweg2NhYDRgwQK1atVLPnj21bNkySdL7779f6Do2m81p2RhTYPtvxcfH6/Tp047XwYMHy6B6AADgqlz+GaDfqlWrllq1aqW9e/cW2B8QEKDMzEyntqNHj6p69eqqV69eodu12+2y2+1lWisAAHBdLn0F6Peys7O1a9cuBQYGFtjfqVMnrVq1yqlt5cqVateundzd3SuiRAAAUAW4dACaNGmSUlNTlZaWpn/9618aOHCgsrKyFBcXJ+nXW1dDhw51jH/wwQe1f/9+TZw4Ubt27dKcOXP03nvvadKkSZV1CAAAwAW59C2wn3/+WX/4wx907Ngx1a9fXx07dtS3336rsLAwSVJGRoYOHDjgGB8REaHly5drwoQJmjVrloKCgvTaa6/xEXgAAODEZvKeEoZDVlaWfHx8dPr0aXl7e1d2OQBwVcKfWFbZJaAcpU+7o7JLcBlX8v7t0rfAAAAAygMBCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWI5L/xo8AFyLyvLHSfkhTJTV+WS1c4krQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHKqV3YBAFDewp9YVibbSZ92R5lspyyV1bEBVsMVIAAAYDkEIAAAYDkEIAAAYDkEIAAAYDkuHYASExPVvn17eXl5yc/PT/3799fu3buLXCclJUU2my3f68cff6ygqgEAgKtz6QCUmpqqRx55RN9++61WrVqly5cvKyYmRufOnSt23d27dysjI8PxatKkSQVUDAAAqgKX/hj8V1995bSclJQkPz8/fffdd+ratWuR6/r5+alOnTrlWB0AAKiqXPoK0O+dPn1akuTr61vs2DZt2igwMFA9evTQ2rVrixybnZ2trKwspxcAALh2VZkAZIzRxIkTdcstt6hly5aFjgsMDNQ777yjRYsWafHixWratKl69OihdevWFbpOYmKifHx8HK/Q0NDyOAQAAOAiXPoW2G+NGTNG27dv1/r164sc17RpUzVt2tSx3KlTJx08eFDTp08v9LZZfHy8Jk6c6FjOysoiBAEAcA2rEleAxo4dq88++0xr165VSEjIFa/fsWNH7d27t9B+u90ub29vpxcAALh2ufQVIGOMxo4dqyVLliglJUURERGl2s7WrVsVGBhYxtUBAICqyqUD0COPPKKPPvpIn376qby8vJSZmSlJ8vHxUY0aNST9evvq0KFDmjt3riRp5syZCg8PV4sWLXTx4kV9+OGHWrRokRYtWlRpxwEAAFyLSwegt956S5LUrVs3p/akpCQNGzZMkpSRkaEDBw44+i5evKhJkybp0KFDqlGjhlq0aKFly5bp9ttvr6iyAQCAi3PpAGSMKXZMcnKy0/LkyZM1efLkcqoIAABcC6rEQ9AAAABliQAEAAAshwAEAAAsx6WfAQJQ9YQ/sazYMenT7qiw7VS0ktQNoPJxBQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFgOAQgAAFhO9couAAULf2JZsWPSp91RAZVUfczltaskf7cAyk5J/81Vhf9TuQIEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAsp0oEoDfffFMRERHy9PRU27Zt9c033xQ5PjU1VW3btpWnp6caNmyo2bNnV1ClAACgKnD5ALRw4UI9+uijeuqpp7R161Z16dJFsbGxOnDgQIHj09LSdPvtt6tLly7aunWrnnzySY0bN06LFi2q4MoBAICrcvkANGPGDI0YMUIjR45Us2bNNHPmTIWGhuqtt94qcPzs2bPVoEEDzZw5U82aNdPIkSM1fPhwTZ8+vYIrBwAArsqlA9DFixf13XffKSYmxqk9JiZGGzZsKHCdjRs35hvfu3dvbd68WZcuXSq3WgEAQNVRvbILKMqxY8eUk5Mjf39/p3Z/f39lZmYWuE5mZmaB4y9fvqxjx44pMDAw3zrZ2dnKzs52LJ8+fVqSlJWVdbWHUGq52b8UO6Yy66tKmMuKVVbzXZLtVLSqWjdQEmV5flfW/6l5+zXGFDvWpQNQHpvN5rRsjMnXVtz4gtrzJCYmaurUqfnaQ0NDr7TUCuUzs7IruHYwlxWrqs53Va0bKImyPL8r+9/KmTNn5OPjU+QYlw5A1113napVq5bvas/Ro0fzXeXJExAQUOD46tWrq169egWuEx8fr4kTJzqWc3NzdeLECdWrV6/IoFUaWVlZCg0N1cGDB+Xt7V2m28b/MM8Vg3muGMxzxWGuK0Z5zbMxRmfOnFFQUFCxY106AHl4eKht27ZatWqV7rrrLkf7qlWr1K9fvwLX6dSpkz7//HOntpUrV6pdu3Zyd3cvcB273S673e7UVqdOnasrvhje3t7846oAzHPFYJ4rBvNccZjrilEe81zclZ88Lv0QtCRNnDhRf//73zVnzhzt2rVLEyZM0IEDB/Tggw9K+vXqzdChQx3jH3zwQe3fv18TJ07Url27NGfOHL333nuaNGlSZR0CAABwMS59BUiSBg8erOPHj+vZZ59VRkaGWrZsqeXLlyssLEySlJGR4fSdQBEREVq+fLkmTJigWbNmKSgoSK+99poGDBhQWYcAAABcjMsHIEl6+OGH9fDDDxfYl5ycnK/t1ltv1ZYtW8q5qtKx2+2aMmVKvltuKFvMc8VgnisG81xxmOuK4QrzbDMl+awYAADANcTlnwECAAAoawQgAABgOQQgAABgOQQgAABgOQSgcpCQkCCbzeb0CggIcPQbY5SQkKCgoCDVqFFD3bp10w8//FCJFVcN69atU9++fRUUFCSbzaalS5c69ZdkXrOzszV27Fhdd911qlWrlu688079/PPPFXgUrq+4eR42bFi+87tjx45OY5jn4iUmJqp9+/by8vKSn5+f+vfvr927dzuN4Zy+eiWZZ87pq/fWW2/phhtucHyxYadOnfTll186+l3xXCYAlZMWLVooIyPD8dqxY4ej78UXX9SMGTP0xhtvaNOmTQoICFCvXr105syZSqzY9Z07d05RUVF64403Cuwvybw++uijWrJkiRYsWKD169fr7Nmz6tOnj3JycirqMFxecfMsSbfddpvT+b18+XKnfua5eKmpqXrkkUf07bffatWqVbp8+bJiYmJ07tw5xxjO6atXknmWOKevVkhIiKZNm6bNmzdr8+bN6t69u/r16+cIOS55LhuUuSlTppioqKgC+3Jzc01AQICZNm2ao+3ChQvGx8fHzJ49u4IqrPokmSVLljiWSzKvp06dMu7u7mbBggWOMYcOHTJubm7mq6++qrDaq5Lfz7MxxsTFxZl+/foVug7zXDpHjx41kkxqaqoxhnO6vPx+no3hnC4vdevWNX//+99d9lzmClA52bt3r4KCghQREaF7771XP/30kyQpLS1NmZmZiomJcYy12+269dZbtWHDhsoqt8orybx+9913unTpktOYoKAgtWzZkrm/QikpKfLz81NkZKRGjRqlo0ePOvqY59I5ffq0JMnX11cS53R5+f085+GcLjs5OTlasGCBzp07p06dOrnsuUwAKgcdOnTQ3LlztWLFCr377rvKzMxU586ddfz4cccv1f/+1+z9/f3z/Yo9Sq4k85qZmSkPDw/VrVu30DEoXmxsrObNm6c1a9bo5Zdf1qZNm9S9e3dlZ2dLYp5LwxijiRMn6pZbblHLli0lcU6Xh4LmWeKcLis7duxQ7dq1Zbfb9eCDD2rJkiVq3ry5y57LVeKnMKqa2NhYx59btWqlTp06qVGjRnr//fcdD9bZbDandYwx+dpw5Uozr8z9lRk8eLDjzy1btlS7du0UFhamZcuW6e677y50Pea5cGPGjNH27du1fv36fH2c02WnsHnmnC4bTZs21bZt23Tq1CktWrRIcXFxSk1NdfS72rnMFaAKUKtWLbVq1Up79+51fBrs94n26NGj+dIxSq4k8xoQEKCLFy/q5MmThY7BlQsMDFRYWJj27t0riXm+UmPHjtVnn32mtWvXKiQkxNHOOV22CpvngnBOl46Hh4caN26sdu3aKTExUVFRUXr11Vdd9lwmAFWA7Oxs7dq1S4GBgYqIiFBAQIBWrVrl6L948aJSU1PVuXPnSqyyaivJvLZt21bu7u5OYzIyMrRz507m/iocP35cBw8eVGBgoCTmuaSMMRozZowWL16sNWvWKCIiwqmfc7psFDfPBeGcLhvGGGVnZ7vuuVwuj1Zb3GOPPWZSUlLMTz/9ZL799lvTp08f4+XlZdLT040xxkybNs34+PiYxYsXmx07dpg//OEPJjAw0GRlZVVy5a7tzJkzZuvWrWbr1q1GkpkxY4bZunWr2b9/vzGmZPP64IMPmpCQELN69WqzZcsW0717dxMVFWUuX75cWYflcoqa5zNnzpjHHnvMbNiwwaSlpZm1a9eaTp06meDgYOb5Cj300EPGx8fHpKSkmIyMDMfrl19+cYzhnL56xc0z53TZiI+PN+vWrTNpaWlm+/bt5sknnzRubm5m5cqVxhjXPJcJQOVg8ODBJjAw0Li7u5ugoCBz9913mx9++MHRn5uba6ZMmWICAgKM3W43Xbt2NTt27KjEiquGtWvXGkn5XnFxccaYks3r+fPnzZgxY4yvr6+pUaOG6dOnjzlw4EAlHI3rKmqef/nlFxMTE2Pq169v3N3dTYMGDUxcXFy+OWSei1fQHEsySUlJjjGc01evuHnmnC4bw4cPN2FhYcbDw8PUr1/f9OjRwxF+jHHNc9lmjDHlc20JAADANfEMEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEIByYbPZtHTp0jLfbnp6umw2m7Zt21bm23Z17733nmJiYiq7DE2aNEnjxo2r7DKAq0IAAqqIYcOGyWazyWazyd3dXf7+/urVq5fmzJmj3Nzcyi4vn4yMDMXGxkqydmgpK9nZ2frLX/6iZ555xtGWkJAgm82m2267Ld/4F198UTabTd26dSvR9o8cOSJ3d3d9+OGHBfaPHj1aN9xwgyRp8uTJSkpKUlpa2pUfCOAiCEBAFXLbbbcpIyND6enp+vLLLxUdHa3x48erT58+unz5cmWX5yQgIEB2u72yy6hQly5dKrdtL1q0SLVr11aXLl2c2gMDA7V27Vr9/PPPTu1JSUlq0KBBibfv7++vO+64Q0lJSfn6zp8/rwULFmjEiBGSJD8/P8XExGj27NmlOBLANRCAgCrEbrcrICBAwcHBuvHGG/Xkk0/q008/1Zdffqnk5GTHuNOnT+tPf/qT/Pz85O3tre7du+v777939CckJKh169Z6++23FRoaqpo1a+qee+7RqVOnHGNyc3P17LPPKiQkRHa7Xa1bt9ZXX33l6L948aLGjBmjwMBAeXp6Kjw8XImJiY7+394Cy/sF7jZt2uS7KpGUlKRmzZrJ09NT119/vd58802nY/73v/+tNm3ayNPTU+3atdPWrVuLnaeLFy9q8uTJCg4OVq1atdShQwelpKQ4+pOTk1WnTh2tWLFCzZo1U+3atR3h8reKqi3vqtbHH3+sbt26ydPTUx9++KEuX76scePGqU6dOqpXr54ef/xxxcXFqX///pKkuXPnql69esrOznba14ABAzR06NBCj2nBggW6884787XnhZH333/f0bZhwwYdO3ZMd9xxR77xRR3TiBEjtHbtWqWnpzut849//EMXLlzQ/fff72i78847NX/+/ELrBVxeuf3KGIAyFRcXZ/r161dgX1RUlImNjTXG/PqjgzfffLPp27ev2bRpk9mzZ4957LHHTL169czx48eNMcZMmTLF1KpVy3Tv3t1s3brVpKammsaNG5shQ4Y4tjljxgzj7e1t5s+fb3788UczefJk4+7ubvbs2WOMMeall14yoaGhZt26dSY9Pd1888035qOPPnKsL8ksWbLEGGPMv//9byPJrF692mRkZDjqeOedd0xgYKBZtGiR+emnn8yiRYuMr6+vSU5ONsYYc/bsWVO/fn0zePBgs3PnTvP555+bhg0bGklm69athc7VkCFDTOfOnc26devMvn37zEsvvWTsdruj9qSkJOPu7m569uxpNm3aZL777jvTrFkzp+Mvrra0tDQjyYSHhzvGHDp0yPz1r381vr6+ZvHixWbXrl3mwQcfNN7e3o6/u19++cX4+PiYjz/+2LGv//73v8bDw8OsWbOm0GOqU6eOWbBggVPblClTTFRUlFm8eLFp3Lixo33EiBFm/PjxZvz48ebWW28t8TFdvnzZBAYGmilTpjjtp1u3bmbQoEFObf/5z3+MJJOenl5ozYArIwABVURRAWjw4MGmWbNmxhhjvv76a+Pt7W0uXLjgNKZRo0bm7bffNsb8+sZZrVo1c/DgQUf/l19+adzc3ExGRoYxxpigoCDz/PPPO22jffv25uGHHzbGGDN27FjTvXt3k5ubW2BNvw1AeWHh96ElNDTUKTQZY8xzzz1nOnXqZIwx5u233za+vr7m3Llzjv633nqryAC0b98+Y7PZzKFDh5zae/ToYeLj440xvwYgSWbfvn2O/lmzZhl/f/8S15Z3TDNnznQa4+/vb1566SXH8uXLl02DBg2c/u4eeughR2A1xpiZM2eahg0bFjqXJ0+eNJLMunXrnNrzAtDFixeNn5+fSU1NNWfPnjVeXl7m+++/zxeAijsmY4x5/PHHTVhYmKOWn376ydhsNrNixQqn9U6fPm0kmZSUlAJrBlxd9cq57gSgLBljZLPZJEnfffedzp49q3r16jmNOX/+vP7v//7PsdygQQOFhIQ4ljt16qTc3Fzt3r1bNWvW1OHDh3XzzTc7bePmm2923EobNmyYevXqpaZNm+q2225Tnz59rugTSv/973918OBBjRgxQqNGjXK0X758WT4+PpKkXbt2KSoqSjVr1nSqsyhbtmyRMUaRkZFO7dnZ2U5zUrNmTTVq1MixHBgYqKNHj5a4tjzt2rVz/Pn06dM6cuSIbrrpJkdbtWrV1LZtW6cH1UeNGqX27dvr0KFDCg4OVlJSkuMh94KcP39ekuTp6Vlgv7u7u+6//34lJSXpp59+UmRkpOOB5TwlPaYRI0bohRde0Jo1a9SjRw/NmTNHISEh6tmzp9P2atSoIUn65ZdfCqwJcHUEIOAasGvXLsdzNrm5uQoMDHR65iVPnTp1Ct1G3pvvb9+Ef/+G/NugdeONNyotLU1ffvmlVq9erUGDBqlnz576xz/+UaKa8wLBu+++qw4dOjj1VatWzbG/K5Wbm6tq1arpu+++c2wnT+3atR1/dnd3d+qz2WyO/ZWktjy1atXKV0NB8/Zbbdq0UVRUlObOnavevXtrx44d+vzzzws9pnr16slms+nkyZOFjhk+fLg6dOignTt3avjw4fn6S3pMTZo0UZcuXZSUlKTo6Gi9//77+uMf/yg3N+dHRk+cOCFJql+/fqE1Aa6MAARUcWvWrNGOHTs0YcIESb8Gk8zMTFWvXl3h4eGFrnfgwAEdPnxYQUFBkqSNGzfKzc1NkZGR8vb2VlBQkNavX6+uXbs61tmwYYPT1Q1vb28NHjxYgwcP1sCBA3XbbbfpxIkT8vX1ddqXh4eHJCknJ8fR5u/vr+DgYP3000+67777CqyxefPm+uCDD3T+/HnHFYdvv/22yPlo06aNcnJydPTo0XyfmCqpktRWEB8fH/n7++vf//63Y985OTnaunWrWrdu7TR25MiReuWVV3To0CH17NlToaGhhW7Xw8NDzZs313/+859Cr7K1aNFCLVq00Pbt2zVkyJCrOqYRI0booYceUr9+/fTzzz/rj3/8Y74xO3fulLu7u1q0aFHktgBXRQACqpDs7GxlZmYqJydHR44c0VdffaXExET16dPH8Qminj17qlOnTurfv79eeOEFNW3aVIcPH9by5cvVv39/xy0bT09PxcXFafr06crKytK4ceM0aNAgBQQESJL+/Oc/a8qUKWrUqJFat26tpKQkbdu2TfPmzZMkvfLKKwoMDFTr1q3l5uamTz75RAEBAQVeZfLz81ONGjX01VdfKSQkRJ6envLx8VFCQoLGjRsnb29vxcbGKjs7W5s3b9bJkyc1ceJEDRkyRE899ZRGjBihp59+Wunp6Zo+fXqRcxQZGan77rtPQ4cO1csvv6w2bdro2LFjWrNmjVq1aqXbb7+9RHNdXG2FGTt2rBITE9W4cWNdf/31ev3113Xy5Ml8V4Xuu+8+TZo0Se+++67mzp1bbD29e/fW+vXr9eijjxY6Zs2aNbp06VKhV/pKekz33HOPxo0bp9GjR6tHjx4FBulvvvlGXbp0cQRToMqpzAeQAJRcXFyckWQkmerVq5v69eubnj17mjlz5picnBynsVlZWWbs2LEmKCjIuLu7m9DQUHPfffeZAwcOGGP+9/Dsm2++aYKCgoynp6e5++67zYkTJxzbyMnJMVOnTjXBwcHG3d3dREVFmS+//NLR/84775jWrVubWrVqGW9vb9OjRw+zZcsWR79+8xC0Mca8++67JjQ01Li5uTk9mDtv3jzTunVr4+HhYerWrWu6du1qFi9e7OjfuHGjiYqKMh4eHqZ169Zm0aJFxX4K7OLFi+Yvf/mLCQ8PN+7u7iYgIMDcddddZvv27caYXx+C9vHxcVpnyZIl5vf/JRZVW2EPdl+6dMmMGTPGeHt7m7p165rHH3/c3HPPPebee+/NV+cDDzxgfH198z2wXpBdu3aZGjVqmFOnTjna8v4eC/P7h6CLO6bf+tOf/mQk5XtoOk9kZKSZP39+sXUDrspmTClusgOo0hISErR06VK+mbkC5ObmqlmzZho0aJCee+45p75evXqpWbNmeu2110q0rUGDBqlNmzaKj48vj1JLbNmyZfrzn/+s7du3q3p1biSgauKLEAGgDO3fv1/vvvuu9uzZox07duihhx5SWlqa03M5J06c0IIFC7RmzRo98sgjJd72Sy+95PQgd2U5d+6ckpKSCD+o0jh7AaAMubm5KTk5WZMmTZIxRi1bttTq1avVrFkzx5gbb7xRJ0+edDyjVVJhYWEaO3ZseZR9RQYNGlTZJQBXjVtgAADAcrgFBgAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALIcABAAALOf/AUdzO7ioryGjAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(shield_energies,bins=np.linspace(50,300,50))\n",
    "plt.title(\"Shield deposited energy\")\n",
    "plt.xlabel(\"Deposited energy (MeV)\")\n",
    "plt.ylabel(\"Counts\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "7456659e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHFCAYAAAAaD0bAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABHxUlEQVR4nO3dd3gU5f7//9eGJBvAJBAkjQQSEJBmQECKUkIJBFBQFDyogKKi0pGDYDkQy4moBzkqivqlWMGCYAEVkATkgEekKCBSPiQEMQFpCTWB5P794S97WNJDyi7zfFzXXhdzzz2z73tm476cuXfXZowxAgAAsBCPyi4AAACgohGAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAAACA5RCAgArw6aefymaz6aOPPsqzLioqSjabTd9++22edQ0aNND1119foud6/fXXtWDBgtKWWqYiIiI0fPjwyi6j2KZPny6bzebUVp7H02azafr06eWybwCFIwABFaBr166y2WxKSEhwaj927Ji2bdum6tWr51n3+++/a9++fYqOji7Rc7lSAHI3999/vzZs2ODUxvEErkyelV0AYAVXX321mjdvrsTERKf2NWvWyNPTUyNGjMgTgHKXSxqAyoMxRufOnVPVqlUru5RyFRYWprCwsMouw6WdP39eNptNnp68fcC9cQUIqCDR0dHatWuXUlNTHW2JiYlq27at+vTpo02bNunkyZNO66pUqaJOnTpJkuLi4tSuXTsFBATIz89P119/vebOnauLf884IiJCO3bs0Jo1a2Sz2WSz2RQREeFYn5GRoUmTJikyMlLe3t6qU6eOxo8fr9OnTzvVarPZNHr0aM2ZM0dNmjSR3W7XO++8U+DYzp8/r8mTJys4OFjVqlXTTTfdpB9//DHfvmlpaRo5cqTCwsLk7e2tyMhIxcXF6cKFC44+ycnJstlseuGFF/Tcc8+pbt268vHxUZs2bfTdd9/l2ee6devUvXt3+fr6qlq1aurYsaOWLVvm1OfMmTOOsfv4+CggIEBt2rTRwoULHX0uvQVWVsczIyNDDzzwgGrVqqWrrrpKvXv31u7duws8npcq6Xl777331KRJE1WrVk1RUVH66quv8uxzz549GjJkiAIDA2W329WkSRPNnj3bqU9iYqJsNpvee+89Pfroo6pTp47sdrv27t0rSXr77bfVqFEj2e12NW3aVB9++KGGDx/uOEbGGDVs2FC9evXK8/ynTp2Sv7+/Ro0aVezjAJQpA6BCLFmyxEgyH374oaOtRYsWZurUqebkyZPG09PTLFu2zLEuMjLStG3b1rE8fPhwM3fuXLNy5UqzcuVK88wzz5iqVauauLg4R5/Nmzeb+vXrm1atWpkNGzaYDRs2mM2bNxtjjDl9+rRp2bKlufrqq83MmTPNqlWrzL///W/j7+9vunXrZnJychz7kWTq1KljrrvuOvPhhx+a1atXm+3btxc4tmHDhhmbzWb+/ve/mxUrVpiZM2eaOnXqGD8/PzNs2DBHv9TUVBMeHm7q1atn3nzzTbNq1SrzzDPPGLvdboYPH+7ol5SUZCSZ8PBwc9NNN5nFixebTz75xLRt29Z4eXmZ9evXO/omJiYaLy8v07p1a/PRRx+ZpUuXmpiYGGOz2cyiRYsc/UaOHGmqVatmZs6caRISEsxXX31lnn/+efPqq686+kybNs1c/J/FsjieOTk5Jjo62tjtdvPcc8+ZFStWmGnTppn69esbSWbatGkFHtfSnLeIiAhzww03mI8//tgsX77cdO3a1Xh6epr/+7//c/TbsWOH8ff3Ny1atDDvvvuuWbFihXn00UeNh4eHmT59uqNfQkKC47Vw++23my+++MJ89dVX5ujRo+bNN980kszAgQPNV199ZT744APTqFEjU69ePVOvXj3HPv79738bm81mdu/e7TSu2bNnG0lmx44dhY4fKC8EIKCCHDt2zHh4eJgHH3zQGGPMkSNHjM1mM998840xxpgbbrjBTJo0yRhjTEpKipFkJk+enO++srOzzfnz583TTz9tatWq5fQm2KxZM9OlS5c828THxxsPDw+zceNGp/ZPP/3USDLLly93tEky/v7+5tixY0WOa+fOnUaSmTBhglP7Bx98YCQ5BaCRI0eaq666yuzfv9+p70svveT0ZpgbgEJDQ83Zs2cd/TIyMkxAQIDp0aOHo619+/YmMDDQnDx50tF24cIF07x5cxMWFuY4Ns2bNzcDBgwodCyXBiBjLv94fv3110aS+fe//+3U77nnnitWACrpeQsKCjIZGRmOtrS0NOPh4WHi4+Mdbb169TJhYWEmPT3daZ+jR482Pj4+jvOeG4A6d+7s1C87O9sEBwebdu3aObXv37/feHl5OQWgjIwM4+vra8aNG+fUt2nTpiY6OrrQsQPliVtgQAWpWbOmoqKiHPOA1qxZoypVqujGG2+UJHXp0sUx7ye/+T+rV69Wjx495O/vrypVqsjLy0v/+Mc/dPToUR0+fLjI5//qq6/UvHlztWzZUhcuXHA8evXqJZvNlmd+Urdu3VSzZs0i95tb61133eXUPmjQoDzzRL766itFR0crNDTUqYbY2FjHMbnYbbfdJh8fH8eyr6+vbr75Zq1du1bZ2dk6ffq0/vvf/+r222/XVVdd5ehXpUoV3XPPPfr999+1a9cuSdINN9ygr7/+WlOmTFFiYqLOnj1b5NgKU9zjWdDxGTJkSJk+T67o6Gj5+vo6loOCghQYGKj9+/dLks6dO6fvvvtOt956q6pVq+a0zz59+ujcuXP64YcfnPY5cOBAp+Vdu3YpLS1NgwYNcmqvW7eu4/Wcy9fXV/fee68WLFjguGW3evVq/frrrxo9enSxjgFQHghAQAWKjo7W7t279ccffyghIUGtW7d2vHF36dJFW7ZsUXp6uhISEuTp6ambbrpJkvTjjz8qJiZG0l/zLv7zn/9o48aNeuKJJySpWG/mhw4d0i+//CIvLy+nh6+vr4wxOnLkiFP/kJCQYo3p6NGjkqTg4GCndk9PT9WqVStPDV9++WWeGpo1ayZJeWq4dJ+5bVlZWTp16pSOHz8uY0y+tYaGhjrV98orr+ixxx7T0qVLFR0drYCAAA0YMEB79uwp1jgvVdzjefTo0XyPRX5ju5znyXXp80iS3W53vEaOHj2qCxcu6NVXX82zzz59+kjKex4uPb65xzQoKCjPc+XXNmbMGJ08eVIffPCBJOm1115TWFiY+vfvX6xjAJQHpvEDFSg6OlozZ85UYmKiEhMTHW84khxhZ+3atY7J0bnhaNGiRfLy8tJXX33ldEVk6dKlxX7uq6++WlWrVtW8efMKXH+xS78PpyC5b7hpaWmqU6eOo/3ChQuON8qLn+O6667Tc889l+++ckNLrrS0tDx90tLS5O3trauuukqenp7y8PBwmlie648//nA8pyRVr15dcXFxiouL06FDhxxXg26++Wb99ttvxRrrpWMpzvGsVauW41hcHE7yG9vlPE9x1axZ03GFrKAJyJGRkU7Ll74Wcsdx6NChPNvmN65rrrlGsbGxmj17tmJjY/XFF18oLi5OVapUKVHtQJmq3DtwgLWkp6ebKlWqmFtvvdXYbDan+RvGGHP99debgQMHGknm8ccfd7RPnDjRXHXVVSYrK8vRdubMGVO3bl0jySQlJTnt44Ybbsjz3M8++6ypVq2a2bdvX5F1SjKjRo0q1ph+/fXXYs8Buv/++01oaGiRc4uKmgPUvXt3R1uHDh1McHCwOXPmjKMtOzvbtGjRwmkOUH7Gjx9vJJnTp08bY/KfA3S5x/Ny5wCVxXmrV6+e03no0aOHiYqKMpmZmYXuL3cO0CeffOLUXpI5QLlWrFhhJJno6Gjj7e1tDh06VOR4gPLEFSCgAuV+fH3p0qXy8PDIM1+iS5cumjVrliTn+T99+/bVzJkzNWTIED344IM6evSoXnrpJdnt9jzP0aJFCy1atEgfffSR6tevLx8fH7Vo0ULjx4/X4sWL1blzZ02YMEHXXXedcnJylJKSohUrVujRRx9Vu3btSjymJk2a6O6779asWbPk5eWlHj16aPv27XrppZfk5+fn1Pfpp5/WypUr1bFjR40dO1aNGzfWuXPnlJycrOXLl2vOnDlO38NTpUoV9ezZUxMnTlROTo5mzJihjIwMxcXFOfrEx8erZ8+eio6O1qRJk+Tt7a3XX39d27dv18KFCx1XL9q1a6d+/frpuuuuU82aNbVz506999576tChg6pVq1bg+C73eMbExKhz586aPHmyTp8+rTZt2ug///mP3nvvvWId3/I4b//+97910003qVOnTnr44YcVERGhkydPau/evfryyy+1evXqQrf38PBQXFycRo4cqdtvv1333XefTpw4obi4OIWEhMjDI+/sip49e6pp06ZKSEjQ3XffrcDAwBLVDJS5yk5ggNVMnjzZSDJt2rTJs27p0qVGkvH29nZclcg1b94807hxY2O32039+vVNfHy8mTt3bp4rQMnJySYmJsb4+voaSU7/N37q1Cnz5JNPmsaNGxtvb2/HR6EnTJhg0tLSHP1UgitAxhiTmZlpHn30URMYGGh8fHxM+/btzYYNG/JceTDGmD///NOMHTvWREZGGi8vLxMQEGBat25tnnjiCXPq1CljzP+uAM2YMcPExcWZsLAw4+3tbVq1amW+/fbbPM///fffm27dupnq1aubqlWrmvbt25svv/zSqc+UKVNMmzZtTM2aNR3HcMKECebIkSOOPvldASqL43nixAlz3333mRo1aphq1aqZnj17mt9++61YV4BK8jwFnbf8zkNSUpK57777TJ06dYyXl5epXbu26dixo3n22WcdfQq6ApTrrbfeMtdcc43x9vY2jRo1MvPmzTP9+/c3rVq1yrf/9OnTjSTzww8/FDlmoLzZjLnoW9QAwAUkJycrMjJSL774oiZNmlTZ5aCYTpw4oUaNGmnAgAF666238qxv06aNbDabNm7cWAnVAc64BQYAKLG0tDQ999xzio6OVq1atbR//369/PLLOnnypMaNG+fol5GRoe3bt+urr77Spk2btGTJkkqsGvgfAhAAoMTsdruSk5P1yCOP6NixY6pWrZrat2+vOXPmOL7WQJI2b97sCEnTpk3TgAEDKq9o4CLcAgMAAJbDFyECAADLIQABAADLIQABAADLYRJ0PnJycvTHH3/I19e32D8HAAAAKpcxRidPnlRoaGi+X8h5MQJQPv744w+Fh4dXdhkAAKAUDhw44PSt8vkhAOXD19dX0l8H8NKv8gcAAK4pIyND4eHhjvfxwhCA8pF728vPz48ABACAmynO9BUmQQMAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMup1AAUHx+vtm3bytfXV4GBgRowYIB27drl1McYo+nTpys0NFRVq1ZV165dtWPHjiL3vXjxYjVt2lR2u11NmzbVkiVLymsYAADAzVRqAFqzZo1GjRqlH374QStXrtSFCxcUExOj06dPO/q88MILmjlzpl577TVt3LhRwcHB6tmzp06ePFngfjds2KDBgwfrnnvu0c8//6x77rlHgwYN0n//+9+KGBYAAHBxNmOMqewicv35558KDAzUmjVr1LlzZxljFBoaqvHjx+uxxx6TJGVmZiooKEgzZszQyJEj893P4MGDlZGRoa+//trR1rt3b9WsWVMLFy4sso6MjAz5+/srPT2dH0MFAMBNlOT926XmAKWnp0uSAgICJElJSUlKS0tTTEyMo4/dbleXLl20fv36AvezYcMGp20kqVevXoVuAwAArMOzsgvIZYzRxIkTddNNN6l58+aSpLS0NElSUFCQU9+goCDt37+/wH2lpaXlu03u/i6VmZmpzMxMx3JGRkapxgAAANyDywSg0aNH65dfftG6devyrLPZbE7Lxpg8bZezTXx8vOLi4kpYMa4kEVOWFdkn+fm+FVAJAKAiuMQtsDFjxuiLL75QQkKCwsLCHO3BwcGSlOfKzeHDh/Nc4blYcHBwibaZOnWq0tPTHY8DBw6UdigAAMANVGoAMsZo9OjR+uyzz7R69WpFRkY6rY+MjFRwcLBWrlzpaMvKytKaNWvUsWPHAvfboUMHp20kacWKFQVuY7fb5efn5/QAAABXrkq9BTZq1Ch9+OGH+vzzz+Xr6+u4auPv76+qVavKZrNp/Pjx+uc//6mGDRuqYcOG+uc//6lq1appyJAhjv0MHTpUderUUXx8vCRp3Lhx6ty5s2bMmKH+/fvr888/16pVq/K9vQYAAKynUgPQG2+8IUnq2rWrU/v8+fM1fPhwSdLkyZN19uxZPfLIIzp+/LjatWunFStWyNfX19E/JSVFHh7/u5jVsWNHLVq0SE8++aSeeuopNWjQQB999JHatWtX7mMCAACuz6W+B8hV8D1A1sMkaABwf277PUAAAAAVgQAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAsp1ID0Nq1a3XzzTcrNDRUNptNS5cudVpvs9nyfbz44osF7nPBggX5bnPu3LlyHg0AAHAXlRqATp8+raioKL322mv5rk9NTXV6zJs3TzabTQMHDix0v35+fnm29fHxKY8hAAAAN+RZmU8eGxur2NjYAtcHBwc7LX/++eeKjo5W/fr1C92vzWbLsy0AAEAut5kDdOjQIS1btkwjRowosu+pU6dUr149hYWFqV+/ftqyZUsFVAgAANyF2wSgd955R76+vrrtttsK7XfttddqwYIF+uKLL7Rw4UL5+Pjoxhtv1J49ewrcJjMzUxkZGU4PAABw5XKbADRv3jzdddddRc7lad++ve6++25FRUWpU6dO+vjjj9WoUSO9+uqrBW4THx8vf39/xyM8PLysywcAAC7ELQLQ999/r127dun+++8v8bYeHh5q27ZtoVeApk6dqvT0dMfjwIEDl1MuAABwcZU6Cbq45s6dq9atWysqKqrE2xpjtHXrVrVo0aLAPna7XXa7/XJKBAAAbqRSA9CpU6e0d+9ex3JSUpK2bt2qgIAA1a1bV5KUkZGhTz75RP/617/y3cfQoUNVp04dxcfHS5Li4uLUvn17NWzYUBkZGXrllVe0detWzZ49u/wHBAAA3EKlBqCffvpJ0dHRjuWJEydKkoYNG6YFCxZIkhYtWiRjjP72t7/lu4+UlBR5ePzvTt6JEyf04IMPKi0tTf7+/mrVqpXWrl2rG264ofwGAgAA3IrNGGMquwhXk5GRIX9/f6Wnp8vPz6+yy0EFiJiyrMg+yc/3rYBKAAClVZL3b7eYBA0AAFCWCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByKjUArV27VjfffLNCQ0Nls9m0dOlSp/XDhw+XzWZzerRv377I/S5evFhNmzaV3W5X06ZNtWTJknIaAQAAcEeVGoBOnz6tqKgovfbaawX26d27t1JTUx2P5cuXF7rPDRs2aPDgwbrnnnv0888/65577tGgQYP03//+t6zLBwAAbsqzMp88NjZWsbGxhfax2+0KDg4u9j5nzZqlnj17aurUqZKkqVOnas2aNZo1a5YWLlx4WfUCAIArg8vPAUpMTFRgYKAaNWqkBx54QIcPHy60/4YNGxQTE+PU1qtXL61fv748ywQAAG6kUq8AFSU2NlZ33HGH6tWrp6SkJD311FPq1q2bNm3aJLvdnu82aWlpCgoKcmoLCgpSWlpagc+TmZmpzMxMx3JGRkbZDAAAALgklw5AgwcPdvy7efPmatOmjerVq6dly5bptttuK3A7m83mtGyMydN2sfj4eMXFxV1+wVe4iCnLiuyT/HzfMtlPcRXn+QAAuJTL3wK7WEhIiOrVq6c9e/YU2Cc4ODjP1Z7Dhw/nuSp0salTpyo9Pd3xOHDgQJnVDAAAXI9bBaCjR4/qwIEDCgkJKbBPhw4dtHLlSqe2FStWqGPHjgVuY7fb5efn5/QAAABXrkq9BXbq1Cnt3bvXsZyUlKStW7cqICBAAQEBmj59ugYOHKiQkBAlJyfr8ccf19VXX61bb73Vsc3QoUNVp04dxcfHS5LGjRunzp07a8aMGerfv78+//xzrVq1SuvWravw8QEAANdUqQHop59+UnR0tGN54sSJkqRhw4bpjTfe0LZt2/Tuu+/qxIkTCgkJUXR0tD766CP5+vo6tklJSZGHx/8uZHXs2FGLFi3Sk08+qaeeekoNGjTQRx99pHbt2lXcwAAAgEur1ADUtWtXGWMKXP/tt98WuY/ExMQ8bbfffrtuv/32yykNAABcwdxqDhAAAEBZIAABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLqdTfAgMuV8SUZUX2SX6+bwVUAgBwJ1wBAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAluNZ2QUA7iJiyrIi+yQ/37cCKgEAXC6uAAEAAMshAAEAAMshAAEAAMup1AC0du1a3XzzzQoNDZXNZtPSpUsd686fP6/HHntMLVq0UPXq1RUaGqqhQ4fqjz/+KHSfCxYskM1my/M4d+5cOY8GAAC4i0oNQKdPn1ZUVJRee+21POvOnDmjzZs366mnntLmzZv12Wefaffu3brllluK3K+fn59SU1OdHj4+PuUxBAAA4IYq9VNgsbGxio2NzXedv7+/Vq5c6dT26quv6oYbblBKSorq1q1b4H5tNpuCg4PLtFYAAHDlcKs5QOnp6bLZbKpRo0ah/U6dOqV69eopLCxM/fr105YtWyqmQAAA4BbcJgCdO3dOU6ZM0ZAhQ+Tn51dgv2uvvVYLFizQF198oYULF8rHx0c33nij9uzZU+A2mZmZysjIcHoAAIArl1sEoPPnz+vOO+9UTk6OXn/99UL7tm/fXnfffbeioqLUqVMnffzxx2rUqJFeffXVAreJj4+Xv7+/4xEeHl7WQwAAAC7E5QPQ+fPnNWjQICUlJWnlypWFXv3Jj4eHh9q2bVvoFaCpU6cqPT3d8Thw4MDllg0AAFyYS/8URm742bNnjxISElSrVq0S78MYo61bt6pFixYF9rHb7bLb7ZdTKgAAcCOVGoBOnTqlvXv3OpaTkpK0detWBQQEKDQ0VLfffrs2b96sr776StnZ2UpLS5MkBQQEyNvbW5I0dOhQ1alTR/Hx8ZKkuLg4tW/fXg0bNlRGRoZeeeUVbd26VbNnz674AQIAAJdUqQHop59+UnR0tGN54sSJkqRhw4Zp+vTp+uKLLyRJLVu2dNouISFBXbt2lSSlpKTIw+N/d/JOnDihBx98UGlpafL391erVq20du1a3XDDDeU7GAAA4DYqNQB17dpVxpgC1xe2LldiYqLT8ssvv6yXX375cksDAABXMJefBA0AAFDWCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByShWANm/erG3btjmWP//8cw0YMECPP/64srKyyqw4AACA8lCqADRy5Ejt3r1bkrRv3z7deeedqlatmj755BNNnjy5TAsEAAAoa6UKQLt373b8PMUnn3yizp0768MPP9SCBQu0ePHisqwPAACgzJUqABljlJOTI0latWqV+vTpI0kKDw/XkSNHyq46AACAclCqANSmTRs9++yzeu+997RmzRr17dtX0l+/5h4UFFSmBQIAAJS1UgWgl19+WZs3b9bo0aP1xBNP6JprrpEkffrpp+rYsWOZFggAAFDWSvVr8FFRUU6fAsv14osvytOzUn9gHgAAoEilugJUv359HT16NE/7uXPn1KhRo8suCgAAoDyVKgAlJycrOzs7T3tmZqZ+//33yy4KAACgPJXoftUXX3zh+Pe3334rf39/x3J2dra+++47RUZGll11AAAA5aBEAWjAgAGSJJvNpmHDhjmt8/LyUkREhP71r3+VWXEAAADloUQBKPe7fyIjI7Vx40ZdffXV5VIUAABAeSrVR7aSkpLKug4AAIAKU+rPrH/33Xf67rvvdPjwYceVoVzz5s277MLgniKmLKvsEipVccef/Hzfcq4EAFCYUgWguLg4Pf3002rTpo1CQkJks9nKui4AAIByU6oANGfOHC1YsED33HNPWdcDAABQ7kr1PUBZWVn85AUAAHBbpQpA999/vz788MOyrgUAAKBClOoW2Llz5/TWW29p1apVuu666+Tl5eW0fubMmWVSHAAAQHkoVQD65Zdf1LJlS0nS9u3bndYxIRoAALi6UgWghISEsq4DAACgwpRqDhAAAIA7K9UVoOjo6EJvda1evbrUBQEAAJS3UgWg3Pk/uc6fP6+tW7dq+/bteX4kFQAAwNWUKgC9/PLL+bZPnz5dp06duqyCAAAAyluZzgG6++67+R0wAADg8so0AG3YsEE+Pj7F7r927VrdfPPNCg0Nlc1m09KlS53WG2M0ffp0hYaGqmrVquratat27NhR5H4XL16spk2bym63q2nTplqyZElJhwIAAK5gpboFdttttzktG2OUmpqqn376SU899VSx93P69GlFRUXp3nvv1cCBA/Osf+GFFzRz5kwtWLBAjRo10rPPPquePXtq165d8vX1zXefGzZs0ODBg/XMM8/o1ltv1ZIlSzRo0CCtW7dO7dq1K9lAAQDAFclmjDEl3ejee+91Wvbw8FDt2rXVrVs3xcTElK4Qm01LlizRgAEDJP0VqkJDQzV+/Hg99thjkqTMzEwFBQVpxowZGjlyZL77GTx4sDIyMvT111872nr37q2aNWtq4cKFxaolIyND/v7+Sk9Pl5+fX6nGcyWKmLKssksoleTn+xbZp6LHVpyaAAAlU5L371JdAZo/f36pCiuJpKQkpaWlOQUqu92uLl26aP369QUGoA0bNmjChAlObb169dKsWbPKs1wAAOBGShWAcm3atEk7d+6UzWZT06ZN1apVq7KqS2lpaZKkoKAgp/agoCDt37+/0O3y2yZ3f/nJzMxUZmamYzkjI6M0JQMAADdRqgB0+PBh3XnnnUpMTFSNGjVkjFF6erqio6O1aNEi1a5du8wKvPQLF40xRf7eWEm3iY+PV1xcXOmLhEtzxVt3xamJ22QAUH5K9SmwMWPGKCMjQzt27NCxY8d0/Phxbd++XRkZGRo7dmyZFBYcHCxJea7cHD58OM8Vnku3K+k2U6dOVXp6uuNx4MCBy6gcAAC4ulIFoG+++UZvvPGGmjRp4mhr2rSpZs+e7TT5+HJERkYqODhYK1eudLRlZWVpzZo16tixY4HbdejQwWkbSVqxYkWh29jtdvn5+Tk9AADAlatUt8BycnLk5eWVp93Ly0s5OTnF3s+pU6e0d+9ex3JSUpK2bt2qgIAA1a1bV+PHj9c///lPNWzYUA0bNtQ///lPVatWTUOGDHFsM3ToUNWpU0fx8fGSpHHjxqlz586aMWOG+vfvr88//1yrVq3SunXrSjNUAABwBSpVAOrWrZvGjRunhQsXKjQ0VJJ08OBBTZgwQd27dy/2fn766SdFR0c7lidOnChJGjZsmBYsWKDJkyfr7NmzeuSRR3T8+HG1a9dOK1ascPoOoJSUFHl4/O9CVseOHbVo0SI9+eSTeuqpp9SgQQN99NFHfAcQAABwKNX3AB04cED9+/fX9u3bFR4eLpvNppSUFLVo0UKff/65wsLCyqPWCsP3AOXPFScTX8mYBA0AJVPu3wMUHh6uzZs3a+XKlfrtt99kjFHTpk3Vo0ePUhUMAABQkUo0CXr16tVq2rSp43tyevbsqTFjxmjs2LFq27atmjVrpu+//75cCgUAACgrJQpAs2bN0gMPPJDvZSV/f3+NHDlSM2fOLLPiAAAAykOJAtDPP/+s3r17F7g+JiZGmzZtuuyiAAAAylOJAtChQ4fy/fh7Lk9PT/3555+XXRQAAEB5KlEAqlOnjrZt21bg+l9++UUhISGXXRQAAEB5KlEA6tOnj/7xj3/o3LlzedadPXtW06ZNU79+/cqsOAAAgPJQoo/BP/nkk/rss8/UqFEjjR49Wo0bN5bNZtPOnTs1e/ZsZWdn64knniivWgEAAMpEiQJQUFCQ1q9fr4cfflhTp05V7nco2mw29erVS6+//nqhPzoKAADgCkr8RYj16tXT8uXLdfz4ce3du1fGGDVs2FA1a9Ysj/oAAADKXKm+CVqSatasqbZt25ZlLQAAABWi1AEI7qG4v9/F704BAKykRJ8CAwAAuBIQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOUQgAAAgOV4VnYBAPIXMWVZkX2Sn+9bAZUAwJWHK0AAAMByCEAAAMByCEAAAMByXD4ARUREyGaz5XmMGjUq3/6JiYn59v/tt98quHIAAOCqXH4S9MaNG5Wdne1Y3r59u3r27Kk77rij0O127dolPz8/x3Lt2rXLrUYAAOBeXD4AXRpcnn/+eTVo0EBdunQpdLvAwEDVqFGjHCsDAADuyuVvgV0sKytL77//vu677z7ZbLZC+7Zq1UohISHq3r27EhISKqhCAADgDlz+CtDFli5dqhMnTmj48OEF9gkJCdFbb72l1q1bKzMzU++99566d++uxMREde7cOd9tMjMzlZmZ6VjOyMgo69IBAIALcasANHfuXMXGxio0NLTAPo0bN1bjxo0dyx06dNCBAwf00ksvFRiA4uPjFRcXV+b1AgAA1+Q2t8D279+vVatW6f777y/xtu3bt9eePXsKXD916lSlp6c7HgcOHLicUgEAgItzmytA8+fPV2BgoPr2LflX/2/ZskUhISEFrrfb7bLb7ZdTHgAAcCNuEYBycnI0f/58DRs2TJ6eziVPnTpVBw8e1LvvvitJmjVrliIiItSsWTPHpOnFixdr8eLFlVE6AABwQW4RgFatWqWUlBTdd999edalpqYqJSXFsZyVlaVJkybp4MGDqlq1qpo1a6Zly5apT58+FVkyAABwYW4RgGJiYmSMyXfdggULnJYnT56syZMnV0BVAADAXbnNJGgAAICyQgACAACWQwACAACW4xZzgFD+IqYsq+wSUArFOW/Jz5f8qyMA4ErHFSAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5npVdAIDyFTFlWZF9kp/vWwGVAIDr4AoQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHJcOQNOnT5fNZnN6BAcHF7rNmjVr1Lp1a/n4+Kh+/fqaM2dOBVULAADchcv/GnyzZs20atUqx3KVKlUK7JuUlKQ+ffrogQce0Pvvv6///Oc/euSRR1S7dm0NHDiwIsoFAABuwOUDkKenZ5FXfXLNmTNHdevW1axZsyRJTZo00U8//aSXXnqJAAQAABxc+haYJO3Zs0ehoaGKjIzUnXfeqX379hXYd8OGDYqJiXFq69Wrl3766SedP3++vEsFAABuwqUDULt27fTuu+/q22+/1dtvv620tDR17NhRR48ezbd/WlqagoKCnNqCgoJ04cIFHTlypMDnyczMVEZGhtMDAABcuVw6AMXGxmrgwIFq0aKFevTooWXLlkmS3nnnnQK3sdlsTsvGmHzbLxYfHy9/f3/HIzw8vAyqBwAArsqlA9ClqlevrhYtWmjPnj35rg8ODlZaWppT2+HDh+Xp6alatWoVuN+pU6cqPT3d8Thw4ECZ1g0AAFyLy0+CvlhmZqZ27typTp065bu+Q4cO+vLLL53aVqxYoTZt2sjLy6vA/drtdtnt9jKtFQAAuC6XvgI0adIkrVmzRklJSfrvf/+r22+/XRkZGRo2bJikv67cDB061NH/oYce0v79+zVx4kTt3LlT8+bN09y5czVp0qTKGgIAAHBBLn0F6Pfff9ff/vY3HTlyRLVr11b79u31ww8/qF69epKk1NRUpaSkOPpHRkZq+fLlmjBhgmbPnq3Q0FC98sorfAQeAAA4sZncWcJwyMjIkL+/v9LT0+Xn51fZ5VyWiCnLKrsEuIHk5/tWdgkAcNlK8v7t0rfAAAAAygMBCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWI5L/xgqANdSkb8tx++TAShPXAECAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACWQwACAACW41nZBQCofBFTllV2CQBQobgCBAAALIcABAAALIcABAAALMelA1B8fLzatm0rX19fBQYGasCAAdq1a1eh2yQmJspms+V5/PbbbxVUNQAAcHUuHYDWrFmjUaNG6YcfftDKlSt14cIFxcTE6PTp00Vuu2vXLqWmpjoeDRs2rICKAQCAO3DpT4F98803Tsvz589XYGCgNm3apM6dOxe6bWBgoGrUqFGO1QEAAHfl0leALpWeni5JCggIKLJvq1atFBISou7duyshIaG8SwMAAG7Epa8AXcwYo4kTJ+qmm25S8+bNC+wXEhKit956S61bt1ZmZqbee+89de/eXYmJiQVeNcrMzFRmZqZjOSMjo8zrBwAArsNtAtDo0aP1yy+/aN26dYX2a9y4sRo3buxY7tChgw4cOKCXXnqpwAAUHx+vuLi4Mq0XAAC4Lre4BTZmzBh98cUXSkhIUFhYWIm3b9++vfbs2VPg+qlTpyo9Pd3xOHDgwOWUCwAAXJxLXwEyxmjMmDFasmSJEhMTFRkZWar9bNmyRSEhIQWut9vtstvtpS0TAAC4GZcOQKNGjdKHH36ozz//XL6+vkpLS5Mk+fv7q2rVqpL+unpz8OBBvfvuu5KkWbNmKSIiQs2aNVNWVpbef/99LV68WIsXL660cQAAANfi0gHojTfekCR17drVqX3+/PkaPny4JCk1NVUpKSmOdVlZWZo0aZIOHjyoqlWrqlmzZlq2bJn69OlTUWUDAAAXZzPGmMouwtVkZGTI399f6enp8vPzq+xyLgu/8g13lfx838ouAYCbKcn7t1tMggYAAChLBCAAAGA5BCAAAGA5Lj0JGgAKU9w5blafT1RWcwGtfhyl4h1LjpN74AoQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHM/KLgD5i5iyrLJLACpVWf4NFGdfyc/3dbvnKu7zAciLK0AAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMBy3CIAvf7664qMjJSPj49at26t77//vtD+a9asUevWreXj46P69etrzpw5FVQpAABwBy4fgD766CONHz9eTzzxhLZs2aJOnTopNjZWKSkp+fZPSkpSnz591KlTJ23ZskWPP/64xo4dq8WLF1dw5QAAwFW5fACaOXOmRowYofvvv19NmjTRrFmzFB4erjfeeCPf/nPmzFHdunU1a9YsNWnSRPfff7/uu+8+vfTSSxVcOQAAcFUuHYCysrK0adMmxcTEOLXHxMRo/fr1+W6zYcOGPP179eqln376SefPny+3WgEAgPvwrOwCCnPkyBFlZ2crKCjIqT0oKEhpaWn5bpOWlpZv/wsXLujIkSMKCQnJs01mZqYyMzMdy+np6ZKkjIyMyx1CqeVknqm05wasqKz+3ovzt1uW/22pyP9WVOZ/E11FRZ9flEzusTfGFNnXpQNQLpvN5rRsjMnTVlT//NpzxcfHKy4uLk97eHh4SUsF4Kb8Z12Zz1WW3LXuisZxqnwnT56Uv79/oX1cOgBdffXVqlKlSp6rPYcPH85zlSdXcHBwvv09PT1Vq1atfLeZOnWqJk6c6FjOycnRsWPHVKtWrUKDVmlkZGQoPDxcBw4ckJ+fX5nu2xVc6eOTrvwxMj73d6WPkfG5v/IaozFGJ0+eVGhoaJF9XToAeXt7q3Xr1lq5cqVuvfVWR/vKlSvVv3//fLfp0KGDvvzyS6e2FStWqE2bNvLy8sp3G7vdLrvd7tRWo0aNyyu+CH5+flfsC1u68scnXfljZHzu70ofI+Nzf+UxxqKu/ORy6UnQkjRx4kT9v//3/zRv3jzt3LlTEyZMUEpKih566CFJf129GTp0qKP/Qw89pP3792vixInauXOn5s2bp7lz52rSpEmVNQQAAOBiXPoKkCQNHjxYR48e1dNPP63U1FQ1b95cy5cvV7169SRJqampTt8JFBkZqeXLl2vChAmaPXu2QkND9corr2jgwIGVNQQAAOBiXD4ASdIjjzyiRx55JN91CxYsyNPWpUsXbd68uZyrKh273a5p06blueV2pbjSxydd+WNkfO7vSh8j43N/rjBGmynOZ8UAAACuIC4/BwgAAKCsEYAAAIDlEIAAAIDlEIAAAIDlEIAq0Ouvv67IyEj5+PiodevW+v777yu7pFKJj49X27Zt5evrq8DAQA0YMEC7du1y6jN8+HDZbDanR/v27Sup4pKbPn16nvqDg4Md640xmj59ukJDQ1W1alV17dpVO3bsqMSKSyYiIiLP+Gw2m0aNGiXJPc/f2rVrdfPNNys0NFQ2m01Lly51Wl+cc5aZmakxY8bo6quvVvXq1XXLLbfo999/r8BRFKyw8Z0/f16PPfaYWrRooerVqys0NFRDhw7VH3/84bSPrl275jmvd955ZwWPJH9Fnb/ivCZd+fxJRY8xv79Jm82mF1980dHHVc9hcd4XXO1vkABUQT766CONHz9eTzzxhLZs2aJOnTopNjbW6TuM3MWaNWs0atQo/fDDD1q5cqUuXLigmJgYnT592qlf7969lZqa6ngsX768kiounWbNmjnVv23bNse6F154QTNnztRrr72mjRs3Kjg4WD179tTJkycrseLi27hxo9PYVq5cKUm64447HH3c7fydPn1aUVFReu211/JdX5xzNn78eC1ZskSLFi3SunXrdOrUKfXr10/Z2dkVNYwCFTa+M2fOaPPmzXrqqae0efNmffbZZ9q9e7duueWWPH0feOABp/P65ptvVkT5RSrq/ElFvyZd+fxJRY/x4rGlpqZq3rx5stlseb7HzhXPYXHeF1zub9CgQtxwww3moYcecmq79tprzZQpUyqporJz+PBhI8msWbPG0TZs2DDTv3//yivqMk2bNs1ERUXluy4nJ8cEBweb559/3tF27tw54+/vb+bMmVNBFZatcePGmQYNGpicnBxjjPufP0lmyZIljuXinLMTJ04YLy8vs2jRIkefgwcPGg8PD/PNN99UWO3Fcen48vPjjz8aSWb//v2Oti5duphx48aVb3FlIL/xFfWadKfzZ0zxzmH//v1Nt27dnNrc5Rxe+r7gin+DXAGqAFlZWdq0aZNiYmKc2mNiYrR+/fpKqqrspKenS5ICAgKc2hMTExUYGKhGjRrpgQce0OHDhyujvFLbs2ePQkNDFRkZqTvvvFP79u2TJCUlJSktLc3pfNrtdnXp0sUtz2dWVpbef/993XfffU4//uvu5+9ixTlnmzZt0vnz5536hIaGqnnz5m55XtPT02Wz2fL8ruEHH3ygq6++Ws2aNdOkSZPc5qqlVPhr8ko7f4cOHdKyZcs0YsSIPOvc4Rxe+r7gin+DbvFN0O7uyJEjys7OzvML9kFBQXl+ud7dGGM0ceJE3XTTTWrevLmjPTY2VnfccYfq1aunpKQkPfXUU+rWrZs2bdrkFt9u2q5dO7377rtq1KiRDh06pGeffVYdO3bUjh07HOcsv/O5f//+yij3sixdulQnTpzQ8OHDHW3ufv4uVZxzlpaWJm9vb9WsWTNPH3f7Oz137pymTJmiIUOGOP3Q5F133aXIyEgFBwdr+/btmjp1qn7++WfHLVBXVtRr8ko6f5L0zjvvyNfXV7fddptTuzucw/zeF1zxb5AAVIEu/r9r6a8XyaVt7mb06NH65ZdftG7dOqf2wYMHO/7dvHlztWnTRvXq1dOyZcvy/EG7otjYWMe/W7RooQ4dOqhBgwZ65513HBMvr5TzOXfuXMXGxio0NNTR5u7nryClOWfudl7Pnz+vO++8Uzk5OXr99ded1j3wwAOOfzdv3lwNGzZUmzZttHnzZl1//fUVXWqJlPY16W7nL9e8efN01113ycfHx6ndHc5hQe8Lkmv9DXILrAJcffXVqlKlSp4Ee/jw4Txp2J2MGTNGX3zxhRISEhQWFlZo35CQENWrV0979uypoOrKVvXq1dWiRQvt2bPH8WmwK+F87t+/X6tWrdL9999faD93P3/FOWfBwcHKysrS8ePHC+zj6s6fP69BgwYpKSlJK1eudLr6k5/rr79eXl5ebnleL31NXgnnL9f333+vXbt2Ffl3KbneOSzofcEV/wYJQBXA29tbrVu3znOJcuXKlerYsWMlVVV6xhiNHj1an332mVavXq3IyMgitzl69KgOHDigkJCQCqiw7GVmZmrnzp0KCQlxXH6++HxmZWVpzZo1bnc+58+fr8DAQPXt27fQfu5+/opzzlq3bi0vLy+nPqmpqdq+fbtbnNfc8LNnzx6tWrVKtWrVKnKbHTt26Pz58255Xi99Tbr7+bvY3Llz1bp1a0VFRRXZ11XOYVHvCy75N1jm06qRr0WLFhkvLy8zd+5c8+uvv5rx48eb6tWrm+Tk5MourcQefvhh4+/vbxITE01qaqrjcebMGWOMMSdPnjSPPvqoWb9+vUlKSjIJCQmmQ4cOpk6dOiYjI6OSqy+eRx991CQmJpp9+/aZH374wfTr18/4+vo6ztfzzz9v/P39zWeffWa2bdtm/va3v5mQkBC3GZ8xxmRnZ5u6deuaxx57zKndXc/fyZMnzZYtW8yWLVuMJDNz5kyzZcsWx6eginPOHnroIRMWFmZWrVplNm/ebLp162aioqLMhQsXKmtYDoWN7/z58+aWW24xYWFhZuvWrU5/l5mZmcYYY/bu3Wvi4uLMxo0bTVJSklm2bJm59tprTatWrVx+fMV9Tbry+TOm6NeoMcakp6ebatWqmTfeeCPP9q58Dot6XzDG9f4GCUAVaPbs2aZevXrG29vbXH/99U4fG3cnkvJ9zJ8/3xhjzJkzZ0xMTIypXbu28fLyMnXr1jXDhg0zKSkplVt4CQwePNiEhIQYLy8vExoaam677TazY8cOx/qcnBwzbdo0ExwcbOx2u+ncubPZtm1bJVZcct9++62RZHbt2uXU7q7nLyEhId/X5bBhw4wxxTtnZ8+eNaNHjzYBAQGmatWqpl+/fi4z7sLGl5SUVODfZUJCgjHGmJSUFNO5c2cTEBBgvL29TYMGDczYsWPN0aNHK3dg/7/Cxlfc16Qrnz9jin6NGmPMm2++aapWrWpOnDiRZ3tXPodFvS8Y43p/g7b/v3AAAADLYA4QAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQgHJhs9m0dOnSMt9vcnKybDabtm7dWub7dnVz585VTExMZZehSZMmaezYsZVdBnBZCECAmxg+fLhsNptsNpu8vLwUFBSknj17at68ecrJyans8vJITU1VbGysJGuHlrKSmZmpf/zjH3rqqaccbdOnT5fNZlPv3r3z9H/hhRdks9nUtWvXYu3/0KFD8vLy0vvvv5/v+pEjR+q6666TJE2ePFnz589XUlJSyQcCuAgCEOBGevfurdTUVCUnJ+vrr79WdHS0xo0bp379+unChQuVXZ6T4OBg2e32yi6jQp0/f77c9r148WJdddVV6tSpk1N7SEiIEhIS9Pvvvzu1z58/X3Xr1i32/oOCgtS3b1/Nnz8/z7qzZ89q0aJFGjFihCQpMDBQMTExmjNnTilGArgGAhDgRux2u4KDg1WnTh1df/31evzxx/X555/r66+/1oIFCxz90tPT9eCDDyowMFB+fn7q1q2bfv75Z8f66dOnq2XLlnrzzTcVHh6uatWq6Y477tCJEyccfXJycvT0008rLCxMdrtdLVu21DfffONYn5WVpdGjRyskJEQ+Pj6KiIhQfHy8Y/3Ft8Byfxm6VatWea5KzJ8/X02aNJGPj4+uvfZavf76605j/vHHH9WqVSv5+PioTZs22rJlS5HHKSsrS5MnT1adOnVUvXp1tWvXTomJiY71CxYsUI0aNfTtt9+qSZMmuuqqqxzh8mKF1ZZ7Vevjjz9W165d5ePjo/fff18XLlzQ2LFjVaNGDdWqVUuPPfaYhg0bpgEDBkiS3n33XdWqVUuZmZlOzzVw4EANHTq0wDEtWrRIt9xyS5723DDyzjvvONrWr1+vI0eOqG/fvnn6FzamESNGKCEhQcnJyU7bfPrppzp37pzuvvtuR9stt9yihQsXFlgv4PLK5RfGAJS5YcOGmf79++e7LioqysTGxhpj/vrBwRtvvNHcfPPNZuPGjWb37t3m0UcfNbVq1XL8aOK0adNM9erVTbdu3cyWLVvMmjVrzDXXXGOGDBni2OfMmTONn5+fWbhwofntt9/M5MmTjZeXl9m9e7cxxpgXX3zRhIeHm7Vr15rk5GTz/fffmw8//NCxvSSzZMkSY4wxP/74o5FkVq1aZVJTUx11vPXWWyYkJMQsXrzY7Nu3zyxevNgEBASYBQsWGGOMOXXqlKldu7YZPHiw2b59u/nyyy9N/fr1jSSzZcuWAo/VkCFDTMeOHc3atWvN3r17zYsvvmjsdruj9vnz5xsvLy/To0cPs3HjRrNp0ybTpEkTp/EXVVvuD5BGREQ4+hw8eNA8++yzJiAgwHz22Wdm586d5qGHHjJ+fn6Oc3fmzBnj7+9vPv74Y8dz/fnnn8bb29usXr26wDHVqFHDLFq0yKlt2rRpJioqynz22WfmmmuucbSPGDHCjBs3zowbN8506dKl2GO6cOGCCQkJMdOmTXN6nq5du5pBgwY5tf36669GkklOTi6wZsCVEYAAN1FYABo8eLBp0qSJMcaY7777zvj5+Zlz58459WnQoIF58803jTF/vXFWqVLFHDhwwLH+66+/Nh4eHiY1NdUYY0xoaKh57rnnnPbRtm1b88gjjxhjjBkzZozp1q2bycnJybemiwNQbli4NLSEh4c7hSZjjHnmmWdMhw4djDF//TJ2QECAOX36tGP9G2+8UWgA2rt3r7HZbObgwYNO7d27dzdTp041xvwVgCSZvXv3OtbPnj3bBAUFFbu23DHNmjXLqU9QUJB58cUXHcsXLlwwdevWdTp3Dz/8sCOwGmPMrFmzTP369Qs8lsePHzeSzNq1a53acwNQVlaWCQwMNGvWrDGnTp0yvr6+5ueff84TgIoakzHGPPbYY6ZevXqOWvbt22dsNpv59ttvnbZLT083kkxiYmK+NQOuzrNyrjsBKEvGGNlsNknSpk2bdOrUKdWqVcupz9mzZ/V///d/juW6desqLCzMsdyhQwfl5ORo165dqlatmv744w/deOONTvu48cYbHbfShg8frp49e6px48bq3bu3+vXrV6JPKP355586cOCARowYoQceeMDRfuHCBfn7+0uSdu7cqaioKFWrVs2pzsJs3rxZxhg1atTIqT0zM9PpmFSrVk0NGjRwLIeEhOjw4cPFri1XmzZtHP9OT0/XoUOHdMMNNzjaqlSpotatWztNVH/ggQfUtm1bHTx4UHXq1NH8+fMdk9zzc/bsWUmSj49Pvuu9vLx09913a/78+dq3b58aNWrkmLCcq7hjGjFihGbMmKHVq1ere/fumjdvnsLCwtSjRw+n/VWtWlWSdObMmXxrAlwdAQi4AuzcudMxzyYnJ0chISFOc15y1ahRo8B95L75XvwmfOkb8sVB6/rrr1dSUpK+/vprrVq1SoMGDVKPHj306aefFqvm3EDw9ttvq127dk7rqlSp4ni+ksrJyVGVKlW0adMmx35yXXXVVY5/e3l5Oa2z2WyO5ytObbmqV6+ep4b8jtvFWrVqpaioKL377rvq1auXtm3bpi+//LLAMdWqVUs2m03Hjx8vsM99992ndu3aafv27brvvvvyrC/umBo2bKhOnTpp/vz5io6O1jvvvKN7771XHh7OU0aPHTsmSapdu3aBNQGujAAEuLnVq1dr27ZtmjBhgqS/gklaWpo8PT0VERFR4HYpKSn6448/FBoaKknasGGDPDw81KhRI/n5+Sk0NFTr1q1T586dHdusX7/e6eqGn5+fBg8erMGDB+v2229X7969dezYMQUEBDg9l7e3tyQpOzvb0RYUFKQ6depo3759uuuuu/KtsWnTpnrvvfd09uxZxxWHH374odDj0apVK2VnZ+vw4cN5PjFVXMWpLT/+/v4KCgrSjz/+6Hju7OxsbdmyRS1btnTqe//99+vll1/WwYMH1aNHD4WHhxe4X29vbzVt2lS//vprgVfZmjVrpmbNmumXX37RkCFDLmtMI0aM0MMPP6z+/fvr999/17333punz/bt2+Xl5aVmzZoVui/AVRGAADeSmZmptLQ0ZWdn69ChQ/rmm28UHx+vfv36OT5B1KNHD3Xo0EEDBgzQjBkz1LhxY/3xxx9avny5BgwY4Lhl4+Pjo2HDhumll15SRkaGxo4dq0GDBik4OFiS9Pe//13Tpk1TgwYN1LJlS82fP19bt27VBx98IEl6+eWXFRISopYtW8rDw0OffPKJgoOD873KFBgYqKpVq+qbb75RWFiYfHx85O/vr+nTp2vs2LHy8/NTbGysMjMz9dNPP+n48eOaOHGihgwZoieeeEIjRozQk08+qeTkZL300kuFHqNGjRrprrvu0tChQ/Wvf/1LrVq10pEjR7R69Wq1aNFCffr0KdaxLqq2gowZM0bx8fG65pprdO211+rVV1/V8ePH81wVuuuuuzRp0iS9/fbbevfdd4usp1evXlq3bp3Gjx9fYJ/Vq1fr/PnzBV7pK+6Y7rjjDo0dO1YjR45U9+7d8w3S33//vTp16uQIpoDbqcwJSACKb9iwYUaSkWQ8PT1N7dq1TY8ePcy8efNMdna2U9+MjAwzZswYExoaary8vEx4eLi56667TEpKijHmf5NnX3/9dRMaGmp8fHzMbbfdZo4dO+bYR3Z2tomLizN16tQxXl5eJioqynz99deO9W+99ZZp2bKlqV69uvHz8zPdu3c3mzdvdqzXRZOgjTHm7bffNuHh4cbDw8NpYu4HH3xgWrZsaby9vU3NmjVN586dzWeffeZYv2HDBhMVFWW8vb1Ny5YtzeLFi4v8FFhWVpb5xz/+YSIiIoyXl5cJDg42t956q/nll1+MMX9Ngvb393faZsmSJebS/yQWVltBE7vPnz9vRo8ebfz8/EzNmjXNY489Zu644w5z55135qnznnvuMQEBAXkmrOdn586dpmrVqubEiROOttzzWJBLJ0EXNaaLPfjgg0ZSnknTuRo1amQWLlxYZN2Aq7IZU4qb7ADc2vTp07V06VK+mbkC5OTkqEmTJho0aJCeeeYZp3U9e/ZUkyZN9MorrxRrX4MGDVKrVq00derU8ii12JYtW6a///3v+uWXX+TpyY0EuCe+CBEAytD+/fv19ttva/fu3dq2bZsefvhhJSUlOc3LOXbsmBYtWqTVq1dr1KhRxd73iy++6DSRu7KcPn1a8+fPJ/zArfHqBYAy5OHhoQULFmjSpEkyxqh58+ZatWqVmjRp4uhz/fXX6/jx4445WsVVr149jRkzpjzKLpFBgwZVdgnAZeMWGAAAsBxugQEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMshAAEAAMv5/wBehucdvINzegAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(water_energies,bins=np.linspace(0,200,50))\n",
    "plt.title(\"Water deposited energy\")\n",
    "plt.xlabel(\"Deposited energy (MeV)\")\n",
    "plt.ylabel(\"Counts\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "1c0721bf",
   "metadata": {},
   "outputs": [],
   "source": [
    "#ratio of deposited energies\n",
    "ratio = np.average(water_energies/shield_energies)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "9eff1f93",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Average ratio of energy deposited in the two objects:  0.23551\n"
     ]
    }
   ],
   "source": [
    "print(\"Average ratio of energy deposited in the two objects: \", round(ratio,5))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "29cec8ad",
   "metadata": {},
   "source": [
    "### Changing the material and thickness of shield"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dda56a61",
   "metadata": {},
   "source": [
    "### Tungsten shield of 1.5cm thickness and water sphere of 50cm radius"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "3b5fd293",
   "metadata": {},
   "outputs": [],
   "source": [
    "shield2 = pd.read_csv(\"output_nt_Shield2.csv\")\n",
    "water2 = pd.read_csv(\"output_nt_Water2.csv\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "d5248e68",
   "metadata": {},
   "outputs": [],
   "source": [
    "shield_energies2 = shield2[4:].values.astype(\"float\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "f6cb0b92",
   "metadata": {},
   "outputs": [],
   "source": [
    "water_energies2 = water2[4:].values.astype(\"float\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "2d543391",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total events:  200\n"
     ]
    }
   ],
   "source": [
    "print(\"Total events: \",shield_energies2.shape[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "47fe7b68",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAHFCAYAAAAHcXhbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA6TklEQVR4nO3deXQUZdr+8auBpEMgCQTIRkIShSCbgMKwqEgQAigIooKDOmFEBpVFBQZFZobgFkVFHHGf1yAKggo6Kjuy6QuMCCiLKPCSKAoBWRMQgiT37w9/6bFNQhZDumvm+zmnz6Ge56mqux4K+jpV1d0uMzMBAAA4VDVfFwAAAPBbEGYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWaACvjXv/6l6667To0aNZLb7VZkZKQ6deqksWPHeo1LSEhQnz59St3eqlWr5HK5tGrVqnLXUp51hwwZooSEhHLvo1BCQoKGDBlS4fWrWlpamlwul1fb888/rxkzZpyX/blcLqWlpZ2XbQMoGWEGKKcFCxaoc+fOysnJ0ZQpU7R06VI988wzuuyyyzR37twKbfOSSy7RunXrdMkll1Rytf/dbr/9dq1bt86r7XyGGQC+UcPXBQBOM2XKFCUmJmrJkiWqUePf/4RuuukmTZkypULbDA0NVceOHSurRPx/sbGxio2N9XUZfu2nn36Sy+XyOpcBp+HKDFBOhw8fVv369Yv9z79ateL/SS1evFiXXHKJatasqYsuukivvvqqV39Jt4o+++wzXXvttQoPD1dQUJDatm2rt956q0x1zpgxQ02bNpXb7VazZs00c+bMsh2gfn6DGz9+vKKiohQcHKzLL79cn376abFjs7OzNXz4cMXGxiowMFCJiYmaPHmyzp496xmTlZUll8ulKVOm6JFHHlGjRo0UFBSkdu3a6aOPPiqyzU8++URXXXWVQkJCFBwcrM6dO2vBggVeY3788UeNGzdOiYmJCgoKUnh4uNq1a6c333zTM+bXt5kSEhK0fft2rV69Wi6XSy6Xy+u2W05OjmebgYGBatiwoe655x6dPHnSa985OTkaNmyY6tWrp9q1a6tXr17auXNnmee3rPtxuVwaOXKkXn/9dTVr1kzBwcFq3bq1PvzwwyLb3LVrlwYPHqyIiAjP3/lzzz3nNabwPHv99dc1duxYNWzYUG63W7t375YkvfLKK0pKSpLb7Vbz5s01e/Zsr1uTZqYmTZqoZ8+eRfZ/4sQJhYWFacSIEWWeB6DSGIByuf32202SjRo1ytavX29nzpwpcWx8fLzFxsZa8+bNbebMmbZkyRK78cYbTZKtXr3aM27lypUmyVauXOlpW7FihQUGBtoVV1xhc+fOtcWLF9uQIUNMkmVkZJxz3YyMDJNk/fr1sw8++MDeeOMNa9y4scXFxVl8fHypx5iammoul8v+/Oc/29KlS23q1KnWsGFDCw0NtdTUVM+4/fv3e7b50ksv2fLly+2hhx4yt9ttQ4YM8YzLzMw0SRYXF2eXX365zZs3z95++21r3769BQQE2Nq1az1jV61aZQEBAXbppZfa3Llz7b333rOUlBRzuVw2Z84cz7jhw4dbcHCwTZ061VauXGkffvihPfbYY/bss896xkyaNMl++d/cpk2b7IILLrC2bdvaunXrbN26dbZp0yYzMzt58qS1adPG6tevb1OnTrXly5fbM888Y2FhYdatWzcrKCgwM7OCggJLTk42t9ttjzzyiC1dutQmTZpkF1xwgUmySZMmnXNuy7ofMzNJlpCQYL/73e/srbfesoULF1rXrl2tRo0a9n//93+ecdu3b7ewsDBr1aqVzZw505YuXWpjx461atWqWVpammdc4bnSsGFDu+GGG+z999+3Dz/80A4fPmwvvfSSSbLrr7/ePvzwQ5s1a5YlJSVZfHy81znzzDPPmMvlsp07d3od13PPPWeSbPv27ec8fuB8IMwA5XTo0CG7/PLLTZJJsoCAAOvcubOlp6dbbm6u19j4+HgLCgqyb775xtN26tQpCw8Pt+HDh3vaigskF110kbVt29Z++uknr2326dPHoqOjLT8/v9h18/PzLSYmxi655BKvN8asrCwLCAgoNczs2LHDJNm9997r1T5r1iyT5BVmhg8fbrVr1/Y6PjOzJ5980uuNrTDMxMTE2KlTpzzjcnJyLDw83Lp37+5p69ixo0VERHjN5dmzZ61ly5YWGxvrOaaWLVta//79z3ksvw4zZmYtWrSwK6+8ssjY9PR0q1atmm3YsMGr/Z133jFJtnDhQjMzW7RokUmyZ555xmvcI488UqYwU9b9mP0cZiIjIy0nJ8fTlp2dbdWqVbP09HRPW8+ePS02NtaOHz/utc2RI0daUFCQHTlyxMz+fa506dLFa1x+fr5FRUVZhw4dvNq/+eabIudMTk6OhYSE2N133+01tnnz5pacnHzOYwfOF24zAeVUr149ffzxx9qwYYMee+wx9evXTzt37tSECRPUqlUrHTp0yGt8mzZt1KhRI89yUFCQkpKS9M0335S4j927d+urr77SzTffLEk6e/as53X11Vdr//79+vrrr4td9+uvv9a+ffs0ePBgr1ss8fHx6ty5c6nHt3LlSkny7LvQwIEDi9xa+/DDD5WcnKyYmBivGnv37i1JWr16tdf4AQMGKCgoyLMcEhKivn37as2aNcrPz9fJkyf1r3/9SzfccINq167tGVe9enXdeuut+u677zzH/bvf/U6LFi3S/fffr1WrVunUqVOlHtu5fPjhh2rZsqXatGnjdSw9e/b0ugVY0vwMHjy4UvdTKDk5WSEhIZ7lyMhIRUREeM6f06dP66OPPtJ1112n4ODgIufK6dOntX79eq9tXn/99V7LX3/9tbKzszVw4ECv9kaNGumyyy7zagsJCdEf//hHzZgxw3NbbMWKFfryyy81cuTIMs0BUNkIM0AFtWvXTvfdd5/efvtt7du3T/fee6+ysrKKPARcr169Iuu63e5zvvkeOHBAkjRu3DgFBAR4ve666y5JKhKaCh0+fFiSFBUVVaSvuLayrl+jRo0ix3LgwAF98MEHRWps0aJFsTWWVNOZM2d04sQJHT16VGam6OjoIuNiYmK86vv73/+u++67T++9956Sk5MVHh6u/v37a9euXaUeY3EOHDigLVu2FDmWkJAQmZnnWA4fPlzsXJRlbsuzn0KlnT+HDx/W2bNn9eyzzxbZ5tVXXy2p6N/Dr+e3cE4jIyOL7Ku4tlGjRik3N1ezZs2SJE2fPl2xsbHq169fmeYAqGw8vg5UgoCAAE2aNElPP/20tm3b9pu3V79+fUnShAkTNGDAgGLHNG3atNj2wje/7OzsIn3FtZ1r/YYNG3raz54963nT+2WdF198sR555JFit1UYQM61/+zsbAUGBqp27dqqUaOGqlWrpv379xcZt2/fPs8+JalWrVqaPHmyJk+erAMHDniu0vTt21dfffVVqcf5a/Xr11fNmjWLPJz9y37p5/kpnItfBo2yzG159lNWdevW9Vy5Kunh28TERK/lX3/3TuFxFIboXyruuBo3bqzevXvrueeeU+/evfX+++9r8uTJql69erlqByoLYQYop/379xd75WDHjh2Sir6BV0TTpk3VpEkTffHFF3r00UfLvW50dLTefPNNjRkzxvPG9c0332jt2rWl1te1a1dJ0qxZs3TppZd62t966y2vTyhJUp8+fbRw4UJdeOGFqlu3bqm1zZ8/X0888YTnVlNubq4++OADXXHFFapevbpq1aqlDh06aP78+XryySdVs2ZNSVJBQYHeeOMNxcbGKikpqch2IyMjNWTIEH3xxReaNm2afvzxRwUHBxdbQ0lXxfr06aNHH31U9erVK/Lm/0vJycmaMmWKZs2apdGjR3vaZ8+eXerxl2c/ZRUcHKzk5GRt3rxZF198sQIDA8u9jaZNmyoqKkpvvfWWxowZ42n/9ttvSzxn7r77bqWkpCg1NVXVq1fXsGHDftNxAL8FYQYop549eyo2NlZ9+/bVRRddpIKCAn3++ed66qmnVLt2bd19992Vsp+XXnpJvXv3Vs+ePTVkyBA1bNhQR44c0Y4dO7Rp0ya9/fbbxa5XrVo1PfTQQ7r99tt13XXXadiwYTp27JjS0tLKdCukWbNmuuWWWzRt2jQFBASoe/fu2rZtm5588kmFhoZ6jX3wwQe1bNkyde7cWaNHj1bTpk11+vRpZWVlaeHChXrxxRe9vuelevXq6tGjh8aMGaOCggI9/vjjysnJ0eTJkz1j0tPT1aNHDyUnJ2vcuHEKDAzU888/r23btunNN9/0hLMOHTqoT58+uvjii1W3bl3t2LFDr7/+ujp16lRikJGkVq1aac6cOZo7d64uuOACBQUFqVWrVrrnnns0b948denSRffee68uvvhiFRQU6Ntvv9XSpUs1duxYdejQQSkpKerSpYvGjx+vkydPql27dvrf//1fvf7666XOraQy76c8nnnmGV1++eW64oordOeddyohIUG5ubnavXu3PvjgA61YseKc61erVk2TJ0/W8OHDdcMNN+i2227TsWPHNHnyZEVHRxf7lQM9evRQ8+bNtXLlSt1yyy2KiIgoV81ApfLxA8iA48ydO9cGDx5sTZo0sdq1a1tAQIA1atTIbr31Vvvyyy+9xsbHx9s111xTZBtXXnml1ydqivs0k5nZF198YQMHDrSIiAgLCAiwqKgo69atm7344oulrvuPf/zDmjRpYoGBgZaUlGSvvvqqpaamlumj2Xl5eTZ27FiLiIiwoKAg69ixo61bt87i4+O9Ps1kZvbDDz/Y6NGjLTEx0QICAiw8PNwuvfRSmzhxop04ccLM/v1ppscff9wmT55ssbGxFhgYaG3btrUlS5YU2f/HH39s3bp1s1q1alnNmjWtY8eO9sEHH3iNuf/++61du3ZWt25dc7vddsEFF9i9995rhw4d8owp7tNMWVlZlpKSYiEhISbJaz5OnDhhf/nLX6xp06YWGBjo+bjzvffea9nZ2Z5xx44ds9tuu83q1KljwcHB1qNHD/vqq6/K9Gmm8uxHko0YMaLI+sX9PWRmZtptt91mDRs2tICAAGvQoIF17tzZHn74Yc+YwnPl7bffLraul19+2Ro3bux1zvTr18/atm1b7Pi0tDSTZOvXry/1mIHzyWVm5rMkBeC/QlZWlhITE/XEE09o3Lhxvi4HZXTs2DElJSWpf//+evnll4v0t2vXTi6XSxs2bPBBdcC/cZsJAKDs7Gw98sgjSk5OVr169fTNN9/o6aefVm5urtet05ycHG3btk0ffvihNm7cqHfffdeHVQM/I8wAAOR2u5WVlaW77rpLR44cUXBwsDp27KgXX3zR81F7Sdq0aZMn8EyaNEn9+/f3XdHA/8dtJgAA4Gh8aR4AAHA0wgwAAHA0wgwAAHA0nz4AnJ6ervnz5+urr75SzZo11blzZz3++ONeX9M+ZMgQvfbaa17rdejQocgPp5WkoKBA+/btU0hISJGv8AYAAP7JzJSbm6uYmJhiv7jxl3waZlavXq0RI0aoffv2Onv2rCZOnKiUlBR9+eWXqlWrlmdcr169lJGR4Vkuz9d179u3T3FxcZVaNwAAqBp79+71+ibx4vg0zCxevNhrOSMjQxEREdq4caO6dOniaXe73WX+RdpfCwkJkfTzZPz6q9gBAIB/ysnJUVxcnOd9/Fz86ntmjh8/LkkKDw/3al+1apUiIiJUp04dXXnllXrkkUdK/B2QvLw85eXleZZzc3MlSaGhoYQZAAAcpiyPiPjN98yYmfr166ejR4/q448/9rTPnTtXtWvXVnx8vDIzM/XXv/5VZ8+e1caNG+V2u4tsJy0tzetH6wodP36cMAMAgEPk5OQoLCysTO/ffhNmRowYoQULFuiTTz45572x/fv3Kz4+XnPmzNGAAQOK9P/6ykzhZSrCDAAAzlGeMOMXt5lGjRql999/X2vWrCn1IZ/o6GjFx8dr165dxfa73e5ir9gAAID/TD4NM2amUaNG6d1339WqVauUmJhY6jqHDx/W3r17FR0dXQUVAgAAf+fTL80bMWKE3njjDc2ePVshISHKzs5Wdna2Tp06JUk6ceKExo0bp3Xr1ikrK0urVq1S3759Vb9+fV133XW+LB0AAPgJnz4zU9ITyhkZGRoyZIhOnTql/v37a/PmzTp27Jiio6OVnJyshx56qMzfHVOee24AAMA/OOaZmdJyVM2aNbVkyZIqqgYAADgRv80EAAAcjTADAAAcjTADAAAcjTADAAAcjTADAAAcjTADAAAcjTADAAAcjTADAAAczS9+aBIAAKdKuH9BqWOyHrumCir578WVGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4Gg1fF0AAAD/6RLuX1Bp28p67JpK29Z/Cq7MAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAAR/NpmElPT1f79u0VEhKiiIgI9e/fX19//bXXGDNTWlqaYmJiVLNmTXXt2lXbt2/3UcUAAMDf+DTMrF69WiNGjND69eu1bNkynT17VikpKTp58qRnzJQpUzR16lRNnz5dGzZsUFRUlHr06KHc3FwfVg4AAPxFDV/ufPHixV7LGRkZioiI0MaNG9WlSxeZmaZNm6aJEydqwIABkqTXXntNkZGRmj17toYPH+6LsgEAgB/xq2dmjh8/LkkKDw+XJGVmZio7O1spKSmeMW63W1deeaXWrl1b7Dby8vKUk5Pj9QIAAP+5/CbMmJnGjBmjyy+/XC1btpQkZWdnS5IiIyO9xkZGRnr6fi09PV1hYWGeV1xc3PktHAAA+JTfhJmRI0dqy5YtevPNN4v0uVwur2UzK9JWaMKECTp+/LjntXfv3vNSLwAA8A8+fWam0KhRo/T+++9rzZo1io2N9bRHRUVJ+vkKTXR0tKf94MGDRa7WFHK73XK73ee3YAAA4Dd8emXGzDRy5EjNnz9fK1asUGJiold/YmKioqKitGzZMk/bmTNntHr1anXu3LmqywUAAH7Ip1dmRowYodmzZ+uf//ynQkJCPM/BhIWFqWbNmnK5XLrnnnv06KOPqkmTJmrSpIkeffRRBQcHa/Dgwb4sHQAA+AmfhpkXXnhBktS1a1ev9oyMDA0ZMkSSNH78eJ06dUp33XWXjh49qg4dOmjp0qUKCQmp4moBAIA/8mmYMbNSx7hcLqWlpSktLe38FwQAABzHbz7NBAAAUBGEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4Gg1fF0AAAD+KuH+Bb4uAWXAlRkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBohBkAAOBoPg0za9asUd++fRUTEyOXy6X33nvPq3/IkCFyuVxer44dO/qmWAAA4Jd8GmZOnjyp1q1ba/r06SWO6dWrl/bv3+95LVy4sAorBAAA/q6GL3feu3dv9e7d+5xj3G63oqKiqqgiAADgNH7/zMyqVasUERGhpKQkDRs2TAcPHvR1SQAAwI/49MpMaXr37q0bb7xR8fHxyszM1F//+ld169ZNGzdulNvtLnadvLw85eXleZZzcnKqqlwAAOADfh1mBg0a5Plzy5Yt1a5dO8XHx2vBggUaMGBAseukp6dr8uTJVVUiAMChEu5f4OsSUEn8/jbTL0VHRys+Pl67du0qccyECRN0/Phxz2vv3r1VWCEAAKhqfn1l5tcOHz6svXv3Kjo6usQxbre7xFtQAADgP49Pw8yJEye0e/duz3JmZqY+//xzhYeHKzw8XGlpabr++usVHR2trKwsPfDAA6pfv76uu+46H1YNAAD8iU/DzGeffabk5GTP8pgxYyRJqampeuGFF7R161bNnDlTx44dU3R0tJKTkzV37lyFhIT4qmQAAOBnfBpmunbtKjMrsX/JkiVVWA0AAHAiRz0ADAAA8GuEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GiEGQAA4GgVCjObNm3S1q1bPcv//Oc/1b9/fz3wwAM6c+ZMpRUHAABQmgqFmeHDh2vnzp2SpD179uimm25ScHCw3n77bY0fP75SCwQAADiXCoWZnTt3qk2bNpKkt99+W126dNHs2bM1Y8YMzZs3rzLrAwAAOKcKhRkzU0FBgSRp+fLluvrqqyVJcXFxOnToUOVVBwAAUIoKhZl27drp4Ycf1uuvv67Vq1frmmuukfTzr15HRkZWaoEAAADnUqEw8/TTT2vTpk0aOXKkJk6cqMaNG0uS3nnnHXXu3LlSCwQAADiXCv1qduvWrb0+zVToiSeeUI0aPv0hbgAA8F+mQldmLrjgAh0+fLhI++nTp5WUlPSbiwIAACirCoWZrKws5efnF2nPy8vTd99995uLAgAAKKty3RN6//33PX9esmSJwsLCPMv5+fn66KOPlJiYWHnVAQAAlKJcYaZ///6SJJfLpdTUVK++gIAAJSQk6Kmnnqq04gAAAEpTrjBT+N0yiYmJ2rBhg+rXr39eigIAACirCn30KDMzs7LrAAAAqJAKf476o48+0kcffaSDBw96rtgUevXVV39zYQAAAGVRoTAzefJkPfjgg2rXrp2io6Plcrkquy4AAIAyqVCYefHFFzVjxgzdeuutlV0PAABAuVQozJw5c4afLQAAwAcS7l9Q6pisx66pgkr8R4W+NO/222/X7NmzK7sWAACAcqvQlZnTp0/r5Zdf1vLly3XxxRcrICDAq3/q1KmVUhwAAEBpKhRmtmzZojZt2kiStm3b5tXHw8AAAKAqVSjMrFy5srLrAAAAqJAKPTMDAADgLyp0ZSY5Ofmct5NWrFhR4YIAAADKo0JhpvB5mUI//fSTPv/8c23btq3ID1ACAACcTxUKM08//XSx7WlpaTpx4sRvKggAAKA8KvWZmVtuuYXfZQIAAFWqUsPMunXrFBQUVJmbBAAAOKcK3WYaMGCA17KZaf/+/frss8/017/+tVIKAwAAKIsKhZmwsDCv5WrVqqlp06Z68MEHlZKSUimFAQAAlEWFwkxGRkZl1wEAAFAhFQozhTZu3KgdO3bI5XKpefPmatu2bWXVBQAAUCYVCjMHDx7UTTfdpFWrVqlOnToyMx0/flzJycmaM2eOGjRoUNl1AgAcLOH+BZWynazHrqmyfcE5KvRpplGjRiknJ0fbt2/XkSNHdPToUW3btk05OTkaPXp0ZdcIAABQogpdmVm8eLGWL1+uZs2aedqaN2+u5557jgeAAQBAlarQlZmCggIFBAQUaQ8ICFBBQcFvLgoAAKCsKhRmunXrprvvvlv79u3ztH3//fe69957ddVVV1VacQAAAKWpUJiZPn26cnNzlZCQoAsvvFCNGzdWYmKicnNz9eyzz1Z2jQAAACWq0DMzcXFx2rRpk5YtW6avvvpKZqbmzZure/fulV0fAADAOZXrysyKFSvUvHlz5eTkSJJ69OihUaNGafTo0Wrfvr1atGihjz/++LwUCgAAUJxyhZlp06Zp2LBhCg0NLdIXFham4cOHa+rUqZVWHAAAQGnKFWa++OIL9erVq8T+lJQUbdy48TcXBQAAUFblCjMHDhwo9iPZhWrUqKEffvjhNxcFAABQVuUKMw0bNtTWrVtL7N+yZYuio6N/c1EAAABlVa4wc/XVV+tvf/ubTp8+XaTv1KlTmjRpkvr06VNpxQEAAJSmXB/N/stf/qL58+crKSlJI0eOVNOmTeVyubRjxw4999xzys/P18SJE89XrQAAAEWUK8xERkZq7dq1uvPOOzVhwgSZmSTJ5XKpZ8+eev755xUZGXleCgUAAChOub80Lz4+XgsXLtTRo0e1e/dumZmaNGmiunXrno/6AAAAzqlC3wAsSXXr1lX79u0rsxYAAIByq9BvMwEAAPgLwgwAAHA0wgwAAHA0wgwAAHA0wgwAAHA0wgwAAHA0wgwAAHA0wgwAAHA0n4aZNWvWqG/fvoqJiZHL5dJ7773n1W9mSktLU0xMjGrWrKmuXbtq+/btvikWAAD4JZ+GmZMnT6p169aaPn16sf1TpkzR1KlTNX36dG3YsEFRUVHq0aOHcnNzq7hSAADgryr8cwaVoXfv3urdu3exfWamadOmaeLEiRowYIAk6bXXXlNkZKRmz56t4cOHV2WpAADAT/ntMzOZmZnKzs5WSkqKp83tduvKK6/U2rVrS1wvLy9POTk5Xi8AAPCfy2/DTHZ2tiQpMjLSqz0yMtLTV5z09HSFhYV5XnFxcee1TgAA4Ft+G2YKuVwur2UzK9L2SxMmTNDx48c9r717957vEgEAgA/59JmZc4mKipL08xWa6OhoT/vBgweLXK35JbfbLbfbfd7rAwAA/sFvr8wkJiYqKipKy5Yt87SdOXNGq1evVufOnX1YGQAA8Cc+vTJz4sQJ7d6927OcmZmpzz//XOHh4WrUqJHuuecePfroo2rSpImaNGmiRx99VMHBwRo8eLAPqwYAAP7Ep2Hms88+U3Jysmd5zJgxkqTU1FTNmDFD48eP16lTp3TXXXfp6NGj6tChg5YuXaqQkBBflQwAAPyMT8NM165dZWYl9rtcLqWlpSktLa3qigIAAI7it8/MAAAAlAVhBgAAOBphBgAAOBphBgAAOBphBgAAOBphBgAAOBphBgAAOBphBgAAOJrf/tAkAMAZEu5f8B+5LzgHV2YAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICj1fB1AQAA/5Rw/wJfl4AKKsvfXdZj11RBJVWDKzMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDR/DrMpKWlyeVyeb2ioqJ8XRYAAPAjfv9zBi1atNDy5cs9y9WrV/dhNQAAwN/4fZipUaMGV2MAAECJ/Po2kyTt2rVLMTExSkxM1E033aQ9e/acc3xeXp5ycnK8XgAA4D+XX1+Z6dChg2bOnKmkpCQdOHBADz/8sDp37qzt27erXr16xa6Tnp6uyZMnV3GlAJzmv+1XhYFf+0/6N+DXV2Z69+6t66+/Xq1atVL37t21YMHPE//aa6+VuM6ECRN0/Phxz2vv3r1VVS4AAPABv74y82u1atVSq1attGvXrhLHuN1uud3uKqwKAAD4kl9fmfm1vLw87dixQ9HR0b4uBQAA+Am/DjPjxo3T6tWrlZmZqX/961+64YYblJOTo9TUVF+XBgAA/IRf32b67rvv9Pvf/16HDh1SgwYN1LFjR61fv17x8fG+Lg0AAPgJvw4zc+bM8XUJAADAz/n1bSYAAIDSEGYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICjEWYAAICj+fWX5gEAyi/h/gWljsl67JoqqASoGlyZAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjkaYAQAAjlbD1wUAAKpewv0LfF0CUGm4MgMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNMAMAAByNX80GAD9Rll+yznrsmiqoBHAWrswAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHI8wAAABHq+HrAgCgrBLuX1CmcVmPXVNl+6usfQGoOK7MAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAAR3NEmHn++eeVmJiooKAgXXrppfr44499XRIAAPATfh9m5s6dq3vuuUcTJ07U5s2bdcUVV6h379769ttvfV0aAADwA34fZqZOnaqhQ4fq9ttvV7NmzTRt2jTFxcXphRde8HVpAADAD/h1mDlz5ow2btyolJQUr/aUlBStXbvWR1UBAAB/4te/mn3o0CHl5+crMjLSqz0yMlLZ2dnFrpOXl6e8vDzP8vHjxyVJOTk5569QAFWiIO/HMo0ry7/3sm6rMvZVVmWpqSqPDfDle2fhvs2s1LF+HWYKuVwur2UzK9JWKD09XZMnTy7SHhcXd15qA+B/wqb9Z+7LF/vDfzd/ON9yc3MVFhZ2zjF+HWbq16+v6tWrF7kKc/DgwSJXawpNmDBBY8aM8SwXFBToyJEjqlevXokBqKJycnIUFxenvXv3KjQ0tFK3jX9jnqsG81w1mOeqwTxXjfM5z2am3NxcxcTElDrWr8NMYGCgLr30Ui1btkzXXXedp33ZsmXq169fseu43W653W6vtjp16pzPMhUaGso/lirAPFcN5rlqMM9Vg3muGudrnku7IlPIr8OMJI0ZM0a33nqr2rVrp06dOunll1/Wt99+qzvuuMPXpQEAAD/g92Fm0KBBOnz4sB588EHt379fLVu21MKFCxUfH+/r0gAAgB/w+zAjSXfddZfuuusuX5dRhNvt1qRJk4rc1kLlYp6rBvNcNZjnqsE8Vw1/mWeXleUzTwAAAH7Kr780DwAAoDSEGQAA4GiEGQAA4GiEGQAA4GiEmVKkpaXJ5XJ5vaKiojz9Zqa0tDTFxMSoZs2a6tq1q7Zv3+7Dip1hzZo16tu3r2JiYuRyufTee+959ZdlXvPy8jRq1CjVr19ftWrV0rXXXqvvvvuuCo/C/5U2z0OGDClyfnfs2NFrDPNcuvT0dLVv314hISGKiIhQ//799fXXX3uN4Zz+7coyz5zTv90LL7ygiy++2PNFeJ06ddKiRYs8/f54LhNmyqBFixbav3+/57V161ZP35QpUzR16lRNnz5dGzZsUFRUlHr06KHc3FwfVuz/Tp48qdatW2v69OnF9pdlXu+55x69++67mjNnjj755BOdOHFCffr0UX5+flUdht8rbZ4lqVevXl7n98KFC736mefSrV69WiNGjND69eu1bNkynT17VikpKTp58qRnDOf0b1eWeZY4p3+r2NhYPfbYY/rss8/02WefqVu3burXr58nsPjluWw4p0mTJlnr1q2L7SsoKLCoqCh77LHHPG2nT5+2sLAwe/HFF6uoQueTZO+++65nuSzzeuzYMQsICLA5c+Z4xnz//fdWrVo1W7x4cZXV7iS/nmczs9TUVOvXr1+J6zDPFXPw4EGTZKtXrzYzzunz5dfzbMY5fb7UrVvX/vGPf/jtucyVmTLYtWuXYmJilJiYqJtuukl79uyRJGVmZio7O1spKSmesW63W1deeaXWrl3rq3IdryzzunHjRv30009eY2JiYtSyZUvmvpxWrVqliIgIJSUladiwYTp48KCnj3mumOPHj0uSwsPDJXFOny+/nudCnNOVJz8/X3PmzNHJkyfVqVMnvz2XCTOl6NChg2bOnKklS5bolVdeUXZ2tjp37qzDhw97fs3717/gHRkZWeSXvlF2ZZnX7OxsBQYGqm7duiWOQel69+6tWbNmacWKFXrqqae0YcMGdevWTXl5eZKY54owM40ZM0aXX365WrZsKYlz+nwobp4lzunKsnXrVtWuXVtut1t33HGH3n33XTVv3txvz2VH/JyBL/Xu3dvz51atWqlTp0668MIL9dprr3keKnO5XF7rmFmRNpRfReaVuS+fQYMGef7csmVLtWvXTvHx8VqwYIEGDBhQ4nrMc8lGjhypLVu26JNPPinSxzldeUqaZ87pytG0aVN9/vnnOnbsmObNm6fU1FStXr3a0+9v5zJXZsqpVq1aatWqlXbt2uX5VNOvk+bBgweLpFaUXVnmNSoqSmfOnNHRo0dLHIPyi46OVnx8vHbt2iWJeS6vUaNG6f3339fKlSsVGxvraeecrlwlzXNxOKcrJjAwUI0bN1a7du2Unp6u1q1b65lnnvHbc5kwU055eXnasWOHoqOjlZiYqKioKC1btszTf+bMGa1evVqdO3f2YZXOVpZ5vfTSSxUQEOA1Zv/+/dq2bRtz/xscPnxYe/fuVXR0tCTmuazMTCNHjtT8+fO1YsUKJSYmevVzTleO0ua5OJzTlcPMlJeX57/n8nl5rPg/yNixY23VqlW2Z88eW79+vfXp08dCQkIsKyvLzMwee+wxCwsLs/nz59vWrVvt97//vUVHR1tOTo6PK/dvubm5tnnzZtu8ebNJsqlTp9rmzZvtm2++MbOyzesdd9xhsbGxtnz5ctu0aZN169bNWrdubWfPnvXVYfmdc81zbm6ujR071tauXWuZmZm2cuVK69SpkzVs2JB5Lqc777zTwsLCbNWqVbZ//37P68cff/SM4Zz+7UqbZ87pyjFhwgRbs2aNZWZm2pYtW+yBBx6watWq2dKlS83MP89lwkwpBg0aZNHR0RYQEGAxMTE2YMAA2759u6e/oKDAJk2aZFFRUeZ2u61Lly62detWH1bsDCtXrjRJRV6pqalmVrZ5PXXqlI0cOdLCw8OtZs2a1qdPH/v22299cDT+61zz/OOPP1pKSoo1aNDAAgICrFGjRpaamlpkDpnn0hU3x5IsIyPDM4Zz+rcrbZ45pyvHbbfdZvHx8RYYGGgNGjSwq666yhNkzPzzXHaZmZ2faz4AAADnH8/MAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMAAAARyPMACiVy+XSe++9V+nbzcrKksvl0ueff17p2/Z3//M//6OUlBRfl6Fx48Zp9OjRvi4D+E0IM4APDBkyRC6XSy6XSwEBAYqMjFSPHj306quvqqCgwNflFbF//37PL8j/NweQypKXl6e//e1v+utf/+ppS0tLk8vlUq9evYqMnzJlilwul7p27Vqm7R84cEABAQF64403iu0fPny4Lr74YknS+PHjlZGRoczMzPIfCOAnCDOAj/Tq1Uv79+9XVlaWFi1apOTkZN19993q06ePzp496+vyvERFRcntdvu6jCr1008/nbdtz5s3T7Vr19YVV1zh1R4dHa2VK1fqu+++82rPyMhQo0aNyrz9yMhIXXPNNcrIyCjSd+rUKc2ZM0dDhw6VJEVERCglJUUvvvhiBY4E8A+EGcBH3G63oqKi1LBhQ11yySV64IEH9M9//lOLFi3SjBkzPOOOHz+uP/3pT4qIiFBoaKi6deumL774wtOflpamNm3a6KWXXlJcXJyCg4N144036tixY54xBQUFevDBBxUbGyu32602bdpo8eLFnv4zZ85o5MiRio6OVlBQkBISEpSenu7p/+VtpsJfKm7btm2RqwUZGRlq1qyZgoKCdNFFF+n555/3OuZPP/1Ubdu2VVBQkNq1a6fNmzeXOk9nzpzR+PHj1bBhQ9WqVUsdOnTQqlWrPP0zZsxQnTp1tGTJEjVr1ky1a9f2BMVfOldthVeb3nrrLXXt2lVBQUF64403dPbsWY0ePVp16tRRvXr1dN999yk1NVX9+/eXJM2cOVP16tVTXl6e176uv/56/eEPfyjxmObMmaNrr722SHthsHjttdc8bWvXrtWhQ4d0zTXXFBl/rmMaOnSoVq5cqaysLK913nnnHZ0+fVq33HKLp+3aa6/Vm2++WWK9gN87b7/6BKBEqamp1q9fv2L7Wrdubb179zazn3/Q7bLLLrO+ffvahg0bbOfOnTZ27FirV6+eHT582MzMJk2aZLVq1bJu3brZ5s2bbfXq1da4cWMbPHiwZ5tTp0610NBQe/PNN+2rr76y8ePHW0BAgO3cudPMzJ544gmLi4uzNWvWWFZWln388cc2e/Zsz/qS7N133zUzs08//dQk2fLly23//v2eOl5++WWLjo62efPm2Z49e2zevHkWHh5uM2bMMDOzEydOWIMGDWzQoEG2bds2++CDD+yCCy4wSbZ58+YS52rw4MHWuXNnW7Nmje3evdueeOIJc7vdntozMjIsICDAunfvbhs2bLCNGzdas2bNvI6/tNoyMzNNkiUkJHjGfP/99/bwww9beHi4zZ8/33bs2GF33HGHhYaGev7ufvzxRwsLC7O33nrLs68ffvjBAgMDbcWKFSUeU506dWzOnDlebZMmTbLWrVvb/PnzrXHjxp72oUOH2t1332133323XXnllWU+prNnz1p0dLRNmjTJaz9du3a1gQMHerV9+eWXJsmysrJKrBnwZ4QZwAfOFWYGDRpkzZo1MzOzjz76yEJDQ+306dNeYy688EJ76aWXzOznN8Hq1avb3r17Pf2LFi2yatWq2f79+83MLCYmxh555BGvbbRv397uuusuMzMbNWqUdevWzQoKCoqt6ZdhpvCN/9cBJC4uzisAmZk99NBD1qlTJzMze+mllyw8PNxOnjzp6X/hhRfOGWZ2795tLpfLvv/+e6/2q666yiZMmGBmP4cZSbZ7925P/3PPPWeRkZFlrq3wmKZNm+Y1JjIy0p544gnP8tmzZ61Ro0Zef3d33nmnJ3yamU2bNs0uuOCCEufy6NGjJsnWrFnj1V4YZs6cOWMRERG2evVqO3HihIWEhNgXX3xRJMyUdkxmZvfdd5/Fx8d7atmzZ4+5XC5bsmSJ13rHjx83SbZq1apiawb8XQ3fXA8CUBIzk8vlkiRt3LhRJ06cUL169bzGnDp1Sv/3f//nWW7UqJFiY2M9y506dVJBQYG+/vprBQcHa9++fbrsssu8tnHZZZd5blcNGTJEPXr0UNOmTdWrVy/16dOnXJ+0+eGHH7R3714NHTpUw4YN87SfPXtWYWFhkqQdO3aodevWCg4O9qrzXDZt2iQzU1JSkld7Xl6e15wEBwfrwgsv9CxHR0fr4MGDZa6tULt27Tx/Pn78uA4cOKDf/e53nrbq1avr0ksv9XpIe9iwYWrfvr2+//57NWzYUBkZGZ4HvItz6tQpSVJQUFCx/QEBAbrllluUkZGhPXv2KCkpyfOwbqGyHtPQoUP1+OOPa8WKFbrqqqv06quvKjY2Vt27d/faXs2aNSVJP/74Y7E1Af6OMAP4mR07dnieSykoKFB0dLTXMyKF6tSpU+I2Ct9If/mG+us311+GpksuuUSZmZlatGiRli9froEDB6p79+565513ylRz4Zv7K6+8og4dOnj1Va9e3bO/8iooKFD16tW1ceNGz3YK1a5d2/PngIAArz6Xy+XZX1lqK1SrVq0iNRQ3b7/Utm1btW7dWjNnzlTPnj21detWffDBByUeU7169eRyuXT06NESx9x2223q0KGDtm3bpttuu61If1mPqUmTJrriiiuUkZGh5ORkvfbaa/rjH/+oatW8H5c8cuSIJKlBgwYl1gT4M8IM4EdWrFihrVu36t5775X0c8jIzs5WjRo1lJCQUOJ63377rfbt26eYmBhJ0rp161StWjUlJSUpNDRUMTEx+uSTT9SlSxfPOmvXrvW66hAaGqpBgwZp0KBBuuGGG9SrVy8dOXJE4eHhXvsKDAyUJOXn53vaIiMj1bBhQ+3Zs0c333xzsTU2b95cr7/+uk6dOuW5ErB+/fpzzkfbtm2Vn5+vgwcPFvnkT1mVpbbihIWFKTIyUp9++qln3/n5+dq8ebPatGnjNfb222/X008/re+//17du3dXXFxcidsNDAxU8+bN9eWXX5Z49atFixZq0aKFtmzZosGDB/+mYxo6dKjuvPNO9evXT999953++Mc/Fhmzbds2BQQEqEWLFufcFuCvCDOAj+Tl5Sk7O1v5+fk6cOCAFi9erPT0dPXp08fzSZju3burU6dO6t+/vx5//HE1bdpU+/bt08KFC9W/f3/PbZGgoCClpqbqySefVE5OjkaPHq2BAwcqKipKkvTnP/9ZkyZN0oUXXqg2bdooIyNDn3/+uWbNmiVJevrppxUdHa02bdqoWrVqevvttxUVFVXs1Z+IiAjVrFlTixcvVmxsrIKCghQWFqa0tDSNHj1aoaGh6t27t/Ly8vTZZ5/p6NGjGjNmjAYPHqyJEydq6NCh+stf/qKsrCw9+eST55yjpKQk3XzzzfrDH/6gp556Sm3bttWhQ4e0YsUKtWrVSldffXWZ5rq02koyatQopaenq3Hjxrrooov07LPP6ujRo0Wu1tx8880aN26cXnnlFc2cObPUenr27KlPPvlE99xzT4ljVqxYoZ9++qnEK3BlPaYbb7xRo0eP1vDhw3XVVVcVG4o//vhjXXHFFZ6QCTiOLx/YAf5bpaammiSTZDVq1LAGDRpY9+7d7dVXX7X8/HyvsTk5OTZq1CiLiYmxgIAAi4uLs5tvvtm+/fZbM/v3g6PPP/+8xcTEWFBQkA0YMMCOHDni2UZ+fr5NnjzZGjZsaAEBAda6dWtbtGiRp//ll1+2Nm3aWK1atSw0NNSuuuoq27Rpk6dfv3gA2MzslVdesbi4OKtWrZrXQ6mzZs2yNm3aWGBgoNWtW9e6dOli8+fP9/SvW7fOWrdubYGBgdamTRubN29eqZ9mOnPmjP3tb3+zhIQECwgIsKioKLvuuutsy5YtZvbzA8BhYWFe67z77rv26//ezlVbSQ81//TTTzZy5EgLDQ21unXr2n333Wc33nij3XTTTUXqvPXWWy08PLzIw9rF2bFjh9WsWdOOHTvmaSv8eyzJrx8ALu2YfulPf/qTSSrywHChpKQke/PNN0utG/BXLrMK3MgG4DfS0tL03nvv8Y28VaCgoEDNmjXTwIED9dBDD3n19ejRQ82aNdPf//73Mm1r4MCBatu2rSZMmHA+Si2zBQsW6M9//rO2bNmiGjW4WA9n4kvzAKAE33zzjV555RXt3LlTW7du1Z133qnMzEyv51iOHDmiOXPmaMWKFRoxYkSZt/3EE094PcTsKydPnlRGRgZBBo7G2QsAJahWrZpmzJihcePGyczUsmVLLV++XM2aNfOMueSSS3T06FHPM01lFR8fr1GjRp2Psstl4MCBvi4B+M24zQQAAByN20wAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDRCDMAAMDR/h+PxXXyAWSz6gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(shield_energies2,bins=np.linspace(50,300,50))\n",
    "plt.title(\"Shield deposited energy\")\n",
    "plt.xlabel(\"Deposited energy (MeV)\")\n",
    "plt.ylabel(\"Counts\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "4492c4e8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAHFCAYAAAAaD0bAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8o6BhiAAAACXBIWXMAAA9hAAAPYQGoP6dpAABHpElEQVR4nO3dd3xUVf7/8feEJBPAJBA0DUISEJBmQECKUkIJRFCwgYsKrKioVJFFsCzE1Y2oi6wCon4pVkClWEAFhICuqFQFRMoSCEIC0jLUBJLz+8NfZhnSQ8oM9/V8PObx4J577p3PmZNx3p65M2MzxhgBAABYiFdFFwAAAFDeCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEBAOfjkk09ks9k0f/78XPtiYmJks9n09ddf59pXt25d3XDDDcW6r+nTp2vOnDklLbVURUVFadCgQRVdRpFNnDhRNpvNpa0sH0+bzaaJEyeWybkBFIwABJSDTp06yWazadWqVS7tx44d05YtW1S1atVc+37//Xft2bNHsbGxxbovdwpAnubBBx/U2rVrXdp4PIErk3dFFwBYwdVXX60mTZooKSnJpX316tXy9vbW4MGDcwWgnO3iBqCyYIzRuXPnVLly5YoupUzVqlVLtWrVqugy3Nr58+dls9nk7c3LBzwbK0BAOYmNjdWOHTuUmprqbEtKSlKrVq10yy23aMOGDTp58qTLvkqVKql9+/aSpISEBLVu3VpBQUEKCAjQDTfcoJkzZ+ri3zOOiorStm3btHr1atlsNtlsNkVFRTn3OxwOjRkzRtHR0fL19VXNmjU1atQonT592qVWm82mYcOGacaMGWrYsKHsdrveeeedfMd2/vx5jR07VqGhoapSpYpuvvlm/fTTT3n2TUtL05AhQ1SrVi35+voqOjpaCQkJunDhgrPP3r17ZbPZ9NJLL+mFF15Q7dq15efnp5YtW+qbb77Jdc7vvvtOXbp0kb+/v6pUqaJ27dppyZIlLn3OnDnjHLufn5+CgoLUsmVLzZ0719nn0rfASuvxdDgceuihh1SjRg1dddVV6tGjh3bu3Jnv43mp4s7be++9p4YNG6pKlSqKiYnRF198keucu3btUv/+/RUcHCy73a6GDRtq2rRpLn2SkpJks9n03nvv6YknnlDNmjVlt9u1e/duSdLbb7+t+vXry263q1GjRvrwww81aNAg52NkjFG9evXUvXv3XPd/6tQpBQYGaujQoUV+HIBSZQCUi0WLFhlJ5sMPP3S2NW3a1IwfP96cPHnSeHt7myVLljj3RUdHm1atWjm3Bw0aZGbOnGmWL19uli9fbv7xj3+YypUrm4SEBGefjRs3mjp16pjmzZubtWvXmrVr15qNGzcaY4w5ffq0adasmbn66qvN5MmTzYoVK8y///1vExgYaDp37myys7Od55Fkatasaa6//nrz4YcfmpUrV5qtW7fmO7aBAwcam81m/va3v5lly5aZyZMnm5o1a5qAgAAzcOBAZ7/U1FQTERFhIiMjzZtvvmlWrFhh/vGPfxi73W4GDRrk7JecnGwkmYiICHPzzTebBQsWmI8//ti0atXK+Pj4mO+//97ZNykpyfj4+JgWLVqY+fPnm8WLF5u4uDhjs9nMvHnznP2GDBliqlSpYiZPnmxWrVplvvjiC/Piiy+a119/3dlnwoQJ5uL/LJbG45mdnW1iY2ON3W43L7zwglm2bJmZMGGCqVOnjpFkJkyYkO/jWpJ5i4qKMjfeeKP56KOPzNKlS02nTp2Mt7e3+e9//+vst23bNhMYGGiaNm1q3n33XbNs2TLzxBNPGC8vLzNx4kRnv1WrVjn/Fu666y7z2WefmS+++MIcPXrUvPnmm0aSufPOO80XX3xhPvjgA1O/fn0TGRlpIiMjnef497//bWw2m9m5c6fLuKZNm2YkmW3bthU4fqCsEICAcnLs2DHj5eVlHn74YWOMMUeOHDE2m8189dVXxhhjbrzxRjNmzBhjjDEpKSlGkhk7dmye58rKyjLnz583zz33nKlRo4bLi2Djxo1Nx44dcx2TmJhovLy8zLp161zaP/nkEyPJLF261NkmyQQGBppjx44VOq7t27cbSebxxx93af/ggw+MJJcANGTIEHPVVVeZffv2ufR95ZVXXF4McwJQeHi4OXv2rLOfw+EwQUFBpmvXrs62Nm3amODgYHPy5Eln24ULF0yTJk1MrVq1nI9NkyZNTJ8+fQocy6UByJjLfzy//PJLI8n8+9//dun3wgsvFCkAFXfeQkJCjMPhcLalpaUZLy8vk5iY6Gzr3r27qVWrlklPT3c557Bhw4yfn59z3nMCUIcOHVz6ZWVlmdDQUNO6dWuX9n379hkfHx+XAORwOIy/v78ZOXKkS99GjRqZ2NjYAscOlCXeAgPKSfXq1RUTE+O8Dmj16tWqVKmSbrrpJklSx44dndf95HX9z8qVK9W1a1cFBgaqUqVK8vHx0d///ncdPXpUhw8fLvT+v/jiCzVp0kTNmjXThQsXnLfu3bvLZrPluj6pc+fOql69eqHnzan13nvvdWnv27dvrutEvvjiC8XGxio8PNylhvj4eOdjcrE77rhDfn5+zm1/f3/deuutWrNmjbKysnT69Gn9+OOPuuuuu3TVVVc5+1WqVEn333+/fv/9d+3YsUOSdOONN+rLL7/UuHHjlJSUpLNnzxY6toIU9fHM7/Hp379/qd5PjtjYWPn7+zu3Q0JCFBwcrH379kmSzp07p2+++Ua33367qlSp4nLOW265RefOndMPP/zgcs4777zTZXvHjh1KS0tT3759Xdpr167t/HvO4e/vr7/+9a+aM2eO8y27lStX6tdff9WwYcOK9BgAZYEABJSj2NhY7dy5UwcPHtSqVavUokUL5wt3x44dtWnTJqWnp2vVqlXy9vbWzTffLEn66aefFBcXJ+nP6y7+85//aN26dXr66aclqUgv5ocOHdIvv/wiHx8fl5u/v7+MMTpy5IhL/7CwsCKN6ejRo5Kk0NBQl3Zvb2/VqFEjVw2ff/55rhoaN24sSblquPScOW2ZmZk6deqUjh8/LmNMnrWGh4e71Pfaa6/pySef1OLFixUbG6ugoCD16dNHu3btKtI4L1XUx/Po0aN5PhZ5je1y7ifHpfcjSXa73fk3cvToUV24cEGvv/56rnPecsstknLPw6WPb85jGhISkuu+8mobPny4Tp48qQ8++ECSNHXqVNWqVUu9e/cu0mMAlAUu4wfKUWxsrCZPnqykpCQlJSU5X3AkOcPOmjVrnBdH54SjefPmycfHR1988YXLisjixYuLfN9XX321KleurFmzZuW7/2KXfh9OfnJecNPS0lSzZk1n+4ULF5wvlBffx/XXX68XXnghz3PlhJYcaWlpufqkpaXJ19dXV111lby9veXl5eVyYXmOgwcPOu9TkqpWraqEhAQlJCTo0KFDztWgW2+9Vb/99luRxnrpWIryeNaoUcP5WFwcTvIa2+XcT1FVr17duUKW3wXI0dHRLtuX/i3kjOPQoUO5js1rXNdee63i4+M1bdo0xcfH67PPPlNCQoIqVapUrNqBUlWx78AB1pKenm4qVapkbr/9dmOz2Vyu3zDGmBtuuMHceeedRpJ56qmnnO2jR482V111lcnMzHS2nTlzxtSuXdtIMsnJyS7nuPHGG3Pd9/PPP2+qVKli9uzZU2idkszQoUOLNKZff/21yNcAPfjggyY8PLzQa4sKuwaoS5cuzra2bdua0NBQc+bMGWdbVlaWadq0qcs1QHkZNWqUkWROnz5tjMn7GqDLfTwv9xqg0pi3yMhIl3no2rWriYmJMRkZGQWeL+caoI8//tilvTjXAOVYtmyZkWRiY2ONr6+vOXToUKHjAcoSK0BAOcr5+PrixYvl5eWV63qJjh07asqUKZJcr//p2bOnJk+erP79++vhhx/W0aNH9corr8hut+e6j6ZNm2revHmaP3++6tSpIz8/PzVt2lSjRo3SggUL1KFDBz3++OO6/vrrlZ2drZSUFC1btkxPPPGEWrduXewxNWzYUPfdd5+mTJkiHx8fde3aVVu3btUrr7yigIAAl77PPfecli9frnbt2mnEiBFq0KCBzp07p71792rp0qWaMWOGy/fwVKpUSd26ddPo0aOVnZ2tSZMmyeFwKCEhwdknMTFR3bp1U2xsrMaMGSNfX19Nnz5dW7du1dy5c52rF61bt1avXr10/fXXq3r16tq+fbvee+89tW3bVlWqVMl3fJf7eMbFxalDhw4aO3asTp8+rZYtW+o///mP3nvvvSI9vmUxb//+97918803q3379nr00UcVFRWlkydPavfu3fr888+1cuXKAo/38vJSQkKChgwZorvuuksPPPCATpw4oYSEBIWFhcnLK/fVFd26dVOjRo20atUq3XfffQoODi5WzUCpq+gEBljN2LFjjSTTsmXLXPsWL15sJBlfX1/nqkSOWbNmmQYNGhi73W7q1KljEhMTzcyZM3OtAO3du9fExcUZf39/I8nl/8ZPnTplnnnmGdOgQQPj6+vr/Cj0448/btLS0pz9VIwVIGOMycjIME888YQJDg42fn5+pk2bNmbt2rW5Vh6MMeaPP/4wI0aMMNHR0cbHx8cEBQWZFi1amKefftqcOnXKGPO/FaBJkyaZhIQEU6tWLePr62uaN29uvv7661z3/+2335rOnTubqlWrmsqVK5s2bdqYzz//3KXPuHHjTMuWLU316tWdj+Hjjz9ujhw54uyT1wpQaTyeJ06cMA888ICpVq2aqVKliunWrZv57bffirQCVJz7yW/e8pqH5ORk88ADD5iaNWsaHx8fc80115h27dqZ559/3tknvxWgHG+99Za59tprja+vr6lfv76ZNWuW6d27t2nevHme/SdOnGgkmR9++KHQMQNlzWbMRd+iBgBuYO/evYqOjtbLL7+sMWPGVHQ5KKITJ06ofv366tOnj956661c+1u2bCmbzaZ169ZVQHWAK94CAwAUW1paml544QXFxsaqRo0a2rdvn1599VWdPHlSI0eOdPZzOBzaunWrvvjiC23YsEGLFi2qwKqB/yEAAQCKzW63a+/evXrsscd07NgxValSRW3atNGMGTOcX2sgSRs3bnSGpAkTJqhPnz4VVzRwEd4CAwAAlsMXIQIAAMshAAEAAMshAAEAAMvhIug8ZGdn6+DBg/L39y/yzwEAAICKZYzRyZMnFR4enucXcl6MAJSHgwcPKiIioqLLAAAAJbB//36Xb5XPCwEoD/7+/pL+fAAv/Sp/AADgnhwOhyIiIpyv4wUhAOUh522vgIAAAhAAAB6mKJevcBE0AACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHAIQAACwHO+KLgCeI2rckkL77H2xZzlUAgDA5WEFCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWA4BCAAAWE6FBqDExES1atVK/v7+Cg4OVp8+fbRjxw6XPsYYTZw4UeHh4apcubI6deqkbdu2FXruBQsWqFGjRrLb7WrUqJEWLVpUVsMAAAAepkID0OrVqzV06FD98MMPWr58uS5cuKC4uDidPn3a2eell17S5MmTNXXqVK1bt06hoaHq1q2bTp48me95165dq379+un+++/Xzz//rPvvv199+/bVjz/+WB7DAgAAbs5mjDEVXUSOP/74Q8HBwVq9erU6dOggY4zCw8M1atQoPfnkk5KkjIwMhYSEaNKkSRoyZEie5+nXr58cDoe+/PJLZ1uPHj1UvXp1zZ07t9A6HA6HAgMDlZ6eroCAgNIZ3BWAX4MHALiz4rx+u9U1QOnp6ZKkoKAgSVJycrLS0tIUFxfn7GO329WxY0d9//33+Z5n7dq1LsdIUvfu3fM9JiMjQw6Hw+UGAACuXN4VXUAOY4xGjx6tm2++WU2aNJEkpaWlSZJCQkJc+oaEhGjfvn35nistLS3PY3LOd6nExEQlJCRcTvluqyirNgAAWI3brAANGzZMv/zyS55vUdlsNpdtY0yutss5Zvz48UpPT3fe9u/fX8zqAQCAJ3GLFaDhw4frs88+05o1a1SrVi1ne2hoqKQ/V3TCwsKc7YcPH861wnOx0NDQXKs9BR1jt9tlt9svZwgAAMCDVOgKkDFGw4YN08KFC7Vy5UpFR0e77I+OjlZoaKiWL1/ubMvMzNTq1avVrl27fM/btm1bl2MkadmyZQUeAwAArKNCV4CGDh2qDz/8UJ9++qn8/f2dqzaBgYGqXLmybDabRo0apX/+85+qV6+e6tWrp3/+85+qUqWK+vfv7zzPgAEDVLNmTSUmJkqSRo4cqQ4dOmjSpEnq3bu3Pv30U61YsULfffddhYwTAAC4lwoNQG+88YYkqVOnTi7ts2fP1qBBgyRJY8eO1dmzZ/XYY4/p+PHjat26tZYtWyZ/f39n/5SUFHl5/W8xq127dpo3b56eeeYZPfvss6pbt67mz5+v1q1bl/mYAACA+3Or7wFyF1fS9wCV96fA+B4gAEBF8djvAQIAACgPBCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA5BCAAAGA53hVdAK4sUeOWFNpn74s9y6ESAADyxwoQAACwHAIQAACwHAIQAACwnAoNQGvWrNGtt96q8PBw2Ww2LV682GW/zWbL8/byyy/ne845c+bkecy5c+fKeDQAAMBTVGgAOn36tGJiYjR16tQ896emprrcZs2aJZvNpjvvvLPA8wYEBOQ61s/PryyGAAAAPFCFfgosPj5e8fHx+e4PDQ112f70008VGxurOnXqFHhem82W61gAAIAcHnMN0KFDh7RkyRINHjy40L6nTp1SZGSkatWqpV69emnTpk0F9s/IyJDD4XC5AQCAK5fHBKB33nlH/v7+uuOOOwrsd91112nOnDn67LPPNHfuXPn5+emmm27Srl278j0mMTFRgYGBzltERERplw8AANyIxwSgWbNm6d577y30Wp42bdrovvvuU0xMjNq3b6+PPvpI9evX1+uvv57vMePHj1d6errztn///tIuHwAAuBGP+Cbob7/9Vjt27ND8+fOLfayXl5datWpV4AqQ3W6X3W6/nBIBAIAH8YgVoJkzZ6pFixaKiYkp9rHGGG3evFlhYWFlUBkAAPBEFboCdOrUKe3evdu5nZycrM2bNysoKEi1a9eWJDkcDn388cf617/+lec5BgwYoJo1ayoxMVGSlJCQoDZt2qhevXpyOBx67bXXtHnzZk2bNq3sBwQAADxChQag9evXKzY21rk9evRoSdLAgQM1Z84cSdK8efNkjNFf/vKXPM+RkpIiL6//LWSdOHFCDz/8sNLS0hQYGKjmzZtrzZo1uvHGG8tuIAAAwKPYjDGmootwNw6HQ4GBgUpPT1dAQEBFl3NZivLr7OWNX4MHAJSF4rx+e8Q1QAAAAKWJAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACynQgPQmjVrdOuttyo8PFw2m02LFy922T9o0CDZbDaXW5s2bQo974IFC9SoUSPZ7XY1atRIixYtKqMRAAAAT1ShAej06dOKiYnR1KlT8+3To0cPpaamOm9Lly4t8Jxr165Vv379dP/99+vnn3/W/fffr759++rHH38s7fIBAICH8q7IO4+Pj1d8fHyBfex2u0JDQ4t8zilTpqhbt24aP368JGn8+PFavXq1pkyZorlz515WvQAA4Mrg9tcAJSUlKTg4WPXr19dDDz2kw4cPF9h/7dq1iouLc2nr3r27vv/++3yPycjIkMPhcLkBAIArV4WuABUmPj5ed999tyIjI5WcnKxnn31WnTt31oYNG2S32/M8Ji0tTSEhIS5tISEhSktLy/d+EhMTlZCQUKq1o3xEjVtSaJ+9L/Ysh0oAAJ7ErQNQv379nP9u0qSJWrZsqcjISC1ZskR33HFHvsfZbDaXbWNMrraLjR8/XqNHj3ZuOxwORUREXEblAADAnbl1ALpUWFiYIiMjtWvXrnz7hIaG5lrtOXz4cK5VoYvZ7fZ8V5QAAMCVx+2vAbrY0aNHtX//foWFheXbp23btlq+fLlL27Jly9SuXbuyLg8AAHiICl0BOnXqlHbv3u3cTk5O1ubNmxUUFKSgoCBNnDhRd955p8LCwrR371499dRTuvrqq3X77bc7jxkwYIBq1qypxMRESdLIkSPVoUMHTZo0Sb1799ann36qFStW6Lvvviv38QEAAPdUoQFo/fr1io2NdW7nXIczcOBAvfHGG9qyZYveffddnThxQmFhYYqNjdX8+fPl7+/vPCYlJUVeXv9byGrXrp3mzZunZ555Rs8++6zq1q2r+fPnq3Xr1uU3MAAA4NYqNAB16tRJxph893/99deFniMpKSlX21133aW77rrrckoDAABXMI+6BggAAKA0EIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlEIAAAIDlVOhPYQDlIWrckkL77H2xZzlUAgBwF6wAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAyyEAAQAAy/Gu6AJgPVHjlhTaZ++LPcuhEgCAVbECBAAALIcABAAALIcABAAALKdCA9CaNWt06623Kjw8XDabTYsXL3buO3/+vJ588kk1bdpUVatWVXh4uAYMGKCDBw8WeM45c+bIZrPlup07d66MRwMAADxFhQag06dPKyYmRlOnTs2178yZM9q4caOeffZZbdy4UQsXLtTOnTt12223FXregIAApaamutz8/PzKYggAAMADVeinwOLj4xUfH5/nvsDAQC1fvtyl7fXXX9eNN96olJQU1a5dO9/z2mw2hYaGlmqtAADgyuFR1wClp6fLZrOpWrVqBfY7deqUIiMjVatWLfXq1UubNm0qsH9GRoYcDofLDQAAXLk8JgCdO3dO48aNU//+/RUQEJBvv+uuu05z5szRZ599prlz58rPz0833XSTdu3ale8xiYmJCgwMdN4iIiLKYggAAMBNeEQAOn/+vO655x5lZ2dr+vTpBfZt06aN7rvvPsXExKh9+/b66KOPVL9+fb3++uv5HjN+/Hilp6c7b/v37y/tIQAAADfi9t8Eff78efXt21fJyclauXJlgas/efHy8lKrVq0KXAGy2+2y2+2XWyoAAPAQJVoB2rhxo7Zs2eLc/vTTT9WnTx899dRTyszMLLXicsLPrl27tGLFCtWoUaPY5zDGaPPmzQoLCyu1ugAAgGcrUQAaMmSIdu7cKUnas2eP7rnnHlWpUkUff/yxxo4dW+TznDp1Sps3b9bmzZslScnJydq8ebNSUlJ04cIF3XXXXVq/fr0++OADZWVlKS0tTWlpaS4ha8CAARo/frxzOyEhQV9//bX27NmjzZs3a/Dgwdq8ebMeeeSRkgwVAABcgUr0FtjOnTvVrFkzSdLHH3+sDh066MMPP9R//vMf3XPPPZoyZUqRzrN+/XrFxsY6t0ePHi1JGjhwoCZOnKjPPvtMkpz3lWPVqlXq1KmTJCklJUVeXv/LcSdOnNDDDz+stLQ0BQYGqnnz5lqzZo1uvPHGkgwVAABcgUoUgIwxys7OliStWLFCvXr1kiRFREToyJEjRT5Pp06dZIwp8H4Kk5SU5LL96quv6tVXXy1yDQAAwHpK9BZYy5Yt9fzzz+u9997T6tWr1bNnT0l/voUVEhJSqgUCAACUthIFoFdffVUbN27UsGHD9PTTT+vaa6+VJH3yySdq165dqRYIAABQ2mymKO8zFdG5c+fk7e0tb2+3/3R9gRwOhwIDA5Wenl7sj927m6hxSyq6BI+w98WeFV0CAOAyFef1u0QrQHXq1NHRo0dztZ87d07169cvySkBAADKTYkC0N69e5WVlZWrPSMjQ7///vtlFwUAAFCWivVeVc7H0iXp66+/VmBgoHM7KytL33zzjaKjo0uvOgAAgDJQrADUp08fSZLNZtPAgQNd9vn4+CgqKkr/+te/Sq04AACAslCsAJTz3T/R0dFat26drr766jIpCgAAoCyV6ONaycnJpV0HAABAuSnx59W/+eYbffPNNzp8+LBzZSjHrFmzLrswAACAslKiAJSQkKDnnntOLVu2VFhYmGw2W2nXBQAAUGZKFIBmzJihOXPm6P777y/tegAAAMpcib4HKDMzk5+8AAAAHqtEAejBBx/Uhx9+WNq1AAAAlIsSvQV27tw5vfXWW1qxYoWuv/56+fj4uOyfPHlyqRQHAABQFkoUgH755Rc1a9ZMkrR161aXfVwQDQAA3F2JAtCqVatKuw4AAIByU6JrgAAAADxZiVaAYmNjC3yra+XKlSUuCEUXNW5JRZeAEirK3O19sWc5VAIA1lSiAJRz/U+O8+fPa/Pmzdq6dWuuH0kFAABwNyUKQK+++mqe7RMnTtSpU6cuqyAAAICyVqrXAN133338DhgAAHB7pRqA1q5dKz8/v9I8JQAAQKkr0Vtgd9xxh8u2MUapqalav369nn322VIpDAAAoKyUKAAFBga6bHt5ealBgwZ67rnnFBcXVyqFAQAAlJUSBaDZs2eXdh0AAADlpkQBKMeGDRu0fft22Ww2NWrUSM2bNy+tugAAAMpMiQLQ4cOHdc899ygpKUnVqlWTMUbp6emKjY3VvHnzdM0115R2nQAAAKWmRJ8CGz58uBwOh7Zt26Zjx47p+PHj2rp1qxwOh0aMGFHaNQIAAJSqEq0AffXVV1qxYoUaNmzobGvUqJGmTZvGRdAAAMDtlWgFKDs7Wz4+PrnafXx8lJ2dfdlFAQAAlKUSBaDOnTtr5MiROnjwoLPtwIEDevzxx9WlS5cin2fNmjW69dZbFR4eLpvNpsWLF7vsN8Zo4sSJCg8PV+XKldWpUydt27at0PMuWLBAjRo1kt1uV6NGjbRo0aIi1wQAAK58JQpAU6dO1cmTJxUVFaW6devq2muvVXR0tE6ePKnXX3+9yOc5ffq0YmJiNHXq1Dz3v/TSS5o8ebKmTp2qdevWKTQ0VN26ddPJkyfzPefatWvVr18/3X///fr55591//33q2/fvvrxxx+LPU4AAHBlshljTEkPXr58uX777TcZY9SoUSN17dq15IXYbFq0aJH69Okj6c/Vn/DwcI0aNUpPPvmkJCkjI0MhISGaNGmShgwZkud5+vXrJ4fDoS+//NLZ1qNHD1WvXl1z584tUi0Oh0OBgYFKT09XQEBAicdU1qLGLanoEq4Ye1/sWa73V5S5K++aAMDTFef1u1grQCtXrlSjRo3kcDgkSd26ddPw4cM1YsQItWrVSo0bN9a3335b8sovkpycrLS0NJeLqu12uzp27Kjvv/8+3+PWrl2b60Ls7t27F3hMRkaGHA6Hyw0AAFy5ihWApkyZooceeijPVBUYGKghQ4Zo8uTJpVJYWlqaJCkkJMSlPSQkxLkvv+OKe0xiYqICAwOdt4iIiMuoHAAAuLtiBaCff/5ZPXr0yHd/XFycNmzYcNlFXcxms7lsG2NytV3uMePHj1d6errztn///pIXDAAA3F6xvgfo0KFDeX783Xkyb2/98ccfl12UJIWGhkr6c0UnLCzM2X748OFcKzyXHnfpak9hx9jtdtnt9susGAAAeIpirQDVrFlTW7ZsyXf/L7/84hJWLkd0dLRCQ0O1fPlyZ1tmZqZWr16tdu3a5Xtc27ZtXY6RpGXLlhV4DAAAsJZirQDdcsst+vvf/674+Hj5+fm57Dt79qwmTJigXr16Ffl8p06d0u7du53bycnJ2rx5s4KCglS7dm2NGjVK//znP1WvXj3Vq1dP//znP1WlShX179/fecyAAQNUs2ZNJSYmSpJGjhypDh06aNKkSerdu7c+/fRTrVixQt99911xhgoAAK5gxQpAzzzzjBYuXKj69etr2LBhatCggWw2m7Zv365p06YpKytLTz/9dJHPt379esXGxjq3R48eLUkaOHCg5syZo7Fjx+rs2bN67LHHdPz4cbVu3VrLli2Tv7+/85iUlBR5ef1vIatdu3aaN2+ennnmGT377LOqW7eu5s+fr9atWxdnqAAA4ApW7O8B2rdvnx599FF9/fXXyjnUZrOpe/fumj59uqKiosqiznLF9wBZD98DBACerziv38X+MdTIyEgtXbpUx48f1+7du2WMUb169VS9evUSFwwAAFCeSvRr8JJUvXp1tWrVqjRrAQAAKBcl+i0wAAAAT0YAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAllPiL0IErIafrwCAKwcrQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHLcPgBFRUXJZrPlug0dOjTP/klJSXn2/+2338q5cgAA4K68K7qAwqxbt05ZWVnO7a1bt6pbt266++67Czxux44dCggIcG5fc801ZVYjAADwLG4fgC4NLi+++KLq1q2rjh07FnhccHCwqlWrVoaVAQAAT+X2b4FdLDMzU++//74eeOAB2Wy2Avs2b95cYWFh6tKli1atWlVg34yMDDkcDpcbAAC4cnlUAFq8eLFOnDihQYMG5dsnLCxMb731lhYsWKCFCxeqQYMG6tKli9asWZPvMYmJiQoMDHTeIiIiyqB6AADgLtz+LbCLzZw5U/Hx8QoPD8+3T4MGDdSgQQPndtu2bbV//3698sor6tChQ57HjB8/XqNHj3ZuOxwOQhAAAFcwjwlA+/bt04oVK7Rw4cJiH9umTRu9//77+e632+2y2+2XUx4AAPAgHvMW2OzZsxUcHKyePXsW+9hNmzYpLCysDKoCAACeyCNWgLKzszV79mwNHDhQ3t6uJY8fP14HDhzQu+++K0maMmWKoqKi1LhxY+dF0wsWLNCCBQsqonQAAOCGPCIArVixQikpKXrggQdy7UtNTVVKSopzOzMzU2PGjNGBAwdUuXJlNW7cWEuWLNEtt9xSniUDAAA35hEBKC4uTsaYPPfNmTPHZXvs2LEaO3ZsOVQFAAA8lcdcAwQAAFBaPGIFyIqixi2p6BIAALhisQIEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAshwAEAAAsx7uiCwCQt6hxSwrts/fFnuVQCQBceVgBAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAlkMAAgAAluPWAWjixImy2Wwut9DQ0AKPWb16tVq0aCE/Pz/VqVNHM2bMKKdqAQCAp3D7X4Nv3LixVqxY4dyuVKlSvn2Tk5N1yy236KGHHtL777+v//znP3rsscd0zTXX6M477yyPcgEAgAdw+wDk7e1d6KpPjhkzZqh27dqaMmWKJKlhw4Zav369XnnlFQIQAABwcuu3wCRp165dCg8PV3R0tO655x7t2bMn375r165VXFycS1v37t21fv16nT9/Pt/jMjIy5HA4XG4AAODK5dYrQK1bt9a7776r+vXr69ChQ3r++efVrl07bdu2TTVq1MjVPy0tTSEhIS5tISEhunDhgo4cOaKwsLA87ycxMVEJCQllMgagLEWNW1Jon70v9iyHSgDAs7j1ClB8fLzuvPNONW3aVF27dtWSJX/+x/6dd97J9xibzeaybYzJs/1i48ePV3p6uvO2f//+UqgeAAC4K7deAbpU1apV1bRpU+3atSvP/aGhoUpLS3NpO3z4sLy9vfNcMcpht9tlt9tLtVYAAOC+3HoF6FIZGRnavn17vm9ltW3bVsuXL3dpW7ZsmVq2bCkfH5/yKBEAAHgAtw5AY8aM0erVq5WcnKwff/xRd911lxwOhwYOHCjpz7euBgwY4Oz/yCOPaN++fRo9erS2b9+uWbNmaebMmRozZkxFDQEAALght34L7Pfff9df/vIXHTlyRNdcc43atGmjH374QZGRkZKk1NRUpaSkOPtHR0dr6dKlevzxxzVt2jSFh4frtdde4yPwAADAhVsHoHnz5hW4f86cObnaOnbsqI0bN5ZRRQAA4Erg1m+BAQAAlAUCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBwCEAAAsBy3/ikMoLxEjVtS0SUAAMoRK0AAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByvCu6AOBKEjVuSUWXAAAoAlaAAACA5RCAAACA5RCAAACA5RCAAACA5bh1AEpMTFSrVq3k7++v4OBg9enTRzt27CjwmKSkJNlstly33377rZyqBgAA7s6tA9Dq1as1dOhQ/fDDD1q+fLkuXLiguLg4nT59utBjd+zYodTUVOetXr165VAxAADwBG79MfivvvrKZXv27NkKDg7Whg0b1KFDhwKPDQ4OVrVq1cqwOgAA4KncegXoUunp6ZKkoKCgQvs2b95cYWFh6tKli1atWlVg34yMDDkcDpcbAAC4cnlMADLGaPTo0br55pvVpEmTfPuFhYXprbfe0oIFC7Rw4UI1aNBAXbp00Zo1a/I9JjExUYGBgc5bREREWQwBAAC4CZsxxlR0EUUxdOhQLVmyRN99951q1apVrGNvvfVW2Ww2ffbZZ3nuz8jIUEZGhnPb4XAoIiJC6enpCggIuKy6S4pvFEZp2ftiz4ouAQDKhcPhUGBgYJFevz1iBWj48OH67LPPtGrVqmKHH0lq06aNdu3ale9+u92ugIAAlxsAALhyufVF0MYYDR8+XIsWLVJSUpKio6NLdJ5NmzYpLCyslKsDAACeyq0D0NChQ/Xhhx/q008/lb+/v9LS0iRJgYGBqly5siRp/PjxOnDggN59911J0pQpUxQVFaXGjRsrMzNT77//vhYsWKAFCxZU2DgAAIB7cesA9MYbb0iSOnXq5NI+e/ZsDRo0SJKUmpqqlJQU577MzEyNGTNGBw4cUOXKldW4cWMtWbJEt9xyS3mVDQAA3JzHXARdnopzEVVZ4SJolBYuggZgFVfcRdAAAAClya3fAgNw+YqymljUVaLSOldp1lSePLVuALmxAgQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACyHAAQAACzHu6ILAFDxosYtcctzXamK+hjtfbFnGVcCWBcrQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHIIQAAAwHI8IgBNnz5d0dHR8vPzU4sWLfTtt98W2H/16tVq0aKF/Pz8VKdOHc2YMaOcKgUAAJ7A7QPQ/PnzNWrUKD399NPatGmT2rdvr/j4eKWkpOTZPzk5Wbfccovat2+vTZs26amnntKIESO0YMGCcq4cAAC4K7cPQJMnT9bgwYP14IMPqmHDhpoyZYoiIiL0xhtv5Nl/xowZql27tqZMmaKGDRvqwQcf1AMPPKBXXnmlnCsHAADuyq0DUGZmpjZs2KC4uDiX9ri4OH3//fd5HrN27dpc/bt3767169fr/PnzZVYrAADwHN4VXUBBjhw5oqysLIWEhLi0h4SEKC0tLc9j0tLS8ux/4cIFHTlyRGFhYbmOycjIUEZGhnM7PT1dkuRwOC53CCWWnXGmwu4bcAcV+fzLT1Gel0Wpu6jPb3d8DAB3lvOcMcYU2tetA1AOm83msm2MydVWWP+82nMkJiYqISEhV3tERERxSwVQSgKnVHQFJVOadXvqYwBUtJMnTyowMLDAPm4dgK6++mpVqlQp12rP4cOHc63y5AgNDc2zv7e3t2rUqJHnMePHj9fo0aOd29nZ2Tp27Jhq1KhRYNAqCYfDoYiICO3fv18BAQGlem53cKWPT7ryx8j4PN+VPkbG5/nKaozGGJ08eVLh4eGF9nXrAOTr66sWLVpo+fLluv32253ty5cvV+/evfM8pm3btvr8889d2pYtW6aWLVvKx8cnz2PsdrvsdrtLW7Vq1S6v+EIEBARcsX/Y0pU/PunKHyPj83xX+hgZn+crizEWtvKTw60vgpak0aNH6//+7/80a9Ysbd++XY8//rhSUlL0yCOPSPpz9WbAgAHO/o888oj27dun0aNHa/v27Zo1a5ZmzpypMWPGVNQQAACAm3HrFSBJ6tevn44eParnnntOqampatKkiZYuXarIyEhJUmpqqst3AkVHR2vp0qV6/PHHNW3aNIWHh+u1117TnXfeWVFDAAAAbsbtA5AkPfbYY3rsscfy3DdnzpxcbR07dtTGjRvLuKqSsdvtmjBhQq633K4UV/r4pCt/jIzP813pY2R8ns8dxmgzRfmsGAAAwBXE7a8BAgAAKG0EIAAAYDkEIAAAYDkEIAAAYDkEoHI0ffp0RUdHy8/PTy1atNC3335b0SWVSGJiolq1aiV/f38FBwerT58+2rFjh0ufQYMGyWazudzatGlTQRUX38SJE3PVHxoa6txvjNHEiRMVHh6uypUrq1OnTtq2bVsFVlw8UVFRucZns9k0dOhQSZ45f2vWrNGtt96q8PBw2Ww2LV682GV/UeYsIyNDw4cP19VXX62qVavqtttu0++//16Oo8hfQeM7f/68nnzySTVt2lRVq1ZVeHi4BgwYoIMHD7qco1OnTrnm9Z577innkeStsPkryt+kO8+fVPgY83pO2mw2vfzyy84+7jqHRXldcLfnIAGonMyfP1+jRo3S008/rU2bNql9+/aKj493+Q4jT7F69WoNHTpUP/zwg5YvX64LFy4oLi5Op0+fdunXo0cPpaamOm9Lly6toIpLpnHjxi71b9myxbnvpZde0uTJkzV16lStW7dOoaGh6tatm06ePFmBFRfdunXrXMa2fPlySdLdd9/t7ONp83f69GnFxMRo6tSpee4vypyNGjVKixYt0rx58/Tdd9/p1KlT6tWrl7KyssprGPkqaHxnzpzRxo0b9eyzz2rjxo1auHChdu7cqdtuuy1X34ceeshlXt98883yKL9Qhc2fVPjfpDvPn1T4GC8eW2pqqmbNmiWbzZbre+zccQ6L8rrgds9Bg3Jx4403mkceecSl7brrrjPjxo2roIpKz+HDh40ks3r1amfbwIEDTe/evSuuqMs0YcIEExMTk+e+7OxsExoaal588UVn27lz50xgYKCZMWNGOVVYukaOHGnq1q1rsrOzjTGeP3+SzKJFi5zbRZmzEydOGB8fHzNv3jxnnwMHDhgvLy/z1VdflVvtRXHp+PLy008/GUlm3759zraOHTuakSNHlm1xpSCv8RX2N+lJ82dM0eawd+/epnPnzi5tnjKHl74uuONzkBWgcpCZmakNGzYoLi7OpT0uLk7ff/99BVVVetLT0yVJQUFBLu1JSUkKDg5W/fr19dBDD+nw4cMVUV6J7dq1S+Hh4YqOjtY999yjPXv2SJKSk5OVlpbmMp92u10dO3b0yPnMzMzU+++/rwceeMDlx389ff4uVpQ527Bhg86fP+/SJzw8XE2aNPHIeU1PT5fNZsv1u4YffPCBrr76ajVu3FhjxozxmFVLqeC/yStt/g4dOqQlS5Zo8ODBufZ5whxe+rrgjs9Bj/gmaE935MgRZWVl5foF+5CQkFy/XO9pjDEaPXq0br75ZjVp0sTZHh8fr7vvvluRkZFKTk7Ws88+q86dO2vDhg0e8e2mrVu31rvvvqv69evr0KFDev7559WuXTtt27bNOWd5zee+ffsqotzLsnjxYp04cUKDBg1ytnn6/F2qKHOWlpYmX19fVa9ePVcfT3uenjt3TuPGjVP//v1dfmjy3nvvVXR0tEJDQ7V161aNHz9eP//8s/MtUHdW2N/klTR/kvTOO+/I399fd9xxh0u7J8xhXq8L7vgcJACVo4v/71r684/k0jZPM2zYMP3yyy/67rvvXNr79evn/HeTJk3UsmVLRUZGasmSJbme0O4oPj7e+e+mTZuqbdu2qlu3rt555x3nhZdXynzOnDlT8fHxCg8Pd7Z5+vzlpyRz5mnzev78ed1zzz3Kzs7W9OnTXfY99NBDzn83adJE9erVU8uWLbVx40bdcMMN5V1qsZT0b9LT5i/HrFmzdO+998rPz8+l3RPmML/XBcm9noO8BVYOrr76alWqVClXgj18+HCuNOxJhg8frs8++0yrVq1SrVq1CuwbFhamyMhI7dq1q5yqK11Vq1ZV06ZNtWvXLuenwa6E+dy3b59WrFihBx98sMB+nj5/RZmz0NBQZWZm6vjx4/n2cXfnz59X3759lZycrOXLl7us/uTlhhtukI+Pj0fO66V/k1fC/OX49ttvtWPHjkKfl5L7zWF+rwvu+BwkAJUDX19ftWjRItcS5fLly9WuXbsKqqrkjDEaNmyYFi5cqJUrVyo6OrrQY44ePar9+/crLCysHCosfRkZGdq+fbvCwsKcy88Xz2dmZqZWr17tcfM5e/ZsBQcHq2fPngX28/T5K8qctWjRQj4+Pi59UlNTtXXrVo+Y15zws2vXLq1YsUI1atQo9Jht27bp/PnzHjmvl/5Nevr8XWzmzJlq0aKFYmJiCu3rLnNY2OuCWz4HS/2yauRp3rx5xsfHx8ycOdP8+uuvZtSoUaZq1apm7969FV1asT366KMmMDDQJCUlmdTUVOftzJkzxhhjTp48aZ544gnz/fffm+TkZLNq1SrTtm1bU7NmTeNwOCq4+qJ54oknTFJSktmzZ4/54YcfTK9evYy/v79zvl588UUTGBhoFi5caLZs2WL+8pe/mLCwMI8ZnzHGZGVlmdq1a5snn3zSpd1T5+/kyZNm06ZNZtOmTUaSmTx5stm0aZPzU1BFmbNHHnnE1KpVy6xYscJs3LjRdO7c2cTExJgLFy5U1LCcChrf+fPnzW233WZq1aplNm/e7PK8zMjIMMYYs3v3bpOQkGDWrVtnkpOTzZIlS8x1111nmjdv7vbjK+rfpDvPnzGF/40aY0x6erqpUqWKeeONN3Id785zWNjrgjHu9xwkAJWjadOmmcjISOPr62tuuOEGl4+NexJJed5mz55tjDHmzJkzJi4uzlxzzTXGx8fH1K5d2wwcONCkpKRUbOHF0K9fPxMWFmZ8fHxMeHi4ueOOO8y2bduc+7Ozs82ECRNMaGiosdvtpkOHDmbLli0VWHHxff3110aS2bFjh0u7p87fqlWr8vy7HDhwoDGmaHN29uxZM2zYMBMUFGQqV65sevXq5TbjLmh8ycnJ+T4vV61aZYwxJiUlxXTo0MEEBQUZX19fU7duXTNixAhz9OjRih3Y/1fQ+Ir6N+nO82dM4X+jxhjz5ptvmsqVK5sTJ07kOt6d57Cw1wVj3O85aPv/hQMAAFgG1wABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABAADLIQABKBM2m02LFy8u9fPu3btXNptNmzdvLvVzu7uZM2cqLi6uosvQmDFjNGLEiIouA7gsBCDAQwwaNEg2m002m00+Pj4KCQlRt27dNGvWLGVnZ1d0ebmkpqYqPj5ekrVDS2nJyMjQ3//+dz377LPOtokTJ8pms6lHjx65+r/00kuy2Wzq1KlTkc5/6NAh+fj46P33389z/5AhQ3T99ddLksaOHavZs2crOTm5+AMB3AQBCPAgPXr0UGpqqvbu3asvv/xSsbGxGjlypHr16qULFy5UdHkuQkNDZbfbK7qMcnX+/PkyO/eCBQt01VVXqX379i7tYWFhWrVqlX7//XeX9tmzZ6t27dpFPn9ISIh69uyp2bNn59p39uxZzZs3T4MHD5YkBQcHKy4uTjNmzCjBSAD3QAACPIjdbldoaKhq1qypG264QU899ZQ+/fRTffnll5ozZ46zX3p6uh5++GEFBwcrICBAnTt31s8//+zcP3HiRDVr1kxvvvmmIiIiVKVKFd199906ceKEs092draee+451apVS3a7Xc2aNdNXX33l3J+Zmalhw4YpLCxMfn5+ioqKUmJionP/xW+B5fwydPPmzXOtSsyePVsNGzaUn5+frrvuOk2fPt1lzD/99JOaN28uPz8/tWzZUps2bSr0ccrMzNTYsWNVs2ZNVa1aVa1bt1ZSUpJz/5w5c1StWjV9/fXXatiwoa666ipnuLxYQbXlrGp99NFH6tSpk/z8/PT+++/rwoULGjFihKpVq6YaNWroySef1MCBA9WnTx9J0rvvvqsaNWooIyPD5b7uvPNODRgwIN8xzZs3T7fddluu9pww8s477zjbvv/+ex05ckQ9e/bM1b+gMQ0ePFirVq3S3r17XY755JNPdO7cOd13333Otttuu01z587Nt17A7ZXJL4wBKHUDBw40vXv3znNfTEyMiY+PN8b8+YODN910k7n11lvNunXrzM6dO80TTzxhatSo4fzRxAkTJpiqVauazp07m02bNpnVq1eba6+91vTv3995zsmTJ5uAgAAzd+5c89tvv5mxY8caHx8fs3PnTmOMMS+//LKJiIgwa9asMXv37jXffvut+fDDD53HSzKLFi0yxhjz008/GUlmxYoVJjU11VnHW2+9ZcLCwsyCBQvMnj17zIIFC0xQUJCZM2eOMcaYU6dOmWuuucb069fPbN261Xz++eemTp06RpLZtGlTvo9V//79Tbt27cyaNWvM7t27zcsvv2zsdruz9tmzZxsfHx/TtWtXs27dOrNhwwbTsGFDl/EXVlvOD5BGRUU5+xw4cMA8//zzJigoyCxcuNBs377dPPLIIyYgIMA5d2fOnDGBgYHmo48+ct7XH3/8YXx9fc3KlSvzHVO1atXMvHnzXNomTJhgYmJizMKFC821117rbB88eLAZOXKkGTlypOnYsWORx3ThwgUTFhZmJkyY4HI/nTp1Mn379nVp+/XXX40ks3fv3nxrBtwZAQjwEAUFoH79+pmGDRsaY4z55ptvTEBAgDl37pxLn7p165o333zTGPPnC2elSpXM/v37nfu//PJL4+XlZVJTU40xxoSHh5sXXnjB5RytWrUyjz32mDHGmOHDh5vOnTub7OzsPGu6OADlhIVLQ0tERIRLaDLGmH/84x+mbdu2xpg/fxk7KCjInD592rn/jTfeKDAA7d6929hsNnPgwAGX9i5dupjx48cbY/4MQJLM7t27nfunTZtmQkJCilxbzpimTJni0ickJMS8/PLLzu0LFy6Y2rVru8zdo48+6gysxhgzZcoUU6dOnXwfy+PHjxtJZs2aNS7tOQEoMzPTBAcHm9WrV5tTp04Zf39/8/PPP+cKQIWNyRhjnnzySRMZGemsZc+ePcZms5mvv/7a5bj09HQjySQlJeVZM+DuvCtm3QlAaTLGyGazSZI2bNigU6dOqUaNGi59zp49q//+97/O7dq1a6tWrVrO7bZt2yo7O1s7duxQlSpVdPDgQd10000u57jpppucb6UNGjRI3bp1U4MGDdSjRw/16tWrWJ9Q+uOPP7R//34NHjxYDz30kLP9woULCgwMlCRt375dMTExqlKlikudBdm4caOMMapfv75Le0ZGhstjUqVKFdWtW9e5HRYWpsOHDxe5thwtW7Z0/js9PV2HDh3SjTfe6GyrVKmSWrRo4XKh+kMPPaRWrVrpwIEDqlmzpmbPnu28yD0vZ8+elST5+fnlud/Hx0f33XefZs+erT179qh+/frOC5ZzFHVMgwcP1qRJk7Ry5Up16dJFs2bNUq1atdS1a1eX81WuXFmSdObMmTxrAtwdAQi4Amzfvt15nU12drbCwsJcrnnJUa1atXzPkfPie/GL8KUvyBcHrRtuuEHJycn68ssvtWLFCvXt21ddu3bVJ598UqSacwLB22+/rdatW7vsq1SpkvP+iis7O1uVKlXShg0bnOfJcdVVVzn/7ePj47LPZrM5768oteWoWrVqrhryetwu1rx5c8XExOjdd99V9+7dtWXLFn3++ef5jqlGjRqy2Ww6fvx4vn0eeOABtW7dWlu3btUDDzyQa39Rx1SvXj21b99es2fPVmxsrN555x399a9/lZeX6yWjx44dkyRdc801+dYEuDMCEODhVq5cqS1btujxxx+X9GcwSUtLk7e3t6KiovI9LiUlRQcPHlR4eLgkae3atfLy8lL9+vUVEBCg8PBwfffdd+rQoYPzmO+//95ldSMgIED9+vVTv379dNddd6lHjx46duyYgoKCXO7L19dXkpSVleVsCwkJUc2aNbVnzx7de++9edbYqFEjvffeezp79qxzxeGHH34o8PFo3ry5srKydPjw4VyfmCqqotSWl8DAQIWEhOinn35y3ndWVpY2bdqkZs2aufR98MEH9eqrr+rAgQPq2rWrIiIi8j2vr6+vGjVqpF9//TXfVbbGjRurcePG+uWXX9S/f//LGtPgwYP16KOPqnfv3vr999/117/+NVefrVu3ysfHR40bNy7wXIC7IgABHiQjI0NpaWnKysrSoUOH9NVXXykxMVG9evVyfoKoa9euatu2rfr06aNJkyapQYMGOnjwoJYuXao+ffo437Lx8/PTwIED9corr8jhcGjEiBHq27evQkNDJUl/+9vfNGHCBNWtW1fNmjXT7NmztXnzZn3wwQeSpFdffVVhYWFq1qyZvLy89PHHHys0NDTPVabg4GBVrlxZX331lWrVqiU/Pz8FBgZq4sSJGjFihAICAhQfH6+MjAytX79ex48f1+jRo9W/f389/fTTGjx4sJ555hnt3btXr7zySoGPUf369XXvvfdqwIAB+te//qXmzZvryJEjWrlypZo2bapbbrmlSI91YbXlZ/jw4UpMTNS1116r6667Tq+//rqOHz+ea1Xo3nvv1ZgxY/T222/r3XffLbSe7t2767vvvtOoUaPy7bNy5UqdP38+35W+oo7p7rvv1ogRIzRkyBB16dIlzyD97bffqn379s5gCnicirwACUDRDRw40Egykoy3t7e55pprTNeuXc2sWbNMVlaWS1+Hw2GGDx9uwsPDjY+Pj4mIiDD33nuvSUlJMcb87+LZ6dOnm/DwcOPn52fuuOMOc+zYMec5srKyTEJCgqlZs6bx8fExMTEx5ssvv3Tuf+utt0yzZs1M1apVTUBAgOnSpYvZuHGjc78uugjaGGPefvttExERYby8vFwuzP3ggw9Ms2bNjK+vr6levbrp0KGDWbhwoXP/2rVrTUxMjPH19TXNmjUzCxYsKPRTYJmZmebvf/+7iYqKMj4+PiY0NNTcfvvt5pdffjHG/HkRdGBgoMsxixYtMpf+J7Gg2vK7sPv8+fNm2LBhJiAgwFSvXt08+eST5u677zb33HNPrjrvv/9+ExQUlOuC9bxs377dVK5c2Zw4ccLZljOP+bn0IujCxnSxhx9+2EjKddF0jvr165u5c+cWWjfgrmzGlOBNdgAebeLEiVq8eDHfzFwOsrOz1bBhQ/Xt21f/+Mc/XPZ169ZNDRs21GuvvVakc/Xt21fNmzfX+PHjy6LUIluyZIn+9re/6ZdffpG3N28kwDPxRYgAUIr27dunt99+Wzt37tSWLVv06KOPKjk52eW6nGPHjmnevHlauXKlhg4dWuRzv/zyyy4XcleU06dPa/bs2YQfeDT+egGgFHl5eWnOnDkaM2aMjDFq0qSJVqxYoYYNGzr73HDDDTp+/LjzGq2iioyM1PDhw8ui7GLp27dvRZcAXDbeAgMAAJbDW2AAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMByCEAAAMBy/h//OvxdkDt3VwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(water_energies2,bins=np.linspace(0,200,50))\n",
    "plt.title(\"Water deposited energy\")\n",
    "plt.xlabel(\"Deposited energy (MeV)\")\n",
    "plt.ylabel(\"Counts\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "05ef2286",
   "metadata": {},
   "outputs": [],
   "source": [
    "#ratio of deposited energies\n",
    "ratio2 = np.average(water_energies2/shield_energies2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "c2ccb9ff",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Average ratio of energy deposited in the two objects:  0.14401\n"
     ]
    }
   ],
   "source": [
    "print(\"Average ratio of energy deposited in the two objects: \", round(ratio2,5))"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5b9c0ecb",
   "metadata": {},
   "source": [
    "The average ratio of deposited energy decreases.\n",
    "\n",
    "We can see that the slightly thinner tungsten shield actually does better than the thicker lead shield.\n",
    "\n",
    "This is because tungsten is a better absorber of radiation than lead and requires less thickness to provide similar levels of protection."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "05587332",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "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.9.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}