{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Dipole elements decay" ] }, { "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": [], "source": [ "state = RydbergState(\"Rb\", n=60, l=0, j=0.5, m=0.5)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of possible spontaneous decay transitions: 275\n", "Number of considered BBR transitions: 430\n" ] } ], "source": [ "states_sp, transition_rates_sp = state.get_spontaneous_transition_rates(unit=\"1/mus\")\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(300, \"kelvin\", unit=\"1/mus\")\n", "print(f\"Number of considered BBR transitions: {len(transition_rates_bbr)}\")" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "relevant_states, energy_differences, electric_dipole_moments = state._get_list_of_dipole_coupled_states( # noqa: SLF001\n", " 1, 120, only_smaller_energy=False\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAG0CAYAAADJpthQAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAg9RJREFUeJzt3Xd4VGX2wPHv9PROKgmBJPQSigEU6UpRENmffTVYV1dUiBVdsYtlBQRZcW3YVpFdsaCiSEepgdBLIAFCSA/pbTJzf38MGQmhZMJMJsmcz/PcJ5l779w5c9NO3nJelaIoCkIIIYQQwkrt7ACEEEIIIVoaSZCEEEIIIc4iCZIQQgghxFkkQRJCCCGEOIskSEIIIYQQZ5EESQghhBDiLJIgCSGEEEKcRevsAFors9nMyZMn8fb2RqVSOTscIYQQQjSCoiiUlpYSHh6OWn3+diJJkJro5MmTREZGOjsMIYQQQjRBRkYG7du3P+9xSZCayNvbG7DcYB8fHydHI4QQQojGKCkpITIy0vp3/HwkQWqium41Hx8fSZCEEEKIVuZiw2NkkLYQQgghxFkkQRJCCCGEOIskSEIIIYQQZ5ExSEIIIUQzMZlMGI1GZ4fRpul0OjQazSVfRxIkIYQQwsEURSE7O5uioiJnh+IS/Pz8CA0NvaQ6hZIgCSGEEA5WlxwFBwfj4eEhBYYdRFEUKioqyM3NBSAsLKzJ15IESQghhHAgk8lkTY4CAwOdHU6b5+7uDkBubi7BwcFN7m6TQdpCCCGEA9WNOfLw8HByJK6j7l5fyngvSZCEEEKIZiDdas3HHvdaEiQhhBBCiLNIgiSEEEIIcRZJkIQQQghhk+HDhzNt2jQAoqOjmTt3rl2vP2XKFCZNmmTXa9pKZrEJIYQQosm2bt2Kp6ens8OwO0mQhBBCCNFk7dq1c3YIDiFdbEIIIUQzUxSFKqPJKZuiKDbFWl5ezh133IGXlxdhYWG89dZb9Y6f3cWmUql49913GTduHO7u7nTq1In//ve/9Z6ze/duRo4cibu7O4GBgdx3332UlZWdNwaz2cysWbPo2LEj7u7u9OnTp8E17U1akIQQQohmVl1r5oaFG53y2kvuH4ybrvHFEx9//HHWrl3Ld999R3BwME8//TTbt28nPj7+vM959tlnee2113j77bf57LPPuPnmm9m9ezfdunWjvLycMWPGMHjwYLZu3Upubi733HMPU6dOZdGiRee83qxZs/j8889ZuHAhcXFxrFu3jr/+9a+0a9eOYcOG2XgHGkcSJCGEEEKcU1lZGR9++CGff/45o0aNAuCTTz6hffv2F3zeDTfcwD333APASy+9xIoVK5g/fz7/+te/+M9//kNVVRWffvqpdezSO++8w4QJE3j99dcJCQmpd63q6mpeffVVfvvtNwYPHgxAp06d2LBhA++9954kSEIIIURbYdCqWXL/YKe9dmMdOXKEmpoaBg4caN0XEBBAly5dLvi8ukTmzMcpKSkA7N+/nz59+tQb2H3FFVdgNps5ePBggwTp8OHDVFRUcNVVV9XbX1NTQ9++fRv9XmwlCZIQQgjRzFQqlU3dXK6sbmzSjz/+SERERL1jBoPBYa8rg7SFEEIIcU4xMTHodDo2b95s3Xfq1CkOHTp0wedt2rSpweNu3boB0K1bN3bu3El5ebn1+O+//45arT5ny1T37t0xGAwcP36c2NjYeltkZOSlvL0LkhYkIYQQQpyTl5cXd999N48//jiBgYEEBwfzzDPPoFZfuH1lyZIlDBgwgCFDhvDFF1+wZcsWPvzwQwBuu+02nnvuORITE3n++efJy8vjoYce4vbbb2/QvQbg7e3NY489xvTp0zGbzQwZMoTi4mJ+//13fHx8SExMdMh7lwRJCCGEEOf15ptvUlZWxoQJE/D29ubRRx+luLj4gs954YUX+Oqrr/j73/9OWFgYX375Jd27dwfAw8ODX375hUceeYTLLrsMDw8P/vKXvzB79uzzXu+ll16iXbt2zJo1i7S0NPz8/OjXrx9PP/20Xd/rmVSKrQURBAAlJSX4+vpSXFyMj4+Ps8MRQgjRQlVVVZGenk7Hjh1xc3NzdjgOp1KpWLp0qVOXCrnQPW/s328ZgySEEEIIcRZJkIQQQgghzuLSCdKyZcvo0qULcXFxfPDBB84ORwghhGj1FEVxaveavbjsIO3a2lqSkpJYvXo1vr6+9O/fn+uvv57AwEBnhyaEEEIIJ3PZFqQtW7bQo0cPIiIi8PLyYty4cfz666/ODksIIYQQLUCrTZDWrVvHhAkTCA8PR6VS8e233zY4Z8GCBURHR+Pm5sbAgQPZsmWL9djJkyfrVeSMiIggMzOzOUIXQgghRAvXahOk8vJy+vTpw4IFC855fPHixSQlJfHcc8+xfft2+vTpw5gxY8jNzW3mSJ2oqgTyLlztVAgBhZWnWJG2mbVpe/kj/Tib0vIpLK9xdlhCCCdqtWOQxo0bx7hx4857fPbs2dx7773ceeedACxcuJAff/yRjz76iKeeeorw8PB6LUaZmZkkJCSc93rV1dVUV1dbH5eUlNjhXTjYyhfhxFa4fiEEd3N2NEK0WLcufZK8qjNbkNWEafuz7I7nnRWSEMLJWm0L0oXU1NSQnJzM6NGjrfvUajWjR49m48aNACQkJLBnzx4yMzMpKyvj559/ZsyYMee95qxZs/D19bVujlz/xW4qCiwf09Y6Nw4hWrDymkprcmRQe6HXaAAzWbXJVNQYnRucEMJp2mSClJ+fj8lkarCmS0hICNnZ2QBotVreeustRowYQXx8PI8++ugFZ7DNmDGD4uJi65aRkeHQ92BXxzc6OwIhWqwD+ccBUCvubLjjK9bd/j80Ki1g5lBelnODE6IFyMvL44EHHiAqKgqDwUBoaChjxozh999/ByA6OhqVSoVKpcLT05N+/fqxZMkSJ0d96VptF5s9TJw4kYkTJzbqXIPBgMFgcHBEDnLqKJRmg3eosyMRosVJLTgBgJc2CLVahRodXlp/io15HCrIJD4iyskRCuFcf/nLX6ipqeGTTz6hU6dO5OTksHLlSgoKCqznvPjii9x7772UlJTw1ltvcdNNNxEREcHll1/uxMgvTZtMkIKCgtBoNOTk5NTbn5OTQ2ioiyYJxzdCj+udHYUQLU5akSVBCjAEW/f56YMoNuZxtOiks8ISokUoKipi/fr1rFmzhmHDhgHQoUOHBmN2vb29CQ0NJTQ0lAULFvD555/zww8/tOoEqU12sen1evr378/KlSut+8xmMytXrmTw4MFOjMyJjm9ydgRCtEgnSizdaCEef/7z1M6jneVYabZTYhIuQFHAWOmczYY16r28vPDy8uLbb7+tN1HpQrRaLTqdjpqa1j0TtNW2IJWVlXH48GHr4/T0dFJSUggICCAqKoqkpCQSExMZMGAACQkJzJ07l/LycuusNpeTuR2MVaBr+ytJC2GLnApLEtTeO9y6L9wrFPIgtzzPWWGJtq62Cj4a65zXvms56NwbdapWq2XRokXce++9LFy4kH79+jFs2DBuvvlmevfu3eD8mpoa3nrrLYqLixk5cqS9I29WrbYFadu2bfTt25e+ffsCkJSURN++fZk5cyYAN910E//85z+ZOXMm8fHxpKSksHz58gYDt12GqQZO7nB2FEK0OKeqLbXROvm1t+6L9rW0JhVUS4IkxF/+8hdOnjzJ999/z9ixY1mzZg39+vVj0aJF1nOefPJJvLy88PDw4PXXX+e1117jmmuucV7QdtBqW5CGDx+OcpFmwqlTpzJ16tRmiqgF820PxScs45A6uGgXoxDnUGGsoMJUCkBc4J+lO2IDLVX2S40FKIqCSqVySnyiDdO6WVpynPXaNnJzc+Oqq67iqquu4tlnn+Wee+7hueeeY8qUKQA8/vjjTJkyBS8vL0JCQtrEz0yrTZCEDaIGw+4llnFIigJt4BtXCHvILD2J0WRGrXgS5e9n3d+lnSVBMlJKYUUlgZ4eTopQtFkqVaO7uVqi7t2711viKygoiNjYWOcF5ACttotN2CC8L2j0UJYDp9KdHY0QLUZqoWUGm04JINBTb90f6O6HTm0p63Eg74RTYhOiJSgoKGDkyJF8/vnn7Nq1i/T0dJYsWcIbb7zBdddd5+zwHOqSWpDKysowm8319vn4+FxSQMIBtG4Q0c/SgnR8EwR0cnZEQrQIqQWWIpE+unao1X+2rKpUKry1gRTWnORwwUmuiO7srBCFcCovLy8GDhzInDlzOHLkCEajkcjISO69916efvppZ4fnUDYnSOnp6UydOpU1a9ZQVVVl3V/XT28ymewaoLCTyIGnE6SNEH+rs6MRokXIOD3Fv51bw8kbAYYgCmtOcqxYqmkL12UwGJg1axazZs067zlHjx5tvoCakc0J0l//+lcUReGjjz5qMwOxXELUYPj9bcjeA2V54NXO2REJ4XQnyy3JT7hXeINjwR7BHC6FTKmFJIRLsjlB2rlzJ8nJyXTp0sUR8QhH8QmDsN6QtQsOLIMBLloPSojTFEUhv8pSbT/Kt2GCFOEdAjmQV5nb3KEJIVoAmwdpX3bZZa1roVbxp+6TLB8P/Ahm6QoVrq3UWEpFbTkAsQERDY539LMkTaeq85s1LiFEy2BzC9IHH3zA/fffT2ZmJj179kSn09U7fq7KmqKF6DgU3P2hPA+O/W55LISLyirLotakoDH7EOHbcHJJ5yBL4chy0ynMZqXeIG4hRNtnc4KUl5fHkSNH6i3ZoVKpZJB2a6DRQddrYMfnsPdbSZCESztZlkWtyYxe8SfYx9DgeFyQpQXJRCWZJUVE+vk3d4hCCCeyOUG666676Nu3L19++aUM0m6Nul4LKV9AZjIUZYBf5MWfI0QbdOTUCRRAqwQQ6NkwQfLSe+Cm8aLKVMaBvBOSIAnhYmxOkI4dO8b333/f5ipmugyfMIgcZJnuv/97GPygsyMSwimOFmUCEGAIRnOe7jMfXSBVpjLSCk8CvZoxOiGEs9k8SHvkyJHs3LnTEbGI5tJjkuXjwZ+httqpoQjhLCdKLVP8g91Dz3tOkJulHMaxYpnqL4SrsbkFacKECUyfPp3du3fTq1evBoO0J06caLfghIO0TwDvMCjNgiOroMs4Z0ckRLNSFIXcCkvSE+nTcIp/nVDPEPYVQVZZTjNFJoRoKWxOkO6//34AXnzxxQbHZJB2K6FWQ7cJsOXfsPu/EDfGsk8IF1FUXURlbRWgItrv/AlShE8IZEJBVV7zBSeEaBFs/qtoNpvPu0ly1Ip0uxZ0HlBwGNLXOjsaIZpVVnkWRrMZjdmXcF/P857X6XTyVFQjCZJwTVOmTEGlUlm3wMBAxo4dy65du6znnHlcq9USFRVFUlIS1dV/DuFYtGiR9Ry1Wk1YWBg33XQTx48fd8bbahRpNnBVbr7Q+0bL59s+ksKRwqVklVtqIGmVAEJ83M57XucgSwHJCvMpamrlZ0S4prFjx5KVlUVWVhYrV65Eq9Vy7bXX1jvn448/Jisri/T0dP71r3/x2Wef8fLLL9c7x8fHh6ysLDIzM/nf//7HwYMHueGGG5rzrdjE5i62c3WtnWnmzJlNDkY0s943wt5voOg4pP4qY5GEyzhZaqmB5GH2J9i74RT/Op0CQlGhQqGWtMJcugaHNWOUQrQMBoOB0FDLZIbQ0FCeeuoprrzySvLy8mjXzjKRwc/Pz3pOZGQk1113Hdu3b693HZVKZT0nLCyMu+++m4cffpiSkhJ8fBoWa3U2mxOkpUuX1ntsNBpJT09Hq9USExMjCVJroveE+Ntg07uQvAhiRoFW7+yohHC4AwWHUQCdEkig1/kTJL1Gj4fWl/LaIg7mZ0qCJOxGURSqTc6ZRWzQGJpcw7CsrIzPP/+c2NhYAgMDz3nOoUOHWLVqFVOmTDnvdXJzc1m6dCkajQaNRtOkWBzN5gRpx44dDfaVlJQwZcoUrr/+ersEJZpRj+th19dQmm1ZxLbnZGdHJIRDFVUVsbdgHwDhbt3OWwOpjq8uiPLaItJPnWyO8ISLqDZVk7g80Smv/cnYT3DTnr9r+WzLli3Dy8sLgPLycsLCwli2bBnqMyb33HLLLWg0Gmpra6murubaa69lxowZ9a5TXFyMl5cXiqJQUVEBwMMPP4yn5/nHATqTXcYg+fj48MILL/Dss8/a43KiOWkN0O92y+c7PgNjlXPjEcLBNmVvosZkQmcKo73P+Wsg1WnnbulCOFEiU/2FaxoxYgQpKSmkpKSwZcsWxowZw7hx4zh27Jj1nDlz5pCSksLOnTtZtmwZhw4d4vbbb693HW9vb1JSUti2bRtvvfUW/fr145VXXmnut9NoNrcgnU9xcTHFxcX2upxoTl2vhZ2LLXWRdi/5M2ESog36PfN3ak1m3Gu7E+x98f+iw7yC2VkIJ8ukWKSwH4PGwCdjP3Haa9vC09Oz3uoZH3zwAb6+vrz//vvWgdihoaHWc7p06UJpaSm33HILL7/8snW/Wq22ft6tWzeOHDnCAw88wGeffWaPt2V3NidI8+bNq/dYURSysrL47LPPGDdOBvm2ShodXHYPrHoJtn8KMSPBN8LZUQlhd3kVeRw6dYhak4K7qdsFZ7DV6RUSw/LjcKzsCGazgvoiXXJCNIZKpbKpm6slqZuqX1lZed5z6sYVXeicp556ipiYGKZPn06/fv3sHuelsjlBmjNnTr3HarWadu3akZiY2KC/UbQisaPg4E+WRWzXvwXXvAWyELFoYzZmbQRAa4pCo3jT3t/9os8ZHzeQt7apqVLy+P1YKld27OzoMIVoUaqrq8nOtrSgnjp1infeeYeysjImTJhgPaeoqIjs7GzMZjOpqam8+OKLdO7cmW7dup33upGRkVx//fXMnDmTZcuWOfx92MrmBCk9Pd0RcQhnU6ngyiRYcqclSUpdAZ2vdnZUQtiVpXtNobYsDjcV9In0u+hzfNy8iPLszNGyA/x06A9JkITLWb58OWFhlhmc3t7edO3alSVLljB8+HDrOXfeeSfw51T+oUOH8uqrr6LVXjjNmD59OoMHD2bLli0kJCQ47D00hd3GIIk2wLc99LsDtn4AG9+BqIGWgpJCtAEny05ytOQoVUYFt9ouxAZ74euuu/gTgcHhCRw9dIDtecnAFIfGKURLsmjRIhYtWnTBcxRFueh1pkyZcs5p/4MGDWrU853B5gSpqqqK+fPns3r1anJzczGbzfWOn10YSjhTE77p+twMh3+DU0dh00IY/qTdoxLCGX4/+TsA7uaOVONB3w7+jX7upG5D+PLQp+TXpHP8VAFR/ueu/yKEaDtsTpDuvvtufv31V/7v//6PhISEJhebEi2URgdDH4fvHrSMSYoZAZEtq9lTCFspisLGkxtBgaqSWFRAvyi/Rj8/NjCCAF0EhcZMvtm/gWmXX+ewWIUQLYPNCdKyZcv46aefuOKKKxwRj2gJQntaCkjuXQqrX4X/+wg8ApwdlRBNdqzkGJllmZjMapTKGDz0GrqEeNt0jT7t+rH6ZCZ/ZG5hGpIgCdHW2VwoMiIiAm9v236xiFZo0N8hMAYqT8Gql+GsrlQhWpMf0n4AwF/dBTUG+kb6odXY9utvfJzln8KjZfuorHHOEhFCiOZjc4L01ltv8eSTT9aroCnaIK0eRj0HWjfLrLad/3F2REI0SVpRGhsyNwCgq7R0F/e1oXutzrDoHhhUPpiUGn5K3WLPEIWLaKmDkdsie9xrmxOkAQMGUFVVRadOnfD29iYgIKDeJtoQ/w4wZJrl860fQvYep4YjhK0UReHz/Z8DkBByOZm5lhXD+0U1foB2HY1GTWffPgCsSPvDfkGKNk+ns8yWrFt/TDhe3b2uu/dNYfMYpFtuuYXMzExeffVVQkJCZJB2W9d57J91kVa+ANe/J+ORRKuxM28newv2olVr6eYxhjVKPhF+7gQ3ooL2uQzvMIjdRevZd2onZsWMWmWX5SxFG6fRaPDz8yM3NxcADw8P+dvpIHUL4ebm5uLn52et6N0UNidIf/zxBxs3bqRPnz5NflHRzC7l51ClgiFJkHcAijLg13/AtXMsi9wK0YKZFbO19Whs9FjSsi2/KPvbML3/bBO7JrAgRU95bTHbMveT0L6HXWIVbV9oqGVh5LokSTiWn5+f9Z43lc0JUteuXS+4topog/QeMGYWfPsA5OyFNa/BqJmyFIlo0dZmrCWjNAMvnReTYibxyJYDAPTr4NfkawZ4ehDu3o0TVTv5cu9ySZBEo6lUKsLCwggODsZoNDo7nDZNp9NdUstRHZsTpNdee41HH32UV155hV69ejXo3/Px8bnkoEQL5BcJV78EPz4KR1aBXxQMuNPZUQlxTpW1lXx96GsAJsVOoqhcQ15pNTqNih7hl1Yd/ppO43hv3042Za8nrzyRdp5B9ghZuAiNRmOXP97C8WzuQB87diwbN25k1KhRBAcH4+/vj7+/P35+fvj7N73pWrQC4X3hykctnycvsoxLEqKFURSFf+/6N4VVhQS7BzM2eiwb0woA6B7ug5vu0v443dH/cjzpQI3JyL+2fW2PkIUQLZDNLUirV692RByiteh6jWUs0s4vYc0s0HtBh8HOjkoIq9+O/8YfJ/9Ao9Iwte9U1CotP+/OAmB45+BLvr6bTsP46IksOTqf346vZFr1bfgaZM1CIdoamxOkYcOGOSIO0Zok3AfluXB4JayYCeNeg4j+zo5KCNKK01i0dxEAt3S9hS4BXfjjcD75ZTX4uGsZ2rmdXV7n3oQRfJ/+XyqMWXycspRpA6fY5bpCiJajUQnSrl276NmzJ2q1ml27dl3w3N69e9slMNGCqdUw/GmorYajG+CXZ2D8Py1LlAjhJOXGcuYkz6HWXMuAkAFc2+laAH7YZWk9GtMjFL3WPtPyA70MXBk6nt9yPuSHIz9zX78b8dB52OXaQoiWoVEJUnx8PNnZ2QQHBxMfH49KpTpnlUqVSoXJZLJ7kKIF0mgtlbZ/eRpObIWfn4RrZ0O7Ls6OTLggk9nEgpQF5FbkEuwezAN9HkClUnE0v5w9mcWoVTCuZ5hdX/P+gVex5tvvKKrK5+v9y5jS+0a7Xl8I4VyN+ncqPT2ddu3aWT9PS0sjPT29wZaWlubQYEULo9XD1S9DWG+oKYNl0yHrwi2MQtibWTHzr53/IjknGa1Ky7T+0/DSewHw4+mxR4M6BdLO2761uzoGeRPvNxqAxfu/o7JWyp8I0ZY0KkHq0KGDtepnhw4dzrtFRUU5NFjRAuncYOzrENYHasrhp8cgQ9apEs3DrJh5b9d7bMjcgBo10/pPI8YvBoDSKiOrD1iK8k3oE+6Q1/9bwng0Zn/yK4pYtFvWKxSiLbG5Q37KlCmUl5c32H/06FGGDh1ql6BEK6P3gPFvQtRgy7ik5TMgba2zoxJtnKIofLTnI9ZkrEGNmof6PsRloZdZj/+2P4fqWjMdAj3oEe6Y+mx9IwPo7nEdZgW+3v8DBwsPOuR1hBDNz+YEaefOnfTu3ZuNGzda933yySf06dOHoCApmOaytAZLIcmYEWCuhd+ehz3fODsq0UbVmmv5YPcHrDi2AhUqHoh/gMsjLrceN5sVftyVDcC1vcMdtu6VSqXiieFj8DD2oqTKyJub38FokirJQrQFNidIW7ZsYfLkyQwfPpynn36aG2+8kalTp/LPf/6TpUuXOiJG0VpodDByJnS7FhQz/P42/DEfzGZnRybakHJjOa9teY3fjv+GChX39r6Xoe3rt16vPphLTkkVXgYtw7vYZ2r/+XQP92Fix1tQK57szTvKkkP/dejrCSGah811kHQ6HW+++SYeHh689NJLaLVa1q5dy+DBUixQYCkBcOVj4B0GW96H3f+FkiwY+Q9LV5wQlyC7PJs3tr5BZlkmBo2hQbcaQGWNiU82HgPg//q3v+TK2Y1x35XdWff5OE7W/pdPd/+Xy8MHE+0b7fDXFUI4js0tSEajkUcffZTXX3+dGTNmMHjwYCZPnsxPP/3kiPhEa6RSQd+/wujnQaOHY7/D9w9ZEiUhmiglN4VnNjxDZlkmAW4BvHj5iw2SI4D/bj/BqfIaQnzcHDY4+2y+7jqmXj4Ot9ouFJRXM2fbfKpN1c3y2kIIx7A5QRowYADff/89a9as4ZVXXmHNmjVMmzaNyZMn8/e//90RMYrWKmYETJgL7v5QcBi+uReOb3J2VKKVqTHV8NGej5i1ZRZlxjJi/WJ5dcir52yhyS2tYun2EwDcNSTaboUhG+OqbiEM8JsMZjdSsg/zbsq756wXJ4RoHZqUIKWkpDBo0CDAMkjxySefZOPGjaxbt87uAYpWLqQHTP43BHeH6lJY/hRs+1jGJYlGOVp8lKfWP8UvR38BYHzH8Tw3+Dn83c69MPYnfxzFaFLoGeHD4E6BzRkqarWKaSP7Elg9mbJqM8uPrOPbw982awxCCPuxOUH68MMP8fT0bLC/b9++JCcn2yUo0cZ4BcOEt6HHJFAUSF4EPz8BFYXOjky0UFW1VXyx/wue3vA0mWWZ+Bn8mJEwg8Qeieg1+nM+Z39WCesO5aNSwT1XdnLYzLUL6Rjkyf2Dh+FbfRX5ZdV8suc/bMve1uxxCCEu3SW1P1dVVVFSUmLdqqulz12ch1YPQ6bDiGcsJQFObIX/3gnHNl78ucJlKIrC5qzNJK1J4vsj32NSTCSEJvDG0DeID44/7/NMZoX311kq+Y/uFkJMO69mirih6+LDuSr6KtyN/cgqrmRu8jwySjOcFo8QomlsTpDKy8uZOnUqwcHBeHp64u/vX28T4oI6Xw3XvweBMVBZZOly+/1tqK1xdmTCydKL03l186vMTp5NQVUB7dzb8fiAx3l0wKP4Gnwv+NzFWzNIzS3DXa/h9kEdminic1OpVEwb1ZkeHhNQGyM5WljEK5teJbs826lxCSFsY3OC9MQTT7Bq1SreffddDAYDH3zwAS+88ALh4eF8+umnjohRtDUBHWHSQuh1g+Xxnm/gf3dDzl7nxiWcIqssi7e3v81T659iV/4utGot18dez1vD32JA6ICLPn9/VgmLtx4H4MERsfh7nrsLrjm56zU8c01Pwkx/wVjtT2p+Fi9ufJHcilxnhyaEaCSVYuM0i6ioKD799FOGDx+Oj48P27dvJzY2ls8++4wvv/zSZab7l5SU4OvrS3FxMT4+jlnG4JItmQKF6XDtbIjo7+xozu34Zlj7mmU8kkoNvW+CAXdauuFEm5ZZlsn3R75nXcY6zJhRoeLy8Mu5scuNhHqGNuoalTUmHvpyBzklVQzv0o5Hr+7i4Khtsz41j1nLt1Pg/gW+3qV0DorgucHP0c7DscUrhRDn19i/3za3IBUWFtKpUycAfHx8KCy0DLQdMmRIq5rFlpGRwfDhw+nevTu9e/dmyZIlzg7JNUUNhBs+gbirLdW3d35paU06meLsyISDHD51mLe2vcWjax5lTcYazJjpF9yP1658jYf7Pdzo5AjgvXVHyCmpItjbwP3DYhwXdBNdGdeOv13Zk8CqWyku9SK14KS0JAnRSthcSbtTp06kp6cTFRVF165d+frrr0lISOCHH37Az8/PASE6hlarZe7cucTHx5OdnU3//v0ZP378OWfoCQdz84GRz0Cn4bD+LSjKgB8egS7jYODfLHWURKtmNBnZlLWJX4/9yqFTh6z7Lwu5jImxE+ns39nma65PzWPl/lzUKph+VWc8DTb/OmsW18VHUFFj4rMtt5DPfzisOsmzvz/Lk5c9SSe/Ts4OTwhxHjb/RrnzzjvZuXMnw4YN46mnnmLChAm88847GI1GZs+e7YgYHSIsLIywsDAAQkNDCQoKorCwUBIkZ4q+AkJ7wdb3Yf8PcPBnOLoBBt4PXcZbljERrUp2eTZrMtaw6vgqimuKAdCqtAyJGMKEmAm0927fpOseySvj7d9SActyIj0jLjyI29luviySihoT/0u5lUKWAPm8sPEFHun3CP1C+jk7PCHEOdg8Bulsx44dIzk5mdjYWHr37m2vuFi3bh1vvvkmycnJZGVlsXTpUiZNmlTvnAULFvDmm2+SnZ1Nnz59mD9/PgkJCTa/VnJyMomJiezZs6fRz5ExSA6WsxfWz7ZU4AYI6gyXPwRh9vseE45RWVvJppObWHtiLfsL91v3B7gFMDpqNKOiRuHn5tfk6xeUVfPokp0UlNXQN8qP5yb0QKNu/ppHtlIUhX+tOcJPe45xyu0bvHxOEujpxp097+Tq6KudHZ4QLqOxf78vuU26Q4cOdOhg/2m15eXl9OnTh7vuuovJkyc3OL548WKSkpJYuHAhAwcOZO7cuYwZM4aDBw8SHBwMQHx8PLW1tQ2e++uvvxIeblmjqbCwkDvuuIP333//gvFUV1fXq/NUUlJyKW+vebTmZQ7qKnDv+QaSP4b8Q5b13GJGWlqUvEOcHaE4g9FkZEfuDv44+QfJOcnUmC1lG1So6N2uNyOjRjIgZABa9aX9yqkymnj5x/0UlNUQGeDOE2O7torkCCzT/x8YFoOHXsP/tt9IsXk5Oeb9fLj7Q46WHOXOHnei0+icHaYQ4rQmtSBt3bqV1atXk5ubi/msJSMc0c2mUqkatCANHDiQyy67jHfeeQcAs9lMZGQkDz30EE899VSjrltdXc1VV13Fvffey+23337Bc59//nleeOGFBvtbdAvS14lw6mjrbEE6U0UhbPsQDvxoSfo0eug5GeJvs4xfEk5RWVvJrrxdbMneQnJOMpW1ldZjYZ5hDGs/jKHthxLobp8lP8xmhdeXH+CPIwV4u2mZfWM8ob5udrl2c/suJZMPNqRRqt2I4v0HIT4Gon07kNQ/yaZB6kII2zmsBenVV1/lH//4B126dCEkJKReOf/mKu1fU1NDcnIyM2bMsO5Tq9WMHj2ajRsbV5lZURSmTJnCyJEjL5ocAcyYMYOkpCTr45KSEiIjI20PXtjOIwCGPg7dJ8Ef8yFrJ+z8ypIwxd9mSZakLECzyK/MJyU3he2529mVtwuj2Wg9FuAWwBXhV3B5+OV09O1o198HiqLw7/Vp/HGkAK1GxTPXdGu1yRFYBm77e+iZ85ua8qJwskw/UGtKZ8b6Gdzf534Ghg10dohCuDybE6S3336bjz76iClTpjggnMbJz8/HZDIRElK/myUkJIQDBw406hq///47ixcvpnfv3nz77bcAfPbZZ/Tq1euc5xsMBgwG+SPsVEFxljXdMjbD5vegMA02L4TdS6DvX6HrtZYlTYTdGE1GDhQeYFf+LlJyUzheerze8WCPYBJCE0gITSDOPw61yv4D6RVF4cMN6fy4KwuVCh4eFUeP8JY9KLsxhnZuR4CnnteXaykoSSTX+B01fvnMTp7NsPbDSOyRiKdOJo0I4Sw2J0hqtZorrrjCEbE0qyFDhjToHhStgEoFUYOgfQIcXgHbPoLSbMtyJSn/gX63Q+dxkig1kclsIr04nX2F+9idt5v9hfvrtRKpURPrH0vf4L4MCBlApHekQ1uOFUXh043H+C7lJAB/Hx7LiC7BDnu95tYzwpe5N8Xz+vID7Mu6lVLjOqr8trMmYy2783fz9z5/p1e7c//TJoRwLJsTpOnTp7NgwQLmzp3rgHAaJygoCI1GQ05OTr39OTk5hIZK/71LUKuh8xjLoO2DP8P2T6E8zzLzbfun0Ptm6HYt6NydHWmLZjQZOVx0mIOnDrK/YD8HCg9QZaqqd46/wZ+eQT3p064P8cHxeOu9my2+L7dk8N/kEwDcPyyGsT3b3s93oJeBV6/vxce/H+X7nSOoKYgjr+ZnjKZ8Xt78MqOjRnNL11vw0jtvAV4hXJHNCdJjjz3GNddcQ0xMDN27d0enqz/r4ptvvrFbcOej1+vp378/K1eutA7cNpvNrFy5kqlTpzr89UULotFB94nQeSwc+AFSvrQkShvfgR2fQa//s4xdksHcABRUFpBalMqhwkOkFqWSVpxGrbn+TE8vnRddA7rSPbA7vdv1pr1X+2YbX1jHbFb4+I+jfLsjE4B7ruzINb3DmjWG5qTVqLl3aCd6hPvwzmotxcWJFFatweC7h9+O/caW7C3c3v12roy4stm/FkK4KpsTpIcffpjVq1czYsQIAgMDHfbDWlZWxuHDh62P09PTSUlJISAggKioKJKSkkhMTGTAgAEkJCQwd+5cysvLufPOOx0ST+vmAr9QtXro+RfoOgFSf7F0t5WchK0fwo4vLFW5e98IPuHOjrTZlNWUcaT4CGlFaRwuOsyR4iOcqjrV4Dw/gx9dArrQxb8LPQJ7EOUT5ZCxRI1VazIzb2Uqqw/mAXD3kI5cFx/htHia0+WxQXQP92H+qsNsSb+a6vyu5Pv8RrVXEQtSFrAmYw139riTSB+ZICKEo9k8zd/b25uvvvqKa665xlExAbBmzRpGjBjRYH9iYiKLFi0C4J133rEWioyPj2fevHkMHNg8sz9aRaFI6zT/ORDhYtV6zSY4stqytltdsUmVGqKHWBKpsD6W8UxtgKIoFFYVcrz0OEeLj5JenE5acRp5lXkNzlWjJsonijj/ODr7d6aLfxeCPYJbTKtEZY2J137ez/bjRahV8MjoOEZ2db2aV4qisHJ/Lv9en0ZFTQ0Vus1o/Lbg56FGq9YwqsMobuh8A76G1j9YXYjm1ti/3zYnSB06dOCXX36ha9eulxxkayYJUiuhKJC5HXYttsx+qxMYY0mUYke3qhIBRpORE2UnOF5ynGMlxyxJUclRSmtKz3l+iEcIMX4xls03ho6+HXHTtszp8fll1bz6435Sc8swaNXMGN+V/h0CnB2WUxWW1/DRhnTWHsqjVlVEjeca9F5p+LjpcNO6cX3c9YzrOA6DpvV8DwvhbA5LkD7++GOWL1/Oxx9/jIeHxyUH2lpJgtQKFabDnv9B6q9Qe7oqusHbMn6p2wTwt39F+KYymU1kV2RzovQEx0uPk1GawYnSE2SVZWGm4exLNWrCvcLp4NOBTn6d6OjbkWif6FYzTfxAdgmv/Lifogoj3m5anpvQgy6hzTcYvKXbmVHEu2uOkFlUSbX6OCbvtXh4FuBh0OBv8GdS3CRGRY6SStxCNILDEqS+ffty5MgRFEUhOjq6wSDt7du3Ny3iVkYSpFasqgQO/gR7v4XSrD/3h8dbail1HNpsrUq15lpyKnI4UXqCzLJMMkozyCzNJLM8s8Hg6TpeOi+ifKLo4N2BSJ9Ion2iifSORK9pnaUNftuXw4I1h6k1KUQFevDsNd1bdRFIR6mpNfPzniy+2pJBaXUNldp9qH024uFRjptOQzv3dkyOm8yV7a9Ep5ZESYjzcVgl7bMXjBWi1XHzgT43Q68b4cRW2PcdHN8IJ1Msm8EbYkdBl2ssxSntMD7HaDJysvwkmaWZnCg7wYmyE2SWZpJdnk2tcu5EyKAx0N6rPe292xPpHWndAtwCWsyYoUtRU2vmo98tBSABBscEMn10Z9z1GidH1jLptWqui49gZNdglmw7wQ+7NBgLu1FakkKZzxaqjNm8t+s9/nvov0yMncjIyJGtNmkWoiVo0lpsQlqQ2pyyXEur0sGfLYUn6wR0hLgxEHcVeAZd9DI1phpLAnS6RajuY055zjm7xgDcNG5EeEXQ3tuSDEV4RRDpHUmQe5BTZ5M50smiSt5YfoAjeeUA3JIQxc2XRaJuJQvPtgQ5JVV8vTWD3w7kYjLXUK7dgcpnG57uNXjoNfjqfRnfaTyjo0ZLDSUhzuCwLjZhIQlSG2U2w8ntlnXejm4Ak2VVelRqyz2MHQ3RV1KrcyOrPIuMkgzrGKGM0gxyK3JROPePlKfOkwivCGsCVJcUBbo5rlxGS7TuUB7vrDpMpdGEt5uWpKs6MyDatQdjX4rs4iq+3pbBygO5mMxGKrQ7MXtuxcOjEi+DFnetGyMiRzC+03iCPdpOFXIhmkoSJAeTBMkFVJfCkdWUH/qRo7l7OKoycUxt4qgaMg3u1Bq8QO/VoAvOW+dtbQ1q792e9l6WViE/g59LJUJnK6+u5d/r0lh1IBeAHuE+PDamC0FeMgPLHnJKqvh2RyYr9uVQVWukUrsPo/tW3NxP4eOuQ6fW0D+kP2Oix9AzqKdLfy8K1+awMUhCtGWVtZWkF6dbCisWHSG9OJ0cbRGEhEB1iWWAt6kGqo24VZcSiY4ov05Ehl1GZPRwIv1j8dH7yB+fs+zMKGLub4fIL6tBpYIb+rfn1oEd0LSULrXaGkBpVSUfzhbi48bfhsVwy8AoftqVxbJdbhSV9qSm4ijZ+i1oPY6zvmYTW3O2EuEVwdXRVzM0YigeOtedjSzEhUgLUhNJC1LbkF+Zz4HCAxwoPMChwkNklGacc6xQsHswHXw6WDZFS4e8I7Q7vhV12RnjlbQGaH8ZdBxmWVBXljehymji041H+WGnZSB2iI8b06+Ko0d4CypwqCjww8OWMhC3LgZ96yiNcDE1tWZ+P5zPsl1ZHMopxajKo0K3A5PbXrzcFbzdtHjo3BgcPphRUaOI84uTxF64hGZpQarLreSHSrQWRVVF7M7fzZ78PezO301BVUGDcwLcAojziyPGL4ZOvpaaQg0GuXbF8oc17yCkrbFspVmWcUtHN4BaY6nW3eEKy+bTdtcRO5/kY6d4d81hckosNafG9gzlris6trxZatm7IWuX5fPSbEsR0TZAr1UzomswI7oGk5pTyk+7s1mfGkJlyTAqK3ZzQrcDvVsRBeW/sfr4GqJ8IhkeOZwrI66UCt1C0MQWpE8//ZQ333yT1NRUADp37szjjz/O7bffbvcAWyppQWodzIqZtKI0tuVsY3vOdo6VHqt3XI2ajr4d6RLQha4BXYnzjyPArQkDhhUFCo7A0XWQttZy38/kHw1Rgy0tSyE9QdN2e7eLK4x8sCGNNafXUmvnbeDBEbH07+Dv5MjOY/WrcOgXy+f/91GbSZDOpaKmlnWH8vh1bw6HcksxqjMp16VQozuAlxt4GbR46HX0C+nLsPbD6BfcT4pPijbHYS1Is2fP5tlnn2Xq1KlcccUVAGzYsIH777+f/Px8pk+f3vSohbADs2Lm0KlDbMjcwNbsrRRVF1mPqVAR7RtNr6Be9AjsQZeALrhr3S/9RVUqCIq1bAPuguJMOPYHHNtgaZ04ddSy7fzSMrA7oh9EDoTIBPBqGzOLzGaFX/Zm8+nGY5RV16JWwYQ+4dw2sEPLazWqU11qaf1zER56LWN7hjG2ZxjHCspZdaA9aw91Iq+shMqq/eRqd2HWZZFb+gcbMrYQ6OHNwLCBXBlxJd0Cu7XZshNCnIvNLUgdO3bkhRde4I477qi3/5NPPuH5558nPT3drgG2VNKC1PJklWWx5sQafs/8vd5CrW4aN+KD4+kf0p/44Hh89M389aoqsRSkPL7Jsh5cVXH94/4dIGIARPS3VPNuhWNgDmSX8N7aNA7nlgHQMciTh0bGEhfSwpcL2fstbJjz5+M23oJ0Lmazwq7MYtYdyuOPI/kU1eRQqd1NpXYfaEvxMmjxNGgJ8wrk8vDBDA4fTJx/nCRLotVyWAtSVlYWl19+eYP9l19+OVlZWed4hhCOY1bM7MjdwS9Hf2Fn3k7rfjeNGwPDBnJ5+OX0COzh3G4CNx9LZe7YUZY6S3n7IWMLnNgGufvg1DHLtud/lnpLwd0tiVJ4Pwjt2aJnVhWUVfPZpmOs3G+Zuu+h1/DXQR0Y3yus5cxQu5ADy5wdgdOp1SriI/2Ij/TjgeExpGQUsf5QNzal51NUfYzK6r1kaw+QXZzJwbylfH3gB9r7BHNFxCAGhQ2SZEm0WTYnSLGxsXz99dc8/fTT9fYvXryYuLg4uwUmxIUYzUbWZqzlu8PfkVtp+eOsQkV8cDzD2g+jf0j/lrnMgloNIT0s24A7La1LJ7dbkqXM7VCSCTl7LNuOz0GjsyRMYb0hLN7yPJ0dugQvUZXRxNIdmfwv+QTVtZZZf6O6BTPl8mj8PFrgfT+XvEOQn2q5x2otGCudHZHT6TRqLosO4LLoAIymOHad6M6mtAQ2puWQU3mISuMBcqsOkVNynL05mXyu/5ZQr0CGRg4iIewyugd2l3XgRJthc4L0wgsvcNNNN7Fu3TrrGKTff/+dlStX8vXXX9s9QCHOVGuuZf2J9XyT+o01MfLSeTEicgRXdbiKEM8QJ0doIzcf6DTcsgGUZFkSppM7LFt5PmTttGx8ZpkdFxgHob0srUshvcAzsNnCNZsVVh7I5YvNxygos1QZ7xrqzT1XdqJLaAvvTjtbXetR9JWWey0JUj06jZr+Hfzp38GfB4bFkJrbmy3pBWxKz+VQ0R4qaw9yqiaVwvIsDuZ+yxf6H/B192JweH+GRg6iT3AfPHWtr7tYiDpNmsWWnJzMnDlz2L9/PwDdunXj0UcfpW/fvnYPsKWSMUjNLyU3hY/3fEx2haX2kJ/Bj0mxkxgZNRKDpuV2QzWZokDxiT8TpKydUJbT8DzvMAju9mfLVEAMaO3biqMoCpvTC/ls4zGOF1YAEOJjIPHyaIbEBrW+Uh/GSvj8L1BTDtfMhlUvQeUplxyD1BR5pdUkHzvF1qO5bD65i2JlP1XaVMwqy9p6Oo0aT72OLv5dGRE9kCvaX0aYZ1jr+z4RbZJD6yD179+fzz//vMnBCWGLoqoiPt33Kb+f/B0AH70P18Vex1UdrmqbiVEdlQr8Ii1bt2st+0qzIXsP5Oy21O8pTLPUXyrNgiOrLOeotRAYC8FdoV1XCOpsKTOgbtpMsl0nivhs4zEOZJcClqngN17Wnmt6haPXttKxJ0dWWZIjn3AId51/7OylnbeBsT1DGdszFKOpJwezS9l+vJANx/ZwqHgHleYjFJny2Vy5k80nd6LXfkCQWzDxwfFc1XEgCRG9cdO6OfttCHFBNidIGo2GrKwsgoPrT00uKCggODgYk8lkt+CEa1MUhTUZa/hs/2eUG8tRo2Zsx7Hc2OVG+0zNb428Qy1b3GjL45pyyD0AuXshZ5/lY1UJ5B2wbHW0BkvXXFCcJXkKjIWAjhccAL73ZDFfbD7O7hOWWXd6rZpJ8eFc3689XoZWXMfp4M9/zlzreo1lXJhoMp1GTc8IX3pG+HLH4I6UVo1l78kS/jh6mI0nt5FZuY9q5Tgna7M5Wbacn9KWo9foaO8ZS3y7PozqeBmXRXRGo5Gvg2hZbP4td74euerqavT6VjI4U7R4NaYaPtz9IWtOrAGgo29H7u11LzF+0v1Rj94T2ve3bGDplivNtsyOy90P+acHIhsr/hz8XUeltrROBXSydMsFdAL/aPaVefLVthPsOF4EgEatYmzPUG7o357A1rywrNkEm/4Fu/9redzhCuh1g3NjaoO83XQM6hTIoE6BwECKKmpIychl9dHt7MjbQXb1AWpMJaSV7CetZD/fHPkKncqLSI8u9ArqxdAO/RgQ2aF1J+GiTWj0d+C8efMAy7IiH3zwAV5efy69YDKZWLduHV27drV/hMLl5FbkMjt5NunF6ahRc2OXG5kYMxFNE7uIXIpKZVnWxCfMUlYALKUFSk5YZm0VHIaCVEvSVFX8Z4mBI6upNJooLK8Bo5qR6hB6aEMJioilX58++IeroDX3iOQfhk0LLDMFAfpPgX6J0nrUDPw89Azv0p7hXdoDEymtMvL70UOsPbaNPQW7yK46glEpI608mbTyZL47tgitOYAgfQxd/buTEN6bPhHhdAzyxKCV3wGi+TQ6QZozx9IkrSgKCxcuRKP58xtVr9cTHR3NwoUL7R+hcCl78vcwJ3kOZcYyvPXePNL3EXq16+XssFo3tRr8oixbXdecokB5PkrhEY6m7ibt4G60JUcJNuegx0g3fQ5XeJxCV3YAfj8920ulAs924NsefCIsm3eoZRyPdxgYvC3ntBTGKstYo/0/WFrUwFIiYcTT0HGoc2NzYd5uOsZ27cHYrj0AqKypZu3R3azPSGZP/m5yqo5jpJDs2kKy87ayJg+0O4JwM3cg0jOOXkE96BEWSkw7LzoGeeKmk6RJOEajE6S6CtkjRozgm2++wd+/ha6rJFqtXXm7eGPrGxjNRmJ8Y0gakESQe5Czw2qTzAr8ngVLthlIz+8GdEPrp+Kqbu24sbOGIGOWpWWpOAOKjlu26lIoy7VsdS0xZ9K5W5ZN8Qq1JFKeQeARaPncIwDc/cHNz+4z7AAw1UJ5nqWFLGevpSsx7yCYLKUIUGss0/n7T7GMvRIthrvewNjOAxjbeQAAFcYKtmXtZv3x7ezM3UN2xQmqagsoM+ezvzqZ/ZmgzQhAb4rEzRxFlHcsXYIiiGnnTcd2nnQK8mw9tbhEi2ZzJ+/q1asdEYdwcbvzdluTo/4h/Zneb7oskukAVUYTv+3P4dsdJ8kpqQLATadmbM8wJsWHnzHGKAaih/z5REWBqiLLGnMlmZbEqTTbUrepJNMyRd5Y+WeX3YXoPCz1nwzelk3vZUmudO6gdbckUGqtZVOdbh1QTJYxRGYj1FRYBqfXlFm6CUuzLcmRYm74Wt5h0G0CdBlnSdJEi+eh82Bo1ECGRg0EoLSmlP0F+9matZvk7F1klGZQZSyhunYXFeadFNbA7kwv9Mcj0ZvaozO3p51bBNGBXkQHelq2IA8iAzyki07YREbBCafbnbeb17e+/mdy1H+6VOO1s6KKGn7anc2yXScpraoFwNtNy4Q+4VzbOwxvt4vcb5XK0gLk7m8pUHm22uo/W5fKTics5flQUWD5vLLIkkSZay0Dxo0VlsTGnjR68I04XQ+qp+Wjb2TL6vYTNvPWe5MQlkBCWAIAZTVlHDx1kP0F+9mVu4+DhYepNNZQXZtKRe1BjCYzBYqOtLwIvDOHozeHAaBWQaivG1EBHpYt0JNIf3fa+3u03nIVwqEkQRJOdaDwgCRHDnS8oILvUjJZfTAXo8kyAzXEx41JfcMZ3S3EfuM3tIY/azadj6JYWn0qi6C6BKrLTn8stSRYxgqorbJ8bm0xqgVUp1uU1JYWJb3n6c3L0gLlFWJpKXL3l0HXLsBL70X/kP70D+kP3aHaVM3hosMcKDjAoVOHOFB4kOKqcqprswjT7yfU3JWj+eWUVtVysqiKk0VVbEortF5PrbL8TEQGeND+dMIUGeBOhJ/7xf9xEG2aJEjCaYqri5mbPBej2Ui/4H6WbjVJji6Z2aywI+MU36ecZPvpqfoAcSFeXN83gstjgpyzkKxK9We3mhB2YtAY6BHYgx6BlkHfZsXMN6nfsOTQEnqH+TCtfy8URaGowsixwgqOFZSTUVjB8dNbebWJrOIqsoqr2JJe/9p+Hjoi/Nxp7+9OuJ8laQr3cyfU1w2d1G1q8yRBEk5hVszM2z6PU9WnaO/Vnkf6PSJjji5RZY2JlQdy+GHnSU4WWcYXqVUwqFMgk/pG0DXUW5Z6EG2eWqXGS+9Vb59KpcLfU4+/p574SD/r/jMTpxOnKjhxqpITpyrIKKyksLyGogojRRVG9p4sOes1LNXEw08nTGG+btaPIT6SPLUVjUqQdu3a1egL9u7du8nBCNex5OAS9hTswU3jxvT+02XZgUuQWVTJT7uyWLE/h8oaSyV7d72Gq7uHcG3vcEJ95d4KcS7nS5zA8g9HZlFd0lTJyaK6rYpKo4mckmpySqqtBVXrqFUQ5GUgxNeNMB83Qn1Pb6c/9zJo5R+VVqJRCVJ8fDwqlQpFUS76hZWlRsTFpOSm8M3hbwC4t/e9tPdu7+SIWh+zWSH5+Cl+3JVF8rFT1v3hfm5M6BPOqK4huOtlxo4QTeWu1xAb7E1scP0uYUVROFVhPCNhqrR20WUVV1JlNJNbWk1uaTW7KW5wXQ+9hhAfN0J8DKc/uhHs/efn8nPbcjQqQaqrgQSwY8cOHnvsMR5//HEGDx4MwMaNG3nrrbd44403HBOlaDOKqop4Z8c7AFzV4SqGRAy5yDPEmYorjfy2L4ef92SRU1INWIb2DOgQwDW9w+gb6YfaGeOLhHARKpWKAE89AZ56ekb41jtW12WXVVxFToklacouriS7pIrskmpOlddQUWMiPb+c9Pzyc17fy6AlxMdA8OnEqZ23gXZeBoJ9DLTzcsPHXVqgmkujEqQOHTpYP7/hhhuYN28e48ePt+7r3bs3kZGRPPvss0yaNMnuQYpL1IJ+mL48+CWlxlI6+HQgsXuis8NpFRRF4UB2KT/vzmLD4XzrbDQvg5ZR3YK5pncYYb4uunivEC3ImV123cN9GhyvMprILakmt7SK7JIqckqqyS2xJFM5JdWUVddatrxajuSdO4HSaVSWpMnbkjAFeesJ8jIQ5GVJpIK89XjoZXixPdh8F3fv3k3Hjg0r0Xbs2JF9+/bZJSjRNqUVp7E2Yy0A9/S6RwZlX0RljYk1B3P5aU82R8/4bzM22IvxvcK4Mi5IllkQohVx02mICvQgKtDjnMcrampPJ1CWJCq3pJq8smrrx1PlNRhNirVcAefowgNL92A7LwOBXnoCPS0fg7z0BHoZCPS07JOWqIuzOUHq1q0bs2bN4oMPPkCvt5Rzr6mpYdasWXTr1s3uAYq2QVEUPt37KQoKV4RfQWf/zs4OqcU6klfG8j3ZrD2YR6XRMqZPp1ExrHMw43uFEhci0+SFaIs89Fqig7REB3me83hNrZnC8hrySqvJK6siv7SGvLJq8kqryS+zbOXVJiprTNYyBuej1agI9NTj76EnwEtPkKfB2nV45uah17hsImVzgrRw4UImTJhA+/btrTPWdu3ahUql4ocffrB7gKJt2JK9hf2F+9Gr9dza7VZnh9PiVNaYWJ+ax/I92aTmlln3R/i5M65XKCO7BkvROiFcnF6rts6KA99znlNZY7ImSwVlNZaE6vTnBeXV1vIFtSbFOhPvQgxaNf6eegI89Ph56gjw0Fsf+3vqLAmWpx4fN12bG/9oc4KUkJBAWloaX3zxBQcOHADgpptu4tZbb8XT89xZr3BtRpORz/d9DsCEmAmyAO0Z0vLKWL43mzUH/mwt0qhVXB4TyLieYfSM8HHZ/96EELZz12uIDLCsPXc+RpOZU+U1FJRbEqj80913heU1FFbUWBOrihoT1bVmsouryC6uuuDrqlXg425JmPw9dPie/ujvocf39Ec/dx1+HrpWk0w1aSSXp6cn9913n71jEW3UT+k/kVuZi7+bPxNiJjg7HKerrDGxLjWPX85qLQr1dWNsj1BGdQuW1ciFEA6j06gts+R8Llwjrcpo4lSFJVk6VW6ksKKGU+U1nLJ+NHKqoobiSiNmBWthzfQLXvXPZMr3dMLk5663JE7uOvys+/X4uuvw99Q5bZHhJiVIn332Ge+99x5paWls3LiRDh06MGfOHDp16sR1111n7xhFK1ZZW8nSw0sBuLXrrbhrXXe21eHcUn7Zm1NvbJFGrWJQp0DG9gyld4Rvq/ivSgjhGtx0GsJ83S86S9ZkViiutCRLRRWWZKoucSosr6Go0kjx6WSqtKq2XjJ1rODCMfxtWCeu7R1ux3fVeDYnSO+++y4zZ85k2rRpvPzyy9bCkP7+/sydO1cSJFHP+hPrqaytJNwz3CVrHlXU1LLukGVs0ZnTdsN83RgjrUVCiDZAo/6zNtTF1JrMlFTVWhKpCiMllUaKK40UVZxOpE4nU3UJl6+788Ze2pwgzZ8/n/fff59Jkybx2muvWfcPGDCAxx57zK7BidZNURR+OfoLAFdHX41a5RrrEymKwqEcy0y09al5VNeaAcuskctjAhnTI5ReEb4ytkgI4XK0GnWjkylFUVCUZgjqPGxOkNLT0+nbt2+D/QaDgfLycxe2Eq5pf+F+TpSdwKAxMLT9UGeH43Bl1bWsOZjLL3tz6tUtigxwZ0yPUIZ3CXbqf0NCCNGaqFQqp9Y5tjlB6tixIykpKfWqawMsX75c6iCJen49+isAV0Zciaeubc5wVBSF/VmlLN+bze+H86k53Vqk06gYEteOMT1C6B4mM9GEEKK1sTlBSkpK4sEHH6SqqgpFUdiyZQtffvmltXikEACnqk6xJXsLYFlzra0prTKy6kAuv+7NqVeMrUOgB2N7hjKsczupWySEEK2YzQnSPffcg7u7O//4xz+oqKjg1ltvJTw8nLfffpubb77ZETGKVmjl8ZWYFBNd/LsQ7Rvt7HDsQlEU9mWV8Mue7Hprohm0aq6Ma8fYnqF0DvGS1iIhhGgDmjTN/7bbbuO2226joqKCsrIygoOD7R2XaMVqzbX8dvw3wDI4u7Urq65l9YFclu/JrtdaFB3kydgeoQzv0g5PgywOKYQQbYnNv9VffPFFhgwZwsiRI/Hw8MDDw1Kts7y8nLfeeouZM2faPUjRuiTnJHOq6hQ+eh8Ghg50djhNoigKh3PL+Gl3NutS86xji/RaNUOltUgIIdo8mxOk559/Hp1Ox6xZs0hKSrLuLysr44UXXpAESbDy+EoARkWNQqdpXeNwqowm1h3K46fdWfXqFkUFejCuZygjugRLa5EQQriAJv2m//TTT3nwwQfZvXs37733Hnq9FLoTFpW1lezN3wvAle2vdHI0jXfiVAU/787mt/05VNRYip/qNCqGxAYxrlcYXUO9pbVICCFcSJMSpBEjRrB582YmTJjA8OHD+fbbb+0clmitduXtolapJdQjlHBP55SHbyyzWWHL0UJ+3JVFSkaRdX+Ijxvje4UyqluI1C0SQggXZXOCVPdfdExMDJs2beLGG2+kf//+LFy40O7BidZne+52APqF9GuxLS4lVUZ+3ZvDT7uzyCutBiyLJw6IDmB8rzD6RvrJmmhCCOHibE6QlDPqfvv4+PDTTz8xbdo0Jk2aZM+4RCtkVszsyNkBQL/gfk6OpqH0/HKW7TzJ6oO51in63m5aru4ewrheYYRcZGVr0Ya10GReCOE8NidIH3/8Mb6+vtbHarWaefPm0bdvX9atW2fX4ETrkl6cTnFNMW4aN7oGdnV2OIClG23r0UK+23mS3SeKrfs7tfNkQu9whnZuh17rGmvECSGEaDybE6TExMRz7r/zzju58847Lzkg0Xol5yQD0KddH3Rq547dqTKa+G1/Dt+nnCSruAqwdKMNjgliYp9wuoXJoGshhBDn16gEad68edx33324ubkxb968856nUql46KGH7BacaF125J7uXgtxXvdaYXkNP+46yU+7symrrgXA06BhTI9Qru0dTjtvg9NiE0II0Xo0KkGaM2cOt912G25ubsyZM+e857XGBKmiooJu3bpxww038M9//tPZ4bRahVWFpBWnAdA3uG+zv35GYQXfbM9kzaFcak+PLwrzdWNifDiju4XgptM0e0xCCCFar0YlSOnp6ef8vC145ZVXGDRokLPDaPXqWo9i/WLxNfhe5Gz7OZhdyv+2n2BTWgF18we6hnpzfb8IBnUMlNloQgghmsSlSwKnpqZy4MABJkyYwJ49e5wdTqvWnLPXFEVhd2Yxi7dmsOuMgdcJHQP4v/7t6Rbm4/AYhBBCtG2NSpDOXFLkYmbPnt3kYM60bt063nzzTZKTk8nKymLp0qUNSgksWLCAN998k+zsbPr06cP8+fNJSEho9Gs89thjvPnmm/zxxx92idlVGU1GdufvBhw7/khRFLYfP8VXWzI4kF0KgFqtYnjndvylX3uiAj0c9tpCCCFcS6MSpB07djTqYvacFVReXk6fPn246667mDx5coPjixcvJikpiYULFzJw4EDmzp3LmDFjOHjwIMHBwQDEx8dTW1vb4Lm//vorW7dupXPnznTu3LlRCVJ1dTXV1dXWxyUlJZfw7tqWvQV7qTJV4e/mT7RPtN2vX5cYfbH5OKk5ZYBlGZCre4QyuW8EwVK/SAghhJ01KkFavXq1o+NoYNy4cYwbN+68x2fPns29995rLS2wcOFCfvzxRz766COeeuopAFJSUs77/E2bNvHVV1+xZMkSysrKMBqN+Pj4nHex3VmzZvHCCy80/Q05hXLxU+xgT76lezK+Xbxdk2RFUUjJKOKLzcc5eLrFSK9VM65nKJP7tSfAU9YAFEKcn9JMvwNF29QqxyDV1NSQnJzMjBkzrPvUajWjR49m48aNjbrGrFmzmDVrFgCLFi1iz549502OAGbMmFGvq7GkpITIyMgmvoO25XDRYQC6+Hex2zUPZpfyycaj1uKOOo2K8b3C+L/+7fHzkMRICCGEYzUpQdq2bRtff/01x48fp6ampt6xb775xi6BXUh+fj4mk4mQkJB6+0NCQjhw4IBDXtNgMGAwSA2ds5kVM+nFlpmNsX6xl3y94wUVfLbpKJvSCgHQalSM6xnKDf0j8ZcWIyGEEM3E5gTpq6++4o477mDMmDH8+uuvXH311Rw6dIicnByuv/56R8TocFOmTHF2CK1WVnkWVaYqDBoD4V7hTb5OUUUNX2w+zq97szErlqrXI7uGcEtCpIwxEkII0exsTpBeffVV5syZw4MPPoi3tzdvv/02HTt25G9/+xthYWGOiLGBoKAgNBoNOTk59fbn5OQQGhraLDEIiyNFRwCI9olGo7a9GGN1rYnvdpzkv8knqDSaABjUKYA7BkcTGSCz0oQQQjiHzat0HjlyhGuuuQYAvV5PeXk5KpWK6dOn8+9//9vuAZ6LXq+nf//+rFy50rrPbDazcuVKBg8e3CwxtC6OK5aYVmSpnh3jF2PT8xRFYeORAh78YjufbTpGpdFEXLAXsyb34plruktyJIQQwqlsbkHy9/entNQyoygiIoI9e/bQq1cvioqKqKiosFtgZWVlHD582Po4PT2dlJQUAgICiIqKIikpicTERAYMGEBCQgJz586lvLxcFsxtZkeKLS1InXw7Nfo5GYUVvL8+jR3HiwAI9NIz5fJohsa1k8rXQgghWgSbE6ShQ4eyYsUKevXqxQ033MAjjzzCqlWrWLFiBaNGjbJbYNu2bWPEiBHWx3UzyBITE1m0aBE33XQTeXl5zJw5k+zsbOLj41m+fHmDgdvCcUxmE0eLjwKNa0GqqTXz9bYM/pt8ApNZQatRMblvBDcMiJS10oQQQrQoNidI77zzDlVVVQA888wz6HQ6/vjjD/7yl7/wj3/8w26BDR8+HEW5cA2LqVOnMnXqVLu9prDNibIT1JhrcNO4Eep54bFfu08U887qVE4WWb53BkT7c++VnQj3c2+OUIUQQgib2JwgBQQEWD9Xq9XWoozC9dSNP+rk1wm16tzD2cqra/lwQzor9lkG1Pt76rl/aCcGxwTataikEEIIYU9NLhSZm5tLbm4uZrO53v7evXtfclCidagbf3S++kcpGUW8/dsh8ssstbLG9QolcXA0noZWWZ9UCCGEC7H5L1VycjKJiYns37+/QReYSqXCZDLZLTjRstW1IHX07Vhvf5XRxKI/jvLjriwAQnzcmDY6jp4Rvs0eoxBCCNEUNidId911F507d+bDDz8kJCREuklclNFs5FjpMQBifP8coJ2WV8Ybyw+SWVQJwPheYUy5PBp3vQzCFkII0XrYnCClpaXxv//9j9jYS19WQrReGSUZ1Jpr8dJ5EewRjKIo/Lwnmw/Wp2E0KQR66XlkVBx9o/ydHaoQQghhM5sTpFGjRrFz505JkFxcWvHpAdq+naioMTFvVSp/HC4ALDPUpo3ujK+7zpkhCiGEEE1mc4L0wQcfkJiYyJ49e+jZsyc6Xf0/ghMnTrRbcKLlqltiJNAQxaNf7ySzqBKNWsWUy6O5Lj5cul6FEEK0ajYnSBs3buT333/n559/bnBMBmm7jiPFR6ioNrFsWy2qqkoCvfQ8Pb4bnUO8nR2aEEIIcclsXovtoYce4q9//StZWVmYzeZ6myRHrqG6tppd2Uc4WVyJuTqErqHezLkxXpIjIYQQbYbNLUgFBQVMnz5dlvRwUWazwj9XbyCvrAq14smYrrE8MDwWvdbmXFsIIYRosWz+qzZ58mRWr17tiFhEC1dTa+aNXw6y6vAhAHqHxPDwqDhJjoQQQrQ5Nrcgde7cmRkzZrBhwwZ69erVYJD2ww8/bLfgRMtRUVPLKz/uZ9eJYsz6IkJ93BgU1UkGYwshhGiTmjSLzcvLi7Vr17J27dp6x1QqlSRIbVBZdS3PfruHw7lluOs09OmqI7VUS4indLOKNuYiC2QLIVyHTQmSoiisWbOG4OBg3N1lFXZXUFlj4rnv9nI4twwfdy3PT+jBJ6k/AhDsEezk6IQQQgjHsGnwiKIoxMXFceLECUfFI1qQKqOJF5ft5VBOKV4GLS9P6kVciDc5FTkAhHhIC5IQQoi2yaYESa1WExcXR0FBgaPiES1ETa2Zl3/cx57MEtz1Gl6a1IOOQZ5U1lZSUlMCSIIkhBCi7bJ5+tFrr73G448/zp49exwRj2gBzGaFN385wM6MYtx0al6Y2IPYYEuNo5xyS+uRt94bD52HM8MUQgghHMbmQdp33HEHFRUV9OnTB71e32AsUmFhod2CE87x+eZjbEorRKdRMfPaHnQL87Eek+410TbJbEwhRH02J0hz5851QBiipVh9MJcl2yxjzB4aFUev9r71jkuCJIQQwhXYnCAlJiY6Ig7RAhzMLmX+ylQA/q9/e0Z0aThLra6LTRIkIYQQbZnNCRKAyWTi22+/Zf/+/QD06NGDiRMnotFo7BqcaD75ZdW8/OM+jCaFhI4B3D6owznPs7YgSQ0kIYQQbZjNCdLhw4cZP348mZmZdOnSBYBZs2YRGRnJjz/+SExMjN2DFI5lNiu89etBiiqMRAV68NjVXVCrzz0mQ1qQhBBCuAKbZ7E9/PDDxMTEkJGRwfbt29m+fTvHjx+nY8eOUkW7lfph10n2ZJbgplPzzPhuuOvP3RJYa64lvzIfkBYkIYQQbZvNLUhr165l06ZNBAQEWPcFBgby2muvccUVV9g1OOF4xwsq+OSPowDcPaQj4X7nr5CeX5mPGTM6tQ5/g38zRSiEEEI0P5tbkAwGA6WlpQ32l5WVodfr7RKUsJOLrCtVazIze8VBjCaF/h38GdMj9ILnnzmDTRapFUK0dAqytp5oOpsTpGuvvZb77ruPzZs3oygKiqKwadMm7r//fiZOnOiIGIWDLN6WwZG8crwMWh4aGXvRpMc6/ki614QQQrRxNidI8+bNIyYmhsGDB+Pm5oabmxtXXHEFsbGxvP32246IUThAWl4ZX2/NAOCB4TEEehku+hypgSSEEMJV2DwGyc/Pj++++47U1FQOHDgAQLdu3YiNjbV7cMJxPv79KGYFLo8JZGjndo16jrQgCSGEcBVNqoMEEBcXR1xcnD1jEc1kx/FTpGQUoVGruGtIx0Y/r64FKdTjwmOVhBBCiNbO5gTJZDKxaNEiVq5cSW5uLmazud7xVatW2S04YX9ms8Ki07PWrukVRoiPW6OepygKuRW5gHSxCSGEaPtsTpAeeeQRFi1axDXXXEPPnj1lNlMrszY1j7S8ctz1Gm68LLLRzyupKaHKVIUKFUHuQQ6MUAghhHA+mxOkr776iq+//prx48c7Ih7hQDW1Zr7YdAywrLXm665r9HOzy7MBCHQPRKdp/POEEEKI1sjmWWx6vV4GZLdSP+/JIqekmkAvPRP7hNv0XJnBJoQQwpXYnCA9+uijvP322ygXKUIoWpYqo4mvtlim9d+aEIWbzraFhSVBEkII4Ups7mLbsGEDq1ev5ueff6ZHjx7odPW7W7755hu7BSfsZ31qPmXVtYT4uDGqm+1JjkzxF0II4UqaVAfp+uuvd0QswoF+3WsZQzSmRwgate0D66UFSQghhCuxOUH6+OOPHRGHcKCs4koOZJtQq1WMbkLrEZzRgiQJkhBCCBdg8xgk0fr8kVYAwMCOAfh72r6gcI2phuKaYgCCPYLtGpsQQtibCik/Iy6dJEhtnBnYln4KsHSvNUVxtSU50ql1eOo87RWaEEII0WJJgtTGlVXXUmk00c7bQN9I/yZdoy5B8tH7SGFQIYQQLkESpDaupNIIwFXdQ1A3YXA2YO1e8zP42SssIYQQokWzS4JUVFRkj8sIO6swmqisMaFSWRKkpiqpLgHAx+Bjr9CEEEKIFs3mBOn1119n8eLF1sc33ngjgYGBREREsHPnTrsGJy5NTkk1AN3CfAnyMjT5OnUtSL56X7vEJYQQQrR0NidICxcuJDLSssjpihUrWLFiBT///DPjxo3j8ccft3uAoumKK2sAGBDtd2nXqRuDJC1IQgghXITNdZCys7OtCdKyZcu48cYbufrqq4mOjmbgwIF2D1A0Xd1qMJ56m7/M9Zw5SFsIIYRwBTa3IPn7+5ORYVnTa/ny5YwePRoARVEwmUz2jU60CCU1ljFIvgbpYhNCCOEabG5amDx5MrfeeitxcXEUFBQwbtw4AHbs2EFsbKzdAxTOVzdIWxIkIYQQrsLmBGnOnDlER0eTkZHBG2+8gZeXFwBZWVn8/e9/t3uAwvlkkLYQQghXY3OCpNPpeOyxxxrsnz59ul0CEi2LWTHLNH8hhBAup0l1kD777DOGDBlCeHg4x44dA2Du3Ll89913dg1OOF+ZsQwzZgC89d5OjkYIIYRoHjYnSO+++y5JSUmMGzeOoqIi68BsPz8/5s6da+/4hJPVzWDz0nmhU+ucHI0QQgjRPGxOkObPn8/777/PM888g0ajse4fMGAAu3fvtmtwwvlkBpsQQghXZHOClJ6eTt++fRvsNxgMlJeX2yUo0XJIDSQhhBCuyOYEqWPHjqSkpDTYv3z5crp162aPmJpNeno6I0aMoHv37vTq1UsSvHOoa0GSAdqiTVM1bSFnIUTbZfMstqSkJB588EGqqqpQFIUtW7bw5ZdfMmvWLD744ANHxOgwU6ZM4eWXX+bKK6+ksLAQg6Hp65W1VcVVMsVfCCGE67E5Qbrnnntwd3fnH//4BxUVFdx6662Eh4fz9ttvc/PNNzsiRofYu3cvOp2OK6+8EoCAgAAnR9QySQuSEEIIV9Skaf633XYbqamplJWVkZ2dzYkTJ7j77rvtGti6deuYMGEC4eHhqFQqvv322wbnLFiwgOjoaNzc3Bg4cCBbtmxp9PVTU1Px8vJiwoQJ9OvXj1dffdWO0bcddWOQ/Ax+zg1ECCGEaEY2J0gjR46kqKgIAA8PD4KDgwEoKSlh5MiRdgusvLycPn36sGDBgnMeX7x4MUlJSTz33HNs376dPn36MGbMGHJzc63nxMfH07NnzwbbyZMnqa2tZf369fzrX/9i48aNrFixghUrVpw3nurqakpKSuptrqCuirYM0hZCCOFKbO5iW7NmDTU1NQ32V1VVsX79ersEBTBu3DjrOm/nMnv2bO69917uvPNOABYuXMiPP/7IRx99xFNPPQVwzsHkdSIiIhgwYACRkZEAjB8/npSUFK666qpznj9r1ixeeOGFJr4b51ChXPI1ZB02IUSrdem/AoULa3SCtGvXLuvn+/btIzs72/rYZDKxfPlyIiIi7BvdedTU1JCcnMyMGTOs+9RqNaNHj2bjxo2NusZll11Gbm4up06dwtfXl3Xr1vG3v/3tvOfPmDGDpKQk6+OSkhJrctWWFVUXATJIWwghhGtpdIIUHx+PSqVCpVKdsyvN3d2d+fPn2zW488nPz8dkMhESElJvf0hICAcOHGjUNbRaLa+++ipDhw5FURSuvvpqrr322vOebzAYXG6WW42phipTFSCDtIUQQriWRidI6enpKIpCp06d2LJlC+3atbMe0+v1BAcH16us3RpcrBvP1dXNYNOqtXhoPZwcjRBCCNF8Gp0gdejQAQCz2eywYBorKCgIjUZDTk5Ovf05OTmEhoY6KaqWrGlF8M6soq2SQnpCCCFciM2DtOvs27eP48ePNxiwPXHixEsO6mL0ej39+/dn5cqVTJo0CbAkbitXrmTq1KkOf31XIeuwCSGEcFU2J0hpaWlcf/317N69G5VKhaJYpgnUtTCYTCa7BFZWVsbhw4etj9PT00lJSSEgIICoqCiSkpJITExkwIABJCQkMHfuXMrLy62z2sSlK6oqAmSAthBCCNdjcx2kRx55hI4dO5Kbm4uHhwd79+5l3bp1DBgwgDVr1tgtsG3bttG3b1/rwrhJSUn07duXmTNnAnDTTTfxz3/+k5kzZxIfH09KSgrLly9vMHBbNJ1U0RZCCOGqbG5B2rhxI6tWrSIoKAi1Wo1arWbIkCHMmjWLhx9+mB07dtglsOHDh1tbp85n6tSp0qXmQFJFWwghhKuyuQXJZDLh7e0NWAZLnzx5ErAM4j548KB9oxNOZW1BkiraQgghXIzNLUg9e/Zk586ddOzYkYEDB/LGG2+g1+v597//TadOnRwRo3CSuhYkGaQthBDC1dicIP3jH/+gvLwcgBdffJFrr72WK6+8ksDAQL766iu7ByicR9ZhE0II4apsTpDGjBlj/Tw2NpYDBw5QWFiIv7+/1MppY+rWYZNB2kIIIVyNzWOQziUgIIDs7GwZMN2GmBXzn3WQZJq/EEIIF2NTC9LevXtZvXo1er2eG2+8ET8/P/Lz83n55Zd57733ZAxSG1JhrMCkWGpaSQuSEEIIV9PoFqTvv/+evn378vDDD3P//fczYMAAVq9eTbdu3Thw4ABLly5l7969joxV2OjCRRIurG78kafOE51aZ5+AhBBCiFai0QnSyy+/zIMPPkhJSQmzZ88mLS2Nhx9+mJ9++only5czduxYR8YpmtmZ67AJIYQQrqbRCdLBgwd58MEH8fLy4qGHHkKtVjNnzhwuu+wyR8YnnESm+AshhHBljU6QSktL8fGxtCZoNBrc3d1lzFEbVtfFJgO0hRCtjQqZUS0unU2DtH/55Rd8fS1/MM1mMytXrmTPnj31zpk4caL9ohNOI1P8hRBCuDKbEqTExMR6j//2t7/Ve6xSqTCZTJcelXA66xR/6WITQgjhghqdIJnNZkfGIVoYGaQthBDCldmlUKRoeyRBEkII4cokQXIBTRmuWG60rLfnpfOybzBCCCFEKyAJkjinKlMVAO5adydHIoQQQjQ/SZDEOVXVnk6QdJIgCSGEcD2SIIkGFEWhwlgBgJvGzcnRCNGcLmWBHiFEW2LTNP8z1dTUkJub22B2W1RU1CUHJZzLaDZixvJ1lS42IYQQrsjmBCk1NZW77rqLP/74o95+RVGkDlIbUVlbaf3cTSstSEIIIVyPzQnSlClT0Gq1LFu2jLCwMFQqKene1tQlSG4aN9Qq6YUVQgjhemxOkFJSUkhOTqZr166OiEe0AHUJknSvCdch/+gJIeqzuXmge/fu5OfnOyIW0UJYW5Cke00IIYSLsjlBev3113niiSdYs2YNBQUFlJSU1NtE61c3xd9D6+HkSIQQQgjnsLmLbfTo0QCMGjWq3n4ZpN12SAuSEEIIV2dzgrR69WpHxCEcoKmjKiRBEkII4epsTpCGDRvmiDhEC1KXIEkXmxBCCFfVpEKRRUVFfPjhh+zfvx+AHj16cNddd+Hr62vX4IRz1K3DJi1IQgghXJXNg7S3bdtGTEwMc+bMobCwkMLCQmbPnk1MTAzbt293RIziEik21qqqW2ZEpvkLIYRwVTa3IE2fPp2JEyfy/vvvo9Vanl5bW8s999zDtGnTWLdund2DFM3LulCtJEhCCCFclM0J0rZt2+olRwBarZYnnniCAQMG2DU44Rx1XWySIAkhhHBVNnex+fj4cPz48Qb7MzIy8Pb2tktQwrnqutjcNDIGSQghhGuyOUG66aabuPvuu1m8eDEZGRlkZGTw1Vdfcc8993DLLbc4IkbRzKwtSDppQRJCCOGabO5i++c//4lKpeKOO+6gtrYWAJ1OxwMPPMBrr71m9wBF85O12IQQQrg6mxMkvV7P22+/zaxZszhy5AgAMTExeHhIzZy2otJ4ulCkdLEJIYRwUU2qgwTg4eFBr1697BmLaCEqTVIoUgghhGtrVII0efJkFi1ahI+PD5MnT77gud98841dAhN2oDTtaXXT/KVQpBBCCFfVqATJ19cX1eligz4+PtbPRdtjNBsxmo2AjEESQgjhuhqVIH388cfWzxctWuSoWEQLUNd6BJIgCSGEcF02T/MfOXIkRUVFDfaXlJQwcuRIe8QknKhuBptOrUOj1jg5GiGEEMI5bE6Q1qxZQ01NTYP9VVVVrF+/3i5BCeepa0GSAdpCCCFcWaNnse3atcv6+b59+8jOzrY+NplMLF++nIiICPtGJ5pdXQuSDNAWQgjhyhqdIMXHx6NSqVCpVOfsSnN3d2f+/Pl2DU40P0mQhBBCCBsSpPT0dBRFoVOnTmzZsoV27dpZj+n1eoKDg9FoZMxKa1eXIEkXmxBCCFfW6ASpQ4cOAJjNZocFI5xPlhkRQgghGpkgff/994wbNw6dTsf3339/wXMnTpxol8CEc0gXmxBCCNHIBGnSpElkZ2cTHBzMpEmTznueSqXCZDLZKzbhBHWz2KQFSQghhCtrVIJ0ZreadLG1bVUmSZCEEEIIm+sgncu5CkeK1qnCWAFIF5sQovWS5bCEPdicIL3++ussXrzY+viGG24gICCAiIgIdu7cadfgRPOTQdpCCCFEExKkhQsXEhkZCcCKFSv47bffWL58OePGjePxxx+3e4CieUkXmxBCCGHDNP862dnZ1gRp2bJl3HjjjVx99dVER0czcOBAuwcomlddF5u7RhIkIYQQrsvmFiR/f38yMjIAWL58OaNHjwZAURSZwdZC2dIdb21B0kmCJIQQwnXZnCBNnjyZW2+9lauuuoqCggLGjRsHwI4dO4iNjbV7gI40Z84cevToQffu3Xn44YdRFMXZITmdTPMXLk1+BwghTrO5i23OnDlER0eTkZHBG2+8gZeXFwBZWVn8/e9/t3uAjpKXl8c777zD3r170el0DB06lE2bNjF48GBnh+ZU1llsGpnFJoQQwnXZnCDpdDoee+yxBvunT59ul4CaU21tLVVVlhYTo9FIcHCwkyNyPhmkLYQQQjQhQQJITU1l9erV5ObmNigcOXPmTLsEtm7dOt58802Sk5PJyspi6dKlDap4L1iwgDfffJPs7Gz69OnD/PnzSUhIaNT127Vrx2OPPUZUVBRarZb777+fmJgYu8TeWpkVM9WmakDqIAkhhHBtNidI77//Pg888ABBQUGEhobWK8ilUqnsliCVl5fTp08f7rrrLiZPntzg+OLFi0lKSmLhwoUMHDiQuXPnMmbMGA4ePGhtCYqPj6e2trbBc3/99Vfc3d1ZtmwZR48exd3dnXHjxrFu3TqGDh1ql/hbo7rxRwAeWg8nRiKEEEI4l80J0ssvv8wrr7zCk08+6Yh4rMaNG2cdAH4us2fP5t577+XOO+8ELPWZfvzxRz766COeeuopAFJSUs77/CVLlhAbG0tAQAAA11xzDZs2bTpvglRdXU11dbX1cUlJia1vqcWrqLWMP9KqtOg0OidHI4QQl0ZBBt2LprN5FtupU6e44YYbHBFLo9XU1JCcnGwtMQCgVqsZPXo0GzdubNQ1IiMj+eOPP6iqqsJkMrFmzRq6dOly3vNnzZqFr6+vdaurBdWW1FXRlu41IYQQrs7mBOmGG27g119/dUQsjZafn4/JZCIkJKTe/pCQELKzsxt1jUGDBjF+/Hj69u1L7969iYmJYeLEiec9f8aMGRQXF1u3ulpQbUldF5t0rwkhhHB1NnexxcbG8uyzz7Jp0yZ69eqFTle/K+bhhx+2W3CO9sorr/DKK6806lyDwYDBYHBwRI7SuEqR0oIkhBBCWNicIP373//Gy8uLtWvXsnbt2nrHVCpVsyRIQUFBaDQacnJy6u3PyckhNDTU4a/fVkmCJIQQQljYnCClp6c7Ig6b6PV6+vfvz8qVK61T/81mMytXrmTq1KnODa4VkyraQgghhEWT6iA1h7KyMg4fPmx9nJ6eTkpKCgEBAURFRZGUlERiYiIDBgwgISGBuXPnUl5ebp3VJmxXabK0IEmCJIQQwtU1KUE6ceIE33//PcePH6empqbesdmzZ9slsG3btjFixAjr46SkJAASExNZtGgRN910E3l5ecycOZPs7Gzi4+NZvnx5g4HbovEqjdLFJoQQQkATEqSVK1cyceJEOnXqxIEDB+jZsydHjx5FURT69etnt8CGDx9+0cVjp06dKl1qdlQ3BklmsQkhhHB1Nk/znzFjBo899hi7d+/Gzc2N//3vf2RkZDBs2DCn10cSl0bWYRNCCCEsbE6Q9u/fzx133AGAVqulsrISLy8vXnzxRV5//XW7ByiaT4XRUklbutiEEEK4OpsTJE9PT+u4o7CwMI4cOWI9lp+fb7/IRLOTWWxCCCGEhc1jkAYNGsSGDRvo1q0b48eP59FHH2X37t188803DBo0yBEximYiXWxCCCGEhc0J0uzZsykrKwPghRdeoKysjMWLFxMXF2e3GWzCOaxdbBrpYhNCCOHabEqQTCYTJ06coHfv3oClu23hwoUOCUw0P2sLkk5akIQQQrg2m8YgaTQarr76ak6dOuWoeIQT1U3zd9dIgiSEEMK12TxIu2fPnqSlpTkiFuFk1gRJxiAJIYRwcTYnSC+//DKPPfYYy5YtIysri5KSknqbaJ0URZFZbEIIIcRpjR6D9OKLL/Loo48yfvx4ACZOnIhKpbIeVxQFlUqFyWSyf5TC4apMVShYKpdLHSQhhBCurtEJ0gsvvMD999/P6tWrHRmPcJK61iM1agwag5OjEUIIIZyr0QlS3bpow4YNc1gwwnnqxh+5ad3qtQwKIYQQrsimMUjyh7PtOjNBEkIIIVydTXWQOnfufNEkqbCw8JICEs4hA7SFEEKIP9mUIL3wwgv4+vo6KhbhRBW1lirakiAJIYQQNiZIN998M8HBwY6KRTiR1EASQggh/tToMUgy/qhtky42IYQQ4k+NTpDqZrGJtsm6DpskSEIIIUTju9jMZrMj4xBOVmG0jEGSWWxCCCFEE5YaEW1TtakaQIpECiGEEEiCJE4zmo0A6DV6J0cihBCXRoWMmRWXThIkAfyZIGlVNk1sFEIIIdokSZBcQGP+mzKaLAmSTqNzdDhCtGAyGUUIYSEJkgCg1lwLgE4tCZIQQgghCZIA/uxikwRJuCSp8yaEOIskSG1a47sLrGOQ1DIGSQjRNkj9PnEpJEESwJ9dbDKLTQghhJAEyUU0YpC2zGITQgghrCRBEsAZY5BkFpsQQgghCVLbZsMYJJMM0hZCCCHqSIIkAJnmL4QQQpxJEiQByDR/IYQQ4kySIAlApvkLIYQQZ5IESQCyWK0QQghxJkmQBPDnIG1pQRJCCCEkQRKAyWzCjBmQMUhCCCEESIIk+LN7DSRBEkIIIUASJMGfU/xBEiQhhBACJEES/NmCpEaNRq1xcjRCCCGE80mCJGSZESGEEOIskiAJKRIphBBCnEUSJCFT/IUQQoizSIIkrIO09WopEimEEEKAJEgCqDHXANKCJIQQQtSRBElYW5BkDJIQQghhIQmSkIVqhRBCiLNIgiRkoVohhBDiLJIgCZnFJoQQQpxFEiQhY5CEEEKIs0iCJKyz2CRBEkIIISwkQRLWFiTpYhNCCCEsJEESshabEEIIcRZJkISsxSaEEEKcxSUSpOuvvx5/f3/+7//+r8GxZcuW0aVLF+Li4vjggw+cEJ3z1c1ikwRJCCGEsHCJBOmRRx7h008/bbC/traWpKQkVq1axY4dO3jzzTcpKChwQoTOJYUihRBCiPpcIkEaPnw43t7eDfZv2bKFHj16EBERgZeXF+PGjePXX391QoTOZV2sVgpFCiGEEEALSJDWrVvHhAkTCA8PR6VS8e233zY4Z8GCBURHR+Pm5sbAgQPZsmWLXV775MmTREREWB9HRESQmZlpl2u3JtbFalXSgiSEEEJAC0iQysvL6dOnDwsWLDjn8cWLF5OUlMRzzz3H9u3b6dOnD2PGjCE3N9d6Tnx8PD179mywnTx5srneRqtmLRQps9iEEEIIAJzeZDBu3DjGjRt33uOzZ8/m3nvv5c477wRg4cKF/Pjjj3z00Uc89dRTAKSkpDTptcPDw+u1GGVmZpKQkHDOc6urq6murrY+Li4uBqCkpKRJr90cyqqMKFUmKCu/YJylJaUYK4xUl1e36PcjhMNUGKGyFkpKwSA/A61deWk5xgojlWWV8jtNNFD3PaEoyoVPVFoQQFm6dKn1cXV1taLRaOrtUxRFueOOO5SJEyfadO3Vq1crf/nLX+rtMxqNSmxsrHLixAmltLRU6dy5s5Kfn3/O5z/33HMKIJtssskmm2yytYEtIyPjgnmD01uQLiQ/Px+TyURISEi9/SEhIRw4cKDR1xk9ejQ7d+6kvLyc9u3bs2TJEgYPHoxWq+Wtt95ixIgRmM1mnnjiCQIDA895jRkzZpCUlGR9bDabKSwsJDAwEJVK1bQ3eA4lJSVERkaSkZGBj4+P3a7bVsn9ajy5V40n96rx5F41ntyrxnPkvVIUhdLSUsLDwy94XotOkOzlt99+O++xiRMnMnHixItew2AwYDAY6u3z8/O71NDOy8fHR36AbCD3q/HkXjWe3KvGk3vVeHKvGs9R98rX1/ei5zh9kPaFBAUFodFoyMnJqbc/JyeH0NBQJ0UlhBBCiLauRSdIer2e/v37s3LlSus+s9nMypUrGTx4sBMjE0IIIURb5vQutrKyMg4fPmx9nJ6eTkpKCgEBAURFRZGUlERiYiIDBgwgISGBuXPnUl5ebp3V1tYYDAaee+65Bt154tzkfjWe3KvGk3vVeHKvGk/uVeO1hHulOj17zGnWrFnDiBEjGuxPTExk0aJFALzzzju8+eabZGdnEx8fz7x58xg4cGAzRyqEEEIIV+H0BEkIIYQQoqVp0WOQhBBCCCGcQRIkIYQQQoizSIIkhBBCCHEWSZBamAULFhAdHY2bmxsDBw5ky5Ytzg7J6WbNmsVll12Gt7c3wcHBTJo0iYMHD9Y7p6qqigcffJDAwEC8vLz4y1/+0qB+lit67bXXUKlUTJs2zbpP7tWfMjMz+etf/0pgYCDu7u706tWLbdu2WY8risLMmTMJCwvD3d2d0aNHk5qa6sSIncNkMvHss8/SsWNH3N3diYmJ4aWXXqq3lpUr36t169YxYcIEwsPDUalUfPvtt/WON+beFBYWctttt+Hj44Ofnx933303ZWVlzfgumseF7pXRaOTJJ5+kV69eeHp6Eh4ezh133NFg4fnmuleSILUgixcvJikpieeee47t27fTp08fxowZQ25urrNDc6q1a9fy4IMPsmnTJlasWIHRaOTqq6+mvLzces706dP54YcfWLJkCWvXruXkyZNMnjzZiVE739atW3nvvffo3bt3vf1yryxOnTrFFVdcgU6n4+eff2bfvn289dZb+Pv7W8954403mDdvHgsXLmTz5s14enoyZswYqqqqnBh583v99dd59913eeedd9i/fz+vv/46b7zxBvPnz7ee48r3qry8nD59+rBgwYJzHm/MvbntttvYu3cvK1asYNmyZaxbt4777ruvud5Cs7nQvaqoqGD79u08++yzbN++nW+++YaDBw82WO2i2e6VTSu+CodKSEhQHnzwQetjk8mkhIeHK7NmzXJiVC1Pbm6uAihr165VFEVRioqKFJ1OpyxZssR6zv79+xVA2bhxo7PCdKrS0lIlLi5OWbFihTJs2DDlkUceURRF7tWZnnzySWXIkCHnPW42m5XQ0FDlzTfftO4rKipSDAaD8uWXXzZHiC3GNddco9x111319k2ePFm57bbbFEWRe3UmqL/oemPuzb59+xRA2bp1q/Wcn3/+WVGpVEpmZmazxd7czr5X57JlyxYFUI4dO6YoSvPeK2lBaiFqampITk5m9OjR1n1qtZrRo0ezceNGJ0bW8hQXFwMQEBAAQHJyMkajsd6969q1K1FRUS577x588EGuueaaevcE5F6d6fvvv2fAgAHccMMNBAcH07dvX95//33r8fT0dLKzs+vdK19fXwYOHOhy9+ryyy9n5cqVHDp0CICdO3eyYcMGxo0bB8i9upDG3JuNGzfi5+fHgAEDrOeMHj0atVrN5s2bmz3mlqS4uBiVSmVd+7Q575XTK2kLi/z8fEwmEyEhIfX2h4SEcODAASdF1fKYzWamTZvGFVdcQc+ePQHIzs5Gr9c3WDw4JCSE7OxsJ0TpXF999RXbt29n69atDY7JvfpTWloa7777LklJSTz99NNs3bqVhx9+GL1eT2JiovV+nOtn0tXu1VNPPUVJSQldu3ZFo9FgMpl45ZVXuO222wDkXl1AY+5NdnY2wcHB9Y5rtVoCAgJc+v5VVVXx5JNPcsstt1gXrG3OeyUJkmhVHnzwQfbs2cOGDRucHUqLlJGRwSOPPMKKFStwc3NzdjgtmtlsZsCAAbz66qsA9O3blz179rBw4UISExOdHF3L8vXXX/PFF1/wn//8hx49epCSksK0adMIDw+XeyUcwmg0cuONN6IoCu+++65TYpAuthYiKCgIjUbTYDZRTk4OoaGhToqqZZk6dSrLli1j9erVtG/f3ro/NDSUmpoaioqK6p3vivcuOTmZ3Nxc+vXrh1arRavVsnbtWubNm4dWqyUkJETu1WlhYWF079693r5u3bpx/PhxAOv9kJ9JePzxx3nqqae4+eab6dWrF7fffjvTp09n1qxZgNyrC2nMvQkNDW0wGae2tpbCwkKXvH91ydGxY8dYsWKFtfUImvdeSYLUQuj1evr378/KlSut+8xmMytXrmTw4MFOjMz5FEVh6tSpLF26lFWrVtGxY8d6x/v3749Op6t37w4ePMjx48dd7t6NGjWK3bt3k5KSYt0GDBjAbbfdZv1c7pXFFVdc0aBcxKFDh+jQoQMAHTt2JDQ0tN69KikpYfPmzS53ryoqKlCr6/+50Gg0mM1mQO7VhTTm3gwePJiioiKSk5Ot56xatQqz2exy647WJUepqan89ttvBAYG1jverPfKrkO+xSX56quvFIPBoCxatEjZt2+fct999yl+fn5Kdna2s0NzqgceeEDx9fVV1qxZo2RlZVm3iooK6zn333+/EhUVpaxatUrZtm2bMnjwYGXw4MFOjLrlOHMWm6LIvaqzZcsWRavVKq+88oqSmpqqfPHFF4qHh4fy+eefW8957bXXFD8/P+W7775Tdu3apVx33XVKx44dlcrKSidG3vwSExOViIgIZdmyZUp6erryzTffKEFBQcoTTzxhPceV71VpaamyY8cOZceOHQqgzJ49W9mxY4d15lVj7s3YsWOVvn37Kps3b1Y2bNigxMXFKbfccouz3pLDXOhe1dTUKBMnTlTat2+vpKSk1Pt9X11dbb1Gc90rSZBamPnz5ytRUVGKXq9XEhISlE2bNjk7JKcDzrl9/PHH1nMqKyuVv//974q/v7/i4eGhXH/99UpWVpbzgm5Bzk6Q5F796YcfflB69uypGAwGpWvXrsq///3vesfNZrPy7LPPKiEhIYrBYFBGjRqlHDx40EnROk9JSYnyyCOPKFFRUYqbm5vSqVMn5Zlnnqn3R8uV79Xq1avP+TsqMTFRUZTG3ZuCggLllltuUby8vBQfHx/lzjvvVEpLS53wbhzrQvcqPT39vL/vV69ebb1Gc90rlaKcUQpVCCGEEELIGCQhhBBCiLNJgiSEEEIIcRZJkIQQQgghziIJkhBCCCHEWSRBEkIIIYQ4iyRIQgghhBBnkQRJCCGEEOIskiAJIYQQQpxFEiQhhBBCiLNIgiSEEEIIcRZJkIRoRYYPH860adNa7PUcdU3heM74uhUUFBAcHMzRo0cb/Zybb76Zt956y3FBCXGaJEhCOMGUKVNQqVSoVCr0ej2xsbG8+OKL1NbWXvB533zzDS+99JLd4rD39RorIyODu+66i/DwcPR6PR06dOCRRx6hoKCg2WMB5yd1dd8Pr732Wr393377LSqVyklROd4rr7zCddddR3R0dKOf849//INXXnmF4uJixwUmBJIgCeE0Y8eOJSsri9TUVB599FGef/553nzzzXOeW1NTA0BAQADe3t52i8He12uMtLQ0BgwYQGpqKl9++SWHDx9m4cKFrFy5ksGDB1NYWNis8bQUbm5uvP7665w6dcrZodhN3fftuVRUVPDhhx9y991323TNnj17EhMTw+eff36p4QlxQZIgCeEkBoOB0NBQOnTowAMPPMDo0aP5/vvvAUuLxtSpU5k2bRpBQUGMGTPGur+upWP48OE8/PDDPPHEEwQEBBAaGsrzzz9f7zXMZjNvvPEGsbGxGAwGoqKieOWVV6zHz245qXvdqVOn4uvrS1BQEM8++yyKoljPWb58OUOGDMHPz4/AwECuvfZajhw50uj3/eCDD6LX6/n1118ZNmwYUVFRjBs3jt9++43MzEyeeeYZ67nR0dHMnTu33vPj4+Prvc+LxXOx+zRlyhTWrl3L22+/bW3Vq+vyaczrDx8+nIceeohp06bh7+9PSEgI77//PuXl5dx55514e3sTGxvLzz//fMH7Mnr0aEJDQ5k1a9YFz3NkTLW1tRf82pvNZmbNmkXHjh1xd3enT58+/Pe//633uuf6vj2Xn376CYPBwKBBg6z7jhw5gkqlYtmyZYwaNQoPDw+6dOnC5s2b6z13woQJfPXVVxe8T0JcKkmQhGgh3N3d6/3H/cknn6DX6/n9999ZuHDhOZ/zySef4OnpyebNm3njjTd48cUXWbFihfX4jBkzeO2113j22WfZt28f//nPfwgJCblgHJ988glarZYtW7bw9ttvM3v2bD744APr8fLycpKSkti2bRsrV65ErVZz/fXXYzabL/oeCwsL+eWXX/j73/+Ou7t7vWOhoaHcdtttLF68uN4f5YtpTDwXuk9vv/02gwcP5t577yUrK4usrCwiIyMb/fp11w8KCmLLli089NBDPPDAA9xwww1cfvnlbN++nauvvprbb7+dioqK815Do9Hw6quvMn/+fE6cOGHT69srpot97WfNmsWnn37KwoUL2bt3L9OnT+evf/0ra9eurXeNi33fAqxfv57+/fvX27dz505UKhWzZ8/m2WefZefOnURFRfHUU0/VOy8hIYEtW7ZQXV19qbdJiPNThBDNLjExUbnuuusURVEUs9msrFixQjEYDMpjjz2mKIqiDBs2TOnbt2+D5w0bNkx55JFHrJ8PGTKk3vHLLrtMefLJJxVFUZSSkhLFYDAo77///nnjOPN6dY+7deummM1m674nn3xS6dat23mvkZeXpwDK7t27z3nNM23atEkBlKVLl57z+OzZsxVAycnJURRFUTp06KDMmTOn3jl9+vRRnnvuOZviudB9ulDMjXn9s69fW1ureHp6Krfffrt1X1ZWlgIoGzduPGfMZ34/DBo0SLnrrrsURVGUpUuXKmf/mnZUTBf72ldVVSkeHh7KH3/8Ue+17777buWWW26xXuNc37fnct1111nfZ52ZM2cq/v7+Sm5urnXfvHnzlB49etQ7b+fOnQqgHD16tFGvJURTSAuSEE6ybNkyvLy8cHNzY9y4cdx00031uknO/u/6XHr37l3vcVhYGLm5uQDs37+f6upqRo0aZVNcgwYNqjcwePDgwaSmpmIymQBITU3llltuoVOnTvj4+FgH2B4/frzRr6FcpIVIr9c3+lqNiedC98kezry+RqMhMDCQXr16WffVtdo15jVff/11PvnkE/bv39/sMV3oa3/48GEqKiq46qqr8PLysm6ffvppvS7NxnzfAlRWVuLm5lZv386dO7nuuuto166ddV96ejqxsbH1zqtrfbxQi5wQl0rr7ACEcFUjRozg3XffRa/XEx4ejlZb/8fR09PzotfQ6XT1HqtUKmvX0tldWPYyYcIEOnTowPvvv094eDhms5mePXtecEBundjYWFQqFfv37+f6669vcHz//v20a9cOPz8/ANRqdYNkymg02hzPhe7ThTTm9c93/TP31SUdjXnNoUOHMmbMGGbMmMGUKVNaREwAZWVlAPz4449ERETUO2YwGKyfN+b7FiAoKKjBgPSdO3cyY8aMevtSUlIYOnRovX11A/nPTKSEsDdpQRLCSTw9PYmNjSUqKqpBcmQPcXFxuLu7s3LlSpued/aA2E2bNhEXF4dGo6GgoICDBw/yj3/8g1GjRtGtWzebZl0FBgZy1VVX8a9//YvKysp6x7Kzs/niiy/qJQXt2rUjKyvL+rikpIT09HTr40uNp45er7e2kJ3pYq/vKK+99ho//PADGzdubNaYLvS17969OwaDgePHjxMbG1tvs3XMFkDfvn3Zt2+f9XFxcTFHjx6lb9++9c5LSUkhPj6+3r49e/bQvn17goKCbH5dIRpLEiQh2ig3NzeefPJJnnjiCWs3yKZNm/jwww8v+Lzjx4+TlJTEwYMH+fLLL5k/fz6PPPIIAP7+/gQGBvLvf/+bw4cPs2rVKpKSkmyK65133qG6upoxY8awbt06MjIyWL58OVdddRWdO3dm5syZ1nNHjhzJZ599xvr169m9ezeJiYloNBrrcXvEA5aZYZs3b+bo0aPk5+dbW1Uu9vqO0qtXL2677TbmzZvX4JgjY7rQ197b25vHHnuM6dOn88knn3DkyBG2b9/O/Pnz+eSTT2x+rTFjxrB3715rQrtr1y60Wm29bsBjx45x6tSpBgnS+vXrufrqq5v+RoVoBOliE6INe/bZZ9FqtcycOZOTJ08SFhbG/ffff8Hn3HHHHVRWVpKQkIBGo+GRRx7hvvvuAyzdO1999RUPP/wwPXv2pEuXLsybN4/hw4c3Oqa4uDi2bt3K888/z4033khubi6KojB58mQ+++wzPDw8rOfOmDGD9PR0rr32Wnx9fXnppZfqtZbYIx6Axx57jMTERLp3705lZSXp6elER0df9PUd6cUXX2Tx4sUN9jsypgt97QFeeukl2rVrx6xZs0hLS8PPz49+/frx9NNP2/xavXr1ol+/fnz99df87W9/Y+fOnXTp0qXeuKQdO3bg5+dXr5BkVVUV3377LcuXL7+k9yrExaiUi42WFEK4jOHDhxMfH9+gzo6jPffcc8yePZsVK1bUq4sj2rYff/yRxx9/nD179qBWN65D491332Xp0qX8+uuvDo5OuDppQRJCON0LL7xAdHQ0mzZtIiEhodF/LEXrds0115CamkpmZmajxzHpdDrmz5/v4MiEkBYkIcQZnNWCJIQQLY0kSEIIIYQQZ5F2bCGEEEKIs0iCJIQQQghxFkmQhBBCCCHOIgmSEEIIIcRZJEESQgghhDiLJEhCCCGEEGeRBEkIIYQQ4iySIAkhhBBCnEUSJCGEEEKIs0iCJIQQQghxlv8Hvq0KKX4l5DcAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots()\n", "\n", "n_list = np.arange(0, np.max([s.n for s in relevant_states]) + 1)\n", "rates_summed = {key: np.zeros(len(n_list)) for key in [\"dipole\", \"SP\", \"BBR\"]}\n", "\n", "for i, s in enumerate(relevant_states):\n", " rates_summed[\"dipole\"][s.n] += electric_dipole_moments[i]\n", "\n", " if i < len(transition_rates_sp):\n", " rates_summed[\"SP\"][s.n] += transition_rates_sp[i]\n", "\n", " if i < len(transition_rates_bbr):\n", " rates_summed[\"BBR\"][s.n] += transition_rates_bbr[i]\n", "\n", "for key, values in rates_summed.items():\n", " values_normalized = np.abs(values) / np.max(np.abs(values))\n", " ax.plot(n_list, values_normalized, label=key, alpha=0.8)\n", "\n", "ax.legend()\n", "\n", "ax.set_xlabel(\"Principal Quantum Number ($n$)\")\n", "ax.set_ylabel(r\"Transition Rates normalized to maximum\")\n", "\n", "ax.set_yscale(\"log\")\n", "ax.set_ylim(1e-10, 1e1)\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 }