{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Generating a sky map and finding sources" ] }, { "cell_type": "raw", "metadata": {}, "source": [ "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**In this tutorial you will learn how to generate a sky map from the selected events and how to find potential sources in the data.**\n", "\n", "As next step you will make a sky map of the selected events to inspect what’s in the data. You do this using the ctskymap tool. \n", "\n", "We start with the usual gammalib, ctools, and cscripts imports." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import gammalib\n", "import ctools\n", "import cscripts" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We will use the matplotlib package to display the results." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We will use the events selected in the previous step. Since the data correspond to observations of the Crab nebula we should centre the sky map on this source. The energy limits of the sky map should be set to the energy interval of the data (0.66 - 100 TeV). " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "obsfile = 'obs_crab_selected.xml'\n", "emin = 0.66\n", "emax = 100.0\n", "skymap = ctools.ctskymap()\n", "skymap['inobs'] = obsfile\n", "skymap['proj'] = 'CAR'\n", "skymap['coordsys'] = 'CEL'\n", "skymap['xref'] = 83.63\n", "skymap['yref'] = 22.01\n", "skymap['binsz'] = 0.02\n", "skymap['nxpix'] = 200\n", "skymap['nypix'] = 200\n", "skymap['emin'] = emin\n", "skymap['emax'] = emax\n", "skymap['bkgsubtract'] = 'NONE'\n", "skymap.run()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can visualise the skymap." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUYAAAEKCAYAAABuTfznAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztvXn8JUV1Nv6cWdgGZGcWQAZlB1mvimAUBzVEQY2Ixi1GURwhBoyiCEQuKsZIQn5GE8lEEY1IIou+OoJgWAReFb0ssjgas4g/wsC44AYKwTnvH933Urdu1alTS9/lO/18Pvdz7+0+dep0d/VTT1VXVxEzo0WLFi1aPI55kw6gRYsWLaYNLTG2aNGihYWWGFu0aNHCQkuMLVq0aGGhJcYWLVq0sNASY4sWLVpYaIwYiWhnIrqOiNYQ0d1EdLKx761E9P16+4eaiqFFixazCSK6gIjWEdFd1vaxcMeCphwDeAzA25n5ViLaAsAtRPRVAIsBvBjA/sz8CBHt0GAMLVq0mE1cCOCjAD7d30BEz8GYuKMxYmTmtQDW1r9/RURrAOwI4E0APsjMj9T71jUVQ4sWLWYTzHwDES23Nr8FY+KOJhXjAPUBHgTgZgDnAvg9IjoHwG8BvIOZvy2n34yBrfKCWLqs+l57X56fEli6TBeH1s60B5o/Rjuu2Dg1Pl3Yfwlwx/15+YwTrmMKHWfONcy6Dj8H88OUmBgAsBsRP6y0XQvcjer+72MVM68KJNsDkdyRisaJkYg2B3AZgFOY+ZdEtADA1gAOBfBUAJ8joiex9W4iEZ0A4ITq35aP/9Sg260+5u+1OUcRyMfMT4qlD00srnTmMbnQxDHa0ByPdD586PuQ0t4R4c/nJ3QOc2Hmu9aRn+8alSinWdc/xElhPAzgzUrbLvBbZu5EZqHijhKgJt+VJqKFAFYDuIqZz6u3fQWVHL6+/v+fAA5l5h/7/SxjNTGO64aIvflTyKIJaMh81uC6tvbxaY435pxMy3ksUa67XaC7DMz3ZSnGHYn4RKXtmcAtIWKsW5qrmXm/+n80d6SiMWIkIgLwKQA/Y+ZTjO0rASxj5vcQ0R4ArgHwRIn1RWIsWTCbKuS133n3n4r15y+aXkIah5oq7X/SxGSiRCzFiC42/apsYtyZiN+mtH17GjFGc0cqmhzHeDiA1wJYQUS3158XALgAwJPqx/D/AuB1TRxYtJqLTROLbhfrl5yrUy2TQqxi0mz3qbZCxzlv5UN646Yr0Vi1GbtPiwmVIQKwUPkJ+iK6GMA3AOxJRPcS0fEYF3eg4aZ0KUQ1pU3YN+U0wEcMTSifcRC+nV/T3QuTVojSOW06ttKEOpQmXzHuQsSnKW1PVCjGSaLxhy8tWrTYMNBXjHMBc/uVwFKd4zGd+KEnyXZz0vg97/5TnWnn3X9qXJ6pscbuT7UPqa4SeaRA0zQuqdpKnt+U81L4XBIqpaX5TDtmgxj7Y7tcKHhxvQRkE1kqMdrpLV/rl5zrdBnsm0ztMgg1tUs/JPHl6zkfSb5jYzD/5/QZ2r+lY/Idv+Q3B+O4vijbxzhpzL0+xlyFGEN8OTe0Jk2M30n3vZmQ+lFd22P9aPMKKP2+QndWSEIFVhw5FUKxGPP7GHcjYu3Ly8dOeR/jbChGE6WaWjnNT5fyC/m3b0rN09QUkvb91/iIhU8N1Z+hrgGf2s4lHN858ik3S6H7VLqTcLTKPHZfaplOaSk0SPJzSTHOPWLU2GhvJp9P303jK6iOm2v9+YvkGDUIHYfWR27+nnMwdIzSuQn5z4GkGrWqXfodq1QDec67/1RnhTLS/5wKVwVVkCgXKD8h+GbXqfe9g4iYiLYrFridx8w1pUvdLD4Fk+Kj/9+1P5S+xrz7T/WrlxZxiL2uMfZS5ZmalzZ/bQtC639oe35Teg8i/gel7fMCTWkiehaAXwP4dH+Ad719ZwAfB7AXgEOY+Sc5MfswNxWj1oemieLb7muqafOwMESKJZqbDTWXZgKxx+4jK61tTBw5TemQykupEPqQHnAqUfKpNDPfAOBnjl1/C+CdABpVdLNHjDY0BS31RgmRnU2KgXxGmkhS3vY2LXIVTaqvcZC3Ri2lokRlZKcvpVq1xJgTR4EZmeYB2FT5AbAdEfWMT/DpKhG9CMD/MPN3soMN5TVzTenS8PQBiU8sI/3PW/nQqB9HniPbE/KaiFKcVL4hpMQ1rcfSOPKb0vsQ8WeUtodEvitNRJsBuA7A85n5F0T0QwCdDbsprZX5oeatJm0N5xPLRPISh4OkxOqLIUZxlMS0kUmoS6NJ5dn3EdNv2AS0zfOCaHiA95MB7ArgOzUp7gTgViJakhm2E7NBjC6ZLzU5fTdDiExChVm6yUL+Q+QUS3SWshXtY2MJISZW7f6ScF2T2GaoMl7n+be7YlKIOLeymUBF1eRwHWa+k5l3YOblzLwcwL0ADmbmRmYung1ibNGixdSjpGL0zK4zNrR9jONEbAe7gKHhPanqonQT2KfiQ4ppmprhwEhMqqFUvq4Rn4osGav5nYz8Psb9ifiLSttdp/zNl7lPjHZBTCmYgWay8+FKir+Y9K7fIXv7Zu8PwPbdsJMmrRI3fYljkM5NSvfLtGAo9nxiPICIr1baLplyYpz7TemU/h2fneu7280nxcJqzzlLj2U/MjGF5Wve/afGTQLri9HhOyp9CVJMPb+u3/a2UL9tycqlNLkWJm0CsHCB7jPtmPvEaCNWbbnSdgOvaMWScWH1qJoUwU5vxTB4Ki/dPCXUkuLYV2iEhURIIfJKjMvrW+tfUxYtG7vcjbxCGMqzQRABCxboPtOOud+UHjekG2PSTVMtNITXhHpxfSYNXxzK+AbjYe21fgLpo/o1iyC/KX3wfOL/u6nOdrOH2qZ0WYQKQkThCxJYirLIUVCl1aTtR6uYfEpGSK9S0Jpzmtq8K0mipeLtdv3KO9D0VnXPaMpe4eayhLmkGBsjRiLamYiuI6I1RHQ3EZ1cb+8S0f9YC2TpESKVUPOk3ibOip2iWJQkHN0ElxRLDDQVhmu7klDFGzm36Ro6/ibUa8ivkuC86WCUwRTy0paJMRPjwo11n2lHk8unLgWwlJlvJaItANwC4CUAXg7g18z812pfyzqMN/fkmly7Tdov2U+qaRdbsEucB41NifMR4yPn2rjIR1spas/LOMtHaqUtVib5TenORsS97XW2dF9wdp0LABwNYJ2xfOq5AI4B8CiA/wTwemb+eU7M3vzH1cdIRP8HwEdRLasaR4yhPkapcDZBChrb2FgaUqgjMdUYGq5j26cghRiaVDMliSrVl02wAZUu9iuWrohGzn0BYtyYuKd8QY9+FD/tGBE9H8C1zPwYEf0VADDzu3Ji9uY/DmKsXwa/AcB+AP4cwJ8A+CWAHoC3M/ODjjQnYMCGWx4CnBKXaaAmHymE2oI3boWgza+JeEJE7vod4zci5gGR5xJU0/CQ/Qru4FrqjdqmlKXYClZVARUixp10tvRfcZNIOPb9IYCXMfOr4yMNo/GHL0S0OYDLAJzCzL8E8DFUL4QfCGAtgL9xpWPmVczcqU7eZvEZW4XBHtYwMoN2qtoJ5BvtO6lQK/Nw2DiHwpRoMmtuXk+F5ToHoXGXI/5d/3MVqqW4TQz+e1oJ11LVFTQ43y71FhmDarudV5MgAPOVn4Rpxyy8AcCVhSIfwYKmHAMAES1ERYoXMfPlAMDMDxj7/wnA6iZjaNGixZjQf1lah5+kDtchojMAPAbgopT0qjwafPhCAD4F4GfMfIqxfSkzr61/vw3A05n5j2RfkeMYm+y7agI5NXpuWvN7XJDy1fbRapqh42xCazDpeMT8CzSlNyPu7amzpdvTmtJE9DoAKwEcycwPJwcbQJNN6cMBvBbACmtozoeI6E4iugPAcwC8rZHcXc3BmGZITv9abrPc5T+mie1rorrsrOah6i0K7fHZ51DTfJSIz5W+RBeDy7aJ/lqNja/chbZJZSgmhhw0PCEjER0F4F0AXtQkKQJz4c0XbS3suqE0SkO6uQti8DAo1n/T9r500jnKQWkVGyLY0shRsL5zLFWOMf5EFFCMmxP3nqKzpW8Gn0pfDOAIANsBeADAWQDeDWBjAD+tzb7JzCszQvZi9t98yVEv9m9NMy2VVHx51TCnECuq2mr0fQ4mhkhR0KmxxF6j0qTo+h3jw/U71kaTr9QqCMU/CYXogv7hiwhmfiUzL2Xmhcy8EzN/gpl3Y+admfnA+tMIKQJzgRib8p1Kitomj+BL/TpYxA3SV6MjU42VQEliLGFj2oYqw5BfbQsj5CeULsZeU4Hn5JOKhpvS48TsEWMT0N7YMU2xJkknoGid047ZviRiiCWoEMnEnKuUyimktEJ5SjFI2zR+Y+xSyC21LDZBlnOIGGe/jzGEuuCkznKizsP+7VIrKfloVYGUR+4xappysXmlxppzHqX/peEj9ZL5FvVXoI9xS+LeM3W2dEU7u04+7FUCY2vlrmOWkxJwNWnM/y7bnDy0tsL2pLn7tDG4KgPf7xj17YtLE5tW7ZcmylJ5xKrImHJX+pjnkGKcDWLsrxIYQxJ9O40C0fgJbZfIcFxNma7nwU2d19CTb4Uv57ZUEtOQpJRvKA+Hz2AloMl7kkg53y0xFsHcbErHFKhStj71GPJXY97Kh4YfjLi+fXH4/lvpR9Z5keKKiT9UQdj2rrQBRaeauDUHISUZuiaSX9fvWPjKV6jyVedZoCm9DXHvuTpbumS6m9JzkxhdiCnQmgJlF06FX+fEFa4YtQjFEFJfMQQYcz6kPD1xjax66LFz5uk7dheZjEsV5p6zHNskFCDGbYl7v6+zpYunmxhnQNS2aNFiJkCohl/PAcxGH2MJ+PqTSjVvFHZetZgSi1YdhhRlSFXWCK48GNs0t/5HTwEXaqJqY4vJL7GbJGijuQYprQlN2v4++wFnCgr2MRLRBUS0jojuMrZtQ0RfJaIf1N9b5wftxtwkxpibVdPclPKIjGNoOnsfWZvx+fKS0sc0qV1pHNtG+vdC5y+miajAEDFrCUo6v5Kdy15LjNL19Nm5YkmtsJXXcwT9B5w5KPvw5UIAR1nbTgNwDTPvDuCa+n8j2HD6GCMxGPfoUnkhokkhDI+qmrfyofiHDrH5mPuk75B9LFLjDNlqrkUoL2tb1sTG5rfWxnVMqedZhQJ9jDsQ916hs6WPxs+uQ0TfB3AEM6+tl065npmV8/nEYW4qRiBNPRkYjHu07Xz+bIXnUxuuPAVyEGet7nZHJ3bNVRuauExbH6HEkEYoDsk+ZKshX4VaHHowJJWDUP6hY7D9SwowhNzrEIvmJ6pd3J+ysP7eoUzgo9gwFGMTtW2CYvAOObGVYUq8mpvLZ+/bF1AtQ6o6ReFo00iVjm1nfmvz9e3X5iVVwrHnXrINVcyh9CIKKMYlxL3X6Gzpb5IU48+ZeStj/4PM3Eg/4+wQY9ca5D1pxBZuM42QTj1erxTZS6pPk9b1u8bIa5ipRB0bl9aHVpWWiiWGKB0xjgxp0pC7en8BYlxK3HuDzpY+MN1N6dkhxpBi1CggrX0IqWl9N2WoALvSpqiO2BvTShOcM9KVj/G/P+VZ9KJWTVw/89vlV3t9pDwkwtPGIW3TxDhOxbiMuPcmnS29N4kYzwXwU2b+IBGdBmAbZn5nTszevGeOGGOVRYqyc/3W5GmnC8Wm2V7DuZDXOIm963gQ5CN567/6AZLvxi+hjFORE1OTqt6C1NLQtUIKEOOOxD3lDIn0nqSJar8A4HMAngjgRwCOY+af5cTszX/miLEPzU0poemmUCosXyP9eKZNV+i3dPk1v33bhFh8PlXE51O8vgqsxLnWVpCaeCw4l3PNPZ8hlCxnIyhAjDsR907S2dLp0/3mS5OLYe0M4NMAlgBYD2AVM3/Y2P8OAOcC2J6ZfyL7qolRuolSUMqHppln33TGt5PcNDey5N+RNut9Y+1NqW0CllJShr+iTXV7n484Q+c5VIGXVMVZwqAQMf6ZzpbeteES41IAS5n5ViLaAsAtAF7CzN+tSfPjAPYCcIiaGCU0WZv6bnbfbw1RuvKQ/sfElmJj55mj1kIk7SIb7flNiScAabSAGItWjZo2JctoboUzZF+AGHcm7imXtqO3TzcxNjaOkZnXMvOt9e9fAVgDYMd6998CeCeA6W/Ht2jRQof+u9Kaz5RjLH2M9dOlGwDsh6pD9UhmPpmIfgig41KM9YDPWiZueQhQL01dstaN6FNKsi+F2Kaxr3kWozAMH+qn0TnIaAZmD3GSlKvG3zi6DwIx5k/LVkAx7kLcO11nSys3UMXYBxFtDuAyVMz2GIAzALwnlI6ZVzFzpzp5mzUcJeJvTPM7Jk2Ez8GTaE8zc6gvzXUTSh/TzhdP/RnpL3PlA4y+z+yJZ4XrfugKE8sGmrTByXdD16vreIPIgGZN7mBeKWXGlZcnvYsUvdejKbQT1SqdEy0EsBrAVcx8HhE9BdXL3/3FsncCcB+ApzHz/V4/yzqMtUfHZT7Ovh3XvtQ8XTeB5UdSByu4g2vPPlq8yZ3DfSTVk3lM3odL9u/S5zVF4frOSWb5sVW3d7C2Nj8F2cedrwKKcVfi3lk6W3r9BqoYiYgAfALAGmY+DwCY+U5m3oGZlzPzcgD3AjhYIsUkuJRRaf+avGMVqMf3QGGZCs6lcrpdXEs9N9H0VabUJLZjcFQuQ7MDCTE78zS391Vmf53rEGKvpcYuRJzSdfQpQs82mwTFLhDJf6gSiSXWkmgVo8Ix0TMB3AjgTlTDdQDgdGa+wrD5ITx9jEO+lnUYb+65L3quKmsyrUYBSfty4oxFbl4+xemzteAcF6jJQ0qbckxaQvQpy1Q0ea1VirSAYnwyce+vdLZ03AaqGJn5JmYmZt6fmQ+sP1dYNstDpDiAqyBqVI9ru1YFxPi14O2zcik9028smRr/1SsAuuKx/WvVi4bIfGTjUlTKmE0MqWBtpRkibc/+eSsfqj5W/928+08d7RP2kPjIvqYrwHFVrnGz60w1ZmPasZRJNK3CsII7VZNUKvjaG0m6AevfzuZj/T3o47P9+/xKRG6SS+hVQY8P58MLXxx2rPW2oSZxSD36jssm6UDcI/99JBwgO5VtbbP+/EXVx3oTaWiaOqEcOUk8hbhCxy+l01ZssSjclCaitxHR3UR0FxFdTESblAs2kPfMvBLYvS98EZuueV151R/nxLamrfnt86XJS/KTghiVaGwzH+IkTaYr+E6y12zTnEfXcfr6ZbvGe+C+ysNln/uuu4TYsjKwLdCU3p2493c6W3pB8F3pHQHcBGAfZv4NEX0OwBXMfGFOjFrMhmJcukx3oWOJI6dgWnnZE9s610hJzcenBjS/7XgkJde3XfmQV0ENdQP0VZQweYEEb0xahJSP73xpylCNEOEPkZwvL0tZihWBppJ0xKne77MvgfJN6QUANiWiBajG7BVYf0GH2VGM9bvSgxpXoxgk+JoU2punJEG70rp8ScdsqheXL5dy8uWdGreGnOztWvWn3efyb/wXB8Vr1FUovxjEpNHa+sp1aF8Jxbgnce9jOls6EvcAMJ8vrGLmVUM2RCcDOAfAbwBczcyvzokvBrNFjGhgIlftDREgpWz1GXuTufJv4kaT8jS3W7DH6QVnKDe2iU+oU8+1I/bBpBOewdEjA+hjyDqG8GO3a2KIPk8FiHEv4t4ndLb0zGBTemtUL4a8AsDPAVwC4FJm/kxOjFrMRlPagLovK3DjAnCPyZN8SQpMG4Ok1Pq/XYRn/3eRlOV76EmplKcvLk0aTwy2ch3qZrAf1Fg+vU+ZY+KUrk+9v/8gxWXnGpAuvtnjikGKyWWrOL/BZWylmHzIqdRNlH348lwA/83MP2bm/wVwOYDDygQahi7EFi1atAihT4xl8CMAhxLRZqia0kcC6BXzHsDMKcYRpCidGtLTwRG11e3639LQNGO1TSWXeoxUiqZqcy756Tpex3Ab8/id6zq7/kvH2Vdqvn5Q+3/Ir+Bj5Fr51HvIr2Eb/bZKIEbxv73dpcTNGDWtFt/+KVSMzHwzgEsB3IrqJZF5AFaJiQpibhBj4OaJagIJ+YgPNnwkJ93omkLZb0KFbnTNjS/kNdKk7HaHKg5xogYNgSmPU0WEUlqTQGIrI9d58x2btM2Xj2Sv+S9Be6zatKko+FSamc9i5r2YeT9mfi0zP1IuUBkz8/Bl3v0/iJ+ZORbam7JELRtJrM4F3+3f1rd3WYSIGKOXdZWI37XP4XNodUFNzJoKSAvpOmiP33euNequRLkKVdBOFHj4sh9x71KdLe29gb4SWBp5c81ZCBVcRD640OzTEoonHmeT2E5n3QCDpnQKKRo+nPnYcXpiMNF/nW7IzmE7FLdPSWqUcAzRSOdIc/58qtM+ztBvjaK194XORZNiwkQ7Ue14MdalDRSkGbTNiUtSO5obQMhLPYbP919zg9t+E2N1prfSOIfU2PmlnH8zncaXltRCsXj2BycKVvqRr1kBxbg/ce/LOlt64nQrxtkhxm496D2pmZCJlBu/Rs64y2BazQ2ggYv0AvDFFj2AWrvdV0nEHr+WkEOkr4lLqtRCRBr6L9klVcQFiPEA4t5XdLa0rCXGbGT1MfoKn/Ymk/xo8s0h7pB6dBFpU+/hxhJRAjkN9S2mxB9SjCWuSSpiz4dLHad0J4WIuiQxHkjc+6rOlnZoiTEbqqa0C5paVKNMNIQpKQJX4YwlE03s9j77d5OEoFVi5rdvm+Q7QtWq/Ln+p+QbuLb9Sss76URMvMa2kclLfOUwiDLE+O1rdbbztm2JMRtD60o3oYRSmiWC+shfmMjhRyrwIZJ12TcQbzB/1z47jSO+4H9Nek3lIClKSY1q/U8C2oq9ADEecjDx12/S2W6yaLqJcWaeSgPQqSOfXajAu37X/72LIXnySm7OWvbBgdCS6ghtM9KuX3Ju3JIFIb8OP74FsEQyctiFXif05e+0t9LaY0XFQe0uv6Hz5aqomiZRDTEWAhPw2Px5qs+0Y3YUY9eacUggsWBTJfHGMfd7FzOy/WlVjL3PF4O936MKk1SgT2HW27zzLvpIw6WmXPm48jfgfCrrSqu91iECbgoR/rNVfEyZ7naB7rJsxXhQZx5/7eaFKtstFzy6YSpGItqZiK4jojX1LLwn19vfR0R3ENHtRHQ1ES1rKoYWLVqMF7+bP1/1mXY0uRjWUgBLmflWItoCwC0AXgLgXmb+ZW3zZ6hm6F0p+5L7GJ1vhWhqZpcCEdI6+/x8abQKJhSPz85OE1JF0n9t/JpmmVZxaY9V48f89uXju9Y+f5KvlP2px6hRfkWQ38d4YGc+/1tPtwb89vTrDVMxMvNaZr61/v0rAGsA7NgnxRqLAOiZ2VM4nGQVSOPd5yFKcX4+RXNQTGPv00DyVWNkgSZFmpF9MURQ/1Yt9O7JQ/W2UQzB+SoOl30ssYX+a/dJiCFwzbluEAzCY5iv+mhARFsR0aVE9L265fmMhg9hgLH0ghLRcgAHAbi5/n8OEf3/AF4N4D2eNCcQUY+IesDDOmXg+h+jYITtwT4ujW/r5vS+dhhzEwo3yMh8gz5/EkkLdoNX/Kzt3hlgHLBXTBRffbRj0hBRiDwD19E7wxAi1sQOxafZHyprIVIfAxiER7Gx6qPEhwF8hZn3AnAAKnE1FjROjES0OaqZeE/pq0VmPoOZdwZwEYA/daVj5lXM3Knk9mb6JpPvfw2VmpHgI8aEZtSAuDS+fPslW9856JNZTOWhJaEaK7gzpFadJNLtjj4oc8QcvTZMjnKzjsUkapvwVRWARtH64jBsgtesCbUaCQbhd5iv+oRARE8A8CwAnwAAZn6UmX/e8CEM0CgxEtFCVKR4ETNf7jD5LIBjo5zG1IoOMomaAdx1o8YoRSk+DRFG3sQaBQR4uh58aYUbWFKk11JvSGWrV090XTPPLNtBxFyrANG4lPGIbai8KK/PiG3/HNS/591/6uPDnzQVXL09ad3xSJQiRgBPAvBjAJ8kotuI6ONE5Gn+lEeTD18IwKcA/IyZTzG2787MP6h/vxXAs5n5ZbIv480Xq2CphjVoSUarBm0iMRSR2N8ZG6Ov0KeQaIzdGNMO9d3GVgZ2mphrHJtHKUiVkMZWk15Tmds+/7EDvq+X9fBlv87GfGlvqcp2b7pHXAyLiDoAvgngcGa+mYg+DOCXzPwXOTFqESRGIpqHqn2/DNUU43cz8wNBx0TPBHAjqtl319ebTwdwPIA96233AFjJzP8j+xImkXDBRSK+m6dUwZeIKyaNFc8K7uBaCszo7iKHGJLV3Ei+fG0fjhs3akyelQ5InHIuRrmVIEBf+tA5jyE3248mfzXyn0rv29mEP9fbSWW7H/1naDGsJQC+yczL6/+/B+A0Zn5hToxaeImRiJ4M4F2oFqX5ASpZuwmAPQA8DOAfAXyKmdc7HZQMclmHsfZonXGsUjPgfQVP8hdTuEP2vvTSzeWyVfg3j1Uc7uQjTvtbqhik8yCpImX8IZiTa4xMvGvFYg4mN99r9g5ut44nuhJwxZGDZB/5xLhPZ1O+qLeryvZgWhMcrkNENwJ4IzN/n4i6ABYxc/P9AZCJ8WIAHwNwI1tGRLQDgFcBeJCZP9V4kK5JJBRqKwoSsUhqIKTUbHuNTSkF6yN3H/GFzqOt5BzLIYiqyfiv7gJxpbcwUNU+e5P4DIIbmeFcIm/JrwLikrC2T1detq+IJ/8jvp3+84lx786m/Onebirbp9FdGmI8EMDHAWwE4L8AvJ6ZH8yJUYvZeCVQqxi1qsOntsYNsaAa20NKMPV4fIrNjsmyMWeKARzrq/iI1ucfgVcsQ8dg+RLJRSA5URUqSVplq6lQDZukfnQtUQ5sShDjZvzJ3p4q22fQ7VM9wFvTx/hSx+ZfALiTmdc1EtVIDHUfo0/xmN8+hG5SyT4ER4Ecml7KRxq+G9faL94YDj/i++IxxxXKS1JmdgxSni6CDsUcQ/5S3j5/oZhsm1Ae2cQl2Nh5hypR2x5ACWLcq7OIP9HbW2X7TLplqolxgcLmeADPAHBd/f8IVE+L9iCi9zLzPzcUW4sWLWYMyqE4Uw/NOMb1APZm5mOZ+VgA+wB4BMDTUT2cGQuC003Zv6XaPVZBnp25AAAgAElEQVQtGSpIzNOw9T7EcaUR/AAYGSs48qqfhSGl6lLW0vnzxepTch7F5F3Myo7F2me+CTM4zhjl5MojkH7o7Rvp26fKPOdGvRa3L/6QAtXu63bd4x4Lo+QA70lD05S+k5mfYvwnVM3o/YjoNmY+qPEg7YcvDV5csckUk0ayNb9T49IQXqy/Ek07A+p+Q885UQ/XiWlWlyw3oWPSNrNd/lLLXEr5BVCiKb1HZwv+u56udfwHdP1UN6U1ivFGIlpNRK8jotcB+CKAG+pR6GN7RWcIUk0b2i7t96kAY5v4bnOM8grtC6k6n1oJ3ag+G8c27wQUDozYYvihjO+VwJHYDH8jy7+60seQiQYhFadNnxtLSvn12TUpJAxMs2Ikoq2JaH+1vUIxEoCXAngmqpVjbwJwmT2Ep0kkr/kiIVJxncmP4v200YgP7fi2od8h8otVcilwEb9rP4ThIS5fVtrgf+H4gsNpfMhVz7n+YtKVUn4ee+9DwBHkK8bdOlvyeb1DVbYvpqsbV4xEdD2AF6F6lnI7qrHYX2PmPw+lDSrGmgB7AL5cv9p3BYDNcwIuAg/BaCYccCqhQAF+P23ktIkiRXu7la/zfVzr2OzZaIIkZPpw5auw9050YfwWZ5oxCHbe/af6ba1zNfQutuPce/vNSqgpjT9fBedL58nH2RdZohLoKvq6C6PktGMFsGU9cc1LAXySmQ9B9cJKEBrF+CZUcm0bZn4yEe0O4HxmPjIzaDVGJqqNqWXNfea3FjZ5aZZOkGLS3EwxSiGUzviOnlPS9u0iWMN+xVmrce3ZR4/mXX9G3ie3/KvWpNZeY995jlV/sfYpCOURW65TYl66Ovtd6Sd3tuIP9p6lsn05fWkcivFOAM9HNWfDGcz8bSK6g5mDTWoNMd4O4GkAbu4/aLEfyDSNYFO6VA2bki6FpB22Q6TST5/iywPvgxANwYcqJCMPAG7yzbm5tcefU/nZ+RYmQxXpSzAr55UPjSrpnJi7XeDyg8HfuTWLGJ/U2Zo/0DtCZftK+sI4iPFlqOZ7vYmZTySiJwE4tx5dI0Lz8OURZn7UyGwBYmbdLoUYhVUir5Sbsd8scqnEgBIcvNImKb9Q/q599XbvioOBY7WHstgTt5r/15+/SB6mI8XuOmfWb9ewmn5zeuShmMKfiaH5Du1rKcWvJCbx/IdgxSUuqVtj5FxJcXa7wB33h+MIgEF4BBupPmPCWmben5lPBABm/i8A52kSahTjh1A9ff5jAG8FcCKA7zLzGVkhRyD64UuoEGjT2z4k8vGpHo0Pl78a2kkJYleVc06W4IhlZNIFx3GsOGv1gNj7ambQbHbN12gSHRQTvrrylQjdIDbVWkA+3x6/we0S8QUqSBEaO9f50aQr0JRe3tmW39M7SmV7PH12HIrxVmY+OLTNmVY57djxqNrqBOAqAB+fyFPp0AUWyMd5E8bcAEIhS37PV4tYn6bigf8BUfSUZor8+mtUrz9/EVactRoAcP0Dz8Eli4/DsWdf4fXpJOdQJROLVALKITOtfaDseisnqQKJyi//qfTyzrZ8Rk83K9gJ9M+NEWO9NsxhAE4B8LfGricA+ENmPiDoYyYmkUgZruMjMqkWTVR1xWxd9j71UzoercoxVaS5zZhYwlSMfTjHIxrKcUCokdN8OY8hcC6C6lqrAFNJukTlGatqg8gnxl062/G7e8eobN9CFzZJjM9G9erySgDnG7t+BeBL/YmyJSwQnN8JoS9R82SnRYsWGw76A7wnDWb+GoCvEdGFzHxPig/p4cvRAI4B8JX68+r6cwWAS1Myy0ZAEajGg3l8qN/JdeWf2ERS5VF/ezvvfXlomqGGYnP6q/ededbpg+2uPscz+dFBX+Lpi88ZpL178b44ffE5OGLxdUM+V5y1erSpb8bSDawR49rmu96OtEm+XTY5ilHKM6fFEmoNNIjq4cvGqs+YsDERrSKiq4no2v5Hk1DTx/h/mfnw0LYmMTIfo6+pIxWWUCEXMDR2MbbJ6uv3kZquQhyAu580dnylpmkuHffQ0+i6GXzJ4uPw9zgRAHDNumNw+Q5/AAC4DQfi6zgM1z/wnMretcaLq8vAOpaoCVp959lhO9TXqilL8PSJ+mLo24fWtkkhWm15DvrOb0rv1FnMJ/deobJ9J31kHA9fvoOqKX0LgN/1tzPzLaG0muE6i+r1W/qZHQYguFoXEe1MRNfVC2XfTUQn19vPrRfQvoOIPk9EWwUjWHvf8H+p5pYKXMjOBZ96UaQTt2nisbY53x02FWXopjN92WMNLdikuMIsw93u4K2Uvp9LFh+HffBdXLPuGFzzwmPwvR12wXM3uRIvXXclvo7DcM26x/ueBv2PDmKUllN1qmZJzdn/PduGSFEqF0Y8Up+pC87jstM5tnnf5Ar5csUesz0RJd+VJqL59QqBqxPDeYyZP8bM32LmW/ofVd4KxXgIgAsAbImqz/EXAN7AzLcG0i0FsJSZbyWiLVCx9ksA7ATgWmZ+jIj+CgCYWZy+TP3wJUXR9X/b+3LyEG7CqLjM/1J6S20lPagJbe92h1VPfx8qortk8XEAgJeuuxJ4NtD9Xm1yDYA1wF+c9G584IFqhJfzwY1r/KPv5jdjk76ltC5Y59HpK5A2Kk9PBRD04TtHrv1q5CvGHTtL+C2916hs/4L+RrO0wZ8D6AB4AjMrF30aSt8FsA7A51FNlQgAYOafhdJq3pW+pX68vT+AA5n5wBAp1unW9u2Y+VcA1gDYkZmvZubHarNvoiLKfPjII1SY7bShG9ORx0j/ZKjghuKQbHy+uo8r26G3ImKURyBP30w3ZhP6ezvsUq0IDKB7Mgak+P6zP4DTF58z6IM0B1Sb8Q69/+yK0T4eH2zytgame9N0u86uClVlaMWovuZaWzsmV3wxFYGJpcvi01ioHr4sUH1CIKKdALwQ1ZovqXgdgFMBfB2VMLsF1bwP4fyFxbBeA+CzvlUA61UElzLzTcFMiJYDuAHAfvVL3f3tXwLwr8z8GUeaEzCQiVseUg1JquEhJ7H/KakWLQChoHpnj9GoRNt/MeUQgKGq+uiPWTwJ/4B98F0s2+QerHkE2JJ3AQDs+8DdI+MYzX69YP+or8KLVYaS8uwrRbsfUDq/paE5zlQ/5j7zu48CA7yXdZbyG3tvUNm+jz4QWlf6UgB/CWALAO9IUYw5kBTjtgBuI6ILiOgkIno5Ef0xEb2XiL4G4EMANOtLbw7gMgCnWKR4BoDHAFzkSsfMq5i5U8ntzYZ3Oi66cwU2340Tql21qtPYbr8eJ6FvO+gzdMWivSFiSdQRizNPB8Gar5mtP38Rjlh8HY5YfB3OPOt0HIav49izr8BFeBWecAXw9I8Ce627BxfhVbh78b74e5yIeSsfGijF/nvh4kgC63/fNviqnu+4fde9f0yhPkzleVbPYxlSidry6UsjVAIjsPvxE8AgPIqNVB8AP+nf3/XHJMWjAazT9gf6UPPVyEeVVupjJKL5AFYAOBzAUgC/QdUkvpKZf6QIbCGA1QCuYubzjO2vQzX48khmfjjsR1gMy4dYYtEitamSipAalOxT94WOsU9U9UOUIxZXywFds+4YzOdf4/TF5+DV+CxOwkedE2PELtalfvNDildDQr4YAkpTjMl3jBpFpznWGCUrnof8PsYlnR35Nb23qGz/hv7C28dIRH8J4LWohNMmqN5YuZyZdR2Yj/v5iPF3EwBHAriVmV8WSiv2MTLz75j5q8zcZeY3M/MpzPyPSlIkAJ8AsMYixaNQrRXzIg0pDiCRQqgg2ra2AoqpWX0F32UXgkR6FvkMbdfcIK5zpL2xhJt4BXec8ylee/bRoI8xLll8HA7C7QNS7I+BNNOsX3Lu0FPuFdwRCcBFiqJqtBRm8Hppz1Uonc/Gdz41ZTUUk4ucpVgaRKk+RmZ+NzPvxMzLAfwRqoe1UaRY+3mr8XkTgIMA3QwWmuE6qTgcFeuvIKLb688LAHwUVb/BV+tt54tebPgKsY8she3rz1/kvsFiClgoHqlwuxSBlS5qMLKkWmLJ3/ruE9m1Zx89eFjS7wq4/oHn4MyzTseKs1bj2LOvwG04cODq/Wd/AMDwwO++n8Hv/sxCEqxYg6/0AdFDarz7tdsC/p0z3mivlcu3poIeM6Z1aYMaDwPYXWM4N96V1tbwLvv6d9Ig7kSoZ5WxISnjmBs6dDN23e9De0nfUxH055gczL5To/gC8lq1JFRK0gxD4jkLnU/hGFZwp6ogYmJ35WXnk4T8pvQOnZ35uN7bVLb/QG8fxwDvL+Hx15rnA9gbwOeY+bRQWlnTtmjRooUSDBrnsgUa/LXx+zEA9zDzvZqEmgHeHwDwIWb+ef1/awBvZ+YzE4ONhksxDpp25gSvPmhqc9f+EFJtU1SpVhWZ3yX820rSmlrMTGvOQm6qRKcyciktX+yGnXqORc3x+ZRvjlrV2Hquk+ucRpcZbRkYsctXjNt1duFjeu9W2V5IbxnL8qlEtBjAU+u/32LmdZp0mj7GP+iTIgAw84MAXhAfYiasgnEt9UZJ0WwaWcNABkNFTDv7t+UjWMhitvtuOim9EKd3xurQjRkLK4138HO3W70LXf82m87eV+4059jc13W8nimRbYiwYoin3reCOyOvR4rlSYrTgrcrJ5b8UyuKTEzb8qlE9HIA3wJwHICXA7i5Xu4gnFahGO8A8FRmfqT+vymAHjPvmxV1BFR9jK7fLhtXLR3qW4whNVcaKS57n0YpaMk1oL40222FeOZZpz++jKwnxpEZeOpv1WTBjvOVNH9i6HxL9hoCT1V0kYidmd2LYHz5inHbznI+qvcele1n6fhxTSLxvL5KJKLtAfxbkYlqieidqNZm/SSqjsw3APgiM38oN3At1DN4x8BHphJx2LYh+5ibVePTjgWOWV5cJB573nzpuv73iEce1thpXJNWpBB3auXkSDvUHdM/jhwi0l7zmIotVLFqKl5VnvnEuE1nV35eLxBHjc/Rn4xllUBz0b56NYLvaBbyUz2VrscePhfV0gZXM/NVGfFGI2kG7z4KkKl5Yw9eZWtKPZr2ru0+fxrCDt2kwv6RQdmSPR4nyqQJGXzwKTfbZ6hCsLadedbp1bCiVBLLsYejb9Hn1/6OPYdimnxi3LrzJD6y9z6V7WX0mnEQ47mo5ni4uN70CgB3hCatAfTjGNcA+Aozvx3AjfVsOZNBSmHISY/h8XBDa6SECrH9P1YpWLZD/aSQBzCPDML2kalL4Tj8SdNmjbzWh8eby87FsCzbPvGaCE4cHCBZ14qCvuPvj7UM5qOBLy7PeQUCfYtWrGJcMfs12yMR+UpgYyCi3YjocGY+FcA/oiLHAwB8A8AqMXHfh6Ip/SZUcm0bZn4yEe0O4HxmPjIr+giMKMaQWrKRSUimTXDVPKlwA6PNSl+MoZi1KjEAr6IzfUox+tLF2KcqSd81S1FTHt/q/k1NWQj5yIGy/DapGLfs7MaH93Q9bFfSsU2u+bIawOnMfIe1vQPgLGYOLkyjIcbbATwNwM3MfFC97U5NO70UsprSgL9AaLdLRKxRh5Jy0UBDkvW3uo9MoRyj5nV0+YuI20bRvj5F/s4VE0MqLiXf1AotYKNa8VH0vawIMR7aUy3bjKvpxU0S413MvJ9nn4q7NE3pR5j5UcPxAgiLZDUC7VxxoZpS0xzR5KEgAfEdZ4k8fXkqY1K/BmekcTa74ZltxoGhZm+3Gs848v6zL27PfnFYTggh5Wtuq7fbD2Ds/cG8YvK098WQuAeqZXATfWsxRcN1NhH2bapxoCHGrxHR6QA2JaLnAbgEwJc0zotBMyVSSKXE7FMQS2jqK9X6K5KqdNlqt/sI23Mz2iRkE50zPjO9NSh5MMY0hFCMPnvHMajyUtqIa4S7YtBUuj4/seo21ybFVokpIsZv112AQyCi41FNVhuEpik9D8DxAJ6P6qn0VQA+zmN8yXpouA4QdVHFCWy1hdKVNrTN3m6osyLrQ8fG4osphYi0CjZ0zFYsqgWmQjFKdua3tC2lbFj2drlzvrHj+i359l1v20aIy+8/vym9RWcPPrD39yrbm+j5TTalF6NazuBRPE6EHVQz6/whM98f8rEgZMDM64noCwC+wMw/zoi3RYsWcxjrMQ+Pjm9pVC+Y+QEAhxHRcwD0+xq/zMyqpVMBoSlNFbpE9BMA3wPwfSL6MRHphrY3BamJUv82X9nqj6UbSaNUA941l301vvS7O7xuSv/1MufQFK0KCp0PSR2Gmn6m6tPOmm2rm27XvVyAI/+RoT3SdUpoPXiVoctOq7RMe2O7rXydx2bH5PIVsg2lK90ED2BKmtIAAGa+jpk/Un/UpAjIfYynoJpT8anMvC0zbwPg6QAOJyLd3EKloWyu2v1b6uaZowBGN39czWdzSn4rzmupV03cepa1QqR0E9h51R/nolzCzeMafzjyv/7WjkfUNk81E826iHSFr/WlrVRChCo1W0O+fOQvEXAsfHml+C9IiMBU9TFmQ1oM6zZU7xn+xNq+Paq3Xw4aQ3x1nsJwndiaUiINB8Rp9V3pfMQkFeKQSggpGx98xycQj73dO+5Sysv132cfUj0hHwH74JAj128pP6V91HCjUMVi2xYmtAr54xg37ezNu/U+rbK9i542ltl1UiEpxoU2KQJA3c+4sLmQIiHVlAqVFSpk4pAVidB8v7XxKm+AoOqLUUaOG9P7GqCVfrAcqhSLKwZfrNZv50Jjtc2I+jR+R5GiZpuPvKzji5qI1zx/imumXjrWEVezKLd86qQhKcZbmfng2H1NIGkxrBg7DTQkk+LT/C6BlFhSlVqfmBxPk71KO+ArSzH6fFh2qqV2pUo1JibXfl+8moox4lpJqnV0X75i3KSzL+/SuzhsCODf6YCpVowSdR9ARL90bCfIAygrI6KdAXwawBIA61GtG/thIjoOQBfVNONPY2bdcP2YWtCnhkL2vv+a9DEoUOiz4jCUh7fvMKISGHrXt680z/fbe6d668cUWvLBp6ykmH3xCsclxpCKkHJOSe+BtF7Q+iUReSrBIDzS8HvQ40Jja74Q0VIAS5n51nrSiVsAvATVWzPrUb3c/Q4NMY4oRhdSVU9q+li4yNrcZ393w/17Z/Kj1dyIPrWjVWFa5eWzF3wMrfci+Q2pSimeEpCujyad+d/81qRpEmZZEpv3+Ypxo85TeHHv8yrbe2n3qVaMja0SyMxrmfnW+vevUM3QsyMzr2Hm70c7zG3C+Pyl3nRW0yw6P3ufg9icr/cZafqkODQzuRmXscB9FGw/viFLoW2wVgC0iGfe/af6nzK74km5VgYxeP3aeeSq9pKti1w/ZlkaA+bKU+kml08dgIiWo1rT9eYiDiXFFWpe+Qqur6nmUgOWn5GmmabQhlRjhP1grKajyRgzh6Lv/Imz70ixBo5p/fmL3MseSETm8yep8e7ouELAUXE4iNuLBJIWx6uGyq2Uf6jca9R4AWwQw3WKZUC0OYCvATiHmS83tl8PoSlNRCegP0Znyycegl+8XpehR32JdrHQpBVuUtHeF3+A3Mymkji1mcPXvJUP6YbkxCJEAtrKQCIPQxFGzQSkuTYa9RdTWcSUxZCCTfFnfo8gvym9oHMAb/WtK1S2P52/01Q3pRslRiJaCGA1gKuY+Txr3/WI6WNMXdpAU0Mn1P4hOOc5jLjRnTd9DJFIRBp7w2sqG60SioFWIfr2K/Mv/j69md73X0pnfucgKvZ8Ypx/8IG82U3/prL99aLtRWL0PbzNiS8GC5pyTEQE4BMA1tik2KJFi7kHZsLvHivWTH4M1TLNg4e3RPRVZv5uqQwkNNnHeDiA1wJYQUS3158XENEfEtG9AJ4B4MtEpF8/JrUGltRPjC9XDJ58vAtU2elccVj/g+MBQ4pKsFP15QV89LePzGPoOb/OZQxc+0LNRQkalVvDOWRJq9y0TVapeSv5jokjFEvBFpETDPzusfmqT9CV5+FtswfwOBrvYywB9Qzekc0X9eLtPr8+klVANV5PgmErNtuFdOa2pCm/JN/apqjmvGnOt0m8sccSe90NBKcXi0Fqs13r2/weQX5Tmg44mOddfaPKdv2Sze8BYL5Zt4qZneux1A9vbwCwHzO7xlYXx+wQY2gcI6ArWJobNVY9uhSVT61oSajETRKjXiR7Y78zdokEfaonZK9RqRK0ZUGKxWcbk19MHjH5xSJI/iWI8RDGlV/XGe+4ierhi+/hbdOYHWIsva40EK0wnWmlNKGbPnTzuxRYShyhvGOhPU+peYVizSkHdVp1a8HeJxGddL2E/1Hr9DSlKJeuBt/XyyPG/Q9hrFaOyNtlYZAYpYe3TWO2iLEENGQYU+A929Sry2n897c7fktN8n5/Xc5Np17V0I7LVpY+FagkWFGpejAUuxRvaHvgmqvOkaZS1BBpTAUXRaQFFON+HcblygW59qTQU2kC8CkAP2PmU3LiSsHcIMbUmlRSW7FqR7LPVU6+7SG1KKkbX36uOPsEYC8da6cL+bW3FVBAWX16NnzK0Laxt5vErVWUvjybUoRB34WI8XNKYtw3SIzPBHAjgDtRDdcBqiVRdQMlMzE7xNitF8RSKgyxQAuFpNEbrVShNxThkCIy9yf6jLJJJXwpzybJIRRHbN4R1zepXEk+6339VyoHr16mnrsSTel9O4yLlcR4gEyMk8bsEGOppnQfpW/KkGI0m12hhxchklPcMMGYfLFo8tDCiqX/TvdQszNF1TpszJmC7GZ31AgAjcoz95nfmjSuuJt6jznqGhZQjPt0GJ9REuMhLTFmI/rhS8pNrVFDBYgz9PqdasbpnHMQeQxJDwbspqVUCTWt3EuVhRIq3FT6GqLOLW+hymZoXwFi3LvDuFBJjIe2xJiNZMUYKliugtgEJJKqv6U5CkdiS2nClmwGl2ru5vrxNWU11z2G+DQqXCBBQLnOeC6yzmcBYtyrw1ilJMZnt8SYjSRiNGrL4LjBUNNVm58vvXQzufKfNOn41Iv57dtWArbSkZRPjNKKqRxKqUYXYkhcyj+mYhfOSXV/7J5PjHt2GB9TEuOR002MCyYdQIsWLeYIGNUbznMAc1cx9pGiyjQKsoTKlPJ3Nd9C9hEosvpfjGLsWg9fIPS3hlR+jsqK9ZHiOzXvUJM9JWa10i7QlN69w/iwUjG+cLoV41gmqi2GEEnFNLU0edkfV34+m1Bskq21zTfRwwruDPIdmQk7cIzi7OC2H+k8RNyo689fNDS8aIQUjXPl7G+17TzwTmYRQuz5cKWTmshWWQhOmBFIH4zR9dvnuwT6ilHzmXLMpmKMLeyuAqat7XOVRUx+kQosKq7E4/A+QfXlUTJ/zY0sEZRBQElvIflsUmPVQKpItZWSZOstYwUU424dxoeUivHY6VaMs0mMGoQKjvmdYpNDmKkkkaLwUuC7sTSq1+XLti0Ro+nHg6j3oQvDS8YhEg4QfXNjTwsQ45M6jPcrifHV002Ms9WUDkFSXapa02MfsnWllbanKLf+glfavGOJR2oOhv77/Lj+C/mu4A5WcEe/oHzg+pivL7rWxAmeI82xGKrU3D+0Dk/s+RMqIpXyHRP5j6BtSo8Xjbz5kgOJ4CS15fpINqEYQuTf5A3iy0vId2RSC+kc+PI0v337JRvomtbehz+aGALk5k0j+W+c8AooxuUdxnuUivH46VaMs0OM9nyMroKUW3BifQRu5qi3NzLUpPpmDfkwbIq/2xtKE6tuBfXm/e/zEapkXD5dZKaBlS5qUt2MMh7uVihAjLt0GGcoifHNLTFmo/ibLxqFI9Xg2hvQQ97qoSiTUg8atSPs877F48vHt83OL0GlmhiZ6TxWlccQqM+Xwz7qtUvX71CeKhQixncpifGklhiz4SRGq5Cox+W5kEA04s2vvens36EbUSIKDSzCcaoIl2/ptyu2EmjSd0lIqtE4VyryK1ThpT10KkCMT+ww3qEkxpM3UGL0LX9IRNsA+FcAywH8EMDLmflB2ZdyaYMcaFWkj7gklePykaM+YmGrVNexho45k4S9fl3/U9Scdr/reOxtodgtX43OkGPGoa1sk9RkAWLcqcN4q5IYT5tuYmzyqXR/+cO9ARwK4CQi2gfAaQCuYebdAVxT/x8/QjdVpN28lQ/pnhxr87W3SWRqYeSpbrc7vBJeSIXa+dr7XPYuRBJM8LfD99Bgblfemv+hYw1UkvbkveJT9ZSKT3ndoyqUJipgBvA75WfK0RgxCssfvhjVlOWov1/SVAwtWrQYM9rhOhGZGMsfAvgRM29l7HuQmbd2pDkBg47FLQ8BrGUfYpoNk4LUnJP+h9Ln7LNtNE3LUnmatnY+IbtYaM+BKx5tGul8Sb5j8gvF4lPxmus5hAJN6WUdxhuVTen3TXdTunFitJc/JKKfa4hx2MeYxjHG3tyxN1RpEh93DKGmZqz/btoiVz6bkbkPtXFpCF/TnA0116WKMrVpLMFH1k4UIMalHcbrlcT4l9NNjI2++VIvf3gZgIuMNWEfIKKl9f6lANY1GUMUNAXIZaO1jySNKBttPNL+Okax367bHX6jQ8pXcQxDfZ81Bn21GjVp+Rr408QSIivNf+k8aM9H7HYD8+4/1f1Wj+1HU0ZzsR7Ab5WfKUeTT6Wdyx8S0bkAfsrMHySi0wBsw8zvlH0pFaNUIEurNU0sktJoQiE0DUmFlYo9R9357M1vlx9VU1MZR0p8KedNpQIjUGIxrMUdxiuUivEjG65iPBzAawGsIKLb688LAHwQwPOI6AcAnlf/j8LIO6h9aGv/EGLUl88+VSn6/MamSVFBLvuQIsxQO9ExpZxLjZrTkoymImiiZZCTVut/7X3pcfTRvis9XhTpY/TdEIGCEzVllbAveaxbipqx44pUYaqZYfr/Hd8j/YYB4ik2+FlTKcScE58vNLy6n4TGWhUF+hi37zBeolSMH59uxTg7xBizrnQsSo8Js9cAABFiSURBVDShfLbmt2J/9CqBGtUWsreJVyJzH0lHxLCCO9U6yFL8ElmXIAdHxeE63sEDHV/f5TggVY6aMqaKuQAxbtdhvEhJjJ9siTEb2UsbSESgVWOSn1CaHLUTqxhjjqUkyWvyiWmyas51yfOvgUstl4DvHDWRhxcFiHHbDuOFSmL855YYsxFFjC6l4ZpC31FInO8/h1SU5maXSM+3PZY8HGlGpvmKySOVzD3nTXOTi0vIaohVey0iCCe6yeyKPRRTU4iqTJflE+M2HcaRSmK8tCXGbDTWxyjZKojU6d/3P9a/Rj24KgHtLC2SYk5RLhpfJdJoUErZCRWk1zaFeFMnP/HFkIQCinHrDuM5SmL8fEuM2Sg67VjpJkpOntqbTWOvISABgxtUikfyJ23TVkqx9naaWDIPQD1LTWp+CjId27yYJYhxqw7jWUpi/FJLjNmY6AzesTdcilpJVTixzTZH/KJiCR2vT03lKExNPiUwrgqyVOw58arSFiDGLTuMw5XEeOV0E+OCSQfQokWLOQIG8L+TDqIM5q5iTOzzCfrTbDfy9S4qL0A9rs/8Dtn6VG+oHy2jv24wLEcTq5R3E03wUN6mX/u3NpYmVGLMOdHYDVBAMT6hw+goFeN1060Y5y4xapDSPyXdRK7tofxj05hpPQTX2GzRMU33GP8p3RWCTXCAeqAyE//HIpbcpe1NdCkMUIAYN+8w9lcS4zdaYszGVKwS2GSfVF9dSp39GcpNyjOZwLS+tCTk8imliYnVRW6J51AcAmXnNQ1Qn7MCxLiow9hHSYy96SbGubWudAxKFlzXTRu62Y3f/dlk7JmgVSQj5aGxi/HXNWbeccWlVZM+/0K+TlvhfJuzqc9b+VBY+SkrnvVLzo17RbTJClXjf5wEPYdm8J5dxZhDHiHb2O2x/jXpzO9SCKk3VzOygRvLVMYjg7q1eUp2keRU/HXDCeUxQJL/Aopxsw5jN6VivLNVjOVR6qbVKJEYO5//WOQqw5AK85D+YG4/TR79/QnHOCBCeFRySN0Z8aoQKC9RbzIllj1nHrZfDbSVxiTQzq4zXgwUo6R4JNg3m/kdwNhnUclRa7lKz6UYUxF7XVL8aUk0Jb9xIiX24sdSQDFu3GHsqFSM/90qxnLQ1rRSAQrd9Nb2pkjRtZKf89veFlKKqYTq8wXEqcgYu4wbe2TW6tzmtyatZlvMfo1dTEUeez5Ko1WM48VQH2NTtf641ESqsikVX6RiFv2EYpQUbK5a9NmZ34H0wSfMOYi9XuM4R6JNAcW4UYexnVIxrp1uxTh7xFgKTdyYhaCeHFdSwik3pv3bQWzRYySbOHclfDZd+YRaACn5NFoOCxDjwg5jKyUx/qQlxmxMxThGE6GmeGwTKrUJnJqvL53r29c8C+SV1TdbQm3F2pZS0lpo8tOqcnt70rkoQIwLOownKInxwekmxsb6GInoAiJaR0R3GdsOIKJvENGdRPQlInpCU/m3aNFizJhDfYxNPny5EMBR1raPAziNmZ8C4PMAlOMtMuCrcWOVgUIheQc/e+yjannXcUTE5/UpNfG0586xbWRZ1Jh4PYpOPTxHgibfCGVcPI5QmZDOve/axMaQivUAfqP8TDkabUoT0XIAq5l5v/r/LwFsycxMRDsDuIqZ9wn7EYbrhBBKMy1Nq3Ej9uYfRz+rNo+Uazbufj2pySv4Vr1ymNvlYscBoEhTmjoMUjaleQNtSntwF4AX1b+PA7Czz5CITiCiHhH1gIerjdqbxvet7WxXqsoVKde1BLnUPlwLrUcNenZ9h+xj9zWBFLWVenypx6Ypaw7fzlcONZWYRl364igJVn6mHOMmxjcAOImIbgGwBYBHfYbMvIqZO1WtspnfY6AGjmpmxDRhgMdXuhPicY5XzCEhY79rfZSoKc60N1xuvFqUIqFJ+UBExRQDTUU9LcQ4RzDWprS1bw8An2Hmp4X9RLwr7cMkmsaeppJXEcTENhcKeInjjjkPJsH4FF2M/xLN2pjKuVEUakpD2ZRG25QegIh2qL/nATgTwPnJzpTN3RF7ab/rt2SnydNCcMW+mDhSm9HTglRilP677DXdBikkl1JOQtfZRcyhGGa9gpxCNDlc52IA3wCwJxHdS0THA3glEf07gO8BuA/AJ1XOli7z7yvZjNPcqLkEGtMsjbwpnEpU0cSaKKHm9pm59rmIRnttY4lX6cfbpRK6RiFSDylaLaR7TI2581h69gZ4S82g0LZS6NZLFpRaX9pn2wTpp8YTm18Jv+NsVsY20WOb2IlNcnGgfNEyXqIpfRADX1Nab9k2pYtCqyS05JOSf7cbfmXPky7YnEuNTUqrbbbF5iPZaI4lZBPyk6LYU5CjXlPta0SXs4li7ozwnj3FmIpctVZagUb6S1rHJXTMEz6miaLpWEuq3VCzuwhKKMYDGLhaab1kqhXj3CfG3MJkq7xcxRWbZwTENaJLwqV+J0GKOXlKFUiDXTDTW3GUIMb9GbhCab3zVBPjgkkH0KJFi7kCxiw8WNFgdvoYU5+8afrbJMWg6OMaee85Fyl9jf1+z5J9sEJeqj7TpuG6btp+SN8xJJSz/hNn15tIWh/F0ESfsRptH+NYQbSM0b2v+lOyzya0LcWXcZMVe6IYsJ3Y8gvjzi/1YVGI7CzfwXWpQ34NjLz7XKqcFU9Xoim9LwMXK60PmOqm9OwQY+okEqnIfVhTCOMbrpHgv7RaHMf1lSpFDfFNSiE33r9dghj3YeAzSutDWmLMBS3rMNYePekwyqC06kn1W/rGSlBTY4e2+S+RotTtkhJLrk0xlCDGvbmabVCDQ6eaGGejj3HtfWnpYomilC8pXTdhXkFNv1dsf1/KudH2nzV1M+f61RJ3zPmO6ZMMnSOXks3FWCuoudPHOBvE2DS0tXdKWse+Iv2BknppQg2a39OE1MqvFAmlEmOMrxyMVXXOnVcC5x4xllYuOc2dJptKRrp5Kx+Kzzsjv0b9xiq5kv2b40RqRVs6r6JgAP+r/Ew3ZqOPsfRiWIHCp37KO619aTZKxTkrxysh1H84Teq45PkO+irRx7g7A+cprV/U9jFOHQKFTZxB2d6e0lcVg5QHNU3FofWTkt+kFamrG6Jhcgz2NTfRHdIo5o5inJvEmFvTum6m0FNKjZ8SsZS0n7Str9uj5DktWQ40NhH5jW3s6djUb0uMU4eh9VdSidGnFHxPIGM705tSUiUeHOTePCmEJVU4vm3T0pyPIX9tzONSzY3lN3eeSi+YdAAtWrSYK+g/lZ59zObDl6ZVg+2/qSeF06B8YjHJuEvk7VP+iX4HD+rGUWZKYiS+Eg9fdmHgNKX1ie3Dl6IYR4HzNIFHJgko2VxPxbQMM5H2lepCaKrPNqNMDb0DHetzksTZSN5zpynd5JovFxDROiK6y9h2IBF9k4hur9eMDq4QOIIJFraR2WtC/Y6pN3fsg45pUSazQIyl8ku9jiX7dacOZR++ENFRRPR9IvoPItJK0SJoUjFeCOAoa9uHAJzNzAcCeE/9vyym4YaS9jXRCT+uGyyHKMbxhN580GFWGKnE3MR5bfoh10TJtpxiJKL5AP4ewB8A2AfVQnr7NBK2K/9xritNRFcBuICZ/5WIXgngGGZ+VdhP4QHeTWOaVJwGhfvdZg4b0rF6UaKPcRkDb1Zad8U+RiJ6BoAuM/9+/f/dAMDMf5kToxbjJsa9AVwFgFCp1cOY+R5P2hPwOBvuCeD7jQVaHtsB+Mmkg2gIc/XYNvTj2oWZt8/JiIi+UuenwSYAfmv8X8XMqwxfLwNwFDO/sf7/WgBPZ+Y/zYlRi3EP13kLgLcx82VE9HIAnwDwXJdhfZJWufZNO4ioN81P3HIwV4+tPa58MLPddZYDl3od2xCacT+Vfh2Ay+vflwCIf/jSokWLDQH3AtjZ+L8TgMT5B+MxbmK8D8Cz698rAPxgzPm3aNFiNvBtALsT0a5EtBGAPwLwxXFl3lhTmoguBnAEgO2I6F4AZwF4E4APE9ECVP0LM/REJQoz2QWgxFw9tva4pgjM/BgR/SmqZxLzUT20vXtc+c/Emy8tWrRoMU7M3psvLVq0aNEwWmJs0aJFCwstMRYAEb2NiO4moruI6GIi2oQqnENE/05Ea4jozyYdZyxcx2Xs+wgR/XqS8aXCc70uql8/u6t+nXXhpONMgefYdiWim4noB0T0r/XDjBYCWmLMBBHtCODPAHTqgezzUT1B+xNUww32Yua9AfzLxIJMgHBcIKIOgK0mGF4yhOO6CMBeAJ4CYFMAb5xYkIkQju2vAPwtM+8O4EEAx08uytlAS4xlsADApvXT9s1QDUt6C4D3MvN6AGDmdROMLxUjx1W/w3ougHdONLI8jBwXM1/BNQB8C9W4uVmEfWxrUQ2Nu7Te/ykAL5lQbDODlhgzwcz/A+CvAfwIVSH8BTNfDeDJAF5RzyJ0JRHtPsk4YyEc158C+CIzr51kfKkQjgsAUDehXwvgK5OJMB2uYwNwC4CfM3N/5oZ7Aew4mQhnBy0xZoKItgbwYgC7AlgGYBERvQbAxgB+W7+O9U8ALphclPHwHNcfAzgOwEcmGVsOhOvVxz8AuIGZb5xEfDlwHRuq2WlstGP0AmiJMR/PBfDfzPxjZv5fVK88HoaqZr6stvk8gP0nFF8qXMd1NoDdAPwHEf0QwGZE9B8TjDEFvusFIjoLwPYA/nyC8eXAd2xb1U1rYMyv1s0qWmLMx48AHEpEmxERATgSwBoAX0DVtwNUr0H++4TiS4XruM5j5iXMvJyZlwN4mJl3m2iU8XBeLyJ6I4DfB/DKfr/wDMJ1bN8FcB2Al9U2rwPwfyYU38ygffOlAIjobACvQDUD522onmhuiupJ5xMB/BrASmb+zsSCTIDruJj5EWP/r5l580nFlwrP9XoIwD0AflWbXc7M751MhOnwHNuOqEZFbFNve415HVuMoiXGFi1atLDQNqVbtGjRwkJLjC1atGhhoSXGFi1atLDQEmOLFi1aWGiJsUWLFi0stMS4gYKIfkdEt9ezsHyJiLyTQhDRQUTERPT7Ef7/PyJ6lmP7EUS0OiPuf6vf8GjRojG0xLjh4jfMfGA9C8vPAJwk2L4SwE31dxBEtA2AQ5n5hvwwR/DPAE5swG+LFgO0xNgCAL4Bz8QC9RsUL0M1jdrzzTkZBbwMxiQMRHQUEX2PiG4C8FJj+6J67sNvE9FtRPTievtmRPQ5Irqjnj/w5nqqM6BaEElF0C1apKIlxg0c9TRiR8K/AtvhqN6//U8A1wN4gcLt4ahmdUFNpP8E4BgAvwdgiWF3BoBrmfmpAJ4D4FwiWoRKET7IzPsDeB+AQ/oJmPlBABsT0bbaY2zRIhYtMW642JSIbgfwU1Svin3VY/dKPD7J7r9Ap9aWAvhx/XsvVMT6g3quw88Yds8HcFodx/UANkH1CuUz+3ky810A7rD8r0M1e0yLFo2gseVTW0w9fsPMBxLRlgBWo+pj/DvToFaTxwJ4ERGdAYAAbEtEWzDzr0Y8Gr5RkVwfvvdOCcCxzPx9K18KxL5JnUeLFo2gVYwbOJj5F6imw3+HY52T5wL4DjPvXM+oswuqqdRCM0CvQTU9GQB8D8CuRPTk+r+pOK8C8NY+ERLRQfX2mwC8vN62D6rlBlD/J1TN8R+qD7JFi0i0xNgCzHwbgO/g8TVdbq93vRLVXJImLgPwqtruCiJyNWm/DOCI2vdvAZwA4Mv1w5d7DLv3AVgI4A4iuqv+D1STxW5PRHcAeBeqpvQv6n2HAPimMSN1ixbF0c6u06IR1CR4NDP/PCHtfAALmfm3tdK8BsAezPwoEX0Y1dIK1xQOuUWLAdo+xhZN4e2oHqREEyOqRZyuq5v2BOAtzPxove+ulhRbNI1WMbZo0aKFhbaPsUWLFi0stMTYokWLFhZaYmzRokULCy0xtmjRooWFlhhbtGjRwsL/AyoiII0p2/R+AAAAAElFTkSuQmCC\n", "text/plain": [ "