{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Lifetimes" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%pip install matplotlib numpy scipy pairinteraction\n", "\n", "from typing import Any\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pairinteraction.real as pi\n", "from scipy.optimize import curve_fit" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "if pi.Database.get_global_database() is None:\n", " pi.Database.initialize_global_database(download_missing=True)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Lifetime at T=0: 241.79 microsecond\n", "Lifetime at T=300K: 101.60 microsecond\n" ] } ], "source": [ "ket = pi.KetAtom(\"Rb\", n=60, l=0, j=0.5, m=0.5)\n", "\n", "T_in_K = 300\n", "lifetime_0 = ket.get_lifetime()\n", "lifetime = ket.get_lifetime(T_in_K, temperature_unit=\"K\")\n", "\n", "print(f\"Lifetime at T=0: {lifetime_0.to('mus'):.2f}\")\n", "print(f\"Lifetime at T={T_in_K}K: {lifetime.to('mus'):.2f}\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of possible spontaneous decay transitions: 275\n", "Number of considered BBR transitions: 430\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj8AAAHcCAYAAAA5lMuGAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAd/BJREFUeJzt3X1cjff/B/DXqZxzulGhdMpSTc09UUpmi2nLzbYyM8zI/WaY1szYyN2sjZn7r2YbsQ1hmGFI7FZCKVIMi4xu3FUrKjqf3x9+XXN0ujnHSXFez8fjPDif6319rvf5nHPVu+t8ruuSCSEEiIiIiIyESW0nQERERPQwsfghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMCosfIqozZs6cCZlMVttpGK2oqCjIZDKcP3/eoLFEdQ2LH6pxJ06cwKuvvgoXFxcolUo0adIEzz//PJYuXVrbqQEADh48iJkzZyI3N7e2U6FqKvvFe++jcePG6N69O37++We9Y7XFm5mZoUmTJhg2bBguXbpUZW73b6uixy+//GKo4ahRdXH/eND3qCJ18bVSzTCr7QTo8Xbw4EF0794dTZs2xejRo6FSqXDx4kUcOnQIixcvxoQJE2o7RRw8eBCzZs3CsGHDYGtrW9vpkA5mz54NNzc3CCGQnZ2NqKgo9O7dGz/99BNefPFFvWPvjS8qKsKhQ4cQFRWFP/74AykpKVAqlRXm9O2332o8X7t2LWJiYsq1t2zZ8gFeec0YMmQIBg4cCIVCIbVVtH9oi33Y9H2PKsKfBcaDxQ/VqLlz58LGxgZHjhwp98MkJyendpKix0avXr3g7e0tPR85ciQcHBywfv36cgWNLrH3x48aNQp2dnb47LPPsH37drz22msV5vTGG29oPD906BBiYmLKtVeksLAQlpaW1Yo1NFNTU5iamho8tqbo+x4R8WsvqlHnzp1D69attf4V1bhxY+n/ZXM9Tp06hddeew3W1tZo1KgRJk6ciKKionLrHjt2DL169YK1tTWsrKzQo0cPHDp0SCOmrM+zZ89Kf8nZ2Nhg+PDhuHnzphTz/vvvAwDc3Nykw+hl8xguXLiAt99+G82bN4e5uTkaNWqE/v37l5vnUJ1t3evSpUsYMWIEHBwcoFAo0Lp1a6xatUqv1zls2DC4urqWW/f++TP//vsvQkND4erqCoVCgcaNG+P5559HYmJiuXXvVd0x0HUc/vjjD3Tq1AlKpRLNmjXDl19+WWke1WFrawtzc3OYmVX9d50usQDwzDPPALj7mTaUsvFKTU3F66+/jgYNGqBr164ADP/Zq877f/88nsr2j4rm/Bhq39SHtvdIl3Gs7GdBdfZZffaxzZs3QyaT4ddffy237Msvv4RMJkNKSore/Zfp3r07nn32WSQmJqJXr16oX78+mjRpgsWLF1e57uOIR36oRrm4uCAuLg4pKSlo06ZNlfGvvfYaXF1dERERgUOHDmHJkiW4ceMG1q5dK8WcPHkSzzzzDKytrTF58mTUq1cPX375Jbp164Zff/0Vvr6+5fp0c3NDREQEEhMT8fXXX6Nx48b47LPP8Morr+Cvv/7C+vXrsXDhQtjZ2QEA7O3tAQBHjhzBwYMHMXDgQDzxxBM4f/48VqxYgW7duiE1NRUWFhbV3laZ7OxsdO7cGTKZDOPHj4e9vT1+/vlnjBw5Evn5+QgNDdXrdVblrbfewubNmzF+/Hi0atUK165dwx9//IG0tDR07NixwvV0HYPqjMOJEyfwwgsvwN7eHjNnzsSdO3cwY8YMODg46PSa8vLycPXqVQghkJOTg6VLl6KgoEDrURZdYrUp+yXYoEEDnXKsjv79+8PDwwOffPIJhBAADP/Z0+f9r2r/uJ8h9019aHuPqjuOlb3W6u6z+oxxnz59YGVlhY0bN8Lf319jWXR0NFq3bi397NR3Hwbu7nNOTk546aWXMHz4cAQHB+Orr77Cu+++i+eeew5t27bVebwfaYKoBu3du1eYmpoKU1NT4efnJyZPniz27NkjSkpKNOJmzJghAIiXX35Zo/3tt98WAERycrLUFhwcLORyuTh37pzUdvnyZVG/fn3x7LPPlutzxIgRGn327dtXNGrUSHo+f/58AUCkp6eXy//mzZvl2uLi4gQAsXbtWp23JYQQI0eOFI6OjuLq1asa7QMHDhQ2NjbSNqv7OkNCQoSLi0u5PMtyKmNjYyPGjRtXLq4q1R2De7dZ1TgEBwcLpVIpLly4ILWlpqYKU1NTUZ0fS6tXrxYAyj0UCoWIiorSO/be+H379okrV66Iixcvis2bNwt7e3uhUCjExYsXq8zvXuPGjavwNZWN16BBg8otM/Rnrzrvf9lrv3dfqGj/0BZr6H2zqjyr8x7p8vmt6LVWd5/Vdx8bNGiQaNy4sbhz547UlpmZKUxMTMTs2bOlNn37v3z5sgAg7O3tNcYmNTVVABBr1qzRuc9HHb/2ohr1/PPPIy4uDi+//DKSk5Mxb948BAYGokmTJti+fXu5+HHjxmk8L5sQvWvXLgBAaWkp9u7di+DgYDz55JNSnKOjI15//XX88ccfyM/P1+jjrbfe0nj+zDPP4Nq1a+XitDE3N5f+f/v2bVy7dg3u7u6wtbXVeqi5qm0JIfDDDz/gpZdeghACV69elR6BgYHIy8tDYmKiXq+zKra2toiPj8fly5d1Wk/XMQAqH4fS0lLs2bMHwcHBaNq0qRTTsmVLBAYG6pTb8uXLERMTg5iYGHz33Xfo3r07Ro0ahS1btjxQLAAEBATA3t4ezs7OePXVV2FpaYnt27fjiSee0CnH6rh/vADDf/b0ff+r62Hvm0D13iN9Pr/3qu4+C+g/xgMGDEBOTo7GGYCbN2+GWq3GgAEDpDZ9+z9x4gQAYMaMGRpjU69ePQCAXC7Xqb/HAYsfqnGdOnXCli1bcOPGDRw+fBhTp07Fv//+i1dffRWpqakasR4eHhrPmzVrBhMTE+lw9pUrV3Dz5k00b9683HZatmwJtVqNixcvarTf+wsW+O+Q+I0bN6rM/datWwgPD4ezszMUCgXs7Oxgb2+P3Nxc5OXllYuvaltXrlxBbm4uVq5cCXt7e43H8OHDAdydCK7P66zKvHnzkJKSAmdnZ/j4+GDmzJn4+++/DT4GVY3DlStXcOvWrXLvNQCtr7cyPj4+CAgIQEBAAAYPHoydO3eiVatWGD9+PEpKSvSOBf4rljZv3ozevXvj6tWrNXZmk5ubW7k2Q3/29H3/q+th75tA9d4jfT6/97+u6uyzgP5j3LNnT9jY2CA6Olpqi46OhqenJ5566impTd/+y4qf4OBgjfZTp04B0H2/exyw+KGHRi6Xo1OnTvjkk0+wYsUK3L59G5s2bap0HUNc8K6iM1LE/8+tqMyECRMwd+5cvPbaa9i4cSP27t2LmJgYNGrUCGq1Wudtla3zxhtvSEch7n88/fTT1X1pACoeo9LSUo3nr732Gv7++28sXboUTk5OmD9/Plq3bq31Wjf30nUMgAcb8wdhYmKC7t27IzMzE2fOnHmg2LJiqV+/fti+fTvatGmD119/HQUFBQbP+96jE2UM/dnT9/2vSQ/6OanOe6TP5/deuuyz+o6xQqFAcHAwtm7dijt37uDSpUv4888/NY76PEj/x48fh0qlQpMmTTTak5OTYWZmhlatWlU5Do8bTnimWlF2empmZqZG+5kzZzT+Cj579izUarV0NpO9vT0sLCxw+vTpcn2eOnUKJiYmcHZ21imXygqszZs3IyQkBAsWLJDaioqK9L4Imr29PerXr4/S0lIEBARUGFdaWlrt19mgQQOt+Vy4cKFcm6OjI95++228/fbbyMnJQceOHTF37lz06tWrwlxqYgzMzc21FhzaXq+u7ty5AwDVKlKqG2tqaoqIiAh0794dy5Ytw5QpUx44z6oYetwB/d7/6v4BUhP7pi4qeo90GUdtr7W6+2wZfcYYuPvV15o1axAbG4u0tDQIIcoVP/r2f+LECbRv375c+/Hjx/HUU0/V6rWaaguP/FCNOnDggNa/4srm8Nx/uHX58uUaz8uuAl22Y5uamuKFF17Ajz/+qHGqanZ2NtatW4euXbvC2tpapxzLrqmi7YehqalpufyXLl1a7qhKdZmamqJfv3744YcfpNNX73XlyhUprrqvs1mzZsjLy8Px48eluMzMTGzdulV6XlpaWu4Qf+PGjeHk5ITi4uIqczb0GAQGBmLbtm3IyMiQ2tPS0rBnzx69+ixz+/Zt7N27F3K5vMqLCOoSCwDdunWDj48PFi1apPXyC4ZmyHF/kPe/sv3j/nwNvW/qStt7pMs4anut1d1nH2SMgbvzlxo2bIjo6GhER0fDx8dH4w9BffsvLS1FWlqa1uInOTkZ7dq1qzK3xxGP/FCNmjBhAm7evIm+ffuiRYsWKCkpwcGDBxEdHQ1XV1fpO/My6enpePnll9GzZ0/ExcXhu+++w+uvv66x43788ceIiYlB165d8fbbb8PMzAxffvkliouLMW/ePJ1z9PLyAgB89NFHGDhwIOrVq4eXXnoJlpaWePHFF/Htt9/CxsYGrVq1QlxcHPbt24dGjRrpPSaffvopDhw4AF9fX4wePRqtWrXC9evXkZiYiH379uH69es6vc6BAwfigw8+QN++ffHOO+/g5s2bWLFiBZ566ilpIua///6LJ554Aq+++irat28PKysr7Nu3D0eOHNH4i1ibmhiDWbNmYffu3XjmmWfw9ttv486dO1i6dClat26tUcRV5eeff5bmLeTk5GDdunU4c+YMpkyZUu4XrS6xFXn//ffRv39/REVFaZ2kbEiGHPcHef8r2j+0MfS+qY/73yNdxrGi11qdffZBxhi4O/n4lVdewYYNG1BYWIjPP/9cY7m+/Z85cwZFRUXlip9bt27h7NmzCAkJqc6wPn5q4QwzMiI///yzGDFihGjRooWwsrIScrlcuLu7iwkTJojs7GwpruzU19TUVPHqq6+K+vXriwYNGojx48eLW7dules3MTFRBAYGCisrK2FhYSG6d+8uDh48qBFT1ueVK1c02rWdojtnzhzRpEkTYWJiorHsxo0bYvjw4cLOzk5YWVmJwMBAcerUKeHi4iJCQkL02pYQQmRnZ4tx48YJZ2dnUa9ePaFSqUSPHj3EypUrdX6dQty9pECbNm2EXC4XzZs3F999953Gqe7FxcXi/fffF+3btxf169cXlpaWon379uJ///tfub7uV90x0HUcfv31V+Hl5SXkcrl48sknRWRkZLnT8yui7fR1pVIpPD09xYoVK4RardYr9t74I0eOlNtuaWmpaNasmWjWrJnGacmVqc6p7vePlxCG/exV9/2v6POqbf+oKNbQ+6Y2urxHunx+K3qtQlS9zz7IPlYmJiZGABAymazcJRX07X/jxo0CgEhJSdFoP3z4sAAgduzYUe38HicyIWp4BiJRNcycOROzZs3ClStXpIuLERER1QTO+SEiIiKjwuKHiIiIjAqLHyIiIjIqnPNDRERERoVHfoiIiMiosPghIiIio8KLHNYharUaly9fRv369Q1yTysiIiJjIYTAv//+CycnJ5iYVH5sh8VPHXL58uUavfcNERHR4+7ixYt44oknKo1h8VOH1K9fH8DdN66m74FDRET0OMnPz4ezs7P0u7QyLH7qkLKvuqytrVn8EBER6aE600Y44ZmIiIiMSp0sfpYvXw5XV1colUr4+vri8OHDlcZv2rQJLVq0gFKpRNu2bbFr1y6N5UIIhIeHw9HREebm5ggICMCZM2ek5efPn8fIkSPh5uYGc3NzNGvWDDNmzEBJSYlGP8ePH8czzzwDpVIJZ2dnrXcprioXIiIiql11rviJjo5GWFgYZsyYgcTERLRv3x6BgYHIycnRGn/w4EEMGjQII0eOxLFjxxAcHIzg4GCkpKRIMfPmzcOSJUsQGRmJ+Ph4WFpaIjAwEEVFRQCAU6dOQa1W48svv8TJkyexcOFCREZG4sMPP5T6yM/PxwsvvAAXFxckJCRg/vz5mDlzJlauXKlTLkRERFS76twVnn19fdGpUycsW7YMwN3Tv52dnTFhwgRMmTKlXPyAAQNQWFiIHTt2SG2dO3eGp6cnIiMjIYSAk5MT3nvvPUyaNAkAkJeXBwcHB0RFRWHgwIFa85g/fz5WrFiBv//+GwCwYsUKfPTRR8jKyoJcLgcATJkyBdu2bcOpU6eqlUtV8vPzYWNjg7y8vArn/AghcOfOHZSWllbZH9HjzNTUFGZmZrwsBBEBqN7v0DJ1asJzSUkJEhISMHXqVKnNxMQEAQEBiIuL07pOXFwcwsLCNNoCAwOxbds2AEB6ejqysrIQEBAgLbexsYGvry/i4uIqLH7y8vLQsGFDje08++yzUuFTtp3PPvsMN27cQIMGDarM5X7FxcUoLi6Wnufn52uNK1NSUoLMzEzcvHmz0jgiY2FhYQFHR0eN/ZKIqCp1qvi5evUqSktL4eDgoNHu4OAgHV25X1ZWltb4rKwsaXlZW0Ux9zt79iyWLl2Kzz//XGM7bm5u5fooW9agQYMqc7lfREQEZs2apXXZ/dRqNdLT02FqagonJyfI5XL+xUtGSwiBkpISXLlyBenp6fDw8KjyomZERGXqVPFTF1y6dAk9e/ZE//79MXr06Brd1tSpUzWOFJVdo0CbkpIS6StACwuLGs2L6FFgbm6OevXq4cKFCygpKYFSqaztlIjoEVGn/lSys7ODqakpsrOzNdqzs7OhUqm0rqNSqSqNL/u3On1evnwZ3bt3R5cuXTQmMle2nXu3UVUu91MoFNI1fap7bR/+dUv0H+4PRKSPOvWTQy6Xw8vLC7GxsVKbWq1GbGws/Pz8tK7j5+enEQ8AMTExUrybmxtUKpVGTH5+PuLj4zX6vHTpErp16wYvLy+sXr263A9VPz8//Pbbb7h9+7bGdpo3b44GDRpUKxciIiKqA0Qds2HDBqFQKERUVJRITU0VY8aMEba2tiIrK0sIIcSQIUPElClTpPg///xTmJmZic8//1ykpaWJGTNmiHr16okTJ05IMZ9++qmwtbUVP/74ozh+/LgICgoSbm5u4tatW0IIIf755x/h7u4uevToIf755x+RmZkpPcrk5uYKBwcHMWTIEJGSkiI2bNggLCwsxJdffqlTLpXJy8sTAEReXl65Zbdu3RKpqalSzkTE/YKI/lPZ79D71bniRwghli5dKpo2bSrkcrnw8fERhw4dkpb5+/uLkJAQjfiNGzeKp556SsjlctG6dWuxc+dOjeVqtVpMnz5dODg4CIVCIXr06CFOnz4tLV+9erUAoPVxr+TkZNG1a1ehUChEkyZNxKefflou96pyqQyLHyLdcL8gojK6FD917jo/xqyyaxQUFRUhPT0dbm5u5SZ2ens/zCyBo0d1i79y5QrCw8Oxc+dOZGdno0GDBmjfvj3Cw8Px9NNP10ySFXB1dUVoaChCQ0Mf6nZry7Bhw7BmzRoAgJmZGRo2bIh27dph0KBBGDZs2CM/Z6ay/YKIjMsje50fejz169cPJSUlWLNmDZ588klkZ2cjNjYW165dq+3UjELPnj2xevVqlJaWIjs7G7t378bEiROxefNmbN++HWZm/DFARMbl0f6zj+q83Nxc/P777/jss8/QvXt3uLi4wMfHB1OnTsXLL78sxclkMqxYsQK9evWCubk5nnzySWzevFmjrxMnTuC5556Dubk5GjVqhDFjxqCgoEBaPmzYMAQHB+Pzzz+Ho6MjGjVqhHHjxkmT1Lt164YLFy7g3XffhUwmk66TdO3aNQwaNAhNmjSBhYUF2rZti/Xr12tsu1u3bnjnnXcwefJkNGzYECqVCjNnziz3WkeNGgV7e3tYW1vjueeeQ3JyskbMihUr0KxZM8jlcjRv3hzffvuttOz8+fOQyWRISkrS6FMmk+GXX34BANy4cQODBw+Gvb09zM3N4eHhgdWrV1f6HigUCqhUKjRp0gQdO3bEhx9+iB9//BE///wzoqKidMr/p59+QqdOnaBUKmFnZ4e+fftKy7799lt4e3ujfv36UKlUeP3116Xb0ggh4O7urnHtLABISkqCTCbD2bNnK30NRESGxOKHapSVlRWsrKywbds2jatZazN9+nT069cPycnJGDx4MAYOHIi0tDQAQGFhIQIDA9GgQQMcOXIEmzZtwr59+zB+/HiNPg4cOIBz587hwIEDWLNmDaKioqRf8Fu2bMETTzyB2bNnIzMzE5mZmQDufnXi5eWFnTt3IiUlBWPGjMGQIUPK3VB3zZo1sLS0RHx8PObNm4fZs2cjJiZGWt6/f3/k5OTg559/RkJCAjp27IgePXrg+vXrAICtW7di4sSJeO+995CSkoI333wTw4cPx4EDB6o9ntOnT0dqaip+/vlnpKWlYcWKFbCzs6v2+mWee+45tG/fHlu2bKl2/jt37kTfvn3Ru3dvHDt2DLGxsfDx8ZHWv337NubMmYPk5GRs27YN58+fx7BhwwDcLW5HjBhRrlBbvXo1nn32Wbi7u+v8GoiI9FbTE5Co+vSd8Ozl9XAfutq8ebNo0KCBUCqVokuXLmLq1KkiOTlZIwaAeOuttzTafH19xdixY4UQQqxcuVI0aNBAFBQUSMt37twpTExMpDMBQ0JChIuLi7hz544U079/fzFgwADpuYuLi1i4cGGVOffp00e899570nN/f3/RtWtXjZhOnTqJDz74QAghxO+//y6sra1FUVGRRkyzZs2kMwK7dOkiRo8erbG8f//+onfv3kIIIdLT0wUAcezYMWn5jRs3BABx4MABIYQQL730khg+fHiV+ZcJCQkRQUFBWpcNGDBAtGzZstr5+/n5icGDB1d720eOHBEAxL///iuEEOLSpUvC1NRUxMfHCyGEKCkpEXZ2diIqKqrafd6PE54fsgf9YUBUg3SZ8MwjP1Tj+vXrh8uXL2P79u3o2bMnfvnlF3Ts2FHjKxcA5a6H5OfnJx35SUtLQ/v27WFpaSktf/rpp6FWq3H69GmprXXr1jA1NZWeOzo6Sl+9VKS0tBRz5sxB27Zt0bBhQ1hZWWHPnj3IyMjQiGvXrp3G83v7Tk5ORkFBARo1aiQd7bKyskJ6ejrOnTsnvYb7J3g//fTT0musjrFjx2LDhg3w9PTE5MmTcfDgwWqvez8hhPTVX3XyT0pKQo8ePSrsLyEhAS+99BKaNm2K+vXrw9/fHwCkcXRyckKfPn2watUqAHe/QisuLkb//v31fg1ERPrgTEd6KJRKJZ5//nk8//zzmD59OkaNGoUZM2ZIX4sYSr169TSey2QyqNXqSteZP38+Fi9ejEWLFqFt27awtLREaGgoSkpKqt13QUEBHB0dpbk597K1ta1W7mVnXol7TsC896KaANCrVy9cuHABu3btQkxMDHr06IFx48aVm0tTHWlpadL96qqTv7m5eYV9lX0tGRgYiO+//x729vbIyMhAYGCgxjiOGjUKQ4YMwcKFC7F69WoMGDCAt2shooeOR36oVrRq1QqFhYUabYcOHSr3vGXLlgCAli1bIjk5WWOdP//8EyYmJmjevHm1tyuXy1FaWqrR9ueffyIoKAhvvPEG2rdvjyeffBJ//fWXTq+nY8eOyMrKgpmZGdzd3TUeZXNyWrZsiT///LPctlu1agUAsLe3BwBpLhIAjcnPZezt7RESEoLvvvsOixYtKncrlurYv38/Tpw4gX79+lU7/3bt2pW7gnmZU6dO4dq1a/j000/xzDPPoEWLFlqPuPXu3RuWlpZYsWIFdu/ejREjRuicOxHRg2LxQzXq2rVreO655/Ddd9/h+PHjSE9Px6ZNmzBv3jwEBQVpxG7atAmrVq3CX3/9hRkzZuDw4cPShObBgwdDqVQiJCQEKSkpOHDgACZMmIAhQ4bAwcGh2vm4urrit99+w6VLl3D16lUAgIeHB2JiYnDw4EGkpaXhzTffLHePtqoEBATAz88PwcHB2Lt3L86fP4+DBw/io48+wtH/vzDS+++/j6ioKKxYsQJnzpzBF198gS1btmDSpEkA7h5Z6dy5Mz799FOkpaXh119/xbRp0zS2Ex4ejh9//BFnz57FyZMnsWPHDqlArEhxcTGysrJw6dIlJCYm4pNPPkFQUBBefPFFDB06tNr5z5gxA+vXr8eMGTOQlpaGEydO4LPPPgMANG3aFHK5HEuXLsXff/+N7du3Y86cOeVyMTU1xbBhwzB16lR4eHjw1i9EVDtqfAYSVdvjeIXnoqIiMWXKFNGxY0dhY2MjLCwsRPPmzcW0adPEzZs3pTgAYvny5eL5558XCoVCuLq6iujoaI2+jh8/Lrp37y6USqVo2LChGD16tDSZVgjtk3snTpwo/P39pedxcXGiXbt2QqFQSFfwvnbtmggKChJWVlaicePGYtq0aWLo0KEaffn7+4uJEydq9B0UFKRxtfH8/HwxYcIE4eTkJOrVqyecnZ3F4MGDRUZGhhTzv//9Tzz55JOiXr164qmnnhJr167V6DM1NVX4+fkJc3Nz4enpKfbu3asx4XnOnDmiZcuWwtzcXDRs2FAEBQWJv//+u8LxDwkJka5WbmZmJuzt7UVAQIBYtWqVKC0t1YitTv4//PCD8PT0FHK5XNjZ2YlXXnlFWrZu3Trh6uoqFAqF8PPzE9u3by83gVsIIc6dOycAiHnz5lWYd3U9qvvFI4sTnqkO4xWeH1H6XuH5cSCTybB161YEBwfXdipUw37//Xf06NEDFy9e1OmonTaP+35R59x/OXldL/dOVIN4hWciqnOKi4tx5coVzJw5E/3793/gwoeISF+c80NED8X69evh4uKC3NxczJs3r7bTISIjxiM/VCfw29fH37Bhwwx+aQOqA+79Koxfg9Ejgkd+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPCU90fB/dfdbWm3Xc6a2rqf////3t06uz8+fNwc3PDsWPH4OnpqX9u96hLV42OiopCaGgocnNzazsVrVxdXREaGorQ0FAAhhu7uvQeEBGV4ZEfqnHDhg2DTCaTHo0aNULPnj1x/Pjx2k5NQ1RUlEaeVlZW8PLywpYtW2o7NQBAt27dpNyUSiWeeuopRERE1Mg1kjIzM9GrV69qx8+cOVNr0aprP0REDwOLH3ooevbsiczMTGRmZiI2NhZmZmZ48cUXazutcqytraU8jx07hsDAQLz22ms4ffp0bacGABg9ejQyMzNx+vRpTJ06FeHh4YiMjDT4dlQqFRQKRZ3ph4jIkFj80EOhUCigUqmgUqng6emJKVOm4OLFi7hy5YrW+NLSUowcORJubm4wNzdH8+bNsXjx4nJxq1atQuvWraFQKODo6Ijx48dXmMOMGTPg6OhY6REnmUwm5enh4YGPP/4YJiYmGuvcuHEDQ4cORYMGDWBhYYFevXrhzJkzGv1ERUWhadOmsLCwQN++fXHt2jVp2fnz52FiYoKj9319uGjRIri4uECtVleYn4WFBVQqFVxcXDB8+HC0a9cOMTEx0vJz584hKCgIDg4OsLKyQqdOnbBv3z6NPnJycvDSSy/B3Nwcbm5u+P7777WOw7Zt26TnH3zwAZ566ilYWFjgySefxPTp03H79m3ptc6aNQvJycnSkamoqCit/Zw4cQLPPfcczM3N0ahRI4wZMwYFBQXS8mHDhiE4OBiff/45HB0d0ahRI4wbN07aFhGRIbD4oYeuoKAA3333Hdzd3dGoUSOtMWq1Gk888QQ2bdqE1NRUhIeH48MPP8TGjRulmBUrVmDcuHEYM2YMTpw4ge3bt8Pd3b1cX0IITJgwAWvXrsXvv/+Odu3aVSvP0tJSrFmzBgDQsWNHqX3YsGE4evQotm/fjri4OAgh0Lt3b+kXdHx8PEaOHInx48cjKSkJ3bt3x8cffyyt7+rqioCAAKxevVpje6tXr8awYcNgYlL1bimEwO+//45Tp05BLpdL7QUFBejduzdiY2Nx7Ngx9OzZEy+99BIyMjI08r948SIOHDiAzZs343//+x9ycnIq3V79+vURFRWF1NRULF68GF999RUWLlwIABgwYADee+89tG7dWjpqNmDAgHJ9FBYWIjAwEA0aNMCRI0ewadMm7Nu3r1zBeuDAAZw7dw4HDhzAmjVrEBUVJRVTREQGIajOyMvLEwBEXl5euWW3bt0Sqamp4tatW+VX9PJ6uI/7nDz530ObkJAQYWpqKiwtLYWlpaUAIBwdHUVCQoIUk56eLgCIY8eOVTg+48aNE/369ZOeOzk5iY8++qjCeABi06ZN4vXXXxctW7YU//zzT4WxQgixevVqAUDK08TERCgUCrF69Wop5q+//hIAxJ9//im1Xb16VZibm4uNGzcKIYQYNGiQ6N27t0bfAwYMEDY2NtLz6Oho0aBBA1FUVCSEECIhIUHIZDKRnp5eYX7+/v6iXr16wtLSUtSrV08AEEqlUiMXbVq3bi2WLl0qhBDi9OnTAoA4fPiwtDwtLU0AEAsXLpTaAIitW7dW2Of8+fOF1z2fhRkzZoj27duXi7u3n5UrV4oGDRqIgoICafnOnTuFiYmJyMrKEkLc/ay4uLiIO3fuSDH9+/cXAwYM0JpHpfsFGZ62nwWV/Gwgepgq+x16Px75oYeie/fuSEpKQlJSEg4fPozAwED06tULFy5cqHCd5cuXw8vLC/b29rCyssLKlSulIxg5OTm4fPkyevToUel23333XcTHx+O3335DkyZNqsyzfv36Up7Hjh3DJ598grfeegs//fQTACAtLQ1mZmbw9fWV1mnUqBGaN2+OtLQ0Kebe5QDg5+en8Tw4OBimpqbYunUrgLtfHXXv3h2urq6V5jd48GAkJSXhzz//RK9evfDRRx+hS5cu0vKCggJMmjQJLVu2hK2tLaysrJCWliaNW1n+Xl5e0jotWrSAra1tpduNjo7G008/DZVKBSsrK0ybNk3jaFJ1pKWloX379rC0tJTann76aajVao05Va1bt4apqan03NHRscojU0REumDxQw+FpaUl3N3d4e7ujk6dOuHrr79GYWEhvvrqK63xGzZswKRJkzBy5Ejs3bsXSUlJGD58OEpKSgAA5ubm1dru888/j0uXLmHPnj3VijcxMZHybNeuHcLCwtCtWzd89tln1Xuh1SSXyzF06FCsXr0aJSUlWLduHUaMGFHlejY2NtIYbty4EcuWLdOY0zNp0iRs3boVn3zyCX7//XckJSWhbdu20rjpIy4uDoMHD0bv3r2xY8cOHDt2DB999NED9VmZevXqaTyXyWSVzoMiItIVix+qFTKZDCYmJrh165bW5X/++Se6dOmCt99+Gx06dIC7uzvOnTsnLa9fvz5cXV0RGxtb6XZefvllrFu3DqNGjcKGDRv0ytXU1FTKs2XLlrhz5w7i4+Ol5deuXcPp06fR6v8vctSyZUuN5QBw6NChcv2OGjUK+/btw//+9z/cuXMHr7zyik55WVlZYeLEiZg0aZJ0uvuff/6JYcOGoW/fvmjbti1UKhXOnz8vrdOiRQvcuXMHCQkJUtvp06crvf7QwYMH4eLigo8++gje3t7w8PAod8ROLpejtLS00nxbtmyJ5ORkFBYWSm1//vknTExM0Lx5cx1eORHRg2HxQw9FcXExsrKykJWVhbS0NEyYMAEFBQV46aWXtMZ7eHjg6NGj2LNnD/766y9Mnz4dR44c0YiZOXMmFixYgCVLluDMmTNITEzE0qVLy/XVt29ffPvttxg+fDg2b95caZ5CCCnP9PR0rFy5Env27EFQUJCUV1BQEEaPHo0//vgDycnJeOONN9CkSRMp5p133sHu3bvx+eef48yZM1i2bBl2795dblstW7ZE586d8cEHH2DQoEHVPpp1rzfffBN//fUXfvjhBym/LVu2ICkpCcnJyXj99dc1jpo0b94cPXv2xJtvvon4+HgkJCRg1KhRlW7bw8MDGRkZ2LBhA86dO4clS5ZIX9eVcXV1RXp6OpKSknD16lUUFxeX62fw4MFQKpUICQlBSkoKDhw4gAkTJmDIkCFwcHDQ+bUTEemtxmcgUbXpPeG5llVnwjMA6VG/fn3RqVMnsXnzZinm/gnPRUVFYtiwYcLGxkbY2tqKsWPHiilTppSbVBsZGSmaN28u6tWrJxwdHcWECROkZbhv0m50dLRQKpXihx9+0Jpn2YTnsodCoRBPPfWUmDt3rsYE3OvXr4shQ4YIGxsbYW5uLgIDA8Vff/2l0dc333wjnnjiCWFubi5eeukl8fnnn2tMeL43DvdNQK6Iv7+/mDhxYrn2N998U7Ru3VqUlpaK9PR00b17d2Fubi6cnZ3FsmXLyq2XmZkp+vTpIxQKhWjatKlYu3atcHFxqXTC8/vvvy8aNWokrKysxIABA8TChQs1Xk9RUZHo16+fsLW1FQCkSeL393P8+HHRvXt3oVQqRcOGDcXo0aPFv//+Ky0PCQkRQUFBGq9v4sSJwt/fX+uY1OX94rHECc9Uh+ky4VkmRA1cHpb0kp+fDxsbG+Tl5cHa2lpjWVFREdLT0+Hm5galUllLGWpniNtbGKs5c+Zg06ZNde5q14+KurxfPJbuv5XO0aOabfddu4roYarsd+j9+LUXUS0oKChASkoKli1bhgkTJtR2OkRERoXFD1EtGD9+PLy8vNCtW7dqneVFRESGw7u6E9UCXrWYiKj28MgPERERGRUWP48Yzk8n+g/3ByLSR50rfpYvXw5XV1colUr4+vri8OHDlcZv2rQJLVq0gFKpRNu2bbFr1y6N5UIIhIeHw9HREebm5ggICCh3B+65c+eiS5cusLCw0HqZ/6ioKOlu1fc/yi67/8svv2hdnpWV9WAD8v/Krnp78+ZNg/RH9Dgo2x/uvyo0EVFl6tScn+joaISFhSEyMhK+vr5YtGgRAgMDcfr0aTRu3Lhc/MGDBzFo0CBERETgxRdfxLp16xAcHIzExES0adMGADBv3jwsWbIEa9asgZubG6ZPn47AwECkpqZKp8aWlJSgf//+8PPzwzfffFNuOwMGDEDPnj012oYNG4aioqJyeZ0+fVrjFDtteevD1NQUtra2UrFlYWEBmUxmkL4f1L13Higqqr08yHgIIXDz5k3k5OTA1tZW415gRERVqVPX+fH19UWnTp2wbNkyAIBarYazszMmTJiAKVOmlIsfMGAACgsLsWPHDqmtc+fO8PT0RGRkJIQQcHJywnvvvYdJkyYBAPLy8uDg4ICoqCgMHDhQo7+oqCiEhoZWeql/ALhy5QqaNGmCb775BkOGDAFw98hP9+7dcePGjSpvElmRqq5RIP7/6sNV5fewZWb+939Hx9rLg4yPra0tVCpVnflD4LHH6/xQHabLdX7qzJGfkpISJCQkYOrUqVKbiYkJAgICEBcXp3WduLg4hIWFabQFBgZi27ZtAID09HRkZWUhICBAWm5jYwNfX1/ExcWVK36qa+3atbCwsMCrr75abpmnpyeKi4vRpk0bzJw5E08//XSF/RQXF2vcBiA/P7/S7cpkMjg6OqJx48a4ffu2XrnXhP+vKwEA/3+XBaIaV69ePR7xISK91Jni5+rVqygtLS13jx8HBwecOnVK6zpZWVla48vm2ZT9W1mMPr755hu8/vrrGvdDcnR0RGRkJLy9vVFcXIyvv/4a3bp1Q3x8PDp27Ki1n4iICMyaNUvn7ZuamtapH/rZ2f/9nxfZJSKiuq7OFD+Piri4OKSlpeHbb7/VaG/evLnGnam7dOmCc+fOYeHCheViy0ydOlXjyFV+fj6cnZ1rJnEiIiICUIfO9rKzs4OpqSmy7z2MACA7OxsqlUrrOiqVqtL4sn916bMqX3/9NTw9PeHl5VVlrI+PD86ePVvhcoVCAWtra40HERER1aw6U/zI5XJ4eXkhNjZWalOr1YiNjYWfn5/Wdfz8/DTiASAmJkaKd3Nzg0ql0ojJz89HfHx8hX1WpqCgABs3bsTIkSOrFZ+UlARHzgAmIiKqU+rU115hYWEICQmBt7c3fHx8sGjRIhQWFmL48OEAgKFDh6JJkyaIiIgAAEycOBH+/v5YsGAB+vTpgw0bNuDo0aNYuXIlgLsThENDQ/Hxxx/Dw8NDOtXdyckJwcHB0nYzMjJw/fp1ZGRkoLS0FElJSQAAd3d3WFlZSXHR0dG4c+cO3njjjXK5L1q0CG5ubmjdujWKiorw9ddfY//+/di7d28NjRYRERHpo04VPwMGDMCVK1cQHh6OrKwseHp6Yvfu3dKE5YyMDJiY/HewqkuXLli3bh2mTZuGDz/8EB4eHti2bZt0jR8AmDx5MgoLCzFmzBjk5uaia9eu2L17t3SNHwAIDw/HmjVrpOcdOnQAABw4cADdunWT2r/55hu88sorWk9lLykpwXvvvYdLly7BwsIC7dq1w759+9C9e3dDDQ8REREZQJ26zo+x0+UaBXUJL/NBZCR4nR+qw3T5HVpn5vwQERERPQwsfoiIiMiosPghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMSp0rfpYvXw5XV1colUr4+vri8OHDlcZv2rQJLVq0gFKpRNu2bbFr1y6N5UIIhIeHw9HREebm5ggICMCZM2c0YubOnYsuXbrAwsICtra2Wrcjk8nKPTZs2KAR88svv6Bjx45QKBRwd3dHVFSUzq+fiIiIaladKn6io6MRFhaGGTNmIDExEe3bt0dgYCBycnK0xh88eBCDBg3CyJEjcezYMQQHByM4OBgpKSlSzLx587BkyRJERkYiPj4elpaWCAwMRFFRkRRTUlKC/v37Y+zYsZXmt3r1amRmZkqP4OBgaVl6ejr69OmD7t27IykpCaGhoRg1ahT27NnzYINCREREBiUTQojaTqKMr68vOnXqhGXLlgEA1Go1nJ2dMWHCBEyZMqVc/IABA1BYWIgdO3ZIbZ07d4anpyciIyMhhICTkxPee+89TJo0CQCQl5cHBwcHREVFYeDAgRr9RUVFITQ0FLm5ueW2JZPJsHXrVo2C514ffPABdu7cqVF4DRw4ELm5udi9e3e1Xn9+fj5sbGyQl5cHa2vraq1TF3h7//f/o0drLw8iqmH37uzA3R2ePwCojtDld2idOfJTUlKChIQEBAQESG0mJiYICAhAXFyc1nXi4uI04gEgMDBQik9PT0dWVpZGjI2NDXx9fSvsszLjxo2DnZ0dfHx8sGrVKtxbN1aVizbFxcXIz8/XeBAREVHNMqvtBMpcvXoVpaWlcHBw0Gh3cHDAqVOntK6TlZWlNT4rK0taXtZWUUx1zZ49G8899xwsLCywd+9evP322ygoKMA777xTaS75+fm4desWzM3Ny/UZERGBWbNm6ZQHERERPZg6U/zUddOnT5f+36FDBxQWFmL+/PlS8aOPqVOnIiwsTHqen58PZ2fnB8qTiIiIKldnvvays7ODqakpsrOzNdqzs7OhUqm0rqNSqSqNL/tXlz6ry9fXF//88w+Ki4srzcXa2lrrUR8AUCgUsLa21ngQERFRzaozxY9cLoeXlxdiY2OlNrVajdjYWPj5+Wldx8/PTyMeAGJiYqR4Nzc3qFQqjZj8/HzEx8dX2Gd1JSUloUGDBlAoFNXKhYiIiOqGOvW1V1hYGEJCQuDt7Q0fHx8sWrQIhYWFGD58OABg6NChaNKkCSIiIgAAEydOhL+/PxYsWIA+ffpgw4YNOHr0KFauXAng7hlaoaGh+Pjjj+Hh4QE3NzdMnz4dTk5OGmdtZWRk4Pr168jIyEBpaSmSkpIAAO7u7rCyssJPP/2E7OxsdO7cGUqlEjExMfjkk0+kM8gA4K233sKyZcswefJkjBgxAvv378fGjRuxc+fOhzN4REREVD2ijlm6dKlo2rSpkMvlwsfHRxw6dEha5u/vL0JCQjTiN27cKJ566ikhl8tF69atxc6dOzWWq9VqMX36dOHg4CAUCoXo0aOHOH36tEZMSEiIAFDuceDAASGEED///LPw9PQUVlZWwtLSUrRv315ERkaK0tJSjX4OHDggPD09hVwuF08++aRYvXq1Tq89Ly9PABB5eXk6rVfbvLz+exDRY+zenb1sh+cPAKojdPkdWqeu82PseJ0fIqrTeJ0fqsMeyev8EBERET0MDzTn5/bt28jKysLNmzdhb2+Phg0bGiovIiIiohqh85Gff//9FytWrIC/vz+sra3h6uqKli1bwt7eHi4uLhg9ejSOHDlSE7kSERERPTCdip8vvvgCrq6uWL16NQICArBt2zYkJSXhr7/+QlxcHGbMmIE7d+7ghRdeQM+ePcvdPZ2IiIiotun0tdeRI0fw22+/oXXr1lqX+/j4YMSIEYiMjMTq1avx+++/w8PDwyCJEhERERmCTsXP+vXrqxWnUCjw1ltv6ZUQERERUU3i2V5ERERkVPQufm7duoWbN29Kzy9cuIBFixZh7969BkmMiIiIqCboXfwEBQVh7dq1AIDc3Fz4+vpiwYIFCAoKwooVKwyWIBEREZEh6V38JCYm4plnngEAbN68GQ4ODrhw4QLWrl2LJUuWGCxBIiIiIkPSu/i5efMm6tevDwDYu3cvXnnlFZiYmKBz5864cOGCwRIkIiIiMiS9ix93d3ds27YNFy9exJ49e/DCCy8AAHJych6p+1IRERGRcdG7+AkPD8ekSZPg6uoKHx8f+Pn5Abh7FKhDhw4GS5CIiIjIkPS+t5evry8yMjKQmZkJT09Pqb1Hjx545ZVXDJEbERERkcHpXfy4uLigYcOGaN++PTw9PaWHqakpPvnkE6xZs8aQeRIREREZhN7FT3p6Oo4dO4akpCQcO3YMGzduxOXLlwGAc36IiIioznqgIz8uLi4IDg6W2uLi4hASEoLZs2cbIjciIiIigzPo7S38/PywePFifP7554bsloiIiMhg9C5+SkpKtLZ7eHjg5MmTeidEREREVJP0/trLysoKrVq1QocOHeDp6YkOHTrAyckJS5cuRUBAgCFzJCIiIjIYvYuf/fv3Izk5GcnJyfj+++8xdepUFBUVAQB69uyJ8PBwtG3bFm3btkWLFi0MljARERHRg9C7+OnatSu6du0qPVer1Th9+jSSkpKQlJSEw4cP46uvvkJOTg5KS0sNkiwRERHRg9K7+LmfiYkJWrZsiZYtW2LQoEFSe3Z2tqE2QURERPTADHq2lzYODg41vQkiIiKiaqvx4oeIiIioLmHxQ0REREaFxQ8REREZFZ2Ln/DwcCQkJNRELkREREQ1Tufi559//kGvXr3wxBNPYOzYsfj5558rvNozERERUV2jc/GzatUqZGVlYf369ahfvz5CQ0NhZ2eHfv36Ye3atbh+/XpN5ElERERkEHrN+TExMcEzzzyDefPm4fTp04iPj4evry++/PJLODk54dlnn8Xnn3+OS5cuGTpfIiIiogdikAnPLVu2xOTJk/Hnn3/i4sWLCAkJwe+//47169cbonsiIiIigzHYFZ7L2NvbY+TIkRg5cqShuyYiIiJ6YDzVnYiIiIwKix8iIiIyKix+iIiIyKjUueJn+fLlcHV1hVKphK+vLw4fPlxp/KZNm9CiRQsolUq0bdsWu3bt0lguhEB4eDgcHR1hbm6OgIAAnDlzRiNm7ty56NKlCywsLGBra1tuG8nJyRg0aBCcnZ1hbm6Oli1bYvHixRoxv/zyC2QyWblHVlaWfgNBRERENcLgxc/FixcxYsQIvdaNjo5GWFgYZsyYgcTERLRv3x6BgYHIycnRGn/w4EEMGjQII0eOxLFjxxAcHIzg4GCkpKRIMfPmzcOSJUsQGRmJ+Ph4WFpaIjAwEEVFRVJMSUkJ+vfvj7Fjx2rdTkJCAho3bozvvvsOJ0+exEcffYSpU6di2bJl5WJPnz6NzMxM6dG4cWO9xoKIiIhqhkwIIQzZYXJyMjp27IjS0lKd1/X19UWnTp2kokKtVsPZ2RkTJkzAlClTysUPGDAAhYWF2LFjh9TWuXNneHp6IjIyEkIIODk54b333sOkSZMAAHl5eXBwcEBUVBQGDhyo0V9UVBRCQ0ORm5tbZa7jxo1DWloa9u/fD+DukZ/u3bvjxo0bWo8eVUd+fj5sbGyQl5cHa2trvfqoDd7e//3/6NHay4OIati9Oztwd4fnDwCqI3T5Harzqe7bt2+vdPnff/+ta5cA7h59SUhIwNSpU6U2ExMTBAQEIC4uTus6cXFxCAsL02gLDAzEtm3bAADp6enIyspCQECAtNzGxga+vr6Ii4srV/zoIi8vDw0bNizX7unpieLiYrRp0wYzZ87E008/XWEfxcXFKC4ulp7n5+frnQ8RERFVj87FT3BwMGQyGSo7YCSTyXRO5OrVqygtLYWDg4NGu4ODA06dOqV1naysLK3xZfNsyv6tLEYfBw8eRHR0NHbu3Cm1OTo6IjIyEt7e3iguLsbXX3+Nbt26IT4+Hh07dtTaT0REBGbNmqV3HkRERKQ7nef8ODo6YsuWLVCr1VofiYmJNZFnnZGSkoKgoCDMmDEDL7zwgtTevHlzvPnmm/Dy8kKXLl2watUqdOnSBQsXLqywr6lTpyIvL096XLx48WG8BCIiIqOmc/Hj5eWFhISECpdXdVSoInZ2djA1NUV2drZGe3Z2NlQqldZ1VCpVpfFl/+rSZ2VSU1PRo0cPjBkzBtOmTasy3sfHB2fPnq1wuUKhgLW1tcaDiIiIapbOxc/777+PLl26VLjc3d0dBw4c0DkRuVwOLy8vxMbGSm1qtRqxsbHw8/PTuo6fn59GPADExMRI8W5ublCpVBox+fn5iI+Pr7DPipw8eRLdu3dHSEgI5s6dW611kpKS4OjoqNN2iIiIqGbpPOfnmWeeqXS5paUl/P399UomLCwMISEh8Pb2ho+PDxYtWoTCwkIMHz4cADB06FA0adIEERERAICJEyfC398fCxYsQJ8+fbBhwwYcPXoUK1euBHD3KFRoaCg+/vhjeHh4wM3NDdOnT4eTkxOCg4Ol7WZkZOD69evIyMhAaWkpkpKSANwt5KysrJCSkoLnnnsOgYGBCAsLk+YLmZqawt7eHgCwaNEiuLm5oXXr1igqKsLXX3+N/fv3Y+/evXqNBREREdUMnYqfjIwMNG3atNrxly5dQpMmTaodP2DAAFy5cgXh4eHIysqCp6cndu/eLU1YzsjIgInJfwerunTpgnXr1mHatGn48MMP4eHhgW3btqFNmzZSzOTJk1FYWIgxY8YgNzcXXbt2xe7du6FUKqWY8PBwrFmzRnreoUMHAMCBAwfQrVs3bN68GVeuXMF3332H7777TopzcXHB+fPnAdw9W+29997DpUuXYGFhgXbt2mHfvn3o3r17tV8/ERER1TydrvPj4OCA4OBgjBo1Cp06ddIak5eXh40bN2Lx4sUYM2YM3nnnHYMl+7jjdX6IqE7jdX6oDqux6/ykpqZi7ty5eP7556FUKuHl5QUnJycolUrcuHEDqampOHnyJDp27Ih58+ahd+/eD/RCiIiIiAxNpwnPjRo1whdffIHMzEwsW7YMHh4euHr1qnSvrMGDByMhIQFxcXEsfIiIiKhO0nnCMwCYm5vj1VdfxauvvmrofIiIiIhqVJ27qzsRERFRTWLxQ0REREaFxQ8REREZFRY/REREZFRY/BAREZFR0bv4uXXrFm7evCk9v3DhAhYtWsTbORAREVGdpnfxExQUhLVr1wIAcnNz4evriwULFiAoKAgrVqwwWIJEREREhqR38ZOYmCjd5HTz5s1wcHDAhQsXsHbtWixZssRgCRIREREZkt7Fz82bN1G/fn0AwN69e/HKK6/AxMQEnTt3xoULFwyWIBEREZEh6V38uLu7Y9u2bbh48SL27NmDF154AQCQk5PzSN2Uk4iIiIyL3sVPeHg4Jk2aBFdXV/j4+MDPzw/A3aNAHTp0MFiCRERERIak1729AODVV19F165dkZmZCU9PT6m9R48e6Nu3ryFyIyIiIjK4B7rOz5kzZ7BgwQI8/fTTuHTpEgDg9OnTuHr1qkGSIyIiIjI0vYufH374AYGBgTA3N0diYiKKi4sBAHl5efjkk08MliARERGRIeld/Hz88ceIjIzEV199hXr16kntTz/9NBITEw2SHBEREZGh6V38nD59Gs8++2y5dhsbG+Tm5j5ITkREREQ1Ru/iR6VS4ezZs+Xa//jjDzz55JMPlBQRERFRTdG7+Bk9ejQmTpyI+Ph4yGQyXL58Gd9//z0mTZqEsWPHGjJHIiIiIoPR+1T3KVOmQK1Wo0ePHrh58yaeffZZKBQKTJo0CRMmTDBkjkREREQGo3fxI5PJ8NFHH+H999/H2bNnUVBQgFatWsHKysqQ+REREREZlN7FT0ZGBpydnSGXy9GqVatyy5o2bfrAyREREREZmt5zftzc3HDlypVy7deuXYObm9sDJUVERERUU/QufoQQkMlk5doLCgqgVCofKCkiIiKimqLz115hYWEA7s75mT59OiwsLKRlpaWliI+P17jXFxEREVFdonPxc+zYMQB3j/ycOHECcrlcWiaXy9G+fXtMmjTJcBkSERERGZDOxc+BAwcAAMOHD8fixYthbW1t8KSIiIiIaoreZ3utXr0aAJCamoqMjAyUlJRoLH/55ZcfLDMiIiKiGqB38ZOeno7g4GCcOHECMpkMQggAkCZBl5aWGiZDIiIiIgPS+2yvd955B25ubsjJyYGFhQVOnjyJ3377Dd7e3vjll18MmCIRERGR4eh95CcuLg779++HnZ0dTExMYGJigq5duyIiIgLvvPOONDGaiIiIqC7R+8hPaWkp6tevDwCws7PD5cuXAQAuLi44ffq0YbIjIiIiMjC9j/y0adMGycnJcHNzg6+vL+bNmwe5XI6VK1fiySefNGSORERERAajd/Ezbdo0FBYWAgBmz56NF198Ec888wwaNWqE6OhogyVIREREZEh6f+0VGBiIV155BQDg7u6OU6dO4erVq8jJycFzzz2nd0LLly+Hq6srlEolfH19cfjw4UrjN23ahBYtWkCpVKJt27bYtWuXxnIhBMLDw+Ho6Ahzc3MEBATgzJkzGjFz585Fly5dYGFhAVtbW63bycjIQJ8+fWBhYYHGjRvj/fffx507dzRifvnlF3Ts2BEKhQLu7u6IiorS+fUTERFRzdK7+NGmYcOGkMlkuHXrll7rR0dHIywsDDNmzEBiYiLat2+PwMBA5OTkaI0/ePAgBg0ahJEjR+LYsWMIDg5GcHAwUlJSpJh58+ZhyZIliIyMRHx8PCwtLREYGIiioiIppqSkBP3798fYsWO1bqe0tBR9+vRBSUkJDh48iDVr1iAqKgrh4eFSTHp6Ovr06YPu3bsjKSkJoaGhGDVqFPbs2aPXWBAREVENEQZUVFQkPv/8c+Hg4KDX+j4+PmLcuHHS89LSUuHk5CQiIiK0xr/22muiT58+Gm2+vr7izTffFEIIoVarhUqlEvPnz5eW5+bmCoVCIdavX1+uv9WrVwsbG5ty7bt27RImJiYiKytLaluxYoWwtrYWxcXFQgghJk+eLFq3bq2x3oABA0RgYGAVr/o/eXl5AoDIy8ur9jp1gZfXfw8ieozdu7OX7fD8AUB1hC6/Q3U+8lNcXIypU6fC29sbXbp0wbZt2wDcveKzm5sbFi1ahHfffVfnIqykpAQJCQkICAiQ2kxMTBAQEIC4uDit68TFxWnEA3e/jiuLT09PR1ZWlkaMjY0NfH19K+yzou20bdsWDg4OGtvJz8/HyZMnq5WLNsXFxcjPz9d4EBERUc3SufgJDw/HihUr4OrqivPnz6N///4YM2YMFi5ciC+++ALnz5/HBx98oHMiV69eRWlpqUaBAQAODg7IysrSuk5WVlal8WX/6tKnLtu5dxsVxeTn51f4NWBERARsbGykh7Ozc7VzIiIiIv3oXPxs2rQJa9euxebNm7F3716Ulpbizp07SE5OxsCBA2FqaloTeT6Wpk6diry8POlx8eLF2k6JiIjosadz8fPPP//Ay8sLwN1r/SgUCrz77rvSPb30ZWdnB1NTU2RnZ2u0Z2dnQ6VSaV1HpVJVGl/2ry596rKde7dRUYy1tTXMzc219qtQKGBtba3xICIiopqlc/FTWloKuVwuPTczM4OVldUDJyKXy+Hl5YXY2FipTa1WIzY2Fn5+flrX8fPz04gHgJiYGCnezc0NKpVKIyY/Px/x8fEV9lnRdk6cOKFx1llMTAysra3RqlWrauVCREREdYPOFzkUQmDYsGFQKBQAgKKiIrz11luwtLTUiNuyZYvOyYSFhSEkJATe3t7w8fHBokWLUFhYiOHDhwMAhg4diiZNmiAiIgIAMHHiRPj7+2PBggXo06cPNmzYgKNHj2LlypUA7t5hPjQ0FB9//DE8PDzg5uaG6dOnw8nJCcHBwdJ2MzIycP36dWRkZKC0tBRJSUkA7l6/yMrKCi+88AJatWqFIUOGYN68ecjKysK0adMwbtw4aRzeeustLFu2DJMnT8aIESOwf/9+bNy4ETt37tR5HIiIiKjm6Fz8hISEaDx/4403DJbMgAEDcOXKFYSHhyMrKwuenp7YvXu3NJE4IyMDJib/Hazq0qUL1q1bh2nTpuHDDz+Eh4cHtm3bhjZt2kgxkydPRmFhIcaMGYPc3Fx07doVu3fvhlKplGLCw8OxZs0a6XmHDh0AAAcOHEC3bt1gamqKHTt2YOzYsfDz84OlpSVCQkIwe/ZsaR03Nzfs3LkT7777LhYvXownnngCX3/9NQIDAw02PkRERPTgZEIIUdtJ0F35+fmwsbFBXl7eIzX/x9v7v/8fPVp7eRBRDbt3Zwfu7vD8AUB1hC6/Qw16hWciIiKiuo7FDxERERkVFj9ERERkVFj8EBERkVFh8UNERERGRedT3e8VGxuL2NhY5OTkQK1WayxbtWrVAyVGREREVBP0Ln5mzZqF2bNnw9vbG46Ojg98ewsiIiKih0Hv4icyMhJRUVEYMmSIIfMhIiIiqlF6z/kpKSlBly5dDJkLERERUY3Tu/gZNWoU1q1bZ8hciIiIiGqc3l97FRUVYeXKldi3bx/atWuHevXqaSz/4osvHjg5IiIiIkPTu/g5fvw4PD09AQApKSkayzj5mYiIiOoqvYufAwcOGDIPIiIiooeCFzkkIiIio/JAFznMzc3FN998g7S0NABAq1atMHLkSNjY2BgkOSIiIiJD0/vIz9GjR9GsWTMsXLgQ169fx/Xr17Fw4UI0a9YMiYmJhsyRiIiIyGD0PvLz7rvv4uWXX8ZXX30FM7O73dy5cwejRo1CaGgofvvtN4MlSURERGQoehc/R48e1Sh8AMDMzAyTJ0+Gt7e3QZIjIiIiMjS9v/aytrZGRkZGufaLFy+ifv36D5QUERERUU3Ru/gZMGAARo4ciejoaFy8eBEXL17Ehg0bMGrUKAwaNMiQORIREREZjN5fe33++eeQyWQYOnQo7ty5AyEE5HI5xo4di08//dSQORIREREZjN7Fj1wux+LFixEREYFz584BAJo1awYLCwuDJUdERERkaDoVP2FhYZgzZw4sLS0RFhZWaSzv7UVERER1kU7Fz7Fjx3D79m3p/xXhvb2IiIiortKp+Ln3fl68txcRERE9ivQ+2ysjIwNCiAqXEREREdVFehc/bm5uuHLlSrn2a9euwc3N7YGSIiIiIqopehc/Qgitc3sKCgqgVCofKCkiIiKimqLzqe5lZ3nJZDJMnz5d49T20tJSxMfHw9PT02AJEhERERmSzsVP2VleQgicOHECcrlcWiaXy9G+fXtMmjTJcBkSERERGZDOxU/ZWV7Dhw/H4sWLYW1tbfCkiIiIiGqK3ld4Xr16tSHzICIiInooeIVnIiIiMiq8wjMREREZFV7hmYiIiIyK3tf5uXXrFm7evCk9v3DhAhYtWoS9e/caJDEiIiKimqB38RMUFIS1a9cCAHJzc+Hj44MFCxYgKCgIK1aseKCkli9fDldXVyiVSvj6+uLw4cOVxm/atAktWrSAUqlE27ZtsWvXLo3lQgiEh4fD0dER5ubmCAgIwJkzZzRirl+/jsGDB8Pa2hq2trYYOXIkCgoKpOUzZ86ETCYr97C0tJRioqKiyi03xgs+entrPoiIiOoSvYufxMREPPPMMwCAzZs3Q6VS4cKFC1i7di2WLFmid0LR0dEICwvDjBkzkJiYiPbt2yMwMBA5OTla4w8ePIhBgwZh5MiROHbsGIKDgxEcHIyUlBQpZt68eViyZAkiIyMRHx8PS0tLBAYGoqioSIoZPHgwTp48iZiYGOzYsQO//fYbxowZIy2fNGkSMjMzNR6tWrVC//79NfKxtrbWiLlw4YLeY0FEREQ1QOjJ3NxcXLhwQQghRP/+/cXMmTOFEEJkZGQIc3NzfbsVPj4+Yty4cdLz0tJS4eTkJCIiIrTGv/baa6JPnz4abb6+vuLNN98UQgihVquFSqUS8+fPl5bn5uYKhUIh1q9fL4QQIjU1VQAQR44ckWJ+/vlnIZPJxKVLl7RuNykpSQAQv/32m9S2evVqYWNjo9sLvkdeXp4AIPLy8vTuozZ4ef33uP95WRsRPQa07dzc2amO0OV3qN5Hftzd3bFt2zZcvHgRe/bswQsvvAAAyMnJ0fvChyUlJUhISEBAQIDUZmJigoCAAMTFxWldJy4uTiMeAAIDA6X49PR0ZGVlacTY2NjA19dXiomLi4OtrS287/mOJiAgACYmJoiPj9e63a+//hpPPfWUdPSrTEFBAVxcXODs7IygoCCcPHmywtdbXFyM/Px8jQcRERHVLL2Ln/DwcEyaNAmurq7w9fWFn58fAGDv3r3o0KGDXn1evXoVpaWlcHBw0Gh3cHBAVlaW1nWysrIqjS/7t6qYxo0bayw3MzNDw4YNtW63qKgI33//PUaOHKnR3rx5c6xatQo//vgjvvvuO6jVanTp0gX//POP1twjIiJgY2MjPZydnbXGERERkeHoXfy8+uqryMjIwNGjR7F7926pvUePHli4cKFBkqurtm7din///RchISEa7X5+fhg6dCg8PT3h7++PLVu2wN7eHl9++aXWfqZOnYq8vDzpcfHixYeRPhERkVHT+/YWAKBSqaBSqTTafHx89O7Pzs4OpqamyM7O1mjPzs4ut517c6gsvuzf7OxsODo6asSU3X1epVKVm1B9584dXL9+Xet2v/76a7z44ovljibdr169eujQoQPOnj2rdblCoYBCoai0DyIiIjKsByp+YmNjERsbi5ycHKjVao1lq1at0rk/uVwOLy8vxMbGIjg4GACgVqsRGxuL8ePHa13Hz88PsbGxCA0NldpiYmKkr+Hc3NygUqkQGxsrFTv5+fmIj4/H2LFjpT5yc3ORkJAALy8vAMD+/fuhVqvh6+ursb309HQcOHAA27dvr/L1lJaW4sSJE+jdu7cuw0BEREQ1SO/iZ9asWZg9eza8vb3h6OhosFtahIWFISQkBN7e3vDx8cGiRYtQWFiI4cOHAwCGDh2KJk2aICIiAgAwceJE+Pv7Y8GCBejTpw82bNiAo0ePYuXKlQDu3mojNDQUH3/8MTw8PODm5obp06fDyclJKrBatmyJnj17YvTo0YiMjMTt27cxfvx4DBw4EE5OThr5rVq1Co6OjujVq1e53GfPno3OnTvD3d0dubm5mD9/Pi5cuIBRo0YZZGyIiIjoweld/ERGRiIqKgpDhgwxZD4YMGAArly5gvDwcGRlZcHT0xO7d++WvmLKyMiAicl/U5W6dOmCdevWYdq0afjwww/h4eGBbdu2oU2bNlLM5MmTUVhYiDFjxiA3Nxddu3bF7t27NS5A+P3332P8+PHo0aMHTExM0K9fv3LXK1Kr1YiKisKwYcNgampaLvcbN25g9OjRyMrKQoMGDeDl5YWDBw+iVatWBh0jIiIi0p9MCCH0WbFRo0Y4fPgwmjVrZuicjFZ+fj5sbGyQl5en9+UCasO9V3E+erT8VZ2PHn24+RBRDdG2c9//A4ColujyO1Tvs71GjRqFdevW6bs6ERERUa3Q+2uvoqIirFy5Evv27UO7du1Qr149jeVffPHFAydHREREZGh6Fz/Hjx+Xzp669z5aAAw2+ZmIiIjI0PQufg4cOGDIPIiIiIgeige6zg8ApKamIiMjAyUlJVKbTCbDSy+99KBdExERERmc3sXP33//jb59++LEiROQyWQoO2ms7Cuv0tJSw2RIREREZEB6n+01ceJEuLm5IScnBxYWFjh58iR+++03eHt745dffjFgikRERESGo/eRn7i4OOzfvx92dnYwMTGBiYkJunbtioiICLzzzjs4duyYIfMkIiIiMgi9j/yUlpaifv36AO7ekPTy5csAABcXF5w+fdow2REREREZmN5Hftq0aYPk5GS4ubnB19cX8+bNg1wux8qVK/Hkk08aMkciIiIig9G7+Jk2bRoKCwsB3L2h54svvohnnnkGjRo1QnR0tMESJCIiIjIkvYufwMBA6f/u7u44deoUrl+/jgYNGvAih0RERFRn6TXn5/bt2+jRowfOnDmj0d6wYUMWPkRERFSn6VX81KtXD8ePHzd0LkREREQ1Tu+zvd544w188803hsyFiIiIqMbpPefnzp07WLVqFfbt2wcvLy9YWlpqLOdd3YmIiKgu0rv4SUlJQceOHQEAf/31l8YyzvshIiKiukrn4mf27NmYNGkS7+pOREREjySd5/zMmjULBQUFNZELERERUY3Tufgpu3s7ERER0aNIr7O9OKeHiIiIHlV6TXh+6qmnqiyArl+/rldCRERERDVJr+Jn1qxZsLGxMXQuRERERDVOr+Jn4MCBaNy4saFzISIiIqpxOs/54XwfIiIiepTxbC8iIiIyKjp/7aVWq2siDyIiIqKHQu8bmxIRERE9ilj8EBERkVFh8UNERERGhcUPERERGRUWP0RERGRUWPwQERGRUWHxQ0REREaFxQ8REREZFRY/REREZFRY/BAREZFRqZPFz/Lly+Hq6gqlUglfX18cPny40vhNmzahRYsWUCqVaNu2LXbt2qWxXAiB8PBwODo6wtzcHAEBAThz5oxGzPXr1zF48GBYW1vD1tYWI0eOREFBgbT8/PnzkMlk5R6HDh3SKRciIiKqXXWu+ImOjkZYWBhmzJiBxMREtG/fHoGBgcjJydEaf/DgQQwaNAgjR47EsWPHEBwcjODgYKSkpEgx8+bNw5IlSxAZGYn4+HhYWloiMDAQRUVFUszgwYNx8uRJxMTEYMeOHfjtt98wZsyYctvbt28fMjMzpYeXl5dOuRAREVHtkok6dpt2X19fdOrUCcuWLQNw90aqzs7OmDBhAqZMmVIufsCAASgsLMSOHTukts6dO8PT0xORkZEQQsDJyQnvvfceJk2aBADIy8uDg4MDoqKiMHDgQKSlpaFVq1Y4cuQIvL29AQC7d+9G79698c8//8DJyQnnz5+Hm5sbjh07Bk9PT625V5XL/YqLi1FcXCw9z8/Ph7OzM/Ly8mBtba374NWS/x8yAMDRo5rPy9qI6DGgbee+/wcAUS3Jz8+HjY1NtX6H1qkjPyUlJUhISEBAQIDUZmJigoCAAMTFxWldJy4uTiMeAAIDA6X49PR0ZGVlacTY2NjA19dXiomLi4Otra1U+ABAQEAATExMEB8fr9H3yy+/jMaNG6Nr167Yvn27TrncLyIiAjY2NtLD2dlZaxwREREZTp0qfq5evYrS0lI4ODhotDs4OCArK0vrOllZWZXGl/1bVUzjxo01lpuZmaFhw4ZSjJWVFRYsWIBNmzZh586d6Nq1K4KDgzUKoKpyud/UqVORl5cnPS5evKg1joiIiAzHrLYTeFTY2dkhLCxMet6pUydcvnwZ8+fPx8svv6xXnwqFAgqFwlApEhERUTXUqSM/dnZ2MDU1RXZ2tkZ7dnY2VCqV1nVUKlWl8WX/VhVz/4TqO3fu4Pr16xVuF7g7P+ns2bPVzoWIiIhqX50qfuRyOby8vBAbGyu1qdVqxMbGws/PT+s6fn5+GvEAEBMTI8W7ublBpVJpxOTn5yM+Pl6K8fPzQ25uLhISEqSY/fv3Q61Ww9fXt8J8k5KS4OjoWO1ciIiIqPbVua+9wsLCEBISAm9vb/j4+GDRokUoLCzE8OHDAQBDhw5FkyZNEBERAQCYOHEi/P39sWDBAvTp0wcbNmzA0aNHsXLlSgCATCZDaGgoPv74Y3h4eMDNzQ3Tp0+Hk5MTgoODAQAtW7ZEz549MXr0aERGRuL27dsYP348Bg4cCCcnJwDAmjVrIJfL0aFDBwDAli1bsGrVKnz99ddS7lXlQkRERLWvzhU/AwYMwJUrVxAeHo6srCx4enpi9+7d0kTijIwMmJj8d8CqS5cuWLduHaZNm4YPP/wQHh4e2LZtG9q0aSPFTJ48GYWFhRgzZgxyc3PRtWtX7N69G0qlUor5/vvvMX78ePTo0QMmJibo168flixZopHbnDlzcOHCBZiZmaFFixaIjo7Gq6++qlMuREREVLvq3HV+jJku1yioS3idHyIjwev8UB32yF7nh4iIiKimsfghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqNS521vQ44lXfSYiorqCR36IiIjIqLD4ISIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqNiVtsJEBFRHeTtrfn86NHayYOoBvDIDxERERkVFj9ERERkVFj8EBERkVFh8UNERERGhcUPERERGRUWP0RERGRUWPwQERGRUWHxQ0REREaFxQ8REREZFRY/REREZFRY/BAREZFRqZPFz/Lly+Hq6gqlUglfX18cPny40vhNmzahRYsWUCqVaNu2LXbt2qWxXAiB8PBwODo6wtzcHAEBAThz5oxGzPXr1zF48GBYW1vD1tYWI0eOREFBgbT8l19+QVBQEBwdHWFpaQlPT098//33Gn1ERUVBJpNpPJRK5QOOBhERERlSnSt+oqOjERYWhhkzZiAxMRHt27dHYGAgcnJytMYfPHgQgwYNwsiRI3Hs2DEEBwcjODgYKSkpUsy8efOwZMkSREZGIj4+HpaWlggMDERRUZEUM3jwYJw8eRIxMTHYsWMHfvvtN4wZM0ZjO+3atcMPP/yA48ePY/jw4Rg6dCh27NihkY+1tTUyMzOlx4ULFww8QkRERPQgZEIIUdtJ3MvX1xedOnXCsmXLAABqtRrOzs6YMGECpkyZUi5+wIABKCws1ChCOnfuDE9PT0RGRkIIAScnJ7z33nuYNGkSACAvLw8ODg6IiorCwIEDkZaWhlatWuHIkSPw/v87Ge/evRu9e/fGP//8AycnJ6259unTBw4ODli1ahWAu0d+QkNDkZubq9drz8/Ph42NDfLy8mBtba1XH7Xh3ps/Hz2q/WbQvEE00SOmujtydX4AED0EuvwOrVNHfkpKSpCQkICAgACpzcTEBAEBAYiLi9O6TlxcnEY8AAQGBkrx6enpyMrK0oixsbGBr6+vFBMXFwdbW1up8AGAgIAAmJiYID4+vsJ88/Ly0LBhQ422goICuLi4wNnZGUFBQTh58mSF6xcXFyM/P1/jYUy8vf97EBERPSx1qvi5evUqSktL4eDgoNHu4OCArKwsretkZWVVGl/2b1UxjRs31lhuZmaGhg0bVrjdjRs34siRIxg+fLjU1rx5c6xatQo//vgjvvvuO6jVanTp0gX//POP1j4iIiJgY2MjPZydnbXGERERkeHUqeLnUXHgwAEMHz4cX331FVq3bi21+/n5YejQofD09IS/vz+2bNkCe3t7fPnll1r7mTp1KvLy8qTHxYsXH9ZLICIiMlp1qvixs7ODqakpsrOzNdqzs7OhUqm0rqNSqSqNL/u3qpj7J1TfuXMH169fL7fdX3/9FS+99BIWLlyIoUOHVvp66tWrhw4dOuDs2bNalysUClhbW2s8iIiIqGbVqeJHLpfDy8sLsbGxUptarUZsbCz8/Py0ruPn56cRDwAxMTFSvJubG1QqlUZMfn4+4uPjpRg/Pz/k5uYiISFBitm/fz/UajV8fX2ltl9++QV9+vTBZ599pnEmWEVKS0tx4sQJODo6VuPVExER0cNgVtsJ3C8sLAwhISHw9vaGj48PFi1ahMLCQmluzdChQ9GkSRNEREQAACZOnAh/f38sWLAAffr0wYYNG3D06FGsXLkSACCTyRAaGoqPP/4YHh4ecHNzw/Tp0+Hk5ITg4GAAQMuWLdGzZ0+MHj0akZGRuH37NsaPH4+BAwdKZ3odOHAAL774IiZOnIh+/fpJc4Hkcrk06Xn27Nno3Lkz3N3dkZubi/nz5+PChQsYNWrUwxxCIiIiqkSdK34GDBiAK1euIDw8HFlZWfD09MTu3bulCcsZGRkwMfnvgFWXLl2wbt06TJs2DR9++CE8PDywbds2tGnTRoqZPHkyCgsLMWbMGOTm5qJr167YvXu3xgUIv//+e4wfPx49evSAiYkJ+vXrhyVLlkjL16xZg5s3byIiIkIqvADA398fv/zyCwDgxo0bGD16NLKystCgQQN4eXnh4MGDaNWqVU0NFxEREemozl3nx5gZ23V+7l+PiOoQXueHHjGP7HV+iIiIiGoaix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqNS567wTMaL10YjIqKHgUd+iIiIyKjwyA8REZU/9Er0GOORHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lr/FCdxqs+ExGRobH4ISKimsW/YqiO4ddeREREZFR45IeIyNjwSAwZOR75ISIiIqPCIz/0yLn3j1b+wUpERLrikR8iIiIyKix+iIiIyKjway8ioscZJzcTlcPihx55/NlO9Aji5D2qRfzai4iIiIwKix8iIiIyKvzaix5L/CqMjNb9H34iKofFDxER1T7+xUIPEYsfMhqcX0mPHRYMRHph8UNGi7836JHDCp7IIFj8EN2DBRHVGfwwcgyoxrD4IaoC/9gmg7v/Q8Vf8tXHHZIMoE6e6r58+XK4urpCqVTC19cXhw8frjR+06ZNaNGiBZRKJdq2bYtdu3ZpLBdCIDw8HI6OjjA3N0dAQADOnDmjEXP9+nUMHjwY1tbWsLW1xciRI1FQUKARc/z4cTzzzDNQKpVwdnbGvHnzdM6FHn3e3pqPitrICFTng8APR83i+JIe6tyRn+joaISFhSEyMhK+vr5YtGgRAgMDcfr0aTRu3Lhc/MGDBzFo0CBERETgxRdfxLp16xAcHIzExES0adMGADBv3jwsWbIEa9asgZubG6ZPn47AwECkpqZCqVQCAAYPHozMzEzExMTg9u3bGD58OMaMGYN169YBAPLz8/HCCy8gICAAkZGROHHiBEaMGAFbW1uMGTOm2rmQ8ajOH/f8g9/A9B10vlmPl+q8x2TU6lzx88UXX2D06NEYPnw4ACAyMhI7d+7EqlWrMGXKlHLxixcvRs+ePfH+++8DAObMmYOYmBgsW7YMkZGREEJg0aJFmDZtGoKCggAAa9euhYODA7Zt24aBAwciLS0Nu3fvxpEjR+D9/zvI0qVL0bt3b3z++edwcnLC999/j5KSEqxatQpyuRytW7dGUlISvvjiC6n4qSoXoqro87tam6NHgVRLzaBWhQ/wS/5RWY+ouvhZNGp1qvgpKSlBQkICpk6dKrWZmJggICAAcXFxWteJi4tDWFiYRltgYCC2bdsGAEhPT0dWVhYCAgKk5TY2NvD19UVcXBwGDhyIuLg42NraSoUPAAQEBMDExATx8fHo27cv4uLi8Oyzz0Iul2ts57PPPsONGzfQoEGDKnO5X3FxMYqLi6XneXl5AO4eZXqUlJb+9//8fM3n1W2rjfX+l+p/X8yvVbZVJ+ZB1rvfg6xXIErva8vHqdT/2lrc97y6bVzv0VrvfoZer9Z35Lq6nv89++ivv2o+r25bTa93v5pcr6IYAyr73SmEqDK2ThU/V69eRWlpKRwcHDTaHRwccOrUKa3rZGVlaY3PysqSlpe1VRZz/1dqZmZmaNiwoUaMm5tbuT7KljVo0KDKXO4XERGBWbNmlWt3dnbWGv8osLHRr60urFcnEtWnb33XqwuDzvX4meJ6tbeePn3ru56+fevh33//hU0Vfdep4sfYTJ06VeNIkVqtxvXr19GoUSPIZDKDbis/Px/Ozs64ePEirK2tDdo3accxf/g45rWD4/7wcczLE0Lg33//hZOTU5Wxdar4sbOzg6mpKbKzszXas7OzoVKptK6jUqkqjS/7Nzs7G46Ojhoxnp6eUkxOTo5GH3fu3MH169c1+tG2nXu3UVUu91MoFFAoFBpttra2WmMNxdramjvKQ8Yxf/g45rWD4/7wccw1VXXEp0ydOtVdLpfDy8sLsbGxUptarUZsbCz8/Py0ruPn56cRDwAxMTFSvJubG1QqlUZMfn4+4uPjpRg/Pz/k5uYiISFBitm/fz/UajV8fX2lmN9++w23b9/W2E7z5s3RoEGDauVCREREdYCoYzZs2CAUCoWIiooSqampYsyYMcLW1lZkZWUJIYQYMmSImDJlihT/559/CjMzM/H555+LtLQ0MWPGDFGvXj1x4sQJKebTTz8Vtra24scffxTHjx8XQUFBws3NTdy6dUuK6dmzp+jQoYOIj48Xf/zxh/Dw8BCDBg2Slufm5goHBwcxZMgQkZKSIjZs2CAsLCzEl19+qVMutSUvL08AEHl5ebWditHgmD98HPPawXF/+DjmD6bOFT9CCLF06VLRtGlTIZfLhY+Pjzh06JC0zN/fX4SEhGjEb9y4UTz11FNCLpeL1q1bi507d2osV6vVYvr06cLBwUEoFArRo0cPcfr0aY2Ya9euiUGDBgkrKythbW0thg8fLv7991+NmOTkZNG1a1ehUChEkyZNxKefflou96pyqS1FRUVixowZoqioqLZTMRoc84ePY147OO4PH8f8wciEqMY5YURERESPiTo154eIiIioprH4ISIiIqPC4oeIiIiMCosfIiIiMiosfozE8uXL4erqCqVSCV9fXxw+fLi2U3psREREoFOnTqhfvz4aN26M4OBgnD59WiOmqKgI48aNQ6NGjWBlZYV+/fqVuyAm6efTTz+FTCZDaGio1MbxrhmXLl3CG2+8gUaNGsHc3Bxt27bF0Xtu4imEQHh4OBwdHWFubo6AgACcOXOmFjN+tJWWlmL69Olwc3ODubk5mjVrhjlz5mjcu4pjrh8WP0YgOjoaYWFhmDFjBhITE9G+fXsEBgaWu6o16efXX3/FuHHjcOjQIcTExOD27dt44YUXUFhYKMW8++67+Omnn7Bp0yb8+uuvuHz5Ml555ZVazPrxcOTIEXz55Zdo166dRjvH2/Bu3LiBp59+GvXq1cPPP/+M1NRULFiwQLrIKwDMmzcPS5YsQWRkJOLj42FpaYnAwEAUFRXVYuaPrs8++wwrVqzAsmXLkJaWhs8++wzz5s3D0qVLpRiOuZ5q9UR7eih8fHzEuHHjpOelpaXCyclJRERE1GJWj6+cnBwBQPz6669CiLsXyKxXr57YtGmTFJOWliYAiLi4uNpK85H377//Cg8PDxETEyP8/f3FxIkThRAc75rywQcfiK5du1a4XK1WC5VKJebPny+15ebmCoVCIdavX/8wUnzs9OnTR4wYMUKj7ZVXXhGDBw8WQnDMHwSP/DzmSkpKkJCQgICAAKnNxMQEAQEBiIuLq8XMHl95eXkAgIYNGwIAEhIScPv2bY33oEWLFmjatCnfgwcwbtw49OnTR2NcAY53Tdm+fTu8vb3Rv39/NG7cGB06dMBXX30lLU9PT0dWVpbGuNvY2MDX15fjrqcuXbogNjYWf/31FwAgOTkZf/zxB3r16gWAY/4g6tSNTcnwrl69itLSUjg4OGi0Ozg44NSpU7WU1eNLrVYjNDQUTz/9NNq0aQMAyMrKglwuL3fTWgcHB2RlZdVClo++DRs2IDExEUeOHCm3jONdM/7++2+sWLECYWFh+PDDD3HkyBG88847kMvlCAkJkcZW288ajrt+pkyZgvz8fLRo0QKmpqYoLS3F3LlzMXjwYADgmD8AFj9EBjRu3DikpKTgjz/+qO1UHlsXL17ExIkTERMTA6VSWdvpGA21Wg1vb2988sknAIAOHTogJSUFkZGRCAkJqeXsHk8bN27E999/j3Xr1qF169ZISkpCaGgonJycOOYPiF97Pebs7Oxgampa7kyX7OxsqFSqWsrq8TR+/Hjs2LEDBw4cwBNPPCG1q1QqlJSUIDc3VyOe74F+EhISkJOTg44dO8LMzAxmZmb49ddfsWTJEpiZmcHBwYHjXQMcHR3RqlUrjbaWLVsiIyMDAKSx5c8aw3n//fcxZcoUDBw4EG3btsWQIUPw7rvvIiIiAgDH/EGw+HnMyeVyeHl5ITY2VmpTq9WIjY2Fn59fLWb2+BBCYPz48di6dSv2798PNzc3jeVeXl6oV6+exntw+vRpZGRk8D3QQ48ePXDixAkkJSVJD29vbwwePFj6P8fb8J5++ulyl3D466+/4OLiAgBwc3ODSqXSGPf8/HzEx8dz3PV08+ZNmJho/po2NTWFWq0GwDF/ILU945pq3oYNG4RCoRBRUVEiNTVVjBkzRtja2oqsrKzaTu2xMHbsWGFjYyN++eUXkZmZKT1u3rwpxbz11luiadOmYv/+/eLo0aPCz89P+Pn51WLWj5d7z/YSguNdEw4fPizMzMzE3LlzxZkzZ8T3338vLCwsxHfffSfFfPrpp8LW1lb8+OOP4vjx4yIoKEi4ubmJW7du1WLmj66QkBDRpEkTsWPHDpGeni62bNki7OzsxOTJk6UYjrl+WPwYiaVLl4qmTZsKuVwufHx8xKFDh2o7pccGAK2P1atXSzG3bt0Sb7/9tmjQoIGwsLAQffv2FZmZmbWX9GPm/uKH410zfvrpJ9GmTRuhUChEixYtxMqVKzWWq9VqMX36dOHg4CAUCoXo0aOHOH36dC1l++jLz88XEydOFE2bNhVKpVI8+eST4qOPPhLFxcVSDMdcPzIh7rlUJBEREdFjjnN+iIiIyKiw+CEiIiKjwuKHiIiIjAqLHyIiIjIqLH6IiIjIqLD4ISIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lih+gx0K1bN4SGhtbZ/mqqT6p5tfG+Xbt2DY0bN8b58+ervc7AgQOxYMGCmkuKHissfojqkGHDhkEmk0Emk0Eul8Pd3R2zZ8/GnTt3Kl1vy5YtmDNnjsHyMHR/1XXx4kWMGDECTk5OkMvlcHFxwcSJE3Ht2rWHngtQ+wVb2efh008/1Wjftm0bZDJZLWVV8+bOnYugoCC4urpWe51p06Zh7ty5yMvLq7nE6LHB4oeojunZsycyMzNx5swZvPfee5g5cybmz5+vNbakpAQA0LBhQ9SvX99gORi6v+r4+++/4e3tjTNnzmD9+vU4e/YsIiMjERsbCz8/P1y/fv2h5lNXKJVKfPbZZ7hx40Ztp2IwZZ9bbW7evIlvvvkGI0eO1KnPNm3aoFmzZvjuu+8eND0yAix+iOoYhUIBlUoFFxcXjB07FgEBAdi+fTuAu0cixo8fj9DQUNjZ2SEwMFBqLztC0a1bN7zzzjuYPHkyGjZsCJVKhZkzZ2psQ61WY968eXB3d4dCoUDTpk0xd+5cafn9RzzKtjt+/HjY2NjAzs4O06dPhxBCitm9eze6du0KW1tbNGrUCC+++CLOnTtX7dc9btw4yOVy7N27F/7+/mjatCl69eqFffv24dKlS/joo4+kWFdXVyxatEhjfU9PT43XWVU+VY3TsGHD8Ouvv2Lx4sXS0biyr2Gqs/1u3bphwoQJCA0NRYMGDeDg4ICvvvoKhYWFGD58OOrXrw93d3f8/PPPlY5LQEAAVCoVIiIiKo2ryZzu3LlT6XuvVqsREREBNzc3mJubo3379ti8ebPGdrV9brXZtWsXFAoFOnfuLLWdO3cOMpkMO3bsQI8ePWBhYYHmzZsjPj5eY92XXnoJGzZsqHSciAAWP0R1nrm5ucZfymvWrIFcLseff/6JyMhIreusWbMGlpaWiI+Px7x58zB79mzExMRIy6dOnYpPP/0U06dPR2pqKtatWwcHB4dK81izZg3MzMxw+PBhLF68GF988QW+/vpraXlhYSHCwsJw9OhRxMbGwsTEBH379oVara7yNV6/fh179uzB22+/DXNzc41lKpUKgwcPRnR0tMYv3KpUJ5/Kxmnx4sXw8/PD6NGjkZmZiczMTDg7O1d7+2X929nZ4fDhw5gwYQLGjh2L/v37o0uXLkhMTMQLL7yAIUOG4ObNmxX2YWpqik8++QRLly7FP//8o9P2DZVTVe99REQE1q5di8jISJw8eRLvvvsu3njjDfz6668afVT1uQWA33//HV5eXhptycnJkMlk+OKLLzB9+nQkJyejadOmmDJlikacj48PDh8+jOLi4gcdJnrcCSKqM0JCQkRQUJAQQgi1Wi1iYmKEQqEQkyZNEkII4e/vLzp06FBuPX9/fzFx4kTp/127dtVY3qlTJ/HBBx8IIYTIz88XCoVCfPXVVxXmcW9/Zc9btmwp1Gq11PbBBx+Ili1bVtjHlStXBABx4sQJrX3e69ChQwKA2Lp1q9blX3zxhQAgsrOzhRBCuLi4iIULF2rEtG/fXsyYMUOnfCobp8pyrs727+//zp07wtLSUgwZMkRqy8zMFABEXFyc1pzv/Tx07txZjBgxQgghxNatW8X9P75rKqeq3vuioiJhYWEhDh48qLHtkSNHikGDBkl9aPvcahMUFCS9zjLh4eGiQYMGIicnR2pbsmSJaN26tUZccnKyACDOnz9frW2R8eKRH6I6ZseOHbCysoJSqUSvXr0wYMAAja8u7v+rWJt27dppPHd0dEROTg4AIC0tDcXFxejRo4dOeXXu3Fljkq2fnx/OnDmD0tJSAMCZM2cwaNAgPPnkk7C2tpYmq2ZkZFR7G6KKIztyubzafVUnn8rGyRDu7d/U1BSNGjVC27Ztpbayo23V2eZnn32GNWvWIC0t7aHnVNl7f/bsWdy8eRPPP/88rKyspMfatWs1vmaszucWAG7dugWlUqnRlpycjKCgINjb20tt6enpcHd314grO2pY2ZE0IgAwq+0EiEhT9+7dsWLFCsjlcjg5OcHMTHM3tbS0rLKPevXqaTyXyWTS1z33f61kKC+99BJcXFzw1VdfwcnJCWq1Gm3atKl0cmsZd3d3yGQypKWloW/fvuWWp6Wlwd7eHra2tgAAExOTcoXS7du3dc6nsnGqTHW2X1H/97aVFRTV2eazzz6LwMBATJ06FcOGDasTOQFAQUEBAGDnzp1o0qSJxjKFQiH9vzqfWwCws7MrN7k7OTkZU6dO1WhLSkrCs88+q9FWNin+3iKJSBse+SGqYywtLeHu7o6mTZuWK3wMwcPDA+bm5oiNjdVpvfsnlx46dAgeHh4wNTXFtWvXcPr0aUybNg09evRAy5YtdTo7qVGjRnj++efxv//9D7du3dJYlpWVhe+//17jF769vT0yMzOl5/n5+UhPT5eeP2g+ZeRyuXRk615Vbb+mfPrpp/jpp58QFxf3UHOq7L1v1aoVFAoFMjIy4O7urvHQdY4UAHTo0AGpqanS87y8PJw/fx4dOnTQiEtKSoKnp6dGW0pKCp544gnY2dnpvF0yLix+iIyMUqnEBx98gMmTJ0tfTRw6dAjffPNNpetlZGQgLCwMp0+fxvr167F06VJMnDgRANCgQQM0atQIK1euxNmzZ7F//36EhYXplNeyZctQXFyMwMBA/Pbbb7h48SJ2796N559/Hk899RTCw8Ol2Oeeew7ffvstfv/9d5w4cQIhISEwNTWVlhsiH+DuGVTx8fE4f/48rl69Kh0NqWr7NaVt27YYPHgwlixZUm5ZTeZU2Xtfv359TJo0Ce+++y7WrFmDc+fOITExEUuXLsWaNWt03lZgYCBOnjwpFavHjx+HmZmZxldzFy5cwI0bN8oVP7///jteeOEF/V8oGQ1+7UVkhKZPnw4zMzOEh4fj8uXLcHR0xFtvvVXpOkOHDsWtW7fg4+MDU1NTTJw4EWPGjAFw9yuXDRs24J133kGbNm3QvHlzLFmyBN26dat2Th4eHjhy5AhmzpyJ1157DTk5ORBC4JVXXsG3334LCwsLKXbq1KlIT0/Hiy++CBsbG8yZM0fjKIch8gGASZMmISQkBK1atcKtW7eQnp4OV1fXKrdfk2bPno3o6Ohy7TWZU2XvPQDMmTMH9vb2iIiIwN9//w1bW1t07NgRH374oc7batu2LTp27IiNGzfizTffRHJyMpo3b64xD+jYsWOwtbXVuAhiUVERtm3bht27dz/QayXjIBNVzTAkIqPXrVs3eHp6lruOTE2bMWMGvvjiC8TExGhc94Uebzt37sT777+PlJQUmJhU7wuKFStWYOvWrdi7d28NZ0ePAx75IaI6a9asWXB1dcWhQ4fg4+NT7V+E9Gjr06cPzpw5g0uXLlV73lC9evWwdOnSGs6MHhc88kNEVaqtIz9ERDWBxQ8REREZFR5DJiIiIqPC4oeIiIiMCosfIiIiMiosfoiIiMiosPghIiIio8Lih4iIiIwKix8iIiIyKix+iIiIyKiw+CEiIiKjwuKHiIiIjMr/ASfxrvpPc184AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "kets_sp, transition_rates_sp = ket.get_spontaneous_transition_rates(unit=\"1/mus\")\n", "print(f\"Number of possible spontaneous decay transitions: {len(transition_rates_sp)}\")\n", "\n", "kets_bbr, transition_rates_bbr = ket.get_black_body_transition_rates(\n", " T_in_K, \"kelvin\", unit=\"1/mus\"\n", ")\n", "print(f\"Number of considered BBR transitions: {len(transition_rates_bbr)}\")\n", "\n", "fig, ax = plt.subplots(figsize=(6, 5))\n", "\n", "n_list = np.arange(0, np.max([s.n for s in kets_bbr]) + 1)\n", "rates_summed = {}\n", "for key, kets, rates in [\n", " (\"BBR\", kets_bbr, transition_rates_bbr),\n", " (\"SP\", kets_sp, transition_rates_sp),\n", "]:\n", " rates_summed[key] = np.zeros(len(n_list))\n", " for i, s in enumerate(kets):\n", " rates_summed[key][s.n] += rates[i]\n", "\n", "ax.bar(n_list, rates_summed[\"SP\"], label=\"Spontaneous Decay\", color=\"blue\", alpha=0.8)\n", "ax.bar(n_list, rates_summed[\"BBR\"], label=\"Black Body Radiation\", color=\"red\", alpha=0.8)\n", "ax.legend()\n", "\n", "ax.set_xlabel(\"Principal Quantum Number ($n$)\")\n", "ax.set_ylabel(r\"Transition Rates (1 / $\\mu$s)\")\n", "ax.set_title(\"Spontaneous and BBR Transition Rates vs $n$\")\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "n_list = list(range(30, 90, 5))\n", "\n", "\n", "def fit_function(x: np.ndarray, a: float, b: float) -> np.typing.NDArray[Any]:\n", " return a * x + b\n", "\n", "\n", "# S states\n", "kets_s = [pi.KetAtom(\"Rb\", n=n, l=0, j=0.5, m=0.5) for n in n_list]\n", "nu_s = [ket.nu for ket in kets_s]\n", "lifetimes_s = [ket.get_lifetime(unit=\"mus\") for ket in kets_s]\n", "popt_s, _ = curve_fit(fit_function, np.log(nu_s), np.log(lifetimes_s))\n", "\n", "# circular states\n", "kets_circular = [pi.KetAtom(\"Rb\", n=n, l=n - 1, j=n - 0.5, m=n - 0.5) for n in n_list]\n", "nu_circular = [ket.nu for ket in kets_circular]\n", "lifetimes_circular = [ket.get_lifetime(unit=\"mus\") for ket in kets_circular]\n", "popt_circular, _ = curve_fit(fit_function, np.log(nu_circular), np.log(lifetimes_circular))" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAigAAAHaCAYAAAAqv7IKAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAe0VJREFUeJzt3XlYlOX+BvD7nWHfl0F2wV1R3AVXyLTUyiWtrF/uZSqYmnlKzym1OpXlqZMVrpWalWmlZnUiy1Q0TdwgDXcRZRFk3xlm5vn9MTA6bLIMzAD357q4gneb78yQc/O8zyIJIQSIiIiITIjM2AUQERERVcSAQkRERCaHAYWIiIhMDgMKERERmRwGFCIiIjI5DChERERkchhQiIiIyOQwoBAREZHJYUAhIiIik8OAQkRERCaHAYWIiIhMDgMKERERmRwGFKJWaMuWLZAkCdevX9dtO3HiBAYPHgxbW1tIkoSYmJhGeywionthQCFqYcoDwcmTJ2t9TmlpKR5//HFkZmbiv//9L7Zt2wY/P79an3/06FGsXLkS2dnZ9aiYiKgyBhSiVmjq1KkoKirShZCrV68iISEBS5YswXPPPYcpU6bA2dm51tc7evQoXnvttSoDSsXHIiKqDTNjF0BETU8ul0Mul+t+TktLAwA4OTk1+mMREdUGW1CIWqG7+4XMmDEDoaGhAIDHH38ckiThvvvu0x2blJSEWbNmwd3dHZaWlujevTs+++wz3f6VK1fiH//4BwCgXbt2kCRJr89JVX1QVq5cCUmScOnSJUyZMgWOjo5wc3PDq6++CiEEbt68ifHjx8PBwQEeHh547733Kj2He9VVlW+//RaSJOHQoUOV9m3YsAGSJOHcuXMAgLy8PCxatAj+/v6wtLREmzZt8MADD+D06dP3fH2HDx+OkJAQnD59GmPGjIG9vT28vb2xZs2ae55LRFpsQSFq5ebMmQNvb2+89dZbWLBgAQYMGAB3d3cAQGpqKgYOHAhJkjB//ny4ubnh559/xjPPPIPc3FwsWrQIEydOxKVLl7B9+3b897//hUKhAAC4ubnd87EnT56Mbt26YdWqVfjpp5/w73//Gy4uLtiwYQPuv/9+vPPOO/jyyy+xZMkSDBgwACEhIbWuqyoPP/ww7OzssHPnTl0oK7djxw50794dPXr0AADMnTsX3377LebPn4+AgABkZGTgyJEjOH/+PPr27Vvj8zp79iy8vLwwduxYzJw5ExMmTMCmTZvwwgsv4P7770dgYOA9XxuiVk8QUYuyefNmAUCcOHHinsfEx8cLIYQ4cOCAACC++eYbveOeeeYZ4enpKdLT0/W2P/nkk8LR0VEUFhYKIYRYvXq13vVqeiwhhFixYoUAIJ577jndNpVKJXx8fIQkSWLVqlW67VlZWcLa2lpMnz69znVV5amnnhJt2rQRKpVKty0lJUXIZDLx+uuv67Y5OjqK8PDwaq9TneTkZAFAuLm5iZs3b+q2x8XFCQBi69atdb4mUWvEWzxEVCUhBL777juMHTsWQgikp6frvkaNGoWcnJxa3e6oybPPPqv7Xi6Xo3///hBC4JlnntFtd3JyQpcuXXDt2jWD1DV58mSkpaXh4MGDum3ffvstNBoNJk+erPe4x48fR3Jycp2e09mzZwEAK1asgI+Pj267ubk5AMDCwqJO1yNqrRhQiKhKt2/fRnZ2NjZu3Ag3Nze9r5kzZwK407m2vtq2bav3s6OjI6ysrHS3ie7enpWVZZC6Ro8eDUdHR+zYsUO3bceOHejduzc6d+6s2/buu+/i3Llz8PX1RVBQEFauXKkLSTUpDygTJkzQ237hwgUAQJcuXao8b8uWLXqhiai1Yx8UIqqSRqMBAEyZMgXTp0+v8piePXs26DGqGt1T3YgfIYRB6rK0tMSECROwe/durF27Fqmpqfjjjz/w1ltv6R33xBNPYNiwYdi9ezf27duH1atX45133sGuXbswZsyYaq//119/wcPDA97e3nrbY2NjYWZmhoCAAL3tX331le45CyHw0UcfISAgACNGjKj2MYhaAwYUIqqSm5sb7O3toVarMXLkyBqPlSSpiaqqW13VmTx5MrZu3Yr9+/fj/PnzEELo3d4p5+npibCwMISFhSEtLQ19+/bFm2++WWNAOXv2LHr16lVp+19//YXOnTvD0tKyUi0ff/wxNm/eDGtra4SFhTGcEIG3eIioGnK5HJMmTcJ3332nG3p7t9u3b+u+t7W1BYAmmUm2LnVVZ+TIkXBxccGOHTuwY8cOBAUFoV27drr9arUaOTk5eue0adMGXl5eKCkpqfa6arUa58+frzKgxMbGVtuyUx7wJEninDFEZdiCQtRCffbZZ4iMjKy0feHChbW+xqpVq3DgwAEEBwdj9uzZCAgIQGZmJk6fPo3ffvsNmZmZAIB+/foBAP71r3/hySefhLm5OcaOHasLLoZW27qqY25ujokTJ+Lrr79GQUEB/vOf/+jtz8vLg4+PDx577DH06tULdnZ2+O2333DixIkq52Qpd/nyZRQXF1cKKEVFRbhy5UqVt6R27NiBNm3aYNGiRfDz88PZs2exf/9+tqJQq8eAQtRCrVu3rsrtM2bMqPU13N3dER0djddffx27du3C2rVr4erqiu7du+Odd97RHTdgwAC88cYbWL9+PSIjI6HRaBAfH99oAaW2ddVk8uTJ+OSTTyBJEp544gm9fTY2NggLC8O+ffuwa9cuaDQadOzYEWvXrsW8efOqvWZ5B9mKAeXcuXNQq9VVtqD83//9H4A7E9otWLCgVvUTtXSSKO95RkRERGQi2AeFiIiITA4DChEREZkcBhQiIiIyOewkW4FGo0FycjLs7e2bdG4HIiKi5k4Igby8PHh5eUEma1gbCANKBcnJyfD19TV2GURERM3WzZs39daiqg8GlArs7e0BaF9cBwcHI1dDRETUfOTm5sLX11f3WdoQDCgVlN/WcXBwYEAhIiKqB0N0kWAn2TIREREICAjAgAEDjF0KERFRq8eJ2irIzc2Fo6MjcnJy2IJCRERUB4b8DGULChEREZkc9kGpB7VajdLSUmOXQS2AXC6HmZkZh7QTEVXAgFJH+fn5SExMBO+MkaHY2NjA09MTFhYWxi6FiMhkMKDUgVqtRmJiImxsbODm5sa/eqlBhBBQKpW4ffs24uPj0alTpwZPbERE1FIwoNRBaWkphBBwc3ODtbW1scuhFsDa2hrm5uZISEiAUqmElZWVsUsiIjIJ/HOtTF2GGbPlhAyJrSZERJXxX8Yy4eHhiIuLw4kTJ4xdChERUavHWzxEREStiUYNJBwF8lMBO3fAbzAgkxu7qkoYUIxArRGIjs9EWl4x2thbIaidC+Qy3jYiIqJGFrcXiHwZyE2+s83BCxj9DhAwznh1VYG3eJpY5LkUDH3ndzy16U8s/DoGT236E0Pf+R2R51Ia7TFv376NefPmoW3btrC0tISHhwdGjRqFP/74o17XW7lyJXr37l3n87Zs2QInJ6d6PSYRETVQ3F5g5zT9cAIAuSna7XF7jVNXNRhQmlDkuRTM++I0UnKK9bbfyinGvC9ON1pImTRpEs6cOYOtW7fi0qVL2Lt3L+677z5kZGQ0yuMREZGJ0ai1LSe4M4eXRgVop/Qq2xa5VHuciWBAaSJqjcBrP8Shqundyre99kMc1BrDTgCXnZ2Nw4cP45133sHw4cPh5+eHoKAgLFu2DOPGVd+cd/DgQQQFBcHW1hZOTk4YMmQIEhISsGXLFrz22muIjY2FJEmQJAlbtmwBALz//vsIDAyEra0tfH19ERYWhvz8fN31Zs6ciZycHN15K1euBACUlJRgyZIl8Pb2hq2tLYKDg3Hw4EFdLQkJCRg7diycnZ1ha2uL7t2743//+59BXyciohYt4aiu5USjkpBxwRZXfnRHfrJl2QECyE3SHmci2AeliUTHZ1ZqObmbAJCSU4zo+EwM6uBqsMe1s7ODnZ0d9uzZg4EDB8LS0vKe56hUKkyYMAGzZ8/G9u3boVQqER0dDUmSMHnyZJw7dw6RkZH47bffAACOjo4AtMNlP/zwQ7Rr1w7Xrl1DWFgYXnrpJaxduxaDBw/GBx98gOXLl+PixYu62gBg/vz5iIuLw9dffw0vLy/s3r0bo0ePxtmzZ9GpUyeEh4dDqVQiKioKtra2iIuL051LRES1kJ8KTamErCu2yLhgC3WJtlNs9jUb2HuX6B1nKhhQykRERCAiIgJqdeM0b6XlVR9O6nNcbZmZmWHLli2YPXs21q9fj759+yI0NBRPPvkkevbsWeU5ubm5yMnJwSOPPIIOHToAALp166bbb2dnBzMzM3h4eOidt2jRIt33/v7++Pe//425c+di7dq1sLCwgKOjIyRJ0jvvxo0b2Lx5M27cuAEvLy8AwJIlSxAZGYnNmzfjrbfewo0bNzBp0iQEBgYCANq3b2+Q14aIqDVQ5xcg66fTyPyhDdRKbTAxt1VB0T0Pjv5F+gfbuRuhwqrxFk+Zxp4HpY197WYIre1xdTFp0iQkJydj7969GD16NA4ePIi+ffvqbs1U5OLighkzZmDUqFEYO3Ys1qxZg5SUe/eP+e233zBixAh4e3vD3t4eU6dORUZGBgoLC6s95+zZs1Cr1ejcubOutcfOzg6HDh3C1atXAQALFizAv//9bwwZMgQrVqzAX3/9Va/XgYioNVHn5SF93TpcHTECt7fshloph7mdCp7BWejwcBqc2hdB0qUACXDw1g45NhEMKE0kqJ0LPB2tUN1gYgmAp6N2yHFjsLKywgMPPIBXX30VR48exYwZM7BixYpqj9+8eTOOHTuGwYMHY8eOHejcuTP+/PPPao+/fv06HnnkEfTs2RPfffcdTp06hYiICACAUqms9rz8/HzI5XKcOnUKMTExuq/z589jzZo1AIBnn30W165dw9SpU3H27Fn0798fH330UT1fCSKilk2dm4vbH0fgyoiRuL3mQ6hzcmDRrh28Fk5Gh4duw6ld8V3BBED5J9PoVSY1HwoDShORyySsGBsAAJVCSvnPK8YGNNl8KAEBASgoKKjxmD59+mDZsmU4evQoevToga+++goAYGFhUelW2KlTp6DRaPDee+9h4MCB6Ny5M5KT9YeyVXVenz59oFarkZaWho4dO+p93X0ryNfXF3PnzsWuXbvw4osvYtOmTQ15+kRELY46Oxu3P/wQV+4fgfSPP4YmNxcWHTrA6z//Qfsff4DjvJWQnvwccPDUP9HBC3jic5ObB4V9UJrQ6B6eWDelL177IU6vw6yHoxVWjA3A6B6eNZxdPxkZGXj88ccxa9Ys9OzZE/b29jh58iTeffddjB8/vspz4uPjsXHjRowbNw5eXl64ePEiLl++jGnTpgHQ9i+Jj49HTEwMfHx8YG9vj44dO6K0tBQfffQRxo4diz/++APr16/Xu66/vz/y8/Oxf/9+9OrVCzY2NujcuTOefvppTJs2De+99x769OmD27dvY//+/ejZsycefvhhLFq0CGPGjEHnzp2RlZWFAwcO6PWJISJqzVRZWcjcshVZX3wBTdkfnpadOkIRFgb7Bx+EJL+rVSRgHND14WYxkywE6cnJyREARE5OTqV9RUVFIi4uThQVFTXoMVRqjTh6JV3sOZMojl5JFyq1pkHXq0lxcbFYunSp6Nu3r3B0dBQ2NjaiS5cu4pVXXhGFhYVVnnPr1i0xYcIE4enpKSwsLISfn59Yvny5UKvVumtOmjRJODk5CQBi8+bNQggh3n//feHp6Smsra3FqFGjxOeffy4AiKysLN21586dK1xdXQUAsWLFCiGEEEqlUixfvlz4+/sLc3Nz4enpKR599FHx119/CSGEmD9/vujQoYOwtLQUbm5uYurUqSI9Pb3RXrOmZqjfKyJqYdQqIa5FCfHXN9r/qlV6u0szMkTqf/4jLvTpK+K6dBVxXbqKq+PGi5zIX4Sm7N/rplbTZ2hdSUIIw0680czl5ubC0dEROTk5cHBw0NtXXFyM+Ph4tGvXDlZWhu/MSq0Tf6+IqJIapqRXtRmMjM82I2v7dogi7Sgcy4BucAsLg93990My4grpNX2G1hVv8RAREZmS8inpK0ztWZqaisxl4ciKd4ZQlgIArLp3hyI8HHbD74Mktaw13RhQiIiITEUVU9KXFsqQccEO2VdtIdQSgFJYBQbCbX44bENCWlwwKceAQkREZCrumpK+tFCGjPP2yL5qA6HRhhBrVyUUPfJg+/J8SO1DjFlpo2NAISIiMhX5qSgtkCP9vB1yrt0VTBQlcOuRBxt3JSQJQEGacetsAgwoZRp7qnsiIqKaKBOTkPHZPmT/0gYoCyY2biVQ9MiDTZuyYFLOhKakbywMKGXCw8MRHh6u64FMRETUFJQ3byJ9wwbk7PkeUKkASLBpow0mtm0qzsQtaUfzmNCU9I2FAYWIiMgIlNevI33DRuTs3QuUtd7bDh4MxcM9YXO2qqVITHNK+sbCgEJERNSESq7FI2PDeuT88COg0QAAbIcOhSI8DDZ9+mgP6taumnlQVpnclPSNhQGFiIioCZRcvYr0deuR+7//3QkmoSFwCwuDda9e+gc3pynpGwkDijFo1Cb5SydJEnbv3o0JEyY02mMcPHgQw4cPR1ZWFpycnBrtcYiITEXJ5ctIX7cOuT9HAmWTt9sNHw5FWBisA3tUf6JMDrQb1kRVmh4GlKZWw/TFjdlsd+vWLbz55pv46aefkJSUhDZt2qB3795YtGgRRowYAQBISUmBs7Nzo9VgTCtXrsSePXsQExNTp/O2bNmCRYsWITs7u1HqIqJmroY/OIsvXkT62nXI++UX3eF2I0fALSwMVgEBxqq42WBAaUrVTF+M3BTt9kZa7vr69esYMmQInJycsHr1agQGBqK0tBS//PILwsPDceHCBQCAh4dHjdcpLS2Fubm5weurK6VSCQsLC2OXQUStXTV/cBZ3WYj0yL+R9+tvus32Dz4IRdg8WHXtaoRCmyfjrSjU2lQxffEdZdsil2qPM7CwsDBIkoTo6GhMmjQJnTt3Rvfu3bF48WL8+eefuuMkScKePXsAaEONJEnYsWMHQkNDYWVlhS+//BIA8Nlnn6F79+6wtLSEp6cn5s+fr3fO3a0U2dnZkCQJBw8erLK2jIwMPPXUU/D29oaNjQ0CAwOxfft2vWPuu+8+zJ8/H4sWLYJCocCoUaOqvNbBgwcRFBQEW1tbODk5YciQIUhISMCWLVvw2muvITY2FpIkQZIkbNmyBQDw/vvvIzAwELa2tvD19UVYWBjy8/N115s5cyZycnJ0561cuRIAUFJSgiVLlsDb2xu2trYIDg7We44JCQkYO3YsnJ2dYWtri+7du+N///tfbd4uImoOyv/gvCucFGWa4+ZPxYh/cY02nEgS7MeMRru938PnwzUMJ3XEFpSmctf0xVUTQG6S9jgD3nPMzMxEZGQk3nzzTdja2lbaf69+IEuXLsV7772HPn36wMrKCuvWrcPixYuxatUqjBkzBjk5Ofjjjz/qXV9xcTH69euHl19+GQ4ODvjpp58wdepUdOjQAUFBQbrjtm7dinnz5lX7WCqVChMmTMDs2bOxfft2KJVKREdHQ5IkTJ48GefOnUNkZCR++037F035XDcymQwffvgh2rVrh2vXriEsLAwvvfQS1q5di8GDB+ODDz7A8uXLcfHiRQCAnZ0dAGD+/PmIi4vD119/DS8vL+zevRujR4/G2bNn0alTJ4SHh0OpVCIqKgq2traIi4vTnUtEzVyFPziLMsyR/rc98pPLVyMXcOggQfHfPbDs3NloZTZ3DChNJT/VsMfV0pUrVyCEQNd6JvdFixZh4sSJup///e9/48UXX8TChQt12wYMGFDv+ry9vbFkyRLdz88//zx++eUX7Ny5Uy+gdOrUCe+++26118nNzUVOTg4eeeQRdOjQAQDQrVs33X47OzuYmZlVuo21aNEi3ff+/v7497//jblz52Lt2rWwsLCAo6MjJEnSO+/GjRvYvHkzbty4AS8vLwDAkiVLEBkZic2bN+Ott97CjRs3MGnSJAQGBgIA2rdvX49Xh4hMUtkfnEXp5rj9tz0KUsqCiSTg0LYIiu75sHRQAeapABhQ6osBpanUdlpiA09fLERVt5Rqr3///rrv09LSkJycrOtUawhqtRpvvfUWdu7ciaSkJCiVSpSUlMDGxkbvuH79+tV4HRcXF8yYMQOjRo3CAw88gJEjR+KJJ56Ap6dnjef99ttvePvtt3HhwgXk5uZCpVKhuLgYhYWFlWood/bsWajVanSu8JdRSUkJXF1dAQALFizAvHnzsG/fPowcORKTJk1Cz5497/VyEFEzUHgyGukHXFCQeieYOPoVQdE9Dxb2d92mN/AfnK0N+6A0Fb/B2tE6qG5ZbAlw8Db49MWdOnWCJEm6jrB1dfdtIWtr6xqPlcm0v053h6LS0tIaz1m9ejXWrFmDl19+GQcOHEBMTAxGjRoFpVJ/eueqbk9VtHnzZhw7dgyDBw/Gjh070LlzZ70+NhVdv34djzzyCHr27InvvvsOp06dQkREBABUevy75efnQy6X49SpU4iJidF9nT9/HmvWrAEAPPvss7h27RqmTp2Ks2fPon///vjoo4/u+RyIyHQVnjiBhBkzkfDqJ9pwIgk4titAh4fS4DUwWz+cAK1ivZzGxIDSVGRy7VBiAJVDSuNNX+zi4oJRo0YhIiICBQUFlfbXZfisvb09/P39sX///ir3u7m5AdAOVy53r2G9f/zxB8aPH48pU6agV69eaN++PS5dulTrmirq06cPli1bhqNHj6JHjx746quvAAAWFhaVFoI8deoUNBoN3nvvPQwcOBCdO3dGcrJ+P6GqzuvTpw/UajXS0tLQsWNHva+7bwX5+vpi7ty52LVrF1588UVs2rSp3s+LiIxDCIGCP48jYeo0JEydhsI//wTMzODUFejw8G14BedUDiaN9Adna8OAUiYiIgIBAQEN6k9xTwHjtEOJHSrcdnDwarQhxgB0qzQHBQXhu+++w+XLl3H+/Hl8+OGHGDRoUJ2utXLlSrz33nv48MMPcfnyZZw+fVrXMmBtbY2BAwdi1apVOH/+PA4dOoRXXnmlxut16tQJv/76K44ePYrz589jzpw5SE2te7NofHw8li1bhmPHjiEhIQH79u3D5cuXdf1Q/P39ER8fj5iYGKSnp6OkpAQdO3ZEaWkpPvroI1y7dg3btm3D+vXr9a7r7++P/Px87N+/H+np6SgsLETnzp3x9NNPY9q0adi1axfi4+MRHR2Nt99+Gz/99BMAbd+WX375BfHx8Th9+jQOHDig1yeGiEybEAIFR48iYepU3JgxA4UnTgDm5nB6cjI6/hIJz7fegYWdGk35B2erI0hPTk6OACBycnIq7SsqKhJxcXGiqKioYQ+iVglxLUqIv77R/letatj1aiE5OVmEh4cLPz8/YWFhIby9vcW4cePEgQMHdMcAELt37xZCCBEfHy8AiDNnzlS61vr160WXLl2Eubm58PT0FM8//7xuX1xcnBg0aJCwtrYWvXv3Fvv27RMAdI9z4MABAUBkZWUJIYTIyMgQ48ePF3Z2dqJNmzbilVdeEdOmTRPjx4/XXTM0NFQsXLiwxud369YtMWHCBOHp6SksLCyEn5+fWL58uVCr1UIIIYqLi8WkSZOEk5OTACA2b94shBDi/fffF56ensLa2lqMGjVKfP7553r1CSHE3LlzhaurqwAgVqxYIYQQQqlUiuXLlwt/f3/d6/Doo4+Kv/76SwghxPz580WHDh2EpaWlcHNzE1OnThXp6elV1m6w3ysiajCNRiPyog6L+CefEnFduoq4Ll3F+R6BIuW114UyOVn/4L+/F+K9rkKscLjz9V437fZWqqbP0LqShGhgL8oWJjc3F46OjsjJyYGDg4PevuLiYsTHx6Ndu3awsrKq5gpEdcPfKyLjE0Kg4PBh3I6IQHHsXwAAydISTk88Addnn4G5ezX9SUx06RJjqekztK44ioeIiFotIQTyDx5E+tp1KD57FgAgWVnBefJkuDwzC+Zt2tR8gVa+Xk5jYkAhIqKW5x4tG0II5P/+O9Ij1qI4Lg4AIFlbw/mpp+A6aybMFApjVU5lGFCIiKhlqWFRVtH1EeT99hvS165DSdn0C5KNDVye/j+4zJgBs7K5jMj4GFCIiKjlqGZRVpGTgrz/PIf0pG4oSbgFAJDZ2sJ5yhS4zJgOsxa6kntzxoBSD+xXTIbE3yciA6liUVahAXJvWiH9b3soc80B3ILMzg4u06bCZdo0yO+xHhkZDwNKHcjl2vuXSqXynrOqEtVWYWEhAMDc3NzIlRA1c3ctyio0QO4Na6TH2ZUFE0BmroFLl3y4LP8v5D1HG7NSqgUGlDowMzODjY0Nbt++DXNzc93U7kT1IYRAYWEh0tLS4OTkpAvARFRP+akQGiAnwRoZcfZQ5mk/4mQWGrh2yYdzpwLILQQg5Ru5UKoNBpQ6kCQJnp6eiI+PR0JCgrHLoRbCycmp0irLRFQ3orQUOUevIP1/bVCar/1ok1toW0ycOxdAbn7XrVSukdMsMKDUkYWFBTp16lTjYnJEtWVubs6WE6IGEEolsr//HhkbNqI0MRGAGeSWarh2LYBTxwrBBJJ2NA/XyGkWGFDqQSaTccZPIiIj0iiVyNm1GxkbN6K0bJFPuasrXB8ZCOeCTyEzE9AfycM1cpobBhQiImo2NCUlyP7uO2Rs+gSqspXT5W4KKJ59Fk5PPAGZtTUQF1LNPCirGm1RVjI8BhQiIjJ5mpISZO/8BhmffAJV2YrnZm3awHX2bDg9/hhkd7dqB4wDuj7MNXKaOQYUIiIyWZqiImTv3ImMTz6F6vZtAICZhwdcZz8Lp8ceg8zSsuoTuUZOs8eAQkREJkdTWIisr3cg47PPoE5PBwCYeXlC8dxzcJw4ETILCyNXSI2NAYWIiEyGpqAAWdu3I+OzzVBnZgIAzL294TrnOThNmACJwaTVYEAhIqKmUcMKw+r8AmR9+SUyN2+GOjsbAGDu6wvF3DlwHDcOEmdabnUYUIiIqPFVs8KwethKZJ3MROaWrVDn5AAAzP3aQjF3HhzHPgLJjB9TrVWLfOf9/f3h4OAAmUwGZ2dnHDhwwNglERG1XlWsMKxWSsg8mofMzcuhKdUuG2LRrh0U8+bC4aGHGEyoZQYUADh69Cjs7OyMXQYRUetWYYVhdYmEzEt2yLxkeyeYOAGKf72rDSacWZnKtNiAQkREJqBshWFViYTMi3bIumQLjUobTCwdS6Hongd7n2JIPZwAhhO6i8ktxxsVFYWxY8fCy8sLkiRhz549lY6JiIiAv78/rKysEBwcjOjoaL39kiQhNDQUAwYMwJdfftlElRMRUUWqpGtIi7XH1R/ckRFnD41KBkunUngPyUS70bfh0LYYkgzajrNEdzG5FpSCggL06tULs2bNwsSJEyvt37FjBxYvXoz169cjODgYH3zwAUaNGoWLFy+iTZs2AIAjR47A29sbKSkpGDlyJAIDA9GzZ88qH6+kpAQlJSW6n3NzcxvniRERtSKq9HRkfLYZWV9+AVFiDwCwdFbCrXs+7LyLIUkVTuAKw1SBJIQQ9z7MOCRJwu7duzFhwgTdtuDgYAwYMAAff/wxAECj0cDX1xfPP/88li5dWuka//jHP9C9e3fMmDGjysdYuXIlXnvttUrbc3Jy4ODgYJDnQUTUWqhu30bGJ58ia8cOiOJiAICVAlB0y4SdVxXBpHyF4UVnORV9C5CbmwtHR0eDfIaa3C2emiiVSpw6dQojR47UbZPJZBg5ciSOHTsGQNsCk5eXBwDIz8/H77//ju7du1d7zWXLliEnJ0f3dfPmzcZ9EkRELVBpahpuvfkWrox8AJlbt0IUF8OqZ0/4blgP/w2rYO9dAqlSOuEKw1Q9k7vFU5P09HSo1Wq4u+s3Bbq7u+PChQsAgNTUVDz66KMAALVajdmzZ2PAgAHVXtPS0hKW1a3lQERENSq9dQsZGzch+9tvIZRKAIB1795QhIfDduiQO6FE+pwrDFOdNKuAUhvt27dHbGysscsgImrRSpOTkb5xI3K+2wVRWgoAsO7XD27hYbAZNKhyawlXGKY6alYBRaFQQC6XIzVVv7d3amoqPDw8GnTtiIgIREREQK1WN+g6REQtmTIxCRkbNyJ7926gLJjYDBgARXg4bIKDqriNcxeuMEx10Kz6oFhYWKBfv37Yv3+/bptGo8H+/fsxaNCgBl07PDwccXFxOHHiREPLJCJqcZQ3byL5lVdwdfRoZO/cCZSWwmbgQLT9fCv8tn0O24HBNYcTojoyuRaU/Px8XLlyRfdzfHw8YmJi4OLigrZt22Lx4sWYPn06+vfvj6CgIHzwwQcoKCjAzJkzjVg1EVHLpLx+HekbNiJn716grIXZdvBgKMLDYNOvn5Gro5bM5ALKyZMnMXz4cN3PixcvBgBMnz4dW7ZsweTJk3H79m0sX74ct27dQu/evREZGVmp4ywREdWghpWFAaDkWjwyNqxHzg8/AhoNAMB22DAowubBpk8fY1VNrYhJz4PSlO7ug3Lp0iXOg0JELVc1Kwtj9DsoseyO9HXrkfu//+mCiV1oKBRh82Ddq5eRCqbmwpDzoDCgVGDIF5eIyORUsbIwABRnmyMjzg65N22Aso8Fu/vvh2LePFgH9jBCodQcGfIz1ORu8RARUSOpsLIwABRnmyH9b3vk3bQu2yJgP3IEFGFhsAoIMEqZRAADChFR61G2sjAAFGeVBZNEa91ue98iKALyYLVoMtCO4YSMiwGFiKi1yE9FUaY50v+2Q37SXS0mvsVQdM+DlZNKdxyRsTGglOFEbUTUkhX99RfS/7MT+dFuZVsEHNoWQdE9H5aOKv2DubIwmQB2kq2AnWSJqCUpionB7Yi1KDh8WLtBAhzaFmqDiUOFYMKVhamB2EmWiIhqVHj6NNI/jkDB0aPaDXI5HMeOhWJkJ1gcWVzFGVxZmEwLAwoRUQtSeOIEbkesReGff2o3mJnBcfw4KObMgUXbttptXgquLEwmjwGFiKiZE0Kg8Hg00iMiUFi+npiZGZwefRSuc56DhY+P/glcWZiaAQYUIqJmSgiBwmPHcHvtWhSdPKXdaG4Op0kToZg9G+be3tWfzJWFycQxoJThKB4iai6EECg48gfS165F0ZkzAADJ3BxOjz8G19mzYe7paeQKiRqOo3gq4CgeIjJVQggUREXh9tq1KI79CwAgWVjA6Ykn4Dr7WZhz0VQyMo7iISJqRYQQyD9wEOlr16L43DkAgGRpCecnJ8PlmWdg3qaNkSskMjwGFCIiEyWEQP7+/bi9di1K4s4DACRrazg/+SRcZ82EmZvbPa5A1HwxoBARGYNGXe0oGqHRIO/X35C+bh1KLlwAAEg2NnD5v6fgMnMmzFxdjVk5UZNgQCEiampxe6uch0Q8+Dbybloife06lFy+DACQ2djAecoUuMycATNnZyMVTNT0GFDKcBQPETWJuL3AzmkA7oxPEBog92wW0r9eAmWuOQBAZmcH56lT4DJtGoMJtUocxVMBR/EQUaPRqIEPeuhaToQGyL1hjfS/7aDMKwsmFoDLs/PgMn065I6OxqyWqM44ioeIqDlKOArkJkNogJwEa2TE2UOZp/1nWGahgUvnfLh0LoB8bB+A4YRaOQYUIqImIrKTkXPNGul/26O0QPvPr9xCA5cu+XDuXAC5eVmDdn6qEaskMg0MKEREjUwolcj+/ntkRHyM0lva/iRySzVcuhbAueNdwaScHSdcI2JAISJqJBqlEjm7diNj40aUJmv7ncitAdcuuXDuUABZxWACSbuqsN/gpi+WyMQwoBARGZimpATZ332HjE2fQJWSAgCQuyng+swzcO5lB9n3z1ZxlqT9z+hVXFWYCAwoREQGoykpQfbOb5DxySdQpWr7kZi5ucF19mw4PfE4ZFZW2gMtLaqcBwWjVwEB44xQOZHpYUApw3lQiKi+NEVFyN65ExmffArV7dsAADN3d7g+NxtOjz0GmaWl/gkB44CuD1c7kywRcR6USjgPChHVlqawEFlf70DGZ59BnZ4OADDz9ITiudlwnDQJMgsLI1dI1LQ4DwoRkRFpCgqQtX07Mj7bDHVmJgDA3MsLrnPmwOnRCZAYTIgajAGFiKiW1PkFyPryS2Ru3gx1djYAwNzHB4q5c+A4fjwkc3PjFkjUgjCgEBHdgzovryyYbIE6JwcAYO7XFoo5c+E49hEGE6JGwIBCRFQNdW4uMrdtQ+bWz6HJzQUAWPj7QzFvLhwefhiSGf8JJWos/L+LiKgCdXY2Mj/fhsxt26DJywMAWLRvD8W8eXB4aAwkOUfbEDU2BhQiaj006hqH9qqyspC5dSuytn0BTUEBAMCyU0co5s2D/ahRDCZETYgBhYhah7i91UyO9g5UHkORuXkLsr78EprCQgCAZefOUISFwf7BByDJZEYqmqj1YkAhopYvbi+wcxoA/WmfVGmpyPhnOLKuu0AUKwEAlt26QRE2D/YjRjCYEBkRAwoRtWwatbbl5K5woiqSIeOCHbKu2ECoZQCUsOoeAEV4OOyGD4ckSUYrl4i0GFDKcKp7ohYq4ajutk5pkQwZ5+2QfdUWQq0NIVYuSih65MFu6UJI7UOMWSkR3YUBpUx4eDjCw8N10/QSUQuRn4rSwruCiUYbTKxdtcHE1qMEkgSgIM24dRKRHgYUImqxSpOTkb7lN+T87H4nmChKoOiRD1v3smBSzs7dOEUSUZUYUIioxVEmJiFj0yZk79oFlJYCkGDjVgJFjzzYtFHqBxNI2tE8foONVC0RVYUBhYhaDOXNm8jYuBHZu/cAKhUAwCY4GIpHesP275VVnFGWVEav0psPhYiMjwGFiJo9ZUIC0jdsRM733wNlHd1tBw+CIiwMNv37aw/q3r6aeVBWAQHjjFA1EdWEAYWImq2S+HhkrN+AnB9/vBNMhg7VBpO+ffQPDhgHdH24xplkich0MKAQUbNTcvUq0tdvQO5PPwEaDQDANjQEbmFhsO7Vq/oTZXKg3bAmqpKIGoIBhYiajZLLl5G+bj1yf/4ZENqJ1+yGD4ciLAzWgT2MXB0RGRIDChGZvOKLl5C+bh3yfvnlTjAZOQJuYWGwCggwcnVE1BgYUIjIZBWfP4/0teuQ9+uvum32Dz4IRdg8WHXtasTKiKixMaAQkckpOvc30tetQ/7+/doNkgT70aOgmDcPVp07G7c4ImoSDChEZDKKzp5FesRa5B88qN0gSXB46CEo5s2FZceORq2NiJoWAwoRGV1RTAxur12LgqjD2g0yGRweeRiKuXNh2b69cYsjIqNgQCnD1YyJml7h6dNIj1iLgj/+0G6Qy+E4diwUc+fAwt/fqLURkXFJQpR1iScA0K1mnJOTAwcHB2OXQ9QiFZ44gdtr16Lw2J/aDXI5HCeMh2LOHFi0bWvc4oio3gz5GcoWFCIyLI26ytlahRAoPB6N9LVrURgdrT3WzAxOjz4K1znPwcLHx7h1E5FJYUAhIsOJ21tpvRth74VC79m4/eNpFJ08pd1obg6nSROhmD0b5t7eRiqWiEwZAwoRGUbcXmDnNADau8ZCAAW3LJH+mxJF6ZsAAJK5OZwefwyus2fD3NPTiMUSkaljQCGihtOotS0nENpgkmKJ23/bozjDAgAgyQScuklw/TgS5p5exq2ViJoFBhQiariEoxA5ychPtkT63/YoziwLJnIB5w4FcOmWD3NrDVB8FQADChHdGwMKETWIEAL5Bw4ifZ8CxVnlwUQD546FcO2aDzNrzZ2D81ONVCURNTcMKERUL0KjQd5vvyF97TqUXLgAwAKSmQYuHQvg0rUAZlaayifZuTd5nUTUPDGgEFGdCI0Gefv2aYPJpUsAAJmNDZw7F8KlXRrMLKsIJpAABy/tkGMiolpgQCGiWhFqNXIjI5G+bh2UV64CAGR2dnCeOgWu06dDnhxVNopHQvlIHi1J+5/RqwCZvKnLJqJmigGFiGokVCrk/vwz0teth/LaNQCAzN4eLtOmwWXaVMgdHbUHOo0Dnvi80jwocPDShpOAcUaonoiaKwYUIqqSUKmQ8+OPyFi3HsqEBACAzNERLtOnwWXqVMjt7SufFDAO6PpwlTPJEhHVBQMKEekRpaXI2fsD0jdsQOmNGwAAuaMjXGbOhPOUpyG3s6v5AjI50G5YE1RKRC0ZAwoRAQCEUons779HxoaNKE1MBADInZ3hMmsmnJ/6P8jtbI1cIRG1JgwoRK2cUCqRvWs3MjZuRGmytu+I3NUVrrNmwfmpJyGzsTFyhUTUGjGgELVSGqUS2d9+i4xNn0CVkgIAkLsp4PrMM3CePBkya2sjV0hErRkDClEroykpQfY33yJj0yaoUrUzu5q5ucF19mw4PfE4ZFZWRq6QiIgBhajV0BQXI3vnTm2Lye3bAAAzd3e4PjcbTo89BpmlpZErJCK6gwGFqIXTFBUh6+sdyPj0U6jT0wEAZp6eUDw3G46TJkFmYWHkComIKmuxAaWwsBDdunXD448/jv/85z/GLoeoyWkKCpD19dfI+Gwz1BkZAABzLy+4zpkDp0cnQGIwISIT1mIDyptvvomBAwcauwyiJqfOL0DWV18hc/NmqLOyAADmPj5QzJ0Dx/HjIZmbG7lCIqJ7a5EB5fLly7hw4QLGjh2Lc+fOGbscoiahzs9H1hdfIHPzFqhzcgAA5n5toZgzF45jH2EwIaJmRWbsAiqKiorC2LFj4eXlBUmSsGfPnkrHREREwN/fH1ZWVggODkZ0dLTe/iVLluDtt99uooqJmoBGDcQfBs5+q/2vRq3bpc7Nxe2ICFy5fwRuf7AG6pwcWPj7w+udVejw009wmvgowwkRNTsNakEpLS3FrVu3UFhYCDc3N7i4uDS4oIKCAvTq1QuzZs3CxIkTK+3fsWMHFi9ejPXr1yM4OBgffPABRo0ahYsXL6JNmzb4/vvv0blzZ3Tu3BlHjx695+OVlJSgpKRE93Nubm6DnwORQcXtrXIBPvXQFcj88zYyt22DJi8PAGDRvj0U8+bB4aExkORc/4aImi9JCCHufdgdeXl5+OKLL/D1118jOjoaSqUSQghIkgQfHx88+OCDeO655zBgwICGFydJ2L17NyZMmKDbFhwcjAEDBuDjjz8GAGg0Gvj6+uL555/H0qVLsWzZMnzxxReQy+XIz89HaWkpXnzxRSxfvrzKx1i5ciVee+21SttzcnLg4ODQ4OdA1CBxe4Gd0wDc+d9UVSIh86Idsi7ZQqPSNoJaduoIxbx5sB81isGEiIwmNzcXjo6OBvkMrVNAef/99/Hmm2+iQ4cOGDt2LIKCguDl5QVra2tkZmbi3LlzOHz4MPbs2YPg4GB89NFH6NSpU/2LqxBQlEolbGxs8O233+qFlunTpyM7Oxvff/+93vlbtmzBuXPnahzFU1ULiq+vLwMKGZ9GDXzQQ9dyoiqRIfOCLbIu3xVMXADFq+9rg4nM5O7YElErY8iAUqdbPCdOnEBUVBS6d+9e5f6goCDMmjUL69evx+bNm3H48OEGBZSK0tPToVar4e7urrfd3d0dFy5cqNc1LS0tYckJqsgUJRwFcpOhKtYGk8wrthDlwcSpFIoeebD3LobU1Q5gOCGiFqZOAWX79u21Os7S0hJz586tV0GGNGPGDGOXQFRvqptXkHHGAVlXbCDU2gBi5ayEokce7LxKIEllB+anGq9IIqJGUu9OskVFRRBCwKZspdOEhATs3r0b3bp1w6hRowxW4N0UCgXkcjlSU/X/QU5NTYWHh0eDrh0REYGIiAio1ep7H0zUiErT0pD56afI2r4dQmkHALByKQsmnncFk3J27pUvQkTUzNW7XXj8+PH4/PPPAQDZ2dkIDg7Ge++9hwkTJmDdunUGK/BuFhYW6NevH/bv36/bptFosH//fgwaNKhB1w4PD0dcXBxOnDjR0DKJ6qU0NRW3/v0mro58AJlbP4dQlsK6DeAbmgn/B9Jh71UxnEiAgzfgN9jgtfj7+6Nnz57o3bs3hg8fXu1xP/74I7p06YJOnTrhk08+AQDcvHkT9913HwICAtCzZ0988803Bq+PiFq+eregnD59Gv/9738BAN9++y3c3d1x5swZfPfdd1i+fDnmzZtXr+vm5+fjypUrup/j4+MRExMDFxcXtG3bFosXL8b06dPRv39/BAUF4YMPPkBBQQFmzpxZ36dCZFSlKSnI2LQJ2d98C1FaCgCw7tsXivAw2DqlQ/pmOgAJd4/k0f4MYPQqQNY4o3aOHj0KOzu7averVCosXrwYBw4cgKOjI/r164dHH30UZmZm+OCDD9C7d2/cunUL/fr1w0MPPQRbW9tGqZOIWqZ6B5TCwkLY29sDAPbt24eJEydCJpNh4MCBSEhIqHdBJ0+e1PuLbfHixQC0I3W2bNmCyZMn4/bt21i+fDlu3bqF3r17IzIyslLHWSJTV5qUhPSNm5C9axdQFkxs+veHYn44bIKDIZU3l0ifVzkPCkavAgLGGaFyrejoaHTv3h3e3t4AgDFjxmDfvn146qmn4OnpCQDw8PCAQqFAZmYmAwoR1Um9A0rHjh2xZ88ePProo/jll1/wwgsvAADS0tIaNLTovvvuw71GPs+fPx/z58+v92NUhX1QqKkoExORsWEDsnfvAVQqAIBNcDAUYWGwDQ6qfELAOKDrw9pRPfmp2j4nfoMbreUE0A7xDw0NhUwmw6JFi/D0009XOiY5OVkXTgDA29sbSUlJesecOnUKarUavr6+jVYrEbVM9Q4oy5cvx//93//hhRdewIgRI3R9QPbt24c+ffoYrMCmEh4ejvDwcN0YbiJDU964gfT1G5Dz/fdAWRC2HTwIirAw2PTvX/PJMjnQblgTVKl15MgReHt7IyUlBSNHjkRgYCB69uxZp2tkZmZi2rRp2LRpUyNVSUQtWb07yT722GO4ceMGTp48icjISN32ESNG6PqmEBFQEh+P5JeX4uqYh5CzaxegVsN26FD4ffUV2n722b3DiYFdvXoVkiThxx9/xIgRI2BjY4MuXbrg+PHjumPKW0Y8PT3x0EMP4fTp05Wu4+XlpddikpSUBC8vLwDaCRAnTJiApUuXYvBgw3fiJaKWr94B5ebNm/Dw8ECfPn0gu2uSqKCgIHTt2tUgxRE1ZyXXriHpHy/h2sOP6FpNbEND4P/1drT9ZBNs+hqmpXHFihXo3LkznnzySeTk5OCHH35Ar1698Omnn1Z5fGxsLCRJwvvvv49XX30VsbGxaNu2LZYuXQpAux5WXtnaPvn5+fj999+rnJwxKCgI586dQ1JSEvLz8/Hzzz9j1KhREEJgxowZuP/++zF16lSDPEcian3qfYvHz88PLi4u6NWrF3r37q37UiqV+PDDD7F161ZD1knUbJRcuYL0teuQ+/PPQFl/Krvhw6EImwfrwECDPtZvv/2GtLQ0nDp1CmvXrsWECROQnZ2Nb775Bp07d67ynNjYWDg5OWHHjh1wc3MDAIwbNw4bNmwAoJ1X6NFHHwUAqNVqzJ49W7e2Vu/evRETEwMAMDMzw3vvvYfhw4dDo9HgpZdegqurK44cOYIdO3agZ8+eutXIt23bhkADP3ciatnqvFhguYSEBJw5cwYxMTE4c+YMTp8+jeRk7SgDBwcHZGVlGbTQxnZ3J9lLly5xLR6qs+KLl5C+bh3yfvnlTjAZOQKKefNgXc3yEA21evVqDB06VNcHLDg4GAsWLKiyU2u5CRMmwNnZGZs3b9ZtW7x4Ma5du6YLFERE9WHItXjqfYvHz88PEyZMwMqVK/H999/j5s2bOHLkCDp06NBoE7U1Jk7URvVVfOECEp9fgPjx45EXGQkIAfsHH0S73bvg+/HHjRZOAKBr167Yu3cvACAyMhKSJOGdd95BSkpKtefExsZWmtgwJiYGvXv31ttW/rfLypUr9X4mImoKBl1hbNCgQVizZk2NqwcTtRRFf/+Nm+HzET/hUeT9+isgSbAfMxrtvv8ePh+ugVW3bo1ew9ixY1FUVIS2bdvirbfewq5du/DSSy9h4MCBVY6eycnJwfXr1yuNtKsqoKxbtw4bN25EQUEBli5diqioqMZ8KkREeup9i0epVMLCwqLS9itXriAwMBBFRUUNLs4YDNk8RS1T0dmzSI9Yi/yDB7UbJAkODz0Exdw5sDTg6t2N4fDhw7j//vuRl5cHKysrANrbtf7+/oiPj4e/v7/e8atXr8Yrr7yC33//HUOGDDFCxUTUnBjyM7TenWTt7OwQEBCAPn36oHfv3ujTpw+8vLzw0UcfYeTIkQ0qisgUFcXG4nZEBAqiDms3yGRweORhKObOhWX79sYtrpZiY2PRpUsXXTgBgDNnzsDJyalSOFm/fj0cHR2xYMEC/PDDD9BoNBg2rOnmYiGi1q3eLShHjhxBbGwsYmNjERMTg3PnzqG4uBgAMHr0aPTv3x+BgYEIDAxsVsOO2YJCFRWePoP0iAgU/PGHdoNcDsexY+E65zlYtmtn3OIakRACkiRh5cqVWLlype5nIqLqGPIztN4BpSKNRoOLFy8iJiYGMTExuvCSlpbWLKaP5ygeqqjw5EncjohA4bE/tRvkcjhOGA/Fc8/Bws/PuMUREZkgkwwo1UlNTW1WC/mxBYUKjkcjPSIChdHR2g1mZnB69FG4znkOFj4+xi2OiMiEGa0Pyo0bN9C2bdtaH5+UlKS3mBiRqRJCoPDPP5EesRaFJ09qN5qbw2nSRChmz4Y5f4+JiJpUnYYZDxgwAHPmzKlxrpCcnBxs2rQJPXr0wHfffdfgAokaTKMG4g8DZ7/V/ldz55ajEAL5R/5AwtNTcGPmLBSePAnJ3BzO//cUOu77BZ4rVzKcEBEZQZ1aUOLi4vDmm2/igQcegJWVFfr16wcvLy9YWVkhKysLcXFx+Pvvv9G3b1+8++67eOihhxqrbqLaidsLRL4M5Cbf2ebgBTFqFQoynJEesRZFsbEAAMnCAk5PPAHX2c/CvBndliQiaonq1QelqKgIP/30E44cOYKEhAQUFRVBoVCgT58+GDVqFHr06NEYtTYJ9kFpQeL2AjunAbjzKy4EkJ9shfS/7VCcqZ3HR7K0hPOTk+HyzDMwb9PGSMUSETV/zaqTbHPBUTwtjEYNfNBD13IiBJCfVBZMssqCiRngPGU6XJ95BmZli+YREVH9MaA0IragtBDxh4Gtj0AIIC/RCul/26Mk2xwAIJlp4NKxAC5dC2A2Zy/QjpOPEREZgknMJEtkykRuCvJulAWTHG0wkZlp4NypLJhYarQH5qcasUoiIqoOAwq1KEKtRm5kJNLXrIfyhgsAQGZeFky65MPMskKDoR07wxIRmSIGFGoRhFqN3P/9D+nr1kN57RoAQGYBuHTKg0uXfMgtKt7JlAAHL8BvcNMXS0RE98SAQs2aUKmQ+9NP2mBy/ToAQOboCJfp0+AywBXyn+ZUcVbZejKjVwEyeZPVSkREtdeggHL48GFs2LABV69exbfffgtvb29s27YN7dq1w9ChQw1VI1ElorQUOT/8iPQN61GacAMAIHd0hMvMmXCe8jTkdnbaA22tq5wHBaNXAQHjjFA5ERHVRr0DynfffYepU6fi6aefxpkzZ1BSUgJAO5PsW2+9hf/9738GK5KonCgtRc733yN9w0aU3rwJAJA7O8Nl1kw4P/V/kNvZ6p8QMA7o+jCQcFTbIdbOXXtbhy0nREQmrd7DjPv06YMXXngB06ZNg729PWJjY9G+fXucOXMGY8aMwa1btwxda5PgMGPTJJRKZO/eg4wNG1CarG0Nkbu6wnXWLDg/ORkyW9t7XIGIiBqbSQwzvnjxIkJCQiptd3R0RHZ2dkNqMoq7J2oj06FRKpHz3XdI37gJqpQUAIDcTQHXZ56B8+TJkFlbG7lCIiJqDPUOKB4eHrhy5Qr8/f31th85cgTt27dvaF1NLjw8HOHh4br0R8alKSlB9jffImPTJqhStXOVmLm5wXX2bDg98ThkVlZGrpCIiBpTvQPK7NmzsXDhQnz22WeQJAnJyck4duwYlixZgldffdWQNVIroikuRvbOncjY9AlUt28DAMzc3eH63Gw4PfYYZJaWRq6QiIiaQr0DytKlS6HRaDBixAgUFhYiJCQElpaWWLJkCZ5//nlD1kitgKaoCFlf70DGp59CnZ4OADDz9ITiudlwnDQJMgsLI1dIRERNqcFr8SiVSly5cgX5+fkICAiAXfnwzmaKnWSblqagAFlff42MzzZDnZEBADD38oLrnDlwenQCJAYTIqJmwyQ6yQJAcXEx/vrrL6SlpUGj0eiN3Bk3jnNMUPXU+QXI+uorZG7eDHVWFgDA3McHirlz4Dh+PCRzcyNXSERExlTvgBIZGYmpU6cio+yv3rtJksTRMFQldX4+sr74Apmbt0CdkwMAMPdrC8WcuXAc+wiDCRERAWhAQHn++efxxBNPYPny5XB354JrVDN1bi4yt21D5tbPocnNBQBY+PtDMW8uHB5+GJIZV10gIqI76v2pkJqaisWLFzOcUI3UOTnI3Po5MrdtgyYvDwBg0b49FPPmweGhMZDknNGViIgqq3dAeeyxx3Dw4EF06NDBkPVQC6HKykLm558ja9sX0OTnAwAsO3WEYt482I8axWBCREQ1qvconsLCQjz++ONwc3NDYGAgzCv0HViwYIFBCmwqd88ke+nSJY7iqSdVVhYyN29B1hdfQFNYCACw7NwZirAw2D/4ACSZzMgVEhFRYzHkKJ56B5RPP/0Uc+fOhZWVFVxdXSFJ0p2LShKuXbvWoMKMhcOM60eVkYHMzZuR+dV2iPJg0q0bFGHzYD9iBIMJEVErYBLDjP/1r3/htddew9KlSyHjh0+rpbp9GxmfbUbW119DFBUBAKwCAqCYHw674cP1gisREVFt1TugKJVKTJ48meGklSpNS0Pmp58i6+sdECUlAACrwEAowsNgFxrKYEJERA1S74Ayffp07NixA//85z8NWQ+ZuNLUVGRs+gTZO3dCKJUAAOtevaCYHw7boUMZTIiIyCDqHVDUajXeffdd/PLLL+jZs2elTrLvv/9+g4sj01GakoKMTZuQ/c23EKWlAADrvn2hCA+D7eDBDCZERGRQ9Q4oZ8+eRZ8+fQAA586d09vHD6tmRKMGEo4C+amAnTvgNxiQ3RkCXJqUhPSNm5C9axdQFkxs+veHYn44bIKD+V4TEVGjqHdAOXDggCHrIGOI2wtEvgzkJt/Z5uAFjH4HSoe+yNiwAdm79wAqFQDAJjgYirAw2AYHGadeIiJqNTi/eGsVtxfYOQ2A/ihzZVIa0l98HjkJdoBGAwCwGTQQbuHhsOnf3wiFEhFRa1SngLJ48WK88cYbsLW1xeLFi2s8ln1QTJhGrW05uSucKPPkSP/bHjkJ1oCQAGhgO2QwFOHhsOnb12ilEhFR61SngHLmzBmUlvVDOHPmTLXHsV+CiUs4qruto1EBKSeckHujPJgAtp7FcOueB+sXZwDtGE6IiKjp1Smg3N3vZOvWrfDx8ak0D4oQAjdv3jRMddQ48lN130pyoLRADggJdl7FUHTPg7VraaXjiIiImlK9+6C0a9cOKSkpaNOmjd72zMxMtGvXDmq1usHFUSOxu7MCtSQB7n1zAQDWLqXVHkdERNSU6j0NbHVL+OTn58PKyqreBVET8BusHa0D7S0da5fSCuFEAhy8tccREREZQZ1bUMo7x0qShOXLl8PGxka3T61W4/jx4+jdu7fBCqRGIJMDo98pG8UjQX8kT1n/odGr9OZDISIiakp1DijlnWOFEDh79iwsLCx0+ywsLNCrVy8sWbLEcBU2kYiICERERLSeW1MB44AnPq9mHpRV2v1ERERGIonq7tXcw8yZM7FmzZoGL6dsagy5VHSzcI+ZZImIiGrLkJ+h9Q4oLVWrCyhEREQGYsjP0Hp3kgWAw4cPY8qUKRg0aBCSkpIAANu2bcORI0caVBQRERG1bvUOKN999x1GjRoFa2trnDlzBiUlJQCAnJwcvPXWWwYrkIiIiFqfegeUf//731i/fj02bdoEc3Nz3fYhQ4bg9OnTBimOiIiIWqd6B5SLFy8iJCSk0nZHR0dkZ2c3pCYiIiJq5eodUDw8PHDlypVK248cOYL27ds3qCgiIiJq3eodUGbPno2FCxfi+PHjkCQJycnJ+PLLL7FkyRLMmzfPkDUSERFRK1PvtXiWLl0KjUaDESNGoLCwECEhIbC0tMSSJUvw/PPPG7JGIiIiamXqNA/KX3/9hR49euitYKxUKnHlyhXk5+cjICAAdnZ2jVJoU+E8KERERPVjtHlQ+vTpg/T0dABA+/btkZGRAQsLCwQEBCAoKKjZhxMiIiIyDXUKKE5OToiPjwcAXL9+HRqNplGKIiIiotatTn1QJk2ahNDQUHh6ekKSJPTv3x9yedXrtly7ds0gBRIREVHrU6eAsnHjRkycOBFXrlzBggULMHv2bNjb2zdWbURERNRK1XkUz+jRowEAp06dwsKFCxlQiIiIyODqPcx48+bNhqyDiIiISKdOAWXx4sV44403YGtri8WLF9d47Pvvv9+gwoiIiKj1qlNAOXPmDEpLS3XfV0eSpIZVRURERK1anQLKgQMHqvz+bomJiXj99dcbVlUDZGdnY+TIkVCpVFCpVFi4cCFmz55ttHqIiIio7uo0k2xtxMbGom/fvlCr1Ya8bK2p1WqUlJTAxsYGBQUF6NGjB06ePAlXV9danc+ZZImIiOrHaDPJNgdyuRw2NjYAgJKSEgghYOAMRkRERI3M5AJKVFQUxo4dCy8vL0iShD179lQ6JiIiAv7+/rCyskJwcDCio6P19mdnZ6NXr17w8fHBP/7xDygUiiaqnoiIiAzB5AJKQUEBevXqhYiIiCr379ixA4sXL8aKFStw+vRp9OrVC6NGjUJaWpruGCcnJ8TGxiI+Ph5fffUVUlNTm6p8IiIiMoA6z4MyceLEGvdnZ2fXtxYAwJgxYzBmzJhq97///vuYPXs2Zs6cCQBYv349fvrpJ3z22WdYunSp3rHu7u7o1asXDh8+jMcee6zK65WUlKCkpET3c25uboPqJyIiooarcwuKo6NjjV9+fn6YNm1aY9QKpVKJU6dOYeTIkbptMpkMI0eOxLFjxwAAqampyMvLAwDk5OQgKioKXbp0qfaab7/9tl79vr6+jVI7ERER1V6dW1CMOYNseno61Go13N3d9ba7u7vjwoULAICEhAQ899xzus6xzz//PAIDA6u95rJly/QmncvNzWVIISIiMrJ6T3VvqoKCghATE1Pr4y0tLWFpadl4BREREVGdmVwn2ZooFArI5fJKnV5TU1Ph4eHRoGtHREQgICAAAwYMaNB1iIiIqOGaVUCxsLBAv379sH//ft02jUaD/fv3Y9CgQQ26dnh4OOLi4nDixImGlklEREQNZHK3ePLz83HlyhXdz/Hx8YiJiYGLiwvatm2LxYsXY/r06ejfvz+CgoLwwQcfoKCgQDeqh4iIiJo/kwsoJ0+exPDhw3U/l3dgnT59OrZs2YLJkyfj9u3bWL58OW7duoXevXsjMjKyUsdZIiKilk6tEYiOz0RaXjHa2FshqJ0L5LKWsWCvwdfiae64Fg8RETUHkedS8NoPcUjJKdZt83S0woqxARjdw9MoNXEtnkbATrJERNRcRJ5LwbwvTuuFEwC4lVOMeV+cRuS5FCNVZjhsQamALShERGTK1BqBoe/8XimclJMAeDha4cjL9zf57R62oBAREbVS0fGZ2nAilcDM7m9YenwHC5dDuv0CQEpOMaLjM41XpAGYXCdZIiIiqlpiXiL2xv8Aa9/fILe5CkmmBgCoS9pAmRmqd2xaXtUtLM0FAwoREZGJUmlUiEmLQVRiFKISo3A15yoAwMxOu1+jdIEqvytU+V2hbTu5c0unjb1V0xdsQAwoZSIiIhAREQG1Wm3sUoiIqAWp61Dg7OJsHEk+gqibUTiSfAR5yjzdPrkkR+82fXD2khey0jtCo3TD3aEEuNMHJaidSyM9o6bBTrIVsJMsEREZSm2GAgshcDn7sq6VJPZ2LDRCozveydIJQ72HItQnFIO9B8PBwkE3igfQtpuUK48q66b0NcpQY0N+hjKgVMCAQkREhlAeIip+yEoAIJViwSMSCuRnEZUYhZQC/WHBnZ07I8QnBKE+oQhUBEIuk1d5/ZY8DwoDSgUMKERE1FBVDQWWzHJgZncBZnbnIbe9CklWqttnKbdEsGcwQrxDEOITAk+72gUMU5tJ1pCfoeyDQkREZGDaocCFkFnfLAslFyC30m8l0ZQ6YrhvKJ7oPgoDPAbA2sy6zo8jl0kY1MHVUGWbFAaUMuwkS0REDZWrzMXR5KP48u9I2Hb6EzKzAt0+ISRoinyhyu8GVX5XaEo88ED/Pgjx8TZixaaLt3gq4C0eIiKqLSEE4nPjcTjxMA4lHsKZ1DNQCdWd/WorqAo6Q5XfBer8LhBqO73zt88e2KJaQHiLh4iIyEiUaiVOpp7UhZKbeTf19rdzbIdh3iHYGWWH9NteEKjcwbWlDAVuTAwoRETUatW2k2l6UToOJx5GVGIUjiYfRaGqULfPXGaO/u79EeobihDvEPg6+AIAAq21o3gkVD0UeMXYAKN2aDV1DChERNQq1TRM98Hu7jifeR5RN6NwKPEQ/s74W+9chbUCw7yHIdQnFAO9BsLW3LbS9Uf38MS6KX0rPYaHkYcCNxfsg1IB+6AQEbV8Vc5RIpXAzPYKzOzPw0VxDXkq/cX2urt2181N0s21G2RS7dbbNbWhwI2JfVCIiIjqSa0ReO2HOO3KNeYZd4YB21zTLb6XpwKszawx2GswQn1CMdR7KNxs3Or1eC15KHBjYkApw2HGREQtX6mmFF/GRCHTchds2l+A3DJNb/+dxfe6IeLJyQjpxNswxsJbPBXwFg8RUcuSVZyFI0lHEJUYhT+S/9BbfE8IGdSFflDld4U6v5ve4ntrnuyN8b05R0ld8BYPERFRNYQQuJR1CYeTDuPQzUP4K/0vvcX3bM0ckJ3eQdtSUtAZ0FQ9g2sbe6umKpmqwIBCREQmpT6dSotVxYi+FY2oRO2om1sFt/T2d3bujFCfUIT4hCDApQdCVx/CrbziKq/FOUpMAwMKERGZjLqs0Hur4BaiEqMQlRiF4ynHUay+c0754nuhPqEY5j2s0uJ7K8YGcI4SE8c+KBWwDwoRkXFUOfQXd0JDxNO94OWergslF7Mu6h3nYeuhayWpzeJ7dQlDVDuG/AxlQKmAAYWIqOmpNQJD3/ldLywAAGRFMLO9DDO78zC3vwTI7yy+J0FCL7de2hlcfULQyakTJKlurR6taY6SpsBOskRE1KJEx2eWhRMBmcVtyHVzk1yHJN3p4Gott0Wo7zCE+IRgqPdQOFs5N+hxOUeJ6WJAKcN5UIiIjEOpVuJo8lFYuv8IM7uLkFlk6O1Xl7SBOr8rVPld8drYcZjYx89IlVJT4i2eCniLh4io8ZUvvnco8RCOJR/TW3xPaORQF7YvmzCtK0TpnRaO7bMHssXDhPEWDxERNSsaocH5jPO6YcBVLb6Xnd4ReZmdoCroCAhLvf0c+tv6MKAQEVGjKCgtwJ/JfyIqSTvqJr0oXW9/D9ceCPENQYhPCLq5dMO+v1M59Jd0GFCIiKhGdRnpcjPvpm4Y8IlbJ1CqKdXtszGzwWCvwQjxCcEwn2FQWCv0zh3dwxPrpvStNPTXg0N/WyUGFCIiqta95gop1ZQiJi1GF0qu5VzTO9/X3lc3N0k/936wkFvU+Hije3jigQAPDv0ldpKtiJ1kiYi0qp04TV4AM9uLCO6Rimv5p5FXemfxPTPJDH3d+yLER3vrxt/Bv85zk1DzxU6yRETUqNQagdd+iCsLJwIyy1swK5ubRGZ9A5IkEJulPdbZ0hnDfIZhmM8wDPYaDAcL/nFHDceAQkRElRy5koI0dQwsPc5rQ4l5jt5+dbEnVPldseL+x/BkryGQy+RGqpRaKgYUIiICcGfxvUOJh3A06U/Y+Cp1+4TGHOqCDlDld4MqvwuEygkAYCd1YDihRsGAUoYzyRJRa6PWqHE2/awulFzKuqS3X1PqBFWedrI0dWEHQJhXukYbe6umKpdaGXaSrYCdZImoJctV5uJo0lEcSjyEI0lHkF2Srdsnk2To5dZLu86N1zDM2HgTqTkllTrJAncmTjvy8v0cYUM67CRLRES1IoRAfG48om5qW0nOpJ2BWtxpKba3sMdQr6EI8Q3BUK+hcLJy0u1bOdaBE6eR0TCgEBE1A3WZLE2pVuLkrZOISorCoZuHkJifqLe/g2MH3TDg3m16w0xW9UcBJ04jY2JAISIycfeaLA0AbhfexuGkw4hKjMLR5KMoUhXpjjWXmSPIIwjDfIYhxCcEvva+tX5sTpxGxsI+KBWwDwoRmZJqJ0uDBjKrZIwfnI1k5ZkqF98L9QnFMJ9hGOQ5CDbmNk1XNLVa7INCRNQK6E+WBkBWAjPby5DbXYCZ3UXIzPLwa8qd4ysuvieTZMYom8ggGFCIiExUdHwmbhUmwdxZO4Or3PYaJOlOB1ehtoCqoBOe7fcIZvQZXWnxPaLmjAGFiMiE3L343k9X9sOu4029/RqlK1T5XaHK6wp1UTtAmKHTsN4MJ9TiMKAQERlZVnEWjiQdKZvB9aje4ntCyKAu9NeGkvxuEEoF7gz01eJkadQSMaAQETUxIQQuZV3SzeD61+2/IO7qBlu++N4Qr2F4fWcpUrNlNU6WFtTOpclqJ2oqDChERE2gSFWEE7dO4NDNQziUeAiphal6+7s4d9HNTRKoCNStbyN7JIWTpVGrxIBCRNRIUvJTEJUYhaikKBxPOY4SdYlun5XcCsGewbpQ4mHrUeU1OFkatVYMKERENajLDK7li+8dSjyEqMSoSovvedp66gJJkEcQrMxq13eEk6VRa8SAQkRUjdrM4FrbxfdCfELQyakTJKl+oUIukzCog2uDng9Rc8KAUiYiIgIRERFQq9X3PpiIWrzqZnC9lVOEsJ2ReOJWHlJVZ+q0+B4R1R6nuq+AU90TkVojMPSd3++0nEgqyG2uwcxOO2GazCJT7/gOjh20M7h617z4HlFLx6nuiYgaUXR8Jm4VpMHc8SLk9udhZnsFkkyp2y80cqgLO2BKz9GY3nsMfOx9jFgtUcvEgEJEBEAjNDifcR6HEg/h+0u/wa7TZf39pfZQ5XeFOr8rVAUdAWGJHkN6w8fe20gVE7VsDChE1GoVlBbgWPIx7VDgxChkFGfo7VcX+ZbN4NoVmmIvcAZXoqbDgEJErcrN3Ju6YcAnUk9ApVHp9tmY2WCI9xAM9Q7BO7uAtCxzzuBKZCQMKETUopUvvlc+g+v13Ot6+9vat0WITwhCfUPRr00/mMvNAQA2D3MGVyJjYkAhohYnszgTR5KOICoxqtLie2aSGfq598Mwn2EI9QmFv6N/ldfgDK5ExsWAQkTNXvnie4cSta0kZ2+f1Vt8z8XKBUO9hyLEJwSDvQbD3sK+VtflDK5ExsOAQkQmoy7TyhepihCdEq3rT1Jx8b2uLl11M7j2cO2hW3yvrjiDK5FxMKAQkUmozbTyyfnJuhE30beiKy2+N9BzIEJ8QzDMe1i1i+8RUfPAgEJERlf9tPKFCP9uFx5OykJiyWlcztKfm6R88b1Qn1AM8BhQ68X3iMj0MaAQkVGpNQKv/RB3J5zICmFmdxlmdudhZnsJklkhDtwq2yXJ0Nutt66Da0enjvVefI+ITBsDChEZ1fFrGUgtSoC5i3adG7lNAiRJo9sv1FZQ5XfB3AFjMaPPKC6+R9RKMKAQUZMrUZfg5K2TOJR4CJHXDsC2wy29/eqSNlDldYM6vyvURW0ByNHOujfDCVErwoBCRE0irTANhxMP41DiIfyZ8ieKVEW6feWL75VPKy9KK8/QymnliVoXBhQiahQaocHf6X8jKikKh24ewvnM83r73azdEOITgqHeIXhlezFSswWnlSciHQYUIjKYfGU+jqVoF987nHhYb/E9CRJ6KHro5ibp5tJN18FV9QinlScifQwoRNQgN3Jv6CZLO5l6Um/xPVtzWwz2GlzWUjIUCmtFldfgtPJEVFGLCyg3b97E1KlTkZaWBjMzM7z66qt4/PHHjV0WUYtRqinFmdQzulBS1eJ7ob6hCPEJ0Vt87144rTwR3U0SQlR127fZSklJQWpqKnr37o1bt26hX79+uHTpEmxtbWt1fm5uLhwdHZGTkwMHB4dGrpaoeShffO/QzUM4mnwU+aX5un3li++V37qpbvE9Imr5DPkZ2uJaUDw9PeHpqW0O9vDwgEKhQGZmZq0DClFLVZd1boQQuJh1EVGJUTUuvhfqE4pBXoNqvfgeEVFtmVxAiYqKwurVq3Hq1CmkpKRg9+7dmDBhgt4xERERWL16NW7duoVevXrho48+QlBQUKVrnTp1Cmq1Gr6+vk1UPZFpqs06N0WqIhxPOa5b66bi4nvdXLrpZnDtoegBmSRr0udARK2LyQWUgoIC9OrVC7NmzcLEiRMr7d+xYwcWL16M9evXIzg4GB988AFGjRqFixcvok2bNrrjMjMzMW3aNGzatKkpyycyOdWvc1OMsK/3Y8r9BUjXxFZafM/azBrBnsEI8eHie0TU9Ey6D4okSZVaUIKDgzFgwAB8/PHHAACNRgNfX188//zzWLp0KQCgpKQEDzzwAGbPno2pU6fW+BglJSUoKbnzj3Jubi58fX3ZB4VaBLVGYOg7v9/VcqKB3PoG5HZl08pb6c/g6mXrpV18z1e7+J6l3LLpiyaiZqvV9kFRKpU4deoUli1bptsmk8kwcuRIHDt2DID23vmMGTNw//333zOcAMDbb7+N1157rdFqJjKm6PhMpORlwsxBu/ie3PYSZGaFuv1CSFAX+eGJgAcxtecYdHDqwMX3iMgkNKuAkp6eDrVaDXd3d73t7u7uuHDhAgDgjz/+wI4dO9CzZ0/s2bMHALBt2zYEBgZWec1ly5Zh8eLFup/LW1CImishBK7lXENUYhR2XfgVdp3/rrD4njVU+V3KppXvDGhs0HdQb3R09jZi1URE+ppVQKmNoUOHQqPR3PvAMpaWlrC0ZDM2NW93L74XlRiFpPwk3T5JAtTF7lDld9VbfO9uXOeGiExNswooCoUCcrkcqan6owtSU1Ph4cEOfNS61LT4noXMAgM8B2CYdyg+/EGOtExbrnNDRM1KswooFhYW6NevH/bv36/rOKvRaLB//37Mnz+/QdeOiIhAREQE1Gq1ASolMrx7Lb7XxrqNbhhwsGcwbMxtAACuaq5zQ0TNj8kFlPz8fFy5ckX3c3x8PGJiYuDi4oK2bdti8eLFmD59Ovr374+goCB88MEHKCgowMyZMxv0uOHh4QgPD9f1QCYyBfdafC9QEaibwbWrS9cqO7hynRsiao5MbpjxwYMHMXz48Erbp0+fji1btgAAPv74Y91Ebb1798aHH36I4OBggzw+p7onY6vN4nuhPqEY6j0Urtautb5uXWaSJSKqD0N+hppcQDE2BhRqavdafM/fwV9366Zvm761XnyPiKiptdp5UBoT+6BQU7rn4nse/RDqo10R2M/Bz4iVEhEZB1tQKmALCjWG8sX3Dt08hKikqCoX3xvmPQyhvqEY5DkIdhZ2RqyWiKh+2IJCZCR16cdRvvhe+a2btMI0vf3dXLppp5X3CUV3RXcuvkdEdBcGFKJaqs2KwMn5yYhKjMKhxEOITomGUqPUHWttZo2BngN1i++527pXegwiItJiQCGqhepXBC5A+He78HBSFm4Wn8KV7Ct6+73tvHXDgLn4HhFR7TGglGEnWaqOWiPw2g9xd8KJrBBmdpdgZncBZraXIJkV4kDZosBySY5ebr0Q6huKEO8QLr5HRFRP7CRbATvJUkVHr6Rjyud7YWZ3AXL785Bb36hy8b15QWMxo88oOFpyoj8iap3YSZaokZWoS3Di1glEJUbh52u/w7aD/vpPdxbf6wZ1kS8AOfytejOcEBEZCAMKUZm0wjREJUYhKjGq0uJ7QmMGdWEHqPK7QpXXFULlXOl8rghMRGQ4DCjUapUvvlc+DLjS4ns2bbQjbrxC8M/tRUjNFk26InB2djZGjhwJlUoFlUqFhQsXYvbs2VUe++OPP+LFF1+ERqPByy+/jGeffbbG7UREpo59UMrc3Un20qVL7IPSQpUvvnfo5iEcTjqMzOJM3T4JEgLdAhHiHYJQ31B0ce6i6+BaPooHqHpF4HVT+hp80T21Wo2SkhLY2NigoKAAPXr0wMmTJ+Hqqr/+jkqlQkBAAA4cOABHR0f069cPR48ehaOjY5XbK55PRGQo7IPSCLiaccuVkJugm8H1VOopvcX37MzttIvv+YZiiNeQahffM8aKwHK5HDY2NgCAkpISCCFQ1d8T0dHR6N69O7y9vQEAY8aMwb59++Dn51fl9qeeesrgtRIRGRoDCrU4pepSnE47jUOJh3A48XCVi++Vz+Dax70PzGW1W3xvdA9PPBDg0aQrAmdnZyM0NBSXL1/G6tWroVAoKh2TnJysCyEA4O3tjaSkJJibm1e5nYioOWBAoRYhoyhDu/he4iEcSz6mv/iezAz93fvrJkxryOJ7cpmEQR2a7haJk5MTYmNjkZqaiokTJ+Kxxx6DuztnoCWilo8BhZolIQQuZF7Qjbo5m1558b3yQGKqi+9dvXoVHTt2xA8//ID//ve/OHbsGHx9ffH5558jODhY71h3d3f06tULhw8fxmOPPaa3z8vLS69lJCkpCUFBQdVuJyJqDthJtgJO1Ga6CksLcTzlOKKSoqpdfC/UNxShPqEIcA1o8sX3VqxYge3bt6Nv377YsGEDoqKi8Morr2DBggV45plnKh2/a9cuPPbYY7jvvvuwfPlyeHt7IywsDCqVCgcOHEBqaipsbGxgb2+PnJwcDBkyBNu3b0dgYKDedVQqFbp164aDBw9W6iRb1XZ2kiWixsJOstRqJOUn6RbfO5FyosrF90J9QjHMZxja2LQxWp2//fYb0tLScOrUKaxduxYTJkxAdnY2vvnmG3Tu3LnKc2JjY+Hk5IQdO3bAzc0NADBu3Dhs2LABAJCQkIDnnntO1zn2+eef14WT3r17IyYmBgBgZmaG9957D8OHD4dGo8FLL72kCyHVbSciMnVsQSnDYcamQaVRIfZ2rO7WTVWL74X6hCLEJwT9PfqbzOJ7q1evxtChQzFo0CAAQHBwMBYsWICnn3662nMmTJgAZ2dnbN68Wbdt8eLFuHbtGvbs2dPYJRMRGRxbUBoBhxkbllojaj3aJackB0eSjiAqMQpHko4gV5mr2yeX5OjdprculLR3bG+Si+917doVe/fuxaBBgxAZGQlJkvDOO+/g/vvvh6dn1UOQY2NjsWzZMr1tMTExCAkJ0dsmhIAkSVi5ciVWrlyp+5mIqCVjQCGDizyXUmm+EM+75gsRQuBq9lXdDK4xt2OgEXcW33O0dMRQ76EI9QnFYK/BzWJ9m7Fjx2L//v1o27Yt/P39sWvXLvz+++8YOHAgXnnllUozwObk5OD69evo06eP3vaYmBgsWLBAb9u6detgZmaGgoICLF26FGPGjEFoaGijPyciImPiLZ4K2Em2YcpnXK34SyVJpZDbXMOIfhm4XngSyQXJevs7OXfSzeDaU9ETcpm86Yo2gsOHD+P+++9HXl4erKy0a/gkJCTA398f8fHx8Pf31zt+9erVeOWVV/D7779jyJAhRqiYiOjeeIuHTJJaI/DaD3G6cCKZ5cDM7iLkdudhZnsFkqwUR29r91nILBDsGazr4Opl52W0uo0hNjYWXbp00YUTADhz5gycnJwqhZP169fD0dERCxYswA8//ACNRoNhw4Y1ccVERE2LLSgVsAWl/v64chtTv/gWZvYXYGZ3AXIr/VYSTakDVPndsHjIeEzrMxLWZtZGqrR5YR8UImou2IJCJiNfmY+jyUdxKPEQ9idEwbZdtm6fEBI0Rb5Q5XeFKr8rNCWeACR4mPdmOKmD8jCycuVKvZ+JiFoyBhSqs+s513XDgE+lnoJK3Fl8T6gtoSroAlV+F6jzu0CoK8/g2sbeqtI2IiKiuzGglLl7HhTSV6ouxam0U7pQkpCboLff38EfoT6hGOo9DIu25iA1p7RSJ1kAkKBd/TeonUuT1E1ERM0X+6BUwD4oWhlFGTicdBhRiVE4mnwUBaUFun3li++Vz03S1qGtbl/5KB4AeiGl/KbEuil9MbpH1fOCEBFR88Y+KGRw5Yvvlc9Nci79nN7ie65WrhjmMwyhPqEY5DUItua2VV5ndA9PrJvSt9I8KB53zYNCRER0LwworVhhaSH+TPkTUYlROJx4GGlF+ovvBbgG6FpJ6rL43ugennggwKPWM8kSERFVxIDSyiTlJ+HQzUOISoqqcvG9QZ6DEOobimHew+Bm41bvx5HLJAzqwIXpiIiofhhQWrjyxfcOJR7C4cTD1S6+F+oTiv4e/WEhtzBSpURERHcwoLRA5YvvHUo8hD+S/qi0+F6fNn0Q4hOCUJ9QtHNsx3k1iIjI5DCgtABCCFzJvqIbBlzV4nvDvO90cG0Oi+8REVHrxoDSTJWoSxCdEq27dVPV4nvlt24CFYEtfvE9IiJqWRhQmpHUglREJUUh6mYUjt86jiJVkW6fpdwSQR5BulE3nnYczktERM0XA0oZU5xJVq1R41zGOd2tmwuZF/T2u9u46wJJkGcQ17chIqIWgzPJVmDomWTVGlGn+UDylHk4mnwUUYlROJJ0BJnFmbp9EiT0dOupCyWdnTuzgysREZkMziTbTESeS6k0o6pnFTOqXs+5rutLUnHxPXtzewz2HoxQn1AM8R4CFyuuY0NERC0fA0ojKV+TpmLz1K2cYsz7IhqLx5mj0OwsohKjcCPvht4x7RzbIcQ7BKG+oejdpjfMZeZNVzgREZEJYEBpBGqNwGs/xOkvlifPg9zuIszsLsDM9jI2Xi7R7TOTmWGA+wCE+IRUWnyPiIioNWJAaQTR8Zm62zpymyuwbPML5NY39Y7RqOwwzDsEj3d7EAO9Bla7+B4REVFrxIDSCNLy7vQ5gZDrwom6yAeq/C5Q5XeDptgLY/r1xQg/byNVSUREZLoYUBpBG3sr3ffqorYoSn4M6oLOECqHao8jIiKiO2TGLqAlCmrnAk9HK2gHAMuhyumvF04kaEfzBLXjiBwiIqKqMKA0ArlMwoqxAQCAirOUlP+8YmxAjfOhEBERtWYMKI1kdA9PrJvSFx6O+rdxPBytsG5KX715UIiIiEgf+6A0otE9PPFAgEedZpIlIiIiBpRGJ5dJGNTB1dhlEBERNSu8xUNEREQmhwGFiIiITA4DSpmIiAgEBARgwIABxi6FiIio1ZOEEBXXs2vVDLlUNBERUWtiyM9QtqAQERGRyWFAISIiIpPDgEJEREQmhwGFiIiITA4DChEREZkcBhQiIiIyOZzqvoLyUde5ublGroSIiKh5Kf/sNMQMJgwoFeTl5QEAfH19jVwJERFR85SXlwdHR8cGXYMTtVWg0WiQnJwMe3t7SFLrWHU4NzcXvr6+uHnzJienM2F8nwyDr+MdLf21aAnPrzk9h/Ja4+Li0KVLF8hkDetFwhaUCmQyGXx8fIxdhlE4ODiY/P8AxPfJUPg63tHSX4uW8Pya03Pw9vZucDgB2EmWiIiITBADChEREZkcBhSCpaUlVqxYAUtLS2OXQjXg+2QYfB3vaOmvRUt4fs3pORi6VnaSJSIiIpPDFhQiIiIyOQwoREREZHIYUIiIiMjkMKAQERGRyWFAaUXWrVuHnj176ib8GTRoEH7++Wfd/uLiYoSHh8PV1RV2dnaYNGkSUlNTjVgxrVq1CpIkYdGiRbptfJ/ubeXKlZAkSe+ra9euuv2t7TVMSkrClClT4OrqCmtrawQGBuLkyZO6/UIILF++HJ6enrC2tsbIkSNx+fJlI1ZcN/7+/pXeb0mSEB4eDsD032+1Wo1XX30V7dq1g7W1NTp06IA33nhDbz0bU3uP8vLysGjRIvj5+cHa2hqDBw/GiRMnDFuvoFZj79694qeffhKXLl0SFy9eFP/85z+Fubm5OHfunBBCiLlz5wpfX1+xf/9+cfLkSTFw4EAxePBgI1fdekVHRwt/f3/Rs2dPsXDhQt12vk/3tmLFCtG9e3eRkpKi+7p9+7Zuf2t6DTMzM4Wfn5+YMWOGOH78uLh27Zr45ZdfxJUrV3THrFq1Sjg6Ooo9e/aI2NhYMW7cONGuXTtRVFRkxMprLy0tTe+9/vXXXwUAceDAASGE6b/fb775pnB1dRU//vijiI+PF998842ws7MTa9as0R1jau/RE088IQICAsShQ4fE5cuXxYoVK4SDg4NITEw0WL0MKK2cs7Oz+OSTT0R2drYwNzcX33zzjW7f+fPnBQBx7NgxI1bYOuXl5YlOnTqJX3/9VYSGhuoCCt+n2lmxYoXo1atXlfta22v48ssvi6FDh1a7X6PRCA8PD7F69WrdtuzsbGFpaSm2b9/eFCUa3MKFC0WHDh2ERqNpFu/3ww8/LGbNmqW3beLEieLpp58WQpjee1RYWCjkcrn48ccf9bb37dtX/Otf/zJYvbzF00qp1Wp8/fXXKCgowKBBg3Dq1CmUlpZi5MiRumO6du2Ktm3b4tixY0astHUKDw/Hww8/rPd+AOD7VAeXL1+Gl5cX2rdvj6effho3btwA0Ppew71796J///54/PHH0aZNG/Tp0webNm3S7Y+Pj8etW7f0Xg9HR0cEBwc3y9dDqVTiiy++wKxZsyBJUrN4vwcPHoz9+/fj0qVLAIDY2FgcOXIEY8aMAWB675FKpYJarYaVlZXedmtraxw5csRg9XKxwFbm7NmzGDRoEIqLi2FnZ4fdu3cjICAAMTExsLCwgJOTk97x7u7uuHXrlnGKbaW+/vprnD59Wu9+brlbt27xfaqF4OBgbNmyBV26dEFKSgpee+01DBs2DOfOnWt1r+G1a9ewbt06LF68GP/85z9x4sQJLFiwABYWFpg+fbruObu7u+ud11xfjz179iA7OxszZswA0Dz+n1m6dClyc3PRtWtXyOVyqNVqvPnmm3j66acBwOTeI3t7ewwaNAhvvPEGunXrBnd3d2zfvh3Hjh1Dx44dDVYvA0or06VLF8TExCAnJwfffvstpk+fjkOHDhm7LCpz8+ZNLFy4EL/++mulv06o9sr/8gSAnj17Ijg4GH5+fti5cyesra2NWFnT02g06N+/P9566y0AQJ8+fXDu3DmsX78e06dPN3J1hvfpp59izJgx8PLyMnYptbZz5058+eWX+Oqrr9C9e3fExMRg0aJF8PLyMtn3aNu2bZg1axa8vb0hl8vRt29fPPXUUzh16pTBHoO3eFoZCwsLdOzYEf369cPbb7+NXr16Yc2aNfDw8IBSqUR2drbe8ampqfDw8DBOsa3QqVOnkJaWhr59+8LMzAxmZmY4dOgQPvzwQ5iZmcHd3Z3vUz04OTmhc+fOuHLlSqv7Xff09ERAQIDetm7duulueZU/54qjWprj65GQkIDffvsNzz77rG5bc3i///GPf2Dp0qV48sknERgYiKlTp+KFF17A22+/DcA036MOHTrg0KFDyM/Px82bNxEdHY3S0lK0b9/eYPUyoLRyGo0GJSUl6NevH8zNzbF//37dvosXL+LGjRsYNGiQEStsXUaMGIGzZ88iJiZG99W/f388/fTTuu/5PtVdfn4+rl69Ck9Pz1b3uz5kyBBcvHhRb9ulS5fg5+cHAGjXrh08PDz0Xo/c3FwcP3682b0emzdvRps2bfDwww/rtjWH97uwsBAymf7HsVwuh0ajAWDa75GtrS08PT2RlZWFX375BePHjzdcvQbr1ksmb+nSpeLQoUMiPj5e/PXXX2Lp0qVCkiSxb98+IYR2KF7btm3F77//Lk6ePCkGDRokBg0aZOSq6e5RPELwfaqNF198URw8eFDEx8eLP/74Q4wcOVIoFAqRlpYmhGhdr2F0dLQwMzMTb775prh8+bL48ssvhY2Njfjiiy90x6xatUo4OTmJ77//Xvz1119i/PjxzWqYsRBCqNVq0bZtW/Hyyy9X2mfq7/f06dOFt7e3bpjxrl27hEKhEC+99JLuGFN7jyIjI8XPP/8srl27Jvbt2yd69eolgoODhVKpNFi9DCityKxZs4Sfn5+wsLAQbm5uYsSIEbpwIoQQRUVFIiwsTDg7OwsbGxvx6KOPipSUFCNWTEJUDih8n+5t8uTJwtPTU1hYWAhvb28xefJkvXk/Wttr+MMPP4gePXoIS0tL0bVrV7Fx40a9/RqNRrz66qvC3d1dWFpaihEjRoiLFy8aqdr6+eWXXwSAKus29fc7NzdXLFy4ULRt21ZYWVmJ9u3bi3/961+ipKREd4ypvUc7duwQ7du3FxYWFsLDw0OEh4eL7Oxsg9YrCXHXVHVEREREJoB9UIiIiMjkMKAQERGRyWFAISIiIpPDgEJEREQmhwGFiIiITA4DChEREZkcBhQiIiIyOQwoREREZHIYUIiIiMjkMKAQERGRyWFAIaqH++67D4sWLTLZ6zXWNanxGeN9y8jIQJs2bXD9+nWDXvfJJ5/Ee++9Z9BrUuvBgEKt2owZMyBJEiRJgoWFBTp27IjXX38dKpWqxvN27dqFN954w2B1GPp6tXXz5k3MmjULXl5esLCwgJ+fHxYuXIiMjIwmrwUwfqgq/31YtWqV3vY9e/ZAkiQjVdX43nzzTYwfPx7+/v4Gve4rr7yCN998Ezk5OQa9LrUODCjU6o0ePRopKSm4fPkyXnzxRaxcuRKrV6+u8lilUgkAcHFxgb29vcFqMPT1auPatWvo378/Ll++jO3bt+PKlStYv3499u/fj0GDBiEzM7NJ6zEVVlZWeOedd5CVlWXsUgym/Pe2KoWFhfj000/xzDPPGPxxe/TogQ4dOuCLL74w+LWp5WNAoVbP0tISHh4e8PPzw7x58zBy5Ejs3bsXgPYv+vnz52PRokVQKBQYNWqUbnv5X/r33XcfFixYgJdeegkuLi7w8PDAypUr9R5Do9Hg3XffRceOHWFpaYm2bdvizTff1O2v2HJQ/rjz58+Ho6MjFAoFXn31Vdy9+HhkZCSGDh0KJycnuLq64pFHHsHVq1dr/bzDw8NhYWGBffv2ITQ0FG3btsWYMWPw22+/ISkpCf/61790x/r7++ODDz7QO7937956z/Ne9dzrdZoxYwYOHTqENWvW6Fq1ym851Obx77vvPjz//PNYtGgRnJ2d4e7ujk2bNqGgoAAzZ86Evb09OnbsiJ9//rnG12XkyJHw8PDA22+/XeNxjVmTSqWq8b3XaDR4++230a5dO1hbW6NXr1749ttv9R63qt/bqvzvf/+DpaUlBg4cqNt29epVSJKEH3/8ESNGjICNjQ26dOmC48eP1/iaVGXs2LH4+uuv63weEQMKUQXW1tZ6f3Fu3boVFhYW+OOPP7B+/foqz9m6dStsbW1x/PhxvPvuu3j99dfx66+/6vYvW7YMq1atwquvvoq4uDh89dVXcHd3r7GOrVu3wszMDNHR0VizZg3ef/99fPLJJ7r9BQUFWLx4MU6ePIn9+/dDJpPh0UcfhUajuedzzMzMxC+//IKwsDBYW1vr7fPw8MDTTz+NHTt26H0o3ktt6qnpdVqzZg0GDRqE2bNnIyUlBSkpKfD19a3145dfX6FQIDo6Gs8//zzmzZuHxx9/HIMHD8bp06fx4IMPYurUqSgsLKz2GnK5HG+99RY++ugjJCYm1unxDVXTvd77t99+G59//jnWr1+Pv//+Gy+88AKmTJmCQ4cO6V3jXr+3AHD48GH069dPb1tsbCwkScL777+PV199FbGxsWjbti2WLl1a5+cfFBSE6OholJSU1PlcauUEUSs2ffp0MX78eCGEEBqNRvz666/C0tJSLFmyRAghRGhoqOjTp0+l80JDQ8XChQt13w8dOlRv/4ABA8TLL78shBAiNzdXWFpaik2bNlVbx93XK/+5W7duQqPR6La9/PLLolu3btVe4/bt2wKAOHv2bJXXvNuff/4pAIjdu3dXuf/9998XAERqaqoQQgg/Pz/x3//+V++YXr16iRUrVtSpnppep5pqrs3jV7y+SqUStra2YurUqbptKSkpAoA4duxYlTXf/fswcOBAMWvWLCGEELt37xYV/7lsrJru9d4XFxcLGxsbcfToUb3HfuaZZ8RTTz2lu0ZVv7dVGT9+vO55llu+fLlwdnYWaWlpum0ffvih6N69u95x5TWWP+e7ay4XGxsrAIjr16/Xqh6icmxBoVbvxx9/hJ2dHaysrDBmzBhMnjxZr5m+4l+XVenZs6fez56enkhLSwMAnD9/HiUlJRgxYkSd6ho4cKBex8xBgwbh8uXLUKvVAIDLly/jqaeeQvv27eHg4KDr4Hjjxo1aP4a4RwuJhYVFra9Vm3pqep0M4e7ry+VyuLq6IjAwULetvNWqNo/5zjvvYOvWrTh//nyT11TTe3/lyhUUFhbigQcegJ2dne7r888/17ulVpvfWwAoKiqClZWV3rbY2FiMHz8ebm5uum3x8fHo2LGj3nHr1q3Dxo0bUVBQgKVLlyIqKqrS9ctb6GpqtSKqipmxCyAytuHDh2PdunWwsLCAl5cXzMz0/7ewtbW95zXMzc31fpYkSXdro+ItFEMZO3Ys/Pz8sGnTJnh5eUGj0aBHjx41dogs17FjR0iShPPnz+PRRx+ttP/8+fNwc3ODk5MTAEAmk1UKM6WlpXWup6bXqSa1efzqrn/3tvIP/do8ZkhICEaNGoVly5ZhxowZJlETAOTn5wMAfvrpJ3h7e+vts7S01H1fm99bAFAoFJU6BMfGxmLZsmV622JiYhASEqK3LSwsDKtXr8aHH36I33//HUOGDKl0/fLO1neHHaLaYAsKtXq2trbo2LEj2rZtWymcGEKnTp1gbW2N/fv31+m8ih0S//zzT3Tq1AlyuRwZGRm4ePEiXnnlFYwYMQLdunWr06gTV1dXPPDAA1i7di2Kior09t26dQtffvml3oeym5sbUlJSdD/n5uYiPj5e93ND6ylnYWGhayG6270ev7GsWrUKP/zwA44dO9akNdX03gcEBMDS0hI3btxAx44d9b7q2mcHAPr06YO4uDjdzzk5Obh+/Tr69Omjd1xMTAx69+6tt239+vVwdHTEggUL8MMPP+Dw4cOVrn/u3Dn4+PhAoVDUuTZq3RhQiBqZlZUVXn75Zbz00ku6Zvg///wTn376aY3n3bhxA4sXL8bFixexfft2fPTRR1i4cCEAwNnZGa6urti4cSOuXLmC33//HYsXL65TXR9//DFKSkowatQoREVF4ebNm4iMjMQDDzyAzp07Y/ny5bpj77//fmzbtg2HDx/G2bNnMX36dMjlct1+Q9QDaEfGHD9+HNevX0d6erquVeFej99YAgMD8fTTT+PDDz+stK8xa6rpvbe3t8eSJUvwwgsvYOvWrbh69SpOnz6Njz76CFu3bq3zY40aNQp///23LlD+9ddfMDMz07sNlZCQgKysrEoBZc6cOXjuuedga2uLVatWYejQoZWuf/jwYTz44IN1rouIt3iImsCrr74KMzMzLF++HMnJyfD09MTcuXNrPGfatGkoKipCUFAQ5HI5Fi5ciOeeew6A9vbC119/jQULFqBHjx7o0qULPvzwQ9x33321rqlTp044ceIEVq5ciSeeeAJpaWkQQmDixInYtm0bbGxsdMcuW7YM8fHxeOSRR+Do6Ig33nhDr7XAEPUAwJIlSzB9+nQEBASgqKgI8fHx8Pf3v+fjN6bXX38dO3bsqLS9MWuq6b0HgDfeeANubm54++23ce3aNTg5OaFv37745z//WefHCgwMRN++fbFz507MmTMHsbGx6NKli16/lDNnzsDJyanSRG7lt6fK+2xVnMyuuLgYe/bsQWRkZJ3rIpLEvXrJEVGTu++++9C7d+9K82w0thUrVuD999/Hr7/+qjcvBrVsP/30E/7xj3/g3LlzkMkM17C+bt067N69G/v27TPYNan1YAsKEem89tpr8Pf3x59//omgoCCDfliR6Xr44Ydx+fJlJCUl1asfS3XMzc3x0UcfGex61LqwBYXIBBmrBYWIyFQwoBAREZHJYfstERERmRwGFCIiIjI5DChERERkchhQiIiIyOQwoBAREZHJYUAhIiIik8OAQkRERCaHAYWIiIhMDgMKERERmRwGFCIiIjI5/w9e0E0UwpKmMAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(6, 5))\n", "\n", "ax.plot(nu_s, lifetimes_s, \"o\", label=\"S states\")\n", "ax.plot(nu_circular, lifetimes_circular, \"o\", label=\"Circular states\")\n", "\n", "fit_s = np.exp(fit_function(np.log(nu_s), *popt_s))\n", "fit_circular = np.exp(fit_function(np.log(nu_circular), *popt_circular))\n", "ax.plot(nu_s, fit_s)\n", "ax.plot(nu_circular, fit_circular)\n", "\n", "ax.text(nu_s[2], fit_s[2], rf\"$\\propto n_*^{{{popt_s[0]:.2f}}}$\", verticalalignment=\"top\")\n", "ax.text(\n", " nu_circular[2],\n", " fit_circular[2],\n", " rf\"$\\propto n_*^{{{popt_circular[0]:.2f}}}$\",\n", " verticalalignment=\"top\",\n", ")\n", "\n", "ax.legend()\n", "ax.set_yscale(\"log\")\n", "ax.set_xscale(\"log\")\n", "ax.set_xlabel(\"Principal Quantum Number ($n_*$)\")\n", "ax.set_ylabel(r\"Lifetime ($\\mu$s)\")\n", "ax.set_title(\"Lifetime vs $n_*$\")\n", "\n", "ax.set_xticks([30, 40, 50, 60, 70, 80, 90])\n", "ax.get_xaxis().set_major_formatter(plt.ScalarFormatter())\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.13.1" } }, "nbformat": 4, "nbformat_minor": 1 }