HHWWgg_hgg_coffea_Plotter.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.offsetbox import AnchoredText\n",
    "import pandas as pd "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {},
   "outputs": [],
   "source": [
    "X = pd.read_parquet(\"b95fc1b0-1c70-11eb-8b03-c22411acbeef_%2FEvents%3B1_0-3597.parquet\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['lead_charge', 'lead_cleanmask', 'lead_cutBased',\n",
       "       'lead_cutBased_Fall17V1Bitmap', 'lead_eCorr', 'lead_electronIdx',\n",
       "       'lead_electronIdxG', 'lead_electronVeto', 'lead_energyErr', 'lead_eta',\n",
       "       'lead_genPartFlav', 'lead_genPartIdx', 'lead_genPartIdxG', 'lead_hoe',\n",
       "       'lead_isScEtaEB', 'lead_isScEtaEE', 'lead_jetIdx', 'lead_jetIdxG',\n",
       "       'lead_mass', 'lead_mvaID', 'lead_mvaID_Fall17V1p1', 'lead_mvaID_WP80',\n",
       "       'lead_mvaID_WP90', 'lead_pdgId', 'lead_pfRelIso03_all',\n",
       "       'lead_pfRelIso03_chg', 'lead_phi', 'lead_pixelSeed', 'lead_pt',\n",
       "       'lead_r9', 'lead_seedGain', 'lead_sieie', 'lead_vidNestedWPBitmap',\n",
       "       'sublead_charge', 'sublead_cleanmask', 'sublead_cutBased',\n",
       "       'sublead_cutBased_Fall17V1Bitmap', 'sublead_eCorr',\n",
       "       'sublead_electronIdx', 'sublead_electronIdxG', 'sublead_electronVeto',\n",
       "       'sublead_energyErr', 'sublead_eta', 'sublead_genPartFlav',\n",
       "       'sublead_genPartIdx', 'sublead_genPartIdxG', 'sublead_hoe',\n",
       "       'sublead_isScEtaEB', 'sublead_isScEtaEE', 'sublead_jetIdx',\n",
       "       'sublead_jetIdxG', 'sublead_mass', 'sublead_mvaID',\n",
       "       'sublead_mvaID_Fall17V1p1', 'sublead_mvaID_WP80', 'sublead_mvaID_WP90',\n",
       "       'sublead_pdgId', 'sublead_pfRelIso03_all', 'sublead_pfRelIso03_chg',\n",
       "       'sublead_phi', 'sublead_pixelSeed', 'sublead_pt', 'sublead_r9',\n",
       "       'sublead_seedGain', 'sublead_sieie', 'sublead_vidNestedWPBitmap', 'pt',\n",
       "       'eta', 'phi', 'mass', 'charge', 'HHWWggTagger_40', 'best_tag', 'event',\n",
       "       'lumi', 'run'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 68,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "X.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {},
   "outputs": [],
   "source": [
    "categories = X['best_tag']%40\n",
    "cats = range(0, 4)\n",
    "objs = [\"Leptons\", \"Jets\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [],
   "source": [
    "for cat in cats:\n",
    "    for obj in objs:\n",
    "        exec(\"cat_%s_n%s = X['n%s'][categories == %s]\"%(cat, obj, obj, cat))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAS50lEQVR4nO3de4wdZ3nH8e8PO05KSEghgMAXnLSRhUHcuk3UgijlVqdAjCgqsYAKsGKQai5qCxhKRfgDEqqqoogU5CYhXIKNCVBsMAQETVPUUGynXBJMKpMm9eIQO6SBhFuU8PSPMyGbrdfe9ZndM7vz/Uir3fOemXeeXVm/HT/zzmyqCknSwvegURcgSZobBr4k9YSBL0k9YeBLUk8Y+JLUEwa+JPWEgS9JPWHga15KclOS50xz2ycn2ZPkZ83nJ8/l8Wcw5/lJPtbmnNJEBr4WtCRLgM8CHwN+E/gw8NlmXOoVA18jl2R5kk8nOZTkR0nen+S3kny1eX1bksuTnNJs/1FgBbAjyV1J3nyE6Z8JLAbeW1W/rKr3AQGedZSazk+yLclHktyZ5PokY9P4Xh6UZFOS7ze1b0vysOa9lUkqyYYkB5LckuQvm/fWAG8DXtp8T99qxh+TZHuS25PsS3LedGtM8pYkP2jeuyHJs49WvxY2A18jlWQR8DngZmAlsBTYyiCULwAeAzwOWA6cD1BVrwD+B3hhVT2kqv72CId4PPDteuAzRL7djB/NOU0tpwDbgfdPY5/XAy8C/qCp/X+BiyZt84fAGcDzgE1JnlNVXwTeDXyi+Z6e1Gy7BRhv5noJ8O5JwX3YGpOsAjYCv1tVJwF/BNw0jfq1gBn4GrUzGYTZm6rqp1X1i6r6WlXtq6ovN2flh4C/ZxCiM/UQ4MeTxn4MnDSNfb9WVTur6l7go8CTjrYD8Brgr6tqvKp+yeCX1EuSLJ6wzTub7/U7wIeAdYebKMly4OnAW5qfyzeBi4FXTKPGe4HjgdVJjquqm6rq+9OoXwuYga9RWw7cXFX3TBxM8sgkW5uWxE8Y9OBPPYb57wJOnjR2MnDnNPb94YSvfwacMCm4D+exwGeS3JHkDmAvg/B91IRt9k/4+mYGv/AO5zHA7VV156Ttlx6txqraB7yRwS+cg83PcqrjqCcMfI3afmDFYYL0AqCAJ1bVycDLGbR57jPdx7xeDzwxycR9n9iMz4b9wNlVdcqEjxOq6gcTtlk+4esVwIHm68nf0wHgYUlOmrT9D5iGqvp4VT2dwS+hAt4zk29EC4+Br1H7BnALcGGSE5OckORpDFoudwF3JFkKvGnSfrcCp09j/qsYnGG/PsnxSTY2419tofbjmnrv+1gMfBB4V5LHAiR5RJK1k/b7myQPTvJ44FXAJ5rxW4GVSR4EUFX7gX8HLmjmfyKwHrj8aIUlWZXkWUmOB34B/JzBz0E9ZuBrpJre8wuB32ZwIXYceCnwTuCpDPrtnwc+PWnXC4C3N62TvzrC/HczuIj6Z8AdwKuBFzXjw9rJIEjv+zgf+AcGF0+/lORO4OvAWZP2+1dgH/AV4O+q6kvN+Cebzz9Kcm3z9ToGF7MPAJ8B3lFVX55GbccDFwK3MWj7PJLBKiD1WPwDKNLcSLIS+G/guMnXLKS54Bm+JPWEga95L8nLmpuVJn8c8cJski9MsZ+tDy1ItnQkqSfm9Ay/WYWxJ8kL5vK4kqTBM0aOWZJLgRcAB6vqCRPG1zBYrbAIuLiqLmzeeguwbbrzn3rqqbVy5cphSpSkXtmzZ89tVfWIw703VOADlzF4dsdH7htono1yEfBcBkvsdiXZzuCuwe8CJ0x38pUrV7J79+4hS5Sk/khy81TvDRX4VXV1s9RsojOBfVV1Y3PwrcBaBs80ORFYDfw8yc6q+tUwx5ckTd+wZ/iHs5QHPitkHDirqjYCJHklcNtUYZ9kA7ABYMWKFbNQniT102wEfg4z9uulQFV12ZF2rqrNSW4BXrhkyZLfabk2Seqt2VilM84DHw61jPsfDjUtVbWjqjY89KEPbbUwSeqz2Qj8XcAZSU7L4M/Incvg2SKSpBEaKvCTbAGuAVYlGU+yvnlGyEbgSgbPAt9WVTN6FG2SFybZ/OMfT/67FZKkY9XpO23HxsbKZZmSNH1J9lTVYf/+ciefpeMZviS1r5OB70VbSWpfJwNfktS+Tga+LR1Jal8nA9+WjiS1r5OBL0lqXycD35aOJLWvk4FvS0eS2tfJwJcktc/Al6SeMPAlqSc6GfhetJWk9nUy8L1oK0nt62TgS5LaZ+BLUk8Y+JLUE50MfC/aSlL7Ohn4XrSVpPZ1MvAlSe0z8CWpJwx8SeoJA1+SesLAl6SeMPAlqSc6Gfiuw5ek9nUy8F2HL0nt62TgS5LaZ+BLUk8Y+JLUEwa+JPWEgS9JPWHgS1JPGPiS1BOL5+pASR4HvAE4FfhKVX1gro7ddSs3fX7UJczITRc+f9QlSDoGQ53hJ7k0ycEk100aX5PkhiT7kmwCqKq9VfVa4E+BsWGOK0mauWFbOpcBayYOJFkEXAScDawG1iVZ3bx3DvA14CtDHleSNENDBX5VXQ3cPmn4TGBfVd1YVXcDW4G1zfbbq+r3gZcNc1xJ0szNRg9/KbB/wutx4KwkzwReDBwP7Jxq5yQbgA0AK1asmIXyJKmfZiPwc5ixqqqrgKuOtnNVbQY2A4yNjVWrlUlSj83GssxxYPmE18uAAzOZwMcjS1L7ZiPwdwFnJDktyRLgXGD7TCbw8ciS1L5hl2VuAa4BViUZT7K+qu4BNgJXAnuBbVV1/Qzn9Qxfklo2VA+/qtZNMb6TI1yYnca8O4AdY2Nj5x3rHJKkB+rkoxU8w5ek9nUy8O3hS1L7Ohn4kqT2dTLwbelIUvs6Gfi2dCSpfZ0MfElS+wx8SeqJTga+PXxJal8nA98eviS1r5OBL0lqn4EvST3RycC3hy9J7etk4NvDl6T2dTLwJUntM/AlqScMfEnqiU4GvhdtJal9nQx8L9pKUvs6GfiSpPYZ+JLUEwa+JPWEgS9JPWHgS1JPGPiS1BOdDHzX4UtS+zoZ+K7Dl6T2dTLwJUntM/AlqScMfEnqCQNfknrCwJeknjDwJaknDHxJ6ok5C/wkL0ryT0k+m+R5c3VcSdLAUIGf5NIkB5NcN2l8TZIbkuxLsgmgqv65qs4DXgm8dJjjSpJmbtgz/MuANRMHkiwCLgLOBlYD65KsnrDJ25v3JUlzaKjAr6qrgdsnDZ8J7KuqG6vqbmArsDYD7wG+UFXXDnNcSdLMzUYPfymwf8Lr8WbsdcBzgJckee1UOyfZkGR3kt2HDh2ahfIkqZ8Wz8KcOcxYVdX7gPcdbeeq2gxsBhgbG6uWa5Ok3pqNM/xxYPmE18uAAzOZwMcjS1L7ZiPwdwFnJDktyRLgXGD7TCbw8ciS1L5hl2VuAa4BViUZT7K+qu4BNgJXAnuBbVV1/Qzn9Qxfklo2VA+/qtZNMb4T2DnEvDuAHWNjY+cd6xySpAfq5KMVPMOXpPZ1MvDt4UtS+zoZ+JKk9nUy8G3pSFL7Ohn4tnQkqX2dDHxJUvsMfEnqiU4Gvj18SWpfJwPfHr4kta+TgS9Jap+BL0k90cnAt4cvSe3rZODbw5ek9nUy8CVJ7TPwJaknDHxJ6olOBr4XbSWpfZ0MfC/aSlL7Ohn4kqT2GfiS1BMGviT1hIEvST1h4EtSTxj4ktQTnQx81+FLUvs6Gfiuw5ek9nUy8CVJ7TPwJaknDHxJ6gkDX5J6wsCXpJ4w8CWpJwx8SeqJOQv8JKcnuSTJFXN1TEnS/YYK/CSXJjmY5LpJ42uS3JBkX5JNAFV1Y1WtH+Z4kqRjN+wZ/mXAmokDSRYBFwFnA6uBdUlWD3kcSdKQhgr8qroauH3S8JnAvuaM/m5gK7B2mONIkoY3Gz38pcD+Ca/HgaVJHp7kg8BTkrx1qp2TbEiyO8nuQ4cOzUJ5ktRPi2dhzhxmrKrqR8Brj7ZzVW0GNgOMjY1Vy7VJUm/Nxhn+OLB8wutlwIGZTODjkSWpfbMR+LuAM5KclmQJcC6wfSYT+HhkSWrfsMsytwDXAKuSjCdZX1X3ABuBK4G9wLaqun6G83qGL0ktG6qHX1XrphjfCewcYt4dwI6xsbHzjnUOSdIDdfLRCp7hS1L7Ohn49vAlqX2dDHxJUvs6Gfi2dCSpfZ0MfFs6ktS+Tga+JKl9Br4k9UQnA98eviS1r5OBbw9fktrXycCXJLXPwJeknuhk4NvDl6T2dTLw7eFLUvs6GfiSpPYZ+JLUEwa+JPVEJwPfi7aS1L5OBr4XbSWpfZ0MfElS+wx8SeoJA1+SesLAl6SeMPAlqScMfEnqiU4GvuvwJal9nQx81+FLUvs6GfiSpPYZ+JLUEwa+JPWEgS9JPWHgS1JPGPiS1BMGviT1xOK5OlCSE4F/BO4Grqqqy+fq2JKkIc/wk1ya5GCS6yaNr0lyQ5J9STY1wy8Grqiq84BzhjmuJGnmhm3pXAasmTiQZBFwEXA2sBpYl2Q1sAzY32x275DHlSTN0FCBX1VXA7dPGj4T2FdVN1bV3cBWYC0wziD0hz6uJGnmZiN4l3L/mTwMgn4p8GngT5J8ANgx1c5JNiTZnWT3oUOHZqE8Seqn2bhom8OMVVX9FHjV0Xauqs3AZoCxsbFquTZJ6q3ZOMMfB5ZPeL0MODCTCXw8siS1bzYCfxdwRpLTkiwBzgW2z2QCH48sSe0bdlnmFuAaYFWS8STrq+oeYCNwJbAX2FZV189wXs/wJallQ/Xwq2rdFOM7gZ1DzLsD2DE2Nnbesc4hSXqgTi6P9AxfktrXycC3hy9J7etk4EuS2tfJwLelI0nt62Tg29KRpPZ1MvAlSe0z8CWpJzoZ+PbwJal9nQx8e/iS1L5OBr4kqX0GviT1RCcD3x6+JLWvk4FvD1+S2tfJwJcktc/Al6SeMPAlqSc6GfhetJWk9nUy8L1oK0nt62TgS5LaZ+BLUk8Y+JLUEwa+JPWEgS9JPWHgS1JPdDLwXYcvSe3rZOC7Dl+S2tfJwJcktc/Al6SeMPAlqScMfEnqCQNfknrCwJeknjDwJaknDHxJ6ok5C/wkpye5JMkVc3VMSdL9phX4SS5NcjDJdZPG1yS5Icm+JJuONEdV3VhV64cpVpJ07BZPc7vLgPcDH7lvIMki4CLgucA4sCvJdmARcMGk/V9dVQeHrlaSdMymFfhVdXWSlZOGzwT2VdWNAEm2Amur6gLgBcdaUJINwIbm5V1JbjjGqU4FbjvWOubYfKqVvGde1TufaoX5Ve98qhXmV73D1PrYqd6Y7hn+4SwF9k94PQ6cNdXGSR4OvAt4SpK3Nr8Y/p+q2gxsHqKu+463u6rGhp1nLsynWmF+1TufaoX5Ve98qhXmV72zVeswgZ/DjNVUG1fVj4DXDnE8SdIQhlmlMw4sn/B6GXBguHIkSbNlmMDfBZyR5LQkS4Bzge3tlNWKodtCc2g+1Qrzq975VCvMr3rnU60wv+qdlVpTNWUX5v6Nki3AMxlcSLgVeEdVXZLkj4H3MliZc2lVvWs2ipQkDW9agS9Jmv98tIIk9cSCC/yZ3P07alPdwdxFSZYn+Zcke5Ncn+QNo67pSJKckOQbSb7V1PvOUdd0NEkWJfnPJJ8bdS1Hk+SmJN9J8s0ku0ddz5EkOSXJFUm+1/z7/b1R1zSVJKuan+l9Hz9J8sbW5l9ILZ3m7t//YsLdv8C6qvruSAubQpJnAHcBH6mqJ4y6niNJ8mjg0VV1bZKTgD3Aizr8sw1wYlXdleQ44GvAG6rq6yMubUpJ/gIYA06uqmO+eXEuJLkJGKuqzt/IlOTDwL9V1cXNApMHV9UdIy7rqJo8+wFwVlXd3MacC+0M/9d3/1bV3cBWYO2Ia5pSVV0N3D7qOqajqm6pqmubr+8E9jK4+a6TauCu5uVxzUdnz26SLAOeD1w86loWkiQnA88ALgGoqrvnQ9g3ng18v62wh4UX+Ie7+7ezoTRfNY/ZeArwHyMu5YiaFsk3gYPAl6uqy/W+F3gz8KsR1zFdBXwpyZ7mcShddTpwCPhQ0y67OMmJoy5qms4FtrQ54UIL/Bnd/auZS/IQ4FPAG6vqJ6Ou50iq6t6qejKDmwLPTNLJtlmSFwAHq2rPqGuZgadV1VOBs4E/b9qTXbQYeCrwgap6CvBToNPX9gCa1tM5wCfbnHehBb53/86iphf+KeDyqvr0qOuZrua/8FcBa0ZbyZSeBpzT9MW3As9K8rHRlnRkVXWg+XwQ+AyDdmoXjQPjE/53dwWDXwBddzZwbVXd2uakCy3wu37377zVXAS9BNhbVX8/6nqOJskjkpzSfP0bwHOA7420qClU1VurallVrWTwb/arVfXyEZc1pSQnNhfuadojzwM6udKsqn4I7E+yqhl6NtDJhQaTrKPldg4M9/C0zqmqe5JsBK7k/rt/rx9xWVOaeAdzknGaO5hHW9WUnga8AvhO0xcHeFtV7RxdSUf0aODDzUqHBwHbqqrzyx3niUcBnxmcA7AY+HhVfXG0JR3R64DLm5PAG4FXjbieI0ryYAYrDV/T+twLaVmmJGlqC62lI0magoEvST1h4EtSTxj4ktQTBr4k9YSBL0k9YeBLUk/8H+joR+MrPoJlAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQQUlEQVR4nO3dcayddX3H8fdnVHSiwAKd01ItrIQMNybkBreROKfOlEmBbEukmW4qA0mG02zTVWMyXWJQ/zDGwJxVGKiMhiEqlTo0cYyZOaVFpzBk6QisV9RWGUV0s8K+++Mc4HLXc7m357n3PKe/9yu56T2/+5zf8z03zadPv+f3e06qCknSoe+nJl2AJGllGPiS1AgDX5IaYeBLUiMMfElqhIEvSY0w8CWpEQa+mpHkniQvW+SxL0iyM8mPhn++oIPz35zkDxd57KJrlRbLwJfmSXI48Gng48DPAFcBnx6OS1PLwNdUSrI2yfVJ9ib5fpJLk/x8ki8MH38vydVJjh4e/zHgucC2JA8lecsC078YWAW8v6p+XFUfAAK85ElqekeSa5N8NMkPktyRZGaB489K8rUkDyT55ySnjKo1ydOSfHz42h5IcmuSZy3hVyYZ+Jo+SQ4DPgPcC6wD1gBbGYTyJcBzgF8A1gLvAKiqVwP/CWysqmdU1XsXOMXzga/XE+878vXh+JM5e1jL0cANwKUjXsNpwBXA64FjgA8BNyR56oha/wA4aviajgEuAv57EfVIjzHwNY1OZxDqb66qH1bV/1TVF6tqV1V9fnhVvhd4H/DrBzH/M4B988b2Ac9cxHO/WFXbq+oR4GPAL4847gLgQ1X15ap6pKquAn4M/MqI43/CIOjXD4/fWVUPLqIe6TEGvqbRWuDeqnp47mCSn02yNcm3kjzIoAd/7EHM/xBw5LyxI4EfLOK535nz/Y+ApyVZdYDjngf86bA980CSBxi8rueMmPdjwE3A1iT3JXlvkqcsoh7pMQa+ptFu4LkHCNJLgAJOqaojgVcxaPM8arG3hr0DOCXJ3OeeMhzvym7gXVV19Jyvp1fVNQeqtap+UlXvrKqTgV8DzgJ+v8N61AADX9PoK8C3gXcnOWL4huYZDFouDwEPJFkDvHne874LnLCI+W8GHgH+OMlTk1w8HP9CJ9UPfBi4KMkLM3BEklckebRt9IRak/xGkl8avn/xIIMWzyMd1qMGGPiaOsP++EZgPYM3N2eBVwLvBE5j0G+/Ebh+3lMvAd4+bKH82QLz7wfOZXAF/QDwOuDc4fjY5Q/PsYNBH/9S4L+AXcBrFqj154DrGIT9ncA/MmhZSYsWPwBFWhlJbgP+sqo+Nela1Cav8KUVkOT5DJaKfnXStahdBr6alOT3hpua5n8t+MZsks+OeN7bFnjOe4DPAX9eVfd2/VqkxbKlI0mN8ApfkhpxoA0hvXHsscfWunXrJl2GJE2NnTt3fq+qVh/oZ70O/HXr1rFjx45JlyFJUyPJyPeJbOlIUiN6GfhJNibZsm/f/PtXSZIOVi8Dv6q2VdWFRx111KRLkaRDRi8DX5LUPQNfkhrR61U6UmvWbb5x0iUs2j3vfsWkS9ASeYUvSY3oZeC7SkeSutfLwHeVjiR1r5eBL0nqnoEvSY0w8CWpEQa+JDXCwJekRvQy8F2WKUnd62XguyxTkrrXy8CXJHXPwJekRnjzNB3SpulmZNNm2n633uzNK3xJaoaBL0mNMPAlqREGviQ1wsCXpEb0MvDdaStJ3etl4LvTVpK618vAlyR1z8CXpEYY+JLUCANfkhrhvXQkNWGa7v2zXPf98Qpfkhph4EtSIwx8SWqEgS9JjTDwJakRBr4kNcLAl6RGGPiS1IgVC/wk5yb5cJJPJ3n5Sp1XkjQwVuAnuSLJniS3zxvfkOSuJLuSbAaoqk9V1QXAa4BXjnNeSdLSjXuFfyWwYe5AksOAy4AzgZOBTUlOnnPI24c/lyStoLECv6puAe6fN3w6sKuq7q6q/cBW4JwMvAf4bFXdNmrOJBcm2ZFkx969e8cpT5I0x3L08NcAu+c8nh2OvQF4GfC7SS4a9eSq2lJVM1U1s3r16mUoT5LatBx3y8wBxqqqPgB8YBnOJ0lahOW4wp8F1s55fBxw31Im8EPMJal7yxH4twInJjk+yeHAecANS5nADzGXpO6NuyzzGuBLwElJZpOcX1UPAxcDNwF3AtdW1R1LnNcrfEnq2Fg9/KraNGJ8O7B9jHm3AdtmZmYuONg5JElP5K0VJKkRBr4kNaKXgW8PX5K618vAd5WOJHWvl4EvSeqegS9Jjehl4NvDl6Tu9TLw7eFLUveW4+ZpOsSt23zjpEuQdBB6eYUvSepeLwPfHr4kda+XgW8PX5K618vAlyR1z8CXpEYY+JLUCANfkhrRy8B3lY4kda+Xge8qHUnqXi8DX5LUPQNfkhph4EtSIwx8SWqEgS9Jjehl4LssU5K618vAd1mmJHWvl4EvSeqegS9JjTDwJakRBr4kNcLAl6RGGPiS1AgDX5IaYeBLUiN6GfjutJWk7vUy8N1pK0nd62XgS5K6Z+BLUiMMfElqhIEvSY0w8CWpEQa+JDXCwJekRhj4ktQIA1+SGmHgS1IjDHxJaoSBL0mNMPAlqRErFvhJTkhyeZLrVuqckqTHjRX4Sa5IsifJ7fPGNyS5K8muJJsBquruqjp/nPNJkg7euFf4VwIb5g4kOQy4DDgTOBnYlOTkMc8jSRrTWIFfVbcA988bPh3YNbyi3w9sBc5Z7JxJLkyyI8mOvXv3jlOeJGmO5ejhrwF2z3k8C6xJckySvwZOTfLWUU+uqi1VNVNVM6tXr16G8iSpTauWYc4cYKyq6vvARYuaINkIbFy/fn2nhUlSy5bjCn8WWDvn8XHAfUuZwM+0laTuLUfg3wqcmOT4JIcD5wE3LMN5JElLMO6yzGuALwEnJZlNcn5VPQxcDNwE3AlcW1V3LHHejUm27Nu3b5zyJElzjNXDr6pNI8a3A9vHmHcbsG1mZuaCg51DkvRE3lpBkhph4EtSI3oZ+PbwJal7vQx8l2VKUvd6GfiSpO4Z+JLUiF4Gvj18SepeLwPfHr4kda+XgS9J6p6BL0mN6GXg28OXpO71MvDt4UtS93oZ+JKk7hn4ktQIA1+SGmHgS1Ijehn4rtKRpO71MvBdpSNJ3etl4EuSumfgS1IjDHxJaoSBL0mNMPAlqRG9DHyXZUpS93oZ+C7LlKTu9TLwJUndM/AlqREGviQ1wsCXpEYY+JLUCANfkhph4EtSIwx8SWrEqkkXcCBJNgIb169fP+lSVsS6zTdOugRJDejlFb47bSWpe70MfElS9wx8SWqEgS9JjTDwJakRBr4kNcLAl6RGGPiS1AgDX5IaYeBLUiMMfElqhIEvSY0w8CWpEQa+JDVixW6PnOQI4K+A/cDNVXX1Sp1bkjTmFX6SK5LsSXL7vPENSe5KsivJ5uHwbwPXVdUFwNnjnFeStHTjtnSuBDbMHUhyGHAZcCZwMrApycnAccDu4WGPjHleSdISjRX4VXULcP+84dOBXVV1d1XtB7YC5wCzDEJ/wfMmuTDJjiQ79u7dO055kqQ5luNN2zU8fiUPg6BfA1wP/E6SDwLbRj25qrZU1UxVzaxevXoZypOkNi3Hm7Y5wFhV1Q+B1y5qgsY+01aSVsJyXOHPAmvnPD4OuG8pE/iZtpLUveUI/FuBE5Mcn+Rw4DzghmU4jyRpCcZdlnkN8CXgpCSzSc6vqoeBi4GbgDuBa6vqjiXOuzHJln379o1TniRpjrF6+FW1acT4dmD7GPNuA7bNzMxccLBzSJKeyFsrSFIjDHxJakQvA98eviR1r5eB77JMSepeLwNfktQ9A1+SGtHLwLeHL0nd62Xg28OXpO71MvAlSd0z8CWpEb0MfHv4ktS9Xga+PXxJ6l4vA1+S1D0DX5IaYeBLUiMMfElqRC8D31U6ktS9Xga+q3QkqXu9DHxJUvcMfElqhIEvSY0w8CWpEQa+JDWil4HvskxJ6l4vA99lmZLUvV4GviSpewa+JDXCwJekRhj4ktQIA1+SGmHgS1IjDHxJaoSBL0mN6GXgu9NWkrrXy8B3p60kda+XgS9J6p6BL0mNMPAlqREGviQ1wsCXpEYY+JLUCANfkhph4EtSIwx8SWqEgS9JjTDwJakRBr4kNcLAl6RGrFjgJzkhyeVJrlupc0qSHreowE9yRZI9SW6fN74hyV1JdiXZvNAcVXV3VZ0/TrGSpIO3apHHXQlcCnz00YEkhwGXAb8JzAK3JrkBOAy4ZN7zX1dVe8auVpJ00BYV+FV1S5J184ZPB3ZV1d0ASbYC51TVJcBZB1tQkguBC4cPH0py10FOdSzwvYOtY4VNU60wXfVOU60wXfVOU60wRfXmPWPV+rxRP1jsFf6BrAF2z3k8C7xw1MFJjgHeBZya5K3Dfxj+n6raAmwZo65Hz7ejqmbGnWclTFOtMF31TlOtMF31TlOtMF31Llet4wR+DjBWow6uqu8DF41xPknSGMZZpTMLrJ3z+DjgvvHKkSQtl3EC/1bgxCTHJzkcOA+4oZuyOjF2W2gFTVOtMF31TlOtMF31TlOtMF31LkutqRrZhXn8oOQa4MUM3vT4LvAXVXV5kt8C3s9gZc4VVfWu5ShSkjS+RQW+JGn6eWsFSWrEIRf4S9n9O2mjdjD3UZK1Sf4hyZ1J7kjyxknXtJAkT0vylST/Oqz3nZOu6ckkOSzJV5N8ZtK1PJkk9yT5RpKvJdkx6XoWkuToJNcl+ebw7++vTrqmUZKcNPydPvr1YJI3dTb/odTSGe7+/Xfm7P4FNlXVv020sBGSvAh4CPhoVf3ipOtZSJJnA8+uqtuSPBPYCZzb499tgCOq6qEkTwG+CLyxqv5lwqWNlORPgBngyKo66M2LKyHJPcBMVfV+I1OSq4B/qqqPDBeYPL2qHphwWU9qmGffAl5YVfd2MeehdoX/2O7fqtoPbAXOmXBNI1XVLcD9k65jMarq21V12/D7HwB3Mth810s18NDw4VOGX729uklyHPAK4COTruVQkuRI4EXA5QBVtX8awn7opcB/dBX2cOgF/oF2//Y2lKbV8DYbpwJfnnApCxq2SL4G7AE+X1V9rvf9wFuA/51wHYtVwOeS7BzeDqWvTgD2An8zbJd9JMkRky5qkc4DrulywkMt8Je0+1dLl+QZwCeAN1XVg5OuZyFV9UhVvYDBpsDTk/SybZbkLGBPVe2cdC1LcEZVnQacCfzRsD3ZR6uA04APVtWpwA+BXr+3BzBsPZ0N/F2X8x5qge/u32U07IV/Ari6qq6fdD2LNfwv/M3AhslWMtIZwNnDvvhW4CVJPj7ZkhZWVfcN/9wDfJJBO7WPZoHZOf+7u47BPwB9dyZwW1V9t8tJD7XA7/vu36k1fBP0cuDOqnrfpOt5MklWJzl6+P1PAy8DvjnRokaoqrdW1XFVtY7B39kvVNWrJlzWSEmOGL5xz7A98nKglyvNquo7wO4kJw2HXgr0cqHBPJvouJ0D4908rXeq6uEkFwM38fju3zsmXNZIc3cwJ5lluIN5slWNdAbwauAbw744wNuqavvkSlrQs4Grhisdfgq4tqp6v9xxSjwL+OTgGoBVwN9W1d9PtqQFvQG4engReDfw2gnXs6AkT2ew0vD1nc99KC3LlCSNdqi1dCRJIxj4ktQIA1+SGmHgS1IjDHxJaoSBL0mNMPAlqRH/B1rP8sXQYJRsAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAARlUlEQVR4nO3df4xlZX3H8ffHXZCKP4iyGl1WFwohbo0WM4VUjbVq7VIFjDWRjZpoCUharKatitZG/UPBpjG2kWpWQPyBrIjQLrqKRmspKVYWqlGKNCuF7rjqrlIVbJWA3/5xz5ZxOrM7s/fM3DPzvF/JZO997jnP+d7J5nPPPOd5zk1VIUla/R4y6QIkScvDwJekRhj4ktQIA1+SGmHgS1IjDHxJaoSBL0mNMPC1IiW5M8nzF7jt1iS3J/lFklct9/EX0efbk3yszz6lmQx8teDrwB8Ct0y6EGmSDHxNXJINSa5Osi/JD5O8L8mvJvlS9/wHSS5PclS3/UeBJwLXJrk3yRsP1H9VXVRVXwR+toia3p7kyiQfSXJPkluTTC1gv4ckOT/Jt7var0zy6O61jUkqyTlJ9iT5bpI/7V7bDLwFeFn3nr7etT8hyfYkdyfZleTshdaY5E1JvtO9dnuS5y30/Wt1MvA1UUnWAJ8G7gI2AuuBbUCAC4AnAE8GNgBvB6iqVwL/CZxWVQ+vqr9covJO72o5CtgOvG8B+/wx8GLgtxjV/l/ARbO2+W3gBOAFwPlJnl9VnwPeBXyie09P67a9Apju+nop8K5ZwT1njUlOBM4DfqOqHgH8LnDnwt62VisDX5N2MqMwe0NV/bSqflZVN1TVrqr6QlX9vKr2Ae9hFKLL6Yaq2lFVDwAfBZ52sB2A1wB/XlXTVfVzRh9SL02ydsY27+je6zeADwFb5uooyQbgWcCbut/L14CLgVcuoMYHgIcCm5IcVlV3VtW3F/i+tUoZ+Jq0DcBdVXX/zMYkj02yrRuS+AnwMeDoZa7tezMe/zdwxKzgnsuTgGuS/CjJj4DbGIXv42Zss3vG47sYfeDN5QnA3VV1z6zt1x+sxqraBbye0QfO3u53Od9x1AgDX5O2G3jiHEF6AVDAU6vqkcArGA3z7DfU27zuBk6tqqNm/BxRVd+Zsc2GGY+fCOzpHs9+T3uARyd5xKztv8MCVNXHq+pZjD6ECnj3Yt6IVh8DX5P2VeC7wIVJjkxyRJJnAo8A7gV+lGQ98IZZ+30fOG4hB0hyeJIjGH1gHNYdo4//+/v72v+zFvgA8M4kT+qOvS7JGbP2+4skD0vya8CrgU/MeE8b99dWVbuBfwYu6Pp/KnAWcPkC3vOJSZ6b5KGMLlb/D6O/NNQwA18T1Y09nwYcz+hC7DTwMuAdwNOBHwOfAa6etesFwFu7oZM/O8hhPs8o8J4BbO0eP7uH8nd0fe3/eTvw14wunn4+yT3AV4BTZu33j8Au4IvAX1XV57v2T3b//jDJ/imkWxhdzN4DXAO8raq+sIDaHgpcCPyA0bDPYxnNAlLD4hegSMsjyUbgP4DDZl+zkJaDZ/iS1AgDXytekpd3i5Vm/9x6kP0+O89+Dn1oVXJIR5Ia4Rm+JDXiYItIepPkycDrGC2e+WJVvf9g+xx99NG1cePGpS5NklaNm2+++QdVtW6u18YK/CSXAi8C9lbVU2a0b2Y0PW0NcHFVXVhVtwHndnOMP7iQ/jdu3MjOnTvHKVGSmpLkrvleG3dI5zJg86yDrWF0s6hTgU3AliSbutdOB25gNP9YkrSMxgr8qroeuHtW88nArqq6o6ruY3QnvzO67bdX1TOAl8/XZ3fr2J1Jdu7bt2+c8iRJMyzFGP56fvnmUNPAKUmeA7yE0QrAHfPtXFVbGa2GZGpqyilEktSTpQj8zNFWVfVl4MtLcDxJ0gIsxbTMaX75boDH8ODdABckyWlJtv74xz/utTBJatlSBP5NwAlJjk1yOHAmo5tJLVhVXVtV5zzqUY9agvIkqU1jBX6SK4AbgROTTCc5q7sp1HnAdYy+/OHKqjrgEvc5+vUMX5J6NuhbK0xNTZXz8CVp4ZLcXFVTc73mrRUkqRGDDHyHdCSpf4MMfC/aSlL/Bhn4kqT+GfiS1IhBBr5j+JLUv0EGvmP4ktS/QQa+JKl/Br4kNWKQge8YviT1b5CB7xi+JPVvkIEvSeqfgS9JjTDwJakRgwx8L9pKUv8GGfhetJWk/g0y8CVJ/TPwJakRBr4kNcLAl6RGGPiS1IhBBr7TMiWpf2snXcBcqupa4NqpqamzD7WPjed/pseKltadF75w0iVIasAgz/AlSf0z8CWpEQa+JDXCwJekRhj4ktQIA1+SGmHgS1IjDHxJasQgA9+VtpLUv0EGvl+AIkn9G2TgS5L6Z+BLUiMMfElqhIEvSY0w8CWpEQa+JDXCwJekRhj4ktQIA1+SGmHgS1Ijli3wk7w4yQeT/H2SFyzXcSVJI2MFfpJLk+xN8s1Z7ZuT3J5kV5LzAarq76rqbOBVwMvGOa4kafHGPcO/DNg8syHJGuAi4FRgE7AlyaYZm7y1e12StIzGCvyquh64e1bzycCuqrqjqu4DtgFnZOTdwGer6pZxjitJWrylGMNfD+ye8Xy6a3st8HzgpUnOnW/nJOck2Zlk5759+5agPElq09ol6DNztFVV/Q3wNwfbuaq2AlsBpqamqufaJKlZS3GGPw1smPH8GGDPYjrwG68kqX9LEfg3ASckOTbJ4cCZwPbFdOA3XklS/8adlnkFcCNwYpLpJGdV1f3AecB1wG3AlVV16/ilSpLGMdYYflVtmad9B7DjUPtNchpw2vHHH3+oXUiSZhnkrRUc0pGk/g0y8CVJ/Rtk4DtLR5L6N8jAd0hHkvo3yMCXJPXPwJekRgwy8B3Dl6T+DTLwHcOXpP4NMvAlSf0z8CWpEYMMfMfwJal/gwx8x/AlqX+DDHxJUv8MfElqhIEvSY0YZOB70VaS+jfIwPeirST1b5CBL0nqn4EvSY0w8CWpEQa+JDXCwJekRgwy8J2WKUn9G2TgOy1Tkvo3yMCXJPXPwJekRhj4ktQIA1+SGmHgS1IjDHxJaoSBL0mNGGTgu/BKkvo3yMB34ZUk9W+QgS9J6p+BL0mNMPAlqREGviQ1wsCXpEYY+JLUCANfkhph4EtSIwx8SWqEgS9JjTDwJakRyxb4SY5LckmSq5brmJKkB40V+EkuTbI3yTdntW9OcnuSXUnOB6iqO6rqrHGOJ0k6dOOe4V8GbJ7ZkGQNcBFwKrAJ2JJk05jHkSSNaazAr6rrgbtnNZ8M7OrO6O8DtgFnLLTPJOck2Zlk5759+8YpT5I0w1KM4a8Hds94Pg2sT/KYJB8ATkry5vl2rqqtVTVVVVPr1q1bgvIkqU1rl6DPzNFWVfVD4NwlOJ4kaQGW4gx/Gtgw4/kxwJ7FdOBXHEpS/5Yi8G8CTkhybJLDgTOB7YvpwK84lKT+jTst8wrgRuDEJNNJzqqq+4HzgOuA24Arq+rW8UuVJI1jrDH8qtoyT/sOYMeh9pvkNOC0448//lC7kCTNMshbKzikI0n9G2TgS5L6N8jAd5aOJPVvkIHvkI4k9W+QgS9J6t8gA98hHUnq3yAD3yEdSerfIANfktQ/A1+SGmHgS1IjBhn4XrSVpP4NMvC9aCtJ/Rtk4EuS+mfgS1IjDHxJasQgA9+LtpLUv0EGvhdtJal/gwx8SVL/DHxJaoSBL0mNMPAlqRGDDHxn6UhS/wYZ+M7SkaT+DTLwJUn9M/AlqREGviQ1wsCXpEYY+JLUCANfkhph4EtSIwx8SWrEIAPflbaS1L9BBr4rbSWpf4MMfElS/wx8SWqEgS9JjTDwJakRBr4kNcLAl6RGGPiS1AgDX5IaYeBLUiMMfElqhIEvSY1Yu1wHSnIk8LfAfcCXq+ry5Tq2JGnMM/wklybZm+Sbs9o3J7k9ya4k53fNLwGuqqqzgdPHOa4kafHGHdK5DNg8syHJGuAi4FRgE7AlySbgGGB3t9kDYx5XkrRIYwV+VV0P3D2r+WRgV1XdUVX3AduAM4BpRqF/wOMmOSfJziQ79+3bN055kqQZluKi7XoePJOHUdCvB64Gfj/J+4Fr59u5qrZW1VRVTa1bt24JypOkNi3FRdvM0VZV9VPg1UtwPEnSAizFGf40sGHG82OAPYvpwK84lKT+LUXg3wSckOTYJIcDZwLbF9OBX3EoSf0bd1rmFcCNwIlJppOcVVX3A+cB1wG3AVdW1a2L7NczfEnq2Vhj+FW1ZZ72HcCOMfq9Frh2amrq7EPtQ5L0y7y1giQ1YpCB75COJPVvkIHvRVtJ6t8gA1+S1D8DX5IaMcjAdwxfkvo3yMB3DF+S+jfIwJck9c/Al6RGDDLwHcOXpP4NMvAdw5ek/g0y8CVJ/TPwJakRBr4kNWKQge9FW0nq3yAD34u2ktS/QQa+JKl/Br4kNcLAl6RGGPiS1AgDX5IaMcjAd1qmJPVvkIHvtExJ6t8gA1+S1D8DX5IaYeBLUiMMfElqhIEvSY0w8CWpEQa+JDXCwJekRgwy8F1pK0n9G2Tgu9JWkvo3yMCXJPXPwJekRhj4ktQIA1+SGmHgS1IjDHxJaoSBL0mNMPAlqREGviQ1wsCXpEYY+JLUiGUL/CTHJbkkyVXLdUxJ0oMWFPhJLk2yN8k3Z7VvTnJ7kl1Jzj9QH1V1R1WdNU6xkqRDt3aB210GvA/4yP6GJGuAi4DfAaaBm5JsB9YAF8za/w+qau/Y1UqSDtmCAr+qrk+ycVbzycCuqroDIMk24IyqugB40aEWlOQc4Jzu6b1Jbj/Ero4GfnCodSynvHvl1NpZSfWupFphZdW7kmqFlVXvOLU+ab4XFnqGP5f1wO4Zz6eBU+bbOMljgHcCJyV5c/fB8P9U1VZg6xh17T/ezqqaGref5bCSaoWVVe9KqhVWVr0rqVZYWfUuVa3jBH7maKv5Nq6qHwLnjnE8SdIYxpmlMw1smPH8GGDPeOVIkpbKOIF/E3BCkmOTHA6cCWzvp6xejD0stIxWUq2wsupdSbXCyqp3JdUKK6veJak1VfOOwjy4UXIF8BxGFxK+D7ytqi5J8nvAexnNzLm0qt65FEVKksa3oMCXJK183lpBkhqx6gJ/Mat/J22+FcxDlGRDkn9IcluSW5O8btI1HUiSI5J8NcnXu3rfMemaDibJmiT/muTTk67lYJLcmeQbSb6WZOek6zmQJEcluSrJt7r/v7856Zrmk+TE7ne6/+cnSV7fW/+raUinW/3778xY/Qtsqap/m2hh80jybOBe4CNV9ZRJ13MgSR4PPL6qbknyCOBm4MUD/t0GOLKq7k1yGHAD8Lqq+sqES5tXkj8BpoBHVtUhL15cDknuBKaqavALmZJ8GPinqrq4m2DysKr60YTLOqguz74DnFJVd/XR52o7w/+/1b9VdR+wDThjwjXNq6quB+6edB0LUVXfrapbusf3ALcxWnw3SDVyb/f0sO5nsGc3SY4BXghcPOlaVpMkjwSeDVwCUFX3rYSw7zwP+HZfYQ+rL/DnWv072FBaqbrbbJwE/MuESzmgbojka8Be4AtVNeR63wu8EfjFhOtYqAI+n+Tm7nYoQ3UcsA/4UDdcdnGSIydd1AKdCVzRZ4erLfAXtfpXi5fk4cCngNdX1U8mXc+BVNUDVfXrjBYFnpxkkMNmSV4E7K2qmyddyyI8s6qeDpwK/FE3PDlEa4GnA++vqpOAnwKDvrYH0A09nQ58ss9+V1vgu/p3CXVj4Z8CLq+qqyddz0J1f8J/Gdg82Urm9Uzg9G5cfBvw3CQfm2xJB1ZVe7p/9wLXMBpOHaJpYHrGX3dXMfoAGLpTgVuq6vt9drraAn/oq39XrO4i6CXAbVX1nknXczBJ1iU5qnv8K8DzgW9NtKh5VNWbq+qYqtrI6P/sl6rqFRMua15Jjuwu3NMNj7wAGORMs6r6HrA7yYld0/OAQU40mGULPQ/nwHg3Txucqro/yXnAdTy4+vfWCZc1r5krmJNM061gnmxV83om8ErgG924OMBbqmrH5Eo6oMcDH+5mOjwEuLKqBj/dcYV4HHDN6ByAtcDHq+pzky3pgF4LXN6dBN4BvHrC9RxQkocxmmn4mt77Xk3TMiVJ81ttQzqSpHkY+JLUCANfkhph4EtSIwx8SWqEgS9JjTDwJakR/wulfpiSPtFaTQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAP2klEQVR4nO3df6zddX3H8efLFsZELQl0TttqYRCybnODNLhJ4n7oTBkUzLZEmummY1SyYTTbdGhMpksM6B/GGJizCuIPRsMQXSt1aOYcM3NKixrpKktHYL2itsoo4qYIvvfHPYy7u57LvT3fc8/39PN8JA093577/b4vgWe//Xy/39NUFZKkY99TJj2AJGl5GHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHw1I8m9SV68yPduS3J3kh8leWVHx/9skj9Y5HsXPau0WAZfOrKvAH8I3DnpQaSuGHxNpSTrktyS5FCS7yS5OslPJfnM4PW3k9yQ5KTB+z8MPAfYmeThJG9YaP9VdU1V/T3w/SXM9JYkNyX5UJLvJtmbZOMC778gyZeTPJjkn5M8b9isSU5I8pHB9/ZgkjuSPHOxs0lg8DWFkqwAPgHcB6wH1gDbgQBXAs8GfhpYB7wFoKpeAfwHsLmqnlZV7xjTeBcOZjkJ2AFcPeR7OBu4Dng1cDLwXmBHkh8bMuvvAasG39PJwGXAf4/pe9AxyuBrGp3DbNRfX1Xfq6rvV9Xnqmp/VX26qn5QVYeAdwK/vMyzfa6qdlXVY8CHgZ8f8r5LgfdW1Req6rGq+iDwA+AXh7z/h8yG/vTB+/dU1UOdT69jmsHXNFoH3FdVj87dmOQnkmxP8vUkDwEfAU5Z5tm+Oefn/wWckGTlEd73XOBPBsszDyZ5kNnv69lD9vth4DZge5L7k7wjyXFdDq5jn8HXNDoAPOcIIb0SKOB5VfUM4OXMLvM8rk8fDXsAeFtVnTTnx1Or6sbBr/+fWavqh1X11qraALwAuAD43WWeWVPO4GsafRH4BnBVkhMHFzTPBZ4OPAw8mGQN8Pp5X/ct4LTFHCDJ8UlOYPY3jOMGx+jy/5f3AZcleX5mnZjk/CRPP9KsSX41yc8Nrl88xOwSz2MdzqMGGHxNncH6+GbgdGYvbs4ALwPeCpwNHAZuBW6Z96VXAm8eLKH86ZMc5lPMXhR9AbBt8PMXdjH+4HvYzew6/tXAfwL7gVcuMOtPAjczG/t9wD8yu2QlLVr8C1Ck5ZHkTuAvqurjk55FbfIMX1oGSX6G2VtFvzTpWdQug68mJfmdwUNN83/sfZKv++SQr3vTAl/zdmaXiP6squ7r+nuRFsslHUlqhGf4ktSIIz0Q0hunnHJKrV+/ftJjSNLU2LNnz7eravWRfq3XwV+/fj27d++e9BiSNDWSDL1O5JKOJDWi12f4UmvWX3HrpEdYtHuvOn/SI2iJPMOXpEYYfElqhMGXpEYYfElqhMGXpEYYfElqhMGXpEYYfElqhMGXpEYYfElqhMGXpEYYfElqxLIFP8lLk7wvyd8meclyHVeSNGuk4Ce5LsnBJHfN274pyd1J9ie5AqCqPl5VlwKvBF42ynElSUs36hn+9cCmuRuSrACuAc4DNgBbkmyY85Y3D35dkrSMRgp+Vd0OPDBv8znA/qq6p6oeAbYDF2XW24FPVtWdw/aZZGuS3Ul2Hzp0aJTxJElzjGMNfw1wYM7rmcG21wAvBn47yWXDvriqtlXVxqrauHr1Ef9aRknSURjH33iVI2yrqno38O4xHE+StAjjOMOfAdbNeb0WuH8pO0iyOcm2w4cPdzqYJLVsHMG/AzgjyalJjgcuBnYsZQdVtbOqtq5atWoM40lSm0a9LfNG4PPAmUlmklxSVY8ClwO3AfuAm6pq7+ijSpJGMdIaflVtGbJ9F7DraPebZDOw+fTTTz/aXUiS5hnHRduRVdVOYOfGjRsvnfQsko5s/RW3TnqEJbn3qvMnPcLE+Vk6ktQIgy9Jjehl8L0tU5K618vge1umJHWvl8GXJHXP4EtSIwy+JDWil8H3oq0kda+XwfeirSR1r5fBlyR1z+BLUiMMviQ1opfB96KtJHWvl8H3oq0kda+XwZckdc/gS1IjDL4kNcLgS1Ijehl879KRpO71MvjepSNJ3etl8CVJ3TP4ktSIlZMeQJKWw/orbp30CIt271Xnj2W/nuFLUiMMviQ1wuBLUiN6GXzvw5ek7vUy+N6HL0nd62XwJUndM/iS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1IheBt/P0pGk7vUy+H6WjiR1r5fBlyR1z+BLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1YtmCn+S0JNcmuXm5jilJesJIwU9yXZKDSe6at31TkruT7E9yBUBV3VNVl4xyPEnS0Rv1DP96YNPcDUlWANcA5wEbgC1JNox4HEnSiEYKflXdDjwwb/M5wP7BGf0jwHbgosXuM8nWJLuT7D506NAo40mS5hjHGv4a4MCc1zPAmiQnJ/kr4Kwkbxz2xVW1rao2VtXG1atXj2E8SWrTyjHsM0fYVlX1HeCyMRxPkrQI4zjDnwHWzXm9Frh/KTtIsjnJtsOHD3c6mCS1bBzBvwM4I8mpSY4HLgZ2LGUHVbWzqrauWrVqDONJUptGvS3zRuDzwJlJZpJcUlWPApcDtwH7gJuqau/oo0qSRjHSGn5VbRmyfRewa5R9S5K61cuPVnANX5K618vgu4YvSd3rZfAlSd3rZfBd0pGk7vUy+C7pSFL3ehl8SVL3DL4kNcLgS1Ijehl8L9pKUvd6GXwv2kpS93oZfElS9wy+JDXC4EtSI3oZfC/aSlL3ehl8L9pKUvd6GXxJUvcMviQ1wuBLUiMMviQ1wuBLUiN6GXxvy5Sk7vUy+N6WKUnd62XwJUndM/iS1AiDL0mNMPiS1AiDL0mNMPiS1IheBt/78CWpe70MvvfhS1L3ehl8SVL3DL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1Ijehl8P1pBkrrXy+D70QqS1L1eBl+S1D2DL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1IiVy3WgJCcCfwk8Any2qm5YrmNLkkY8w09yXZKDSe6at31TkruT7E9yxWDzbwI3V9WlwIWjHFeStHSjLulcD2yauyHJCuAa4DxgA7AlyQZgLXBg8LbHRjyuJGmJRgp+Vd0OPDBv8znA/qq6p6oeAbYDFwEzzEZ/weMm2Zpkd5Ldhw4dGmU8SdIc47hou4YnzuRhNvRrgFuA30ryHmDnsC+uqm1VtbGqNq5evXoM40lSm8Zx0TZH2FZV9T3gVWM4niRpEcZxhj8DrJvzei1w/xiOI0lagnEE/w7gjCSnJjkeuBjYsZQdJNmcZNvhw4fHMJ4ktWnU2zJvBD4PnJlkJsklVfUocDlwG7APuKmq9i5lv1W1s6q2rlq1apTxJElzjLSGX1VbhmzfBewaZd+SpG718qMVXNKRpO71Mvgu6UhS93oZfElS9wy+JDWil8F3DV+SutfL4LuGL0nd62XwJUndM/iS1IheBt81fEnqXi+D7xq+JHWvl8GXJHXP4EtSIwy+JDWil8H3oq0kda+XwfeirSR1r5fBlyR1z+BLUiMMviQ1wuBLUiMMviQ1opfB97ZMSepeL4PvbZmS1L1eBl+S1D2DL0mNMPiS1AiDL0mNMPiS1AiDL0mN6GXwvQ9fkrrXy+B7H74kda+XwZckdc/gS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1Ijehl8P0tHkrrXy+D7WTqS1L1eBl+S1D2DL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1IhlC36S05Jcm+Tm5TqmJOkJiwp+kuuSHExy17ztm5LcnWR/kisW2kdV3VNVl4wyrCTp6K1c5PuuB64GPvT4hiQrgGuAXwdmgDuS7ABWAFfO+/rfr6qDI08rSTpqiwp+Vd2eZP28zecA+6vqHoAk24GLqupK4IKjHSjJVmDr4OXDSe4+yl2dAnz7aOdYZtM0K0zXvNM0K0zXvNM0K0zRvHn7SLM+d9gvLPYM/0jWAAfmvJ4Bnj/szUlOBt4GnJXkjYPfGP6fqtoGbBthrsePt7uqNo66n+UwTbPCdM07TbPCdM07TbPCdM07rllHCX6OsK2GvbmqvgNcNsLxJEkjGOUunRlg3ZzXa4H7RxtHkjQuowT/DuCMJKcmOR64GNjRzVidGHlZaBlN06wwXfNO06wwXfNO06wwXfOOZdZUDV2FeeJNyY3ArzB70eNbwJ9X1bVJfgN4F7N35lxXVW8bx5CSpNEtKviSpOnnRytIUiOOueAv5enfSRv2BHMfJVmX5B+S7EuyN8lrJz3TQpKckOSLSb4ymPetk57pySRZkeRLST4x6VmeTJJ7k3w1yZeT7J70PAtJclKSm5N8bfDf7y9NeqZhkpw5+Hf6+I+Hkryus/0fS0s6g6d//405T/8CW6rqXyc62BBJXgg8DHyoqn520vMsJMmzgGdV1Z1Jng7sAV7a43+3AU6sqoeTHAd8DnhtVf3LhEcbKskfAxuBZ1TVUT+8uByS3AtsrKreP8iU5IPAP1XV+wc3mDy1qh6c8FhPatCzrwPPr6r7utjnsXaG/79P/1bVI8B24KIJzzRUVd0OPDDpORajqr5RVXcOfv5dYB+zD9/1Us16ePDyuMGP3p7dJFkLnA+8f9KzHEuSPAN4IXAtQFU9Mg2xH3gR8O9dxR6OveAf6enf3kZpWg0+ZuMs4AsTHmVBgyWSLwMHgU9XVZ/nfRfwBuBHE55jsQr4VJI9g49D6avTgEPABwbLZe9PcuKkh1qki4Ebu9zhsRb8JT39q6VL8jTgo8DrquqhSc+zkKp6rKp+gdmHAs9J0stlsyQXAAeras+kZ1mCc6vqbOA84I8Gy5N9tBI4G3hPVZ0FfA/o9bU9gMHS04XA33S532Mt+D79O0aDtfCPAjdU1S2TnmexBn+E/yywabKTDHUucOFgXXw78GtJPjLZkRZWVfcP/nkQ+Bizy6l9NAPMzPnT3c3M/gbQd+cBd1bVt7rc6bEW/L4//Tu1BhdBrwX2VdU7Jz3Pk0myOslJg5//OPBi4GsTHWqIqnpjVa2tqvXM/jf7map6+YTHGirJiYML9wyWR14C9PJOs6r6JnAgyZmDTS8CenmjwTxb6Hg5B0b78LTeqapHk1wO3MYTT//unfBYQ819gjnJDIMnmCc71VDnAq8AvjpYFwd4U1XtmtxIC3oW8MHBnQ5PAW6qqt7f7jglngl8bPYcgJXAX1fV3012pAW9BrhhcBJ4D/CqCc+zoCRPZfZOw1d3vu9j6bZMSdJwx9qSjiRpCIMvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUiP8BwT61C1loe9YAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAARFUlEQVR4nO3db4xldX3H8fcHFrSigMpqdFkc6NJNqdFiVkyDsa1SuxQXbGsim+oDSlxNitXYqqu1UR8o2BirRqqhgPgHQcR/i6yCqbXUVgsLxShFmpUu3WGVXaAq2CqFfvtgzpbJdGd3Zu+ZuWfm934lk5175t5zvjPZvOfO754zk6pCkrT8HTLuASRJi8PgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLga0lKsiPJaXO43y8l+WKSPUnuT3JdkrWLdfx57vMdST7Z5z6l6Qy+lrujgS3AWuCpwI3AF8c5kDQuBl9jl2R1ks91z8LvS/KhJL+Y5Gvd7XuTXJ7k6O7+nwCOA65J8mCSN82276q6saouqar7q+q/gb8E1iZ58gFmekeSq5J8PMkDSW5Lsm4On8shSTYn+X43+1VJntR9bCJJJdmUZFeSHyT5k+5j64G3Ai/vPqdvd9ufnmRL99PJ9iSvmuuMSd6c5O7uY3ckedGB5tfyZvA1VkkOBb4E3AVMAKuAK4EA5wNPB34ZWA28A6CqXgn8O7Chqh5fVX8xj0O+APhhVd03h/ue2c1yNFM/JXxoDo/5Y+ClwK93s/8HcOGM+/wmcCLwYmBzktOq6ivAu4FPd5/Ts7v7XgFMdvt6GfDuGeHe54zdstV5wHOr6gnAbwM75jC/ljGDr3E7hamYvbGqflpVP6uqb1TV9qr6alX9vKr2AO9jKqIHLcmxTMX3DXN8yDeqamtVPQJ8Anj2gR4AvBr4s6qarKqfM/VN6mVJVky7zzu7z/U7wEeBjbPMuxp4PvDm7utyK3Ax8Mo5zPgI8BjgpCSHVdWOqvr+HD9vLVMGX+O2Grirqh6evjHJU5Jc2S1J/AT4JHDMwR4kyUrgeuCvquqKOT7sh9Pe/0/gsTPCvS/PAD6f5EdJfgTczlR8nzrtPjunvX8XU9/w9uXpwP1V9cCM+6860IxVtR14PVPfcHZ3X8vZjqNGGHyN207guH2E9HyggGdV1ZHAK5ha5tlrzr/mNckTmYr9lqp614jzHshO4PSqOnra22Or6u5p91k97f3jgF3d+zM/p13Ak5I8Ycb972YOqupTVfV8pr4JFfCe+XwiWn4MvsbtRuAHwAVJjkjy2CSnAk8AHgR+lGQV8MYZj7sHOOFAO09yJHAd8A9Vtbnf0Tmsm3fv2wrgI8C7kjyjO/7KJGfNeNyfJ3lckl8BzgE+3W2/B5hIcghAVe0E/hE4v9v/s4BzgcsPNFiStUlemOQxwM+A/2LqJw01zOBrrLq15w3AGqZeiJ0EXg68E3gO8GPgWuBzMx56PvC2bunkT/dziN8Fnguc0539svftuB7G38pUSPe+vQP4AFMvnl6f5AHgW8DzZjzu74DtwN8A762q67vtn+n+vS/JLd37G5l6MXsX8Hng7VX11TnM9hjgAuBeppZ9nsLUWUBqWPwDKNLiSDIB/Btw2MzXLKTF4DN8SWqEwdeSl+QPZizX7H277QCP+/Isj3PpQ8uSSzqS1Aif4UtSIw50EclYHXPMMTUxMTHuMSRpybj55pvvraqV+/rYoIM/MTHBtm3bxj2GJC0ZSe6a7WMu6UhSIwy+JDVikMFPsiHJRT/+8Y/HPYokLRuDDH5VXVNVm4466qhxjyJJy8Yggy9J6p/Bl6RGDPq0TA3TxOZrxz3CnO244IxxjyANxiCf4fuirST1b5DB90VbSerfIIMvSeqfwZekRhh8SWrEIIPvi7aS1L9BBt8XbSWpf4MMviSpfwZfkhph8CWpEQZfkhph8CWpEYMMvqdlSlL/Bhl8T8uUpP4NMviSpP4ZfElqhMGXpEYYfElqhMGXpEYYfElqhMGXpEYYfElqxCCD75W2ktS/QQbfK20lqX+DDL4kqX8GX5IaYfAlqREGX5IaYfAlqREGX5IaYfAlqREGX5IaYfAlqREGX5IaYfAlqREGX5IasWjBT/LSJH+d5ItJXrxYx5UkTVkxyoOTXAq8BNhdVc+ctn098AHgUODiqrqgqr4AfCHJE4H3AtePcmxpLiY2XzvuEeZlxwVnjHsELWOjPsO/DFg/fUOSQ4ELgdOBk4CNSU6adpe3dR+XJC2ikYJfVTcA98/YfAqwvarurKqHgCuBszLlPcCXq+qW2faZZFOSbUm27dmzZ5TxJEnTLMQa/ipg57Tbk9221wKnAS9L8prZHlxVF1XVuqpat3LlygUYT5LaNNIa/iyyj21VVR8EPrgAx5MkzcFCPMOfBFZPu30ssGs+O/Bv2kpS/xYi+DcBJyY5PsnhwNnAlvnswL9pK0n9Gyn4Sa4AvgmsTTKZ5Nyqehg4D7gOuB24qqpuG31USdIoRlrDr6qNs2zfCmw92P0m2QBsWLNmzcHuQpI0wyB/tYJLOpLUv0EGX5LUv0EG37N0JKl/gwy+SzqS1L9BBl+S1D+DL0mNGGTwXcOXpP4NMviu4UtS/wYZfElS/wy+JDVikMF3DV+S+jfI4LuGL0n9G2TwJUn9M/iS1AiDL0mNMPiS1IhBBt+zdCSpf4MMvmfpSFL/Bhl8SVL/DL4kNcLgS1IjDL4kNcLgS1IjBhl8T8uUpP4NMvielilJ/Rtk8CVJ/TP4ktQIgy9JjTD4ktQIgy9JjTD4ktQIgy9JjTD4ktSIQQbfK20lqX8rxj3AvlTVNcA169ate9W4Z5EW08Tma8c9wpztuOCMcY+geRrkM3xJUv8MviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMWLfhJTkhySZKrF+uYkqRHjRT8JJcm2Z3kuzO2r09yR5LtSTYDVNWdVXXuKMeTJB28UZ/hXwasn74hyaHAhcDpwEnAxiQnjXgcSdKIRgp+Vd0A3D9j8ynA9u4Z/UPAlcBZoxxHkjS6hVjDXwXsnHZ7EliV5MlJPgKcnOQtsz04yaYk25Js27NnzwKMJ0ltWog/gJJ9bKuqug94zYEeXFUXARcBrFu3rnqeTZKatRDP8CeB1dNuHwvsms8O/BOHktS/hQj+TcCJSY5PcjhwNrBlPjuoqmuqatNRRx21AONJUptGPS3zCuCbwNokk0nOraqHgfOA64Dbgauq6rbRR5UkjWKkNfyq2jjL9q3A1oPdb5INwIY1a9Yc7C4kSTMM8lcruKQjSf0bZPAlSf0bZPA9S0eS+jfI4LukI0n9G2TwJUn9M/iS1IhBBt81fEnq3yCD7xq+JPVvkMGXJPXP4EtSIwy+JDVikMH3RVtJ6t8gg++LtpLUv0EGX5LUP4MvSY0w+JLUiEEG3xdtJal/gwy+L9pKUv8GGXxJUv8MviQ1wuBLUiMMviQ1wuBLUiMGGXxPy5Sk/g0y+J6WKUn9G2TwJUn9M/iS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1IhBBt8rbSWpf4MMvlfaSlL/Bhl8SVL/DL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNWLFYh0oyRHAXwEPAV+vqssX69iSpBGf4Se5NMnuJN+dsX19kjuSbE+yudv8e8DVVfUq4MxRjitJmr9Rl3QuA9ZP35DkUOBC4HTgJGBjkpOAY4Gd3d0eGfG4kqR5GmlJp6puSDIxY/MpwPaquhMgyZXAWcAkU9G/lf18o0myCdgEcNxxx40ynqQFNLH52nGPMC87Ljhj3COM3UK8aLuKR5/Jw1ToVwGfA34/yYeBa2Z7cFVdVFXrqmrdypUrF2A8SWrTQrxom31sq6r6KXDOAhxPkjQHC/EMfxJYPe32scCu+ezAv2krSf1biODfBJyY5PgkhwNnA1vmswP/pq0k9W/U0zKvAL4JrE0ymeTcqnoYOA+4DrgduKqqbht9VEnSKEY9S2fjLNu3AlsPdr9JNgAb1qxZc7C7kCTNMMhfreCSjiT1b5DBlyT1b5DB9ywdSerfIIPvko4k9W+QwZck9c/gS1IjBhl81/AlqX+DDL5r+JLUv0EGX5LUP4MvSY0YZPBdw5ek/g0y+K7hS1L/Bhl8SVL/DL4kNcLgS1IjDL4kNWKQwfcsHUnq3yCD71k6ktS/QQZfktQ/gy9JjTD4ktQIgy9JjTD4ktSIQQbf0zIlqX+DDL6nZUpS/wYZfElS/wy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDVikMH3SltJ6t8gg++VtpLUv0EGX5LUP4MvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUCIMvSY1YtOAnOSHJJUmuXqxjSpIeNafgJ7k0ye4k352xfX2SO5JsT7J5f/uoqjur6txRhpUkHbwVc7zfZcCHgI/v3ZDkUOBC4LeASeCmJFuAQ4HzZzz+D6tq98jTSpIO2pyCX1U3JJmYsfkUYHtV3QmQ5ErgrKo6H3jJwQ6UZBOwqbv5YJI7DnJXxwD3Huwci2wpzQpLa96lNCssrXmX0qzkPUtq3lFmfcZsH5jrM/x9WQXsnHZ7EnjebHdO8mTgXcDJSd7SfWP4f6rqIuCiEebae7xtVbVu1P0shqU0KyyteZfSrLC05l1Ks8LSmnehZh0l+NnHtprtzlV1H/CaEY4nSRrBKGfpTAKrp90+Ftg12jiSpIUySvBvAk5McnySw4GzgS39jNWLkZeFFtFSmhWW1rxLaVZYWvMupVlhac27ILOmatZVmEfvlFwB/AZTLyTcA7y9qi5J8jvA+5k6M+fSqnrXQgwpSRrdnIIvSVr6/NUKktSIZRf8+Vz9O26zXcE8RElWJ/nbJLcnuS3J68Y90/4keWySG5N8u5v3neOe6UCSHJrkn5N8adyzHEiSHUm+k+TWJNvGPc/+JDk6ydVJvtf9//21cc80myRru6/p3refJHl9b/tfTks63dW//8q0q3+BjVX1L2MdbBZJXgA8CHy8qp457nn2J8nTgKdV1S1JngDcDLx0wF/bAEdU1YNJDgO+Abyuqr415tFmleQNwDrgyKo66IsXF0OSHcC6qhr8hUxJPgb8fVVd3J1g8riq+tGYxzqgrmd3A8+rqrv62Odye4b/f1f/VtVDwJXAWWOeaVZVdQNw/7jnmIuq+kFV3dK9/wBwO1MX3w1STXmwu3lY9zbYZzdJjgXOAC4e9yzLSZIjgRcAlwBU1UNLIfadFwHf7yv2sPyCv6+rfwcbpaWq+zUbJwP/NOZR9qtbIrkV2A18taqGPO/7gTcB/zPmOeaqgOuT3Nz9OpShOgHYA3y0Wy67OMkR4x5qjs4Gruhzh8st+PO6+lfzl+TxwGeB11fVT8Y9z/5U1SNV9atMXRR4SpJBLpsleQmwu6puHvcs83BqVT0HOB34o255cohWAM8BPlxVJwM/BQb92h5At/R0JvCZPve73ILv1b8LqFsL/yxweVV9btzzzFX3I/zXgfXjnWRWpwJnduviVwIvTPLJ8Y60f1W1q/t3N/B5ppZTh2gSmJz2093VTH0DGLrTgVuq6p4+d7rcgj/0q3+XrO5F0EuA26vqfeOe50CSrExydPf+LwCnAd8b61CzqKq3VNWxVTXB1P/Zr1XVK8Y81qySHNG9cE+3PPJiYJBnmlXVD4GdSdZ2m14EDPJEgxk20vNyDoz2y9MGp6oeTnIecB2PXv1725jHmtX0K5iTTNJdwTzeqWZ1KvBK4DvdujjAW6tq6/hG2q+nAR/rznQ4BLiqqgZ/uuMS8VTg81PPAVgBfKqqvjLekfbrtcDl3ZPAO4FzxjzPfiV5HFNnGr66930vp9MyJUmzW25LOpKkWRh8SWqEwZekRhh8SWqEwZekRhh8SWqEwZekRvwvIDJC30vuI6EAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQO0lEQVR4nO3df6zddX3H8edLkDlRygKd07ZaWFkztrlBOtxG4n7ITJkU3I9EmummY1aSYTRuumpMpksI6h/OGNBZheEPRsMQtUgdmjjH3JxQ0EU6ZOkIjCtqK4wibhNh7/1xD3J37bm9t+d77vme+3k+kobeb8/5ft+3gWe//Xy/30OqCknSyvekSQ8gSVoeBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl/NSHJ3krMW8bqfSPKJJAeSPJDkxiQbOzj+55L8YZezSkth8KUfdDywC9gIPAO4GfjEJAeSumDwNZWSrEty3eAs/P4klyb58SSfHXz9rSRXJTl+8PoPA88Grk/ycJI3DNt3Vd1cVZdX1QNV9T3gL4CNSU44zExvSXJNkg8l+XaSvUk2LfD6c5J8OcmDSf4pyXOHzZrkKUk+MvjeHkxyS5JnLPk3Tk0z+Jo6SY4CPgncA6wH1gA7gQCXAM8CfhJYB7wFoKpeBvwHsKWqnlZV71jCIZ8PfKOq7l/Ea88dzHI8s39LuHTI93A6cAXwKuAE4H3AriQ/NGTW3wdWDb6nE4ALgf9ewvcgGXxNpTOYjfrrq+o7VfU/VfX5qtpXVZ+pqu9W1QHgncAvj3KgJGuBy4DXLfItn6+q3VX1GPBh4GeHvO6VwPuq6otV9VhVfRD4LvALQ17/PWZDv2Hw+lur6qHFfyeSwdd0WgfcU1WPzt2Y5EeT7EzytSQPAR8BTjzSgyRZDXwaeE9VXb3It31jzs//C3hKkqMP8brnAH88WJ55MMmDzH5fzxqy3w8DNwI7k9yX5B1JnrzImSTA4Gs63Qs8+xAhvQQo4LlVdRzwUmaXeR636I+GTfIjzMZ+V1VdPOK8h3IvcHFVHT/nx1Pn/MHy/2atqu9V1Vur6lTgl4BzgN8bw1xawQy+ptHNwNeBtyU5dnBB80zg6cDDwINJ1gCvn/e+bwInH27nSY5j9mz6H6tqe7ejf9/7gQuTPC+zjk3yoiRPP9SsSX41yc8Mrl88xOwSz2Njmk0rlMHX1Bmsj28BNjB7cXMGeAnwVuB04CBwA3DdvLdeArx5sITyJwsc4jeBnwdeMbhL5vEfz+5i/MH3sIfZdfxLgf8E9gEvX2DWHwOuZTb2dwB/z+ySlbRo8X+AIi2PJLcBf15VH5/0LGqTZ/jSMkjyU8zeKvqlSc+idhl8NSnJ785brnn8x97DvO9TQ973pgXe83ZmLwD/aVXd0/X3Ii2WSzqS1AjP8CWpEYd6IKQ3TjzxxFq/fv2kx5CkqXHrrbd+q6pWH+rXeh389evXs2fPnkmPIUlTI8nQ60Qu6UhSIwy+JDXC4EtSIwy+JDXC4EtSI3p9l47UmvXbb5j0CIt299teNOkRtES9PMNPsiXJjoMHD056FElaMXoZ/Kq6vqq2rVq1atKjSNKK0cvgS5K6Z/AlqREGX5IaYfAlqREGX5Ia4X34ko7IND0zAD43AJ7hS1IzDL4kNcLgS1IjDL4kNcLgS1IjvEtHUhOm6a6icd1RZPC1ok3Tf+TSuLmkI0mNWLbgJ3lxkvcn+USSFy7XcSVJs0YKfpIrkuxPcvu87ZuT3JlkX5LtAFX18ap6JfBy4CWjHFeStHSjnuFfCWyeuyHJUcBlwNnAqcDWJKfOecmbB78uSVpGIwW/qm4CHpi3+QxgX1XdVVWPADuB8zLr7cCnquq2YftMsi3JniR7Dhw4MMp4kqQ5xrGGvwa4d87XM4NtrwbOAn4nyYXD3lxVO6pqU1VtWr169RjGk6Q2jeO2zBxiW1XVu4F3j+F4kqRFGMcZ/gywbs7Xa4H7lrKDJFuS7Dh48GCng0lSy8YR/FuAU5KclOQY4Hxg11J2UFXXV9W2VatWjWE8SWrTqLdlXg18AdiYZCbJBVX1KHARcCNwB3BNVe0dfVRJ0ihGWsOvqq1Dtu8Gdh/pfpNsAbZs2LDhSHchSZqnlx+t4JKOJHWvl8GXJHXP4EtSI3oZfG/LlKTu9TL4ruFLUvd6GXxJUvcMviQ1opfBdw1fkrrXy+C7hi9J3etl8CVJ3TP4ktQIgy9Jjehl8L1oK0nd62XwvWgrSd3rZfAlSd0z+JLUCIMvSY0w+JLUiF4G37t0JKl7vQy+d+lIUvd6GXxJUvcMviQ1wuBLUiMMviQ1wuBLUiMMviQ1opfB9z58SepeL4PvffiS1L1eBl+S1D2DL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1IheBt+PVpCk7vUy+H60giR1r5fBlyR1z+BLUiMMviQ1wuBLUiMMviQ1wuBLUiOOnvQAmj7rt98w6REkHQHP8CWpEQZfkhph8CWpEQZfkhph8CWpEQZfkhph8CWpEcsW/CQnJ7k8ybXLdUxJ0hNGCn6SK5LsT3L7vO2bk9yZZF+S7QBVdVdVXTDK8SRJR27UM/wrgc1zNyQ5CrgMOBs4Fdia5NQRjyNJGtFIwa+qm4AH5m0+A9g3OKN/BNgJnDfKcSRJoxvHGv4a4N45X88Aa5KckOQvgdOSvHHYm5NsS7InyZ4DBw6MYTxJatM4Pjwth9hWVXU/cOHh3lxVO4AdAJs2baqOZ5OkZo3jDH8GWDfn67XAfUvZQZItSXYcPHiw08EkqWXjCP4twClJTkpyDHA+sGspO6iq66tq26pVq8YwniS1adTbMq8GvgBsTDKT5IKqehS4CLgRuAO4pqr2jj6qJGkUI63hV9XWIdt3A7tH2bckqVu9/GgF1/AlqXu9DL5r+JLUvV4GX5LUvV4G3yUdSepeL4Pvko4kda+XwZckdc/gS1IjDL4kNaKXwfeirSR1r5fB96KtJHWvl8GXJHXP4EtSIwy+JDWil8H3oq0kda+XwfeirSR1r5fBlyR1z+BLUiMMviQ1wuBLUiN6GXzv0pGk7vUy+N6lI0nd62XwJUndM/iS1AiDL0mNMPiS1AiDL0mNMPiS1IheBt/78CWpe70MvvfhS1L3ehl8SVL3DL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNcLgS1Ijjp70AIeSZAuwZcOGDZMeZVms337DpEeQ1IBenuH70QqS1L1eBl+S1D2DL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNWLbPw09yLPAe4BHgc1V11XIdW5I04hl+kiuS7E9y+7ztm5PcmWRfku2Dzb8FXFtVrwTOHeW4kqSlG3VJ50pg89wNSY4CLgPOBk4FtiY5FVgL3Dt42WMjHleStEQjBb+qbgIemLf5DGBfVd1VVY8AO4HzgBlmo7/gcZNsS7InyZ4DBw6MMp4kaY5xXLRdwxNn8jAb+jXAdcBvJ3kvcP2wN1fVjqraVFWbVq9ePYbxJKlN47hom0Nsq6r6DvCKMRxPkrQI4zjDnwHWzfl6LXDfUnaQZEuSHQcPHux0MElq2TiCfwtwSpKTkhwDnA/sWsoOqur6qtq2atWqMYwnSW0a9bbMq4EvABuTzCS5oKoeBS4CbgTuAK6pqr2jjypJGsVIa/hVtXXI9t3A7iPdb5ItwJYNGzYc6S4kSfP08qMVXNKRpO71MviSpO4ZfElqRC+D722ZktS9XgbfNXxJ6l4vgy9J6p7Bl6RG9DL4ruFLUvd6GXzX8CWpe70MviSpewZfkhph8CWpEb0MvhdtJal7vQy+F20lqXu9DL4kqXsGX5IaYfAlqREGX5Ia0cvge5eOJHWvl8H3Lh1J6l4vgy9J6p7Bl6RGGHxJaoTBl6RGGHxJaoTBl6RG9DL43ocvSd3rZfC9D1+SutfL4EuSumfwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRvQy+n6UjSd3rZfD9LB1J6l4vgy9J6p7Bl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGLFvwk5yc5PIk1y7XMSVJT1hU8JNckWR/ktvnbd+c5M4k+5JsX2gfVXVXVV0wyrCSpCN39CJfdyVwKfChxzckOQq4DPh1YAa4Jcku4Cjgknnv/4Oq2j/ytJKkI7ao4FfVTUnWz9t8BrCvqu4CSLITOK+qLgHOOdKBkmwDtg2+fDjJnUe4qxOBbx3pHMtsmmaF6Zp3mmaF6Zp3mmaFKZo3bx9p1ucM+4XFnuEfyhrg3jlfzwDPG/biJCcAFwOnJXnj4A+GH1BVO4AdI8z1+PH2VNWmUfezHKZpVpiueadpVpiueadpVpiuecc16yjBzyG21bAXV9X9wIUjHE+SNIJR7tKZAdbN+XotcN9o40iSxmWU4N8CnJLkpCTHAOcDu7oZqxMjLwsto2maFaZr3mmaFaZr3mmaFaZr3rHMmqqhqzBPvCi5GvgVZi96fBP4s6q6PMlvAO9i9s6cK6rq4nEMKUka3aKCL0mafn60giQ1YsUFfylP/07asCeY+yjJuiR/l+SOJHuTvGbSMy0kyVOS3JzkXwbzvnXSMx1OkqOSfCnJJyc9y+EkuTvJV5J8OcmeSc+zkCTHJ7k2yVcH//7+4qRnGibJxsHv6eM/Hkry2s72v5KWdAZP//4bc57+BbZW1b9OdLAhkjwfeBj4UFX99KTnWUiSZwLPrKrbkjwduBV4cY9/bwMcW1UPJ3ky8HngNVX1zxMebagkrwM2AcdV1RE/vLgcktwNbKqq3j/IlOSDwD9U1QcGN5g8taoenPBYhzXo2deA51XVPV3sc6Wd4X//6d+qegTYCZw34ZmGqqqbgAcmPcdiVNXXq+q2wc+/DdzB7MN3vVSzHh58+eTBj96e3SRZC7wI+MCkZ1lJkhwHPB+4HKCqHpmG2A+8APj3rmIPKy/4h3r6t7dRmlaDj9k4DfjihEdZ0GCJ5MvAfuAzVdXned8FvAH43wnPsVgFfDrJrYOPQ+mrk4EDwF8Nlss+kOTYSQ+1SOcDV3e5w5UW/CU9/aulS/I04KPAa6vqoUnPs5Cqeqyqfo7ZhwLPSNLLZbMk5wD7q+rWSc+yBGdW1enA2cAfDZYn++ho4HTgvVV1GvAdoNfX9gAGS0/nAn/T5X5XWvB9+neMBmvhHwWuqqrrJj3PYg3+Cv85YPNkJxnqTODcwbr4TuDXknxksiMtrKruG/xzP/AxZpdT+2gGmJnzt7trmf0DoO/OBm6rqm92udOVFvy+P/07tQYXQS8H7qiqd056nsNJsjrJ8YOf/zBwFvDViQ41RFW9sarWVtV6Zv+d/WxVvXTCYw2V5NjBhXsGyyMvBHp5p1lVfQO4N8nGwaYXAL280WCerXS8nAOjfXha71TVo0kuAm7kiad/9054rKHmPsGcZIbBE8yTnWqoM4GXAV8ZrIsDvKmqdk9upAU9E/jg4E6HJwHXVFXvb3ecEs8APjZ7DsDRwF9X1d9OdqQFvRq4anASeBfwignPs6AkT2X2TsNXdb7vlXRbpiRpuJW2pCNJGsLgS1IjDL4kNcLgS1IjDL4kNcLgS1IjDL4kNeL/AIki2k4n0eHoAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAR2ElEQVR4nO3df5BdZX3H8feHhB+KIFWio0k0UJiMqaPF2UJbHGuV2qBCHOuMpOpMlQGZKVanrRqtHfEPBTsdxzpSnQiIPzAREdqgUXS0lNKiEqhWaKSNFJolagIUBFtlwG//uIeybrPJ3dyze8/ueb9mdnLvc895zvfuZD737HOe59xUFZKkxe+gcRcgSZofBr4k9YSBL0k9YeBLUk8Y+JLUEwa+JPWEgS9JPWHga0FKckeSU4bY7ugk/5jkniT3JbkhycnzdfxZ9nlekk+32ac01dJxFyDNsQeBNwD/DhSwDrg6yVOq6uGxVibNM8/wNXZJVia5Msme5kz8w0l+OcnXm+d3J7ksyVHN9p8CnsEguB9M8raZ+q6qn1bVbVX1cyDAI8AvAU/aT03nJbk8ySeTPJDk1iQTQ7yXg5JsSPL9pvbLkzypeW1VkkpydpJdSX6Q5E+a19YC7wRe3byn7zTtT0+yJcm9SXYkOWvYGpO8PcldzWu3JXnx/urX4mbga6ySLAG+ANwJrAKWA5sZhPP5wNOBZwErgfMAqup1wH8Cp1XVE6rqL4Y4zr8APwW2ABdV1e4hyju9qeWoZr8PD7HPHwGvAH6rqf2/gAunbfPbwPHAS4ANSU6pqi8D7wM+27yn5zbbbgImm75eBbxvWnDvtcYkq4FzgV+rqiOA3wXuGKJ+LWIGvsbtRAZh9taq+klzRn59Ve2oqq9W1c+qag/wAQYhekCq6jnAkcDvA9cPudv1VbW1qh4BPgU8d387AG8E/qyqJqvqZww+pF6VZOrw6Xua9/pd4OPA+r11lGQl8Hzg7c3v5dvARcDrhqjxEeBQYE2Sg6vqjqr6/pDvW4uUga9xWwncOX08PclTkmxuhiR+DHwaOHqUAzWhuYnBWfUw4f3DKY//GzhsWnDvzTOBq5oLxPcB2xmE71OnbLNzyuM7GXzg7c3TgXur6oFp2y/fX41VtQN4C4MPnN3N73Km46gnDHyN207gGXsJ0vMZXGR9TlUdCbyWwTDPo0a5zevBwLEj7L8vO4FTq+qoKT+HVdVdU7ZZOeXxM4BdzePp72kX8KQkR0zb/i6GUFWfqarnM/gQKuD9s3kjWnwMfI3bt4AfABckOTzJYc20ySMYzLC5L8ly4K3T9vsRQ4R2kl9P8vwkhyR5XJK3Mzjb/mYLtR/c1Pvoz1Lgo8B7kzyzOf6yJOum7ffnSR6f5FeA1wOfnfKeViU5CKCqdgL/BJzf9P8c4EzgsiHe9+okL0pyKINrF//D4C8N9ZiBr7Fqxp5PA45jcCF2Eng18B7gecD9wBeBK6ftej7wrmbo5E/3cYhDGVw0vYfBmfFLgZdV1a597DOsrQyC9NGf84C/YnDx9CtJHgC+AZw0bb+/B3YAXwP+sqq+0rR/rvn3niQ3N4/XM7iYvQu4Cnh3VX11iNoOBS4A7mYw7PMUBrOA1GPxC1Ck+ZFkFfAfwMGuAdA4eIYvST1h4GvBS/KaZrHS9J9b97Pfl2bYz6EPLUoO6UhST3iGL0k9MW83T0vyLODNDBbPfK2qPrK/fY4++uhatWrVXJcmSYvGTTfddHdVLdvbayMFfpJLgJcDu6vq2VPa1zKYnraEwX1LLqiq7cA5zRzjjw3T/6pVq9i2bdsoJUpSryS5c6bXRh3SuRRYO+1gSxjMez4VWAOsT7Kmee10Bvcx+dqIx5UkzdJIgV9V1wH3Tms+EdhRVbdX1UMM7uS3rtl+S1X9JvCamfpsbh27Lcm2PXv2jFKeJGmKuRjDX84v3hxqEjgpyQuBVzJYAbh1pp2raiOwEWBiYsIpRJLUkrkI/OylrarqWuDaOTieJGkIczEtc5JfvBvgCh67G+BQkpyWZOP999/famGS1GdzEfg3AscnOSbJIcAZDG4mNbSqurqqzn7iE584B+VJUj+NFPhJNgE3AKuTTCY5s7kp1LnANQy+/OHyqtrnEve99OsZviS1rNO3VpiYmCjn4UvS8JLcVFUTe3vNWytIUk90MvAd0pGk9nUy8L1oK0nt62TgS5LaZ+BLUk90MvAdw5ek9nUy8B3Dl6T2dTLwJUntM/AlqSc6GfiO4UtS+zoZ+I7hS1L7Ohn4kqT2GfiS1BMGviT1RCcD34u2ktS+Tga+F20lqX2dDHxJUvsMfEnqCQNfknrCwJeknjDwJaknOhn4TsuUpPZ1MvCdlilJ7etk4EuS2mfgS1JPGPiS1BMGviT1hIEvST1h4EtSTywddwFzZdWGL467hKHdccHLxl2CpB7wDF+SeqKTge9KW0lqXycD35W2ktS+Tga+JKl9Br4k9YSBL0k9YeBLUk8Y+JLUEwa+JPWEgS9JPWHgS1JPGPiS1BMGviT1xLwFfpJXJPlYkr9N8pL5Oq4kaWCkwE9ySZLdSW6Z1r42yW1JdiTZAFBVf1NVZwF/ALx6lONKkmZv1DP8S4G1UxuSLAEuBE4F1gDrk6yZssm7mtclSfNopMCvquuAe6c1nwjsqKrbq+ohYDOwLgPvB75UVTePclxJ0uzNxRj+cmDnlOeTTdubgFOAVyU5Z6adk5ydZFuSbXv27JmD8iSpn+biKw6zl7aqqg8BH9rfzlW1EdgIMDExUS3XJkm9NRdn+JPAyinPVwC7ZtOB33glSe2bi8C/ETg+yTFJDgHOALbMpgO/8UqS2jfqtMxNwA3A6iSTSc6sqoeBc4FrgO3A5VV16+ilSpJGMdIYflWtn6F9K7D1QPtNchpw2nHHHXegXUiSpunkrRUc0pGk9nUy8CVJ7etk4DtLR5La18nAd0hHktrXycCXJLXPwJeknuhk4DuGL0nt62TgO4YvSe3rZOBLktpn4EtST3Qy8B3Dl6T2dTLwHcOXpPZ1MvAlSe0z8CWpJwx8SeqJTga+F20lqX2dDHwv2kpS+zoZ+JKk9hn4ktQTBr4k9YSBL0k9YeBLUk90MvCdlilJ7etk4DstU5La18nAlyS1z8CXpJ4w8CWpJwx8SeoJA1+SesLAl6SeMPAlqSc6GfguvJKk9nUy8F14JUnt62TgS5LaZ+BLUk8Y+JLUEwa+JPWEgS9JPWHgS1JPGPiS1BMGviT1hIEvST1h4EtSTxj4ktQT8xb4SY5NcnGSK+brmJKkx4wU+EkuSbI7yS3T2tcmuS3JjiQbAKrq9qo6c5TjSZIO3Khn+JcCa6c2JFkCXAicCqwB1idZM+JxJEkjGinwq+o64N5pzScCO5oz+oeAzcC6YftMcnaSbUm27dmzZ5TyJElTzMUY/nJg55Tnk8DyJE9O8lHghCTvmGnnqtpYVRNVNbFs2bI5KE+S+mnpHPSZvbRVVd0DnDMHx5MkDWEuzvAngZVTnq8Ads2mA7/iUJLaNxeBfyNwfJJjkhwCnAFsmU0HfsWhJLVv1GmZm4AbgNVJJpOcWVUPA+cC1wDbgcur6tbRS5UkjWKkMfyqWj9D+1Zg64H2m+Q04LTjjjvuQLuQJE3TyVsrOKQjSe3rZOBLktrXycB3lo4kta+Tge+QjiS1r5OBL0lqXycD3yEdSWpfJwPfIR1Jal8nA1+S1D4DX5J6wsCXpJ7oZOB70VaS2tfJwPeirSS1r5OBL0lqn4EvST1h4EtST3Qy8L1oK0nt62Tge9FWktrXycCXJLXPwJeknjDwJaknDHxJ6olOBr6zdCSpfZ0MfGfpSFL7Ohn4kqT2GfiS1BMGviT1hIEvST1h4EtSTxj4ktQTBr4k9YSBL0k90cnAd6WtJLWvk4HvSltJal8nA1+S1D4DX5J6wsCXpJ4w8CWpJwx8SeoJA1+SesLAl6SeMPAlqScMfEnqCQNfknrCwJeknlg6XwdKcjjw18BDwLVVddl8HVuSNOIZfpJLkuxOcsu09rVJbkuyI8mGpvmVwBVVdRZw+ijHlSTN3qhDOpcCa6c2JFkCXAicCqwB1idZA6wAdjabPTLicSVJszRS4FfVdcC905pPBHZU1e1V9RCwGVgHTDII/X0eN8nZSbYl2bZnz55RypMkTTEXF22X89iZPAyCfjlwJfB7ST4CXD3TzlW1saomqmpi2bJlc1CeJPXTXFy0zV7aqqp+Arx+Do4nSRrCXJzhTwIrpzxfAeyaTQd+xaEktW8uAv9G4PgkxyQ5BDgD2DKbDvyKQ0lq36jTMjcBNwCrk0wmObOqHgbOBa4BtgOXV9Wts+zXM3xJatlIY/hVtX6G9q3A1hH6vRq4emJi4qwD7UOS9Iu8tYIk9UQnA98hHUlqXycD34u2ktS+Tga+JKl9Br4k9UQnA98xfElqXycD3zF8SWpfJwNfktQ+A1+SeqKTge8YviS1r5OB7xi+JLWvk4EvSWqfgS9JPWHgS1JPdDLwvWgrSe3rZOB70VaS2tfJwJcktc/Al6SeMPAlqScMfEnqCQNfknqik4HvtExJal8nA99pmZLUvk4GviSpfQa+JPWEgS9JPWHgS1JPGPiS1BMGviT1hIEvST1h4EtST3Qy8F1pK0nt62Tgu9JWktrXycCXJLXPwJeknjDwJaknDHxJ6gkDX5J6wsCXpJ4w8CWpJwx8SeoJA1+SesLAl6SeMPAlqSfmLfCTHJvk4iRXzNcxJUmPGSrwk1ySZHeSW6a1r01yW5IdSTbsq4+qur2qzhylWEnSgVs65HaXAh8GPvloQ5IlwIXA7wCTwI1JtgBLgPOn7f+Gqto9crWSpAM2VOBX1XVJVk1rPhHYUVW3AyTZDKyrqvOBlx9oQUnOBs5unj6Y5LYD7Opo4O4DrWM+5f0Lp9bGQqp3IdUKC6vehVQrLKx6R6n1mTO9MOwZ/t4sB3ZOeT4JnDTTxkmeDLwXOCHJO5oPhv+nqjYCG0eo69HjbauqiVH7mQ8LqVZYWPUupFphYdW7kGqFhVXvXNU6SuBnL20108ZVdQ9wzgjHkySNYJRZOpPAyinPVwC7RitHkjRXRgn8G4HjkxyT5BDgDGBLO2W1YuRhoXm0kGqFhVXvQqoVFla9C6lWWFj1zkmtqZpxFOaxjZJNwAsZXEj4EfDuqro4yUuBDzKYmXNJVb13LoqUJI1uqMCXJC183lpBknpi0QX+bFb/jttMK5i7KMnKJH+XZHuSW5O8edw17UuSw5J8K8l3mnrfM+6a9ifJkiT/nOQL465lf5LckeS7Sb6dZNu469mXJEcluSLJ95r/v78x7ppmkmR18zt99OfHSd7SWv+LaUinWf37b0xZ/Qusr6p/HWthM0jyAuBB4JNV9exx17MvSZ4GPK2qbk5yBHAT8IoO/24DHF5VDyY5GLgeeHNVfWPMpc0oyR8DE8CRVXXAixfnQ5I7gImq6vxCpiSfAP6hqi5qJpg8vqruG3NZ+9Xk2V3ASVV1Zxt9LrYz/P9b/VtVDwGbgXVjrmlGVXUdcO+46xhGVf2gqm5uHj8AbGew+K6TauDB5unBzU9nz26SrABeBlw07loWkyRHAi8ALgaoqocWQtg3Xgx8v62wh8UX+Htb/dvZUFqomttsnAB8c8yl7FMzRPJtYDfw1arqcr0fBN4G/HzMdQyrgK8kuam5HUpXHQvsAT7eDJddlOTwcRc1pDOATW12uNgCf1arfzV7SZ4AfB54S1X9eNz17EtVPVJVv8pgUeCJSTo5bJbk5cDuqrpp3LXMwslV9TzgVOAPm+HJLloKPA/4SFWdAPwE6PS1PYBm6Ol04HNt9rvYAt/Vv3OoGQv/PHBZVV057nqG1fwJfy2wdryVzOhk4PRmXHwz8KIknx5vSftWVbuaf3cDVzEYTu2iSWByyl93VzD4AOi6U4Gbq+pHbXa62AK/66t/F6zmIujFwPaq+sC469mfJMuSHNU8fhxwCvC9sRY1g6p6R1WtqKpVDP7Pfr2qXjvmsmaU5PDmwj3N8MhLgE7ONKuqHwI7k6xuml4MdHKiwTTraXk4B0a7eVrnVNXDSc4FruGx1b+3jrmsGU1dwZxkkmYF83irmtHJwOuA7zbj4gDvrKqt4ytpn54GfKKZ6XAQcHlVdX664wLxVOCqwTkAS4HPVNWXx1vSPr0JuKw5CbwdeP2Y69mnJI9nMNPwja33vZimZUqSZrbYhnQkSTMw8CWpJwx8SeoJA1+SesLAl6SeMPAlqScMfEnqif8FPLu+6TV5Jc4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQHklEQVR4nO3df6zddX3H8efLVvyBUjZQ59pqYRCybnODNOjG5n7oTBEKZlsi3XT+YFaSYTTLVHQm0yUG3R/GGJhageEP1oYhuiJ1aOIcY2NKQZx0la0jsF5RW2EWcVME3/vjHOTurqe9t+d77/mefp6P5Iaeb8/5ft+3Ic9+7+d8v6epKiRJR77HTXoASdLSMPiS1AiDL0mNMPiS1AiDL0mNMPiS1AiDL0mNMPhqRpK7k7xwHs87Psk/JrkvybeT3JzkjA6O//kkf9DlrNJCLJ/0AFIPPQi8Gvh3oIBzgeuSPL2qHp7oZNIYPMPXVEqyOsm1SfYNz8QvSfJTST43fPytJFclOXb4/I8Cz2IQ7geTvGnUvqvqe1V1Z1X9EAjwCPBjwI8fYqa3J7k6yUeSfCfJziTrDvL8s5PcPvwp4p+SPGfUrEmemORjs37quCXJMxb4x6bGGXxNnSTLgE8B9wBrgJXAVgZxvhj4SeCngdXA2wGq6uXAfwIbquopVfXn8zjOvwDfA7YBl1XV3nmMd85wlmOHr7tkxL5PA64AXgscB3wQ2JbkCSNmfQWwYvg9HQdcAPzPPOaRfsTgaxqdziDqb6yq7w7PyG+qqt1V9dmq+n5V7QPeA/zq4R6kqp4DHAP8LnDTPF92U1Vtr6pHgI8CPz/iea8BPlhVX6iqR6rqw8D3geeNeP4PGIT+pOHzb62qB+b9zUi4hq/ptBq4Z+56epKnA+8DfgV4KoMTmv8a50BV9T1gS5JdSW6vqi8f4iXfmPXr/waemGT5Adb+nw28IsnrZm07isFfZAfyUQbf99bhMtXHgD+pqh/M93uRPMPXNNoDPCvJ3BOWixm8yfqcqjoGeBmDZZ5HjfPRsI8HThzj9XPtAd5ZVcfO+npyVW0Z/v7/mbWqflBV76iqtcAvAWcDv9/hPGqAwdc0+iLwdeBdSY4evqF5BoOz+geBbydZCbxxzuu+yTyineR5SX45yVFJnpTkzcAzgC90+D18CLggyXMzcHSSs5I89UCzJvn1JD83fP/iAQZLPI90OI8aYPA1dYbr4xuAkxi8uTkDvBR4B3AasB+4Hrh2zksvBt42vMrljw9yiCcAlwL3AV8DXgycVVX3djH+8HvYwWAd/xIGy067gVceZNafAK5hEPtdwN8zWNaR5i3+AyjS0khyG/BnVfXJSc+iNnmGLy2BJD/D4FLRL016FrXL4KtJSX5veFPT3K+dh3jdp0e87q0Hec27gc8Ab66qe7r+XqT5cklHkhrhGb4kNaLXN14df/zxtWbNmkmPIUlT49Zbb/1WVT3tQL/X6+CvWbOGHTt2THoMSZoaSUa+T+SSjiQ1opfBT7Ihyeb9+/dPehRJOmL0MvhVdV1VbVqxYsWkR5GkI0Yvgy9J6p7Bl6RG9PoqHWlcay66ftIjLMjd7zpr0iPoCOYZviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1wuBLUiMMviQ1YsmCn+QlST6U5G+SvGipjitJGhgr+EmuSLI3yR1ztq9PcmeS3UkuAqiqT1bVa4BXAi8d57iSpIUb9wz/SmD97A1JlgGXAmcCa4GNSdbOesrbhr8vSVpCYwW/qm4E7p+z+XRgd1XdVVUPAVuBczPwbuDTVXXbOMeVJC3cYqzhrwT2zHo8M9z2OuCFwO8kuWDUi5NsSrIjyY59+/YtwniS1KbF+AdQcoBtVVXvA953qBdX1WZgM8C6deuq49kkqVmLcYY/A6ye9XgVcO9CdpBkQ5LN+/fv73QwSWrZYgT/FuDkJCckOQo4D9i2kB1U1XVVtWnFihWLMJ4ktWncyzK3ADcDpySZSXJ+VT0MXAjcAOwCrq6qneOPKkkax1hr+FW1ccT27cD2cfYtSepWLz9awTV8SepeL4PvGr4kda+XwZckdc/gS1Ijehl81/AlqXu9DL5r+JLUvV4GX5LUPYMvSY3oZfBdw5ek7vUy+K7hS1L3ehl8SVL3DL4kNcLgS1Ijehl837SVpO71Mvi+aStJ3etl8CVJ3TP4ktQIgy9JjTD4ktSIXgbfq3QkqXu9DL5X6UhS93oZfElS9wy+JDVi+aQH0PRZc9H1kx5B0mHwDF+SGmHwJakRBl+SGmHwJakRvQy+N15JUvd6GXxvvJKk7vUy+JKk7hl8SWqEwZekRhh8SWqEwZekRhh8SWqEwZekRhh8SWqEwZekRvQy+H60giR1r5fB96MVJKl7vQy+JKl7Bl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGmHwJakRBl+SGrFkwU9yYpLLk1yzVMeUJD1mrOAnuSLJ3iR3zNm+PsmdSXYnuQigqu6qqvPHOZ4k6fCNe4Z/JbB+9oYky4BLgTOBtcDGJGvHPI4kaUxjBb+qbgTun7P5dGD38Iz+IWArcO44x5EkjW8x1vBXAntmPZ4BViY5LskHgFOTvGXUi5NsSrIjyY59+/YtwniS1Kbli7DPHGBbVdV9wAWHenFVbQY2A6xbt646nk2SmrUYZ/gzwOpZj1cB9y7CcSRJC7AYwb8FODnJCUmOAs4Dti1kB0k2JNm8f//+RRhPkto07mWZW4CbgVOSzCQ5v6oeBi4EbgB2AVdX1c6F7LeqrquqTStWrBhnPEnSLGOt4VfVxhHbtwPbx9m3JKlbvfxoBZd0JKl7vQy+SzqS1L1eBl+S1D2DL0mN6GXwXcOXpO71Mviu4UtS93oZfElS9wy+JDWil8F3DV+SutfL4LuGL0nd62XwJUndM/iS1AiDL0mN6GXwfdNWkrrXy+D7pq0kda+XwZckdc/gS1IjDL4kNcLgS1IjDL4kNaKXwfeyTEnqXi+D72WZktS9XgZfktQ9gy9JjTD4ktQIgy9JjTD4ktQIgy9Jjehl8L0OX5K618vgex2+JHWvl8GXJHXP4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDVi+aQHOJAkG4ANJ5100qRHWRJrLrp+0iNIakAvz/D9aAVJ6l4vgy9J6p7Bl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGGHxJaoTBl6RGLNnn4Sc5GvgL4CHg81V11VIdW5I05hl+kiuS7E1yx5zt65PcmWR3kouGm38LuKaqXgOcM85xJUkLN+6SzpXA+tkbkiwDLgXOBNYCG5OsBVYBe4ZPe2TM40qSFmis4FfVjcD9czafDuyuqruq6iFgK3AuMMMg+mMfV5K0cIsR3pU8diYPg9CvBK4FfjvJ+4HrRr04yaYkO5Ls2Ldv3yKMJ0ltWow3bXOAbVVV3wVedagXV9VmYDPAunXrquPZJKlZi3GGPwOsnvV4FXDvIhxHkrQAixH8W4CTk5yQ5CjgPGDbQnaQZEOSzfv371+E8SSpTeNelrkFuBk4JclMkvOr6mHgQuAGYBdwdVXtXMh+q+q6qtq0YsWKccaTJM0y1hp+VW0csX07sH2cfUuSutXLyyNd0pGk7vUy+C7pSFL3ehl8SVL3DL4kNaKXwXcNX5K618vgu4YvSd3rZfAlSd0z+JLUiF4G3zV8SepeL4PvGr4kda+XwZckdc/gS1IjDL4kNaKXwfdNW0nqXi+D75u2ktS9XgZfktQ9gy9JjTD4ktQIgy9JjTD4ktSIXgbfyzIlqXu9DL6XZUpS93oZfElS9wy+JDXC4EtSIwy+JDXC4EtSIwy+JDWil8H3OnxJ6l4vg+91+JLUvV4GX5LUPYMvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUCIMvSY0w+JLUiF4G349WkKTu9TL4frSCJHWvl8GXJHXP4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDXC4EtSIwy+JDViyYKf5MQklye5ZqmOKUl6zLyCn+SKJHuT3DFn+/okdybZneSig+2jqu6qqvPHGVaSdPiWz/N5VwKXAB95dEOSZcClwG8CM8AtSbYBy4CL57z+1VW1d+xpJUmHbV7Br6obk6yZs/l0YHdV3QWQZCtwblVdDJx9uAMl2QRsGj58MMmdh7mr44FvHe4cS2yaZoXpmneaZiXvnqp5p2lWmK55x5n12aN+Y75n+AeyEtgz6/EM8NxRT05yHPBO4NQkbxn+xfD/VNVmYPMYcz16vB1VtW7c/SyFaZoVpmveaZoVpmveaZoVpmvexZp1nODnANtq1JOr6j7ggjGOJ0kawzhX6cwAq2c9XgXcO944kqTFMk7wbwFOTnJCkqOA84Bt3YzVibGXhZbQNM0K0zXvNM0K0zXvNM0K0zXvosyaqpGrMI89KdkC/BqDNxK+CfxpVV2e5MXAexlcmXNFVb1zMYaUJI1vXsGXJE0/P1pBkhpxxAV/IXf/TtqoO5j7KMnqJH+XZFeSnUleP+mZDibJE5N8McmXh/O+Y9IzHUqSZUm+lORTk57lUJLcneQrSW5PsmPS8xxMkmOTXJPkq8P/f39x0jONkuSU4Z/po18PJHlDZ/s/kpZ0hnf//huz7v4FNlbVv050sBGSPB94EPhIVf3spOc5mCTPBJ5ZVbcleSpwK/CSHv/ZBji6qh5M8njgJuD1VfXPEx5tpCR/BKwDjqmqw755cSkkuRtYV1W9v5EpyYeBf6iqy4YXmDy5qr494bEOadizrwHPrap7utjnkXaG/6O7f6vqIWArcO6EZxqpqm4E7p/0HPNRVV+vqtuGv/4OsIvBzXe9VAMPDh8+fvjV27ObJKuAs4DLJj3LkSTJMcDzgcsBquqhaYj90AuA/+gq9nDkBf9Ad//2NkrTavgxG6cCX5jwKAc1XCK5HdgLfLaq+jzve4E3AT+c8BzzVcBnktw6/DiUvjoR2Af85XC57LIkR096qHk6D9jS5Q6PtOAv6O5fLVySpwAfB95QVQ9Mep6DqapHquoXGNwUeHqSXi6bJTkb2FtVt056lgU4o6pOA84E/nC4PNlHy4HTgPdX1anAd4Fev7cHMFx6Ogf46y73e6QF37t/F9FwLfzjwFVVde2k55mv4Y/wnwfWT3aSkc4Azhmui28FfiPJxyY70sFV1b3D/+4FPsFgObWPZoCZWT/dXcPgL4C+OxO4raq+2eVOj7Tg9/3u36k1fBP0cmBXVb1n0vMcSpKnJTl2+OsnAS8EvjrRoUaoqrdU1aqqWsPg/9nPVdXLJjzWSEmOHr5xz3B55EVAL680q6pvAHuSnDLc9AKglxcazLGRjpdzYLwPT+udqno4yYXADTx29+/OCY810uw7mJPMMLyDebJTjXQG8HLgK8N1cYC3VtX2yY10UM8EPjy80uFxwNVV1fvLHafEM4BPDM4BWA78VVX97WRHOqjXAVcNTwLvAl414XkOKsmTGVxp+NrO930kXZYpSRrtSFvSkSSNYPAlqREGX5IaYfAlqREGX5IaYfAlqREGX5Ia8b+eTtqQLVGQJAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "for cat in cats:\n",
    "    for obj in objs:\n",
    "        quantity = \"cat_%s_n%s\"%(cat, obj)\n",
    "        exec(\"vals = %s\"%(quantity))\n",
    "        fig, ax = plt.subplots()\n",
    "        plt.hist(vals, bins = np.linspace(0, 7, 8))\n",
    "        plt.yscale('log')\n",
    "        ymin, ymax = ax.get_ylim()\n",
    "        ax.set_ylim(0.1, ymax)\n",
    "        plt.title(quantity)\n",
    "        plt.show()\n",
    "        plt.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}