{ "cells": [ { "cell_type": "markdown", "id": "878619db-6670-431a-ae94-b36492cc0565", "metadata": {}, "source": [ "# Tools" ] }, { "cell_type": "code", "execution_count": 1, "id": "ad68153a-c8e9-44ca-b6e8-a3cd4c54e51b", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkoAAAGwCAYAAABWwkp7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABagElEQVR4nO3deVxU9f4/8NfMMDOAbCKyo7grLuACpOWWJGJZaXXLvOZSdvNqtyK76bfVbmW/vJl5natlubTc9LZotyxTSUUTl1ByQVEMAWUTkG2AGZg5vz8GJonFmWGGMzO8no/HPIBzzpzznvE4vPh8PudzJIIgCCAiIiKiZqRiF0BERERkrxiUiIiIiFrBoERERETUCgYlIiIiolYwKBERERG1gkGJiIiIqBUMSkREREStcBG7AEej1+uRl5cHT09PSCQSscshIiIiEwiCgMrKSgQHB0MqNb2diEHJRCqVCiqVClqtFpcuXRK7HCIiIrJAbm4uQkNDTd5ewpm5zVNeXg4fHx/k5ubCy8tL7HKIiIjIBBUVFQgLC0NZWRm8vb1Nfh5blMzU2N3m5eXFoERERORgzB02w8HcRERERK1gUCIiIiJqBYMSERERUSsYlIiIiIhawaBERERE1AoGJSIiIqJWMCgRERERtYJBiYiIiKgVDEpEREREreiUQem7777DgAED0K9fP3z44Ydil0NERER2qtPdwqS+vh6JiYnYt28fvL29MXLkSEyfPh3dunUTuzQiIiKyM52uRenYsWMYPHgwQkJC4OHhgYSEBOzevVvssoiIiMgOOVxQSk5OxrRp0xAcHAyJRIIdO3Y020alUiE8PByurq6IjY3FsWPHjOvy8vIQEhJi/DkkJARXr17tiNJNoqnXoV6nF7sMIiIiggN2vanVakRGRmL+/PmYMWNGs/Xbtm1DYmIi1q9fj9jYWKxevRrx8fHIyMiAv7+/2cfTaDTQaDTGnysqKtpVf1uOZZXi0S3H4aF0gb+nEr8Vq3FrHz+8PC0CwT5uNjsuERERtczhWpQSEhLw+uuvY/r06S2uX7VqFRYsWIB58+YhIiIC69evh7u7OzZu3AgACA4ObtKCdPXqVQQHB7d6vBUrVsDb29v4CAsLs+4LusH6A5dQWVuP/PJa/HqlHJW19dh1tgBTVidj15l8mx2XiIiIWuZwQaktWq0WqampiIuLMy6TSqWIi4tDSkoKACAmJgZnzpzB1atXUVVVhR9++AHx8fGt7nPZsmUoLy83PnJzc21Wf2ZRFQBgRA8fLEsYiM8X3ILIUG9U1NbjiU9P4OENR5B0rhA6vWCzGoiIiOh3Dtf11pbi4mLodDoEBAQ0WR4QEIDz588DAFxcXPDOO+9g4sSJ0Ov1+Pvf/97mFW9KpRJKpdKmdQNAbZ0OuderAQDrZ4+Ev6crAODLhWPw7p4L2HDwNxy+VILDl0oQ4uOGOWN64s+39IS7wqn+CYmIiOxKp/wte/fdd+Puu+826zkqlQoqlQo6nc4mNeWUVkMQAE+lC7p7/B7M5DIp/j5lIGbG9MDHKZfxReoVXC2rwZvfn8cHyVlYNLEPZsb0gKtcZpO6iIiIOjOn6nrz8/ODTCZDYWFhk+WFhYUIDAxs174XLVqE9PR0HD9+vF37aU1xlWHAeIC3KyQSSbP1Yb7ueOHOCBxZNglv3zcMYb5uKK7SYPm36Zj4z/34z9Ec1PFqOSIiIqtyqqCkUCgwcuRIJCUlGZfp9XokJSVh9OjRIlZ2c+XVdQAAHzd5m9u5ymX4U3QYfnp2At6YPgSBXq7IL6/F/20/jUnvHMD2k1c4homIiMhKHC4oVVVVIS0tDWlpaQCArKwspKWlIScnBwCQmJiIDRs2YMuWLTh37hwWLlwItVqNefPmteu4KpUKERERiI6Obu9LaFFZTUNQcm87KDWSy6SYFdsT+5+bgJfvioCfhwI5pdV4ZtuvmLI6GTtP5UPPwERERNQuEkEQHOq36f79+zFx4sRmy+fMmYPNmzcDANauXYuVK1eioKAAUVFRWLNmDWJjY61y/IqKCnh7e6O8vBxeXl5W2ScArNt/Cf9v13ncNyIU7/wp0uznqzX12Hz4Mt4/cAkVtfUAgIGBnnhqUj/EDw6EVNq8O4+IiKizsPT3t8MFJbHZKiit+OEc3j/wGx69rRdeuivC4v2U19Rh46EsbDyUhUqNITANCvLCc/H9MXGAf4vjn4iIiJydpb+/Ha7rTSy27nozdYzSzXi7yfHMHf1x6Pnb8bfb+8JD6YJz+RWYv/kX3L8+BUd+K7FGuURERJ0CW5TMZKsWpSc+ScWuswX4xz2DMXt0uNX2W1atxboDl7D558vQ1BuuihvXvzuemzwAQ0O9rXYcIiIie8YWJQdXVqMFAHi7K6y6Xx93BZYlDELy3yfiz7f0gItUguQL1zBt7SH89bNUZBZVWvV4REREzoQtSmayVYtSRkElCipqMSjQE/5erlbb7x/llFTj3b0XsCPtKgQBkEqA+0aE4qm4fgjt6m6z4xIREYmJg7lt7MaZuS9cuGD1oNTRMgoq8c7uDOxON0zOqZBJ8XBsDzx5e19087D9LVuIiIg6EoNSB7FVi5JYTuZcx8ofM3D4kmGQdxeFDAvG9cZjY3vDQ9kp73BDREROiEGpgzhbUGp06GIx/t+u8zh9tRwA0K2LAotv74uHY3tA6cL7yBERkWNjUOogzhqUAECvF/DDmQL8c3cGsorVAIDQrm54dnJ/3BMZwkkriYjIYTEodRBnDkqN6nR6/PeXXLy39yKKKg036x0Y6Im/TxnASSuJiMghMSjZmLMN5jZFjVaHjT9nYf2BS6hsuC1KbC9fvHRXBIaEcA4mIiJyHAxKHaQztCj9UVm1Fuv2X8Kmw5ehrddDIgGmDw/Bc/EDEOTtJnZ5REREN8Wg1EE6Y1BqdLWsBit3nceOtDwAgKtcigVje+OJ8X3QhVfIERGRHWNQ6iCdOSg1+jW3DK/vTMfxy9cBAH4eSiyZ3B8PjAqDjAO+iYjIDjEodRAGJQNBEPDj2QKs+OE8skuqARgGfP/f1EEY17+7yNURERE1xaBkY51xMLcptPV6fHIkG2uSLqK8pg4AML5/d7xw5yD0D/AUuToiIiIDBqUOwhallpVVa7EmKROfHLmMOp0AqQR4KKYHnonrj+6evCUKERGJi0GpgzAote1ysRpv/XAeu84WADDcEuWvE/vi0dt6wVXOGb6JiEgcDEodhEHJNMeySvHGznT8esVwS5Rgb1c8N2UAZ/gmIiJRMCh1EAYl0+n1Av73ax7e3nUeeeW1AIDIMB8sv3swosJ8xC2OiIg6FQalDsKgZL7aOh0+OpSFf+/LhFqrAwA8MDIUz00ZAH9PV5GrIyKizoBBqYMwKFmuqKIWb+06j69PXAUAeChd8NSkfpgzJhwKF6nI1RERkTOz9Pc3fzuZSKVSISIiAtHR0WKX4rD8vVyx6k9R+PqvYzAs1BtVmnq88f05THkvGfszisQuj4iIqBm2KJmJLUrWodcL+PLEFby96zyKq7QAgLhB/njxzgiE+3URuToiInI27HrrIAxK1lVRW4c1ey9i8+HLqNcLUMikeHxcbyya2BduCk4nQERE1sGg1EEYlGwjs6gKr32XjuQL1wAAoV3d8Oq0wYiLCBC5MiIicgYco0QOra+/B7bMi8b6P49EsLcrrlyvwWMf/4LHtvyC3NJqscsjIqJOikGJ7IZEIsGUIYHY++x4PDG+D1ykEuw9V4g73j0A1b5MaOp1YpdIRESdDIMS2R13hQuWJgzED0+NxS29fVFbp8fKHzOQ8N5B/JxZLHZ5RETUiTAokd3qF+CJzxfcgtUPRsHPQ4nfrqkx68OjePLzkyisqBW7PCIi6gQYlMiuSSQS3Ds8BEnPjsfcMeGQSoBvf83DpHcOYOOhLNTr9GKXSEREToxXvZmJV72J68zVcryw4wx+zS0DAAwK8sKKGUN57zgiImoTr3qzMc7MbR+GhHhj+8IxWDFjKLzd5DiXX4EZ//4Zr32bDrWmXuzyiIjIybBFyUxsUbIfJVUavL7zHLafNNw7LsTHDa9PH4KJA/xFroyIiOwNW5So0+nmocS7D0Zhy/wYhHZ1w9WyGszbdBxPbT2JkiqN2OUREZETYFAihze+f3fsfmYcHrutF6QS4Ju0PExadQBfpl4BG0yJiKg9GJTIKbgrXPDiXRHYsehWDAryQll1HZZ88Stmf3QMOSWc2ZuIiCzDoEROZVioD/63+FY8P2UglC5SHMosxuTVB/BB8iVOJUBERGZjUCKnI5dJsXBCH/z49DiM7t0NtXV6vPn9eUz/92FkFFSKXR4RETkQBiVyWuF+XfCfBbF4+/5h8HaT4/TVckz71yGo9mWydYmIiEzCoEROTSKR4E+jwrDnmXGIG+QPrc5w37j71h3GxUK2LhERUdsYlKhT8PdyxYZHRuGdByLh5eqCX6+U485/HcL6A5eg0/PKOCIiahmDEnUaEokE940Mxe5nxmPigO7Q1uvx1g/ncf/6w8gsqhK7PCIiskMMStTpBHq7YuPcaLx9/zB4Kl1wMqcMU9ccxIbk39i6RERETXTKoDR9+nR07doV999/v9ilkEgaxy79+Mw4jOtvaF164/tz+NP7KfjtGluXiIjIoFMGpaeeegoff/yx2GWQHQj2ccOWedF4a8ZQeChdkJp9HVPXHMSnR7I5qzcREXXOoDRhwgR4enqKXQbZCYlEgodieuDHZ8bh1r6GeZde3HEGj235BcW8ZxwRUadmd0EpOTkZ06ZNQ3BwMCQSCXbs2NFsG5VKhfDwcLi6uiI2NhbHjh3r+ELJ6YT4uOGT+bF48c5BUMikSDpfhCmrk5F0rlDs0oiISCR2F5TUajUiIyOhUqlaXL9t2zYkJibilVdewYkTJxAZGYn4+HgUFRUZt4mKisKQIUOaPfLy8syuR6PRoKKiosmDnJdUKsFjY3vjf0/eioGBniiu0uLRLb/ghe2nUa2tF7s8IiLqYBLBjgdiSCQSbN++Hffee69xWWxsLKKjo7F27VoAgF6vR1hYGJ588kksXbrU5H3v378fa9euxZdfftnmdq+++iqWL1/ebHl5eTm8vLxMPh45nto6Hf75YwY+PJQFAOjt1wWrH4rCsFAfcQsjIiKzVVRUwNvb2+zf33bXotQWrVaL1NRUxMXFGZdJpVLExcUhJSXFJsdctmwZysvLjY/c3FybHIfsj6tchhfvisBnj8Ui0MsVvxWrMePfh6Hal8lpBIiIOgmHCkrFxcXQ6XQICAhosjwgIAAFBQUm7ycuLg4PPPAAvv/+e4SGhrYZspRKJby8vJo8qHO5ta8fdj09FncODUK9XsDKHzPwyMajKKqsFbs0IiKyMYcKStayd+9eXLt2DdXV1bhy5QpGjx590+eoVCpEREQgOjq6Ayoke+PjrsDah4dj5f3D4CaX4efMEkx97yAOXrwmdmlERGRDDhWU/Pz8IJPJUFjY9CqkwsJCBAYG2vTYixYtQnp6Oo4fP27T45D9kkgkeGBUGL69YaD3IxuP4e1d51Gn04tdHhER2YBDBSWFQoGRI0ciKSnJuEyv1yMpKcmkViEia+jr74kdi27FrNgeEATg3/sv4aEPjuBqWY3YpRERkZXZXVCqqqpCWloa0tLSAABZWVlIS0tDTk4OACAxMREbNmzAli1bcO7cOSxcuBBqtRrz5s2zaV3seqMbucpleGP6UKgeHgHPxhm93zuIH8+aPlaOiIjsn91ND7B//35MnDix2fI5c+Zg8+bNAIC1a9di5cqVKCgoQFRUFNasWYPY2NgOqc/SywvJeeWUVOPJz0/g1yvlAIC5Y8KxbOpAKF1kIldGRESNLP39bXdByd4xKFFLtPV6rPzxPDYcNMy5FBnmg3/PGoEQHzeRKyMiIqCTzKMkJna9UVsULlK8cGcENs4dBW83OX7NLcNdaw4i+QKviiMicmRsUTITW5ToZnJLq/HXz07g9NVySCTAM3H9sXhiX0ilErFLIyLqtNiiRGQnwnzd8cUTozEzJgyCAKzacwGPbjmOsmqt2KUREZGZGJSIbMBVLsOKGcPw9v3DoHSRYl/GNdz1r0M4c7Vc7NKIiMgMDEom4hglssSfRoXh67+OQQ9fd1y5XoMZ6w7jq9QrYpdFREQm4hglM3GMElmivLoOif9NQ9L5IgDAY7f1wtKEgXCR8W8VIqKOwDFKRHbM212ODY+Mwt9u7wsA+PBQFuZt5rglIiJ7x6BE1EGkUgkSJw/Av2eNgJtchoMXi3GP6mdcKKwUuzQiImoFgxJRB5s6NAhf/3UMQru6IbukGtNVP2M3b31CRGSXGJRMxMHcZE2Dgrzwv8W34ZbevlBrdXj8k1So9mWCQwaJiOwLB3ObiYO5yZrqdHq8/l06tqRkAwDuGxGKFTOGQuHCv2GIiKyJg7mJHJBcJsXye4bgH/cOgUwqwVcnrmD2R0c5yJuIyE4wKBHZgdm39MTGudHwULrgaFYppv/7MLKK1WKXRUTU6TEoEdmJ8f2748uFoxHi44asYjWm//tnHMsqFbssIqJOjUHJRBzMTR1hYKAXti8ag8hQb5RV12HWh0fw9QnO5E1EJBYO5jYTB3NTR6jR6pD43zT8cMYwbcDfpwzAwvF9IJFIRK6MiMgxcTA3kRNxU8igengEHh/XGwDw9q4MLP82HXo9/64hIupIDEpEdkoqleD/pg7Ci3cOAgBsPnwZT249CU29TuTKiIg6DwYlIjv32NjeWDNzOOQyCXaeysfcjcdRUVsndllERJ0CgxKRA7g7Mhib5sagi0KGlN9K8OD7R1BUUSt2WURETo9BichB3NbPD9v+Mhp+Hgqcy6/AjHWHcZlzLRER2RSDkok4PQDZgyEh3vhq4Rj07OaOK9dr8MD7KbhQWCl2WURETovTA5iJ0wOQPbhWqcHsj47ifEElurrL8fH8WAwN9Ra7LCIiu8XpAYg6ke6eSmx9/BZEhnrjenUdHt5wBL9c5izeRETWxqBE5KB83BX49LFYxIT7olJTj9kfHcPPmcVil0VE5FQYlIgcmKerHFvmx2BsPz/U1Okwb/Nx7E0vFLssIiKnwaBE5ODcFDJ8OGcUJkcEQFuvxxOfpmJXw61PiIiofRiUiJyA0kUG1awRuDsyGPV6AYv/cwJ72LJERNRuDEpETkIuk2LVnyIxrSEs/fWzVPx0nmGJiKg9GJSInIiLTIp3/xSJO4cGoU4n4IlPTmB/RpHYZREROSwGJSIn4yKTYvVDUUgYEgitTo/HP0lF8oVrYpdFROSQGJRMxJm5yZHIZVKsmTncOMB7wce/4NBFTh1ARGQuzsxtJs7MTY5EW6/HXz9Lxd5zRXCVS/Hpo7EYFe4rdllERB2OM3MTUTMKFylUs0ZgwoDuqK3TY/7m4ziXXyF2WUREDoNBicjJKV1kWDdrJEb17IqKWsMM3tklarHLIiJyCAxKRJ2Am0KGj+ZGY1CQF4qrNPjzR0dRWFErdllERHaPQYmok/B2k2PL/Gj07OaO3NIazP7oKMqqtWKXRURk1xiUiDoRf09XfPpoLAK8lLhQWIV5m4+jWlsvdllERHaLQYmokwnzdcfH82Ph7SbHyZwy/O3zk9DpefErEVFLGJSIOqEBgZ7YODcaShcp9p4rwvJvz4IzhRARNcegRNRJjezZFasfjIJEAnycko2PDmWJXRIRkd1hUCLqxBKGBuH/EgYBAN74/hx+OJ0vckVERPaFQYmok3tsbC88MronBAF4elsaUrOvi10SEZHd6HRBKTc3FxMmTEBERASGDRuGL774QuySiEQlkUjw8l0RmDTQH5qG+8LlllaLXRYRkV3odEHJxcUFq1evRnp6Onbv3o2nn34aajVnKabOzaXhJrpDQrxQqtZiwce/QK3htAFERJ0uKAUFBSEqKgoAEBgYCD8/P5SWlopbFJEd6KJ0wYZHRsHPQ4nzBZVY8sWvvBKOiDo9uwtKycnJmDZtGoKDgyGRSLBjx45m26hUKoSHh8PV1RWxsbE4duyYRcdKTU2FTqdDWFhYO6smcg5B3m5Y/+cRkMsk+OFMAf71U6bYJRERicrugpJarUZkZCRUKlWL67dt24bExES88sorOHHiBCIjIxEfH4+ioiLjNlFRURgyZEizR15ennGb0tJSPPLII/jggw/arEej0aCioqLJg8iZjQr3xev3DgEArNpzAT+eLRC5IiIi8UgEO25bl0gk2L59O+69917jstjYWERHR2Pt2rUAAL1ej7CwMDz55JNYunSpSfvVaDS44447sGDBAsyePbvNbV999VUsX7682fLy8nJ4eXmZ/mKIHMyr/zuLzYcvo4tChu2LbkX/AE+xSyIislhFRQW8vb3N/v1tdy1KbdFqtUhNTUVcXJxxmVQqRVxcHFJSUkzahyAImDt3Lm6//fabhiQAWLZsGcrLy42P3Nxci+snciQv3DkIo3t3g1qrw8JPUzm4m4g6JYcKSsXFxdDpdAgICGiyPCAgAAUFpnUP/Pzzz9i2bRt27NiBqKgoREVF4fTp061ur1Qq4eXl1eRB1BnIZVKsfXg4Ar1ccemaGku/Ps3B3UTU6biIXUBHu+2226DX681+nkqlgkqlgk6ns0FVRPapm4cSax8ejoc+OIJvf81DTHhXzB4dLnZZREQdpl0tSnV1dcjNzUVGRkaHXGLv5+cHmUyGwsLCJssLCwsRGBho02MvWrQI6enpOH78uE2PQ2RvRoX7YmnCQADAa9+lIy23TNyCiIg6kNlBqbKyEuvWrcP48ePh5eWF8PBwDBo0CN27d0fPnj2xYMECm4UJhUKBkSNHIikpybhMr9cjKSkJo0ePtskxiQh49LZeiB8cgDqdgEWfnUBZtVbskoiIOoRZQWnVqlUIDw/Hpk2bEBcXhx07diAtLQ0XLlxASkoKXnnlFdTX12Py5MmYMmUKLl68aHZBVVVVSEtLQ1paGgAgKysLaWlpyMnJAQAkJiZiw4YN2LJlC86dO4eFCxdCrVZj3rx5Zh/LHCqVChEREYiOjrbpcYjskUQiwcoHItGzmzuultVgyRenOF6JiDoFs6YHmDlzJl588UUMHjy4ze00Gg02bdoEhUKB+fPnm1XQ/v37MXHixGbL58yZg82bNwMA1q5di5UrV6KgoABRUVFYs2YNYmNjzTqOpSy9vJDIGZzNK8d01WFodXq8OX0oHo7tIXZJREQmsfT3t13Po2SPGJSos/vw4G94fec5uMql+O7Jsejr7yF2SUQtqq3Toay6DuU1dVBr61Gj1UGtqUdNnQ7VWsOjRluPaq0Omno9dHoBdbrGrwJ0ej3q9AJ0OgH1egESCSCVAFKJBFKJpOFnye/LpBIoXKRQukihdJEZvspv+N5FCqVcBne5DB6uLvBQusCz4auHqwuULjKx3zKnZunvb4uvequpqYEgCHB3dwcAZGdnY/v27Rg0aBDi4+Mt3a3d4lVvRAbzb+2F/RnXcCizGE9vO4mvF94KhYtDzTRCDkqvF1Ci1qKwohbXKjUorKhFYYUGRZW1KKnS4nq1FuU1dSirrkNZjRa1deZf4SwmuUxiDE0eSjm8XF3Q1V2Brl0U6Oouh28XRcPPcsPXhnVeri6QSCRil++0LG5Rmjx5MmbMmIEnnngCZWVlGDhwIORyOYqLi7Fq1SosXLjQ2rXaBbYoEQGFFbWIX52Msuo6PDG+j/GqOKL20OsFFFVqkHu9Gjkl1YavpdW4UlqD3OvVKKrUQKc371eWi1QCLzc5uihlcJe7wE0hg7tCBneFS8NXGdwUMihdZJDLJJBJJXCRSuAik8JF2vCzTApZQxDRCwIEQYBeMHyvF9Dws6HVqa5egKbe0EKlqddBU6f//ft6PTR1elRr61GlaXjU1kOtbd8f4HKZBH4eSvh7KtHd0xX+Xobv/T1dDV+9DN/7eSjgIuu8f9R0eIvSiRMn8O677wIAvvzySwQEBODkyZP46quv8PLLLzttUCIiIMDLFW/NGIYnPk3F+8mXML5/d4zu003ssshB6PQCckqrcbGwEheLqoxfL12rumkrkEQCYygI8HJFgJchHHT3UMDHXQEfd0Nri7ebHD7ucngo7b+1RacXoNYaQlOVph6VDV8NrWNalKq1KKuuQ6na0Gp2vVqL6+o6XK/WolqrQ51OQH55LfLLawGUt3ociQTo7qFEsI8bQnzcEOzjimAfN+PPIT5u8HGX2/371dEsDkrV1dXw9DTc+2n37t2YMWMGpFIpbrnlFmRnZ1utQCKyT1OGBOKh6DBsPZ6LZ/+bhh+fGQdPV7nYZZGdqa3TIaOgEqeuluPMlXKcvlqOzGtV0Na3HIhkUgmCvF3Rw9cdYV3d0aObO0K7uiHM1x3B3m5O2Soik0rg5SqHlwX/f2rrdChRa3GtUoOiiloUVWpQVKnBtcpaFFU0fq/BtSpDa1zj+tbmQ3OTy4wBKrSrO3p2c0dPX3f07NYFPbu5o4uy081TbXlQ6tu3L3bs2IHp06fjxx9/xDPPPAMAKCoqcsouKY5RImrupbsicPhSCXJKq/HWD+fxxvShYpdEIhIEQ0vRsaxSnMi5jlNXypFRUIn6FrrLlC5S9PX3QP8AT/T190A/fw/09fdAmK875E4WhGzJVS4ztga1Ra8XUFqtRX5ZLa6W1SCv8VFeg6tltcgrq8G1Sg1q6nS4dE2NS9fULe7Hz0NhCE2+hhAb3q0LejSEKd8uCqdsjbJ4jNKXX36Jhx9+GDqdDpMmTcLu3bsBACtWrEBycjJ++OEHqxZqLzhGiaipw5eK8fCGowCA/yyIxZg+fiJXRB1FEARkFFbiWFYpjmWV4vjlUhRWaJpt19VdjqGhPhga4oWhId4YFOSF0K7ukEmd75eqI9PU65DfEJqulNXgyvUaZJeokV1iGCtWqm57ollvNzn6+nugT/cuDV8N4dde/q1FmR6goKAA+fn5iIyMhFRq+Avg2LFj8PLywsCBzjm4k0GJqLkXtp/GZ0dz0MPXHbueHgt3Redrnu8srqu1OJRZjAMXriH5wjUUVTYNRnKZBENDvBHdyxfDw3wwJMQbIT5uTtnS0NlU1NYhp6Qa2SXVyC5VI7vY8DWnpBr5FbVoLU0oXKTo7dcFfbp7oM8fgpSrvOOmROA8Sh2EQYmoucraOsS/m4y88lrMv7UXXp4WIXZJZCWCIOBCYRV2nSnAvowinLpShht70lzlUozq6YuYXr6IDvdFVJgP3BScD6izqa3TIatYjcyGQfmGr2r8dq0KmlbGo0kkQA9f94ZuV0/08/dAvwBDK5Qt/tgSJSglJSUhKSkJRUVF0OubvhEbN260dLd2jUGJqGX7M4owd9NxSCTAF38ZjVHhvmKXRBYSBAFn8yrw/el87DpTgN+Km45XGRDgifEDumNcv+4YFd61Q1sFyLHo9ALyymr+EKAMX69X17X6vH/PGoGpQ4OsWkuHTw+wfPlyvPbaaxg1ahSCgoKcvlmVg7mJ2jZhgD/uHxmKL1Ov4P+2n8bOv43loFwHk1WsxtcnrmBH2lXkltYYlytcpBjXzw93RARgXP/uCPJue+AwUSOZVIIwX3eE+bpj4kD/JuuKqzS4WFiFi0WVTb6WqLXo4esuUsXNWdyiFBQUhLfffhuzZ8+2dk12jS1KRK27rtZi0qoDKFVrsSxhIP4yvo/YJdFNlNfUYeepfHx14gpSs68bl7vKpZg4wB8JQ4Nw+0B/eHTCy8JJHCVVGni5ya3+h1aHtyhptVqMGTPG0qcTkRPq2kWBpQkD8fcvT2H13ouYFhmM4JtctkziOHO1HB+nXMb/fs0zTvIolQDj+nfHjBGhiBvkz0H5JIpuHkqxS2jC4hal559/Hh4eHnjppZesXZNdY4sSUdv0egEPfpCC45evI35wAN6fPUrskqiBtl6Pnafz8HFKNk7mlBmX9w/wwH0jQnHv8BAEeLmKVyCRDXV4i1JtbS0++OAD7N27F8OGDYNc3nRG0VWrVlm6ayJyYFKpBK/fOxRT1xzEj2cL8dP5Qtw+MEDssjq1Kk09th7LwYcHs1BQUQvAcBl/wpAgPDK6J0b27Or040yJLGVxUDp16hSioqIAAGfOnGmyzhn/w3EwN5HpBgR64tHbeuGD5N+w/Nt03NrXD0oXXhnV0UrVWmz+OQtbUrJRXmO4wsjfU4nZt/TEQzE90N3Tvro4iOwR51EyE7veiEyj1tRj4j/3o6hSgxemDsKCcb3FLqnTqNLU48ODv+HDg1mo0tQDAHr5dcFfxvXG9BEhDK3UKVn6+9vsIeUvv/wyUlNTzX0aEXUyXZQuWBI/AACw5qeLKKlqfmsLsq7aOh0+PPgbxr29D6v3XkSVph4RQV5YN2sE9iaOx0MxPRiSiMxkdtfblStXkJCQAIVCgWnTpuHuu+/GpEmToFAobFEfETmw+0aEYvPPl5GeX4HVey/iH/cOEbskpyQIAn44U4A3dp7D1TLD/Ee9/bogcXJ/TB0SBKkd3GeLyFFZ1PWm1+vx888/49tvv8U333yD/Px83HHHHbjnnntw1113wdfXeWfkZdcbkXlSLpVg5oYjkEkl2PXUWPQL8BS7JKeSUVCJ5d+exeFLJQCAIG9XPB3XD/eNCIULJ/wkMhL1Xm/nzp0zhqbU1FTExMTg7rvvxsyZMxESEtLe3dsVBiUi8z3+8S/YnV6I8f27Y8v8GLHLcQo1Wh1W7cnAxp8vQ6cXoHSR4onxffDE+D681xpRCzpsjNLjjz+OqqqqJssGDRqEv//97/j555+Rm5uLOXPm4ODBg/j888/N3T0ROaFlUwfBRSrBgQvXcOS3ErHLcXgpl0ow5b1kbDiYBZ1ewJTBgdibOB7P3NGfIYnIysxuUZLJZMjPz4e/v//NN3YiN04PcOHCBbYoEZnphe2n8dnRHIzq2RVfPDHaKacRsbUqTT3e2HkOnx/LAWDoZntj+hDOU0Vkgg7repNKpSgoKOh0QakRu96ILFNQXovxK/dBU6/HpnnRmDigc36GWOr0lXI8+fkJXC6pBgD8+ZYeeH7KQHi6ym/yTCICOrDrjYjIEoHernhkdE8AwD9/zIBezyncTKHXC/jw4G+Yse5nXC6pRrC3Kz5fcAtev3coQxJRBzA7KEkkEjaZE5FFFk7oiy4KGc7mVWDX2QKxy7F7ZdVaPLrlOF7feQ51OsNYpB+eGofRfbqJXRpRp2H2PEqCIGDu3LlQKtue+v7rr7+2uCgick6+XRR49LZeWPNTJt7ZnYH4wYGQcY6fFl0orMSCj39Bdkk1lC5SvHRXBGbF9uAfqkQdzOygNGfOHFvUQUSdxGPjemNLSjYuXVPj+9P5mBYZLHZJdufHswVI3JYGtVaH0K5ueH/2SAwO9ha7LKJOyeygtGnTJlvUQUSdhJerHPNuDcfqvReh2peJu4YFsZWkgSAI2HDwN7z5/XkAwJg+3bD24RHw7cI7HxCJhYO5iajDzR0TDneFDOcLKrEvo0jscuyCXi/g9Z3njCFp7phwfDw/hiGJSGRmBaWcnByzdn716lWztieizsHHXYE/32K4Am7tT5mwwg0CHJqmXoentqXho0NZAIAXpg7Cq3cP5i1IiOyAWf8Lo6Oj8Ze//AXHjx9vdZvy8nJs2LABQ4YMwVdffdXuAonIOT12Wy8oZFKcyCnD0axSscsRjaZehyc+ScW3v+ZBLpNg9YNRWDCut9hlEVEDs8Yopaen44033sAdd9wBV1dXjBw5EsHBwXB1dcX169eRnp6Os2fPYsSIEXj77bcxdepUW9Xd4W6cmZuI2s/fyxUPjArFZ0dzoNqXiVt6d75L3mvrdFj4aSr2ZVyDq1yKDY+Mwth+3cUui4huYNFNcWtqarBz504cOnQI2dnZqKmpgZ+fH4YPH474+HgMGTLEFrXaBc7MTWQ9OSXVmPDPfdALwO5nxqF/gKfYJXWY2jodnvg0FfsbQtLGOdEY09dP7LKInFaH3cKks2NQIrKuv3zyC348W4iZMT2wYsZQscvpEPU6PZ749AT2nis0hKS50RjThyGJyJZ4CxMickjzb+0FAPj6xBVcV2tFrsb2BEHAS9+cwd5zhVC6MCQR2TsGJSISVUwvXwwO9oKmXo/Pj5t3Za0jei/pIj4/lgupBFgzczhDEpGdY1AiIlFJJBJjq9LHh7NRp9OLXJHtbD2Wg9V7LwIAlt8zBPGDA0WuiIhuhkGJiER3V2QQ/DyUKKioxQ9nnPNmuceySvHijjMAgMUT+2J2wzxSRGTfLA5Kc+bMQXJysjVrIaJOSukiw6zYHgCAz486X/dbfnkN/vpZKur1Au4aFoRnJ/cXuyQiMpHFQam8vBxxcXHo168f3nzzTc7CTUTt8qfoMEgkQMpvJcgqVotdjtXU1hkmlCyu0mJQkBfevn8Y721H5EAsDko7duzA1atXsXDhQmzbtg3h4eFISEjAl19+ibq6OmvWSESdQIiPG8b3N0y2uNWJBnX/47t0/HqlHD7ucnwweyTcFWbfi5yIRNSuMUrdu3dHYmIifv31Vxw9ehR9+/bF7NmzERwcjGeeeQYXL160Vp1E1Ak8FG3ofvsq9Qq09Y4/qPvHswX4rKErcc1DwxHm6y5yRURkLqsM5s7Pz8eePXuwZ88eyGQyTJ06FadPn0ZERATeffddaxyCiDqBSYP84eehRHGVFknnCsUup13yy2vw/FenAAB/Gdcb4/rz1iREjsjioFRXV4evvvoKd911F3r27IkvvvgCTz/9NPLy8rBlyxbs3bsX//3vf/Haa69Zs952Kysrw6hRoxAVFYUhQ4Zgw4YNYpdERA3kMikeGBUKAPj8eK7I1VhOrxeQuO1XlFXXYWiIN56dPEDskojIQhZ3lgcFBUGv12PmzJk4duwYoqKimm0zceJE+Pj4tKM86/P09ERycjLc3d2hVqsxZMgQzJgxA926db4bchLZo4eiw7Bu/yUcvHgNV8tqEOLjJnZJZvvkSDZSfiuBu0KG9x6KgsKFM7EQOSqL//e+++67yMvLg0qlajEkAYCPjw+ysrIsPYRNyGQyuLsbxgloNBoIggDe7o7IfvTs1gUxvXwhCMD/0vLELsdsV8tq8Pau8wCApQkD0bu7h8gVEVF7WByUxo8fD6VS2Wy5IAjIybH8ipXk5GRMmzYNwcHBkEgk2LFjR7NtVCoVwsPD4erqitjYWBw7dsysY5SVlSEyMhKhoaF47rnn4OfHWwgQ2ZPpw0MAAN+kOda0I4Ig4MXtp6HW6jCqZ1f8OZaTShI5OouDUq9evXDt2rVmy0tLS9GrVy+LC1Kr1YiMjIRKpWpx/bZt25CYmIhXXnkFJ06cQGRkJOLj41FUVGTcpnH80R8feXmGv059fHzw66+/IisrC//5z39QWOjYg0aJnM3UIUFQyKQ4X1CJc/kVYpdjsv/9mod9GdegkEnx1n3DIJVyviQiR2fxGCVBEFqcNK2qqgqurq4WF5SQkICEhIRW169atQoLFizAvHnzAADr16/Hzp07sXHjRixduhQAkJaWZtKxAgICEBkZiYMHD+L+++9vcRuNRgONRmP8uaLCcT60iRyVt7sctw/0x66zBdhx8ioGBXmJXdJNVdbW4R/fnQMAPHl7X/T1Z5cbkTMwOyglJiYCMNzI8qWXXjKO9wEAnU6Ho0ePtjpmqb20Wi1SU1OxbNky4zKpVIq4uDikpKSYtI/CwkK4u7vD09MT5eXlSE5OxsKFC1vdfsWKFVi+fHm7ayci89w7PBi7zhbgm7Q8PD9loN23zqj2XUJxlQa9/LrgL+P7iF0OEVmJ2UHp5MmTAAwtSqdPn4ZCoTCuUygUiIyMxJIlS6xX4Q2Ki4uh0+kQEBDQZHlAQADOnz9v0j6ys7Px+OOPGwdxP/nkkxg6dGir2y9btswYDgFDi1JYWJhlL4CITDZxoD+8XF1QUFGLY5dLcUtv+70yNbtEjY2HDBeuvHjnIF7lRuREzA5K+/btAwDMmzcP7733Hry87L9J/EYxMTEmd80BgFKphFKphEqlgkqlgk6ns11xRGSkdJFh8uBAfJl6BbvOFNh1UFrx/XlodXqM7eeH2wf6i10OEVmRxX/2bNq0qcNDkp+fH2QyWbPB14WFhQgMDLTpsRctWoT09HQcP37cpschot9NGWz4f/3j2QLo9fY5jcfJnOvYdbYAUgnw0l0RvOEtkZMxq0UpMTER//jHP9ClS5cm3VEtWbVqVbsKa4lCocDIkSORlJSEe++9FwCg1+uRlJSExYsXW/14RCSu2/r5oYtChvzyWpy6Wo6oMB+xS2pm1Z4LAIAZI0LRP8BT5GqIyNrMCkonT55EXV2d8fvWtOcvqqqqKmRmZhp/zsrKQlpaGnx9fdGjRw8kJiZizpw5GDVqFGJiYrB69Wqo1WrjVXC2wq43oo7nKpdh4kB/fHcqHz+cybe7oHTktxIcvFgMuUyCpyb1E7scIrIBiWBn01Lv378fEydObLZ8zpw52Lx5MwBg7dq1WLlyJQoKChAVFYU1a9YgNja2Q+qrqKiAt7c3ysvLHW58FpEj+u5UHhb/5yTCu7lj35IJdtO1JQgC/vR+Co5fvo7Zt/TEP+4dInZJRNQGS39/WxyUampqIAiCcXqA7OxsbN++HREREZg8ebIlu3QIDEpEHatKU48R/9gDbb0ePzw11m7mVDryWwke+uAIFC5SHPz7RAR4WT5/HBHZnqW/vy0ezH3PPffg448/BmC4JUhMTAzeeecd3HPPPVi3bp2lu7VbKpUKERERiI6OFrsUok7FQ+mCcf0MtxlKOmc/s+i/f+ASAOBPo0IZkoicmMVB6cSJExg7diwA4Msvv0RgYCCys7Px8ccfY82aNVYr0F7wqjci8UwYYLjkfn9G89smieF8QQX2ZVyDVAI8dltvscshIhuyOChVV1fD09Nwhcfu3bsxY8YMSKVS3HLLLcjOzrZagUREEwZ0BwCcyLmO8uo6kasBPkj+DQAwZUggwv26iFwNEdmSxUGpb9++2LFjB3Jzc/Hjjz8axyUVFRVx7A4RWVVoV3f08/eAXgCSL4rbqlRQXov/pRlusP2XcbxVCZGzszgovfzyy1iyZAnCw8MRGxuL0aNHAzC0Lg0fPtxqBdoLjlEiEldjq5LY3W9bj+egXi8gOrwrIu1sugIisr52TQ9QUFCA/Px8REZGQio1ZK5jx47By8sLAwcOtFqR9oRXvRGJ43BmMR7+8Cj8PJQ49n+TRLlJbr1Oj9v+3z4UVNTivYeicE9USIfXQESWsfT3t9n3ertRYGBgs1uHxMTEtGeXREQtGhXuiy4KGYqrNDibV4Ghod4dXsO+jGsoqKiFbxcFpgyx7W2TiMg+tCsoJSUlISkpCUVFRdDr9U3Wbdy4sV2FERHdSOEixZi+ftiTXogDF4pECUqfHTVcqPLAyFAoXWQdfnwi6ngWj1Favnw5Jk+ejKSkJBQXF+P69etNHkRE1ja2YT6llN9KOvzYV65X48AFw/iomTE9Ovz4RCQOi1uU1q9fj82bN2P27NnWrMdu8V5vROIb06cbAOCXy9ehqdd1aKvO9hNXIQiGGjglAFHnYXGLklarxZgxY6xZi13jhJNE4uvT3QPdPZXQ1OtxMqesw44rCAJ2pF0FAEwfzgHcRJ2JxUHpsccew3/+8x9r1kJE1CaJRILRvQ2tSimXOq777WxeBS5dU0PhIuUgbqJOxuKut9raWnzwwQfYu3cvhg0bBrlc3mT9qlWr2l0cEdEfje7TDf/7NQ8pl0rwzB0dc8z//WqYYDJukD88XeU32ZqInInFQenUqVOIiooCAJw5c6bJOomk4+c3IaLOoXGc0snc66jR6uCmsO04JZ1eMM7EzXmTiDofi4PSvn37rFmH3eNgbiL70MPXHQFeShRWaHDqShliG7ribOVYVikKKmrh5epinB2ciDoPi8coAcDBgwfx5z//GWPGjMHVq4aBjp988gkOHTpkleLsCQdzE9kHiUSCUT19AQCpObafiuSHM/kADDfA5dxJRJ2PxUHpq6++Qnx8PNzc3HDixAloNBoAQHl5Od58802rFUhE9EcjenYFAKRetm1QEgQBe9MLAQDxgzmIm6gzsjgovf7661i/fj02bNjQZCD3rbfeihMnTlilOCKiloxsDEo519GO21Xe1Nm8CuSV18JNLsOtff1sdhwisl8WB6WMjAyMGzeu2XJvb2+UlZW1pyYiojYNDvaCq1yKsuo6/Fasttlx9jS0Jo3r7wdXObvdiDoji4NSYGAgMjMzmy0/dOgQevfu3a6iiIjaIpdJMSzUBwCQmm277rfGoBQ3KMBmxyAi+2ZxUFqwYAGeeuopHD16FBKJBHl5efjss8+wZMkSLFy40Jo1EhE1M9LG45QKymuRnl8BiQS4faC/TY5BRPbP4ukBli5dCr1ej0mTJqG6uhrjxo2DUqnEkiVL8OSTT1qzRrvA6QGI7EtUmA8A4NTVcpvs/1BmMQBgaIg3unkobXIMIrJ/EqGdIyG1Wi0yMzNRVVWFiIgIeHh4WKs2u1RRUQFvb2+Ul5fDy8tL7HKIOq28shqMeesnuEglOLM83upjiJ7eehI70vLw1wl98PcpA626byLqeJb+/ra46y0nJweCIEChUCAiIgIxMTHGkJSTk2PpbomITBLk7YpuXRSo1ws4X1Bp1X0LgoBDmYZ7yd3Wj1e7EXVmFgelXr164dq1a82Wl5SUoFevXu0qiojoZiQSCQaHeAMATlu5++18QSWKqzRwk8uMY6GIqHOyOCgJgtDiPd2qqqrg6urarqKIiEwxNMTQfH7minWD0qGLhvFJMb18ORs3USdn9mDuxMREAIa/5l566SW4u7sb1+l0Ohw9etR4s1wiIlsaaqMWpSO/Gbrdbu1r2/vIEZH9MzsonTx5EoChRen06dNQKBTGdQqFApGRkViyZIn1KiQiasWQhqB0obASmnqdVVp/9HrBeA+5mF4MSkSdndlBad++fQCAefPmYc2aNfD09LR6UUREpgjxcUNXdzmuV9cho6DSOAlle1y6VoWy6jq4yqUYHMwrW4k6O4vnUdq0aROSkpKQlJSEoqIi6PX6Jus3btzY7uKIiNoikUgQEeyFnzNLcN5KQel4wwSWw8O6Qi6zeBgnETkJiz8FXnvtNUyePBlJSUkoLi7G9evXmzyIiDpC/wBDq/YFK00R8Et2KQAgOpxXuxFRO1qU1q1bh82bN2P27NnWrMducWZuIvtkDEpFVVbZ3y8NLUojw32tsj8icmwWtyhptVqMGTPGmrXYtUWLFiE9PR3Hjx8XuxQiuoE1W5SKKmqRU1oNqQQY0cOn3fsjIsdncVB67LHH8J///MeatRARma1/gOGOAAUVtSivrmvXvtJyyxr26QlPV3l7SyMiJ2Bx11ttbS0++OAD7N27F8OGDYNc3vRDZdWqVe0ujojoZjxd5QjxccPVshpcKKpEdDu6zM7kVQD4fX4mIiKLg9KpU6eME0ueOXOmybqWZuwmIrKVfgEeuFpWg4yCdgalhokrhzAoEVEDi4NS43xKRERiGxDgif0Z13CxsH3jlH4PSpw/iYgMOEkIETm8xgHd59sxoLuoohZFlRpIJcCgIAYlIjIwu0VpxowZJm339ddfm10MEZEl+jUM6L50TW3xPs7kGVqT+nT3gLvC4sZ2InIyZn8aeHuz756I7Eu4XxcAQHGVBpW1dRZdsXbmKgdyE1FzZgelTZs22aIOIiKLebnK4eehRHGVBpeLqzE01Pywc7phfNJgBiUiugHHKBGRU+jd0Kr0W7FlM3SfLzC0KEVwfBIR3YBBiYicQrifOwAgq9j8cUq1dTpcuV4D4PfxTkREQCcOStXV1ejZsyeWLFkidilEZAW9/AwB57IFQenStSoIAuDjLke3Lgprl0ZEDqzTBqU33ngDt9xyi9hlEJGV9GroerOkRSmz4Ya6/fw9OGEuETXRKYPSxYsXcf78eSQkJIhdChFZSWNQulxSbfZzG4NSX392uxFRU3YXlJKTkzFt2jQEBwdDIpFgx44dzbZRqVQIDw+Hq6srYmNjcezYMbOOsWTJEqxYscJKFRORPQjt6gYAKK+pQ0WteTfHbQxKfbozKBFRU3YXlNRqNSIjI6FSqVpcv23bNiQmJuKVV17BiRMnEBkZifj4eBQVFRm3iYqKwpAhQ5o98vLy8M0336B///7o379/R70kIuoAXZQuxvFFV0przHouW5SIqDV2N/1sQkJCm11iq1atwoIFCzBv3jwAwPr167Fz505s3LgRS5cuBQCkpaW1+vwjR45g69at+OKLL1BVVYW6ujp4eXnh5ZdfbnF7jUYDjUZj/LmiosKCV0VEHSG0qxtK1FrkXq9GRLBpl/nX6fS4XGIY18SgRER/ZHctSm3RarVITU1FXFyccZlUKkVcXBxSUlJM2seKFSuQm5uLy5cv45///CcWLFjQakhq3N7b29v4CAsLa/frICLbCPU1TBGQW2r6OKXskmrU6QS4K2QI9nazVWlE5KAcKigVFxdDp9MhICCgyfKAgAAUFBTY5JjLli1DeXm58ZGbm2uT4xBR+4V1NQSlxjmRTHHj+CSplFe8EVFTdtf11pHmzp17022USiWUSiVUKhVUKhV0Op3tCyMiizQO6L5y3fQWpcZut8ar5oiIbuRQLUp+fn6QyWQoLCxssrywsBCBgYE2PfaiRYuQnp6O48eP2/Q4RGS5MF/zW5Qau+l6dnO3SU1E5NgcKigpFAqMHDkSSUlJxmV6vR5JSUkYPXq0iJURkT1obFHKLa2GIAgmPSenISg1hiwiohvZXddbVVUVMjMzjT9nZWUhLS0Nvr6+6NGjBxITEzFnzhyMGjUKMTExWL16NdRqtfEqOFth1xuR/WscjK3W6lCpqYeXq/ymz2lsUerBoERELbC7oPTLL79g4sSJxp8TExMBAHPmzMHmzZvx4IMP4tq1a3j55ZdRUFCAqKgo7Nq1q9kAb2tbtGgRFi1ahIqKCnh7e9v0WERkGTeFDN5ucpTX1KGwvPamQUmnF4zddAxKRNQSuwtKEyZMuGmT+eLFi7F48eIOqoiIHEmQtyvKa+qQX16LfgGebW6bX16Der0AhUyKAC/XDqqQiByJQ41REpNKpUJERASio6PFLoWI2tAYeArKa2+6beP4pNCubpBxagAiagGDkol41RuRYwjybghKFTcPSo23OuFAbiJqDYMSETmVxhalfBNalK6UGYJSSFfOyE1ELWNQIiKn0tiiVGhCi1J+Q1AK9ub4JCJqGYOSiThGicgxBHib3qLU2D0XyHu8EVErGJRMxDFKRI7BOEap/OazczeGKbYoEVFrGJSIyKkENoxRul5dh9q6tieIbbwyLpBBiYhawaBERE7F200OucxwqX+pWtvqdhW1dajS1ANgUCKi1jEoEZFTkUgk6NZFCQAortK0ul1ja5K3mxzuCrube5eI7ASDkok4mJvIcXTzUAAASqpab1FqHJ8UxNYkImoDg5KJOJibyHH4eRhalK610aLUODUAgxIRtYVBiYicjjktShyfRERtYVAiIqfT3ePmY5Qa13X3ZFAiotYxKBGR0/m9RcmEoNSwLRFRSxiUTMTB3ESOw8/YotR611tjt1y3hm2JiFrCoGQiDuYmchzdzOh682NQIqI2MCgRkdPxa+hOa6tFqdjYosSuNyJqHYMSETmdxlaiUrUGer3QbH1tnc44KzdblIioLQxKROR0fLsYWon0AlBWU9dsfWO3m0ImhZcrZ+UmotYxKBGR05HLpPBUGgJQWXXz7reSG7rdJBJJh9ZGRI6FQYmInJKXmxxA2y1K7HYjopthUDIRpwcgciw+7oagVN5GUOJAbiK6GQYlE3F6ACLH4t3QolRe3VJQMnS9sUWJiG6GQYmInFJbLUolnBqAiEzEoERETqmxRamshRalxgHeXd0ZlIiobQxKROSUvN0MIailFqXGZY1hioioNQxKROSUjC1KNc2nB2BQIiJTMSgRkVNqHKNUwRYlImoHBiUickptjVFiUCIiUzEoEZFT8nFr/ao3BiUiMhWDEhE5pdZm5q6t00FTrwcAeLszKBFR2xiUTMSZuYkcy43zKAmCYFze2JoklQAeCt4Ql4jaxqBkIs7MTeRYGrvVtPV61Nbpjcsbg5KXmxxSKW+IS0RtY1AiIqfkoXSBrCEI3ThOieOTiMgcDEpE5JQkEgk8XQ1daxW1NwSlagYlIjIdgxIROa0uDWOQ1Jp64zK2KBGRORiUiMhpeSgbg5LOuOzGMUpERDfDoERETstdKQMAqLVsUSIiyzAoEZHT+r1FqXlQ8mFQIiITMCgRkdNqaYxSBbveiMgMDEpE5LS6NLYoaX8fo9TYDdfY2kRE1BYGJSJyWl0axyjd0KJU3RCaGtcREbWFQYmInFZji1LVDUGpMTS58/YlRGSCTvlJER4eDi8vL0ilUnTt2hX79u0TuyQisoHG7rXqG6YHMLYoMSgRkQk67SfF4cOH4eHhIXYZRGRD7gpD91rVDdMDNI5RcmfXGxGZgF1vROS0urQwPUBj6xJblIjIFHYXlJKTkzFt2jQEBwdDIpFgx44dzbZRqVQIDw+Hq6srYmNjcezYMbOOIZFIMH78eERHR+Ozzz6zUuVEZG9a6noztigp2KJERDdnd39SqdVqREZGYv78+ZgxY0az9du2bUNiYiLWr1+P2NhYrF69GvHx8cjIyIC/vz8AICoqCvX19c2eu3v3bgQHB+PQoUMICQlBfn4+4uLiMHToUAwbNszmr42IOpax662hRUmnF1Bbpwfwe2sTEVFb7O6TIiEhAQkJCa2uX7VqFRYsWIB58+YBANavX4+dO3di48aNWLp0KQAgLS2tzWOEhIQAAIKCgjB16lScOHGi1aCk0Wig0WiMP1dUVJjzcohIRMYWpYZWpOobxiqxRYmITGF3XW9t0Wq1SE1NRVxcnHGZVCpFXFwcUlJSTNqHWq1GZWUlAKCqqgo//fQTBg8e3Or2K1asgLe3t/ERFhbWvhdBRB3m9+kBDF1vjVe8yaQSKF0c6uOPiETiUJ8UxcXF0Ol0CAgIaLI8ICAABQUFJu2jsLAQt912GyIjI3HLLbfgkUceQXR0dKvbL1u2DOXl5cZHbm5uu14DEXWcP97r7fc5lGSQSCSi1UVEjsPuut5srXfv3vj1119N3l6pVEKpVEKlUkGlUkGn0938SURkFxq712rqdNDpBc6hRERmc6gWJT8/P8hkMhQWFjZZXlhYiMDAQJsee9GiRUhPT8fx48dtehwisp4bB2xXa+t/b1HiHEpEZCKHCkoKhQIjR45EUlKScZler0dSUhJGjx4tYmVEZI+ULlK4SA1dbGqNji1KRGQ2u/u0qKqqQmZmpvHnrKwspKWlwdfXFz169EBiYiLmzJmDUaNGISYmBqtXr4ZarTZeBWcr7HojcjwSiQTuChkqauuh1tZzDiUiMpvdBaVffvkFEydONP6cmJgIAJgzZw42b96MBx98ENeuXcPLL7+MgoICREVFYdeuXc0GeFvbokWLsGjRIlRUVMDb29umxyIi6/FQuhiCkqb+91m5OYcSEZnI7j4tJkyYAEEQ2txm8eLFWLx4cQdVRESOzLWh9ai2Ts8WJSIym0ONUSIiMpery+9XvnGMEhGZi0HJRCqVChEREW3OuURE9sfN2KKk41VvRGQ2BiUTcXoAIsfkKjd8zNWyRYmILMCgREROrbHrjS1KRGQJBiUTseuNyDHdOJibLUpEZC4GJROx643IMTUdzG1oUXLjVW9EZCIGJSJyajeOUaqt0zcsY1AiItMwKBGRU3OT/971Vltv6HpTuvCjj4hMw08LInJqrvLfB3Nr2KJERGZiUDIRB3MTOaYb51FqbFFyZYsSEZmInxYm4mBuIsfU2M1WwxYlIrIAgxIRObUbu95q6xrGKMn50UdEpuGnBRE5tcbB3DV1emjqG1qUXNiiRESmYVAiIqfWUosSu96IyFQMSibiYG4ix+SmMHzMqTX1qNcLADg9ABGZjp8WJuJgbiLH1NjNVlZd9/sytigRkYkYlIjIqSkbQlF5ze9BiS1KRGQqfloQkVNrHMxdpTHc500hk0IqlYhZEhE5EAYlInJqrn+YCoBTAxCROfiJQURO7Y/jkTg+iYjMwaBERE7N7Q/BiOOTiMgc/MQgIqfGFiUiag8GJRNxHiUix/THFqQ/jlkiImoLPzFMxHmUiByTVCppEpZ4+xIiMgeDEhE5vRu723jVGxGZg58YROT0bhzQzRYlIjIHgxIROb0bxyVxMDcRmYNBiYicXpOuN04PQERm4CcGETm9pmOU2KJERKZjUCIip9e0640fe0RkOn5iEJHTc2vS9cYWJSIyHYMSETm9G7ve2KJERObgJ4aJODM3keNqGpTYokREpmNQMhFn5iZyXE2CEq96IyIz8BODiJzejd1tvOqNiMzBoERETs+NY5SIyEL8xCAip3dj15sbW5SIyAwMSkTk9G5sRermoRSxEiJyNAxKROT0bmxR8vdkUCIi0zEoEZHTq9cJxu/9PV1FrISIHA2DEhE5vYraOuP3bgqOUSIi0zEoEZHT0ws334aIqCUMSkTk9B4Z3RN9unfBc/EDxC6FiByMi9gFiCErKwvz589HYWEhZDIZjhw5gi5duohdFhHZiJ+HEknPThC7DCJyQJ0yKM2dOxevv/46xo4di9LSUiiVvAqGiIiImut0Qens2bOQy+UYO3YsAMDX11fkioiIiMhe2d0YpeTkZEybNg3BwcGQSCTYsWNHs21UKhXCw8Ph6uqK2NhYHDt2zOT9X7x4ER4eHpg2bRpGjBiBN99804rVExERkTOxuxYltVqNyMhIzJ8/HzNmzGi2ftu2bUhMTMT69esRGxuL1atXIz4+HhkZGfD39wcAREVFob6+vtlzd+/ejfr6ehw8eBBpaWnw9/fHlClTEB0djTvuuMPmr42IiIgci90FpYSEBCQkJLS6ftWqVViwYAHmzZsHAFi/fj127tyJjRs3YunSpQCAtLS0Vp8fEhKCUaNGISwsDAAwdepUpKWltRqUNBoNNBqN8eeKigpzXxIRERE5KLvremuLVqtFamoq4uLijMukUini4uKQkpJi0j6io6NRVFSE69evQ6/XIzk5GYMGDWp1+xUrVsDb29v4aAxYRERE5PwcKigVFxdDp9MhICCgyfKAgAAUFBSYtA8XFxe8+eabGDduHIYNG4Z+/frhrrvuanX7ZcuWoby83PjIzc1t12sgIiIix2F3XW8d4WbdezdSKpVQKpVQqVRQqVTQ6XQ2ro6IiIjshUO1KPn5+UEmk6GwsLDJ8sLCQgQGBtr02IsWLUJ6ejqOHz9u0+MQERGR/XCooKRQKDBy5EgkJSUZl+n1eiQlJWH06NEiVkZERETOyO663qqqqpCZmWn8OSsrC2lpafD19UWPHj2QmJiIOXPmYNSoUYiJicHq1auhVquNV8HZCrveiIiIOh+JIAh2dV/t/fv3Y+LEic2Wz5kzB5s3bwYArF27FitXrkRBQQGioqKwZs0axMbGdkh9FRUV8Pb2Rnl5Oby8vDrkmERERNQ+lv7+trugZO8YlIiIiByPpb+/HWqMEhEREVFHsrsxSvaqcYxS461ROEM3ERGR42j8vW1uRxq73sx05coVzs5NRETkoHJzcxEaGmry9gxKZtLr9cjLy4OnpyckEonV9ltRUYGwsDDk5uZ2+rFPfC8M+D4Y8H34Hd8LA74PBnwfDEx9HwRBQGVlJYKDgyGVmj7yiF1vZpJKpWYlUXN5eXl16hP+RnwvDPg+GPB9+B3fCwO+DwZ8HwxMeR+8vb3N3i8HcxMRERG1gkGJiIiIqBUMSnZCqVTilVdegVKpFLsU0fG9MOD7YMD34Xd8Lwz4PhjwfTCw9fvAwdxERERErWCLEhEREVErGJSIiIiIWsGgRERERNQKBiUiIiKiVjAo2QmVSoXw8HC4uroiNjYWx44dE7skm1qxYgWio6Ph6ekJf39/3HvvvcjIyGiyzYQJEyCRSJo8nnjiCZEqto1XX3212WscOHCgcX1tbS0WLVqEbt26wcPDA/fddx8KCwtFrNh2wsPDm70XEokEixYtAuC850NycjKmTZuG4OBgSCQS7Nixo8l6QRDw8ssvIygoCG5uboiLi8PFixebbFNaWopZs2bBy8sLPj4+ePTRR1FVVdWBr6L92nof6urq8Pzzz2Po0KHo0qULgoOD8cgjjyAvL6/JPlo6h956660OfiXtd7NzYu7cuc1e55QpU5ps4+znBIAWPy8kEglWrlxp3MYa5wSDkh3Ytm0bEhMT8corr+DEiROIjIxEfHw8ioqKxC7NZg4cOIBFixbhyJEj2LNnD+rq6jB58mSo1eom2y1YsAD5+fnGx9tvvy1SxbYzePDgJq/x0KFDxnXPPPMMvv32W3zxxRc4cOAA8vLyMGPGDBGrtZ3jx483eR/27NkDAHjggQeM2zjj+aBWqxEZGQmVStXi+rfffhtr1qzB+vXrcfToUXTp0gXx8fGora01bjNr1iycPXsWe/bswXfffYfk5GQ8/vjjHfUSrKKt96G6uhonTpzASy+9hBMnTuDrr79GRkYG7r777mbbvvbaa03OkSeffLIjyreqm50TADBlypQmr/Pzzz9vst7ZzwkATV5/fn4+Nm7cCIlEgvvuu6/Jdu0+JwQSXUxMjLBo0SLjzzqdTggODhZWrFghYlUdq6ioSAAgHDhwwLhs/PjxwlNPPSVeUR3glVdeESIjI1tcV1ZWJsjlcuGLL74wLjt37pwAQEhJSemgCsXz1FNPCX369BH0er0gCJ3jfAAgbN++3fizXq8XAgMDhZUrVxqXlZWVCUqlUvj8888FQRCE9PR0AYBw/Phx4zY//PCDIJFIhKtXr3ZY7db0x/ehJceOHRMACNnZ2cZlPXv2FN59913bFtfBWnov5syZI9xzzz2tPqeznhP33HOPcPvttzdZZo1zgi1KItNqtUhNTUVcXJxxmVQqRVxcHFJSUkSsrGOVl5cDAHx9fZss/+yzz+Dn54chQ4Zg2bJlqK6uFqM8m7p48SKCg4PRu3dvzJo1Czk5OQCA1NRU1NXVNTk3Bg4ciB49ejj9uaHVavHpp59i/vz5TW4+3RnOhxtlZWWhoKCgyTng7e2N2NhY4zmQkpICHx8fjBo1yrhNXFwcpFIpjh492uE1d5Ty8nJIJBL4+Pg0Wf7WW2+hW7duGD58OFauXIn6+npxCrSx/fv3w9/fHwMGDMDChQtRUlJiXNcZz4nCwkLs3LkTjz76aLN17T0neFNckRUXF0On0yEgIKDJ8oCAAJw/f16kqjqWXq/H008/jVtvvRVDhgwxLn/44YfRs2dPBAcH49SpU3j++eeRkZGBr7/+WsRqrSs2NhabN2/GgAEDkJ+fj+XLl2Ps2LE4c+YMCgoKoFAomv0iCAgIQEFBgTgFd5AdO3agrKwMc+fONS7rDOfDHzX+O7f0+dC4rqCgAP7+/k3Wu7i4wNfX12nPk9raWjz//POYOXNmk5ug/u1vf8OIESPg6+uLw4cPY9myZcjPz8eqVatErNb6pkyZghkzZqBXr164dOkS/u///g8JCQlISUmBTCbrlOfEli1b4Onp2WxogjXOCQYlEt2iRYtw5syZJmNzADTpTx86dCiCgoIwadIkXLp0CX369OnoMm0iISHB+P2wYcMQGxuLnj174r///S/c3NxErExcH330ERISEhAcHGxc1hnOB7q5uro6/OlPf4IgCFi3bl2TdYmJicbvhw0bBoVCgb/85S9YsWKFU93m46GHHjJ+P3ToUAwbNgx9+vTB/v37MWnSJBErE8/GjRsxa9YsuLq6NllujXOCXW8i8/Pzg0wma3YlU2FhIQIDA0WqquMsXrwY3333Hfbt24fQ0NA2t42NjQUAZGZmdkRpovDx8UH//v2RmZmJwMBAaLValJWVNdnG2c+N7Oxs7N27F4899lib23WG86Hx37mtz4fAwMBmF37U19ejtLTU6c6TxpCUnZ2NPXv2NGlNaklsbCzq6+tx+fLljilQJL1794afn5/x/0JnOicA4ODBg8jIyLjpZwZg2TnBoCQyhUKBkSNHIikpybhMr9cjKSkJo0ePFrEy2xIEAYsXL8b27dvx008/oVevXjd9TlpaGgAgKCjIxtWJp6qqCpcuXUJQUBBGjhwJuVze5NzIyMhATk6OU58bmzZtgr+/P+688842t+sM50OvXr0QGBjY5ByoqKjA0aNHjefA6NGjUVZWhtTUVOM2P/30E/R6vTFMOoPGkHTx4kXs3bsX3bp1u+lz0tLSIJVKm3VDOZsrV66gpKTE+H+hs5wTjT766COMHDkSkZGRN93WonOiXUPBySq2bt0qKJVKYfPmzUJ6errw+OOPCz4+PkJBQYHYpdnMwoULBW9vb2H//v1Cfn6+8VFdXS0IgiBkZmYKr732mvDLL78IWVlZwjfffCP07t1bGDdunMiVW9ezzz4r7N+/X8jKyhJ+/vlnIS4uTvDz8xOKiooEQRCEJ554QujRo4fw008/Cb/88oswevRoYfTo0SJXbTs6nU7o0aOH8PzzzzdZ7sznQ2VlpXDy5Enh5MmTAgBh1apVwsmTJ41Xc7311luCj4+P8M033winTp0S7rnnHqFXr15CTU2NcR9TpkwRhg8fLhw9elQ4dOiQ0K9fP2HmzJlivSSLtPU+aLVa4e677xZCQ0OFtLS0Jp8ZGo1GEARBOHz4sPDuu+8KaWlpwqVLl4RPP/1U6N69u/DII4+I/MrM19Z7UVlZKSxZskRISUkRsrKyhL179wojRowQ+vXrJ9TW1hr34eznRKPy8nLB3d1dWLduXbPnW+ucYFCyE//617+EHj16CAqFQoiJiRGOHDkidkk2BaDFx6ZNmwRBEIScnBxh3Lhxgq+vr6BUKoW+ffsKzz33nFBeXi5u4Vb24IMPCkFBQYJCoRBCQkKEBx98UMjMzDSur6mpEf76178KXbt2Fdzd3YXp06cL+fn5IlZsWz/++KMAQMjIyGiy3JnPh3379rX4f2HOnDmCIBimCHjppZeEgIAAQalUCpMmTWr2/pSUlAgzZ84UPDw8BC8vL2HevHlCZWWlCK/Gcm29D1lZWa1+Zuzbt08QBEFITU0VYmNjBW9vb8HV1VUYNGiQ8OabbzYJD46irfeiurpamDx5stC9e3dBLpcLPXv2FBYsWNDsD2tnPycavf/++4Kbm5tQVlbW7PnWOickgiAIprc/EREREXUeHKNERERE1AoGJSIiIqJWMCgRERERtYJBiYiIiKgVDEpERERErWBQIiIiImoFgxIRERFRKxiUiIiIiFrBoEREDmP//v2QSCTNbhTcUZKSkjBo0CDodDqr7G/p0qV48sknrbIvIrINzsxNRHZpwoQJiIqKwurVq43LtFotSktLERAQAIlE0uE1jRw5EomJiZg1a5ZV9ldcXIzevXsjLS0NvXv3tso+ici62KJERA5DoVAgMDBQlJB06NAhXLp0Cffdd5/V9unn54f4+HisW7fOavskIutiUCIiuzN37lwcOHAA7733HiQSCSQSCS5fvtys623z5s3w8fHBd999hwEDBsDd3R33338/qqursWXLFoSHh6Nr167429/+1qS7TKPRYMmSJQgJCUGXLl0QGxuL/fv3t1nT1q1bcccdd8DV1dW47NVXX0VUVBQ++eQThIeHw9vbGw899BAqKyuN23z55ZcYOnQo3Nzc0K1bN8TFxUGtVhvXT5s2DVu3brXOG0dEVsegRER257333sPo0aOxYMEC5OfnIz8/H2FhYS1uW11djTVr1mDr1q3YtWsX9u/fj+nTp+P777/H999/j08++QTvv/8+vvzyS+NzFi9ejJSUFGzduhWnTp3CAw88gClTpuDixYut1nTw4EGMGjWq2fJLly5hx44d+O677/Ddd9/hwIEDeOuttwAA+fn5mDlzJubPn49z585h//79mDFjBm4c8RATE4MrV67g8uXLFr5bRGRLLmIXQET0R97e3lAoFHB3d0dgYGCb29bV1WHdunXo06cPAOD+++/HJ598gsLCQnh4eCAiIgITJ07Evn378OCDDyInJwebNm1CTk4OgoODAQBLlizBrl27sGnTJrz55pstHic7O9u4/Y30ej02b94MT09PAMDs2bORlJSEN954A/n5+aivr8eMGTPQs2dPAMDQoUObPL9xn9nZ2QgPDzf9TSKiDsGgREQOzd3d3RiSACAgIADh4eHw8PBosqyoqAgAcPr0aeh0OvTv37/JfjQaDbp169bqcWpqapp0uzUKDw83hiQACAoKMh4rMjISkyZNwtChQxEfH4/Jkyfj/vvvR9euXY3bu7m5ATC0jBGR/WFQIiKHJpfLm/wskUhaXKbX6wEAVVVVkMlkSE1NhUwma7LdjeHqj/z8/HD9+nWTjt94LJlMhj179uDw4cPYvXs3/vWvf+GFF17A0aNH0atXLwBAaWkpAKB79+6mvFwi6mAco0REdkmhUFhtvqIbDR8+HDqdDkVFRejbt2+TR1vdfMOHD0d6errZx5NIJLj11luxfPlynDx5EgqFAtu3bzeuP3PmDORyOQYPHmzR6yEi22KLEhHZpfDwcBw9ehSXL1+Gh4cHfH19rbLf/v37Y9asWXjkkUfwzjvvYPjw4bh27RqSkpIwbNgw3HnnnS0+Lz4+Hlu2bDHrWEePHkVSUhImT54Mf39/HD16FNeuXcOgQYOM2xw8eBBjx441dsERkX1hixIR2aUlS5ZAJpMhIiIC3bt3R05OjtX2vWnTJjzyyCN49tlnMWDAANx77704fvw4evTo0epzZs2ahbNnzyIjI8Pk43h5eSE5ORlTp05F//798eKLL+Kdd95BQkKCcZutW7diwYIF7Xo9RGQ7nJmbiMhEzz33HCoqKvD+++9bZX8//PADnn32WZw6dQouLmzgJ7JHbFEiIjLRCy+8gJ49exoHa7eXWq3Gpk2bGJKI7BhblIiIiIhawRYlIiIiolYwKBERERG1gkGJiIiIqBUMSkREREStYFAiIiIiagWDEhEREVErGJSIiIiIWsGgRERERNQKBiUiIiKiVvx/0/t6dIrqn5IAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "iron = nx.Material.Template(nx.lib.material.SS_enriched)\n", "\n", "layer_Fe = nx.Layer(id = \"Fe layer\",\n", " material = iron,\n", " thickness = 3000)\n", "\n", "site = nx.Hyperfine(magnetic_field = 0,\n", " isotropic = True)\n", "\n", "iron.hyperfine_sites = [site]\n", "\n", "sample = nx.Sample(layers = [layer_Fe])\n", "\n", "beam = nx.Beam(polarization = 0)\n", "\n", "air = nx.Air(2)\n", "\n", "exp = nx.Experiment(beam = beam,\n", " objects = [air, sample],\n", " isotope = nx.lib.moessbauer.Fe57)\n", "\n", "# define the time spectrum measurement object\n", "time_spectrum = nx.TimeSpectrum(experiment = exp,\n", " time_length = 170, # in ns\n", " time_step = 0.2) # in ns\n", "\n", "# calculate the time spectrum from the object\n", "# a calculation can also be performed by the () operator instead of .Calculate()\n", "# the time spectrum will return the time axis and the intensity\n", "time_axis, intensity = time_spectrum()\n", "\n", "plt.semilogy(time_axis, intensity)\n", "plt.xlabel('time (ns)')\n", "plt.ylabel(r'Intensity ($\\Gamma$/ns)')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 2, "id": "35750cfb-b02a-4834-90b8-d7f1241efb35", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkYAAAGwCAYAAABM/qr1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABQiUlEQVR4nO3deXxU5d0+/mtmksm+kH0hK1sIhCQsiWBRkcimFkQqUqq44VMbtDbSFmvFolV8akXU5gfWbwH1UXGpYtWKQAQBDWENWwiEJCQh+76SSWbm/v0xyUDMwmQykzPL9X698mozc3LO52Y45PI+9yITQggQEREREeRSF0BERERkKRiMiIiIiLowGBERERF1YTAiIiIi6sJgRERERNSFwYiIiIioC4MRERERURcHqQuwNlqtFmVlZfDw8IBMJpO6HCIiIjKAEALNzc0ICQmBXN5/vxCD0SCVlZUhLCxM6jKIiIjICCUlJRg5cmS/7zMYGSg9PR3p6elQq9UAdH+wnp6eEldFREREhmhqakJYWBg8PDwGPE7GLUEGp6mpCV5eXmhsbGQwIiIishKG/v7m4GsiIiKiLgxGRERERF0YjIiIiIi6MBgRERERdWEwIiIiIurCYERERETUhcGIiIiIqAuDEREREVEXBiMiIiKiLgxGRERERF0YjAyUnp6O2NhYTJs2TepSiIiIyEy4V9ogca80IiIi68O90oiIiMiqtHdqcL6iGe2dGslqcJDsykRERGR3tFqB0oYrKKhpRWF1i+5/a1pRUN2KssYrEAL47DczMDl8hCT1MRgRERGRybWq1CiobkVBTQvyq1qQX92K/OoWFNa0QqXW9vtzns4OqG/tGMZKe2IwIiIiIqMIIVDVrOoKPi24eE0AKm9s7/fnHBUyRPi6IcrPDdH+bhjl544ofzdE+7nBx00JmUw2jK3oicGIiIiIBqTRCpTUtSGvShd+Lla14GJ1CwqqWtCsUvf7c37uSkT7u2OUvxtG+bsj2t8N0X7uGDnCBQ4KyxzmzGBEREREAIBOjRZFta3Iq2xBXlXXV2UzCmpa0dHP4y+5DIjwddOFnwB3jPLv/nKDt6tymFswdAxGREREdkat0eJSbRsuVDbjQmWzPgAV1rSiU9P3Kj5ODnJE+7tjTIA7Rl/zFeHrCicHxTC3wHwYjIiIiGxU9wyw8xXNOF/ZjPMVuiBUUN2KDk3fPUCuSkVX+PHAmMCrQWjkCFco5NKN/RkuDEZEREQ2oK61A7kVTboQVNGM3Ipm5FU2o7Wj7zWBugPQmEAPjA3U/e+YAHeEeLlAbgcBqD8MRkRERFakQ63FxaoWfQg6V9GM3PImVDWr+jxeqZBjVIA7xnWFn5ggD4wN9ECot30HoP4wGBEREVmo2hYVzpU3I6e8Ebnlzcgpb0J+dUu/44DCfVwxLkgXfrr/N9LXzWJngFkiBiMiIiKJabUCRXVtOFvWiJyyJuSUN+FceRMqm/ruBfJwdkBMkAfGB3siJsgTMcEeGBfoATcn/lofKv4JEhERDSOVWoO8yhacLWvE2bImnC3ThaC2fsYCRfq6Ynyw5zVfusdgUi6CaMsYjIiIiMykrUONc+VNOFPahLNljThT2oQLlc1Qa3s/CnNykCMmyAOxIboANCHEE+OCPOHOXqBhxT9tIiIiE2hRqZFT1oTTpY04U9qI06WNyK9ugehjOJCXiyMmhHh2fXlhQognovw4FsgSMBgRERENUluHLgSduqwLQKcuN6CgprXPEBTg4YSJoV6YGOKJ2K4QNHIEH4VZKgYjA6WnpyM9PR0aTd/PgImIyDap1Brkljfj1OUGnLrciFOXG5FX1Yw+noYhyNMZE0O9EBfqhbiRnpgY4oUAT+fhL5qMJhOir3xL/WlqaoKXlxcaGxvh6ekpdTlERGRCWq1AQU0rTpY04OTlBpwsacC58uY+V4kO9HRCXKg3Jo3UBaGJoV7w93CSoGoyhKG/v9ljREREdqu2RYXskgacKG5AdlcYam7vvVv8CFdHTBrpjfiRXogbqQtDgewJskkMRkREZBc6NVrklDXhRHE9TnSFoeK6tl7HOTnIERfqhfgwb8SHeSNhpDfCfDgmyF4wGBERkU2qblbheHE9jhfV43hxPU5dboRK3fuR2JgAd10A6voaF+QBR84Os1sMRkREZPW0WoELVc04ekkXhI4V16OotndvkJeLIxLDvTE5fAQSw70xaaQ3vFwcJaiYLBWDERERWZ0rHRpklzTg6KU6HO3qEfrp2CCZDBgb4IHJEd5IDB+ByeEjEO3nxo1TaUAMRkREZPFqW1Q4WlSPo5fqcORSPc6UNvZaPdpVqUBiuDemRPhgSsQIJISxN4gGj8GIiIgsTmnDFRwurMXhwnocuVSHi1UtvY4J8nTG1MgRmBoxAlMjfRAT5MGVo2nIGIyIiEhSQghcqm3D4cJaZBXUIauwDqUNV3odNybAHdOifDAtcgSmRfpwI1UyCwYjIiIaVkIIFNa04lBBHQ4V1CKrsBaVTaoexyjkMkwM8URSlA+mReq+RrgpJaqY7AmDERERmZUQAkW1bcgsqEVmfi0OFdSiqrlnEFIq5IgP80JylC+SonwwOWIEd5UnSfBvHRERmVxFYzt+uFiDH/NrkZlfg7LG9h7vKx3kSAzzxg3RvkiO9sHk8BFwdlRIVC3RVQxGREQ0ZI1tncgsqMHBizX48WItCmpae7zvqJAhMWwEbhjli+nRvkgM92YQIovEYERERIOmUmtwrKgeB/Nq8MPFGpwubeyx27xcBsSN9MaMUb6YMcoXUyN84KJkECLLx2BERETXJYTAxaoW7M+rwYG8amQV1OFKp6bHMaMD3PGz0X64cbQfkqN94OnMNYTI+jAYERFRnxraOnDwYg32X6jGgbwalP9knJCfuxNmjvHTh6EgL+42T9aPwYiIiADo9hs7VdqI789X4/sLVcguaejxeMzJQY6kKB/cNMYfM8f6YVygB9cRIpvDYEREZMfqWzuwP68ae3OrsD+vBnWtHT3eHxvojpvG+OOmsf5IivLhgGmyeQxGRER2RAiBs2VN2He+Ct/l9u4V8nBywI2j/XDLOF0YCvF2ka5YIgkwGBER2bgrHRr8cLEGGblV2JtbhYqmnmOFYoI8MCsmALeM9cfkiBFw5H5jZMcYjIiIbFBFYzsyciuRca4KP1ysgUqt1b/n4qjAjaP9cGtMAG4Zx14homsxGBER2QAhBHLKm7A7pxJ7zlXiTGlTj/dDvV0we3wAbo0JwA3RvhwrRNQPBiMiIivVqdHiSGEdduVUYndOZY8d6WUyIDHMG7PHByJlfCDGBrpzBhmRARiMiIisyJUODfbnVePbsxXIOFeFxiud+vecHeWYOcYft8UG4taYAPi5O0lYKZF1srtg1NDQgJSUFKjVaqjVavz2t7/FypUrpS6LiKhfTe2d2JtbhW9OV2DfhSq0d14dL+TjpsTsmADMmRCEn43247YbRENkd8HIw8MD+/fvh6urK1pbWzFx4kQsXrwYvr6+UpdGRKTX0NaBXWcr8c2ZcvxwsRYdmqthKNTbBXMnBGHuhEBMjfSBQs5HZESmYnfBSKFQwNXVFQCgUqkghIAQ4jo/RURkfnWtHdh1tgJfny5HZn4t1NcsMDTK3w3zJwZj3sQgTAjx5HghIjOxuGC0f/9+vPLKKzh27BjKy8vx+eefY9GiRT2OSU9PxyuvvIKKigrEx8fjzTffRFJSksHXaGhowM0334y8vDy88sor8PPzM3EriIgM09jWiW/PVuDLU2X4Mb8WmmvCUEyQBxbEBWP+xCCMCfSQsEoi+2Fxwai1tRXx8fF46KGHsHjx4l7vf/TRR0hLS8PmzZuRnJyMjRs3Yu7cuTh//jwCAgIAAAkJCVCr1b1+dteuXQgJCYG3tzdOnjyJyspKLF68GEuWLEFgYKDZ20ZEBAAtKjV251Tgq5Pl2J9XjU7N1TAUG+yJ2yfpwlC0v7uEVRLZJ5mw4OdIMpmsV49RcnIypk2bhn/84x8AAK1Wi7CwMDz++ONYs2bNoK/xm9/8BrfeeiuWLFnS5/sqlQoqlUr/fVNTE8LCwtDY2AhPT89BX4+I7JNKrcH356vxn5Nl2HOusscA6nGBHrhjUjDuiA9BlJ+bhFUS2a6mpiZ4eXld9/e3xfUYDaSjowPHjh3D008/rX9NLpcjJSUFmZmZBp2jsrISrq6u8PDwQGNjI/bv34/HHnus3+PXr1+PdevWDbl2IrI/Wq3A4Ut1+CK7FF+fKkdT+9We7EhfV/w8PgR3xIdgLB+TEVkMqwpGNTU10Gg0vR57BQYGIjc316BzFBUV4dFHH9UPun788ccRFxfX7/FPP/000tLS9N939xgREfXnQmUzPj9Rii9OlKKs8eq+ZIGeTrhzUgh+nhCCuFAvDqAmskBWFYxMISkpCdnZ2QYf7+TkBCcnLpJGRAOraVHhP9ll+OzE5R7bcXg4OWBBXDAWJoYgOcqXU+uJLJxVBSM/Pz8oFApUVlb2eL2yshJBQUESVUVE9kql1uC7c1X49Nhl7LtQrZ9R5qiQ4eaxAVg8ORS3xgRwXzIiK2JVwUipVGLKlCnIyMjQD8jWarXIyMjAqlWrzHrt9PR0pKenQ6PRmPU6RGT5zpQ24tNjl7EjuxQNbVe35Igf6YW7p4zEHZNC4OOmlLBCIjKWxQWjlpYWXLx4Uf99YWEhsrOz4ePjg/DwcKSlpWHFihWYOnUqkpKSsHHjRrS2tuLBBx80a12pqalITU3Vj2onIvtS39qBHdml+PjoZZwrv/qoLNDTCYsnj8Tdk0MxOoCDqImsncUFo6NHj2LWrFn677sHPq9YsQLbtm3D0qVLUV1djbVr16KiogIJCQnYuXMn1yEiIpPTagV+yK/BR0dKsOtspX5bDqVCjtsmBOIXU0Zi5hh/jhsisiEWvY6RJTJ0HQQisl6VTe345GgJth8pweX6K/rXY4M9sXRaGBYmhMDblY/KiKyJTa5jRERkLhqtwP4L1fjgcDG+y63SD6T2cHbAooRQLJ0WhomhfIxOZOsYjAzEwddEtqmquR2fHL2MD7KKUdpwtXdoWuQI3DstHAviguGi5KwyInvBR2mDxEdpRNZPCIHMglq8f6gY356t0O9i7+XiiLsnj8SypDBu2kpkY/gojYjoJ1pUanx+/DLezSxCXlWL/vXJ4d5YnhyB2ycFc80hIjvHYERENi+/ugXv/ngJ/z5eihaVbr8yV6UCixJD8avkCMSGsPeXiHQYjIjIJmm1AvvzqrH1h0v4/kK1/vVoPzfcNz0Cd08ZCU9nRwkrJCJLxGBERDalrUONfx8vxdYfClFQ3QoAkMmA2TEBWDEjEjeO8oOc6w4RUT8YjAzEWWlElq2yqR3v/HgJ72cVo/GKbpsODycH/GJqGFbMiECEr5vEFRKRNeCstEHirDQiy3KuvAlv7y/Al6fK0KnR/XMW4euKB2dEYsnUMLg78b//iIiz0ojIhnVPt3/r+4Ie44eSIn3w8MwopIwP5DYdRGQUBiMishoarcDOMxV4a38+Tl1uBADIZcD8uGCsnBmNhDBvaQskIqvHYEREFq9DrcWOE6XY/H0+Cmp0A6qdHeW4Z2oYHvlZNMJ9XSWukIhsBYMREVmsKx0afHi4GG8fKEB5YzsA3erUK2ZE4oEZkfBx40auRGRaDEZEZHFaVWr836EivH2gADUtHQCAAA8nrJwZjWXJ4RxQTURmw39dDMTp+kTm19zeiXczi/D/DhSgvk035T7MxwWP3Twad08JhZMDt+sgIvPidP1B4nR9ItNrUamx7YdCvH2gUL8GUaSvK1bdOgYLE0LgqJBLXCERWTtO1ycii9eqUuOdzEt4e//VHqJofzc8cesY3DEpGA4MREQ0zBiMiGjYtXdq8H+HirBpXz5qW3VjiKL93PDE7DG4Mz6EaxARkWQYjIho2HRqtPj02GW8kZGnn2UW4euKJ7oembGHiIikxmBERGan1Qp8dbocG3adx6XaNgBAiJcznpg9BndPGckxRERkMRiMiMisDubV4OWd53CmtAkA4OumROqs0fhlcjicHTnLjIgsC4MREZnFmdJG/O/OXBzIqwEAuDs54H9uisZDP4uCG9chIiILxX+dDMR1jIgMU9ZwBa98ex6fnygFADgqZPjVDRFYNWs0fN2dJK6OiGhgXMdokLiOEVHfWlRqbN6Xj7cPFECl1gIAFiaEYPWccQjz4V5mRCQtrmNERMNCoxX4+GgJXt11ATUtKgBAUpQP/nz7eEwa6S1tcUREg8RgRERGyyqoxbovc5BTrhtYHeXnhjXzYzAnNhAyGdciIiLrw2BERIN2ub4N67/JxdenygEAns4OeDJlLH51QwSUDpx6T0TWi8GIiAzW3qnB5u/zsWlfPlRqLeQy4JfJ4Ui7bRx83JRSl0dENGQMRkRkkIxzlVj3ZQ6K63QLNN4Q7YO1d0xAbAgnIRCR7WAwIqIBldS1Yd2XZ7HnXBUAINDTCX++PRZ3TArmOCIisjkMRkTUp06NFm8fKMDre/KgUmvhIJfh4Z9F4fHZY+DOBRqJyEbxXzci6uXopTr86fPTuFDZAgCYHu2LFxZNwOgAD4krIyIyLwYjA3Hla7IHjVc68fI35/Dh4RIAgI+bEn++fTzuSgzlYzMisgtc+XqQuPI12apvz1bg2R1nUNWsW6Rx6dQwrJkfgxGcbUZENoArXxORQaqbVfjLf87i69O6NYmi/dywfnEckqN9Ja6MiGj4MRgR2SkhBHZkl2LdlzloaOuEQi7D/9wUjSdmj4Gzo0Lq8oiIJMFgRGSHqptV+NPnp7E7pxIAEBvsib8tmYSJoV4SV0ZEJC0GIyI78/Wpcvx5x2nUt3XCUSHDE7eOwa9vGQVHBbfyICJiMCKyEw1tHXj2i7P48mQZAGB8sCde/UU8V64mIroGgxGRHfjhYg2e+vgkKpraoZDLkHrLKKy6dQw3fCUi+gkGIyIbplJr8Pdvz+PtA4UAdDPOXluagPgwb2kLIyKyUAxGRDbqQmUznvjwBHIrmgEAy5PD8czt4+Gq5G1PRNQf/gtJZGOEENh+pAR/+c9ZqNRa+Lgp8b93T8JtsYFSl0ZEZPEYjIhsSHN7J57+7DS+OqVbrHHmGD+8ek88AjycJa6MiMg6MBgZiHulkaU7dbkBqz44geK6NijkMvx+7jg8OjMacjn3OCMiMhT3Shsk7pVGlkYIgXczi/DXr3PQqREI9XbBG8sSMSVihNSlERFZDO6VRmQH2jrUWPPv0/hP19pEcycE4m93x8PL1VHiyoiIrBODEZGVyq9uwWP/dwwXKlugkMvw9PwYPPyzKMhkfHRGRGQsBiMiK/TN6XKs/uQkWjs0CPBwwj9+ORlJUT5Sl0VEZPUYjIisiFYr8Oru80jfmw8ASI7ywZu/TOSsMyIiE2EwIrISTe2d+N32bGTkVgEAVs6Mwh/nxcCBm78SEZkMgxGRFcivbsHKd4+ioLoVSgc5/vfuONyVOFLqsoiIbA6DEZGF23+hGqkfHEdzuxrBXs54674pmDTSW+qyiIhsEoMRkQV771AR/vKfs9BoBaZEjMCmX03meCIiIjNiMCKyQBqtwF+/zsHWHy4BABZPDsX6xXFwclBIWxgRkY1jMCKyMC0qNX774Qn9IOvfzx2H39wyiusTERENAwYjIgtS1dSOB7YeQU55E5wc5Hj1nnjcMSlE6rKIiOwGgxGRhbhY1YIVWw6jtOEK/NyVePv+qUgM535nRETDicGIyAIcK6rHw+8cQUNbJ6L83PDOg0kI93WVuiwiIrvDYEQksd05lVj1wXGo1FrEh3ljy4qp8HV3krosIiK7xGBEJKGPj5RgzWenoBXA7JgAvPnLRLgqeVsSEUmF/wITSeT/HSjAX78+BwBYOjUML941kdt7EBFJjMHIQOnp6UhPT4dGo5G6FLJyQgi8ticPb2TkAQD+5+ZorJkXw+n4REQWQCaEEFIXYU2amprg5eWFxsZGeHp6Sl0OWRmtVuCFaxZu5BpFRETDw9Df3+wxIhomGq3A05+dwsdHLwMAnl84AfdPj5S2KCIi6mFIwaizsxMVFRVoa2uDv78/fHx8TFUXkU3RaAV+/+lJfHa8FAq5DK8smYTFk0dKXRYREf3EoEd6Njc3Y9OmTbj55pvh6emJyMhIjB8/Hv7+/oiIiMDKlStx5MgRc9RKZJU0WoHff3I1FL1xbyJDERGRhRpUMNqwYQMiIyOxdetWpKSkYMeOHcjOzsaFCxeQmZmJ5557Dmq1GnPmzMG8efOQl5dnrrqJrIJGK/DUx9n47EQpHOQy/GNZIm6fFCx1WURE1I9BDb5etmwZ/vznP2PChAkDHqdSqbB161YolUo89NBDQy7SknDwNRlKrdHiqU9O4ovsMl0o+mUi5k1kKCIikoKhv785K22QGIzIEFqtwOquMUW6UDQZ8yYGSV0WEZHdMvT3t9GryV25cgVtbW3674uKirBx40Z8++23xp6SyCYIIbD2P2f0Y4oYioiIrIfRwWjhwoV49913AQANDQ1ITk7Gq6++ikWLFmHTpk0mK5DImggh8PLOXPzfoWLIZMCGe+IZioiIrIjRwej48eOYOXMmAODTTz9FYGAgioqK8O677+KNN94wWYFE1iR970W89X0BAODFRXFYmBAqcUVERDQYRgejtrY2eHh4AAB27dqFxYsXQy6X44YbbkBRUZHJCiSyFlt/KMTfd10AAPz59vH4ZXK4xBUREdFgGR2MRo8ejR07dqCkpATffvst5syZAwCoqqrioGSyO/85WYZ1X+YAAJ5MGYNHZkZLXBERERnD6GC0du1arF69GpGRkUhOTsb06dMB6HqPEhMTTVYgkaX78WINnvo4GwDwwIxI/Hb2GGkLIiIiow1pun5FRQXKy8sRHx8PuVyXsQ4fPgxPT0/ExMSYrEhLwun6dK1z5U24Z3MmmlVqLIgLwj+WTYZczg1hiYgszbBsIhsUFISgoJ4zbpKSkoZySiKrUdpwBQ9sPYxmlRpJkT7YcE8CQxERkZUbUjDKyMhARkYGqqqqoNVqe7y3ZcuWIRVGZMka2zrxwJbDqGxSYUyAO96+fyqcHRVSl0VERENkdDBat24dnn/+eUydOhXBwcGQyfhfymQfOjVa/OaDY8irakGQpzPeeSgJXq6OUpdFREQmYHQw2rx5M7Zt24b77rvPlPUQWbznv8zBDxdr4apUYMsD0xDi7SJ1SUREZCJGz0rr6OjAjBkzTFkLkcV7N/MS3jtUBJkMeP3eRMSGcAA+EZEtMToYPfLII/jggw9MWQuRRTuQV61fq+gPc2NwW2ygxBUREZGpGf0orb29Hf/85z+xZ88eTJo0CY6OPcdYbNiwYcjFEVmKguoWpL5/HBqtwOLJofj1zVzAkYjIFhkdjE6dOoWEhAQAwJkzZ3q8x4HYZEua2zvxyDtH0dSuxuRwb7x0Vxz/jhMR2Sijg9HevXtNWQeRRRJC4PefnEJBTSuCvZzx1n2clk9EZMsGPcZo7dq1OHbsmDlqIbI4bx8owM6zFXBUyLDpV1Pg7+EkdUlERGRGgw5Gly9fxvz58zFy5Eg89thj+Oabb9DR0WGO2syqra0NERERWL16tdSlkIU6VFCL/915HgCw9s4JSAjzlrYgIiIyu0EHoy1btqCiogIffvghPDw88OSTT8LPzw9333033n33XdTV1ZmjTpN78cUXccMNN0hdBlmoyqZ2rPrgBDRagbsSQ/Gr5HCpSyIiomFg1HR9uVyOmTNn4m9/+xvOnz+PrKwsJCcn46233kJISAhuuukm/P3vf0dpaamp6zWJvLw85ObmYv78+VKXQhaoU6PFqg+Oo6ZFhZggDw62JiKyI4MORo8++ihaWlp6vDZ+/Hj84Q9/wA8//ICSkhKsWLECBw4cwIcffjjogvbv348777wTISEhkMlk2LFjR69j0tPTERkZCWdnZyQnJ+Pw4cODusbq1auxfv36QddG9uHvu87jyKV6eDg5YNOvpsBFycHWRET2YtDB6F//+hfa2tr6fd/f3x8PP/wwvvjiC6PG77S2tiI+Ph7p6el9vv/RRx8hLS0Nzz33HI4fP474+HjMnTsXVVVV+mMSEhIwceLEXl9lZWX44osvMHbsWIwdO3bQtZHtO5hXg7e+LwAA/G3JJET5uUlcERERDadBT9cXQpijDr358+cP+Ihrw4YNWLlyJR588EEAuj3bvv76a2zZsgVr1qwBAGRnZ/f784cOHcL27dvxySefoKWlBZ2dnfD09MTatWv7PF6lUkGlUum/b2pqMqJVZA3qWjuQ9nE2AGBZUjjmxwVLWxAREQ07o7cEkUJHRweOHTuGlJQU/WtyuRwpKSnIzMw06Bzr169HSUkJLl26hL///e9YuXJlv6Go+3gvLy/9V1hY2JDbQZZHCIE//vsUqppVGOXvhmfvGC91SUREJIFBByOZTCbZQNSamhpoNBoEBvbcoyowMBAVFRVmuebTTz+NxsZG/VdJSYlZrkPSej+rGLtzKqFUyPH6vYlwVRq99ikREVkxox6lPfDAA3ByGnihu88++8zooobLAw88cN1jnJycrttWsm55lc144auuzWHnjcPEUC+JKyIiIqkMOhitWLHCHHUYxM/PDwqFApWVlT1er6ysRFBQkERVkTXrUGvx+IcnoFJrcdNYfzx0Y5TUJRERkYQGHYy2bt1qjjoMolQqMWXKFGRkZGDRokUAAK1Wi4yMDKxatcqs105PT0d6ejo0Go1Zr0PD6x/f5SG3ohk+bkr8/ReTIJdzvSIiIntmcQMpWlpacPHiRf33hYWFyM7Oho+PD8LDw5GWloYVK1Zg6tSpSEpKwsaNG9Ha2qqfpWYuqampSE1NRVNTE7y8+KjFFpwpbUT6vnwAwF8XTUSAh7PEFRERkdQGFYyKi4sRHm741gilpaUIDQ0dVEFHjx7FrFmz9N+npaUB0D3C27ZtG5YuXYrq6mqsXbsWFRUVSEhIwM6dO3sNyCYaSIdai9WfnIRGK3B7XDAWcGo+EREBkIlBLEwUGBiIRYsW4ZFHHsG0adP6PKaxsREff/wxXn/9dTz66KN44oknTFasJejuMWpsbISnp6fU5ZCRXtt9Aa9n5MHHTYldv7sJfu4cYE9EZMsM/f09qB6jnJwcvPjii7jtttvg7OyMKVOmICQkBM7Ozqivr0dOTg7Onj2LyZMn429/+xsWLFgw5IYQmdrZskak79U9rn1+4QSGIiIi0htUj1G3K1eu4Ouvv8bBgwdRVFSEK1euwM/PD4mJiZg7dy4mTpxojlotAnuMrFunRouF//gBOeVNmD8xCP/f8sncIJaIyA6Ypceom4uLC5YsWYIlS5YYXaC14aw02/DP/QXIKW/CCFdHPL9wIkMRERH1YFSPkT1jj5H1KqlrQ8qG76FSa7FxaQIWJQ5uYgAREVkvQ39/W9VeaUTGEkLguf+chUqtxYxRvliYECJ1SUREZIEYjMgu7M6pxHe5VXBUyPgIjYiI+sVgRDavrUONdV/q9kJ79KZojA5wl7giIiKyVEYHoxUrVmD//v2mrMWipaenIzY2tt/1m8hyvfndRZQ2XEGotwtWzRojdTlERGTBjA5GjY2NSElJwZgxY/DSSy+htLTUlHVZnNTUVOTk5ODIkSNSl0KDkFfZjLf3FwAA1v18AlyUCokrIiIiS2Z0MNqxYwdKS0vx2GOP4aOPPkJkZCTmz5+PTz/9FJ2dnaaskcgoQgg8+8UZqLUCKeMDkRLLbWOIiGhgQxpj5O/vj7S0NJw8eRJZWVkYPXo07rvvPoSEhOB3v/sd8vLyTFUn0aB9e7YChwrq4Owox3N3xkpdDhERWQGTDL4uLy/H7t27sXv3bigUCixYsACnT59GbGwsXnvtNVNcgmhQOtRavPxNLgDg0ZnRCPNxlbgiIiKyBkYHo87OTvz73//GHXfcgYiICHzyySd48sknUVZWhnfeeQd79uzBxx9/jOeff96U9RIZ5P2sIlyqbYOfuxMevXmU1OUQEZGVMGpLEAAIDg6GVqvFsmXLcPjwYSQkJPQ6ZtasWfD29h5CeUSD19jWidczdI9x024bC3cno/+aExGRnTH6N8Zrr72GX/ziF3B2du73GG9vbxQWFhp7CYvCvdKsR/q+i2ho68TYQHfcM3Wk1OUQEZEVMfpR2s033wwnJ6derwshUFxcPKSiLBGn61uHkro2bPvhEgDg6QXj4aDgGqZERGQ4o39rREVFobq6utfrdXV1iIqKGlJRRMb627fn0aHR4mej/XDLWH+pyyEiIitjdDASQvS531RLS8uAj9eIzOVEcT2+PFkGmQz404Lx3A+NiIgGbdBjjNLS0gAAMpkMzz77LFxdr06D1mg0yMrK6nMgNpG5/X3XeQDA3ZNHIjbEU+JqiIjIGg06GJ04cQKArsfo9OnTUCqV+veUSiXi4+OxevVq01VIZIDDhXX44WItHBUyPJnC/dCIiMg4gw5Ge/fuBQA8+OCDeP311+Hpyf8yJ+m9tvsCAOAXU8MwcgQXcyQiIuMYPV1/69atpqyDyGiHCmqRWaDrLUqdNVrqcoiIyIoNKhilpaXhhRdegJubm36sUX82bNgwpMIsDdcxslzdvUVLp4Uh1NtF4mqIiMiaDSoYnThxAp2dnfr/3x9bnA2UmpqK1NRUNDU1wcvLS+pyqMuP+TXIKqyDUiFnbxEREQ3ZoIJR9/iin/5/IikIIbBxj27rj3uTwhDsxd4iIiIaGqPXMbpy5Qra2tr03xcVFWHjxo3YtWuXSQojup7M/Foc7uot+s0t7C0iIqKhMzoYLVy4EO+++y4AoKGhAUlJSXj11VexcOFCbNq0yWQF2hO1Rit1CVZDCIHX9ujGFi1LCkOQFxcVJSKioTM6GB0/fhwzZ84EAHz66acICgpCUVER3n33XbzxxhsmK9AeCCHwvztzccP671BY0yp1OVbhaFE9jlyqh9JBjsfYW0RERCZidDBqa2uDh4cHAGDXrl1YvHgx5HI5brjhBhQVFZmsQHsgk8lwrrwJNS0qbD9iexvwmsPb+wsAAHdPDmVvERERmYzRwWj06NHYsWMHSkpK8O2332LOnDkAgKqqKi76aIRlSeEAgE+PXkaHmo/UBlJQ3YLd5yoBAA//LFriaoiIyJYYHYzWrl2L1atXIzIyEsnJyZg+fToAXe9RYmKiyQq0F7fGBCDAwwm1rR3YnVMpdTkW7V8HCyEEMDsmAKMD3KUuh4iIbIjRwWjJkiUoLi7G0aNHsXPnTv3rs2fPxmuvvWaS4uyJo0KOpdPCAAAfHubjtP7Utqjw6bHLAICVN7G3iIiITMvoYAQAQUFBSExMhFx+9TRJSUmIiYkZcmGWJj09HbGxsZg2bZrZrnHP1DDIZMDBizW4xEHYffq/Q8VQqbWIC/VCcpSP1OUQEZGNMXqvNADIyMhARkYGqqqqoNX2HBezZcuWIRVmaYZj5eswH1fcNMYf31+oxvYjJVgz3/YC5lC0d2rwbuYlALreIltcYZ2IiKRldI/RunXrMGfOHGRkZKCmpgb19fU9vsg4+kHYx0o4CPsnPj9RitrWDoR6u2DBxCCpyyEiIhtkdI/R5s2bsW3bNtx3332mrMfuzR4fAH8PJ1Q3q7DnXCUWxAVLXZJF0GoF3j6gm6L/4I2RcFAM6SkwERFRn4z+7dLR0YEZM2aYshaCbhD2PVNHAuAg7GvtPV+FgupWeDg56AepExERmZrRweiRRx7BBx98YMpaqMu903SP0w7k1aC4tu06R9uHdzJ1i4YuSw6Hh7OjxNUQEZGtMvpRWnt7O/75z39iz549mDRpEhwde/6y2rBhw5CLs1dhPq6YOcYPB/Jq8OGRYvxxnn0Pwi6pa8OBvGoAwK+SIySuhoiIbJnRwejUqVNISEgAAJw5c6bHe5wtNHTLk8NxIK8Gnx67jNVzxkEht98/04+OlEAIYOYYP4T7ukpdDhER2TCjg9HevXtNWQf9xK0xgfBycUR1swqHCmpx42g/qUuSRKdGi4+OlgC4OmOPiIjIXIY0tefAgQP41a9+hRkzZqC0tBQA8N577+HgwYMmKc6eKR3k+hlpO06USlyNdDLOVaG6WQU/dyVSxgdKXQ4REdk4o4PRv//9b8ydOxcuLi44fvw4VCoVAKCxsREvvfSSyQq0ZwsTQgAAO89UoL1TI3E10uiembdkShiUDpyiT0RE5mX0b5q//vWv2Lx5M95+++0eA69vvPFGHD9+3CTF2bukSB8EezmjWaXGvvNVUpcz7Erq2rC/a9D1vZyiT0REw8DoYHT+/HncdNNNvV738vJCQ0PDUGqiLnK5DD+P1/UafZFdJnE1w+/jo7pB1zeO9kWkn5vU5RARkR0wOhgFBQXh4sWLvV4/ePAgoqNtb9fz4dhEti8/73qclpFbhab2zmG9tpTUGi0+OqIbdP3LJE7RJyKi4WF0MFq5ciV++9vfIisrCzKZDGVlZXj//fexevVqPPbYY6as0SKkpqYiJycHR44cGdbrxgZ7YnSAOzrUWuw8UzGs15bSd7lVqGpWwddNidtiOeiaiIiGh9HT9desWQOtVovZs2ejra0NN910E5ycnLB69Wo8/vjjpqzRrslkMixKCMHfd13Af7LLcM9U+xhrox90PXUkB10TEdGwMfo3jkwmwzPPPIO6ujqcOXMGhw4dQnV1NV544QVT1kcAfh4fCgD4Mb8GVU3tEldjfhWN7dh3oXvQNdcuIiKi4WN0MCouLoYQAkqlErGxsUhKSoK7u7v+PTKdcF9XJIZ7QyuAL0+VS12O2X11qgxCAFMjRiCKg66JiGgYGR2MoqKiUF1d3ev12tpaREVFDako6m1Rgq7X6D/Ztr/Y439O6mbgdQ88JyIiGi5GByMhRJ97orW0tMDZ2XlIRVFvC+KCoZDLcPJyIwprWqUux2wKa1px6nIjFHKZfuVvIiKi4TLowddpaWkAdGOMnn32Wbi6Xt3UU6PRICsrS7+5LJmOv4cTZozyxYG8Gvz3dDlSZ42WuiSz+LKrt+jG0X7wc3eSuBoiIrI3gw5GJ06cAKDrMTp9+jSUSqX+PaVSifj4eKxevdp0FZLevIlBOJBXgz3nKm0yGAkh8EXXo8LuhS2JiIiG06CD0d69ewEADz74IN544w14eHiYvCjqW8r4QDzz+RlklzSgqrkdAR629cgyp7wJ+dWtUDrIMXcC1y4iIqLhZ/Q6Rlu3bkVGRgYyMjJQVVUFrVbb4/0tW7YMuTjqKdDTGfEjvXDyciO+O1eFe5Nsayp796DrW8cFwMPZ8TpHExERmZ7Rg6+ff/55zJkzBxkZGaipqUF9fX2PLzKPlPG6npQ95yolrsS0tFqBr07qliLgbDQiIpKK0T1GmzZtwrZt23DfffeZsh66jpTYQLy6+wIO5NWgrUMNV6XRH6FFOV5cj9KGK3B3csCtMQFSl0NERHbK6B6jjo4OzJgxw5S1kAFigjwwcoQLVGotDubVSF2OyXQ/RpsTGwhnR4XE1RARkb0yOhg98sgj+OCDD0xZCxlAJpPpH6ftzrGNx2lqjRb/Pa17jHYnH6MREZGEjH4O097ejn/+85/Ys2cPJk2aBEfHnoNlN2zYMOTiqG9zYgOx7cdL+C63ChqtgELee6FNa/Jjfi1qWjowwtURPxvtJ3U5RERkx4wORqdOndIv5HjmzJke7/W1IjaZzrQoH3g6O6C2tQMniusxNdJH6pKG5Ouu/d8WxAXDUWF0JyYREdGQGR2Mutczshfp6elIT0+HRqORuhQ4KuSYFROAL7LLsPtcpVUHI61W4LvzVQCA+RO5BQgREUmL/3luoNTUVOTk5ODIkSNSlwLgmmn7Vj7O6HRpI6qbVXBTKpAUZb0Bj4iIbMOge4wWL15s0HGfffbZoIshw908zh+OChnyq1tRUN2CaH93qUsySkaurrfoprH+UDowpxMRkbQGHYy8vLzMUQcNkqezI26I9tXvnfaolQaj73J1PV5cu4iIiCzBoIPR1q1bzVEHGSFlfCAO5NVgd04lHr1plNTlDFpFYzvOlDZBJgNuGcdgRERE0uOzCyvW3ctyvLgBze2dElczeHu7Bl3Hj/SGv4eTxNUQERExGFm1MB9XRPi6QqMVOFxYJ3U5g5ZxTheMUsazt4iIiCwDg5GVmzFKtyDiDxdrJa5kcNo7Nfjhom5Lk1tjAiWuhoiISIfByMrdONoXAPBjvnXtm5ZZUIsrnRoEezljfLCH1OUQEREBYDCyetOjdcEot6IZNS0qiasxXMa5q7PRuFI6ERFZCgYjK+fr7oSYIF2PS2a+dTxOE0Lgu67xRbM5voiIiCwIg5ENuLFr41VreZyWW9GMssZ2ODvK9WOkiIiILAGDkQ3oHmdkLQOwv+ta7fpno/3g7KiQuBoiIqKrGIxsQFKULxzkMhTXtaGkrk3qcq7r6vgizkYjIiLLwmBkA9ydHBAf5g3A8h+n1bd24ERJAwBuA0JERJaHwchG3DjKOh6nZRXWQghgbKA7grycpS6HiIioBwYjGzFDPwC7FkIIiavp36EC3QrdN3QtM0BERGRJGIxsRGK4N5wd5ahpUSGvqkXqcvp1qEDXo8VgRERElojByEY4OSgwLdIHAPRbbVia+tYO5FY0AwCSonwkroaIiKg3BiMbYun7pmV1bXQ7NtAdfu5OEldDRETUG4ORDelezyiroBZqjVbianrrfoyWHMXHaEREZJkYjGzIhBAveDo7oFmlxunSRqnL6YXji4iIyNIxGNkQhVymDx2Hux5bWYprxxclR3N8ERERWSYGIxszJWIEACC7axFFS9E9vmhMAMcXERGR5WIwsjGJ4bpgdKK4QdpCfiKrsGt8EXuLiIjIgjEY2Zi4UC8o5DJUNLWjvPGK1OXocWFHIiKyBg5SFyCFyMhIeHp6Qi6XY8SIEdi7d6/UJZmMi1KBmCAPnC1rwoniBgTHuUhdEhraOpBb0QSAM9KIiMiy2WUwAoAff/wR7u7uUpdhFonh3l3BqB4L4oKlLgdZhXUQAhgd4A5/D44vIiIiy8VHaTYoMcyyBmBfXb+I44uIiMiyWVww2r9/P+68806EhIRAJpNhx44dvY5JT09HZGQknJ2dkZycjMOHDw/qGjKZDDfffDOmTZuG999/30SVW46EcG8AwKnLjei0gIUeOb6IiIishcU9SmttbUV8fDweeughLF68uNf7H330EdLS0rB582YkJydj48aNmDt3Ls6fP4+AgAAAQEJCAtRqda+f3bVrF0JCQnDw4EGEhoaivLwcKSkpiIuLw6RJk/qsR6VSQaVS6b9vamoyUUvNJ8rXDV4ujmi80onc8mbEjfSSrJYe44s4I42IiCycxQWj+fPnY/78+f2+v2HDBqxcuRIPPvggAGDz5s34+uuvsWXLFqxZswYAkJ2dPeA1QkNDAQDBwcFYsGABjh8/3m8wWr9+PdatW2dES6Qjl8uQEOaN7y9U40RJvaTB6HDX+KJR/m4I8HCWrA4iIiJDWNyjtIF0dHTg2LFjSElJ0b8ml8uRkpKCzMxMg87R2tqK5mbdCswtLS347rvvMGHChH6Pf/rpp9HY2Kj/KikpGVojhkli1+O0bInXM+JjNCIisiYW12M0kJqaGmg0GgQGBvZ4PTAwELm5uQado7KyEnfddRcAQKPRYOXKlZg2bVq/xzs5OcHJyfpmUiWEeQMATkg8APvkZd31u1fkJiIismRWFYxMITo6GidPnpS6DLPrDkaFNa2ob+3ACDflsNeg1mhxtky3me2kkd7Dfn0iIqLBsqpHaX5+flAoFKisrOzxemVlJYKCgsx67fT0dMTGxg7Yu2RJvF2ViPZ3AyDdtP28qha0d2rh7uSAaD83SWogIiIaDKsKRkqlElOmTEFGRob+Na1Wi4yMDEyfPt2s105NTUVOTg6OHDli1uuYUvd6RlI9Tjt9WddbNDHUE3K5TJIaiIiIBsPiHqW1tLTg4sWL+u8LCwuRnZ0NHx8fhIeHIy0tDStWrMDUqVORlJSEjRs3orW1VT9Lja5KCPfGv49fxoniekmuf6q0AQAfoxERkfWwuGB09OhRzJo1S/99WloaAGDFihXYtm0bli5diurqaqxduxYVFRVISEjAzp07ew3IJiCxa5xRdkkDtFox7L02p7p6jOJCpVsugIiIaDAsLhjdcsstEEIMeMyqVauwatWqYarIesUEecDZUY7mdjUKalowOsBj2K6tUmtwrly3sOMkCddRIiIiGgyrGmNEg+OgkOsfY50Y5vWMLlS0oFMj4OXiiHAf12G9NhERkbEYjAxkbbPSuiVKtJ7R1fFFXpDJOPCaiIisA4ORgaxxVhpwdQXs4e4xOlXC8UVERGR9GIxsXPejtLzKZqjUmmG77qlSLuxIRETWh8HIxgV7OcPT2QFqrUB+VeuwXLO9U4MLlbr96DjwmoiIrAmDkY2TyWSICfYEAORWNA3LNc+WNUGjFfBzVyLYy3lYrklERGQKDEZ2ICZIN03/fEXzsFzvdNfGsZNGenPgNRERWRUGIwNZ66w0AIgJ0vUYnRumYNQ9vogDr4mIyNowGBnIWmelAUBMsK7HKLd8eB6lda94zfFFRERkbRiM7MC4QF0wqmpWoa61w6zXalGpkV/dAgCIYzAiIiIrw2BkB9ycHBDhq1t92twDsM+WNkII3Wy4AA8OvCYiIuvCYGQnunuNcsvNO87oNMcXERGRFWMwshPDNWX/ZNf4oviurUiIiIisCYORgax5VhoAjO+asp9r5plp3VP12WNERETWiMHIQNY8Kw242mN0obIZGq0wyzUar3TiUm0bAAYjIiKyTgxGdiLcxxXOjnK0d2pRVGuerUG6F5AM9XbBCDelWa5BRERkTgxGdkIhl10dgG2mx2nd+6ONCXQ3y/mJiIjMjcHIjnSvgG2uhR7zuoLR2K4ARkREZG0YjOyIfgVss/UY6RZ2HBPAHiMiIrJODEZ2ZJyZZ6blVXU/SmOPERERWScGIzvS/SituK4NLSq1Sc9d19qBmhbddiPsMSIiImvFYGQga1/HCAB83JQI9HQCcHUGmal0jy8K9XaBm5ODSc9NREQ0XBiMDGTt6xh16+41MnUwulClG180ljPSiIjIijEY2ZkY/Tgj085M44w0IiKyBQxGdkY/M83Em8leXcOIwYiIiKwXg5Gd6X6Udq6iCUKYbmuQvEo+SiMiIuvHYGRnRvm7w0EuQ3O7GmWN7SY5Z22LCrWtuhlpozkjjYiIrBiDkZ1ROsgxyl8XXs6baJxR98KOYT4ucFVyRhoREVkvBiM71L3Q4/mKFpOcr3thx7EBHF9ERETWjcHIDnX3GBVUmyYYceA1ERHZCgYjOxTt7wYAKKhpNcn5LnDgNRER2QgGIwPZwsrX3bqDUb4JeoyEEFzDiIiIbAaDkYFsZeVrAIj20/XsNLR1oq5rNpmxals7UN/WCZns6iM6IiIia8VgZIdclAqEersAGPo4o+7xReE+rnBRKoZcGxERkZQYjOyUfpxR9dDGGXUv7DiGM9KIiMgGMBjZqe7HXkMdZ3RBP76Ij9GIiMj6MRjZqasDsE3TY8SB10REZAsYjOxU9wDsghrje4yEELhQ1b2GEXuMiIjI+jEY2anuHqPi2jZ0arRGnaO6RYWGtk7IOSONiIhsBIORnQrydIarUgG1VqC4rs2oc3Q/RovwdYOzI2ekERGR9WMwslNyuQxRfkObmabfCiSAvUVERGQbGIzsWPQQZ6blVXVN1ef4IiIishEMRnYsWt9jZFwwuljFNYyIiMi2MBjZsVFdj8CMfZSW3xWMRvNRGhER2QgGIwPZ0iay3bp7jIx5lFbf2oHarn3Wume4ERERWTsGIwPZ0iay3boDTX1bJ+oHuZlsd5gK9XaBq9LB5LURERFJgcHIjrkqHRDi5Qxg8As9do8vGsXHaEREZEMYjOzc1Zlpgxtn1N1jNIqP0YiIyIYwGNm5q3umGddjxIHXRERkSxiM7Fy0kYs8dvcwcSsQIiKyJQxGdu7qlH3De4zaOzUoqddtI8IeIyIisiUMRnaue4xR0SA2ky2saYUQgJeLI3zdlOYsj4iIaFgxGNm5YE9nODvKodYKlBi4mey144tkMpk5yyMiIhpWDEZ2TreZ7OBWwOaMNCIislUMRqQPOIbOTOOMNCIislUMRqQfZ2Roj5F+cUfOSCMiIhvDYET6HiNDVr/WaAUKa3QBij1GRERkaxiMSN/zk1fVAiHEgMeW1l+BSq2F0kGOkSNch6M8IiKiYcNgRBgd4A4HuQwNbZ0obbgy4LHd45Ci/dygkHNGGhER2RYGI4KzowJjAz0AAGdKGwc8luOLiIjIljEYEQBg0kgvAMCpywMHI/1UfY4vIiIiG8RgRACAuK5gdNrgHiOuYURERLaHwchA6enpiI2NxbRp06QuxSwmhXoD0PUYDTQAu7vHiDPSiIjIFjEYGSg1NRU5OTk4cuSI1KWYxdggdygVcjRe6cTl+r4HYNe2qFDf1gmZDIj2YzAiIiLbw2BEAAAnBwVignUDsPsbZ5TftQBkqLcLXJSKYauNiIhouDAYkd7E0K4B2KUNfb7PGWlERGTrGIxIb1JXMDrdT48R90gjIiJbx2BEetfOTOtrALZ+qj57jIiIyEYxGJHe2EAPKB3kaG5Xo6i2rdf77DEiIiJbx2BEeo4KOWKDPQEAp36yntGVDo1+uxCuYURERLaKwYh6iNOPM2ro8fqunAoAQICHE3zdnYa7LCIiomHBYEQ9xPWxNYgQApv25QMA7rshQpK6iIiIhgODEfXQvWfa2bImaLW6Adj7LlQjt6IZbkoF7p8eKWF1RERE5sVgRD2M9neHs6McLSo1Cmt1Czp29xb9MjkcXq6OUpZHRERkVgxG1IODQo4JIVfXMzpWVIfDhXVwVMjw8M+iJa6OiIjIvBiMqJfuAdinLjdi074CAMBdiaEI8nKWsiwiIiKzc5C6ALI83cHo27MVKG24ApkMePSmURJXRUREZH7sMaJeugdgd69bNDc2iIs6EhGRXWAwol6i/d3hqlTov//1LewtIiIi+8BgRL0o5DJM7BqAPWOULxLCvKUtiIiIaJhwjBH1afkN4ahuUWHN/BipSyEiIho2DEbUp4UJoViYECp1GURERMOKj9KIiIiIujAYEREREXVhMCIiIiLqwmBERERE1IXBiIiIiKiLXQajwsJCzJo1C7GxsYiLi0Nra6vUJREREZEFsMvp+g888AD++te/YubMmairq4OTk5PUJREREZEFsLtgdPbsWTg6OmLmzJkAAB8fH4krIiIiIkthcY/S9u/fjzvvvBMhISGQyWTYsWNHr2PS09MRGRkJZ2dnJCcn4/DhwwafPy8vD+7u7rjzzjsxefJkvPTSSyasnoiIiKyZxfUYtba2Ij4+Hg899BAWL17c6/2PPvoIaWlp2Lx5M5KTk7Fx40bMnTsX58+fR0BAAAAgISEBarW618/u2rULarUaBw4cQHZ2NgICAjBv3jxMmzYNt912W5/1qFQqqFQq/fdNTU0maikRERFZGosLRvPnz8f8+fP7fX/Dhg1YuXIlHnzwQQDA5s2b8fXXX2PLli1Ys2YNACA7O7vfnw8NDcXUqVMRFhYGAFiwYAGys7P7DUbr16/HunXrjGwNERERWROLe5Q2kI6ODhw7dgwpKSn61+RyOVJSUpCZmWnQOaZNm4aqqirU19dDq9Vi//79GD9+fL/HP/3002hsbNR/lZSUDLkdREREZJksrsdoIDU1NdBoNAgMDOzxemBgIHJzcw06h4ODA1566SXcdNNNEEJgzpw5uOOOO/o93snJibPWiIiI7IRVBSNTud7jOiIiIrJPVhWM/Pz8oFAoUFlZ2eP1yspKBAUFDUsNQggAHIRNRERkTbp/b3f/Hu+PVQUjpVKJKVOmICMjA4sWLQIAaLVaZGRkYNWqVWa9dnp6OtLT09HR0QEA+sHbREREZD2am5vh5eXV7/sWF4xaWlpw8eJF/feFhYXIzs6Gj48PwsPDkZaWhhUrVmDq1KlISkrCxo0b0draqp+lZi6pqalITU2FVqtFWVkZPDw8IJPJTHb+pqYmhIWFoaSkBJ6eniY7ryWx9TayfdbP1tto6+0DbL+NbJ/xhBBobm5GSEjIgMdZXDA6evQoZs2apf8+LS0NALBixQps27YNS5cuRXV1NdauXYuKigokJCRg586dvQZkm4tcLsfIkSPNdn5PT0+b/Mt+LVtvI9tn/Wy9jbbePsD228j2GWegnqJuFheMbrnllus+/1u1apXZH50RERGR/bGqdYyIiIiIzInByEI4OTnhueees+k1k2y9jWyf9bP1Ntp6+wDbbyPbZ34ycb3nVkRERER2gj1GRERERF0YjIiIiIi6MBgRERERdWEwIiIiIurCYGQGGo0Gzz77LKKiouDi4oJRo0bhhRde6LE+kxACa9euRXBwMFxcXJCSkoK8vLzrnjs9PR2RkZFwdnZGcnIyDh8+bM6m9Ol67evs7MQf//hHxMXFwc3NDSEhIbj//vtRVlY24Hn/8pe/QCaT9fiKiYkZjib1Yshn+MADD/Sqd968edc9tzV8hgB6ta3765VXXun3vJb0GTY3N+PJJ59EREQEXFxcMGPGDBw5ckT/vjXfg90GaqMt3IfX+wyt+R4Ert8+a7sH9+/fjzvvvBMhISGQyWTYsWNHj/cNuefq6uqwfPlyeHp6wtvbGw8//DBaWloGvG57eztSU1Ph6+sLd3d33H333b32VB0UQSb34osvCl9fX/HVV1+JwsJC8cknnwh3d3fx+uuv6495+eWXhZeXl9ixY4c4efKk+PnPfy6ioqLElStX+j3v9u3bhVKpFFu2bBFnz54VK1euFN7e3qKysnI4mqV3vfY1NDSIlJQU8dFHH4nc3FyRmZkpkpKSxJQpUwY873PPPScmTJggysvL9V/V1dXD0aReDPkMV6xYIebNm9ej3rq6ugHPay2foRCiR7vKy8vFli1bhEwmE/n5+f2e15I+w3vuuUfExsaK77//XuTl5YnnnntOeHp6isuXLwshrPse7DZQG23hPrzeZ2jN96AQ12+ftd2D//3vf8UzzzwjPvvsMwFAfP755z3eN+SemzdvnoiPjxeHDh0SBw4cEKNHjxbLli0b8Lq//vWvRVhYmMjIyBBHjx4VN9xwg5gxY4bR7WAwMoPbb79dPPTQQz1eW7x4sVi+fLkQQgitViuCgoLEK6+8on+/oaFBODk5iQ8//LDf8yYlJYnU1FT99xqNRoSEhIj169ebuAUDu177+nL48GEBQBQVFfV7zHPPPSfi4+NNVeaQGNLGFStWiIULFw7qvNb8GS5cuFDceuutA57XUj7DtrY2oVAoxFdffdXj9cmTJ4tnnnnG6u9BIa7fxr5Y031oSPus+R405vOzpnvwp8HIkHsuJydHABBHjhzRH/PNN98ImUwmSktL+7xOQ0ODcHR0FJ988on+tXPnzgkAIjMz06ja+SjNDGbMmIGMjAxcuHABAHDy5EkcPHgQ8+fPB6DbGLeiogIpKSn6n/Hy8kJycjIyMzP7PGdHRweOHTvW42fkcjlSUlL6/RlzuV77+tLY2AiZTAZvb+8Bz52Xl4eQkBBER0dj+fLlKC4uNmXpBjO0jfv27UNAQADGjRuHxx57DLW1tf2e05o/w8rKSnz99dd4+OGHr3tuS/gM1Wo1NBoNnJ2de7zu4uKCgwcPWv09CFy/jX2xpvvQ0PZZ6z042M/P2u7BnzLknsvMzIS3tzemTp2qPyYlJQVyuRxZWVl9nvfYsWPo7Ozscd6YmBiEh4cb/Zla3F5ptmDNmjVoampCTEwMFAoFNBoNXnzxRSxfvhwAUFFRAQC9Nr4NDAzUv/dTNTU10Gg0ff5Mbm6uGVrRv+u176fa29vxxz/+EcuWLRtwU8Dk5GRs27YN48aNQ3l5OdatW4eZM2fizJkz8PDwMFdz+mRIG+fNm4fFixcjKioK+fn5+NOf/oT58+cjMzMTCoWi1zmt+TN855134OHhgcWLFw94Xkv5DD08PDB9+nS88MILGD9+PAIDA/Hhhx8iMzMTo0ePtvp7ELh+G3/K2u5DQ9pnzffgYD8/a7sHf8qQe66iogIBAQE93ndwcICPj0+/92VFRQWUSmWvsD/QvXw9DEZm8PHHH+P999/HBx98gAkTJiA7OxtPPvkkQkJCsGLFCqnLG7LBtK+zsxP33HMPhBDYtGnTgOe9trdi0qRJSE5ORkREBD7++GOD/ivJlAxp47333qs/Pi4uDpMmTcKoUaOwb98+zJ49e1jrHazB/h3dsmULli9f3uu/bn/Kkj7D9957Dw899BBCQ0OhUCgwefJkLFu2DMeOHRvWOszJ0DZa6314vfZZ8z0IDO7vqDXeg9aKj9LM4Pe//z3WrFmDe++9F3Fxcbjvvvvwu9/9DuvXrwcABAUFAUCvUfOVlZX6937Kz88PCoViUD9jLtdrX7fuf4yLioqwe/fuAf8rtS/e3t4YO3YsLl68aMryDWJoG68VHR0NPz+/fuu1xs8QAA4cOIDz58/jkUceGfR1pPwMR40ahe+//x4tLS0oKSnB4cOH0dnZiejoaKu/B7sN1MZu1nwfGtK+a1nTPQgY3j5rvQevZcg9FxQUhKqqqh7vq9Vq1NXV9fv5BAUFoaOjAw0NDf2ed7AYjMygra0NcnnPP1qFQgGtVgsAiIqKQlBQEDIyMvTvNzU1ISsrC9OnT+/znEqlElOmTOnxM1qtFhkZGf3+jLlcr33A1X+M8/LysGfPHvj6+g76Oi0tLcjPz0dwcPCQax4sQ9r4U5cvX0ZtbW2/9VrbZ9jtX//6F6ZMmYL4+PhBX0fKz7Cbm5sbgoODUV9fj2+//RYLFy60+nvwp/pqI2D992G3/tr3U9Z0D17reu2z9nsQMOz33vTp09HQ0NCjx+y7776DVqtFcnJyn+edMmUKHB0de5z3/PnzKC4uNv4zNWrINg1oxYoVIjQ0VD8V+rPPPhN+fn7iD3/4g/6Yl19+WXh7e4svvvhCnDp1SixcuLDXtMVbb71VvPnmm/rvt2/fLpycnMS2bdtETk6OePTRR4W3t7eoqKiwqPZ1dHSIn//852LkyJEiOzu7x7RRlUrVb/ueeuopsW/fPlFYWCh++OEHkZKSIvz8/ERVVdWwts+QNjY3N4vVq1eLzMxMUVhYKPbs2SMmT54sxowZI9rb2/tto7V8ht0aGxuFq6ur2LRpU5/nseTPcOfOneKbb74RBQUFYteuXSI+Pl4kJyeLjo4OIYR134PdBmqjLdyHA7XP2u/B67WvmzXdg83NzeLEiRPixIkTAoDYsGGDOHHihH4WpCH33Lx580RiYqLIysoSBw8eFGPGjOkxXf/y5cti3LhxIisrS//ar3/9axEeHi6+++47cfToUTF9+nQxffp0o9vBYGQGTU1N4re//a0IDw8Xzs7OIjo6WjzzzDM9/jHSarXi2WefFYGBgcLJyUnMnj1bnD9/vsd5IiIixHPPPdfjtTfffFOEh4cLpVIpkpKSxKFDh4ajST1cr32FhYUCQJ9fe/fu1Z/np+1bunSpCA4OFkqlUoSGhoqlS5eKixcvDnPrdK7Xxra2NjFnzhzh7+8vHB0dRUREhFi5cmWvf1yt9TPs9tZbbwkXFxfR0NDQ53ks+TP86KOPRHR0tFAqlSIoKEikpqb2aIc134PdBmqjLdyHA7XP2u9BIa7/d1QI67oH9+7d2+fftxUrVgghDLvnamtrxbJly4S7u7vw9PQUDz74oGhubta/3/33+tq/w1euXBG/+c1vxIgRI4Srq6u46667RHl5udHtkAlxzVK3RERERHaMY4yIiIiIujAYEREREXVhMCIiIiLqwmBERERE1IXBiIiIiKgLgxERERFRFwYjIiIioi4MRkRERERdGIyIyGrs27cPMpms14aRwyUjIwPjx4+HRqMxyfnWrFmDxx9/3CTnIiLT4MrXRGSRbrnlFiQkJGDjxo361zo6OlBXV4fAwEDIZLJhr2nKlClIS0vD8uXLTXK+mpoaREdHIzs7u98d44loeLHHiIishlKpRFBQkCSh6ODBg8jPz8fdd99tsnP6+flh7ty52LRpk8nOSURDw2BERBbngQcewPfff4/XX38dMpkMMpkMly5d6vUobdu2bfD29sZXX32FcePGwdXVFUuWLEFbWxveeecdREZGYsSIEXjiiSd6PP5SqVRYvXo1QkND4ebmhuTkZOzbt2/AmrZv347bbrsNzs7O+tf+8pe/ICEhAe+99x4iIyPh5eWFe++9F83NzfpjPv30U8TFxcHFxQW+vr5ISUlBa2ur/v0777wT27dvN80fHBENGYMREVmc119/HdOnT8fKlStRXl6O8vJyhIWF9XlsW1sb3njjDWzfvh07d+7Evn37cNddd+G///0v/vvf/+K9997DW2+9hU8//VT/M6tWrUJmZia2b9+OU6dO4Re/+AXmzZuHvLy8fms6cOAApk6d2uv1/Px87NixA1999RW++uorfP/993j55ZcBAOXl5Vi2bBkeeughnDt3Dvv27cPixYtx7QiGpKQkXL58GZcuXTLyT4uITMlB6gKIiH7Ky8sLSqUSrq6uCAoKGvDYzs5ObNq0CaNGjQIALFmyBO+99x4qKyvh7u6O2NhYzJo1C3v37sXSpUtRXFyMrVu3ori4GCEhIQCA1atXY+fOndi6dSteeumlPq9TVFSkP/5aWq0W27Ztg4eHBwDgvvvuQ0ZGBl588UWUl5dDrVZj8eLFiIiIAADExcX1+PnucxYVFSEyMtLwPyQiMgsGIyKyaq6urvpQBACBgYGIjIyEu7t7j9eqqqoAAKdPn4ZGo8HYsWN7nEelUsHX17ff61y5cqXHY7RukZGR+lAEAMHBwfprxcfHY/bs2YiLi8PcuXMxZ84cLFmyBCNGjNAf7+LiAkDX80VE0mMwIiKr5ujo2ON7mUzW52tarRYA0NLSAoVCgWPHjkGhUPQ47tow9VN+fn6or6836Prd11IoFNi9ezd+/PFH7Nq1C2+++SaeeeYZZGVlISoqCgBQV1cHAPD39zekuURkZhxjREQWSalUmmy9oGslJiZCo9GgqqoKo0eP7vE10GO7xMRE5OTkDPp6MpkMN954I9atW4cTJ05AqVTi888/179/5swZODo6YsKECUa1h4hMiz1GRGSRIiMjkZWVhUuXLsHd3R0+Pj4mOe/YsWOxfPly3H///Xj11VeRmJiI6upqZGRkYNKkSbj99tv7/Lm5c+finXfeGdS1srKykJGRgTlz5iAgIABZWVmorq7G+PHj9cccOHAAM2fO1D9SIyJpsceIiCzS6tWroVAoEBsbC39/fxQXF5vs3Fu3bsX999+Pp556CuPGjcOiRYtw5MgRhIeH9/szy5cvx9mzZ3H+/HmDr+Pp6Yn9+/djwYIFGDt2LP785z/j1Vdfxfz58/XHbN++HStXrhxSe4jIdLjyNRGRgX7/+9+jqakJb731lknO98033+Cpp57CqVOn4ODADnwiS8AeIyIiAz3zzDOIiIjQD64eqtbWVmzdupWhiMiCsMeIiIiIqAt7jIiIiIi6MBgRERERdWEwIiIiIurCYERERETUhcGIiIiIqAuDEREREVEXBiMiIiKiLgxGRERERF0YjIiIiIi6/P8YqDuwn3stZAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "t_m, i_m = nx.Mask(time_axis, intensity, 80, 100)\n", "\n", "plt.semilogy(t_m, i_m)\n", "plt.xlabel('time (ns)')\n", "plt.ylabel(r'Intensity ($\\Gamma$/ns)')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 3, "id": "19766082-5db5-4a82-92c2-d773fc4a5151", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "34065.619827558825\n" ] } ], "source": [ "count_rate = nx.CountRateEstimator(5e9, time_axis, intensity, nx.lib.moessbauer.Fe57, time_start=12, time_end=170, detector_efficiency=0.65, attenuation=0.25)\n", "print(count_rate)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.12.0" } }, "nbformat": 4, "nbformat_minor": 5 }