201 lines
76 KiB
Plaintext
201 lines
76 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 35,
|
|
"id": "34599fff-eb66-402d-922b-4cd54ff6d76c",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"%matplotlib widget\n",
|
|
"\n",
|
|
"from collections import defaultdict\n",
|
|
"import numpy as np\n",
|
|
"\n",
|
|
"from matplotlib import pyplot as plt\n",
|
|
"\n",
|
|
"\n",
|
|
"import cristallina as cr\n",
|
|
"\n",
|
|
"from sfdata import SFScanInfo, SFDataFiles\n",
|
|
"\n"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 36,
|
|
"id": "66210aab-28d1-4372-9b07-97286a92716d",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": [
|
|
"scan = SFScanInfo(\"/tmp/scan_info/scan0005_test_bs_scan_info.json\")\n",
|
|
"\n",
|
|
"res = defaultdict(list)\n",
|
|
"res['x'] = scan.readbacks\n",
|
|
"\n",
|
|
"for step in scan:\n",
|
|
" intensity = step[\"SARFE10-PBIG050-EVR0:CALCI\"]\n",
|
|
" res['y'].append(np.mean(intensity.data))\n",
|
|
" "
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 37,
|
|
"id": "fcdc4072-28c0-41cb-b32c-b3107c76b997",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"[<matplotlib.lines.Line2D at 0x7fbb4c6a2210>]"
|
|
]
|
|
},
|
|
"execution_count": 37,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
},
|
|
{
|
|
"data": {
|
|
"application/vnd.jupyter.widget-view+json": {
|
|
"model_id": "6edc7ab4200a47dbaf8a0ba842f922f7",
|
|
"version_major": 2,
|
|
"version_minor": 0
|
|
},
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABl6ElEQVR4nO3deXhU5f3+8fdk3yYTQjJZSICwhbCHRQggomJAsYq2LsWmWhWXAorU1qLtT7sZa1u1/VpxaaVVBGkramo1GgTjAgEFIggh7DshgSSTkJB1zu+PkJGwL8mcSeZ+Xddcl0zOnPkcJJk7z3mez2MxDMNARERERLyGj9kFiIiIiIh7KQCKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl7Gz+wC2jOn08n+/fuxWq1YLBazyxEREZFzYBgGlZWVxMfH4+PjnWNhCoAXYf/+/SQmJppdhoiIiFyAPXv2kJCQYHYZplAAvAhWqxVo+gcUHh5ucjUiIiJyLioqKkhMTHR9jnsjBcCL0HzbNzw8XAFQRESknfHm6VveeeNbRERExIspAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVExONU1tTzu/9tZFNRhdmliHRICoAiIuJx/v75Dl75bAf3z19DfaPT7HJEOhwFQBER8TgfbTgIwI5DVSxYudvkakQ6HgVAERHxKHvLqtl44Ntbv3/+eAsVNfUmViTS8SgAioiIR1mysWn0b2jXCHpEh1JaVcfcT7aZXJVIx6IAKCIiHiWnoCkAXj0gjp9P6gvAq5/vYH/5UTPLEulQFABFRMRjOI7Ws3J7KQBX9Yvhqn4xXNI9ktoGJ3/8qNDk6kQ6DgVAERHxGJ8UFtPgNOhtD6N7VCgWi4VHJ6cA8PbafXyzz2FyhSIdgwKgiIh4jI+Ozf+7ql+M67khiRF8Z3A8hgGZHxRgGIZZ5Yl0GAqAIiLiEWobGsktLAFaBkCAn01MJsDXhy+2HuaTzSVmlCfSoSgAioiIR8jbXsqR2gbs1kAGJ0S0+FpiZAi3j+4GQOb7BTSoObTIRVEAFBERj5CzsQiAK1Ni8PGxnPT1GZf3xhbsz+aDR/jP6r3uLk+kQ1EAFBER0xmGwZKNxQCkn3D7t5ktxJ+ZV/QC4JmczVTVNritPpGORgFQRERMt36fg6KKGkICfEnr2fm0x2WkdSMxMpjiylpe+Wy7GysU6VgUAEVExHQ5x1b/XtYnmiB/39MeF+jnyyPHmkO//Ol2iitq3FKfSEfjcQEwMzOTESNGYLVasdvtTJkyhcLCls0/Fy9ezMSJE4mKisJisZCfn3/SecaPH4/FYmnxuPXWW1scU1ZWRkZGBjabDZvNRkZGBuXl5W14dSIicio5p2j/cjqTB8YxJDGC6rpGnl2yua1LE+mQPC4A5ubmMn36dPLy8sjJyaGhoYH09HSqqqpcx1RVVTFmzBieeuqpM55r2rRpHDhwwPV46aWXWnx96tSp5Ofnk52dTXZ2Nvn5+WRkZLTJdYmIyKntKa1mU1Elvj4WruhrP+vxFouFx441h1705R42H6xs6xJFOhw/sws4UXZ2dos/z5s3D7vdzurVqxk3bhyAK6Tt3LnzjOcKCQkhNjb2lF8rKCggOzubvLw8Ro4cCcArr7xCWloahYWFJCcnX+SViIjIuWhu/jyieyciQgLO6TUjukcysX8MH244yFMfbOLVO0a0ZYkiHY7HjQCeyOFo2vYnMjLyvF/7xhtvEBUVRf/+/Xn44YeprPz2t8QVK1Zgs9lc4Q9g1KhR2Gw2li9ffsrz1dbWUlFR0eIhIiIXp7n9y1X9Tv0L++k8Mqkvfj4Wlm4qZvnWQ21RmkiH5dEB0DAMZs+ezdixYxkwYMB5vfa2225j4cKFfPLJJ/zyl7/krbfe4sYbb3R9vaioCLv95FsNdrudoqKiU54zMzPTNV/QZrORmJh4fhckIiItlFfX8eXOMuD07V9Op0d0GFNHdgXgd+8X4HRqiziRc+Vxt4CPN2PGDNatW8fnn39+3q+dNm2a678HDBhA7969GT58OGvWrGHo0KFA0zySExmGccrnAebMmcPs2bNdf66oqFAIFBG5CMsKi2l0GvSNtZIYGXLer3/wyt4sXrOPDfsrePfrfdyQmtAGVYp0PB47Ajhz5kyysrJYtmwZCQkX/w09dOhQ/P392bJlCwCxsbEcPHjwpONKSkqIiTn1b6GBgYGEh4e3eIiIyIU7n9W/p9I5LJD7x/cE4A/ZhdTUN7ZabSIdmccFQMMwmDFjBosXL2bp0qUkJSW1ynk3bNhAfX09cXFxAKSlpeFwOFi1apXrmJUrV+JwOBg9enSrvKeIiJxebUMjuYUlwIUHQIC7xiYRZwtiv6OGeV/sbKXqRDo2jwuA06dPZ/78+SxYsACr1UpRURFFRUUcPXrUdUxpaSn5+fls3LgRgMLCQvLz811z97Zt28avf/1rvvrqK3bu3Mn777/PTTfdRGpqKmPGjAEgJSWFSZMmMW3aNPLy8sjLy2PatGlce+21WgEsIuIGy7cdpqqukdjwIAZ2sV3weYL8fXk4venn9gvLtlJaVddaJYp0WB4XAOfOnYvD4WD8+PHExcW5HosWLXIdk5WVRWpqKpMnTwbg1ltvJTU1lRdffBGAgIAAPv74YyZOnEhycjIPPPAA6enpLFmyBF/fbzvMv/HGGwwcOJD09HTS09MZNGgQr7/+unsvWETESzXf/p3Qz37audfn6obULvSLC6eytoG/fLylNcoT6dAshmFo2dQFqqiowGaz4XA4NB9QROQ8OJ0GozI/priyln/eeQmX9Ym+6HN+sfUQt/1tJX4+FnJmX0ZSVGgrVCodkT6/PXAEUEREOr51+xwUV9YSFujHqB7n3+f1VMb0imJ8cjQNToPff7CpVc4p0lEpAIqIiNs1N3++LDmaQD/fsxx97uZcnYKPBbI3FPHVztJWO69IR6MAKCIibtc8/+98mz+fTXKslZuHN/Vn/d37BWiWk8ipKQCKiIhb7TpcxeaDR/DzsTA++eQdmS7W7Kv6EOzvy9rd5by//tQ7O4l4OwVAERFxq+bRv5E9IrEF+7f6+e3hQdwzrgcAT3+4iboGZ6u/h0h7pwAoIiJu9VHz7h8prXv793j3jOtBtDWQXYermZ+3q83eR6S9UgAUERG3Ka2qcy3OmNDK8/+OFxrox0MT+gDwl6VbcBytb7P3EmmPFABFRMRtlm4qxmlAv7hwEjqFtOl73Tw8gd72MMqr63lh2dY2fS+R9kYBUERE3Ka5/Utbjv418/P1Yc41fQGY98VO9pRWt/l7irQXCoAiIuIWNfWNfLr5END67V9O5/JkO6N7dqau0ckfPyp0y3uKtAcKgCIi4hZfbD3E0fpG4m1B9I93z/ZbFouFR69JAeDd/P2s21vulvcV8XQKgCIi4hbN7V8m9IvBYrG47X0HdLFxQ2oXAH73PzWHFgEFQBERcQOn02BJQTEAV7np9u/xHp6YTICfDyt3lPLxsTpEvJkCoIiItLm1e8o5dKQWa6AfI5M6u/39u0QEc+eYJAAyPyigoVHNocW7KQCKiEiba779O76vnQA/cz56fnx5TzqF+LOtpIo3v9xjSg0inkIBUERE2lxz+xczbv82Cw/y58ErewPw3JLNHKltMK0WEbMpAIqISJvaXnKEbSVV+PtaGJ8cbWotU0d2o3vnEA4dqePl3G2m1iJiJgVAERFpU823f0f16Ex4kL+ptQT4+fDzq5uaQ7/82XaKHDWm1iNiFgVAERFpU80B0Mzbv8eb2D+W4d06UVPv5JkcNYcW76QAKCIibebQkVpW7y4DYEKKZwRAi8XCnGPNof+9ei8FBypMrkjE/RQARUSkzSwtKMYwYECXcOIjgs0ux2VYt05MHhiHYUDmB5vMLkfE7RQARUSkzXzUfPs3JdbkSk72s0nJ+Pta+HRzCZ9uLjG7HBG3UgAUEZE2cbSukc+3NgUrT5n/d7xunUPJGNUdgCffL6DRqS3ixHsoAIp4oJLKWt7N30e9diuQduyzLSXU1DvpEhFMSpzV7HJOaeYVvbAG+bGpqJLFa/aaXY6I2ygAinig37y3kQffzOcpzU2Sduz41b8Wi8Xkak6tU2gAMy7vBcCfPtrM0bpGkysScQ8FQBEPYxgGK7YfBmDeFzv4Zp/D5IpEzl+j02DppmIA0j3w9u/xbh/dnS4RwRRV1PD3z7ebXY6IWygAiniYvWVHKamsBcBpwJzF67VxvbQ7a3aXcbiqjvAgP0YkRZpdzhkF+fvys0nJALyYu51DR2pNrkik7SkAiniYNcd6pvWICiU8yI/1+xy8tmKXyVWJnJ/m279X9LXj7+v5HzXfGRTPoAQbR2ob+POSLWaXI9LmPP+7UsTLrNnVFADH9Ynm51c3Nav900eF7C8/amZZIudliWv+n+e1fzkVHx8Ljx5rDr1g1W62Fh8xuSKRtqUAKOJhmndNGNatE7eOSGR4t05U1TXyeNYGkysTOTdbi4+w/VAVAb4+XJYcbXY552xUj85MSLHT6DT4fbYWYEnHpgAo4kGq6xooOFAJwNBunfDxsfDkjQPx87GQs/Eg2d8UmVyhyNk13/5N69mZsEA/k6s5Pz+/ui++x77fVh5bjCXSESkAiniQr/c4aHQaxIYHEW8LAqBPjJV7L+sBwBNZG6isqTezRJGzytnY9IuKJzZ/Pptediu3jkgEmppDO9UcWjooBUARD9K8AGRot4gWfdNmXtGbbp1DKKqo4U8fbTarPJGzKqmsZe2ecqB9BkCAWRP6EBrgy9d7Hby3/oDZ5Yi0CQVAEQ/SvABkaNdOLZ4P8vfld1MGAvDPFTv5+tgHrIin+bjgIIYBgxNsxIQHmV3OBYm2BnLfZT0BeDp7E7UNag4tHY8CoIiHMAzjuBHATid9fWzvKG5I7YKh3oDiwY7f/aM9u/vSHsSEB7K37CivLVcbJul4PC4AZmZmMmLECKxWK3a7nSlTplBYWNjimMWLFzNx4kSioqKwWCzk5+ef9nyGYXD11VdjsVh45513WnytrKyMjIwMbDYbNpuNjIwMysvLW/+iRM7BjkNVlFXXE+DnQ//48FMe89jkFCJC/Nl4oIJ5X+x0b4EiZ1Fd18DnWw8B7af9y+kEB/jyk/Sm5tD/t3QL5dV1Jlck0ro8LgDm5uYyffp08vLyyMnJoaGhgfT0dKqqqlzHVFVVMWbMGJ566qmznu+555477R6UU6dOJT8/n+zsbLKzs8nPzycjI6PVrkXkfKzZXQ7AwC42Av18T3lMVFggjx7rDfhMzmb2lFa7qzyRs/p08yFqG5x0jQyhT0yY2eVctO8OTaBvrJWKmgb+b+lWs8sRaVUetz4/Ozu7xZ/nzZuH3W5n9erVjBs3DsAV0nbu3HnGc3399dc888wzfPnll8TFxbX4WkFBAdnZ2eTl5TFy5EgAXnnlFdLS0igsLCQ5ObmVrkjk3Kw5rv/fmdw0PIH/rNnLqh2l/L93v+HVO0ac9pccEXc6/vZvR/g36etjYc41Kdz+6ipeW7GT29O607VziNllibQKjxsBPJHD4QAgMvL89pKsrq7m+9//Ps8//zyxsSffilixYgU2m80V/gBGjRqFzWZj+fLlpzxnbW0tFRUVLR4ireXbBSARZzzOYrHw5A0DCfD1YVlhCe+vV29AMV9Do5Olm5oC4ISU9j3/73iX9Ynm0t5R1Dca/P5DNYeWjsOjA6BhGMyePZuxY8cyYMCA83rtQw89xOjRo7n++utP+fWioiLsdvtJz9vtdoqKTv2BmpmZ6ZovaLPZSExMPK+aRE6nsqaewoPHGkB3PfMIIEAvexj3j29apfjEfzdQod6AYrLVu8ooq64nIsSfEd3P/m+4PXn0mhQsFvjfugOsPTZSL9LeeXQAnDFjBuvWrWPhwoXn9bqsrCyWLl3Kc889d8bjTnWLwjCM0966mDNnDg6Hw/XYs2fPedUlcjr5e8oxDEjoFIz9HFtn3D++Jz2iQimprOVpbVslJmu+/XtFsh0/X4/+aDlvKXHhfG9oAtDUHNow1Bxa2j+P/S6dOXMmWVlZLFu2jISEhPN67dKlS9m2bRsRERH4+fnh59c01fG73/0u48ePByA2NpaDBw+e9NqSkhJiYk59+yIwMJDw8PAWD5HWsGZXOXD2+X/HC/L35Xc3NPUGfGPlblbv0siEmMMwDHIKOkb7l9P5SXoyQf4+fLmzjA83nPzZIdLeeFwANAyDGTNmsHjxYpYuXUpSUtJ5n+PnP/8569atIz8/3/UAePbZZ5k3bx4AaWlpOBwOVq1a5XrdypUrcTgcjB49ulWuReRcrd596gbQZ5PWszPfG5aAYcCji9dTr96AYoItxUfYdbiaAD8fxvWJNrucNhFrC+LusU1bMv4+e5O+16Td87hVwNOnT2fBggW8++67WK1W13w8m81GcHAwAKWlpezevZv9+/cDuPoExsbGtnicqGvXrq5AmZKSwqRJk5g2bRovvfQSAPfccw/XXnutVgCLWzmdhmte0fmMADZ77JoUlm4qpvBgJX/7bIdrbqCIuzTf/h3TszOhgR73sdJq7hvfkze/3M2OQ1UsXLWbH6Z1N7skkQvmcSOAc+fOxeFwMH78eOLi4lyPRYsWuY7JysoiNTWVyZMnA3DrrbeSmprKiy++eF7v9cYbbzBw4EDS09NJT09n0KBBvP766616PSJns7XkCJU1DQT7+9I31nrer+8UGsBj1zT1Bvzzx5vZfVi9AcW9PnK1f2nfzZ/PJizQjwcn9AHguSVbtPhK2jWP+1XtXCbX3nHHHdxxxx0Xfd7IyEjmz59/XucRaW3N7V8GJ9ouePL8jUO78NaavSzfdphfvPsN//yRegOKexysqHHtTT0h5eTOCh3NrSMSmffFDraXVPHiJ9v42aS+ZpckckE8bgRQxNus3nVh8/+OZ7FY+N0NAwnw8+HTzSVkfb2/tcoTOaMlxxZ/DEmMOOcV7O2Zv68Pc47txvP3z3ewv/yoyRWJXBgFQBGTnesOIGeTFBXKzMt7AfCb9zbiqNbtKWl7x+/+4S0mpNi5JCmS2gYnf/yo8OwvEPFACoAiJiqvrmNbSdM+16kXMQLY7N7LetLLHsahI3U8lV1w0ecTOZMjtQ0s33oYgHQvCoAWi8U17/bttfv4Zp/D5IpEzp8CoIiJ1u4uB6BHVCiRoQEXfb4APx+ePNYbcOGqPXy5s/SizylyOp9uLqGu0Un3ziH0soeZXY5bDU6M4LrB8RgGZH6g5tDS/igAipioef5fa4z+NbskKZJbRzRtU/jo4vXUNahfmbSN42//euOio59OTCbA14cvth7mk80lZpcjcl4UAEVM1Frz/07086v7EhUWwJbiI7z86bZWPbcIQH2jk6WbioGO3/7ldBIjQ7hjTHcAMt8voEHNoaUdUQAUMUlDo5P8Y+0zhnaLaNVzR4QE8Mtr+wHwl6Vb2XGoqlXPL/LlzlIcR+uJDA1o9V9g2pPp43thC/Zn88Ej/Gf1XrPLETlnCoAiJik8WEl1XSPWQD9628+/AfTZXDc4nkt7R1HX4OQX76zXHCVpVc23f6/oa8fXx/tu/zazhfjzwJW9AXgmZzPVdQ0mVyRybhQARUzS3AB6SNeINvkAtVgs/HbKAAL9muYovb12X6u/h3gnwzC8sv3L6WSM6kbXyBCKK2t55dMdZpcjck4UAEVMsubYCuCLaQB9Nt06h/LghKbRid/+r4Cyqro2ey/xHpuKKtlbdpRAPx8u7R1ldjmmC/Dz4WeTmvaQf+nTbRRX1phckcjZKQCKmMS1A0gbz5+admkPkmOslFbV8eT76g0oF6959O/S3lGEBHjcjqKmmDwwjiGJEVTXNfJszhazyxE5KwVAEROUVNayu7Qai6VpC6225O/rw5M3NvUG/PfqvazYdrhN3086Pt3+PZnFYuEXk5uaQy/6cjdbDlaaXJHImSkAipiguf1Lb3sYtmD/Nn+/Yd06cdvIrgA89vZ6ahsa2/w9pWM64DjK+n0OLBa4oq8C4PGGd49kYv8YnAZkfrDJ7HJEzkgBUMQEbdX/70x+Nqkv0dZAth+q4oVl6g0oF2bJsdG/oV07EW0NNLkaz/PIpL74+VhYuqmY5VsPmV2OyGkpAIqYYE0b7AByNrZgfx7/TlNvwLmfbGNr8RG3vbd0HB/p9u8Z9YgOc422/+79ApxOtV8Sz6QAKOJmdQ1O1u1t2jze3Q10Jw+M4/LkaOoanTz2tnoDyvmprKknb3vTHFIFwNN74MreWAP92LC/gne/Vvsl8UwKgCJutvFABbUNTiJC/OkRFerW97ZYLPz6+gEE+/uyckcp/9bOBXIecjeXUN9o0CM6lJ7RYWaX47E6hwVy/+U9Afjjh5upqdecW/E8CoAibtZ8+3do105YLO7fQSExMoSHrmrqDfjk+wUcPlLr9hqkfdLq33N355gk4m1B7Cs/yrwvdppdjshJFABF3Gz17uYAGGFaDT8ak0RKXDjl1fX87n/qDShnV9/oZNmmYgDSFQDPKsjfl4cnNjWHfmHZVkrVhF08jAKgiJutdVMD6DPx9/Uh88aBWCyweO0+Pt+i1YpyZqt2lFJR00BUWABDEs37t9ueTBnShf7x4VTWNvCXj9UcWjyLAqCIGx1wHGW/owYfCwxOiDC1liGJEfxwVDcAfvHOes1TkjNqvv17Zd+YNtm7uiPy8bHw2DVNzaHn5+1ix6EqkysS+ZYCoIgbrdlVDkBKXDihgeZvofXwxGRiwgPZebia55duNbsc8VCGYbgC4ATd/j0vo3tFcXlyNA1Og6ez1RxaPIcCoIgbrT5uAYgnsAb586vr+gNNm9hv1vZVcgobD1Swr/woQf4+jO0VZXY57c6ca1LwscAH3xTx1c5Ss8sRARQARdzKjB1AzmZi/1gmpMRQ32jw6OL1alwrJ2ke/bu0dzTBAb4mV9P+9ImxcsuIRKBp5b36b4onUAAUcZOa+kY27G9qAO0pI4DQ1BvwV9f3JyTAl692lbHoqz1mlyQeRu1fLt5DE/oQEuDLmt3lfPBNkdnliCgAirjLN/sc1DcaRIUFkhgZbHY5LXSJCOYn6U0tKzLfL6C4ssbkisRT7Cs/yob9FfhY4Mq+drPLabfs4UFMu7QHAL/P3kRdg9PkisTbKQCKuMm38/8iTGkAfTa3p3VjQJdwKmoa+O176g0oTZYcG/0b1q0TncMCTa6mfbtnXA+irYHsOlzN/LxdZpcjXk4BUMRNPHH+3/H8fH3IvGEQPhbI+no/nxQWm12SeADd/m09oYF+zL6qDwB/WboFx9F6kysSb6YAKOIGhmGw+lgLGDMbQJ/NwAQbd4xOAuCX737D0Tr1BvRmjqP15G0/DMBV/WJNrqZjuGlYAn1iwiivrueFZWq9JOZRABRxg71lRzl0pBZ/XwsDu9jMLueMfpLeh3hbEHtKj/Jn7V7g1T4pLKbBadDLHkZSVKjZ5XQIfr4+zLm6qTn0vOU72VtWbXJF4q0UAEXcoHn+X794G0H+nt1GIzTQj19fPwCAv322nU1FFSZXJGbR7d+2MT45mtE9O1PX4OSPHxaaXY54KQVAETdwzf/zoPYvZzKhXwyT+sfS4DSYo96AXqmuwUluYQmgANjaLBYLj16TgsUC7+TvZ93ecrNLEi+kACjiBq4VwN0izC3kPDxxXX/CAv1Yu7ucN1btNrsccbO87YeprG0g2hrIEJP3re6IBnSxccOQLoCaQ4s5FABF2lhVbQObipq2WPPUFcCnEmsL4qcTm3oDPv3BJg5WqDegN3Ht/Ztix8fH89oWdQQ/mZhMoJ8PedtL+bhAq+7FvRQARdrY13vLaXQaxNmCiLN5VgPos/nBqG4MToygsraBX/93o9nliJsYhsGSAs3/a2tdIoK5c2zTqvvMDwpoaFRzaHEfjwuAmZmZjBgxAqvVit1uZ8qUKRQWtpwku3jxYiZOnEhUVBQWi4X8/PyTznPvvffSs2dPgoODiY6O5vrrr2fTpk0tjikrKyMjIwObzYbNZiMjI4Py8vI2vDrxRmt3lwOe3f7ldHx9LGTeMBBfHwv/W3+ApZsOml2SuME3+yo44KghJMCX0T2jzC6nQ7t/fE8iQwPYVlKlbRjFrTwuAObm5jJ9+nTy8vLIycmhoaGB9PR0qqqqXMdUVVUxZswYnnrqqdOeZ9iwYcybN4+CggI+/PBDDMMgPT2dxsZv+5pNnTqV/Px8srOzyc7OJj8/n4yMjDa9PvE+3+4A0v4CIEC/+HDuPjZK8ct3NlBd12ByRdLWcjY27VU7rne0x69ab+/Cg/x58MreADybs5kjtfr+EvewGB4+87SkpAS73U5ubi7jxo1r8bWdO3eSlJTE2rVrGTJkyBnPs27dOgYPHszWrVvp2bMnBQUF9OvXj7y8PEaOHAlAXl4eaWlpbNq0ieTk5LPWVlFRgc1mw+FwEB4efsHXKB2XYRik/iaH8up63pk+hiGJEWaXdEGq6xq46plP2Vd+lGmXJvHY5H5mlyRtaNJzn7KpqJI/3TSY7w5LMLucDq++0Un6s5+y41AVD1zRi9npZ//8kYujz28PHAE8kcPhACAyMvKCz1FVVcW8efNISkoiMTERgBUrVmCz2VzhD2DUqFHYbDaWL19+yvPU1tZSUVHR4iFyJtsPVVFeXU+gnw/94trvD5mQAD9+e0NTb8BXv9jJN/scJlckbWVPaTWbiirx9bFwRV+72eV4BX9fHx6Z1BT6Xv5sO0UOLbiStufRAdAwDGbPns3YsWMZMGDAeb/+hRdeICwsjLCwMLKzs8nJySEgIACAoqIi7PaTf7jZ7XaKiopOeb7MzEzXfEGbzeYKkyKns+bY7d9BCTYC/Dz62+2sLk+2M3lQHI1Og0ffXk+jegN2SM2rf4d360Sn0ACTq/EeE/vHMrxbJ2rqnTyTo+bQ0vY8+hNpxowZrFu3joULF17Q62+77TbWrl1Lbm4uvXv35uabb6am5tvfrCyWk1sbGIZxyucB5syZg8PhcD327NGEXTmz5gbQ7XX+34kev7Yf1iA/1u118PqKnWaXI21Au3+Yw2Kx8Ojkpi3i/r16r3bgkTbnsQFw5syZZGVlsWzZMhISLmwOis1mo3fv3owbN47//Oc/bNq0ibfffhuA2NhYDh48eUVjSUkJMTGn/sEXGBhIeHh4i4fImazZVQ60zxXAp2IPD+KRSX0B+MOHhRxwHDW5ImlN5dV1rNpZCkB6v1iTq/E+Q7t2YvKgOAwDMt/fdPYXiFwEjwuAhmEwY8YMFi9ezNKlS0lKSmrVc9fW1gKQlpaGw+Fg1apVrq+vXLkSh8PB6NGjW+09xXtV1NSzubipAXRHGQEEmHpJV4Z2jaCqrpEnsjaYXY60omWFxTQ6DZJjrHTtHGJ2OV7pkYl98fe1kLu5hM+2lJhdjnRgHhcAp0+fzvz581mwYAFWq5WioiKKioo4evTbkYbS0lLy8/PZuLGpMW1hYSH5+fmuuXvbt28nMzOT1atXs3v3blasWMHNN99McHAw11xzDQApKSlMmjSJadOmkZeXR15eHtOmTePaa689pxXAImeTv7scw4CukSFEWwPNLqfV+PhYyLxxEH4+Fj7ccJCPNpx6zqy0P7r9a76unUPIGNUdgN/9r0BzbaXNeFwAnDt3Lg6Hg/HjxxMXF+d6LFq0yHVMVlYWqampTJ48GYBbb72V1NRUXnzxRQCCgoL47LPPuOaaa+jVqxc333wzoaGhLF++vMXCjzfeeIOBAweSnp5Oeno6gwYN4vXXX3fvBUuH9W3/vwhzC2kDybFW7hnXA4DHszaod1kHUNvQSG5h04iTAqC5Zl7Ri/AgPzYVVbJ4zV6zy5EOyuP7AHoy9RGSM8n4+0o+23KI31zfn4y07maX0+pq6htJf/ZTdpdW86Mx3Xn8O/3NLkkuwrLCYn4070tiwgNZ8fMrtf+vyV7+dBtPvr8JuzWQd2eMaXfbSHo6fX574AigSEfgdBrkH9sCLrUDzf87XpC/L7+d0tSe6Z/Ld7Jub7m5BclFab79OyElRuHPA/wwrTs9okIprqxl6isr1RtQWp0CoEgb2FJ8hMraBkICfOkbazW7nDYzrk801w+Jx2nAnMXrtZl9O+V0GizR/D+PEuTvy+t3jyQxMpgdh6qY+koeBysUAqX1KACKtIHm+X+DEyLw8+3Y32a/vLYftmB/Nuyv4B/Ld5pdjlyA9fscFFfWEhboR1rPzmaXI8d0iQhm4bRRJHQKZvuhKr7/Sh7FCoHSSjr2J5OISZobQA/rIP3/ziQqLJA5Vzf1BnwmZzP7ytUbsL1pvv17WZ9oAv18Ta5GjpfQKYSF00bRJSKY7SXHQmClQqBcPAVAkTbQvAXc0G4R5hbiJjcPT+SS7pFU1zXy/975Bq0ta1/U/sWzJUY2hcB4WxDbSqqY+spKSiprzS5L2jkFQJFWVlpVx/ZDVQCkJnb8EUBo6g345I0D8Pe18PGmYrK/UW/A9mL34WoKD1bi62NhfHK02eXIaXTtHMLCe0YRZwtia/ERbvtbHoeOKATKhVMAFGlla4/d/u0RHUqn0ACTq3GfXnYr91/WE2jqDVhRU29yRXIuPtrYFNYv6R5JRIj3/Httj7p1DmXhtFHEhgex+eARbntlJYcVAuUCKQCKtDLX/L8O2v7lTH58eS+SjrWu+OOHhWaXI+dAt3/bl+5RoSy8ZxQx4YEUHqzktr+tpLSqzuyypB1SABRpZa4dQLxgAciJgvx9+d2x3oCv5+1yhWHxTGVVdXy5sxRQAGxPkqKaRgLt1kA2FVUy9ZU8yhQC5TwpAIq0ooZGJ1/vcQDesQL4VEb3iuLGoV0wDHh08Xrq1RvQYy3dVIzTgL6xVhIjQ8wuR85Dj+gwFkwbRfSxEHjb31ZSXq0QKOdOAVCkFW0qquRofSPWID96RYeZXY5pfjG5H51C/NlUVMnfP99hdjlyGs23f9M1+tcu9bKHsXDaSKLCAtl4oEIhUM6LAqBIK2q+5ZnatZNXb6cVGRrAo9ekAPDcks3sKa02uSI5UU19I59uKQHgqn6xJlcjF6qX3XosBAawYX8FGX9fhaNaC7Dk7BQARVqRa/5f1whzC/EA3xuWwKgekdTUO/mFegN6nOXbDlFd10icLYgBXcLNLkcuQu8YKwumjaJzaADr9znIeHUljqMKgXJmCoAircibdgA5G4vFwu9uGEiArw+5m0t4b90Bs0uS4zTf/p2QEoPF4r2j1R1Fn2MhMDI0gHV7Hfzw7yvViknOSAFQpJUUV9awp/QoFgsMSYwwuxyP0DM6jOmX9wLgV//dqFtTHsLpNFhSUAxo9W9HkhxrZcG0kXQK8efrvQ5++PdVCoFyWgqAIq1kza5yAJJjrFiD/M0txoPcN74HPaNDOXSklt9/uMnscgTI31tOSWUt1kA/RvXobHY50or6xobzxt2jiAjxJ39PObe/uopKhUA5BQVAkVZy/AIQ+Vagny9P3jAQgAUrd/PVsb5zYp7m27+XJUcT4KePgY6mX3w4b9w9EluwP2t3l3PHvC85UttgdlniYfSdL9JK1uzS/L/TGdmjMzcPTwDg0bfXU9eg3oBm0u4fHV//eJsrBK7eVcYdr65SCJQWFABFWkFdg5N1+5oaQGsF8Kk9ek0KnUMD2HzwCK98tt3scrzWjkNVbC0+gp+PhfHJdrPLkTY0oIuN+XeNJDzIj692lXHnvC+pUgiUYxQARVrBhv0O6hqcdArxJykq1OxyPFJESAC/uLapN+BfPt7CzkNVJlfknXI2FgEwqkdnbMGaq9rRDUywMf/ukViD/Fi1s5Qf/eNLqusUAkUBUKRVfNv/r5NaapzBlCFdGNsritoG9QY0i27/ep9BCRG8ftdIrIF+rNpRyp3/+JKjdY1mlyUmUwAUaQVrd5cDMFTz/87IYrHw2ykDCPTz4fOth3g3f7/ZJXmVw0dqXb+sTFAA9CpDEiN47a5LCAv0I297KXf9UyHQ2ykAirSC40cA5cy6R4XywJW9AfjNexu1d6kbfbypGKcB/ePD6RIRbHY54mapXTvxzzubQuDybYe5+7UvqalXCPRWCoAiF2l/+VGKKmrw9bEwONFmdjntwrRLe9AnJozDVXVkvq/egO6i278yrFsn/nnnCEIDfPli62GmvfaVQqCXUgAUuUjNo38pcVZCAvxMrqZ9CPDzcfUGXPTVHlZuP2xyRR3f0bpGPttSAigAerth3SL5x52XEBLgy2dbDikEeikFQJGL5Nr/V7d/z8vw7pFMHdkVaOoNWNugD6C29PnWQ9TUO+kSEUy/uHCzyxGTjegeyT9+9G0IvPf11QqBXkYBUOQiNTeA1gKQ8/fIxL5EhQWyraSKFz9Rb8C21Nz+5ap+MVqpLgBckhTJq3eMINjfl9zNJdw/f7V+EfMiCoAiF6GmvpEN+ysALQC5ELYQfx7/Tj8A/rpsK9tLjphcUcfU6DT4uKAY0O1faWlUj868escIgvx9WFZYwv3z1ygEegkFQJGLsG6vgwanQbQ1kIROWlV5Ia4dFMdlfaKpa3Ty2NvqDdgW1u4u43BVHeFBflySFGl2OeJh0np25tXbm0Lg0k3FTH9jjbZr9AIKgCIXoXn+39CuEbqtdoGaewMG+fuwYvth3lqzz+ySOpzm1b+X97Xj76sf+3Ky0b2i+PvtIwj082FJQTHTFygEdnT6SSByEZpXAA/T/L+LkhgZwqwJfQD43f82Ulql3oCtSe1f5FyM6RXF324fToCfDzkbDzJjwRrqGxUCOyoFQJELZBgGa3erAXRruWtsEn1jrZRV1/O7/xWYXU6HsbX4CNsPVeHva+GyPtFmlyMe7tLe0bzyw6YQ+NHGg8xcsFYhsINSABS5QLtLqzl0pA5/XwsDuqgB9MXy9/Uh88aBWCzw1pq9LN96yOySOoTm0b+0nlFYg/xNrkbag8v6RPNyxjACfH3I3lDEg28qBHZECoAiF6h5/l//eBtB/r4mV9MxpHbtRMaobgA89s436kvWCo5v/yJyrsYn23npWAh8f30Rsxbl06AQ2KEoAIpcIM3/axsPT0zGbg1kx6EqnsnZrFXBF6Gkspa1e8oBmJBiN7cYaXcu72vnxYyh+Pta+N+6Azz0r68VAjsQjwuAmZmZjBgxAqvVit1uZ8qUKRQWFrY4ZvHixUycOJGoqCgsFgv5+fktvl5aWsrMmTNJTk4mJCSErl278sADD+BwOFocV1ZWRkZGBjabDZvNRkZGBuXl5W18hdJRrNlVDmj+X2sLD/LnV9f1B+DlT7fz6/c20uhUCLwQHxccxDBgYBcbcTa1KZLzd0XfGObeNgx/Xwv//Xo/sxUCOwyPC4C5ublMnz6dvLw8cnJyaGhoID09naqqKtcxVVVVjBkzhqeeeuqU59i/fz/79+/nj3/8I+vXr+cf//gH2dnZ3HXXXS2Omzp1Kvn5+WRnZ5OdnU1+fj4ZGRlten3SMVTVNrCp6FgD6G4R5hbTAV09MI7HrkkBYN4XO5mxYI1uB1+AJQVa/SsXb0K/GP46dSh+Phayvt7Pw//+Wr+UdQAWw8Pvr5SUlGC328nNzWXcuHEtvrZz506SkpJYu3YtQ4YMOeN5/v3vf/ODH/yAqqoq/Pz8KCgooF+/fuTl5TFy5EgA8vLySEtLY9OmTSQnJ5+1toqKCmw2Gw6Hg/Bw7a3pTZZvPcTUv60k3hbE8jlXml1Oh5X19X4e/tfX1DU6Gd6tE6/8cDidQgPMLqtdqK5rIPXXOdQ2OPngwUtJ0f6/cpGyvylixoI1NDgNbkztwh9uGoyvT/vsf6rPbw8cATxR823byMiL617f/D/Zz88PgBUrVmCz2VzhD2DUqFHYbDaWL19+ynPU1tZSUVHR4iHeydUAWvP/2tR1g+N57a5LCA/y46tdZXz3xeXsKa02u6x24bMth6htcJLQKZi+sVazy5EOYNKAWJ6fmoqvj4XFa/fxs/+s00hgO+bRAdAwDGbPns3YsWMZMGDABZ/n8OHD/OY3v+Hee+91PVdUVITdfvKkaLvdTlFR0SnPk5mZ6ZovaLPZSExMvOCapH1rXgCi+X9tb1SPzvzn/tHE24LYXlLFDS8sZ/1ex9lf6OWOb/6sXWqktUwaEMf/fb8pBL61Zi8/f2sdToXAdsmjA+CMGTNYt24dCxcuvOBzVFRUMHnyZPr168fjjz/e4mun+qFoGMZpf1jOmTMHh8PheuzZs+eC65L2y+k0XCsrtQLYPfrEWFn84zH0jbVy6Egtt7y8gk8Ki80uy2M1Og2Wbmr6+9H8P2lt1wyM48+3DsHXx8K/V+9lzuL1CoHtkMcGwJkzZ5KVlcWyZctISEi4oHNUVlYyadIkwsLCePvtt/H3/7YJamxsLAcPHjzpNSUlJcTEnPoHZmBgIOHh4S0e4n22H6qivLqeQD8fzatyo1hbEP++L42xvaKormvkrn9+xb++1C9hp7J6VxmlVXXYgv25pPvFTZ8ROZVrB8Xz7C1D8LHAoq/28OjbCoHtjccFQMMwmDFjBosXL2bp0qUkJSVd0HkqKipIT08nICCArKwsgoKCWnw9LS0Nh8PBqlWrXM+tXLkSh8PB6NGjL+oapGNrnv83OCGCAD+P+xbq0KxB/rx6xwhuTO1Co9PgZ2+t47kl6hV4oubmz1f0tePnq3+j0jauG/xtCHzzyz384t1vFALbET+zCzjR9OnTWbBgAe+++y5Wq9U1H89msxEc3NTHqrS0lN27d7N//34AV5/A2NhYYmNjqaysJD09nerqaubPn99iwUZ0dDS+vr6kpKQwadIkpk2bxksvvQTAPffcw7XXXntOK4DFe605Nv8vVe1fTBHg58Ofbh5MfEQwzy/bynNLtrC//Ci/u2Eg/go7GIbRYv6fSFu6fkgXDANm/yufBSt342OB31w/QPNO2wGP+2k5d+5cHA4H48ePJy4uzvVYtGiR65isrCxSU1OZPHkyALfeeiupqam8+OKLAKxevZqVK1eyfv16evXq1eI8x8/be+ONNxg4cCDp6emkp6czaNAgXn/9dfdesLQ7zSOAw7QAxDQWi4WHJybzuxsG4GOBf321l7v/+RVVtQ1ml2a6rcVH2Hm4mgBfH8b1iTa7HPECU1K78MebBmOxwPy83TyetUGj8u2Ax/cB9GTqI+R9HEfrGfyrjwD46hcTiAoLNLkiWbLxIDMWrqGm3smALuG8escI7Nags7+wg/rrsq384cNCxidH848fXWJ2OeJF/rN6Lz/9z9cYBtwxujuPf6efx44E6vPbA0cARTxZ/rHVv906hyj8eYgJ/WJ48540OocG8M2+Cm58YTnbSo6YXZZpdPtXzPK9YQn8/ruDsFjgH8t38uv3Nmok0IMpAIqcB/X/80xDEiN46/7RdO8cwt6yo3x37nK+2llqdlluV1xR4/olZUKKAqC4383DE3nqxoFA0zaOv3mvQCHQQykAipyHtdoBxGN1jwrlrftHMyQxgvLqeqb+bSUfrD9gdllutaSgqfff4MQIYsK99za4mOuWEV3JPBYCX/1iB0++rxDoiRQARc5Ro9Ng7e5yAIZ2jTC1Fjm1zmGBLJw2igkpMdQ1OPnxgjXM+2KH2WW5TXP7l3Td/hWTff+Srjx5Q1MIfOWzHTz1wSaFQA+jAChyjrYUV3KktoHQAF+SY7S3qqcKDvDlxR8M5QejumIY8Kv/buR3/9vY4fuTVdU28MW2w4Dm/4lnmDqyK7+Z0rSN60ufbuf32YUKgR5EAVDkHDXP/xucGKHmuh7Oz9eH31w/gJ9Naurp+cpnO3jgzbXUNjSaXFnb+XRzCXUNTrp1DqG3PczsckQAyBjVjV9f3x+AF3O38YcPFQI9hT7FRM7Rml3lgPb/bS8sFgs/Ht+LZ28ZjL+vhffWHeCHf1+Fo7re7NLahGv1b0qMx7beEO/0w7TuPPGdfgC88Mk2/vSRdu/xBAqAIueouQG0VgC3LzekJvCPH12CNdCPlTtK+d6Ly9lXftTsslpVQ6OTpYVNC0B0+1c80R1jkvh/1zaFwOeXbeXZJVtMrkgUAEXOQWlVHTsOVQGQqgUg7c6YXlH86740YsOD2FJ8hBv++gUb91eYXVar+XJnGeXV9XQK8dcItXisO8cm8YvJKQD85eMtPLdks8kVeTcFQJFz0Lz/b8/oUCJCAkyuRi5ESlw4i388mj4xYRRX1nLzSyv4fMshs8tqFc23f6/oG6P5qeLR7r60B49d0xQCn1uyhb98rJFAs+gnhcg5cO3/q9GVdi0+Iph/3zeaUT0iOVLbwB3zVrF4zV6zy7oohmGQU9DU/kW3f6U9mDauB3Ou7gvAMzmbeX6pQqAZFABFzoF2AOk4bMH+/PPOS7hucDwNToPZ//qavy7b2m4npRcerGRP6VEC/XwY1yfK7HJEzsm9l/XkkUlNIfCPH23mr8u2mlyR91EAFDmL+kYn6/Y6AI0AdhSBfr48d8sQ7r2sBwB/+LCQx975hoZGp8mVnb+cDU23f8f2iiIkwM/kakTO3f3je/LTiU2tmv7wYSFzP9lmckXeRQFQ5Cw2HajkaH0j4UF+9IxWf7WOwsfHwpyrU/jVdf2xWGDByt3c+/pqqusazC7tvOQUHGv/otu/0g5Nv7wXD6f3AeD32Zt4+VOFQHdRABQ5i+b5f6ldO+Hjo/5qHc3to7sz97ZhBPr58PGmYr7/ykoOHak1u6xzUuSoYd1eBxYLXJmiACjt04wrejP7qqYQ+OT7m/jbZ9tNrsg7KACKnIXm/3V8kwbEsmDaKDqF+PP1nnK+O3e5q+2PJ2se/UtNjCDaGmhyNSIX7oEre/Pglb0B+O3/Cvj7596zh7dZFABFzkIrgL3DsG6deOv+0SRGBrPrcDXfnbvc9f/eU7l2/+gXa3IlIhdv1oTePHBFLwB+895G5n2hENiWFABFzqC4ooa9ZUexWGBwos3scqSN9YgOY/H9YxiUYKO0qo6pr+S5QpanqaypZ8W2pj6Gmv8nHYHFYuGhq/ow4/KmEPir/27kn8t3mltUB6YAKHIGzSNAyTFWrEH+Jlcj7hBtDWThtFFcnhxNTb2Te1//itfzdpld1klyN5dQ32iQFBVKz+hQs8sRaRUWi4WfpPfhx+N7AvB41gZeW7HT3KI6KAVAkTNwzf/T7V+vEhroxys/HM73L0nEacAv3/mG32dvwun0nF6B397+jcFi0eIk6TgsFgs/nZjMfZc1hcC6hvbXnqk9UNMokTNYs7scgGFaAOJ1/Hx9ePKGgcTbgvlTzmbmfrKNA+VHefp7gwnwM/d35/pGJ8s2FQO6/Ssdk8Vi4ZFJyVzR184lSZFml9MhaQRQ5DRqGxpZf6wBtEYAvZPFYmHmlb35w/cG4edj4Z38/dwxbxUVNfWm1rVqRykVNQ10Dg3Q6nTpsCwWi8JfG1IAFDmNDfsrqGt0EhkaQPfOIWaXIya6aXgir94xgtAAX5ZvO8zNL67ggOOoafU03/69oq8dX/WmFJELoAAochprXP3/IjTHShjXJ5pF96YRbQ1kU1ElN76wnMKiSrfXYRhGi/l/IiIXQgFQ5DSaVwDr9q80G9DFxts/Hk0vexgHHDV878XlLD/WisVdCg5Usq/8KEH+PlzaO9qt7y0iHYcCoMgpGIahHUDklBI6hfCf+9K4pHsklTUN3P7qKt7N3+e2928e/RvbK5rgAF+3va+IdCwKgCKnsN9Rw8GKWnx9LAxOiDC7HPEwESEBvHbXJUweGEd9o8GDb+bzYu42DKPt28TkFBQBkK7bvyJyERQARU6hefSvX1y4RlnklIL8ffm/76dy19gkAJ76YBOPZ22gsQ17Be4vP8o3+yqwWOCKFHubvY+IdHwKgCKn0LwARPv/ypn4+Fj45bX9+MXkFCwWeG3FLu6fv5qa+sY2eb8lBU23f4d17URUWGCbvIeIeAcFQJFTaF4Akto1wtxCpF24+9IePP/9oQT4+fDRxoNMfSWP0qq6Vn8frf4VkdaiAChygqN1jWzcXwFoBFDO3eRBccy/ayS2YH/W7C7nu3OXs/twdaudv6KmnrzthwEFQBG5eAqAIidYt7ecBqeB3RpIl4hgs8uRduSSpEjeuj+NLhHB7DhUxY1zv+DrPeWtcu5PCkuobzToGR1Kj+iwVjmniHgvBUCRE7j2/+3WSQ2g5bz1slt5+8ej6R8fzqEjddz6cp5r396L8e3t39iLPpeIiAKgyAnU/08ulj08iEX3pjGuTzRH6xu5+7WvWLhq9wWfr67BySfHQqRu/4pIa1AAFDmOYRis1Q4g0grCAv34++3DuWlYAo1OgzmL1/PMR4UX1Ctw5Y7DVNY2EBUWSGpiROsXKyJex+MCYGZmJiNGjMBqtWK325kyZQqFhYUtjlm8eDETJ04kKioKi8VCfn7+Sed5+eWXGT9+POHh4VgsFsrLy086pqysjIyMDGw2GzabjYyMjFMeJ95j1+FqDlfVEeDrw4Au4WaXI+2cv68PT39vEA9e2RuAvyzdysP/Xkd9o/O8ztN8+3dCih0fH01LEJGL53EBMDc3l+nTp5OXl0dOTg4NDQ2kp6dTVVXlOqaqqooxY8bw1FNPnfY81dXVTJo0iUcfffS0x0ydOpX8/Hyys7PJzs4mPz+fjIyMVr0eaV+a278M6BJOoJ8aQMvFs1gsPHRVH566cSC+PhbeWrOXO//xJUdqG87p9YZhsETtX0SklfmZXcCJsrOzW/x53rx52O12Vq9ezbhx4wBcIW3nzp2nPc+sWbMA+OSTT0759YKCArKzs8nLy2PkyJEAvPLKK6SlpVFYWEhycvLFXYi0S5r/J23l1ku6EhMexI/fWMNnWw5x84sr+MePRmAPDzrj6zbsr2C/o4Zgf1/G9IpyU7Ui0tF53AjgiRwOBwCRkZGtet4VK1Zgs9lc4Q9g1KhR2Gw2li9ffsrX1NbWUlFR0eIhHcvxK4BFWtvlfe0suncUUWEBbDxQwQ0vLGdrceUZX/PRsdG/cX2iCPLXqLSItA6PDoCGYTB79mzGjh3LgAEDWvXcRUVF2O0n76Vpt9spKio65WsyMzNd8wVtNhuJiYmtWpOY60htA4VFTaFeC0CkrQxKiGDx/WPoERXKvvKj3PjCclbtKD3t8Wr/IiJtwaMD4IwZM1i3bh0LFy5sk/OfqsebYRin7f02Z84cHA6H67Fnz542qUvM8fWecpwGdIkIJuYst+VELkbXziG8df9ohnXrREVNAz/420r+t+7AScftKa2m4EAFPha4ou/Jv7CKiFwojw2AM2fOJCsri2XLlpGQkNDq54+NjeXgwYMnPV9SUkJMzKknWgcGBhIeHt7iIR2Ha/6fRv/EDTqFBvDG3SOZ2D+GukYnMxau4W+fbW9xzJKCpp9Rw7tHEhkaYEaZItJBeVwANAyDGTNmsHjxYpYuXUpSUlKbvE9aWhoOh4NVq1a5nlu5ciUOh4PRo0e3yXuKZ2teATysa4S5hYjXCPL35YXbhnF7WjcMA377vwJ+/d+NOJ1NvQKbb/+ma/WviLQyj1sFPH36dBYsWMC7776L1Wp1zcez2WwEBzfty1paWsru3bvZv38/gKtPYGxsLLGxTfNkioqKKCoqYuvWrQCsX78eq9VK165diYyMJCUlhUmTJjFt2jReeuklAO655x6uvfZarQD2Qk6nwRqNAIoJfH0sPHFdf7p0CubJ9zfx6hc7KKo4yhPf6c/KY3MD1f5FRFqbx40Azp07F4fDwfjx44mLi3M9Fi1a5DomKyuL1NRUJk+eDMCtt95KamoqL774ouuYF198kdTUVKZNmwbAuHHjSE1NJSsry3XMG2+8wcCBA0lPTyc9PZ1Bgwbx+uuvu+lKxZNsP3SEipoGgvx9SInTrX1xL4vFwj3jevKX76cS4OvD++uLmPx/n9PoNOgTE0a3zqFmlygiHYzFuJB9iQSAiooKbDYbDoejVecDllXV8cpn27lxaBd62a2tdl45vUVf7uaRt9ZzSVIk/7o3zexyxIut2HaYe17/isqapkbR0y/vyU8n9jW5KpGOpa0+v9sTjxsBFPjlu9/wwifb+PPHW80uxWus2VUOqP+fmC+tZ2feun80cbYgfCxw7aB4s0sSkQ5IAdAD3T++JwDvrdtPYdGZm8RK61i9WzuAiOfoE2MlZ/ZlLJl9maYkiEibUAD0QP3jbVw9IBbDgD9/vNnscjo8R3U9W4uPADBUK4DFQ4QF+tEjOszsMkSkg1IA9FCzJvTBYoH31xexcb+2nGtLa/Y0jf517xxC57BAk6sRERFpewqAHio51srkgXEAPLdEo4Btaa3av4iIiJdRAPRgsyb0xmJp2gx+/V6H2eV0WJr/JyIi3kYB0IP1slu5fnDTCkCNAraNRqdB/u5yQCuARUTEeygAergHruyNjwU+3lRM/p5ys8vpcAqLKqmqayQs0I8+Meq5KCIi3kEB0MP1iA7jhtQEAJ7N0Shga2ve/3dIYgS+PhaTqxEREXEPBcB24IEre+HrYyF3cwmrd5WaXU6H4tr/V+1fRETEiygAtgPdOofyvaHNo4BbTK6mY2keAdQKYBER8SYKgO3EjCt64e9r4fOth1i5/bDZ5XQIh47UsvNwNQCpiQqAIiLiPRQA24nEyBBuHp4IwDM5mzEMw+SK2r+1x1b/9raHYQvxN7cYERERN1IAbEemX96LAF8fVu4oZcU2jQJerNW71P9PRES8kwJgOxIfEcz3L9EoYGtpnv+n/n8iIuJtFADbmR9f3osAPx++2lXGZ1sOmV1Ou1Xf6GTd3nIAhnaLMLUWERERd1MAbGdiwoP4wchugEYBL0bBgQpq6p3Ygv3pERVmdjkiIiJupQDYDt03vgdB/j7k7ynnk8ISs8tpl5rn/6V2jcBHDaBFRMTLKAC2Q3ZrED9M6w5oFPBCrWne/1cLQERExAspALZT947rQUiAL+v3OVhSUGx2Oe2OawcQLQAREREvpADYTnUOC+T20d2BplFAp1OjgOeqyFHDvvKj+FhgcGKE2eWIiIi4nQJgO3bPpT0IDfCl4EAFH20sMrucdqO5/UtybDhhgX4mVyMiIuJ+CoDtWKfQAO4cmwQ07RGsUcBz03z7d5jav4iIiJdSAGzn7h7bA2ugH4UHK3n/mwNml9MurN6tHUBERMS7KQC2c7YQf+66tGkU8LklW2jUKOAZ1dQ3smFfBaAdQERExHspAHYAd45Nwhbsz9biI7y3br/Z5Xi0Dfsd1DU66RwaQNfIELPLERERMYUCYAcQHuTPPeN6APDnJVtoaHSaXJHnWrOrHGhq/2KxqAG0iIh4JwXADuL20d3pFOLP9kNVvJOvUcDTad4BRPP/RETEmykAdhBhgX7ce1lPAP7y8RbqNQp4EsMwXAtANP9PRES8mQJgB/LDtG50Dg1gd2k1i9fsNbscj7O37CgllbX4+VgYlGAzuxwRERHTKAB2ICEBftw/vnkUcCt1DRoFPF5zA+j+8eEE+fuaXI2IiIh5FAA7mNtGdiPaGsi+8qP8e/Ues8vxKM0NoFM1/09ERLycAmAHExzgy4+PjQI+v3QrtQ2NJlfkOdbsLgc0/09EREQBsAP6/iVdiQkP5ICjhkVfahQQoLqugY0HmhpAD1UAFBERL+dxATAzM5MRI0ZgtVqx2+1MmTKFwsLCFscsXryYiRMnEhUVhcViIT8//6Tz1NbWMnPmTKKioggNDeW6665j796WCyPKysrIyMjAZrNhs9nIyMigvLy8Da/OPYL8fZlxeS8A/rpsKzX1GgVct9dBo9MgNjyIeFuQ2eWIiIiYyuMCYG5uLtOnTycvL4+cnBwaGhpIT0+nqqrKdUxVVRVjxozhqaeeOu15Zs2axdtvv82bb77J559/zpEjR7j22mtpbPw2DE2dOpX8/Hyys7PJzs4mPz+fjIyMNr0+d7l5RCLxtiAOVtSyYOVus8sxXfMCkKHdItQAWkREvJ7FMAyP3jy2pKQEu91Obm4u48aNa/G1nTt3kpSUxNq1axkyZIjreYfDQXR0NK+//jq33HILAPv37ycxMZH333+fiRMnUlBQQL9+/cjLy2PkyJEA5OXlkZaWxqZNm0hOTj5rbRUVFdhsNhwOB+Hh4a130a1kwcrdPPr2eqLCAvnsZ5cTHOC9K1/v/ueXLCko5heTU7j70h5mlyMiIiby9M9vd/C4EcATORwOACIjI8/5NatXr6a+vp709HTXc/Hx8QwYMIDly5cDsGLFCmw2myv8AYwaNQqbzeY6pr373rAEEjoFc+hILfPzdpldjmkMw3AtANH8PxEREQ8PgIZhMHv2bMaOHcuAAQPO+XVFRUUEBATQqVPLD/uYmBiKiopcx9jt9pNea7fbXcecqLa2loqKihYPTxbg58MDV/YG4MXcbVTVNphckTl2Hq6mtKqOAD8f+sd75296IiIix/PoADhjxgzWrVvHwoULW+V8hmG0mP91qrlgJx5zvMzMTNeCEZvNRmJiYqvU1ZZuTO1Ct84hHK6q47UV3jkK2Nz/b2AXG4F+3nsbXEREpJnHBsCZM2eSlZXFsmXLSEhIOK/XxsbGUldXR1lZWYvni4uLiYmJcR1z8ODBk15bUlLiOuZEc+bMweFwuB579nh+ixU/Xx8ePDYK+NKn26isqTe5IvfT/r8iIiIteVwANAyDGTNmsHjxYpYuXUpSUtJ5n2PYsGH4+/uTk5Pjeu7AgQN88803jB49GoC0tDQcDgerVq1yHbNy5UocDofrmBMFBgYSHh7e4tEeXDc4nh7RoZRX1/OPL3aaXY7bNY8ADu0aYW4hIiIiHsLjAuD06dOZP38+CxYswGq1UlRURFFREUePHnUdU1paSn5+Phs3bgSgsLCQ/Px819w9m83GXXfdxU9+8hM+/vhj1q5dyw9+8AMGDhzIhAkTAEhJSWHSpElMmzaNvLw88vLymDZtGtdee+05rQBuT44fBXzls+04jnrPKGBlTT2FBysBGKot4ERERAAPDIBz587F4XAwfvx44uLiXI9Fixa5jsnKyiI1NZXJkycDcOutt5KamsqLL77oOubZZ59lypQp3HzzzYwZM4aQkBD++9//4uv77RywN954g4EDB5Kenk56ejqDBg3i9ddfd9/FutG1g+LpbQ+joqaBVz/fYXY5bvP1HgeGAQmdgrGHqwG0iIgItIM+gJ6svfUR+t+6A0xfsAZroB+fPXI5ESEBZpfU5v68ZAvPLtnM9UPi+fOtqWaXIyIiHqC9fX63BY8bAZS2c/WAWPrGWqmsbeBvn3nHKKBrBxDd/hUREXFRAPQiPj4WZk3oA8C8L3ZQWlVnckVty+k0XAFQK4BFRES+pQDoZSb2j6F/fDhVdY28/Ol2s8tpU9tKjlBZ00Cwvy99Y61mlyMiIuIxFAC9jMVi4aFjo4D/XL6TQ0dqTa6o7aw+1v5lcKINP1/9UxcREWmmT0UvdGWKnUEJNo7WN/JS7jazy2kzmv8nIiJyagqAXshisfDQVU2jgK+t2EVxRY3JFbWN5hFAzf8TERFpSQHQS43vE01q1whqG5zM7YCjgOXVdWwrqQIgVSOAIiIiLSgAeimLxcJPrmra8eSNlbspcnSsUcC1u8sB6BEVSmRox+93KCIicj4UAL3YmF6duaR7JHUNTl74ZKvZ5bSq5vl/Gv0TERE5mQKgFzt+LuCbq/awr/zoWV7Rfmj+n4iIyOkpAHq5tJ6dSevRmbpGJ88v7RijgA2NTr7eUw7A0G4RptYiIiLiiRQAxTUK+O+v9rCntNrkai5e4cFKquoasQb60duuBtAiIiInUgAULkmK5NLeUTQ4Df5v6Razy7loa44tABnSNQJfH4u5xYiIiHggBUABcO0R/Naafew8VGVyNRdnzS41gBYRETkTBUABmhZLjE+OptFp8Jd2Pgro2gFEC0BEREROSQFQXJr3CH5n7T62lRwxuZoLc+hILbsOV2OxwJDECLPLERER8UgKgOIyODGCCSl2nAb85eP2OQrYfPu3tz0MW7C/ydWIiIh4JgVAaaF5LmDW1/vZcrDS5GrO3+rd6v8nIiJyNgqA0sKALjYm9o/BMOC5djgKuHZXOaAdQERERM5EAVBO0jwK+L91B9hUVGFyNeeursHJ13vLAY0AioiInIkCoJwkJS6cyYPiAHgup/2MAhYcqKC2wUlEiD89okLNLkdERMRjKQDKKc26sjcWC2RvKOKbfQ6zyzknq4/r/2exqAG0iIjI6SgAyin1jrFy3eB4AJ5b0j5GAV39/7pGmFuIiIiIh1MAlNN64Mre+FhgScFB1h2bW+fJXDuAaP6fiIjIGSkAymn1jA5jSmoXAJ7J2WxyNWd2wHGU/Y4afCwwOCHC7HJEREQ8mgKgnNEDV/TG18fCJ4Ulrjl2nmjNsfYvKXHhhAb6mVuMiIiIh1MAlDPqHhXKd4c2jQI+t8RzRwG/nf+n278iIiJnowAoZzXzit74+Vj4bMshVu0oNbucU2oenVT/PxERkbNTAJSzSowM4abhiQA864FzAWvqG9mwv6lVjUYARUREzk4BUM7JjCt64e9rYcX2wyzfdsjsclr4Zp+D+kaDqLBAEiODzS5HRETE4ykAyjnpEhHMrSO6Ak27gxiGYXJF3zq+/58aQIuIiJydAqCcsx9f3pMAPx9W7Szli62HzS7HRfP/REREzo8CoJyzOFswUy9pGgV8JqfQI0YBDcNgze5yQA2gRUREzpUCoJyXH1/ekyB/H9bsLid3c4nZ5bC37CgllbX4+1oY2MVmdjkiIiLtggKgnBe7NYiMUd2AphXBZo8CNs//6xdvI8jf19RaRERE2guPC4CZmZmMGDECq9WK3W5nypQpFBYWtjjGMAyeeOIJ4uPjCQ4OZvz48WzYsKHFMdu2beOGG24gOjqa8PBwbr75Zg4ePNjimLKyMjIyMrDZbNhsNjIyMigvL2/rS2z37r2sJ8H+vny918HSTcWm1uKa/6f2LyIiIufM4wJgbm4u06dPJy8vj5ycHBoaGkhPT6eqqsp1zNNPP80zzzzD888/z5dffklsbCxXXXUVlZWVAFRVVZGeno7FYmHp0qV88cUX1NXV8Z3vfAen0+k6z9SpU8nPzyc7O5vs7Gzy8/PJyMhw+zW3N1Fhgdw+ujvQtEewmaOArhXA3SJMq0FERKS9sRhm38M7i5KSEux2O7m5uYwbNw7DMIiPj2fWrFk88sgjANTW1hITE8Pvf/977r33Xj766COuvvpqysrKCA8PB5pG+yIjI8nJyWHChAkUFBTQr18/8vLyGDlyJAB5eXmkpaWxadMmkpOTz1pbRUUFNpsNh8Pheh9vUVpVx6W/X0pVXSMvZQxjYv9Yt9dQXdfAwCc+otFpsPznVxAfoR6AIiJydt78+d3M40YAT+RwNO3wEBkZCcCOHTsoKioiPT3ddUxgYCCXXXYZy5cvB5oCocViITAw0HVMUFAQPj4+fP755wCsWLECm83mCn8Ao0aNwmazuc4jpxcZGsCPxiQBTXMBnU73/x7x9R4HjU6DOFuQwp+IiMh58OgAaBgGs2fPZuzYsQwYMACAoqIiAGJiYlocGxMT4/raqFGjCA0N5ZFHHqG6upqqqip++tOf4nQ6OXDggOs8drv9pPe02+2u85yotraWioqKFg9vdvelSVgD/dhUVMkH35z676wtfdsAWvP/REREzodHB8AZM2awbt06Fi5ceNLXTtzxwTAM13PR0dH8+9//5r///S9hYWGuYd6hQ4fi6+t72nOceJ4TZWZmuhaM2Gw2EhMTL+by2r2IkADuHNs0Cvjcks00unkUcM2u5vl/CoAiIiLnw2MD4MyZM8nKymLZsmUkJCS4no+NbZprduIoXXFxcYtRwfT0dLZt20ZxcTGHDh3i9ddfZ9++fSQlJbnOc+KqYGiac3ji6GKzOXPm4HA4XI89e/Zc9HW2d3eOTSI8yI8txUd4b91+t71vUwPob7eAExERkXPncQHQMAxmzJjB4sWLWbp0qSuwNUtKSiI2NpacnBzXc3V1deTm5jJ69OiTzhcVFUVERARLly6luLiY6667DoC0tDQcDgerVq1yHbty5UocDscpzwNNcw3Dw8NbPLydLdifaZf2AODPS7bQ0Og8yytax45DVZRV1xPg50P/eDWAFhEROR8eFwCnT5/O/PnzWbBgAVarlaKiIoqKijh69CjQdNt21qxZPPnkk7z99tt888033HHHHYSEhDB16lTXeebNm0deXh7btm1j/vz53HTTTTz00EOu1b0pKSlMmjSJadOmkZeXR15eHtOmTePaa689pxXA8q07xnQnIsSf7YeqyPraPaOAzf3/BnWxEeDncf+MRUREPJqf2QWcaO7cuQCMHz++xfPz5s3jjjvuAOBnP/sZR48e5cc//jFlZWWMHDmSjz76CKvV6jq+sLCQOXPmUFpaSvfu3Xnsscd46KGHWpzzjTfe4IEHHnCtKL7uuut4/vnn2+7iOihrkD/3jOvB09mF/PnjLVw3OB4/37YNZc37/w7T/D8REZHz5vF9AD2Z+gh9q6q2gUufXkZpVR1Pf28QNw9v2wUyE5/9lMKDlbz4g2FMGuD+HoQiItJ+6fPbA28BS/sUGujHfZc1zQX8y8dbqG/DuYAVNfVsLm7a9UU7gIiIiJw/BUBpNRmjuhMVFsjesqP8Z/XeNnuf/N3lGAYkRgZjtwa12fuIiIh0VAqA0mqCA3z58fieADy/dCu1DY1t8j7N7V+GqQG0iIjIBVEAlFY1dWRXYsID2Vd+lH991TajgKvVAFpEROSiKABKqwry92X65b0A+OvSrdTUt+4ooNNpkH9sBbC2gBMREbkwCoDS6m4ZkUicLYiiihreXLW7Vc+9pfgIlbUNhAT40jfWevYXiIiIyEkUAKXVBfr5MuOKY6OAn2xr1VHA5vl/gxMi2rzXoIiISEelT1BpEzcNS6RLRDAllbXMz9vVauf9dv5fRKudU0RExNsoAEqbCPDz4YErm0YB536yjeq6hlY5r2sFsBaAiIiIXDAFQGkzNw5NoGtkCIer6nhtxcWPApZV1bG9pAqA1EQFQBERkQulAChtxt/Xhweu7A3AS7nbOFJ7caOAa/c0jf71iA6lU2jARdcnIiLirRQApU1NGRJPUlQoZdX1/HP5zos6l2v+n9q/iIiIXBQFQGlTfr4+PHhsFPDlT7dTUVN/wedas6sc0Pw/ERGRi6UAKG3uO4Pj6WUPw3G0nnmf77ygczQ0OsnfUw5oBFBERORiKQBKm/P1sTBrQtMo4N8+346j+vxHATcVVXK0vhFroB+97WGtXaKIiIhXUQAUt7hmQBzJMVYqaxr4++fbz/v1ze1fhnSNwMfH0trliYiIeBUFQHELHx8LD13VNAr46hc7KauqO6/Xr9ml/n8iIiKtRQFQ3Ca9Xyz94sI5UtvAK5+d3yjg6t1aASwiItJaFADFbZpGAfsA8I/lOzl8pPacXldcWcOe0qNYLE23gEVEROTiKACKW01IsTOwi43qukZe/vTcRgGb27/0sVsJD/Jvw+pERES8gwKguJXFYmH2sVHAf67YSUnl2UcB1zbf/tX8PxERkVahAChuNz45miGJEdTUO5n7ybazHv/tDiARbVyZiIiId1AAFLc7fhRw/spdHKyoOe2xdQ1O1u1zAFoBLCIi0loUAMUUl/aOYni3TtQ1OHlh2dbTHrdhv4O6BiedQvxJigp1Y4UiIiIdlwKgmOL4UcCFq/awv/zoKY9bs7scaGr/YrGoAbSIiEhrUAAU06T17MzIpEjqGp389TSjgM0NoLUAREREpPUoAIppjh8F/NdXe9hTWn3SMWvUAFpERKTVKQCKqUb26MzYXlHUNxonjQLuLz/KAUcNvj4WBifaTKpQRESk41EAFNM17xH879V72XW4yvV88+hfSpyVkAA/U2oTERHpiBQAxXTDukVyWZ9oGp0G/7f021HAb/v/6faviIhIa1IAFI/QvEfw4jV72XGoaRSweQWw+v+JiIi0LgVA8QhDEiO4sq8dpwF/+XgLNfWNbDjWAFojgCIiIq1LAVA8RvMo4Lv5+3h77T4anAbR1kASOgWbXJmIiEjHogAoHmNAFxvp/WJwGvDr/24Emvb/VQNoERGR1qUAKB5l1oSmUcCj9Y2A5v+JiIi0BY8LgJmZmYwYMQKr1YrdbmfKlCkUFha2OMYwDJ544gni4+MJDg5m/PjxbNiwocUxRUVFZGRkEBsbS2hoKEOHDuU///lPi2PKysrIyMjAZrNhs9nIyMigvLy8rS9RzqBffDjXDIx1/Vnz/0RERFqfxwXA3Nxcpk+fTl5eHjk5OTQ0NJCenk5V1bf94Z5++mmeeeYZnn/+eb788ktiY2O56qqrqKysdB2TkZFBYWEhWVlZrF+/nhtvvJFbbrmFtWvXuo6ZOnUq+fn5ZGdnk52dTX5+PhkZGW69XjnZrAl98LGANdCPAV3UAFpERKS1WQzDMMwu4kxKSkqw2+3k5uYybtw4DMMgPj6eWbNm8cgjjwBQW1tLTEwMv//977n33nsBCAsLY+7cuS0CXefOnXn66ae56667KCgooF+/fuTl5TFy5EgA8vLySEtLY9OmTSQnJ5+1toqKCmw2Gw6Hg/Dw8Da4eu+1elcZgX4+CoAiItLq9PntgSOAJ3I4mlqBREZGArBjxw6KiopIT093HRMYGMhll13G8uXLXc+NHTuWRYsWUVpaitPp5M0336S2tpbx48cDsGLFCmw2myv8AYwaNQqbzdbiPGKOYd06KfyJiIi0EY/eX8swDGbPns3YsWMZMGAA0DS3DyAmJqbFsTExMezatcv150WLFnHLLbfQuXNn/Pz8CAkJ4e2336Znz56u89jt9pPe0263u97jRLW1tdTW1rr+XFFRcXEXKCIiImICjx4BnDFjBuvWrWPhwoUnfe3E1iCGYbR47he/+AVlZWUsWbKEr776itmzZ3PTTTexfv36057jVOc5XmZmpmvBiM1mIzEx8UIvTURERMQ0HhsAZ86cSVZWFsuWLSMhIcH1fGxs0wrRE0fpiouLXaOC27Zt4/nnn+fVV1/lyiuvZPDgwTz++OMMHz6cv/71r67zHDx48KT3LSkpOWl0sdmcOXNwOByux549e1rlWkVERETcyeMCoGEYzJgxg8WLF7N06VKSkpJafD0pKYnY2FhycnJcz9XV1ZGbm8vo0aMBqK6uBsDHp+Xl+fr64nQ6AUhLS8PhcLBq1SrX11euXInD4XCd50SBgYGEh4e3eIiIiIi0Nx43B3D69OksWLCAd999F6vV6hrps9lsBAcHY7FYmDVrFk8++SS9e/emd+/ePPnkk4SEhDB16lQA+vbtS69evbj33nv54x//SOfOnXnnnXfIycnhvffeAyAlJYVJkyYxbdo0XnrpJQDuuecerr322nNaASwiIiLSXnlcG5jTzb+bN28ed9xxB9A0SvirX/2Kl156ibKyMkaOHMlf//pX10IRgC1btvDzn/+czz//nCNHjtCrVy8efvjhFm1hSktLeeCBB8jKygLguuuu4/nnnyciIuKcatUychERkfZHn98eGADbE/0DEhERaX/0+e2BcwBFREREpG0pAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl/G4RtDtSXMHnYqKCpMrERERkXPV/LntzZ3wFAAvQmVlJQCJiYkmVyIiIiLnq7KyEpvNZnYZplAj6IvgdDrZv38/Vqv1tDuYXKiKigoSExPZs2ePVzap1PV79/WD/g68/fpBfwe6/ra7fsMwqKysJD4+Hh8f75wNpxHAi+Dj40NCQkKbvkd4eLhXfuM30/V79/WD/g68/fpBfwe6/ra5fm8d+WvmnbFXRERExIspAIqIiIh4GQVADxUYGMjjjz9OYGCg2aWYQtfv3dcP+jvw9usH/R3o+r37+tuaFoGIiIiIeBmNAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqAHeuGFF0hKSiIoKIhhw4bx2WefmV2S23z66ad85zvfIT4+HovFwjvvvGN2SW6VmZnJiBEjsFqt2O12pkyZQmFhodlluc3cuXMZNGiQq/FrWloaH3zwgdllmSYzMxOLxcKsWbPMLsVtnnjiCSwWS4tHbGys2WW53b59+/jBD35A586dCQkJYciQIaxevdrsstyie/fuJ/0bsFgsTJ8+3ezSOhQFQA+zaNEiZs2axWOPPcbatWu59NJLufrqq9m9e7fZpblFVVUVgwcP5vnnnze7FFPk5uYyffp08vLyyMnJoaGhgfT0dKqqqswuzS0SEhJ46qmn+Oqrr/jqq6+44ooruP7669mwYYPZpbndl19+ycsvv8ygQYPMLsXt+vfvz4EDB1yP9evXm12SW5WVlTFmzBj8/f354IMP2LhxI3/605+IiIgwuzS3+PLLL1v8/8/JyQHgpptuMrmyjkVtYDzMyJEjGTp0KHPnznU9l5KSwpQpU8jMzDSxMvezWCy8/fbbTJkyxexSTFNSUoLdbic3N5dx48aZXY4pIiMj+cMf/sBdd91ldiluc+TIEYYOHcoLL7zAb3/7W4YMGcJzzz1ndllu8cQTT/DOO++Qn59vdimm+fnPf84XX3zhVXd/zmTWrFm89957bNmyBYvFYnY5HYZGAD1IXV0dq1evJj09vcXz6enpLF++3KSqxEwOhwNoCkHeprGxkTfffJOqqirS0tLMLsetpk+fzuTJk5kwYYLZpZhiy5YtxMfHk5SUxK233sr27dvNLsmtsrKyGD58ODfddBN2u53U1FReeeUVs8syRV1dHfPnz+fOO+9U+GtlCoAe5NChQzQ2NhITE9Pi+ZiYGIqKikyqSsxiGAazZ89m7NixDBgwwOxy3Gb9+vWEhYURGBjIfffdx9tvv02/fv3MLstt3nzzTdasWeN1I/7NRo4cyWuvvcaHH37IK6+8QlFREaNHj+bw4cNml+Y227dvZ+7cufTu3ZsPP/yQ++67jwceeIDXXnvN7NLc7p133qG8vJw77rjD7FI6HD+zC5CTnfhbjmEY+s3HC82YMYN169bx+eefm12KWyUnJ5Ofn095eTlvvfUWt99+O7m5uV4RAvfs2cODDz7IRx99RFBQkNnlmOLqq692/ffAgQNJS0ujZ8+e/POf/2T27NkmVuY+TqeT4cOH8+STTwKQmprKhg0bmDt3Lj/84Q9Nrs69/v73v3P11VcTHx9vdikdjkYAPUhUVBS+vr4njfYVFxefNCooHdvMmTPJyspi2bJlJCQkmF2OWwUEBNCrVy+GDx9OZmYmgwcP5s9//rPZZbnF6tWrKS4uZtiwYfj5+eHn50dubi5/+ctf8PPzo7Gx0ewS3S40NJSBAweyZcsWs0txm7i4uJN+4UlJSfGaxYDNdu3axZIlS7j77rvNLqVDUgD0IAEBAQwbNsy14qlZTk4Oo0ePNqkqcSfDMJgxYwaLFy9m6dKlJCUlmV2S6QzDoLa21uwy3OLKK69k/fr15Ofnux7Dhw/ntttuIz8/H19fX7NLdLva2loKCgqIi4szuxS3GTNmzEntnzZv3ky3bt1Mqsgc8+bNw263M3nyZLNL6ZB0C9jDzJ49m4yMDIYPH05aWhovv/wyu3fv5r777jO7NLc4cuQIW7dudf15x44d5OfnExkZSdeuXU2szD2mT5/OggULePfdd7Fara7RYJvNRnBwsMnVtb1HH32Uq6++msTERCorK3nzzTf55JNPyM7ONrs0t7BarSfN9wwNDaVz585eMw/04Ycf5jvf+Q5du3aluLiY3/72t1RUVHD77bebXZrbPPTQQ4wePZonn3ySm2++mVWrVvHyyy/z8ssvm12a2zidTubNm8ftt9+On5+iSpswxOP89a9/Nbp162YEBAQYQ4cONXJzc80uyW2WLVtmACc9br/9drNLc4tTXTtgzJs3z+zS3OLOO+90/duPjo42rrzySuOjjz4yuyxTXXbZZcaDDz5odhluc8sttxhxcXGGv7+/ER8fb9x4443Ghg0bzC7L7f773/8aAwYMMAIDA42+ffsaL7/8stkludWHH35oAEZhYaHZpXRY6gMoIiIi4mU0B1BERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TL/H+lHdNEfpEu4AAAAAElFTkSuQmCC",
|
|
"text/html": [
|
|
"\n",
|
|
" <div style=\"display: inline-block;\">\n",
|
|
" <div class=\"jupyter-widgets widget-label\" style=\"text-align: center;\">\n",
|
|
" Figure\n",
|
|
" </div>\n",
|
|
" <img src='data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy81sbWrAAAACXBIWXMAAA9hAAAPYQGoP6dpAABl6ElEQVR4nO3deXhU5f3+8fdk3yYTQjJZSICwhbCHRQggomJAsYq2LsWmWhWXAorU1qLtT7sZa1u1/VpxaaVVBGkramo1GgTjAgEFIggh7DshgSSTkJB1zu+PkJGwL8mcSeZ+Xddcl0zOnPkcJJk7z3mez2MxDMNARERERLyGj9kFiIiIiIh7KQCKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl7Gz+wC2jOn08n+/fuxWq1YLBazyxEREZFzYBgGlZWVxMfH4+PjnWNhCoAXYf/+/SQmJppdhoiIiFyAPXv2kJCQYHYZplAAvAhWqxVo+gcUHh5ucjUiIiJyLioqKkhMTHR9jnsjBcCL0HzbNzw8XAFQRESknfHm6VveeeNbRERExIspAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVExONU1tTzu/9tZFNRhdmliHRICoAiIuJx/v75Dl75bAf3z19DfaPT7HJEOhwFQBER8TgfbTgIwI5DVSxYudvkakQ6HgVAERHxKHvLqtl44Ntbv3/+eAsVNfUmViTS8SgAioiIR1mysWn0b2jXCHpEh1JaVcfcT7aZXJVIx6IAKCIiHiWnoCkAXj0gjp9P6gvAq5/vYH/5UTPLEulQFABFRMRjOI7Ws3J7KQBX9Yvhqn4xXNI9ktoGJ3/8qNDk6kQ6DgVAERHxGJ8UFtPgNOhtD6N7VCgWi4VHJ6cA8PbafXyzz2FyhSIdgwKgiIh4jI+Ozf+7ql+M67khiRF8Z3A8hgGZHxRgGIZZ5Yl0GAqAIiLiEWobGsktLAFaBkCAn01MJsDXhy+2HuaTzSVmlCfSoSgAioiIR8jbXsqR2gbs1kAGJ0S0+FpiZAi3j+4GQOb7BTSoObTIRVEAFBERj5CzsQiAK1Ni8PGxnPT1GZf3xhbsz+aDR/jP6r3uLk+kQ1EAFBER0xmGwZKNxQCkn3D7t5ktxJ+ZV/QC4JmczVTVNritPpGORgFQRERMt36fg6KKGkICfEnr2fm0x2WkdSMxMpjiylpe+Wy7GysU6VgUAEVExHQ5x1b/XtYnmiB/39MeF+jnyyPHmkO//Ol2iitq3FKfSEfjcQEwMzOTESNGYLVasdvtTJkyhcLCls0/Fy9ezMSJE4mKisJisZCfn3/SecaPH4/FYmnxuPXWW1scU1ZWRkZGBjabDZvNRkZGBuXl5W14dSIicio5p2j/cjqTB8YxJDGC6rpGnl2yua1LE+mQPC4A5ubmMn36dPLy8sjJyaGhoYH09HSqqqpcx1RVVTFmzBieeuqpM55r2rRpHDhwwPV46aWXWnx96tSp5Ofnk52dTXZ2Nvn5+WRkZLTJdYmIyKntKa1mU1Elvj4WruhrP+vxFouFx441h1705R42H6xs6xJFOhw/sws4UXZ2dos/z5s3D7vdzurVqxk3bhyAK6Tt3LnzjOcKCQkhNjb2lF8rKCggOzubvLw8Ro4cCcArr7xCWloahYWFJCcnX+SViIjIuWhu/jyieyciQgLO6TUjukcysX8MH244yFMfbOLVO0a0ZYkiHY7HjQCeyOFo2vYnMjLyvF/7xhtvEBUVRf/+/Xn44YeprPz2t8QVK1Zgs9lc4Q9g1KhR2Gw2li9ffsrz1dbWUlFR0eIhIiIXp7n9y1X9Tv0L++k8Mqkvfj4Wlm4qZvnWQ21RmkiH5dEB0DAMZs+ezdixYxkwYMB5vfa2225j4cKFfPLJJ/zyl7/krbfe4sYbb3R9vaioCLv95FsNdrudoqKiU54zMzPTNV/QZrORmJh4fhckIiItlFfX8eXOMuD07V9Op0d0GFNHdgXgd+8X4HRqiziRc+Vxt4CPN2PGDNatW8fnn39+3q+dNm2a678HDBhA7969GT58OGvWrGHo0KFA0zySExmGccrnAebMmcPs2bNdf66oqFAIFBG5CMsKi2l0GvSNtZIYGXLer3/wyt4sXrOPDfsrePfrfdyQmtAGVYp0PB47Ajhz5kyysrJYtmwZCQkX/w09dOhQ/P392bJlCwCxsbEcPHjwpONKSkqIiTn1b6GBgYGEh4e3eIiIyIU7n9W/p9I5LJD7x/cE4A/ZhdTUN7ZabSIdmccFQMMwmDFjBosXL2bp0qUkJSW1ynk3bNhAfX09cXFxAKSlpeFwOFi1apXrmJUrV+JwOBg9enSrvKeIiJxebUMjuYUlwIUHQIC7xiYRZwtiv6OGeV/sbKXqRDo2jwuA06dPZ/78+SxYsACr1UpRURFFRUUcPXrUdUxpaSn5+fls3LgRgMLCQvLz811z97Zt28avf/1rvvrqK3bu3Mn777/PTTfdRGpqKmPGjAEgJSWFSZMmMW3aNPLy8sjLy2PatGlce+21WgEsIuIGy7cdpqqukdjwIAZ2sV3weYL8fXk4venn9gvLtlJaVddaJYp0WB4XAOfOnYvD4WD8+PHExcW5HosWLXIdk5WVRWpqKpMnTwbg1ltvJTU1lRdffBGAgIAAPv74YyZOnEhycjIPPPAA6enpLFmyBF/fbzvMv/HGGwwcOJD09HTS09MZNGgQr7/+unsvWETESzXf/p3Qz37audfn6obULvSLC6eytoG/fLylNcoT6dAshmFo2dQFqqiowGaz4XA4NB9QROQ8OJ0GozI/priyln/eeQmX9Ym+6HN+sfUQt/1tJX4+FnJmX0ZSVGgrVCodkT6/PXAEUEREOr51+xwUV9YSFujHqB7n3+f1VMb0imJ8cjQNToPff7CpVc4p0lEpAIqIiNs1N3++LDmaQD/fsxx97uZcnYKPBbI3FPHVztJWO69IR6MAKCIibtc8/+98mz+fTXKslZuHN/Vn/d37BWiWk8ipKQCKiIhb7TpcxeaDR/DzsTA++eQdmS7W7Kv6EOzvy9rd5by//tQ7O4l4OwVAERFxq+bRv5E9IrEF+7f6+e3hQdwzrgcAT3+4iboGZ6u/h0h7pwAoIiJu9VHz7h8prXv793j3jOtBtDWQXYermZ+3q83eR6S9UgAUERG3Ka2qcy3OmNDK8/+OFxrox0MT+gDwl6VbcBytb7P3EmmPFABFRMRtlm4qxmlAv7hwEjqFtOl73Tw8gd72MMqr63lh2dY2fS+R9kYBUERE3Ka5/Utbjv418/P1Yc41fQGY98VO9pRWt/l7irQXCoAiIuIWNfWNfLr5END67V9O5/JkO6N7dqau0ckfPyp0y3uKtAcKgCIi4hZfbD3E0fpG4m1B9I93z/ZbFouFR69JAeDd/P2s21vulvcV8XQKgCIi4hbN7V8m9IvBYrG47X0HdLFxQ2oXAH73PzWHFgEFQBERcQOn02BJQTEAV7np9u/xHp6YTICfDyt3lPLxsTpEvJkCoIiItLm1e8o5dKQWa6AfI5M6u/39u0QEc+eYJAAyPyigoVHNocW7KQCKiEiba779O76vnQA/cz56fnx5TzqF+LOtpIo3v9xjSg0inkIBUERE2lxz+xczbv82Cw/y58ErewPw3JLNHKltMK0WEbMpAIqISJvaXnKEbSVV+PtaGJ8cbWotU0d2o3vnEA4dqePl3G2m1iJiJgVAERFpU823f0f16Ex4kL+ptQT4+fDzq5uaQ7/82XaKHDWm1iNiFgVAERFpU80B0Mzbv8eb2D+W4d06UVPv5JkcNYcW76QAKCIibebQkVpW7y4DYEKKZwRAi8XCnGPNof+9ei8FBypMrkjE/RQARUSkzSwtKMYwYECXcOIjgs0ux2VYt05MHhiHYUDmB5vMLkfE7RQARUSkzXzUfPs3JdbkSk72s0nJ+Pta+HRzCZ9uLjG7HBG3UgAUEZE2cbSukc+3NgUrT5n/d7xunUPJGNUdgCffL6DRqS3ixHsoAIp4oJLKWt7N30e9diuQduyzLSXU1DvpEhFMSpzV7HJOaeYVvbAG+bGpqJLFa/aaXY6I2ygAinig37y3kQffzOcpzU2Sduz41b8Wi8Xkak6tU2gAMy7vBcCfPtrM0bpGkysScQ8FQBEPYxgGK7YfBmDeFzv4Zp/D5IpEzl+j02DppmIA0j3w9u/xbh/dnS4RwRRV1PD3z7ebXY6IWygAiniYvWVHKamsBcBpwJzF67VxvbQ7a3aXcbiqjvAgP0YkRZpdzhkF+fvys0nJALyYu51DR2pNrkik7SkAiniYNcd6pvWICiU8yI/1+xy8tmKXyVWJnJ/m279X9LXj7+v5HzXfGRTPoAQbR2ob+POSLWaXI9LmPP+7UsTLrNnVFADH9Ynm51c3Nav900eF7C8/amZZIudliWv+n+e1fzkVHx8Ljx5rDr1g1W62Fh8xuSKRtqUAKOJhmndNGNatE7eOSGR4t05U1TXyeNYGkysTOTdbi4+w/VAVAb4+XJYcbXY552xUj85MSLHT6DT4fbYWYEnHpgAo4kGq6xooOFAJwNBunfDxsfDkjQPx87GQs/Eg2d8UmVyhyNk13/5N69mZsEA/k6s5Pz+/ui++x77fVh5bjCXSESkAiniQr/c4aHQaxIYHEW8LAqBPjJV7L+sBwBNZG6isqTezRJGzytnY9IuKJzZ/Pptediu3jkgEmppDO9UcWjooBUARD9K8AGRot4gWfdNmXtGbbp1DKKqo4U8fbTarPJGzKqmsZe2ecqB9BkCAWRP6EBrgy9d7Hby3/oDZ5Yi0CQVAEQ/SvABkaNdOLZ4P8vfld1MGAvDPFTv5+tgHrIin+bjgIIYBgxNsxIQHmV3OBYm2BnLfZT0BeDp7E7UNag4tHY8CoIiHMAzjuBHATid9fWzvKG5I7YKh3oDiwY7f/aM9u/vSHsSEB7K37CivLVcbJul4PC4AZmZmMmLECKxWK3a7nSlTplBYWNjimMWLFzNx4kSioqKwWCzk5+ef9nyGYXD11VdjsVh45513WnytrKyMjIwMbDYbNpuNjIwMysvLW/+iRM7BjkNVlFXXE+DnQ//48FMe89jkFCJC/Nl4oIJ5X+x0b4EiZ1Fd18DnWw8B7af9y+kEB/jyk/Sm5tD/t3QL5dV1Jlck0ro8LgDm5uYyffp08vLyyMnJoaGhgfT0dKqqqlzHVFVVMWbMGJ566qmznu+555477R6UU6dOJT8/n+zsbLKzs8nPzycjI6PVrkXkfKzZXQ7AwC42Av18T3lMVFggjx7rDfhMzmb2lFa7qzyRs/p08yFqG5x0jQyhT0yY2eVctO8OTaBvrJWKmgb+b+lWs8sRaVUetz4/Ozu7xZ/nzZuH3W5n9erVjBs3DsAV0nbu3HnGc3399dc888wzfPnll8TFxbX4WkFBAdnZ2eTl5TFy5EgAXnnlFdLS0igsLCQ5ObmVrkjk3Kw5rv/fmdw0PIH/rNnLqh2l/L93v+HVO0ac9pccEXc6/vZvR/g36etjYc41Kdz+6ipeW7GT29O607VziNllibQKjxsBPJHD4QAgMvL89pKsrq7m+9//Ps8//zyxsSffilixYgU2m80V/gBGjRqFzWZj+fLlpzxnbW0tFRUVLR4ireXbBSARZzzOYrHw5A0DCfD1YVlhCe+vV29AMV9Do5Olm5oC4ISU9j3/73iX9Ynm0t5R1Dca/P5DNYeWjsOjA6BhGMyePZuxY8cyYMCA83rtQw89xOjRo7n++utP+fWioiLsdvtJz9vtdoqKTv2BmpmZ6ZovaLPZSExMPK+aRE6nsqaewoPHGkB3PfMIIEAvexj3j29apfjEfzdQod6AYrLVu8ooq64nIsSfEd3P/m+4PXn0mhQsFvjfugOsPTZSL9LeeXQAnDFjBuvWrWPhwoXn9bqsrCyWLl3Kc889d8bjTnWLwjCM0966mDNnDg6Hw/XYs2fPedUlcjr5e8oxDEjoFIz9HFtn3D++Jz2iQimprOVpbVslJmu+/XtFsh0/X4/+aDlvKXHhfG9oAtDUHNow1Bxa2j+P/S6dOXMmWVlZLFu2jISEhPN67dKlS9m2bRsRERH4+fnh59c01fG73/0u48ePByA2NpaDBw+e9NqSkhJiYk59+yIwMJDw8PAWD5HWsGZXOXD2+X/HC/L35Xc3NPUGfGPlblbv0siEmMMwDHIKOkb7l9P5SXoyQf4+fLmzjA83nPzZIdLeeFwANAyDGTNmsHjxYpYuXUpSUtJ5n+PnP/8569atIz8/3/UAePbZZ5k3bx4AaWlpOBwOVq1a5XrdypUrcTgcjB49ulWuReRcrd596gbQZ5PWszPfG5aAYcCji9dTr96AYoItxUfYdbiaAD8fxvWJNrucNhFrC+LusU1bMv4+e5O+16Td87hVwNOnT2fBggW8++67WK1W13w8m81GcHAwAKWlpezevZv9+/cDuPoExsbGtnicqGvXrq5AmZKSwqRJk5g2bRovvfQSAPfccw/XXnutVgCLWzmdhmte0fmMADZ77JoUlm4qpvBgJX/7bIdrbqCIuzTf/h3TszOhgR73sdJq7hvfkze/3M2OQ1UsXLWbH6Z1N7skkQvmcSOAc+fOxeFwMH78eOLi4lyPRYsWuY7JysoiNTWVyZMnA3DrrbeSmprKiy++eF7v9cYbbzBw4EDS09NJT09n0KBBvP766616PSJns7XkCJU1DQT7+9I31nrer+8UGsBj1zT1Bvzzx5vZfVi9AcW9PnK1f2nfzZ/PJizQjwcn9AHguSVbtPhK2jWP+1XtXCbX3nHHHdxxxx0Xfd7IyEjmz59/XucRaW3N7V8GJ9ouePL8jUO78NaavSzfdphfvPsN//yRegOKexysqHHtTT0h5eTOCh3NrSMSmffFDraXVPHiJ9v42aS+ZpckckE8bgRQxNus3nVh8/+OZ7FY+N0NAwnw8+HTzSVkfb2/tcoTOaMlxxZ/DEmMOOcV7O2Zv68Pc47txvP3z3ewv/yoyRWJXBgFQBGTnesOIGeTFBXKzMt7AfCb9zbiqNbtKWl7x+/+4S0mpNi5JCmS2gYnf/yo8OwvEPFACoAiJiqvrmNbSdM+16kXMQLY7N7LetLLHsahI3U8lV1w0ecTOZMjtQ0s33oYgHQvCoAWi8U17/bttfv4Zp/D5IpEzp8CoIiJ1u4uB6BHVCiRoQEXfb4APx+ePNYbcOGqPXy5s/SizylyOp9uLqGu0Un3ziH0soeZXY5bDU6M4LrB8RgGZH6g5tDS/igAipioef5fa4z+NbskKZJbRzRtU/jo4vXUNahfmbSN42//euOio59OTCbA14cvth7mk80lZpcjcl4UAEVM1Frz/07086v7EhUWwJbiI7z86bZWPbcIQH2jk6WbioGO3/7ldBIjQ7hjTHcAMt8voEHNoaUdUQAUMUlDo5P8Y+0zhnaLaNVzR4QE8Mtr+wHwl6Vb2XGoqlXPL/LlzlIcR+uJDA1o9V9g2pPp43thC/Zn88Ej/Gf1XrPLETlnCoAiJik8WEl1XSPWQD9628+/AfTZXDc4nkt7R1HX4OQX76zXHCVpVc23f6/oa8fXx/tu/zazhfjzwJW9AXgmZzPVdQ0mVyRybhQARUzS3AB6SNeINvkAtVgs/HbKAAL9muYovb12X6u/h3gnwzC8sv3L6WSM6kbXyBCKK2t55dMdZpcjck4UAEVMsubYCuCLaQB9Nt06h/LghKbRid/+r4Cyqro2ey/xHpuKKtlbdpRAPx8u7R1ldjmmC/Dz4WeTmvaQf+nTbRRX1phckcjZKQCKmMS1A0gbz5+admkPkmOslFbV8eT76g0oF6959O/S3lGEBHjcjqKmmDwwjiGJEVTXNfJszhazyxE5KwVAEROUVNayu7Qai6VpC6225O/rw5M3NvUG/PfqvazYdrhN3086Pt3+PZnFYuEXk5uaQy/6cjdbDlaaXJHImSkAipiguf1Lb3sYtmD/Nn+/Yd06cdvIrgA89vZ6ahsa2/w9pWM64DjK+n0OLBa4oq8C4PGGd49kYv8YnAZkfrDJ7HJEzkgBUMQEbdX/70x+Nqkv0dZAth+q4oVl6g0oF2bJsdG/oV07EW0NNLkaz/PIpL74+VhYuqmY5VsPmV2OyGkpAIqYYE0b7AByNrZgfx7/TlNvwLmfbGNr8RG3vbd0HB/p9u8Z9YgOc422/+79ApxOtV8Sz6QAKOJmdQ1O1u1t2jze3Q10Jw+M4/LkaOoanTz2tnoDyvmprKknb3vTHFIFwNN74MreWAP92LC/gne/Vvsl8UwKgCJutvFABbUNTiJC/OkRFerW97ZYLPz6+gEE+/uyckcp/9bOBXIecjeXUN9o0CM6lJ7RYWaX47E6hwVy/+U9Afjjh5upqdecW/E8CoAibtZ8+3do105YLO7fQSExMoSHrmrqDfjk+wUcPlLr9hqkfdLq33N355gk4m1B7Cs/yrwvdppdjshJFABF3Gz17uYAGGFaDT8ak0RKXDjl1fX87n/qDShnV9/oZNmmYgDSFQDPKsjfl4cnNjWHfmHZVkrVhF08jAKgiJutdVMD6DPx9/Uh88aBWCyweO0+Pt+i1YpyZqt2lFJR00BUWABDEs37t9ueTBnShf7x4VTWNvCXj9UcWjyLAqCIGx1wHGW/owYfCwxOiDC1liGJEfxwVDcAfvHOes1TkjNqvv17Zd+YNtm7uiPy8bHw2DVNzaHn5+1ix6EqkysS+ZYCoIgbrdlVDkBKXDihgeZvofXwxGRiwgPZebia55duNbsc8VCGYbgC4ATd/j0vo3tFcXlyNA1Og6ez1RxaPIcCoIgbrT5uAYgnsAb586vr+gNNm9hv1vZVcgobD1Swr/woQf4+jO0VZXY57c6ca1LwscAH3xTx1c5Ss8sRARQARdzKjB1AzmZi/1gmpMRQ32jw6OL1alwrJ2ke/bu0dzTBAb4mV9P+9ImxcsuIRKBp5b36b4onUAAUcZOa+kY27G9qAO0pI4DQ1BvwV9f3JyTAl692lbHoqz1mlyQeRu1fLt5DE/oQEuDLmt3lfPBNkdnliCgAirjLN/sc1DcaRIUFkhgZbHY5LXSJCOYn6U0tKzLfL6C4ssbkisRT7Cs/yob9FfhY4Mq+drPLabfs4UFMu7QHAL/P3kRdg9PkisTbKQCKuMm38/8iTGkAfTa3p3VjQJdwKmoa+O176g0oTZYcG/0b1q0TncMCTa6mfbtnXA+irYHsOlzN/LxdZpcjXk4BUMRNPHH+3/H8fH3IvGEQPhbI+no/nxQWm12SeADd/m09oYF+zL6qDwB/WboFx9F6kysSb6YAKOIGhmGw+lgLGDMbQJ/NwAQbd4xOAuCX737D0Tr1BvRmjqP15G0/DMBV/WJNrqZjuGlYAn1iwiivrueFZWq9JOZRABRxg71lRzl0pBZ/XwsDu9jMLueMfpLeh3hbEHtKj/Jn7V7g1T4pLKbBadDLHkZSVKjZ5XQIfr4+zLm6qTn0vOU72VtWbXJF4q0UAEXcoHn+X794G0H+nt1GIzTQj19fPwCAv322nU1FFSZXJGbR7d+2MT45mtE9O1PX4OSPHxaaXY54KQVAETdwzf/zoPYvZzKhXwyT+sfS4DSYo96AXqmuwUluYQmgANjaLBYLj16TgsUC7+TvZ93ecrNLEi+kACjiBq4VwN0izC3kPDxxXX/CAv1Yu7ucN1btNrsccbO87YeprG0g2hrIEJP3re6IBnSxccOQLoCaQ4s5FABF2lhVbQObipq2WPPUFcCnEmsL4qcTm3oDPv3BJg5WqDegN3Ht/Ztix8fH89oWdQQ/mZhMoJ8PedtL+bhAq+7FvRQARdrY13vLaXQaxNmCiLN5VgPos/nBqG4MToygsraBX/93o9nliJsYhsGSAs3/a2tdIoK5c2zTqvvMDwpoaFRzaHEfjwuAmZmZjBgxAqvVit1uZ8qUKRQWtpwku3jxYiZOnEhUVBQWi4X8/PyTznPvvffSs2dPgoODiY6O5vrrr2fTpk0tjikrKyMjIwObzYbNZiMjI4Py8vI2vDrxRmt3lwOe3f7ldHx9LGTeMBBfHwv/W3+ApZsOml2SuME3+yo44KghJMCX0T2jzC6nQ7t/fE8iQwPYVlKlbRjFrTwuAObm5jJ9+nTy8vLIycmhoaGB9PR0qqqqXMdUVVUxZswYnnrqqdOeZ9iwYcybN4+CggI+/PBDDMMgPT2dxsZv+5pNnTqV/Px8srOzyc7OJj8/n4yMjDa9PvE+3+4A0v4CIEC/+HDuPjZK8ct3NlBd12ByRdLWcjY27VU7rne0x69ab+/Cg/x58MreADybs5kjtfr+EvewGB4+87SkpAS73U5ubi7jxo1r8bWdO3eSlJTE2rVrGTJkyBnPs27dOgYPHszWrVvp2bMnBQUF9OvXj7y8PEaOHAlAXl4eaWlpbNq0ieTk5LPWVlFRgc1mw+FwEB4efsHXKB2XYRik/iaH8up63pk+hiGJEWaXdEGq6xq46plP2Vd+lGmXJvHY5H5mlyRtaNJzn7KpqJI/3TSY7w5LMLucDq++0Un6s5+y41AVD1zRi9npZ//8kYujz28PHAE8kcPhACAyMvKCz1FVVcW8efNISkoiMTERgBUrVmCz2VzhD2DUqFHYbDaWL19+yvPU1tZSUVHR4iFyJtsPVVFeXU+gnw/94trvD5mQAD9+e0NTb8BXv9jJN/scJlckbWVPaTWbiirx9bFwRV+72eV4BX9fHx6Z1BT6Xv5sO0UOLbiStufRAdAwDGbPns3YsWMZMGDAeb/+hRdeICwsjLCwMLKzs8nJySEgIACAoqIi7PaTf7jZ7XaKiopOeb7MzEzXfEGbzeYKkyKns+bY7d9BCTYC/Dz62+2sLk+2M3lQHI1Og0ffXk+jegN2SM2rf4d360Sn0ACTq/EeE/vHMrxbJ2rqnTyTo+bQ0vY8+hNpxowZrFu3joULF17Q62+77TbWrl1Lbm4uvXv35uabb6am5tvfrCyWk1sbGIZxyucB5syZg8PhcD327NGEXTmz5gbQ7XX+34kev7Yf1iA/1u118PqKnWaXI21Au3+Yw2Kx8Ojkpi3i/r16r3bgkTbnsQFw5syZZGVlsWzZMhISLmwOis1mo3fv3owbN47//Oc/bNq0ibfffhuA2NhYDh48eUVjSUkJMTGn/sEXGBhIeHh4i4fImazZVQ60zxXAp2IPD+KRSX0B+MOHhRxwHDW5ImlN5dV1rNpZCkB6v1iTq/E+Q7t2YvKgOAwDMt/fdPYXiFwEjwuAhmEwY8YMFi9ezNKlS0lKSmrVc9fW1gKQlpaGw+Fg1apVrq+vXLkSh8PB6NGjW+09xXtV1NSzubipAXRHGQEEmHpJV4Z2jaCqrpEnsjaYXY60omWFxTQ6DZJjrHTtHGJ2OV7pkYl98fe1kLu5hM+2lJhdjnRgHhcAp0+fzvz581mwYAFWq5WioiKKioo4evTbkYbS0lLy8/PZuLGpMW1hYSH5+fmuuXvbt28nMzOT1atXs3v3blasWMHNN99McHAw11xzDQApKSlMmjSJadOmkZeXR15eHtOmTePaa689pxXAImeTv7scw4CukSFEWwPNLqfV+PhYyLxxEH4+Fj7ccJCPNpx6zqy0P7r9a76unUPIGNUdgN/9r0BzbaXNeFwAnDt3Lg6Hg/HjxxMXF+d6LFq0yHVMVlYWqampTJ48GYBbb72V1NRUXnzxRQCCgoL47LPPuOaaa+jVqxc333wzoaGhLF++vMXCjzfeeIOBAweSnp5Oeno6gwYN4vXXX3fvBUuH9W3/vwhzC2kDybFW7hnXA4DHszaod1kHUNvQSG5h04iTAqC5Zl7Ri/AgPzYVVbJ4zV6zy5EOyuP7AHoy9RGSM8n4+0o+23KI31zfn4y07maX0+pq6htJf/ZTdpdW86Mx3Xn8O/3NLkkuwrLCYn4070tiwgNZ8fMrtf+vyV7+dBtPvr8JuzWQd2eMaXfbSHo6fX574AigSEfgdBrkH9sCLrUDzf87XpC/L7+d0tSe6Z/Ld7Jub7m5BclFab79OyElRuHPA/wwrTs9okIprqxl6isr1RtQWp0CoEgb2FJ8hMraBkICfOkbazW7nDYzrk801w+Jx2nAnMXrtZl9O+V0GizR/D+PEuTvy+t3jyQxMpgdh6qY+koeBysUAqX1KACKtIHm+X+DEyLw8+3Y32a/vLYftmB/Nuyv4B/Ld5pdjlyA9fscFFfWEhboR1rPzmaXI8d0iQhm4bRRJHQKZvuhKr7/Sh7FCoHSSjr2J5OISZobQA/rIP3/ziQqLJA5Vzf1BnwmZzP7ytUbsL1pvv17WZ9oAv18Ta5GjpfQKYSF00bRJSKY7SXHQmClQqBcPAVAkTbQvAXc0G4R5hbiJjcPT+SS7pFU1zXy/975Bq0ta1/U/sWzJUY2hcB4WxDbSqqY+spKSiprzS5L2jkFQJFWVlpVx/ZDVQCkJnb8EUBo6g345I0D8Pe18PGmYrK/UW/A9mL34WoKD1bi62NhfHK02eXIaXTtHMLCe0YRZwtia/ERbvtbHoeOKATKhVMAFGlla4/d/u0RHUqn0ACTq3GfXnYr91/WE2jqDVhRU29yRXIuPtrYFNYv6R5JRIj3/Httj7p1DmXhtFHEhgex+eARbntlJYcVAuUCKQCKtDLX/L8O2v7lTH58eS+SjrWu+OOHhWaXI+dAt3/bl+5RoSy8ZxQx4YEUHqzktr+tpLSqzuyypB1SABRpZa4dQLxgAciJgvx9+d2x3oCv5+1yhWHxTGVVdXy5sxRQAGxPkqKaRgLt1kA2FVUy9ZU8yhQC5TwpAIq0ooZGJ1/vcQDesQL4VEb3iuLGoV0wDHh08Xrq1RvQYy3dVIzTgL6xVhIjQ8wuR85Dj+gwFkwbRfSxEHjb31ZSXq0QKOdOAVCkFW0qquRofSPWID96RYeZXY5pfjG5H51C/NlUVMnfP99hdjlyGs23f9M1+tcu9bKHsXDaSKLCAtl4oEIhUM6LAqBIK2q+5ZnatZNXb6cVGRrAo9ekAPDcks3sKa02uSI5UU19I59uKQHgqn6xJlcjF6qX3XosBAawYX8FGX9fhaNaC7Dk7BQARVqRa/5f1whzC/EA3xuWwKgekdTUO/mFegN6nOXbDlFd10icLYgBXcLNLkcuQu8YKwumjaJzaADr9znIeHUljqMKgXJmCoAircibdgA5G4vFwu9uGEiArw+5m0t4b90Bs0uS4zTf/p2QEoPF4r2j1R1Fn2MhMDI0gHV7Hfzw7yvViknOSAFQpJUUV9awp/QoFgsMSYwwuxyP0DM6jOmX9wLgV//dqFtTHsLpNFhSUAxo9W9HkhxrZcG0kXQK8efrvQ5++PdVCoFyWgqAIq1kza5yAJJjrFiD/M0txoPcN74HPaNDOXSklt9/uMnscgTI31tOSWUt1kA/RvXobHY50or6xobzxt2jiAjxJ39PObe/uopKhUA5BQVAkVZy/AIQ+Vagny9P3jAQgAUrd/PVsb5zYp7m27+XJUcT4KePgY6mX3w4b9w9EluwP2t3l3PHvC85UttgdlniYfSdL9JK1uzS/L/TGdmjMzcPTwDg0bfXU9eg3oBm0u4fHV//eJsrBK7eVcYdr65SCJQWFABFWkFdg5N1+5oaQGsF8Kk9ek0KnUMD2HzwCK98tt3scrzWjkNVbC0+gp+PhfHJdrPLkTY0oIuN+XeNJDzIj692lXHnvC+pUgiUYxQARVrBhv0O6hqcdArxJykq1OxyPFJESAC/uLapN+BfPt7CzkNVJlfknXI2FgEwqkdnbMGaq9rRDUywMf/ukViD/Fi1s5Qf/eNLqusUAkUBUKRVfNv/r5NaapzBlCFdGNsritoG9QY0i27/ep9BCRG8ftdIrIF+rNpRyp3/+JKjdY1mlyUmUwAUaQVrd5cDMFTz/87IYrHw2ykDCPTz4fOth3g3f7/ZJXmVw0dqXb+sTFAA9CpDEiN47a5LCAv0I297KXf9UyHQ2ykAirSC40cA5cy6R4XywJW9AfjNexu1d6kbfbypGKcB/ePD6RIRbHY54mapXTvxzzubQuDybYe5+7UvqalXCPRWCoAiF2l/+VGKKmrw9bEwONFmdjntwrRLe9AnJozDVXVkvq/egO6i278yrFsn/nnnCEIDfPli62GmvfaVQqCXUgAUuUjNo38pcVZCAvxMrqZ9CPDzcfUGXPTVHlZuP2xyRR3f0bpGPttSAigAerth3SL5x52XEBLgy2dbDikEeikFQJGL5Nr/V7d/z8vw7pFMHdkVaOoNWNugD6C29PnWQ9TUO+kSEUy/uHCzyxGTjegeyT9+9G0IvPf11QqBXkYBUOQiNTeA1gKQ8/fIxL5EhQWyraSKFz9Rb8C21Nz+5ap+MVqpLgBckhTJq3eMINjfl9zNJdw/f7V+EfMiCoAiF6GmvpEN+ysALQC5ELYQfx7/Tj8A/rpsK9tLjphcUcfU6DT4uKAY0O1faWlUj868escIgvx9WFZYwv3z1ygEegkFQJGLsG6vgwanQbQ1kIROWlV5Ia4dFMdlfaKpa3Ty2NvqDdgW1u4u43BVHeFBflySFGl2OeJh0np25tXbm0Lg0k3FTH9jjbZr9AIKgCIXoXn+39CuEbqtdoGaewMG+fuwYvth3lqzz+ySOpzm1b+X97Xj76sf+3Ky0b2i+PvtIwj082FJQTHTFygEdnT6SSByEZpXAA/T/L+LkhgZwqwJfQD43f82Ulql3oCtSe1f5FyM6RXF324fToCfDzkbDzJjwRrqGxUCOyoFQJELZBgGa3erAXRruWtsEn1jrZRV1/O7/xWYXU6HsbX4CNsPVeHva+GyPtFmlyMe7tLe0bzyw6YQ+NHGg8xcsFYhsINSABS5QLtLqzl0pA5/XwsDuqgB9MXy9/Uh88aBWCzw1pq9LN96yOySOoTm0b+0nlFYg/xNrkbag8v6RPNyxjACfH3I3lDEg28qBHZECoAiF6h5/l//eBtB/r4mV9MxpHbtRMaobgA89s436kvWCo5v/yJyrsYn23npWAh8f30Rsxbl06AQ2KEoAIpcIM3/axsPT0zGbg1kx6EqnsnZrFXBF6Gkspa1e8oBmJBiN7cYaXcu72vnxYyh+Pta+N+6Azz0r68VAjsQjwuAmZmZjBgxAqvVit1uZ8qUKRQWFrY4ZvHixUycOJGoqCgsFgv5+fktvl5aWsrMmTNJTk4mJCSErl278sADD+BwOFocV1ZWRkZGBjabDZvNRkZGBuXl5W18hdJRrNlVDmj+X2sLD/LnV9f1B+DlT7fz6/c20uhUCLwQHxccxDBgYBcbcTa1KZLzd0XfGObeNgx/Xwv//Xo/sxUCOwyPC4C5ublMnz6dvLw8cnJyaGhoID09naqqKtcxVVVVjBkzhqeeeuqU59i/fz/79+/nj3/8I+vXr+cf//gH2dnZ3HXXXS2Omzp1Kvn5+WRnZ5OdnU1+fj4ZGRlten3SMVTVNrCp6FgD6G4R5hbTAV09MI7HrkkBYN4XO5mxYI1uB1+AJQVa/SsXb0K/GP46dSh+Phayvt7Pw//+Wr+UdQAWw8Pvr5SUlGC328nNzWXcuHEtvrZz506SkpJYu3YtQ4YMOeN5/v3vf/ODH/yAqqoq/Pz8KCgooF+/fuTl5TFy5EgA8vLySEtLY9OmTSQnJ5+1toqKCmw2Gw6Hg/Bw7a3pTZZvPcTUv60k3hbE8jlXml1Oh5X19X4e/tfX1DU6Gd6tE6/8cDidQgPMLqtdqK5rIPXXOdQ2OPngwUtJ0f6/cpGyvylixoI1NDgNbkztwh9uGoyvT/vsf6rPbw8cATxR823byMiL617f/D/Zz88PgBUrVmCz2VzhD2DUqFHYbDaWL19+ynPU1tZSUVHR4iHeydUAWvP/2tR1g+N57a5LCA/y46tdZXz3xeXsKa02u6x24bMth6htcJLQKZi+sVazy5EOYNKAWJ6fmoqvj4XFa/fxs/+s00hgO+bRAdAwDGbPns3YsWMZMGDABZ/n8OHD/OY3v+Hee+91PVdUVITdfvKkaLvdTlFR0SnPk5mZ6ZovaLPZSExMvOCapH1rXgCi+X9tb1SPzvzn/tHE24LYXlLFDS8sZ/1ex9lf6OWOb/6sXWqktUwaEMf/fb8pBL61Zi8/f2sdToXAdsmjA+CMGTNYt24dCxcuvOBzVFRUMHnyZPr168fjjz/e4mun+qFoGMZpf1jOmTMHh8PheuzZs+eC65L2y+k0XCsrtQLYPfrEWFn84zH0jbVy6Egtt7y8gk8Ki80uy2M1Og2Wbmr6+9H8P2lt1wyM48+3DsHXx8K/V+9lzuL1CoHtkMcGwJkzZ5KVlcWyZctISEi4oHNUVlYyadIkwsLCePvtt/H3/7YJamxsLAcPHjzpNSUlJcTEnPoHZmBgIOHh4S0e4n22H6qivLqeQD8fzatyo1hbEP++L42xvaKormvkrn9+xb++1C9hp7J6VxmlVXXYgv25pPvFTZ8ROZVrB8Xz7C1D8LHAoq/28OjbCoHtjccFQMMwmDFjBosXL2bp0qUkJSVd0HkqKipIT08nICCArKwsgoKCWnw9LS0Nh8PBqlWrXM+tXLkSh8PB6NGjL+oapGNrnv83OCGCAD+P+xbq0KxB/rx6xwhuTO1Co9PgZ2+t47kl6hV4oubmz1f0tePnq3+j0jauG/xtCHzzyz384t1vFALbET+zCzjR9OnTWbBgAe+++y5Wq9U1H89msxEc3NTHqrS0lN27d7N//34AV5/A2NhYYmNjqaysJD09nerqaubPn99iwUZ0dDS+vr6kpKQwadIkpk2bxksvvQTAPffcw7XXXntOK4DFe605Nv8vVe1fTBHg58Ofbh5MfEQwzy/bynNLtrC//Ci/u2Eg/go7GIbRYv6fSFu6fkgXDANm/yufBSt342OB31w/QPNO2wGP+2k5d+5cHA4H48ePJy4uzvVYtGiR65isrCxSU1OZPHkyALfeeiupqam8+OKLAKxevZqVK1eyfv16evXq1eI8x8/be+ONNxg4cCDp6emkp6czaNAgXn/9dfdesLQ7zSOAw7QAxDQWi4WHJybzuxsG4GOBf321l7v/+RVVtQ1ml2a6rcVH2Hm4mgBfH8b1iTa7HPECU1K78MebBmOxwPy83TyetUGj8u2Ax/cB9GTqI+R9HEfrGfyrjwD46hcTiAoLNLkiWbLxIDMWrqGm3smALuG8escI7Nags7+wg/rrsq384cNCxidH848fXWJ2OeJF/rN6Lz/9z9cYBtwxujuPf6efx44E6vPbA0cARTxZ/rHVv906hyj8eYgJ/WJ48540OocG8M2+Cm58YTnbSo6YXZZpdPtXzPK9YQn8/ruDsFjgH8t38uv3Nmok0IMpAIqcB/X/80xDEiN46/7RdO8cwt6yo3x37nK+2llqdlluV1xR4/olZUKKAqC4383DE3nqxoFA0zaOv3mvQCHQQykAipyHtdoBxGN1jwrlrftHMyQxgvLqeqb+bSUfrD9gdllutaSgqfff4MQIYsK99za4mOuWEV3JPBYCX/1iB0++rxDoiRQARc5Ro9Ng7e5yAIZ2jTC1Fjm1zmGBLJw2igkpMdQ1OPnxgjXM+2KH2WW5TXP7l3Td/hWTff+Srjx5Q1MIfOWzHTz1wSaFQA+jAChyjrYUV3KktoHQAF+SY7S3qqcKDvDlxR8M5QejumIY8Kv/buR3/9vY4fuTVdU28MW2w4Dm/4lnmDqyK7+Z0rSN60ufbuf32YUKgR5EAVDkHDXP/xucGKHmuh7Oz9eH31w/gJ9Naurp+cpnO3jgzbXUNjSaXFnb+XRzCXUNTrp1DqG3PczsckQAyBjVjV9f3x+AF3O38YcPFQI9hT7FRM7Rml3lgPb/bS8sFgs/Ht+LZ28ZjL+vhffWHeCHf1+Fo7re7NLahGv1b0qMx7beEO/0w7TuPPGdfgC88Mk2/vSRdu/xBAqAIueouQG0VgC3LzekJvCPH12CNdCPlTtK+d6Ly9lXftTsslpVQ6OTpYVNC0B0+1c80R1jkvh/1zaFwOeXbeXZJVtMrkgUAEXOQWlVHTsOVQGQqgUg7c6YXlH86740YsOD2FJ8hBv++gUb91eYXVar+XJnGeXV9XQK8dcItXisO8cm8YvJKQD85eMtPLdks8kVeTcFQJFz0Lz/b8/oUCJCAkyuRi5ESlw4i388mj4xYRRX1nLzSyv4fMshs8tqFc23f6/oG6P5qeLR7r60B49d0xQCn1uyhb98rJFAs+gnhcg5cO3/q9GVdi0+Iph/3zeaUT0iOVLbwB3zVrF4zV6zy7oohmGQU9DU/kW3f6U9mDauB3Ou7gvAMzmbeX6pQqAZFABFzoF2AOk4bMH+/PPOS7hucDwNToPZ//qavy7b2m4npRcerGRP6VEC/XwY1yfK7HJEzsm9l/XkkUlNIfCPH23mr8u2mlyR91EAFDmL+kYn6/Y6AI0AdhSBfr48d8sQ7r2sBwB/+LCQx975hoZGp8mVnb+cDU23f8f2iiIkwM/kakTO3f3je/LTiU2tmv7wYSFzP9lmckXeRQFQ5Cw2HajkaH0j4UF+9IxWf7WOwsfHwpyrU/jVdf2xWGDByt3c+/pqqusazC7tvOQUHGv/otu/0g5Nv7wXD6f3AeD32Zt4+VOFQHdRABQ5i+b5f6ldO+Hjo/5qHc3to7sz97ZhBPr58PGmYr7/ykoOHak1u6xzUuSoYd1eBxYLXJmiACjt04wrejP7qqYQ+OT7m/jbZ9tNrsg7KACKnIXm/3V8kwbEsmDaKDqF+PP1nnK+O3e5q+2PJ2se/UtNjCDaGmhyNSIX7oEre/Pglb0B+O3/Cvj7596zh7dZFABFzkIrgL3DsG6deOv+0SRGBrPrcDXfnbvc9f/eU7l2/+gXa3IlIhdv1oTePHBFLwB+895G5n2hENiWFABFzqC4ooa9ZUexWGBwos3scqSN9YgOY/H9YxiUYKO0qo6pr+S5QpanqaypZ8W2pj6Gmv8nHYHFYuGhq/ow4/KmEPir/27kn8t3mltUB6YAKHIGzSNAyTFWrEH+Jlcj7hBtDWThtFFcnhxNTb2Te1//itfzdpld1klyN5dQ32iQFBVKz+hQs8sRaRUWi4WfpPfhx+N7AvB41gZeW7HT3KI6KAVAkTNwzf/T7V+vEhroxys/HM73L0nEacAv3/mG32dvwun0nF6B397+jcFi0eIk6TgsFgs/nZjMfZc1hcC6hvbXnqk9UNMokTNYs7scgGFaAOJ1/Hx9ePKGgcTbgvlTzmbmfrKNA+VHefp7gwnwM/d35/pGJ8s2FQO6/Ssdk8Vi4ZFJyVzR184lSZFml9MhaQRQ5DRqGxpZf6wBtEYAvZPFYmHmlb35w/cG4edj4Z38/dwxbxUVNfWm1rVqRykVNQ10Dg3Q6nTpsCwWi8JfG1IAFDmNDfsrqGt0EhkaQPfOIWaXIya6aXgir94xgtAAX5ZvO8zNL67ggOOoafU03/69oq8dX/WmFJELoAAochprXP3/IjTHShjXJ5pF96YRbQ1kU1ElN76wnMKiSrfXYRhGi/l/IiIXQgFQ5DSaVwDr9q80G9DFxts/Hk0vexgHHDV878XlLD/WisVdCg5Usq/8KEH+PlzaO9qt7y0iHYcCoMgpGIahHUDklBI6hfCf+9K4pHsklTUN3P7qKt7N3+e2928e/RvbK5rgAF+3va+IdCwKgCKnsN9Rw8GKWnx9LAxOiDC7HPEwESEBvHbXJUweGEd9o8GDb+bzYu42DKPt28TkFBQBkK7bvyJyERQARU6hefSvX1y4RlnklIL8ffm/76dy19gkAJ76YBOPZ22gsQ17Be4vP8o3+yqwWOCKFHubvY+IdHwKgCKn0LwARPv/ypn4+Fj45bX9+MXkFCwWeG3FLu6fv5qa+sY2eb8lBU23f4d17URUWGCbvIeIeAcFQJFTaF4Akto1wtxCpF24+9IePP/9oQT4+fDRxoNMfSWP0qq6Vn8frf4VkdaiAChygqN1jWzcXwFoBFDO3eRBccy/ayS2YH/W7C7nu3OXs/twdaudv6KmnrzthwEFQBG5eAqAIidYt7ecBqeB3RpIl4hgs8uRduSSpEjeuj+NLhHB7DhUxY1zv+DrPeWtcu5PCkuobzToGR1Kj+iwVjmniHgvBUCRE7j2/+3WSQ2g5bz1slt5+8ej6R8fzqEjddz6cp5r396L8e3t39iLPpeIiAKgyAnU/08ulj08iEX3pjGuTzRH6xu5+7WvWLhq9wWfr67BySfHQqRu/4pIa1AAFDmOYRis1Q4g0grCAv34++3DuWlYAo1OgzmL1/PMR4UX1Ctw5Y7DVNY2EBUWSGpiROsXKyJex+MCYGZmJiNGjMBqtWK325kyZQqFhYUtjlm8eDETJ04kKioKi8VCfn7+Sed5+eWXGT9+POHh4VgsFsrLy086pqysjIyMDGw2GzabjYyMjFMeJ95j1+FqDlfVEeDrw4Au4WaXI+2cv68PT39vEA9e2RuAvyzdysP/Xkd9o/O8ztN8+3dCih0fH01LEJGL53EBMDc3l+nTp5OXl0dOTg4NDQ2kp6dTVVXlOqaqqooxY8bw1FNPnfY81dXVTJo0iUcfffS0x0ydOpX8/Hyys7PJzs4mPz+fjIyMVr0eaV+a278M6BJOoJ8aQMvFs1gsPHRVH566cSC+PhbeWrOXO//xJUdqG87p9YZhsETtX0SklfmZXcCJsrOzW/x53rx52O12Vq9ezbhx4wBcIW3nzp2nPc+sWbMA+OSTT0759YKCArKzs8nLy2PkyJEAvPLKK6SlpVFYWEhycvLFXYi0S5r/J23l1ku6EhMexI/fWMNnWw5x84sr+MePRmAPDzrj6zbsr2C/o4Zgf1/G9IpyU7Ui0tF53AjgiRwOBwCRkZGtet4VK1Zgs9lc4Q9g1KhR2Gw2li9ffsrX1NbWUlFR0eIhHcvxK4BFWtvlfe0suncUUWEBbDxQwQ0vLGdrceUZX/PRsdG/cX2iCPLXqLSItA6PDoCGYTB79mzGjh3LgAEDWvXcRUVF2O0n76Vpt9spKio65WsyMzNd8wVtNhuJiYmtWpOY60htA4VFTaFeC0CkrQxKiGDx/WPoERXKvvKj3PjCclbtKD3t8Wr/IiJtwaMD4IwZM1i3bh0LFy5sk/OfqsebYRin7f02Z84cHA6H67Fnz542qUvM8fWecpwGdIkIJuYst+VELkbXziG8df9ohnXrREVNAz/420r+t+7AScftKa2m4EAFPha4ou/Jv7CKiFwojw2AM2fOJCsri2XLlpGQkNDq54+NjeXgwYMnPV9SUkJMzKknWgcGBhIeHt7iIR2Ha/6fRv/EDTqFBvDG3SOZ2D+GukYnMxau4W+fbW9xzJKCpp9Rw7tHEhkaYEaZItJBeVwANAyDGTNmsHjxYpYuXUpSUlKbvE9aWhoOh4NVq1a5nlu5ciUOh4PRo0e3yXuKZ2teATysa4S5hYjXCPL35YXbhnF7WjcMA377vwJ+/d+NOJ1NvQKbb/+ma/WviLQyj1sFPH36dBYsWMC7776L1Wp1zcez2WwEBzfty1paWsru3bvZv38/gKtPYGxsLLGxTfNkioqKKCoqYuvWrQCsX78eq9VK165diYyMJCUlhUmTJjFt2jReeuklAO655x6uvfZarQD2Qk6nwRqNAIoJfH0sPHFdf7p0CubJ9zfx6hc7KKo4yhPf6c/KY3MD1f5FRFqbx40Azp07F4fDwfjx44mLi3M9Fi1a5DomKyuL1NRUJk+eDMCtt95KamoqL774ouuYF198kdTUVKZNmwbAuHHjSE1NJSsry3XMG2+8wcCBA0lPTyc9PZ1Bgwbx+uuvu+lKxZNsP3SEipoGgvx9SInTrX1xL4vFwj3jevKX76cS4OvD++uLmPx/n9PoNOgTE0a3zqFmlygiHYzFuJB9iQSAiooKbDYbDoejVecDllXV8cpn27lxaBd62a2tdl45vUVf7uaRt9ZzSVIk/7o3zexyxIut2HaYe17/isqapkbR0y/vyU8n9jW5KpGOpa0+v9sTjxsBFPjlu9/wwifb+PPHW80uxWus2VUOqP+fmC+tZ2feun80cbYgfCxw7aB4s0sSkQ5IAdAD3T++JwDvrdtPYdGZm8RK61i9WzuAiOfoE2MlZ/ZlLJl9maYkiEibUAD0QP3jbVw9IBbDgD9/vNnscjo8R3U9W4uPADBUK4DFQ4QF+tEjOszsMkSkg1IA9FCzJvTBYoH31xexcb+2nGtLa/Y0jf517xxC57BAk6sRERFpewqAHio51srkgXEAPLdEo4Btaa3av4iIiJdRAPRgsyb0xmJp2gx+/V6H2eV0WJr/JyIi3kYB0IP1slu5fnDTCkCNAraNRqdB/u5yQCuARUTEeygAergHruyNjwU+3lRM/p5ys8vpcAqLKqmqayQs0I8+Meq5KCIi3kEB0MP1iA7jhtQEAJ7N0Shga2ve/3dIYgS+PhaTqxEREXEPBcB24IEre+HrYyF3cwmrd5WaXU6H4tr/V+1fRETEiygAtgPdOofyvaHNo4BbTK6mY2keAdQKYBER8SYKgO3EjCt64e9r4fOth1i5/bDZ5XQIh47UsvNwNQCpiQqAIiLiPRQA24nEyBBuHp4IwDM5mzEMw+SK2r+1x1b/9raHYQvxN7cYERERN1IAbEemX96LAF8fVu4oZcU2jQJerNW71P9PRES8kwJgOxIfEcz3L9EoYGtpnv+n/n8iIuJtFADbmR9f3osAPx++2lXGZ1sOmV1Ou1Xf6GTd3nIAhnaLMLUWERERd1MAbGdiwoP4wchugEYBL0bBgQpq6p3Ygv3pERVmdjkiIiJupQDYDt03vgdB/j7k7ynnk8ISs8tpl5rn/6V2jcBHDaBFRMTLKAC2Q3ZrED9M6w5oFPBCrWne/1cLQERExAspALZT947rQUiAL+v3OVhSUGx2Oe2OawcQLQAREREvpADYTnUOC+T20d2BplFAp1OjgOeqyFHDvvKj+FhgcGKE2eWIiIi4nQJgO3bPpT0IDfCl4EAFH20sMrucdqO5/UtybDhhgX4mVyMiIuJ+CoDtWKfQAO4cmwQ07RGsUcBz03z7d5jav4iIiJdSAGzn7h7bA2ugH4UHK3n/mwNml9MurN6tHUBERMS7KQC2c7YQf+66tGkU8LklW2jUKOAZ1dQ3smFfBaAdQERExHspAHYAd45Nwhbsz9biI7y3br/Z5Xi0Dfsd1DU66RwaQNfIELPLERERMYUCYAcQHuTPPeN6APDnJVtoaHSaXJHnWrOrHGhq/2KxqAG0iIh4JwXADuL20d3pFOLP9kNVvJOvUcDTad4BRPP/RETEmykAdhBhgX7ce1lPAP7y8RbqNQp4EsMwXAtANP9PRES8mQJgB/LDtG50Dg1gd2k1i9fsNbscj7O37CgllbX4+VgYlGAzuxwRERHTKAB2ICEBftw/vnkUcCt1DRoFPF5zA+j+8eEE+fuaXI2IiIh5FAA7mNtGdiPaGsi+8qP8e/Ues8vxKM0NoFM1/09ERLycAmAHExzgy4+PjQI+v3QrtQ2NJlfkOdbsLgc0/09EREQBsAP6/iVdiQkP5ICjhkVfahQQoLqugY0HmhpAD1UAFBERL+dxATAzM5MRI0ZgtVqx2+1MmTKFwsLCFscsXryYiRMnEhUVhcViIT8//6Tz1NbWMnPmTKKioggNDeW6665j796WCyPKysrIyMjAZrNhs9nIyMigvLy8Da/OPYL8fZlxeS8A/rpsKzX1GgVct9dBo9MgNjyIeFuQ2eWIiIiYyuMCYG5uLtOnTycvL4+cnBwaGhpIT0+nqqrKdUxVVRVjxozhqaeeOu15Zs2axdtvv82bb77J559/zpEjR7j22mtpbPw2DE2dOpX8/Hyys7PJzs4mPz+fjIyMNr0+d7l5RCLxtiAOVtSyYOVus8sxXfMCkKHdItQAWkREvJ7FMAyP3jy2pKQEu91Obm4u48aNa/G1nTt3kpSUxNq1axkyZIjreYfDQXR0NK+//jq33HILAPv37ycxMZH333+fiRMnUlBQQL9+/cjLy2PkyJEA5OXlkZaWxqZNm0hOTj5rbRUVFdhsNhwOB+Hh4a130a1kwcrdPPr2eqLCAvnsZ5cTHOC9K1/v/ueXLCko5heTU7j70h5mlyMiIiby9M9vd/C4EcATORwOACIjI8/5NatXr6a+vp709HTXc/Hx8QwYMIDly5cDsGLFCmw2myv8AYwaNQqbzeY6pr373rAEEjoFc+hILfPzdpldjmkMw3AtANH8PxEREQ8PgIZhMHv2bMaOHcuAAQPO+XVFRUUEBATQqVPLD/uYmBiKiopcx9jt9pNea7fbXcecqLa2loqKihYPTxbg58MDV/YG4MXcbVTVNphckTl2Hq6mtKqOAD8f+sd75296IiIix/PoADhjxgzWrVvHwoULW+V8hmG0mP91qrlgJx5zvMzMTNeCEZvNRmJiYqvU1ZZuTO1Ct84hHK6q47UV3jkK2Nz/b2AXG4F+3nsbXEREpJnHBsCZM2eSlZXFsmXLSEhIOK/XxsbGUldXR1lZWYvni4uLiYmJcR1z8ODBk15bUlLiOuZEc+bMweFwuB579nh+ixU/Xx8ePDYK+NKn26isqTe5IvfT/r8iIiIteVwANAyDGTNmsHjxYpYuXUpSUtJ5n2PYsGH4+/uTk5Pjeu7AgQN88803jB49GoC0tDQcDgerVq1yHbNy5UocDofrmBMFBgYSHh7e4tEeXDc4nh7RoZRX1/OPL3aaXY7bNY8ADu0aYW4hIiIiHsLjAuD06dOZP38+CxYswGq1UlRURFFREUePHnUdU1paSn5+Phs3bgSgsLCQ/Px819w9m83GXXfdxU9+8hM+/vhj1q5dyw9+8AMGDhzIhAkTAEhJSWHSpElMmzaNvLw88vLymDZtGtdee+05rQBuT44fBXzls+04jnrPKGBlTT2FBysBGKot4ERERAAPDIBz587F4XAwfvx44uLiXI9Fixa5jsnKyiI1NZXJkycDcOutt5KamsqLL77oOubZZ59lypQp3HzzzYwZM4aQkBD++9//4uv77RywN954g4EDB5Kenk56ejqDBg3i9ddfd9/FutG1g+LpbQ+joqaBVz/fYXY5bvP1HgeGAQmdgrGHqwG0iIgItIM+gJ6svfUR+t+6A0xfsAZroB+fPXI5ESEBZpfU5v68ZAvPLtnM9UPi+fOtqWaXIyIiHqC9fX63BY8bAZS2c/WAWPrGWqmsbeBvn3nHKKBrBxDd/hUREXFRAPQiPj4WZk3oA8C8L3ZQWlVnckVty+k0XAFQK4BFRES+pQDoZSb2j6F/fDhVdY28/Ol2s8tpU9tKjlBZ00Cwvy99Y61mlyMiIuIxFAC9jMVi4aFjo4D/XL6TQ0dqTa6o7aw+1v5lcKINP1/9UxcREWmmT0UvdGWKnUEJNo7WN/JS7jazy2kzmv8nIiJyagqAXshisfDQVU2jgK+t2EVxRY3JFbWN5hFAzf8TERFpSQHQS43vE01q1whqG5zM7YCjgOXVdWwrqQIgVSOAIiIiLSgAeimLxcJPrmra8eSNlbspcnSsUcC1u8sB6BEVSmRox+93KCIicj4UAL3YmF6duaR7JHUNTl74ZKvZ5bSq5vl/Gv0TERE5mQKgFzt+LuCbq/awr/zoWV7Rfmj+n4iIyOkpAHq5tJ6dSevRmbpGJ88v7RijgA2NTr7eUw7A0G4RptYiIiLiiRQAxTUK+O+v9rCntNrkai5e4cFKquoasQb60duuBtAiIiInUgAULkmK5NLeUTQ4Df5v6Razy7loa44tABnSNQJfH4u5xYiIiHggBUABcO0R/Naafew8VGVyNRdnzS41gBYRETkTBUABmhZLjE+OptFp8Jd2Pgro2gFEC0BEREROSQFQXJr3CH5n7T62lRwxuZoLc+hILbsOV2OxwJDECLPLERER8UgKgOIyODGCCSl2nAb85eP2OQrYfPu3tz0MW7C/ydWIiIh4JgVAaaF5LmDW1/vZcrDS5GrO3+rd6v8nIiJyNgqA0sKALjYm9o/BMOC5djgKuHZXOaAdQERERM5EAVBO0jwK+L91B9hUVGFyNeeursHJ13vLAY0AioiInIkCoJwkJS6cyYPiAHgup/2MAhYcqKC2wUlEiD89okLNLkdERMRjKQDKKc26sjcWC2RvKOKbfQ6zyzknq4/r/2exqAG0iIjI6SgAyin1jrFy3eB4AJ5b0j5GAV39/7pGmFuIiIiIh1MAlNN64Mre+FhgScFB1h2bW+fJXDuAaP6fiIjIGSkAymn1jA5jSmoXAJ7J2WxyNWd2wHGU/Y4afCwwOCHC7HJEREQ8mgKgnNEDV/TG18fCJ4Ulrjl2nmjNsfYvKXHhhAb6mVuMiIiIh1MAlDPqHhXKd4c2jQI+t8RzRwG/nf+n278iIiJnowAoZzXzit74+Vj4bMshVu0oNbucU2oenVT/PxERkbNTAJSzSowM4abhiQA864FzAWvqG9mwv6lVjUYARUREzk4BUM7JjCt64e9rYcX2wyzfdsjsclr4Zp+D+kaDqLBAEiODzS5HRETE4ykAyjnpEhHMrSO6Ak27gxiGYXJF3zq+/58aQIuIiJydAqCcsx9f3pMAPx9W7Szli62HzS7HRfP/REREzo8CoJyzOFswUy9pGgV8JqfQI0YBDcNgze5yQA2gRUREzpUCoJyXH1/ekyB/H9bsLid3c4nZ5bC37CgllbX4+1oY2MVmdjkiIiLtggKgnBe7NYiMUd2AphXBZo8CNs//6xdvI8jf19RaRERE2guPC4CZmZmMGDECq9WK3W5nypQpFBYWtjjGMAyeeOIJ4uPjCQ4OZvz48WzYsKHFMdu2beOGG24gOjqa8PBwbr75Zg4ePNjimLKyMjIyMrDZbNhsNjIyMigvL2/rS2z37r2sJ8H+vny918HSTcWm1uKa/6f2LyIiIufM4wJgbm4u06dPJy8vj5ycHBoaGkhPT6eqqsp1zNNPP80zzzzD888/z5dffklsbCxXXXUVlZWVAFRVVZGeno7FYmHp0qV88cUX1NXV8Z3vfAen0+k6z9SpU8nPzyc7O5vs7Gzy8/PJyMhw+zW3N1Fhgdw+ujvQtEewmaOArhXA3SJMq0FERKS9sRhm38M7i5KSEux2O7m5uYwbNw7DMIiPj2fWrFk88sgjANTW1hITE8Pvf/977r33Xj766COuvvpqysrKCA8PB5pG+yIjI8nJyWHChAkUFBTQr18/8vLyGDlyJAB5eXmkpaWxadMmkpOTz1pbRUUFNpsNh8Pheh9vUVpVx6W/X0pVXSMvZQxjYv9Yt9dQXdfAwCc+otFpsPznVxAfoR6AIiJydt78+d3M40YAT+RwNO3wEBkZCcCOHTsoKioiPT3ddUxgYCCXXXYZy5cvB5oCocViITAw0HVMUFAQPj4+fP755wCsWLECm83mCn8Ao0aNwmazuc4jpxcZGsCPxiQBTXMBnU73/x7x9R4HjU6DOFuQwp+IiMh58OgAaBgGs2fPZuzYsQwYMACAoqIiAGJiYlocGxMT4/raqFGjCA0N5ZFHHqG6upqqqip++tOf4nQ6OXDggOs8drv9pPe02+2u85yotraWioqKFg9vdvelSVgD/dhUVMkH35z676wtfdsAWvP/REREzodHB8AZM2awbt06Fi5ceNLXTtzxwTAM13PR0dH8+9//5r///S9hYWGuYd6hQ4fi6+t72nOceJ4TZWZmuhaM2Gw2EhMTL+by2r2IkADuHNs0Cvjcks00unkUcM2u5vl/CoAiIiLnw2MD4MyZM8nKymLZsmUkJCS4no+NbZprduIoXXFxcYtRwfT0dLZt20ZxcTGHDh3i9ddfZ9++fSQlJbnOc+KqYGiac3ji6GKzOXPm4HA4XI89e/Zc9HW2d3eOTSI8yI8txUd4b91+t71vUwPob7eAExERkXPncQHQMAxmzJjB4sWLWbp0qSuwNUtKSiI2NpacnBzXc3V1deTm5jJ69OiTzhcVFUVERARLly6luLiY6667DoC0tDQcDgerVq1yHbty5UocDscpzwNNcw3Dw8NbPLydLdifaZf2AODPS7bQ0Og8yytax45DVZRV1xPg50P/eDWAFhEROR8eFwCnT5/O/PnzWbBgAVarlaKiIoqKijh69CjQdNt21qxZPPnkk7z99tt888033HHHHYSEhDB16lTXeebNm0deXh7btm1j/vz53HTTTTz00EOu1b0pKSlMmjSJadOmkZeXR15eHtOmTePaa689pxXA8q07xnQnIsSf7YeqyPraPaOAzf3/BnWxEeDncf+MRUREPJqf2QWcaO7cuQCMHz++xfPz5s3jjjvuAOBnP/sZR48e5cc//jFlZWWMHDmSjz76CKvV6jq+sLCQOXPmUFpaSvfu3Xnsscd46KGHWpzzjTfe4IEHHnCtKL7uuut4/vnn2+7iOihrkD/3jOvB09mF/PnjLVw3OB4/37YNZc37/w7T/D8REZHz5vF9AD2Z+gh9q6q2gUufXkZpVR1Pf28QNw9v2wUyE5/9lMKDlbz4g2FMGuD+HoQiItJ+6fPbA28BS/sUGujHfZc1zQX8y8dbqG/DuYAVNfVsLm7a9UU7gIiIiJw/BUBpNRmjuhMVFsjesqP8Z/XeNnuf/N3lGAYkRgZjtwa12fuIiIh0VAqA0mqCA3z58fieADy/dCu1DY1t8j7N7V+GqQG0iIjIBVEAlFY1dWRXYsID2Vd+lH991TajgKvVAFpEROSiKABKqwry92X65b0A+OvSrdTUt+4ooNNpkH9sBbC2gBMREbkwCoDS6m4ZkUicLYiiihreXLW7Vc+9pfgIlbUNhAT40jfWevYXiIiIyEkUAKXVBfr5MuOKY6OAn2xr1VHA5vl/gxMi2rzXoIiISEelT1BpEzcNS6RLRDAllbXMz9vVauf9dv5fRKudU0RExNsoAEqbCPDz4YErm0YB536yjeq6hlY5r2sFsBaAiIiIXDAFQGkzNw5NoGtkCIer6nhtxcWPApZV1bG9pAqA1EQFQBERkQulAChtxt/Xhweu7A3AS7nbOFJ7caOAa/c0jf71iA6lU2jARdcnIiLirRQApU1NGRJPUlQoZdX1/HP5zos6l2v+n9q/iIiIXBQFQGlTfr4+PHhsFPDlT7dTUVN/wedas6sc0Pw/ERGRi6UAKG3uO4Pj6WUPw3G0nnmf77ygczQ0OsnfUw5oBFBERORiKQBKm/P1sTBrQtMo4N8+346j+vxHATcVVXK0vhFroB+97WGtXaKIiIhXUQAUt7hmQBzJMVYqaxr4++fbz/v1ze1fhnSNwMfH0trliYiIeBUFQHELHx8LD13VNAr46hc7KauqO6/Xr9ml/n8iIiKtRQFQ3Ca9Xyz94sI5UtvAK5+d3yjg6t1aASwiItJaFADFbZpGAfsA8I/lOzl8pPacXldcWcOe0qNYLE23gEVEROTiKACKW01IsTOwi43qukZe/vTcRgGb27/0sVsJD/Jvw+pERES8gwKguJXFYmH2sVHAf67YSUnl2UcB1zbf/tX8PxERkVahAChuNz45miGJEdTUO5n7ybazHv/tDiARbVyZiIiId1AAFLc7fhRw/spdHKyoOe2xdQ1O1u1zAFoBLCIi0loUAMUUl/aOYni3TtQ1OHlh2dbTHrdhv4O6BiedQvxJigp1Y4UiIiIdlwKgmOL4UcCFq/awv/zoKY9bs7scaGr/YrGoAbSIiEhrUAAU06T17MzIpEjqGp389TSjgM0NoLUAREREpPUoAIppjh8F/NdXe9hTWn3SMWvUAFpERKTVKQCKqUb26MzYXlHUNxonjQLuLz/KAUcNvj4WBifaTKpQRESk41EAFNM17xH879V72XW4yvV88+hfSpyVkAA/U2oTERHpiBQAxXTDukVyWZ9oGp0G/7f021HAb/v/6faviIhIa1IAFI/QvEfw4jV72XGoaRSweQWw+v+JiIi0LgVA8QhDEiO4sq8dpwF/+XgLNfWNbDjWAFojgCIiIq1LAVA8RvMo4Lv5+3h77T4anAbR1kASOgWbXJmIiEjHogAoHmNAFxvp/WJwGvDr/24Emvb/VQNoERGR1qUAKB5l1oSmUcCj9Y2A5v+JiIi0BY8LgJmZmYwYMQKr1YrdbmfKlCkUFha2OMYwDJ544gni4+MJDg5m/PjxbNiwocUxRUVFZGRkEBsbS2hoKEOHDuU///lPi2PKysrIyMjAZrNhs9nIyMigvLy8rS9RzqBffDjXDIx1/Vnz/0RERFqfxwXA3Nxcpk+fTl5eHjk5OTQ0NJCenk5V1bf94Z5++mmeeeYZnn/+eb788ktiY2O56qqrqKysdB2TkZFBYWEhWVlZrF+/nhtvvJFbbrmFtWvXuo6ZOnUq+fn5ZGdnk52dTX5+PhkZGW69XjnZrAl98LGANdCPAV3UAFpERKS1WQzDMMwu4kxKSkqw2+3k5uYybtw4DMMgPj6eWbNm8cgjjwBQW1tLTEwMv//977n33nsBCAsLY+7cuS0CXefOnXn66ae56667KCgooF+/fuTl5TFy5EgA8vLySEtLY9OmTSQnJ5+1toqKCmw2Gw6Hg/Dw8Da4eu+1elcZgX4+CoAiItLq9PntgSOAJ3I4mlqBREZGArBjxw6KiopIT093HRMYGMhll13G8uXLXc+NHTuWRYsWUVpaitPp5M0336S2tpbx48cDsGLFCmw2myv8AYwaNQqbzdbiPGKOYd06KfyJiIi0EY/eX8swDGbPns3YsWMZMGAA0DS3DyAmJqbFsTExMezatcv150WLFnHLLbfQuXNn/Pz8CAkJ4e2336Znz56u89jt9pPe0263u97jRLW1tdTW1rr+XFFRcXEXKCIiImICjx4BnDFjBuvWrWPhwoUnfe3E1iCGYbR47he/+AVlZWUsWbKEr776itmzZ3PTTTexfv36057jVOc5XmZmpmvBiM1mIzEx8UIvTURERMQ0HhsAZ86cSVZWFsuWLSMhIcH1fGxs0wrRE0fpiouLXaOC27Zt4/nnn+fVV1/lyiuvZPDgwTz++OMMHz6cv/71r67zHDx48KT3LSkpOWl0sdmcOXNwOByux549e1rlWkVERETcyeMCoGEYzJgxg8WLF7N06VKSkpJafD0pKYnY2FhycnJcz9XV1ZGbm8vo0aMBqK6uBsDHp+Xl+fr64nQ6AUhLS8PhcLBq1SrX11euXInD4XCd50SBgYGEh4e3eIiIiIi0Nx43B3D69OksWLCAd999F6vV6hrps9lsBAcHY7FYmDVrFk8++SS9e/emd+/ePPnkk4SEhDB16lQA+vbtS69evbj33nv54x//SOfOnXnnnXfIycnhvffeAyAlJYVJkyYxbdo0XnrpJQDuuecerr322nNaASwiIiLSXnlcG5jTzb+bN28ed9xxB9A0SvirX/2Kl156ibKyMkaOHMlf//pX10IRgC1btvDzn/+czz//nCNHjtCrVy8efvjhFm1hSktLeeCBB8jKygLguuuu4/nnnyciIuKcatUychERkfZHn98eGADbE/0DEhERaX/0+e2BcwBFREREpG0pAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl/G4RtDtSXMHnYqKCpMrERERkXPV/LntzZ3wFAAvQmVlJQCJiYkmVyIiIiLnq7KyEpvNZnYZplAj6IvgdDrZv38/Vqv1tDuYXKiKigoSExPZs2ePVzap1PV79/WD/g68/fpBfwe6/ra7fsMwqKysJD4+Hh8f75wNpxHAi+Dj40NCQkKbvkd4eLhXfuM30/V79/WD/g68/fpBfwe6/ra5fm8d+WvmnbFXRERExIspAIqIiIh4GQVADxUYGMjjjz9OYGCg2aWYQtfv3dcP+jvw9usH/R3o+r37+tuaFoGIiIiIeBmNAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqAHeuGFF0hKSiIoKIhhw4bx2WefmV2S23z66ad85zvfIT4+HovFwjvvvGN2SW6VmZnJiBEjsFqt2O12pkyZQmFhodlluc3cuXMZNGiQq/FrWloaH3zwgdllmSYzMxOLxcKsWbPMLsVtnnjiCSwWS4tHbGys2WW53b59+/jBD35A586dCQkJYciQIaxevdrsstyie/fuJ/0bsFgsTJ8+3ezSOhQFQA+zaNEiZs2axWOPPcbatWu59NJLufrqq9m9e7fZpblFVVUVgwcP5vnnnze7FFPk5uYyffp08vLyyMnJoaGhgfT0dKqqqswuzS0SEhJ46qmn+Oqrr/jqq6+44ooruP7669mwYYPZpbndl19+ycsvv8ygQYPMLsXt+vfvz4EDB1yP9evXm12SW5WVlTFmzBj8/f354IMP2LhxI3/605+IiIgwuzS3+PLLL1v8/8/JyQHgpptuMrmyjkVtYDzMyJEjGTp0KHPnznU9l5KSwpQpU8jMzDSxMvezWCy8/fbbTJkyxexSTFNSUoLdbic3N5dx48aZXY4pIiMj+cMf/sBdd91ldiluc+TIEYYOHcoLL7zAb3/7W4YMGcJzzz1ndllu8cQTT/DOO++Qn59vdimm+fnPf84XX3zhVXd/zmTWrFm89957bNmyBYvFYnY5HYZGAD1IXV0dq1evJj09vcXz6enpLF++3KSqxEwOhwNoCkHeprGxkTfffJOqqirS0tLMLsetpk+fzuTJk5kwYYLZpZhiy5YtxMfHk5SUxK233sr27dvNLsmtsrKyGD58ODfddBN2u53U1FReeeUVs8syRV1dHfPnz+fOO+9U+GtlCoAe5NChQzQ2NhITE9Pi+ZiYGIqKikyqSsxiGAazZ89m7NixDBgwwOxy3Gb9+vWEhYURGBjIfffdx9tvv02/fv3MLstt3nzzTdasWeN1I/7NRo4cyWuvvcaHH37IK6+8QlFREaNHj+bw4cNml+Y227dvZ+7cufTu3ZsPP/yQ++67jwceeIDXXnvN7NLc7p133qG8vJw77rjD7FI6HD+zC5CTnfhbjmEY+s3HC82YMYN169bx+eefm12KWyUnJ5Ofn095eTlvvfUWt99+O7m5uV4RAvfs2cODDz7IRx99RFBQkNnlmOLqq692/ffAgQNJS0ujZ8+e/POf/2T27NkmVuY+TqeT4cOH8+STTwKQmprKhg0bmDt3Lj/84Q9Nrs69/v73v3P11VcTHx9vdikdjkYAPUhUVBS+vr4njfYVFxefNCooHdvMmTPJyspi2bJlJCQkmF2OWwUEBNCrVy+GDx9OZmYmgwcP5s9//rPZZbnF6tWrKS4uZtiwYfj5+eHn50dubi5/+ctf8PPzo7Gx0ewS3S40NJSBAweyZcsWs0txm7i4uJN+4UlJSfGaxYDNdu3axZIlS7j77rvNLqVDUgD0IAEBAQwbNsy14qlZTk4Oo0ePNqkqcSfDMJgxYwaLFy9m6dKlJCUlmV2S6QzDoLa21uwy3OLKK69k/fr15Ofnux7Dhw/ntttuIz8/H19fX7NLdLva2loKCgqIi4szuxS3GTNmzEntnzZv3ky3bt1Mqsgc8+bNw263M3nyZLNL6ZB0C9jDzJ49m4yMDIYPH05aWhovv/wyu3fv5r777jO7NLc4cuQIW7dudf15x44d5OfnExkZSdeuXU2szD2mT5/OggULePfdd7Fara7RYJvNRnBwsMnVtb1HH32Uq6++msTERCorK3nzzTf55JNPyM7ONrs0t7BarSfN9wwNDaVz585eMw/04Ycf5jvf+Q5du3aluLiY3/72t1RUVHD77bebXZrbPPTQQ4wePZonn3ySm2++mVWrVvHyyy/z8ssvm12a2zidTubNm8ftt9+On5+iSpswxOP89a9/Nbp162YEBAQYQ4cONXJzc80uyW2WLVtmACc9br/9drNLc4tTXTtgzJs3z+zS3OLOO+90/duPjo42rrzySuOjjz4yuyxTXXbZZcaDDz5odhluc8sttxhxcXGGv7+/ER8fb9x4443Ghg0bzC7L7f773/8aAwYMMAIDA42+ffsaL7/8stkludWHH35oAEZhYaHZpXRY6gMoIiIi4mU0B1BERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TIKgCIiIiJeRgFQRERExMsoAIqIiIh4GQVAERERES+jACgiIiLiZRQARURERLyMAqCIiIiIl1EAFBEREfEyCoAiIiIiXkYBUERERMTLKACKiIiIeBkFQBEREREvowAoIiIi4mUUAEVERES8jAKgiIiIiJdRABQRERHxMgqAIiIiIl5GAVBERETEyygAioiIiHgZBUARERERL6MAKCIiIuJlFABFREREvIwCoIiIiIiXUQAUERER8TL/H+lHdNEfpEu4AAAAAElFTkSuQmCC' width=640.0/>\n",
|
|
" </div>\n",
|
|
" "
|
|
],
|
|
"text/plain": [
|
|
"Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …"
|
|
]
|
|
},
|
|
"metadata": {},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"fig, ax = plt.subplots()\n",
|
|
"ax.plot(res['x'], res['y'])"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 48,
|
|
"id": "6c1df388-6e5f-4fd6-8282-b319a43af001",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"scan0000_test_bs_scan_info.json scan0005_test_bs_scan_info.json\n",
|
|
"scan0001_test_bs_scan_info.json scan0006_test_bs_scan_info.json\n",
|
|
"scan0002_test_bs_scan_info.json scan0007_test_bs_scan_info.json\n",
|
|
"scan0003_test_bs_scan_info.json scan0008_test_bs_scan_info.json\n",
|
|
"scan0004_test_bs_scan_info.json\n"
|
|
]
|
|
}
|
|
],
|
|
"source": [
|
|
"!ls /tmp/scan_info/\n"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 47,
|
|
"id": "3dcac2e6-84ee-4f6d-95f2-20ce1f94d593",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"name": "stdout",
|
|
"output_type": "stream",
|
|
"text": [
|
|
"scan0000_test_bs scan0003_test_bs scan0006_test_bs scan0008_test_bs\n",
|
|
"scan0001_test_bs scan0004_test_bs scan0007_test_bs scan0008_test_pv\n",
|
|
"scan0002_test_bs scan0005_test_bs scan0007_test_pv\n"
|
|
]
|
|
}
|
|
],
|
|
"source": [
|
|
"!ls /tmp/scan_data/"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 39,
|
|
"id": "8e2a3f25-106a-4684-a6f6-441948a09e31",
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"ename": "FileNotFoundError",
|
|
"evalue": "[Errno 2] No such file or directory: '/tmp/scan_info/scan0007_test_pv_scan_info.json'",
|
|
"output_type": "error",
|
|
"traceback": [
|
|
"\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
|
|
"\u001b[0;31mFileNotFoundError\u001b[0m Traceback (most recent call last)",
|
|
"Cell \u001b[0;32mIn[39], line 1\u001b[0m\n\u001b[0;32m----> 1\u001b[0m scan \u001b[38;5;241m=\u001b[39m \u001b[43mSFScanInfo\u001b[49m\u001b[43m(\u001b[49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[38;5;124;43m/tmp/scan_info/scan0007_test_pv_scan_info.json\u001b[39;49m\u001b[38;5;124;43m\"\u001b[39;49m\u001b[43m)\u001b[49m\n\u001b[1;32m 3\u001b[0m res \u001b[38;5;241m=\u001b[39m defaultdict(\u001b[38;5;28mlist\u001b[39m)\n\u001b[1;32m 4\u001b[0m res[\u001b[38;5;124m'\u001b[39m\u001b[38;5;124mx\u001b[39m\u001b[38;5;124m'\u001b[39m] \u001b[38;5;241m=\u001b[39m scan\u001b[38;5;241m.\u001b[39mreadbacks\n",
|
|
"File \u001b[0;32m/sf/cristallina/applications/conda/envs/analysis_edge/lib/python3.11/site-packages/sfdata-0.4.0-py3.11.egg/sfdata/sfscaninfo.py:13\u001b[0m, in \u001b[0;36mSFScanInfo.__init__\u001b[0;34m(self, fname)\u001b[0m\n\u001b[1;32m 11\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mfname \u001b[38;5;241m=\u001b[39m fname\n\u001b[1;32m 12\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mfs \u001b[38;5;241m=\u001b[39m FileStatus(fname)\n\u001b[0;32m---> 13\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39minfo \u001b[38;5;241m=\u001b[39m info \u001b[38;5;241m=\u001b[39m \u001b[43mjson_load\u001b[49m\u001b[43m(\u001b[49m\u001b[43mfname\u001b[49m\u001b[43m)\u001b[49m\n\u001b[1;32m 15\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mfiles \u001b[38;5;241m=\u001b[39m info[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mscan_files\u001b[39m\u001b[38;5;124m\"\u001b[39m]\n\u001b[1;32m 16\u001b[0m \u001b[38;5;28mself\u001b[39m\u001b[38;5;241m.\u001b[39mparameters \u001b[38;5;241m=\u001b[39m info[\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mscan_parameters\u001b[39m\u001b[38;5;124m\"\u001b[39m]\n",
|
|
"File \u001b[0;32m/sf/cristallina/applications/conda/envs/analysis_edge/lib/python3.11/site-packages/sfdata-0.4.0-py3.11.egg/sfdata/utils/json.py:5\u001b[0m, in \u001b[0;36mjson_load\u001b[0;34m(filename, *args, **kwargs)\u001b[0m\n\u001b[1;32m 4\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21mjson_load\u001b[39m(filename, \u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs):\n\u001b[0;32m----> 5\u001b[0m \u001b[38;5;28;01mwith\u001b[39;00m \u001b[38;5;28mopen\u001b[39m(filename, \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mr\u001b[39m\u001b[38;5;124m\"\u001b[39m) \u001b[38;5;28;01mas\u001b[39;00m f:\n\u001b[1;32m 6\u001b[0m \u001b[38;5;28;01mreturn\u001b[39;00m json\u001b[38;5;241m.\u001b[39mload(f, \u001b[38;5;241m*\u001b[39margs, \u001b[38;5;241m*\u001b[39m\u001b[38;5;241m*\u001b[39mkwargs)\n",
|
|
"\u001b[0;31mFileNotFoundError\u001b[0m: [Errno 2] No such file or directory: '/tmp/scan_info/scan0007_test_pv_scan_info.json'"
|
|
]
|
|
}
|
|
],
|
|
"source": [
|
|
"scan = SFScanInfo(\"/tmp/scan_info/scan0007_test_pv_scan_info.json\")\n",
|
|
"\n",
|
|
"res = defaultdict(list)\n",
|
|
"res['x'] = scan.readbacks\n",
|
|
"\n",
|
|
"for step in scan:\n",
|
|
" intensity = step[\"SARFE10-PBIG050-EVR0:CALCI\"]\n",
|
|
" res['y'].append(np.mean(intensity.data))"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"id": "03eb2976-ae1f-43b9-8169-c9bf1eb73550",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"id": "980f9e86-c2d0-49fd-802f-d04d9791bae6",
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
}
|
|
],
|
|
"metadata": {
|
|
"kernelspec": {
|
|
"display_name": "Python [conda env:conda-analysis_edge] *",
|
|
"language": "python",
|
|
"name": "conda-env-conda-analysis_edge-py"
|
|
},
|
|
"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
|
|
}
|