Skip to content
Snippets Groups Projects
reshape_and_compress_ascii_with_WF.ipynb 84.9 KiB
Newer Older
{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Script to reshape and compress original ascii files \n",
    "```bash\n",
    "#!/bin/bash\n",
    "if [ $# -ne 1 ];then\n",
    "    echo \"ERROR: use 1 argument (run number)\"\n",
    "    exit 1\n",
    "fi\n",
    "\n",
    "nrun=$( printf \"%6.6d\" \"$1\")\n",
    "ascii_dir=\"/data/insudaq/ascii_asa_2021\"\n",
    "new_ascii_compact_dir=\"/home/insudaq/tlktest2011/maroc3_2012/vme14/ana_insulab2021_asa/python/compact_data_WF\"\n",
    "awk '{for(i=1;i<=NF;i++)printf(\"%s \",$i);if(NR%3==0)printf(\"\\n\")}' $ascii_dir/run$nrun*wave*.dat | gzip > $new_ascii_compact_dir\"/run\"$nrun\"_WF.gz\"\n",
    "\n",
    "```\n",
    "The script can be run from here: `/home/insudaq/valerio/bar-fiber-sipm-2021`\n",
    "\n",
    "It creates a `run260158_WF.gz` here `/home/insudaq/tlktest2011/maroc3_2012/vme14/ana_insulab2021_asa/python/compact_data_WF` folder.\n",
    "\n",
    "The gz files contains **1 row per event**, each row consisting of a total of 23 + 512 + 512 = 1047 columns:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Data import"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = np.loadtxt('compact_data_WF/run260158_WF.gz')"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
Valerio Mascagna's avatar
Valerio Mascagna committed
       "(11943, 1047)"
Valerio Mascagna's avatar
Valerio Mascagna committed
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.shape"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# reorganizing data\n",
    "\n",
    "sili     = data[:, 0  : 4]                # 4 coordinates from the silicon track detector [cm]\n",
    "digiph   = data[:, 4  : 4+8]              # 8 pulse heigth of the digitizer [ADC, fullscale=16384]\n",
    "digitime = data[:, 4+8: 4+8+8]            # 8 time of the maximum of the digitizer [2ns, range 0-511]\n",
    "dtime    = data[:, 4+8+8]                 # 1 time from the start of run [s]\n",
    "abstime  = data[:, 4+8+8+1]               # 1 absolute time from 1970-01-01 [s]\n",
    "nevent   = data[:, 4+8+8+2]               # 1 event number\n",
    "wave1    = data[:, 4+8+8+3: 4+8+8+3+512]  # 512 sampling points of 1st waveform [ADC, fullscale=16384]\n",
    "wave2    = data[:, 4+8+8+3+512:]          # 512 sampling points of 2nd waveform [ADC, fullscale=16384]"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
Valerio Mascagna's avatar
Valerio Mascagna committed
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAHSCAYAAAA5eGh0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAqSUlEQVR4nO3df7BlZX3n+/dn/DVV6oQmdHo6QOeI1bEKUzNIThFmYizmkihgSvRWioE/lCiV1grUaI1VN62ZipYpq0hu1NKaXHLbgQKqFDFBIxUw2mHMcK0KmAOD/BCRhjRFdzXdJ+IFHKecgN/5Y6+Dm8M5fX7tH8/e+/2q2nXWftbaZ3179fnu73qe9ey1U1VIkqQ2/LNxByBJkn7KwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDXjruAABOOumkmpubG3cYUvPuuuuuf6yq7eOO43jMZ2l9VsvnNQtzklOB64EdQAH7qurTST4K/A6w2G364aq6tXvNh4DLgOeA/1BVXzvePubm5lhYWNjAP0eaTUkeG3cMazGfpfVZLZ/X02N+FvhgVd2d5NXAXUn2d+s+VVV/smxHpwMXA68Hfh74myS/WFXPbT58SZJmw5rXmKvqSFXd3S0/AzwInHycl1wIfKGqflxV/wAcAM4aRLCSJE27DU3+SjIHvAG4s2u6Ism9Sa5Jsq1rOxl4vO9lhzh+IZckSZ11F+YkrwJuAj5QVU8DVwGvBc4AjgCf2MiOk+xJspBkYXFxce0XSJI0A9ZVmJO8jF5R/lxVfQmgqo5W1XNV9RPgs/x0uPowcGrfy0/p2l6gqvZV1XxVzW/f3vQkU0mSRmbNwpwkwNXAg1X1yb72nX2bvQO4v1u+Gbg4ySuSvAbYDXxrcCFL2qwkpyb5RpLvJHkgyfu79o8mOZzknu5xQd9rPpTkQJKHkrxlfNFLs2E9s7J/FXgncF+Se7q2DwOXJDmD3keoDgLvBaiqB5J8EfgOvRndlzsjW2qGn7KQGrdmYa6qbwJZYdWtx3nNx4GPbyEurWBu7y0AHLzyrWOORJOqqo7QmxNCVT2TZN2fsgD+IcnSpyz+bujB6nnm/mzxlpzSjPJTFlKbLMzSDPJTFlK7LMxTbm7vLc8Pg0ngpyyk1lmYpRnipyyk9jXx7VKSRsZPWUiNszBLM8RPWUjtcyh7hnn9WZLaY2GeMhZbSZpsFmZJkhpiYZYkPc9Rt/GzME8wE0iSpo+FecZYzKXpZG5PDz8uNYHWk3yDTNC5vbd483xpQviFF5PPHrMkzRh7122zMOtFTFpJGh8LsyRJDbEwN2y9Pddh9W7tOUvS6FmYJUlqiIV5RtkTlqQ2+XGpCTCI4WwLsSRNBnvMktSorczz8GR8ctljlqTGbKSoWoCnjz1mreus3Bna0mwy90fPHrOeZ/JJ08Xbc04me8yNaf3stPX4pElkXqmfPWZJGpPNFGN7wdPPHnOjPHuWpNlkYZYkqSEOZUvSjBjESJxD6cNnj3lMpnWyx7T+u6SNmvZcmOZ/27itWZiTnJrkG0m+k+SBJO/v2k9Msj/Jw93PbV17knwmyYEk9yY5c9j/CEmSpsV6eszPAh+sqtOBs4HLk5wO7AVuq6rdwG3dc4Dzgd3dYw9w1cCjVjOmvVcgzQJzuC1rXmOuqiPAkW75mSQPAicDFwLndJtdB/wt8Htd+/VVVcAdSU5IsrP7PZoSJrI0+czjNm3oGnOSOeANwJ3Ajr5i+wSwo1s+GXi872WHurblv2tPkoUkC4uLixuNW5JmjoV0Nqx7VnaSVwE3AR+oqqeTPL+uqipJbWTHVbUP2AcwPz+/oddK2pwkpwLX0zuRLmBfVX06yYnAjcAccBC4qKp+kF6ifxq4APgR8NtVdfc4Yp82rRXZ5fG0Ft8sWVePOcnL6BXlz1XVl7rmo0l2dut3Ase69sPAqX0vP6VrkzR+zhmRGrdmj7k7Y74aeLCqPtm36mbgUuDK7udX+tqvSPIF4FeAp7y+vLpp+UygZ9eTwTkj42euaC3rGcr+VeCdwH1J7unaPkyvIH8xyWXAY8BF3bpb6Q17HaA39PXuQQYsaTC2OGfEwjxBPBmYLOuZlf1NIKusPneF7Qu4fItxSRqiQc8ZSbKH3lA3u3btGmSo0szxlpyN8IxWo3K8OSNVdWQzc0aczCkNjrfklGbIOuaMwIvnjLyru6Pf2ThnZMO8CY82yh6zNFucM6KBm5ZJrK2wMI+YZ84aJ+eMjM96ct/3B4FD2RqS473BOLQnbZ15NL0szJIkNcShbEnSpthjHw4L8wj5RyxpUg3r/cuJYy/mULYkSQ2xMEuSBsIJaYNhYdbYmMSaNv49axAszJIkNcTJX5K0RfaUNUj2mIfAIVpJ0mZZmDV0nqhI0vo5lK2hsRhL0sbZY5YkDZUn6Rtjj3kE/KPs8ThIs8Fc3xoLsyQNkEVJW2VhlqRNsghrGLzGLElSQyzMkiQ1xMIsSVJDvMY8RF5/kiRtlIVZTfLL09UyT7o1TA5lS5LUEHvMA+DZsyRpUOwxqwl+0YUk9ViYJWkNnjhqlNYszEmuSXIsyf19bR9NcjjJPd3jgr51H0pyIMlDSd4yrMAlSZpG6+kxXwuct0L7p6rqjO5xK0CS04GLgdd3r/l/krxkUMFKkjTt1izMVXU78OQ6f9+FwBeq6sdV9Q/AAeCsLcQnSdJM2co15iuS3NsNdW/r2k4GHu/b5lDXJkmS1mGzhfkq4LXAGcAR4BMb/QVJ9iRZSLKwuLi4yTDGywkhmjTOGZHat6nCXFVHq+q5qvoJ8Fl+Olx9GDi1b9NTuraVfse+qpqvqvnt27dvJgxJG3ctzhmRmrapwpxkZ9/TdwBLZ983AxcneUWS1wC7gW9tLURJg+KcEal9a975K8kNwDnASUkOAR8BzklyBlDAQeC9AFX1QJIvAt8BngUur6rnhhK5pEG6Ism7gAXgg1X1A3rzQ+7o28Y5I9IIrFmYq+qSFZqvPs72Hwc+vpWgpLXM7b3FL7gYnKuAP6R3ov2H9OaMvGcjvyDJHmAPwK5duwYdn6bA8vk45u/qvFf2OvltR5pWVXV0aTnJZ4G/6p5uaM4IsA9gfn6+hhPp+DnZU6NgYVZTlr/xeUI0fEl2VtWR7unyOSOfT/JJ4Odxzog0EhbmTfLMWZPIOSNS+yzM0gxxzojUPgvzBtlTliQNk1/7KElSQyzMkiQ1xKFsjZ2XB9Qq/zaHx2O7OnvMktTHL6cZD4/7T1mYNRFMWkmzwsIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQ/wcsyaKM7MlTTt7zJIkNcTCLEkrcHRG42JhliSpIV5jliTsIasd9pglSWqIhVnSTPM+7G2a5f8Th7IlSc2Y5YK8xB6zJEkNsTBr4jkUKWmaWJglSWqIhXkV9sIkSeNgYZYkqSEWZkmSGmJhliSpIesqzEmuSXIsyf19bScm2Z/k4e7ntq49ST6T5ECSe5OcOazgJUmaNuvtMV8LnLesbS9wW1XtBm7rngOcD+zuHnuAq7YepiRp1szqJNx1Feaquh14clnzhcB13fJ1wNv72q+vnjuAE5LsHECskrbI0a+fmtU3/Uk0a/9XW7nGvKOqjnTLTwA7uuWTgcf7tjvUtb1Akj1JFpIsLC4ubiEMSRtwLY5+SU0byOSvqiqgNviafVU1X1Xz27dvH0QYktbg6JfUvq18icXRJDur6kiXrMe69sPAqX3bndK1NW9u7y0cvPKtL2qTptxGR7+OMMHMabVuKz3mm4FLu+VLga/0tb+ruz51NvBUX9JLQzNr16GGYTOjX+ClKWmQ1vtxqRuAvwNel+RQksuAK4HfSPIw8Ovdc4BbgUeBA8Bngd8deNSSBuno0hD1Zke/vDQlDc66hrKr6pJVVp27wrYFXL6VoCSN1NLo15W8ePTriiRfAH4FR7+kkdjKNWZJE6Yb/ToHOCnJIeAj9AryF7uRsMeAi7rNbwUuoDf69SPg3SMPWOqzdKlq+VygaWNhlmbIrI1+zcobuaaLhVkTy4lekqaRX2IhSVJD7DEvYy9MkjRO9pglTT0/465JYo8Ze8mSpHbYY5YkqSEWZkmSGmJhliRNlGmfM2BhljQzpvnNXNPDwixJUkMszJoJ0z70JWl6WJg1tSzGkiaRhVmSpIbM9A1G7E1J0uSa1m8Ps8csSVJDLMySJDXEwqyp5yULSZPEwixJUkNmevKXppM9ZEmTzB6zJEkNsTBLktSQmRzKdqhTktQqe8ySJDXEwixJUkMszJIkNcTCLElSQyzMkqaGX/WpabDlWdlJDgLPAM8Bz1bVfJITgRuBOeAgcFFV/WCr+5Kk9bA4a5INqsf876rqjKqa757vBW6rqt3Abd1zSZK0hmENZV8IXNctXwe8fUj7kTQgSQ4muS/JPUkWurYTk+xP8nD3c9u445Sm3SAKcwFfT3JXkj1d246qOtItPwHsGMB+JA2fo1+aONN26WIQd/56Y1UdTvJzwP4k3+1fWVWVpJa/qCviewB27do1gDAkDcGFwDnd8nXA3wK/N65gpNUsFeeDV751zJFs3ZZ7zFV1uPt5DPgycBZwNMlOgO7nsRVet6+q5qtqfvv27VsNQ9LWbXr0K8meJAtJFhYXF0cRqzS1tlSYk7wyyauXloE3A/cDNwOXdptdCnxlK/uRNBJvrKozgfOBy5O8qX9lVRW94v0inmhLg7PVoewdwJeTLP2uz1fVXyf5e+CLSS4DHgMu2uJ+pIGapmGvQekf/UrygtGvqjqy2uiXpMHaUmGuqkeBf71C+/eBc7fyuyWNTjfi9c+q6pm+0a+P8dPRrytx9EsaiZn82kdJL+Lol9QIC7MkR7+khliYJU28afscq2abX2IhSVJDLMySJDXEoWzNFIc8JbXOHrMkSQ2xMEuS1JCZGsp2GFOSpts03NXPHrMkSQ2xMEuS1JCZGsqWNF28PKVpZI9ZkjR15vbeMrEnbhZmzbRJTl5J02kmhrJ949VapmEmp6TpYI9Z0kTyhFvTysIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSZpak3gTIQuzJEkNsTBLfSbtzFrS9LEwS5LUEAuzJEkNmYkvsZA0PbzcoGlnj1laZhJncUqaHhZmSdLUm6QT7qEV5iTnJXkoyYEke4e1H0nDZS5LozWUwpzkJcCfAucDpwOXJDl9GPuSNDzmsjR6w+oxnwUcqKpHq+p/AV8ALhzSvqShmKShryEyl6URG1ZhPhl4vO/5oa5N0mQxl6URG9vHpZLsAfZ0T3+Y5KEh7/Ik4B+HvI/NaDUuaDe2kcaVP9rQ5sOO7ReG+Ls3bcT53OrfJbQbW6txQbv5PIq4VsznYRXmw8Cpfc9P6dqeV1X7gH1D2v+LJFmoqvlR7W+9Wo0L2o2t1big7dg2ac1chtHmc8vHuNXYWo0L2o1tnHENayj774HdSV6T5OXAxcDNQ9qXpOExl6URG0qPuaqeTXIF8DXgJcA1VfXAMPYlaXjMZWn0hnaNuapuBW4d1u/fhJENm29Qq3FBu7G1Ghe0HdummMsb0mpsrcYF7cY2trhSVePatyRJWsZbckqS1JCpK8xr3T4wySuS3NitvzPJ3AhiOjXJN5J8J8kDSd6/wjbnJHkqyT3d4w+GHVffvg8mua/b78IK65PkM90xuzfJmSOI6XV9x+KeJE8n+cCybUZ2zJJck+RYkvv72k5Msj/Jw93Pbau89tJum4eTXDqsGKdNi7nc7bfZfG4xl7v9NpPPE5HLVTU1D3qTUx4BTgNeDnwbOH3ZNr8L/Fm3fDFw4wji2gmc2S2/GvjeCnGdA/zVmI7bQeCk46y/APgqEOBs4M4x/L8+AfzCuI4Z8CbgTOD+vrY/BvZ2y3uBP1rhdScCj3Y/t3XL28bx/zxJj1ZzudtXs/ncei73/d+OLZ8nIZenrce8ntsHXghc1y3/BXBukgwzqKo6UlV3d8vPAA8yWXdPuhC4vnruAE5IsnOE+z8XeKSqHhvhPl+gqm4HnlzW3P+3dB3w9hVe+hZgf1U9WVU/APYD5w0rzinSZC7DxOfzuHMZxpzPk5DL01aY13P7wOe3qapngaeAnx1JdEA33PYG4M4VVv+bJN9O8tUkrx9VTEABX09yV3p3cFpu3LdlvBi4YZV14zpmADuq6ki3/ASwY4Vtxn3sJlXzuQxN5nPruQxt5nNTuTy2W3LOoiSvAm4CPlBVTy9bfTe9oZ0fJrkA+Etg94hCe2NVHU7yc8D+JN/tzirHLr2bWrwN+NAKq8d5zF6gqiqJH3GYIY3mc7O5DJORzy3k8rT1mNdz+8Dnt0nyUuBngO8PO7AkL6OXxJ+rqi8tX19VT1fVD7vlW4GXJTlp2HF1+zvc/TwGfJneMGK/dd2WcUjOB+6uqqPLV4zzmHWOLg0Ddj+PrbDNOI/dJGs2l7v9NZnPjecytJvPTeXytBXm9dw+8GZgaTbdbwH/tbor+8PSXfe6Gniwqj65yjb/cun6WJKz6P3fjOKE4ZVJXr20DLwZuH/ZZjcD7+pmdJ4NPNU37DNsl7DKsNe4jlmf/r+lS4GvrLDN14A3J9nWzfR8c9em42syl6HdfJ6AXIZ287mtXB7WzLdxPejNOvwevRmdv9+1fQx4W7f8z4E/Bw4A3wJOG0FMb6R37ede4J7ucQHwPuB93TZXAA/Qm316B/BvR3S8Tuv2+e1u/0vHrD+2AH/aHdP7gPkRxfZKeon5M31tYzlm9N5MjgD/RO/a0mX0rmfeBjwM/A1wYrftPPBf+l77nu7v7QDw7lEcu2l4tJjL3X6bzOeWc7nbdxP5PAm57J2/JElqyLQNZUuSNNEszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDXjruAABOOumkmpubG3cYUvPuuuuuf6yq7eOO43jMZ2l9VsvnJgrz3NwcCwsL4w5Dal6Sx7b4+lOB64EdQAH7qurTSU4EbgTmgIPARVX1gyQBPg1cAPwI+O2quvt4+zCfpfVZLZ8dypZmy7PAB6vqdOBs4PIkpwN7gduqajdwW/cc4Hxgd/fYA1w1+pCl2WJhlmZIVR1Z6vFW1TPAg8DJwIXAdd1m1wFv75YvBK6vnjuAE5LsHG3U0myxMEszKskc8AbgTmBHVR3pVj1Bb6gbekX78b6XHeraJA2JhVmaQUleBdwEfKCqnu5fV1VF7/rzRn7fniQLSRYWFxcHGKk0e9YszElOTfKNJN9J8kCS93ftJybZn+Th7ue2rj1JPpPkQJJ7k5w57H+EpPVL8jJ6RflzVfWlrvno0hB19/NY134YOLXv5ad0bS9QVfuqar6q5rdvb3rSuNS89fSYnSwiTYlulvXVwINV9cm+VTcDl3bLlwJf6Wt/V3fCfTbwVN+Qt6QhWPPjUl0SHumWn0nSP1nknG6z64C/BX6PvskiwB1JTkiy02SWmvCrwDuB+5Lc07V9GLgS+GKSy4DHgIu6dbfS+6jUAXofl3r3SKOVZtCGPse8xckiFmZpzKrqm0BWWX3uCtsXcPlQg5L0Auue/OVkkbXN7b2Fub23rLtdkmaF74Prt67C7GSRtvnHLknTYz2zsp0sMiAWUEnSWtZzjdnJIpIkjch6ZmU7WUSSpBHxzl+SJDXEwixJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMyN8HZ1kiSwMI/NuAqxJwCS1DYLsyRJDbEwS5LUEAuzNEOSXJPkWJL7+9puTHJP9zi4dE/8JHNJ/mffuj8bW+Bak5eppsd6vsRC0vS4FvjPwPVLDVX175eWk3wCeKpv+0eq6oxRBaeVLRXcg1e+dcyRaBQszBNk2Mlp8k+/qro9ydxK67qveL0I+D9GGpQmlu8Zw+FQ9gxz6EvL/BpwtKoe7mt7TZL/nuS/Jfm1cQWmzTHHJ9OaPeYk1wC/CRyrql/q2m4EXtdtcgLw/1fVGd2Z+IPAQ926O6rqfYMOepxaOENsIQZNpUuAG/qeHwF2VdX3k/wy8JdJXl9VTy9/YZI9wB6AXbt2jSTYaTBJuTxJsU669fSYrwXO62+oqn9fVWd0155uAr7Ut/qRpXXTVpRbM86zYc/Ep0uSlwL/J3DjUltV/biqvt8t3wU8AvziSq+vqn1VNV9V89u3bx9FyDNpszm3Vr6ay21ZszBX1e3Akyut67smdcNK6yVNjF8HvltVh5YakmxP8pJu+TRgN/DomOIbOYuVxmWrk79WvSYFPA38p6r6/7a4j6nWcvI7dDV9ktwAnAOclOQQ8JGquhq4mBefYL8J+FiSfwJ+AryvqlY8Sdf6rSevzL3ZttXC7DUpaYJU1SWrtP/2Cm030btUpSnVcsdglm26MPddk/rlpbaq+jHw4275riRL16QWlr++qvYB+wDm5+drs3FofUxAaePsuWoctvJxKa9JzRAne0mD1Z9Pk5hfkxbvJFmzMHfXpP4OeF2SQ0ku61atdk3q3u6Wfn+B16RGZhITW5L0YmsOZXtNSpJmj8P44+OdvyRpgCZ19GpS455G3itbA9XiWfbc3luaikezadBFb7Vca6m4tvh+MAkszMvMyh9SS8krTTNzbWNm5T34eBzKbozDSZI02+wx92mtIE7SmeOwj92wjsUkHWNJs8Ee8ya1VsRXYw9cOj5zRK2xx6wtGdYbmj1ZzaLN5NOknVQsz+3V4p/lSZv2mCVJaoiFWSPn0KE0HubeZHAoW00nasuxSRquWb2kZWEegmn+Y7JQSlrNsN/7pvm9tZ+FeQvW+iPZShGbhgK4mSSahn93y5JcA/wmcKyqfqlr+yjwO8Bit9mHq+rWbt2HgMuA54D/UFVfG3nQ0iZMchG3MEuz5VrgPwPXL2v/VFX9SX9DktPpfYvc64GfB/4myS9W1XOjCFRt8uR5+CzMGhkTevyq6vYkc+vc/ELgC1X1Y+AfkhwAzqL3NbDS88ztwXJWtiSAK5Lcm+SaJNu6tpOBx/u2OdS1SRqiNQtzl6jHktzf1/bRJIeT3NM9Luhb96EkB5I8lOQtwwpc0sBcBbwWOAM4Anxio78gyZ4kC0kWFhcX137BlJrGnuMkfMSq9fg2aj1D2dfiNanjmrY/iuWm/d8366rq6NJyks8Cf9U9PQyc2rfpKV3bSr9jH7APYH5+voYT6cZN8gQg/dSsvQetWZi9JiVNtyQ7q+pI9/QdwNLo2M3A55N8kt6J9m7gW2MIcShm7c1+vSb9uEx6/LC1yV9XJHkXsAB8sKp+QO/60x1923hNasSm4Y9Sw5PkBuAc4KQkh4CPAOckOQMo4CDwXoCqeiDJF4HvAM8Cl0/76FdrzOf1m6ZjtdnCfBXwh/QS+Q/pXZN6z0Z+QZI9wB6AXbt2bTIMSRtRVZes0Hz1cbb/OPDx4UXUnuO9wTs0rlHYVGGe5mtSwzZNZ3XH47dOadha+vahWcnrQfKYrW5ThXlWr0nNinEkzKj2aWGfbv7/ahqsWZi9JtXj2d3wbPbY+iaslvmeoc1az6zsmbwm5Zv+5PONUcPm35iGwVtyrsHEa996TqL8f9RW+PejUbIwa+hWe1Mb9Judb55azTSNgPl3Pv0szJIm1moF1+KlSWZhlqQNsvBrmPx2KUkzw4KqSWCPWZKwaE+rSZxfYI9ZkqSGWJglSRNlEr4jeiscytZMmuakljTZLMySpIk0rSfYDmVLfaY10aVRmvah5mGzMEszJMk1SY4lub+v7f9O8t0k9yb5cpITuva5JP8zyT3d48/GFrg0QyzM0my5FjhvWdt+4Jeq6l8B3wM+1Lfukao6o3u8b0QxSjPNwizNkKq6HXhyWdvXq+rZ7ukdwCkjD0zS89YszA59qUVevxqa9wBf7Xv+miT/Pcl/S/Jr4wpKmiXr6TFfi0Nf0tRL8vvAs8DnuqYjwK6qegPwH4HPJ/kXq7x2T5KFJAuLi4ujCViaUmt+XKqqbk8yt6zt631P7wB+a8BxSRqhJL8N/CZwblUVQFX9GPhxt3xXkkeAXwQWlr++qvYB+wDm5+drRGFviqMto+cx35hBXGN26EuaYEnOA/4v4G1V9aO+9u1JXtItnwbsBh4dT5TS7NjSDUaOM/T1/SS/DPxlktdX1dMrvHYPsAdg165dWwlD0joluQE4BzgpySHgI/QuRb0C2J8E4I7uMtSbgI8l+SfgJ8D7qurJFX+xpIHZdGGepaEvTb9ZGWqrqktWaL56lW1vAm4abkQah1n5e59UmxrKduhLkqThWLPH7NCXJEmjs55Z2Q59SZI0It75S5KkhliYJUlqiIVZkqSGWJglSWrIlm4wIkkt8HO5mib2mCVJashMFua5vbd4hi1JatJMFmZJklplYZYkqSEWZknS1JukS5gzXZgn6T9KkjQbZrowS5LUGguzJEkNsTBLMyTJNUmOJbm/r+3EJPuTPNz93Na1J8lnkhxIcm+SM8cXuTQ7LMzSbLkWOG9Z217gtqraDdzWPQc4H9jdPfYAV40oRmmmraswe5YtTYequh14clnzhcB13fJ1wNv72q+vnjuAE5LsHEmg0gxbb4/5WjzLlqbVjqo60i0/Aezolk8GHu/b7lDX9iJJ9iRZSLKwuLg4vEilGbCuwuxZtjQbqqqA2sTr9lXVfFXNb9++fQiRSbNjK98utdGz7CM0ys8ya8YdTbKzqo50J9HHuvbDwKl9253StUkaooFM/trMWbZDX1IzbgYu7ZYvBb7S1/6ubt7I2cBTfSfjkoZkK4X56NIQ9WbOsh36kkYvyQ3A3wGvS3IoyWXAlcBvJHkY+PXuOcCtwKPAAeCzwO+OIWRp5mxlKHvpLPtKXnyWfUWSLwC/gmfZUjOq6pJVVp27wrYFXD7ciCQtt67C3J1lnwOclOQQ8BF6BfmL3Rn3Y8BF3ea3AhfQO8v+EfDuAccsSdLUWldh9ixbs2RpMuDBK9865kgkzSLv/CVJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMySJDXEwixJUkO2cucvSRo5v3RG084esyRJDbEwS5LUEAuzJEkNsTBLktQQC7MkSQ1xVrYkkrwOuLGv6TTgD4ATgN8BFrv2D1fVraONTpotmy7MJrI0ParqIeAMgCQvAQ4DX6b3feqfqqo/GV900mzZdGE2kaWpdS7wSFU9lmTcsUgzZ1DXmJ9P5AH9PknjczFwQ9/zK5Lcm+SaJNtWekGSPUkWkiwsLi6utImkdRpUYd5wIktqT5KXA28D/rxrugp4Lb3RsSPAJ1Z6XVXtq6r5qprfvn37KEKVptaWC/NmE9kzbKlJ5wN3V9VRgKo6WlXPVdVPgM8CZ401OmkGDKLHvKlE9gxbatIl9I1+JdnZt+4dwP0jj0iaMYP4uNSLErmqjnRPTeQB8Kb9GoUkrwR+A3hvX/MfJzkDKODgsnWShmBLhdlElqZHVf0P4GeXtb1zTOFIM2tLhdlEliRpsLwlpyRJDbEwS5LUEAuzJEkNsTBLktSQmfp2KT92JElqnT1mSZIaYmGWJKkhFmZJkhpiYZYkqSEWZkmSGmJhliSpIRZmSZIaYmGWJKkhFmZJkhoyU3f+krS6JAeBZ4DngGeraj7JicCNwBy971e/qKp+MK4YpVmw5R5zkoNJ7ktyT5KFru3EJPuTPNz93Lb1UCWNwL+rqjOqar57vhe4rap2A7d1zyUN0aCGsk1maTpdCFzXLV8HvH18oUizYVjXmE1mafIU8PUkdyXZ07XtqKoj3fITwI7xhCbNjkFcY15K5gL+36rah8ksTaI3VtXhJD8H7E/y3f6VVVVdnr9IV8j3AOzatWv4kUpTbBA95jdW1ZnA+cDlSd7Uv7Kqil7xfoEke5IsJFlYXFwcQBiStqKqDnc/jwFfBs4CjibZCdD9PLbKa/dV1XxVzW/fvn1UIUtTacuFebPJbCJL7UjyyiSvXloG3gzcD9wMXNptdinwlfFEKM2OLRVmk1maGjuAbyb5NvAt4Jaq+mvgSuA3kjwM/Hr3XNIQbfUa8w7gy0mWftfnq+qvk/w98MUklwGPARdtcT+ShqiqHgX+9Qrt3wfOHX1E0uzaUmE2mSVJGixvySlJUkMszJIkNcTCLElSQyzMkiQ1xMIsSVJDLMySJDXEwixJUkMszNIq5vbewtzeW8YdhqQZY2GWJKkhFmZJkhpiYZYkqSEWZkmSGmJhliSpIRZmSZIaYmGWRJJTk3wjyXeSPJDk/V37R5McTnJP97hg3LFK027ThdlElqbKs8AHq+p04Gzg8iSnd+s+VVVndI9bxxeiNBteuoXXLiXy3UleDdyVZH+37lNV9SdbD0/SKFTVEeBIt/xMkgeBk8cblTSbNt1jrqojVXV3t/wMYCJLUyDJHPAG4M6u6Yok9ya5Jsm28UUmzYaBXGPeTCIn2ZNkIcnC4uLiIMKQtEVJXgXcBHygqp4GrgJeC5xBr0f9iVVeZz5LA7LlwrzZRK6qfVU1X1Xz27dv32oYkrYoycvo5fLnqupLAFV1tKqeq6qfAJ8FzlrpteazNDhbKsxbSWRJ7UgS4Grgwar6ZF/7zr7N3gHcP+rYpFmz6clfx0vkbiIJmMjSpPhV4J3AfUnu6do+DFyS5AyggIPAe8cRnDRLtjIr20SWpkRVfRPICqv8eJQ0YpsuzCayJEmD552/JElqiIVZkqSGWJglSWqIhVmSpIZYmCVJaoiFWZKkhliYJUlqyFZuMDIx5vbeMu4QJElaF3vMkiQ1xMIsSVJDLMySJDXEwixJUkMszJIkNcTCLElSQ4ZWmJOcl+ShJAeS7B3WfiQNl7ksjdZQCnOSlwB/CpwPnA5ckuT0YexL0vCYy9LoDavHfBZwoKoerar/BXwBuHBI+5I0POayNGLDKswnA4/3PT/UtUmaLOayNGJjuyVnkj3Anu7pD5M8NORdngT845D3sRmtxgXtxjbSuPJHG9p82LH9whB/96aNOJ9b/buEdmNrNS5oN59HEdeK+TyswnwYOLXv+Sld2/Oqah+wb0j7f5EkC1U1P6r9rVercUG7sbUaF7Qd2yatmcsw2nxu+Ri3GlurcUG7sY0zrmENZf89sDvJa5K8HLgYuHlI+5I0POayNGJD6TFX1bNJrgC+BrwEuKaqHhjGviQNj7ksjd7QrjFX1a3ArcP6/ZswsmHzDWo1Lmg3tlbjgrZj2xRzeUNaja3VuKDd2MYWV6pqXPuWJEnLeEtOSZIaMnWFea3bByZ5RZIbu/V3JpkbQUynJvlGku8keSDJ+1fY5pwkTyW5p3v8wbDj6tv3wST3dftdWGF9knymO2b3JjlzBDG9ru9Y3JPk6SQfWLbNyI5ZkmuSHEtyf1/biUn2J3m4+7ltldde2m3zcJJLhxXjtGkxl7v9NpvPLeZyt99m8nkicrmqpuZBb3LKI8BpwMuBbwOnL9vmd4E/65YvBm4cQVw7gTO75VcD31shrnOAvxrTcTsInHSc9RcAXwUCnA3cOYb/1yeAXxjXMQPeBJwJ3N/X9sfA3m55L/BHK7zuRODR7ue2bnnbOP6fJ+nRai53+2o2n1vP5b7/27Hl8yTk8rT1mNdz+8ALgeu65b8Azk2SYQZVVUeq6u5u+RngQSbr7kkXAtdXzx3ACUl2jnD/5wKPVNVjI9znC1TV7cCTy5r7/5auA96+wkvfAuyvqier6gfAfuC8YcU5RZrMZZj4fB53LsOY83kScnnaCvN6bh/4/DZV9SzwFPCzI4kO6Ibb3gDcucLqf5Pk20m+muT1o4oJKODrSe5K7w5Oy437towXAzessm5cxwxgR1Ud6ZafAHassM24j92kaj6Xocl8bj2Xoc18biqXx3ZLzlmU5FXATcAHqurpZavvpje088MkFwB/CeweUWhvrKrDSX4O2J/ku91Z5dild1OLtwEfWmH1OI/ZC1RVJfEjDjOk0XxuNpdhMvK5hVyeth7zem4f+Pw2SV4K/Azw/WEHluRl9JL4c1X1peXrq+rpqvpht3wr8LIkJw07rm5/h7ufx4Av0xtG7Leu2zIOyfnA3VV1dPmKcR6zztGlYcDu57EVthnnsZtkzeZyt78m87nxXIZ287mpXJ62wrye2wfeDCzNpvst4L9Wd2V/WLrrXlcDD1bVJ1fZ5l8uXR9Lcha9/5tRnDC8Msmrl5aBNwP3L9vsZuBd3YzOs4Gn+oZ9hu0SVhn2Gtcx69P/t3Qp8JUVtvka8OYk27qZnm/u2nR8TeYytJvPE5DL0G4+t5XLw5r5Nq4HvVmH36M3o/P3u7aPAW/rlv858OfAAeBbwGkjiOmN9K793Avc0z0uAN4HvK/b5grgAXqzT+8A/u2Ijtdp3T6/3e1/6Zj1xxbgT7tjeh8wP6LYXkkvMX+mr20sx4zem8kR4J/oXVu6jN71zNuAh4G/AU7stp0H/kvfa9/T/b0dAN49imM3DY8Wc7nbb5P53HIud/tuIp8nIZe985ckSQ2ZtqFsSZImmoVZkqSGWJglSWqIhVmSpIZYmCVJaoiFWZKkhliYJUlqiIVZkqSG/G8PcwZO43IALAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 576x576 with 4 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(2,2, figsize=(8,8))\n",
    "ax = ax.flatten()\n",
    "for i in range(4):\n",
    "    ax[i].hist(sili[:,i], 100, range=[-1,11])\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
Valerio Mascagna's avatar
Valerio Mascagna committed
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABBQAAAFpCAYAAAAsm3skAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABRR0lEQVR4nO39f9RdZX3g/b8/JcWZtigEE8g3IQVJrCSIETNgVllZ0BgKtoIWaomsMWlSGTvyfLVM52t8OmtZpg+P6FpoZcaxxQeWlCVBH6dTojKxAWS5yjJKpJhCaCYp0JUwaQgEqFStCf18/zj7Dic3949z7vNj7332+7XWve5z9tnnnOucz7X3vvbnXNe1IzORJEmSJEnqxs+UXQBJkiRJklQ/JhQkSZIkSVLXTChIkiRJkqSumVCQJEmSJEldM6EgSZIkSZK6ZkJBkiRJkiR1zYSCJEmSJEnqmgkFSZIkSZLUNRMKkvoqIs6KiD+JiK9GxO+WXR5JkiRJg2FCQdK0IuK2iHgmIh4dt/ySiNgVEXsiYiNAZj6emR8E3gv8chnllSRJkjR4kZlll4HXv/71efrpp5ddDHXg+9///rOZOacfrxUR7wLedcIJJ3zgjW98Yz9eUgPywx/+kOOOO47HH3/85cycBRARxwH/C1gN7AMeAtZk5s6IuAz4XeCOzLxzutd3H1AP/dz+2xn/+vAY0GyD2Ae4/deHx4BmM/6arA7MKqMwY8YaE4sWLWL79u1lFkUdioi/79drZebXgK8tX778A8a/+p566inOOOOMw22LzgP2ZOYTABFxF3A5sDMzNwObI+IbwIQJhYi4BrgGYOHChe4DaqCf23+7008/3fjXhMeAZutn/G0D1k+/jwHWgXox/pqsDpQ65CEzv5aZ17zuda8rsxiSZmY+sLft/j5gfkRcGBE3R8SfAvdM9uTMvCUzl2fm8jlz+p7wliRVmG1AWQeazfiPjlJ7KEgaPZn5APBAJ+u2Z6fVPMa/2Yy/JEn156SMKk1EvCsibnnxxRfLLopm5mngtLb7C4plHTM73WzGv9mMf7PZBpB1oNmM/+goNaFgRWo2G5O19xCwOCLOiIjjgauAzSWXSZJUA7YBZB1oNuM/OpxDoRMRZZdAZTL+rFmzhhUrVgC8JiL2RcSGzDwCXAt8E3gc+EpmPtbN63acVIwwDk1m7EfSlNu/MdewWec0xrrQHMa6LxzyIGlamzZtYv/+/QAPZ+aCzLwVIDPvycw3ZuaZmXlDt69bm6SipL5z+5dUGk8k1SnryrRMKFhJpNI47KnZjL8kSVK9mVBQaTyZUN9/oTRBWCv+Qt1sHgOazfjLOtBsxn90OCmjSuPJhLo2VcLAZIJUKx4Dms34yzrQbMZ/dDgpo6T6M5kgSZIkDZ1DHiSVZka9lKZLHphcGC3Gs7mMvSRJlWdCQVJp7KUkSaoFE1ySNCETCpIkaeicR0mSpPozodAvZq67ZmNSM64D7dub215tuQ9oNnsoSaol2x3SMbzKg0ozUo1JDy4z0nMd8HuvtZHaB0jqim1AWQeazfiPDq/yIKmevIRksxhTaaTYBpR1oNlGJv62TxzyoBHnRi5JkiQ1i+cAQ2NCQZIkSZIkdc2EgqTSOH5Oai63f0mS6s+EgqTSDHT8nF3dpEobmfGzkqrF4780VCYUJElSNXliIElSpY1+QsHGiCQNTUScFRF/EhFfjYjfLbs8kiRJGpzRTyho6DyhkEbL+vXrmTt3LsDS9uURcUlE7IqIPRGxESAzH8/MDwLvBX55+KWVVCbbAM1m/JvN+DdT3xMKVqTR5AmF1Fzr1q1jy5YtxyyLiOOAzwGXAkuANRGxpHjsMuAbwD09vbE9zKRKqFUboKz9xgjvr2oVf/Wd8dd0OkooRMRtEfFMRDw6brkVqSFKO6HQSHOW93pYuXIls2fPHr/4PGBPZj6RmT8F7gIuB8jMzZl5KXD1ZK8ZEddExPaI2H7w4MGJVuhX8SX1yDZAsxn/ZjP+mk6nPRS+CFzSvsCK1Cz9PqGY9mRCjdDRLO+9nFiOf64nqf00H9jbdn8fMD8iLoyImyPiT5niGJCZt2Tm8sxcPmfOnEGXVVIPbAM0WylJZVWG8dd0ZnWyUmZ+OyJOH7f4aEUCiIixirQzMzcDmyPiG8CdE71mRFwDXAOwcOHCmZVeZZvohOL8iLgQ+A3gNUxyQpGZtwC3ACxfvjwHWkpJQ5OZDwAPdLJuRLwLeNeiRYsGWSQNUUScBXwYeD1wX2Z+vuQiaXBsAzTbjOMP1oERYPx1VEcJhUlYkTShTk8oPJlQKSIg3eX0wdPAaW33FxTLOpaZXwO+tnz58g/0s2Dqr/Xr1/P1r399bAztURFxCfBZ4Djg/8nMGzPzceCDEfEzwJ8BJhQaxjZAs5lUbjbj30x9n5QxMx/IzP9vZv67zPzcVOs6frr2ej6haBS725fH734QHgIWR8QZEXE8cBWwuZsX8BhQD46f1SR6agN0NORNVdaXpPLI1oHRb3cYfx3VS0LBk0n1dELhjkRDNfoH94FZs2YNK1asAHhNROyLiA2ZeQS4Fvgm8Djwlcx8rJvX7XofYAxL4fhZTaKnNoAJxdozqdxsxr8bI95+6SWh0HNF8oSyPgZ1QiEN3IjvxIdh06ZN7N+/H+DhzFyQmbcCZOY9mfnGzDwzM2/o9nUb1ZgYPU7K2SCDaAPYBqyPyiSVVQrjr+l0NIdCRGwCLgReHxH7gI9n5q0RMVaRjgNu82RydG3atAmAiHg4M5ePLc/Me5hht1bHTknN5hwKo8fxs6PJNkCzDSL+xetZB2rA+Gs6HfVQyMw1mTkvM3/WX6fUL2Ym5T5Aqi3Hz6onxl/WgWYz/qOj75MydsOK1GyeTKq0fYDDICrBfUCtOX5WkiSVm1CwMdFsJpSkZnMfUA9j42d37doFcM7Qx8+aABxJtgFlHWg24z867KEgSZImNTYp5+HDhwF2OOxR/WAbUNaBZjP+o6PUhIKazcak1GzuA5rNxqQkSfXnkAeVxsak1GzuA5rNNkCzGX9ZB5rN+I8OhzxIaibHZUulsg3QbMZf1oFmM/6jwyEPkqRS+OuEJElSvZlQUGk8mZCazV8nms1jQLMZf1kHms34jw7nUFBpPJlQ6Rz2UD/GbGR4DGg24y/rQLMZ/9HhHAqSJEnSoJgIlTTCZpVdAEmjJSLeDfwa8Frg1sz8y3JLJEmSJGkQnENB0rTWr1/P3LlzAZa2L4+ISyJiV0TsiYiNAJn5F5n5AeCDwG/19Mb+qjPSHPbWbMZfkqT6M6Gg0vTcmPRkc2jWrVvHli1bjlkWEccBnwMuBZYAayJiSdsq/6l4vHvGthEc9tZsxr/ZTCjJOtBsPcXfdmKlOCmjSmNjcpwK7xxXrlzJ7Nmzxy8+D9iTmU9k5k+Bu4DLo+WTwP/MzIcne82IuCYitkfE9oMHDw6u8NOJqPR3L0mjyDaArAPNZvxHh5MyVoknNaqX+cDetvv7imX/B/AO4MqI+OBkT87MWzJzeWYunzNnzmBLKkmSVBW2+TVejeuEkzJK6qvMvBm4uZN1I+JdwLsWLVo02EJJqpyutv8IyBx4mSRJUnecQ0HSTD0NnNZ2f0GxrGP2UpKay+1fkqT6M6EgaaYeAhZHxBkRcTxwFbC5mxeo5DwqNe5yJknSyPM4LVWKCQWVppInk5rQmjVrWLFiBcBrImJfRGzIzCPAtcA3gceBr2TmY928rr9QNpv7AEmSpHrzKg8qjSeT9bFp0yb2798P8HBmLsjMWwEy857MfGNmnpmZN5RbStWN+wCpuWwDyjrQbMZ/dHiVB0mlqdzBxG6UkjQUtgFlHWg24z866j3kwca/VGuVPpi4f5EkSZKmVO+EgiQNgskEaeAq10NJkiR1zYSCpNJU6oTCJII0VJXuoSRJkjpiQkFSaTyhkCR1xeSvJFWKCQVJ1WSjUVI79wmSJFXOrLILIEkaHRHxbuDXgNcCt2bmX5ZbIkmSJA3KQHooRMS7I+ILEfHliLh4EO+h6jL+6lSl5lDQpNavX8/cuXMBlrYvj4hLImJXROyJiI0AmfkXmfkB4IPAbw2/tJLKZBugSyPW88b4N5vxb6aOEwoRcVtEPBMRj45bboOyATyh0CA4h0I9rFu3ji1bthyzLCKOAz4HXAosAdZExJK2Vf5T8bikmrMN0GzGv9mMv6bTTQ+FLwKXtC+wQVmCkjLZnlBIzbVy5Upmz549fvF5wJ7MfCIzfwrcBVweLZ8E/mdmPjzsskrqP9sAzWb8m834azodJxQy89vAoXGLbVA2RL9PKCLimojYHhHbDx48ONjCSxqE+cDetvv7imX/B/AO4MqI+OBkT3YfMLrs8jp6bAM02yCSytaB+jD+mk6vcyjMuEFpRRoJM45/Zt6Smcszc/mcOXMGX1JVknMojJ7MvDkz35aZH8zMP5livVuA64GHjz/++O7eZMTGHNeBXV41AdsAzdZTUtk6UHvGX0cN5CoPmXkzcPM069wSEfuBdx1//PFvG0Q5VI5O4g+tk0ngXYsWLRp8oVRJmfk14GvLly//QNllmVAEZJZdiqp6Gjit7f6CYlnHKh9/HbVu3TquvfZa3va2Vw7XbV1eV9NqTD4UEZszc2exil1eG8g2QLN1Gn+wDowi499MvfZQ6KlB6YRstWf8peZ6CFgcEWdExPHAVcDmbl7AHir1YZd3TcA2QLP1JalsHagt46+jek0o9NSgtDFZe8Z/2OzqXY6Gf+9r1qxhxYoVAK+JiH0RsSEzjwDXAt8EHge+kpmPdfO6NiZqb/hd3hu+LVaMbYBmM6ncbMZfR3Vz2chNwHeAX+pXg9LGZH0M4oTC+KtWGnwis2nTJvbv3w/wcGYuyMxbATLznsx8Y2aemZk3dPu6NiZGU6fzaMw4/g3eFstiG6DZTCo3m/HXdDqeQyEz10yy/B7gnpm8uWNn6mPTpk0ARMTDmbl8bLnxlzRTzqFQez13eVU92AZotkHEv3g968BEKjZ/k/HXdHod8tATM1PNZvxVi1+oJ/o11F9IJeixy6vHgGYz/rIONJvxn4GKtj9LTSio2WpxMqmB8mDSbO4D6mNQXV4lSVK9lZpQsDHZbJ5MSs3W0z6goln6UTWIeTRsAzSb8Zd1oNmM/+hwyIMkdcqTWKlvbAM0m/GXdaDZjP/ocMiDSmNmUrVngqEn7gOazfhLklR/DnlQacxMSs3mPqDZeo6/Cb1asw0o60CzGf/R4ZAHSeqVJzbScI1tc257tWUbUNaBZjP+o8MhD5LUC09opBnx1ylJkurPhIJKY2NSta4DJhJ6Vuv4q2f+OtVsbv+yDjSb8R8dzqGg0tiYVGPqgMmHCTmGXn1jXaidxuz/NSnrQLMZ/9HhHAqSJGno/FFBkqT6c8iDJPWLv5JKHfNHBUmS6s+EgiT1U4SJBUmSJDWCcyioNMZfaraR2weMTyRNd19qsJHb/tU160CzGf/R4RwKKo3xVy11c1Job4UpVX4fMF3sOo2tdWCw/H5rqfLbf1WMcP22DjSb8R8dDnmQJKnfpuqdYKIJGMCvU36nkiQNnQkFSRqWmZ7weKI0PBN91932ShnEuiNoIL9ONfw7lSRp2EwoSFI/9HIi40lQPTjEQZIk6RgmFCRJajfRkITJ7veaPOjX60iSJJXAhIKkvoqIN0TErRHx1bLLUqpBzvA/IieftZzhudshEV7pQZKk7nisrBUvGzmKarIRGv8BGUD8169fz9y5cwGWHvtWcUlE7IqIPRGxESAzn8jMDX0vhEbOSM3wXJP9bmMYj8qzDSDrQA9GYB9n/EeHl41UaYx/faxbt44tW7YcsywijgM+B1wKLAHWRMSSEopXL+Nn+x//2DAaCSPQEJlUnT+bwx+Gx++4dLYBZB1oNuM/OhzyIGlaK1euZPbs2eMXnwfsKXok/BS4C7i809eMiGsiYntEbD948GAfS1sD/TiZ8YRodDQ0lqX9OtXQ71uSNOJKOr6ZUJA0U/OBvW339wHzI+LkiPgT4K0R8bHJnpyZt2Tm8sxcPmfOnEGXVaNqWL061HdD+3WqvY5YVyRJ6qtZZRdA0mjJzOeAD5ZdDmlaYyeXmeWWQ5JmIsL9l6TS2UNB0kw9DZzWdn9BsaxjTsjTppPLEvrrancm+74mu1JDFb7fKpRBkiSpQyYUJM3UQ8DiiDgjIo4HrgI2d/MCTsijvhqlk/GqJDjqrJPvz+9YkqSemFCQNK01a9awYsUKgNdExL6I2JCZR4BrgW8CjwNfyczHunldeygM0GQ9HiRJkqQ+GUhCISLeEBG3RsRXB/H6qjbjP3o2bdrE/v37AR7OzAWZeStAZt6TmW/MzDMz84ZuX9ceCqPH7X8ATAZ1x++rVO4Dmq0v8Xcbri23/2bqOKEQEbdFxDMR8ei45ZdExK6I2BMRGwGKy8ht6HdhVZ7169czd+5cgKXty42/NPpqv/1P1zid6vHJ5lvQhGxMjqba7wPUE+PfbMZf0+mmh8IXgUvaF0TEccDngEuBJcCaiFjSt9KpMtatW8eWLVuOWWb81SuHPEyj1wkZ+3Ti6/bfbGONybPPPvuY5TYmm8N9QLMZ/2Yz/ppOxwmFzPw2cGjc4vOAPUUD4qfAXcDlnbxeRFwTEdsjYvvBgwc7LrDKsXLlSmbPnj1+8eDi7y+AjeCQh3ro9/YPJRwD6rpPqUC5bUxq6G0AVcpIHAM0Y8Zf0+l1DoX5wN62+/uA+RFxckT8CfDWiPjYRE/MzFsyc3lmLp8zZ06PxVBJjL/UXDPe/sF9QJ0MojGpkWAboNk8BjSb8ddRswbxopn5HPDB6daLiHcB71q0aNFED0Jm/wungetL/NUIja4D3Y7bn2y9fu4n+/B6nW7/rbcbQvxH9WoXY5+jvOPkRI3J8yPiZOAGisZkZn5ioidHxDXANQALFy4cdFmnZnujr2wDNFvljgEaKuPfTL32UHgaOK3t/oJiWUfs7lx7xl89sQ7UWk/bPxj/UZSZz2XmB4srv0yYTCjWuwW4Hnj4+OOPH07hRiWZVB22AZrNY0CzGX8d1WtC4SFgcUScERHHA1cBmzt9shOy1Z7xl5qrp+0fStgH1P2EslrltzEp2wDNVr9jgPrJ+Ouobi4buQn4DvBLEbEvIjZk5hHgWuCbwOPAVzLzsU5f08ZEfaxZs4YVK1YAvMb4S80yiO0f3AfUXD0ak90MH9KkbAM0m8eAZjP+mk5kieMG28bOfGD37t3jH5x+TGM/1hnW+3SiSp950qfG9zNz+Yye/OrXKj/+nahSXEr+XvoZ/+L1Jq4DNu471x7H9nH148fYt8d8svH309SLocW/vYydav/ME92fbv2qGx/TIc6hsGbNGh544AGeffZZjhw5chj43cy8NSLeCfwxcBxwW2beMJPXX758eW7fvv3YhcOIzUR1RlMaWhtg6icOpo3T6+tO9fhMH+v1uX021DbAoGLRyTrDaK/1GtcS5n4ZehugLm3mupWlB5PVgV6HPPTEzFSzGX9ZB5rN+NfDpk2b2L9/P4cPHwbYkZm3AmTmPZn5xmK+hK6TCZXs7lqnJFPNuf3LOtBsxn90lJpQULNVsjEpaWgGug8YtRPDyT5PxCt/NWNjUpKk+is1oeAJZbPZmJSazX2A1Fy2AWUdaDbjPzoc8iCpNB5M+mCiX6in+jV7svs1/IVb9Vb69m+dL5VtQFkHms34jw6HPEgqjQeTZiv9hFKlcvuXJKn+HPKg0hh/qdk8oZSayzaArAPNZvxHh0MeVBrjL0nNVbnG5GRDgJo4NGIIn9k2gKwDzWb8R4dDHiRJ0tDZmJQkqf4c8qDSGH+p2fq6D+j2F9VR/NV5FD+TRtaU2791uRFsBzab8R8dDnlQaYy/PJgMyVSN806vCDEA7gOarbLb/1TbxEy2C0+OJ+T2L+tAsxn/0eGQB0ml8WAiNZfbvyRJ9WdCQZIkSZIkdc2EgiRJkiRJ6pqTMqo0xl9qNvcBMzTdmHzH7KsG3P5lHRigGhwHjP/ocFJGlcb4S83mPqDZatuYLHEi01Hi9i/rQLMZ/9HhkAdJkjR0NiYlSao/EwqSJEmSJKlrJhQklaa2XZ7rZHw37E66Zdt1W5IkSR0woSCpNHZ5liRJkurLhIIkSZIkSeqal41UaYy/1GzuAwao12ErQxj2Uov4z/R7qOuVICKGVsZaxF8DZR1oNuNfkgHs471spEpj/KVmcx/QbMa/2Yy/rAPNZvxHh0MeJEmSJElS10woSJIkSZKkrplQkCRJkiRJXTOhIEmSJEmSumZCQZIkSZIkdc2EgiRJkiRJ6tqsfr9gRPw88N+AnwIPZOaX+v0eqjbrQLMZ/2Yz/s1m/JvN+Ms60GzGv5k66qEQEbdFxDMR8ei45ZdExK6I2BMRG4vFvwF8NTM/AFzW5/KqJOvXr2fu3LkAS9uXWweawfg3m/FvtrH4n3322ccsN/7NYPzlMaDZjL+m0+mQhy8Cl7QviIjjgM8BlwJLgDURsQRYAOwtVnu5P8VU2datW8eWLVuOWWYdaA7j32zGv9mMf7MZf1kHms34azodJRQy89vAoXGLzwP2ZOYTmflT4C7gcmAfrco05etHxDURsT0ith88eLD7kmuwIo65u3LlSmbPnj1+rRnXAeNfcQOOf+strAOV0R7vcbEH418rE8Tv6LKIV/660Nj4T7VdtH+X7d/p+O+3fXm/ytLPdTvQ2Pg3wWR1xTaA2uqA8W+gLo8jvUzKOJ9XMlDQqkDzgT8HroiIzwNfm+zJmXlLZi7PzOVz5szpoRgq0YzrgPEfCe4Dms34N5vxbzbjL+tAsxl/HdX3SRkz85+A3+5k3Yh4F/CuRYsW9bsYKlGndcD4jyb3Ac1m/JvN+Deb8Zd1oNmMfzP10kPhaeC0tvsLimUdy8yvZeY1r3vd63oohkrUUx0w/rXnPqDZjH+zGf9mM/6yDjSb8ddRvSQUHgIWR8QZEXE8cBWwuZsXiIh3RcQtL774Yg/FUIl6qgPGv/Z63geo1jwGNJvxbzbjL+tAsxl/HdXpZSM3Ad8Bfiki9kXEhsw8AlwLfBN4HPhKZj7WzZubmaqPNWvWsGLFCoDX9KsOGP/6GET8wYNJXQwq/u4D6mEs/rt27QI4x/g3i/GXx4BmM/6aTmRmeW/+ytiZD+zevXv8gzBd2fqxzrDepxM1+MwR8f3MXD71kztTifh3ogZxGVZZ+hn/dsuXL8/t27cf+/4q37g60O/4T7sPUG/G4je2Pbd/p3U4BtTdRN/5TI9J3Tyvl+PeNHWkEm2AfhxLB/G6My3zoF63E5M9f0htgEnrwCBj0ck6w2iv9RrXQW0HUzx/6G2AmreZR7Esk9WBUhMKRwsRcRD4+3GLXw88W0JxZqop5f3FzOzrdKzGvzQzKXPf4w8T1oE6fp/DUPb3Mqz4Q/mfdZDq/Nk8BjS7bE2PfxXLNcwy2QZ4RdXLOIjy2QaYXFPKO2EdqERCYSIRsX0Qv4QOiuXtr6qXb7y6lReqXeYql61MTfpeRvmzjvJn65cqf0eWbfCq+jmqWK4qlqlXdfhMVS9j1cs3nbqVv+nl7WVSRkmSJEmS1FAmFCRJkiRJUteqnFC4pewCdMny9lfVyzde3coL1S5zlctWpiZ9L6P8WUf5s/VLlb8jyzZ4Vf0cVSxXFcvUqzp8pqqXserlm07dyt/o8lZ2DgVJkiRJklRdVe6hIEmSJEmSKsqEgiRJkiRJ6lolEwoR8VRE/E1EPBIR28suz3gRcVtEPBMRj7Ytmx0RWyNid/H/pDLL2G6S8v5hRDxdfMePRMQ7yyxjO+PfXzWM/yURsSsi9kTExrLLM0gRcVpEfCsidkbEYxHx4WL5hPUpWm4uvpsdEXFu22utLdbfHRFry/pMvapr/CfabzU5jjNVVvyrFr9ujjMzKU9EvK34vHuK58ZMy9pvVWgDVPU4X7fj+UxUIf7jVbU+TFO+WtaLKsZ/vKrXh/GGUj8ys3J/wFPA68suxxTlWwmcCzzatuxTwMbi9kbgk2WXc5ry/iHw+2WXzfgb/3FlPQ74O+ANwPHAD4AlZZdrgJ93HnBucfsE4H8BSyarT8A7gf8JBPB24LvF8tnAE8X/k4rbJ5X9+ZoU/4n2W02NYx3jX7X4dXOcmUl5gO8V60bx3EvLjv9UsSihDJU8zk9Srj+kgsfzOse/LvVhFOtFFeNft/pQRv2oZA+FqsvMbwOHxi2+HLi9uH078O5hlmkqk5RXM2T8B+o8YE9mPpGZPwXuovXdjqTM3J+ZDxe3fwg8Dsxn8vp0OfBn2bINODEi5gG/CmzNzEOZ+TywFbhkeJ+kb0Yt/k2N40xVLf6lxa/L40xX5Skee21mbstWy/LPqNAxqwqqepyv2fF8ZFS1PoyxXgxX1evDeMOoH1VNKCTwlxHx/Yi4puzCdOiUzNxf3P4H4JQyC9Oha4vukbdVqWsOxn9Yqhj/+cDetvv7imUjLyJOB94KfJfJ69Nk38+ofG91/hwT7beaGseZKvPz1yF+/SrP/OL2oMrZq6q2Aap8nK/i8Xymqhr/8apcH8bUsV7UJf7j1aE+jNe3+lHVhMIFmXkucCnwoYhYWXaBulFk/Kt+Pc7PA2cCy4D9wE2lluZYxn/wqhz/xomIXwD+O/CRzPzH9sdqUp80zX7LOFZereJXtfL0WeXbABX7/kfteF75+I9Xsfowpq71onbxH6+i9WG8vtaPSiYUMvPp4v8zwP+g1Q2y6g4U3Qgp/j9TcnmmlJkHMvPlzPwX4AtU6Ds2/oNX4fg/DZzWdn9BsWxkRcTP0komfCkz/7xYPFl9muz7GZXvrbafY5L9VlPjOFOlff6axK9f5Xm6uD2ocvakwm2ASh7nK3w8n5EKx3+8StaHMXWtFzWK/3iVrg/j9bt+VC6hEBE/HxEnjN0GLgYenfpZlbAZGJtBeS1wd4llmdZYpS+8h4p8x8Z/OKoaf+AhYHFEnBERxwNX0fpuR1JEBHAr8Hhmfrrtocnq02bg/dHyduDFoovdN4GLI+KkotvaxcWyuqll/KfYbzU1jjNVSvxrFL++lKd47B8j4u3FPuj9VOSYVfE2QCWP8xU+nnet4vEfr5L1YUwd60XN4j9epevDeH2vH5PN1ljWH63ZnX9Q/D0G/EHZZZqgjJtodQ85TGvs4QbgZOA+YDdwLzC77HJOU947gL8BdtDaCOaVXU7jb/yL8r6T1tUO/q6K8e/zZ72AVre4HcAjxd87J6tPtGZk/1zx3fwNsLzttdYDe4q/3y77szUp/pPtt5ocxzrFv4rx6+Y4M5PyAMtpNSD/DvivQJQd+6liUUI5Knmcn6RclT2e1zX+dakPo1Yvqhr/utWHMupHFG8kSZIkSZLUscoNeZAkSZIkSdVnQkGSJEmSJHXNhIIkSZIkSeqaCQVJkiRJktQ1EwqSJEmSJKlrJhQkSZIkSVLXTChIkiRJkqSumVCQJEmSJEldM6EgSZIkSZK6ZkJBkiRJkiR1zYSCJEmSJEnqmgkFSZIkSZLUNRMKkiRJkiSpa7PKLgDA61//+jz99NPLLoY68P3vf//ZzJzTz9c0/vUxiPiDdaAujL88BjSb8W82jwHNZvw1WR2oRELh9NNPZ/v27WUXQx2IiL/v92sa//oYRPzBOlAXxl8eA5rN+Debx4BmM/6arA445EGSJEmSJHXNhIIkSQ23fv165s6dy9lnn3102X/8j/+RN73pTZxzzjm85z3v4YUXXjj6WER8LCL2RMSuiPjVtuWXFMv2RMTGoX4I9d3evXu56KKLWLJkCUuXLgWYCxARfxgRT0fEI8XfO8eeM1ndkCSNJhMKkiQ13Lp169iyZcsxy1avXs2jjz7Kjh07eOMb38gnPvEJACJiCXAVsBS4BPhvEXFcRBwHfA64FFgCrCnWVU3NmjWLm266iZ07d7Jt2zaAuW0x/UxmLiv+7oHJ60YphZckDYUJBUmSGm7lypXMnj37mGUXX3wxs2a1plp6+9vfzr59+8Yeuhy4KzP/OTOfBPYA5xV/ezLzicz8KXBXsa5qat68eZx77rkAnHDCCQA/BuZP8ZTJ6oYkaUSZUJAkSVO67bbbuPTSS8fuzgf2tj28r1g22fJXiYhrImJ7RGw/ePDgAEqsfnvqqacAfg74brHo2ojYERG3RcRJxbKO6oDxl6TRUY+Ewp1RdglUFdYFjbEuaCasN1274YYbmDVrFldffXXfXjMzb8nM5Zm5fM6cvl+FTL0at5289NJLXHHFFQB7M/Mfgc8DZwLLgP3ATd28/LTxdzttNuMvsB7USD0SCpIkaei++MUv8vWvf50vfelLRBxt3D0NnNa22oJi2WTLVWOHDx/miiuuGEsovQCQmQcy8+XM/BfgC7wyrME6IKl3JhNqxYSCJEl6lS1btvCpT32KzZs383M/93PtD20GroqI10TEGcBi4HvAQ8DiiDgjIo6nNTnf5qEXXH2TmWzYsIGzzjqL66677ujyiJjXttp7gEeL25PVDUnSiJpVdgEkSVK51qxZwwMPPMCzzz7LggULuP766/nEJz7BP//zP7N69WqgNTEjQGY+FhFfAXYCR4APZebLABFxLfBN4Djgtsx8rIzPo/548MEHueOOO3jzm9/MsmXLAJYUl4hcExHLgASeAv4dTF03JEmjyYSCqu/OgPdl2aWQRt7pp5/OCSecwHHHHcesWbPYvn07hw4d4rd+67fGJmRbHBEnZebz0er//lngncCPgHWZ+TBARKwF/lPxsv9XZt5ewsdRFzZt2vSqZRs2bHjVsj/90z8FIDNvAG4Y/3hx+cB7ei6Q+/1KuOCCC8h8JQ4RsXO6GE9WNyRJo8khD5Kko771rW/xyCOPsH37dgBuvPFGVq1axe7duwF+CGwsVr2UVnfmxcA1tCZpIyJmAx8Hzqc1rvrjbTPAS5I0PcfQS7VhQkGSNKm7776btWvXjt19Dnh3cfty4M+yZRtwYjGu+leBrZl5KDOfB7YClwy52JIkSRoCEwqSJAAigosvvpi3ve1t3HLLLQAcOHCAefOOzr92GDiluD3Z9eY7ug69JKna9u7dy0UXXcSSJUtYunQpwFxo9USLiK0Rsbv4f1KxPCLi5ojYExE7IuLcngsxUU8Fey9IleIcCpIkAP7qr/6K+fPn88wzz7B69Wre9KY3TbRaXwa2R8Q1tIZKsHDhwn68pCSpj2bNmsVNN93Eueeeyw9/+ENe+9rXzo2IJcA64L7MvDEiNtIaCvdRjh0Kdz6toXDnl1R8jTrn2qkMeyhIkgCYP7/VkWDu3Lm85z3v4Xvf+x6nnHIK+/fvH1vlZ4FnituTXW++o+vQZ+Ytmbk8M5fPmTOnvx9EktSzefPmce65rU4GJ5xwAsCPafU4uxwYm2z3dqYfCtc7eyWUor2XCrA0Ij4MQ+6losqbNqEQEadFxLciYmdEPFZqRXJnIg3M+vXrmTt3LmefffbRZYcOHWL16tUsXrx47NJxx8HU23lErC32C7uL2f5VA//0T//ED3/4w6O3//Iv/5Kzzz6byy67jNtvP3qRhpOBu4vbm4H3F3Xh7cCLmbmf1iUDL46Ik4rjwsXFMklSTRVX+vk54LvAKcX+HuAfmH4o3KtExDURsT0ith88eHAgZVbvxnqp7Ny5E+Bx4ENFL5WNtHqpLAbuY5oJmzXaOumhcAT4D5m5BHg7ViRpJK1bt44tW7Ycs6x9hv9Vq1YBnFo85Az/I+bAgQNccMEFvOUtb+G8887j137t17jkkkvYuHEjW7duZfHixQCvBW4snnIP8ASwB/gC8O8BMvMQ8EfAQ8Xffy6WSZJq6KWXXuKKK64A2JuZ/9j+WLauK9p1v3N7qdVDey8V4F9oJRXK6aWiypp2DoUiA7m/uP3DiGivSBcWq90OPEBr/NTRigRsi4gTI2JeWyazc46NkYZm5cqVY79AHHX33XfzwAMPALB27Vo+9rGPjSUHJtzOae0Tto6dQEbE2Az/r77IvSrlDW94Az/4wQ9etfzkk0/mvvvuAyAi/tdYbIvYf2ii18rM24DbBldaSdIwHD58mCuuuIKrr76ahx9++IVi8YGxtn1x7J9uKJxGw/HAW5lZL5VjzgN7mkfJ88PK6WoOhYg4nZlXpPGvZVcnqeLaZ/g/9dRT4ZUkZM8z/LsPkCSpujKTDRs2cNZZZ3Hddde1P7QZGBvSuJbph8Kp5l566SWAM4GP9KOXij1URkvHCYWI+AXgvzOsiuR8CVKlRPR3m/RgIkmalO3A0j344IPccccd3H///SxbtgxgSUS8k9bQt9URsRt4B9MMhVO9jfVSAQ5l5p8Xiw+MDWWwl4o6SihExM/SSiZ8yYrULJW4BrFK0z7Df/H/SPFQTzP8S5KkarvgggvITHbs2MEjjzwCsDMz78nM5zJzVWYuzsx3tA+Fy8wPZeaZmfnmzNxe6gdQz9p7qQAH2h6yl4qO6uQqDwHcCjyemZ9ue8iK1ADts7tu27YNYOwaxE7K2QDtM/wX/18oHnKGf0mSpBHW3kuFVg+VR+ylovGmnZQR+GXg3wJ/ExGPFMv+T1oV5ysRsQH4e+C9xWP3AO+kVZF+BPx2Pwus4Zo3b97RMfQTXIP4wmK1wUzKqaFas2YNDzzwAM8++ywLFizg+uuvZ+PGjbz3ve/l1ltv5Rd/8RfhlUl1JtzOM/NQRIzN8A/O8C9JkobFCfv6aqyXCkBE7MzM5W0Prxq//lQTNg+F8S9FJ1d5+CtgsoFs1atIGpger0Hcv9ld27nj6JtNmya+EMPYDP8AEfEyOMO/JEmSpC6v8lAr4yfzcXKfnvT7GsQdTchnzCRpKNavX8/cuXM5++yzjy47dOgQq1evZvHixaxevZrnn38emHqunIhYW8ytszsi1r76nSRJ6pLnBJU2ugkF9U37NYh5ZQz9YCflnG7H4Y5Fkvpm3bp1bNmy5ZhlN954I6tWrWL37t2sWrWKG28cGyI78Vw5ETEb+DhwPnAe8PGxCXslSSqd5w8DYUJBU6rFNYjdOUhST1auXMns2bOPWXb33Xezdm1rN7927Vr+4i/+Yuyho3PlZOY24MQisfyrwNbMPJSZzwNbgUuG9BEkSaPM9n5lmVDQlLwGsSQ104EDB45Oynvqqady4MDRK4ZNNlfOZMtfJSKuiYjtEbH94MGD/S66RoknEWo3WX2wnkil6eQqD2qw9tld4egMr/cUd52UU5IaICJoXUW6PzLzFuAWgOXLlzuzriRJNWUPBUmS9CqnnHIK+/e3Rqzt37+fuXPnjj002Vw5/ZlDR5Wxd+9eLrroIpYsWcLSpUsB5kJrvoyI2FpMvrl1bK6MqSbslCSNJhMKkiTpVS677DJuv/12AG6//XYuv/zysYcmmyvnm8DFEXFScYJ5cbFMNTVr1ixuuukmdu7cybZt2wDmRsQSYCNwX2YuBu4r7sMkE3ZKM+IwBqkWTChIkprBxumk1qxZw4oVK9i1axcLFizg1ltvZePGjWzdupXFixdz7733snHj2DnjxHPlZOYh4I+Ah4q//1wsU03NmzePc89tdTI44YQTAH5Ma16My4Hbi9VuB95d3J5swk5pMPqxX/fYUH/GsFT1m0PhzoD3OdxSkqR+2bRp04TL77vvvlctm2qunMy8Dbitn2VTNTz11FMAPwd8Fzil7QpO/wCcUtyebGLOY672FBHX0OrBwMKFCwdWZkkN5jnj0IxGDwWzUpIkSQPx0ksvccUVVwDszcx/bH+sSDB11WrPzFsyc3lmLp8zZ04fSypJGrbRSCioGWaSODLZJEnSjB0+fJgrrriCq6++GuCFYvGBsaEMxf9niuXDmZjTY7skVUY9Ewp3xuAPJh6sRovxlCSpK5nJhg0bOOuss7juuuvaH9oMrC1urwXubls+0YSdUv/ZtpMqoZ4JBUlS37388su89a1v5dd//dcBePLJJzn//PNZtGgRv/VbvwUQABHxmoj4cnFpuO9GxOljrxERHyuW74qIXy3jc0jqjwcffJA77riD+++/n2XLlgEsiYh3AjcCqyNiN/CO4j5MMmGnJA2ESaVKMKEAVkZJAj772c9y1llnHb3/0Y9+lN/7vd9jz549nHTSSQCvLx7aADyfmYuAzwCfBCguJ3cVsBS4BPhvEXHcED+CpD664IILyEx27NjBI488ArAzM+/JzOcyc1VmLs7Md4xdzaO4usOHMvPMzHxzZm4v9QNotPTSXret3yzGe6hMKGh0uPOQZmzfvn184xvf4Hd+53eAVlfn+++/nyuvvBKAtWvXApxYrN5+ybivAqsiIorld2XmP2fmk7R+pTxvaB9CkqR+sV0pdcSEgkaTBwGpKx/5yEf41Kc+xc/8TOuw8Nxzz3HiiScya1br6sILFiwAOL5Y/eil4TLzCPAicDKTXzLuVSLimojYHhHbDx482P8PJEm2BUbDoOJo/WgOYz1Q1U8oTFUBJnqsmwpj5ZI68pnPfIalS5cCLI2ITRHxryLijGL8/J5iPP3xMPX4elXT17/+debOncvb3va2ob2nl42TJElT8gpvtVD9hEKnrDyjzfiW5umnn+bmm29m+/btAI8Bx9EaJ/9J4DPFOPrnaY2rh0nG16u6HnzwQTZv3szpp5/OVVddxf3338+HP/xhXnjhBY4cOQK0hkQAPy2ecvTScBExC3gd8BzDumScJIFtA0mqgNFJKEB/DiwenKRXOXLkCD/+8Y/H7v4csB/4FVrj56E1nv7dxe3Jxteroj7xiU+wb98+nnrqKe666y5+5Vd+hS996UtcdNFFfPWrrRDffvvt8Mo16NsvGXclcH9mZrH8qqKXyhnAYuB7w/wskiRJU/J8r69GK6Egqe/mz5/P7//+77Nw4UKAt9AaL/994IVi/DwcO1Z+svH1x3AMffV98pOf5NOf/jSLFi3iueeeA3i2eOhW4OSI2ANcB2wEyMzHgK8AO4EtwIcy8+Xhl1zSyPEEQDNl3ZEGalbZBSiNOxepI88//zx33303Tz75JHPnzt0B/DytSwL2JDNvAW4BWL58efb6euqPCy+8kAsvvBCAN7zhDXzve690MIiIBMjMnwC/OdHzM/MG4IaBF1RS89wZ8L5pDhedrCOpPjxnq7z691CwkkkDde+993LGGWdQTJyXwJ8DvwycWIyfh2PHyk82vl6SJNXI+vXrmTt3LmefffbRZRHxhxHxdEQ8Uvy9s+2xjxWTMu+KiF8tpdCeG6gT1pO+qX9CYTr9rCxWPDXQwoUL2bZtGz/60Y/GFq2i1aX9W7TGz0NrPP3dxe3JxtdLkqQaWbduHVu2bJnooc9k5rLi7x6AiFhCa9LmpbR6Mv63iDhueKVVv1UuoTR2LuY5WaVMm1CIiNsi4pmIeLRtWbUzk+NZ6aQZO//887nyyis599xzodVI+BlaQxU+ClxXjKM/mda4ephkfL0kSUNhu69vVq5cyezZsztd/XLgrsz858x8EtgDnDewwk2n20vP61VMKKkTnfRQ+CITj5e2IqlZGnzwuf766/nbv/1bgMcy898WjYUnMvO8zFyUmb+Zmf8MrfH1xf1FxeNPlFt6Sb34zGc+w9KlS8d+oTojIv5VRJwREd8tfkD4ckQcD1Bc4ePLxfLvRsTppRZe0qBcGxE7ih8eTyqWHZ2UudA+YfMxKj0xc4Pbe+PVOqE0E8Z+RqZNKGTmt4FDHb5e9SuSFWX0GWNJ6ounn36am2++me3bt/Poo48CBK0fDj5J64eFRcDzwIbiKRuA54vlnynWkzRaPg+cCSyjdRnpm7p9gcy8JTOXZ+byYo4m1UtPCaWh8rxg4HqZQ2F0M5PTsWI2i/GW1GBHjhzhxz/+MUeOHIFWu2E/8CvAV4tVbgfeXdy+vLhP8fiqiHAnqpnx+FtJmXkgM1/OzH8BvsArPx4enZS50D5hc/V1Wt+slz0nlGp9HqhXmWlCoTmZSXca1ZuQpV+MrSRNaf78+fz+7/8+CxcuZN68eQAvA98HXsjMI8Vq7T8eHP1hoXj8RVpzrBzDxqS60snxuttjum2AGYuIeW133wOMzbO2GbiqGPp0BrAY+N7459ea9aYvCaW+ngcak9LNKKEwsplJTWjkJ2RxRyRJE3r++ee5++67efLJJ/nf//t/Q6vdMNG8Sl2pzY8KUsOtWbOGFStWsGvXLoBzImID8KmI+JuI2AFcBPweQGY+BnyF1pWgtgAfysyXSyq6BmRkE0qeD8zYrOlXebWImJeZ+4u74yvSnRHxaeD/Q9kVyWx1X6xcuZKnnnqq09WPzqMBPFnM9H8e8J0BFU+SNCD33nsvZ5xxBm0n/S8AvwycGBGzil4I7T8ejP2wsC8iZgGvA54bbqkl9cumTZuO3o6IHZl5K69c1elVMvMG4IYhFE1DsGbNGh544AGeffZZeCWhdGFELAMSeAr4d9BKKEXEWELpCKOSULoz4H1e/Xwq0yYUImITcCHw+ojYB3ycplUkTebaiHg/sB34D5n5PK3urtva1plwHo2IuAa4BmDhwoVDKCruECSpSwsXLmTbtm386Ec/4l//638NcAKtY/y3gCuBu4C1wN3FUzYX979TPH5/ZrrjlTR4/jDYdyaU1IlpEwqZuWaCxVYkfR74I1pJpT+iNY/G+k6fnJm3ALcALF++vLqNzckOTiYnJDXA+eefz5VXXsm5557LrFmzoHWVh1uAbwB3RcT/Bfw1r7QLbgXuKHqnHaI1BE6SJI2oGQ15kDLzwNjtiPgC8PXirvNoSNIIuf7667n++usBiIgniyFtTzDBZaEz8yfAbw63hGoUZ+JXJ4y/NDS9XDZyuMrcMbhTepVaT8hiPCVJGg2DuAKEpOZw/9Cz+iQUVBpn+JUkqXlG6rLRnjQIpq8H4x8fu2/9kSY1mgkFN/q+2rRpE/v37+fw4cMAOzLz1sz8t5n55sw8JzMva7vqB5l5Q2aemZm/lJn/s+8Fqmp8q1ouSZJmYOQvGy1ptNk2H4rRTCiUod8V1g2gdzP9Dv3uJUli5cqVzJ49u9PVj142OjOfBMYuG11tHvMFDq2WemBCYdjcaUiSpHq7NiJ2RMRtEXFSsWw+sLdtnQkvGw2tS0dHxPaI2H7w4MGZlaCX9tRMnmv7rbqmik37Y8ZQ4011NTd1zISCqm3YG7Q7EEmSpvJ54ExgGbCf1mWju5KZt2Tm8sxcPmfOnD4Xb8BsJ0jNdWc4EewEmplQaFiQG2XYv1hII+InP/kJ5513Hm95y1tYunQpH//4xwF48sknOf/881m0aBHAGyLieIDiSi5fLiZg+25EnD72WpWemE1STzLzQGa+nJn/AnyBV4Y1eNloDY9tNnXLOjMwzUwoaDQMMnkwWQbSnZFG1Gte8xruv/9+fvCDH/DII4+wZcsWtm3bxkc/+lF+7/d+jz179gAcATYUT9kAPJ+Zi4DPAJ8EJ2aTRl2tLxs9nclm+JdUfZ1ur2W170d4f2JCoVcjXDlqpapxqGq5pHEigl/4hV8A4PDhwxw+fJiI4P777+fKK68cW+054N3F7cuB24vbXwVWRURQ14nZJL1K+2WjFyxYAPB6RvWy0R6vm2NYsbZOqSFMKEjDVNODywsvvDB2Urk0Ih6PiBURMTsitkbE7uL/SQDRcnPR5X1HRJxbbunVqZdffplly5Yxd+5cVq9ezZlnnsmJJ57IrFmzxlb5Ka9MsnZ0ArbMPAK8CJxMhxOz9WVSNkkD1X7Z6H379gE8W+plo6UqqWmbTuo3Ewr95I6l+oZ5ec8Rqg8f/vCHueSSSwAeA94CPA5sBO7LzMXAfcV9gEtpdXVdDFxDawIv1cBxxx3HI488wr59+/je977H3/7t3w7svWo9KZskSd3yihPN0qAYm1CQNKUXX3yRb3/722zY0Bo6n5k/zcwXOLbL++0c2xX+z7JlG3DiuDG3qrgTTzyRiy66iO985zu88MILHDlyZOyh43llkrWjE7BFxCzgdbSGRDgxmyRpdDToxFCaCRMKVeIOazSMWByffPJJ5syZw2//9m8DLImI/ycifh44pa2r6z8ApxS37fJeQwcPHuSFF14A4Mc//jFbt27lrLPO4qKLLuKrX/3q2GonA3cXtzcDa4vbVwL3Z2YyChOzSaqXETvuqsYGXRet68Pl990REwpqtjIm5qnZLNJHjhzh4Ycf5nd/93ehNdnWP/HK8AYAihPJ7OZ17fJeLfv37+eiiy7inHPO4d/8m3/D6tWr+fVf/3U++clP8ulPf3rsspGzgFuLp9wKnBwRe4DrKOpE7Sdmk9QcFT/+qmb6VZ+sl/XTyxUmRsCs6VdRT+4MeN8051mdrCOVZMGCBSxYsIDzzz9/bNFXaZ08HoiIeZm5vxjS8EzxuF3ea+icc87hr//6r1+1/A1veAPf+16rg0FEPJGZ/wyQmT8BfnOi18rMG4AbBldaSaoQ23Eab0RPHFXo1zY/IvuOZvVQ6OfGPd1ruSPRiDj11FM57bTT2LVr19iiVbR+fW7v8r6WY7vCv7+42sPbgRfbZwGXVC9jV3l505veBK0rvXiVF0kaBM8fmmcEYt6shILUDxNt+COwM5jKf/kv/4Wrr74aYAmwDPi/gRuB1RGxG3hHcR/gHuAJYA/wBeDfD7u8kvpn7CovxVU/duJVXlRVg/ixZ8SP75rEsLqwe+WH0dJLDGt85TgTCoPghCwaMcuWLWP79u0AOzPz3Zn5fGY+l5mrMnNxZr4jMw9Baz6FzPxQcS3yN2fm9nJLL2mmxl/lhdYm/gJe5UWamG00DZt1TiUzoTARJ1VRmaw3kiqi/Sovb33rWwF+0au8aKTdGbWbPFkl66Z+WJfUiZnWk5LqlwmFKnJnI0mqgParvBSTdv4LXuVFmp5tOWm0jOI23afPZEKhTkaxIjeZ8ZRUcRNc5eV54FyKq7wAeJUXaXStX7+euXPncvbZZx9dNlKTsg66V3LZbb2y31+NYEJBkiRNaIKrvLwWr/IizUwNT+7WrVvHli1bxi92UlY1SxW23QpP4DltQiEibouIZyLi0bZlo5OZHNOv2VyHNSusOtfJd133eNS9/JIqa+wqL+eccw7Av8arvEivNqJXgFq5ciWzZ88ev7iZk7KOQDy7NfI9VPqhgfVivE56KHwRuGTcMjOTM2GFkyTVzNhVXnbs2AHwd17lRWq8niZlBSdmrYtG91Cp+nlbhco3bUIhM78NHBq3uJmZyfGGEcix9yix0pidHJIK7RgkSWqMbnsXeLw+aiaTshbPc2LWGrCHyoAMax8ypPeZ6RwKZiYbpNHZyTE2HiRJqp8mHL+H/xmdlLUbNbsEYAdG+zywn8PXqxvDyc2gzD1Pylj7zGQdAz1kZicrzDk7JEnqXD+Ph809tjopq4AROA/sVYUnShymmSYUzEz2Q68Vr9yTyZ6yk5XOTFbdoHZYDd4RSpI0I/1OUIz9VcSaNWtYsWLF2JVezomIDTgpa/dDYnrppVCxOoHngd3pJXZlxb3L951pQsHMZLdGbKxMu5lkJ0cmMwn1j221DlKSJFVfQ46dmzZtYv/+/Rw+fBhgR2be6qSsDTB1/fY8sBfDnINvSDq5bOQm4DvAL0XEvkZlJofRjaWqk/5M/95mJwdpfEZ6EHWhIY0hSdLMOTFziTxOq1t9rjP2UFEnOrnKw5rMnJeZP5uZC8xMzkAVDwi9l8ns5LAMa2hMk/idSFJHnJi5Aibr8t7PrvXSBOyhMgM1GabQTz1PyiiqsfMeYBnMTkqS1ExOzNzG3oIalJnUg9G7eoQGaYBxN6HQFD1UIrOTk3CmaI2IvXv3ctFFF7FkyRKWLl3KZz/7WQAOHTrE6tWrWbx4MatXrwY4Dqbu1hwRa4tu0LsjYu2Ebyip7pyYeSLdHsuHeRKp0VfHyf80EkwolMGNVg3sDqXqmjVrFjfddBM7d+5k27ZtfO5zn2Pnzp3ceOONrFq1it27d7Nq1SqAU4unTNitOSJmAx8HzgfOAz4+NrZa0mhq/MTMM+Eky5ImM3477nVOvyEMjTKhUAfD6NJUvUvSNNN0MTBGGoB58+Zx7rmtTgYnnHACZ511Fk8//TR33303a9e2OhkU/8eSA5N1a/5VYGtmHsrM54GtwCVD/jiSBs+JmfvN47umU/erimlkmVCom6myVp0+1x3F4FT7yhydrTvBYy+//DJvfetbARYBRMQZEfHdosv7lyPi+GL5a4r7e4rHT+/mI6h8Tz31FH/913/N+eefz4EDB5g3rzX8+dRTTwWYVaw2Wbfmjro7wwh3eZaawYmZuzHRL4y2xdSLbn6AmqjOWf/URyYUqsYNfDSM2DVmP/vZz3LWWWe1L/ok8JnMXAQ8D2wolm8Ani+Wf6ZYTzXx0ksvccUVV/DHf/zHvPa1rz3msYj+1rfGd3mWasKJmXtUx3ZdHctcJ36/Gq+XYQn9qk89vI4JBamOBn0wanv9ffv28Y1vfIPf+Z3fAVoT8gG/Any1WGX8DN9jM39/FVgV/T4T1UAcPnyYK664gquvvprf+I3fAOCUU05h//7Wj4vF/yPF6pN1a7a7szRinJh5BqpyOUdPXJvJS4xqyEwo6FjuXDTORz7yET71qU/xMz9zdHdxMvBCZo6dXLZ3az/a5b14/MVi/Vexy3t1ZCYbNmzgrLPO4rrrrju6/LLLLuP221v5oeL/C8VDk3Vr/iZwcUScVEzGeHGxTDXnsCepQmyrjYYqXPbRulRvFblSjAmFUeTOQX3y9a9/nblz5/K2t72t769tl/fqePDBB7njjju4//77WbZsGcuWLeOee+5h48aNbN26lcWLF3PvvfcCjI2FnrBbc/Er5R8BDxV//3nsl0vVm8OepAEq46TA8fSjxRhqTAn7k1nTr6Ja6qRijK3zvpx4uertznh1bLv04IMPsnnzZu655x5+8pOfAJwAfJbWrP6zil4I7d3ax7q874uIWcDrgOd6KoQG7oILLqB15bdXu++++47ejoiX4ehl4j400fqZeRtwW/9LqbKMDXv6gz/4AzZt2tQ+7Ol9xSq3A39I6/Khlxe3oTXs6b9GRORkFUzS4PWhPaAhsy2uGrGHQlW441AFfeITn2Dfvn089dRT3HXXXQA/zMyrgW8BVxarjZ/he2zm7yuB+z2RkOptEMOeHPIklWSy9qbtUGlw6jJMZYavbUJhlHgwqKZhXiLqzhhWPfgocF1E7KF1snBrsfxW4ORi+XXAxmEURtJgDGrYk0OepBJ4yUpNZXhtSFVJH2LukIdR546h2gYVnwG87oUXXgitMfNk5hPAeePXycyfAL/Z9zcfBLuAStNy2JNUsmFfOk5SeQY1UeeA9w/2UJAkSRNy2JM0Ikw4aDLWDfXIhILUFFMdMOwGKak7DnuS6m6iY77tgHrpNl7GtzoGHYshxtohD6PAnYMkacBGbtiTVGX9OtkvuSu0pNFnD4W6qstsoZIkSRoe23WShsiEgiRJktR0JiIkzYAJBanuemkA2HiQJGn0eHyXNCQmFCRJkiRJUtdMKDSNGWtJkqTRYLtOUslMKEiSJEmS1BR9TEZ62Uipafw1Q5Ik9SgingJ+CLwMHMnM5RExG/gycDrwFPDezHy+rDJKajOgc4CeeihExFMR8TcR8UhEbC+WzY6IrRGxu/h/Un+KKkmSJKlCLsrMZZm5vLi/EbgvMxcD9xX3VYYB/4DkeaDG9GPIgzuSBnNnIkmSpMLlwO3F7duBd5dXFA2B54EayBwK7kiax52JJElSsyTwlxHx/Yi4plh2SmbuL27/A3DKRE+MiGsiYntEbD948OAwyqrh8DywgXpNKLgj0UTcmUiSNMLsoSjggsw8F7gU+FBErGx/MDOT1rnCq2TmLZm5PDOXz5kzZwhF1QB4Hiig94SCOxLNaGfijmRInIBRkjQ49lBssMx8uvj/DPA/gPOAAxExD6D4/0x5JdSAeR4ooMeEgjsSMcOdiTsSqVrWr1/P3LlzOfvss1sL7gwOHTrE6tWrWbx4MatXrwY4DiBabo6IPRGxIyLOHXudiFhb/Dq5OyLWlvFZJJXGHooNERE/HxEnjN0GLgYeBTYDY/v+tcDd5ZRQg+Z5oMbMOKHgjkTgzqQJ9u7dy0UXXcSSJUsAlkbEh2Hyrq1TnWyqutatW8eWLVuOWXbjjTeyatUqdu/ezapVqwBOLR66FFhc/F0DfB5adQL4OHA+rX3Bx+3yLI0suzs32ynAX0XED4DvAd/IzC3AjcDqiNgNvKO4rxHjeaDa9dJDwR1Jw7kzaYZZs2Zx0003sXPnToDHafVEWcLkXVsnPNlUta1cuZLZs2cfs+zuu+9m7drWplz8H0sOXA78WbZsA04skoe/CmzNzEPFdce3ApcM6SOoiRzWVSa7OzdYZj6RmW8p/pZm5g3F8ucyc1VmLs7Md2TmobLLqoHwPFBHzZrpEzPzCeAtEyx/DljVS6FUG6cA/yMioFWX7szMLRHxEPCViNgA/D3w3hLLqB7NmzePefPmjd39F1pJhfm0TiovLJbfDjwAfJS2k01gW0ScGBHz2n61Uk0cOHDgaOxPPfVUeOWYMR/Y27bqvmLZZMtfpfhF8xqAhQsX9rPY6qO9e/fy/ve/nwMHDlDs6+fC0d4oXwZOB54C3puZz0drpc8C7wR+BKzLzIfLKLsGq72HYkQc00MxM/fbQ1EaXZ4Hqt0gLhuphjA73UjHA28FvsvkXVs7OqmsTJdXf+HsSHEy2Tf+QlkP7T2Utm3bBjDXHkqyh6IkaYwJBUkdeemllwDOBD6Smf/Y/thUXVsn4wll9Z1yyins39/KGRX/jxQPPQ2c1rbqgmLZZMtVU/PmzePcc1vToJxwwgkAP+aVHkoTTb432XAYjRa7O0uSABMKkjpw+PBhrrjiCoBDmfnnxeLJJt/0pHJEXHbZZdx+e+ucsfj/QvHQZuD9xQScbwdeLHqrfBO4OCJOKiZjvLhYVi32SpmRp556CuDnGKUeSpoReyhKksaYUJA0pcxkw4YNnHXWWQAH2h6arGvrZCebqrA1a9awYsUKdu3axYIFC7j1Adi4cSNbt25l8eLF3HvvvQBjcbwHeALYA3wB+PcAxcnDHwEPFX//2ROK0fDSSy+NJRX32kNJkiSNmfGkjJKa4cEHH+SOO+7gzW9+M8CSiHgE+D9pdWWdaPLNe2hNyLaH1qRsvz30Qk/FX6cntGnTpmMX3Blw8sncd999RxdFxMtw9ATyQxO9TmbeBtw2sIJq6MZ6KF199dU8/PDDLxSLJ5t8zx5KkiQ1iD0UJE3pggsuIDPZsWMHwM7MXJaZ90zWtbUYO/2hzDwzM9+cmdvL/QSSZqq9h9J1113X/pA9lCRJkj0UJEnSxNp7KC1btgxavZTeSV17KGl69uKSJHXBhIKk+roz4H1dDd0efX4n6qOxHkpjImJnZt5T3H3VtcanGg4jSZJGj0MeJAkG/6tcP16/19fo5Pn+Oqm6qUKdrUIZJEkqgQkFSRqv25ODqp5MtJfLZIJGzajV11H7PJKkRjChIEndmqrhP9Fj45dN9vyZJDLGntPpcyd7TrfJB0mdq8s2VZdySpIqw4SCpNEwqIZwP4cZTJZYmGj5+L/xj3XznlMlOTp5PWlU9aPnzvhtaaLttpv369VE+wxJkgbEhIKk5prsRH6ydbvtmdBLWaZ7fKITBk8ipKnNdLvoZl8x0/UGwX2CJGnATChI0kQ6GQ4w3dCFmfxCaaNfqo7ptuGZbK/d7EN6MVnScbIeTO57JEkzYEJBUr3NtEdBp68x3XNmekJgI35qfi/q1vhtqpttc7JhC4Mo31Rl6+S5Ey136JIkqSQmFCTVX1VOzst+f0kt3QwJ6ma7HdYcCL3cn8n7uO+SJM3QrLILIEl9N1FD+X1ZTlkkDd5UE5xOtHyqdad7j37ptifF+P2a+zRJUgWYUJA0WgY5JlmSBmlYvSXcH0qS+sQhD5IkSZIkqWsmFCRVj7+eSVL33HdKkobMhIIkSZIkSeqaCQVJ0ujzl1tJkqS+M6EgSZIkSZK6NrCEQkRcEhG7ImJPRGwc1Puomox/sxn/ZjP+zWb8ZR1oNuPfbMa/eQaSUIiI44DPAZcCS4A1EbFkEO+l6jH+zWb8m22o8R/GMAaHSnTF7V/WgWYz/s1m/Jtp1oBe9zxgT2Y+ARARdwGXAzsH9H6qFuPfbMa/2UYn/mPJhJkmFd6Xrzz3ffnKa40tH1s2WkYn/pop60CzGf9mM/4NNKiEwnxgb9v9fcD57StExDXANcXdlyJi14DKMlOvB54tuxAdGFw5r56wEf2LHTzT+A9PLeMP09aB13N1VPX7r2rd6H+5Xl0HhhV/6LQOTFxPpzfT53XzPbe/x0S3Z16GiQyrXk5XB4Yb/+Fr8vYPtgGqZlTbAGUz/lB2G6BMxh+62gcMKqEwrcy8BbilrPefTkRsz8zlZZdjOnUp53jGvz/qUs6JTFUHqvy5qlq2qpZrMtPtA6r6eSxXfxj//qpquSZjG6A/6lLOiVS5DtTle61LOSdi/HtXpXIOalLGp4HT2u4vKJapGYx/sxn/ZjP+zWb8ZR1oNuPfbMa/gQaVUHgIWBwRZ0TE8cBVwOYBvZeqx/g3m/FvNuPfbMZf1oFmM/7NZvwbaCBDHjLzSERcC3wTOA64LTMfG8R7DVAlu+FMoHLlNP5DVbly9in+lftcbapatkqUq4/bfyU+zwQs1xSMf2kqUy7bAENVuXIa/6GqXDmN/1BVppyROZKzTEuSJEmSpAEa1JAHSZIkSZI0wkwoSJIkSZKkrjUyoRARp0XEtyJiZ0Q8FhEfLpbPjoitEbG7+H9SsTwi4uaI2BMROyLi3CGX97iI+OuI+Hpx/4yI+G5Rni8Xk54QEa8p7u8pHj99mOWsC+NfbRFxSUTsKj7HxiG/d6XrRhPqgvGfsnzGf7DvbfwboOpxnqC8xr2PjH+zGf/BaGRCATgC/IfMXAK8HfhQRCwBNgL3ZeZi4L7iPsClwOLi7xrg80Mu74eBx9vufxL4TGYuAp4HNhTLNwDPF8s/U6ynVzP+FRURxwGfo/WdLwHWFLEZlqrXjZGuC8Z/WsZ/sIx/M1Q9zuMZ9/4y/s1m/AchMxv/B9wNrAZ2AfOKZfOAXcXtPwXWtK1/dL0hlG0BrYr9K8DXgQCeBWYVj68Avlnc/iaworg9q1gvyv5+q/5n/Kvz1/55ivsfAz5m3WhGXTD+xt/4Nzf+xtm4G3/jb/zrGf+m9lA4qugS8lbgu8Apmbm/eOgfgFOK2/OBvW1P21csG4Y/Bv5/wL8U908GXsjMIxOU5Wg5i8dfLNbXJIx/5ZT5XR+jgnXjjxn9umD8J/fHGP+hMf7NUME4j/fHGPeBMf7NZvz7p9EJhYj4BeC/Ax/JzH9sfyxbKZ5Sr6kZEb8OPJOZ3y+zHKPK+GsyVasb1oXhMv7NZvyboWpxHs+4D5bxbzbj31+zyi5AWSLiZ2lVpC9l5p8Xiw9ExLzM3B8R84BniuVPA6e1PX1BsWzQfhm4LCLeCfwr4LXAZ4ETI2JWkYFqL8tYOfdFxCzgdcBzQyhn7Rj/yirruz6qonWjKXXB+E/M+A+J8W+GisZ5POM+IMa/2Yx//zWyh0JEBHAr8Hhmfrrtoc3A2uL2WlrjasaWv7+Y6fPtwItt3WIGJjM/lpkLMvN04Crg/sy8GvgWcOUk5Rwr/5XF+qVm2KrI+FfaQ8DiYhbb42l97s3DevOq1o0G1QXjPwHjPxzGvxmqGufxjPtgGP9mM/6DK3Dj/oALaHVl2QE8Uvy9k9ZYk/uA3cC9wOxi/aA18/TfAX8DLC+hzBcCXy9uvwH4HrAH+H+B1xTL/1Vxf0/x+BvK/q6r+Gf8q/1XxOJ/Fd/3H1g3mlUXjL/xN/7Njb9xNu7G3/gb//rFP4pCSJIkSZIkdayRQx4kSZIkSVJvTChIkiRJkqSumVCQJEmSJEldM6EgSZIkSZK6ZkJBkiRJkiR1zYSCJEmSJEnqmgkFSZIkSZLUtf8/vh/gm5R74ZEAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1296x432 with 16 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(2,8, figsize=(18,6))\n",
    "for i in range(8):\n",
Valerio Mascagna's avatar
Valerio Mascagna committed
    "    ax[0][i].hist(digiph[:,i], 100, log=True, color='red')\n",
    "    ax[1][i].hist(digitime[:,i], 100, color='orange')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
Valerio Mascagna's avatar
Valerio Mascagna committed
   "execution_count": 25,
   "metadata": {},
Valerio Mascagna's avatar
Valerio Mascagna committed
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA3cAAAE/CAYAAADlpzo+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABylUlEQVR4nO3de7xcVX3//9fnnNzIhSRAjJgQwyVeEG8QAbUPUSk3L6XtTyjaarR8S/v9YrVqRWj7LRW1X8RWSr9trVioUBUEqw/5KRRTFG1/CpKAlZtKCCBJCQRyISGQk3PO+v2x15q9ZmbPzJ6ZPff38/E4jzOzL2uvfTn7zJr12Z9lzjlERERERERksI31ugIiIiIiIiLSPjXuREREREREhoAadyIiIiIiIkNAjTsREREREZEhoMadiIiIiIjIEFDjTkREREREZAiocdenzOxWM/sfva5HYGb3mtkbe12PevrtmAVm5szsiF7XQ6Qo/fa3pvtT63R/kmHSb39nuje1Tvem1qlx10Nm9rCZPWtmu83scTP7opnN71Fddkc/01G9dpvZbzvnXuacu7UXdRtmZvZXZvaAme0ys5+Z2Xt6XScR0P1JwMwuMbNHzexpM3vEzP6k13US0b1JAjM7wMy2mtl/9rou/USNu957u3NuPnA0sBr4s15Uwjk3P/wAvwz18j9f7kWdRsQzwNuBhcAa4DIze11vqyRSovvTaLsCeIlzbn/gdcBvm9lv9rhOIqB7kyQ+Ddzf60r0GzXu+oRzbjNwE3BUNPmFZvb/+V6d75jZQVnrmtn9Zva26P0M/03G0WY2x8y+ZGZPmdkOM7vDzJY2Wz//Tdmv+td/YWbX+3J3mdndZvYiM7vAzJ7w3/SeHK270MyuMLPHzGyzmX3SzMabrUNU3ulm9hP/bfKDZnZqNLvmMfN13mJmO83sB2b2smjeF83s783s237d283s8Gi+M7M/8L1sO/yyFs3/XX8etpvZzWb2wjz74py70Dn3M+fctHPuduA/gNe2emxEOkH3p6bqMkz3p587556JJk0DCpOSvqF7U1N1GZp7k1/3dSTn/Z9bPSbDSo27PmFmhwBvAe6KJr8LeB/wPGAW8Mc1Vr8GeGf0/hTgSefcnSS9QQuBQ4ADgT8Ani2gym8H/gVY7Ot8M8n1tAy4CPh8tOwXgUmSDwWvBk4GWorvNrNjgauBjwKLgDcAD0eL1DtmNwGr/Lw7gcpv1c4CPu73aQPwqYr5bwNeA7wCOJPkOGNmpwN/AvwmsISkgXZNC/u2ny//3mbXFekk3Z/yGcb7k5mdb2a7gU3APOAredcV6TTdm/IZtnuTb+T+HfB+wOVZZ6Q45/TTox+SP6zdwA7gEeAfgP38vFuBP4uW/V/Av9Uo5whgFzDXv/8y8Of+9e8CPwRe0WS9frXWNOAvgLXRvLf7/Rj37xeQ/LEtApYCe8N++fnvBL7X4jH7PHBpjXnNHLNFvo4L/fsvAv8UzX8L8LPovQN+JXp/HXC+f30TcHY0bwzYA7wwWveIHPt2FfBvgPX62tSPfnR/aumYDeX9CTCSD5cfBxb0+trUz2j/6N7U0jEbqnsT8CHgc/71e4H/7PV12U8/6rnrvV93zi1yzr3QOfe/nHPxN0Nbotd7gMwHhp1zG0hijt9uZnOBXyP9dvVfSL4ZutbM/tuSB+RnFlDvx6PXz5J82zUVvcfX94XATOAx3yW/g+Qm87ysQq384eQVGYscAjxYp16Zx8zMxs3sYh+K8DTpN1YHNVo3x/wXkjwrF/ZvG8mHoWV16lnGzD5DEl5wpvN3K5E+oPtTZFTvTy5xF8mx+3je9UQ6SPemyCjdm8zsBcAHgD+tt9wom9HrCkhhQnjBGHCfv2nhnNtH8s/442a2ErgR+DnJg/Ld8CjJt08HOecmGy3skoeSG5V3eINlsrwLOB34VZKb00JgO8mNpF2PAp9yLT48bWYfB04DTnDOPV1AfUT6je5P9fXt/anCDFrbP5F+pXtTff14bzoWOBi4zz++tx+wn5ltAZZFjeWRpZ674XEtSTz2/yR6JsLM3mRmL/fxyU8D+0geiu8K59xjwHeAvzaz/c1szMwON7MTWizyCuB9ZnaiL2uZmb0kx3oLSG6UTwFzgb9scftZ/hG4IDxkbMlD0GfkWdHMLiC5ef6qc+6pAusk0k90f6qv7+5Pvv6/b2aLLXEscC5wS4F1E+k13Zvq67t7E0k450rgVf7nz0meX3yVGnYJNe6GhL8R/IgkXfVXo1nPB75GcnO6H/g+SbhBN72H5AHd+0i+8fkaybcuTXPO/Zjkod9LgZ0k+5Mnu9LVJLH5m309bmtl+zXq9A2SdLzX+rCFe0h64vL4S2AFsCEKqdBYUjJUdH9qqF/vT79BEsq1C/gS8H/9j8hQ0L2pob67Nznn9jrntoQfkv3Z518LPnGDiIiIiIiIDDb13ImIiIiIiAwBNe5ERERERESGgBp3IiIiIiIiQ0CNOxERERERkSGgxp2IiIiIiMgQ6OtBzA866CC3cuXKwst94P7/Lr12e54rvX7RMYcVvi0RqbZ+/fonnXNLel2PdnTq/iSD5RfrN5ZeV/4PiedlzZf+o3uTNPLAPZvK3q86annmcnn//ou4T+heM/yauTf1deNu5cqVrFu3rvByTznmwtLr6bvuK71eu+76wrclItXM7JFe16Fdnbo/yWA5aSwdc7fyf0g8L2u+9B/dm6SR01adV/b+pnWXZC6X9++/iPuE7jXDr5l7U1837kREREREuqnsS5vp7jeUer19GWwj37gbe/WR3Lz+472uhoiIDKB6H7z0oUxk+Nz0QHZPXaW8f/9F3Cd0r5GYEqqIiIiIiIgMgZHvuRMRERFphZldCbwNeMI5d5Sf9hng7cAE8CDwPufcDj/vAuBsYAr4gHPuZj/9VOAyYBz4J+fcxX76ocC1wIHAeuDdzrmJru3gCKh6Xq3LvWC93r4Mn5Fs3MVJVERERERa9EXg74Cro2lrgQucc5Nm9mngAuBjZnYkcBbwMuAFwL+b2Yv8On8PnARsAu4wsxucc/cBnwYudc5da2b/SNIw/FwX9muk9bqB1evty2BTWKaIiIhIC5xzPwC2VUz7jnNu0r+9DQi58k8HrnXO7XXOPQRsAI71Pxuccxt9r9y1wOlmZsCbga/59a8Cfr2T+yMig0+NOxEREZHO+F3gJv96GfBoNG+Tn1Zr+oHAjqihGKaLiNQ0kmGZj3/wdb2ugoiIiAwxM/tTYBL4che2dQ5wDsCKFSs6vbmh0usQyF5vX4aPeu5IHmatfKBVREREpBVm9l6SRCu/7ZxzfvJm4JBoseV+Wq3pTwGLzGxGxfQqzrnLnXOrnXOrlyxZUth+iMjgGcmeu8W/2Fd6PeumO3pYExERERkmPvPlecAJzrk90awbgK+Y2WdJEqqsAn4MGLDKZ8bcTJJ05V3OOWdm3wPeQfIc3hrgm93bExEZROq5ExEREWmBmV0D/Ah4sZltMrOzSbJnLgDWmtlPfJZLnHP3AtcB9wH/BpzrnJvyz9S9H7gZuB+4zi8L8DHgw2a2geQZvCu6uHsiMoBGsudOREREpF3OuXdmTK7ZAHPOfQr4VMb0G4EbM6ZvJMmmKSKSy0g27r7/7fNKr/WsnYiIiIiIDAOFZYqIiIiIiAyBkey5O+WYC3tdBRERERERkUKp505ERERERGQIqHEnIiIiIiIyBNS4ExERERERGQIj2bjb+dKFpZ9d73ota6evZ+309b2ulog0ycw+ZGb3mtk9ZnaNmc0xs0PN7HYz22BmXzWzWX7Z2f79Bj9/ZVTOBX76z83slJ7tkIiIiEgbRrJxJyKDz8yWAR8AVjvnjgLGgbOATwOXOueOALYDZ/tVzga2++mX+uUwsyP9ei8DTgX+wczGu7kvIiIiIkUYycbdxDwr+xGRgTUD2M/MZgBzgceANwNf8/OvAn7dvz7dv8fPP9HMzE+/1jm31zn3ELABDRosIiIiA2gkG3ciMvicc5uBvwJ+SdKo2wmsB3Y45yb9YpuAZf71MuBRv+6kX/7AeHrGOiIiIiIDYyTHuat00tgZAHruTmSAmNlikl63Q4EdwPUkYZWd3OY5wDkAK1as6OSmpMfC/wXQ/wYRERkcI9lz9/Rh6c+BX/hhr6sjIq35VeAh59xW59w+4OvA64FFPkwTYDmw2b/eDBwC4OcvBJ6Kp2esU8Y5d7lzbrVzbvWSJUuK3h8RERGRtjRs3JnZIWb2PTO7z2el+6Cf/hdmttnMfuJ/3hKtk5l5zsxO9dM2mNn5ndklERkRvwSON7O5/tm5E4H7gO8B7/DLrAG+6V/f4N/j53/XOef89LN8Ns1DgVXAj7u0DyIiIiKFyROWOQl8xDl3p5ktANab2Vo/71Ln3F/FC1dknnsB8O9m9iI/+++Bk0ieabnDzG5wzt1XxI6IyGhxzt1uZl8D7iS5T90FXA58G7jWzD7pp13hV7kC+Bcz2wBsI7lP4Zy718yuI2kYTgLnOuemuroz0jficEwREZFB07Bx55x7jCRZAc65XWZ2P/WTDZQyzwEP+Q9SIfPcBufcRgAzu9Yv2/XG3eHX7UjfvPpIbl7/8W5XQUQK4Jy7ELiwYvJGMrJdOueeAzI/uTvnPgV8qvAKykDTs3YiIjJomkqo4gf9fTVwO8mzLe83s/cA60h697aTNPxui1aLM89VZqQ7rrVqi4iISD+r7AVVY1lEpPNyN+7MbD7wr8AfOeeeNrPPAZ8AnP/918DvtluhbmSj2/nShR0pV0REBpsaICIiMshyZcs0s5kkDbsvO+e+DuCce9w5N+Wcmwa+QBoGVSvzXK6MdMpGJyIiIiIi0rw82TKNJBHB/c65z0bTD44W+w3gHv+6Vua5O4BVZnaomc0iSWZwQzG7ISIiIiIiMtryhGW+Hng3cLeZ/cRP+xPgnWb2KpKwzIeB34f6mefM7P3AzcA4cKVz7t7C9qQJOw63XmxWRERkZCjEVUSk+/Jky/xPIKs1dGOddTIzzznnbqy3noiIiIiIiLSmqWyZIiIiMjqU8VJEZLCMZONu6bp9va6CiIiIiIhIoXJlyxQREREREZH+psadiIiIiIjIEBjJsMw5W57pdRVERET6np6xExEZLOq5ExERERERGQIj2XO36aRFva6CiIiIiIhIodRzJyIiIiIiMgRGsudORESGXz+O0daPdRIRkeExko27qdm9roGIiIiIiEixFJYpIiIiIiIyBEay505ERKSeOHyy1dDJVkIwi9iuiIiMrpFs3C160PW6CiIi0mH92DjqxzqJiMjwGMnGnYiIyKBSUhYRaUT3idE1ko27iXnW6yqIiEgfK+KDUCtl6APYYDGzK4G3AU84547y0w4AvgqsBB4GznTObTczAy4D3gLsAd7rnLvTr7MG+DNf7Cedc1f56ccAXwT2A24EPuicU/iRiNSkhCoiIiIirfkicGrFtPOBW5xzq4Bb/HuA04BV/ucc4HNQagxeCBwHHAtcaGaL/TqfA34vWq9yWyIiZUay505EhoOZvZjkG/LgMODPgasp6Jtz6R/1ko20E4Kk8CVplXPuB2a2smLy6cAb/eurgFuBj/npV/uet9vMbJGZHeyXXeuc2wZgZmuBU83sVmB/59xtfvrVwK8DN3Vuj0Rk0KlxJyIDyzn3c+BVAGY2DmwGvkH6zfnFZna+f/8xyr85P47kW/Hjom/OVwMOWG9mNzjntnd3j0RkCCx1zj3mX28BlvrXy4BHo+U2+Wn1pm/KmC7SkL6kGl0KyxSRYXEi8KBz7hGSb8hDz9tVJN92Q/TNuf82PHxzfgr+m3PfoFuLwp9EpE2+l67jz8iZ2Tlmts7M1m3durXTmxORPqaeOxEZFmcB1/jXRX1zLkOiMvQyaPTtdlEhmxq/bqQ8bmYHO+ce818ePeGnbwYOiZZb7qdtJg3jDNNv9dOXZyxfxTl3OXA5wOrVq5VwRWSEjWTjbtYzuu+JDBMzmwX8GnBB5TznnDOzQv7ozewckkQIrFixoogipQn1GkXtNJhmHLay5XV7QY3DvncDsAa42P/+ZjT9/WZ2LUlY+E7fALwZ+MsoicrJwAXOuW1m9rSZHQ/cDrwH+L/d3BERGTwKyxSRYXAacKdz7nH//nH/jTlNfHOeNb2Mc+5y59xq59zqJUuWFLwLIjJozOwa4EfAi81sk5mdTdKoO8nMHgB+1b+HZCiDjcAG4AvA/wLwiVQ+Adzhfy4KyVX8Mv/k13kQJVMRkQZGsudORIbOO0lDMqGgb867UnPJrdXQShnNsNBuZEF1zr2zxqwTM5Z1wLk1yrkSuDJj+jrgqHbqKCKjRY07ERloZjYPOAn4/WjyxcB1/lv0R4Az/fQbSYZB2EAyFML7IPnm3MzCN+dQ/s25DLmbNnym5ryiGgSj0qASEZHeUuNORAaac+4Z4MCKaU9R0DfnMnhq9fB1elt5GnAaU09ERDpJjTsRERkI7TaERrkhNYr7Por7LCKihCoiIiIiIiJDQD13IiIiFdpJQNJs4heFahZrFJPHiIgEI9m427PEel0FERHpgk5/uG+2fDU2RESkkxSWKSIiIiIiMgRGsudORET602lHfLT0+qYNnyl7H6a1U14j3cy0KSIiUrSRbNzteX6vayAiIp0SQh8rG4atltNpCtUslo6niIyykWzciYhI/4h7y2YctrLt8tpt1MX6MdlJP9ZpkKm3VkSGyUg27qbmTve6CiIikmFy48NljZVmwzD7gRpbIiLSK0qoIiIiIiIiMgRGsudORESGQxFjmnUjLE9jr4mISDeMZONu8T3qsBQR6RdFN3baDeWM69OPDbF+rNMgi4+nmcbBFZHBplaOiIiIiIjIEBjJnjsRERleeUMgOx2OqSyM+dTL/qnMoCIizWnYc2dmh5jZ98zsPjO718w+6KcfYGZrzewB/3uxn25m9rdmtsHMfmpmR0dlrfHLP2Bmazq3W/XtPqT8R0REpNPWTl+vxomIiHRUnrDMSeAjzrkjgeOBc83sSOB84Bbn3CrgFv8e4DRglf85B/gcJI1B4ELgOOBY4MLQIBQREREREZH2NAzLdM49BjzmX+8ys/uBZcDpwBv9YlcBtwIf89Ovds454DYzW2RmB/tl1zrntgGY2VrgVOCaAvenJSHsQ9+oioj0VqeySnYqRLKd+iqDZm/p+IvIMGrqmTszWwm8GrgdWOobfgBbgKX+9TLg0Wi1TX5areldN+ep9PXSy37YiyqIiEiHtPJBvRMf7tVgyKfecdIxFBFpTu7GnZnNB/4V+CPn3NNxumDnnDMzV0SFzOwcknBOVqxYUUSRIiIyQk474qMATG58uLcVEREZYur97k+5GndmNpOkYfdl59zX/eTHzexg59xjPuzyCT99MxCnKVnup20mDeMM02+t3JZz7nLgcoDVq1cX0mCs9PRLJ9PX/3gsD5/z0U5sRkREmjTjsJVNLd/MB4q109cXnn2xnfX1Yai3dPxFZBjlyZZpwBXA/c65z0azbgBCxss1wDej6e/xWTOPB3b68M2bgZPNbLFPpHKynyYiIiIiIiJtytNz93rg3cDdZvYTP+1PgIuB68zsbOAR4Ew/70bgLcAGYA/wPgDn3DYz+wRwh1/uopBcRURERtcpr/l46XWeFM5ZQihmK2olW8nTs6OwpM7q5Dh3GodQRIZRnmyZ/wlYjdknZizvgHNrlHUlcGUzFeyE/e/X2O0iw8DMFgH/BBwFOOB3gZ8DXwVWAg8DZzrntvsohMtIvnzaA7zXOXenL2cN8Ge+2E86567q3l5IbHrxfG6+48K2yoiftatsDKgBJiJSDN1P+1OrX5KKiPSDy4B/c869BHglcD8ag1NERERGlLqwRGQgmdlC4A3AewGccxPAhJkNzRicki0OwQy9dHkTsRQR5tfrcL5Ohiq2qld1iq+FmzZ8pivbFBHpZyPZuHvuwF7XQEQKcCiwFfhnM3slsB74IAM8BueoajcME5JGXlENimbL6YfG1bDq5LHNKjse5klEZBCNZONORIbCDOBo4A+dc7eb2WWkIZhAsWNwgsbh7JSsXp9WkpzU68WpVV6z07O2H5ZttiGiZCy9peMvIsNoJBt3M3f3ugYiUoBNwCbn3O3+/ddIGncdGYMTujMOpzSWFX7XTLbM8EG+1Q/38bK9CNHsx4ZIr+qkgepFRMopoYqIDCTn3BbgUTN7sZ90InAfGoNTRERERtRI9tw1qx8fXhcRAP4Q+LKZzQI2koyrOYbG4Bxqee7J9UIl2xkTrxmV22kl4UeekNVa0/pJv4RA9joZjohIp41k425qdq9rICJFcM79BFidMWsgx+AcVY0+7OdpDMQNp6wP8KGMvA27vA2Qdhoq/dYAG0TtHkOdAxEZNgrLFBERERERGQIj2XOXl8I3RES665TXfLzmvFq9bmF6vbDHRj12tXr78oY7dissst3/S/0evikiIu0Zycbd3C3Nr6N/gCIi/WHGYStLDbmsRluj+3VYt9Nf4BUxqHatzJyV+9jv/6P6vX6dYGYfAv4H4IC7SZ7zPRi4FjiQZGzOdzvnJsxsNnA1cAzwFPBbzrmHfTkXAGcDU8AHnHNK+CQiNSksU0RERKRAZrYM+ACw2jl3FDAOnAV8GrjUOXcEsJ2k0Yb/vd1Pv9Qvh5kd6dd7GXAq8A9mNt7NfRGRwTKSPXciItI/4h6psWOOqjl/xmErc02L9arHqNWBzSvXb6eMTunH0M48GUR7YAawn5ntA+YCjwFvBt7l518F/AXwOeB0/xqSMTv/zszMT7/WObcXeMjMNgDHAj/q0j6IyIAZycbdzhfnG3u4H/5hiYiMkun195ReV35An9z4cK5pWbIGu84azDyeXvm6nm4MbF7E/yT9X+vOMXDObTazvwJ+CTwLfIckDHOHc27SL7YJWOZfLwMe9etOmtlOktDNZcBtUdHxOiIiVUaycSciIpJnWIR2x5Sr7IFrp9er2d68RtvqZu9gq9tq9Xj1utfOzBaT9LodCuwAricJq+zU9s4BzgFYsWJFpzYjIgNAjTsREempOLQyj6wP+O2MlRcnaOkX/dzD1o91y5t4pot+FXjIObcVwMy+DrweWGRmM3zv3XJgs19+M3AIsMnMZgALSRKrhOlBvE6Jc+5y4HKA1atX5wtPEpGhpIQqIiIiIsX6JXC8mc31z86dCNwHfA94h19mDfBN//oG/x4//7vOOeenn2Vms83sUGAV8OMu7YOIDCD13ImISNcUkYwjb6hkHHYZnrlrtpdwmLSS5KWV89Wo56xPetY6yjl3u5l9DbgTmATuIulZ+zZwrZl90k+7wq9yBfAvPmHKNpIMmTjn7jWz60gahpPAuc65qa7ujIgMlJFs3I3vsV5XQUREvKxkJ0FW9sxubDdoJQQ0yGpIdrMx0yjJSy/q0uyzcK3UsV8ajM65C4ELKyZvJMl2Wbnsc0DmwXHOfQr4VOEVFJGhpLBMERERERGRITCSPXd5jULoiIhIvzpp7Iyqnrusnp9a01oJwSwyu2Un9er/U7vj9+Utv+ht6P+5iIyKkWzc7Vs43esqiIiMpLxhjuHD+PT6ezJD+vJMC2PgVU4Lpk84uuX9qKfVEMROlttqg6aV8fvqbUsNKxGRzlJYpoiIiIiIyBAYyZ47ERHpH1lZLWN5QzGzBiXP29uUd912whI7HdLYrHqhikWEMfY6pLXWuVeIpogMs5Fs3C2+J1+HpW76IiK9U+uZt0bDHoRByePl4ufvKgcsz2rYVZYP7WV6LDJEs1f/m7qx3U5uI+vaEREZNiPZuBOR7uj1N/cyGrKSp9RrsMXzKht6tcrPM61e71cz5bc6Dl2//H0NWgNq0OorIlLPSDbunj6s1zUQEZGgiHHu4jJCQy/04MUa9dJVllFZt3Z6f/ql8RV0I/FJL/e5nTEKRUQGlRKqiIiIiIiIDIGR7LkTEZH+l3ecu3oahWDm7cnLo1Hd8oRy5i1z2HugujHeXdFli4j0g5Fs3M3d0usaiIwGfXCSPLIGGw/hkGPbd1dNq6Uy8UlWudOL55dex2XnefaucjuxRo2weg2/dv5O+vFvrB/rVE9cXzPrYU1ERNqnsEwRGWhm9rCZ3W1mPzGzdX7aAWa21swe8L8X++lmZn9rZhvM7KdmdnRUzhq//ANmtqZX+yMiIiLSqpHsuRORofMm59yT0fvzgVuccxeb2fn+/ceA04BV/uc44HPAcWZ2AHAhsBpwwHozu8E5t72bOzEKOhkWd8prPt7yunl61vopVLAXWTL7JTOnQitFRGobycbdnuf3ugYi0mGnA2/0r68CbiVp3J0OXO2cc8BtZrbIzA72y651zm0DMLO1wKnANd2t9miKQyIrn4GLQzHD83dxOGVYP27YjR1zFDffcWHD7U4vns/Nd1xYSCr8Rg2Mfhsrr991qsGmhqCIDLuRbNxJsfrl21wZWQ74jpk54PPOucuBpc65x/z8LcBS/3oZ8Gi07iY/rdZ0GVChkdjMs3T9ItxTe9EjV7ndTjYo9b9DRHolz/1nUKME1LgTkUH3K865zWb2PGCtmf0snumcc77h1zYzOwc4B2DFihVFFDly2vnnOL3+HtZOX5+Z4XJ6/T2l11mJVII8PXpBZV2L+Mder27N6NWHjEH5cCMiMqqUUEVEBppzbrP//QTwDeBY4HEfbon//YRffDNwSLT6cj+t1vTKbV3unFvtnFu9ZMmSondFREREpC3quRtSCneRUWBm84Ax59wu//pk4CLgBmANcLH//U2/yg3A+83sWpKEKjudc4+Z2c3AX4asmr6cC7q4KyMjK8yl3bHm6oUOVg6LEA+FAM315HVbq+GZw3T/H5RkMyIi/WIkG3eTC6Z6XYWhon+e0kNLgW/4salmAF9xzv2bmd0BXGdmZwOPAGf65W8E3gJsAPYA7wNwzm0zs08Ad/jlLgrJVaS36o2BF5s+4ej0zaPbai7Xa1nPAPY6aUq9e3gn7++tlq3/OSLSrjz3kUG914xk405EhoNzbiPwyozpTwEnZkx3wLk1yroSuLLoOoqIiIh0ixp3dD8zmfTWoGY/EpH2xrJrVuX/hsoQz6CyN25QwgL7uZ61ejAVqioivTBI95CGjTszuxJ4G/CEc+4oP+0vgN8DtvrF/sQ5d6OfdwFwNjAFfMA5d7OffipwGTAO/JNz7uJidyW/GbvGS68PveCHvapGR/X7hScio6nRvSk8ExfGqoufx6ucB+k/3LHv35kuF8bDyxjvLqvRkFWnXo0pV288vLz39WG7/zc7RuCw7b+ISDPyZMv8IslgvpUudc69yv+Eht2RwFnAy/w6/2Bm42Y2Dvw9cBpwJPBOv6yIiIiIiIgUoGHPnXPuB2a2Mmd5pwPXOuf2Ag+Z2QaStOQAG/zzMfhMdacD9zVf5e5TGF//6adzMkhd9SLd1O7fab0smnkzbDbq7Tlp7AzGfE8fwNj23fkqJ2WyznWv7tO6J4s0pr+TRKN7V7MZevOu00ntPHP3fjN7D7AO+IhzbjuwDLgtWmaTnwbwaMX049rYdlsmnz9Rev3AP6/m4TUf61VVpAd6/UcnIuUmDjmg9DqEV2ZlyRzbvpubNnymrGE347CV3LThM2X/XOMBzbPkCfML5WZplIWziHtMN+5T/XwvzKpbO/Xt530Vkf43SPeQVht3nwM+ATj/+6+B3y2iQmZ2DnAOwIoVK4oosiElVBERGQwnjZ1R1fDr5PNxecuutVy9/yul5wWjnkNobey9Rt8c9+O39P1YJxEZbb163rpILTXunHOPh9dm9gXgW/7tZuCQaNHlfhp1pleWfTlwOcDq1atdK/VraG/6qOGqP/hxw8X1T6f/9NM56ae6iPSTvH8bcTKUyvWy/tFObny4bu9b3Pir1fsWm15/T+6kHVnhhlJ8T1s7dE8WaUx/J4la/2+aOT79dizzJFSpYmYHR29/AwgxMDcAZ5nZbDM7FFgF/JhkYOBVZnaomc0iSbpyQ+vVFhERERERkVieoRCuAd4IHGRmm4ALgTea2atIwjIfBn4fwDl3r5ldR5IoZRI41zk35ct5P3AzyVAIVzrn7i16ZySlcJfa+u3BV5F+1Mt7SHimrl6vW7O9Z7VCIItQry6VxzFvvePnCuuNoxfKHRSNxq+D2vvTbOhrZdKdPL24IiKDLk+2zHdmTL6izvKfAj6VMf1G4Mamatch47vTce7GXq0RGURE+kXl83TxB/cwgHllwpTKRlOjhCf1VDYUshK75F03T0OuqIZZo3L6tQGoMFcR6Rf9ep9sVkthmSIiIiIiItJf2hkKQUREpGMqx6BrtYx21mum564I9babN/wTmh+bqTKhTLNJaYqqU6MyRESkvpFs3E3PTJNw/uK9+/PQH36kh7XpjGHpWu4EHRuR/tUopDLvMAGV4X6lxoofKw/yNRzizJz1tpVX3sZio+0O6n2s3ayazSyrZ+xEZBSNZONu0JxyTPmHmZvXf7xHNRER6a3TjvhoVQOw2YZOMz1gldPaaVRljdE3LIpK8tKJxCp5yxWR0dPOvatf7ysj2bibnjfV6yqIiIhXb1y7OHlKvYZR1j/WePl2kqy0Kmu/GtWj1d66Vj5YdHpsuiLK6qcPTCIig0AJVUREREQKZmaLzOxrZvYzM7vfzF5rZgeY2Voze8D/XuyXNTP7WzPbYGY/NbOjo3LW+OUfMLM1vdsjERkEI9lzVyl8o6pvCIvTStiMiAy/To3RFo9pFnrHBiEEMs+YfnnkPa6DPEbeALoM+Dfn3DvMbBYwF/gT4Bbn3MVmdj5wPvAx4DRglf85DvgccJyZHUAyvvBqkrGF15vZDc657d3fHZHRMqj3y5Fs3NmklV6v+oMf97Am+egZOxHptukTjm68UIfVG78uzz/ZOClJs184xVk6K5O4tPIPPiv7ZOUg25JPUSGonWRmC4E3AO8FcM5NABNmdjrwRr/YVcCtJI2704GrnXMOuM33+h3sl13rnNvmy10LnApc0619ERlm7dwb+rWxp7BMERERkWIdCmwF/tnM7jKzfzKzecBS59xjfpktwFL/ehnwaLT+Jj+t1vQyZnaOma0zs3Vbt24teFdEZJCMZM+dFKtfswUNi0EMC9A1MVziHqZuppc/5TXtRy3kGe6g0TJxPSp78WqNSzfoafh7dc6LkPf+0+Ex9GYARwN/6Jy73cwuIwnBLHHOOTNzmWs3yTl3OXA5wOrVqwspU2SYDfPnlJFs3LkZ0Th3//QaHvnd83pYm+E0bH8oIqNmxu4JAMa27y603Eb3htBIml48v+a8LFnDI2R9MZInFDKuYxENTMhuIBXZaMp7z+33e3O/168Jm4BNzrnb/fuvkTTuHjezg51zj/mwyyf8/M3AIdH6y/20zaRhnGH6rR2st4h4g3o/UlimiIiISIGcc1uAR83sxX7SicB9wA1AyHi5Bvimf30D8B6fNfN4YKcP37wZONnMFvvMmif7aSIimUay5y7LMHfPVhrEML9h0M1rrJWwpEG9DsxsHFgHbHbOvc3MDgWuBQ4E1gPvds5NmNls4GrgGOAp4Leccw/7Mi4AzgamgA845/ThaUB0OLSuZZX1qtXr2Kj3Lm8iGN3X+9IfAl/2mTI3Au8j+VL9OjM7G3gEONMveyPwFmADsMcvi3Num5l9ArjDL3dRSK4iIsWo/CyU9dlokD4vjWTjbr9HZva6CkOl3y/yQTeIx7fLdf4gcD+wv3//aeBS59y1ZvaPJI22z/nf251zR5jZWX653zKzI4GzgJcBLwD+3cxe5Jyb6uZO9JvwbFmvGk9j23dz04bPVIVR1msM5bnumg2FjAdRrydPZs5B0uxx6qf7VDshqmaWsWRrnHM/IRnCoNKJGcs64Nwa5VwJXFlYxURkqO7XlUaycTcMivyWVuP89c4gfRPUilOOKU8+UfSwHma2HHgr8Cngw5Z8Mnsz8C6/yFXAX5A07k73ryF5/uXv/PKnA9c65/YCD5nZBuBY4EeFVnZA9OIfXbzNuIfrpLEzmhqrrlZyk8r57f6tFX2Msu7njbYxCPeOIv5PqUdSpHsG4b4ijY1k427ixc9WTRulizje12H8xqJfdfMaG5bkCjn8DXAesMC/PxDY4Zyb9O/jtOGllOLOuUkz2+mXXwbcFpWZmWpcuiP0gOW9N4WGXGVClTAv9EDVK6/e30E/jKnWKASzE9sUERkFWffOvNP6lRKqiMhAMrO3AU8459Z3cZsaS0pERET61kj23FVSWGJ/GqWkIEXoZvhSn4RKvR74NTN7CzCH5Jm7y4BFZjbD996FdOKQphrfZGYzgIUkiVVqpSCvorGkuqPTEQWt3jO6FenQ6+00Oibt/P0Xfe/Q/V9EWjHM946RbNxN7xsvvV71vnU9rEnriroQh+2CHjTDfvyLfsYu5py7ALgAwMzeCPyxc+63zex64B0kGTMrU42vIXmW7h3Ad/0gwjcAXzGzz5IkVFkF/LhjFR9Anb5Oa4WKr52+vuE4cyHsMk68UitEs9a8VoQ65xk3L29ZkO8Dx6DcN4qo56Dsq8gw0N/bcBjJxp2IDLWPAdea2SeBu4Ar/PQrgH/xCVO2kWTIxDl3r5ldRzIG1SRw7qhnyhQREZHBpMadDJVh7mYfdJWZM4vknLsVuNW/3kiS7bJymeeAzDg059ynSDJuSoa4d6rZFPmVmskMedoRH4XF8+uuHzTKkNlM1s08ssodpvtPr/el1ZBREZFW9fq+V5SRbNzNmjfR6ypIDs1mfBz1zJ+DfCOS/rHvlNcAMPPmZMzkohtFeYXtNjPOXOW68bR4HLpWZd1rms3umXcb3dbMdtupY6cziYqI5DHM946RbNzVk+db5W5fEJ18AL3V8np9TCRVxDdNnR6PTqQfjPIQAu3cJ+qt26lyRUS6YRg7BkaycTf1y3ml1xv/6rU8+OEP97A2IoOnlQ9icYPR7KIiqyMF2u8BP8RDh3rsGo0fFBKoTAM333FhS18kZTUaivwHPuhjIAV5xvfr1QefQTyeIjL4huHeo3HuREREREREhsBI9twNm06FtnQjjKeVsosqoxWNegSKOE4hRHL6rvsKLbdIjY6/wq2klnb/divHJe3Ftd7N+0/lUAvtJrSJ9eLvdBhDoERkOA3qI0gj2bibmt9clvNBOZl5DevYQ/1Yp24oYr/1jJ0EE4ccAMCM3UniqbwJTYrSqe3FiWGKbCANilbvE43WKypJjYhILwzjPUhhmSIiIiIiIkNgJHvumjWo3bIiIqMkK+SvU+PcFUn/YzpPoeIi0sw4q4N8n1DjrkI/nsxOh8V0otw8Wdja0U9jQRV9nIoIkez08enVNSmdN/b9O5PfvjE0dtjKQsMYG10blY2wePl6Qxk0+qfd6jh3vc6MGY/VV7Re/J3q3iAig2JQ71cKyxQRERERERkC6rmjOvta5fR+lKdunc4cV0T5ox4qc8JbLwFg1k13lE2vdSxG/XhJb8QZGzuRjKQyI2Teef1EoZWDR/dTkdHW6LN0KyGb/XBfGfnG3ao/+HGu5XTjF5FRMHbMUQDsXj4PgPl3P97L6pSZ0WSI6KDctzsdxi46jiJSO8x/2O4PI9+467R+/za3U2PZtVuXeiq3FcaFg+6l9O/GOG/9cuNp9M1Vv9RThs9JY2c0TITSToRF0X+nRdSjnbr027ZERLqtG59r+/1zz2g27iat1zVoW3wx9aJxU5RayRBGxfe/fV5Ty/fjTUSG2/Ti+R3fRtwbd+Kb/g+QJnZpNRFKLwxKPSWh8yUi9WTdIxp9Xu2H+4oSqoiIiIiIiAyB0ey5q6MbPUjNjrMR98xB//bO5X3wtN9DkPq9u71SK2EEsUHYRxl89a7TsmvzhKMbljUoSVYGhUI1RUTK1Uq2OAhGsnE3tmgi97LtntRuXBTtNPaarV8/jqnXi8buIP6xt0pj2o2WXYcmYZjj+xwAe5bP4z++8cdd2/7Y9++s+4VXpWaTrMSKuHaLGBy9m39D+nsVkWGW9x5XbwzTev9/BuEeqrBMERERERGRIdCw587MrgTeBjzhnDvKTzsA+CqwEngYONM5t93MDLgMeAuwB3ivc+5Ov84a4M98sZ90zl1V7K5IlnbCbbLCqFoJWWx2nX5LEJMnjLbTGTGbPSbdzNbZblio9L9Wr6dWz3m710orIeL99m1ss3XrVWhlp45hP58bkWEw6n9jw/yZJE9Y5heBvwOujqadD9zinLvYzM737z8GnAas8j/HAZ8DjvONwQuB1YAD1pvZDc657UXtSDNmbNyvF5stqdcVnKUfGjitqBVa1e8hSEXUr8h9HHv1kYWUM+qZSSWfyblJNuHxnUlY5pxv3l5Y2bX+Lpode6gTA6m3qp/q0qpR/GAnIjKsGjbunHM/MLOVFZNPB97oX18F3ErSuDsduNo554DbzGyRmR3sl13rnNsGYGZrgVOBa9rfhc5q5oHKTn1zWmRClVr7Myjf4IR6VjZ42u3NKuLc9WtSgkbnNs81Pixjv8hw6Ne/NRER6Q+j/H+i1YQqS51zj/nXW4Cl/vUy4NFouU1+Wq3pPTF52LO92nTX1fpGvNl1601rdfla+q2nsqj9KmK7eZer10vXzjVRb/vNjP2SRHBLP1p49Y8A2PWu1wIwp42y+r3nvJ//2fdLsqtebbefz43IMNDfWGrYjkXbCVV8L50roC4AmNk5ZrbOzNZt3bq1qGJFZAiZ2Rwz+7GZ/ZeZ3WtmH/fTDzWz281sg5l91cxm+emz/fsNfv7KqKwL/PSfm9kpPdolERERkZa12nP3uJkd7Jx7zIddPuGnbwYOiZZb7qdtJg3jDNNvzSrYOXc5cDnA6tWrC2s0FuWEt14CwKyb7mi4bCvhbnm7kdvtbi4yfK5TyTu6EeJXr7ep6DoVmRSl1SQvRT1vl7ecLoRF7AXe7JzbbWYzgf80s5uADwOXOueuNbN/BM4meQb4bGC7c+4IMzsL+DTwW2Z2JHAW8DLgBcC/m9mLnHNTRVdY6mvmuq83r5nQ80Z/T8Pwre6w7U+zRn3/RaSxYblPtNq4uwFYA1zsf38zmv5+M7uWJKHKTt8AvBn4SzNb7Jc7Gbig9Wq3x8bTNuPGa17Fg7/1p6X3RY1rV2TSiqLCFIcpocb0XfcVHrY0CGMadkKeJBftltUpPnJgt3870/844M3Au/z0q4C/IGncne5fA3wN+Duf5fd04Frn3F7gITPbABwL/Kjze9FfQjjm2ERynxw75qheVkdERKRpozxGb56hEK4h6XU7yMw2kWS9vBi4zszOBh4BzvSL30gyDMIGkqEQ3gfgnNtmZp8AQnfXRSG5iohIO8xsHFgPHAH8PfAgsMM5N+kXiZ/xLT3/65ybNLOdwIF++m1RsT19LlhERESkFXmyZb6zxqwTM5Z1wLk1yrkSuLKp2nVJrXCezB6u017TjSqVaSbcLszLCgFstseuUVhinrq2ktWySPFxmL7rvrrLttqjGbZRWf7a6eurMp22KoQDA3z/2+c1Va9mheMwEV3rebfZCz508lVmtgj4BvCSTm3LzM4BzgFYsWJFpzbTdWXXvu+56+Q2mgm37EejmoVtWEKWRGS05B37tNJpq8o/+9z0wCU1luwvrYZlDrSpifHS61Vr1je17qyb7qg5flus09na2s14WLndIkIS2/mA1svsdrXqXtT4crHKxm2e/ZnI+YVCXFZRjco826onb92L4JzbYWbfA14LLDKzGb73Ljz7C+lzwZvMbAawEHiK2s8LV26jr58JLsLCn+8C4KlXLABgcv6srtehlVDhdkJwepURt5uGbX+a1av995EF64DNzrm3mdmhwLUkEQPrgXc75ybMbDbJeMLHkNyTfss597Av4wKS54WngA84527u/p6IyKBoO1umiEivmNkS32OHme0HnATcD3wPeIdfrPK54DX+9TuA7/qIgxuAs3w2zUOBVcCPu7ITIjLMPkhyTwo+TZLs6QhgO0mjDaJkT8Clfjkqkj2dCvyDbzCKiGQayZ67TupmhsdOlN9uaGI728gT+pm1nVrhlkUen1KW1GharQHVB02711O9sNRayxboYOAq/2FnDLjOOfctM7sPuNbMPgncBVzhl78C+BefMGUbyYcmnHP3mtl1wH3AJHCuMmX2zmlHfBSAmzZ8JnN61rwsoxo+Kf3BzJYDbwU+BXzYJ29SsieRnFq5h8ehlI3CKDv5ebqXRrJxN2u/faXX8QfzZkLmOpVmvlF4UKNBqYOsBk8R2Q/zfmAPx/W5589repuNtteM+Jw22xgLGTnj593qLVdrXqvaCUEL+115biqv8XrXUxHXc/ysXidCRZ1zPwVenTF9I8kHoMrpzwGZO+ac+xTJh7CRk3WPmTz+dQDMerSY3Fd5r+e4Adevhu2DQF6jut9t+BvgPGCBf38gSvYk0nV5P7NXLjcoz9hVGsnGXbPyjCnWSjlQv3FZRPlQ3bMSb7PZxsfxv/PZ0uvbvvThZqvYkni/mukly/uNTMOewhwN1HCMi0heE8ubeKKZMRLjxmrW01RZvZRZvXN5Pui1eu6ku3qd0CTe/ozDVvauItI0JVnJZmZvA55wzq03szd2YXtDmexJJI+y/yGHH9rW+pUafeHdyfteq/+bR7JxNxklVJm+66eFlt3pf27NJKhoJVtjKx/AQ+9c3BgIjYBZNJ+8pd78rF6yvBd/s72t4VjM2fJM3YZ3p5OXtKNR433itNfUzIQZ5jXav2YbfNLfnnnH8QCMTfS4IhnyhGLGdD1KD70e+DUzewswB9gfuAwlexLJrZV7+OSDD+Vab5jGfq6khCoiIiIiBXLOXeCcW+6cW0nybO93nXO/jZI9iUiHjWTPXbvytPLbCdts9puKWj0rJ42dUUgoXLtjYOUJ26kVvteJcfF6/S1NM9tv52HfvNupWq4HYznK8FP4nggAH0PJnkS6qpVxm7PWHRQj2bibMSu9L+YJc2x3TLlGZTZT7qyb7qg5b+zVR1Y985UVllcvtLFRGN+Cr0QJur6U/AphffE+7PKNwLLlO6SorvXKY5I33DIraUw7yWvyhIrWW79RGXEYZtZyYSzHoFEDO+++xted2UW51pHuKBsj8TXhfM8HYHrx/K5uXwaLzl1jzrlbgVv9ayV7EilYEWM1Q+PPkd2+38XbS5Ln5qOwTBERERERkSEwkj139dTqrckbHlhUlsR2ymx2u/W+iai37VbrVcTYU50Yv6reUAmV26vV41t5XFsJwWy2bs1q5ZrL+larXk9nrR7gRsNKiIiIiDSjmc+EWZ+BTj52uKKJRrJxt++xuaXXs266rc6S1Zq5YJoNvczTQMn6cJ/V8Gz2wm4mk2S9RktVOV+qWY2a67XacGtnuVr7c/P6j+dqDH3/2+c11ZDLk1o3a17WcBZ5t9VundrVTKZX6Z29B+0HwH5bJ8ret0vheyIi0gtZn29s9curlrPVL+c7P/7zrtWrU0aycVdLnmesmk0KUC/ZSda0dsagi4UekjlbnimbntUQbDZpxwlvvSRzfLSsMuNyazUaGo4zR+Oe02bHmQt1anTO2xnmIM9+tbJuq7LGr2tHGPNw9vZkPN45BZUrIiIi0qp2ktGdfOxFjG/fXXo/iAOZj2Tjbmxf9UOJoVG186ULS9MW3r+zqXI7+c10kb0qeRNvZG0z9L5U9hx1Iqtlvd6pogfEjpPRBKEBFs8L07IS2zQ6rll1DtNqNeSbGZw87zL1xrYL84PK5UK58WD2QVliGV9GfJwqv2iQ/jQ1p/pR7Ftvqn29iIiI9FLez7XxAOeT6+6ump/VmzeIlFBFRERERERkCIxkz10rssLo4mmd6LlqpHL7ZXXMSM8f1qnVS9ROWGje8MlW5T3WecMosxKU1EtakjeMsuje21rhu3HvWtwjVnlsytaP1qlMbFLWE1nn2bjScYh6uGX4nXbER0uvb9rwmR7WREREJNHoUakixletl6eiX8dvHcnG3dSc6dLrB/7uOB7+X39ceh+fqLBU/GE/hOjlbUTEH7bzPIcWjwVW74JqtP3wYX26iXXyCGOg1Rp0vLKeWdOysirGjZVS3f2xyKp3fFyzni/MCrMMys6x30YzyUvqNqCi7VcuX3rmLWpItRKq2apw7hplrKwc5y5LPH5h1j5kHbui90c6a2LBOAAL79zS45qIiIi0L35+LmTIjJ+vY/tubnrgEt584sXJvO+t72r9iqKwTBERERERkSEwkj13leqFKmapDF8M4h6RuHemmdT7cTlZWQ3rZd+slWSkiB67ym3VUrmtyuNaRLf1CW+9JFcYYd4MmkUen6yyujm2W9h+kePMNeqlFWlXJ8atFBERiTUaCqG0zJuOKXRb0N3/ayPZuBtbNFF6veqv91bNrwwfzPqgHIdP5lUr3DIuJ/7wXS+rYZyZMCtzY73w0ThksTJUMGu/ag1OXbk/eT/4x2VkZWasPD5ZIZix0v630PAIdc7KiFkvtDOWN9w0hI1OR+s2Ot5xHWNztjxT1WitF+IJ6TUTh682GypZKiMqN2SYjbPLNqoLdza1Wemi8eeSK3TfvOTfw+TGh3tYGxERkWz1HgFqxPlsmZNUf56tFY5ZWXa/fhGpsEwREREREZEhMJI9d0Ur9TTVCA/MSqQRq5dpsl64ZZ56NQzVzMiqGcYwC0kzapXRyph7za5TRBbOo3//UgAObLBcvQyaQa1kJ+18e1MZZpr3fHdioPNY0YOei4iIiBSlXrbK01alkW+TDz5Ut5x42Tzbq/eZr9OfzfIYycbd1HPpbk/fdRdQP+NfpVontbLxVi8zYa0MlnlC5irDLafrLJ8n7A9g17tem7mtShMZg1MH8QDwWQ3DrPrVGww8q55ZoY31tv/U772uarm8GRzbye4YjlOt8xn2Mev6qMw0mRXGWuuLglplQdpIizOoZmUarRQfrxAOGp+bOByzXmbQsmtMYZl9a+bNyTnb77CVQBKy0mn9GtoiIiLDIx7EfO+KAwCY/cttVfPDvHazZfbqf9tINu4aydPqzpOE46SxMxr25s1pssxWVD472K9JMSobBHnHoOsHzZ67biZZacUpx1xYc6zEE956Sdl1G2vUiy0iIiIyyCp7+hr1DOZV1Lh5I9m4G5+Tfhdd2Vhodyyuyocys3pY4iQblQ2uOFEKGb0/QeVy3//2eaXwQ4AFm/bVXDfex12+Z+vAL/ywrB5AKUFJvHxIpJHVczZ7e3pcQzKYuNETEpRkNZ4nMhoDYR/jY1iaFi0X6hxvP1iybkfV/mQdu7ielduo1dPX7PiBWec4K6tqZSP38Q+mvY9LL0vOU9wLFvZ77+L0z3nBXdl1qLxZhOsjrlOYVuolzeg5zApLjc/rky+fCcDyLemx27PEsislfWX6hKMBeOqFswFYqIQqIiLSh+o1gLIaXJMPPlT1OT3+5Djx1mNLr7///3607HPNlM+gGff0tVKvblBCFRERERERkSEwkj137aqXqKQbD1JmhXueNHYGZDxf1m86FXrajuN/57MsbLxYlaL3JSSymV1oqeVaTd7SyPG/81kWFFJSfmZ2CHA1sBRwwOXOucvM7ADgq8BK4GHgTOfcdjMz4DLgLcAe4L3OuTt9WWuAP/NFf9I5d1U390VERKTT+n1M0aLCEus5+diLmq5LKzTOXZfNmDVVer119SIgDWPMCvdrFKqZleQiK8wwK0QyJPwI289aftNJi0rTXnCJD8tbPrM0LYRUHviFH5YuntBQmHVT9phyQdhuvQ/5cVjgwT/YUbV81vhpeeVJEBJvKywX70scjliSUW4oJ2v5+Bzv9SGP2/1+h1DIeLuNwl1L11HGeHAAt33pw2XrhvMVC8d98S/SayMrQczC+5O6xOc6yArfjFVe/7GsMe3Cscsa5y6+ThZscmVlVO5HQSaBjzjn7jSzBcB6M1sLvBe4xTl3sZmdD5wPfAw4DVjlf44DPgcc5xuDFwKrSRqJ683sBufc9qIrPAhm7E7GAZ29M7nHPHf6cfzHN/64l1USERFpSlZjKqtht+c3jyu9nvv128vmxQlYZvxyGzc9cElHG21FlaWwTBEZSM65x0LPm3NuF3A/sAw4HQg9b1cBv+5fnw5c7RK3AYvM7GDgFGCtc26bb9CtBU7t3p6IiIiIFGMke+6GWVYPUK+84kOXcnCvK9Gn4pDOdsbwk4SZrQReDdwOLHXOPeZnbSEJ24Sk4fdotNomP63WdBEREemByjHluhGymbdO/W4kG3cTz6TBl7OeScLHQpjjrCiMbWKez+63+nVVy0N1aF38gT0rnX8pbC8KKazcfgixBEohfVPRQ1ihvDgL5FYfqhfXLYw99ly0razMiPWEcLsQihlUZryckxHS+dgbFlVtKyuUMeu4h7qH5R9526LSvKXPT/anLDOkH9Puv89Lz1MIVdz6e9nnLgihivF5CtNm+6C8OBQyHPdGobqhvDgUE5JrJutZvayQycrjDmnoZbyvpWNcI1Q1XKchC2dWeG4s1Ll0DUXnMBzr+LiGkM6JeWmocNb+hAyaANxYtwpNMbP5wL8Cf+Scezp5tC7hnHNmVn3iW9/WOcA5ACtWrCiq2L4yvf4eAKYOPb7HNRERkSL12zN2lbIackVy6+6umjb367dXbTcOx7zpgUuq6tTvx1FhmSIysMxsJknD7svOua/7yY/7cEv87yf89M3AIdHqy/20WtOrOOcud86tds6tXrJkSXE7IiIiIlKAkey5K0poxYekFa1kXBxEJ7z1krLEM5JfVmbO43/nsx3NkDmsA4v77JdXAPc75+J45BuANcDF/vc3o+nvN7NrSRKq7HTOPWZmNwN/aWaL/XInAxd0Yx8GxSB9YykiIsOt8n9SJ5Oc1OtFzPrfmHf5WL11W9mXkWzcHbBkV+n17O37A7DgK9VhlPWyQGaF5cXTamXarDxJYeDxEI4Zh8yF0M4XUh3aGS8XQgXj8LmQBbFsEO+KcmPPZWX8XF2dmTFsIyu7ZryvSzMG0S6FiGbUHdJnzyoH8Y73P4QMhvDA2EF3p9kYQz3jLJD1MnLGWS3D+c7KJJk1iHn4Y8wK34zrGb4EiM9TVubOMG3v4mRf4xDHEAK7fO2O2vtCdubM7Gym1SG6IRwzK4w4DE4fzlFcbqhvrW3te+2iunVuweuBdwN3m9lP/LQ/IWnUXWdmZwOPAGf6eTeSDIOwgWQohPcBOOe2mdkngHCBXOScazxK6ZCacdjKsvfzvnZbbyoiIiIjKatBU0TIZhxuCWnIZWCrXw7AZEb4ZuW6/WwkG3ciRQjJa7o9vluvnXLMhZDzuc1Ocs79J2A1Zp+YsbwDzq1R1pXAlcXVTkREZPj0+1h5UDsZS2UD7eRjLyp7Di807vKW3w2tbGskG3fbtqYfx2eFnpOM8ctCD8qu5YvS5X1SjtnPT3szQs/KREbykkZKCTpC+Vm9avEHad9zUp4VM53//W+fB6QNj8qEHgALol610JsUepjiHrnQYxX3Vj3pE3ks35KWUUrAEfWwhSQoh3wiShATthmN0Rf39gThGIRx02bsieqeMS5fVoKYUEY8zlsIJY2nZfUAVvb6xecynKe4jCDuhSwd95dWLxeXN2dL0mMZJ1mZs6W8HpV1+/63z8tMyhKS2AAcHJ2L0CMabhBZ4xbG26rsnczqpSu7Tv25KDuWGT27h35jR+n1A1VzpV9ML57f6yqIiIhkqterB/UbQ1MN/r+Nb98NJIPoVpp88KHMENBadWklxLJyfquNSCVUERERERERGQIj2XPXS6ccc2H58AAtlkFGjxykz6uR8SxXvxuU8UOKdsJbL2nqOmh2+axevqbW7YMQTBEREZFmwkJ7FT7Z61DVwWsBFOGZ8dLLrHDEyjC2rLHqsqaVh+pVjz2WJzFL1hhkcQhc2O7C+3dW1TNed/b2Sb7/7fPKLrYwP14uJM8gY/+zxkULSUvK6v3S6kQhIRwzXjdrLL/KsNB42tytrmpeVqKaECr6gkuqQ0CzZIVixuZseYab13+cV3woSXazmPKwWIiOG5SOXdaYdnH4bAjfDYlaoDxZSxD2LYSRZl0TWdsqu4llhEWGhC/huMbbqkyaAun1Go4HZDcUw3mflTEtDsHNSu4j/Wu/rRNAdYIVEREZXb1uuDRSKyyyVr1PW3Ve2fs4yUrlvFrldKoRGW8rHsO3EYVlioiIiIiIDIG2eu7M7GFgFzAFTDrnVpvZAcBXgZXAw8CZzrntfkyqy0hSke8B3uucu7Od7XdKO2Fs3dSonq2OR1cqd4TD8U546yXwopmNF2xURpfCYzt9zYYeyFEZy1FERIZbP4XRdUvcEzX54EOl150aF67Vcnt5bk4aO6PmsAeVWTjrldHquo3KyaOIT55vcs49Gb0/H7jFOXexmZ3v338MOA1Y5X+OAz7nf3fd7Oen6RfjkLfKLIRZoZphftlzcxmhkqUPw1H4XlYIXMhcGLIWxqGbpfC9jPHTYqVB1KNtZY1Ht903VsK2IA3vCyGDWaGj8VhsWWGc8X4HYf8n5pV3I9/5+Q9x/DOOSnFIa9hGVrbIkIVzbjQGX2nMt6hOpQyeOcMn44yQweJfVGfLzKpTEI9HF2Sd/yUVIZOVxy5cYyH7ZahHbM+S9LiG85kVArvn+ek6S9cl5cTnMyvks/IaiI9Xaf8z1hvLGCswHm9v03nRuIl9+bWOAOw9aD8AnjswuU7Gvv9wD2sjIiLSvLyNqDgMM6sxFea30tAqYmy+VnSiW+F04I3+9VXArSSNu9OBq/1YU7eZ2SIzO9g591gH6tCSk8bOyBycG9LnzxolsiidwIxU+cPspLEzyhsKNXqsyodw6A9H//6lPR2r7hUfurTsOTj1jomIiPSnQRjnrVXN9pgV3UvXTjl5ttFoEPNWy21lfqfqBO037hzwHTNzwOedc5cDS6MG2xZgqX+9DHg0WneTn1bWuDOzc4BzAFasWNFm9bLtfXK/0us5W5IxLUKvSmWCk8pxvip7uiDtYVl6WZrQIyTtmCYVeuni5RbzmrIypmany4/vTX7H47xlbT/0DsU9LGmPWdrTFMqJlwvlhN6ceNy5MC/ufQq9iXEikKyx+bLqFOxanvY6hQQtceKNh34j2cYBP0saPHOiBndW71PoTdr8hrSMw6/bUbXdsFzW/sTbn1XRszgnGtMvLBcfp1De9iiMM/S2xfsfxhfcdNKi0rRwjmNhndDge3z1zKrlw3GD9NqJe2TD2Ixx+XszGtzh3D0V9Q5XjrMX976G/V+wKT1GYf7jUaKYWU8nv2dEPawyGPZ7YKv/7ScooYqIyFAoouHQDb1qMBcxRl0IfY3DXvOsn2XqTcekb777tdzrtZtQ5Vecc0eThFyea2ZviGf6XrrqGLw6nHOXO+dWO+dWL1mypM3qiYiIiHSXmR1iZt8zs/vM7F4z+6CffoCZrTWzB/zvxX66mdnfmtkGM/upmR0dlbXGL/+Ama3p1T6JyGBoq+fOObfZ/37CzL4BHAs8HsItzexg4Am/+GbgkGj15X5aX3vFhy5lca8rISIiIoNkEviIc+5OM1sArDeztcB7aSIvgU9SdyGwmuTL8vVmdoNzbnvX90gGSqNwx0FIKFNvmII863Szp7JWmGXl/G4c65Ybd2Y2Dxhzzu3yr08GLgJuANYAF/vf3/Sr3AC838yuJblx7ezV83bjCydKr0NI21w/Vlsc7hjC7OIEEyFBRPxsXhgHLSs5RRwKF4djBk++vDwj4zMrptJtbUvG44uf81uwabKsbgATPhxubhQ+GMQhfaUy1qX7mIbZ7cuozwy/zKLSlMm55esBTGaEmz7nj0VlkpHvf/u8zOfu4sQfh35jR1KGDwvMCiPNsv/G9HUIfcwKH9z+onR/QuhjHHpYGaoaj0sXEtqExC4Ae5Ykryf2T7cfrpn42D1yUXKcZu6urvuOw9Pth/MYwkPj0Mowlt+j/zsNd8wK7dztv0aJtxX2Kw5fXeiv4zgUNfxNhP2Jr6vZ/uNEfLyykvyEY7Y7+jpnam5TnfjSI9OL55e/X39Pj2oiMrj855vH/OtdZnY/yaMoTeUl8Muudc5tA/ANxFOBa7q2M0Oq1Q/Z/doQKkKz+5a1fBFlFLl8q+s0U27cUJ588KG6Del661aasWui5rx62um5Wwp8ww+qNwP4inPu38zsDuA6MzsbeAQ40y9/I8kwCBtIhkJ4XxvbFhEREel7ZrYSeDVwO83nJag1XUQkU8uNO+fcRuCVGdOfAk7MmO6Ac1vdngyHl1x4KYt6XQkREZEuMLP5wL8Cf+Sce9p/IQ4kn4t8QroittPxZHRSrdehjZ3afr1yiwr3bLXuzWTXDMva6pcD4NbdXbVuM/vTSn0r1Uq4krXNUO9mdWeE5T5jD88tva4My5sdRbGHeXFY3iwfvhmH0YXMgHGmyzD2WhyWlzUeXbD7RUlY5PxfpCFzISPiphPT0zRjTzL/uQPTdec8VVVcScguCWmIZlZmyBDmGYf4hTC++JiEdSfTQ1ja75ANFLJD+kLoZXzsxvdWbzeEY4b9nvdouvzjxyWhqgs2ptMm5iVl7HxxNJzAz83XN10uDdFMp217SfJ6OopGXbDJl+HP04bfjjN+Lqqqb7hOpman5cZhm0EIkSyNyxct98JvpdNCWOTk3PIxAiENx4yPfxCPwTfff88bLxfO8Vg0bN6DZybbWvaDdGIIzQ3ZVeMynn5pck3sf396TYYQ2D3L0tywM/0lXh6CWr0/0n8m588CYNtLk4Dw529f2cPaiAwuM5tJ0rD7snPu635ys3kJNpOGcYbpt1Zuy2crvxxg9erVioGX3Nkd+1m9BlAt7YxLFzQ7Pl2zmTSbXb9Z7WbLFBEREZGIJV10VwD3O+fiB81DXgKozkvwHp8183jSvAQ3Ayeb2WKfWfNkP01EJNNI9tyJiIiIdNDrgXcDd5vZT/y0PyFJNpc7L4FzbpuZfQIIA8peFJKrSDFq9c4025vS6TDDdnqiOqmbWSCzjkE7YZN5y6t08rEXlV5/58d/XrZOnC2zV+MKjmTjbvIFaUzdni2zy+ZlhSzujcZCCKFyIRQT0sGj41C8EO4WZ3AMYX5x+Fwwc2t1VssQHjceDWIe6hKH1pUGJ4/CEm1fCIFLyw0hmvE+bvUj6Sy5M1k3DpkMYXxxCOjEAUk2z3m/HC9NC+GOT748epZgZjJt8uk4BLP8d1z3OKtiyFIZwjHjTJpzN49V1ekAHxa55M502p6MIRJDiGa8/YnnJWGGM7anfwph+788NYktXfFv1RlUy0NVq8MNw/UR1z1sNx7EfNGDSd3ja2ffQj94+xPVocJxeUEIKY3Dd0NY5HT1ZcW+hWn45Ownk+MZZ1+dXJCc48Ov21VV3zmbfQbVw9IyFmxMylh8TxoIEP/NBMqWORgmFiXXwvzHkutgcuPDPayNyGByzv0ntWPRm8pL4Jy7EriyuNpJUZoN3+tkHbpVZjyvnYZVnuVqbSvPcW9UfrPlheXihl2z4jDTWkMmxPYu3a+l7Yxk405EREREBkc/JS8J8nxAr1dGMw2QZsttdltFa7VHsp+fBTxp7IyqJCe1nqOrXK+d+c0aycbdgoXPll5P7J/03IXEH3GPw84XJ7/nRUmIF/8i6jLzQq9L6IWBtJcuJMeA7LHXQm9O6C2ZuyXtESsldDk63dZ835sVj0u2Z0kyLe2ty96f0IsXJ8gIvV0hsUjcu7LoQV/+89NyZ+xK6xfEY54Fk4uS3q6JnWmPUFYvUki8csDP0u2G3sOwfDg2APtvTLYfJ68J24+PSegRjMeeCz2X8biFc7cmyVK2vSRdLmw/7rGrFB/XcD7jXs8wP647VB+72duTdeJrJ/QEhh67MAYhwOTc5KCE3j2AsT3J8gdEPZfx/gQh8c7YvrSHLe6BDkKvbHztBuEYxuMnhnMRn5N988u3CTA1Pz4W0q/+4xt/DMCJb/o/AEyfcDS3fO+CXlZJRKRj+rkxUSmr0dZL3Tx2RYyp1+w2so61W3d3VYKUGYcf2nQIZhgPr15Y8LFrqseFzkMJVURERERERIbASPbciYhItted+dcAxJH+vQ6HEhGpp6ikKFnltVJGZX2aDd/shjzHLCtRSOXwBEX/T2jmObp6obLN1LOI/3Gt9KjWW+e0VefB6zISLeQwko27p59ME5rY0iRUbNbTSShaCCeDdKy0OAQvJJ6Iw/1CaFsIbYR0fLs4eUcIW4sTj4TwwVnbqkP2wnILf169DyEUE9LkIjP2VC/3zIo0FG7f/GQb86Mw0xBKGJJsjEWhnZXJVuI6xSF4216ZbOOF356OlvPhg9HxPPhH1SF9WUlewrTHXptM2/+B9NiEpDHje6pDQeNjEsTj7JWShrxoUbqty34IwMTvpWP0pYl0Zvh9qQ53jfcrSxg3LuwDpPsf72sIuYyPSbBk3Y6kjCg8MpzjeP/DGHlx4pOwrThRSlg3vk5CgpZ4PMTKRDpxyGyYF4eRhrDY+O8kXGNxkh9mp9dHUczsSuBtwBPOuaP8tAOArwIrgYeBM51z231q8stIMtLtAd7rnLvTr7MG+DNf7Cedc1cVXtkB88yy5I98wVd+1OOaiIj05xdLIbQuVu8De1b4XrP7lTc8sdfHq5XEKt0MOc2TRGXqTceUvf/uLeeXvW9mTL2109eXPZ8XNz4r6xIaqOP7WktEp7BMERlkXwROrZh2PnCLc24VcIt/D3AasMr/nAN8DkqNwQuB44BjgQv9eFIiIiIiA2Uke+5EZDg4535gZisrJp8OvNG/vgq4FfiYn361Tzl+m5ktMrOD/bJrw9hRZraWpMF4Tafr3y/KvnV8x/G9q4iISAc0k/a+lz1eRYQHxr1Dlb2EneoZa6bcdo9vq/tQOR5fPyWmKdrIN+5C9scQqjY/Cll7+rDkdxzaF8IRs8It4/C1rFC17AyWiRDaGYd7huXjsc1CBsl4rLIXfbE65WEI5YvHowuhhGXjxvltvPBbSRk7X7qwunKRg+6u7iJe+K1kPLR4/L6Q6TEO3wxhjnGGxpC5MoSxxtMO/UbyPg43DGPfxeWG8MWtq9PlQmhlPAbdQXdXh0U++r+TcMz4HM/ePlm2XDymYMiWuvkNM6uWn/V0Oi2EpcbZIkN5cbbScEzi62n52vL9CeMoQnp+pmem+xXOdZyFNYRPxuuG5eKwzHA9xRlEZ/hzMbZ/qHe6rbCNrBDYeB/CNb78ljTj6COndy1QYKlz7jH/eguw1L9eBkRByWzy02pNH2lznkrO3dgxR/W4JiIi9XW7QZZ3HLhObqddeQbb7mVDt93x8opocI5/b33T69droFc+B1irnFDGlJ65ExEp55xzZlbY6Olmdg5JSCcrVqwoqliRpgzSuFAigy6rkVD5Ab5yrLM8afFb7aXL2lZW4pNa26pXZtZzYL1S2dPW7nJ5ysiT+KZyEPPv/PjPm16n00aycWfj1eOhhd6suFdpam7SO3bwD9KJoXcq9LgAbHldstz43rSXLPQshd4iSHtuQq8WwIFfSBJ6/Pd5r/NlpPPCOHfxOHOHX7cDKE+yEXpz4t6UsP24NyckT4l7qeJeJEjH54PynrBKce/jxLxkudCDBvC4XzcrUUfc6xa2EdcjjBt4yCeSYzP+hjTZSTgmmWPrxeP3rQv7kfZEZvWchWMRemkBFmwqLy/ufQsJUuIezFBuqBukvXnxdRLOTzxWYujhi89TOJ+h9/HBMxdV1TdcB5Bek3HylNADHPfEhvrFPXEh8UnWmHYhaUp8Xe1ZYvz00g9x9O9fWpoWjlO8X7NuugOAxz+Ynrvx3YW1sRp53MwOds495sMun/DTNwOHRMst99M2k4Zxhum3ZhXsnLscuBxg9erVXduhbtqzNOlh3XFE8vuF123tZXVERNrWbMKRWqF7/TbOXJa8vUNBnkG4m9UPx6myoVYrkU3R9WymkVlr2TBQ+rMHVn/WzUMJVURk2NwArPGv1wDfjKa/xxLHAzt9+ObNwMlmttgnUjnZTxMREREZKCPZcycirTvlmAuhTq9uN5nZNSS9bgeZ2SaSrJcXA9eZ2dnAI8CZfvEbSYZB2EAyFML7AJxz28zsE8AdfrmLQnIVERGRZmX1BrUSNthKCGCjejSjVur+IuWtY7uhjSeNnVHqEQMY37677TIgPUb1ktdUrtNo3XaNZONufE6a5GF65iwgTZoRJ6CYuzUJaYtD1uLQuyCM7/bky9OwzBCO2WhMtxC29sxhSZ3mbK4O48savy6eFsIt47HHwrTYogcdt33pw8mHcy+E7WWFm4aQ0YnTXlOaFuZPzq0uP06oEuoSlxfEoadhnLd4uanZVrbdOMnMxP7JvDC2G6QJV+LQxum77gNgTxTSGcJB47DYEF4aTwvnJ4QZxucrJNyJw1ezktCE/QnhmfE24uQlCxYv9NuoTloSpsUhoGG78TbDcvHYeyHcMg6pDMcxDj2uDMuFKFQ5I4w01D1cGwBP+TECy8bUW/66srIgHTeySM65d9aYdWLGsg44t0Y5VwJXFli1gfXskuS387ezyY0P96wukq3Xz8KIjLKs57LiD+ftNKjy/G230hCo9yxZEY22rPIbhUU2KifvftY6ZnnPSdbYg/Wyq7p1d2fWtx637u6mwmHDNma85nU1l6lHYZkiIiIiIiJDYCR77kbZCW+9hDm9roSIiIhID3U74Ue97dWalxUCWFR4ZKNEKieNnZG7Z6qVdRsd/3azZRZ1fotKOFN57lrJaprXSDbupp5K0zUu9OFrIdNjHFoYMjLGYWkh3G7BV35UmhbC0uJww5At8IVbjixNq8yCCGl4YwjH3HtQOn7d+N6xqnJDdsl4rLoQ7hkyL8bT4pDCEAKZlQVzz5LyUMR4v8qW89uNQ0BDqGB87EKoYBwqufSyH1aVm3U89/hQ1TCWXJzdM4i3FcaFC6GYALve9VqgPCww7Ft8nMKxjcNxH/qNRX5fk+XiMMYQ2hiPvRfE+xo88ta0c/yA/0ri3PYuTkMqQ9jmxLy0Tmk4ZnVo77aXVG83hBTPjMLHtx+VXEeL70m3Xy/TaJZwTOKw5DAG3qyXvrY0LYS2xhlUS2HOP0qvp/hYSP/atyC5TqbnJNfQMxrUXERGWGWjoZmwyDwDptf7QB+H8zVqZBQVrh3XafLBh5ouN163mW3G62aJwyGbzXRZakhFz9nVq1+9Y5433LRR47bRsV07fT2v/MClNefXo09bIiIiIiIiQ2Ake+5EREREpL9ljTPX7Drxuq2EvdUbbLzVgciz1m+0TLthjnnrV1lO5T6++cSLS+/Hv7c+d/n1eh5bCVmtVX47YatFDSxfhNNWnQenHdzSuqPZuJuRhs+FULJ6A3aHEEuAOa9Owix3visNSwthlvFg4yHTY5wZMYTgxdNCmN3UqiR+0KbSMmbfsx9QPsB2CCOMQwXj0MsgDLY9Hc3a+eLkdB/2x2mmwzB4+p5lIRw0XSGE5cXHJs5SGYQQxTjjZFq3dH8W+mMXh6WGY7IrOp4hfDBkzYxDK3f7Iajj7JohzHFi9euqlosHIA/rxGGuoew49DCEFMahp0HI1hnva9hW1gDws7alGVSDOJPp9MzkeMchlSGrZ6hvXI8QUhkfw20vSabFIb2zn0w65eNsmNtemcxf9f7bS9M2/tVrfT2q9/Xpw5J9mFwwVZoWQkvjLJzheMbnKdQ5XIfJAmnmUOkvMw5bWXo9fdAEADN9VuHnFs/LWkVEZCjFH+qLegYqKxV+KSNi1Gir1zAJYXxxncK6rTRoQjhhreEOWm0M5w0jDVkqm9lGrQZXK9k68wytkHVc8x7ryuXyNBYrj138uEszRrNxJyIiIiJ9pdmek7y9Ws04+diLyp7vqrXNPNNLDZzF86vmdVtlYys0Poo6du32YubR6NzUkzVGXZDVEC26Ny4ryUzR49sFI9+4C8k4QoKOuFctJOgIvVvx8lkO+UTaIxZ6ouZEyUBCD1s8bfKwZwF42QuSiT9//HmleTtf7HuwlkykK2xJunWeXpX2poRx9hZsShcLSTgO/UZa39A79ej/TvcnJF7ZvjdZPk7UEpafSvPPZPbShUQi8ZhuoQcw7mEK5cXj1oXl4t6p0CsUthHXKfSqbT2aKkvuTMs44GfJ77jnMLx+7nnpcnOeSKbFvXnhddZYgaEnNK5T6HV78MxFpWmh7jMzxsjMStASb79S3EtWmhb1EofeyUO/kRYczlPcm7b/A0mvW3w9T8/0187G9PHbcL7DuHlx7+O2V06VlRVvIz6v4VjHPacT6aUtfebZVUtKr5cf/CQA82Ym950HD1fPnYj0RjuNhLyNllof+PPUIavnLGu7WclA8sjqXYrrO0UysHmj3rxWx7nLqut3bzm/alpWw7ZWuZVl1qpb3p62tdPXl80L5e1dcUBpWhw+WppWZxDzRtdO5T5k9dK10+sb9v3lH1ZCFRERERERkZE18j13IiL9KG+ISzdCYUREitbsvSsOk2slnC3PM1b1tvv0K5Pwk7n1Fm5BvR6eRr2KjdZvV1ZCmzxJburVKRzPKR+qWp2ZIC1j6k3HlN7XarCctuq8roe9hn3Y85Ik6mVWxfxmksBk9XCWeg7froQquc16It3tR962CEjDE+PkIU+fmYSvHXR3dbKJOAQxJFx5/IPV48JNzEvLC2Ppxaanks7T8w65CYA7l6wszbtq7VsB2HZgGoIXQubiULknXz7u6xSV68MH4/0J4XMhYQnAphOTY7H8lmT/D/5BGpb6i/cmcXZzN6cdvFnhgyEEL4yVB2nY4rIfxKGa4bin+5M19tpzBya/w2Dr8TbD/odwSkjDB3ccni4X1olDSvctTLY1c2e6bnhYdf+N1fsV1o3H2Qvjx2WFMcayjlM4//E5CccpTnwSQk+DcDwA5vtxGeNwy1BGGJ8PYGpuCJWNkvxkhICGxCtxSGUIkX38uOR32b7OHauuk69LnLznuWXJeX/+D9LrdOfLq/+OpD9se0l68l63+DEAlsxKQlZ+dsCyntRJREZPs19QhcbPZMVzWGunry9rDIaGRFa4XbxuKaSx4n0s67mpWs9nhUZoqw2wrNDB8Yyx2uJG4GRGopa4LpVqJVRpV96w2KmMhlnWc3U7Xpost3Bd/W2EBtfcn20tTdvrG4lTc9PPJHMzju3us5JxXRdseKZqXjjGWeckKylOPC3sYwijhXzP9E1nfL7MQ2GZIiIiIiIiQ2Ake+5ERAZZJ8NwREQ6pVZvRa0QzaxEIWXl1AjHe/07/or9M+cU780nXkytDpasUNBWw0OhvQyOYbt5e46y6vnmEy+uCqM84e2fqQpLlNrCca0VjlqEkWzchZA1gBf+8Y8AeOr3kpDKBZvSmLg4++T3v538EZxyzIVA+dh3wdLLflj2vjK8IKwbMnMCLFiYdP2+fk7SiXrIjHtL876+6eTkxZJdpWl7NibxcHHoXshWGRvbZ1X7E7J1xuOxzd5eniU0Dhmcuzn5fdDdFXGClI8zd+AXkmMYxvYDmLs1mf/ky9Nwr1DPOKRvxp7qes56Jlk3ZImM9y+Eb8bLB/G2gnhcvrBvcQhkCEGMM2OGMMwQMhpnq1x8T3Kewth2kGaEjMM3Q5hjCAWFNMx1el5a9wNvT+ocslACzHk8+ZMP+xpfryFLZrwPcUjpnZ//EABHfPrSmvuaNUbdM4ek2whhtiF8Nc4MuvDnfky9aOyVMA5jHCqahj6n5dpzChRoRjPhSe0+a/dcmiyT/+eAJOZljiXX6TeXNn7mQ0SkWc2Oc1Z6RisKixvfvruqjD2/eVzy+yVLqDS5Iv2nOGOvf3whaiDuXeCbKT7T4p7os868Xyb/5Ha8Zmlp2n6PV38WyQozDM+Ozdg1UTVv4q3HpvuzZ6p8uToZHWNu3d1VY6TFIZZ5nt2DtAEd9qEs3NMfm7iRHY5xvFwIqYyP/+zHk8zwIbQytuj+8vMJ8LQ/h9Oz0s8r8zclxyQ+XlP+eMXhlpV1gzREc+sb0g8043uSczz7l9tK08L5DPu64/ei7PLfTR5ZiLNwTixIPteMR9NCQz8+/jvf/dqqfS09r5cRChu2MbGwalYuI9m4a0W4cYz5gbibXa+VdUWadcJbfUz9G6obujKY8nzwCcsooYqI9EJlL1CtRlvlM1K1xhfLaiDVsvrsz9LiWM9D5fXv+CuAhj2Wx675LAAHNlhOBtdINu7G90RJJnxv04FfSHrd4gZYGOculjUtr7Bu3H09/fLkW4FzXpD8vn97+o3QfN+b9uQP0j/BF/px9na+NG3Oh56YuOdwtt+vuIdxid+3uOcw9NiFuh14V1q3MFZfVi9lvA/hmMW9eaF3MPQMxttaeH9VceW9mV9JegLx2y+9Jx4/MH3YNRyLrDEIy8r1vX3xugA3r/84R/9+OpZIuBZK61VXlwVRL2XY77hHNIivl9A7vGTd01XLzXomPZ8LvpJsP+zrEV+uLjfUG6IGXeTw63YA5ftfOo++3NjU7PRvIu4pve1LH+b43/lsVLcfldUN4KC7J8vLp7wXN/R6x/V8JHOPpBmhIVdEiObEAWnP8XGzk7+P+WNJSqNjD/5l2+WLiEDSm1HZw1Qv8UacLTHYdUT6OWnffsn/ru3vS/8nzXsiuZ8987y0N2f2ziQZ2rMHptP2eypZLg6pfHZp8sXoc4uTHpk526ej7Sb/T/cuSKNQ9u3nPw2tTHuE9n846WHa/pK05H2+w3DOtvoBjFNzks8TM3cny43vS/+HP7M02e509Ml95p4km+Lcx9N7eKlnMephnLXL7/8B6f6H7J9Z9s2rjrSJe51KdfLHeNcL0k9KYy9PzsVYlE9vYsE8vz/p54twjJ+JelP3+Z69ffOrx/cdm5hRVbdnjkqmzU4735jxXLKN+DjNfnx+VZ2efEVyfhYsjsaXPjTZnzAe8d5F6fL/fVpyrGdG4zdPzUnqufv56bme7a+FsddFYU/etpenX1rMfDYp5+n/WZ2McebuZJ4dka/ntpLipERERERERIaAGncy8gY1OcVJY2d0vO6DemxERKR36iXrePOJF/PmEy/uYm1ERstIhmXuW5h2s4dQsjgcM4S7BY0GbKz3AbjRMzAh5O17P3gFUD5+246T4N6LP5QZdheHAC74Shr6F7YXkrfw6iOr9qc0j/phpnE4ZJ5neeJygzhUb9qXU7k/IWwvCMcza/vheO186UJu+9KHM9cr237F+1rnrjIUc+309aX9qTxGtcqYhlzLPReFLIbjs+CuskXLBgoNxy2rrFBG5cCioS6zqA7fi6+dsG9ZoadBWbipv54aXUNZobzSmhmHrSy9ntz4cNX8op+zu39fEpZyyIwkHGT7RNFD9orIsInHUKsXbhknowgJJeJQwSCE7+1alvZBzPRJuyaix/HGfd6ReKzWZw9KPto+e3AaPjfvl8k2pvaL6rIombZjVfofcF946UK5Y1XzpmZFSez8R7b4s9uulUmI3vhz6WITi3xCuSgScvy56tDDSZ9Azfw/cWfpMiGkME6ot9c/tfPMsvQYhvnx2LPT48l+xInPxifGy/YLYMwnWQvL7YuO6775yYI7V6UnwHydZkUfE9yMpM775qXHab/Hx/y8dLmw7pbXpsd49rZk3fFnfb2jmNnn/HUyFuWkmfIf8uLwyT1zQ1K4dFoIh9y1Mkoo58/PUy9Lj91Mn7/w2eeVLwPpmMZ7o7Gnwz64qKvMjfv9b/Dw48ynqxPkTc5PjllIinjM8kdL8x6oX1wZ9dyJiIiIiIgMgZHsuZPm1RqDpsiy65Vb2n5GMpBmZfUwNiNvqOKwZjA8aewMZX7tgtVnJ73U1cm8y3Xyb1NEJMi615zw9s8AsOPQ9OPkC7pbLRGpMJKNuxd9McpWWCccM2j0gamdD1Qh9HDH4Um2nDDuWCyMT5e8XsSdn/9Qw1DRWvsC5WF0leGoWeXmbcw0Cl8M4jDMrLLrhSAu+MqPch/vrBDQSvXCccP7PPufJ1S3Xr2bbXBmlZX3Om3mObpa5z/vuY4pVLN9Mw5byU0bPlN4ufv/PA1L+e4bXgrAQh8Xc++W6oxfIiK1bH5bklUwDpUM4XVTc9JQvbGJJPRsekY6LYTtOX9LGo/CCPf5cLznlqVxieO7fajes2kg2sRiny1xdvqwxO5DrWybkIY+Ts5Plxv35YQwu8kFUQjmlFWV6+Ykr/dOpp/Tpn3Y5t4l01Xruv3SrJaTz/qdjPbf/Di403uTebY33a8pv7ibmZZbmh9Vc2zC78N4dFznJ9t1T6Uf+59b4MtJq8T4nmTdUgjq/unMiUXJPIsPyZIkLnH62fQhmH0+tNDNTBd89mBf32hb4Vy4KCp374HJOjOeqQ5ZLe1LFHMYwmxDOCOkIY27D42O/0PJSpPRuMH7fHTp+N50W/v29/MtlJ+WMePpcb+tdNrMXb7ceem0SR9ua1NRuYv8jkdV2uezrsbXE/56Yiypx/GL0oGMryW/rodlmtmpZvZzM9tgZud3e/siIll0bxKRfqX7k4jk1dWeOzMbB/4eOAnYBNxhZjc45zKzejz27AY+de/bAPjTl32ra/UMavVwtNpT1+nwqcren0aJYYoySBkVmwkpLHq/4vNTr2e13rbzXkP9Eqo3KNdGs/cmgMef+wWX3n8yAB966Xe6Uk8Y3nBfEcnW9P1p3z1Mb3kRAGPP/0VH6vSHd/62f/WqjpQvIq3rdljmscAG59xGADO7FjgdaH1k8BY8eOaisvcbPvahbm6+TPiAdvhnk+drFj1YvUxlJkc+/6G2P9jVWr+VcL9Wl21l+WbLygoFnL7rvra220yoal5xY6/TDaJ43/M0Ajtxjsxqh1z0SF/cm2LbXun87+ez8QMfAZoP981r4cNprMwPnzocgANmJ/FQE4/Oz1xHRLqmP+5P310OwEv2f6I0adZH/huAPVsPLE0bH09Cy57bnYbqjflQwvjWPznpw9fG0lC56aeTdZwPVSwLbTvIpyacTqeFkM44DG3mLp+tMJoapk0sSkPgZm9L5o/tq87WOeOZ5PfeNLlnKYxzxrPp8uP+9XSUmtucz/i4d7xqXXakdZr02SRdFJXnfIpLC/sYHRv88Zr5dPrRfWo/H4IahXtOh3DIfdHBngghgHGmTx8CuSsqb25SmZDJc9aTUTPBrzq5IKrwM+N++WixMV/P3dG58/9GLN5Xf3hmP1WdwbKUQTI6sSFEt6yMUL191Z8p5j8UnX9/PkOGyricOMxz1k5/7fiw1HmPpufwOf8A/KzoHIZ9mPFMOm2+T3C5J3qiYe6j1ZlJw2NYITwTwIX9WJLMXBTHJTeh22GZy4BHo/eb/DQRkV7SvUlE+pXuTyKSmznnGi9V1MbM3gGc6pz7H/79u4HjnHPvj5Y5BzjHvz0KuKdrFeyMg4Ane12JNqj+vTWs9X+hc65RIsiuyXNv8tOH6f40rNfWIBn0fRjG+vfVvQlG8rPToF9XMPj7oPr3Vlv3pm6HZW4GDoneL/fTSpxzlwOXA5jZOufc6u5Vr3iDvg+qf2+p/l3T8N4Ew3V/Uv17b9D3QfXvmpH67DTo9YfB3wfVv7farX+3wzLvAFaZ2aFmNgs4C7ihy3UQEamke5OI9Cvdn0Qkt6723DnnJs3s/cDNwDhwpXPu3m7WQUSkku5NItKvdH8SkWZ0fRBz59yNwI05F7+8k3XpkkHfB9W/t1T/Lmny3gQDtG81qP69N+j7oPp3yYh9dhr0+sPg74Pq31tt1b+rCVVERERERESkM7r9zJ2IiIiIiIh0QN827szsVDP7uZltMLPze12fLGZ2pZk9YWb3RNMOMLO1ZvaA/73YTzcz+1u/Pz81s6N7V/NSXQ8xs++Z2X1mdq+ZfdBPH4h9MLM5ZvZjM/svX/+P++mHmtntvp5f9Q+gY2az/fsNfv7KXtY/MLNxM7vLzL7l3w9a/R82s7vN7Cdmts5PG4hrqBWDcG8C3Z96vQ+6P/VN/XV/6jO6N+neVATdm2rry8admY0Dfw+cBhwJvNPMjuxtrTJ9ETi1Ytr5wC3OuVXALf49JPuyyv+cA3yuS3WsZxL4iHPuSOB44Fx/nAdlH/YCb3bOvRJ4FXCqmR0PfBq41Dl3BLAdONsvfzaw3U+/1C/XDz4I3B+9H7T6A7zJOfeqKHXvoFxDTRmgexPo/tTrfdD9qX/o/tRfvojuTb2ke1P/6My9yTnXdz/Aa4Gbo/cXABf0ul416roSuCd6/3PgYP/6YODn/vXngXdmLdcvP8A3gZMGcR+AucCdwHEkAz/OqLyWSDKNvda/nuGXsx7Xe7n/A34z8C3ABqn+vi4PAwdVTBu4ayjnvg7MvcnXT/en/qi77k+92wfdn/qgbhl11b2pP+que1Pv9qFj96a+7LkDlgGPRu83+WmDYKlz7jH/eguw1L/u633y3dSvBm5ngPbBd8v/BHgCWAs8COxwzk36ReI6lurv5+8EDuxqhav9DXAeMO3fH8hg1R/AAd8xs/Vmdo6fNjDXUJMGvf4DeV50f+qZv0H3p775O8hhkOs/kOdE96ae+Rt0b6p5DXV9KIRR4pxzZtb36UjNbD7wr8AfOeeeNrPSvH7fB+fcFPAqM1sEfAN4SW9rlJ+ZvQ14wjm33sze2OPqtONXnHObzex5wFoz+1k8s9+voVE1KOdF96fe0P1JemVQzonuTb2he1Nj/dpztxk4JHq/3E8bBI+b2cEA/vcTfnpf7pOZzSS5OX3ZOfd1P3mg9gHAObcD+B5JV/wiMwtfXMR1LNXfz18IPNXdmpZ5PfBrZvYwcC1JeMFlDE79AXDObfa/nyD5J3EsA3gN5TTo9R+o86L7k+5P7dL9aWDqP1DnRPcm3Zva1cl7U7827u4AVvnMN7OAs4AbelynvG4A1vjXa0hiscP09/iMN8cDO6Ou156w5GumK4D7nXOfjWYNxD6Y2RL/rRNmth9JzPv9JDeqd/jFKusf9usdwHedD17uBefcBc655c65lSTX+Hedc7/NgNQfwMzmmdmC8Bo4GbiHAbmGWjDI9yYYoPOi+5PuT+3S/Wmg7k8Dc050b9K9qV0dvze1+0Bgp36AtwC/IIkD/tNe16dGHa8BHgP2kcS/nk0Sx3sL8ADw78ABflkjyWL1IHA3sLoP6v8rJDG/PwV+4n/eMij7ALwCuMvX/x7gz/30w4AfAxuA64HZfvoc/36Dn39Yr89BtC9vBL41aPX3df0v/3Nv+FsdlGuoxX3u+3uTr6fuT72tv+5Pva+37k99UKeMOure1Nv6697U+3p39N5kfiUREREREREZYP0alikiIiIiIiJNUONORERERERkCKhxJyIiIiIiMgTUuBMRERERERkCatyJiIiIiIgMATXuREREREREhoAadyIiIiIiIkNAjTsREREREZEh8P8DJmLGe4m2ZWgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1080x360 with 3 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from matplotlib.colors import LogNorm \n",
    "\n",
    "fig, ax = plt.subplots(1,3, figsize=(15,5))\n",
    "for i in range(3):\n",
    "    ax[i].set_title('Ph vs Time - channel {}'.format(2+i))\n",
    "    ax[i].hist2d(digitime[:,2+i], digiph[:,2+i], norm=LogNorm(), bins=100)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}