Skip to content

Instantly share code, notes, and snippets.

@alexander-held
Last active August 21, 2023 13:59
Show Gist options
  • Save alexander-held/9eb02d00986ebfbc908a887d8df64ef9 to your computer and use it in GitHub Desktop.
Save alexander-held/9eb02d00986ebfbc908a887d8df64ef9 to your computer and use it in GitHub Desktop.
model predictions for custom parameter values
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "a3425b9c",
"metadata": {},
"source": [
"## model predictions for custom parameter values"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "a9ffd410",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[{'figure': <Figure size 600x600 with 2 Axes>, 'region': 'singlechannel'}]"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmMAAAJjCAYAAABAyBZ9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACP1UlEQVR4nOzdd1hT1/8H8HdANoIiICAoiKIiKg7QWlfdinuvWidVcba2aq3FUVetq4paR9G6LVXrrFvrqtvWgcW6CogiKiCyyfn9wY98jQElIeRC8n49T54295x77uckknxy77nnyIQQAkREREQkCSOpAyAiIiIyZEzGiIiIiCTEZIyIiIhIQkzGiIiIiCTEZIyIiIhIQkzGiIiIiCTEZIyIiIhIQiWkDkAfyOVyPH78GCVLloRMJpM6HCIiItKQEAKvXr2Ci4sLjIx0c86KyZgWPHz4EJ6enlKHQURERFpy7949VKxYUSfHYjKmBSYmJgCA27dvo1y5cvnax8/PD5cuXVLrOOruU9jHSExMhJubGyIjI2FjY1NocenitVJ3H33qu7r1Nel7UeyHuvvo6j3XZB996XtRfK0M9d87YLh9j46Ohre3t+K7XReYjGlBzqXJkiVL5vsfrLGxsVofaprso4tjAICNjU2hxqWrfhhq3zU5BqBe34tqP4rie67JPvrS96L6WgGG++8dMLy+JyYmAoBOhx1xAL9EgoKCCn0fXRxDE0W1H4ba96LYb0320Zf3XJN99KXvRfW1UldR7Qf7XnTJuFB4wUVFRSlO5bq6ukodjs4kJibC1tYWCQkJGv3aKs7Yd8Pru6H2G2Df2XfD6rsU3+k8M6YFZmZmSv81FGZmZggODja4fgPsuyH23VD7DbDv7Lth9V2K73SeGdMCQ/31QEREpG+k+E4vFmfGkpKSEBwcjLZt28LOzg4ymQzr16/Pta5cLsfKlSvh6+sLCwsLlClTBs2bN8dff/2lUu+7776Dh4cHzM3NUbNmTWzdulUHvSEiIiL6n2KRjMXFxWHmzJkIDw9HrVq13ll3yJAhGDt2LOrWrYtly5bhm2++Qfny5REbG6tUb+rUqZg0aRJatWqFZcuWoXz58ujXrx+2bdtWmF0hIiIiUlIsLlOmpaXh5cuXcHJywuXLl+Hn54fQ0FAMGjRIqd6OHTvQu3dv7Ny5E127ds2zvejoaHh4eCAwMBDLly8HkD3jbtOmTfHgwQM8fPgQxsbG+Y6PlymJiIj0Ay9T5sHMzAxOTk7vrbdo0SL4+/uja9eukMvleP36da71fvvtN2RkZGDUqFGKbTKZDCNHjkRUVBTOnz+vtdiJiIiI3qVYJGP5kZiYiIsXL8LPzw9fffUVbG1tYW1tjYoVK2LHjh1Kda9duwYrKytUq1ZNabu/v7+iXNMY3nykpaVp1hkiIiLSibS0NJXvb13Tm2Ts3r17EEJg27Zt+Omnn/Ddd99h8+bNcHBwQJ8+ffD7778r6sbExKBs2bIqs+s6OzsDAB4/fqxRDG5ubrC1tVU85s6dq3mHiIiIqNDNnTtX6bvbzc1N5zHozXJISUlJAIDnz5/jzz//RP369QEAnTp1goeHB7799lu0bdsWAJCSkpLr/CHm5uaKck28vX6Xoc3NQkREVNxMmTIFn332meJ5zpqcuqQ3yZiFhQUAwMPDQ5GIAYC1tTU6duyITZs2ITMzEyVKlICFhUWulxBTU1OV2lKXuuu2ERERkbTMzMwkP3miN5cpXVxcAABly5ZVKXN0dERGRoZiQL+zszOePHmCt28kjYmJUWqLiIiIqLDpVTLm5OSE6OholbLHjx/D3NwcJUuWBAD4+voiOTkZ4eHhSvUuXLigKC+KBg0aBHd390I/jkwmw/Tp0wul7YcPH0Imk+H7778vlPZ1Lac/eU1CTESUF119pueXu7u7ypRRpBt6k4wBQO/evREZGYkjR44otsXFxeG3335D8+bNYWSU3d3OnTvDxMQEK1asUNQTQmDVqlUoV64cGjZsqPPYiYiIyDAVmzFjy5cvR3x8vOJOx7179yIqKgoAMGbMGNja2mLKlCnYsWMHunfvjs8++wy2trZYtWoVMjIyMGfOHEVbrq6uGD9+PBYsWICMjAz4+flh9+7dOH36NDZv3qzWhK+6tGbNGsjlcqnDICIqdLLQiVKHAAAQgwvvLD4/0ylHsUnGvv/+ezx69EjxfOfOndi5cycAYMCAAbC1tUXZsmVx5swZTJw4EYsXL0ZGRgY++OADbNq0SWUZpXnz5qF06dL48ccfsX79elSuXBmbNm1Cv379dNovdZiYmEgdAhERaQk/0ylHsblM+fDhQwghcn28ec29YsWK2LlzJxISEpCcnIxjx47Bz89PpT0jIyNMmTIFDx8+RFpaGm7evIn+/fvrsEeqXr16hfHjx8Pd3R1mZmZwdHREq1atcPXqVQCq4wveHH+1evVqeHp6wszMDH5+frh06ZJK+7/88gu8vb1hbm4OHx8f7Nq1K99jFqKjozFkyBCULVsWZmZmqF69On766SeVeqmpqZg+fTq8vLxgbm4OZ2dndOvWDffu3VOp+76Y//77bwwaNAgVK1aEubk5nJycMGTIEDx//lyp3vTp0yGTyfDvv/9i0KBBKFWqFGxtbTF48GAkJycr1ZXJZBg9ejR2794NHx8fRV/enIdO3T4TEeVG3c90IHt6po8//hg2NjYoVaoUPvnkE/z1118qY1MHDRoEa2trREdHo0uXLrC2toaDgwMmTpyIrKwspTa///57NGzYEGXKlIGFhQXq1q2LsLCwwu4+qaHYnBkzBCNGjEBYWBhGjx4Nb29vPH/+HGfOnEF4eDjq1KmT535btmzBq1ev8Omnn0Imk+G7775Dt27dcP/+fcUvr/3796N3796oUaMG5s6di5cvX2Lo0KEoV67ce+N6+vQpGjRooEhkHBwccPDgQQwdOhSJiYkYP348ACArKwsdOnTAsWPH0KdPH4wbNw6vXr3CkSNHcPPmTXh6eqoV85EjR3D//n0MHjwYTk5OuHXrFlavXo1bt27hzz//VJm0t1evXvDw8MDcuXNx9epVrF27Fo6Ojpg/f75SvTNnzmDnzp0YNWoUSpYsiR9++AHdu3fHf//9hzJlyqjVZyKivKj7mS6Xy9GxY0dcvHgRI0eORNWqVfHbb7/hk08+ybX9rKwstGnTBvXr18f333+Po0ePYuHChfD09MTIkSMV9ZYuXYpOnTqhf//+SE9Px7Zt29CzZ0/s27cPAQEBhdZ/yj8mY0XI/v37MXz4cCxcuFCx7csvv3zvfv/99x/u3r2L0qVLAwCqVKmCzp0749ChQ+jQoQOA7EntypUrh7Nnz8La2hoA0KJFCzRr1gwVKlR4Z/tTp05FVlYWbty4oUhWRowYgb59+2L69On49NNPYWFhgZ9//hnHjh3DokWLMGHCBMX+kydPVplGJD8xjxo1Cp9//rnSfg0aNEDfvn1x5swZNG7cWKmsdu3aWLduneL58+fPsW7dOpVkLDw8HLdv31Ykhx999BFq1aqFrVu3YvTo0Wr1mYgoL+p+pu/evRvnz5/HkiVLMG7cOADAyJEj0apVq1zrp6amonfv3pg2bRqA7M+oOnXqYN26dUrJWEREhNLn1ejRo1GnTh0sWrSIyVgRUWwuUxqCUqVK4cKFC2ovx9S7d29FUgNAkaTcv38fQPbUHjdu3MDAgQMViRgANG3aFDVq1Hhn20II/Prrr+jYsSOEEIiLi1M82rRpg4SEBMUp919//RX29vYYM2aMSjtvn8V6X8yA8uS7qampiIuLQ4MGDQBAccw3jRgxQul548aN8fz5c5V1xlq2bKl0lq5mzZqwsbFRHFudPhMR5UXdz/Tff/8dJiYmGD58uGKbkZERgoKC8twnt8+9Nz9HAeXP0pcvXyIhIQGNGzfm51gRwmRMi/z8/ODt7Y2QkBCN9v/uu+9w8+ZNuLm5wd/fH9OnT1f5o8pN+fLllZ7nJDkvX74EAMWND5UqVVLZN7dtb3r27Bni4+OxevVqODg4KD0GDx4MAIiNjQWQvT5olSpVUKLE+0+4vi9mAHjx4gXGjRuHsmXLwsLCAg4ODvDw8AAAJCQkaNRmbvVy6ubUU6fPRER5Ufcz/dGjR3B2doalpaXS9rw+p83NzeHg4KC07c3Pshz79u1DgwYNYG5uDjs7Ozg4OGDlypW5fo4aspCQEHh7e+c6zryw8TKlFl26dKlAyyH16tULjRs3xq5du3D48GEsWLAA8+fPx86dO9GuXbs898trKo63Lw1qIue26wEDBuQ5bqFmzZpqt5ufmHv16oVz587hiy++gK+vL6ytrSGXy9G2bdtcbwfP7+vwvnqF1WciMiyafqbnV36mYTp9+jQ6deqEJk2aYMWKFXB2doaJiQlCQ0OxZcuWAsegT4KCghAUFITExETY2trq9NhMxooYZ2dnjBo1CqNGjUJsbCzq1KmD2bNnF+gPN2dM2L///qtSltu2Nzk4OKBkyZLIyspCy5Yt31nX09MTFy5cQEZGRoFv2X758iWOHTuGGTNm4JtvvlFsv3v3boHazQ91+kxE9C7qfKZXqFABJ06cQHJystLZsfd9Tr/Lr7/+CnNzcxw6dEhp/cXQ0FCN2yTt42XKIiIrK0vllLGjoyNcXFxyXdRcHS4uLvDx8cHPP/+MpKQkxfZTp07hxo0b79zX2NgY3bt3x6+//oqbN2+qlD979kzx/927d0dcXByWL1+uUk/ds3Q5v/je3m/JkiVqtaMJdfpMRJQbTT7T27Rpg4yMDKxZs0axTS6Xazz0Bcj+PJPJZErTXTx8+BC7d+/WuE3SPp4ZKyJevXoFV1dX9OjRA7Vq1YK1tTWOHj2KS5cuKd2Jo6k5c+agc+fO+PDDDzF48GC8fPkSy5cvh4+Pj1KClpt58+bhxIkTqF+/PoYPHw5vb2+8ePECV69exdGjR/HixQsAwMCBA/Hzzz/js88+w8WLF9G4cWO8fv0aR48exahRo9C5c+d8x2tjY4MmTZrgu+++Q0ZGBsqVK4fDhw/jwYMHBXod8iu/fSYiyo0mn+ldunSBv78/Pv/8c/z777+oWrUq9uzZo/i8eftGqPwICAjAokWL0LZtW/Tr1w+xsbEICQlBpUqV8Pfffxeoj6Q9TMaKCEtLS4waNQqHDx/Gzp07IZfLUalSJaxYsULpFmVNdezYEVu3bsX06dMxefJkVK5cGevXr8eGDRtw69atd+5btmxZXLx4ETNnzsTOnTuxYsUKlClTBtWrV1eaNsLY2BgHDhzA7NmzsWXLFvz6668oU6YMGjVq9N67NnOzZcsWjBkzBiEhIRBCoHXr1jh48CBcXFzUbktd+e0zERWOwlyGSBc0+Uw3NjbG/v37MW7cOGzYsAFGRkbo2rUrgoOD8eGHH8Lc3FztOJo3b45169Zh3rx5GD9+PDw8PDB//nw8fPiQyVgRIhPaGOVt4HIG+yUkJBRoAL8UfH194eDgoLS4OhERFR27d+9G165dcebMGXz44YdSh6P3pPhO55gxA5GRkYHMzEylbSdPnsRff/2FZs2aSRMUEREpSUlJUXqelZWFZcuWwcbG5p0rsVDxxsuUBiI6OhotW7bEgAED4OLigjt37mDVqlVwcnJSmTSQiIikMWbMGKSkpOCDDz5AWloadu7ciXPnzmHOnDlc9UOPMRkzEKVLl0bdunWxdu1aPHv2DFZWVggICMC8efMUy/0QEZG0mjdvjoULF2Lfvn1ITU1FpUqVsGzZMsVSbaSfOGZMC4rzmDEiIiL6H44ZIyIiIjIwTMa0qKBrUxIREZE0pFybkpcptYCXKYmIiPQDL1MSERERGRgmY0REREQSYjJGREREJCEmY0REREQSYjJmoJKTkzF9+nScPHlS6lDoHbZs2YIlS5ZIHQYRaWD69OmQyWSIi4vT2TGbNWsGHx8fnR2vqHF3d8egQYOkDkNtTMYMVHJyMmbMmMFkrIhjMkZEpP+4HBIZvNTUVJiamsLIyDB+m8jlcqSnp8Pc3FzqUIjyJPt8r9QhAADEwo5Sh0AGwDC+fYqR6OhoDB06FC4uLjAzM4OHhwdGjhyJ9PR0xSnvt61fvx4ymQwPHz5UbLt8+TLatGkDe3t7WFhYwMPDA0OGDAEAPHz4EA4ODgCAGTNmQCaTQSaTYfr06Yr9jx8/jsaNG8PKygqlSpVC586dER4ernTcnHgiIiIwYMAA2NrawsHBAdOmTYMQApGRkejcuTNsbGzg5OSEhQsXqv16uLu7o0OHDjh8+DB8fX1hbm4Ob29v7Ny5U6neixcvMHHiRNSoUQPW1tawsbFBu3bt8NdffynVO3nyJGQyGbZt24avv/4a5cqVg6WlJRITE9VuY8eOHZgxYwbKlSuHkiVLokePHkhISEBaWhrGjx8PR0dHWFtbY/DgwUhLS1Pp26ZNm1C3bl1YWFjAzs4Offr0QWRkpKK8WbNm2L9/Px49eqR4j9zd3RXlaWlpCA4ORqVKlWBmZgY3Nzd8+eWXKseSyWQYPXo0Nm/ejOrVq8PMzAy///47AGDbtm2oW7cuSpYsCRsbG9SoUQNLly5V+30iIsPx+vVrqUPQO0zGipDHjx/D398f27ZtQ+/evfHDDz/g448/xqlTp5CcnJzvdmJjY9G6dWs8fPgQkydPxrJly9C/f3/8+eefAAAHBwesXLkSANC1a1ds3LgRGzduRLdu3QAAR48eRZs2bRAbG4vp06fjs88+w7lz5/Dhhx8qJXw5evfuDblcjnnz5qF+/fr49ttvsWTJErRq1QrlypXD/PnzUalSJUycOBF//PGH2q/L3bt30bt3b7Rr1w5z585FiRIl0LNnTxw5ckRR5/79+9i9ezc6dOiARYsW4YsvvsCNGzfQtGlTPH78WKXNWbNmYf/+/Zg4cSLmzJkDU1NTtduYO3cuDh06hMmTJ2PIkCHYuXMnRowYgSFDhiAiIgLTp09Ht27dsH79esyfP19p39mzZ2PgwIGoXLkyFi1ahPHjx+PYsWNo0qQJ4uPjAQBTp06Fr68v7O3tFe9RziVLuVyOTp064fvvv0fHjh2xbNkydOnSBYsXL0bv3r1VYj1+/DgmTJiA3r17Y+nSpXB3d8eRI0fQt29flC5dGvPnz8e8efPQrFkznD17Vu33iIjyFhcXh169esHGxgZlypTBuHHjkJqaqigPDQ1F8+bN4ejoCDMzM3h7eys+o9928OBBNG3aVPEDys/PD1u2bHnn8Q8fPgxLS0v07dsXmZmZAICUlBSMHTsW9vb2KFmyJDp16oTo6GiVH+Y5P7pv376Nfv36oXTp0mjUqBEAIDMzE7NmzYKnpyfMzMzg7u6Or776KtcfhG+2mePt8V05JxbOnj2Lzz77DA4ODrCyskLXrl3x7NkzpX2FEPj222/h6uoKS0tLfPTRR7h169Y7X4ciTVCBJSQkCAAiISGhQO0MHDhQGBkZiUuXLqmUyeVyERwcLHJ7y0JDQwUA8eDBAyGEELt27RIAcm0nx7NnzwQAERwcrFLm6+srHB0dxfPnzxXb/vrrL2FkZCQGDhyo2JYTT2BgoGJbZmamcHV1FTKZTMybN0+x/eXLl8LCwkJ88skn73oJVFSoUEEAEL/++qtiW0JCgnB2dha1a9dWbEtNTRVZWVlK+z548ECYmZmJmTNnKradOHFCABAVK1YUycnJSvXVbcPHx0ekp6crtvft21fIZDLRrl07pTY++OADUaFCBcXzhw8fCmNjYzF79mylejdu3BAlSpRQ2h4QEKC0b46NGzcKIyMjcfr0aaXtq1atEgDE2bNnFdsACCMjI3Hr1i2luuPGjRM2NjYiMzNTpX0iqeGzPUXiURA5n5E1atQQHTt2FMuXLxcDBgwQAMTHH3+sqOfn5ycGDRokFi9eLJYtWyZat24tAIjly5crtRcaGipkMpnw8fERs2fPFiEhIWLYsGFKbTVt2lRUr15d8Xzv3r3CzMxMDBw4UOlvvVevXoo4QkJCRK9evUStWrVUvhdy+uDt7S06d+4sVqxYIUJCQoQQQnzyyScCgOjRo4cICQkRAwcOFABEly5dlN/LPL5rKlSooPSdkPNdVrt2bdG8eXOxbNky8fnnnwtjY2PRq1cvpX2//vprAUC0b99eLF++XAwZMkS4uLgIe3t7tb9n3qat73R1cMyYFvn5+cHY2BhBQUEICgpSa1+5XI7du3ejY8eOqFevnkp5bpcn81KqVCkAwL59+1CrVi2YmJjke9+YmBhcv34dX375Jezs7BTba9asiVatWuHAgQMq+wwbNkzx/8bGxqhXrx6ioqIwdOhQpZiqVKmC+/fv5zuWHC4uLujataviuY2NDQYOHIj58+fjyZMncHJygpmZmaI8KysL8fHxsLa2RpUqVXD16lWVNj/55BNYWFgobVO3jYEDByq9tvXr18fWrVsVl4Pf3P7DDz8gMzMTJUqUwM6dOyGXy9GrVy+lu6ycnJxQuXJlnDhxAl999dU7X5NffvkF1apVQ9WqVZXaaN68OQDgxIkTaNiwoWJ706ZN4e3trdRGqVKl8Pr1axw5cgRt27Z95/GISHMeHh747bffAABBQUGwsbHBihUrMHHiRNSsWROnTp1S+jwaPXo02rZti0WLFim+SxISEjB27Fj4+/vj5MmTSmM+RR6rGu7cuRN9+vTBoEGDsGrVKsW42KtXr2LHjh0YP348Fi9eDAAYNWoUBg8erDIsI0etWrWUzsD99ddf2LBhA4YNG4Y1a9Yo2nB0dMT333+PEydO4KOPPtLo9SpTpgwOHz6s+N6Ty+X44YcfkJCQAFtbWzx79gzfffcdAgICsHfvXkW9qVOnYs6cORodE8hemzIkJARZWVkat6EpXqbUokuXLuH27dtqJ2IA8OzZMyQmJmrlluSmTZuie/fumDFjBuzt7dG5c2eEhobmOm7pbY8ePQIAVKlSRaWsWrVqiIuLUxkvUL58eaXntra2MDc3h729vcr2ly9fqtsdVKpUSSUZ9fLyAgDFZVO5XI7FixejcuXKMDMzg729PRwcHPD3338jISFBpU0PDw+Vbeq2kVu/AcDNzU1lu1wuV7Rx9+5dCCFQuXJlODg4KD3Cw8MRGxv73tfk7t27uHXrlsr+Oa/L223k1t9Ro0bBy8sL7dq1g6urK4YMGaIYS0ZE2vP2d8KYMWMAQPHj9s1ELCEhAXFxcWjatCnu37+v+Nw4cuQIXr16hcmTJ6vcfJPbj/WtW7eid+/e+PTTT/Hjjz8q3aCU83c+atSoXOPKzYgRI5Se58T+2WefKW3//PPPAQD79+/Ps633CQwMVOpT48aNkZWVpfh+Onr0KNLT0zFmzBileuPHj9f4mED2+3T79m1cunSpQO1ogmfGipG8zo69ncXLZDKEhYXhzz//xN69e3Ho0CEMGTIECxcuxJ9//glra2utxmVsbJyvbUDev+AKas6cOZg2bRqGDBmCWbNmwc7ODkZGRhg/fjzkcrlK/bfPimnSRl59fF/f5XI5ZDIZDh48mGvd/Lw/crkcNWrUwKJFi3ItfzshzK2/jo6OuH79Og4dOoSDBw/i4MGDCA0NxcCBA7Fhw4b3xkBE+VO5cmWl556enjAyMlL8mDx79iyCg4Nx/vx5lfHBOWeD7t27BwD5+sH+4MEDDBgwAD179sSyZctUyh89egQjIyOVH2mVKlXKs8236+a08fY+Tk5OKFWqlCJx0sTbP3RLly4NAIof8zltv/26Ojg4KOoWN0zGiggHBwfY2Njg5s2bedbJ+UcWHx+vuBQJIM9/9A0aNECDBg0we/ZsbNmyBf3798e2bdswbNiwPBO7ChUqAAD++ecflbI7d+7A3t4eVlZW+e2WVvz7778QQijFHBERAQCKuwvDwsLw0UcfYd26dUr7xsfHq5yhy4s22sgPT09PCCHg4eGhOJOVl7zeJ09PT/z1119o0aKFWpew32ZqaoqOHTuiY8eOkMvlGDVqFH788UdMmzbtnR/MRKS5N/9m7927hxYtWqBq1apYtGgR3NzcYGpqigMHDmDx4sW5/hB8H2dnZzg7O+PAgQO4fPlyrkNf1JXbDzpAvSE0b8vrcqCuf8wXBbxMWUQYGRmhS5cu2Lt3Ly5fvqxSLoSAp6cnACjdkfj69WuVsxgvX75U+Ufr6+sLAIpLlZaWlgCguHMvh7OzM3x9fbFhwwalsps3b+Lw4cNo3769Rv0riMePH2PXrl2K54mJifj555/h6+sLJycnANl/vG/3+ZdffkF0dHS+j6ONNvKjW7duMDY2xowZM1SOJ4TA8+fPFc+trKxyvUTaq1cvREdHK8ZqvCklJSVft56/eRwg+99gzZo1ASBfl7SJKH/u3r2r9Pzff/+FXC6Hu7s79u7di7S0NOzZsweffvop2rdvj5YtW6okPzmf/+/6wZ7D3Nwc+/btQ+XKldG2bVuVuwwrVKgAuVyOBw8eqMSVXzltvN23p0+fIj4+XvHDHsg+kfD2d016ejpiYmLyfby3jw2ovq7Pnj3TaChMUcAzY0XInDlzcPjwYTRt2hSBgYGoVq0aYmJi8Msvv+DMmTNo3bo1ypcvj6FDh+KLL76AsbExfvrpJzg4OOC///5TtLNhwwasWLECXbt2haenJ169eoU1a9bAxsZGkUxZWFjA29sb27dvh5eXF+zs7ODj4wMfHx8sWLAA7dq1wwcffIChQ4ciJSUFy5Ytg62tba63Jxc2Ly8vDB06FJcuXULZsmXx008/4enTpwgNDVXU6dChA2bOnInBgwejYcOGuHHjBjZv3oyKFSvm+zjaaCM/PD098e2332LKlCl4+PAhunTpgpIlS+LBgwfYtWsXAgMDMXHiRABA3bp1sX37dnz22Wfw8/ODtbU1OnbsiI8//hg7duzAiBEjcOLECXz44YfIysrCnTt3sGPHDhw6dOi9v4aHDRuGFy9eoHnz5nB1dcWjR4+wbNky+Pr6olq1alrtM5EhCwkJQevWrRXPcy4dtmvXTvHj+s0fZgkJCUqfbwDQunVrlCxZEnPnzkXbtm1VBvC/fYbK1tYWhw4dQpMmTdCqVSucPn1akdC1adMGU6dOxYoVKxQD+N+MKz/at2+Pr776CkuWLMGPP/6o2J4zdCIgIECxzdPTU2Vao9WrV2s8UL5ly5YwMTHBsmXL0Lp1a0Xfi/NqJcUiGUtKSsKCBQtw4cIFXLx4ES9fvkRoaOg715/KyMhArVq1EB4ejgULFii+3HLI5XJ8//33WLlyJWJiYuDl5YUpU6agb9++hdybvJUrVw4XLlzAtGnTsHnzZiQmJqJcuXJo164dLC0tYWJigl27dmHUqFGYNm0anJycMH78eJQuXRqDBw9WtNO0aVNcvHgR27Ztw9OnT2Frawt/f39s3rxZ6br/2rVrMWbMGEyYMAHp6ekIDg6Gj48PWrZsid9//x3BwcH45ptvYGJigqZNm2L+/Pm5DgQvbJUrV8ayZcvwxRdf4J9//oGHhwe2b9+ONm3aKOp89dVXeP36NbZs2YLt27ejTp062L9/PyZPnpzv42ijjfyaPHkyvLy8sHjxYsyYMQNA9jiv1q1bo1OnTop6o0aNwvXr1xEaGorFixejQoUK6NixI4yMjLB7924sXrwYP//8M3bt2gVLS0tUrFgR48aNe+/lTwAYMGAAVq9ejRUrViA+Ph5OTk7o3bs3pk+fbjCrERDpwoMHD9CpUye0bdsW58+fx6ZNm9CvXz/UqlUL5ubmiuECn376KZKSkrBmzRo4OjoqnTmysbHB4sWLMWzYMPj5+Snm/Prrr7+QnJyc6zhPe3t7HDlyBI0aNULLli1x5swZlCtXDnXr1kX37t2xZMkSPH/+HA0aNMCpU6cUwz/yc+mxVq1a+OSTT7B69WrEx8crvnc2bNiALl26KN1JOWzYMIwYMQLdu3dHq1at8Ndff+HQoUMaD/9wcHDAxIkTMXfuXHTo0AHt27fHtWvXcPDgQa0OKdEpnU2iUQAPHjwQAET58uVFs2bNBAARGhr6zn0WLlworKysBACxYMEClfLJkycLAGL48OFi9erVIiAgQAAQW7duVTs+KeYkMRQVKlQQAQEBUodBRDom9fxi2pxn7Pbt26JHjx6iZMmSonTp0mL06NEiJSVFUW/Pnj2iZs2awtzcXLi7u4v58+eLn376SWn+yDfrNmzYUFhYWAgbGxvh7++v9L319jxjQgjx77//CmdnZ1GtWjXx7NkzIYQQr1+/FkFBQcLOzk5YW1uLLl26iH/++UcAUJojMqcPOfu9KSMjQ8yYMUN4eHgIExMT4ebmJqZMmSJSU1OV6mVlZYlJkyYJe3t7YWlpKdq0aSP+/fffPOcZe3uOzJy5HU+cOKHU5owZM4Szs7OwsLAQzZo1Ezdv3lRpUxNSfKcXi2QsNTVVxMTECCGEuHTp0nuTsadPnwpbW1sxc+bMXJOxqKgoYWJiIoKCghTb5HK5aNy4sXB1dVV7EkwmY4WHyRgRkW5cu3ZNABCbNm2SOhRJSfGdXiyuRZiZmSkGaufH5MmTUaVKFQwYMCDX8t9++w0ZGRlKc6zIZDKMHDkSUVFROH/+fIFjpnd79uwZnjx5kufjxYsXUodIRKS3UlJSVLYtWbIERkZGaNKkiQQRGbZiMWZMHTnXrM+cOZPnde9r167ByspKZZCyv7+/ojxn7S11JCYmKj03MzNTmtWd/sfPz++d89A0bdoUJ0+e1F1AREQG5LvvvsOVK1fw0UcfoUSJEoq5BgMDA1XmKdR3aWlpSneQv/1drgt6lYwJITBmzBj07t0bH3zwQa6LWgPZS/6ULVtWJVlzdnYGgFwXhc6Pt/8BBwcHS3L3YXGwefPmXH+Z5ciZUy2v95CIiDTXsGFDHDlyBLNmzUJSUhLKly+P6dOnY+rUqVKHpnNz585V3EglFb1KxtavX48bN24gLCzsnfVSUlJyPWOVc6vwu5KEd4mMjISNjY3iOc+K5e3DDz+UOgQiIoPVqlUrtGrVSuowioQpU6YoLeuUmJio87ODepOMJSYmYsqUKfjiiy/e+yJaWFjkOqllamqqolwTNjY2SskYERERFW1FYUhRsRjAnx/ff/890tPT0bt3bzx8+BAPHz5EVFQUgOwZ6R8+fIj09HQA2Zcjnzx5ojL7ec6cLi4uLroNnoiIiAyW3iRj//33H16+fInq1avDw8MDHh4eaNy4MYDsme09PDxw+/ZtANlLAyUnJyM8PFypjQsXLijKiYiIiHRBby5Tjh07Fl26dFHaFhsbi08//RSDBg1C586dFbPHd+7cGRMmTMCKFSuwfPlyANmD/1etWoVy5cqhYcOGug6fiIiIDFSxScaWL1+O+Ph4xZ2Oe/fuVVyGHDNmDOrUqYM6deoo7ZNzJ1716tWVEjVXV1eMHz8eCxYsQEZGBvz8/LB7926cPn0amzdvznPFeCIiIiJtKzbJ2Pfff680L9XOnTuxc+dOANlr7Nna2qrV3rx581C6dGn8+OOPWL9+PSpXrqxYL4yIiIhIV2Ti7VHspLbExETY2toiISGBd1MSEREVY1J8p+vNAH4iIiKi4ojJGBERERWKZs2aoVmzZlKHUeQxGSMiomLl2rVrmDFjBvbt26cyX+TbIiIi8O2332L79u3Iysp6Z93o6GjMmzcPP/30U64Tg6tr+vTpkMlkiIuLy7Xcx8enSCcq586dw/Tp0xEfHy/J8bds2YIlS5ZIcmxdYzKmRX5+fvD29kZISIjUoRAR6aVr165hz549qFu3LgICAlTWGH5TREQEduzYgcqVK6NHjx7vvFM+OjoaGzduhKOjI/r37y/5jOxFwblz5zBjxowCJWOHDx/G4cOHNdpX18lYSEgIvL294efnp7Nj5ig2d1MWB5cuXeIAfiKiQsJETDdev34NKysrrbRlamqqlXZ0ISgoCEFBQYoB/LrEM2NERFTk6XsidvLkSchkMuzYsQOzZ8+Gq6srzM3N0aJFC/z7778q9S9cuID27dujdOnSsLKyQs2aNbF06VKlOnfu3EGPHj1gZ2cHc3Nz1KtXD3v27FGqs379eshkMpw6dQqjRo2Co6MjXF1dMX36dHzxxRcAAA8PD8hkMshkMsX8naGhoWjevDkcHR1hZmYGb29vrFy5UiXOt8eM5befzZo1w/79+/Ho0SPFsd3d3ZGUlAQrKyuMGzdO5VhRUVEwNjbG3Llz8/26FxU8M0ZEREWavidib5o3bx6MjIwwceJEJCQk4LvvvkP//v0Vy/UBwJEjR9ChQwc4Oztj3LhxcHJyQnh4OPbt26dIUm7duoUPP/wQ5cqVw+TJk2FlZYUdO3agS5cu+PXXX9G1a1el444aNQoODg745ptv8Pr1a7Rr1w4RERHYunUrFi9eDHt7ewCAg4MDAGDlypWoXr06OnXqhBIlSmDv3r0YNWoU5HI5goKCCtzPqVOnIiEhAVFRUVi8eDEAwNraGtbW1ujatSu2b9+ORYsWKb2/W7duhRAC/fv3L8A7IBFBBZaQkCAAiISEBKlDISLSK1evXhXTp08Xe/fuFXK5/J11//nnHzFr1iyxbds2kZmZ+c66UVFRYu7cuWLdunUiNTVVmyErBAcHCwDi2bNnuZZXr15dNG3aVAghxIkTJwQAUa1aNZGWlqaos3TpUgFA3LhxQwghRGZmpvDw8BAVKlQQL1++VGrvzdenRYsWokaNGkp9k8vlomHDhqJy5cqKbaGhoQKAaNSokcprtmDBAgFAPHjwQCX25ORklW1t2rQRFStWVNrWtGlTRR/V6acQQgQEBIgKFSqoHOfQoUMCgDh48KDS9po1ayodS1NSfKfzMiURERVJhnRGLMfgwYOVxlk1btwYAHD//n0A2a/JgwcPMH78eJQqVUpp35zX58WLFzh+/Dh69eqFV69eIS4uDnFxcXj+/DnatGmDu3fvIjo6Wmnf4cOHq7UUoIWFheL/ExISEBcXh6ZNm+L+/ftISEgocD/fpWXLlnBxccHmzZsV227evIm///4bAwYMyHcfihJepiQioiKnKCRiz58/R5kyZQrUj/d5u1/ly5dXel66dGkAwMuXLwEA9+7dA5A9LUZe/v33XwghMG3aNEybNi3XOrGxsShXrpziuYeHh1pxnz17FsHBwTh//jySk5OVyhISEt47AP59/XwXIyMj9O/fHytXrkRycjIsLS2xefNmmJubo2fPnmr1o6hgMkZEREVOUUjENmzYgM8++0zjPpibmwMAUlJSci1PTk5W1MmRV/xCjZUL5XI5AGDixIlo06ZNrnUqVaqk9PzNM13vc+/ePbRo0QJVq1bFokWL4ObmBlNTUxw4cACLFy9WHP9dCtrPgQMHYsGCBdi9ezf69u2LLVu2oEOHDjq/C1JbmIwREVGRUxQSsYJevqxQoQIA4J9//oGbm5tSWXJyMiIjI9G6dWu12vT09ASQfVmuZcuWudapWLEiAMDExCTPOvmR12u/d+9epKWlYc+ePUpnuE6cOKHxsdQ5PpB9ZrB27drYvHkzXF1d8d9//2HZsmVaPb4uccwYEREVOUUhEfvkk08K1IcWLVrA1NQUK1euVDlbtHr1amRmZqJdu3ZqtVmnTh14eHhgyZIlKpOx5pxVcnR0RLNmzfDjjz8iJiZGpY1nz57l61g5c429fZyc1/rNs1gJCQkIDQ3Nbzfyffx3jT/7+OOPcfjwYSxZsgRlypRR+7UsSnhmjIiIipyikIhZW1sXqA+Ojo745ptv8PXXX6NJkybo1KkTLC0tce7cOWzduhWtW7dGx44d1WrTyMgIK1euRMeOHeHr64vBgwfD2dkZd+7cwa1bt3Do0CEA2bPJN2rUCDVq1MDw4cNRsWJFPH36FOfPn0dUVBT++uuv9x6rbt26ALKnmejTpw9MTEzQsWNHtG7dGqampujYsSM+/fRTJCUlYc2aNXB0dMw1+dNU3bp1sX37dnz22Wfw8/ODtbW10uvVr18/fPnll9i1axdGjhwJExMTrR1b53R236Ye49QWRESUl02bNokGDRoIKysrYWZmJqpWrSpmzJihNO1EzpQPv/zyi9K+Dx48EABEaGio0vYzZ86IVq1aiZIlSworKytRs2ZNsWzZMqU69+7dEwMHDhROTk7CxMRElCtXTnTo0EGEhYUp6uRMbXHp0qVcY581a5YoV66cMDIyUprmYs+ePaJmzZrC3NxcuLu7i/nz54uffvpJZSqMvKa2yE8/k5KSRL9+/USpUqUEgFynuWjfvr0AIM6dO5dr/JqQ4jtdJoQaowIpVzlLJ3h5ecHY2FixpAIREREVnq5du+LGjRu5rlKgrpCQEISEhCArKwsRERFISEjQ2RKHTMa0ICcZ0+UbR0REZMhiYmJQoUIFTJ06FcHBwVprV4rvdI4ZIyIiomLjwYMHOHv2LNauXQsTExN8+umnUodUYLybkoiIiIqNU6dO4eOPP8aDBw+wYcMGODk5SR1SgfEypRbwMiUREZF+kOI7nWfGiIiIiCTEZIyIiIhIQkzGiIiIiCTEZIyIiIhIQkzGiIiIiCTEZIyIiIhIQkzGiIhIr7x+/RoymQwymQyvX7+WOhyi92IypkV+fn7w9vZGSEiI1KEQERGRGkJCQuDt7Q0/Pz+dH5uTvmoBJ30lIio6Xr9+DWtrawBAUlISrKysJIlj/fr1GDx4sOK5mZkZ7OzsUKNGDQQEBGDw4MEoWbKk2u2eO3cOhw8fxvjx41GqVCktRkwAJ30lIiLSOzNnzsTGjRuxcuVKjBkzBgAwfvx41KhRA3///bfa7Z07dw4zZsxAfHy8liMlqXChcCIi0itZWVmK///jjz/QunVrGBsbSxZPu3btUK9ePcXzKVOm4Pjx4+jQoQM6deqE8PBwWFhYSBYfSY9nxoiISG/s3LkT3t7eiuft27eHu7s7du7cKWFUqpo3b45p06bh0aNH2LRpEwDg77//xqBBg1CxYkWYm5vDyckJQ4YMwfPnzxX7TZ8+HV988QUAwMPDQ3GjwsOHDwEAoaGhaN68ORwdHWFmZgZvb2+sXLlS5/0j9fDMGBER6YWdO3eiR48eeHsodHR0NHr06IGwsDB069ZNouhUffzxx/jqq69w+PBhDB8+HEeOHMH9+/cxePBgODk54datW1i9ejVu3bqFP//8EzKZDN26dUNERAS2bt2KxYsXw97eHgDg4OAAAFi5ciWqV6+OTp06oUSJEti7dy9GjRoFuVyOoKAgKbtL7yKKgVevXolvvvlGtGnTRpQuXVoAEKGhoUp1srKyRGhoqOjYsaNwdXUVlpaWonr16mLWrFkiJSUl13bXrl0rqlatKszMzESlSpXEDz/8oFF8CQkJAoBISEjQaH8iIiqYzMxM4erqKgDk+pDJZMLNzU1kZmbqLKbQ0FABQFy6dCnPOra2tqJ27dpCCCGSk5NVyrdu3SoAiD/++EOxbcGCBQKAePDggUr93Npo06aNqFixogY9MExSfKcXi8uUcXFxmDlzJsLDw1GrVq1c6yQnJ2Pw4MF49uwZRowYgSVLlsDf3x/BwcFo166dyi+lH3/8EcOGDUP16tWxbNkyfPDBBxg7dizmz5+viy4REZEWnT59GlFRUXmWCyEQGRmJ06dP6zCq97O2tsarV68AQGncWGpqKuLi4tCgQQMAwNWrV/PV3pttJCQkIC4uDk2bNsX9+/eRkJCgxchJm4rFZUpnZ2fExMTAyckJly9fznUOEFNTU5w9exYNGzZUbBs+fDjc3d0RHByMY8eOoWXLlgCAlJQUTJ06FQEBAQgLC1PUlcvlmDVrFgIDA1G6dGnddI6IiAosJiZGq/V0JSkpCY6OjgCAFy9eYMaMGdi2bRtiY2OV6uU3kTp79iyCg4Nx/vx5JCcnq7Rha2urncBJq4rFmTEzMzM4OTm9s46pqalSIpaja9euAIDw8HDFthMnTuD58+cYNWqUUt2goCC8fv0a+/fv10LURESkK87OzlqtpwtRUVFISEhApUqVAAC9evXCmjVrMGLECOzcuROHDx/G77//DgCQy+Xvbe/evXto0aIF4uLisGjRIuzfvx9HjhzBhAkT8t0GSaNYnBkriCdPngCAYpAjAFy7dg0AlG41BoC6devCyMgI165dw4ABA3QXJBERFUjjxo3h6uqK6OholWEpACCTyeDq6orGjRtLEF3uNm7cCABo06YNXr58iWPHjmHGjBn45ptvFHXu3r2rsp9MJsu1vb179yItLQ179uxB+fLlFdtPnDih5chJ24rFmbGC+O6772BjY4N27doptsXExMDY2FhxajiHqakpypQpg8ePH2t0rMTERKVHWlpagWInIqL8MTY2xtKlSwGoJis5z5csWSLpfGNvOn78OGbNmgUPDw/0799fEdfbieSSJUtU9s1ZUeDtSV9zayMhIQGhoaFajFz/pKWlqXx/65penxmbM2cOjh49ihUrVigtGZGSkgJTU9Nc9zE3N0dKSopGx3Nzc1N6HhwcjOnTp2vUFhERqadbt24ICwvD2LFjER0drdju6uqKJUuWSDatxcGDB3Hnzh1kZmbi6dOnOH78OI4cOYIKFSpgz549MDc3h7m5OZo0aYLvvvsOGRkZKFeuHA4fPowHDx6otFe3bl0AwNSpU9GnTx+YmJigY8eOaN26NUxNTdGxY0d8+umnSEpKwpo1a+Do6FjkxsoVJXPnzsWMGTOkDUJn921qyaVLl3Kd2uJt27ZtEzKZTAwdOlSlLCgoSBgbG+e6n4ODg+jTp49aMeXcBhsZGSkSEhIUj9TUVLXaISKigsv5TAYgDhw4oNPpLN6UM7VFzsPU1FQ4OTmJVq1aiaVLl4rExESl+lFRUaJr166iVKlSwtbWVvTs2VM8fvxYABDBwcFKdWfNmiXKlSsnjIyMlKa52LNnj6hZs6YwNzcX7u7uYv78+eKnn37KcyoMEiI1NVXpuzsyMlLnU1sUu4XCc+6mDA0NxaBBg3Ktc+TIEXTo0AGtW7fGrl27UKKE8gnA2bNn4+uvv8bTp0+VLlWmp6fDwsIC48ePx8KFC/MdExcKJyIqOorKQuFUPHGhcC24cOECunbtinr16mHHjh0qiRgA+Pr6AshO7N50+fJlyOVyRTkRERFRYdOrZCw8PBwBAQFwd3fHvn378lx4tXnz5rCzs1NZr2vlypWwtLREQECALsIlIqJCYGVlBSEEhBA8K0bFQrEZwL98+XLEx8cr7nTcu3evYrblMWPGwMjISHF78BdffKEyV5inpyc++OADANkzFM+aNQtBQUHo2bMn2rRpg9OnT2PTpk2YPXs27OzsdNs5IiIiMljFZsyYu7s7Hj16lGtZzt0mHh4eee7/ySefYP369Urb1qxZg4ULF+LBgwdwc3PD6NGjMW7cuDzncMkLx4wRERHpBym+04tNMlaUMRkjIiLSDxzAT0RERGRgmIwRERERSYjJGBEREZGEmIxpkZ+fH7y9vRESEiJ1KERERKSGkJAQeHt7w8/PT+fH5gB+LeAAfiIiIv3AAfxEREREBobJGBEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYjJGBEREZGESkgdgD7x8/ODsbExgoKCEBQUJHU4RKRlMv8qQEam1GEQUWF4lgjEJQASrBLJZEyLLl26xLUpifRZRiYwtqvUURBRYUpJA4KW6/SQvExJREREJCEmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQS4nJIWsS1KYmIiIqpY9eA49cBOdemLNa4NiUREVEx1aJ29oNrUxIREREZFiZjRERERBJiMkZEREQkISZjRERERBJiMkZEREQkoWKRjCUlJSE4OBht27aFnZ0dZDIZ1q9fn2vd8PBwtG3bFtbW1rCzs8PHH3+MZ8+eqdSTy+X47rvv4OHhAXNzc9SsWRNbt24t5J4QERERKSsWyVhcXBxmzpyJ8PBw1KpVK896UVFRaNKkCf7991/MmTMHEydOxP79+9GqVSukp6cr1Z06dSomTZqEVq1aYdmyZShfvjz69euHbdu2FXZ3iIiIiBSKxTxjzs7OiImJgZOTEy5fvgw/P79c682ZMwevX7/GlStXUL58eQCAv78/WrVqhfXr1yMwMBAAEB0djYULFyIoKAjLl2fPJTJs2DA0bdoUX3zxBXr27AljY2PddI6IiIgMWrE4M2ZmZgYnJ6f31vv111/RoUMHRSIGAC1btoSXlxd27Nih2Pbbb78hIyMDo0aNUmyTyWQYOXIkoqKicP78ee12gIiIiCgPxSIZy4/o6GjExsaiXr16KmX+/v64du2a4vm1a9dgZWWFatWqqdTLKddEYmKi0iMtLU2jdoiIiEhHMjKzZ91XPNLfv4+WFYvLlPkRExMDIPuS5tucnZ3x4sULpKWlwczMDDExMShbtixkMplKPQB4/PixRjG4ubkpPQ8ODsb06dM1aouIiIh0YP9FYI+0V8T0JhlLSUkBkH1J823m5uaKOmZmZor/vqueJiIjI5XWpsztGERERFSEBPgDber+73lKOjBxtU5D0JtkzMLCAgByvTSYmpqqVMfCwiJf9dRlY2PDhcKJiIiKE5MS2Q8J6c2YsZxLjDmXK98UExMDOzs7xZkqZ2dnPHnyBEIIlXoA4OLiUsjREhEREWXTm2SsXLlycHBwwOXLl1XKLl68CF9fX8VzX19fJCcnIzw8XKnehQsXFOVEREREuqA3yRgAdO/eHfv27UNkZKRi27FjxxAREYGePXsqtnXu3BkmJiZYsWKFYpsQAqtWrUK5cuXQsGFDncZNREREhqvYjBlbvnw54uPjFXc67t27F1FRUQCAMWPGwNbWFl999RV++eUXfPTRRxg3bhySkpKwYMEC1KhRA4MHD1a05erqivHjx2PBggXIyMiAn58fdu/ejdOnT2Pz5s2c8JWIiIh0RibeHjhVRLm7u+PRo0e5lj148ADu7u4AgFu3buGzzz7DmTNnYGpqioCAACxcuBBly5ZV2kcul2P+/Pn48ccfERMTg8qVK2PKlCno37+/2rElJibC1tYWCQkJHMBPpMdktT2BsV2lDoOIClNKGhC0XKff6cUmGSvKmIwRGQYmY0QGQIJkTK/GjBEREREVN0zGiIiIiCTEZIyIiIhIQsXmbsriwM/PD8bGxggKCkJQUJDU4RAREVF+HbsGHL8OyHU/lJ7JmBZdunSJA/iJiIiKoxa1sx//P4Bfl3iZkoiIiEhCTMaIiIiIJKSTZOz+/fsq60ASERERkZaTsR9++AF9+vRR2jZ48GBUrlwZPj4+qFevHmJjY7V5SCIiIqJiTavJ2Nq1a5WWHTp06BA2bNiAwMBALFu2DPfv38eMGTO0eUgiIiKiYk2rd1M+evQI1apVUzzfsWMHPDw8sHLlSgDAkydPsHHjRm0ekoiIiKhY0+qZsbeXuTx8+DDatWuneO7u7o4nT55o85BERERExZpWkzEvLy/s2rULQPYlysePHyslY1FRUShVqpQ2D0lERERUrGn1MuXEiRPRr18/lC5dGq9fv0a1atXQpk0bRfnx48fh6+urzUMSERERFWtaTcb69OmDMmXK4MCBAyhVqhRGjRqFEiWyD/HixQvY2dlhwIAB2jxkkcLlkIiIiIopCZdDkom3B3qR2hITE2Fra4uEhAQuh0Skx2S1PYGxXaUOg4gK0/8vh6TL73StjhkzNjbGli1b8izfvn07jI2NtXlIIiIiomKtUO+mfFtWVhZkMpk2D0lERERUrGl9OaS8kq3ExEQcOnQI9vb22j4kERERUbFV4GRsxowZMDY2hrGxMWQyGQYMGKB4/uajdOnS2Lhxo8pySURERESGrMB3U/r7+2PUqFEQQmDFihVo1aoVvLy8lOrIZDJYWVmhbt266NatW0EPSURERKQ3CpyMtWvXTjGx6+vXrzFixAjUr1+/wIERERERGQKtzjMWGhqqzeaIiIiI9J5WkzEg+47JQ4cO4f79+3j58qXKHZYymQzTpk3T9mGJiIiIiiWtJmOXL19G9+7dERUVlec0F0zGiIiIiP5Hq1NbjBo1CikpKdi9ezdevHgBuVyu8sjKytLmIYmIiIiKNa2eGfv7778xe/ZsdOzYUZvNFhtcm5KIiKiYknBtSq0mY66uru+dhV+fXbp0iWtTEhERFUctamc//n9tSl3S6mXKSZMmYc2aNUhMTNRms0RERER6S6tnxl69egVra2tUqlQJffr0gZubm8rC4DKZDBMmTNDmYYmIiIiKLZnQ4nVFI6P3n2iTyWR6N4g/MTERtra2SEhI4GVKIj0mq+0JjO0qdRhEVJj+/zKlLr/TtXpm7MGDB9psjoiIiEjvaTUZq1ChgjabIyIiItJ7Wh3AnyM6Ohpbt27F0qVLERUVBSB7Zv4XL14U+iXKu3fvok+fPnB1dYWlpSWqVq2KmTNnIjk5WaneuXPn0KhRI1haWsLJyQljx45FUlJSocZGRERE9DatnhkTQuDzzz/H8uXLkZmZCZlMhho1asDV1RVJSUlwd3fHzJkzMX78eG0eViEyMhL+/v6wtbXF6NGjYWdnh/PnzyM4OBhXrlzBb7/9BgC4fv06WrRogWrVqmHRokWIiorC999/j7t37+LgwYOFEhsRERFRbrSajC1YsABLly7FpEmT0KJFC7Rq1UpRZmtri27duuHXX38ttGRs48aNiI+Px5kzZ1C9enUAQGBgIORyOX7++We8fPkSpUuXxldffYXSpUvj5MmTisF57u7uGD58OA4fPozWrVsXSnxEREREb9PqZco1a9Zg4MCBmDNnDnx9fVXKa9asiYiICG0eUknO/GZly5ZV2u7s7AwjIyOYmpoiMTERR44cwYABA5Tukhg4cCCsra2xY8eOQouPiIiI6G1aTcYiIyPRsGHDPMutrKwKdULYZs2aAQCGDh2K69evIzIyEtu3b8fKlSsxduxYWFlZ4caNG8jMzES9evWU9jU1NYWvry+uXbum8fETExOVHmlpaQXpDhERERW2jMzs6SwUj3Sdh6DVy5SOjo6IjIzMs/zKlSsoX768Ng+ppG3btpg1axbmzJmDPXv2KLZPnToV3377LQAgJiYGQPbZsrc5Ozvj9OnTGh/fzc1N6XlwcDCmT5+ucXtERERUyPZfBPaclzQErSZj3bp1w6pVqzBo0CDY2toCyJ7kFQAOHz6M9evX48svv9TmIVW4u7ujSZMm6N69O8qUKYP9+/djzpw5cHJywujRo5GSkgIAMDMzU9nX3NxcUa6JyMhIpUufuR2DiIiIipAAf6BN3f89T0kHJq7WaQhaTcZmzJiBEydOwNfXF40bN4ZMJsP8+fMxbdo0nD9/HrVr18ZXX32lzUMq2bZtGwIDAxEREQFXV1cA2QmiXC7HpEmT0LdvX1hYWABArpcQU1NTFeWasLGx4Qz8RERExYlJieyHhLQ6ZszW1hZ//vknvvzyS0RHR8Pc3BynTp1CfHw8goODcfr0aVhaWmrzkEpWrFiB2rVrKxKxHJ06dUJycjKuXbumuDyZc7nyTTExMXBxcSm0+IiIiIjepvVJXy0sLPD111/j+vXreP36NVJSUnDz5k188803BTrrlB9Pnz7NdVLZjIwMAEBmZiZ8fHxQokQJXL58WalOeno6rl+/nutdoERERESFRavJ2IoVK/Ds2TNtNqkWLy8vXLt2TWX6jK1bt8LIyAg1a9aEra0tWrZsiU2bNuHVq1eKOhs3bkRSUhJ69uyp67CJiIjIgMmEEEJbjRkZGaFEiRJo2rQp+vTpg65du8LOzk5bzb/XH3/8gebNm6NMmTIYPXo0ypQpg3379uHgwYMYNmwY1qxZAwC4evUqGjZsCG9vbwQGBiIqKgoLFy5EkyZNcOjQIbWPm5iYCFtbW52u8E5Euier7QmM7Sp1GERUmFLSgKDlOv1O1+qZsTt37uDrr79GTEwMhg8fDmdnZ7Rv3x4bN24s1PnFcjRp0gTnzp1D3bp1sWLFCowfPx737t3D7NmzsXLlSkW9OnXq4OjRo7CwsMCECROwevVqDB06FGFhYYUeIxEREdGbtHpm7E23bt3Ctm3b8MsvvyAiIgJmZmZo06YN+vTpgz59+hTGISXDM2NEhoFnxogMgARnxgotGXvTX3/9hW3btmHFihV4/fo1MjMzC/uQOsVkjMgwMBkjMgDF/TJlbv7++2/s2LEDYWFhePXqFSdCJSIiInpDoSRjt2/fRnBwMKpVq4batWtj4cKF8Pb2xqZNm/D06dPCOCQRERFRsaTVKWdnzZqFHTt24Pbt2zA2NkaLFi0wefJkdOnSRbE8kj6zdSgDyGSAvS3gwMuVRHpH4lm6iagQHbsGHL8OyAt99JYKrY4ZMzExQdOmTdG7d29069YNZcqU0VbTRVrOmDGEjAYseBmWiIio2JJgzJhWf+ZFR0fD0dERaWlpuHr1KmJjY/Hhhx/C3t5em4chIiIi0htaHTPm6OiIH374Ac7OzmjUqBG6deuGv//+GwAQFxcHe3t7/PTTT9o8JBEREVGxptVkLDQ0FOPHj0fbtm2xbt06vHkF1N7eHs2bN8e2bdu0eUgiIiKiYk2rydjChQvRuXNnbNmyBR07dlQpr1u3Lm7duqXNQxIRERWOtAxgyMLsR1qG1NGQHtNqMvbvv/+iXbt2eZbb2dnh+fPn2jwkERERUbGm1WSsVKlSiIuLy7P89u3bcHJy0uYhiYiIiIo1rSZj7du3x+rVqxEfH69SduvWLaxZswadOnXS5iGJiIiIijWtJmPffvstsrKy4OPjg6+//hoymQwbNmzAgAEDUK9ePTg6OuKbb77R5iGJiIiIijWtJmMuLi64cuUK2rZti+3bt0MIgY0bN2Lv3r3o27cv/vzzT845RkRERPQGra/t4ejoiLVr12Lt2rV49uwZ5HI5HBwcYGRU6GuSExERERU7hbrQmoODQ2E2X/TM3AwYyYDmvkCL2lJHQ0RERPkl4dqUXPVWm77pz7UpiYiIiqMWtbMf/782pS7x2iERERGRhJiMEREREUmIyRgRERGRhJiMEREREUmIyRgRERGRhJiMERER5UYu/9//R0QpPyfSIk5tQURE9LYrd4HNx//3fPFOoLQ10K85ULeydHGRXuKZMSIiojdduQuE7AHik5S3v0zK3n7lrjRxkd5iMkZERJRDLge2HH93na0neMmStIrJGBERUY6I6OwzYO/y4lV2PSIt4ZgxbeLalERExdvblyYLWo+KD65NqSe4NiURUfFWylq79aj44NqURERERYBXuey7Jt/FrmR2PSItYTJGRESUw8goe/qKd+n7UXY9Ii3hvyYiIqI31a0MBHVSvRRpVzJ7O+cZIy3jmDEiIqK31a0MeJf/39ihCd2A6hV4RowKhV7+q7p69So6deoEOzs7WFpawsfHBz/88INSnXPnzqFRo0awtLSEk5MTxo4di6Qk3h1DRET/783Ey8uViRgVGr07M3b48GF07NgRtWvXxrRp02BtbY179+4hKipKUef69eto0aIFqlWrhkWLFiEqKgrff/897t69i4MHD0oYPRERERkavUrGEhMTMXDgQAQEBCAsLAxGefyK+eqrr1C6dGmcPHkSNjY2AAB3d3cMHz4chw8fRuvWrXUZNhERERkwvTrnumXLFjx9+hSzZ8+GkZERXr9+DflbS1YkJibiyJEjGDBggCIRA4CBAwfC2toaO3bs0HXYREREZMD0Khk7evQobGxsEB0djSpVqsDa2ho2NjYYOXIkUlNTAQA3btxAZmYm6tWrp7SvqakpfH19ce3aNc0DSEnPniwu55GRWZDuEBERUWHLyFT+7k5J13kIenWZ8u7du8jMzETnzp0xdOhQzJ07FydPnsSyZcsQHx+PrVu3IiYmBgDg7Oyssr+zszNOnz6teQATVys/7/QB0KWh5u1R0ZSWAYz8/xtCVo4FzEykjYd0Z94FIFP3H9QkEXnW//5//lnAyFi6WKjwvIwCXj6WNAS9SsaSkpKQnJyMESNGKO6e7NatG9LT0/Hjjz9i5syZSElJAQCYmakuW2Rubq4o18j3gYCF6f+el+AfLpFeyUwHunwtdRSkKxmpwLJe2f/fcRJgYi5tPFQ4MjOArIz/PU9PBtYM0WkIenWZ0sLCAgDQt29fpe39+vUDAJw/f15RJy0tTWX/1NRURblmAZhmr02Z8zDRq1yXiIhI/5QwAcws//cwtdR5CHqVjLm4uAAAypYtq7Td0dERAPDy5UvF5cmcy5VviomJUbRBREREpAt6lYzVrVsXABAdHa20/fHj7GvBDg4O8PHxQYkSJXD58mWlOunp6bh+/Tp8fX11EisRERERoGfJWK9e2df2161bp7R97dq1KFGiBJo1awZbW1u0bNkSmzZtwqtXrxR1Nm7ciKSkJPTs2VOnMRMREZFh06tBTbVr18aQIUPw008/ITMzE02bNsXJkyfxyy+/YMqUKYpLkLNnz0bDhg3RtGlTBAYGIioqCgsXLkTr1q3Rtm1biXtBREREhkSvkjEAWLVqFcqXL4/Q0FDs2rULFSpUwOLFizF+/HhFnTp16uDo0aOYNGkSJkyYgJIlSyqmwiAiIiLSJb1LxkxMTBAcHIzg4OB31mvUqBHOnj2ro6iIiKjYMTEHPtsjdRRkAPRqzBgRERFRccNkjIiIiEhCTMaIiIiIJKR3Y8YkNXMzYCQDmvsCLWpLHQ0VFrn8f/8fEQVUrwAY8XcNEVGxdn0/cP0AIOTvr6tlTMa06Zv+2csgkf66chfYfPx/zxfvBEpbA/2aA3UrSxcXEREVjG9A9iMtGQjpo9ND8+c8UX5duQuE7AHik5S3v0zK3n7lrjRxERFRscZkjCg/5HJgy/F319l6QvkSJhERUT4wGSPKj4jo7DNg7/LiVXY9IiIiNTAZI8qPty9NFrQeERHR/2MyRpQfpay1W4+IiOj/MRkjyg+vctl3Tb6LXcnsekRERGpgMkaUH0ZG2dNXvEvfjzjfGBERqY3fHET5VbcyENRJ9VKkXcns7ZxnjIiINMBJX4nUUbcy4F0eCFqe/XxCN87AT0REBcJvECJ1vZl4ebkyESMiogLhmTFtGrcGkMkA27KATVmpo6HCIs/63//PPwsYGUsXC+lWCVOpIyCiwiLh2pQyIYTQ+VH1TGJiImxtbYGgbYCZpdThUGHLSAWW9cr+/zE7ABNzaeMhIiLt+f+1KRMSEmBjY6OTQ/L6ChEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYjJGBEREZGEmIwRERERSYgz8BOpy8Qc+GyP1FEQEZGe4JkxIiIiIgnxzJg2bfkckBkBvu0B3wCpoyEiIqL84tqUxRvXpiQiItITXJuSiIiIyLAwGSMiIiKSEJMxIiIiIgkxGSMiIiKSkF4nY7Nnz4ZMJoOPj49K2blz59CoUSNYWlrCyckJY8eORVJSkgRREhERkSHT26ktoqKiMGfOHFhZWamUXb9+HS1atEC1atWwaNEiREVF4fvvv8fdu3dx8OBBCaIlIiIiQ6W3ydjEiRPRoEEDZGVlIS4uTqnsq6++QunSpXHy5EnFbavu7u4YPnw4Dh8+jNatW0sRMhERERkgvbxM+ccffyAsLAxLlixRKUtMTMSRI0cwYMAApflDBg4cCGtra+zYsUOHkRIREZGh07tkLCsrC2PGjMGwYcNQo0YNlfIbN24gMzMT9erVU9puamoKX19fXLt2TVehEhEREenfZcpVq1bh0aNHOHr0aK7lMTExAABnZ2eVMmdnZ5w+fVrzg6cnKz83NgFKmGjeHhERERWuzAwgK+N/z9/+LtcBvUrGnj9/jm+++QbTpk2Dg4NDrnVSUlIAAGZmZipl5ubminKNrBmi/LxBH6BhP83bIyIiosJ18Rfgz22ShqBXydjXX38NOzs7jBkzJs86FhYWAIC0tDSVstTUVEW5Rob/BJi+sTalMc+KERERFWn+PYG6Xf73PD1Z9eRKIdObZOzu3btYvXo1lixZgsePHyu2p6amIiMjAw8fPoSNjY3i8mTO5co3xcTEwMXFRfMgTC25UDgREVFxUkL6IUV6M4A/OjoacrkcY8eOhYeHh+Jx4cIFREREwMPDAzNnzoSPjw9KlCiBy5cvK+2fnp6O69evw9fXV5oOEBERkUHSmzNjPj4+2LVrl8r2r7/+Gq9evcLSpUvh6ekJW1tbtGzZEps2bcK0adNQsmRJAMDGjRuRlJSEnj176jp0IiIiMmAyIYSQOojC1KxZM8TFxeHmzZuKbVevXkXDhg3h7e2NwMBAREVFYeHChWjSpAkOHTqk9jESExNha2sLBG3jZUoiIqLiLC0ZCOmDhIQEpflIC5PeXKZUR506dXD06FFYWFhgwoQJWL16NYYOHYqwsDCpQyMiIiIDozeXKfNy8uTJXLc3atQIZ8+e1W0wRERERG8xyDNjREREREUFkzEiIiIiCTEZIyIiIpKQ3o8Z06ktnwMyI8C3PeAbIHU0RERElF/X9wPXDwBCrvND6/3UFrrAqS2IiIj0BKe2ICIiIjIsTMaIiIiIJMRkjIiIiEhCTMaIiIiIJMRkjIiIiEhCTMaIiIiIJMRkjIiIiEhCTMaIiIiIJMRkjIiIiEhCTMaIiIiIJMS1KbWJa1MSEREVT1ybsnjj2pRERER6gmtTEhERERkWJmNEREREEmIyRkRERCQhJmNEREREEmIyRkRERCQhJmNEREREEmIyRkRERCQhJmNEREREEmIyRkRERCQhJmNEREREEuLalNrEtSmJiIiKJ65NWbxxbUoiIiI9wbUpiYiIiAwLkzEiIiIiCTEZIyIiIpIQkzEiIiIiCTEZIyIiIpKQXiVjly5dwujRo1G9enVYWVmhfPny6NWrFyIiIlTqhoeHo23btrC2toadnR0+/vhjPHv2TIKoiYiIyJDp1Txj8+fPx9mzZ9GzZ0/UrFkTT548wfLly1GnTh38+eef8PHxAQBERUWhSZMmsLW1xZw5c5CUlITvv/8eN27cwMWLF2FqaipxT4iIiMhQ6FUy9tlnn2HLli1KyVTv3r1Ro0YNzJs3D5s2bQIAzJkzB69fv8aVK1dQvnx5AIC/vz9atWqF9evXIzAwUJL4iYiIyPDo1WXKhg0bqpzVqly5MqpXr47w8HDFtl9//RUdOnRQJGIA0LJlS3h5eWHHjh06i5eIiIhIr5Kx3Agh8PTpU9jb2wMAoqOjERsbi3r16qnU9ff3x7Vr1zQ/WHpy9sy9OY/MDM3bIiIiosKXmaH83Z2erPMQ9OoyZW42b96M6OhozJw5EwAQExMDAHB2dlap6+zsjBcvXiAtLQ1mZmbqH2zNEOXnDfoADfup3w4RERHpxsVfgD+3SRqCXidjd+7cQVBQED744AN88sknAICUlBQAyDXZMjc3V9TRKBkb/hNg+sbalMYm6rdBREREuuPfE6jb5X/P05NVT64UMr1Nxp48eYKAgADY2toiLCwMxsbGAAALCwsAQFpamso+qampSnXUZmrJhcKJiIiKkxIm2Q8pQ5D06IUkISEB7dq1Q3x8PE6fPg0XFxdFWc7lyZzLlW+KiYmBnZ2dZmfFiIiIiDSgd8lYamoqOnbsiIiICBw9ehTe3t5K5eXKlYODgwMuX76ssu/Fixfh6+uro0iJiIiI9OxuyqysLPTu3Rvnz5/HL7/8gg8++CDXet27d8e+ffsQGRmp2Hbs2DFERESgZ8+eugqXiIiICDIhhJA6CG0ZP348li5dio4dO6JXr14q5QMGDAAAREZGonbt2ihVqhTGjRuHpKQkLFiwAK6urrh06ZLalykTExNha2sLBG3jmDEiIqLiLC0ZCOmDhIQE2NjY6OSQepWMNWvWDKdOncqz/M2u3rp1C5999hnOnDkDU1NTBAQEYOHChShbtqzax2UyRkREpCckSMb0aszYyZMn8123evXqOHToUOEFQ0RERJQPejVmjIiIiKi4YTJGREREJCEmY0REREQS0qsxY5Lb8jkgMwJ82wO+AVJHQ0RERPl1fT9w/QAg5Do/tF7dTSkV3k1JRESkJyS4m5KXKYmIiIgkxGSMiIiISEJMxoiIiIgkxGSMiIiISEJMxoiIiIgkxGSMiIiISEJMxoiIiIgkxGSMiIiISEJMxoiIiIgkxOWQtInLIRERERVPXA6peONySERERHqCyyERERERGRYmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQSYjJGREREJCEmY0REREQS4tqU2sS1KYmIiIonrk1ZvHFtSiIiIj3BtSmJiIiIDAuTMSIiIiIJMRkjIiIikhCTMSIiIiIJMRkjIiIikpDBJmNpaWmYNGkSXFxcYGFhgfr16+PIkSMatwUAyMzQYoTFQGYGcG6L4fUbYN8Nse+G2m+AfWffpY5Et/6/v4rvdh0w2GRs0KBBWLRoEfr374+lS5fC2NgY7du3x5kzZ9RuS/GGyQ3sH2xWBvDntuz/Ghr23fD6bqj9Bth39l3qSHRLzmRMJy5evIht27Zh7ty5WLBgAQIDA3H8+HFUqFABX375pW6CuL6/8PfRxTE0UVT7Yah9L4r91mQffXnPNdlHX/peVF8rdRXVfrDvRZZBJmNhYWEwNjZGYGCgYpu5uTmGDh2K8+fPIzIysvCDuH6g8PfRxTE0UVT7Yah9L4r91mQffXnPNdlHX/peVF8rdRXVfrDvRZZBLod07do1eHl5qcys6+/vDwC4fv063Nzc8t2eYhGDpBfKBUYmQAmTPHaSZ8/yqw519ynsY6QnK/+3sOLSxWul7j761Hd162vS96LYD3X30dV7rsk++tL3ovhaGeq/d8Bw+p6ZoTzM6P+/y3W5QJFBLofk4+ODsmXL4tixY0rbb9++jerVq2PVqlX49NNP893e/fv34enpqe0wiYiISCL37t1DxYoVdXIsgzwzlpKSAjMzM5Xt5ubminJ1uLu74969ezAxMYFMJlNsNzMzy/U4REREVDSkpaUpDdYXQiAjIwPu7u46i8EgkzELC4tc75JITU1VlKvDyMhIZ9kzERER6ReDHMDv7OyMmJgYle0521xcXHQdEhERERkog0zGfH19ERERgcTERKXtFy5cUJQTERER6YJBJmM9evRAVlYWVq9erdiWlpaG0NBQ1K9fX607KYmIiIgKwiDHjNWvXx89e/bElClTEBsbi0qVKmHDhg14+PAh1q1bJ3V4REREZEAMcmoLIHuw/rRp07Bp0ya8fPkSNWvWxKxZs9CmTRupQyMiIiIDYrDJGBEREVFRYJBjxoiIiIiKCiZjRERERBJiMkZEREQkISZjRERERBJiMkZEREQkISZjRERERBJiMkZEREQkISZjRERERBJiMkZEREQkIYNcm1Lb5HI5Hj58CBMTE8hkMsV2MzMzmJmZSRgZERERvUtaWhrS0tIUz4UQyMjIgLu7O4yMdHPOismYFjx8+BCenp5Sh0FERERacu/ePVSsWFEnx2IypgUmJiYAgIsXL8LZ2Vmx/V1nxvz8/HDp0iW1jqPuPoV9jMTERLi5uSEyMhI2NjaFFpcuXit199GnvqtbX5O+F8V+qLuPrt5zTfbRl74XxdfKUP+9A4bT97fPjMXExMDf31/x3a4LTMa0IOfSpLOzM1xdXfO1j7GxsVofaprso4tjAICNjU2hxqWrfhhq3zU5BqBe34tqP4rie67JPvrS96L6WgGG++8dMNy+vznsqLBxAL9EgoKCCn0fXRxDE0W1H4ba96LYb0320Zf3XJN99KXvRfW1UldR7Qf7XnTJhBBC6iCKu6ioKMWp3PyeGdMHiYmJsLW1RUJCgka/tooz9t3w+m6o/QbYd/bdsPouxXc6z4xpQc64MEO7c9LMzAzBwcEG12+AfTfEvhtqvwH2nX03rL5L8Z3OM2NaYKi/HoiIiPSNFN/pPDOmRX5+fvD29kZISIjUoRAREZEaQkJC4O3tDT8/P50fm2fGtIBnxoiIiPQDz4wRERERGRgmY0REREQSYjJGREREJCEmY0REREQSYjKmRbybkoiIqHji3ZTFHO+mJCIi0g+8m5KIiIjIwJTIT6UhQ4ao3bBMJsO6devU3o+IiIjIkOQrGTt+/DhkMplaDatbn4iIiMgQ5SsZe/jwYSGHQURERGSYOGaMiIiISEIFSsb+/PNPzJ07FxMmTMDdu3cBAMnJybh69SqSkpK0EmBxwqktiIiIiqdiN7VFeno6+vTpg99++w1CCMhkMhw5cgTNmzdHamoqXF1dMWHCBEydOrUwYi5yOLUFERGRfig2U1tMmzYN+/btw8qVK/HPP//gzXzO3NwcPXv2xG+//aa1IImIiIj0lUbJ2NatWzFy5EgEBgbCzs5OpbxatWq4f/9+gYMjIiIi0ncaJWOxsbGoUaNGnuXGxsZITk7WOCgiIiIiQ6FRMubm5oY7d+7kWX727FlUqlRJ46CIiIiIDIVGyVi/fv3w448/4vz584ptOZO8rlmzBjt27MDAgQO1E2ExwrspiYiIiqdieTdlx44dcfz4cVSrVg23bt1CjRo18OLFC0RFRaF9+/b47bffYGxsXBgxFzm8m5KIiEg/FJu7KU1NTfH7778jNDQUFStWRNWqVZGWloaaNWti/fr12Lt3r0aJWFJSEoKDg9G2bVvY2dlBJpNh/fr1+dr32LFjGDJkCLy8vGBpaYmKFSti2LBhiImJUanbrFkzyGQylUfbtm3VjpmIiIioIPK1HFJuZDIZBgwYgAEDBmgtmLi4OMycORPly5dHrVq1cPLkyXzvO2nSJLx48QI9e/ZE5cqVcf/+fSxfvhz79u3D9evX4eTkpFTf1dUVc+fOVdrm4uKijW4QERER5ZvGyVhhcHZ2RkxMDJycnHD58mW1rtsuWrQIjRo1gpHR/072tW3bFk2bNsXy5cvx7bffKtW3tbXVaiJJREREpIl8JWPNmzdXu2GZTIZjx46ptY+ZmZnKGaz8atKkSa7b7OzsEB4enus+mZmZSE1NhbW1tUbHJCIiIiqofCVjcrlccbdkjsjISNy/fx+2traoWLEiAODBgweIj4+Hp6cn3NzctB+tmpKSkpCUlAR7e3uVsoiICFhZWSE9PR1ly5bF8OHD8c0338DExETj4yUmJio9NzMzg5mZmcbtERERUeFKS0tDWlqa4vnb3+W6kK9k7O2xW2fOnEGnTp2wZs0afPLJJyhRIruZzMxMhIaGYtKkSfkeeF+YlixZgvT0dPTu3Vtpu6enJz766CPUqFEDr1+/RlhYGL799ltERERg+/btGh/v7QQ0ODgY06dP17g9IiIiKlxz587FjBkzJI1Bo6ktGjRogA8//BALFy7Mtfzzzz/HmTNncOHCBY0DyxkzFhoaikGDBqm9/x9//IEWLVqgW7du+UqwAgMDsWbNGpw/fx4NGjRQ61g5t8FGRkYq3QbLM2NERERFW25nxtzc3Ir+1BZ///234tJkbjw8PHDjxg2NgyqoO3fuoGvXrvDx8cHatWvztc/nn38OADh69KjGx7WxsVF6MBEjIiIq2szMzFS+v3VNo2TMxcUF27dvR2ZmpkpZZmYmtm/fLtk0EZGRkWjdujVsbW1x4MABlCxZMl/75VxifPHiRWGGR0RERKREo6ktvvzyS4wYMQINGjTAiBEjFOtQ3r17F6tWrcL169exYsUKrQaaH8+fP0fr1q2RlpaGY8eOwdnZOd/73r9/HwDg4OBQWOERERERqdAoGQsMDISxsTGmTp2KwMBAxZ2WQgg4ODhg1apVGD58uFYDfVNMTAwSEhLg6empuPvx9evXaN++PaKjo3HixAlUrlw5130TExNVxnIJIRTzkLVp06bQ4iYiIiJ6m8aTvg4dOhSffPIJLl++jEePHgEAKlSogHr16inurtTE8uXLER8fj8ePHwMA9u7di6ioKADAmDFjYGtriylTpmDDhg148OAB3N3dAQD9+/fHxYsXMWTIEISHhyvNLWZtbY0uXboAAK5evYq+ffuib9++qFSpElJSUrBr1y6cPXsWgYGBqFOnjsax+/n5wdjYGEFBQQgKCtK4HSIiItKtkJAQhISEICsrS+fH1uhuysLk7u6uSO7elpN8DRo0SCUZe9d+FSpUwMOHDxVtTJo0CZcuXcKTJ09gZGSEatWqYfjw4Upn+dTBhcKJiIj0gxTf6RonY1lZWdi0aRP279+vdGasQ4cO6N+/v0YLhRdXTMaIiIj0Q7FJxhISEtCmTRtcunQJJUuWVJqBPzExEf7+/jh06JDBJCZMxoiIiPSDFN/pGk1tMXXqVFy5cgXLli3Ds2fPcPXqVVy9ehWxsbFYvnw5Ll++jKlTp2o7ViIiIiK9o9GZsXLlyqFHjx5YunRpruVjx45FWFiYYhC+vuOZMSIiIv1QbM6MPX/+HFWqVMmzvGrVqpw8lYiIiCgfNErGKlWqhD179uRZvmfPHnh6emocVHHl5+cHb29vhISESB0KERERqSEkJATe3t7w8/PT+bE1uky5YsUKjB49Gm3btsX48ePh5eUFAPjnn3/www8/4Pfff8fy5csxcuRIrQdcFPEyJRERkX6Q4jtdo9lZR40ahdjYWMybNw+HDh1SKjMxMcE333xjMIkYERERUUEUaNLXuLg4HD16VGmesZYtW8Le3l5rARYHPDNGRESkH4rNmbEc9vb26NOnj7ZiISIiIjI4BUrGMjIyEB0djZcvXyK3E2wFWeeRiIiIyBBolIzFx8dj4sSJ2Lx5M9LT01XKhRCQyWSSLLYpJS4UTkREVDwVu4XCu3Tpgr1796JPnz6oX78+bG1tc633ySefFDjA4oBjxoiIiPRDsRkzdvjwYYwdOxaLFy/WdjxEREREBkWjSV/LlCmDSpUqaTsWIiIiIoOjUTIWGBiIbdu2QS6XazseIiIiIoOi0WXKadOmIS0tDfXq1cPHH38MV1dXGBsbq9Tr1q1bgQMkIiIi0mcaJWPR0dE4fvw4rl+/juvXr+daxxDvpiQiIiJSl0bJ2JAhQ3D16lVMmTLlnXdTGhpObUFERFQ8FbupLaysrDBx4kTMmDGjMGLKl6SkJCxYsAAXLlzAxYsX8fLlS4SGhmLQoEH52j8+Ph5ffvkldu3aheTkZPj7+2PhwoUaTVTLqS2IiIj0gxTf6RoN4HdycoKdnZ22Y1FLXFwcZs6cifDwcNSqVUutfeVyOQICArBlyxaMHj0a3333HWJjY9GsWTPcvXu3kCImIiIiUqVRMvb5559j7dq1SEpK0nY8+ebs7IyYmBg8evQICxYsUGvfsLAwnDt3DuvXr0dwcDCCgoJw8uRJGBsbIzg4uJAiJiIiIlKl0Zix1NRUmJiYoFKlSujVqxfc3NxU7qaUyWSYMGGCVoLMjZmZGZycnDTaNywsDGXLllW629PBwQG9evXCpk2bkJaWBjMzM22FSkRERJQnjZKxiRMnKv5/+fLludYp7GSsIK5du4Y6derAyEj5xKC/vz9Wr16NiIgI1KhRQ+12z549CysrK8VzR0dHVK1aFampqbh9+7ZK/Zzxaf/88w9ev36tVObu7g47Ozs8e/YMkZGRSmUlS5ZE5cqVkZWVhb/++kul3Ro1asDExAT37t1DQkKCUlm5cuVQtmxZvHz5Eg8ePFAqs7CwQLVq1QBkv0ZvDyesVq0aLCws8OjRIzx//lyprGzZsihXrhxevXqlcqnXxMRE8XreuHEDGRkZSuWVK1dGyZIlER0djadPnyqVlSlTBhUqVEBKSgrCw8OVymQyGWrXrg0ACA8PR0pKilK5h4cHSpcujadPnyI6OlqpzNbWFp6ensjIyMCNGzfwtlq1asHY2Bh3797Fq1evlMrc3Nzg4OCAFy9e4OHDh0plVlZWqFKlCgDg6tWrKu16e3vD3NwcDx48wMuXL5XKnJ2d4ezsjMTERPz7779KZWZmZqhevToA4O+//0ZmZqZSuZeXF6ytrREVFYXY2FilMnt7e5QvXx7Jycm4c+eOUpmRkRF8fX0BALdv30ZqaqpSecWKFVGqVCk8efIEjx8/ViorVaoUKlasiPT0dNy8eVOlr76+vjAyMkJERITKWfTy5cvD3t4ecXFx+O+//5TKrK2t4eXlBblcnuvd2j4+PjA1NcX9+/cRHx+vVObi4gInJyfEx8fj/v37SmXm5ubw9vYGAFy/fl1lnsSqVavC0tIS//33H+Li4pTKHB0d4erqiqSkJERERCiVlShRAjVr1gQA3Lp1C2lpaUrllSpVgo2NDWJiYhATE6NUVrp0aXh4ePAzgp8RCvyMyKbLz4j09HSkp6fD3NwcVatWRWJiokqbhU5o4OHDh/l66MqlS5cEABEaGpqv+lZWVmLIkCEq2/fv3y8AiN9//12t4yckJAgAKo8aNWoIIYS4e/duruU5GjRooFK2ceNGIYQQy5cvVylr3br1O48bGxsrhBCiY8eOKmULFy4UQgixY8cOlbLatWsrYjI1NVUpv3nzphBCiKFDh6qUTZ48WQghxIkTJ1TKypUrp2i3XLlyKuUnTpwQQggxefJklbKhQ4cKIYS4efOmSpmpqami3dq1a6uU79ixQwghxMKFC1XKOnbsKIQQIjY2NtfXMCEhQQghROvWrVXKli9fLoQQYuPGjSplDRo0UMSUW7t3794VQgjRv39/lbLg4GAhhBC///67Spmnp6eiXXt7e5Xyc+fOCSGEmDBhgkrZqFGjhBBCXLlyRaWsZMmSina9vb1Vyn/77TchhBBz5sxRKevRo4cQQojIyMhc+5qamiqEEKJp06YqZWvWrBFCCLFmzRqVsqZNmwohhEhNTc213cjISCGEED169FApmzNnjhBCiN9++02lzNvbW9HXkiVLqpRfuXJFCCHEqFGjVMomTJgghBDi3LlzKmX29vaKdj09PVXKcz5PgoODVcr69+8vhOBnBD8j+BnxdpkUnxF5vce6oNHdlEXN5cuX4efnl++7KY2NjfHpp59ixYoVStuPHz+OFi1aYNeuXejSpUu+j59z58WBAwd4Zoy/evmr9/8VlV+9OXhm7H/4GZGNnxHZDP0zIrczY25ubjq9m9IgkzFra2v07t0b69atU9p+4MABBAQE4Pfff0ebNm3yfXxObUFERKQfis3UFsVdzp2Yb8vZ5uLiouuQiIiIyEAZZDLm6+uLq1evqlymuHDhAiwtLeHl5SVRZERERGRo9D4Zi4mJwZ07d5TGIPTo0QNPnz7Fzp07Fdvi4uLwyy+/oGPHjpzWgoiIiHQm31NbZGVlqcwlJrXly5cjPj5eMWhw7969iIqKAgCMGTMGtra2mDJlCjZs2IAHDx7A3d0dQHYy1qBBAwwePBi3b9+Gvb09VqxYgaysLEmXeCIiIiLDk+8B/KVKlUKbNm3QoUMHtG3bFg4ODoUd23u5u7vj0aNHuZblJF+DBg1SScYA4OXLl/jiiy+we/dupKSkwM/PD99//z3q1aundhw5g/28vLy4UDgREVEx9OZC4REREUXzbsply5bhwIEDOHXqFNLT01G3bl0EBAQgICAAdevWLew4izTeTUlERKQfpPhOV3tqi5SUFBw9ehQHDx7EgQMH8N9//8HJyQlt27ZFhw4d0KpVK5QsWbKw4i2SmIwRERHph2KRjL3t5s2b2L9/Pw4ePIhz585BJpOhUaNGaN++PQICAlC1alVtxVpkMRkjIiLSD8UyGXtTQkICDh06hAMHDuD333/Hs2fPMGfOHEyaNElbhyiSmIwRERHph2KfjL3t0qVLAAA/P7/COkSRwGSMiIhIP+jdDPx+fn56n4i9yc/PD97e3ggJCZE6FCIiIlJDSEgIvL29Jclb9GJtSqnxzBgREZF+0LszY0RERET0bkzGiIiIiCTEZIyIiIhIQkzGiIiIiCSU74XC3/bkyROsW7cOV69eRUJCAuRyuVK5TCbDsWPHChwgERERkT7TKBn7+++/0axZM6SkpKBKlSq4ceMGvL29ER8fj+joaHh6esLNzU3bsRZ5fn5+XCiciIioGHpzoXBd02hqi/bt2+PmzZs4c+YMLC0t4ejoiKNHj6J58+b45ZdfMHLkSBw4cAD+/v6FEXORw6ktiIiI9EOxmdri7Nmz+PTTT1G+fHkYGWU3kXOZsmfPnujfvz+++OIL7UVJREREpKc0SsbkcjnKli0LAChVqhSMjY3x4sULRXmNGjVw5coV7URIREREpMc0SsY8PDzw4MGD7AaMjODh4YGjR48qys+dO4dSpUppJUAiIiIifaZRMta6dWv88ssviucjR47E2rVr0bJlS7Ro0QIbNmxAv379tBYkERERkb7SaAD/y5cvcf/+fdSsWRMmJiYQQmD27Nn49ddfYWxsjA4dOuCrr76CqalpYcRc5OQM9vPy8uLdlAbg9evXsLa2BgAkJSXByspK4oiIiKig3rybMiIiQqcD+LlQuBbwbkrDwmSMiEh/FZu7KZs3b/7OCV1PnDiB5s2baxRQWloaJk2aBBcXF1hYWKB+/fo4cuTIe/dzd3eHTCbL9VG5cmWlunnVmzdvnkYxExEREWlKo0lfT548iWHDhuVZHhsbi1OnTmkU0KBBgxAWFobx48ejcuXKWL9+Pdq3b48TJ06gUaNGee63ZMkSJCUlKW179OgRvv76a7Ru3VqlfqtWrTBw4EClbbVr19YoZiIiIiJNabwckkwmy7Ps33//RcmSJdVu8+LFi9i2bRsWLFiAiRMnAgAGDhwIHx8ffPnllzh37lye+3bp0kVl27fffgsA6N+/v0qZl5cXBgwYoHaMRERERNqU72Rsw4YN2LBhg+L5t99+izVr1qjUi4+Px99//4327durHUxYWBiMjY0RGBio2GZubo6hQ4fiq6++QmRkpFrLLG3ZsgUeHh5o2LBhruUpKSmQyWQwNzdXO1YiIiIibcj3mLHk5GQ8e/YMz549AwC8evVK8TznERcXBzMzM4wYMQJr165VO5hr167By8tLZcBczrJK169fV6ut8PDwPKfYWL9+PaysrGBhYQFvb29s2bJF7XjflpiYqPRIS0srcJtERERUeNLS0lS+v3Ut32fGRo4ciZEjRwLInvR16dKl6NSpk1aDiYmJgbOzs8r2nG2PHz/Od1ubN28GkPslyoYNG6JXr17w8PDA48ePERISgv79+yMhIUHRR028fdYuODgY06dP17g9KpreXET2jz/+QOvWrWFsbCxhREREpKm5c+dixowZksZQpKa28PT0RJUqVXDgwAGl7ffv34enpycWL16M8ePHv7cduVyO8uXLw9HREVevXn1v/fT0dNStWxdRUVF4/PgxLCws1Io75zbYyMhIpbN6ZmZmMDMzU6stKtp27tyJsWPHIjo6WrHN1dUVS5cuRbdu3SSMjIiINJGWlqZ0JSsxMRFubm5Ff2qLN7169QpRUVH477//VB7qsrCwyPXSXmpqqqI8P06dOoXo6Ohcz4rlxtTUFKNHj0Z8fHyB1tS0sbFRejAR0y87d+5Ejx49lBIxAIiOjkaPHj2wc+dOiSIjIiJNmZmZqXx/65rGd1OuXLkSixYtwv379/Os8+blnPxwdnZW+aIDsi9fAoCLi0u+2tm8eTOMjIzQt2/ffB875xLjmwueE+XIysrCuHHjkNuJZCEEZDIZxo8fj86dO/OSJRERqUWjM2OrVq1CUFAQKlWqhG+//RZCCIwfPx6TJ0+Gk5MTatWqhXXr1qndrq+vLyIiIlQGz124cEFR/j5paWn49ddf0axZs3wnbwAUSaWDg0P+AyaDcfr0aURFReVZLoRAZGQkTp8+rcOoiIhIH2iUjC1btgxt2rTBwYMHFdNQBAQEYPbs2bh9+zZevXqF58+fq91ujx49kJWVhdWrVyu2paWlITQ0FPXr11ecvfrvv/9w586dXNs4cOAA4uPj87xEmXM36JtevXqFJUuWwN7eHnXr1lU7btJ/OWdntVWPiIgoh0aXKe/du6dYCNvExARA9iB4ALC1tcWwYcOwYsUKfP7552q1W79+ffTs2RNTpkxBbGwsKlWqhA0bNuDhw4dKZ9oGDhyIU6dO5XrJaPPmzTAzM0P37t1zPUZISAh2796Njh07onz58oiJicFPP/2E//77Dxs3bizQ4uZ+fn5cKFxP5XaXb0HqERFR0fLmQuG6plEyZmtri8zMTADZg9YtLS0RGRmpKC9ZsiSePHmiUUA///wzpk2bho0bN+Lly5eoWbMm9u3bhyZNmrx338TEROzfvx8BAQGwtbXNtc6HH36Ic+fOYe3atXj+/DmsrKzg7++Pn376SeP1NHNcunSJC4XrqcaNG8PV1RXR0dG5/giQyWRwdXVF48aNJYiOiIgKKudESs4MCbqk0dQWLVq0gLu7u+JsVcuWLfHixQvs3bsXcrkcHTp0gJGREa5du6b1gIsiKVZ4J93LuZsSgFJClrM0WFhYGKe3ICIq5qT4TtdozNiAAQNw8+ZNxTQUM2bMQHh4OMqXLw93d3f8888/inUhifRFt27dEBYWpnJjiKurKxMxIiLSmNYmfb1//z727t0LY2NjtG7dGl5eXtpotljgmTHD8uYp7AMHDnAGfiIiPSLFd3qRmoG/uGIyZlhev34Na2trAEBSUhKsrKwkjoiICgP/1g1TsblMSURERETaka9kzMjICMbGxmo/DI2fnx+8vb0REhIidShERESkhpCQEHh7e8PPz0/nx87X1BbffPON4o6xHLt27cKtW7fQpk0bVKlSBQBw584dHD58GD4+PujSpYvWgy3qOLUFERFR8STl1Bb5SsamT5+u9Hz16tWIjY3FzZs3FYlYjvDwcDRv3lytpYiIiIiIDJVGY8YWLFiA0aNHqyRiAFCtWjWMHj0a3333XYGDIyIiksqbM7H/8ccfkszMToZBo2QsKipKsQxSbkxMTN65qDIREVFRtnPnTnh7eyuet2/fHu7u7ti5c6eEUZG+0igZ8/HxwYoVKxAdHa1SFhUVhRUrVqBGjRoFDo6oKLKysoIQAkII3upOpIdyVtt4+zsuOjoaPXr0YEJGWqfRPGNnzpxBmzZtAABdu3ZFpUqVAAB3797F7t27IYTA4cOH0ahRI+1GW0TlDPbz8vLiQuFERMVYVlYW3N3d87y6k7MO7YMHDwxy1gB99uZC4REREcVj0tebN29i2rRpOHz4MFJSUgAAFhYWaNOmDWbMmGFQZ8Y46SsRkX44efIkPvroo/fWO3HiBJo1a1b4AZHOSfGdnq+7KXPj4+ODXbt2QS6X49mzZwAABwcHGBlxHlkiIiqeYmJitFqPKD80TsZyGBkZoWzZstqIhYiISFLOzs5arUeUHzyNRURE9P8aN24MV1dXlYnOc8hkMri5uaFx48Y6joz0GZMxIiKi/2dsbIylS5cCgEpClvN8yZIlHLxPWsVkjIiI6A3dunVDWFiYykoyrq6uCAsLQ7du3SSKjPQVkzEt4kLhRET6oVu3brh9+7bi+YEDB/DgwQMmYnpMyoXCi20ylpaWhkmTJsHFxQUWFhaoX78+jhw58t79pk+fDplMpvIwNzcvUCxA9vxrt2/fNpg5xtLS0jB9+nRF/w0J+254fTfUfgOG2/c3L0XWr1/f4C5NGtr7HhQUhNu3b+PMmTMAoNN+azzPmNT69u2LsLAwjB8/HpUrV8b69etx6dIlnDhx4p2TzU6fPh0zZszAypUrYW1trdhubGyMvn37ahRLVFQU3NzcEBkZCVdXV43aKI4MeX419t3w+m6o/QYMt++vX79WfE88fvzY4O6gNNT3XYrv9GJ5ZuzixYvYtm0b5s6diwULFiAwMBDHjx9HhQoV8OWXX+arjR49emDAgAGKh6aJmKY0uZSp7j66OIYmimo/DLXvRbHfmuyjL++5JvvoS9+L6mulrqLaD/a9CBPF0BdffCGMjY1FQkKC0vY5c+YIAOK///7Lc9/g4GABQMTGxoqEhAQhl8sLHE9kZKQAICIjI/O9T7Vq1dQ+jrr7FPYxEhISBACV90Gbx9Ckvi720ae+q1tfk74XxX6ou4+u3nNN9tGXvhfF18pQ/70LYbh91+Q7vaAKPOmrFK5duwYvLy+V06b+/v4AgOvXr8PNze2dbVSsWBFJSUmwsrJCly5dsHDhQo0nrxX/f6X37RmZzczMYGZmlus+WVlZSExMVOs46u5T2MfIqVfc+6HJPvrUd3Xra9L3otgPdffR1XuuyT760vei+FoZ6r93wHD6npaWpjQ+LOe7XOhwFFexHDPm4+ODsmXL4tixY0rbb9++jerVq2PVqlX49NNPc9136dKl+Pfff/HBBx/AzMwMp0+fRkhICDw8PHD58mWNrovfv38fnp6eGvWFiIiIip579+6hYsWKOjlWsTwzlpKSkusZp5w7InMWLs/NuHHjlJ53794d/v7+6N+/P1asWIHJkyerHY+7uzvu3bsHExMTpUkC33VmjIiIiKT39pkxIQQyMjLg7u6usxiKZTJmYWGR6y2nqampinJ19OvXD59//jmOHj2qUTJmZGSks+yZiIiI9EuxvJvS2dlZZXwW8L/rvG/Pmpwfbm5uePHiRYFjIyIiIlJHsUzGfH19ERERoTIY78KFC4pydQgh8PDhQzg4OGgrRCIiIqJ8KZbJWI8ePZCVlYXVq1crtqWlpSE0NBT169dX3En533//4c6dO0r7Pnv2TKW9lStX4tmzZ2jbtm3hBk5ERET0lmJ5NyUA9OrVC7t27cKECRNQqVIlbNiwARcvXsSxY8fQpEkTAECzZs1w6tQppdtTLS0t0bt3b9SoUQPm5uY4c+YMtm3bhlq1auHs2bOwtLSUqktERERkgIrlAH4A+PnnnzFt2jRs3LgRL1++RM2aNbFv3z5FIpaX/v3749y5c/j111+RmpqqmLV/6tSpTMSIiIhI54rtmTEiIiIifVAsx4wRERER6QsmY0REREQSYjJGREREJCEmY0REREQSYjKG7DnKJk2aBBcXF1hYWKB+/fo4cuRIvvaNjo5Gr169UKpUKdjY2KBz5864f/9+rnXXrVuHatWqwdzcHJUrV8ayZcu02Q2NaNr3nTt3onfv3qhYsSIsLS1RpUoVfP7554iPj1ep6+7uDplMpvIYMWJEIfQo/zTt+/Tp03PtT87aqG/Tp/c9r/dSJpOhcuXKSnXzqjdv3rzC6tZ7JSUlITg4GG3btoWdnR1kMhnWr1+f7/3j4+MRGBgIBwcHWFlZ4aOPPsLVq1dzrbtnzx7UqVMH5ubmKF++PIKDg5GZmamlnqivIH0/duwYhgwZAi8vL1haWqJixYoYNmxYriuhNGvWLNf3Xap5HAvS7/Xr1+f57/jJkycq9fXpPc/rfZTJZDAxMVGqW9Q+4y9duoTRo0ejevXqsLKyQvny5dGrVy9ERETka38p/s6L7dQW2jRo0CCEhYVh/PjxqFy5MtavX4/27dvjxIkTaNSoUZ77JSUl4aOPPkJCQgK++uormJiYYPHixWjatCmuX7+OMmXKKOr++OOPGDFiBLp3747PPvsMp0+fxtixY5GcnIxJkybpopu50rTvgYGBcHFxwYABA1C+fHncuHEDy5cvx4EDB3D16lWV9UF9fX3x+eefK23z8vIqlD7ll6Z9z7Fy5UpYW1srnhsbG6vU0bf3fcmSJUhKSlLa9ujRI3z99ddo3bq1Sv1WrVph4MCBSttq166tnU5oIC4uDjNnzkT58uVRq1YtnDx5Mt/7yuVyBAQE4K+//sIXX3wBe3t7rFixAs2aNcOVK1eUktGDBw+iS5cuaNasGZYtW4YbN27g22+/RWxsLFauXFkIPXu/gvR90qRJePHiBXr27InKlSvj/v37WL58Ofbt24fr16/DyclJqb6rqyvmzp2rtE2TZeq0oSD9zjFz5kx4eHgobStVqpTSc317z6dOnYphw4YpbXv9+jVGjBiR6996UfqMnz9/Ps6ePYuePXuiZs2aePLkCZYvX446dergzz//hI+PT577SvZ3LgzchQsXBACxYMECxbaUlBTh6ekpPvjgg3fuO3/+fAFAXLx4UbEtPDxcGBsbiylTpii2JScnizJlyoiAgACl/fv37y+srKzEixcvtNQb9RSk7ydOnFDZtmHDBgFArFmzRml7hQoVVPoutYL0PTg4WAAQz549e2c9fXzfczNr1iwBQJw9e1ZpOwARFBRU4Hi1KTU1VcTExAghhLh06ZIAIEJDQ/O17/bt2wUA8csvvyi2xcbGilKlSom+ffsq1fX29ha1atUSGRkZim1Tp04VMplMhIeHF7wjGihI30+dOiWysrJUtgEQU6dOVdretGlTUb16da3ErA0F6XdoaKgAIC5duvTeuvr2nudm48aNAoDYvHmz0vai9hl/9uxZkZaWprQtIiJCmJmZif79+79zX6n+zg3+MmVYWBiMjY0RGBio2GZubo6hQ4fi/PnziIyMfOe+fn5+8PPzU2yrWrUqWrRogR07dii2nThxAs+fP8eoUaOU9g8KCsLr16+xf/9+LfYo/wrS92bNmqls69q1KwAgPDw8133S09Px+vXrggWtJQXpew4hBBITE5VWeHiTPr7vudmyZQs8PDzQsGHDXMtTUlKQmppaoJi1xczMTOUsTn6FhYWhbNmy6Natm2Kbg4MDevXqhd9++w1paWkAgNu3b+P27dsIDAxEiRL/u/gwatQoCCEQFhZWsE5oqCB9b9KkCYyMjFS22dnZ5fn3npmZqXIWVQoF6febXr16haysrFzL9PE9z82WLVtgZWWFzp0751peVD7jGzZsCFNTU6VtlStXRvXq1fP895pDqr9zg0/Grl27Bi8vL9jY2Cht9/f3BwBcv3491/3kcjn+/vtv1KtXT6XM398f9+7dw6tXrxTHAKBSt27dujAyMlKU65qmfc9LzhgKe3t7lbLjx4/D0tIS1tbWcHd3x9KlSzULWku00feKFSvC1tYWJUuWxIABA/D06VOVYwD6/b5fu3YN4eHh6NevX67l69evh5WVFSwsLODt7Y0tW7ZoHLfUrl27hjp16qgkJf7+/khOTlaMR8nrfXdxcYGrq6tk77u2JSUlISkpKde/94iICFhZWaFkyZJwcnLCtGnTkJGRIUGU2vHRRx/BxsYGlpaW6NSpE+7evatUbgjv+bNnz3DkyBF06dIFVlZWKuVF7TP+bUIIPH36NNd/r2+S6u/c4MeMxcTEwNnZWWV7zrbHjx/nut+LFy+Qlpb23n2rVKmCmJgYGBsbw9HRUameqakpypQpk+cxCpumfc/L/PnzYWxsjB49eihtr1mzJho1aoQqVarg+fPnWL9+PcaPH4/Hjx9j/vz5mnegAArS99KlS2P06NH44IMPYGZmhtOnTyMkJAQXL17E5cuXFUmOIbzvmzdvBpC9zNjbGjZsiF69esHDwwOPHz9GSEgI+vfvj4SEBIwcOVLD6KUTExOT63Jrb75uNWrUUAxqz+s1lup917YlS5YgPT0dvXv3Vtru6emJjz76CDVq1MDr168RFhaGb7/9FhEREdi+fbtE0WrG0tISgwYNUiRjV65cwaJFi9CwYUNcvXoVbm5uAGAQ7/n27duRmZmZ6996UfyMf9vmzZsRHR2NmTNnvrOeVH/nBp+MpaSkwMzMTGV7zp1xKSkpee4HIF/7pqSkqJwyfbNuXscobJr2PTdbtmzBunXr8OWXX6rcVbdnzx6l54MHD0a7du2waNEijBkzBq6urhpEXzAF6fu4ceOUnnfv3h3+/v7o378/VqxYgcmTJyva0Of3XS6XY9u2bahduzaqVaumUn727Fml50OGDEHdunXx1VdfYdCgQSo3eRR1+X3d3vfZkJiYWIhR6sYff/yBGTNmoFevXmjevLlS2bp165Sef/zxxwgMDMSaNWswYcIENGjQQJehFkivXr3Qq1cvxfMuXbqgTZs2aNKkCWbPno1Vq1YBMIz3fMuWLXBwcECrVq1UyoriZ/yb7ty5g6CgIHzwwQf45JNP3llXqr9zg79MaWFhobgG/KacMS55fWHkbM/PvhYWFkhPT8+1ndTUVMm+lDTt+9tOnz6NoUOHok2bNpg9e/Z768tkMkyYMAGZmZka3dmkDdrqe45+/frByckJR48eVTqGPr/vp06dQnR0dK6/lHNjamqK0aNHIz4+HleuXMl/wEVEfl+39302FLck9G137txB165d4ePjg7Vr1+Zrn5y77N78+yiuGjVqhPr166v8rQP6+57fv38f58+fR+/evZXGR+WlKHzG53jy5AkCAgJga2urGC/7LlL9nRt8Mubs7JzrXDk52/K6HdvOzg5mZmb52tfZ2RlZWVmIjY1Vqpeeno7nz59Ldsu3pn1/019//YVOnTrBx8cHYWFh+fpDBaA4vf/ixQs1ItYebfT9bW5ubkr90ef3Hcg+7W9kZIS+ffvm+9hSv+8Fkd/XLeeyRV51pXrftSEyMhKtW7eGra0tDhw4gJIlS+Zrv+L8vucmt791QD/fcwCKsZ75/eEFFI33PCEhAe3atUN8fDx+//33fL0PUv2dG3wy5uvri4iICJVTihcuXFCU58bIyAg1atTA5cuXVcouXLiAihUrKj6octp4u+7ly5chl8vzPEZh07TvOe7du4e2bdvC0dERBw4cUJpz631yJsZ1cHBQL2gtKWjf3yaEwMOHD5X6o6/vO5D9a/DXX39Fs2bN1PrQkfp9LwhfX19cvXoVcrlcafuFCxdgaWmpmFMpr/f98ePHiIqKkux9L6jnz5+jdevWSEtLw6FDh3IdK5OX4vy+5+b+/fv5+lsv7u95ji1btsDT01OtS8xSv+epqano2LEjIiIisG/fPnh7e+drP8n+ztWaCEMP/fnnnypzLqWmpopKlSqJ+vXrK7Y9evRIZd6QefPmqcxBc+fOHWFsbCwmTZqk2JacnCzs7OxEhw4dlPYfMGCAsLS0FM+fP9d2t/KlIH2PiYkRFStWFC4uLuLBgwd5HuP58+ciMzNTaVt6err48MMPhampqWIOHF0rSN9jY2NV2gsJCREAxKJFixTb9PF9z7Fz504BQKxbty7X8txeo8TEROHp6Sns7e1V5gCSwrvmXXr8+LEIDw8X6enpim3btm1TmX/o2bNnolSpUqJ37/9r515Cquj/OI5/x44eMy+FFJJ5CY2ioiNFEV0QgqIijUyNWpR0OSB5obIMMkwkoygE5YQZUaCb2kiIEUJEi2hjQZsojC60CZLKKHCh5/MsQv/PefT4LzHH9P2CWcyc38z8vsztM+M4e0LmX7JkiXw+X8i+X1VVJcdx9OLFi/Ev5jf9bu3fv3/XmjVrFBcXp66urrDL7e3tVV9fX8i0YDCoPXv2yMz09OnTcathLH637pH2446ODpmZysrKQqZPtW0+6NmzZzIznT17dsTlTsZzfH9/v3Jzc+XxeNTR0RG23WQ6zqd9GJOkgoICeTwenTx5UteuXdO6devk8Xj06NGjoTbZ2dn6b3YdvLjMmzdPly5dUn19vVJSUjR//vxhB/HgxTo/P1/Xr1/X/v37ZWY6f/78hNQYzlhr9/l8MjOdOnVKLS0tIUNnZ+dQu5s3byojI0OVlZVqampSXV2dli9fLjNTXV3dhNU5krHWPnPmTBUVFenKlSsKBALau3evHMdRVlaWfvz4EdJ2qm33Qbt375bX69XXr19H/L26ulo+n09VVVVqbm5WTU2N0tLS5DiOWltb/0hNv6qxsVG1tbUqLi6WmSkvL0+1tbWqra0dqufAgQMys5Abjf7+fq1du1axsbGqqalRIBDQsmXLFBcXp5cvX4aso729XY7jaNOmTWpublZZWZkiIiJ05MiRiSx1mLHWvnPnTpmZDh48OOx4b2trG2r38OFDJSUl6dixYwoEArp8+bLWr18vM5Pf75/gav9nrHVnZmaqoKBAFy9eVFNTk/x+vzwej1JSUvTx48eQdUy1bT7oxIkTMrNh+/igyXiOLy8vl5kpJydn2P7a0tIy1G4yHeeEMf38+nhFRYWSkpLk9Xq1evVq3b9/P6RNuAvThw8flJ+fr/j4eMXGxmrHjh3q7u4ecT3Nzc1avHixoqKilJGRofr6egWDwT9S068aa+1mFnbIzs4eatfV1aWcnBwlJycrKipKsbGx2rBhg+7cuTMR5Y1qrLUfPnxYS5cuVVxcnCIjI5WZmanKykp9+/ZtxPVMpe0u/Xz6ER0drby8vLDL7+zs1ObNm5WUlKTIyEjNnj1bW7Zs0YMHD8a9lt+VlpYWdt8dPCmHuzh9/vxZhw4dUmJiomJiYpSdnR326+xtbW3KysqS1+vVggULVFVVNeKTh4k01tpHmy8tLW2o3Zs3b1RQUKD09HRFR0crJiZGq1atUlNTk6v7/FjrPnPmjLKyspSQkKDIyEilpqaquLh4WBAbNJW2uSQNDAwoOTlZK1euDLv8yXiOHzx3hRsGTabj3JHCfD4cAAAAf9y0f4EfAADATYQxAAAAFxHGAAAAXEQYAwAAcBFhDAAAwEWEMQAAABcRxgAAAFxEGAMAAHARYQwAAMBFhDEA09q5c+fMcRzr6ekZtV16eroVFRVNTKcATCuEMQAAABd53O4AAPwNXr16ZRER3L8CGH+EMQD4BV6v1+0uAJiiuM0DADPr6emxwsJCi4+Pt8TERCsvL7e+vr6h3//7ztitW7fMcRx7/PixHT9+3ObOnWuzZs2yXbt22adPn1yoAMDfijAGAGZWWFhofX19duHCBdu+fbs1NDSY3+//v/OVlpba8+fPrbq62oqLi629vd1KSkomoMcApgr+TAkAZrZw4UK7e/eumZkdPXrU4uPj7erVq1ZRUWErVqwIO19iYqJ1dnaa4zhmZhYMBq2hocF6e3stISFhQvoO4O/GkzEAsJ8B7N9KS0vNzOzevXujzuf3+4eCmJnZxo0bbWBgwN6/fz/+nQQwJRHGAMDMFi1aFDKekZFhERER9u7du1HnS01NDRmfM2eOmZl9+fJlXPsHYOoijAHACP79tGs0M2bMGHG6pPHsDoApjDAGAGbW3d0dMv769WsLBoOWnp7uTocATBuEMQAws0AgEDLe2NhoZmbbtm1zozsAphH+mxIAzOzt27eWm5trW7dutSdPnlhra6vt27fPfD6f210DMMXxZAwAzOz27dvm9Xrt9OnT1tHRYSUlJXbjxg23uwVgGnDEW6YAAACu4ckYAACAiwhjAAAALiKMAQAAuIgwBgAA4CLCGAAAgIsIYwAAAC4ijAEAALiIMAYAAOAiwhgAAICLCGMAAAAuIowBAAC4iDAGAADgon8AtBy9tPOz02QAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 600x600 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"import json\n",
"import cabinetry\n",
"import pyhf\n",
"import numpy as np\n",
"\n",
"# parse model spec\n",
"modelspec = json.load(open(\"simple-model.json\"))\n",
"\n",
"# construct a workspace from a background-only model and a signal hypothesis\n",
"workspace = pyhf.Workspace(modelspec)\n",
"\n",
"# construct the probability model and observations\n",
"model, data = cabinetry.model_utils.model_and_data(workspace)\n",
"\n",
"# attempt to change model parameters\n",
"params = {\"mu\": 5, \"uncorr_bkguncrt[0]\": 1}\n",
"\n",
"\n",
"def model_prediction_without_uncertainties(model, params):\n",
" # start with default pre-fit parameter values for everything\n",
" param_values = cabinetry.model_utils.asimov_parameters(model)\n",
"\n",
" # update central values of parameters as specified in params dict\n",
" labels = model.config.par_names\n",
" for k, v in params.items():\n",
" idx = cabinetry.model_utils._parameter_index(k, labels) # index of parameter\n",
" param_values[idx] = v\n",
"\n",
" # get model prediction per bin and slice into desired format\n",
" yields_combined = pyhf.tensorlib.to_numpy(\n",
" model.main_model.expected_data(param_values, return_by_sample=True)\n",
" )\n",
" model_yields = [\n",
" yields_combined[:, model.config.channel_slices[ch]].tolist()\n",
" for ch in model.config.channels\n",
" ]\n",
"\n",
" # need placeholders for total model uncertainty per bin and channel (set all to zero)\n",
" total_stdev_model_bins = [np.zeros(shape=(len(model.config.samples) + 1, model.config.channel_nbins[ch])).tolist() for ch in model.config.channels]\n",
" total_stdev_model_channels = [np.zeros(shape=(len(model.config.samples) + 1)).tolist() for ch in model.config.channels]\n",
" return cabinetry.model_utils.ModelPrediction(model, model_yields, total_stdev_model_bins, total_stdev_model_channels, \"custom_parameters\")\n",
" \n",
"model_prediction = model_prediction_without_uncertainties(model, params)\n",
"cabinetry.visualize.data_mc(model_prediction, data)"
]
},
{
"cell_type": "markdown",
"id": "cdc0142b",
"metadata": {},
"source": [
"## with proper uncertainties"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "73d778db",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[{'figure': <Figure size 600x600 with 2 Axes>, 'region': 'singlechannel'}]"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAmMAAAJjCAYAAABAyBZ9AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACjE0lEQVR4nOzdd1hTVx8H8G9ACEtQFAQBBVGruHBgrXXVgQOxal1176qoVWuHVYujrrqrqHUUrdtStc66t9aJW4tbQRRRhqwwct8/8iYhBpDEkJDw/TxPHs295957Tg6555dz7z1HJAiCACIiIiIyCDNDZ4CIiIioKGMwRkRERGRADMaIiIiIDIjBGBEREZEBMRgjIiIiMiAGY0REREQGxGCMiIiIyICKGToDpkAqleL58+coXrw4RCKRobNDREREWhIEAW/fvkXZsmVhZqafPisGYzrw+PFjeHt7GzobREREpCMPHjxAhQoV9HIsBmM6YGFhAQC4ffs23Nzc8rWNn58fLl68qNFxNN2moI+RmJgIDw8PPHv2DPb29gWWL318VppuY0pl1zS9NmUvjOXQdBt91bk225hK2QvjZ1VU/96Bolv2qKgo+Pj4KNp2fWAwpgPyS5PFixfP9x+subm5Ric1bbbRxzEAwN7evkDzpa9yFNWya3MMQLOyF9ZyFMY612YbUyl7Yf2sgKL79w4UvbInJiYCgF5vO+IN/AYSFBRU4Nvo4xjaKKzlKKplL4zl1mYbU6lzbbYxlbIX1s9KU4W1HCx74SXiROEfLjIyUtGV6+7ubujs6E1iYiIcHByQkJCg1a8tY8ayF72yF9VyAyw7y160ym6INp09YzogFotV/i0qxGIxgoODi1y5AZa9KJa9qJYbYNlZ9qJVdkO06ewZ04Gi+uuBiIjI1BiiTTeKnrGkpCQEBwejTZs2cHR0hEgkwtq1a3NMK5VKsXz5cvj6+sLa2hqlSpVC8+bNce3aNbV0v/zyC7y8vGBlZYWaNWti8+bNeigNERERkZJRBGOxsbGYNm0a7ty5g1q1auWZduDAgRg9ejTq1q2LJUuW4KeffkK5cuUQExOjkm7ixIn4/vvv0apVKyxZsgTlypVDz549sWXLloIsChEREZEKo7hMKZFIEBcXBxcXF1y6dAl+fn4IDQ1F//79VdJt27YN3bt3x/bt29GpU6dc9xcVFQUvLy8MHToUS5cuBSAbcbdp06Z49OgRHj9+DHNz83znj5cpiYiITAMvU+ZCLBbDxcXlvekWLFiA+vXro1OnTpBKpUhOTs4x3d9//42MjAyMGDFCsUwkEmH48OGIjIzEuXPndJZ3IiIiorwYRTCWH4mJibhw4QL8/Pzw448/wsHBAXZ2dqhQoQK2bdumkjY8PBy2traoWrWqyvL69esr1mubh+wviUSiXWGIiIhILyQSiVr7rW8mE4w9ePAAgiBgy5Yt+P333/HLL79g48aNcHJyQo8ePfDPP/8o0kZHR6NMmTJqo+u6uroCAJ4/f65VHjw8PODg4KB4zZo1S/sCERERUYGbNWuWStvt4eGh9zyYzHRISUlJAIDXr1/j33//xccffwwA6NChA7y8vPDzzz+jTZs2AIDU1NQcxw+xsrJSrNfGu/N3FbWxWYiIiIzNhAkTMG7cOMV7+Zyc+mQywZi1tTUAwMvLSxGIAYCdnR0CAwOxYcMGZGZmolixYrC2ts7xEmJaWprKvjSl6bxtREREZFhisdjgnScmc5mybNmyAIAyZcqorXN2dkZGRobihn5XV1e8ePEC7z5IGh0drbIvIiIiooJmUsGYi4sLoqKi1NY9f/4cVlZWKF68OADA19cXKSkpuHPnjkq68+fPK9YXRv3794enp2eBH0ckEmHKlCkFsu/Hjx9DJBJh3rx5BbJ/fZOXJ7dBiImIcqOvc3p+eXp6qg0ZRfphMsEYAHTv3h3Pnj3DoUOHFMtiY2Px999/o3nz5jAzkxX3888/h4WFBZYtW6ZIJwgCVqxYATc3NzRs2FDveSciIqKiyWjuGVu6dCni4+MVTzru3r0bkZGRAIBRo0bBwcEBEyZMwLZt2/DFF19g3LhxcHBwwIoVK5CRkYGZM2cq9uXu7o4xY8Zg7ty5yMjIgJ+fH3bu3IlTp05h48aNGg34qk+rVq2CVCo1dDaIiAqcKHS8obMAABAGFFwvPs/pJGc0wdi8efPw5MkTxfvt27dj+/btAIDevXvDwcEBZcqUwenTpzF+/HgsXLgQGRkZ+OSTT7Bhwwa1aZRmz56NkiVL4rfffsPatWtRqVIlbNiwAT179tRruTRhYWFh6CwQEZGO8JxOckZzmfLx48cQBCHHV/Zr7hUqVMD27duRkJCAlJQUHDlyBH5+fmr7MzMzw4QJE/D48WNIJBLcvHkTvXr10mOJ1L19+xZjxoyBp6cnxGIxnJ2d0apVK1y5cgWA+v0F2e+/WrlyJby9vSEWi+Hn54eLFy+q7f/PP/+Ej48PrKysUL16dezYsSPf9yxERUVh4MCBKFOmDMRiMapVq4bff/9dLV1aWhqmTJmCypUrw8rKCq6urujcuTMePHiglvZ9eb5+/Tr69++PChUqwMrKCi4uLhg4cCBev36tkm7KlCkQiUS4f/8++vfvjxIlSsDBwQEDBgxASkqKSlqRSISRI0di586dqF69uqIs2ceh07TMREQ50fScDsiGZ+rTpw/s7e1RokQJ9OvXD9euXVO7N7V///6ws7NDVFQUOnbsCDs7Ozg5OWH8+PHIyspS2ee8efPQsGFDlCpVCtbW1qhbty7CwsIKuvikAaPpGSsKhg0bhrCwMIwcORI+Pj54/fo1Tp8+jTt37qBOnTq5brdp0ya8ffsWX331FUQiEX755Rd07twZDx8+VPzy2rt3L7p3744aNWpg1qxZiIuLw6BBg+Dm5vbefL18+RINGjRQBDJOTk7Yv38/Bg0ahMTERIwZMwYAkJWVhfbt2+PIkSPo0aMHvv76a7x9+xaHDh3CzZs34e3trVGeDx06hIcPH2LAgAFwcXHBrVu3sHLlSty6dQv//vuv2qC93bp1g5eXF2bNmoUrV65g9erVcHZ2xpw5c1TSnT59Gtu3b8eIESNQvHhx/Prrr/jiiy/w9OlTlCpVSqMyExHlRtNzulQqRWBgIC5cuIDhw4ejSpUq+Pvvv9GvX78c95+VlYXWrVvj448/xrx583D48GHMnz8f3t7eGD58uCLd4sWL0aFDB/Tq1Qvp6enYsmULunbtij179iAgIKDAyk/5x2CsENm7dy+GDBmC+fPnK5Z99913793u6dOnuHfvHkqWLAkA+Oijj/D555/jwIEDaN++PQDZoHZubm44c+YM7OzsAAAtWrRAs2bNUL58+Tz3P3HiRGRlZeHGjRuKYGXYsGH48ssvMWXKFHz11VewtrbGH3/8gSNHjmDBggUYO3asYvsffvhBbRiR/OR5xIgR+Oabb1S2a9CgAb788kucPn0ajRs3VllXu3ZtrFmzRvH+9evXWLNmjVowdufOHdy+fVsRHH722WeoVasWNm/ejJEjR2pUZiKi3Gh6Tt+5cyfOnTuHRYsW4euvvwYADB8+HK1atcoxfVpaGrp3747JkycDkJ2j6tSpgzVr1qgEYxERESrnq5EjR6JOnTpYsGABg7FCwmguUxYFJUqUwPnz5zWejql79+6KoAaAIkh5+PAhANnQHjdu3EDfvn0VgRgANG3aFDVq1Mhz34Ig4K+//kJgYCAEQUBsbKzi1bp1ayQkJCi63P/66y+ULl0ao0aNUtvPu71Y78szoDr4blpaGmJjY9GgQQMAUBwzu2HDhqm8b9y4MV6/fq02z1jLli1Veulq1qwJe3t7xbE1KTMRUW40Paf/888/sLCwwJAhQxTLzMzMEBQUlOs2OZ33sp9HAdVzaVxcHBISEtC4cWOexwoRBmM65OfnBx8fH4SEhGi1/S+//IKbN2/Cw8MD9evXx5QpU9S+VDkpV66cynt5kBMXFwcAigcfKlasqLZtTsuye/XqFeLj47Fy5Uo4OTmpvAYMGAAAiImJASCbH/Sjjz5CsWLv73B9X54B4M2bN/j6669RpkwZWFtbw8nJCV5eXgCAhIQErfaZUzp5Wnk6TcpMRJQbTc/pT548gaurK2xsbFSW53aetrKygpOTk8qy7OcyuT179qBBgwawsrKCo6MjnJycsHz58hzPo0VZSEgIfHx8crzPvKDxMqUOXbx48YOmQ+rWrRsaN26MHTt24ODBg5g7dy7mzJmD7du3o23btrlul9tQHO9eGtSG/LHr3r1753rfQs2aNTXeb37y3K1bN5w9exbffvstfH19YWdnB6lUijZt2uT4OHh+P4f3pSuoMhNR0aLtOT2/8jMM06lTp9ChQwc0adIEy5Ytg6urKywsLBAaGopNmzZ9cB5MSVBQEIKCgpCYmAgHBwe9HpvBWCHj6uqKESNGYMSIEYiJiUGdOnUwY8aMD/riyu8Ju3//vtq6nJZl5+TkhOLFiyMrKwstW7bMM623tzfOnz+PjIyMD35kOy4uDkeOHMHUqVPx008/KZbfu3fvg/abH5qUmYgoL5qc08uXL49jx44hJSVFpXfsfefpvPz111+wsrLCgQMHVOZfDA0N1XqfpHu8TFlIZGVlqXUZOzs7o2zZsjlOaq6JsmXLonr16vjjjz+QlJSkWH7ixAncuHEjz23Nzc3xxRdf4K+//sLNmzfV1r969Urx/y+++AKxsbFYunSpWjpNe+nkv/je3W7RokUa7UcbmpSZiCgn2pzTW7dujYyMDKxatUqxTCqVan3rCyA7n4lEIpXhLh4/foydO3dqvU/SPfaMFRJv376Fu7s7unTpglq1asHOzg6HDx/GxYsXVZ7E0dbMmTPx+eef49NPP8WAAQMQFxeHpUuXonr16ioBWk5mz56NY8eO4eOPP8aQIUPg4+ODN2/e4MqVKzh8+DDevHkDAOjbty/++OMPjBs3DhcuXEDjxo2RnJyMw4cPY8SIEfj888/znV97e3s0adIEv/zyCzIyMuDm5oaDBw/i0aNHH/Q55Fd+y0xElBNtzukdO3ZE/fr18c033+D+/fuoUqUKdu3apTjfvPsgVH4EBARgwYIFaNOmDXr27ImYmBiEhISgYsWKuH79+geVkXSHwVghYWNjgxEjRuDgwYPYvn07pFIpKlasiGXLlqk8oqytwMBAbN68GVOmTMEPP/yASpUqYe3atVi3bh1u3bqV57ZlypTBhQsXMG3aNGzfvh3Lli1DqVKlUK1aNZVhI8zNzbFv3z7MmDEDmzZtwl9//YVSpUqhUaNG731qMyebNm3CqFGjEBISAkEQ4O/vj/3796Ns2bIa70tT+S0zERWMgpyGSB+0Oaebm5tj7969+Prrr7Fu3TqYmZmhU6dOCA4OxqeffgorKyuN89G8eXOsWbMGs2fPxpgxY+Dl5YU5c+bg8ePHDMYKEZGgi7u8izj5zX4JCQkfdAO/Ifj6+sLJyUllcnUiIio8du7ciU6dOuH06dP49NNPDZ0dk2eINp33jBURGRkZyMzMVFl2/PhxXLt2Dc2aNTNMpoiISEVqaqrK+6ysLCxZsgT29vZ5zsRCxo2XKYuIqKgotGzZEr1790bZsmVx9+5drFixAi4uLmqDBhIRkWGMGjUKqamp+OSTTyCRSLB9+3acPXsWM2fO5KwfJozBWBFRsmRJ1K1bF6tXr8arV69ga2uLgIAAzJ49WzHdDxERGVbz5s0xf/587NmzB2lpaahYsSKWLFmimKqNTBPvGdMBY75njIiIiJR4zxgRERFREcNgTIc+dG5KIiIiMgxDzk3Jy5Q6wMuUREREpoGXKYmIiIiKGAZjRERERAbEYIyIiIjIgBiMERERERkQg7EiKiUlBVOmTMHx48cNnRXKw6ZNm7Bo0SJDZ4OItDBlyhSIRCLExsbq7ZjNmjVD9erV9Xa8wsbT0xP9+/c3dDY0xmCsiEpJScHUqVMZjBVyDMaIiEwfp0OiIi8tLQ2WlpYwMysav02kUinS09NhZWVl6KwQ5Ur0zW5DZwEAIMwPNHQWqAgoGq2PEYmKisKgQYNQtmxZiMVieHl5Yfjw4UhPT1d0eb9r7dq1EIlEePz4sWLZpUuX0Lp1a5QuXRrW1tbw8vLCwIEDAQCPHz+Gk5MTAGDq1KkQiUQQiUSYMmWKYvujR4+icePGsLW1RYkSJfD555/jzp07KseV5yciIgK9e/eGg4MDnJycMHnyZAiCgGfPnuHzzz+Hvb09XFxcMH/+fI0/D09PT7Rv3x4HDx6Er68vrKys4OPjg+3bt6uke/PmDcaPH48aNWrAzs4O9vb2aNu2La5du6aS7vjx4xCJRNiyZQsmTZoENzc32NjYIDExUeN9bNu2DVOnToWbmxuKFy+OLl26ICEhARKJBGPGjIGzszPs7OwwYMAASCQStbJt2LABdevWhbW1NRwdHdGjRw88e/ZMsb5Zs2bYu3cvnjx5oqgjT09PxXqJRILg4GBUrFgRYrEYHh4e+O6779SOJRKJMHLkSGzcuBHVqlWDWCzGP//8AwDYsmUL6tati+LFi8Pe3h41atTA4sWLNa4nIio6kpOTDZ0Fk8NgrBB5/vw56tevjy1btqB79+749ddf0adPH5w4cQIpKSn53k9MTAz8/f3x+PFj/PDDD1iyZAl69eqFf//9FwDg5OSE5cuXAwA6deqE9evXY/369ejcuTMA4PDhw2jdujViYmIwZcoUjBs3DmfPnsWnn36qEvDJde/eHVKpFLNnz8bHH3+Mn3/+GYsWLUKrVq3g5uaGOXPmoGLFihg/fjxOnjyp8edy7949dO/eHW3btsWsWbNQrFgxdO3aFYcOHVKkefjwIXbu3In27dtjwYIF+Pbbb3Hjxg00bdoUz58/V9vn9OnTsXfvXowfPx4zZ86EpaWlxvuYNWsWDhw4gB9++AEDBw7E9u3bMWzYMAwcOBARERGYMmUKOnfujLVr12LOnDkq286YMQN9+/ZFpUqVsGDBAowZMwZHjhxBkyZNEB8fDwCYOHEifH19Ubp0aUUdyS9ZSqVSdOjQAfPmzUNgYCCWLFmCjh07YuHChejevbtaXo8ePYqxY8eie/fuWLx4MTw9PXHo0CF8+eWXKFmyJObMmYPZs2ejWbNmOHPmjMZ1RES5i42NRbdu3WBvb49SpUrh66+/RlpammJ9aGgomjdvDmdnZ4jFYvj4+CjO0e/av38/mjZtqvgB5efnh02bNuV5/IMHD8LGxgZffvklMjMzAQCpqakYPXo0SpcujeLFi6NDhw6IiopS+2Eu/9F9+/Zt9OzZEyVLlkSjRo0AAJmZmZg+fTq8vb0hFovh6emJH3/8MccfhNn3Kffu/V3yjoUzZ85g3LhxcHJygq2tLTp16oRXr16pbCsIAn7++We4u7vDxsYGn332GW7dupXn51CoCfTBEhISBABCQkLCB+2nb9++gpmZmXDx4kW1dVKpVAgODhZyqrLQ0FABgPDo0SNBEARhx44dAoAc9yP36tUrAYAQHBysts7X11dwdnYWXr9+rVh27do1wczMTOjbt69imTw/Q4cOVSzLzMwU3N3dBZFIJMyePVuxPC4uTrC2thb69euX10egpnz58gIA4a+//lIsS0hIEFxdXYXatWsrlqWlpQlZWVkq2z569EgQi8XCtGnTFMuOHTsmABAqVKggpKSkqKTXdB/Vq1cX0tPTFcu//PJLQSQSCW3btlXZxyeffCKUL19e8f7x48eCubm5MGPGDJV0N27cEIoVK6ayPCAgQGVbufXr1wtmZmbCqVOnVJavWLFCACCcOXNGsQyAYGZmJty6dUsl7ddffy3Y29sLmZmZavsnMjSM21UoXh9Cfo6sUaOGEBgYKCxdulTo3bu3AEDo06ePIp2fn5/Qv39/YeHChcKSJUsEf39/AYCwdOlSlf2FhoYKIpFIqF69ujBjxgwhJCREGDx4sMq+mjZtKlSrVk3xfvfu3YJYLBb69u2r8l3v1q2bIh8hISFCt27dhFq1aqm1C/Iy+Pj4CJ9//rmwbNkyISQkRBAEQejXr58AQOjSpYsQEhIi9O3bVwAgdOzYUbUuc2lrypcvr9ImyNuy2rVrC82bNxeWLFkifPPNN4K5ubnQrVs3lW0nTZokABDatWsnLF26VBg4cKBQtmxZoXTp0hq3M+/SVZuuCd4zpkN+fn4wNzdHUFAQgoKCNNpWKpVi586dCAwMRL169dTW53R5MjclSpQAAOzZswe1atWChYVFvreNjo7G1atX8d1338HR0VGxvGbNmmjVqhX27dunts3gwYMV/zc3N0e9evUQGRmJQYMGqeTpo48+wsOHD/OdF7myZcuiU6dOivf29vbo27cv5syZgxcvXsDFxQVisVixPisrC/Hx8bCzs8NHH32EK1euqO2zX79+sLa2Vlmm6T769u2r8tl+/PHH2Lx5s+JycPblv/76KzIzM1GsWDFs374dUqkU3bp1U3nKysXFBZUqVcKxY8fw448/5vmZ/Pnnn6hatSqqVKmiso/mzZsDAI4dO4aGDRsqljdt2hQ+Pj4q+yhRogSSk5Nx6NAhtGnTJs/jEZH2vLy88PfffwMAgoKCYG9vj2XLlmH8+PGoWbMmTpw4oXI+GjlyJNq0aYMFCxYo2pKEhASMHj0a9evXx/Hjx1Xu+RRymdVw+/bt6NGjB/r3748VK1Yo7ou9cuUKtm3bhjFjxmDhwoUAgBEjRmDAgAFqt2XI1apVS6UH7tq1a1i3bh0GDx6MVatWKfbh7OyMefPm4dixY/jss8+0+rxKlSqFgwcPKto9qVSKX3/9FQkJCXBwcMCrV6/wyy+/ICAgALt371akmzhxImbOnKnVMQHZ3JQhISHIysrSeh/a4mVKHbp48SJu376tcSAGAK9evUJiYqJOHklu2rQpvvjiC0ydOhWlS5fG559/jtDQ0BzvW3rXkydPAAAfffSR2rqqVasiNjZW7X6BcuXKqbx3cHCAlZUVSpcurbY8Li5O0+KgYsWKasFo5cqVAUBx2VQqlWLhwoWoVKkSxGIxSpcuDScnJ1y/fh0JCQlq+/Ty8lJbpuk+cio3AHh4eKgtl0qlin3cu3cPgiCgUqVKcHJyUnnduXMHMTEx7/1M7t27h1u3bqltL/9c3t1HTuUdMWIEKleujLZt28Ld3R0DBw5U3EtGRLrzbpswatQoAFD8uM0eiCUkJCA2NhZNmzbFw4cPFeeNQ4cO4e3bt/jhhx/UHr7J6cf65s2b0b17d3z11Vf47bffVB5Qkn/PR4wYkWO+cjJs2DCV9/K8jxs3TmX5N998AwDYu3dvrvt6n6FDh6qUqXHjxsjKylK0T4cPH0Z6ejpGjRqlkm7MmDFaHxOQ1dPt27dx8eLFD9qPNtgzZkRy6x17N4oXiUQICwvDv//+i927d+PAgQMYOHAg5s+fj3///Rd2dnY6zZe5uXm+lgG5/4L7UDNnzsTkyZMxcOBATJ8+HY6OjjAzM8OYMWMglUrV0r/bK6bNPnIr4/vKLpVKIRKJsH///hzT5qd+pFIpatSogQULFuS4/t2AMKfyOjs74+rVqzhw4AD279+P/fv3IzQ0FH379sW6devemwciyp9KlSqpvPf29oaZmZnix+SZM2cQHByMc+fOqd0fLO8NevDgAQDk6wf7o0eP0Lt3b3Tt2hVLlixRW//kyROYmZmp/UirWLFirvt8N618H+9u4+LighIlSigCJ228+0O3ZMmSAKD4MS/f97ufq5OTkyKtsWEwVkg4OTnB3t4eN2/ezDWN/I8sPj5ecSkSQK5/9A0aNECDBg0wY8YMbNq0Cb169cKWLVswePDgXAO78uXLAwD+++8/tXV3795F6dKlYWtrm99i6cT9+/chCIJKniMiIgBA8XRhWFgYPvvsM6xZs0Zl2/j4eLUeutzoYh/54e3tDUEQ4OXlpejJyk1u9eTt7Y1r166hRYsWGl3CfpelpSUCAwMRGBgIqVSKESNG4LfffsPkyZPzPDETkfayf2cfPHiAFi1aoEqVKliwYAE8PDxgaWmJffv2YeHChTn+EHwfV1dXuLq6Yt++fbh06VKOt75oKqcfdIBmt9C8K7fLgfr+MV8Y8DJlIWFmZoaOHTti9+7duHTpktp6QRDg7e0NACpPJCYnJ6v1YsTFxan90fr6+gKA4lKljY0NACie3JNzdXWFr68v1q1bp7Lu5s2bOHjwINq1a6dV+T7E8+fPsWPHDsX7xMRE/PHHH/D19YWLiwsA2Zf33TL/+eefiIqKyvdxdLGP/OjcuTPMzc0xdepUteMJgoDXr18r3tva2uZ4ibRbt26IiopS3KuRXWpqar4ePc9+HED2N1izZk0AyNclbSLKn3v37qm8v3//PqRSKTw9PbF7925IJBLs2rULX331Fdq1a4eWLVuqBT/y839eP9jlrKyssGfPHlSqVAlt2rRRe8qwfPnykEqlePTokVq+8ku+j3fL9vLlS8THxyt+2AOyjoR325r09HRER0fn+3jvHhtQ/1xfvXql1a0whQF7xgqRmTNn4uDBg2jatCmGDh2KqlWrIjo6Gn/++SdOnz4Nf39/lCtXDoMGDcK3334Lc3Nz/P7773BycsLTp08V+1m3bh2WLVuGTp06wdvbG2/fvsWqVatgb2+vCKasra3h4+ODrVu3onLlynB0dET16tVRvXp1zJ07F23btsUnn3yCQYMGITU1FUuWLIGDg0OOjycXtMqVK2PQoEG4ePEiypQpg99//x0vX75EaGioIk379u0xbdo0DBgwAA0bNsSNGzewceNGVKhQId/H0cU+8sPb2xs///wzJkyYgMePH6Njx44oXrw4Hj16hB07dmDo0KEYP348AKBu3brYunUrxo0bBz8/P9jZ2SEwMBB9+vTBtm3bMGzYMBw7dgyffvopsrKycPfuXWzbtg0HDhx476/hwYMH482bN2jevDnc3d3x5MkTLFmyBL6+vqhatapOy0xUlIWEhMDf31/xXn7psG3btoof19l/mCUkJKic3wDA398fxYsXx6xZs9CmTRu1G/jf7aFycHDAgQMH0KRJE7Rq1QqnTp1SBHStW7fGxIkTsWzZMsUN/NnzlR/t2rXDjz/+iEWLFuG3335TLJffOhEQEKBY5u3trTas0cqVK7W+Ub5ly5awsLDAkiVL4O/vryi7Mc9WYhTBWFJSEubOnYvz58/jwoULiIuLQ2hoaJ7zT2VkZKBWrVq4c+cO5s6dq2jc5KRSKebNm4fly5cjOjoalStXxoQJE/Dll18WcGly5+bmhvPnz2Py5MnYuHEjEhMT4ebmhrZt28LGxgYWFhbYsWMHRowYgcmTJ8PFxQVjxoxByZIlMWDAAMV+mjZtigsXLmDLli14+fIlHBwcUL9+fWzcuFHluv/q1asxatQojB07Funp6QgODkb16tXRsmVL/PPPPwgODsZPP/0ECwsLNG3aFHPmzMnxRvCCVqlSJSxZsgTffvst/vvvP3h5eWHr1q1o3bq1Is2PP/6I5ORkbNq0CVu3bkWdOnWwd+9e/PDDD/k+ji72kV8//PADKleujIULF2Lq1KkAZPd5+fv7o0OHDop0I0aMwNWrVxEaGoqFCxeifPnyCAwMhJmZGXbu3ImFCxfijz/+wI4dO2BjY4MKFSrg66+/fu/lTwDo3bs3Vq5ciWXLliE+Ph4uLi7o3r07pkyZUmRmIyDSh0ePHqFDhw5o06YNzp07hw0bNqBnz56oVasWrKysFLcLfPXVV0hKSsKqVavg7Oys0nNkb2+PhQsXYvDgwfDz81OM+XXt2jWkpKTkeJ9n6dKlcejQITRq1AgtW7bE6dOn4ebmhrp16+KLL77AokWL8Pr1azRo0AAnTpxQ3P6Rn0uPtWrVQr9+/bBy5UrEx8cr2p1169ahY8eOKk9SDh48GMOGDcMXX3yBVq1a4dq1azhw4IDWt384OTlh/PjxmDVrFtq3b4927dohPDwc+/fv1+ktJXqlt0E0PsCjR48EAEK5cuWEZs2aCQCE0NDQPLeZP3++YGtrKwAQ5s6dq7b+hx9+EAAIQ4YMEVauXCkEBAQIAITNmzdrnD9DjElSVJQvX14ICAgwdDaISM8MPb6YLscZu337ttClSxehePHiQsmSJYWRI0cKqampinS7du0SatasKVhZWQmenp7CnDlzhN9//11l/MjsaRs2bChYW1sL9vb2Qv369VXarXfHGRMEQbh//77g6uoqVK1aVXj16pUgCIKQnJwsBAUFCY6OjoKdnZ3QsWNH4b///hMAqIwRKS+DfLvsMjIyhKlTpwpeXl6ChYWF4OHhIUyYMEFIS0tTSZeVlSV8//33QunSpQUbGxuhdevWwv3793MdZ+zdMTLlYzseO3ZMZZ9Tp04VXF1dBWtra6FZs2bCzZs31fapDUO06UYRjKWlpQnR0dGCIAjCxYsX3xuMvXz5UnBwcBCmTZuWYzAWGRkpWFhYCEFBQYplUqlUaNy4seDu7q7xIJgMxgoOgzEiIv0IDw8XAAgbNmwwdFYMyhBtulFcixCLxYobtfPjhx9+wEcffYTevXvnuP7vv/9GRkaGyhgrIpEIw4cPR2RkJM6dO/fBeaa8vXr1Ci9evMj19ebNG0NnkYjIZKWmpqotW7RoEczMzNCkSRMD5KhoM4p7xjQhv2Z9+vTpXK97h4eHw9bWVu0m5fr16yvWy+fe0kRiYqLKe7FYrDKqOyn5+fnlOQ5N06ZNcfz4cf1liIioCPnll19w+fJlfPbZZyhWrJhirMGhQ4eqjVNo6iQSicoT5O+25fpgUsGYIAgYNWoUunfvjk8++STHSa0B2ZQ/ZcqUUQvWXF1dASDHSaHz490/4ODgYIM8fWgMNm7cmOMvMzn5mGq51SEREWmvYcOGOHToEKZPn46kpCSUK1cOU6ZMwcSJEw2dNb2bNWuW4kEqQzGpYGzt2rW4ceMGwsLC8kyXmpqaY4+V/FHhvIKEvDx79gz29vaK9+wVy92nn35q6CwQERVZrVq1QqtWrQydjUJhwoQJKtM6JSYm6r130GSCscTEREyYMAHffvvtez9Ea2vrHAe1TEtLU6zXhr29vUowRkRERIVbYbilyChu4M+PefPmIT09Hd27d8fjx4/x+PFjREZGApCNSP/48WOkp6cDkF2OfPHihdro5/IxXcqWLavfzBMREVGRZTLB2NOnTxEXF4dq1arBy8sLXl5eaNy4MQDZyPZeXl64ffs2ANnUQCkpKbhz547KPs6fP69YT0RERKQPJnOZcvTo0ejYsaPKspiYGHz11Vfo378/Pv/8c8Xo8Z9//jnGjh2LZcuWYenSpQBkN/+vWLECbm5uaNiwob6zT0REREWU0QRjS5cuRXx8vOJJx927dysuQ44aNQp16tRBnTp1VLaRP4lXrVo1lUDN3d0dY8aMwdy5c5GRkQE/Pz/s3LkTp06dwsaNG3OdMZ6IiIhI14wmGJs3b57KuFTbt2/H9u3bAcjm2HNwcNBof7Nnz0bJkiXx22+/Ye3atahUqZJivjAiIiIifREJ797FThpLTEyEg4MDEhIS+DQlERGRETNEm24yN/ATERERGSMGY0RERFQgmjVrhmbNmhk6G4UegzEiIjIq4eHhmDp1Kvbs2aM2XuS7IiIi8PPPP2Pr1q3IysrKM21UVBRmz56N33//PceBwTU1ZcoUiEQixMbG5ri+evXqhTpQOXv2LKZMmYL4+HiDHH/Tpk1YtGiRQY6tbwzGdMjPzw8+Pj4ICQkxdFaIiExSeHg4du3ahbp16yIgIEBtjuHsIiIisG3bNlSqVAldunTJ80n5qKgorF+/Hs7OzujVq5fBR2QvDM6ePYupU6d+UDB28OBBHDx4UKtt9R2MhYSEwMfHB35+fno7ppzRPE1pDC5evMgb+ImICggDMf1ITk6Gra2tTvZlaWmpk/3oQ1BQEIKCghQ38OsTe8aIiKjQM/VA7Pjx4xCJRNi2bRtmzJgBd3d3WFlZoUWLFrh//75a+vPnz6Ndu3YoWbIkbG1tUbNmTSxevFglzd27d9GlSxc4OjrCysoK9erVw65du1TSrF27FiKRCCdOnMCIESPg7OwMd3d3TJkyBd9++y0AwMvLCyKRCCKRSDF+Z2hoKJo3bw5nZ2eIxWL4+Phg+fLlavl8956x/JazWbNm2Lt3L548eaI4tqenJ5KSkmBra4uvv/5a7ViRkZEwNzfHrFmz8v25FxbsGSMiokLN1AOx7GbPng0zMzOMHz8eCQkJ+OWXX9CrVy/FdH0AcOjQIbRv3x6urq74+uuv4eLigjt37mDPnj2KIOXWrVv49NNP4ebmhh9++AG2trbYtm0bOnbsiL/++gudOnVSOe6IESPg5OSEn376CcnJyWjbti0iIiKwefNmLFy4EKVLlwYAODk5AQCWL1+OatWqoUOHDihWrBh2796NESNGQCqVIigo6IPLOXHiRCQkJCAyMhILFy4EANjZ2cHOzg6dOnXC1q1bsWDBApX63bx5MwRBQK9evT6gBgxEoA+WkJAgABASEhIMnRUiIpNy5coVYcqUKcLu3bsFqVSaZ9r//vtPmD59urBlyxYhMzMzz7SRkZHCrFmzhDVr1ghpaWm6zLJCcHCwAEB49epVjuurVasmNG3aVBAEQTh27JgAQKhataogkUgUaRYvXiwAEG7cuCEIgiBkZmYKXl5eQvny5YW4uDiV/WX/fFq0aCHUqFFDpWxSqVRo2LChUKlSJcWy0NBQAYDQqFEjtc9s7ty5AgDh0aNHanlPSUlRW9a6dWuhQoUKKsuaNm2qKKMm5RQEQQgICBDKly+vdpwDBw4IAIT9+/erLK9Zs6bKsbRliDadlymJiKhQKko9YnIDBgxQuc+qcePGAICHDx8CkH0mjx49wpgxY1CiRAmVbeWfz5s3b3D06FF069YNb9++RWxsLGJjY/H69Wu0bt0a9+7dQ1RUlMq2Q4YM0WgqQGtra8X/ExISEBsbi6ZNm+Lhw4dISEj44HLmpWXLlihbtiw2btyoWHbz5k1cv34dvXv3zncZChNepiQiokKnMARir1+/RqlSpT6oHO/zbrnKlSun8r5kyZIAgLi4OADAgwcPAMiGxcjN/fv3IQgCJk+ejMmTJ+eYJiYmBm5ubor3Xl5eGuX7zJkzCA4Oxrlz55CSkqKyLiEh4b03wL+vnHkxMzNDr169sHz5cqSkpMDGxgYbN26ElZUVunbtqlE5CgsGY0REVOgUhkBs3bp1GDdunNZlsLKyAgCkpqbmuD4lJUWRRi63/AsazFwolUoBAOPHj0fr1q1zTFOxYkWV99l7ut7nwYMHaNGiBapUqYIFCxbAw8MDlpaW2LdvHxYuXKg4fl4+tJx9+/bF3LlzsXPnTnz55ZfYtGkT2rdvr/enIHWFwRgRERU6hSEQ+9DLl+XLlwcA/Pfff/Dw8FBZl5KSgmfPnsHf31+jfXp7ewOQXZZr2bJljmkqVKgAALCwsMg1TX7k9tnv3r0bEokEu3btUunhOnbsmNbH0uT4gKxnsHbt2ti4cSPc3d3x9OlTLFmyRKfH1yfeM0ZERIVOYQjE+vXr90FlaNGiBSwtLbF8+XK13qKVK1ciMzMTbdu21WifderUgZeXFxYtWqQ2GKu8V8nZ2RnNmjXDb7/9hujoaLV9vHr1Kl/Hko819u5x5J919l6shIQEhIaG5rcY+T5+Xvef9enTBwcPHsSiRYtQqlQpjT/LwoQ9Y0REVOgUhkDMzs7ug8rg7OyMn376CZMmTUKTJk3QoUMH2NjY4OzZs9i8eTP8/f0RGBio0T7NzMywfPlyBAYGwtfXFwMGDICrqyvu3r2LW7du4cCBAwBko8k3atQINWrUwJAhQ1ChQgW8fPkS586dQ2RkJK5du/beY9WtWxeAbJiJHj16wMLCAoGBgfD394elpSUCAwPx1VdfISkpCatWrYKzs3OOwZ+26tati61bt2LcuHHw8/ODnZ2dyufVs2dPfPfdd9ixYweGDx8OCwsLnR1b7/T23KYJ49AWRESUmw0bNggNGjQQbG1tBbFYLFSpUkWYOnWqyrAT8iEf/vzzT5VtHz16JAAQQkNDVZafPn1aaNWqlVC8eHHB1tZWqFmzprBkyRKVNA8ePBD69u0ruLi4CBYWFoKbm5vQvn17ISwsTJFGPrTFxYsXc8z79OnTBTc3N8HMzExlmItdu3YJNWvWFKysrARPT09hzpw5wu+//642FEZuQ1vkp5xJSUlCz549hRIlSggAchzmol27dgIA4ezZsznmXxuGaNNFgqDBXYGUI/nUCZUrV4a5ubliSgUiIiIqOJ06dcKNGzdynKVAUyEhIQgJCUFWVhYiIiKQkJCgtykOGYzpgDwY02fFERERFWXR0dEoX748Jk6ciODgYJ3t1xBtOu8ZIyIiIqPx6NEjnDlzBqtXr4aFhQW++uorQ2fpg/FpSiIiIjIaJ06cQJ8+ffDo0SOsW7cOLi4uhs7SB+NlSh3gZUoiIiLTYIg2nT1jRERERAbEYIyIiIjIgBiMERERERkQgzEiIiIiA2IwRkRERGRADMaIiIiIDIjBGBERmZTk5GSIRCKIRCIkJycbOjtE78VgTIf8/Pzg4+ODkJAQQ2eFiIiINBASEgIfHx/4+fnp/dgc9FUHOOgrEVHhkZycDDs7OwBAUlISbG1tDZKPtWvXYsCAAYr3YrEYjo6OqFGjBgICAjBgwAAUL15c4/2ePXsWBw8exJgxY1CiRAkd5pgADvpKRERkcqZNm4b169dj+fLlGDVqFABgzJgxqFGjBq5fv67x/s6ePYupU6ciPj5exzklQ+FE4UREZFKysrIU/z958iT8/f1hbm5usPy0bdsW9erVU7yfMGECjh49ivbt26NDhw64c+cOrK2tDZY/Mjz2jBERkcnYvn07fHx8FO/btWsHT09PbN++3YC5Ute8eXNMnjwZT548wYYNGwAA169fR//+/VGhQgVYWVnBxcUFAwcOxOvXrxXbTZkyBd9++y0AwMvLS/GgwuPHjwEAoaGhaN68OZydnSEWi+Hj44Ply5frvXykGfaMERGRSdi+fTu6dOmCd2+FjoqKQpcuXRAWFobOnTsbKHfq+vTpgx9//BEHDx7EkCFDcOjQITx8+BADBgyAi4sLbt26hZUrV+LWrVv4999/IRKJ0LlzZ0RERGDz5s1YuHAhSpcuDQBwcnICACxfvhzVqlVDhw4dUKxYMezevRsjRoyAVCpFUFCQIYtLeRGMwNu3b4WffvpJaN26tVCyZEkBgBAaGqqSJisrSwgNDRUCAwMFd3d3wcbGRqhWrZowffp0ITU1Ncf9rl69WqhSpYogFouFihUrCr/++qtW+UtISBAACAkJCVptT0REHyYzM1Nwd3cXAOT4EolEgoeHh5CZmam3PIWGhgoAhIsXL+aaxsHBQahdu7YgCIKQkpKitn7z5s0CAOHkyZOKZXPnzhUACI8ePVJLn9M+WrduLVSoUEGLEhRNhmjTjeIyZWxsLKZNm4Y7d+6gVq1aOaZJSUnBgAED8OrVKwwbNgyLFi1C/fr1ERwcjLZt26r9Uvrtt98wePBgVKtWDUuWLMEnn3yC0aNHY86cOfooEhEZoQULFiAkJARJSUl5pktKSkJISAgWLFigcokpJxKJBL///jtmz56NqKioPNNmZWVh69at+PnnnxEREZFnWkEQsGfPHkydOhXh4eF5pgVk91ZNnToVJ0+efG/a8PBwTJ06FXv27FE7t74rIiICP//8M7Zu3apyL1dOoqKiMHv2bPz++++QSCR5pn39+rVKfZw6dQqRkZG5phcEAc+ePcOePXvy3K++68POzg5v374FAJX7xg4dOoTvvvsOIpEIAHDlypVc9529PqysrBTLExISEBsbi6ZNm+Lhw4dISEjQW33kxRi/HwXNKC5Turq6Ijo6Gi4uLrh06VKOY4BYWlrizJkzaNiwoWLZkCFD4OnpieDgYBw5cgQtW7YEAKSmpmLixIkICAhAWFiYIq1UKsX06dMxdOhQlCxZUj+FIyKjIRaL0a9fP8WwCTlJSkrCunXrIJFI0K9fP5QqVSrXtBKJBBs3bkRMTAz69OkDNze3XNNmZWUhLCwM9+7dQ7du3VC5cuVc0wqCgL179+Ly5cvo0KEDateunWe5Tp48iWPHjuGzzz5DkyZN8kwbHh6OXbt2oW7duggICFAECzmJiIjAtm3bUKlSJXTp0iXPm+ijoqKwfv16ODs7o1evXhCLxbmmff36NdatW6dSH9HR0XnmWy4lJSXXdYaoj6SkJDg7OwMA3rx5g6lTp2L9+vWIi4tT2UdCQkKO+363Ps6ePYvg4GCcO3dOrazh4eE4ffq0XuojN8b4/dAHo+gZE4vFcHFxyTONpaWlSiAm16lTJwDAnTt3FMuOHTuG169fY8SIESppg4KCkJycjL179+og10RkakyxoTGFQAyQ/WjPj9zSFVR9AMDp06dzrI/IyEgkJCSgYsWKAIBu3brht99+Q82aNTF9+nQcPHgQ//zzDwBAKpWq7ffWrVsq9fHw4UO0aNECsbGxWLBgAfbu3YtDhw5h7NixAIDdu3czECuEgRhgJD1jH+LFixcAoLjJEYCiSzL7o8YAULduXZiZmSE8PBy9e/fWXyaJyCiYWkNjKoEYADRu3Bju7u6IiorK8dKpSCSCu7s7GjdurLauIOsDkHUGDB8+XK0+1q9fDwBo3bo14uLicOTIETRr1gxTp05V1Me9e/dyLAsgu5TZsmVLRX3s3r0bEokEu3btQrly5RTp//rrLwCAp6cnA7FCGIgBRtIz9iF++eUX2Nvbo23btopl0dHRMDc3V3QNy1laWqJUqVJ4/vy5VsdKTExUeb3vGjsRmQZjbGhMKRADAHNzcyxevBgA1PInf79o0SK1/BRkfVy7dg0A0LRpU7X6OHr0KKZPnw4vLy/06tUL//77LwDZcBXZ62PRokVq+37z5g0AoEKFCir1IS9b9mD08uXL2LhxIwDZMB8MxNS/HxKJRK391jeT7hmbOXMmDh8+jGXLlqlMGZGamgpLS8sct7GyskJqaqpWx/Pw8FB5HxwcjClTpmi1LyIyDoW9ocmJqQVicp07d8b69esRFBSkco+Vu7s7Fi1apDasRUHXh3zsr4cPH2LDhg3IzMzEy5cvcfToURw6dAjly5fHrl27cOHCBVy4cAG1atXC1q1bUbZsWbi5ueHgwYN49OiRyr7Dw8Px6tUrAMCBAwdQsmRJWFhYIDAwEP7+/rC0tERgYCC++uorPHr0CKGhoXB0dMTbt28ZiOXy/Zg1axamTp2a6/b6YLLB2NatWzFp0iQMGjQIw4cPV1lnbW2N9PT0HLdLS0vTeiTkZ8+eqcxjldcfMxEZP2NoaN5lqoEYIKuP+Ph4fPvtt5g0aRIAYN++fTmOwK+P+qhduzb+/vtv/PTTTwBkV1/kc1MuWrQIAwYMQHh4uKI+Bg8ejFGjRiEkJASCIMDf3x/79+9H2bJlASjro0OHDvDw8MBvv/2Gf/75B1KpFI8ePcJHH32EsLAwTJo0Cd988w2sra3Rrl07NG/eHIMHD841z0U5EANkMyKMGzcOAHDmzBkcOnQICxcuzHWfBcHoJgqXP00ZGhqK/v3755jm0KFDaN++Pfz9/bFjxw4UK6Yac86YMQOTJk3Cy5cvVS5Vpqenw9raGmPGjMH8+fPznSdOFE5U9BhLQ5OdqQdi8vro0qULypcvr1j+7kThrA+loh6IZSevj/r166Ndu3acKPxDnD9/Hp06dUK9evWwbds2tUAMAHx9fQHIArvsLl26BKlUqlhPRJQTY25o2PCzPuRYH0rZ6+PTTz/NM21BMKlg7M6dOwgICICnpyf27NmT6+XG5s2bw9HRUW2+ruXLl8PGxgYBAQH6yC4RGSFjb2iKQsNva2sLQRAgCIJKrxjrQ4mBmJIm9VFQjOaesaVLlyI+Pl7xpOPu3bsVoy2PGjUKZmZmiseDv/32W7Wxwry9vfHJJ58AkN0zNn36dAQFBaFr165o3bo1Tp06hQ0bNmDGjBlwdHTUb+GIyCiYekPDhl+G9aHE+tAPo7lnzNPTE0+ePMlxnfxpEy8vr1y379evH9auXauybNWqVZg/fz4ePXoEDw8PjBw5El9//XWef/A54T1jREVDSEiIyTY0bPhlWB9KRbU+DNGmG00wVpgxGCMqGhYsWGD0DU1O2PDLsD6UinJ9GKJNN6l7xoiICpIpNDTvYsMvw/pQYn3oH4MxIqJ8MrWGhg2/DOtDifVhGAzGiIg+kDE2NGz4ZVgfSqwPw2EwpkN+fn7w8fFBSEiIobNCRHpijA0NG34Z1ocS60P2gI6Pjw9q1qyZ574KAm/g1wHewE9UNBlTQyPHhl+G9aHE+lAKDw/Htm3bMHv2bN7AT0RU2BlrQ8OGn/WRHetDSV4f79tnQTCaQV+JiAoLY25o2PCzPuRYH0rZ66Nx48Z5pi0I7BkjItKAsTc0bPhZHwDrIztN6qOgMBgjIsonU29o2PDLsD6UWB/6wWCMiCifTLmhYcMvw/pQYn3oD4MxIqJ8MtWGhg2/DOtDifWhXwzGiIjyyRQbGjb8MqwPJdaH/jEYIyLKJ1NraNjwy7A+lFgfhsFgjIjoAxljQ8OGX4b1ocT6MBwGY0REH8AYGxo2/DKsDyXWh9L9+/fz3FdB4HRIOiCfDqly5cowNzdHUFAQgoKCDJ0tItIxUf2PgIxM5QIBKJtuDtssEZ6Ls5BsnsfpVACcM8xRIlOEF5ZSJBaT5nksxwwzlM4wQ6yFFG8s8k5rn2kGl3QzxBcTEGORBeTxg982S4SyEnMkmwt4bpl3WqssEdzTzSERCYgUZ0HII62FVAQPiTmkIgHPxFJkiXL/LMwFETwkZjATRHgmzkKGWe5pRQLgLjGHWBAh0jILae/5jFkfMqwPpXzXx6tEiGITUEwKZGRk6HU6JI7Ar0MXL17k3JREpiwjExjdCQBgLgBdYq1QKbUYtjqlItk6K9fNRALQLk6MekkW2OWYhkS7zFzTAkCTBEs0TxDjqIMEEQ7peaatnVQMHd5Y4ZJdBvaVlOTZmFdKNUf3V9a4Z52JsNJpeaZ1k5ihT4wNYiyk2OCcAiGP6yilMkToF2MDiUjAujKpyMqj0bXLEqHfS2uIBRHWOacgwyL3tJZSoHeMDZwzzLDeOQVp4twbXdaHEutDSZv6uC1KxvbgeXnuV9cYjBERaejdhuaeBg1NuAYNzUkNG5q8ekrebfizNGj40zVo+JM0aPhfa9DwR2nQ8LM+WB+A9vWx0zYtz/0WBN4zRkSkAWNvaNjwsz4A1kd279aH1AD39TMYIyLKLxNoaNjwsz5YH0qa1EdBYjBGRJRPZdPNTbahYcOvxPqQYX3oD4MxIqJ8ss0SmWRDw4ZfifUhw/rQL97AT0SUT8/FWRo9FWYMDQ0bfiXWhwzrQ//YM0ZElE95jZNkjA0NG34l1ocM68MwGIwREX0gY2xo2PArsT5kWB+Gw2CMiOgDGGNDw4ZfifUhw/pQcpXoPzRiMEZEpCVjbGjY8CuxPmRYH0puEjP0eGWT57ELAm/g1yE/Pz/OTUlURBhrQ8OGX4b1IcP6UHI4cBVpR6/hd2nec10WBAZjOsS5KYmKBmNsaNjwK7E+ZFgfSm4SM/TxaYSYWg2xvngcMGppnvnQNV6mJCLSgNE2NGz4AbA+5FgfSu/WR4YBIiMGY0RE+WUCDQ0bftYH60NJk/ooSAzGiIjyyTnD3GQbGjb8SqwPGdaH/jAYIyLKpxKZIpNsaNjwK7E+ZFgf+sUb+ImI8umFpRSJJtbQsOFXYn3IsD70rxBk4f2SkpIQHByMNm3awNHRESKRCGvXrs0x7Z07d9CmTRvY2dnB0dERffr0watXr9TSSaVS/PLLL/Dy8oKVlRVq1qyJzZs3F3BJiMiYJRbL+5F3Y2to2PArsT5kWB+GUUiykbfY2FhMmzYNd+7cQa1atXJNFxkZiSZNmuD+/fuYOXMmxo8fj71796JVq1ZIT1et+IkTJ+L7779Hq1atsGTJEpQrVw49e/bEli1bCro4RGSCjK2hYcOvxPqQYX0YjlFcpnR1dUV0dDRcXFxw6dIl+Pn55Zhu5syZSE5OxuXLl1GuXDkAQP369dGqVSusXbsWQ4cOBQBERUVh/vz5CAoKwtKlsrFEBg8ejKZNm+Lbb79F165dYW5urp/CEZHRM7aGhg2/EutDhvWhVDJDhLg8c6l7hSguzJ1YLIaLi8t70/31119o3769IhADgJYtW6Jy5crYtm2bYtnff/+NjIwMjBgxQrFMJBJh+PDhiIyMxLlz53RbACIyWcbW0LDhV2J9yLA+lEpliNDTANMhGUUwlh9RUVGIiYlBvXr11NbVr18f4eHhivfh4eGwtbVF1apV1dLJ12sjMTFR5SWRSLTaDxEZB2NsaNjwy7A+ZFgfADIygVQJSiamo8dTcySlp+WZ14JgFJcp8yM6OhqA7JLmu1xdXfHmzRtIJBKIxWJER0ejTJkyEIlEaukA4Pnz51rlwcPDQ+V9cHAwpkyZotW+iKhwM5qG5v/Y8CuxPmRYH/+39wKw6xziAITkmdOCYzLBWGpqKgDZJc13WVlZKdKIxWLFv3ml08azZ89U5qbM6RhEZPyMqqEBG/7sWB8yrA8lR/+P8WWNJkgXCdjknIrkdAkwfmWe+dY1kwnGrK2tASDHS4NpaWkqaaytrfOVTlP29vacKJzIxBlbQ8OGX4n1IcP6UCqVIUK/OHtIrAVsLpOKZHNLIDX3z7mgmMw9Y/JLjPLLldlFR0fD0dFR0VPl6uqKFy9eQBAEtXQAULZs2QLOLREZI8cMM+NraNjwA2B9yLE+lDSpj4JmMsGYm5sbnJyccOnSJbV1Fy5cgK+vr+K9r68vUlJScOfOHZV058+fV6wnInpX6Qwzk2xo2PArsT5kWB/6JRLe7R4q5OTjjIWGhqJ///4q64YPH45169bhv//+U9xMf+TIEbRs2RLLly/HsGHDAMgGh61QoQKGDh2qGGdMEAQ0bdoUDx8+xJMnTzQaZywxMREODg5ISEjgZUoiEyaq7Q2M7mTobBBRQUqVAEFL9dqmG809Y0uXLkV8fLziScfdu3cjMjISADBq1Cg4ODjgxx9/xJ9//onPPvsMX3/9NZKSkjB37lzUqFEDAwYMUOzL3d0dY8aMwdy5c5GRkQE/Pz/s3LkTp06dwsaNGzngKxEREemN0fSMeXp64smTJzmue/ToETw9PQEAt27dwrhx43D69GlYWloiICAA8+fPR5kyZVS2kUqlmDNnDn777TdER0ejUqVKmDBhAnr16qVx3tgzRlQ0sGeMqAgwQM+Y0QRjhRmDMaKigcEYURFggGDMZG7gJyIiIjJGDMaIiIiIDIjBGBEREZEBGc3TlMbAz88P5ubmCAoKQlBQkKGzQ0RERPl1JBw4ehWQ6v9WegZjOnTx4kXewE9ERGSMWtSWvf5/A78+8TIlERERkQExGCMiIiIyIL0EYw8fPlSbB5KIiIiIdByM/frrr+jRo4fKsgEDBqBSpUqoXr066tWrh5iYGF0ekoiIiMio6TQYW716tcq0QwcOHMC6deswdOhQLFmyBA8fPsTUqVN1eUgiIiIio6bTpymfPHmCqlWrKt5v27YNXl5eWL58OQDgxYsXWL9+vS4PSURERGTUdNoz9u40lwcPHkTbtm0V7z09PfHixQtdHpKIiIjIqOk0GKtcuTJ27NgBQHaJ8vnz5yrBWGRkJEqUKKHLQxIREREZNZ1ephw/fjx69uyJkiVLIjk5GVWrVkXr1q0V648ePQpfX19dHpKIiIjIqOk0GOvRowdKlSqFffv2oUSJEhgxYgSKFZMd4s2bN3B0dETv3r11echChdMhERERGSkDTockEt690Ys0lpiYCAcHByQkJHA6JCITJqrtDYzuZOhsEFFB+v90SPps03V6z5i5uTk2bdqU6/qtW7fC3Nxcl4ckIiIiMmoF+jTlu7KysiASiXR5SCIiIiKjpvPpkHILthITE3HgwAGULl1a14ckIiIiMlofHIxNnToV5ubmMDc3h0gkQu/evRXvs79KliyJ9evXq02XRERERFSUffDTlPXr18eIESMgCAKWLVuGVq1aoXLlyippRCIRbG1tUbduXXTu3PlDD0lERERkMj44GGvbtq1iYNfk5GQMGzYMH3/88QdnjIiIiKgo0Ok4Y6GhobrcHREREZHJ02kwBsiemDxw4AAePnyIuLg4tScsRSIRJk+erOvDEhERERklnQZjly5dwhdffIHIyMhch7lgMEZERESkpNOhLUaMGIHU1FTs3LkTb968gVQqVXtlZWXp8pBERERERk2nPWPXr1/HjBkzEBgYqMvdGg3OTUlERGSkDDg3pU6DMXd39/eOwm/KLl68yLkpiYiIjFGL2rLX/+em1CedXqb8/vvvsWrVKiQmJupyt0REREQmS6c9Y2/fvoWdnR0qVqyIHj16wMPDQ21icJFIhLFjx+rysERERERGSyTo8Lqimdn7O9pEIpHJ3cSfmJgIBwcHJCQk8DIlkQkT1fYGRncydDaIqCD9/zKlPtt0nfaMPXr0SJe7IyIiIjJ5Og3Gypcvr8vdEREREZk8nd7ALxcVFYXNmzdj8eLFiIyMBCAbmf/NmzcFfony3r176NGjB9zd3WFjY4MqVapg2rRpSElJUUl39uxZNGrUCDY2NnBxccHo0aORlJRUoHkjIiIiepdOe8YEQcA333yDpUuXIjMzEyKRCDVq1IC7uzuSkpLg6emJadOmYcyYMbo8rMKzZ89Qv359ODg4YOTIkXB0dMS5c+cQHByMy5cv4++//wYAXL16FS1atEDVqlWxYMECREZGYt68ebh37x72799fIHkjIiIiyolOg7G5c+di8eLF+P7779GiRQu0atVKsc7BwQGdO3fGX3/9VWDB2Pr16xEfH4/Tp0+jWrVqAIChQ4dCKpXijz/+QFxcHEqWLIkff/wRJUuWxPHjxxU353l6emLIkCE4ePAg/P39CyR/RERERO/S6WXKVatWoW/fvpg5cyZ8fX3V1tesWRMRERG6PKQK+fhmZcqUUVnu6uoKMzMzWFpaIjExEYcOHULv3r1VnpLo27cv7OzssG3btgLLHxEREdG7dBqMPXv2DA0bNsx1va2tbYEOCNusWTMAwKBBg3D16lU8e/YMW7duxfLlyzF69GjY2trixo0byMzMRL169VS2tbS0hK+vL8LDw7U+fmJiospLIpF8SHGIiIiooGVkyoazULzS9Z4FnV6mdHZ2xrNnz3Jdf/nyZZQrV06Xh1TRpk0bTJ8+HTNnzsSuXbsUyydOnIiff/4ZABAdHQ1A1lv2LldXV5w6dUrr43t4eKi8Dw4ORvPmzXHs2DF89tlnaNKkSZ7bh4eHY9euXahbty4CAgIgEolyTRsREYFt27ahUqVK6NKli9rgutlFRUVh/fr1cHZ2Rq9evSAWi3NN+/r1a6xbtw5isRj9+vWDnZ1drmmTkpKwbt06SCQS9OvXD6VKlco1rUQiwcaNGxETE4M+ffrAzc0t17RZWVkICwvDvXv30K1bN1SuXDnXtIIgYO/evbh8+TI6dOiA2rVr55oWAE6ePMn6AOsjO03qg4hM0N4LwK5zBs2CToOxzp07Y8WKFejfvz8cHBwAQHHCPHjwINauXYvvvvtOl4dU4+npiSZNmuCLL75AqVKlsHfvXsycORMuLi4YOXIkUlNTASDHE66VlZVivTaePXumcunzwoULRtXQsOFXYn3IsD6UXr9+nWceichIBdQHWtdVvk9NB8av1GsWdBqMTZ06FceOHYOvry8aN24MkUiEOXPmYPLkyTh37hxq166NH3/8UZeHVLFlyxYMHToUERERcHd3ByALEKVSKb7//nt8+eWXsLa2BoAcLyGmpaUp1mvD3t5eEYydPHkSZ86cMaqGhg2/DOtDhvWhJK8PIjJBFsVkLwPS6T1jDg4O+Pfff/Hdd98hKioKVlZWOHHiBOLj4xEcHIxTp07BxsZGl4dUsWzZMtSuXVsRiMl16NABKSkpCA8PV1yelF+uzC46Ohply5b94HwYa0PDhp/1Icf6UMpeH0REBUHng75aW1tj0qRJuHr1KpKTk5GamoqbN2/ip59++qBep/x4+fJljoPKZmRkAAAyMzNRvXp1FCtWDJcuXVJJk56ejqtXr+b4FKgmjLmhYcPP+gBYH9m9Wx9ERAVBp8HYsmXL8OrVK13uUiOVK1dGeHi42vAZmzdvhpmZGWrWrAkHBwe0bNkSGzZswNu3bxVp1q9fj6SkJHTt2lXr4585c8aoGxo2/KwP1oeSJvVBRPQhdBqMjRw5Em5ubmjVqhXWrFmDN2/e6HL37/Xtt98iKysLjRs3xvTp07Fs2TK0a9cOO3fuxMCBAxWXIGfMmIE3b96gadOmWLFiBSZNmoSRI0fC398fbdq00fr4J0+eNMmGhg2/EutDhvVBRKQ7Og3G7t69i0mTJiE6OhpDhgyBq6sr2rVrh/Xr1xfo+GJyTZo0wdmzZ1G3bl0sW7YMY8aMwYMHDzBjxgwsX75cka5OnTo4fPgwrK2tMXbsWKxcuRKDBg1CWFjYBx/f1BoaNvxKrA8Z1gcRkW6JBEEQCmLHt27dwpYtW/Dnn38iIiICYrEYrVu3Ro8ePdCjR4+COKTBJCYmwsHBAQkJCSpDW7zL2BoaNvxKrA+Zol4fotrewOhOeZaDiIxcqgQIWvreNl2XCiwYy+7atWvYsmULli1bhuTkZGRmZhb0IfUqP8GYMTQ02bHhV2J9yLA+GIwRFQkGCMZ0/jTlu65fv45t27YhLCwMb9++LZKPhxtLQyPHhl+J9SHD+iAiKjgFEozdvn0bwcHBqFq1KmrXro358+fDx8cHGzZswMuXLwvikIWWsTU0bPiVWB8yrA8iooKl08uU06dPx7Zt23D79m2Ym5ujRYsW6NGjBzp27KiYHskUyS9TwrIYIBIBpR0AJ3vYZ5rBJd0M8cUExFhkAbm3M7DNEqGsxBzJ5gKeW+ad1ipLBPd0c0hEAiLFWRDySGshFcFDYg6pSMAzsRRZotyr21wQwUNiBjNBhGfiLGSY5Z5WJADuEnOIBREiLbOQZp7Hn5EAlE03h22WCM/FWUh+T1rnDHOUyBThhaUUicWkuacF4JhhhtIZZoi1kOKNRd5pWR//x/pQ0LQ+sizNgeGBeZaLiIzUkXDg6FVAKgAv44z3njELCws0bdoU3bt3R+fOnfP8FW9KFMFYyEjAWvYLvHZSMXR4Y4VLdhnYV1KSZ4NQKdUc3V9Z4551JsJKpyErj7RuEjP0ibFBjIUUG5xTkJ5H32apDBH6xdhAIhKwrkwqkvJodO2yROj30hpiQYR1zil4bZF7Wksp0DvGBs4ZZljvnIIoce6NrrkAdIm1QqXUYtjqlIp71uqD8sqJBKBdnBj1kiywyzEN4XZ531vYJMESzRPEOOogwUmH9DzTsj5kWB9K2tRHSNmUPMtFRCbAAPeM6XQypqioKDg7O0MikeDKlSuIiYnBp59+itKlS+vyMIWesTY0bPhZH3KsD6Xs9UFEVBB0es+Ys7Mzfv31V7i6uqJRo0bo3Lkzrl+/DgCIjY1F6dKl8fvvv+vykIWOMTc0bPhZHwDrI7t364OIqCDoNBgLDQ3FmDFj0KZNG6xZswbZr4CWLl0azZs3x5YtW3R5yEKlppE3NGz4WR+sDyVN6oOI6EPoNBibP38+Pv/8c2zatAmBgeo3udatWxe3bt3S5SELlXZxptnQsOFXYn0osT7I5EkygIHzZS9JhqFzQyZMp8HY/fv30bZt21zXOzo64vXr17o8ZKFyxQQbGjb8SqwPJdYHEZHu6DQYK1GiBGJjY3Ndf/v2bbi4uOjykIXKwRKm1dCw4VdifSixPoiIdEunwVi7du2wcuVKxMfHq627desWVq1ahQ4dOujykIWLCTU0bPiVWB9KrA8iIt3TaTD2888/IysrC9WrV8ekSZMgEomwbt069O7dG/Xq1YOzszN++uknXR7SKBhbQ8OGX4n1ocT6ICIqGDoNxsqWLYvLly+jTZs22Lp1KwRBwPr167F79258+eWX+Pfff4vcmGPG1tCw4VdifSixPoiICo5OB30FZGONrV69GqtXr8arV68glUrh5OQEM7MCn5O80DG2hoYNvxLrQ4n1IWMpRZ7HJSLSlk6nQyqqFNMhlSkJmImA5r5Ai7wnUyYiokJOkgEM/1X2/+WjAbGFYfNDBcuAc1PqvGesSPupl2JuSiIiIjIiLWrLXv+fm1Kf2OlOREREZEAMxoiIiIgMiMEYERERkQExGCMiIiIyIAZjRERERAbEYIyIiCgn0mxj+0VEqr4n0iEObUFERPSuy/eAjUeV7xduB0raAT2bA3UrGS5fZJLYM0ZERJTd5XtAyC4gPkl1eVySbPnle4bJF5ksBmNERERyUimw6WjeaTYf4yVL0ikGY0RERHIRUbIesLy8eStLR6QjvGdMl6Zt5NyURETG7N1Lkx+ajoxH9rkp9YzBmC5xbkoiIuNWwk636ch4cG5KIiKiQqCym+ypybw4FpelI9IRBmNERERyZmay4Svy8uVnsnREOsK/JiIiouzqVgKCOqhfinQsLlvOccZIx3jPGBER0bvqVgJ8yinvHRrbGahWnj1iVCBM8q/qypUr6NChAxwdHWFjY4Pq1avj119/VUlz9uxZNGrUCDY2NnBxccHo0aORlMSnY4iI6P+yB16V3RmIUYExuZ6xgwcPIjAwELVr18bkyZNhZ2eHBw8eIDIyUpHm6tWraNGiBapWrYoFCxYgMjIS8+bNw71797B//34D5p6IiIiKGpMKxhITE9G3b18EBAQgLCwMZrn8ivnxxx9RsmRJHD9+HPb29gAAT09PDBkyBAcPHoS/v78+s01ERERFmEn1uW7atAkvX77EjBkzYGZmhuTkZEjfmbIiMTERhw4dQu/evRWBGAD07dsXdnZ22LZtm76zTUREREWYSQVjhw8fhr29PaKiovDRRx/Bzs4O9vb2GD58ONLS0gAAN27cQGZmJurVq6eyraWlJXx9fREeHq59BlLTZYPFyV8ZmR9SHCIiIipoGZmqbXdqut6zYFKXKe/du4fMzEx8/vnnGDRoEGbNmoXjx49jyZIliI+Px+bNmxEdHQ0AcHV1Vdve1dUVp06d0j4D41eqvu/wCdCxofb7o8JJkgEM//8DIctHA2ILw+aH9Gf2eSBT/ydqMhBplvL/c84AZuaGywsVnLhIIO65QbNgUsFYUlISUlJSMGzYMMXTk507d0Z6ejp+++03TJs2DampqQAAsVh92iIrKyvFeq3MGwpYWyrfF+MXl8ikZKYDHScZOhekLxlpwJJusv8Hfg9YWBk2P1QwMjOArAzl+/QUYNVAvWbBpC5TWltbAwC+/PJLleU9e/YEAJw7d06RRiKRqG2flpamWK9dBixlc1PKXxYmFesSERGZnmIWgNhG+bK00XsWTCoYK1u2LACgTJkyKsudnZ0BAHFxcYrLk/LLldlFR0cr9kFERESkDyYVjNWtWxcAEBUVpbL8+XPZtWAnJydUr14dxYoVw6VLl1TSpKen4+rVq/D19dVLXomIiIgAEwvGunWTXdtfs2aNyvLVq1ejWLFiaNasGRwcHNCyZUts2LABb9++VaRZv349kpKS0LVrV73mmYiIiIo2k7qpqXbt2hg4cCB+//13ZGZmomnTpjh+/Dj+/PNPTJgwQXEJcsaMGWjYsCGaNm2KoUOHIjIyEvPnz4e/vz/atGlj4FIQERFRUWJSwRgArFixAuXKlUNoaCh27NiB8uXLY+HChRgzZowiTZ06dXD48GF8//33GDt2LIoXL64YCoOIiIhIn0wuGLOwsEBwcDCCg4PzTNeoUSOcOXNGT7kiIiKjY2EFjNtl6FxQEWBS94wRERERGRsGY0REREQGxGCMiIiIyIBM7p4xg5q2ETATAc19gRa1DZ0bKihSqfL/EZFAtfKAGX/XEBEZtat7gav7AEH6/rQ6xmBMl37qJZsGiUzX5XvAxqPK9wu3AyXtgJ7NgbqVDJcvIiL6ML4BspckBQjpoddD8+c8UX5dvgeE7ALik1SXxyXJll++Z5h8ERGRUWMwRpQfUimw6WjeaTYfU72ESURElA8MxojyIyJK1gOWlzdvZemIiIg0wGCMKD/evTT5oemIiIj+j8EYUX6UsNNtOiIiov9jMEaUH5XdZE9N5sWxuCwdERGRBhiMEeWHmZls+Iq8fPkZxxsjIiKNseUgyq+6lYCgDuqXIh2Ly5ZznDEiItICB30l0kTdSoBPOSBoqez92M4cgZ+IiD4IWxAiTWUPvCq7MxAjIqIPwp4xXfp6FSASAQ5lAPsyhs4NFRRplvL/c84AZuaGywvpVzFLQ+eAiAqKAeemFAmCIOj9qCYmMTERDg4OQNAWQGxj6OxQQctIA5Z0k/1/1DbAwsqw+SEiIt35/9yUCQkJsLe318sheX2FiIiIyIAYjBEREREZEIMxIiIiIgNiMEZERERkQAzGiIiIiAyIwRgRERGRATEYIyIiIjIgBmNEREREBsQR+Ik0ZWEFjNtl6FwQEZGJYM8YERERkQGxZ0yXNn0DiMwA33aAb4Chc0NERET5xbkpjRvnpiQiIjIRnJuSiIiIqGhhMEZERERkQAzGiIiIiAyIwRgRERGRAZl0MDZjxgyIRCJUr15dbd3Zs2fRqFEj2NjYwMXFBaNHj0ZSUpIBcklERERFmckObREZGYmZM2fC1tZWbd3Vq1fRokULVK1aFQsWLEBkZCTmzZuHe/fuYf/+/QbILRERERVVJhuMjR8/Hg0aNEBWVhZiY2NV1v34448oWbIkjh8/rnhs1dPTE0OGDMHBgwfh7+9viCwTERFREWSSlylPnjyJsLAwLFq0SG1dYmIiDh06hN69e6uMH9K3b1/Y2dlh27ZteswpERERFXUmF4xlZWVh1KhRGDx4MGrUqKG2/saNG8jMzES9evVUlltaWsLX1xfh4eH6yioRERGR6V2mXLFiBZ48eYLDhw/nuD46OhoA4OrqqrbO1dUVp06d0v7g6Smq780tgGIW2u+PiIiIClZmBpCVoXz/bluuByYVjL1+/Ro//fQTJk+eDCcnpxzTpKamAgDEYrHaOisrK8V6rawaqPq+QQ+gYU/t90dEREQF68KfwL9bDJoFkwrGJk2aBEdHR4waNSrXNNbW1gAAiUSiti4tLU2xXitDfgcss81Nac5eMSIiokKtflegbkfl+/QU9c6VAmYywdi9e/ewcuVKLFq0CM+fP1csT0tLQ0ZGBh4/fgx7e3vF5Un55crsoqOjUbZsWe0zYWnDicKJiIiMSTHD31JkMjfwR0VFQSqVYvTo0fDy8lK8zp8/j4iICHh5eWHatGmoXr06ihUrhkuXLqlsn56ejqtXr8LX19cwBSAiIqIiyWR6xqpXr44dO3aoLZ80aRLevn2LxYsXw9vbGw4ODmjZsiU2bNiAyZMno3jx4gCA9evXIykpCV27dtV31omIiKgIEwmCIBg6EwWpWbNmiI2Nxc2bNxXLrly5goYNG8LHxwdDhw5FZGQk5s+fjyZNmuDAgQMaHyMxMREODg5A0BZepiQiIjJmkhQgpAcSEhJUxiMtSCZzmVITderUweHDh2FtbY2xY8di5cqVGDRoEMLCwgydNSIiIipiTOYyZW6OHz+e4/JGjRrhzJkz+s0MERER0TuKZM8YERERUWHBYIyIiIjIgBiMERERERmQyd8zplebvgFEZoBvO8A3wNC5ISIiovy6uhe4ug8QpHo/tMkPbaEPHNqCiIjIRHBoCyIiIqKihcEYERERkQExGCMiIiIyIAZjRERERAbEYIyIiIjIgBiMERERERkQgzEiIiIiA2IwRkRERGRADMaIiIiIDIjBGBEREZEBcW5KXeLclERERMaJc1MaN85NSUREZCI4NyURERFR0cJgjIiIiMiAGIwRERERGRCDMSIiIiIDYjBGREREZEAMxoiIiIgMiMEYERERkQExGCMiIiIyIAZjRERERAbEYIyIiIjIgDg3pS5xbkoiIiLjxLkpjRvnpiQiIjIRnJuSiIiIqGhhMEZERERkQAzGiIiIiAyIwRgRERGRATEYIyIiIjIgkwrGLl68iJEjR6JatWqwtbVFuXLl0K1bN0RERKilvXPnDtq0aQM7Ozs4OjqiT58+ePXqlQFyTUREREWZSY0zNmfOHJw5cwZdu3ZFzZo18eLFCyxduhR16tTBv//+i+rVqwMAIiMj0aRJEzg4OGDmzJlISkrCvHnzcOPGDVy4cAGWlpYGLgkREREVFSYVjI0bNw6bNm1SCaa6d++OGjVqYPbs2diwYQMAYObMmUhOTsbly5dRrlw5AED9+vXRqlUrrF27FkOHDjVI/omIiKjoManLlA0bNlTr1apUqRKqVauGO3fuKJb99ddfaN++vSIQA4CWLVuicuXK2LZtm97yS0RERGRSwVhOBEHAy5cvUbp0aQBAVFQUYmJiUK9ePbW09evXR3h4uPYHS0+Rjdwrf2VmaL8vIiIiKniZGaptd3qK3rNgUpcpc7Jx40ZERUVh2rRpAIDo6GgAgKurq1paV1dXvHnzBhKJBGKxWPODrRqo+r5BD6BhT833Q0RERPpx4U/g3y0GzYJJB2N3795FUFAQPvnkE/Tr1w8AkJqaCgA5BltWVlaKNFoFY0N+ByyzzU1pbqH5PoiIiEh/6ncF6nZUvk9PUe9cKWAmG4y9ePECAQEBcHBwQFhYGMzNzQEA1tbWAACJRKK2TVpamkoajVnacKJwIiIiY1LMQvYyZBYMevQCkpCQgLZt2yI+Ph6nTp1C2bJlFevklyfllyuzi46OhqOjo3a9YkRERERaMLlgLC0tDYGBgYiIiMDhw4fh4+Ojst7NzQ1OTk64dOmS2rYXLlyAr6+vnnJKREREZGJPU2ZlZaF79+44d+4c/vzzT3zyySc5pvviiy+wZ88ePHv2TLHsyJEjiIiIQNeuXfWVXSIiIiKIBEEQDJ0JXRkzZgwWL16MwMBAdOvWTW197969AQDPnj1D7dq1UaJECXz99ddISkrC3Llz4e7ujosXL2p8mTIxMREODg5A0BbeM0ZERGTMJClASA8kJCTA3t5eL4c0qWCsWbNmOHHiRK7rsxf11q1bGDduHE6fPg1LS0sEBARg/vz5KFOmjMbHZTBGRERkIgwQjJnUPWPHjx/Pd9pq1arhwIEDBZcZIiIionwwqXvGiIiIiIwNgzEiIiIiA2IwRkRERGRAJnXPmMFt+gYQmQG+7QDfAEPnhoiIiPLr6l7g6j5AkOr90Cb1NKWh8GlKIiIiE2GApyl5mZKIiIjIgBiMERERERkQgzEiIiIiA2IwRkRERGRADMaIiIiIDIjBGBEREZEBMRgjIiIiMiAGY0REREQGxGCMiIiIyIA4HZIucTokIiIi48TpkIwbp0MiIiIyEZwOiYiIiKhoYTBGREREZEAMxoiIiIgMiMEYERERkQExGCMiIiIyIAZjRERERAbEYIyIiIjIgBiMERERERkQgzEiIiIiA2IwRkRERGRAnJtSlzg3JRERkXHi3JTGjXNTEhERmQjOTUlERERUtDAYIyIiIjIgBmNEREREBsRgjIiIiMiAGIwRERERGVCRDcYkEgm+//57lC1bFtbW1vj4449x6NAhrfcFAMjM0GEOjUBmBnB2U9ErN8CyF8WyF9VyAyw7y27onOjX/8uraNv1oMgGY/3798eCBQvQq1cvLF68GObm5mjXrh1Onz6t8b4UFSYtYn+wWRnAv1tk/xY1LHvRK3tRLTfAsrPshs6JfkkZjOnFhQsXsGXLFsyaNQtz587F0KFDcfToUZQvXx7fffedfjJxdW/Bb6OPY2ijsJajqJa9MJZbm21Mpc612cZUyl5YPytNFdZysOyFVpEMxsLCwmBubo6hQ4cqlllZWWHQoEE4d+4cnj17VvCZuLqv4LfRxzG0UVjLUVTLXhjLrc02plLn2mxjKmUvrJ+VpgprOVj2QqtITocUHh6OypUrq42sW79+fQDA1atX4eHhke/9KSYxSHqjusLMAihmkctGUtkov5rQdJuCPkZ6iuq/BZUvfXxWmm5jSmXXNL02ZS+M5dB0G33VuTbbmErZC+NnVVT/3oGiU/bMDNXbjP7flutzgqIiOR1S9erVUaZMGRw5ckRl+e3bt1GtWjWsWLECX331Vb739/DhQ3h7e+s6m0RERGQgDx48QIUKFfRyrCLZM5aamgqxWKy23MrKSrFeE56ennjw4AEsLCwgEokUy8VicY7HISIiosJBIpGo3KwvCAIyMjLg6emptzwUyWDM2to6x6ck0tLSFOs1YWZmprfomYiIiExLkbyB39XVFdHR0WrL5cvKli2r7ywRERFREVUkgzFfX19EREQgMTFRZfn58+cV64mIiIj0oUgGY126dEFWVhZWrlypWCaRSBAaGoqPP/5YoycpiYiIiD5Ekbxn7OOPP0bXrl0xYcIExMTEoGLFili3bh0eP36MNWvWGDp7REREVIQUyaEtANnN+pMnT8aGDRsQFxeHmjVrYvr06WjdurWhs0ZERERFSJENxoiIiIgKgyJ5zxgRERFRYcFgjIiIiMiAGIwRERERGRCDMSIiIiIDYjBGREREZEAMxoiIiIgMiMEYERERkQExGCMiIiIyIAZjRERERAZUJOem1DWpVIrHjx/DwsICIpFIsVwsFkMsFhswZ0RERJQXiUQCiUSieC8IAjIyMuDp6QkzM/30WTEY04HHjx/D29vb0NkgIiIiHXnw4AEqVKigl2MxGNMBCwsLAMCFCxfg6uqqWJ5Xz5ifnx8uXryo0XE03aagj5GYmAgPDw88e/YM9vb2BZYvfXxWmm5jSmXXNL02ZS+M5dB0G33VuTbbmErZC+NnVVT/3oGiU/Z3e8aio6NRv359RduuDwzGdEB+adLV1RXu7u752sbc3Fyjk5o22+jjGABgb29foPnSVzmKatm1OQagWdkLazkKY51rs42plL2wflZA0f17B4pu2bPfdlTQeAO/gQQFBRX4Nvo4hjYKazmKatkLY7m12cZU6lybbUyl7IX1s9JUYS0Hy154iQRBEAydCWMXGRmp6MrNb8+YKUhMTISDgwMSEhK0+rVlzFj2olf2olpugGVn2YtW2Q3RprNnTAfk94UVtScnxWIxgoODi1y5AZa9KJa9qJYbYNlZ9qJVdkO06ewZ04Gi+uuBiIjI1BiiTWfPmA75+fnBx8cHISEhhs4KERERaSAkJAQ+Pj7w8/PT+7HZM6YD7BkjIiIyDewZIyIiIipiGIwRERERGRCDMSIiIiIDYjBGREREZEAMxnSIT1MSEREZJz5NaeT4NCUREZFp4NOUREREREVMsfwkGjhwoMY7FolEWLNmjcbbERERERUl+QrGjh49CpFIpNGONU1PREREVBTlKxh7/PhxAWeDiIiIqGjiPWNEREREBvRBwdi///6LWbNmYezYsbh37x4AICUlBVeuXEFSUpJOMmhMOLQFERGRcTK6oS3S09PRo0cP/P333xAEASKRCIcOHULz5s2RlpYGd3d3jB07FhMnTiyIPBc6HNqCiIjINBjN0BaTJ0/Gnj17sHz5cvz333/IHs9ZWVmha9eu+Pvvv3WWSSIiIiJTpVUwtnnzZgwfPhxDhw6Fo6Oj2vqqVavi4cOHH5w5IiIiIlOnVTAWExODGjVq5Lre3NwcKSkpWmeKiIiIqKjQKhjz8PDA3bt3c11/5swZVKxYUetMERERERUVWgVjPXv2xG+//YZz584plskHeV21ahW2bduGvn376iaHRoRPUxIRERkno3yaMjAwEEePHkXVqlVx69Yt1KhRA2/evEFkZCTatWuHv//+G+bm5gWR50KHT1MSERGZBqN5mtLS0hL//PMPQkNDUaFCBVSpUgUSiQQ1a9bE2rVrsXv3bq0CsaSkJAQHB6NNmzZwdHSESCTC2rVr87XtkSNHMHDgQFSuXBk2NjaoUKECBg8ejOjoaLW0zZo1g0gkUnu1adNG4zwTERERfYh8TYeUE5FIhN69e6N37946y0xsbCymTZuGcuXKoVatWjh+/Hi+t/3+++/x5s0bdO3aFZUqVcLDhw+xdOlS7NmzB1evXoWLi4tKend3d8yaNUtlWdmyZXVRDCIiIqJ80zoYKwiurq6Ijo6Gi4sLLl26pNF12wULFqBRo0YwM1N29rVp0wZNmzbF0qVL8fPPP6ukd3Bw0GkgSURERKSNfAVjzZs313jHIpEIR44c0WgbsVis1oOVX02aNMlxmaOjI+7cuZPjNpmZmUhLS4OdnZ1WxyQiIiL6UPkKxqRSqeJpSblnz57h4cOHcHBwQIUKFQAAjx49Qnx8PLy9veHh4aH73GooKSkJSUlJKF26tNq6iIgI2NraIj09HWXKlMGQIUPw008/wcLCQuvjJSYmqrwXi8UQi8Va74+IiIgKlkQigUQiUbx/ty3Xh3wFY+/eu3X69Gl06NABq1atQr9+/VCsmGw3mZmZCA0Nxffff5/vG+8L0qJFi5Ceno7u3burLPf29sZnn32GGjVqIDk5GWFhYfj5558RERGBrVu3an28dwPQ4OBgTJkyRev9ERERUcGaNWsWpk6datA8aDW0RYMGDfDpp59i/vz5Oa7/5ptvcPr0aZw/f17rjMnvGQsNDUX//v013v7kyZNo0aIFOnfunK8Aa+jQoVi1ahXOnTuHBg0aaHQs+WOwz549U3kMlj1jREREhVtOPWMeHh6Ff2iL69evKy5N5sTLyws3btzQOlMf6u7du+jUqROqV6+O1atX52ubb775BgBw+PBhrY9rb2+v8mIgRkREVLiJxWK19lvftArGypYti61btyIzM1NtXWZmJrZu3WqwYSKePXsGf39/ODg4YN++fShevHi+tpNfYnzz5k1BZo+IiIhIhVZDW3z33XcYNmwYGjRogGHDhinmobx37x5WrFiBq1evYtmyZTrNaH68fv0a/v7+kEgkOHLkCFxdXfO97cOHDwEATk5OBZU9IiIiIjVaBWNDhw6Fubk5Jk6ciKFDhyqetBQEAU5OTlixYgWGDBmi04xmFx0djYSEBHh7eyuefkxOTka7du0QFRWFY8eOoVKlSjlum5iYqHYvlyAIinHIWrduXWD5JiIiInqX1oO+Dho0CP369cOlS5fw5MkTAED58uVRr149xdOV2li6dCni4+Px/PlzAMDu3bsRGRkJABg1ahQcHBwwYcIErFu3Do8ePYKnpycAoFevXrhw4QIGDhyIO3fuqIwtZmdnh44dOwIArly5gi+//BJffvklKlasiNTUVOzYsQNnzpzB0KFDUadOHa3z7ufnB3NzcwQFBSEoKEjr/RAREZF+hYSEICQkBFlZWXo/tlZPUxYkT09PRXD3Lnnw1b9/f7VgLK/typcvj8ePHyv28f333+PixYt48eIFzMzMULVqVQwZMkSll08TnCiciIjINBiiTdc6GMvKysKGDRuwd+9elZ6x9u3bo1evXlpNFG6sGIwRERGZBqMJxhISEtC6dWtcvHgRxYsXVxmBPzExEfXr18eBAweKTGDCYIyIiMg0GKJN12poi4kTJ+Ly5ctYsmQJXr16hStXruDKlSuIiYnB0qVLcenSJUycOFHXeSUiIiIyOVr1jLm5uaFLly5YvHhxjutHjx6NsLAwxU34po49Y0RERKbBaHrGXr9+jY8++ijX9VWqVOHgqURERET5oFUwVrFiRezatSvX9bt27YK3t7fWmTJWfn5+8PHxQUhIiKGzQkRERBoICQmBj48P/Pz89H5srS5TLlu2DCNHjkSbNm0wZswYVK5cGQDw33//4ddff8U///yDpUuXYvjw4TrPcGHEy5RERESmwRBtulajs44YMQIxMTGYPXs2Dhw4oLLOwsICP/30U5EJxIiIiIg+xAcN+hobG4vDhw+rjDPWsmVLlC5dWmcZNAbsGSMiIjINRtMzJle6dGn06NFDV3khIiIiKnI+KBjLyMhAVFQU4uLikFMH24fM80hEZGzCw8Oxa9cu1K1bFwEBAXlOrxYREYFt27ahUqVK6NKlS56zlkRFRWH9+vVwdnZGr169IBaLc037+vVrrFu3DmKxGP369YOdnV2uaZOSkrBu3TpIJBL069cPpUqVyjWtRCLBxo0bERMTgz59+sDNzS3XtFlZWQgLC8O9e/fQrVs3xX3FOREEAXv37sXly5fRoUMH1K5dO9e0AHDy5EkcO3YMn332GZo0aZJnWtaHDOtDSZP60CetgrH4+HiMHz8eGzduRHp6utp6QRAgEokMMtmmIXGicKKiy9gaGjb8SqwPmaJeH0Y3UXjHjh2xe/du9OjRAx9//DEcHBxyTNevX78PzqAx4D1jREWbMTQ02bHhV2J9yLA+lJ4/fw43N7fCf8/YwYMHMXr0aCxcuFDX+SEiMirG1tCw4VdifciwPpSSkpKwadOmPMtUELQa9LVUqVKoWLGirvNCRGRUjLGhYcMvw/qQYX0oyesjp9uvCppWwdjQoUOxZcsWSKVSXeeHiMgoGGtDw4af9SHH+lDKXh89e/bMs2wFQavLlJMnT4ZEIkG9evXQp08fuLu75/jhde7c+YMzaAr4RVLiiU2J9SHD+lBifSixPmRYH0r6qg8LC4s8y1cQtArGoqKicPToUVy9ehVXr17NMU1RfJoyJ/wiKfHEpsT6kGF9KLE+lFgfMqwPJX3WR2JiYp5lLAhaBWMDBw7ElStXMGHChDyfpixq3h3agl8kJZ7YlFgfMqwPJdaHEutDhvWhpK/62LJli3ENbWFra4vx48dj6tSpBZGnfElKSsLcuXNx/vx5XLhwAXFxcQgNDUX//v3ztX18fDy+++477NixAykpKahfvz7mz5+v1UC1OQ1twS+SEk9sSqwPGdaHEutDifUhw/pQMkR9GGK4Kq1u4HdxcYGjo6Ou86KR2NhYTJs2DXfu3EGtWrU02lYqlSIgIACbNm3CyJEj8csvvyAmJgbNmjXDvXv3Pjhv/CIp8cSmxPqQYX0osT6UWB8yrA+lwlAf+qJVMPbNN99g9erVSEpK0nV+8s3V1RXR0dF48uQJ5s6dq9G2YWFhOHv2LNauXYvg4GAEBQXh+PHjMDc3R3Bw8Afli18kpcLwRWJ9KLE+lFgfMqwPJdaHEutD/7S6ZywtLQ0WFhaoWLEiunXrBg8PD7UPVSQSYezYsTrJZE7EYjFcXFy02jYsLAxlypRRedrTyckJ3bp1w4YNGyCRSPKs+NxIJBKEhYXxi4TC8UXiiU2J9aHE+pBhfSixPpRYH4ahVTA2fvx4xf+XLl2aY5qCDsY+RHh4OOrUqQMzM9WOwfr162PlypWIiIhAjRo1NN7vzJkzkZKSAn9/fzx48ABv375FlSpVkJaWhtu3byvSZWVl4ciRI8jMzES3bt0gCAKuXLmisi9PT084OjoiJiYGmzZtwp07d9C0aVMIgoB79+6hUqVKyMrKwrVr11S2u3LlCp4+fYqWLVvCzc1Nbb9ubm4oU6YM4uLisH//fpw4cQJVq1aFq6sr7t69i6pVqyo+o+y3Ez558gRXr16Fj48P/Pz81I5bpkwZuLm54e3btzhz5gz27t0LR0dH+Pj4qHyeN27cQEZGhmK7+Ph4/Pvvv3BwcECrVq0QERGhst9SpUqhfPnySE1NxeXLl7Fnzx6kp6ejffv2ePr0qeLLdOfOHaSmpiq2S09Px8WLF5GcnIw2bdrg5cuXePnypWK9g4MDvL29kZGRgatXr+LIkSN4+vQpWrVqhaSkJGRlZcHc3Bz37t3D27dvFdsJgoC7d+/i/v37ivrI/hnb2trio48+UtTFlStXcPHiRfj5+cHOzg5paWmwsrLCo0ePEBcXp1LW2NhYnDt3TlEf4eHhinVisRjVqlUDAFy/fh0PHjzAoUOHUK5cOVSoUAGpqamws7NDZGQkYmJiVPabmZmJo0ePwsHBAT4+Prh165ZinZmZGXx9fQEAt2/fxosXL7Bnzx5YWlqiffv2yMzMBAC8ePECz58/V9mvhYUFTp48ieTkZNSrVw9PnjzBkydPFOt9fX1hZmaGiIgIvHnzBvv378ebN28QEBCgOLnGxsbi6dOnKvu1trbG9evX8d9//6FGjRpISkpS+YyrV68OS0tLPHz4EHFxcTh16pTi++Hq6gpA9nf18OFDlf1aWVkhNjYWx44dg7u7O+zs7FT2W6VKFdjY2ODp06eIjY3F3bt3Fd+PWrVqQSQSISkpSe1vtFixYrCyssK2bdtgZWWFChUqqHw/KlasCHt7e0RHRyM6OhoxMTGK70fDhg0hFovVzhFy5cuXx7p16xAfHw9/f3+VY8vPEa9evcKzZ8+QkpKi+H706NEDpUqVyvEcAQCVK1fGtm3bcOfOHbRo0ULl+5H9HPHo0SPF+erp06cIDAxUNPzvniPkn+HRo0dx5MgR+Pn5qXw/sp8j5LeDyL8fDRo0UDT8754jAKBSpUq4f/8+NmzYAA8PD5XvR/ZzxJ07dwDIzlfy74e84X/3HAEAXl5eSElJwbJly1CsWDGV70f2c8SNGzcAyP6u5N+PCRMmwM7OTu0cAQAeHh6wtrbGsmXL8OLFC7Rv317x/Xj3HAHIzlfy78eYMWPg5uaW4znC1dUVzs7O+OOPP3Dx4kXF+erKlStq54jMzEwIgqD4fvTt2xe1a9fO8RxRunRplCtXDgcPHsRff/2lOF9duXJF7RyRlpYGAIrvx2effYaAgAC8fPlS7RxRokQJVKhQATdv3sSyZcsU5yv532T2c4T8Spv8+1GxYkX06tULb9++VTlnAYCdnR0qV66MV69eYdasWYrzlfz7kf0cER8fDwCK74eVlRVGjRoFc3Nzlc8/PT0dVlZWqFKlCmJjY6F3ghYeP36cr5e+XLx4UQAghIaG5iu9ra2tMHDgQLXle/fuFQAI//zzj0bHT0hIEACovWrUqCEIgiDcu3cvx/X//fefIAiC0KBBA7V169evF6RSqTBs2DC1df7+/nke9++//xYEQRACAwPV1s2fP18QBEGYM2eO2rratWsrymRpaam2fu7cuUJmZqYwaNAgtXU//PCDIAiCsG3bNrV1bm5uiv26ubmprR89erTw9u1b4YcfflBbN2jQIEEQBOHChQtq6ywtLRX7rV27ttr6nj17CpGRkcL8+fPV1gUGBgqCIAjR0dE5foYJCQmCIAiCv7+/2rp27doJV65cEdavX6+2rkGDBoo85bTfe/fuCYIgCL169VJb17RpU2H37t3C/v371dZ5e3sr9luyZEm19WfPnhUEQRDGjh2bY57WrFkjnDt3Tm1d8eLFFfv96KOPcv1bmjlzZo5/L/PnzxeuX7+eY1nT0tIEQRCExo0bq61btWqVIAiCsGrVKrV1VatWFaZPny7cuHEjx/0+e/ZMEARB6NKli9q6mTNnCoIgCH///bfaOk9PT2HKlCnCiRMnhOLFi6utv3z5siAIgjBixAi1dWPHjhUEQRDOnj2rtq5kyZLC9OnThS1btgje3t5q6+Xnk+DgYLV1vXr1yvMcMX/+fGHp0qWCn5+f2rr169cLgiAIS5cuVVv3vnPEokWLhFmzZgmtWrXK8Zi5fZffd45YunSpMGXKFKFjx45q6+TniGPHjqmte985YuXKlcKUKVNyrHP5OeLmzZtq6953jlixYoUwa9YsoVu3bmrr5OeImJiYHD/DvM4R8nrr2bOn2roPOUdMnjxZ2LJli9C3b1+1ddnPEaVLl1Zbn9c5YsSIEcKJEyeEoUOHqq3Lfo7w8fFRW79z585czxFdunQR/vvvP2H8+PE5llV+jmjatKnaumXLluV6jmjatKkQGxubYzsG5H2OmDRpUq7niNzqWB+0epqysLl06RL8/Pzy/TSlubk5vvrqKyxbtkxl+dGjR9GiRQvs2LEDHTt2zPfx5U9e/PHHHyhfvrxiubOzs0rPWPZfmK1atUKnTp0AAP/99x+Sk5NV9lm+fHmcO3cOJ0+ehK+vL6pUqaJYV7x4cbWesew9MPJB6x48eICEhASV/bq5ueH58+fYunUrXFxc0LhxY0XXsrW1tVrPWPZfmCNHjoSdnR2ePHmC169fq+y3TJkyAIDVq1dDKpWibdu2sLS0BCDrQXm3Zyz7L8yvv/4arq6uiIqKUum9AmS/ekuVKoVVq1bh2bNnaN++PUqUKAFA1vsq73KX/+rN/gtz2LBhqFatGl6+fImoqCiV/To4OMDT0xNbtmzBmTNn0KpVK5W6q1WrlkrPmJDtF2b37t3RvHlzvHnzBo8fP1bZr/xX78mTJ7Fx40b4+fmpPKHr4+Oj1jMm/4X56aefok+fPnj79i3u37+vsl/5r96IiAj8+uuvcHd3R4sWLRRd/ZUrV1brGZP/wixXrhxGjx6NrKws3L17V2W/8l+9r1+/xpw5c2BmZob27dvDxsYGAFChQgWUKFFCpWdM/gvT3Nwc48aNQ/HixXHz5k28y9fXFxkZGZg/fz6io6MREBAAZ2dnAEC5cuVQunRplZ4x+ffjxYsXGDFiBCpWrJjjOIbVq1eHhYUF1qxZg4sXL6Jp06aK70fZsmXh4uKi1jN25coVhIeHo3v37mjSpAmuXr2qNoOIvGds3759+Pvvv1G1alXF98PZ2Rnu7u5qPWNPnjzBkSNH0LhxY3Tp0gV3796FRCJR2a+8Z+zKlSsIDQ2Fo6Oj4vtRsmRJeHl5qfWMyb8f3t7e6NevH6KiotTOEfKesUePHuG3335T9BiXKFEix3MEoOyBsbS0RP/+/ZGWlpbjOaJMmTKIjY3F6tWrFeer8uXL53iOAKD4fsTGxuKLL76Ao6NjjucIec/Y1q1bFeerOnXq5HiOkLt79y5u3ryJhg0bwtfXV61XJ3vP2D///KM4X7Vo0QLFihVTO0fIxcTE4Ny5cyhfvjxatmyp1huSvWfs1KlTKj3GNjY2aucIuZSUFJw+fRqWlpbo0KGD2phV2XvG/v33X5UeY2dn5xzPEYDs+3H58mXExsaiXbt2eFf2nrFr167h2LFjih7jKlWq5HiOkHvw4AFu376NTz75BKVLl1ZZ927P2NWrVxU9xo0bN4a3t7faOUIuLi4O586dg5eXFypXrqx2aTJ7z9jDhw8VPcZt27ZFxYoV1c4RchkZGThz5gwsLCxQu3ZtxflKLnvP2PPnz1WuqPj4+KidI+Q9Y2ZmZrh27RoiIyMxe/ZsvT5NWSSDMTs7O3Tv3h1r1qxRWb5v3z4EBATgn3/+QevWrfN9/Pw8Bstr/Eq850KJ9SHD+lBifSixPmRYH0r6qI+OHTuiatWqeg3GtLpnzNjJn8R8l3xZ2bJldXo8fpGUeGJTYn3IsD6UWB9KrA8Z1oeSvuqjePHiuaYtKFoNbWHsfH19ceXKFbXLFOfPn4eNjU2ef+ya4hdJiSc2JdaHDOtDifWhxPqQYX0oFYb6KEgmH4xFR0fj7t27KvcgdOnSBS9fvsT27dsVy2JjY/Hnn38iMDBQq2EtcsIvklJh+CKxPpRYH0qsDxnWhxLrQ4n1oR/5vkwpf9S/MFm6dCni4+MVNw3u3r0bkZGRAIBRo0bBwcEBEyZMwLp16/Do0SN4enoCkAVjDRo0wIABA3D79m2ULl0ay5YtQ1ZWls6meOIXSakwfJFYH0qsDyXWhwzrQ4n1ocT60J9838BfokQJtG7dGu3bt0ebNm3g5ORU0Hl7L09PT5VxjbKTB1/9+/dXC8YA2VMe3377LXbu3InU1FT4+flh3rx5qFevnsb5kN/AL39aZPjw4XB2duYXCYXji8QTmxLrQ4n1IcP6UGJ9KBXF+ggJCVFMFB4REVE4n6ZcsmQJ9u3bhxMnTiA9PV1RQQEBAahbt25B57NQy/40pa2tLb9I/8cTmxLrQ4b1ocT6UGJ9yLA+lAxZH4aYKFzjoS1SU1Nx+PBh7N+/H/v27cPTp0/h4uKCNm3aoH379mjVqpVBnkQwJENUHBEVLsbS0Mix4VdifSixPowkGHvXzZs3sXfvXuzfvx9nz56FSCRCo0aN0K5dOwQEBKgMVmqqGIwRFW3G1NAAbPizY30osT5k4uLi4OjoqNc2/YOfpqxevTq+//57HD9+HK9evcL69evh4eGBuXPnolq1apgzZ44u8klEVCgZW0PDhl+J9aHE+pDJysrCzp07c11fUHQ66KuDgwO6deuGbt26AQAuXryoy90TERUqxtjQsOGXYX0osT5k5PXx4MGDXNMUlAIdgd/Pz68gd1/o+Pn5wdzcHEFBQQgKCjJ0doioABlrQ8OGn/WRHetDJisrC0OGDMGePXvyPH5BKZLTIRWUixcv8p4xoiLAGBsaNvwyrA8l1oeMvD4qVKiA06dPw8XFBQ4ODrmmLwgmMVG4ofEGfiIiItNgiDbd5KdDIiIiIirMGIwRERERGRCDMSIiIiID0voG/hcvXmDNmjW4cuUKEhISIJVKVdaLRCIcOXLkgzNIREREZMq0CsauX7+OZs2aITU1FR999BFu3LgBHx8fxMfHIyoqCt7e3vDw8NB1Xgs9Dm1BRERknLJPFK5vWj1N2a5dO9y8eROnT5+GjY0NnJ2dcfjwYTRv3hx//vknhg8fjn379qF+/foFkedCh09TEhERmQajeZryzJkz+Oqrr1CuXDmYmcl2Ib9M2bVrV/Tq1Qvffvut7nJJREREZKK0CsakUinKlCkDAChRogTMzc3x5s0bxfoaNWrg8uXLuskhERERkQnTKhjz8vLCo0ePZDswM4OXlxcOHz6sWH/27FmUKFFCJxkkIiIiMmVaBWP+/v74888/Fe+HDx+O1atXo2XLlmjRogXWrVuHnj176iyTRERERKZKqxv44+Li8PDhQ9SsWRMWFhYQBAEzZszAX3/9BXNzc7Rv3x4//vgjLC0tCyLPhY78Zr/KlSvzacoiIDk5WTF/WlJSEmxtbQ2cIyIi+lDZn6aMiIjQ6w38nJtSB/g0ZdHCYIyIyHQZzdOUzZs3z3NA12PHjqF58+ZaZUgikeD7779H2bJlYW1tjY8//hiHDh1673aenp4QiUQ5vipVqqSSNrd0s2fP1irPRERERNrSatDX48ePY/Dgwbmuj4mJwYkTJ7TKUP/+/REWFoYxY8agUqVKWLt2Ldq1a4djx46hUaNGuW63aNEiJCUlqSx78uQJJk2aBH9/f7X0rVq1Qt++fVWW1a5dW6s8ExEREWlL6+mQRCJRruvu37+P4sWLa7zPCxcuYMuWLZg7dy7Gjx8PAOjbty+qV6+O7777DmfPns11244dO6ot+/nnnwEAvXr1UltXuXJl9O7dW+M8EhEREelSvoOxdevWYd26dYr3P//8M1atWqWWLj4+HtevX0e7du00zkxYWBjMzc0xdOhQxTIrKysMGjQIP/74I549e6bRNEubNm2Cl5cXGjZsmOP61NRUiEQiWFlZaZxXIiIiIl3I9z1jKSkpePXqFV69egUAePv2reK9/BUbGwuxWIxhw4Zh9erVGmcmPDwclStXVrthTj6t0tWrVzXa1507d3IdYmPt2rWwtbWFtbU1fHx8sGnTJo3z+67ExESVl0Qi+eB9EhERUcGRSCRq7be+5btnbPjw4Rg+fDgA2aCvixcvRocOHXSamejoaLi6uqotly97/vx5vve1ceNGADlfomzYsCG6desGLy8vPH/+HCEhIejVqxcSEhIUZdTGu712wcHBmDJlitb7o8Ip+ySyJ0+ehL+/P8zNzQ2YIyIi0tasWbMwdepUg+ahUA1t4e3tjY8++gj79u1TWf7w4UN4e3tj4cKFGDNmzHv3I5VKUa5cOTg7O+PKlSvvTZ+eno66desiMjISz58/h7W1tUb5lj8G++zZM5VePbFYDLFYrNG+qHDbvn07Ro8ejaioKMUyd3d3LF68GJ07dzZgzoiISBsSiUTlSlZiYiI8PDwK/9AW2b19+xaRkZF4+vSp2ktT1tbWOV7aS0tLU6zPjxMnTiAqKirHXrGcWFpaYuTIkYiPj/+gOTXt7e1VXgzETMv27dvRpUsXlUAMAKKiotClSxds377dQDkjIiJticVitfZb37R+mnL58uVYsGABHj58mGua7Jdz8sPV1VWtoQNkly8BoGzZsvnaz8aNG2FmZoYvv/wy38eWX2LMPuE5kVxWVha+/vpr5NSRLAgCRCIRxowZg88//5yXLImISCNa9YytWLECQUFBqFixIn7++WcIgoAxY8bghx9+gIuLC2rVqoU1a9ZovF9fX19ERESo3Tx3/vx5xfr3kUgk+Ouvv9CsWbN8B28AFEGlk5NT/jNMRcapU6cQGRmZ63pBEPDs2TOcOnVKj7kiIiJToFUwtmTJErRu3Rr79+9XDEMREBCAGTNm4Pbt23j79i1ev36t8X67dOmCrKwsrFy5UrFMIpEgNDQUH3/8saL36unTp7h7926O+9i3bx/i4+NzvUQpfxo0u7dv32LRokUoXbo06tatq3G+yfTJe2d1lY6IiEhOq8uUDx48UEyEbWFhAUB2EzwAODg4YPDgwVi2bBm++eYbjfb78ccfo2vXrpgwYQJiYmJQsWJFrFu3Do8fP1bpaevbty9OnDiR4yWjjRs3QiwW44svvsjxGCEhIdi5cycCAwNRrlw5REdH4/fff8fTp0+xfv36D5rc3M/PjxOFm6icnvL9kHRERFS4ZJ8oXN+0CsYcHByQmZkJQHbTuo2NDZ49e6ZYX7x4cbx48UKrDP3xxx+YPHky1q9fj7i4ONSsWRN79uxBkyZN3rttYmIi9u7di4CAADg4OOSY5tNPP8XZs2exevVqvH79Gra2tqhfvz5+//13refTlLt48SInCjdRjRs3hru7O6KionL8ESASieDu7o7GjRsbIHdERPSh5B0p8hES9EmroS1atGgBT09PRW9Vy5Yt8ebNG+zevRtSqRTt27eHmZkZwsPDdZ7hwsgQM7yT/smfpgSgEpDJpwYLCwvj8BZEREbOEG26VveM9e7dGzdv3lQMQzF16lTcuXMH5cqVg6enJ/777z/FvJBEpqJz584ICwtTezDE3d2dgRgREWlNZ4O+Pnz4ELt374a5uTn8/f1RuXJlXezWKLBnrGjJ3oW9b98+jsBPRGRCDNGmF6oR+I0Vg7GiJTk5GXZ2dgCApKQk2NraGjhHRFQQ+F0vmozmMiURERER6Ua+gjEzMzOYm5tr/Cpq/Pz84OPjg5CQEENnhYiIiDQQEhICHx8f+Pn56f3Y+Rra4qefflI8MSa3Y8cO3Lp1C61bt8ZHH30EALh79y4OHjyI6tWro2PHjjrPbGHHoS2IiIiMkyGHtshXMDZlyhSV9ytXrkRMTAxu3rypCMTk7ty5g+bNm2s0FRERERFRUaXVPWNz587FyJEj1QIxAKhatSpGjhyJX3755YMzR0REZCjZR2I/efKkQUZmp6JBq2AsMjJSMQ1STiwsLPKcVJmIiKgw2759O3x8fBTv27VrB09PT2zfvt2AuSJTpVUwVr16dSxbtgxRUVFq6yIjI7Fs2TLUqFHjgzNHVBjZ2tpCEAQIgsBH3YlMkHy2jXfbuKioKHTp0oUBGemcVuOMnT59Gq1btwYAdOrUCRUrVgQA3Lt3Dzt37oQgCDh48CAaNWqk29wWUvKb/SpXrsyJwomIjFhWVhY8PT1zvbojn4f20aNHRXLUAFOWfaLwiIgI4xj09ebNm5g8eTIOHjyI1NRUAIC1tTVat26NqVOnFqmeMQ76SkRkGo4fP47PPvvsvemOHTuGZs2aFXyGSO8M0abn62nKnFSvXh07duyAVCrFq1evAABOTk4wM+M4skREZJyio6N1mo4oP7QOxuTMzMxQpkwZXeSFiIjIoFxdXXWajig/2I1FRET0f40bN4a7u7vaQOdyIpEIHh4eaNy4sZ5zRqaMwRgREdH/mZubY/HixQCgFpDJ3y9atIg375NOMRgjIiLKpnPnzggLC1ObScb9f+3df1RUZf4H8Pfwa2D4mYgh8RvRFBTSFZJcwfZomvkLYTiKW6y6rAb+yhQ39ICaGtumlaKky5HW5GhLeSplW81Y11hDEXU7B1hM1EgwEAICGRD4fP/wO7NeZgaHmZELM5/XOfxxn/s8c5/PPHPvfLg/nvH0RF5eHqKjo0XqGTNVnIwZEf9QOGOMmYbo6GiUlpaqlvPz83Hjxg1OxEyYmD8UPmiTsfb2dqSkpMDDwwN2dnYIDw/H6dOnH9kuPT0dEolE7c/W1tagvgAP5l8rLS01mznG2tvbkZ6erorfnHDs5he7ucYNmG/sD1+KDA8PN7tLk+Y27klJSSgtLcU333wDAP0at97zjIlt4cKFyMvLw5o1axAYGIicnBxcvHgRBQUFvU42m56eji1btmD//v1wcHBQlVtaWmLhwoV69eXHH3+El5cXqqqq4OnpqddrDEbmPL8ax25+sZtr3ID5xt7a2qr6nqiurja7JyjNddzF+E4flGfGLly4gKNHj2Lnzp14++23kZiYiK+//ho+Pj7YsGGDTq8RExODxYsXq/70TcT0pc+lzL626Y9t6GOgxmGusQ/EuPVpYypjrk8bU4l9oL5XfTVQ4+DYBzAahNavX0+WlpbU1NQkKN+xYwcBoB9++EFr27S0NAJAtbW11NTURN3d3Qb3p6qqigBQVVWVzm1Gjx7d5+30tc3j3kZTUxMBUBsHY25Dn/r90caUYu9rfX1iH4hx9LVNf425Pm1MJfaB+F6Z6+edyHxj1+c73VAGT/oqhsuXL2PkyJFqp03DwsIAAFeuXIGXl1evr+Hv74+WlhbY29tj3rx5eOedd/SevJb+/0pvzxmZpVIppFKpxjZdXV1obm7u03b62uZxb0NZb7DHoU8bU4q9r/X1iX0gxtHXNv015vq0MZXYB+J7Za6fd8B8Ym9vbxfcH6b8Lqd+vItrUN4zFhwcjCeffBJnzpwRlJeWliIoKAhZWVn4wx/+oLHte++9h++//x6TJk2CVCrFuXPnkJmZCT8/PxQXF+t1XbyyshIBAQF6xcIYY4yxgef69evw9/fvl20NyjNjbW1tGs84KZ+IVP5wuSarV68WLC9YsABhYWGIj4/Hvn37sHHjxj73x9fXF9evX4e1tbVgksDezowxxhhjTHw9z4wREe7fvw9fX99+68OgTMbs7Ow0PnKqUChU6/ti0aJFWLduHb766iu9kjELC4t+y54ZY4wxZloG5dOUw4cPV7s/C/jfdd6esybrwsvLCw0NDQb3jTHGGGOsLwZlMhYaGoqKigq1m/GKiopU6/uCiHDz5k24ubkZq4uMMcYYYzoZlMlYTEwMurq6cODAAVVZe3s7Dh06hPDwcNWTlD/88APKy8sFbevq6tReb//+/airq8OMGTMeb8cZY4wxxnoYlE9TAoBcLsfx48exdu1ajBgxAh9++CEuXLiAM2fOYMqUKQCAqKgonD17VvB4qkwmQ1xcHMaOHQtbW1t88803OHr0KEJCQlBYWAiZTCZWSIwxxhgzQ4PyBn4A+Otf/4rNmzfj8OHD+PnnnzFu3DicOHFClYhpEx8fj3//+9/45JNPoFAoVLP2p6amciLGGGOMsX43aM+MMcYYY4yZgkF5zxhjjDHGmKngZIwxxhhjTEScjDHGGGOMiYiTMcYYY4wxEXEyhgdzlKWkpMDDwwN2dnYIDw/H6dOndWp7+/ZtyOVyuLi4wMnJCXPnzkVlZaXGutnZ2Rg9ejRsbW0RGBiIPXv2GDMMvegb+6effoq4uDj4+/tDJpNh1KhRWLduHRobG9Xq+vr6QiKRqP0tX778MUSkO31jT09P1xiP8rdRezKlcdc2lhKJBIGBgYK62uq99dZbjyusR2ppaUFaWhpmzJiBIUOGQCKRICcnR+f2jY2NSExMhJubG+zt7TF16lSUlJRorPv5559j/PjxsLW1hbe3N9LS0tDZ2WmkSPrOkNjPnDmDJUuWYOTIkZDJZPD398eyZcs0/hJKVFSUxnEXax5HQ+LOycnR+jm+c+eOWn1TGnNt4yiRSGBtbS2oO9CO8RcvXkRycjKCgoJgb28Pb29vyOVyVFRU6NRejP180E5tYUwJCQnIy8vDmjVrEBgYiJycHLz44osoKCjA5MmTtbZraWnB1KlT0dTUhDfeeAPW1tbYvXs3IiMjceXKFbi6uqrqfvDBB1i+fDkWLFiA1157DefOncOqVatw7949pKSk9EeYGukbe2JiIjw8PLB48WJ4e3vju+++w969e5Gfn4+SkhK13wcNDQ3FunXrBGUjR458LDHpSt/Ylfbv3w8HBwfVsqWlpVodUxv3d999Fy0tLYKyW7duYdOmTZg+fbpa/WnTpuHll18WlD3zzDPGCUIPd+/exdatW+Ht7Y2QkBD885//1Lltd3c3Zs2ahatXr2L9+vUYOnQo9u3bh6ioKFy6dEmQjP7973/HvHnzEBUVhT179uC7777Dm2++idraWuzfv/8xRPZohsSekpKChoYGxMbGIjAwEJWVldi7dy9OnDiBK1euwN3dXVDf09MTO3fuFJTp8zN1xmBI3Epbt26Fn5+foMzFxUWwbGpjnpqaimXLlgnKWltbsXz5co37+kA6xmdkZKCwsBCxsbEYN24c7ty5g71792L8+PH49ttvERwcrLWtaPs5mbmioiICQG+//baqrK2tjQICAmjSpEm9ts3IyCAAdOHCBVVZWVkZWVpa0h//+EdV2b1798jV1ZVmzZolaB8fH0/29vbU0NBgpGj6xpDYCwoK1Mo+/PBDAkAHDx4UlPv4+KjFLjZDYk9LSyMAVFdX12s9Uxx3TbZt20YAqLCwUFAOgJKSkgzurzEpFAqqqakhIqKLFy8SADp06JBObY8dO0YA6G9/+5uqrLa2llxcXGjhwoWCumPGjKGQkBC6f/++qiw1NZUkEgmVlZUZHogeDIn97Nmz1NXVpVYGgFJTUwXlkZGRFBQUZJQ+G4MhcR86dIgA0MWLFx9Z19TGXJPDhw8TADpy5IigfKAd4wsLC6m9vV1QVlFRQVKplOLj43ttK9Z+bvaXKfPy8mBpaYnExERVma2tLZYuXYrz58+jqqqq17YTJ07ExIkTVWVPP/00fvOb3+Djjz9WlRUUFKC+vh6vvvqqoH1SUhJaW1tx8uRJI0akO0Nij4qKUiubP38+AKCsrExjm46ODrS2thrWaSMxJHYlIkJzc7PgFx4eZorjrklubi78/PwQERGhcX1bWxsUCoVBfTYWqVSqdhZHV3l5eXjyyScRHR2tKnNzc4NcLsdnn32G9vZ2AEBpaSlKS0uRmJgIK6v/XXx49dVXQUTIy8szLAg9GRL7lClTYGFhoVY2ZMgQrft7Z2en2llUMRgS98N++eUXdHV1aVxnimOuSW5uLuzt7TF37lyN6wfKMT4iIgI2NjaCssDAQAQFBWn9vCqJtZ+bfTJ2+fJljBw5Ek5OToLysLAwAMCVK1c0tuvu7sZ//vMf/OpXv1JbFxYWhuvXr+OXX35RbQOAWt0JEybAwsJCtb6/6Ru7Nsp7KIYOHaq27uuvv4ZMJoODgwN8fX3x3nvv6ddpIzFG7P7+/nB2doajoyMWL16Mn376SW0bgGmP++XLl1FWVoZFixZpXJ+TkwN7e3vY2dlhzJgxyM3N1bvfYrt8+TLGjx+vlpSEhYXh3r17qvtRtI27h4cHPD09RRt3Y2tpaUFLS4vG/b2iogL29vZwdHSEu7s7Nm/ejPv374vQS+OYOnUqnJycIJPJMGfOHFy7dk2w3hzGvK6uDqdPn8a8efNgb2+vtn6gHeN7IiL89NNPGj+vDxNrPzf7e8ZqamowfPhwtXJlWXV1tcZ2DQ0NaG9vf2TbUaNGoaamBpaWlhg2bJigno2NDVxdXbVu43HTN3ZtMjIyYGlpiZiYGEH5uHHjMHnyZIwaNQr19fXIycnBmjVrUF1djYyMDP0DMIAhsT/xxBNITk7GpEmTIJVKce7cOWRmZuLChQsoLi5WJTnmMO5HjhwB8OBnxnqKiIiAXC6Hn58fqqurkZmZifj4eDQ1NWHFihV69l48NTU1Gn9u7eH3bezYsaqb2rW9x2KNu7G9++676OjoQFxcnKA8ICAAU6dOxdixY9Ha2oq8vDy8+eabqKiowLFjx0TqrX5kMhkSEhJUydilS5ewa9cuREREoKSkBF5eXgBgFmN+7NgxdHZ2atzXB+IxvqcjR47g9u3b2Lp1a6/1xNrPzT4Za2trg1QqVStXPhnX1tamtR0Andq2tbWpnTJ9uK62bTxu+sauSW5uLrKzs7Fhwwa1p+o+//xzwfLvfvc7zJw5E7t27cLKlSvh6empR+8NY0jsq1evFiwvWLAAYWFhiI+Px759+7Bx40bVa5jyuHd3d+Po0aN45plnMHr0aLX1hYWFguUlS5ZgwoQJeOONN5CQkKD2kMdAp+v79qhjQ3Nz82PsZf/417/+hS1btkAul+P5558XrMvOzhYs//a3v0ViYiIOHjyItWvX4tlnn+3PrhpELpdDLperlufNm4cXXngBU6ZMwfbt25GVlQXAPMY8NzcXbm5umDZtmtq6gXiMf1h5eTmSkpIwadIkvPLKK73WFWs/N/vLlHZ2dqprwA9T3uOi7QtDWa5LWzs7O3R0dGh8HYVCIdqXkr6x93Tu3DksXboUL7zwArZv3/7I+hKJBGvXrkVnZ6deTzYZg7FiV1q0aBHc3d3x1VdfCbZhyuN+9uxZ3L59W+N/yprY2NggOTkZjY2NuHTpku4dHiB0fd8edWwYbEloT+Xl5Zg/fz6Cg4Pxl7/8Rac2yqfsHt4/BqvJkycjPDxcbV8HTHfMKysrcf78ecTFxQnuj9JmIBzjle7cuYNZs2bB2dlZdb9sb8Taz80+GRs+fLjGuXKUZdoexx4yZAikUqlObYcPH46uri7U1tYK6nV0dKC+vl60R771jf1hV69exZw5cxAcHIy8vDyddlQAqtP7DQ0Nfeix8Rgj9p68vLwE8ZjyuAMPTvtbWFhg4cKFOm9b7HE3hK7vm/Kyhba6Yo27MVRVVWH69OlwdnZGfn4+HB0ddWo3mMddE037OmCaYw5Ada+nrv94AQNjzJuamjBz5kw0Njbiyy+/1GkcxNrPzT4ZCw0NRUVFhdopxaKiItV6TSwsLDB27FgUFxerrSsqKoK/v7/qQKV8jZ51i4uL0d3drXUbj5u+sStdv34dM2bMwLBhw5Cfny+Yc+tRlBPjurm59a3TRmJo7D0REW7evCmIx1THHXjw3+Ann3yCqKioPh10xB53Q4SGhqKkpATd3d2C8qKiIshkMtWcStrGvbq6Gj/++KNo426o+vp6TJ8+He3t7fjHP/6h8V4ZbQbzuGtSWVmp074+2MdcKTc3FwEBAX26xCz2mCsUCsyePRsVFRU4ceIExowZo1M70fbzPk2EYYK+/fZbtTmXFAoFjRgxgsLDw1Vlt27dUps35K233lKbg6a8vJwsLS0pJSVFVXbv3j0aMmQIvfTSS4L2ixcvJplMRvX19cYOSyeGxF5TU0P+/v7k4eFBN27c0LqN+vp66uzsFJR1dHTQc889RzY2Nqo5cPqbIbHX1taqvV5mZiYBoF27dqnKTHHclT799FMCQNnZ2RrXa3qPmpubKSAggIYOHao2B5AYept3qbq6msrKyqijo0NVdvToUbX5h+rq6sjFxYXi4uIE7Z9++mkKCQkRfPY3bdpEEomESktLjR9MH/U19paWFgoLCyNHR0cqLi7W+rpNTU2kUCgEZd3d3RQXF0cA6NKlS0aLQR99jVvT5/jkyZMEgFatWiUoN7UxVyopKSEAtHnzZo2vOxCP8Z2dnTRnzhyysrKikydPaq03kPZzs0/GiIhiY2PJysqK1q9fTx988AFFRESQlZUVnT17VlUnMjKSeuauyi+XYcOG0Z/+9CfavXs3eXl5kYeHh9pOrPyyjomJoYMHD9LLL79MAGj79u39EqM2+sYeEhJCAGjDhg10+PBhwd+pU6dU9Q4dOkQBAQGUkpJCWVlZtGPHDgoODiYAtGPHjn6LUxN9Y7ezs6OEhAR65513KDMzkxYuXEgSiYRCQ0OptbVVUNfUxl1pwYIFJJVKqbGxUeP6tLQ0CgkJoU2bNtGBAwdoy5Yt5OPjQxKJhD766KPHEpOu9uzZQ9u2baMVK1YQAIqOjqZt27bRtm3bVPG88sorBEDwj0ZnZyc9++yz5ODgQFu2bKHMzEwKCgoiR0dHKi8vF2zjiy++IIlEQs8//zwdOHCAVq1aRRYWFvT73/++P0NVo2/sc+fOJQC0ZMkStf39+PHjqnoFBQXk7u5Oa9eupczMTPrzn/9Mzz33HAGgxMTEfo72f/SNe8SIERQbG0sZGRmUlZVFiYmJZGVlRV5eXnTnzh3BNkxtzJXWrVtHANQ+40oD8Ri/evVqAkCzZ89W+7wePnxYVW8g7eecjNGD2cdff/11cnd3J6lUShMnTqQvv/xSUEfbF1NVVRXFxMSQk5MTOTg40EsvvUTXrl3TuJ0DBw7QqFGjyMbGhgICAmj37t3U3d39WGLSlb6xA9D6FxkZqapXXFxMs2fPpqeeeopsbGzIwcGBJk+eTB9//HF/hNcrfWNftmwZjRkzhhwdHcna2ppGjBhBKSkp1NzcrHE7pjTuRA/Oftja2lJ0dLTW1z916hRNmzaN3N3dydramlxcXGj69Ol05swZo8fSVz4+Plo/u8qDsrYvp4aGBlq6dCm5urqSTCajyMhIrbOzHz9+nEJDQ0kqlZKnpydt2rRJ45mH/qRv7L218/HxUdWrrKyk2NhY8vX1JVtbW5LJZDRhwgTKysoS9TOvb9ypqakUGhpKzs7OZG1tTd7e3rRixQq1REzJlMaciKirq4ueeuopGj9+vNbXH4jHeOWxS9uf0kDazyVEWqYPZ4wxxhhjj53Z38DPGGOMMSYmTsYYY4wxxkTEyRhjjDHGmIg4GWOMMcYYExEnY4wxxhhjIuJkjDHGGGNMRJyMMcYYY4yJiJMxxhhjjDERcTLGGGOMMSYiTsYYY2YtPT0dEokEd+/e7bWer68vEhIS+qdTjDGzwskYY4wxxpiIrMTuAGOMDQb//e9/YWHB/78yxoyPkzHGGNOBVCoVuwuMMRPF/+YxxhiAu3fvQi6Xw8nJCa6urli9ejUUCoVqfc97xnJyciCRSFBYWIjXXnsNbm5usLe3x/z581FXVydCBIyxwYqTMcYYAyCXy6FQKLBz5068+OKLeP/995GYmPjIditXrsTVq1eRlpaGFStW4IsvvkBycnI/9JgxZir4MiVjjAHw8/PDZ599BgBISkqCk5MT9u3bh9dffx3jxo3T2s7V1RWnTp2CRCIBAHR3d+P9999HU1MTnJ2d+6XvjLHBjc+MMcYYHiRgD1u5ciUAID8/v9d2iYmJqkQMAH7961+jq6sLt27dMn4nGWMmiZMxxhgDEBgYKFgOCAiAhYUFbt682Ws7b29vwfITTzwBAPj555+N2j/GmOniZIwxxjR4+GxXbywtLTWWE5Exu8MYM2GcjDHGGIBr164Jlr///nt0d3fD19dXnA4xxswGJ2OMMQYgMzNTsLxnzx4AwMyZM8XoDmPMjPDTlIwxBuDGjRuYM2cOZsyYgfPnz+Ojjz7CokWLEBISInbXGGMmjs+MMcYYgGPHjkEqlWLjxo04efIkkpOTkZ2dLXa3GGNmQEJ8lyljjDHGmGj4zBhjjDHGmIg4GWOMMcYYExEnY4wxxhhjIuJkjDHGGGNMRJyMMcYYY4yJiJMxxhhjjDERcTLGGGOMMSYiTsYYY4wxxkTEyRhjjDHGmIg4GWOMMcYYExEnY4wxxhhjIuJkjDHGGGNMRP8HvghrkIc/06sAAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 600x600 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"def get_fit_results(model, params):\n",
" # start with default pre-fit parameter values for everything\n",
" param_values = cabinetry.model_utils.asimov_parameters(model)\n",
" # corresponding default pre-fit uncertainties\n",
" param_uncertainty = cabinetry.model_utils.prefit_uncertainties(model)\n",
" # diagonal correlation matrix\n",
" corr_mat = np.zeros(shape=(len(param_values), len(param_values)))\n",
" np.fill_diagonal(corr_mat, 1.0)\n",
" \n",
" # update central values of parameters as specified in params dict\n",
" labels = model.config.par_names\n",
" for k, v in params.items():\n",
" idx = cabinetry.model_utils._parameter_index(k, labels) # index of parameter\n",
" param_values[idx] = v\n",
" \n",
" return cabinetry.fit.FitResults(param_values, param_uncertainty, labels, corr_mat, np.nan)\n",
"\n",
"\n",
"# get a custom fit_results object that contains all relevant information\n",
"custom_parameters = get_fit_results(model, params)\n",
"\n",
"# produce visualizations of the pre-fit model and observed data\n",
"prefit_model = cabinetry.model_utils.prediction(model, fit_results = custom_parameters, label=\"custom_parameters\")\n",
"cabinetry.visualize.data_mc(prefit_model, data)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "951b3992",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "stats",
"language": "python",
"name": "stats"
},
"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.16"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
{
"channels": [
{
"name": "singlechannel",
"samples": [
{
"name": "signal",
"data": [
5.0,
10.0
],
"modifiers": [
{
"name": "mu",
"type": "normfactor",
"data": null
}
]
},
{
"name": "background",
"data": [
50.0,
60.0
],
"modifiers": [
{
"name": "uncorr_bkguncrt",
"type": "shapesys",
"data": [
5.0,
12.0
]
}
]
}
]
}
],
"observations": [
{
"name": "singlechannel",
"data": [
53.0,
65.0
]
}
],
"measurements": [
{
"name": "Measurement",
"config": {
"poi": "mu",
"parameters": []
}
}
],
"version": "1.0.0"
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment