{ "cells": [ { "cell_type": "markdown", "id": "bb4539da-7b99-4fbe-8db9-a15334e8dbea", "metadata": {}, "source": [ "# Energy spectrum time gate" ] }, { "cell_type": "code", "execution_count": 1, "id": "0e08d6e9", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAGwCAYAAABSN5pGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABj3ElEQVR4nO3deXgUVdo28Lu7syckLIEsECBAWGQLi8REFJVoEFziOEx0UGJePlxRMQIjyOI2E0RBQJ1B5xUYfQdRZhwclYkyAdyIYQsq+2IwCHTCYhYSsnWf749OVVd1dUISutJJ9/27rlwk1dVVp7pi6vE5zznHIIQQICIiIiKZ0d0NICIiImprGCAREREROWCAREREROSAARIRERGRAwZIRERERA4YIBERERE5YIBERERE5MDH3Q1or6xWK06fPo0OHTrAYDC4uzlERETUBEIIlJeXIzo6GkZjw3kiBkgtdPr0acTExLi7GURERNQCJ0+eRI8ePRp8nQFSC3Xo0AGA7QMODQ11c2uIiIioKcrKyhATEyM/xxvCAKmFpG610NBQBkhERETtzOXKY9pEkfabb76J3r17IyAgAAkJCdixY0ej+2/YsAEDBw5EQEAAhg4dik2bNqleF0Jg4cKFiIqKQmBgIJKTk3H06FHVPr1794bBYFB9LV682OXXRkRERO2P2wOkDz74AJmZmVi0aBH27NmD4cOHIyUlBcXFxU733759O+69915MmzYN+fn5SE1NRWpqKvbt2yfvs2TJEqxcuRKrVq1CXl4egoODkZKSgqqqKtWxXnjhBZw5c0b+evzxx3W9ViIiImofDEII4c4GJCQk4Oqrr8Ybb7wBwDY6LCYmBo8//jieeeYZzf5paWmoqKjAp59+Km+75pprEB8fj1WrVkEIgejoaDz99NOYNWsWAKC0tBQRERFYu3Yt7rnnHgC2DNLMmTMxc+bMJrWzuroa1dXV8s9SH2ZpaSm72IiIiNqJsrIyhIWFXfb57dYMUk1NDXbv3o3k5GR5m9FoRHJyMnJzc52+Jzc3V7U/AKSkpMj7FxQUwGw2q/YJCwtDQkKC5piLFy9Gly5dMGLECLzyyiuoq6trsK1ZWVkICwuTvziCjYiIyHO5tUj73LlzsFgsiIiIUG2PiIjAoUOHnL7HbDY73d9sNsuvS9sa2gcAnnjiCYwcORKdO3fG9u3bMXfuXJw5cwbLli1zet65c+ciMzNT/lnKIBEREZHn8dpRbMpgZ9iwYfDz88NDDz2ErKws+Pv7a/b39/d3up2IiIg8j1u72MLDw2EymVBUVKTaXlRUhMjISKfviYyMbHR/6d/mHBOw1ULV1dXhxIkTzb0MIiIi8jBuDZD8/PwwatQo5OTkyNusVitycnKQmJjo9D2JiYmq/QFg8+bN8v6xsbGIjIxU7VNWVoa8vLwGjwkAe/fuhdFoRLdu3a7kkoiIiMgDuL2LLTMzE+np6Rg9ejTGjBmD5cuXo6KiAhkZGQCAqVOnonv37sjKygIAPPnkkxg3bhyWLl2KSZMmYf369di1axfefvttALaJn2bOnImXXnoJcXFxiI2NxYIFCxAdHY3U1FQAtkLvvLw83HjjjejQoQNyc3Px1FNP4b777kOnTp3c8jkQERFR2+H2ACktLQ1nz57FwoULYTabER8fj+zsbLnIurCwULWYXFJSEtatW4f58+dj3rx5iIuLw8aNGzFkyBB5nzlz5qCiogIPPvggSkpKMHbsWGRnZyMgIACArZ5o/fr1eO6551BdXY3Y2Fg89dRTqrokIiIi8l5unwepvWrqPApERETUdrSLeZCIiIiI2iIGSETUoFqLFbUWq7ubQa3kUo3F3U0gajMYIBGRUxarwNiXtyBp8RZYrOyJ93SbfjyDQQuzsfqbAnc3hahNYIBERE6VVNagqKwaZ8urUVJZ4+7mkM6eeD8fAPDCpwfc3BKitoEBEhEREZEDBkhE5BQ71YjImzFAIqLLMhgM7m4C6Yy3mEiNARIRXRanSyMib8MAiYicUsZEDI+IyNswQCIip6yKCMnKDBIReRkGSETklDIoYnzk+QxgERKREgMkInJKOTckM0hE5G0YIBGRU1arsovNjQ0hInIDBkhE5JSqBokREhF5GQZIROQUu9iIyJsxQCIip9Sj2NzYECIiN2CAREROqWuQGCERkXdhgERETimzRpxJ2wtwlD+RCgMkInKKXWxE5M0YIBGRU5xJm4i8GQMkInLKanX+PRGRN2CAREROMYPkXViCRKTGAImInOJabETkzRggEZFTzCARkTdjgERETnEmbe9iYB8bkQoDJCJyiovVEpE3Y4BERE5xokgi8mYMkIjIKWW3moUpJCLyMgyQiMgpzqRNRN6MARIROcUuNiLyZgyQiMgpFmkTkTdjgERETnEeJO9i4FzaRCoMkIjIKc6DRETejAESETnFpUaIyJsxQCIip9Q1SIyQiMi7MEAiIqfUXWzuawe1Di41QqTGAImInGKRNhF5MwZIROSUugaJARIReRcGSETkFGfS9i7sYSNSY4BERE5ZrYrvmUEiIi/DAImInGIGiYi8GQMkInJKFSAxQiIiL8MAiYic4kzaROTNGCARkVPsYiMib8YAiYic4kzaROTNGCARkVPKrBHnQSIib8MAiYicYhebdzFwrREiFQZIROQUi7SJyJsxQCIip9Q1SG5sCLU6dqkSMUAiogZwLTbvouxg4+0mYoBERA1QdbExheRV2KVKxACJiBrAIm3vxftNxACJiBrAeZC8F+83EQMkImqAeh4k97WDWh8DJCIGSETUAAFmkLwVbzcRAyQiaoDyIWnhE9PjKe8wA2IiBkhE1ADlI5LPS8+nmtbBje0gaisYIBGRU8q5jzjM3/MJ1pwRqbSJAOnNN99E7969ERAQgISEBOzYsaPR/Tds2ICBAwciICAAQ4cOxaZNm1SvCyGwcOFCREVFITAwEMnJyTh69KjTY1VXVyM+Ph4GgwF79+511SURtXvsYvNivN1E7g+QPvjgA2RmZmLRokXYs2cPhg8fjpSUFBQXFzvdf/v27bj33nsxbdo05OfnIzU1Fampqdi3b5+8z5IlS7By5UqsWrUKeXl5CA4ORkpKCqqqqjTHmzNnDqKjo3W7PqL2SlmkzfjIu7AGiagNBEjLli3D9OnTkZGRgauuugqrVq1CUFAQVq9e7XT/FStWYMKECZg9ezYGDRqEF198ESNHjsQbb7wBwJY9Wr58OebPn48777wTw4YNw7vvvovTp09j48aNqmP95z//wRdffIFXX31V78skandUXS7uawa5Ae83kZsDpJqaGuzevRvJycnyNqPRiOTkZOTm5jp9T25urmp/AEhJSZH3LygogNlsVu0TFhaGhIQE1TGLioowffp0vPfeewgKCrpsW6urq1FWVqb6IvJkqockMwpehWvvEbk5QDp37hwsFgsiIiJU2yMiImA2m52+x2w2N7q/9G9j+wgh8MADD+Dhhx/G6NGjm9TWrKwshIWFyV8xMTFNeh9Re8UMkndRdam6sR1EbYXbu9jc4fXXX0d5eTnmzp3b5PfMnTsXpaWl8tfJkyd1bCGR+7EGybso7zFrkIjcHCCFh4fDZDKhqKhItb2oqAiRkZFO3xMZGdno/tK/je2zZcsW5Obmwt/fHz4+PujXrx8AYPTo0UhPT3d6Xn9/f4SGhqq+iDyaKoPEB6anEw3+QOSd3Bog+fn5YdSoUcjJyZG3Wa1W5OTkIDEx0el7EhMTVfsDwObNm+X9Y2NjERkZqdqnrKwMeXl58j4rV67E999/j71792Lv3r3yNAEffPAB/vjHP7r0GonaK04U6WXYpUqk4uPuBmRmZiI9PR2jR4/GmDFjsHz5clRUVCAjIwMAMHXqVHTv3h1ZWVkAgCeffBLjxo3D0qVLMWnSJKxfvx67du3C22+/DQAwGAyYOXMmXnrpJcTFxSE2NhYLFixAdHQ0UlNTAQA9e/ZUtSEkJAQA0LdvX/To0aOVrpyobVNODskHpudjlyqRmtsDpLS0NJw9exYLFy6E2WxGfHw8srOz5SLrwsJCGI32RFdSUhLWrVuH+fPnY968eYiLi8PGjRsxZMgQeZ85c+agoqICDz74IEpKSjB27FhkZ2cjICCg1a+PqL3i2lzexcouVSIVg+B4zhYpKytDWFgYSktLWY9EHumFTw5g9bcFAICHru+DuRMHublFpKc+cz+Tg6Rvn7kJ3TsGurdBRDpp6vPbK0exEdHlcdi3d1HXnPGOEzFAIiKn1IuX8oHp6bhYLZEaAyQiuiw+MInI2zBAIiKnlFkjxkeezTFDyKJ8IgZIRNQAzoPkPRzvL+83EQMkImqA4LBvr+F4d3m3iRggEVEDlN0szCh4NscuNRblEzFAIqIGcNi393C8vVbebiIGSETknLqLjTyZtguVd5yIARIRNYBdbN6CRdpEWgyQiMgpFml7L95tIgZIRNQAzqzsPbQ1SLzhRAyQiMgprsXmPRwzhIyPiBggEVEDmEHyHqxBItJigERETnGYv/fQzIPEnCERAyQico4TRXoPzSB/3m8iBkhE1ACOYvMa7GIj0mKAREROcbFaL+IYIDEgJmKARETOKeuO+Lj0bBzFRqTFAImInGIGyXtwHiQiLQZIROQUZ9L2HlyJjUiLARIROSUa/IE8jeM0DkwgETFAIqIGsAbJe2jvL+84EQMkInJK+YhkTYpn00wUydtNxACJiJwTnCjSe2iKtN3TDKK2hAESETmlLtImT6adSZt3nIgBEhE5pV6slg9MT6aZSds9zSBqUxggEZFTyqH9fGB6Nk4USaTFAImInFI9JPnA9Gjatdh4w4kYIBGRU+r4iA9MT8aJIom0GCARkVPqGiT3tYP0x4kiibQYIBFRA+xPSc6D5Nm0Rdq830QMkIjIKWaQvId2sVr3tIOoLWGAREROWbnUiNfQjmLjHSdigERETqmKtPm89GicB4lIiwESETmlfmjykenJNHeXt5uIARIROccMkvdw7FJjUT4RAyQiaoBgDZLX0K7F5pZmELUpDJCI6LJYtOvZWINEpMUAiYicUj40Oezbs2kniuQNJ2KAREROcbFa7+F4fxkQEzFAIqIGqCeK5BPTk2lvL+83EQMkInKKI5m8h3aiSDc1hKgNYYBERE5xqRHvwSJtIi0GSETklGoeJD4yPZp2LTbebyIGSETkHDNIXoNdbERaDJCIyCnVKDY+MD0au9iItJodIBUVFeH+++9HdHQ0fHx8YDKZVF9E5BlUNUh8ZHoVjlokAnya+4YHHngAhYWFWLBgAaKiomAwGPRoFxG5mfIRyXlxPBtrjoi0mh0gffPNN/j6668RHx+vQ3OIqK0Q6hQSeTAWaRNpNbuLLSYmhulXIi/AUWzeg4vVEmk1O0Bavnw5nnnmGZw4cUKH5hBRW2HlKDavoV2LzU0NIWpDmt3FlpaWhsrKSvTt2xdBQUHw9fVVvX7hwgWXNY6I3EjxlOTz0rNpMkhuaQVR29LsAGn58uU6NIOI2hpVFxtTCh6NNUhEWs0OkNLT0/VoBxG1MazR9iacCInIUbMDJACwWCzYuHEjDh48CAAYPHgw7rjjDs6DRORBOFGk99BOFMkbTtTsAOnYsWOYOHEiTp06hQEDBgAAsrKyEBMTg88++wx9+/Z1eSOJqPWpF6vlA9OTOc5zxdtN1IJRbE888QT69u2LkydPYs+ePdizZw8KCwsRGxuLJ554Qo82EpEbsIvNezgGwJwYlKgFAdKXX36JJUuWoHPnzvK2Ll26YPHixfjyyy9b1Ig333wTvXv3RkBAABISErBjx45G99+wYQMGDhyIgIAADB06FJs2bVK9LoTAwoULERUVhcDAQCQnJ+Po0aOqfe644w707NkTAQEBiIqKwv3334/Tp0+3qP1EnkhdpO22ZlAr0I5i4w0nanaA5O/vj/Lycs32ixcvws/Pr9kN+OCDD5CZmYlFixZhz549GD58OFJSUlBcXOx0/+3bt+Pee+/FtGnTkJ+fj9TUVKSmpmLfvn3yPkuWLMHKlSuxatUq5OXlITg4GCkpKaiqqpL3ufHGG/Hhhx/i8OHD+Oc//4njx4/jt7/9bbPbT+SphGqYPx+YnkxTg8TbTdT8AOm2227Dgw8+iLy8PAghIITAd999h4cffhh33HFHsxuwbNkyTJ8+HRkZGbjqqquwatUqBAUFYfXq1U73X7FiBSZMmIDZs2dj0KBBePHFFzFy5Ei88cYbAGx/1JcvX4758+fjzjvvxLBhw/Duu+/i9OnT2Lhxo3ycp556Ctdccw169eqFpKQkPPPMM/juu+9QW1vb7Gsg8kTqGiT3tYP05xgA83YTtSBAWrlyJfr27YvExEQEBAQgICAA1157Lfr164cVK1Y061g1NTXYvXs3kpOT7Q0yGpGcnIzc3Fyn78nNzVXtDwApKSny/gUFBTCbzap9wsLCkJCQ0OAxL1y4gL///e9ISkrSTHwpqa6uRllZmeqLyJNxFJsX0WSQeMOJmj2KrWPHjvj4449x9OhRHDp0CAAwaNAg9OvXr9knP3fuHCwWCyIiIlTbIyIi5GM7MpvNTvc3m83y69K2hvaR/OEPf8Abb7yByspKXHPNNfj0008bbGtWVhaef/75pl0YkQdgkbb34FpsRFrNziBJ4uLicPvtt+P2229vUXDUFsyePRv5+fn44osvYDKZMHXq1Ab/z2nu3LkoLS2Vv06ePNnKrSVqXZxJ23toa5B4v4malEHKzMzEiy++iODgYGRmZja677Jly5p88vDwcJhMJhQVFam2FxUVITIy0ul7IiMjG91f+reoqAhRUVGqfeLj4zXnDw8PR//+/TFo0CDExMTgu+++Q2Jioua8/v7+8Pf3b/K1EbV3qiJtPi89muPSIrzdRE0MkPLz8+Xi5fz8/Ab3MxgMzTq5n58fRo0ahZycHKSmpgIArFYrcnJyMGPGDKfvSUxMRE5ODmbOnClv27x5sxzUxMbGIjIyEjk5OXJAVFZWhry8PDzyyCMNtsVqtQKw1RoRkUMGiY9Mj+Z4dzkPElETA6StW7c6/d4VMjMzkZ6ejtGjR2PMmDFYvnw5KioqkJGRAQCYOnUqunfvjqysLADAk08+iXHjxmHp0qWYNGkS1q9fj127duHtt98GYAvSZs6ciZdeeglxcXGIjY3FggULEB0dLQdheXl52LlzJ8aOHYtOnTrh+PHjWLBggVx8TkRQPTWZQfJsjl1q7GIjauFabEplZWXYsmULBg4ciIEDBzb7/WlpaTh79iwWLlwIs9mM+Ph4ZGdny0XWhYWFMBrtpVJJSUlYt24d5s+fj3nz5iEuLg4bN27EkCFD5H3mzJmDiooKPPjggygpKcHYsWORnZ2NgIAAAEBQUBA++ugjLFq0CBUVFYiKisKECRMwf/58dqMR1RMNfE+eh/eXSMsgmvm/Cr/73e9w/fXXY8aMGbh06RKGDx+OEydOQAiB9evX4+6779arrW1KWVkZwsLCUFpaitDQUHc3h8jlbnhlK06crwQA9O0ajJynb3Bvg0g3Ww8VI2PtTvnnZycOwvTr+7ixRUT6aerzu9mj2L766itcd911AIB//etfEEKgpKQEK1euxEsvvdTyFhNRm2LlMH+v4Vhj5li0TeSNmh0glZaWyuuwZWdn4+6770ZQUBAmTZqkWe+MiNov1UOTz0uPphnm755mELUpzQ6QYmJikJubi4qKCmRnZ+OWW24BAPz6669yjQ8RtX+cKNJ7cC02Iq1mF2nPnDkTU6ZMQUhICHr16oUbbrgBgK3rbejQoa5uHxG5iXotNj4xPZlmJm2GxETND5AeffRRjBkzBidPnsTNN98sjzDr06cPa5CIPBTnxfFsmokieb+JWjbMf/To0Rg9ejQAwGKx4Mcff0RSUhI6derk0sYRkfuoZtJmRsGjcakRIq1m1yDNnDkT77zzDgBbcDRu3DiMHDkSMTEx2LZtm6vbR0Ruol6LzW3NoFbBDBKRo2YHSP/4xz8wfPhwAMAnn3yCgoICHDp0CE899RSeffZZlzeQiNxDXYPkvnaQ/jiKjUir2QHSuXPn5AVhN23ahMmTJ6N///74n//5H/z4448ubyARuQe71byHpkibt56o+QFSREQEDhw4AIvFguzsbNx8880AgMrKSphMJpc3kIjcw8pRbF7D8fZyokiiFhRpZ2Rk4He/+x2ioqJgMBiQnJwMwLYAbEvWYiOitonzIHkPx2wh7zdRCwKk5557DkOGDMHJkycxefJkeXFXk8mEZ555xuUNJCJ3UYxi4xPTo2nuL284UcuG+f/2t7/VbEtPT7/ixhBR26F8RrLLxbNp5kFyUzuI2pImBUgrV67Egw8+iICAAKxcubLRfZ944gmXNIyI3Es08D15PgbERE0MkF577TVMmTIFAQEBeO211xrcz2AwMEAi8hCqiSL5vPRoXIuNSKtJAVJBQYHT74nIc7HTxXuwSJtIq9nD/InIO3CiSO/BDBKRVrOLtIUQ+Mc//oGtW7eiuLgYVqtV9fpHH33kssYRkfuo12IjT8a12Ii0mh0gzZw5E2+99RZuvPFGREREwGAw6NEuInIzdQaJD0xPphnl75ZWELUtzQ6Q3nvvPXz00UeYOHGiHu0hojaCo9i8h2MAzICYqAU1SGFhYejTp48ebSGiNkT5kLRa+cD0ZKxBItJqdoD03HPP4fnnn8elS5f0aA8RtRHMIHkPx1FsjIeJWtDF9rvf/Q7vv/8+unXrht69e8PX11f1+p49e1zWOCJyHysXY/MamgwSbzhR8wOk9PR07N69G/fddx+LtIk8GOMj78Gl2Ii0mh0gffbZZ/j8888xduxYPdpDRG2EqouNT0yPxttLpNXsGqSYmBiEhobq0RYiakNURdp8gHo0x7XXuBYbUQsCpKVLl2LOnDk4ceKEDs0horZC3cXGB6YnYxcbkVazu9juu+8+VFZWom/fvggKCtIUaV+4cMFljSMi97Eyg+Q1NPMgMSAman6AtHz5ch2aQURtjWjwB/I0nAeJSKtFo9iIyPMpH5KsSfFs2hokNzWEqA1pdg3Snj178OOPP8o/f/zxx0hNTcW8efNQU1Pj0sYRkXtou1zIk2njX95xomYHSA899BCOHDkCAPjpp5+QlpaGoKAgbNiwAXPmzHF5A4mo9Tk+MJlB8myO95e3m6gFAdKRI0cQHx8PANiwYQPGjRuHdevWYe3atfjnP//p6vYRkRvwgendeL+JWhAgCSFgtVoBAP/9738xceJEALb5kc6dO+fa1hGRWzh7PnKySM+lCYjZxUbU/ABp9OjReOmll/Dee+/hyy+/xKRJkwAABQUFiIiIcHkDiaj1OetSY3zkubRdqu5pB1Fb0uwAafny5dizZw9mzJiBZ599Fv369QMA/OMf/0BSUpLLG0hErc9ZMMQ6JM/lGBDxVhO1YJj/sGHDVKPYJK+88gpMJpNLGkVEbQ+fmZ7LsUuNXWxELQiQJDU1NSguLpbrkSQ9e/a84kYRkXs5yxYxg+S5NLeWt5qo+QHSkSNHMG3aNGzfvl21XQgBg8EAi8XissYRkXs4i4UYH3kuxwJ8BsNELQiQMjIy4OPjg08//RRRUVEwGAx6tIuI3Ki1i7RfzzmKHScu4J30q+Hn0+zSSI+25VARVuQcw9LJw9GvW4gu59DUIOlyFqL2pdkB0t69e7F7924MHDhQj/YQURvgdJi/jo/NpZttk89uOVSMCUMidTtPe/Q/a3cBAB7+v934b+Y4Xc7BtdiItJr9v2pXXXUV5zsi8nDOR7Hpc66aOnsdoz+zRw06VnxRt2Nr50Eiomb/NXr55ZcxZ84cbNu2DefPn0dZWZnqi4jaP2eTQuo1UeT5imr5+w4BLR43QlfA8c6yBomoBV1sycnJAIDx48ertrNIm8hztGYG6Wy5PUDiY9k9NMEvbwRR8wOkrVu36tEOImpDnGYQdHponrtoD5AsnMJZIzTAB2VVdbqeQ7rdJqMBFqvgPEhEaEGANG6cPkWCRNR2OHs86tXtcrHannW2MkDSCPbXP0CS7q3JYIAFgkXaRLiCiSIrKytRWFiImpoa1fZhw4ZdcaOIyL1aMYGkCoosfDK7hfSpG40ALKxBIgJaECCdPXsWGRkZ+M9//uP0ddYgEbV/Uk2KwWAPlvR6aCqPywSSe0j3wFg/rx3jI6IWjGKbOXMmSkpKkJeXh8DAQGRnZ+Nvf/sb4uLi8O9//1uPNhJRK5MzCoqJYPV6aCrrjtjFptUqU/FKNUhSgNQa5yRq45qdQdqyZQs+/vhjjB49GkajEb169cLNN9+M0NBQZGVlYdKkSXq0k4hakZRRMAAwGmyZHb2G+SsPyyJt95AzSEZmkIgkzc4gVVRUoFu3bgCATp064ezZswCAoUOHYs+ePa5tHRG5hfSANBoM8nJCej0zlXVHrEFyD+UoNtvPvA9EzQ6QBgwYgMOHDwMAhg8fjrfeegunTp3CqlWrEBUV5fIGElHrkx+PBlsGCWidGiQ+mN3DKgfEtn95F4ha0MX25JNP4syZMwCARYsWYcKECfj73/8OPz8/rF271tXtIyI3kGqBbN1rBkDHod+qUWzWRnYk3UjzHtmLtBkiETU7QLrvvvvk70eNGoWff/4Zhw4dQs+ePREeHu7SxhGRexlggMFge1jql0Gyf88uNvfQdLG5sS1EbUWzuthqa2vRt29fHDx4UN4WFBSEkSNHMjgi8iBWxTB/aSAbR7F5LuEwzJ+3gaiZGSRfX19UVVXp1RYiaiOURdqO21xNmZniKDb3kGuQ6v+XmV1sRC0o0n7sscfw8ssvo65O36nvich9pMejAfZ5ePRan0s9USQfzI3RbaoF2JcaISKbZtcg7dy5Ezk5Ofjiiy8wdOhQBAcHq17/6KOPXNY4InIPZRebEfp2uyiPywCpcVYBmHSIYawOGUPeBqIWBEgdO3bE3XffrUdbiKiNkB6QBoNBzi7olb2wcBRbk1msQi6kdiW5S1Uu0maERNTsLrY1a9Y0+tUSb775Jnr37o2AgAAkJCRgx44dje6/YcMGDBw4EAEBARg6dCg2bdqkel0IgYULFyIqKgqBgYFITk7G0aNH5ddPnDiBadOmITY2FoGBgejbty8WLVqkWXiXyFvZi3b1L9wVnCiyyfTKsEn3QOpiszJQJWp+gHTTTTehpKREs72srAw33XRTsxvwwQcfIDMzE4sWLcKePXswfPhwpKSkoLi42On+27dvx7333otp06YhPz8fqampSE1Nxb59++R9lixZgpUrV2LVqlXIy8tDcHAwUlJS5ALzQ4cOwWq14q233sL+/fvx2muvYdWqVZg3b16z20/kieQaJIMB9rIUvTJI9u85iq1x+gVItn+ZQSKya3aAtG3bNqeZlqqqKnz99dfNbsCyZcswffp0ZGRk4KqrrsKqVasQFBSE1atXO91/xYoVmDBhAmbPno1BgwbhxRdfxMiRI/HGG28AsP2f0PLlyzF//nzceeedGDZsGN59912cPn0aGzduBABMmDABa9aswS233II+ffrgjjvuwKxZsxqtn6qurkZZWZnqi8hTyV1s0D+DxCLtptP7Hhh1ntKBqD1pcoD0ww8/4IcffgAAHDhwQP75hx9+QH5+Pt555x107969WSevqanB7t27kZycbG+Q0Yjk5GTk5uY6fU9ubq5qfwBISUmR9y8oKIDZbFbtExYWhoSEhAaPCQClpaXo3Llzg69nZWUhLCxM/oqJiWnSNRK1R/YibYN9FBuH+budXp+PdFROFElk1+Qi7fj4+Pp0u8FpV1pgYCBef/31Zp383LlzsFgsiIiIUG2PiIjAoUOHnL7HbDY73d9sNsuvS9sa2sfRsWPH8Prrr+PVV19tsK1z585FZmam/HNZWRmDJPJY9iJtAPIoNg7zdze9CuXtGSQuNUIkaXKAVFBQACEE+vTpgx07dqBr167ya35+fujWrRtMJpMujdTTqVOnMGHCBEyePBnTp09vcD9/f3/4+/u3YsuI3EfZ5SKVCOk3k7bz70lLtwyb41IjjI+Imh4g9erVCwBgdeHwhvDwcJhMJhQVFam2FxUVITIy0ul7IiMjG91f+reoqAhRUVGqfeLj41XvO336NG688UYkJSXh7bffvtLLIfI4Bhhg1HktNsEMUqPU80TpdQ6HGiR9TkPUrjS7SNuV/Pz8MGrUKOTk5MjbrFYrcnJykJiY6PQ9iYmJqv0BYPPmzfL+sbGxiIyMVO1TVlaGvLw81TFPnTqFG264AaNGjcKaNWtgNLr1oyBqU1RrsUHf2ZW5FlvjWqMLUjoqu9iI7Jo9UaSrZWZmIj09HaNHj8aYMWOwfPlyVFRUICMjAwAwdepUdO/eHVlZWQCAJ598EuPGjcPSpUsxadIkrF+/Hrt27ZIzQAaDATNnzsRLL72EuLg4xMbGYsGCBYiOjkZqaioAe3DUq1cvvPrqqzh79qzcnoYyV0TeRLkWmzTkW78aJPv3nAdJqzUCJKtDFxvjVKI2ECClpaXh7NmzWLhwIcxmM+Lj45GdnS0XWRcWFqqyO0lJSVi3bh3mz5+PefPmIS4uDhs3bsSQIUPkfebMmYOKigo8+OCDKCkpwdixY5GdnY2AgAAAtozTsWPHcOzYMfTo0UPVHv6fE5G6i8Wg8/ITqgCAT2YNVQCp1yg2aaJIjmIjkrk9QAKAGTNmYMaMGU5f27Ztm2bb5MmTMXny5AaPZzAY8MILL+CFF15w+voDDzyABx54oCVNJfIKck2K0R4YtcYoNmaQtJRBkV4fj3JpGV1PRNSOsPCGiDTsE0XaZ9LW65GpqkHic1mjNeaJkrpRTSzSJpK5LEBKT09v0VIjRNT2OFuLTb85eBTfM0LSsKoCSJ3uQf3gZHsNEu8Dkcu62Lp3786RYEQeQrkWm5RO0q0Gyap/hqQ9Uw/z1zeDZNS53oyoPXFZgPSnP/3JVYciIjdTrsVm0HktNgtrkBql+nx0mkhTurcMkIjsmPIhIg3lPEiQ63b1L9Lmg1mrNSbSlA7LUWxEdk3KICnXILucZcuWtbgxRNQ2KEc1SdNE6jaLM7vYGmVphc9HrjkzcqJIIkmTAqT8/PwmHUweIkpE7ZqySFt6Vgqd8gqcKLJxys9Ht2H+9f+aOMqfSNakAGnr1q16t4OI2hC5SFs5zL8VapA4ik3N8fPQK4C0r8UmdbHxPhC1uAbp2LFj+Pzzz3Hp0iUATMkSeRJ7Fxt0D5CEYBdbQxxrjvSuQbJ3selyGqJ2pdkB0vnz5zF+/Hj0798fEydOxJkzZwAA06ZNw9NPP+3yBhJR67MXaRsUo9j0eWpyosiGOWaM9MqwSffWpPO9JmpPmh0gPfXUU/D19UVhYSGCgoLk7WlpacjOznZp44jIPexdbPIgNt06XVpjnp/2yvHj0DuANHIUG5Gs2fMgffHFF/j88881i7zGxcXh559/dlnDiMh9lGuxOW5z+bk4iq1Bjp+5Xp+PVVGUD4AREhFaECBVVFSoMkeSCxcuwN/f3yWNIiI3U6zF5rjN1bhYbcMcAyK9aj05DxKRVrO72K677jq8++678s8GgwFWqxVLlizBjTfe6NLGEZF72JeesGcVdKtB4lpsDXL8OFprFBu7OolakEFasmQJxo8fj127dqGmpgZz5szB/v37ceHCBXz77bd6tJGIWpm0eCkUc5u1xig2PpjVHANGveJHTQaJt4Go+RmkIUOG4MiRIxg7dizuvPNOVFRU4De/+Q3y8/PRt29fPdpIRK1MWaStewZJVYOkyynaLc0wf91m0rb9K91rzoNE1MwMUm1tLSZMmIBVq1bh2Wef1atNRORmmqJd6DmKjRmkhmiG+etVgyR1qTKDRCRrVoDk6+uLH374Qa+2EFEb4WwtNt0Wq1VkjTiKTc3xI9dvFJvtX2keJAZIRC3oYrvvvvvwzjvv6NEWImoz7BkkvWfSZgapYY4BkX41SPUTRXKxWiJZs4u06+rqsHr1avz3v//FqFGjEBwcrHp92bJlLmscEbmH9CA2wCDPFKnXw9nCAKlBrbXUiFWRMQQ4zJ8IaEGAtG/fPowcORIAcOTIEdVrBsWIFyJqv5RrscnbdHpsKgMvdrGpWR2K1vWrQbJhFxuRXbMDpK1bt+rRDiJqQ+xrsdkni9QrdlGOzHIMCLxda82kbe9iq/+ZOSSi5tcgEZHnsw/zNyhqkPSdpBDgTNqOHD8P/eaisv0rjWJjIo+oBRkkIvJ8UjBkNNozSHo9nJVZEdYgqTkGpfqvxcYuNiIJM0hEpCHXICkzSDp1uygfxlxqRM1x4kzdapAchvmzTJuIARIROSEFQ8oibb3qgyzsYmtQ641iU99v3gYidrERkRNSMKSaKFKvcykDJBZpq7TWPEgSk5GL1RJJmEEiIg3lWmwGnddiU3arcYJCtdabSdthokhdzkLUvjCDREQayrXYDPYiJJ3OZf+e8yCpaUex6VuDZGCRNpGMGSQi0nKyFpteGSRlUMQaJLXWmgdJziAZuNQIkYQZJCLSEIq12KDzRJHKhzFHsak5fh663YP6f+WJInkbiJhBIiIt+4PYUB8k6ZhB4ii2BjkGRHrXgZmMtkcC7wIRM0hE5IS9aNc+eaDeC6UCXGrEkXYUm75Bqq+Jo9iIJMwgEZGGlFEwGgzy8hO61b9wJu0GaWfS1uc8UmDqW9/HxmJ5IgZIROSE9IA0Gg1y4a5eD01VFxsfzCqOXY56F8pLARIDVSIGSETkhBSnmAwG+9w4eg3zZwapQZoaJJ1HsfmYuFgtkYQBEhFpqOdBsm3Tq4Ca8yA1zDEg0u8e1GeQjOxiI5IwQCIijdbsYlNmjfhcVtOuxabPeaR7K2WQAE65QMQAiYg0pGejUdHFptcDUxl48aGsphnFpvM9kGqQAHZ3EjFAIiIN5czKRqPeE0Xav+c8SGqtNQ+SdFg/RYDEe0HejgESEWlYnXWxtcZEkcwgqWiWGtH5Hii72BgfkbdjgEREGhZFkbY8k3Yr1CDxoazmGCDp9fnYu9gMmm1E3ooBEhFp2JeeUEwUqcPTWQjBLrZGOMYoehfK+xjZxUYkYYBERBqqIm2DfkXarRUAtFfaxWp1LtL2sT8SBJd9IS/HAImINOxdbIpRbDo8nJ0FRBzJZqcZ5q9bBsn2r6qLjRkk8nIMkIhIw97FBsVabDqcx8lDmA9mO+1ita4/hzLoUnWxMVAlL8cAiYg0rIoMklykrUPg4uyYnH/HzvGj0CN4VB5TWZTvuFAukbdhgEREGlK2SO+ZtJ0d0sraF5ljQKRH0KK8r0ajvUuVmTzydgyQiEjD2USRumQvnERIfDDbaeZB0iVItR/TZDDAqPPSMkTtBQMkItJQLlYrZZD0yF44OyYfzHbaUWw6nENxTJPRHiAxTiVvxwCJiDSUi9Xai7RbJ4PE2hc7zVIjOt8D5ahFBqrk7RggEZGG1ckwf31GsWm38cFspx3FpkMXm+IcJqMB9QkkdnWS12OAREQaUqG0rculfpuOo9j4YHZOuxab68/hOIpNCoiZySNvxwCJiDQsqmH++k0UqSwGt8/Y7fLTtFuaiSJ1vAcGA2AwKEctuvxURO0KAyQi0lAVabdCDZLBAPvoKWYuZK1RgyRnC+s/f4OOATFRe8IAiYg0lIvV6rnUiHRI26K46nNT69QgydnC+vtsMjo/N5G38XF3A4io7ZFqXQw6z4sjj5YzGCCtAsbMhZ1jHZAuhfJWezcnAF27VInaE2aQiEjDXhuk6PrScS02owG6TifQXkmfubSIrJ4zaUvF+PYAyeWnImpX3B4gvfnmm+jduzcCAgKQkJCAHTt2NLr/hg0bMHDgQAQEBGDo0KHYtGmT6nUhBBYuXIioqCgEBgYiOTkZR48eVe3zxz/+EUlJSQgKCkLHjh1dfUlE7Z66i822TY+Hs1XRvaNnV157JX0W0iKyetRnWTVdbAxUiQA3B0gffPABMjMzsWjRIuzZswfDhw9HSkoKiouLne6/fft23HvvvZg2bRry8/ORmpqK1NRU7Nu3T95nyZIlWLlyJVatWoW8vDwEBwcjJSUFVVVV8j41NTWYPHkyHnnkEd2vkag9shdPG3QtnpaeweolLlx+mnZLDpBM+mV1lFMtANB1Wgei9sStAdKyZcswffp0ZGRk4KqrrsKqVasQFBSE1atXO91/xYoVmDBhAmbPno1BgwbhxRdfxMiRI/HGG28AsP0f7vLlyzF//nzceeedGDZsGN59912cPn0aGzdulI/z/PPP46mnnsLQoUNb4zKJ2h05cDHqO7Oys0CMD2Y76bPwrU/j6TOTtu1fuQZJyuQxg0Rezm0BUk1NDXbv3o3k5GR7Y4xGJCcnIzc31+l7cnNzVfsDQEpKirx/QUEBzGazap+wsDAkJCQ0eMymqq6uRllZmeqLyFOp5ifSsevLnr3g6ClnpODFR8d7oFxWBrAHSpxugbyd2wKkc+fOwWKxICIiQrU9IiICZrPZ6XvMZnOj+0v/NueYTZWVlYWwsDD5KyYm5oqOR9SWKecnMug4ik2ag8eonCiSD2aZcMgg6XIPFIXytn85YScR0AaKtNuLuXPnorS0VP46efKku5tEpBtlXYpJx1FNyjXfOIpNS/osfORRbK4/hzJbCCi62Biokpdz2zxI4eHhMJlMKCoqUm0vKipCZGSk0/dERkY2ur/0b1FREaKiolT7xMfHX1F7/f394e/vf0XHIGov1IvV1m/To/5FHkEFSDMhMT6yq6v/MPxM+o1ic+xikzJJ7GIjb+e2DJKfnx9GjRqFnJwceZvVakVOTg4SExOdvicxMVG1PwBs3rxZ3j82NhaRkZGqfcrKypCXl9fgMYlIS+76Muo8is3aOrVO7ZVjkbaeXWwmh2H+XKyWvJ1bZ9LOzMxEeno6Ro8ejTFjxmD58uWoqKhARkYGAGDq1Kno3r07srKyAABPPvkkxo0bh6VLl2LSpElYv349du3ahbfffhuArVZi5syZeOmllxAXF4fY2FgsWLAA0dHRSE1Nlc9bWFiICxcuoLCwEBaLBXv37gUA9OvXDyEhIa36GRC1RRZnRdq6PJxt/xoNBhikzAVTSDIpg+TrUz+KTZcMku1fx5m0Od0CeTu3BkhpaWk4e/YsFi5cCLPZjPj4eGRnZ8tF1oWFhTAa7UmupKQkrFu3DvPnz8e8efMQFxeHjRs3YsiQIfI+c+bMQUVFBR588EGUlJRg7NixyM7ORkBAgLzPwoUL8be//U3+ecSIEQCArVu34oYbbtD5qonaPqtidmW5NkiHh3NdfarKZLQHSBxebmepX/PF30e/DFKDXWy8D+Tl3L4W24wZMzBjxgynr23btk2zbfLkyZg8eXKDxzMYDHjhhRfwwgsvNLjP2rVrsXbt2uY2lchrKGdX1nNUk70I2Sh36bD2xU7KIOkZIAmHUWzsYiOy4Sg2ItKwKLq+9Bx+LwUAPspAjM9lmfSZ+/uYANg/L1eyKAryAcW0DgyQyMu5PYNERG2PfS022wgzQKfuHYu6QFh5bmqdDJLFqr4HJh3nvSJqTxggEZGGapi/jhkFqQbJx2iAdHQ+mO0s9Z+PFCDV6lA53fAoNpefiqhdYRcbEWnIhbsGgzxJoR6BS50ie2Hi/Dsa0mfu76tfBqnOIYvH0YRENgyQiEhDmVWQ5uCRHqSuJD2EfU1GuQaGxcF2coCkYw1SneIeAPZAiYEqeTsGSESkoZyfyKe+CKlGh+4dZfbCvtSIy0/TbrVGDZLUbedrUtcgMVAlb8cAiYg0lPMg+fnYHph61L9YFKPYuIq8lsUhQNLjHtTUSQGS7RzSKDY9slVE7QkDJCLSqJWKp01GXbvYahUTReo5Y3d7JQVIfjpmkBy72Hx1rDkjak8YIBGRRm2dfZFUH5N+XWzK1ertXWx8MEvqWqEGybGLTQqUpMwSkbdigEREGvLwe5NBfnDq0b0jZaV8jEZ5JmcuVmsnBYsBvlIWz/X3oNaiziBJoxbZxUbejgESEWkoH5q+9UXaQrg+u+OsBokBkp1jF5tVuL4LUgp8pWJ8P5N+wRhRe8IAiYg0lN0u0kryyu2uopwHiaPYtByH+QOuL2Kvre9Kk4rxpQxSjQ41Z0TtCQMkItKQu75MRvgolgFxdR2SlKXwMXEUmzOOw/wB12fxaq32bk7lv8wgkbdjgEREKkIIeXSZr8k+USTg+pFsdYqHM0exaUmfhTSTNuD6LJ49W1jfxVYfjLEGibwdAyQiUrFYhbwOl2994CIFL65+OCsXSjWwSFujTl6LTdHF5uLARcoU+UpdbPX3mqPYyNsxQCIiFWXmQKo/0uuhWacs0uYwfw3HIm3A9ZkduSDfKM2DJGWQGCCRd2OAREQqyiyRFBjJI5t0yl6YTBzF5oxFsSaeXgFkjUMXmzStgx4TgxK1JwyQiEilVvFgdJwbR69RbD4cxeaUxaLNsLn8HigK5W3/6jcxKFF7wgCJiFSkB6bRYF/ZXQqU9KpB4kSRzimnQfDVKYMkBcR+JocuNmaQyMsxQCIilRo5o2D/82APkPQaxca12JxRFrFLn4/ra5AclxrRb+Z0ovaEARIRqdQ5ZBQA/R6aFqu9BsnIeZA0pM/Cx2iQA1bXZ5DUAbFewTBRe8MAiYhUlOuwSfTqYqtzUmOjVwap1mJFxpodWJlz1CXHKy6vwuRV2/Gv/F9ccjxnpBokk2KeKL3ugRQQ+8iZKmaQyLsxQCIilZo69czKgD27oFcXm8lo1D2DtO3wWWw9fBbLNh9xyfGyNh3CzhO/4qkPvnfJ8ZyRPx+DQQ5c9BrF5mPSt96MqL1hgEREKlLmwE+RQZK+d/U8SNLD3lfRxdYaJUjlVbVXfIyCcxUuaEnj5GH+JoMcwOhXg8QuNiIlBkhEpFKrWIdNEuBrm8m5qtbi0nNJwZitCNm2Ta8uNuXouDOlVVd8PFcEWZdjURSxSxk9V2eQLtXa7kFg/T2WAzFmkMjLMUAiIhXHUU0AEOhne3hecnGAJGWkbEua6LsG2MWqOvn70yWXrvh45Yrj6UEIIQdDRoNBvh+uzuJV1djuqXSP/ZhBIgLAAImIHEhFu8pFaoOkAKnGtQFSdf3D3t/XKC+noVftS0WNPaC5UFFzxcfTeyJF5fH9fY3yemzVda69B5W1ts9FCpD0mhSUqL1hgEREKrVORrFJXWyuziBJXXYBvib41wdI1bX6PJiVGZ8qF5xDKpoGXB+02I5pb6OfyYgAX9vn44q2K12qUXexsUibyIYBEhGpSAGKcgV56eGpWwbJxyh37eiVmblYbQ+QXBHoSUXlAFBR7foASdmV5u+jXwbpUn1mTbrHUqDKpUbI2zFAIiKVS/VdLlK3mvJ7XTNI9RkSPbIxAFBRrcwgXfk5lAGE8tiuIgWPfiYjDAaDnEFyZYZNCCHfU+keB+rUnUrU3jBAIiIVqctF6lYDWjmD5OIiZImySNsV16E8hh4F2zWKz8b2r+tHElbXWeVpFQKkAEmne03U3jBAIiIVKaMQqAyQ/HxUr7mKVE+jziDpEyAp236lQYbVKlTtdPXnAtgzaVLxuh6fj/JzkO53oKLeTHDZF/JiDJCISKXKWYBU/3B2fQbJdjxbBkmqsdE/QLrSgMaxja6eHwpoLIPkus9H+hx8TQa5OFvqYrMK1iGRd2OAREQqlTXqYd/K712dKalWZpB89M0gVbkwQHJ8vx4BklyDVP+5BOhQo1VZY68Bkyi/ZzcbeTMGSESkIg/79tN2sVXWuK7WxmoVcobC38c+D5JeNUiXFJmXKw1oHAMkPbrY7Bkkk+pfl2aQarTZQl+TUZ6UUo/rImovGCARkYqzGqQOAbYAqeyS6wIkZaZInUHS56FcpciGXGlmxPH9emRa5O5HX/0ySGWXbMulhAb6qrYHsFCbiAESEak5q0EKq3+All5y3fpjyixO62SQlEXaV3YOxwyUnjVI0ug+PTJI0v0McwiQ9JrWgag9YYBERCpS1iDATxsglbkwQJIySD5GA3xMyokQ234NkjZAcn2blcuwAPpkkBoKkAJ1WpyYqD1hgEREKpVOMkgd6x+g5dV1LlvlXXr4Sl1rrZtBavs1SNWtmEHq2EAXWyW72MiLMUAiIhVpVuhgRQZJWaNS5qJJEaXJFTsE2I6tew2SK0exOdYg6RggSYGRVAdWXuW6LF5pAzVIwf62c+kxQzhRe8EAiYhU5G6XIPtD09dklAMmV9UhlVVJD2fbw1g5zN/VExTWWqyotdiPecVF2o4ZJB0yLTUOw/ylIMZVASrQcBebHjVnRO0NAyQiUmnoodkxyA8A8GtljUvOI4+gqs8gSVkLIfRb862hn5vLcT00PWp1pKVRpM8lVB5J6LqgpaSyvostyOFeM0AiYoBERHZCCJTKD00/1WvdQv0BAMVlVS45l5RBkgKxID8TTEbb/DuuXttMO7HjldXxOB5Pj1odqStNCoz0KJQvqr+XEaEBqu1S9lAKoIi8EQMkIpJdqrXIkzc6Fu5GdLA9RIvKql1yLmlOJanryGAwIMTf9XU2gHqhWsC2hMaVFJs7Bkh61OrYa7TqM0iKQnmL1TVdkEXlUoDkr9ouBWMlzCCRF2OAREQyqUvFx2iQ58KRRIZJAZJrM0hShgSw1yOVunBCSgC4WB/AdA62Z8WqrmC0nFRzJB2vXIcAyV6jZQtWOig+J8eAryWEEHKw262DOoPELjYiBkhEpCB1qYQF+sJgMKhek7rYzC4KkJyNoOrgX58lcXEGScrGdFEESFdSWC3VHIWH+KmO70qOGSR/H5M8F1LJpSuvA/u1slYuBO/mmEGSu9hcU29G1B4xQCIiWXG5LaPQtYO/5rWYTkEAgMLzla45V5n2XFIGyZUjtQB1sOGKSRCl90ptv1jt+kyLPcNmDyCl850tv/Juzp/PVwCwda9JUwlIwkNcdx6i9ooBEhHJzKWXANi705Riw4MBAAXnKlxyrjPSuRQFwtKcSK7PINXKx5cW4b2SAEmqQZK6pvTNINkDpKjQQADA6dIrz+JJ91G6r0pR9ff/jAvOQ9ReMUAiIpm51JYxiAzVBki96x+k5ytqXNL1Ij18o8IC5W1S7cuvFa7t2pFqkEIUGaQrmUrgUv0oODmDVFXn8rmbpOkUlNMtSIGrFMheCXuAFKJ5TRrVVl5Vx8kiyWsxQCIimbms4QxSiL8Pena2dbN9/0vpFZ2n1mLF2Yu2YCyqo/1ccgDgojoniVTUHBrgI9fxXMnQ/Es1tuN1re+KqrMKl64hd6nGIteDKe+F9FmdLrnyz2f/6TIAwIAIbYDUIcBXHlHo6ntB1F4wQCIi2YlztvqiHvX1Ro5G9+oEANh94sIVnqcCQtiWM1EWTkvZpDMuCACUpILwDgG+clH4lYzQkt4bGRaA+qmbXDriSwpKgvxMqlF+vTrbsng/XWE3pxACP56yBblDe3R0uk+PTrZ7IdUqEXkbBkhEJDt29iIAoF83bVYBAEb1rg+QCn+9ovMcMpcDAPpHdlCNltOr9qWoXBrO7o9O9RNgXkk34a/12Z0uIX5yN5urpj8AFPVZYQGqz2dApO2+HKn//Fqq4FwFzpZXw89kxFVRoU736Vv/O3Cs+OIVnYuovWKAREQAgNLKWnnUUt+u2sJdABjdqzMAIL+w5IqKnO3dOx1U27vXZy1OXqh0aU2PsvhcCpAuVLQ84yPVSHUK8kOklPVyYVB38oItk9e9Y6Bqe//6z8tcVoVzF1s+wuzLI2cBAKN6dZKL1h3162oLkI4WMUAi78QAiYgAAHt/KQEA9OwcpBo5pRTXLQTdOwaissaCz/ebW3yuvILzAGwPaKXeXYLhazKgvLoOp0quvBBZInVZRYYGoNMVzvFjtQp5hunOwX6IDHV9BungmfoMm0MA2SHAFwMjbdvyfmp5N+dHe04BAG6+KqLBfQZH2zJL39f/XhB5GwZIRAQA2FVfVzTaIWhRMhoNuHtUDwDABztPtug8xeVV+P5kCQAgqV+46jU/HyP61mcupCzTlaqzWFFUPzovIjQAnYKlDFLLAqTyKvtSHx2DfOW6qVO/ui6gO1B/7YOcdH9dW/+Zbdp3psXH/vFUKXxNBqSO6N7gfiPrfw+OFF1s8WdF1J4xQCIiAMB/DxYDABL7dml0v8mjesBkNGD78fP4+ujZZp/nH7t/gVUA8TEdNV1IgD2rtP3YuWYf25mjxRdRY7EixN8H3TsGyl1s51v40D9fYQu2gv1M8PcxyVkeVwV0FdV1yD9pq/FyzLABwN0jbQHqF/vNze5mE0JgcfYhAMAtgyNVS684Cg/xx5DutgDtsx9bFowRtWcMkIgI+0+X4uCZMvgYDRg/qOFuFwCI6RyEqYm9AADP/PPHZnUtFZdV4e2vfgIA3H9NL6f7jOvfFQDw+f6iK1pQVrKnvqB8aPcwGI0G9O5iG6HX0gkvfzpre1+vLrY6reExYQCA70+WoNYF7f3vwSLUWgR6dg6S26p0VXQo4mM6otYisPSLw8069ru5P+OrI2fhZzJi9i0DLrv/ncNtGaaN+aeadR4iT8AAicjLCSGw9IsjAICUIY1nFSRP3dwffcKDcarkEu5841tszD912ckdT16oxNTVO1BSWYtBUaG4Mz7a6X7jBnRF52A/mMuq8MGulnXjKX3y/WkAwLX9bJmxfvXz/vx8vqJFheZHim31QXH1xxkQ0QHhIX4or67Dl4ebn1FTsliFHEDeNaK7Zj08yR8mDITBALy/4yTe2HK00UDSahU4WlSO+Rt/xKJ/7wcAZN7SX574szG3D4+GyWjA7p9/xdbDxS24IqL2y+fyu+jvzTffxCuvvAKz2Yzhw4fj9ddfx5gxYxrcf8OGDViwYAFOnDiBuLg4vPzyy5g4caL8uhACixYtwl//+leUlJTg2muvxV/+8hfExcXJ+1y4cAGPP/44PvnkExiNRtx9991YsWIFQkKcD28m8kQWq8DL2Yew5VAxfE0GzBwfd/k3wbY+2N/+ZwzSV+/AT+cqMPODvQBsa511DPJFpyA/dAjwQZCfbeZqc2kV9hT+ijqrQHiIH966bxR8TM7//8zfx4THbuyHFz89gMWbDmF4j44Y0j2sRdf38d5T+O6nC/Ax2uttuob4o3OwHy5U1CC/sOSyXYqO9vxsy0hJXWs+JiPuGtEdf/26AIuzD+Hq3p3lxV6bQwiBFz89gP2ny9DB30fO0jmT2LcLnhwfh+X/PYpXvziC//uuEEN7hCHIz4SLVXUoq6pFSWUtLlTU4NfKGlgVAwIfHtcXD13fp0ltigwLwP9c2xt//boAszd8j7UZY1p8L4jaG4Nw9fz4zfTBBx9g6tSpWLVqFRISErB8+XJs2LABhw8fRrdu3TT7b9++Hddffz2ysrJw2223Yd26dXj55ZexZ88eDBkyBADw8ssvIysrC3/7298QGxuLBQsW4Mcff8SBAwcQEGCbZ+XWW2/FmTNn8NZbb6G2thYZGRm4+uqrsW7duia1u6ysDGFhYSgtLUVoqPN5RIjaGqtV4OcLldh3qhT7TpXi0x/OyKPFsn4zFPeO6dms41VU1+GdbwqwMf9UkyYvTOrbBYt/Mww9nXQdKVmsAmlv5WLXz78iwNeI9KTeuG1oNOIiQhDg63xYuvK9R4vL8Y9dv2DN9hOwWAUeGtcHc28dJO8ze8P32LD7F9w6JBJ/njKywUyNo5/PVyB52ZeotQh88dT1cpB0/mI1Jqz4GmfLq9EnPBiPj++HGwd0Q8egxrNxQgicu1iDnScu4G/bTyCvwFYovzwtvtECasn6HYVY8vnhyxZR+/sYcU2fLnjo+j6awvjLqaiuw29X5eLgmTIYDcANA7phXP+uGBwdioFRofKM20TtRVOf324PkBISEnD11VfjjTfeAABYrVbExMTg8ccfxzPPPKPZPy0tDRUVFfj000/lbddccw3i4+OxatUqCCEQHR2Np59+GrNmzQIAlJaWIiIiAmvXrsU999yDgwcP4qqrrsLOnTsxevRoAEB2djYmTpyIX375BdHRzlP/SnoFSNn7zDhdcgm+Pkb4mQzwNRnlL38f6XtD/esN95A63lUB0eDrjr8Ajr8SQvWa5kzNOG7TzqE9juZimtw+5Xu1rzmeUzT4mvacjRy30fbp99lahIDFaoXFClisVtRYBOosVpRdqsMvv1bixPkKHDxTLq9LJgkL9MVzd1yFu0b0cGxAs5RV2eZRKqmsRUllDS5W16Gi2oLKmjp07eCPwdFhDU5A6Ux5VS0e/fsefH1UXaztYzTAx2SAj9EIgwEwGgwwGmyfiRBAZU0dai32T+i3o3pgyd3DYDTag6B9p0px+xvfQAhbzVPqiGj07hKMYH8fBPiYYBWi/vO0fZ2/WIMfT5Vi9bcFOFtejaS+XbBu+jWqdh0yl+H+d3bIc0kBtuVZfEwGGGBrn9Uq5HtVZ7XdH2V2J8DXiJdSh+K3o5p+L6pqLdhRcAGFFypxqcaCDgE+CA30RVigLzoH+6FLsB86B/s1mLFritLKWjzz0Q/4zz7t1A69uwShT9cQ9OwchO4dAxHoZ4K/jxE+JgOMBgNMRgN8jLbvfX2M8DcZ5b9fpvp7IsWnBhjs3xvsPxsU+9h+Uv6s/G/D2ePMvqP9PIpXFQcyOOwntclxmytc7smr+bvX7Pdf7vyXOf5lzy1gFbbvrULAKoTie9T/XP+9VdrX/lpj+9fUWVFRU4fKagsmDYtCtJPBHFeiqc9vt4b+NTU12L17N+bOnStvMxqNSE5ORm5urtP35ObmIjMzU7UtJSUFGzduBAAUFBTAbDYjOTlZfj0sLAwJCQnIzc3FPffcg9zcXHTs2FEOjgAgOTkZRqMReXl5uOuuuzTnra6uRnW1/Q9fWZlrRqw4WrejEF8dubI6BqLG+PkYMSiyA66KDsM1fTojZXDkZbMyTREa4IvQBuZPaokOAb5493/G4IsDRfjH7l/w3U/nUV5VZwssrAJAw3U3fj5GXB8Xjt8n9MRNA7VF50O6h2Hxb4bi2X/tw5dHzsoTJzZFn67BeHXycM32gZGh+G/mOKz+pgCf/HAaP52t0ASjzhgMQGyXYNw8OAL3X9OrwWVeGhLga8L19YXtegkL8sVf7huFo0Xl+Hy/GXsKS3DgdBnMZVU4cb4SJ85X6np+8l6DokJdHiA1lVsDpHPnzsFisSAiQv0HLCIiAocOHXL6HrPZ7HR/s9ksvy5ta2wfx+47Hx8fdO7cWd7HUVZWFp5//vkmXlnLXdOnM8ICfVFbZ0WtxYoaixU19d/XWoS8rbZ+uyODk/9bsr+m+L6R/x3SvM/J/0215ByGBn9o+vu052v4etXtbMYxm3EO9TU1p90tbBsMzneELbNiqv8/dR+jLfvoYzIg2M8HPToFIqZzEAZFhaJv1+Aryia0JoPBgJTBkUgZHAkhBEov1aKq1vb7r/w/T2t9GsZgAAL9fBAZGiBnJxqSdnVPjOzZCR/sPInvfynB6ZIqVNbUoarWWv85Aqb6zzQs0Bex4SG4aWA33D2qO/x9nAeUYYG+eOrm/njq5v6oqK6DuaxK/r9io5ztsmVFfExG+BoNCA30dUmA2hriIjogTjF55fmL1ThkLseJ8xUoPF8Jc1kVqmotqK6zytk36csW2Nr+btn+ptl+BuzZEClrIWUohJAyg0J+XQj1z/asT2OZJfV77Bu13zrLJNvbp84cN/o3x/E/UOVrjb6vkdca+7vdyPsae7HF51Nkbw31/yp/v42KbQZnrxml15TvtX1uvj62v1sh/j7oFOy6/+lqLnYeN9HcuXNVmauysjLExMS4/DyP3tDP5cck8gQGg+GyNT3NFRfRAfNvu8qlx5QE+/vIk156qi4h/ri2n788eSWRJ3Hr/0aGh4fDZDKhqKhItb2oqAiRkZFO3xMZGdno/tK/l9unuFg9ZLWurg4XLlxo8Lz+/v4IDQ1VfREREZFncmuA5Ofnh1GjRiEnJ0feZrVakZOTg8TERKfvSUxMVO0PAJs3b5b3j42NRWRkpGqfsrIy5OXlyfskJiaipKQEu3fvlvfZsmULrFYrEhISXHZ9RERE1D65vYstMzMT6enpGD16NMaMGYPly5ejoqICGRkZAICpU6eie/fuyMrKAgA8+eSTGDduHJYuXYpJkyZh/fr12LVrF95++20AtjT8zJkz8dJLLyEuLk4e5h8dHY3U1FQAwKBBgzBhwgRMnz4dq1atQm1tLWbMmIF77rmnSSPYiIiIyLO5PUBKS0vD2bNnsXDhQpjNZsTHxyM7O1susi4sLITRaE90JSUlYd26dZg/fz7mzZuHuLg4bNy4UZ4DCQDmzJmDiooKPPjggygpKcHYsWORnZ0tz4EEAH//+98xY8YMjB8/Xp4ocuXKla134URERNRmuX0epPaKE0USERG1P019frePsb5ERERErYgBEhEREZEDBkhEREREDhggERERETlggERERETkgAESERERkQMGSEREREQOGCAREREROWCAREREROTA7UuNtFfSBORlZWVubgkRERE1lfTcvtxCIgyQWqi8vBwAEBMT4+aWEBERUXOVl5cjLCyswde5FlsLWa1WnD59Gh06dIDBYHDZccvKyhATE4OTJ0967Bpvnn6NvL72z9Ov0dOvD/D8a+T1tZwQAuXl5YiOjobR2HClETNILWQ0GtGjRw/djh8aGuqRv/RKnn6NvL72z9Ov0dOvD/D8a+T1tUxjmSMJi7SJiIiIHDBAIiIiInLAAKmN8ff3x6JFi+Dv7+/upujG06+R19f+efo1evr1AZ5/jbw+/bFIm4iIiMgBM0hEREREDhggERERETlggERERETkgAESERERkQMGSG5y4sQJTJs2DbGxsQgMDETfvn2xaNEi1NTUqPb74YcfcN111yEgIAAxMTFYsmSJ5lgbNmzAwIEDERAQgKFDh2LTpk2tdRmN+uMf/4ikpCQEBQWhY8eOTvcxGAyar/Xr16v22bZtG0aOHAl/f3/069cPa9eu1b/xTdSUaywsLMSkSZMQFBSEbt26Yfbs2airq1Pt05avUal3796a+7V48WLVPk35nW3L3nzzTfTu3RsBAQFISEjAjh073N2kFnnuuec092rgwIHy61VVVXjsscfQpUsXhISE4O6770ZRUZEbW3x5X331FW6//XZER0fDYDBg48aNqteFEFi4cCGioqIQGBiI5ORkHD16VLXPhQsXMGXKFISGhqJjx46YNm0aLl682IpX0bDLXd8DDzyguacTJkxQ7dOWry8rKwtXX301OnTogG7duiE1NRWHDx9W7dOU38um/E11CUFu8Z///Ec88MAD4vPPPxfHjx8XH3/8sejWrZt4+umn5X1KS0tFRESEmDJliti3b594//33RWBgoHjrrbfkfb799lthMpnEkiVLxIEDB8T8+fOFr6+v+PHHH91xWSoLFy4Uy5YtE5mZmSIsLMzpPgDEmjVrxJkzZ+SvS5cuya//9NNPIigoSGRmZooDBw6I119/XZhMJpGdnd1KV9G4y11jXV2dGDJkiEhOThb5+fli06ZNIjw8XMydO1fep61fo1KvXr3ECy+8oLpfFy9elF9vyu9sW7Z+/Xrh5+cnVq9eLfbv3y+mT58uOnbsKIqKitzdtGZbtGiRGDx4sOpenT17Vn794YcfFjExMSInJ0fs2rVLXHPNNSIpKcmNLb68TZs2iWeffVZ89NFHAoD417/+pXp98eLFIiwsTGzcuFF8//334o477hCxsbGqvykTJkwQw4cPF9999534+uuvRb9+/cS9997bylfi3OWuLz09XUyYMEF1Ty9cuKDapy1fX0pKilizZo3Yt2+f2Lt3r5g4caLo2bOn6m/I5X4vm/I31VUYILUhS5YsEbGxsfLPf/7zn0WnTp1EdXW1vO0Pf/iDGDBggPzz7373OzFp0iTVcRISEsRDDz2kf4ObaM2aNY0GSI5/BJTmzJkjBg8erNqWlpYmUlJSXNjCK9fQNW7atEkYjUZhNpvlbX/5y19EaGiofF/byzUKYQuQXnvttQZfb8rvbFs2ZswY8dhjj8k/WywWER0dLbKystzYqpZZtGiRGD58uNPXSkpKhK+vr9iwYYO87eDBgwKAyM3NbaUWXhnHvx1Wq1VERkaKV155Rd5WUlIi/P39xfvvvy+EEOLAgQMCgNi5c6e8z3/+8x9hMBjEqVOnWq3tTdFQgHTnnXc2+J72dH1CCFFcXCwAiC+//FII0bTfy6b8TXUVdrG1IaWlpejcubP8c25uLq6//nr4+fnJ21JSUnD48GH8+uuv8j7Jycmq46SkpCA3N7d1Gu0Cjz32GMLDwzFmzBisXr0aQjE1V3u/vtzcXAwdOhQRERHytpSUFJSVlWH//v3yPu3pGhcvXowuXbpgxIgReOWVV1Sp7ab8zrZVNTU12L17t+peGI1GJCcnt9l7cTlHjx5FdHQ0+vTpgylTpqCwsBAAsHv3btTW1qqudeDAgejZs2e7vdaCggKYzWbVNYWFhSEhIUG+ptzcXHTs2BGjR4+W90lOTobRaEReXl6rt7kltm3bhm7dumHAgAF45JFHcP78efm19nZ9paWlACA/95rye9mUv6muwsVq24hjx47h9ddfx6uvvipvM5vNiI2NVe0n/VKYzWZ06tQJZrNZ9Ysi7WM2m/VvtAu88MILuOmmmxAUFIQvvvgCjz76KC5evIgnnngCABq8vrKyMly6dAmBgYHuaHaTNdR+6bXG9mmL1/jEE09g5MiR6Ny5M7Zv3465c+fizJkzWLZsGYCm/c62VefOnYPFYnF6Lw4dOuSmVrVcQkIC1q5diwEDBuDMmTN4/vnncd1112Hfvn0wm83w8/PT1M21p78djqR2N/b30Gw2o1u3bqrXfXx80Llz53Zx3RMmTMBvfvMbxMbG4vjx45g3bx5uvfVW5ObmwmQytavrs1qtmDlzJq699loMGTIEAJr0e9mUv6muwgDJxZ555hm8/PLLje5z8OBBVbHkqVOnMGHCBEyePBnTp0/Xu4lXpCXX15gFCxbI348YMQIVFRV45ZVX5ADJHVx9jW1dc643MzNT3jZs2DD4+fnhoYceQlZWlscuedBe3XrrrfL3w4YNQ0JCAnr16oUPP/ywTQXd1HT33HOP/P3QoUMxbNgw9O3bF9u2bcP48ePd2LLme+yxx7Bv3z5888037m5KgxggudjTTz+NBx54oNF9+vTpI39/+vRp3HjjjUhKSsLbb7+t2i8yMlJTvS/9HBkZ2eg+0uuu1tzra66EhAS8+OKLqK6uhr+/f4PXFxoaqtsfeVdeY2RkpGYUVFPvoZ7XqHQl15uQkIC6ujqcOHECAwYMaNLvbFsVHh4Ok8nUqv89taaOHTuif//+OHbsGG6++WbU1NSgpKRE9X/r7flapXYXFRUhKipK3l5UVIT4+Hh5n+LiYtX76urqcOHChXZ53X369EF4eDiOHTuG8ePHt5vrmzFjBj799FN89dVX6NGjh7w9MjLysr+XTfmb6ioMkFysa9eu6Nq1a5P2PXXqFG688UaMGjUKa9asgdGoLglLTEzEs88+i9raWvj6+gIANm/ejAEDBshdFYmJicjJycHMmTPl923evBmJiYmuuSAHzbm+lti7dy86deokZyMSExM10xboeX2Aa68xMTERf/zjH1FcXCynvjdv3ozQ0FBcddVV8j6tfY1KV3K9e/fuhdFolK+tKb+zbZWfnx9GjRqFnJwcpKamArB1A+Tk5GDGjBnubZwLXLx4EcePH8f999+PUaNGwdfXFzk5Obj77rsBAIcPH0ZhYWGr/d65WmxsLCIjI5GTkyMHRGVlZcjLy8MjjzwCwPb7WVJSgt27d2PUqFEAgC1btsBqtSIhIcFdTW+xX375BefPn5cDwrZ+fUIIPP744/jXv/6Fbdu2abrjm/J72ZS/qa5sMLnBL7/8Ivr16yfGjx8vfvnlF9WwTUlJSYmIiIgQ999/v9i3b59Yv369CAoK0gzz9/HxEa+++qo4ePCgWLRoUZsZ5v/zzz+L/Px88fzzz4uQkBCRn58v8vPzRXl5uRBCiH//+9/ir3/9q/jxxx/F0aNHxZ///GcRFBQkFi5cKB9DGgI/e/ZscfDgQfHmm2+2qSHwl7tGaUjqLbfcIvbu3Suys7NF165dnQ7zb6vXKNm+fbt47bXXxN69e8Xx48fF//3f/4muXbuKqVOnyvs05Xe2LVu/fr3w9/cXa9euFQcOHBAPPvig6Nixo2rETHvx9NNPi23btomCggLx7bffiuTkZBEeHi6Ki4uFELbh1D179hRbtmwRu3btEomJiSIxMdHNrW5ceXm5/N8YALFs2TKRn58vfv75ZyGEbZh/x44dxccffyx++OEHceeddzod5j9ixAiRl5cnvvnmGxEXF9dmhsE3dn3l5eVi1qxZIjc3VxQUFIj//ve/YuTIkSIuLk5UVVXJx2jL1/fII4+IsLAwsW3bNtUzr7KyUt7ncr+XTfmb6ioMkNxkzZo1AoDTL6Xvv/9ejB07Vvj7+4vu3buLxYsXa4714Ycfiv79+ws/Pz8xePBg8dlnn7XWZTQqPT3d6fVt3bpVCGEbfhofHy9CQkJEcHCwGD58uFi1apWwWCyq42zdulXEx8cLPz8/0adPH7FmzZrWv5gGXO4ahRDixIkT4tZbbxWBgYEiPDxcPP3006K2tlZ1nLZ8jZLdu3eLhIQEERYWJgICAsSgQYPEn/70J9UfZyGa9jvblr3++uuiZ8+ews/PT4wZM0Z899137m5Si6SlpYmoqCjh5+cnunfvLtLS0sSxY8fk1y9duiQeffRR0alTJxEUFCTuuusu1f+gtUVbt251+t9benq6EMI21H/BggUiIiJC+Pv7i/Hjx4vDhw+rjnH+/Hlx7733ipCQEBEaGioyMjLk/6Fxt8aur7KyUtxyyy2ia9euwtfXV/Tq1UtMnz5dE7y35etr6Jmn/HvXlN/LpvxNdQVDfaOJiIiIqB7nQSIiIiJywACJiIiIyAEDJCIiIiIHDJCIiIiIHDBAIiIiInLAAImIiIjIAQMkIiIiIgcMkIiIiIgcMEAioia54YYbVGv+udPatWtVi1nq7frrr8e6deta7XzuUlNTg969e2PXrl3ubgqR2zFAIiJd6BnEpKWl4ciRI7oc29G///1vFBUV4Z577lFtz8/PR1paGqKiouDv749evXrhtttuwyeffIL2ukCBn58fZs2ahT/84Q/ubgqR2zFAIqJ2JzAwUF7JW28rV65ERkYGjEb7n8uPP/4Y11xzDS5evIi//e1vOHjwILKzs3HXXXdh/vz5KC0tbZW26WHKlCn45ptvsH//fnc3hcitGCARkUZFRQWmTp2KkJAQREVFYenSpZp9qqurMWvWLHTv3h3BwcFISEjAtm3bAADbtm1DRkYGSktLYTAYYDAY8NxzzwEADAYDNm7cqDpWx44dsXbtWgDAiRMnYDAY8NFHH+HGG29EUFAQhg8fjtzcXHl/x+zUc889h/j4eLz33nvo3bs3wsLCcM8996C8vFzep7y8HFOmTEFwcDCioqLw2muvXbbb8OzZs9iyZQtuv/121Wczbdo0TJo0CZ999hluueUW9OnTB4MGDcK0adPw/fffIywsDABgsVgwbdo0xMbGIjAwEAMGDMCKFStU53jggQeQmpqKP/3pT4iIiEDHjh3xwgsvoK6uDrNnz0bnzp3Ro0cPrFmzRn6P9Bl9+OGHuO666xAYGIirr74aR44cwc6dOzF69GiEhITg1ltvxdmzZ+X37dy5EzfffDPCw8MRFhaGcePGYc+ePar2dOrUCddeey3Wr1/f4OdC5A0YIBGRxuzZs/Hll1/i448/xhdffIFt27ZpHqQzZsxAbm4u1q9fjx9++AGTJ0/GhAkTcPToUSQlJWH58uUIDQ3FmTNncObMGcyaNatZbXj22Wcxa9Ys7N27F/3798e9996Lurq6Bvc/fvw4Nm7ciE8//RSffvopvvzySyxevFh+PTMzE99++y3+/e9/Y/Pmzfj666811+Tom2++QVBQEAYNGiRv++KLL3D+/HnMmTOnwfcZDAYAgNVqRY8ePbBhwwYcOHAACxcuxLx58/Dhhx+q9t+yZQtOnz6Nr776CsuWLcOiRYtw2223oVOnTsjLy8PDDz+Mhx56CL/88ovqfYsWLcL8+fOxZ88e+Pj44Pe//z3mzJmDFStW4Ouvv8axY8ewcOFCef/y8nKkp6fjm2++wXfffYe4uDhMnDhRFUgCwJgxY/D11183+tkQeTxBRKRQXl4u/Pz8xIcffihvO3/+vAgMDBRPPvmkEEKIn3/+WZhMJnHq1CnVe8ePHy/mzp0rhBBizZo1IiwsTHN8AOJf//qXaltYWJhYs2aNEEKIgoICAUD87//+r/z6/v37BQBx8OBBp8detGiRCAoKEmVlZfK22bNni4SEBCGEEGVlZcLX11ds2LBBfr2kpEQEBQXJ1+TMa6+9Jvr06aPatnjxYgFAXLhwQd62Y8cOERwcLH998sknDR7zscceE3fffbf8c3p6uujVq5ewWCzytgEDBojrrrtO/rmurk4EBweL999/v8HP6P333xcARE5OjrwtKytLDBgwoMG2WCwW0aFDB017V6xYIXr37t3g+4i8gY87gzMianuOHz+OmpoaJCQkyNs6d+6MAQMGyD//+OOPsFgs6N+/v+q91dXV6NKli0vaMWzYMPn7qKgoAEBxcTEGDhzodP/evXujQ4cOqvcUFxcDAH766SfU1tZizJgx8uthYWGqa3Lm0qVLCAgIaFJb9+7dCwCIi4tTZbrefPNNrF69GoWFhbh06RJqamoQHx+vev/gwYNVNU4REREYMmSI/LPJZEKXLl3k61GeV/keABg6dKhqm/I9RUVFmD9/PrZt24bi4mJYLBZUVlaisLBQddzAwEBUVlZe9rqJPBkDJCJqtosXL8JkMmH37t0wmUyq10JCQhp9r8Fg0Izyqq2t1ezn6+ureg9g67JqiHJ/6T2N7d8U4eHh+PXXX1Xb4uLiAACHDx/GNddcAwDw9/dHv379NO9fv349Zs2ahaVLlyIxMREdOnTAK6+8gry8vMu2vSnX4+wzctymfE96ejrOnz+PFStWoFevXvD390diYiJqampUx71w4QK6du3q5BMh8h6sQSIilb59+8LX11f1EP/1119Vw+pHjBgBi8WC4uJi9OvXT/UVGRkJwDZk3GKxaI7ftWtXnDlzRv756NGjumcr+vTpA19fX+zcuVPeVlpaetmpAkaMGAGz2awKkm655RZ07twZL7/88mXP++233yIpKQmPPvooRowYgX79+uH48eMtv5Ar9O233+KJJ57AxIkTMXjwYPj7++PcuXOa/fbt24cRI0a4oYVEbQcDJCJSCQkJwbRp0zB79mxs2bIF+/btwwMPPKDqAurfvz+mTJmCqVOn4qOPPkJBQQF27NiBrKwsfPbZZwBsXV4XL15ETk4Ozp07JwdBN910E9544w3k5+dj165dePjhhzXZElfr0KED0tPTMXv2bGzduhX79+/HtGnTYDQa5cyLMyNGjEB4eDi+/fZbeVtISAj+93//F5999hkmTZqEzz//HD/99BN++OEHLFmyBADkrFpcXBx27dqFzz//HEeOHMGCBQtUQVpri4uLw3vvvYeDBw8iLy8PU6ZMQWBgoGa/r7/+GrfccosbWkjUdjBAIiKNV155Bddddx1uv/12JCcnY+zYsRg1apRqnzVr1mDq1Kl4+umnMWDAAKSmpmLnzp3o2bMnACApKQkPP/ww0tLS0LVrVzl4WLp0KWJiYnDdddfh97//PWbNmoWgoCDdr2nZsmVITEzEbbfdhuTkZFx77bUYNGhQozVGJpMJGRkZ+Pvf/67aftddd2H79u0ICgrC1KlTMWDAANx0003YsmUL1q9fj9tuuw0A8NBDD+E3v/kN0tLSkJCQgPPnz+PRRx/V9Tob88477+DXX3/FyJEjcf/99+OJJ57QzCeVm5uL0tJS/Pa3v3VTK4naBoNwLAYgIvICFRUV6N69O5YuXYpp06Y1uJ/ZbMbgwYOxZ88e9OrVqxVb6B5paWkYPnw45s2b5+6mELkVM0hE5BXy8/Px/vvv4/jx49izZw+mTJkCALjzzjsbfV9kZCTeeecdzUgvT1RTU4OhQ4fiqaeecndTiNyOGSQi8gr5+fn4f//v/+Hw4cPw8/PDqFGjsGzZMtWweCIiCQMkIiIiIgfsYiMiIiJywACJiIiIyAEDJCIiIiIHDJCIiIiIHDBAIiIiInLAAImIiIjIAQMkIiIiIgcMkIiIiIgc/H+gb3KjKv1+iAAAAABJRU5ErkJggg==", "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.Fe)\n", "\n", "layer_Fe = nx.Layer(id = \"Fe layer\",\n", " material = iron,\n", " thickness = 3000)\n", "\n", "site = nx.Hyperfine(magnetic_field = 33,\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", "exp = nx.Experiment(beam = beam,\n", " objects = [sample],\n", " isotope = nx.lib.moessbauer.Fe57)\n", "\n", "detuning = np.linspace(-200, 200, 1001)\n", "\n", "energy_spectrum = nx.EnergySpectrum(experiment = exp,\n", " detuning = detuning,\n", " electronic = True,\n", " scaling = \"auto\",\n", " background = 0,\n", " resolution = 1,\n", " distribution_points = 1,\n", " fit_weight = 1,\n", " kernel_type = \"Gauss\",\n", " residual = None,\n", " time_gate = [100, np.inf],\n", " id = \"my energy spectrum\")\n", "\n", "intensity = energy_spectrum.Calculate()\n", "\n", "plt.plot(detuning, intensity)\n", "plt.xlabel('detuning (Gamma)')\n", "plt.ylabel('rel. transmission')\n", "plt.show()" ] } ], "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 }