diff --git a/Code in Progress.ipynb b/Code in Progress.ipynb index 076666f..67d3b09 100644 --- a/Code in Progress.ipynb +++ b/Code in Progress.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 2, + "execution_count": 1, "metadata": { "tags": [] }, @@ -5479,7 +5479,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 2, "metadata": {}, "outputs": [ { @@ -5555,7 +5555,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": 3, "metadata": {}, "outputs": [ { @@ -5565,7 +5565,7 @@ " )" ] }, - "execution_count": 4, + "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, @@ -5594,7 +5594,7 @@ }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 4, "metadata": {}, "outputs": [], "source": [ @@ -5609,22 +5609,20 @@ }, { "cell_type": "code", - "execution_count": 67, + "execution_count": 8, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ - "1946\n", - "584\n", "[3.17786029e-01 1.22745332e-32 3.91216727e-01 8.77867643e-26]\n" ] } ], "source": [ - "print(len(Shuffle_Params))\n", - "print(len(Shuffle_Scorr))\n", + "#print(len(Shuffle_Params))\n", + "#print(len(Shuffle_Scorr))\n", "\n", "Scorr = np.genfromtxt(string)\n", "print(Scorr)" @@ -5632,7 +5630,7 @@ }, { "cell_type": "code", - "execution_count": 76, + "execution_count": 9, "metadata": {}, "outputs": [ { @@ -5750,7 +5748,7 @@ "[578 rows x 2 columns]" ] }, - "execution_count": 76, + "execution_count": 9, "metadata": {}, "output_type": "execute_result" } @@ -5797,16 +5795,16 @@ }, { "cell_type": "code", - "execution_count": 77, + "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "" + "" ] }, - "execution_count": 77, + "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, @@ -5837,13 +5835,40 @@ }, { "cell_type": "code", - "execution_count": 129, + "execution_count": 14, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([0.55316554, 0.5341635 , 0.59764036, 0.54541506, 0.53308938,\n", + " 0.51051648, 0.49712333, 0.46845947, 0.44224636, 0.42824524,\n", + " 0.39121673])" + ] + }, + "execution_count": 14, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "array" + ] + }, + { + "cell_type": "code", + "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "df_max = df.groupby('sigma').max()\n", "df_min = df.groupby('sigma').min()\n", "\n", + "# Need to plot 1 , 3 ,5 ... as didn't generate shuffle error bars for 0.1, 0.3, etc.\n", + "# Honestly, these aren't error bars per say, more like \"shuffle ranges\" \n", + "array = np.diag(df_pivot)\n", + "array = array[0:11:2]\n", + "\n", "# Error bar above. \n", "yerr_1 = df_max['FCFC'] - array\n", "# Error bar below\n", @@ -5852,12 +5877,12 @@ }, { "cell_type": "code", - "execution_count": 130, + "execution_count": 31, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnMAAAFLCAYAAAC9cQpTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deZyddX33/9dnZjKTSWYyeybJZJLJZA9bIGENiIARtBaolrr0Vlur0sVWaYuFn0ot2mqLP1Fr74re9aZVW6kUQUVlC4gGCQRZE8hONiBk39eZ7/3HORmSYUJmMsuZM/N6Ph7ncTLfc11XPicXc86b6/oukVJCkiRJ+akg1wVIkiTpxBnmJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLyWFGuC8il2tra1NTUlOsyJEmSjuuJJ57YlFKqa98+qMNcU1MTCxcuzHUZkiRJxxURqztq9zarJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxw5wkSVIeM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh4zzEmSJOWxQb2cl6TcuPm+pXz1gWXH3e7jl0zmmrlT+qAiScpfkVLKdQ05M3v27OTarFLuvfuWXwNw29Xn5rgSSeq/IuKJlNLs9u3eZpUkScpjhjlJkqQ8ZpiTJEnKYw6AkCQ5KEXKY4Y5SRLXzJ1yVEhzUIqUP7zNKkmSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zAEQkiT1U44yVmcY5iRJ6qccZazO8DarJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxw5wkSVIeM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh5zBQgNaC6FI0ka6AxzGtBcCkeSNNB5m1WSJCmPGeYkSZLymGFOkiQpj+U8zEVEY0TcHhHbI2JHRNwREeO6sP/0iPhBRGyKiL0RsSQiPt6bNUuSJPUXOR0AERHDgHnAfuCDQAI+DzwYEaemlHYfZ//Z2f0fAj4MbAcmA2W9WLYkSVK/kevRrB8BmoGpKaXlABHxDLAMuBr48rF2jIgC4N+BB1JKv3PESw/2XrmSJEn9S65vs14OPHo4yAGklFYB84ErjrPvm4EZvEHgkyRJGuhyfWXuJOCuDtoXAVcdZ9/zs89DI+JRYBawFfg+8Dcppb09VqUGhDufXM+Ta7ZxoKWVOV+cx7WXTuXK0xtyXZYkKY/0x8nocx3mqskEsPa2AFXH2XdM9vk24OvAdcBs4EagEfidY+ynQejOJ9dz/R3PcqClFYD12/Zy/R3PAhjoJEmd1h8no8/1bVbIDHpoLzqx3+Hav5tSuiGl9FBK6UvA3wFXRsSMjnaKiI9GxMKIWLhx48YTLFn55qZ7lrD3YMtRbXsPtnDTPUtyVJEkST0j12FuK5mrc+1V0fEVuyNtzj7f16793uzzzI52Sil9M6U0O6U0u66urtOFKr+klFi7ZQ8/fvolPveTxazf1vFd9/Xb9rJm854+rk6SpJ6T69usi8j0m2tvBrC4E/vC66/sHb6q19qNupRntu05wNPrtvPUmm08vW4bT6/dxubdBwAoKSqguLCg7RZre2+66UGm1pczd0Y9c2fUc0pDBQUFnbk4LElS7uU6zP0I+FJENKeUVgJERBMwh0wfuDfyMzLz010G/OSI9kuzzwt7tFL1G/sPtbD4pR08tTYT2p5et51VmzJTEkbApLoyLpo2kpmNlcxsrGTqqHLufuZlrr/j2aNutZYOKeSv3prp93Df4g3874eW8/UHl1M/ooS3TM8Eu3Mn1lBSVJiT9ylJUmfkOsx9C/gYcFdEfJrMVbbPAWuBWw5vFBHjgRXAjSmlGwFSSpsj4gvAZyJiB5nJg2cDNwD/fuR0J8pfra2JVZt3H3XFbfHLOzjYkrkgO7K8hJmNlVw1eywzx1ZyytgKyocOed1xDg9y+OTtz3CgpZWGytKjRrN++IJmtu4+wLwXXuW+xRv44ZPr+d6CNQwvLuTCqXXMnVHPxVPrqRj2+mNLkpRLOQ1zKaXdEXExcDPwHTK3SB8APpFS2nXEpgEU8vo+fjcCO4E/Bf4aeBm4iUwgVB7auHP/EVfcMs879h0CYHhxIaeOreSPzm9mZmMFMxurGFUxtNPHvvL0Bv7rsTVAx6OOqoYX865ZY3nXrLHsO9jCIys2cd/iDdz//Kv89NlXKCwIzmqqbrsd21g9rGfetCRJ3ZDrK3OklNYA7zrONi/SwQjXlFIiM2mwEwfnoT0HDvHsuu3Z0Ladp9ZuaxuoUFgQTBtVzjtOG8PMsZXMHFfJxLoyCvuoL9vQIYVcPK2ei6fV8/etiafXbeO+xRu4b/EGbvzJYm78yWKmjSrnrTPqmTtjFCc3jCDCfnaSpL6X8zCnwaGlNbF0w06eXruNp7KPpRt20podvtJYXcrp4yr5wzlNzGys5KQxFZQW94++agUFwenjqjh9XBWfvGwaqzbt5v5ssPv6g8v52rzljK4Y2tbP7pzmGoqLcj1QXJI0WBjm1ONSSry0fd9Rwe259dvZcyAz+KCidAinNVby1pNGMbOxglPHVlJbVpLjqjtvQu1wPvKmZj7ypmY279rf1s/u9ifW8Z1HV1NeUtTWz+7NU0dSUWo/O0lS7zHMqdu27z3Is+u289TarTy1NnPbdOPO/QAUFxYwY8wIfm92IzMbKzmtsZKmmmED5pZkTVkJV81u5KrZjew72MKvlmX62T3wwgZ+8szLFBUEZzdXM3d6PXNPGkVDZWmuS5YkDTCGOXXJgUOtvPDKjrYrbk+v3caKjbvbXm+uG84Fk2vbpgWZNmrEoLnlOHRIIW+ZUc9bZtTT0pp4au1W7s3ejv3sjxfz2R8vZsboEW0DKE4aYz87SVL3GeZ0TCklVm/e0xbcnlq7jcUv7WibfLe2LDMtyO+c3sDMxipOGVvhLcWswoJg1vhqZo2v5vq3TWfFxl1t/ey+Nm8ZX31gGQ2Vpbxl+kjmzhjF2c3VDCkcHKFXktSzDHNqs3nXfp5et42nsiNLn167je17DwKZCXZPGVvBH2QHKJzWWMmYiqFeWeqkiXVlTLywjKsvnMimXfuZ9/yr3Lt4A99/fC3//uvVlA8t4qKpI7P97Oo6nCtPkqSOGOYGqb0HWlj00vbXbpeu28baLZlpQQoCptSX87aTR7UFt8kjyyjyylGPqC0r4ffObOT3zmxk74EWfrlsY7af3av86OmXGFIYnNNcw1tn1HPJ9HrG2M9OkvQGDHODQEtrYsXGXUf1c3vhlZ20ZOcFaags5bTGCt5/znhOy66iMKzY/zT6QmlxIW89aRRvPWkULa2J36zZ2jaf3WfuWsRn7lrEyQ0jmDt9FHNn1DN9dLlXQyVJR/EbewB6Zfu+o4Lbs+u3s2t/ZhWF8qFFzGys5E8unMhpjZWc1ljByPLOr6Kg3lNYEJzZVM2ZTdVc/7ZprNi4q20AxVceWMrN9y9lbFUpb5lez1tn1HPmBPvZSZIMc3lv576DPLv+tT5uT63dxoYdmWlBhhQGM0aP4J1nNHBadhWFCTXDKeijVRR04iKCSSPLmTSynD998yRe3bmPB55/lfsXb+A/H1vDrY+8SEXpEC6aWsfcGaO4cGodZSX+OkvSYOSnfx452NLKkld2HhXclm/cRcquojChdjjnNtdwWnZakOmjRzB0SP9YRUHdM7J8KO89axzvPWscew4c4uGlmfns5r2wgTufeoniwgLOnVjTNu1J/QivtkrSYGGY66dSSqzbupcnjwhuz63fzv5DmWlBqocXM7Oxkt8+bUzmdunYCiqHFee4avWFYcVFXHbyKC47eRSHWlp5YnW2n93zG/j0nc/x6Tuf47SxFdlgN4op9WX2s5OkAcww109s23Mge8Uts5LC0+u2s2X3AQBKigo4pSE7QCF71W1sValf0KKosICzm2s4u7mGT/3WdJa9uov7Fm/g3sUb+NK9S/nSvUsZVz2sbd3YM5uqHJUsSQOMYa6X3Pnkem66ZwkvbdvLmMpSrr10Klee3gDAvoMtLH55R9sVt6fXbuPFzXsAiIDJI8t4y/SRbcFtSn25Hd11XBHBlPpyptSX82cXTWLDjn3c/3xmAMV3H13Nt+evonLYEC7Ozmf3pil1DLefnSTlPT/Je8GdT67n+jueZe/BzMLy67ft5drbn+YHC9eyc/8hnn95BwdbMh3dRlcM5bSxlbz7zHHMbMxMC2JHdvWE+hFD+f2zx/P7Z49n1/5DPLz0tfns7nhyPcVFBcyZWMPcGaN4y/SRjLSfnSTlJVNDL7jpniVtQe6wgy2JR1Zs5tyJNXz4gua2tUvtqK6+UFZSxNtPGc3bTxnNoZZWHn/xcD+7V3jwh8/y//0QZjZWMndGZtqTSSPtZydJ+cIw14PefcuvWbBqyzFfT8AjKzbT0pr4m8um9V1h0hGKsiNfz51Yw2feMZ0lG3Zy36LMAIqb7lnCTfcsoalmWNsAilnjqyh0OhtJ6rcMcz3otqvPBWDOF+exftve173eUFnK/Osu7uuypGOKCKaNGsG0USP480sm8/L2vdz//Kvct3gDtz7yIt/65Sqqhg3h4mn12X52ta4OIkn9jJ/KveDaS6ce1WcOMgvVX3vp1BxWJR3f6IpS3n/OeN5/znh27jvIL7L97O5b/Ar/85t1lBQVcP6kWuZm142tKy/JdcmSNOgZ5nrB4VGrxxrNKuWD8qFDeMepY3jHqWM42NLKY6u2tK0b+8ALrxLxLKc3VjJ3Rmbd2Ekjy3JdsiQNSoa5XnLl6Q2GNw0YQwoLmDOpljmTavnb357B8y/vbBtA8Y8/f4F//PkLNNcOb1uB4vRx9rOTpL5imJPUJRHBjDEjmDFmBB9/y2Re2ra3bT67f/vVKm55eCU1w4u5ZPpI5s4YxfmTaiktdlk5SeothjlJ3TKmspQPnNvEB85tYse+gzy0JNPP7mfPvsJ/L1zH0CEFnD+pjrfOqOfi6SOpLbOfnST1JMOcpB4zYugQLj9tDJefNoYDh1pZsGoz9y3ewP2LN3D/8xuIgFnjqtpuxzbXlXHnk+t5cs02DrS0MueL8+xfKkldZJiT1CuKiwq4YHIdF0yu4+8uP4lFL+3g3uwAii/87AW+8LMXGFlezJbdBznUmlkRZf22vVx/x7MABjpJ6iTDnKReFxGc3FDByQ0V/OXcKazdsocHns+EusNB7rC9B1u44a7nGFFaxJT6choqS12NQpLegGFOUp9rrB7GH8yZwN/9eHGHr+/Yd4gP3boQyCxFNrm+jKn15UypL2fqqMxzbVmxIU+SMMxJyqExlaUdrpYyumIo//ze01myYSdLX9nJkg07uWfRK3z/8bVt21QPL2bK4ZA3qpyp9eVMri+nonRIX74FSco5w5yknDnWail/c9k0ZjdVM7upuq09pcSmXQdYumEnS17ZmXnesJPbn1jH7gOv7T+6YuhRV/Cm1pczaWSZ06NIGrAMc5Jy5vAgh0/e/gwHWlppeIPVUiKCuvIS6spLmDOptq09pcT6bXuzIW9XW9j79crNHDjUmt0XxlcPOzrkjSpnQu1whhQW9M2blaReYpiTlFNXnt7Afz22BoDbrj63y/tHBGOrhjG2ahgXT6tvaz/U0srqLXtYdmTI27CTB154lZbsoIshhUFzbVn2Nm1ZW8hrrBpGgStYSMoThjlJA1JRYQET68qYWFfGZSe/1r7/UAsrN+4+6nbtU2u38uOnX2rbpnRIIZMPh7sj+uTVjyhx0IWkfscwJ2lQKSkqZProEUwfPeKo9t37D7Hs1V1tAy6WbtjJw0s3cvsT69q2GTG06KjbtFOyI2yrhxf39duQpDaGOUkChpcUMbOxkpmNlUe1b92dGXRx+Dbt0ld28eOnX+J7Cw61bVNXXsKUdlfyptSXU1biR6yk3ucnjSS9garhxZzdXMPZzTVtbSklXt25/7VRtdnn7z+29qiRuQ2VpUdcycuEvYl1ZQwd4shaST3HMCdJXRQR1I8YSv2IobxpSl1be2trYt3WvW23aQ+HvF8u28jBlsygi4KAptrhr5sEualmGEWOrJV0AnIe5iKiEbgZmAsEcD/wiZTSmk7sm47x0ukppad6rkpJOr6CgmBczTDG1Qxj7ozXRtYebGnlxU27j5oE+YVXdvLzRa+Qsp9ixYUFTBxZlhlVmx1wcXg5M0fWSnojOQ1zETEMmAfsBz4IJODzwIMRcWpKaXcnDnMrcEu7tqU9WackdceQwgImZ1eo4NTX2vceaGHFxl1HTYL82Kot3PnUayNrhxcXMrndqNopo8qoK3NkraSMXF+Z+wjQDExNKS0HiIhngGXA1cCXO3GM9SmlR3uvREnqHaXFhZzcUMHJDRVHte/Yd/Do+fFe2cl9z2/gtoWvLWdWNWzI6yZBnjKynIphLmcmDTa5DnOXA48eDnIAKaVVETEfuILOhTlJGlBGDB3CrPHVzBpffVT7pl37j5o6ZckrO7njN+vZtf+1kbWjRgx93STIk0aWMaw41x/3knpLrn+7TwLu6qB9EXBVJ4/xJxFxLdACPAr8bUrplz1UnyT1G7VlJdROKuG8dsuZvbR932shL/v8Hys3s/+I5czGHV7O7IjbtRNqh1Nc9PpBF3c+uZ4n12zjQEsrc74475hLrEnqH3Id5qqBrR20bwGqOrH/d4GfAC8B44FrgXkRMTel9FBPFSlJ/VVE0FBZSkNlKRdNG9nW3tKaWLNlz1H98Za+spN5RyxnVlQQTKgdftSAi7Vb9vDl+5ZwoCUTBNdv28v1dzwLYKCT+qlchznIDHpor1O9elNK7z/ix19GxF3Ac2QGUZzf0T4R8VHgowDjxo3rWqWSlCcKs0FtQu1wLjt5VFv7/kMtrNq0+4g58nbx7Lrt3P3My8c81t6DLdx0zxLDnNRP5TrMbSVzda69Kjq+YveGUko7I+Ju4I/eYJtvAt8EmD179rGmNpGkAamkqJBpo0YwbdTrlzNb/uourviX+R3ut37bXm646znOm1jLuc01DrSQ+pFch7lFZPrNtTcDWHyCxww6vtonSTqG4SVFnNZYSUNlKeu37X3d6yVFBfxg4Tr+49erKQg4paGC8ybVMmdiLbObqlzVQsqhXIe5HwFfiojmlNJKgIhoAuYA13X1YBExAvgtYEEP1ihJg8a1l07l+juePWpZstIhhXzhnafw9lNG8+SarcxfsZlHlm/iWw+v5F8fWkFxUQGzxlUxZ1IN502q5dSGClezkPpQrsPct4CPAXdFxKfJXFH7HLCWIyYCjojxwArgxpTSjdm2vwamAg/y2gCIvwZGAb/fh+9BkgaMw/3iPnn7MxxoaaWhsvSo0ayH16n9y7lT2LX/EI+v2sL85ZuYv2IzX7p3Kdy7lPKSIs5urmbOpFrmTKpl8sgyJziWelFOw1xKaXdEXExmOa/vkLlF+gCZ5bx2HbFpAIXAkf+rtwT4neyjAtgBzAf+KKX0WB+UL0kD0pWnN/Bfj2VWVLzt6nOPuV1ZSREXTRvZNop28679/HrlZuYv38z85Zu4//lXAagrL+G8iTXMmVjLeZNqGFs1rPffhDSI5PrKHNk1WN91nG1epN0I15TSj4Ef915lkqSuqCkr4R2njuEdp44BYO2WPTyyYlM23G3mruwyZeNrhnHexNrMbdmJtVQPL85l2VLey3mYkyQNTI3Vw3h39TjefeY4Ukos3bCL+cs38ciKTfz46Zfarv7NGD2irb/dWU3VDC/xq0nqCn9jJEm9LiKYOiqzvNiHzp/AoZZWnlm/nUeWZ67c/fsjq/nWL1dRVBCcPq4ye+WulpmNlR2uUiHpNYY5SVKfKyos4IxxVZwxroqPXTyZvQdaWLh6C/OXb+aRFZv42rxlfPWBZQwrLuTMpmrOn5Tpbzd91AgKChxMIR3JMCdJyrnS4kIumFzHBZPrANi+5yC/Xrk52+duE3//0+cBqB5ezLnNNZw3KTOgYnzNMEfKatAzzEmS+p2KYUO47ORRbUuRvbJ9X3YKlE08snwzdz+bWX6sobI0M1I2e+VuZPnQXJYt5YRhTpLU742qGMq7Zo3lXbPGklJi5abdbf3t7l28gR88sQ6AKfVlbf3tzm6uZsRQlx3TwGeYkyTllYhgYl0ZE+vKeP+5TbS0Jha/tIP52Vuy3398Dbc+8iIFAaeOrWRO9pbsGeNddkwDk2FOkpTXCguCU8ZWcMrYCv74wonsP9TCb1Zva+tv941frORfHlxBSVEBs5uqMitTTKzl5IYKCh1MoQHAMCdJGlBKigo5d2IN506s4a/eOpWd+w7y2KrXRsr+08+XAEsYMbSIc5prssuO1TCxzmXHlJ8Mc5KkAa186BAumV7PJdPrAdi4cz+PZAdSzF+xiXsXbwCgfkQJ502sbRtQMaayNJdlS51mmJMkDSp15SVcMbOBK2Y2ALBm8562/nYPL93ID59cD0Bz7fC2KVDOnVhD5TCXHVP/ZJiTJA1q42qGMa5mHO89axytrYklG3Zmlx3bzA9/s57vPrqGCDhpzAjmTKxtW3astNjBFOofDHOSJGUVFATTR49g+ugRfPiCZg62tPL02m3Mz96S/fb8Vdzy8EqGFAanj6vi/Gx/u1PHVjKk0GXHlBuGOUmSjmFIYQGzm6qZ3VTNx98ymT0HDvH4i1szc9yt2MTN9y/ly/fB8OJCzm6uaetvN7W+3GXH1Gc6FeYiogT4JbATuCyldPAY2xUDPwOGAxccaztJkvLRsOIiLpxSx4VTMsuObd19gEdXbs72udvMvBdeBaBmeDHnZoPd+ZNqaawelsuyNcB19src7wOzgN9+o4CWUjoQETcBP83uc2u3K5QkqZ+qGl7M204ZzdtOGQ3AS9v2tvW3m798Ez95JrPsWGN1aVt/u/Mm1lBbVpLLsjXAdDbMvRNYmVL66fE2TCn9PCKWAVdhmJMkDSJjKku5anYjV81uJKXEio27Mv3tlm/i7mdf5vuPrwVg2qjy7LJjNZzdXENZib2edOI6+1/P6WSutnXWw8Dbu16OJEkDQ0QwaWQ5k0aW88HzMsuOPbd+O79avolHVmziewtW8+35qygsCGY2VjJnYg3nTarl9HGVlBQ5Ulad19kwVwts6MJxNwA1XS9HkqSBqbAgOK2xktMaK/mziyax72ALv1m9ta2/3dcfXM7X5i1n6JACzmyqblt2bMaYERQWBHc+uZ4n12zjQEsrc744j2svncqVpzfk+m2pH+hsmNsLlHXhuGXAvq6XI0nS4DB0SGGmD92kWq69FHbsO8iClVuyfe428cWfvQBARekQxleX8vwrOznYkgBYv20v19/xLICBTp0Oc2uBM7tw3NnAmq6XI0nS4DRi6BDmzqhn7ozMsmOv7tjXNpDijt+spyWlo7bfe7CFf/jp81wxc4xryg5ynZ3h8CHgnIiYfbwNI2IWcB7wYDfqkiRpUBs5YihXnt7ATVedRmu7IHfYqzv3M+vz93P1dxby7V+t4rn122lp7XhbDVydvTL3deBPgB9ExNtTSs93tFFETAN+ALQA/7tnSpQkaXAbU1nK+m17X9deOWwIF08byYJVm7lnUaZre/nQIs5squasCdWcPaGakxsqXJ1igOtUmEspLYmIG4HPAk9GxO3APGAdkICxwCXAu4AS4IaU0pJeqViSpEHm2kuncv0dz7L3YEtbW+mQQj772ye19Zl7adteHn9xC4+u3MJjq16bwHhYcSGzxldxVjbgndZYydAhjpYdSDo9sU1K6caIOAT8LfA+4L3tNgngIPCplNIXeq5ESZIGt8OB7ZO3P8OBllYaKktfN5p1TGUpV8xs4IqZmbaNO/fz+ItbeGzVFh5duZkv37+UlKC4qICZjZWcM6GasybUcMb4SoYVO89dPuvS2Usp/UNEfA/4EDAHGE0mxL0E/Ar4vyml1T1epSRJg9yVpzfwX49lxhbedvW5x92+rryEt58ymrdnV6fYtucAj7+4lcdWbWbBqi18/cHltM5bTlFBcMrYCs6aUM05E2qY1VTFiKFDevW9qGd1OYpnw9rf9kItkiSpl1QOKz5qtOyu/Yd4YvVWFqzczGOrtvDtX63ill+spCBg+ugRnD2hhrMmZG7NVg8vznH1eiNeV5UkaRAqKyniwil1XDilDiAzifGarTy2agsLVm5pW6ECYEp9WXZARQ1nT6hm5IihuSxd7XQqzEXESuArKaWvHdE2DmhKKT3cW8VJkqS+MXRIIedNrOW8ibUAHDjUyrPrt2UHVGzhh79Zz3cfzdzmnVA7vG1AxdnN1YytGpbL0ge9zl6ZawIq27X9IXAD4JAYSZIGmOKiAmaNr2bW+Gr+7CI41NLK4pd3sGDlFhas2sLPF73CbQvXAtBQWcrZ2VuyZ02oZkLtcCcy7kPeZpUkScdVVFjAqWMrOXVsJR95UzOtrYklG3Zmbsuu2szDyzZyx5Prgczgi8Pz3J09oYbJI8soKDDc9RbDnCRJ6rKCgmD66BFMHz2CD57XREqJlZt2syA7z92CVVu4+5mXAagaNqRtIuNzmmuYPnoEhYa7HmOYkyRJ3RYRTKwrY2JdGe87exwpJdZt3cuCVVsyI2Zf3MK9i7OrVJQUMaupqm3E7KljXaWiOwxzkiSpx0UEjdXDaKwexu/OGgvAK9v3sWDV5uyt2S08tOQFILOaxRnjKzmrqYazm6uZ6SoVXdKVMPfmdp0Z3wwQEZ8hM3Fweyml9LkTL02SJA0koyqGHrVKxaZd+3k8G+weW7WFrzywlHQ/FBdmVqk4PKBi1vgqhpd4/elYuhTmso/2/q7dz4lMuEvAccNcRDQCNwNzs/vdD3wipbSmC7UREdcD/wDMTymd35V9JUlS36stK+Ftp4zmbdlVKrbvOcjC1Zlwt2DVFv71Fyv4+oPLKSwITm6oyC5BVs3spmoqSl2l4rDOhrn2ga1HRMQwYB6wH/ggmQD4eeDBiDg1pbS7k8dpBj4FvNobdUqSpN5XMWwIl0yv55Lpr61S8ZvVW9tGzP7f+S9yy8MriYDpo0ZkB1RUc2ZTNTVlJTmuPnc6FeZSSr0S5oCPAM3A1JTScoCIeAZYBlwNfLmTx/lX4HvAVOwHKEnSgFBWUsSbptTxpiNWqXhyzTYeW7WFx+CYh1gAABQ1SURBVF7czPcfX8Otj7wIwKSRZW1z3Z3TXEP9IFqlItfB53Lg0cNBDiCltCoi5gNX0IkwFxHvA84A3gvc0VuFSpKk3Bo6pJBzJ9Zw7sQaYHJ2lYrtbYMq7nrqJb63INNLa3zNsGy4yyxBNraqdMBOZNzpMBcR7wCqgO+nlA4eY5ti4N3AlpTS3Z047EnAXR20LwKu6kRNVWT6230ypbRloJ4kSZL0eplVKqqYNb6KP31zZpWK51/e2Rbu7l28gf9euA6AMRVDswMqMiNmmwfQKhWdXZv1FOBO4F9TSt851nYppQMRMRv404iYmVJadJxDVwNbO2jfQiY4Hs9NwFLg1k5sK0mSBrCiwgJOGVvBKWMr+PAFmVUqlr26iwXZSYx/tXwzdz71EpAZfHHkEmRT68vzdpWKzl6Z+xBwgM4NhPgc8OHs45pObJ86aDvuv2ZEXAB8ADgjpdTRMY6130eBjwKMGzeus7tJkqQ8U1AQTB1VztRR5Xzg3MwqFas27W6b527Bys3c/WxmlYqK0swqFec0Z8LdjNEjKMqTiYw7G+beDDyUUtp0vA1TSpsi4iHg4k4cdyuZq3PtVdHxFbsj3QL8G7AuIiqzbUVAYfbnvSml/R3U903gmwCzZ8/udAiUJEn5LSJoriujua6M95yVuaCzdsuezICK7IjZ+5/PrFJRVlLErPFVbSNmT2mopLjo6HB355PreXLNNg60tDLni/O49tKpXHl6Q5+/r86GuYnAA1047vNAZ+Z6W0Sm31x7M4DFx9l3evbxxx28tpXMVcGvdKIGSZI0SB1epeJd2VUqNuzYl53EeDMLVm7hpnuWADB0SAFnjKtquy27futebrhrEQdaWgFYv20v19/xLECfB7rOhrliMrdZO+sA0JkJX34EfCkimlNKKwEiogmYA1x3nH0v6qDtK0Ah8OfA8g5elyRJOqb6EUO5/LQxXH7aGAA279rP4y9ubRtU8dUHlnGszl17D7Zw0z1L+m2Y2wJ0pYPZOGBzJ7b7FvAx4K6I+DSvrRqxlsxtVAAiYjywArgxpXQjQErpofYHi4htQFFHr0mSJHVVTVkJl508istOHgXA9r0HeWL1Fj5068IOt39p296+LA+Azvbs+w0wNyKOe7UtIoaSWZrryeNtm13h4WIyI1K/Q2bi31XAxSmlXUcelswVt/zoiShJkgakitIhXDytnobK0g5fH3OM9t7U2XD0Q6AO+PtObHsjUAv8T2cOnFJak1J6V0ppREqpPKV0ZUrpxXbbvJhSipTSZ49zrDe7LqskSept1146ldIhhUe1lQ4p5NpLp/Z5LZ0Nc/8BvABcExHfiYjJ7TeIiMkR8R/AX5EZAHHM+egkSZLy2ZWnN/CFd55CcXb6kobKUr7wzlP672jWlNLBiLiSzIjW3wfeFxHrgXVk+rmNzT4i23ZlSulQ75QsSZKUe1ee3sB/PZZZPuy2q8/NWR2d7oOWUloKzCIzt9sBMuHtHOBcoDHb9n+A2UeutSpJkqTe0+m1WQFSSq8CH4mIPwdmA6PJXI17CViYUtrX8yVKkiTpWLoU5g7LhrZf9XAtkiRJ6qJO3WaNiA9ExKm9XYwkSZK6prN95m4FrjyyISI+GBHzerwiSZIkdVp3JuFtAi7soTokSZJ0AlxRQZIkKY8Z5iRJkvKYYU6SJCmPdSXMpV6rQpIkSSekK/PMfTYiPtu+MSJajrF9Simd0Dx2kiRJ6pyuhK3o4rG7ur0kSZK6qFNhLqVk3zpJkqR+yJAmSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUx4pyXYDUm26+bylffWDZ69qbrrv7qJ8/fslkrpk7pa/KkiSpxxjmNKBdM3eKIU2SNKB5m1WSJCmPGeYkSZLymGFOkiQpjxnmJEmS8phhTpIkKY8Z5iRJkvJYzsNcRDRGxO0RsT0idkTEHRExrhP7jY+IuyJidUTsjYhNEfFQRLytL+qWJEnqD3Ia5iJiGDAPmAZ8EHg/MBl4MCKGH2f3MmAT8Gng7cAfAbuAn0bEO3utaEmSpH4k15MGfwRoBqamlJYDRMQzwDLgauDLx9oxpbSITIBrExF3A6uAPwTu6KWaJUmS+o1c32a9HHj0cJADSCmtAuYDV3T1YCmlQ8B24GCPVShJktSP5TrMnQQ810H7ImBGZw4QEQURURQRoyLiM8AU4F96sEZJkqR+K9e3WauBrR20bwGqOnmMfwL+KvvnXcB7UkoP9EBtkiRJ/V6ur8wBpA7aogv7fwU4E/ht4GfAf0bEO461cUR8NCIWRsTCjRs3dq1SSZKkfibXYW4rmatz7VXR8RW710kprUspLUwp/SSl9HvAo8CX3mD7b6aUZqeUZtfV1Z1Q0ZIkSf1Frm+zLiLTb669GcDiEzzmQuATJ1yRpF53831L+eoDy17X3nTd3Uf9/PFLJnPN3Cl9VZYk5aVch7kfAV+KiOaU0kqAiGgC5gDXdfVgEVEAnA+s6MEaJfWwa+ZOMaRJUg/JdZj7FvAx4K6I+DSZ/nOfA9YCtxzeKCLGkwloN6aUbsy2fZbMLdr5wCvAKDLzzp0FvK/v3oIkSVLu5DTMpZR2R8TFwM3Ad8gMfHgA+ERKadcRmwZQyNF9/H5D5nbqe4AKMoHuaeCClNL8PihfkiQp53J9ZY6U0hrgXcfZ5kXajXBNKf2IzG1aSZKkQSvXo1klSZLUDYY5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxWlOsCJEm5d/N9S/nqA8te19503d1H/fzxSyZzzdwpfVWWpE4wzEmSuGbuFEOalKe8zSpJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMdyHuYiojEibo+I7RGxIyLuiIhxndhvdkR8MyJeiIg9EbEmIr4XERP6om5JkqT+IKdhLiKGAfOAacAHgfcDk4EHI2L4cXZ/D3AS8DXgbcB1wBnAwoho7LWiJUmS+pGiHP/9HwGagakppeUAEfEMsAy4GvjyG+z7jymljUc2RMR8YFX2uDf0SsWSJEn9SK5vs14OPHo4yAGklFYB84Er3mjH9kEu27Ya2Ag09HCdkiRJ/VKuw9xJwHMdtC8CZnT1YBExHRgJPN/NuiRJkvJCrsNcNbC1g/YtQFVXDhQRRcA3yFyZ+7fulyZJktT/5TrMAaQO2uIEjvN14Dzgf6WUOgqImQNHfDQiFkbEwo0bX3enVpIkKa/kOsxtJXN1rr0qOr5i16GI+ALwUeBDKaV732jblNI3U0qzU0qz6+rqulSsJElSf5Pr0ayLyPSba28GsLgzB4iIT5GZluQvUkrf6cHaJEmS+r1cX5n7EXBORDQfboiIJmBO9rU3FBF/AXwe+FRK6Z97qUZJkqR+K9dh7lvAi8BdEXFFRFwO3AWsBW45vFFEjI+IQxFxwxFt7wG+AvwcmBcR5xzx6PJIWEmSpHyU09usKaXdEXExcDPwHTIDHx4APpFS2nXEpgEUcnT4vCzbfln2caRfAG/upbIlSZL6jVz3mSOltAZ413G2eZF2I1xTSn8A/EFv1SVJkpQPcn2bVZIkSd1gmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPJbz0aySJKljN9+3lK8+sOx17U3X3X3Uzx+/ZDLXzJ3SV2WpnzHMSZLUT10zd4ohTcflbVZJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjTk0iSZLUSf1x7r9IKfXJX9QfzZ49Oy1cuDDXZUiSJB1XRDyRUprdvt3brJIkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJeWxQr80aERuB1b3819QCm3r571DXeE76J89L/+M56Z88L/1PX52T8SmluvaNgzrM9YWIWNjRorjKHc9J/+R56X88J/2T56X/yfU58TarJElSHjPMSZIk5THDXO/7Zq4L0Ot4Tvonz0v/4znpnzwv/U9Oz4l95iRJkvKYV+YkSZLymGHuBEREY0TcHhHbI2JHRNwREeM6ue/QiLgpIl6OiL0R8euIeFNv1zwYnOh5iYjZEfHNiHghIvZExJqI+F5ETOiLugey7vyutDvO9RGRIuJXvVHnYNPd8xIR0yPiBxGxKfs5tiQiPt6bNQ903fxeGRcR/5797NoTEUsj4vMRMby36x7IImJsRPxz9nt6T/YzqKmT+/bpd71hrosiYhgwD5gGfBB4PzAZeLCTvzj/BnwEuAF4B/AycE9EzOydigeHbp6X9wAnAV8D3gZcB5wBLIyIxl4reoDrgd+Vw8dpBj4FvNobdQ423T0vETEbWACUAB8G3g78/0Bhb9U80HXnnGRfvx94E/AZ4LeA/wP8FfDtXix7MJgE/B6wFfhlF/ft2+/6lJKPLjyAjwMtwKQj2iYAh4C/PM6+pwEJ+MMj2oqAJcCPcv3e8vnRzfNS10HbeKAVuDHX7y1fH905J+2Ocw9wC/AQ8Ktcv698f3Tzd6UAWAT8MNfvYyA9unlO3pr9Xnlru/YvZvcfluv3l68PoOCIP384++/c1In9+vy73itzXXc58GhKafnhhpTSKmA+cEUn9j0I3HbEvoeA7wOXRkRJz5c7aJzweUkpbeygbTWwEWjo4ToHk+78rgAQEe8jc5X0+l6pcHDqznl5MzAD+HKvVTc4deecFGefd7Rr30YmfEdPFTnYpJRaT3DXPv+uN8x13UnAcx20LyLzIXe8fVellPZ0sG8xmUu6OjHdOS+vExHTgZHA892sazDr1jmJiCrgZuCTKaUtPVzbYNad83J+9nloRDwaEQcj4tWI+FpElPZolYNLd87J/cAy4B8jYkZElEXExWSu9n0jpbS7Z0tVJ/T5d71hruuqydw/b28LUNWNfQ+/rhPTnfNylIgoAr5B5srcv3W/tEGru+fkJmApcGsP1qTunZcx2efbgHuBucA/kbkF9Z89VeAgdMLnJKW0j0zIPnwLfCfwAPAT4GM9W6Y6qc+/64t6+oCDREeT83XmUnZ0Y18dX0/9234dOA/4rZRSR7+Q6rwTOicRcQHwAeCMlO1woh51or8rhy8AfDeldEP2zw9FRCHwxYiYkVJa3CMVDj4n+rsylEy4Hklm4MQa4CwyHe8PAX/SgzWqc/r8u94w13Vb6ThVV9FxEj/SFqCjoeZVR7yuE9Od89ImIr4AfBT4YErp3h6qbbDqzjm5hcxV0XURUZltKwIKsz/vTSnt77FKB5funJfN2ef72rXfS6bD/UzAMNd13Tknf0SmL+OklNKKbNvDEbEd+GZEfCOl9HSPVarO6PPvem+zdt0iMvfD25vB8T/EFgETssPQ2+97AFj++l3USd05LwBExKfITEvy8ZTSd3qwtsGqO+dkOvDHZL7IDj/mAOdk/+zVhhPX3c8weP1Vh8NXHE60w/hg151zcgqw9Yggd9hj2efp3axNXdfn3/WGua77EXBOdu4rALKTCM7Jvna8fYcAVx2xbxHwbuBerzR0S3fOCxHxF8DngU+llP65l2ocbLpzTi7q4PE0mU7iFwG393y5g0Z3zsvPgP3AZe3aL80+L+yZEged7pyTV4CqiGjfqf7s7PP6HqpRndf33/W5nscl3x7AcDKp+lkyQ8YvJ/MlsxIoO2K78WT6K9zQbv/vk7my8GHgEjJfSvvI9A3K+fvL10d3zguZSYNbyXxRndPuMSPX7y1fH939XengeA/hPHM5Py/A32bb/wF4C5mr2XuBW3P93vL10c3PryYy05IsJTPh8EXAtdm2hRwxV5qPEzo3v5t9/CuZK9J/kv35wmOdk2x7n37X5/wfKh8fZO6F/0/2l2UncCftJhLM/oIl4LPt2kvJzNH0SvbELgDenOv3NBAeJ3peyIyWTMd4PJTr95XPj+78rnRwLMNcPzgvZG6p/mU2fBwAVgM3AkNy/b7y+dHNczID+G9gLZlgvRT4ElCV6/eV74/jfTf0l+/6yP6lkiRJykP2mZMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYk6RuioiyiLghIp6MiJ0RkY7xGJnrWiUNPEW5LkCS8lk2oP0CmAY8A3wDKAGuAkYBB4E1wKaU0qu5qlPSwBUppVzXIEl5KyLuBy4B/gm4LmU/VCOiEVgGFAKjU0qbclelpIHM26ySdIIiYi6ZIDcfuD4d8X/HKaW1wC/J3AGZmZsKJQ0GhjlJOnH/K/t8c0qptYPXt2ef/ayV1Gv8gJGkE3cB0Ar8/Bivj80+L++bciQNRvaZk6QTEBEFZAY3vJpSGt3B6/XAS8DqlFJzX9cnafDwypwknZjD/ydcng127X2SzGfsLX1XkqTByDAnSScgO9jhaWA48N4jX4uI3wU+AbwAfLXvq5M0mHibVZJOUERcCdwBHAJuA9YCZwJvITMtyaUppVW5q1DSYGCYk6RuiIgrgL8BTgMCWAr8N/C1lNKuXNYmaXAwzEmSJOUx+8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxw5wkSVIeM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh77f/qANvlHjNKaAAAAAElFTkSuQmCC", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnMAAAFLCAYAAAC9cQpTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deZyddX33/9dn1kwyk9kz2SaZTPawBRLWgAgYQWuBaqlLb7W1Kl1skbZY+KnUoq22+BO19q7oXW9atdVKEVRUtoBokCXIEhLIQgJZgOz7npnv/cc5GZLJhMxkljNn5vV8PM7jZL7nuq58Ti7mnDfX9V0ipYQkSZLyU0GuC5AkSdKJM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh4zzEmSJOWxolwXkEt1dXWpqakp12VIkiQd15NPPrkxpVTfvn1Qh7mmpiYWLFiQ6zIkSZKOKyJe7qjd26ySJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcG9XJeknLjlvuW8pUHlh13u2sumcy1c6f0QUWSlL8ipZTrGnJm9uzZybVZpdx7962/BuD7V5+b40okqf+KiCdTSrPbt3ubVZIkKY8Z5iRJkvKYYU6SJCmPOQBCkuSgFCmPGeYkSVw7d8oRIc1BKVL+8DarJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxB0BIktRPOcpYnWGYkySpn3KUsTrD26ySJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwVIDSguRSOJGmgM8xpQHMpHEnSQOdtVkmSpDxmmJMkScpjhjlJkqQ8lvMwFxGNEXF7RGyLiO0RcUdEjOvC/tMj4gcRsTEi9kTEkoi4pjdrliRJ6i9yOgAiIoYC84B9wAeBBHwOeDAiTk0p7TrO/rOz+z8EfBjYBkwGynuxbEmSpH4j16NZPwI0A1NTSssBIuJZYBlwNfClY+0YEQXAvwMPpJR+57CXHuy9ciVJkvqXXN9mvRx49FCQA0gprQTmA1ccZ983AzN4g8AnSZI00OX6ytxJwF0dtC8CrjrOvudnn4dExKPALGAL8D3gb1JKe3qsSg0Idz61lqdWbWV/SytzvjCP6y6dypWnj8l1WZKkPNIfJ6PPdZirIRPA2tsMVB9n39HZ5+8DXwOuB2YDNwGNwO8cYz8NQnc+tZYb7ljI/pZWANZu3cMNdywEMNBJkjqtP05Gn+vbrJAZ9NBedGK/Q7V/J6V0Y0rpoZTSF4G/A66MiBkd7RQRH42IBRGxYMOGDSdYsvLNzfcsYc+BliPa9hxo4eZ7luSoIkmSekauw9wWMlfn2qum4yt2h9uUfb6vXfu92eeZHe2UUvpGSml2Sml2fX19pwtVfkkpsXrzbn78zCt89ieLWbu147vua7fuYdWm3X1cnSRJPSfXt1kXkek3194MYHEn9oWjr+wduqrX2o26lGe27t7PM2u28fSqrTyzZivPrN7Kpl37ASgtKqCksKDtFmt7b7r5QaY2VDB3RgNzZzRwyphKCgo6c3FYkqTcy3WY+xHwxYhoTimtAIiIJmAOmT5wb+RnZOanuwz4yWHtl2afF/Ropeo39h1sYfEr23l6dSa0PbNmGys3ZqYkjIBJ9eVcNG0EMxurmNlYxdSRFdz97KvccMfCI261lhUX8ldvzfR7uG/xOv73Q8v52oPLaRheylumZ4LduRNrKS0qzMn7lCSpM3Id5r4JfAy4KyI+ReYq22eB1cCthzaKiPHAi8BNKaWbAFJKmyLi88CnI2I7mcmDZwM3Av9++HQnfaU/jnDJd62tiZWbdh1xxW3xq9s50JK5IDuiopSZjVVcNXssM8dWccrYSiqGFB91nEODHD5x+7Psb2llTFXZEaNZP3xBM1t27WfeC+u5b/E6fvjUWr772CqGlRRy4dR65s5o4OKpDVQOPfrYkiTlUk7DXEppV0RcDNwCfJvMLdIHgI+nlHYetmkAhRzdx+8mYAfwp8BfA68CN5MJhH2uP45wyTcbduw77Ipb5nn73oMADCsp5NSxVfzR+c3MbKxkZmM1IyuHdPrYV54+hv96fBXQ8TmpHlbCu2aN5V2zxrL3QAuPvLiR+xav4/7n1/PTha9RWBCc1VTTdju2sWZoz7xpSZK6IddX5kgprQLedZxtXqKDEa4ppURm0mAnDs5Du/cfZOGabdnQto2nV29tG6hQWBBMG1nBO04bzcyxVcwcV8XE+nIK+6gv25DiQi6e1sDF0xr4+9bEM2u2ct/iddy3eB03/WQxN/1kMdNGVvDWGQ3MnTGSk8cMJ8J+dpKkvpfzMKfBoaU1sXTdDp5ZvZWns4+l63bQmh2+0lhTxunjqvjDOU3MbKzipNGVlJX0j75qBQXB6eOqOX1cNZ+4bBorN+7i/myw+9qDy/nqvOWMqhzS1s/unOZaSopyPVBckjRYGObU41JKvLJt7xHB7bm129i9PzP4oLKsmNMaq3jrSSOZ2VjJqWOrqCsvzXHVnTehbhgfeVMzH3lTM5t27mvrZ3f7k2v49qMvU1Fa1NbP7s1TR1BZZj87SVLvMcyp27btOcDCNdt4evUWnl6duW26Ycc+AEoKC5gxeji/N7uRmY1VnNZYRVPt0AFzS7K2vJSrZjdy1exG9h5o4VfLMv3sHnhhHT959lWKCoKzm2uYO72BuSeNZExVWa5LliQNMIY5dcn+g6288Nr2tituz6zeyosbdrW93lw/jAsm17VNCzJt5PBBc8txSHEhb5nRwFtmNNDSmnh69Rbuzd6O/cyPF/OZHy9mxqjhbQMoThptPztJUvcZ5nRMKSVe3rS7Lbg9vXori1/Z3jb5bl15ZlqQ3zl9DDMbqzllbKW3FLMKC4JZ42uYNb6GG942nRc37GzrZ/fVecv4ygPLGFNVxlumj2DujJGc3VxDceHgCL2SpJ5lmFObTTv38cyarTydHVn6zOqtbNtzAMhMsHvK2Er+IDtA4bTGKkZXDvHKUidNrC9n4oXlXH3hRDbu3Me859dz7+J1fO+J1fz7r1+mYkgRF00dke1nV9/hXHmSJHXEMDdI7dnfwqJXtr1+u3TNVlZvzkwLUhAwpaGCt508si24TR5RTpFXjnpEXXkpv3dmI793ZiN79rfwy2Ubsv3s1vOjZ16huDA4p7mWt85o4JLpDYy2n50k6Q0Y5gaBltbEixt2HtHP7YXXdtCSnRdkTFUZpzVW8v5zxnNadhWFoSX+p9EXykoKeetJI3nrSSNpaU38ZtWWtvnsPn3XIj591yJOHjOcudNHMndGA9NHVXg1VJJ0BL+xB6DXtu09IrgtXLuNnfsyqyhUDCliZmMVf3LhRE5rrOK0xkpGVHR+FQX1nsKC4MymGs5squGGt03jxQ072wZQfPmBpdxy/1LGVpfxlukNvHVGA2dOsJ+dJMkwl/d27D3AwrWv93F7evVW1m3PTAtSXBjMGDWcd54xhtOyqyhMqB1GQR+toqATFxFMGlHBpBEV/OmbJ7F+x14eeH499y9ex38+vorbHnmJyrJiLppaz9wZI7lwaj3lpf46S9Jg5Kd/HjnQ0sqS13YcEdyWb9hJyq6iMKFuGOc213JadlqQ6aOGM6S4f6yioO4ZUTGE9541jveeNY7d+w/y8NLMfHbzXljHnU+/QklhAedOrG2b9qRhuFdbJWmwMMz1Uykl1mzZw1OHBbfn1m5j38HMtCA1w0qY2VjFb582OnO7dGwlVUNLcly1+sLQkiIuO3kkl508koMtrTz5craf3fPr+NSdz/GpO5/jtLGV2WA3kikN5fazk6QBzDDXT2zdvT97xS2zksIza7axedd+AEqLCjhlTHaAQvaq29jqMr+gRVFhAWc313J2cy2f/K3pLFu/k/sWr+Pexev44r1L+eK9SxlXM7Rt3dgzm6odlSxJA4xhrpfc+dRanlq1lf0trcz5wjyuu3QqV54+BoC9B1pY/Or2tituz6zeykubdgMQAZNHlPOW6SPagtuUhgo7uuu4IoIpDRVMaajgzy6axLrte7n/+cwAiu88+jLfmr+SqqHFXJydz+5NU+oZZj87Scp7fpL3gjufWssNdyxsWylh7dY9XHf7M/xgwWp27DvI869u50BLpqPbqMohnDa2inefOY6ZjZlpQezIrp7QMHwIv3/2eH7/7PHs3HeQh5e+Pp/dHU+tpaSogDkTa5k7YyRvmT6CEfazk6S8ZGroBTffs4Q9B1qOaDvQknjkxU2cO7GWD1/Q3LZ2qR3V1RfKS4t4+ymjePspozjY0soTLx3qZ/caD/5wIf/fD2FmYxVzZ2SmPZk0wn52kpQvDHO94JWte4752n9+5Jw+rEQ6WlF25Ou5E2v59Dums2TdDu5blBlAcfM9S7j5niU01Q5tG0Axa3w1hU5nI0n9lmGuF4yuKmNtB4HOZZnU30QE00YOZ9rI4fz5JZN5ddse7n9+PfctXsdtj7zEN3+5kuqhxVw8rSHbz67O1UEkqZ/xU7kXXHfpVG64Y+ERt1rLigu57tKpOaxKOr5RlWW8/5zxvP+c8ezYe4BfZPvZ3bf4Nf7nN2soLSrg/El1zM2uG1tfUZrrkiVp0DPM9YJDo1Y/cfuz7G9pZUxV2RGjWaV8UDGkmHecOpp3nDqaAy2tPL5yc9u6sQ+8sJ6IhZzeWMXcGZl1YyeNKM91yZI0KBnmesmVp4/hvx5fBcD3rz43x9VI3VNcWMCcSXXMmVTH3/72DJ5/dUfbAIp//PkL/OPPX6C5bljbChSnj7OfnST1FcOcpC6JCGaMHs6M0cO55i2TeWXrnrb57P7tVyu59eEV1A4r4ZLpI5g7YyTnT6qjrMRl5SSptxjmJHXL6KoyPnBuEx84t4ntew/w0JJMP7ufLXyN/16whiHFBZw/qZ63zmjg4ukjqCu3n50k9STDnKQeM3xIMZefNprLTxvN/oOtPLZyE/ctXsf9i9dx//PriIBZ46rbbsc215e/4WopkqTjM8xJ6hUlRQVcMLmeCybX83eXn8SiV7Zzb3YAxed/9gKf/9kLjKgoYfOuAxxszayIsnbrHm64YyGAgU6SOskwJ6nXRQQnj6nk5DGV/OXcKazevJsHns+EukNB7pA9B1q48a7nGF5WxJSGCsZUlbkahSS9AcOcpD7XWDOUP5gzgb/78eIOX9++9yAfum0BkFmKbHJDOVMbKpjSUMHUkZnnuvISQ54kYZiTlEPHWi1lVOUQ/vm9p7Nk3Q6WvraDJet2cM+i1/jeE6vbtqkZVsKUQyFvZAVTGyqY3FBBZVlxX74FSco5w5yknDnWail/c9k0ZjfVMLuppq09pcTGnftZum4HS17bkXlet4Pbn1zDrv2v7z+qcsgRV/CmNlQwaUS506NIGrAMc5JypiurpUQE9RWl1FeUMmdSXVt7Som1W/dkQ97OtrD36xWb2H+wNbsvjK8ZemTIG1nBhLphFBcW9M2blaReYpiTlFPdXS0lIhhbPZSx1UO5eFpDW/vBllZe3rybZYeHvHU7eOCF9bRkB10UFwbNdeXZ27TlbSGvsXooBa5gISlPGOYkDUhFhQVMrC9nYn05l538evu+gy2s2LDriNu1T6/ewo+feaVtm7LiQiYfCneH9clrGF7qoAtJ/Y5hTtKgUlpUyPRRw5k+avgR7bv2HWTZ+p1tAy6WrtvBw0s3cPuTa9q2GT6k6IjbtFOyI2xrhpX09duQpDaGOUkChpUWMbOxipmNVUe0b9mVGXRx6Dbt0td28uNnXuG7jx1s26a+opQp7a7kTWmooLzUj1hJvc9PGkl6A9XDSji7uZazm2vb2lJKrN+x7/VRtdnn7z2++oiRuWOqyg67kpcJexPryxlS7MhaST3HMNeDbrlvKV95YNlR7U3X333Ez9dcMplr507pq7Ik9bCIoGH4EBqGD+FNU+rb2ltbE2u27Gm7TXso5P1y2QYOtGQGXRQENNUNO2oS5KbaoRQ5slbSCch5mIuIRuAWYC4QwP3Ax1NKqzqxbzrGS6enlJ7uuSo759q5Uwxp0iBWUBCMqx3KuNqhzJ3x+sjaAy2tvLRx1xGTIL/w2g5+vug1UvZTrKSwgIkjyjOjarMDLg4tZ+bIWklvJKdhLiKGAvOAfcAHgQR8DngwIk5NKe3qxGFuA25t17a0J+uUpO4oLixgcnaFCk59vX3P/hZe3LDziEmQH1+5mTuffn1k7bCSQia3G1U7ZWQ59eWOrJWUkesrcx8BmoGpKaXlABHxLLAMuBr4UieOsTal9GjvlShJvaOspJCTx1Ry8pjKI9q37z1w5Px4r+3gvufX8f0Fry9nVj20+KhJkKeMqKByqMuZSYNNrsPc5cCjh4IcQEppZUTMB66gc2FOkgaU4UOKmTW+hlnja45o37hz3xFTpyx5bQd3/GYtO/e9PrJ25PAhR02CPGlEOUNLcv1xL6m35Pq3+yTgrg7aFwFXdfIYfxIR1wEtwKPA36aUftlD9UlSv1FXXkrdpFLOa7ec2Svb9r4e8rLP/7FiE/sOW85s3KHlzA67XTuhbhglRUcPurjzqbU8tWor+1tamfOFecdcYk1S/5DrMFcDbOmgfTNQ3Yn9vwP8BHgFGA9cB8yLiLkppYd6qkhJ6q8igjFVZYypKuOiaSPa2ltaE6s27z6iP97S13Yw77DlzIoKggl1w44YcLF6826+dN8S9rdkguDarXu44Y6FAAY6qZ/KdZiDzKCH9jrVqzel9P7DfvxlRNwFPEdmEMX5He0TER8FPgowbty4rlUqSXmiMBvUJtQN47KTR7a17zvYwsqNuw6bI28nC9ds4+5nXz3msfYcaOHme5YY5qR+KtdhbguZq3PtVdPxFbs3lFLaERF3A3/0Btt8A/gGwOzZs481tYkkDUilRYVMGzmcaSOPXs5s+fqdXPEv8zvcb+3WPdx413OcN7GOc5trHWgh9SO5DnOLyPSba28GsPgEjxl0fLVPknQMw0qLOK2xijFVZazduueo10uLCvjBgjX8x69fpiDglDGVnDepjjkT65jdVO2qFlIO5TrM/Qj4YkQ0p5RWAEREEzAHuL6rB4uI4cBvAY/1YI2SNGhcd+lUbrhj4RHLkpUVF/L5d57C208ZxVOrtjD/xU08snwj33x4Bf/60IuUFBUwa1w1cybVct6kOk4dU+lqFlIfynWY+ybwMeCuiPgUmStqnwVWc9hEwBExHngRuCmldFO27a+BqcCDvD4A4q+BkcDv9+F7kKQB41C/uE/c/iz7W1oZU1V2xGjWQ+vU/uXcKezcd5AnVm5m/vKNzH9xE1+8dyncu5SK0iLObq5hzqQ65kyqY/KIcic4lnpRTsNcSmlXRFxMZjmvb5O5RfoAmeW8dh62aQCFwOH/q7cE+J3soxLYDswH/iil9HgflC9JA9KVp4/hvx7PrKj4/avPPeZ25aVFXDRtRNso2k079/HrFZuYv3wT85dv5P7n1wNQX1HKeRNrmTOxjvMm1TK2emjvvwlpEMn1lTmya7C+6zjbvES7Ea4ppR8DP+69yiRJXVFbXso7Th3NO04dDcDqzbt55MWN2XC3ibuyy5SNrx3KeRPrMrdlJ9ZRM6wkl2VLeS/nYU6SNDA11gzl3TXjePeZ40gpsXTdTuYv38gjL27kx8+80nb1b8ao4W397c5qqmFYqV9NUlf4GyNJ6nURwdSRmeXFPnT+BA62tPLs2m08sjxz5e7fH3mZb/5yJUUFwenjqrJX7uqY2VjV4SoVkl5nmJMk9bmiwgLOGFfNGeOq+djFk9mzv4UFL29m/vJNPPLiRr46bxlfeWAZQ0sKObOphvMnZfrbTR85nIICB1NIhzPMSZJyrqykkAsm13PB5HoAtu0+wK9XbMr2udvI3//0eQBqhpVwbnMt503KDKgYXzvUkbIa9AxzkqR+p3JoMZedPLJtKbLXtu3NToGykUeWb+LuhZnlx8ZUlWVGymav3I2oGJLLsqWcMMxJkvq9kZVDeNessbxr1lhSSqzYuKutv929i9fxgyfXADClobytv93ZzTUMH+KyYxr4DHOSpLwSEUysL2difTnvP7eJltbE4le2Mz97S/Z7T6zitkdeoiDg1LFVzMnekj1jvMuOaWAyzEmS8lphQXDK2EpOGVvJH184kX0HW/jNy1vb+tt9/Rcr+JcHX6S0qIDZTdWZlSkm1nHymEoKHUyhAcAwJ0kaUEqLCjl3Yi3nTqzlr946lR17D/D4ytdHyv7Tz5cASxg+pIhzmmuzy47VMrHeZceUnwxzkqQBrWJIMZdMb+CS6Q0AbNixj0eyAynmv7iRexevA6BheCnnTaxrG1Axuqosl2VLnWaYkyQNKvUVpVwxcwxXzBwDwKpNu9v62z28dAM/fGotAM11w9qmQDl3Yi1VQ112TP2TYU6SNKiNqx3KuNpxvPescbS2Jpas25FddmwTP/zNWr7z6Coi4KTRw5kzsa5t2bGyEgdTqH8wzEmSlFVQEEwfNZzpo4bz4QuaOdDSyjOrtzI/e0v2W/NXcuvDKyguDE4fV8352f52p46torjQZceUG4Y5SZKOobiwgNlNNcxuquGat0xm9/6DPPHSlswcdy9u5Jb7l/Kl+2BYSSFnN9e29beb2lDhsmPqM50KcxFRCvwS2AFcllI6cIztSoCfAcOAC461nSRJ+WhoSREXTqnnwimZZce27NrPoys2ZfvcbWLeC+sBqB1WwrnZYHf+pDoaa4bmsmwNcJ29Mvf7wCzgt98ooKWU9kfEzcBPs/vc1u0KJUnqp6qHlfC2U0bxtlNGAfDK1j1t/e3mL9/IT57NLDvWWFPW1t/uvIm11JWX5rJsDTCdDXPvBFaklH56vA1TSj+PiGXAVRjmJEmDyOiqMq6a3chVsxtJKfHihp2Z/nbLN3L3wlf53hOrAZg2siK77FgtZzfXUl5qryeduM7+13M6mattnfUw8PaulyNJ0sAQEUwaUcGkERV88LzMsmPPrd3Gr5Zv5JEXN/Ldx17mW/NXUlgQzGysYs7EWs6bVMfp46ooLXKkrDqvs2GuDljXheOuA2q7Xo4kSQNTYUFwWmMVpzVW8WcXTWLvgRZ+8/KWtv52X3twOV+dt5whxQWc2VTTtuzYjNHDKSwI7nxqLU+t2sr+llbmfGEe1106lStPH5Prt6V+oLNhbg9Q3oXjlgN7u16OJEmDw5Diwkwfukl1XHcpbN97gMdWbM72udvIF372AgCVZcWMrynj+dd2cKAlAbB26x5uuGMhgIFOnQ5zq4Ezu3Dc2cCqrpcjSdLgNHxIMXNnNDB3RmbZsfXb97YNpLjjN2tpSemI7fccaOEffvo8V8wc7Zqyg1xnZzh8CDgnImYfb8OImAWcBzzYjbokSRrURgwfwpWnj+Hmq06jtV2QO2T9jn3M+tz9XP3tBXzrVyt5bu02Wlo73lYDV2evzH0N+BPgBxHx9pTS8x1tFBHTgB8ALcD/7pkSJUka3EZXlbF2656j2quGFnPxtBE8tnIT9yzKdG2vGFLEmU01nDWhhrMn1HDymEpXpxjgOhXmUkpLIuIm4DPAUxFxOzAPWAMkYCxwCfAuoBS4MaW0pFcqliRpkLnu0qnccMdC9hxoaWsrKy7kM799UlufuVe27uGJlzbz6IrNPL7y9QmMh5YUMmt8NWdlA95pjVUMKXa07EDS6YltUko3RcRB4G+B9wHvbbdJAAeAT6aUPt9zJUqSNLgdCmyfuP1Z9re0Mqaq7KjRrKOryrhi5hiumJlp27BjH0+8tJnHV27m0RWb+NL9S0kJSooKmNlYxTkTajhrQi1njK9iaInz3OWzLp29lNI/RMR3gQ8Bc4BRZELcK8CvgP+bUnq5x6uUJGmQu/L0MfzX45mxhd+/+tzjbl9fUcrbTxnF27OrU2zdvZ8nXtrC4ys38djKzXztweW0zltOUUFwythKzppQwzkTapnVVM3wIcW9+l7Us7ocxbNh7W97oRZJktRLqoaWHDFadue+gzz58hYeW7GJx1du5lu/Wsmtv1hBQcD0UcM5e0ItZ03I3JqtGVaS4+r1RryuKknSIFReWsSFU+q5cEo9QGYS41VbeHzlZh5bsblthQqAKQ3l2QEVtZw9oYYRw4fksnS106kwFxErgC+nlL56WNs4oCml9HBvFSdJkvrGkOJCzptYx3kT6wDYf7CVhWu3ZgdUbOaHv1nLdx7N3OadUDesbUDF2c01jK0emsvSB73OXplrAqratf0hcCPgkBhJkgaYkqICZo2vYdb4Gv7sIjjY0sriV7fz2IrNPLZyMz9f9BrfX7AagDFVZZydvSV71oQaJtQNcyLjPuRtVkmSdFxFhQWcOraKU8dW8ZE3NdPamliybkfmtuzKTTy8bAN3PLUWyAy+ODTP3dkTapk8opyCAsNdbzHMSZKkLisoCKaPGs70UcP54HlNpJRYsXEXj2XnuXts5WbufvZVAKqHFrdNZHxOcy3TRw2n0HDXYwxzkiSp2yKCifXlTKwv531njyOlxJote3hs5ebMiNmXNnPv4uwqFaVFzGqqbhsxe+pYV6noDsOcJEnqcRFBY81QGmuG8ruzxgLw2ra9PLZyU/bW7GYeWvICkFnN4ozxVZzVVMvZzTXMdJWKLulKmHtzu86MbwaIiE+TmTi4vZRS+uyJlyZJkgaSkZVDjlilYuPOfTyRDXaPr9zMlx9YSrofSgozq1QcGlAxa3w1w0q9/nQsXQpz2Ud7f9fu50Qm3CXguGEuIhqBW4C52f3uBz6eUlrVhdqIiBuAfwDmp5TO78q+kiSp79WVl/K2U0bxtuwqFdt2H2DBy5lw99jKzfzrL17kaw8up7AgOHlMZXYJshpmN9VQWeYqFYd0Nsy1D2w9IiKGAvOAfcAHyQTAzwEPRsSpKaVdnTxOM/BJYH1v1ClJknpf5dBiLpnewCXTX1+l4jcvb2kbMft/57/ErQ+vIAKmjxyeHVBRw5lNNdSWl+a4+tzpVJhLKfVKmAM+AjQDU1NKywEi4llgGXA18KVOHudfge8CU7EfoCRJA0J5aRFvmlLPmw5bpeKpVVt5fOVmHs/hPh8AABQ1SURBVH9pE997YhW3PfISAJNGlLfNdXdOcy0Ng2iVilwHn8uBRw8FOYCU0sqImA9cQSfCXES8DzgDeC9wR28VKkmScmtIcSHnTqzl3Im1wOTsKhXb2gZV3PX0K3z3sUwvrfG1Q7PhLrME2djqsgE7kXGnw1xEvAOoBr6XUjpwjG1KgHcDm1NKd3fisCcBd3XQvgi4qhM1VZPpb/eJlNLmgXqSJEnS0TKrVFQza3w1f/rmzCoVz7+6oy3c3bt4Hf+9YA0AoyuHZAdUZEbMNg+gVSo6uzbrKcCdwL+mlL59rO1SSvsjYjbwpxExM6W06DiHrgG2dNC+mUxwPJ6bgaXAbZ3YVpIkDWBFhQWcMraSU8ZW8uELMqtULFu/k8eykxj/avkm7nz6FSAz+OLwJcimNlTk7SoVnb0y9yFgP50bCPFZ4MPZx7Wd2D510Hbcf82IuAD4AHBGSqmjYxxrv48CHwUYN25cZ3eTJEl5pqAgmDqygqkjK/jAuZlVKlZu3NU2z91jKzZx98LMKhWVZZlVKs5pzoS7GaOGU5QnExl3Nsy9GXgopbTxeBumlDZGxEPAxZ047hYyV+faq6bjK3aHuxX4N2BNRFRl24qAwuzPe1JK+zqo7xvANwBmz57d6RAoSZLyW0TQXF9Oc3057zkrc0Fn9ebdmQEV2RGz9z+fWaWivLSIWeOr20bMnjKmipKiI8PdnU+t5alVW9nf0sqcL8zjukuncuXpY/r8fXU2zE0EHujCcZ8HOjPX2yIy/ebamwEsPs6+07OPP+7gtS1krgp+uRM1SJKkQerQKhXvyq5SsW773uwkxpt4bMVmbr5nCQBDigs4Y1x1223ZtVv2cONdi9jf0grA2q17uOGOhQB9Hug6G+ZKyNxm7az9QGcmfPkR8MWIaE4prQCIiCZgDnD9cfa9qIO2LwOFwJ8Dyzt4XZIk6Zgahg/h8tNGc/lpowHYtHMfT7y0pW1QxVceWMaxOnftOdDCzfcs6bdhbjPQlQ5m44BNndjum8DHgLsi4lO8vmrEajK3UQGIiPHAi8BNKaWbAFJKD7U/WERsBYo6ek2SJKmrastLuezkkVx28kgAtu05wJMvb+ZDty3ocPtXtu7py/IA6GzPvt8AcyPiuFfbImIImaW5njrettkVHi4mMyL122Qm/l0JXJxS2nn4YclcccuPnoiSJGlAqiwr5uJpDYypKuvw9dHHaO9NnQ1HPwTqgb/vxLY3AXXA/3TmwCmlVSmld6WUhqeUKlJKV6aUXmq3zUsppUgpfeY4x3qz67JKkqTedt2lUykrLjyiray4kOsundrntXQ2zP0H8AJwbUR8OyImt98gIiZHxH8Af0VmAMQx56OTJEnKZ1eePobPv/MUSrLTl4ypKuPz7zyl/45mTSkdiIgryYxo/X3gfRGxFlhDpp/b2Owjsm1XppQO9k7JkiRJuXfl6WP4r8czy4d9/+pzc1ZHp/ugpZSWArPIzO22n0x4Owc4F2jMtv0fYPbha61KkiSp93R6bVaAlNJ64CMR8efAbGAUmatxrwALUkp7e75ESZIkHUuXwtwh2dD2qx6uRZIkSV3UqdusEfGBiDi1t4uRJElS13S2z9xtwJWHN0TEByNiXo9XJEmSpE7rziS8TcCFPVSHJEmSToArKkiSJOUxw5wkSVIeM8xJkiTlsa6EudRrVUiSJOmEdGWeuc9ExGfaN0ZEyzG2TymlE5rHTpIkSZ3TlbAVXTx2V7eXJElSF3UqzKWU7FsnSZLUDxnSJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLymGFOkiQpjxnmJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLymGFOkiQpjxnmJEmS8lhRrguQetMt9y3lKw8sO6q96fq7j/j5mksmc+3cKX1VliRJPcYwpwHt2rlTDGmSpAHN26ySJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHch7mIqIxIm6PiG0RsT0i7oiIcZ3Yb3xE3BURL0fEnojYGBEPRcTb+qJuSZKk/iCnYS4ihgLzgGnAB4H3A5OBByNi2HF2Lwc2Ap8C3g78EbAT+GlEvLPXipYkSepHcj1p8EeAZmBqSmk5QEQ8CywDrga+dKwdU0qLyAS4NhFxN7AS+EPgjl6qWZIkqd/I9W3Wy4FHDwU5gJTSSmA+cEVXD5ZSOghsAw70WIWSJEn9WK7D3EnAcx20LwJmdOYAEVEQEUURMTIiPg1MAf6lB2uUJEnqt3J9m7UG2NJB+2agupPH+Cfgr7J/3gm8J6X0QA/UJkmS1O/l+socQOqgLbqw/5eBM4HfBn4G/GdEvONYG0fERyNiQUQs2LBhQ9cqlSRJ6mdyHea2kLk61141HV+xO0pKaU1KaUFK6Scppd8DHgW++AbbfyOlNDulNLu+vv6EipYkSeovcn2bdRGZfnPtzQAWn+AxFwAfP+GKJPW6W+5bylceWHZUe9P1dx/x8zWXTObauVP6qixJyku5DnM/Ar4YEc0ppRUAEdEEzAGu7+rBIqIAOB94sQdrlNTDrp07xZAmST0k12Hum8DHgLsi4lNk+s99FlgN3Hpoo4gYTyag3ZRSuinb9hkyt2jnA68BI8nMO3cW8L6+ewuSJEm5k9Mwl1LaFREXA7cA3yYz8OEB4OMppZ2HbRpAIUf28fsNmdup7wEqyQS6Z4ALUkrz+6B8SZKknMv1lTlSSquAdx1nm5doN8I1pfQjMrdpJUmSBq1cj2aVJElSNxjmJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLymGFOkiQpjxnmJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLymGFOkiQpjxnmJEmS8phhTpIkKY8Z5iRJkvKYYU6SJCmPGeYkSZLymGFOkiQpjxnmJEmS8phhTpIkKY8Z5iRJkvJYUa4LkCTl3i33LeUrDyw7qr3p+ruP+PmaSyZz7dwpfVWWpE4wzEmSuHbuFEOalKe8zSpJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMcMc5IkSXnMMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJecwwJ0mSlMdyHuYiojEibo+IbRGxPSLuiIhxndhvdkR8IyJeiIjdEbEqIr4bERP6om5JkqT+IKdhLiKGAvOAacAHgfcDk4EHI2LYcXZ/D3AS8FXgbcD1wBnAgoho7LWiJUmS+pGiHP/9HwGagakppeUAEfEssAy4GvjSG+z7jymlDYc3RMR8YGX2uDf2SsWSJEn9SK5vs14OPHooyAGklFYC84Er3mjH9kEu2/YysAEY08N1SpIk9Uu5DnMnAc910L4ImNHVg0XEdGAE8Hw365IkScoLuQ5zNcCWDto3A9VdOVBEFAFfJ3Nl7t+6X5okSVL/l+swB5A6aIsTOM7XgPOA/5VS6iggZg4c8dGIWBARCzZsOOpOrSRJUl7JdZjbQubqXHvVdHzFrkMR8Xngo8CHUkr3vtG2KaVvpJRmp5Rm19fXd6lYSZKk/ibXo1kXkek3194MYHFnDhARnyQzLclfpJS+3YO1SZIk9Xu5vjL3I+CciGg+1BARTcCc7GtvKCL+Avgc8MmU0j/3Uo2SJEn9Vq7D3DeBl4C7IuKKiLgcuAtYDdx6aKOIGB8RByPixsPa3gN8Gfg5MC8izjns0eWRsJIkSfkop7dZU0q7IuJi4Bbg22QGPjwAfDyltPOwTQMo5MjweVm2/bLs43C/AN7cS2VLkiT1G7nuM0dKaRXwruNs8xLtRrimlP4A+IPeqkuSJCkf5Po2qyRJkrrBMCdJkpTHDHOSJEl5zDAnSZKUxwxzkiRJeSzno1klSVLHbrlvKV95YNlR7U3X333Ez9dcMplr507pq7LUzxjmJEnqp66dO8WQpuPyNqskSVIeM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxpyaRJEnqpP4491+klPrkL+qPZs+enRYsWJDrMiRJko4rIp5MKc1u3+5tVkmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPDao12aNiA3Ay73819QBG3v571DXeE76J89L/+M56Z88L/1PX52T8Sml+vaNgzrM9YWIWNDRorjKHc9J/+R56X88J/2T56X/yfU58TarJElSHjPMSZIk5THDXO/7Rq4L0FE8J/2T56X/8Zz0T56X/ien58Q+c5IkSXnMK3OSJEl5zDB3AiKiMSJuj4htEbE9Iu6IiHGd3HdIRNwcEa9GxJ6I+HVEvKm3ax4MTvS8RMTsiPhGRLwQEbsjYlVEfDciJvRF3QNZd35X2h3nhohIEfGr3qhzsOnueYmI6RHxg4jYmP0cWxIR1/RmzQNdN79XxkXEv2c/u3ZHxNKI+FxEDOvtugeyiBgbEf+c/Z7enf0Maurkvn36XW+Y66KIGArMA6YBHwTeD0wGHuzkL86/AR8BbgTeAbwK3BMRM3un4sGhm+flPcBJwFeBtwHXA2cACyKisdeKHuB64Hfl0HGagU8C63ujzsGmu+clImYDjwGlwIeBtwP/P1DYWzUPdN05J9nX7wfeBHwa+C3g/wB/BXyrF8seDCYBvwdsAX7ZxX379rs+peSjCw/gGqAFmHRY2wTgIPCXx9n3NCABf3hYWxGwBPhRrt9bPj+6eV7qO2gbD7QCN+X6veXrozvnpN1x7gFuBR4CfpXr95Xvj27+rhQAi4Af5vp9DKRHN8/JW7PfK29t1/6F7P5Dc/3+8vUBFBz25w9n/52bOrFfn3/Xe2Wu6y4HHk0pLT/UkFJaCcwHrujEvgeA7x+270Hge8ClEVHa8+UOGid8XlJKGzpoexnYAIzp4ToHk+78rgAQEe8jc5X0hl6pcHDqznl5MzAD+FKvVTc4deeclGSft7dr30omfEdPFTnYpJRaT3DXPv+uN8x13UnAcx20LyLzIXe8fVemlHZ3sG8JmUu6OjHdOS9HiYjpwAjg+W7WNZh165xERDVwC/CJlNLmHq5tMOvOeTk/+zwkIh6NiAMRsT4ivhoRZT1a5eDSnXNyP7AM+MeImBER5RFxMZmrfV9PKe3q2VLVCX3+XW+Y67oaMvfP29sMVHdj30Ov68R057wcISKKgK+TuTL3b90vbdDq7jm5GVgK3NaDNal752V09vn7wL3AXOCfyNyC+s+eKnAQOuFzklLaSyZkH7oFvgN4APgJ8LGeLVOd1Off9UU9fcBBoqPJ+TpzKTu6sa+Or6f+bb8GnAf8Vkqpo19Idd4JnZOIuAD4AHBGynY4UY860d+VQxcAvpNSujH754ciohD4QkTMSCkt7pEKB58T/V0ZQiZcjyAzcGIVcBaZjvcHgT/pwRrVOX3+XW+Y67otdJyqq+k4iR9uM9DRUPPqw17XienOeWkTEZ8HPgp8MKV0bw/VNlh155zcSuaq6JqIqMq2FQGF2Z/3pJT29Vilg0t3zsum7PN97drvJdPhfiZgmOu67pyTPyLTl3FSSunFbNvDEbEN+EZEfD2l9EyPVarO6PPvem+zdt0iMvfD25vB8T/EFgETssPQ2++7H1h+9C7qpO6cFwAi4pNkpiW5JqX07R6sbbDqzjmZDvwxmS+yQ485wDnZP3u14cR19zMMjr7qcOiKw4l2GB/sunNOTgG2HBbkDnk8+zy9m7Wp6/r8u94w13U/As7Jzn0FQHYSwTnZ1463bzFw1WH7FgHvBu71SkO3dOe8EBF/AXwO+GRK6Z97qcbBpjvn5KIOHs+Q6SR+EXB7z5c7aHTnvPwM2Adc1q790uzzgp4pcdDpzjl5DaiOiPad6s/OPq/toRrVeX3/XZ/reVzy7QEMI5OqF5IZMn45mS+ZFUD5YduNJ9Nf4cZ2+3+PzJWFDwOXkPlS2kumb1DO31++PrpzXshMGtxK5ovqnHaPGbl+b/n66O7vSgfHewjnmcv5eQH+Ntv+D8BbyFzN3gPcluv3lq+Pbn5+NZGZlmQpmQmHLwKuy7Yt4LC50nyc0Ln53ezjX8lckf6T7M8XHuucZNv79Ls+5/9Q+fggcy/8f7K/LDuAO2k3kWD2FywBn2nXXkZmjqbXsif2MeDNuX5PA+FxoueFzGjJdIzHQ7l+X/n86M7vSgfHMsz1g/NC5pbqX2bDx37gZeAmoDjX7yufH908JzOA/wZWkwnWS4EvAtW5fl/5/jjed0N/+a6P7F8qSZKkPGSfOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SZKkPGaYkyRJymOGOUmSpDxmmJMkScpjhjlJkqQ8ZpiTJEnKY4Y5SeqmiCiPiBsj4qmI2BER6RiPEbmuVdLAU5TrAiQpn2UD2i+AacCzwNeBUuAqYCRwAFgFbEwprc9VnZIGrkgp5boGScpbEXE/cAnwT8D1KfuhGhGNwDKgEBiVUtqYuyolDWTeZpWkExQRc8kEufnADemw/ztOKa0GfknmDsjM3FQoaTAwzEnSiftf2edbUkqtHby+LfvsZ62kXuMHjCSduAuAVuDnx3h9bPZ5ed+UI2kwss+cJJ2AiCggM7hhfUppVAevNwCvAC+nlJr7uj5Jg4dX5iTpxBz6P+GKbLBr7xNkPmNv7buSJA1GhjlJOgHZwQ7PAMOA9x7+WkT8LvBx4AXgK31fnaTBxNusknSCIuJK4A7gIPB9YDVwJvAWMtOSXJpSWpm7CiUNBoY5SeqGiLgC+BvgNCCApcB/A19NKe3MZW2SBgfDnCRJUh6zz5wkSVIeM8xJkiTlMcOcJElSHjPMSZIk5THDnCRJUh4zzEmSJOUxw5wkSVIeM8xJkiTlMcOcJElSHjPMSZIk5bH/Bw5mQ/LVCjaKAAAAAElFTkSuQmCC", "text/plain": [ "
" ] @@ -5869,11 +5894,6 @@ } ], "source": [ - "# Need to plot 1 , 3 ,5 ... as didn't generate shuffle error bars for 0.1, 0.3, etc.\n", - "# Honestly, these aren't error bars per say, more like \"shuffle ranges\" \n", - "array = np.diag(df_pivot)\n", - "array = array[0:11:2]\n", - "\n", "plt.errorbar([0,0.2,0.4,0.6,0.8,1.0],array,yerr=(yerr_2,yerr_1),fmt='-o',capsize=5)\n", "plt.xlabel(\"$\\sigma$\",fontsize=20)\n", "plt.ylabel(\"FCFC\",fontsize=20)\n", @@ -5884,7 +5904,7 @@ }, { "cell_type": "code", - "execution_count": 131, + "execution_count": 17, "metadata": {}, "outputs": [ { @@ -5915,6 +5935,46 @@ "plt.savefig(\"do-not-track\\\\Het_sweep_reviewer2.pdf\",bbox_inches='tight')" ] }, + { + "cell_type": "code", + "execution_count": 98, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAn0AAAFLCAYAAACjuDrgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdeVxVdf7H8deXHVwQFVBAwBXEFUXTcAHUbFUrayrbm6yZalqtzBZbdbJyaprJcipb+FWj2b664b4E7rgvqGjuu7ggfH9/XHSQUC6KHpb38/G4j8v9nu8553O6ph++q7HWIiIiIiKVm4fTAYiIiIjI+aekT0RERKQKUNInIiIiUgUo6RMRERGpApT0iYiIiFQBSvpEREREqgAvpwOoCOrWrWujo6OdDkNERESkRBkZGTuttcFFy5X0uSE6Opr09HSnwxAREREpkTFmQ3Hl6t4VERERqQKU9ImIiIhUAUr6RERERKoAx8f0GWMaACOBXoABJgIPWWs3unl+c+AFIBmoBmwE/m2tfbNQHQ/gCeAeoB6wEnjBWvtlGT6KiIhIuZWbm0t2djZHjhxxOhQpI35+fkRERODt7e1WfUeTPmNMADAZOArcBljgJWCKMaa1tfZQCecnFJyfBvwZ2Ac0BaoXqfoi8BgwBMgAbgDGGmOutNb+WGYPJCIiUk5lZ2dTo0YNoqOjMcY4HY6cI2stu3btIjs7m4YNG7p1jtMtfXcDjYAYa+0aAGPMYmA1rla5N053YkHr3UfAJGvt1YUOTSlSLwRXwjfcWvvaiTrGmCbAcEBJn4iIVHpHjhxRwleJGGOoU6cOO3bscPscp8f09QHmnEj4AKy164GZQN8Szk0C4jhDYligN+ADfFqk/FOglTHGvfRYRESkglPCV7mU9vt0OulrASwtpjwTV0J3Jl0K3v2MMXOMMbnGmO3GmLeMMf5F7nEUWFPk/MyC95LuI1XQ1ws2kzh8Mg2f/IHE4ZP5esFmp0MSEanwjDHccsstJz8fP36c4OBgrrzyylJdJzo6mp07d55znarG6aSvNrCnmPLdQFAJ54YVvH8B/IprIsiruMb2/V+Re+y11tpi7nHi+B8YYwYaY9KNMemlaTqViu/rBZsZPH4Jm/cexgKb9x5m8PglSvxERM5RtWrVWLp0KYcPHwZgwoQJhIeHOxxV1eF00geuyRtFudNeeSL2T621z1pr0wrG7D0P9DPGnGjBM2dzD2vte9baBGttQnDwH3YykUpsxC8rOZybd0rZ4dw8Rvyy0qGIREQuvPPV43HZZZfxww8/APDZZ59x4403njy2e/du+vXrR+vWrenUqROLFy8GYNeuXVxyySXEx8dzzz33ULgd59NPP6Vjx460bduWe+65h7y8U//+lv9xOunbQ/EtbUEU3wJY2K6C9wlFyn8teG9b8L4bCDJ/7PgOKnRchGPH8/klcyub9x4u9viW05SLiFQ257PH44YbbuDzzz/nyJEjLF68mIsuuujkseeee474+HgWL17MK6+8wq233grA888/T5cuXViwYAF9+vRh40bXqm7Lly/niy++YObMmSxcuBBPT09SU1PPOcbKyunZu5m4xtwVFQcsc+Nc+GMr3onkLr9QPV+gMaeO6zvREljSfaSSW7ZlP+Mysvl64WZ2HzqGh4H8YtqGgwK8sdZqILSIVHjPf5fJsi37T3t8wca9HMvLP6XscG4ej49bzGfzil9GNy6sJs9dVdw/6adq3bo1WVlZfPbZZ1x++eWnHJsxYwZffulaQjclJYVdu3axb98+pk2bxvjx4wG44oorCApytdtMmjSJjIwMOnTo4Irx8GFCQkJKjKGqcjrp+xZ4zRjTyFq7DsAYEw0kAk+WcO5PuCZoXAp8X6i8d8F7esH7z8AxYACurt8TbgaWFswWlipm96FjfLNwM+Myssncsh9vT0PP5qFclxDB3kPHGPJ15ildvAbYnZPLNe/M4qnLm9MhutihoCIilULRhK+k8tLq06cPjz32GGlpaezatetk+R+H3/9vhmpxv3Bba7ntttsYNmxYmcRV2Tmd9I0G7ge+McY8javV7kVgE/DuiUrGmChgLa5dNF4AsNbuMsYMA54xxuzHtUhzAvAs8NGJZWCstduNMSOBwcaYA8B84E9ACiUvCyOVyPG8fNJW7mBcRjaTVmwjN8/SMrwmz/dpQZ82YQRV8zlZ18PDgxG/rGTL3sOE1fLn0V5Nyc23vP7rKq4bNZtL4kJ5/NJYmoQUXQdcRKT8K6lFLnH45GKHuoTX8ueLezqf8/3vvPNOAgMDadWqFWlpaSfLu3XrRmpqKs888wxpaWnUrVuXmjVrnix/+umn+emnn9izxzUCrEePHvTt25eHH36YkJAQdu/ezYEDB4iKijrnGCsjR5M+a+0hY0wKrm3YPsHVoDIJ1zZsBwtVNYAnfxyD+AJwAPgrrgWYfwdG4EocCxsCHAQe5H/bsF1vrf2uTB9IyqVV2w4wNn0TXy3Yws6DR6lTzYdbO0fTv30EzevXLPacfvHh9Iv/44yyq9qE8cGM9Yyauo7e/5jGDR0a8GDPpoTU8DvfjyEicsEM6h3D4PFLTunx8Pf2ZFDvmDK5fkREBA8++OAfyocOHcodd9xB69atCQgI4KOPPgJcY/1uvPFG2rVrR/fu3YmMjAQgLi6Ol156iUsuuYT8/Hy8vb3517/+paTvNExxTalyqoSEBJuenl5yRSk39uYc47tFWxiXkc2i7H14eRhSYkPo3z6C5NgQvD3PbQ7TzoNH+eek1aTO3YiPlwcDuzXi7q6NqObrdOO5iEjxli9fTvPmzd2u//WCzaf0eAzqHVPsL8PirOK+V2NMhrU2oWhdJX1uUNJXMeTlW6atdnXfTsjcxrG8fGLr1eC6hAb0bRtG3eq+ZX7P9TsPMeKXFfy4ZCt1q/vycK+m/CmhAV7nmFSKiJS10iZ9UjGUJulTs4RUeGu2H2RcRjZfLchm2/6j1Arw5qaLIunfPoIWYTXP62zbhnWr8e8B7cnYsIdhPy5nyFdL+WDGep64NJZecaGa6SsiIuWGkj6pkPYfyeX7Rb8zNmMTCzbuxdPDkNQsmKFXRZDSPARfL88LGk/7qCDG3tuZCcu2MfznFQz8JIMO0UEMvrw57SJL2lxGRETk/FPSJxVGXr5l1tqdjE3P5pfMrRw9nk/TkOo8dXks/eLDHZ9MYYzhkhb1SIkN4Yv0TYycsJpr/j2Ly1vVY1DvWBrWreZofCIiUrUp6ZNyL2vnIcZlZDN+fjZb9h2hpp8X1yc0oH/7CFpHBJa7LlQvTw8GXBRFv7bhvDdtHaOnr+PXzG0MuCiSB3o0PS9jC0VEREqipE/KpYNHj/PDYtfs29+y9uBhoGvTYJ66ojk9m4fi531hu2/PRjVfLx7u1YwBF0Xyj0mr+XTuRr6cv5l7uzfiri6N8Pcp/88gIiKVh6YYSrmRX9B9+8gXC+nw0kSe+HIJuw4d4/FLY5j1ZA8+urMjV7YOqxAJX2EhNf145epW/PJQVzo3rsNrv64i6bUpfPHbRvKK2+9NRMQpQ4eCMSW/hg49q8sbY3j00UdPfn7ttdcYWsK1srKyaNmyJQBpaWlceeWVf6iTlpZGYGAg8fHxxMTE0K1bN77//n+bdY0aNYqPP/74tPdIS0tj1qxZpz3+7bffMnz4cABuv/12xo0bd8aYi3rllVdO+XzxxReX6vyyopY+cdym3TmMy8jmy/nZZO85TA1fL/rFh9O/fQTtImuVu+7bs9UkpAajb01g3vrdvPLjcp74cgnvz1jPk5fFkhwTUmmeU0QqsKFDT03okpJc74V2zTgXvr6+jB8/nsGDB1O3bt0yueYJXbt2PZnoLVy4kH79+uHv70+PHj249957z3huWloa1atXLzYZO378OH369KFPnz5nHdsrr7zCU089dfLzmRLM80ktfeKInGPHGZeRzQ3vzabrq1N4a/JqGtatxps3tOW3p3sy7JpWtI8KqpSJUMeGtfnqrxfz7wHtOHY8nzvHpHPj6Dkszt7rdGgiIv+Tmgpz5sDUqRAd7fp8jry8vBg4cCAjR478w7GiLWjVq5/9Npdt27bl2Wef5e233wZcO3289tprALz11lvExcXRunVrbrjhBrKyshg1ahQjR46kbdu2TJ8+ndtvv51HHnmE5ORknnjiCcaMGcP9999/8voTJ06ka9euNGvW7GSiWbTOlVdeSVpaGk8++SSHDx+mbdu2DBgw4JRns9YyaNAgWrZsSatWrfjiiy8AVxKalJRE//79iY2NZcCAAcXuS1xaaumTC8Zay29Zexibvokfl/zOoWN5RNcJ4LFLmnF1uwjCa/k7HeIFY4zh8lb16dk8lM/mbeTNSavp8/ZMrmoTxqBLYoisE+B0iCJSlaWmwsCBcPSo6/OGDa7PAAWJy9m67777aN26NY8//vg5Bnlm7dq1Y8SIEX8oHz58OOvXr8fX15e9e/dSq1Yt7r33XqpXr85jjz0GwPvvv8+qVauYOHEinp6ejBkz5pRrZGVlMXXqVNauXUtycjJr1qw5bRzDhw/n7bffZuHChX84Nn78eBYuXMiiRYvYuXMnHTp0oFu3bgAsWLCAzMxMwsLCSExMZObMmXTp0uUc/oso6ZMLYPPew4zPyGbc/Gw27Mqhmo8nV7Suz3UJDUiopK157vLx8uC2i6O5pl04705dx39mrOPnpb9zS6doHkhpQlA1H6dDFJGqaMgQyMk5tSwnx1V+jklfzZo1ufXWW3nrrbfw9z9/v+yfrmWsdevWDBgwgH79+tGvX7/Tnn/dddfh6Vn8GPLrr78eDw8PmjZtSqNGjVixYsVZxThjxgxuvPFGPD09CQ0NpXv37vz222/UrFmTjh07EhERAbhaLrOyspT0Sfl0+Fgev2RuZVxGNjPX7sRa6NyoDg/2aMqlLesR4KM/eoXV8PPmsd4x3NwpipETVjFm1nrGZmzir0lNuCMxusJNXhGRCm7jxtKVl9JDDz1Eu3btuOOOO06WeXl5kZ+fD7gStmPHjp3TPRYsWFDstnM//PAD06ZN49tvv+XFF18kMzOz2POrVTv92qpFGyuMMafED3DkyJESYzxTl62v7/+W9/L09OT48eMlXq8kGtMnZcZaS8aG3Qwev5iOL0/koS8WkrXrEA/2aMr0x5P5bGAnrmkXoYTvDOoF+vH3/q356cFudIiuzd9/XkHya2mMy8jWTF8RuXAiI0tXXkq1a9fm+uuv5/333z9ZFh0dTUZGBgDffPMNubm5Z339xYsX8+KLL3LfffedUp6fn8+mTZtITk7m1VdfZe/evRw8eJAaNWpw4MABt68/duxY8vPzWbt2LevWrSMmJobo6GgWLlx48h7z5s07Wd/b27vY5+nWrRtffPEFeXl57Nixg2nTptGxY8ezfu6S6F9fOWdb9x3hy/nZfJmRzbqdh/D39uTyVvXp3z6CixrWxsOj6nbfnq2YejX44PYOzFq7k+E/reCxsYv4z/R1DL68Od2bBTsd3oUzdCg8/3zJ9Z577qyXkBCRYrz8smsMX+Eu3oAAV3kZefTRR09OtAC4++676du3Lx07dqRHjx5nbGkrzvTp04mPjycnJ4eQkBDeeustevTocUqdvLw8br75Zvbt24e1locffphatWpx1VVX0b9/f7755hv++c9/lnivmJgYunfvzrZt2xg1ahR+fn4kJibSsGFDWrVqRcuWLWnXrt3J+gMHDqR169a0a9eO1EITYq6++mpmz55NmzZtMMbw6quvUq9evbPuLi6JKYvZIJVdQkKCTU9PdzqMcuVIbh4Tlm1jbEY2M1bvIN9Cx+ja9E+I4PJW9anuq98nykp+vuX7Jb8z4pcVbNp9mK5N6/LEpbG0DA90OrQLr4yXjxCpSpYvX15sd+dppabCXXe5JnNERbkSvnMczydlr7jv1RiTYa1NKFpX/zKL26y1LMrex7iMTXy7cAv7jxwnLNCP+5KbcG27CKIr+t6y5bRVycPD0KdNGL1bhPLpnI38c/Jqrnp7Ble3DeeRS5oREaSZviJyHgwYAKNHu37WL1qVgpI+KdH2/Uf4asFmxmVks3r7QXy9PLisZT2uS2hA50Z1Kk/37XlelPRc+Xp5cleXhvRvH8E7aWv5YOZ6vl/yO3dcHM1fk5oQGODtdIgiUtGd7pffoqssaEhFhaSkT4p19Hgek5ZvZ1xGNlNX7SAv39I+Kohh17Tiitb1qemnBMMpgf7ePHlZLLd0juKNX1fx3vR1fP7bJh5IacItnaPw9dJMXxE5S0V/+ZVKRUmfnGStJXPLfsamb+KbRVvYm5NLvZp+3NOtEde2j6Bx8Nmvji5lL7yWP69f34a7ujRk+M8reOmH5YyZlcWg3jFc1Tqs8rTAikiZsdZW6bVRK5vSzstQ0ifsPHiUrwu6b1dsPYCPlweXxIVyXUIDujSpi6eSh3ItLqwmH9/ZkemrdzDsxxU8+PlCRk9fx1OXNefiJmW7t6WIVFx+fn7s2rWLOnXqKPGrBKy17Nq1Cz8/P7fPUdJXReXm5TN5hav7dsqK7RzPt7RpUIsX+7WkT+swjQ+rgLo2DSbxgbp8s2gzr/2yipv+M5ekmGCevCyW2Ho1nQ5PRBwWERFBdnY2O3bscDoUKSN+fn4nd+1wh5K+Kmb57/sZm57NNws3s+vQMYJr+J6cHNA0tIbT4ck58vAwXB0fwWUt6/Px7CzenryGy96czrXtInj0kmbUD6w6+xuLyKm8vb1p2LCh02GIg5T0VQG7Dx3j24WbGZuRTeaW/Xh7GnrFhdK/fQTdmgbj5amNWSobP29PBnZrzPUJDfjXlDV8NGsD3y3awp1dGvKXpMaaiCMiUgUp6aukjuflM3XVDsamZzNpxTZy8yytwgN5vk8L+rQJI6iaj9MhygVQK8CHIVfEcWvnaF7/dSXvpK3l83kb+VuPpgy4KAofLyX8UkrldD1LESmZduRwQ0XakWPVtgOMy8hm/PzN7Dx4lDrVfOgXH07/9hE0r69xXaVSztbpKwtLsvcx7KflzFq7i8jaATx+aQxXtKpfcQZ1V8LvpMLTdyJS7mhHjkpsX04u3y5yzb5dlL0PLw9DSmwI/dtHkBwbgre6b6VAq4hAUv98EVNX7WD4Tyu4//8WMDrCtadvp0Z1nA5PRETOIyV9FVRevmX66h2MzchmQuY2juXlE1uvBs9cGUfftmHUre7rdIhSThljSIoJoWvTYMbPz+b1X1dxw3tz6Nk8hCcujdWEHhGRSkpJXwWzdsfBgu7bbLbtP0pQgDc3XRRJ//YRtAwPdDo8qUA8PQzXJTTgytZhfDhrPe9MWUvvf0zjTx0a8FDPZoTWdH/tJxERKf+U9FUA+4/k8v2i3xmXsYn5G/fi6WFIahbM831c3bfadkvOhb+PJ39NasINHSL55+TVfDpnA18v2MLdXRsysHtjqvvqrwkRkcpAf5uXU3n5lllrdzIuI5ufl27l6PF8moZU56nLY+kXH05IDbXCSNmqXc2H565qwe0XR/PqLyt5a/IaUudu5KGeTbmhY6TGhoqIVHCOJ33GmAbASKAXYICJwEPW2o1unHu6qcfx1tqFheplAVHF1LvaWvt1qYMuQ1+//V9GrDzCloAgwnL2cHtjf/Y1b834+dls2XeEmn5eXJ/QgP7tI2gdEVhxZllKhRVVpxr/uqkdd3fdyys/LueZbzL5YGYWT1waQ+8W9fRnUESkgnI06TPGBACTgaPAbYAFXgKmGGNaW2sPuXGZMcC7RcpWFVPvF2BokbKVpYm3rH399n8ZvN6Tw9VcsyY3V6vDy1uBratJignhqSua07N5KH7e6r6VC69tg1p8MbATk1dsZ9hPK7j30/m0i6zFU5c3JyG6ttPhiYhIKTnd0nc30AiIsdauATDGLAZWA/cAb7hxjc3W2jlu1NvpZr0LZsTKIycTvsLq5exlzB1XOhCRyKmMMfRoHkr3ZsGMy8jmjQmr6D9qNr1bhPL4pbE0Dq7udIgiIuImpwfp9AHmnEj4AKy164GZQF/HojrfkpLAGLYEBBV7eJt/IBjzv0VPRRzm5enBDR0jSRuUxKO9mjFj9U4uGTmNp79ewo4DR50OT0RE3OB00tcCWFpMeSYQ5+Y1/mKMOWqMyTHGTDbGdD1NvasK6hw1xswxxvQ7q4jLQloaWEtYzp5iD4fl7AFrtcK9lDsBPl480KMpUx9PZsBFkXw+bxNJI6bw5sTVHDp63OnwRETkDJxO+moDxWU+u4Him8FO9SnwV6AnMBCoA0w2xiQVqfcd8ADQGxgAHAG+MsbcfLoLG2MGGmPSjTHpO3bscCOU0hsU44d/7pFTyvxzjzAoRjNzpXyrW92XF/q25NeHu9GtWTAjJ64i6bU0/m/uRo7n5TsdnoiIFMPppA9ckzeKcmt6oLX2FmvtF9ba6dbaT4EuwBZck0EK13vAWvtxQb1xQA8gHRh2hmu/Z61NsNYmBAcHu/0wpdHv/usZ1jCP8EO7MDaf8EO7GNYwj373X39e7idS1hoFV+edm9vz5V86E1k7gKe+WkLvf0xjwrJtaF9vEZHyxemJHHtwtfYVFUTxLYBnZK09YIz5AbirhHp5xpixwN+NMfWttb+X9l5lpd/91+NcP7NI2WgfVZtx93bm12Xb+PtPK7j743Q6Rtdm8OWxxEe602gvIiLnm9MtfZm4xvUVFQcsO8trGopvPSyuHm7WFZESGGPo3aIevzzcjRf7tWTdzoNc/e9Z3Jc6n6yd7qy+JCIi55PTSd+3QCdjTKMTBcaYaCCx4FipGGNqAlcAc0uo5wVcB2y01m4t7X1E5PS8PT24pVMUaYOS+VuPpkxesZ2eb0xl6LeZ7Dqomb4iIk5xOukbDWQB3xhj+hpj+gDfAJsotOCyMSbKGHPcGPNsobLHjDGjjTE3GWOSjDG34VrqpR7wdKF6NxpjPjfG3GqMSTbG3ABMAdoDT1yIhxSpiqr7evFIr2ZMHZTEdQkN+Hh2Ft1HpPGvKWs4fCzP6fBERKocR5O+gh03UnDtoPEJkAqsB1KstQcLVTWAJ6fGuxJXN/BbwARcCzmvB7pYa6cXqrceCAFGAL/iSiaPApdaaz8/D48lIoWE1PRj2DWt+PXhbnRqVIcRv6wk6bUp/Pe3TeTla3SFiMiF4vREDgr22L22hDpZFJnRa639DtdSLCVdfw6uxFJEHNQkpAb/uS2Buet28cpPK3j8y8W8P2M9T14WS1JMsPb0FRE5z5zu3hWRKuaiRnX4+q8X86+b2nHkeB53jPmNm0bPZXH2XqdDExGp1JT0icgFZ4zhitb1mfBwd4ZeFcfKbQfo8/ZM/vbZAjbtznE6PBGRSklJn4g4xsfLg9sTG5I2KIn7khvzS+ZWerw+lRe/X8aeQ8ecDk+k4hg61LVne0mvoUOdjlQcpKRPRBxX08+bQb1jSRuURL/4MD6YuZ5uI6YwaupajuRqpq9IiYYOde3ZfuLVvbvrVbjMWiV9VZySPhEpN+oH+vNq/zb89GBXEqKCGP7TClJeS+PLjGy+ysgmMX4gDS96jMThk/l6wWanwxURqVAcn70rIlJUbL2afHhHR2at2cmwn1bw6NhFrq12fAMB2Lz3MIPHLwGgX3y4g5GKiFQcaukTkXLr4iZ1+ea+RIICvP+wX+Lh3DxG/LLCkbhERCoiJX0iUq55eBj25uQWe2zz3iM8/10ms9bsJDcv/wJHJiJSsah7V0TKvbBa/mzee/gP5b5eHqTO3ciHM7Oo4edFckwIPeNC6d4smEB/bwciFREpv5T0iUi5N6h3DIPHL+FwoZm8/t6eDLumFZe0CGX66p1MXLaNySu28+2iLXh5GC5qVJuezUPp2TyUBrUDHIxeRKR8UNInIuXeickaIz6eyhafmoQFBTCod8zJ8t4t6tG7RT3y8i0LN+1hwrLtTFy+jee/W8bz3y0jtl4NVwIYF0rr8EA8PLTlm4hUPUr6RKRC6BcfTr+H33N9SEsrto6nh6F9VG3aR9XmyctiWb/zEJOWb2PCsm38O20Nb09ZQ3ANX3o2D6Fn81ASm9TFz9vzwj2EiIiDlPSJSKXVsG41/ty1EX/u2og9h46Rtmo7E5dt59uFW/hs3ib8vD3o2jSYXs1DSY4NIbiGr9Mhi4icN0r6RKRKCKrmw9XxEVwdH8HR43nMXbebicu3MXGZqyXQGIhvUIuecaH0ah5Kk5DqGKNuYBGpPJT0iUiV4+vlSbdmwXRrFszzfVqw7Pf9TCwYB/jqzyt59eeVRNUJODkRpEN0EF6eWuFKRCo2JX0iUqUZY2gRFkiLsEAe7NmU3/cdZtJyVwL4yewNvD9jPYH+3iTHBJ9cDqaGn5aDEZGKR0mfiEgh9QP9ublTFDd3iuLg0ePMWL2DCcu2M3nFNr5euAVvT0OnRnXoFRdKj+ahhNfydzpkERG3KOkTETmN6r5eXNqyPpe2rE9evmX+xj2uMYDLt/HsN5k8+00mcfVrnhwH2DK8psYBiki5paRPRMQNnh6GDtG16RBdm8GXN2ftjoNMWr6Nicu28/bk1bw1aTWhNX1PrgfYuVEdLQcjIuWKkj4RkbPQOLg6jYOrM7BbY3YfOsaUFa5xgF8t2Ezq3I0E+HjSrWkwPZqHkBIbQp3qWg5GRJylpE+kOKmpMGcOHD0K0dHw8sswYIDTUUk5VbuaD9e2j+Da9hEcyc1jzrpdBcvBbOfnzK0YA+0jg+gZ55oN3CSkutMhi0gVpKRPpKjUVBg40JXwAWzY4PoMSvykRH7eniTFhJAUE8KLfS2ZW/YzYdk2Ji7fxvCfVjD8pxU0qlvtZALYLrKWloMRkQtCSZ9IUUOGQE7OqWU5Oa5yJX1SCsYYWoYH0jI8kId7NWPz3sNMXr6NX5dt48OZ63lv2jqCArxJjg2hV/NQujYLprqv/loWkfNDf7uIFLVxY+nKRdwUXsufWzpHc0vnaA4cyWXaqp1MXL6NScu3M37+Znw8PejcuE5BK2AI9QO1HIyIlB0lfSJFRUa6unSLKxcpIzX8vLmidX2uaF2f43n5pG/433Iwz3y9lGe+hpbhNU/uCtIiTMvBiMi5UdInUtTLL7vG8BXu4g0IcJWLnAdenobjzJ8AACAASURBVB50alSHTo3qMOQK13IwEwq2hXtz0mr+MXE19QP9Ti4H06lRbXy9tByMiJSOkj6Rok6M27vrLtdkjqgozd6VC8YYQ5OQGjQJqcFfkhqz8+BRJq/YzsRl2xiXkc0nczZQzceT7jHB9GweSnJMCEHVfJwOW0QqACV9IsUZMABGj3b9nJbmaChStdWt7sv1CQ24PqEBR3LzmLV2JxOWbWfS8m38uGQrHgYSomvTq6AVsGHdak6HLCLllJI+EZEKws/bk5TYUFJiQ8nPb8mSzfuYuHwbE5Zt4+Ufl/Pyj8tpHFzt5LZw8ZFBeHpoHKCIuDie9BljGgAjgV6AASYCD1lrS5wqaYyxpzkUb61dWKieB/AEcA9QD1gJvGCt/fIcwxcRcYSHh6FNg1q0aVCLRy+JYdPuHNe2cMu38/709bw7dR21q/mQEhtCz+ahdG1al2paDkakSnP0bwBjTAAwGTgK3AZY4CVgijGmtbX2kBuXGQO8W6RsVZHPLwKPAUOADOAGYKwx5kpr7Y9n/wQiIuVDg9oB3J7YkNsTG7L/SC5TV+5g4vJt/Jq5lXEZ2fh4eZB4cjmYUEJr+jkdskjlNnQoPP98yfWee85V9wJw+te+u4FGQIy1dg2AMWYxsBpXq9wbblxjs7V2zukOGmNCcCV8w621rxUUTzHGNAGGA0r6RKRSqennzVVtwriqTRi5efn8lrWbicu2M2H5VqZ8tYMhXy2ldUTgyeVgmtevoeVgRMra0KGnJnNJSa53B8eJO5309QHmnEj4AKy1640xM4G+uJf0laQ34AN8WqT8U+ADY0xDa+36MriPiEi54+3pwcWN63Jx47o8c2VzVm8/eHJbuJETV/HGhFWE1/KnZ/MQesaFclHDOvh4aVs4kcrI6aSvBfBNMeWZwHVuXuMvxphBQB4wB3jOWju9yD2OAmuKnJdZ8B4HKOkTkUrPGEOz0Bo0C63BfclN2H7gCFNWbGfCsu18kb6Jj2ZvoIavF91igunVPJSkmGBqBWg5GJHKwumkrzawp5jy3UCQG+d/CnwPbAGigEHAZGNML2ttWqF77LXWFp30sbvQcRGRKiekhh9/6hDJnzpEcvhYHjPXuLaFm7h8Oz8s/h1PD0OH6CB6Ng+lV1woUXW0HIxIReZ00geuyRtFuTW4xFp7S6GP040x3wBLcU0G6VLoWqW+hzFmIDAQIFLbb4lIJefv4+ma5BEXSn6+ZVH2XlcCuGw7L/2wnJd+WE7TkOonJ4LEN6jFt4u2MCJ+IFt8ahI2fDKDesfQLz7c6UcRkdNwOunbQ/EtbUEU3wJ4RtbaA8aYH4C7ChXvBoKMMaZIa19QoePFXes94D2AhISE0y0NIyJS6Xh4GOIjg4iPDGJQ71g27sopaAHcxnvT1vFO2lqq+3pyODefPN9AADbvPczg8UsAlPiJlFNOj9bNxDXmrqg4YNlZXrNoy14m4As0LuYenMN9RESqhMg6AdzZpSH/d3cn5j/dizdvaEtePuTln/r78OHcPEb8stKhKEWkJE4nfd8CnYwxjU4UGGOigcSCY6VijKkJXAHMLVT8M3AMKLpx6s3AUs3cFakgUlNhzhyYOhWio12f5YILDPCmb9twjuTmFXt8897DLP99/wWOSkTc4XTSNxrIAr4xxvQ1xvTBNZt3E4UWXDbGRBljjhtjni1U9pgxZrQx5iZjTJIx5jZgJq4dN54+Uc9aux3Xjh+DjTGPFNR9B0gBnroAzygi5yo1FQYOhKNHXZ83bHB9VuLnmLBa/sWWG+CyN6dz+4fzmLtuF3+cQyciTnE06SvYcSMF1w4anwCpuJZPSbHWHixU1QCenBrvSlxdtG8BE3Ct6bce6FJkyRZw7cTxEvAg8AuulsTrrbXflfUzich5MGQI5OScWpaT4yoXRwzqHYO/t+cpZf7enrx8dUse7dWMJdn7+NN7c7jmnVn8krmV/HwlfyJOc3oiBwV77F5bQp0sisy2LUjY3ErarLV5uJK+l84uShFx1MbTbMV9unI5705M1hjx8VTX7N2ggFNm797drRFj0zfx3vR13PNJBo2Dq3FPt8b0jQ/D18vzTJcWkfPE8aRPRKREkZGuLt3iysUx/eLD6ffwe64PRbaW8vP25JbO0dzYMZIfl25lVNpaHv9yMa9PWMldXRpyY8dIavh5X/igRaowp8f0iYiU7OWXISDg1LKAAFe5lGtenh70aRPGD3/rwsd3dqRxcHVe+XEFFw+fzKs/r2DHgaNOhyhSZailT0TKvwEFk+/vuss1mSMqypXwDSg6KV/KK2MM3ZoF061ZMIs27eXdaWt5Z+pa/jNjPf3bRzCwayOi62rHD5HzSUmfiFQMAwbA6NGun4t0JUrF0qZBLf49oD3rdx7ivWnrGJeezefzNnJZy/rc270xrSICnQ5RpFJS0iciIo5oWLcaw65pxcO9mvLhzCw+nb2BH5b8TmKTOtzbvTFdmtTFGLd25RQRN2hMn4iIOCqkhh9PXBrLrMEpDL4sltXbDnLL+/O46u0ZfLdoC8fz8p0OUaRSUNInIiLlQg0/b+7p3pjpTyTz92tbkXMsjwc+W0DK61P5ZM6G0+4CIiLuUdInIiLliq+XJ3/qEMnEh7sz6ub2BFXz4Zmvl9Ll75N5e/Jq9uXkOh2iSIWkMX0iIlIueXgYLm1Zj94tQpmzbjejpq7ltV9X8U7aWm7sGMldXRtSP7D47eBE5I+U9ImISLlmjKFz4zp0blyHZVv28+60tXw4K4uPZmfRt20493ZvRJOQGk6HKVLuqXtXREQqjLiwmrx5QzxpjyVxU8dIvl+8hZ5vTOPPH6WTsWG30+GJlGtK+kREpMJpUDuA5/u2ZOYTKfytR1PSN+zm2ndmc92oWUxesQ1rrdMhipQ7biV9xhhfY8w8Y8wkY8xpN0s0xvgU1JlzpnoiIiJloU51Xx7p1YyZT6Tw7JVxbN5zmDvHpHPpP6Yzfn42uVruReQkd1v6BgDtgdettaedNmWtPQaMADoWnCMiInLeVfP14s4uDZn6eDJvXN8GgEf+u4ikEWl8MGM9OceOOxyhiPPcTfquAdZZa38sqaK19mdgNXDduQQmIiJSWt6eHlzTLoKfH+rKB7cnEF7Lnxe+X8bFwyfzxoRV7D50zOkQRRzj7uzdeKDEhK+QacDlpQ9HRETk3BljSIkNJSU2lIwNexg1dS1vTVrNe9PW8qeEBvy5ayMa1A5wOkyRC8rdpK8usK0U190G1Cl9OCIiImWrfVQQo29NYM32A7w7dR3/N28jn87dyJWt63NPt8bEhdV0OkSRC8Ld7t3DQPVSXLc6cKT04YiIiJwfTUJqMOK6Nkx7PJk7E6OZuGwbl781nds+mMfstbs041cqPXeTvk1Ah1JcNwHYWPpwRESkwkhNhTlzYOpUiI52fa4A6gf6M+SKOGY92YNBvWPI3LKPG0fPod+/Z/Hz0t/Jz1fyJ5WTu0lfGtDJGJNQUkVjTHvgYmDKOcQlIiLlWWoqDBwIR4+6Pm/Y4PpcQRI/gMAAb+5LbsKMJ1J4sV9L9hw6xr2fzqfnG1P5fN5Gjh7PczpEkTLlbtL3NmCBscaY5qerZIyJBcYCecC/zz08EREpl4YMgZycU8tyclzlFYyftye3dIpiymNJvH1TPP4+njw5fgld/z6FUVPXsv/IaVcqE6lQ3JrIYa1daYx5ARgKLDDGjAMmA9m4ksEIoAdwLeALPGutXXleIhYREedtPM0IntOVVwCeHoYrW4dxRav6zFizk1FT1zL8pxX8a/IaBnSK4s7EaEJq+jkdpshZc3f2LtbaF4wxx4HngJuAG4tUMUAuMMRaO6zsQhQRkXInMtLVpVtceQVnjKFr02C6Ng1mSfY+Rk1dy3vT1vLBjPVc2z6cgd0a07BuNafDFCk1t5M+AGvtK8aYVOBOIBGojyvZ2wLMAD601hbzt4CIiFQqL7/sGsNXuIs3IMBVXom0igjkXwPakbXzEO9NX8e4jGw+/20Tl7aox73dG9OmQS2nQxRxW6mSPoCCpO658xCLiIhUFAMKdtq86y7XZI6oKFfCN6By7sAZXbcar1zdiod6NmXMzCw+mbOBn5ZupXOjOtyb1JhuTetijHE6TJEzKnXSJyIiArgSvNGjXT+npTkayoUSUsOPxy+N5S9Jjfls3kben7Ge2z6YR1z9mtyb1JjLW9bDy9PdOZIiF5ZbfzKNMeuMMX8rUhZpjOl2fsISEREpv2r4eTOwW2OmPZ7Mq9e25sjxPP722QKSX0/jk9lZHMnVci9S/rj760g0UHTgwh1oLT4REanCfL08ub5DAyY+3J13b2lP3eq+PPNNJonDJ/PWpNXszTnmdIgiJ6l7V0RE5Bx5eBh6t6jHJXGhzFu/m1FT1/LGhFWMmrqWGztGcleXhoTV8nc6TKnilPSJiIiUEWMMFzWqw0WN6rBi637enbqOMbOy+GhWFn3bhnNv90Y0Da3hdJhSRTk+2tQY08AYM84Ys88Ys98YM94YU+qFnowxg40x1hgzo5hjWQXHir76lc1TiIiInCq2Xk1G/qktUwclcXOnKH5c8ju9Rk7jzx/9RnrWbqfDkyrI0ZY+Y0wArp09jgK34drd4yVgijGmtbX2kJvXaQQMAbafodovuHYUKUy7hoiIyHkVERTA0D4t+FuPpnw8O4sxs7LoP2o2CVFB3Nu9MSmxIXh4aLkXOf9Kk/QlFVmDKAnAGPMMrgWai7LW2hdLuObdQCMgxlq7puB6i4HVwD3AG27G9g6QCsRw+mfaaa2d4+b1pCoaOhSef/6P5UXX3nruOVddEZFSqF3Nh4d6NmNgt0Z88dsm/jN9PX/+OJ2mIdW5p3tj+rQJw8fL8Q44qcSMtbbkSsbkl+KaFlcSaK21niVcdxLgZ61NLFI+FdcFursR203Am7gSvvGAl7W2S5E6WcAMa+3NpXiOkxISEmx6evrZnCoiZSkpyfVeRdaEqxD0nZy13Lx8flj8O6OmrmXF1gPUD/Tjri4NubFjJNV8z7EjTt9L+XMBvxNjTIa1NqFoubt/qopp/igTLYBviinPBK4r6WRjTBAwEnjcWru7hNXQrzLG5ACewAJguLX269KHLCIicu68PT3oFx9O37ZhpK3awTtpa3nph+X8c/Iabu0cxW0XR1O3uq/TYUol4lbSZ609X0lfbWBPMeW7gSA3zh8BrALGlFDvO+A3YD0QCtwPfGWMucVa+6nb0YqIiJQxYwzJMSEkx4Qwf+MeRqWt5e0pa3hv2jquT2jA3V0bEVknwOkwpRIoD0u2FNe/XOKIVmNMV+BWoJ0toY/aWvtAkXO/AuYAw4Bikz5jzEBgIEBkZKknE4uIiJRau8gg3rs1gTXbD/LetLV8/ttGUudu4IrWYdzbvREtwgKdDlEqMLdHjBpjrjTG3GKM8T5DHZ+COle4edk9uFr7igqi+BbAwt4F3geyjTG1jDG1cCWxngWfT9smbq3NA8YCEcaY+qep8561NsFamxAcHOzOs4iIiJSJJiHVebV/G6Y/nsKfuzZiyortXPHWDG55fy6z1uzEnfH4IkW5u/duK+BroKO1Nvd09ay1x4AE4GtjTAs3Lp2Ja1xfUXHAshLObQ7ciys5PPFKBDoV/PyXEs4/0Zqo/3NERKRcqhfox1OXN2fmkykM6h3D8t8PcNN/5tL3XzP5ccnv5OXrnzBxn7stfXcCx3BvQseLBXX/7Ebdb4FOBevsAWCMicaVvH1bwrnJxbwWAUsLfh53uhONMV64JopstNZudSNOERERxwT6e3NfchNmPJHMy1e3ZN/hXP6aOp+eb0zl/+Zu5EhuntMhSgXg7pi+JCDNWruzpIrW2p3GmDQgxY3rjsY1qeIbY8zTuFrdXgQ24eq+BcAYEwWsBV6w1r5QcJ+0ohczxuzFtWRLWqGyG4G+wI8F1w0F7gPaAze6EaOIiEi54OftyYCLorihQyQ/L93KqKlreeqrJYycuIo7EqO5uVMUk//zFSPibmNLjbqEPfAxg2L86Hf/9U6HLuWAu0lfY2BSKa67HOhSUiVr7SFjTAquZVc+wdXlOgl4yFp7sFBVg2uplbNZtXI9EIJrpm9tIAfXTN5LrbW/nMX1REREHOXpYbiidX0ub1WPWWt3MWrqWl79eSVv/ryCvDxfjtesBsDmanUYvP4IvP1fJX7idtLng6vL1l3HALcWF7LWbgSuLaFOFm7M6LXWJhVTNgf3Wh1FREQqFGMMiU3qktikLks37+O6kZM46nXqP7+Hvf0YsXIX2mxe3G052w2UZt2SSGBX6cMRERGRs9EyPJAjnsUvsLE5oDbH80qzuZZURu4mffOBXmdaBuUEY4wf0AvXrhciIiJygYTlnGa1M2PoPiKNd6euZV/OaRfhkErO3aTvKyAYeNmNui8AdYEvzzYoERERKb1BMX745x45pcwv9yh31D5CRJA/w35aQadhkxjy1RJWbzvgUJTiFHfH9H0MPAo8bIwJxTWLdnXhCsaYpsAzwABcEzk+KctARURE5Mz63X89vP1fRizb7pq9m7PnlNm7mVv2MWZmFmMzskmdu5GuTetyR2I0Sc1C8PAocei8VHDG3VW9jTHNcM2sDce1tMpmILvg54iClykoS7HWrjkfATshISHBpqenOx2GiCQlud7T0pyMQgrTd1I+lfC97Dp4lM/mbeSTORvYtv8o0XUCuO3iaPq3j6CG32k33pJzcQH/XzHGZFhrE4qWu70EirV2Fa617d7HNTs3AtfuF52BBgVl/wESKlPCJyIiUtnUqe7L/SlNmfFECm/e0Jagaj48/90yOg+bzPPfZZK185DTIcp54G73LgDW2u3A3caYB3Btt1YfV+veFiDdWnvkTOeLiIhI+eHt6UHftuH0bRvOwk17+XDmej6ZvYExs7JIiQnhjsSGJDapgzHq+q0MSpX0nVCQ3M0o41hERETEIW0b1OLNG+J56vLmpM7ZQOrcjdz8/lyahlTn9sRoromPwN/H0+kw5Ry41b1rjLnVGNP6fAcjIiIizgqt6ccjl8Qw88kUXruuDT5eHgz5aimdhk1i2E/Lyd6T43SIcpbcHdM3Bk5dzNsYc5sxZnKZRyQiIiKO8/P2pH/7CL5/oAv/vacziU3qMHraOrq9OoW/fJrB3HW7cHcyqJQPZ9W9WyAa6F5GcYiIiEg5ZIyhY8PadGxYm817D/Px7Cw+n7eJn5ZuJa5+Te5IjOaqNmH4eavrt7xze/auiIiIVG3htfwZfFlz5gzuwStXt+J4fj6Dxi0mcfhkXv91Jdv2az5neXYuLX0iIiJSBfn7eHLTRZHc2LEBs9bu4sOZ63l7yhreSVvL5a3qc0diNPGRQU6HKUUo6RMREZGzYowhsUldEpvUJWvnIT6evYH/pm/i20VbaNugFnckRnNZy/r4eKljsTwozbeg0ZoiIiJSrOi61Xj2qjjmPNWDoVfFse9wLg9+vpAuf5/MPyetZtfBo06HWOWVpqVvqDFmaNFCY0zeaepba61aEkVERKqQ6r5e3J7YkFs7RzN11Q4+mLme1yes4p9T1tCnTRh3JEbTIizQ6TCrpNIkZaVdjlvLd4uIiFRRHh6G5NgQkmNDWLP9AGNmZfFlxmbGZWTTsWFt7rg4ml5xoXh5quv3QnEr6bPW6hsRERGRs9IkpAYv9WvFoEti+SJ9Ix/N2sBfUucTXsufWztHcUOHSAIDvJ0Os9JTMiciIiIXRGCANwO7NWba48mMurk9EUH+DPtpBZ2GTeKpr5awetsBp0Os1DTmTkRERC4oTw/DpS3rcWnLemRu2cdHs7IYl5HN/83dSNemdbn94miSY0Lw8NBIsbKklj4RERFxTIuwQF7t34bZT6bw2CXNWLXtAHd9lE7K62l8OHM9B47kOh1ipaGkT0RERBxXp7ov96c0ZcYTKbx1Yzy1q/nw/HfL6DxsMkO/zWT9zkNOh1jhqXtXREREyg1vTw/6tAmjT5swFm3ay4cz15M6dwMfzc4iOSaEOxKj6dKkLsao67e01NInIiIi5VKbBrX4xw3xzHwihQdSmrI4ey+3vD+PXiOn8emcDeQcO+50iBWKkj4REREp10Jq+vFIr2bMfDKF165rg6+XB09/vZROr0xi2I/Lyd6T43SIFYK6d0VERKRC8PXypH/7CK5tF076hj18OHM9/5mxntHT13FJXD3uSIymY8Pa6vo9DSV9IiIiUqEYY+gQXZsO0bXZvPcwn8zewGfzNvJz5lbi6tfk9sRo+rQJw8/b0+lQXVJTYc4cOHoUoqPh5ZdhwIALHoa6d0VERKTCCq/lz5OXxTJncA+GXdOK4/n5PD5uMYnDJ/P6ryvZtv+IswGmpsLAga6ED2DDBtfn1NQLHoqSPhEpv4YOBWP+95o61fUqXGaMq56IVGn+Pp7c2DGSXx7qRuqfLyI+shZvT1lD4vDJ/O2zBSzYuMeZwIYMgZwiYw5zclzlF5ix1l7wm1Y0CQkJNj093ekwRETKn6Qk13tampNRSFH6XgDYsOsQH83awNj0TRw4epw2DWpxZ2I0l7Wsj4/XBWr38vCA4nItYyA//7zc0hiTYa1N+EMo5+VupWCMaWCMGWeM2WeM2W+MGW+MiTyL6ww2xlhjzIxijnkUHM8yxhwxxiwyxlxbNk8gIiIi5VFUnWo8e1Ucs5/qwfN9WrD/cC4Pfr6QLn+fzFuTVrPz4NHzH0TkaVKa05WfR44mfcaYAGAyEAvcBtwCNAWmGGOqleI6jYAhwPbTVHkRGAq8DVwGzAHGGmMuP+vgRUREpEKo7uvFbRdHM+mR7nx4ewdi6tXgjQmruHjYZB797yKWbt53/m7+8ssQEHBqWUCAq/wCc3r27t1AIyDGWrsGwBizGFgN3AO84eZ13gFSgRiKPJMxJgR4DBhurX2toHiKMaYJMBz48VwfQkRERMo/Dw9DcmwIybEhrNl+gDGzsvgyYzNfzs+mY3Rt7kiMpldcKF6eZdgmdmKW7l13uSZzREVV2dm7fYA5JxI+AGvtemAm0NedCxhjbgLaAYNPU6U34AN8WqT8U6CVMaZhaYMWERGRiq1JSA1e6teKOYN7MOTy5mzZd5i/pM6n+4g0Rk1dy96cY2V3swEDoFMn6N4dsrIcSfjA+aSvBbC0mPJMIK6kk40xQcBI4HFr7e4z3OMosKZIeWbBe4n3ERERkcopMMCbu7s1YuqgZN69pT0Navsz/KcVdBo2icHjl7Bq2wGnQywzTnfv1gaKm0O9Gwhy4/wRwCpgTAn32Gv/OE15d6HjIiIiUoV5ehh6t6hH7xb1WLZlPx/NyuLL+dl8Nm8jXZrU5faLo0mJDcHDo+Lu9uF00gdQ3JoxJf4XNcZ0BW4F2hWT0BW9VqnvYYwZCAwEiHRgho2IiIg4Iy6sJn/v35onLovls3kb+WT2Bv78cTpRdQK4rXM01yVEUMPP2+kwS83p7t09FN/SFkTxLYCFvQu8D2QbY2oZY2rhSmI9Cz77FtTbDQSZP27EF1To+B9Ya9+z1iZYaxOCg4PdeRYRERGpRGpX8+G+5CZMfyKZf94YT51qPrzw/TI6vTKJod9msn7nIadDLBWnW/oycY25KyoOWFbCuc0LXvcWc2wP8DDwj4J7+AKNOXVc34mxfCXdR0RERKowb08PrmoTxlVtwli0aS9jZmWROncDH83OIjkmhNsvjqZr07r8sX2pfHG6pe9boFPBOnsAGGOigcSCY2eSXMxrEa6JIcnAuIJ6PwPHgKJTZW4GlhbMFhYREREpUZsGtRj5p7bMfCKFB1Kasjh7L7d+MI9eI6fx6ZwN5Bw77nSIp+V0S99o4H7gG2PM07jG3r0IbMLVfQuAMSYKWAu8YK19AcBam1b0YsaYvYBX4WPW2u3GmJHAYGPMAWA+8CcgBTeXhREREREpLKSmH4/0asZ9yY35ftHvfDhrPU9/vZRXf17BDR0jubVzFBFBASVf6AJyNOmz1h4yxqTgWnblE1yTKyYBD1lrDxaqagBPzr5lcghwEHgQqAesBK631n53trGLiIiI+Hp5cm37CK5pF076hj2MmZnF+zPW85/p67gkrh63J0bz+97DvBY/kC0+NQkbPplBvWPoFx9+wWN1uqUPa+1G4Iz74Fprs3BjRq+1Nuk05XnASwUvERERkTJljKFDdG06RNdm897DfDJ7A5//tpGfM7e6lhHxDQRg897DDB6/BOCCJ35Oj+kTERERqVTCa/nz5GWxzH6yB7X8vf+wbtzh3DxG/LLygselpE9ERETkPPD38WTf4dxij23Ze/gCR6OkT0REROS8CavlX6ry80lJn4iIiMh5Mqh3DP7enqeU+Xt7Mqh3zAWPRUmfiIiIyHnSLz6cYde0IvzoPoy1hNfyZ9g1rarm7F0RERGRyqxffDj9Hn7P9SEtzbE41NInIiIiUgUo6RMRERGpApT0/X979x5sV1mfcfz7hHhDrI3TILZCgOKFOF7GRssUsd4RL+BURbwyDop1qoXR1kpVegaZShEv9TIFWjt2pA4qxYqIQkVINdPYpnZiiVXAhlsrA2kwyFWQX/9Y6+hhZ8ezc/bZZ51z1vcz884O716Xd+Vl7/XkXWu9W5IkqQcMfZIkST1g6JMkSeoBQ58kSVIPGPokSZJ6wNAnSZLUA4Y+SZKkHjD0SZIk9YChT5IkqQcMfZKk0U1NQfKLsn59U2bWJc1ykhaVlV03QJK0hExNGeikJcqRPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk9YOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHOg99SfZNcl6SHUluTXJ+kv1GWG9Nki8luTbJnUm2Jbk8yRFDlr0mSQ0pL5vMUUmSJC0uK7vceZI9gW8AdwPHAgWcClyW5ElVdfsvWX0vYBvwXuAG4FeANwMXJXl5VZ0/sPzFwNRA3Q/GYuvEAgAADblJREFUPghJkqQloNPQRxPSDgQeV1VXAyT5LnAV8Bbgw7tasaq2AMfNrEvyFWAr8EZgMPRtq6qN89d0SZKkpaPry7tHAhunAx9AVW0FNgBH7e7GqupeYAdwz7y1UJIkaRnoOvQ9AbhiSP0WYO0oG0iyIsnKJPskeR/wWOCTQxZ9aZI7ktydZKP380mSpD7pOvQ9ArhlSP12YNWI2zidZmTvR8C7gGOq6tKBZb4MvB04HHgtcBfwxSSvm0ujJUmSlpqu7+mD5uGNQdmN9T8KnAvsA7wB+GySV1TVhT/fQdXb77fx5IvARuADwDnDNprkeOB4gP32m/VhYkmSpEWt65G+W2hG+watYvgI4E6q6oaq2lRVF1bV0TRh7oxZ1vkZ8AXg0UketYtlzq6qdVW1bvXq1aM0RZIkadHqOvRtobmvb9Ba4Htz3OYm4KARlpseTRw20ihJkrSsdB36LgAOSXLgdEWS/YFD2/d2S5IVwDOAH86y3ErglcB1VXXj7u5HkqRFZWoKkl+U9eubMrMuaZZTb3V9T99fA28DvpTkvTSjbu8HrgfOml4oyRqaIHdKVZ3S1k3RXBreANxIc0/fccDTgdfMWPfVNNO/XNRu95HAHwC/Bbx6okcnSdJCmJoy0GlWnYa+qro9yXOAjwCfobnkeilwYlXdNmPRAHtw/5HJ7wAnAscAD6cJfpuBw6pqw4zltgJ7Ax+kCYl3AP8GvLCqLp7EcUmSJC02XY/0UVXXAS+fZZlrGHiit6ouYIRLwO2vcDxnjCZKkiQteV3f0ydJkqQFYOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk9YOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk9YOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk9YOiTJEnqAUOfJElSDxj6JEmS5tvUFCS/KOvXN2VmXdIst0BSVQu2s6Vq3bp1tWnTpq6bIUmSNKsk/15V6wbrHemTJEnqgc5DX5J9k5yXZEeSW5Ocn2S/EdZbk+RLSa5NcmeSbUkuT3LEkGVXJDkpyTVJ7kqyOcnLJ3NEkiRJi0+noS/JnsA3gMcDxwKvBx4DXJbkobOsvhewDXgv8CLgOOA24KIkvzew7PuBKeATwBHARuALSV40P0ciSZK0uK3seP9vBg4EHldVVwMk+S5wFfAW4MO7WrGqttAEvZ9L8hVgK/BG4Py2bm/gj4DTquqMdtHLkhwEnAZcNJ8HJEmStBh1fXn3SGDjdOADqKqtwAbgqN3dWFXdC+wA7plRfTjwQOCcgcXPAZ6Y5IDd3Y8kSdJS03XoewJwxZD6LcDaUTbQ3q+3Msk+Sd4HPBb45MA+7gauHlh1S/s60n4kSZKWsq4v7z4CuGVI/XZg1YjbOB14Z/vn24BjqurSgX38uHaem2b7jPclSZKWta5H+gCGTRSY3Vj/o8DTgJcCXwU+m+QlA9va7X0kOT7JpiSbbr755t1ojiRJ0uLTdei7heEjbasYPgK4k6q6oao2VdWFVXU0zZO5Z8xYZDuwKslgyFs14/1h2z27qtZV1brVq1eP0hRJkqRFq+vQt4XmnrtBa4HvzXGbm4CDBvbxIOA3h+yDMfYjSZK0ZHQd+i4ADkly4HRFkv2BQ9v3dkuSFcAzgB/OqP4a8FPgtQOLvw64on1aWJIkaVnr9Ld32wmYNwN30kyyXDQTKT8MeFJV3dYut4YmyJ1SVae0dVM0l4Y3ADcC+9DM2/c84DVVde6M/ZwGnAj8KfAd4FU08wAeVVVfHqGdNwPXjn/Ev9Sv0Uw2rcXFfll87JPFxz5ZnOyXxWeh+mRNVe10b1qnT+9W1e1JngN8BPgMzcMVlwInTge+VoA9uP/I5HdogtwxwMNpgt9m4LCq2jCwq/fQPNl7Ak04/AFw9CiBr23nxG/qS7Jp2I8jq1v2y+Jjnyw+9sniZL8sPl33SddTtlBV1wG/9Hdwq+oaBp62raoLGPEScFX9DDi1LZIkSb3T9T19kiRJWgCGvsXj7K4boKHsl8XHPll87JPFyX5ZfDrtk04f5JAkSdLCcKRPkiSpBwx9E5Rk3yTnJdmR5NYk5yfZb8R1H5zkg0l+lOTOJP+S5JmTbnMfzLVfkqxLcnaS7ye5I8l1Sf4+yQEL0e7lbJzPysB2TkpSSb41iXb2zbj9kuTgJF9Isq39HvtBkhMm2eblbszzyn5J/q797rojyZVJTm2nT9McJXl0ko+35+k72u+g/Udcd0HP9Ya+CUmyJ/AN4PHAscDrgccAl434AfsU8GbgZOAlwI+Ai5M8ZTIt7ocx++UYml+Q+RhwBPBu4KnApiT7TqzRy9w8fFamt3MgzfRMN02inX0zbr8kWQd8m+YXkd4EvAj4EM30W5qDcfqkff/rwDOB9wEvBv4GeCfwtxNsdh8cBBxN8/Ox39zNdRf2XF9VlgkUmjkBfwYcNKPuAOBe4B2zrPtkmomq3zijbiXN/IIXdH1sS7mM2S+rh9StAe6jmTi88+NbimWcPhnYzsXAWcDlwLe6Pq6lXsb8rKyg+QnML3Z9HMupjNknL2jPKy8YqD+tXX/Pro9vqRZgxYw/v6n9e95/hPUW/FzvSN/kHAlsrKqrpyuq+cm3DcBRI6x7D/C5GeveC5wLHJ7kQfPf3N6Yc79U1c1D6q4FbgZ+Y57b2SfjfFYASPIamlHXkybSwn4ap1+eRfP75h+eWOv6aZw+eWD7eutA/Y9pQnrQnFTVfXNcdcHP9Ya+yXkCcMWQ+i00X4azrbu1qu4Ysu4DaYaSNTfj9MtOkhwM7A3815jt6rOx+iTJKppf9XlXVW2f57b12Tj98oz29cFJNia5J8lNST6W5CHz2sp+GadPvg5cBfxFkrVJ9mp/EesE4Myqun1+m6oRLPi53tA3OY+gub4/aDuwaox1p9/X3IzTL/eTZCVwJs1I36fGb1pvjdsnHwSuBD49j23SeP3y6+3r54BLgOcDp9Nc+vrsfDWwh+bcJ1V1F00Yn770/hOanz29EHjb/DZTI1rwc33nP8O2zA2bBHGUIfSMsa5mN19/t58Afgd4cVUN++BqdHPqkySHAW8AnlrtDTGaV3P9rEwPKJxTVSe3f748yR7AaUnWVtX35qWF/TPXz8qDaUL43jQPgFwHPJ3mAYJ7gbfOYxs1mgU/1xv6JucWhqf0VQxP9jNtB4Y9gr9qxvuam3H65eeSfAA4Hji2qi6Zp7b11Th9chbNKOsNSX61rVsJ7NH+951Vdfe8tbRfxumX/2tf/2mg/hKaBweeAhj6dt84fXIczb2WB1XVD9u6f06yAzg7yZlVtXneWqpRLPi53su7k7OF5nr9oLXM/mW3BTigfTx/cN2fAlfvvIpGNE6/AJDkPTTTtZxQVZ+Zx7b11Th9cjDw+zQnvOlyKHBI+2dHL+Zu3O8w2HkUY3oEY643vvfdOH3yROCWGYFv2r+2rweP2TbtvgU/1xv6JucC4JB27jAA2skaD23fm23dBwCvnLHuSuBVwCWOXIxlnH4hyR8CpwLvqaqPT6iNfTNOnzx7SNlMc7P7s4Hz5r+5vTFOv3wVuBt44UD94e3rpvlpYu+M0yc3AquSDD4c8Nvt6//MUxs1uoU/13c9v81yLcBDaVL6f9I8Sn8kzcnov4G9Ziy3huZ+ipMH1j+XZqTiTcBzaU5ed9Hcu9T58S3VMk6/0EzOfB/NCe2QgbK262NbqmXcz8qQ7V2O8/R13i/An7X1fw48j2Z0/E7g010f21ItY35/7U8zXcuVNBM7Pxv447ZuEzPmmrPMqW9e0Za/ohnhfmv737+7qz5p6xf0XN/5X9RyLjTX6v+h/VD9BPhHBiZsbD+IBUwN1D+EZo6rG9v/Ab4NPKvrY1oOZa79QvN0aO2iXN71cS3lMs5nZci2DH2LoF9oLuW+ow0pPwWuBU4BHtD1cS3lMmafrAU+D1xPE8CvBM4AVnV9XEu9zHZuWCzn+rQ7lSRJ0jLmPX2SJEk9YOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk9YOiTJEnqAUOfJC2QJHslOTnJfyT5SZLaRdm767ZKWn5Wdt0ASeqDNsitBx4PfBc4E3gQ8EpgH+Ae4DpgW1Xd1FU7JS1fqaqu2yBJy16SrwPPBU4H3l3tl2+SfYGrgD2AR1XVtu5aKWk58/KuJE1YkufTBL4NwEk141/bVXU98E2aKy9P6aaFkvrA0CdJk/e69vUjVXXfkPd3tK9+J0uaGL9gJGnyDgPuA762i/cf3b5evTDNkdRH3tMnSROUZAXNQxo3VdWjhrz/SOB/gWur6sCFbp+k/nCkT5Ima/pf1g9rA+Cgd9F8F5+1cE2S1EeGPkmaoPahjc3AQ4FXz3wvySuAE4HvA3+58K2T1Cde3pWkCUvyMuB84F7gc8D1wNOA59FM13J4VW3troWS+sDQJ0kLIMlRwJ8ATwYCXAl8HvhYVd3WZdsk9YOhT5IkqQe8p0+SJKkHDH2SJEk9YOiTJEnqAUOfJElSDxj6JEmSesDQJ0mS1AOGPkmSpB4w9EmSJPWAoU+SJKkHDH2SJEk98P82xJCm3lpTIwAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "# Instead want:\n", + "# 1. FC-FC as Blue circles\n", + "# 2. Mean +/- std for null distribution for each shuffle in red. No need for sigma = 0.0 \n", + "\n", + "mean = df.groupby('sigma').mean()\n", + "std = df.groupby('sigma').std()\n", + "\n", + "# Plot FC-FC blue\n", + "plt.plot([0,0.2,0.4,0.6,0.8,1.0],array,'-o',zorder=2)\n", + "# Null dist in red (shuffle dist)\n", + "plt.errorbar([0,0.2,0.4,0.6,0.8,1.0],mean[\"FCFC\"],yerr=std[\"FCFC\"],fmt='ro',capsize=5,zorder=1)\n", + "\n", + "plt.legend(('Model','Null Distribution'))\n", + "\n", + "plt.xlabel(\"$\\sigma$\",fontsize=20)\n", + "plt.ylabel(\"FCFC\",fontsize=20)\n", + "plt.xticks(fontsize=16)\n", + "plt.yticks(fontsize=16)\n", + "plt.savefig(\"do-not-track\\\\Het_sweep_reviewer3.pdf\",bbox_inches='tight')" + ] + }, { "cell_type": "code", "execution_count": 7,