{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial PAMtools" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import importlib\n", "import numpy as np\n", "import pandas as pd\n", "import pamtools.pypam as pp\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Extract 2D histogram" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "scrolled": true }, "outputs": [], "source": [ "h = pp.pamplot.Hist2d.fromfile('data/hist2d.json')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "h.fit_histogram(0, [0.8],[0.1])\n", "h.fit_histogram(1, [0.8],[0.1])" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "scrolled": true }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAATMAAAEiCAYAAACV0jVgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAylklEQVR4nO3de1yMed8H8M+UjiqJKLLLFqW2IpWyIRXrcVhTuMntcEdax5bbYnncWG6xOZa2Za1jIqxD6WHXjXVodXR+dEAWoROpKJWa6/mj18zTNDM1M00z0zXf9+vV66Xr+s0135m7+7O/33X9ruvHYRiGASGEtHFaqi6AEEIUgcKMEMIKFGaEEFagMCOEsAKFGSGEFSjMCCGsQGFGCGEFCjNCCCtQmBFCWIHCjBDCChRmhBBWaKfqAghRB5cvX8a//vUvZGdnw9nZGaNGjYKdnR0AYMyYMSqujkiDQzeaE0137NgxTJkyBTweT7BNR0cHK1euhKurK4VZG0HDTKLRMjMzERQUBB6PhyFDhiA8PBzu7u74+PEjtmzZgsLCQlWXSKREPTOisRiGgZeXF27cuAEXFxesWbMGWlpaqK2txXfffYecnBy4ubkhLS1N1aUSKVDPjGisU6dO4caNGzA0NERoaCi0tOr/79CuXTssWrQI2traSE9Px+XLl1VcKZEGhRnRSAzDYNOmTQCAJUuWoFOnTkL7rays8F//9V8AgA0bNii9PiI7CjOikZKTk5GRkQEdHR3Y2tqKbRMQEABtbW1cvnwZGRkZSq6QyIrCjGikiIgIAIC3tzc6dOggto25uTkGDx4MAPj555+VVhuRD4UZ0Th5eXk4efIkAGDs2LFNtv3yyy8BAEePHsX79+9bvTYiPwozonGOHj2Kuro6ODg4oFevXk22dXBwQO/evfH+/XscP35cSRUSeVCYEY3DD6WhQ4c225bD4WDQoEEAgM2bNyMxMbFVayPyozAjGuXJkye4efMmtLS0BCHVHF9fX2hpaSE7OxuvXr1q5QqJvCjMiEY5ceIEAMDHx0fiif/GOnbsCCcnJwDAn3/+2Wq1kZahMCMahT/E/Nvf/ibT67y8vAAASUlJCq+JKAaFGdEYjx8/xq1bt6CtrQ1/f3+ZXuvh4QEtLS08efIEjx8/bqUKSUtQmBGNwR9i+vr6onPnzjK9tkOHDnB2dhY6DlEvFGZEY/CHmBMnTpTr9V988YXQcYh6oTAjGuHRo0e4c+cOtLS0YGhoKNcUC/5Q886dOzTUVEMUZkQj8IeG/fr1g4mJiVzHoKGmeqMwIxqBPzTkX5WUFw011ReFGWG9nJwc3L17F+3atYOHh0eLjuXh4QFtbW0aaqohCjPCevwh4fDhw2FsbNyiY3Xo0AE+Pj5CxyXqgcKMsB4/dGSdKCsJ/2oohZl6oTUACKtlZ2ejb9++0NHRQWFhoUJuR/Lw8ICFhQXq6urw+PFjWFtbK6BS0lLUMyOs9u9//xsA4OzsrLD7KlNSUuDo6AiAemfqhMKMsBr/XsqWXsVsjH9Vk8JMfVCYEdbKysrCs2fP0K5dOwwcOFChx/b09ISWlhZu3bqFJ0+eKPTYRD4UZoS14uLiANRPlDUyMlLosTt06EBDTTVDYUZYiWEYHD58GED9oiWtgYaa6oXCjLBSSkoKnjx5AgMDA4UPMfn4Q82bN28iNze3Vd6DSI/CjLBSbGwsgPrA0dfXb5X3MDU1ha+vLwDg4MGDrfIeRHoUZoR1Pn78KDhfJs2iJS0xc+ZMAPVhxuPxWvW9SNMozAjrxMfH482bN7CwsBA85aK1cLlcmJqa4vnz57h8+XKrvhdpGoUZYZ3169cDAIYMGQJtbe1WfS99fX1MmTIFALBv375WfS/SNAozwio5OTm4d+8etLS0BKuRt7agoCAAwKlTp1BYWKiU9ySiKMwIq+zatQsA4Orqii5durT6+yUmJiI/Px+2traorq5GdHR0q78nEY/CjLBGaWkp9u/fDwAYNWqU0t6Xw+Fg3LhxAIDo6Gh8+PBBae9N/h+FGWGNyMhIlJWV4ZNPPkH//v2V+t6DBg1Cly5d8Pr1a8TExCj1vUk9CjPCCmVlZdi+fTsAYPLkydDSUu6ftra2NsaOHQsACAsLQ1VVlVLfn1CYERZITEzEvHnzUFpaih49emDQoEEqqWPkyJHo1q0bnj17hqioKJXUoMkozEibV1BQgF9//RVAfa+stadjSKKvry94ftqGDRtQUlKikjo0FYUZadMYhsFPP/2EmpoaODk5YfDgwSqtZ/r06XByckJpaSkWL16s0lo0DT02m7RJ/EV8//Of/yAyMhLt2rVDVFQUunfvruLKgMzMTKxYsQI8Hg9xcXGYNGmSqkvSCNQzI21WVlaWYF5XYGCgWgQZANjb2wsWPZkzZw6ysrJUXJFmaKfqAgiRBr8nxvf06VOEhYWhtrYWnp6emDBhgooqE2/y5Ml4/vw5kpOT4efnh2vXrtHCJ62MhplEbTUOML67d+8iLCwMlZWV6NWrF3744QcYGBgoubrmeXh4wNvbGw8ePEC3bt0QFxen8nN6bEZhRtRW4zB7+/YtDh8+jAsXLgAAHBwcsGrVKoU/EluR3r59i//+7/9GXl4etLS0sHjxYixfvhzm5uaqLo11KMyI2mgcXgzDoKSkBA8ePEBqaipu3LiB2tpaAMCXX36JkJAQ6OrqqqJUmXz48AG7du0SPCKoffv28Pf3R0BAAAYNGoSuXbuquEJ2oDADsHPnTlRXVwtta/y1iPua5Gmj7NepW02PHz+W+LoPHz6gsrISFRUVKC0txcuXL0Xuc7S1tUVQUBAcHBxE3kfdZWRkIDY2VvAd8Jmbm6NHjx6wsrKCqakp2rdvj/bt20NfXx9aWlpN/ixZskRFn0b9UJih/kZhop44HA569uwJZ2dnDB06FDY2NqouqUUYhkFWVhaSkpJw9+5d5OXlif0PgCzHI/UozAghrEDzzAghrEBhRghhBQozQggraPwdAK6urqipqaF5P6TNKS4uhq6uLjIyMlRdilrQ+J5ZdXW1YO4SIW1JbW2tyJQiPk18YofG98z4i15cunRJxZUQIhv+auriaOKCxBrfMyOEsAOFGSGEFSjMCCGsQGFGCGEFCjNCCCtQmBFCWIHCjBDCChRmhBBWUKsw2717N6ZNm9Zkm7dv32LJkiVwc3ODu7s7vv/+e5EH+J0/fx6jRo2Ck5MTuFwukpOTW7NsQogaUJswi42NxY4dO5ptFxoaimfPnuHAgQOIiIjA1atXsXbtWsH+lJQULF26FJMnT8bp06fh6emJkJAQ5Obmtl7xhBCVU3mYFRYWYs6cOdiyZQt69uzZZNvbt28jLS0NP/zwAxwcHODp6Yl169YhPj4ehYWFAIA9e/bAz88P06dPh7W1NZYvXw4HBwccPHhQCZ+GEKIqKg+zBw8eQEdHBwkJCXB2dm6ybUZGBszNzYXWH3R3dweHw8HNmzfB4/Fw69YteHp6Cr1u4MCBSE9Pb5X6CSHqQeU3mvv4+MDHx0eqtoWFhbC0tBTapqurC1NTU+Tn56O8vByVlZWwsLAQatOlSxcUFBRIPG5+fn6TN+3STeiEqD+V98xk8eHDB7FLi+np6aG6uhpVVVUAINKGv58Qwl4q75nJQl9fHzU1NSLbq6urYWhoCD09PQAQaVNdXd3kiteWlpbU+yKkjWtTPTMLCwsUFRUJbaupqUFpaSm6dOkCU1NTGBoairQpKiqihVYJYbk2FWZubm4oKCjAs2fPBNvS0tIAAAMGDACHw4GLi4tgG19qaipcXV2VWishRLnUOszq6upQXFwsOBfm7OwMFxcXLF68GPfu3UNKSgpWr14NLpcr6HkFBQXhf/7nf7B//37k5uYiPDwcWVlZmDFjhio/CiGklal1mOXn58PLywvnzp0DUL+6dVRUFKysrDBjxgwsWrQIQ4YMEZo06+XlhbCwMBw9ehT+/v5ISUnBrl27hKZzEELYR60uAGzatEnodysrK+Tk5Aht69SpEyIjI5s8DpfLBZfLVXR5hBA1ptY9M0IIkRaFGSGEFSjMCCGsQGFGCGEFCjNCCCtQmBFCWIHCjBDCChRmhBBWoDAjhLAChRkhhBUozAghrEBhRghhBQozQggrUJgRQliBwowQwgoUZoQQVqAwI4SwglxhNmvWLJw7d07ssm+EEKIKcj02u66uDt9++y2MjIwwatQoBAQEwMnJSdG1EUKI1OQKswMHDqCgoABnzpzBmTNnEBcXB2tra/j7+2PcuHEwNzdXdJ2EENIkuc+ZWVhYYM6cOfjtt99w9OhReHl54dixYxg2bBjmzJmDP/74Q5F1EkJIkxRyAYBhGPB4PNTW1oJhGBQVFWHBggUYO3YsHj58qIi3IISQJsm91FxeXh7i4+ORkJCAvLw89OjRA5MmTYK/vz+6du2KwsJCzJ49G0uWLMHZs2cVWTMhhIiQK8wmT56Mu3fvQk9PDyNGjMC///1vuLu7C7Xp2rUrRowYgQMHDiiiTkIIaZJcYVZbW4s1a9ZgzJgxMDIyktjOz88PgwcPbvJYPB4PUVFROHHiBN69ewc3NzesXr0aPXr0EGm7c+dOREVFiT1OQEAANm7cCAAICgrCjRs3hPa7u7sjJiamuY9GCCswDKPqEpROrjAzMTHBp59+2mSQAYCdnV2zx4qOjsaRI0ewadMmWFhYYPPmzQgODsbZs2ehq6sr1HbmzJmYPHmy0Lb9+/fj6NGj+Mc//iHYlpOTg7Vr18LPz0+wTUdHR4pPRgg78Hg8VZegdHKF2e3bt8HhcFr85jU1Ndi3bx++/fZbeHt7AwC2b9+OwYMH48KFCxgzZoxQ+/bt26N9+/aC3zMzM3Ho0CGsX78etra2AIA3b97gzZs3cHZ2pikiRGNpYpjJdTVz8ODBSEhIwMePH1v05tnZ2aioqICnp6dgm4mJCezt7ZGent7s69etWwdXV1f4+/sLtuXk5IDD4aBXr14tqo2QtkwTw0yunpmenh4SEhJw/vx5WFtbw9DQUGg/h8PBwYMHmz1OQUEBAMDS0lJoe5cuXQT7JPnjjz9w+/ZtnDlzRmj7w4cPYWxsjHXr1uHPP/+EoaEhRo4ciXnz5okMW/ny8/Ph6+sr8b0uXbrU7GchRJ1QmEmpoKAA/fv3F/ze+GSjtCcfP3z4AAAiIaOnp4eysrImX7t//34MGzYMffv2Fdr+8OFDVFdXw8nJCUFBQcjKykJ4eDhevXqF8PBwqeoipK2jCwBSUtRVQX19fQD15874/waA6upqGBgYSHzdq1evkJqaip9//llk37p167B8+XJ06NABANCnTx/o6Ohg8eLFWLZsGTp37izyGktLS+p9EVbRxJ6ZXOfMpk+fjtzcXLH7srOzMXbsWKmOwx9eFhUVCW0vKipC165dJb7u4sWLMDMzwxdffCGyr127doIg4+vduzcANDt0JYQtWno+uy2SumeWkZEh6LqmpaUhPT0dJSUlIu3++OMP5OXlSXVMOzs7GBkZITU1FZ988gkAoLy8HJmZmZg6dWqTtbi7u6NdO9Hyp02bBisrK8GcMwC4f/8+dHR00LNnT6nqIqStq6ysVHUJSid1mJ04cQLx8fHgcDjgcDj4/vvvRdrww67xlApJdHV1MXXqVGzZsgVmZmbo3r07Nm/eDAsLC4wYMQJ1dXUoKSmBsbGx0DA0MzMT48ePF3vML7/8EmFhYXBycoKXlxfu37+P8PBwzJo1q9l5cYSwRUVFhapLUDqpw2zVqlUYP348GIbBjBkzsHr1atjY2Ai10dLSgomJiWBYJ43Q0FDU1tZi1apVqKqqgpubG/bu3QsdHR28ePECvr6+2LhxIwICAgSvKS4uhqmpqdjjTZ06FRwOBzExMQgLC4O5uTn+8Y9/ICQkROqaCGnr3r17p+oSlI7DyHHZIy0tDfb29qzo6fCnZNAFANLWNPW3a21tLfG8NlvJdTXT3d0dNTU1OHLkCG7cuIHi4mKEhYUhLS0NDg4O9NRZQlQsLy8PPB4PWlqas8yHXJ+0pKQE48ePx4YNG/Ds2TPcu3cPVVVVuHLlCqZNm4bbt28ruk5CiAw+fvyIp0+fqroMpZIrzMLDw1FRUYFz587h9OnTghP/kZGRcHR0RGRkpEKLJITILjk5WdUlKJVcYfbHH3/gm2++waeffip0w7menh5mzpyJBw8eKKxAQoh8EhISVF2CUskVZtXV1RKvJmpra2vkhD1C1E18fDyKi4tVXYbSyBVmjo6OOHLkiNh9Z8+exeeff96iogghLWNjY4Pq6mqx80HZSq4w++abb/Dnn39i3LhxiIiIAIfDQWJiomC1pvnz5yu6TkKUauDAgUI/bU1QUBAAYNeuXRoz7UiuMHN1dcX+/fthYGCAX375BQzD4MCBAyguLsbu3bvh4eGh6DoJaRXiwkpceLW1UHNycoKPjw/q6uowceJEZGdnq7qkVif36kxubm6Ii4tDVVUVysrKYGRkJPQUWELUlaSwkvW1qampCqupNcyfPx8vX75ETk4O3N3dkZiYiCFDhqi6rFbTohl179+/R0lJCerq6lBWVoZXr14JfghRB42Hi4rsXal7T01XVxf/+te/0KdPH7x79w5+fn7YsWMH6urqVF1aq5CrZ5adnY2lS5fi8ePHEttkZWXJXRQhLaWsoBk4cKBa99A6dOiADRs2YPv27bhx4wYWL16Mo0ePYs+ePay7U0euMFu9ejXevn2LZcuWSZyiQYgqqKK3xH9PdQ01fX19LF++HAUFBVi6dCnS0tLg4uKCv//971i1apVMD4ZQZ3LdaN6vXz9s374dw4YNa42alIpuNGcPeYJM0rqu169fF9um4XZJlBVqTf3tJiYmin3NmzdvsHv3bsHdAVpaWpgyZQqWL1/e5qdUyXXOrEePHoLn9xOiDmQNssGDBze5QDV/f+M2zS1qLU8tytSpUyesXLkSW7duhaurK3g8Hg4fPgxHR0cMHjwYhw8fRlVVlarLlItcYfbPf/4TERERSEtLa7MfnLCHLOHRXIgp6xiq1qdPH6xZswbbtm1DQEAAtLW1kZSUhGnTpqF79+5YvHix0NOl2wK5zpn16tVL8JBGcTgcDjIzM1tUGCHSkDbImgofLy8vqd8vKSlJqnbqfmGAr3fv3ujduze++uorXLx4Eb/99htev36NHTt2YMeOHbC0tERwcDACAwNFVkJTN3KdMwsMDERubi5GjRoldrUjAFiwYEGLi1MGOmfWdkkTZJJCrLkAc3R0BFC/fkRjDQOtqXNorR1m8pwza05dXR1u3bqFK1euICUlBTU1NYJ9Tk5O8Pf3B5fLhbOzs9BDJtSBXGHm7OyMjRs3YtSoUa1Rk1JRmLU9LemNiQsxfnA1pXGoSRNobTHMGvrw4QNSU1ORk5OD3377DbW1tYJ9PXv2BJfLBZfLxRdffCF2cSFlk+ucWZcuXZpc15IQRZN10mvjIPPy8hIJMkdHR7FBZv3ZZ0I/4to2PFZbP38miYGBAby9vXH27FkUFhbi4MGD4HK5MDAwwNOnT7Fjxw54e3vD0tISM2fOREJCgkovDMrVMzt+/DhiY2MRERHR5pdvo56Z+pP1BH9j4kKMjx9Wzcl98kTw74a9tKZ6aG29ZyZJVVUV7ty5g+TkZKSnpwstntK+fXuMHDkS/v7+GD16tFLnocoVZsHBwbh9+zYqKythYmIisrAJh8PBxYsXFVZka6IwU2/yTLloSJYgq8u6KvaY2n2HCv7NDzVpAo2tYdZQXV0dMjMzkZKSguTkZKHnp+no6MDHxwf+/v4YN24cLCwsWrUWuQa65ubmGDFihKJrIUSIIoNMUohJCrDyq78BAEyGjhS00e47FNaffYbcJ0/g6Ogo9uJA4/rbwhXNltDW1hYMwYODg9GtWzecPn0ap0+fRmZmJn7//Xf8/vvvmDt3Ljw9PREYGIjAwEB06tRJ4bXI1TNjE+qZqa+WhFlzQcYPKH5oNcdk6EgA/99La9hDU9VQUx16Zk158eIFkpOTkZKSgocPHwq2t2vXDm5ubvD19cWAAQPA5XIV8n4tugRx7do1pKWloby8HB07doSrqytrT4YS9aaIIHsUnyb22L3HuaP86m+CXlrDYSf//aSdf6ZJrKysMHHiREycOBGvX7/GjRs3cOnSJTx58gTJyclITk6Gqakp7t69i3nz5sHc3LxF7yfX1cyamhoEBwcjJCQE+/fvx+XLl/HLL78gJCQEQUFBQnNTmsPj8RAZGYnBgwejX79+mD17NvLy8iS2T0hIgK2trcjPixcvBG3Onz+PUaNGwcnJCVwuV+NWqSGimjrR/yg+TWKQ8fcD0vfiiKjOnTvjq6++QkREBCIjI8HlcmFqaorS0lKsXbsWn376KVatWtWildjlCrOdO3fi5s2bCA8Px71795CUlIS7d+9i48aNuHPnDn766SepjxUdHY0jR45g/fr1iIuLA4/HQ3BwsMRA5D9oLikpSejH0tISAJCSkoKlS5di8uTJOH36NDw9PRESEqJxqztrEkmjAX6vTNw5Mn4wNQyx5BsvRX4kkRSOjWtR5/s0VaVXr16YNWsW9u/fj2XLlsHGxgYfPnzAhg0bYGNjgy1btsg1TJYrzBITE7FgwQJ89dVX0NbWBlA/DuZyuViwYAHOnj0r1XFqamqwb98+hIaGwtvbG3Z2dti+fTsKCgpw4cIFsa95+PAhbG1tYW5uLvTDr2PPnj3w8/PD9OnTYW1tjeXLl8PBwQEHDx6U56MSNdfcaQ1pp15ICi7+9oahJ+miAZFNu3btMHjwYGzbtg0rV65E165dUVRUhJUrV8q1XKXcK5rb29uL3Wdvb4/CwkKpjpOdnY2Kigp4enoKtpmYmMDe3h7p6eliX5OTkwNra2ux+3g8Hm7duiV0PKD+v46Sjkc0hyJCiIaaisfhcODp6YnIyEi4urri48ePiIyMlPkmd7kuAHzyySe4efOmSGgAQHp6umDI15yCggIAEGnfpUsXwb6GysrKUFhYiIyMDBw5cgRv376Fk5MTli5dil69eqG8vByVlZUi81kkHY8vPz9fcGVIHLrSqb6uX7/eZO8s98kTQe9Mu+/QFgca/6omUTxDQ0MsXboUkyZNwqtXr/DixQv06NFD6tfL1TObPHkydu/ejV9++QX5+fn4+PEj8vPzsWfPHuzZswfjx4+X6jj8Wx90dXWFtuvp6aG6ulqk/aNHjwAADMNg48aN2LFjB6qrqzFlyhS8fv1a8DgiaY9H2EGaBybyNb4Syec5qLtcx5Dmvk46bya9K1euAKh/Oq6sVzfl6pkFBgYiMzMTW7ZswdatWwXbGYaBv78/QkJCpDqOvr4+gPpzZ/x/A/Urpou799PV1RXJycno2LGj4I79qKgoeHt749SpU5g4caLgeA1JOh6fpaUl9b7auIY9tKSkJHh5eeH+/ftwdHQU6p01xXNQ9yZP+pPWU1lZiZiYGMGJ/1WrVgllgjTkCjMtLS1s2LABQUFBgnlmHTp0gLu7u8TzWeLwh5dFRUX45JNPBNuLiopga2sr9jVmZmZCvxsYGMDKygqFhYUwNTWFoaEhioqKhNoUFRWha9euUtdF2KdxoJkMHSn2/JekQOs9zl3sELO5uwBI0z58+IDz58/j9OnTKC0tBQAsW7YMK1askPlYLVpqzsbGBlOmTMGcOXMQGBgoU5ABgJ2dHYyMjIRmSZeXlyMzMxNubm4i7Y8dO4aBAweisrJSsO39+/d4+vQpbGxswOFw4OLigrQ04TlDqampcHV1lfHTEVVTxOz5hmHDn7XfcJjYe5y7yGs8B3UX+hEXZA1vPJcGDTX/35gxY9CzZ09cuHABwcHB2L9/P0pLS2FjY4MLFy7ghx9+gJaW7NEkdc/M19cXP/74I+zs7ODj49Pkg9mkvdFcV1cXU6dOxZYtW2BmZobu3btj8+bNsLCwwIgRI1BXV4eSkhIYGxtDX18fQ4YMwZYtW7Bs2TJ88803qKqqwrZt22BmZoaAgAAA9cvSh4SEwN7eHkOGDMHJkyeRlZWFDRs2SPtRCQvwh5oABMNNQLiHxu+d8QNN3MTZxmEn6ZwbaV5RURGuXbuGlStXCv1Hpk+fPlixYgX+/ve/Q0dHR+7jSx1m7u7ughXL3d3dFfaUydDQUNTW1mLVqlWoqqqCm5sb9u7dCx0dHbx48QK+vr7YuHEjAgICYGlpiQMHDmDr1q0IDAwEwzD44osvcOjQIejp6QGov7UkLCwM0dHR2L59O2xsbLBr1y6Ze42EXRoGGiD+ymbvce54FJ8mEmB0BVN+r1+/RmpqKq5fvy40d0xHRwejR4/G3Llz4efnJ1dPrDG60ZxuNFd7LXmeWUtuOG8cYuJuMgdU99RZdbzRnGEY5OXlCW4wb7xQ+Oeff47Q0FBMmDABHTt2VOh7t+hG84qKCrx79w48Hk9kX7du3VpyaEIUTtxwk99Da6731TjIxJFligib8Hg85OTkoLi4GGfOnBFMoQLqTzkNGjQIXC4XkyZNkmnemKzkCrPnz5/jn//8Z5O3HGRlZcldFCENpaamSt07azyJtuG5s8YaBhqfuEm14oJMXK+sqfrZ5uPHj7h//76gB8a/EgnUnwsfPnw4uFwuxo4dq7SZBHKF2ffff4+8vDzMmTMHVlZWChnvEqIoTd0V0PjcWcOAahxsDUl6bHbD92S7qqoq3Lx5EykpKUhPT0dFRYVgX/v27QULnHz55ZcwNjZWen1yhdmtW7ewZs0ahT1UjZDmyNI7A0Qn0QIQurrJJynYxJG0QpMql5trTWPGjEF5eTnOnDmDU6dO4ffffxda9NvCwgJcLhf+/v7w9vYWufNG2eQKs/bt27f4QWqEyKolgQaIH3LKM+lV2nUz26ra2lrcvn0bsbGxiI+PF1px6bPPPoO/vz8CAgLg4eGhVqMyucJs3LhxOHToEDw8PASP3iFEGRQRaIBsq5g3fF3jY7NJXl4ezp8/j6tXr6K8vFyw3dbWFpMmTUJAQACcnJzUbvFfPqnDrOHtBbW1tbh+/TqGDx8OJycnkfseORwOwsLCFFclIQ20NNAA8aEmy6OvpQmytjDErKurQ0ZGBhITE3Hnzh3BdlNTU8yYMQNTp07FgAED1DbAGpI6zBr/D8N/zM69e/dE2raFD07aNnkCDRCdhybPs/vZEGR1dXW4cuUK4uLiBI/H4nA4GDt2LObMmYPhw4erxSrlspC62suXL7dmHYQoRXPPP2vutW0dwzBISUlBTEyMYK0NIyMjzJ07F3PnzkWvXr1UXKH85I5e/iXa+fPnAwAyMzOxe/duzJ49G59//rnCCiREEll7Z3ySemmS2slTlzp68+YNIiMjcevWLQD1ITZhwgSMHj0aEyZMUHF1LSdXmF29ehXz58+Ho6OjIMw4HA6ePn2KKVOmYN++ffSUCqIU8gYawI6elrSuXbuG6OhoVFRUQF9fH0uWLMG3334LU1NTVZemMHKF2c6dOzF69Ghs2rRJsK1v376Ij4/H8uXLsW3bNhw5ckRhRRLSFHE9IXrkTj2GYXD48GEcP34cQP1juxISEtC3b18VV6Z4ck0Syc3NBZfLFXuin8vlIjs7u8WFEdIS6jrUU7Zjx44JgmzChAnYvHkzK4MMkDPMjI2N8ddff4ndl5eXB0NDwxYVRYgiqCLQ1ClEL1++jNjYWABASEgIZsyY0eauUMpCrk82fPhwREREwNLSEsOGDRNsv379OiIiIjBixAiFFUiIulGnwJKkpKQEu3fvBlA/R1QT5n3KFWaLFy/G/fv3MXfuXOjo6AiWWa+trYWzszOWLFmi6DoJkQs/eCSdQ2tuv6T26u7gwYOorKxE7969sX79elWXoxRyhZmRkRHi4uJw9epV3Lx5E2VlZTA2Noarqyu8vb3V6n4tQoDmQ6jx/obh1lYCrKFr164BAL7++muNueVQ7gG0lpYWhg0bJjTMJIQt2mKANVRbW4s+ffpo1ChJpnsz582bhx49ejS7DBTdm0mI6sl6M31bJ9O9mTNmzBD8uyl0byYhqifNautsIte9mXSfJiHqjcPhtOrz9tVRiyadlJeX486dO3j37h3MzMzg6OgIIyMjRdVGCJGTmZmZYPlFTSF3mP3888+Ijo4Weoyurq4uvv76a8H9moQQ1TAxMVF1CUonV5idPHkS27Ztw4QJE/DVV1+hc+fOKC4uRnx8PKKiotCtWzf4+/srulZCiJRUsaCIqskVZgcOHEBgYCDWrFkj2PbZZ59h4MCB0NfXx6FDh2QKMx6Ph6ioKJw4cQLv3r2Dm5sbVq9eLXHM/+jRI2zevBl3796FlpYW3Nzc8N133wnW6qyrq0P//v1RXV0t9LoFCxZg4cKFcnxiQtoWGxsbjBkzRtVlKJVcs1ufPXsGPz8/sft8fX3xpJlVbhqLjo7GkSNHsH79esTFxYHH4yE4OBg1NTUibd++fYugoCDo6+sjJiYGe/bsQUlJCYKDgwXh9fTpU1RXVyM+Ph5JSUmCn5kzZ8r+YQlpg1S9UpIqyBVmXbt2xatXr8Tue/HihUwXAWpqarBv3z6EhobC29sbdnZ22L59OwoKCnDhwgWR9hcvXkRlZSXCw8PRp08ffP7559i8eTNyc3MFD53LycmBkZER7OzsYG5uLvhp3769PB+XENIGyBVmPj4+iIiIEHn+/927d7Fz5074+PhIfazs7GxUVFTA09NTsM3ExAT29vZIT08Xae/p6Yno6Gjo6+sLtvFvn+KvKJOTkwNra2uZPhMhpG2T65zZwoULcePGDUyaNAndu3dH586d8fr1a7x8+RLW1tYy3ULBX0zB0tJSaHuXLl0E+xqysrKClZWV0Laff/4Z+vr6cHNzAwA8fPgQtbW1mDVrFrKzs9G1a1fMmDED48aNE1tDfn4+fH19JdZ46dIlqT8PIepAEyeuy32j+a+//oqTJ08iPT0dZWVlcHR0xMyZMxEQECDUa2oOf4HRxmN8PT09lJWVNfv6mJgYHD58GKtWrYKZmRmA+gsEPB4PoaGhsLCwwNWrV7FixQp8/PiRFc86J4SIkivMzpw5g6FDh2LKlCmYMmWK0L7i4mLExMRg9uzZUh2LH3w1NTVCIVhdXS2yHmdDDMMgIiICP/30E+bOnYtp06YJ9iUmJqKurk5wjszOzg6vXr3C3r17xYaZpaUl9b4Iq2hiz0yuc2YrVqwQLFPVWFZWFiIjI6U+Fn94WVRUJLS9qKgIXbt2Ffuajx8/YunSpdi1axdWrFiBRYsWCe3X19cXOdnfp08fscNWQgg7SN0zCwkJQW5uLoD6XtH8+fPFXv598+YNPvnkE6kLsLOzg5GREVJTUwWvKy8vR2ZmJqZOnSr2NcuWLcN//vMfbN26FaNHjxbaV15eDj8/P3z33XcICAgQbL9//z569+4tdV2EtGWa2DOTOszmzJmDEydOgGEYvHz5Evb29oJzVHxaWlowMTERCpHm6OrqYurUqdiyZQvMzMzQvXt3bN68GRYWFhgxYgTq6upQUlICY2Nj6Ovr49SpUzh37hyWLVsGd3d3FBcXC45lbGwMExMTeHh4YPv27ejUqRM+/fRTXLhwAQkJCYLHCBNC2EfqMHNxcYGLiwuA+tTnP9sMAMrKyvD8+XP07NlTrtsoQkNDUVtbi1WrVqGqqgpubm7Yu3cvdHR08OLFC/j6+mLjxo0ICAhAYmIiACA8PBzh4eFCx+G3CQsLw86dO7FmzRq8efMG1tbWiIyMlHsla0KI+uMwDMNI2/jevXuIjo7GyJEjweVyAQCHDx/G5s2bUVNTAz09PSxcuBCzZs1qrXoVjj8lgy4AkLamqb/dadOmISYmRtklqZTUFwCys7Mxbdo0ZGVlCZaSu3//PjZs2IAePXpg586dmDdvHrZv346LFy+2WsGEECKO1MPM3bt3w87ODgcOHBBMmTh06BAAYMuWLbCzswMAvH79GjExMRLv3SSEtD5NvAAgdc8sPT0d06ZNE5r7lZSUhB49egiCDKh/7nhmZqZiqySEkGZIHWalpaWwsLAQ/J6bm4u3b9+KrDdoYGAg9mkXhBDSmqQOM1NTU7x580bwe0pKCjgcjtAN4kB9yDWeskEIUS4aZjbB3d0dx48fB8MwqK2txcmTJ6Gnpyc03aGmpgaxsbGCKRyEEKIsUl8AmDt3LiZNmgQ/Pz8wDINXr15h/vz5gnllJ0+eRGxsLP766y+R+V+EEOXSxJ6Z1GHWu3dvHD9+HPv27cObN28we/ZsBAYGCvbv2LED7dq1w48//oi+ffu2SrGEECKJTE/NsLGxkbhS+a+//gpzc3PBgxIJIapDPbMWkPSEC0IIUQbqRhFCWIHCjBAW0sRhJoUZIYQVKMwIYSHqmRFCSBtFYUYIC1HPjBBC2igKM0IIK1CYEcJCNMwkhJA2isKMEBainhkhhLRRFGaEEFagMCOEhWiYqQI8Hk+w2ni/fv0we/Zs5OXlSWz/9u1bLFmyBG5ubnB3d8f333+PDx8+CLU5f/48Ro0aBScnJ3C5XCQnJ7f2xyCEqJjKwyw6OhpHjhzB+vXrERcXBx6Ph+DgYIkrPIWGhuLZs2c4cOAAIiIicPXqVaxdu1awPyUlBUuXLsXkyZNx+vRpeHp6IiQkBLm5uUr6RISoHvXMlKympgb79u1DaGgovL29YWdnh+3bt6OgoAAXLlwQaX/79m2kpaXhhx9+gIODAzw9PbFu3TrEx8ejsLAQALBnzx74+flh+vTpsLa2xvLly+Hg4ICDBw8q++MRQpRIpWGWnZ2NiooKoeXqTExMYG9vj/T0dJH2GRkZMDc3h7W1tWCbu7s7OBwObt68CR6Ph1u3boksfzdw4ECxxyOErTSxZ6awx2bLo6CgAABgaWkptL1Lly6CfQ0VFhaKtNXV1YWpqSny8/NRXl6OyspKocWKmzoeX35+Pnx9fSXuv3TpUrOfhRCiWirtmfFP3Ovq6gpt19PTQ3V1tdj2jds2bF9VVSXT8Qgh7KHSnpm+vj6A+nNn/H8DQHV1NQwMDMS2F3dhoLq6GoaGhtDT0xMcr/F+ccfjs7S0pN4XYRVNHGaqtGfGHzIWFRUJbS8qKhK72pOFhYVI25qaGpSWlqJLly4wNTWFoaGh1McjhLCHSsPMzs4ORkZGSE1NFWwrLy9HZmYm3NzcRNq7ubmhoKAAz549E2xLS0sDAAwYMAAcDgcuLi6CbXypqalwdXVtpU9BiPrRxJ6ZSoeZurq6mDp1KrZs2QIzMzN0794dmzdvhoWFBUaMGIG6ujqUlJTA2NgY+vr6cHZ2houLCxYvXoy1a9eisrISq1evBpfLFfS8goKCEBISAnt7ewwZMgQnT55EVlYWNmzYoMqPSghpZSqfNBsaGooJEyZg1apVCAwMhLa2Nvbu3QsdHR3k5+fDy8sL586dA1D/X5uoqChYWVlhxowZWLRoEYYMGSI0adbLywthYWE4evQo/P39kZKSgl27dglN5yCEsA+HYRhG1UWoEn9KBl0AIG1NU3+7CxYsQFRUlLJLUimV98zaMl9f3ybnpxH6jppD34/iUJgRwkKaeAGAwowQwgoUZoSwEPXMCCGkjaIwI4SwAoUZISykicNMld4BoA6KiopQV1cn1+Xx/Px8AKBL602g76hpLfl+8vPzoa2tLXZfUw9WYCuNnzTr6uqKmpoamJubq7oUQmRSXFwMXV1dZGRkqLoUtaDxYUYIYQc6Z0YIYQUKM0IIK1CYEUJYgcKsCa2xQDHbyPodJSQkwNbWVuTnxYsXSqxaNXbv3o1p06Y12UYT/4YUhcKsCYpeoJiNZP2OcnJy4O7ujqSkJKGfxqtusU1sbCx27NjRbDtN/BtSGIaIVV1dzfTv35+JjY0VbCsrK2OcnJyYs2fPirS/desW06dPH+bx48eCbdevX2dsbW2ZgoICpdSsbLJ+RwzDMMHBwcz69euVVaLKFRQUMF9//TXTr18/ZuTIkczUqVMlttXEvyFFop6ZBIpeoJiNZP2OgPqemSY99ffBgwfQ0dFBQkICnJ2dm2yriX9DiqTxdwBIougFitlI1u+orKwMhYWFyMjIwJEjR/D27Vs4OTlh6dKl6NWrl1JqVjYfHx/4+PhI1VYT/4YUiXpmEih6gWI2kvU7evToEQCAYRhs3LgRO3bsQHV1NaZMmYLXr1+3fsFqThP/hhSJemYSKHqBYjaS9TtydXVFcnIyOnbsKLgROioqCt7e3jh16hRCQkKUU7ia0sS/IUWinpkEil6gmI1k/Y4AwMzMTOiJDgYGBrCyskJhYWHrFdpGaOLfkCJRmEmg6AWK2UjW7+jYsWMYOHAgKisrBdvev3+Pp0+fwsbGRik1qzNN/BtSJAozCRouUHzp0iVkZ2dj8eLFQgsUFxcXo6qqCgCEFii+d+8eUlJSRBYoZhtZv6MhQ4aAx+Nh2bJlePToEe7fv4+FCxfCzMwMAQEBKv40ykd/Qwqm6rkh6qy2tpYJDw9nPDw8mH79+jGzZ89m8vLyGIZhmLy8PKZPnz7MyZMnBe1fv37NLFy4kOnXrx8zcOBAZs2aNUxVVZWqylcKWb+j//3f/2WCgoKYAQMGMC4uLszChQuZV69eqap8pVq+fLnQPDP6G1IsegQQIYQVaJhJCGEFCjNCCCtQmBFCWIHCjBDCChRmhBBWoDAjhLAChRkhhBXoRnMi5LvvvsPp06cl7o+IiMCVK1fEtjE0NISVlRUCAgIQFBQk9TE7d+6MP//8Ezt37kRUVFST9XXv3h2XL1+W4pMQTUNhRkSYm5tLDJWePXviypUrIm0YhsHr168RFxeHTZs2QU9PD1OmTJHqmDo6OgCAiRMnYvDgwYLtJ06cwK+//opjx44Jtol7RA4hAIUZEUNXVxf9+vWTq423tzf8/Pxw6tQpoTCT5pgWFhawsLAQ/H79+nUAaPZ1hAB0zowomI6ODgwMDIQe80OIMlDPjIhVW1srsk1bW1sopBq24fF4KCoqQkxMDP766y+sXr1aqmOKOy4h8qAwIyJevnwJBwcHke1LliwRPA1WUpuePXtizZo1CAwMlOqYALBs2TLMmjVLAZUTTUZhRkSYm5vjp59+Etne8HxWwzbl5eWIjo7G8+fPsWnTJvTv31/qYwKiC6IQIg8KMyJCV1cXjo6OMrVxcXHB+PHjMXv2bJw4cUJktSVpjklIS9AFAKIQBgYG2LRpEyoqKrBixQrQY/KIslGYEYVxcnLC3/72N9y+fRtnzpxRdTlEw9AwkyjUokWLcP78eWzduhXDhw+HkZERgPpVhu7cuSPxdba2tmKXpyNEWhRmRKE6duyIb775BuvWrcOPP/6I5cuXAwCKi4sxadIkia87c+YM+vbtq6wyCQvRGgCEEFagc2aEEFagMCOEsAKFGSGEFSjMCCGsQGFGCGEFCjNCCCtQmBFCWIHCjBDCChRmhBBWoDAjhLAChRkhhBUozAghrPB/r75oHa0sse0AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h.plot2Dhist(style='contour', imgOffset=0, PAMtools_fit=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Export FCS fit results" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "fcs = pp.pamplot.FCS.fromfile('data/fcs.txt')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQwAAADjCAYAAABw6SpwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA4kElEQVR4nO3deVxU5f7A8c8wwLANArKpiIqyiIqCLJpoSmq5lbZaomYulXa9eX+lLeY1SzO1xfJqWlqWa5qmphZiKuIVRCV3FneQXfYdZs7vj7lMjoACggzwvF8vXuk5zznz5QRfn+08j0ySJAlBEIQaMGjsAARBaDpEwhAEocZEwhAEocZEwhAEocZEwhAEocZEwhAEocZEwhAEocZEwhAEocYMGzsAfePr60tpaSl2dnaNHYog1Ep6ejrGxsacPHmywT5DJIy7lJSUoFKpGjsMQai18vJyGnritkgYd7G3twfg4MGDjRyJINTOY4891uCfIfowBEGoMZEwBEGoMZEw7iMjO49zlxPIyM5r7FAEodGJPox7OBB5npXbDyJJEjKZjOnPPsaQgO6NHZYgNBpRw6hGRnaeNlkASJLEqu0HRU1DaNFEwqhGzPWkSkNUakkiOSO7cQISBD0gEkYVikpKWbZhf6XjBjIZbWytHn5AgqAnRMKoQm5BcaVjBjIZrz/7GLZWykaISBD0g0gYVao8W27q6EdxbN3qvn0YYlRFaM7EKEmVZJWOrNl5GAnuOVoiRlWE5k6vaxirV69m/Pjx9yyTlZXF//3f/+Hn54e/vz8ffvghRUVFOmX279/P8OHD8fLyYvTo0Rw/fvye97Q0N8FApps0KuockiSxcltopRqEGFURWgK9TRgbN27kyy+/vG+5mTNncuPGDX744QeWL1/OkSNHmD9/vvZ8REQEb7/9NmPHjmXnzp307duXadOmceXKlWrvaaowZs37r/BW8LAqz0vA97vDtMkgIzuPY2fixKiK0OzpXZMkNTWVf//730RGRtKxY8d7lo2OjubEiRPs27ePzp07A7BgwQKmTJnCv/71LxwcHPj2228ZPHgwEyZMAGDOnDlER0ezfv16FixYUO29ba2UtLIwq/b8sbPx/PfcZQb29uDwqZhq3xJMSM2kja2V6CwVmgW9q2FcuHABIyMjdu/eTc+ePe9Z9uTJk9jZ2WmTBYC/vz8ymYxTp06hVqs5ffo0ffv21bkuICCAqKio+8aiMDZCVrk7Q0uSJA6dvHTPV4rX7DzE1I/XciDy/H0/TxD0nd7VMIKCgggKCqpR2dTUVNq0aaNzzNjYGCsrK5KTk8nNzaWwsBBHR0edMvb29qSkpFR73+TkZNy7eWFoYYOZlR2mVnaYWdujMG+F7F4ZpBoSsGr7QbzdO4iahtCk6V3CqI2ioiKMjY0rHVcoFJSUlFBcrJlPcXeZivPVkSTIz7kN2bfJTozXHpcbKbCwa4fSzgkLOyfMbRwxkMtrFGtFf4ZIGEJT1qQThomJCaWlpZWOl5SUYGZmhkKhAKhUpqSkBFNT02rva2tnR7s+o8hPv0VRdjqF2WkU5dxGVVZCTtJVcpKuAmBgaISlYwes2nbBqq0LRqbm1d5ThqaJk5GdR8z1JAA8OrYVCURoUpp0wnB0dCQ0NFTnWGlpKdnZ2djb22NlZYWZmRlpaWk6ZdLS0nBwcKj2voZyOUrbtli0/ru5I6lVFGalk5eeSF56IvnpiZSXFJGdeJnsxMsAWNg5YdupGzbO7siNFDr3lIDZX23ROSYDpj83WMzVEJoMvev0rA0/Pz9SUlK4ceOG9tiJEycA6N27NzKZDB8fH+2xCpGRkfj6+lZ7XwMDzaSrirkYBjIZb7zwONMnPIujhy+u/UfTa8wMPB+fQNvuj2Bmo0k++emJXD/xB3/tXMmV//5GXlrCPTtEK/o2xFwNoaloUjUMlUpFZmYmSqUSExMTevbsiY+PD7NmzWL+/PkUFhYyb948Ro8era1BTJo0iWnTpuHp6cmAAQP45ZdfuHTpEgsXLrznZw0J6I63eweSM7IrDYuu3xuOTCbD3MYBcxsH2vXoR2lhHrevXyTj2nmKczPJvHGJzBuXMLN2wNHDFxtnD2QGlfOz6NsQmpImVcNITk4mMDCQffv2AZpp2itWrMDJyYmJEyfy5ptvMmDAAJ2JW4GBgSxatIjNmzczZswYIiIi+Oabb3SGYqtja6WkR5f2Or/MYwb58t3cybwVPIxHvd21w64KMyUvvzKZHiMm03VoMHadvZDJDSnMSuXq8b2c3/8DmTdjK9U4xBuwQlMikxp6XfImpmLl5ZquGp6RnadTCzkQeZ6V20KRgLKSQtLjz5ASexJVqWbExszKHqdeA2jVppO2D6O6mowg1EZtf3brQiSMu9THQ4+7maLTwVleWkJq7ElSYk6iLteM2Fg5ufLpJ4tIyCwkJOL8fV9sE4T7eRgJo0k1SZoKN2dHnnrUR/t3Q2MF7Xr0o+eT03Bw7w0yGdmJ8bz+8ov8sG4tarUaEC+sCfpPJIwGMqq/d6VZoUYKU776fBnrftqC0r49alU5iWfCuBS6ieLcTEC8sCboN5EwGoitlbLS0Oz05wYT2MuNpHw17kEv0ClgGHIjYwpuJ3Ph9/Wkxp1GkiTRCSrorSY1rNrUVDU0m5Gdx64jp5HJZNi6dEfp4Mz1yN/JTb3BzVMHyUtL4NeD3fHo7AyI2aCCfhEJo4HZWil1fuHvbm4ozC1xG/QcafHRJEQfIishjkXvvUmXwKcws9bs8zpxRCBjBlU/0UwQHhbRJHnI2thaVerbkMlkOLj54DH4JYzNLCnJz+bigY1kXNO8Er9+bzg7D51sjHAFQYdIGA/Z3X0bd7Jo3YZuT0ygVZtOSKpyrkXsJ+GvI0iSxPq94cTdrP6VfEF4GMQ8jLs8jLFs0J3wdTQ6lvV7w7XnJEni1rljJF/QrD1q5dQFlz4jkBsZM0O8rCZUQ8zDaMbunHY+ZpAvE0cEapsqMpkMJ69AXPqOQGYgJzvxMjGhmykpyOU/VSxALAgPi+j01BNjBvnS39ud5IxsFMZG/HwggpOAwqIV8WG/UpidxqWQDXQZMIbPN+5n0YznGztkoQUSNQw9UlHrcHN25LVnNNVLC9t2eD4ejGkrW8qKC4g9uIVjYYf55Ps9oqYhPHQiYegpWyslM54bDIDCvBVdh4yjVZtOqFXlXA7fxa5ftjD5o+/4fk+YSBzCQyOaJHpsSEB3WlmYsej73ciNjHEd8DQ3Tx8kLf4vEqIPU5yXza+SxK4jpxka0J3nhwSISV5Cg6pzwsjKyiI0NJTjx4+TmJhIXl4e1tbWtG3blgEDBjBw4EAsLS3rM9YWyb+bC329unD87GVkBgY49x6MwsKahOhDpF/+i9KCHDr3G0VI5HlCIs/z1KM+jOrvLRKH0CBqPayamZnJqlWr2L59OyqVis6dO9OuXTtMTU3Jzc0lJSWF+Ph4jI2NGTt2LFOnTqV169YNFX+9e1jDqrW1Zuch9h07o/17VmI8V/+7F7WqDFMrO1wHPI3CXJOgxWvyLZPerYexf/9+PvroI7y8vBg1ahRBQUFVrr6dn59PWFgYP//8MzExMcybN4/hw4fXa+ANRV8TBsDOQyd15msU3E4hPmwHZcUFGJma4zrgGcz/t76oDPh27mRR02hB9C5hjB8/nvfee4+uXbvW+APOnTvH4sWL2bhxY50CfNj0OWEA2m0Koi5c40h0DCUFucQf+YWinAwM5EZ07jcSq3ZdAOjn5cqkJweIpNFC6F3CqCBJUp12AGsK9D1h3KmixqEqK+Fy+G5yU64D4OwTpFmoB9E8aUn0dqanj48PR44cqe9YhFoaM8iXJTPHIjdS4Pro09h10exFe/P0n9w4GYparUKSJFaK2aFCPalTwigqKiI3N7e+YwFArVbz1Vdf0b9/f3r16sXUqVNJSEiosuzXX3+Nu7t7lV/vvvuuttykSZMqnR8/fnyDxP+wuTk7MnFEIAYGcjr4DqF9r4EApMVHE3twK6WFeUjAttDIRo1TaB4aZB5GZmYmJSUllTZKromVK1eyadMmFi9ejKOjI0uXLmXKlCns2bOn0h6pr7zyCmPHjtU59v3337N582Zefvll7bHY2Fjmz5/P4MGDtceMjIxqHZu+qlgrY/3ecBy7+qFQWnMtYh/5Gbe48PuPdH5kJH9EQFFxGf7dXcSiPEKd1XmmZ2xsLDdu3KhyZ69Dhw7VeAf2O5WWlrJu3TpmzpzJwIED8fDw4IsvviAlJYWQkJBK5c3NzbGzs9N+paen8+OPPzJv3jzc3d0BuH37Nrdv36Znz546Za2srGodnz4bM8iXt4KHAWDt1AXPx8djamVHeUkhsYe3kXQhgiPRMSzbsJ+pH6/lQOT5Ro5YaIrqXMNYu3Yta9euRaFQ4OLigqurKy4uLtjY2PDLL7/g6OhY63vGxMRQUFBA3759tccsLS3x9PQkKiqKkSNH3vP6BQsW4Ovry5gxY7THYmNjkclkdOrUqdbxNDUeHdsik8mQJAkTpTWeQ8Zx42QoGdfOc+vsUfIzbuHSZziGClP+sy2UsnIV/t1cRG1DqLE6J4wZM2bg5OREXFwc8fHxREZGsmvXLgDkcjnz5s2r9T1TUjQLxNzdlLG3t9eeq86hQ4eIjo7m119/1TkeFxeHUqlkwYIFHDt2DDMzM5544gmmT59eqYnT1FUszrNq+0HUkoSBoRGd+gzDwq4dN06GkpN0lQu/r6dzvyexsG3Lmp2HWLPzkFgCUKixOieMDh06MGrUKJ1jhYWF3L59G2traywsLGp9z6KiIoBKv8gKhYKcnJx7Xvv9998zaNCgSnNE4uLiKCkpwcvLi0mTJnHp0iWWLFlCUlISS5YsqfJeycnJ2iGqqujzkOudCw+fib/J9oNR2HX2wtzGgcvhuynJzyYmdDNOvQbg4O6LTCZj/d5wriWlM2FEoKhtCPdUr52eZmZmmJmZ1fl6ExMTQNOXUfFngJKSkipnlFZISkoiMjKSNWvWVDq3YMEC5syZQ6tWrQBwc3PDyMiIWbNmMXv2bGxtbescr76qWHi4R5f2PNHXi9+ORvPrkdN0e2IC1078QdbNWBKiD5OXlkinPsMwNDYhLDqWo9GxTBcregn3UKeEMX/+/BptZlxbFU2RtLQ0nJ2dtcfT0tK0nZhVCQ0NxcbGhn79+lU6Z2hoqE0WFVxdXQFNE6iqhNGmTRu9rkXUhq2VkpdHDeCRnm7M/moLnR8ZRbp9e26ePkT2rcuaUZR+o7Bo3QYJWLX9IN7uHURNQ6hSnUZJxo4di6enZ33HgoeHBxYWFkRG/j1nIDc3l4sXL+Ln51ftdSdPnsTf3x9Dw8r5b/z48TpzMkAzXd3IyIiOHTvWW+z6rmL7RplMhr2rN12HvITCohWlBTnEhG7SbqIkdl4T7qVWCWP27NlkZGTU6gNSUlL4v//7vxqVNTY2Jjg4mGXLlnHw4EFiYmKYNWsWjo6ODB06FJVKRXp6OsXFxTrXXbx4EQ8Pjyrv+fjjj7Nr1y42b95MQkIC+/btY8mSJUyePLlO/SxN2Z3bN5rbOOL5+ASsnFyR1GpunjrIjZOhSGo1lxNSGzlSQV/Vqkni4eHByJEjGTlyJE8++SReXl7Vlj179ixbt24lJCSEV199tcafMXPmTMrLy5k7dy7FxcX4+fmxdu1ajIyMSExM5LHHHuOTTz7h6aef1l6Tnp5e7byK4OBgZDIZP/30E4sWLcLOzo6XX36ZadOm1Tim5uLuURRDYxO6BD5FSkwUiX8d0ayvUZjLOlU5RSWl+Hq64OZc++Fxofmq9ctnV65cYdmyZRw+fBh7e3t69OiBk5MTpqam5OXlkZycTHR0NFlZWQwcOJBZs2bh5ubWUPHXu6b08lld3bnFQcz1JJZt2E9mQhxXj+9FUpVjZm2P66PPYGxqwSDfrvxz7OONHbJQAw/jZ7fWnZ6dO3dm1apVxMXFsWfPHiIjIzl16pR2xa127drx4osvMnTo0Ht2VAqN587tGz06tgXApr0bxqYWxIftpDArjZiDW/AIeoFDJy8x7JGeoqYhAGIjo0paQg3jbt/vCWPXkdMAFOdnE/vnVkoLclFYWOHx2Fj69e7B+5NHN26Qwn3p7evtQvNyZ2eoiYUVHkFjMTZvpZnkdXALx06dZ8O+Y40cpaAP6jQPY//+/QAMGzYMtVrNkCFDdM6PGjWKN99884GDEx6OuztDFRat8HjsBWIPbqUkP5u4I9vZYqh5uzd4eOW5LkLLUasahkqlYsaMGfzrX/8iLCwM+N8+oLdu4erqir+/P46Ojnz33XfcvHmzQQIWGsaQgO6sef8V7RuvCvNWuAc9j5GJOUXZ6Vw+upOfDxwXNY0WrlY1jJ9//pmwsDCWL1/O0KFDdc794x//oFu3bhQXF/P444+zZcsWZs+eXa/BCg3L1kpJYC930rPyWL83HIWFFa6PPkPMwS3kpSVwLfJ3thnIAVHTaKlqVcPYtWsXL7zwQqVkcScTExOeeeYZjh0T/xI1VWMG+fJskGZmrbmNA136PwUyGZk3LpEac5Ltf0ax89DJRo5SaAy1ShiXL19mwIAB9y3n4+MjmiRNXPDwftqk0cqxI84+mgWREs4cISf5Oj/uOybWCW2BapUwysvLK701KpfLCQkJ0b7QVXHMwEAMwDR1dyYNe1dvbDt1B0niyn/3UJyXLd45aYFq9Vvt4ODAtWvXKh13dnbWWcMiLi6Otm3bPnh0QqOrSBoymYwOfkMwb90GVWkxV4/v5fSlyj8LQvNWq4QRGBjI1q1bUavV1ZYpKytj+/btDBo06IGDE/RDRdIwkBvS+ZGRGBgak59xixUrvhajJi1MrRLGuHHjuHLlCm+++SZZWVmVzhcWFjJnzhySk5N58cUX6y1IofH1dNOsT6KwsKKjn2beTdKF43y/dafoAG1BajWs6uLiwqJFi3jvvfd47LHH6Nu3r3ZNiVu3bhEeHk55eTlLliyp0xYDgv5qY2ul/XPrjp7kJF/j9vWLXD2+l7WtbOnv7S4W3WkBaj3Tc/jw4Xh4ePDtt9/y559/auetm5qaEhQUxKuvvtqk3k4VasbWSsnEEYHazaA7+A4mP+MWJfk5JEQfZluoL68/O/g+dxGaujpNDXdxceGTTz4BNCtiqdXqZrfPh1DZmEG+FBSVsP3PKORGCjr1GU5M6GYyrp5j2697UZqZigldzdwDj31aWlqKZNGCBA/vx9D/LRKstHPC3rUXANej/mBryH9Ff0YzJyZLCLX2/JAA7Z+deg7AyExJSX4Ot86Fs35vOHE3772HjNB0iYQh1FpFfwaA3EhBR1/NqElq7Cnybycz+6stYivGZkokDKFO7nzfxKpdZ2w6dAVJ4saJECS1mpXbQsXU8WZIJAyhzu7sz3D2CUJupKAwO430K2eQgN+ORjdugEK9EwlDeCAV/RlGJma089I0UxLPhlNWUsiuI6dFLaOZaZCEUVZWRlJSUp2uVavVfPXVV/Tv359evXoxdepUEhISqi2/e/du3N3dK30lJiZqy+zfv5/hw4fj5eXF6NGjOX78eJ1iEyqztVIy4znN/Av7Lr0wtbJDVVrMrbPhopbRDDVIwrh48eI9NzO+l5UrV7Jp0yY++ugjtmzZglqtZsqUKZSWllZZPjY2Fn9/f8LDw3W+KmaaRkRE8PbbbzN27Fh27txJ3759mTZtGleuXKnz9yfoGhLQnSUzxyIzMKBDb83/9/TLZyjITBG1jGamQRKGo6MjM2bMqPV1paWlrFu3jpkzZzJw4EA8PDz44osvSElJISQkpMpr4uLicHd3x87OTudLLtesDPXtt98yePBgJkyYQOfOnZkzZw7dunVj/fr1D/Q9CroqtmJU2rfXdIACN06GopYkUctoRhokYTg4OPDGG2/U+rqYmBgKCgro27ev9pilpSWenp5ERUVVeU1sbGy1G0Or1WpOnz6tcz+AgICAau8n1N2o/t4AtO/1KAaGRhTcTibzRoyoZTQjdZoaDpCTk8NXX33F6dOnyc3NrXReJpMRGhpaq3umpGgm/Nz94pq9vb323N0xpKamcvLkSTZt2kRWVhZeXl68/fbbdOrUidzcXAoLC3F0dKzR/SokJyffs0nVkvYsqQ1bKyVPPerDriOnadM1gFvnwkk8G4Z1e1d+OxrNy6Puv1qboN/qnDA++OADDh48SP/+/avdCLm2ioqKAHQW4wFQKBTk5ORUKh8fHw9oVi7/5JNPKC4uZtWqVbz00kvs2bOH8vLyau9XUlJSLzELukb192bXkdM4ePiSdvkvSgtySY09xS65ISP7e4s3Wpu4OieM//73v8ydO7de170wMTEBNH0ZFX8GKCkpqbQ0IICvry/Hjx/H2tpauxHPihUrGDhwIDt27OC5557T3u9O1d2vQps2bUQtoo7urGU4efXnWuR+ki9GYNu5B9tCI8UbrU1cnfswzM3NcXJyqs9YtE2RtLQ0neNpaWk4ODhUeY2NjY02WYDmNXsnJydSU1OxsrLCzMysVvcTHlxFX0brTt0ws7JHVVZK0vnj/BFxXryc1sTVOWGMGzeOtWvXUlBQUG/BeHh4YGFhQWRkpPZYbm4uFy9exM/Pr1L5rVu3EhAQQGFhofZYfn4+169fp0uXLshkMnx8fDhx4oTOdZGRkfj6+tZb3IKuilqGTCajvfdAANLj/6I4N5P1e8NFB2gTVucmSXBwMDt37uTRRx+lU6dOlar4Mpms1kOXxsbGBAcHs2zZMmxsbGjXrh1Lly7F0dGRoUOHolKpyMzMRKlUYmJiwoABA1i2bBmzZ8/mn//8J8XFxXz++efY2Njw9NNPAzBp0iSmTZuGp6cnAwYM4JdffuHSpUssXLiwrt+6UAMVfRmWjh1o1daFnKSrJJwJw7X/aNE0acLqXMOYN28e165dw97eHhMTEyRJ0vm610LB9zJz5kyeffZZbf+IXC5n7dq1GBkZkZycTGBgIPv27QM0TZgffviBwsJCXnzxRV5++WWUSiU//vgjCoUC0CxcvGjRIjZv3syYMWOIiIjgm2++qXYoVqgfd77R2r7XoyCTkZ0YT15aAn9EiM2dmyqZJElSXS708fHh9ddfZ+rUqfUdU6OqGE4VnZ71Y+W2UEIiz3M9KoT0y2cwt3Gk69BgZDIZE0cEMmaQaBrWl4fxs1vnGoaxsTHdu3evz1iEZqji5bR23ftpJnNlppB5IwZA9Gc0QXVOGE899RSbN2+uc9NDaBkqmiZGpua08dQkj8SzYahVmjkyYtp401LnTk+lUsn27dsJCgrCy8sLc3NznfMymYxFixY9cIBC01exePDW8jLS4v83mSvuNG26+rPryGkxoasJqXPC2LFjB61atQLg/PnKy7HdOTdCEIKH9yMtM5esm/+bzHUhAluX7hgpzIi6cIVh/Xo1dohCDdQ5Yfz555/1GYfQAvh3d+HIaU9SY09RmJ1G0vnjdOj9GGt2HsbQ0JAhAaJPTN/VuQ/jgw8+4MyZM/UZi9DMeXRsi4GBge5krrwsJGDV9oOiA7QJqHPC2L17d73O8hSaP1srJdOfG0yr/03mkiQ1iX8dAUAtSSRnZDdugMJ91TlheHt760zhFoSaGBLQnU9njtVO5spKjCcvXbOc4pm4m40cnXA/de7DcHd3Z+3atfz+++94eHhgZmamc16MkgjVcXN25LWXRvNh7CnSr5wlIfowXYeMY/ufUZibKsRkLj1W54Rx4MAB7O3tKSsr49y5c5XOi1ES4V66tHegXY9+3L5xSbMy181YWnfwYP3ecLETvB4ToyRCo2hja4WRqQVtuvpz69wxEs+EYe3UBQO5oVidS4/VOWFUCAsL48SJE+Tm5mJjY0Pv3r3p379/fcQmNGMVM0DXlZeSdvkMpQU5pMScpG23PmIylx6rc8IoLS1l+vTphIeHI5fLsba2Jisri9WrV9OnTx9Wr15daWk8QbjTmEG+JGdkk514mavH95J84TitO3qiMLcUtQw9VedRkq+//ppTp06xZMkSzp49S3h4OGfOnOGTTz7hr7/+YtWqVfUZp9BMPT8kAJsOXbGwc0KtKich+hCAWGlcT9U5Yfz222+88cYbPPnkk9o9QAwNDRk9ejRvvPEGe/bsqbcghebL1krJ6IG9NRsgyWRkJcSRk3IdCdgWKobt9U2dE0ZmZiaenp5VnvP09CQ1NbXOQQkty6j+3phZ22PvqlkL9Oapg6hVKrHQjh6qc8Jwdnbm1KlTVZ6LioqqtLeIIFSnYg3Qdj36Yagwozg3k9Q4zc/W9j+jxMLBeqTOCWPs2LGsXr2a7777juTkZMrKykhOTubbb7/l22+/5ZlnnqnPOIVmblR/bwyNTXDqpenoTDr/X0oKNBtkiYV29EedR0lefPFFLl68yLJly/jss8+0xyVJYsyYMUybNq1eAhRahoph1h9+k8i4co78jFvcOBmK64AxyGQyMWqiJ+qcMAwMDFi4cCGTJk0iKiqKnJwcWrVqhb+/v1hgV6iTioV2fsq9zYXffyQn6Yp2BqiYm6EfHnjilp2dHYMGDdJZqi8pKQmAtm3b1vp+arWaFStWsG3bNvLy8vDz82PevHm0b9++yvLx8fEsXbqUM2fOYGBggJ+fH++88472s1UqFd7e3pW2RnzjjTf4xz/+Uev4hIZVsdBOZkIcSeeOcfPUQVo5dsBQYSoW2tEDdU4YN27cYM6cOfdcE+PSpUu1vu/KlSvZtGkTixcvxtHRkaVLlzJlyhT27NlTaSJYVlYWkyZNwsfHh59++onS0lIWL17MlClT2LlzJwqFguvXr1NSUsKuXbto3bq19tq7X5YT9Id/dxcOdw0g62YsRTkZXI86QOd+o8RCO3qgzgnjo48+4vr167zxxhs4OjpiYFDn/lOt0tJS1q1bx1tvvcXAgQMB+OKLL+jfvz8hISGMHDlSp3xoaCiFhYUsWbJEuxfr0qVLGThwIKdPn6Zv377ExsZiYWFRbxtGCw3Po2Nb5HI5nQKGcenARrISYrl9vTO2nbqxclso3u4dRNOkkdQ5YURFRbFw4cJKv8QPIiYmhoKCAvr27as9ZmlpiaenJ1FRUZU+q2/fvqxcuVJn4+aKxJWbq+lhj42NFX0qTUzFQjsrt4XStvsj3DoXzs1ToSjt26MwtxQ7pzWiOicMCwsL7SLA9SUlJQWg0hwOe3t77bk7OTk5VdoQes2aNZiYmGj3Yo2Li6O8vJzJkycTExODg4MDEydO5Kmnnqo2juTkZO2mMFURmxw1vCEB3enQxpa3v1STnXyVgowkrkXswz3oBf6IOI9jayuxbkYjeKB9STZu3EgdN06rUlFREUClvgqFQlGp07IqP/30Exs2bOCtt97CxsYG0HSKZmdnM378eNauXcvjjz/Ou+++y/bt2+stbqFhuDk7MnqQLy59hmNgaEReWgKpMZpJXGJuRuOocw3D1NSUU6dOMWTIEHr06KHTLIC6rbhVcY/S0lKd+5WUlFTa7PlOkiSxfPlyVq1axeuvv8748eO153777TdUKpV23xQPDw+SkpJYu3Ytzz77bJX3a9OmjahF6ImKTZ2dvQdxPSqExLNHUdq3x7y1o2iaNII6J4ydO3eiVCpRq9VVjpTUZcWtiqZIWloazs7O2uNpaWm4u7tXeU1ZWRnvvvsuv/32G++++y4vv/yyzvm7ExmAm5sbu3fvrnV8wsN354Su7ORrZCfGc/nYLro9PoE/Is6jNDMleHi/xg6zxdCrFbc8PDywsLAgMjJSmzByc3O5ePEiwcHBVV4ze/ZsDhw4wGeffcaIESN0zuXm5jJ48GDeeecdnn76ae3xc+fO4erqWu/xCw2jYt0MVVkJF7PTKcnP5mrEPlwHPM32P6MARNJ4SB544lZ9MjY2Jjg4mGXLlmFjY0O7du1YunQpjo6ODB06FJVKRWZmJkqlEhMTE3bs2MG+ffuYPXs2/v7+pKena++lVCqxtLSkT58+fPHFF7Ru3ZoOHToQEhLC7t27Wb16dSN+p0JtPT8kgJDI83QJfJKLBzaRk3SV5IsRtO3WVySNh+jBJ0/Us5kzZ/Lss88yd+5cXnzxReRyOWvXrsXIyIjk5GQCAwPZt28foOmfAFiyZAmBgYE6XxVlFi1axPDhw/n3v//NqFGj2LdvH1999ZVYRrCJqWiamFk70KG3pt/i1tlwsm9dAcRbrQ+LTKrPYY5moGI4VXR66qcN+46x/c8orkeFkH75DAaGRnQd/BJm1vYAfDd3coud1PUwfnb1roYhCPcSPLwfzwb54dz7MZQOzqjLy4gP20lZkWYXPrFKV8MSCUNocoKH92Ngb0+6BD6FidKa0sJc4o/uRF1exh8R54m7WXmSn1A/RMIQmqQJIwIxNDbB9dFnkBubUHA7mcvH9qBWq5j91RbRn9FARMIQmqSKTlATpTWuA8ZgIDckJ+kK1yL2I0kS6/eGi/VAG4BIGEKTNWaQLxNHBKK0c6Jz4FPIZAZk3rjEzVMHkSRJjJw0AJEwhCZtzCBflswci1VbFzr1GQ5AWnw0t86Fa2sa4p2T+iMShtDkuTk7MnFEIK07dqWD7xAAki9EcOvsUSRJ4ptfDoqkUU/0aqanINRVxXqg20G7g1ryxUjUqnIATl5ay8QRgeKV+AckEobQbAQP70duQREhgIFczo2ToaTGnkJdXkYH3yGs3xtOQVGJmEL+AESTRGhWnh8SAIC9qzcd/Z8AIP3KWeKP7kRVVsr2P6PE6MkDEAlDaFYqhlsB7Dr30IyeyA3JSbpKTOhmSgvzRNJ4ACJhCM1OxXArgE17NzyCXsBQYUZhdhoXQzaQn3GL7X9GsXJbqOgMrSWRMIRmacwgX76bOxmvLu2xsG2L59BxmFi2pqwon5jQLaTGnuKPiHNM+XitmKtRCyJhCM2WrZWSmWOHAqCwsMJz6Dis27sjSWpunv6Ty0d/payogPV7w/l8435R26gBkTCEZs3WSsmM5zTrZ8iNFHTuNwpnnyBkBgZk37rM+X3fc/tGDGHRsUz5eK1optyHWA/jLmI9jOYpIzuPbaGR/BFxHoDCrDSuRuyjKFuzSpuVUxfaew/CxMIKgAHe7kwYEdik1tZ4GD+7ImHcRSSM5u3OxKFWqUi+eJykCxEgScgM5Di496Zttz7IjRQATWqy18P42RUTt4QWxdZKyevPDkZpZsr2P6No1yMQ6/buJJw+RG7qDVIunSDj6nkcPXyxd+3F+r3hXLuVjn93Fzw6tm1SNY6GIBKG0CIFD++HuamC9XvDMbOyw23Qc+QkXeVm9CFK8rJIPBNG8sVI7N18OFiUT9hfsQAM6OVO105tUJqbtsgEIpokdxFNkpbl7r4NtVpF5vVLJF+KpDg3U1NIJsOqXWfsOveklWNHZHdsPD6gl7ve1D5EH0YjEAmjZcrIzuOLTfu5cDUJ0Oyml5UYT0pMFAUZSdpyhgozrJ26YO3khtLBGQO5XHtuaJ/udGpjq/37w66FtMiEoVarWbFiBdu2bSMvLw8/Pz/mzZtH+/btqyyflZXFxx9/TFhYGDKZjBEjRjB79mydrRX379/P119/TWJiIi4uLsyZM0dnh/g7iYTRssXdTGHbgUiiLl3THivKySD9ylluX7tIeWmR9riBoRFK+/ZY2jujdHDGzMpOp/ZRIaCbC73cNBtzKc1NsbdpRVpmDkC9JpQWmTBWrFjBhg0bWLx4MY6OjixdupTExET27NlTaZNmgPHjx1NUVMSHH35Ibm4u77//Pn5+fnz66acAREREMGXKFGbPnk2/fv3Yvn07GzZs4Ndff6Vz586V7icShgCaGsdvR6P59chp7TG1WkVeWgJZCfFkJ8ZTVlygc41MboiZlR1m1g6Y2zhg2qo1CqU1hsam99w61KuzE91d26M0Vegcr0guVxJSuJWejaWFaaUyd/r4nZkYGxpy6FD970pYQa8SRmlpKX369OGtt97ipZdeAjTbHfbv35+FCxcycuRInfLR0dGMHTuWffv2aX/5w8PDmTJlCkeOHMHBwYHJkyejVCr58ssvtdeNHTsWNzc3FixYUCkGkTCEO93dx1FBkiSKstPITb1JbmoC+ekJqMpKq7yH3EiBidIahYUVRqbmGJlaYGxqgZGpBUam5hgamyA3MtFp3tTF5YMbADgf3XBT3fVqlCQmJoaCggKd5oKlpSWenp5ERUVVShgnT57Ezs5Op6bg7++PTCbj1KlTPPHEE5w+fZp33nlH57qAgABCQkIa9psRmoWKYdjnBgcQcz2JqAtXORIdi0wmw8zaATNrBxw9/JAkiZL8bAoyUyjMTKUwK5WivCzKCvNQlZVQkJlCQea9tz8wkBsiN1YgNzJBbmSMgdwIA0NDDAzkyOSGGBgaYWAgx0BuiIGhITIDOTKZATIDA2QGclTlZcgNjRr0eehVwkhJ0TzQil3cK9jb22vP3Sk1NbVSWWNjY6ysrEhOTiY3N5fCwkIcHR1rdL8KycnJ2ppGVfSl9iFqQxoP4znYWikJ7OVOYC93xo8IJOZ6EnkFxQBcT84gJOIcJkprTJTWtO7QVXudqryMkvxsivOyKC3IoayogNKifMqKCigr1vxXVVYCaFYKUxeVazdlqi0zMzPkFq0e/Ju9B71KGEVFmg6lu/sqFAoFOTk5VZavql9DoVBQUlJCcXFxtfcrKSmpr7CFFqYiedzpucH+Oknk0KmLxN1MRW5opOnXsLKr9n6SWo2qvBRVaTHlpSWoyjRfalU5alU50v/+e+eXpCpHrVYjqVVIajWSpKY8J7VBv2/Qs4RhYmICaPoyKv4MUFJSojPqcWf50tLK7caSkhLMzMxQKBTa+919vqr7VWjTpk2L/1dbqJ27k8iwfj2Ju5nCqYvXMDQ0wMJU8/Os6ci05OCJ89p+EZmBAYbGJhgam1B9l+b9VfRhNCS9ShgVzYu0tDScnZ21x9PS0nB3d69U3tHRkdDQUJ1jpaWlZGdnY29vj5WVFWZmZqSlpemUSUtLw8HBoQG+A0H4m5uzI27OjtWee25wAFEXrnArPZtWFqbapFKhIrlUjJJUVeZOb/2+DklS1ev3cDe9ShgeHh5YWFgQGRmpTRi5ublcvHiR4ODgSuX9/PxYtmwZN27coEOHDgCcOHECgN69eyOTyfDx8eHEiRM899xz2usiIyPx9W0aLxQJzZetlZJh/Xrdt1x1Sedu/6cuf8CI7k+vEoaxsTHBwcEsW7YMGxsb2rVrx9KlS3F0dGTo0KGoVCoyMzNRKpWYmJjQs2dPfHx8mDVrFvPnz6ewsJB58+YxevRobQ1i0qRJTJs2DU9PTwYMGMAvv/zCpUuXWLhwYSN/t4LQ9OjdAjozZ87k2WefZe7cubz44ovI5XLWrl2LkZERycnJBAYGsm/fPgBkMhkrVqzAycmJiRMn8uabbzJgwADmz5+vvV9gYCCLFi1i8+bNjBkzhoiICL755psqJ20JgnBvejVxSx/06NEDlUpVabhWHyUnJwOVh6FbGvEcNBITEwGIjY1tsM/QqyaJPlAoFFWOvOijlv4LUkE8Bw2FQlHlNIP6JGoYgiDUmN71YQiCoL9EwhAEocZEwhAEocZEp2cLMWHCBLKysjD43wIv69ato3Xr1o0c1cO1bNkyDh06hIGBAa+//jrDhw9v7JAaVWpqKsHBwRw4cKDG14iE0QJIksStW7cIDQ2950IuzVlERATnz59nz5495ObmMnz4cAYPHtzgowr66vjx43z44YdkZGTU6jrRJGkBrl69ikqlYuLEiYwZM6ZFrgXSp08fvvvuOwwMDEhLS8PY2Bj5Ay5Y05T98ssvOotK1ZSoYbQAubm59OnTh/nz55OVlcW4cePo2rVrteukNleGhoZ88sknbNy4kddee61FJ4xly5bV6ToxD6MZ+e2331iyZInOsWHDhvHuu+/qHPv4449xd3fXeSGvuajJM8jOzmb8+PH8+9//btYvIdbkWXh7exMdHV3je4qE0QKcOnWKsrIy+vTpA2gSRrdu3RgzZkwjR/bwXLt2jdLSUu0yCZ9++int27fXrh3bUtU2YYg+jBYgPz+fzz77jNLSUjIzMzl8+DCPPPJIY4f1UN28eZMFCxZQXl5Ofn4+x44do1evXo0dVpMj+jBagEcffZTTp08zevRo1Go1//rXv1rcAkIVz+DJJ59ELpcTHByMp6dnY4fV9EiCXvrmm2+k4OBgnWMqlUpavny5FBgYKPXs2VOaMmWKdPPmzUaKsOGJZ/A3fXkWokmihzZu3FjlkNfKlSvZtGkTH330EVu2bEGtVjNlypQm83ZtbYhn8De9ehYNmo6EWklJSZFeffVVqVevXtITTzyh8y9KSUmJ5O3tLW3cuFF7LCcnR/Ly8pL27NnTGOE2CPEM/qaPz0LUMPTIhQsXMDIyYvfu3fTs2VPn3P02eWouxDP4mz4+C9HpqUeCgoIICgqq8lxtN3lqqsQz+Js+PgtRw2gi7rXJU0vZlEk8g7811rMQCaOJuHOTpzvdb1Om5kQ8g7811rMQCaOJuHOTpzu1pE2ZxDP4W2M9C5Ewmog7N3mqULHJk5+fXyNG9vCIZ/C3xnoWotOzibjfJk8tgXgGf2usZyESRhMyc+ZMysvLmTt3LsXFxfj5+Wk3eWopxDP4W2M8C/G2qiAINSb6MARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMARBqDGRMIR7Em8OCHcSCaMJcnd35+uvv27wzzl48CBz5szRORYUFMQ777zT4J/dUNatW8dbb71V4/LBwcHs27evASNqWsTbqkK1fvjhh0rHVqxYgYWFxcMPph5cuXKF1atXs3v37hpf89577zF58mQCAgJo3bp1A0bXNIiEIdRKY+4WtnnzZubPn1/teVNTU06dOlXtruxLly5l5MiRtVqRytPTEy8vL1atWsXcuXNrG3KzIxJGE1dcXMx//vMf/vjjD5KSkjA2NqZnz57Mnj2brl27AlBWVsby5cvZvXs32dnZ+Pn5MWrUKObMmcPBgwdxcnKqdN/x48dz4sQJQNME+vHHHwkICCAoKAh/f38WL14MaJooTz/9NLm5uezatYvS0lKCgoJYsGABGzduZMOGDRQUFPDII4+wYMECrK2tAdi2bRs//PADN27cwNbWlmeeeYbp06dX+8sO0KNHD+1u5Ddu3OA///kPr7zyCh4eHoBmmf3qro+Li+Pw4cNs3rxZe+z8+fMsXbqU8+fPo1ar6dmzJ2+++WalPVdHjRrF3LlzmT59OjY2NjX539J8NdiOJ0KDcXNzk7766itJkiTpH//4h9S3b19p27ZtUmRkpPTzzz9L/fr1k4YNGyap1WpJkiTpnXfekbp37y6tXr1aCgsLk2bPni11795dcnNzkxISEqr8jPj4eGn06NHS6NGjpejoaCkvL0+SJEkaNGiQNGfOHG25QYMGSd7e3tIbb7whHTt2TFq9erXk5uYmPf7449L48eOlw4cPSxs2bJC6du0qzZ8/X5IkzbZ/7u7u0kcffSQdPXpUWrNmjdSjRw/p3XffrfEz2LFjh+Tm5iZdvXq1RuWXLVsmBQYGap9JXl6eFBAQIP3zn/+Ujh07Jh06dEh6/vnnJR8fHyk3N1fn2oKCAql79+7Sli1bahxfcyVqGE1YaWkpBQUFzJ07l+HDhwPg7+9Pfn4+ixcvJiMjg6KiInbu3MmcOXOYNGkSAP379ycjI4Pw8PBq792lSxdtX8X9djm3sLDgiy++wNDQkEceeYSdO3eSmprKtm3bUCqVABw9epTTp0+Tl5fHypUreeGFF7RV/MDAQKysrJg7dy6TJk3C1dX1vt97XFwcCoUCZ2fn+5YFiIiIoEePHshkMgAuX75MVlYWEyZMwMfHBwAXFxe2bt1KQUGBNm4AMzMzOnfuzPHjx3nhhRdq9HnNlUgYTZixsTFr164FIDU1lWvXrnH9+nUOHToEaBJKZGQkkiTxxBNP6Fw7cuTIeyaM2vDy8sLQ8O8fJVtbW8zMzHR+6aysrIiLiyM6Opri4mKCgoIoLy/Xnq/YsOfYsWM1Shjx8fG4uLjcswlzp4SEBLy9vbV/d3V1xcbGhtdee40nnniC/v37069fP95+++0qr2/Xrh2JiYk1+qzmTCSMJu7o0aMsWrSIq1evYm5ujoeHB2ZmZoBmDkVmZiZApR7++uzxr2rUpCKGu2VnZwMwbdq0Ks/fvWx+deLj4/H3969ZgEB+fr7Ofh3m5uZs3LiRVatWsX//frZu3YqJiQlPPfUUc+fOrbRBkKmpKXl5eTX+vOZKJIwm7ObNm8yYMYPBgwezevVq2rdvj0wmY+PGjRw9ehRAOyKQkZFB27ZttddWJJKHzdLSEoBly5bRsWPHSudtbW3ve4+CggJSUlJwcXGp8edaWVlV+oV3cXFh6dKlqFQqzp49y65du9i8eTPOzs5MmTJFp2xubq62w7YlExO3mrDz589TUlLCtGnTcHZ21rbPK5KFJEn07t0buVzOgQMHdK4NCQm57/0NDOr/x6Nnz54YGRmRmppKjx49tF+GhoZ8/vnnNar23759G6BWIxbt2rUjOTlZ+/fff/+dPn36kJ6ejlwux9vbm/nz52NpaUlSUlKl61NSUmjXrl2NP6+5EjWMJqxbt24YGhqydOlSXnnlFUpLS9mxYweHDx8GoLCwEHd3d5555hk+//xzysrK8PDw4MCBA9p+jnslBUtLS6Kjozl+/Dienp60atXqgWO2trZmypQpLF++nPz8fAICAkhNTWX58uXIZDLtEOn97iGXy9mxYwdyuZxBgwbdt4nVr18/Nm3ahCRJyGQyfHx8UKvVzJgxg2nTpmFubs7+/fvJy8urtK9HXl4e8fHxvPLKKw/0vTcHoobRhHXo0IHPPvuM1NRUXn/9debNmwfATz/9hEwm4+TJkwB88MEHjB07lnXr1jF9+nRSUlJ4/fXXger7GgDGjRuHkZERU6dOJSwsrN7ifvPNN3nnnXc4cOAAU6dOZenSpfTu3ZsNGzbodJRWR6lU8sEHH5CSksL7779PcXHxfa8ZOnQoWVlZnD17FtDscv7dd9+hVCp5//33efXVV7lw4QJff/01ffr00bn26NGjGBkZMXDgwDp9v82J2JekmcvOziYsLIz+/fvrtME//fRTduzYobPVXnP32muvYW1tzSeffFKr6yZOnIibmxvvv/9+A0XWdIgaRjNnamrKwoULmTVrFocOHSIyMpLVq1ezYcMGxo8f39jhPVSzZs0iJCSkyj6K6pw7d46YmJhqR3VaGlHDaAEuXbrEl19+yV9//UVRURHOzs6MHTuWcePGaTtKW4o1a9YQExPD559/XqPyL730Ei+99BIjR45s4MiaBpEwBEGoMdEkEQShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxkTCEAShxv4fOX7PlePDsHYAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fcs.plotFCS(time_limits=[5*10**-7, 0.1], normalization='G0')" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "N 0.258003\n", "D 44.769663\n", "w_xy 0.313000\n", "w_z 1.745000\n", "tau_t 6.794145\n", "A_t 0.202325\n", "tau_t2 194.279258\n", "A_t2 0.344981\n", "y0 0.000000\n", "c_nM 0.450065\n", "Rh(298K)_nm 5.478213\n", "Name: 0, dtype: float64" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fcs.parameters" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Extract timetrace" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "trace = pp.pamplot.Timetrace.fromfiles('data/timetrace_GG.json', 'data/timetrace_GR.json')" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdUAAADVCAYAAAABzmpDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAABklklEQVR4nO2dd1gUxxvHv1c4ekeKYkEUEBUEEcXeYxJj7FFjjT2WaFTUmGj8WRM1mlijxtiNvRtj7xUbNqyIoDTp9eDu9vfHccfu3V7ljjtwPs+TJ7I7u/ve7Oy8M++87zsciqIoEAgEAoFAKDNcUwtAIBAIBEJlgShVAoFAIBAMBFGqBAKBQCAYCKJUCQQCgUAwEESpEggEAoFgIIhSJRAIBALBQBClSiAQCASCgSBKlUAgEAgEA8E3tQDlTVhYGIqKilClShVTi0IgEAgEE5KamgqBQICoqCiD3fOjU6pCoRBisdjUYhAIBALBxIhEIhg6qeBHp1Td3d0BAGfPnjWxJAQCgUAwJR06dDD4PcmaKoFAIBAIBoIoVQKBQCAQDARRqgQCgUAgGAiiVAkEAoFAMBBEqRIIBAKBYCCIUiUQCAQCwUAQpUogEAgEgoEgSpVAIBAIBANR4ZRqbGwsQkJCcODAAfmxp0+fYuDAgWjUqBHat2+PrVu3mlBCAoFAIHysVCilWlxcjKlTpyI/P19+LCMjA8OGDUONGjWwf/9+jBs3DkuXLsX+/ftNKCmBQCAQPkYqVJrClStXws7OjnFsz549sLCwwP/+9z/w+Xz4+voiLi4O69evR69evUwkKYFAIBA+RirMTPX27dvYvXs3Fi9ezDgeFRWF8PBw8Pml44NmzZrhzZs3+PDhQ3mLSSAQCISPmAoxU83OzkZkZCR+/PFHeHl5Mc4lJSXBz8+PcUyWND8xMRFubm5K90tMTFSbSJkk2ycQCASCPlSImerPP/+MkJAQfPHFF0rnCgsLIRAIGMcsLS0BSLd5IxAIBAKhvDD7meqhQ4cQFRWFo0ePsp63srJCUVER45hMmdrY2LBe4+XlRWajBAKBQDA4Zq9U9+/fj7S0NLRt25ZxfM6cOThx4gQ8PT2RkpLCOCf728PDo7zEJBAIBALB/JXq0qVLUVhYyDjWuXNnTJw4Ed26dcPhw4fxzz//QCwWg8fjAQBu3LgBHx8fuLq6mkJkAoFAIHykmP2aqoeHB2rWrMn4DwBcXV3h4eGBXr16ITc3F7NmzcLLly9x4MABbN68GaNHjzax5AQCgUD42DB7paoJV1dXbNy4EbGxsejRowdWrVqFyMhI9OjRw9SiEQgEAuEjw+zNv2w8e/aM8XdQUBB2795tImkIBAKBQJBS4WeqBAKBQCCYC0SpEggEAoFgIIhSJRAIBALBQBClSiAQCASCgdDLUSk+Ph5FRUXw9fVFTk4OVqxYgXfv3qFLly7o3r27gUUkEAgEAqFioPNM9eLFi/j000+xb98+AMDs2bPxzz//IDk5GTNnzsTevXsNLiSBQCAQCBUBnZXq2rVr0bJlS4wbNw7Z2dk4ffo0Ro0ahYMHD2LUqFHYunWrMeQkEAgEAsHs0VmpxsTEYMiQIbCzs8OlS5cgFovxySefAABatGiBuLg4gwtJIBAIBEJFQGelamlpCZFIBAC4cuUKXF1dERAQAAD48OEDHBwcDCshgUAgEAgVBJ0dlUJDQ7Fp0yZkZ2fjv//+k6cDfPToEVatWoXQ0FCDC0kgEAgEQkVA55nqDz/8gKSkJEyZMgXVqlXD2LFjAQCjR4+GUCjE1KlTDS4kgUAgEAgVAZ1nqtWrV8eJEyeQlpYGNzc3+fHVq1cjMDAQAoHAoAISCAQCoXwoFAlhxbc0tRgVGp1nqoMHD8br168ZChUAGjVqhNevX+OLL74wmHAEAoFAKB9uJ0SjzYYB2Bi1x9SiVGi0mqlGRUWBoigAwK1bt3D79m2kp6crlTt//jzi4+MNKyGBQCCYIUJREXY/PI4WNRrD17WGqcUpM79cWg8A2HB7N0aE9TWxNBUXrZTq3r17cfjwYXA4HHA4HMydO1epjEzpdu3a1bASEio0FEWBw+FU2ucRPl623T+EDbd3Y/WN7bg5dr+pxSkzPC7P1CJUCrRSqj/++CN69eoFiqIwZMgQzJ49G3Xq1GGU4XK5cHBwQN26dY0iKKHisfrGNhx/dgFbey+Fm62z0Z934+09/Hh6OWa1+xbtajcz+vMIHzdPUl6aWgSDcfrlVbzJSDC1GJUCrZSqvb09wsPDAQBbt25FYGAg7OzsjCoYoWKTWZCNrfcOAQC23z+ESS2GGf2Z3x2fDwCY8d+SSjFzIJg3lcke8uPp30wtQqVBZ+/f8PBw5OTk4PTp08jPz5ebfemQpPqEqf8ulv+bwyGbIREqH6RdE9jQWalevnwZEydORGFhIatC5XA4RKkS8DD5mfzfXLLGSaiEkHZNYENnpbps2TLUrl0bM2fOhIeHB7hc44/WMjMz8dtvv+HChQvIzc2Fv78/pkyZgrCwMADA9evXsWTJErx69QpeXl6YMGECPv/8c6PLRdAOLhnREyohHB0MwIUiIW68vY8m3kGwFVgbUSqCqdFZqb569Qpr1qyRK7Ty4Pvvv0dqaip+++03uLq6Ytu2bRg+fDgOHjwIiqIwevRoDBs2DEuWLMGFCxcQGRkJFxcXRERElJuMBNWQET3hY2fZ5b9wJOYsIqqHYEXXH00tDsGI6KxUq1atitzcXGPIwkpcXByuXr2KnTt3onHjxgCAn376CZcvX8bRo0eRlpYGf39/TJ48GQDg6+uLJ0+eYOPGjUSpmglkpkqojOjSro/EnAUAXI+/ZyxxCGaCzr3d6NGjsXr1aiQklI/7tbOzM9avX4+GDRvKj8niZbOzsxEVFaWkPJs1a4Y7d+6wrvkSyh+iVAmVEWKAIbCh80z16NGjSE5ORqdOneDi4gIrKyvGeQ6HgzNnzhhMQAcHB7Rp04Zx7L///kNcXBx++OEHHDx4EJ6enozz7u7uKCgoQEZGBlxcXJTumZiYiA4dOqh85tmzZw0jPAEASDIGQqWEeP8S2NBZqXp6eiopsfLk7t27mDlzJjp37oy2bduisLBQKYm/7O+ioiJTiEhQgEc6H5Ox68FRvE6Pxw9tx5LBjYHhVqpIVYKh0FmpLlq0yBhyaMWZM2cwdepUhIaGYunSpQCkm6YrKk/Z39bW7F52Xl5eH/VsNDk3DbYWVrCztC2X531sI3qKovA28z2qOXqCb+LUbyuubQYAdK7bEk28g0wqS0XkQ14GBHwLOFgqJ7upzIMUCSUhyzZ6UmFqbfv27ZgwYQLatWuHdevWwdJSuj2Rl5cXUlJSGGVTUlJgY2MDe3t7U4hq1qTlZ6DbtlHo9PfQcnvmxzZTPfHsAvr+M9GsstQUFAtNLUKFI1eYh8+3jkCnTUNYz+sSUlPRKBIXm1qECovOM9WAgACNI7SnT5/qLRAbO3fuxLx58zBo0CDMmjWL8fywsDDcunWLUf7GjRsIDQ0tlxjaisaj5BcApCPR8kKxvbxIe4Pd0ccxJLQndj04iuY1QtGyVvmFaBmbrfcOAgDOv75hYklKqcSTKqMRl/lefQED1Wl8ViI23z2AwY26o6ZzNcPctIyIxCKA7KuqFzor1XHjxil1knl5ebh79y7evn2LqVOnGkw4AIiNjcXChQvRqVMnjB49Gh8+fJCfs7KywqBBg9CjRw8sXboUPXr0wMWLF3Hy5Els3LjRoHJUFkwxulacqQ7bNx3FEhGOxpwDAOx//F+lytVLwTy8zune7x+bCd7QsJlDuQYy9E06Ph8JWUm4/vYuTgz5yyD3LCsiidjUIlRYdFaqEyZMUHkuMjISjx49Qq9evcokFJ3//vsPxcXFOH36NE6fPs0416NHDyxevBhr1qzBkiVLsGXLFnh7e2PJkiUkRtWMUByEFUtEJpLk40JEq2cyUS0bYokYXB5TiRpqTTUhKwkAkJafaZD7GQKiVPVHZ6Wqjh49emDSpEmYM2eOwe45ZswYjBkzRm2Z1q1bo3Xr1gZ7JsGwaDM7PvXiCjrXbVkO0hgfcwmPpneMlXn9z1jQLQ4iiRgWPAvGebpS1bSPLwccs7FgaIOIDHz1xqA2obdv30IkIi/DnKF/97J1VYqiUKyHY4KEkmh1nUQLLfPTmeU6P99cMZfOk65USarIsiFmmbnRQ2rEGnwUKpqnsJgiM1V90XmmumrVKqVjEokESUlJOHHiBNq1a2cQwQjGgT5jkZm0pv27GHcTn+DQwLWsoQOq+O7YPMRmJGBf/1WwslDt1ECh/JyizAFzUarFYtoAt4J16uYGm9KkK0qxRGzy8ClDIhITpaovBlGqAGBnZ4eOHTti5syZZRaKUD6IJRJY8IDLcVEAgAuvb6JbPdWZpuhQFIVbCdEAgPuJT9CsRoj8nKJnsViinVLVZELTFYqisO/Rv6jv4YdA9zoGu6/mB5ffo9RBN+FJtHwHhFLo7ViTOVTTzK6iDWkq05rqh7wM/Pv8Ir4IaA8nawejP09npRoTE2MMOUyOzFOyoplpyoK0oyidYeqSK1koLk24wecymxFjhgTtZ25iiRh8nuGW+U+9uIKlV6TelOXpXWyKmSrbgISuCIg5T3foikXEMiih17eiEtJ3gGjogaW+VKY11cknFuD5h1g8SIrB0k9nGP15eq+pZmdn49KlSzh+/DiuX79erjvXGBqKojD28GxMODq38ifh56heB1p4cS3eaorNKyGvKJ92T+Y5xXVWbWdJIgN3/NHJphkAlncbmn9+NfrumohChQQP9I5eW2sBoRQxo/40rKnSzs85+zv6756EQhHtfWipKLts/gZ33j3SQ1rDUplmqs8/xAIArr29Wy7P00uprl+/Hq1bt8aoUaMwZcoUDBs2DM2bN8fq1asNLV+5kJKXhnuJT3D73UOmsqgA5AjzcCzmHHKFedpdQDE9GhVZc3OHVrfJFZbWk1DETBOpGDKjyYlDnTxlgS6jKckV5uFozDnkaPuOSsgqzMGRp2eRX1ygttzRmHN4m/UeR2LO4sjTs8grkpYnM1UmeUUFOPL0LLIKc7QqzxiUsNSfhGaRoA9aTj6/hNiMBFyNuyM/psr7Ojk3jfF3ZmE2pp/8VSv5jImmb/FR8nNcf2ue29gVioQ4FnNOKUTJ1qJ8NofX2da2f/9+/Pbbb+jduze6desGNzc3pKam4vDhw1i1ahWqVq2KHj16GENWo0FfO5GYy4KYlvx0ejmux9/DxdhbWKKFaYPeObB9ONpmWsqlDT4YI3IARYrmXy1nbmyzgbKQW6SbEjMUir92ztnfcSXuDs6/voHfPvtB6/vMOfM7rsffQ9S7h/hfx0kayy8rMXXfiL+HhZ2nMszwlWnmoS+/Xl6Pk88v4fzrG1j++SyN5emDEtZvhaZIiyXKXvDaDKKG7JumdMwcHN00mX+HH5D6zhwauBZe9u7lIZLWrLmxA7sfHkctp2rY3f8P+XE7gU25PF9npbp582b079+fEYtau3ZtNG3aFFZWVti6dWuFU6r0UaY+oSV0KIpC5MlfYG1hpVVHWFZkmx5fenNbq/L0tSGxRKykyHgc7TwY84pLlWpBcSHjnGIdajtLMrSJMsdkM1Vmp3ilZMZCn7log+zd/vfisk5t6eyr6xi8dxo87Fzlxww9YDE2FEVh1ulliEl9jaWfzkBtlxplut+pF1dw8vklANqbAcUK34rSedoAVJZbmT6AlFkM1JFRkKV0rFgswuC9U9HBtzmGhPbUSlZDo24QRv+NiTmpZqdUZelB32S+Yxy3LSelqrP5Ny4uDh07dmQ916FDB7x+/brMQpU39OTRik42iuQXF+C/55fxPjsFtxOilcw3CdlJuPTmNv57cRnvs5ONIm+OMA9X3kRJ83PqiFjBpKX48fC0zJdMN60WKM1UmUpV65mqgU2UdGeq8qQ811RVWRaefXjNGGhpa4I3F/KLC3H21XW8y07G5TdRet0jOilG7iOgTxw0w3zOomQeJj2T/1tmoqcvfdCXknRxPRKKi/DsQ6zWSzFlhe23qZupVlSvcku+QHMhA6DzTNXDwwPv37M7syQkJMDOTvs4R3OBrgQ07c7w561/8E/0McYxumcp3RTaY8e3RvE6/eHUUtxKiMbIJl/pfC1dcYklbEpVu5kqw/yr4CCjPFPVdk3VsB6HpnI6M9RTna0ckFGYrbaMpkGgjIo2Uy2iDYj02fwhISsJIw9KTbz6foNM719m/SXnfsDbrNJ+UKZU6f1HJu3dmYNHryqOxZxXOqZupkpf3jG0c6EhoJvPTdEH6DxTbd++PX7//XdER0czjj948AArV65E+/btDSZceUH/gDV17HffP1Z7XlHBGANZfOiuB0d1vpY+shRJxEq/V5X59977xxh96Ee8TIsDwByFK66pKjoqbb13EHsf/atRNk3m3/ziAkw89j8cePyfxnsB2q1NpednodX6fmi6theepRrGymKoD1kxLR4b2i5XVDylWtqGisUi/HlrF2b+t1Tr3/EqPU7t+Vmnf0NhsRCTjs3H3ocnWMsoWnXoxGcmMv7OLyoskbX0fTxOeamVrKaGbTcldUr1WMlGGKrK/X5tC+aeXWkWkRT/vbgs/3d5yaOzUp0wYQKcnZ3x1VdfoWPHjujXr5/8/05OTpgyZYox5DQqRSL6TJVdqablZ+B1ejwkbKYSsQgPk55BJBEzZnDGRp9nKc9Umb/37vvHEEnEyCrMwYu0N/LjYw7Pxv3Ep5j672KlZ8sUbGGxEI+Sn7PO9pde3oinqa/Uyqypw9z36CRuxj/AL5fWqy0ngz7DUbz3u+wkvM9Owfrb/8jl/f7EQq3ua06oaq+KmKP5V0JJ8Cj5OetAtFjBerTpzj6ce30dtxOilcpqgq1dnXl5FQsvrsX1+HvyWGZFmCFJzHsobit5+100CouFjLYfmxGvs1yGRtbO1WHBEhuubnLx29VN8n8XKXj+iyVi7HxwBCeeX1Ba0wSApJxUvMtO0iS2wZhz9vdye5YMnc2/dnZ22LdvH/bv34/bt28jKysLDRs2xDfffIOePXvCysrKGHIaFframypz2mdbRgBgX+z+/foW7Hl4AoNDuqOOay2jyGgoGI5KlERppJmYk4I/b+3EgcenkFuUj219lsLPzUd+PiVXuvUeXTlmC6UxypH//YKb8Q/QsU4L1mcP3ReJWk6q94vU5KGqc4gMbWRaLBHJTdtF4mL03DEOANCqVhN5mQ/5GbrdX9Vjy9F7s0jLdWNznKkeizmPBRfWoGn1YPzRdTbjHP135dC8uPV5R4Ui9jrKKFBvWlfnKa8YIrPv0UkkZCVhaqsR8mNCUREKigthbWHFGlJj7AF4oUgob+dXR+1WmViFTalq214UrVT0vlQx1E4sEePL7dLNUS6M2AFrC+PoClNvHqFXnKpQKIS3tzeWL1+OTZs2YcqUKdKk7MUVc7d4+ssXsbjG08+zxbHuKTEfbb13SCle1NAmB7o5Qx+YHYWIdRCx9d4h+Qd/LY7dU5KhVAtzceLZBdyMfwBAOgtQBdvolU02NrR1opJBr/mpJxYhPV/qaUlXzukKsWzTT/6qcWSvCboDmabO6Vnqa0Se/AWzTv+GVde36fwsmYewRpnMQKmuubED62/9I/9736OTAICb8Q+U6pw+A6e/o6xC7ZLM0JcSVMWeK+7zeyn2NmadWib/hvfTlhkU3yNbW70Rf1/JHC/37lXo56OTYjD2MHMgwcZTNSZksUSMeedX48DjU6zn6Z7FOSrCy16lvcXZV9eVji+5vBFb7h4AILUa/Hz2DxyNOacUM63o+V9AszqIJCJQFIVlV/7CtnuHkE8rS5ftacpLTD/5q9aJZ2QIRUX48fRvmHPmd42+MED5WWt0VqqvXr3C559/jp9//ll+LD4+HosWLUKvXr1UOjGZK9nCXCTlpsr/LhKLkFuUzwgcTspJZbmSnQQFj1+hqAgpuWlKjY9OZkG2PCD9fXaK2nWy2WdWaJThfXYKCkVC5ArzkKYwspcohAloWkNm64wpikIM7WPPEuZg7rmVGuXShCYFRHei0mawQi9z+91DrLm5HQAgFJd++I9TXjCuuRB7EzP++xXZwlx8yMvAi7Q3ciVJURTisxLVOs6IJWKGc5Fi/SXlpDLawvijc3Ex9hbOvLyKbfcPyZVLWn4GowNjew8f8jKw+OKfqiuARoKWJjfZb5TVnVgiRkJWUpkHh5kF2dhy7wD+urMXT1JeIj6LuSa54MIahlmQ/g3Qv0VtZ6r0GVN0EntmLcW2P+3kYpx5dQ2b7uwDALn/AACIFN45fcmITmxGAuNv2WxYcfY08uAsvEp/q+4nAACG7p+uILNYvv/qzfgHOBZzDr9cKm0D9PdHb2eqZsXfHGCPbc8vLsCamztAURT+iT6Of59fxPzzq5Gam84oR7cCvM9OQY6wdNBTUFyIF2lvsOfhCay6sU1le/7mwExciL2JH0//prIe2Ih69xCnX17FyReXcPjJaaTmpastX9ZwSW3RWakuWbIEHh4e2LVrl/xYREQELl68CCcnJ/z6q+mzgWgLRVEYvHcafr+2RX6sWCJC1y0j8dmW4XKzZmKu9kr1QeJTxt+xmQn4Ytso9N89ibV8sbgYn2wehs5/D0XUu4fosWMsIsuQUeVlWhx67BiLyccXoPeuCfhsywhGEDpzTVWiccNwNqX7T/RxvKR1CNlazh40wZZflQ6f5kSlTbiM4m97kyGdJSuapRR59iEWnTYNwedbR2DgnimYVfKx73t8Er13jsfqG9tVXrvu1i7G34rK8MvtYxhrt9lCZt2l5qUjqzAHn20ZwegIFTuE+KxE9Nz5rdrfQefQk9N4kKg5beP2+4fRe+d4/BW1FwCw/Orf6LVzHP59flHrZ7FRICrt4Iftn47eO8fj2YdSx7Codw/Rc8c4eUwvfeaRXpAp/7e2mcPoM6ZZKjpr5rJP6fOSWL53xQGfqpnRD6eWMf421JKCjIUX1qDXznE49+o6o23J2sqOB0fQe+d4/H13PyNOW1W9KZpvFRGKirDqRqkFJV0hrlamuKOTYtBjx1iMO/Kz/Fx+cQEjVpc+OKIrY9kgVZ0Vi433OaXWjaVX/kLXrSORXLI8xYY2s1lDoLNSvXv3LiZMmAAPDw/GcVdXV4wZMwY3bih7kpkbReJiTD/5K36/thmJOUyz05QTC+UdwKs0qeLIKtAurRmPw0VKHjNu9eLrmwCkQdKAtBH/emkDbidE48Dj/7Di2mZ5WZm55drbu9j78ATmnl2JD3kZeJedhEUX12llHrkQK33e3feP5SYWWe7LzXcPMAYQ447+jNMvVJtqAWkKvGU0Rw4xJcGKa38zyuj6MXA57M3uj+tb1Jq36WEJMm9LALgYewvLr/6tVF7VyFSTUlVEVqcrrm4GIFU8qth67yDjb7ZOS50H+ahDs3CKpQ4UO4RjMed0/h27otV7i++OPi7vQDdE7cb1t/fkXtuydnoz/gFmnVrGiNFUBUVRWHtzB4YfmIk1N7SLudxRUrf0d/chr1QxqTJjKlIoUm0ZkhFN+w2KsdaKM3O6AisSF2PhxbVayfE4WWoJsbe0lR/TNV0lnePPLgAA/rqzl3Fc1o+tvL4VAPDnrV2MWWNybhp+vbRe57zCCy6sYfz9O62/AqRtedHFdfLwJfogIq+ogDEopyu8rMJsLLm8ATfoqQ4paVawy2oS2dx59wi/XPwT+cUFrMs03baNVuqDZcRnJcrfh0y+TC1TVuqCzo5KHA4HBQXsmUJEIlGFWFfNEebJO0p18Lk8CEVFyChUznrChp2lrZIJgt6pFhYLse3+Iex/fBL7H59kuUOp0pB5JLrZOuPf5xeRmpeuVUYetvyWGQVZSMvPwFqWYPIt9w6ovd+H/Az5mrGhsODyWWea0UkxiE6KQUSNENgLbJFTlMfY35W+/puc9wEuNo4AgMiTvyjdK7coX6mjLJaIUFBcqDH2kw3ZfpmyTiKrMAeOVvYQS8QoFotgZWHJmlNWVd7ezIJsOFixx3SzeaPSf3teEXuHoglZW8wV5sFGYI3MghzYCqxhyRfgVdpbhlcnAEw6Pl/+76zCHBSLizH33B9Iy8/E6/R47PjqNwhFRSodTp5/iMXmkoHio+TnWskoUz6MkBpax5yalw6xRIyC4kJYWVihUCQERVEMpQUwZ6raQC+fI8xXGgzRZ6rb7h3SOn/wg6SnEIlFsOSVJh6gz9DpjAnvr2TpAKTKyMPOjXmQohgz+Nfp8Upl6AOQtbd24k1GAvY//g83xuxDfnGhVvu/nnp5hfH309RXjL/vvH+EO+/ZFXVqXjpjcBKXUTr43nL3AG6/eyhfVwekloM9D09gz8MTOPj1Wnjau8kH4BRFIVuYi2+PSDP58Xl8vFEwtWvDNwdm4ObY/cgR5uHsq2vIL8qHk5W9zvdRh85KtUmTJli9ejXCw8Ph4uIiP56ZmYl169YhPDzcoAIaA22DyeOzEjHx2P8YC+zqYPvQ6OaSjMIsxCg0Sk1y0Wc+mtYMAHaPxnfZyajtUl3jteWFpvrvtGkIfJy9EZuRgD+7z0cjr3oAmGuhQ/dFYmOPhWjo6c96jw5/DVI6FpP6Cm03fq2XzN8emQMLLh+FkMrQ+e+hmNN+Aq6/vYdrb++if/AX2HB7t9J1+SpS1X2yeRi61G2t9fNlM9VsYS56bv9W6xkbU5ZCJGQl4at/vmPMIM5+s1WrHX0G7Z0qN+G9zojHpqh92Bi1Bxt6LGB9D/rMyGT5Wdly6QJS5Txg92Qk5qTCxcZJPkPb3PtX1KviKy+nzoeBDfqs7lbCAxx+eoZxXlZf19/ew/rb/0Bb7ic+xQ+nlzGcZOgmUjpWfPbBSbdtozGr7Vh0q1eaye5l+lssurhO/jebrwV9WYaugHY8OIKV17didrvxWv4K/VDMCLX21k75v2+/e6j22h47xqJbQAfMaidd4vjz9j/4u2StG0CZBvqv099iwO7vjealz6F09ECIjY1F3759IRKJ0KhRI7i4uCAjIwP379+HQCDArl274OPjo/lGBkYikWDVqlXYu3cvcnJy0KRJE8yePRvVqzOVSYcOHZCU+wF2QzVvWh1WrSGiNLx8XVj1xRyMPzrXYPdT5OroPfj10nqlDgEAPO3ckKRmvcGc6RfUFRMiBuP3a5uVPqYbY/ah2breJpJMM53qtMBpNd7Q2jIkpCc6+EZg/NG5SuuwZUWbzE2aGBzSAyPC+iItPwPHn11Ai5qNseD8asbau7YcGbQeow/9qLQ0o4kVn/8If7fa2Hx3P3Y/PK7TtW42zhrXP79q+DlOvrik1Sy1jmtNhqOTNszv9L1KZx1/Nx/83vUndNn8jU73rOh42bsj2CtAnrdZET6XDwkl0SvrFgBkb3qGqg4eOHv2bFnEZMqk6wU+Pj44duwYNm/ejLt37+L9+/ewt7dH3759MXToUHh6ehpMOF1Ys2YNdu7cicWLF8PT0xNLlizBiBEjcPToUQgEzJyP2r4AdYve+iBLnGAsEnNS8E5FvmFzVqg1HKsyUr4p8k/0MfA4XNbZx30FxzBzwxAKFZCa6Y88PWNwhQqgzAoVkFpU3O1csfzq3xBLxNgYtUfve3XbNkqv6yYdn48v63VkHVTS+arh50pKVxuHIjZFbSuwYQ3ZqeOiu1J1tnJQee5tZiJ2G3gZRlfa+ITjYuytcn1mYk6K2sGVs7UDsgtzTZbnmw294lQ9PDwwffp07N69G//99x/27duHGTNmmEyhFhUVYdOmTZg4cSLatm2LgIAALF++HElJSTh1iiWGi2Vyvu7L/ykdU3T793fzwdRWI+Bm4yw/FuTpj4Wdp2DZpzMZZdm2GdLkaacrka1GMoLNe+8cr/fMujx21GGjafVgOFs7aiy348ERHKWlR5Mx5vBPGq+d1XasxjKDQ3pge99lGssZi1+7RGoso075KbY/OpNbDNNLJl1ZenmjzkkmJkQMhrUKs6cifRt+prGMJoUKAP2Du2r1PG2gr5XS0SfXrzrFUCAqZJg/jcG6L+dhYecp2N2PPQvRmPAB8n9XNZOdaaz4liodH02FeUmjJzExMcjLy0NERIT8mIODAwIDA3H7trInmWJO1VY1wxBStT5WfTEHtYssMPBaAjo9SoWNUIRGcVnwypDOkEaH90efBp8ynCK+b/ENOvg2h++zBHSWSBtanwafYmHnqXDPFqLF83RUzShAp0ep+PRBCqpmFMC2UITWMWno8DgVzrlF4FAUIl6k49MHKfBNZl+Hcs8SovnzdHAkFEBRmFvoi8Aj59HNKwwjwvqqrBtZdiOrIjFaPkuDfUExGsVlYZlvL7jZOKPeuxwEvM+Bz8v3aJguncH7Jueh7dMPaBSXVfrsF+kY0EDaGXk7eqIJ3x0jEgSYkeKCzg9T4JpThF71P8EQ+4YIfZMJGwtr9PPvjDYxaWj79AMax2ZiwPUEuOVIBxahVevDy94d8zpOlid1CH2TiZof1GeZqZFWgC7RKej0KBUeWaoHKY1jpfcK9qyHrnXbYnC6PTyyhIioHsJavnfd9rC+cBlD3MLweVwROrvVByAdHHX0bc4oWzslT143lsVitHqWji7RKfDMlLYTr8xCNHuZAQ7L4E1QLEHL5+nwyixEh8epGHwzGbvazUQbn6YAAI5E2hY8MwtRJzkPDePZFalfUi4aJOTAWWCPebxgBAutMTTDEW7Z0jqxKhKj5fN02BcU4zP/tvLrAt/loEt0irxNh77JROeHKfjsQTKavcxAi1dZaCxxRFi1hirrFgC806Vt+pPoFHSJTsHgK/FoFJeFVs/S0PFxKrwyS60K9d3rymWuU9K+O/o2h5+bDzoVOGCZ1yfy390lOgUtn6WBK5HW3Wd+bRGRWIyFkvr40s4fYbGZCInLQoOEbLR6lgYLkQRNX2XALUeIZi8z4Jiv2VHSy95dbVuzEEnQ7skHfHPpLerkADU/5KP9kw8Yejkewy++hXd6AayFYrR6loZJDXuitq0HvohORbcHqRAUS2AvsMXwsD5Y1Hkqqti64NcukajpVBU2NCdC/8RcBL7LgQA8DM5wQLjIHtVuP8I4nj86JxSjRloBvrn0FnWTcsEXSfBJdAq63k8GXyyBtVD6LdsVlqyLl/Qf/W8nw6Hk99sXFKPlszRYFZUOcuwKRWj5PB3WQukxWZ17ZBWiY+0IeN5/ijChLWo5e8uvCYnLQq3UfExtNQK1Xaqjd61W6J0ArGsXCUueAByKQtOXGQh+m4VPHqagw+NU2BaK0I7nhRZvcuBY4mhoVyhCq2dpsBFKZQ54n4OOj1MR8kb6HX3XfAhmtRmDtnH5qJUjhp3ABp6ZhYh4UdLnlcATU2j1LA3tn3yAtVAMZ2tHOFuXzvA/K3JFp0el91XEqkiML+4l4ZOHKdI2ZoSNDnQ2/5ojSUnSYGgvLy/GcXd3d/k5OqJsIfhviyCqIR1l3jpwCR2WSUe4P0mks9OGyEHomyy4lDTS6X3rIbjEYaZZ9UaIzUiQKiX3Osh/G4tni39ABwBzT0eDy5dW68zbuShKY5ou2j5LA79xE4iipa7dYbFZ2BPuhe73pGbbZq8z8HN3P1AcDnycvbHrqxVotq43Rl2Ig71QDEdYYFSHkXi+bA6SAWTHPMLIrcdx/e09eSKDtj5NsfgT6ebHHA4H8ztOxh9DIhAcn43PolPAowDRzfnw+qYjBl6VOjC8vzIXAwHM6e6HURdL18FqVvND8/8uAADqNOdj4ph94HA4uDtxILIfSjModQDQxyEQQa1H4UK7QAQCaDRwCW7ePomQaObvD47PwfS+9bDyizngcbjgcDiwt7RFjbQCfHVLWvdtzj3G6ZdXWbfrGnf2jfzfQfHZWPZpqXPK/gGr0WvnOPgl5qLvbem92v70LxIObEf9M7dQH0DbGccw69QynHl1DVVsXeTOX7kH/kH8zo0ILLmXXYYV5v25HxRFSeuQ+h4cDgfiIiEufyJVzCsmbEL8vq14/0DqgNEuJg3T+9bDpFPSEKYCARcPajBn4Z8/SEaz15mMY2/HDkbt80+k7eFNlrwtyFjoao0sm9KBIFdCYfglaV7Z2kH98XrdEtzZ9Q/qAagHaVv98m4SQt9mIyw2E3aTbeDvVhtx719gSMn7rpOSh53NqsnrnMGd92h7fhNyi/IZDl8NPPywscdC3Et8gqRhX8GmiLmMUv99qWm60+MPmN5X+r1Maz0SY3ZOkcs898u6WNB5CqMud4+YhNh7K+TXz+j9E1xbtIMw6T1uzl8L4CUS9u1DHwVRuz5gtq8msVY40qcF4jLfo7Zzdbxmyb+b/TRa/rtz165C14D2mHJioTxGNuJlBro8kobA+f+rvLzgl5yHV1Vs4JuajypOZ/A/V2+8i7kAAOhdpy38J0qtXtUdvdDeVzrQb+PTFBdjbyLy5K/giSl8c1kqV83vhiBuzwLUB/AGt1ADQA1IvykA8E/Kw/kAV7SLkYaJ2ApFcJLwUTs+HQGJudjYtiY21uqDtL3zAADt3AIR9MufuDikCxCfAs8sIfaFVwUAfH0tAbU/FCA8xwK/NbZHxKsMfFnS1hr3XII7o6W+CW3OPYaNhTWcUtLR76Z0Wabt7E+l7/X0A2Tee4oPFr/j1KzNGL7gS/S8y+xjeRIKHZ5K+6KZn0/DSfs8cP83F9UyhQgU26Lb9N/waNAX8vIuPv4YENwNyWePw+Gm1GTe9vwTXGgn/RpFPC5u13YCAAQk5jDeeZNOrZEtzMPamzvgyLdGmz2l3sq5Vjy88GR62Ld/+gEtX0hN/cU8Li7bOCm937Kis6OSOXL48GFERkbi6dOnjETXkZGRSElJwebNm+XHOnTogNzEdxjiZYMdjasi+H02+r6Trjcu5XhgKsW+JmmzZSfCazQCULKn6t+L4HbyPAK++wmPf57EKOvdewjEBflIPL5X+UZaULBoIXI3rkF1K1c0/WMH7ic9RdbX0u6E5+oGu2o1kRVdGl5Ta9gEvPm7NKOR7/ZDqF7Nj3FPWQOlY9WuEwrPn2YcKwpqAEG0+lg2p9BmyLyrHI9s5VkNhUlSt3m/yXOQn56KhC1rlMpxbO0QvuEArL2kI+K4jHeI7dmp9LyFBdzadkHytXNAYQG4YmnnXehgB6ts5pri+0U/wyY9Cy5/bUXd4ZPw9OAWcF+WZntqc+4xHs4cg/Sb0thPt1Yd8eGydADlu/0QnuQlw8naAfxFi5DzlJmsvW2JoovdvAofLp2GTQ0fpF4sXU6o2n0A3h/aCW0RWwrAE6o28dmNGoe0+zdheUt5/1AJnwfqhx8Rby1BzYUrwMmROsu4tmiPtKtMs3j2wP5w2F4amsGzsYU4n2kBEVtZwmbGTxD+/COrLO4duyJw1q+4FHsb2cIcZBXmIiKDg7f/m85aXhUUjwdLZ1eIBRYQv5e2Deeps9CgQ0/cGNAZxRlSZeFQPwTZj0tjFmuPmoIa/YcjalRv5L54otMzPWbMRfJiaeiFReRM5N6PguWp05DwuPBYtAJu2UI8nS8ddFp6VIUwWao4kiInwe30RfDvaZf6URXND17BgynDkff6GQRu7mj2zxlkPYjCkwWRSOzeBbZ+9WE9Q7W53tBwHJ1AZWUyjtl06Yr8k6VbWDb6fSvufzcYANBk81G8unQC6ZtUx+JyeHy0ORONiztXgNrA3OCC37I1RFekjkU1B4+FW8euuDP4c7Uy2gc0hL1fIN4fkXrQuzZvh7RrpVvSUX+uR0xGHPwevAavJPFQdmgQui7dCYqisOfhv/A9ehb5Z5k7WL1ZOg95xQVofj8eWYf2AEKmdeunLAEEzq6mdVQyR2RJ/IuKihgJ/YVCIaytleM2+aAQKCnEme/W4ObALvLjv7QOZnSadIJtqsr/bWNhDedd+yAGlBQqACTs26J0TBea8qog+lUCCpGAh7PGwbN1J8iMGeK0DyhQGAbRFSoAvBrYHba/rIdLeEsAqlP6KSpUALCNewdNBjQ2hQpArlABIPHEfuTFsYcPUXm5iNv+J+zr1EOVdl3gKeIjln6+uBipp48qrU0oKlQAaPE6C8mnDiM/PQ3PlvykdI0oLwei3NLrZAoVAHD1Kjr3GoyUi/8hNovdUSU75iHiSgYGebHMlIa6KFQAahUqAOSuXw0bG1uwrUpyRWLY796P8KlzcSdHvfcpXaECUFKoAMArFKLKo+dQFemXcuYYAqYvQGuf0g0HrnzZXEVp1XDEYhR9YM4mczeswbN79+UKFQBDoQLStlTw7q3OChWAXKECQPGvi2BZ8m+uWIKM+T9BHNhIfl6mUAHA89cVOj+LjdiNK5D3WppYouhDCl7+sRDvj0jDcNw2bUf9eX9A/QaShkVRoQJgKFQAjH6vIP6NWoUKANbeNSEpLoLDrXtQNLTKFCoAcLhcZN/VnBMgJ+YhcmJK/UHoChUAat9+guqZ6ch88ACyr7l6MR9cDhfCtBS0S5EgRkGhAsCn6QJ4dOqFCzOVJxUAIClkD3krC5VCqcrMvikpKahRo4b8eEpKCvz92eMYAeBJyWhVhiqFCkg/cku38lmcz6F1JOk3LiL9BjNFXFG6Zk/e6Omj0GTLMdjWqM1orJooVqFc1NFoxVa8XLUQuS9LYx1znqmf7Sad2I8kSE1xbi3034M3/p+/QIlVO8cUZ6SjOJM9wwqHb4HEfw/g+bI5rOcpiQR3x+q+EXxZYFOAMgpTEiFRTLyiIbWjOhL2bVV7XpiSBOuq9JA0wxi1irMykHJWfchLYVICXq5VTtfJ4fNh5eWNgvg3ej1blJ2l9D0ZGsVvSKZQZTz+aaJRn68P7w6Upt4sytQcD8+3tUPiiQPIeiD1WbHzq49aQ77Fo1njFEpyGINtfYnduELpWGGK1IT/+OdJyH58n/W6pwunw7pG+YZ46qxUKYrC3r17cf78eRQUFDAStAPSNbwtW8o2U9OVgIAA2NnZ4ebNm3Klmp2djSdPnmDgwIEqr9NF2dyb8DWsqlaHa0Rb1B2vn+nGpmZt5Mdp3gg7doPyWqI2WDi5oJj2QeTHvoRtjdq4+20/ne9lVyeAoSRV4TdlLpyCw2BTw1e78lP/h+dLS3fnSD59BMmnj+gsnwxRjvpwkKK0VBS8Y4+VfLlS9f6pjg1DUZSmfc7n8qA4Iw33JjLbM6Wjt61X1z7IengX+QpWBK8v+oISiZD0b2mGrRd/LEBxVgYc6gUh7eYliLJVZxZrtucc4rauReIx/ZY8FEm/Vbo2VnPwWMRtlc6cXJu3R1H6B7VK1b39Z0g5Z7rwkw9XDGdKNAWqBpn29YLkSyTZTx4g+8kD+TmepSX4tiyZibhcgyhVNoTJiazLWor9rKYBvqHR2ft32bJlmD17Nl68eAGRSLq1D/0/RSVbHggEAgwcOBBLly7F2bNnERMTg8mTJ8PT0xOdO3fW+77W3jUZfxe+j8e7/dtAURS4AksVV6kmbP0B2NUJUHne45PuOt+TTtAv69Fkc2l+16RThyEq0C/PaMNF6yCgzcwdAoNRpW0XRplqvQahalfpWq/wA/tatCICZ1eV52x96motX42vR4FrpWzaB6RrspwSZzFFk622iIVC43yMCqnhmu9nD2rXFtlasTZ49x0K/ylzUWccc2cSrsAS/t//jIDI+Yzj6TcvISfmId4d3IHC96o33HYKaQqrKp7wmzwHVlUNn7mr5sDRcAwKg8C1CvynzgVXwB7GAgBNd/yHwJ+WwqE+u5e3ttjU9GU9zrU0zX7RHBV7ocqwrl5L/m+31tr1ed59h2r9fPeOXREwQ/UgVCwUws5X2SrI4XAg/KAcZ+rR+UutnmtZRfcwzcbr9zPaoS7fiCHQWakeOnQIw4YNw5kzZ7B9+3Zs27ZN6T9TMHHiRPTu3Rs//vgj+vfvDx6Ph7/++gsWFhaaL1ZB2MaDcGUxTYrzchmNWIZV1eqwrq7a1MCxsICFU6lS4dkw85U6BgbrLSsAcAUC2Nb0hVeJoku7dh6Jx/SLbbNwdkH41uNosuUYmu46jUYrtiDwxyVwblwatkQ3DVpX9Wa7jRI8FYoQANxad0LLo5rXXwCg5qAxELi4KR0XuLihxaFrcGokDVEpTNVuyzMZjkGNAQASYQHidqzXUFozrc9EI2I/3dxIwWfEJPlfAhc3+E+bV+bnaKL+zyvgO0a63GHlydwonr7mbumueydmXVVqHeJwuWjy1yF5HRoC/8j54FoIELzsLzTdfhIW9o5qB7TcEp+KRsv/hkP9Rno9065uIIJ/28R6jv7uykLNwWPhM1w7M3DTnacQse8C/BUGPXTqfFvqQBY461eErCxNEVilTWfY1KytdI3AyUXpmCpsa9aGbY3aaPwne38iERaCb2ePxuuULRVsyy/2foEI+IGZDMe2tp9Suapf6mZl41nbgCewZMQJK67PNvn7SOkAXsc9mrVB5zvm5uaibdu2BhekrPB4PEybNg3Xr1/HvXv3sH79enh7q+joaRXOs1ZO0iA/Z2kF/6nKSSGyntxH3iupI4JDg1D5cUmREJauVVTej8PhwN6vnvxvpVmbDi/Yo1M3pWOyzsbKq/R3v1pTmmzeu/cQ1P1Oc7IEnrUNuDw++Na2sK1RG9ae1cC1EIDD48G792B5Ob59aXyYzzffaSW3OqUqcHIF3449uXXgz0yTONdCAC7LgIlnYwu+ja38PvG7NmollyISobBMMWwWzq4IXfMPuDw+LF1obUIiQfU+Q1H1y/5ouFBq0mQ1mylQ42v9sgzZ+dVHjQEj4da6k7yjsVacTdIyjAX/przbDx36DKDupNnw+rwPfEZOkh/jWVmj3izmWqjnZ72k3p0B6mNfZTg0CIVL09bw6toHHh2ksdFcvoW87dCVqkODUIaFh2cl/Z65FgJYODhp9TxFBK5VlPoFno0tAmYsAk8HC5Vd3Xrw7NIDVbspKwY73wBYuldluUoZay9vCJxcGL/HKaQpo33a1PRFrW8mou53P4ErYP52K89qELgo90sWjs5Kx1Th3WeoVO669VBj4GjY1Q2ES7M28vOSYqkTnnW1GozrJMVFKM5U9tPg8PgQOJcOigVuHmi4YLVSOcWJhyZ4Gr4l9/afw7ZWHTSYtxKen/WCpavh/WR0VqqNGzfG3bt3DS5IuaKlUgWkozmvz5nRcQ+nj5b/23f0FPm/JUVCcCzYTVOWHtIPyD4gSH5M0QTC4XJRa6iGJNccDqyr+6DeD4vR5G/mWqSss7GwY093VmfcdFTr3h+OwaUenfV/XqFUjq+mM3Jt1gb2/g3AtbSCc0gz+XHLKh7w7NJDuXyL9qgz4QcAQLWeA1WabKXPlcZ0BsyUjmBlH1StoePg3uYT2PnVl5flcLkAlJWez3CpclelnDUh+w1iYSGoYv1Tn9Wb9Qsc6pW+6yrtpHF+Vbv1A1cggN+kn+AaIe2UODzNu4XUHjFJXj/a0vrUfYT9uRe1R05mjNw5PB48OpXGCVK04HqbajUZsxw6fHsHhPxR6tDiFBIO/6lzIVDonK3cvRiKLmDaPDReuxsB0xdoJXe1L/shaPE6+E9hN/XSjwX++CvsaDMcHs08K1sC0AR9cAhI60dxNhwwYxE8P/kSHB0sX44NQhEwfQH8Js9G6Dpm2kYLJ2et5PPuPaRUTtp3zbezh9+kUt8EroUAtQaNQbXu/aW/gVZHnJIBsSIWWs5Uq381TF6vHA4HtYd/h7D1+xC0aC08P+0JAKg1WJr4nmfLjAv9cOUsRLnKfg8cHg/2daXroXwHRzTfex5WntXg2LB0kuLcpKXG/lkRfkmfwbVk72dqDpFmV7OuVgMB0+Zp3UZ0kkHXC0aMGIFp06ZBJBIhODiYNWSlSZMmLFeaJ+o6eXkZgeoPif7xSYqKlGZPLY/ehCg/V/7h0jt7W5+68O47FAl7NgMAOFweag4agzebV7E+i2/vgLANB8C3l3auiqM4WWfDV6FUZQQv2YC82JcQuFVBYZJyzl1F86DS9b/9DXFhPnMGBsDSw0upLM/KGtV6fA2X8JawrlpDrdOC7IPw7NwNjg1DYVnFA4WJ7+Rr20qzXIWZZPi2E7CuJi0rqyNVhP11CC9XLULmvVJzc8DMxXJFKCkSykffgHTm6RzSDCnnmF6r1XoORPW+Q3GjX0fGccU1sHozFqF6nyGw81N2rNC0Xia/x6xfGQM6dTRYsBpcFQM8APCfOg/Jp0vW3xVyYTs2CGGNfw3beBAWNMWuzgzLZrGxrVUHTXeeQvyezerDkTRYCLh8WiIMSytGx0hXHhwttjYDpO09l+bwxuFywVV4J7LOnv5sTdC/dcXrBE6uKM5Q9rJ1adoatrV8Eb9bajHwHVsaoUC/H4dvwWg3ioMP+rvn8i1KBqFMtDX/qluf9p8yFzUGjJB/dxwOB7VHfY/X66UbA+S9Zt/yj8Pnw8LRCRF7LzC+64a//AlRViYoiQSWVTyQdu2CVjKW3ljadpzDIuRhTXQU36sx0PkJw4ZJ84iuXi2dqtNHwLLsM0+fmneSc7rMqsyRDeaXKja1azi0xsymVPl29oyPga4IuXwLWHvSTNRcLmvjlyHKy4OVR6nJSFmpSuXk2bHv0ykvZyGAfUnnLmKJYdOkVPklJlZFanz1DfLjXsM5pCmeL59b8iwLcDgc2HjXksqmZhBD/z2yxBA2tLVrnsIAjl5X/tPmyZ8BqJ+pBv2yHna1/WBby5ehVD06fSF3uJIICyEpkipVC2dXhPy+DUknSz1ja3w9CgUJb+Az/DvWuuAo5CPlCgSMmSsDnnYGI56K0Tcb6hzCZPLIYYljVjIRA7B0dQeHx0OtoeMhysuBtZp2UmPASBS8i4N7e2bQv7WXt8aZuWLdqYNraalytqHtLKTG16OQeHw/Mm5fUVlGpoB0manS27PieqFVVW/kx8cqXgK+vQNqDZuAwpQkuEW0ZbRxdUpacQDF6Ie4XIaDXM1BYyEuzGc4IqqiStsucGveTuV5Do/H+O4AwLvPELlSVX2d9N0ohinyrW3Bty6tNzbzb+3RU/D6T/Y83VSJo2ytoePwbv92ULR9eQHt20RZ0PkJW7eqj22rEHA4kMXcuTZvxzqaosdOqjLpAgoKVyIGR8NIlt5gOAILRocq+4Acg5vI478YKIRPKJpGZB+WLqMxtkZr7aVeqaq8l7UN6s/5DZKiIrlSVUSdZUDT+omyQi4dHHl91otxRp1SlSXFoHsW1ho6HhwOB3wb6YCEEolQXBJCErT4T9hUr8UwKXl0+Fytt7KNDrFxmpRMjYHS2akuHYK2s1+V17O0Y5mctYZ8q/F6vq0d69ICALhGtMG7/aodGu3rqV97lRSVZsXhCazg1KhkD2eFAam2dWDh4IQG8/7A5S6hGsvSlZlnl+5IOnlIZVl6e+ZwOAhauhHRU0fArk6AdPbI8j7d230KnqUV6s9WVhrMtXeKsR8oR81MFRQFr896If3GRdj7N4DPNxMAAMJ0zSFj9WYu1mp5gvFsvgV8v53O8OfgWlpBIizNB63tPRX7BEsPL9ToN1ylUpUIpXHcfGtbNFy8DtHTRjDO6zIo0hedv7yKsAm5JiiRCID0pToEBiNswwFEjewpP99sD9Pspc6MRleqfAcntWUBpiLk8gWMUbnMXNVg3h/SODweD0VpqSwB1SXXK3Qasoaqi9s/myKz8tBPqcrloDVcxVk+T41s9AEHG0rrK2qshJpM4ADg0KDUrCWTk29rB56tPcR5ORDn5ZSck75TulJXVDo8WzuI86S5XppsPqpxpkiHrfNvsukwuFbWKEpLlSsZXTo3dRYPbeAacUTvHBqh8lyDBavlVgpV0JOfcHg82NaojbCNB5XqXNUgpHq/4ZAIC/DuoNQEzeVbKNSt6oZFv6dTcLiSUpUpTkC533AOjUDo6l3yZAT0e3l06oZqPQaodebi0rLFSUTMGZhSX6CgZN1adkDo6l2MUCEuT1nB8Gxs4dgwVB6Goq8SUvztdIUKaD9A5FqW9h/W1X0Qulp9FrNiWiw1W9/G9psNjV5fXmxsLCZPnowWLVqgYcOGaN26Nb7//nu8esWels6c4VlZK8WOWinERqlVqhYChKzaCYf6IQj65U+NjZD+oilKwlz34Uo/Zgt7RzgEBsPBv4Fa0wsA1tAeh8BgVGmjXawamyOAJvOvJjgcDqr3Gw67OgHw+oK5g47Ug3iI1HtRC1nocK2Y59Vtr6WNo5J9nVJPbPrsR3GmLnv/TKXK7BQaLfsbDvVDEPLHdtiqiHFUBZuytPWpC2svbzg2CJF3mJqUKt2RCzrOLpRlMp5SVffeHLTwEGbLKGbn668UYqVqYMC1tIRzk5al8ggEWq+/0r9vSw8vePcdCsegMPkxa69Ss7liG+FwOHAIDJY7EtI7eAtHJzjUC1JbN/RZMiUSgaNG+TOUbMmynENgMPg0RyI2xSZwcUPtEZNhVycA1b8aptcWdoBmZaz1TFVQOpBwCgqDhQZfCXraQba2VB7mX52V6suXL9GnTx9cu3YNrVu3xtChQ9G8eXNcvXoVffr0qXCKVTZDkcWjsnnfqgs251pbw7F+I4Su2gGHgIYaHRnoioMqLlYw/+reEdLXfkvvw1VpelOETV6ramUP4PcdPQVhGw4wFJeMOuOm6+U+X+otW/JhqFl702amSn8X9I7a0p3pcCUb9XNpa7qK9WbvXx+hq3YwvBe1ResOXY2iq9ZzIOP3aNNpubWSOld5fd5b+Vnl0PmwoY2Vxb2jNMyGHs7Ghqo64HC4DKuIKkceVvn4NAcgCwHqjI1Eg/krGcfo91UrnwoHK21QXCvUB7Z3LBEKYVcnAGEbDsjjmvVBcSIie2fyZ2s5aKPPVBV9KtigR2pwuFzU+2kp87nmaP5dunQpvL29sW3bNtjbl84GcnJyMGTIECxfvhyrVrF7r5ojMg/PgMj5yHwQxT4zVPHBBf60TCluja6AQ1bvUryE+fFSFKND1cdkZ1ujNhqt2AqBq3IiBH1RnKkbA561jdLOJJo6VNfwVgj6ZT1sapXMBMs4U6VDzx/MVwgLkM9ULVXPVMuCth2qumdKiosYHZk2bSlg5iJkftIdzmHKifJNp1Q1x4HW6PcN7HzqwrGR+igDjgpTH4fHZToMKlqi1LQrxtJGyXWMQRHtWk2RBQxPZV2VqqhYBwXBnrOZ7R0batMy+oDC6/M+8PlmAlLOlCbx1/b30vsEdZMbQJrhrbZCcg5Fs3h5eP/q3Ivfvn0bY8aMYShUALC3t8eoUaNYNwU3ZxxL1tUsHJxQpVVH1pctWytTxD5QhTdnCTYsplkAqNbjazg2DIVTaFNm56enyc4pOEzJA08RxYB8OtW/+gYCFzdYeXmjjp55jfXB3r/UXOnZpbtWpiaX8JawKplJqjX/apFQAQDqTp4Nm5q+qDGwNLmCKgcweidmEqVKK+cc1gIWtDVEiVDI6HS0mQnwrW3h1qI96zq3PlaTsuLesatWIStcCwHcWnZQGY9digoFweGCZ0Ob1SsoJ3XNkC6fzILBt7VDlTafwLVFewhcq6BajwFwbBgKlyYt1ErHnKnqNoOixGJUafsJ7OsFofpX36gvq0JRsr9jwyhVDi0M0b1dF6VYYG2VG6Nt0uStPWoKbGv7Q0AL3arSqpOSMyPj2+JwdB686IPOPQOfz4elitGkQCBAUZH+AfPlhawx1Ro2QauPWKRi5xC2hqFpoRwA6k6cpSQLoF5JlAkuFx4K5hc6vmOmwnfMVOM8Ww2ODULku2METFedV1QVrhFtkPPsEWsQu+JHrIpq3fqhmkLGG56Cw5Tcq5oeH6lDvKImtDWF0cu5t+2C4CUb5AnFJUVCRhhFWR2VypJNSh+qdvsKfpPZE7nri6IzjwwOl8d4xzp5VTNmqqX/rk/L+FV3IvsetUr3oitoHTt7iUgEnsASjdf8o7mwKqVqxHfMGHxYWCi1ce3Nv7RkHjR5a/QfLt1vd2TP0o0vuMq/hzFwMeA3qw6dlWrDhg2xc+dOtG3bVilGdceOHWjQoIFBBTQGAjd3+EdOZGSVUYeqzDps5hdhamliea1GY2VcU9UGYzqdlIUqbbsgQFQMB3/t0tcpUqP/SFi6e7GaL3kCS9jW9lMZfK4OxcGQ7D2beqZKb0+KbU9SJGSGcBi6LRkhR6oM/2nz4N7+M4PflypWsTMwl8O0RuiwCQj9HWjy9NflXrqbf7VfUy3zAEsPmN+K8pq11tYZxnVsSpM+kFS+J73vK0uGNF3QuWf47rvv0L9/f3Tr1g1dunRBlSpVkJqaipMnTyI2NhZ//60+d6g5wOHx4PVpT80FS/DuOxQp507A1qcuMu5cL70PS8daa/BY5Dx7iNqjtZv5MRqClkkAtCXw5+V4/tvPCGSJeTMHOBwOPFlyGGsLVyBQ+x6dQ5sZRqmWDB4ZnaAhZ6paOyqVllOcKUuKhIxcrmU1cyl6luqzK5NWcHlKMcaGQiJiV6ocDhc8axs4hTaDpLCAZScU7UJqylonTPOvbl2xNo5K3n2GIOX8SVTtPkCHGxtoTVWDw5Zeg1K2mSj9m2QZPJSHuVcRvWaqGzduxLJly7Bq1Sp5FqUGDRpgw4YNFSpFobZYVfFExL6LyHn+mKFU2Waijg1D0eLwda1NK/SGoEsmGW1wb/MJqrTubFQzj1mj5whdVbwsPQONIetU+zVVWntTGIBJioTMjsrAnYm6+OIyQWk/S9T51qosTFweOBwOgpf+Jf1bz3dZ1pkq0zSp/RIAJRbBUYPnMyDducZ3bKRuv88IjkpsFj19rGds/SOj3tiUqgkc7vR6YrNmzbB3714UFBQgOzsbDg4O8hzASUlJ8PQ0vvdoecPhcJQDrFV43+nSiJnev4YfVX20ChX6D1LoTiyuNG9way9v1PvhF63Xa7VFH0clxd8mKSpijtoNrFS18crVCwN14mxI1Jh/AT2/DZq8ZQ3PYChVLQeATTYfwYcrZ1Hty/7aPUPDbwxe+hcKkt7h+dLZasvpCj35BOtMVZ/2yfJbmLme2WaqpXVMj002Jjr3OvXq1UN0tHT3d2tra3h4eMgValRUFD799FPDSmhGKDYEg6xV0mccKsy/slirgBmLyv68jwk9BxQ8WpIJl/BWjHMenb6AK23LK0OgtaMSfdRd8ttqDZVm26ozfqZeyeRVP4z+b45ejmSmhlJj/tUXS8+qsPTwgk2N2mU2/1rYlSYyYNvYgg0b71qo0W+4zru3qMK5cQSqssQplxWNM1UdZpDu7T8H19KKNZ6ayxiYsK2plh7T5I1tKLT6ZZs2bUJ+fj4AqUPS3r17cenSJaVy9+7dg0BDLFFFRnEdzRCzQKb3L/vH7tH+M7i17KDTXo4E6K1U6eEA5WE+Uhxhq4pvZM5US5TqkHGo3n8EeAJLpF46pfKeugtVWnetTt6tkG2PvsuQZRVPCGUb1mvyXVCX1YjHR9PtJ8HhcMv8/dNDoHTd1s9YGCxOVdOaqg6Dvno//gqquJg1TpXRJ2tYUy0vU7BWTxEKhfKEDhwOB3v3Ku/uzuVyYW9vj7FjxxpWQjPCGC+FsaaqpiOsiJ1aRYXRIZSD5zR9puravB18VTi5MWa09J2WStoGY4BWVvMvy/0rGnTzb8APi/Fg8lAAOsxUORxW87Qhw6ka/7kPKedOwLvH1wa7p1lAV2ZsM1UdBn0cDkcpl7H8nCZHJSOFwalDqx5j7NixcmUZEBCAPXv2IChIfeKDyogxPMmY3r/l76lGUIaRmagc3gn9Gb5jpqlMGsLoNDSNyk2QvEEXfIZ/h9i/ftc6plMfag4ajegHt+HR+UtYutKczLSsG78pc/F86WzU1GJXHn2x9wuUb8NYmWCsHrB9Qwby9eDoYP41q5kqnZiYGPm/CwoKkJubCycnJ1gYMadiYmIilixZgps3b6KoqAhBQUGYMWMG6tYt3Xrr33//xcqVK5GQkIDatWtj+vTpiIhQvRuGPhjlpdDcxA3t/UvQD4ZSLQ/zL73TYQkbYL1Gg9NGWeNK1SVrNwQ1B46G1+e9ddrNR1dcGjdH8wOXYeHkgsL38aUnNNVxSd1W/bw33CLaMjJXVXoMZf6lbz5hxIEpV4OzlymUql5fXlRUFPr27YvGjRujdevWCAoKwldffYUbN24YWj4UFRVh1KhRSE1Nxbp167Bz507Y2tpiyJAhSE9PBwDcuHED06ZNQ79+/XDw4EFERERg1KhRBk/ub4wkCmXN/UtQjb6DYboDSrkkzqB/+NoqM1alSutgyjoTKAevcWMqVPozOBwOs451GLwKXNw+Lg96AylVgZMLfEZMQu0xUzVu6VgWGN8nS//JSNZRTuZfnXvxu3fvYujQocjJycG3336LOXPmYOzYscjMzMSIESNw7949zTfRgaioKDx//hxLly5Fw4YNUbduXSxZsgT5+fk4d0667+mGDRvQsWNHDB48GL6+vpg+fTrq16+PLVu2GFQWY6yvGSL3L8Gw0NdvymWrKNrHTs9lqg77usomQ0POCJybSLNU8bTMoWzuML8z9m5PtpepV5ce5SFSpafm16NQQ0VeYsVt+vRFU1gSx0gJW9Shc4+xYsUKhIWF4a+//gKP9hGPHz8ew4cPx8qVK7Fp0yaDCVi3bl2sX78eHh4e8mPcksrLzs6GRCLB3bt3MWPGDMZ1TZs2xalTp2BI6C+wwbyVakrqAItHJ8G0lLv5l8tF84NXQEnESgnBFWlx9AbE+XmsnZIh11FtqtVE012nYeHgZLB7mhJtvOxDft+GwuT3Kte0CYbBrXVnje1cWzSnKawAa6oPHz7EsmXLGAoVkCq6gQMHYvr06QYTDgCqVKmCNm2YcYHbtm1DYWEhWrRogezsbOTn5yslnHB3d0dSUhLrPRMTE9GhQweVzzx79izrcfoL4mvYLFdbmN6/ZKZqDnDLYc9FRQQsmwKwYWHnoHJ3FkOvXVmXcbN6c0Kb74wrEHz0CtVQITXqMKTZn6spo5IG72BjoLNStbW1hUhFMmeRSKTzS0lISFCr4K5fvw4Xl9IO5/Tp01i2bBmGDh0Kf39/ueJUjI+1tLSEUCjUSRZNMEdFhnlBxsz9S9APRvo5sfHS6BkaU+Q5rTCQutES4ylVvoMjRNlZcG3e1mD31BxSwx6GZkx0VqqhoaFYv349WrVqJc+kBAD5+flYv349wsLCdLqfh4cHTpw4ofK8o2PpjHDXrl2YN28eunXrhsjISACQb0OnuOWcUChkyEfHy8tL5WxUHXTzrKYNtbW+pxFz/xJoHxGXB0jEWnnF0tdUKSPmpjU0xNKhGvp3RhlRcRBU03Tbv8iPfwOHwGCD3ZOjYZcvjj5OgGVEZ6U6ZcoU9OzZEx06dEDbtm3lu9RcuHABhYWFWLBggU73s7CwgK+vr8ZyS5YswcaNGzFs2DBMnz5druCcnJxgY2ODlJQURvmUlBTGOqyh8O49GMLUJNjVrWeYG2rIXUkwDCF/bMPLVYu02oSd4ZCmw7ZgpobMVFXD8BItBxNnhcWIdWPh4ATH+o0Me1O6ImXdxYa5SXl5oLNSrVmzJvbs2YOVK1fi4sWLyMrKgqOjI8LDwzF+/HjUqVPH4ELKFOr06dPxzTdMbzIOh4PQ0FDcunULffr0kR+/efOmzrNmbagzbobmQjpAvH/LB8f6jdB47W6dr6MqkFKVea8SlCEDVi2pYOMNTQlPVGUhMyZ6uUP5+vpixYoVBhaFnZs3b2Ljxo0YNGgQvvjiC6SmpsrP2djYwNbWFsOGDcOoUaMQGBiI1q1bY//+/Xj69KnOs2ZToI1XIsGEVCCPbHu/QAQt3Qgr96qmFsXsYHS4ZKZaaeDT9j5mXVNlLK+ZsVJNT0/HgwcPkJ2dzeqY1L1797LKJefYsWMApB6/27ZtY5wbP348JkyYgJYtW2LhwoVYs2YNli9fjjp16mDdunVamZVNjba5fwnlS7UeA5D9JBquTVubWhSdcGnc3NQimCd0J0CiVFVSHt6/hoRHTyyhqf80V6V66dIlfPfddygsLGR9ARwOx6BKdd68eZg3b57Gct27dzfoc8sN+osm5l+zwZg5aQnlD3HiqpzwNMxU6Vh6eBlbHAB6KNVly5ahZs2amD59Ory9veWJGAhlh8xUDUwFMt0SjAvd9Ee8f9VRseqGb2Mn/7eqgVPIH9tRlJkOm2o1y0cmXS+IjY3FypUrDZ6snkDWVAmEcqFi6Q2CGugzVVXmX8eGoeUkTYkYul7g5eWFgoICY8hCIDMrAoFgSirYmirdpGsuaV51nqmOGTMGf/zxBwICAlCrVi0jiPQRU8EaNIFQISHfWaXB2rMaAmYsAs/WeDvh6IpWSrV9+/aMUUBiYiI+/fRTODs7K2Ut4nA4OHPmjGGlrMQI3Nxh4eQCDo/HNGUQCARCeVMBBxyen3xpahEYaKVUw8PDzWZqXdng8viI2HMO4HCIoxKBUB5UQMVRXlS0kBpzRCulunjxYmPL8VHDSOBOMBxkIEhggXj/EoyJ3hvMXblyBbdv30ZmZibc3NwQERFhlLSABAKBQCBUFHRWqllZWRg9ejTu378PPp8PJycnZGZmYs2aNWjdujVWrlyptA0bgUAgmAsWBtoLmUBgQ+dFvIULFyI2NharVq3Cw4cPceXKFURHR+P333/H/fv3sXz5cmPISSDoDPEDINDxn74AVb/sD9eItqYWxXwha6plRmeleuHCBUydOhUdO3aUd1pcLhedO3fG5MmTcfToUYMLSSDog2uL9gAAC2dXE0tCMAe8uvSA36SfiEMgwajobP6lKApubm6s57y8vJCfn19moQgEQ+Dg3wDhW45D4OZualEIhIoBmamWGZ2HbD169MDatWuRl5fHOC4SibB9+3b06NHDYMIRCGXFpoYPY3soAoGgGhJSU3Z0nqlaW1vjzZs36NChAzp06AAPDw9kZGTg4sWLSEpKgqOjI2bOnAlAuqa1cOFCgwtNIBAIBII5orNSPXLkCOzspDsDXL9+nXHO09MTd+/elf9NHEUIBAKB8DGhs1I9d+6cMeQgEAgEgskh5t+yQtzgCAQCgSCFrKmWGaJUCQQCgUAwEESpEggEAoFgICqcUo2KikK9evVw8+ZNxvHr16+jZ8+eCA4ORpcuXXD8+HETSUggEAgVC9vafgAA1+btTSxJxUfvhPqmICcnB5GRkZBIJIzjr169wujRozFs2DAsWbIEFy5cQGRkJFxcXBAREWEiaQkEAqFiELRkA1IvnoJnp26mFqXCU6GU6s8//4zq1avj3bt3jONbtmyBv78/Jk+eDADw9fXFkydPsHHjRqJUCQQCQQOWLlXg3eNrU4tRKagw5t/Dhw/j3r17+OGHH5TORUVFKSnPZs2a4c6dOyRDCIFAIBDKjQoxU01ISMCCBQuwZs0a2Noqp5xLSkqCp6cn45i7uzsKCgqQkZEBFxcXxrnExER06NBB5fPOnj1rGMEJBAKB8FFhcqWakJCgVsFdvXoV06ZNw1dffYWwsDAkJCQolSksLFTaw1X2d1FRkWEFJhAIBAJBBSZXqh4eHjhx4oTK8//88w8KCgowYcIElWUsLS2VlKfsb2tra6XyXl5eZDZKIBAIBINjcqVqYWEBX19flecPHDiAlJQUNG3aFEDpLgojR45E9+7d8b///Q9eXl5ISUlhXJeSkgIbGxvY29sbT3gCgUAgEGiYXKlqYtu2bRCJRPK/k5OTMWjQIMyfPx8tWrQAAISFheHWrVuM627cuIHQ0FBwyYbEBAKBQCgnzF6pVqtWjfE3j8cDIDUbu7q6AgAGDRqEHj16YOnSpejRowcuXryIkydPYuPGjeUuL4FAIBA+XirFNK5u3bpYs2YNLl68iO7du2Pv3r1YsmQJiVElEAgEQrli9jNVRby9vfHs2TOl461bt0br1q1NIBGBQCAQCFIqxUyVQCAQCARzoMLNVMtKYmKiqUWokMhiiUkokvaQOtMPUm+6Q+pMP4yhD8hMlUAgEAgEA0GUKoFAIBAIBoIoVQKBQCAQDARRqgQCgUAgGAiiVAkEAoFAMBBEqRIIBAKBYCA41Ee2i7e/vz8AaRIJgvbIXM+9vLxMLEnFgdSZfpB60x1SZ/oh20qULaGQvnx0car29vZkj1U9IB+r7pA60w9Sb7pD6kw/LC0tlfbiLisf3UyVQCAQCARjQdZUCQQCgUAwEESpEggEAoFgIIhSJRAIBALBQFQ6pSqRSPDHH3+gVatWaNSoEUaOHIn4+HiV5TMyMjBlyhQ0adIE4eHhmDt3LgoKCspRYvNA13p78eIFRo0ahaZNmyIiIgITJ07E+/fvy1Fi06NrndE5cuQI/P395d6HHxO61ltxcTGWLVsmLz9w4EA8ffq0HCU2PbrWWVpaGqZMmYJmzZqhadOmmDx5MpKTk8tRYvPjzz//xKBBg9SWMYg+oCoZK1eupJo2bUqdP3+eevr0KfXNN99QnTt3poRCIWv5gQMHUr169aIePXpEXbt2jWrXrh0VGRlZzlKbHl3qLT09nWrRogU1YcIE6tmzZ9TDhw+pr7/+mvr000+pwsJCE0hvGnRtazISEhKoxo0bU35+flR8fHw5SWs+6FpvP/zwA9W8eXPq0qVL1MuXL6kJEyZQLVq0oLKzs8tZctOhT7/Wr18/6smTJ9Tjx4+pvn37Ur169Spnqc2H7du3UwEBAdTAgQPVljOEPqhUSlUoFFIhISHUjh075MeysrKooKAg6ujRo0rl7969S/n5+VEvX76UH7t8+TLl7+9PJSUllYvM5oCu9bZnzx4qJCSEKigokB97//495efnR127dq1cZDY1utaZDLFYTPXv358aPHjwR6lUda23t2/fUv7+/tT58+cZ5du1a0famoo6y8rKovz8/KizZ8/Kj505c4by8/OjMjIyykNksyEpKYkaPXo01ahRI6pLly5qlaqh9EGlMv/GxMQgLy8PERER8mMODg4IDAzE7du3lcpHRUWhSpUq8PX1lR8LDw8Hh8PBnTt3ykVmc0DXeouIiMCaNWtgZWUlP8blSptSdna28QU2A3StMxnr1q1DcXExRo8eXR5imh261tvVq1dhb2+P1q1bM8qfO3eOcY/KjK51ZmVlBVtbWxw6dAi5ubnIzc3F4cOH4ePjAwcHh/IU3eQ8fvwYFhYWOHLkCIKDg9WWNZQ+qFTJH5KSkgAoB0K7u7vLz9FJTk5WKisQCODk5PRRbWaua715e3srZaRav349rKys0KRJE+MJakboWmcAEB0djU2bNmHfvn0f7fqWrvUWGxuL6tWr49SpU1i/fj2Sk5MRGBiIGTNmMDq/yoyudSYQCLB48WLMnj0bYWFh4HA4cHd3x/bt2+WD34+F9u3bo3379lqVNZQ+qFQ1LFtQVsyQYWlpCaFQyFqeLZuGqvKVFV3rTZFt27Zh+/btmDp1KlxcXIwio7mha53l5+dj6tSpmDp1KmrVqlUeIpolutZbbm4u4uLisGbNGnz//fdYu3Yt+Hw+BgwYgLS0tHKR2dToWmcUReHp06cICQnBjh07sGXLFlStWhXffvstcnNzy0Xmioih9EGlUqoyc6RiGkKhUAhra2vW8mwpC4VCIWxsbIwjpBmia73JoCgKK1aswPz58zF27FiNnnWVCV3rbP78+fDx8UG/fv3KRT5zRdd64/P5yM3NxfLly9GyZUsEBQVh+fLlAICDBw8aX2AzQNc6+/fff7F9+3YsWbIEjRs3Rnh4ONatW4d3795h37595SJzRcRQ+qBSKVXZ1D0lJYVxPCUlBR4eHkrlPT09lcoWFRUhMzMT7u7uxhPUzNC13gBpmMO0adOwbt06zJw5E5MmTTK2mGaFrnW2f/9+XLt2DSEhIQgJCcHIkSMBAF27dsW6deuML7CZoM83yufzGaZeKysrVK9e/aMJR9K1zqKiouDj4wM7Ozv5MUdHR/j4+CAuLs64wlZgDKUPKpVSDQgIgJ2dHW7evCk/lp2djSdPnrCu9TVp0gRJSUmMhnbr1i0AQOPGjY0vsJmga70BQGRkJE6ePIlly5Zh6NCh5SSp+aBrnZ06dQrHjh3DoUOHcOjQIcyfPx+AdC36Y5q96vONikQiPHz4UH6ssLAQ8fHxqFmzZrnIbGp0rTNPT0/ExcUxTJb5+flISEj4qJceNGEofVCpHJUEAgEGDhyIpUuXwsXFBdWqVcOSJUvg6emJzp07QywWIz09Hfb29rCyskJwcDBCQ0MxefJk/Pzzz8jPz8fs2bPRvXt3lTO0yoiu9XbgwAGcOHECkZGRCA8PR2pqqvxesjKVHV3rTFEByBxMqlatCicnJxP8AtOga72FhYWhefPmmD59Ov73v//ByckJf/zxB3g8Hr788ktT/5xyQdc66969O/766y9MmjQJ3333HQBgxYoVsLS0RM+ePU38a8wHo+mDMoQAmSUikYj69ddfqWbNmlGNGjWiRo4cKY8FjI+Pp/z8/Kj9+/fLy3/48IGaMGEC1ahRI6pp06bUnDlzPqoEBjJ0qbdhw4ZRfn5+rP/R67ayo2tbo3Pjxo2PMk6VonSvt5ycHGrOnDlU06ZNqeDgYGrYsGHUixcvTCW+SdC1zl6+fEmNHj2aCg8Pp5o1a0aNHz/+o2xrdKZPn86IUzWWPiBbvxEIBAKBYCAq1ZoqgUAgEAimhChVAoFAIBAMBFGqBAKBQCAYCKJUCQQCgUAwEESpEggEAoFgIIhSJRAIBALBQBClSiCYGSTKjUCouBClSiCYEWfPnsX06dPlf9+8eRP+/v6MFHWmIDo6Gp988glrwnE2Bg4ciBMnThhZKgLB/CDJHwgEM0K208+2bdsASLc+e/nyJerUqcNIkF6eCIVCdO/eHVOmTEHHjh21uubJkycYPnw4jh07BldXVyNLSCCYD2SmSiCYMXZ2dmjUqJHJFCoA7Ny5E3w+X2uFCgCBgYEICgrC2rVrjSgZgWB+EKVKIJgJgwYNwq1bt3Dr1i25yVfR/Lty5Up06dIFp0+fRteuXdGwYUN8+eWXuHfvHu7fv48+ffogKCgIXbt2xfXr1xn3f/78OUaPHo3Q0FCEhoZi3LhxiI+PVytTUVER/v77b3Tt2pVx/NixY+jWrRuCgoLQrFkzTJ06FcnJyYwyX3zxBfbt24f09HQD1A6BUDEgSpVAMBPmzJmDwMBABAYGYvfu3ahfvz5ruaSkJCxevBhjxozB77//juzsbEycOBHff/89+vTpg9WrV4OiKEyePBmFhYUAgNjYWPTr1w9paWn45ZdfsGDBAsTHx6N///5IS0tTKdPNmzeRnJyMzp07y4/duXMHkZGR6Ny5MzZs2ICZM2fixo0bmDJlCuPa9u3bQywW4/Tp0waoHQKhYlCptn4jECoy9HXTRo0aqSxXUFCAOXPmoHXr1gCAly9fYtmyZViwYAF69+4NQLp/5sSJExEbG4t69eph1apVsLa2xubNm+XPiIiIQMeOHbFx40aGcxSdGzduwMHBAT4+PvJjd+7cgZWVFUaNGgWBQAAAcHJywsOHD0FRFDgcDgDAxsYGvr6+uH79Or766quyVQ6BUEEgM1UCoQISGhoq/7ebmxsAIDg4WH5MtkdrdnY2AKlyDA8Ph5WVFUQiEUQiEezs7BAWFoZr166pfE58fDyqVavGONakSRMUFBSga9euWLZsGaKiotCyZUuMHz9erlBlVKtWDQkJCWX6rQRCRYLMVAmECgib45K1tbXK8pmZmThx4gRrmIuLi4vK63Jzc5XuGxISgvXr12Pz5s34+++/sX79eri5uWHMmDFy72W6TDk5OZp+DoFQaSBKlUD4CLC3t0fz5s0xbNgwpXN8vupuwNnZGSkpKUrHW7VqhVatWqGgoAA3btzA1q1bMX/+fAQHByMoKEheLjs7G87Ozob5EQRCBYCYfwkEM4LLNc4nGR4ejpcvX6JevXpo2LAhGjZsiAYNGmDz5s1qHYmqVq2KpKQkRpanX375Bb169QJFUbC2tka7du3ka7Lv379nXJ+UlKRkPiYQKjNEqRIIZoSDgwNiY2Nx/fp1ZGVlGey+3377Ld6+fYvRo0fjzJkzuHz5MiZMmIDjx48jICBA5XUtWrRATk4Onj9/Lj/WrFkzPH78GDNmzMDVq1dx4cIFzJ8/H05OTmjWrJm8XE5ODl68eIFWrVoZ7HcQCOYOUaoEghnx9ddfw8LCAiNHjsSlS5cMdt+AgADs2LEDHA4HkZGRmDhxIlJTU7F69WpGuIwiYWFhcHV1xcWLF+XH2rRpg6VLl+LFixcYP348vv/+e1hbW2Pr1q1yBykAuHz5MiwsLNC2bVuD/Q4CwdwhaQoJBIJaNm3ahF27duHUqVNK3r3qGDJkCPz8/DBr1iwjSkcgmBdkpkogENQyYMAASCQSnDx5UutrHj58iJiYGIwaNcqIkhEI5geZqRIIBI3cvXsXM2bMwLFjx+QJH9QxYMAADBgwQCm9IYFQ2SFKlUAgEAgEA0HMvwQCgUAgGAiiVAkEAoFAMBBEqRIIBAKBYCCIUiUQCAQCwUAQpUogEAgEgoEgSpVAIBAIBANBlCqBQCAQCAaCKFUCgUAgEAwEUaoEAoFAIBiI/wPpN8aFCrGqHwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "trace.plotTrace(time_limits=[0,1])" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "fret", "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.9.15" } }, "nbformat": 4, "nbformat_minor": 2 }