{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Special cases" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Geant4 analysis tools output" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from physt.compat.geant4 import load_csv" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Histogram1D('Edep in absorber', bins=(100,), total=10000.0, dtype=float64)" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h = load_csv(\"../tests/data/geant-h1.csv\")\n", "h" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\janpi\\Documents\\code\\physt\\src\\physt\\statistics.py:36: RuntimeWarning:\n", "\n", "divide by zero encountered in scalar divide\n", "\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABiE0lEQVR4nO3deVxU5eI/8M8wMMPmgIIwoIAoJm64JlGmliQqWpr3tlxUTNNUqNRyuz/Xule+Wd20NG0T6qbXtNRKSyPc0hCNwl1cUrFkcYORdWDm+f2Bc2BiUXSGGWY+79eL13XOeebMc4409+OzyoQQAkRERETU5DlYugJEREREZBoMdkREREQ2gsGOiIiIyEYw2BERERHZCAY7IiIiIhvBYEdERERkIxjsiIiIiGwEgx0RERGRjWCwIyIiIrIRDHZEZFIXLlyATCZDUlJSo3/2okWLIJPJGv1zDff81ltvNfpn16VNmzYYNmyYpatBRI2MwY6IjCQlJUEmk9X5c+DAAUtXkYiI6uBo6QoQkXV67bXXEBwcXON4SEiIBWpzZ+bNm4c5c+ZYuhpERBbDYEdEtRoyZAh69+5t6Wo0iKOjIxwd7ftrrbi4GK6uro3yWRUVFdDr9VAoFI3yeUR0e+yKJaK7lp+fj3HjxsHDwwOenp6IjY1Ffn5+rWVPnTqFv/3tb2jRogWcnZ3Ru3dvfPPNN0ZlDN3Ae/fuxQsvvAAvLy+oVCqMHTsWN27cuG19ahtjJ5PJEB8fjy1btqBLly5QKpXo3Lkztm/fftvrabVaLFiwAL169YKHhwfc3Nzw8MMPY9euXXW+55133kFQUBBcXFzQv39/HDt2zOh8Tk4OnnvuObRu3RpKpRJ+fn544okncOHCBaNy77//Pjp37gylUgl/f3/ExcXVeLYDBgxAly5dkJ6ejn79+sHV1RX//Oc/jcr88MMP6N69O5ydndGpUyds2rSpRp3z8/Mxbdo0BAQEQKlUIiQkBG+88Qb0er1Upvo4wmXLlqFdu3ZQKpU4ceLEbZ8jETUe+/6nLRHVqaCgAFevXjU6JpPJ4OXlBQAQQuCJJ57Avn37MHnyZHTs2BGbN29GbGxsjWsdP34cDz30EFq1aoU5c+bAzc0NGzZswIgRI/DVV19h5MiRRuXj4+Ph6emJRYsWITMzE6tWrcLFixexe/fuu5ocsW/fPmzatAlTp05Fs2bN8O6772LUqFHIysqS7qc2Go0GH3/8MZ599llMnDgRN2/exCeffIKoqCgcPHgQ3bt3Nyr/2Wef4ebNm4iLi0NpaSmWL1+ORx99FEePHoWvry8AYNSoUTh+/DhefPFFtGnTBnl5eUhOTkZWVhbatGkDoDKgLl68GJGRkZgyZYr0DA4dOoT9+/fDyclJ+sxr165hyJAheOaZZzB69GjpcwDgzJkzePrppzF58mTExsYiMTERf//737F9+3Y89thjACpb+Pr3748///wTL7zwAgIDA/Hzzz9j7ty5yM7OxrJly4zuMTExEaWlpZg0aRKUSiVatGjR4L8PIjIjQURUTWJiogBQ649SqZTKbdmyRQAQS5culY5VVFSIhx9+WAAQiYmJ0vGBAweKrl27itLSUumYXq8XDz74oGjfvn2Nz+7Vq5fQarXS8aVLlwoA4uuvv6637gsXLhR//VoDIBQKhTh79qx07PDhwwKAeO+99+q9XkVFhSgrKzM6duPGDeHr6yvGjx8vHTt//rwAIFxcXMQff/whHU9LSxMAxPTp06X3AhBvvvlmnZ+Zl5cnFAqFGDRokNDpdNLxFStWCABizZo10rH+/fsLAGL16tU1rhMUFCQAiK+++ko6VlBQIPz8/ESPHj2kY6+//rpwc3MTp0+fNnr/nDlzhFwuF1lZWUb3qFKpRF5eXp31JyLLYlcsEdVq5cqVSE5ONvr5/vvvpfPfffcdHB0dMWXKFOmYXC7Hiy++aHSd69evY+fOnXjqqadw8+ZNXL16FVevXsW1a9cQFRWFM2fO4M8//zR6z6RJk4xapaZMmQJHR0d89913d3UvkZGRaNeunfQ6LCwMKpUKv//+e73vk8vl0vgxvV6P69evo6KiAr1798avv/5ao/yIESPQqlUr6XWfPn0QHh4u1dvFxQUKhQK7d++us2v5xx9/hFarxbRp0+DgUPUVPXHiRKhUKmzbts2ovFKpxHPPPVfrtfz9/Y1aQw3d2r/99htycnIAABs3bsTDDz+M5s2bS383V69eRWRkJHQ6Hfbu3Wt0zVGjRqFly5Z1PjMisix2xRJRrfr06VPv5ImLFy/Cz88P7u7uRsc7dOhg9Prs2bMQQmD+/PmYP39+rdfKy8szCkTt27c3Ou/u7g4/P78a49DuVGBgYI1jzZs3v6Nxe59++inefvttnDp1CuXl5dLx2mYM/7XeAHDfffdhw4YNACpD2BtvvIFXXnkFvr6+eOCBBzBs2DCMHTsWarUaQOVzBWo+R4VCgbZt20rnDVq1alXn5IWQkJAaXdf33XcfgMoxc2q1GmfOnMGRI0fqDGt5eXlGr2u7byKyHgx2RGRWhgH4r776KqKiomotY+4lVORyea3HhRD1vu/zzz/HuHHjMGLECMycORM+Pj6Qy+VISEjAuXPn7qou06ZNw/Dhw7Flyxbs2LED8+fPR0JCAnbu3IkePXo0+HouLi53VQ8DvV6Pxx57DLNmzar1vCEImurziMi8GOyI6K4EBQUhJSUFhYWFRq12mZmZRuXatm0LAHByckJkZOQdXfvMmTN45JFHpNeFhYXIzs7G0KFDTVDzO/fll1+ibdu22LRpk1HL18KFC2stf+bMmRrHTp8+LU2KMGjXrh1eeeUVvPLKKzhz5gy6d++Ot99+G59//jmCgoIAVD5Hw7MDKmfonj9//o6fIVDVWlq97qdPnwYAqU7t2rVDYWFhg65LRNaLY+yI6K4MHToUFRUVWLVqlXRMp9PhvffeMyrn4+ODAQMG4IMPPkB2dnaN61y5cqXGsQ8//NCo23PVqlWoqKjAkCFDTHgHt2do6avespeWlobU1NRay2/ZssVovODBgweRlpYm1bu4uBilpaVG72nXrh2aNWuGsrIyAJXjARUKBd59912jz/3kk09QUFCA6OjoO67/5cuXsXnzZum1RqPBZ599hu7du0tdv0899RRSU1OxY8eOGu/Pz89HRUXFHX8eEVkeW+yIqFbff/89Tp06VeP4gw8+iLZt22L48OF46KGHMGfOHFy4cEFaI62goKDGe1auXIm+ffuia9eumDhxItq2bYvc3Fykpqbijz/+wOHDh43Ka7VaDBw4EE899RQyMzPx/vvvo2/fvnj88cfNdr+1GTZsGDZt2oSRI0ciOjoa58+fx+rVq9GpUycUFhbWKB8SEoK+fftiypQpKCsrw7Jly+Dl5SV1c54+fVq6r06dOsHR0RGbN29Gbm4unnnmGQBAy5YtMXfuXCxevBiDBw/G448/Lj2D+++/H6NHj77j+t93332YMGECDh06BF9fX6xZswa5ublITEyUysycORPffPMNhg0bhnHjxqFXr14oKirC0aNH8eWXX+LChQvw9va+xydJRI3GonNyicjq1LfcCf6yjMm1a9fEmDFjhEqlEh4eHmLMmDHit99+q1FOCCHOnTsnxo4dK9RqtXBychKtWrUSw4YNE19++WWNz96zZ4+YNGmSaN68uXB3dxcxMTHi2rVrt617XcudxMXF1SgbFBQkYmNj672eXq8XS5YsEUFBQUKpVIoePXqIrVu3itjYWBEUFCSVMywF8uabb4q3335bBAQECKVSKR5++GFx+PBhqdzVq1dFXFycCA0NFW5ubsLDw0OEh4eLDRs21PjsFStWiNDQUOHk5CR8fX3FlClTxI0bN4zK9O/fX3Tu3LnWugcFBYno6GixY8cOERYWJpRKpQgNDRUbN26sUfbmzZti7ty5IiQkRCgUCuHt7S0efPBB8dZbb0nLzlS/RyKyXjIhbjN6mIiokSQlJeG5557DoUOHmtx2ZkRE1oBj7IiIiIhsBIMdERERkY1gsCMiIiKyERxjR0RERGQj2GJHREREZCMY7IiIiIhshM0uUKzX63H58mU0a9asxibYRERERJYghMDNmzfh7+8PBwfTt6/ZbLC7fPkyAgICLF0NIiIiohouXbqE1q1bm/y6NhvsmjVrBqDywalUKgvXhoiIiKhyz+aAgAApp5iazQY7Q/erSqVisCMiIiKrYq5hYpw8QURERGQjGOyIiIiIbASDHREREZGNYLAjIiIishEMdkREREQ2gsGOiIiIyEYw2BERERHZCAY7IiIiIhvBYEdERERkIxjsiIiIiGwEgx0RERGRjWCwIyIiIrIRDHZERERENoLBjoiIiMhGMNgRERER2QgGOyIiIiIbwWBHREREZCMY7Jo4mUyGLVu2WLoaREREZAUY7ExEJpPV+7No0aI633vhwgXIZDJkZGQ0Sl3//e9/48EHH4Srqys8PT1rLZOVlYXo6Gi4urrCx8cHM2fOREVFhVGZ3bt3o2fPnlAqlQgJCUFSUlKN66xcuRJt2rSBs7MzwsPDcfDgQaPzpaWliIuLg5eXF9zd3TFq1Cjk5ubWW38hBBYsWAA/Pz+4uLggMjISZ86cadAzICIiskUMdiaSnZ0t/Sxbtgwqlcro2KuvvmrpKkq0Wi3+/ve/Y8qUKbWe1+l0iI6Ohlarxc8//4xPP/0USUlJWLBggVTm/PnziI6OxiOPPIKMjAxMmzYNzz//PHbs2CGV+eKLLzBjxgwsXLgQv/76K7p164aoqCjk5eVJZaZPn45vv/0WGzduxJ49e3D58mU8+eST9dZ/6dKlePfdd7F69WqkpaXBzc0NUVFRKC0tvccnQ0RE1MQJG1VQUCAAiIKCgkb/7MTEROHh4SG91ul0YvHixaJVq1ZCoVCIbt26ie+//146D8Dop3///kIIIQ4ePCgiIyOFl5eXUKlUol+/fiI9Pd3oswCIzZs3m6SeBt99951wcHAQOTk50rFVq1YJlUolysrKhBBCzJo1S3Tu3NnofU8//bSIioqSXvfp00fExcUZPQd/f3+RkJAghBAiPz9fODk5iY0bN0plTp48KQCI1NTUWuus1+uFWq0Wb775pnQsPz9fKJVK8b///a8Bd09ERNT4zJ1P2GLXCJYvX463334bb731Fo4cOYKoqCg8/vjjUvehoXvyxx9/RHZ2NjZt2gQAuHnzJmJjY7Fv3z4cOHAA7du3x9ChQ3Hz5s06P2vAgAEYN27cPdU3NTUVXbt2ha+vr3QsKioKGo0Gx48fl8pERkYavS8qKgqpqakAKlsF09PTjco4ODggMjJSKpOeno7y8nKjMqGhoQgMDJTK/NX58+eRk5Nj9B4PDw+Eh4fX+R4iIiJ74WjpCtiDt956C7Nnz8YzzzwDAHjjjTewa9cuLFu2DCtXrkTLli0BAF5eXlCr1dL7Hn30UaPrfPjhh/D09MSePXswbNiwWj8rMDAQfn5+91TfnJwco1AHQHqdk5NTbxmNRoOSkhLcuHEDOp2u1jKnTp2SrqFQKGqM8/P19ZU+p7a6Va/PnbyHiIjIXjDYmZlGo8Hly5fx0EMPGR1/6KGHcPjw4Xrfm5ubi3nz5mH37t3Iy8uDTqdDcXExsrKy6nzPZ599ZpJ6ExERUdPDrlgrFhsbi4yMDCxfvhw///wzMjIy4OXlBa1Wa9bPVavVNWamGl4bWhTrKqNSqeDi4gJvb2/I5fJay1S/hlarRX5+fp1laqtb9frcyXuIiIjsBYOdmalUKvj7+2P//v1Gx/fv349OnToBABQKBYDK2ah/LfPSSy9h6NCh6Ny5M5RKJa5evWr2OkdERODo0aNGs1eTk5OhUqmkOkdERCAlJcXofcnJyYiIiABQeU+9evUyKqPX65GSkiKV6dWrF5ycnIzKZGZmIisrSyrzV8HBwVCr1Ubv0Wg0SEtLq/M9RERE9oJdsY1g5syZWLhwIdq1a4fu3bsjMTERGRkZWLt2LQDAx8cHLi4u2L59O1q3bg1nZ2d4eHigffv2+O9//4vevXtDo9Fg5syZcHFxqfezxo4di1atWiEhIaHOMllZWbh+/TqysrKg0+mk9fNCQkLg7u6OQYMGoVOnThgzZgyWLl2KnJwczJs3D3FxcVAqlQCAyZMnY8WKFZg1axbGjx+PnTt3YsOGDdi2bZv0OTNmzEBsbCx69+6NPn36YNmyZSgqKsJzzz0HoHLSw4QJEzBjxgy0aNECKpUKL774IiIiIvDAAw9I1wkNDUVCQgJGjhwJmUyGadOm4V//+hfat2+P4OBgzJ8/H/7+/hgxYsTd/PUQERHZDrPMtbUC1rbcyaJFi0SrVq2Ek5NTjeVOhBDio48+EgEBAcLBwUFa7uTXX38VvXv3Fs7OzqJ9+/Zi48aNIigoSLzzzjvS+/CX5U769+8vYmNj661fbGxsjSVWAIhdu3ZJZS5cuCCGDBkiXFxchLe3t3jllVdEeXm50XV27dolunfvLhQKhWjbtq1ITEys8VnvvfeeCAwMFAqFQvTp00ccOHDA6HxJSYmYOnWqaN68uXB1dRUjR44U2dnZRmUAGF1br9eL+fPnC19fX6FUKsXAgQNFZmZmvfdMRERkDcydT2RCCGGxVGlGGo0GHh4eKCgogEqlsnR1iIiIiMyeTzjGjoiIiMhGMNgRERER2QgGOyIiG/DN4cuY/kUGSst1ty9MRDaLwY6IyAa8tSMTm3/7E4cuXLd0VYjIghjsCADQpk0bLFu2rEHv2b9/P7p27QonJycuNUJkQaXlOly6UQwAKCqrsHBtiMiSGOxMaNy4cZDJZJg8eXKNc3FxcZDJZBg3blzjV+wOHDp0CJMmTWrQe2bMmIHu3bvj/PnzSEpKMk/FiOi2fr9SBMP6BiXsiiWyawx2JhYQEID169ejpKREOlZaWop169YhMDDQgjWrX8uWLeHq6tqg95w7dw6PPvooWrduDU9PT/NUjIhu6+yVQunPxVoGOyJ7xmBnYj179kRAQAA2bdokHdu0aRMCAwPRo0cPo7J6vR4JCQkIDg6Gi4sLunXrhi+//FI6r9PpMGHCBOl8hw4dsHz5cqNrjBs3DiNGjMBbb70FPz8/eHl5IS4uDuXl5Q2q91+7YmUyGT7++GOMHDkSrq6uaN++Pb755hsAwIULFyCTyXDt2jWMHz8eMpmMLXZEFnQuryrYlTDYEdk1BjszGD9+PBITE6XXa9askbbRqi4hIQGfffYZVq9ejePHj2P69OkYPXo09uzZA6Ay+LVu3RobN27EiRMnsGDBAvzzn//Ehg0bjK6za9cunDt3Drt27cKnn36KpKQko6C1aNEitGnTpsH3sXjxYjz11FM4cuQIhg4dipiYGFy/fh0BAQHIzs6GSqXCsmXLkJ2djaeffrrB1yci06jeYsdZsUT2jcHODEaPHo19+/bh4sWLuHjxIvbv34/Ro0cblSkrK8OSJUuwZs0aREVFoW3bthg3bhxGjx6NDz74AADg5OSExYsXo3fv3ggODkZMTAyee+65GsGuefPmWLFiBUJDQzFs2DBER0cjJSVFOu/t7Y127do1+D7GjRuHZ599FiEhIViyZAkKCwtx8OBByOVyqNVqyGQyeHh4QK1W33YPWyIyn+otduyKJbJvjpaugC1q2bIloqOjkZSUBCEEoqOj4e3tbVTm7NmzKC4uxmOPPWZ0XKvVGnXZrly5EmvWrEFWVhZKSkqg1WrRvXt3o/d07twZcrlceu3n54ejR49Kr+Pj4xEfH9/g+wgLC5P+7ObmBpVKhby8vAZfh4jMR6cX+P1qkfSakyeI7BuDnZmMHz9eClMrV66scb6wsPJf2Nu2bUOrVq2MzimVSgDA+vXr8eqrr+Ltt99GREQEmjVrhjfffBNpaWlG5Z2cnIxey2Qy6PX6e74Hc12XiEznjxvF0FZU/XfJrlgi+8ZgZyaDBw+GVquFTCZDVFRUjfOdOnWCUqlEVlYW+vfvX+s19u/fjwcffBBTp06Vjp07d85sdSaipudstW5YgF2xRPauQWPsEhIScP/996NZs2bw8fHBiBEjkJmZaVRmwIABkMlkRj9/XdctKysL0dHRcHV1hY+PD2bOnImKCuNFNXfv3o2ePXtCqVQiJCSkyc26lMvlOHnyJE6cOGHUTWrQrFkzvPrqq5g+fTo+/fRTnDt3Dr/++ivee+89fPrppwCA9u3b45dffsGOHTtw+vRpzJ8/H4cOHWpwXVasWIGBAwfe8z0RkfU5d8U42HFWLJF9a1Cw27NnD+Li4nDgwAEkJyejvLwcgwYNQlFRkVG5iRMnIjs7W/pZunSpdE6n0yE6OhparRY///yzNItzwYIFUpnz588jOjoajzzyCDIyMjBt2jQ8//zz2LFjxz3ebuNSqVRQqVR1nn/99dcxf/58JCQkoGPHjhg8eDC2bduG4OBgAMALL7yAJ598Ek8//TTCw8Nx7do1o9a7O3X16lW29BHZKEOLXZBX5TqUHGNHZN9kQhjWK2+4K1euwMfHB3v27EG/fv0AVLbYde/evc7tqb7//nsMGzYMly9fhq+vLwBg9erVmD17Nq5cuQKFQoHZs2dj27ZtOHbsmPS+Z555Bvn5+di+ffsd1U2j0cDDwwMFBQX1hisioqbsyff349esfER39cO2o9noHdQcX0550NLVIqI6mDuf3NNyJwUFBQCAFi1aGB1fu3YtvL290aVLF8ydOxfFxcXSudTUVHTt2lUKdQAQFRUFjUaD48ePS2UiIyONrhkVFYXU1NQ661JWVgaNRmP0Q0Rky4QQOHelsseka2sPAGyxI7J3dz15Qq/XY9q0aXjooYfQpUsX6fg//vEPBAUFwd/fH0eOHMHs2bORmZkp7cSQk5NjFOoASK9zcnLqLaPRaFBSUlLrmmkJCQlYvHjx3d4OEVGTc7VQi4KScshkQCe/yn/5M9gR2be7DnZxcXE4duwY9u3bZ3S8+kbyXbt2hZ+fHwYOHIhz587d1SK5d2ru3LmYMWOG9Fqj0SAgIMBsn0dEZGmG8XUBzV3h6Vq5PBEnTxDZt7vqio2Pj8fWrVuxa9cutG7dut6y4eHhACoX5AUAtVqN3NxcozKG12q1ut4yKpWqzh0OlEqlNFnhdpMWiIhsgWFGbIiPO1ycKmffs8WOyL41KNgJIRAfH4/Nmzdj586d0uzN+mRkZACo3A0BACIiInD06FGjHQySk5OhUqnQqVMnqUz1LbEMZSIiIhpS3SZt0aJFNXaYuJ3du3dDJpMhPz/fLHUiIutiaLEL8XGHi+JWsGOLHZFda1Cwi4uLw+eff45169ahWbNmyMnJQU5ODkpKSgBULp77+uuvIz09HRcuXMA333yDsWPHol+/ftL2VIMGDUKnTp0wZswYHD58GDt27MC8efMQFxcn7bgwefJk/P7775g1axZOnTqF999/Hxs2bMD06dNNfPumdeXKFUyZMgWBgYFQKpVQq9WIiorC/v37pTIymQxbtmyxXCWJyGYYWuzatXSTWuzKKvTQ6e96sQMiauIaNMZu1apVACqXNKkuMTER48aNg0KhwI8//ohly5ahqKgIAQEBGDVqFObNmyeVlcvl2Lp1K6ZMmYKIiAi4ubkhNjYWr732mlQmODgY27Ztw/Tp07F8+XK0bt0aH3/8ca07OFiTUaNGQavV4tNPP0Xbtm2Rm5uLlJQUXLt2zdJVIyIbdK6WFjugclsxNyU3FiKyS8JGFRQUCACioKCgUT7vxo0bAoDYvXt3nWWCgoIEAOknKChIOpeQkCB8fHyEu7u7GD9+vJg9e7bo1q1bvZ+5bds20b59e+Hs7CwGDBggEhMTBQBx48YNqcxPP/0k+vbtK5ydnUXr1q3Fiy++KAoLC4UQQsydO1f06dOnxnXDwsLE4sWL7+i+ExMThYeHh9i+fbsIDQ0Vbm5uIioqSly+fFkqc/DgQREZGSm8vLyESqUS/fr1E+np6UbXASA++ugjMWLECOHi4iJCQkLE119/fUd1ILJHhaXlImj2VhE0e6u4UVQmdDq99DpPU2rp6hFRHcydT+5pHTuq4u7uDnd3d2zZsgVlZWW1ljFsB5aYmIjs7Gzp9YYNG7Bo0SIsWbIEv/zyC/z8/PD+++/X+3mXLl3Ck08+ieHDhyMjIwPPP/885syZY1Tm3LlzGDx4MEaNGoUjR47giy++wL59+xAfHw8AiImJwcGDB412pTh+/DiOHDmCf/zjHwAqu45vt51bcXEx3nrrLfz3v//F3r17kZWVhVdffVU6f/PmTcTGxmLfvn04cOAA2rdvj6FDh+LmzZtG11m8eDGeeuopHDlyBEOHDkVMTAyuX79e72cT2StDN6y3uwKergo4OMjg7FT5lV7KCRRE9ssscdEKNHaLnRBCfPnll6J58+bC2dlZPPjgg2Lu3Lni8OHDRmUAiM2bNxsdi4iIEFOnTjU6Fh4eXm+L3dy5c0WnTp2Mjs2ePduoxW7ChAli0qRJRmV++ukn4eDgIEpKSoQQQnTr1k289tprRtcNDw+XXnfo0EFs2rSpznoYWgnPnj0rHVu5cqXw9fWt8z06nU40a9ZMfPvtt9IxAGLevHnS68LCQgFAfP/993Veh8iebfr1kgiavVU8tfpn6ViP134QQbO3iswcjQVrRkT1YYtdEzJq1ChcvnwZ33zzDQYPHozdu3ejZ8+et23xOnnypLQsjMHtZgDfyXsOHz6MpKQkqTXR3d0dUVFR0Ov1OH/+PIDKVrt169YBqJz1/L///Q8xMTHSNU6dOoWRI0fWWxdXV1ejNQr9/PyMZj3n5uZi4sSJaN++PTw8PKBSqVBYWIisrCyj6xgm2ACAm5sbVCqV0XWIqIphRmw7H3fpmGECRTFnxhLZLY6uNTFnZ2c89thjeOyxxzB//nw8//zzWLhwIcaNG9fodSksLMQLL7yAl156qca5wMBAAMCzzz6L2bNn49dff0VJSQkuXbqEp59+ukGf4+TkZPRaJpNBVNuCODY2FteuXcPy5csRFBQEpVKJiIgIaLXa215Hr9c3qC5E9kJa6qRlVbAzdMVyyRMi+8VgZ2adOnUyWt7EyckJOp3xl27Hjh2RlpaGsWPHSscOHDhQ73U7duyIb775xujYX9/Ts2dPnDhxAiEhIXVep3Xr1ujfvz/Wrl2LkpISPPbYY/Dx8bndbTXI/v378f7772Po0KEAKscHXr161aSfQWRvDHvEhlRrsXNVVH6lc4wdkf1iV6yJXLt2DY8++ig+//xzHDlyBOfPn8fGjRuxdOlSPPHEE1K5Nm3aICUlBTk5Obhx4wYA4OWXX8aaNWuQmJiI06dPY+HChTh+/LjR9Tdv3ozQ0FDp9eTJk3HmzBnMnDkTmZmZWLduXY0u39mzZ+Pnn39GfHw8MjIycObMGXz99dfS5AmDmJgYrF+/Hhs3bjTqhgWA0NBQbN68+Z6eTfv27fHf//4XJ0+eRFpaGmJiYurcQYSIbq9cp8eFq5XBjl2xRFQdg52JuLu7Izw8HO+88w769euHLl26YP78+Zg4cSJWrFghlXv77beRnJyMgIAA9OjRAwDw9NNPY/78+Zg1axZ69eqFixcvYsqUKUbXLygoQGZmpvQ6MDAQX331FbZs2YJu3bph9erVWLJkidF7wsLCsGfPHpw+fRoPP/wwevTogQULFsDf39+o3N/+9jdcu3YNxcXFGDFihNG5zMxMFBQU3NOz+eSTT3Djxg307NkTY8aMwUsvvWTyVkEie3LxWjEq9AKuCjn8PZyl484KbitGZO9kovpgKBui0Wjg4eGBgoIC7htLRDZlx/EcvPDfdHRt5YFvX+wrHZ/833RsP56D10d0wZgHgixYQyKqi7nzCVvsiIiaGGlGbEs3o+NV+8VWNHqdiMg6MNgRETUx1bcSq87ZyRDsOJucyF4x2BERNTGGXSf+GuxcOcaOyO4x2BERNSFCCGmpk3YtjYOdixO7YonsHYOdBQ0YMADTpk2zys85deoUHnjgATg7O6N79+5mqRcRNVyOphSFZRWQO8gQ5FXHGDu22BHZLQY7Exo3bhxkMhkmT55c41xcXBxkMpnRDhSbNm3C66+/bvZ63c3nLFy4EG5ubsjMzERKSoqZakZEDXUur7K1LsjLFQpH469wrmNHRAx2JhYQEID169ejpKREOlZaWop169ZJ23gZtGjRAs2aNTN7ne7mc86dO4e+ffsiKCgIXl5eZqoZETVU3s1SAIC/R81Fvg0tdtx5gsh+MdiZWM+ePREQEIBNmzZJxzZt2oTAwEBpQWKDv3aRtmnTBkuWLMH48ePRrFkzBAYG4sMPP7znOjX0c2QyGdLT0/Haa69BJpNh0aJF91wHIjKNolutce7KmjtCcvIEETHYmcH48eORmJgovV6zZg2ee+65O3rv22+/jd69e+O3337D1KlTMWXKFKMdJ/4qKSkJMpmswXWs73Oys7PRuXNnvPLKK8jOzsarr77a4OsTkXkUlVVOjHBVymucc2ZXLJHdY7Azg9GjR2Pfvn24ePEiLl68iP3792P06NF39N6hQ4di6tSpCAkJwezZs+Ht7Y1du3bVWd7DwwMdOnRocB3r+xy1Wg1HR0e4u7tDrVbD3d39NlcjosZiCHa1tdhVzYplsCOyVzW/GeietWzZEtHR0UhKSoIQAtHR0fD29r6j94aFhUl/lslkUKvVyMvLq7P8yJEjMXLkyAbXsaGfQ0TWofBWsHOrpyuWY+yI7BeDnZmMHz8e8fHxAICVK1fe8fucnJyMXstkMuj1pl9FvrE+h4hMq74WO3bFEhGDnZkMHjwYWq0WMpkMUVFRlq4OEdkIw+QJN0XNMXZcx46IOMbOTORyOU6ePIkTJ05ALq/5BXy3Bg4ciBUrVkivN2/ejNDQUJNdn4isW9XkCXbFElFNbLEzI5VKZfJrnjt3DlevXpVeFxQU1Dtrlohsy51MnijXCZTr9HCS89/uRPZGJoQQlq6EOWg0Gnh4eKCgoMAsAYuIyBKGLP8JJ7M1+HR8H/S/r6XRubIKHTrM2w4AOLJoEFTOTrVdgogsyNz5hP+cIyJqQqpa7GoO8VDIHeBwa1nLUk6gILJLDHZERE1Isbbu5U5kMhn3iyWycwx2RERNiLSOnaL2IdIut45zZiyRfWKwIyJqIip0epSWV643WVuLHQC4KCq/1tliR2SfGOwIADBu3DiMGDGiQe/JycnBY489Bjc3N3h6epqlXkRUpahaWHOrZYwdUDUzlkueENknBjsTGjduHGQyGSZPnlzjXFxcHGQyGcaNG9f4FbsDy5cvR1JSUoPe88477yA7OxsZGRk4ffq0eSpGRBLD+DonuQxKxzqCnaErli12RHaJwc7EAgICsH79epSUlEjHSktLsW7dOgQGBlqwZvXz8PBocKvbuXPn0KtXL7Rv3x4+Pj7mqRgRSaTFiesYXwcALk63umLZYkdklxjsTKxnz54ICAjApk2bpGObNm1CYGAgevToYVR2+/bt6Nu3Lzw9PeHl5YVhw4bh3Llz0vnPPvsM7u7uOHPmjHRs6tSpCA0NRXFxsUnr/deu2AEDBuCll17CrFmz0KJFC6jVaixatEg636ZNG3z11Vf47LPPrLolksiWFJZVhrXaFic2kLpi2WJHZJcY7Mxg/PjxSExMlF6vWbMGzz33XI1yRUVFmDFjBn755RekpKTAwcEBI0eOhF5fOTh67NixGDp0KGJiYlBRUYFt27bh448/xtq1a+Hq6lrrZ1+4cAEymQy7d+++5/v49NNP4ebmhrS0NCxduhSvvfYakpOTAQCHDh3C4MGD8dRTTyE7OxvLly+/588jovoZWuzqGl8HVLXmcVYskX3ilmJmMHr0aMydOxcXL14EAOzfvx/r16+vEbZGjRpl9HrNmjVo2bIlTpw4gS5dugAAPvjgA4SFheGll17Cpk2bsGjRIvTq1avOz3ZyckKHDh3qDH4NERYWhoULFwIA2rdvjxUrViAlJQWPPfYYWrZsCaVSCRcXF6jV6nv+LCK6PWmpk3pa7Jy5jh2RXWOLnRm0bNkS0dHRSEpKQmJiIqKjo+Ht7V2j3JkzZ/Dss8+ibdu2UKlUaNOmDQAgKytLKtO8eXN88sknWLVqFdq1a4c5c+bU+9mtWrXCqVOn0KdPn1rPZ2Vlwd3dXfpZsmRJndcKCwszeu3n54e8vLx6P5+IzMcweaLerthby52wxY7IPrHFzkzGjx+P+Ph4AMDKlStrLTN8+HAEBQXho48+gr+/P/R6Pbp06QKtVmtUbu/evZDL5cjOzkZRURGaNWt21/Xy9/dHRkaG9LpFixZ1lnVyMt5nUiaTSd3ERNT4DGPsXBW374rlcidE9oktdmYyePBgaLValJeXIyoqqsb5a9euITMzE/PmzcPAgQPRsWNH3Lhxo0a5n3/+GW+88Qa+/fZbuLu7S2Hxbjk6OiIkJET6qS/YEZF1KWpQV2xFo9SJiKwLg52ZyOVynDx5EidOnIBcXvNf182bN4eXlxc+/PBDnD17Fjt37sSMGTOMyty8eRNjxozBSy+9hCFDhmDt2rX44osv8OWXX0pl5s6di7Fjx0qv//zzT4SGhuLgwYPmuzkisghDsKuvK9bQmleiZes6kT1isDMjlUoFlUpV6zkHBwesX78e6enp6NKlC6ZPn44333zTqMzLL78MNzc3aRxc165dsWTJErzwwgv4888/AQDZ2dlGY/LKy8uRmZlp8uVQiMjyim51xdbXYsedJ4jsm0wIISxdCXPQaDTw8PBAQUFBneGKiKgpmf3lEXzxyyXMjOqAuEdCai2z4dAlzPrqCB7p0BKJz9U+iYqILMfc+YQtdkRETUSh1rDzRN2TJ1wMXbFssSOySwx2RERNxJ1MnjB0xXKvWCL7xGBHRNRE3MnkCbbYEdk3BjsioibijiZPMNgR2TUGOyKiJqJI2nminjF27IolsmsMdkRETYShK9awu0RtGOyI7BuDHRFRE1HYkAWKy3Ww0dWsiKgeDHZERE1AhU6P0vLK3STq3VLsVrDTC6CsgrtPENkbBjsioiaguNpkCLc7GGMHcPcJInvEYEdE1AQYxtc5OsigkNf91e0kd4CTXAaAM2OJ7BGDHRFRE1B9cWKZTFZvWedbrXbFnEBBZHcY7IiImoDCW2vY1TdxwkCaQMFgR2R3GOyIiJqAYqnFru7xdQaGcXYcY0dkfxjsiIiagMI72CfWgF2xRPaLwY6IqAkw7DrhVs/ixAau3FaMyG4x2BERNQGF0j6xd9AVyzF2RHarQcEuISEB999/P5o1awYfHx+MGDECmZmZRmVKS0sRFxcHLy8vuLu7Y9SoUcjNzTUqk5WVhejoaLi6usLHxwczZ85ERUWFUZndu3ejZ8+eUCqVCAkJQVJS0t3dIRGRDShqQFestK0YW+yI7E6Dgt2ePXsQFxeHAwcOIDk5GeXl5Rg0aBCKioqkMtOnT8e3336LjRs3Ys+ePbh8+TKefPJJ6bxOp0N0dDS0Wi1+/vlnfPrpp0hKSsKCBQukMufPn0d0dDQeeeQRZGRkYNq0aXj++eexY8cOE9wyEVHTU3wH24kZuNzqrmWLHZH9uf03RDXbt283ep2UlAQfHx+kp6ejX79+KCgowCeffIJ169bh0UcfBQAkJiaiY8eOOHDgAB544AH88MMPOHHiBH788Uf4+vqie/fueP311zF79mwsWrQICoUCq1evRnBwMN5++20AQMeOHbFv3z688847iIqKMtGtExE1HVVdsXfSYlf5b3a22BHZn3saY1dQUAAAaNGiBQAgPT0d5eXliIyMlMqEhoYiMDAQqampAIDU1FR07doVvr6+UpmoqChoNBocP35cKlP9GoYyhmvUpqysDBqNxuiHiMhWSF2xituPsXNlix2R3brrYKfX6zFt2jQ89NBD6NKlCwAgJycHCoUCnp6eRmV9fX2Rk5Mjlake6gznDefqK6PRaFBSUlJrfRISEuDh4SH9BAQE3O2tERFZnUJtw5c7YYsdkf2562AXFxeHY8eOYf369aasz12bO3cuCgoKpJ9Lly5ZukpERCZzN5MnuI4dkf1p0Bg7g/j4eGzduhV79+5F69atpeNqtRparRb5+flGrXa5ublQq9VSmYMHDxpdzzBrtnqZv86kzc3NhUqlgouLS611UiqVUCqVd3M7RERWr/guthTjzhNE9qdBLXZCCMTHx2Pz5s3YuXMngoODjc736tULTk5OSElJkY5lZmYiKysLERERAICIiAgcPXoUeXl5Upnk5GSoVCp06tRJKlP9GoYyhmsQEdmbBu08wXXsiOxWg1rs4uLisG7dOnz99ddo1qyZNCbOw8MDLi4u8PDwwIQJEzBjxgy0aNECKpUKL774IiIiIvDAAw8AAAYNGoROnTphzJgxWLp0KXJycjBv3jzExcVJLW6TJ0/GihUrMGvWLIwfPx47d+7Ehg0bsG3bNhPfPhFR01C188Sd7xVbzBY7IrvToBa7VatWoaCgAAMGDICfn5/088UXX0hl3nnnHQwbNgyjRo1Cv379oFarsWnTJum8XC7H1q1bIZfLERERgdGjR2Ps2LF47bXXpDLBwcHYtm0bkpOT0a1bN7z99tv4+OOPudQJEdmthoyxk7pi2WJHZHca1GInhLhtGWdnZ6xcuRIrV66ss0xQUBC+++67eq8zYMAA/Pbbbw2pHhGRzSpqwBg77jxBZL+4VywRkZXT6YUU0hqy3EmxtuI2JYnI1jDYERFZuaJqAc31jhYoNsyK1ZutTkRknRjsiIisnGF8naODDErH239tuyjYYkdkrxjsiIisXPWJEzKZ7LblOcaOyH4x2BERWbmGTJwAqlrsSsv10OtvP+mNiGwHgx0RkZWrarG7/fg6oKrFDgBKK9hqR2RPGOyIiKycYdcJV8UdtthVC3bcfYLIvjDYERFZOcOs2DvtinWoNsmC4+yI7AuDHRGRlTOMsbvTrligapwdW+yI7AuDHRGRlWvIdmIGrpwZS2SXGOyIiKycIdjdaVcsADizxY7ILjHYERFZucJbXbF3OnkCqJpAUcwWOyK7wmBHRGTlqlrs7nyMnbStGFvsiOwKgx0RkZUzzIptyBg7Z46xI7JLDHZERFburiZPSPvFMtgR2RMGOyIiK9fQLcWAqjF2pWyxI7IrDHZERFauaueJhq9jxxY7IvvCYEdEZOUauvMEALg4VZblGDsi+8JgR0Rk5ap2nmhAsFPc2lKMLXZEdoXBjojIyt3NAsWGMXYMdkT2hcGOiMiK6fRC6k5t2Bg7dsUS2SMGOyIiK2YYXwc0sCuW69gR2SUGOyIiK1Z8a3ydo4MMSsc7/8rmGDsi+8RgR0RkxQqrLU4sk8nu+H2cFUtknxjsiIis2N1MnACq1rFjix2RfWGwIyKyYkV3sTgxwDF2RPaKwY6IyIoV3sU+sUBVEGSLHZF9YbAjIrJihi3BGtoV68wWOyK7xGBHRGTFqlrsGtYVyxY7IvvEYEdEZMWK7rIr1jDGTqvTo0KnN3m9iMg6MdgREVkxKdgp7m5WLMDuWCJ7wmBHRGTFCm8tUNzQFjulowMMy94x2BHZDwY7IiIrVqw1rGPXsDF2MplM6o4t1bIrlsheMNgREVmxu13uBKgaZ1dcXnGbkkRkKxjsiIis2N1OngC4+wSRPWKwIyKyYkWGMXYNnDwBcPcJInvEYEdEZMWKtHe3jh3AFjsie8RgR0RkxQxdsQ3deQJgix2RPWKwIyKyYne73AnAFjsie8RgR0Rkxe52gWKg2rZibLEjshsMdkREVkqnF1Iou5sxds5ObLEjsjcMdkREVsqwODFwb+vYscWOyH4w2BERWSnDUieODjIoHRv+de3KMXZEdofBjojISlXfdUJm2Pi1AdhiR2R/GOyIiKxU1cSJho+vAwDnW+8rZosdkd1gsCMislJVixM3fHwdALiyxY7I7jDYERFZqaJ7WMMOqFrHrpQtdkR2g8GOiMhK3cuuE0DVcifsiiWyHwx2RERWqmryxN2NsXO9tahxMbtiiewGgx0RkZW6l10nAMDZqfIrvozBjshuMNgREVmpIu29jbEzdMWWVehNVicism4MdkREVqqo7N5mxTo73po8wRY7IrvBYEdEZKWqJk/c5Tp2t7piGeyI7AeDHRGRlSq81xY7J0OLHbtiiewFgx0RkZW618kTSkOLXYUOQgiT1YuIrBeDHRGRlTLV5AkhAK2OrXZE9qDBwW7v3r0YPnw4/P39IZPJsGXLFqPz48aNg0wmM/oZPHiwUZnr168jJiYGKpUKnp6emDBhAgoLC43KHDlyBA8//DCcnZ0REBCApUuXNvzuiIiasKJ7XMfOMHkCYHcskb1ocLArKipCt27dsHLlyjrLDB48GNnZ2dLP//73P6PzMTExOH78OJKTk7F161bs3bsXkyZNks5rNBoMGjQIQUFBSE9Px5tvvolFixbhww8/bGh1iYiarHvdecJJLoODrPLPXMuOyD40+NtiyJAhGDJkSL1llEol1Gp1redOnjyJ7du349ChQ+jduzcA4L333sPQoUPx1ltvwd/fH2vXroVWq8WaNWugUCjQuXNnZGRk4D//+Y9RACQismU3S28FO+e7C3YymQzOTnIUa3VssSOyE2YZY7d79274+PigQ4cOmDJlCq5duyadS01NhaenpxTqACAyMhIODg5IS0uTyvTr1w8KhUIqExUVhczMTNy4ccMcVSYisjqGYNfM2emuryHNjK1gix2RPbi7fwbWY/DgwXjyyScRHByMc+fO4Z///CeGDBmC1NRUyOVy5OTkwMfHx7gSjo5o0aIFcnJyAAA5OTkIDg42KuPr6yuda968eY3PLSsrQ1lZmfRao9GY+taIiBpNablOmvDQ7C5b7ADA2ZFr2RHZE5MHu2eeeUb6c9euXREWFoZ27dph9+7dGDhwoKk/TpKQkIDFixeb7fpERI3JsIadTAa43+VyJwDXsiOyN2Zf7qRt27bw9vbG2bNnAQBqtRp5eXlGZSoqKnD9+nVpXJ5arUZubq5RGcPrusbuzZ07FwUFBdLPpUuXTH0rRESNRhpfp3CEg2EGxF1QOnFbMSJ7YvZg98cff+DatWvw8/MDAERERCA/Px/p6elSmZ07d0Kv1yM8PFwqs3fvXpSXl0tlkpOT0aFDh1q7YYHKCRsqlcroh4ioqbpZWvn9dy/dsAC3FSOyNw0OdoWFhcjIyEBGRgYA4Pz588jIyEBWVhYKCwsxc+ZMHDhwABcuXEBKSgqeeOIJhISEICoqCgDQsWNHDB48GBMnTsTBgwexf/9+xMfH45lnnoG/vz8A4B//+AcUCgUmTJiA48eP44svvsDy5csxY8YM0905EZEVu9cZsQZKwxi7CnbFEtmDBge7X375BT169ECPHj0AADNmzECPHj2wYMECyOVyHDlyBI8//jjuu+8+TJgwAb169cJPP/0EpVIpXWPt2rUIDQ3FwIEDMXToUPTt29dojToPDw/88MMPOH/+PHr16oVXXnkFCxYs4FInRGQ3qlrs7n5GLFB9jB1b7IjsQYP/KThgwIB69xzcsWPHba/RokULrFu3rt4yYWFh+OmnnxpaPSIim6CRljq5x67YW7tPcIFiIvvAvWKJiKyQKdawA6qPsWNXLJE9YLAjIrJChaZqsWNXLJFdYbAjIrJCppsVy50niOwJgx0RkRWSumKV9zgrll2xRHaFwY6IyArdLDPRrFhHdsUS2RMGOyIiK3TT5GPs2GJHZA8Y7IiIrJDG1LNiOcaOyC4w2BERWSFTT57gOnZE9oHBjojICpluuRNOniCyJwx2RERWyDDGTsXJE0TUAAx2RERWplynR8mtIOZ+j8udcB07IvvCYEdEZGUM3bAA4H6PXbFcx47IvjDYERFZGUM3rIuTHE7ye/ua5pZiRPaFwY6IyMpoTDQjFqg+xo4tdkT2gMGOiMjKmGpxYqBqViyXOyGyDwx2RERWprDMNIsTA5w8QWRvGOyIiKyMqRYnBqqCXblOQKcX93w9IrJuDHZERFbGHF2xACdQENkDBjsiIisjtdgpTdAVe2vyBMBgR2QPGOyIiKyMKVvsHBxkUNxaMqW0gjNjiWwdgx0RkZXRlJpu8gRQfZFittgR2ToGOyIiK2PKyRMAFykmsicMdkREVqZquRNTBTtuK0ZkLxjsiIiszE0Td8UaJlCUcS07IpvHYEdEZGXM1RVbxhY7IpvHYEdEZGVMOSsWqN4VyxY7IlvHYEdEZGVM3hXLbcWI7AaDHRGRFdHphcknTygdDbNi2RVLZOsY7IiIrIgh1AHsiiWihmOwIyKyIoZgp3B0kFra7lXVOnZssSOydQx2RERWxDAjVmWi1jqALXZE9oTBjojIihgmTrgrTRjsHDl5gsheMNgREVmRqjXsTDMjFuA6dkT2hMGOiMiKmHoNO4BdsUT2hMGOiMiKaMwS7AyTJxjsiGwdgx0RkRUxR1eskrNiiewGgx0RkRUxS1es462uWE6eILJ5DHZERFak0MTbiQHsiiWyJwx2RERWROqKNeVyJ+yKJbIbDHZERFaEs2KJ6F4w2BERWZGbZuyKLatgix2RrWOwIyKyIhppVqwZdp5gix2RzWOwIyKyIuboilWyK5bIbjDYERFZEbNsKebIyRNE9oLBjojISgghUFhmxskTFToIIUx2XSKyPgx2RERWolirg/5W7jJtV2xli50QgFbHVjsiW8ZgR0RkJQzj6+QOMrjcCmOmYGixA9gdS2TrGOyIiKzEzWozYmUymcmuq5A7wHC5Mk6gILJpDHZERFZCY4YZsQAgk8k4gYLITjDYERFZiartxEw3I9ag+gQKIrJdDHZERFbCHGvYGVTtF8tgR2TLGOyIiKyEOZY6MagKduyKJbJlDHZERFbCHIsTGygdufsEkT1gsCMishLsiiWie8VgR0RkJcwb7AyTJ9gVS2TLGOyIiKyExoxdsWyxI7IPDQ52e/fuxfDhw+Hv7w+ZTIYtW7YYnRdCYMGCBfDz84OLiwsiIyNx5swZozLXr19HTEwMVCoVPD09MWHCBBQWFhqVOXLkCB5++GE4OzsjICAAS5cubfjdERE1IWZtsbu1jh0XKCaybQ0OdkVFRejWrRtWrlxZ6/mlS5fi3XffxerVq5GWlgY3NzdERUWhtLRUKhMTE4Pjx48jOTkZW7duxd69ezFp0iTpvEajwaBBgxAUFIT09HS8+eabWLRoET788MO7uEUioqbBnJMnpK5YzoolsmkN/mfhkCFDMGTIkFrPCSGwbNkyzJs3D0888QQA4LPPPoOvry+2bNmCZ555BidPnsT27dtx6NAh9O7dGwDw3nvvYejQoXjrrbfg7++PtWvXQqvVYs2aNVAoFOjcuTMyMjLwn//8xygAEhHZEmm5EyUnTxDR3THpGLvz588jJycHkZGR0jEPDw+Eh4cjNTUVAJCamgpPT08p1AFAZGQkHBwckJaWJpXp168fFAqFVCYqKgqZmZm4ceOGKatMRGQ1GmVWLHeeILJpJv32yMnJAQD4+voaHff19ZXO5eTkwMfHx7gSjo5o0aKFUZng4OAa1zCca968eY3PLisrQ1lZmfRao9Hc490QETWuqmBnhnXs2BVLZBdsZlZsQkICPDw8pJ+AgABLV4mI6I4JIaqNsTPf5Al2xRLZNpMGO7VaDQDIzc01Op6bmyudU6vVyMvLMzpfUVGB69evG5Wp7RrVP+Ov5s6di4KCAunn0qVL935DRESNpKxCj3KdAMAtxYjo7pk02AUHB0OtViMlJUU6ptFokJaWhoiICABAREQE8vPzkZ6eLpXZuXMn9Ho9wsPDpTJ79+5FeXm5VCY5ORkdOnSotRsWAJRKJVQqldEPEVFTYVjDTiYD3BTmXKCYLXZEtqzBwa6wsBAZGRnIyMgAUDlhIiMjA1lZWZDJZJg2bRr+9a9/4ZtvvsHRo0cxduxY+Pv7Y8SIEQCAjh07YvDgwZg4cSIOHjyI/fv3Iz4+Hs888wz8/f0BAP/4xz+gUCgwYcIEHD9+HF988QWWL1+OGTNmmOzGiYisiWF8nbvSEQ4OMpNf39Bix3XsiGxbg/9Z+Msvv+CRRx6RXhvCVmxsLJKSkjBr1iwUFRVh0qRJyM/PR9++fbF9+3Y4OztL71m7di3i4+MxcOBAODg4YNSoUXj33Xel8x4eHvjhhx8QFxeHXr16wdvbGwsWLOBSJ0RkswpLzbfUCVDVYlfGLcWIbFqDv0EGDBgAIUSd52UyGV577TW89tprdZZp0aIF1q1bV+/nhIWF4aeffmpo9YiImiRzzogFOHmCyF7YzKxYIqKmzJwzYgFOniCyFwx2RERWwJyLEwPV17Fjix2RLWOwIyKyAhoz7hMLcOcJInvBYEdEZAXM3WJXNcaOXbFEtozBjojICph98gS7YonsAoMdEZEVKCxrnMkTZWyxI7JpDHZERFbA7F2xt4KdVqeHTl/3klVE1LQx2BERWQHzB7uqr/syTqAgslkMdkREVkBax05p3gWKAU6gILJlDHZERFbA3C12Dg4yKOScQEFk6xjsiIisgMbMs2IBLlJMZA8Y7IiIrIC5txQDuK0YkT1gsCMisjBthR5lFZVhy7zB7laLHSdPENksBjsiIgsrLKuQ/uyuNGOwk3afYLAjslUMdkREFmbohnVVyOEoN9/XMhcpJrJ9DHZERBZm7hmxBtxWjMj2MdgREVmYRpo4Yb4ZsUC1yRMcY0dksxjsiIgsrLFa7JSOnBVLZOsY7IiILMwQ7Mw5cQJgVyyRPWCwIyKysGuFZQAAb3elWT+H69gR2T4GOyIiC7sqBTuFWT9H6cgWOyJbx2BHRGRhV242cosdJ08Q2SwGOyIiC7taqAUAtGxm7mBX+ZXPdeyIbBeDHRGRhV1trDF23HmCyOYx2BERWVijd8Uy2BHZLAY7IiILqtDpcb24cbtiOSuWyHYx2BERWdD1Yi2EABxkQAs3M8+K5eQJIpvHYEdEZEGGbtgWbgrIHWRm/Sx2xRLZPgY7IiILMsyINff4OgBwdmRXLJGtY7AjIrKgq7da7Mw9vg5gix2RPWCwIyKyoMZa6gSoCnZlFWyxI7JVDHZERBZUtdSJeSdOANVnxbLFjshWMdgREVmQocWOXbFEZAoMdkREFtS4kycMwY5dsUS2isGOiMiCGmvXCaBaV2yFDkIIs38eETU+BjsiIgtqzK5YwwLFQgBaHVvtiGwRgx0RkYVU306sMVvsAHbHEtkqBjsiIgtpzO3EAEAhd4Ds1uYWZZxAQWSTGOyIiCykMbcTAwCZTMYJFEQ2jsGOiMhCGnNGrIGhO7asgi12RLaIwY6IyEIaczsxg6q17NhiR2SLGOyIiCzkSiNuJ2YgBTu22BHZJAY7IiILsUSLndKR24oR2TIGOyIiC7la2Hj7xBqwK5bItjHYERFZiGW6YtliR2TLGOyIiCzk6s3KWbGWmTzBYEdkixjsiIgs5KolWuwM69hVsCuWyBYx2BERWUBjbydmIK1jxxY7IpvEYEdEZAHXixp3OzEDdsUS2TYGOyIiCzBMnGjhpmyU7cQMOCuWyLYx2BERWUDVdmKN11oHAErOiiWyaQx2REQWcMUCixMD1SdPMNgR2SIGOyIiCzDMiG3ZiBMnAHbFEtk6BjsiIgswbCfm3dgtduyKJbJpDHZERBZgie3EALbYEdk6BjsiIgswzIpt9DF2hnXsOMaOyCaZPNgtWrQIMpnM6Cc0NFQ6X1pairi4OHh5ecHd3R2jRo1Cbm6u0TWysrIQHR0NV1dX+Pj4YObMmaioqDB1VYmILMawnVhjLk4MVJs8wa5YIpvkaI6Ldu7cGT/++GPVhzhWfcz06dOxbds2bNy4ER4eHoiPj8eTTz6J/fv3AwB0Oh2io6OhVqvx888/Izs7G2PHjoWTkxOWLFlijuoSETU6S2wnBrArlsjWmSXYOTo6Qq1W1zheUFCATz75BOvWrcOjjz4KAEhMTETHjh1x4MABPPDAA/jhhx9w4sQJ/Pjjj/D19UX37t3x+uuvY/bs2Vi0aBEUisYdj0JEZGqW2k4M4Dp2RLbOLGPszpw5A39/f7Rt2xYxMTHIysoCAKSnp6O8vByRkZFS2dDQUAQGBiI1NRUAkJqaiq5du8LX11cqExUVBY1Gg+PHj5ujukREjcpS24kB1VrsOMaOyCaZvMUuPDwcSUlJ6NChA7Kzs7F48WI8/PDDOHbsGHJycqBQKODp6Wn0Hl9fX+Tk5AAAcnJyjEKd4bzhXF3KyspQVlYmvdZoNCa6IyIi07LUdmJA9TF27IolskUmD3ZDhgyR/hwWFobw8HAEBQVhw4YNcHFxMfXHSRISErB48WKzXZ+IyFQMu0409lInANexI7J1Zl/uxNPTE/fddx/Onj0LtVoNrVaL/Px8ozK5ubnSmDy1Wl1jlqzhdW3j9gzmzp2LgoIC6efSpUumvREiIhMx7BPb2EudAFVdsWVssSOySWYPdoWFhTh37hz8/PzQq1cvODk5ISUlRTqfmZmJrKwsREREAAAiIiJw9OhR5OXlSWWSk5OhUqnQqVOnOj9HqVRCpVIZ/RARWSNLbScGVAU7rU4PnV40+ucTkXmZvCv21VdfxfDhwxEUFITLly9j4cKFkMvlePbZZ+Hh4YEJEyZgxowZaNGiBVQqFV588UVERETggQceAAAMGjQInTp1wpgxY7B06VLk5ORg3rx5iIuLg1LZ+F+CRESmZqntxICqrligcpFiV4VZFkcgIgsx+X/Rf/zxB5599llcu3YNLVu2RN++fXHgwAG0bNkSAPDOO+/AwcEBo0aNQllZGaKiovD+++9L75fL5di6dSumTJmCiIgIuLm5ITY2Fq+99pqpq0pEZBFXLNhip7w1eQKonEDhyhWkiGyKyYPd+vXr6z3v7OyMlStXYuXKlXWWCQoKwnfffWfqqhERWQVpceJmjZ+q5A4yOMllKNcJTqAgskHcK5aIqJFZajsxA24rRmS7GOyIiBrZFQttJ2ag5LZiRDaLwY6IqBFV6PS4UWy55U6AamvZcfcJIpvDYEdE1IiqbyfW3EIzF6RtxdgVS2RzGOyIiBqRJbcTMzC02HGRYiLbw2BHRNSIDNuJWaobFgA8XJwAAPklWovVgYjMg8GOiKgRGbYTs8Q+sQa+KmcAQHZBqcXqQETmwWBHRNSILLmdmIGfR2Wwy2GwI7I5DHZERI3IGrpi1R4uABjsiGwRgx0RUSO6auE17ADA71ZXbI6GwY7I1jDYERE1IktuJ2ag9uAYOyJbxWBHRNSIDF2xlmyxMwS7q4VlKNdxyRMiW8JgR0TUSIQQuJxf2UqmvtUdagktXBVQyB0gBJB3K2gSkW1gsCMiaiR/3ChBYVkFnOQytPF2s1g9HBxk8FFVthjmFJRYrB5EZHoMdkREjeRUzk0AQIhPMzjJLfv168dxdkQ2icGOiKiRZOZoAACh6mYWrgmXPCGyVQx2RESN5OStFjtrCHZcpJjINjHYERE1ksxbwa6DFQQ7aVsxrmVHZFMY7IiIGkFpuQ6/XykEAHT0U1m4NmyxI7JVDHZERI3gbF4h9AJo7uoEHwtuJ2agZrAjskkMdkREjeBUtW5YmUxm4dpUraOXqymFXi8sXBsiMhUGOyKiRnAq2zAj1vLdsADQspkSDjKgQi9wtYiLFBPZCgY7IqJGkJlrPTNiAcBJ7oCWzQyLFLM7lshWMNgRETWCk9m3gp0VTJww4Fp2RLaHwY6IyMyuFpbhamEZZDLgPl93S1dHojZsK8YlT4hsBoMdEZGZGdavC2rhCleFo4VrU8XvVosdtxUjsh0MdkREZnbSyiZOGBiWPMllsCOyGQx2RERmdsqKdpyozrDkCVvsiGwHgx0RkZkZumI7+llZsDMsUswxdkQ2g8GOiMiMdHqB09JSJ9bVFWvYViy7oARCcJFiIlvAYEdEZEYXrhWhrEIPFyc5Alu4Wro6RnxvdcWWluuhKamwcG2IyBQY7IiIzOjUrfXr7lM3g4OD5bcSq87ZSY7mrk4AgGxNiYVrQ0SmwGBHRGRGp3JuzYj1ta7xdQZqLnlCZFMY7IiIzMgwIzbUyiZOGBjG2XH3CSLbwGBHRGRGUoudlU2cMFAz2BHZFAY7IiIzKSyrwKXrlWPXQq1sDTsDw1p2DHZEtoHBjojITAzr1/mqlGjuprBwbWpnaLHL5lp2RDaBwY6IyEwM3bAdrLQbFqgaY8dtxYhsA4MdEZGZSDtOWGk3LFB9WzEud0JkCxjsiIjMxLCGnbXOiAWqumI1pRUoKuMixURNHYMdEZEZCCGqumJ9rbcrtpmzE9yVjgC4ZyyRLWCwIyIyg+yCUmhKK+DoIEM7HzdLV6deao6zI7IZDHZERGZgaK1r29INSke5hWtTv6pxdgx2RE0dgx0RkRn8N/UiACCstadlK3IHpEWK2RVL1OQx2BERmVjKyVzsyrwCJ7kMUwe0s3R1bsuw5AlnxhI1fQx2REQmVFquw2tbTwAAxvcNRtuW7hau0e35SrtPlFm4JkR0rxjsiIhM6JN953HxWjF8minx4qPtLV2dO+IndcWyxY6oqWOwIyIykeyCEqzYeRYAMHdoqLSMiLWTxthx8gRRk8dgR0RkIgnfnUJJuQ69gppjRPdWlq7OHfPzcAEAXC3UQluht3BtiOheMNgREZlA2u/X8M3hy5DJgMWPd4ZMJrN0le5Yc1cnKBwr/+8glzNjiZo0BjsiontUodNj4TfHAQDP9glEl1YeFq5Rw8hkMmktOy55QtS0MdgREd2jdQezcCrnJlTOjnh1UAdLV+euqD24SDGRLWgaI3uJiKxQQUk5/vNDJv57oHIx4lcGdUALN4WFa3V3DC123FaMqGljsCMiaiAhBL7OuIx/bTuJq4WVa7892bMVYsIDLVyzu+fHFjsim8BgR0TUAGdyb2L+18dw4PfrACr3gn39iS54KMTbwjW7N2quZUdkExjsiIhuQ1NajuTjudh65DJ+OnMVFXoBZycHvPhoezz/cDCUjnJLV/GeGVrsTlzWQFNaDpWzk4VrRER3w6onT6xcuRJt2rSBs7MzwsPDcfDgQUtXiYjsgBAC1wrL8HXGn5j42S/o/fqPeGXjYezKvIIKvUBkR18kT++PuEdCbCLUAUDPwOZwVchx4Vox/rbqZ1y6XmzpKhHRXZAJIYSlK1GbL774AmPHjsXq1asRHh6OZcuWYePGjcjMzISPj89t36/RaODh4YGCggKoVKpGqDERNSUlWh0uXCvC+atFuHCtCH/cKMGfN0rwZ37l/5aU64zKh/i4Y1iYH4aF+SPEx/r3f70bx/4swIRPDyFXUwZvdwU+HNsbPQObW7paRDbF3PnEaoNdeHg47r//fqxYsQIAoNfrERAQgBdffBFz5sy57fsZ7IiorEKHS9eLcf5qMc5fLcT5q8W4cCvI3ckkgbbebhja1Q/Duvmhg2+zJrXo8N3KLijBhKRfcCJbA4WjA97+ezcM7+Zv6WoR2Qxz5xOrHGOn1WqRnp6OuXPnSsccHBwQGRmJ1NTUBl3r+J8FOPzbFTjJHaBwdIBCLoOT3AFOcgc4yGQo1+mh1ekr/7dCjwqdgNxBBqfblC2v0KNcJ+DgIKtxXUcHGf76/S8EUFKuQ5FWh+KyChSVVaBIq0OFTg8XhSPclXK4KhzhppTDTekIR4d7+z8QvYB0T+U6UVlnXeVWQUpHB6muTo4OcKqjvjohbr1fD61OoLxCjwp97dsNOcgqn4PhugrHyudQ+bmiWl300NfxbwnpvXIHODlWPk+5rPa6aXXG92XYBsnw2ZXXkt36u7vzZ1ai1aFYW4HCssr/LSrTQavT1Vrexcn478xN4Qilo0ON+tZFp0etv1MCNZ+PDLLK+7l1b0pHBzg6OEAvjJ+tViegq+PvqDZ/va7i1t+B4bq3+7sXAiir0KNIe+t3ukyHorIKlJTrUHarTpX3Vnk9uYOs8lkp5XBTOMJN6QgXhbzG35EQgF6IW9eo/Hxttd/hms9SSHU1PIsSrQ7ZBSXQ1/NPV09XJ7TxckMbL1cEtnBFq+YuaOVZ+b9+Hs5wdrKNbtaG8PNwwcbJEXh5/W/48WQeXvzfbziZrUFY66a16DKRNRrQ4fY9jvfKKoPd1atXodPp4Ovra3Tc19cXp06dqvU9ZWVlKCsrk14XFBQAAPafvIR3f/rTfJUlIqvmppQjsIUrglq4oo2XGwK9XBHo5YqgFm5oXueaczpoS4qgteMJom8+cR/echH4b+pFrNhx1NLVIbIJu17tD6XQAqgcy2sOVhns7kZCQgIWL15c4/jMUQ9ZoDZEZE1q/+cgEVHjCllW9edr167Bw8P0LeFWGey8vb0hl8uRm5trdDw3NxdqtbrW98ydOxczZsyQXufn5yMoKAhZWVlmeXBNkUajQUBAAC5dusRxh7fwmdTEZ1ITn0lNfCY18ZnUxGdSU0FBAQIDA9GiRQuzXN8qg51CoUCvXr2QkpKCESNGAKicPJGSkoL4+Pha36NUKqFUKmsc9/Dw4C/TX6hUKj6Tv+AzqYnPpCY+k5r4TGriM6mJz6QmBwfzrDhnlcEOAGbMmIHY2Fj07t0bffr0wbJly1BUVITnnnvO0lUjIiIiskpWG+yefvppXLlyBQsWLEBOTg66d++O7du315hQQURERESVrDbYAUB8fHydXa+3o1QqsXDhwlq7Z+0Vn0lNfCY18ZnUxGdSE59JTXwmNfGZ1GTuZ2K1CxQTERERUcNY9V6xRERERHTnGOyIiIiIbASDHREREZGNsMlgt3LlSrRp0wbOzs4IDw/HwYMHLV0ls9m7dy+GDx8Of39/yGQybNmyxei8EAILFiyAn58fXFxcEBkZiTNnzhiVuX79OmJiYqBSqeDp6YkJEyagsLCwEe/CtBISEnD//fejWbNm8PHxwYgRI5CZmWlUprS0FHFxcfDy8oK7uztGjRpVY0HsrKwsREdHw9XVFT4+Ppg5cyYqKioa81ZMZtWqVQgLC5PWkoqIiMD3338vnbe35/FX//d//weZTIZp06ZJx+zxmSxatAgymczoJzQ0VDpvj88EAP7880+MHj0aXl5ecHFxQdeuXfHLL79I5+3te7ZNmzY1fk9kMhni4uIA2OfviU6nw/z58xEcHAwXFxe0a9cOr7/+utG2YY32eyJszPr164VCoRBr1qwRx48fFxMnThSenp4iNzfX0lUzi++++078v//3/8SmTZsEALF582aj8//3f/8nPDw8xJYtW8Thw4fF448/LoKDg0VJSYlUZvDgwaJbt27iwIED4qeffhIhISHi2WefbeQ7MZ2oqCiRmJgojh07JjIyMsTQoUNFYGCgKCwslMpMnjxZBAQEiJSUFPHLL7+IBx54QDz44IPS+YqKCtGlSxcRGRkpfvvtN/Hdd98Jb29vMXfuXEvc0j375ptvxLZt28Tp06dFZmam+Oc//ymcnJzEsWPHhBD29zyqO3jwoGjTpo0ICwsTL7/8snTcHp/JwoULRefOnUV2drb0c+XKFem8PT6T69evi6CgIDFu3DiRlpYmfv/9d7Fjxw5x9uxZqYy9fc/m5eUZ/Y4kJycLAGLXrl1CCPv8Pfn3v/8tvLy8xNatW8X58+fFxo0bhbu7u1i+fLlUprF+T2wu2PXp00fExcVJr3U6nfD39xcJCQkWrFXj+Guw0+v1Qq1WizfffFM6lp+fL5RKpfjf//4nhBDixIkTAoA4dOiQVOb7778XMplM/Pnnn41Wd3PKy8sTAMSePXuEEJXPwMnJSWzcuFEqc/LkSQFApKamCiEqA7ODg4PIycmRyqxatUqoVCpRVlbWuDdgJs2bNxcff/yxXT+Pmzdvivbt24vk5GTRv39/KdjZ6zNZuHCh6NatW63n7PWZzJ49W/Tt27fO8/yeFeLll18W7dq1E3q93m5/T6Kjo8X48eONjj355JMiJiZGCNG4vyc21RWr1WqRnp6OyMhI6ZiDgwMiIyORmppqwZpZxvnz55GTk2P0PDw8PBAeHi49j9TUVHh6eqJ3795SmcjISDg4OCAtLa3R62wOBQUFACDty5eeno7y8nKj5xIaGorAwECj59K1a1ejBbGjoqKg0Whw/PjxRqy96el0Oqxfvx5FRUWIiIiw6+cRFxeH6Ohoo3sH7Pt35MyZM/D390fbtm0RExODrKwsAPb7TL755hv07t0bf//73+Hj44MePXrgo48+ks7b+/esVqvF559/jvHjx0Mmk9nt78mDDz6IlJQUnD59GgBw+PBh7Nu3D0OGDAHQuL8nVr1AcUNdvXoVOp2uxu4Uvr6+OHXqlIVqZTk5OTkAUOvzMJzLycmBj4+P0XlHR0e0aNFCKtOU6fV6TJs2DQ899BC6dOkCoPKeFQoFPD09jcr+9bnU9twM55qio0ePIiIiAqWlpXB3d8fmzZvRqVMnZGRk2OXzWL9+PX799VccOnSoxjl7/R0JDw9HUlISOnTogOzsbCxevBgPP/wwjh07ZrfP5Pfff8eqVaswY8YM/POf/8ShQ4fw0ksvQaFQIDY21u6/Z7ds2YL8/HyMGzcOgP3+tzNnzhxoNBqEhoZCLpdDp9Ph3//+N2JiYgA07v8f21SwI/qruLg4HDt2DPv27bN0VSyuQ4cOyMjIQEFBAb788kvExsZiz549lq6WRVy6dAkvv/wykpOT4ezsbOnqWA1D6wIAhIWFITw8HEFBQdiwYQNcXFwsWDPL0ev16N27N5YsWQIA6NGjB44dO4bVq1cjNjbWwrWzvE8++QRDhgyBv7+/patiURs2bMDatWuxbt06dO7cGRkZGZg2bRr8/f0b/ffEprpivb29IZfLa8y+yc3NhVqttlCtLMdwz/U9D7Vajby8PKPzFRUVuH79epN/ZvHx8di6dSt27dqF1q1bS8fVajW0Wi3y8/ONyv/1udT23AznmiKFQoGQkBD06tULCQkJ6NatG5YvX26XzyM9PR15eXno2bMnHB0d4ejoiD179uDdd9+Fo6MjfH197e6Z1MbT0xP33Xcfzp49a5e/JwDg5+eHTp06GR3r2LGj1EVtz9+zFy9exI8//ojnn39eOmavvyczZ87EnDlz8Mwzz6Br164YM2YMpk+fjoSEBACN+3tiU8FOoVCgV69eSElJkY7p9XqkpKQgIiLCgjWzjODgYKjVaqPnodFokJaWJj2PiIgI5OfnIz09XSqzc+dO6PV6hIeHN3qdTUEIgfj4eGzevBk7d+5EcHCw0flevXrBycnJ6LlkZmYiKyvL6LkcPXrU6D+y5ORkqFSqGl/yTZVer0dZWZldPo+BAwfi6NGjyMjIkH569+6NmJgY6c/29kxqU1hYiHPnzsHPz88uf08A4KGHHqqxXNLp06cRFBQEwH6/ZwEgMTERPj4+iI6Olo7Z6+9JcXExHByMI5VcLoderwfQyL8n9zAJxCqtX79eKJVKkZSUJE6cOCEmTZokPD09jWbf2JKbN2+K3377Tfz2228CgPjPf/4jfvvtN3Hx4kUhROX0ak9PT/H111+LI0eOiCeeeKLW6dU9evQQaWlpYt++faJ9+/ZNdhq+EEJMmTJFeHh4iN27dxtNyS8uLpbKTJ48WQQGBoqdO3eKX375RURERIiIiAjpvGE6/qBBg0RGRobYvn27aNmyZZOdjj9nzhyxZ88ecf78eXHkyBExZ84cIZPJxA8//CCEsL/nUZvqs2KFsM9n8sorr4jdu3eL8+fPi/3794vIyEjh7e0t8vLyhBD2+UwOHjwoHB0dxb///W9x5swZsXbtWuHq6io+//xzqYw9fs/qdDoRGBgoZs+eXeOcPf6exMbGilatWknLnWzatEl4e3uLWbNmSWUa6/fE5oKdEEK89957IjAwUCgUCtGnTx9x4MABS1fJbHbt2iUA1PiJjY0VQlROsZ4/f77w9fUVSqVSDBw4UGRmZhpd49q1a+LZZ58V7u7uQqVSieeee07cvHnTAndjGrU9DwAiMTFRKlNSUiKmTp0qmjdvLlxdXcXIkSNFdna20XUuXLgghgwZIlxcXIS3t7d45ZVXRHl5eSPfjWmMHz9eBAUFCYVCIVq2bCkGDhwohToh7O951Oavwc4en8nTTz8t/Pz8hEKhEK1atRJPP/200Xpt9vhMhBDi22+/FV26dBFKpVKEhoaKDz/80Oi8PX7P7tixQwCocZ9C2OfviUajES+//LIIDAwUzs7Oom3btuL//b//Z7R8S2P9nsiEqLYsMhERERE1WTY1xo6IiIjInjHYEREREdkIBjsiIiIiG8FgR0RERGQjGOyIiIiIbASDHREREZGNYLAjIiIishEMdkREREQ2gsGOiOg22rRpg2XLllm6GkREt8VgR0R0G4cOHcKkSZPuqKwQAgsWLICfnx9cXFwQGRmJM2fOmLmGRESVGOyIiG6jZcuWcHV1vaOyS5cuxbvvvovVq1cjLS0Nbm5uiIqKQmlpqZlrSUTEYEdEdmD79u3o27cvPD094eXlhWHDhuHcuXMAgM8++wzu7u5GrWpTp05FaGgoiouLARh3xQohsGjRIgQGBkKpVMLf3x8vvfSSdG7ZsmWYN28ennjiCYSFheGzzz7D5cuXsWXLlka9ZyKyTwx2RGTzioqKMGPGDPzyyy9ISUmBg4MDRo4cCb1ej7Fjx2Lo0KGIiYlBRUUFtm3bho8//hhr166ttZXuq6++wjvvvIMPPvgAZ86cwZYtW9C1a1cAwPnz55GTk4PIyEipvIeHB8LDw5Gamtpo90tE9svR0hUgIjK3UaNGGb1es2YNWrZsiRMnTqBLly744IMPEBYWhpdeegmbNm3CokWL0KtXr1qvlZWVBbVajcjISDg5OSEwMBB9+vQBAOTk5AAAfH19jd7j6+srnSMiMie22BGRzTtz5gyeffZZtG3bFiqVCm3atAFQGdIAoHnz5vjkk0+watUqtGvXDnPmzKnzWn//+99RUlKCtm3bYuLEidi8eTMqKioa4zaIiG6LwY6IbN7w4cNx/fp1fPTRR0hLS0NaWhoAQKvVSmX27t0LuVyO7OxsFBUV1XmtgIAAZGZm4v3334eLiwumTp2Kfv36oby8HGq1GgCQm5tr9J7c3FzpHBGROTHYEZFNu3btGjIzMzFv3jwMHDgQHTt2xI0bN4zK/Pzzz3jjjTfw7bffwt3dHfHx8fVe08XFBcOHD8e7776L3bt3IzU1FUePHkVwcDDUajVSUlKkshqNBmlpaYiIiDDL/RERVccxdkRk05o3bw4vLy98+OGH8PPzQ1ZWllFX682bNzFmzBi89NJLGDJkCFq3bo37778fw4cPx9/+9rca10tKSoJOp0N4eDhcXV3x+eefw8XFBUFBQZDJZJg2bRr+9a9/oX379ggODsb8+fPh7++PESNGNOJdE5G9YosdEdk0BwcHrF+/Hunp6ejSpQumT5+ON998Uzr/8ssvw83NDUuWLAEAdO3aFUuWLMELL7yAP//8s8b1PD098dFHH+Ghhx5CWFgYfvzxR3z77bfw8vICAMyaNQsvvvgiJk2ahPvvvx+FhYXYvn07nJ2dG+eGiciuyYQQwtKVICIiIqJ7xxY7IiIiIhvBYEdERERkIxjsiIiIiGwEgx0RERGRjWCwIyIiIrIRDHZERERENoLBjoiIiMhGMNgRERER2QgGOyIiIiIbwWBHREREZCMY7IiIiIhsBIMdERERkY34/zDgNc1eSrmAAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h.plot(\"line\", show_stats=True);" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Histogram1D('Edep in absorber', bins=(100,), total=10000.0, dtype=float64)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h2 = load_csv(\"../tests/data/geant-h2.csv\")\n", "h" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHWCAYAAACWppaPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8g+/7EAAAACXBIWXMAAA9hAAAPYQGoP6dpAABRjklEQVR4nO3de1hU1f4/8DcoDCDOAAqMJCBe0VAxMxwvZUriJZMyU/PnpbwUgamYqSfvlZhapqaYHYNuHtNKS0UMMfWkSMoR88pRIzF1sFRAMAaB9fvDL/s4ctvAMNf363nm0dlrzd5rsWX2x8/aa207IYQAERERkQ2zN3UDiIiIiEyNARERERHZPAZEREREZPMYEBEREZHNY0BERERENo8BEREREdk8BkRERERk8xgQERERkc1jQEREREQ2jwERUT0YP348WrRoobctPz8fEydOhFqthp2dHaZNm1avbdi/fz/s7OzwzTff1OtxaiI+Ph52dnY4duyYqZtCRKSHARER/nehLns5OTnBx8cHYWFhWL16NW7fvl3nYyxZsgTx8fGIiIjAF198gTFjxuDw4cNYuHAhcnJyarSv/fv347nnnoNarYajoyO8vLwwZMgQfPfdd3Vup7XLz8/HggULMGDAAHh4eMDOzg7x8fGyPiuEQO/eveHp6YkbN26UK3/11Vfh4OCA9PR0wzZaps2bN8POzg4ff/xxheURERFwcHDAiRMnjNwyIgsgiEjExcUJAGLx4sXiiy++EJ9++qlYsmSJ6N+/v7CzsxP+/v7ixIkTsvdXVFQkCgsL9baFhISInj176m1bvny5ACAyMzNl73v+/PkCgGjTpo2YP3++2Lhxo1i2bJno06ePACC++uorIYQQP/30kwAgtm7dKnvf9a3s53z06FGTtSEzM1MAEH5+ftLPLC4uTvbnT58+LRwcHMT48eP1th8+fFjY2dmJGTNmGLjFNTNgwADh5uYmtFqt3vbU1FRhb28vZs6caaKWEZk3BkREouoLdXJysnB2dhb+/v7izp07Ve4nPz+/0rKAgAAxePBgvW01DYi2bt0qAIjnn39eFBUVlStPTEwUO3bsEELYdkBU1XkoLCwU165dE0IIcfTo0RoHREII8Y9//EMAEPv37xdC3AuAg4KChJ+fX5XHNobMzEzh4uIiRo0aJW0rLi4WwcHBokWLFqKgoMCErSMyXxwyI6pG3759MW/ePFy6dAlffvmltH38+PFwdXXFxYsXMWjQIDRu3BijR4+WysruISq7lyczMxO7du2ShuXGjx+PmTNnAgACAgKk7b///nulbZk3bx48PDzw6aefwsHBoVx5WFgYnn76ab1tpaWlePfdd9G8eXM4OTmhX79+uHDhgl6df//73xg+fDj8/PygUCjg6+uL6dOn4++//9arV9bnrKwsPP3003B1dcVDDz2EtWvXAgBOnjyJvn37olGjRvD398emTZsq7MedO3fwyiuvoEmTJlAqlRg7dixu3bpVrt7u3bvRu3dvNGrUCI0bN8bgwYNx+vTpCttU0XmoiEKhgFqtrrRcjnnz5qFVq1Z45ZVXUFRUhPfffx+nTp3CRx99hEaNGlX52aCgIDz55JPltpeWluKhhx7C888/L23bvHkzunbtisaNG0OpVKJjx45YtWpVlftv0aIFFi5ciH/9619ISkoCAKxevRrp6emIjY2Fi4tLLXpMZP0YEBHJMGbMGADAjz/+qLe9uLgYYWFh8PLywooVKzBs2LByn23fvj2++OILNG3aFMHBwfjiiy/wxRdfYPr06Rg1ahQAYOXKldJ2T0/PCttw/vx5nDt3DuHh4WjcuLHsti9duhTbtm3DG2+8gTlz5uDIkSPlAoatW7fizp07iIiIwJo1axAWFoY1a9Zg7Nix5fZXUlKCgQMHwtfXF8uWLUOLFi0QFRWF+Ph4DBgwAI8++ijee+89NG7cGGPHjkVmZma5fURFReHs2bNYuHAhxo4di6+++grh4eEQQkh1vvjiCwwePBiurq547733MG/ePJw5cwa9evUqFzTKOQ+G5OTkhHXr1iEjIwOvvfYaFi9ejGeffRZDhgyp9rMjRozAwYMHodVq9bb//PPPuHr1KkaOHAkASEpKwqhRo+Du7o733nsPS5cuRZ8+fXDo0KFqjzF9+nR07twZERERuHDhAubPn4+RI0diwIABteswkS0wdYqKyBzIGcpRqVSiS5cu0vtx48YJAGL27Nnl6o4bN074+/vrbfP396/TkNn3338vAIiVK1dWW1eI/w2ZtW/fXuh0Omn7qlWrBABx8uRJaVtFQ4ExMTHCzs5OXLp0Sa9fAMSSJUukbbdu3RLOzs7Czs5ObN68Wdp+7tw5AUAsWLBA2lb2c+7atavekN+yZcsEAPH9998LIYS4ffu2cHNzE5MmTdJrk1arFSqVSm97VeehOrUdMiszatQoAUA0btxYXL58WdZnMjIyBACxZs0ave2vvfaacHV1lc7F1KlThVKpFMXFxbVqW9k9Qx4eHhXeU0RE+pghIpLJ1dW1wtlmERERRjl+Xl4eANQoOwQAL730EhwdHaX3vXv3BgD89ttv0jZnZ2fp7wUFBfjrr7/Qo0cPCCFw/PjxcvucOHGi9Hc3Nze0a9cOjRo1wgsvvCBtb9euHdzc3PSOU2by5Ml6Q34RERFo2LAhEhISANzLjuTk5GDUqFH466+/pFeDBg0QEhKCn376qdw+jXUe7te0aVMAQIcOHdC8eXNZn2nbti2Cg4Px9ddfS9tKSkrwzTffYMiQIdK5cHNzQ0FBgTTsVVOPPfYYXn31Vdy8eRMxMTHw9vau1X6IbAUDIiKZ8vPzywUjDRs2lH0hrCulUgkANV4CwM/PT++9u7s7AOjds5OVlYXx48fDw8MDrq6u8PT0xBNPPAEAyM3N1fu8k5NTuWE9lUqF5s2bw87Ortz2iu4NatOmjd57V1dXNGvWTBoKO3/+PIB79295enrqvX788Udcv35d7/PGPA9ljh07hrVr1yIoKAipqal695dVZ8SIETh06BCuXLkC4N59ZtevX8eIESOkOq+99hratm2LgQMHonnz5nj55ZeRmJhYozZ269YNAPDoo4/W6HNEtqihqRtAZAn++OMP5ObmonXr1nrbFQoF7O2N8/+KwMBAAPduXK6JBg0aVLhd/N/9OiUlJXjqqadw8+ZNzJo1C4GBgWjUqBGuXLmC8ePHo7S0VNb+qjtOTZQd84svvqjwBuiGDfW/uox5HoB7P7PJkyfDx8cHhw4dQv/+/TFjxgw8/fTTcHNzq/bzI0aMwJw5c7B161ZMmzYNW7ZsgUql0rvHx8vLC+np6dizZw92796N3bt3Iy4uDmPHjsVnn31Wj70jsk3MEBHJ8MUXXwC4N4vLkB7MqFSlbdu2aNeuHb7//nvk5+cbrA0nT57Ef//7X7z//vuYNWsWhg4ditDQUPj4+BjsGA8qywCVyc/Px7Vr16SZea1atQJwLygIDQ0t9+rTp0+9tU2O1atX4/jx41izZg2USiXWr1+PGzduYPbs2bI+HxAQgMceewxff/01iouL8d133yE8PBwKhUKvnqOjI4YMGYJ169bh4sWLeOWVV/D555+XmyVIRHXHgIioGvv27cPbb7+NgICAKqdz10bZFG25K1UvWrQIN27cwMSJE1FcXFyu/Mcff8TOnTtr1IayzM79mRwhRLXTu+tiw4YNuHv3rvQ+NjYWxcXFGDhwIIB7gadSqcSSJUv06pX5888/661t1bl8+TLmz5+PZ555BuHh4QCA4OBgvP766/jkk0+Qmpoqaz8jRozAkSNH8Omnn+Kvv/7SGy4DUG4lbHt7e3Tq1AkAoNPp6t4RItLDITOi++zevRvnzp1DcXExsrOzsW/fPiQlJcHf3x8//PADnJycDHq8rl27AgDeeustjBw5Eg4ODhgyZEila9mMGDECJ0+exLvvvovjx49j1KhR8Pf3x40bN5CYmIjk5ORK1/6pTGBgIFq1aoU33ngDV65cgVKpxLffflvhvT+GUlRUhH79+uGFF15ARkYG1q1bh169euGZZ54BcO9+qdjYWIwZMwaPPPIIRo4cCU9PT2RlZWHXrl3o2bMnPvroo1of/6OPPkJOTg6uXr0KANixYwf++OMPAMCUKVOgUqkq/eyUKVMghMCaNWv0ti9atAhbtmzBq6++imPHjlU6hFjmhRdewBtvvIE33ngDHh4eCA0N1SufOHEibt68ib59+6J58+a4dOkS1qxZg+DgYLRv37423SaiqphyihuRuSibDl72cnR0FGq1Wjz11FNi1apVIi8vr9xnxo0bJxo1alTh/uROuxdCiLfffls89NBDwt7eXvYU/OTkZDF06FDh5eUlGjZsKDw9PcWQIUOkaetCVL5SddmjK+6fan7mzBkRGhoqXF1dRdOmTcWkSZPEiRMnytWrrM9PPPGEePjhh8ttf7DPZT/nAwcOiMmTJwt3d3fh6uoqRo8eLW7cuFHu8z/99JMICwsTKpVKODk5iVatWonx48eLY8eOVdumqvj7++ud7/tfVf38t23bJgCIFStWVFj+zTffCADigw8+kNWOnj17CgBi4sSJFe6rf//+wsvLSzg6Ogo/Pz/xyiuvSKtsy2EOj0ohshR2QtTijkciIiIiK8J7iIiIiMjmMSAiIiIim8eAiIiIiGyexQREsbGx6NSpE5RKJZRKJTQaDXbv3i2VFxYWIjIyEk2aNIGrqyuGDRuG7OxsvX1kZWVh8ODBcHFxgZeXF2bOnFnh1GUiIiIyDy1atICdnV25V2RkJAB51385LCYgat68OZYuXYq0tDQcO3YMffv2xdChQ3H69GkA957uvGPHDmzduhUHDhzA1atX8dxzz0mfLykpweDBg1FUVITDhw/js88+Q3x8PObPn2+qLhEREVE1jh49imvXrkmvsuf7DR8+HED113/ZTD3NrS7c3d3FP//5T5GTkyMcHBz0phefPXtWABApKSlCCCESEhKEvb293hOfY2NjhVKp1HsSOBEREZmvqVOnilatWonS0lJZ13+5LHJhxpKSEmzduhUFBQXQaDRIS0vD3bt39RY2CwwMhJ+fH1JSUtC9e3ekpKSgY8eOek98DgsLQ0REBE6fPo0uXbpUeCydTqe3KmxpaSlu3ryJJk2a1OixC0RERIYkhMDt27fh4+NjlGf5FRYWoqioyCD7EkKUu4YqFIpyj695UFFREb788ktER0fDzs5O1vVfLosKiE6ePAmNRoPCwkK4urpi27Zt6NChA9LT0+Ho6FjuoYre3t7QarUAAK1WqxcMlZWXlVUmJiYGixYtMmxHiIiIDOTy5cto3rx5vR6jsLAQzs7OBtufq6truWcyLliwAAsXLqzyc9u3b0dOTg7Gjx8P4N71u7rrv1wWFRC1a9cO6enpyM3NxTfffINx48bhwIED9XrMOXPmIDo6Wnqfm5sLPz8/XL58GUqlsl6PTURE9WPPnj3o379/peU//vijwR/mbGh5eXnw9fVF48aN6/1YhsoMlcnPzy93Ha0uOwQAGzduxMCBA+vl4dMWFRA5OjqidevWAO49A+ro0aNYtWoVRowYgaKiIuTk5OhFidnZ2VCr1QAAtVqNX375RW9/ZXehl9WpSGUpvLLZbkREZHlcXFyq/A6vrtycGPv2jboeT/zfAzJqeh29dOkS9u7di++++07aplarq73+y2Uxs8wqUlpaCp1Oh65du8LBwQHJyclSWUZGBrKysqDRaAAAGo0GJ0+exPXr16U6SUlJUCqV6NChg9HbTkREZGkqmv5em1dtxMXFwcvLC4MHD5a2ybn+y2UxGaI5c+Zg4MCB8PPzw+3bt7Fp0ybs378fe/bsgUqlwoQJExAdHQ0PDw8olUpMmTIFGo1GuqGqf//+6NChA8aMGYNly5ZBq9Vi7ty5iIyMlJWmIyIiItMoLS1FXFwcxo0bh4YN/xe6yLn+y2UxAdH169cxduxYXLt2DSqVCp06dcKePXvw1FNPAQBWrlwJe3t7DBs2DDqdDmFhYVi3bp30+QYNGmDnzp2IiIiARqNBo0aNMG7cOCxevNhUXSIiIrIodcnw3E/U8Lnye/fuRVZWFl5++eVyZdVd/+Xi0+5rKC8vDyqVCrm5uRYzvkxERPp27dqFQYMGVVqekJCgNzRjjox5PSo7VsOGDQ1yD1FxcbHZXUct+h4iIiIiIkOwmCEzIiIiQ1EoFEhISKiynMoz1JCZOWJARERENuf+lY1JPmsOiDhkRkRERDaPGSIiIiKSxZozRAyIiIiISBZrDog4ZEZEREQ2jxkiIiIiksWaM0QMiIiIiEgWe3t7gz3c1dxwyIyIiIhsHjNEREREJAuHzIiIiMjmWXNAxCEzIiIisnnMEBEREZEs1pwhYkBEREREslhzQMQhMyIiIrJ5zBARERGRLNacIWJARERERLLY2dnB3r5ug0ulpaUGao1hcciMiIiIbB4zRERERCSLIYbMzHXIjQERERERyWLNARGHzIiIiMjmMUNEREREslhzhogBEREREclizQERh8yIiIjI5jFDRERERLJYc4aIARERERHJYm9vX+eFGc2VdfaKiIiIqAaYISIiIiJZOGRGRERENs+aAyIOmREREZHNY4aIiIiIZLHmDBEDIiIiIpLFmgMiDpkRERGRzWOGiIiIiGSx5gwRAyIiIiKShQszEhEREVkxZoiIiIhIFg6ZERERkc2z5oCIQ2ZERERk85ghIiIiIlmsOUPEgIiIiIhkM9eApq44ZEZEREQ2jxkiIiIiksUQ6xAJIQzUGsNihoiIiIhkKbuHqK6vmrpy5Qr+3//7f2jSpAmcnZ3RsWNHHDt2TCoXQmD+/Plo1qwZnJ2dERoaivPnz9foGAyIiIiIyGzdunULPXv2hIODA3bv3o0zZ87g/fffh7u7u1Rn2bJlWL16NdavX4/U1FQ0atQIYWFhKCwslH0cDpkRERGRLKaYZfbee+/B19cXcXFx0raAgADp70IIfPjhh5g7dy6GDh0KAPj888/h7e2N7du3Y+TIkbKOwwwRERERyVJ2D1FdXwCQl5en99LpdBUe84cffsCjjz6K4cOHw8vLC126dMEnn3wilWdmZkKr1SI0NFTaplKpEBISgpSUFPl9q+XPhIiIiKjWfH19oVKppFdMTEyF9X777TfExsaiTZs22LNnDyIiIvD666/js88+AwBotVoAgLe3t97nvL29pTI5OGRGREREshhyyOzy5ctQKpXSdoVCUWH90tJSPProo1iyZAkAoEuXLjh16hTWr1+PcePG1akt92OGiIiIiGQx5JCZUqnUe1UWEDVr1gwdOnTQ29a+fXtkZWUBANRqNQAgOztbr052drZUJqtvsmsSERERGVnPnj2RkZGht+2///0v/P39Ady7wVqtViM5OVkqz8vLQ2pqKjQajezjcMiMiIiIZDHFLLPp06ejR48eWLJkCV544QX88ssv2LBhAzZs2CDtb9q0aXjnnXfQpk0bBAQEYN68efDx8UF4eLjs41hMhigmJgbdunVD48aN4eXlhfDw8HIRY2FhISIjI9GkSRO4urpi2LBh5VJoWVlZGDx4MFxcXODl5YWZM2eiuLjYmF0hIiKySKZYmLFbt27Ytm0b/vWvfyEoKAhvv/02PvzwQ4wePVqq8+abb2LKlCmYPHkyunXrhvz8fCQmJsLJyUl+34S5rqH9gAEDBmDkyJHo1q0biouL8Y9//AOnTp3CmTNn0KhRIwBAREQEdu3ahfj4eKhUKkRFRcHe3h6HDh0CAJSUlCA4OBhqtRrLly/HtWvXMHbsWEyaNEm6Was6eXl5UKlUyM3N1bsZjIiIyJiMeT0qO1bHjh3RoEGDOu2rpKQEJ0+eNLvrqMUERA/6888/4eXlhQMHDuDxxx9Hbm4uPD09sWnTJjz//PMAgHPnzqF9+/ZISUlB9+7dsXv3bjz99NO4evWqND1v/fr1mDVrFv788084OjpWe1wGRLZr7969la6TUUahUOithUFEVF9MERB17tzZIAHRiRMnzO46arH3EOXm5gIAPDw8AABpaWm4e/eu3sUoMDAQfn5+UkCUkpKCjh076q1VEBYWhoiICJw+fRpdunQpdxydTqd3EczLy6uvLpGZ0+l0GDRoUJV1EhISjNQaIiLjM8U9RMZiMfcQ3a+0tBTTpk1Dz549ERQUBODewkyOjo5wc3PTq3v/wkxarbbChZvKyioSExOjt3CUr6+vgXtDREREpmaRAVFkZCROnTqFzZs31/ux5syZg9zcXOl1+fLlej8mERGROTLkOkTmxuKGzKKiorBz504cPHgQzZs3l7ar1WoUFRUhJydHL0t0/8JMarUav/zyi97+ymahVbZ4k0KhqHSxKCIiIlvCITMzIIRAVFQUtm3bhn379uk96RYAunbtCgcHB72FmTIyMpCVlSUtzKTRaHDy5Elcv35dqpOUlASlUlluFUwiIiKyHRaTIYqMjMSmTZvw/fffo3HjxtI9PyqVCs7OzlCpVJgwYQKio6Ph4eEBpVKJKVOmQKPRoHv37gCA/v37o0OHDhgzZgyWLVsGrVaLuXPnIjIyklkgIiKiatjZ2dV5yKu0tNRArTEsiwmIYmNjAQB9+vTR2x4XF4fx48cDAFauXAl7e3sMGzYMOp0OYWFhWLdunVS3QYMG2LlzJyIiIqDRaNCoUSOMGzcOixcvNlY3iIiILJY1D5lZTEAkZ7kkJycnrF27FmvXrq20jr+/P6dGU60oFIpq/+0w00hEZJksJiAiMjUuuEhEto4ZIiIiIrJ5hpg2b67T7s2zVURERERGxAwRERERycIhMyIiIrJ5HDIjIiIismLMEBEREZEsHDIjIiIim2fNARGHzIiIiMjmMUNEREREsljzTdUMiIiIiEgWDpkRERERWTFmiIiIiEgWDpkRERGRzbPmITMGRERUK3v37oVOp6u0XKFQIDQ01IgtIiKqPQZERFQrOp0OgwYNqrQ8ISHBiK0hImNghoiIiIhsnp2dXZ3vATLXgMg872wiIiIiMiJmiIiIiEgWDpkRERGRzbPmaffm2SoiIiIiI2KGiIiIiGThkBkR0QMUCkWVU+sVCoURW0NExmDNQ2YMiIioVrjoIhFZEwZEREREJAuHzIiIiMjmWXNAZJ4DeURERERGxAwRERERyWLNGSIGRERERCSLNQdEHDIjIiIim8cMEREREclizRkiBkREREQkizUHRBwyIyIiIpvHDBERERHJYs0ZIgZEREREJIs1B0QcMiMiIiKztXDhQikQK3sFBgZK5YWFhYiMjESTJk3g6uqKYcOGITs7u8bHYYaIiIiIZDHV0+4ffvhh7N27V3rfsOH/wpfp06dj165d2Lp1K1QqFaKiovDcc8/h0KFDNToGAyIiIqJa2Lt3L3Q6XaXlCoUCoaGh9XaMO3fu1GnftWGqIbOGDRtCrVaX256bm4uNGzdi06ZN6Nu3LwAgLi4O7du3x5EjR9C9e3f5x6hxq4iIiAg6nQ6DBg2qtDwhIaFej5GXl1fn/ZvSg+1XKBRQKBQV1j1//jx8fHzg5OQEjUaDmJgY+Pn5IS0tDXfv3tULPAMDA+Hn54eUlJQaBUS8h4iIiIhke/B+npq+yvj6+kKlUkmvmJiYCo8XEhKC+Ph4JCYmIjY2FpmZmejduzdu374NrVYLR0dHuLm56X3G29sbWq22Rv1ihoiIiIhkMeSQ2eXLl6FUKqXtlWWHBg4cKP29U6dOCAkJgb+/P7Zs2QJnZ+c6teV+zBARERGR0SmVSr1XZQHRg9zc3NC2bVtcuHABarUaRUVFyMnJ0auTnZ1d4T1HVWFARERERLLUdbjMEBmm/Px8XLx4Ec2aNUPXrl3h4OCA5ORkqTwjIwNZWVnQaDQ12i+HzIiIiEgWU8wye+ONNzBkyBD4+/vj6tWrWLBgARo0aIBRo0ZBpVJhwoQJiI6OhoeHB5RKJaZMmQKNRlOjG6oBBkRERERkxv744w+MGjUKN27cgKenJ3r16oUjR47A09MTALBy5UrY29tj2LBh0Ol0CAsLw7p162p8HAZERER1ZIz1aMj8KBSKKqfWy70nprbHsJV1iDZv3lxluZOTE9auXYu1a9fWpVkMiIiI6soY69GQ+TFGkFvVMUyxDhGfZUZERERkxZghIiIiIlmsOUPEgIiIiIhkseaAiENmREREZPOYISIiIiJZrDlDxICIiKyaMabEG2P6NZE5YEBkJg4ePIjly5cjLS0N165dw7Zt2xAeHi6VCyGwYMECfPLJJ8jJyUHPnj0RGxuLNm3aSHVu3ryJKVOmYMeOHdJCTqtWrYKrq6sJekRE9c0YU+K5xhCR5bOoe4gKCgrQuXPnShdfWrZsGVavXo3169cjNTUVjRo1QlhYGAoLC6U6o0ePxunTp5GUlISdO3fi4MGDmDx5srG6QEREZLHM4Vlm9cWiMkQDBw7EwIEDKywTQuDDDz/E3LlzMXToUADA559/Dm9vb2zfvh0jR47E2bNnkZiYiKNHj+LRRx8FAKxZswaDBg3CihUr4OPjY7S+EBERWRprHjKzqAxRVTIzM6HVavVS1yqVCiEhIUhJSQEApKSkwM3NTQqGgHupbnt7e6Smpla4X51Oh7y8PL0XERERWRerCYi0Wi0AwNvbW2+7t7e3VKbVauHl5aVX3rBhQ3h4eEh1HhQTEwOVSiW9fH1966H1RERE5s+ah8ysJiCqL3PmzEFubq70unz5sqmbREREZBIMiCyAWq0GAGRnZ+ttz87OlsrUajWuX7+uV15cXIybN29KdR6kUCigVCr1XkRERGRdLOqm6qoEBARArVYjOTkZwcHBAO49CTg1NRUREREAAI1Gg5ycHKSlpaFr164AgH379qG0tBQhISGmajoR1SOuEURkONZ8U7VFBUT5+fm4cOGC9D4zMxPp6enw8PCAn58fpk2bhnfeeQdt2rRBQEAA5s2bBx8fH2mtovbt22PAgAGYNGkS1q9fj7t37yIqKgojR47kDDMiK8U1gogMhwGRmTh27BiefPJJ6X10dDQAYNy4cYiPj8ebb76JgoICTJ48GTk5OejVqxcSExPh5OQkfearr75CVFQU+vXrJy3MuHr1aqP3hYiIiMyHnRBCmLoRliQvLw8qlQq5ubm8n4iIiEzGmNejsmNNnjwZjo6OddpXUVERNmzYYHbXUYvKEBEREZHpWPOQmdXMMiMiIiKqLWaIiIiISDZzzfDUFQMiIiIikoVDZkRERERWjBkiIiIiksWaM0QMiIiIiEgWaw6IOGRGRERENo8ZIiIiIpLFmjNEDIiIiIhIFmsOiDhkRkRERDaPGSIioirs3bsXOp2uyjoKhQKhoaFGahGR6djb28Pevm65lLp+vr4wICIiqoJOp8OgQYOqrJOQkGCk1hCZFofMiIiIiKwYM0REREQkizVniBgQERERkSzWHBBxyIyIiIhsHjNEREREJIs1Z4gYEBEREZEsDIiIiGyUQqGodlq9QqEwUmuIqL4wICIiqgIXXCT6H2aIiIiIyOZZc0DEWWZERERk85ghIiIiIlmsOUPEgIiIiIhkseaAiENmREREZPMYEBEREZEsZRmiur5qa+nSpbCzs8O0adOkbYWFhYiMjESTJk3g6uqKYcOGITs7u8b75pAZEZGN2Lt3L3Q6XaXlCoWCywxQlezs7GBvX7dcSm0DoqNHj+Ljjz9Gp06d9LZPnz4du3btwtatW6FSqRAVFYXnnnsOhw4dqtH+GRAREdkInU6HQYMGVVpe3QKURKaSn5+P0aNH45NPPsE777wjbc/NzcXGjRuxadMm9O3bFwAQFxeH9u3b48iRI+jevbvsY3DIjIiIiGQx5JBZXl6e3quq7GVkZCQGDx5cLoOZlpaGu3fv6m0PDAyEn58fUlJSatQ3BkRERERkdL6+vlCpVNIrJiamwnqbN2/Gf/7znwrLtVotHB0d4ebmprfd29sbWq22Ru3hkBkRERHJYshp95cvX4ZSqZS2V/RMwMuXL2Pq1KlISkqCk5NTnY5bHWaIiIiISBZDDpkplUq9V0UBUVpaGq5fv45HHnkEDRs2RMOGDXHgwAGsXr0aDRs2hLe3N4qKipCTk6P3uezsbKjV6hr1jRkiIiIiMkv9+vXDyZMn9ba99NJLCAwMxKxZs+Dr6wsHBwckJydj2LBhAICMjAxkZWVBo9HU6FgMiIgMiNOajau6nzfAn/n9FApFlTPJKvofuinw98h8GXul6saNGyMoKEhvW6NGjdCkSRNp+4QJExAdHQ0PDw8olUpMmTIFGo2mRjPMAAZERAbFac3GVd3PG+DP/H6WEkTw98h8meOjO1auXAl7e3sMGzYMOp0OYWFhWLduXY33w4CIiIiILMb+/fv13js5OWHt2rVYu3ZtnfbLgIiIiIhkMccMkaEwICIiIiJZrDkg4rR7IiIisnnMEBEREZEs1pwhYkBEREREsjAgkik7Oxsff/wx5s+fb8jdElkMS1nnxVpU9/Muq2NqXC+pZmzp96i6fxv5+flwdXWtsOzOnTv11SybZNCASKvVYtGiRQyIyGbxgmZclvLz5npJNWMp59UQqvu3sXXr1krL8/Ly6qtZlbK3t4e9fd1uP67r5+tLjQKiX3/9tcryjIyMOjWGiIiIzBeHzP5PcHAw7OzsIIQoV1a23Vw7SkRERFSZGgVEHh4eWLZsGfr161dh+enTpzFkyBCDNIyIiIjMCzNE/6dr1664evUq/P39KyzPycmpMHtERERElo8B0f959dVXUVBQUGm5n58f4uLi6twoIiIiImOqUUD07LPPVlnu7u6OcePG1alBREREZJ6YIapGSUkJTp48CX9/f7i7uxtil0Rk5qpbPyUzMxMlJSWVlhcXF6Nhw6q/gjw9PfHiiy/Wuo2GYIg1hCxlvSQyvur+bQghKi03xTpEDIgeMG3aNHTs2BETJkxASUkJnnjiCRw+fBguLi7YuXMn+vTpY+BmEpG5qW79lJiYGCxfvrzS8qlTp+L999+v8hjvvPNOrdtnKIZYQ8iW1tWhmqnLvw1TrENkzWq1OtI333yDzp07AwB27NiBzMxMnDt3DtOnT8dbb71l0AYSERGR+SjLEtX2Za5qFRD99ddfUKvVAO79z2j48OFo27YtXn75ZZw8edKgDawva9euRYsWLeDk5ISQkBD88ssvpm4SERGRWatrMGTOQVGtAiJvb2+cOXMGJSUlSExMxFNPPQXg3nhmgwYNDNrA+vD1118jOjoaCxYswH/+8x907twZYWFhuH79uqmbRkRERCZQq4DopZdewgsvvICgoCDY2dlJY6CpqakIDAw0aAPrwwcffIBJkybhpZdeQocOHbB+/Xq4uLjg008/NXXTiIiIzJY1Z4hqdVP1woULERQUhMuXL2P48OHS7IgGDRpg9uzZBm2goRUVFSEtLQ1z5syRttnb2yM0NBQpKSkmbBkREZF54yyzCjz//PPltlnCGkR//fUXSkpK4O3trbfd29sb586dK1dfp9PpTbnlXf1ERETWR3ZAtHr1akyePBlOTk5YvXp1lXVff/31OjfMXMTExGDRokWmbgaRUclZe+fWrVvYunVrpeX29vaYNm1apeUNGjRAdHR0lccoLi6ustwYuIYQ0f/Y29vD3r5Wd9vo7cMcyQ6IVq5cidGjR8PJyQkrV66stJ6dnZ1ZB0RNmzZFgwYNkJ2drbc9Oztbmjl3vzlz5uh9aefl5cHX17fe20lkSnLW3tmyZQuGDh1aafnKlSvx+eefV1o+evRovP3221Ue44MPPqi6oUbANYSI/odDZri36mxFf7c0jo6O6Nq1K5KTkxEeHg4AKC0tRXJyMqKiosrVVygU/N8fERGRlavVPUSFhYVwcnKqsOzatWto1qxZnRpV36KjozFu3Dg8+uijeOyxx/Dhhx+ioKAAL730kqmbRkREZLasOUNUq4G8Rx55BOnp6eW2f/vtt+jUqVNd21TvRowYgRUrVmD+/PkIDg5Geno6EhMTy91oTURERP9jzdPuaxUQ9enTB927d8d7770HACgoKMD48eMxZswY/OMf/zBoA+tLVFQULl26BJ1Oh9TUVISEhJi6SURERGQitRoyW7duHQYPHoyJEydi586duHbtGlxdXfHLL78gKCjI0G0kIiIiM8BZZhUYOHAgnnvuOcTGxqJhw4bYsWMHgyEiKyFnqjkAJCUlVVrWuHFjTJ8+vdJynU5X7SwyLy+vattgCeQsY6BQKDijjcyeNd9DVKuA6OLFi3jxxReh1WqxZ88eHDhwAM888wymTp2Kd999Fw4ODoZuJxEZES/MhiVnGQM5ASgR1Z9a5a2Cg4MREBCAEydO4KmnnsI777yDn376Cd999x0ee+wxQ7eRiIiIzABvqn7AunXrsHnzZri5uUnbevTogfT0dDzyyCOGahsRERGZEWsOiGo1ZDZmzBgAwJkzZ5CVlYWioiKprKqVa4mIiIjMUa0Cot9++w3PPvssTp48CTs7OwghAPzvRqmSkhLDtZCIiIjMgjXfVF2rIbOpU6ciICAA169fh4uLC06fPo2DBw/i0Ucfxf79+w3cRCIiIjIHZdPu6/oyR7XKEKWkpGDfvn1o2rSp1LlevXohJiYGr7/+Oo4fP27odhIRERHVm1oFRCUlJWjcuDGAe0+Pv3r1Ktq1awd/f39kZGQYtIFERJZOzrpOfIi0eeHaURWz5iGzWgVEQUFBOHHiBAICAhASEoJly5bB0dERGzZsQMuWLQ3dRiIii2ZrF01rwLWjKmeuAU1d1Sogmjt3LgoKCgAAixcvxtNPP43evXujSZMm+Prrrw3aQCIiIqL6VquAKCwsTPp769atce7cOdy8eRPu7u5WGzkSERHZOg6ZyeDh4WGoXREREZEZsuaHu5pnq4iIiIiMyGAZIiIiIrJu1jxkxgwRERERyWKKZ5nFxsaiU6dOUCqVUCqV0Gg02L17t1ReWFiIyMhINGnSBK6urhg2bBiys7Nr3DdmiIiIiB7AtaPMR/PmzbF06VK0adMGQgh89tlnGDp0KI4fP46HH34Y06dPx65du7B161aoVCpERUXhueeew6FDh2p0HAZERERED+DaURUzxZDZkCFD9N6/++67iI2NxZEjR9C8eXNs3LgRmzZtQt++fQEAcXFxaN++PY4cOYLu3bvLPg6HzIiIiEgWUwyZ3a+kpASbN29GQUEBNBoN0tLScPfuXb0ANjAwEH5+fkhJSanRvpkhIiIiIqPLy8vTe69QKCodhjx58iQ0Gg0KCwvh6uqKbdu2oUOHDkhPT4ejoyPc3Nz06nt7e0Or1daoPcwQERERkSyGfNq9r68vVCqV9IqJian0uO3atUN6ejpSU1MRERGBcePG4cyZMwbtGzNEREREJIsh7yG6fPkylEqltL2qm9QdHR3RunVrAEDXrl1x9OhRrFq1CiNGjEBRURFycnL0skTZ2dlQq9U1ahczRERERGR0ZdPoy141mbVXWloKnU6Hrl27wsHBAcnJyVJZRkYGsrKyoNFoatQeZoiIiIhIFlPMMpszZw4GDhwIPz8/3L59G5s2bcL+/fuxZ88eqFQqTJgwAdHR0fDw8IBSqcSUKVOg0WhqNMMMYEBEREREMpkiILp+/TrGjh2La9euQaVSoVOnTtizZw+eeuopAMDKlSthb2+PYcOGQafTISwsDOvWratxuxgQERERkdnauHFjleVOTk5Yu3Yt1q5dW6fjMCAiIiIiWaz5afcMiIiIiEgWPtyViIiIyIoxQ0RERESyWHOGiAERERERyWLNARGHzIiIiMjmMUNEREREslhzhogBEREREcliZ2dX52nz5hoQcciMiIiIbB4zRERERCQLh8yIiIjI5llzQMQhMyIiIrJ5zBARERGRLNacIWJARGREe/fuhU6nq7KOQqFAaGiokVpERCQfH+5KRAah0+kwaNCgKuskJCQYqTVERFSGARERERHJYs1DZuaZtyIiIiIyIgZEREREZPM4ZEZERESyWPOQGQMiIiIikoUBEREZhEKhqHYWmUKhMFJriIioDAMiIiPi+kJEZMmYISIiIiKbZ80BEWeZERERkc1jhoiIiIhkseYMEQMiIiIiksWaAyKLGTJ799130aNHD7i4uMDNza3COllZWRg8eDBcXFzg5eWFmTNnori4WK/O/v378cgjj0ChUKB169aIj4+v/8YTERGRWbOYgKioqAjDhw9HREREheUlJSUYPHgwioqKcPjwYXz22WeIj4/H/PnzpTqZmZkYPHgwnnzySaSnp2PatGmYOHEi9uzZY6xuEBERWayyDFFdX+bIYobMFi1aBACVZnR+/PFHnDlzBnv37oW3tzeCg4Px9ttvY9asWVi4cCEcHR2xfv16BAQE4P333wcAtG/fHj///DNWrlyJsLAwY3WFiIiIzIzFZIiqk5KSgo4dO8Lb21vaFhYWhry8PJw+fVqq8+A6MGFhYUhJSTFqW4mIiMi8WEyGqDparVYvGAIgvddqtVXWycvLw99//w1nZ+dy+9XpdNDpdNL7vLw8QzediIjIIvCm6noye/bsascZz507Z8omIiYmBiqVSnr5+vqatD1ERESmwnuI6smMGTMwfvz4Kuu0bNlS1r7UajV++eUXvW3Z2dlSWdmfZdvur6NUKivMDgHAnDlzEB0dLb3Py8tjUERERGRlTBoQeXp6wtPT0yD70mg0ePfdd3H9+nV4eXkBAJKSkqBUKtGhQwepzoMP1kxKSoJGo6l0vwqFgg/bJCIiAofMzEJWVhbS09ORlZWFkpISpKenIz09Hfn5+QCA/v37o0OHDhgzZgxOnDiBPXv2YO7cuYiMjJQCmldffRW//fYb3nzzTZw7dw7r1q3Dli1bMH36dFN2jYiIyGJY43AZYEE3Vc+fPx+fffaZ9L5Lly4AgJ9++gl9+vRBgwYNsHPnTkRERECj0aBRo0YYN24cFi9eLH0mICAAu3btwvTp07Fq1So0b94c//znPznlnoiIyMbZCSGEqRthSfLy8qBSqZCbmwulUmnq5hARkY0y5vWo7Fj79u2Dq6trnfaVn5+Pvn37mt111GIyRERERGRavIeIiIiIyIoxQ0RERESyMENEREREZMUYEBEREZHN45AZERERyWLNQ2YMiIiIiEgWaw6IOGRGREREZismJgbdunVD48aN4eXlhfDwcGRkZOjVKSwsRGRkJJo0aQJXV1cMGzas3LNLq8OAiIiIiGQxxdPuDxw4gMjISBw5cgRJSUm4e/cu+vfvj4KCAqnO9OnTsWPHDmzduhUHDhzA1atX8dxzz9XoOBwyIyIiIllMMWSWmJio9z4+Ph5eXl5IS0vD448/jtzcXGzcuBGbNm1C3759AQBxcXFo3749jhw5gu7du8s6DjNEREREZDFyc3MBAB4eHgCAtLQ03L17F6GhoVKdwMBA+Pn5ISUlRfZ+mSEiIiIiWQyZIcrLy9PbrlAooFAoqvxsaWkppk2bhp49eyIoKAgAoNVq4ejoCDc3N7263t7e0Gq1stvFDBERERHJYsh7iHx9faFSqaRXTExMtcePjIzEqVOnsHnzZoP3jRkiIiIiMrrLly/rPe2+uuxQVFQUdu7ciYMHD6J58+bSdrVajaKiIuTk5OhlibKzs6FWq2W3hwERkZnZu3cvdDpdpeUKhUJvrJyIyFgMOWSmVCr1AqLKCCEwZcoUbNu2Dfv370dAQIBeedeuXeHg4IDk5GQMGzYMAJCRkYGsrCxoNBrZ7WJARGRmdDodBg0aVGl5QkKCEVtDRPQ/pphlFhkZiU2bNuH7779H48aNpfuCVCoVnJ2doVKpMGHCBERHR8PDwwNKpRJTpkyBRqORPcMMYEBEREREZiw2NhYA0KdPH73tcXFxGD9+PABg5cqVsLe3x7Bhw6DT6RAWFoZ169bV6DgMiIiIiMhsCSGqrePk5IS1a9di7dq1tT4OAyIiIiKShc8yIyIiIrJizBARERGRLNacIWJARERERLIwICIio1EoFFVOra9u8TIiIqo5BkREZoaLLhKRuWKGiIiIiGyeNQdEnGVGRERENo8ZIiIiIpKFGSIiIiIiK8aAiIiIiGweh8yIiIhINnMd8qorBkREREQkC+8hIiIiIrJiDIiIiIjI5nHIjIiIiGThkBkRERGRFWOGiIiIiGSx5gwRAyIiIiKSxZoDIg6ZERERkc1jQEREREQ2j0NmREREJAuHzIiIiIisGDNEREREJAszRERERERWjAERERER2TwOmREREZEs1jxkxoCI6P/s3bsXOp2u0nKFQoHQ0FAjtoiIyLwwICKyATqdDoMGDaq0PCEhwYitISIiY2JARERERLJYc4aIN1UTERGRzWNARERERDaPQ2ZEREQkizUPmTEgIiIiIlmsOSDikBkRERHZPIsIiH7//XdMmDABAQEBcHZ2RqtWrbBgwQIUFRXp1fv111/Ru3dvODk5wdfXF8uWLSu3r61btyIwMBBOTk7o2LEjp1KTRKFQICEhodKXQqEwdROJiEyqLENU15c5soghs3PnzqG0tBQff/wxWrdujVOnTmHSpEkoKCjAihUrAAB5eXno378/QkNDsX79epw8eRIvv/wy3NzcMHnyZADA4cOHMWrUKMTExODpp5/Gpk2bEB4ejv/85z8ICgoyZRfJDHDRRSIi22UnhBCmbkRtLF++HLGxsfjtt98AALGxsXjrrbeg1Wrh6OgIAJg9eza2b9+Oc+fOAQBGjBiBgoIC7Ny5U9pP9+7dERwcjPXr18s6bl5eHlQqFXJzc6FUKg3cKyIiInmMeT0qO9Z///tfNG7cuE77un37Ntq2bWt211GLGDKrSG5uLjw8PKT3KSkpePzxx6VgCADCwsKQkZGBW7duSXUezAKEhYUhJSXFOI0mIiKyYBwyMzMXLlzAmjVrpOEyANBqtQgICNCr5+3tLZW5u7tDq9VK2+6vo9VqKz2WTqfTe75VXl6eIbpARERkcTjLrJ7Mnj272iiybLirzJUrVzBgwAAMHz4ckyZNqvc2xsTEQKVSSS9fX996PyYREREZl0kzRDNmzMD48eOrrNOyZUvp71evXsWTTz6JHj16YMOGDXr11Go1srOz9baVvVer1VXWKSuvyJw5cxAdHS29z8vLY1BERERkZUwaEHl6esLT01NW3StXruDJJ59E165dERcXB3t7/eSWRqPBW2+9hbt378LBwQEAkJSUhHbt2sHd3V2qk5ycjGnTpkmfS0pKgkajqfS4CoWC062JiIjAITOTu3LlCvr06QM/Pz+sWLECf/75J7Rard69Py+++CIcHR0xYcIEnD59Gl9//TVWrVqll92ZOnUqEhMT8f777+PcuXNYuHAhjh07hqioKFN0i4iIiMyERQRESUlJuHDhApKTk9G8eXM0a9ZMepVRqVT48ccfkZmZia5du2LGjBmYP3++tAYRAPTo0QObNm3Chg0b0LlzZ3zzzTfYvn071yAiIiKSwRSzzA4ePIghQ4bAx8cHdnZ22L59u165EALz589Hs2bN4OzsjNDQUJw/f77GfbOIgGj8+PEQQlT4ul+nTp3w73//G4WFhfjjjz8wa9ascvsaPnw4MjIyoNPpcOrUKQwaNMhY3SAiIrJopgiICgoK0LlzZ6xdu7bC8mXLlmH16tVYv349UlNT0ahRI4SFhaGwsLBGx7HIafdERERkGwYOHIiBAwdWWCaEwIcffoi5c+di6NChAIDPP/8c3t7e2L59O0aOHCn7OBaRISIiIiJ6UGZmJrRard6iyyqVCiEhITVedJkZIiIiIpLNULPEHlzouDazussmV9V00eWKMENERERERufr66u38HFMTIxJ28MMEREREcliyHWILl++rPdw19qs+Ve2sHJ2drbezPPs7GwEBwfXaF/MEBEREZHRKZVKvVdtAqKAgACo1WokJydL2/Ly8pCamlrlossVYYaIiIiIzFZ+fj4uXLggvc/MzER6ejo8PDzg5+eHadOm4Z133kGbNm0QEBCAefPmwcfHB+Hh4TU6DgMiIiIiksUUj+44duwYnnzySel92RMoxo0bh/j4eLz55psoKCjA5MmTkZOTg169eiExMRFOTk41a5d4cHVDqlJeXh5UKhVyc3P1xj6JiIiMyZjXo7JjZWVl1flYeXl58PPzM7vrKO8hIiIiIpvHgIiIiIhsHu8hIiIiIllMcQ+RsTBDRERERDaPARERERHZPA6ZERERkSwcMiMiIiKyYgyIiIiIyOZxyIyIiIhk4ZAZERERkRVjQEREREQ2j0NmREREJAuHzIiIiIisGAMiIiIisnkcMiMiIiJZOGRGREREZMUYEBEREZHN45AZERERycIhMyIiIiIrxoCIiIiIbB6HzIiIiEgWDpkRERERWTEGRERERGTzGBARERGRzWNARERERDaPN1UTERGRLLypmoiIiMiKMSAiIiIim8chMyIiIpKFQ2ZEREREVowBEREREdk8DpkRERGRLBwyIyIiIrJiDIiIiIjI5nHIjIiIiGThkBkRERGRFWNARERERDaPQ2ZEREQkC4fMiIiIiKwYAyIiIiKyeRwyIyIiIlk4ZEZERERkxRgQERERkc3jkBkRERHJwiEzM/DMM8/Az88PTk5OaNasGcaMGYOrV6/q1fn111/Ru3dvODk5wdfXF8uWLSu3n61btyIwMBBOTk7o2LEjEhISjNUFIiIiqoW1a9eiRYsWcHJyQkhICH755ReDH8NiAqInn3wSW7ZsQUZGBr799ltcvHgRzz//vFSel5eH/v37w9/fH2lpaVi+fDkWLlyIDRs2SHUOHz6MUaNGYcKECTh+/DjCw8MRHh6OU6dOmaJLREREVI2vv/4a0dHRWLBgAf7zn/+gc+fOCAsLw/Xr1w16HDshhDDoHo3khx9+QHh4OHQ6HRwcHBAbG4u33noLWq0Wjo6OAIDZs2dj+/btOHfuHABgxIgRKCgowM6dO6X9dO/eHcHBwVi/fr2s4+bl5UGlUiE3NxdKpdLwHSMiIpLBmNcjQx6rpvsKCQlBt27d8NFHHwEASktL4evriylTpmD27Nl1asv9LCZDdL+bN2/iq6++Qo8ePeDg4AAASElJweOPPy4FQwAQFhaGjIwM3Lp1S6oTGhqqt6+wsDCkpKQYr/FEREQkS1FREdLS0vSu3fb29ggNDTX4tduibqqeNWsWPvroI9y5cwfdu3fXy/RotVoEBATo1ff29pbK3N3dodVqpW3319FqtZUeU6fTQafTSe9zc3MB3ItwiYiITKXsOmTMgR5DXPvK9vHgvhQKBRQKhd62v/76CyUlJRVeu8tGfwzFpAHR7Nmz8d5771VZ5+zZswgMDAQAzJw5ExMmTMClS5ewaNEijB07Fjt37qzXO9ZjYmKwaNGictt9fX3r7ZhERERy3bhxAyqVql6P4ejoCLVabbBrn6ura7l9LViwAAsXLjTI/mvDpAHRjBkzMH78+CrrtGzZUvp706ZN0bRpU7Rt2xbt27eHr68vjhw5Ao1GA7VajezsbL3Plr1Xq9XSnxXVKSuvyJw5cxAdHS29z8nJgb+/P7Kysur9H6Cx5OXlwdfXF5cvX7aK+6KsrT+A9fXJ2voDWF+frK0/gPX1KTc3F35+fvDw8Kj3Yzk5OSEzMxNFRUUG2Z8Qolwy48HsEHDvut+gQYMaX7trw6QBkaenJzw9PWv12dLSUgCQhrM0Gg3eeust3L17V7qvKCkpCe3atYO7u7tUJzk5GdOmTZP2k5SUBI1GU+lxKkrhAYBKpbKKX6j7KZVKq+qTtfUHsL4+WVt/AOvrk7X1B7C+PtnbG+d2YCcnJzg5ORnlWGUcHR3RtWtXJCcnIzw8HMC9639ycjKioqIMeiyLuKk6NTUVH330EdLT03Hp0iXs27cPo0aNQqtWraRg5sUXX4SjoyMmTJiA06dP4+uvv8aqVav0sjtTp05FYmIi3n//fZw7dw4LFy7EsWPHDP5DJSIiIsOIjo7GJ598gs8++wxnz55FREQECgoK8NJLLxn0OBZxU7WLiwu+++47LFiwAAUFBWjWrBkGDBiAuXPnStkblUqFH3/8EZGRkejatSuaNm2K+fPnY/LkydJ+evTogU2bNmHu3Ln4xz/+gTZt2mD79u0ICgoyVdeIiIioCiNGjMCff/6J+fPnQ6vVIjg4GImJieVutK4riwiIOnbsiH379lVbr1OnTvj3v/9dZZ3hw4dj+PDhtW6LQqHAggULKhxGs1TW1idr6w9gfX2ytv4A1tcna+sPYH19srb+VCUqKqreR3MsdmFGIiIiIkOxiHuIiIiIiOoTAyIiIiKyeQyIiIiIyOYxIPo/7777Lnr06AEXFxe4ublVWCcrKwuDBw+Gi4sLvLy8MHPmTBQXF+vV2b9/Px555BEoFAq0bt0a8fHx5fazdu1atGjRAk5OTggJCcEvv/xSDz3St3//ftjZ2VX4Onr0KADg999/r7D8yJEjevvaunUrAgMD4eTkhI4dOyIhIaHe21+ZFi1alGvv0qVL9er8+uuv6N27N5ycnODr64tly5aV24859On333/HhAkTEBAQAGdnZ7Rq1QoLFizQWwjNEs9RRUzxO1AbMTEx6NatGxo3bgwvLy+Eh4cjIyNDr06fPn3KnY9XX31Vr46c7w5jWLhwYbm2lj0JAAAKCwsRGRmJJk2awNXVFcOGDSu3IJ659KVMRd8BdnZ2iIyMBGD+5+fgwYMYMmQIfHx8YGdnh+3bt+uVCyEwf/58NGvWDM7OzggNDcX58+f16ty8eROjR4+GUqmEm5sbJkyYgPz8fL06cr4HbZ4gIYQQ8+fPFx988IGIjo4WKpWqXHlxcbEICgoSoaGh4vjx4yIhIUE0bdpUzJkzR6rz22+/CRcXFxEdHS3OnDkj1qxZIxo0aCASExOlOps3bxaOjo7i008/FadPnxaTJk0Sbm5uIjs7u177p9PpxLVr1/ReEydOFAEBAaK0tFQIIURmZqYAIPbu3atXr6ioSNrPoUOHRIMGDcSyZcvEmTNnxNy5c4WDg4M4efJkvba/Mv7+/mLx4sV67c3Pz5fKc3Nzhbe3txg9erQ4deqU+Ne//iWcnZ3Fxx9/LNUxlz7t3r1bjB8/XuzZs0dcvHhRfP/998LLy0vMmDFDqmOJ5+hBpvodqI2wsDARFxcnTp06JdLT08WgQYOEn5+f3r+xJ554QkyaNEnvfOTm5krlcr47jGXBggXi4Ycf1mvrn3/+KZW/+uqrwtfXVyQnJ4tjx46J7t27ix49ephlX8pcv35drz9JSUkCgPjpp5+EEOZ/fhISEsRbb70lvvvuOwFAbNu2Ta986dKlQqVSie3bt4sTJ06IZ555RgQEBIi///5bqjNgwADRuXNnceTIEfHvf/9btG7dWowaNUoql/M9SEIwIHpAXFxchQFRQkKCsLe3F1qtVtoWGxsrlEql0Ol0Qggh3nzzTfHwww/rfW7EiBEiLCxMev/YY4+JyMhI6X1JSYnw8fERMTExBu5J1YqKioSnp6dYvHixtK3sYnv8+PFKP/fCCy+IwYMH620LCQkRr7zySn01tUr+/v5i5cqVlZavW7dOuLu7S+dICCFmzZol2rVrJ703tz7db9myZSIgIEB6b4nn6EHm8jtQG9evXxcAxIEDB6RtTzzxhJg6dWqln5Hz3WEsCxYsEJ07d66wLCcnRzg4OIitW7dK286ePSsAiJSUFCGEefWlMlOnThWtWrWS/qNnSefnwYCotLRUqNVqsXz5cmlbTk6OUCgU4l//+pcQQogzZ84IAOLo0aNSnd27dws7Oztx5coVIYS870ESgkNmMqWkpKBjx456C0GFhYUhLy8Pp0+fluqEhobqfS4sLAwpKSkAgKKiIqSlpenVsbe3R2hoqFTHWH744QfcuHGjwpU+n3nmGXh5eaFXr1744Ycf9Mqq66MpLF26FE2aNEGXLl2wfPlyvVR3SkoKHn/8cTg6OkrbwsLCkJGRgVu3bkl1zK1PZXJzcyt8TpGlnaMy5vQ7UBu5ubkAUO6cfPXVV2jatCmCgoIwZ84c3LlzRyqT891hTOfPn4ePjw9atmyJ0aNHIysrCwCQlpaGu3fv6p2bwMBA+Pn5SefG3PryoKKiInz55Zd4+eWX9Z6TZUnn536ZmZnQarV650SlUiEkJETvnLi5ueHRRx+V6oSGhsLe3h6pqalSneq+B8lCFmY0B1qtttyqmGXvtVptlXXy8vLw999/49atWygpKamwzrlz5+qx9eVt3LgRYWFhaN68ubTN1dUV77//Pnr27Al7e3t8++23CA8Px/bt2/HMM88AqLyPZT8DY3v99dfxyCOPwMPDA4cPH8acOXNw7do1fPDBBwDutTcgIEDvM/efN3d3d7PrU5kLFy5gzZo1WLFihbTNEs/R/f766y+z+R2oqdLSUkybNg09e/bUW93+xRdfhL+/P3x8fPDrr79i1qxZyMjIwHfffQdA3neHsYSEhCA+Ph7t2rXDtWvXsGjRIvTu3RunTp2CVquFo6NjuXso7/+3Y059qcj27duRk5Oj99BwSzo/Dyo7flW/z1qtFl5eXnrlDRs2hIeHh16d6r4HycoDotmzZ+O9996rss7Zs2f1biq0NLXp4x9//IE9e/Zgy5YtevWaNm2q9+y3bt264erVq1i+fLl0sTWGmvTp/vZ26tQJjo6OeOWVVxATE2M2q7fW5hxduXIFAwYMwPDhwzFp0iRpu7mcI1sUGRmJU6dO4eeff9bbfv/jgTp27IhmzZqhX79+uHjxIlq1amXsZlZp4MCB0t87deqEkJAQ+Pv7Y8uWLXB2djZhywxj48aNGDhwIHx8fKRtlnR+yLSsOiCaMWOG3v8UKtKyZUtZ+1Kr1eVmwpTNvlCr1dKfD87IyM7OhlKphLOzMxo0aIAGDRpUWKdsHzVVmz7GxcWhSZMmsi6gISEhSEpKkt5X1sfatr8idTlvISEhKC4uxu+//4527dpV2l6g+vNmqD7VtD9Xr17Fk08+iR49emDDhg3V7t8U56i2mjZtavDfAWOIiorCzp07cfDgQb2sakVCQkIA3MvwtWrVStZ3h6m4ubmhbdu2uHDhAp566ikUFRUhJydHL0t0/7kx575cunQJe/fulTI/lbGk81N2/OzsbDRr1kzanp2djeDgYKnO9evX9T5XXFyMmzdvVvsdd/8xCJxl9qDqbqq+fybMxx9/LJRKpSgsLBRC3LupOigoSO9zo0aNKndTdVRUlPS+pKREPPTQQ0a7obS0tFQEBATozVyqysSJE0WXLl2k9y+88IJ4+umn9epoNBqzuWH3yy+/FPb29uLmzZtCiP/dTHj/LKw5c+aUu6naXPr0xx9/iDZt2oiRI0eK4uJiWZ+xtHNk6t+BmigtLRWRkZHCx8dH/Pe//5X1mZ9//lkAECdOnBBCyPvuMJXbt28Ld3d3sWrVKumm6m+++UYqP3fuXIU3VZtjXxYsWCDUarW4e/dulfXM+fygkpuqV6xYIW3Lzc2t8KbqY8eOSXX27NlT4U3VVX0PEmeZSS5duiSOHz8uFi1aJFxdXcXx48fF8ePHxe3bt4UQ/5ua2b9/f5Geni4SExOFp6dnhdPuZ86cKc6ePSvWrl1b4bR7hUIh4uPjxZkzZ8TkyZOFm5ub3gyH+rR3714BQJw9e7ZcWXx8vNi0aZM4e/asOHv2rHj33XeFvb29+PTTT6U6hw4dEg0bNhQrVqwQZ8+eFQsWLDDZlO7Dhw+LlStXivT0dHHx4kXx5ZdfCk9PTzF27FipTk5OjvD29hZjxowRp06dEps3bxYuLi7lpt2bQ5/++OMP0bp1a9GvXz/xxx9/6E0TLmNp56gipv4dqImIiAihUqnE/v379c7HnTt3hBBCXLhwQSxevFgcO3ZMZGZmiu+//160bNlSPP7449I+5Hx3GMuMGTPE/v37RWZmpjh06JAIDQ0VTZs2FdevXxdC3Jt27+fnJ/bt2yeOHTsmNBqN0Gg0ZtmX+5WUlAg/Pz8xa9Ysve2WcH5u374tXW8AiA8++EAcP35cXLp0SQhxb9q9m5ub+P7778Wvv/4qhg4dWuG0+y5duojU1FTx888/izZt2uhNu5fzPUgMiCTjxo0TAMq9ytayEEKI33//XQwcOFA4OzuLpk2bihkzZpT738hPP/0kgoODhaOjo2jZsqWIi4srd6w1a9YIPz8/4ejoKB577DFx5MiReu7d/4waNUpvXZH7xcfHi/bt2wsXFxehVCrFY489pjcFt8yWLVtE27ZthaOjo3j44YfFrl276rvZFUpLSxMhISFCpVIJJycn0b59e7FkyZJy/6s7ceKE6NWrl1AoFOKhhx4SS5cuLbcvc+hTXFxchf8G70/kWto5qowpfwdqorLzUfZ7nZWVJR5//HHh4eEhFAqFaN26tZg5c6beOjdCyPvuMIYRI0aIZs2aCUdHR/HQQw+JESNGiAsXLkjlf//9t3jttdeEu7u7cHFxEc8++6xeQC6E+fTlfnv27BEAREZGht52Szg/P/30U4X/xsaNGyeEuJclmjdvnvD29hYKhUL069evXD9v3LghRo0aJVxdXYVSqRQvvfSS9J/5MnK+B20dn3ZPRERENo/rEBEREZHNY0BERERENo8BEREREdk8BkRERERk8xgQERERkc1jQEREREQ2jwERERER2TwGRERERGTzGBARkUm0aNECH374oambQUQEgAEREZnI0aNHMXnyZFl1hRCYP38+mjVrBmdnZ4SGhuL8+fP13EIisiUMiIjIJDw9PeHi4iKr7rJly7B69WqsX78eqampaNSoEcLCwlBYWFjPrSQiW8GAiIhqJTExEb169YKbmxuaNGmCp59+GhcvXgQAfP7553B1ddXL4rz22msIDAzEnTt3AOgPmQkhsHDhQvj5+UGhUMDHxwevv/66VPbhhx9i7ty5GDp0KDp16oTPP/8cV69exfbt243aZyKyXgyIiKhWCgoKEB0djWPHjiE5ORn29vZ49tlnUVpairFjx2LQoEEYPXo0iouLsWvXLvzzn//EV199VWFW6Ntvv8XKlSvx8ccf4/z589i+fTs6duwIAMjMzIRWq0VoaKhUX6VSISQkBCkpKUbrLxFZt4ambgARWaZhw4bpvf/000/h6emJM2fOICgoCB9//DE6deqE119/Hd999x0WLlyIrl27VrivrKwsqNVqhIaGwsHBAX5+fnjssccAAFqtFgDg7e2t9xlvb2+pjIiorpghIqJaOX/+PEaNGoWWLVtCqVSiRYsWAO4FNwDg7u6OjRs3IjY2Fq1atcLs2bMr3dfw4cPx999/o2XLlpg0aRK2bduG4uJiY3SDiAgAAyIiqqUhQ4bg5s2b+OSTT5CamorU1FQAQFFRkVTn4MGDaNCgAa5du4aCgoJK9+Xr64uMjAysW7cOzs7OeO211/D444/j7t27UKvVAIDs7Gy9z2RnZ0tlRER1xYCIiGrsxo0byMjIwNy5c9GvXz+0b98et27d0qtz+PBhvPfee9ixYwdcXV0RFRVV5T6dnZ0xZMgQrF69Gvv370dKSgpOnjyJgIAAqNVqJCcnS3Xz8vKQmpoKjUZTL/0jItvDe4iIqMbc3d3RpEkTbNiwAc2aNUNWVpbekNjt27cxZswYvP766xg4cCCaN2+Obt26YciQIXj++efL7S8+Ph4lJSUICQmBi4sLvvzySzg7O8Pf3x92dnaYNm0a3nnnHbRp0wYBAQGYN28efHx8EB4ebsReE5E1Y4aIiGrM3t4emzdvRlpaGoKCgjB9+nQsX75cKp86dSoaNWqEJUuWAAA6duyIJUuW4JVXXsGVK1fK7c/NzQ2ffPIJevbsiU6dOmHv3r3YsWMHmjRpAgB48803MWXKFEyePBndunVDfn4+EhMT4eTkZJwOE5HVsxNCCFM3goiIiMiUmCEiIiIim8eAiIiIiGweAyIiIiKyeQyIiIiIyOYxICIiIiKbx4CIiIiIbB4DIiIiIrJ5DIiIiIjI5jEgIiIiIpvHgIiIiIhsHgMiIiIisnkMiIiIiMjm/X/RyiAxhEBK6AAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h2.plot(show_zero=False)" ] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.5" } }, "nbformat": 4, "nbformat_minor": 1 }