{ "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": "iVBORw0KGgoAAAANSUhEUgAAAkgAAAGwCAYAAABSN5pGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAABgl0lEQVR4nO3deXgUVdo28Lu6sweSAIEsECBAWGRfJAZhUIkGwdE4DkYHZXn5UEcZZSIwgixuI4qCgPoOOqMw4zuIMqPoIEaZAC4Qwy77KhC2JEDIQgJZus/3R6eqq7o6IQld6XT3/buuXEmqq7tOdYWuh+c85xxJCCFARERERAqTuxtARERE1NQwQCIiIiJywACJiIiIyAEDJCIiIiIHDJCIiIiIHDBAIiIiInLAAImIiIjIgZ+7G+CprFYrzp07h+bNm0OSJHc3h4iIiOpACIGSkhLExsbCZKo5T8QAqYHOnTuHuLg4dzeDiIiIGuD06dNo165djY8zQGqg5s2bA7C9wWFhYW5uDREREdVFcXEx4uLilPt4TRggNZDcrRYWFsYAiYiIyMNcrzyGRdpEREREDhggERERETlggERERETkgAESERERkQMGSEREREQOGCAREREROWCAREREROSAARIRERGRAwZIRERERA4YIBERERE5YIBERERE5IABEhEREZEDBkhEVKNKixWVFqu7m0GN5GqFxd1NIGoyGCARkVMWq8DQ1zdgyGsbYLEKdzeHDLZu73n0mJuBD3884e6mEDUJDJCIyKnCsgrkFZfjQkk5Cssq3N0cMtjTH+8CALy09oCbW0LUNDBAIiIiInLAAImInGKnGhH5MgZIRERERA4YIBERESTJ3S0galoYIBGRU0LVx8buNiLyNQyQiMgpqypCUv9MROQLGCARkVPqoIjxkfeTwD42IjUGSETklHpuSGaQiMjXMEAiIqesVnUXmxsbQkTkBgyQiMgpTQ0SIyQi8jEMkIjIKXaxEZEvY4BERE5pR7G5sSFERG7AAImInNLWIDFCIiLfwgCJiJxSZ40EAyTvx1H+RBoMkIjIKXaxEZEvY4BERE5xJm0i8mUMkIjIKavV+c9ERL6AARIROcUMkm9hCRKRFgMkInKKa7ERkS9jgERETjGDRES+jAESETnFmbSJyJcxQCIip7hYrW+RWIREpMEAiYic4kSRROTLGCARkVOcKJKIfBkDJCJyikXavkXiQH8iDQZIROSUpkibKSSvJ8BrTKTGAImInGIXGxH5MgZIROSUdhQbIyQi8i0MkIjIKc6D5FtYg0SkxQCJiJziUiNE5MsYIBGRU+xiIyJfxgCJiJzSdrG5rx3UODiTNpEWAyQicorzIBGRL2OAREROaWuQGCARkW9hgERETnEeJCLyZU0iQHr33XfRsWNHBAUFITExEVu3bq11/9WrV6N79+4ICgpC7969sW7dOs3jQgjMnTsXMTExCA4ORnJyMo4ePer0tcrLy9GvXz9IkoTdu3e76pSIPJ7VqvqZGSSvxxIkIi23B0iffPIJ0tPTMW/ePOzcuRN9+/ZFSkoK8vPzne6/ZcsWPPzww5g0aRJ27dqF1NRUpKamYt++fco+CxYswNKlS7Fs2TJkZ2cjNDQUKSkpuHbtmu71ZsyYgdjYWMPOj8hTMYNERL7M7QHSokWLMHnyZEycOBE33XQTli1bhpCQEHz44YdO91+yZAlGjhyJ6dOno0ePHnj55ZcxYMAAvPPOOwBs2aPFixdj9uzZuO+++9CnTx/84x//wLlz57BmzRrNa3399df49ttv8eabb163neXl5SguLtZ8EXkzddKINUhE5GvcGiBVVFRgx44dSE5OVraZTCYkJycjKyvL6XOysrI0+wNASkqKsv+JEyeQm5ur2Sc8PByJiYma18zLy8PkyZPx0UcfISQk5LptnT9/PsLDw5WvuLi4ep0rkaexcBQbEfkwtwZIFy9ehMViQVRUlGZ7VFQUcnNznT4nNze31v3l77XtI4TAhAkT8MQTT2DQoEF1auvMmTNRVFSkfJ0+fbpOzyPyVJouNmstO5JXkDgREpGGn7sb4A5vv/02SkpKMHPmzDo/JzAwEIGBgQa2iqhpUdcdWZhB8nrsRiXScmsGKTIyEmazGXl5eZrteXl5iI6Odvqc6OjoWveXv9e2z4YNG5CVlYXAwED4+fmhS5cuAIBBgwZh/PjxN35iRF5AvdQIb55E5GvcGiAFBARg4MCByMzMVLZZrVZkZmYiKSnJ6XOSkpI0+wPA+vXrlf3j4+MRHR2t2ae4uBjZ2dnKPkuXLsXPP/+M3bt3Y/fu3co0AZ988gn+/Oc/u/QciTwVR7H5FnaxEWm5vYstPT0d48ePx6BBgzB48GAsXrwYpaWlmDhxIgBg3LhxaNu2LebPnw8AeOaZZzB8+HAsXLgQo0ePxqpVq7B9+3a8//77AGz/yKdOnYpXXnkFCQkJiI+Px5w5cxAbG4vU1FQAQPv27TVtaNasGQCgc+fOaNeuXSOdOVHTpl2LjRESEfkWtwdIaWlpuHDhAubOnYvc3Fz069cPGRkZSpF1Tk4OTCZ7omvIkCFYuXIlZs+ejVmzZiEhIQFr1qxBr169lH1mzJiB0tJSPPbYYygsLMTQoUORkZGBoKCgRj8/Ik8lmEEiIh8mCRYXNEhxcTHCw8NRVFSEsLAwdzeHyOXe++445n99CADw0n09MS6po3sbRIbqNe8bXCmvAgD88uoomEzsciPvVNf7t9sniiSipknTxcYUktdTh0PsUiVigERENWCRtu/i9SZigERENVBnjZhR8C283kQMkIioBlbNWmzuawc1PgZIRAyQiKgGAswg+RL1FeblJmKAREQ1EJp5kNzXDmp8vNxEDJCIqAbqmyQzSN5PO+8VrzcRAyQick5wLTZfwi42Ii0GSETklLpbzWJ1XzuocWiCIgZIRAyQiMg5Fmn7FvX1FoyQiBggEZFz6piIt0vvJ6nm0mZRPhEDJCKqgeYeyQyST2HNGREDJCKqATNIvkXbxUZEDJCIyCnNDZN3TK+nCYh5vYkYIBFRDTQZJN4xvZ12mD+vNxEDJCJyivPi+Bh2qRJpMEAiIqfUWQTeML0fu1SJtBggEZFTVtak+BTBLlUiDQZIROSUtmiXN0xvZ9WsxebGhhA1EQyQiMgpDvv2LSzSJtJigERETjGD5Fs4zJ9IiwESEV0Xb5hE5GsYIBGRUxzF5jscM4RcnJiIARIR1YDzIPkOx+vL603EAImIasBh377D8eryahMxQCKiGnDiQN/h2MXGonwiBkhEVAMrR7H5DMd5j3i1iRggEVENtF1s5M0cu1AZEBMxQCKiGrGLzVewSJtIjwESETnFIm3fxatNxACJiGrAmZV9BzNIRHoMkIjIKa7F5jscM4ScKJKIARIR1YAZJN/BDBKRHgMkInJK1PIbeRf9RJG83kQMkIjIKWaQfId+okg3NYSoCWGAREROqW+arEnxbrqJInm5iRggEZFzXKzWh+hm0uYFJ2KAREROqTNIvF16N/1M2m5qCFETwgCJiJxiBsl36EaxuacZRE0KAyQicoozafsOx6vLmjMiBkhEVANR4y/kbTiKjUiPARIROcUaJN+hv7684kQMkIjIKW0NEm+Y3owzaRPpMUAiIudUN0nHeXLIuzgGwLzeRAyQiKgGVnax+QzdUiNMIRExQCIi57RLjfCG6c04zJ9IjwESETmlHtrPG6Z340SRRHoMkIjIKc1NkjdMr6Yv0uYFJ2KAREROaeMj3jC9ma4GyS2tIGpaGCARkVPaGiT3tYOMx4kiifQYIBFRDVQ1SLxhejV9kTYvOBEDJCJyimux+S4GxEQMkIioBup7JCcO9G6Oi9NysVoiBkhEVAPNRJG8X3o1zoNEpMcAiYic0t40ecv0Zrqry8tNxACJiJzTLlbrtmZQI9CNYmOERMQAiYicE1yLzWc4Xl+r1S3NIGpSGCAR0XVxZmXvxhokIj0GSETklHaYP3k3x4kiecWJGCARkVOCE0X6DGaQiPQYIBGRU+qbJufF8W6O81wxg0TEAImIasB7pO9wHLXGa0/EAImIasAuNt/BLjYivSYRIL377rvo2LEjgoKCkJiYiK1bt9a6/+rVq9G9e3cEBQWhd+/eWLduneZxIQTmzp2LmJgYBAcHIzk5GUePHtXsc++996J9+/YICgpCTEwMHn30UZw7d87l50bkqayaIm3eMr2ZLkDi5SZyf4D0ySefID09HfPmzcPOnTvRt29fpKSkID8/3+n+W7ZswcMPP4xJkyZh165dSE1NRWpqKvbt26fss2DBAixduhTLli1DdnY2QkNDkZKSgmvXrin73H777fj0009x+PBh/Pvf/8bx48fx29/+1vDzJfIY6gCJN0yv5hgAs+aMCJCEm6vxEhMTcfPNN+Odd94BAFitVsTFxeEPf/gDnnvuOd3+aWlpKC0txdq1a5Vtt9xyC/r164dly5ZBCIHY2Fg8++yzmDZtGgCgqKgIUVFRWLFiBR566CGn7fjyyy+RmpqK8vJy+Pv76x4vLy9HeXm58ntxcTHi4uJQVFSEsLCwG3oPiJqiMcu2YNvJywCApE6t8PFjt7i5RWSUfWeLcM/bPyq/L324P+7tG+vGFhEZp7i4GOHh4de9f9c7g5SXl4dHH30UsbGx8PPzg9ls1nzVR0VFBXbs2IHk5GR7g0wmJCcnIysry+lzsrKyNPsDQEpKirL/iRMnkJubq9knPDwciYmJNb5mQUEB/vnPf2LIkCFOgyMAmD9/PsLDw5WvuLi4ep0rkacR7GLzWRzFRgT41fcJEyZMQE5ODubMmYOYmBhIktTgg1+8eBEWiwVRUVGa7VFRUTh06JDT5+Tm5jrdPzc3V3lc3lbTPrI//elPeOedd1BWVoZbbrlFk5VyNHPmTKSnpyu/yxkkIm/Ftdh8B68vkV69A6Qff/wRP/zwA/r162dAcxrX9OnTMWnSJJw6dQovvvgixo0bh7Vr1zoN+gIDAxEYGOiGVhK5h2YtNt5AvZpjzRFrkIgaECDFxcW5LP0aGRkJs9mMvLw8zfa8vDxER0c7fU50dHSt+8vf8/LyEBMTo9nHMaiLjIxEZGQkunbtih49eiAuLg4//fQTkpKSbvTUiDyeJoPELjav5nh1GR8RNaAGafHixXjuuedw8uTJGz54QEAABg4ciMzMTGWb1WpFZmZmjUFKUlKSZn8AWL9+vbJ/fHw8oqOjNfsUFxcjOzu71sDHWr18tboQm8iXCY5i8xmO/+nl9SZqQAYpLS0NZWVl6Ny5M0JCQnRFzQUFBfV6vfT0dIwfPx6DBg3C4MGDsXjxYpSWlmLixIkAgHHjxqFt27aYP38+AOCZZ57B8OHDsXDhQowePRqrVq3C9u3b8f777wMAJEnC1KlT8corryAhIQHx8fGYM2cOYmNjkZqaCgDIzs7Gtm3bMHToULRo0QLHjx/HnDlz0LlzZ2aPiKqJGn4m76PLILmlFURNS70DpMWLF7u0AWlpabhw4QLmzp2L3Nxc9OvXDxkZGUqRdU5ODkwme6JryJAhWLlyJWbPno1Zs2YhISEBa9asQa9evZR9ZsyYgdLSUjz22GMoLCzE0KFDkZGRgaCgIABASEgIPvvsM8ybNw+lpaWIiYnByJEjMXv2bNYZEVXT1iDxlunNHC8va5CImsA8SJ6qrvMoEHmqX7/9I/aeLQIA9G8fgc+fvNXNLSKj7DhVgAf+Yp8GZcEDffDgzRylS96prvfvemeQAMBisWDNmjU4ePAgAKBnz56499576z0PEhE1XVyLzXfo12LjBSeqd4B07NgxjBo1CmfPnkW3bt0A2CZRjIuLw1dffYXOnTu7vJFE1Pi0E0WSN+MoNiK9eo9ie/rpp9G5c2ecPn0aO3fuxM6dO5GTk4P4+Hg8/fTTRrSRiNxAc5PkHdOr6TNIRFTvDNJ3332Hn376CS1btlS2tWrVCq+99hpuvZU1CkTeQn2TtPKO6dU4USSRXr0zSIGBgSgpKdFtv3LlCgICAlzSKCJyP80oNuYUvJoug8TLTVT/AOmee+7BY489huzsbAghIITATz/9hCeeeAL33nuvEW0kIjfjDdO7OQbAvNxEDQiQli5dqkyoGBQUhKCgINx6663o0qULlixZYkQbicgNOJO2D2GVNpFOvWuQIiIi8MUXX+Do0aM4dOgQAKBHjx7o0qWLyxtHRO5j1XSxkTdzvL6sOSNq4DxIAJCQkICEhARXtoWImhDtIDbeMb2ZvgaJ15uoTgFSeno6Xn75ZYSGhiI9Pb3WfRctWuSShhGRe/Em6TtYg0SkV6cAadeuXaisrFR+rokkSa5pFRG5HadB8h0cxUakV6cAaePGjU5/JiIvprpJcl4c78Z5kIj06j2KzVFxcTHWrFmjFGwTkXcQNfxM3ofXl0iv3gHSgw8+iHfeeQcAcPXqVQwaNAgPPvggevfujX//+98ubyARuYdmokhmFLwbu9iIdOodIH3//fcYNmwYAODzzz+HEAKFhYVYunQpXnnlFZc3kIjcgxkk36Ev0uYVJ6p3gFRUVKSsw5aRkYEHHngAISEhGD16NI4ePeryBhKRewhGSD7DMWPEeZCIGhAgxcXFISsrC6WlpcjIyMBdd90FALh8+TKCgoJc3kAicg9OFOk7OIqNSK/eE0VOnToVY8eORbNmzdChQwfcdtttAGxdb71793Z1+4jITbRLjfCO6c10K40wJCaqf4D05JNPYvDgwTh9+jTuvPNOmEy2JFSnTp1Yg0TkpXi79G6OATDjYaIGLjUyaNAgDBo0CABgsViwd+9eDBkyBC1atHBp44jIfbSj2NzYEDKcfq1aXnCietcgTZ06FR988AEAW3A0fPhwDBgwAHFxcdi0aZOr20dEbqK+RXLiQO/GDBKRXr0DpH/961/o27cvAOA///kPTpw4gUOHDuGPf/wjnn/+eZc3kIjcQ1uD5L52kPF0RdruaQZRk1LvAOnixYuIjo4GAKxbtw5jxoxB165d8T//8z/Yu3evyxtIRO7BQl3foe9ic0sziJqUegdIUVFROHDgACwWCzIyMnDnnXcCAMrKymA2m13eQCJyD45i8x36eZB4vYnqXaQ9ceJEPPjgg4iJiYEkSUhOTgYAZGdno3v37i5vIBG5B+eJ9B36mbSJqN4B0gsvvIBevXrh9OnTGDNmDAIDAwEAZrMZzz33nMsbSETuwVFsvkN3fXnBiRo2zP+3v/2tbtv48eNvuDFE1HRoutiYU/Bq+okiiahOAdLSpUvx2GOPISgoCEuXLq1136efftolDSMi99J0sfGO6dU4zJ9Ir04B0ltvvYWxY8ciKCgIb731Vo37SZLEAInIS6hvmly81LuxSJtIr04B0okTJ5z+TETei2W7voNF2kR69R7mT0S+gRNF+g7dRJG83kT1L9IWQuBf//oXNm7ciPz8fFitVs3jn332mcsaR0TuoxnF5sZ2kPH0M2nzihPVO0CaOnUq3nvvPdx+++2IioqCJElGtIuI3ExbpM0bpjfjKH8ivXoHSB999BE+++wzjBo1yoj2EFFTIZz+SF5IP4qNV5yo3jVI4eHh6NSpkxFtIaImxMqJIn0GM0hEevUOkF544QW8+OKLuHr1qhHtIaImgl1sPkRXg0RE9e5ie/DBB/Hxxx+jTZs26NixI/z9/TWP79y502WNIyL3Eexi8xmORdmcB4moAQHS+PHjsWPHDjzyyCMs0ibyYuxi8x2OE4HyehM1IED66quv8M0332Do0KFGtIeImgh2sfkOXl4ivXrXIMXFxSEsLMyIthBRE8J5kHyHbiZtRkxE9Q+QFi5ciBkzZuDkyZMGNIeImgr1PZI1Kd7NsYuNa+8RNaCL7ZFHHkFZWRk6d+6MkJAQXZF2QUGByxpHRO6j7WJzWzOoMTjOg8ScIVH9A6TFixcb0AwiampYpO07WKRNpNegUWxE5P20w/x5x/Rmupm03dQOoqak3jVIO3fuxN69e5Xfv/jiC6SmpmLWrFmoqKhwaeOIyD0cb5isSfFu+gwSLzhRvQOkxx9/HEeOHAEA/PLLL0hLS0NISAhWr16NGTNmuLyBRNT4dKu784bp1bjUCJFevQOkI0eOoF+/fgCA1atXY/jw4Vi5ciVWrFiBf//7365uHxG5geP9kRkk76ZfrNZNDSFqQuodIAkhYLVaAQD//e9/MWrUKAC2+ZEuXrzo2tYRkVs4G9bPLJL30mUMWYVEVP8AadCgQXjllVfw0Ucf4bvvvsPo0aMBACdOnEBUVJTLG0hEjc9ZLMT4yHs5BsS81kQNCJAWL16MnTt3YsqUKXj++efRpUsXAMC//vUvDBkyxOUNJKLG5zSD5IZ2UONglyqRXr2H+ffp00czik32xhtvwGw2u6RRRNT0WIWAGVyc2hvpMkgMh4nqHyDJKioqkJ+fr9Qjydq3b3/DjSIi92IXm2/RXVtea6L6B0hHjhzBpEmTsGXLFs12IQQkSYLFYnFZ44jIPZx1sXE9Nu/FiSKJ9OodIE2cOBF+fn5Yu3YtYmJiIElMuRN5m8a+QW47WYCjeVfwu0RmoB2dLijDur3nMfaWDmgW2OCkf60cY18Gw0QNCJB2796NHTt2oHv37ka0h4iagMbOII1ZlgUA6BbdDAM7tDTsOJ7o4b/+hDOXr+JwbgkWpfUz5Bhci41Ir96j2G666SbOd0Tk5RqzBsmqujtfKOFyRY7OXL4KAFi797xhx3AsymZ8RNSAAOn111/HjBkzsGnTJly6dAnFxcWaLyLyAk7ukEZlkAqvVio/twjxN+QY3qCiynr9nRqIa7ER6dW7iy05ORkAMGLECM12FmkTeY/GnAfpQkm58rPZxJpGt+BEkUQ69Q6QNm7caEQ7iKgJcXZ/FAYlMNQBkoUzFOoE+ZtwrdK47BFgzyCZTRIsVsF5kIjQgABp+PDhRrSDiJoQ5xkko7rY7HVHFqYudCKCA5Bbec3QY8jX1iQBFjCDRATcwESRZWVlyMnJQUWFtqiyT58+N9woInIvZzdIo5I76qyR1dhEiUdqjJlU5EtgkiQAggESERoQIF24cAETJ07E119/7fRx1iAReT45oyBJ9mDJqMJddbaK8++4h1B1sQG8DkRAA0axTZ06FYWFhcjOzkZwcDAyMjLw97//HQkJCfjyyy+NaCMRNTKhySjYGJdBUv3MG7NOY5Sty8Gvufp68yoQNSBA2rBhAxYtWoRBgwbBZDKhQ4cOeOSRR7BgwQLMnz+/QY1499130bFjRwQFBSExMRFbt26tdf/Vq1eje/fuCAoKQu/evbFu3TrN40IIzJ07FzExMQgODkZycjKOHj2qPH7y5ElMmjQJ8fHxCA4ORufOnTFv3jxddyGRr5LjFAm2uhTAuBokTQaJRdpuIb/rcjzMOJWoAQFSaWkp2rRpAwBo0aIFLly4AADo3bs3du7cWe8GfPLJJ0hPT8e8efOwc+dO9O3bFykpKcjPz3e6/5YtW/Dwww9j0qRJ2LVrF1JTU5Gamop9+/Yp+yxYsABLly7FsmXLkJ2djdDQUKSkpODaNVuh46FDh2C1WvHee+9h//79eOutt7Bs2TLMmjWr3u0n8kZy0CJJUJYTaoyJIjmKzT3ka2CfZoHXgajeAVK3bt1w+PBhAEDfvn3x3nvv4ezZs1i2bBliYmLq3YBFixZh8uTJmDhxIm666SYsW7YMISEh+PDDD53uv2TJEowcORLTp09Hjx498PLLL2PAgAF45513ANiyR4sXL8bs2bNx3333oU+fPvjHP/6Bc+fOYc2aNQCAkSNHYvny5bjrrrvQqVMn3HvvvZg2bRo+++yzerefyBvZMwqSPYNkVIAk1D/zxuwO8rtur0FyX1uImop6B0jPPPMMzp+3TXk/b948fP3112jfvj2WLl2KV199tV6vVVFRgR07diiTTwKAyWRCcnIysrKynD4nKytLsz8ApKSkKPufOHECubm5mn3Cw8ORmJhY42sCQFFREVq2rHkNqPLycs4aTj5DzihIACQYW7hr0RRpG3IIug752pqUbCEvBFG9R7E98sgjys8DBw7EqVOncOjQIbRv3x6RkZH1eq2LFy/CYrEgKipKsz0qKgqHDh1y+pzc3Fyn++fm5iqPy9tq2sfRsWPH8Pbbb+PNN9+ssa3z58/Hiy++WPsJEXkZkyTBWp1fMOqWqb4Zs4vNPRyL8nkViOqZQaqsrETnzp1x8OBBZVtISAgGDBhQ7+CoqTh79ixGjhyJMWPGYPLkyTXuN3PmTBQVFSlfp0+fbsRWEjUupUhbshfuGlVArZkHiZkLt1BGsZmMrTcj8iT1CpD8/f2VQmdXiIyMhNlsRl5enmZ7Xl4eoqOjnT4nOjq61v3l73V5zXPnzuH222/HkCFD8P7779fa1sDAQISFhWm+iLyVustFMniguTruYgbJPeR33VR9R2CgStSAGqSnnnoKr7/+Oqqqqm744AEBARg4cCAyMzOVbVarFZmZmUhKSnL6nKSkJM3+ALB+/Xpl//j4eERHR2v2KS4uRnZ2tuY1z549i9tuuw0DBw7E8uXLYTLV+60g8lpKkTbsw/yNumlyFFvdGT1Zp6kxpu0m8hD1rkHatm0bMjMz8e2336J3794IDQ3VPF7fkWDp6ekYP348Bg0ahMGDB2Px4sUoLS3FxIkTAQDjxo1D27ZtlTmWnnnmGQwfPhwLFy7E6NGjsWrVKmzfvl3JAEmShKlTp+KVV15BQkIC4uPjMWfOHMTGxiI1NRWAPTjq0KED3nzzTWWqAgA1Zq6IfIlVNRGSnEEybhSb/YWZuKidEMYsPaLMpG3wlA5EnqTeAVJERAQeeOABlzUgLS0NFy5cwNy5c5Gbm4t+/fohIyNDKbLOycnRZHeGDBmClStXYvbs2Zg1axYSEhKwZs0a9OrVS9lnxowZKC0txWOPPYbCwkIMHToUGRkZCAoKAmDLOB07dgzHjh1Du3btNO3h6A0ibdGuXKTdGKPYOJN27SxCwGRAl6dVVXMGGDcpKJEnqXeAtHz5cpc3YsqUKZgyZYrTxzZt2qTbNmbMGIwZM6bG15MkCS+99BJeeuklp49PmDABEyZMaEhTiXyEaqJIod7ieuxiqzuLVcDfbMQrs0ibyFG9C2/uuOMOFBYW6rYXFxfjjjvucEWbiMjN1Ku7m0zGzo3DiSLrzrjZzG3f5RokXgeiBgRImzZtcrpm2bVr1/DDDz+4pFFE5F7qtdgkh22uphnmzwxSrYzqghTMIBHp1LmLbc+ePcrPBw4c0Ey6aLFYkJGRgbZt27q2dUTkFuq12EwwdvkJzUSRvDHXyrCRhJwokkinzgFSv379IEkSJEly2pUWHByMt99+26WNIyL3sE8UKQHKTNqNsNQIM0i1Mur9UYrylYX3DDkMkUepc4B04sQJCCHQqVMnbN26Fa1bt1YeCwgIQJs2bWA2G1I9SESNTA6GJEAZ2iTXqbiaZqJI9u3UyugsntngOa+IPEmdA6QOHToAsE3kSETeTTPMXxibQeIottqp3xGj3h/5VZUaJEOOQuRZ6j3Mn4i8n3otNlOjThTJW7Ojxnh/7DVnxo5YJPIkXF+DiHQ0a7HJZSmGjWJz/jPZaN4fo0axOc6kbchRiDwLAyQi0lHfII2eG8eqGcXGW7MjdTbHqB5I+RrIXWzs6SRigERETgjVMH9lm0HHsnIUW60aY5Sf/Kr2UWy8DkSsQSIiHfW8OEIydi02ZpBqpw6KjLoGQulSrf7dkKMQeRaXZZDGjx/PpUaIvIZqLTaDi7TVNTYcXq6nmQbB4HmQlBokXgYi12WQ2rZtC5OJPXZE3kCTQZInijQ4ewGwi80ZayPWIJlMXIuNSOayAOnVV1911UsRkZup12KDwSObLJp5kAw6iAezNEoXm+07M0hEdqxBIiIdZ0XaRmV31C/LzIWeaIT3R8kYVncC8CoQ1TFASk9Pr/MLLlq0qMGNIaKmQb5hSpKk3KEbYxQbZ9LWszTK+2Of9wrgRJFEQB0DpF27dtXpxST1fzeJyGOp12KTGnEeJGaQ9LQzaRt1DNt3ZakRXgaiugVIGzduNLodRNSEqNdisxdpG3Osxqix8VRCCM37btwoNnmxWrnejNeBqMHDzo4dO4ZvvvkGV69eBcCULJE30azFZnDhbmMEAJ7K8e0wugZJYpE2kaLeAdKlS5cwYsQIdO3aFaNGjcL58+cBAJMmTcKzzz7r8gYSUeNTuthU3eZG3Zw5iq1mjgGj0d2cZhZpEynqHSD98Y9/hL+/P3JychASEqJsT0tLQ0ZGhksbR0TuoWQUoMogGXYs41er91SOAZHRCTYz50EiUtR7mP+3336Lb775Bu3atdNsT0hIwKlTp1zWMCJyH6fD/LnUSKNzfM+N6oJUJoqUuNYIkazeAVJpaakmcyQrKChAYGCgSxpFRO6lLtK2bzTmWI2xlIanaqwaJMfrzatA1IAutmHDhuEf//iH8rskSbBarViwYAFuv/12lzaOiNxDqNdiq46RGqMGiV07WroaJINqtOw1SJwHiUhW7wzSggULMGLECGzfvh0VFRWYMWMG9u/fj4KCAmzevNmINhJRI5NvxOoibePm4OFEkTVxDFSYQSJqPPXOIPXq1QtHjhzB0KFDcd9996G0tBS/+c1vsGvXLnTu3NmINhJRI5NvkFL1F9BYE0UacgiP5RgwGlWjpazFVn1HYCaPqJ4ZpMrKSowcORLLli3D888/b1SbiMjNhFK0q9pm0LE0XWyMkDR0NUhGTRSpW2rEkMMQeZR6BUj+/v7Ys2ePUW0hoibCceJAwLi6FE2RNu/MGo01zN++WC0DJCJZvbvYHnnkEXzwwQdGtIWImgz7WmxyFsmwGiQra5Bq0ljD/B2XGiGiBhRpV1VV4cMPP8R///tfDBw4EKGhoZrHFy1a5LLGEZF7KBkFzUzaRh2Lo9hq4hgQGZ3FM3GiSCJFvQOkffv2YcCAAQCAI0eOaB6T+L8PIq+g3B8le5G2UQuYWlQva9Qwdk/lGKcYVqRd/d3obCGRJ6l3gLRx40Yj2kFETYi9aNe+zagMkjorwhokLf1abMYcx7GLzahgmMiT1DtAIiLvZ1+LTVImijSqe4ej2GqmK9I2rAbJ9p1F2kR29S7SJiLvp16LTTK6SJuj2GqkH8Vm7FxU9sVqDTkMkUdhBomIdNQzK9vXLzXo5qxZasSQQ3gsx/fDuFFstu/2LlVeCCJmkIhIR70Wm8zodcBsx+CNWU0/is2Y48jXgBNFEtkxg0REOvINUpIk1Sg2Y6i71TgPkpZuHiSDIxdlsVpDj0LkGZhBIiIde5G2PYtk9EKpRh7DUzlm7RqvBonXgYgZJCLSUa/FJtmLkAxh4UzaNWrsUWwSu9iIFMwgEZGOsy42o7MXRh7DUzl2qRk9m7kyDxKvAxEDJCLSE6q12OSsgmE3Z45iq5FjoGLYKLbq7yZjk4VEHoUBEhHp2DNI6pumseuAAexic2RppBok+WXNnCiSSMEAiYh0LKph3/LQb6PqXziKrWaNNVGk/L77m02GHofIk7BIm4h0rKqMghwgGTdJIWuQauIYlBoVP8rX1s/MUWxEMmaQiEhHvjGbJElZn8vSCKPYeGPWaqyZtOX3XckgGTQpKJEnYYBERDryjdhkkmCW50Ey7OasPq4hh/BYjgGj0QsGB1QHSFwTj4gBEhE5YR/2bV/h3bBh/swg1cgxUDEqgJQvAbvYiOwYIBGRjnptLnluHKOyClYWadfIMWNk9FxUfiZT9XE5FxIRAyQi0pEzFSZT445iY+ZCq7GG+dtHsUm6bUS+igESEenYM0j2Ljaju3cA44IwT9VYw/ytDsP8bccy5FBEHoMBEhHpyDdMs0mCfM80rItNvRYbM0gajgGjUUGqRTgLkHgtyLcxQCIiHYuTGiTDRlCpu9g4ik3DMYtjdA0Su9iI7BggEZGOfG/UzINkwA1TCKFZ1oIZJC3dYrVGTbVQHZj6MYNEpGCAREQ66i42k4Gj2BorQ+KpdIvVGpzF8zPZM0jM5pGvY4BERDpyoCJJ9gVMjche6CdC5PByNcesnVFvjTJRpB8zSEQyBkhEpCNnFMySei02A47jJOhi7YtdYyw1og581RkkdneSr2OAREQ6zkaxGZFRcPaSvDHb6RerNTaLZ+tSdX5sIl/DAImIdOR7o6QaxWbEzdlZMMTaF7vGmAdJfQ00E4MyPiIfxwCJiHQsSgbJFiSpt7mSsxs+a1/s9KPYXH8M9WtqRi3yOpCPY4BERDpWVQ2S2cDFap114/DGbKerQTI4g2RWZwyZQiIfxwCJiHTso9jsAZIxGSQn23hjVjRGDZL6uppMsNcgMVAlH8cAiYh05BFrZoNrUjiKrXa6GiSDJuuUmQ2eGJTIkzBAIiId+UZskmAfxWbgzdmsnqCQ92WFY5BidJCqDYh5Ici3uT1Aevfdd9GxY0cEBQUhMTERW7durXX/1atXo3v37ggKCkLv3r2xbt06zeNCCMydOxcxMTEIDg5GcnIyjh49qtnnz3/+M4YMGYKQkBBERES4+pSIPJ58czQZPJO2fc03GFrr5Kkc3wojr4EkabtUGaiSr3NrgPTJJ58gPT0d8+bNw86dO9G3b1+kpKQgPz/f6f5btmzBww8/jEmTJmHXrl1ITU1Famoq9u3bp+yzYMECLF26FMuWLUN2djZCQ0ORkpKCa9euKftUVFRgzJgx+P3vf2/4ORJ5IucTRRpXg6ReFJddO3aOAZERs4zLo9jk62zk9SbyJG4NkBYtWoTJkydj4sSJuOmmm7Bs2TKEhITgww8/dLr/kiVLMHLkSEyfPh09evTAyy+/jAEDBuCdd94BYPvwWLx4MWbPno377rsPffr0wT/+8Q+cO3cOa9asUV7nxRdfxB//+Ef07t27MU6TyOPYu9jsGQUjEjvq45iqP414Y7ZzzKYZOdWCWQmQjDsWkSdxW4BUUVGBHTt2IDk52d4YkwnJycnIyspy+pysrCzN/gCQkpKi7H/ixAnk5uZq9gkPD0diYmKNr1lX5eXlKC4u1nwReSsls2MytmjXqqpBYu2Lnn4Um+uPIV9XOUA1MiAm8iRuC5AuXrwIi8WCqKgozfaoqCjk5uY6fU5ubm6t+8vf6/OadTV//nyEh4crX3FxcTf0ekRNmb2LzZ5ZMKT+xapaFJczOOvI74VfIywYbHbsYmOERD7O7UXanmLmzJkoKipSvk6fPu3uJhEZRqiKtI0cxSa/pNngTJWnkt8LP7OBy71Y7dfa9h2GHYvIk7gtQIqMjITZbEZeXp5me15eHqKjo50+Jzo6utb95e/1ec26CgwMRFhYmOaLyFtZVLVBkoEZBWW0nMEzdnsqOVD1r45aLEbUgamuAQDOpE1UzW0BUkBAAAYOHIjMzExlm9VqRWZmJpKSkpw+JykpSbM/AKxfv17ZPz4+HtHR0Zp9iouLkZ2dXeNrEpGeRTWyyciuL/XNmaOn9OSgVM4gGTKKTZXFAziKjUjm586Dp6enY/z48Rg0aBAGDx6MxYsXo7S0FBMnTgQAjBs3Dm3btsX8+fMBAM888wyGDx+OhQsXYvTo0Vi1ahW2b9+O999/H4BtDo+pU6filVdeQUJCAuLj4zFnzhzExsYiNTVVOW5OTg4KCgqQk5MDi8WC3bt3AwC6dOmCZs2aNep7QNQU2SdwVM1PZMAN056p0m8jVQ1SdT+nEe+NOlsI2LvaeBnI17k1QEpLS8OFCxcwd+5c5Obmol+/fsjIyFCKrHNycmAy2ZNcQ4YMwcqVKzF79mzMmjULCQkJWLNmDXr16qXsM2PGDJSWluKxxx5DYWEhhg4dioyMDAQFBSn7zJ07F3//+9+V3/v37w8A2LhxI2677TaDz5qo6bOoMzsG1gYJh+yFehvZg1J/A7sf5esq15qZOZqQCICbAyQAmDJlCqZMmeL0sU2bNum2jRkzBmPGjKnx9SRJwksvvYSXXnqpxn1WrFiBFStW1LepRD5DnVVojFFscvbCqON4KqvSxWaLXuRJHY04hnyd5UvBTB75OrcHSETU9KgzO8aOYrPPwSOBtS+OquQMUiOOYmOxPJENAyQi0lHPTyQZ2OXirEibN2Y7ixIgyaPYjF3uRf2d14F8HQMkItJRJoo0qbvYDDhOdZeRWZIAdu3oyO9FoJ/cxWbsbOYAVDVnLj8UkUdhgEREOkJVl9Ioo9hMcgcbMxdq8vsT4GdcBslxJGF1bx6vA/k8BkhEpGPvYjN2FJsyU7RqFJsRhcieqkrJIJltvxuQxrNaHTJInCiSCACXGiEiJ+T7sLqLzYiMQlV1NORnlrgGmBNWhy42Q4JUh5m0lYCY14F8HDNIRKSjnijSZGCXi30OHhMk2H5m5sJOySD5mzS/u5JjkTYXDSayYQaJiHQaq4ut0mLvYjMyEPNUFscuNgP6H3VdbAZO60DkSZhBIiIdZX0uSVLqg4yuQbIKzoPkSA6IlC42A2qQqmqqQWKgSj6OGSQi0rFa7XUp8hw8lYbcnO01SFziQk8eai8HSEZ0sVVVHySg+jpzsVoiGwZIRKRjUc1wLa8kX2nAxDjqGiTOv6NnqQ4gAwwMkCos9iAV4EzaRDIGSESko84gBSgZJNdHLlWqGiQj13zzVPph/q6/BnJm0F+XQXL5oYg8CgMkItKxrwFmUm6cRszBo65/kYuDBQMkhZzFMXKYvxx0yddZXvfNwgmpyMcxQCIinUrlpikpXS8VhnSx2Y/D2hc9OSg1cpi/+loDgF91oFRhxNoyRB6EARIR6diH35uM7WJT1SCZDRwt56n0w/yNm2pBySBVXwcjuvOIPAkDJCLSqVKKgyUlo2AVrg9enNUgsTjYzl6DZO9ic3UXZKVDkbbSpcpAlXwcAyQi0qmsqr5pmkxK1wvg+iyStgaJxcGOlBokf/tHtauD1EqHYf5Kl2oVLwT5NgZIRKRTKRdp+9mLtAHXB0jaGiTbNmaQ7JQapOouNsD1mR1dF5uSQWKARL6NARIR6SiFuyZJEyC5eiSbOoPE+Xf0LA5dbIARAZJjF5s87xWvA/k2BkhEpGGrc7H97Ge2FU/L2R3XZ5DsxeAcxaYnzwkVoAqQXL3ciBxw2bvYjCvKJ/IkDJCISEN9Y9QP/XbtTVPOUqgzSAyQ7CwOwQvg+q4vudZIV6TNDBL5OAZIRKShDZBsHxEBBt005RokP5N9HiT2sNnJwZC/n0nJ4hnVxeY4zJ8ZJPJ1DJCISEMdBDnOrmzUKDY/9USRjJAUcneaWbJPt+DqAKnKsUjbz7jFiYk8CQMkItKorM5aSJJ94VKjutgsmokitdvIHiyaTRL85C5IFwcuupm05YkiOYqNfBwDJCLScBz2DRjXxVapnihSHsXGAElhUWXYzAYFLsqUDnJ3qh+LtIkABkhE5KBKNcRf5mdQF5tcg2Q2sYvNGWUaBElSZXZcHKQqRdrVo9hM7GIjAhggEZED+7w49o8Hf4O62KqU7IU9QDIqgWS1Cry89gA+23nGJa93pbwKf/rXHmw+dtElr+eMRTVPlJ9hWTx5Jm25O5VF2kQAAyQicuCsi82ood/aGiRju9i2HL+ED348gfRPf3bJ6735zWF8sv00xv4t2yWv54x6nig/g6ZB0HWxcZg/EQAGSETkwD6qyd7F5m/Q+lzqxWqN7mIrrahSfr5aYbnh19tzpvCGX+N65AybyQTjapAcu9iYQSICwACJiBxUOMyLAwBB1WuBXau68cBCrUpVgyQfzqgMkqT6+VzR1Rt+vZJrVdff6QZZnWSQXF2DdLXSdk2D/LRrsTFAIl/HAImINKoc1uYCgOAAW4DkisyLWqUqW2UyeCZtdQbpfOG1G349owMkIYQmg6QELi7O4l2rDpBCAvwA2LOFrg7EiDwNAyQi0lAKp032j4dg/+oMUqVrA6Ty6oxUkL9ZqX0xKnNxRRXQXLxSfsOv5+psmiP1KLJAPzMC/W3vT7mLA6Sy6qA3OICj2IjUGCARkYbSxeZnzyCFVGeQylycQbpWaTtWoJ9ZWbHe1QGA7Eq5ve1XXRDomSX7+2NEUFeuCsAC/UxKN2e5iwMz+b0I9pczSOxiIwIYIBGRA7kQW1ODJHexGZRBCvQ3IbA6AHB1IbjsSnml8rMrugpNqnmiSstd392mfh8CzCYlgyQHla5yTckg2d7/ACVQNTZDRtTUMUAiIg17TYpZ2Rbib0yAJN/sg/zMqhuz8V1srjiPKlWGxYh6pHIlULXVZxmRQRJCoMzhesvfXR2IEXkaBkhEpKHUpFR3uQDGFWlrM0iN18VW7oIASR1AXDEwgyRn1ozIIFVahFIUH1QdBMv1Zq6+1kSehgESEWlcdehyAew3T1ffNJ1nkIy5MV+tdF0GSQiheQ1XZ9YAe6Aovy9GZJDU11MOjIJV3amCy76QD2OAREQa9qJd+8dDiFE1SJX6GiSjMkjqYOBGz8OxjdcMyLYo2bXqAMmIDJL8PviZJHsgVh0oWayCI9nIpzFAIiINJYPkb88gGdXtcq060Ajyt2eQjCrSVgcWVytu7BiO74MRQ/7tXWzVAZIBGaSy6rmh1NdaXXvGbjbyZQyQiEhDySAF6GuQXDnM32oVmiDA6BokddboRudzcsxA3WjA5YxjF5sRGST5egapgiJ/s33WbiO6Dok8BQMkItKwd7HZb5rNg2zBUolqqPyNqlCNAtNmkIy5KV9zYc2QLkAyIJBwLNI2ogZJHn0nX19ZsEGjFok8CQMkItK45jCzMgCEBwcAAIquui5AUgcsnpZBcny+MUXattc0MoMkX8+IYH/NdqNGLRJ5EgZIRKRR5qQGKbz6BlpY5roASQ6EzCYJ/mZTI9QguS6D5BgguWLaAEflDjVIQQYUsRdXB0jhNQVIlcYvyEvUVDFAIiINZXV3JwFSybUqly0mKwcZ+iLkRhjFdoOZEceaIyMyLTXXILnuWEU1BUhKUT4niyTfxQCJiDTkkU2hgfa6FPUNtOSaa7JIcv1Ls+rjBDbiKDaXF2kbWoNke1/k9+mKC2ftrilAkkeyGTEBJpGnYIBERBpyN5r6phngZ1Jumq6qQyq+pj1OoGqiSFdPUFhlsWqKwj2hSNueQbK97/L75Mo6sMKrFbbXDgnQbJePVezCYxF5GgZIRKRRU1ZBLuQtKK1wyXHkm29Y9evKdS9W4fputmtVru0Sc5wY0pXdXrKy6uxNaPX7Ir9PxS7K4AHOg2EAiAhxfVE+kadhgEREGoXyyKYQ7U2zdVgQACC/pNwlxym+agsAwqqHmIcG+KF6+h2XZy70EzveWADmODGkETVIciAkB0ZhQa7P6uQX265lVFigZrtSlH/VNcEwkSdigEREimuVFqX2JcKh2yWque0mml98zSXHcuxiM5kkpc6m2IV1NoC+lqaiynpDxeaOAZERtTrKHEXV74lSKF9eBauLCuVzq69lVHXwKzOiO4/I0zBAIiKF3OViNklK144sOtx2E80rdlUGSZshUf/sym4kwF7Y3EKVFbuRbjG55kjOspW4OKAD9BkkeTJHIWxB0o0SQiBPDpCaOw+QXDmtA5GnYYBERAqlaDfYH5IkaR6Tswzni1yTQZKzE3LXEQA0DzIm4JBH3rVqZu9KupHCavm5ratfz4gAyXGW6yB/s1LI7oputqKrlUqtVxuHLjY58GOARL6MARIRKS5U1xe1bhaoe6xdi2AAwOmCMpccS65limxm78qT65FcXYMkZ1yaB/khqHo+oRupG5KLtFtXdzsa0cVW7CSAjKy+Lq6oAzt1yXYdWzcP1Mx5pT7OxSuuyRYSeSIGSESkkLNDUeFBusfiI0MBAL9cLHXpsWIigpVtxmWQ5ADJX5kE0RVdbHKA5Kq5odScrZMWG2G7LrkuyOKdvGS7jvJ1VZO7U12VLSTyRAyQiEgh33hjwvQBUsfqG+nFK+UuKd5VjqUKxuTal8tlrh09daU6gGke6KcKkBo+kk1+rpxpu1Je5fK5m+T3IFxVNxUdbgsmzxddveHXP37BFiB1qiVAKrpayfXYyGcxQCIihTKqyUkGKSzIH22rsz17zxTd0HGqLFbkl9iOFa06VnS4LeDIc9FIOZk6GxNUXXwuzxjeEGUOXWyVFuHSuZuuVVpwubr+J1oVrMa4MLNz4FwxACAhqrnuseaBfsrEoK4Ixog8EQMkIlKcqu52keuNHA3q2AIAsP1UwQ0d5+SlMliFbUmLyFB7vVNMdYbkXKFrAyR5RFizQD+lpudGsmByfVBMRDAkA+ZukrNrwf5mzSSO7VuGAAB+uXDlhl5fCIG9ZwsBAH3bheselyRJ+Rs45aKaMyJPwwCJiBTH8m033i5tmjl9fFAHW4C049TlGzrO4dwSALbshclkHy2n1NgUuzZrkadMiBikDPW/kRFacvdXZGiA0s3mqukPAFV9VniQZjRht2hbtkd+/xoqp6AMecXl8DdLuCk2zOk+8t/A8fwbC8aIPBUDJCICYMuoyDf5mgKkgR1aAgB2nrp8Q0XO+8/Zuui6RWmPE1vdhZdzqcylNT1yRiY6PAgtQm2j5gpuoM5JDpAiQgJU3V6uC+pOX7ZlbWIitF2dXau7w84VXbuhJV++P3IBANC/fQuEBPg53adLa9u1OZrHAIl8EwMkIgIA7DlTCMDWvaYeWq7WPbo5YsKDUFphwX8P5jX4WNknbF10gzq21GyPjwyFn0lC8bUql46gkmurosOD0KJ6hvCGFoILIZT6oJahAcr8UK6smzp43lYf1C1Km90JD/ZH1+qgcuuJSw1+/bV7zgMA7ujepsZ9boq1db39XP13QeRrGCAREQBg20lbt5ncjeaMySThgQHtAACfbDvdoONcvFKO3acLAQBJnVppHgv0MyvZq/3VRcQ3ymIV9gxSWBBaVmeQLjcwA1N8rUpZpiQixF/JIJ0pdF0GST73HjH6AuohnSMBAOv25jbwtYuQfaIAJgm4p09MjfsNrP47OJxX0uD3isiTMUAiIgBAZnVGSL4B12TMoHYwScAPRy9iy7GL9T7Ov3ecgcUq0KddOOKqi47V+rePAABsbsBrO3M0vwQVFitCA8yIjQhWMkiXrjTspn+pevLEkAAzgvzNyiiwAy4K6ErLq7A7pxCArQvMkRygZuzLVdpSV0IIvLX+KADg131j0a6F/v2XtW4eiB4xYRAC+Hpfw4IxIk/GAImIcPB8MfafK4bZJGFEj5q7XQCgQ6tQPHJLBwDAjH/vUYbr18WFknK89/0vAKC8hqPhXW3H/3Z/7g0tKCvbVR1s9G4XDrNJQodWtqDgRAMnvJSfJ48o6xcXAQD4+XQhqiw3PtR/w6F8VFisaNciGJ1b6+co6t0uHH3ahaPCYsWi9Ufq9dr/99Mp/PdgHvxMEp68rct190/tFwsA+HzXmXodh8gbMEAi8nFCCCz89jAAIKVnlGa9spo8e2c3dGgVgjOXr+L+d7fgy5/Poeg6o8LOXC7DhOVbUVBage7RzXF//7ZO97utW2u0CPHHuaJrWL29Yd14amv3nANgz4wlVHfhnbxU2qBC8yPVRctywXT36OZoGRqA4mtV+P7ohRtqq9Uq8NcfbAHkb/q31a2HJ5uR0h0A8M/sHPxl0/FaA0khBI7kleCFL/djzhf7AQDpd3VVRsTV5t5+sTBJtu7XTYfz63s6RB7N+fCFRvbuu+/ijTfeQG5uLvr27Yu3334bgwcPrnH/1atXY86cOTh58iQSEhLw+uuvY9SoUcrjQgjMmzcPf/3rX1FYWIhbb70Vf/nLX5CQkKDsU1BQgD/84Q/4z3/+A5PJhAceeABLlixBs2bOR+8QeSOLVeCNbw7jvwfz4WeS8MyIrnV6XniIP/7xP4Mx7sOtOHWpDE9/vAuArSYnItgfESEBaB7kh9AA29pnucXXsOPUZVRaBFqFBuAvjwyEv9n5/8+C/M148rYu+PO6g3h13UH0aRdR41D06/lqz3lsPnYJZpOkBGStmweiZWgACkorsPt0IW5xqIO6HnmKAznA8DObkNqvLT7cfAKvfX0IAzu01MxdVFdCCLy67iD2nClCs0A/jBvSscZ9hyZEYsrtXfDOxmN4PeMQPso6iT7tIhASYEZJeRWKr1ai6GolLpVW4HJpBapUAdTkYfF44led69SmmPBgTBgSjw83n8C01XuwYuLN6NVWP28SkTeShKvnx6+nTz75BOPGjcOyZcuQmJiIxYsXY/Xq1Th8+DDatNGn+rds2YJf/epXmD9/Pu655x6sXLkSr7/+Onbu3IlevXoBAF5//XXMnz8ff//73xEfH485c+Zg7969OHDgAIKCbAWVd999N86fP4/33nsPlZWVmDhxIm6++WasXLmyTu0uLi5GeHg4ioqKEBbWsA9vosZmtQrkFJRh/7li7D1bhLV7zuHMZVtx8SupvWrs9qrJlfIq/PX7X/DF7rM4een6EwomxrfEaw/0cbr+l1qVxYoH38vCzpxCBPubMeHWjhjdOwYJUc0Q6Geu9bkWq8Cx/Cv4984z+PDHE6iyCjz2q06YNaqHss+01T/jXzvOYFTvaLz7uwE1Zmoc5VwqQ/Ki71BhsSJj6jB0j7b92794pRwjF/+Ai1fK0bl1KJ4ekYDburbRLBPijBACBaUV2HayAH/fcgpZv9hGpi0c0xcPDGx33ef+MzsHb357+LpzOgX6mZDUuRUmDY3HsITWdTpXWWl5FR74yxYcyi2BSbKNfPtV19boGRuG7tFhCA1sEv/PJqqzut6/3R4gJSYm4uabb8Y777wDALBarYiLi8Mf/vAHPPfcc7r909LSUFpairVr1yrbbrnlFvTr1w/Lli2DEAKxsbF49tlnMW3aNABAUVERoqKisGLFCjz00EM4ePAgbrrpJmzbtg2DBg0CAGRkZGDUqFE4c+YMYmNjr9tuowKkjH25OF90Ff5mEwLMJvj7SfA3m2y/+1VvM5vgb7ZtV3+u13Yl1Y8JiFoec3yeqOUx3VFqfEz9q+6xOh7Dsd2otd01P6+u56t7XHf4Ws7XJe9pfdoNWISAxWqFxQpYrFZUWgQqLVaUXKvCmctlOHmxDAfPFysr28vCgvww99c98dvr3JCvp6isEvkl11B4tRKXSytwpbwKpRUWXK2oQmSzQPRuG+50WYsaX+9qJZ785w5sPqYdzh5gNsFskuBnkiBJtpF1JkmCELYrUlZuQYWqFug3/dvijTF9YVZNSLn3TBHuffdHCAEM79oaqf1j0bFVKEID/RDkZ4ZViOr30/Z16UoF9p4twvLNJ5BfUo6kTq2wcnKiJrA6cK4Y4z7ciouqwunmgX7wM0uQqttnFYBVCFitAlXVX+rusUA/E16+rxcevDmuzu9TWUUVsk8U4ExBGcoqLGge5I/mQX4ID/ZHy9AAtGoWgFahgQjwa3hFRVFZJWb8+2d8s187tYMkAR1bhaJTZCjiWoagbUQwggPMCPQzwc9suy7ytTJJEvz9TAis/izzr76O8lsoQVJeU/4uQar+DtVnnXY/QP1vw9mHoH1H+7HU5yDp9lS/tmO7XOV6d17d5129n3+941/n9a97bNvfs5D/poVQ/Yzq36t/tsr7qv4N1LJ/RZUVpRVVKCu3YHSfGGV+NFep6/3braF/RUUFduzYgZkzZyrbTCYTkpOTkZWV5fQ5WVlZSE9P12xLSUnBmjVrAAAnTpxAbm4ukpOTlcfDw8ORmJiIrKwsPPTQQ8jKykJERIQSHAFAcnIyTCYTsrOzcf/99+uOW15ejvJy+wdfcbFrRqw4Wrk1R5nEjcgIAX4mdI9ujp6xYRgc3xIje8YgOKD2rExdhIf4XzdjUq/XC/bH/01KxDf7c/GvHWeR/csllJRX2YKf65QOBfiZMKxLJB4e3B4jerTRZYh6twvH/Pt7Y/aaffjuyAV8V49/c50iQ/Hmg311r3lTbBgy04fjg80nsPbnc/jlYqkuGHVGkoAOLUOQ0jMaj9zSwenIvtqEBPjh9m61F9bfqPAQf7z36CAczi3BN/tzsSvnMg6cL0ZecTlOXCxtcME70fV0j2nu8gCprtwaIF28eBEWiwVRUVGa7VFRUTh06JDT5+Tm5jrdPzc3V3lc3lbbPo7dd35+fmjZsqWyj6P58+fjxRdfrOOZNdwtnWz1C5VVVlRarKiwWFFR/bOcEaiwVD/mZHFMycn/luyPOfxey3+JdM918j+quhzH8RhSjb/U3j79Y47Prfm8a3rN671ubcfQ7as/UI3HafD7pTum9lGzSVL+1+5X/bOf2YTQADPatQhBXMtg9IgJQ+fWzWqs/2lqJEnCyF4xGNkrRpmg8VqlBVUWofwvVP7fJ2B7v4IDzIgJD9ZkjJx5aHB7DOjQAp9uO42fzxTiXOE1lFVU4VqlFWaTBJNke0/NJgnhwf6Ij2yG27u3xm8Htquxmy88xB/pd3ZF+p1dcaW8CrlF15T/FZtU2S4JsF0fkwkRIf4I8r/xALUxdIturinuvnilHIfOl+BUQSlyLpXhfNE1lFdZUF5lVbJv8pctY2b73LJ9ptl+B+zZEDlrIWcohLBtk7Mdovox9e/2rE9tmSXtc+wb9T86y/Da26fNGtf6WaP/VLA/Vuvzanmsts/sWp5X24MNPp4EmCTbvxOp+rvtd0nzmPZ31c8m+TH1c23vm7+fhNAAP4QG+inzlrkDO4/raObMmZrMVXFxMeLi6p4Gr6u6DL0l8kWSJLn8w7JrVHPMvucml76mrFmgX41LtniLyGaBGJoQiKGofe4sIk/k1v9GRkZGwmw2Iy9P26+dl5eH6Ohop8+Jjo6udX/5+/X2yc/XDlmtqqpCQUFBjccNDAxEWFiY5ouIiIi8k1sDpICAAAwcOBCZmZnKNqvViszMTCQlJTl9TlJSkmZ/AFi/fr2yf3x8PKKjozX7FBcXIzs7W9knKSkJhYWF2LFjh7LPhg0bYLVakZiY6LLzIyIiIs/k9i629PR0jB8/HoMGDcLgwYOxePFilJaWYuLEiQCAcePGoW3btpg/fz4A4JlnnsHw4cOxcOFCjB49GqtWrcL27dvx/vvvA7Cl4adOnYpXXnkFCQkJyjD/2NhYpKamAgB69OiBkSNHYvLkyVi2bBkqKysxZcoUPPTQQ3UawUZERETeze0BUlpaGi5cuIC5c+ciNzcX/fr1Q0ZGhlJknZOTA5PJnugaMmQIVq5cidmzZ2PWrFlISEjAmjVrlDmQAGDGjBkoLS3FY489hsLCQgwdOhQZGRnKHEgA8M9//hNTpkzBiBEjlIkily5d2ngnTkRERE2W2+dB8lScKJKIiMjz1PX+7RljfYmIiIgaEQMkIiIiIgcMkIiIiIgcMEAiIiIicsAAiYiIiMgBAyQiIiIiBwyQiIiIiBwwQCIiIiJywACJiIiIyIHblxrxVPIE5MXFxW5uCREREdWVfN++3kIiDJAaqKSkBAAQFxfn5pYQERFRfZWUlCA8PLzGx7kWWwNZrVacO3cOzZs3hyRJLnvd4uJixMXF4fTp0167xpu3nyPPz/N5+zl6+/kB3n+OPL+GE0KgpKQEsbGxMJlqrjRiBqmBTCYT2rVrZ9jrh4WFeeUfvZq3nyPPz/N5+zl6+/kB3n+OPL+GqS1zJGORNhEREZEDBkhEREREDhggNTGBgYGYN28eAgMD3d0Uw3j7OfL8PJ+3n6O3nx/g/efI8zMei7SJiIiIHDCDREREROSAARIRERGRAwZIRERERA4YIBERERE5YIDkJidPnsSkSZMQHx+P4OBgdO7cGfPmzUNFRYVmvz179mDYsGEICgpCXFwcFixYoHut1atXo3v37ggKCkLv3r2xbt26xjqNWv35z3/GkCFDEBISgoiICKf7SJKk+1q1apVmn02bNmHAgAEIDAxEly5dsGLFCuMbX0d1OcecnByMHj0aISEhaNOmDaZPn46qqirNPk35HNU6duyou16vvfaaZp+6/M02Ze+++y46duyIoKAgJCYmYuvWre5uUoO88MILumvVvXt35fFr167hqaeeQqtWrdCsWTM88MADyMvLc2OLr+/777/Hr3/9a8TGxkKSJKxZs0bzuBACc+fORUxMDIKDg5GcnIyjR49q9ikoKMDYsWMRFhaGiIgITJo0CVeuXGnEs6jZ9c5vwoQJums6cuRIzT5N+fzmz5+Pm2++Gc2bN0ebNm2QmpqKw4cPa/apy99lXT5TXUKQW3z99ddiwoQJ4ptvvhHHjx8XX3zxhWjTpo149tlnlX2KiopEVFSUGDt2rNi3b5/4+OOPRXBwsHjvvfeUfTZv3izMZrNYsGCBOHDggJg9e7bw9/cXe/fudcdpacydO1csWrRIpKeni/DwcKf7ABDLly8X58+fV76uXr2qPP7LL7+IkJAQkZ6eLg4cOCDefvttYTabRUZGRiOdRe2ud45VVVWiV69eIjk5WezatUusW7dOREZGipkzZyr7NPVzVOvQoYN46aWXNNfrypUryuN1+ZttylatWiUCAgLEhx9+KPbv3y8mT54sIiIiRF5enrubVm/z5s0TPXv21FyrCxcuKI8/8cQTIi4uTmRmZort27eLW265RQwZMsSNLb6+devWieeff1589tlnAoD4/PPPNY+/9tprIjw8XKxZs0b8/PPP4t577xXx8fGaz5SRI0eKvn37ip9++kn88MMPokuXLuLhhx9u5DNx7nrnN378eDFy5EjNNS0oKNDs05TPLyUlRSxfvlzs27dP7N69W4waNUq0b99e8xlyvb/LunymugoDpCZkwYIFIj4+Xvn9f//3f0WLFi1EeXm5su1Pf/qT6Natm/L7gw8+KEaPHq15ncTERPH4448b3+A6Wr58ea0BkuOHgNqMGTNEz549NdvS0tJESkqKC1t442o6x3Xr1gmTySRyc3OVbX/5y19EWFiYcl095RyFsAVIb731Vo2P1+VvtikbPHiweOqpp5TfLRaLiI2NFfPnz3djqxpm3rx5om/fvk4fKywsFP7+/mL16tXKtoMHDwoAIisrq5FaeGMcPzusVquIjo4Wb7zxhrKtsLBQBAYGio8//lgIIcSBAwcEALFt2zZln6+//lpIkiTOnj3baG2vi5oCpPvuu6/G53jS+QkhRH5+vgAgvvvuOyFE3f4u6/KZ6irsYmtCioqK0LJlS+X3rKws/OpXv0JAQICyLSUlBYcPH8bly5eVfZKTkzWvk5KSgqysrMZptAs89dRTiIyMxODBg/Hhhx9CqKbm8vTzy8rKQu/evREVFaVsS0lJQXFxMfbv36/s40nn+Nprr6FVq1bo378/3njjDU1quy5/s01VRUUFduzYobkWJpMJycnJTfZaXM/Ro0cRGxuLTp06YezYscjJyQEA7NixA5WVlZpz7d69O9q3b++x53rixAnk5uZqzik8PByJiYnKOWVlZSEiIgKDBg1S9klOTobJZEJ2dnajt7khNm3ahDZt2qBbt274/e9/j0uXLimPedr5FRUVAYBy36vL32VdPlNdhYvVNhHHjh3D22+/jTfffFPZlpubi/j4eM1+8h9Fbm4uWrRogdzcXM0firxPbm6u8Y12gZdeegl33HEHQkJC8O233+LJJ5/ElStX8PTTTwNAjedXXFyMq1evIjg42B3NrrOa2i8/Vts+TfEcn376aQwYMAAtW7bEli1bMHPmTJw/fx6LFi0CULe/2abq4sWLsFgsTq/FoUOH3NSqhktMTMSKFSvQrVs3nD9/Hi+++CKGDRuGffv2ITc3FwEBAbq6OU/67HAkt7u2z8Pc3Fy0adNG87ifnx9atmzpEec9cuRI/OY3v0F8fDyOHz+OWbNm4e6770ZWVhbMZrNHnZ/VasXUqVNx6623olevXgBQp7/LunymugoDJBd77rnn8Prrr9e6z8GDBzXFkmfPnsXIkSMxZswYTJ482egm3pCGnF9t5syZo/zcv39/lJaW4o033lACJHdw9Tk2dfU53/T0dGVbnz59EBAQgMcffxzz58/32iUPPNXdd9+t/NynTx8kJiaiQ4cO+PTTT5tU0E1199BDDyk/9+7dG3369EHnzp2xadMmjBgxwo0tq7+nnnoK+/btw48//ujuptSIAZKLPfvss5gwYUKt+3Tq1En5+dy5c7j99tsxZMgQvP/++5r9oqOjddX78u/R0dG17iM/7mr1Pb/6SkxMxMsvv4zy8nIEBgbWeH5hYWGGfci78hyjo6N1o6Dqeg2NPEe1GznfxMREVFVV4eTJk+jWrVud/mabqsjISJjN5kb999SYIiIi0LVrVxw7dgx33nknKioqUFhYqPnfuiefq9zuvLw8xMTEKNvz8vLQr18/ZZ/8/HzN86qqqlBQUOCR592pUydERkbi2LFjGDFihMec35QpU7B27Vp8//33aNeunbI9Ojr6un+XdflMdRUGSC7WunVrtG7duk77nj17FrfffjsGDhyI5cuXw2TSloQlJSXh+eefR2VlJfz9/QEA69evR7du3ZSuiqSkJGRmZmLq1KnK89avX4+kpCTXnJCD+pxfQ+zevRstWrRQshFJSUm6aQuMPD/AteeYlJSEP//5z8jPz1dS3+vXr0dYWBhuuukmZZ/GPke1Gznf3bt3w2QyKedWl7/ZpiogIAADBw5EZmYmUlNTAdi6ATIzMzFlyhT3Ns4Frly5guPHj+PRRx/FwIED4e/vj8zMTDzwwAMAgMOHDyMnJ6fR/u5cLT4+HtHR0cjMzFQCouLiYmRnZ+P3v/89ANvfZ2FhIXbs2IGBAwcCADZs2ACr1YrExER3Nb3Bzpw5g0uXLikBYVM/PyEE/vCHP+Dzzz/Hpk2bdN3xdfm7rMtnqisbTG5w5swZ0aVLFzFixAhx5swZzbBNWWFhoYiKihKPPvqo2Ldvn1i1apUICQnRDfP38/MTb775pjh48KCYN29ekxnmf+rUKbFr1y7x4osvimbNmoldu3aJXbt2iZKSEiGEEF9++aX461//Kvbu3SuOHj0q/vd//1eEhISIuXPnKq8hD4GfPn26OHjwoHj33Xeb1BD4652jPCT1rrvuErt37xYZGRmidevWTof5N9VzlG3ZskW89dZbYvfu3eL48ePi//7v/0Tr1q3FuHHjlH3q8jfblK1atUoEBgaKFStWiAMHDojHHntMREREaEbMeIpnn31WbNq0SZw4cUJs3rxZJCcni8jISJGfny+EsA2nbt++vdiwYYPYvn27SEpKEklJSW5ude1KSkqUf2MAxKJFi8SuXbvEqVOnhBC2Yf4RERHiiy++EHv27BH33Xef02H+/fv3F9nZ2eLHH38UCQkJTWYYfG3nV1JSIqZNmyaysrLEiRMnxH//+18xYMAAkZCQIK5du6a8RlM+v9///vciPDxcbNq0SXPPKysrU/a53t9lXT5TXYUBkpssX75cAHD6pfbzzz+LoUOHisDAQNG2bVvx2muv6V7r008/FV27dhUBAQGiZ8+e4quvvmqs06jV+PHjnZ7fxo0bhRC24af9+vUTzZo1E6GhoaJv375i2bJlwmKxaF5n48aNol+/fiIgIEB06tRJLF++vPFPpgbXO0chhDh58qS4++67RXBwsIiMjBTPPvusqKys1LxOUz5H2Y4dO0RiYqIIDw8XQUFBokePHuLVV1/VfDgLUbe/2abs7bffFu3btxcBAQFi8ODB4qeffnJ3kxokLS1NxMTEiICAANG2bVuRlpYmjh07pjx+9epV8eSTT4oWLVqIkJAQcf/992v+g9YUbdy40em/t/HjxwshbEP958yZI6KiokRgYKAYMWKEOHz4sOY1Ll26JB5++GHRrFkzERYWJiZOnKj8h8bdaju/srIycdddd4nWrVsLf39/0aFDBzF58mRd8N6Uz6+me576864uf5d1+Ux1Bam60URERERUjfMgERERETlggERERETkgAESERERkQMGSEREREQOGCAREREROWCAREREROSAARIRERGRAwZIRERERA4YIBFRndx2222aNf/cacWKFZrFLI32q1/9CitXrmy047lLRUUFOnbsiO3bt7u7KURuxwCJiAxhZBCTlpaGI0eOGPLajr788kvk5eXhoYce0mzftWsX0tLSEBMTg8DAQHTo0AH33HMP/vOf/8BTFygICAjAtGnT8Kc//cndTSFyOwZIRORxgoODlZW8jbZ06VJMnDgRJpP94/KLL77ALbfcgitXruDvf/87Dh48iIyMDNx///2YPXs2ioqKGqVtRhg7dix+/PFH7N+/391NIXIrBkhEpFNaWopx48ahWbNmiImJwcKFC3X7lJeXY9q0aWjbti1CQ0ORmJiITZs2AQA2bdqEiRMnoqioCJIkQZIkvPDCCwAASZKwZs0azWtFRERgxYoVAICTJ09CkiR89tlnuP322xESEoK+ffsiKytL2d8xO/XCCy+gX79++Oijj9CxY0eEh4fjoYceQklJibJPSUkJxo4di9DQUMTExOCtt966brfhhQsXsGHDBvz617/WvDeTJk3C6NGj8dVXX+Guu+5Cp06d0KNHD0yaNAk///wzwsPDAQAWiwWTJk1CfHw8goOD0a1bNyxZskRzjAkTJiA1NRWvvvoqoqKiEBERgZdeeglVVVWYPn06WrZsiXbt2mH58uXKc+T36NNPP8WwYcMQHByMm2++GUeOHMG2bdswaNAgNGvWDHfffTcuXLigPG/btm248847ERkZifDwcAwfPhw7d+7UtKdFixa49dZbsWrVqhrfFyJfwACJiHSmT5+O7777Dl988QW+/fZbbNq0SXcjnTJlCrKysrBq1Srs2bMHY8aMwciRI3H06FEMGTIEixcvRlhYGM6fP4/z589j2rRp9WrD888/j2nTpmH37t3o2rUrHn74YVRVVdW4//Hjx7FmzRqsXbsWa9euxXfffYfXXntNeTw9PR2bN2/Gl19+ifXr1+OHH37QnZOjH3/8ESEhIejRo4ey7dtvv8WlS5cwY8aMGp8nSRIAwGq1ol27dli9ejUOHDiAuXPnYtasWfj00081+2/YsAHnzp3D999/j0WLFmHevHm455570KJFC2RnZ+OJJ57A448/jjNnzmieN2/ePMyePRs7d+6En58ffve732HGjBlYsmQJfvjhBxw7dgxz585V9i8pKcH48ePx448/4qeffkJCQgJGjRqlCSQBYPDgwfjhhx9qfW+IvJ4gIlIpKSkRAQEB4tNPP1W2Xbp0SQQHB4tnnnlGCCHEqVOnhNlsFmfPntU8d8SIEWLmzJlCCCGWL18uwsPDda8PQHz++eeabeHh4WL58uVCCCFOnDghAIi//e1vyuP79+8XAMTBgwedvva8efNESEiIKC4uVrZNnz5dJCYmCiGEKC4uFv7+/mL16tXK44WFhSIkJEQ5J2feeust0alTJ8221157TQAQBQUFyratW7eK0NBQ5es///lPja/51FNPiQceeED5ffz48aJDhw7CYrEo27p16yaGDRum/F5VVSVCQ0PFxx9/XON79PHHHwsAIjMzU9k2f/580a1btxrbYrFYRPPmzXXtXbJkiejYsWONzyPyBX7uDM6IqOk5fvw4KioqkJiYqGxr2bIlunXrpvy+d+9eWCwWdO3aVfPc8vJytGrVyiXt6NOnj/JzTEwMACA/Px/du3d3un/Hjh3RvHlzzXPy8/MBAL/88gsqKysxePBg5fHw8HDNOTlz9epVBAUF1amtu3fvBgAkJCRoMl3vvvsuPvzwQ+Tk5ODq1auoqKhAv379NM/v2bOnpsYpKioKvXr1Un43m81o1aqVcj7q46qfAwC9e/fWbFM/Jy8vD7Nnz8amTZuQn58Pi8WCsrIy5OTkaF43ODgYZWVl1z1vIm/GAImI6u3KlSswm83YsWMHzGaz5rFmzZrV+lxJknSjvCorK3X7+fv7a54D2LqsaqLeX35ObfvXRWRkJC5fvqzZlpCQAAA4fPgwbrnlFgBAYGAgunTponv+qlWrMG3aNCxcuBBJSUlo3rw53njjDWRnZ1+37XU5H2fvkeM29XPGjx+PS5cuYcmSJejQoQMCAwORlJSEiooKzesWFBSgdevWTt4RIt/BGiQi0ujcuTP8/f01N/HLly9rhtX3798fFosF+fn56NKli+YrOjoagG3IuMVi0b1+69atcf78eeX3o0ePGp6t6NSpE/z9/bFt2zZlW1FR0XWnCujfvz9yc3M1QdJdd92Fli1b4vXXX7/ucTdv3owhQ4bgySefRP/+/dGlSxccP3684SdygzZv3oynn34ao0aNQs+ePREYGIiLFy/q9tu3bx/69+/vhhYSNR0MkIhIo1mzZpg0aRKmT5+ODRs2YN++fZgwYYKmC6hr164YO3Ysxo0bh88++wwnTpzA1q1bMX/+fHz11VcAbF1eV65cQWZmJi5evKgEQXfccQfeeecd7Nq1C9u3b8cTTzyhy5a4WvPmzTF+/HhMnz4dGzduxP79+zFp0iSYTCYl8+JM//79ERkZic2bNyvbmjVrhr/97W/46quvMHr0aHzzzTf45ZdfsGfPHixYsAAAlKxaQkICtm/fjm+++QZHjhzBnDlzNEFaY0tISMBHH32EgwcPIjs7G2PHjkVwcLBuvx9++AF33XWXG1pI1HQwQCIinTfeeAPDhg3Dr3/9ayQnJ2Po0KEYOHCgZp/ly5dj3LhxePbZZ9GtWzekpqZi27ZtaN++PQBgyJAheOKJJ5CWlobWrVsrwcPChQsRFxeHYcOG4Xe/+x2mTZuGkJAQw89p0aJFSEpKwj333IPk5GTceuut6NGjR601RmazGRMnTsQ///lPzfb7778fW7ZsQUhICMaNG4du3brhjjvuwIYNG7Bq1Srcc889AIDHH38cv/nNb5CWlobExERcunQJTz75pKHnWZsPPvgAly9fxoABA/Doo4/i6aef1s0nlZWVhaKiIvz2t791UyuJmgZJOBYDEBH5gNLSUrRt2xYLFy7EpEmTatwvNzcXPXv2xM6dO9GhQ4dGbKF7pKWloW/fvpg1a5a7m0LkVswgEZFP2LVrFz7++GMcP34cO3fuxNixYwEA9913X63Pi46OxgcffKAb6eWNKioq0Lt3b/zxj390d1OI3I4ZJCLyCbt27cL/+3//D4cPH0ZAQAAGDhyIRYsWaYbFExHJGCAREREROWAXGxEREZEDBkhEREREDhggERERETlggERERETkgAESERERkQMGSEREREQOGCAREREROWCAREREROTg/wNTpRfNAJqZQAAAAABJRU5ErkJggg==", "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", " distribution_points = 1,\n", " fit_weight = 1,\n", " # resolution = 1, # only in version 1\n", " # kernel_type = \"Gauss\", # only in version 1\n", " instrument = nx.lib.instrument.Gaussian(1.0), # since version 2\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()" ] }, { "cell_type": "code", "execution_count": null, "id": "b6b59f0a-d539-4eda-be0e-f6934339686b", "metadata": {}, "outputs": [], "source": [] } ], "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.11.9" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }