{ "cells": [ { "cell_type": "markdown", "id": "3e6f4cef-d07e-4b44-85af-5a17674e271d", "metadata": {}, "source": [ "# FindPeaks" ] }, { "cell_type": "code", "execution_count": 1, "id": "0de772f9-4d9e-4ade-a04e-b20f1f04207e", "metadata": {}, "outputs": [], "source": [ "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "ec8317f2-ce43-419c-a3d4-058d1478ca36", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[ 7.395 13.572 19.923 24.534 30.711 36.975]\n", "[ 656.10337734 968.54481784 1472.52126398 1429.91065595 1031.57776675\n", " 642.49844298]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAACpMUlEQVR4nO2deVwU9f/HX8sqiAool6iL4pVHHqmZXywKjzy+VuaKpJlZmWbhgZqWv0zUDktNscOsvqWVmqKsWXaJBoaplRZZZqYmKiieCaKJsMzvD5ppZndmd2Z39pjZ9/Px4KE7Mzv7meszr8/7+hgYhmFAEARBEAShIYJ83QCCIAiCIAilkIAhCIIgCEJzkIAhCIIgCEJzkIAhCIIgCEJzkIAhCIIgCEJzkIAhCIIgCEJzkIAhCIIgCEJzkIAhCIIgCEJz1PJ1AzxFdXU1Tp06hbCwMBgMBl83hyAIgiAIGTAMg8uXL6NJkyYICpK2s+hWwJw6dQrx8fG+bgZBEARBEC5w8uRJmEwmyfW6FTBhYWEAak5AeHi4j1sToFy/DrzySs3/p08HgoN92x5vEajHzUfL50DLbXeHQD1ulkA/fj+irKwM8fHx3HtcCt0KGNZtFB4eTgLGV1y/DoSE1Pw/PDxwOoRAPW4+Wj4HWm67OwTqcbME+vH7Ic7CPyiIlyAIgiAIzUEChiAIgiAIzUEChiAIgiAIzUEChiAIgiAIzUEChiAIgiAIzUEChiAIgiAIzUEChiAIgiAIzUEChiAIgiAIzaHbQnZEYGO1WpGfn4/Tp0+jcePGSEpKgtFo9HWzCIIgCJUgAUPoDovFgilTpqCoqIhbZjKZsGzZMpjNZh+2jCCUQUKcIKQhFxKhKywWC1JSUgTiBQCKi4uRkpICi8Xio5YRhDIsFgsSEhLQu3dv3H///ejduzcSEhLoHiaIfyABQ+gGq9WKKVOmgGEYu3XssvT0dFitVm83jSAUQUKcIJxDAobQDfn5+XYdPh+GYXDy5Enk5+d7sVUEoQwS4gQhDxIwhG44ffq0qtsRhC8gIU4Q8iABQ+iGxo0bq7odQfgCEuIEIQ/KQiJ0Q1JSEkwmE4qLi0XN7waDASaTCUlJST5onXMo44QASIgThFzIAkPoBqPRiGXLlgGoESt82M+ZmZl+KQoo44RgYYW47T3MYjAYEB8f77dCnCC8BQkYQleYzWZs3LgRTZs2FSw3mUzYuHGjX9aBoYwTgo+WhThBeBMSMITuMJvNKCws5D43atQIx44d80vxQhknhBisELd1E/mzECcIb0MChtAl/NFpSEiI345WKeOEkMJsNuOnn37iPi9cuNBvhTjhX1itVuTl5eGjjz5CXl6ebgdAFMRL+D3uBrdWV1d7sHXuQRknhCP4927Hjh39VogT/kMgTaVCAobwa9R4GP1ZwFDGCeGIyspK7v/Xr1/3YUscQxl03kXqfLPxdLYuaTaeTm/uR0UupAULFqBHjx4ICwtDbGws7r33Xhw6dEiwTXJyMgwGg+BvwoQJgm1OnDiBwYMHo27duoiNjcWMGTNQVVUl2CYvLw/dunVDSEgIWrdujVWrVrl2hIRmUSu41Z8FDGWcEI7QgoDxdgZdoLhHxLBarZg/fz5iY2MF5zs2NhZz584NuHg6RQJmx44dSEtLw549e5CTk4PKykr0798fV65cEWw3btw4nD59mvtbuHAht85qtWLw4MG4fv06du3ahffffx+rVq3CnDlzuG2OHTuGwYMHo3fv3igoKEB6ejoeffRRfPXVV24eLqEV1Axu9ecHlp9xYgtlnBD+LmC8nUEXyOUGLBYLGjVqhIyMDFy8eFGw7uLFi5g3b17AxdMpciF9+eWXgs+rVq1CbGws9u3bh9tvv51bXrduXcTFxYnuY+vWrfjtt9+wbds2NGrUCDfddBOee+45PPXUU5g7dy6Cg4OxYsUKtGjRAq+88goAoH379ti5cyeWLl2KAQMGKD1Gv4RMro5REtyanJzscF/+bIEB/s04sTX9mkwmZGZmumTy1cr9pZV2+gp/FjDOBhkGgwHp6ekYMmSIKtc00NwjfCwWC4YNG6bKvvQUT+dWFlJpaSkAIDIyUrB8zZo1iI6ORseOHTFr1ixcvXqVW7d792506tQJjRo14pYNGDAAZWVlOHDgALdNv379BPscMGAAdu/eLdmWiooKlJWVCf78lUAeRchFzeBWsQ7W3zCbzahXrx73OTc31+WME63cX1pppy/xZwHjzQy6QC43wB67Wugpns5lAVNdXY309HTceuut6NixI7f8/vvvx+rVq5Gbm4tZs2bhww8/xAMPPMCtLykpEYgXANznkpISh9uUlZXh77//Fm3PggULEBERwf3Fx8e7emgehYqWyUPN4FYtdmrJyckujVq1cn9ppZ2+xp8FjDcz6AK53ICzY5eLHuPpXM5CSktLw6+//oqdO3cKlo8fP577f6dOndC4cWP07dsXR48eRatWrVxvqRNmzZqFadOmcZ/Lysr8TsR42+SqZc6dOwej0SgpPpzNa8T/nr+7kFjcbadW7i+ttNMf4Cc3+JuAcXWQ4YrbMJDLDahxTO7E0/mzm9clC8zEiROxZcsW5ObmwmQyOdy2Z8+eAIAjR44AAOLi4nDmzBnBNuxnNm5Gapvw8HCEhoaK/k5ISAjCw8MFf/5GII8ilGCxWHDfffc5tZw4ehj5I1etCBh30cr9pZV2qomrmTP+bIFxJYPOVbdhIJcbUOOYXK3g7O9uXkUChmEYTJw4EZs2bcLXX3+NFi1aOP1OQUEBgH8vQmJiIn755RecPXuW2yYnJwfh4eHo0KEDt8327dsF+8nJyUFiYqKS5voNbOeVnZ0ta3tPjCK0knroaHTOYjQakZWV5fBhrKioEOxTC7gbq6OVUapW2qkW7rwE/FnAKM2gc8dtGMjlBpwduy2227kST8emaw8bNsyv3byKBExaWhpWr16NtWvXIiwsDCUlJSgpKeHiUo4ePYrnnnsO+/btQ2FhIT755BM8+OCDuP3229G5c2cAQP/+/dGhQweMHj0aP//8M7766ivMnj0baWlpCAkJAQBMmDABf/75J2bOnInff/8dy5cvR1ZWFqZOnary4Xsefuf1+uuvy/qO2qMIf1fRfOT4e61WK6Kjox1uw+/sA0XAaGWUqpV2OkLugMDdWB++gOH/319gM+hsLaG2I353g3DlTnAJQBMDNSU4Eoq2iIkcpfF0FosFzZs3R0ZGhuh6vwqaZhQAQPRv5cqVDMMwzIkTJ5jbb7+diYyMZEJCQpjWrVszM2bMYEpLSwX7KSwsZAYNGsSEhoYy0dHRzPTp05nKykrBNrm5ucxNN93EBAcHMy1btuR+Qy6lpaUMALvf9ibZ2dmMwWCQPG+2fwaDgYmPj2eqqqo83gaDwcAYDAYmOztbtd+yo6KCYTIyav4qKmR9Ze3atbLO1dq1ax3up6ioSLB9dXW1+8cjFxeOm2EYJiQkhGuvK1RVVTEmk0nynvPE/SWJg3PgV+0Uw8n1y87OZkwmk6DNJpPJ7llij9Od533z5s3c9tOnT2dyc3OZtWvXMrm5ueqfHxfvW4ZhmLi4OK6dYm3Lzc2V9Vzn5uY6/J3s7GymSZMmgu/Ex8cz2dnZsq+LJ47fG2RnZzN16tRxeP7Yc8FfpvQ3lLyznF0vV5H7/natp9QAvhYwzjovbwgKVzrQqqoq9TpJFzoEtTq6P//8U7D91atXXT8OpbjYEQYHB7slYBjm3w7IthPyimDlI0ME+EU7xXDQdiUDAjXu5Y0bN3Lb1a9f3/WXs5vH7Yzo6GiH965aAxOGEQ5O5s6dy1RVVakzUPNzAcMwDHPHHXdInrvBgwczVVVVTHV1tWB5hZNjYfv81atXMzExMbLfWXKvlyvIfX/TbNQeQmnqm6tBVu60gbEJlvQHV5Mcf68cX7dtvMC1a9dUaZ+/wro0KioqMHfuXLsyBJ64v9yBdT00bdpUsNzf2slHqRtEjVgfvtuovLxcsM6fYhFsp4KxxVNuQzazNVBqxFy+fFlyXVxcHIxGo12pEdv7hg+/z3/ggQdw7tw5Re3xtZuXJnNUAbE0MyUBiOvWrUNKSorqqWlKOlB/qXLJ+ntTUlJgMBhEOyU5qYC2Aubvv/9Gw4YNVW2r2ogdqxzEJrzkC5hx48bhzTff9JvURxaz2YwhQ4agVq1/u6Fjx475XTtZlFaHViPN+Oeff3b4e/6Scu5MwLADk+LiYtH73FlZBKnfunr1qqpVu/0dRwVaq6qqkJeXh4MHDwqWl5eX2xWbBaQrG8tByfXyJCRg3ERqtuRx48bJ3kf37t090vnI7UBjY2Px0EMP+U1dDnZ0bnteAeCpp56SJaT4WUiANiwwrnQkUp0QP8svIiLCb0WBbbv8tZ2AcouKKy9tsf7EEf7ychYTMKxVMC8vDwAwduxYzJ8/325gorRGia2ACZSsNqvVivPnz0uu37BhA95//3275WzFfNt9jR8/3q3EAX+Yo41cSG7gKMMgIyMDUVFRslLfPJVdICf10GQyoaCgwO/qcpjNZhQWFtot79Gjh6zv21pgvvnmG783IcvpTPjZL9u3b3dqOgeA77//XjcZGb5EqUVFbuaMszRjOfj65WwrYNiJB/v164fnn38ezz//PObNm4d69erZWUKVug1tBYwestpYpLLbWFfPpUuXJL/Ln7KHT2JiIjZs2CBY9sILL+DChQsutdGf3LwkYFxETjVRKWzX8QWMmvVanHWgDMPg77//xpNPPilrf97uJMXUPf98ODpX27ZtE3zv4Ycf9tvUcRZnAsY2Rqlfv36yXnbffPONX6fOawVXapGw1sSYmBjBtkrSjOXg65czX1SwEw+KvSDLy8sFMym//PLLimuU8H+rvLwcVqtV1EXCopUaMVIxiDNnznRZ2ALAlStXkJqaipkzZwKoudfkpmXbMn78eBQWFvqFeAHgRrqDn+PpLCS5GQbz5s2zS32Lj49nZs+ezX3eu3cvwzDy0zOVIrbfqKgoRdHmkJH5Y4cKUf22bVizZo3kMbHnyjaNkP3zWoaLi8fNz6KwRWl6o0+Pn2EUnQN+G5WiatYci0pZSPz2zZkzh9t29erVLqcZi/1FRUV5/LidwW9P06ZNZbf9888/V9zMgoIC7vu2mVlu3fM+zEJS4/mW87dhwwa37rX8/HyvnA/KQvIwcq0Rbdq0weDBg7nPn3/+OY4dO8ZNsQDUWGA8ObmdrTsmJSVFckoGMfxpBHPt2jWHFSKHDRsmmI+LD+PnGQkMb/TN/7+7o3Pbfap1/L6s7sxWCo2NjfVq1hxrUbFFzKzOH1HPnz+fW37+/Hk766I71s0LFy5g8+bNLn9fbYqLi2Vv6yz419l3HGXYAP7l7pBCredbDk888YSi6wMAderU4f5vG1voayiI10WU+F2Dgv7ViZ07d4bRaBS4ja5du+axye2uX7+O5cuX4+jRo9wyo9GoyBzJMAweffRRRb+rBhaLxS7gb+zYsZLbs9s58u0yEkGP/jZhmdVq5bJz1JqNFpA+fqVIBa8vW7bMrZfFRx99JHr++dfn8OHDePXVV0Wvszey5mz3m5uba9deRxke6enpiI+PF+zHXReQP2QiuYIr8X9yX6Lz5s3DM8884/fnRM3n2xnnzp1TnCrNv4f9LRmCLDAuosQfzg+uYiPC+Q/uvn37PBJEO3PmTNStWxdTp04VTGOwd+9eRfsBgIyMDK/GULiT4icH/ojXH+rf2MK/PzwRe+TOPj1pLRQ7/7bXJyMjQ1Kk+sLKZluqXc6I2rZ9Sue7sUWrk1+6YoHZt2+frO2CgoL8XrwA3o8tjImJcXqvRUVFcf/n36f+ZoEhAeMiSiYy4wsYNoqc/4LauXOnrN+UutGvX7+OzMxMTJo0CZmZmbh+/TpmzpyJRYsWiXbifGuMErxVOMsbJlV2xLtx40a/nLCMf394IkDT1X26O6eNLVLnlz3/rgQwuir41ULOiPrkyZOYO3cu53pzFHAvF39xI9kWJ3SEKxYYuRYERynH/oS3A7CbNm3q9F5bsWIF93++yCQBoyNYf3j9+vUFy239rleuXOHWbd68mauYyvL111/L+r3Dhw/DarVi+/btePbZZ/Hss8/ivvvuE1hZpk6dipCQECxevFiFIxTirdGtJ02qfMvYhg0bMGLECNHtvD2StxUE/I69V69edlksfNh0eDm4G88kt2iYHFHOiiGp/TAMgyVLlrgsZH2VWiz3d59//nmBtUmqOrFc1qxZ4xexXa+++qrsbV2xwERERMjaLiwsTPG+fYG71jclGI1GnD9/XvJei42NRXZ2Nv773/+Kfp8EjM4wm81IS0vjPq9fv94uLZDfoS1cuBC9e/cWpC47qq7IJyMjAw0aNBDUVsjKyhLttDxlvfDG6NZTLx6+ZWzz5s1ITU112OF7cyRv25GzAsZisaBVq1aSo072mOSkRSotGMZiW3tGDiUlJU63kTvzuKucOXPGJy90pSNqvrWPDbiPj49X/Lvnzp3zidXJ9t41m83YsGEDgoOD7baNiorCnDlzJL8rh/bt28varmXLlor37WnEAt/lzDZ91113ITw8XJXfT01N5e61I0eOCNZ/8MEHMJvNdtMRsFAMjA7hj5bbtWtnF8x38uRJu+/89ddfLv2Ws6h7b+HJ0a2nTKqsZWzIkCGSI38xvDGSty285ygzjY9UlsUjjzxit210dDSmTJmCyMhIRS4efvzJ888/L+t7cXFxTrfx9HmdOnWqT2KZ2BG1XGytfUajEdXV1S79ti+sTrYvtY0bN2Lq1Kl29/STTz6JM2fO4Oabb+aWueJC+vbbbyXX8a0Y/mCNAv4VLVOnTkXjxo1FY+2krPnsu2TLli2yB7pyYO8123PEfpYqiudvFhjKQlIB/gPMFyaOTORax5N+W2cl2F3llVdegdlsRl5eniIXlTd81GKTTzqLAwoODsaRI0dER7rt2rUTfI6KisK5c+eQmZmJzMxMWRlDrgRSs+6s2267DXDwogG8c169PZcXUOMmlhrBSsG39l28eBGnTp1y6bd9UdDOVsCkpqaK3jOLFy/GsWPHBMscWWBspyJITk7GX3/9hRdffFHyOxEREVycoe0z5Qs2btyIJ554QtKCant/5uTkCOJPPCHC+Pda165dBevYcyZ1/5KA0Qn8tE7+Q/nJJ59wD6XcF2X9+vVx5coVr9QBcBdvTOLFn9BRTaZPnw6z2axolOqt+je2I9Hdu3c7vXeuX7+OXbt2iaZD23betlk7/I5zyJAhdinkgPQMv1IodVGxQtWTKaTuliFQirvZc5s3b8ayZcsUf9+Xk+vZvuwctT07O1vwWSqjyGKxYPz48YL79vnnnxeUpBCjdu3a3P99LWDYRApH2N6falpZnLF9+3a7gQ4rULTiQiIB4wKOJlxbsmQJlixZomh/N9xwA3788Ue1mucxXI2hcAWz2Yw1a9bg/vvvV22f7KhDySjVWxOW2Xa2cmJIAGmXgbOREttxjh8/3u5ejo6ORlpammJhYTKZkJmZWWPpkPnyGDduHDIyMhT9jlLYEedrr72GSZMmeex6qpE9995777kkXgDfTa7HT1JQyqpVq3D33XfbFQAcNmyY6PbOXGt8S4en5pizRayGlMVicSpeWGytb97i+eefx7vvvitYxvZDUi6kv//+G9u3bxdYxWzLCHgTioFRiDsTrknhr+LFNire21Ut+RWMe/fuLVhXr149O/OnHE6fPi0r6j8oKAhz585FRUWFV6rM2goY2wnvpGDFmO1LT87ok2EYXLhwwe5ePn/+PObNmyfr91lyc3MVzWnDxtZIiRdPdIiejolRI3vOlRG4yWTC+vXrERkZ6fWqyBaLxe1ZsKdMmcK112q1YvLkySq0zDsWGLEaUs2bN3ep8GdxcbHoBLaexHagtH//fnz00UfYvXu36PaZmZmCJJJ+/fqhUaNGPis3QRYYBXiz5LM/ULduXW50lZubi169emHXrl2S1VLVhj+C6tixI3Jzc7nPV69exU8//aR4n40bNxa4qGwr/bI0bNgQc+fO5T6rUWXWEbadbadOnZzGAdWtW5dzGfC/bzAYvG4+V/ISc+ZmueGGG/DHH3+o1DIhnoyJ8XYArcFgQE5ODv766y9MnTpV9arIzlCr2GRRURFXGTo/P19xqXspPP0MSB2/q+2fOnWq4iq57mLbdmflN8RcSBcuXMCwYcOQnZ3t9SkbyAKjAG+WfPYkTZo0kbUd399cUlKCVq1aebVaLStgDAYDQkJCBOtc6TT58Sxs1H90dLTotlIxI546XtvOtrq62mlq5c0338wJSL7JNygoyOPBdgaDwaU6G3IGAcePH3enaQ7xZH2f2NhYVffnDIZh8O233yI1NdWuXyoqKsKwYcMwdepUj1hk1B7MseLPHRFoMBgEz7MnXUieGMyqLV5uvfVWVffnDL4lzVuQgFGArwpjqY3cCpXsXDxATYl3T1WrlZoUkO2AateurYqlxzZGwGw246233pL1XU8XtrPtbCsrKzmRxZ9MjU/z5s25//MFjNVq9aiAYV1vffv2VfxdOYMAT4svT9T3sVgsGDNmjGr7k4uzgN/MzEyPDDbUHsyxrlB3sqgYhsGDDz7IffakBUYLg9lOnTp59fdYS5o3IQGjAF+kKHqC69evo27duk6347+o1Sodb4ujeYj4AoYvplyhSZMmnHmTL5i+++472fvwZGE7sTowQI3ISkxM5JZ//PHH3P/5MTy2WQOXL19WrW2RkZGCz2ws1A033KB4X3IHAZGRkYork3bu3NnOUqdGW5zBuhLcdX0YDAbBHDRykBv0qbYFUc3BHFuXyGq1IikpyeVKxHPmzBHExXnSAqPW8TuqsO0uDRo08Ni+pfD2IJ8EjAK8WfLZ03Ts2NHpNnJqWbjzUnc2KeBnn30GwDUBw9ZGGTBgAIB/3WG2gunll19W3G5PPKS255rva+ZbJHr06MH9n5+RYZs1wE4a6gi593FWVhb3//DwcC5Y1xXzudxBgCv1k/744w9F1pvDhw8r/g0x1HAlsNfi7bffFlSqVQu1LYhyr6NtYTYxLl68iH79+iEhIQGbN29GZmamS2266aabBAMBT1pg1BjMvvDCC1i6dKkKrRHHFwLG24N8EjAKUGPCNX9BzuhZyQhG6UtdzqSAbMGqWrVqKRIw48ePx6BBgwD8W3a8tLQUU6dOFZ24USlqP6QWiwX33XefYNlDDz3EjZb54oYvbPgdtK2AEav+bIucl258fLwgQLd+/fqi7jy5lWPlDALq1KmDZ555Bhs3bpS1TxalNSrmzp2rikVCDVdCZGQkF1jcpUsXt/cnhpoWxKSkJFnWonbt2sl2NRYXF2PYsGGYMGGCS226cuWK1wSMGoPZpk2bumxtkoO7Vmul+KIOEQkYhTgL/tQCwcHBqr+ElY5m5UwKyKb4KY2B6dOnDye+2Iqmly9fdnlkx+LuRIhisFYo27ikCxcucCZ/voDhu4b4HbSty0iOBQaomWPFEbZxQ/yqv3wBJHdOGznzvrRs2RJGo9ErGQ2uWCRcmb/HGWzQuMViwbhx4xxu66yYmzPUsCBu3rzZLtBdjPDwcNx4442y9sneT3L2K0Zubi4OHDjAffakC0mNwWx6ejrOnTvnMau+3Gk/1GLZsmVerwdDAkYhVqsVkZGRgmAxrRETE6P6w/3OO+8oehEo6UStVquiTvvq1avcy53v/lCCWIfCMAyGDRuG/Px8VczwzjIZGIbBhAkTBNkJ/Doh7DFu3LjRrmqx3PPLd0nJgV/plI+Slzo7CJDaV7169RS1yVVctUh4YmTPFhVMSUlxGtfi6jxJLO4OXuRMkcJew/r168uKt1OD9957D6+//jr32dNp1O7OHn7p0iXcd999GDlypMot+3f/3iAqKsonKdQACRhF8OMnXnnlFbf2JdZ5r1mzxq19KvlttTM9lEagK+lEz549iwULFsjevry83K2S10aj0S5FmB1ZqJnVISeT4dy5c4IR6eeff879v7KyEjNnzsTw4cPtCqA5e8mxpefffvtth9ukp6cLxIm7FhgWs9mMFi1aiK5zp7KrKyi1SLhadNARbFFBT9aYYq+51Wp1q+CdnPuWvYZhYWE+q5tlO++SJ2BnD3cVhmHw0UcfYd26dYoC0P2FTz/9FGfOnPGJeAGokJ1s1CraxCJmAenWrZsq+3aGwWDwSKqqkheB0gkbxVwi4eHhopVLy8vL3eq8rFar3X5tO3o1CqK5Ysp/6aWXuP8XFRUhJyfHpd9mGAZt2rQRFAcU2+bkyZPYvn07t0wtAQNIp0tfuXIFeXl5Ds+P0WhULZ1dqUXCtt2uzizvbRiGwd9//41+/fpxy1wpeKfkvj137hw2b96sqJ1qsX//fm52b0/i7v6Lioowbtw4v5soUQ7enLtJDLLAyMBbFXi9dQMHBwd7ZFIuJS8CNXzI7MNj63IoLS2VPZeQq6iR1eFu4bNDhw65/N369es7FC98tm7dyv2fHxjIfx5cOQdS9+DJkye5LDEp1BAvrsY0+eqFrAZqFGhU8pxv3boV5eXlsrdXk2vXrnmk5IFt3Sr+YHTAgAGyMq9s8bUQcJVRo0Z5vKCpI0jAyMBbRYukJtBSm8rKStXFkisvAtaHLLcysBS2I6A//vjDKxUh3U0hT01Ndev33YmFUPJS4U9OWlBQwHVWfKuLKxYYKQHjboyHHFydANFqtXp88klv4ooQl5OB4y/uELXFpljdKr4r9IYbbsDQoUO5z/z53PSKp6uUO4IEjAy8VZxn7969quwnPT0dJpNJcv3x48ft/LbuRsGPGDHCJVOq2WzGvn373Ppt25gJV7MYXEXp/cHOtuvNmWfVorKykuus3BUwnrY4RkdHC6oV83F1YtK8vDyPW/e8DSvEd+7cKWt7OZlkrlxbV6amcMaaNWtUG8xI1a1iMx2BmnPDdyk66of1AsMwYBjGY1XKHUECRgbeKs7DdozupkkOGTIEhYWFkjMKW61WO9Hg7jEuXrzYZQXubkaU7UNz9uxZt/anFCXnTk4Gh7/DZkf99ttv3DI1LTBqsWLFCtFpGObPn69o5mwWNaxm/owSYcZaT8PDw1X7/UmTJqlusTh37pwqbiQ5dauAGuHGj79zZcJZreKpKuWOIAEjA29V4JWTbmg0Gh22g+/Keeedd2T/9sGDB/HAAw/I3l4MVxU4vzy+Gvz++++q7k8KV2IotDCHihzOnTuHHTt2cJ/Zqsn+hNVqFRVJ7du3V2wtZEffWrSaySUuLk7R9mazGSNGjFDt91988UWP3EdqWNDlPrcrV64U1KL5/vvv3f5tLeHt+DBFAmbBggXo0aMHwsLCEBsbi3vvvVcQSHjx4kVMmjQJbdu2RWhoKJo1a4bJkyfbZZAYDAa7v3Xr1gm2ycvLQ7du3RASEoLWrVtj1apVrh+lm8gxmarBmTNnADiOAbBarZLBxAaDgfPpK31RhoaGonPnzsoazMPVeBCr1er1gktq4GoMhV4mBLVl8uTJPgvkk6KiokJUwCh1b3griN+XmEwm3HbbbYq/98cff3igNeqihgVd7nPraauiv6Omy04OigTMjh07kJaWhj179iAnJweVlZXo378/F4Nw6tQpnDp1CosXL8avv/6KVatW4csvv8TYsWPt9rVy5UqcPn2a+7v33nu5dceOHcPgwYPRu3dvFBQUID09HY8++ii++uor947WDViTaUREhMd+48MPP3Tr+0uWLOHM4kpflLVq1VLFHKz0d/Pz8znhpiVcjaHQy4SgYvjCB+6I69evqyJg9GI1c8S4ceOcCnHb7Jvr1697pdaKO6hVOVvPz62aqOWyk4uiOjBffvml4POqVasQGxuLffv24fbbb0fHjh2RnZ3NrW/VqhVeeOEFPPDAA6iqqhKkYDZo0EDSZLlixQq0aNGCKxbXvn177Ny5E0uXLuUm5/MVcku0u4LSyokxMTEoKyvjOmR2/h9A+QNnMBhUETBKf1eLFolt27YhOTnZpaBlpfVvbFGz/oma8C1w/LmTpPBGxl1FRYWoWFE6SpZ7jwYHB3u8+qunaNOmjcP1FosFU6ZMEQi5oKAgr2SMuYNSC6kU7j63atKhQwdB/JkUERERHn1fSeHNPt2tGBj25ERGRjrcJjw83G5iqbS0NERHR+OWW27Be++9J7gpdu/eLSi2BNTk1+/evVvydyoqKlBWVib4UxN/DL48d+6coLAYv7qvK3E77ggYV2tqyBU8jmqCeJuqqipkZWW5VMnUXXekuynnUgQHB3vNAmexWJy+MNVAyoV04MABRZVo5d6jLVu2VNxGf8HRMUpl37gjXp588kmXvyuXiIgI1SrEeiuMQA5yK//yB7TexJvWKpcFTHV1NdLT03HrrbeiY8eOotucP38ezz33HMaPHy9YPn/+fGRlZSEnJwfDhg3DE088gddee41bX1JSgkaNGgm+06hRI5SVlQkmtuOzYMECREREcH/x8fGuHpooL7zwgl+akfkZS3wB40qhOHdeYAzDuDTakTOrbZ06ddyeLyQxMdGt7/MZOHAgVwPClSJOrDtSzmy+tsiZZdoV3nzzTVy8eNHtLBtnnRf7MuSnnnoKqVo3y5cvV3T9kpKSZM13405JeSnCwsIwY8YMj6XjOht4eCr+x5WYxpiYGEXbX758WdW4LPa5tR2wq5mJJQe51st169ZhyJAhHm7Nv3hisltnuCxg0tLS8Ouvv9oF37KUlZVh8ODB6NChA+bOnStY9+yzz+LWW29F165d8dRTT2HmzJlYtGiRq00BAMyaNQulpaXcn5odvcVi8dviVfyCUbbzKzmabExsqnV3JtFz5WUM1AitMWPGONwmNjbWYWEsR5Uv2YfKU5YLV4s4mc1mvPzyy4JlUnUwXD23SmCrAtu6iZUQExPjsPPydjCsHPO5nOtnNBrtBmFieCKA8/Lly/jPf/6DwsJCzJ49W9V9ywlE91T8j+3s645gn+GioiKMHj1a9veqq6sxbNgwbNiwwZUmimI2m7Fw4ULu89KlS/Hqq6+qtn+1+fHHH7F+/XpER0d79HdcTWpwF5cEzMSJE7Flyxbk5uaKjgwuX76MgQMHIiwsDJs2bZKcdZalZ8+eKCoq4vzVcXFxdoGdZ86cQXh4OEJDQ0X3ERISgvDwcMGfGvij6wj496Hmp16LnWd2srHc3Fx88MEH3HLbuh0Wi8WtGgwXL150uRqjs+ynEydOYNeuXZLr+/bt6/D7r7zyCrZt26a4XXJwZ0oB2xfs7Nmz8X//93+CZZ988onkPa8mx48fx/PPP++W63XUqFEOOy9vB8PKORa5188bLi8x2Ak1Aef3uVLkBKL7Q4waa90NDg4WWPvlTqI5cuRIbNy4UbX28IVqQkKC4jpW7k4hooSTJ08iNjYWJSUlknXB1MDVpAZ3USRgGIbBxIkTsWnTJnz99deis8mWlZWhf//+CA4OxieffCJaSMqWgoICNGzYkBtlJyYmCiaQA4CcnBxV3QBy8ccMBL7a5buHpISi0WhEcnKyQzdMSkqKWxVG3XmRyzGJvvHGG6LLDQYDOnXqJPm9jRs3IiYmxqPBbEpTyNlsDtvtGzRogJtvvlmwrG7dul65/yZOnOh2B+fMXO3tl6FcMSbn+vkqC4XfNrlxbWIW1zFjxqBPnz7c5w4dOsgq5ucP2Tfp6elcO/lWyiVLlmDp0qVOv2+1WjF8+HDV3En8/uqvv/5SbHnzpoABap47o9GIOXPmIDs7W3V35Lx581wqDKkGigRMWloaVq9ejbVr1yIsLAwlJSUoKSnh4lJY8XLlyhW8++67KCsr47ZhX2qffvop/ve//+HXX3/FkSNH8Oabb+LFF1/EpEmTuN+ZMGEC/vzzT8ycORO///47li9fjqysLEydOlXFQ5eHGp1uTEyMYtP87NmzBRYTPny1y+/M+AG9tlgsFlHByaKGWd/VWjD8wk9SnDt3TnR5cHCwqDsMqLHKmc1mr7045QawsnOp2BbwCw8PF7jxQkJCFFUVZu+FlJQU2d/h4849EBERgeLiYoeBsd5+GSq1Jjm6fqx48BXsS8hZXNt9990nuu7ChQu4fPky97lu3bqyTP3eKuLpCL4w5guY4OBgu1hJR6iV5s+fuuTSpUuicZmOzpfSbFN34T93fIv8I4884va+DQYD/ve//7m9H1dRJGDefPNNlJaWIjk5GY0bN+b+1q9fD6DG3/bdd9/hl19+QevWrQXbsDEptWvXxhtvvIHExETcdNNNeOutt7BkyRJBjEmLFi3w2WefIScnB126dMErr7yC//3vfz5JoVbS6dretGyRvhUrViiu4JmcnCyojMumnL/wwgsCtSvHAiOVReAplAoGJf5wW0JCQiQ7Ytaid/jwYZf3rwS5AaxS12Hjxo0Cl2DdunUV3X8mkwnZ2dlIS0uT/R21KC0txQMPPMAFxm7atMluG2+9DFn3sVIB4+hc+zoLhW2bo7i2tm3bIisrS/T+2rJlC4qLi7nPctO9fX3ctnFV/Hi3WrVqKXo+1Cp1z7fA5Ofn23kLAMeDAbn9sG3QstIgZqmgWtYi7+pAh487E9qqAqNTSktLGQBMaWmpW/upqqpiTCYTYzAYGACifyaTidmwYQNjMpkEy+Pj45ns7GyGYRgmNzdX8vtS+8zOzmaMRiMDgImOjmYAMDk5OYL2tW7dmvuOo/Yr+W13/3Jzc2t+vKKCYTIyav4qKiTP8ZgxY1z+rejoaObFF1+UXJedne10HwaDgYmKinK5DQaDgYmPj2eqqqokj1vudViwYIHgHpDzvZCQEGbbtm3c78u5Zz35ZzAYmGCA+S011e7aZ2dne7xdzZo1YwAwHTt2VO36sTz44INeP5eCtvGe69zcXOaee+7htq1Xr57DfQUFBXH/b9++vf2D6OC4s7OzBd/31t+GDRsE7di2bRu3Ljs7m6mqqmKaNGkie3/p6emS/ZDc/mrAgAEeP+74+HimoqKCyc3NZdauXcvk5uYyq1evVnzvsO8fW7Kzs5mmTZs63UdkZKSs31q7dq30eXUBue9vmgvJCY7MtqyFZdmyZUhJSeFMc2vXrkVubq7AUqJ09MlmR7DbsymhtgGdzvbnzRgeV9PolI4s+C4jRxaY4OBg2QHYb7/9NrKzsyUzfqSsW3Kj7+VeB34QL2vmd5ZC3rVrV/Tt25f7fVdS6NWE+Wf0+eWXX9rVCmEtCI4yx9yFLS0gxwKjNHuibdu27jVOAY7axo6ie/bsyS2znZXdFv61UFpwz2w2CzL5HNX+EsOV+/DJJ5+0sxLwXUi1atWC0WgUZAU5IzMz061YGIvF4pWK8EuWLEFwcDAXu5icnCwrlZ/PO++8IxqXwlqC+RY5KeT2n76KlSIBIwMps61t5DU/WNa2UqvSlwr7EmB9tmygmG1QtLN9eSv+w500OqWdKd/NEhwcLPl71dXVskTDs88+C7PZDLPZjDNnzmDbtm2YPXs2/vvf/3LbSGUayI2+l3sdGJ7puW7durBarfjoo48cfmfPnj2cG5fFkavBVRzFWNnCACgtK8OJEyfs1pnNZsHUISxr164VnEeptHJnsPVYWNeko0QCpdkTalTalSvYo6OjnbbN1XPkynHwBw7sRJ5169Z1ejxRUVF292GDBg0cfqd27dqipTVsXUgAMHz4cIf74sNmdbkSC+PNjFSxtOekpCRFwvGee+6xW6aklMGsWbPwzDPPOBx4+6L2Cx8SMDLhBz+JWVjk7kPJS4VhGLsbTakFxhPKOCoqyq7jdCeNTmkUPz9FvqqqSlDMj4/cUV+7du24/xuNRvTt2xddu3bF559/7vB73bp1k30PuHIdgoODZVtuRowYgZkzZwqW8e9Zd6lfv75LVhN+4CgfsUDGkSNHonXr1tzn9PR0PPXUU4p/k4WNVZB6yW/YsEHxM7x//36X2wPUiBd2ihRnLF261GnbXBUwRUVFii0RfAHDBq7WrVtXMsCe5cKFC1i1apWg7/zpp58cfkcqRZo/eDlw4ACsVqvTMh18GDdiNrxpzRYb8BiNRkUCSiy5Qckx9O7d26kHAvB+7Rc+JGAU4MjCIhf+S2XixImKv69UwCQlJalaxGjGjBk4c+YMN+eVyWRyScyxWK1WxRPC8TNzTp48ifnz54tuJyeFH7AXF1arFZMnT3b6PbYDlYMr1+Hy5cuKpqdftGiRXb0L9p51lyeeeMJhMUEpxF6wjq45//4+evSoIveAUrp166boGbZarcjJyXHrN8+dO4cLFy7I2lbOQId/fhs0aCBbtDMMo7huE/+FyFq35IqHs2fPCvpO/r6mTp2K1atXC1KixZ5di8WCW2+9lfs8c+ZMyYBxZyi1TFutVtFgXU8hNeB55plnZNeFErs2So6bvQZyPRC+gASMD2BfKsOGDVP8XbkvZf5vLV++XPHvSBEXFwej0ci9zOrVq+eymGNTiuW8FPgds635WyotsWHDhrLijmzNn/n5+bL8wxUVFTCZTLJeAq5ch4qKCqxevVrRd5544gk7UaXEXC5VALJ3796KRroGABHh4WjWrJlgOXvNpVLn+ff3Z599pkp6v5SFQKnlLz8/X3J6AiXExMTISsmWk53HFzAPPfSQonYwDKPIncJ/jthCo3KtcvwXssViEcTuLF26FE8//bTgetj2c1JTULCxgkpRkpnI3rPPP/+84t9xBUcuGaPRKFlvybafE7PAKLEE86+BGh4IT0ACxoc4qy1hMBjshIFSCwxQ4yOeMWOGa420gQ3kY9tlW9FXLkpTu5XUe2AJCQmRFXdke46VjFLOnTsneySr9DpUVVUpTjEXm85errk8Ojoa586d4zqp7t27c+vq1KkjOwaGPdcDBw4UuPecXXOLxSLoND09k66nZqV2RtOmTbFkyRKn202bNs2puOC7VGJiYpCVlaUoTkKJO4V/vtiil7GxsYpiJBwJEX4AO/8+cBS3IbZMjjjMyMiQ9cy6U4JCyuJqW6zSlldeecXhgFCqcrnJZBLElokJGCXJJLYWVzU8EGpDAsaHOKqxwN5gtp2RKwIGABYuXIgNGzY4nGJBbmEr4N+Hw51gODmj67Vr12LevHkuBR1eunQJkZGRyMrKUhTM6kq8ityR7MKFCyXnD7NFicWDj+2LVu6L94EHHhBkPvBfBCEhIbLbYzKZsG7dOrRv355bJueap6enu+SmcpVdu3Z5ZFZqR0RFRSEpKUlWIK8zccG+XFmeeeYZTJ06VVAUVA5y7w9+/RNWwNSvX192jIQSIcIXMM7iNmy/K7duzfjx4x1ed3fn7srMzOQGA/xib3feeSc2bNggeQ9MmzbNobjiC5MVK1YILCL8LDmx/lxOTAuLUmu/LyAB42PMZjPuuOMOu+WsfzEiIoJbZjAY7EbBStITU1JScPHiRS7L5v/+7/+wePFirF69Grm5uVi3bh2XGi4FO+JjHyJXLDBKAskOHz6MuXPnKi4ECAAHDx5E7969MXXqVCxZskRg/nSEs6BEW5QGBt53332yrBlK51hhsX3Ryn3x2pqm+VWB69SpI0vAsB3p0KFDBcvlXPOTJ0/i+PHjstqqBpMmTVI8KzXf4uEKkydPhtFolC0apLZjxYvtnHHFxcWYN2+eoglA5d4ffAHDBjPXrVtXdoyEEiHCf3kqtXyZzWaHVcdZLly4gLy8PMn17gbtNm3alBsM9OrVi1seHByMlJQUSZeyswlG+cLE1iLCP4dS/bij69W7d2/uMwkYQhb8zAugJn+f9S/y1XadOnWcqmZnsFk2zz33HF544QVMnz4do0aN4iozOsuSYkfI7riQlHRIGRkZbsdBFBcX47777sPFixe5h10Kq9WKadOmufQ7So6LHfnZvjz4ndOJEycUm2nF/OdyUhzFvsd/Wefn5+PgwYNO9yNlWnZllB8VFeXw/lbTWiN3VuouXbq4/BtRUVF45plnAMgXDWLbObNkyO0TlKTAWiwWQUYZOwhgA5LlxEgoeT7419YVy5dca6EjAeOOy7B+/fqC88rvx4ODg2G1WiWnxmGvq5RVl/98uTppsdT1atWqFbeNN62hrkICxg+wHY3zXwL8B1Es+lztQmXsjS2Vnsne1O64kLxd9MhZh8DHnVGX3MDAqqoqrh0///wzcnNzuRmH3QnANRgMkkXPHKWaS32PL2DS09NdtggB8q85P9bpsccec7htnz59HBYfVILce4Q1+8udCZnPihUruHPsLBbBkbiQY8m4cOGCLEuMnBRYW1cVn2+++YYTfc5iJFydlkXOubJFan40JbjTT9keP78fl1MewZFVt6Kigvu/qwIGcH69yAJDyMJWwPBN97YWGEfI9eU7w2g0SgYCsje5Oy4kX0wQJ9fN486oa+7cubICA/nBkPXq1UNSUpJd+rMttp2LrSCJj493mNIo9ZJylArJFzDuWsHkXPP4+Hh069aN+9ynTx9s3LhR8jsdO3bkig8uXrzYrfYB/94jO3fulNyGfXksXrxYcQVpflCnO/U15N6jbdq0wZkzZ+wywVhWrlzpNItEbuySnH7HFSECKIvbYGnevLnT9gBwaI11ZwLPG2+8UfDZVsC440LkZ8J9//33gnPvbp/Kr9ZMAoaQhe1ogX/j8NfZWmAsFgt+//137rNcX74zLBaLXZaALe64kHxZ6v706dMOrRzuWofkdOZ8ARMSEiLL6mO1WrF06VLO3Pv333+7nNIo93tqmpDlTAqYmZlpF3djNpsFqbr80gN8d2Z6erpqopgNUBWDDSYvLCxUHCtl+zJytb6GEveT0WhEx44duWX82KSBAwc63Yfc2CU58V+uCBEWZ+fKlpYtWzptT1RUlEMBo3QiS/6gwrZPtxUwrroQLRaLoF5Onz59BH2+uwMNvoBRUnnbV5CA8TEWiwVvv/22YFnnzp25G1LKAsOadW3N+nJ8+c7aI7Zf/np+u1y1+Hii1L0cDh8+jISEBMEyfgfgjnVIzgge+HcUX7t2bUUBnY0aNeLMvbbzpDhzA/CPR+731BaXQ4YMwdy5cyUDYc1ms0CkswKF3/l36NCB+z//eVBz1mR25ncxWAEjVWHYEWIvLVfqayh1P/HdDLfeeiv3YpKT2edusLEtcoWI2LEpOVd8ISzF22+/7fS5MZvNmDdvntN9AcKX/2uvvSbog/n3cFBQkEsuRLZvtq1FxO/z+W1wxSKvpjXHG5CA8SHsDWnbGfJvSDELjJx0RFfm+1BiLnbHhcRiNpsFowlPYjAYEBUVhblz59qNKPnnWw3rkKMRPGA/r5U7AZ2eRGweI1dhi4FlZGRwgbpiFh5+R//LL78I7jVAWDjN9vtmsxlPPvmkW+2Mj4/HbbfdJrmefenHxsbK3qezYFml9TWUup/4AiYiIoI7x3Jimjxxb8oRIlJ9kNi5ErOoOsoWM5lMyM7Olm2xZOcDUkJZWZlgIMm/r9n6XkquoZw+f/z48VixYgW33BWLvO3kq/4OCRgfIVeE8Dvvy5cvIy8vD3l5eS4HgDlCibnY3UJ2ABxG4qsNe07liD6pUaKcUR3geAQP/Ctg2BewOwGdclEqxiwWi9NJJOWyadMm0WJg/GBE9jcHDRrEfX7kkUeQkJAg2I4fXG7ro5cz8aUzRowYISogrFYr8vLyuLTlzp07y7LUeWq+GCXuJ/45i4iIUGSBkRMHYjQaFRdcVKsoGiuM+SQkJODo0aOi2/fr1w+FhYWKKsgajUaXB1psn2IrYABl11Bu4LYj64wcSMAQspAbhc4fBR86dAi9e/dGamqqrN9QGpCqxFzMCqvq6mqX/a5KM34yMjIUVRnlEx4e7nAOGlvRJzZKfPrppx3+Bis0HI3gAXsLjDcmTFMiYFjLIDthnxJsO8rq6mpMnz5d1j2SkpJiZ70qLi4WVOTlv4xtLTBqTLa3ePFiu7l12Jdk7969uefx4YcfxsiRIwE4PreenC9GrkuFb4EJDw9XZIExGo3ccUphtVqRmprqduwd35Jy8eJFpxZkqSq5xcXF+PDDD0W/06BBA5eeo8GDByv+Dr9PERMwgPxr6GpygVKLPP+eUCspxJOQgPERcm9IsZGE3KJuSl0OSszFfMuQqze53HMQFRWF7OxszJ07F++8847TYntilJWVKW4Tf5R48eJFh7MIKxEatgIG8PyEaXwB4ahjcrf6qG1HeeLECRTJmFfKto1SyxxZYNQq9f/kk09yI1Ep69G5c+ewePFiPPnkk3bXLCYmBunp6V6ZL0aOJUPKhSTHAmOxWGRleCmdV0nsd/iWlL179zp0fyidXoBF6fQRUt+rW7euZGkCW06fPi1o06FDhwTnSc41dMd9LNcib7FY8Nlnn3Gf1UoK8SQkYHyEJ+MZXHU5yE11TUpKEjxkrnZacs/B+vXruZeAp4N/xdrEjvSkJo0ElAkNMQEDeG7CNIvFInDDOOqY3LVi2HaUrgS7OoJ/zmwtMGo8UwzD4GRREU6cOCHLerRu3TocPXpUcM1Onz6NpUuXuuUaURO+6/P333+XbYFRKmZdcVsDji0pUu4PV+/TnJwcl17Itu7Oq1evyna32CYOLFq0SLEwUKP0hCOBL2V1dTcpxOMwOqW0tJQBwJSWlvq6KaJUVVUxJpOJMRgMDADV/gwGA2MwGJjs7GyX2pWdnc3tQ2z/7H7Ly8u5ZeXl5eI7q6hgmIyMmr+KCslz4OhY4uPjmaqqKtHv5ubmMmvXrmUmTZok69zExMRIHpfUbzlrIwCmTp06TAX/+Jwc92effcYAYG6++WaZV8V12Osp9z5Zu3at2/fg2rVruXOQl5zM1Fbx/t62bRv3/82bNwvartYzVRtg9pvNzLExY2S1PTc31+PX0VWys7OZyMhIQXtr1arFAGC++eYb4cY2921ubq5r114BrvYB7tynkv2jg+f22LFjLv1OVFSUoufPEc76ZlfvU3f6YU8h9/1NFhgfwcY9MG7m7dvGhLjrcnBk4QgODub2y3chuRrIazQakZmZKbrOmUuGb3a97777HP4Oa5Fi5x5REmciZ6R37do17Nq1y+E2ttsDni/V7Uq2mhpWDHYfBw8eRK6DUu2uwI8lEJstV61U6rCwMNnWI7VcV2rDjqptXc7s87p9+3aH33fluJTeP65WpPVGvSY+thYYOTAKEgfk4MjN7GjaDWcWeXeqAvsaEjAaJysrS3WXg60rg4X/8lDDhQTYTyDIokSIOSq4xBcnUnM9OfottetgANIuJLVxpWNy11RtMpmQlJQEq9WKL774wqV9OMJZZWq2k3e1gqrBYEC8yYRmzZpJTqdhi7fT2+Ugx/2zfPlyh8+u0uOKiYlR7LZ29flSo16TkheyUgETFhaGefPmKUockIOYm7mwsJCrJeZKEoAn+jhv4f6EEYRLsB2MI8RqHLAYDAaYTCaP+dlZCwcA3H///QCElSb5v+lOKjXfD79lyxaUlZWhcePGdnE2ztrKEhkZKRhxmkwmZGZmCmJohgwZgvz8fJw+fdrpb3miDoa3BIwrHRNrxUhJSYHBYLCb3Zb9bLuOZf78+TAajdixcyfKXIh/cfSbgGMLDAv/Gm/evBlr1qyRXTWXYRgsXrwYQQcPolmzZjA1bYrCU6dEj5V9Bt1Jb/cUciyH586dQ35+vmQ1WlYkyI01Wb58ueK+yNXny9F9KhdXBh1ymTNnjuw4PaXCgN83s7DCfcqUKYLrZdv/ieGvdajkQBYYHyG3fLwYnqot4Qy+gGGLMQHqCZg+ffq4VBeCv+3y5cudWqSU1KCQM9Jj5zOSi7cEjKsdkyNTdXZ2NrKzsyU752eeeQYWi8VpMT8xMjIyRH+TD98C48gFx17jpUuX4vTp04IJM+USFBTEZZ55Kr3dU6gxqmZFghwrx4wZMyQnfHSEOzWQ3A3oV/JCVmqB6dGjh9eFgatJAN6oQ+UpSMD4CLkdTKdOneyWebK2hCNs0wb5VSJdhS9g+KNrJfBfIEpL7MvZt7PKvD169FD0O2xn6GkB4+7LQaozdFRBuaSkBCkpKbJn5uYzb948MAyDefPmCX6Tj5LJTVmMRqOsCTOBmnPCT6MeOnQoNm7caDdxY1RUlCA7zt9Q6+XJigR+BWQ+MTExyMrKwsKFCxW3EXC/BhL/Pp09e7bs31X6QhYTMI6E3cGDB9GrVy+vCwNXCgR6ow6VpyAB4yPkdjCtWrXi/j9x4kSv1JaQwlbAqDGdAF/AuPqA8F9qnnjInI30pGb7lcJbFhh3OyapztBRBWXWlP/ee+8hPCwMSiMUTp06hblz5yIkJES0A+bfgz/99JNs8Sw37ZbhpVGzmM1mTJ48WbDd+fPnMW3aNL9NL5VjOQwLC0OvXr2c7stsNnPFM1NSUrB06VKsXr2aSxkfPny4W211twYSe5/OnTtXdnVkpS9kWxfSLbfcYtde/r35+OOPo1WrVpLFDv1NGHi6DpWnIAHjI+SOjvkvx/79+3u9tgT/BVFVVSX4rKaAqV27tsuBo/zzYTsLrFrwR3qvvfaaYJ0S/7vVauVmED9//rzHK116omOSExxcVFyM7t27A4AiEeMsO6Nfv37c/0eNGiW7nobSOAN+BpLFYhEd3ftzjQw5lsPLly+jVatWstrPvsB79eqF9PR0jBo1StW+SI0aSHLnMXPlvre1wLRo0YJrL+uatK0LU1xcLFns0B+FgafqUHkSEjA+Qu7omJ9hI2XG9RS21TFLS0sFLww1XUiuuo/47bD9v9qwI71hw4YJlsspyQ78ez4/+OADAMDmzZu9UulS7Y5JrhiIjIxEamqqw5GqGI6yM8SmGpAjIpTGGbAZSI6C7V1JhfUmcmJE5J6/K1euAJA/H5grqDE/ktQxN2jQAEBNP+PKfW8rYIKDg526Jtn7Q6zYob8KA7XmqPIWJGB8iJzRsdRMvJ5GTnVMtS0wruINCwwf29+QE+DnSrVRNVGzY5IrBsLCwtC+fXu7eBi5Ljc5QkmuiJAzKSEgTKMGgJ07d2q2RgZQ08ccPXrULoaHRe75YwWMo1me/QUxwf7bb78BqOlvXMlYsnUhsQHkcksV7Nq1S1PCQCuQgPExzkbHvhAwcgugqZmF5I6A8XQMjKPfA5wLGFcKyvkzstyfPBFge02U1FeRc07kiAglGTWLFy/mrESnTp2S1VZ/rJHBsmvXLodp5Oz527lzp+Q23rDAqImtYOcLLzlzQNli+4yzAkbLNVT0AAkYP8DR6Jj/YveWgJE7qmDFx2effebyzKWs+NGyBaakpAQfffSR5DnQcqVLMeTEGvBFgC1NmjSRFf91/vx5gQvTGc5eEs6K3MXHx2Pjxo0YOnQogJpMkunTp8v6bX+skcEi9+XpKPVdawLGFn7KvStVdaUsMFquoaIHSMD4OfzAMG8JGLkd3l9//QUAmD17tsszl6rtQpI7Q6w72AqYH3/8Effffz93DjZt2iRYr8dRmpT7MyYmRiACxKhXr57T+K8RI0YgNTVV0YR9cl4SfIvn6tWrMWHCBAA1QZl8y+fBgweRlZWF8w4qqbLt9dcaGSxyX55xcXGS67QuYPixhK4IGLEYGEDbNVT0AAkYP4f/4Hir85Db4YlF3SuN51BbwLgzW6tcHFl5iouLMWLECBw8eJBbptdRGl8MtG7dGgCwYsUK0eBE/vWtV6+ew/iv9evX46OPPpIdq6D0JcFaPEeNGoWxY8cCqLkP+UHpX3zxBZz9ur+lwkoh9yV72223Se7j6tWrALQrYIKCgrh7MC8vz6HFVAwpF5KWa6joARIwfg7fdOlozh81cXWeEVfiOdSOgfE17Dn48ssvOYGn51EaKwaaNGkCQDojLTQ0lPs/G48gFf8VExMj2/Li7kuCzU65dOkSt2ynzGkQoqOj/S4VVgw1XrJat8AA/w50UlNTBRZTJWnkLHyXlFZrqOgBRQJmwYIF6NGjB8LCwhAbG4t7770Xhw4dEmxz7do1pKWlISoqCvXr18ewYcNw5swZwTYnTpzA4MGDUbduXcTGxmLGjBl2gaB5eXno1q0bQkJC0Lp1a6xatcq1I9Q47MgHgMtxJkqRW09BDKXxHKyAcUeEeHt08+233zpczwAoLSvjiqEFwiiNFddSAZJ8AWMbdG0b/6XElebuS6Jhw4YAgPLycu5elDsNwtKlSzXzcnLnJcswjKaykMSwWCyi8xnJtRpLuZBYtFhDRQ8oEjA7duxAWloa9uzZg5ycHFRWVqJ///7czQ0AU6dOxaeffooNGzZgx44dOHXqlOAiWq1WDB48GNevX8euXbvw/vvvY9WqVZgzZw63zbFjxzB48GD07t0bBQUFSE9Px6OPPoqvvvpKhUPWDhaLBatXr+Y+uxpn4gpSHV5kZKSs78t9CalhgeGjpDqrq8g9Nn4xND2P0qxWK3es+/fvFz3//BHriRMnVJkFeenSpW6/JPhxZU8//TS2b98umXJsi6tz8PgK9iU7btw4AMDAgQNlnb/r169z10uLFhh3a/lYrVb8+eefgmVi83BprYaKLmDc4OzZswwAZseOHQzDMMylS5eY2rVrMxs2bOC2OXjwIAOA2b17N8MwDPP5558zQUFBTElJCbfNm2++yYSHhzMVFRUMwzDMzJkzmRtvvFHwW/fddx8zYMAA2W0rLS1lADClpaUuH58vyc7OZgwGA4OaAT33ZzAYGIPBwGRnZ3ulHVVVVUxubi6zdu1aJjc3l9m2bZtdm8T+cnNzGaaigmEyMmr+/rm2tvudPn06A4Dp0aOHS+3Lzs5mmjZtKvhtk8nk0fOTm5vr8NhrA0wGwBwbM0b0uJs1a8YAYJYsWcJUVVV5rJ3eIDs7mzGZTHbn37JuHXftLevWMUajUfY1qqqqYkwmk+j9zz4D8fHxbp+77OxsJioqym7/jRo2ZF4ODWXm/nMtPfX7vmLJkiUMAOb+++8XrpB4Xtl+HgCTk5OjueN29ryyf3lbt9odv9j9DYCZMGGCT49J78h9f7slYA4fPswAYH755ReGYRhm+/btDADmr7/+EmzXrFkzZsmSJQzDMMyzzz7LdOnSRbD+zz//ZAAwP/74I8MwDJOUlMRMmTJFsM17773HhIeHS7bl2rVrTGlpKfd38uRJzQoYtgOXetB82YEqerlIdIhinUJwcLBi0eErkVdVVeWwIwwGmCXh4Yz12WftBAzDMExCQoJA1GsVR+c/GGB+S01lfktNZYIl7hNH14jdt+3+1bq22dnZTgVoxj/XUqzt3hpAeII33niDAcCYzWbhCpHnNTs7m2ncuLFXBwhqk56eLkvArPvgA8HxS93f7J+WzoHWkCtgXA7ira6uRnp6Om699VZ07NgRQI3vODg4mAuMY2nUqBHnVy4pKUGjRo3s1rPrHG1TVlaGv//+W7Q9CxYsQEREBPcXHx/v6qH5HH+uG8KP57BFTjyHVEXa69evK8pg8mVxOEemYfYcDBw4UDKlm72H+XEhWkPO+f/iiy8ks3mcXSNPutysVqvd5Ixi1A0N5YKTWRo0aKB5lx87iaizdGL2WbV1mfrzHFC2WCwWZGZmytqWn0bu6P5m0VLxSb3isoBJS0vDr7/+inXr1qnZHpeZNWsWSktLub+TJ0/6ukku4+91Q9iXi20gm7OXi5qdgr+KPJPJhHXr1qF9+/aS2+hBwDg9/wDKLl92mM3j7Bp5KjAyPz8fxcXFTre7+vffdhN3Llq0SNPiBfg3fsORgNFD9WhHsS98xNLInU0hAUBTxSf1iksCZuLEidiyZQtyc3MFVS3j4uJw/fp1QUoiAJw5c4ZTt3FxcXZZSexnZ9uEh4dLdvohISEIDw8X/GkVLdQNMZvN6NWrF/e5du3aWLBgASIjIyU7NTVFhz+KvA8//BDHjh1zWMQN0IeAUfO8OtqXJwIjlbTdVugcP37cr1/acpAjYLQ+BxTgvL9hYRjGzmosNxNNS8Un9YgiAcMwDCZOnIhNmzbh66+/RosWLQTru3fvjtq1a2P79u3cskOHDuHEiRNITEwEACQmJuKXX37B2bNnuW1ycnIQHh6ODh06cNvw98Fuw+5D72ilbgg/c6iyshIPPPCAw0wpNUWHP4q8Ll26OH3BWq1WLvNKywJGzfPqbSGu5Pf++OMPwefnn3/ea5mAnkKOgNHDC1xu29LT0+2sao6qEvPRWvFJvaFIwKSlpWH16tVYu3YtwsLCUFJSgpKSEm5EGRERgbFjx2LatGnIzc3Fvn378PDDDyMxMRH/+c9/AAD9+/dHhw4dMHr0aPz888/46quvMHv2bKSlpXEP1oQJE/Dnn39i5syZ+P3337F8+XJkZWVh6tSpKh++f6KVuiFSE8SxPnLbkvpqig5/FHl33nmn0xcbP4ZLywLG6fkHEB4WhvCwMEhVEvKVEE9KSpKVAl03NBSvvf663XItxYCIIUfA6OEFLrdtQ4YMsVt22223OS3maTQacf78eZfbR6iAkshgSERjr1y5ktvm77//Zp544gmmYcOGTN26dZmhQ4cyp0+fFuynsLCQGTRoEBMaGspER0cz06dPZyorKwXb5ObmMjfddBMTHBzMtGzZUvAbctB6GjXDiGfrxMfH+0X0e1VVFVOnTh2HmVItTaaaTJx/ovrVTo/1dKaKo991dNz8FGLbLCR+SqrVavVI+7yFo/PvbhaSN9oudQ3ZLKSXQ0NF06hduVf9iby8PAYA07ZtW+EKXhZS1dWrXkll9yTO+hu2P+WOwSYLy1kWEnse/KE/1hteSaP2Z/QgYBjGvg6Lv3QYcmor1MY/tVBs0jIddYqupFJ7U+TJSXG3FW58jh8/zgA1aeN6QOr880XchrVrRV8cvu74perAxEVGMnnJyUwGxOvA8P9yc3N9egyusGfPHgYA06hRI2G/IvMF7mvxqQQpkc3+CY5BJI08KyvLroaR1oScFvF4GjXhHfy1uqMr1WiBmuDfrKwsu+3q1avnUnqqt0t4ywpELiriphGwRQ8BvHzY87906VIAQHx8vF0g812DB3P/f/fdd/2mzLrZbMaZM2ewbds2tGnTBgAwffp0FBUVqV5x2p/YuXMngJrECEezqLPZhux0Cyxaqh4tlY4PAOPGjXN6DDExMQ6DthkNBDPrGRIwhEvI9S+HhYUJPlssFtFYpsTERJc7RG+KPFeFG4veBAzw7/kHhLM6s7BBywD8TogbjUb07dsXN910EwCgRYsWMBqNdvetFP4cAyKGxWLBjBkz7JaLzaIO1AgAdpqX//znP34jPpXAiuyZM2cKljsqdcDij9mOxL+QgCFcIikpyeFL2GAwIN5kQrNmzbhlUkXsAGDbtm2aCIp0Vbix6FHAAOKzOrPwJ2pVa84rtWGLu7ET/jVr1gzhDkSMv2QCKkFObRf+LOos7ASdbdu29SvxqQSj0WiXxSpVFJWPP2Y7Ev9CAoZwCaPRyI1abWEj9xcvXsxVo9VLZUs2+0YKMeHGR+8C5tq1a3az/vItMP768rMVMEFBQRg4cKDotv6UCagEOcUH+bOos7DnhD1HWsXWFSY2O7Utsp53jQlZPUEChnCZ2NhY0eWsj5wfB6GXypZGo5GL95CCL9xs0auA4c9S/PnnnwuEKGuBMRqNDtNSfYmtgDl48CC+/PJL0W21FAPCx1X3p14EjG1x06tXrzr9jrtTpxCehQQM4TL8kTVLamqqqI9cD4WxAOkYHgCIjo62E2626FHAWCwWtGzZkvs8bNgwtGnThounYO8Tf3UfAf/WRrl27Ro2bdqErKwsyWkQXnnlFc2JF8B196ceBIzFYsFgXjA5ALz11luy3NZmsxlRUVF2y7UqZPUECRjCZcQETPPmzUVHI3ILYx0+fNjtdnkKRzE8APDggw867cz0JmCkzklxcTGysrJw8OBBzgJTq1YtXzRRFuzL+erVq5g+fbroBJRAzah7+vTpfu/qFENO8cGI8HA796fWBYzUpJTl5eWyCxKy87698847Xsl2JORBAoZwGTa4Tw5yKlsCwNy5c/0ymFdODM/777/v9MV25coVAMCFCxeQl5enyRchi8Og0H/+/fLLL7mKr/5sgWFfzoWFhShyMNGjltNm5bhDxGZR17KAUSv2jhXhiYmJfpdJF8iQgCFcRswCIyVQ2M7TUUfC4o/BvHImhrtw4YLDF5vFYsGsWbMAAD/++KPDeaO0gNyg0B9++AGANiwwpaWlsrb3d1enFGxdFFsczaLOChgtWg3VmkCW7ev8+R4OREjAEC6jRMAANZ3nvHnzHO7TX0e47taDYM3YtmnGWp5XR+45YeOftGCBkTuq1nLarK3bY8KECQ5nUdeyBUatOi6sBcaf7+FAhAQM4TJiLiRnlhO24qkz/G2E6049CDn1N/zR6uQMueeETV/159Er+3IOCwuDqWlTv5uA0pM0btzYoXDTsoBRq46LFuK4AhESMITLiAkYZ3ExWi0M5SwAEqh5UYu92NQyY/sbcoNCb7zxRgD+PXplX84VFRV45ZVXAMBOxOglbdZWKEtVjWbRsoBRa9Z6EjD+CQkYwiUsFotoxtDvv//u8HtJSUlc0TMx/HWEyw+AlOoMhw0bJvpi02s5ckfnhP00cOBArrKrP3f+/DowQ4cORWpqql06sR7SZi0WCxISEgTLVqxY4dCFyWbOaVHAyHlunQlShmFIwPgpJGAIxbDxHPwS8Sxbt2512BkajUakpKSIrvP3Ea7UxHB169YFAHTs2FH0e1q1OslB6pyYTCakpqaiffv2muj82Zfz2bNnsX79eoSGhmLKlCm4sUMHAMC8efM0nzYrlfLOphPbTubIomULDCB9j9aqVUuWIOVbrPz5Hg5ESMAQilAjLbFTp06iy7UwwhWb/ZoNfhQLagbUM2P7K+w5SU1NBVBTzPCPP/7gMlq0UMiOzZT6448/MPrBB7Hq/fexbNkyLgX8pptu8ktRLRc5z+2TTz5pNw8SoH0BAwif24yMDAA1019ERkbKTqEG/PseDkRIwBCKkJNO7CyeQ8xys3DhQs2McG1nv2aLXIkdF7u93suRG41G3HDDDQBqppjgH4u/W2AsFotodtzly5dx5OhRANp/ccmKwyoqspsHCdCHgAFq7tGLFy9ixYoVAIDz58/LKmXAf6799R4OVEjAEIpQI55DbMRzyy23aPYFzr7cpAQM8K8Zm3U3sWjB6iQX2/mEWPw5BZW1TIjBt1X46xxOcnF1HiRAPwKGdaGdOXNGsJwtZSDlQuNbVknA+BckYAhFqBHPISZg2LlotAjbqTkSMECNiLn77rsBAKNHj9ZdOXK20Bkb9Mniz0XA5FgUAefB6f6Oq/MgAfoQMHJKGYi50KxWK3bs2MF91rqQ1RskYAhFyEkndhbPISZgtNw5yhUwwL/H3rNnT92VI5eywLDH7I8WGLmWCdsChFpDVhyWyWQ3DxKgDwHjiguNzdjiF/hr0aKFJotO6hUSMIQi1EhLDFQLDKCNgFZXkRIw/myBkWuZiI2N9XBLPIuc53bx4sV28yAxDKMLAaPUhbZp0ybJSUq1Wjlbj5CAIRQjlZYIAK1bt1aUlshCAkb7aDEGJikpCVFRUU6369Gjhxda41mknluDwYC5c+finnvusfsO3x34ww8/aK5aNIsSF1p1dXXNjOQ6q5ytR0jAEC7BT0vkU79+faffFXvRk4DRPmwMjJYsMJs3b8aFCxecbqfl+5MP+9zOmzePKyjJMAwyMjLQpk0bHDx4kNvWYrEIpv64++67NTv5qBIX2okTJ3Q7I7neIAFDuAybTszH2VQCQGDHwPjzy9xd2GtoG8Trr2nUjjKQbNGT4Ny8eTPmzp0rOrFoVlYWDh48yLlQTp06ZbeNFl0oDqtG//OZdaE5m1qBRWuVs/UICRhCVVwVMFoe4bIvZqlCdnz82Z3iLs5iYPztmOVmIAH+13ZXcZiN88+/X3zxhS5dKFIutOjoaKxfv54L1hXLxBJDi5Wz9QYJGEJVAlnABLoLyVbAVFdXo7CwEN9//z0A4MSJE8jLy/ObF5+SEbS/WY9cxWk2DoCyy5d160Ixm81YunSpoB7TuXPnMG3aNK4OTLNmzWpmJNdp5Ww9QQKGUBVXBYyWXxByCtmx6FnA8GNgNm3ahMzMTKx6/32sz8oCAOzatUtW5VNvoWQErZfrpabbQ4suFIvFgtTUVFy9elWwvLi4GCNGjMDBgwcRFBTEzUhui14qZ+sFEjCEqsgRMGIvei0XiCILTA2sBaa0tBQjRoxAmUQsQVFRkV/EUTgN7OT9X8sCm4+abg+tuVDkFLP78ssvUV1djaFDh2Ljxo12aeV6qpytB0jAEG5ha02xjX+Q8x2tQwKmBlbAXLlyBdJTBtbAMIzP4yjk1EZh0cv1kiPawsPCdOlCkeM+Ky0r44rZmc1mxMXFceubNWumq8rZeoAEDOEybKVKPlevXnU6sg5kAeOvGTlqoDSTzB/iKKQCO+vVrSt4UelFwDjMxvnn30GDBunSheLKfFD85zQiIkJzx6x3SMAQLsFOjCY2onHmHghkAaNnCwwbA6MEf4ij4Nc0MplM3DJ2dm1AX4JTSrSZTCakpqaiffv2nAvFtsifll0orswHxXch6eke0AskYAjFOPIlszhyD5CA0aeAcaWWj7/EUbA1jWJiYgDUuLj4E/vp7Xqxoi0+Ph5AjUXljz/+QPv27QXbvPrqqwCADh06aH7yUTnus4jwcMF8UHyLCwkY/4MEDKEYWROjnTyJnTt3iq4nAaO/FyIAxeb1mJgYv4ujCA4OBlBzj/LvUz26DoxGI1eN98YbbxQ9RjYov3nz5pqffFROMbuBAwcKrC78/+vxmdU6igXMN998g7vvvhtNmjSBwWDAxx9/LFhvMBhE/xYtWsRtk5CQYLf+pZdeEuxn//79SEpKQp06dRAfH4+FCxe6doSE6sg1+3/66aeiy+W86LWEkkJ2ehUwFosFLVq0UPSd5cuX+90Lka1HZLVaOQtMrVq1NJ0l5whWsEllD+phIkc+jtxn69atE1igALLA+DuKBcyVK1fQpUsXvPHGG6LrT58+Lfh77733YDAYMGzYMMF28+fPF2w3adIkbl1ZWRn69++P5s2bY9++fVi0aBHmzp2Lt99+W2lzCQ8g1+z/6muvCeZWYQlEC4zVakVeXh5XZl9PL0RH8VBSzJgxAykpKR5slWvwLTCsgNGb2OTDCrZAETDAv+4z9phWr16NY8eOcZV4+fAFjG1KNeF7FEvKQYMGYdCgQZLr+WlnQM28G71790bLli0Fy8PCwuy2ZVmzZg2uX7+O9957D8HBwbjxxhtRUFCAJUuWYPz48UqbTKgM60t29sIyoKauQtu2bQVKWW8CxlkhO4vFgilTpgjOV3JyMl5//XXNxhOwyImH4hMTE4M33ngDw4cP93DLXEPMhaTnkXegWWBYjEYjwsLCcO3aNXTu3LlGqIj0S3zRord+Sw94VFKeOXMGn332GcaOHWu37qWXXkJUVBS6du2KRYsWCTr/3bt34/bbb+ceLgAYMGAADh06hL/++kv0tyoqKlBWVib4IzwD35fsCNu6Cix66wgcWWDYSfFsxV5JSYlfFHNzF7nzCc2ePRu5ubk4ffq034oX4N8XelVVVUBYYAJVwABAvXr1ANR4FaTgW2D05vrWAx4VMO+//z7CwsLsRpmTJ0/GunXrkJubi8ceewwvvvgiZs6cya0vKSlBo0aNBN9hP5eUlIj+1oIFCxAREcH9sdH1hGcwm81IT0+Xta3t7K5iAsaf5shRipSAqa6u1uWkeHzkxkN16NBBE0GggeZCYo+3oqJCdD0rYFxJkfd3WAFjO60AH/6zKSfGjfAuHhUw7733HkaNGmWn3qdNm4bk5GR07twZEyZMwCuvvILXXntN8iGSw6xZs1BaWsr9nTx50t3mE04YMmSIrO1sZ3cVG8n40xw5SpESMCdOnNDtpHgscuOh/CVd2hliAoZcSPq0wLATOjqywPDPC1lg/A+PCZj8/HwcOnQIjz76qNNte/bsiaqqKhQWFgKoiaM5c+aMYBv2s1TcTEhICMLDwwV/hGdxpa4CID1qLy4u1qRbRUrA2FqepPCHYm6u4uweAIC6oaF+ly4tBd+FxLrGqqqqNG0lc0QgCxg5LiT+eSELjP/hMQHz7rvvonv37ujSpYvTbQsKChAUFITY2FgAQGJiIr755hvBDZOTk4O2bduiYcOGnmoyoRBX6ipYrVYcOnRIdH9adauwAqasrAwfffQRduzYgerqajvLkxRasU6IIWc+oVtuucXvXUcs7Av922+/xWeffw6gZvCkVeugM5xlIbFZc3oWMI5cSGSB8W8UC5jy8nIUFBSgoKAAAHDs2DEUFBQIAjXLysqwYcMGUevL7t27kZmZiZ9//hl//vkn1qxZg6lTp+KBBx7gxMn999+P4OBgjB07FgcOHMD69euxbNkyTJs2zcXDJDyF0roK+fn5Dl2FWnSr7NixA0CNBen+++/Hnf37IzMzE1euXNHlpHi2SN0Ddf55ObZq1coXzXIJ1vVcYfNC16p10BmBbIEhF5L2USxg9u7di65du6Jr164AauJZunbtijlz5nDbrFu3DgzDYOTIkXbfDwkJwbp163DHHXfgxhtvxAsvvICpU6cKarxERERg69atOHbsGLp3747p06djzpw5lELtp7B1FTp37gwAyMjIkKyrINddohW3ysyZM/HMM8/YLS+7fBkbN27EfffdJ/o9LU+KJwZ7DyQmJgIApk2dip49ewLQTgyJ1WrlxKgtWrUOOiOQBQxZYLSP4p4lOTnZac2H8ePHS4qNbt26Yc+ePU5/p3PnzpoahQc6RqMRcXFx2L9/P1q1aiVZV0FPQZ8bNmwQVJi2hQGQlZWFrKwsjB07VpDabzKZkJmZqfk6MHyMRiNnhWnevDmYvXsBaEfA5Ofno7y8HFI5R3zrYHJysjeb5jHkZiHpUcAotcBQDIz/QaUFCdVgOzm20xMjKSnJYVqqVtwqVqsVTzzxhNPtThYVITo6WmClCQkJ0fSkeI5gYyquXbumuTRkvVkH5UAWGMcChi/srl69qivrmx4gAUOoBlsrgg38E8NoNMJkMomu05JbJT8/H+fPn5e17enTpwV1NBiG8fvjcxVWwFRUVGguDVlP1kG5BLKAYZ/Jn3/+WbQO1caNGwUzkl+6dEm3wdxahQQMoRpyLDDAv3VhoqOjBctNJhM2btyoCcuEklF448aNBVYIqZeFHmDvgevXr2tOwCQlJSEiIkJyvVasg0oIxLmQgJrpPdjsue3bt6N3795o06YNN3fbpk2bkJqaavc9vQZzaxVt9CyEJpArYNiRztq1a1G7dm2cPn0ajRs3RlJSkmYsE0pG4efPnxeM5PSMli0wRqMRKSkp+ODdd+3Wack6qIRAtMCwk4/axnIWFxcjKysLKSkpmP6//0lW0DYYDEhPT8eQIUN0dS9oEbLAEKqhVMAEBwcjOTkZI0eO1ESZeT5sATc5TJs2LWACAMUEjJau6y233CK6XEvWQSU4EjBWqxWXLl0CABw4cEAX8R+OJh9ll3z++ee6r6CtF0jAEKohJwYG+DcdUSsjczHkTmgJ1NQW+eOPPzzcIv9AyxYYAIIJZFkSExN1G3QtlYW0adMmJCQk4OzZswCAxx9/XBfxH84mH2UAXHGQVs1HT8HcWoUEDKEaSi0wWhqZi6FkQkt2JKt3xLKQtCRg2PbzCQsL0/y9KoWYBebgwYMYMWKE3YteD/EfaooOPQVzaxUSMIRqBJqAAeRPaMmmbOodPVpgHBU60zrs8ZaUlGD9+vX4888/8cUXX0Cs0pceivnJFR0x0dG6r6CtB0jAEKrBChhnLiQ9CRg5E1rGx8cjPj7euw3zEVrOQgLE70lHdUK0DjslzN69ezH6wQfxwYcfoszBJKRaj/+QOwHtq6++Kr5ep8HcWoUEDKEagWiBcRQLw3aRmZmZTqtX6wUtW2AsFgvGjRtnt/zMmTM+aI3nsVgsWLJkiUvf1Wr8h5zndeDAgRg2bBhefvllu230GsytVUjAEKrBBvE6EzB6COLlw05mGBUVJVgeHh6OdevWwWw2a9bkrhS+gGGPWQvXmU2tFStOeOrUKWzYsMEHrfIcbDaOq2g5/oN9XmNiYgTLTSYTUlNTuQlob7vtNgBAkyZNsHbtWuTm5uo2mFur+H/PQmiGQLTAsJjNZgQFBWHo0KFo1bIlHkpKQrNmzRD0z4SWthPB5eXlaarujVy0GMTrKLWWZeTIkTAYDEhJSfFiyzyHs2wcKQwGA0wmk+bjP8xmM6Kjo3HHHXcgLi4OH330EZJ69oSRZ3Xhu8K1VqcqUCALDKEagRgDw4cNiIyIiEBCQgKCgv59vH777TfBtr1799ZFWqotWnQhyXmZW61WDB8+XDfXyxUXkN7iP9h7NTQ01K4O1aZNm7iZ5E+dOqXb51XrkIAhVCOQLTDAvy9qW2uLxWLBxx9/bLe9HtJSbWHvAS0JGCUvcy1n4PBxxQWkt/gPqeeVTSO3dSfq8XnVOiRgCNWQGwMTSALGUayBHtJSbeHPrcNV4vVzAaPkZa7lDBw+zipJGwAupqtVq1a6jP9g+x/+81pdXa3rNHK9QQKGUA1nFhir1Yq8vDyuaJZUKqNWYQUMf96jnTt3Oq78qfG0VFtEXUh+LlSVTAsBaDcDh4/RaMTIkSMdbjN69GgAQEJCguam+pCD2IDjxIkTuk4j1xskYAjVYGdc/uuvv+ymp2dLk/fu3Zt7sf3nP//RlTlWrEMsKSmR9V09vBQBbcbAKJkWAtB2Bg6LxWLB4sWLJddPmzYN7dq1AwDUr1/fW83yKux9ye+nLjsQL3z08rxqHRIwhCpYLBYMGjQIQE0nwJ+eXqo0eUlJia58ymIm6bi4OFnf1cNLEdBmFhJQk5WSlZUlCLy2RS8VWOVkXWVlZaG0tBRAzVQKekRswCH3WPXyvGodEjCE27A1NGytDcXFxViflYVPP/00IHzKYh3ibbfd5jjWQCcvRRYtWmBYhg8fjvfee090nZ4ycORkXZ0sKsKBAwcABJaAadasGcLDwiDl3Nbb86p1SMAQbiFnevqrDtKq9eRTFusQHVb+1NFLkUWLWUh8hg8fLrpcTxk4ct0f586dAxBYAiYoKIizJNuix+dV65CAIdzC1YJYtujBpyzmUwdq3BO333673fZ6eimycC8Fq5UL5taSgLGdzHHmjBm6y8CR6/5gX9h6FTBiLl8AaN++PdatW4cGDRoIluvxedU6JGAIt1BLeOjBpyxVVwIAmjZtCgBIS0vTbVlyi8WCbt26cZ8vl5cDAL799ltfNUkxRqNR4D7o16+f7jJw5ExoaGralJvEsqSkRBcuXlv4WYO2FuShQ4di1qxZAGrcwHp8XvUACRjCLdwVHnryKTsSMOyytm3bYuTIkbp7KbJxUKdOnbJb9+ycOZoJ1DYYDFwcD6At65Fc+G5NWxFjQI3r9++//8aOHTsAAG+88YYuq9Dyr62YQGNdoG3atNHd86oXSMAQbiFnNFc3NBQGiHSWOvMpS5mkAf0W7wOcZ7UYoK1Abb4biS0NoDfYCQ1ZyyALK94uXLwoWK7HKrR8ASP2zFZWVgLQ7z2gB0jAEG7hbDQHAHfffTfWrVsnOvurnnzKjiwwWpqZWSnO4qAYaKuCLV/A6PF6sZjNZhQWFiI3Nxd169YFIH28essYBJwLGHaZnu8BrUMChnAbqdEcf3r6oUOHYtOmTYL1evMpy3Eh6dECIzcOSiuB2oFggWExGo1ITk5GvXr1AADl/8S9iKGnjEGALDB6gAQMoQrsaC4iIgIA8O677+KPP/5A+/btuW1sR256e5mzHSID2LlT9OxCkhsHpZVA7UCxwPCxzb5yhFaEqDP4zyJZYLQJCRhCNYxGI1d2vEuXLnYv64qKCl80y2vwOzr+fEiAvl1IcuKgtBSoHYgCRomVQStC1BlBQUHcPUsWGG1CAoZQFbbzZx9+Puwkjiy28yVpHb5gsxUwenYhOYqDYtFSoDb/hRUoLy/2ua1Xt25AVaGVqt0EkAVGC5CAIVSF7fBtxQoALi2TpXfv3rpKzwxUCwwgHQcFACtWrNBUrFMgW2B69OgBAHYiRm8ZgyyO4tbIAuP/kIAhVEXKArNp0yYsXLjQbns9pWfKETB66vxtYeOg7rjjDsHyu+++20ctcg1+HZhAeXmxx8kG3osF5OspY5BFTuB9oIhYLUIChlAVMQtMdXU1pk+fLrq9ntIz+eLENohXzy4kPkajEa1btxYs09oLgH+NtNZ2V2EHHlarFe3bt8fhw4e5Z/mjjz7SXcYgi6PaTWSB8X9IwBCqImaBOXHiBIqKiyW/o5f0zKCgIAQF1TxSgeZC4mM7d46SY/7mm29w9913o0mTJjAYDPj4448lt50wYQIMBgMyMzO5ZXl5eTAYDKJ/P/zwA7fNkCFD0LhxY9SrVw833XQT1qxZw+2DvYaA8OUlts9169bJPjZ/hj1O9r7li7jbbrtNt8KbLDDaRrGAcdbBPPTQQ3YP+cCBAwXbXLx4EaNGjUJ4eDgaNGiAsWPHovyfeVNY9u/fj6SkJNSpUwfx8fGi7gfC/2AFDN8Cc/nyZVnf1UN6JtvRB6ILiSU8PFzwWckL4MqVK+jSpQveeOMNh9tt2rQJe/bsQZMmTQTLe/XqhdOnTwv+Hn30UbRo0QI333wzAGDXrl3o3LkzsrOzsX//fjz88MN48MEHsWXLFgDCQGTbtq9cuVKw73vvvVf2sfkzfAsMSyC8wOUIGLLA+C+K70y2g3nkkUckTYoDBw7EypUruc98nzIAjBo1CqdPn0ZOTg4qKyvx8MMPY/z48Vi7di0AoKysDP3790e/fv2wYsUK/PLLL3jkkUfQoEEDjB8/XmmTCS8i5kKSO5utHtIza9WqharKyoDKQrLFHQvMoEGDMGjQIIfbFBcXY9KkSfjqq68wePBgwbrg4GDExcVxnysrK7F582ZMmjSJEyb/93//J/jOlClTsHXrVlgsFtx1112SFhgAaNCggWD/eoE9TlbAMAzDuUH1fM86ykJirch6FnBaR/GVkdPBhISESD7kBw8exJdffokffviBGxG99tpr+O9//4vFixejSZMmWLNmDa5fv4733nsPwcHBuPHGG1FQUIAlS5aQgPFzxFxIzZo1g6lpUxyTcCMZDAaYTCZdpGfWqlULVQhsF5KtBUbNF2B1dTVGjx6NGTNm4MYbb3S6/SeffIILFy7g4YcfdrhdaWkpV3TRkQUmLS0Njz76KFq2bIkJEybg4Ycflkwd1xK2Aob/Qg8EAUMWGG3ikRiYvLw8xMbGom3btnj88cdx4cIFbt3u3bvRoEEDTrwANVPWBwUF4bvvvuO2uf322wXpjAMGDMChQ4fw119/if5mRUUFysrKBH+E9xGzwAQFBeGVV14R3V5v6ZlshxjILiS+BSboHzeyWrz88suoVasWJk+eLGv7d999FwMGDIDJZJLcJisrCz/88AMncvgWGP7/58+fj6ysLOTk5GDYsGF44okn8Nprr7l4JP4F29ey9y3/ha5n0S0njVrPx691VBcwAwcOxAcffIDt27fj5Zdfxo4dOzBo0CCuAy8pKUFsbKzgO7Vq1UJkZCRKSkq4bRo1aiTYhv3MbmPLggULEBERwf3Fx8erfWiEDKTSqIcOHYphw4bZba+39EwpARNILiR2Xh2Wr7/+WpUMs3379mHZsmVYtWqVLFFUVFSEr776CmPHjpXcJjc3Fw8//DDeeecdzqIjte9nn30Wt956K7p27YqnnnoKM2fOxKJFi1w7GD8jUC0wjrKQyALj/6guYEaMGIF77rkHnTp1wr333ostW7bghx9+QF5ento/JWDWrFkoLS3l/k6ePOnR3yPEEQviZWnZsiUAIDU1FWvXrkVubq7u0jOdWWD0PpqzWCwYPXo097maYTBw0CA0atTI7Vo/+fn5OHv2LJo1a4ZatWqhVq1aOH78OKZPn46EhAS77VeuXImoqCjcc889ovvbsWMH7r77bixduhQPPvggt5xvdXFEz549UVRUpIspMgJVwJAFRtt4/Mq0bNkS0dHROHLkCPr27Yu4uDicPXtWsE1VVRUuXrzIxc3ExcXhzJkzgm3Yz1KxNSEhIXbBwoT3YTtCsakE2I6+VatWGDlypFfb5S0COQvJYrFwVjbbMeuFCxcwbNgwZGdnuyxYR48ejX79+gmWDRgwAKNHj7aLcWEYBitXrsSDDz4oOoLOy8vDXXfdhZdfftkurk6ugCkoKEDDhg110e8EugvJ0VQCZIHxXzxeB6aoqAgXLlzgMkwSExNx6dIl7Nu3j9vm66+/RnV1NXr27Mlt88033whegjk5OWjbti0aNmzo6SYTbuDIAsMu00OHL0WgupCsVqusuJQpU6Y4dCeVl5ejoKAABQUFAIBjx46hoKAAJ06cQFRUFDp27Cj4q127NuLi4tC2bVvBfr7++mscO3YMjz76qN1v5ObmYvDgwZg8eTKGDRuGkpISlJSU4OLFiwBqXEgHAbzO+86nn36K//3vf/j1119x5MgRvPnmm3jxxRcxadIkAO7XryksLMTYsWPRokULhIaGolWrVsjIyLB7jrKysnDTTTehbt26aN68uSwX1o8//og777wTDRo0QFRUFMaPH29XtoIsMGSB0SKKBYyjDqa8vBwzZszAnj17UFhYiO3bt2PIkCFo3bo1BgwYAABo3749Bg4ciHHjxuH777/Ht99+i4kTJ2LEiBFcTYf7778fwcHBGDt2LA4cOID169dj2bJlmDZtmnpHTngER3Mhscv4wdl6g+3sbCvx6t2FlJ+fj2IHxQpZioqKHBYs3Lt3L7p27YquXbsCAKZNm4auXbtizpw5itrz7rvvolevXmjXrp3duvfffx9Xr17FggUL0LhxY+6PtQwFBQWhAsAF3ndq166NN954A4mJibjpppvw1ltvYcmSJcjIyADgfv2a33//HdXV1Xjrrbdw4MABLF26FCtWrBCkfH/xxRcYNWoUJkyYgF9//RXLly/H0qVL8frrr9v+DMepU6fQr18/tG7dGt999x2+/PJLHDhwAA899JBgO9s6MOy/BpWDsP0NykLSNop7071796J3797cZ1ZUjBkzBm+++Sb279+P999/H5cuXUKTJk3Qv39/PPfcc4JR95o1azBx4kT07dsXQUFBGDZsGF599VVufUREBLZu3Yq0tDR0794d0dHRmDNnDqVQawBHs1GzLqRAtMDo3YWkpAiho22Tk5PtxJ8jCgsLRZezNaXEWLVqFVatWiW53mAwoAuAm3jLBg4caFeQk4+79Wts99+yZUscOnQIb775JhYvXgwA+PDDD3HvvfdiwoQJ3DazZs3Cyy+/jLS0NFGhsWXLFk58sa6xFStWoHPnzjhy5Ag37YNtJd5AKGIHkAVG6yi+Ms46mK+++srpPiIjIx12MADQuXNnzZeWD0TkuJACwQIj5ULSa2eopAihvxcs9ITFQWn9GqCmNk1kZCT3uaKiAnXr1hVsExoaiqKiIhw/flw0kLmiogLBwcGCuJ7Q0FAAwM6dO+0EjK0FRq+Cm4WykLQNzYVEqIqcIF49W2ACNYg3KSnJbgZjMbRQsJAvYHbs2KFKCrjS+jVHjhzBa6+9hscee4xbNmDAAFgsFmzfvh3V1dX4448/uPpKUlatPn36oKSkBIsWLcL169fx119/4emnn7b7jpQLSa/3KwtZYLQNCRhCVcgCE5guJKPRKHADS7Fs2TK/PgcWiwVbt27lPt/Zvz8SEhLcSgFXWr+muLgYAwcOxPDhwzFu3Dhu+bhx4zBx4kTcddddCA4Oxn/+8x+MGDECgHTm1I033oj3338fr7zyCurWrYu4uDi0aNECjRo1Ep0yIVBdSJSFpE1IwBCqQkG8gelCAgCz2Yzs7GxERUXZrYuKinIrhdobWCwWpKSk4JpNXZfi4mKkpKS4LGKU1K85deoUevfujV69euHtt98WrDMYDHj55ZdRXl6O48ePo6SkBLfccguAf2ssiXH//fejpKQExcXFuHDhAubOnYtz584JvkMWGLLAaBESMISqUBBvYFpgWMxmM86cOYMvv/gCtycl4fakJHz5xRc4c+aMX4sXq9WKKVOmiMb3scvS09NdcieNHj0a+/fv57I3CwoK0KRJE8yYMUMQM1hcXIzk5GR0794dK1eulLSqGI1GNG3aFMHBwfjoo4+QmJiImJgYp+1o1KgR6tevj/Xr16NOnTq48847uXW2MTDs/av3+5WykLQNSUtCVciFZC9gGIYJmBcCUHOMffr0AfbsqVnQpw/g58edn5+PoqIiyfUMw+DkyZPIz89HcnKy3fry8nIcOXKE+8yWl4iMjESzZs3srFK29WtY8dK8eXMsXrwY586d47Zli3eeP38eGzduRHJyMq5du4aVK1diw4YN2LFjB7ft999/jwcffBDbt2/nYpJef/119OrVC/Xr10dOTg5mzJiBl156CQ0aNBC0BwhcFxJZYLQJXRlCVSiI1z6Ilz9qp87QP5GbBi61naPyEo5StllycnJw5MgRHDlyxG7iSb5V6P3338eTTz4JhmGQmJiIvLw8zo0EAFevXsWhQ4cEz9/333+PjIwMlJeXo127dnjrrbcE0z0A/96XFy9eRGFhITfY0LvgpiwkbUO9KaEqUhYYq9XKzSR+4MAB9OnTR5edo5gFJlCqmmoZuandUtu5W7/moYcesisuZ0t0dDR2797tcBuxdnzwwQcOv2OxWDB79mwAwKnTp7Hq/ffxv88/B6D/+5UsMNqGYmAIVRGzwBw8eBBt2rThJticOnWq25kd/goJGG2SlJQEk8kkmSVkMBgQHx/v9yngSmEDl9nBBcvZf1xY165d80WzvAZlIWkbEjCEqthaYA4ePIisrCwU2ZSZdzezw18Rm0ogUCbG0zJGoxHLli0DANhKGFbUZGZm6kqAOgpcZrlw4YIqdXD8FSkLDMMwZIHRACRgCFXhCxir1YovvvgCYt2ju5kd/gpZYLSL2WzGxo0b7QrymUwmbNy40a+zqFzBWeAyUPNi13NFdCkBw39+yQLjv5C0JFSF70LauXMnyi5fltzWWWaHFhETMPzOkQSMf2M2mzFk0CCcfPxxXL58Gb0nTMBtOo3XcjdwWQ9IBfHyXeBkgfFf6MoQqsK3wJSUlMj6jp46SEdZSAaDQbK2B+E/GI3GfwvM3XGH36eAu4q7gct6QMoCQ25fbUC9KaEq7Av6/PnzOHPmjKzv6KmDdORC0uMontAuzgKXgRqLqt4Cl/lICRi+BYZcSP4LCRhCNSwWC0aNGgUAKCoqwpMzZiDIQeeox8wORy4kEjCEP8EPXJYiISFB1/etVBYSWWC0AQkYQhXYdEx+BVEAqJbIcNBrZoeYgOHXxMnLy9NV0DKhbdjAZdupCKIiIwFA1hQFWsaZBcZoNMqagJPwDSRgCLeRk45pa4nRa2aHrYDZtGkTbr/9dgA1lYh79+6t2xo4hDYxm834+OOPAQD16tbFQ2PGcDOL62lwIYazGBiyvvg3JGAIt5GTjslaYoKDg5Gbm4tjx47pTrwAwiDegwcPYsSIEXbBzHqtgUNoF3Z6D0EAM/T/AmetKwcPHsSOHTu4gQdrgaH4F/+GBAzhNkqyiOrXr4/k5GTdjuzYIOYzZ85gy5YtAVUDh9AuUpM56vU5BWrc3mwM0LZt23Bn//7IzMzEwYMHyQKjEUjAEG6jJIuoTp06HmyJb7FYLHjvvfcAAIf++ANXrl6V3JZfA4cgfI2t61PvmXNszF5ZWZlg+eXLl5GVlYUvvvgCQM15oLg1/4UEDOE2TueRARD7TzBgaGioF1vmPdgOsby8XNH39FQDh9AuthYY9oWtRwuEo5g95p+/p59+GkCNoKG4Nf+FBAzhNoJ5ZGxEDPvpscceA6BPC4ycIGYp9FQDh9AutunEenYhKYnZY6G4Nf+EBAyhCmw6ZpMmTQTLw8PDkZqaiptvvhmAPgWMnA7RFj3WwCG0i5QFRo8CxhWrJ8Wt+SckYAjVMJvNKCws5Kww69etw5QpU9C+fXv8/fffAPTpQlLaIeq1Bg6hXaRiYPToQnLV6klxa/4HCRhCVWrVqsWJlK5du3JZORUVFQD0aYFR2iHqtQYOoV1YCwyDmhe1nl1IcqZQcATFrfkPJGAI1WFrSly7do1bxv5fjwJGbod45513Ytu2bbqtgUNoF76lxWq16tqFJGcKBUdQ3Jr/QAKGUB1WpLBWF0DfAkZuh5iTk4OHHnoImzdv9kKrCEI+/IJt1dXVunYhAf/G7EVFRQmWB9rcbVqHBAyhOqxI4Vtg9BwDA/zbIUb+M4eMFJTNQPgjfKFSXV2taxcSi9lsxocffggAaNasGXK2bkVKSgrEJAzFrfknJGAI1RETMHq2wLCYzWa88cYbAKRHcpTNQPgjgeRC4lOvXj0ANQOrO+64Ax06dEBqaioaNmgg2I7i1vwTEjCE6ogJGD0H8fJhLUxSs3ADlM1A+B9BQUGc6A4EFxILG6/Hd3e3b98e//d//wcASExM1PXcbVqHBAyhOoHoQmJhO0Q5UDYD4U/wa8EEigVGrK8CgOvXrwMA2rZtq+u527QOCRhCdQItiJePEgFD2QyEP8GvBRMoAkbMAgMETn+ldfRtHyR8QqDGwAD/dohBBoPk1AIGgwEmk4myGQi/grXA8GNg9O5CkrLAsIJG7xZjraPYAvPNN9/g7rvvRpMmTWAwGPDxxx9z6yorK/HUU0+hU6dOqFevHpo0aYIHH3wQp06dEuwjISEBBoNB8PfSSy8Jttm/fz+SkpJQp04dxMfHY+HCha4dIeF1xCwwrAspUASMlCWGshkIf4VvgQmELCRA2FfxBxyBMuDSOooFzJUrV9ClSxcu24LP1atX8eOPP+LZZ5/Fjz/+CIvFgkOHDuGee+6x23b+/Pk4ffo09zdp0iRuXVlZGfr374/mzZtj3759WLRoEebOnYu3335baXMJH+AoiFfvIxpWuBgMBqSmpiKyYUPBespmIPyVQIyBYZ9XvmgDAmfApXUU2wcHDRqEQYMGia6LiIhATk6OYNnrr7+OW265BSdOnECzZs245WFhYYiLixPdz5o1a3D9+nW89957CA4Oxo033oiCggIsWbIE48ePV9pkwsvYVuKtrq5GcXExAKCwsBBWq1W3HSN77FarFe3bt0dGr154Ij0dN998MxYtWoSkpCTdHjuhbcRiYALFhQTU9Fe1ef+3XU/4Hx4P4i0tLYXBYEADm7z6l156CVFRUejatSsWLVokUL+7d+/G7bffjuDgYG7ZgAEDcOjQIfz111+ebjLhJnwLzG+//YZXXnkFP+zdCwB4+eWXkZCQoNtCbuyxs/dzWVkZAKBz586UzUD4NfwYmEBxIfFdvYGYdKB1PCqvr127hqeeegojR45EeHg4t3zy5Mno1q0bIiMjsWvXLsyaNQunT5/GkiVLAAAlJSVo0aKFYF+NGjXi1jW0McsDNTcf/wZkXxyE92Ef+i+//BLWf4QLn6KiIqSkpOjSlcJZYKqrwTAMLl26BAB2Ap4g/A2BBeafZXoXMEFBQahduzYqKysD0uWtdTwmYCorK5GamgqGYfDmm28K1k2bNo37f+fOnREcHIzHHnsMCxYsUJSGymfBggWYN2+eW20m1IEVMD/s3Ys7JbZhGAbp6ekYMmSIrjpJ/v1rtVpJwBCagRUwp06dwuF/YkBcnbFZS4SEhKCyspIsMBrEIy4kVrwcP34cOTk5AuuLGD179kRVVRUKCwsBAHFxcThz5oxgG/azVNzMrFmzUFpayv2dPHnS/QMhXILv+nOEHqvR8gVMZWUlDh8+DAA4d+4cTR1A+DVXrlwBAHy1dSu++ee5XLRokW7dvSyBnDWpdVQXMKx4OXz4MLZt22Y326cYBQUFCAoKQmxsLICa8s3ffPMNKisruW1ycnLQtm1bUfcRUPPiCA8PF/wRvkFJhVm9VaPli7fly5dzL4LXXntN17E/hLaxWCw4fuKE3fLS0lLdTz4qVsyOLDDaQLGAKS8vR0FBAQoKCgAAx44dQ0FBAU6cOIHKykqkpKRg7969WLNmDaxWK0pKSlBSUsKVZt69ezcyMzPx888/488//8SaNWswdepUPPDAA5w4uf/++xEcHIyxY8fiwIEDWL9+PZYtWyZwPRH+i21RKEforRptUFAQjEE1j9Xl8nLBOpqJmvBHrFYrpkyZ4nAbPU8+GsiFN7WOYgGzd+9edO3aFV27dgVQE8/StWtXzJkzB8XFxfjkk09QVFSEm266CY0bN+b+du3aBaBG7a5btw533HEHbrzxRrzwwguYOnWqoMZLREQEtm7dimPHjqF79+6YPn065syZQynUGkGO1Q0AYmJidFeN1mq1orq6WnQdzURN+CP5+fkoKiqSXK/3yUfFXEgUxKsNFAfxJicnS5ZIB+BwHQB069YNe/bscfo7nTt31u0Do2csFgtWrVola9vly5frKoAXqHkZOHoC+C+D5ORkbzWLICSR68bVm7uXxbZuFUAxMFpB31WKCK9isViQkpLiVMQCwIwZM5CSkuKFVnmXQH8ZENpDrhtXb+5elkCefFbr0GzUhCqwfnRn4iU6OhpZWVm6ndsq0F8GhPZISkqCyWSSXG8wGBAfH687dy8LBfFqFxIwhCo486OzrFu3DsOHD/dCi3xDUlISajlwi+n9ZUBoD6PRiGXLljncRs+Tj9oG8VZXV+Pq1asAgJ9++oni1fwYEjCEKsh1iZw9e9bDLfEtRqMRTZs2FV1HM1ET/orZbMZ/eva0Wx4WFob169frrmI2i9Vq5aq2//LLL/jtt9+QmZkJ6z+B+KmpqVT+wI8hAUOoArlO/oWd9iLUxvxMM1ET/gx/sl2Wy5cvY9q0abp8gVssFiQkJODbb78FACx79VVkbdiAssuXBdtR+QP/hQQMoQqsH91R6XE9pk2LwfrU+/Tpwy3btm0bjh07RuKF8FtKSkpEl+vxBc4mHMhxe1P5A/+FBAyhCnw/upSImTRpku5dJ1arlSvJzk5/YTAY0KdPH90fO6FdrFYr9opMvAro7wUuN+GAj95r4WgVEjCEapjNZmzcuJFzobAE/SNoevfu7YtmeQ3WJP3jTz8BAPbu2wegRtwFwqR4hHbJz8/H1X9qn4ihpxe43IQDMaj8gX9BAoZQFbPZLBjJjXnwQURERACAri0QjkzSVVVVujK/E/ojkOoXuXMMgRDDpyVIwBCqU69ePe7/zZo140y1QUH6vN3kmKT1Yn4n9EkgBeG7cgxU/sA/0ecbhfAp/OJPVVVV3NxAerXAyDFJ68X8TuiTpKQk1OcNPGzR0wtcTsIBHyp/4L+QgCFUh83CAWoEDGuZ0OvDH0jmd0KfGI1G9OvXT3Sd3l7gchIO+FD5A/+FBAyhOgaDAcHBwQACwwITSOZ3Qr+0a9dOdLkeX+BswoFt0cmI8HDc+Y+Qa9iwIXJzc6n8gR9DAobwCKwbKRAsMHJM0noxvxP6xGKxYPny5YJlDRs2xLx583T7AjebzSgsLMSkSZMAAMl33IEpU6YgISEBQE0V4uTkZN32W3qABAzhEfgCRu8WGDkmab2Y3wn9wWbQ2VagvXTpEubOnYvNmzf7qGWex2g0okOHDgCABg0aICgoSPcDLj1BAobwCIFkgQGkTdIA0KJFC12OYAnt4yiDTm8F7KRg3d3Xr18HAG7AVatWLZ+1iZAHCRjCIwSSBYaFNUmPHzcOABBcuzYAIDY21pfNIghJnGXQ6amAnRSsgKmsrASAgOmv9AAJGMIjBJoFhoVvkmY7xNr/CBmC8Dcog+7f55O1wARSf6V1SMAQHiEQLTAsbBo5a5RnR3gE4W9QBh25kLQMCRjCIwSqBQYQ1sEByAJD+C/OMuj0VMBOCnIhaRcSMIRHEAiYf5YFSodAAobQCvwMOlv0VsBOCnIhaRcSMIRHYAUMO6oBAqdDsHUZkQuJ8GfYDLo4m1nk9VjATgwpCwy5kPwfEjCER2CtEOyoBggcAUMWGEJrmM1m7Nq1i/vc9oYbZBew++abb3D33XejSZMmMBgM+PjjjwXrGYbBnDlz0LhxY4SGhqJfv344fPiwYJt77rkHzZo1Q506ddC4cWOMHj0ap06d4tZfu3YNDz30EDp16oRatWrh3nvvlX1sn332GXr27InQ0FA0bNjQ7rtSMTCB0l9pGRIwhEcIZAuMrYAhCwyhBSIjI7n/h4aGyn5er1y5gi5duuCNN94QXb9w4UK8+uqrWLFiBb777jvUq1cPAwYMwLVr17htevfujaysLBw6dAjZ2dk4evQoUlJSuPVWqxWhoaGYPHmy5JxNYmRnZ2P06NF4+OGH8fPPP+Pbb7/F/fffL9iGHWCwfRW5kLQD2cgIj0AC5l/IAkNogdDQUO7/fHHhjEGDBmHQoEGi6xiGQWZmJmbPno0hQ4YAAD744AM0atQIH3/8MUaMGAEAmDp1Kved5s2b4+mnn8a9996LyspK1K5dG/Xq1cObb74JAPj2229x6dIlp+2qqqrClClTsGjRIowdO5ZbzpY5YCELjHYhCwzhEVgBQy4kssAQ2oCfiaREwDji2LFjKCkpEVhNIiIi0LNnT+zevVv0OxcvXsSaNWvQq1cvt8T/jz/+iOLiYgQFBaFr165o3LgxBg0ahF9//VWwHaVRaxcSMIRHYDuFs2fPcssCVcCQBYbQGmoJmJKSEgBAI5sA4UaNGnHrWJ566inUq1cPUVFROHHihNtzMP35558AgLlz52L27NnYsmULGjZsiOTkZFy8eJHbjlxI2oUEDKE6FosFb7/9NgCg8PhxbrltcJ9esbW4kIAhtAB/vqPLly97ff6jGTNm4KeffsLWrVthNBrx4IMPis7RJBfWkvLMM89g2LBh6N69O1auXAmDwYANGzZw25ELSbuQgCFUhZ3Z9rLNzLYAkJqaCovF4oNWeRfWfcZCLiTC37FYLGjTpg33+crVq0hISHD7eY2LiwMAnDlzRrD8zJkz3DqW6Oho3HDDDbjzzjuxbt06fP7559izZ4/Lv81WD+bHvISEhKBly5Y4ceIEt8xWwLCiiVxI/g8JGEI1HM1sy6L3mW0BciER2oIddBQVFwuWFxcXIyUlxS0R06JFC8TFxWH79u3csrKyMnz33XdITEyU/B5rBamoqHD5t7t3746QkBAcOnSIW1ZZWYnCwkI0b96cW8Y+n8w/v0sWGO1AAoZQDZrZtgYK4iW0gqNBB7vM2aCjvLwcBQUFKCgoAFATuFtQUIATJ07AYDAgPT0dzz//PD755BP88ssvePDBB9GkSROuHst3332H119/HQUFBTh+/Di+/vprjBw5Eq1atRKInN9++w0FBQW4ePEiSktLBb8JAN9//z3atWuH4n+EWHh4OCZMmICMjAxs3boVhw4dwuOPPw4AGD58OPc9/vNptVpJwGgIspERqkEz29ZAFhhCKygZdCQnJ4tus3fvXvTu3Zv7PG3aNADAmDFjsGrVKsycORNXrlzB+PHjcenSJdx222348ssvOVdr3bp1YbFYkJGRgStXrqBx48YYOHAgZs+eLXiW/vvf/+I4L6aua9euXBsB4OrVqzh06JCgdMOiRYtQq1YtjB49Gn///Td69uyJr7/+Gg0bNuS2sRUwFMSrHUjAEKpBM9vWQBYYQiuoMehITk526DY2GAyYP38+5s+fL7q+U6dO+Prrr522obCw0OF6sXbUrl0bixcvxuLFiyW/xx9g8F1IFAPj/yh2IalRNvrixYsYNWoUwsPD0aBBA4wdOxbl5eWCbfbv34+kpCTUqVMH8fHxWLhwofKjI7wKzWxbQ1BQEIJ450CuBeby5ctIT09H8+bNERoail69euGHH34QbHPw4EHcc889iIiIQL169dCjRw9BQGJJSQlGjx6NuLg41KtXD926dUN2drbD3507dy4MBoPgr127dgqOmNAqNOioeV5Zawu5kLSFYgGjRtnoUaNG4cCBA8jJycGWLVvwzTffYPz48dz6srIy9O/fH82bN8e+ffuwaNEizJ07l0vNJfwTRzPbsuh9ZlsW/uhNroB59NFHkZOTgw8//BC//PIL+vfvj379+nE+/aNHj+K2225Du3btkJeXh/379+PZZ58VZD09+OCDOHToEBdvYDabkZqaip9++snhb9944404ffo097dz504XjprQGjToqIG1kpILSWMwbgCA2bRpE/e5urqaiYuLYxYtWsQtu3TpEhMSEsJ89NFHDMMwzG+//cYAYH744Qdumy+++IIxGAxMcXExwzAMs3z5cqZhw4ZMRUUFt81TTz3FtG3bVnbbSktLGQBMaWmpq4dHuEh2djYTFRXF1AaYjH/+QgwGJjs729dN8w4VFcxLdeowGQBTG2Dee+89p1+5evUqYzQamS1btgiWd+vWjXnmmWcYhmGY++67j3nggQcc7qdevXrMBx98IFgWGRnJvPPOO5LfycjIYLp06eK0jYqoqGCYjIyaP95zrAm03HYXyM7OZgwGAxPMe15rA4zBYGAMAfLchoeHM7UB5vzEiUxecjJTG2Aee+wxXzcrYJH7/lY1C0lO2ejdu3ejQYMGuPnmm7lt+vXrh6CgIHz33XfcNrfffrsgdmDAgAE4dOgQ/vrrL9HfrqioQFlZmeCP8A1msxkrVqwAANT6ZxQTExMja2ZbvaDUAlNVVQWr1WpXQyY0NBQ7d+5EdXU1PvvsM9xwww0YMGAAYmNj0bNnTzsXbq9evbB+/XpcvHgR1dXVWLduHa5duyYZgMly+PBhNGnSBC1btsSoUaMEbilC35jNZmzcuBFNmzYVLDeZTNi4cWNAPLd8Cwy5kLSDqgJGTtnokpISxMbGCtbXqlULkZGRgm3E9sH/DVsWLFiAiIgI7i8+Pt79AyJchu0QmAAtCsXv/A4fPuy09k1YWBgSExPx3HPP4dSpU7BarVi9ejV2796N06dP4+zZsygvL8dLL72EgQMHYuvWrRg6dCjMZjN27NjB7ScrKwuVlZWIiopCSEgIHnvsMWzatAmtW7eW/O2ePXti1apV+PLLL/Hmm2/i2LFjSEpKEi1GSOgTs9mMw4cP46ExYzDMbEbO1q04duxYQIgXgFxIWkU3dWBmzZqF0tJS7u/kyZO+blJAwwXFBeBoJjs7W2ABnD9/vqyqph9++CEYhkHTpk0REhKCV199FSNHjkRQUBA3KhwyZAimTp2Km266CU8//TTuuusuztoFAM8++ywuXbqEbdu2Ye/evZg2bRpSU1Pxyy+/SP7uoEGDMHz4cHTu3BkDBgzA559/jkuXLiErK8vNM0FoCaPRiISEBHTq1Al33HFHQD2zrJWUCtlpC1UFjJyy0XFxcYIJ/oAa8/nFixcF24jtg/8btoSEhCA8PFzwR/gOW4tLoHQGM2fOxMj770e1TTpnUVGR06qmrVq1wo4dO1BeXo6TJ0/i+++/R2VlJVq2bIno6GjUqlVLUBYdANq3b8+5e44ePYrXX38d7733Hvr27YsuXbogIyMDN998s2TQvRgNGjTADTfcgCNHjig4coLQLmIupECzGmsRVQWMnLLRiYmJuHTpEvbt28dt8/XXX6O6uho9e/bktvnmm28EBYlycnLQtm1bQQEiwn+xFSyBIGA2bNiARYsWSa5nGEbWVAr16tVD48aN8ddff+Grr77CkCFDEBwcjB49egjKogPAH3/8wZVFv3r1KoCatFA+RqOR65TlUF5ejqNHj+o6dZYg+JALSZsoFjDulo1u3749Bg4ciHHjxuH777/Ht99+i4kTJ2LEiBFo0qQJAOD+++9HcHAwxo4diwMHDmD9+vVYtmwZV+GR8H8CzQJjtVrxxBNPON3O0VQKX331Fb788kscO3YMOTk56N27N9q1a4eHH34YQM1svevXr8c777yDI0eO4PXXX8enn37K/W67du3QunVrPPbYY/j+++9x9OhRvPLKK8jJyeGePwDo27cvXn/9de7zk08+iR07dqCwsBC7du3C0KFDYTQaMXLkSDfOCEFoB9aFREG82kKxjczdstEAsGbNGkycOBF9+/ZFUFAQhg0bhldffZVbHxERga1btyItLQ3du3dHdHQ05syZI6gVQ/g3gWaByc/Px/nz52VtK1XVtLS0FLNmzUJRUREiIyMxbNgwvPDCC1znOnToUKxYsQILFizA5MmT0bZtW2RnZ+O2224DUNMJf/7553j66adx9913o7y8HK1bt8b777+P//73v9zvHD16VNDWoqIijBw5EhcuXEBMTAxuu+027NmzBzExMa6eDoLQFORC0iaKr5C7ZaMBIDIyEmvXrnX4O507d9b9pH96xvbht3Vr6A0l8ztJuWZSU1ORmprq8LuPPPIIHnnkEcn1bdq0cVp517Yk+7p16xxuTxB6xmq1cu7X06dPkwVGQ+j7rUL4DNuHX+8CRm68SExMjO6rmhKEVrBYLEhISMCvv/4KAMjbsYPL2CMB4//o+61C+IxAi4FhS7I7Y/ny5bo/FwShBSwWC1JSUuxm477+T/LI77//7otmEQogAUN4hECLgWHngZKaUwaoCcJNSUnxYqsIghDDarViypQpDsMhtmzZ4jRjkPAtJGAIjxBoFhjg35LsJpuS7DExMcjKyqIZ1QnCT8jPz7ezvNhy6dIlisP0cyjMmvAIgWaBYTGbzRgyaBBOPv44Ll++jN4TJuC2Pn0C5vgJQgvIDbpXEpxPeB8SMIRHCEQLDAtbkh0AcMcdQAAdO0FoAblB91TM0b8hFxLhEQLVAkMQhP/DBt07illr2LAhZQz6OSRgCI9AAoYgCH+FDboHICliRowYQf2Wn0MChvAIgexCIgjC/2GD7pvaBN0b/6lZdcstt/iiWYQCSMAQHoEsMARB+DtmsxlLly4VTBJs/acS748//uirZhEyIQFDeASywBAE4e9YLBakpqbir7/+slv32muvwWKx+KBVhFxIwBAegSwwBEH4M3KK2aWnp1MxOz+GBAzhEcgCQxCEPyOnmN3JkyepmJ0fQwKG8AiBNpkjQRDagorZaR96qxAegSwwBEH4M1TMTvuQgCE8AsXAEAThz8gpZhcfH0/F7PwYEjCERwgKCgK/WyABQxCEPyGnmF1mZib1XX4MCRjCY/AffOoECILwN6SK2QFARkYGzGazD1pFyIUEDOEx+KKFgngJgvBHzGYzCgsLkXTbbYLlycnJvmkQIRt6qxAegx/ISxYYgiD8FaPRaBesa5uIQPgfJGAIj0EuJIIgtEJwcLDgM/VZ/g8JGMJjkAWGIAitULt2bcFn6rP8HxIwhMcgCwxBEFrBVsCQC8n/IQFDeAyywBAEoRXIhaQ9SMAQHoMsMARBaAVyIWkPEjCEx+B3AEVFRTSrK0EQfoutBYZcSP4PCRjCI1gsFsFMr+vWr0dCQgIsFosPW0UQBCEOWWC0BwkYQnUsFgtSUlJgra4WLC8uLkZKSgqJGIIg/A5bi8sPP/xAVmM/hwQMoSpWqxVTpkwBwzB269hl6enp1DEQBOFXHD58WPB51KhRZDX2c0jAEKqSn58vcB3ZwjAMTp48ifz8fC+2iiAIQhqLxYI1a9faLSersX9DAoZQldOnT6u6HUEQhCdhrcZikNXYvyEBQ6iK7Xwi7m5HEAThSchqrF1IwBCqkpSUBJPJBIPBILreYDAgPj4eSUlJXm4ZQRCEPWQ11i6qC5iEhAQYDAa7v7S0NAA1U5TbrpswYYJgHydOnMDgwYNRt25dxMbGYsaMGaiqqlK7qYQHMBqNWLZsmeg6VtRkZmZSiiJBEH4BWY21i+oC5ocffsDp06e5v5ycHADA8OHDuW3GjRsn2GbhwoXcOqvVisGDB+P69evYtWsX3n//faxatQpz5sxRu6mEhzCbzdi4cSOCbeoqmEwmbNy4EWaz2UctIwiCEMJajaUgq7H/orqAiYmJQVxcHPe3ZcsWtGrVCnfccQe3Td26dQXbhIeHc+u2bt2K3377DatXr8ZNN92EQYMG4bnnnsMbb7yB69evq91cwkOYzWb06tWL+/zySy/h2LFjJF4IgvAryGqsXTwaA3P9+nWsXr0ajzzyiCAmYs2aNYiOjkbHjh0xa9YsXL16lVu3e/dudOrUCY0aNeKWDRgwAGVlZThw4IDkb1VUVKCsrEzwR/gWfmXL//znP9QBEAThl5jNZjz91FN2y8lq7N94VMB8/PHHuHTpEh566CFu2f3334/Vq1cjNzcXs2bNwocffogHHniAW19SUiIQLwC4zyUlJZK/tWDBAkRERHB/8fHx6h4MoRi+YKlbt64PW0IQBOEYvosoJjoaubm5ZDX2czw6W9W7776LQYMGoUmTJtyy8ePHc//v1KkTGjdujL59++Lo0aNo1aqVy781a9YsTJs2jftcVlZGIsbH8Etzh4aG+rAlBEEQjuFP5hgTE4Pk5GTfNYaQhccEzPHjx7Ft2zanFQx79uwJADhy5AhatWqFuLg4fP/994Jtzpw5AwCIi4uT3E9ISAhCQkLcbDWhJmSBIQhCK/Bd3jTg0gYecyGtXLkSsbGxGDx4sMPtCgoKAPybopaYmIhffvkFZ8+e5bbJyclBeHg4OnTo4KnmEh6GBAxBEP4MX8BQf6UNPGKBqa6uxsqVKzFmzBiBG+Ho0aNYu3Yt/vvf/yIqKgr79+/H1KlTcfvtt6Nz584AgP79+6NDhw4YPXo0Fi5ciJKSEsyePRtpaWlkYdEY/Kwx6hAIgvBn+C4kssBoA48ImG3btuHEiRN45JFHBMuDg4Oxbds2ZGZm4sqVK4iPj8ewYcMwe/Zsbhuj0YgtW7bg8ccfR2JiIurVq4cxY8Zg/vz5nmgq4UEqKiq4/5P4JAjCnyELjPbwiAupf//+YBgGN9xwg2B5fHw8duzYgQsXLuDatWs4fPgwFi5cKKgDAwDNmzfH559/jqtXr+LcuXNYvHixwJJDaAO+BUZqagExiouL8cADDyAqKgqhoaHo1KkT9u7dy61/6KGH7Ko5Dxw40G4/n332GXr27InQ0FA0bNgQ9957r8PftVgs6N+/P6KiomAwGDj3JkEQ+odiYLQHqQLCY/AtMHL566+/cOutt6J379744osvEBMTg8OHD6Nhw4aC7QYOHIiVK1dyn20tPNnZ2Rg3bhxefPFF9OnTB1VVVfj1118d/vaVK1dw2223ITU1FePGjVPcdoIgtAsJGO1BAobwGK5UTn755ZcRHx8vECctWrSw2y4kJEQyK62qqgpTpkzBokWLMHbsWG65syDw0aNHAwAKCwsVt5sgCG1DLiTtQbNREx7DFQHzySef4Oabb8bw4cMRGxuLrl274p133rHbLi8vD7GxsWjbti0ef/xxXLhwgVv3448/ori4GEFBQejatSsaN26MQYMGObXAEAQRuPCDeOvUqSP7e85c3nwmTJgAg8GAzMxMbllhYSHGjh2LFi1aIDQ0FK1atUJGRobT/lPOxMh6hywwhMdwRcD8+eefePPNNzFt2jT83//9H3744QdMnjwZwcHBGDNmDIAa95HZbEaLFi1w9OhR/N///R8GDRqE3bt3w2g04s8//wQAzJ07F0uWLEFCQgJeeeUVJCcn448//kBkZKSqx0kQhPbhW2Bq20xEK4VclzcAbNq0CXv27BEUdgWA33//HdXV1XjrrbfQunVr/Prrrxg3bhyuXLmCxYsXO/z9cePGCRJcAs1yRAKG8Bj8GJgdO3bgtj59nM6HVF1djZtvvhkvvvgiAKBr16749ddfsWLFCk7AjBgxgtu+U6dO6Ny5M1q1aoW8vDz07dsX1dXVAIBnnnkGw4YNA1BTl8hkMmHDhg147LHHVD1OgiC0D1+0yE0akevyLi4uxqRJk/DVV1/Z1UYbOHCgIAmhZcuWOHToEN58802nAoadGDlQIRcS4REsFgtOnTrFfb6zf38kJCQ4rczcuHFju1iV9u3b48SJE5LfadmyJaKjo3HkyBFuH4Aw5iUkJAQtW7Z0uB+CIAIXvgtJroCR4/Kurq7G6NGjMWPGDNx4442y9ltaWirLUuxoYuRAgAQMoToWiwUpKSmoZhjB8uLiYqSkpDgUMbfeeisOHTokWPbHH3+gefPmkt8pKirChQsXOOHSvXt3hISECPZTWVmJwsJCh/shCCJw4Zd6OHHiBKxWq9PvsC7vNm3a4KuvvsLjjz+OyZMn4/333+e2efnll1GrVi1MnjxZVjuOHDmC1157zaml2NnEyIEACRhCVaxWK6ZMmQLGRrwA4Jalp6dLdg5Tp07Fnj178OKLL+LIkSNYu3Yt3n77baSlpQEAysvLMWPGDOzZsweFhYXYvn07hgwZgtatW2PAgAEAgPDwcEyYMAEZGRnYunUrDh06hMcffxwAMHz4cO632rVrh02bNnGfL168iIKCAvz2228AgEOHDqGgoMDhLOgEQWgfi8UiqFv27nvvybIYV1dXo1u3bnjxxRfRtWtXjB8/HuPGjcOKFSsAAPv27cOyZcuwatUqWbWwiouLMXDgQAwfPtxpKYfx48djwIAB6NSpE0aNGoUPPvgAmzZtwtGjR2UcsT4gAUOoSn5+PoqKiiTXMwyDkydPIj8/X3R9jx49sGnTJnz00Ufo2LEjnnvuOWRmZmLUqFEAaio179+/H/fccw9uuOEGjB07Ft27d0d+fr6gFsyiRYswYsQIjB49Gj169MDx48fx9ddfC4LrDh06hNLSUu7zJ598gq5du3I+6hEjRqBr165cZ0QQhP5gLcZFxcWC5XIsxs5c3vn5+Th79iyaNWuGWrVqoVatWjh+/DimT5+OhIQEwfdOnTqF3r17o1evXnj77bcVHwd/YuRAgYJ4CVU5ffq029vddddduOuuu0TXhYaG4quvvnK6/9q1a2Px4sUOg+BsrUQPPfQQHnroIaf7JghCHzizGBsMBqSnp2PIkCGiCQjOXN6jR49Gv379BOsHDBiA0aNH4+GHH+aWFRcXo3fv3ujevTtWrlyJoCDltgXbiZEDARIwhKrIfXgC6SEjCMI/UWIxTk5Otls/depU9OrVCy+++CJSU1Px/fff4+233+YsKFFRUYiKihJ8p3bt2oiLi0Pbtm0B1IiX5ORkNG/eHIsXL8a5c+e4bdkMo+LiYvTt2xcffPABbrnlFlkTIwcCJGAIVUlKSoLJZEJxcTEgMqoxGAwwmUxISkryQesIgiD+xV2LMevynjVrFubPn48WLVoIXN5yyMnJwZEjR3DkyBGYTCbBOtYyVFlZiUOHDnFZRnImRg4ESMAQqmI0GrFs2TKkpKTANmSNDWLLzMx0Wg+GIAjC06hhMXbk8hbDdqoSOa7rhIQEgZuLnRg50KEgXkJ1zGYzNm7ciKZNmwqWm0wmbNy4EWaz2UctIwiC+BfWYiyVIWQwGBAfH08WYz+FBAzhEcxmMw4fPoyHxozBMLMZOVu34tixYyReCILwG1iLMQCyGGsQEjCExzAajUhISECnTp1wxx13UCdAEITfQRZj7UIxMARBEERAYzabMWTQIJx8/HFcvnwZvSdMkDV3G+FbSMAQBEEQAQ9rMQYA3HEHQOLF7yEXEkEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmoMEDEEQBEEQmkO3lXjZqcfLysp83JIA5vp1oKKi5v9lZUBwsG/b4y0C9bj5aPkcaLnt7hCox80S6MfvR7DvbfY9LoWBcbaFRikqKkJ8fLyvm0EQBEEQhAucPHkSJpNJcr1uBUx1dTVOnTqFsLAwblp0NSgrK0N8fDxOnjyJ8PBw1fZLuAddF/+Erov/QdfEP6Hr8i8Mw+Dy5cto0qQJgoKkI11060IKCgpyqNzcJTw8POBvMn+Erot/QtfF/6Br4p/QdakhIiLC6TYUxEsQBEEQhOYgAUMQBEEQhOYgAaOQkJAQZGRkICQkxNdNIXjQdfFP6Lr4H3RN/BO6LsrRbRAvQRAEQRD6hSwwBEEQBEFoDhIwBEEQBEFoDhIwBEEQBEFoDhIwBEEQBEFoDhIwCnnjjTeQkJCAOnXqoGfPnvj+++993aSA4ptvvsHdd9+NJk2awGAw4OOPPxasZxgGc+bMQePGjREaGop+/frh8OHDvmlsgLBgwQL06NEDYWFhiI2Nxb333otDhw4Jtrl27RrS0tIQFRWF+vXrY9iwYThz5oyPWhwYvPnmm+jcuTNXGC0xMRFffPEFt56uie956aWXYDAYkJ6ezi2j6yIfEjAKWL9+PaZNm4aMjAz8+OOP6NKlCwYMGICzZ8/6umkBw5UrV9ClSxe88cYbousXLlyIV199FStWrMB3332HevXqYcCAAbh27ZqXWxo47NixA2lpadizZw9ycnJQWVmJ/v3748qVK9w2U6dOxaeffooNGzZgx44dOHXqFMxmsw9brX9MJhNeeukl7Nu3D3v37kWfPn0wZMgQHDhwAABdE1/zww8/4K233kLnzp0Fy+m6KIAhZHPLLbcwaWlp3Ger1co0adKEWbBggQ9bFbgAYDZt2sR9rq6uZuLi4phFixZxyy5dusSEhIQwH330kQ9aGJicPXuWAcDs2LGDYZiaa1C7dm1mw4YN3DYHDx5kADC7d+/2VTMDkoYNGzL/+9//6Jr4mMuXLzNt2rRhcnJymDvuuIOZMmUKwzD0rCiFLDAyuX79Ovbt24d+/fpxy4KCgtCvXz/s3r3bhy0jWI4dO4aSkhLBNYqIiEDPnj3pGnmR0tJSAEBkZCQAYN++faisrBRcl3bt2qFZs2Z0XbyE1WrFunXrcOXKFSQmJtI18TFpaWkYPHiw4PwD9KwoRbeTOarN+fPnYbVa0ahRI8HyRo0a4ffff/dRqwg+JSUlACB6jdh1hGeprq5Geno6br31VnTs2BFAzXUJDg5GgwYNBNvSdfE8v/zyCxITE3Ht2jXUr18fmzZtQocOHVBQUEDXxEesW7cOP/74I3744Qe7dfSsKIMEDEEQqpGWloZff/0VO3fu9HVTCABt27ZFQUEBSktLsXHjRowZMwY7duzwdbMClpMnT2LKlCnIyclBnTp1fN0czUMuJJlER0fDaDTaRYOfOXMGcXFxPmoVwYe9DnSNfMPEiROxZcsW5ObmwmQyccvj4uJw/fp1XLp0SbA9XRfPExwcjNatW6N79+5YsGABunTpgmXLltE18RH79u3D2bNn0a1bN9SqVQu1atXCjh078Oqrr6JWrVpo1KgRXRcFkICRSXBwMLp3747t27dzy6qrq7F9+3YkJib6sGUES4sWLRAXFye4RmVlZfjuu+/oGnkQhmEwceJEbNq0CV9//TVatGghWN+9e3fUrl1bcF0OHTqEEydO0HXxMtXV1aioqKBr4iP69u2LX375BQUFBdzfzTffjFGjRnH/p+siH3IhKWDatGkYM2YMbr75Ztxyyy3IzMzElStX8PDDD/u6aQFDeXk5jhw5wn0+duwYCgoKEBkZiWbNmiE9PR3PP/882rRpgxYtWuDZZ59FkyZNcO+99/qu0TonLS0Na9euxebNmxEWFsb56iMiIhAaGoqIiAiMHTsW06ZNQ2RkJMLDwzFp0iQkJibiP//5j49br19mzZqFQYMGoVmzZrh8+TLWrl2LvLw8fPXVV3RNfERYWBgXG8ZSr149REVFccvpuijA12lQWuO1115jmjVrxgQHBzO33HILs2fPHl83KaDIzc1lANj9jRkzhmGYmlTqZ599lmnUqBETEhLC9O3blzl06JBvG61zxK4HAGblypXcNn///TfzxBNPMA0bNmTq1q3LDB06lDl9+rTvGh0APPLII0zz5s2Z4OBgJiYmhunbty+zdetWbj1dE/+An0bNMHRdlGBgGIbxkXYiCIIgCIJwCYqBIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc5CAIQiCIAhCc/w/8R57z6VofW4AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data = np.loadtxt(\"Fe_alpha_spectrum.txt\")\n", "\n", "# create an arbitrary x axis for this example\n", "x_axis = np.arange(len(data))*0.087\n", "\n", "x, y, indices, n = nx.data.FindPeaks(x_axis,\n", " data,\n", " n=6,\n", " neg=True)\n", "\n", "print(x)\n", "print(y)\n", "\n", "plt.plot(x_axis, data, color = 'black', marker=\"o\")\n", "\n", "for i, j in zip(x,y):\n", " plt.axvline(i, color = 'red', alpha=0.5)\n", " plt.text(i+1, j, f\"{j:.1f}\")\n", "\n", "plt.savefig(\"fig_findpeaks.png\")\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.9" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }