{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Statistics in physt" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Necessary import evil\n", "import physt\n", "import numpy as np\n", "\n", "np.random.seed(42) # Repeatable results" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "data = np.random.normal(150, 25, 1000000)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": "Statistics(sum=149960006.08864358, sum2=23113238494.210587, min=29.26409978722107, max=266.9737275094181, weight=1000000.0)" }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAABVIUlEQVR4nO3deXxOd/7//+cVZCGuBNkriKWW2oM02mo7UkFqUJ3RVlu7oUmRTG2thrYz1eooOpQuQ5hSywxatGnT2L4IKpWxpyjS4qKtJiFIIjm/P/xyPq7GkigSp4/77XZuN+e8X+dcr3Od0uftXGexGYZhCAAAAHc8l7JuAAAAADcHwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIuoWNYNlGeFhYU6fvy4qlatKpvNVtbtAACA3yHDMHTmzBkFBQXJxeXa5+QIdtdw/PhxBQcHl3UbAAAA+v7771WzZs1r1hDsrqFq1aqSLn2Rdru9jLsBAAC/R9nZ2QoODjZzybUQ7K6h6OdXu91OsAMAAGWqJJeFcfMEAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsIEmy2WxasWJFWbcBAAB+g1IFu1mzZql58+ay2+2y2+0KDw/X559/bo4/9NBDstlsTtPQoUOdtpGRkaGoqChVrlxZfn5+GjVqlC5evOhUs27dOrVu3Vpubm6qX7++EhISivUyc+ZM1alTR+7u7goLC9O2bducxi9cuKDo6GjVqFFDnp6e6tWrl06ePFma3S0Tv/7+fj1NnDjxquseOXJENptNaWlpt6XXv//972rfvr0qV64sb2/vK9aUp+N9s3oBAKDcMkrh008/NVavXm18++23Rnp6uvHiiy8alSpVMnbv3m0YhmE8+OCDxuDBg40TJ06YU1ZWlrn+xYsXjaZNmxoRERHGjh07jM8++8zw8fExxo0bZ9Z89913RuXKlY24uDhj7969xj//+U+jQoUKRmJiolmzaNEiw9XV1ZgzZ46xZ88eY/DgwYa3t7dx8uRJs2bo0KFGcHCwkZycbGzfvt249957jfbt25dmd42srCxDktM+3GqXf3fTpk0z7Ha707IzZ85cdd3Dhw8bkowdO3aU+nMlGcuXLy/VOvHx8cbbb79txMXFGV5eXsXGy9Pxvlm9AABwu5Umj5Qq2F1JtWrVjA8//NAwjEvBbsSIEVet/eyzzwwXFxfD4XCYy2bNmmXY7XYjNzfXMAzDGD16tHHPPfc4rde7d28jMjLSnG/Xrp0RHR1tzhcUFBhBQUHGpEmTDMMwjMzMTKNSpUrG0qVLzZp9+/YZkoyUlJQS71tZBLvLzZ071ykwFRQUGK+88opx1113Ga6urkaLFi2Mzz//3ByX5DQ9+OCDhmEYxrZt24yIiAijRo0aht1uNzp06GCkpqY6fdaNBLur9VmkPB3vm9ULAAC3W2nyyA1fY1dQUKBFixYpJydH4eHh5vIFCxbIx8dHTZs21bhx43Tu3DlzLCUlRc2aNZO/v7+5LDIyUtnZ2dqzZ49ZExER4fRZkZGRSklJkSTl5eUpNTXVqcbFxUURERFmTWpqqvLz851qGjVqpFq1apk1V5Kbm6vs7GynqTyZPn26pkyZon/84x/auXOnIiMj9cc//lEHDhyQJPPnya+++konTpzQsmXLJElnzpxR3759tXHjRm3ZskUNGjRQ165ddebMmat+1kMPPaR+/fr9pn7L0/G+Gb0AAFDeVSztCrt27VJ4eLguXLggT09PLV++XE2aNJEkPfXUU6pdu7aCgoK0c+dOjRkzRunp6WbAcDgcTv9jlWTOOxyOa9ZkZ2fr/Pnz+uWXX1RQUHDFmv3795vbcHV1LXbdl7+/v/k5VzJp0iS98sorpfxGbp9//OMfGjNmjJ544glJ0ptvvqm1a9dq2rRpmjlzpnx9fSVJNWrUUEBAgLneH/7wB6ftvP/++/L29tb69ev16KOPXvGzatWqpcDAwN/Ub3k63jejFw8Pj1LtPwAAt1upg13Dhg2VlpamrKws/ec//1Hfvn21fv16NWnSREOGDDHrmjVrpsDAQHXs2FGHDh1SvXr1bmrjt8K4ceMUFxdnzmdnZys4OLgMO/o/2dnZOn78uO677z6n5ffdd5/+97//XXPdkydPavz48Vq3bp1OnTqlgoICnTt3ThkZGVddZ/78+TelbwAAcPuU+qdYV1dX1a9fX6GhoZo0aZJatGih6dOnX7E2LCxMknTw4EFJUkBAQLE7FYvmi84wXa3GbrfLw8NDPj4+qlChwhVrLt9GXl6eMjMzr1pzJW5ubuYdv0WTFfTt21dpaWmaPn26Nm/erLS0NNWoUUN5eXm39HPL0/G+Gb0AAFDe/ebn2BUWFio3N/eKY0WP3Sj6SS88PFy7du3SqVOnzJqkpCTZ7Xbz59zw8HAlJyc7bScpKcm8js/V1VWhoaFONYWFhUpOTjZrQkNDValSJaea9PR0ZWRkOF0PeCex2+0KCgrSpk2bnJZv2rTJ/O5cXV0lXbr+8dc1w4cPV9euXXXPPffIzc1NP/300y3vuTwd75vRCwAA5V5p7soYO3assX79euPw4cPGzp07jbFjxxo2m8348ssvjYMHDxqvvvqqsX37duPw4cPGJ598YtStW9fo0KGDuX7RIyc6depkpKWlGYmJiYavr+8VHzkxatQoY9++fcbMmTOv+PgLNzc3IyEhwdi7d68xZMgQw9vb2+mOx6FDhxq1atUy1qxZY2zfvt0IDw83wsPDS7O75e6u2KlTpxp2u91YtGiRsX//fmPMmDFGpUqVjG+//dYwDMPIz883PDw8jL/97W+Gw+EwMjMzDcMwjFatWhmPPPKIsXfvXmPLli3GAw88YHh4eBhTp041t61f3RX7zDPPGGPHjr1mf0ePHjV27NhhvPLKK4anp6exY8cOY8eOHeYjWcrT8b5ZvQAAcLvdssedDBgwwKhdu7bh6upq+Pr6Gh07djS+/PJLwzAMIyMjw+jQoYNRvXp1w83Nzahfv74xatSoYk0cOXLE6NKli+Hh4WH4+PgYf/3rX438/HynmrVr1xotW7Y0XF1djbp16xpz584t1ss///lPo1atWoarq6vRrl07Y8uWLU7j58+fN5577jmjWrVqRuXKlY2ePXsaJ06cKM3ulrtgV1BQYEycONG46667jEqVKhV73IlhGMYHH3xgBAcHGy4uLubjTr755hujTZs2hru7u9GgQQNj6dKlRu3ata8Z7B588EGjb9++1+yvb9++xR6xIslYu3atWVOejvfN6gUAgNupNHnEZhiGUUYnC8u97OxseXl5KSsryzLX2wEAgDtLafII74oFAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIso9SvFAAClU2fs6rJu4ZY58kZUWbcA4DKcsSuH+vXrJ5vNpqFDhxYbi46Ols1mU79+/W5/YyWwbNkyderUSTVq1JDNZjPfPnIlhmGoS5custlsWrFihdNYcnKy2rdvr6pVqyogIEBjxozRxYsXr/nZhw4dUs+ePeXr6yu73a4///nPxV4RJkmrV69WWFiYPDw8VK1aNfXo0eMG9hQAgPKHYFdOBQcHa9GiRTp//ry57MKFC1q4cKFq1apVhp1dW05Oju6//369+eab162dNm2abDZbseX/+9//1LVrV3Xu3Fk7duzQ4sWL9emnn2rs2LHX/NxOnTrJZrNpzZo12rRpk/Ly8tStWzcVFhaadf/973/1zDPPqH///vrf//6nTZs26amnnrqxnQUAoJwh2JVTrVu3VnBwsJYtW2YuW7ZsmWrVqqVWrVo51RYWFmrSpEkKCQmRh4eHWrRoof/85z/meEFBgQYOHGiON2zYUNOnT3faRr9+/dSjRw/94x//UGBgoGrUqKHo6Gjl5+eXqu9nnnlG8fHxioiIuGZdWlqapkyZojlz5hQbW7x4sZo3b674+HjVr19fDz74oCZPnqyZM2fqzJkzV9zepk2bdOTIESUkJKhZs2Zq1qyZ5s2bp+3bt2vNmjWSpIsXL2rEiBF66623NHToUN19991q0qSJ/vznP5dqHwEAKK8IduXYgAEDNHfuXHN+zpw56t+/f7G6SZMmaf78+Zo9e7b27Nmj2NhYPf3001q/fr2kS8GvZs2aWrp0qfbu3av4+Hi9+OKLWrJkidN21q5dq0OHDmnt2rWaN2+eEhISlJCQYI5PnDhRderU+c37de7cOT311FOaOXOmAgICio3n5ubK3d3daZmHh4cuXLig1NTUK24zNzdXNptNbm5u5jJ3d3e5uLho48aNkqRvvvlGx44dk4uLi1q1aqXAwEB16dJFu3fv/s37BABAeUCwK8eefvppbdy4UUePHtXRo0e1adMmPf300041ubm5ev311zVnzhxFRkaqbt266tevn55++mm99957kqRKlSrplVdeUZs2bRQSEqI+ffqof//+xYJdtWrVNGPGDDVq1EiPPvqooqKilJycbI77+PioXr16v3m/YmNj1b59e3Xv3v2K45GRkdq8ebM+/vhjFRQU6NixY3r11VclSSdOnLjiOvfee6+qVKmiMWPG6Ny5c8rJydELL7yggoICc53vvvtO0qWAOn78eK1atUrVqlXTQw89pNOnT//m/QIAoKwR7MoxX19fRUVFKSEhQXPnzlVUVJR8fHycag4ePKhz587pkUcekaenpznNnz9fhw4dMutmzpyp0NBQ+fr6ytPTU++//74yMjKctnXPPfeoQoUK5nxgYKBOnTplzsfExDgFvRvx6aefas2aNZo2bdpVazp16mT+XOrm5qa7775bXbt2lSS5uFz5P1lfX18tXbpUK1eulKenp7y8vJSZmanWrVub6xRda/fSSy+pV69eCg0N1dy5c2Wz2bR06dLftF8AAJQHPO6knBswYIBiYmIkXQpnv3b27FlJl+70vOuuu5zGin6WXLRokV544QVNmTJF4eHhqlq1qt566y1t3brVqb5SpUpO8zabzenGg5thzZo1OnTokLy9vZ2W9+rVSw888IDWrVsnSYqLi1NsbKxOnDihatWq6ciRIxo3bpzq1q171W136tRJhw4d0k8//aSKFSvK29tbAQEB5jqBgYGSpCZNmpjruLm5qW7dusVCLgAAdyKCXTnXuXNn5eXlyWazKTIysth4kyZN5ObmpoyMDD344INX3MamTZvUvn17Pffcc+ayy8/m3U5jx47VoEGDnJY1a9ZMU6dOVbdu3ZyW22w2BQUFSZI+/vhjBQcHq3Xr1tf9jKKzmmvWrNGpU6f0xz/+UZIUGhoqNzc3paen6/7775ck5efn68iRI6pdu/Zv3jcAAMoawa6cq1Chgvbt22f++deqVq2qF154QbGxsSosLNT999+vrKwsbdq0SXa7XX379lWDBg00f/58ffHFFwoJCdG///1vff311woJCSlVLzNmzNDy5cuv+XPs6dOnlZGRoePHj0uS0tPTJUkBAQFO06/VqlXLqZ+33npLnTt3louLi5YtW6Y33nhDS5YsMb+DY8eOqWPHjpo/f77atWsnSZo7d64aN24sX19fpaSkaMSIEYqNjVXDhg0lSXa7XUOHDtWECRMUHBys2rVr66233pIk/elPfyrVdwEAQHlEsLsD2O32a46/9tpr8vX11aRJk/Tdd9/J29tbrVu31osvvihJ+stf/qIdO3aod+/estlsevLJJ/Xcc8/p888/L1UfP/3003XP9H366adOd+4+8cQTkqQJEyZo4sSJJf6szz//XH//+9+Vm5urFi1a6JNPPlGXLl3M8fz8fKWnp+vcuXPmsvT0dI0bN06nT59WnTp19NJLLyk2NtZpu2+99ZYqVqyoZ555RufPn1dYWJjWrFmjatWqlbg3AADKK5thGEZZN1FeZWdny8vLS1lZWdcNVwBwNbxSDMBvUZo8wl2xAAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYGdxEydOVMuWLUu1zrp162Sz2ZSZmXlLegIAALcGwa4c+vHHHzVs2DDVqlVLbm5uCggIUGRkpDZt2mTW2Gw2rVixouyavMnWrVun7t27KzAwUFWqVFHLli21YMECp5qEhATZbDanyd3dvUTbbt26tdzc3FS/fn0lJCQUq5k5c6bq1Kkjd3d3hYWFadu2bTdr1wAAuG14jl051KtXL+Xl5WnevHmqW7euTp48qeTkZP38889l3dots3nzZjVv3lxjxoyRv7+/Vq1apWeffVZeXl569NFHzTq73W4+9Fi6FHCv5fDhw4qKitLQoUO1YMECJScna9CgQQoMDDTf5LF48WLFxcVp9uzZCgsL07Rp0xQZGan09HT5+fndmh0GAOAW4IxdOZOZman/9//+n9588009/PDDql27ttq1a6dx48aZr8aqU6eOJKlnz56y2WzmvCS98cYb8vf3V9WqVTVw4EBduHDhup/52Wef6e6775aHh4cefvhhHTlypFjNxo0b9cADD8jDw0PBwcEaPny4cnJyJEkvvviiwsLCiq3TokULvfrqqyXa7xdffFGvvfaa2rdvr3r16mnEiBHq3Lmzli1b5lRns9mc3mDh7+9/ze3Onj1bISEhmjJliho3bqyYmBg9/vjjmjp1qlnz9ttva/Dgwerfv7+aNGmi2bNnq3LlypozZ06JegcAoLwg2JUznp6e8vT01IoVK5Sbm3vFmq+//lrSpVdonThxwpxfsmSJJk6cqNdff13bt29XYGCg3n333Wt+3vfff6/HHntM3bp1U1pamgYNGqSxY8c61Rw6dEidO3dWr169tHPnTi1evFgbN25UTEyMJKlPnz7atm2b01sp9uzZo507d+qpp56SdCmQXekn0GvJyspS9erVnZadPXtWtWvXVnBwsLp37649e/ZccxspKSmKiIhwWhYZGamUlBRJUl5enlJTU51qXFxcFBERYdYAAHCnINiVMxUrVlRCQoLmzZsnb29v3XfffXrxxRe1c+dOs8bX11eS5O3trYCAAHN+2rRpGjhwoAYOHKiGDRvqb3/7m5o0aXLNz5s1a5bq1aunKVOmqGHDhurTp4/69evnVDNp0iT16dNHI0eOVIMGDdS+fXu98847mj9/vi5cuKB77rlHLVq00MKFC811FixYoLCwMNWvX1+S1LBhQ3l5eZX4e1iyZIm+/vprp9eTNWzYUHPmzNEnn3yijz76SIWFhWrfvr1++OGHq27H4XAUO6vn7++v7OxsnT9/Xj/99JMKCgquWONwOErcLwAA5QHBrhzq1auXjh8/rk8//VSdO3c2L/6/3hmvffv2FftJNDw8/Dev87///U8JCQnm2URPT09FRkaqsLBQhw8flnTprF1RsDMMQx9//LH69OljbmP//v3q2bPnNXspsnbtWvXv318ffPCB7rnnHqe+nn32WbVs2VIPPvigli1bJl9fX7333nsl2i4AAFbHzRPllLu7ux555BE98sgjevnllzVo0CBNmDCh2Nm02+Hs2bP6y1/+ouHDhxcbq1WrliTpySef1JgxY/TNN9/o/Pnz+v7779W7d+9Sf9b69evVrVs3TZ06Vc8+++w1aytVqqRWrVrp4MGDV60JCAjQyZMnnZadPHlSdrtdHh4eqlChgipUqHDFmoCAgFL3jxtn5fepAsDtwhm7O0STJk3MmxWkS6GmoKDAqaZx48baunWr07ItW7Zcc7uNGzcu9miPX6/TunVr7d27V/Xr1y82ubq6SpJq1qypBx98UAsWLNCCBQv0yCOPlPqO0nXr1ikqKkpvvvmmhgwZct36goIC7dq1S4GBgVetCQ8PV3JystOypKQk86ykq6urQkNDnWoKCwuVnJx83bOdAACUNwS7cubnn3/WH/7wB3300UfauXOnDh8+rKVLl2ry5Mnq3r27WVenTh0lJyfL4XDol19+kSSNGDFCc+bM0dy5c/Xtt99qwoQJxW4uWL58uRo1amTODx06VAcOHNCoUaOUnp6uhQsXFvvJd8yYMdq8ebNiYmKUlpamAwcO6JNPPjFvnijSp08fLVq0SEuXLnX6GVaSGjVqpOXLl191v9euXauoqCgNHz5cvXr1ksPhkMPh0OnTp82aV199VV9++aW+++47ffPNN3r66ad19OhRDRo0yKwZN26c05m+oUOH6rvvvtPo0aO1f/9+vfvuu1qyZIliY2PNmri4OH3wwQeaN2+e9u3bp2HDhiknJ8fp+j4AAO4EBLtyxtPTU2FhYZo6dao6dOigpk2b6uWXX9bgwYM1Y8YMs27KlClKSkpScHCwWrVqJUnq3bu3Xn75ZY0ePVqhoaE6evSohg0b5rT9rKwsp+fA1apVS//973+1YsUKtWjRQrNnz9brr7/utE7z5s21fv16ffvtt3rggQfUqlUrxcfHKygoyKnu8ccf188//6xz586pR48eTmPp6enKysq66n7PmzdP586d06RJkxQYGGhOjz32mFnzyy+/aPDgwWrcuLG6du2q7Oxsbd682ekGkRMnTigjI8OcDwkJ0erVq5WUlKQWLVpoypQp+vDDD81n2BV9b//4xz8UHx+vli1bKi0tTYmJidd9lAoAAOWNzTAMo6ybKK+ys7Pl5eWlrKws2e32sm4HsDSusbszHXkjqqxbACyvNHmEM3YAAAAWQbADAACwCIIdAACARRDsAAAALIJgZwEPPfSQRo4cWdZtAACAMkawK4f69esnm82moUOHFhuLjo6WzWZzegPFsmXL9Nprr93yviZNmqS2bduqatWq8vPzU48ePZwenSJJhw4dUs+ePeXr6yu73a4///nPxd7qcCPblaSUlBT94Q9/UJUqVWS329WhQwedP3/+pu4jAAB3MoJdORUcHKxFixY5BZcLFy5o4cKF5mu8ilSvXl1Vq1a95T2tX79e0dHR2rJli5KSkpSfn69OnTqZb8TIyclRp06dZLPZtGbNGm3atEl5eXnq1q2bCgsLb3i70qVQ17lzZ3Xq1Enbtm3T119/rZiYGLm48J8wAABFeFdsOdW6dWsdOnRIy5YtM9/isGzZMtWqVUshISFOtQ899JBatmypadOmSbr0VoohQ4bo4MGDWrp0qapVq6bx48eX6DVd15KYmOg0n5CQID8/P6WmpqpDhw7atGmTjhw5oh07dpjP2Zk3b56qVaumNWvWKCIi4oa2K0mxsbEaPny4xo4da9Y1bNjwN+0PAABWw+mOcmzAgAGaO3euOT9nzpwSv+ZqypQpatOmjXbs2KHnnntOw4YNu+LPm0USEhJks9lK1V/RmySqV68uScrNzZXNZpObm5tZ4+7uLhcXF23cuPGGt3vq1Clt3bpVfn5+at++vfz9/fXggw+WapsAAPweEOzKsaefflobN27U0aNHdfToUW3atElPP/10idbt2rWrnnvuOdWvX19jxoyRj4+P1q5de9V6Ly+vUp0BKyws1MiRI3XfffepadOmkqR7771XVapU0ZgxY3Tu3Dnl5OTohRdeUEFBgU6cOHHD2/3uu+8kSRMnTtTgwYOVmJio1q1bq2PHjjpw4ECJewYAwOpKFexmzZql5s2by263y263Kzw8XJ9//rk5fuHCBUVHR6tGjRry9PRUr169il04n5GRoaioKFWuXFl+fn4aNWqULl686FSzbt06tW7dWm5ubqpfv36xl9JL0syZM1WnTh25u7srLCxM27ZtcxovSS/lna+vr6KiopSQkKC5c+cqKipKPj4+JVq3efPm5p9tNpsCAgJ06tSpq9b37NlT+/fvL3Fv0dHR2r17txYtWuTU79KlS7Vy5Up5enrKy8tLmZmZat26dYmvhbvSdouuz/vLX/6i/v37q1WrVpo6daoaNmyoOXPmlLhnAACsrlTBrmbNmnrjjTeUmpqq7du36w9/+IO6d++uPXv2SLp0HdTKlSu1dOlSrV+/XsePH3d6iXtBQYGioqKUl5enzZs3a968eUpISFB8fLxZc/jwYUVFRenhhx9WWlqaRo4cqUGDBumLL74waxYvXqy4uDhNmDBB33zzjVq0aKHIyEin4HK9Xu4UAwYMUEJCgubNm6cBAwaUeL1KlSo5zdtstmvewFAaMTExWrVqldauXauaNWs6jXXq1EmHDh3SqVOn9NNPP+nf//63jh07prp1697wdgMDAyVJTZo0capv3LixMjIybsIeAQBgDaUKdt26dVPXrl3VoEED3X333fr73/8uT09PbdmyRVlZWfrXv/6lt99+W3/4wx8UGhqquXPnavPmzdqyZYsk6csvv9TevXv10UcfqWXLlurSpYtee+01zZw5U3l5eZKk2bNnKyQkRFOmTFHjxo0VExOjxx9/XFOnTjX7ePvttzV48GD1799fTZo00ezZs1W5cmXz7E1JerlTdO7cWXl5ecrPz1dkZGSZ9mIYhmJiYrR8+XKtWbOm2E0cl/Px8ZG3t7fWrFmjU6dO6Y9//OMNb7dOnToKCgoqdo3gt99+q9q1a/+2nQIAwEJu+Bq7goICLVq0SDk5OQoPD1dqaqry8/Od7nxs1KiRatWqpZSUFEmXHlnRrFkz+fv7mzWRkZHKzs42z/qlpKQUu3syMjLS3EZeXp5SU1OdalxcXBQREWHWlKSXO0WFChW0b98+7d27VxUqVLhp2+3YsaNmzJhhzi9fvlyNGjW65jrR0dH66KOPtHDhQlWtWlUOh0MOh8PpkSxz587Vli1bdOjQIX300Uf605/+pNjYWKfr93792dfbrs1m06hRo/TOO+/oP//5jw4ePKiXX35Z+/fv18CBA2/WVwIAwB2v1I872bVrl8LDw3XhwgV5enpq+fLlatKkidLS0uTq6ipvb2+nen9/fzkcDkmSw+FwCnVF40Vj16rJzs7W+fPn9csvv6igoOCKNUXXiDkcjuv2ciW5ubnKzc0157Ozs6/zbdweRY8OuZkOHTqkn376yZzPysq65l2z0qVrLKVLj1e53Ny5c80HJqenp2vcuHE6ffq06tSpo5deekmxsbHX/OySbHfkyJG6cOGCYmNjdfr0abVo0UJJSUmqV69eSXcZAADLK3Wwa9iwodLS0pSVlaX//Oc/6tu3r9avX38rervtJk2apFdeeaWs27jizSKXW7FihdP8unXrnOaPHDlSbJ20tLRr1vTr18/pbRZXYhjGNccl6Y033tAbb7xxzZpff3ZJtitJY8eOdXqOHQAAcFbqn2JdXV1Vv359hYaGatKkSWrRooWmT5+ugIAA5eXlKTMz06n+5MmTCggIkCQFBAQUuzO1aP56NXa7XR4eHvLx8VGFChWuWHP5Nq7Xy5WMGzdOWVlZ5vT999+X7EsBAAAoB37zc+wKCwuVm5ur0NBQVapUScnJyeZYenq6MjIyFB4eLkkKDw/Xrl27nO5eTUpKkt1uN+94DA8Pd9pGUU3RNlxdXRUaGupUU1hYqOTkZLOmJL1ciZubm/kol6IJAADgTlGqn2LHjRunLl26qFatWjpz5owWLlyodevW6YsvvpCXl5cGDhyouLg4Va9eXXa7Xc8//7zCw8N17733Srr0KIwmTZromWee0eTJk+VwODR+/HhFR0ebbysYOnSoZsyYodGjR2vAgAFas2aNlixZotWrV5t9xMXFqW/fvmrTpo3atWunadOmKScnx3wrQ0l6AQAAsJpSBbtTp07p2Wef1YkTJ+Tl5aXmzZvriy++0COPPCJJmjp1qlxcXNSrVy/l5uYqMjJS7777rrl+hQoVtGrVKg0bNkzh4eGqUqWK+vbtq1dffdWsCQkJ0erVqxUbG6vp06erZs2a+vDDD50e9dG7d2/9+OOPio+Pl8PhUMuWLZWYmOh0Q8X1egEAALAam1HSK9d/h7Kzs+Xl5aWsrCx+lgVusTpjV1+/COXOkTeiyroFwPJKk0d4V2w51K9fP9lsNg0dOrTYWHR0tGw223XvYL0dli1bpkceeUS+vr7mK+Yuf0NIkWPHjunpp59WjRo15OHhoWbNmmn79u1ONfv27dMf//hHeXl5qUqVKmrbtu013yqRn5+vV199VfXq1ZO7u7tatGihxMREp5o6derIZrMVm6Kjo2/OFwAAQDlDsCungoODtWjRIqeH/164cEELFy5UrVq1yrCz/7NhwwY98sgj+uyzz5SamqqHH35Y3bp1044dO8yaX375Rffdd58qVaqkzz//XHv37tWUKVNUrVo1s+bQoUO6//771ahRI61bt047d+7Uyy+/LHd396t+9vjx4/Xee+/pn//8p/bu3auhQ4eqZ8+eTp/99ddf68SJE+aUlJQkSfrTn/50C74NAADKHsGunGrdurWCg4O1bNkyc9myZctUq1YttWrVyqk2MTFR999/v7y9vVWjRg09+uijOnTokDk+f/58eXp66sCBA+ay5557To0aNdK5c+duuMdp06Zp9OjRatu2rRo0aKDXX39dDRo00MqVK82aN998U8HBwZo7d67atWunkJAQderUyenBwi+99JK6du2qyZMnq1WrVqpXr57++Mc/ys/P76qf/e9//1svvviiunbtqrp162rYsGHq2rWrpkyZYtb4+voqICDAnFatWqV69erpwQcfvOF9BgCgPCPYlWMDBgzQ3Llzzfk5c+aYd/5eLicnR3Fxcdq+fbuSk5Pl4uKinj17qrCwUJL07LPPqmvXrurTp48uXryo1atX68MPP9SCBQtUuXLlK372kSNHZLPZij38+FoKCwt15swZVa9e3Vz26aefqk2bNvrTn/4kPz8/tWrVSh988IHTOqtXr9bdd9+tyMhI+fn5KSwsrNhDmH8tNze32Bk9Dw8Pbdy48Yr1eXl5+uijjzRgwADZbLYS7xMAAHcSgl059vTTT2vjxo06evSojh49qk2bNunpp58uVterVy899thjql+/vlq2bKk5c+Zo165d2rt3r1nz3nvv6cSJExo+fLgGDhyoiRMnKjQ09KqfXalSJTVs2PCqwe9K/vGPf+js2bP685//bC777rvvNGvWLDVo0EBffPGFhg0bpuHDh2vevHmSLt1pffbsWb3xxhvq3LmzvvzyS/Xs2VOPPfbYNd9oEhkZqbffflsHDhxQYWGhkpKStGzZMp04ceKK9StWrFBmZma5uDYRAIBbpdSvFMPt4+vrq6ioKCUkJMgwDEVFRcnHx6dY3YEDBxQfH6+tW7fqp59+Ms/UZWRkqGnTppKkatWq6V//+pciIyPVvn37676a66677jLfvVsSCxcu1CuvvKJPPvnE6SfUwsJCtWnTRq+//rokqVWrVtq9e7dmz56tvn37mr12797dfKdsy5YttXnzZs2ePfuqP5tOnz5dgwcPVqNGjWSz2VSvXj31799fc+bMuWL9v/71L3Xp0kVBQUEl3icAAO40nLEr5wYMGKCEhATNmzdPAwYMuGJNt27ddPr0aX3wwQfaunWrtm7dKunSz4+X27BhgypUqKATJ04oJyfnpvW4aNEiDRo0SEuWLFFERITTWGBgoPlWkSKNGzc273j18fFRxYoVr1lzJb6+vlqxYoVycnJ09OhR7d+/X56enqpbt26x2qNHj+qrr77SoEGDbnQXAQC4IxDsyrnOnTsrLy9P+fn5Tg9pLvLzzz8rPT1d48ePV8eOHdW4cWP98ssvxeo2b96sN998UytXrpSnp6diYmJuSn8ff/yx+vfvr48//lhRUcWfZ3XfffcpPT3dadm3336r2rVrS7r0iri2bdtes+Za3N3dddddd+nixYv673//q+7duxermTt3rvz8/K7YHwAAVsJPseVchQoVtG/fPvPPv1atWjXVqFFD77//vgIDA5WRkVHsZ9YzZ87omWee0fDhw9WlSxfVrFlTbdu2Vbdu3fT4449LuvS6uGPHjmn+/PmSLj17rmPHjpo/f77atWt3xd4WLlyovn37avr06QoLC5PD4ZB06SYGLy8vSVJsbKzat2+v119/XX/+85+1bds2vf/++3r//ffN7YwaNUq9e/dWhw4d9PDDDysxMVErV650unHj2Wef1V133aVJkyZJkrZu3apjx46pZcuWOnbsmCZOnKjCwkKNHj3aqcfCwkLNnTtXffv2VcWK/OcOALA2ztjdAex2+1WfNO3i4qJFixYpNTVVTZs2VWxsrN566y2nmhEjRqhKlSrmdW7NmjXT66+/rr/85S86duyYJOnEiRNOP33m5+crPT39mo9Def/993Xx4kVFR0crMDDQnEaMGGHWtG3bVsuXL9fHH3+spk2b6rXXXtO0adPUp08fs6Znz56aPXu2Jk+erGbNmunDDz/Uf//7X91///1mTUZGhtONERcuXND48ePVpEkT9ezZU3fddZc2btwob29vpx6/+uorZWRkXPVnbAAArIRXil0DrxQDbh9eKXZn4pViwK3HK8UAAAB+hwh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIkoV7CZNmqS2bduqatWq8vPzU48ePZSenu5U89BDD8lmszlNQ4cOdarJyMhQVFSUKleuLD8/P40aNUoXL150qlm3bp1at24tNzc31a9fXwkJCcX6mTlzpurUqSN3d3eFhYVp27ZtTuMXLlxQdHS0atSoIU9PT/Xq1UsnT54szS4DAADcMUoV7NavX6/o6Ght2bJFSUlJys/PV6dOnZSTk+NUN3jwYJ04ccKcJk+ebI4VFBQoKipKeXl52rx5s+bNm6eEhATFx8ebNYcPH1ZUVJQefvhhpaWlaeTIkRo0aJC++OILs2bx4sWKi4vThAkT9M0336hFixaKjIzUqVOnzJrY2FitXLlSS5cu1fr163X8+HE99thjpf6SAAAA7gQ2wzCMG135xx9/lJ+fn9avX68OHTpIunTGrmXLlpo2bdoV1/n888/16KOP6vjx4/L395ckzZ49W2PGjNGPP/4oV1dXjRkzRqtXr9bu3bvN9Z544gllZmYqMTFRkhQWFqa2bdtqxowZkqTCwkIFBwfr+eef19ixY5WVlSVfX18tXLhQjz/+uCRp//79aty4sVJSUnTvvfded/+ys7Pl5eWlrKws2e32G/2aAJRAnbGry7oF3IAjb0SVdQuA5ZUmj/yma+yysrIkSdWrV3davmDBAvn4+Khp06YaN26czp07Z46lpKSoWbNmZqiTpMjISGVnZ2vPnj1mTUREhNM2IyMjlZKSIknKy8tTamqqU42Li4siIiLMmtTUVOXn5zvVNGrUSLVq1TJrfi03N1fZ2dlOEwAAwJ2i4o2uWFhYqJEjR+q+++5T06ZNzeVPPfWUateuraCgIO3cuVNjxoxRenq6li1bJklyOBxOoU6SOe9wOK5Zk52drfPnz+uXX35RQUHBFWv2799vbsPV1VXe3t7Faoo+59cmTZqkV155pZTfBAAAQPlww8EuOjpau3fv1saNG52WDxkyxPxzs2bNFBgYqI4dO+rQoUOqV6/ejXd6G4wbN05xcXHmfHZ2toKDg8uwIwAAgJK7oZ9iY2JitGrVKq1du1Y1a9a8Zm1YWJgk6eDBg5KkgICAYnemFs0HBARcs8Zut8vDw0M+Pj6qUKHCFWsu30ZeXp4yMzOvWvNrbm5ustvtThMAAMCdolTBzjAMxcTEaPny5VqzZo1CQkKuu05aWpokKTAwUJIUHh6uXbt2Od29mpSUJLvdriZNmpg1ycnJTttJSkpSeHi4JMnV1VWhoaFONYWFhUpOTjZrQkNDValSJaea9PR0ZWRkmDUAAABWUqqfYqOjo7Vw4UJ98sknqlq1qnmtmpeXlzw8PHTo0CEtXLhQXbt2VY0aNbRz507FxsaqQ4cOat68uSSpU6dOatKkiZ555hlNnjxZDodD48ePV3R0tNzc3CRJQ4cO1YwZMzR69GgNGDBAa9as0ZIlS7R69f/dNRcXF6e+ffuqTZs2ateunaZNm6acnBz179/f7GngwIGKi4tT9erVZbfb9fzzzys8PLxEd8QCAADcaUoV7GbNmiXp0iNNLjd37lz169dPrq6u+uqrr8yQFRwcrF69emn8+PFmbYUKFbRq1SoNGzZM4eHhqlKlivr27atXX33VrAkJCdHq1asVGxur6dOnq2bNmvrwww8VGRlp1vTu3Vs//vij4uPj5XA41LJlSyUmJjrdUDF16lS5uLioV69eys3NVWRkpN59991SfUEAAAB3it/0HDur4zl2wO3Dc+zuTDzHDrj1bttz7AAAAFB+EOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEaUKdpMmTVLbtm1VtWpV+fn5qUePHkpPT3equXDhgqKjo1WjRg15enqqV69eOnnypFNNRkaGoqKiVLlyZfn5+WnUqFG6ePGiU826devUunVrubm5qX79+kpISCjWz8yZM1WnTh25u7srLCxM27ZtK3UvAAAAVlGqYLd+/XpFR0dry5YtSkpKUn5+vjp16qScnByzJjY2VitXrtTSpUu1fv16HT9+XI899pg5XlBQoKioKOXl5Wnz5s2aN2+eEhISFB8fb9YcPnxYUVFRevjhh5WWlqaRI0dq0KBB+uKLL8yaxYsXKy4uThMmTNA333yjFi1aKDIyUqdOnSpxLwAAAFZiMwzDuNGVf/zxR/n5+Wn9+vXq0KGDsrKy5Ovrq4ULF+rxxx+XJO3fv1+NGzdWSkqK7r33Xn3++ed69NFHdfz4cfn7+0uSZs+erTFjxujHH3+Uq6urxowZo9WrV2v37t3mZz3xxBPKzMxUYmKiJCksLExt27bVjBkzJEmFhYUKDg7W888/r7Fjx5aol+vJzs6Wl5eXsrKyZLfbb/RrAlACdcauLusWcAOOvBFV1i0AlleaPPKbrrHLysqSJFWvXl2SlJqaqvz8fEVERJg1jRo1Uq1atZSSkiJJSklJUbNmzcxQJ0mRkZHKzs7Wnj17zJrLt1FUU7SNvLw8paamOtW4uLgoIiLCrClJLwAAAFZS8UZXLCws1MiRI3XfffepadOmkiSHwyFXV1d5e3s71fr7+8vhcJg1l4e6ovGisWvVZGdn6/z58/rll19UUFBwxZr9+/eXuJdfy83NVW5urjmfnZ19va8BAACg3LjhM3bR0dHavXu3Fi1adDP7KVOTJk2Sl5eXOQUHB5d1SwAAACV2Q8EuJiZGq1at0tq1a1WzZk1zeUBAgPLy8pSZmelUf/LkSQUEBJg1v74ztWj+ejV2u10eHh7y8fFRhQoVrlhz+Tau18uvjRs3TllZWeb0/fffl+DbAAAAKB9KFewMw1BMTIyWL1+uNWvWKCQkxGk8NDRUlSpVUnJysrksPT1dGRkZCg8PlySFh4dr165dTnevJiUlyW63q0mTJmbN5dsoqinahqurq0JDQ51qCgsLlZycbNaUpJdfc3Nzk91ud5oAAADuFKW6xi46OloLFy7UJ598oqpVq5rXqnl5ecnDw0NeXl4aOHCg4uLiVL16ddntdj3//PMKDw8370Lt1KmTmjRpomeeeUaTJ0+Ww+HQ+PHjFR0dLTc3N0nS0KFDNWPGDI0ePVoDBgzQmjVrtGTJEq1e/X93zcXFxalv375q06aN2rVrp2nTpiknJ0f9+/c3e7peLwAAAFZSqmA3a9YsSdJDDz3ktHzu3Lnq16+fJGnq1KlycXFRr169lJubq8jISL377rtmbYUKFbRq1SoNGzZM4eHhqlKlivr27atXX33VrAkJCdHq1asVGxur6dOnq2bNmvrwww8VGRlp1vTu3Vs//vij4uPj5XA41LJlSyUmJjrdUHG9XgAAAKzkNz3Hzup4jh1w+/AcuzsTz7EDbr3b9hw7AAAAlB8EOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEVy7oBAMCdq87Y1WXdwi1z5I2osm4BKDXO2AEAAFgEwQ4AAMAiSh3sNmzYoG7duikoKEg2m00rVqxwGu/Xr59sNpvT1LlzZ6ea06dPq0+fPrLb7fL29tbAgQN19uxZp5qdO3fqgQcekLu7u4KDgzV58uRivSxdulSNGjWSu7u7mjVrps8++8xp3DAMxcfHKzAwUB4eHoqIiNCBAwdKu8sAAAB3hFIHu5ycHLVo0UIzZ868ak3nzp114sQJc/r444+dxvv06aM9e/YoKSlJq1at0oYNGzRkyBBzPDs7W506dVLt2rWVmpqqt956SxMnTtT7779v1mzevFlPPvmkBg4cqB07dqhHjx7q0aOHdu/ebdZMnjxZ77zzjmbPnq2tW7eqSpUqioyM1IULF0q72wAAAOWezTAM44ZXttm0fPly9ejRw1zWr18/ZWZmFjuTV2Tfvn1q0qSJvv76a7Vp00aSlJiYqK5du+qHH35QUFCQZs2apZdeekkOh0Ourq6SpLFjx2rFihXav3+/JKl3797KycnRqlWrzG3fe++9atmypWbPni3DMBQUFKS//vWveuGFFyRJWVlZ8vf3V0JCgp544onr7l92dra8vLyUlZUlu91+I18RgBKy8kX4uDNx8wTKi9LkkVtyjd26devk5+enhg0batiwYfr555/NsZSUFHl7e5uhTpIiIiLk4uKirVu3mjUdOnQwQ50kRUZGKj09Xb/88otZExER4fS5kZGRSklJkSQdPnxYDofDqcbLy0thYWFmza/l5uYqOzvbaQIAALhT3PRg17lzZ82fP1/Jycl68803tX79enXp0kUFBQWSJIfDIT8/P6d1KlasqOrVq8vhcJg1/v7+TjVF89eruXz88vWuVPNrkyZNkpeXlzkFBweXev8BAADKyk1/jt3lP3E2a9ZMzZs3V7169bRu3Tp17NjxZn/cTTVu3DjFxcWZ89nZ2YQ7AABwx7jljzupW7eufHx8dPDgQUlSQECATp065VRz8eJFnT59WgEBAWbNyZMnnWqK5q9Xc/n45etdqebX3NzcZLfbnSYAAIA7xS0Pdj/88IN+/vlnBQYGSpLCw8OVmZmp1NRUs2bNmjUqLCxUWFiYWbNhwwbl5+ebNUlJSWrYsKGqVatm1iQnJzt9VlJSksLDwyVJISEhCggIcKrJzs7W1q1bzRoAAAArKXWwO3v2rNLS0pSWlibp0k0KaWlpysjI0NmzZzVq1Cht2bJFR44cUXJysrp376769esrMjJSktS4cWN17txZgwcP1rZt27Rp0ybFxMToiSeeUFBQkCTpqaeekqurqwYOHKg9e/Zo8eLFmj59utPPpCNGjFBiYqKmTJmi/fv3a+LEidq+fbtiYmIkXbpjd+TIkfrb3/6mTz/9VLt27dKzzz6roKAgp7t4AQAArKLU19ht375dDz/8sDlfFLb69u2rWbNmaefOnZo3b54yMzMVFBSkTp066bXXXpObm5u5zoIFCxQTE6OOHTvKxcVFvXr10jvvvGOOe3l56csvv1R0dLRCQ0Pl4+Oj+Ph4p2fdtW/fXgsXLtT48eP14osvqkGDBlqxYoWaNm1q1owePVo5OTkaMmSIMjMzdf/99ysxMVHu7u6l3W0AAIBy7zc9x87qeI4dcPvwHDuUNzzHDuVFmT/HDgAAALcfwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFhExbJuAEDp1Bm7uqxbAACUU5yxAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsIhSB7sNGzaoW7duCgoKks1m04oVK5zGDcNQfHy8AgMD5eHhoYiICB04cMCp5vTp0+rTp4/sdru8vb01cOBAnT171qlm586deuCBB+Tu7q7g4GBNnjy5WC9Lly5Vo0aN5O7urmbNmumzzz4rdS8AAABWUepgl5OToxYtWmjmzJlXHJ88ebLeeecdzZ49W1u3blWVKlUUGRmpCxcumDV9+vTRnj17lJSUpFWrVmnDhg0aMmSIOZ6dna1OnTqpdu3aSk1N1VtvvaWJEyfq/fffN2s2b96sJ598UgMHDtSOHTvUo0cP9ejRQ7t37y5VLwAAAFZhMwzDuOGVbTYtX75cPXr0kHTpDFlQUJD++te/6oUXXpAkZWVlyd/fXwkJCXriiSe0b98+NWnSRF9//bXatGkjSUpMTFTXrl31ww8/KCgoSLNmzdJLL70kh8MhV1dXSdLYsWO1YsUK7d+/X5LUu3dv5eTkaNWqVWY/9957r1q2bKnZs2eXqJfryc7OlpeXl7KysmS322/0awJuqjpjV5d1C8DvwpE3osq6BUBS6fLITb3G7vDhw3I4HIqIiDCXeXl5KSwsTCkpKZKklJQUeXt7m6FOkiIiIuTi4qKtW7eaNR06dDBDnSRFRkYqPT1dv/zyi1lz+ecU1RR9Tkl6AQAAsJKKN3NjDodDkuTv7++03N/f3xxzOBzy8/NzbqJiRVWvXt2pJiQkpNg2isaqVasmh8Nx3c+5Xi+/lpubq9zcXHM+Ozv7OnsMAABQfnBX7GUmTZokLy8vcwoODi7rlgAAAErspga7gIAASdLJkyedlp88edIcCwgI0KlTp5zGL168qNOnTzvVXGkbl3/G1WouH79eL782btw4ZWVlmdP3339fgr0GAAAoH25qsAsJCVFAQICSk5PNZdnZ2dq6davCw8MlSeHh4crMzFRqaqpZs2bNGhUWFiosLMys2bBhg/Lz882apKQkNWzYUNWqVTNrLv+copqizylJL7/m5uYmu93uNAEAANwpSh3szp49q7S0NKWlpUm6dJNCWlqaMjIyZLPZNHLkSP3tb3/Tp59+ql27dunZZ59VUFCQeeds48aN1blzZw0ePFjbtm3Tpk2bFBMToyeeeEJBQUGSpKeeekqurq4aOHCg9uzZo8WLF2v69OmKi4sz+xgxYoQSExM1ZcoU7d+/XxMnTtT27dsVExMjSSXqBQAAwEpKffPE9u3b9fDDD5vzRWGrb9++SkhI0OjRo5WTk6MhQ4YoMzNT999/vxITE+Xu7m6us2DBAsXExKhjx45ycXFRr1699M4775jjXl5e+vLLLxUdHa3Q0FD5+PgoPj7e6Vl37du318KFCzV+/Hi9+OKLatCggVasWKGmTZuaNSXpBQAAwCp+03PsrI7n2KE84jl2wO3Bc+xQXpTZc+wAAABQdgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEXc9GA3ceJE2Ww2p6lRo0bm+IULFxQdHa0aNWrI09NTvXr10smTJ522kZGRoaioKFWuXFl+fn4aNWqULl686FSzbt06tW7dWm5ubqpfv74SEhKK9TJz5kzVqVNH7u7uCgsL07Zt22727gIAAJQbt+SM3T333KMTJ06Y08aNG82x2NhYrVy5UkuXLtX69et1/PhxPfbYY+Z4QUGBoqKilJeXp82bN2vevHlKSEhQfHy8WXP48GFFRUXp4YcfVlpamkaOHKlBgwbpiy++MGsWL16suLg4TZgwQd98841atGihyMhInTp16lbsMgAAQJmzGYZh3MwNTpw4UStWrFBaWlqxsaysLPn6+mrhwoV6/PHHJUn79+9X48aNlZKSonvvvVeff/65Hn30UR0/flz+/v6SpNmzZ2vMmDH68ccf5erqqjFjxmj16tXavXu3ue0nnnhCmZmZSkxMlCSFhYWpbdu2mjFjhiSpsLBQwcHBev755zV27NgS7Ut2dra8vLyUlZUlu93+W74W4KapM3Z1WbcA/C4ceSOqrFsAJJUuj9ySM3YHDhxQUFCQ6tatqz59+igjI0OSlJqaqvz8fEVERJi1jRo1Uq1atZSSkiJJSklJUbNmzcxQJ0mRkZHKzs7Wnj17zJrLt1FUU7SNvLw8paamOtW4uLgoIiLCrLmS3NxcZWdnO00AAAB3ipse7MLCwpSQkKDExETNmjVLhw8f1gMPPKAzZ87I4XDI1dVV3t7eTuv4+/vL4XBIkhwOh1OoKxovGrtWTXZ2ts6fP6+ffvpJBQUFV6wp2saVTJo0SV5eXuYUHBx8Q98BAABAWah4szfYpUsX88/NmzdXWFiYateurSVLlsjDw+Nmf9xNNW7cOMXFxZnz2dnZhDsAAHDHuOWPO/H29tbdd9+tgwcPKiAgQHl5ecrMzHSqOXnypAICAiRJAQEBxe6SLZq/Xo3dbpeHh4d8fHxUoUKFK9YUbeNK3NzcZLfbnSYAAIA7xS0PdmfPntWhQ4cUGBio0NBQVapUScnJyeZ4enq6MjIyFB4eLkkKDw/Xrl27nO5eTUpKkt1uV5MmTcyay7dRVFO0DVdXV4WGhjrVFBYWKjk52awBAACwmpse7F544QWtX79eR44c0ebNm9WzZ09VqFBBTz75pLy8vDRw4EDFxcVp7dq1Sk1NVf/+/RUeHq57771XktSpUyc1adJEzzzzjP73v//piy++0Pjx4xUdHS03NzdJ0tChQ/Xdd99p9OjR2r9/v959910tWbJEsbGxZh9xcXH64IMPNG/ePO3bt0/Dhg1TTk6O+vfvf7N3GQAAoFy46dfY/fDDD3ryySf1888/y9fXV/fff7+2bNkiX19fSdLUqVPl4uKiXr16KTc3V5GRkXr33XfN9StUqKBVq1Zp2LBhCg8PV5UqVdS3b1+9+uqrZk1ISIhWr16t2NhYTZ8+XTVr1tSHH36oyMhIs6Z379768ccfFR8fL4fDoZYtWyoxMbHYDRUAAABWcdOfY2clPMcO5RHPsQNuD55jh/KizJ9jBwAAgNuPYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBAABYBMEOAADAIgh2AAAAFlGxrBsAAKA8qjN2dVm3cEsceSOqrFvALcQZOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBE8oBiWZNUHiwIAcC2csQMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAsgmAHAABgEQQ7AAAAiyDYAQAAWATBDgAAwCIIdgAAABZBsAMAALAIgh0AAIBFEOwAAAAs4ncR7GbOnKk6derI3d1dYWFh2rZtW1m3BAAAcNNZPtgtXrxYcXFxmjBhgr755hu1aNFCkZGROnXqVFm3BgAAcFPZDMMwyrqJWyksLExt27bVjBkzJEmFhYUKDg7W888/r7Fjx15z3ezsbHl5eSkrK0t2u/12tIubpM7Y1WXdAgDgNjvyRlRZt3BLlCaPVLxNPZWJvLw8paamaty4ceYyFxcXRUREKCUlpQw7Kx8IPwAAWIulg91PP/2kgoIC+fv7Oy339/fX/v37i9Xn5uYqNzfXnM/KypJ0KSlbUWHuubJuAQCAm8aq/78u2q+S/Mhq6WBXWpMmTdIrr7xSbHlwcHAZdAMAAErDa1pZd3BrnTlzRl5eXtessXSw8/HxUYUKFXTy5Emn5SdPnlRAQECx+nHjxikuLs6cLyws1OnTp1WjRg3ZbLZb3u/Nlp2dreDgYH3//fdcI3iH41haC8fTOjiW1lGej6VhGDpz5oyCgoKuW2vpYOfq6qrQ0FAlJyerR48eki6FteTkZMXExBSrd3Nzk5ubm9Myb2/v29DprWW328vdf6S4MRxLa+F4WgfH0jrK67G83pm6IpYOdpIUFxenvn37qk2bNmrXrp2mTZumnJwc9e/fv6xbAwAAuKksH+x69+6tH3/8UfHx8XI4HGrZsqUSExOL3VABAABwp7N8sJOkmJiYK/70anVubm6aMGFCsZ+XcefhWFoLx9M6OJbWYZVjafkHFAMAAPxeWP6VYgAAAL8XBDsAAACLINgBAABYBMHuDjdx4kTZbDanqVGjRub4hQsXFB0drRo1asjT01O9evUq9sBmlJ0NGzaoW7duCgoKks1m04oVK5zGDcNQfHy8AgMD5eHhoYiICB04cMCp5vTp0+rTp4/sdru8vb01cOBAnT179jbuBaTrH8t+/foV+7vauXNnpxqOZfkwadIktW3bVlWrVpWfn5969Oih9PR0p5qS/NuakZGhqKgoVa5cWX5+fho1apQuXrx4O3fld68kx/Khhx4q9ndz6NChTjV30rEk2FnAPffcoxMnTpjTxo0bzbHY2FitXLlSS5cu1fr163X8+HE99thjZdgtLpeTk6MWLVpo5syZVxyfPHmy3nnnHc2ePVtbt25VlSpVFBkZqQsXLpg1ffr00Z49e5SUlKRVq1Zpw4YNGjJkyO3aBfz/rncsJalz585Of1c//vhjp3GOZfmwfv16RUdHa8uWLUpKSlJ+fr46deqknJwcs+Z6/7YWFBQoKipKeXl52rx5s+bNm6eEhATFx8eXxS79bpXkWErS4MGDnf5uTp482Ry7446lgTvahAkTjBYtWlxxLDMz06hUqZKxdOlSc9m+ffsMSUZKSspt6hAlJclYvny5OV9YWGgEBAQYb731lrksMzPTcHNzMz7++GPDMAxj7969hiTj66+/Nms+//xzw2azGceOHbttvcPZr4+lYRhG3759je7du191HY5l+XXq1ClDkrF+/XrDMEr2b+tnn31muLi4GA6Hw6yZNWuWYbfbjdzc3Nu7AzD9+lgahmE8+OCDxogRI666zp12LDljZwEHDhxQUFCQ6tatqz59+igjI0OSlJqaqvz8fEVERJi1jRo1Uq1atZSSklJW7aKEDh8+LIfD4XT8vLy8FBYWZh6/lJQUeXt7q02bNmZNRESEXFxctHXr1tveM65t3bp18vPzU8OGDTVs2DD9/PPP5hjHsvzKysqSJFWvXl1Syf5tTUlJUbNmzZwehh8ZGans7Gzt2bPnNnaPy/36WBZZsGCBfHx81LRpU40bN07nzp0zx+60Y/m7eECxlYWFhSkhIUENGzbUiRMn9Morr+iBBx7Q7t275XA45OrqWux9t/7+/nI4HGXTMEqs6Bj9+i0plx8/h8MhPz8/p/GKFSuqevXqHONypnPnznrssccUEhKiQ4cO6cUXX1SXLl2UkpKiChUqcCzLqcLCQo0cOVL33XefmjZtKkkl+rfV4XBc8e9u0RhuvysdS0l66qmnVLt2bQUFBWnnzp0aM2aM0tPTtWzZMkl33rEk2N3hunTpYv65efPmCgsLU+3atbVkyRJ5eHiUYWcALvfEE0+Yf27WrJmaN2+uevXqad26derYsWMZdoZriY6O1u7du52uXcad6WrH8vLrWJs1a6bAwEB17NhRhw4dUr169W53m78ZP8VajLe3t+6++24dPHhQAQEBysvLU2ZmplPNyZMnFRAQUDYNosSKjtGv77S7/PgFBATo1KlTTuMXL17U6dOnOcblXN26deXj46ODBw9K4liWRzExMVq1apXWrl2rmjVrmstL8m9rQEDAFf/uFo3h9rrasbySsLAwSXL6u3knHUuCncWcPXtWhw4dUmBgoEJDQ1WpUiUlJyeb4+np6crIyFB4eHgZdomSCAkJUUBAgNPxy87O1tatW83jFx4erszMTKWmppo1a9asUWFhofmPE8qnH374QT///LMCAwMlcSzLE8MwFBMTo+XLl2vNmjUKCQlxGi/Jv63h4eHatWuXU1hPSkqS3W5XkyZNbs+O4LrH8krS0tIkyenv5h11LMv67g38Nn/961+NdevWGYcPHzY2bdpkREREGD4+PsapU6cMwzCMoUOHGrVq1TLWrFljbN++3QgPDzfCw8PLuGsUOXPmjLFjxw5jx44dhiTj7bffNnbs2GEcPXrUMAzDeOONNwxvb2/jk08+MXbu3Gl0797dCAkJMc6fP29uo3PnzkarVq2MrVu3Ghs3bjQaNGhgPPnkk2W1S79b1zqWZ86cMV544QUjJSXFOHz4sPHVV18ZrVu3Nho0aGBcuHDB3AbHsnwYNmyY4eXlZaxbt844ceKEOZ07d86sud6/rRcvXjSaNm1qdOrUyUhLSzMSExMNX19fY9y4cWWxS79b1zuWBw8eNF599VVj+/btxuHDh41PPvnEqFu3rtGhQwdzG3fasSTY3eF69+5tBAYGGq6ursZdd91l9O7d2zh48KA5fv78eeO5554zqlWrZlSuXNno2bOnceLEiTLsGJdbu3atIanY1LdvX8MwLj3y5OWXXzb8/f0NNzc3o2PHjkZ6errTNn7++WfjySefNDw9PQ273W7079/fOHPmTBnsze/btY7luXPnjE6dOhm+vr5GpUqVjNq1axuDBw92enyCYXAsy4srHUdJxty5c82akvzbeuTIEaNLly6Gh4eH4ePjY/z1r3818vPzb/Pe/L5d71hmZGQYHTp0MKpXr264ubkZ9evXN0aNGmVkZWU5bedOOpY2wzCM23d+EAAAALcK19gBAABYBMEOAADAIgh2AAAAFkGwAwAAsAiCHQAAgEUQ7AAAACyCYAcAAGARBDsAAACLINgBwC1Wp04dTZs2razbAPA7QLADgFvs66+/1pAhQ0pUaxiG4uPjFRgYKA8PD0VEROjAgQO3uEMAVkGwA4BbzNfXV5UrVy5R7eTJk/XOO+9o9uzZ2rp1q6pUqaLIyEhduHDhFncJwAoIdgBwHYmJibr//vvl7e2tGjVq6NFHH9WhQ4ckSfPnz5enp6fTWbXnnntOjRo10rlz5yQ5/xRrGIYmTpyoWrVqyc3NTUFBQRo+fLg5Nm3aNI0fP17du3dX8+bNNX/+fB0/flwrVqy4rfsM4M5EsAOA68jJyVFcXJy2b9+u5ORkubi4qGfPniosLNSzzz6rrl27qk+fPrp48aJWr16tDz/8UAsWLLjiWbr//ve/mjp1qt577z0dOHBAK1asULNmzSRJhw8flsPhUEREhFnv5eWlsLAwpaSk3Lb9BXDnqljWDQBAederVy+n+Tlz5sjX11d79+5V06ZN9d5776l58+YaPny4li1bpokTJyo0NPSK28rIyFBAQIAiIiJUqVIl1apVS+3atZMkORwOSZK/v7/TOv7+/uYYAFwLZ+wA4DoOHDigJ598UnXr1pXdbledOnUkXQppklStWjX961//0qxZs1SvXj2NHTv2qtv605/+pPPnz6tu3boaPHiwli9frosXL96O3QDwO0CwA4Dr6Natm06fPq0PPvhAW7du1datWyVJeXl5Zs2GDRtUoUIFnThxQjk5OVfdVnBwsNLT0/Xuu+/Kw8NDzz33nDp06KD8/HwFBARIkk6ePOm0zsmTJ80xALgWgh0AXMPPP/+s9PR0jR8/Xh07dlTjxo31yy+/ONVs3rxZb775plauXClPT0/FxMRcc5seHh7q1q2b3nnnHa1bt04pKSnatWuXQkJCFBAQoOTkZLM2OztbW7duVXh4+C3ZPwDWwjV2AHAN1apVU40aNfT+++8rMDBQGRkZTj+1njlzRs8884yGDx+uLl26qGbNmmrbtq26deumxx9/vNj2EhISVFBQoLCwMFWuXFkfffSRPDw8VLt2bdlsNo0cOVJ/+9vf1KBBA4WEhOjll19WUFCQevTocRv3GsCdijN2AHANLi4uWrRokVJTU9W0aVPFxsbqrbfeMsdHjBihKlWq6PXXX5ckNWvWTK+//rr+8pe/6NixY8W25+3trQ8++ED33Xefmjdvrq+++korV65UjRo1JEmjR4/W888/ryFDhqht27Y6e/asEhMT5e7ufnt2GMAdzWYYhlHWTQAAAOC344wdAACARRDsAAAALIJgBwAAYBEEOwAAAIsg2AEAAFgEwQ4AAMAiCHYAAAAWQbADAACwCIIdAACARRDsAAAALIJgBwAAYBEEOwAAAIv4/wC+rQNucIktjAAAAABJRU5ErkJggg==\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h = physt.h1(data)\n", "h.plot(show_stats=True)\n", "h.statistics" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA0U0lEQVR4nO3deXhU1f3H8c8QYBLIBgmQhIQQkFUWUWgaWVuiEBBxqUWKGJFClaQqFAupZbNCqLaIrYCoRaRFU/nJoqJYWQKCIbIKggaIwVghUKFZ2AZIzu8Pf8zPIQEyQBI4eb+e5z4P955z7/3eOZPh89x7547DGGMEAACA616Nqi4AAAAAVwfBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAElUW7NatW6cBAwYoIiJCDodDS5cu9Wr9yZMny+FwlJrq1q1bMQUDAABc46os2B0/flwdO3bUrFmzLmv9sWPH6uDBgx5T27Ztdd99913lSgEAAK4PVRbsEhIS9Mwzz+juu+8us93lcmns2LFq3Lix6tatq9jYWKWnp7vb/f39FRYW5p4OHTqk3bt3a/jw4ZV0BAAAANeWa/Yeu+TkZGVkZCgtLU07duzQfffdp759+2rv3r1l9n/11VfVsmVLde/evZIrBQAAuDZck8EuNzdXr732mhYtWqTu3burefPmGjt2rLp166bXXnutVP9Tp05p4cKFnK0DAADVWs2qLqAsO3fuVHFxsVq2bOmx3OVyKSQkpFT/JUuWqKioSImJiZVVIgAAwDXnmgx2x44dk4+Pj7Zs2SIfHx+PNn9//1L9X331Vd1xxx1q1KhRZZUIAABwzbkmg12nTp1UXFysw4cPX/KeuZycHK1Zs0bvvPNOJVUHAABwbaqyYHfs2DHt27fPPZ+Tk6Pt27erfv36atmypYYMGaIHH3xQf/7zn9WpUyf95z//0apVq9ShQwf179/fvd68efMUHh6uhISEqjgMAACAa4bDGGOqYsfp6en6yU9+Ump5YmKi5s+frzNnzuiZZ57RggUL9O233yo0NFQ//vGPNWXKFLVv316SVFJSoujoaD344IOaOnVqZR8CAADANaXKgh0AAACurmvycScAAADwHsEOAADAEpX+5YmSkhIdOHBAAQEBcjgclb17AACA64oxRkVFRYqIiFCNGhc/J1fpwe7AgQOKioqq7N0CAABc17755htFRkZetE+lB7uAgABJ3xcXGBhY2bsHAAC4rhQWFioqKsqdoS6m0oPducuvgYGBBDsAAIByKs8tbHx5AgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBLsK5nA4tHTp0qouAwAAVAPVJtg5HI6LTpMnT77guvv375fD4dD27dsrpdapU6fq1ltvVZ06dRQcHHzRvkeOHFFkZKQcDofy8/M92mbNmqU2bdrIz89PrVq10oIFCzzaz5w5o6efflrNmzeXr6+vOnbsqBUrVpTax7fffqsHHnhAISEh8vPzU/v27bV582Z3uzFGEydOVHh4uPz8/BQfH6+9e/d6bOPo0aMaMmSIAgMDFRwcrOHDh+vYsWMefXbs2KHu3bvL19dXUVFRevbZZ0vVkp+fr6SkJIWHh8vpdKply5Z6//33L/oalWe7M2fOVKtWreTn56eoqCiNHj1ap06dcrcXFxdrwoQJiomJkZ+fn5o3b64//OEPMsZ4bOeLL77QnXfeqaCgINWtW1ddunRRbm5uqf0ZY5SQkFAq+B85ckR9+/ZVRESEnE6noqKilJycrMLCQo/1XS6XnnrqKUVHR8vpdKpp06aaN2/eRV8HAEA1YCpZQUGBkWQKCgoqdb8HDx50TzNnzjSBgYEey4qKii64bk5OjpFktm3b5vV+JZklS5Z4tc7EiRPNjBkzzJgxY0xQUNBF+w4cONAkJCQYSea///2ve/ns2bNNQECASUtLM9nZ2ebNN980/v7+5p133nH3+e1vf2siIiLM8uXLTXZ2tpk9e7bx9fU1W7dudfc5evSoiY6ONg899JDJzMw0X331lfnwww/Nvn373H2mT59ugoKCzNKlS81nn31m7rzzThMTE2NOnjzp7tO3b1/TsWNHs3HjRvPxxx+bG264wQwePNjdXlBQYBo1amSGDBliPv/8c/Pmm28aPz8/M3fuXHcfl8tlOnfubPr162fWr19vcnJyTHp6utm+ffsFX5/ybHfhwoXG6XSahQsXmpycHPPhhx+a8PBwM3r0aHefqVOnmpCQEPPee++ZnJwcs2jRIuPv729eeOEFd599+/aZ+vXrmyeffNJs3brV7Nu3zyxbtswcOnSoVF0zZsxwj9sP3x9Hjx41s2fPNps2bTL79+83K1euNK1atfJ4rYwx5s477zSxsbHmo48+Mjk5OeaTTz4x69evv+DrAAC4fnmTnapNsPuh1157zSMwFRcXmylTppjGjRub2rVrm44dO5oPPvjA3S7JY+rZs6cxxphPP/3UxMfHm5CQEBMYGGh69OhhtmzZ4rGvywl2F6rzfLNnzzY9e/Y0q1atKhXs4uLizNixYz36jxkzxnTt2tU9Hx4ebl588UWPPvfcc48ZMmSIe37cuHGmW7duF6yhpKTEhIWFmeeee869LD8/3zidTvPmm28aY4zZvXu3kWQ2bdrk7vPBBx8Yh8Nhvv32W/ex1KtXz7hcLo99t2rVyj0/Z84c06xZM3P69OkL1nO+8mw3KSnJ/PSnP/VY7/zXqn///ubhhx/26HP+azVo0CDzwAMPXLKmbdu2mcaNG5uDBw+W6/3xwgsvmMjISPf8Bx98YIKCgsyRI0cuuS8AwPXPm+xUbS7FXswLL7ygP//5z/rTn/6kHTt2qE+fPrrzzjvdlxM//fRTSdLKlSt18OBBLV68WJJUVFSkxMRErV+/Xhs3blSLFi3Ur18/FRUVXXBfvXr10kMPPXTFNe/evVtPP/20FixYoBo1Sg+jy+WSr6+vxzI/Pz99+umnOnPmzEX7rF+/3j3/zjvvqHPnzrrvvvvUsGFDderUSa+88oq7PScnR3l5eYqPj3cvCwoKUmxsrDIyMiRJGRkZCg4OVufOnd194uPjVaNGDWVmZrr79OjRQ7Vr13b36dOnj7KysvTf//7XXUtcXJySkpLUqFEjtWvXTtOmTVNxcfEFX6fybPfWW2/Vli1b3OP81Vdf6f3331e/fv3c69x6661atWqV9uzZI0n67LPPtH79eiUkJEiSSkpKtHz5crVs2VJ9+vRRw4YNFRsbW+r+yhMnTugXv/iFZs2apbCwsAvWfc6BAwe0ePFi9ezZ073s3Jg8++yzaty4sVq2bKmxY8fq5MmTl9weAMBylRA0PVyLZ+wiIiLM1KlTPfp06dLFjBo1yhhT/kuxxcXFJiAgwLz77rvuZTrvjMzQoUPN+PHjL6vOc06dOmU6dOhg/v73vxtjjFmzZk2pM3YpKSkmLCzMbN682ZSUlJhNmzaZRo0aGUnmwIEDxhhjBg8ebNq2bWv27NljiouLzb/+9S/j5+dnateu7d6O0+k0TqfTpKSkmK1bt5q5c+caX19fM3/+fGOMMRs2bPDY5jn33Xef+fnPf26M+f4yZsuWLUsdR4MGDczs2bONMcbcdtttZuTIkR7tu3btMpLM7t27jTHGtGrVyjidTvPwww+bzZs3m7S0NFO/fn0zefLkC76G5dmuMd+fFatVq5apWbOmkWQeeeQRj3WKi4vNuHHjjMPhMDVr1jQOh8NMmzbN3X7u7FudOnXMjBkzzLZt20xqaqpxOBwmPT3d3W/kyJFm+PDh7vnz3x/n3H///cbPz89IMgMGDPC4rN2nTx/jdDpN//79TWZmplm+fLn7cjkAwD6csfNCYWGhDhw4oK5du3os79q1q7744ouLrnvo0CGNGDFCLVq0UFBQkAIDA3Xs2LEyb5Y/Z8GCBUpNTb2imlNSUtSmTRs98MADF+wzYcIEJSQk6Mc//rFq1aqlgQMHKjExUZLcZ/heeOEFtWjRQq1bt1bt2rWVnJysYcOGeZwBLCkp0c0336xp06apU6dOGjlypEaMGKGXXnrpio7hcpSUlKhhw4Z6+eWXdcstt2jQoEF66qmnrriW9PR0TZs2TbNnz9bWrVu1ePFiLV++XH/4wx/cfd566y0tXLhQb7zxhrZu3arXX39df/rTn/T666+7a5OkgQMHavTo0brppps0fvx43XHHHe763nnnHa1evVozZ868ZE3PP/+8tm7dqmXLlik7O1tjxozxeB0cDocWLlyoH/3oR+rXr59mzJih119/nbN2AFDNVftgdyUSExO1fft2vfDCC/rkk0+0fft2hYSE6PTp0xW639WrV2vRokWqWbOmatasqd69e0uSQkNDNWnSJEnfX1KdN2+eTpw4of379ys3N1dNmzZVQECAGjRoIElq0KCBli5dquPHj+vrr7/Wl19+KX9/fzVr1sy9r/DwcLVt29Zj/23atHGH13OXEw8dOuTR59ChQ+62sLAwHT582KP97NmzOnr0qEefsrbxw32Eh4erZcuW8vHx8aglLy/vgq95ebY7YcIEDR06VL/85S/Vvn173X333Zo2bZpSU1Pdge3JJ5/U+PHjdf/996t9+/YaOnSoRo8e7Q7poaGhqlmz5kVfq9WrVys7O1vBwcHusZOke++9V7169SpVd+vWrXXnnXdq7ty5mjNnjg4ePOh+HRo3bqygoCCP/Rhj9O9//7vM1wEAUD1U+2AXGBioiIgIbdiwwWP5hg0b3P9Jn7s/6/x7uTZs2KDHHntM/fr104033iin06nvvvuuwmt+++239dlnn2n79u3avn27Xn31VUnSxx9/rKSkJI++tWrVUmRkpHx8fJSWlqY77rij1D15vr6+aty4sc6ePau3335bAwcOdLd17dpVWVlZHv337Nmj6OhoSVJMTIzCwsK0atUqd3thYaEyMzMVFxcnSYqLi1N+fr62bNni7rN69WqVlJQoNjbW3WfdunXu+/8k6aOPPlKrVq1Ur149dy379u1zh61ztYSHh3vcQ/dD5dnuiRMnSr0m58Kj+b/HmVyoz7laateurS5dulz0tRo/frx27NjhHrdzj895/vnn9dprr5VZv/T/ZwNdLpf7dThw4IDH42L27NmjGjVqKDIy8oLbAQBUAxV+Yfg81+I9ds8//7wJDAw0aWlp5ssvvzTjxo0ztWrVMnv27DHGGHPmzBnj5+dnnnnmGZOXl2fy8/ONMcZ06tTJ3HbbbWb37t1m48aNpnv37sbPz888//zz7m3rMu6x+/rrr822bdvMlClTjL+/v9m2bZvZtm3bBR/JUtY9dllZWebvf/+72bNnj8nMzDSDBg0y9evXNzk5Oe4+GzduNG+//bbJzs4269atMz/96U9NTEyMx3Y+/fRTU7NmTTN16lSzd+9es3DhQlOnTh3zj3/8w91n+vTpJjg42Cxbtszs2LHDDBw4sMzHnXTq1MlkZmaa9evXmxYtWng8wiM/P980atTIDB061Hz++ecmLS3N1KlTx+OxJLm5uSYgIMAkJyebrKws895775mGDRuaZ555xt3nr3/9q8c3XMuz3UmTJpmAgADz5ptvmq+++sr861//Ms2bN3ffI2iMMYmJiaZx48bux50sXrzYhIaGmt/+9rfuPosXLza1atUyL7/8stm7d6/561//anx8fMzHH398oaEu9f5Yvny5mTdvntm5c6fJyckx7733nmnTpo3HN3SLiopMZGSk+dnPfmZ27dpl1q5da1q0aGF++ctfXnA/AIDrF487uYSyHncyefJk07hxY1OrVq1SjzsxxphXXnnFREVFmRo1argfd7J161bTuXNn4+vra1q0aGEWLVpkoqOjLxrsevbsaRITEy9aX2JiYqlHrEgya9asKbN/WcFu9+7d5qabbjJ+fn4mMDDQDBw40Hz55Zce66Wnp5s2bdoYp9NpQkJCzNChQ92PH/mhd99917Rr1844nU7TunVr8/LLL3u0l5SUmAkTJphGjRoZp9NpevfubbKysjz6HDlyxAwePNj4+/ubwMBAM2zYsFJB9bPPPjPdunUzTqfTNG7c2EyfPr1ULZ988omJjY01TqfTNGvWzEydOtWcPXvW3T5p0iQTHR3t1XbPnDljJk+ebJo3b258fX1NVFSUGTVqlMfrWVhYaB5//HHTpEkT4+vra5o1a2aeeuopj8eoGGPM3/72N3PDDTcYX19f07FjR7N06dJSx/BD578/Vq9ebeLi4kxQUJD7fTVu3DiPWowx5osvvjDx8fHGz8/PREZGmjFjxpgTJ05cdF8AgOuTN9nJYcx5j86vYIWFhQoKClJBQYECAwMrc9cAAADXHW+yU7W/xw4AAMAWBDsAAABLEOwAAAAsQbADAACwRM2qLgAArhVNxy+v6hIqzf7p/au6BAAVoFqdsXvooYfkcDj0yCOPlGpLSkqSw+HQQw89VPmFlcPixYt1++23KyQkRA6Hw/1w2x/q1auXHA6Hx1TWsUrSkSNHFBkZKYfDofz8/Ivue+rUqbr11ltVp04dBQcHl7mtvn37KiIiQk6nU1FRUUpOTlZhYeFlHCkAALhc1SrYSVJUVJTS0tI8flPz1KlTeuONN9SkSZMqrOzijh8/rm7duumPf/zjRfuNGDFCBw8edE/PPvtsmf2GDx+uDh06lGvfp0+f1n333adHH320zPYaNWpo4MCBeuedd7Rnzx7Nnz9fK1euvGCoBAAAFaPaBbubb75ZUVFRWrx4sXvZ4sWL1aRJE3Xq1Mmjb0lJiVJTUxUTEyM/Pz917NhR//M//+NuLy4u1vDhw93trVq10gsvvOCxjYceekh33XWX/vSnPyk8PFwhISFKSkry+Imr8hg6dKgmTpyo+Pj4i/arU6eOwsLC3FNZz7uZM2eO8vPzNXbs2HLte8qUKRo9erTat29fZnu9evX06KOPqnPnzoqOjlbv3r01atQoffzxx+XaPgAAuDqqXbCTpIcfftjjtznnzZunYcOGleqXmpqqBQsW6KWXXtKuXbs0evRoPfDAA1q7dq2k74NfZGSkFi1apN27d2vixIn63e9+p7feestjO2vWrFF2drbWrFmj119/XfPnz9f8+fPd7ZMnT1bTpk2vyrEtXLhQoaGhateunVJSUnTixAmP9t27d+vpp5/WggULSv326dVy4MABLV68WD179qyQ7QMAgLJVyy9PPPDAA0pJSdHXX38tSdqwYYPS0tKUnp7u7uNyuTRt2jStXLnS/WP2zZo10/r16zV37lz17NlTtWrV0pQpU9zrxMTEKCMjQ2+99ZZ+/vOfu5fXq1dPL774onx8fNS6dWv1799fq1at0ogRIyRJoaGhat68+RUf1y9+8QtFR0crIiJCO3bs0Lhx45SVleU+O+lyuTR48GA999xzatKkib766qsr3ucPDR48WMuWLdPJkyc1YMAAvfrqq1d1+wAA4OKqZbBr0KCB+vfvr/nz58sYo/79+ys0NNSjz759+3TixAnddtttHstPnz7tccl21qxZmjdvnnJzc3Xy5EmdPn1aN910k8c6N954o3x8fNzz4eHh2rlzp3s+OTlZycnJV3xcI0eOdP+7ffv2Cg8PV+/evZWdna3mzZsrJSVFbdq00QMPPHDF+yrL888/r0mTJmnPnj1KSUnRmDFjNHv27ArZFwAAKK1aBjvp+8ux58LUrFmzSrUfO3ZMkrR8+XI1btzYo83pdEqS0tLSNHbsWP35z39WXFycAgIC9NxzzykzM9Ojf61atTzmHQ6HSkpKrtqxXEhsbKyk70Nq8+bNtXr1au3cudN9n+C5nwkODQ3VU0895XH28XKcu6+vdevWql+/vrp3764JEyYoPDz8yg4EAACUS7UNdn379tXp06flcDjUp0+fUu1t27aV0+lUbm7uBe8V27Bhg2699VaNGjXKvSw7O7vCavbWuUeinAtWb7/9tse3gTdt2qSHH35YH3/88VW5FPxD54Kry+W6qtsFAAAXVm2DnY+Pj7744gv3v88XEBCgsWPHavTo0SopKVG3bt1UUFCgDRs2KDAwUImJiWrRooUWLFigDz/8UDExMfr73/+uTZs2KSYmxqtaXnzxRS1ZskSrVq26YJ+jR48qNzdXBw4ckCRlZWVJ+v+zZNnZ2XrjjTfUr18/hYSEaMeOHRo9erR69OjhfqzJ+eHtu+++kyS1adPG/Xy6Tz/9VA8++KBWrVrlPlOZm5vr3n9xcbE7MN5www3y9/fX+++/r0OHDqlLly7y9/fXrl279OSTT6pr165X7UshAADg0qptsJNU5qNAfugPf/iDGjRooNTUVH311VcKDg7WzTffrN/97neSpF/96lfatm2bBg0aJIfDocGDB2vUqFH64IMPvKrju+++u+SZvnfeecfjm7v333+/JGnSpEmaPHmyateurZUrV2rmzJk6fvy4oqKidO+99+r3v/+9V7WcOHFCWVlZHo9jmThxol5//XX3/Ll7DNesWaNevXrJz89Pr7zyikaPHi2Xy6WoqCjdc889Gj9+vFf7BgAAV8Zhzt1oVUkKCwsVFBSkgoKCSwYrAKhM/KQYgGuRN9mpWj7HDgAAwEYEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAAS3gV7CZPniyHw+ExtW7duqJqAwAAgBdqervCjTfeqJUrV/7/Bmp6vQkAAABUAK9TWc2aNRUWFlYRtQAAAOAKeH2P3d69exUREaFmzZppyJAhys3NvWh/l8ulwsJCjwkAAABXn1fBLjY2VvPnz9eKFSs0Z84c5eTkqHv37ioqKrrgOqmpqQoKCnJPUVFRV1w0AAAASnMYY8zlrpyfn6/o6GjNmDFDw4cPL7OPy+WSy+VyzxcWFioqKkoFBQUKDAy83F0DwFXXdPzyqi6h0uyf3r+qSwBQToWFhQoKCipXdrqibz4EBwerZcuW2rdv3wX7OJ1OOZ3OK9kNAAAAyuGKnmN37NgxZWdnKzw8/GrVAwAAgMvkVbAbO3as1q5dq/379+uTTz7R3XffLR8fHw0ePLii6gMAAEA5eXUp9t///rcGDx6sI0eOqEGDBurWrZs2btyoBg0aVFR9AAAAKCevgl1aWlpF1QEAAIArxG/FAgAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJWpWdQEArm1Nxy+v6hIAAOXEGTsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsMQVBbvp06fL4XDoiSeeuErlAAAA4HJddrDbtGmT5s6dqw4dOlzNegAAAHCZLivYHTt2TEOGDNErr7yievXqXe2aAAAAcBkuK9glJSWpf//+io+Pv2Rfl8ulwsJCjwkAAABXX01vV0hLS9PWrVu1adOmcvVPTU3VlClTvC4MAAAA3vHqjN0333yjxx9/XAsXLpSvr2+51klJSVFBQYF7+uabby6rUAAAAFycV2fstmzZosOHD+vmm292LysuLta6dev04osvyuVyycfHx2Mdp9Mpp9N5daoFAADABXkV7Hr37q2dO3d6LBs2bJhat26tcePGlQp1AAAAqDxeBbuAgAC1a9fOY1ndunUVEhJSajkAAAAqF788AQAAYAmvvxV7vvT09KtQBgAAAK4UZ+wAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsIRXwW7OnDnq0KGDAgMDFRgYqLi4OH3wwQcVVRsAAAC84FWwi4yM1PTp07VlyxZt3rxZP/3pTzVw4EDt2rWrouoDAABAOdX0pvOAAQM85qdOnao5c+Zo48aNuvHGG69qYQAAAPCOV8Huh4qLi7Vo0SIdP35ccXFxV7MmAAAAXAavg93OnTsVFxenU6dOyd/fX0uWLFHbtm0v2N/lcsnlcrnnCwsLL69SAAAAXJTX34pt1aqVtm/frszMTD366KNKTEzU7t27L9g/NTVVQUFB7ikqKuqKCgYAAEDZHMYYcyUbiI+PV/PmzTV37twy28s6YxcVFaWCggIFBgZeya4BVIKm45dXdQmoAPun96/qEgCUU2FhoYKCgsqVnS77HrtzSkpKPILb+ZxOp5xO55XuBgAAAJfgVbBLSUlRQkKCmjRpoqKiIr3xxhtKT0/Xhx9+WFH1AQAAoJy8CnaHDx/Wgw8+qIMHDyooKEgdOnTQhx9+qNtuu62i6gMAAEA5eRXs/va3v1VUHQAAALhC/FYsAACAJQh2AAAAliDYAQAAWIJgBwAAYIkrfo4dAOD6U90ePM0DmVFdcMYOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAAS3gV7FJTU9WlSxcFBASoYcOGuuuuu5SVlVVRtQEAAMALXgW7tWvXKikpSRs3btRHH32kM2fO6Pbbb9fx48crqj4AAACUU01vOq9YscJjfv78+WrYsKG2bNmiHj16XNXCAAAA4J0ruseuoKBAklS/fv2rUgwAAAAun1dn7H6opKRETzzxhLp27ap27dpdsJ/L5ZLL5XLPFxYWXu4uAQAAcBGXfcYuKSlJn3/+udLS0i7aLzU1VUFBQe4pKirqcncJAACAi7isYJecnKz33ntPa9asUWRk5EX7pqSkqKCgwD198803l1UoAAAALs6rS7HGGP3617/WkiVLlJ6erpiYmEuu43Q65XQ6L7tAAAAAlI9XwS4pKUlvvPGGli1bpoCAAOXl5UmSgoKC5OfnVyEFAgAAoHy8uhQ7Z84cFRQUqFevXgoPD3dP//znPyuqPgAAAJST15diAQAAcG3it2IBAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAAS9Ss6gKA61HT8curugQAAErhjB0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACW8DrYrVu3TgMGDFBERIQcDoeWLl1aAWUBAADAW14Hu+PHj6tjx46aNWtWRdQDAACAy1TT2xUSEhKUkJBQEbUAAADgCnCPHQAAgCW8PmPnLZfLJZfL5Z4vLCys6F0CAABUSxV+xi41NVVBQUHuKSoqqqJ3CQAAUC1VeLBLSUlRQUGBe/rmm28qepcAAADVUoVfinU6nXI6nRW9GwAAgGrP62B37Ngx7du3zz2fk5Oj7du3q379+mrSpMlVLQ4AAADl53Ww27x5s37yk5+458eMGSNJSkxM1Pz5869aYQAAAPCO18GuV69eMsZURC0AAAC4AjzHDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEsQ7AAAACxBsAMAALAEwQ4AAMASBDsAAABL1KzqAgAAqGhNxy+v6hIqzf7p/au6BFQhztgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAliDYAQAAWIJgBwAAYAmCHQAAgCUIdgAAAJYg2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJWpWdQGwR9Pxy6u6BAAAqjXO2AEAAFiCYAcAAGAJgh0AAIAlCHYAAACWINgBAABYgmAHAABgCYIdAACAJQh2AAAAlrisYDdr1iw1bdpUvr6+io2N1aeffnq16wIAAICXvP7liX/+858aM2aMXnrpJcXGxmrmzJnq06ePsrKy1LBhw4qoEQAAlFN1+hWg/dP7V3UJ1xyvz9jNmDFDI0aM0LBhw9S2bVu99NJLqlOnjubNm1cR9QEAAKCcvAp2p0+f1pYtWxQfH///G6hRQ/Hx8crIyLjqxQEAAKD8vLoU+91336m4uFiNGjXyWN6oUSN9+eWXZa7jcrnkcrnc8wUFBZKkwsJCb2u9LrWb9GFVlwAAgJWqS5Y4d5zGmEv29foeO2+lpqZqypQppZZHRUVV9K4BAIDFgmZWdQWVq6ioSEFBQRft41WwCw0NlY+Pjw4dOuSx/NChQwoLCytznZSUFI0ZM8Y9X1JSoqNHjyokJEQOh8Ob3V9zCgsLFRUVpW+++UaBgYFVXQ4qAWNevTDe1QvjXf1cL2NujFFRUZEiIiIu2derYFe7dm3dcsstWrVqle666y5J3we1VatWKTk5ucx1nE6nnE6nx7Lg4GBvdnvNCwwMvKbfELj6GPPqhfGuXhjv6ud6GPNLnak7x+tLsWPGjFFiYqI6d+6sH/3oR5o5c6aOHz+uYcOGeV0kAAAArh6vg92gQYP0n//8RxMnTlReXp5uuukmrVixotQXKgAAAFC5LuvLE8nJyRe89FqdOJ1OTZo0qdSlZtiLMa9eGO/qhfGufmwcc4cpz3dnAQAAcM27rN+KBQAAwLWHYAcAAGAJgh0AAIAlCHaXMHnyZDkcDo+pdevW7vZTp04pKSlJISEh8vf317333lvqAc64tq1bt04DBgxQRESEHA6Hli5d6tFujNHEiRMVHh4uPz8/xcfHa+/evR59jh49qiFDhigwMFDBwcEaPny4jh07VolHgfK61Hg/9NBDpf7m+/bt69GH8b5+pKamqkuXLgoICFDDhg111113KSsry6NPeT7Hc3Nz1b9/f9WpU0cNGzbUk08+qbNnz1bmoaCcyjPmvXr1KvV3/sgjj3j0uV7HnGBXDjfeeKMOHjzontavX+9uGz16tN59910tWrRIa9eu1YEDB3TPPfdUYbXw1vHjx9WxY0fNmjWrzPZnn31Wf/nLX/TSSy8pMzNTdevWVZ8+fXTq1Cl3nyFDhmjXrl366KOP9N5772ndunUaOXJkZR0CvHCp8Zakvn37evzNv/nmmx7tjPf1Y+3atUpKStLGjRv10Ucf6cyZM7r99tt1/Phxd59LfY4XFxerf//+On36tD755BO9/vrrmj9/viZOnFgVh4RLKM+YS9KIESM8/s6fffZZd9t1PeYGFzVp0iTTsWPHMtvy8/NNrVq1zKJFi9zLvvjiCyPJZGRkVFKFuJokmSVLlrjnS0pKTFhYmHnuuefcy/Lz843T6TRvvvmmMcaY3bt3G0lm06ZN7j4ffPCBcTgc5ttvv6202uG988fbGGMSExPNwIEDL7gO4319O3z4sJFk1q5da4wp3+f4+++/b2rUqGHy8vLcfebMmWMCAwONy+Wq3AOA184fc2OM6dmzp3n88ccvuM71POacsSuHvXv3KiIiQs2aNdOQIUOUm5srSdqyZYvOnDmj+Ph4d9/WrVurSZMmysjIqKpycRXl5OQoLy/PY4yDgoIUGxvrHuOMjAwFBwerc+fO7j7x8fGqUaOGMjMzK71mXLn09HQ1bNhQrVq10qOPPqojR4642xjv61tBQYEkqX79+pLK9zmekZGh9u3bezyIv0+fPiosLNSuXbsqsXpcjvPH/JyFCxcqNDRU7dq1U0pKik6cOOFuu57H/LIeUFydxMbGav78+WrVqpUOHjyoKVOmqHv37vr888+Vl5en2rVrl/rt20aNGikvL69qCsZVdW4cz/9llR+OcV5enho2bOjRXrNmTdWvX5/3wXWob9++uueeexQTE6Ps7Gz97ne/U0JCgjIyMuTj48N4X8dKSkr0xBNPqGvXrmrXrp0kletzPC8vr8zPgHNtuHaVNeaS9Itf/ELR0dGKiIjQjh07NG7cOGVlZWnx4sWSru8xJ9hdQkJCgvvfHTp0UGxsrKKjo/XWW2/Jz8+vCisDUBHuv/9+97/bt2+vDh06qHnz5kpPT1fv3r2rsDJcqaSkJH3++ece90nDbhca8x/eE9u+fXuFh4erd+/eys7OVvPmzSu7zKuKS7FeCg4OVsuWLbVv3z6FhYXp9OnTys/P9+hz6NAhhYWFVU2BuKrOjeP535D74RiHhYXp8OHDHu1nz57V0aNHeR9YoFmzZgoNDdW+ffskMd7Xq+TkZL333ntas2aNIiMj3cvL8zkeFhZW5mfAuTZcmy405mWJjY2VJI+/8+t1zAl2Xjp27Jiys7MVHh6uW265RbVq1dKqVavc7VlZWcrNzVVcXFwVVomrJSYmRmFhYR5jXFhYqMzMTPcYx8XFKT8/X1u2bHH3Wb16tUpKStwfFrh+/fvf/9aRI0cUHh4uifG+3hhjlJycrCVLlmj16tWKiYnxaC/P53hcXJx27tzpEeg/+ugjBQYGqm3btpVzICi3S415WbZv3y5JHn/n1+2YV/W3N651v/nNb0x6errJyckxGzZsMPHx8SY0NNQcPnzYGGPMI488Ypo0aWJWr15tNm/ebOLi4kxcXFwVVw1vFBUVmW3btplt27YZSWbGjBlm27Zt5uuvvzbGGDN9+nQTHBxsli1bZnbs2GEGDhxoYmJizMmTJ93b6Nu3r+nUqZPJzMw069evNy1atDCDBw+uqkPCRVxsvIuKiszYsWNNRkaGycnJMStXrjQ333yzadGihTl16pR7G4z39ePRRx81QUFBJj093Rw8eNA9nThxwt3nUp/jZ8+eNe3atTO333672b59u1mxYoVp0KCBSUlJqYpDwiVcasz37dtnnn76abN582aTk5Njli1bZpo1a2Z69Ojh3sb1POYEu0sYNGiQCQ8PN7Vr1zaNGzc2gwYNMvv27XO3nzx50owaNcrUq1fP1KlTx9x9993m4MGDVVgxvLVmzRojqdSUmJhojPn+kScTJkwwjRo1Mk6n0/Tu3dtkZWV5bOPIkSNm8ODBxt/f3wQGBpphw4aZoqKiKjgaXMrFxvvEiRPm9ttvNw0aNDC1atUy0dHRZsSIER6PPDCG8b6elDXWksxrr73m7lOez/H9+/ebhIQE4+fnZ0JDQ81vfvMbc+bMmUo+GpTHpcY8NzfX9OjRw9SvX984nU5zww03mCeffNIUFBR4bOd6HXOHMcZU3vlBAAAAVBTusQMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwA4D/07RpU82cObOqywCAy0awA4D/s2nTJo0cObJcfY0xmjhxosLDw+Xn56f4+Hjt3bu3gisEgIsj2AHA/2nQoIHq1KlTrr7PPvus/vKXv+ill15SZmam6tatqz59+ujUqVMVXCUAXBjBDoA1VqxYoW7duik4OFghISG64447lJ2dLUlasGCB/P39Pc6qjRo1Sq1bt9aJEyckeV6KNcZo8uTJatKkiZxOpyIiIvTYY4+522bOnKnf//73GjhwoDp06KAFCxbowIEDWrp0aaUeMwD8EMEOgDWOHz+uMWPGaPPmzVq1apVq1Kihu+++WyUlJXrwwQfVr18/DRkyRGfPntXy5cv16quvauHChWWepXv77bf1/PPPa+7cudq7d6+WLl2q9u3bS5JycnKUl5en+Ph4d/+goCDFxsYqIyOj0o4XAM5Xs6oLAICr5d577/WYnzdvnho0aKDdu3erXbt2mjt3rjp06KDHHntMixcv1uTJk3XLLbeUua3c3FyFhYUpPj5etWrVUpMmTfSjH/1IkpSXlydJatSokcc6jRo1crcBQFXgjB0Aa+zdu1eDBw9Ws2bNFBgYqKZNm0r6PqRJUr169fS3v/1Nc+bMUfPmzTV+/PgLbuu+++7TyZMn1axZM40YMUJLlizR2bNnK+MwAOCyEewAWGPAgAE6evSoXnnlFWVmZiozM1OSdPr0aXefdevWycfHRwcPHtTx48cvuK2oqChlZWVp9uzZ8vPz06hRo9SjRw+dOXNGYWFhkqRDhw55rHPo0CF3GwBUBYIdACscOXJEWVlZ+v3vf6/evXurTZs2+u9//+vR55NPPtEf//hHvfvuu/L391dycvJFt+nn56cBAwboL3/5i9LT05WRkaGdO3cqJiZGYWFhWrVqlbtvYWGhMjMzFRcXVyHHBwDlwT12AKxQr149hYSE6OWXX1Z4eLhyc3M9LrUWFRVp6NCheuyxx5SQkKDIyEh16dJFAwYM0M9+9rNS25s/f76Ki4sVGxurOnXq6B//+If8/PwUHR0th8OhJ554Qs8884xatGihmJgYTZgwQREREbrrrrsq8agBwBNn7ABYoUaNGkpLS9OWLVvUrl07jR49Ws8995y7/fHHH1fdunU1bdo0SVL79u01bdo0/epXv9K3335banvBwcF65ZVX1LVrV3Xo0EErV67Uu+++q5CQEEnSb3/7W/3617/WyJEj1aVLFx07dkwrVqyQr69v5RwwAJTBYYwxVV0EAAAArhxn7AAAACxBsAMAALAEwQ4AAMASBDsAAABLEOwAAAAsQbADAACwBMEOAADAEgQ7AAAASxDsAAAALEGwAwAAsATBDgAAwBIEOwAAAEv8L2cfHGn4c/vaAAAAAElFTkSuQmCC\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h = physt.h1(data, weights=data)\n", "h.plot(show_stats=[\"total\", \"mean\"]);" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": "
", "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAABkF0lEQVR4nO3deVhO6f8H8PdTadcirbTZKqQSmrKPVDT2GVu2bIMalDFkTJaZL2MbmWHsirEbZBuRUgYJ0djSkCVLxUhFqNT5/eHXGY+WKUvp9H5d13NdnnM+5z73ueeZ4+2sMkEQBBARERFRladQ2R0gIiIioveDwY6IiIhIIhjsiIiIiCSCwY6IiIhIIhjsiIiIiCSCwY6IiIhIIhjsiIiIiCSCwY6IiIhIIpQquwOVqaCgAPfv30fNmjUhk8kquztEREQkAYIg4MmTJzAxMYGCQsUeQ6vWwe7+/fswNTWt7G4QERGRBN25cwd169at0HVW62BXs2ZNAK8GXktLq5J7Q0RERFKQlZUFU1NTMWdUpGod7ApPv2ppaTHYERER0XtVGZd58eYJIiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSCAY7IiIiIolgsCMiIiKSiHIFu+XLl6NZs2bQ0tKClpYWnJ2dcfDgQXH+ixcv4OPjAz09PWhqaqJPnz5IS0uTayM5ORmenp5QV1eHgYEBJk+ejJcvX8rVREVFoXnz5lBRUUGDBg0QEhJSpC/Lli2DhYUFVFVV4eTkhNOnT5dnU4iIiIgkp1zBrm7duvjxxx8RFxeHs2fP4tNPP0WPHj1w+fJlAICfnx/27duHHTt2IDo6Gvfv30fv3r3F5fPz8+Hp6Ync3FycPHkS69evR0hICAIDA8WamzdvwtPTEx07dkR8fDwmTpyIkSNH4tChQ2LNtm3b4O/vjxkzZuDcuXOws7ODu7s7Hjx48K7jITkymQyhoaGV3Q0iIiKqCMI70tXVFdasWSNkZGQINWrUEHbs2CHOS0hIEAAIMTExgiAIwh9//CEoKCgIqampYs3y5csFLS0tIScnRxAEQfjmm2+EJk2ayK2jX79+gru7u/i9VatWgo+Pj/g9Pz9fMDExEebOnVuuvmdmZgoAhMzMzHIt9zYAlPqZMWNGicvevHlTACCcP3/+rda7e/fuMtffvHlTGD58uGBhYSGoqqoK9erVEwIDA8X/PoIgCEePHhW6d+8uGBkZCerq6oKdnZ2wceNGuXZyc3OFWbNmCfXq1RNUVFSEZs2aCQcPHpSrycrKEiZMmCCYmZkJqqqqgrOzs3D69Oki/S/uM3/+fLEvJdW82RYREVFFqMh88aa3vsYuPz8fW7duRXZ2NpydnREXF4e8vDy4urqKNdbW1jAzM0NMTAwAICYmBra2tjA0NBRr3N3dkZWVJR71i4mJkWujsKawjdzcXMTFxcnVKCgowNXVVawpSU5ODrKysuQ+FSUlJUX8BAUFQUtLS27a119/XWF9Kc3Vq1dRUFCAlStX4vLly1i8eDFWrFiBadOmiTUnT55Es2bNsHPnTly4cAHe3t4YMmQI9u/fL9ZMnz4dK1euxC+//IIrV65gzJgx6NWrF86fPy/WjBw5EuHh4fjtt99w8eJFuLm5wdXVFffu3RNrXh+jlJQUrFu3DjKZDH369AEAuLi4FKkZOXIkLC0t0aJFiwoYMSIioo9IeZPghQsXBA0NDUFRUVHQ1tYWDhw4IAiCIGzatElQVlYuUt+yZUvhm2++EQRBEEaNGiW4ubnJzc/OzhYACH/88YcgCILQsGFDYc6cOXI1Bw4cEAAIz549E+7duycAEE6ePClXM3nyZKFVq1al9n3GjBnFHtmp6EQdHBwsaGtri9/z8/OFWbNmCXXq1BGUlZUFOzs7uaNbb/a3ffv2giAIwunTpwVXV1dBT09P0NLSEtq1ayfExcXJrQvlPGJXnPnz5wuWlpal1nTt2lXw9vYWvxsbGwtLly6Vq+ndu7fg5eUlCIIgPHv2TFBUVBT2798vV9O8eXPh22+/LXE9PXr0ED799NMS5+fm5gr6+vrC7NmzS+0vERHRh1KljthZWVkhPj4esbGxGDt2LIYOHYorV668n5T5gQUEBCAzM1P83Llzp7K7BABYsmQJFi1ahIULF+LChQtwd3dH9+7dce3aNQAQbww5cuQIUlJSsGvXLgDAkydPMHToUBw/fhynTp1Cw4YN0bVrVzx58qTEdXXo0AHDhg0rV/8yMzNRq1atctXk5ORAVVVVrkZNTQ3Hjx8HALx8+RL5+fml1rwpLS0NBw4cwIgRI0rsx969e/Ho0SN4e3uX2l8iIiIpKnewU1ZWRoMGDeDo6Ii5c+fCzs4OS5YsgZGREXJzc5GRkSFXn5aWBiMjIwCAkZFRkbtkC7//V42WlhbU1NRQu3ZtKCoqFltT2EZJVFRUxDt6Cz8fg4ULF2LKlCno378/rKysMG/ePNjb2yMoKAgAoK+vDwDQ09ODkZGRGKA+/fRTDBo0CNbW1rCxscGqVavw7NkzREdHl7guMzMzGBsbl7lv169fxy+//IIvv/yyxJrt27fjzJkzcmHK3d0dP/30E65du4aCggKEh4dj165dSElJAQDUrFkTzs7O+P7773H//n3k5+dj48aNiImJEWvetH79etSsWVPuhpw3rV27Fu7u7qhbt26Zt5GIiEgq3vk5dgUFBcjJyYGjoyNq1KiBiIgIcV5iYiKSk5Ph7OwMAHB2dsbFixfl7l4NDw+HlpYWGjduLNa83kZhTWEbysrKcHR0lKspKChARESEWFOVZGVl4f79+2jdurXc9NatWyMhIaHUZdPS0jBq1Cg0bNgQ2tra0NLSwtOnT5GcnFziMhs2bMDcuXPL1Ld79+7Bw8MDX3zxBUaNGlVszdGjR+Ht7Y3Vq1ejSZMm4vQlS5agYcOGsLa2hrKyMnx9feHt7Q0FhX9/cr/99hsEQUCdOnWgoqKCn3/+GQMGDJCred26devg5eVV5Chfobt37+LQoUOlHtEjIiKSMqXyFAcEBKBLly4wMzPDkydPsHnzZkRFReHQoUPQ1tbGiBEj4O/vj1q1akFLSwtfffUVnJ2d8cknnwAA3Nzc0LhxYwwePBjz589Hamoqpk+fDh8fH6ioqAAAxowZg6VLl+Kbb77B8OHDERkZie3bt+PAgQNiP/z9/TF06FC0aNECrVq1QlBQELKzs6vd6behQ4fi0aNHWLJkCczNzaGiogJnZ2fk5ua+c9v3799Hx44d4eLiglWrVhVbEx0djW7dumHx4sUYMmSI3Dx9fX2EhobixYsXePToEUxMTDB16lTUq1dPrKlfvz6io6ORnZ2NrKwsGBsbo1+/fnI1hf78808kJiZi27ZtJfY5ODgYenp66N69+1tuNRERURVXngvyhg8fLpibmwvKysqCvr6+0KlTJ+Hw4cPi/OfPnwvjxo0TdHV1BXV1daFXr15CSkqKXBu3bt0SunTpIqipqQm1a9cWJk2aJOTl5cnVHD16VLC3txeUlZWFevXqCcHBwUX68ssvvwhmZmaCsrKy0KpVK+HUqVPl2RRBECrv4sY3b54wMTER/ve//8nVtGzZUnykS+ENI2fPnpWr0dTUFDZs2CB+T05OFgAIixcvFqfhLW6euHv3rtCwYUOhf//+wsuXL4utOXr0qKChoVHkBomS5ObmCvXr1xcCAgJKrElPTxe0tbWFlStXFpk3dOhQwdHRscRlCwoKBEtLS2HSpEll6g8REdGHUpk3T7zzc+yqso8l2C1evFjQ0tIStm7dKly9elWYMmWKUKNGDeHvv/8WBEEQ8vLyBDU1NeGHH34QUlNThYyMDEEQBMHBwUHo3LmzcOXKFeHUqVNC27ZtBTU1tVKD3eDBg4WpU6eW2Le7d+8KDRo0EDp16iTcvXtXSElJET+FIiMjBXV1dSEgIEBu/qNHj8SaU6dOCTt37hSSkpKEY8eOCZ9++qlgaWkpPH78WKwJCwsTDh48KNy4cUM4fPiwYGdnJzg5OQm5ublyfcrMzBTU1dWF5cuXl9jvI0eOCACEhISEEmuIiIgqAoNdJflYgl1+fr4wc+ZMoU6dOkKNGjWKPO5EEARh9erVgqmpqaCgoCA+7uTcuXNCixYtBFVVVaFhw4bCjh07BHNz81KDXfv27YWhQ4eW2jeU8MDfQkOHDi12fmG/BEEQoqKiBBsbG0FFRUXQ09MTBg8eLNy7d09uXdu2bRPq1asnKCsrC0ZGRoKPj48YWl+3cuVKQU1Nrdh5hQYMGCC4uLiUOJ+IiKiiVGawkwmCIFTgmd+PSlZWFrS1tZGZmfnR3CFLREREVVtl5ot3viuWiIiIiD4ODHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQR5XrzBH3cLCwsMHHiREycOPGd2hEEAV9++SV+//13PH78GOfPn8fEiRPl3l9LRET0vlhMPfDfRVVIQc6zSls3j9hVoA4dOhQbukJCQqCjo1Ph/SlJWFgYQkJCsH//fqSkpKBp06aV3SUcO3YM3bp1g4mJCWQyGUJDQ4vU7Nq1C25ubtDT04NMJkN8fHyRmg4dOkAmk8l9xowZI85/9OgRPDw8YGJiAhUVFZiamsLX1xdZWVnF9uvEiRNQUlKCvb19kXn37t3DoEGDoKenBzU1Ndja2uLs2bNvOwRERET/icGORIXvmE1KSoKxsTFcXFxgZGQEJaXKP7CbnZ0NOzs7LFu2rNSaNm3aYN68eaW2NWrUKKSkpIif+fPni/MUFBTQo0cP7N27F3///TdCQkJw5MgRufBXKCMjA0OGDEGnTp2KzHv8+DFat26NGjVq4ODBg7hy5QoWLVoEXV3dcmw1ERFR+TDYfYSGDRuGnj17YuHChTA2Noaenh58fHyQl5cn1jx48ADdunWDmpoaLC0tsWnTpiLtZGRkYOTIkdDX14eWlhY+/fRT/PXXX+L8mTNnwt7eHmvWrIGlpSVUVVUxbNgwfPXVV0hOToZMJoOFhUWxfXz8+DGGDBkCXV1dqKuro0uXLrh27RqAV6dy9fX18fvvv4v19vb2MDY2Fr8fP34cKioqePasbIeru3Tpgh9++AG9evUqsWbw4MEIDAyEq6trqW2pq6vDyMhI/Lz+8EhdXV2MHTsWLVq0gLm5OTp16oRx48bhzz//LNLOmDFjMHDgQDg7OxeZN2/ePJiamiI4OBitWrWCpaUl3NzcUL9+/TJtLxER0dtgsPtIHT16FElJSTh69CjWr1+PkJAQhISEiPOHDRuGO3fu4OjRo/j999/x66+/4sGDB3JtfPHFF3jw4AEOHjyIuLg4NG/eHJ06dUJ6erpYc/36dezcuRO7du1CfHw8lixZgtmzZ6Nu3bpISUnBmTNniu3fsGHDcPbsWezduxcxMTEQBAFdu3ZFXl4eZDIZ2rVrh6ioKACvQmBCQgKeP3+Oq1evAgCio6PRsmVLqKurIyQkBDKZ7P0OYCk2bdqE2rVro2nTpggICCg1XN6/fx+7du1C+/bt5aYHBwfjxo0bmDFjRrHL7d27Fy1atMAXX3wBAwMDODg4YPXq1e91O4iIiN5U+efYqFi6urpYunQpFBUVYW1tDU9PT0RERGDUqFH4+++/cfDgQZw+fRotW7YEAKxduxY2Njbi8sePH8fp06fx4MEDqKioAAAWLlyI0NBQ/P777xg9ejSAV6dfN2zYAH19fXHZmjVrQlFREUZGRsX27dq1a9i7dy9OnDgBFxcXAK/CkqmpKUJDQ/HFF1+gQ4cOWLlyJYBX18c5ODjAyMgIUVFRsLa2RlRUlBiWtLW1YWVl9Z5HsHgDBw6Eubk5TExMcOHCBUyZMgWJiYnYtWuXXN2AAQOwZ88ePH/+HN26dcOaNWvktn/q1Kn4888/SzxNfePGDSxfvhz+/v6YNm0azpw5g/Hjx0NZWRlDhw79oNtIRETVF4/YfaSaNGkCRUVF8buxsbF4RC4hIQFKSkpwdHQU51tbW8vdgPHXX3/h6dOn0NPTg6ampvi5efMmkpKSxDpzc3O5UFcWhet3cnISp+np6cHKygoJCQkAgPbt2+PKlSt4+PAhoqOj0aFDB3To0AFRUVHIy8vDyZMn0aFDBwBAr169xCN5H9ro0aPh7u4OW1tbeHl5YcOGDdi9e7fcmADA4sWLce7cOezZswdJSUnw9/cHAOTn52PgwIGYNWsWGjVqVOJ6CgoK0Lx5c8yZMwcODg4YPXo0Ro0ahRUrVnzQ7SMiouqNR+wqkJaWFjIzM4tMz8jIgLa2tty0GjVqyH2XyWQoKCgo87qePn0KY2Nj8XTo614PgBoaGmVuszxsbW1Rq1YtREdHIzo6Gv/73/9gZGSEefPm4cyZM8jLyxOP9lWmwnB6/fp1uevfCq+/s7a2Rq1atdC2bVt89913UFNTw9mzZ3H+/Hn4+voCeBXiBEGAkpISDh8+jE8//RTGxsZo3Lix3LpsbGywc+fOits4IiKqdhjsKpCVlRUOHz5cZPq5c+dKPfrzJmtra7x8+RJxcXHiqdjExERkZGSINc2bN0dqaiqUlJRKvAHibdnY2ODly5eIjY0Vw9mjR4+QmJgohhmZTIa2bdtiz549uHz5Mtq0aQN1dXXk5ORg5cqVaNGixQcLleVR+EiU12/seFNhoM7JyYGhoSEuXrwoN//XX39FZGQkfv/9d1haWgIAWrdujcTERLm6v//+G+bm5u+x90RERPIY7CrQ2LFjsXTpUowfPx4jR46EiooKDhw4gC1btmDfvn1lbsfKygoeHh748ssvsXz5cigpKWHixIlQU1MTa1xdXeHs7IyePXti/vz5aNSoEe7fv48DBw6gV69eaNGixVtvR8OGDdGjRw+MGjUKK1euRM2aNTF16lTUqVMHPXr0EOs6dOiASZMmoUWLFtDU1AQAtGvXDps2bcLkyZPFut27dyMgIKDU07FPnz7F9evXxe83b95EfHw8atWqBTMzMwBAeno6kpOTcf/+fQAQg1Xh0bekpCRs3rwZXbt2hZ6eHi5cuAA/Pz+0a9cOzZo1AwD88ccfSEtLQ8uWLaGpqYnLly9j8uTJaN26tRiQ33yun4GBAVRVVeWm+/n5wcXFBXPmzEHfvn1x+vRprFq1CqtWrSr3eBMREZUVr7GrQPXq1cOxY8dw9epVuLq6wsnJCdu3b8eOHTvg4eFRrraCg4NhYmKC9u3bo3fv3hg9ejQMDAzE+TKZDH/88QfatWsHb29vNGrUCP3798ft27dhaGj4ztsSHBwMR0dHfPbZZ3B2doYgCPjjjz/kTiG3b98e+fn54rV0wKuw9+a0zMzMIke33nT27Fk4ODjAwcEBAODv7w8HBwcEBgaKNXv37oWDgwM8PT0BAP3794eDg4N4XZuysjKOHDkCNzc3WFtbY9KkSejTp49cqFZTU8Pq1avRpk0b2NjYwM/PD927d8f+/fvLNT4tW7bE7t27sWXLFjRt2hTff/89goKC4OXlVa52iIiIykMmCIJQ2Z2oLFlZWdDW1kZmZqbcs8yIiIio4kjxlWJ3gvpWSr7gETsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBrtqKDQ0FA0aNICioiImTpyIkJAQuYcWExERUdXE59hVsDt37mDGjBkICwvDP//8A2NjY/Ts2ROBgYHQ09OrkD58+eWX8Pb2xvjx41GzZs2P4m0IgiBg0aJFWLVqFW7fvo3atWtj3Lhx+PbbbwEAUVFR6NixY5HlUlJSin2n7Y8//oiAgABMmDABQUFBcvNiYmLw7bffIjY2FoqKirC3t8ehQ4egpqZW4noAiO/mnTlzJmbNmlVkvrq6OrKzs99i64mIykZqd4/S+8dgV4Fu3LgBZ2dnNGrUCFu2bIGlpaX4ANyDBw/i1KlTqFWr1gdbf15eHnJycvDgwQO4u7vDxMTkg62rvCZMmIDDhw9j4cKFsLW1RXp6OtLT04vUJSYmyt06/vqz+wqdOXMGK1euFB86/LqYmBh4eHggICAAv/zyC5SUlPDXX39BQeHVwWsXFxekpKTILfPdd98hIiJCfKjz119/jTFjxsjVdOrUSXwLCBERUWXhqdgK5OPjA2VlZRw+fBjt27eHmZkZunTpgiNHjuDevXvi0alp06aJ7zB9nZ2dHWbPni1+X7NmDWxsbKCqqgpra2v8+uuv4rxbt25BJpNh27ZtaN++PVRVVbFp0ybUrFkTAPDpp59CJpMV+y5ZAFi+fDnq168PZWVlWFlZ4bfffhPnff311/jss8/E70FBQZDJZAgLCxOnNWjQAGvWrCnTuCQkJGD58uXYs2cPunfvDktLSzg6OqJz585Fag0MDMQ3SRgZGYmBrNDTp0/h5eWF1atXQ1dXt8jyfn5+GD9+PKZOnYomTZrAysoKffv2hYqKCoBXDzF+vX09PT3s2bMH3t7ekMlkAABNTU25mrS0NFy5cgUjRowo0/YSERF9KAx2FSQ9PR2HDh3CuHHj5F79Bbx65ZWXlxe2bdsGQRDg5eWF06dPIykpSay5fPkyLly4gIEDBwIANm3ahMDAQPzvf/9DQkIC5syZg++++w7r16+Xa3vq1KmYMGECEhIS0LFjR/ENDzt37kRKSor4rtfX7d69GxMmTMCkSZNw6dIl8dTt0aNHAbx6o8Tx48eRn58PAIiOjkbt2rXFkHjv3j0kJSWJb5eQyWQICQkpcWz27duHevXqYf/+/bC0tISFhQVGjhxZ7BE7e3t7GBsbo3Pnzjhx4kSR+T4+PvD09ISrq2uReQ8ePEBsbCwMDAzg4uICQ0NDcVtKsnfvXjx69Aje3t4l1qxZswaNGjVC27ZtS6whIiKqCAx2FeTatWsQBAE2NjbFzrexscHjx4/x8OFDNGnSBHZ2dti8ebM4f9OmTXByckKDBg0AADNmzMCiRYvQu3dvWFpaonfv3vDz88PKlSvl2p04caJYY25uLp66rFWrFoyMjKCsrFykLwsXLsSwYcMwbtw4NGrUCP7+/ujduzcWLlwIAGjbti2ePHmC8+fPQxAEHDt2DJMmTRKDXVRUFOrUqSP21crKCtra2iWOzY0bN3D79m3s2LEDGzZsQEhICOLi4vD555+LNcbGxlixYgV27tyJnTt3wtTUFB06dMC5c+fEmq1bt+LcuXOYO3duiesBgJkzZ2LUqFEICwtD8+bN0alTJ1y7dq3YZdauXQt3d3fUrVu32PkvXrzApk2beLSOiIg+CrzGroKV9Q1uXl5eWLduHb777jsIgoAtW7bA398fAJCdnY2kpCSMGDECo0aNEpd5+fJlkQBVeF1YeSQkJGD06NFy01q3bo0lS5YAAHR0dGBnZ4eoqCgoKytDWVkZo0ePxowZM/D06VNER0ejffv24rJXr14tdX0FBQXIycnBhg0b0KhRIwCvApWjoyMSExNhZWUlfgq5uLggKSkJixcvxm+//YY7d+5gwoQJCA8Ph6qqaonrAf69eQQAHBwcEBERgXXr1hUJhHfv3sWhQ4ewffv2Evu+e/duPHnyBEOHDi11G4mIiCoCg10FadCgAWQyGRISEtCrV68i8xMSEqCrqwt9fX0AwIABAzBlyhScO3cOz58/x507d9CvXz8Ar64jA4DVq1cXuRZPUVFR7ruGhsaH2Bx06NABUVFRUFFRQfv27VGrVi3Y2Njg+PHjiI6OxqRJk8rclrGxMZSUlMRQB0A8spmcnCwX6F7XqlUr8TRqXFwcHjx4gObNm4vz8/PzcezYMSxduhQ5OTkwNjYGADRu3FiuHRsbGyQnJxdpPzg4GHp6eujevXuJfV+zZg0+++wzGBoalnFriYiIPhyeiq0genp66Ny5M3799Vc8f/5cbl5qaio2bdqEfv36iRfo161bF+3bt8emTZuwadMmdO7cWTyNamhoCBMTE9y4cQMNGjSQ+1haWr5zX21sbIpcv3bixAm5QFR4bVpERIR4LV2HDh2wZcsW/P333+K0smjdujVevnwpd03h33//DQAwNzcvcbn4+HgxrHXq1AkXL15EfHy8+GnRogW8vLwQHx8PRUVFWFhYwMTERLzO8PV1vbkeQRAQHByMIUOGoEaNGsWu/+bNmzh69ChPwxIR0UeDR+wq0NKlS+Hi4gJ3d3f88MMPco87qVOnDv73v//J1Xt5eWHGjBnIzc3F4sWL5ebNmjUL48ePh7a2Njw8PJCTk4OzZ8/i8ePH4inbtzV58mT07dsXDg4OcHV1xb59+7Br1y4cOXJErGnXrh2ePHmC/fv348cffwTwKth9/vnnMDY2ljv6Zm1tjblz5xZ7pBIAXF1d0bx5cwwfPhxBQUEoKCiAj48POnfuLLYTFBQES0tLNGnSBC9evMCaNWsQGRmJw4cPAwBq1qyJpk2byrWroaEBPT09cbpMJsPkyZMxY8YM2NnZwd7eHuvXr8fVq1fx+++/yy0bGRmJmzdvYuTIkSWO07p162BsbIwuXbqUdWiJiIg+KAa7CtSwYUOcPXsWM2bMQN++fZGeng4jIyP07NkTM2bMKPIMu88//xy+vr5QVFREz5495eaNHDkS6urqWLBgASZPngwNDQ3Y2tpi4sSJ79zPnj17YsmSJVi4cCEmTJgAS0tLBAcHyx2F09XVha2tLdLS0mBtbQ3gVdgrKCiQu74OePXsuczMzBLXp6CggH379uGrr75Cu3btoKGhgS5dumDRokViTW5uLiZNmoR79+5BXV0dzZo1w5EjR0p8mHBJJk6ciBcvXsDPzw/p6emws7NDeHg46tevL1e3du1auLi4iNv2poKCAoSEhGDYsGFFTn8TERFVFplQ1qv5JSgrKwva2trIzMyUe+gtERHRx4hvnqgaCnKe4U5Q30rJFzxiR0REksMARNUVb54gIiIikggGOyIiIiKJYLAjIiIikggGOyIiIiKJYLCrQMOGDYNMJsOYMWOKzPPx8YFMJsOwYcMqvmNlsGvXLri5uUFPTw8ymQzx8fEl1gqCgC5dukAmkyE0NFRuXkREBFxcXFCzZk0YGRlhypQpePnyZanrTkpKQq9evaCvrw8tLS307dsXaWlpReoOHDgAJycnqKmpQVdXt8gjYoiIiKSOwa6CmZqaYuvWrXJvn3jx4gU2b94MMzOzSuxZ6bKzs9GmTRvMmzfvP2uDgoLEN2i87q+//kLXrl3h4eGB8+fPY9u2bdi7dy+mTp1a6nrd3Nwgk8kQGRmJEydOIDc3F926dRPf/QoAO3fuxODBg+Ht7Y2//voLJ06cwMCBA99uY4mIiKooBrsK1rx5c5iammLXrl3itF27dsHMzAwODg5ytQUFBZg7dy4sLS2hpqYGOzs7uTck5OfnY8SIEeJ8KysrLFmyRK6NYcOGoWfPnli4cCGMjY2hp6cHHx8f5OXllavfgwcPRmBgIFxdXUuti4+Px6JFi7Bu3boi87Zt24ZmzZohMDAQDRo0QPv27TF//nwsW7YMT548Kba9EydO4NatWwgJCYGtrS1sbW2xfv16nD17FpGRkQCAly9fYsKECViwYAHGjBmDRo0aoXHjxujbt2+5tpGIiKiqY7CrBMOHD0dwcLD4fd26dfD29i5SN3fuXGzYsAErVqzA5cuX4efnh0GDBiE6OhrAq+BXt25d7NixA1euXEFgYCCmTZuG7du3y7Vz9OhRJCUl4ejRo1i/fj1CQkIQEhIizp85cyYsLCzeebuePXuGgQMHYtmyZTAyMioyPycnB6qqqnLT1NTU8OLFC8TFxRXbZk5ODmQyGVRUVMRpqqqqUFBQwPHjxwEA586dw71796CgoAAHBwfxNV+XLl16520iIiKqShjsKsGgQYNw/Phx3L59G7dv38aJEycwaNAguZqcnBzMmTMH69atg7u7O+rVq4dhw4Zh0KBBWLlyJQCgRo0amDVrFlq0aAFLS0t4eXnB29u7SLDT1dXF0qVLYW1tjc8++wyenp6IiIgQ59euXbvIK7Xehp+fH1xcXNCjR49i57u7u+PkyZPYsmUL8vPzce/ePcyePRsAkJKSUuwyn3zyCTQ0NDBlyhQ8e/YM2dnZ+Prrr5Gfny8uc+PGDQCvAur06dOxf/9+6OrqokOHDkhPT3/n7SIiIqoqGOwqgb6+Pjw9PRESEoLg4GB4enqidu3acjXXr1/Hs2fP0LlzZ2hqaoqfDRs2ICkpSaxbtmwZHB0doa+vD01NTaxatQrJyclybTVp0kTufabGxsZ48OCB+N3X11cu6L2NvXv3IjIyEkFBQSXWuLm5iadLVVRU0KhRI3Tt2hXAq/fFFkdfXx87duzAvn37oKmpCW1tbWRkZKB58+biMoXX2n377bfo06cPHB0dERwcDJlMhh07drzTdhEREVUlfKVYJRk+fDh8fX0BvApnb3r69CmAV3d61qlTR25e4WnJrVu34uuvv8aiRYvg7OyMmjVrYsGCBYiNjZWrr1Gjhtx3mUwmd+PB+xAZGYmkpCTo6OjITe/Tpw/atm2LqKgoAIC/vz/8/PyQkpICXV1d3Lp1CwEBAahXr16Jbbu5uSEpKQn//PMPlJSUoKOjAyMjI3EZY2NjAEDjxo3FZVRUVFCvXr0iIZeIiEjKGOwqiYeHB3JzcyGTyeDu7l5kfuPGjaGiooLk5GS0b9++2DZOnDgBFxcXjBs3Tpz2+tG8ijR16lSMHDlSbpqtrS0WL16Mbt26yU2XyWQwMTEBAGzZsgWmpqZo3rz5f66j8KhmZGQkHjx4gO7duwMAHB0doaKigsTERLRp0wYAkJeXh1u3bsHc3Pydt42IiKiqYLCrJIqKikhISBD//KaaNWvi66+/hp+fHwoKCtCmTRtkZmbixIkT0NLSwtChQ9GwYUNs2LABhw4dgqWlJX777TecOXMGlpaW5erL0qVLsXv37lJPx6anpyM5ORn3798HACQmJgIAjIyM5D5vMjMzk+vPggUL4OHhAQUFBezatQs//vgjtm/fLo7BvXv30KlTJ2zYsAGtWrUCAAQHB8PGxgb6+vqIiYnBhAkT4OfnBysrKwCAlpYWxowZgxkzZsDU1BTm5uZYsGABAOCLL74o11gQERFVZQx2lUhLS6vU+d9//z309fUxd+5c3LhxAzo6OmjevDmmTZsGAPjyyy9x/vx59OvXDzKZDAMGDMC4ceNw8ODBcvXjn3/++c8jfXv37pW7c7d///4AgBkzZmDmzJllXtfBgwfxv//9Dzk5ObCzs8OePXvQpUsXcX5eXh4SExPx7NkzcVpiYiICAgKQnp4OCwsLfPvtt/Dz85Nrd8GCBVBSUsLgwYPx/PlzODk5ITIyErq6umXuGxERUVUnEwRBqOxOVJasrCxoa2sjMzPzP0MWERFVHRZTD1R2F6gaK8h5hjtBfSslX/CIHRFRNccQRCQdfNwJERERkUSUK9jNnTsXLVu2RM2aNWFgYICePXuKF9EX6tChA2QymdznzZfeJycnw9PTE+rq6jAwMMDkyZOLvAg+KioKzZs3h4qKCho0aCD3poRCy5Ytg4WFBVRVVeHk5ITTp0+XZ3OIiIiIJKVcwS46Oho+Pj44deoUwsPDkZeXBzc3N2RnZ8vVjRo1CikpKeJn/vz54rz8/Hx4enoiNzcXJ0+eFF9xFRgYKNbcvHkTnp6e6NixI+Lj4zFx4kSMHDkShw4dEmu2bdsGf39/zJgxA+fOnYOdnR3c3d3lHrxLREREVJ28080TDx8+hIGBAaKjo9GuXTsAr47Y2dvbl/gGgoMHD+Kzzz7D/fv3YWhoCABYsWIFpkyZgocPH0JZWRlTpkzBgQMH5N712b9/f2RkZCAsLAwA4OTkhJYtW2Lp0qUAXr19wNTUFF999RWmTp1apv5X5ZsnZs6cidDQUMTHx5d5maioKHTs2BGPHz8u8iBhIqq+eI0d0ftVmTdPvNM1dpmZmQCAWrVqyU3ftGkTateujaZNmyIgIEDu0RUxMTGwtbUVQx3w6h2iWVlZuHz5sljj6uoq16a7uztiYmIAALm5uYiLi5OrUVBQgKurq1hTnJycHGRlZcl9KtLDhw8xduxYmJmZQUVFBUZGRnB3d8eJEyfEGplMhtDQ0Art14cUFRWFHj16wNjYGBoaGrC3t8emTZvkakJCQoqcvldVVS1T2zxdT0RE9K+3viu2oKAAEydOROvWrdG0aVNx+sCBA2Fubg4TExNcuHABU6ZMQWJiInbt2gUASE1NlQt1AMTvqamppdZkZWXh+fPnePz4MfLz84utuXr1aol9njt3LmbNmvW2m/zO+vTpg9zcXKxfvx716tVDWloaIiIi8OjRo0rr04d28uRJNGvWDFOmTIGhoSH279+PIUOGQFtbG5999plYp6WlJXe9pkwmK7XdwtP1Y8aMwaZNmxAREYGRI0fC2NhYfJNH4en6FStWwMnJCUFBQXB3d0diYiIMDAw+zAYTERFVorc+Yufj44NLly5h69atctNHjx4Nd3d32NrawsvLCxs2bMDu3bsr7VVXrwsICEBmZqb4uXPnToWtOyMjA3/++SfmzZuHjh07wtzcHK1atUJAQID4aiwLCwsAQK9evSCTycTvAPDjjz/C0NAQNWvWxIgRI/DixYv/XOcff/yBRo0aQU1NDR07dsStW7eK1Bw/fhxt27aFmpoaTE1NMX78ePGayWnTpsHJyanIMnZ2dpg9e3aZtnvatGn4/vvv4eLigvr162PChAnw8PAQg34hmUwm9waLN0P7m1asWAFLS0ssWrQINjY28PX1xeeff47FixeLNT/99BNGjRoFb29vNG7cGCtWrIC6ujrWrVtXpr4TERFVNW8V7Hx9fbF//34cPXoUdevWLbW2MBhcv34dwKtXUKWlpcnVFH4vfCVVSTVaWlpQU1ND7dq1oaioWGxNca+1KqSiogItLS25T0XR1NSEpqYmQkNDkZOTU2zNmTNnALx6hVZKSor4ffv27Zg5cybmzJmDs2fPwtjYGL/++mup67tz5w569+6Nbt26IT4+HiNHjixy7WFSUhI8PDzQp08fXLhwAdu2bcPx48fh6+sLAPDy8sLp06flQvnly5dx4cIFDBw4EMCrQFbcKdDSZGZmFjl9//TpU5ibm8PU1BQ9evQQT8uX5EOdriciIqrKyhXsBEGAr68vdu/ejcjIyDK9k7Tw4n5jY2MAgLOzMy5evCh392p4eDi0tLTQuHFjsebN95aGh4fD2dkZAKCsrAxHR0e5moKCAkRERIg1HxslJSWEhIRg/fr10NHRQevWrTFt2jRcuHBBrNHX1wcA6OjowMjISPweFBSEESNGYMSIEbCyssIPP/wgjlVJli9fjvr162PRokWwsrKCl5cXhg0bJlczd+5ceHl5YeLEiWjYsCFcXFzw888/Y8OGDXjx4gWaNGkCOzs7bN68WVxm06ZNcHJyQoMGDQAAVlZW0NbWLvM4bN++HWfOnJF7PZmVlRXWrVuHPXv2YOPGjSgoKICLiwvu3r1bYjv/dbr+n3/+KfF0feEpfyIiIqkpV7Dz8fHBxo0bsXnzZtSsWROpqalITU3F8+fPAbw6AvT9998jLi4Ot27dwt69ezFkyBC0a9cOzZo1AwC4ubmhcePGGDx4MP766y8cOnQI06dPh4+PD1RUVAAAY8aMwY0bN/DNN9/g6tWr+PXXX7F9+3a594P6+/tj9erVWL9+PRISEjB27FhkZ2fLBYaPTZ8+fXD//n3s3bsXHh4e4sX//3XEKyEhocgp0f8KsGVZ5q+//kJISIh4NFFTUxPu7u4oKCjAzZs3Abw6alcY7ARBwJYtW+Dl5SW2cfXqVfTq1avUvhQ6evQovL29sXr1ajRp0kSuX0OGDIG9vT3at2+PXbt2QV9fHytXrixTu0RERPRKuW6eWL58OYBXjzR5XXBwMIYNGwZlZWUcOXIEQUFByM7OhqmpKfr06YPp06eLtYqKiti/fz/Gjh0LZ2dnaGhoYOjQoXLXbFlaWuLAgQPw8/PDkiVLULduXaxZs0a8KB4A+vXrh4cPHyIwMBCpqamwt7dHWFjYf16bVdlUVVXRuXNndO7cGd999x1GjhyJGTNmFDmaVhGePn2KL7/8EuPHjy8yz8zMDAAwYMAATJkyBefOncPz589x584d9OvXr9zrio6ORrdu3bB48WIMGTKk1NoaNWrAwcFBPH1fnP86Xa+oqPhWp+uJiIiqsnIFu/965J2pqSmio6P/sx1zc3P88ccfpdZ06NAB58+fL7XG19dXvB6sqmrcuLHc401q1KiB/Px8uRobGxvExsbKBaJTp06V2q6NjQ327t0rN+3NZZo3b44rV66Ip1WLU7duXbRv3x6bNm3C8+fP0blz53LfURoVFYXPPvsM8+bNw+jRo/+zPj8/HxcvXkTXrl1LrHF2di7yGyrpdH3Pnj0B/Hu6vqr/ZoiIiErCd8VWkEePHuHTTz/Fxo0bceHCBdy8eRM7duzA/Pnz0aNHD7HOwsICERERSE1NxePHjwEAEyZMwLp16xAcHIy///4bM2bMKHJzwe7du2FtbS1+HzNmDK5du4bJkycjMTERmzdvLnLKd8qUKTh58iR8fX0RHx+Pa9euYc+ePUWCj5eXF7Zu3YodO3bInYYFAGtra+zevbvE7T569Cg8PT0xfvx49OnTRzx9n56eLtbMnj0bhw8fxo0bN3Du3DkMGjQIt2/fxsiRI8WagIAAuWAr1dP1RERE74LBroJoamrCyckJixcvRrt27dC0aVN89913GDVqlPj2DABYtGgRwsPDYWpqCgcHBwCvTjt/9913+Oabb+Do6Ijbt29j7Nixcu1nZmbKPQfOzMwMO3fuRGhoKOzs7LBixQrMmTNHbplmzZohOjoaf//9N9q2bQsHBwcEBgbCxMREru7zzz/Ho0eP8OzZM/HoV6HExETxQdXFWb9+PZ49e4a5c+fC2NhY/PTu3Vusefz4MUaNGgUbGxt07doVWVlZOHnypNwNIikpKUhOTha/F56uDw8Ph52dHRYtWlTs6fqFCxciMDAQ9vb2iI+PrxKn64mIiN7WO71SrKqryq8UIyJ6X/hKMaL3q8q+UoyIiIiIPh4MdkREREQSwWBHREREJBHletwJEVF1xmvRiOhjxyN2H7EOHTpg4sSJld0NIiIiqiIY7CrQsGHDIJPJMGbMmCLzfHx8IJPJ5N5AsWvXLnz//fcfvF9z585Fy5YtUbNmTRgYGKBnz55yj04BXr0urlevXtDX14eWlhb69u1b5K0Ob9MuAMTExODTTz+FhoYGtLS00K5dO/E1dURERFR2DHYVzNTUFFu3bpULLi9evMDmzZvF13gVqlWrFmrWrPnB+xQdHQ0fHx+cOnUK4eHhyMvLg5ubG7KzswEA2dnZcHNzg0wmQ2RkJE6cOIHc3Fx069YNBQUFb90u8CrUeXh4wM3NDadPn8aZM2fg6+sLBQX+NImIiMqL19hVsObNmyMpKQm7du0S3+Kwa9cumJmZwdLSUq62Q4cOsLe3R1BQEIBXb6UYPXo0rl+/jh07dkBXVxfTp08v02u6ShMWFib3PSQkBAYGBoiLi0O7du1w4sQJ3Lp1C+fPnxefx7N+/Xro6uoiMjISrq6ub9UuAPj5+WH8+PGYOnWqWGdlZfVO20NERFRd8bBIJRg+fDiCg4PF7+vWrSvza64WLVqEFi1a4Pz58xg3bhzGjh1b7OnNQiEhIZDJZOXqX+GbJGrVqgUAyMnJgUwmg4qKilijqqoKBQUFHD9+/K3bffDgAWJjY2FgYAAXFxcYGhqiffv25WqTiIiI/sVgVwkGDRqE48eP4/bt27h9+zZOnDiBQYMGlWnZrl27Yty4cWjQoAGmTJmC2rVr4+jRoyXWa2trl+sIWEFBASZOnIjWrVujadOmAIBPPvkEGhoamDJlCp49e4bs7Gx8/fXXyM/PR0pKylu3e+PGDQDAzJkzMWrUKISFhaF58+bo1KkTrl27VuY+ExER0SsMdpVAX18fnp6eCAkJQXBwMDw9PVG7du0yLdusWTPxzzKZDEZGRnjw4EGJ9b169cLVq1fL3DcfHx9cunQJW7dulevvjh07sG/fPmhqakJbWxsZGRlo3rx5ma+FK67dwuvzvvzyS3h7e8PBwQGLFy+GlZUV1q1bV+Y+ExER0Su8xq6SDB8+HL6+vgCAZcuWlXm5GjVqyH2XyWSl3sBQHr6+vti/fz+OHTuGunXrys1zc3NDUlIS/vnnHygpKUFHRwdGRkaoV6/eW7drbGwMAGjcuLFcvY2NDZKTk9/DFhEREVUvPGJXSTw8PJCbm4u8vDy4u7tXal8EQYCvry92796NyMjIIjdxvK527drQ0dFBZGQkHjx4gO7du791uxYWFjAxMSlyjeDff/8Nc3Pzd9soIiKiaojBrpIoKioiISEBV65cgaKi4ntrt1OnTli6dKn4fffu3bC2ti51GR8fH2zcuBGbN29GzZo1kZqaitTUVLlHsgQHB+PUqVNISkrCxo0b8cUXX8DPz0/u+r031/1f7cpkMkyePBk///wzfv/9d1y/fh3fffcdrl69ihEjRryvISEiIqo2eCq2EhU+OuR9KjxdWigzM7PUu2YBYPny5QBePV7ldcHBweIDkxMTExEQEID09HRYWFjg22+/hZ+fX6nrLku7EydOxIsXL+Dn54f09HTY2dkhPDwc9evXL+smExER0f+TCYIgVHYnKktWVha0tbWRmZn5QUIWEUkL3xVLRGVRkPMMd4L6Vkq+4KlYIiIiIolgsCMiIiKSCF5jR0TvHU9ZEhFVDh6xIyIiIpIIBjsiIiIiiWCwIyIiIpIIBrsKNGzYMMhkMowZM6bIPB8fH8hkMvH5bpVp165d6Ny5M/T19aGlpQVnZ2ccOnSoSN29e/cwaNAg6OnpQU1NDba2tjh79qxcTUJCArp37w5tbW1oaGigZcuWpb4uLC8vD7Nnz0b9+vWhqqoKOzs7hIWFydVYWFhAJpMV+fj4+LyfASAiIqqiGOwqmKmpKbZu3Sr3VocXL15g8+bNMDMzq8Se/evYsWPo3Lkz/vjjD8TFxaFjx47o1q0bzp8/L9Y8fvwYrVu3Ro0aNXDw4EFcuXIFixYtgq6urliTlJSENm3awNraGlFRUbhw4QK+++47qKqqlrju6dOnY+XKlfjll19w5coVjBkzBr169ZJb95kzZ5CSkiJ+wsPDAQBffPHFBxgNIiKiqoMPKK7ABxQPGzYMGRkZSEpKwtSpU+Hl5QUA2Lx5M+bNmwdLS0vo6OggJCQEABAWFoYffvgBly5dgqKiIpydnbFkyRLxrQwbNmzAuHHjcP78eTRs2BAAMG7cOERGRuLcuXNQV1d/b31v0qQJ+vXrh8DAQADA1KlTceLECfz5558lLtO/f3/UqFEDv/32W5nXY2Jigm+//Vbu6FufPn2gpqaGjRs3FrvMxIkTsX//fly7dg0ymazM66IPh3fFElF1xgcUVzPDhw9HcHCw+H3dunXw9vYuUpednQ1/f3+cPXsWERERUFBQQK9evVBQUAAAGDJkCLp27QovLy+8fPkSBw4cwJo1a7Bp06YSQ92tW7cgk8kQFRVV5v4WFBTgyZMnqFWrljht7969aNGiBb744gsYGBjAwcEBq1evllvmwIEDaNSoEdzd3WFgYAAnJyeEhoaWuq6cnJwiR/TU1NRw/PjxYutzc3OxceNGDB8+nKGOiIiqPQa7SjBo0CAcP34ct2/fxu3bt3HixAkMGjSoSF2fPn3Qu3dvNGjQAPb29li3bh0uXryIK1euiDUrV65ESkoKxo8fjxEjRmDmzJlwdHQscd01atSAlZVVuY7mLVy4EE+fPkXfvn3FaTdu3MDy5cvRsGFDHDp0CGPHjsX48eOxfv16AMCDBw/w9OlT/Pjjj/Dw8MDhw4fRq1cv9O7dG9HR0SWuy93dHT/99BOuXbuGgoIChIeHY9euXUhJSSm2PjQ0FBkZGR/FtYlERESVjQ8orgT6+vrw9PRESEgIBEGAp6cnateuXaTu2rVrCAwMRGxsLP755x/xSF1ycjKaNm0KANDV1cXatWvh7u4OFxcXTJ06tdR116lTB1evXi1zXzdv3oxZs2Zhz549MDAwEKcXFBSgRYsWmDNnDgDAwcEBly5dwooVKzB06FCxrz169ICfnx8AwN7eHidPnsSKFSvQvn37Yte3ZMkSjBo1CtbW1pDJZKhfvz68vb2xbt26YuvXrl2LLl26wMTEpMzbREREJFU8YldJhg8fjpCQEKxfvx7Dhw8vtqZbt25IT0/H6tWrERsbi9jYWACvTj++7tixY1BUVERKSgqys7PfWx+3bt2KkSNHYvv27XB1dZWbZ2xsjMaNG8tNs7GxEe94rV27NpSUlEqtKY6+vj5CQ0ORnZ2N27dv4+rVq9DU1ES9evWK1N6+fRtHjhzByJEj33YTiYiIJIXBrpJ4eHggNzcXeXl5cHd3LzL/0aNHSExMxPTp09GpUyfY2Njg8ePHRepOnjyJefPmYd++fdDU1ISvr+976d+WLVvg7e2NLVu2wNPTs8j81q1bIzExUW7a33//DXNzcwCAsrIyWrZsWWpNaVRVVVGnTh28fPkSO3fuRI8ePYrUBAcHw8DAoNj+ERERVUc8FVtJFBUVkZCQIP75Tbq6utDT08OqVatgbGyM5OTkIqdZnzx5gsGDB2P8+PHo0qUL6tati5YtW6Jbt274/PPPAQABAQG4d+8eNmzYAODVs+c6deqEDRs2oFWrVsX2bfPmzRg6dCiWLFkCJycnpKamAnh1E4O2tjYAwM/PDy4uLpgzZw769u2L06dPY9WqVVi1apXYzuTJk9GvXz+0a9cOHTt2RFhYGPbt2yd348aQIUNQp04dzJ07FwAQGxuLe/fuwd7eHvfu3cPMmTNRUFCAb775Rq6PBQUFCA4OxtChQ6GkxJ8xERERwCN2lUpLS6vE26AVFBSwdetWxMXFoWnTpvDz88OCBQvkaiZMmAANDQ3xOjdbW1vMmTMHX375Je7duwcASElJkTv1mZeXh8TERDx79qzEfq1atQovX76Ej48PjI2Nxc+ECRPEmpYtW2L37t3YsmULmjZtiu+//x5BQUHiI1wAoFevXlixYgXmz58PW1tbrFmzBjt37kSbNm3EmuTkZLkbI168eIHp06ejcePG6NWrF+rUqYPjx49DR0dHro9HjhxBcnJyiaexiYiIqiM+x64Cn2NHVF3wOXZEVJ3xOXZERERE9M4Y7IiIiIgkgsGOiIiISCIY7IiIiIgkgsGOiIiISCIY7IiIiIgkgsGOiIiISCIY7IiIiIgkgsGOiIiISCIY7IiIiIgkgm9PJ6pkfP0WERG9LzxiR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREElGuYDd37ly0bNkSNWvWhIGBAXr27InExES5mhcvXsDHxwd6enrQ1NREnz59kJaWJleTnJwMT09PqKurw8DAAJMnT8bLly/laqKiotC8eXOoqKigQYMGCAkJKdKfZcuWwcLCAqqqqnBycsLp06fLszlEREREklKuYBcdHQ0fHx+cOnUK4eHhyMvLg5ubG7Kzs8UaPz8/7Nu3Dzt27EB0dDTu37+P3r17i/Pz8/Ph6emJ3NxcnDx5EuvXr0dISAgCAwPFmps3b8LT0xMdO3ZEfHw8Jk6ciJEjR+LQoUNizbZt2+Dv748ZM2bg3LlzsLOzg7u7Ox48ePAu40FERERUZckEQRDeduGHDx/CwMAA0dHRaNeuHTIzM6Gvr4/Nmzfj888/BwBcvXoVNjY2iImJwSeffIKDBw/is88+w/3792FoaAgAWLFiBaZMmYKHDx9CWVkZU6ZMwYEDB3Dp0iVxXf3790dGRgbCwsIAAE5OTmjZsiWWLl0KACgoKICpqSm++uorTJ06tUz9z8rKgra2NjIzM6GlpfW2w0D0TviAYiIiaSnIeYY7QX0rJV+80zV2mZmZAIBatWoBAOLi4pCXlwdXV1exxtraGmZmZoiJiQEAxMTEwNbWVgx1AODu7o6srCxcvnxZrHm9jcKawjZyc3MRFxcnV6OgoABXV1expjg5OTnIysqS+xARERFJxVsHu4KCAkycOBGtW7dG06ZNAQCpqalQVlaGjo6OXK2hoSFSU1PFmtdDXeH8wnml1WRlZeH58+f4559/kJ+fX2xNYRvFmTt3LrS1tcWPqalp+TeciIiI6CP11sHOx8cHly5dwtatW99nfz6ogIAAZGZmip87d+5UdpeIiIiI3hult1nI19cX+/fvx7Fjx1C3bl1xupGREXJzc5GRkSF31C4tLQ1GRkZizZt3rxbeNft6zZt30qalpUFLSwtqampQVFSEoqJisTWFbRRHRUUFKioq5d9gIiIioiqgXEfsBEGAr68vdu/ejcjISFhaWsrNd3R0RI0aNRARESFOS0xMRHJyMpydnQEAzs7OuHjxotzdq+Hh4dDS0kLjxo3FmtfbKKwpbENZWRmOjo5yNQUFBYiIiBBriIiIiKqbch2x8/HxwebNm7Fnzx7UrFlTvJ5NW1sbampq0NbWxogRI+Dv749atWpBS0sLX331FZydnfHJJ58AANzc3NC4cWMMHjwY8+fPR2pqKqZPnw4fHx/xaNqYMWOwdOlSfPPNNxg+fDgiIyOxfft2HDjw792D/v7+GDp0KFq0aIFWrVohKCgI2dnZ8Pb2fl9jQ0RERFSllCvYLV++HADQoUMHuenBwcEYNmwYAGDx4sVQUFBAnz59kJOTA3d3d/z6669iraKiIvbv34+xY8fC2dkZGhoaGDp0KGbPni3WWFpa4sCBA/Dz88OSJUtQt25drFmzBu7u7mJNv3798PDhQwQGBiI1NRX29vYICwsrckMFERERUXXxTs+xq+r4HDv6GPA5dkRE0lJln2NHRERERB+Pt7orlqgy8MgWERFR6XjEjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgiGOyIiIiIJILBjoiIiEgilCq7A/T+WUw9UNldICIiokrAI3ZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEsFgR0RERCQRDHZEREREEqFU2R34GDSdcQgKKuqV3Q0iIiKid8IjdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBHlDnbHjh1Dt27dYGJiAplMhtDQULn5w4YNg0wmk/t4eHjI1aSnp8PLywtaWlrQ0dHBiBEj8PTpU7maCxcuoG3btlBVVYWpqSnmz59fpC87duyAtbU1VFVVYWtriz/++KO8m0NEREQkGeUOdtnZ2bCzs8OyZctKrPHw8EBKSor42bJli9x8Ly8vXL58GeHh4di/fz+OHTuG0aNHi/OzsrLg5uYGc3NzxMXFYcGCBZg5cyZWrVol1pw8eRIDBgzAiBEjcP78efTs2RM9e/bEpUuXyrtJRERERJIgEwRBeOuFZTLs3r0bPXv2FKcNGzYMGRkZRY7kFUpISEDjxo1x5swZtGjRAgAQFhaGrl274u7duzAxMcHy5cvx7bffIjU1FcrKygCAqVOnIjQ0FFevXgUA9OvXD9nZ2di/f7/Y9ieffAJ7e3usWLGiTP3PysqCtrY2TCduh4KK+luMABEREZG8gpxnuBPUF5mZmdDS0qrQdX+Qa+yioqJgYGAAKysrjB07Fo8ePRLnxcTEQEdHRwx1AODq6goFBQXExsaKNe3atRNDHQC4u7sjMTERjx8/FmtcXV3l1uvu7o6YmJgPsUlEREREHz2l992gh4cHevfuDUtLSyQlJWHatGno0qULYmJioKioiNTUVBgYGMh3QkkJtWrVQmpqKgAgNTUVlpaWcjWGhobiPF1dXaSmporTXq8pbKM4OTk5yMnJEb9nZWW907YSERERfUzee7Dr37+/+GdbW1s0a9YM9evXR1RUFDp16vS+V1cuc+fOxaxZsyq1D0REREQfygd/3Em9evVQu3ZtXL9+HQBgZGSEBw8eyNW8fPkS6enpMDIyEmvS0tLkagq//1dN4fziBAQEIDMzU/zcuXPn3TaOiIiI6CPywYPd3bt38ejRIxgbGwMAnJ2dkZGRgbi4OLEmMjISBQUFcHJyEmuOHTuGvLw8sSY8PBxWVlbQ1dUVayIiIuTWFR4eDmdn5xL7oqKiAi0tLbkPERERkVSUO9g9ffoU8fHxiI+PBwDcvHkT8fHxSE5OxtOnTzF58mScOnUKt27dQkREBHr06IEGDRrA3d0dAGBjYwMPDw+MGjUKp0+fxokTJ+Dr64v+/fvDxMQEADBw4EAoKytjxIgRuHz5MrZt24YlS5bA399f7MeECRMQFhaGRYsW4erVq5g5cybOnj0LX1/f9zAsRERERFVPuYPd2bNn4eDgAAcHBwCAv78/HBwcEBgYCEVFRVy4cAHdu3dHo0aNMGLECDg6OuLPP/+EioqK2MamTZtgbW2NTp06oWvXrmjTpo3cM+q0tbVx+PBh3Lx5E46Ojpg0aRICAwPlnnXn4uKCzZs3Y9WqVbCzs8Pvv/+O0NBQNG3a9F3Gg4iIiKjKeqfn2FV1fI4dERERvW+Se44dEREREVU8BjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpIIBjsiIiIiiWCwIyIiIpKIcge7Y8eOoVu3bjAxMYFMJkNoaKjcfEEQEBgYCGNjY6ipqcHV1RXXrl2Tq0lPT4eXlxe0tLSgo6ODESNG4OnTp3I1Fy5cQNu2baGqqgpTU1PMnz+/SF927NgBa2trqKqqwtbWFn/88Ud5N4eIiIhIMsod7LKzs2FnZ4dly5YVO3/+/Pn4+eefsWLFCsTGxkJDQwPu7u548eKFWOPl5YXLly8jPDwc+/fvx7FjxzB69GhxflZWFtzc3GBubo64uDgsWLAAM2fOxKpVq8SakydPYsCAARgxYgTOnz+Pnj17omfPnrh06VJ5N4mIiIhIEmSCIAhvvbBMht27d6Nnz54AXh2tMzExwaRJk/D1118DADIzM2FoaIiQkBD0798fCQkJaNy4Mc6cOYMWLVoAAMLCwtC1a1fcvXsXJiYmWL58Ob799lukpqZCWVkZADB16lSEhobi6tWrAIB+/fohOzsb+/fvF/vzySefwN7eHitWrChT/7OysqCtrQ3TiduhoKL+tsNAREREJCrIeYY7QX2RmZkJLS2tCl33e73G7ubNm0hNTYWrq6s4TVtbG05OToiJiQEAxMTEQEdHRwx1AODq6goFBQXExsaKNe3atRNDHQC4u7sjMTERjx8/FmteX09hTeF6iIiIiKobpffZWGpqKgDA0NBQbrqhoaE4LzU1FQYGBvKdUFJCrVq15GosLS2LtFE4T1dXF6mpqaWupzg5OTnIyckRv2dlZZVn84iIiIg+atXqrti5c+dCW1tb/JiamlZ2l4iIiIjem/ca7IyMjAAAaWlpctPT0tLEeUZGRnjw4IHc/JcvXyI9PV2uprg2Xl9HSTWF84sTEBCAzMxM8XPnzp3ybiIRERHRR+u9BjtLS0sYGRkhIiJCnJaVlYXY2Fg4OzsDAJydnZGRkYG4uDixJjIyEgUFBXBychJrjh07hry8PLEmPDwcVlZW0NXVFWteX09hTeF6iqOiogItLS25DxEREZFUlDvYPX36FPHx8YiPjwfw6oaJ+Ph4JCcnQyaTYeLEifjhhx+wd+9eXLx4EUOGDIGJiYl456yNjQ08PDwwatQonD59GidOnICvry/69+8PExMTAMDAgQOhrKyMESNG4PLly9i2bRuWLFkCf39/sR8TJkxAWFgYFi1ahKtXr2LmzJk4e/YsfH19331UiIiIiKqgct88cfbsWXTs2FH8Xhi2hg4dipCQEHzzzTfIzs7G6NGjkZGRgTZt2iAsLAyqqqriMps2bYKvry86deoEBQUF9OnTBz///LM4X1tbG4cPH4aPjw8cHR1Ru3ZtBAYGyj3rzsXFBZs3b8b06dMxbdo0NGzYEKGhoWjatOlbDQQRERFRVfdOz7Gr6vgcOyIiInrfJPMcOyIiIiKqPAx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEQx2RERERBLBYEdEREQkEe892M2cORMymUzuY21tLc5/8eIFfHx8oKenB01NTfTp0wdpaWlybSQnJ8PT0xPq6uowMDDA5MmT8fLlS7maqKgoNG/eHCoqKmjQoAFCQkLe96YQERERVSkf5IhdkyZNkJKSIn6OHz8uzvPz88O+ffuwY8cOREdH4/79++jdu7c4Pz8/H56ensjNzcXJkyexfv16hISEIDAwUKy5efMmPD090bFjR8THx2PixIkYOXIkDh069CE2h4iIiKhKUPogjSopwcjIqMj0zMxMrF27Fps3b8ann34KAAgODoaNjQ1OnTqFTz75BIcPH8aVK1dw5MgRGBoawt7eHt9//z2mTJmCmTNnQllZGStWrIClpSUWLVoEALCxscHx48exePFiuLu7f4hNIiIiIvrofZAjdteuXYOJiQnq1asHLy8vJCcnAwDi4uKQl5cHV1dXsdba2hpmZmaIiYkBAMTExMDW1haGhoZijbu7O7KysnD58mWx5vU2CmsK2yAiIiKqjt77ETsnJyeEhITAysoKKSkpmDVrFtq2bYtLly4hNTUVysrK0NHRkVvG0NAQqampAIDU1FS5UFc4v3BeaTVZWVl4/vw51NTUiu1bTk4OcnJyxO9ZWVnvtK1EREREH5P3Huy6dOki/rlZs2ZwcnKCubk5tm/fXmLgqihz587FrFmzKrUPRERERB/KB3/ciY6ODho1aoTr16/DyMgIubm5yMjIkKtJS0sTr8kzMjIqcpds4ff/qtHS0io1PAYEBCAzM1P83Llz5103j4iIiOij8cGD3dOnT5GUlARjY2M4OjqiRo0aiIiIEOcnJiYiOTkZzs7OAABnZ2dcvHgRDx48EGvCw8OhpaWFxo0bizWvt1FYU9hGSVRUVKClpSX3ISIiIpKK9x7svv76a0RHR+PWrVs4efIkevXqBUVFRQwYMADa2toYMWIE/P39cfToUcTFxcHb2xvOzs745JNPAABubm5o3LgxBg8ejL/++guHDh3C9OnT4ePjAxUVFQDAmDFjcOPGDXzzzTe4evUqfv31V2zfvh1+fn7ve3OIiIiIqoz3fo3d3bt3MWDAADx69Aj6+vpo06YNTp06BX19fQDA4sWLoaCggD59+iAnJwfu7u749ddfxeUVFRWxf/9+jB07Fs7OztDQ0MDQoUMxe/ZsscbS0hIHDhyAn58flixZgrp162LNmjV81AkRERFVazJBEITK7kRlycrKgra2NkwnboeCinpld4eIiIgkoCDnGe4E9UVmZmaFX/bFd8USERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSQSDHREREZFEMNgRERERSUSVD3bLli2DhYUFVFVV4eTkhNOnT1d2l4iIiIgqRZUOdtu2bYO/vz9mzJiBc+fOwc7ODu7u7njw4EFld42IiIiowlXpYPfTTz9h1KhR8Pb2RuPGjbFixQqoq6tj3bp1ld01IiIiogpXZYNdbm4u4uLi4OrqKk5TUFCAq6srYmJiKrFnRERERJVDqbI78Lb++ecf5Ofnw9DQUG66oaEhrl69WuwyOTk5yMnJEb9nZmYCAApynn24jhIREVG1UpgrBEGo8HVX2WD3NubOnYtZs2YVmX5v+bCK7wwRERFJ2qNHj6CtrV2h66yywa527dpQVFREWlqa3PS0tDQYGRkVu0xAQAD8/f3F7xkZGTA3N0dycnKFD/zHJisrC6amprhz5w60tLQquzuVimPxL47FvzgW/+JY/Itj8S+Oxb8yMzNhZmaGWrVqVfi6q2ywU1ZWhqOjIyIiItCzZ08AQEFBASIiIuDr61vsMioqKlBRUSkyXVtbu9r/CAtpaWlxLP4fx+JfHIt/cSz+xbH4F8fiXxyLfykoVPytDFU22AGAv78/hg4dihYtWqBVq1YICgpCdnY2vL29K7trRERERBWuSge7fv364eHDhwgMDERqairs7e0RFhZW5IYKIiIiouqgSgc7APD19S3x1Ot/UVFRwYwZM4o9PVvdcCz+xbH4F8fiXxyLf3Es/sWx+BfH4l+VORYyoTLuxSUiIiKi967KPqCYiIiIiOQx2BERERFJBIMdERERkURUyWB37NgxdOvWDSYmJpDJZAgNDRXn5eXlYcqUKbC1tYWGhgZMTEwwZMgQ3L9/X66N9PR0eHl5QUtLCzo6OhgxYgSePn1a6npfvHgBHx8f6OnpQVNTE3369CnygOSK9q5jcevWLYwYMQKWlpZQU1ND/fr1MWPGDOTm5pa63g4dOkAmk8l9xowZ86E2s8zex2/DwsKiyLb9+OOPpa5Xir+NqKioIuNQ+Dlz5kyJ6/0YfxuljQUAzJw5E9bW1tDQ0ICuri5cXV0RGxsrV1Md9hnAf4+FlPYZ7+N3UR32F8B/j0V12l+8bsyYMZDJZAgKCpKbXpn7iyoZ7LKzs2FnZ4dly5YVmffs2TOcO3cO3333Hc6dO4ddu3YhMTER3bt3l6vz8vLC5cuXER4ejv379+PYsWMYPXp0qev18/PDvn37sGPHDkRHR+P+/fvo3bv3e9228nrXsbh69SoKCgqwcuVKXL58GYsXL8aKFSswbdq0/1z3qFGjkJKSIn7mz5//XrftbbyP3wYAzJ49W27bvvrqq1LXK8XfhouLi9wYpKSkYOTIkbC0tESLFi1KXffH9tsobSwAoFGjRli6dCkuXryI48ePw8LCAm5ubnj48KFYUx32GcB/j4WU9hnv43cBSH9/Afz3WFSn/UWh3bt349SpUzAxMSkyr1L3F0IVB0DYvXt3qTWnT58WAAi3b98WBEEQrly5IgAQzpw5I9YcPHhQkMlkwr1794ptIyMjQ6hRo4awY8cOcVpCQoIAQIiJiXn3DXkP3mYsijN//nzB0tKy1Hbat28vTJgw4S16WXHedjzMzc2FxYsXl3k91eW3kZubK+jr6wuzZ88utZ2P/bdRlrHIzMwUAAhHjhwRBKF67zPeHIviSGGf8bZjUV33F//1u5D6/uLu3btCnTp1hEuXLhX5DVT2/qJKHrErr8zMTMhkMujo6AAAYmJioKOjI/evCFdXVygoKBQ5zF4oLi4OeXl5cHV1FadZW1vDzMwMMTExH7T/79ObY1FSTVneb7dp0ybUrl0bTZs2RUBAAJ49e/Yee1oxShqPH3/8EXp6enBwcMCCBQvw8uXLEtuoLr+NvXv34tGjR2V6s0tV/m3k5uZi1apV0NbWhp2dHYDqu88obiyKUx32GaWNRXXbX5TldyHl/UVBQQEGDx6MyZMno0mTJkXmV/b+oso/oPi/vHjxAlOmTMGAAQPEd9elpqbCwMBArk5JSQm1atVCampqse2kpqZCWVm5yF96hoaGJS7zsSluLN50/fp1/PLLL1i4cGGpbQ0cOBDm5uYwMTHBhQsXMGXKFCQmJmLXrl0fousfREnjMX78eDRv3hy1atXCyZMnERAQgJSUFPz000/FtlNdfhtr166Fu7s76tatW2pbVfW3sX//fvTv3x/Pnj2DsbExwsPDUbt2bQDVb59R2li8Ser7jP8ai+q0vyjP70LK+4t58+ZBSUkJ48ePL3Z+Ze8vJB3s8vLy0LdvXwiCgOXLl1d2dypVWcbi3r178PDwwBdffIFRo0aV2t7r1wrY2trC2NgYnTp1QlJSEurXr/9e+/4hlDYe/v7+4p+bNWsGZWVlfPnll5g7d64kn6helt/G3bt3cejQIWzfvv0/26uqv42OHTsiPj4e//zzD1avXo2+ffsiNja2yA66OijrWFSHfcZ/jUV12l+U9Xch5f1FXFwclixZgnPnzkEmk1V2d4ol2VOxhX9Z3b59G+Hh4XJHIYyMjPDgwQO5+pcvXyI9PR1GRkbFtmdkZITc3FxkZGTITU9LSytxmY9FaWNR6P79++jYsSNcXFywatWqcq/DyckJwKt/vX/syjIer3NycsLLly9x69atYudL/bcBAMHBwdDT0yv2RpP/UlV+GxoaGmjQoAE++eQTrF27FkpKSli7di2A6rfPKG0sClWXfUZZxuJ1Ut5flHUspLy/+PPPP/HgwQOYmZlBSUkJSkpKuH37NiZNmgQLCwsAlb+/kGSwK/zL6tq1azhy5Aj09PTk5js7OyMjIwNxcXHitMjISBQUFIg/qjc5OjqiRo0aiIiIEKclJiYiOTkZzs7OH2ZD3oP/Ggvg1b+6O3ToAEdHRwQHB0NBofw/i/j4eACAsbHxu3b5gyrLeLwpPj4eCgoKJR65kfJvAwAEQUBwcDCGDBmCGjVqlHs9VeW38aaCggLk5OQAqF77jOK8PhZA9dpnvOnNsXiTVPcXxSluLKS+vxg8eDAuXLiA+Ph48WNiYoLJkyfj0KFDAD6C/UW5brX4SDx58kQ4f/68cP78eQGA8NNPPwnnz58Xbt++LeTm5grdu3cX6tatK8THxwspKSniJycnR2zDw8NDcHBwEGJjY4Xjx48LDRs2FAYMGCDOv3v3rmBlZSXExsaK08aMGSOYmZkJkZGRwtmzZwVnZ2fB2dm5Qrf9Te86Fnfv3hUaNGggdOrUSbh7965cTaE3x+L69evC7NmzhbNnzwo3b94U9uzZI9SrV09o165dpYzB6951PE6ePCksXrxYiI+PF5KSkoSNGzcK+vr6wpAhQ8R1VJffRqEjR44IAISEhIQi66gqv43SxuLp06dCQECAEBMTI9y6dUs4e/as4O3tLaioqAiXLl0S26gO+4yyjIWU9hnvOhbVZX9R1v9HBEH6+4viFHdndGXuL6pksDt69KgAoMhn6NChws2bN4udB0A4evSo2MajR4+EAQMGCJqamoKWlpbg7e0tPHnyRJxf2M7ryzx//lwYN26coKurK6irqwu9evWS25lVhncdi+Dg4BJrCr05FsnJyUK7du2EWrVqCSoqKkKDBg2EyZMnC5mZmZUwAvLedTzi4uIEJycnQVtbW1BVVRVsbGyEOXPmCC9evBDXUV1+G4UGDBgguLi4FLuOqvLbKG0snj9/LvTq1UswMTERlJWVBWNjY6F79+7C6dOn5dqoDvuMsoyFlPYZ7zoW1WV/Udb/RwRB+vuL4hQX7CpzfyETBEEo9ZAeEREREVUJkrzGjoiIiKg6YrAjIiIikggGOyIiIiKJYLAjIiIikggGOyIiIiKJYLAjIiIikggGOyIiIiKJYLAjIiIikggGOyKiMrCwsEBQUFBld4OIqFQMdkREZXDmzBmMHj26TLWCICAwMBDGxsZQU1ODq6srrl279oF7SETEYEdEVCb6+vpQV1cvU+38+fPx888/Y8WKFYiNjYWGhgbc3d3x4sWLD9xLIqruGOyIqFoICwtDmzZtoKOjAz09PXz22WdISkoCAGzYsAGamppyR9XGjRsHa2trPHv2DID8qVhBEDBz5kyYmZlBRUUFJiYmGD9+vDgvKCgI06dPR48ePdCsWTNs2LAB9+/fR2hoaIVuMxFVPwx2RFQtZGdnw9/fH2fPnkVERAQUFBTQq1cvFBQUYMiQIejatSu8vLzw8uVLHDhwAGvWrMGmTZuKPUq3c+dOLF68GCtXrsS1a9cQGhoKW1tbAMDNmzeRmpoKV1dXsV5bWxtOTk6IiYmpsO0loupJqbI7QERUEfr06SP3fd26ddDX18eVK1fQtGlTrFy5Es2aNcP48eOxa9cuzJw5E46OjsW2lZycDCMjI7i6uqJGjRowMzNDq1atAACpqakAAENDQ7llDA0NxXlERB8Kj9gRUbVw7do1DBgwAPXq1YOWlhYsLCwAvAppAKCrq4u1a9di+fLlqF+/PqZOnVpiW1988QWeP3+OevXqYdSoUdi9ezdevnxZEZtBRFQqBjsiqha6deuG9PR0rF69GrGxsYiNjQUA5ObmijXHjh2DoqIiUlJSkJ2dXWJbpqamSExMxK+//go1NTWMGzcO7dq1Q15eHoyMjAAAaWlpcsukpaWJ84iIPhQGOyKSvEePHiExMRHTp09Hp06dYGNjg8ePH8vVnDx5EvPmzcO+ffugqakJX1/fUttUU1NDt27d8PPPPyMqKgoxMTG4ePEiLC0tYWRkhIiICLE2KysLsbGxcHZ2/iDbR0RUiNfYEZHk6erqQk9PD6tWrYKxsTGSk5PlTrU+efIEgwcPxvjx49GlSxfUrVsXLVu2RLdu3fD5558XaS8kJAT5+flwcnKCuro6Nm7cCDU1NZibm0Mmk2HixIn44Ycf0LBhQ1haWuK7776DiYkJevbsWYFbTUTVEY/YEZHkKSgoYOvWrYiLi0PTpk3h5+eHBQsWiPMnTJgADQ0NzJkzBwBga2uLOXPm4Msvv8S9e/eKtKejo4PVq1ejdevWaNasGY4cOYJ9+/ZBT08PAPDNN9/gq6++wujRo9GyZUs8ffoUYWFhUFVVrZgNJqJqSyYIglDZnSAiIiKid8cjdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBEMdkREREQSwWBHREREJBH/B+Ga5tJE8J/0AAAAAElFTkSuQmCC\n" }, "metadata": {}, "output_type": "display_data" } ], "source": [ "h = physt.h1(data, range=(120, 140))\n", "h.plot(show_stats=True);" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": "Statistics(sum=149960006.08864358, sum2=23113238494.210587, min=29.26409978722107, max=266.9737275094181, weight=1000000.0)" }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h._stats" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": "229977" }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "h.total" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.8.5" } }, "nbformat": 4, "nbformat_minor": 1 }