{ "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": "iVBORw0KGgoAAAANSUhEUgAAAkoAAAGwCAYAAABWwkp7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAABabklEQVR4nO3deVxU5f4H8M+ZgZlhR0BWQdwVF3BBUnNL3CpbLOuWt1zKytAWslvebpndbvarm5k5N9vUdm21sjKVVDQNcCEXFEURUHaQbYAZmDm/PwZGCVGWmTnDzOf9es1LOHPmnO+MxzMfn+c5zxFEURRBRERERM3IpC6AiIiIyFYxKBERERG1gEGJiIiIqAUMSkREREQtYFAiIiIiagGDEhEREVELGJSIiIiIWuAkdQGdjcFgQG5uLjw8PCAIgtTlEBERUSuIoojKykoEBwdDJmt9OxGDUhvl5uYiNDRU6jKIiIioHXJyctCtW7dWr8+g1EYeHh4AjB+0p6enxNUQERFRa1RUVCA0NNT0Pd5aDEpt1Njd5unpyaBERETUybR12AwHcxMRERG1gEGJiIiIqAUMSkREREQtYFAiIiIiagGDEhEREVELGJSIiIiIWsCgRERERNQChwxKW7ZsQb9+/dCnTx988MEHUpdDRERENsrhJpysr69HfHw8du7cCS8vLwwfPhy33347fH19pS6NiIiIbIzDtSglJydj4MCBCAkJgbu7O6ZPn45t27ZJXRYRERHZoE4XlBITEzFjxgwEBwdDEARs3ry52TpqtRrh4eFQqVSIiYlBcnKy6bnc3FyEhISYfg8JCcGFCxesUToRERF1Mp0uKGk0GkRGRkKtVl/x+U2bNiE+Ph7Lli3DoUOHEBkZialTp6KwsNDKlRIREVFn1+mC0vTp0/Hyyy/j9ttvv+LzK1euxIIFCzBv3jxERERg7dq1cHV1xbp16wAAwcHBTVqQLly4gODg4Bb3p9VqUVFR0eRhSdp6Per1BgCAKIoW3RcRERFdnV0N5tbpdDh48CCWLl1qWiaTyRAbG4v9+/cDAEaOHIljx47hwoUL8PLywi+//ILnn3++xW2uWLECy5cvt3jtAJByrhTzN6TAVSGHh8oZZ4uqcF1PX7x06yD09ne3Sg1ERER0SadrUbqa4uJi6PV6BAQENFkeEBCA/Px8AICTkxPeeOMNTJw4EVFRUXjqqaeuesXb0qVLUV5ebnrk5ORYrP53dp1BZW09Ciq0yCisgkEE9p0pwY2r92DD75mmliYiIiKyDrtqUWqtW265Bbfcckur1lUqlVAqlVCr1VCr1dDr9RarK6OwCgAwsocPpg0MxMgePnj913TsPlWEF39Mw0f7szA7Jgwzh3WDj5vCYnUQERGRkV21KPn5+UEul6OgoKDJ8oKCAgQGBnZo23FxcUhLS0NKSkqHttOS2jo9ci5WAwDU9w7D/Ot7YFCIFzbMi8ZLtw6Et6szMos1ePmnE7julQTEb0o1BSsiIiKyDLsKSgqFAsOHD0dCQoJpmcFgQEJCAkaNGiVhZdeWXVoNUQQ8VE7wc7/UWiQIAu4fFY69z9yAl28bhEEhntDpDfj28AVMfnM3HvviMDIKKyWsnIiIyH51uq63qqoqZGRkmH7PzMxEamoqfHx8EBYWhvj4eMyZMwcjRozAyJEjsWrVKmg0GsybN69D+7V011txlRYAEOCpgiAIzZ53Vzrh79d1x9+v644/c8rw9m8Z2HGiAD/8mYsfj+TipsFBeHxSH/QJ8LBIfURERI5IEDvZNei7du3CxIkTmy2fM2cONmzYAABYs2YNXn/9deTn5yMqKgqrV69GTEyMWfZfUVEBLy8vlJeXw9PT0yzbBIBfjuZh4WeHMKJ7F3y9cHSrXnPsQjlWJ5zGtjRjV6MgADcNDsJjk/qgLwMTERGRSXu/vztdUJKapYLSF8nZWPrtUcQO8McHc6Lb9NrjucbA9OvxS4Fp+qBAxE3sjYHBXmarkYiIqLNq7/e3XY1RsiS1Wo2IiAhER7ctxLTWxWodAMDLpe1Xsw0M9sK7943Az4+NxbSBgRBF4Oej+bhp9V48sCEFh7MvmrtcIiIih8AWpTayVIvSip9P4N3Es3jw+h74180RHdpWen4l1DszsOVILgwNf7vj+3bFk5P7IirUu+PFEhERdTJsUerkyqrrAADers4d3la/QA+svmcodsSPx53Du0EuE7D7VBFuU/+OeeuTceR8WYf3QURE5AgYlGxEWU1D15ur+SaS7NnVHf+dFYmE+PG4Y1g3yARgZ3oRblnzOx7YkIKj58vNti8iIiJ7xKDUSpYeo2RqUXLpeIvSX4X7ueGNuyKR8NQEzBwaApkAJJwsxIw1e/HgRwdw7AIDExER0ZVwjFIbWWqMUk5pNYqqtOju4wpfd6XZtnslZ4qqsOa3DHyfesE0hmlKRACeiO2LiGDzvSciIiJbwekBrMRSQUkKGYVVePu30/jhz1w0HgXTBwXiycl9OQ8TERHZFQYlK7GnoNQoo7ASbyUYr5ITReM8TLdEBuOJ2L7o4ecmdXlEREQdxqBkJfYYlBqdKqjEm9tP4Zdj+QAAuUzAzKEheGxSH4T6uEpcHRERUfsxKFnY5fd6O3XqlF0GpUbHLpTjze2nkHCyEADgLBdw14hQLLqhN4K8XCSujoiIqO0YlKzEnluU/upQ9kW8uf0U9pwuBgAonGSYOzocC8f3Qhc3801jQEREZGkMSlbiSEGpUdLZEryx7RSSz5UCADyUTnh4fE/Mv74HXBVOEldHRER0bQxKVuKIQQkARFHErlNFeG1rOk7kVQAA/NyVeHxSb/xtZBic5ZySi4iIbBeDkpU4alBqZDCI+PFILt7YdgrZpdUAgO6+roif3BczhgRDJhMkrpCIiKg5BiULc6TB3K2hqzdgU0o23krIQHGVFgAQEeSJf0zrh/F9u0IQGJiIiMh2MChZiaO3KP2VRluP9b9n4t3dZ1GprQcAxPTwwTPT+2NYWBeJqyMiIjJiULISBqUru6jR4X+7MvDR/izo6g0AjLdFeXpqP/ThLN9ERCQxBiUrYVC6utyyGqzacQpfHzwPgwjIBODO4d3w1JR+CPBUSV0eERE5KAYlK2FQap2Mwkq8/ms6fj1eAABwcZZjwbieeHhcT7gpOaUAERFZF4OSlTAotc2h7Iv4z08ncDDrIgCgq4cST03ui1kjQiHnFXJERGQlDEpWwqDUdqIoYuuxfLy69SSySoxTCvQL8MA/bxqA8X27SlwdERE5AgYlC+P0AB2nqzfgkz+ysDrhNMpr6gAAY/v44bmbBqB/ID9LIiKyHAYlK2GLUseVVeuw5rcMfLT/HOr0ImQCMGt4KOKn9OWAbyIisggGJSthUDKfrBINXtuajp+O5gEwDvh+eHxPPDSuJ+8hR0REZsWgZCUMSuZ3MKsUL/90AoezywAA/h5KLJnSD3cM78YB30REZBYMSlbCoGQZoiji56P5eHXrCeSU1gAABoV4YtmMgYgO95G4OiIi6uwYlKyEQcmytPV6fLzPOOC78ZYot0QG49np/RHs7SJxdURE1FkxKFkJg5J1FFdp8ca2dGxMyYEoGscvLZzQCw+N6wmVs1zq8oiIqJNhULISBiXrOnahHMt/PI6Uc8YJK0O8XfDcTQMwfVAgBIHjl4iIqHUYlKyEQcn6RFHEj0fysOLnE8grrwUAXNfTB8tmDMSAIP4dEBHRtTEoWQmDknSqdfVYu/ss3t19Btp6A2QCMDumO5ZM6QcvV2epyyMiIhvW3u9vmQVrsitqtRoRERGIjo6WuhSH5apwQvzkvkh4ajxuGhwEgwh88kcWbnhjF74+eB7M/EREZG5sUWojtijZjn1nivHC98eRUVgFAIgO74J/3zaIt0MhIqJm2KJEDmd0Lz/8/NhYPDOtP1yc5Ug5dxE3rd6Lf29JQ1XD1AJEREQdwaBEnZrCSYaFE3phx1PjMW1gIPQGER/uzcSkN3bhhz9z2R1HREQdwqBEdiHE2wVr7xuODfOi0d3XFQUVWjz2xWH8/cMkU9ccERFRWzEokV2Z0M8fvz4xDk/G9oXCSYbfM0ow/a1EvLb1JKp17I4jIqK2YVAiu6NyluPx2D7Y8eR43NDfH3V6Ef/bdQaTVyYi4USB1OUREVEnwqBEdivM1xUfzhmB9+4bjhBvF1woq8EDHx3Aos8PoahSK3V5RETUCTAokV0TBAFTBgZie/w4PDSuJ2QCsOVIHmJX7saXKTkc7E1ERFfFoEQOwVXhhH/eOAA/LLoeA4M9UV5Th398cwT3vp+EzGKN1OUREZGNYlAihzIoxAvfx43BP2/sD5WzDPvPlmDqqkSod2agTm+QujwiIrIxDErkcJzkMjw0rhe2PTEeY/v4QVdvwOu/pmPG23uRmlMmdXlERGRDHDIo3X777ejSpQvuvPNOqUshCYX5uuLj+SOx8q5IdHF1xsn8Stz+v9+x/Mfj0HBmbyIigoMGpccffxwff/yx1GWQDRAEATOHdcOO+PGYOTQEogis//0cpryZiL2ni6Uuj4iIJOaQQWnChAnw8PCQugyyIb7uSqy8Owofzx+Jbl2MUwn8/cMkPPfdUd43jojIgdlcUEpMTMSMGTMQHBwMQRCwefPmZuuo1WqEh4dDpVIhJiYGycnJ1i+U7NK4vl3x6xPjcP+o7gCAz5KyMW1VIvZlsHWJiMgR2VxQ0mg0iIyMhFqtvuLzmzZtQnx8PJYtW4ZDhw4hMjISU6dORWFhoWmdqKgoDBo0qNkjNze3zfVotVpUVFQ0eZB9c1M64aVbB+HzB2PQrYsLzl+swb0fJOH5zcc4domIyMEIog3PuCcIAr777jvcdtttpmUxMTGIjo7GmjVrAAAGgwGhoaFYvHgxnn322VZve9euXVizZg2+/vrrq6734osvYvny5c2Wl5eXw9PTs9X7o86pSluPFT+fwGdJ2QCAUB8XvHZHJEb18pW4MiIiaouKigp4eXm1+fvb5lqUrkan0+HgwYOIjY01LZPJZIiNjcX+/fstss+lS5eivLzc9MjJybHIfsg2uSud8J/bB+OzB2MQ4u2CnNIa3PP+H1j2/THeZJeIyAF0qqBUXFwMvV6PgICAJssDAgKQn5/f6u3ExsZi1qxZ+Pnnn9GtW7erhiylUglPT88mD3I8Y3r7YesTY3HPyDAAwEf7szD9rT04mHVR4sqIiMiSnKQuQAo7duxo82vUajXUajX0er0FKqLOwEPljBUzB+PGwYF45usjyCqpxqy1+7BoYm8sntQHzvJO9f8OIiJqhU51Zvfz84NcLkdBQUGT5QUFBQgMDLTovuPi4pCWloaUlBSL7ods39g+XbH1yXG4fWgIDCKw+rcM3PnOPpwpqpK6NCIiMrNOFZQUCgWGDx+OhIQE0zKDwYCEhASMGjVKwsrI0XiqnPHm3VF4+56h8FQ54c/z5bhp9R588kcWbPj6CCIiaiObC0pVVVVITU1FamoqACAzMxOpqanIzjZedRQfH4/3338fH330EU6cOIGFCxdCo9Fg3rx5Fq1LrVYjIiIC0dHRFt0PdS4zIoPx65PjMKa3L2rrDHh+8zHM35CCwspaqUsjIiIzsLnpAXbt2oWJEyc2Wz5nzhxs2LABALBmzRq8/vrryM/PR1RUFFavXo2YmBir1NfeywvJvhkMItbvO4f/23oSunoDfNwUWDFzMKYOtGyXMBERtU57v79tLijZOgYlupr0/Eo8sSkVJ/KME5Ped113PHfTAKic5RJXRkTk2BxiHiUpseuNWqNfoAc2x43GQ+N6AgA++SMLt6l/R0ZhpcSVERFRe7BFqY3YokSttSu9EE99+SdKNDq4OMux/NaBmDW8GwRBkLo0IiKHwxYlIhszoZ8/fnl8LMb09kVNnR7/+PoIHt+YisraOqlLIyKiVmJQIrIgf08VPpkfg6en9oNcJuCHP3Nx89t7ceR8mdSlERFRKzAotRLHKFF7yWQC4ib2xpcPj0KItwuySqpxxzv7sP73TM65RERk4zhGqY04Rok6orymDs9+cwS/HDPem/DmIUH4vzuGwE3pkHcTIiKyGo5RIuoEvFyc8b/Zw7BsRgScZAK2HMnDrbwqjojIZjEoEVmZIAiYN6YHNj18HQI8lcgorMIta37HD3/mSl0aERH9BYNSK3GMEpnb8O4++OmxsRjdyxfVOj0e++IwXvzhOHT1BqlLIyKiBhyj1EYco0TmpjeIWLk9HeqdZwAAQ8O88b/ZwxDk5SJxZURE9oNjlIg6KblMwNNT++OD+0fAQ+WEw9lluGn1Xuw7Uyx1aUREDo9BichGxEYE4KfFYzEw2BOlGh3u+zAZH+8/xykEiIgkxKBEZEPCfF3xzcLRuDUqGHqDiBe+P45/fneU45aIiCTCoERkY1TOcqy6OwpLp/eHIABfJOfg3vf/QFGlVurSiIgcDoNSK/GqN7ImQRDw8PheWD83Gh4qJxzIuohb1uzF0fPlUpdGRORQeNVbG/GqN7K2s0VVePDjAzhbpIHSSYbX7hyCW6NCpC6LiKhT4VVvRHaqZ1d3bI4bgxv6+0Nbb8DjG1Oxcls6B3kTEVkBgxJRJ+Cpcsb794/Awgm9AACrf8vA4xtTUVunl7gyIiL7xqBE1EnIZQKemdYfr905BE4yAT/8mYu/f5CEkioO8iYishQGJaJO5q4Rofh4/kh4NgzynvnOPpwpqpK6LCIiu8SgRNQJje7th28fHYNQHxdklVRj5v/2Yf+ZEqnLIiKyOwxKrcTpAcjW9PZ3x+ZHx2BYmDfKa+pw/7ok/PBnrtRlERHZFU4P0EacHoBsTW2dHk99+Sd+OpoHAHjh5gjMv76HxFUREdkWTg9A5KBUznK8fc9QzB0dDgB4aUsaXv3lJKcPICIyAwYlIjsgkwlYNiMC/5jWDwCwdvcZPPXVn6jT8x5xREQdwaBEZCcEQcCjE3rj9TuHQC4T8O2hC3jwowPQaOulLo2IqNNiUCKyM7NGhOKD+0fAxVmO3aeKcO8HSSir1kldFhFRp8SgRGSHJvb3x+cLYuDt6ow/c8rwt/f+QFElJ6YkImorBiUiOzU0rAu+fHgUunoocTK/Ene/ux+5ZTVSl0VE1KkwKBHZsb4BHvjq4VEI8XbB2WINZq3dj6wSjdRlERF1GgxKRHYu3M8NXz4yCj383HChrAaz1u7H6YJKqcsiIuoUGJRaiTNzU2cW4u2CTQ9fh34BHiis1OKud/cjLbdC6rKIiGweZ+ZuI87MTZ3ZRY0Oc9Yn48j5cnRxdcbnC67DgCAex0Rk/zgzNxFdUxc3BT55IAaR3bxwsboOsz9Iwsl8tiwREbWEQYnIwXi5OOPjB2IwpJsXSjU63Pt+EtLzOWaJiOhKGJSIHJCXizM+mR+DwSGNYekPhiUioitgUCJyUF6uzvj0gRgMCvFESUNY4tVwRERNMSgRObDGsDQw2BiWZn+QhOySaqnLIiKyGQxKRA7O21WBTx+IQd8AdxRWavH3D5NQUFErdVlERDaBQYmI0MXNGJbCfFyRXVqN+z5MwkUNb6RLRMSgREQAAH9PFT57MAYBnkqcKqjCnPXJqKytk7osIiJJMSgRkUmojys+fSAGXVydceR8OR786ABq6/RSl0VEJBkGJSJqok+ABz6eHwN3pROSMkvxxMZU6A2cwJ+IHJPDBaWcnBxMmDABERERGDJkCL766iupSyKyOYO7eeH9+0dAIZdh6/F8/HtLGni3IyJyRA4XlJycnLBq1SqkpaVh27ZteOKJJ6DRaKQui8jmjOrli//eFQkA2LDvHD7cmylxRURE1udwQSkoKAhRUVEAgMDAQPj5+aG0tFTaoohs1C2Rwfjnjf0BAC//dAI//pkrcUVERNZlc0EpMTERM2bMQHBwMARBwObNm5uto1arER4eDpVKhZiYGCQnJ7drXwcPHoRer0doaGgHqyayXwvG9sTc0eEAgKe+/BN/nC2RtiAiIiuyuaCk0WgQGRkJtVp9xec3bdqE+Ph4LFu2DIcOHUJkZCSmTp2KwsJC0zpRUVEYNGhQs0du7qX/DZeWluL+++/He++9d9V6tFotKioqmjyIHIkgCHj+5ghMHRgAnd6Ahz4+gLNFVVKXRURkFYJowyM0BUHAd999h9tuu820LCYmBtHR0VizZg0AwGAwIDQ0FIsXL8azzz7bqu1qtVpMnjwZCxYswH333XfVdV988UUsX7682fLy8nJ4enq2/s0QdXK1dXrc+/4fOJRdhp5+bvju0THwcnWWuiwiolapqKiAl5dXm7+/ba5F6Wp0Oh0OHjyI2NhY0zKZTIbY2Fjs37+/VdsQRRFz587FDTfccM2QBABLly5FeXm56ZGTk9Pu+ok6M5WzHGvvG45gLxXOFmuw6ItDqNcbpC6LiMiiOlVQKi4uhl6vR0BAQJPlAQEByM/Pb9U2fv/9d2zatAmbN29GVFQUoqKicPTo0RbXVyqV8PT0bPIgclT+Hiq8d/8IuDjLsed0MV75+aTUJRERWZST1AVY2/XXXw+Doe3/C1ar1VCr1dDrOUsxObZBIV54465IPPrZIaz7PRP9At1xd3SY1GUREVlEp2pR8vPzg1wuR0FBQZPlBQUFCAwMtOi+4+LikJaWhpSUFIvuh6gzuHFwEJ6I7QMA+NfmY0g5xyk2iMg+daqgpFAoMHz4cCQkJJiWGQwGJCQkYNSoURJWRuR4HruhD24cHIg6vYiFnx5EQUWt1CUREZldh4JSXV0dcnJykJ6ebrZJG6uqqpCamorU1FQAQGZmJlJTU5GdnQ0AiI+Px/vvv4+PPvoIJ06cwMKFC6HRaDBv3jyz7L8larUaERERiI6Otuh+iDoLmUzAf2dFon+gB4qrdFj0+SHUcXA3EdmZNk8PUFlZiU8//RQbN25EcnIydDodRFGEIAjo1q0bpkyZgoceeqjdgWLXrl2YOHFis+Vz5szBhg0bAABr1qzB66+/jvz8fERFRWH16tWIiYlp1/7aqr2XFxLZq8xiDWa8vRdV2nosGNsDz90UIXVJRETNtPf7u01BaeXKlfjPf/6DXr16YcaMGRg5ciSCg4Ph4uKC0tJSHDt2DHv27MHmzZsRExODt99+G3369GnXG7JVDEpEzW09lodHPj0EAFj792GYNihI4oqIiJqySlC655578K9//QsDBw686nparRbr16+HQqHA/PnzW12MLbv8qrdTp04xKBH9xctb0vDB3kx4KJ3ww+Lr0cPPTeqSiIhMrBKUiC1KRC2p0xtw7/t/IOXcRfQP9MDmuDFQOculLouo3QwGERpdPTRaPaq09ait00OnN0BXb3zUNf6sb/xdhL5h+pnGL9bGb9jLv2plMgFOMhmc5QIUTjI4y2VwkglwdpJBIW/4XS5A6SSDq8IJbgo5XJVOcHGWQy4TrPwp2I/2fn+3ex6lmpoaiKIIV1dXAEBWVha+++47DBgwAFOnTm3vZomok3KWy7Dm3mG4afUenMyvxKu/nMSLt1y99ZnIGgwGEaXVOhRVanFRo0Nptc74p6YOF6t1KNXocLFah/KaOlRp61FVWw+Nth4ane3Nm6dyNoYnV4W84eEEd6UTPF2c4OXiDE8XZ3i18PB2UcDTxQmCwLDVFu1uUZoyZQpmzpyJRx55BGVlZejfvz+cnZ1RXFyMlStXYuHCheau1SawRYno6nalF2LueuN8Y+vmjsAN/QOu8QqijqmorUN2STWyS6uRW1aD/PJa5FfUmv4sqKhFnb79nSdymQA3hRwuCjkUDa0+Cic5FJe1CF3eMnR5DhHQ5BcAxuBWpxdRbzC2StXVi6i7ws+1dQbU6PSo1tXDYKa+HyeZAF93BXzdlPB1V8DPXQk/dwV83ZXwdTP+3ri8q4cSzvJONYvQVVm9RenQoUN48803AQBff/01AgICcPjwYXzzzTd44YUX7C4ocWZuotaZ0M8f88f0wLrfM7HkqyPY+vhY+HuqpC6LOrkanR5niqpwqqAS54o1yCqtRlZJNbJKNLhYXXfN1wsC4OOqgI+bAl3cFPBxbfjTzRldGpZ7uTjDXekEN6WxlcZdZfxT6SSTtBVGFEVo6w3QaOtRrdM3PIw/a7T1qNLWo6KmDuU19SivqTM9Ki7/ubYO1To96g0iCiq0KKjQXnO/ggD4uikR6KVEoKdLw58qBHq5NPypRICnCh4q+745drtblFxdXXHy5EmEhYXhrrvuwsCBA7Fs2TLk5OSgX79+qK6uNnetNoEtSkTXpq3X4zb1PpzIq8DYPn74aN5IyDi2glrBYBBxtliD47nlSM+vxKmCKpwurER2aTWu9m3l565AmI8rQrq4ItDT+AUe5NXw5e7lAn87ax1pD229HqUaHYordSjWaFFSpUNJlRYlGh2Kqxp+12hRXGn8s7WtcO5KJwR5qdCtiwu6dXFFqE/Dn11c0a2LC7xdnW2iu8/qLUq9e/fG5s2bcfvtt+PXX3/Fk08+CQAoLCxkgCBycEonOVb/LQoz1uzFntPFWPd7Jh4c21PqssjGiKKIrJJqHLlQjqPny3DkfDmOXShvcWxQF1dn9A3wQM+u7uju64pwX1eE+bghzNcV7kqHu3Vpmymd5AjyckGQl8s1120c15Vf3rQLM6/c+GfjsspaY4vW6cIqnC6suuK23JVOphDVrYsLQn1cEdrFBWG+rgjzcYWrwrb/7trdovT111/j3nvvhV6vx6RJk7Bt2zYAwIoVK5CYmIhffvnFrIXaCrYoEbXeZ0lZeO67Y1A4yfDT4uvRJ8BD6pJIQnV6A47nViAlsxTJ50px4FzpFbvNVM4yRAR5on+QJ/r6u6NvgAf6BHjAz11hEy0TdIlGW4/8ilrkltXg/MUa5JRWG/+8aPyzqPLaXXxdPZQI83FFdx9XhPm6oruvK0b38kOAmbvsJZkeID8/H3l5eYiMjIRMZmzSTE5OhqenJ/r379/ezdokzqNE1HaiKGL+hhTsTC9CZDcvfLNwNJwcvPvDkRgMItLyKrD7VBH2nynBoeyLqP5La5HCyRiKhnTzwuAQLwzp5o1eXd14nNiJ2jp9k+B0viFIZZcaB9+X11x5fNkH949AbIR5LwThPEpWwhYlorbJL6/F5Dd3o7K2Hv+Y1g+PTugtdUlkQUWVWiSeKkLi6SLsPV2MEo2uyfOeKieM7OGD6HAfRPfwwaBgLyicGIocVXl1HbJKNchuGJzfePXiKzMHm33SWkmCUkJCAhISElBYWAiDoenNMNetW9fezdo0BiWitvvqQA6e/voIFHIZtjx2PfqyC86uZJdU49fj+fj1eD4OZl9sMujaTSHHqF5+GNvHDzE9fdDX34MD+0kSVh/MvXz5crz00ksYMWIEgoKC2G9MRC26c3g3/HIsH7+dLMSSr/7Et+yC6/TOFWvwfWouth7Px4m8iibPDQz2xPi+XTGub1cMC+vCFiPq1NrdohQUFITXXnsN9913n7lrsmlsUSJqn4KKWkxeuRsV7ILrtC5qdNhyJBffHr6Aw9llpuVymYCR4T6YOjAAUwYGItj72ldVEVmb1VuUdDodRo8e3d6XE5GDCfBU4YUZA7Hkqz+xOuE0ZgwJRqiPq9Rl0TXoDSJ2pRdiU0oOdqYXmubWkQnA9X264uYhQYgdEAAfN4XElRJZRruD0oMPPojPP/8czz//vDnrsVmcmZuo4+4YFoJvDp7H/rMleP77Y1g/N5rd9jaqpEqLTQdy8Nkf2bhQVmNaHhHkiZnDQnBLZDBnXCeH0O6ut8cffxwff/wxhgwZgiFDhsDZuekU5itXrjRLgbaGXW9EHXOmqArTV+2BTm/A/2YPw42Dg6QuiS5z7EI5PthzFj8fzYdOb7xIx8vFGbOGd8OsEaHoF8iB+NQ5Wb3r7ciRI4iKigIAHDt2rMlz/B8iEbWkV1d3PDK+J1b/loHlPx7H2D5+dn+vKFsniiL2nSnB2t1nsOd0sWl5ZDcv/P267pgRGQyVs1zCComk0+6gtHPnTnPWQUQO5NGJvfHDn7k4V1KNldtPYdmMgVKX5JBEUcSOE4V4+7fTOHK+HIBxYPbNQ4Iwf0wPRIZ6S1sgkQ1o8zWbL7zwAg4ePGiJWojIQaic5fj3bYMAAB/vz8LpgkqJK3Isoihiz+ki3Pa/fVjw8QEcOV8OlbMMc0Z1x64lE/DW34YyJBE1aHOL0vnz5zF9+nQoFArMmDEDt9xyCyZNmgSFglc8EFHrje3TFVMiArAtrQAvbUnDx/NHstveCg5lX8RrW0/ij7OlAAAXZznmjgnHg9f3gK+7UuLqiGxPuwZzGwwG/P777/jxxx/x/fffIy8vD5MnT8att96Km2++GT4+Ppao1SZwMDeR+WSVaDB5ZSJ0eoNF7u1ElxRU1OLVX07iu8MXAAAKuQyzrwvDoxN6o6sHAxLZP0nv9XbixAlTaDp48CBGjhyJW265Bffccw9CQkI6unmbwJviElnG/209iXd2nUG4ryt+fXIclE4cNGxO2no9PtybiTW/ZaBap4cgAHcO64YnJvdFCCeGJAditaD00EMPYeXKlXB3d7/i80VFRfjhhx/www8/YOzYsViyZElbNm/z2KJEZF5V2npM/O8uFFVq8c8b++Ohcb2kLslu7DtTjH9+exTnSqoBAMPCvPHiLQMxpJu3tIURScBqQUkulyMvLw/+/v5tLtIeMCgRmV/jTXM9lE5I/MdEdOEszx1SWVuHFb+cxOdJ2QAAfw8lnp3eH7dFhfCGtOSw2vv93ear3szQU0dE1MQdw7phQJAnKrX1eGf3GanL6dR2nyrClDcTTSHp3pgwJDw1HjOHdWNIImoH3tKZiCQnkwn4x9R+AICP9p1DfnmtxBV1Ptp6Pf69JQ1z1iUjr7wWYT6u+GLBdXjl9sGc0JOoA9oclARB4CW8RGR2E/p1RXR4F2jrDXgr4bTU5XQqZ4qqMPN/+/Dh3kwAwJxR3bH1ibEY1ctX4sqIOr82z6MkiiLmzp0LpfLql5N+++237S6KiByPIAj4x7T+mLV2P748kIOHxvVEDz83qcuyeZsPX8DSb4+ipk6PLq7OeP3OSE6zQGRGbQ5Kc+bMsUQdRESIDvfBxH5dsTO9CG9sS8eae4dJXZLN0htE/N/Wk3gv8SwAYFRPX7x5dxQCvVQSV0ZkX8wyj5Ij4VVvRJZ1PLccN63eCwDY+sRY9A/kv7O/Kq+pw2NfHMbuU0UAgEcn9MJTU/pBzsHaRC2y2lVvRESWNDDYCzcODgQArPktQ+JqbE9WiQa3q3/H7lNFUDnL8PY9Q/GPaf0ZkogspE1BKTs7u00bv3DhQpvWt2VqtRoRERGIjo6WuhQiu7doYh8AwE9H83CmqEriamzHsQvluOOdfThbrEGwlwpfPzIaMyKDpS6LyK61KShFR0fj4YcfRkpKSovrlJeX4/3338egQYPwzTffdLhAWxEXF4e0tLSrvnciMo+IYE/EDgiAKALqnWxVAoB9GcX423t/oLhKh4ggT2xeNAaDQrykLovI7rVpMHdaWhr+85//YPLkyVCpVBg+fDiCg4OhUqlw8eJFpKWl4fjx4xg2bBhee+013HjjjZaqm4js3OIbemPHiQJ8n5qLJyb1RZivq9QlSebX4/lY/Plh6PQGjOrpi/fuH865kYispF2DuWtqavDTTz9h7969yMrKQk1NDfz8/DB06FBMnToVgwYNskStNoGDuYms5/51yUg8VYR7RoZixcwhUpcjiW3H8/HoZ4dQbxAxfVAg3rw7Cipn3jiYqK2sdq83R8egRGQ9KedKMWvtfiicZNj37A3wc7/6/G32ZkdaARZ+dhB1ehG3RgVj5V1RHLRN1E686o2I7M6I7l0QGeoNXb0Bn/6RJXU5VrXzZKEpJM2IDMYbsyIZkogkwKBERDZLEAQ8eH0PAMAn+7NQW6eXuCLrOJR90RSSbhochDfvioSTnKdrIinwXx4R2bTpgwIR7KVCiUaHH1JzpS7H4s4WVeGBDSmorTNgQr+uWPW3KIYkIgnxXx8R2TQnuQxzx4QDAD7Yexb2PKyyqFKLOeuTcbG6DpHdvPC/2cPgzJBEJKl2/wucM2cOEhMTzVkLEdEV3R0dBjeFHKcKqrDndLHU5ViErt6AhZ8eRE5pDcJ9XfHh3Gi4Ktp8O04iMrN2B6Xy8nLExsaiT58+eOWVV+xqFm4isi1eLs6YNSIUAOx2UPeLPx7HgayL8FA54cO50Q53hR+RrWp3UNq8eTMuXLiAhQsXYtOmTQgPD8f06dPx9ddfo66uzpw1mlVZWRlGjBiBqKgoDBo0CO+//77UJRFRK8yOCQMAJJwsREFFrcTVmNenf2Th86RsCAKw+m9D0auru9QlEVGDDnV+d+3aFfHx8fjzzz+RlJSE3r1747777kNwcDCefPJJnD592lx1mo2HhwcSExORmpqKpKQkvPLKKygpKZG6LCK6hj4BHogO7wK9QcSmlBypyzGbo+fLsfzH4wCAp6f2w8T+/hJXRESXM8sowby8PGzfvh3bt2+HXC7HjTfeiKNHjyIiIgJvvvmmOXZhNnK5HK6uxlshaLVaiKJo14NDiezJ7JjuAICNydnQGzr/v9sqbT0Wf3EIdXoRUwcGYOH4XlKXRER/0e6gVFdXh2+++QY333wzunfvjq+++gpPPPEEcnNz8dFHH2HHjh348ssv8dJLL7Vpu4mJiZgxYwaCg4MhCAI2b97cbB21Wo3w8HCoVCrExMQgOTm5TfsoKytDZGQkunXrhqeffhp+fn5tej0RSWPaoEB0cXVGbnktdqUXSl1Oh73w/TGcK6lGsJcK/3fHEAgCJ5QksjXtDkpBQUFYsGABunfvjuTkZBw4cACPPPJIk2nBJ06cCG9v7zZtV6PRIDIyEmq1+orPb9q0CfHx8Vi2bBkOHTqEyMhITJ06FYWFl06ajeOP/vrIzTXOweLt7Y0///wTmZmZ+Pzzz1FQUND2D4CIrE7lLMcdw7oBAD5Pypa4mo75PvUCvj10ATIBWPW3ofB2VUhdEhFdQbvv9fbJJ59g1qxZUKlU5q7JRBAEfPfdd7jttttMy2JiYhAdHY01a9YAAAwGA0JDQ7F48WI8++yzbd7Ho48+ihtuuAF33nnnFZ/XarXQarWm3ysqKhAaGsp7vRFJ5ExRFSa9sRsyAdj37CQEelnuHGQphZW1mLwyEeU1dXh8Uh88Obmv1CUR2T2r3+tt/PjxUCqbX74qiiKysy3zPz2dToeDBw8iNjbWtEwmkyE2Nhb79+9v1TYKCgpQWVkJwDjFQWJiIvr169fi+itWrICXl5fpERoa2rE3QUQd0qurO6LDu8AgGltlOqMXfziO8po6DArxxOIbektdDhFdRbuDUo8ePVBUVNRseWlpKXr06NGholpSXFwMvV6PgICAJssDAgKQn5/fqm1kZWVh7NixiIyMxNixY7F48WIMHjy4xfWXLl2K8vJy0yMnx36utiHqrG4faux++/bQhU53McbWY3n4+Wg+5DIB/3fHEN6ehMjGtXvaV1EUrzjwsKqqyqLdcR01cuRIpKamtnp9pVIJpVIJtVoNtVoNvd4xbspJZMtuGhyEF384jvSCSqTlVWBgsJfUJbVKeXUdnv/eOBXAI+N7dpq6iRxZm4NSfHw8AOP4oeeff950qT0A6PV6JCUlISoqymwFXs7Pzw9yubzZ4OuCggIEBgZaZJ+N4uLiEBcXZ+rjJCLpeLk6Y9IAf/xyLB/fHbrQaQLHyu3pKKrUoldXNyy+oY/U5RBRK7S5zffw4cM4fPgwRFHE0aNHTb8fPnwYJ0+eRGRkJDZs2GCBUgGFQoHhw4cjISHBtMxgMCAhIQGjRo2yyD6JyDbNbLj67fs/c1GvN0hczbWdzK/Apw1X6v371kFQOcslroiIWqPNLUo7d+4EAMybNw9vvfWW2a/8qqqqQkZGhun3zMxMpKamwsfHB2FhYYiPj8ecOXMwYsQIjBw5EqtWrYJGo8G8efPMWsdfseuNyLaM79sVXVydUVSpxR9nS3F9H9udD00URSz/IQ16g4hpAwMxurft1kpETbV7egBL2bVrFyZOnNhs+Zw5c0wtVWvWrMHrr7+O/Px8REVFYfXq1YiJibFKfe29vJCIzG/pt0fwRXIO7o0Jwyu3t3xRhtS2HsvHI58ehNJJhh3x4xHq43rtFxGRWbX3+7tNQSk+Ph7//ve/4ebmZhqr1JKVK1e2uojOhEGJyHYknirC/euS4eumQPJzsZDLbG9ma71BxJQ3d+NMkQaLJvbGkqktT0dCRJbT3u/vNnW9HT58GHV1daafW8Jp+InIGkb18oW3qzNKNDokZ5ZiVC9fqUtq5rvDF3CmSANvV2c8NL6n1OUQURu1KSg1jk/668+OgGOUiGyPs1yGKREB+PLAefx8NM/mgpK2Xo83t58CACwc3wueKmeJKyKitmr3TGc1NTWorq42/Z6VlYVVq1Zh27ZtZinM1sTFxSEtLQ0pKSlSl0JEl5k+OAgAsPV4PvQGmxpyiS9TcnChrAb+HkrcPypc6nKIqB3aHZRuvfVWfPzxxwCAsrIyjBw5Em+88QZuvfVWvPPOO2YrkIjoasb08oOnyglFlVoczLoodTkmdXoD1u4+CwBYdENvuCg4HQBRZ9TuoHTo0CGMHTsWAPD1118jMDAQWVlZ+Pjjj7F69WqzFWgr1Go1IiIiEB0dLXUpRHQZhZMMkwYYb2uUcKLgGmtbz5YjubhQVgM/dwXuGsF7RBJ1Vu0OStXV1fDw8AAAbNu2DTNnzoRMJsN1112HrKwssxVoK9j1RmS7Jvb3BwD8drJQ4kqMDAYR7+w6AwCYN6YHJ5ck6sTaHZR69+6NzZs3IycnB7/++iumTJkCACgsLORl80RkVeP7dIVcJuB0YRVySquv/QIL25leiFMFVXBXOuHv13WXuhwi6oB2B6UXXngBS5YsQXh4OGJiYky3ENm2bRuGDh1qtgKJiK7Fy9UZw7t3AWAbrUrv7zGOTZp9XRi8XHilG1Fn1u6gdOeddyI7OxsHDhzA1q1bTcsnTZqEN9980yzFERG11iQb6X47XVCJP86WQiYAc3ilG1Gn1+6gBACBgYEYOnQoZLJLmxk5ciT69+/f4cJsDQdzE9m2GxqC0v6zJajW1UtWx6d/GMdoxg4IQLC3i2R1EJF5tPmmuJdLSEhAQkICCgsLYTA0vXv3unXrOlSYrYmLi0NcXJxpCnQisi29/d3RrYsLzl+swe8ZJZgcEWD1GjTaenxz6AIA4L5RHJtEZA/a3aK0fPlyTJkyBQkJCSguLsbFixebPIiIrEkQBEzsZ2xV2nO6SJIaNqdeQJW2Hj383DCml58kNRCRebW7RWnt2rXYsGED7rvvPnPWQ0TUbmN6++KTP7Lwe0axJPv/PCkbADA7JgwyG7xBLxG1XbtblHQ6HUaPHm3OWoiIOuS6nr4QBOBMkQb55bVW3fepgkocz62Ak0zAzGHdrLpvIrKcdgelBx98EJ9//rk5a7FpHMxNZPu8XRUYFGwcQ7jvjHVblTYfNo5NmtDPHz5uCqvum4gsp91db7W1tXjvvfewY8cODBkyBM7OTecKWblyZYeLsyUczE3UOYzu7YujF8rxe0aJ1Vp2DAYR36fmAgBuHxpilX0SkXW0OygdOXIEUVFRAIBjx441eU4Q2DdPRNIY08sP7+4+i31niiGKolXORynnSnGhrAYeSidMGuBv8f0RkfW0Oyjt3LnTnHUQEZlFdLgPnOUC8sprcf5iDUJ9XC2+z82pxm63aYMCeV83IjvToQkn9+zZg7///e8YPXo0Llwwnig++eQT7N271yzFERG1lYtCjoEN45RSzpVafH91egN+PpoPALiN3W5EdqfdQembb77B1KlT4eLigkOHDkGr1QIAysvL8corr5itQCKitooON9737UCW5ed0SzpbivKaOvi6KXBdT1+L74+IrKvdQenll1/G2rVr8f777zcZyD1mzBgcOnTILMUREbXHiHAfAMABK7Qo/Xrc2Jo0OSIAcs6dRGR32h2U0tPTMW7cuGbLvby8UFZW1pGabBKnByDqPEZ0N7YonSqoQlm1zmL7MRhEbE8rAABMHRhosf0QkXTaHZQCAwORkZHRbPnevXvRs2fPDhVli+Li4pCWloaUlBSpSyGia/B1V6JnVzcAwEELdr8duVCO/IpauCnkGN2b3W5E9qjdQWnBggV4/PHHkZSUBEEQkJubi88++wxLlizBwoULzVkjEVGbRXc3dr+lnLNcUNqeZux2m9DfH0onXu1GZI/aPT3As88+C4PBgEmTJqG6uhrjxo2DUqnEkiVLsHjxYnPWSETUZsO6e2PTgRyk5lguKO0+Zbz57qT+nDuJyF61OygJgoDnnnsOTz/9NDIyMlBVVYWIiAi4u7ubsz4ionYZ0s0bAHDsQgUMBtHsN6ktrtLi2IUKAMDYPl3Num0ish3t7nrLzs6GKIpQKBSIiIjAyJEjTSEpOzvbbAUSEbVHH393qJxlqNLW42yxxuzb33vaeC+5iCBPdPVQmn37RGQb2h2UevTogaKiombLS0pK0KNHjw4VRUTUUU5ymekGuUfOl5l9+4kN3W7j+rI1icietTsotXQPpaqqKqhUqg4VRURkDoO7NQalcrNu12AQkdjQojSur59Zt01EtqXNY5Ti4+MBGMcoPf/883B1vXQfJb1ej6SkJNPNcomIpBTZME7pTzO3KJ0urEJxlRYuznKMaLi6jojsU5uD0uHDhwEYW5SOHj0KhUJhek6hUCAyMhJLliwxX4VERO00pKFFKS23AnV6A5zlHbq9pUlyZgkAYHj3LlA4mWebRGSb2hyUdu7cCQCYN28eVq9eDQ8PD7MXZYvUajXUajX0er3UpRBRK4X7usFD5YTK2nqcKqg03Sy3o5Ib5mYa2YOtSUT2rt3TA6xfvx4JCQlISEhAYWEhDAZDk+fXrVvX4eJsSVxcHOLi4lBRUQEvL/OcbInIsmQyAQOCPJGcWYqTeeYJSqIomlqUosMZlIjsXbvbjF966SVMmTIFCQkJKC4uxsWLF5s8iIhsQf9AY6t3ekGlWbaXU1qDggotnOUChoZ5m2WbRGS72t2i9M4772DDhg247777zFkPEZFZ9Q/0BACczDdPUEo+VwoAGBziBZUzb1tCZO/a3aKk0+kwevRoc9ZCRGR2/RpalE7mVZhle6ZuN45PInII7Q5KDz74ID7//HNz1kJEZHaNQamwUotSja7D20vNKQMATgtA5CDa3fVWW1uL9957Dzt27MCQIUPg7Ozc5PmVK1d2uDgioo5yVzoh1McFOaU1OJlfgdG92j9BZLWuHhmFVQAuTT1ARPat3UHpyJEjpokljx071uS5K83YTUQklX4BnsgprUF6fmWHglJabgUMIhDgqUSAJ+9AQOQI2h2UGudTIiKydQOCPLDjRAFO5nVsQHfjrVAGh3iboSoi6gw4pSwR2T3TgO4OThFw9IIxKLHbjchxtLlFaebMma1a79tvv21zMUREltDb3x0AkFlU1eINvVvjSMM94wYzKBE5jDYHJc5KTUSdTbivGwCgorYeF6vr4OOmuMYrmqusrcPZYg0A4xxKROQY2hyU1q9fb4k6rK66uhoDBgzArFmz8N///lfqcojIglTOcgR7qZBbXovM4ir4uLX90v7juRUQRSDE2wV+7koLVElEtshhxyj95z//wXXXXSd1GURkJT26GluVzhZp2vX6xgkrBwR5mq0mIrJ9DhmUTp8+jZMnT2L69OlSl0JEVtLDzxiUzpW0Lyidapg/qW+Au9lqIiLbZ3NBKTExETNmzEBwcDAEQcDmzZubraNWqxEeHg6VSoWYmBgkJye3aR9LlizBihUrzFQxEXUGPfwaBnQXty8onW64Yq5vgIfZaiIi22dzQUmj0SAyMhJqtfqKz2/atAnx8fFYtmwZDh06hMjISEydOhWFhYWmdaKiojBo0KBmj9zcXHz//ffo27cv+vbta623REQ2oKdf+7veRFHEqQJji1IftigROZR2TzhpKdOnT79ql9jKlSuxYMECzJs3DwCwdu1a/PTTT1i3bh2effZZAEBqamqLr//jjz+wceNGfPXVV6iqqkJdXR08PT3xwgsvXHF9rVYLrVZr+r2iwjw31iQi67q8681gECGTtX6KgKJKLcpr6iATgF5dGZSIHInNtShdjU6nw8GDBxEbG2taJpPJEBsbi/3797dqGytWrEBOTg7OnTuH//73v1iwYEGLIalxfS8vL9MjNDS0w++DiKyvWxcXOMkE1NYZUFBZ26bXNrYmdfd1g8pZbonyiMhGdaqgVFxcDL1ej4CAgCbLAwICkJ+fb5F9Ll26FOXl5aZHTk6ORfZDRJblJJchzMcVAJDZxu63Uw3jk/r4szWJyNHYXNebNc2dO/ea6yiVSiiVSqjVaqjVauj1essXRkQWEe7nhrPFGmSWaDC6d+tvjnu6kAO5iRxVp2pR8vPzg1wuR0FBQZPlBQUFCAwMtOi+4+LikJaWhpSUFIvuh4gsJ7SLCwDgwsWaNr2OA7mJHFenCkoKhQLDhw9HQkKCaZnBYEBCQgJGjRolYWVE1BmENASl820MSmeKGoKSP1uUiByNzXW9VVVVISMjw/R7ZmYmUlNT4ePjg7CwMMTHx2POnDkYMWIERo4ciVWrVkGj0ZiugrMUdr0RdX7duhjHKF0oa31QKqvWoay6DgAQ7udqkbqIyHbZXFA6cOAAJk6caPo9Pj4eADBnzhxs2LABd999N4qKivDCCy8gPz8fUVFR2Lp1a7MB3uYWFxeHuLg4VFRU8MbARJ1UiHdji1J1q1+TVWJc199DCVeFzZ0yicjCbO5f/YQJEyCK4lXXWbRoERYtWmSliojIXnRr6HorrNRCW6+H0unal/pnlRqDUrivm0VrIyLb1KnGKElJrVYjIiIC0dHRUpdCRO3k46aAylkGUQTyylo3l1JWwy1PwnzZ7UbkiBiUWolXvRF1foIgtHmc0rmSxhYlBiUiR8SgREQOpa3jlLJLG1uU2PVG5IgYlIjIoQQ3BKXcVna9sUWJyLExKLUSxygR2YdATxUAoLAV93urrdOjqNJ4U+zQLgxKRI6IQamVOEaJyD4EeCoBAPnl1w5Kjeu4OMvh7eps0bqIyDYxKBGRQwnwMrYo5Vdor7lubsOA72BvFQRBsGhdRGSbGJSIyKGYut4qrt2idMEUlFwsWhMR2S4GJSJyKAENQalEo4O2/uq3JGoc8B3sxaBE5KgYlFqJg7mJ7EMXV2conIynvsJrdL/lskWJyOExKLUSB3MT2QdBEEwDuq915Vtu+aUxSkTkmBiUiMjhBHg0DOgub12LUghblIgcFoMSETmcS1e+tdyiJIqiaYxSEIMSkcNiUCIih9PYonS1rreK2nrU1BkHewd5seuNyFExKLUSB3MT2Q9fdwUAoLRK1+I6jTNye6icoHKWW6UuIrI9DEqtxMHcRPbD180YlEo01w5KXT2UVqmJiGwTgxIRORxfd2P4uWpQqmoISu4MSkSOjEGJiByOT2OLUlXLV72xRYmIAAYlInJAfo1jlNj1RkTXwKBERA6nsUWpWqdHje7KtzFhUCIigEGJiByQu9LJdBuTEs2Vu984RomIAAalVuP0AET2QxCES1e+tTBFQGOLkh9blIgcGoNSK3F6ACL74nuNcUqmrje2KBE5NAYlInJIPm7GAFR8hSvf9AYRpQ1dcv5sUSJyaAxKROSQ/NxablEq1ehgEAFBuDTwm4gcE4MSETkkn6vMzl1WbVzmqXKGk5ynSSJHxjMAETmkLg1BqTEUXa6sps64jquzVWsiItvDoEREDsnTxRiCKmrqmz13saGVycuV3W5Ejo5BiYgckqfKCQBQ3tB6dDm2KBFRIwYlInJIphal2uZBqbzauMzbhUGJyNExKBGRQ/JqCEFXalG62DBuyZtdb0QOj0GplTgzN5F98VQ1jlFquevNm11vRA6PQamVODM3kX1pbFGq1NbDYBCbPNd4JVwXtigROTwGJSJySJ4uxsHcomgMS5crq2aLEhEZMSgRkUNSOsmhcjaeAv/a/XbRFJTYokTk6BiUiMhhNY5T+uuA7vLGwdy86o3I4TEoEZHD8mphioDGFiWOUSIiBiUicliXZue+FJRq6/SoqdMDALw4RonI4TEoEZHDapyd+/LbmDSGJrlMMD1PRI6LQYmIHNaVut4axyt5qJwgCIIkdRGR7WBQIiKH5XmF2bkbpwpwV7I1iYgYlIjIgXldYYyShkGJiC7DoEREDst0G5PaS2OUGJSI6HIOeSYIDw+Hp6cnZDIZunTpgp07d0pdEhFJwK0hDFVdNjN3lVbf5DkicmwOeybYt28f3N3dpS6DiCTkppQDuNSKBABVDQO72aJERAC73ojIgTWGocuDkkbX2KIkl6QmIrItNheUEhMTMWPGDAQHB0MQBGzevLnZOmq1GuHh4VCpVIiJiUFycnKb9iEIAsaPH4/o6Gh89tlnZqqciDobV8WVut6MP7PrjYgAG+x602g0iIyMxPz58zFz5sxmz2/atAnx8fFYu3YtYmJisGrVKkydOhXp6enw9/cHAERFRaG+vr7Za7dt24bg4GDs3bsXISEhyMvLQ2xsLAYPHowhQ4ZcsR6tVgutVmv6vaKiwkzvlIikdqlFSW9a1ti65MGgRESwwaA0ffp0TJ8+vcXnV65ciQULFmDevHkAgLVr1+Knn37CunXr8OyzzwIAUlNTr7qPkJAQAEBQUBBuvPFGHDp0qMWgtGLFCixfvrwd74SIbN2VxyixRYmILrG5rrer0el0OHjwIGJjY03LZDIZYmNjsX///lZtQ6PRoLKyEgBQVVWF3377DQMHDmxx/aVLl6K8vNz0yMnJ6dibICKbYWpR0tVDFEUA7HojoqY61ZmguLgYer0eAQEBTZYHBATg5MmTrdpGQUEBbr/9dgCAXq/HggULEB0d3eL6SqUSSqWy/UUTkc1qDEMGEaitM8BFIYdGx3mUiOgShzsT9OzZE3/++WebX6dWq6FWq6HX66+9MhF1Cq4KOQQBEEVjS5KLQm6aR4lBiYiATtb15ufnB7lcjoKCgibLCwoKEBgYaNF9x8XFIS0tDSkpKRbdDxFZjyAIcFM0nSKgcR4ldr0REdDJgpJCocDw4cORkJBgWmYwGJCQkIBRo0ZJWBkRdVYuCuOA7uqG+ZM0bFEiosvY3JmgqqoKGRkZpt8zMzORmpoKHx8fhIWFIT4+HnPmzMGIESMwcuRIrFq1ChqNxnQVnKWw643IPrk4G4NSbX1jUGoczM0JJ4nIBoPSgQMHMHHiRNPv8fHxAIA5c+Zgw4YNuPvuu1FUVIQXXngB+fn5iIqKwtatW5sN8Da3uLg4xMXFoaKiAl5eXhbdFxFZjyko6fQQRRFVHMxNRJexuTPBhAkTTJfptmTRokVYtGiRlSoiInumcjaOQKip06Nap0fj6cddZXOnRyKSQKcaoyQltVqNiIiIq04lQESdj6qxRanOYOp2kwmXWpqIyLExKLUSr3ojsk+Ng7lr6vSXJptUOEEQBCnLIiIbwaBERA5N5XQpKDVe8capAYioEYMSETm0xhYlbZ0eNXXGoOSqYLcbERkxKLUSxygR2SfTYG6dHrUNQUnJ8UlE1IBBqZU4RonIPjUO5q65rEWpMTwREfFsQEQOzeWyq94aW5Qaxy0RETEoEZFDu7xFSVtnaFjGUyMRGfFs0Eoco0RknxpblLR1etNtTFQco0REDRiUWoljlIjs0+Uzc5u63hiUiKgBgxIRObTLu95qTV1vDEpEZMSgREQOrXEepdomLUo8NRKREc8GROTQLs3MbWCLEhE1w6BERA7N1KKku2weJU4PQEQNGJRaiVe9Edmnxm622no9tOx6I6K/4NmglXjVG5F9Mg3m1nF6ACJqjkGJiByayxWveuOpkYiMeDYgIoemMk04aeA8SkTUDIMSETm0xhYlnd4AjY5BiYiaYlAiIofWeNUbAJRV6wAwKBHRJQxKROTQlE6XToMXNQ1ByYmnRiIy4tmglTg9AJF9EgTBNHi7orYeAFuUiOgSBqVW4vQARPbrr8GIQYmIGjEoEZHDc2kWlHhqJCIjng2IyOGxRYmIWsKgREQOT/mXwdsMSkTUiEGJiByekl1vRNQCng2IyOFdPh2AIAAKOU+NRGTEswERObzLW5RUTnIIgiBhNURkSxiUiMjhXd6ixG43IroczwhE5PCatChxIDcRXYZBqZU4MzeR/WraosSgRESXMCi1EmfmJrJfysu62/46VQAROTaeEYjI4amcLrUiuSrYokRElzAoEZHDu7xFyVXhJGElRGRrGJSIyOEpL2tRcmGLEhFdhkGJiByeqkmLEoMSEV3CoEREDk/JMUpE1AIGJSJyeJe3KLk4c4wSEV3CoEREDo8tSkTUEgYlInJ4HqpLrUgczE1El2NQIiKHF+ilMv3swpm5iegyDEpE5PCCvVxMPzvJBQkrISJb45BBKTMzExMnTkRERAQGDx4MjUYjdUlEJCFvV2fTzxU1dRJWQkS2xiEv75g7dy5efvlljB07FqWlpVAqlVKXREQSEoRLrUj1BlHCSojI1jhci9Lx48fh7OyMsWPHAgB8fHzg5OSQeZGILrNkSl/06uqG+67rLnUpRGRDbC4oJSYmYsaMGQgODoYgCNi8eXOzddRqNcLDw6FSqRATE4Pk5ORWb//06dNwd3fHjBkzMGzYMLzyyitmrJ6IOqtFN/RBwlMT4OvOFmYiusTmmlI0Gg0iIyMxf/58zJw5s9nzmzZtQnx8PNauXYuYmBisWrUKU6dORXp6Ovz9/QEAUVFRqK+vb/babdu2ob6+Hnv27EFqair8/f0xbdo0REdHY/LkyVesR6vVQqvVmn6vqKgw0zslIiIiW2dzQWn69OmYPn16i8+vXLkSCxYswLx58wAAa9euxU8//YR169bh2WefBQCkpqa2+PqQkBCMGDECoaGhAIAbb7wRqampLQalFStWYPny5e18N0RERNSZ2VzX29XodDocPHgQsbGxpmUymQyxsbHYv39/q7YRHR2NwsJCXLx4EQaDAYmJiRgwYECL6y9duhTl5eWmR05OToffBxEREXUONteidDXFxcXQ6/UICAhosjwgIAAnT55s1TacnJzwyiuvYNy4cRBFEVOmTMHNN9/c4vpKpZJXxRERETmoThWUzOVa3XtXolaroVarodfrLVQVERER2ZpO1fXm5+cHuVyOgoKCJssLCgoQGBho0X3HxcUhLS0NKSkpFt0PERER2Y5OFZQUCgWGDx+OhIQE0zKDwYCEhASMGjVKwsqIiIjIHtlc11tVVRUyMjJMv2dmZiI1NRU+Pj4ICwtDfHw85syZgxEjRmDkyJFYtWoVNBqN6So4S2HXGxERkeMRRFG0qfn6d+3ahYkTJzZbPmfOHGzYsAEAsGbNGrz++uvIz89HVFQUVq9ejZiYGKvUV1FRAS8vL5SXl8PT09Mq+yQiIqKOae/3t80FJVvHoERERNT5tPf7u1ONUZKSWq1GREQEoqOjpS6FiIiIrIQtSm3EFiUiIqLOhy1KRERERGbGoERERETUApubHsBWNU4PUF9fD8DYhEdERESdQ+P3dltHHHGMUhudP38eoaGhUpdBRERE7ZCTk4Nu3bq1en0GpTYyGAzIzc2Fh4cHBEEw23YrKioQGhqKnJwchx8kzs/CiJ+DET+HS/hZGPFzMOLnYNTaz0EURVRWViI4OBgyWetHHrHrrY1kMlmbkmhbeXp6OvQBfzl+Fkb8HIz4OVzCz8KIn4MRPwej1nwOXl5ebd4uB3MTERERtYBBiYiIiKgFDEo2QqlUYtmyZVAqlVKXIjl+Fkb8HIz4OVzCz8KIn4MRPwcjS38OHMxNRERE1AK2KBERERG1gEGJiIiIqAUMSkREREQtYFAiIiIiagGDko1Qq9UIDw+HSqVCTEwMkpOTpS7JolasWIHo6Gh4eHjA398ft912G9LT05usM2HCBAiC0OTxyCOPSFSxZbz44ovN3mP//v1Nz9fW1iIuLg6+vr5wd3fHHXfcgYKCAgkrtpzw8PBmn4UgCIiLiwNgv8dDYmIiZsyYgeDgYAiCgM2bNzd5XhRFvPDCCwgKCoKLiwtiY2Nx+vTpJuuUlpZi9uzZ8PT0hLe3Nx544AFUVVVZ8V103NU+h7q6OjzzzDMYPHgw3NzcEBwcjPvvvx+5ublNtnGlY+jVV1+18jvpuGsdE3Pnzm32PqdNm9ZkHXs/JgBc8XwhCAJef/110zrmOCYYlGzApk2bEB8fj2XLluHQoUOIjIzE1KlTUVhYKHVpFrN7927ExcXhjz/+wPbt21FXV4cpU6ZAo9E0WW/BggXIy8szPV577TWJKracgQMHNnmPe/fuNT335JNP4scff8RXX32F3bt3Izc3FzNnzpSwWstJSUlp8jls374dADBr1izTOvZ4PGg0GkRGRkKtVl/x+ddeew2rV6/G2rVrkZSUBDc3N0ydOhW1tbWmdWbPno3jx49j+/bt2LJlCxITE/HQQw9Z6y2YxdU+h+rqahw6dAjPP/88Dh06hG+//Rbp6em45ZZbmq370ksvNTlGFi9ebI3yzepaxwQATJs2rcn7/OKLL5o8b+/HBIAm7z8vLw/r1q2DIAi44447mqzX4WNCJMmNHDlSjIuLM/2u1+vF4OBgccWKFRJWZV2FhYUiAHH37t2mZePHjxcff/xx6YqygmXLlomRkZFXfK6srEx0dnYWv/rqK9OyEydOiADE/fv3W6lC6Tz++ONir169RIPBIIqiYxwPAMTvvvvO9LvBYBADAwPF119/3bSsrKxMVCqV4hdffCGKoiimpaWJAMSUlBTTOr/88osoCIJ44cIFq9VuTn/9HK4kOTlZBCBmZWWZlnXv3l188803LVuclV3ps5gzZ4546623tvgaRz0mbr31VvGGG25osswcxwRblCSm0+lw8OBBxMbGmpbJZDLExsZi//79ElZmXeXl5QAAHx+fJss/++wz+Pn5YdCgQVi6dCmqq6ulKM+iTp8+jeDgYPTs2ROzZ89GdnY2AODgwYOoq6trcmz0798fYWFhdn9s6HQ6fPrpp5g/f36Tm087wvFwuczMTOTn5zc5Bry8vBATE2M6Bvbv3w9vb2+MGDHCtE5sbCxkMhmSkpKsXrO1lJeXQxAEeHt7N1n+6quvwtfXF0OHDsXrr7+O+vp6aQq0sF27dsHf3x/9+vXDwoULUVJSYnrOEY+JgoIC/PTTT3jggQeaPdfRY4I3xZVYcXEx9Ho9AgICmiwPCAjAyZMnJarKugwGA5544gmMGTMGgwYNMi2/99570b17dwQHB+PIkSN45plnkJ6ejm+//VbCas0rJiYGGzZsQL9+/ZCXl4fly5dj7NixOHbsGPLz86FQKJp9EQQEBCA/P1+agq1k8+bNKCsrw9y5c03LHOF4+KvGv+crnR8an8vPz4e/v3+T552cnODj42O3x0ltbS2eeeYZ3HPPPU1ugvrYY49h2LBh8PHxwb59+7B06VLk5eVh5cqVElZrftOmTcPMmTPRo0cPnDlzBv/85z8xffp07N+/H3K53CGPiY8++ggeHh7NhiaY45hgUCLJxcXF4dixY03G5gBo0p8+ePBgBAUFYdKkSThz5gx69epl7TItYvr06aafhwwZgpiYGHTv3h1ffvklXFxcJKxMWh9++CGmT5+O4OBg0zJHOB7o2urq6nDXXXdBFEW88847TZ6Lj483/TxkyBAoFAo8/PDDWLFihV3d5uNvf/ub6efBgwdjyJAh6NWrF3bt2oVJkyZJWJl01q1bh9mzZ0OlUjVZbo5jgl1vEvPz84NcLm92JVNBQQECAwMlqsp6Fi1ahC1btmDnzp3o1q3bVdeNiYkBAGRkZFijNEl4e3ujb9++yMjIQGBgIHQ6HcrKypqsY+/HRlZWFnbs2IEHH3zwqus5wvHQ+Pd8tfNDYGBgsws/6uvrUVpaanfHSWNIysrKwvbt25u0Jl1JTEwM6uvrce7cOesUKJGePXvCz8/P9G/BkY4JANizZw/S09Ovec4A2ndMMChJTKFQYPjw4UhISDAtMxgMSEhIwKhRoySszLJEUcSiRYvw3Xff4bfffkOPHj2u+ZrU1FQAQFBQkIWrk05VVRXOnDmDoKAgDB8+HM7Ozk2OjfT0dGRnZ9v1sbF+/Xr4+/vjpptuuup6jnA89OjRA4GBgU2OgYqKCiQlJZmOgVGjRqGsrAwHDx40rfPbb7/BYDCYwqQ9aAxJp0+fxo4dO+Dr63vN16SmpkImkzXrhrI358+fR0lJienfgqMcE40+/PBDDB8+HJGRkddct13HRIeGgpNZbNy4UVQqleKGDRvEtLQ08aGHHhK9vb3F/Px8qUuzmIULF4peXl7irl27xLy8PNOjurpaFEVRzMjIEF966SXxwIEDYmZmpvj999+LPXv2FMeNGydx5eb11FNPibt27RIzMzPF33//XYyNjRX9/PzEwsJCURRF8ZFHHhHDwsLE3377TTxw4IA4atQocdSoURJXbTl6vV4MCwsTn3nmmSbL7fl4qKysFA8fPiwePnxYBCCuXLlSPHz4sOlqrldffVX09vYWv//+e/HIkSPirbfeKvbo0UOsqakxbWPatGni0KFDxaSkJHHv3r1inz59xHvuuUeqt9QuV/scdDqdeMstt4jdunUTU1NTm5wztFqtKIqiuG/fPvHNN98UU1NTxTNnzoiffvqp2LVrV/H++++X+J213dU+i8rKSnHJkiXi/v37xczMTHHHjh3isGHDxD59+oi1tbWmbdj7MdGovLxcdHV1Fd95551mrzfXMcGgZCPefvttMSwsTFQoFOLIkSPFP/74Q+qSLArAFR/r168XRVEUs7OzxXHjxok+Pj6iUqkUe/fuLT799NNieXm5tIWb2d133y0GBQWJCoVCDAkJEe+++24xIyPD9HxNTY346KOPil26dBFdXV3F22+/XczLy5OwYsv69ddfRQBienp6k+X2fDzs3Lnziv8W5syZI4qicYqA559/XgwICBCVSqU4adKkZp9PSUmJeM8994ju7u6ip6enOG/ePLGyslKCd9N+V/scMjMzWzxn7Ny5UxRFUTx48KAYExMjenl5iSqVShwwYID4yiuvNAkPncXVPovq6mpxypQpYteuXUVnZ2exe/fu4oIFC5r9x9rej4lG7777ruji4iKWlZU1e725jglBFEWx9e1PRERERI6DY5SIiIiIWsCgRERERNQCBiUiIiKiFjAoEREREbWAQYmIiIioBQxKRERERC1gUCIiIiJqAYMSERERUQsYlIio09i1axcEQWh2o2BrSUhIwIABA6DX682yvWeffRaLFy82y7aIyDI4MzcR2aQJEyYgKioKq1atMi3T6XQoLS1FQEAABEGwek3Dhw9HfHw8Zs+ebZbtFRcXo2fPnkhNTUXPnj3Nsk0iMi+2KBFRp6FQKBAYGChJSNq7dy/OnDmDO+64w2zb9PPzw9SpU/HOO++YbZtEZF4MSkRkc+bOnYvdu3fjrbfegiAIEAQB586da9b1tmHDBnh7e2PLli3o168fXF1dceedd6K6uhofffQRwsPD0aVLFzz22GNNusu0Wi2WLFmCkJAQuLm5ISYmBrt27bpqTRs3bsTkyZOhUqlMy1588UVERUXhk08+QXh4OLy8vPC3v/0NlZWVpnW+/vprDB48GC4uLvD19UVsbCw0Go3p+RkzZmDjxo3m+eCIyOwYlIjI5rz11lsYNWoUFixYgLy8POTl5SE0NPSK61ZXV2P16tXYuHEjtm7dil27duH222/Hzz//jJ9//hmffPIJ3n33XXz99dem1yxatAj79+/Hxo0bceTIEcyaNQvTpk3D6dOnW6xpz549GDFiRLPlZ86cwebNm7FlyxZs2bIFu3fvxquvvgoAyMvLwz333IP58+fjxIkT2LVrF2bOnInLRzyMHDkS58+fx7lz59r5aRGRJTlJXQAR0V95eXlBoVDA1dUVgYGBV123rq4O77zzDnr16gUAuPPOO/HJJ5+goKAA7u7uiIiIwMSJE7Fz507cfffdyM7Oxvr165GdnY3g4GAAwJIlS7B161asX78er7zyyhX3k5WVZVr/cgaDARs2bICHhwcA4L777kNCQgL+85//IC8vD/X19Zg5cya6d+8OABg8eHCT1zduMysrC+Hh4a3/kIjIKhiUiKhTc3V1NYUkAAgICEB4eDjc3d2bLCssLAQAHD16FHq9Hn379m2yHa1WC19f3xb3U1NT06TbrVF4eLgpJAFAUFCQaV+RkZGYNGkSBg8ejKlTp2LKlCm488470aVLF9P6Li4uAIwtY0RkexiUiKhTc3Z2bvK7IAhXXGYwGAAAVVVVkMvlOHjwIORyeZP1Lg9Xf+Xn54eLFy+2av+N+5LL5di+fTv27duHbdu24e2338Zzzz2HpKQk9OjRAwBQWloKAOjatWtr3i4RWRnHKBGRTVIoFGabr+hyQ4cOhV6vR2FhIXr37t3kcbVuvqFDhyItLa3N+xMEAWPGjMHy5ctx+PBhKBQKfPfdd6bnjx07BmdnZwwcOLBd74eILIstSkRkk8LDw5GUlIRz587B3d0dPj4+Ztlu3759MXv2bNx///144403MHToUBQVFSEhIQFDhgzBTTfddMXXTZ06FR999FGb9pWUlISEhARMmTIF/v7+SEpKQlFREQYMGGBaZ8+ePRg7dqypC46IbAtblIjIJi1ZsgRyuRwRERHo2rUrsrOzzbbt9evX4/7778dTTz2Ffv364bbbbkNKSgrCwsJafM3s2bNx/PhxpKent3o/np6eSExMxI033oi+ffviX//6F9544w1Mnz7dtM7GjRuxYMGCDr0fIrIczsxNRNRKTz/9NCoqKvDuu++aZXu//PILnnrqKRw5cgROTmzgJ7JFbFEiImql5557Dt27dzcN1u4ojUaD9evXMyQR2TC2KBERERG1gC1KRERERC1gUCIiIiJqAYMSERERUQsYlIiIiIhawKBERERE1AIGJSIiIqIWMCgRERERtYBBiYiIiKgFDEpERERELfh/VsSv8CuawCwAAAAASUVORK5CYII=", "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/qr1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAABW0klEQVR4nO3deVzUdf4H8NfMwHAz3KfcnigCIuCRlkleZRq5W64Zaelua+1htmvbZtdW++swO0irTa220g51d61MJc9UPMEDRS4B5b5vBma+vz+AUURggBm+c7yejwePR8zx/b6/jV948TklgiAIICIiIiJIxS6AiIiIyFAwGBERERG1YzAiIiIiasdgRERERNSOwYiIiIioHYMRERERUTsGIyIiIqJ2FmIXYGzUajUKCgrg4OAAiUQidjlERESkBUEQUFtbCx8fH0il3bcLMRj1UUFBAfz8/MQug4iIiPohPz8fQ4YM6fZ5BqM+cnBwAND2P9bR0VHkaoiIiEgbNTU18PPz0/we7w6DUR91dJ85OjoyGBERERmZ3obBcPA1ERERUTsGIyIiIqJ2DEZERERE7RiMiIiIiNoxGBERERG1YzAiIiIiasdgRERERNSOwYiIiIioHYMRERERUTsGIy0lJiYiNDQU0dHRYpdCREREeiIRBEEQuwhjUlNTA4VCgerqam4JQkREZCS0/f3NFiMiIiKidgxGREREZBCUrWqcv1YNZatatBosRDszERERma1WlRqXi+tw7loVzl6txrlr1bhUWAulSo3/PjEZY4c4iVIXgxERERHplSAIuFLegNT8KqRebQtCFwqq0dTStWXIydYSJTXNIlTZhsGIiIiIdKq0tlkTglLy24JQdWNLl9fZW1kgzFeBsUMUCBuiQPgQJwxxtoFEIhGh6jYMRkRERNRvTS0qnLtWjZS8thCUkl+Fa1WNXV4nt5BitI8jwoc4YewQBcYOcUKwmx2kUvFC0K0wGBEREZFWBEFAbnkDTudV4kx7ELpYWINWdeeVfyQSYJiHPSL8nDB2iBMi/JwwwssBljLDn/PFYERERES3VN/citT8KpzOq8TpvCqcyatEZUPXLjF3BytE+jkhwt9JE4bsrYwzYhhn1URERKRTgiDgamUjTuVWar4uFdXgpsYgyGVSjPF1xDh/Z0T4OyHS3xk+CmtRxwXpEoMRERGRGWpRqZFWUIMTVyo0QaiktutsMB+FNSIDnDHO3xnj/J0Q6uMIKwuZCBUPDgYjIiIiM1DX3IrTuZU4eaUCJ65UIiW/Co0tqk6vsZRJMNpHgagAZ0S1hyEvhbVIFYuDwYiIiMgEldc148SVChzPqcSJKxW4UFDdpVtMYWOpCUHjA5wR7ucEa0vTbQ3SBoMRERGRCSiuacKx7HIcz6lAck4FMkvqurzG18kGMUEuGB/ojOhAFwx1tze46fJiYzAiIiIyQoXVjUjOrsCx7HIcyy7HlfKGLq8Z7mmP6EAXxAS5IDrQBT5ONiJUalwYjIiIiIxASU0TjraHoKNZXYOQVAKE+jgiJtAVscFtQcjFTi5StcaLwYiIiMgAVdYrcSy7HEeyynEkqwxZpfWdnpdKgDG+CkwIdsWEYBeMD3SBo7WlSNWaDgYjIiIiA9CgbMXxnAr8klmGXzLLcbGoBsINg6UlEiDU2xETg10xMcQV0UEMQvrAYERERCSCVpUaqVer8UtmGQ5nluFMXiVaVJ2njQ33tMekEDdNq5CTLbvG9I3BiIiIaBAIgoAr5Q04lFGKQxllOJZVjtrm1k6v8XWywW1D3TBpaFurkIeDea0hZAgYjIiIiPSkuqEFR7LKcDCjDIcySnG1svOu8wobS0wKccVtw9xw21A3+LvYmszWGsaKwYiIiEhHVGoBZ69W4eDlMhy4XIKU/KpOiyrKZVJEBThrgtAYXwVkXEfIoDAYERERDUBpbTMOXi7F/sulOJRRiqqbdp8f6mGPKcPcMHWYO2KDXWAr569eQ8ZPh4iIqA9UagEp+VXYn16C/emlOHetutPzDtYWuG2oG24f7o4pw93hy0UVjQqDERERUS+qGpQ4cLkU+9NLceByKSrqlZ2eH+PriDuGe+D2Ee6I9HOChUwqUqU0UAxGRERENxEEAVmldUi6WIKkSyU4lVsJ1Q2DhRytLTBluDumjfDA1OFunD1mQswuGFVVVSEuLg6tra1obW3FH//4RyxbtkzssoiISGStKjVOXKnE3ovF2HuxGLk3bbkxwtMB00Z6YNoId0QFOLNVyESZXTBycHDAwYMHYWtri/r6eowZMwbx8fFwdXUVuzQiIhpkdc2tOHi5FHvSivHzpRJUN14fOC2XSTEhxBXTR3rgzpEe8HOxFbFSGixmF4xkMhlsbdv+cTc3N0MQBAiC0Mu7iIjIVJTWNiPpYjF2pxXjcGYZlK1qzXPOtpa4c6Qn4kZ5YMpwd9hbmd2vSbNncJ/4wYMH8cYbb+DUqVMoLCzE9u3bMX/+/E6vSUxMxBtvvIGioiKEh4fjvffeQ0xMjNbnqKqqwu23346MjAy88cYbcHNz0/FVEBGRIcmvaMBPF4rw04UinMyt7LQHWaCrLe4K9cRdoV6ICnDmukJmzuCCUX19PcLDw7F06VLEx8d3eX7r1q1YuXIlNmzYgNjYWKxbtw4zZ85Eeno6PDw8AAARERFobW3t8t7du3fDx8cHTk5OSE1NRXFxMeLj47FgwQJ4enresp7m5mY0Nzdrvq+pqdHRlRIRkT5ll9bhx/NF+PF8Ic5f6/yze+wQBWaEemLGaC8M87DnatOkIREMuB9JIpF0aTGKjY1FdHQ03n//fQCAWq2Gn58fnnzySaxevbrP5/j973+PO++8EwsWLLjl8y+88AJefPHFLo9XV1fD0dGxz+cjIiL9ySiuxQ/n2sLQpaJazeNSCRAd6IJZY7wwY7QX1xYyQzU1NVAoFL3+/ja4FqOeKJVKnDp1Cs8884zmMalUiri4OBw9elSrYxQXF8PW1hYODg6orq7GwYMH8fjjj3f7+meeeQYrV67UfF9TUwM/P7/+XwQREelURnEtvj9XiO/PFiKjpE7zuIVUgslD3TBrjBfuCvWEm72ViFWSsTCqYFRWVgaVStWl28vT0xOXLl3S6hi5ublYvny5ZtD1k08+ibCwsG5fb2VlBSsr3kxERIYkq7QOO1ML8f25Alwuvh6G5DIppgxzw+wwb9w1yhMKW0sRqyRjZFTBSBdiYmKQkpIidhlERNRH+RUN2Hm2EP9LLUBa4fUxQ5YyCaYOc8fdY70RF+oJR2uGIeo/owpGbm5ukMlkKC4u7vR4cXExvLy8RKqKiIj0pbS2Gd+fLcB/UwtwOq9K87iFVILbhrnhnrE+uCvUEwobhiHSDaMKRnK5HFFRUUhKStIMyFar1UhKSsITTzyh13MnJiYiMTERKpVKr+chIjJ3tU0t+OlCMf6Tcg2/ZJahYycOiQSYGOyKueE+mDXaC852cnELJZNkcMGorq4OmZmZmu9zcnKQkpICFxcX+Pv7Y+XKlUhISMD48eMRExODdevWob6+HkuWLNFrXStWrMCKFSs0o9qJiEh3WlRqHLxciu1nrmFPWjGab1h0McLPCfeG++Cesd7wcOSeZKRfBheMTp48iWnTpmm+75gRlpCQgM2bN+OBBx5AaWkp1qxZg6KiIkRERGDXrl3drkNERESGSRAEnL1ajW2nr+K/qQWobLi+HUeIux3mR/ji3ggfBLjaiVglmRuDXsfIEGm7DgIREd1aQVUjtp+5hm2nryKrtF7zuJu9FeZF+GB+hC/G+Dpy0UXSKZNcx4iIiIxTU4sKP10owrenruJwZplmSw5rSylmjvbCfZG+uG2oG3esJ9ExGGmJg6+JiPpGEASk5Ffh65NXsTO1ALXN17dqig1ywf1RQzB7jBccOL2eDAi70vqIXWlERD2rqFdi+5lr2Hoir9Pii75ONlgQNQT3jxsCf1dbESskc8SuNCIiGjRqtYCj2eX48nge9lwohlLVNqvMykKKOWHe+NX4IZgQ5Aopd64nA8dgRERE/VZa24xvT13FlhN5yC1v0Dwe5qvAr6P9cG+4DxdfJKPCYERERH0iCG2tQ18k5+Gn80VobV+B0cHKAvMjffFAtB/G+HK9NzJODEZERKSV6sYWfHfqKr5Izu00zT7S3wkLY/xxz1hv2Mr5a4WMG/8Fa4mz0ojIXF0srMFnR69g+5lraGppGztkK5dhfqQvHooNQKgPJ6KQ6eCstD7irDQiMgctKjV+ulCEz47k4viVCs3jIzwd8NAEf8yP9OU0ezIqnJVGRER9VlGvxFfH8/D50VwU1TQBAGRSCWaN9sLDEwMQE+TCFanJpDEYERERLhXVYNPhK9iRck2zgaubvRV+E+uP38T4w0vBzVvJPDAYERGZKbVawIGMUnxyKAeHM8s0j4f5KrBkciDuHusNKwuZiBUSDT4GIyIiM9PUosL2M9fwyeEcZJa0rUwtlQCzx3hjyeRARAU4s7uMzBaDERGRmahqUOLzo7nYfOQKyuuVAAB7Kws8GO2HhEmB8HPhNh1EDEZa4nR9IjJW+RUN+ORwDraeyEdjS9vPMF8nGyyZHIgHov04u4zoBpyu30ecrk9ExuJSUQ027M/C/84WQtW+OvVoH0csnxqMu8O8YSGTilwh0eDhdH0iIjN18koFPtifhZ8vlWgemzLMDb+dGoLJQ105foioBwxGREQmQBAEHM4sw3s/Z+J4TtuCjBIJMCfMG4/fHsK9y4i0xGBERGTEBEFA0sUSvLcvE6n5VQAAuUyK+6N8sXxqCILc7MQtkMjIMBgRERkhtVrA7rRivJuUgbTCGgCAtaUUC2P88dupIVyQkaifGIyIiIxIRyB6JykDF9sDkZ1chsUTA/HobUFwd7ASuUIi48ZgRERkBAShLRCt23s9ENlbWWDJ5LZA5GQrF7lCItPAYKQlrmNERGIQBAEHLpdi7Z7LOHu1GgADEZE+cR2jPuI6RkQ0WI5ll+Ot3ek4caUSAGArl+GRSYFYPjWYgYioj7iOERGRkbpQUI3Xd6XjwOVSAIDcQoqHJwTgd3eEwM2eY4iI9InBiIjIQOSW1+Ot3Zfx39QCAICFVIIHY/zwxLRhnGVGNEgYjIiIRFZRr8S7SRn497FctLZv3XFvuA+emjEcAa5ch4hoMDEYERGJpKlFhU2/XMEH+zJR29wKAJg63B1/mTmCK1UTiYTBiIhokKnVAv6bWoA3fkrHtapGAECotyOevXsUJg91E7k6IvPGYERENIhO51Xixf+labbv8FZY4+mZIzA/whdSKTd3JRIbgxER0SAoqGrE/+26hP+ktA2stpPL8PtpQ/HobUGwtpSJXB0RdWAwIiLSo6YWFT46mI0P9meiqUUNiQRYMG4Inp45Ah6OnGlGZGgYjIiI9GRvWjFe2pmGvIoGAEB0oDPW3DMaYUM4sJrIUDEYaYlbghCRtq6U1eOlnWn4+VIJAMDT0Qp/mzMK94b7QCLhOCIiQ8YtQfqIW4IQUXeaWlTYcCALH+zPgrJVDUuZBI/eFown7xwKOyv+HUokJm4JQkQ0iA5llOK5Hedxpbyt22zKMDe8cO9ohLjbi1wZEfUFgxER0QCU1DbhHzsvarbx8HCwwvNzR2NOmBe7zYiMEIMREVE/CIKAr0/m45XvL6KmqRVSCfDwxEA8NWM4HKwtxS6PiPqJwYiIqI9yyurxzLazOJZdAQAI81Xg1fvCONuMyAQwGBERaalVpcZHh7Kxbm8GlK1qWFtKsWrGCDwyKRAWMqnY5RGRDjAYGYj65lYcyijDrDFeYpdCRLdwqagGT39zFueuVQNoG1z9yvww+LvailwZEekSg5EBULaqsfzzk/glsxx/mzMSy6eGiF0SEbVrUamxfn8W3vs5Ay0qAQobS6y5JxTx43w5uJrIBDEYGQBLmQRjfBX4JbMcr/5wCVUNLXh65gj+0CUS2cXCGqz6JhUXCmoAAHGjPPHqfWO4lQeRCWMwMgASiQTPzB4FJxs5/m/XJXywPwtVjS14ed4YyLjbNtGgU6kFfHgwC2/vuYwWlQAnW0u8eO9orlxNZAYYjAzI43eEwMnWEn/bfg5fJuehprEFa38dAbkFB3USDZbc8no89XUqTuZWAmhvJYofAw8HthIRmQMGIwOzMMYfDtYW+PPWFOw8W4i65lZseCgK1pYysUsjMmmCIGDLiXy8vDMNDUoV7K0s8PzcUCyIGsJWIiIzwqYILSUmJiI0NBTR0dF6P9c9Y33wr4RoWFtKsT+9FEs3n0CDslXv5yUyVxX1Siz//BSe2XYODUoVYoNc8OMfp+BX4/0YiojMDDeR7aPB3ET2eE4Flmw6jnqlCtGBztj4SDRX1CXSsV8yy/DnrSkoqW2GXCbF0zNH4NHbgiDl+D4ik6Lt72+2GBmwmCAXfP5YLBysLXDiSiUWf3Ic1Y0tYpdFZBKUrWq89sNFPPRJMkpqmxHiboftKyZh2dRghiIiM8ZgZODG+Tvjq2UT4GRriZT8Kvzm42OorFeKXRaRUcsrb8CCDUfw4cFsCALwm1h/7HxyCkb7cEsPInPHYGQExvgqsGX5BLjayXGhoAaL/pWMqgaGI6L+2HW+EHe/dwhnr1ZDYWOJDxdH4dX7wmAj5wQHImIwMhojvRyx9bcT4GZvhbRChiOivmpuVeGF/17A7/59GrVNrRjn74Qf/jgFM0dzGx4iuo7ByIgM9XDAV8tiNS1Hiz85juoGjjki6k1+RQMWrD+KzUeuAAB+OzUYW387Eb5ONuIWRkQGh8HIyAzzdMCXyybAxU6Oc9eqsXhjMgdkE/XgwOVSzH3/MM5dq4aTrSU2PjIez8wZBUsZf/wRUVf8yWCERng54MtlsXC2tcTZq9V4eONx1DVznSOiG6nVAt5LysAjm46jqqEFY4co8P0fpuDOkZ5il0ZEBozByEiN9HLEl+2z1VLzq/Do5hNoVKrELovIIFQ3tmDZZyfx1p7LEIS2FeW/ZtcZEWmBwciIjfJ2xGdLY2BvZYHknAo8/sUpKFvVYpdFJKrMkjrMT/wFSZdKILeQ4vX7x+K1+DBuq0NEWmEwMnJjhzhh4yPXtw/545YzaFUxHJF52pdegvsSf0FOWT18nWzw3e8m4dfRfmKXRURGhMHIBMQEueCjxeMhl0nx4/ki/OW7s1CrudMLmQ9BEPDxwWw8uvkEaptbER3ojP88MRlhQ7hgIxH1DYORiZg63B3v/yYSMqkE205fw0s708Bt8MgcNLeq8NQ3qXjlh4tQC8CD0X744rG2Nb+IiPqKwciEzBjthbd+FQ4A2HzkCtYfyBK5IiL9qqxXYvG/jmPb6WuQSSV4YW4oXosPg9yCP9qIqH8sxC6AdGt+pC/K65V4eWcaXt+VDjd7K/x6PMdYkOm5UlaPJZtPIKesHg5WFvjgoXGYMsxd7LKIyMgxGJmgR28LQmltMzYcyMIz287BxVaOuFCu3UKm4+SVCiz77CQqG1rg62SDTUuiMdzTQeyyiMgEsL3ZRP111gjcP24IVGoBK748jVO5FWKXRKQTO88W4Df/SkZlQwvChyiwfcUkhiIi0hkGIy0lJiYiNDQU0dHRYpeiFYlEgn/eH4Y7R3qguVWNRz89iezSOrHLIhqQzb/k4MmvzkDZqsaMUE9sWT4RHg7WYpdFRCZEInDqUp/U1NRAoVCguroajo6OYpfTq0alCg9+fAyp+VUIdLXFtt9PhoudXOyyiPpEEAS8tfsy3t+XCQBImBiANXNHQyaViFwZERkLbX9/s8XIxNnIZfjXw+MxxNkGV8obsOyzk2hq4dYhZDxaVWqs/u6cJhStmjEcL9zLUERE+sFgZAbcHayw6ZFoOFhb4FRuJVZ9k8oFIMkoNLWo8PgXp7H1ZD6kEuC1+DA8cecwSCQMRUSkHwxGZmKYpwM+fCgKFlIJdp4txFt70sUuiahH9c2tWLr5BPakFUNuIcX6h6KwMMZf7LKIyMQxGJmRSUPd8Fp8GAAgcV8WvjmZL3JFRLdW09SChzcex5GsctjJZfhsaQxmjvYSuywiMgMMRmbmV+P98MS0oQCAZ7efx6ncSpErIuqssl6JRR8n41RuJRytLfDvx2IxIdhV7LKIyEwwGJmhlXcNx4xQTyhVavzu36dQVN0kdklEAICS2iY8+NExnLtWDRc7Ob5aPgGR/s5il0VEZoTByAxJpRK8/UAERng6oLS2Gcs/50w1El9JbRMWfnQM6cW18HCwwte/nYDRPgqxyyIiM8NgZKbsrCzw8cPj4WRribNXq7H6u7PgklYkltLaZvzm42RkldbDR2GNr387EUM9uJo1EQ0+BiMz5u9qiw9+Mw4yqQQ7Ugrw0cFssUsiM1RW14xF/zqGzJI6eCus8dXyCQh0sxO7LCIyUwxGZm7SUDesuScUAPB/uy7hSFaZyBWROSmva8aij5NxubgOno5W+GrZBAS4MhQRkXgYjAgPTwzAgqghUAvAH746w8HYNCgq65VY9K9kzZiiLcsnsqWIiETHYESQSCR4ed4YjPRyQFmdEk98eRotKrXYZZEJq2tuxSObjuNSUS3cHazw1fIJCGIoIiIDwGBEANr2VNvwUBQcrCxwMrcS//zxktglkYlqalFh2acnkXq1Gs62lvjisViEuNuLXRYREQAGI7pBoJsd3vx1OADgk8M5+P5socgVkalpVanx5FdncDS7bUXrzUtiMNyTs8+IyHAwGFEnM0d74be3BwMA/vJtKrJK60SuiEyFWi3gL9+e1ex99q+EaIT7OYldFhFRJwxG1MXTM0YgNsgF9UoVnvzyDJpbufgjDYwgCHj5+zRsO3MNMqkEH/xmHCaGcJsPIjI8AwpGLS0tyM/PR3p6OioqKnRVE4nMQibFuwsj4WxribTCGvzfj+lil0RG7uND2dj0yxUAwFu/CkdcqKe4BRERdaPPwai2thbr16/H7bffDkdHRwQGBmLUqFFwd3dHQEAAli1bhhMnTuijVhpEno7WePNXbeONNv6Sg32XSkSuiIzVf1ML8OoPbYP5/373KMyP9BW5IiKi7vUpGK1duxaBgYHYtGkT4uLisGPHDqSkpODy5cs4evQonn/+ebS2tmLGjBmYNWsWMjIy9FU3DYLpozzxyKRAAMCqb1JRUsP1jahvjmWXY9XXqQCAJZMD8ehtQSJXRETUM4nQhw2yFi5ciL///e8YPXp0j69rbm7Gpk2bIJfLsXTp0gEXaUhqamqgUChQXV0NR0dHscvRu6YWFe774AguFtbgtqFu+GxpDKRSidhlkRHIKK7F/euPoKapFbNGeyFxUdv2M0REYtD293efghGZXzACgMySWtzz3mE0tajx11kj8fgdIWKXRAaupKYJ931wBNeqGhEV4IwvHouFtaVM7LKIyIxp+/u734OvGxsb0dDQoPk+NzcX69atw08//dTfQ5KBGurhgBfmtrUSrt2TjgsF1SJXRIasqUWFZZ+fwrWqRgS72eFfD49nKCIio9HvYDRv3jx89tlnAICqqirExsbirbfewvz587F+/XqdFUiG4YFoP8wI9USLSsBTX6dyCj/dkiAIWP3dWaTmV8HJ1hIbH4mGs51c7LKIiLTW72B0+vRpTJkyBQDw7bffwtPTE7m5ufjss8/w7rvv6qxAMgwSiQSv3BcGFzs5LhXV4t0kDqynrj7Yn4UdKQWwkErwwaJx3BSWiIxOv4NRQ0MDHBzalvLfvXs34uPjIZVKMWHCBOTm5uqsQDIc7g5W+Mf8MQCA9fuzcCavUuSKyJDsvlCEN35qW/Pq+XtHY1KIm8gVERH1Xb+D0dChQ7Fjxw7k5+fjp59+wowZMwAAJSUlRjEouaGhAQEBAVi1apXYpRiVOWHemBfhA7UAPPV1KhqV7FIj4GJhDf60NQUAsHhCABZPCBC3ICKifup3MFqzZg1WrVqFwMBAxMbGYuLEiQDaWo8iIyN1VqC+vPLKK5gwYYLYZRill+4dA09HK2SX1WtaCMh8VTUoseyzk2hQqjApxBVr5oaKXRIRUb/1OxgtWLAAeXl5OHnyJHbt2qV5fPr06Xj77bd1Upy+ZGRk4NKlS5g9e7bYpRglha0l/nn/WABtq2Ifyy4XuSISi1ot4M9bU3C1shF+Ljb4YNE4WMq4BSMRGa8B/QTz8vJCZGQkpNLrh4mJicHIkSP7fcyDBw9i7ty58PHxgUQiwY4dO7q8JjExEYGBgbC2tkZsbCyOHz/ep3OsWrUKr732Wr9rJGDaCA88GO0HAPjbtnNoamGXmjlK3JeJfemlsLKQYv2iKDjZcgYaERk3i4G8OSkpCUlJSSgpKYFare703MaNG/t1zPr6eoSHh2Pp0qWIj4/v8vzWrVuxcuVKbNiwAbGxsVi3bh1mzpyJ9PR0eHh4AAAiIiLQ2tra5b27d+/GiRMnMHz4cAwfPhxHjhzptZ7m5mY0Nzdrvq+pqenXdZmiZ+aMws+XSpBdVo/EfZl4asYIsUuiQXQooxRr914GALw8fwzG+CpEroiIaOD6vfL1iy++iJdeegnjx4+Ht7c3JJLOS/1v37594MVJJNi+fTvmz5+veSw2NhbR0dF4//33AQBqtRp+fn548sknsXr16l6P+cwzz+Df//43ZDIZ6urq0NLSgqeeegpr1qy55etfeOEFvPjii10eN6eVr3vy47lCPP7FaVhIJfj+D1MwwstB7JJoEFyrasQ97x5CZUMLHoz203StEhEZKr1vCeLt7Y3XX38dixcv7neRvbk5GCmVStja2uLbb7/tFJYSEhJQVVWF//znP306/ubNm3H+/Hm8+eab3b7mVi1Gfn5+DEbtBEHA8s9PYU9aMSL9nfDt7yZxPywT19yqwq8/PIbU/CqM8XXEt7+bxJWticjg6X1LEKVSiUmTJvX37f1SVlYGlUoFT0/PTo97enqiqKhIL+e0srKCo6Njpy+6TiKR4KV5o2FvZYEzeVX4IplrWJm613elIzW/CgobS6xfFMVQREQmpd/B6LHHHsOXX36py1oG3SOPPNJjaxFpx1thg7/Mahtf9PqudBRWN4pcEenL/vQSfHI4BwDw1q/C4ediK3JFRES61e/B101NTfjoo4+wd+9ejB07FpaWlp2eX7t27YCLu5mbmxtkMhmKi4s7PV5cXAwvLy+dn4+091BsAHacuYbTeVVY858L+Pjh8WKXRDpWWtuMVd+kAgASJgYgLtSzl3cQERmffrcYnT17FhEREZBKpTh//jzOnDmj+UpJSdFhidfJ5XJERUUhKSlJ85harUZSUpJmgUl9SUxMRGhoKKKjo/V6HmMllUrwWvxYWEgl2JNWjH2XSsQuiXRIrRaw6ptUlNUpMcLTAc/MGSV2SUREetHvFqN9+/bpsg6Nuro6ZGZmar7PyclBSkoKXFxc4O/vj5UrVyIhIQHjx49HTEwM1q1bh/r6eixZskQv9XRYsWIFVqxYoRm8RV2N8HLAksmB+PhQDl783wVMGuoKKwuOPzEFG3/JwYHLbesVvfebSI4rIiKT1ecWozVr1uDUqVP6qAUAcPLkSURGRmq2FVm5ciUiIyM10+kfeOABvPnmm1izZg0iIiKQkpKCXbt2dRmQTeL4w/RhcHewwpXyBs1YFDJu569V4/92XQIA/P2eUAz35JIMRGS6+jxdf+nSpdi5cyfkcjnmzp2Le++9F9OnT4dcbh4r3mo73c+cbTt9FSu/ToWtXIafn7oDXgprsUuifmpqUeGe9w4js6QOM0I98eHiqC5rlhERGQO9TdffuHEjioqK8NVXX8HBwQF/+tOf4Obmhvvvvx+fffYZKioqBlQ4Gb/7In0RFeCMBqUKr/5wUexyaADe3nMZmSV1cHewwv/dP5ahiIhMXr8GX0ulUkyZMgWvv/460tPTkZycjNjYWHz44Yfw8fHB1KlT8eabb+LatWu6rlc0HHytPYlEghfvHQ2JBPhvagGSucmsUTqdV4mPD2UDAF67LwzOdubRKkxE5q3PwWj58uWoq6vr9NioUaPwl7/8Bb/88gvy8/ORkJCAQ4cO4auvvtJZoWJbsWIF0tLScOLECbFLMQpjfBVYGOMPAHj+vxfQqlL38g4yJE0tKqz6JhVqAYiP9OXUfCIyG30eYySTyVBYWKjZsNXccIyR9irrlbjjzf2obmzBy/PHYPGEALFLIi29+sNFfHQwGx4OVtjz59uhsLXs/U1ERAZMb2OM+rm1GpkhZzs5Vt41HADwzt7LqGtuFbki0sap3Bu60OLDGIqIyKz0e4FHIm38JtYfQW52KKtT4uOD2WKXQ71oalHh6W9SIQhA/DhfTB/FLjQiMi99DkYSiYQzU0hrljIpnp7Zto/ax4eyUVLTJHJF1JN3kzKQXVYPT0crPH/PaLHLISIadH1e+VoQBDzyyCOwsrLq8XXbtm3rd1GGKDExEYmJiVCpVGKXYnRmj/FChJ8TUvKrsC4pA6/eFyZ2SXQLmSW1mi60l+aNYRcaEZmlPrcYJSQkwMPDAwqFoscvU8NZaf0nkUjwt/a9tbaeyEdmSV0v76DBJggC/r7jPFpUAuJGeWAGZ6ERkZnqc4vRpk2b9FEHmbiYIBfEjfLE3ovFeH3XJXz08HixS6IbbD9zDceyK2BtKcXzc0ezu5yIzBYHX9OgWT17BKQSYHdaMU5e4QrphqKqQYlXvm9bofwP04fBz8VW5IqIiMTTp2CUl5fXp4Ob0srXNHBDPRzwQHTboo+v/nCRSz8YiNd/Skd5vRLDPOzx2G3BYpdDRCSqPgWj6Oho/Pa3v+1xnE11dTU+/vhjjBkzBt99992ACyTT8ue4YbC2lOJ0XhUOXC4VuxyzdzqvEl8db/uD5x/zx0BuwUZkIjJvfRpjlJaWhldeeQV33XUXrK2tERUVBR8fH1hbW6OyshJpaWm4cOECxo0bh9dffx1z5szRV91kpDwcrbF4QgA+PpSDdXszcPtwd45nEYlKLeC5HechCMD944YgNthV7JKIiETXpz8PXV1dsXbtWhQWFuL999/HsGHDUFZWhoyMDADAokWLcOrUKRw9etTkQhE3kdWd5VNDYG0pRUo+W43E9N3pq7hQUAMHaws8M2ek2OUQERmEPu+VZu64V5pudOzFFe7nhB2/n8RWo0HWoGzFHW/sR0ltM56dMwrLpnJsERGZNr3tlUakC8unBsPaUorU/CrsZ6vRoPvwQDZKapvh72KLhydxc18iog4MRiQKN3srPDwxEACwbm8GZ6gNoqLqJnx4MAsAsHr2SFhZyESuiIjIcDAYkWg6tRqls9VosLy5Ox1NLWqMD3DG7DFeYpdDRGRQGIxINJ1bjS6z1WgQnL9Wje9OXwUAPHv3KI7tIiK6Sb+DUUJCAg4ePKjLWsgMLZ8aDBtLGVKvVrPVSM8EQcAr31+EIADzInwQ6e8sdklERAan38GouroacXFxGDZsGF599VWuck394mZvhcUT2wb/rj+QJXI1pu3nSyU4ml0OuYUUT88cIXY5REQGqd/BaMeOHbh27Roef/xxbN26FYGBgZg9eza+/fZbtLS06LJGg8B1jPRn6eQgWMokOJ5TgZT8KrHLMUlqtYA3d18GACyZHIghztwPjYjoVgY0xsjd3R0rV65EamoqkpOTMXToUCxevBg+Pj7485//rFn40RSsWLECaWlpPW6HQv3jpbDGveG+AICPDrLVSB92pxXhYmEN7K0s8LupIWKXQ0RksHQy+LqwsBB79uzBnj17IJPJMGfOHJw7dw6hoaF4++23dXEKMnHL2xcY3HW+CLnl9SJXY1rUagFv72n7I2Xp5EA428lFroiIyHD1Oxi1tLTgu+++wz333IOAgAB88803+NOf/oSCggJ8+umn2Lt3L77++mu89NJLuqyXTNQILwfcPtwdagH45HCO2OWYlO/PFSK9uBYO1hZ49DaucE1E1JM+bSJ7I29vb6jVaixcuBDHjx9HREREl9dMmzYNTk5OAyiPzMlvpwbjwOVSfH0yH3+OG86WDR1QqQWs29s2tuix24KhsLUUuSIiIsPW7xajt99+GwUFBUhMTLxlKAIAJycn5OTwr3/SzsQQV4zxdURTixqfH8sVuxyT8L/UAmSV1kNhY4kltwWKXQ4RkcHrdzC6/fbbYWVl1eVxQRCQl5c3oKLIPEkkEiyb0tbV8+mRK2hqUYlckXFrVanxTlLb2KLlU4PhaM3WIiKi3vQ7GAUFBaG0tOuCfBUVFQgKChpQUWS+7g7zhq+TDcrrldh2mmtjDcT2M9eQU1YPFzs5EiYFil0OEZFR6HcwEgThltsJ1NXVwdraekBFkfmykEnx6G1twfpfh7KhVnObkP5oVanx3s+ZANrGbtlb9Xs4IRGRWenzT8uVK1cCaOv2eO6552Bre32hOJVKheTk5G7HHBFp44FoP7y99zKyy+pxOLMMU4e7i12S0dl1oQh5FQ1wtrXUrCxORES963MwOnPmDIC2FqNz585BLr8+c0gulyM8PByrVq3SXYVkduysLHD/uCHYfOQK/n0sl8GojwRBwEcHswEAiycGwlbO1iIiIm31+Sfmvn37AABLlizBO++8A0dHR50XZYgSExORmJgIlYoDggfDolh/bD5yBXsvFqOwuhHeChuxSzIax7IrcPZqNawspEhgaxERUZ/0e4zRpk2bzCYUAdwSZLAN83RAbJAL1ALwVTJnOfZFx7YqC6KGwNW+68xRIiLqXp9ajFauXImXX34ZdnZ2mrFG3Vm7du2ACiN6aEIAknMqsOVEPp6cPgyWMp3sYGPSLhfXYl96KSQS4LEpXOWaiKiv+hSMzpw5g5aWFs1/d+dWs9WI+mrmaC+42VuhpLYZe9KKMSfMW+ySDF7H2KKZoV4IcrMTuRoiIuPTp2DUMb7o5v8m0ge5hRQPRvvh/X2Z+PexXAajXhRVN+E/KW1rPy2/na1FRET90e++icbGRjQ0NGi+z83Nxbp167B7926dFEYEAAtj/SGVAEeyypFZUid2OQZt05EctKgERAc6Y5y/s9jlEBEZpX4Ho3nz5uGzzz4DAFRVVSEmJgZvvfUW5s2bh/Xr1+usQDJvvk42uHOkBwDgi2Tun9ad2qYWfHmsbZD68qkhIldDRGS8+h2MTp8+jSlTpgAAvv32W3h5eSE3NxefffYZ3n33XZ0VSLRoQtuU8+9OXUWjkssl3Mo3J6+itrkVIe52mN4eJImIqO/6HYwaGhrg4OAAANi9ezfi4+MhlUoxYcIE5ObyL3vSnduHucPPxQY1Ta34X2qB2OUYHEEQNK1pj0wOglTKyQ9ERP3V72A0dOhQ7NixA/n5+fjpp58wY8YMAEBJSYlZrW9E+ieVSvBgtD8A4NvTV0WuxvAk51Qgq7QednIZ7ov0FbscIiKj1u9gtGbNGqxatQqBgYGIjY3FxIkTAbS1HkVGRuqsQCIAuC/SFxIJcDynAvkVDb2/wYx80b4A5rxIX24WS0Q0QP0ORgsWLEBeXh5OnjyJXbt2aR6fPn063n77bZ0UR9TBx8kGE4NdAQA7zlwTuRrDUVrbjF3nCwEAv4nxF7kaIiLjN6ClhL28vBAZGQmp9PphYmJiMHLkyAEXRnSz+HFDAADbzlyDIAgiV2MYvjmVjxaVgAg/J4zxVYhdDhGR0RtQu3tSUhKSkpJQUlICtVrd6bmNGzcOqDBDw01kxTdrjBf+vuMccsrqkZJfhUgzX6tHrRbwZXs32qJYthYREelCv1uMXnzxRcyYMQNJSUkoKytDZWVlpy9Tw01kxWdvZYFZo70AANtOszvtYEYprlY2wtHaAveM9RG7HCIik9DvFqMNGzZg8+bNWLx4sS7rIerRfeOGYEdKAf53tgDP3RMKuYX5bizbMej6/qghsJHLRK6GiMg09Pu3ilKpxKRJk3RZC1GvJoe4wsPBClUNLdiXXiJ2OaIprG5E0sViAOxGIyLSpX4Ho8ceewxffvmlLmsh6pWFTIr57Wv1bDfj7rQtx/OhFoDYIBcM9XAQuxwiIpPR7660pqYmfPTRR9i7dy/Gjh0LS0vLTs+vXbt2wMUR3cp9kb746GA2ki4Vo6pBCSdbudglDSqVWsDWE/kArm+XQkREutHvYHT27FlEREQAAM6fP9/pOYmEWxKQ/ozydsQob0dcLKzBzrOFeMjMwkFydjmKapqgsLHEzNGeYpdDRGRS+h2M9u3bp8s6iPokPtIXrxTWYNvpq2YXjLa3L3A5J8wbVhYcdE1EpEsDmtJz6NAhPPTQQ5g0aRKuXWv7Yf3555/j8OHDOimOqDvzInwglQCn86pwtdJ8tghpalFh1/kiAOC+aEREetDvYPTdd99h5syZsLGxwenTp9Hc3AwAqK6uxquvvqqzAoluxcPRGuMDXQAAuy8Ui1zN4Em6WILa5lb4OtlgfIB5L3BJRKQP/Q5G//jHP7BhwwZ8/PHHnQZeT548GadPn9ZJcUQ96VjscdeFIpErGTw7UtpaZu+N8IFUyrF8RES61u9glJ6ejqlTp3Z5XKFQoKqqaiA1EWll5pi2YHTiSgVKa5tFrkb/qhqU2N++dhO70YiI9KPfwcjLywuZmZldHj98+DCCg4MHVBSRNnydbDB2iAKCAOy9aPrdad+fK0SLSsAob0cM9+TaRURE+tDvYLRs2TL88Y9/RHJyMiQSCQoKCvDFF19g1apVePzxx3VZI1G3ZrZ3p/1kBt1pO9pno82P4L5oRET60u/p+qtXr4Zarcb06dPR0NCAqVOnwsrKCqtWrcKTTz6pyxqJujVztBfe+Ckdv2SWoaapBY7Wlr2/yQjlVzTgxJVKSCRt44uIiEg/+t1iJJFI8Oyzz6KiogLnz5/HsWPHUFpaipdfflmX9RH1aKiHPYZ62KNFJWDfJdPdO+2/qQUAgAlBrvBW2IhcDRGR6ep3MMrLy4MgCJDL5QgNDUVMTAzs7e01zxENFs3stPOm2Z0mCIJmUUcOuiYi0q9+B6OgoCCUlpZ2eby8vBxBQUEDKoqoL2a1z07bn16KRqVK5Gp070JBDTJL6iC3kGJWmJfY5RARmbR+ByNBEG65J1pdXR2sra0HVJQhSkxMRGhoKKKjo8UuhW4y2scRvk42aGxR4WBG17Bu7L4/VwgAmD7Sw2THUBERGYo+D75euXIlgLYxRs899xxsbW01z6lUKiQnJ2s2lzUlK1aswIoVK1BTUwOFQiF2OXQDiUSCWWO88MnhHPx0vkgzU81U7ElrW4qgo2WMiIj0p8/B6MyZMwDaWozOnTsHuVyueU4ulyM8PByrVq3SXYVEWugIRnsvFqNFpYalbEDbABqMnLJ6ZJbUwUIqwR0jPMQuh4jI5PU5GO3btw8AsGTJErz77rtwcOBCcyS+cf7OcLOXo6xOiWPZ5ZgyzF3sknRiT1rbgPIJwa5Q2LAbjYhI3/q9jtGmTZuQlJSEpKQklJSUQK1Wd3p+48aNAy6OSFsyqQR3hXrhq+N52JNWbDLBqGOD3LtCPUWuhIjIPPS7v+Gll17CjBkzkJSUhLKyMlRWVnb6Ihpsd4xoC0OHMspErkQ3yuqacSqv7V6KYzAiIhoU/W4xWr9+PTZv3ozFixfrsh6ifpsU4gqZVIKcsnrkVzTAz8W29zcZsJ8vlkAQrs+6IyIi/et3i5FSqcSkSZN0WQvRgDhYW2KcvxMA02g12p3GbjQiosHW72D02GOP4csvv9RlLUQD1jG26OBl417PqFGpwuHMtmtgMCIiGjz97kpramrCRx99hL1792Ls2LGwtOw8Y2bt2rUDLo6or6YOd8faPZfxS1YZWlVqWBjptP1DGaVoalHD18kGod6OYpdDRGQ2+h2Mzp49q1nI8fz5852eu9WK2ESDIcxXASdbS1Q1tCD1ajWiApzFLqlf9tzQjcb7iYho8PQ7GHWsZ0RkSGRSCSYPdcP3Zwtx8HKpUQYjlVrAz5dKALAbjYhosBlnPwNRD6YOcwPQ1h1ljE7nVaK8XglHawvEBLmIXQ4RkVnpc4tRfHy8Vq/btm1bn4sh0oWOAdgp+VWobmiBwta4Vozu6EabNtLDZLY2ISIyFn0ORtxAlQydj5MNhnrYI7OkDkeyyjA7zFvskvpkL6fpExGJps/BaNOmTfqog0inpgxzQ2ZJHQ5mGFcwyq9oQHZZPWRSCW4fbhrbmhARGRO205NJmnrDekaCIIhcjfaOZLUtTBk+RAEHa+PqAiQiMgUMRmSSYoNdIJdJca2qETll9WKXo7VfMssBAJOHuolcCRGReWIwIpNkK7fA+MC2qfrGsgq2IAg4ksVgREQkJgYjMllT28foGMu+aZeL61BW1wxrSyki2/d8IyKiwcVgRCZrSvt6Rkezy9GiUotcTe9+yWwLcNGBLrCykIlcDRGReWIwIpM1yssRTraWaFCqcKGgRuxyetUx8JrdaERE4mEwIpMllUoQ5d82zujklQqRq+lZq0qNY9ltNU4OYTAiIhILgxGZtKj2AdincitFrqRnqVerUdfcCoWNJUJ9HMUuh4jIbDEYkUmLDmzba+zElUqDXs/oSPv4oonBrpBJJSJXQ0RkvhiMyKSF+Sogl0lRVteMvIoGscvp1i8d44uGsRuNiEhMfd4SxBQEBgbC0dERUqkUzs7O2Ldvn9glkZ5YW8owxtcRp/OqcPJKJQJc7cQuqYtGpQqnc6sAAJNDXMUthojIzJllMAKAI0eOwN7eXuwyaBCMD3RpC0a5lbg/aojY5XRxMrcCSpUa3gprBLkZXnAjIjIn7Eojkzc+wLBnpnVsAzIpxA0SCccXERGJyeCC0cGDBzF37lz4+PhAIpFgx44dXV6TmJiIwMBAWFtbIzY2FsePH+/TOSQSCW6//XZER0fjiy++0FHlZKii2oNRRkkdqhqUIlfTVcfCjpOHshuNiEhsBteVVl9fj/DwcCxduhTx8fFdnt+6dStWrlyJDRs2IDY2FuvWrcPMmTORnp4ODw8PAEBERARaW1u7vHf37t3w8fHB4cOH4evri8LCQsTFxSEsLAxjx47V+7WROFztrRDsZofssnqczqvEnSM9xS5Jo6pBifMF1QC4sCMRkSEwuGA0e/ZszJ49u9vn165di2XLlmHJkiUAgA0bNuD777/Hxo0bsXr1agBASkpKj+fw9fUFAHh7e2POnDk4ffp0t8GoubkZzc3Nmu9ragx/BWXqKirAGdll9ThxxbCC0bHscggCEOJuB09Ha7HLISIyewbXldYTpVKJU6dOIS4uTvOYVCpFXFwcjh49qtUx6uvrUVtbCwCoq6vDzz//jNGjR3f7+tdeew0KhULz5efnN7CLIFF0rGd06ophLfR4PKetnomcjUZEZBCMKhiVlZVBpVLB07PzX/yenp4oKirS6hjFxcW47bbbEB4ejgkTJuDhhx9GdHR0t69/5plnUF1drfnKz88f0DWQODpWwE69WgVlq+FsKHsmvy0YdYyDIiIicRlcV5q+BQcHIzU1VevXW1lZwcrKSo8V0WAIdrODi50cFfVtY3rG+YsfRJpbr29uG+knfj1ERGRkLUZubm6QyWQoLi7u9HhxcTG8vLxEqoqMgUQi0YQhQ5m2f7GwFspWNZxtLRHgait2OUREBCMLRnK5HFFRUUhKStI8plarkZSUhIkTJ4pYGRmD6MCOYGQY44zO5LXVEeHnxPWLiIgMhMF1pdXV1SEzM1PzfU5ODlJSUuDi4gJ/f3+sXLkSCQkJGD9+PGJiYrBu3TrU19drZqnpS2JiIhITE6FSqfR6HtKf8e3B6FRu24ayYoeRlPwqAECkAXTrERFRG4MLRidPnsS0adM0369cuRIAkJCQgM2bN+OBBx5AaWkp1qxZg6KiIkRERGDXrl1dBmTr2ooVK7BixQrU1NRAoVDo9VykH2N8FZBbSFFer8SV8gbRt984k1cFAIj0dxK1DiIius7ggtEdd9wBQRB6fM0TTzyBJ554YpAqIlNhZSFD+BAFTlypxIkrFaIGo/K6ZuRVNEAiAcL9nESrg4iIOjOqMUZEAzWufVp8RzeWWDrOH+JuD0drS1FrISKi6xiMyKyM8WnrBu2YJi8WTTcaW4uIiAwKg5GWEhMTERoa2uNikGT4Rvs4AgAuFdagVSXeQo8dCzty4DURkWFhMNLSihUrkJaWhhMnTohdCg1AoKsd7OQyNLeqkV1WL0oNKrWA1Py2jWM58JqIyLAwGJFZkUolGOXd1mp0oX1X+8GWVVqHuuZW2MplGO7pIEoNRER0awxGZHY6utPSRBpn1LGw49ghCsikXNiRiMiQMBiR2Rkt8gDs6+sXcXwREZGhYTAisxPq09GVVtPrmln6oFnxmjPSiIgMDoORljgrzXQM87SHhVSC6sYWXKtqHNRz1zW3Ir24FgAQwYHXREQGh8FIS5yVZjqsLGQY1j7oebC7087mV0EQAF8nG3g4WA/quYmIqHcMRmSWRt/QnTaYzmg2jnUa1PMSEZF2GIzILF2fmTa4U/Y58JqIyLAxGJFZEmNmmiAISGlf8TqCA6+JiAwSgxGZpVHebWOMCqubUFGvHJRzltY2o6xOCankeosVEREZFgYjMksO1pYIdLUFMHgrYF8sapuNFuRmB2tL2aCck4iI+obBSEucrm96OrrTBmsF7PSitvOM9GZrERGRoWIw0hKn65ue0EGemXapsK3FaCT3RyMiMlgMRmS2rk/ZH9yuNLYYEREZLgYjMlsdXWnZZfVoULbq9VwtKjWySuoAACO92GJERGSoGIzIbLk7WMHdwQqCAFxs7+bSl5yyeihVathbWcDXyUav5yIiov5jMCKzNlgLPV4sbBvHNMLLAVKpRK/nIiKi/mMwIrM2WFuDpLePLxrBbjQiIoPGYERmbbBWwL7UHoxGMRgRERk0BiMtcR0j09TRYpReVIsWlVpv50nnjDQiIqPAYKQlrmNkmvxdbGErl0GpUiOvokEv56hubMG1qkYAwHCuYUREZNAYjMisSSQSBLvbAYBmOr2udbQW+TrZQGFjqZdzEBGRbjAYkdkb6m4PAMgs1Vcwat8KhOOLiIgMHoMRmb2Q9mCUVVKvl+Nf5Iw0IiKjwWBEZm+oh35bjC4VcvNYIiJjwWBEZi+kPRhll9RBEASdHlutFnC5uC1wcao+EZHhYzAisxfgaguZVILa5laU1Dbr9NjXqhpR19wKuUyKQDc7nR6biIh0j8GIzJ6VhQz+LrYAdD8zrWMrkKEe9rCU8XYjIjJ0/ElNhOsDsHU9zkizsCO70YiIjAKDERGAEA/9rGV0SbPiNYMREZExYDDSErcEMW36ajG6qFnDiDPSiIiMAYORlrgliGnrmLKvy7WMmlpUuFLWdjx2pRERGQcGIyJcbzEqqmlCXXOrTo6ZUVwHtQC42Mnh7mClk2MSEZF+MRgRAVDYWGrCi67GGV28YSsQiUSik2MSEZF+MRgRtQvp2ExWR+OMMtsD1nBPdqMRERkLBiOidpqtQXTUYtTR8tQRuIiIyPAxGBG102wmq6MWo+z2gdfB7cclIiLDx2BE1E6XLUbKVjXyKhoAAMFsMSIiMhoMRkTtOlqMcssb0KJSD+hYeRUNUKkF2Mpl8HK01kV5REQ0CBiMiNp5K6xhK5ehVS0gt7xhQMfKbu+OC3Kz44w0IiIjwmBE1E4ikehsnBHHFxERGScGI6IbdMwgG+g4o44Wo2A3ji8iIjImDEZEN9BsDTLQFqPSjhYjBiMiImPCYKQlbiJrHq53pQ1sz7SOrrQQdqURERkVBiMtcRNZ83B9M9k6CILQr2NUNShRUa8E0Db4moiIjAeDEdENAlztIJNKUNfcipLa5n4do6O1ycvRGnZWFrosj4iI9IzBiOgGcgspAlxsAfR/ALZm4DXHFxERGR0GI6KbBA9wyv71qfoMRkRExobBiOgmIR4Dm7J/fao+B14TERkbBiOim4S0B5rsfs5M41R9IiLjxWBEdJOOFqPsfnSlqW7YToRT9YmIjA+DEdFNOrrACqqb0KBs7dN7r1Y2QKlSQ24hhY+TjT7KIyIiPWIwIrqJs50czraWAICcsr51p3V0owW1T/snIiLjwmBEdAsdM9P6Os4oi1P1iYiMGoMR0S10bP7a12DEqfpERMaNwYjoFjQtRmV9G4DNqfpERMaNwYjoFjpafPrcYsSp+kRERo3BiOgWQtyvT9nXdjPZ2qYWzf5qwZyqT0RklBiMiG7B36VtVlm9UqX1ZrIdM9jc7OVQ2FjqszwiItITBiOiW5BbSOHn3LYOkbZ7pmm60Ti+iIjIaDEYEXWjr1P2szlVn4jI6DEYaSkxMRGhoaGIjo4WuxQaJH2dsp/FqfpEREaPwUhLK1asQFpaGk6cOCF2KTRI+jpln11pRETGj8GIqBt9mbKvVgvIKWNXGhGRsWMwIupGR8C5WtmA5lZVj6/NKa9HU4sa1pZS+LvYDkZ5RESkBwxGRN1wt7eCg5UF1AKQW97Q42vPX6sGAIzydoSFjLcVEZGx4k9wom5IJJIbutN6Hmd0oaAGADDGR6H3uoiISH8YjIh60DEAO6uXcUYdLUZjfB31XhMREekPgxFRD7SZsi8IgiYYjWaLERGRUWMwIuqBNlP2r1Y2oqapFZYyCYZ7OgxWaUREpAcMRkQ9CPG43mLU3WayHa1FI7wcILfgLUVEZMz4U5yoB4GudpBIgOrGFlTUK2/5mvMF7eOL2I1GRGT0GIyIemBtKYOvU9tmstlltx5ndP5a24y00b4MRkRExo7BiKgX1zeT7TrO6MaB12N8OCONiMjYMRgR9aJjZtqtpuwX1zSjvF4JmVSCUd4MRkRExo7BiKgXIT0s8nihfXzRUHd7WFvKBrUuIiLSPQYjol5c70rr2mJ0fXwRW4uIiEwBgxFRLzq2BcmraEBTS+fNZDkjjYjItDAYEfXCy9Eavk42aFUL+Op4XqfnLmi2AmEwIiIyBQxGRL2QSCT4/bQQAMAH+7M0rUbldc0oqG4CAIRyRhoRkUlgMCLSwq+i/ODrZIPS2mb8+1guAOBCQdv4omA3O9hbWYhZHhER6QiDEZEW5BZS/GH6UADAhgNZaFC2asYXcWFHIiLTwWBEpKX4cUPg72KLsjolPj+aiwvtM9K4sCMRkelgMCLSkqVMij9MHwagrdXodF4lAA68JiIyJQxGRH0wP8IHQW52qGxoQWH7wOvRbDEiIjIZDEZEfWAhk+KP7a1GAODrZAMnW7mIFRERkS4xGBH10dxwH802IWO44jURkUkxy2CUk5ODadOmITQ0FGFhYaiv77rVA1F3ZFIJXp43BgGutngg2k/scoiISIckgiAIYhcx2G6//Xb84x//wJQpU1BRUQFHR0dYWGi3Dk1NTQ0UCgWqq6vh6MjWAiIiImOg7e9vs1uV7sKFC7C0tMSUKVMAAC4uLiJXRERERIbC4LrSDh48iLlz58LHxwcSiQQ7duzo8prExEQEBgbC2toasbGxOH78uNbHz8jIgL29PebOnYtx48bh1Vdf1WH1REREZMwMrsWovr4e4eHhWLp0KeLj47s8v3XrVqxcuRIbNmxAbGws1q1bh5kzZyI9PR0eHh4AgIiICLS2tnZ57+7du9Ha2opDhw4hJSUFHh4emDVrFqKjo3HXXXfp/dqIiIjIsBlcMJo9ezZmz57d7fNr167FsmXLsGTJEgDAhg0b8P3332Pjxo1YvXo1ACAlJaXb9/v6+mL8+PHw82sbNDtnzhykpKR0G4yam5vR3Nys+b6mpqavl0RERERGwuC60nqiVCpx6tQpxMXFaR6TSqWIi4vD0aNHtTpGdHQ0SkpKUFlZCbVajYMHD2LUqFHdvv61116DQqHQfHUEKiIiIjI9RhWMysrKoFKp4Onp2elxT09PFBUVaXUMCwsLvPrqq5g6dSrGjh2LYcOG4Z577un29c888wyqq6s1X/n5+QO6BiIiIjJcBteVNhh66667kZWVFaysrPRcERERERkCo2oxcnNzg0wmQ3FxcafHi4uL4eXlJVJVREREZCqMKhjJ5XJERUUhKSlJ85harUZSUhImTpwoYmVERERkCgyuK62urg6ZmZma73NycpCSkgIXFxf4+/tj5cqVSEhIwPjx4xETE4N169ahvr5eM0tNXxITE5GYmAiVSqXX8xAREZF4DG5LkP3792PatGldHk9ISMDmzZsBAO+//z7eeOMNFBUVISIiAu+++y5iY2MHpT5uCUJERGR8tP39bXDByNAxGBERERkfbX9/G9UYIyIiIiJ9YjAiIiIiamdwg68NVcfg64492Lg1CBERkfHo+L3d2wgijjHqo6tXr3JbECIiIiOVn5+PIUOGdPs8g1EfqdVqFBQUwMHBARKJRGfHrampgZ+fH/Lz8012ULepXyOvz/iZ+jWa+vUBpn+NvL7+EwQBtbW18PHxgVTa/UgidqX1kVQq7TFpDpSjo6NJ/mO/kalfI6/P+Jn6NZr69QGmf428vv5RKBS9voaDr4mIiIjaMRgRERERtWMwMhBWVlZ4/vnnYWVlJXYpemPq18jrM36mfo2mfn2A6V8jr0//OPiaiIiIqB1bjIiIiIjaMRgRERERtWMwIiIiImrHYERERETUjsFID1QqFZ577jkEBQXBxsYGISEhePnllzvtzyIIAtasWQNvb2/Y2NggLi4OGRkZvR47MTERgYGBsLa2RmxsLI4fP67PS7ml3q6vpaUFf/3rXxEWFgY7Ozv4+Pjg4YcfRkFBQY/HfeGFFyCRSDp9jRw5cjAuqQttPsNHHnmkS72zZs3q9djG8BkC6HJtHV9vvPFGt8c1pM+wtrYWf/rTnxAQEAAbGxtMmjQJJ06c0DxvzPdgh56u0RTuw94+Q2O+B4Her8/Y7sGDBw9i7ty58PHxgUQiwY4dOzo9r809V1FRgUWLFsHR0RFOTk549NFHUVdX1+N5m5qasGLFCri6usLe3h73338/iouL+38hAuncK6+8Iri6ugo7d+4UcnJyhG+++Uawt7cX3nnnHc1r/vnPfwoKhULYsWOHkJqaKtx7771CUFCQ0NjY2O1xt2zZIsjlcmHjxo3ChQsXhGXLlglOTk5CcXHxYFyWRm/XV1VVJcTFxQlbt24VLl26JBw9elSIiYkRoqKiejzu888/L4wePVooLCzUfJWWlg7GJXWhzWeYkJAgzJo1q1O9FRUVPR7XWD5DQRA6XVdhYaGwceNGQSKRCFlZWd0e15A+w1//+tdCaGiocODAASEjI0N4/vnnBUdHR+Hq1auCIBj3Pdihp2s0hfuwt8/QmO9BQej9+oztHvzhhx+EZ599Vti2bZsAQNi+fXun57W552bNmiWEh4cLx44dEw4dOiQMHTpUWLhwYY/n/d3vfif4+fkJSUlJwsmTJ4UJEyYIkyZN6vd1MBjpwd133y0sXbq002Px8fHCokWLBEEQBLVaLXh5eQlvvPGG5vmqqirByspK+Oqrr7o9bkxMjLBixQrN9yqVSvDx8RFee+01HV9Bz3q7vls5fvy4AEDIzc3t9jXPP/+8EB4erqsyB0Sba0xISBDmzZvXp+Ma82c4b9484c477+zxuIbyGTY0NAgymUzYuXNnp8fHjRsnPPvss0Z/DwpC79d4K8Z0H2pzfcZ8D/bn8zOme/DmYKTNPZeWliYAEE6cOKF5zY8//ihIJBLh2rVrtzxPVVWVYGlpKXzzzTeaxy5evCgAEI4ePdqv2tmVpgeTJk1CUlISLl++DABITU3F4cOHMXv2bABATk4OioqKEBcXp3mPQqFAbGwsjh49estjKpVKnDp1qtN7pFIp4uLiun2PvvR2fbdSXV0NiUQCJyenHo+dkZEBHx8fBAcHY9GiRcjLy9Nl6VrT9hr3798PDw8PjBgxAo8//jjKy8u7PaYxf4bFxcX4/vvv8eijj/Z6bEP4DFtbW6FSqWBtbd3pcRsbGxw+fNjo70Gg92u8FWO6D7W9PmO9B/v6+RnbPXgzbe65o0ePwsnJCePHj9e8Ji4uDlKpFMnJybc87qlTp9DS0tLpuCNHjoS/v3+/P1NuIqsHq1evRk1NDUaOHAmZTAaVSoVXXnkFixYtAgAUFRUBADw9PTu9z9PTU/PczcrKyqBSqW75nkuXLunhKrrX2/XdrKmpCX/961+xcOHCHjcFjI2NxebNmzFixAgUFhbixRdfxJQpU3D+/Hk4ODjo63JuSZtrnDVrFuLj4xEUFISsrCz87W9/w+zZs3H06FHIZLIuxzTmz/DTTz+Fg4MD4uPjezyuoXyGDg4OmDhxIl5++WWMGjUKnp6e+Oqrr3D06FEMHTrU6O9BoPdrvJmx3YfaXJ8x34N9/fyM7R68mTb3XFFRETw8PDo9b2FhARcXl27vy6KiIsjl8i5hv6d7uTcMRnrw9ddf44svvsCXX36J0aNHIyUlBX/605/g4+ODhIQEscsbsL5cX0tLC379619DEASsX7++x+Pe2FoxduxYxMbGIiAgAF9//bVWfyXpkjbX+OCDD2peHxYWhrFjxyIkJAT79+/H9OnTB7Xevurrv9GNGzdi0aJFXf66vZkhfYaff/45li5dCl9fX8hkMowbNw4LFy7EqVOnBrUOfdL2Go31Puzt+oz5HgT69m/UGO9BY8WuND14+umnsXr1ajz44IMICwvD4sWL8ec//xmvvfYaAMDLywsAuoyaLy4u1jx3Mzc3N8hksj69R196u74OHT+Mc3NzsWfPnh7/Sr0VJycnDB8+HJmZmbosXyvaXuONgoOD4ebm1m29xvgZAsChQ4eQnp6Oxx57rM/nEfMzDAkJwYEDB1BXV4f8/HwcP34cLS0tCA4ONvp7sENP19jBmO9Dba7vRsZ0DwLaX5+x3oM30uae8/LyQklJSafnW1tbUVFR0e3n4+XlBaVSiaqqqm6P21cMRnrQ0NAAqbTz/1qZTAa1Wg0ACAoKgpeXF5KSkjTP19TUIDk5GRMnTrzlMeVyOaKiojq9R61WIykpqdv36Etv1wdc/2GckZGBvXv3wtXVtc/nqaurQ1ZWFry9vQdcc19pc403u3r1KsrLy7ut19g+ww6ffPIJoqKiEB4e3ufziPkZdrCzs4O3tzcqKyvx008/Yd68eUZ/D97sVtcIGP992KG767uZMd2DN+rt+oz9HgS0+703ceJEVFVVdWox+/nnn6FWqxEbG3vL40ZFRcHS0rLTcdPT05GXl9f/z7RfQ7apRwkJCYKvr69mKvS2bdsENzc34S9/+YvmNf/85z8FJycn4T//+Y9w9uxZYd68eV2mLd55553Ce++9p/l+y5YtgpWVlbB582YhLS1NWL58ueDk5CQUFRUZ1PUplUrh3nvvFYYMGSKkpKR0mjba3Nzc7fU99dRTwv79+4WcnBzhl19+EeLi4gQ3NzehpKRkUK9Pm2usra0VVq1aJRw9elTIyckR9u7dK4wbN04YNmyY0NTU1O01Gstn2KG6ulqwtbUV1q9ff8vjGPJnuGvXLuHHH38UsrOzhd27dwvh4eFCbGysoFQqBUEw7nuwQ0/XaAr3YU/XZ+z3YG/X18GY7sHa2lrhzJkzwpkzZwQAwtq1a4UzZ85oZkFqc8/NmjVLiIyMFJKTk4XDhw8Lw4YN6zRd/+rVq8KIESOE5ORkzWO/+93vBH9/f+Hnn38WTp48KUycOFGYOHFiv6+DwUgPampqhD/+8Y+Cv7+/YG1tLQQHBwvPPvtspx9GarVaeO655wRPT0/ByspKmD59upCent7pOAEBAcLzzz/f6bH33ntP8Pf3F+RyuRATEyMcO3ZsMC6pk96uLycnRwBwy699+/ZpjnPz9T3wwAOCt7e3IJfLBV9fX+GBBx4QMjMzB/nq2vR2jQ0NDcKMGTMEd3d3wdLSUggICBCWLVvW5YersX6GHT788EPBxsZGqKqquuVxDPkz3Lp1qxAcHCzI5XLBy8tLWLFiRafrMOZ7sENP12gK92FP12fs96Ag9P5vVBCM6x7ct2/fLf+9JSQkCIKg3T1XXl4uLFy4ULC3txccHR2FJUuWCLW1tZrnO/5d3/hvuLGxUfj9738vODs7C7a2tsJ9990nFBYW9vs6JIJww1K3RERERGaMY4yIiIiI2jEYEREREbVjMCIiIiJqx2BERERE1I7BiIiIiKgdgxERERFROwYjIiIionYMRkRERETtGIyIyGjs378fEomky4aRgyUpKQmjRo2CSqXSyfFWr16NJ598UifHIiLd4MrXRGSQ7rjjDkRERGDdunWax5RKJSoqKuDp6QmJRDLoNUVFRWHlypVYtGiRTo5XVlaG4OBgpKSkdLtjPBENLrYYEZHRkMvl8PLyEiUUHT58GFlZWbj//vt1dkw3NzfMnDkT69ev19kxiWhgGIyIyOA88sgjOHDgAN555x1IJBJIJBJcuXKlS1fa5s2b4eTkhJ07d2LEiBGwtbXFggUL0NDQgE8//RSBgYFwdnbGH/7wh07dX83NzVi1ahV8fX1hZ2eH2NhY7N+/v8eatmzZgrvuugvW1taax1544QVERETg888/R2BgIBQKBR588EHU1tZqXvPtt98iLCwMNjY2cHV1RVxcHOrr6zXPz507F1u2bNHN/zgiGjAGIyIyOO+88w4mTpyIZcuWobCwEIWFhfDz87vlaxsaGvDuu+9iy5Yt2LVrF/bv34/77rsPP/zwA3744Qd8/vnn+PDDD/Htt99q3vPEE0/g6NGj2LJlC86ePYtf/epXmDVrFjIyMrqt6dChQxg/fnyXx7OysrBjxw7s3LkTO3fuxIEDB/DPf/4TAFBYWIiFCxdi6dKluHjxIvbv34/4+HjcOIIhJiYGV69exZUrV/r5f4uIdMlC7AKIiG6mUCggl8tha2sLLy+vHl/b0tKC9evXIyQkBACwYMECfP755yguLoa9vT1CQ0Mxbdo07Nu3Dw888ADy8vKwadMm5OXlwcfHBwCwatUq7Nq1C5s2bcKrr756y/Pk5uZqXn8jtVqNzZs3w8HBAQCwePFiJCUl4ZVXXkFhYSFaW1sRHx+PgIAAAEBYWFin93ccMzc3F4GBgdr/TyIivWAwIiKjZmtrqwlFAODp6YnAwEDY29t3eqykpAQAcO7cOahUKgwfPrzTcZqbm+Hq6trteRobGzt1o3UIDAzUhCIA8Pb21pwrPDwc06dPR1hYGGbOnIkZM2ZgwYIFcHZ21rzexsYGQFvLFxGJj8GIiIyapaVlp+8lEsktH1Or1QCAuro6yGQynDp1CjKZrNPrbgxTN3Nzc0NlZaVW5+84l0wmw549e3DkyBHs3r0b7733Hp599lkkJycjKCgIAFBRUQEAcHd31+ZyiUjPOMaIiAySXC7X2XpBN4qMjIRKpUJJSQmGDh3a6aunbrvIyEikpaX1+XwSiQSTJ0/Giy++iDNnzkAul2P79u2a58+fPw9LS0uMHj26X9dDRLrFFiMiMkiBgYFITk7GlStXYG9vDxcXF50cd/jw4Vi0aBEefvhhvPXWW4iMjERpaSmSkpIwduxY3H333bd838yZM/Hpp5/26VzJyclISkrCjBkz4OHhgeTkZJSWlmLUqFGa1xw6dAhTpkzRdKkRkbjYYkREBmnVqlWQyWQIDQ2Fu7s78vLydHbsTZs24eGHH8ZTTz2FESNGYP78+Thx4gT8/f27fc+iRYtw4cIFpKena30eR0dHHDx4EHPmzMHw4cPx97//HW+99RZmz56tec2WLVuwbNmyAV0PEekOV74mItLS008/jZqaGnz44Yc6Od6PP/6Ip556CmfPnoWFBRvwiQwBW4yIiLT07LPPIiAgQDO4eqDq6+uxadMmhiIiA8IWIyIiIqJ2bDEiIiIiasdgRERERNSOwYiIiIioHYMRERERUTsGIyIiIqJ2DEZERERE7RiMiIiIiNoxGBERERG1YzAiIiIiavf/iBhbyiPRGfgAAAAASUVORK5CYII=", "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": [ "33358.68514589628\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.11.9" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }