Deane_CERNHSTWG5_Sample2.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Deane_CERNHSTWG5_Sample2\n",
    "\n",
    "Deane, J.K.\n",
    "\n",
    "2018-07-09\n",
    "\n",
    "An example Jupyter notebook to upload to Github."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# This notebook runs on Python 3, so we need to import some modules.\n",
    "\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from scipy.stats import norm\n",
    "\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "muons = pd.read_csv('http://opendata.cern.ch/record/545/files/Dimuon_DoubleMu.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "100000"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# This shows the number of rows in the data.\n",
    "\n",
    "len(muons)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA5sAAAJ7CAYAAABtWIUmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XmYZVV9L/zvT3DWCCgQAmjjK9GgN059UeIYiYBgBI0DGV7RS4I3ksHEvIIxV3BKMPqKMY6oRMwgDnHAYFQikKBGpHEWMLSKAmFoRgdEQdf9Y6+CQ1FVfbrZ1c0pPp/nOU/tvfZwfvuc06frW2vtvau1FgAAABjT7TZ3AQAAAKw8wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphE1iRqurIqmpVde4iy8/ty4/cxKVtVlW1pqretRHb/UFVtYn5x/fX70FTbr9df09WbcBztqr6g4n5U6vqAxtS9xL7fmZVPWeB9tGeY0xV9dKqurCqfrYx79/mUlV7VdULNtNz3+S9XKyWqnpXVa3ZiP3/Yv9Mb3VLawVYqYRNYCW7NskuVbV6srGq/meSVX05G+cLSfZI8s0p198uyREZXvdp7ZHk/RtW1tSemeQ5C7Q/P8mLl+k5N0r//L4syRuTPCrJKzZvRRtkrySbJWzm5u/l2LX8YobPtLAJsIgtN3cBAMvohxlC0YFJJnsuDkxycpKHb46iNlZV3T7Jz1prP93ctbTWvpfkc8ux76q6c2vtR621Zdn/UlprZ23q55zCA/rPN/XXnSncSt9LgNsUPZvASnd8kmdWVSVJ//nM3n4zVfWYqvr3qrqmqi6vqrdX1d0nlu9QVcdW1beq6kdV9V9V9cqqusO8/by4qtZW1bVVdUlVfbyqfr4ve04fInq3educV1WvnZg/tao+UFWHVNU3M/TE/kJf9qCqOrGqvt8f75/b/8T2D6qqz/Qazq6qp0zzglXVHavqjVV1VVVdUVVHJ7n9vHVuNoy2qg6uqrP663JZfx0f2IfOfrWvdkrfrs3bz95VdUJV/SBDD97NhtFOPM8h/bX6UX8NdlyqrsnXsk+/K8lvJHncXC1zw6kXGkZbVU+oqtMn3ss3T753E8/5+P4+/KB/Pp4/xWu9RR+K+d2q+nFVfb2qfmti+buS/H2fvXrueZbY36Kfi6q6a1X9sKoOXWC7M6rqHybm711Vx/f3/5qq+kRV3X9i+apeyzOr6m1VdXVVXVBVL6uq2/V1jkzywiT3mXid39WXPbD/m7ii13T2QnVNPN9xVfXJifn79/19cKLt4b1t1z4/+Z4vWsvE9k+sqq/0ej5dVQ9cop7HJ/lon/123995VXWn/j5Ovod/1Zc/ZaLtb6vqMxPz9+rHeHl/vU+teSMyFqhh7j04sKr+rqq+19+D3+nLX1RV/11V66rq1XPvS192s6HDE/t78kTbXarqDVV1cf/8n1FVe83bbu576rdq+M77XlX9a1XtNG+9Rb8TgZVL2ARWug8m2T7Jo/v8Y5Js29tvoqoeleTfklyc5OkZhtztm+TvJla7V5Irkvxpkn2SvCbJc5P87cR+np3kz5O8LsneSX4/ydokd92I+h/Vtz8sya9nCBz3S/KZJHdK8jsZhoM+MMlHq24I1XdO8okkd0vyW0lemeT1Se49xXMeleR3MwzX/O0k98nwi/qiquqxSd6aIRg9Kcn/SvLZJPdIclHfT5IcmmF47B7zdvHOJF9O8pQ+vZg9kvxhhtf/4CS/nOTDUxzTpFckOSXJFydqecdCK/bA8fEkl2UIqEdkeD0XOq/z7f0Ynprk1CRvqqrd11PLy5O8JMkxGY79M0n+sap+c6LWV/bpJ/Rav7BIrUt+LlprP0zyLxn+2DK53X2TrE7/A0xVbZPk00nun+R/9/XvmuTf+udq0l8n+UGGfy//kOSlfToZXtN/yvDvae51nhsC/NEkP+11PiXDv58b/qizgNOS7FFVW/T5x2b448ujJ9Z5bJJLWmsLnae9VC3J8O/iNUleleQ3Mwz7fu/cv6cFfCHJn/Xpp/X9PbW1dm2SMzJ8z0zWde0CbadNzH84w3fFnyV5Vobfz07p7+n6vDrDv7Hf6Ps8rqr+/yS7Z/h3+PokL8q8931Kb8/w/faqDJ/r85OcWFWPnrfeI5L8QYbviUOSPCzDZzrJ6N+JwCxprXl4eHisuEeSI5Nc1qc/kmEIYpK8OcmH+/RlSY6c2Oa0JKfM288TkrQkD1rkebbMED6uTXKH3vbGJP+8RG3P6fu827z285K8dmL+1CQ/SrL9vPX+Psk35p6vt+2a4Zf3/fr885Ncl2SniXUe1Z/3XUvUds/+nIdNtN0uyTnDfxk3tD1+8nXJ8EvymUvs90F9/cfPa5/bz9ELbNOS/MG81+O6JPde4Jj2Waiuedt+YGL+A0lOXeA55693fJJzk2wx0fbM/hx7zHvOl0+sc/sk65IctcRrsk2God5HzGv/WJJvrO/zssD+pvlcPLXP/8LEOi/O8AeU2/f5VyS5PMk2E+tsneTqJIf2+VW9pnfPq+FLSY6fmH9tkvPmrXOvvu3/2IB/z/fv26zu8+/O8MeN65I8oLd9MMn7l3gvb1ZLb39XkuuT7DrRdkB/vgcsUdOT+zqr5rX/VZKv9ek7Jflxhu+Ez/W2rea9J/v0/TxuYh937Z+fty3x/HPvwd9NtP1cf03mf2Y/n+S98455zSL7e3Kf/6UkP0ty0MQ6t0vytSSfmPc6X51k64m2F/R93bnPL/md6OHhsXIfejaB24Ljkzy9qu6YodflZkNoq+ouGXon3ldVW849MvTwXJd+fmcNXlB9uGhf9o9J7pgbew2/lGTfPqRw94nemI1xZmvtknltv5bkQ0l+NlHntzOE1bmhd7v3bS+Y26i19pkkl67n+f5Hhl+QPzKx3c8m5xfxpSQPraqjq+qxNW9Y8RROnHK9L7TWvjtR29wxra8HcWPtnuRD7abnyf5zhnAyv3fnhmGerbW5X/h3yuIelOQuuflFkN6b5BeratsNrHWaz8W/ZuiJfMbEds/KcIzXTeznpCTfm9jP95OcObGfOZ+cN39Wlj7mZAi25yd5a1U9q6q2W9+Btda+keF9nusdfGw/li9MtD06N+0t3BDntZv2iM6d77m+Y1nIfyTZrfcQPzLD6/2WJA/r3zNzn5u5YbS7J7m0tfbvcztoN/ZCz/+MLeRTE9t9L0NI/fd5n9m1SXacv+F6/M8klYnPZ/8ueP8CdZ3RWrtyYn7u9Zt7zjG/E4EZImwCtwUnZBhO+qoMPQYfXWCdrZNskaHn87qJx48z9FLt3Nd7QYYekg8l2T/DL4pz55rdqf88NsOQsWcmOT3JJTWc17kxv2DND5rJ0DN02Lw6r0ty34k6fz4LB8v1hc25c6jmr7fkdq21f8sw3O6xGXo6LquqN1XVtMPkFjrOhSx2TDtMuf2G2iHzauu/xF+eoWdy0lXz5n+SGz8Ti+078/c/MT9//+uz3s9FG4Z5fiRDwEwN52E+ODf9A8y9+vL5+/nV3Pj5mrOhxzwXWPbKMKT12CQXV9VpVfXQ9RzfaUkeU1U7Z/jDzqcn2n4pw/D4jQ2bCx1Hsp5jWcRnM/TqPTpDEP5MhvB1dYbw+ZgMPZ9zz7lDFv5cX5LpPgML1b7B78sCdkjyg9baNQvUdZf+x7ulakiW5zsRmCGuRguseK21H1bVvyT5kwzD7H64wGpXZfgF8cgMwxjn++/+8xkZhua9ZG5BVe027/l+luToJEf3X4x/O0PQvSDD0L+5W67M7/3beqHyF2i7IkPYXeg8w8v6z4tz41VMJ62vF+niifWu2IDt0lo7LsP5YttmOI/t6Aw9Yoevb9ssfJwLWaiO7TKcs5Ys/dpelg130fzn7L8g3zM3fX02xlzN22UIr3O27z83dP/TfC6Soef0o1V17wyhcl2GqzNP7ueELHyLle9vYE0Laq2dk+Q3arjC8mMynHd4YlXt1P/9LOS0DOe3PjbJWa21y6vqtAznJH4myfcynDO7WbXWrq6qr2Q4rodkGHLaqurTvW3++Zo3+4x12+eWf8YWc23W//1zUZK7VdVd5gXO7ZNc01r78bRPNsV3IrBC6dkEbivekqFHc8FfbHoA/VyS+7fW1izwmAubd87Q2znpt7OI1tr5rbWjMgxjmwulc0Nbf2luvap6RIbzrabxqQwXfjlzgTrP6+uckeThk1eE7BdAWl9o/GqGX0T3n9judpPz69NaW9dae1uGX6jnjvmW9BRNelgPSXO1zR3T53vTQq/tzrl58J62p+f0JE+d1wPztAx/rP30hpV+M19Lck1uOqQ1GXp//qu1tm4D9zfN5yIZhr5e1Z/nWRn+ePLTBfbz9QX2840NrGnJ17m1dl1r7eQMF47ZIUvfs/I/MvReHtKnk+Ezdp8M501/domgut5aNsJSn+n/yI0XdPqPiba9MwzJnwybpyfZrl9kK8kNw/r3yy3/jC3mgiSrqmqy9r3mrXNGhj8CzV3wae5q3k+/JXUt8p0IrFB6NoHbhNbaqRmGdy7lRUk+VVU/y3ABme9nGK63X5KXtNb+K8O5bH9UVacn+WaGoHmTK0ZW1dsy9Eh8LsPQuV/NcKGWw/oqn09yYZI3VNX/yTBU7kUZemamcWTfx4lVdWyGXqsdkzwxw8V/Ts1wBd2/6OscmSEkvyLr6d3rvUXHJHlZVV2f5OtJfi/DMORFVdXL+nGc2p/joUkelxt7Nb+b4cJDB1XV1Umua62tWWBX67OuH9MRGX7Jf3WG8zg/3uu/oIZbOryiqq7J8EfVP8/Ne4jOSbJ/VR2Q4Rfv/574g8KkV2a4au2Hq+otGc7he3WG3qr/3Ij6b9Bau6KqXp/kL/prvSZDkN03wxVRN9SRWf/nIq2162q4ZcifZgh482/R8roMV4k9uar+NsNndfsM7+enW2vv2YCazkmyfVU9J0O4vizDH1Vem6GH9VsZetQOS/Ll1tpSPXlfzvBv5LEZ/ng09xqe1dtessS2C9YyL4RvqLng/byqOj5Db9/cLX5OS/JHGc7X/MJE2+smptOP4RNV9dkMV789PEMv959l+Df7mltQ31I+nOFKyO+o4RYwD81w5dobtNbOrqr3JHljDbd/+maG74IHZLia7NSm+E4EVig9mwBda+3TGX5p3TbDlT0/miEEnp8bz6N7eZL3ZAgh78nQu/FH83b1n30/f5dhSO5Tk/xea+3D/Xl+0tvmQu0LM/zydmWm0EPvIzP0ih2T4UIpL8vQ47q2r3NNhl6UH2Y4H++I/jzfmeIpXpThHKuX9mP879z4S/JizsjQS/HWDLdc+f0M4edvej3XZvhF9eFJ/r2vvzE+m+RNGYZOvjNDaDhg3jq/mSHc/kOSv8zwns3vkXtzhh6+Y3sthyz0ZK21r2e4lct2Ga52Ove+P32h9TfCSzNcvfT3M1wQ5rFJfqe1tuB9YJcyzediwvEZguZ/Z955jq21y/p+zskw9PGTGW5xco8kX9nAst6X4cqnf53hdT4yw1DtSzKEw3/N8F6cneEWKEsd388yvP/Jjb2Fmah/fb1tC9Wy0Vpr38kQCp+WYRjv5LngczX9Z2vt+j79xQzh89sL/GHjgAx/yHp9hgvwVJIntNbmv2+jaK19LUO43CPDkOnHZTjner7fS3Jchs/pRzL0Ij+5f1duiCW/E4GVq1qb9jQZAAAAmI6eTQAAAEYnbAIAADA6YRMAAIDRCZsAAACMTtgEAABgdMImAAAAoxM2AQAAGJ2wCQAAwOiETQAAAEYnbAIAADA6YRMAAIDRCZsAAACMTtgEAABgdMImAAAAoxM2AQAAGJ2wCQAAwOiETQAAAEYnbAIAADA6YRMAAIDRCZsAAACMTtgEAABgdMImAAAAoxM2AQAAGN2Wm7uAWXOve92rrVq1anOXAQAAsFmceeaZl7XWtl3fesLmBlq1alXWrFmzucsAAADYLKrqO9OsZxgtAAAAoxM2AQAAGJ2wCQAAwOiETQAAAEYnbAIAADA6YRMAAIDRCZsAAACMTtgEAABgdMImAAAAoxM2AQAAGJ2wCQAAwOiETQAAAEYnbAIAADA6YRMAAIDRCZsAAACMbibDZlX9SVV9vaq+VlXvqao7VdUuVXV6Va2tqvdW1R36unfs82v78lUT+3lxb/9GVe29uY4HAABgpZm5sFlVOyb5oySrW2sPSrJFkgOTvDrJ0a21+yW5MsnBfZODk1zZ24/u66WqduvbPTDJPkneXFVbbMpjAQAAWKlmLmx2Wya5c1VtmeQuSS5K8oQkH+jLj0tyQJ/ev8+nL9+zqqq3H99a+3Fr7dtJ1ibZfRPVDwAAsKLNXNhsrV2Y5LVJvpshZF6d5MwkV7XWru+rXZBkxz69Y5Lz+7bX9/XvOdm+wDYAAADcAjMXNqtq6wy9krsk+YUkd80wDHY5n/OQqlpTVWvWrVu3nE8FAACwIsxc2Ezya0m+3Vpb11q7LskHkzwqyVZ9WG2S7JTkwj59YZKdk6Qvv0eSyyfbF9jmJlprx7TWVrfWVm+77bZjHw8AAMCKM4th87tJHllVd+nnXu6Z5KwkpyR5el/noCQf6dMn9Pn05Se31lpvP7BfrXaXJLsm+fwmOgYAAIAVbcv1r3Lr0lo7vao+kOQLSa5P8sUkxyQ5McnxVfXK3vbOvsk7k/x9Va1NckWGK9Cmtfb1qnpfhqB6fZJDW2s/3aQHAwAAsELV0MnHtFavXt3WrFmzucsY1arDT7xh+ryj9tuMlQAAALd2VXVma231+tabxWG0AAAA3MoJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMLqZC5tVdf+q+tLE43tV9YKq2qaqTqqqc/vPrfv6VVVvqKq1VfWVqnrYxL4O6uufW1UHbb6jAgAAWFlmLmy21r7RWntIa+0hSR6e5JokH0pyeJJPtdZ2TfKpPp8kT0qya38ckuQtSVJV2yQ5Iskjkuye5Ii5gAoAAMAtM3Nhc549k3yztfadJPsnOa63H5fkgD69f5J3t8HnkmxVVTsk2TvJSa21K1prVyY5Kck+m7Z8AACAlWnWw+aBSd7Tp7dvrV3Upy9Osn2f3jHJ+RPbXNDbFmsHAADgFprZsFlVd0jylCTvn7+stdaStBGf65CqWlNVa9atWzfWbgEAAFasmQ2bGc7F/EJr7ZI+f0kfHpv+89LefmGSnSe226m3LdZ+M621Y1prq1trq7fddtsRDwEAAGBlmuWw+Zu5cQhtkpyQZO6Ksgcl+chE+7P7VWkfmeTqPtz2E0n2qqqt+4WB9uptAAAA3EJbbu4CNkZV3TXJE5M8b6L5qCTvq6qDk3wnyTN7+8eS7JtkbYYr1z43SVprV1TVK5Kc0dd7eWvtik1QPgAAwIo3k2GztfbDJPec13Z5hqvTzl+3JTl0kf0cm+TY5agRAADgtmyWh9ECAABwKyVsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARjeTYbOqtqqqD1TVOVV1dlXtUVXbVNVJVXVu/7l1X7eq6g1VtbaqvlJVD5vYz0F9/XOr6qDNd0QAAAAry0yGzSR/k+TjrbUHJHlwkrOTHJ7kU621XZN8qs8nyZOS7NofhyR5S5JU1TZJjkjyiCS7JzliLqACAABwy8xc2KyqeyR5bJJ3Jklr7SettauS7J/kuL7acUkO6NP7J3l3G3wuyVZVtUOSvZOc1Fq7orV2ZZKTkuyzCQ8FAABgxZq5sJlklyTrkvxdVX2xqt5RVXdNsn1r7aK+zsVJtu/TOyY5f2L7C3rbYu0AAADcQrMYNrdM8rAkb2mtPTTJD3PjkNkkSWutJWljPWFVHVJVa6pqzbp168baLQAAwIo1i2HzgiQXtNZO7/MfyBA+L+nDY9N/XtqXX5hk54ntd+pti7XfTGvtmNba6tba6m233Xa0AwEAAFipZi5sttYuTnJ+Vd2/N+2Z5KwkJySZu6LsQUk+0qdPSPLsflXaRya5ug+3/USSvapq635hoL16GwAAALfQlpu7gI30h0n+sarukORbSZ6bITi/r6oOTvKdJM/s634syb5J1ia5pq+b1toVVfWKJGf09V7eWrti0x0CAADAyjWTYbO19qUkqxdYtOcC67Ykhy6yn2OTHDtudQAAAMzcMFoAAABu/YRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARjdV2Kyq7apql4n5qqpDqur1VfXry1ceAAAAs2jans13JfmTifmXJ3lzkn2SfKiqnjNuWQAAAMyyacPmw5KcnCRVdbsk/zvJn7fWHpDkVUlesDzlAQAAMIumDZv3SHJ5n354km2S/GOfPznJ/UauCwAAgBk2bdi8IMlufXq/JOe01i7s8/dIcu3YhQEAADC7tpxyvWOT/HVV/VqGsPniiWWPTHL22IUBAAAwu6bq2Wyt/VWSP0xycf/5honF2yR5x/ilLa6qzquqr1bVl6pqTW/bpqpOqqpz+8+te3tV1Ruqam1VfaWqHjaxn4P6+udW1UGb8hgAAABWsql6Nqvq3kne01p79wKL/zDJDqNWNZ1fba1dNjF/eJJPtdaOqqrD+/xhSZ6UZNf+eESStyR5RFVtk+SIJKuTtCRnVtUJrbUrN+VBAAAArETTnrP57SQPXWTZL/flm9v+SY7r08clOWCi/d1t8LkkW1XVDkn2TnJSa+2KHjBPynArFwAAAG6hacNmLbHsTkl+PEItG6Il+WRVnVlVh/S27VtrF/Xpi5Ns36d3THL+xLYX9LbF2gEAALiFFh1GW1W/nOQhE037VtUD5q12pyTPTPJfy1DbUh7dWruwqrZLclJVnTO5sLXWqqqN9WQ90B6SJPe+973H2i0AAMCKtdQ5m0/NcE5jMvQkvnSR9b6d5HljFrU+c7ddaa1dWlUfSrJ7kkuqaofW2kV9mOylffULk+w8sflOve3CJI+f137qIs93TJJjkmT16tWjhVgAAICVaqlhtH+Z5O5Jfi7DMNon9PnJxx1ba/9Pa+3flrvQOVV116q6+9x0kr2SfC3JCUnmrih7UJKP9OkTkjy7X5X2kUmu7sNtP5Fkr6raul+5dq/eBgAAwC20aM9ma+26JNf12WnP7dwUtk/yoapKhvr/qbX28ao6I8n7qurgJN/JMLw3ST6WZN8ka5Nck+S5SdJau6KqXpHkjL7ey1trV2y6wwAAAFi5prr1yZyq+sUMw03vNH9Za+1jYxW1lNbat5I8eIH2y5PsuUB7S3LoIvs6NsmxY9cIAABwWzftfTZ3S3J8kgdm4SvTtiRbjFgXAAAAM2zans23JbljkqclOSvJT5atIgAAAGbetGHzoUkObK39y3IWAwAAwMow7YV/vpkFztMEAACAhUwbNl+Y5M+r6r7LWQwAAAArw7TDaP8qyY5Jzqmq85JcNX+F1truI9YFAADADJs2bH6tPwAAAGC9pgqbrbXnLnchAAAArBzTnrOZJKnBzlX1K1V11+UqCgAAgNk2ddisqucnuTDJd5KcluT+vf2DVfWC5SkPAACAWTRV2Kyq/y/J65K8PckTktTE4lOTPGv0ygAAAJhZ014g6NAkL22t/XVVbTFv2TeS/OK4ZQEAADDLph1G+/NJzlxk2c+S3GmccgAAAFgJpg2ba5M8bpFlj01y1jjlAAAAsBJMO4z29UneXFU/SfKB3rZdVR2c5E+T/N5yFAcAAMBsmvY+m++oqq2TvDTJy3rzx5Jck+TI1to/LVN9AAAAzKBpezbTWntNVb01yR5J7pXkiiT/2Vq7ermKAwAAYDZNFTar6r6ttW+11r6f5JPLXBMAAAAzbuoLBFXV56vqT6pqp2WtCAAAgJk3bdj89SRnJzkiyXlVdVpHBGBiAAAgAElEQVRVHVpV2y9faQAAAMyqqcJma+3E1tpBSbZL8vQk5yc5KskFVfWpqvrdZawRAACAGTNtz2aSpLX2k9bah1trv5UheB6U5AFJ3rYcxQEAADCbpr4a7Zyqul2SJyR5VpKnJtk6yWdHrgsAAIAZNnXPZlU9rqrenOSiDFekfXCSv0xyn9baY5apPgAAAGbQtLc+uSjDsNmvJnl9kve21r61nIUBAAAwu6YdRvvWDAHznOUsBgAAgJVhqrDZWnvZ3HRVVZIdklzaWrt+uQoDAABgdm3IOZv7VtXpSa5N8t0kv9zbj6mq31mm+gAAAJhBU4XNqnp2khOSnJPkkHnbnZvk4PFLAwAAYFZN27P5kiSvaa0dlOQf5i37epLdRq0KAACAmTZt2LxPkpMWWXZtkp8bpxwAAABWgmnD5vlJHrrIstVJ1o5TDgAAACvBtGHznUmO6BcCunNvq6raM8mLkrx9OYoDAABgNk17n81XJ9k5yXFJftrbPptkiyRva629YRlqAwAAYEZNe5/NluTQqnpdkj2T3CvJFUlObq391zLWBwAAwAyatmczSdJa+2aSby5TLQAAAKwQ056zCQAAAFMTNgEAABidsAkAAMDoFg2bVXXvqrr9piwGAACAlWGpns1vJ3loklTVyVX1gE1TEgAAALNuqbD5oyR36dOPT/Jzy14NAAAAK8JStz75YpK/qaqT+vwfVtVFi6zbWmuHjVsaAAAAs2qpsPl7SV6TZP8kLcmeSX68yLotibAJAABAkiXCZmvtnCS/niRV9bMkB7TWPr+pCgMAAGB2LdWzOWmXJIsNoQUAAICbmCpstta+U1VbVtWzkjw6yTZJrkhyWpIPttauX8YaAQAAmDFThc2q2i7JJ5P8cpLzklySZI8khyb5clXt1Vpbt1xFAgAAMFuWuvXJpNcluWeSR7bW7tta26O1dt8kj+jtr1uuAgEAAJg904bNfZMcNv8CQa21M5K8OMl+YxcGAADA7Jo2bN4xyfcXWfb9JHcYpxwAAABWgmnD5ueSHFZVd51s7POH9eUAAACQZPpbn7wwySlJzq+qT2a4QNB2SfZOUkkevyzVAQAAMJOm6tlsrX0pya5JjkmybZInZgibb02ya2vty8tWIQAAADNn2p7NtNYuS3L4MtYCAADACjHtOZsAAAAwtZkNm1W1RVV9sar+pc/vUlWnV9XaqnpvVd2ht9+xz6/ty1dN7OPFvf0bVbX35jkSAACAlWdmw2aSP05y9sT8q5Mc3Vq7X5Irkxzc2w9OcmVvP7qvl6raLcmBSR6YZJ8kb66qLTZR7QAAACvaTIbNqtopyX5J3tHnK8kTknygr3JckgP69P59Pn35nn39/ZMc31r7cWvt20nWJtl90xwBAADAyrbesNmHob6kqh68KQqa0uuTvCjJz/r8PZNc1Vq7vs9fkGTHPr1jkvOTpC+/uq9/Q/sC2wAAAHALrDdsttZ+nOQlSbZa/nLWr6qenOTS1tqZm/A5D6mqNVW1Zt26dZvqaQEAAGbWtMNoT0/ysOUsZAM8KslTquq8JMdnGD77N0m2qqq5W7nslOTCPn1hkp2TpC+/R5LLJ9sX2OYmWmvHtNZWt9ZWb7vttuMeDQAAwAo0bdh8UZLnV9UfVNV9q+quVXWXycdyFjmptfbi1tpOrbVVGS7wc3Jr7beTnJLk6X21g5J8pE+f0OfTl5/cWmu9/cA+THiXJLsm+fwmOgwAAIAVbcv1r5Jk6NlMkjdk6EVcyOa+kuthSY6vqlcm+WKSd/b2dyb5+6pam+SKDAE1rbWvV9X7kpyV5Pokh7bWfrrpywYAAFh5pg2b/ytJW85CNkZr7dQkp/bpb2WBq8m21q5N8oxFtn9VklctX4UAAAC3TVOFzdbau5a5DgAAAFaQaXs2kyRVtVuSh2e4sM6xrbWLq+p+SS5prX1/OQoEAABg9kwVNqvqbkmOzXCBnev6dh9PcnGSv0zy3SR/tkw1AgAAMGOmvRrt65L8SpI9k9w9SU0s+1iSfUauCwAAgBk27TDapyX549baKVU1/6qz30lyn3HLAgAAYJZN27N55ySXL7Ls7kncMgQAAIAbTBs2z0jy7EWWPT3JZ8cpBwAAgJVg2mG0/yfJSVX1b0nen+Gem/tW1Z9kCJuPXab6AAAAmEFT9Wy21k7LcHGgOyZ5Y4YLBL0syX2T/Fpr7YxlqxAAAICZM/V9Nltrn0nymKq6c5Ktk1zVWrtm2SoDAABgZk17zuakazPca/NHI9cCAADACjF12KyqfavqsxnC5sVJrq2qz1bVfstWHQAAADNpqrBZVc9L8tEkP0jyx0me0X/+IMkJfTkAAAAkmf6czT9P8rbW2vPntb+1qt6a5CVJ3jZqZQAAAMysaYfR3jPJhxZZ9s9JthmnHAAAAFaCacPmKUket8iyxyX5j3HKAQAAYCVYdBhtVe02MfuGJO+oqnsm+XCSS5Nsl+SpSZ6U5HeXs0gAAABmy1LnbH4tSZuYryTP64/W5+d8PMkWo1cHAADATFoqbP7qJqsCAACAFWXRsNla+/dNWQgAAAArx7S3PrlBVW2Z5A7z21tr14xSEQAAADNvqqvRVtU9qurNVXVRkmuTfH+BBwAAACSZvmfzXRlucfL2JGuT/GS5CgIAAGD2TRs290zyvNbae5azGAAAAFaGqYbRJvluEudkAgAAMJVpw+aLkvxFVd17OYsBAABgZZhqGG1r7WNV9WtJ1lbVeUmuWmCd3UeuDQAAgBk1VdisqtcmeUGSM+ICQQAAAKzHtBcI+t0kL2mt/dVyFgMAAMDKMO05m9ckOXM5CwEAAGDlmDZs/k2SQ6qqlrMYAAAAVoZph9HeK8kjknyjqk7NzS8Q1Fprh41ZGAAAALNr2rD59CTXJ7l9kicusLwlETYBAABIMv2tT3ZZ7kIAAABYOaY9ZxMAAACmNu19Np+/vnVaa2++5eUAAACwEkx7zuYbl1jW+k9hkw2y6vATb5g+76j9NmMlAADA2KYaRttau938R5Jtkvxmki8n2W05iwQAAGC2TNuzeTOttauSvLeq7pHkbUkeP1ZRAAAAzLYxLhD07SSrR9gPAAAAK8QtCptVtUOSF2YInAAAAJBk+qvRrsuNFwKac4ckd09ybZKnjVwXAAAAM2zaczbflJuHzWuTXJDk4621y0etCgAAgJk2VdhsrR25zHWwDOZuLeK2IgAAwKa20Vej5dZp8t6Vs8i9NwEAYGVYNGxW1ckbsJ/WWttzhHoAAABYAZbq2ZzmPMwdkvxKbn4+JwAAALdhi4bN1tozFltWVfdOcliSJye5LMnR45cGAADArNqgczar6n5JXpzkd5Jc2qff1lr70TLUBgAAwIya9j6bD0zykiTPSHJ+kj9Ocmxr7SfLWBsAAAAz6nZLLayqh1fVB5N8JcnDkvxukl1ba28VNAEAAFjMUlej/dckeyX5apIDW2vv32RVAQAAMNOWGka7d/+5U5I3VdWbltpRa2270aoCAABgpi0VNl+2yaoAAABgRVnq1ifCJgAAABtlyQsEAQAAwMaYubBZVXeqqs9X1Zer6utV9bLevktVnV5Va6vqvVV1h95+xz6/ti9fNbGvF/f2b1TV3gs/IwAAABtq5sJmkh8neUJr7cFJHpJkn6p6ZJJXJzm6tXa/JFcmObivf3CSK3v70X29VNVuSQ5M8sAk+yR5c1VtsUmPBAAAYIWaubDZBj/os7fvj5bkCUk+0NuPS3JAn96/z6cv37Oqqrcf31r7cWvt20nWJtl9ExwCAADAijdzYTNJqmqLqvpSkkuTnJTkm0muaq1d31e5IMmOfXrHJOcnSV9+dZJ7TrYvsA0AAAC3wEyGzdbaT1trD8lwD9DdkzxgOZ+vqg6pqjVVtWbdunXL+VQAAAArwlL32bzVa61dVVWnJNkjyVZVtWXvvdwpyYV9tQuT7JzkgqraMsk9klw+0T5ncpv5z3NMkmOSZPXq1W05juW2btXhJ27uEgAAgBHNXM9mVW1bVVv16TsneWKSs5OckuTpfbWDknykT5/Q59OXn9xaa739wH612l2S7Jrk85vmKAAAAFa2WezZ3CHJcf3KsbdL8r7W2r9U1VlJjq+qVyb5YpJ39vXfmeTvq2ptkisyXIE2rbWvV9X7kpyV5Pokh7bWfrqJjwUAAGBFmrmw2Vr7SpKHLtD+rSxwNdnW2rVJnrHIvl6V5FVj1wgAAHBbN3PDaAEAALj1EzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxu5u6zyYZbdfiJN0yfd9R+m7ESAADgtkLPJgAAAKMTNgEAABidsAkAAMDohE0AAABGJ2wCAAAwOmETAACA0QmbAAAAjE7YBAAAYHTCJgAAAKMTNgEAABidsAkAAMDohE0AAABGJ2wCAAAwOmETAACA0QmbAAAAjE7YvI1ZdfiJWXX4iZu7DAAAYIUTNgEAABjdlpu7AFjMZA/seUftt+jyhZYBAACbl55NAAAARidsAgAAMDphEwAAgNE5Z5OZ4kq6AAAwG/RsAgAAMDphkxXFfUQBAODWQdgEAABgdMImAAAAoxM2AQAAGJ2r0TITnIcJAACzRdhk5gmiAABw62MYLQAAAKMTNgEAABidsAkAAMDohE0AAABGJ2wCAAAwOmETAACA0bn1CSvS5O1Qzjtqv81YCQAA3Dbp2QQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARidsAgAAMDq3PrmNmrw1CAAAwNj0bLLirTr8ROEaAAA2MT2b3GZMBs7zjtpvM1YCAAArn55NAAAARidsAgAAMDphEwAAgNEJmwAAAIxu5sJmVe1cVadU1VlV9fWq+uPevk1VnVRV5/afW/f2qqo3VNXaqvpKVT1sYl8H9fXPraqDNtcxAQAArDQzFzaTXJ/kha213ZI8MsmhVbVbksOTfKq1tmuST/X5JHlSkl3745Akb0mGcJrkiCSPSLJ7kiPmAioAAAC3zMyFzdbaRa21L/Tp7yc5O8mOSfZPclxf7bgkB/Tp/ZO8uw0+l2Srqtohyd5JTmqtXdFauzLJSUn22YSHAgAAsGLNXNicVFWrkjw0yelJtm+tXdQXXZxk+z69Y5LzJza7oLct1g4AAMAtNLNhs6ruluSfk7ygtfa9yWWttZakjfhch1TVmqpas27durF2CwAAsGLNZNisqttnCJr/2Fr7YG++pA+PTf95aW+/MMnOE5vv1NsWa7+Z1toxrbXVrbXV22677XgHAgAAsELNXNisqkryziRnt9ZeN7HohCRzV5Q9KMlHJtqf3a9K+8gkV/fhtp9IsldVbd0vDLRXbwMAAOAW2nJzF7ARHpXk/03y1ar6Um/78yRHJXlfVR2c5DtJntmXfSzJvknWJrkmyXOTpLV2RVW9IskZfb2Xt9au2DSHAAAAsLLNXNhsrX06SS2yeM8F1m9JDl1kX8cmOXa86gAAAEhmcBgtjGHV4Sdm1eEnbu4yAABgxRI2AQAAGJ2wCQAAwOiETQAAAEYnbAIAADA6YRMAAIDRCZsAAACMTtgEAABgdMImAAAAoxM2AQAAGN2Wm7sA2JxWHX7iDdPnHbXfZqwEAABWFj2bAAAAjE7YBAAAYHTCJnSrDj/xJsNqAQCAjSdsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARrfl5i4Abm1WHX7iDdPnHbXfZqwEAABml55NAAAARidsAgAAMDrDaGEJhtQCAMDG0bMJAADA6IRNAAAARidsAgAAMDphEwAAgNEJmwAAAIxO2AQAAGB0wiYAAACjEzYBAAAYnbAJAADA6IRNAAAARrfl5i6Acaw6/MTNXQIAAMAN9GwCAAAwOmETAACA0QmbAAAAjE7YBAAAYHTCJkxp1eEnuhATAABMydVo2eQENgAAWPn0bAIAADA6PZuwgSZ7Zs87ar/NWAkAANx66dkEAABgdMIm/N/27j3csrq+7/j7w6AY0qQDYpAyEyGBakCtF7xQmtRExUHUodUCxiRotDxtSRN9bOMYWzX65ClGGy+1qCjjiKWgRROhgxeCREXFMIgiVyWIYUZuMgEvJFzkmz/WOrBms8+cc2bWnLXP2e/X85znnPVba//29+zf8ztnf/fvsiRJkiT1zmRTkiRJktQ7k01JkiRJUu9MNiVJkiRJvTPZlHbCAes2et9QSZIkaQyTTUmSJElS77zPptQD770pSZIkbcuRTUmSJElS70w2JUmSJEm9M9mUJEmSJPVuSSabSdYnuTXJFZ2yvZOcn+Q77fe92vIkeU+S65JcnuQpncec0F7/nSQnDPG7SJIkSdJytCSTTWADsGakbB1wQVUdDFzQHgMcBRzcfp0IvA+a5BR4E/AM4OnAm2YSVGlneDsUSZIkaYkmm1X1RWDrSPFa4CPtzx8BjumUn16Ni4GVSfYDngecX1Vbq+rvgPN5aAIrSZIkSdoBSzLZnMW+VXVT+/PNwL7tz/sDN3au29yWzVY+1RyVkyRJktSH5ZRsPqCqCqi+6ktyYpJNSTbddtttfVUrSZIkScvW7kMH0KNbkuxXVTe102Rvbcu3AKs7161qy7YAzxop/6txFVfVqcCpAIcddlhvSayWt+4I8Q0nHz1gJJIkSdLiW04jm+cAMzvKngB8qlP+O+2utM8E7myn234WODLJXu3GQEe2ZZIkSZKknbQkRzaTnEkzKrlPks00u8qeDHw8ySuB7wHHtpefBzwfuA64C3gFQFVtTfJW4JL2urdU1eimQ5IkSZKkHbAkk82qeuksp5495toCTpqlnvXA+h5Dk8aamVLrdFpJkiRNi+U0jVaSJEmSNCFMNiVJkiRJvTPZlCRJkiT1zmRTkiRJktQ7k01JkiRJUu+W5G600lI1syttlzvUSpIkaTlyZFOSJEmS1DuTTUmSJElS70w2JUmSJEm9M9mUJEmSJPXODYK0KMZtjCNJkiRp+XJkU5ogB6zbaGIuSZKkZcGRTWlgJpeSJElajhzZlCRJkiT1zmRTkiRJktQ7k01JkiRJUu9csylNoO46zhtOPnrASCRJkqQd48imJEmSJKl3jmxKE85RTkmSJC1FjmxKkiRJknrnyKZ2Ke8h2a+Z19MRTkmSJE06k02N5dRNSZIkSTvDabSSJEmSpN6ZbEqSJEmSeuc0Ws1p3LpLp9ZOHqc+S5IkaZKYbGqnjEtw3BRIkiRJksmmdsi4hNIkc/H4WkuSJGnSuWZTWoYOWLfRhFSSJEmDMtmUJEmSJPXOZFOSJEmS1DvXbEqSpKni7t2StDhMNqVlzDdUkiRJGorTaCVJkiRJvXNkU5oSM6OcjnBKmlbbu22XfxslqX8mm9KUcWqtpGnibaAkaThOo5UkSZIk9c6RTWmKjfvE39FOSZIk9cFkU5IkLTtOn5Wk4TmNVpIkSZLUO0c2JY3lRkKSlpqdGc30b54k9c9kU9I2tndrgC7fjEmSJGl7nEYrSZIkSeqdI5uSJEkdTqmVpH6YbEraIb4ZkyRJ0vaYbEraadvblMNEVNKu5m1OJGkymWxK2qUcAZUkSZpObhAkadEcsG6jIxCSJElTwpFNSYvOW6lIkiQtf45sSpo4joBKmhT+PZKkHefIpiRJWnJMACVp8plsSpoI4944Ot1WkiRp6TLZlLSkzDWaYTIqLW9DjWi6s7YkLZzJpqRlxTeE0vIzaVNmZ+Lxb4wkbZ8bBElattzYQ5IkaTiObEpa9nYk4ZwZsXCkVBrOpH9Y5N8HSdo+k01JGmO+GxZ1+WZT2jGTnlRKknbM1CebSdYA7wZWAB+qqpMHDknSErW9N8zdRHTcei/XgElLm31Ykh4qVTV0DINJsgL4NvBcYDNwCfDSqrpqtsccdthhtWnTpkWKcP78VFhSl294NST/J23L/ihpuUlyaVUdNtd10z6y+XTguqq6HiDJWcBaYNZkU5KWgr7f7I97s7yQ9WpL8Z6pu3I93s7UPddj+xphM2Hsz3xfy0nvE5K0UNM+svkSYE1Vvao9/m3gGVX1e7M9xpFNSZI0yfpKWp0aLGk2jmz2KMmJwInt4Y+TXDtkPLPYB/jB0EFoG7bJZLJdJo9tMplsl8kzrzbJ2/p90r7rW2bsJ5PJdtn1HjOfi6Y92dwCrO4cr2rLtlFVpwKnLlZQOyLJpvl8uqDFY5tMJttl8tgmk8l2mTy2yeSxTSaT7TI5dhs6gIFdAhyc5MAkDweOB84ZOCZJkiRJWvKmemSzqu5L8nvAZ2lufbK+qq4cOCxJkiRJWvKmOtkEqKrzgPOGjqMHEz3Nd0rZJpPJdpk8tslksl0mj20yeWyTyWS7TIip3o1WkiRJkrRrTPuaTUmSJEnSLmCyucQlWZPk2iTXJVk3dDzTKsnqJBcmuSrJlUn+oC3fO8n5Sb7Tft9r6FinTZIVSS5L8v/b4wOTfK3tMx9rNwfTIkqyMsnZSa5JcnWSw+0rw0rymvZv1xVJzkzyCPvK4kuyPsmtSa7olI3tG2m8p22fy5M8ZbjIl69Z2uTt7d+vy5P8eZKVnXOvb9vk2iTPGybq5W9cu3TOvTZJJdmnPbavDMhkcwlLsgL438BRwCHAS5McMmxUU+s+4LVVdQjwTOCkti3WARdU1cHABe2xFtcfAFd3jt8GvLOqDgL+DnjlIFFNt3cDn6mqxwH/gqZ97CsDSbI/8PvAYVX1eJoN847HvjKEDcCakbLZ+sZRwMHt14nA+xYpxmmzgYe2yfnA46vqicC3gdcDtP/3jwcObR9zSvteTf3bwEPbhSSrgSOBv+0U21cGZLK5tD0duK6qrq+qe4CzgLUDxzSVquqmqvp6+/OPaN4870/THh9pL/sIcMwwEU6nJKuAo4EPtccBfgM4u73ENllkSf4p8GvAaQBVdU9V3YF9ZWi7Az+TZHdgT+Am7CuLrqq+CGwdKZ6tb6wFTq/GxcDKJPstTqTTY1ybVNXnquq+9vBimvu0Q9MmZ1XV3VX1XeA6mvdq6tksfQXgncAfAt1NaewrAzLZXNr2B27sHG9uyzSgJAcATwa+BuxbVTe1p24G9h0orGn1Lpp/Ove3x48E7ui8SbDPLL4DgduAD7fTmz+U5GexrwymqrYA76AZCbgJuBO4FPvKpJitb/geYDL8LvDp9mfbZEBJ1gJbquqbI6dslwGZbEo9SvJPgE8Ar66qH3bPVbP1s9s/L5IkLwBurapLh45F29gdeArwvqp6MvATRqbM2lcWV7sGcC3NBwH/DPhZxkxP0/DsG5MlyRtoltGcMXQs0y7JnsAfAW8cOhZty2RzadsCrO4cr2rLNIAkD6NJNM+oqk+2xbfMTNVov986VHxT6AjgRUluoJli/hs0awVXtlMFwT4zhM3A5qr6Wnt8Nk3yaV8ZznOA71bVbVV1L/BJmv5jX5kMs/UN3wMMKMnLgRcAL6sH7yNomwznl2k+MPtm+39/FfD1JI/GdhmUyebSdglwcLtj4MNpFqWfM3BMU6ldC3gacHVV/Vnn1DnACe3PJwCfWuzYplVVvb6qVlXVATR94/NV9TLgQuAl7WW2ySKrqpuBG5M8ti16NnAV9pUh/S3wzCR7tn/LZtrEvjIZZusb5wC/0+60+Uzgzs50W+1CSdbQLNF4UVXd1Tl1DnB8kj2SHEizIc1fDxHjtKmqb1XVL1TVAe3//c3AU9r/OfaVAeXBD2O0FCV5Ps26tBXA+qr6k4FDmkpJ/hXwJeBbPLg+8I9o1m1+HPhF4HvAsVU1bkG7dqEkzwL+S1W9IMkv0Yx07g1cBvxWVd09ZHzTJsmTaDZtejhwPfAKmg8/7SsDSfLHwHE0UwIvA15Fs6bJvrKIkpwJPAvYB7gFeBPwF4zpG+0HA++lmfJ8F/CKqto0RNzL2Sxt8npgD+D29rKLq+o/tNe/gWYd5300S2o+PVqndt64dqmq0zrnb6DZYfsH9pVhmWxKkiRJknrnNFpJkiRJUu9MNiVJkiRJvTPZlCRJkiT1zmRTkiRJktQ7k01JkiRJUu9MNiVJkiRJvTPZlCRJkiT1zmRTkiQta0lWJ7kgydVJrkzyp+2N3iVJu5DJpiRJWu7uA15XVb8CPBl4BvBvhw1JkpY/k01JUq+SvDnJD4aOY1SSDUk27aK6j03y8l1R99DS+EaSE8acOybJ55LcnuSeJFuSnJ1kzQKf49wk39rO+fcmuSPJHp2yY5PcPJ8Ryqq6qao2tT/fA1wOrB7zHKctJG5J0vaZbEqSpsVbgZfvorqP3YV1D+1YYG/g/3YLk7wT+ASwBXgV8BxgHfAzwKeT/PICnuNM4PFJDhk9kWQF8BLgk1V1d+fU0cB5VVULeB6SPBI4BvjsyKl3AC9LctBC6pMkzW73oQOQJGlXaROVFVV1T1X9zdDxLFG/D3y0qu6dKUiyFng18Iqq2jBy/UeTvBD4+wU8x6eAu4CXAv995NyvA/vSJKQzz78bcBTwHxfwHLQjo2cD76qqq7vnquqGJBe1db52IfVKksZzZFOStEvNTF9N8twklyf5SZKLkhzann95OwVz5cjjDk1SSZ7THh+e5JwkN7V1fCPJy2Z5rmOSXAn8A836vIdMo11gfWNjn7kGeDHwr9t4K8mb53gtjk5yVZK7kmxMsneSg5Jc2D7HpiRPHHnsfOI9NMlnkmxtr7k6yUnzPT8m3oOAf0mToHW9GrhkTKIJQFWdW1XfH6nrV5N8of2db0/ywSQ/117/E+Bc4Lgx1R0P3Ap8vlP2NGAlcH6n/hVJXtO20z+0r9MZM9Ns2w8ezgAuq6r/Ocuv/Ama0U3fH0lSDxzZlCQthl8E3g78Cc2I1zuAjyV5AvAXwAeAfwN8uPOY44BbgAvb48cAXwbeT5NEHgF8OMn9VXVm53EHAH8KvAW4GfjuLDHNt75ZY2+ncL61vWYl8J/ax2ye47V4C/DfgAAVSJ0AAAXgSURBVD2B/wWc2sb9wTb2/wGcleTQzjTR+cR7LnA18FvA3cBjgZ/vPPdc50c9G/gJ8M2ZgiS7A4e3r8O8JDkC+Euatn4J8EjgZGCv9hiakcvjkjy1qi5tH/cwmo18zqiqn3aqPBr4UlX9sL1uN5qE+Fk0r93XadZkHtV5/T4A/Ijtj1p+hWYU9Qnd31mStGNMNiVJi2Fv4Iiq+g48kBz8OfDYqromyWdoksvRZPPsmSSjqs6aOdGOVn0RWAX8ezpTLGkSmedU1Te2F9AC6ps1duCaqvqbJFuB3arq4nm+FofPTOttRzD/K3BCVZ3eiWcj8Dia5HDOeJPsAxwIrK2qmc12Lug8ZrvnZ/FU4Oqqur9T9khgD+DG7oVtTCs6RT/tJHonA1+pquM6128BLkjy+Kq6Avg0cAfNSOal7WXPo0lIu+0BTbJ5Ruf4NcBzgaeNTI/9cPtcRwCvBK4ALmsHO9dX1XtG6r0S+CnwdEw2JWmnOU1EkrQYbphJ1lpXtd9Xtd8/Bjw7zeYtJHkS8M/bctqyvZK8J8n3gHvbrxPb67q2zJVoLrC+uWJfqBtG1o9e137//Jiy/RcQ71aaBPD9SY5L8gsjzzvX+XEeDcy2s/Doxjyv7cR1L3BSG/eeNCOhH0+y+8wXcFF73VPhgV1iPwkcOzP1leYDh+8BX+28DvvR3L5kY3u8G/A64L2j6zAfCLTqy1WVqnpCVT2p/RpNNKmq+2gS3kdv5zWRJM2TyaYkaTHcMXJ8T/v9Ee33c2gSjxe3x8fRTEW9qPOYDW3524Ejadbtre/UMeOWecY03/rmin2hZqvvjjFl3efYwHbibUcfj6SZOrweuDnJl5I8eT7nZ/EImum2Xbe3ZaPJ9kfbmJ42Ur4XzYjnKWybjN4NPIxtb0FyJs0048OTPAJYC5w1suPs84Hrq+ra9viJwKNoRpv7cDc73raSpA6n0UqSBldVP06ykSaZOpXmdhv/bybJaBOPFwAnVdX7Zx43y0Yuc94KY4H1DW6+8VbVNcCL27WOvwq8DdiYZFVV3T/X+TFPvZWRUb6qui/JV2kS1zd2ym+hTfSz7a0v76BpkzcD5415ju5GQhe2dRwP7Af8HOOn0G7sHO/Xfr95TN07YiXN7y1J2kkT+U9VkjSVzqLZ0fWFwC+1xzP2oPmf9cAoW7uT6Yt28Ln6ru8edu1o2ILirap7q+rzwJ/RJGMrF3K+41qadZ6j3gU8I8lvzxV4u9PsxTTrczeN+fp+59qfAh8H/h3wmzTrRbubEz2cZm1mN9mcSTIfN1csc0nyKJpNm769s3VJkhzZlCRNjvNo7rX4AeC7VfXXMyeq6s4klwBvTPJD4H5gHXAn299Nday+6wOuAdYmOYZm+u/3R2/9sTPmE2+70dA7aNa5Xk8zffV1wDerautc52d56i+3z/moqrqtE8+nkrwL2JDk12l2uf0BzeZBR7aX/bhTzx/SbAZ0P82usT+imS57NPCGquomd2cC/5lmd+I3jcTza0CAL3TKvkWzsc/7kryRZqT0EGB1Vb1ult9rNofRjMJ+ZYGPkySN4cimJGkiVNXf06zd3I/OxkAdv0mTJJ0OvJvmnoin78RT9lnfKcDnaNZCXkKzcU/f5or3ZpopqG+g2dn1FJqdbF80z/Pj/BXNlNI1oyeq6jU0ty1ZDZxGs8HRKTS3Dnl+9x6cVXURTaL4KJq1nefSJKA3MrLGtqq+CtxAk1SOm0L7l1V1d+f6+4AXApcB72zrflV7vFBrgC9U1e078FhJ0ohsu+ZekiTpQUneDRxUVUdPQCzfBt5eVR/cBXWvoNn5dl1V/Z++65ekaWSyKUmSZpVkFc0axieNTHddVpIcD7wV+JV2tFSStJOcRitJkmZVVZuB3+XBXV+XqwCvNNGUpP44silJkiRJ6p0jm5IkSZKk3plsSpIkSZJ6Z7IpSZIkSeqdyaYkSZIkqXcmm5IkSZKk3plsSpIkSZJ6Z7IpSZIkSeqdyaYkSZIkqXf/CI6sSZoMfptaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1080x720 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure(figsize=(15, 10))\n",
    "\n",
    "plt.hist(muons.M, 300, range = (0,150))\n",
    "\n",
    "plt.xlabel('Invariant mass (GeV/$c^2$)', fontsize = 15)\n",
    "plt.ylabel('Number of events \\n', fontsize = 15)\n",
    "plt.title('Measured distribution of events with two muons \\n', fontsize = 15)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA44AAAJ7CAYAAACs38+KAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XmYbFddL/zvj4QZhAROYiSEAy8IAlcZjpFcRhNAJAgBkcEpeJF4BVEEXxLEC0FFg/ACImMY41UJg0DAIIKBSAANSRBkSJDpAMGMYJhCIIH1/rF3h0qne50653R3VfX5fJ6nnq491q+qdlfXt9fea1VrLQAAALCaa8y6AAAAAOab4AgAAECX4AgAAECX4AgAAECX4AgAAECX4AgAAECX4AjMvao6tqpaVX1mleWfGZcfu8GlzVRVnVlVr9uF7X6nqtrE9H3G1++OU26/3/iebN2Jx2xV9TsT06dW1Zt3pu7Ovh9RVY9ZYf6aPcZaqqpnVNVXquoHu/L+zUpV3b+qnjSjx77Ke7laLVX1uqo6cxf2/+PjMX3j3a0VYLMSHIFFcVmSW1bVtsmZVfXTSbaOy9k1H0lySJLPTbn+fkmemeF1n9YhSd60c2VN7RFJHrPC/Mcnedo6PeYuGY/fZyV5cZK7J/mT2Va0U+6fZCbBMVd/L9e6lh/PcEwLjgCr2HvWBQBM6dsZAs6jkky2KDwqyXuT3HUWRe2qqrpmkh+01r4/61paa99I8m/rse+qum5r7TuttXXZf09r7VMb/ZhTuN348yXj684U5vS9BNijaHEEFsmJSR5RVZUk489HjPOvpqruWVX/UlWXVtVXq+qVVXXDieUHVNVrqurzVfWdqvrPqvrTqrrWsv08rao+W1WXVdUFVfWuqvrRcdljxtMwb7Bsm+1V9byJ6VOr6s1VdVRVfS5DC+mPjcvuWFUnV9U3x9ublvY/sf0dq+qDYw1nV9WDp3nBquraVfXiqrqkqr5WVS9Ics1l61ztVNWqemxVfWp8XS4eX8c7jKenfnxc7X3jdm3Zfn6uqt5eVd/K0LJ2tVNVJx7nqPG1+s74GtysV9fkaznef12SX0xy76Valk5ZXulU1ao6tKpOn3gvXzr53k085n3G9+Fb4/Hx+Cle673G0x2/VFXfrapPVtUvTyx/XZL/O05+felxOvtb9bioqutX1ber6gkrbHdGVf3NxPRBVXXi+P5fWlX/VFW3nVi+dazlEVX1iqr6elWdW1XPqqprjOscm+QpSW4x8Tq/blx2h/F34mtjTWevVNfE451QVe+emL7tuL+3TMy76zjvNuP05Hu+ai0T29+vqv5jrOcDVXWHTj33SfKOcfIL4/62V9V1xvdx8j3883H5gyfm/VVVfXBi+qbjc/zq+HqfWsvOlFihhqX34FFV9dqq+sb4HvzquPypVfVfVXVRVT1n6X0Zl13t9NyJ/T1oYt71qupFVXX+ePyfUVX3X7bd0ufUL9fwmfeNqvrHqjpw2XqrfiYCm5fgCCyStyTZP8k9xul7Jtkyzr+Kqrp7kn9Ocn6Sh2c4re2BSV47sdpNk3wtyZOTPCDJc5P8RpK/mtjPryf5wyTPT/JzSX47yWeTXH8X6r/7uP3RSX4hQ3i4dZIPJrlOkl/NcMrlHZK8o+rKgHzdJP+U5AZJfjnJnyZ5YZKDpnjM45L8ZoZTIn8lyS0yfOleVVXdK8nLM4Scn0/yv5J8KMmNkpw37idJnpDhFNRDlu3i1Uk+luTB4/3VHJLkiRle/8cm+ckkb5viOU36kyTvS/LvE7W8aqUVx/DwriQXZwibz8zweq50HeQrx+fw0CSnJnlJVR28g1r+OMnTkxyf4bl/MMnfVtWjJ2r90/H+oWOtH1ml1u5x0Vr7dpJ/yPCPk8ntbpVkW8Z/plTVvkk+kOS2Sf73uP71k/zzeFxN+osk38rw+/I3SZ4x3k+G1/TvMvw+Lb3OS6fZviPJ98c6H5zh9+fKf9Cs4LQkh1TVXuP0vTL8I+UeE+vcK8kFrbWVrmvu1ZIMvxfPTfLsJI/OcGr1G5Z+n1bwkSR/MN5/2Li/h7bWLktyRobPmcm6Llth3mkT02/L8FnxB0kemeG71vvG93RHnpPhd+wXx32eUFX/X5KDM/wevjDJU7PsfZ/SKzN8vj07w3H95SQnV9U9lq33M0l+J8PnxFFJ7pLhmE6y5p+JwCJprbm5ubnN9S3JsUkuHu+flOE0vyR5aZK3jfcvTnLsxDanJXnfsv0cmqQlueMqj7N3hiBxWZJrjfNenOTvO7U9ZtznDZbN357keRPTpyb5TpL9l633f5N8eunxxnm3yfBF/PBx+vFJLk9y4MQ6dx8f93Wd2m4yPubRE/OukeSc4eP/ynn3mXxdMnzhPauz3zuO699n2fyl/bxghW1akt9Z9npcnuSgFZ7TA1aqa9m2b56YfnOSU1d4zOXrnZjkM0n2mpj3iPExDln2mH88sc41k1yU5LjOa7JvhtOpn7ls/juTfHpHx8sK+5vmuHjoOP1jE+s8LcM/Q645Tv9Jkq8m2XdinX2SfD3JE8bprWNNf72sho8mOXFi+nlJti9b56bjtv9jJ36fbztus22c/usM/6i4PMntxnlvSfKmznt5tVrG+a9LckWS20zMO2J8vNt1anrQuM7WZfP/PMknxvvXSfLdDJ8J/zbOu/Gy9+QB437uPbGP64/Hzys6j7/0Hrx2Yt6PjK/J8mP2w0nesOw5n7nK/h40Tv9Ekh8kOXJinWsk+USSf1r2On89yT4T85407uu643T3M9HNzW3z3rQ4AovmxCQPr6prZ2gNudppqlV1vQytBm+sqr2XbhlaXi7PeD1kDZ5U4ymZ47K/TXLt/LA176NJHjietnfwRCvJrjirtXbBsnn3TfLWJD+YqPMLGYLn0ultB4/bnru0UWvtg0ku3MHj/Y8MX3ZPmtjuB5PTq/hokjtX1Quq6l617NTdKZw85Xofaa19aaK2pee0o5a9XXVwkre2q15X+vcZgsbyVpcrT6VsrS19eT8wq7tjkuvl6h0AvSHJj1fVlp2sdZrj4h8ztBD+0sR2j8zwHC+f2M97knxjYj/fTHLWxH6WvHvZ9KfSf87JEFK/nOTlVfXIqtpvR0+stfbpDO/zUqvdvcbn8pGJeffIVVvxdsb2dtWWyqXrI3f0XFby/iS3H1tu75bh9X5ZkruMnzNLx83SqaoHJ7mwtfYvSztoP2wdXn6MreSUie2+kSFw/suyY/azSW62fMMd+OkklYnjc/wseNMKdZ3RWvvvieml12/pMdfyMxFYIIIjsGjenuGUzWdn+E/+O1ZYZ58ke2Vokbx84vbdDK1HNx/Xe1KGlou3JnlIhi99S9dmXWf8+ZoMp2U9IsnpSS6o4TrIXfmytDw0JkOLzdHL6rw8ya0m6vzRrBwSdxQcl645Wr5ed7vW2j9nOKXtXhlaIC6uqpdU1bSnoq30PFey2nM6YMrtd9YBWVbb+IX8qxlaDCddsmz6e/nhMbHavrN8/xPTy/e/Izs8LtpwKuVJGcJiarhu8ady1X+m3HRcvnw/P5sfHl9LdvY5L4WP+2c4bfQ1Sc6vqtOq6s47eH6nJblnVd08wz9pPjAx7ycynIK+q8FxpeeR7OC5rOJDGVrb7pEh1H4wQ5D6eoYgec8MLZJLj3lAVj6uL8h0x8BKte/0+7KCA5J8q7V26Qp1XW/8R1yvhmR9PhOBBaJXVWChtNa+XVX/kOT3M5zK9u0VVrskw5e9YzOcKrjcf40/fynD6W9PX1pQVbdf9ng/SPKCJC8Yv+T+SobQem6G0+uWhgFZ3iq3z0rlrzDvaxmC60rX5V08/jw/P+yNc9KOWnfOn1jvazuxXVprJ2S4vmpLhuu+XpChpeqYHW2blZ/nSlaqY78M13gl/df24uy885Y/5vhl9ya56uuzK5Zq3i9DEF2y//hzZ/c/zXGRDC2a76iqgzIExIsy9DI8uZ+3Z+VhP765kzWtqLV2TpJfrKGn4HtmuE7v5Ko6cPz9WclpGa4HvVeST7XWvlpVp2W4hu+DSb6R4RrTmWqtfb2q/iPD87pThtM6W1V9YJy3/PrGqx1jo/2z+8fYai7Ljj9/zktyg6q63rLwuH+SS1tr3532wab4TAQ2KS2OwCJ6WYaWxhW/pIxh8t+S3La1duYKt6XgeN0MrZCTfiWraK19ubV2XIZTxZYC5tLpoz+xtF5V/UyG65OmcUqGTk/OWqHO7eM6ZyS562TPhmPnPzsKgB/P8KXyIRPbXWNyekdaaxe11l6R4cvx0nPenRacSXcZA89SbUvP6cPjrJVe25vn6iF62haY05M8dFnLyMMy/BP1AztX+tV8Ismlueppo8nQKvOfrbWLdnJ/0xwXyXB66SXj4zwywz9Cvr/Cfj65wn4+vZM1dV/n1trlrbX3Zug05YD0x0R8f4ZWxaPG+8lwjN0iw3XGH+qEzh3Wsgt6x/T788POjN4/Me/nMpz2PhkcT0+y39jBVJIrT50/PLt/jK3m3CRbq2qy9vsvW+eMDP/QWersaKlX6ofvTl2rfCYCm5QWR2DhtNZOzXAKZc9Tk5xSVT/I0HnKNzOcEnd4kqe31v4zw7Vfv1tVpyf5XIbQeJWeD6vqFRlaCv4tw+lpP5uhk5Kjx1U+nOQrSV5UVf8nw+loT83QYjKNY8d9nFxVr8nQmnSzJPfL0PHNqRl6gv2jcZ1jMwTeP8kOWt3GVpzjkzyrqq5I8skkj8twqu+qqupZ4/M4dXyMOye5d37Y2vilDJ3uHFlVX09yeWvtzBV2tSMXjc/pmRm+sD8nw3WP7xrrP7eGYQb+pKouzfDPzj/M1VtuzknykKo6IsOX6P+a+OfApD/N0Pvq26rqZRmueXtOhlakf92F+q/UWvtaVb0wyR+Nr/WZGULpAzP07Lmzjs2Oj4u01i6vYRiLJ2cIa8uHDXl+ht5O31tVf5XhWN0/w/v5gdba63eipnOS7F9Vj8kQlC/O8A+S52Vo+fx8hpauo5N8rLXWa2H7WIbfkXtl+EfQ0mv4qXHe0zvbrljLskC9s5ZC9G9V1YkZWuGWhp05LcnvZri+8SMT854/cT/jc/inqvpQhl5cj8nQ+vwHGX5nn7sb9fW8LUOPvq+qYViSO2fogfVKrbWzq+r1SV5cw5BEn8vwWXC7DL2iTm2Kz0Rgk9LiCGxKrbUPZPgCuiVDD5XvyBDovpwfXnf2x0lenyFQvD5Dq8PvLtvVv477eW2G014fmuRxrbW3jY/zvXHeUkB9SoYvYv+dKYwB9m4ZWquOz9BJyLMytIR+dlzn0gytG9/OcP3aM8fH+eIUD/HUDNckPWN8jv+VH37hXc0ZGVoPXp5hGJDfzhBk/nKs57IMXzrvmuRfxvV3xYeSvCTD6YmvzhAAjli2zqMzBNW/SfJnGd6z5S1lL83Q8vaasZajVnqw1tonMwwvsl+GXjuX3veHr7T+LnhGhl44fztDZyj3SvKrrbUVxxntmea4mHBihtD4X1l2XWBr7eJxP+dkOL3w3RmG3bhRkv/YybLemKEHz7/I8Dofm+F06AsyBL1/zPBenJ1hWI7e8/tBhvc/+WErXibq31Er2Eq17LLW2hczBLyHZThVdvLa6aWa/rW1dsV4/98zBMkvrPBPiiMy/FPqhRk6n6kkh7bWlr9va6K19okMQfGQDKcl3zvDNcrLPS7JCRmO05MytO4+aPys3Bndz0Rg86rWpr0UBQAAgD2RFkcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC69p51AbN005vetG3dunXWZQAAAMzEWWeddXFrbcuO1tujg+PWrVtz5plnzroMAACAmaiqL06znlNVAQAA6BIcAQAA6BIcAQAA6BIcAQAA6BIcAQAA6BIcAQAA6FrI4FhVv19Vn6yqT1TV66vqOlV1y6o6vao+W1VvqKprzbpOAACAzWDhgmNV3SzJ7ybZ1lq7Y5K9kjwqyXOSvKC1dusk/53ksbOrEgAAYPNYuOA42jvJdatq7yTXS3JekkOTvHlcfkKSI2ZUGwAAwKaycMGxtfaVJM9L8qUMgfHrSc5Kcklr7YpxtXOT3Gw2FQIAAGwuCxccq2qfJA9JcsskP5bk+kkesBPbH1VVZ1bVmRdddNE6VQkAALB5LFxwTHLfJF9orV3UWrs8yVuS3D3JjcdTV5PkwCRfWWnj1trxrbVtrbVtW7Zs2ZiKAQAAFtgiBscvJblbVV2vqirJYUk+leR9SR4+rnNkkpNmVB8AAMCmsnDBsbV2eoZOcD6S5OMZnsPxSY5O8uSq+mySmyR59cyKBAAA2ET23vEq86e19swkz1w2+/NJDp5BOQAAAJvawrU4AgAAsLEERwAAALoERwAAALoERwAAALoERwAAALoERwAAALoERwAAALoERwAAALoERwBg3Ww95uRsPebkWZcBwG4SHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOgSHAEAAOjae9YFAAAs2XrMyVfe337c4TOsBIBJWhwBAADoEhwBAADoWrjgWFW3raqPTty+UVVPqqp9q+o9VfWZ8ec+s64VAABgM1i44Nha+3Rr7U6ttTsluWuSS5O8NckxSU5prd0mySnjNAAAALtp4YLjMocl+Vxr7YtJHpLkhHH+CUmOmFlVAAAAm8iiB8dHJXn9eH//1tp54/3zk+y/0gZVdVRVnVlVZ1500UUbUSMAAMBCW9jgWFXXSvLgJG9avqy11pK0lbZrrR3fWtvWWtu2ZcuWda4SAABg8S1scEzy80k+0lq7YJy+oKoOSJLx54UzqwwAAGATWeTg+Oj88DTVJHl7kiPH+0cmOWnDKwIAANiEFjI4VtX1k9wvyVsmZh+X5H5V9Zkk9x2nAQAA2E17z7qAXdFa+3aSmyyb99UMvawCAACwhhayxREAAICNs5AtjgDAxth6zMlX3t9+3OEzrASAWdLiCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAB7kK3HnHyVITY2alsAFpvgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQNfesy4AAFgsW485+cr72487fIaVALBRtDgCAADQJTgCAADQJTgCAADQJTgCAADQJTgCAADQJTgCAADQZTgOANjkJofPAIBdocURAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACALsERAACAroUMjlV146p6c1WdU1VnV9UhVbVvVb2nqj4z/txn1nUCAABsBgsZHJP8ZZJ3tdZul+Snkpyd5Jgkp7TWbpPklHEaAACA3bRwwbGqbpTkXklenSStte+11i5J8pAkJ4yrnZDkiNlUCAAAsLksXHBMcsskFyV5bVX9e1W9qqqun2T/1tp54zrnJ9l/ZhUCAABsIosYHPdOcpckL2ut3TnJt7PstNTWWkvSVtq4qo6qqjOr6syLLrpo3YsFAABYdIsYHM9Ncm5r7fRx+s0ZguQFVXVAkow/L1xp49ba8a21ba21bVu2bNmQggEAABbZwgXH1tr5Sb5cVbcdZx2W5FNJ3p7kyHHekUlOmkF5AAAAm87esy5gFz0xyd9W1bWSfD7Jb2QIwW+sqscm+WKSR8ywPgAAgE1jIYNja+2jSbatsOiwja4FAABgs1u4U1UBAADYWAvZ4ggA7J6tx5x85f3txx2+x9YAwHS0OAIAANAlOAIAANAlOAIAANAlOAIAANAlOAIAANAlOAIAANBlOA4AWCPrObzE0r4NWwHALGhxBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoGvvWRcAAMzW1mNOTpJsP+7wNdl2aR4Am4cWRwAAALoERwAAALoERwAAALoERwAAALoERwAAALoERwAAALoMxwEAXM1aD6kxub9dGfZjvfcHQJ8WRwAAALoERwAAALoERwAAALoERwAAALoERwAAALoERwAAALoMxwEAC8QwFADMghZHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAuvSqCgAbbEc9oy4t12sqAPNCiyMAAABdgiMAAABdgiMAAABdgiMAAABdgiMAAABdgiMAAABdhuMAYI8178NeTA7b0Vs+r/Xvrs3+/AAWiRZHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugzHAQBzYEdDb+zKfgxjAcBa0eIIAABA11QtjlW1X5Lrt9a+ME5XkscluX2SU1pr71i/ElesZ3uSbyb5fpIrWmvbqmrfJG9IsjXJ9iSPaK3990bWBQAAsBlN2+L4uiS/PzH9x0lemuQBSd5aVY9Z27Km8rOttTu11raN08dkCLG3SXLKOA0AAMBumjY43iXJe5Okqq6R5H8n+cPW2u2SPDvJk9anvJ3ykCQnjPdPSHLEDGsBAADYNKYNjjdK8tXx/l2T7Jvkb8fp9ya59RrXtSMtybur6qyqOmqct39r7bzx/vlJ9t/gmgAAADalaXtVPTfD9YynJTk8yTmtta+My26U5LJ1qK3nHq21r4zXXr6nqs6ZXNhaa1XVVtpwDJpHJclBBx20/pUCMFfWqvdSANiTTNvi+Jokf1FVb0ry1CTHTyy7W5Kz17qwnqXQ2lq7MMlbkxyc5IKqOiBJxp8XrrLt8a21ba21bVu2bNmokgEAABbWVMGxtfbnSZ6Y4RTQJyZ50cTifZO8au1LW1lVXb+qbrh0P8n9k3wiyduTHDmudmSSkzaqJgAAgM1s2uE4Dkry+tbaX6+w+IlJDljTqvr2z9CTazLU/3ettXdV1RlJ3lhVj03yxSSP2MCaAAAANq1pr3H8QpJDknx4hWU/Oc7fa62K6mmtfT7JT60w/6tJDtuIGgAAAPYk017jWJ1l10ny3TWoBQAAgDm0aotjVf1kkjtNzHpgVd1u2WrXyXBK6H+uQ20AAADMgd6pqg9N8szxfkvyjFXW+0KS31rLogBg0S0N+7H9uMOnWm/RbdTz2CyvF8Ci6Z2q+mdJbpjkRzKcqnroOD15u3Zr7f9prf3zehcKAADAbKza4thauzzJ5ePktNdCAgAAsMlM26tqkqSqfjzJgRmubbyK1to716ooAAAA5se04zjePsmJSe6QlXtYbdmg4TgAAADYWNO2OL4iybWTPCzJp5J8b90qAgAAYK5MGxzvnORRrbV/WM9iAAAAmD/TBsfPZYXrGgGA2VurISo2+1AXk89vR8OkAHBV0/aW+pQkf1hVt1rPYgAAAJg/07Y4/nmSmyU5p6q2J7lk+QqttYPXsC4AAADmxLTB8RPjDQAAgD3MVMGxtfYb610IAAAA82naaxyTJDW4eVX9z6q6/noVBQAAwPyY9lTVVNXjk/xRkh9N0pL8dJKPVNVbkry/tfbC9SkRADaHzd5r6bRWeh16r80sX7elx9YLK7Cnm6rFsar+3yTPT/LKJIcmqYnFpyZ55JpXBgAAwFyYtsXxCUme0Vr7i6raa9myTyf58bUtCwAAgHkx7TWOP5rkrFWW/SDJddamHAAAAObNtMHxs0nuvcqyeyX51NqUAwAAwLyZ9lTVFyZ5aVV9L8mbx3n7VdVjkzw5yePWozgAAABmb9pxHF9VVfskeUaSZ42z35nk0iTHttb+bp3qAwAAYMamHo6jtfbcqnp5kkOS3DTJ15L8a2vt6+tVHADA7jAECsDamCo4VtWtWmufb619M8m717kmAAAA5sjUneNU1Yer6ver6sB1rQgAAIC5Mm1w/IUkZyd5ZpLtVXVaVT2hqvZfv9IAAACYB1MFx9baya21I5Psl+ThSb6c5Lgk51bVKVX1m+tYIwAAADM0bYtjkqS19r3W2ttaa7+cIUQemeR2SV6xHsUBAAAwe1P3qrqkqq6R5NAkj0zy0CT7JPnQGtcFADM12Rvn9uMOX3UeAOwJpm5xrKp7V9VLk5yXoWfVn0ryZ0lu0Vq75zrVBwAAwIxNOxzHeRlOTf14khcmeUNr7fPrWRgAAADzYdpTVV+eISyes57FAAAAMH+mCo6ttWct3a+qSnJAkgtba1esV2EAAADMh525xvGBVXV6ksuSfCnJT47zj6+qX12n+gAAAJixqYJjVf16krcnOSfJUcu2+0ySx659aQAAAMyDaa9xfHqS57bWnlZVeyV57cSyTyb5gzWvDADYLZPDh7AyQ6wATGfaU1VvkeQ9qyy7LMmPrE05AAAAzJtpg+OXk9x5lWXbknx2bcoBAABg3kwbHF+d5JljJzjXHedVVR2W5KlJXrkexQEAADB7017j+JwkN09yQpLvj/M+lGSvJK9orb1oHWoDAABgDkw7jmNL8oSqen6Sw5LcNMnXkry3tfaf61gfAAAAMzZti2OSpLX2uSSfW6cwP1N1AAAgAElEQVRaAAAAmEM7FRwBYLObdggLQ10AsCeZtnMcAAAA9lCCIwAAAF2rBseqOqiqrrmRxQAAADB/ei2OX0hy5ySpqvdW1e02piQAAADmSS84fifJ9cb790nyI+teDQAAAHOn16vqvyf5y6p6zzj9xKo6b5V1W2vt6LUtDQAAgHnQC46PS/LcJA9J0pIcluS7q6zbkgiOAAAAm9CqwbG1dk6SX0iSqvpBkiNaax/eqMIAAACYD70Wx0m3TLLaaaoAAABsYlMFx9baF6tq76p6ZJJ7JNk3ydeSnJbkLa21K9axRgAAAGZoquBYVfsleXeSn0yyPckFSQ5J8oQkH6uq+7fWLlqvIgEAAJid3nAck56f5CZJ7tZau1Vr7ZDW2q2S/Mw4//nrVSAAAACzNe01jg9M8jvLO8dprZ1RVU9L8ldrXhkAbAJbjzl51iXsMSZf6+3HHT7DSgA2n2lbHK+d5JurLPtmkmutTTkAAADMm2mD478lObqqrj85c5w+elwOAADAJjTtqapPSfK+JF+uqndn6BxnvyQ/l6SS3Gddquuoqr2SnJnkK621B1XVLZOcmOGay7OS/Fpr7XsbXRcAAMBmM1WLY2vto0luk+T4JFuS3C9DcHx5ktu01j62bhWu7veSnD0x/ZwkL2it3TrJfyd57AxqAgAA2HSmbXFMa+3iJMesYy1Tq6oDkxye5NlJnlxVleTQJL88rnJCkmOTvGwmBQIAAGwi017jOG9emOSpSX4wTt8kySWttSvG6XOT3GwWhQEAAGw2U7c4zouqelCSC1trZ1XVfXZh+6OSHJUkBx100BpXB8AiMmQGAPQtYovj3ZM8uKq2Z+gM59Akf5nkxlW1FIQPTPKVlTZurR3fWtvWWtu2ZcuWjagXAABgoS1ccGytPa21dmBrbWuSRyV5b2vtVzL0+vrwcbUjk5w0oxIBAAA2lR0Gx6q6dlU9vap+aiMK2g1HZ+go57MZrnl89YzrAQAA2BR2eI1ja+27VfX0JB/YgHp2Smvt1CSnjvc/n+TgWdYDAACwGU17qurpSe6ynoUAAAAwn6btVfWpSf6uqi5P8s4kFyRpkyu01i5d49oAYC7odRWAPd20wfH08eeLMvRgupK9dr8cAAAA5s20wfF/ZVkLIwAAAHuGqYJja+1161wHAAAAc2raFsckSVXdPsldk9w8yWtaa+dX1a2TXNBa++Z6FAgAAMBsTRUcq+oGSV6T5OFJLh+3e1eS85P8WZIvJfmDdaoRAACAGZp2OI7nJ/mfSQ5LcsMkNbHsnUkesMZ1AQAAMCemPVX1YUl+r7X2vqpa3nvqF5PcYm3LAgDYdYZQAVhb07Y4XjfJV1dZdsMk31+bcgAAAJg30wbHM5L8+irLHp7kQ2tTDgAAAPNm2lNV/0+S91TVPyd5U4YxHR9YVb+fITjea53qAwAAYMamanFsrZ2WoWOcayd5cYbOcZ6V5FZJ7ttaO2PdKgQAAGCmph7HsbX2wST3rKrrJtknySWttUvXrTIAAADmwrTXOE66LMNYjt9Z41oAAACYQ1O3OFbVA5P8UZK7jttdUVVnJXl2a02f1wCsi6VhFbYfd/jV5i2fD+tl2mPOsQlsVlO1OFbVbyV5R5JvJfm9JL80/vxWkrePywEAANiEpm1x/MMkr2itPX7Z/JdX1cuTPD3JK9a0MgAAAObCtNc43iTJW1dZ9vdJ9l2bcgAAAJg30wbH9yW59yrL7p3k/WtTDgAAAPNm1VNVq+r2E5MvSvKqqrpJkrcluTDJfkkemuTnk/zmehYJAADA7PSucfxEkjYxXUl+a7y1cXrJu5LstebVAcAuWKknVlgrji9gT9QLjj+7YVUAAAAwt1YNjq21f9nIQgAAAJhP0w7HcaWq2jvJtZbPb61duiYVAQAAMFem6lW1qm5UVS+tqvOSXJbkmyvcAAAA2ISmbXF8XYZhN16Z5LNJvrdeBQEAADBfpg2OhyX5rdba69ezGAAAAObPtMHxS0lcwwjAhlga7gAAmA9TXeOY5KlJ/qiqDlrPYgAAAJg/U7U4ttbeWVX3TfLZqtqe5JIV1jl4jWsDAABgDkwVHKvqeUmelOSM6BwHAABgjzLtNY6/meTprbU/X89iAAAAmD/TXuN4aZKz1rMQAAAA5tO0LY5/meSoqnpPa62tZ0EAADtjV3rh1XMvwM6ZNjjeNMnPJPl0VZ2aq3eO01prR69lYQAAAMyHaYPjw5NckeSaSe63wvKWRHAEAADYhKYdjuOW610IAAAA82naznEAAADYQ007juPjd7ROa+2lu18OAAAA82baaxxf3Fm21Muq4AgAALAJTXWqamvtGstvSfZN8ugkH0ty+/UsEgAAgNmZtsXxalprlyR5Q1XdKMkrktxnrYoCAABgfqxF5zhfSLJtDfYDAADAHNqt4FhVByR5SobwCAAAwCY0ba+qF+WHneAsuVaSGya5LMnD1rguAAAA5sS01zi+JFcPjpclOTfJu1prX13TqgAAAJgbUwXH1tqx61wHAAAAc2qXe1UFgLW29ZiTZ10CALCCVYNjVb13J/bTWmuHrUE9AAAAzJlei+M01y0ekOR/5urXPwIAALBJrBocW2u/tNqyqjooydFJHpTk4iQvWPvSAAAAmAc7dY1jVd06ydOS/GqSC8f7r2itfWcdagMAAGAOTDuO4x2SPD3JLyX5cpLfS/Ka1tr31rE2AAAA5sA1egur6q5V9ZYk/5HkLkl+M8ltWmsvFxoBmAdbjzlZb6ysCccSwOp6var+Y5L7J/l4kke11t60YVUBAAAwN3qnqv7c+PPAJC+pqpf0dtRa22/NqgIAAGBu9ILjszasip1QVddJ8v4k185Q/5tba8+sqlsmOTHJTZKcleTXnE4LAACw+3rDccxlcEzy3SSHtta+VVXXTPKB8bTaJyd5QWvtxKp6eZLHJnnZLAsFAADYDLqd48yjNvjWOHnN8daSHJrkzeP8E5IcMYPyAAAANp2FC45JUlV7VdVHM4wl+Z4kn0tySWvtinGVc5PcbFb1AQAAbCZTjeM4b1pr309yp6q6cZK3JrndtNtW1VFJjkqSgw46aH0KBGDNGSYBAGZnIVscl7TWLknyviSHJLlxVS0F4QOTfGWVbY5vrW1rrW3bsmXLBlUKAACwuBYuOFbVlrGlMVV13ST3S3J2hgD58HG1I5OcNJsKAQAANpdFPFX1gCQnVNVeGYLvG1tr/1BVn0pyYlX9aZJ/T/LqWRYJAACwWSxccGyt/UeSO68w//NJDt74igAAADa3hTtVFQAAgI0lOAIAANC1cKeqAsBKVhquY3Le9uMO38hy2AM4voA9iRZHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAuvaedQEAsBG2HnPyrEtgE3N8AZudFkcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6DMcBADBho4bWWHqc7ccdviGPB7A7tDgCAADQJTgCAADQJTgCAADQJTgCAADQJTgCAADQpVdVADbMZG+VepIEgMWhxREAAIAuwREAAIAuwREAAIAuwREAAIAuwREAAIAuwREAAIAuw3EAAKyjXRmGZq2HrjEUDrC7tDgCAADQJTgCAADQJTgCAADQJTgCAADQJTgCAADQpVdVAHbKevb2CHu6aX8fltbTQyqwUbQ4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0GU4DoA9xFoPozG5T0MCwK4zJA2wCLQ4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0CU4AgAA0GU4DgB2264M9WEIAjY7xziwmWhxBAAAoGvhgmNV3byq3ldVn6qqT1bV743z962q91TVZ8af+8y6VgAAgM1g4YJjkiuSPKW1dvskd0vyhKq6fZJjkpzSWrtNklPGaQAAAHbTwgXH1tp5rbWPjPe/meTsJDdL8pAkJ4yrnZDkiNlUCAAAsLksXHCcVFVbk9w5yelJ9m+tnTcuOj/J/jMqCwAAYFNZ2F5Vq+oGSf4+yZNaa9+oqiuXtdZaVbVVtjsqyVFJctBBB21EqQC7ZFd6Kl2Lx5t8rJXm7cp+gPnk9xWY1kK2OFbVNTOExr9trb1lnH1BVR0wLj8gyYUrbdtaO761tq21tm3Lli0bUzAAAMACW7jgWEPT4quTnN1ae/7EorcnOXK8f2SSkza6NgAAgM1oEU9VvXuSX0vy8ar66DjvD5Mcl+SNVfXYJF9M8ogZ1QcAALCpLFxwbK19IEmtsviwjawFAABgT7Bwp6oCAACwsQRHAAAAuhbuVFUAFs/k0CIAwOLR4ggAAECX4AgAAECX4AgAAECX4AgAAECX4AgAAECX4AgAAECX4TgAFtxKQ11sP+7wNd3fLPcD7Jzd/d1b2n6tP0d2Z3/A7GlxBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEtwBAAAoEuvqgAAezi9IAM7osURAACALsERAACALsERAACALsERAACALsERAACALsERAACALsNxAOyBlrre337c4Tu1PrB7duV3aWd/XwHWgxZHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAugzHAeyRJrvE18U9sKh25bNsrYbX2Z1hQgwxAotHiyMAAABdgiMAAABdgiMAAABdgiMAAABdgiMAAABdelUFmJKeWIE91Z7SC6rPeVidFkcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6DMcBbHrr2b36rux7I7q1n6wL2Bzm9ffaEBawZ9DiCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJfgCAAAQJdeVYF1tQi97a3Uy+mO6u71brgez3lHNQLMO59ZsNi0OAIAANAlOAIAANAlOAIAANAlOAIAANAlOAIAANAlOAIAANBlOA6AGVmpa/qNHrJE9/jAaubt82GlYYmAjaPFEQAAgC7BEQAAgK6FDI5V9ZqqurCqPjExb9+qek9VfWb8uc8sawQAANgsFjI4Jnldkgcsm3dMklNaa7dJcso4DQAAwG5ayODYWnt/kq8tm/2QJCeM909IcsSGFgUAALBJLWRwXMX+rbXzxvvnJ9l/lsUAAABsFptyOI7WWquqttKyqjoqyVFJctBBB21oXcDOW+vu13ele/mN7JJ+R481b93jA+yOyc80w2zAfNtMLY4XVNUBSTL+vHCllVprx7fWtrXWtm3ZsmVDCwQAAFhEmyk4vj3JkeP9I5OcNMNaAAAANo2FDI5V9fok/5rktlV1blU9NslxSe5XVZ9Jct9xGgAAgN20kNc4ttYevcqiwza0EAAAgD3AQrY4AgAAsHEWssWRzWcz9qq2o+e0GZ/zjqx1D6kbte+1oDdUYL35nFkfe+Lfa1iJFkcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6BEcAAAC6DMexh5j3oQp2x57UTfbuvI872ra3fNrXeD3fi0V9n3WPD2wm8/qZtit/IzbzdyNYD1ocAQAA6BIcAQAA6BIcAQAA6BIcAQAA6BIcAQAA6BIcAQAA6DIcxwJaz+6jN0vX1Kt1y73Wz2/a/a119+U72t9GvH9rNTzGWgwxAsB82Ii/d7szRNTO1Lezz2Wl9Tf6+9RaD121O8OcrFUNzA8tjgAAAHQJjgAAAHQJjgAAAHQJjgAAAHQJjgAAAHTpVXWTmFVPqxvdc9Z6Pt6u9NK21na3B7i1rmFW+1mrbXf2GNnd566XV4CNtaNeV+fdrHqz1/Mpu0KLIwAAAF2CIwAAAF2CIwAAAF2CIwAAAF2CIwAAAF2CIwAAAF2G49jEdtRF9e4MVTDtcBXr2cXzRg8b0VtPt9a7b2fez0Xqah2A6czqb+lG/U3ZlaE3pt1mVn8XN+J75SzNariUeaXFEQAAgC7BEQAAgC7BEQAAgC7BEQAAgC7BEQAAgC7BEQAAgC7DcezB5nVIg5W6Pp7XWlcyD7WuZ/fdOxrmZVf3sR7m4b0AYPbW+u/BWg8lsSv7W+thyTZqyIl5fy9YnRZHAAAAugRHAAAAugRHAAAAugRHAAAAugRHAAAAuvSqusDWonfLtXzsHa23Uk9Xvf2s53PZ0b7noYextXj+szxGAGBPt1Z/h9e6B9VF2HZHvez3vleutM1afafb0f5Weq6bpbdXLY4AAAB0CY4AAAB0CY7w/7d371F3VOUdx78/EglidSVAxEjQglAQxIXcL4XWghBACAolQUQQKKsW2kptBRctUlmuyq0gVeRSIIKYgHiDJspdLipKKAESEyRAFMJNiICIJiQ8/WPvQyeHmXnPyXuZnMPvs9as953L3rPP88475zxn9uwxMzMzM7NaThzNzMzMzMyslhNHMzMzMzMzq+XE0czMzMzMzGr5cRw2Ykb6ERft+x3psmZmZtZfRuIxV0199hiO/a7K49uGou6hemxcXX1vxM+IvuJoZmZmZmZmtZw4mpmZmZmZWS0njmZmZmZmZlbLiaOZmZmZmZnVcuJoZmZmZmZmtTyq6mpodRilaTjbsDq8PjMzM1t9+bPC4I10DHvpbzbSo6Q29WSBoeYrjmZmZmZmZlbLiaOZmZmZmZnV6rvEUdIkSQ9KWijppKbbY2ZmZmZm1uv6KnGUNAr4KrAPsAVwqKQtmm2VmZmZmZlZb+urxBHYAVgYEY9ExDJgBjC54TaZmZmZmZn1tH5LHDcAHivMP56XmZmZmZmZ2SpSRDTdhiEj6WBgUkQck+cPB3aMiOML2xwLHJtnNwMeHPGGDmw94NmmG/EG5dg3x7FvluPfHMe+OY59cxz75jj2zVldY//uiBg/0Eb99hzHxcCGhfmJedlrIuIi4KKRbFS3JM2OiO2abscbkWPfHMe+WY5/cxz75jj2zXHsm+PYN6fXY99vXVXvBjaVtJGkNYGpwLUNt8nMzMzMzKyn9dUVx4hYLul44HpgFHBpRMxruFlmZmZmZmY9ra8SR4CImAXMarodg7Rad6Xtc459cxz7Zjn+zXHsm+PYN8exb45j35yejn1fDY5jZmZmZmZmQ6/f7nE0MzMzMzOzIebEcYRJOkHSPElzJU2XtFYezOdnkhZKuioP7FNW9nN5mwcl7T3Sbe91FbG/MsdzrqRLJb2pouwKSXPy5AGXulQR+2mSHi3EdeuKskdIeihPR4x023tdRezvKMT9CUnfqyjr434QJP1jjvs8SZ/Oy9aRdGM+nm+UNK6irI/7QaiI/ZmSFki6X9J3JY2tKLtI0gP5uJ89si3vfRWxP1XS4sL5ZN+KspPye/JCSSeNbMv7Q0X8ryrEfpGkORVlfex3IX9ufEbS3MKy0nO8kvPysX2/pG0q6tw2/w0W5u01Uq+nIxHhaYQmYAPgUeDNef5q4Mj8c2pedgHwqZKyWwD3AWOAjYCHgVFNv6ZemWpivy+gPE0vi33e/qWmX0OvTjWxnwYcPEDZdYBH8s9x+fdxTb+mXpmqYt+2zbeBT1SU93G/6rF/HzAXWJs0nsBNwCbAGcBJeZuTgNNLyvq4H57Y7wWMztucXhb7vG4RsF7Tr6MXp5rYnwr88wBlR+XPNhsDa+bPPFs0/Zp6aaqKf9s2ZwOnVJT3sd9dvHcHtgHmFpaVnuPz580f5M+bOwE/q6jz53m98vb7NP06i5OvOI680cCbJY0m/WM/CfwVcE1e/3XgwJJyk4EZEbE0Ih4FFgI7jEB7+0l77J+IiFmRkf5ZJzbawv71uth3WG5v4MaIWBIRvwVuBCYNUxv7VWXsJb2NdP4pveJog/Je0geDlyNiOXAb8FHSufzreZuq872P+8EpjX1E3JDnAe7C5/vhUHXcd2IHYGFEPBIRy4AZpP8X61xt/PPVq0NIX5TbIEXE7cCStsVV5/jJwOX5I+ddwFhJE4oF8/zbIuKu/Ln0csrfIxrjxHEERcRi4Czg16SE8QXgHuD5wpvZ46SrBO02AB4rzFdtZyXKYh8RN7TW5y6qhwM/rKhiLUmzJd0labX6J17dDRD7L+YuG+dIGlNS3Mf9IAx03JPekG6OiBcrqvBxv+rmArtJWlfS2qRvmzcE1o+IJ/M2TwHrl5T1cT84VbEvOor0bX6ZAG6QdI+kY4exnf2oLvbH5/P9pRVdtH3cD95Ax/5uwNMR8VBFeR/7g1d1ju/k+N4gL6/bplFOHEdQPlFOJnU1fSfwFvwt8ogoi72kjxc2OR+4PSLuqKji3RGxHfAx4FxJ7xnWBveRmth/Dtgc2J7UJe/ExhrZpzo47g+l/ptnH/erKCLmk7pD3kD6QmoOsKJtmyB9ULMhNFDsJZ0MLAeurKjizyNiG2Af4DhJuw9vi/tHTey/BrwH2Jr0JdbZTbWxn3Vw3hnonO9jfwj14zneiePI2hN4NCJ+ExGvAN8BdiVdrm49U3MisLik7GJW/taoajsrVxb7XQAkfR4YD/xTVeF85YaIeAT4EfCB4W5wHymNfUQ8mbtsLAUuo7zrtY/7wak77tcjxXxmVWEf94MTEZdExLYRsTvwW+CXwNOt7kn55zMlRX3cD1JF7JF0JPBh4LD8oa6sbOu4fwb4Lr4tpCtlsY+IpyNiRUS8ClyMz/fDpubYH03qtnpVTVkf+4NXdY7v5PhezMpd6Fe7/wEnjiPr18BOktbO/cz3AH4B3AocnLc5Avh+SdlrgamSxkjaCNiUdE+edaYs9vMlHUO6n+jQ/Ib2OpLGtbpR5g/bu5L+btaZqti3TqwidZmcW1L2emCv/DcYRxrc4voRanc/KI19Xncw8D8R8ceygj7uB0/S2/PPd5E+sH2TdC5vjZJadb73cT9IZbGXNAn4LHBARLxcUe4tkt7a+p0U+7Jzk1WoiH3xXq6PUB7Tu4FNlUaaXxOYSvp/sS5UnHcgfZG4ICIeryjnY39oVJ3jrwU+kUdX3Yl068iTxYJ5/kVJO+X37E9Q/h7RnKEaZcdTZxPw78AC0j/jFaRRUjcmJYELgW8BY/K2BwBfKJQ9mTTi2IOsZqMs9cJUEfvlOaZz8nRK3nY74L/z77sAD5BGeHsAOLrp19JrU0Xsb8nxnAt8A/iT9tjn+aPy/8ZC4JNNv5Zem8pin5f/CJjUtq2P+6GN/R2kZPs+YI+8bF3gZuAh0oiH67THPs/7uB/62C8k3WPUOt9fkJe/E5iVf984l7kPmAec3PRr6bWpIvZX5PPI/aQP0BPaY5/n9yVdIXvYsR+6+Ofl04C/bdvWx/7gYj2d1PX6FdL9iEfXnOMFfDUf2w8A2xXqmVP4fbv8fv0w8BVATb/O4qTcSDMzMzMzM7NS7qpqZmZmZmZmtZw4mpmZmZmZWS0njmZmZmZmZlbLiaOZmZmZmZnVcuJoZmZmZmZmtZw4mpmZmZmZWS0njmZmZmZmZlbLiaOZmZn1DEkbSrpZ0nxJ8ySdIUlNt8vMrN85cTQzM7Neshw4MSLeC3wA2BH4aLNNMjPrf04czcyskqRTJT3bdDvaSZomafYw1X2IpCOHo+6mKZkj6YiSdQdKukHSc5KWSVos6RpJk7rcx3WSHqhZ/xVJz0saU1h2iKSnOrlyGBFPRsTs/Psy4H5gw5J9XNJNu83MrJ4TRzMz60WnAUcOU92HDGPdTTsEWAf4ZnGhpHOAbwOLgWOAPYGTgDcDP5D0ni72MR14n6Qt2ldIGgUcDHwnIpYWVu0HzIqI6GI/SFoXOBC4vm3VWcBhkjbppj4zM6s2uukGmJmZdSInHaMiYllEPNx0e3rUPwBXRMQrrQWSJgOfBj4ZEdPatr9C0v7AH7rYx/eBl4FDgX9rW/dBYH1Sctna/xrAPsCnutgH+YrlNcC5ETG/uC4iFkm6M9f5mW7qNTOzcr7iaGZmHWt1EZX0IUn3S/q9pDslbZnXH5m7OY5tK7elpJC0Z57fWdK1kp7MdcyRdFjFvg6UNA/4I+l+ttd1Ve2yvtK2t7YBDgL+Irc3JJ06QCz2k/QLSS9LmilpHUmbSLo172O2pPe3le2kvVtK+qGkJXmb+ZKO63R9SXs3AXYhJVtFnwbuLkkaAYiI6yLiiba6dpN0W37Nz0m6WNJb8/a/B64DppRUNxV4BrilsGx7YCxwY6H+UZJOyH+nP+Y4Xdnqypq/RLgSuDcizq54yd8mXXX0Zx0zsyHgK45mZtatdwFnAl8kXYk6C7hK0lbA94ALgY8AlxXKTAGeBm7N8+8GfgxcQEoIdwUuk/RqREwvlPtT4AzgC8BTwKMVbeq0vsq2526Sp+VtxgJ/l8s8PkAsvgD8K7A28F/ARbndF+e2/wcwQ9KWha6YnbT3OmA+8HFgKbAZ8LbCvgda324P4PfAfa0FkkYDO+c4dETSrsBNpL/1wcC6wJeAcXke0hXFKZK2jYh7crk3kQaxuTIiVhSq3A+4IyJezNutQUpu/5IUu/8l3cO4TyF+FwK/o/5q4k9IVze3Kr5mMzNbNU4czcysW+sAu0bEQ/DaB/3vAptFxAJJPyQliu2J4zWthCEiZrRW5KtItwMTgb+h0I2RlJTsGRFz6hrURX2VbQcWRMTDkpYAa0TEXR3GYudW19l8ZfFfgCMi4vJCe2YCm5MSvQHbK2k9YCNgckS0Bpq5uVCmdn2FbYH5EfFqYdm6wBjgseKGuU2jCotWFJK2LwE/iYgphe0XAzdLel9EzAV+ADxPusJ4T95sb1JyWfx7QEocryzMnwB8CNi+rQvqZXlfuwJHA3OBe/NFyEsj4ry2eucBK4AdcOJoZjZo7r5hZmbdWtRKvLJf5J8T88+rgD2UBi5B0tbAn+Xl5GXjJJ0n6VfAK3k6Nm9XtHigpLHL+gZqe7cWtd1vuTD/vKVk2QZdtHcJKZm7QNIUSW9v2+9A68u8A6gaIbd9UJrPFNr1CnBcbvfapCuUV0sa3ZqAO/N228Jro51+B9YBk2gAAAP0SURBVDik1b2U9OXBr4CfFuIwgfRIjZl5fg3gROAr7fctvtbQiB9HhCJiq4jYOk/tSSMRsZyUvL6jJiZmZtYhJ45mZtat59vml+Wfa+Wf15KSiIPy/BRSd887C2Wm5eVnAnuR7nO7tFBHy9MdtqnT+gZqe7eq6nu+ZFlxH9OoaW++KrgXqXvupcBTku6Q9IFO1ldYi9Sltei5vKw9cb4it2n7tuXjSFciz2flxHIp8CZWfizGdFJX3p0lrQVMBma0jZy6L/BIRDyY598PjCddBR4KS1n1v62ZmRW4q6qZmQ2piHhJ0kxSYnQR6REQ32olDDmJ+DBwXERc0CpXMYjJgI9n6LK+xnXa3ohYAByU7w3cDTgdmClpYkS8OtD6kl0voe3qW0Qsl/RTUhJ6SmH50+SkXSs/WvF50t/kVGBWyT6Kg+jcmuuYCkwA3kp5N9WZhfkJ+edTJXWvirGk121mZoO0Wr6pmplZz5tBGpl0f2DjPN8yhvT+89rVrzwi5wGruK+hrm8Zw3uVqqv2RsQrEXEL8J+kxGpsN+sLHiTdF9nuXGBHSYcP1PA8YupdpPtZZ5dMTxS2XQFcDfw18DHS/ZXFgXnWJN3LWEwcWwnj5gO1ZSCSxpMGLPrlYOsyMzNfcTQzs+Exi/QsvwuBRyPi560VEfGCpLuBUyS9CLxKetj8C9SPClpqqOsDFgCTJR1I6mL7RPvjKAajk/bmQXbOIt0X+gipi+iJwH0RsWSg9RW7/nHe5/iI+E2hPd+XdC4wTdIHSaO1PksaOGevvNlLhXo+SxoI51XS6Ke/I3VJ3Q84OSKKidp04O9Jo+x+vq09uwMCbisse4A0qM3XJJ1CuoK5BbBhRJxY8bqqbEe6OvqTLsuZmVkJX3E0M7MhFxF/IN3rOIHCoDgFHyMlPJcDXyY9c+/yQexyKOs7H7iBdO/g3aRBa4baQO19itTN82TSCKXnk0ZkPaDD9WV+ROq2Oal9RUScQHqUxobAJaTBfc4nPc5i3+IzHiPiTlLSN550L+R1pGTyMdruSY2InwKLSAliWTfVmyJiaWH75cD+wL3AObnuY/J8tyYBt0XEc6tQ1szM2mjle9TNzMysX0n6MrBJROy3GrTll8CZEXHxMNQ9ijSC60kR8Y2hrt/M7I3IiaOZmdkbhKSJpHv+tm7rUtpXJE0FTgPem69impnZILmrqpmZ2RtERDwOHMX/j17arwQc7aTRzGzo+IqjmZmZmZmZ1fIVRzMzMzMzM6vlxNHMzMzMzMxqOXE0MzMzMzOzWk4czczMzMzMrJYTRzMzMzMzM6vlxNHMzMzMzMxqOXE0MzMzMzOzWk4czczMzMzMrNb/AYl3rpHEkwDaAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1080x720 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "bump = muons[(muons.M >= 80) & (muons.M <= 100)]\n",
    "\n",
    "fig = plt.figure(figsize=(15, 10))\n",
    "\n",
    "plt.hist(bump.M, 300, range = (80,100))\n",
    "\n",
    "plt.xlabel('Invariant mass (GeV/$c^2$)', fontsize = 15)\n",
    "plt.ylabel('Number of events \\n', fontsize = 15)\n",
    "plt.title('Measured distribution of events with two muons \\n', fontsize = 15)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA44AAAJ7CAYAAACs38+KAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3X28rfWcP/7Xu446DaZI9e0rqZlJFEmd+tbQzRTyZdwbhmbEIDOIDIrJd6yYGQy/ZIyb0I2ZMTJTlLsxbs+IDBUhhUQRIZkSiXPy+f1xrX2ss9v7Ouuc9j5r732ez8djP/a67t/Xta69Huu1P5/ruqq1FgAAAJjNZpMuAAAAgIVNcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY7AgldVg6pqVXX5LNMvH04fbOTSJqqqLqyqMzZguedUVRsZPnR4/O495vLbD9+TXdZjm62qnjMyvLKqzlqfunvW/fiqesoM4+dsG3Opqv66qr5XVb/ekPdvUqrqwVV17IS2vdZ7OVstVXVGVV24Aeu/x/Cc3ua21gqwVAmOwGJxc5Jdq2rF6Miq2i/JLsPpbJgvJDkwyRVjzr99kpelO+7jOjDJv69fWWN7fJKnzDD+WUleMk/b3CDD8/fEJP+Y5P5JXjHZitbLg5NMJDjm1u/lXNdyj3TntOAIMItlky4AYEw/Txdw/jjJaIvCHyf5RJJ9J1HUhqqq2yX5dWvtlknX0lr7aZL/no91V9VWrbVftNbmZf19WmuXbuxtjuGew99vHB53xrBA30uATYoWR2AxOTPJ46uqkmT4+/HD8bdSVQdV1X9V1U1VdV1Vva2q7jgyfceqOq2qvlVVv6iqb1TV31TVFtPW85Kq+mZV3VxVP6yqD1fV/xpOe8qwG+Ydpi1zZVW9dmR4ZVWdVVVHV9UV6VpI//dw2r2r6oNVdePw59+n1j+y/L2r6jPDGi6rqkeMc8Cqasuq+sequr6qflJVr0tyu2nz3KqralU9raouHR6XHw+P457D7qlfGc72yeFybdp6jqiq91XVz9K1rN2qq+rIdo4eHqtfDI/BXfvqGj2Ww9dnJHlskkOmapnqsjxTV9WqOqyqPjfyXr5p9L0b2eahw/fhZ8Pz41ljHOvNh90dv1NVv6yqr1bVk0amn5Hkn4eDN0xtp2d9s54XVXX7qvp5VT17huUuqKp/GRneuarOHL7/N1XVf1bV7iPTdxnW8viqOqWqbqiqq6vqxKrabDjPIMkLktx95DifMZy25/Bv4ifDmi6bqa6R7b2jqj4yMrz7cH3vGRm373DcbsPh0fd81lpGln9QVX15WM+nq2rPnnoOTfL+4eC3h+u7sqqWD9/H0ffwlcPpjxgZ94aq+szI8F2G+3jd8HivrGk9JWaoYeo9+OOqOr2qfjp8D/5kOP24qvp+VV1bVa+eel+G027VPXdkfX84Mu63quofquoHw/P/gqp68LTlpj6nnlTdZ95Pq+o/qmqnafPN+pkILF2CI7CYvCfJDkkeMBw+KMl2w/Frqar7J/lYkh8keVy6bm0PTXL6yGx3SfKTJH+Z5CFJXpPkqUneMLKeJyf5qyQnJTkiyV8k+WaS229A/fcfLn98koenCw+/l+QzSZYn+ZN0XS73TPL+qjUBeask/5nkDkmelORvkpycZOcxtvmqJE9P1yXyyCR3T/ele1ZVdXCSt6QLOf83yZ8lOT/J1kmuGa4nSZ6drgvqgdNWcWqSLyV5xPD1bA5Mcky64/+0JHslOWeMfRr1iiSfTPLFkVrePtOMw/Dw4SQ/Thc2X5bueM50HeTbhvvw6CQrk7yxqvZfRy0vT3JCkrem2/fPJHlnVT1xpNa/Gb4+bFjrF2aptfe8aK39PMkH0v3jZHS530myIsN/plTVnZN8OsnuSf58OP/tk3xseF6N+vskP0v39/IvSf56+Drpjum/pvt7mjrOU91s35/klmGdj0j397PmHzQzOC/JgVW1+XD44HT/SHnAyDwHJ/lha22m65r7akm6v4vXJPnbJE9M17X63VN/TzP4QpIXDl8/Zri+R7fWbk5yQbrPmdG6bp5h3Hkjw+ek+6x4YZInpPuu9cnhe7our073N/bY4TrfUVX/X5L90/0dnpzkuEx738f0tnSfb3+b7rz+bpIPVtUDps33f5I8J93nxNFJ9kl3TieZ889EYDFprfnx48fPgv5JMkjy4+Hrc9N180uSNyU5Z/j6x0kGI8ucl+ST09ZzWJKW5N6zbGdZuiBxc5IthuP+McnZPbU9ZbjOO0wbf2WS144Mr0zyiyQ7TJvvn5N8fWp7w3G7pfsi/rDh8LOSrEqy08g89x9u94ye2rYdbvP4kXGbJfla9/G/Ztyho8cl3Rfei3rWe+/h/IdOGz+1ntfNsExL8pxpx2NVkp1n2KeHzFTXtGXPGhk+K8nKGbY5fb4zk1yeZPORcY8fbuPAadt8+cg8t0tybZJX9RyTO6frTv2yaeM/lOTr6zpfZljfOOfFo4fD/3tknpek+2fI7YbDr0hyXZI7j8xzpyQ3JHn2cHiXYU3/NK2Gi5OcOTL82iRXTpvnLsNl77Mef8+7D5dZMRz+p3T/qFiV5J7Dce9J8u897+WtahmOPyPJ6iS7jYx71HB79+yp6Q+H8+wybfwrk1wyfL08yS/TfSb893DcNtPek4cM13PIyDpuPzx/TunZ/tR7cPrIuN8eHpPp5+znk7x72j5fOMv6/nA4fK8kv05y1Mg8myW5JMl/TjvONyS508i4Y4fr2mo43PuZ6MePn6X7o8URWGzOTPK4qtoyXWvIrbqpVtVvpWs1+LeqWjb1k67lZVWG10NW59gadskcTntnki3zm9a8i5M8dNhtb/+RVpINcVFr7YfTxj0wyXuT/Hqkzm+nC55T3dv2Hy579dRCrbXPJPnROrZ3n3Rfds8dWe7Xo8OzuDjJ/arqdVV1cE3rujuGD4453xdaa98ZqW1qn9bVsreh9k/y3rb2daVnpwsa01td1nSlbK1NfXnfKbO7d5Lfyq1vAPTuJPeoqu3Ws9Zxzov/SNdC+Ecjyz0h3T6uGlnPR5P8dGQ9Nya5aGQ9Uz4ybfjS9O9z0oXU7yZ5S1U9oaq2X9eOtda+nu59nmq1O3i4L18YGfeArN2Ktz6ubGu3VE5dH7mufZnJp5LsMWy5PSDd8X5zkn2GnzNT581UV9X9k/yotfZfUytov2kdnn6OzeTjI8v9NF3g/K9p5+w3k9x1+oLrsF+Sysj5Ofws+PcZ6rqgtfY/I8NTx29qm3P5mQgsIoIjsNi8L12Xzb9N95/8988wz52SbJ6uRXLVyM8v07Ue3W0437HpWi7em+SR6b70TV2btXz4+7R03bIen+RzSX5Y3XWQG/JlaXpoTLoWm+On1bkqye+M1Pm/MnNIXFdwnLrmaPp8vcu11j6WrkvbwelaIH5cVW+sqnG7os20nzOZbZ92HHP59bVjptU2/EJ+XboWw1HXTxv+VX5zTsy27kxf/8jw9PWvyzrPi9Z1pTw3XVhMddct3jdr/zPlLsPp09fzB/nN+TVlffd5Knw8OF230dOS/KCqzquq+61j/85LclBV3S3dP2k+PTLuXum6oG9ocJxpP5J17Msszk/X2vaAdKH2M+mC1A3pguRB6Vokp7a5Y2Y+r3+Y8c6BmWpf7/dlBjsm+Vlr7aYZ6vqt4T/i+mpI5uczEVhE3FUVWFRaaz+vqg8keX66rmw/n2G269N92Ruk6yo43feHv/8oXfe3E6YmVNUe07b36ySvS/K64ZfcI9OF1qvTda+begzI9Fa5O81U/gzjfpIuuM50Xd6Ph79/kN/cjXPUulp3fjAy30/WY7m01t6R7vqq7dJd9/W6dC1VL17Xspl5P2cyUx3bp7vGK+k/tj/O+rtm+jaHX3a3zdrHZ0NM1bx9uiA6ZYfh7/Vd/zjnRdK1aL6/qnZOFxCvTXeX4dH1vC8zP/bjxvWsaUatta8leWx1dwo+KN11eh+sqp2Gfz8zOS/d9aAHJ7m0tXZdVZ2X7hq+zyT5abprTCeqtXZDVX053X7tna5bZ6uqTw/HTb++8Vbn2NAOue3n2Gxuzro/f65Jcoeq+q1p4XGHJDe11n457sbG+EwEligtjsBi9OZ0LY0zfkkZhsn/TrJ7a+3CGX6mguNW6VohRx2ZWbTWvttae1W6rmJTAXOq++i9puarqv+T7vqkcXw83U1PLpqhziuH81yQZN/ROxsOb/6zrgD4lXRfKh85stxmo8Pr0lq7trV2Srovx1P7fFtacEbtMww8U7VN7dPnh6NmOrZ3y61D9LgtMJ9L8uhpLSOPSfdP1E+vX+m3ckmSm7J2t9Gka5X5Rmvt2vVc3zjnRdJ1L71+uJ0npPtHyC0zrOerM6zn6+tZU+9xbq2taq19It1NU3ZM/zMRP5WuVfHo4eukO8funu464/N7Quc6a9kAfef0p/Kbmxl9amTcEem6vY8Gx88l2X54g6kka7rOPyy3/RybzdVJdqmq0dofPG2eC9L9Q2fqZkdTd6V+3G2pa5bPRGCJ0uIILDqttZXpulD2OS7Jx6vq1+lunnJjui5xD0tyQmvtG+mu/XpuVX0uyRXpQuNadz6sqlPStRT8d7ruaX+Q7iYlxw9n+XyS7yX5h6r6f+m6ox2XrsVkHIPhOj5YVaela026a5IHpbvxzcp0d4J96XCeQbrA+4qso9Vt2Irz1iQnVtXqJF9N8ox0XX1nVVUnDvdj5XAb90tySH7T2viddDfdOaqqbkiyqrV24QyrWpdrh/v0snRf2F+d7rrHDw/rv7q6xwy8oqpuSvfPzr/KrVtuvpbkkVX1qHRfor8/8s+BUX+T7u6r51TVm9Nd8/bqdK1In92A+tdorf2kqk5O8tLhsb4wXSh9aLo7e66vQdZ9XqS1tqq6x1j8ZbqwNv2xISelu9vpJ6rqDenO1R3SvZ+fbq29az1q+lqSHarqKemC8o/T/YPktelaPr+VrqXr+CRfaq31tbB9Kd3fyMHp/hE0dQwvHY47oWfZGWuZFqjX11SIfmZVnZmuFW7qsTPnJXluuusbvzAy7qSR1xnuw39W1fnp7uL64nStzy9M9zf7mttQX59z0t3R9+3VPZbkfunuwLpGa+2yqnpXkn+s7pFEV6T7LLhnuruijm2Mz0RgidLiCCxJrbVPp/sCul26O1S+P12g+25+c93Zy5O8K12geFe6VofnTlvVZ4frOT1dt9dHJ3lGa+2c4XZ+NRw3FVBfkO6L2P9kDMMAe0C61qq3prtJyInpWkK/OZznpnStGz9Pd/3ay4bbuWqMTRyX7pqkvx7u4/fzmy+8s7kgXevBW9I9BuQv0gWZ1w/ruTndl859k/zXcP4NcX6SN6brnnhqugDwqGnzPDFdUP2XJH+X7j2b3lL2pnQtb6cNazl6po211r6a7vEi26e7a+fU+/64mebfAH+d7i6cf5HuZigHJ/mT1tqMzxntM855MeLMdKHx+5l2XWBr7cfD9XwtXffCj6R77MbWSb68nmX9W7o7eP59uuM8SNcd+ofpgt5/pHsvLkv3WI6+/ft1uvc/+U0rXkbqX1cr2Ey1bLDW2lXpAt5j0nWVHb12eqqmz7bWVg9ffzFdkPz2DP+keFS6f0qdnO7mM5XksNba9PdtTrTWLkkXFA9M1y35kHTXKE/3jCTvSHeenpuudfcPh5+V66P3MxFYuqq1cS9FAQAAYFOkxREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAAoJfgCAAAQC/BEQAAgF6CIwAAAL0ERwAAAHoJjgAAAPQSHAEAAOglOAIAANBLcAQAAKDXskkXMEl3uctd2i677DLpMgAAACbioosu+nFrbbt1zbdJB8dddtklF1544aTLAAAAmIiqumqc+XRVBQAAoJfgCAAAQC/BEQAAgF6b9DWOAADAxrFq1apcffXVufnmmyddyiZp+fLl2WmnnXK7291ug5YXHAEAgHl39dVX5453vGN22WWXVNWky9mktNZy3XXX5eqrr86uu+66QevQVRUAAJh3N998c7bddluhcQKqKttuu+1tau0VHAEAgI1CaJyc23rsBUcAAGCTcIc73GGt4TPOOCPPec5zkiRvectb8k//9E+9y4/Ov6lxjSMAALDRDQaDBbW+P//zP5+bQubJ6tWrs2zZ5OKbFkcAAGCTNxgM8trXvjZJcsEFF2SvvfbK3nvvnRe96EW5973vvWa+73//+3nIQx6S3XbbLccdd9yM67roootyyCGHZN99980RRxyRa665Jkly6KGH5vjjj8/++++fe9zjHjnvvPOSJLfcckte9KIXZb/99stee+2VU045JUmycuXKHHTQQXnEIx6RPfbYI0nyile8Irvvvnse8IAH5IlPfGJe+9rX5oorrsg+++yzZvuXX375WsNzQYsjAACwSfjFL36Rvffee83wT37ykzziEY+41XxPfepT87a3vS0HHnhgXvziF6817eKLL84Xv/jFbLnlltl9991zzDHH5G53u9ua6atWrcoxxxyTc889N9ttt13e/e5354QTTshpp52WpGs5/PznP58PfehDOfHEE/Oxj30sp556arbeeutccMEF+eUvf5n73//+efCDH5wk+cIXvpBLLrkku+66ay644IKcffbZ+dKXvpRVq1Zln332yb777pvf/d3fzdZbb52LL744e++9d04//fQ89alPndNjJzgCAACbhK222ioXX3zxmuEzzjgjF1544VrzXH/99bnxxhtz4IEHJkme9KQn5QMf+MCa6Ycffni23nrrJMkee+yRq666aq3g+PWvfz2XXHJJHvSgByXpWhN33HHHNdMf85jHJEn23XffXHnllUmSj3zkI/nyl7+cs846K0lyww035PLLL88WW2yR/ffff80jND7zmc/kkY98ZJYvX57ly5fn4Q9/+Jr1Pv3pT8/pp5+ek046Ke9+97vz+c9//rYdrGkERwAAgDFtueWWa15vvvnmWb169VrTW2vZc88989nPfrZ3+dFlW2t5wxvekCOOOGKteVeuXJnb3/72Y9X12Mc+NieeeGIOO+yw7Lvvvtl2223H3qdxuMYRAABgaJtttskd73jHfO5zn0uSnHnmmeu1/O67755rr712TXBctWpVvvrVr/Yuc8QRR+TNb35zVq1alST5xje+kZ///Oe3mu/+979/3v/+9+fmm2/Oz372s7VaQpcvX54jjjgif/EXfzHn3VQTLY4AAABrOfXUU/OMZzwjm222WQ455JA1XVPHscUWW+Sss87Kc5/73Nxwww1ZvXp1jj322Oy5556zLvP0pz89V155ZfbZZ5+01rLddtvlnHPOudV8++23Xx7xiEdkr732yg477JD73Oc+a9V25JFH5r3vfe+a6yPnUrXW5nyli8WKFSva9D7NAADA3Lvssstyr3vda9JljOVnP/vZmmc+vupVr8o111yT17/+9ROuqjNV20033ZSDDz44b33rW9fcQfW1r31tbrjhhrziFa+YcdmZ3oOquqi1tmJd29XiCAAAMOKDH/xgXvnKV2b16tW5+93vnjPOOGPSJa1x9NFH59JLL83NN9+co446ak1ofPSjH50rrrgin/jEJ+Zlu1octTgCAMC8W0wtjkvVbWlxdHMcAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAANglTj9hg/XkcBwCwYAwGgxlfA0vQXP+NL8HPjNWrV2fZsoUR2bQ4AgAAm6wrr7wyhx12WPbaa68cfvjh+c53vpNbbrklu+66a1pruf7667P55pvnU5/6VJLk4IMPzuWXX77WOm655Za86EUvyn777Ze99torp5xySpJk5cqVOfTQQ/O4xz0u97znPXPkkUdm6nGIF110UQ455JDsu+++OeKII3LNNdckSQ499NAce+yxWbFiRV7/+tfniiuuyAEHHJD73Oc+eelLX7qm1fTJT35yzjnnnDU1HHnkkTn33HPn7TgJjgAAwCbrmGOOyVFHHZUvf/nLOfLII/Pc5z43m2++eXbfffdceuml+fSnP5199tkn5513Xn75y1/mu9/9bnbbbbe11nHqqadm6623zgUXXJALLrggb3vb2/Ltb387SfLFL34xJ598ci699NJ861vfymc+85msWrUqxxxzTM4666xcdNFF+bM/+7OccMIJa9b3q1/9KhdeeGFe8IIX5HnPe16e97zn5Stf+Up22mmnNfM87WlPyxlnnJEkueGGG3L++efnYQ972Lwdp4XR7gkAADABn/3sZ/Oe97wnSfKnf/qnOe6445IkBx10UD71qU/l29/+dl7ykpfkbW97Ww455JDst99+t1rHRz7ykXz5y1/OWWedlaQLcpdffnm22GKL7L///msC3957750rr7wy22yzTS655JI86EEPStK1WO64445r1veEJzxhrfqmWhaf9KQn5YUvfGGS5JBDDsmznvWsXHvttTn77LPz2Mc+dl67tQqOAAAA0xx88MF585vfnO9///t5+ctfnte85jVZuXJlDjrooFvN21rLG97whhxxxBFrjV+5cmW23HLLNcObb755Vq9endZa9txzz3z2s5+dcdu3v/3tx6rxyU9+cv7lX/4lZ555Zk4//fT12Lv1p6sqAACwyfr93//9nHnmmUmSd77znWuC4f7775/zzz8/m222WZYvX5699947p5xySg4++OBbreOII47Im9/85qxatSpJ8o1vfCM///nPZ93m7rvvnmuvvXZNcFy1alW++tWvzjjvAQcckLPPPjtJ1tQ55SlPeUpOPvnkJMkee+yxPru93gRHAABgk3DTTTdlp512WvNz0kkn5Q1veENOP/307LXXXvnnf/7nvP71r0+SbLnllrnb3e6WAw44IEnXdfXGG2/Mfe5zn1ut9+lPf3r22GOP7LPPPrn3ve+dZz7zmVm9evWsdWyxxRY566yzcvzxx+e+971v9t5775x//vkzznvyySfnpJNOyl577ZVvfvOb2XrrrddM22GHHXKve90rT33qU2/LYRlLTd3VZ1O0YsWKduGFF066DABgyOM4YOm67LLLcq973WvSZSw6N910U7baaqtUVc4888y8613vWnP31Jtuuin3uc998oUvfGGtQDmbmd6DqrqotbZiXcu6xhEAAGCBuuiii/Kc5zwnrbVss802Oe2005IkH/vYx/K0pz0tz3/+88cKjbeV4AgAALBAHXTQQfnSl750q/EPfOADc9VVV220OlzjCAAAQC/BEQAA2Cg25furTNptPfaCIwAAMO+WL1+e6667TnicgNZarrvuuixfvnyD1+EaRwAAYN7ttNNOufrqq3PttddOupRN0vLly7PTTjtt8PKCIwAAMO9ud7vbZdddd510GWwgXVUBAADoJTgCAADQS3AEAACgl+AIAABALzfHAQBmNRgMZnwNwKZFiyMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHANiEDAaDDAaDSZcBwCIjOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHAAAAei2bdAEboqquTHJjkluSrG6traiqOyd5d5JdklyZ5PGttf+ZVI0AAABLxWJucfyD1trerbUVw+EXJ/l4a223JB8fDgMAAHAbLebgON0jk7xj+PodSR41wVoAAACWjMUaHFuSj1TVRVV19HDcDq21a4avf5Bkh5kWrKqjq+rCqrrw2muv3Ri1AsAmazAYZDAYTLoMAG6jRXmNY5IHtNa+V1XbJ/loVX1tdGJrrVVVm2nB1tpbk7w1SVasWDHjPAAAAPzGomxxbK19b/j7R0nem2T/JD+sqh2TZPj7R5OrEAAAYOlYdMGxqm5fVXecep3kwUkuSfK+JEcNZzsqybmTqRAAAGBpWYxdVXdI8t6qSrr6/7W19uGquiDJv1XV05JcleTxE6wRAABgyVh0wbG19q0k951h/HVJDt/4FQEAACxti66rKgAAABuX4AgAAEAvwREAAIBegiMAAAC9Ft3NcQCAyRoMBjO+BmDp0uIIAABAL8ERAACAXoIjAAAAvQRHAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAAoNeySRcAACwOg8Fg0iUAMCFaHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADotWzSBQAA82swGEy6hBmN1rVQawSgo8URAACAXoIjAAAAvQRHAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAAoJfgCAAAQC/BEQDYYIPBIIPBYNJlADDPBEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHAAAAegmOALAIDAaDDAaDSZcBwCZKcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQa9mkCwCApWIwGMz4eiFaV60LvX4ANi4tjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHAAAAegmOAAAA9BIcAQAA6LVs0gUAAHNnMBjM+BoAbgstjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6LZt0AQBAMhgMZny9VCz1/QNY6rQ4AgAA0EtwBAAAoJfgCAAAQC/BEQAAgF6CIwAAAL0WbXCsqs2r6otV9YHh8K5V9bmq+mZVvbuqtph0jQAAAEvBog2OSZ6X5LKR4VcneV1r7feS/E+Sp02kKgAAgCVmUQbHqtopycOSvH04XEkOS3LWcJZ3JHnUZKoDAABYWhZlcExycpL/Hc5PAAAgAElEQVTjkvx6OLxtkutba6uHw1cnueskCgMAAFhqlk26gPVVVX+Y5EettYuq6tANWP7oJEcnyc477zzH1QHA3BkMBpMuYcEaPTaOE8D8W4wtjvdP8oiqujLJmem6qL4+yTZVNRWEd0ryvZkWbq29tbW2orW2YrvtttsY9QIAACxqiy44ttZe0lrbqbW2S5I/TvKJ1tqRST6Z5HHD2Y5Kcu6ESgQAAFhSFl1w7HF8kr+sqm+mu+bx1AnXAwAAsCQsumscR7XWViZZOXz9rST7T7IeAACApWgptTgCAAAwDwRHAAAAegmOAAAA9BIcAQAA6LWob44DAEtR3wPtN6UH30/t31LfT4DFQIsjAAAAvQRHAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAAoJfgCAAAQC/BEQAAgF7LJl0AALBhBoPBWr/pN3qcHDOA9aPFEQAAgF6CIwAAAL0ERwAAAHoJjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBeyyZdAABsKgaDwVq/N1W3Zf839WMHMClaHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0WjbpAgBgYxp9gPykHiY/19ud1H7MZqHVA8Btp8URAACAXoIjAAAAvcYKjlW1fVXtOjJcVXV0VZ1cVQ+fv/IAAACYtHFbHM9I8vyR4ZcneVOShyR5b1U9ZW7LAgAAYKEYNzjuk+QTSVJVmyX58yR/1Vq7Z5K/TXLs/JQHAADApI0bHLdOct3w9b5J7pzkncPhTyT5vTmuCwAAgAVi3OB4dZI9hq8fluRrrbXvDYe3TnLzXBcGAADAwjDucxxPS/L3VfXAdMHxJSPTDkhy2VwXBgAAwMIwVnBsrb2yqr6XZL8kx6QLklPunOTt81AbAAAAC8BYwbGqdk7yrtbaP80w+ZgkO85pVQCwEQwGg7V+Mznz9R7c1vU6RwA6417j+O0k95tl2l7D6QAAACxB4wbH6pm2PMkv56AWAAAAFqBZu6pW1V5J9h4Z9dCquue02ZYneXySb8xDbQAAACwAfdc4PjrJy4avW5K/nmW+byd55lwWBQAAwMLR11X175LcMclvp+uqethwePRny9ba77bWPjbfhQIAADAZs7Y4ttZWJVk1HBz3WkgAAACWmLEexzGlqu6RZKd01zaupbX2obkqCgAAgIVj3Oc47pHkzCR7ZuY7rLYkm89hXQAAACwQ47Y4npJkyySPSXJpkl/NW0UAsACMPvDdw98B2NSNGxzvl+SPW2sfmM9iAAAAWHjGvenNFZnhukYAAACWvnGD4wuS/FVV/c58FgMAAMDCM25X1VcmuWuSr1XVlUmunz5Da23/OawLAACABWLc4HjJ8AcAAIBNzFjBsbX21PkuBAAAgIVp3GsckyTVuVtV/X5V3X6+igIAAGDhGDs4VtWzknwvyVVJzkuy+3D8e6rq2PkpDwAAgEkbq6tqVb0oySuSvDrJJ5N8YmTyyiRPTHLyXBcHAAvJYDCY8fViN7Uvt2Wf5vN4bKx1L6X3FGCujXtznGcn+evW2t9X1ebTpn09yT3mtiwAAAAWinG7qv6vJBfNMu3XSZbPTTkAAAAsNOMGx28mOWSWaQcnuXRuygEAAGChGber6slJ3lRVv0py1nDc9lX1tCR/meQZ81EcAAAAkzfucxzfXlV3SvLXSU4cjv5QkpuSDFpr/zpP9QEAADBh47Y4prX2mqp6S5IDk9wlyU+SfLa1dsN8FQcAAMDkjfs4jt9prX2rtXZjko/Mc00AAAAsIGPfHKeqPl9Vz6+qnea1IgAAABaUcYPjw5NcluRlSa6sqvOq6tlVtcP8lQYAAMBCMO7NcT6Y5INVtUWShyZ5fJJXJTm5qj6V5F2ttbfPX5kAMH8Gg8GMr2/LMlPjxl0fACxk47Y4Jklaa79qrZ3TWntSku2THJXknklOmY/iAAAAmLyx76o6pao2S3JYkickeXSSOyU5f47rAgAAYIEYu8Wxqg6pqjcluSbdnVXvm+Tvkty9tXbQPNUHAADAhI37OI5r0nVN/UqSk5O8u7X2rfksDAAAgIVh3K6qb0kXFr82n8UAAACw8Ix7V9UTp15XVSXZMcmPWmur56swAAAAFob1ucbxoVX1uSQ3J/lOkr2G499aVX8yT/UBAAAwYWMFx6p6cpL3JflakqOnLXd5kqfNfWkAAAAsBONe43hCkte01l5SVZsnOX1k2leTvHDOKwOACRgMBnO6zIasDwAWmnG7qt49yUdnmXZzkt+em3IAAABYaMYNjt9Ncr9Zpq1I8s25KQcAAICFZtzgeGqSlw1vgrPVcFxV1eFJjkvytvkoDgAAgMkb9xrHVye5W5J3JLllOO78JJsnOaW19g/zUBsAAAALwLjPcWxJnl1VJyU5PMldkvwkySdaa9+Yx/oAAACYsHFbHJMkrbUrklwxT7WMpaqWJ/lUki3T1X9Wa+1lVbVrkjOTbJvkoiR/2lr71eQqBQAAWBrGvcZxIfllksNaa/dNsneSh1TVAem6076utfZ7Sf4nni0JAAAwJxZdcGydnw0Hbzf8aUkOS3LWcPw7kjxqAuUBAAAsOYsuOCZJVW1eVRcn+VG650tekeT61trq4SxXJ7nrpOoDAABYSma9xrGqdk5yTWtt1UasZyyttVuS7F1V2yR5b5J7jrtsVR2d5Ogk2XnnneenQAAmajAYrPWbTc+k3vvR7Tr/gKWkr8Xx20nulyRV9YmqGjucbSytteuTfDLJgUm2qaqpILxTku/NssxbW2srWmsrtttuu41UKQAAwOLVFxx/keS3hq8PTfLb817NGKpqu2FLY6pqqyQPSnJZugD5uOFsRyU5dzIVAgAALC19j+P4YpLXV9VHh8PHVNU1s8zbWmvHz21ps9oxyTuqavN0wfffWmsfqKpLk5xZVX+TrvZTN1I9AAAAS1pfcHxGktckeWS6u5Yenu5RGDNpSTZKcGytfTnDLrTTxn8ryf4bowYAAIBNyazBsbX2tSQPT5Kq+nWSR7XWPr+xCgMAAGBh6GtxHLVrktm6qQIAALCEjRUcW2tXVdWyqnpCkgckuXOSnyQ5L8l7Rp6fCAAAwBIzVnCsqu2TfCTJXkmuTPLDdI/AeHaSL1XVg1tr185XkQAAAEzOuF1VT0qybZIDRq9zrKr9kpw9nP6nc18eALCxLPYH1m+s+qe2s9iPF8D66HuO46iHJjl++s1xWmsXJHlJkofNdWEAAAAsDOMGxy2T3DjLtBuTbDE35QAAALDQjBsc/zvJ8VV1+9GRw+Hjh9MBAABYgsa9xvEFST6Z5LtV9ZF0N8fZPskRSSrJofNSHQAAABM3Votja+3iJLsleWuS7ZI8KF1wfEuS3VprX5q3CgEAAJiocVsc01r7cZIXz2MtAAAALEDjXuMIAADAJkpwBAAAoJfgCAAAQK+xr3EEABaXwWAw6RIAWCLW2eJYVVtW1QlVdd+NURAAAAALyzqDY2vtl0lOSLLN/JcDAADAQjPuNY6fS7LPfBYCAADAwjTuNY7HJfnXqlqV5ENJfpikjc7QWrtpjmsDAABgARg3OH5u+Psfkrx+lnk2v+3lAAAAsNCMGxz/LNNaGAEAANg0jBUcW2tnzHMdAAAALFDr9RzHqtojyb5J7pbktNbaD6rq95L8sLV243wUCAAAwGSNFRyr6g5JTkvyuCSrhst9OMkPkvxdku8keeE81QjAJmzqIfaTfJj9JLfNwuAcADZ14z6O46Qkv5/k8CR3TFIj0z6U5CFzXBcAAAALxLhdVR+T5HmttU9W1fS7p16V5O5zWxYAAAALxbgtjlsluW6WaXdMcsvclAMAAMBCM25wvCDJk2eZ9rgk589NOQAAACw043ZV/X9JPlpVH0vy7+me6fjQqnp+uuB48DzVBwAAwISN1eLYWjsv3Y1xtkzyj+lujnNikt9J8sDW2gXzViEAAAATNfZzHFtrn0lyUFVtleROSa5vrd00b5UBAACwIIx7jeOom9M9y/EXc1wLAAAAC9DYLY5V9dAkL02y73C51VV1UZK/ba19cJ7qAwCYc4PBYNIlACwqY7U4VtUzk7w/yc+SPC/JHw1//yzJ+4bTAQAAWILGbXH8qySntNaeNW38W6rqLUlOSHLKnFYGAADAgjDuNY7bJnnvLNPOTnLnuSkHAACAhWbc4PjJJIfMMu2QJJ+am3IAAABYaGbtqlpVe4wM/kOSt1fVtknOSfKjJNsneXSS/5vk6fNZJAAAAJPTd43jJUnayHAleebwpw2Hp3w4yeZzXh0AAAAT1xcc/2CjVQEAAMCCNWtwbK3918YsBAAAgIVp3MdxrFFVy5JsMX18a+2mOakIAACABWWsu6pW1dZV9aaquibJzUlunOEHAACAJWjcFscz0j12421JvpnkV/NVEAAAAAvLuMHx8CTPbK29az6LAQAAYOEZq6tqku8kcQ0jAADAJmjc4HhckpdW1c7zWQwAAAALz1hdVVtrH6qqByb5ZlVdmeT6GebZf45rAwAAYAEYKzhW1WuTHJvkgrg5DgAAwCZl3JvjPD3JCa21V85nMQAAACw8417jeFOSi+azEAAAABamcVscX5/k6Kr6aGutzWdBADBXBoPBjK9hLji/gE3JuMHxLkn+T5KvV9XK3PrmOK21dvxcFgYAAMDCMG5wfFyS1Ulul+RBM0xvSQRHAACAJWjcx3HsOt+FAAAAsDCNe3McAAAANlHjPsfxWeuap7X2ptteDgAAAAvNuNc4/mPPtKm7rAqOAAAAS9BYXVVba5tN/0ly5yRPTPKlJHvMZ5EAAABMzrgtjrfSWrs+yburauskpyQ5dK6KAgAAYOGYi5vjfDvJijlYDwAAAAvQBrc4JklV7ZjkBenCIwBsFIPBYKxxsCGcSwC3Nu5dVa/Nb26CM2WLJHdMcnOSx8xxXQAAACwQ47Y4vjG3Do43J7k6yYdba9fNaVUAAAAsGGMFx9baYJ7rAAAAYIGai5vjAAAAsITN2uJYVZ9Yj/W01trhc1APAAAAC0xfV9VxrlvcMcnv59bXPwIAALBEzBocW2t/NNu0qto5yfFJ/jDJj5O8bu5LAwAAYCFYr+c4VtXvJXlJkj9J8qPh61Naa7+Yh9oAAABYAMZ9juOeSU5I8kdJvpvkeUlOa639ah5rA2AT5QHsALCw9N5Vtar2rar3JPlykn2SPD3Jbq21twiNAAAAm4a+u6r+R5IHJ/lKkj9urf37RqsKAACABaOvq+oRw987JXljVb2xb0Wtte3nrCoAAAAWjL7geOJGqwIAAIAFq+9xHIIjAAAA/TfHAQAAAMERAACAXoIjAAAAvfpujgMAi8ZgMLhN08E5AjA7LY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAr0UXHKvqblX1yaq6tKq+WlXPG46/c1V9tKouH/6+06RrBQAAWAoWXXBMsjrJC1preyQ5IMmzq2qPJC9O8vHW2m5JPj4cBgAA4DZadMGxtXZNa+0Lw9c3JrksyV2TPDLJO4azvSPJoyZTIQAAwNKy6ILjqKraJcn9knwuyQ6ttWuGk36QZIcJlQUAALCkLJt0ARuqqu6Q5Owkx7bWflpVa6a11lpVtVmWOzrJ0Umy8847b4xSARjTYDCYdAkwUVN/A/4WgIVmUbY4VtXt0oXGd7bW3jMc/cOq2nE4fcckP5pp2dbaW1trK1prK7bbbruNUzAAAMAituiCY3VNi6cmuay1dtLIpPclOWr4+qgk527s2gAAAJaixdhV9f5J/jTJV6rq4uG4v0ryqiT/VlVPS3JVksdPqD4AAIAlZdEFx9bap5PULJMP35i1AAAAbAoWXVdVAAAANi7BEQAAgF6CIwAAAL0ERwAAAHotupvjALBhRh8oPqmHiy+EGlh6Fvq5NB/n/dR6Fvq+A0uHFkcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXoIjAAAAvQRHAAAAei2bdAEAMI7BYDDpEgBgk6XFEQAAgF6CIwAAAL0ERwAAAHoJjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAECvZZMuAIClZTAYrPV73PlhU7O+fytLtQZgcdDiCAAAQC/BEQAAgF6CIwAAAL0ERwAAAHoJjgAAAPQSHAEAAOglOAIAANBLcAQAAKDXskkXAACw2A0Gg3lbdmr6uNsYne+21DVX1rd+YGHS4ggAAEAvwREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAr2WTLgAAgLUNBoN5W99M617XdAAtjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAEAvwREAAIBegiMAAAC9BEcAAAB6LZt0AQAsfLM9HHzq9boeGL6uB47DpmIxnPfj/l2Psw5g6dDiCAAAQC/BEQAAgF6CIwAAAL0ERwAAAHoJjgAAAPQSHAEAAOglOAIAANBLcAQAAKCX4AgAAECvZZMuAICNbzAYrPV7+utxll3XONjULfW/i5k+R4ClS4sjAAAAvQRHAAAAegmOAAAA9BIcAQAA6CU4AgAA0EtwBAAAoJfgCAAAQC/BEQAAgF7LJl0AAEuTh4IDwNKhxREAAIBegiMAAAC9BEcAAAB6CY4AAAD0EhwBAADoJTgCAADQS3AEAACgl+AIAABAL8ERAACAXssmXQDAJAwGgxlfL3ZT+7JU9w9Y27r+1ufq73+mz5Zxa1hMltK+wFzT4ggAAEAvwREAAIBegiMAAAC9BEcAAAB6LcrgWFWnVdWPquqSkXF3rqqPVtXlw993mmSNAAAAS8WiDI5JzkjykGnjXpzk46213ZJ8fDgMAADAbbQog2Nr7VNJfjJt9COTvGP4+h1JHrVRiwIAAFiiFmVwnMUOrbVrhq9/kGSHSRYDAACwVCybdAHzobXWqqrNNK2qjk5ydJLsvPPOG7UuYOnoexj2UtgezJVDV65Mkqw89NCJ1sH6m+vPm9H1zfSZ5vMNFral1OL4w6raMUmGv38000yttbe21la01lZst912G7VAAACAxWgpBcf3JTlq+PqoJOdOsBYAAIAlY1EGx6p6V5LPJtm9qq6uqqcleVWSB1XV5UkeOBwGAADgNlqU1zi21p44y6TDN2ohAAAAm4BF2eL4/7d379Gy1NWBx7/bexUko4uHxBAgjoqjomah4JPRORkNXiAKiQygRiHquCajM6OTmUAWM+SMrqyJr9EYQ1AjEggBjI8IAyqIHgWVBBx5XAT08kjk8jBK0KiR554/6tdYt6mq7j7dp+ucc7+ftXr16ar6/WrX7urq2qeqqyRJkiRJ82PhKEmSJEnqZOEoSZIkSepk4ShJkiRJ6rQmL44jSZOY9w2mR82v6cbX0lqwsLQEwNLCwljTaXVaDdue1RCDpMl4xFGSJEmS1MnCUZIkSZLUycJRkiRJktTJwlGSJEmS1MnCUZIkSZLUycJRkiRJktTJwlGSJEmS1MnCUZIkSZLUycJRkiRJktRpY98BSNJasbi42Dmsafys+2saP2o6SeuDn2tJffKIoyRJkiSpk4WjJEmSJKmThaMkSZIkqZOFoyRJkiSpk4WjJEmSJKmThaMkSZIkqZOFoyRJkiSpk4WjJEmSJKnTxr4DkKY1zQ3YNZ5Z3eTe92e2ZpHP7f09WVhaAmBpYaGX+a70vCddvnHjmlf8UpN5fa809b29bzO1ffOIoyRJkiSpk4WjJEmSJKmThaMkSZIkqZOFoyRJkiSpk4WjJEmSJKmThaMkSZIkqZOFoyRJkiSpk4WjJEmSJKmThaMkSZIkqdPGvgOQNB+Li4vbPM9rfrPqZyXjbppH29/TzmNW02n7sLC09ODfSwsLM+9z1n1PGkPTfJviq5t3rAMr8V7M2jy2H/PeRo27LW4atxa2p/P+bpam4RFHSZIkSVInC0dJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUqeNfQeg2Rh1A9lpbjA7zxuwj5xPfdwqWpam+EfdjHjWN5hfT8bN3axv+Dzq86PZaroJfNeN4ZvajjPtJLHU+xt32LRxjbvMTfOYdF7L0TS/lehnmvlMmsOVNKt1c9b9zCs3s9oGz3r7Pu48VmvffcWw3vOlyXjEUZIkSZLUycJRkiRJktTJwlGSJEmS1MnCUZIkSZLUycJRkiRJktTJwlGSJEmS1MnCUZIkSZLUycJRkiRJktRpY98BaHVZzs3rV6V6fCt4w9p5WM4NkZcV46DNMtqupvWhz1hmPe9xbwA+qxuF96Xpxu8ruUwrcUP7WcQ4q7iW0/dKznueRr0nq229GrRZi5/btWDc78+V+N6YVd/j9jPpfGa179CU4zW1vzjCSq4ja5FHHCVJkiRJnSwcJUmSJEmdLBwlSZIkSZ0sHCVJkiRJnSwcJUmSJEmdLBwlSZIkSZ0sHCVJkiRJnSwcJUmSJEmdLBwlSZIkSZ029h2AZmtxcbHx73HbzGJ+y4mhqc9RbaeJuy9tMTcNb5x2MGyK3LSNW1haqp7H7GdpoTZl+XvQx0PGj2Fe72fTfAZxTxrzcue3WjXloWvY8PCVVp/vctqMG2vTfEbNezltJjXr/lbb/Ma1FvOwEp+ZcT+vKxnDNGa13Z3Hvs6sdO3fTLLvNKv9tnlayfmOyt24+5WzNu3+8GrkEUdJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUicLR0mSJElSp419B6DJLTXcbL1rumrihbbJtjHujY7HvpFp/eank7adwlJtORZqwwfLt9gQFw03ah1189bGGxhPcYP5ke/poJ853Ui3a/zCMvoeLH89D6Pm15S7sd+LMfqdZHzTjbbHNe57P68bc497o/CuYcvpZyWWb965W436vLn7pJ+FWbUdt++2fHTNe96f13HzsJL5aprPNMs3r3Wyr8//uDd3bxw35vfsyH5mbNQyLSeGSdtMMt+J8z5iuknbLKftWucRR0mSJElSJwtHSZIkSVInC0dJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUicLR0mSJElSJwtHSZIkSVInC0dJkiRJUqeNfQegBouLP/tz6LluaWHhZy/qf3dYapquaVgthoWlpdZ+msZtM2xEXA/G0zDd4ogYBi0al2kCS6Xvpdr8xo2hq7+mZar3MU3cTf10vU/18U3DltP3NBrnURvflJtp4ppV/LPOwyTzWM77N7BS7+Oo+U4z3Ur0PY/3bzXYXpZzNVgN69xKfqaa2oz67hp3uzRSw3fyuDF0zbdpG7rYMK/2sB46bdOwzmWe4Pu/a5lH7VNMslyzaDtos5x9nXHnN2rfd9x+xn0flzt927TTvCergUccJUmSJEmdLBwlSZIkSZ0sHCVJkiRJnSwcJUmSJEmd1l3hGBGbIuL6iNgSEcf3HY8kSZIkrXXrqnCMiA3AnwAHA/sCr4yIffuNSpIkSZLWtnVVOALPAbZk5o2ZeQ9wFnBYzzFJkiRJ0pq23grHPYHv1F7fUoZJkiRJkpYpMrPvGGYmIo4ANmXmG8rr1wDPzcw316Z5I/DG8vLJwPVzD3S0xwDf6zuI7ZS574+575f574+574+574+574+5789qzf3jMnP3URNtnEckc7QV2Lv2eq8y7EGZ+SHgQ/MMalIRcXlmHtB3HNsjc98fc98v898fc98fc98fc98fc9+ftZ779Xaq6mXAkyLi8RHxCOBo4JyeY5IkSZKkNW1dHXHMzPsi4s3A54ANwCmZeU3PYUmSJEnSmrauCkeAzDwfOL/vOKa0qk+lXefMfX/Mfb/Mf3/MfX/MfX/MfX/MfX/WdO7X1cVxJEmSJEmzt95+4yhJkiRJmjELxzmLiLdGxDURsTkizoyIHcvFfP4mIrZExNnlwj5NbX+vTHN9RLx03rGvdS25P6Pkc3NEnBIRD29pe39EXFEeXnBpQi25PzUibqrldb+WtsdExLfL45h5x77WteT+4lreb42Iv25p63o/hYj4LyXv10TEW8qwXSPiwrI+XxgRu7S0db2fQkvu3xUR10XEVRHxqYjYuaXtzRFxdVnvL59v5GtfS+4XI2JrbXtySEvbTeU7eUtEHD/fyNeHlvyfXcv9zRFxRUtb1/0JlP3G70bE5tqwxm18VN5f1u2rIuJZLX3uX96DLWX6mNfyjCUzfczpAewJ3AQ8srz+GHBseT66DDsZ+O2GtvsCVwI7AI8HbgA29L1Ma+XRkftDgCiPM5tyX6b/Ud/LsFYfHbk/FThiRNtdgRvL8y7l7136Xqa18mjL/dA0nwBe29Le9X75uX86sBnYiep6Ap8H9gHeCRxfpjkeeEdDW9f7lcn9QcDGMs07mnJfxt0MPKbv5ViLj47cLwL/bUTbDWXf5gnAI8o+z759L9NaerTlf2ia9wAntrR33Z8s3y8CngVsrg1r3MaX/c3PlP3N5wF/09Ln35bxUaY/uO/lrD884jh/G4FHRsRGqg/2bcC/BT5exv85cHhDu8OAszLz7sy8CdgCPGcO8a4nw7m/NTPPz4Lqw7pXrxGuXw/J/ZjtXgpcmJl3ZuY/AhcCm1YoxvWqNfcR8Wiq7U/jEUdN5alUOwY/ycz7gC8Bv0G1Lf/zMk3b9t71fjqNuc/MC8prgEtxe78S2tb7cTwH2JKZN2bmPcBZVJ8Xja8z/+Xo1ZFU/yjXlDLzy8CdQ4PbtvGHAaeVXc5LgZ0jYo96w/L60Zl5adkvPY3m74jeWDjOUWZuBd4N/D1VwfgD4OvAXbUvs1uojhIM2xP4Tu1123Rq0JT7zLxgML6covoa4LMtXewYEZdHxKURsao+xKvdiNz/QTll470RsUNDc57TyP0AAAyLSURBVNf7KYxa76m+kC7KzB+2dOF6v3ybgRdGxG4RsRPVf5v3Bh6bmbeVaW4HHtvQ1vV+Om25r3sd1X/zmyRwQUR8PSLeuIJxrkdduX9z2d6f0nKKtuv99Eat+y8E7sjMb7e0d92fXts2fpz1e88yvGuaXlk4zlHZUB5GdarpLwI/h/9Fnoum3EfEb9YmOQn4cmZe3NLF4zLzAOBVwPsi4okrGvA60pH73wOeAjyb6pS843oLcp0aY71/Jd3/eXa9X6bMvJbqdMgLqP4hdQVw/9A0SbWjphkalfuIOAG4DzijpYt/nZnPAg4G3hQRL1rZiNePjtz/KfBEYD+qf2K9p68Y17Mxtjujtvmu+zO0HrfxFo7z9RLgpsz8h8y8F/gkcCDV4erBPTX3ArY2tN3Ktv81aptOzZpy/wKAiPh9YHfgv7Y1LkduyMwbgSXgmSsd8DrSmPvMvK2csnE38FGaT712vZ9O13r/GKqcn9fW2PV+Opn5kczcPzNfBPwj8C3gjsHpSeX5uw1NXe+n1JJ7IuJY4NeAV5eduqa2g/X+u8Cn8GchE2nKfWbekZn3Z+YDwIdxe79iOtb9jVSnrZ7d0dZ1f3pt2/hx1u+tbHsK/ar7DFg4ztffA8+LiJ3KeeYvBr4JfBE4okxzDPDphrbnAEdHxA4R8XjgSVS/ydN4mnJ/bUS8ger3RK8sX2gPERG7DE6jLDvbB1K9bxpPW+4HG9agOmVyc0PbzwEHlfdgF6qLW3xuTnGvB425L+OOAP5vZv60qaHr/fQi4ufL8y9R7bD9JdW2fHCV1Lbtvev9lJpyHxGbgN8FXp6ZP2lp93MR8ajB31S5b9o2qUVL7uu/5fp1mnN6GfCkqK40/wjgaKrPiybQst2B6h+J12XmLS3tXPdno20bfw7w2nJ11edR/XTktnrD8vqHEfG88p39Wpq/I/ozq6vs+BjvAfwv4DqqD+PpVFdJfQJVEbgF+CtghzLty4G31dqeQHXFsetZZVdZWguPltzfV3J6RXmcWKY9APiz8vcLgKuprvB2NfD6vpdlrT1acv+Fks/NwF8A/2I49+X168pnYwvwW30vy1p7NOW+DF8CNg1N63o/29xfTFVsXwm8uAzbDbgI+DbVFQ93Hc59ee16P/vcb6H6jdFge39yGf6LwPnl7yeUNlcC1wAn9L0sa+3RkvvTy3bkKqod6D2Gc19eH0J1hOwGcz+7/JfhpwL/YWha1/3pcn0m1anX91L9HvH1Hdv4AP6krNtXAwfU+rmi9vcB5fv6BuADQPS9nPVHlCAlSZIkSWrkqaqSJEmSpE4WjpIkSZKkThaOkiRJkqROFo6SJEmSpE4WjpIkSZKkThaOkiRJkqROFo6SJEmSpE4WjpIkac2IiL0j4qKIuDYiromId0ZE9B2XJK13Fo6SJGktuQ84LjOfCjwTeC7wG/2GJEnrn4WjJKlVRCxGxPf6jmNYRJwaEZevUN9HRsSxK9F336JyRUQc0zDu8Ii4ICK+HxH3RMTWiPh4RGyacB7nRsTVHeM/EBF3RcQOtWFHRsTt4xw5zMzbMvPy8vc9wFXA3g3z+MgkcUuSulk4SpLWorcDx65Q30euYN99OxLYFfjL+sCIeC/wCWAr8AbgJcDxwCOBz0TEEyeYx5nA0yNi3+EREbEBOAL4ZGbeXRt1KHB+ZuYE8yEidgMOBz43NOrdwKsjYp9J+pMktdvYdwCSJI2jFB0bMvOezLyh73jWqP8MnJ6Z9w4GRMRhwFuA38rMU4emPz0iXgb88wTz+DTwE+CVwP8cGvcrwGOpisvB/B8GHAz89gTzoByx/Djwvsy8tj4uM2+OiEtKn78zSb+SpGYecZQkjW1wimhE/GpEXBURP46ISyLiaWX8seU0x52H2j0tIjIiXlJePz8izomI20ofV0TEq1vmdXhEXAP8lOr3bA85VXXC/hpjH0wDvAL4NyXejIjFEbk4NCK+GRE/iYjzImLXiNgnIr5Y5nF5RPzyUNtx4n1aRHw2Iu4s01wbEW8ad3xDvPsAL6AqtureAlzWUDQCkJnnZuatQ329MCK+VJb5+xHx4Yh4VJn+x8C5wFEN3R0NfBf4Qm3Ys4GdgQtr/W+IiLeW9+mnJU9nDE5lLf9EOAP4Rma+p2WRP0F11NF9HUmaAY84SpIm9UvAu4A/oDoS9W7g7Ih4BvDXwAeBXwc+WmtzFHAH8MXy+nHAV4CTqQrCA4GPRsQDmXlmrd2/BN4JvA24HbipJaZx+2uNvZwm+fYyzc7AfyxtbhmRi7cB/wPYCfhj4EMl7g+X2P83cFZEPK12KuY48Z4LXAv8JnA38GTg0bV5jxo/7MXAj4ErBwMiYiPw/JKHsUTEgcDnqd7rI4DdgD8EdimvoTqieFRE7J+ZXy/tHk51EZszMvP+WpeHAhdn5g/LdA+jKm4XqHL3/6h+w3hwLX8fBP6J7qOJX6U6uvmM+jJLkpbHwlGSNKldgQMz89vw4I7+p4AnZ+Z1EfFZqkJxuHD8+KBgyMyzBiPKUaQvA3sB/57aaYxURclLMvOKroAm6K81duC6zLwhIu4EHpaZl46Zi+cPTp0tRxb/O3BMZp5Wi+c84ClUhd7IeCPiMcDjgcMyc3ChmYtqbTrHt9gfuDYzH6gN2w3YAfhOfcIS04baoPtrRdsfAl/NzKNq028FLoqIp2fmZuAzwF1URxi/XiZ7KVVxWX8/oCocz6i9fivwq8Czh05B/WiZ14HA64HNwDfKQchTMvP9Q/1eA9wPPAcLR0mamqdvSJImdfOg8Cq+WZ73Ks9nAy+O6sIlRMR+wL8qwynDdomI90fE3wH3lscby3R1W0cVjRP2Nyr2Sd089HvLLeX5Cw3D9pwg3jupirmTI+KoiPj5ofmOGt/kF4C2K+QOX5Tmd2px3Qu8qcS9E9URyo9FxMbBA7ikTLc/PHi1008CRw5OL6X658HfAV+r5WEPqltqnFdePww4DvjA8O8WHww08yuZGZn5jMzcrzyGi0Yy8z6q4vUXOnIiSRqThaMkaVJ3Db2+pzzvWJ7PoSoiXlFeH0V1uucltTanluHvAg6i+p3bKbU+Bu4YM6Zx+xsV+6Ta+rurYVh9HqfSEW85KngQ1em5pwC3R8TFEfHMcca32JHqlNa675dhw4Xz6SWmZw8N34XqSORJbFtY3g08nG1vi3Em1am8z4+IHYHDgLOGrpx6CHBjZl5fXv8ysDvVUeBZuJvlv7eSpBpPVZUkzVRm/igizqMqjD5EdQuIvxoUDKWI+DXgTZl58qBdy0VMRt6eYcL+ejduvJl5HfCK8tvAFwLvAM6LiL0y84FR4xtmfSdDR98y876I+BpVEXpibfgdlKI9tr214l1U78kicH7DPOoX0fli6eNoYA/gUTSfpnpe7fUe5fn2hr6XY2eq5ZYkTWlVfqlKkta8s6iuTPoy4Anl9cAOVN8/Dx79KlfkfPky5zXr/u5hZY9STRRvZt6bmV8A/g9VYbXzJONrrqf6XeSw9wHPjYjXjAq8XDH1Uqrfs17e8Li1Nu39wMeAfwe8iur3lfUL8zyC6reM9cJxUDA+ZVQso0TE7lQXLPrWtH1JkjziKElaGedT3cvvg8BNmfm3gxGZ+YOIuAw4MSJ+CDxAdbP5H9B9VdBGs+4PuA44LCIOpzrF9tbh21FMY5x4y0V23k31u9AbqU4RPQ64MjPvHDW+ZdZfKfPcPTP/oRbPpyPifcCpEfErVFdr/R7VhXMOKpP9qNbP71JdCOcBqquf/hPVKamHAidkZr1QOxP4T1RX2f39oXheBATwpdqwq6kuavOnEXEi1RHMfYG9M/O4luVqcwDV0dGvTthOktTAI46SpJnLzH+m+q3jHtQuilPzKqqC5zTgj6juuXfaFLOcZX8nARdQ/XbwMqqL1szaqHhvpzrN8wSqK5SeRHVF1pePOb7JEtVpm5uGR2TmW6lupbE38BGqi/ucRHU7i0Pq93jMzEuoir7dqX4LeS5VMfkdhn6TmplfA26mKhCbTlP9fGbeXZv+PuBlwDeA95a+31BeT2oT8KXM/P4y2kqShsS2v1GXJEnrVUT8EbBPZh66CmL5FvCuzPzwCvS9geoKrsdn5l/Mun9J2h5ZOEqStJ2IiL2ofvO339AppetKRBwNvB14ajmKKUmakqeqSpK0ncjMW4DX8bOrl65XAbzeolGSZscjjpIkSZKkTh5xlCRJkiR1snCUJEmSJHWycJQkSZIkdbJwlCRJkiR1snCUJEmSJHWycJQkSZIkdbJwlCRJkiR1snCUJEmSJHX6/5BOFSFsUtXyAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1080x720 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "threshold = 30\n",
    "\n",
    "highE = bump[(bump.pt1 >= threshold) & (bump.pt2 >= threshold)]\n",
    "lowE = bump[(bump.pt1 < threshold) & (bump.pt2 < threshold)]\n",
    "\n",
    "fig = plt.figure(figsize=(15, 10))\n",
    "\n",
    "plt.hist(highE.M, 300, range = (80,100), alpha = 0.5 , color = 'black', label = 'High energy')\n",
    "plt.hist(lowE.M, 300, range = (80,100), alpha = 0.5 , color = 'red', label = 'Low energy')\n",
    "\n",
    "plt.xlabel('Invariant mass (GeV/$c^2$)', fontsize = 15)\n",
    "plt.ylabel('Number of events \\n', fontsize = 15)\n",
    "plt.title('Measured distribution of events with two muons \\n', fontsize = 15)\n",
    "plt.legend()\n",
    "\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.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}