{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Calculating Lifetimes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "from ryd_numerov.rydberg import RydbergState"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Lifetime at T=0: 25.50 microsecond\n",
      "Lifetime at T=300K: 15.80 microsecond\n"
     ]
    }
   ],
   "source": [
    "method = \"exact\"\n",
    "T_in_K = 300\n",
    "state = RydbergState(\"Rb\", n=30, l=0, j=0.5, m=0.5)\n",
    "\n",
    "lifetime_0 = state.get_lifetime(method=method)\n",
    "lifetime = state.get_lifetime(T_in_K, temperature_unit=\"K\", method=method)\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": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of possible spontaneous decay transitions: 125\n",
      "Number of considered BBR transitions: 280\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlUAAAHNCAYAAAA3244BAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAb7BJREFUeJzt3XlcVFX/B/DPsAz7ooAMKAIpuQsKgribJJobbqGZIm4t7mSmpeBSYW65JumTopVJ+CiPlaJIZouIC4uiaGoopiyukKigcH5/+OPmyIAwXAT083695lVzzrnnfu+ZGebruXfOVQghBIiIiIioUnSqOwAiIiKi5wGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiKqFebNmweFQlHdYbywwsPDoVAocPHiRVnbEj1PmFRRtTp58iSGDBkCR0dHGBoaon79+nj11VexevXq6g4NAHDo0CHMmzcPt2/fru5QqJyKv9Aff9SrVw/du3fHnj17tG6rqb2enh7q16+P0aNH48qVK0+N7cl9lfb45Zdf5BqOKlUTPx+VfY1KUxOPlWoeveoOgF5chw4dQvfu3dGwYUOMHz8eKpUKly9fxuHDh7Fy5UpMnjy5ukPEoUOHMH/+fIwePRqWlpbVHQ5VwIIFC+Ds7AwhBLKyshAeHo7XXnsNP/zwA/r27at128fb379/H4cPH0Z4eDh+//13pKSkwNDQsNSYvv76a7XnW7ZsQUxMTInyZs2aVeLIq8bIkSMxbNgwGBgYSGWlfT40tX3WtH2NSsO/BVQeTKqo2nzyySewsLDA0aNHS/yRys7Orp6g6LnRu3dveHh4SM/Hjh0LW1tbfPfddyUSpYq0fbL9uHHjYG1tjc8++wy7du3C66+/XmpMb775ptrzw4cPIyYmpkR5afLy8mBiYlKutnLT1dWFrq6u7G2riravEVFl8PQfVZsLFy6gRYsWGv/VV69ePen/i6+lOXPmDF5//XWYm5vDysoKU6dOxf3790tsm5iYiN69e8Pc3Bympqbo0aMHDh8+rNamuM/z589L//K0sLBAYGAg7t69K7V5//33AQDOzs7S6YTi60QuXbqEd999F02aNIGRkRGsrKwwdOjQEteRlGdfj7ty5QrGjBkDW1tbGBgYoEWLFti4caNWxzl69Gg4OTmV2PbJ65P++ecfTJs2DU5OTjAwMEC9evXw6quvIiEhocS2jyvvGFR0HH7//Xe0a9cOhoaGaNSoEb788ssy4ygPS0tLGBkZQU/v6f+WrEhbAOjcuTOAR+9puRSP1+nTp/HGG2+gTp066NSpEwD533vlef2fvE6qrM9HaddUyfXZ1Iam16gi41jW34LyfGa1+Yxt374dCoUCBw8eLFH35ZdfQqFQICUlRev+i3Xv3h1dunRBQkICevfuDTMzM9SvXx8rV6586rakjjNVVG0cHR0RFxeHlJQUtGzZ8qntX3/9dTg5OSE0NBSHDx/GqlWrcOvWLWzZskVqc+rUKXTu3Bnm5uaYOXMm9PX18eWXX6Jbt244ePAgvLy8SvTp7OyM0NBQJCQk4D//+Q/q1auHzz77DIMGDcKff/6J7777Dp9//jmsra0BADY2NgCAo0eP4tChQxg2bBgaNGiAixcvYt26dejWrRtOnz4NY2Pjcu+rWFZWFtq3bw+FQoFJkybBxsYGe/bswdixY5Gbm4tp06ZpdZxP8/bbb2P79u2YNGkSmjdvjhs3buD3339Hamoq2rZtW+p2FR2D8ozDyZMn0bNnT9jY2GDevHl4+PAhQkJCYGtrW6FjysnJwfXr1yGEQHZ2NlavXo07d+5onBWqSFtNir9c69SpU6EYy2Po0KFwcXHBp59+CiEEAPnfe9q8/k/7fDxJzs+mNjS9RuUdx7KOtbyfWW3GuE+fPjA1NcX333+Prl27qtVFRESgRYsW0t9ObT/DwKPPnL29Pfr164fAwED4+flhw4YNmD59Ol555RW0atWqwuP9whJE1WTfvn1CV1dX6OrqCm9vbzFz5kyxd+9eUVBQoNYuJCREABD9+/dXK3/33XcFAJGcnCyV+fn5CaVSKS5cuCCVXb16VZiZmYkuXbqU6HPMmDFqfQ4cOFBYWVlJz5csWSIAiLS0tBLx3717t0RZXFycACC2bNlS4X0JIcTYsWOFnZ2duH79ulr5sGHDhIWFhbTP8h5nQECAcHR0LBFncUzFLCwsxMSJE0u0e5ryjsHj+3zaOPj5+QlDQ0Nx6dIlqez06dNCV1dXlOdP1qZNmwSAEg8DAwMRHh6uddvH2+/fv19cu3ZNXL58WWzfvl3Y2NgIAwMDcfny5afG97iJEyeWekzF4zV8+PASdXK/98rz+hcf++OfhdI+H5rayv3ZfFqc5XmNKvL+Le1Yy/uZ1fYzNnz4cFGvXj3x8OFDqSwjI0Po6OiIBQsWSGXa9n/16lUBQNjY2KiNzenTpwUAsXnz5gr3+SLj6T+qNq+++iri4uLQv39/JCcnY/HixfD19UX9+vWxa9euEu0nTpyo9rz4Qvbdu3cDAAoLC7Fv3z74+fnhpZdektrZ2dnhjTfewO+//47c3Fy1Pt5++2215507d8aNGzdKtNPEyMhI+v8HDx7gxo0baNy4MSwtLTVOuT9tX0II/Pe//0W/fv0ghMD169elh6+vL3JycpCQkKDVcT6NpaUl4uPjcfXq1QptV9ExAMoeh8LCQuzduxd+fn5o2LCh1KZZs2bw9fWtUGxr165FTEwMYmJi8M0336B79+4YN24cduzYUam2AODj4wMbGxs4ODhgyJAhMDExwa5du9CgQYMKxVgeT44XIP97T9vXv7ye9WcTKN9rpM3793Hl/cwC2o+xv78/srOz1X4Run37dhQVFcHf318q07b/kydPAgBCQkLUxkZfXx8AoFQqK9Tfi45JFVWrdu3aYceOHbh16xaOHDmC2bNn459//sGQIUNw+vRptbYuLi5qzxs1agQdHR1pWv/atWu4e/cumjRpUmI/zZo1Q1FRES5fvqxW/vgXN/DvqYFbt249NfZ79+4hODgYDg4OMDAwgLW1NWxsbHD79m3k5OSUaP+0fV27dg23b9/G+vXrYWNjo/YIDAwE8OgCfm2O82kWL16MlJQUODg4wNPTE/PmzcNff/0l+xg8bRyuXbuGe/fulXitAWg83rJ4enrCx8cHPj4+GDFiBH766Sc0b94ckyZNQkFBgdZtgX+TsO3bt+O1117D9evXq+yXbs7OziXK5H7vafv6l9ez/mwC5XuNtHn/Pnlc5fnMAtqPca9evWBhYYGIiAipLCIiAm5ubnj55ZelMm37L06q/Pz81MrPnDkDoOKfuxcdkyqqEZRKJdq1a4dPP/0U69atw4MHDxAZGVnmNnIsBFnaL5TE/1+7UpbJkyfjk08+weuvv47vv/8e+/btQ0xMDKysrFBUVFThfRVv8+abb0qzJk8+OnbsWN5DA1D6GBUWFqo9f/311/HXX39h9erVsLe3x5IlS9CiRQuNazU9rqJjAFRuzCtDR0cH3bt3R0ZGBs6dO1eptsVJ2ODBg7Fr1y60bNkSb7zxBu7cuSN73I/PphST+72n7etflSr7PinPa6TN+/dxFfnMajvGBgYG8PPzw86dO/Hw4UNcuXIFf/zxh9osVWX6P3HiBFQqFerXr69WnpycDD09PTRv3vyp40D/4oXqVOMU/ww6IyNDrfzcuXNq/2o/f/48ioqKpF+32djYwNjYGGfPni3R55kzZ6CjowMHB4cKxVJW4rZ9+3YEBARg2bJlUtn9+/e1XhzQxsYGZmZmKCwshI+PT6ntCgsLy32cderU0RjPpUuXSpTZ2dnh3Xffxbvvvovs7Gy0bdsWn3zyCXr37l1qLFUxBkZGRhoTGU3HW1EPHz4EgHIlP+Vtq6uri9DQUHTv3h1r1qzBrFmzKh3n08g97oB2r395/2FTFZ/NiijtNarIOGo61vJ+ZotpM8bAo1OAmzdvRmxsLFJTUyGEKJFUadv/yZMn4erqWqL8xIkTePnll6t1rbHaiDNVVG0OHDig8V+dxddIPTntvHbtWrXnxauuF//B0NXVRc+ePfG///1P7SfRWVlZ2Lp1Kzp16gRzc/MKxVi8JpCmP7K6urol4l+9enWJWaDy0tXVxeDBg/Hf//5X+pn0465duya1K+9xNmrUCDk5OThx4oTULiMjAzt37pSeFxYWljjVUa9ePdjb2yM/P/+pMcs9Br6+voiKikJ6erpUnpqair1792rVZ7EHDx5g3759UCqVT11csyJtAaBbt27w9PTEihUrNC7zITc5x70yr39Zn48n45X7s1lRml6jioyjpmMt72e2MmMMPLo+rG7duoiIiEBERAQ8PT3V/oGpbf+FhYVITU3VmFQlJyejdevWT42N1HGmiqrN5MmTcffuXQwcOBBNmzZFQUEBDh06hIiICDg5OUnXJBRLS0tD//790atXL8TFxeGbb77BG2+8ofYH4eOPP0ZMTAw6deqEd999F3p6evjyyy+Rn5+PxYsXVzhGd3d3AMBHH32EYcOGQV9fH/369YOJiQn69u2Lr7/+GhYWFmjevDni4uKwf/9+WFlZaT0mixYtwoEDB+Dl5YXx48ejefPmuHnzJhISErB//37cvHmzQsc5bNgwfPDBBxg4cCCmTJmCu3fvYt26dXj55ZelC2j/+ecfNGjQAEOGDIGrqytMTU2xf/9+HD16VO1f8JpUxRjMnz8f0dHR6Ny5M9599108fPgQq1evRosWLdSSw6fZs2ePdF1IdnY2tm7dinPnzmHWrFklvsAr0rY077//PoYOHYrw8HCNF5fLSc5xr8zrX9rnQxO5P5vaePI1qsg4lnas5fnMVmaMgUcXjQ8aNAjbtm1DXl4eli5dqlavbf/nzp3D/fv3SyRV9+7dw/nz5xEQEFCeYaXHVcMvDomEEELs2bNHjBkzRjRt2lSYmpoKpVIpGjduLCZPniyysrKkdsU/sT59+rQYMmSIMDMzE3Xq1BGTJk0S9+7dK9FvQkKC8PX1FaampsLY2Fh0795dHDp0SK1NcZ/Xrl1TK9f0U/CFCxeK+vXrCx0dHbW6W7duicDAQGFtbS1MTU2Fr6+vOHPmjHB0dBQBAQFa7UsIIbKyssTEiROFg4OD0NfXFyqVSvTo0UOsX7++wscpxKOlK1q2bCmUSqVo0qSJ+Oabb9SWVMjPzxfvv/++cHV1FWZmZsLExES4urqKL774okRfTyrvGFR0HA4ePCjc3d2FUqkUL730kggLCyuxDERpNC2TYGhoKNzc3MS6detEUVGRVm0fb3/06NES+y0sLBSNGjUSjRo1Uvv5e1nKs6TCk+MlhLzvvfK+/qW9XzV9PkprK/dnU5OKvEYVef+WdqxCPP0zW5nPWLGYmBgBQCgUihJLd2jb//fffy8AiJSUFLXyI0eOCADixx9/LHd89IhCiCq+OpSokubNm4f58+fj2rVr0qJ7RERENQ2vqSIiIiKSAZMqIiIiIhkwqSIiIiKSAa+pIiIiIpIBZ6qIiIiIZMCkioiIiEgGXPzzGSoqKsLVq1dhZmYmy33riIiIqOoJIfDPP//A3t4eOjqlz0cxqXqGrl69WqX3tyIiIqKqc/nyZTRo0KDUeiZVz5CZmRmARy9KVd/nioiIiOSRm5sLBwcH6Xu8NEyqnqHiU37m5uZMqoiIiGqZp126wwvViYiIiGTApIqIiIhIBkyqiIiIiGTAa6qIiF5QRUVFKCgoqO4wiKqdvr4+dHV1K90PkyoiohdQQUEB0tLSUFRUVN2hENUIlpaWUKlUlVpHkkkVEdELRgiBjIwM6OrqwsHBoczFDImed0II3L17F9nZ2QAAOzs7rftiUkVE9IJ5+PAh7t69C3t7exgbG1d3OETVzsjICACQnZ2NevXqaX0qkP88ISJ6wRQWFgIAlEplNUdCVHMU/wPjwYMHWvfBpIqI6AXFe5AS/UuOzwOTKiIiIiIZMKkiIiIikgEvVCciIgCAh8ez3d+xYxVrf+3aNQQHB+Onn35CVlYW6tSpA1dXVwQHB6Njx45VE2QpnJycMG3aNEybNu2Z7re6jB49Gps3bwYA6OnpoW7dumjdujWGDx+O0aNH8xek/49JFRER1QqDBw9GQUEBNm/ejJdeeglZWVmIjY3FjRs3qju0F0KvXr2wadMmFBYWIisrC9HR0Zg6dSq2b9+OXbt2QU+PKUWNSy3Xrl0LJycnGBoawsvLC0eOHCmzfWRkJJo2bQpDQ0O0atUKu3fvVqvfsWMHevbsCSsrKygUCiQlJanVX7x4EQqFQuMjMjJSaqepftu2bbIdNxERle727dv47bff8Nlnn6F79+5wdHSEp6cnZs+ejf79+0vtFAoF1q1bh969e8PIyAgvvfQStm/frtbXyZMn8corr8DIyAhWVlaYMGEC7ty5I9WPHj0afn5+WLp0Kezs7GBlZYWJEydKvwrr1q0bLl26hOnTp0vfBwBw48YNDB8+HPXr14exsTFatWqF7777Tm3f3bp1w5QpUzBz5kzUrVsXKpUK8+bNK3Gs48aNg42NDczNzfHKK68gOTlZrc26devQqFEjKJVKNGnSBF9//bVUV/y99vj33e3bt6FQKPDLL78AAG7duoURI0bAxsYGRkZGcHFxwaZNm8p8DQwMDKBSqVC/fn20bdsWH374If73v/9hz549CA8Pr1D8P/zwA9q1awdDQ0NYW1tj4MCBUt3XX38NDw8PmJmZQaVS4Y033pDWkBJCoHHjxli6dKlaf0lJSVAoFDh//nyZx1DValRSFRERgaCgIISEhCAhIQGurq7w9fWVBvNJhw4dwvDhwzF27FgkJibCz88Pfn5+SElJkdrk5eWhU6dO+OyzzzT24eDggIyMDLXH/PnzYWpqit69e6u13bRpk1o7Pz8/2Y6diIhKZ2pqClNTU0RFRSE/P7/MtnPnzsXgwYORnJyMESNGYNiwYUhNTQXw6DvB19cXderUwdGjRxEZGYn9+/dj0qRJan0cOHAAFy5cwIEDB7B582aEh4dLicOOHTvQoEEDLFiwQPo+AID79+/D3d0dP/30E1JSUjBhwgSMHDmyxOTA5s2bYWJigvj4eCxevBgLFixATEyMVD906FBkZ2djz549OH78ONq2bYsePXrg5s2bAICdO3di6tSpeO+995CSkoK33noLgYGBOHDgQLnHc+7cuTh9+jT27NmD1NRUrFu3DtbW1uXevtgrr7wCV1dX7Nixo9zx//TTTxg4cCBee+01JCYmIjY2Fp6entL2Dx48wMKFC5GcnIyoqChcvHgRo0ePBvAoaR4zZkyJBHDTpk3o0qULGjduXOFjkJWoQTw9PcXEiROl54WFhcLe3l6EhoZqbP/666+LPn36qJV5eXmJt956q0TbtLQ0AUAkJiY+NQ43NzcxZswYtTIAYufOnU8/iDLk5OQIACInJ6dS/dQW7u6aH0RUve7duydOnz4t7t27p1Ze2me2qh4VtX37dlGnTh1haGgoOnToIGbPni2Sk5PV2gAQb7/9tlqZl5eXeOedd4QQQqxfv17UqVNH3LlzR6r/6aefhI6OjsjMzBRCCBEQECAcHR3Fw4cPpTZDhw4V/v7+0nNHR0fx+eefPzXmPn36iPfee0963rVrV9GpUye1Nu3atRMffPCBEEKI3377TZibm4v79++rtWnUqJH48ssvhRBCdOjQQYwfP16tfujQoeK1114TQmj+vrt165YAIA4cOCCEEKJfv34iMDDwqfEXCwgIEAMGDNBY5+/vL5o1a1bu+L29vcWIESPKve+jR48KAOKff/4RQghx5coVoaurK+Lj44UQQhQUFAhra2sRHh5e7j41Ke1zIUT5v79rzExVQUEBjh8/Dh8fH6lMR0cHPj4+iIuL07hNXFycWnsA8PX1LbV9eRw/fhxJSUkYO3ZsibqJEyfC2toanp6e2LhxI4QQZfaVn5+P3NxctQcREWln8ODBuHr1Knbt2oVevXrhl19+Qdu2bdVOPQGAt7d3iefFM1WpqalwdXWFiYmJVN+xY0cUFRXh7NmzUlmLFi3UVtW2s7Mr9axJscLCQixcuBCtWrVC3bp1YWpqir179yI9PV2tXevWrdWeP953cnIy7ty5AysrK2l2ztTUFGlpabhw4YJ0DE9emN+xY0fpGMvjnXfewbZt2+Dm5oaZM2fi0KFD5d72SUII6RRoeeJPSkpCjx49Su3v+PHj6NevHxo2bAgzMzN07doVAKRxtLe3R58+fbBx40YAj04l5ufnY+jQoVofg1xqzFVl169fR2FhIWxtbdXKbW1tcebMGY3bZGZmamyfmZmpdRxfffUVmjVrhg4dOqiVL1iwAK+88gqMjY2xb98+vPvuu7hz5w6mTJlSal+hoaGYP3++1rEQEZE6Q0NDvPrqq3j11Vcxd+5cjBs3DiEhIdLpIbno6+urPVcoFE+9+fSSJUuwcuVKrFixAq1atYKJiQmmTZuGgoKCcvd9584d2NnZSdc+Pc7S0rJcsRf/Eu/xf/g/uUp47969cenSJezevRsxMTHo0aMHJk6cWOJapfJITU2Fs7NzueMvviWMJsWnZ319ffHtt9/CxsYG6enp8PX1VRvHcePGYeTIkfj888+xadMm+Pv714hbLtWYmaqa4N69e9i6davGWaq5c+eiY8eOaNOmDT744APMnDkTS5YsKbO/2bNnIycnR3pcvny5qkInInohNW/eHHl5eWplhw8fLvG8WbNmAIBmzZohOTlZbZs//vgDOjo6aNKkSbn3q1Qqpdv9PN7PgAED8Oabb8LV1RUvvfQS/vzzzwodT9u2bZGZmQk9PT00btxY7VF8zVOzZs3wxx9/lNh38+bNAQA2NjYAIF3rBaDEj7SK2wUEBOCbb77BihUrsH79+grFCgA///wzTp48icGDB5c7/tatWyM2NlZjf2fOnMGNGzewaNEidO7cGU2bNtU4Q/jaa6/BxMQE69atQ3R0NMaMGVPh2KtCjUmqrK2toauri6ysLLXyrKwsqFQqjduoVKoKtX+a7du34+7duxg1atRT23p5eeHvv/8u84JJAwMDmJubqz2IiKjibty4gVdeeQXffPMNTpw4gbS0NERGRmLx4sUYMGCAWtvIyEhs3LgRf/75J0JCQnDkyBHpQvQRI0bA0NAQAQEBSElJwYEDBzB58mSMHDmyxJmPsjg5OeHXX3/FlStXcP36dQCAi4sLYmJicOjQIaSmpuKtt94q8R31ND4+PvD29oafnx/27duHixcv4tChQ/joo49w7P8X9nr//fcRHh6OdevW4dy5c1i+fDl27NiBGTNmAHg0E9S+fXssWrQIqampOHjwIObMmaO2n+DgYPzvf//D+fPncerUKfz4449S4lma/Px8ZGZm4sqVK0hISMCnn36KAQMGoG/fvtL3ZnniDwkJwXfffYeQkBCkpqbi5MmT0o/JGjZsCKVSidWrV+Ovv/7Crl27sHDhwhKx6OrqYvTo0Zg9ezZcXFxKnPKtLjUmqVIqlXB3d1fLXouKihAbG1vqYHl7e5fIdmNiYrQe3K+++gr9+/eXsvyyJCUloU6dOjAwMNBqX0REVH6mpqbw8vLC559/ji5duqBly5aYO3cuxo8fjzVr1qi1nT9/PrZt24bWrVtjy5Yt+O6776RZHGNjY+zduxc3b95Eu3btMGTIEPTo0aNEH0+zYMECXLx4EY0aNZK+M+bMmYO2bdvC19cX3bp1g0qlqvCvxBUKBXbv3o0uXbogMDAQL7/8MoYNG4ZLly5JSZ+fnx9WrlyJpUuXokWLFvjyyy+xadMmdOvWTepn48aNePjwIdzd3TFt2jR8/PHHavtRKpWYPXs2WrdujS5dukBXV/epywRFR0fDzs4OTk5O6NWrFw4cOIBVq1bhf//7n3T9WXni79atGyIjI7Fr1y64ubnhlVdekX4haWNjg/DwcERGRqJ58+ZYtGhRqackx44di4KCAgQGBlZojKtUpS6Vl9m2bduEgYGBCA8PF6dPnxYTJkwQlpaW0i8yRo4cKWbNmiW1/+OPP4Senp5YunSpSE1NFSEhIUJfX1+cPHlSanPjxg2RmJgofvrpJwFAbNu2TSQmJoqMjAy1fZ87d04oFAqxZ8+eEnHt2rVLbNiwQZw8eVKcO3dOfPHFF8LY2FgEBwdX6Pj46z/++o+oJijrV061HWT4pTbVDr/++qvQ19eXcoTKkuPXfzXmQnUA8Pf3l25DkJmZCTc3N0RHR0vZbXp6utpS+B06dMDWrVsxZ84cfPjhh3BxcUFUVBRatmwptdm1a5daFjts2DAAj6YfH19wbePGjWjQoAF69uxZIi59fX2sXbsW06dPlxYeW758OcaPHy/3EBAREVEZ8vPzce3aNcybNw9Dhw6t0GnbqqYQ4inrApBscnNzYWFhgZycnBfi+qrS7iNW0ft9EZG87t+/j7S0NDg7O8PQ0LC6w5GVQqHAzp07uTjzcyw8PBxjx46Fm5sbdu3ahfr168vSb1mfi/J+f9eomSoiIqLK4DzB82/06NGyL6EhlxpzoToRERFRbcakioiIiEgGTKqIiIiIZMCkioiIiEgGTKqIiIiIZMCkioiIiEgGTKqIiOi5cfHiRSgUCo03ENaWQqFAVFSUbP1VRnh4OCwtLas7jFI5OTlhxYoV0nO5xq4mvQZl4TpVRET0SGkr9laVCq4EPHr0aGzevFl6XrduXbRr1w6LFy9G69at5Y5Oa+Hh4Wp38jAxMUGTJk3w0UcfYdCgQdUY2SPdunXDwYMHAQAGBgZo2LAhAgMDMWvWLCgUCln3lZGRgTp16pS7/bx58xAVFVUiKa5oP9WFM1VERFRr9OrVCxkZGcjIyEBsbCz09PTQt2/f6g6rBHNzcynOxMRE+Pr64vXXX8fZs2erOzQAwPjx45GRkYGzZ89i9uzZCA4ORlhYmOz7UalUMDAwqDH9VDUmVUREVGsYGBhApVJBpVLBzc0Ns2bNwuXLl3Ht2jWN7QsLCzF27Fg4OzvDyMgITZo0wcqVK0u027hxI1q0aAEDAwPY2dlh0qRJpcYQEhICOzs7nDhxotQ2CoVCitPFxQUff/wxdHR01La5desWRo0ahTp16sDY2Bi9e/fGuXPn1PoJDw9Hw4YNYWxsjIEDB+LGjRtS3cWLF6Gjo4NjT8z4rVixAo6OjigqKio1PmNjY6hUKjg6OiIwMBCtW7dGTEyMVH/hwgUMGDAAtra2MDU1Rbt27bB//361PrKzs9GvXz8YGRnB2dkZ3377rcZxePy03QcffICXX34ZxsbGeOmllzB37lw8ePBAOtb58+cjOTkZCoUCCoUC4eHhGvs5efIkXnnlFRgZGcHKygoTJkzAnTt3pPrRo0fDz88PS5cuhZ2dHaysrDBx4kRpX1WFSRUREdVKd+7cwTfffIPGjRvDyspKY5uioiI0aNAAkZGROH36NIKDg/Hhhx/i+++/l9qsW7cOEydOxIQJE3Dy5Ens2rULjRs3LtGXEAKTJ0/Gli1b8Ntvv5X7lGNhYaF02rJt27ZS+ejRo3Hs2DHs2rULcXFxEELgtddek7744+PjMXbsWEyaNAlJSUno3r07Pv74Y2l7Jycn+Pj4YNOmTWr727RpE0aPHg0dnad/xQsh8Ntvv+HMmTNQKpVS+Z07d/Daa68hNjYWiYmJ6NWrF/r164f09HS1+C9fvowDBw5g+/bt+OKLL5CdnV3m/szMzBAeHo7Tp09j5cqV2LBhAz7//HMAgL+/P9577z20aNFCmuXz9/cv0UdeXh58fX1Rp04dHD16FJGRkdi/f3+JRPjAgQO4cOECDhw4gM2bNyM8PFxK0qqMoGcmJydHABA5OTnVHcoz4e6u+UFE1evevXvi9OnT4t69e+oVpX1oq+pRQQEBAUJXV1eYmJgIExMTAUDY2dmJ48ePS23S0tIEAJGYmFhqPxMnThSDBw+Wntvb24uPPvqo1PYARGRkpHjjjTdEs2bNxN9//11mnJs2bRIApDh1dHSEgYGB2LRpk9Tmzz//FADEH3/8IZVdv35dGBkZie+//14IIcTw4cPFa6+9pta3v7+/sLCwkJ5HRESIOnXqiPv37wshhDh+/LhQKBQiLS2t1Pi6du0q9PX1hYmJidDX1xcAhKGhoVosmrRo0UKsXr1aCCHE2bNnBQBx5MgRqT41NVUAEJ9//rlUBkDs3Lmz1D6XLFki3B97L4SEhAhXV9cS7R7vZ/369aJOnTrizp07Uv1PP/0kdHR0RGZmphDi0XvF0dFRPHz4UGozdOhQ4e/vX2ospX4uRPm/vzlTRUREtUb37t2RlJSEpKQkHDlyBL6+vujduzcuXbpU6jZr166Fu7s7bGxsYGpqivXr10szLtnZ2bh69Sp69OhR5n6nT5+O+Ph4/Prrr6hfv/5T4zQzM5PiTExMxKeffoq3334bP/zwAwAgNTUVenp68PLykraxsrJCkyZNkJqaKrV5vB4AvL291Z77+flBV1cXO3fuBPDoFFr37t3h5ORUZnwjRoxAUlIS/vjjD/Tu3RsfffQROnToINXfuXMHM2bMQLNmzWBpaQlTU1OkpqZK41Ycv7u7u7RN06ZNn/rLxIiICHTs2BEqlQqmpqaYM2eO2uxXeaSmpsLV1RUmJiZSWceOHVFUVKR2zVqLFi2gq6srPbezs3vqTFplMakiIqJaw8TEBI0bN0bjxo3Rrl07/Oc//0FeXh42bNigsf22bdswY8YMjB07Fvv27UNSUhICAwNRUFAAADAyMirXfl999VVcuXIFe/fuLVd7HR0dKc7WrVsjKCgI3bp1w2effVa+Ay0npVKJUaNGYdOmTSgoKMDWrVsxZsyYp25nYWEhjeH333+PNWvWqF0zNWPGDOzcuROffvopfvvtNyQlJaFVq1bSuGkjLi4OI0aMwGuvvYYff/wRiYmJ+OijjyrVZ1n09fXVnisUijKvM5MDkyoiIqq1FAoFdHR0cO/ePY31f/zxBzp06IB3330Xbdq0QePGjXHhwgWp3szMDE5OToiNjS1zP/3798fWrVsxbtw4bNu2TatYdXV1pTibNWuGhw8fIj4+Xqq/ceMGzp49i+bNm0ttHq8HgMOHD5fod9y4cdi/fz+++OILPHz4sMLLNpiammLq1KmYMWMGHp1pezRuo0ePxsCBA9GqVSuoVCpcvHhR2qZp06Z4+PAhjh8/LpWdPXsWt2/fLnU/hw4dgqOjIz766CN4eHjAxcWlxAyjUqlEYWFhmfE2a9YMycnJyMvLk8r++OMP6OjooEmTJhU4cvkxqSIiolojPz8fmZmZyMzMRGpqKiZPnow7d+6gX79+Gtu7uLjg2LFj2Lt3L/7880/MnTsXR48eVWszb948LFu2DKtWrcK5c+eQkJCA1atXl+hr4MCB+PrrrxEYGIjt27eXGacQQoozLS0N69evx969ezFgwAAprgEDBmD8+PH4/fffkZycjDfffBP169eX2kyZMgXR0dFYunQpzp07hzVr1iA6OrrEvpo1a4b27dvjgw8+wPDhw8s9+/a4t956C3/++Sf++9//SvHt2LEDSUlJSE5OxhtvvKE2y9OkSRP06tULb731FuLj43H8+HGMGzeuzH27uLggPT0d27Ztw4ULF7Bq1SrptGUxJycnpKWlISkpCdevX0d+fn6JfkaMGAFDQ0MEBAQgJSUFBw4cwOTJkzFy5EjY2tpW+NjlxKSKiIhqjejoaNjZ2cHOzg5eXl7Sr7+6deumsf1bb72FQYMGwd/fH15eXrhx4wbeffddtTYBAQFYsWIFvvjiC7Ro0QJ9+/YtsbRBsSFDhmDz5s0YOXIkduzYUWqcubm5UpzNmjXDsmXLsGDBAnz00UdSm02bNsHd3R19+/aFt7c3hBDYvXu3dNqqffv22LBhA1auXAlXV1fs27cPc+bM0bi/sWPHoqCgoFyn/jSpW7cuRo0ahXnz5qGoqAjLly9HnTp10KFDB/Tr1w++vr5qv1wsjt/e3h5du3bFoEGDMGHCBNSrV6/UffTv3x/Tp0/HpEmT4ObmhkOHDmHu3LlqbQYPHoxevXqhe/fusLGxwXfffVeiH2NjY+zduxc3b95Eu3btMGTIEPTo0QNr1qzR6tjlpBDFc31U5XJzc2FhYYGcnByYm5tXdzhVrrTFmSu4iDIRyez+/ftIS0uDs7MzDA0NqzscksHChQsRGRlZ5tpZVLayPhfl/f7mTBUREVEtdefOHaSkpGDNmjWYPHlydYfzwmNSRUREVEtNmjQJ7u7u6Natm9an/kg+vKEyERFRLfVMVgmncuNMFREREZEMmFQREb2g+Dslon/J8XlgUkVE9IIpvnVHVa1kTVQb3b17F0DJldgrgtdUERG9YPT09GBsbIxr165BX18fOjr89zW9uIQQuHv3LrKzs2Fpaal2v8CKYlJFRPSCUSgUsLOzQ1paWpk3IiZ6kVhaWkKlUlWqDyZVREQvIKVSCRcXF54CJMKjU36VmaEqxqSKiOgFpaOjwxXViWTEE+lEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCSDGpdUrV27Fk5OTjA0NISXlxeOHDlSZvvIyEg0bdoUhoaGaNWqFXbv3q1Wv2PHDvTs2RNWVlZQKBRISkoq0Ue3bt2gUCjUHm+//bZam/T0dPTp0wfGxsaoV68e3n//fTx8+LDSx0tERETPhxqVVEVERCAoKAghISFISEiAq6srfH19kZ2drbH9oUOHMHz4cIwdOxaJiYnw8/ODn58fUlJSpDZ5eXno1KkTPvvsszL3PX78eGRkZEiPxYsXS3WFhYXo06cPCgoKcOjQIWzevBnh4eEIDg6W58CJiIio1lMIIUR1B1HMy8sL7dq1w5o1awAARUVFcHBwwOTJkzFr1qwS7f39/ZGXl4cff/xRKmvfvj3c3NwQFham1vbixYtwdnZGYmIi3Nzc1Oq6desGNzc3rFixQmNce/bsQd++fXH16lXY2toCAMLCwvDBBx/g2rVrUCqV5Tq+3NxcWFhYICcnB+bm5uXapjbz8NBcfuzYs42DiIioMsr7/V1jZqoKCgpw/Phx+Pj4SGU6Ojrw8fFBXFycxm3i4uLU2gOAr69vqe3L8u2338La2hotW7bE7NmzcffuXbX9tGrVSkqoiveTm5uLU6dOldpnfn4+cnNz1R5ERET0fNKr7gCKXb9+HYWFhWqJCwDY2trizJkzGrfJzMzU2D4zM7NC+37jjTfg6OgIe3t7nDhxAh988AHOnj2LHTt2lLmf4rrShIaGYv78+RWKhYiIiGqnGpNUVacJEyZI/9+qVSvY2dmhR48euHDhAho1aqR1v7Nnz0ZQUJD0PDc3Fw4ODpWKlYiIiGqmGnP6z9raGrq6usjKylIrz8rKgkql0riNSqWqUPvy8vLyAgCcP3++zP0U15XGwMAA5ubmag8iIiJ6PtWYpEqpVMLd3R2xsbFSWVFREWJjY+Ht7a1xG29vb7X2ABATE1Nq+/IqXnbBzs5O2s/JkyfVfoUYExMDc3NzNG/evFL7IiIioudDjTr9FxQUhICAAHh4eMDT0xMrVqxAXl4eAgMDAQCjRo1C/fr1ERoaCgCYOnUqunbtimXLlqFPnz7Ytm0bjh07hvXr10t93rx5E+np6bh69SoA4OzZswAezTCpVCpcuHABW7duxWuvvQYrKyucOHEC06dPR5cuXdC6dWsAQM+ePdG8eXOMHDkSixcvRmZmJubMmYOJEyfCwMDgWQ4RERER1VA1Kqny9/fHtWvXEBwcjMzMTLi5uSE6Olq6KDw9PR06Ov9OrnXo0AFbt27FnDlz8OGHH8LFxQVRUVFo2bKl1GbXrl1SUgYAw4YNAwCEhIRg3rx5UCqV2L9/v5TAOTg4YPDgwZgzZ460ja6uLn788Ue888478Pb2homJCQICArBgwYKqHhIiIiKqJWrUOlXPO65T9QjXqSIiotqk1q1TRURERFSbMakiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikgGTKiIiIiIZMKkiIiIikoFeZTZ+8OABMjMzcffuXdjY2KBu3bpyxUVERERUq1R4puqff/7BunXr0LVrV5ibm8PJyQnNmjWDjY0NHB0dMX78eBw9erQqYiUiIiKqsSqUVC1fvhxOTk7YtGkTfHx8EBUVhaSkJPz555+Ii4tDSEgIHj58iJ49e6JXr144d+5cVcVNREREVKNU6PTf0aNH8euvv6JFixYa6z09PTFmzBiEhYVh06ZN+O233+Di4iJLoEREREQ1mUIIIao7iBdFbm4uLCwskJOTA3Nz8+oOp8p5eGguP3bs2cZBRERUGeX9/uav/4iIiIhkoHVSde/ePdy9e1d6funSJaxYsQL79u2TJTAiIiKi2kTrpGrAgAHYsmULAOD27dvw8vLCsmXLMGDAAKxbt062AImIiIhqA62TqoSEBHTu3BkAsH37dtja2uLSpUvYsmULVq1aJVuARERERLWB1knV3bt3YWZmBgDYt28fBg0aBB0dHbRv3x6XLl2SLUAiIiKi2kDrpKpx48aIiorC5cuXsXfvXvTs2RMAkJ2d/UL8so2IiIjocVonVcHBwZgxYwacnJzg6ekJb29vAI9mrdq0aSNbgERERES1gdb3/vPy8kJ6ejoyMjLg5uYmlffo0QODBg2SIzYiIiKiWkPrmSpHR0e0bNkSM2bMwIwZM/D111/j5MmT0NXVxaeffqp1QGvXroWTkxMMDQ3h5eWFI0eOlNk+MjISTZs2haGhIVq1aoXdu3er1e/YsQM9e/aElZUVFAoFkpKS1Opv3ryJyZMno0mTJjAyMkLDhg0xZcoU5OTkqLVTKBQlHtu2bdP6OImIiOj5ovVMVVpaGhITE5GUlITExER8//33uHr1KgBofU1VREQEgoKCEBYWBi8vL6xYsQK+vr44e/Ys6tWrV6L9oUOHMHz4cISGhqJv377YunUr/Pz8kJCQgJYtWwIA8vLy0KlTJ7z++usYP358iT6uXr2Kq1evYunSpWjevDkuXbqEt99+G1evXsX27dvV2m7atAm9evWSnltaWmp1nERERPT8kfU2NXFxcQgICMCCBQswbNiwCm/v5eWFdu3aYc2aNQCAoqIiODg4YPLkyZg1a1aJ9v7+/sjLy8OPP/4olbVv3x5ubm4ICwtTa3vx4kU4OzsjMTFR7XSlJpGRkXjzzTeRl5cHPb1HeadCocDOnTvh5+dX4eMqxtvUPMLb1BARUW1SLbep8fb2xsqVK7F06dIKb1tQUIDjx4/Dx8fn3+B0dODj44O4uDiN28TFxam1BwBfX99S25dX8aAVJ1TFJk6cCGtra3h6emLjxo14Wj6an5+P3NxctQcRERE9n7Q+/VdQUAClUlmi3MXFBadOnapwf9evX0dhYSFsbW3Vym1tbXHmzBmN22RmZmpsn5mZWeH9Px7HwoULMWHCBLXyBQsW4JVXXoGxsTH27duHd999F3fu3MGUKVNK7Ss0NBTz58/XOhYiIiKqPbROqkxNTdG8eXO0adMGbm5uaNOmDezt7bF69eoSs0e1RW5uLvr06YPmzZtj3rx5anVz586V/r9NmzbIy8vDkiVLykyqZs+ejaCgILX+HRwcZI+biIiIqp/Wp/9+/vlnjB8/Hvr6+vj222/Rq1cvvPzyy1i9ejUKCwsRHByMyMjIUmeZnmRtbQ1dXV1kZWWplWdlZUGlUmncRqVSVah9Wf755x/06tULZmZm2LlzJ/T19cts7+Xlhb///hv5+fmltjEwMIC5ubnag4iIiJ5PWs9UderUCZ06dZKeFxUV4ezZs0hKSkJSUhKOHDmCDRs2IDs7G4WFhU/tT6lUwt3dHbGxsdLF4EVFRYiNjcWkSZM0buPt7Y3Y2FhMmzZNKouJiZEWIi2v3Nxc+Pr6wsDAALt27YKhoeFTt0lKSkKdOnVgYGBQoX0RERHR80nrpOpJOjo6aNasGZo1a4bhw4dL5U/OJJUlKCgIAQEB8PDwgKenJ1asWIG8vDwEBgYCAEaNGoX69esjNDQUADB16lR07doVy5YtQ58+fbBt2zYcO3YM69evl/q8efMm0tPTpeUezp49C+DRLJdKpUJubi569uyJu3fv4ptvvlG7oNzGxga6urr44YcfkJWVhfbt28PQ0BAxMTH49NNPMWPGjMoNGhERET03ZEuqSvPkheRl8ff3x7Vr1xAcHIzMzEy4ubkhOjpa6iM9PR06Ov+esezQoQO2bt2KOXPm4MMPP4SLiwuioqKkNaoAYNeuXVJSBkBa6iEkJATz5s1DQkIC4uPjATy6n+Hj0tLS4OTkBH19faxduxbTp0+HEAKNGzfG8uXLNa57RURERC8mWdeporJxnapHuE4VERHVJtWyThURERHRi4pJFREREZEMKpxUBQcH4/jx41URCxEREVGtVeGk6u+//0bv3r3RoEEDvPPOO9izZw8KCgqqIjYiIiKiWqPCSdXGjRuRmZmJ7777DmZmZpg2bRqsra0xePBgbNmyBTdv3qyKOImIiIhqNK2uqdLR0UHnzp2xePFinD17FvHx8fDy8sKXX34Je3t7dOnSBUuXLsWVK1fkjpeIiIioRpLlQvVmzZph5syZ+OOPP3D58mUEBATgt99+w3fffSdH90REREQ1Htepeoa4TtUjXKeKiIhqE65TRURERPQMMakiIiIikgGTKiIiIiIZMKkiIiIikoHsSdXly5cxZswYubslIiIiqtFkT6pu3ryJzZs3y90tERERUY2mV9ENdu3aVWb9X3/9pXUwRERERLVVhZMqPz8/KBQKlLW8lUKhqFRQRERERLVNhU//2dnZYceOHSgqKtL4SEhIqIo4iYiIiGq0CidV7u7uOH78eKn1T5vFIiIiInoeVfj03/vvv4+8vLxS6xs3bowDBw5UKigiIiKi2qbCSVXnzp3LrDcxMUHXrl21DoiIiIioNqrQ6b/09PQKdX7lypUKtSciIiKqrSqUVLVr1w5vvfUWjh49WmqbnJwcbNiwAS1btsR///vfSgdIREREVBtU6PTf6dOn8cknn+DVV1+FoaEh3N3dYW9vD0NDQ9y6dQunT5/GqVOn0LZtWyxevBivvfZaVcVNREREVKMohBY/1bt37x5++ukn/P7777h06RLu3bsHa2trtGnTBr6+vmjZsmVVxFrr5ebmwsLCAjk5OTA3N6/ucKqch4fm8mPHnm0cRERElVHe7+8KX6gOAEZGRhgyZAiGDBmidYBEREREzxPZ7/1HRERE9CJiUkVEREQkAyZVRERERDJgUkVEREQkAyZVRERERDLQOqm6d+8e7t69Kz2/dOkSVqxYgX379skSGBEREVFtonVSNWDAAGzZsgUAcPv2bXh5eWHZsmUYMGAA1q1bJ1uARERERLWB1klVQkKCdHPl7du3w9bWFpcuXcKWLVuwatUq2QIkIiIiqg20Tqru3r0LMzMzAMC+ffswaNAg6OjooH379rh06ZJsARIRERHVBlonVY0bN0ZUVBQuX76MvXv3omfPngCA7OzsF+IWLERERESP0zqpCg4OxowZM+Dk5ARPT094e3sDeDRr1aZNG9kCJCIiIqoNtLr3HwAMGTIEnTp1QkZGBtzc3KTyHj16YODAgXLERkRERFRrVGqdqnPnzmHZsmXo2LEjrly5AgA4e/Ysrl+/LktwRERERLWF1knVf//7X/j6+sLIyAgJCQnIz88HAOTk5ODTTz+VLUAiIiKi2kDrpOrjjz9GWFgYNmzYAH19fam8Y8eOSEhI0DqgtWvXwsnJCYaGhvDy8sKRI0fKbB8ZGYmmTZvC0NAQrVq1wu7du9Xqd+zYgZ49e8LKygoKhQJJSUkl+rh//z4mTpwIKysrmJqaYvDgwcjKylJrk56ejj59+sDY2Bj16tXD+++/j4cPH2p9nERERPR80TqpOnv2LLp06VKi3MLCArdv39aqz4iICAQFBSEkJAQJCQlwdXWFr68vsrOzNbY/dOgQhg8fjrFjxyIxMRF+fn7w8/NDSkqK1CYvLw+dOnXCZ599Vup+p0+fjh9++AGRkZE4ePAgrl69ikGDBkn1hYWF6NOnDwoKCnDo0CFs3rwZ4eHhCA4O1uo4iYiI6DkktOTs7CxiYmKEEEKYmpqKCxcuCCGE2Lx5s2jWrJlWfXp6eoqJEydKzwsLC4W9vb0IDQ3V2P71118Xffr0USvz8vISb731Vom2aWlpAoBITExUK799+7bQ19cXkZGRUllqaqoAIOLi4oQQQuzevVvo6OiIzMxMqc26deuEubm5yM/PL/fx5eTkCAAiJyen3NvUZu7umh9ERES1SXm/v7WeqRo/fjymTp2K+Ph4KBQKXL16Fd9++y1mzJiBd955p8L9FRQU4Pjx4/Dx8ZHKdHR04OPjg7i4OI3bxMXFqbUHAF9f31Lba3L8+HE8ePBArZ+mTZuiYcOGUj9xcXFo1aoVbG1t1faTm5uLU6dOldp3fn4+cnNz1R5ERET0fNJ6SYVZs2ahqKgIPXr0wN27d9GlSxcYGBhgxowZmDx5coX7u379OgoLC9USFwCwtbXFmTNnNG6TmZmpsX1mZma595uZmQmlUglLS8tS+yltP8V1pQkNDcX8+fPLHQsRERHVXlrPVCkUCnz00Ue4efMmUlJScPjwYVy7dg0LFy6UM75abfbs2cjJyZEely9fru6QiIiIqIpoPVOVnp4OBwcHKJVKNG/evERdw4YNK9SftbU1dHV1S/zqLisrCyqVSuM2KpWqQu1L66OgoAC3b99Wm616vB+VSlXiV4jF+y1rXwYGBjAwMCh3LERERFR7aT1T5ezsjGvXrpUov3HjBpydnSvcn1KphLu7O2JjY6WyoqIixMbGSrfAeZK3t7daewCIiYkptb0m7u7u0NfXV+vn7NmzSE9Pl/rx9vbGyZMn1X6FGBMTA3Nz8xIJJREREb2YtJ6pEkJAoVCUKL9z5w4MDQ216jMoKAgBAQHw8PCAp6cnVqxYgby8PAQGBgIARo0ahfr16yM0NBQAMHXqVHTt2hXLli1Dnz59sG3bNhw7dgzr16+X+rx58ybS09Nx9epVAI8SJuDRDJNKpYKFhQXGjh2LoKAg1K1bF+bm5pg8eTK8vb3Rvn17AEDPnj3RvHlzjBw5EosXL0ZmZibmzJmDiRMnciaKiIiIAGiRVAUFBQF4dE3V3LlzYWxsLNUVFhYiPj5e7V6AFeHv749r164hODgYmZmZcHNzQ3R0tHRReHp6OnR0/p1c69ChA7Zu3Yo5c+bgww8/hIuLC6KiotCyZUupza5du6SkDACGDRsGAAgJCcG8efMAAJ9//jl0dHQwePBg5Ofnw9fXF1988YW0ja6uLn788Ue888478Pb2homJCQICArBgwQKtjpOIiIiePwohhKjIBt27dwcAHDx4EN7e3lAqlVKdUqmEk5MTZsyYARcXF3kjfQ7k5ubCwsICOTk5MDc3r+5wqpyHh+byY8eebRxERESVUd7v7wrPVB04cAAAEBgYiJUrV74QyQERERHR02h9TdWmTZsAAKdPn0Z6ejoKCgrU6vv371+5yIiIiIhqEa2TqrS0NPj5+eHkyZNQKBQoPotYfPF6YWGhPBESERER1QJaL6kwZcoUODs7Izs7G8bGxjh16hR+/fVXeHh44JdffpExRCIiIqKaT+uZqri4OPz888+wtraGjo4OdHR00KlTJ4SGhmLKlClITEyUM04iIiKiGk3rmarCwkKYmZkBeLQaevE6UI6OjtJaUEREREQvCq1nqlq2bInk5GQ4OzvDy8sLixcvhlKpxPr16/HSSy/JGSMRERFRjad1UjVnzhzk5eUBABYsWIC+ffuic+fOsLKyQkREhGwBEhEREdUGWidVvr6+0v83btwYZ86cwc2bN1GnTh2Nt68hIiIiep5pfU2VJnXr1oVCocC9e/fk7JaIiIioxpM1qcrPz8eyZcvg7OwsZ7dERERENV6Fk6r8/HzMnj0bHh4e6NChA6KiogA8WmHd2dkZK1aswPTp0+WOk4iIiKhGq/A1VcHBwfjyyy/h4+ODQ4cOYejQoQgMDMThw4exfPlyDB06FLq6ulURKxEREVGNVeGkKjIyElu2bEH//v2RkpKC1q1b4+HDh0hOTuYF6kRERPTCqvDpv7///hvu7u4AHq1VZWBggOnTpzOhIiIiohdahZOqwsJCKJVK6bmenh5MTU1lDYqIiIiotqnw6T8hBEaPHg0DAwMAwP379/H222/DxMRErd2OHTvkiZCIiIioFqhwUhUQEKD2/M0335QtGCIiIqLaqsJJ1aZNm6oiDiIiIqJaTdbFP4mIiIheVFrf+4+IiAB4eJQsO3bs2cdBRNWOM1VEREREMmBSRURERCQDnv4jIiqLptN7AE/xEVEJlUqqYmNjERsbi+zsbBQVFanVbdy4sVKBEREREdUmWidV8+fPx4IFC+Dh4QE7OzvepoaIiIheaFonVWFhYQgPD8fIkSPljIeIiIioVtL6QvWCggJ06NBBzliIiIiIai2tk6px48Zh69atcsZCREREVGtpffrv/v37WL9+Pfbv34/WrVtDX19frX758uWVDo6IiIiottA6qTpx4gTc3NwAACkpKWp1vGidiIiIXjRaJ1UHDhyQMw4iIiKiWo0rqhMRERHJoFKLf96+fRtfffUVUlNTAQDNmzfH2LFjYWFhIUtwRERERLWF1jNVx44dQ6NGjfD555/j5s2buHnzJj7//HM0atQICQkJcsZIREREVONpPVM1ffp09O/fHxs2bICe3qNuHj58iHHjxmHatGn49ddfZQuSiIiIqKbTOqk6duyYWkIFAHp6epg5cyY8SrsBKREREdFzSuvTf+bm5khPTy9RfvnyZZiZmVUqKCIiIqLaRuukyt/fH2PHjkVERAQuX76My5cvY9u2bRg3bhyGDx8uZ4xERERENZ7WSdXSpUsxaNAgjBo1Ck5OTnB0dMTo0aMxZMgQfPbZZ5UKau3atXBycoKhoSG8vLxw5MiRMttHRkaiadOmMDQ0RKtWrbB79261eiEEgoODYWdnByMjI/j4+ODcuXNS/S+//AKFQqHxcfToUQDAxYsXNdYfPny4UsdKREREzwetkyqlUomVK1fi1q1bSEpKQnJysvQLQAMDA60DioiIQFBQEEJCQpCQkABXV1f4+voiOztbY/tDhw5h+PDhGDt2LBITE+Hn5wc/Pz+1Vd4XL16MVatWISwsDPHx8TAxMYGvry/u378PAOjQoQMyMjLUHuPGjYOzs3OJ68P279+v1s7d3V3rYyUiIqLnh0IIIcrbOCgoCAsXLoSJiQmCgoLKbKvtvf+8vLzQrl07rFmzBgBQVFQEBwcHTJ48GbNmzSrR3t/fH3l5efjxxx+lsvbt28PNzQ1hYWEQQsDe3h7vvfceZsyYAQDIycmBra0twsPDMWzYsBJ9PnjwAPXr18fkyZMxd+5cAI9mqpydnZGYmCjdnqeicnNzYWFhgZycHJibm2vVR21S2u8Vjh17tnEQVcrT3sia6vkmJ3qulPf7u0K//ktMTMSDBw+k/y+Ntvf+KygowPHjxzF79mypTEdHBz4+PoiLi9O4TVxcXIkEz9fXF1FRUQCAtLQ0ZGZmwsfHR6q3sLCAl5cX4uLiNCZVu3btwo0bNxAYGFiirn///rh//z5efvllzJw5E/379y/1ePLz85Gfny89z83NLbUtERER1W4VSqoev99fVdz77/r16ygsLIStra1aua2tLc6cOaNxm8zMTI3tMzMzpfristLaPOmrr76Cr68vGjRoIJWZmppi2bJl6NixI3R0dPDf//4Xfn5+iIqKKjWxCg0Nxfz588s4YiIiInpeaL1OVXp6OhwcHDTOSqWnp6Nhw4aVCqy6/P3339i7dy++//57tXJra2u1GbF27drh6tWrWLJkSalJ1ezZs9W2yc3NhYODQ9UETkRERNVK6wvVnZ2dce3atRLlN27cgLOzs1Z9WltbQ1dXF1lZWWrlWVlZUKlUGrdRqVRlti/+b3n73LRpE6ysrMo8rVfMy8sL58+fL7XewMAA5ubmag8iIiJ6PmmdVAkhNM5S3blzB4aGhlr1qVQq4e7ujtjYWKmsqKgIsbGx8Pb21riNt7e3WnsAiImJkdo7OztDpVKptcnNzUV8fHyJPoUQ2LRpE0aNGgV9ff2nxpuUlAQ7O7tyHx+p8/DQ/CAiIqqNKnz6r/h0lkKhwNy5c2FsbCzVFRYWIj4+XutfxxX3HxAQAA8PD3h6emLFihXIy8uTLhofNWoU6tevj9DQUADA1KlT0bVrVyxbtgx9+vTBtm3bcOzYMaxfv16Kc9q0afj444/h4uICZ2dnzJ07F/b29vDz81Pb988//4y0tDSMGzeuRFybN2+GUqlEmzZtAAA7duzAxo0b8Z///EfrYyUiIqLnR4WTquJf/QkhcPLkSSiVSqlOqVTC1dVVWrpAG/7+/rh27RqCg4ORmZkJNzc3REdHSxeap6enQ0fn3wm2Dh06YOvWrZgzZw4+/PBDuLi4ICoqCi1btpTazJw5E3l5eZgwYQJu376NTp06ITo6usSM2ldffYUOHTqgadOmGmNbuHAhLl26BD09PTRt2hQREREYMmSI1sdKREREz48KrVP1uMDAQKxcuZLXCVUA16l6pKzlfR6vJ6oRuE4V0QuvStapetymTZu03ZSIiIjouVOhpOpZrKhOREREVBvVqBXViYiIiGqrGrWiOhEREVFtpfU6Vffu3cPdu3el55cuXcKKFSuwb98+WQIjIiIiqk20TqoGDBiALVu2AABu374NT09PLFu2DAMGDMC6detkC5CIiIioNtA6qUpISEDnzp0BANu3b4dKpcKlS5ewZcsWrFq1SrYAiYiIiGoDrZOqu3fvwszMDACwb98+DBo0CDo6Omjfvj0uXbokW4BEREREtYHWSVXjxo0RFRWFy5cvY+/evejZsycAIDs7+4VY2JKIiIjocVonVcHBwZgxYwacnJzg5eUl3Zx437590v3xiIiIiF4UWq+oPmTIEHTq1AkZGRlwdXWVynv06IGBAwfKEhwRERFRbaF1UgUAKpUKKpVKrczT07NSARERERHVRpVKqmJjYxEbG4vs7GwUFRWp1W3cuLFSgRERERHVJlonVfPnz8eCBQvg4eEBOzs73pqGiIiIXmhaJ1VhYWEIDw/HyJEj5YyHiIiIqFbS+td/BQUF6NChg5yxEBEREdVaWidV48aNw9atW+WMhYiIiKjW0vr03/3797F+/Xrs378frVu3hr6+vlr98uXLKx0cERERUW2hdVJ14sQJuLm5AQBSUlLU6njROhEREb1otE6qDhw4IGccRERERLVapdapAoDTp08jPT0dBQUFUplCoUC/fv0q2zURERFRraF1UvXXX39h4MCBOHnyJBQKBYQQAP499VdYWChPhERERES1gNa//ps6dSqcnZ2RnZ0NY2NjnDp1Cr/++is8PDzwyy+/yBgiERERUc2n9UxVXFwcfv75Z1hbW0NHRwc6Ojro1KkTQkNDMWXKFCQmJsoZJxEREVGNpvVMVWFhIczMzAAA1tbWuHr1KgDA0dERZ8+elSc6IiIiolpC65mqli1bIjk5Gc7OzvDy8sLixYuhVCqxfv16vPTSS3LGSERERFTjaZ1UzZkzB3l5eQCABQsWoG/fvujcuTOsrKwQEREhW4BEREREtYHWSZWvr6/0/40bN8aZM2dw8+ZN1KlTh4t/EhER0QtHq2uqHjx4gB49euDcuXNq5XXr1mVCRURERC8krZIqfX19nDhxQu5YiIiIiGotrX/99+abb+Krr76SMxYiIiKiWkvra6oePnyIjRs3Yv/+/XB3d4eJiYla/fLlyysdHBEREVFtoXVSlZKSgrZt2wIA/vzzT7U6XldFREREL5oKJ1ULFizAjBkzcODAgaqIh4iIiKhWqvA1VfPnz8edO3eqIhYiIiKiWqvCSZUQoiriICIiIqrVtPr1H6+ZIiIiIlKn1YXqL7/88lMTq5s3b2oVEBEREVFtpFVSNX/+fFhYWMgdCxEREVGtpdXpv2HDhiEgIKDMR2WsXbsWTk5OMDQ0hJeXF44cOVJm+8jISDRt2hSGhoZo1aoVdu/erVYvhEBwcDDs7OxgZGQEHx+fErfYcXJygkKhUHssWrRIrc2JEyfQuXNnGBoawsHBAYsXL67UcRIREdHzo8JJVVVfTxUREYGgoCCEhIQgISEBrq6u8PX1RXZ2tsb2hw4dwvDhwzF27FgkJibCz88Pfn5+SElJkdosXrwYq1atQlhYGOLj42FiYgJfX1/cv39fra8FCxYgIyNDekyePFmqy83NRc+ePeHo6Ijjx49jyZIlmDdvHtavX181A0FERES1So379d/y5csxfvx4BAYGonnz5ggLC4OxsTE2btyosf3KlSvRq1cvvP/++2jWrBkWLlyItm3bYs2aNVK8K1aswJw5czBgwAC0bt0aW7ZswdWrVxEVFaXWl5mZGVQqlfR4fJX4b7/9FgUFBdi4cSNatGiBYcOGYcqUKVw5noiIiABokVQVFRWhXr16VRELCgoKcPz4cfj4+EhlOjo68PHxQVxcnMZt4uLi1NoDgK+vr9Q+LS0NmZmZam0sLCzg5eVVos9FixbBysoKbdq0wZIlS/Dw4UO1/XTp0gVKpVJtP2fPnsWtW7c0xpafn4/c3Fy1BxERET2ftL5NTVW4fv06CgsLYWtrq1Zua2uLM2fOaNwmMzNTY/vMzEypvristDYAMGXKFLRt2xZ169bFoUOHMHv2bGRkZEgzUZmZmXB2di7RR3FdnTp1SsQWGhqK+fPnP/W4iYiIqParUUlVdQoKCpL+v3Xr1lAqlXjrrbcQGhoKAwMDrfqcPXu2Wr+5ublwcHCodKxERERU82j167+qYm1tDV1dXWRlZamVZ2VlQaVSadxGpVKV2b74vxXpEwC8vLzw8OFDXLx4scz9PL6PJxkYGMDc3FztQURERM+nGpVUKZVKuLu7IzY2ViorKipCbGwsvL29NW7j7e2t1h4AYmJipPbOzs5QqVRqbXJzcxEfH19qnwCQlJQEHR0d6foxb29v/Prrr3jw4IHafpo0aaLx1B8RERG9WGpUUgU8Og23YcMGbN68GampqXjnnXeQl5eHwMBAAMCoUaMwe/Zsqf3UqVMRHR2NZcuW4cyZM5g3bx6OHTuGSZMmAXi0BMS0adPw8ccfY9euXTh58iRGjRoFe3t7+Pn5AXh0EfqKFSuQnJyMv/76C99++y2mT5+ON998U0qY3njjDSiVSowdOxanTp1CREQEVq5cqXZ6j4iIiF5cNe6aKn9/f1y7dg3BwcHIzMyEm5sboqOjpYvC09PToaPzby7YoUMHbN26FXPmzMGHH34IFxcXREVFoWXLllKbmTNnIi8vDxMmTMDt27fRqVMnREdHw9DQEMCj03Tbtm3DvHnzkJ+fD2dnZ0yfPl0tYbKwsMC+ffswceJEuLu7w9raGsHBwZgwYcIzGhkiIiKqyRSiqheeIklubi4sLCyQk5PzQlxf5eGhufzYsfLVE9UI2ryR+SYmeq6U9/u7xp3+IyIiIqqNmFQRERERyYBJFREREZEMmFQRERERyYBJFREREZEMmFQRERERyYBJFREREZEMmFQRERERyYBJFREREZEMmFQRERERyYBJFREREZEMmFQRERERyYBJFREREZEMmFQRERERyUCvugMg0sTDQ3P5sWPPNg4iIqLy4kwVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJoEYmVWvXroWTkxMMDQ3h5eWFI0eOlNk+MjISTZs2haGhIVq1aoXdu3er1QshEBwcDDs7OxgZGcHHxwfnzp2T6i9evIixY8fC2dkZRkZGaNSoEUJCQlBQUKDWRqFQlHgcPnxY3oMnIiKiWqnGJVUREREICgpCSEgIEhIS4OrqCl9fX2RnZ2tsf+jQIQwfPhxjx45FYmIi/Pz84Ofnh5SUFKnN4sWLsWrVKoSFhSE+Ph4mJibw9fXF/fv3AQBnzpxBUVERvvzyS5w6dQqff/45wsLC8OGHH5bY3/79+5GRkSE93N3dq2YgiIiIqFZRCCFEdQfxOC8vL7Rr1w5r1qwBABQVFcHBwQGTJ0/GrFmzSrT39/dHXl4efvzxR6msffv2cHNzQ1hYGIQQsLe3x3vvvYcZM2YAAHJycmBra4vw8HAMGzZMYxxLlizBunXr8NdffwF4NFPl7OyMxMREuLm5aXVsubm5sLCwQE5ODszNzbXqozbx8NBcfuzY0+ufti3RM6PNG5lvVKLnSnm/v2vUTFVBQQGOHz8OHx8fqUxHRwc+Pj6Ii4vTuE1cXJxaewDw9fWV2qelpSEzM1OtjYWFBby8vErtE3iUeNWtW7dEef/+/VGvXj106tQJu3btKvN48vPzkZubq/YgIiKi51ONSqquX7+OwsJC2NraqpXb2toiMzNT4zaZmZllti/+b0X6PH/+PFavXo233npLKjM1NcWyZcsQGRmJn376CZ06dYKfn1+ZiVVoaCgsLCykh4ODQ6ltiYiIqHbTq+4AaporV66gV69eGDp0KMaPHy+VW1tbIygoSHrerl07XL16FUuWLEH//v019jV79my1bXJzc5lYERERPadqVFJlbW0NXV1dZGVlqZVnZWVBpVJp3EalUpXZvvi/WVlZsLOzU2vz5LVRV69eRffu3dGhQwesX7/+qfF6eXkhJiam1HoDAwMYGBg8tR+qGF5vRURENVGNOv2nVCrh7u6O2NhYqayoqAixsbHw9vbWuI23t7daewCIiYmR2js7O0OlUqm1yc3NRXx8vFqfV65cQbdu3eDu7o5NmzZBR+fpQ5OUlKSWqBEREdGLq0bNVAFAUFAQAgIC4OHhAU9PT6xYsQJ5eXkIDAwEAIwaNQr169dHaGgoAGDq1Kno2rUrli1bhj59+mDbtm04duyYNNOkUCgwbdo0fPzxx3BxcYGzszPmzp0Le3t7+Pn5Afg3oXJ0dMTSpUtx7do1KZ7ima7NmzdDqVSiTZs2AIAdO3Zg48aN+M9//vOshoaIiIhqsBqXVPn7++PatWsIDg5GZmYm3NzcEB0dLV1onp6erjaL1KFDB2zduhVz5szBhx9+CBcXF0RFRaFly5ZSm5kzZyIvLw8TJkzA7du30alTJ0RHR8PQ0BDAo5mt8+fP4/z582jQoIFaPI+vOLFw4UJcunQJenp6aNq0KSIiIjBkyJCqHA4iIiKqJWrcOlXPM65T9Uhl16niNVX0THGdKqIXXq1cp4qIiIiotmJSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMmBSRURERCQDJlVEREREMqhxN1QmqizeG5CIiKoDZ6qIiIiIZMCZKiIiTm8SkQw4U0VEREQkAyZVRERERDJgUkVEREQkAyZVRERERDJgUkVEREQkAyZVRERERDLgkgr0QuEv5+mZ4huO6IXCmSoiIiIiGTCpIiIiIpIBkyoiIiIiGTCpIiIiIpIBkyoiIiIiGfDXf0SP4Y+1iIhIW5ypIiIiIpIBkyoiIiIiGfD0HxE9/3hel4ieASZVROXE72UiIioLT/8RERERyYAzVUQy4CwWERExqSIiqg7MxImeO0yqiJ4Bfn8+AxxkIqpmvKaKiIiISAacqSKqZpxgISJ6PjCpIqrhmHT9vxdtIF604yV6DtTI039r166Fk5MTDA0N4eXlhSNHjpTZPjIyEk2bNoWhoSFatWqF3bt3q9ULIRAcHAw7OzsYGRnBx8cH586dU2tz8+ZNjBgxAubm5rC0tMTYsWNx584dtTYnTpxA586dYWhoCAcHByxevFieAybSkoeH5kdZ9TXa0w6oGpxO1fwgInpSjZupioiIQFBQEMLCwuDl5YUVK1bA19cXZ8+eRb169Uq0P3ToEIYPH47Q0FD07dsXW7duhZ+fHxISEtCyZUsAwOLFi7Fq1Sps3rwZzs7OmDt3Lnx9fXH69GkYGhoCAEaMGIGMjAzExMTgwYMHCAwMxIQJE7B161YAQG5uLnr27AkfHx+EhYXh5MmTGDNmDCwtLTFhwoRnN0BEMilrIqRKJ0k4A1N5HEOiGqnGJVXLly/H+PHjERgYCAAICwvDTz/9hI0bN2LWrFkl2q9cuRK9evXC+++/DwBYuHAhYmJisGbNGoSFhUEIgRUrVmDOnDkYMGAAAGDLli2wtbVFVFQUhg0bhtTUVERHR+Po0aPw+P8/VqtXr8Zrr72GpUuXwt7eHt9++y0KCgqwceNGKJVKtGjRAklJSVi+fDmTKnrhnDbR/KXePO9YmXVA6bM8zStRV5l+n7ZtrcSki6ha1KikqqCgAMePH8fs2bOlMh0dHfj4+CAuLk7jNnFxcQgKClIr8/X1RVRUFAAgLS0NmZmZ8PHxkeotLCzg5eWFuLg4DBs2DHFxcbC0tJQSKgDw8fGBjo4O4uPjMXDgQMTFxaFLly5QKpVq+/nss89w69Yt1KlTR44hIJLVllRNX67Hyqh7VF9WHdVylZmiZLJGVKYalVRdv34dhYWFsLW1VSu3tbXFmTNnNG6TmZmpsX1mZqZUX1xWVpsnTy3q6emhbt26am2cnZ1L9FFcpympys/PR35+vvQ8JycHwKNTiS+CwkLN5cWHX1a9tnWV6bcqY/ridNdS6g+WWVfZbe+IkoEVv/801RXXl1VXXdvWxJgAaH7xq/PNWJltK9NvV83vRRx89F4ss16bbauq36retibGRE9V/HkXQpTZrkYlVc+b0NBQzJ8/v0S5g4NDNURTc1hYaF+vbV11bfu0fmtkULVtW8ZU9dsypqrftibGRCX8888/sChjzGpUUmVtbQ1dXV1kZWWplWdlZUGlUmncRqVSldm++L9ZWVmws7NTa+Pm5ia1yc7OVuvj4cOHuHnzplo/mvbz+D6eNHv2bLVTk0VFRbh58yasrKygUCg0blNZubm5cHBwwOXLl2Fubl4l+3gecJzKh+NUfhyr8uE4lQ/HqXye1TgJIfDPP//A3t6+zHY1KqlSKpVwd3dHbGws/Pz8ADxKRGJjYzFp0iSN23h7eyM2NhbTpk2TymJiYuDt7Q0AcHZ2hkqlQmxsrJRE5ebmIj4+Hu+8847Ux+3bt3H8+HG4u7sDAH7++WcUFRXBy8tLavPRRx/hwYMH0NfXl/bTpEmTUq+nMjAwgIGBgVqZpaVlhcdFG+bm5vwglgPHqXw4TuXHsSofjlP5cJzK51mMU1kzVBJRw2zbtk0YGBiI8PBwcfr0aTFhwgRhaWkpMjMzhRBCjBw5UsyaNUtq/8cffwg9PT2xdOlSkZqaKkJCQoS+vr44efKk1GbRokXC0tJS/O9//xMnTpwQAwYMEM7OzuLevXtSm169eok2bdqI+Ph48fvvvwsXFxcxfPhwqf727dvC1tZWjBw5UqSkpIht27YJY2Nj8eWXXz6DUSm/nJwcAUDk5ORUdyg1GsepfDhO5cexKh+OU/lwnMqnpo1TjZqpAgB/f39cu3YNwcHByMzMhJubG6Kjo6WLwtPT06Gj8++apR06dMDWrVsxZ84cfPjhh3BxcUFUVJS0RhUAzJw5E3l5eZgwYQJu376NTp06ITo6WlqjCgC+/fZbTJo0CT169ICOjg4GDx6MVatWSfUWFhbYt28fJk6cCHd3d1hbWyM4OJjLKRAREdEj1Z3Vkbzu378vQkJCxP3796s7lBqN41Q+HKfy41iVD8epfDhO5VPTxkkhxFN+H0hERERET1Uj7/1HREREVNswqSIiIiKSAZMqIiIiIhkwqSIiIiKSAZOq58zatWvh5OQEQ0NDeHl54ciRI9UdUrX69ddf0a9fP9jb20OhUEg32i4mhEBwcDDs7OxgZGQEHx8fnDt3rnqCrUahoaFo164dzMzMUK9ePfj5+eHs2bNqbe7fv4+JEyfCysoKpqamGDx4cIm7DDzv1q1bh9atW0sLDXp7e2PPnj1SPcdIs0WLFkGhUKgt0syxAubNmweFQqH2aNq0qVTPMfrXlStX8Oabb8LKygpGRkZo1aoVjj12I++a8recSdVzJCIiAkFBQQgJCUFCQgJcXV3h6+tb4hY8L5K8vDy4urpi7dq1GusXL16MVatWISwsDPHx8TAxMYGvry/u37//jCOtXgcPHsTEiRNx+PBhxMTE4MGDB+jZsyfy8vKkNtOnT8cPP/yAyMhIHDx4EFevXsWgQYOqMepnr0GDBli0aBGOHz+OY8eO4ZVXXsGAAQNw6tQpABwjTY4ePYovv/wSrVu3VivnWD3SokULZGRkSI/ff/9dquMYPXLr1i107NgR+vr62LNnD06fPo1ly5ap3c2kxvwtr94VHUhOnp6eYuLEidLzwsJCYW9vL0JDQ6sxqpoDgNi5c6f0vKioSKhUKrFkyRKp7Pbt28LAwEB899131RBhzZGdnS0AiIMHDwohHo2Lvr6+iIyMlNqkpqYKACIuLq66wqwR6tSpI/7zn/9wjDT4559/hIuLi4iJiRFdu3YVU6dOFULw/VQsJCREuLq6aqzjGP3rgw8+EJ06dSq1vib9LedM1XOioKAAx48fh4+Pj1Smo6MDHx8fxMXFVWNkNVdaWhoyMzPVxszCwgJeXl4v/Jjl5OQAAOrWrQsAOH78OB48eKA2Vk2bNkXDhg1f2LEqLCzEtm3bkJeXB29vb46RBhMnTkSfPn3UxgTg++lx586dg729PV566SWMGDEC6enpADhGj9u1axc8PDwwdOhQ1KtXD23atMGGDRuk+pr0t5xJ1XPi+vXrKCwslG7nU8zW1haZmZnVFFXNVjwuHDN1RUVFmDZtGjp27Cjd7ikzMxNKpbLEDcFfxLE6efIkTE1NYWBggLfffhs7d+5E8+bNOUZP2LZtGxISEhAaGlqijmP1iJeXF8LDwxEdHY1169YhLS0NnTt3xj///MMxesxff/2FdevWwcXFBXv37sU777yDKVOmYPPmzQBq1t/yGnfvPyKqXhMnTkRKSoratR30ryZNmiApKQk5OTnYvn07AgICcPDgweoOq0a5fPkypk6dipiYGLV7rJK63r17S//funVreHl5wdHREd9//z2MjIyqMbKapaioCB4eHvj0008BAG3atEFKSgrCwsIQEBBQzdGp40zVc8La2hq6urolfhmSlZUFlUpVTVHVbMXjwjH716RJk/Djjz/iwIEDaNCggVSuUqlQUFCA27dvq7V/EcdKqVSicePGcHd3R2hoKFxdXbFy5UqO0WOOHz+O7OxstG3bFnp6etDT08PBgwexatUq6OnpwdbWlmOlgaWlJV5++WWcP3+e76fH2NnZoXnz5mplzZo1k06V1qS/5UyqnhNKpRLu7u6IjY2VyoqKihAbGwtvb+9qjKzmcnZ2hkqlUhuz3NxcxMfHv3BjJoTApEmTsHPnTvz8889wdnZWq3d3d4e+vr7aWJ09exbp6ekv3Fg9qaioCPn5+Ryjx/To0QMnT55EUlKS9PDw8MCIESOk/+dYlXTnzh1cuHABdnZ2fD89pmPHjiWWePnzzz/h6OgIoIb9LX+ml8VTldq2bZswMDAQ4eHh4vTp02LChAnC0tJSZGZmVndo1eaff/4RiYmJIjExUQAQy5cvF4mJieLSpUtCCCEWLVokLC0txf/+9z9x4sQJMWDAAOHs7Czu3btXzZE/W++8846wsLAQv/zyi8jIyJAed+/eldq8/fbbomHDhuLnn38Wx44dE97e3sLb27sao372Zs2aJQ4ePCjS0tLEiRMnxKxZs4RCoRD79u0TQnCMyvL4r/+E4FgJIcR7770nfvnlF5GWlib++OMP4ePjI6ytrUV2drYQgmNU7MiRI0JPT0988skn4ty5c+Lbb78VxsbG4ptvvpHa1JS/5UyqnjOrV68WDRs2FEqlUnh6eorDhw9Xd0jV6sCBAwJAiUdAQIAQ4tFPcefOnStsbW2FgYGB6NGjhzh79mz1Bl0NNI0RALFp0yapzb1798S7774r6tSpI4yNjcXAgQNFRkZG9QVdDcaMGSMcHR2FUqkUNjY2okePHlJCJQTHqCxPJlUcKyH8/f2FnZ2dUCqVon79+sLf31+cP39equcY/euHH34QLVu2FAYGBqJp06Zi/fr1avU15W+5Qgghnu3cGBEREdHzh9dUEREREcmASRURERGRDJhUEREREcmASRURERGRDJhUEREREcmASRURERGRDJhUEREREcmASRURERGRDJhUEREREcmASRXRc65bt26YNm1aje2vqvqkqlcdr9uNGzdQr149XLx4sdzbDBs2DMuWLau6oIj+H5Mqolpi9OjRUCgUUCgUUCqVaNy4MRYsWICHDx+Wud2OHTuwcOFC2eKQu7/yunz5MsaMGQN7e3solUo4Ojpi6tSpuHHjxjOPBaj+RLD4/bBo0SK18qioKCgUimqKqup98sknGDBgAJycnMq9zZw5c/DJJ58gJyen6gIjApMqolqlV69eyMjIwLlz5/Dee+9h3rx5WLJkica2BQUFAIC6devCzMxMthjk7q88/vrrL3h4eODcuXP47rvvcP78eYSFhSE2Nhbe3t64efPmM42npjA0NMRnn32GW7duVXcosil+32py9+5dfPXVVxg7dmyF+mzZsiUaNWqEb775prLhEZWJSRVRLWJgYACVSgVHR0e888478PHxwa5duwA8mjmZNGkSpk2bBmtra/j6+krlxTMq3bp1w5QpUzBz5kzUrVsXKpUK8+bNU9tHUVERFi9ejMaNG8PAwAANGzbEJ598ItU/OUNTvN9JkybBwsIC1tbWmDt3Lh6/V3t0dDQ6deoES0tLWFlZoW/fvrhw4UK5j3vixIlQKpXYt28funbtioYNG6J3797Yv38/rly5go8++khq6+TkhBUrVqht7+bmpnacT4vnaeM0evRoHDx4ECtXrpRmD4tPR5Vn/926dcPkyZMxbdo01KlTB7a2ttiwYQPy8vIQGBgIMzMzNG7cGHv27ClzXHx8fKBSqRAaGlpmu6qM6eHDh2W+9kVFRQgNDYWzszOMjIzg6uqK7du3q+1X0/tWk927d8PAwADt27eXyi5cuACFQoEff/wRPXr0gLGxMZo0aYL4+Hi1bfv164dt27aVOU5ElcWkiqgWMzIyUvuX/ebNm6FUKvHHH38gLCxM4zabN2+GiYkJ4uPjsXjxYixYsAAxMTFS/ezZs7Fo0SLMnTsXp0+fxtatW2Fra1tmHJs3b4aenh6OHDmClStXYvny5fjPf/4j1efl5SEoKAjHjh1DbGwsdHR0MHDgQBQVFT31GG/evIm9e/fi3XffhZGRkVqdSqXCiBEjEBERofZF/jTliaescVq5ciW8vb0xfvx4ZGRkICMjAw4ODuXef3H/1tbWOHLkCCZPnox33nkHQ4cORYcOHZCQkICePXti5MiRuHv3bql96Orq4tNPP8Xq1avx999/V2j/csX0tNc+NDQUW7ZsQVhYGE6dOoXp06fjzTffxMGDB9X6eNr7FgB+++03uLu7q5UlJydDoVBg+fLlmDt3LpKTk9GwYUPMmjVLrZ2npyeOHDmC/Pz8yg4TUekEEdUKAQEBYsCAAUIIIYqKikRMTIwwMDAQM2bMEEII0bVrV9GmTZsS23Xt2lVMnTpV+v9OnTqp1bdr10588MEHQgghcnNzhYGBgdiwYUOpcTzeX/HzZs2aiaKiIqnsgw8+EM2aNSu1j2vXrgkA4uTJkxr7fNzhw4cFALFz506N9cuXLxcARFZWlhBCCEdHR/H555+rtXF1dRUhISEViqescSor5vLs/8n+Hz58KExMTMTIkSOlsoyMDAFAxMXFaYz58fdD+/btxZgxY4QQQuzcuVM8+ae9qmJ62mt///59YWxsLA4dOqS277Fjx4rhw4dLfWh632oyYMAA6TiLBQcHizp16ojs7GypbNWqVaJFixZq7ZKTkwUAcfHixXLti0gbnKkiqkV+/PFHmJqawtDQEL1794a/v7/aKZwn/xWvSevWrdWe29nZITs7GwCQmpqK/Px89OjRo0JxtW/fXu3iaG9vb5w7dw6FhYUAgHPnzmH48OF46aWXYG5uLl1knJ6eXu59iKfMRCmVynL3VZ54yhonOTzev66uLqysrNCqVSuprHh2sDz7/Oyzz7B582akpqY+85jKeu3Pnz+Pu3fv4tVXX4Wpqan02LJli9rp1vK8bwHg3r17MDQ0VCtLTk7GgAEDYGNjI5WlpaWhcePGau2KZznLmvkjqiy96g6AiMqve/fuWLduHZRKJezt7aGnp/4RNjExeWof+vr6as8VCoV02uvJ02ty6devHxwdHbFhwwbY29ujqKgILVu2LPOi5GKNGzeGQqFAamoqBg4cWKI+NTUVNjY2sLS0BADo6OiUSMAePHhQ4XjKGqeylGf/pfX/eFlxolKefXbp0gW+vr6YPXs2Ro8eXSNiAoA7d+4AAH766SfUr19frc7AwED6//K8bwHA2tq6xEX5ycnJmD17tlpZUlISunTpolZW/GOGx5MvIrlxpoqoFjExMUHjxo3RsGHDEgmVHFxcXGBkZITY2NgKbffkRcGHDx+Gi4sLdHV1cePGDZw9exZz5sxBjx490KxZswr9Ws3KygqvvvoqvvjiC9y7d0+tLjMzE99++61aImFjY4OMjAzpeW5uLtLS0qTnlY2nmFKplGbiHve0/VeVRYsW4YcffkBcXNwzjams17558+YwMDBAeno6GjdurPao6DVoANCmTRucPn1aep6Tk4OLFy+iTZs2au2SkpLg5uamVpaSkoIGDRrA2tq6wvslKi8mVUQkMTQ0xAcffICZM2dKp2gOHz6Mr776qszt0tPTERQUhLNnz+K7777D6tWrMXXqVABAnTp1YGVlhfXr1+P8+fP4+eefERQUVKG41qxZg/z8fPj6+uLXX3/F5cuXER0djVdffRUvv/wygoODpbavvPIKvv76a/z22284efIkAgICoKurK9XLEQ/w6Bd18fHxuHjxIq5fvy7N3jxt/1WlVatWGDFiBFatWlWiripjKuu1NzMzw4wZMzB9+nRs3rwZFy5cQEJCAlavXo3NmzdXeF++vr44deqUlASfOHECenp6aqcoL126hFu3bpVIqn777Tf07NlT+wMlKgee/iMiNXPnzoWenh6Cg4Nx9epV2NnZ4e233y5zm1GjRuHevXvw9PSErq4upk6digkTJgB4dOpp27ZtmDJlClq2bIkmTZpg1apV6NatW7ljcnFxwdGjRzFv3jy8/vrryM7OhhACgwYNwtdffw1jY2Op7ezZs5GWloa+ffvCwsICCxcuVJuVkSMeAJgxYwYCAgLQvHlz3Lt3D2lpaXBycnrq/qvSggULEBERUaK8KmMq67UHgIULF8LGxgahoaH466+/YGlpibZt2+LDDz+s8L5atWqFtm3b4vvvv8dbb72F5ORkNGnSRO06q8TERFhaWqotDnr//n1ERUUhOjq6UsdK9DQK8bSrP4mIytCtWze4ubmVWAepqoWEhGD58uWIiYlRW7eInm8//fQT3n//faSkpEBHp3wnW9atW4edO3di3759VRwdveg4U0VEtdL8+fPh5OSEw4cPw9PTs9xfsFS79enTB+fOncOVK1fKfV2Wvr4+Vq9eXcWREXGmiogqqbpmqoiIahomVUREREQy4Hw5ERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQyYVBERERHJgEkVERERkQz+Dw8bZeFdHFzQAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "states_sp, transition_rates_sp = state.get_spontaneous_transition_rates(unit=\"1/mus\", method=method)\n",
    "print(f\"Number of possible spontaneous decay transitions: {len(transition_rates_sp)}\")\n",
    "\n",
    "states_bbr, transition_rates_bbr = state.get_black_body_transition_rates(T_in_K, \"kelvin\", unit=\"1/mus\", method=method)\n",
    "print(f\"Number of considered BBR transitions: {len(transition_rates_bbr)}\")\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "\n",
    "n_list = np.arange(0, np.max([s.n for s in states_bbr]) + 1)\n",
    "rates_summed = {}\n",
    "for key, states, rates in [(\"BBR\", states_bbr, transition_rates_bbr), (\"SP\", states_sp, transition_rates_sp)]:\n",
    "    rates_summed[key] = np.zeros(len(n_list))\n",
    "    for i, s in enumerate(states):\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()"
   ]
  }
 ],
 "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": 2
}