{ "cells": [ { "cell_type": "markdown", "id": "410e24b9", "metadata": {}, "source": [ "# Example Usage" ] }, { "cell_type": "code", "execution_count": 1, "id": "9b15c2b8", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import matplotlib.ticker as mticker\n", "\n", "import ppops" ] }, { "cell_type": "code", "execution_count": 2, "id": "e6acbf03", "metadata": {}, "outputs": [], "source": [ "# Create PPOPS instance\n", "ops = ppops.OpticalParticleSpectrometer(laser_wavelength=0.405, laser_power=70)" ] }, { "cell_type": "code", "execution_count": null, "id": "e9a1607b", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk4AAAGGCAYAAACNCg6xAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAVA1JREFUeJzt3QdY1WX7B/AveyiigIADELc4UNx7ZO5VVmrlyizL/g1taHub1WtpqWmZs2HmyFJTK/deqLhRFEQUQUWW7P91P3COrAMHBc76fq7rvJzzO79zzsMPX7i7n/u5H6usrKwsEBEREVGxrIs/hYiIiIgYOBERERGVADNORERERHpi4ERERESkJwZORERERHpi4ERERESkJwZORERERHpi4ERERESkJwZORERERHpi4ERk4RYtWgQrKyudt61bt+r9XleuXMH777+P4ODgAs/JcXk/Qzh58qT6/IsXL5b6e+v7fckmDb/++is6d+4MT09PODo6ombNmujduzd++OGHPOfK+8n7GgtjGw+RIdka9NOJyGgsXLgQDRs2LHA8ICCgRIHTBx98gFq1aqF58+Z5nnv66afRp08fGCpwknF169ZNjc0Qpk6diunTp2P8+PF47bXX4OLigkuXLuG///7DH3/8oa6Pxp49e1RQRUTGh4ETESlNmjRBq1atyuxqSCBgqcFAcnIyvv76a4waNQrz58/P89yYMWOQmZmZ51i7du3KeYREpC9O1RGR3lasWIG2bdvC1dUVzs7OqF27Np566in1nEzptW7dWt0fO3asdqpPM8VT2JSWZH8GDBiAv/76Cy1atICTkxMaNWqkHmumEeVxhQoV0KZNGxw8eDDP6+Xx8OHD1fvIa+XriBEjVCZHQ97j0UcfVfe7d++uHZcc1/jnn3/wwAMPoFKlSur76tixI/79998C3/+6detUJs3BwQH+/v748ssv9bpuiYmJSElJQbVq1Qp93trautipsZ07d6J9+/Zqiq9GjRp455131BSfnJt7ClJzTf/++28EBQWp6yKZxB9//DHP+12/fh3PP/+8yihWrFhRTR/26NEDO3bs0Ot7IrJUDJyISMnIyEB6enqemxzLPX00bNgwFSxJrY4EEe+++646T8gfaZnuE2+//bY6X265p6AKc/ToUTWN9cYbb2DVqlUqKHv44Yfx3nvvqcDg008/xU8//YS4uDgVEEj2RkMChgYNGqhszsaNG9VUWFRUlArgYmJi1Dn9+/dX7yFmz56tHZccF8uWLUOvXr1U0LR48WL89ttvcHNzU7VHuYMnuT948GA1xSbf/xdffKHO1XzPRfHw8EDdunUxZ84czJgxA6dPn1Y1T/o6duwYHnzwQSQlJakxfvfddzh8+DA++eQTndd08uTJeOWVV9Q0YLNmzTBu3Dhs375de86NGzfUV7nO8rOU70N+tjKdWZK6NiKLk0VEFm3hwoXyF7zQm42Njfa8L7/8Uh27deuWzvc6cOCAOkfeM7/33ntPPZebn59flpOTU9bly5e1x4KDg9V51apVy0pMTNQeX7NmjTq+du1anZ+fnp6elZCQkFWhQoWsmTNnao+vWLFCvXbLli15zpf3d3Nzyxo4cGCe4xkZGVmBgYFZbdq00R5r27ZtVvXq1bOSk5O1x27fvq1er8+v0v3792f5+vpqr62Li0vWgAEDspYsWZKVmZmZ51x5Xq6XxqOPPqq+p+vXr+cZY0BAgDo3LCwszzV1dHTMunTpkvaYjFnG+eyzzxZ57dLS0rIeeOCBrIceeqjI8RBZMmaciEhZsmQJDhw4kOe2b98+7dXRTMM99thjKtMSGRlZKldOpr5k6klDpuaEZD5k2iz/8dzTcAkJCSpTJdkcW1tbdZNpJ5kaO3XqVLGfvXv3bpV5GT16dJ5Mm9QcSSG7XAN5L7nJfcmEyVSZhmSfBg4cqNf3KdcvNDRUTaG9+eabatpNslhS9zRo0KAiM1Dbtm1T02iSuco9vSc/C13X1NfXV/tYxly/fv08105I5koyhfK8XDs7Ozs1Jn2uHZGlYnE4EWkDk6KKw7t06YI1a9Zg1qxZ6o+91Ow0btwYb731lqorulcyLZabvb19kcfv3LmjPfb444+rP/RS7yOBiUy3Sc1Pv3798kzp6XLt2jX19ZFHHtF5jgRW8p4STHl7exd4vrBjukhgIlOAchOxsbHqs6Wma8OGDWrchZHzvLy8Chwv7Jhwd3cvcEzqsnJfE5kylOm8CRMm4KOPPlJBmY2NjbqWDJyIdGPgRER6kxofuUnQtHfvXkybNk0FL1KQLBmU8iQ1TxJwSI3OlClTtMdlbJr6neJoMjjffPONzpVsEpykpaWp4Onq1asFni/smL4kwHn55ZdVTVFISIjOwEnO0wR5pfXZUtslWb25c+fmOR4fH3/P70lkCRg4EVGJSfaia9euqFy5sirKPnLkiAqc5LjQJ9tzvySQkektzWdqSEF57qJ2zXgLG5esnpPvQfo8vfDCCzo/S7JdsqpPitelKFwzXSdBxp9//lnsWCXwun37dqGZIE12p3r16jpfL9d6/fr1quBdE+xJBkxWOd7P9ct/7aQIXQrnfXx87vl9icwdAyciUiTjoVkhl1udOnVQtWpVtYLu8uXLatm+9GO6desWZs6cqaaf5A+75lxZ/i6r4GTqT+qNJCAoKii4VzItJ9OHEshIMCFZL6kFWrBggQqG8veoEtJDSeqSJPCRdgISyEi2SWqcJEsl02ayLF+W6svKNPmqycjIdJbUPcnqNpnikuBMVvFJq4TiMlySHZPxSVuEnj17qsBE6rMk0yTXUK6V1E/pItOhEqDJtZf7co2lPklqrwprZ6APWaEo35Nk7OTnd+bMGXz44YfquhT274CIchi6Op2IjHdVndy+//57dd5ff/2V1bdv36waNWpk2dvbZ3l6emb169cva8eOHXne75dffslq2LBhlp2dXZ7VWLpW1fXv37/AmOS8iRMn5jkmK8fk+BdffKE9Jqvxhg4dmlWlShW1Sq1Pnz5ZISEh6n1Hjx6d5/Vff/11lr+/v1opmH/l37Zt29Q4ZOWZjFu+R3ksq/FykxV9zZo1U9+/rJD77LPPCv2+8ktJSVGrEuX6yescHBzUyrdGjRplvf7661mxsbHFrmKT6ywr++S13t7eWa+99lrW9OnTC6x01HVNu3btqm65x/Tqq6+q71XGEhQUpFYuynWT9yhuPESWykr+RxNEERGR6ZD+U9LL6uzZs4YeCpHF4FQdEZEJmDRpkuquLtN8MjUo06GbN29WU5NEVH4YOBERmQCpqZI6M1lJJ4XdslXK0qVL8eSTTxp6aEQWhVN1RERERHpi53AiIiIiPZl94CR9VqSjsGxB0LRpU3z//feGHhIRERGZKLOfqpO6AOkkLHteyc7i0s9F9pwqrBEdERERkUUXh8veS5qNQmWPKwmkShIrSnfeK1euqKZ5UpBJRERE5kXiApmhkma9xTaUNXQjKWk8N2DAgKxq1aqpJmurV68ucM7s2bOzatWqpRq/SZO27du3l+gzbt68qZrWOTk5ZX377bclem1ERESRzQF54zXgvwH+G+C/Af4b4L8BmMU1kL/5xTF4xkm2DAgMDMTYsWMxdOjQAs8vX75cbYI5Z84cta/UvHnz0LdvX7W3lK+vrzqnZcuWajouv02bNqnoUbZfkO0TZJNM2dZAtlXQtat4fpJpEhEREWqLByIiIjIvspek9EjT/M03mRonmQpbvXo1hgwZoj3Wtm1bBAUF5dnBW/Z1knNkZ/aSeu6559CjRw+1Z1RhJADLHYRpLqbsNcXAiYiIyPzI33pXV1e9/tYb9aq61NRUHDp0SG0rkJs83r17t17vIVkmuSBCvm7fvh0NGjTQeb4EY3LxNDfuEk5EREQmETjFxMSoYu7802ryWLrn6kN2c5cd1GU6sFOnTnjhhRfQrFkznedPnTpVRZyam0zREREREQmD1zjpI/9qNpld1HeFm9Q/BQcH6/1ZDg4O6kZERERkUhknDw8P1U4gf3YpOjpa7+JuIiIiIosInOzt7VXGSHYAz00ed+jQoUw/e/bs2WoTTek6TkRERGQUU3UJCQkIDQ3VPg4LC1NTa25ubqrdwKRJkzBy5Ei0atUK7du3x/z58xEeHo4JEyaU6bgmTpyobppKeyIiIiKDB04HDx5E9+7dtY8lUBKjR4/GokWLMGzYMMTGxuLDDz9EVFSU2jJl/fr18PPzM+CoiYiIyBIZVR8nU+/tQERERKbHbPo4GRJrnIiIiCg/ZpyKwYwTERGRebvNjBMREREZm61novHcskOITSi4v6ypMHhxOBEREVmGMQsPqK8VHGzx5aOBMEWscSIiIqJyFXEjyWSvOAMnHVgcTkREVDayTPjCMnDSQZpfnjx5EgcOZKcViYiIqHSYcickBk5ERERUrKV7L+Hx7/ciISX9vq+WCcdNDJyIiIioeO+sCcHu87FYtCvsvi9XZgkip6TUdKPKUDHjRERERHqLv1MKGSc9z9sVGoOAdzdi/vYLMBYMnHRgcTgREVHpycqVNcrUM3KasuqY+jptw2mj+VEwcNKBxeFERESlJyU98+4DPafeUnO9pjRqq0oDAyciIiIqUkp6xn1foeTUjBJlnOQzbySmah8Hh9+CMWDgREREREVKSrkb9GTd47VKTssoEIiFxyZhwtJDOHM1vsD5IZFxSMu4+2n/nY42ip8SAyciIiIqUu5psju5AqB7DZwScwKxF389gr9PXMWQ2bsKnP/n0Sj11cUhe3e4tUevID0j13SfgTBwIiIiMlOh0fHqdr+SUgsGPfczVacJok5fva19LMXjcs7SPRcRFpOIP49eUc/NGNYc7hXsEZOQgvUhV2FoDJx04Ko6IiIyZRKE9JyxXd3ut0Ypd8YpsYgi7ZuJqfh1fzji76QVeC53pkrzHpUc7bTHLsUmYfKKYLzzxwl0/3IrYhNTUd3VEd0aVMWo9rXUOd/+dw6Z+i7JKyMMnHTgqjoiIjJl127f0d6/nXx/K9KkCaVGYq77+T330yFMWXUcX248U0yNU6Z6z+sJKdpj3++4gPXH82aUJnSrAzsba4zpWAsujra4mZSGcANvEMzAiYiIyAzJ1JZGYRkgjT3nYzFm4X6cv56g85xEPTNOey/cUF//yJlm0zVVJ05FxefpSvDTvvACrxnW2kd9dXWyw+Kn2mDH691Ry6MCDImBExERkRmKjs8dOOkOdj748wS2nrmOB/63Tec5CbnqmpLyBUB3j9/9jJpVnIrMOIkd566rr54uDnmOt/CtjNpVK2DFhPZwsLXRHg/yrQJHu7uPDSW7VJ2IiIjMd6quiIxT7qBKsk/t67gXGRQl6Mg4Hbscp71vY303LxMdfwdOdjYFVuN9/c859XVcJ3/sD7uBf3PaDbw7IAAtfKvAWDHjREREZCIkAHrgf1vx0V8nSy3j5GR/N4uzLyy20HNyB0tJOjJOx3MFTjE5ny19mnp8uQ1P/LCvwFSdkOLv0R1q4dOHm6KaqyN83JzQuLorjBkzTkRERCZCumefv56IuORIvDMgoMhzo2/fDZxuJ+vOOOXuzp37vq4GmIk6Mk6Rt5K196XoW9oLfP3vWRV0STaqTS23Aq/5/JFANf0mt38nd4W1lRXsbY07p8PAiYiIyERczZl+i0lIVfu4FRVkyBRZcRmnjMws3Ey6GyxJC4DiMk4p6ZmqEaWtTd7Pvhp39/NkbDLW9cezm1iKAxezC8cfblEDXRtURRVne3Sq56F93tneNEIS4w7rDIh9nIiIyNhE56pbyh0YFeZ6nqm6wjNOEjTlXtl2I0FHxilfC4LEnGm33w9dxpYz2bVJUbnGJlYdjsSdtLudvo/mTOW5V7TH4OY10KV+VZgiBk46sI8TEREZa8ZJXMs1FVd8cXjhGaf8U3O6purydwtPTElX+8u9uuIoxi48oDJM13JlnMTPOe0FrK3yvpeniyNMGQMnIiIiE5E7WModGOUngYw0iyxuVV1svgyTrqm6/E0vk1LTceLK3WLwk1G3tc0sa7k756l5+mBQ4zyvHdqyJkwZAyciIiIDkZVmzy07hOUHCjZ/LEzuYCl3TVF++afxdHUOj03MDnb8coIdmborbEuT/MXlcclpCInM3mdO/Hc6WtVL2Vhb4YFGXtrjD7WogSfb+WHm8OaoYG+DN/s1hFsFe5gy06jEIiIiMkN/n4jChpCr6iY1P9VcCzaO1BU4XSuixunyzbsr3IqqcdJMzdWtWlHtFSfBj2SnKjvfDW4kkDp3LUHbPuBK3B0cCb+FkMi7GafNJ6+pr1UrOuCtfo3QKyA7eGrj7wYrKytV0zQosLq6b+qYcSIiIjKQyzfuBjhztpwv8lwJanIXfOevKcotIt9+brpqnGR1nvB2dYSLg612uu7ctXhM23BKFaPL3nDxKelqBd+T7f3UOdvOXs8zVXcqKjv75OXqCGtrK7St7a5uuQMlcwiaBAMnIiKiUnQzMbXYFW8aF2IStfeXH4go0F07/95zuWfRcheK5xeRk3Fq6O1STMYpOxBzr+gAt4r22izU9L9PY962C3jwq+3YeyFW+15d6mWvhNtxLkatrJNgKndLhAZeFWHuGDgRERGVEqk76vG/rXhwxvYiN8PVyL2xbmpGZp6mlfnlLwYv6tzLORknTRduXQ0wL8YkafeL09QeSVbrn1PR2lqmKauO57xXJTSqVgnuuWqUPnu4KV7uWU/dl4zVpAcbwNyxxomIiKgUyFTaS78e0a5mC4tJRJMaurcPkc7aF67fzTiJmMQU+OYUauenKQZ3dbJTAY1knOQ9CpsCi7ipCZwqYeXh7AaWUqsk02gakt3an9OUsl1tN2zN6ce0+3xMoZ8vQZiNtRV+GN1K9WjyqOigir/TpSjcygqt/d3UlJ+5Y+BERERUCjafvIp9YdmBiJDaoKICJ8nsSEAjsUwD70qqTih/e4DCpt+a+1RWNUayZ5y83sXRrsC58tmawEnIFJ+0FJBzYxNScCoqHhlZWaptgewRV6dqRfi4ZQdsy/Zmr/DrVNcDF64nqGJwyUYNbFZdHZcNeHNvwmtnY4Vnu9aBpeBUnQ7sHE5ERCWx98LdoKmwAu38ZM85IQGLrFbT1DHpcjqnAFsCJ1nan31+wUBLMkmafk91PSvCIacGSYIyyVCNW3wQTy7Yh9E/7lfHO9fzUFmr0e1r5Xkf2Q5F9pIb3Lw6/pjYEa7OBQM0S8TASQd2DiciopI4En5TffXNydxosj7F1TfV9qigtiERkg3S5fTVePW1UTUXeFZyLLAFS/5WBM72NipTJMGTOHMtXtUuBUfcynN+v6bV1NdaHhUwvrO/uv9ggBdGtfdTwdPM4S202Shi4ERERFRo1mb8koOYtv6U3uefuJKdEZIMTe6pNV2k27aQwEbqhXRlkIRsqiuBj5ACbemXJDTdunPbee66dhWcZJIaemdP18lU4Pztd1seSCbqoyFN0K2Bp/bYm/0aYfeUHpg/sqXJbLpb3nhViIiI8llxMELb1HFC1zqoUky36+ORcapIWgKg9nXc8c1/ocVO1e08l12E3dbfHZdyztU1VXcxNlHVI0kWyaeKM6q65AROufo6afx5LEp97Z9TkyQZKrHnfCwOXsrOiu2Z2gPuFRzytBIQEmhVr1x0E05Lx6k6IiKifNmdedsvaB9r+hgV5XBOQBLkW1k7VXf5ZnYn7sKExyapqTxbayu0q+MOD+1UXeEZJynmFg28XdTKOF2Bk3TzPnTpJmSh3YBm1bQZKiGF61lZ2Zko6VCeP2gi/fCqERER5bLueFSeLUt2ny8+cNqfs5ouyK+KCkokIErLyNK5Ee+O0OzptBa+lVHRwVY7VafZOy4/TV2SZtotf+C070Isun6xBUNm71KPu9avCq+cOihNE0wNeY7uHQMnIiKiXFYejlRfm9V0LbKvkYYUdEt7ANGtQVXV66hGFaciC8S3nck+v3NOJ25NcbiuVXKrDl/Wvn+ewClnak8yZLLXnEwXSmD0v0cDta+XruC1cvWG6pmzjxzdG9Y4ERER5UhOzdBOzb07IACPztuj2gZI5kiTwclv7dErKmBpWsNVmxGS6ToJZCRwalfbPc/5ssfb5lPZ9VMPNMouzNZknG4mpaqpQlubu3mNNUciVVPNGpWd0LORV4GMk2zKuyOnIHzR2NYqcMrfFHPh2DbqHE8XR7Su5caf931gxomIiCjHngsxqghbgpSWflXQJGfLEims1uX3Q9nZoEda1tQe8/eooL6ezWkhoCE1Tx/9dVLVGkkNkmZLlCrO9qouSY7fSLqbddpwPArv/3lC3R/dwU9ls4RmVV20bI9y8pqaFqznWVGtkCusk7iMZ1T7WujTxJs/6/vEwImIiMySTLG9syZEZWT0tTVnCq1rg+ysTRv/7OyMFFwX5s+jV1QbAnsbawwKzF7FJjQBl6y205AtT17//ZhqlCmtAKb0bah9TgIiN+e8BeIyPTfx58O4k5apskgS+GjI3nLZ56bgr5xVdH1z+jFR2WLgREREZkeW9T+37DCW7r2E2VtC9XqNdNXWBE7dcgqoJeskDuc0t8xNpu/e+SNE3Z/QLW/LgqY59VESVEnAJKRFwcrDl1WQ9NWw5qhZJW9TSY9cWSRpZfDa78fUVinDW/vgxzGt4WiX3S1cU7ckySd5/r/T2XvM9WvKbFJ5YOBERERm59N1p9RGuOKnveGIy9l4t7iGlFKTJHuvdajroY4F5ezJJs0jE1PStefKHnHjFh/AraQ0NKlRCf/Xo26e95JpM8kqyXnSg0lWvc3896x6btrDTbXdunOr65Xd4Ts4/BZWHLqspvXa+rvh04eaaqfoNORx7m7e0n28gVfe1XNUNhg4ERGRWTlw8QZWHYlUNUOyga0EL0v2XCz2dYt2ZZ/Tq7G3ahEgvF0d1T5yktk5ejm7JYAENBN/OoyQyNtqS5NvRgTBLlcxt5Dibk3/JHndW2tC1HsMDaqJx1r5FPr5HepkF5HvDL2OlTl1U4+39VV9mwrzeu+7U30y5sJqm6j0MXAiIiKzsmTPJfV1WCsfbR3Rot0X1Wo1XaRW6I+jV9T9pzrm3ey2Rc503ZHw7MDpm//OqfYDjnbWWDimtbYQPD9ZZSembziD0OgEVHG2w7sDA3SOoWOd7CzXgYs3EXkrGS6OtujdWPf0W/9m1fBEW1+4V7DHiDaFB2NU+hg46TB79mwEBASgdevWZXDZiYioLNxMTMXGkKvq/pPt/NC/aTW4OtkhNjEVRy/fLdQuLNiS1XTSu0kzPaeheSxNLrecjsbMf8+pxzKFFuhTWed7agKnqzlNMP+vRz01Fl383J3VNKGGjD93XVNhPnmoKQ6+3RN+7oUHb1T6GDjpMHHiRJw8eRIHDhwog8tORERlYfWRSKRmZKJx9UpoUsNVTZl1qpedydl2JruIOr+TV25j7tbszW+f7ly7wJRXu9rZK+sky/TM0oOqZcCINr54OOhu+4HCSLaoS/2qqvlkz0aeeKKdb5Hny+cOCqyh7ku26eWe9fT6njlFV77YAJOIiIyKFHKvOBQBz0qOeZb467Mq7reDEeq+rETTkBVy645FqcBnUq8GeV6TlJqOV5YHq2BLmksOzNnfLTfptfRqr/r4ctNZ1S9JapE+HNy42PG4OtthyVNtUBKvPFgPbhXsMLpDLTjYFp1tIsNg4EREREZjwc4wfLX5rCro1jR6bJ9TNF2cLWeicfpqvKo9GtQ8O3OTe2+2Y5FxqpZJlvILmZqbsOwwzlyLV60Apg9tqjN780KPeiqA2hsWi+e71S1QDF5apEXBW/1110GR4XGqjoiIjIJsXSJdtSVokqkq8ebq42qvNn2yTTM2Zy/3H92+Vp5aIslcBVSrpKbYtudsTSK9lV5dcRTbz16Hk50N5o9qqQ2odOne0BNT+zYqsk6JzB8DJyIiMrgzV+Pxxu/H1P1nu9TGzjd6wKuSA8JiEjErpxi7KJtOXlPtAZztbfBMl9oFntdsjjtny3kVmH3w5wkVqNlaW2Huk0EFCsKJdGHgREREBiWNKp9dehDJaRnoXM8Dr/dpqLI6Hwxqom0lUNS2KfJ6yVSJMR1qFZo5GtvRX21Tci46AR0/+w+Lc1oW/O+xQLW/G5G+GDgREZHByJTZ5N+CcTE2SW2sO3N4C22X7N6NvVDXsyKSUjOw5kikzim6N1cdx+WbyahZxUltfVKYqi4O+PbxIPXeEmhJg8vPH2mGwblqoYj0wcCJiIgMQoKeLzadwT+nomFva415I1uqTtwaUqgtDR7Fsr2X1Pn5zd9+AeuOR6kpt29GtEAlR931R7Jh77wnW+KVnvWx9bVuOjt4ExWFgRMREZU7CYI+33hG2z/pkyFNVN+l/KRXkhRvn72WgIOX8m60+/O+cEzbcFrdn9qvEVroUafUM8ALL/Wsp91Ql6ikGDgREdF9uXwzSU2lzd4SigvXE4o9X9oAvLk6RBs0vTcwAI/qyP5IrZOml9PH604hOTVDBV3fbTuvVtyJ57rVwbhO/vwpUrmwyios90lat2/fhqurK+Li4lCpUvaGjURElB0wvb0mBFvPZC/xF1KeJIHOG30bopqrU4HLdO32Hbzw82G1H5u0TPpwcBOMbOdX5OW8FJuIwbN34VZSmuribW9ro9oICAmY3u7fiN2zqdz+1jNwKsWLSURkCTIys7Bkz0V8sfGMKtyWgmuZZnNxsMXO0Bh1jtQqfT2sudpyRFMEvvLwZbX67fad7D5Ns4a3UL2R9HHw4g088cM+pKRnb9QrQde7AwLUajmi+8XAqRQxcCIiyttvacqqYzgSfks9blPLDdOGNkWdqhXV45DIOLyx8hhOXLmtghvpj1TPs6LKEF2Jy97sVjbSlaCqds5r9CVbpkiReHOfymp7lAbeLvzRUKlg4FSKGDgREQG3klLx9T/nsHTvJZVxkuX8U/o2xONtfGGd0z5AQzp9S2bpp33heY5LRur57nUxvrO/2nyXyFgwcDLQxSQiMjfSZXvx7otq2b/0P9L0V5LmlN6ujkW+NuJGEvZeiMXZa/FqxVuPhp5wtOPGtWTaf+u5yS8RERVwPT4FS/dcxJK9l1RRtqjvVRHvDmiMTvU89LpiPm7O6kZkThg4ERGRtoB7z4VY/HYwAhtCrqq2AaK2RwW8+EA9DAysru3qTWSpLCZwSkpKQqNGjfDoo4/iyy+/NPRwiIiMJljaF3YDG0KisPHEVVy7naJ9Toqwx3eurabmWJNEZGGB0yeffIK2bdsaehhEREZT7L10zyX8digCETeStccrOdpiUPPqeLSlDwJ9Kht0jETGyCICp3PnzuH06dMYOHAgQkJCDD0cIiKDkZ7HfwRfUaveYhNTtavd+jb1Rt8m1dChrjscbFnATaSLwdeDbt++XQU01atXV51f16xZU+CcOXPmwN/fH46OjmjZsiV27NhRos949dVXMW3atFIcNRGR6ZFVceMWH8TLy4NV0CT9lWY8Foj9b/XE548EqmaUDJqIjDzjlJiYiMDAQIwdOxZDhw4t8Pzy5cvx8ssvq+CpY8eOmDdvHvr27YuTJ0/C1zd712wJplJS7s7La2zatAkHDhxA/fr11W337t3l8j0RERkb2bbkqUUHcP56IuxtrfFij7p4pksddZ+I9GdUW65Ixmn16tUYMmSI9pjUJQUFBWHu3LnaY1LkLefok0WaOnUqli1bBhsbGyQkJCAtLQ2TJ0/Gu+++W+j5EoDlDsKkt4OPjw/7OBGRydofdgPPLj2Im0lp8K7kiB9Gt1JbpBBRyfs4GfV/aqSmpuLQoUPo1atXnuPyWN/skQRXERERuHjxolpNN378eJ1Bk+Z8uXiamwRNRESmasXBCDzxw14VNMlWJ3+80JFBE9F9MOrAKSYmBhkZGfDy8spzXB5fvXq1TD5TMlQScWpuEnQREZka6cH0ybqTeO33Y0jLyEK/pt5Y/kx7eFUquts3ERl5jZO+U3i5yexi/mP6GDNmTLHnODg4qBsRkamSLU5e+jUYp6Juq8f/16MuXulZv8CeckRkZoGTh4eHqk3Kn12Kjo4ukIUiIrJ0smruu23nsWBnmMo4VXG2w7SHm6FPE29DD43IbBh14GRvb69WzG3evBkPPfSQ9rg8Hjx4cJl+9uzZs9VNpgqJiIxZbEIKfj0Qge93XNDuK9etQVV8PrQZPDk1R2RegZOsdAsNDdU+DgsLQ3BwMNzc3FS7gUmTJmHkyJFo1aoV2rdvj/nz5yM8PBwTJkwo03FNnDhR3TSV9kRExkRKFg5cvIlley/hb9lXLiN7XznpzfRGn4Z4oJHnPZU0EJGRB04HDx5E9+7dtY8lUBKjR4/GokWLMGzYMMTGxuLDDz9EVFQUmjRpgvXr18PPz8+AoyYiMoz4O2lYfSQSP+0Nx5lr8drjgTVdMap9LQxuXp37yhFZSh8nU+/tQERUVs5cjcei3RfxR3AkklKzSwgc7awxOLAGnmznh6Y1mRknKo+/9QbPOBkr1jgRkTE4GnEL324JxeaT17TH6npWxJNtffFQUE24OtkZdHxEloYZp2Iw40REhhB9+w4+/Osk/joWpR5LuVKfxt4Y3aEW2vq7sX6JqBQx40REZKKkeuKnfeGY/vdpxN9Jh7ReGtKiBp7vVldlmojIsDhVR0RkJG7fScPk345qp+Wa1nDFtIebcosUIiPCwEkH1jgRUXmSLt/PLTuEi7FJsLexxht9G2JMh1qwYbdvIqPCGqdisMaJiMra3yFReGX5USSnZaBGZSfMfTIIzWpW5oUnKiescSIiMpF6pjlbz+OLjWfU4871PDBreAtUqWBv6KERkQ6cqiMiMoA7aRl4c9VxrDoSqR7LtNzb/RuxeSWRkWPgRERUzi7GJOK5nw6ruiapYXp/YABGtq/FnwORCWDgRERUjmRfudd+P6paDbhXsMesES3Qsa4HfwZEJoKBkw5cVUdEpSkqLhkfrD2Jv09cVY9b+lXB7MeD4O3qyAtNZEK4qq4YXFVHRPcjKTUdS/dcwqx/zyExNUNNzY3vXBuTe9WHnY01Ly6REeCqOiIiA4tLSsOSPRfx464w3ExKU8eCfCvj04eboqE3NwwnMlWcqiMiKsX2Akcvx+GXfeFYe/SK6ssk/NydMbF7XTwSVBPWbGhJZNIYOBER3af4O2n4I/gKft4XjpNRt7XHG3q74LluddC/aTW2GSAyEwyciIju0fnrCfhhRxj+CI5EUmp2dsne1loFSo+39UUrvyqwsrLi9SWy1MApLi4Oq1evxo4dO3Dx4kUkJSWhatWqaNGiBXr37o0OHTrAXHBVHRHpEhIZh7lbz2N9SBSysrKP1a5aAY+38cXQoJrs/E1k6avqoqKi8O677+Knn36Ct7c32rRpgxo1asDJyQk3btxASEgIDh06BD8/P7z33nsYNmwYzAVX1RGRRsSNJEzbcArrj2e3FBA9G3nh6c7+aOvvxuwSkYkq9VV1gYGBGDVqFPbv348mTZoUek5ycjLWrFmDGTNmICIiAq+++uq9jZ6IyMgkpqSrDNP8HReQmp4Jqe8eGFhd1S9xhRyRZdEr43T9+nU1Jaevkp5vzJhxIrJcmZlZWBMciel/n8a12ynqWPva7nh3YAAaVWNLASJzUeoZp5IGQeYSNBGR5ToSfhMf/HkSwRG31GNfN2e82a8Rejf24pQckQUrcdvaxYsXY926ddrHr7/+OipXrqwKwy9dulTa4yMiKlfHLt/CuEUH8NCc3SpoqmBvgzf6NMTmSV3Qp4k3gyYiC1fiwOnTTz9VReFiz549+Pbbb/H555/Dw8MDr7zySlmMkYioTKVlZGLb2esYu3A/Bn27C/+ejlZ1TI+2rIktr3VTtUwOtjb8KRBRyfs4SeF33bp11X0pBn/kkUfwzDPPoGPHjujWrRsvKREZPSntvBJ3R7UV2HzymrrFJWdviyJ7yQ1uXh0vdK+L2lUrGnqoRGTqgVPFihURGxsLX19fbNq0SZtlcnR0VCvrzAX7OBGZJln1diMxFTEJKeoWm5B9PzbnWHhsEs5ci0f8nfQ8r/OoaI++TaphXCd/1PKoYLDxE5GZBU4PPvggnn76adX08uzZs+jfv786fuLECdSqVQvmYuLEieqmqbQnIuORkZmFC9cTEHIlDiev3EbEjWRExSUj8tYdFRzpw9baSjWt7FDHA32beKNVLTeVbSIiKtXASTIxb7/9tpqyW7lyJdzd3dVxaYA5YsSIkr4dEZHe25tsOnENW05H43hknHYD3cJIAORewR7uFR1UJsmjooN67OHigGqujmjg7YLaHhXV9ihERKXex8mSsY8TkWE7dS8/EIG/T1xFaHRCnuec7W0QUK0SGlevBH+PCqhW2Qk1KjupwKiKsz2smT0iIkP1ccrv5s2bWLBgAU6dOqWW5jZs2BBPPfUU3Nzc7uXtiIjyOBx+Ewt2hGFDSBQys+5OrbWv447ejb3Rrra7CpY4tUZERp9x2rZtGwYNGqQis1atWmmn6W7duoW1a9eia9euMCfMOBGVH2kJMOvfczh06ab2WOd6HnikZU10a+AJVyc7/jiIyKB/60scOMleddLscu7cubCxye5rkpGRgeeffx67du1SG/6aEwZORGXvVNRtfLr+FHaci1GP7WykJUANtcKNW5sQkUkHTtL8Mjg4GA0aNMhz/MyZM2jevLlZtSQQDJyIys6123fwv01nsOLQZchvIgmYRrWvhWe71IZnJUdeeiIy/RqnoKAgVduUP3CSYxI4EREVJyk1HfO2XcD87Re0q+P6N62G1/s0gJ87eygRkfEqceD04osv4qWXXkJoaCjatWunju3du1e1Kfjss89w7Ngx7bnNmjUr3dESkcn3X/r9UAT+t+ksouOz+y0F+VbGW/0D0NKviqGHR0RUrBJP1VlbF933RFbZyVvKV6l9MofO4dLoU5/0HREVLjMzC5tPXcNXm8/i9NV4dczXzRlT+jZUzSfl9wURkVnWOF26dEnvc/38/GDqWONEdH8ZJmkp8O1/odqASVbG/V+PuhjZ3o8b5xKR+dc4mUMwRERlKzk1A+uOR2Hu1lCcv56ojlV0sMXoDn4Y37k2Kjvb80dARCbpnhpgRkZGqtYD0dHRyMzMLFADRUSWR5LX0n/p90OX8dexKCSkZG+iW8nRFk918sfYDv5wdWYfJiKysMBp4cKFmDBhAuzt7dU+dblrE+Q+Ayciy3EnLUN1+d4VGoN1x6JwMTZJ+5yPmxOGt/bFqPZ+cHFkwERE5qHENU4+Pj4qcJo6dWqxheLmgDVORHlrlk5euY2doTHYfT4G+8NuICU9M8/+cf2aVlOdvtvUcuN+cURkEsq0xikpKQnDhw+3iKCJyNLJf1ddik3SBkq7z8fiVlJannOqujigYx13dKlfVe0jV8HhnioAiIhMQol/w40bNw4rVqzAlClTymZERGRQ1+NTVJAk02+7QmMReSvvbgBS5N2uths61PFAp3oeqOdZke0EiMhilHiqTvoaDRgwQG2t0rRpU9jZ5a1dmDFjBswJp+rIEqbfpE5p88lr2H72urZtgIZsg9LCtwo61fVAx7ruaFazMuxsmHEmIvNRplN1n376KTZu3KjdciV/cTgRmUZR985zMdh08ir+PRWN2MTUPM8HVKukgqSOdT3Qxt8NzvacfiMiEiX+bSgZpR9//BFjxozhFSQyIZJclszSr/sjVI+lpNS7nf2lZUCPhp7o0chL1Su5V3Qw6FiJiMwmcHJwcEDHjh3LZjREVOpiE1Kw+kgkfj0QgdDoBO3xGpWd8GCAF3oFeKG1vxun34iIyiJwkg1+v/nmG8yaNaukLyWicnT2WjxmbwnF+uNRSMvILmV0srPBgGbVMLyND4J8q3B6nYiorAOn/fv347///sNff/2Fxo0bFygOX7VqVUnfkohK0Zmr8Zj17zmsD4mCZulHs5quGNbaB4MCq7MZJRFReQZOlStXxsMPPwxzN3v2bHWTVYREpuBU1G0VMG0Iuao91reJNyZ2r4smNVwNOjYiIottR2Bp2I6AjN2123fw8bpT+PPoFe2x/k2r4f8eqIuG3kUvqyUiIpRtO4LC3Lx5E8uWLcOCBQsQHBzMnwFROUjPyMSSPZcwY/NZtaGudAOR7U5e7FEPDbxd+DMgIioD9xU4/fPPPypYWrNmDTw8PCxiCo/IGEhbgbdXh+Bk1G31uLlPZXw8pAmn5IiIjC1wCg8Px8KFC9UtISFBZZt+++03DB06tGxGSERat5JSMf3vM/j1QLgq/HZ1ssMbfRpieGsfbqhLRFQO9N43QYKjXr16oVGjRggJCcHMmTNx5coVtdmvHCOisiXboTz41Xb8sj87aHqkZU38O7krHm/ry6CJiMjYMk6PP/44Xn/9daxcuRIuLqyfICovKekZ+HLjGXy/I0w9rlO1AqY93ExthUJEREYaOD311FOYM2cOtm3bhpEjR2LYsGGoUqVK2Y6OyMJJp++Xfj2CE1eya5mebOeLt/oFwMnextBDIyKySHpP1c2fPx9RUVF45pln8Msvv6BatWoYPHiw2v8qMzOzbEdJZGHk/1cyJTfgmx0qaKribIfvR7XCx0OaMmgiIjKFwEk4OTlh9OjRKut0/PhxBAQEwMvLS+1dJ1N57BpOdP9u30nDc8sOY+qq47iTlolOdT3w98td1L5yRERk4g0wJdu0bt061ZZgw4YNSElJgTlhA0wq7/3lnl16CGExibCzscJrvRvg6U61WfxNRGQkf+tLtXN4dHQ0PD09YU4YOFF5WXcsCq/9fhRJqRmoUdkJc54IQqBPZf4AiIiM6G+9XlN1e/bs0euDJWhKTEzEiRMn9BspEakO4NPWn8LEnw+roKljXXesfaEjgyYiIiOkV+A0atQoPPjgg6qXkzS9LMzJkyfx5ptvom7dujh8+HBpj5PILMUmpGDUj/sxb/sF9fjZrrWxeGwbuFd0MPTQiIjoXtsRSFA0b948vPvuu3jiiSdQv359VK9eHY6Ojqpz+OnTp1WmSbZc2bx5M5o0aaLP2xJZtOOX4zBh2SFE3kqGs70NvngkEP2bVTP0sIiIqDRrnCSbtGPHDly8eBHJyclqj7oWLVqge/fucHMzv4Z8rHGisvDXsSt4dcVRtWrO36MC5o1sifpebCxLRGTsf+tLvFddUFCQuhFRycl/p3z9zznM/PecetytQVXMGtEClRzteDmJiExAiQMnIro3yakZePX3o2r1nHi6kz+m9msEG2srXlIiIhNhEYGTra2ttu6qVatW+OGHHww9JLIwV+PuYPySgzgeGaf6M308pAmGtfY19LCIiKiELCJwqly5MoKDgw09DLJQRyNuqaApOj5FbZ3y3ZMt0ba2u6GHRURE98AiAiciQ/nzaHYReEp6Jup7VcSC0a3h4+bMHwgRkSXsVVcWtm/fjoEDB6r2BlZWVlizZk2Bc+bMmQN/f3/V/qBly5ZqVV9Jq+XldZ06dVL77BGVtczMLMzYfBb/98sRFTT1aOiJlc91YNBERGQJGadZs2bp/YYvvvhiiQYg/Z8CAwMxduxYDB06tMDzy5cvx8svv6yCJ9lMWPpJ9e3bV/WW8vXNrhGRoKiwPfI2bdqkAjJpnSBfQ0JC0L9/f7VBcXHLDYnupwh88opgrD9+VT0e39kfU/qyCJyIyGL6OEm2R683s7LChQsX7n0wVlZYvXo1hgwZoj3Wtm1b1f5g7ty52mONGjVS50ybNq3EnyFB10cffaSKxAsjAVjuIEyyVT4+Pnr1diCKiktW9UwhkbdVEfgnDzXFY618eGGIiCypj1NYWBgMITU1FYcOHcKUKVPyHO/Vqxd2796t13tIZ3NnZ2c4ODjg8uXLKlNVu3ZtnedLMPbBBx/c99jJ8gRH3MIzOUXgbhXsVRF4G3/zawpLRGTJjLo4PCYmBhkZGfDy8spzXB5fvZo9DVKcU6dO4dlnn4W1tbXKaM2cObPIDudTp07FpEmTCmSciIqy9ugVvJZTBN7AywU/jG7FeiYiIjN0T4GTZG7Wrl2L8PBwlRXKbcaMGShtEvDkJrOL+Y/p0qFDB1XTpC/JTMmNSB/yb/Grf85hVk4ncCkCnzm8OVzYCZyIyCyVOHD6999/MWjQIFX3dObMGdVYUoqv5Q9IaW/FIvvg2djYFMguRUdHF8hCEZW3O2kZeP33YyrbJJ7pUhtv9GnITuBERGasxO0IZCpr8uTJaoWatAdYuXIlIiIi0LVrVzz66KOlOjh7e3u1Ym7z5s15jstjySSVpdmzZyMgIACtW7cu088h0xSbkIInftingiZbaytMH9oUb3L7FCIis6fXqrrcXFxcVBfuOnXqoEqVKti5cycaN26Mo0ePYvDgwSr7VBIJCQkIDQ1V91u0aKGm+rp3767qkKTdgLQjGDlyJL777ju0b98e8+fPx/fff48TJ07Az88PxlRpT5YhNDoeYxcdQMSNZFRytFVF4B3qehh6WEREZCyr6nKrUKGCdrm+9EY6f/68Cpw0xdwldfDgQRUoaWgKs0ePHo1FixZh2LBhiI2NxYcffoioqCg1Nbh+/fpyCZqI8tsVGoMJyw4h/k46fN2c8eOY1qjrWZEXiojIQpQ4cGrXrh127dqlprGkmaRM20nx9apVq9RzJdWtWzdVH1WU559/Xt2IDGn5gXC8tToE6ZlZaOlXBfNHtoR7RS4kICKyJCUOnGQqTabXxPvvv6/uy3Ra3bp18dVXX8FcSI2T3KQdAlk22T7l841n8N228+rxoMDq+PyRZnC0szH00IiIyNhrnCwNa5wsm2yfMum3YGwIyV7Z+eID9fBKz3p6t8MgIiILr3HSkP5N0hYgMzMzz3HN/nFEpi46/g7GLz6Io5fjYG9jjemPNMVDLWoaelhERGRAJQ6czp49i3HjxhXY8kTTlJJTW2QOzlyNx1OLDiDyVjKqONth3shW3D6FiIhKHjiNHTsWtra2+Ouvv1CtWjWznbJgjZPl2nb2Oib+dBgJKemo7VFBrZyr5VHB0MMiIiJTrHGSdgSy8W7Dhg1hCVjjZFmW7r2E99eeQEZmFtrVdlM9mio72xt6WEREZKo1TtKG4F76NREZMwmUPl1/Cgt2hqnHQ4NqYtrDTWFvW+Lm+kREZMZK/Fdh+vTpeP3117F161bVmFKitNw3IlOTlJqOZ5ce0gZNr/VugC8fbcagiYiI7n+qzto6O9bKX9tkrsXhnKozbzcSU1UReHDELRUozXgsEAOaVTf0sIiIyFym6rZs2XI/YyMyGhE3kjD6x/24EJOIys52WDC6FVr6uRl6WEREZMRKHDh17doVloCr6szbySu3MXrhflyPT0GNyk5Y/JTsOedi6GEREZG5TdUdO3as8DeysoKjo6NqgOngYD77d3GqzvzsPh+DZ5ccQnxKOhp6u2DR2DbwdnU09LCIiMgcp+qaN29eZO8mOzs7DBs2DPPmzVOBFJEx+evYFUxafhSpGZlo6++G+aNawdXJztDDIiIic11Vt3r1atSrVw/z589HcHAwjhw5ou43aNAAP//8MxYsWID//vsPb7/9dtmMmOgeLdoVhv/75YgKmvo28cbip9owaCIiohIpccbpk08+wcyZM9G7d2/tsWbNmqFmzZp45513sH//ftUkc/Lkyfjyyy9L+vZEpU5moz/feAZzt55Xj0e198N7AxvDxto8u94TEZERBU7Hjx+Hn59fgeNyTJ7TTOdFRUWVzgiJ7rOx5ZurjmP5wQhtj6bnu9Ux262CiIjIyKbqZKuVzz77DKmpqdpjaWlp6phmG5bIyEh4eXnB1FfVSZf01q1bG3oodI/SMjLx0q9HVNAkyaXpQ5tiYve6DJqIiKj8VtXt3r0bgwYNUo0wZYpO/stdVtpJ40vZ+Lddu3ZYunQprl69itdeew2mjqvqTNOdtAy88PNh/HMqGnY2Vpg5vAX6Na1m6GEREZGJ/60vceAkEhISsGzZMpw9e1bVj0im6fHHH4eLi/n1wWHgZJpbqIxfchC7QmPhYGutNurt3tDT0MMiIiJLbEcgKlasiAkTJtzr+IjKzO07aXhq4QEcvHQTzvY2+GF0K3So48ErTkREpUKvwGnt2rXo27ev6tEk94si03hEhtp3btSP+xASeRuVHG2x6Kk2CPKtwh8GERGVGr2m6qSeSWqWPD09tZv8Fvpm3OSXDCT69h088cM+nItOgHsFeywZ1waNq7vy50FEROU/VZeZmVnofSJjcPlmkgqaLsUmwauSA356uh3qelY09LCIiMgM3VONE5GxuBiTiMe/34srcXfg4+aEn8a1g6+7s6GHRURElt7Had++fdiwYUOeY0uWLIG/v7+awnvmmWeQkpICc8E+TsbvUmwiRuQETbWrVsBvz7Zn0ERERMYROL3//vuqX5OGdAkfN24cevbsiSlTpuDPP//EtGnTYC4mTpyIkydP4sCBA4YeChUi4kYSRszfi6i4O6hTtQKWP9Me1VydeK2IiMg4AifZ0PeBBx7QPv7111/Rtm1bfP/995g0aRJmzZqF3377razGSZSnpmn4/JxMk0cF/DK+Haq6OPAKERGR8QRON2/ezLONyrZt29CnTx/tY9maJCIiez8worISeStZTc/JV38Jmp5pB89KjrzgRERkXIGTBE1hYWHqvuxTd/jwYbRv3177fHx8vOrzRFRWouKS1fRcxI1k+Lk7q0yTF4MmIiIyxsBJsktSy7Rjxw5MnToVzs7O6Ny5s/Z5qX+qU6dOWY2TLNy123dU0BR+I0mtnpOgyduVmSYiIjLSdgQff/wxHn74YXTt2lVtubJ48WLY29trn//xxx/Rq1evshonWXhzSwmaLsYmoWaV7KCpemUWghMRUfkr8Sa/0lVTAicbG5s8x2/cuKGO5w6mzAE3+TWs6PjsoOn89UTUqOyEX59pBx839mkiIiIT2eRX3rgwbm5uJX0roiLFJKTgie/3qaCpmqujyjQxaCIiIpOocSIqT3HJaXgyZ+8570rZQRM7ghMRkaExcNKBncMNJzk1A+MWHcDpq/HwqOiAn8e3RS2PCgYcERER0T3WOFka1jiVr9T0TDyz9CC2nrkOF0dbtY1Ko2pFzzcTERGV1996ZpzIaGRkZmHyiqMqaHK0s8bCMa0ZNBERkVFh4ERGQRKf760NwZ9Hr8DOxgrfPdkSrWpxwQERERkXBk5kFGZsPotle8NhZQXMeKw5ujXwNPSQiIiICmDgRAb3w44L+Oa/UHX/o8FNMDCwuqGHREREVCgGTmRQKw5G4ON1p9T913o3wJPt/PgTISIio8XAiQxm44mreGPlMXV/fGd/PN+Nex0SEZFxY+BEBrE7NAb/9/MRZGYBj7asiTf7NYKVFDgREREZMQZOVO6OXb6F8UsOIjUjE70be2Haw00ZNBERkUlg4ETlKjw2CWMXHkBiagY61nXHzOEtYGvDf4ZERGQa+BeLys3NxFSMWbgfsYmpaFy9EuaNbAVHOxv+BIiIyGQwcKJycSctA08vOYgLMYmoUdlJdQWv6GDLq09ERCaFgROVy1YqrywPxqFLN1HJ0RaLxraGZyVHXnkiIjI5DJx0mD17NgICAtC6devy/YmYoU/Xn8KGkKuwt7HG/FGtUM/LxdBDIiIiuidWWbJJGJXKjslU0I87w/DhXyfV/VkjWmAQu4ITEZEJ/61nxonKzIbjUfhoXXbQNKVvQwZNRERk8hg4UZk4dOkGXl4eDMlnPtnOF892qc0rTUREJo+BE5W6C9cT8PTig0hJz0TPRp54f2BjNrgkIiKzwMCJSlVMQgrGLDyAm0lpCKzpquqa2OCSiIjMBQMnKjVJqekYt/ggwm8kwcfNCQvGtIazPXs1ERGR+WDgRKXWq+nFX4JxNOIWKjvbYdHYNvCo6MCrS0REZoWBE5Var6Z/Tl2Dva01fhjVCnWqVuSVJSIis8PAie7bL/vDsWBnmLr/1WPN0aqWG68qERGZJQZOdF92n4/BO2tC1P1JD9ZH/2bVeEWJiMhsMXCiexYWk4jnlh1GemaWam75fz3q8moSEZFZY+BE9yQuKQ3jFh9AXHIamvtUxuePNGOvJiIiMnsMnKjE0jIyMfHnw7hwPRHVXR0xf1RLONrZ8EoSEZHZY+BEJfbRXyexMzQGzvY2+GF0a3i6OPIqEhGRRWDgRCWyZM9FLNlzCVZWwFfDmiOgetG7SBMREZkTBk6kt+1nr+ODP0+q+6/3bojejb159YiIyKJYROAUFhaG7t27IyAgAE2bNkViYqKhh2RyQqMTVF2TdAgfGlQTE7rWNvSQiIiIyp1FbCQ2ZswYfPzxx+jcuTNu3LgBBwduBVISNxNT1Qq6+DvpaF2rCj59uAlX0BERkUUy+8DpxIkTsLOzU0GTcHNjV+uSSE3PxIRlh3ApNgk1qzjhuydbwsGWK+iIiMgyGXyqbvv27Rg4cCCqV6+ushhr1qwpcM6cOXPg7+8PR0dHtGzZEjt27ND7/c+dO4eKFSti0KBBCAoKwqefflrK34H5ysrKwntrQ7Av7AYqOthiwejWcOfGvUREZMEMnnGSeqPAwECMHTsWQ4cOLfD88uXL8fLLL6vgqWPHjpg3bx769u2LkydPwtfXV50jwVRKSkqB127atAlpaWkq0AoODoanpyf69OmD1q1b48EHHyyX78+U/bjrIn7ZHwFrK+CbES3QwNvF0EMiIiKy7MBJgiC56TJjxgyMGzcOTz/9tHr89ddfY+PGjZg7dy6mTZumjh06dEjn62vWrKkCJR8fH/W4X79+Kohi4FS0nedi8Mm67BV0b/ZrhO4NPe/hp0tERGReDD5VV5TU1FQVFPXq1SvPcXm8e/duvd5DgqZr167h5s2byMzMVFODjRo10nm+ZK5u376d52ZpwmOT8MIvh5GZBTzSsibGdfI39JCIiIiMglEHTjExMcjIyICXl1ee4/L46tWrer2Hra2tqmvq0qULmjVrhnr16mHAgAE6z5cslqurq/amyVRZisSUdDyz9CBuJaUh0KcyPh7CFXRERERGM1WnDykaz1+0nP/Y/UwH5jZ16lRMmjRJ+1gyTpYSPMl1fe33ozh9NR5VXRww70nuQUdERGQygZOHhwdsbGwKZJeio6MLZKFKi/R4stQ+T3O2nsf641dhZ2OF754Mgrcr96AjIiIymak6e3t7tWJu8+bNeY7L4w4dOpTpZ8+ePVt1GpcaKUvw3+lr+HLTGXX/w8FN0NKP/a6IiIiMLuOUkJCA0NDQPNujyKo3aVQp7QZk2mzkyJFo1aoV2rdvj/nz5yM8PBwTJkwo03FNnDhR3WSqTmqdzNn56wl46ZdgZGUBT7bzxYg22W0eiIiIyMgCp4MHD6p95DQ09UWjR4/GokWLMGzYMMTGxuLDDz9EVFQUmjRpgvXr18PPz8+AozYft++kYfySg4hPSUebWm54d0BjQw+JiIjIaFllSUUw6aTJOMXFxaFSpUpmdaUyM7NU0PTv6WhUc3XE2hc6qaJwIiIiS3K7BH/rjbrGyZAsocbpq3/OqqDJwdYa80e2YtBERERUDGacLDTjtOF4FJ776bC6/9WwQDzUoqahh0RERGQQzDhRkU5fvY3JK46q+0938mfQREREpCdO1VmYW0mpeGbJISSlZqBTXQ9M6dvQ0EMiIiIyGQycLEh6Rib+75cjCL+RBB83J3wzogVsbfhPgIiISF/8q2lBxeHT/z6NHedi4GRno4rBq1SwN/SQiIiITAqLwy2kOHzt0St48Zcj6v6cJ4LQr2k1Qw+JiIjIKLA4nPI4ey0eb/x+TN1/rlsdBk1ERET3iFN1FtAZ/Nmlh5Ccll0M/mqvBoYeEhERkcli4GTGpDP45N+OIiwmETUqO2HWiBawsbYy9LCIiIhMFgMnMy4On7vtPDafvAZ7G2tV1+TGYnAiIqL7wuJwMy0O33HuOkb/uB+ZWcBnDzfF8Da+hh4SERGRUWJxuIW7fDNJraCToGl4ax8GTURERKWEU3Vm5k5aBp5bdhg3k9LQrKYr3h/U2NBDIiIiMhsMnMzM+2tP4HhkHKo426m6Jkc7G0MPiYiIyGwwcDIjv+4Px68HIiAL574ZEYSaVZwNPSQiIiKzwsDJTByNuIV3/zih7k/u1QCd6nkYekhERERmh4GTGbQjuJGYiueWHUJqRiZ6BXjh+W51DD0kIiIis8R2BCbejiAjM0u1HdgZGoPaHhWw5oWOqORoZ+hhERERmQy2I7Ag/9t0RgVNzvY2+G5kSwZNREREZYhTdSZs44mrmLP1vLo/fWgz1PdyMfSQiIiIzBoDJxN1/nqC2odOjOvkj4GB1Q09JCIiIrPHwMkEJaakY8LSQ0hISUcbfzdM6dvQ0EMiIiKyCAycTExWVhamrDqOc9EJ8KrkgNmPB8HOhj9GIiKi8sC/uCbWjmDp3kv48+gV2FpbqaCpqouDoYdERERkMdiOwITaERwJv4nH5u1BWkYW3hkQoGqbiIiI6P6wHYEZupmYiok/HVZBU7+m3niqYy1DD4mIiMjicKrOBGRmZuHl5cG4EncH/h4VVOsBKysrQw+LiIjI4jBwMgHfbgnFtrPX4WhnjblPBsGFncGJiIgMgoGTkdtx7jq++uesuv/JkKZo6G18274QERFZCgZORiwqLhkv/RqMrCxgRBsfDG1Z09BDIiIismgMnIxUanqmKga/kZiKxtUr4b2BjQ09JCIiIovHwMlIfbbhNA6H34KLoy3mPtESjnY2hh4SERGRxWPgZITWH4/Cj7vC1P0ZjzWHr7uzoYdEREREnKozvs7hF64n4PXfj6n7E7rWwYMBXuX6+URERKQbO4cbUefw5NQMDJm9C2euxaOtvxt+erotbLkPHRERUZli53AT3bz3rTXHVdAk+89983gLBk1ERERGhjVORuLXAxFYdTgS1lbANyNawNPF0dBDIiIionwYOBmBkMg4vLf2hLr/Wu+GaFfb3dBDIiIiokIwcDKwuKQ0PPfTIdW3qWcjLzzbpbahh0REREQ6MHAy8Oa9k1cEI+JGMnzcnPC/RwNhLXN1REREZJQYOBnQvO0X8M+paNjbWqsml67OdoYcDhERERWDgZOB7Dkfiy82nlb3PxjUGE1quBpqKERERKQnW31PpNJVvbIjGnpXQsNqLhje2oeXl4iIyAQwcDIQP/cKWPV8B2RlAVZWrGsiIiIyBQycDIgb9xIREZkW1jgRERER6YmBExEREZGeGDjpMHv2bAQEBKB169b6XksiIiIyc1ZZsrsslcqOyURERGTef+uZcSIiIiLSEwMnIiIiIj0xcCIiIiJi4ERERERUuphxIiIiItITAyciIiIiPXHLlWJoujXIUkUiIiIyP5q/8fp0aGLgVIz4+Hj11cfHpzR+NkRERGTEf/Oln1NR2ACzGJmZmbhy5QpcXFxgZWVV4HnpLH7gwIF7jnAlIIuIiGBzTRNwPz9rU2aK37exjLm8x1HWn1fa71+a78ffxZajdRn8O5dMkwRN1atXh7V10VVMzDgVQy5gzZo1dT5vY2Nz30GPvJ5dyY1fafysTZEpft/GMubyHkdZf15pv39pvh9/F1sOmzL6d15cpkmDxeH3aeLEiff7FmQiLPVnbYrft7GMubzHUdafV9rvX5rvZyw/cyp7hv5Zc6rOgLgPHhGR4fF3MZUEM04G5ODggPfee099JSIi/i4m48eMExEREZGemHEiIiIi0hMDJyIiIiI9MXAiIiIi0hMDp1I0Z84c+Pv7w9HRES1btsSOHTt0nhsVFYXHH38cDRo0UL2iXn755ULPW7lyJQICAlQBuXxdvXp1aQ6ZiMisbN++HQMHDlSNDKVp8Zo1a4p9zbZt29TvbPndXbt2bXz33XcFzuHvYtJg4FRKli9froKft956C0eOHEHnzp3Rt29fhIeHF3p+SkoKqlatqs4PDAws9Jw9e/Zg2LBhGDlyJI4ePaq+PvbYY9i3b19pDZuIyKwkJiaq36nffvutXueHhYWhX79+6ne2/O5+88038eKLL6pASYO/iyk3rqorJW3btkVQUBDmzp2rPdaoUSMMGTIE06ZNK/K13bp1Q/PmzfH111/nOS5Bk/QX2bBhg/ZYnz59UKVKFfzyyy+lNXQiIrMkGSfJ0svvYV3eeOMNrF27FqdOndIemzBhgvqPVQmYBH8XU27MOJWC1NRUHDp0CL169cpzXB7v3r37nt9X/k+b/z179+59X+9JRETF/549ePAg0tLS+LuYCmDgVApiYmKQkZEBLy+vPMfl8dWrV+/5feW1pf2eRERU/O/Z9PR09budv4spPwZOpZwWzr/bcv5jxvCeRERU9O/Z/Mf5u5g0GDiVAg8PD7Vbc/5MUHR0dIH/kikJb2/vUn9PIiIq/vesra0t3N3d+buYCmDgVArs7e3VUtbNmzfnOS6PO3TocM/v2759+wLvuWnTpvt6TyIiKv73bKtWrWBnZ8ffxVSAbcFDdC8mTZqk2gXI/9nk/4jz589XrQhkdYaYOnUqIiMjsWTJEu1rgoOD1deEhARcv35dPZYgTPo1iZdeegldunTB9OnTMXjwYPzxxx/4559/sHPnTv6QiIgKIb9PQ0ND87QbkN+tbm5u8PX1LfC7WH5HS+sC+R0+fvx4VSy+YMGCPCuX+buY8siiUjN79uwsPz+/LHt7+6ygoKCsbdu2aZ8bPXp0VteuXfOcL5c//01en9uKFSuyGjRokGVnZ5fVsGHDrJUrV/InRkSkw5YtWwr93Sq/g3X9Lt66dWtWixYt1O/uWrVqZc2dO7fA+/J3MWmwjxMRERGRnljjRERERKQnBk5EREREemLgRERERKQnBk5EREREemLgRERERKQnBk5EREREemLgRERERKQnBk5EREREemLgREQGIbvNr1mzxmKuvmzJ9Omnn5b557z66qt48cUXy/xziCwVO4cTUakZM2YMFi9erO7L7vKyP1izZs0wYsQI9Zy19d3/VpMd6atUqQIHBweD/gTef/99FcBp9o4sC8eOHUO3bt1w6dIluLi4oCxFR0ejTp066jP9/f3L9LOILBEzTkRUqvr06YOoqChcvHgRGzZsQPfu3dUmqQMGDEB6err2PG9vb4MHTaUpLS1N53Oyieyjjz5a5kGT8PT0RK9evfDdd9+V+WcRWSIGTkRUqiQYkqCoRo0aCAoKwptvvok//vhDBVGLFi3SOVX3xhtvoH79+nB2dkbt2rXxzjvv5AlGJDPUvHlz/Pjjj2qX+4oVK+K5555DRkYGPv/8c/WZEjR88sknecYTFxeHZ555Rj1XqVIl9OjRA0ePHlXPyXg++OAD9VjGIzfNGIt6Xf7xyHjl+87euzuvzMxMrFixAoMGDSp2qrJy5craz5fAU8757bff0LlzZzg5OaF169Y4e/YsDhw4gFatWqlrIIHq9evX87yPfNYvv/xSwp8cEenDVq+ziIjugwQdgYGBWLVqFZ5++ulCz5FsjAQN1atXx/HjxzF+/Hh17PXXX9eec/78eRWA/f333+r+I488grCwMBVwbdu2Dbt378ZTTz2FBx54AO3atVOBTP/+/dWU4fr16+Hq6op58+ap5yUAGTZsGEJCQtT7/fPPP+oz5JziXifHRWhoqApsVq5cCRsbm0K/L5kyu3Xrlgp07sV7772Hr7/+WgWL8r3JtKcEcjNnzlRB5mOPPYZ3330Xc+fO1b6mTZs2iIiIUFODfn5+9/S5RFQ4Bk5EVC4aNmyogghd3n77be39WrVqYfLkyVi+fHmewEmyN5LhkYAqICBATQOeOXNGBTdSP9WgQQNMnz4dW7duVYHTli1bVBAmdT+aacEvv/xSZXp+//13lVGSrI3UY0nGSuO///4r9nUiNTUVS5cuRdWqVXV+X5I5kqBKMlf3Wuzdu3dvdV+mPCVw+vfff9GxY0d1bNy4cXkyeUKyfZrPZuBEVLoYOBFRuZAsjkw96SIBiWRWJIuTkJCg6qEks5KbBFS564S8vLxUUJK76FyOScAjDh06pN7L3d09z/skJyerjJUu+r5OgpKigibNayT4Kup7L4oU1+f+3kTTpk0L/X41ZFpPJCUl3dNnEpFuDJyIqFycOnVK5yqvvXv3Yvjw4areSLIrMjX266+/4n//+1+e8+zs7PI8lmCksGOSmRLytVq1aioDlZ/UE+mi7+sqVKiA4nh4eKgARrJT9vb2RZ4r9Vr55f7+NMFX/mOa71fjxo0b6mtxQR0RlRwDJyIqc5qpr1deeaXQ53ft2qWyN2+99Zb2mNTn3C8pTpe2BzIVJ9mqwkgwkz9g0ed1+pICcnHy5EntfQ35DA0p8JYsV2mQui0Jrho3blwq70dEd3FVHRGVqpSUFBUQREZG4vDhw6rp4+DBg1U7glGjRhX6mrp16yI8PFxlmWQqbNasWVi9evV9j6Vnz55o3749hgwZgo0bN6qaHykgl3qqgwcPqnMkMJICc+njFBMTo8avz+v0JVkfCcR27txZ4DmZmpRsm2TjZIWgkJqt2NjY+/q+d+zYoV2JR0Sli4ETEZUqWaEm01wSkMhSeSnQlkBIWhLoWnkmgZVko1544QWVlZEgRdoR3C+ZxpLC8S5duqgVabL6TqYEJRDS1AsNHTpUjVMKzSXIkWX8+ryuJKSY/KeffipwXD5X3ldW3Pn4+OD555/H7NmzVcbofsj3IKsSiaj0sXM4EVEZu3PnjlrxJxk1yWSpX75WViqrJlmt0rRu3Tq89tpragWjTDUSUelixomIqIw5OjpiyZIlaiqwrCUmJmLhwoUMmojKCP9zhIioHHTt2rVcrrM0xCSissPAiYjIAArbnoWIjB+n6oiIiIj0xMCJiIiISE8MnIiIiIj0xMCJiIiISE8MnIiIiIj0xMCJiIiISE8MnIiIiIj0xMCJiIiISE8MnIiIiIign/8HTk79PYpMDhsAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "### Estimate the signal response across a diameter range ###\n", "# Define diameters from 100 nm to 2000 nm for signal estimation\n", "diameters = np.arange(0.1, 2.0, 0.01) # in micrometers\n", "\n", "signal, noise = ops.estimate_signal_noise(ri=1.5, diameters=diameters)\n", "\n", "fig, ax = plt.subplots(figsize=(6, 4))\n", "ax.plot(diameters, signal)\n", "ax.set_yscale(\"log\")\n", "ax.set_xscale(\"log\")\n", "ax.xaxis.set_major_formatter(mticker.ScalarFormatter())\n", "ax.set_xlabel(\"Diameter (µm)\")\n", "ax.set_ylabel(\"Signal (Amps)\")\n", "ax.set_title(\"Estimated Signal\")\n", "fig.tight_layout()" ] }, { "cell_type": "code", "execution_count": 4, "id": "1f48ccbc", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk4AAAGGCAYAAACNCg6xAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAVidJREFUeJzt3QdYleX7B/AvG0EUAUFURBQXKk6cmavcpllq+ctRjiz7m+3MtLJh06zcLTPLzExbjly5B6K4cCIqKoiAiuz5v+4HzxEQ8DAOZ30/13XinPec856HF8Pb+7mf+7HKycnJARERERHdk/W9X0JEREREDJyIiIiISoAZJyIiIiIdMXAiIiIi0hEDJyIiIiIdMXAiIiIi0hEDJyIiIiIdMXAiIiIi0hEDJyIiIiIdMXAisjBLliyBlZVVkbf//vtP53NduXIFb7/9NkJDQ+96To7L+QwhLCxMff758+fL/dy6fl+yKcMvv/yCLl26wNPTE46OjqhduzZ69+6Nb775Jt9r5XxyXmNhbOMhMia2hh4AERnG999/j8aNG991PCAgoESB0zvvvIO6deuiZcuW+Z4bN24c+vTpA0MFTjKubt26qbEZwtSpU/HRRx9h/PjxeOWVV+Di4oILFy5gy5Yt+OOPP9T10dizZ48KqojI+DFwIrJQzZo1Q9u2bfV2fgkELDUYSElJwZw5czBq1CgsXrw433NjxoxBdnZ2vmMdOnSo4BESUWlxqo6IirRy5Uq0b98eVatWhZOTE+rVq4ennnpKPSdTekFBQer+k08+qZ3q00zxFDalJdmfAQMG4O+//0arVq1QqVIlNGnSRD3WTCPKY2dnZ7Rr1w4HDhzI9355/Nhjj6nzyHvl6+OPP64yORpyjqFDh6r73bt3145Ljmts2rQJPXv2RJUqVdT31blzZ2zevPmu7/+ff/5RmTQHBwf4+fnh008/1elPS1JSEtLS0uDt7V3o89bW1vecGtu5cyc6duyopvhq1aqF6dOnqyk+eW3eKUjNNV2/fj1at26trotkEr/77rt857t27RqeffZZlVGsXLmymj7s0aMHduzYodP3RES5GDgRWaisrCxkZmbmu8mxvNNHw4cPV8GS1OpIEDFjxgz1OiF/Sct0n3jzzTfV6+WWdwqqMIcPH1bTWK+99hp+//13FZQNGTIEb731lgoMPvjgA/z000+4efOmCggke6MhAUOjRo1UNmfDhg1qKiwqKkoFcLGxseo1/fv3V+cQ8+bN045Ljotly5ahV69eKmj64Ycf8Ouvv8LNzU3VHuUNnuT+oEGD1BSbfP+ffPKJeq3mey6Oh4cH/P39MX/+fMyePRsnT55UNU+6OnLkCB588EEkJyerMS5cuBAHDx7E+++/X+Q1femll/DCCy+oacDAwECMHTsW27dv174mPj5efZXrLD9L+T7kZyvTmSWpayOyeDlEZFG+//57+Ru80JuNjY32dZ9++qk6duPGjSLPFRwcrF4j5yzorbfeUs/l5evrm1OpUqWcS5cuaY+Fhoaq13l7e+ckJSVpj69Zs0Yd//PPP4v8/MzMzJzExMQcZ2fnnC+++EJ7fOXKleq9W7duzfd6Ob+bm1vOwIED8x3PysrKadGiRU67du20x9q3b59Ts2bNnJSUFO2xhIQE9f6C31dh9u/fn1OnTh3ttXVxcckZMGBAztKlS3Oys7PzvVael+ulMXToUPU9Xbt2Ld8YAwIC1GsjIiLyXVNHR8ecCxcuaI/JmGWcTz/9dLHXLiMjI6dnz545Dz/8cLHjIaI7mHEislBLly5FcHBwvtu+ffu0z2um4YYNG6YyLZcvXy6Xz5WpL5l60pCpOSGZD5k2K3g87zRcYmKiylRJNsfW1lbdZNpJpsZOnDhxz8/evXu3yryMHj06X6ZNao6kkF2ugZxLbnJfMmEyVaYh2aeBAwfq9H3K9Tt79qyaQnvjjTfUtJtksaTu6aGHHio2A7Vt2zY1jSaZq7zTe/KzKOqa1qlTR/tYxtywYcN8105I5koyhfK8XDs7Ozs1Jl2uHRHlYnE4kYWSwKS44vD7778fa9aswZdffqn+speanaZNm2LatGmqrqi0ZFosL3t7+2KPp6amao+NGDFC/UUv9T4SmMh0m9T89OvXL9+UXlGuXr2qvj766KNFvkYCKzmnBFM1atS46/nCjhVFAhOZApSbiIuLU58tNV3r1q1T4y6MvM7Ly+uu44UdE+7u7ncdk7qsvNdEpgxlOm/ixIl49913VVBmY2OjriUDJyLdMXAioiJJjY/cJGjau3cvZs2apYIXKUiWDEpFkponCTikRuf111/XHpexaep37kWTwfnqq6+KXMkmwUlGRoYKnqKjo+96vrBjupIAZ8qUKaqm6NixY0UGTvI6TZBXXp8ttV2S1VuwYEG+47du3Sr1OYksEQMnIronyV507doVrq6uqij70KFDKnCS40KXbE9ZSSAj01uaz9SQgvK8Re2a8RY2Llk9J9+D9Hl67rnnivwsyXbJqj4pXpeicM10nQQZf/311z3HKoFXQkJCoZkgTXanZs2aRb5frvXatWtVwbsm2JMMmKxyLMv1K3jtpAhdCud9fHxKfV4iS8PAichCScZDs0Iur/r166N69epqBd2lS5fUsn3px3Tjxg188cUXavpJ/mLXvFaWv8sqOJn6k3ojCQiKCwpKS6blZPpQAhkJJiTrJbVA3377rQqGCvaoEtJDSeqSJPCRdgISyEi2SWqcJEsl02ayLF+W6svKNPmqycjIdJbUPcnqNpnikuBMVvFJq4R7ZbgkOybjk7YIDzzwgApMpD5LMk1yDeVaSf1UUWQ6VAI0ufZyX66x1CdJ7VVh7Qx0ISsU5XuSjJ38/E6dOoWZM2eq61LYnwMiKkKeQnEisvBVdXL7+uuv1ev+/vvvnL59++bUqlUrx97ePsfT0zOnX79+OTt27Mh3vuXLl+c0btw4x87OLt9qrKJW1fXv3/+uMcnrJk2alO+YrByT45988on2mKzGe+SRR3KqVaumVqn16dMn59ixY+q8o0ePzvf+OXPm5Pj5+amVggVX/m3btk2NQ1aeybjle5THshovL1nRFxgYqL5/WSH34YcfFvp9FZSWlqZWJcr1k/c5ODiolW9NmjTJefXVV3Pi4uLuuYpNrrOs7JP31qhRI+eVV17J+eijj+5a6VjUNe3atau65R3Tyy+/rL5XGUvr1q3VykW5bnKOe42HiHJZyX+KCqqIiMh4SP8p6WV1+vRpQw+FyGJxqo6IyAi9+OKLqru6TPPJ1KBMh27cuFFNTRKR4TBwIiIyQlJTJXVmspJOCrtlq5Qff/wRTzzxhKGHRmTROFVHREREpCN2DiciIiLSEQMnIiIiIh0xcCIiIiLSkcUXh0s33itXrqgmeVKASURERJYlJydH7QwgzXvv1WDW4gMnCZq43QARERFFRkaqnRIYOBVDMk2aiyVbOhAREZFlSUhIUEkUTUxQHIvPOGmm5yRoYuBERERkuax0KNlhcTgRERGRjhg4EREREemIgRMRERGRjhg4EREREemIgRMRERGRjhg4EREREemIgRMRERGRjhg4EREREemIgRMRERGRjhg4EREREemIgRMRERHp3a3UDEz66SDWH4s26attsYHTvHnzEBAQgKCgIEMPhYiIyOzN/y8c/xyNwsRlITBlFhs4TZo0CWFhYQgODjb0UIiIiMze1ZupMAcWGzgRERFRBbIyj6vNwImIiIj0H3BYmUfkxMCJiIiIirX+WBSGLtyNS9eTS32lzCNsYuBERERE9zBx2UEEn7+Ot/88Xi4Zp6zsHJ3ek5mVjdSMLKP6+TDjRERERDq5lphe+oAjT8SRomMw9PyKUAS9vwmXb6TAWDBwIiIiIr1Pt2Vn37mfnJZ5z9enZWbhnyNRuJWaieX7LsJYMHAiIiIivUvNvJNlSkq/d8bpZNQt7X1mnIiIiMiipOQJlpLT751xOnLphvb+gQvxMBbMOBEREZHepeSpa0q+HUR9sekMvtx8ptDXSzG6RmR8Ck5fvZOBMiQGTkRERFTsyrbykJoncEpKy8TVhFR8vuk0Zm88jcj4/G0O5PlNJ66q+x6V7dXXVQcvGcVPiYETERGRGQq5cB0xCWXf5kSXeiRdaLJMmmm78JhE7eP9EfHaMf8RelntaSevr+vuhPcGN1PPrTxwSQVUhmZr6AEQERFR+Tp08ToeWbAbTvY2CJvZp0znyluPlFFM9iknJwerDl5GW99qqOvhXOxUnQRjsYlp2sf7IuLQqo6rGnNew4Pq4IEmXiqAOh+XjKV7LuCZbvVhSMw4ERERmZmtJ2PuyvKUVt4sT3Ix5/slOBIvrzyMRxfuKfT51ALF4eHXkrSPd52Nw0frT+Z7vauTHUZ29IWtjTX+r0eDuwrGDYUZJyIiIjOja4PJ7aev4btdEXjnoabwdb87SyQS0+6cK7GYqbLfQnJrkPJmkooakwRg4dcS87UbKNhy4PmeDVDZITdMGdSyJnzcnNDOzw2GxowTERGRmUnNyM7XSLIob6w+iv9OXUPfL3YU+Zq8zSqTigmc8hZ4y7RdcYFTfFI6jl9JUPfruDlpjzvYWqNrw+ro07QGRnbw1R6XrJMxBE1mFTglJyfD19cXL7/8sqGHQkREZFBJeeqSpPN2UTTPSQbowPnCeyXlzTIlp2chu5B95tIzsxFz606m6dbt90gReNTNFPWevMHc4u3nVPBUy7USPn40UHu8fvXK+OGpdlg4so0KloyRcY6qFN5//320b9/e0MMgIiLSi/XHotDhg83Yey7unq+NzbOnXHGBk4vjnYqdfbdXthUXhInkQqYBC/ZYunYrTWWdnloSjC4fbdVmlwp6pXcjdKjnjhceaKgej+50J8tkrMwicDpz5gxOnjyJfv36GXooREREerHuWDSiE1Kx+XZ/o+LkbUOQkJJR5Ovi8gRYkgG6V41TUfvMXbqeclfgtOdcnLplZufg7yNX7npPrwAvVbsknn+gAUJnPIhhbX1g7AweOG3fvh0DBw5EzZo1YWVlhTVr1tz1mvnz58PPzw+Ojo5o06YNduzIPxcr03OzZs2qwFETERFVrOibucFQdELhxdd5SXNJjYTUwgMnWdmWt+4oroii7oKBUmIhgVP0zbsDpx92n9c+Do28oa1hWvVMJ7w7uBk+H95S/b2v4epkn++xsTJ44JSUlIQWLVpg7ty5hT6/YsUKTJkyBdOmTcOhQ4fQpUsX9O3bFxcv5u6U/Mcff6Bhw4bqRkREZK40wdDV2wFUUaQY/HrynWApISXzntkm9biIjFPBgvCk2xmoHWeu4dfgyEKDuYvxydh2+tpd04CV7G3QxreaKvx2vr1iztQYfNQSBMmtKLNnz8bYsWMxbtw49XjOnDnYsGEDFixYoLJMe/fuxS+//IKVK1ciMTERGRkZqFKlCmbMmFGB3wUREZH+SL3Q1dvBiUzXFUeyPXndKiLjVLBtQMFAqqipuqT0TGRl52Dkt/vV4wZele/KOP1+8JIqBq9kZ5Mvq+Xp4gBTZ/CMU3HS09MREhKCXr165Tsuj3fvzu0uKsFTZGQkzp8/j08//RTjx48vNmhKS0tDQkJCvhsREZExS0i9M60mgVNhy/018q5uy31v4YGTpqZJMztWVI1T3s7hmgxU3h5MUvgddTsL1riGi/qqaW45PMgHgbWral87ulNdmDqjDpxiY2ORlZUFLy+vfMflcXR0dKnOKYFW1apVtTcfH+MvRCMiIvMj01zP/XywyMCmqJolWfp/I89U3F2vLTCVd6+puoaeucFOXFLuSriCbhYoLr+RnIEjl25qHx+7fFM7vl5Na2iP16vujIld62PhE23QwLMyWtSuiqFtTP/vXINP1emiYLGY/GALKyAbM2bMPc81depUvPjii9rHknFi8ERERBXt1VVH1Fd7W2vMHtZSp8Jw7eOEVFRzti/0tQVXuBUVmMUm5Wam/L0q49TVW8jIylH9l6o42uV7nTwn/DycERGbhEOR12FrbZ2v8FuTcXq0dW081KImLsYnob2fu7aO6d8X7jeJwm+Tzzh5eHjAxsbmruxSTEzMXVkoXTk4OKgaqLw3IiKiipR3ZdrvBy8XuaJNo2BdU3F1TlKYLexvN5Asqo+TJuNU27USnO1ttMckOPpg7QlcuZGiximPxYT766mvu8/G4eDF69rznIy+hbTM3OaWnlUc4O9ZGT0ae+Ur/jaXoMnoAyd7e3vVfmDjxo35jsvjTp06lenc8+bNQ0BAAIKCgso4SiIiIiAjK1sbZNzLuTw1QmLl7X3edJ1+K25lnSZwalqrSrF9nDQ1TW7O9nCrnJu9ik9Kw+yNp1Vn726f/oedZ2Ihs3feVR3Rr7k3rK2Ac7FJaqrO1tpKvVdDMlKOdrkBmDkzeOAkK+FCQ0PVTURERKj7mnYDMq32zTff4LvvvsOJEyfwwgsvqOcmTpxYps+dNGkSwsLCEBwcXC7fBxERWXbQNOLrvej+6X8IuXAnG1OUvMXVBfd5K2vGKfL67cCpZpVip+okABJeVRzh7uyg7TguHco1tVQTl4XcPldVVK1kh5Y+rtr3j+tSD3NHtNI+fvuhprAEBq9xOnDgALp37659rKk/Gj16NJYsWYLhw4cjLi4OM2fORFRUFJo1a4a1a9eqfemIiIiMwecbTyP4fG7AtD8iXvUqKs6526vOimoNUJCm+FqW88uqubzF4nnJnnCX4lO0wU5RU3XXk9Jx5FJuU8r29dzw1+Hczt6yX53UOhXU7Hb26tOhLfDzvouqG/jzPRuovkyfD2+hap5kc15LYPDAqVu3bsUuqxTPPvusuhERERkbKdxesC1c+1g2tdU14xRUt5oKuPLuLVeYyNvBUOs61bD+eLS2p1NBV2+lIj0rW02jNfRyuWuqTuqWpHhcAi/5q7eRlwu8q1aCj5uTev7rHRHqa3s/N7WaTuqX7GystFuh1KteGW8OCMj3mQ+3qg1LYvDAyVCkxklu0u6AiIiotILPx6sgRONCXPHTbiI8JjfjJCvPcgOnojNOqRlZ2kDrvgYeKnAq2ORS4+Ltz67pWgket+uWpCO4NKzMzM7GY4v3amugxP0NPdTXMZ3q4oc957XfR2d/D/VZ3+86j8k9/NX5yEhqnAyFNU5ERFQeNDVNkr3RpV5JghhNEXk7P7diu3aLszGJamrM1ckOLWrn1hjF3Cp8qu7C7c+u4+aEWq6V1N5wsuLtQlwSVgRH5guahBR8i7oezpjQJXfVXM/GnnjqPj+V3frq8VZocPv7IgsPnIiIiArb523ijyF49+8wnS+OZmn+4Fa11FeZCpPgqChnYm6p6TTZjkQTCMmyf8ksFSYsKneHiyY1qqjl/kKm9qSeqaAdZ2K1Hbxtbay1nbzlHAv/uzOdKO0Hvny8FVrVuVOL9Xrfxtj1eg98M7otKpvoPnIVgYETERHRbatCLqupsG93RhSZ1Sm4HYlsOSIGBHqr2iIJiooq3hayxF+TbapSyVY1wBRFTb+d0ARO3lXg7myvtkiRwCw+Of2urVA2huX2PRzQoqb2PZpeUVdupqpg7eS7fXBoRi/VqDIv6bUkWSpz6rmkDwyciIiIAGRmZWPR9jtZmT3hcfe8LtLPSIKYGlUcUbtaJXUTBafE8tp+O3Dq0sBDBSket3shFVXndCdwys0iSfAkYgoUiMtUnGysW9fdSW1vkjdw2nIyRn3t7O+uei1pgjUqOYu9cmyASUREea09Fp2vsHvX2dwApzj7zsWrr9J+QIIgzeq0ogInmY7bdy43ILv/9vJ9D5c702+FTR1qMlqaIMijcu7rr90OtBZtC0fQ+5sw8/b04tC2PtqsUcDtXk4a3Rp58odeRhYbOLE4nIiI8lp5IFLbIkDsOhtXbLsceW71odyO390aVdcWZRdXIL73XJwq1pYMlWx8mzcQKmzblb8OR6k+TPJ6Tb2SZxVH9TUmIVUFYl9sPqOd5nu6az3t1iiaYMvRLvevehtrK3RvzMCprFj9RUREFk9qlTTZoxkDmmLIgl24fCNFZaBkxVlhpI3A+bhkONnbaFenaQKnwloSSDG3BDniwQAvbVbIvYipOnn9dztz+yqN6uSrpulE9TwZp/9OXUNyepbqtbRnak9tEKYhRd6/Pt0Rhy7eUIGa1DBR2TBwIiIiiycb10pRt49bJdUlW1abSQfwXeGxRQZOmgyVFIVrNrSVBpHi9NVbd71+1cFLKoCRFW3/18Nfe7ywqToJ5F5YEapWw0nG6PGgOtrnNCvrpMbpRNQtbR+mgkGTRmBtV3Wj8mGxU3VERGSeZEuQLzadKXS5flG2nMotnu7eyFNlgjrUc1ePQ25vo1LQ8Ss3sSb0sramSCPwdlG2BE4p6XfaC2w5eRXTVh9T9yf18NdOtxVWsyStCcZ8F4wNx6/C3sYaHz0SiGp5NtOVbVc004GbT1xV9/sH5l8hR/pjsRkndg4nIjI/shLujdVH1f3G3i7o3bTGPd8jtUr/nbwTOAnNXnMht3s05SV1RZINkj3degV4oW2efelkw1zNfnJhUTfRxtdNrYp7ZtlBldHq39xb22hSQ9Ph+9rtVXKfrD+J/efj4eJgi++fDELburlNMjWq3w6cNt8es0y/aVbRkf5ZbMaJxeFEROZFAhpN0CTmbT17z71QRWjkDdXjSLpsazJNLX1cVb8kqVXK219JWhZMXn4Ip68mqkzRrCHN7+p71LxWVW2rAumtNOmng6ogXDbBnfNYS22tkoZmT7mjl2+q/eGk55KQBpUFgybhV2DqsH+gN3svVSCLDZyIiMi8LN5+Tm1lIhkZafQogYumZ1JxZD82IdmgSvY26n7VSnZo6OmSrzO4eOevMPwbdlX1QZo7ohXcC6kran47+yOB0IL/wnEuNkmtivt8eEvYFQiaNFu1SIF4SkYWPlx3ArfSMlWtlQRahWlas6oKljRk3FRxGDgREZHJk0zQsr0X1P1p/Zrg8Xa5xdTf7DhX7PuibqZg7dEodV/2Z8ur9e0puIO396L7LeQSftx7QWWivnyslTY7VZCmzkk6hH99+/PfGdQUbnnqlPKytrZCJ//czXaX788tOB/axkcdL8p7g5qhoVdl1dBS83lUMRg4ERGRyZPO2FJXJJkbaQ0wsqOvtubpVmpGke/7dkeE2kC3vZ8bmt2eYtPQ1Dnti4jH4cgbeHNN7jTg8z0boE+zomunmtfSbMSbpqboZGsVqYUqzn3+7vlaCAwPulNwXhgpFt8w5X78NK4Dp+kqGAMnIiIyecv3X1RfH21TW02jSR2Q3CQokkaWhTlwPh7f7YrQNo4sSIIpyS5JDdTQhXvUdibS6HJyjwbFjkWmCp++vx78PSurabiZg5reM7jpfDvjJD4Y0lwVmd8L95QzDItdVUdERMZDirglO7TnXBxGdayrXTmmiys3UrDt9DV1P2+mRmqEpObpv1Mxd2WIbiSnY8qKUEjHgiGtaqFH47szQrJ9ysyHmuLtv8LUijjpwj13ROtip9A0pvZrom66ql3NSRWai4Kb75JxsdjAie0IiIiMw4W4JEz+JVRNh4kD56/j5/Htdc6oyF5tEgB1rOeubUApZHuRJbvPq+7aEphpzicr3cZ8H4xL11NUEbbUHxVlZMe6atXbxrCraisTmUbTF01dFhk3i52qYzsCIiLDk/qjp5YEq6BJ2gHINJtknX4Jzi2SvpdL15Px8+1pusk9G9w11SZdt6MTUnHqdidv2TR3wo8H1PSbq5Mdvh0dBBdHu2I/o309d7w5ICBf00qyXBYbOBERkWFJZ++Xfj2M8Gu5y/X/e6UbXu3dSD33wT8ncDUh9Z7n+HLzGdWIUlaXdayff5Wbo50NOtfPrR36ZP0pFTQ9vzxU1TzJtidLnmyn7aFEpCsGTkREZBBzt57N7YlkY42FI9vAu2olPNnZDy18XFUvI03hdlF2h8eqFgHipV65AVdBLzzYUGWxpMt2+w82Y/3xaPX461FtVZNLopJi4ERERBVO9lj7fNNpdf+9wc20QYyNtRUmdauv7v8aHKm6gRdGunk//0tucfewtrXRus6dbU/ykhYDcn5xIzlDdfte9EQbbd8kopKy2OJwIiIyjGOXb2LKL6GQ3VBGdvDFsAI9i3o28VL7r12+kYJ/jkThkTa18z0vm+c++1OICp6kCeQ7D+UGRkUZ1tZHFYffSs1UBdjOeizwJvPHjBMREVWYQxev4/Gv96qpuHZ13TB9QMBdr5Gs04j2uSvMpFN3Xpri7uDz19UmuPNGtNZuk1Kc4UF1MK5LPQZNVGbWltyOICAgAEFBQYYeChGRyYlNTMP6Y9GYu+UMQi7E6/Sef49H44lv9qnMT1vfavh2TFtVb1QY6cdkZ2OlVr/JlJ2IT0rHyG/3Y8eZWDhJcfdTQWjA4m6qYFY5umwdbcYSEhJQtWpV3Lx5E1WqVDH0cIiIjJoEL+//cwKrDuYWZWtID6Wp/RojsLZrofvIyeq3L7ecVY9lBdzikW3vmf2Z/e8p9R7JQD0W5KP6Mcn0nWSaFo1qg063V8wRVWQswMCJgRMR0T3Jv7HXhF7Gu3+fUMGTkO1EpLv2ttMxqiWArbUVXu/bGGPv89M2mzx+5SZeW3UExy4nqMdjOtXFtP5NYGdjrdNnvrzySL4grXa1SvhuTBDbCFC5YuCkp4tFRGSpnb3fXHNMTZFpAqZZjzTXrmSTLNB7f4dh3bFo9bhedWcE+bqpppMy1SaqVrJTe7YNalmrRJ+dkZWND9edxPWkdLUSrldTL1S5R8NKopJi4KSni0VEZGldvaXX0vc7z6u92qQe6fmeDdTWIwUzRpIdWrbvIt79Owzpmdna47KtW7/m3pgxMACeLuy8TaYfC3BNJhER5SN7uclqtsXbz2mn5bo08MDMQc3g5+Fc6NWSqTlpLTCguTeCz8erTJNXFUcVNJVkw14iY8fAiYiIlJhbqVi254IKmq4nZ6hj9TycVU1Sj8aeOm26W83ZHr2a1lA3InPEwImIyILJircdZ3O3LpF2AVLkLeq6O+H/ejTAoJY1YatDITeRpWDgRERkicHSmVjVh2nTiauIuz0dJ1rXccXY++qhd1MvBkxEhWDgRERkIWT129I95/H7wctquxKNak52arXbo21qq73diKhotpbcOVxuWVmFbyBJRGQuZEn/dzsjMGfTGaTc3jTX3dleFW73bloD7eu56dRXiYjYAJPtCIjIrF2MS8bEZSEIi8ptQClbnYy/vx66N/IscrsTIkuj93YEkZGROH/+PJKTk1G9enU0bdoUDg5cbkpEZEz2hMfhmZ9CcCM5Q03HvdGviZqO02V1HBGVMXC6cOECFi5ciOXLl6vAKe8Wd/b29ujSpQsmTJiARx55BNbW/FcMEZEh/bzvImb8cQyZ2TloUbsqFo9qq/oqEVHZ6BThPP/882jevDnOnDmDmTNn4vjx4yqdlZ6ejujoaKxduxb33Xcfpk+fjsDAQAQHB5dxWEREVNoVc2//eRxvrD6qgqaHWtTEiqc7MmgiqsiMk2SUwsPD1bRcQZ6enujRo4e6vfXWWyqIkuxUUFBQeY2RiIh0EJOQihd/PYydZ3P3lHuldyM8260+p+aIypFVTt45NwvEveqIyBz8cyQK09YcVfVMTvY2mD2sJfo0Y/duIoMXh6ekpKj6JicnJ/VYskurV69GkyZN0Lt375KejoiIyuBU9C18suEkNp2IUY+b1aqCz4e1RAMvF15XIj0oceA0aNAgDBkyBBMnTsSNGzfQvn172NnZITY2FrNnz8Yzzzyjj3ESEVEeRy7dwPe7zmNN6GXIvIGNtRUmdauP53o0YJsBImMKnA4ePIjPP/9c3f/tt9/g5eWFQ4cOYdWqVZgxYwYDJyIiPUlKy8Sfh6+oFXNHL9/UHu/f3Bsv9mqI+tUr89oTGVvgJL2bXFxyU8D//vuvyj5J+4EOHTqoaTsiIipfx6/cVMHSH6FXkJiWqY7Z21ijb/MaGHufHwJru/KSExlr4OTv7481a9bg4YcfxoYNG/DCCy+o4zExMfcsqCIiIt2kZWbhz9Ar+GnfRYRG3tAe9/Nwxoh2dfBIm9pwc7bn5SQy9sBJpuNGjBihAqaePXuiY8eO2uxTq1at9DFGIiKLmo5bvv8ivt5xDlcTcjfitbOxQq+mNfC/dnXQsb472wsQmVo7Aml6GRUVhRYtWmi7hO/fv19lnBo3bgxTwnYERGQM5FexFHq//88JxCamq2NeVRwwulNdDG3jg+ou3NaKyGT3qqtRo4a65dWuXTuYknnz5qlbVlbuTuFERIYSfi0Rb64+hj3n4tRjX3cnPNO1Ph5uXQsOtjb8wRCZcsYpKSkJH374ITZv3qzqmrKzs/M9f+7cOZgSZpyIyJDWHLqMqb8fRUpGFhxsrTG5ZwOM71KPLQWIzCXjNG7cOGzbtg0jR46Et7c359qJiEohNSML7/4dpoq/RWd/d3w4JBA+brnNhYnIOJU4cFq3bh3++ecfdO7cWT8jIiKygD3lxv8YgsORN2BlBUzu0UBlmqSJJRGZWeBUrVo1uLm56Wc0RERm7uilmxi/9ACiE1Lh6mSHOcNbolsjT0MPi4h0lLskrgTeffdd1ZJAGmESEZHu/gi9jKGLdqugyd+zMv6Y1JlBE5G5Z5w+++wzhIeHq61W6tatq/apK7glCxER5a9neu+fMCzbm1vP1L1RdXzxeCtUccz/+5OIzDBwGjx4sH5GQkRkhiJik/B/yw/i2OUE9fi57v544cGGrGcisqQGmOaE7QiISF9ZpoXbwjH/v3CkZ2ajmpMdPmc9E5FlNsAkIqLCZWfn4N+waHy47iTOx+XWgnZp4IGPHglETddKvGxEJk6nwElW0Z0+fRoeHh5qVZ2VrJ8tQnx8fHmOj4jIJEhWSbZMkSzTuWtJ2i1TZgxoin7Na7DnHZElBU6ff/45XFxc1P05c+boe0xERCZVw/RL8EWsCrmk3WOuiqMtRnWsi6e71oMLC8CJzAprnEowr0lEJNIys7D+WDR+2R+p3V9OeLo4YFwXPzzerg4DJiITovcaJ9kYd/Xq1Thx4oRKPzdp0gSDBg2CrS1LpojIfF2+kYIluyLwW8glXE/OUMekcqFbw+p4rF0d9GjsCTubErfHIyITUuJI59ixYypIio6ORqNGjdQxqX+qXr06/vzzTzRv3lwf4yQiMpgzV29h4bZzqoFlZnbuQmTvqo4Y1tYHw4J8UItF30QWo1Sb/DZt2hQHDhxQheLi+vXrGDNmDCZMmIA9e/boY5xERBXuyo0UfLz+JNaEXtEe61TfHWPv81Mdv7m3HJHlKXHgdPjw4XxBk5D777//PoKCgsp7fEREFS4pLROLtoVj8Y5zSM3IVsf6NK2Bid3qo6WPK38iRBasxIGTTM9dvXpVZZ3yiomJgb+/f3mOjYiownswrTp4CZ9sOIWYW2nqWLu6bpg+IADNa1flT4OIdAucpNpc44MPPsDkyZPx9ttvo0OHDurY3r17MXPmTHz00UcVfklv3bqFHj16ICMjQxWty9jGjx9f4eMgItO271wc3v0nTLs1Sh03J0zt2xh9mrEHExGVsB2BtbV1vuZtmrdojuV9LMFLRZLPS0tLg5OTE5KTk9GsWTMEBwfD3d1dp/ezHQGRZZOA6astZ7HzbKx67OJgi+d6+GNM57pwsLUx9PCIyBTbEWzduhXGysbGRgVNIjU1VQVSFr79HhHdQ0p6FradjsF3O89j//nc3Q5sra0wPMhHbcDrUdmB15CISh84de3aFfqyfft2fPLJJwgJCUFUVJTqDzV48OB8r5k/f756jTwvtVXSvbxLly7a52/cuKHGeObMGfU62RqGiEgjKzsHF+KScPjSDfx7/Cr+O3UNKRm52XF7G2sMbVsbE7vWh49b7j/CiIiKYvCOlUlJSWjRogWefPJJPPLII3c9v2LFCkyZMkUFT507d8aiRYvQt29fhIWFoU6dOuo1rq6uarWfFK0PGTIEjz76KLy8vAzw3RBRRZLscnJ6FuIS0xGblKa+xiWmIS4pHbGJaWoLlIjYRJy5moi0zNzVcRrSe2lAoLeakvOuys13icgEt1yRGqmCGaf27dujdevWWLBggfaYdCqX18yaNeuuczzzzDOqWHzo0KGFfobUQ8kt77ymj48Pt1whMmI3UzIQdiUBx6/cxNmYRFy5mYqoGymIupmKxLRMnc7haGeNBp4uuL+hB/o09UazWlW48S4RVcyWKxUlPT1dTeG9/vrr+Y736tULu3fvVvcly1SpUiX1jco3LlN/EjwVRYKtd955R+9jJ6Ky1SBJsfbGsGjsPRePi/HJ9wyKpC7JvbIDPJzt4V5Zbg5wd7ZH7WpOaFzDRU3DsWElEZVViQInSU5dvHgRnp6eKljRt9jYWFXsXXDaTR7Lli/i0qVLGDt2rBqb3J577jkEBgYWec6pU6fixRdfvCvjRESGlZGVjbVHo/D3kSjsOHNN23gy79Ra05pV0Ni7Cmq7VoK3q6OaYpOtT5wdjPrfgERkyYFTgwYNcPz4cfW1ouRthaAZh+ZYmzZtEBoaqvO5HBwc1I2IjMON5HT8tO8ilu45j6sJafkCpQcDvNC9sScCa1VFNWd7g46TiKjEgZP0c5KAKS4urkICJ1kdJ+0GNNmlvF3KWfxNZNqib6Zi3taz+C3kknaFW3UXBzwe5IPezWogwJs1SERkfEqc3/7444/xyiuvqGJtaTapT/b29iqjtHHjRjz88MPa4/J40KBBZTr3vHnz1K2iG3YSWbpbqRlYtO0cvtl5Zx84CZJk49wBLbzZdJKIzGtVnWzoKx26MzMzVWBTsNYpPj63mZyuEhMTcfbsWXW/VatWmD17Nrp37w43NzfVbkDaEYwcORILFy5Ex44dsXjxYnz99ddqutDX1xdlxc7hRBUjPTMby/dfxBebzyA+KV0da+tbDS/2aoiO9dy5wo2IzHNVnTSfLE8HDhxQgZKGpnB79OjRWLJkCYYPH66mBmUvPGmAKVmutWvXlkvQRET6J/82W3csGh+vP4nzcbmr4+pVd8brfRqrGqaCNYxERMbMqPo4GQIzTkT6sz8iHh+sPYHQyBvqsbQMeOHBBhje1ge2Nta89ERkGX2cwsPD8f3336uvX3zxhWpPsH79erWsX7ZEMQWscSLSn2OXb2LOpjPYdOKqeuxkb4PxXephwv312DqAiCwr47Rt2za15YlsfyLNJk+cOIF69eqpovH9+/fjt99+gylhxomo/Bw4H4+vtpzFttPX1GNpOPlYkA+ef6ABPF0ceamJyPIyTtLF+7333lO1SC4uLtrjUqck2ScisiyZWdnYfuaaWim3LyJ3cYi1FTCoZS1M6u4Pf8/Khh4iEVG5KXHgdPToUfz88893Ha9evboq4iYiyyB7x/1+8BLWhF5RG+oKOxsrPNqmNiZ2rQ9fd2dDD5GIyPCBk6urq1rd5ufnl+/4oUOHUKtWrfIcGxEZWTuBw5duYPfZOKw/Ho0TUQna59yc7fFwq1qqF1NNV/1vx0REZDKB04gRI/Daa69h5cqVahlxdnY2du3ahZdffhmjRo2CqWBxOFHxsrNzEBaVgN3hsdh1Ng7B5+ORnH6nYay9jTV6NvHEkNa10a1RddhxlRwRWYASF4dnZGRgzJgx+OWXX1R/FltbW9V9WwIq6bskW6SYEhaHE+WS/5/PxSZh99lY7A6Pw55zcbiRnJHv8khmqWN9d9zn74G+zWrA1Yn7xxGR6StJLFDqPk7SikCm5yTjJB2/K3LT3/LEwIks2ZUbKSpI0gRL0Qmp+Z53trdB+3ru6FRfbh5oXMMF1lL5TURkRvS6qk5aEDRu3Bj169dXt7yZqD179uD+++8v3aiJSO/SMrOwJzwOm0/EYOfZWETEJuV73t7WGm3qVENnf3d0rO+BwNpVOQVHRFSWwKlbt27w8vLC77//rvaOy7tHnbQk4Ka5RMZF9oXbejJGNaPcfvoakvLUKUnyKLC2q8oodfb3QBvfanC0M63pdiKiilSqzuGPPfYYevbsifnz56t6Jw1T2r2FxeFk7ivgJFBaERyJHWeuITvP/5peVRzwQBMvdGvkifb13FDF0c6QQyUiMiklrnGS4m9pR7Bz5061Ee/48ePx2WefISYmBjVr1jS5jBNrnMicnIq+pYKlNaGXVaZJQ2qTegV44YEALzSrWZV1SkREFVXjpImzhgwZono5DRo0CGFhYewaTmTAtgEbjkdj0fZz2s10haeLg2pGObStD/w82IySiMhgU3UasppO9qcbPHiwmrojooqTlZ2Df45GYe6WMzh9NVEds7W2UtNww4Jq4/4G1WHL3kpERIYNnGR6rlKlO52Ba9SooTb+nTBhglpxR0T63xvuj9ArmPffWZy7lrsqzsXBFmM618XoTnXhUdmBPwIiIj0pdR8nc8EaJzIlm8Ku4t1/wnAhLlk9rlrJTm1zIgGT3CciIiOocTpy5AiaNWsGa2trdb84gYGBMAVcVUem5GJcMt756zg2n4xRj92d7TGuSz2M7OiLyg5lmnEnIqLyzjhJwBQdHQ1PT091X/aoy/s2zWP5ylV1ROUnNSMLC7eFY/5/4arFgJ2NlQqYnuvuD2cGTERExplxioiIQPXq1bX3iUj/tpy8irf/DMPF+NxpOdkf7u2HmsLfszIvPxGRgegUOPn6+mrvX7hwAZ06dVKb++aVmZmJ3bt353stEZVcQmoG3v7jOH4/dFk9rlHFEdMHBKBf8xoqq0tERIZT4uII2VZFGmDKtF1ekt7ilitEZRN8Ph5TfgnF5RspajsUmZZ7vmcDTssRERmJUjXALOxfvXFxcXB2ZpM9otLIyMrGF5vOYP5/Z9X2KD5ulTBneEu08XXjBSUiMsXASTqFCwmaZH86B4c7vWKkIFxW28kUHhGVzLlriXhhRSgOX7qpHku377cGBsCFe8gREZlu4CTV5pqMk4uLS74mmPb29ujQoYPat46IdPdrcCTe+vM4UjKyVB+mWUOao19zb15CIiJTD5y+//579bVu3bp4+eWXTX5ajn2cyNBtBmb8cQy/HrikHnf2d8enQ1vAu+qdf5AQEZHxYefwEvRuICoPkfHJmLgsBMevJKgC8Jd6NcIzXevDWh4QEZHp93Fq3bo1Nm/ejGrVqqmNfYtbEn3w4MGSj5jIQmw9GYMpK0JxMyUDbs72+OrxVujs72HoYRERkY50CpwGDRqkLQYfPHiwrucmotuysnPwxeYz+GrLGUjT/ZY+rpj/v9ao6cqpOSIiU8KpOk7VkZ7dSE7H87+EYtvpa+rxyA6+eHNAEzjY2vDaExGZ41QdEZXOmau3MG7pAVyIS4ajnTU+eLg5hrSuzctJRGSiShw4SZ1TYTVOcszR0RH+/v6qz9OTTz5ZXmMkMkmbwq6qeqbEtEzUrlYJi0e2RUBNLkAgIrKowGnGjBl4//330bdvX7Rr1071dQoODsb69esxadIktQnwM888o/auY18nskTy/8SCbeH4ZMMpVc/U3s8NC55oo4rBiYjIwgKnnTt34r333sPEiRPzHV+0aBH+/fdfrFq1CoGBgfjyyy8ZOJHFSUnPwmurjuDPw1e09UwzBgbAzsba0EMjIiJDFIdXrlwZoaGhakour7Nnz6Jly5ZITExEeHi4Cp6SkpJg7FgcTuUl6mYKJiwNwdHLN2FrbYW3H2qKJzr48gITEZlRLFDifwa7ubnhr7/+uuu4HJPnhARMsi2LsXcODwgIQFBQkKGHQmYg5MJ1DPxqlwqaZEpu2bj2DJqIiMxQiafqpk+frmqYtm7dqmqcpCh8//79WLt2LRYuXKhes3HjRnTt2hXGTOqx5KaJMolKa+WBSExbfQzpWdloXMMFX49qCx83J15QIiIzVKo+Trt27cLcuXNx6pQUv+agcePG+L//+z906tQJpoZTdVRamVnZmLXuJL7dGaEe927qhdnDWsLZgV0+iIjMNRZgA0w2wKRSuJmcgeeWH8SOM7Hq8fM9G6gb95sjIjI95d4AU06oOZHcLw43yiVzdzbmFsYvDUFEbBIq2dngs2Et0K+5t6GHRUREFcBW16aXUVFR8PT0hKura6ENMGXKTo5nZWXpY5xERrNJ7+Tlh3ArLRO1XCth8ag2aFqTNXJERJZCp8Bpy5Yt2hVzUhROZGnkHwZf7zinapqkKrBdXTfMf6I1PCrnbn5NRESWgTVOrHGie0jPzMa01UexMuSSevx4uzp456GmsLdlU0siInNQ7jVOR44c0fnDpfElkbmIT0rHxGUh2B8RD2srYPqAAIzpVLfQ6WoiIjJ/OgVO0hFc/qLQ1DEVhzVOZE5F4E8tOYCL8cmo7GCLuSNaoVsjT0MPi4iIjD1wko17NQ4dOoSXX34Zr7zyCjp27KiO7dmzB5999hk+/vhj/Y2UqAJtP30Nk34+iFupmfBxq4RvRwehoZdxd8MnIiIjCZx8fe/stzV06FC1gW+/fv3yTc/5+PioruKDBw/Wz0iJKsjSPefxzl9hyMrOQVDdalj4RBu4swiciIhKs+XK0aNH4efnd9dxORYWFsaLSibdCXzm32FYuueCevxI69r4YEgzONjaGHpoRERkJEq8LKhJkyZ47733kJqaqj2WlpamjslzRKboZkoGnlwSrIImKeN7rU9jfDo0kEETERGVLeMkG/kOHDhQTc21aNFCHTt8+LAqGv/7779hKubNm6duLGanC3FJeGpJMMKv5XYC/3x4S/RpVoMXhoiIyqePU3JyMpYtW4aTJ6UZYA4CAgIwYsQIODs7w9Swj5Nl23cuTrUbuJ6cgRpVHPHN6LZoVoudwImILElCefdxKsjJyQkTJkwo7fiIjMKvByJVY8uMrBwE1q6Kr0e1hVcVR0MPi4iITL3GSdoN6CopKQnHjx8vy5iI9EpWy81aewKv/nZEBU39A72xYkJHBk1ERFQ+gdOoUaPw4IMP4tdff0ViYmKhr5EVdW+88Qb8/f1x8OBBXU5LVOGS0jLx9I8hWLT9nHo8uWcDfPVYK1Sy58o5IiK6N52m6iQoWrRoEWbMmIH//e9/aNiwIWrWrAlHR0dcv35d1TpJpmnIkCHYuHEjmjVrpstpiSpU1M0U1Qn8RFSC2mfuk0cDMahlLf4UiIhIf8Xhkk3asWMHzp8/j5SUFHh4eKBVq1bo3r073NzcYGpYHG4ZTkYnYMx3wYhOSIVHZXssGtkWbXyrGXpYRERk7sXhrVu3VjciU7E7PBZPLw3BrbRM+HtWxpIng1C7mpOhh0VERCaoVKvqiEzFH6GX8fLKw6oIvF1dNywe1QauTvaGHhYREZkoBk5klmQGevH2c5i17qR63L+5Nz4b1gKOdiwCJyKi0mPgRGbZbuDdv8OwZPd59fipzn54s38TWFtbGXpoRERk4hg4kVlJzcjC878cwobjV9Wec9P6NcG4LvUMPSwiIjITDJzIbFxPSse4pQcQcuE67G2s1Z5z0tySiIioQgOnL7/8UucTTp48uSzjISp1j6YnvtmnNuqt4mirtk9pX8+dV5OIiCq+j5Ofn59uJ7OywrlzuR2ZTQX7OJm+iNgkFTRdvpEC76qO+OGpdmjo5WLoYRERkaX2cYqIiCivsRGVK+kCPvLb/YhNTEM9D2f8OK49arlW4lUmIiK9YI0TmSypZXry+/1ISM1EE+8qWPpUO1R3cTD0sIiIyIyVKnC6dOkS/vzzT1y8eBHp6en5nps9e3Z5jY2oSDvOXMOEpSFIychSW6d8NyYIVSvZ8YoREZFxBU6bN2/GQw89pOqeTp06pTb0lX3rpFTKEFuxREZGYuTIkYiJiYGtrS2mT5+OoUOHVvg4qOKsPxaNycsPIT0rG10aeGDRyDZwsmfylIiI9M+6pG+YOnUqXnrpJRw7dgyOjo5YtWqVCl66du1qkIBFgqU5c+YgLCwMmzZtwgsvvICkpKQKHwdVjJUHIvHsTyEqaOrbrAa+Gd2WQRMRERlv4HTixAmMHj1aG7SkpKSgcuXKmDlzJj766CNUNG9vb7Rs2VLd9/T0hJubG+Lj4yt8HKR/3+2MwCu/HUF2DjC0TW189XgrONhyCxUiIjLiwMnZ2RlpaWnqfs2aNREeHq59LjY2tsQD2L59OwYOHKjOJe0M1qxZc9dr5s+fr6YGJcPVpk0b7Nixo9BzHThwANnZ2fDx8SnxOMi4fbn5DGb+Habuj73PDx89EghbmxL/8SUiIiqTEv/N06FDB+zatUvd79+/v5q2e//99/HUU0+p50pKptVatGiBuXPnFvr8ihUrMGXKFEybNg2HDh1Cly5d0LdvX1WYnldcXBxGjRqFxYsXl3gMZNw+33gaszeeVvdffLAh950jIiLjboCZlzS4TExMRGBgIJKTk/Hyyy9j586d8Pf3x+effw5fX9/SD8bKCqtXr8bgwYO1x9q3b6+KzhcsWKA91qRJE/WaWbNmqceSAXvwwQcxfvx4VSheHHmtJmOmaXolGSpdml6RYYKmLzafUfen9m2Mp7vW54+BiIiMuwFmXvXq3dkw1cnJSU2j6Yu0OggJCcHrr7+e73ivXr2we/dudV/ivjFjxqBHjx73DJqEBFvvvPOO3sZM5UN+rp9vOqOm6MQb/Rpjwv0MmoiIyLCsyxLUSD8nmTLLeytPUjOVlZUFLy+vfMflcXR0tLov04YynSe1UVIkLrejR48WuypQIkrNTVYEknEHTdP6NWHQRERERqHEGafTp09j7Nix2oxP3r/sZKpNAp3yJuct7LPEfffdpwrCdeXg4KBuZMRB08bT+HLLWfX4zf5NMK7LnSwnERGRSQVOTz75pGpD8Pfff6tWAAWDmvLk4eEBGxsbbXZJQ5pdFsxCkXkETVIE/hWDJiIiMpfAKTQ0VNUdNW7cGPpmb2+v2g9s3LgRDz/8sPa4PB40aFCZzj1v3jx100eGjEqOQRMREZll4BQQEFCqfk1FkRV6Z8/mTsuIiIgIFZxJI8s6dergxRdfVEXfbdu2RceOHVW7AamlmjhxYpk+d9KkSeqmqaQnwwZNn/17GnO35v45mD4gQPVqIiIiMvnASbqDv/rqq/jggw/QvHlz2Nnl31i1pEv6pWll9+7dtY8lUBLSnXzJkiUYPny46tEkncmjoqLU3nhr164tU9sDMh4MmoiIyKz7OFlbWxdbsG1qU18l6d1A5e+LTWfw+abc5pYzBgTgKWaaiIjInPo4bd26FeaANU6G9+3OCG3QJNNzDJqIiMjsMk7mhhknw/g1OBKvrjqi7r/0YEP8X88GBhoJERFZugR9ZpyOHMn9y64gmaaTTXiloJt9kqg4/xyJwuu/5/45mnB/PTzXw58XjIiITEKJAyfpzF1c7yYpFpeC7kWLFqlAiiivrSdjMGXFIWTnAI+3q6P2n9NnLzAiIiKDbrkim/A2aNBAtQWQtgGHDh1S9xs1aoSff/4Z3377LbZs2YI333wTxl7jJK0VgoKCDD0Ui7H3XBwmLgtBRlYOHmpRE+8NbsagiYiIzLvGqV27dnj33XfRu3fvfMc3bNiA6dOnY//+/WrfuJdeegnh4eEwdqxxqhiHI2/gf9/sQ2JaJh5o4okFT7SBnU2pt0okIiIySCxQ4r+5ZAPdwnooyTHN5roynSc9l4jEqehbGP39fhU0daznjrkjWjNoIiIik1TiwEm2Wvnwww+Rnp6uPZaRkaGOabZhuXz5MveSI+VCXBKe+HYfbiRnoKWPK74e3RaOdja8OkREZBnF4VIb9NBDD6F27doIDAxUNSqy0k4aX8rGv+LcuXN49tln9TFeMiFRN1PU9Ny1W2loXMMFS54MQmWHEv+RIyIiMu0+TrK/3LJly3D69GnVMVwyTSNGjICLiwtMDWuc9CMuMQ3DFu1B+LUk1HV3wq8TO8LThassiYjIwvo4icqVK5d5k11DY+dw/UlOz8RTPxxQQVPNqo5YNq49gyYiIrKcjNOff/6Jvn37qh5Ncr84Mo1nSphxKl+ZWdmY8GMItpyMgauTHX6b2An+npXL+VOIiIgMEwvoFDjJxr7R0dHw9PTUbvJb6Mm4ya9Fkz9Kr686ihUHIuFga42fx3dAG99qhh4WERFRxU7VZWdnF3qfKK85m86ooMnaCqrlAIMmIiIyN+xASOXi530X8cXmM+r+u4Ob4cEAL15ZIiKy3MBp3759WLduXb5jS5cuhZ+fn5rCmzBhAtLS0vQxRjJym8Ku4s01uc1PJ/fwx//a390glYiIyKICp7ffflv1a9KQLuFjx47FAw88gNdffx1//fUXZs2aBVPBverKx8GL1/Hc8oNq095hbWvjhQcbltOZiYiITLiPk7e3twqO2rZtqx5PmzYN27Ztw86dO9XjlStX4q233kJYWBhMCVfVlV74tUQ8umA3ridnoHuj6lg8qi23UiEiIpOjl73qrl+/nm8bFQma+vTpo30cFBSEyMjI0o6ZTEzMrVSM/m6/Cppa1K6Kef/j/nNERGT+dA6cJGiKiIhQ92WfuoMHD6Jjx47a52/duqX6PJH5k816n/w+GJeup6iu4N+OCYKTPbdSISIi86dz4CTZJall2rFjB6ZOnQonJyd06dJF+7zUP9WvX19f4yQjkZ6ZjWeWheD4lQS4O9vjh6fawaOyg6GHRUREVCF0ThO89957GDJkCLp27aq2XPnhhx9gb2+vff67775Dr1699DVOMpoGl0ew40wsnOxt8P2TQfB1dzb0sIiIiIwvcKpevbrKNknhlARONjY2+Z6X4nA5Tubr801n8Puhy7CxtsL8/7VGYG1XQw+JiIioQpW4MEWqzgvj5uZWHuMhI7X60CV8ebvB5ayHm6NbI09DD4mIiKjCWWzncPZx0t3+iHi89ltug8tnutXHsCAfvf1ciIiIzKKPk7liH6finY9NwsPzd6m2A32b1cC8Ea1hLZvRERERmQm99HEiy3MzOQNP/RCs7dU0e1hLBk1ERGTRGDhRkW0HJi4LwblrSahZ1RFfj2qLSvb5FwQQERFZGgZOdBeZvZ2+5hj2nIuDs72NanDpWcWRV4qIiCweAye6y6Lt57DiQCSklGnuiNZo4l38fC8REZGlYOBE+aw/FoUP151U998a2BTdG7PtABERkQYDJ9I6HHkDU1aEqvtjOtXF6E51eXWIiIjyYOBEyuUbKRi39ABSM7LRrVF1vNm/Ca8MERFRAQycCIlpmRi7JBjXbqWhcQ0XfPV4K9ja8I8GERFRQRb7tyM7h+fKzMrG//18ECejb8GjsoNaQefiaGfgnw4REZFxYufwEnQLNUfv/HUc3+86D0c7a6yY0BEtfLhxLxERWZYEdg4nXfyy/6IKmsTnw1oyaCIiIroHi52qs3TB5+Mx/Y9j6v5LDzZE3+behh4SERGR0WPgZKEr6Cb+GIKMrBz0b+6N53r4G3pIREREJoGBk4VJTs/E+B8OIC4pHQHeVfDJ0EBYWVkZelhEREQmgYGThe1B9/LKwwiLSoC7sz2+Ht0WTva2hh4WERGRyWDgZEG+2nIWa49Gw87GCgtHtkEt10qGHhIREZFJYeBkIdYfi8bsjafV/fcGN0NQXTdDD4mIiMjkMHCyACejE/Dir3f2oBseVMfQQyIiIjJJDJzMXHxSOsb9cADJ6Vno7O/OPeiIiIjKgIGTGcvIysYzy0Jw6XoKfN2dMPfx1tyDjoiIqAwYOJn5dir7IuJR2cEWX49qi2rO9oYeEhERkUlj4GSmlu29gGV7L0JaNM0Z3hINvVwMPSQiIiKTZ7GB07x58xAQEICgoCCYm73n4vD2n8fV/Vd6N8IDAV6GHhIREZFZsMqRrogWrCQ7IpuCyPhkPDR3J64nZ+ChFjXxxWMt2RmciIionGIBi804maOktEyMX3pABU3Na1XFx49yOxUiIqLyxMDJTEji8JXfDuNk9C14VHbA4lFt4GhnY+hhERERmRUGTmZi4bZz2u1UFo1sDe+q3E6FiIiovDFwMgPbTl/DxxtOqvtvP9QUbXy5nQoREZE+MHAycRfjkjF5+SFIif9jQT4Y0Y7bqRAREekLAycTlpyeiQk/HsDNlAy09HHFO4OacgUdERGRHjFwMuFi8NdXHdUWgy98og0cbFkMTkREpE8MnEzUtzsj8OfhK7C1tsL8/7VGjaqOhh4SERGR2WPgZIJ2n43FB2tPqPvTBwSgnR+LwYmIiCoCAycTc+l6Mib9fBDZOcAjrWtjVEdfQw+JiIjIYjBwMiGpGVl4+scQbWfw9x9uxmJwIiKiCsTAyYSKwaf+fhTHryTAzdkeC0eyMzgREVFFY+BkIpbsPo/Vhy7DxtoKc0e0Qi1XdgYnIiKqaAycTMDec3F475/cYvA3+jVBp/oehh4SERGRRWLgZOSu3EjBpJ8OIis7B4Na1sRTnesaekhEREQWi4GTkReDP7MsBHFJ6WjiXQUfDglkMTgREZEBmUXg9PDDD6NatWp49NFHYU7F4NPXHMPhSzfh6mSHxSPboJI9O4MTEREZklkETpMnT8bSpUthTpbtu4iVIZdgbQV89Xgr+Lg5GXpIREREFs8sAqfu3bvDxcUF5iLkwnXM/Ou4uv9an8bo0qC6oYdERERExpBx2r59OwYOHIiaNWuq+p01a9bc9Zr58+fDz88Pjo6OaNOmDXbs2AFzde1WGp79KQQZWTno39wbE+6vZ+ghERERkbEETklJSWjRogXmzp1b6PMrVqzAlClTMG3aNBw6dAhdunRB3759cfHiRZibzKxs/N/yg7iakIb61Z3x0aMsBiciIjImtoYegARBcivK7NmzMXbsWIwbN049njNnDjZs2IAFCxZg1qxZJf68tLQ0ddNISEiAsfhkwynsPRcPZ3sbLBrZBpUdDP7jISIiImPKOBUnPT0dISEh6NWrV77j8nj37t2lOqcEW1WrVtXefHx8YAzWHY3Cou3n1P1PhraAv6f51GwRERGZC6MOnGJjY5GVlQUvL698x+VxdHS09nHv3r0xdOhQrF27FrVr10ZwcHCR55w6dSpu3rypvUVGRsLQwq8l4pXfjqj747v4oV9zb0MPiYiIiAphEnNBUjResMdR3mMydacrBwcHdTMWSWmZmPhjCBLTMtHez02toiMiIiLjZNQZJw8PD9jY2OTLLomYmJi7slCmSALA11YdwZmYRHhVccDcEa1ha2PUPxIiIiKLZtR/S9vb26v2Axs3bsx3XB536tSpTOeeN28eAgICEBQUBEP5btd5/H0kCrbWVpj/v9ao7mI8mTAiIiIywqm6xMREnD17Vvs4IiICoaGhcHNzQ506dfDiiy9i5MiRaNu2LTp27IjFixerVgQTJ04s0+dOmjRJ3WRVnRSJV7T9EfGYtfaEuv9m/yZo4+tW4WMgIiIiEwucDhw4oDp/a0igJEaPHo0lS5Zg+PDhiIuLw8yZMxEVFYVmzZqpInBfX1+YqpiEVEz6+SAys3MwqGVNjO5U19BDIiIiIh1Y5UihjQXTZJxkhV2VKlX0/nkZWdkY8fVeBJ+/jkZeLlg9qROc7A0evxIREVmshBLEAkZd46RPhqpx+nDdSRU0uTjYYsETrRk0ERERmRBmnCow4/T3kSt47udD6r50Bu/dtIZeP4+IiIjujRknI3Tm6i28ervJ5TPd6jNoIiIiMkEWO1VXkW6lZuDpZSFITs9CZ393vPRgQ0MPiYiIiErBYgOniqpxktr7V1YewblrSfCu6ogvH2vFJpdEREQmijVOeq5xWrw9HB+sPQk7Gyv8+nRHtKpTrdw/g4iIiEqPNU5GYnd4rFpFJ94a2JRBExERkYmz2Kk6fYu+mYrJyw8hOwcY0roW/te+jqGHRERERGXEzot6UtnRFu383BARm4z3BzeHlZWVvj6KiIiIKggDJz2p7GCLeSNaIyElE5XsbfT1MURERFSBLHaqriJW1UmWqaqTnd7OT0RERBWLq+oqeK86IiIiMi5cVUdERESkBxY7VUdERERUUgyciIiIiHTEwImIiIhIRxYbOFXUXnVERERkPriqjqvqiIiILFpCCWIBi804EREREZUUAyciIiIiHTFwIiIiItKRxe9Vl5OTo53fJCIiIsuTcDsG0MQExbH4wOnWrVvqQvj4+Oj9B0NERETGHRNIkXhxLH5VXXZ2Nq5cuQIXFxe1KW9B0q4gODi41BGsBGSRkZHcB8+ElOVnbupM8Xs3pjEbYiz6/kx9nL+8zlnW8/B3tOkJ0tOfd8k0SdBUs2ZNWFsXX8Vk8RknuUC1a9cu8gLZ2NiUOeiR93MDYdNRHj9zU2WK37sxjdkQY9H3Z+rj/OV1zvI6D39Hmw4bPf55v1emSYPF4fcwadKk8vh5kAmx5J+5KX7vxjRmQ4xF35+pj/OX1zmN6WdPFcMYfuYWP1WnT2yuSURkvPg7mkqDGSc9cnBwwFtvvaW+EhGRceHvaCoNZpyIiIiIdMSMExEREZGOGDgRERER6YiBExEREZGOGDiVwPz58+Hn5wdHR0e0adMGO3bsKPK1UVFRGDFiBBo1aqR6RU2ZMqXQ161atQoBAQGqSFG+rl69uiRDIiIiANu3b8fAgQNVA0NpZrxmzZp7Xpdt27ap3+XyO71evXpYuHAhf0fTPTFw0tGKFStU8DNt2jQcOnQIXbp0Qd++fXHx4sVCX5+Wlobq1aur17do0aLQ1+zZswfDhw/HyJEjcfjwYfV12LBh2Ldvn67DIiIiAElJSep37dy5c3W6HhEREejXr5/6XS6/09944w1MnjxZ/WOWv6OpOFxVp6P27dujdevWWLBggfZYkyZNMHjwYMyaNavY93br1g0tW7bEnDlz8h2XoEn6iKxbt057rE+fPqhWrRqWL1+u69CIiCgPyThJ9l5+Pxfltddew59//okTJ05oj02cOFH9I1b+Ucvf0VQUZpx0kJ6ejpCQEPTq1SvfcXm8e/dulJb8z1nwnL179y7TOYmIqPS/fw8cOICMjAz+jqYiMXDSQWxsLLKysuDl5ZXvuDyOjo5Gacl7y/ucRERU+t+/mZmZ6nc+f0dTURg4lTD9W3A35YLHjOGcRERUut+/BY/zdzQVxMBJBx4eHmpH5oKZoJiYmLv+xVISNWrUKPdzEhFR6X//2trawt3dnb+jqUgMnHRgb2+vlqxu3Lgx33F53KlTJ5RWx44d7zrnv//+W6ZzEhFR6X//tm3bFnZ2dvwdTUWyLfopyuvFF19U7QLkfyr5H27x4sWqFYGswhBTp07F5cuXsXTpUu17QkND1dfExERcu3ZNPZYgTPo1ieeffx73338/PvroIwwaNAh//PEHNm3ahJ07d/LiExGVgPyePXv2bL52A/I7183NDXXq1Lnrd7T87pbWBfK7ffz48apY/Ntvv823opm/o6lQOaSzefPm5fj6+ubY29vntG7dOmfbtm3a50aPHp3TtWvXfK+Xy1vwJu/Pa+XKlTmNGjXKsbOzy2ncuHHOqlWr+BMhIiqhrVu3Fvo7V343F/U7+r///stp1aqV+p1et27dnAULFtx1Xv6OpoLYx4mIiIhIR6xxIiIiItIRAyciIiIiHTFwIiIiItIRAyciIiIiHTFwIiIiItIRAyciIiIiHTFwIiIiItIRAyciIiIiHTFwIqIKIbvMr1mzxmKutmzR9MEHH+j9c15++WVMnjxZ759DRLnYOZyISm3MmDH44Ycf1H3ZVV72BQsMDMTjjz+unrO2vvNvM9mJvlq1anBwcDDoFX/77bdVAKfZS1Ifjhw5gm7duuHChQtwcXGBPsXExKB+/frqM/38/PT6WUTEjBMRlVGfPn0QFRWF8+fPY926dejevbvaHHXAgAHIzMzUvq5GjRoGD5rKU0ZGRpHPyeaxQ4cO1XvQJDw9PdGrVy8sXLhQ759FRAyciKiMJBiSoKhWrVpo3bo13njjDfzxxx8qiFqyZEmRU3WvvfYaGjZsCCcnJ9SrVw/Tp0/PF4xIZqhly5b47rvv1O72lStXxjPPPIOsrCx8/PHH6jMlaHj//ffzjefmzZuYMGGCeq5KlSro0aMHDh8+rJ6T8bzzzjvqsYxHbpoxFve+guOR8cr3nbuXd37Z2dlYuXIlHnrooXtOVbq6umo/XwJPec2vv/6KLl26oFKlSggKCsLp06cRHByMtm3bqmsggeq1a9fynUc+a/ny5SX8yRFRadiW6l1ERMWQoKNFixb4/fffMW7cuEJfI9kYCRpq1qyJo0ePYvz48erYq6++qn1NeHi4CsDWr1+v7j/66KOIiIhQAde2bduwe/duPPXUU+jZsyc6dOigApn+/furKcO1a9eiatWqWLRokXpeApDhw4fj2LFj6nybNm1SnyGvudf75Lg4e/asCmxWrVoFGxubQr8vmTK7ceOGCnRK46233sKcOXNUsCjfm0x7SiD3xRdfqCBz2LBhmDFjBhYsWKB9T7t27RAZGammBn19fflnk0iPGDgRkV40btxYBRFFefPNN7X369ati5deegkrVqzIFzhJ9kYyPBJQBQQEqGnAU6dOqeBG6qcaNWqEjz76CP/9958KnLZu3aqCMKn70UwLfvrppyrT89tvv6mMkmRtpB5LMlYaW7Zsuef7RHp6On788UdUr169yO9LMkcSVEnmqrTF3r1791b3ZcpTAqfNmzejc+fO6tjYsWPzZfKEZPs0n83AiUi/GDgRkV5IFkemnooiAYlkViSLk5iYqOqhJLOSlwRUeeuEvLy8VFCSt+hcjknAI0JCQtS53N3d850nJSVFZayKouv7JCgpLmjSvEeCr+K+9+JIcX3e7000b9680O9XQ6b1RHJycqk+k4h0x8CJiPTixIkTRa7y2rt3Lx577DFVbyTZFZka++WXX/DZZ5/le52dnV2+xxKMFHZMMlNCvnp7e6sMVEFST1QUXd/n7OyMe/Hw8FABjGSn7O3ti32t1GsVlPf70wRfBY9pvl+N+Ph49fVeQR0RlR0DJyIqd5qprxdeeKHQ53ft2qWyN9OmTdMek/qcspLidGl7IFNxkq0qjAQzBQMWXd6nKykgF2FhYdr7GvIZGlLgLVmu8iB1WxJcNW3atFzOR0RFYwNMIiqTtLQ0FRBcvnwZBw8eVE0fBw0apNoRjBo1qtD3+Pv74+LFiyrLJFNhX375JVavXl3mn8QDDzyAjh07YvDgwdiwYYOq+ZECcqmnOnDggHqNBEZSYC59nGJjY9X4dXmfriTrI4HYzp0773pOpiYl2ybZOFkhKKRmKy4urkzf944dO7Qr8YhIvxg4EVGZyAo1meaSgESWykuBtgRC0pKgqJVnElhJNuq5555TWRkJUqQdQVnJNJYUjt9///1qRZqsvpMpQQmENPVCjzzyiBqnFJpLkCPL+HV5X0lIMflPP/1013H5XDmvrLjz8fHBs88+i3nz5qmMUVnI9yCrEolI/9g5nIionKWmpqoVf5JRk0yW+mVrZaWyapLVKk///PMPXnnlFbWCUaYaiUi/mHEiIipnjo6OWLp0qZoK1LekpCR8//33DJqIKgj/eUJEpAddu3atkOsqDTGJqOIwcCIiqgCFbc9CRKaHU3VEREREOmLgRERERKQjBk5EREREOmLgRERERKQjBk5EREREOmLgRERERKQjBk5EREREOmLgRERERKQjBk5ERERE0M3/AyjXF0khN7F7AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "### Convert the estimated signal to a digitized signal for direct comparison with a POPS instrument ###\n", "digitized_signal = ppops.digitize_signal(signal)\n", "\n", "fig, ax = plt.subplots(figsize=(6, 4))\n", "ax.plot(diameters, digitized_signal)\n", "ax.set_yscale(\"log\")\n", "ax.set_xscale(\"log\")\n", "ax.xaxis.set_major_formatter(mticker.ScalarFormatter())\n", "ax.set_xlabel(\"Diameter (µm)\")\n", "ax.set_ylabel(\"Signal (digitizer bins)\")\n", "ax.set_title(\"Estimated Signal\")\n", "fig.tight_layout()" ] } ], "metadata": { "kernelspec": { "display_name": "ppops_dev", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.14.2" } }, "nbformat": 4, "nbformat_minor": 5 }