Skip to content

Instantly share code, notes, and snippets.

@andrzejnovak
Created February 9, 2024 09:48
Show Gist options
  • Save andrzejnovak/b5531edb78c153df7d1a9280fde68320 to your computer and use it in GitHub Desktop.
Save andrzejnovak/b5531edb78c153df7d1a9280fde68320 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 46,
"id": "b189bc69-10e0-4aea-b40c-fd3dab00c82e",
"metadata": {},
"outputs": [],
"source": [
"import ROOT as r\n",
"import plt\n",
"import np"
]
},
{
"cell_type": "code",
"execution_count": 52,
"id": "aacb6655-f43d-4317-b2ff-92946343e63f",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PolyCollection at 0x7fdc61f38d50>"
]
},
"execution_count": 52,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGdCAYAAABO2DpVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABOaklEQVR4nO3de5yWc/7H8dfMpCnUEDqZknOSCtGGxIoklk2OIac9hpJj64e1i8Qip81pV3bXoUU57WITHRwiJeSQENJRqpmKppr7/v3xua7ue2oO9z1zXfd1uN/Px2Meu9fdPff9bUrX+/5+P9/PtyCZTCYRERER8UBh0AMQERGR+FCwEBEREc8oWIiIiIhnFCxERETEMwoWIiIi4hkFCxEREfGMgoWIiIh4RsFCREREPNMo12+YSCRYtGgRzZo1o6CgINdvLyIiIvWQTCZZvXo1bdu2pbCw5nmJnAeLRYsW0a5du1y/rYiIiHhgwYIFlJaW1vjrOQ8WzZo1A2xgzZs3z/Xbi4iISD2Ul5fTrl27TffxmuQ8WLjLH82bN1ewEBERiZi6yhhUvCkiIiKeUbAQERERzyhYiIiIiGcULERERMQzChYiIiLiGQULERER8YyChYiIiHhGwUJEREQ8o2AhIiIinlGwEBEREc8oWIiIiIhnFCxERETEMwoWIiIi4hkFCxERkVioAJ4EzgNWBTYKBQsREZFI+wq4GmgDnAGMBT4IbDSNAntnERERqaeNwH+A+4BXsXmCykBH5FKwEBERiYyFwMPAGGApUAQkCUuoAAULERGRkEtgsxJ/BV5IewzCFChcChYiIiKh9D1WL3Ef8A02O5Go7RtCQcFCREQkNJLAm9jsxFNYkAjv7ER1FCxEREQCVwb8E7gXmIvdnjcGOqL6UrAQEREJzEysEPMxrA+FK5qhAhQsREREcmwtMA6bnXifKM9OVEfBQkREJCc+wWYnxgJrSPWojE+oAAULERERH1UA47GdHW9SdXYi/Ds86kPBQkRExHNfAQ86XyuJ6+xEdRQsREREPOG22b4Xa2hVRGqLaDxnJ6qjYCEiItIg1bXZhqj0nfCagoWIiEjWotVmO5cULERERDIWzTbbuaRgISIiUqvot9nOJQULERGRasWnzXYuKViIiIhUEb8227mkYCEiIkIl8A/gbmA2mp2oPwULERHJc2uBU4H/kk+NrPyiYCEiInlsCXAsMMe51g6PhlKwEBGRPPUJcAwWLrS7wyuFdT9FREQkbl4HeqBQ4T0FCxERyTP/wmYqfkShwnsKFiIikieSwJ+Bs7HiTNVT+EE1FiIikgc2AL8BHgl6ILGnYCEiIjFXBgzA6irEbwoWIiISYwuAvsDn2FKI+E01FiIiElOzge7APFSkmTsKFiIiEkMvA4cAP6AumrmlYCEiIjHzENAfWIdmKnIvq2Dxxz/+kYKCgipfHTt29GtsIiIiWUgAI4BfO/9fNRVByLp4c9999+XVV19NvUAj1X+KiEjQKoBzgScDHkewyoBngHPYGNjujKzft1GjRrRu3dqPsYiIiNTDCuAXwNtBDyRQLwPnYU3K9+QTenFUIOPIusZi3rx5tG3blt12241Bgwbx7bff1vr8iooKysvLq3yJiIh44yvgYGA6+dpJcxVwPtAPWOo8lgjwZ5FVsOjRowdjx47l5ZdfZsyYMcyfP59evXqxevXqGr9n5MiRlJSUbPpq165dgwctIiIC7wIHAd+Qr0WaLwEdgX8412GoKilIJpP1HseqVavYZZdduOOOO7jggguqfU5FRQUVFRWbrsvLy2nXrh1lZWU0b968vm8tIiJ57VngdKxVd/7NVKwChmMNygvZ8icwmdH0Zqin71leXk5JSUmd9+8G1XZst9127LXXXnzxxRc1Pqe4uJji4uKGvI2IiEiau4BLnf8fhs/oufUSVkux3LkOW6xqUB+LNWvW8OWXX9KmTRuvxiMiIlKDSmAoMAwLFPkVKlZhtRTHAd8T3sWfrILF5ZdfzpQpU/j666956623+OUvf0lRURFnnHGGX+MTEREBfgROBu4JeiCB2LyWImyzFOmyWgr57rvvOOOMM/jhhx/YaaedOOyww5g+fTo77bSTX+MTEZG8twz7nP4++ThLcSkwluprKcIoq2Dx5JP53XhERERybS5wDLCQaNxWvRP2Woqa6KwQEREJqWlADyxUhLWiwHursEAR9lqKmqgft4iIhNA44Czsc3pUPqs3XFRnKdJpxkJEREIkCYzCelRsJJq31uytItqzFOk0YyEiIiGxERgCPBj0QHLqv9g20ijPUqRTsBARkRBYDQwEJgY9kJxZhXXkeBQoID77XRQsREQkYIuAvsCnxOf2WrvNZyni9LtWjYWIiAToI+BA4DOiXVmQmVXAuUB/rDtHHH/HmrEQEZGAvAqcBKwjnrfYqv6LFWj+4FzHaZYinWYsREQkAGOBY7FW3fEOFatIzVJEfcdHJhQsREQkh5LAddhn90ri+7nd/BfYG/iXcx3v363RUoiIiOTIeuACUrfZ+FpFPHd8ZELBQkREcmAVVk8xLdhh5EC+1FLUREshIiLis2+AnwFvEP32TzVbCQzGaimWE/9aippoxkJERHw0C+tRsZI432r/g/WlcGcp4huf6qYZCxER8cmLwKHEOVS4sxTHk9+zFOkULERExAdjgF8AFcT1dvsfoCPwmHOdz7MU6RQsRETEQwngcuD3WNli/EoXNUtRO9VYiIiIR9YBZwHPBD0Q36iWom4KFiIi4oHl2Gf4GUEPxBcrgaHAP7GpfgWKmilYiIhIA30BHAN8SxxvuZqlyI5qLEREpJ4+A24BDgIWELdqg5XAOaiWIluasRARkQwlgOnAc8BTwHzs82n8ijRfxGYpVjjXmqXInIKFiIjUYh0wCZgAPIstCDQCNjq/Hq9b7krgEuw0E9VS1I+ChYiIbGYFVlkwAXgJCxfpYWJjDd8XbZql8IaChYiIYOd5PIdtFXXP9CgiVVkQzzCRwALFKOAtNEvhBQULEZG8lAQ+wJY3ngHmsGW9RHzLFddhW0dHAV9iEQoUKrygYCEikjc2YMeWPwc8DSyi6qxE/G+rP2DNxu/EljwKnMfjG6FyT8FCRCTW1gCvYPUSzwOrqVovkR+31K+AO4CHsXjlRqh47WUJBwULEZHYWQK8gIWJV7FbafyLL6vzLnArMB5b6MmPGBUsBQsRkViYS2qJ4z3ss3h6KWL+hIkEtqdlFPAmdqNLolCRKwoWIiKRlMA+jz+LhYkv2bL4Mv41E+nWYf0nRmFNxt2CzPyJVOGgYCEiEhnrgNewMDEBazQd32ZVmVpBqiDzB1SQGTQFCxGRUFsJ/BcLEv8FfiJf6yU2N59UQeZ6VJAZFgoWIiKh8y1WLzEe2x5aST40q8rUDKwg8xlUkBlGChYiIoFLAh+SKr78CJvQLyD1OTy/b58JbL7mFlSQGXYKFiIigUgAU0l1vvyOqrMS8TsxtD7WAY9hBZnzUEFmFChYiIjkXBL4FfB38rFZVSZWAPdjNRQqyIwWBQsRkZy7DQsVoM/eVc3Hdnc8hAoyo0rBQkQkpyYAVwc9iNB5DyvIfBoVZEadgoWISM7MAs4MehChkQBewuonpqGCzLhQsBARyYmFQD/s3I78ntivIFWQ+TkqyIwbBQsREd+txULFCvL58/hKUgWZy1FBZlwpWIiI+CqBLX98Qr7eQr/GCjIfRAWZ+UDBQkTEV1djR5jn3210Jrb/5SlsdiI/Y1UOVWLFKovh/TbzOKxXJUVFRXV9l+cKc/6OIiJ542/YrTV/QoXbIbM30B1r/ZVAocJ344EOwJHAmXDpkffRoUMHxo8fn/OhKFiIiPhiMvCboAeRMxXAI8A+QH+s7TaoIDMnxgMDseataRYuXMjAgQNzHi4ULEREPPc5cCL5MFPxE3Z+RzvgfKztNmiGImcqgaFU+1ctmbQHhw0bRmVl7v5EFCxERDy1AjgW2wmSqOO50fYjttflGuB757H4R6mQmcYWMxXpkskkCxYsYNq0aTkbkoKFiIhn1gMnYceex/szu7uBdhpxj08htzjDpy3O8IkeaFCwuOWWWygoKGDYsGEeDUdEJKqSwG+BN8iHUHEc9jtVqAjYl5k9rU2bNv6OI029t5vOmDGDBx54gC5dung5HhGRiPoLVr4Yb+5MxZsoVARqNVZbUcdfuYKCAkpLS+nVq1cuRgXUc8ZizZo1DBo0iIceeojtt9/e6zGJiETMs8BVQQ/Cd26oeAuFikC9AXTFQkUBtvpWQKqVqaOgwB4YPXp0TvtZ1CtYDBkyhP79+9OnT586n1tRUUF5eXmVLxGR+JgFnBH0IHy3FitJfYu4L/SE2HpgBHA4dr78Ltiu5gnYsbA7V316aWkpTz/9NAMGDMjpMLNeCnnyySeZNWsWM2bMyOj5I0eO5IYbbsh6YCIi4ZcfB4utwULF22imIjAfA2cBs53rwcDdQHPnegC2w9npvHlnmyFc3Ouu8HfeXLBgAUOHDuWxxx6jSZMmGX3PiBEjKCsr2/S1YMGCeg1URCRc3BLGH4jzZ/g1QF9gOgoVgUgAdwEHYqFiB2x2YiypUOEqAo4AzoD9j9gzkFABWc5YzJw5k2XLlnHAAQdseqyyspKpU6dy7733UlFRscVvpLi4mOLiYm9GKyISCu7BYnOI8+3WDRXvEOfoFGLfAecBrzrXxwJ/B3K3waNesgoWRx11FB999FGVx8477zw6duzIVVddFVg6EhHJrRHA80EPwlersVDxLgoVgRiH7V5eBTQFbneuC2r5npDIKlg0a9aMzp07V3lsm222YYcddtjicRGReHoEuDXoQfhqNXAMMAOFipxbCVwEPO5cHwT8E9g7sBFlTZ03RUQyNhn4VdCD8JVCRYBeA7pgoaIIuA5rGBKhUAENaJDlmjx5sgfDEBEJu3lYw4D47v4ox0LFeyhU5NQ64A/Anc71Htgsxc8CG1GDNDhYiIjE3wqs4mANcS3WLAeOBmaiUJFTs7FtpB8717/BmrhuG9SAGk5LISIitVoP/JI4HyxWDvRBoSKnKoFRwMFYqGgJvAjcT6RDBWjGQkSkFkngd1jXoXgugZRhMxWzUKjIma+Bc7C/VmCNrR4CdgpqQN7SjIWISI1uxxoHxDdU9EGhImeSwKNYgeY0bGbib1hL7piECtCMhYhIDZ4Frgx6EL5xQ8X7KFTkxHKsD8UzzvWhwD+A3QIbkW80YyEisoX3ifPBYquAn6NQkTMvAfthoaIRcDMwhViGCtCMhYjIZuJ9sNgq4CjgAxQqfPcjcAXwV+d6H+BfwAE1fkcsaMZCRGSTeB8stgqFipyZAexPKlQMxbbdxDxUgIKFiIgjgTUUmANsDHgs3luFLX8oVPhsI/AnoCfwOdAW+B8wGjvzIw9oKUREBLCDxZ4NehC+WImFio9QqPDVPOBs7DhYgNOwGYsWgY0oEJqxEBGJ8cFiK4EjUajwVRJ4AOiGhYoS4DHgCfIuVIBmLEQk700Bfh30IHzhhoo5KFT4ZilwIdY1E+wHPhZoH9SAgqcZCxHJY/OwtofxO/9jBXAEChW+eg7ojIWKYuAO4FXyOlSAZixEJG+tAI4ljgeLuaHiExQqfLEauBTrmgnWSfMxLGSIZixEJB9twA4W+4a43XoVKnz2JtAVCxUFWHPWd1GoSKMZCxHJM/E9WOwHLFR8ikKF59YDNwC3YBNcu2AtuQ8PclDhpGAhInnmDlJz2PGhUOGjT7EWJ7Oc68HAXdjuD9mClkJEJI88h/VYjpflQG8UKjyXAO7GumXOwraOPo3t+lCoqJFmLEQkT7wPnB70IDznhoq5KFTUWyW2MrYYaAP0ApYA5wETneccC/zd+XWplYKFiOSBRcTxYDGFCg+Mx87x+C7tsRZABXZ0TFPgL1hZTkHORxdJChYiEnPxPFjseyxUfE6cflc5Nh4YyJZZc4Xzv7sD/wH2zuWgok/BQkRizD1Y7CPi1KtCocIDldhMRW0TWBXAHrkZTpyoeFNEYuwP2MFi8QkVy7AdjgoV9bARq3B9Ajss7Lvan853WO2FZEUzFiISU2OBUUEPwlNuqPgChYo6rQI+xM6Jd7/mAOuyfJ3F3g4rHyhYiEgMTQF+FfQgPKVQUYMEMJ9UeJjt/O83NTx/G2A/YCfghQxeX7tAsqZgISIx8wVxO1hMocKxFiuXSZ+F+BA77qU67bH22+lfu2NFAJVAB2Ah1ddZFACl2NZTyYqChYjEyEqgL3E6WGwpFiq+JI9CRRKrb/hgs695VB8CioF9seDQzfnfLsD2tbxHEdY9cyAWItJf191WOtp5nmRFwUJEYiJ+B4uFPlRU11gq2xtxBXZiWvoyxoektnxurjVbzkLsTf3uZgOwTpqb97EoxULFgHq8pihYiEgcJIHfA1OJSwOsJVio+IqQhorqGkuVYrMANd2Ql7LlLMRn2G6NzTUCOrJliGjlwdjTDcBWzhoakGQTBQsRiYE7gIeDHoRnlmD3tvmEOFRU11hqofP4OKATVQPEbCxYVKcFWwaITtgSRy4UYSe4iScULEQk4p4nTgeLLcZmKr4mpKGitsZS7mOn1vC9BcCeVA0Q3YCdUbvsGFGwEJEImw2cFvQgPLMYm6n4hupXB0JhGnU3lgI7Y+MAqoaIzth2T4k1BQsRiahF2JGT8ThYLBKhAjJvGPUQMMjPgUhYqaW3iETQj0B/7HzPUC4YZGURcBgRCBUALTN83s6+jkJCTDMWIhIx7sFiHxKHXhWLsJmKb4lAqFiNHSFeGzWWynuasRCRCFkLnAdMIA6hYiE2UxGJULEAG+zLQGPnsc0LLtVYSlCwEJHIeB+rAPxX0APxxKfYh/oFRCBUzAR6YJNErYA3gGfYcrmjFGs4pcZSeU1LISIScgms69KVWJFmdGcqksCb2JmrL2If6kNfIfIccCZW1tIZG/guwEGosZRUS8FCREJsKXA2MDHogTRIJfAsMBL78N8o7fHQSmJLGpc5//8Y4N9ASdpz1FhKqqFgISIh9TK2X7Es6IHU24/AI8CtWB2Fu/Yc+qWPjcAlwBjn+jfAveiOIRnRXxMRCZkK4Grs4/Lmx05GwzLsPnw3UE7qdxCJRZxyrOfYy9iP/zZgOOqMKRlTsBCREPkMOAX42LmOVqj4HNuNORZb5ohEkEj3LXA88BHWOfNx4KQgByRRpGAhIiGQBP4GXEwUO2mmF2QWEvLaiZq8B5yAnYDWGngB6B7oiCSiFCxEJGArgQuxIzOjoxLbMHELMAOrY0wS0VDxLLbz4ydgPywhtQ9yQBJlChYiEqBp2IL+sqAHkrEfgUexgsyvSRVkRjJQJLET569w/v+x2JHnzYMclESdgoWIBGAj8CfgRqwqMPzVCN8D92EdNdL3qYR/5DXYCFwEPOBc/w6rNtVdQRpIf4VEJMe+AU4H3sE+Joe7nmIecDu2bXQjEQ4S6cqBU4FXsFx3OzAM7fwQTyhYiEgOjcPqKdYR9kDxFrbc8TwRLsiszjfYzo85wNbYzo8TAx2RxIyChYjkwBpsx8dYwtybohILErcA72L/QEa2ILM6M7CdH0uxFtwvAAcGOiLxWBFQDOxO28DGkNUhZGPGjKFLly40b96c5s2b07NnT1566SW/xiYisTALOzzsH851+ELFT8D9wJ7Y+VnvOY+HvkNmNsYDvbFQ0QVbiVKoiA33Zn4SMBcopWXgY8lIaWkpt9xyCzNnzuS9997j5z//OSeeeCIff/xx3d8sInkmgS3e98Dm38NXnbAcuAFoC/wemO88Hr6RNkAS6555MpagjsNOJ20X5KDESwXYuXCvYIfLlgY7HAqSyWSDPj60aNGC2267jQsuuCCj55eXl1NSUkJZWRnNm2tPk0g8LQHOAiYFPZBqfUGqIHMDMQsS6TYAQ4CHnOuLgDvRInhMFGF/lNdhZ8UVV/nVydgUlXcyvX/X+69XZWUlTz31FGvXrqVnz541Pq+iooKKiooqAxOROPsvFirC99/6dKxD5nPErCCzOmVYd/SJ2Efa0djBYhJ5hVgYPg7bIdwh0NFsKetg8dFHH9GzZ0/WrVvHtttuy4QJE+jUqVONzx85ciQ33HBDgwYpIlGwDrgK+6fO/acveAmsRvEWLFjEriCzOl8D/YFPgG2AJ7CiTYm8AmBn7ODZ/gGPpSZZL4WsX7+eb7/9lrKyMp5++mkefvhhpkyZUmO4qG7Gol27dloKEYmVT7GPx58SlkDxE1YueivwFTZtHOsw4XoH+AXWzLQt1p57/0BHJB4owuL6H7D43rTO75hMUEshDa6x6NOnD7vvvjsPPPBA3U/OYmAiEgVJbAH/EmwPRfC37uXAX7GZ/5WEeXOrD54GzsYmj7phUzVBV/JJg7h/f4/BOr/ukfF3TiZyNRauRCJRZUZCRPLFCuAC7ASr4H2JHXvxN6oWZOZFqEhiUzNXO9f9gSeBbQMbkXigEGiFBYqTiE5j1KyCxYgRI+jXrx/t27dn9erVPP7440yePJlXXnnFr/GJSChNwdpyfx/0QHgHK8h8ljwoyKzOBuycj78515dgCasosBFJAxVhIeIy4FqsTCZKsgoWy5Yt45xzzmHx4sWUlJTQpUsXXnnlFY4++mi/xiciobIR+CNwM0EeHpbASgdGYa2386IgszqrgIHYrt5CbP3n4gDHIw3iLnschjVs6xjscOotq2Dxt7/9re4niUhMzceOOH+PIA8PexE7L+tLUh/KY9UhM1PzsSWPT7GPtOMI7zYBqVMhsCO2p+pUorPsUZ2sOm+KSL56AtgPa88dTKCoBP4P2zX5VdpjeWk61tD0U2zv4RsoVESUu9tjKNa47TSiHSpA/ddEpFarsXaN/yDI/RUrsH9w3T6eeVGQWZOngHOwnR/7Yzs/dg50RFIP7n9NPbBlj/2CHY6nNGMhIjV4Dzut6jHnOpjb+fvYEWavBzaCkEgCI7F58nXY1M1UFCoiqAjYHovrbxCvUAEKFiKyhQS2d/FnwAKCXHD4hzOKxYGOIgTWAxdi3ZHAikwmoO2kEePu9vgtViN0NtFf9qiOlkJEJM1iYBA2PxCc9di9c0ygowiJldjOj9ewj4L3YEexSuTsDzwAHBD0QHymYCEijhexz1BrAh3FIuCXwIxARxESX2FFmZ9hsxP/BvoFOiLJUhH2R3c7cB75sUyQD79HEanVOqz5wQnYkZjBbd6cilV1BLf3JETextaBPsPacr+BQkWEuMse52PLHheQPzdczViI5LWPscPD5jrXwdzOk9j+/eHOdTiOMQvQOGAwUIHNm7+AHSgmkbEvdorOwUEPJAD5EqBEpIoktsntAOBzgryVrwXOxGoqEoGOJASSwE1Yt/QK4ERsGkehIhKKgGbYIXizyM9QAZqxEMlDP2ATtM8HPRC+wE74nlvXE/PBeuDXwKPO9XBsc47O/Ai9ImzX0lnYH1nLYIcTOAULkbzyOnAGdrh4sP6DfTD/iTycpagEpmGbcNpg8+anYiddF2E7P34X1OAkW3tjuz0OC3ogIaFgIZIXktjBYdcS5OFhOO98A/AnguzlGaDxWP/m79Iea4TVzDbDdn4cG8C4JCtFQDH2X9UQdDNNp5+FSOwlsc5Kt6RdB2MlVk/xcuAjCch4rCfF5r9xdyPODShUhJy77HEKdjp9m2CHE0oq3hSJtSRwGalQEZwPgW7AxIDHEZhKbKaipjRVANxJnrcYDbcCYDesV9kTKFTURMFCJLYSWH+KO4MeCI9hFfILyeP75lSqLn9sLol1UJ+Wm+FI5oqAJlg8nwMcGexwQk9LISKxlAB+Azwc6Cg2YPMl9wQ6ioB9hvWleDDD5y/2cSySlULsv6QTgdFAu0BHEx0KFiKxU4n1+Xu0rif6ajFwMjA90FEE5CssTIwDPsjyezW/HgoFwC7YeTV9Ax5L1ChYiMTKRuy8j3GBjuIt4CRgBXlUoLkA29ExjqoHnTQCjsGq/a7BEld1P5QCrHV3L3+HKbVrhM1UXAtcge38kOwoWIjExgasR8V4gmzNfR/WRTNJHvSnWAw8DTyJpSlXIfBz4DTsRLUdnMebY7tCNt9n656dPRo1xArAVth/PWCbcu4Gdg1uOJGnYCESCxXYR+IXCSpU/Ig1jnwskHfPoeXAM1iYmELqx12AzTachq0BtarmewdgQWTzPhalWKgY4MuIxdEIWyh0/8h2wA696wp0xo41j/uR5rmgYCESeeuwj8X/I6hQ8RXWmvvTQN49B1YCE7BljklU3dryM6yF6EBg5wxeawBWDZjeebMXmqnwkPujdP+YtsGam+6PBYjOzvVOuR9aXlCwEIm0H7HjzicT1MLDS9h99cfARuCTcuw4lXHAK6TmygEOxGYmTsUq/LJVBBzRwPEJhc6X219sK6y9thsg9nP+t5TUapP4T8FCJLLWAMcBbxLELT2BtTO+zrmORZHmWuwQk3HO/1ak/dp+WJg4Ddgj90PLZwXYzcrNdoVYDUQ3UuFhP6x5lW5qwdOfgUgklWOb4GYQRKhYhZ3k+J+cv7MP1mE9xp8EXsCmXlx7kwoTnXI/tHyUXkgJtlLUDauFcJcxOmINqyScFCxEImcV0AeYTRB9LOdgiy8Lcv7OHloPvIqFieewnObaFQsSp2N3M82h+6IRFondWLwdqUJKdxaiE1ASxOCkQRQsRCLlB2wf48cEESrGAedinygj15p7I3Zq/DhsR+7KtF8rJTUz0R2FCQ9tXkjZFCuc7EbVOoiWOR+Z+EXBQiQyvscq/uaS69v6BuBKbEdkpCSw3RfjsG2e36f9Wmtsh+5pQE90clIDFWAhwi2kbATsRdWdGJ2B9uhHHXcKFiKRsAToDXxJrkPFUmwn5Zs5fdcaVFL3Ns0k1kd8HPAUsCjt13bEekycBhxezfdK1gqwH+uhVK2D2AOrl5D8o2AhEnoLsbvgN+Q6VEzHWnMvJwS7PsZTfWOpu7A2HrNInc/xbdpztnN+/XTsWErd7TxRhM08XANchYopJUXBQiTUvsVCxXfkMlQkgQewQ9fTC+wCMx6bNtk83SzEZiBaY5M6rm2xRHQadk5HY/+HmC/cbuRHY+3bdwt2OBJCChYiofUVFiqWkMtQ8RPwO4I+GzVNJTZTUd2UifvYEuwj8y+wMNEPqxIUTxViK1B/xXYGqcZVqqNgIRJK87BQ8T25DBVfY92m5+TsHTMwjarLHzUZjwUK8VwRFiKuAv4AbB3scCTkFCxEQudTrFBzBbkMFROxTRJrCMHSR7rFGT5vlZ+DyE/usseR2CzFnsEORyJCu35EQmUOcBi5DBVJYCTWx7M8Z++aoQQwNcPntvFzIPmnECtdeQY73k6hQjKlGQuR0Hgfa361mlzd3suBs7GztkJnATAYa2pVmwJsd0gv30eUF9xlj8uAa7GTQUWyoRkLkVCYgS1/5C5UfAIcQEjP+/g31hThdWxB/7fY3W7zakH3ejTqSdFA7o+yF/ARcAsKFVI/ChYigXsbW8VeS65CxdPYyd9f5+wdM1SOzVKchtVMHIRN5IzBBr3zZs8vdR4fkLshxlEhsBPWAuQ17JAvkfrSUohIoKYCx2Lnc/tfMrkRGAH8hVRhXmi8ia3LzMfudH/AzmR3G1oNwLas1NV5UzLm/uiGAdcDzYIbisSIgoVIYCYBx2NHbfofKr7Hdn24tZChCRUbgD8BN2M/hg7AP7Ea1s0VYcelSIO4ofIQbDJo32CHIzGjYCESiFewbk4byUWomAScge01CU2gAGvXcRbwrnN9DnA3OivbR4VAC6wT+hmoyZV4TzUWIjn3AjZTsQG/Q8V64Aqs/fIPhKieIgk8jB19+S52nsc4rN2nQoUv3LM9LsGOsjsThQrxh2YsRHJqPHAqFij8nTv43HmnD513Cs1MxXLgV8CzzvWRWKBoF9SA4s1d9jgYuB/bbCPiJ81YiOTMOKzKwd9QkQT+DnTF2m2FJlCArQDth4WKrYBbgVdRqPBJEbA9ltveRKFCckMzFiI58U9sH6W/t/mV2GTAM76+Sz38BFyN1U8A7AM8DnQLakDxVoTF198CN2IrTSK5omAh4ru/Axfid6iYhrV/WObru9TDB8Ag4GPn+iJspkKnj/rmAOzY+/2DHojkJS2FiPhqDHABfoaKDcD/YX07lxKiAs0EcDu2uP8x0Ar4L3APChU+KMJmJv4GTEehQoKjGQsR39yFtR7yz1fA6cB7hKxA8zvgXGyfK9jO2oex9o7iKXfZ41fATdhWUpEgacZCxBe34XeoeByrg3yfEAUKsBbbXbBQsTU2J/8sChU+6YLt2B2DQoWEg4KFiOduBK707dXLsZKFQcCPWIutUCgHzsM2vqwEumOp59eoYYLHirB2Hw9is1Xdgx2OSBVZBYuRI0dy0EEH0axZM1q2bMlJJ53E3Llz/RqbSMQksYOmr/XtHaYDnbGNq6HyFrbDYyz2r8o1zmN7BTekOHLP9jgP+AJb/tCnQwmbrP5OTpkyhSFDhjB9+nQmTpzIhg0bOOaYY1i7dq1f4xOJiCRwFTZb4b1KbP38UGARISrQ3ICdXtULOzxsF2Ay9mPYquZvk/rZFwuXDwE7BjwWkZpkVbz58ssvV7keO3YsLVu2ZObMmRx++OGeDkwkOpJYPcXddTyvfhZgZzq86curN8AX2Dkf7zjXZwH3opbcHivCSlVGYatKOsxVwq5Bu0LKysoAaNGi5pKhiooKKioqNl2Xl5c35C1FQiYBDMGaJXvvaeB8rJYiNNzWnkOBtViQuB/bniKeKcJmps7GQkXLYIcjkrF6L88lEgmGDRvGoYceSufOnWt83siRIykpKdn01a6devdKXCSwz5Deh4q1WPeLU4A1hGjp4wfgZKzf11rsCPMPUajwQUfgDeARFCokWuodLIYMGcKcOXN48skna33eiBEjKCsr2/S1YMGC+r6lSIhUYi26/+75K8/CthCOda5Ds5X0f9j+1glY/cQo7JyP9kEOKl6KgG2wRbXZWE2NSNTUaynkoosu4sUXX2Tq1KmUlpbW+tzi4mKKi4vrNTiRcNqAFRQ8hZe3/QRwBzDCeVV/D1TPwjrsnI+7nOuOwGNY32jxhLvscTrwF6B1sMMRaZCsgkUymeTiiy9mwoQJTJ48mV133dWvcYmE1HrsRI7n8DJULMb6Urzu2St65ENsYHOc6yHYOR9bBzai2CkA9sT6iKkEXuIgq2AxZMgQHn/8cZ577jmaNWvGkiVLACgpKaFpUzX/l7irwAoMXsLLUPECcA5WSxEaCWA0Nn2yHlvkfwQ4LsAxxUwhUIxtI74I7c6V+ChIJpMZ/wtZUFB9+7xHHnmEc889N6PXKC8vp6SkhLKyMpo3b57pW4sE7CfgRKxPtTeLFD8BlwN/xT61hqaWYiF2zserzvXx2MlWqiBssELsb09LbIPyr1A/CvHLZOxoQu9kev/OeilEJP+sxe6uU/EqVHyE7fiY51yH5r+sZ7CNLiuwE0jvRC25PdAIa72+PxYmT0YzFBJfOt1UpFbzsQgwGy9CRRLrIXWZ82qhKdBcDVxCaivKgViB5t5BDSge3GZWJ2MzFD8LbigSW+7fMndT+tZYf/2OgYwGFCxEavEk1rBhHV50kliGbVB9ua4n5trb2CaXr7CZiauBPwKNAxxThLnLWiXA752v2vfOiWSiALtlb3Cui4Ddse1Z+6V97ULQU4wKFiJbWANcjH1896b64RVsc8WqBr9SPVUC07DtJ22wsz2S2JkeNzq/3h74J9qaUE/ultE9seWOQWjzjNTXVtjimftvT2tsFqIrqQDRkbCmfwULkSpmYUsfXzvXDQsVFdjGijtJFe7l3His/fZ3aY+1BpqRKvIYhK3RbJfTkcWCGyiOAS4F+hD050WJjs2XMbbFQkM3UgGiM1H7D1PBQgSwW/5dwJV41Z7qMyyifJL2Djk3HhjIlvloifO1NfAwdsqZZKUA2y56AVaeohPipWYFWIjY6Fw3wua2Nl/GKCUOsVTBQoRl2FFP//Pk1ZLYsdaXYP+MBFagWYnNVNQ26VICnJqb4cSBOztRis1OnE/UPkuK/7YiVQcBsDNbLmPsRZz3BSlYSJ77H3AmXlU//IB9gn3Ok1droGlUXf6ozmLneUf4PppIc7eL/gzb0fMLdHy5NKLq3q7m2Ck/3ai6jNEsiMEFSsFC8tR64BrsZAZvqh9ex1YUljf4lTyyMMPnLfZ1FJFW6HydiW0X3T/Q0Ugw3L8F7jLGVtg+7PRljC5Y4VL0lzG8oGAheWgeNv//gXPdsFCxAbgOO+yzoMGv5oFKYBxWNZqJNj6OJYLcmNkCW876LdAq0BFJ7rhzU2D/NbfH4mQXbPZhP2APdOusnX46kmf+CfwGiwMN30b6BXYk2fvOqwXaQbMS+DfwJ6xyFGrfLVuAFQv08n9oUeDeUvbFtouehhVnSpy5VTNgPSGOJRUkOmGH2Eu2FCwkT5Rjnz2f8OTVksA/gN/hVURpgAR2gvsNwKfOY9tjxQC7YCecQdVBujO2o8n7YoEi7Ed4PFaQ2QtNaMdXAal5xRKgr/N1NNAuwHHFi4KF5IF3sY2fmRYd1G4VNufxb09erQESwNNYoHD3tG4HDMfm8Eucx7Zmyz4WpVioGJCDcYaQGxy2wY5CuRjoENhoxF/uXFQj4BBsVuIYbGaiMMBxxZeChcRYArgNK9IEL9pyv4VVZyxp8Cs1QALrT3EDMMd5rAQLFENJBQrXAOxg1s07b+bhTIU78d0Bm9AZjLUkkjhJX97YCzgOCxKHo6WN3FCwkJhajLWTfN2TV9sI3ITdywvxIqLUQwJ41hnEh85jzbH5+2HU3lChiLzeUup+Zu2N5a9+6LNqfKQvb2yPzUj0xXqg7hzguPKXgoXE0H+whlerPXm1r7FtpO9gZQo5DxVJUoHC3cjSDAsTl2L/lkq1CrF/5AZjkzn7Bjsc8YwbFbcCDiO1vNEFRcbgKVhIjFRgLbnvxovDwyqBe7Bdm+nHAeVMEngeO2l0tvNYM+wOeSm2H1K24G4XbYVlr18BOwQ5IPFA+jbQfUgtb/QCmgY1KKmBgoXExGfYoRjutoiGxYBZWLvmD+p6oh+SwAtYoHjfeWxbrCBzOLpL1sC99RyA1U+cTJybJsdd+vLGjlRd3mgd4LgkEwoWEnFJ4O/ARdjGz4a1p1oDXIsdR5bzCdUktorzR2Cm89g2WKC4DAWKGrg1qAOxGYoewQ1FGsSNho2xQkt3eaMz2gAcLQoWEmGrgAuBZzx5tRexbaRLyHEtRRJ4CQsUM5zHtsGy0uXYBzapwl3oKgGGAL9HZXrRk7680ZnU8sahQJOgBiUeULCQiPJu4+cirI/BeLw6NSRDSeBlLFC86zy2NalAsVOuBhId7kbCvbFJnEFohT063DnABNCSVJA4yrmWuFCwkIipBG7G7sYFNGReIQHcj5V7rkt7zHdJ7FDV67GtJmB3xyHAFejf2M24n2vbA2dhrbb3Q5Pj0eD+6TXB9jv3xcLEPuhPML4ULCRCvsM2fr5JQ4szP8KON59R1xO9lARexQLF285jTbG+4Feik67SuLejNtisxGnAgehWFH7un1wB0BXrGHIM0BOdvJI/FCwkIp4FzgXW0pBQ8SPWDuJ2L4aUqSQwCZtkedN5rAmpQKEidyB1S9oRi4+nAz9DXQnCrZDU8XutseWNvsDPUXFQ/lKwkJD7CdtjeT8N7U3xP6zUcyE5XPJ4HZuheMN5rBg7C+0qdFw5qTCxHTYrcTp52208QjZf3nBnJfZGc0oCChYSanOwTYTznOv6hYql2DbEJ8lhcebr2AzFVOe6GNtychXQNhcDCC/3ttQM+9M9Hft8q3+Mwir97I1O2DmwfbHdG1rekC3pv2UJoSQ2QzEM+wetflEggXW4uAxbQHEf89UUbIZiinPdGDs+82ryej+kGya2Bn6JhYljsB+PhE0BqRNxtsdmJI7FjhbXup3UTcFCQuYHrOfl8w16lU+xZY+3PBhRRqZhgcI986wx1kv6auyI8jzkfs4tBn6BhYl+aHtoOKUfLd6T1FbQbqjKRbKlYCEhMgVbaV9e71dYh21GHYkHZ3tUUvdR429gSx6TnOutsEQzAmjX0AFEjxsmtsJCxJnYxLkOqw6bImz+LgnsQmp54whskUqk/hQsJAQ2Yns1biJ1PkD2Xsfu6fPxIFSMxw77+i7tsVKs1/cAbCrkemz7KNid9HzgD1jDhTzi7gsoxCbLzwROxE50lzBJX5DqQ6rocrcgByUxpGAhAfsamyR/l9S2tewsx+oo/kHqJtcg47Gqws1faCF2slVXUqeTNSIVKHZp6BtHR3rtf2+s18Qv0YGr4ZLeU6Ib0B+bleiBjmcTPylYSID+jbWp+on6xIEk8E9sYmG181iDizMrnResbjjuYx9gM8nnAdcAHRr6ptHgnjeZxPYDnInlLDUKDYv0E0HdltnuiaDqKSG5o2AhAViLHdn5d+rbm2Iettlicr1foQbTqLr8UZN/Yl2c8oC7Rbc7NjNxCnm/YzZE3FmJrbATQfthYWJf1FNCgqJgITk2G1tnmO9cZxcJ1gO3An8iNTvhWagAK9SUTWGiKxYmTiWvVnpCLL2nxF6kljcOR/ttJCwULCRHksA9WDVEkvosWryBLZzMw+MwkS7TE0Vj2DXTvWXtQ+p8jj0CHZEY90+mORYijsWKLvN0H7OEnoKF5MD3wGDgpXp990qsYeVD2D+xvoWKV7D6itoUYP+e9/JrELnl3rJ2JxUmOgU6IkktbxQCB5GalTgQNTuXKFCwEJ9NwooRVmT9nUlgHHARsMp5rP6HpNdiHjaR8oJz3QyrBt28eMNdsh5NpP99d29b7bAwcTrQBa3IByf9IK+dsSBxLHAkdoqKSLQoWIhP1gPXArdRn94U87GjNSbicXFmunLgRiwobMD+a7gYuA54jer7WIzG+lhEjBsmWmO7OU7HijEVJnKtCPupb3Sum2AnpRyLzUrsif5UJOoULMQHH2O3ro/JtjfFBuBO7N7u/tPreahIAI9i3TGXOo8d67xxR+d6ANblqa7OmyHm3p6KsG2hF2K3MDVozgX3L4o7x9YYCw1dscUm92t39M+wxI3+RouHEsDdwJWk2gVn7h2s19SnWX9nFt7CdrrOdK73wgLFcdU8twjrcBwx7uxEJ2zWZxBqXOWfRtjfdXdGrgl2fPjmAaIDkUqlIg2gYCEeWQCcTepYz8yVY5MHY/Coc2Z1vsMqQB93rptj0yIXE4sjNt3todsCZ2GzEwegSXXvbB4gtsb2z3TBgsO+zv+2Q3NCku8ULMQDT2CfjX/K6ruSwATgd1hb7iQ+FGf+BNyOnUr2I3anvQCrrWjl9Zvlnrur41CsYdjJqJtBw2yFzfe48bYZFhjcAOF+7Yxim0j1FCykAVYCv8Vac2dXYrkA+D3wYtbfmaEkdubH5dhxJGB337uwXXsR5oaJVtjMxHnYSr1kY/MAUQJ0BvYjNfvQCfspK0CIZEPBQurpVWzS3T3iPLNoUIm1yRqBFWpm/p1Z+BDb0THZuS7FNqecRqTvEYXY8E8AfoXtIdCqfV22IvU3DWAHUgEifQYi085oIlIXBQvJ0k/A1ViRpruyn5lZWHHmB3U9sb6WYztcH8SG1QSrI70S2MavN/WXW4i5J7bYdDY69Kt6mweIVlh46EwqPOyDylhF/KdgIVl4H9tG+oVznVmoWIPVSY7Gp7K2DVjl5/WkOmmdih0qEsEDLty41hRrLXYh8DMiPdnisfRA25bqA0RJMEMTEQULyUQldpe+jmzP+XgR+6S9BJ+KMycCw4BPnOuuWB1Fb6/fyH9u7cRB2M/sFGyXh7jzNgXAwViDkaOwEKGfkEjYKFhIHeZjvRqnZ/Vdy7DizGfIdsEkQ19gbbifd653BG7CdnxEqPDADRMtsKGfT6pHV35zw0RzrMX18VhVyQ5BDkpEMqBgITVIAmOBIVh77sw9hW19XO1cexoqVmMB4k5nWI2ww0SuA7b38o385fbr6IsVYvbHqgTyV/p5GZ2wWYn+2CJQhJKiiGS/5D116lROOOEE2rZtS0FBAc8++6wPw5JgfQ+chH1+/olMFzC+BwZi5Q1lGX9XhhJYztkLGIWFimOwHSB3EolQ4ab4Dlgbje+A/2A/6fwMFe5PpAnQDyuU+RZrBX8ztj9YoUIkarKesVi7di1du3bl/PPPZ8CACJ7GJHX4D3bE+aqsvusZ7JN3uXPt6RbS6Vgb7hnO9R5YmOhP6Csa3eE1xgLXhdiRIyEftk8KSB1IV4pFqv5Y3/QmgY1KRLyVdbDo168f/fr182MsEqi1wHBsr2bmLauWY7UUT2X1XRlaiO1s/Zdz3QzbTnoJUOzlG3nPrZ3ohhVink6+7lNwayWKgEOAX2BhoiP5Gq9E4k41FoId/3U6Ng0NmcaDCdgn8LKsvisD64A7sNnwtdj95zystqK1V2/iPTdMlGBzPhdiGyHzj/uT2B5r53U8tm6Vn9FKJN/4HiwqKiqoqKjYdF1eXl7LsyW3NgB/xu7Y7hR13X7ASjrH4fEshXt4yGWk2nAfgm0f7e7Vm3jPLTs8AitaPZHQT6h4zA0SYPt93cLL7uhALpH843uwGDlyJDfccIPfbyNZm4u1X5pNNtHgOWxb5Crn2rNQ8RHWj+I153pnrHXGGYRyxtyd4N8Zqy05l0j24moA9yfQFNvbcgJWgNkmyEGJSAj4/nFixIgRlJWVbfpasGCB328ptUoCf8U+WX5IptFgBdbN4iTn/3u248Od/uiGhYpi4P+w3HMmoQoVbulhI+wU0YnY4tH15EOocH/3ALti582/ikXMCdgOIoUKEcnBjEVxcTHFxfk1MRxei7HV/4lZfdfz2CzFSufak1mKjcD9WP8J94UHYoeFdfDiDbzjfjbvhBViDiJfTpxwf+eNgMNJFV7uEeSgRCTksg4Wa9as4Ysvvth0PX/+fGbPnk2LFi1o3769p4MTL43HPlWuyfg7VmKfSx/D41qKV7Flj4+d6y5YHcURXr1Bw7m31BZYkBgMHECoJlB84tZL7IQFieOBPqh1tohkKutg8d5773HkkUduuh4+fDgAgwcPZuzYsZ4NTLxShu3P/AfZxIMXsRiywrnOKlRUAtOwCZI2WOOGIuBLrDDzOed5O2Cdoi4kFPuT3FtqE2AAdpJoH0IxNB+5v+sC4EBShZfdyIcYJSLey/rfzCOOOIJk0tNuBeKbqVihwhLnuu4/t1VYDPkn9TzjYzwwFGsr6WqLdWZ+EeuYWYTVVfyRwDtmpt86f47NTJxI3D+fu/Mx2wLHYbMSx2KzFCIiDRPvD2N5qwLrJPUXstlG+l+sXcQPznW9QsVAtswvi5xfAzgaOz+9U7Yv7q30BlbnAqcBrQIcj7/Sz+HYEyvBPR7by6t/AkTEW/pXJXbmYM2uPiF1M6ndKuBS7CiOetdSVGIzFbV9845Yegnob537Ob0dFiYGAXsHM5QccH+3WwFHkiq87BDgmEQkHyhYxEYCmwq4ikwDBcDL2E12uXNd70WuaVRd/qjOcuANclqk6d5eS7BFobOxVZl4Vg+48zCtsAWd47EFnm2CHJSI5BkFi1j4FrtlTs34O8qwk0H+Tj1rKTZ/sYcyfO7ihrxRZtzfT2Ns0v9srKF0Y//fOsfSCy8PIlV42YW4RicRCT8Fi0hLAo8Dv8UO2MjM/7Aixe+d63qHim+xbaIPAasz/B4feyi5lQS9gXOwnR3N/Xu7gLhzMM2oWni5Y5CDEhHZRMEislZggeKpjL+jHJul+BsNnKWYCdwO/JtUC859sM0nq6h+PaUAOym7V33ftHruZ/bO2JLO6Vib7fjYvPDyl9ishAovRSSc9C9TJE3EJviX1/XEKt8xGFjmXGcdKhJY4eXtwOS0x48CLseOi5iA7QrZvALUnZUfjSWBBnI/s7fFfk+DgH0b/rIhkl54eQSpJY4OwQ1JRCRDChaR8hNWnHkPmc45rMZ6Uj2U8XdsZh3W1OIO4DPnsUbY1MBl2H5N1wDgabbsY1GKhYoB2b55ijsz0cx567OAw4jT2ZnphZe/wA71UuGliESPgkVkzMJuqV8613VHhEnYJ3q3XjKrUPE9dlbZfaSKMZpjh2VcgoWF6gzAPmBX13kzS+4iQBFWSXAOVlUQj5Nn0gsvu5PaxaHCSxGJNgWL0KsERmGndUGmsxRXAA9Qj1mKucCdwKOk6kHbY40uLsCmDOpSRIO2lLpjPgQLRicTeINOj6QXXvZDHS9FJI4ULELtK6z7wrtk2mHiNexmvMi5zihUJLEZhtuBF9LeqjtWP3Eyvv9NcW+5e2PdP8/A8ky0qeOliOQf/esWSkmsw8TF2OEadYeKNcCVwBiymKXYCDyDBYoZaY//Aquf6IWvs/JumGiJLXOcRRwWAjYvvHQ7Xu4a4JhERHJHwSJ0yrA5h+fqeuImk7Eb80Lnus5QsRrbczoa+MZ5rInztpfia59rt7Jga+BUbG9LbzzZLBKgzQsvj8e2y6jwUkTyj4JFqMzDyhPnZ/TstdgekfvIcJbiO2xDyQNYfgFb3h8C/B7flvrdBYFCrLLgHOzW29Sft8uB6gov+wNdifp8i4hIQylYhMZErJjhR1Jdp2o2Ffu07+7qrDVUfIAtdzyBzdKDzUoMd17Epzu8e/s9CJsMOYUo94fUUeMiIplQsAhcEuuLPTztumZrgRFk0MkiCbyCBYpX0x7vjRVkHocvTSDc2++uWBHmIGA3798mB9zfCVhb0eNJdbzcKqhBiYiEnoJFoCqwxhCPZvTsadgEwwLnutpQUYEdH3I78LHzWBE2XXAZNnPvMXdmojlWgHkOcDBRWxRI38HREktefZyvVgGOS0QkWhQsArMEK/SbWeczfwT+ANxNLbMUK4D7samMJc5j2wK/wjph7tLQ8Vbldu1uhPWIHIzVT0TrBFE3EjXFulweAxwNdCRqsUhEJCwULALxHja1vpy6Si7fxGYpvsFu5FtUX3yJNbR6BEsgYKdwDcVCxXbejNjl3ooPxpY6TiVKzavc5Y1CbOqmHzYj0QMtb4iIeEPBIucex27JldRWpPkTcA22I7SAauLHW9hyxwRSZRndsPqJU/H0PunejjtgIz+LqNRNuDEIbMTu8sYRQElAYxIRiTcFi5ypxBY0bq35GZUwbRpMWQwPtIElvSBZlFbOWQk8iwWKt9O+8TisfuJIPJvBT6+bOBOrm/iZdy/vkwJSMWx77MjVo7EwEf0+niIiUaBgkRNlWJPql2t8xvjxcMlQWLj5qaB3YffHR7Alj6+cX2uMrZEMBzp5M8r0Q7/6A+dimSXcdRPufEpjrFXosViQ6EKczj4VEYkKBQvfpTe9qn4r6fjxcPLAan55IdbaYhtsnylAC6yZ1RCgtTcjdGcnupOqm2jhzUv7wA0SBVh4cOskDsXah4qISJAULHz1P2AgtTW9Wl4JZw+l+szhPrYW2B2bnRiMJ52i3dtze1J1E3s0/GV9UIQtbSSxqtR+2O6NI4lyuy0RkbhSsPBFEiu7vCzteksTgHOnwY/fVfvLVT2AHT/RAOl1E2dgdRM9CVvdRAG2hFGJHS/eBwsSfbB0Fa7RiohIVQoWnqu76dUSbCVjPMDiDF92Wf1G496Gi7AFmcFY/URx/V7OJ+78SSOsRNStk+hO1I8nExHJNwoWnqq96VUS+AdwCU7JxEZgeoYv3Sa7kbizEwcC52N1Eztk9xI+Sm+XvTcWeY7Gii+3DWpQIiLiAQULz9Te9Oob4ELSju14CVsp+bSOly3Adof0qnsE7u26Ham6iT3r/rYccHdnJLBDu/phQeIosk5MIiISagoWnqi56VUC+CtwBbAB7PyOy7ADwsDqD38JPOxcp5djuOsYo6lxRcCdmWgGnI7VTRxK0JUI6XUSTbBCy77Y8kYngh6diIj4R8GiQWpvevUZFjemA3wPXAc8iKWNrbC229dgbbePda4372MxGhhQ9XXd23Ih9tn/XKxuIrjNlukzEk2Aw7Aw0Rs7ND3cnTBERMQ7Chb1VnPTqw3AbcD1QKICOz3sRqDcecIALIvsnvZNA4ATsSNMF2MrBL2oMlPhzk4cgAWW0whqw2V6kNgGG6gbJA5A526IiOQvBYt6qbnp1Sxs58WcJLbt40pS3TL3x7pn9q7hZYuwYyzSuHUTpdjMxNnAXg0cffbSg0Qz4HDsNNDe2AEl2rkhIiJGwSJr1Te9+gm4AZupYCbWzGqq84ttgJuxAogMuky7MxPbUrVuIncNqtObUpVgsxFHYEFiPxQkRESkJgoWGau56dU0bDZh/iJIXoO1sEhi5QZXYLMWGeyidGcnegO/BU4gV3UT6UGiBTYbcYQzkk7ozA0REcmUgkVGqm96tRq4ChjzIxTcDslRpM70GASMxPZ+1sItxNwa6zcxBOvs4K/0PhI7YUHCrZHYG+3aEBGR+lKwqNNirKqyatOrl4DzE7D0CeBqSLq7OXpidRQ9an9V99a+FzAM6znhX2uo9CDRGusfcQQWJPZAQUJERLyiYFGrGdiCRKrp1Q/YrtDH3gIuBd51ntoeGIVt1ajlPl2ELTicgHXg7F370+spPUiUYv0jejtfHXx5RxEREVCwqEXVpldJ4CngV9/A6quAcc7TtgVGYCGjafWvVIjFkhbA77BFlTpWSLKUHiR2wWYk3KUNb99JRESkNgoWW9iy6dUi4MLV8NJI4A6s5KIAK4q4EVtdqIZ7uz8AW+4YiFeHf22F08cTa4aRPiPR1pN3EBERqQ8FiyqqNr1KAg9VwiWPQMX/AUudpx2JBYxu1b9KIbbkcSZwEXZGZ8NsfmiXGyQOB1o1+NVFRES8omCxSdWmV18Bv3wNPhwOfOA8ZQ/gL9gBppuVKbi9J9pitRMX0JCumG6QKAD2oWqQCKbXpoiISCYULID0pleVVHLtPBh1OSSed355O+ycjyFsceyFGyh6Y0Wd/alP+6j0GYmO2MEhR2KtsrfP+tVERESCkufBIr3pVQFvrkww4E+w7F7sPl+EVVteT5WJAneyoimp3hMds3rf9IZUbbEg0QfrJ6GlDRERia48DhappldrN8Ap9yd56Y/ACueXj8OWPfZJfYc7r7AHVox5NnZyRt0KnK8E1iL7aCxI9AF2Q9s/RUQkLvI0WFjTq2TyPe78L4y4HNZ/5vzSvlhh5jGpZ7vzC8dj9RNHkEkUcGNIMbakcQwWJLqiFtkiIhJXeRgsrOnVjDnLGDg8ybcTnYd3BP4MXAg0SvWe2A74PTa30b7W13WDRCFwILa8cRTwM7zaZCoiIhJ2eRYsHmfZsnO58LoNvPAQlhwaY1WX1wAlqXjQDVvuOIWaDgJLL7jcC+iHBYneQHPffgciIiJhlifBopKKiisYdfed3HgjbCh3Hj4Za8O9e6r3xBlY74mDtniN9ILL1lQtuGzj++9AREQkCmIRLCorK5k2bRqLFy+mTZs29OrVi6Ii2/SZTK5i/PgjuejK2Sz5yvmGA4A7oKh3qvfExVjviZ02vWoBFjcqsRmIPljR5VHo4C4REZHq1StY3Hfffdx2220sWbKErl27cs8993DwwQd7PbaMjB8/nqFDh/Ldd99teqy0tJS77rqLXXaBi4efydtTK+wX2gA3Q9E5UFloJZVDsaJM+0G4yxuNgcOAvliQ6EZ9ulOIiIjkm6yDxbhx4xg+fDj3338/PXr0YPTo0fTt25e5c+fSsmVLP8ZYo/HjxzNw4ECSyWSVxxcuXMjJJ5+ceqApcDlwJTTdNtV7Yp8qBZcHkAoSPampskJERERqVpDc/K5chx49enDQQQdx7733ApBIJGjXrh0XX3wxV199dZ3fX15eTklJCWVlZTRvXv8ix8rKSjp06FBlpqJaZwCjYI92MIxCzibhlFbuSapOojfWX0JERESqk+n9O6sZi/Xr1zNz5kxGjBix6bHCwkL69OnD22+/Xe33VFRUUFFRUWVgXpg2bVrdoQI45Nfwp3bwc1pSQD9SBZc6BVRERMRrWXVqWr58OZWVlbRqVbXtdKtWrViyZEm13zNy5EhKSko2fbVr167+o02zePHijJ530eLzOIrPKGAJMBY4C4UKERERf/jeAnLEiBGUlZVt+lqwYIEnr9umTWZbPNu0OQc7aly7OERERPyW1VLIjjvuSFFREUuXLq3y+NKlS2ndunW131NcXExxsfedJ3v16kVpaSkLFy7congToKCggNLSUnr16uX5e4uIiEj1spqxaNy4MQceeCCTJk3a9FgikWDSpEn07NnT88HVpqioiLvuuguwEJHOvR49evSmfhYiIiLiv6yXQoYPH85DDz3Eo48+yqeffsrvfvc71q5dy3nnnefH+Go1YMAAnn76aXbeeecqj5eWlvL0008zYMCAnI9JREQkn2Xdx+K0007j+++/57rrrmPJkiV069aNl19+eYuCzlwZMGAAJ554Yo2dN0VERCR3su5j0VBe9bEQERGR3Mn0/u37rhARERHJHwoWIiIi4hkFCxEREfGMgoWIiIh4RsFCREREPKNgISIiIp5RsBARERHPKFiIiIiIZxQsRERExDNZt/RuKLfRZ3l5ea7fWkREROrJvW/X1bA758Fi9erVALRr1y7Xby0iIiINtHr1akpKSmr89ZyfFZJIJFi0aBHNmjXb4rjzhigvL6ddu3YsWLBAZ5D4SD/n3NHPOjf0c84N/Zxzw8+fczKZZPXq1bRt25bCwporKXI+Y1FYWEhpaalvr9+8eXP9pc0B/ZxzRz/r3NDPOTf0c84Nv37Otc1UuFS8KSIiIp5RsBARERHPxCZYFBcXc/3111NcXBz0UGJNP+fc0c86N/Rzzg39nHMjDD/nnBdvioiISHzFZsZCREREgqdgISIiIp5RsBARERHPKFiIiIiIZ2ITLO677z46dOhAkyZN6NGjB++++27QQ4qVkSNHctBBB9GsWTNatmzJSSedxNy5c4MeVuzdcsstFBQUMGzYsKCHEjsLFy7krLPOYocddqBp06bst99+vPfee0EPK1YqKyu59tpr2XXXXWnatCm77747f/7zn+s8a0LqNnXqVE444QTatm1LQUEBzz77bJVfTyaTXHfddbRp04amTZvSp08f5s2bl5OxxSJYjBs3juHDh3P99dcza9YsunbtSt++fVm2bFnQQ4uNKVOmMGTIEKZPn87EiRPZsGEDxxxzDGvXrg16aLE1Y8YMHnjgAbp06RL0UGJn5cqVHHrooWy11Va89NJLfPLJJ9x+++1sv/32QQ8tVkaNGsWYMWO49957+fTTTxk1ahS33nor99xzT9BDi7y1a9fStWtX7rvvvmp//dZbb+Xuu+/m/vvv55133mGbbbahb9++rFu3zv/BJWPg4IMPTg4ZMmTTdWVlZbJt27bJkSNHBjiqeFu2bFkSSE6ZMiXoocTS6tWrk3vuuWdy4sSJyd69eyeHDh0a9JBi5aqrrkoedthhQQ8j9vr37588//zzqzw2YMCA5KBBgwIaUTwByQkTJmy6TiQSydatWydvu+22TY+tWrUqWVxcnHziiSd8H0/kZyzWr1/PzJkz6dOnz6bHCgsL6dOnD2+//XaAI4u3srIyAFq0aBHwSOJpyJAh9O/fv8rfa/HO888/T/fu3TnllFNo2bIl+++/Pw899FDQw4qdQw45hEmTJvH5558D8MEHH/DGG2/Qr1+/gEcWb/Pnz2fJkiVV/v0oKSmhR48eObkv5vwQMq8tX76cyspKWrVqVeXxVq1a8dlnnwU0qnhLJBIMGzaMQw89lM6dOwc9nNh58sknmTVrFjNmzAh6KLH11VdfMWbMGIYPH84f/vAHZsyYwSWXXELjxo0ZPHhw0MOLjauvvpry8nI6duxIUVERlZWV3HTTTQwaNCjoocXakiVLAKq9L7q/5qfIBwvJvSFDhjBnzhzeeOONoIcSOwsWLGDo0KFMnDiRJk2aBD2c2EokEnTv3p2bb74ZgP333585c+Zw//33K1h46N///jePPfYYjz/+OPvuuy+zZ89m2LBhtG3bVj/nGIv8UsiOO+5IUVERS5curfL40qVLad26dUCjiq+LLrqIF198kddff53S0tKghxM7M2fOZNmyZRxwwAE0atSIRo0aMWXKFO6++24aNWpEZWVl0EOMhTZt2tCpU6cqj+2zzz58++23AY0onq644gquvvpqTj/9dPbbbz/OPvtsLr30UkaOHBn00GLNvfcFdV+MfLBo3LgxBx54IJMmTdr0WCKRYNKkSfTs2TPAkcVLMpnkoosuYsKECbz22mvsuuuuQQ8plo466ig++ugjZs+evemre/fuDBo0iNmzZ1NUVBT0EGPh0EMP3WK79Oeff84uu+wS0Iji6ccff6SwsOptpqioiEQiEdCI8sOuu+5K69atq9wXy8vLeeedd3JyX4zFUsjw4cMZPHgw3bt35+CDD2b06NGsXbuW8847L+ihxcaQIUN4/PHHee6552jWrNmmdbqSkhKaNm0a8Ojio1mzZlvUrWyzzTbssMMOqmfx0KWXXsohhxzCzTffzKmnnsq7777Lgw8+yIMPPhj00GLlhBNO4KabbqJ9+/bsu+++vP/++9xxxx2cf/75QQ8t8tasWcMXX3yx6Xr+/PnMnj2bFi1a0L59e4YNG8aNN97Innvuya677sq1115L27ZtOemkk/wfnO/7TnLknnvuSbZv3z7ZuHHj5MEHH5ycPn160EOKFaDar0ceeSToocWetpv644UXXkh27tw5WVxcnOzYsWPywQcfDHpIsVNeXp4cOnRosn379skmTZokd9ttt+Q111yTrKioCHpokff6669X+2/y4MGDk8mkbTm99tprk61atUoWFxcnjzrqqOTcuXNzMjYdmy4iIiKeiXyNhYiIiISHgoWIiIh4RsFCREREPKNgISIiIp5RsBARERHPKFiIiIiIZxQsRERExDMKFiIiIuIZBQsRERHxjIKFiIiIeEbBQkRERDyjYCEiIiKe+X9wcZMyR2NZQQAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"green = r.gROOT.GetColor(r.kGreen).AsHexString()\n",
"yellow = r.gROOT.GetColor(r.kYellow).AsHexString()\n",
"plt.plot(np.linspace(0, 10, 10), np.linspace(0, 10, 10)*0.4*np.random.normal(1, 0.05, 10), color='k', marker='o')\n",
"plt.fill_between(np.linspace(0, 10, 10), \n",
" np.linspace(0, 10, 10)*0.6*np.random.normal(1, 0.05, 10), \n",
" np.linspace(0, 10, 10)*0.2*np.random.normal(1, 0.05, 10), color=yellow)\n",
"plt.fill_between(np.linspace(0, 10, 10), \n",
" np.linspace(0, 10, 10)*0.5*np.random.normal(1, 0.05, 10), \n",
" np.linspace(0, 10, 10)*0.3*np.random.normal(1, 0.05, 10), color=green)"
]
},
{
"cell_type": "code",
"execution_count": 53,
"id": "f9fa48de-e27f-4340-a9f1-08e5c2b695c8",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PolyCollection at 0x7fdc61f9ff50>"
]
},
"execution_count": 53,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGdCAYAAABO2DpVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAABfGklEQVR4nO3deXycVdn/8c89S/a1SZOmTbqyryLbj7XQlqUCD4vgBoKAqAiyFBBwQ0QoUEFWUcAH0UdAtgIuiAiyCciuIDttoVknaZbJOst9n98fM5OFpm2STuaezHzfr1demGmSuRrTzrfnXOc6ljHGICIiIpIEHrcLEBERkcyhYCEiIiJJo2AhIiIiSaNgISIiIkmjYCEiIiJJo2AhIiIiSaNgISIiIkmjYCEiIiJJ40v1EzqOQ2NjI8XFxViWleqnFxERkQkwxtDd3c3MmTPxeDa8LpHyYNHY2EhdXV2qn1ZERESSYO3atdTW1m7w11MeLIqLi4FYYSUlJal+ehEREZmAYDBIXV3d4Ov4hqQ8WCS2P0pKShQsREREpphNtTGoeVNERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREQkU/S/DE2ngd3lWgkKFiIiIlOd3QHNp8PHe0LX7RB6w7VSfK49s4iIiGweYyD4f9ByDjhdgHG7IgULERGRKSn0NjR/E/qfAyzSIVSAgoWIiMjU4vRC20+g/ZphD6ZHqIAJ9Fg0NDRwwgknUFFRQX5+PjvuuCOvvPLKZNQmIiIiCcZA90Owaito/xlgx9/Sy7hWLDo6Othnn3048MADefTRR5k+fToffPAB5eXlk1WfiIiIhFdDy5nQ+xdiawKO2xVt0LiCxVVXXUVdXR133HHH4GPz5s1LelEiIiICOKHY6kTbTxgKE+kbKmCcWyGPPPIIu+22G8cddxxVVVXssssu3HbbbRv9nFAoRDAYHPEmIiIim9D7BKzeHtp+CISBqNsVjcm4gsWqVau45ZZb2HLLLXnsscc4/fTTOeuss7jzzjs3+DnLly+ntLR08K2urm6zixYREclY0SZo+BKsXQKR1aRTY+ZYWMaYMVeck5PDbrvtxvPPPz/42FlnncXLL7/MCy+8MOrnhEIhQqHQ4PvBYJC6ujq6urooKSnZjNJFREQyiIlCxy3QejGYATarMXP2U1CwMFmVAbHX79LS0k2+fo+rx6KmpobttttuxGPbbrstDzzwwAY/Jzc3l9zc3PE8jYiISHbp/xc0nwahN92uZLONK1jss88+vPfeeyMee//995kzZ05SixIREckKdjsELoau28iUWzbGFSzOPfdc9t57b6644gq+8IUv8NJLL3Hrrbdy6623TlZ9IiIimccY6LoTAueC002sjyL9ZlJMxLji0e67787KlSu5++672WGHHbjsssu47rrrOP744yerPhERkcwSegs+3geaTwank0wJFAnjHul9+OGHc/jhh09GLSIiIpnL6YHWH0PHz4nd7ZGZdFeIiIjIZDIGuh+MTc60A6T7gKvNpWAhIiIyWcIfQfO3oe9vpNMNpJNJwUJERCTZnAFovxrafsrQCkXmhwpQsBAREUmu3seh+ZsQWUO2hInhFCxERESSIdIAgXOg+35ihy6zL1SAgoWIiMjmMVHouBFavw8mHH8wsxs0N0bBQkREZKL6no+N4g6/7XYlaUPBQkREZLzsdRC4ELp+DXjdriatKFiIiIiMlXGg6w4InBcbeAVk2uTMzaVgISIiMhYD/4ltewy85HYlaU3BQkREZGPsbmj7EXTcQCaP4k4WBQsREZHRGAPd90HLd8BuI5tPeoyHgoWIiMinhT+A5tOh7wmyZRR3sihYiIiIJDj9sO5KWHcFQ2FCoWI8FCxEREQAev4Kzd+C6CcoTEycgoWIiEjnHdB8Ctk8ijtZPG4XICIi4qr+V2KXhgFq0Nx8ChYiIpK9om3QcCQKFMmjYCEiItnJ2ND4BYi2oOmZyaMeCxERyU6t34e+p1BPRXIpWIiISPbpXgntV7ldRUbSVoiIiGSX0LvQeAIazz05FCxERCR72N1QfwSYENoCmRwKFiIikh2MgaaTILIaNWtOHvVYiIhIdmj/GfSsdLuKjKcVCxERyXy9T0LrhW5XkRUULEREJLNFPoGGz6NmzdRQsBARkczlDED9UeB0o+maqaEeCxERyVwt34HQv1GoSB0FCxERyUydv4au292uIutoK0RERDJP/8vQfLrbVWQlBQsREcks0Vao142lblGwEBGRzGGisRtL7QAaguUO9ViIiEjmaP0+9D2NxnW7R8FCREQyQ/ABaL/a7SqynrZCRERk6gu9A01fRUOw3KdgISIiU5vdDfX/AyaMtkDcp2AhIiJTl24sTTvqsRARkamrfYVuLE0zWrEQEZGpqfcJaL3I7SrkUxQsRERk6ol8Ag3HombN9KNgISIiU4szEJus6fSg6ZrpRz0WIiIytbScCaH/oFCRnsa1YvHjH/8Yy7JGvG2zzTaTVZuIiMhInbdB169RqEhf416x2H777fn73/8+9AV8WvQQEZEU6H8Jmr/tdhWyCeNOBT6fjxkzZkxGLSIiIqOLtkL9UWgAVvobd/PmBx98wMyZM5k/fz7HH388n3zyyUY/PhQKEQwGR7yJiIiMmYlCw3G6sXSKGFew2HPPPfnNb37DX//6V2655RZWr17NfvvtR3d39wY/Z/ny5ZSWlg6+1dXVbXbRIiKSRVq/B/3PoFAxNVjGmAmvK3V2djJnzhyuvfZaTj311FE/JhQKEQqFBt8PBoPU1dXR1dVFSUnJRJ9aRESyQfB+aDzO7SqmntlPQcHCpH7JYDBIaWnpJl+/N6vzsqysjK222ooPP/xwgx+Tm5tLbm7u5jyNiIhko9A70HQisSFY6q2YKjZrQFZPTw8fffQRNTU1yapHREQE7CDUH6EbS6egcQWL888/n6effpo1a9bw/PPPc/TRR+P1evnyl788WfWJiEi2GbyxdA3qq5h6xrUVUl9fz5e//GXWrVvH9OnT2XfffXnxxReZPn36ZNUnIiLZpv1q6HnI7SpkgsYVLO65557JqkNERAR6/w6tF7tdhWwGXUImIiLpIfJxbF6Fbiyd0hQsRETEfbqxNGPoog8REXGXMdDybQi9iULF1KdgISIi7uq6DbrucLsKSRJthYiIiHv6X4LmM9yuQpJIwUJERNwRDcT6KjQAK6MoWIiISOoN3ljaioZgZRb1WIiISOoFLoL+Z9FqReZRsBARkdQK3gcd17hdhUwSbYWIiEjqhN4edmOpZCIFCxERSY3BG0sjaAskc2krREREJp9xoPGrsbHdatbEMRDBIoxFFA8RE/vfkfj7YRP737G32K8PvW8RMZ5hnx97P/H5PgynOSHyXfq9KViIiMjka78Keh9xu4oNihp43i6lG2/sBdrEXrDDeOIv5PEX+MEX9qEX/XA8FISJvbiHh73oR/AQNYnAMPTmbPZWkMFiaNvBiq8AGcDGwz4Da9m1aDOfYoIULEREZHL1Pg6t33e7ig2yDfwgPJfH7WmDj1nxF26LoRft4Qyxl3ZDYlMn1T0jVjxEDL2fLhQsRERk8gzeWJqebAOXhOfwd7t8xONmRJxInxftqUDBQkREJofTD/X/A04v6dis6Ri4LDybR+1pKDwkj4KFiIgknzHQ/G0IvUU63lhqDCwP1/FHuwKFiuRSsBARkeTrvBWCv3G7ilEZAyvCtTxoT3e7lIykORYiIpJc/f+CljPdrmJUxsDPI7P4g13ldikZS8FCRESSJ9qStjeWGgM3Rmby+2i126VkNAULERFJjsEbS9tIxyFYv4rUcGd0httlZDz1WIiISHIELoT+50jH1YrbIzO4LVrjdhlZQcFCREQ2X/Be6LjW7SpG9ZtINbdEZrpdRtbQVoiIiGye0H+h6STS8djm7yNV3BiZ5XYZWUXBQkREJi70XmwIVhreWPqHyHSujdS6XUbW0VaIiIiMT+hd6L4PgndD+B3AS7o1az4QqeTqSJ3bZWQlBQsREdm00H+h+37ouhsi7zEyTKRXqHg4WsEVkdlul5G1FCxERGR9xsTGcXffH1uZiHxAOoeJhD9Hp/GTsEKFmxQsREQkxhgI/WdYmPiIqRAmEh6LlnNJeI7bZWQ9BQsRkWxmDITeiIeJuyCyhqkUJhKeiJbxg/DcePto+p1OySYKFiIi2cYYCL0GwXgDZvQTpmKYSHgqWspF4XkKFWlCwUJEJBsYAwOvDJ3miNYzlcNEwnN2Cd8Nz49fzK5QkQ4ULEREMpUxMPBSPEzcA9EGMiFMJLxgF3NeSKEi3ShYiIhkEuPAwL/i2xz3gN1E7K/6aPwDpnaYSHjJLuLc0AJsLIxCRVpRsBARmeqMA/3Pxxsw/wB2MyPDRHQjnzz1vGYXcXZoC6IKFWlJwUJEZCoyNvT/M7Yy0X0v2AEyOUwk/Nsu5MzQFkQUKtKWgoWIyFRhbOh7NtYz0X0v2G1kQ5hIeMsu4NsKFRtkbEPP6z1E2iK8sv27fOZwG6/Xm/I6LGNMSm+NCQaDlJaW0tXVRUlJSSqfWkRk6jFR6Hsm3oB5HzjrGBkmssM7Tj6nDWxFCA+OQsV6Op7soH5FPZFAZPCx2tparr/+eo455pikPMdYX791u6mISLoxUej9OzR9Az6ogrWLofP2eKiAbAsV7zv5fEuhYoM6nuxg9QWrR4QKgIaGBo499lgefPDBlNajYCEikg5MBHr/Bk1fhw8qYe1B0HUHOB3xD8iuMJHwoZPHNwa2pE+hYlTGNtSvqB/91+IbEueccw62nbrTQOqxEBFxiwlD75OxfonuB8AJkk09E5uy2snlGwNb0YtXoWIDel7vWW+lYjhjDGvXruXZZ5/lgAMOSElNChYiIm4I/gGav6EwsQGfOLmcNrAV3QoVG2R327Q+2Dqmj21qaprkaoYoWIiIpFrfs9B4AgoTo6t3cvj6wFYE8SlUjCIcCBO4O0DbA204vc6YPqempmaSqxqiYCEikkrhD6D+CGBsLwjZpikeKjrxYStUjNC/up/A7wK0/7kdE431T+TOyyW6LoodHL2HwrIsamtr2W+//VJW52Y1b1555ZVYlsU555yTpHJERDKYvQ7WHgJODwoW62tx/Jwa2op2/AoVw/S80cNHyz7inWPfYd3D6zBRQ9EuRSy4bgHb3bsds384e9TPs6zY9/C6665L6TyLCa9YvPzyy/zqV79ip512SmY9IiKZyQlB/ZEQ+YRMua8jmVodP18PbUWbUagAMI6h69kuWu5sofffvbEHLSg9oJTqE6sp2qlo8GPLF5XDCkadY3HdddclbY7FWE0oWPT09HD88cdz22238dOf/jTZNYmIZBZjoPnU2H0epHQm4ZSwzvj4emhLWkxO1ocKJ+zQ/mg7Lb9tIbQmBIDlt5h22DSqv1pN3ty8UT+vfFE5ZQvLBidvXrT9sXz98Ctdmbw5oWBxxhlncNhhh7FkyZJNBotQKEQoFBp8PxgMTuQpRUSmrrbLIPh7t6tISx3Gx2kDW9FkcrM6VNjdNm0r2wjcFSDSGlt18BR6mH7sdKq+XIV/un+TX8PyWhTvVgzAbrO3cSVUwASCxT333MNrr73Gyy+/PKaPX758OZdeeum4CxMRyQjBu2HdJW5XkZY6jZdvDGxJfRaHikhrhMDdAVrvbx084eGf7qfqK1VUHlOJt8idcLA5xhUs1q5dy9lnn83jjz9OXt7oyzGfdvHFF7Ns2bLB94PBIHV1deOrUkRkKur7JzSe6HYVaSlovHxzYEs+NnlZGSoGVg/Q8rsW2v/SjonEtsfy5uVRfWI15UvL8fin7mDscQWLV199lUAgwGc/+9nBx2zb5plnnuGmm24iFAqtt/SSm5tLbm5ucqoVEZkqwh9B/eHo9Mf6uo2H0we2ZJXJz7o5FT3/6aHlzha6nu4abLcp/Ewh1SdVU7pvKZZn6n8/xhUsFi9ezJtvvjnisZNPPpltttmGCy+80LX9HBGRtGK3x4+VdqNgMVKv8XBGaEvez6JQYRxD8LkgzXc20/tG7+DjpQtLqT6pmqKdizby2VPPuIJFcXExO+yww4jHCgsLqaioWO9xEZGsZMJQfxRE1qBjpSP1Gw9nhrbgHacgK0KFE3HoeLSDlt+1MLBqAADLN+yEx7yxtRRMNZq8KSKSLMZA02nQ/0+0UjFSv7E4K7SAt5zCjA8Vds+wEx6BkSc8pn95OjnTc1yucHJtdrB46qmnklCGiEgGWHcFBH/rdhVpJ2Qszg0t4A2nKKNDRaQ1QuCeAG33t2H3xFar/JV+qo6vovLoSrzF2dEuoBULEZFkCP4B2n7gdhVpJ2wszgvN51WnOGNDxcCa+AmPPw+d8Midm0v1idVMWzoNT87UPeExEQoWIiKbq+95aPyq21WknYix+G5oHi86JZgMDBW9b/bSfGczXU8NO+Gxc/yEx36ZccJjIhQsREQ2R3hV/FipGjWHixi4ODSX55zSjAoVxjEE/xmk5bct9LzWM/h46cL4HR6fyawTHhOhYCEiMlF2h46VjiJq4IfhuTzllGVMqHAiDh2PddDy2xYGPhp2wuNz06g6sYr8efkuV5g+FCxERCbChKH+aIisRqsVQ2wDPw7P4e92+ZQJFcY2g5d3+Sv9FO1ShOWN1W73Djvh0TJ0wqPy85VUfbmKnKrMPuExEQoWIiLjZQw0fwv6n0G3lQ5xDFwWns2j9jSYIqGi48mO9a4b91f5qflWDaH6EG33tWF3x4Kjr8IXu8Pj85X4ivXyuSH6zoiIjFf7VdB1h9tVpBVj4IpwHX+0K5hKoWL1BavXezwSiPDJTz4ZfD93TvyEx+ey74THRChYiIiMR/A+aL3Y7SrSijFwdbiWlfZ0t0sZM2Mb6lfUb/RjLL/F3MvnUnZgWdae8JgIRS8RkbHqfxEaT3C7irRiDFwbnsW9dpXbpYxLz+s9I7Y/RmMiBl+pT6FinLRiISIyFuE1sPYwIOp2JWmh23j5S3Qa90UrWW2m3omISNvGQ8V4P06GKFiIiGyK3Qn1h4DTRTYfKzUG/u0U8mC0kr/Z5USxpmTrqnEMPf/u2fQHEhvJLeOjYCEisjEmAg3HQPgjsvVYaZfx8ufoNO6LTucTk4cXgz1FGjQ/LRwI8/GPP6b7X92b/Fh/dezoqYyPgoWIyIYYA82nQ99TZNuxUmPgdaeIB6KV/N0uGxEkpmqo6Hi8g0+u+AQ7aGPlWkxbOo11D60b/YMtqD2/dnCehYydgoWIyIa0r4CuX7tdRUp1GC9/ilZwf7SS+im+OpFgd9usXbGW9j+3A1CwbQFzL5tL3rw8SvYpWX+ORbWf2vNrKV9U7lbJU5qChYjIaIIPQOuFbleREo6BV5wiHoxW8qRdjsPQ+sxUDxXdr3Xz8Y8+JtwUBg/MOHkGNafVYPljv6/yReWULSzb4ORNGT8FCxGRT+t/CRq/QmzQU+ZugbQbH3+MVnBftJImk5sRqxMJTtih6ZdNtPy2BQzkzMph7mVzKdp5/Z4Jy2tRvFuxC1VmJgULEZHhIh9D/eeIHSvNvFDhGHjJKeaBSGX8krDMWZ1I6P+onzU/WEP/+/0AVBxZQe15tXgLvS5Xlh0ULEREEuyu2G2ldieZdqy0Nb46cX90Oi0mBy8GJ0OCRIJxDK33tNJwYwMmbPCV+Zj9g9mUHVjmdmlZRcFCRATix0o/D+EPyZRjpbaBF50SHohU8qxTCjAYJjJldSLh08dIS/YpYc6P5mgOhQsULEREjIHmM6HvSTJh+yPg+HnYruCBSCWtZObqxHAdf+vgk+VDx0hrz62l8thKLCtzf8/pTMFCRKT9Wui61e0qNkvUwPNOCQ9EpvNPpwSLzF2dSLC7bdZevZb2v8SPkW4XP0Y6N8/lyrKbgoWIZLfuh6D1ArermLAmx8/D0UoeiFbSjh8vBjNFR22PR/er8WOkzfFjpKfMoObrQ8dIxT0KFiKSvfpfgcYvuV3FuEUMPGeX8kC0khezZHUiwQk7NN7SSOB3ATCQW5vLnMvmULSTRm+nCwULEclOkU+gfimYqXOstMHJ4eFoBQ9GK+nAjydLVicS+j/sZ80Phx0jPSp+jLRAx0jTiYKFiGQfOwhrDwW7g3Q/ARIx8Ixdxv3RSl5yivEwtDqRyQ2ZwxnHELg7QONNjUPHSH84m7IDytwuTUahYCEi2cVEoeFYCL9POoeKtU4uK6MVrIxWEsSHBwNYGTZdY9PCLfFjpC/Fj5HuGz9GWqFjpOlKwUJEsocx0HIW9P2ddNz+CBuLp+wy7otW8ppTjGfYMdFsWZ0Yrv2xdtYuX4vdbePJ8zBr2Swqj9Ex0nSnYCEi2aPjeui8xe0q1rPGyWVltJKHoxV0D65OZGeYAIh2R1l75Vo6/toBQMH28WOkc3SMdCpQsBCR7ND9MASWuV3FoIix+Ltdxv3R6bzhFI24ACxbAwVA98vdrLlkDZGWCHjjx0hP1THSqUTBQkQy38CraXWstMnJ4fzQPN41hYOrE5l+THRTnLBD4y8aCfzf0DHSuT+dS+GOhW6XJuOkYCEimS2yFtYujd0FkgZ9Fc/bJVwcmks/sSOS2bw6kdD/QfwY6Qc6RpoJFCxEJHPZ3fFjpe24fQLENnBrpIbbozOwAKNAETtGelf8GGnE4CuPHyNdWOZ2abIZFCxEJDOZKDR+AcLv4XaoaDc+LgrN41WnCLJooNXGhJvjx0hf1jHSTKNgISKZqeUc6H0Mt7c/3rALuSA0ny58oFUKANr/2s7aK3WMNFMpWIhI5mm/ATpvdrUEY+CuaBXXRWYB6qUAiAbjx0gf0zHSTKZgISKZpedPEDjH1RK6jYcfh+bylFPmah3ppPulbtb8eNgx0lNnUHOKjpFmIgULEckcA69Dw3GulvC+k8+y0HxaTI6rdaQLJ+TQeHMjgd8HAMity2XuZTpGmskULEQkM0QaYidAXDxW+kh0GleEZ2NjaesD6PugjzU/WMPAhwMAVB5TyaxzZ+kYaYZTsBCRqc/piV2B7tKx0gFjcWW4jj/alRC/LCybGccQ+L8Ajb8Ydoz0R7Mp27/M7dIkBRQsRGRqMzY0fAlCb+NGqFjr5HJeaD6rTaIBMXtChbENPa/3EGmL4K/0U7RLEZFAhDWXrKHn1R4ASvcvZfYPZ+OfpmOk2ULBQkSmtsAy6P0Lbmx/PBkt5UfhuYTxZN3WR8eTHdSvqCcSiAw+5i3xYsIGZ8DBk+eh9rxaKo6u0DHSLKNgISJTV8fN0HFDyp82YuCmyCz+L1pNNm59dDzZweoLVq/3uB2MrRjl1OWwxQ1bkDdbx0izkYKFiExNPX+Blu+k/GkDjp/vhufxlpM41ZBdocLYhvoV9Rv/mLAhd1ZuiiqSdONxuwARkXEb+Dc0HJvyp33JLuZLA9vytlOYtXd99LzeM2L7YzSRlgg9r/ekqCJJN+MKFrfccgs77bQTJSUllJSUsNdee/Hoo49OVm0iIuuL1MePlYZJVV+FY+DXkRl8O7QFQbxZfcV5uCU8po+LtG08fEjmGlewqK2t5corr+TVV1/llVdeYdGiRRx55JH897//naz6RERinH5YdyWs2hrsVlJ1AqTTePlOaAG/iNRgsLJ2pQKg+7VuGn/ZOKaP9VfqFEi2GlePxRFHHDHi/csvv5xbbrmFF198ke233z6phYmIAGAcCN4Fge+C3UwqT3+8ZRdwfmg+7fjJtl6K4SKtERpuaKD9L+2xB2L3vm+Qvzp29FSy04SbN23b5r777qO3t5e99tprgx8XCoUIhUKD7weDwYk+pYhkm96nYvd+hP7NJl/NksgYuC9ayc8idRiy9wIxEzW03ttK4y8bcXodsKDy6EoKdyrk4x9/PPonWVB7fi2WNzu/Z27yxNfUbDx4cW+66biDxZtvvslee+3FwMAARUVFrFy5ku22226DH798+XIuvfTSzSpSRLJM6F0IXAC9f4LBvyBTEyr6jIefhGfzuD0tJc+Xrrpf62btlWsZ+Cg2jrtg+wLqLqyjcPvYaRhPoWe9ORb+aj+159dSvqjclZozn8GHwf7UllwJUWZaIWZ7Qsyywsz1DLB9/mzXqrSMMeP60xoOh/nkk0/o6uri/vvv5/bbb+fpp5/eYLgYbcWirq6Orq4uSkpKNq96Ecks0QC0/Rg6f0VshSK1kzRXOXksC82nweRm7SpFpDVC/fX1dDwau9rcW+pl1pmzqDiqAssz8nsy2uRNrVRsHi8OYI1oEM7FocYKU2cNUOsJM8uKBYiZnhAzrTAFlrP+F5r9FBQsTGptwWCQ0tLSTb5+jztYfNqSJUtYsGABv/rVr5JamIhkEacfOq6Dtp+CCeHGaO5Ho+X8JDwH+1N/qWcLEzEE7g3Q9KumoW2PYyqZ+e2Z+Mo08ihZPBg8GKLDzk54MFRZEWqtAeqsMLPiKw8zrRCzPGHKiDLu4aUuBovN/mlxHGfEioSIyJgZB4K/h8CFKW/MTAgbi5+Fa3nAnk42TtEE6H61m7VXfWrb46I6CrfT1ebjZWHwjrJdUUaEWZ4ws61QbMUhvtowywpTZYXJpIWecQWLiy++mKVLlzJ79my6u7u56667eOqpp3jssccmqz4RyVS9/4g3Zv6HVDZmDtfo5HB+aD7vm/z4Ixn0t/sYhFvDNFzXQMdfh217fGcWFUeuv+0hCbE+BwdrxHZZPjY1Vpg51gCz4tsVMz2x4FBjhcizUv/z7ZZxBYtAIMCJJ55IU1MTpaWl7LTTTjz22GMcdNBBk1WfiGSa0DsQOD9+cVhqGzOHe9Yu4fuheQzEe+mziYkYAn+Ib3v0Ddv2OGMmvlJte3gxED9dkeDDocqKMNsKUeuJrTrMHLZtUWKlfvsuXY3rJ+jXv/71ZNUhIpkuGoDWH0HXbQytDKT+L+OogV9GZnJHdAYWJutCRfcr3ay9Wtsen+bBkIfDlp5+6kZsVcSCQ6UVQYs4Y6NoKiKTy+mD9utg3eXxxsxROthTZJ3xcVFoHq85seFN2RQqwq1hGn7eQMdjw7Y9zppFxf9k97aHN95MeaKvhZP8LRSOdsJCxkXBQkQmh3Eg+H/xxswW3NjuGO41u4jvhuYRxEc29VKYiCFwT4CmW4dte3w+ftoji7c9Yg2WcKi3nTP8jVR7dLdJsmTvT5WITJ7eJ+ONmW/iVmNmgjHwu2gVN0ZmAdk1RbP75fi2x6r4tscOBcy+aDYF2xa4XJl7PPHGy509PZyXU882nn63S8o4ChYikjyht+ONmY/iZmNmQrfx8qPQHJ5xylyrwQ3hQPy0R3zbw1fmY+ZZM6k4Ipu3PWI/h7VWiGU59ezrCY5/NoSMiYKFiGy+aMuwxsxEJ727XfLvOvmcF1pAwGTPLZsmYgjcHaDptmHbHsdWMvP07N72sDAUY3OGv5GjfG34FCgmVfb+pInI5nP6oP3n8cbMMLF/FbobKIyBh+wKrgzXrTdrIJN9etujcMdC6i6sy+ptj0Rj5ld9AU7yN1OkxsyUULAQkfEzDgR/F2/MDOB2Y2ZCv7FYHp7Nn+0Kt0tJmdG2PWadPYtph0/L2m2PxOTLQ+KNmTPUmJlSChYiMj69T0DLORB+C7cbM4f72MnlvNB8PjZ5bpeSEk7EofXu1qFtDw9MP3Y6NafX4CvJzr/aE3NJdoo3Zm6rxkxXZOdPn4iMX+i/EDgPeh9jqI8iPULF36NlXBKeQwRPVmx9BF8KsvaqtYTWxO5pKtwpvu2xTbZue8R+DmdZIc7zN7Cft0uNmS5SsBCRjYs2xxszb2coUKTHXnXEWFwfnsnddnVWTNEMt4Sp/3k9nY93AuAr9zHrrOze9vBgKMLm2/HGTH92fhvSioKFiIzO6YP2a2HdFWnTmDlci+PngvB83nZi/0rP5FDhRBwCdwVovq0Zpz++7XHcdGq+lb3bHonGzBN8Ab6mxsy0kp0/kSKyYcaGrt9B64Vgt5Iu2x3DvWAXc3FoHn14MzpQwCjbHjvHtz22zs5tj0Rj5sHxxswaNWamHQULERnS+/d4Y+Z/SafGzIQPnHzuiUznIbsiXl3mhgpte4yU2Ora0dPL+Tlr1ZiZxhQsRARCb0HLedD3N9KtMTNi4B92GXdHq/iPUxS/0tpKk+qSz4k4BH4foPn2T217nF6Drzgb/8oeasxc5m9gfzVmpr1s/CkVkYRoM7T+ELp+zdAI7vTYq241PlZGK7k3Mp0O/HjiLzB2BqxSGNvQ83oPkbYI/ko/RbsUYXktgv+Kb3t8rG0PiDVmFsYbM49WY+aUoWAhko2cXmi/BtYtBxMh9q/CqNtVYQy84RRyT3Q6T9rlwNClYZlyjLTjyQ7qV9QTCQz1BvgqfeTMzKHvP32x96fFh1wdNg0rC/95HuucGWrMLFZj5pSiYCGSTYwNXXdC68Vgt5EuqxP9xsOjdjl3R6pYZfLxxm+gzDQdT3aw+oLV6z0ebYsSbYuCBdO/GD/tkYXbHonGzIO8HZzpb6TGE3a7JJmA7PvJFclGxkDv3yCwDMJvky6NmWudXO6LVrIyWkkfnsEokQnbHZ9mbEP9ivqNfoxvmo/aZbVY3sz7/W9MojFzB08v5/vr2c7b53ZJshkULEQyXf8LEPgu9D9HOlxl7hh43inh7sh0XnRK8DIUJNyPOpOn5/WeEdsfo4mui9Lzeg/FuxWnqCq3xf4fn2mFWeavZ6EaMzOCgoVIphr4D7R+D3r/zFCgcG/AVZfx8ki0gnuiVTSbnHgzppVGI7cmV6ghNKaPi7Rlx1wGD4YCbL7tb+IYX6saMzOIgoVIpgl/BK0/gO4/kA6B4l0nn3sj0/mLPY3osGOimdhDMRrjGNr/1E79dRvfBknwV/onuSJ3JRozj/cFOFmNmRlJwUIkU0Qaoe0n8Ts9wM2THhFj8UR89sRbTuFgU1626Xmjh/qf1dP3TrxnILbvs0H+6tjR00yU+BlYEm/MnKnGzIylYCEy1dntsO5K6LgBTBQ3VycCjp8HopXcF51OF76Mmj0xHuGmMA03NtDxWAcAnkIPNafV4K/2s+biNaN/kgW152de42aiMXP7eGPm9mrMzHgKFiJTldMD7dfHQoXpw62jo8bAa04R90Sm85RTFistw2ZPjJXT79D822Za7mzBhAxYUHFUBTO/PRP/tNgWh+Wz1ptj4a/2U3t+LeWLyt0qfZIYauKNmQeoMTNrKFiITDVOCDpvhbYfg9OJW4Giz3h4NDqNu6LTWZPBsyfGwhhDx986aLi+gUhLLDAUfbaI2vNqKdhm5NTM8kXllC0sG3XyZqZINGae7m/i82rMzDoKFiJTReLW0bYfQLQRtw5nfuzkcl/8IrD+DJ89MRZ97/SxdsVaev/dC0DOjBxmnTOLsiVlG5yaaXmtjDxS6o3/TH7FF+BUfzPFVrac+ZHhFCxE0p0x0LMSAhdB5APcGG5lG/inXcrd0em85JSMaMbM5NkTGxNpi9B4cyPr/rgODHjyPFSfXE31CdV48jyb/gIZxBNfrVrs7eA7aszMegoWIums9++x4Vah13Hj1tFO4+XhaCX3RKcTGJw9kb2rEwBO2CFwV4Dm/23G6Y1tQ0373DRmnjmTnOocl6tLneH3y+7m6ebb/kZ2VGOmoGAhkp76/wWBC6H/ady4dfQdJ597IlU8Zpdn5eyJ0Rhj6Hqqi4brGgjVx4ZdFWxfQO35tRTtlJlHRNdn4utlFtt7+vict53Fvg4qLfcvsJP0oWAhkk5C/4XA96D3EVI93CpsLP5ul3NXZDrvmOydPTGa/g/7qb+mnu6XuoHYEKuZ35nJtM9Nw/Jk+vfI4CEWKrey+lnqa+cgbwc1nuyYECrjp2Ahkg7Cq6HtRxD8PUNbHqkJFM2Onwei07k/WkkQ3+D6hEIFRDujNP6ykbYH2sABK8ei6oQqZpw8A2+Bd9NfYApL9E3MtUJ8ztfOwd4O6jxjG0su2U3BQsRN0WZouyx2fBQTf5v8QGEMvOIUcU+kiqedUiyGtjmMAgUmYmi9v5WmW5uwg7H/P8oWlzHr7Fnkzsp1ubrJk1ilqrFCfM7bzsG+DrbwDLhdlkwxChYibrA7of1qaP85mAipWp3oNR7+HJ3G3dEqPjF5eONTEbP1ZMdogs8Hqb+2noHVsRfU/C3zqT2/NiOPh8JQmKgkzFJfB4f42tnG6tcwK5kwBQuRVHL6YqO3110BTi+pasg0Bu6IVPPr6AxCg2c7tN0x3MDHA9T/vJ7gs0EAfGU+ar5dQ+VRlRk1vAqGwkQJUQ71tnOIr4OdPL1kfLuIpISChUgqmDB03g5tl4C9jlQeGQ0bi0vDs/mrXZGy55xK7G6bptubaL2nFRM14IWqL1Ux47QZ+Ioz56/IRJgoxOYgbweH+Dr4rKcbn8KEJFnm/KkRSUfGhuBd0Pp9iK5N+dN3GS/nhhbwH6cw5c+d7oxtWPfwOhp/0Ui0I3ZcsmSfEmrPrSVvXp7L1SVHIkzk4rAoHib+n6cbv6XNL5k8ChYik8EY6HkEWi+C8LvgwpbDWieXM0Jb0GRy1JD5Kd2vdVP/s3r63+sHIHdOLrXLaindt9TlyjZf4jSHH4f9vV0c6u1gb28XeQoTkiIKFiLJ1vsPaL0QBl7GjWmZAG/YhZwTWkAf3qweavVpocYQDdc30Pn3TgC8RV5qvlnD9OOmY03hm7JiYSI2+WQvT5BDfe3s7+2i0HLngjrJbgoWIsnS/0osUPQ9iRvTMhMei5bzo/AcHCyFiji7z6blNy20/K4FEzbggcqjK6k5vQZ/ud/t8iZk+FmeXT3dLPV1cKC3kxJd/CUuU7AQ2Vyhd2I9FD0rSfW0zOGMgTui1dwcmUVshUShwjiG9r+203hDI5HW+HXmuxVRe34tBVsWbOKz09HQSO0dPb0s9baz2NdJhUZqSxpRsBCZqMgn0HoJBO8k1dMy1yvFWFweruOPdmX8EYWK3jd7WfuztfS9FbsYK2dW/DrzAzd8nXl6GhqpvbXVz+d87SzxdjBDI7UlTSlYiIxXNADrLoeOX5DKaZkb0m28nBeaz2tOtlyEtXHh1jCNNzbS/ud2ADz5HmacMoOq46vw5E6d68wTTZjzrAE+F7+fo07XkcsUoGAhMlZ2F7T/LPaWwmmZG9Pg5HBmaAvqTW7Wn/xwQg6B38evM++PX2d+ePw68+lT4zrzxPHQmVaIw+IjtedrpLZMMeMKFsuXL+fBBx/k3XffJT8/n7333purrrqKrbfeerLqE3Gf0w8dN8VWKZxu3GjIHM2bdgFnh7agJ4tOfhjb0PN6D5G2CP5KP0W7FIEHOp/spOG6BsKNsX/RF+5YSO35tRTukP7zO4aP1D7MFwsTW2uktkxh4woWTz/9NGeccQa777470WiU733vexx88MG8/fbbFBam/x9gkXEbeBMajoTIGlJ9ZHRjnoiW8f3wXOwsOvnR8WQH9SvqiQSGegt803z4ynwMrIr9q95f5WfWd2ZRvrQ8rfsorPgdLaVEWRofqb2DRmpLhhhXsPjrX/864v3f/OY3VFVV8eqrr7L//vsntTAR13X9Hpq/Ht/2SI9QYQz8LlrF9ZHawRenbNDxZAerL1i93uPR9ijR9ij4YMbJM6g+qRpvfvpeZ57om6izQpzsb2apt0NTMCXJYq2+bnY6bNYzd3V1ATBt2rQNfkwoFCIUCg2+HwwGN+cpRSafCUPLMui82e1KRogYuCpcx0p7OjC+681H20JI94u1TNRg99pEu6KsXb7xcei+Mh81p9Wk7e8psd2xnaeXU30t7Ovt0uqEJIkVf3PAKoLCA6HwEMjfw7WKJhwsHMfhnHPOYZ999mGHHXbY4MctX76cSy+9dKJPI5JakQZoOBoGXnG7khG6jYfvhubzsjP+q7tH20LwV/mpvaCW8kXlySxzkIkY7B576K3XHvn+p96cHmf9x/rH3ssSbYvS83pP2l1tnggUe3m6ONnfwme8vW6XJBnBC9hg5UL+flC4BAoWQ94uYLm/YmcZYya0Dnf66afz6KOP8txzz1FbW7vBjxttxaKuro6uri5KSkom8tQik6P3H9BwLDhdpMOJj4Qmx8+ZoS34xOSNu59iQ1sICfNWzFsvXDhhZ5Mv+pt6M6EkLu/7gDHMf5p7+VymHbrh1dNU8sTXk5Z62znJ36KTHbKZEn8IfJC/JxQeBAWLYv/bSt2Jp2AwSGlp6SZfvye0YnHmmWfypz/9iWeeeWajoQIgNzeX3NzciTyNSGoYA+0rYheGJZYU08Q7Tj5nDmxBN75xhwpjG+pX1G/0Y9Z8fw3NC5pxeofCgwknLxR48jx4i7x4imL/Hevb8I/v/XcvH3zzg00+l7/S/dHcFoYcDMf5WvmKL0C1hljJhCSChAdyPwOFB0PhIsjfBzzpPzF2XMHCGMN3vvMdVq5cyVNPPcW8efMmqy6R1LCD0HQS9DwUfyB9GumeipZycXge0Qme/Oh5vWfE9sdoTNjQ/07/qL/mKfTgLdz0C/8G3wq8SbnYq2iXIvxV/o3+XvzV8aOnroj9zJRic4K/hWN9bbqvQ8Zp2LJczraxHomCRVCwP3in3o274woWZ5xxBnfddRcPP/wwxcXFNDc3A1BaWkp+fv6kFCgyaUL/hfr/gcjHblcygjFwd3Q610QSq4ETe3GOtI3tX8tVX62ibGHZyLBQ4E2bRkjLa1F7Qe2Gt3QsqD2/NuX1Jk541FhhTva1cLhvHbk64SFjMixI+OdBwcFQuBgKDgDfdDcLS4pxBYtbbrkFgAMOOGDE43fccQdf+9rXklWTyOQL3gNNJ6fNBM2EqIFrwrXca1dt1tcxxtDzWs+YPrZ031IX/7U/NuWLymEF6zehVvupPX/ymlBHkwgUW1r9nOJv5kBvJ2mSwSRteYltsRrwzoDCQ2NbGwUHgn/j7QRT0bi3QkSmNBOBwPnQcYPblayn13i4MDSPF53Na2p2Qg6fXP7J4F0ZG+PuFsL4lC8qp2xhmWvHZhMnPHbzdHOKv5ndPD2ajikb4GHwHiHPtFizZeHi2PaGfz6Z/oOju0Ike0QaoeHzMPAvtytZTyB+8mO1ydusoVfhQJhV56+i77994IWKwypY98i60T/YpS2EzWF5rZQfKfViMMASbwdf87ewlWf0nhTJZqPMkig4KLYqkbNdxgeJT1OwkOzQ93QsVNidpFODJsB78ZMfXRM4+TFc75u9rDp/FZG2CN4SL/OumkfJHiWU7FeSFlsIU40HgxfDMd42TvAHmKmbRWWExCyJPMjfd+gIaJrMknCTgoVkNmOg/Vpo/W78gfQ5SgrwrF3ChaH5RDbzzo91f1zHJ5d/gokY8hbkseCaBeTWxY55u72FMLXE1osKsfmKr5Uv+AOU64SHACNnSfy/YbMk9kjpLImpQMFCMpfdDU1fg54H3a5kVPdGKrk6UgeMbzz3cCZqaLi+gcBdAQBKF5Yy97K5eAtH/ovJjS2EqSRx78p0InzN38KRvnXkW+kVQiXVhs+S2AWKDo4Hib2nxCwJNylYSGYKvRM/SrrhqZNusQ1cF57FXXb1Zn2daFeU1Revpvtf3QDMOG0GNd+owdIlFGOWOOExzxrgFH8zS7wdJGH0hkxJicu7iPVFFB4SH0q135ScJeEmBQvJPMF7Y0Ov0uwoKUC/8fC90FyecTbvL6r+Vf2sOncVofoQnjwPcy6dQ/kS9UuMVeKEx2c8PZzsb2EvTzDb+utkULxXIn8/KPt67Cior9LtoqY0BQvJHCYCge9Cx3XEOrTTq0mz1fg4a2ALPjT5THToFUDnM52s+cEanF6HnJoc5l87n4KttDQ7Fl4MDrDQ08nX/C1s7+1zuyRxRTxMeGfEwkTpyZAz3+2iMoaChWSGaDPUfx4GXog/kF6h4gMnjzMHtqAD/4SbNI0xtNzRQuMvGsFA0a5FzLtqHv5y9+/ISHceDB4MR3jXcaI/wGxPaNOfJBnIE3srOgLKTovdwZHlJzgmg4KFTH19z8WuOrc7SLdAAfC8XcwFofmE8Uw4VNj9Nh9f+jGdj3cCUHlcJXXn1SXlLo7MFWuJzcPhS75WvuQPUGmN4ZpUyTDxJkz/VlD+TSj5akaMzU5nChYydRkDHdfHJmliSLejpAAPRitYHp4dn8E3sRAQbgrz0Xkf0f9eP3hh9kWzqTxGe8AbkjjhMY0oJ/pbONrXRpFOeGSZeCOmlQ8lX4ltd+TtmXWDqtyiYCFTk9MDTadA931uVzIqx8CNkZn8Njpjs75Oz+s9rLpgFdGOKL5yH/NXzJ8yI7hTLXHCo84KcbK/maXeDvy6FCzLxHsn8vaEsm9AybHg0Z+XVFOwkKkn9C40HAnhj9yuZFQDxuIHobn8wynbrK/T+kAra69aG2tY3zqf+dfMJ7cmNzlFZpDECY/tPL2c4mthP28XOnGbTeJhwlMBZadC6SmQu7XbRWU1BQuZWoIPQNNXwYRJt6OkAOuMj7MHFvCeKWCiJz9MxLD2Z2tpu78NgLKDyphzyRy8+WoyGy4RKPbyBDnZ38xnvL1ulyQp5Yn9p/DQWCNm0efAUiNzOlCwkKnBRCFwEXRcQzoeJQVY5eRxRmgL1pmJn/yIdERY/d3VsSvPLZj57ZlUn1yNpb1hINY/AbGXlKXedk7ytzDfM+BuUZJCiUbMefGtjhPBP9PtouRTFCwk/UVboOE46H8u/kD6hYqX7GKWheYT2oyTH33v97Fq2SrCTWE8hR7m/nQuZfuXJbfQKSgRJgwW23t6OdTbwRJfB9N1wiNLxP88WTlQ/IVYI2b+fmrETGMKFpLe+p6PHyVdRzoGCoBHotO4LDxns05+dDzRwcc/+hhnwCG3Npf5P59P/vz85BY6hYwME30c4m1nsbeTak9kE58pmSO+OpH7GSj7JpR8SaO1pwgFC0lPxkDHTRA4l3Q9SuoYuCVSw/9Gayb8NYxjaPpVE823NwNQvGcx85bPw1eafX80FSZkqBGzBEq/BqVfh7wd3S5Kxin7/vaS9Of0QtPXofsetyvZoJCxuCQ8h8ftaRP+GnavzZofraHrqS4Aqo6vYtZZs7B82bTEa+IdMxbbefo41NvOIm8nMxQmsowXcKDgwHgj5pHg0QmoqUrBQtJL+P3YraThD9yuZIM6jJdzQgv4r1M44a8RWhvio/M+YuCjASy/xezvz6biiIokVpnOPh0mOljs7VCYyDrxrQ7frFgjZunXwD/b7aIkCRQsJH10r4TGE8CESMetD4CPnVzOCG1BwORMuJ8i+FKQ1Retxu6y8Vf6mf+z+RTuOPGQMjUMCxNWH4f4YmGiRmEiC1mAF4o/H2vELFgElsftoiSJFCzEfSYKrd+H9qtJ16OkAK/ZRZwTWsAAHuwJhApjDK33tFL/83qwoWD7AuZfM5+c6TmTUG06GAoT21p9HKowkcXivRM528caMUuPB+/EtxElvSlYiLuigfhR0mfjD6RnqPhLdBo/Ds/BYWInP5yww9or17Lu4XUATDtsGrO/PxtPbqb9S20oTGwTDxNLvJ3UeMJuFyYpl7ivowhKT4itTuR+VsdEs4CChbin/0WoPwrsNtI1UKx2cvlDZDr32VXEahz/X4qRtgirLlhF7396wQOzzppF1QlVGTT0aihMbG31c6ivnSXeTmYqTGSp+OpE/j6x3oniY8BT4HZRkkIKFpJ6xkDnL6DlHGK9FOnVT9FrPPzdLufBaCVvOYV4BkPP+INA79u9rDpvFZFABG+Rl3nL51Gyd0lyC3aFif97dChMLPZ2MkthIkvFw4S3KrYyUXoy5GzhdlHiEgULSS2nF5q/AcG73K5kBGPg304hD0Ur+Js9jRDWYIyY6CTN9kfb+fiyjzEhQ+7cXBZcu4C8OXnJKzrlhsLEVlY/SxUmsoyXWLgeNvHUyo+N187dKbbdUXgIWHpZyXb6CZDUCX8YP0r6ntuVDGozPv4creDBaAX1Jm/wYiuY+OaMsQ2NNzfScmcLACX7ljDvp/PwFk/FS8RGholEA2atwkSG8hH7yR92wZ+nDPxbQM7WkLMgthLhXxD7394q9UzIehQsJDW6H4HG48H04/bWR8TAP+1SVkYred4piY/ijpnIaY/hot1R1nx/DcF/BgGoPrmamafPxPJOpb98h8LElvEwsURhIoP4iQWHxJ9DC7zVkLPlUHjwDwsQ3kzYupNUUrCQyWVsaP0htC/H7aOkq51cHolW8nC0gi58eDET3uYYzcCaAT5a9hGhj0NYuRZzLpnDtEOmypG6kWHikHiYqFOYmIIsBodPDf5584G/Nh4cthpacfAviG1leKbyFp2kGwULmTzRVmj8IvQ9FX8g9aGi13h4PN6I+V+ncMRWx+auTgzX9c8u1nxvDXaPjb/az4JrFlCwbbp3wg+FiS2sAQ7xtXOQt5M6T8jtwmSTRut3yAP//Hh42GJoxcG/APx1YE3FrTiZihQsZHL0PQ8Nnwe7lVQHCmPgDaeQh0dpxExmmIg9lyHwuwANNzSAgcKdC5m/Yj7+Cn9Snyd5hsLEAmuAQxUm0tiG+h0WQM428X6HBfH+B/U7SPpQsJDkMgbar4XWC0n1raStg42YlTSY3KQ0Ym6MM+Dw8U8/puPRDgAqjq6g7sI6PP70Gno1/NbQ+dYAS+NzJmYrTKSBsfY7xP+ra8NlClCwkOSxO6DxJOj9Y8qeMmLgObuUhyahEXNjwi1hVp2/ir63+8ALdefXUXlcZdoMvUqEKg+GXTw9HOjtZH9vl46Gui7R+wDkbAcFB8R6HtTvIBlEwUKSo/8VaDgaok0pebrVTh4PRyt4OFpBcBIaMSF2bLTn9R4ibRH8lX6KdinC8lr0/LuHVResIrouirfUy/yr51O8W3FSn3siEmEiH5v9vF0c4O1ib2+QYsve9CfLJBkeJLaBwoNjYaJgf/Bmy222km0ULGTzGAOdN0PLuay3H5xkvcbD3+xyVk5yIyZAx5Md1K+oJxIYujDLX+WndGEp6x5ah4kY8rbIY8G1C8idlZvU5x67oX6J6YRZ7OtkobeLXTw9+K30HJGe+YYFCf+W8SBxYCxI+Ka7WplIqihYyMTZQWg+BbofmLSnSDRiPhSt5G92OeFJbMRM6Hiyg9UXrF7v8UggQtt9bQCUHljK3J/MxVuQ2k57Dya+3RMbpb3Y18n+3k62sAbUt+eK4UFiARQeFA8SC8FX7WplIm5RsJCJGfg3NBwFkbWT8uVbjY8/xRsxG1PQiJlgbEP9ivqNfoyn0MO85fNS1qTpxcHGgw+HPTzdHOjrZD9PkOm6ftwFw4PEPChYAoWL4kGixtXKRNKFgoWMjzHQdTs0n0Gskz15Wx8RA8/GGzFfSGEj5nA9r/eM2P4YjdPr0Pvv3kntq0gEqWKiHODtYn9vJ3t5u8m30uvCtsw3LEj4Zg9bkTgA/LPcLEwkbSlYyNg5PdD8zaRfILYq3oj5SLwR0zMJjZhj1f9h/5g+LtKW3NWC4UdCa60BFsdPcezo6WVKTQOf8oYHidphKxIHxIZMicgmKVjI2IT+C/VHQ+SjpHy5nkQjZqSSt83IRsxUhgpjDP3v9tP5dCddT3fR//7YgoW/cvMHYMUCVGx+4o6eXg70xpov52i+RArFr/sG8NYMO7VxAOTMda8skSlMwUI2retOaP4WmAibM/DKGHjdKeKhaAWPp6gRczRO2KHnlZ5YmHima+TWhwWWz8JENtzJ4a+OHT2diESAysVhL08XB/q62MfbRbmOhKbI8CBRFd/aSKxIzNPkSpEkULCQDXP6oOVM6Lpjs75MwPHzJ3saD0YraUphI+Zw0a4oXf/souvpLoIvBHF6hwKSJ99DyV4llC4spXSfUrpf7x71VAgAFtSeXzuu20oTv99yIizydrLQ18Vunm5ydSQ0BbzEwrAB7/T41saBsT4J/wIFCZFJoGAhowu9Fxt4FX5vQp9uDDzvlHBPZDovOCVYDK11pGp1IlQfouuZLjqf6qTnjZ4Rfab+ythMitL9SynevRhP7tAJj/JF5bCC9edYVPupPb829usbMXKEdj+LvZ0s9HWyjdWv17FJNyxIeCqgcHFsRaLwwNhcCf0fIDLpFCxkfcF7oOkUMGEmsvWx1snhqnAdLzil8bkLVkpWJoxj6PtvH53PxPolBj4aGPHreVvkUbawjNKFpRRsW4Dl2fCLTPmicsoWlo06eXM0iVUJL4bPero50NvF/t4uajRCe5INDxLl8SARX5HI2UZBQsQFChYyxBmAwLnQ+csJffqAsfhNZAZ3RKsxKWrEdAYcul/uHuyXiK4bdo20F4p2KYqFif1Lya0d34RMy2tt9EhpIkwUDI7Q7mQvb5BiHQmdRB5IHET2lMZXIxbFg8R2ChIiaWDcweKZZ55hxYoVvPrqqzQ1NbFy5UqOOuqoSShNUiq8Krb1EXprQp/+jF3CleHZBIx/MFRMlmhHlK5nu+h8ppPuF7pxBob1SxR6KNm7hLKFZZTsU4KvJJnZeWiEdrUVZrG3g/29XXzG04Nfr2eTJLFF5YCnOH5iI74qkbsDWOl1k6yITCBY9Pb2svPOO3PKKadwzDHHTEZNkmrdK6HxRDADjHfro8HJ4epwHc85pVjxbY/JMPDxAF1Pd9H5dCe9/+kdUaa/2j+4KlG0W1HSJ2Im5mpsa/WxKH4fx3yN0J4kVvzNAasgNtGyYHGsRyJ3Z7BSO0JdRMZv3MFi6dKlLF26dDJqkVQzYQh8FzquJ/aX+dg7IULG4rfRan4dmTG43ZHMUGFsQ++bvbHmy6c7Ca0ZOdshf+t8SheWUrawjPyt8yflunILQx4OX/EFOM7XphHak8Iitiphg5UH+fsN9UnkfRYs7daKTDWT/qc2FAoRCg29KASDwcl+ShmLyCfQcAwMvBZ/YOyh4p92CcvDdTSbnKSGCaffIfivIF1Pd9H1bBfRjqF+CctnUbRrEWUHlFG6Xyk5NTlJe96RYt+HUmxO9LdwrK+VIvVMJFl8loSVC/l7D21t5O8O1uYPHhMRd016sFi+fDmXXnrpZD+NjEfPn6HxeHB6GU+gaHJyuDpcyzNO2Ua3PYxtxnyaIrIuQtczXXQ900XwX0FMaKgeb5GXkn3j8yX2KsVbPHnL4Infz3QinOJv5n9868jTnIkkSYzJ9kH+XlC4JL4isQd43LpyXkQmy6QHi4svvphly5YNvh8MBqmr08x9V5gotH4f2q9mPFsfYWPxf9Eqbo3UbHLbo+PJjvXnP1T5qb0gNv/BGMPA6li/RNfTXfS+1TuijJyZOZTuH9viKNqlCGuSuyIT/ROzrBBf9zdzqLddjZibbViQyNt9KEjk/z/w5LtdnIhMskkPFrm5ueTm6l8lros0QMMXYOCF+ANjCxUv2MVcEZ5N0xi2PTqe7Bh1YmUkEGH1BatZt986QmtChNaO7Jco2L5gMEzkbZE3Kf0Sn5YIFPOtAb7hb+IAb6cu+5qwRJDwxPoiEjeA5u8NnkK3ixORFFNnVDbofRwavghON2MNFM2On5+Fa/mHUz445GpjjG2oX1G/0Y8JPhvrr7H8FsV7FMe2OPYrJadqsvol1peYPbGDp5fT/M3s5QnqdMe4JYKEFTupMRgk9gXv5F0lLyJTw7iDRU9PDx9++OHg+6tXr+aNN95g2rRpzJ49O6nFyWYyNrRdCut+yuARvk2IGIvfR6v4VaRmXLeN9rzeM/Iyrw2Y8c0ZVB9fjbcwtccGE4Fid0+Qr/ub2cXbm9Lnn9qGBYmc7eI3gB4IBfuBt8zl2kQk3Yw7WLzyyisceOCBg+8n+idOOukkfvOb3yStMNlM0ZbYKkX/MwxOKtyEl+xiLg/X0WByx33aY6B+YNMfBOTNzktpqIgFCjjA28kp/ma28YztWvTslggSQM7Ww4LEQvBOc7UyEUl/4w4WBxxwAMaoWz6t9T4FjceB3cFYAkXA8fOzSC1P2GPb9hgu2hElcHeAlt+3jOnj/ZWpOU7oxWCApd52TvY3M9cT2uTnZK9hQcK/YFiQOAB8090sTESmIPVYZBLjwLoroe0HjGXrI2Lg7mgVv4zMJDrOuz0irRFa/q+FtvvbhkZqx8cTbIi/Onb0dDJ5MHgxHO1t40R/QJeAjWpYkPDNifVIFC6KB4kaNwsTkQygYJEpom3Q+BXoezz+wMZXKl62i7giPJtPTC6MY4Ui1BSi5c4W1j28DhOOPUf+NvnUnFqDcQyrL1z/VAjEnqL2/NoNzrPYXBaGXBy+5GvlK/4AFVZ005+U0RJ/tD/1ffDVQsGSoYu7/LUpr0xEMpuCRSboex4aPg926yY/tNXxc21kFn+zp+HBMNZQMbB2gJY7Wlj3p3WDqxKFOxcy49QZlOxdMnREdAXrz7Go9lN7fmyORXLFNm2KsPmqP8AXfK0UWxtZMskYXmJjsD/VLGsVgm8W5MwD/xzw1YG/Dnyz4/+tBU+eGwWLSBZRsJjKjIH2a6D1ImIrFBve+ogYuDdaxS8iNUTiN0aOZduj/6N+mv+3mY6/dQx++eI9iplx6gyKdi1ab+ZE+aJyyhaWjXny5sTEVkqmEeVkfzNH+9aRnzFjtz3EgkOUEatOVg54a8A/N/5WB/7Zw8JDHXhLXKlYRGQ4BYupyu6I3Uja+6dNfuhrdhGXh+tYY/IY6wpF3zt9NP9vM51Pdg4+VrJvSSxQ7LTxPgnLa1G8W/LnGSTGbtdYYU71N3OYt52cKTV22yL2R85mZAj0gLc6FhT88+L/rftUaKhEAzdEZCpQsJiK+l+OXSAWbdroh7UZH9eFZ/GoXTHmbY+e//TQfHszwX8OXRZXtriMGafMoGCbgs2tfEISUzLnWgOc5m9msbcDX1q+xvqJrTJ8qq/BMy2+wjB/5CpDYpvCN0PXgYtIxlCwmEqMgY6bILCM2AvY6P0EUQP3Radzc2Qm4TFsexhj6Hmlh+ZfN9P9cnfsQQ+UH1LOjJNnkL/AnfsdEoFiG6uP03Ka2M/VKZkbaIa0iuIhYd6nQkP8f/tm6aItEckqChZThd0FTadAz4Mb/bA37EIuD89m1Ri2PYwxBP8ZpPnXzfT+JzaJ0vJZTDt8GtVfqyavzp1Gv8SUzF093Xzd38yunh4XAoXF4JFd/wIo2D8eGj61TeGZ3OOzIiJTjYLFVDDwBjQcDZG1G/yQduPj+vAs/jSGbQ/jGDqf6qT59mb634tNorRyLCqPrqT6q9Xk1KTu7o7hEoFiX08Xp/qb2d7bl+IKEvMdvLFx1UVHQtHhkLNFiusQEZm6FCzSmTHQdRs0n0ms2W/9rQ/bwAPRSm6MzGJgE9seJmroeLyD5v9tZmBVbAS3J99D5bGVVB9fjX96aqZiflpiSuZB3nZO8bewwDO28eDJEQ8TnlIoOiL2VngIeEtTWIOISOZQsEhXTg80fQO6797gh7xpF/DT8Gw+NBtvqnQiDu1/aafljpbBK8s9hR6qvlRF1Zer8JW782Pgic+hOMrbxon+FmpTMiXTE/+vA/6toPjoWJjI/39qoBQRSQIFi3QU+i/UHwWRVaP+cofxcUN4Jo/YlfFtj9E5IYd1D6+j5c4Wws2xF21vqZfqE6qZftx0vMXuvJB6MPgxfMHXyvG+ANM9m74VdfMktjh8sV6JwS2O+ZP8vCIi2UfBIp0YA52/hMC5YKJ8euCVbWBltJIbIrPo38i2h91n0/ZAGy2/ayG6LnaKwVfho/rEaiqPqcRb4E6gsDAU4HCCv4Uv+Fopm9QpmYktjjIo+p/4FsfBGiIlIjLJFCzSRbQFmk6G3kdH/eX/2gVcHp7Ne6YANtCcaXfbBO4NEPh9ALsr9qKdMyOH6pOqqTiyAk+uZ73PmWxWvH+ijCgn+Vv4vK+NwkmZkulh8Hr4nK2h+JhYmMjbQ1scIiIppGCRDnr+BI0ngdO13i91Gi83hWey0q5kKBaMDBWJq8tb/9CK3RMLFLl1uVSfXM20z03D4099oEic8JhuRTjV18wRvnXkJn1K5vAtjgOg+EgoPBxy5ib5eUREZKwULNzk9EHgvNj2BxbD74ZwDDxkV3B9eBZ9eAFrvZtAIm3Dri7vj/1q3oI8Zpw8g/KDyrFSPJ4yESbKibDE28mBvk529XQneUpm/G52z7TYFkfxEVBwEHiTP0JcRETGT8HCLQOvQcMXhzVoDoWKd5x8fhqazbumkNG2PcJNYVp+20LbQ23rXV1eekAplidVgcLgIdbnUWOFONjbwYHeTrb39JG8EoZvcWw3dIojb3ewUr8SIyIiG6dgkWrGhvafQesP+PSNpJ1RD999cRqPBwrJrbQp2MWMuBV01KvLdypkxtc/dXX5JLLiAchgsYU1wBJfLEwssAaSOB0zscXhh4IDY1scRYfHpl6KiEhaU7BIpcgn0Hg89D834mFj4AePFXHN8hChQBvQBoC/yk/tBbXkzc2LXV3+2LCry3ePX12+2/pXlydb7M6O2LrJjp5elng7OcDbyaykzp1IbHFUxIPEEVC4RCOzRUSmGAWLVAneA02ngekf8fBHTh7feLScJ89f/6bSSCDC6gtWj3isZJ/41eU7T+4LbqJfwothd0+Qxb5OFnq7qLCim/7kMRm+xbHDsC2OXbXFISIyhSlYTDa7C1q+DcG7GN6g2Wc83BaZwe9CVfzn6v9u8suULiql5pQaCradvKvLE2EiF4d9vV0s8nayj7eL4qQdD01sceRA4eL4fInDwV+bpK8vIiJuU7CYTH3PQuOXIdocf8BgDPzDLuOqcC3r8BN8vYdIYNOTJ6u+WDUpoSIRJoqIcqC3k0XeLvbwBslL2tHQ+BaHd3p84uURsVDhKUzS1xcRkXSiYDEZTARaL4H2Kxm8ehtY6+RwZbiOF53S+OAoi0jb2MZZj/XjxiIRJiqIDDZf7uLpSdKx0MSqhAW5u8RWJIqOgLzPaotDRCQLKFgkW/h9aPgShN4g0UMQMhZ3Rqv5dWRG/Nqt2KmKvnf6CNwdGNOX9Vduzs2jQ8dCa60BDvJ2cqC3k+08fUk4yRFfkQDwzYGiz8WaLgsOBG/55n5xERGZYhQsksUY6LodWs6KrVjEeymet4u5IjybZpMzGCr63uuj6VdNdD29/qTN0fir/RTtMr5mzeHHQrey+jnI18GB3i7mbfaV5BaDqzCe0tj9G4UHQ8ESTbwUEREFi6SItkHTKdD7x8GHWhw/PwvX8qRTjie+7dH/QT9NtzXR+URn7IM8MO3QaRTsUED91fWjf20Las+vHTHPYkOGHwv9jKeHxfFjoTWbfXvo8LkS+0DhobFVidxdtL0hIiIjKFhsrp6/QtNXwe4AIGLgnmgVt0RmEo2vUPSuGqDp1iY6H++MfY4F5QeXU3NaDXnz8gDwT/dTv6J+RCOnv9pP7fm1lC/a8JZCol/Ch8Oenm4W+zrZ39tF+WYdC00ECWJHQYuWQuFBkL8vePI34+uKiEimU7CYKKcfAt+Fzpsg3sHwul3IT8OzWWPyAIuBNQM03dYUG2wVP2RRtqSMmm/UkL9g5At0+aJyyhaW0fN6D5G2CP7K2PbHaCsVXhxsPORhs7+3iwO9nezjDW7GraHDZkp4Z8SDxMFQsBh80yf4NUVEJBspWEzEwL+h8YsQ/gCAduPhunAdf7Yr8GAYWBui+bZm2h9tH5yUWXpgKTXfrKFgyw0fGbW8FsW7jX6ZVmJlooQoi+IXfO3h6SZnwsdC402XVmHs+GfhwbFVCf+WJHE2t4iIZBkFi/EwDrT/HFovAgy2cVgZreT6yCwG8BBqCNH86+YRd3mU7l9KzTcmNthq8Opxwhzk62CRt5OdPL2Mod1iFIntDS/k7RFflVgSv8xLPwYiIpIcekUZq0gDNJ0AfU8BsRtILwvN4T1TQKgpHigeGQoUJfuUUPPNGgq3H98gqESYmGMNsMTbwSJfJ1tb/RNYRBjWJ+HfcugYaP5CXTEuIiKTRsFiLIL3Q/Op4PTSbbzcHK7hPns60ZYwjf/7CeseWoeJxrYkiv9fMTO/NZPCHcceKBJhYoYV5jDvOg71dTB/3MdCE6czHPBWQuEhQ30S/lnj/FoiIiITo2CxMXY3tHwHgndijMVf7HKuCdeyrtWh6Y562h5sw0TigWL3Ymq+WTPmeRPDeyaWettZ6mtnh3EPrEr0SeTFViKKDon1SeRsrz4JERFxhYLFhvS/EJugGW3gIyePK8KzeTmQS8udTbQ+0IYJxQJF0WeLqPlWDcW7bnp7IREm8rBZ7O1kqa+d3T3d4xilPWxcdt5n4/MkDoL8vcDKmejvVEREJGkULD7NRKHtp7DuMvqMl9si1dzROo2WOwO03Ns6GCgKdy5k5rdmUrR7EdZGVgcSQ6u8wL6eLj7na2dfb9cYL/nSuGwREZlaFCyGC38EjV/G9L/CP+xSfhKo4b3ftRP4w9s4/bFzowU7FDDzWzMp/n/FGwkUieHdsQmYh/naWeTtpMSyN1HAsD4JjcsWEZEpSMEC4vd8/AZazmCtbfHjtjn85bcDBO5+D6cvHii2K6DmmzWU7FOywUARW52I3c1xuK+dg70dTN/oOO1hQcIqhsKFULAICg6A3J01LltERKYcBQt7HTSdRqj7IW5pn86K31q0/P4T7N5YoMjfOp+ab9VQul/pqIEi0Tcx0wpxhG8dh3g7mOMJbeDJhgeJIihYGBtOVXAA5O4ElndSfosiIiKpkt3Bovfv0Hg8j7eHOPO30/nod+uwu2PbFXlb5DHzmzMpPaAUyzMyUCTCRDkRPudr51BvB9uOeqJjtCAxfEVCQUJERDJLdgYLZwBaL2bVJzfxtd+V8Pxve7C7YleY583Po+YbNZQtLhsRKBJhIh+bg70dHOprZ1dPz6emYA4PEoWxIFGwCAoPgNzPKEiIiEjGy75gEXqLzveP5Vt3NPDAHYZoZzsAuXNyqflmDeVLygcv/kqc6PBj2N/bxVJvO/t4g8Pu5xgeJAqGgkTBAZC3i4KEiIhknewJFsahv/EafnDDj7j59gih9tiWR25dLjO+MYNph0zD8lpYiVs+gd083XzO186B3k6KLAew4m8mHiT2jweJAyHvM7pzQ0REsl5GvBLats2zzz5LU1MTNTU17Lfffni9Q6sFAz2rue5nh3D5L1bR0xoLFDmzcqj5eg3TPjcNy5e4+By2tfo4zNfOQb4OKiybWJBwwMqH/P2Hmi3zdlGQEBER+ZQJvTLefPPNrFixgubmZnbeeWduvPFG9thjj2TXNiYPPvggZ599NvX19YOP1dbWcv3113PYYYdx203f4ZIV/0t7SzxQzMhhxtdnUHF4BT4/2FjUWQMc7mvnEG8HdZ4IQ0Fi8VCzZd5nFSREREQ2YdyvlH/4wx9YtmwZv/zlL9lzzz257rrrOOSQQ3jvvfeoqqqajBo36MEHH+TYY4/FmJFTLBsaGvj85z9PWXkOnR1hAPzVfmacMoOqI6dh/F4qiHCYr51DfZ1sZfVieeL3bYwIEv6U/n5ERESmOst8+lV5E/bcc0923313brrpJgAcx6Guro7vfOc7XHTRRZv8/GAwSGlpKV1dXZSUlEysamLbH3Pnzh2xUjEaX6WPmlNnUHFUJcU5hkO8HSz1tfMZbxRPwb7DtjZ2VZAQERHZgLG+fo9rxSIcDvPqq69y8cUXDz7m8XhYsmQJL7zwwqifEwqFCIWGBkYFg8HxPOUGPfvss5sMFQALLpnN0fs5LPXVs1fRzviLvhgPErspSIiIiCTZuIJFW1sbtm1TXV094vHq6mrefffdUT9n+fLlXHrppROvcAOamprG9HEX+T/L17ZapiAhIiKSApN+GcXFF19MV1fX4NvatWuT8nVramrG9HFztz4jfq24QoWIiMhkG9eKRWVlJV6vl5aWlhGPt7S0MGPGjFE/Jzc3l9zc3IlXuAH77bcftbW1NDQ0rNe8CWBZFrW1tey3335Jf24REREZ3bhWLHJycth111154oknBh9zHIcnnniCvfbaK+nFbYzX6+X6668HWO9ysMT711133Yh5FiIiIjK5xr0VsmzZMm677TbuvPNO3nnnHU4//XR6e3s5+eSTJ6O+jTrmmGO4//77mTVr1ojHa2truf/++znmmGNSXpOIiEg2G/cciy9+8Yu0trbyox/9iObmZj7zmc/w17/+db2GzlQ55phjOPLIIzc6eVNERERSY9xzLDZXsuZYiIiISOqM9fV70k+FiIiISPZQsBAREZGkUbAQERGRpFGwEBERkaRRsBAREZGkUbAQERGRpFGwEBERkaRRsBAREZGkUbAQERGRpBn3SO/NlRj0GQwGU/3UIiIiMkGJ1+1NDexOebDo7u4GoK6uLtVPLSIiIpupu7ub0tLSDf56yu8KcRyHxsZGiouL17vufHMEg0Hq6upYu3at7iCZRPo+p46+16mh73Nq6PucGpP5fTbG0N3dzcyZM/F4NtxJkfIVC4/HQ21t7aR9/ZKSEv3QpoC+z6mj73Vq6PucGvo+p8ZkfZ83tlKRoOZNERERSRoFCxEREUmajAkWubm5XHLJJeTm5rpdSkbT9zl19L1ODX2fU0Pf59RIh+9zyps3RUREJHNlzIqFiIiIuE/BQkRERJJGwUJERESSRsFCREREkiZjgsXNN9/M3LlzycvLY8899+Sll15yu6SMsnz5cnbffXeKi4upqqriqKOO4r333nO7rIx35ZVXYlkW55xzjtulZJyGhgZOOOEEKioqyM/PZ8cdd+SVV15xu6yMYts2P/zhD5k3bx75+fksWLCAyy67bJN3TcimPfPMMxxxxBHMnDkTy7J46KGHRvy6MYYf/ehH1NTUkJ+fz5IlS/jggw9SUltGBIs//OEPLFu2jEsuuYTXXnuNnXfemUMOOYRAIOB2aRnj6aef5owzzuDFF1/k8ccfJxKJcPDBB9Pb2+t2aRnr5Zdf5le/+hU77bST26VknI6ODvbZZx/8fj+PPvoob7/9Ntdccw3l5eVul5ZRrrrqKm655RZuuukm3nnnHa666iquvvpqbrzxRrdLm/J6e3vZeeedufnmm0f99auvvpobbriBX/7yl/zrX/+isLCQQw45hIGBgckvzmSAPfbYw5xxxhmD79u2bWbOnGmWL1/uYlWZLRAIGMA8/fTTbpeSkbq7u82WW25pHn/8cbNw4UJz9tlnu11SRrnwwgvNvvvu63YZGe+www4zp5xyyojHjjnmGHP88ce7VFFmAszKlSsH33ccx8yYMcOsWLFi8LHOzk6Tm5tr7r777kmvZ8qvWITDYV599VWWLFky+JjH42HJkiW88MILLlaW2bq6ugCYNm2ay5VkpjPOOIPDDjtsxM+1JM8jjzzCbrvtxnHHHUdVVRW77LILt912m9tlZZy9996bJ554gvfffx+Af//73zz33HMsXbrU5coy2+rVq2lubh7x90dpaSl77rlnSl4XU34JWbK1tbVh2zbV1dUjHq+urubdd991qarM5jgO55xzDvvssw877LCD2+VknHvuuYfXXnuNl19+2e1SMtaqVau45ZZbWLZsGd/73vd4+eWXOeuss8jJyeGkk05yu7yMcdFFFxEMBtlmm23wer3Yts3ll1/O8ccf73ZpGa25uRlg1NfFxK9NpikfLCT1zjjjDN566y2ee+45t0vJOGvXruXss8/m8ccfJy8vz+1yMpbjOOy2225cccUVAOyyyy689dZb/PKXv1SwSKJ7772X3//+99x1111sv/32vPHGG5xzzjnMnDlT3+cMNuW3QiorK/F6vbS0tIx4vKWlhRkzZrhUVeY688wz+dOf/sQ//vEPamtr3S4n47z66qsEAgE++9nP4vP58Pl8PP3009xwww34fD5s23a7xIxQU1PDdtttN+Kxbbfdlk8++cSlijLTBRdcwEUXXcSXvvQldtxxR7761a9y7rnnsnz5crdLy2iJ1z63XhenfLDIyclh11135Yknnhh8zHEcnnjiCfbaay8XK8ssxhjOPPNMVq5cyZNPPsm8efPcLikjLV68mDfffJM33nhj8G233Xbj+OOP54033sDr9bpdYkbYZ5991jsu/f777zNnzhyXKspMfX19eDwjX2a8Xi+O47hUUXaYN28eM2bMGPG6GAwG+de//pWS18WM2ApZtmwZJ510Ervttht77LEH1113Hb29vZx88slul5YxzjjjDO666y4efvhhiouLB/fpSktLyc/Pd7m6zFFcXLxe30phYSEVFRXqZ0mic889l7333psrrriCL3zhC7z00kvceuut3HrrrW6XllGOOOIILr/8cmbPns3222/P66+/zrXXXsspp5zidmlTXk9PDx9++OHg+6tXr+aNN95g2rRpzJ49m3POOYef/vSnbLnllsybN48f/vCHzJw5k6OOOmryi5v0cycpcuONN5rZs2ebnJwcs8cee5gXX3zR7ZIyCjDq2x133OF2aRlPx00nxx//+Eezww47mNzcXLPNNtuYW2+91e2SMk4wGDRnn322mT17tsnLyzPz58833//+900oFHK7tCnvH//4x6h/J5900knGmNiR0x/+8Iemurra5ObmmsWLF5v33nsvJbXp2nQRERFJminfYyEiIiLpQ8FCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJJGwUJERESSRsFCREREkkbBQkRERJLm/wOpyOw0LnkkzgAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(np.linspace(0, 10, 10), np.linspace(0, 10, 10)*0.4*np.random.normal(1, 0.05, 10), color='k', marker='o')\n",
"plt.fill_between(np.linspace(0, 10, 10), \n",
" np.linspace(0, 10, 10)*0.6*np.random.normal(1, 0.05, 10), \n",
" np.linspace(0, 10, 10)*0.2*np.random.normal(1, 0.05, 10), color='gold')\n",
"plt.fill_between(np.linspace(0, 10, 10), \n",
" np.linspace(0, 10, 10)*0.5*np.random.normal(1, 0.05, 10), \n",
" np.linspace(0, 10, 10)*0.3*np.random.normal(1, 0.05, 10), color='limegreen')"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "b531317b-a7a5-4e32-a482-68ad124757d1",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.6"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment