MS3_MS2_Comparison.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "329c7b6d-de9f-4c68-9744-f9b54e09f2e6",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.24/00\n"
     ]
    }
   ],
   "source": [
    "import ROOT as R"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "c24d3d7c-648f-43ce-81b7-8a2e015a4f4a",
   "metadata": {},
   "outputs": [],
   "source": [
    "files={}\n",
    "files['MS3']=R.TFile(\"MS3/prodAnalysis_PMT.root\")\n",
    "files['MS2']=R.TFile(\"MS2/prodAnalysis_PMT.root\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "id": "5289818f-7c34-4210-a0c1-e758fbcf6f72",
   "metadata": {},
   "outputs": [],
   "source": [
    "g={}\n",
    "g['lyDT_MS3']=files['MS3'].Get(\"lyOverDT_byProd_geo2_PREIRR\")\n",
    "g['lyAbs_MS3']=files['MS3'].Get(\"lyAbs_byProd_geo2_PREIRR\")\n",
    "g['lyNorm_MS3']=files['MS3'].Get(\"lyNorm_byProd_geo2_PREIRR\")\n",
    "g['tau_MS3']=files['MS3'].Get(\"dt_byProd_geo2_PREIRR\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "id": "5435bbaf-5f01-4250-baf3-22f792894a14",
   "metadata": {},
   "outputs": [],
   "source": [
    "g['lyDT_MS2']=R.TGraphErrors()\n",
    "g['lyAbs_MS2']=R.TGraphErrors()\n",
    "g['lyNorm_MS2']=R.TGraphErrors()\n",
    "g['tau_MS2']=R.TGraphErrors()\n",
    "for i,vendor in enumerate([1,2,3,4,5,6,8,9,10]):\n",
    "    prod='prod%d'%vendor\n",
    "    if vendor in [4,5,6,10]:\n",
    "        prod='prod%d_ms2'%vendor\n",
    "    g['lyAbs_MS2'].SetPoint(i,i,files['MS2'].Get(\"lyAbsAvgQECorr_%s\"%prod).GetMean())\n",
    "    g['lyAbs_MS2'].SetPointError(i,0,files['MS2'].Get(\"lyAbsAvgQECorr_%s\"%prod).GetRMS())\n",
    "    g['lyNorm_MS2'].SetPoint(i,i,files['MS2'].Get(\"lyNorm_%s\"%prod).GetMean())\n",
    "    g['lyNorm_MS2'].SetPointError(i,0,files['MS2'].Get(\"lyNorm_%s\"%prod).GetRMS())\n",
    "    g['tau_MS2'].SetPoint(i,i,files['MS2'].Get(\"dt_%s\"%prod).GetMean())\n",
    "    g['tau_MS2'].SetPointError(i,0,files['MS2'].Get(\"dt_%s\"%prod).GetRMS())\n",
    "    g['lyDT_MS2'].SetPoint(i,i,files['MS2'].Get(\"lyAbsOverDt_%s\"%prod).GetMean())\n",
    "    g['lyDT_MS2'].SetPointError(i,0,files['MS2'].Get(\"lyAbsOverDt_%s\"%prod).GetRMS())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "id": "dc91f34f-bdc3-4db3-8831-6e00b132ca62",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning in <TCanvas::Constructor>: Deleting canvas with same name: c1\n"
     ]
    }
   ],
   "source": [
    "c1=R.TCanvas(\"c1\",\"c1\",800,600)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "id": "28b2d471-42f1-4f5d-9f97-3f9427ec7ab3",
   "metadata": {},
   "outputs": [],
   "source": [
    "l=R.TLegend(0.8,0.8,0.9,0.88)\n",
    "l.SetFillColorAlpha(0,0)\n",
    "l.SetBorderSize(0)\n",
    "l.SetTextSize(0.04)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "id": "3f0a32ce-2be4-4355-9f6d-2fa7dcd59e99",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dMXLbyKI2UEDlZFzlkjSxPWXtwF4BQe/lyqrr0PMH9gtIRvOCUTjjuvLbyhPBHXgFtgPGT1akbPgHPcbAIAVRYhMkgHMCFwmQUrcAEp+7G93pYrFIAADYzMGuCwAA0AVCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABI2GqjzPsyzLsmw8Hld2pUsqrxmPx7e99869AADb9qix3zQejyeTSZIkg8FgMplMJpPFYhF25Xle/94sy2azWXg8m83yPC+/pX4vAEADGmqpyvN8MpmMRqPFYpHn+XQ6TZIky7Lya6bT6aKkaHMaj8ez2Sy8d7FYjEajkJzW2QsA0Iy0aC7aqtCYVP5dITMV/5YbrqpFTNMkScp70zQdDAYhOdXvBQBoRkPdf7PZbDAYlLeUBz8VASg8qLRgJUlSee9gMCj6++7cCwDQgOYGqodR5FmWpWmaZdlyS1KapsPhcDgchhdU3vvgpwAATVhsXxhBFQwGg9FoFB4Xg6jC02JcVGh5Ck/De4tdxcCpouT1eyua/dMCQN9FCxNt0Nzdf0lp5NN4PA7tUmHLdDotNy+FmRcmk0lo2dpSGQCArQrjnvujie6/EIwqI5+K9qpkVYdd2HLbYPP6QeiGqAMAzWt0TNWD3yJFAQB7rqFQFSb8LG8p3/G3PH96OSct381Xvpewfi8AQDMaClUhMxUtT8WMnWFjiFzL83mufG94UJ4atGYvAEAzVkz+mX+XfM8oYVm9DX9TsUxNUJmfszKWbTQaVSayGg6HxdPlge01e8vStKHJTgGAvl12/6ltnuehiei2l1aCzsPcNr1n8j3M1QS4mvfeuTfo29EFgB3q22X379qGZWQGg8HKWQyKtqvQK9fezrW+HV0A2KG+XXbTxWKRpun6Ual+soM917ejCwA71LfLbs9q27OjCwA71LfL7kHiXjkAaETTy6bs2q7/3k37u/svSZLbBlR1Sd8iMwB7ZfNr0J13ZRXjc1aOkC5uRxsMBg1c8Q8ODvp12V0sFqPRqDxb5mAwKJY67pikf/9LAGB//LWZYgbH0Wh022uKq3l54+Xl5coMcHl5WXn75eVleXrt5RfcS98uu393/+V5vviermaz2XA4TNM0y7KWDkgHgA6rLFJSWHnVzvP81atXSZKMRqMQki4vL0M+e/XqVfkt4ZVFa9ZsNnv16lV5JkjusDJqTafTTrZd3VZfAGjAJq0+f5VaqpJVjUx//fVX+dpd2bj8+vDTBoNBuWGp8sqa37WOvl12Vy9TE9qoFt/TVWi7ihvmAIAHCGFo5WTdK1e/Da9cHkFV+Tmhyaoy1ir0G9ZMDE5Z3dp/d86xDgDsxHIPYNhSbs2601+l4Vbhcr/y7cYCrWlFqAprxaRpOhwOQ+YNfX/NF24b0lq7Lh0A/XVwUNfSURaiTyXr3HZjYGi7Ojg4WDkSq3h9GPxeeXt4S7dnBojon+O3nKVGo9FisQjbd1fCyOp7Q3ddOgD6Zf0gVRZCVSUkrez7S5JkOp2GB5PJ5ODgIKSr+sanPM8nk0k5DzygkD30KEmS8XhcPjCtXt0PAFrkr7/+Ojg4+Ks0FULlaY3y+Jz6vr+//vor5KTwluKiPxqNVr4l3C1YvGadwpCElqriSIRuPokKABoTclV4vH6iqvQA3jkpaJZl0+n0r9J8CkmShOaolUUqT7tw2wwOVKTd6+CrYUZ1AHao5hoUclV9oppMJpPJpGheOjg4COOelx/f+aOS0vxVl5eXNTFgzZ9223t7ddk9SG4JtnmeG+0PANtwsOS27fWDrsrdeSv76cI4qpUX9CzLyrMquO5vbsWhKoarhybBMLV60+UCgO5aniSzaApaOYXmSkUPYE3fXxi6vk7/3WQyKQ+l4gGqoSpN09lsVl4NMEz+KVcBwDaUx1GVx1fdqbgH8Lb7/pLvSWs2m61shSo3cYVXVuKXtqt7+eHIhSHq0+l0PB4XKSrP89FoZApQAIhueWT6fUcv1UzaGYTpPcMqfmEyhXAnYEhv5UWak+/DtorXFIOu7ler3irP0hRybngc/rjlCZw6sPxf0rNFiADYK3eulFf/gnBpDrN0lrdU3ri85bZUVP5Rt73swQv//dW/tf9+uBsuy7LZbBa2hMmrir1pmk6n07Z3Arr7D4Ad2u01KM/zot+pssbfypfVvGZNfbv774eQkef5cDgM4akcqsphq9WEKgB2qG/XoL6FqkflJ1mWDQaD4XBYjHcrJls3oSoAQI0VLTeVVWuSJOlAx1+gpQqAHerbNahvLVV1IaN7M60LVQDsUN+uQUJVlwlVAOxQ365BfQtVj5La9RfLzAAGAHCbdLFYpGm6zks7EDa1VAGwQ5tfg2pWpCm/YOVr1pxSIaK+tVT9EKoGg0F5LvXuuTM+9urYA9CwDa8yYbrzJElGo9Ftd+UXq9z8VZqZPc/zlev6XV5eLl/0i98Sgtcmt//3LVQdJKUp1Gez2XA4DCsod7Wzr34u1F2XDgDudtsCySuv3UWiGo1GYXr0y8vLEJVevXpVfkue5wcHB8UPn81mxWo2rOOfv9R4PO5VugKA9qpZIHnlxhCkQrtUaH8qFvsrXllkr2KdmTBv5XA43EYVumdF/LwtXTVeNgCgKoShYnRU2Ww2K6bvLm9MVg2xqvycIqWV+/um0+ltv4tldW16IV3VHDwAIJZ7XWqXG6UesAJKeRHl8NuX3x5Smm6rddwaqsLMn2maWqYGALZnPp//13/9188//zwcDn/++ef379/P5/P6t4SLciXo3HZjYEhF5cFSZcXrw4ir5cv9bQ1dLKuGqiJLDYfD2Ww2Go1CV+B4PN5F8QCgy+bz+dnZ2dXV1f/+7//+3//93+Xl5fX19evXr+tz1fJwqOSWvr/kexdeeP3BwUFIVytbnpaTUxhNtfLHsuzvULWcpabTqSwFAFv14cOHX3755ffff3/58uXR0dGLFy9+//3358+f//nnn3e+t9xdWN+tFLr5imw0mUxevXp1W9tVIUSCpBTLqHeQJMnKLKWhDwC27cOHD//6178eP35cbPnpp5/+9a9/ffjwof6NlR7AOycFzbJsOp2W51NIkmQymay8sy+0aYWmr/J8V9QzozoANKS4Bt1r8qci1oRpOYvZEA4ODgaDQWhGqjxOfpz8c6Vi/qryFKDlaUJXTg16L32b/PNRoq8UAJpVJJ6ff/758vLyxYsX5b2fPn0aDodXV1d3/pzQPVfT9xdy2Mp4FGarmkwms9ks7C1PE+oGtQd4lLhPEgB25Ozs7OPHj7///vtPP/0Uttzc3PzP//zP2dnZne8NkSjP85q+v9BuMplM1mlzWm644l7+HlMVulp3XRgA6Jd///vfX79+/fXXXz99+nR9ff3p06f/9//+3+fPn//973/f+d7iHsDb7vtLviet2WxWMwN7eXhWMes6D/D32n8h7YZ05Y4/AGjG06dP//Of/xweHg6Hw+Pj4+Fw+OTJk//85z/Pnj1b8yfcNmlnIUzv+erVq+FwWLRsFYv6lQetJ6VpFyo2rGZfVFYULm65HAwGxSRVnbFcXwBozF+1wmrHNUIAKq/NV0Si8suWtxTTplcsL/N3m/qC3aZvl91b74abTqdheerwR+9G85W7/wDYod1eg/I8L6a2GgwGDXTz9e3uv7tDRpGuLi8v275OtVAFwA717RokVP0jTLPeYGG2TqgCYIf6dg3qW6haMfQsrFcTpllPvt8b2HS5AABapRqq0jQN69UUA9YGg0ExLVgHpLV2XToAoK1+CFVhNPp0Oi2v/Zfn+Wg0Kq/a2Gr14/Z3XToAoK1+CFW3TckawpaJ1wFgE/W9Jd2z679300znBQAQwYruv+UWqdB21ZlhVQAA0T0qP8mybDAYDIfDYpT6eDyuWfsaAIBgxbxNRZAqmFEdALivvl12Tf4JAGxF3y67j1ZuDUtYF087Fq0AAKKrRsg8z5cX+BsMBt2YT6FvkRkAdqhvl93qlAohUY1Go+l0ulgsptNpx2ZUh37q+eQxAA24Y0b1LMuKGdW70VgF/VSsGWDxAIAtuceM6gAA3OaHUHVbH582KgCAetURZGmaLs9KlaZpN8aq923EHFT4CABN6tt3zg+1zfN8PB7PZrMkSYpJ1StPkyTJsqylHYJ9O7pQ4SMANKlv3zk/1HZ5LvXbtPRv1LejCxU+AkCT+vad07Pa9uzoQoWPANCkvn3nrJ5RvcPqZ+jp1bEHACI6SJIkTdP1x0gVU1i11KLWrksHALTVQZIki8ViMpmkaRqm+lz5upCl0jStLAsIAEBSdP8tFotw69/ywn+FwWCgLQcAYKUVI8jy78LT7LuGS7YNfRsxBxU+AkCT+vad07Pa9uzoQoWPANCkvn3nHNz9EgAA7iJUAQBEIFQBAEQgVAEARCBUAQBEsGKZmvJkCk0WBQCgvf5pqcrzPMyZPvwuTdMwhXqsXxZ+RZZlK1fFCZO2b2MvAMC2/T2BxHg8nkwmSZIMBoOigSrP89lsliTJaDTaPKyUf0X4seW5K7IsCxuDwWBQDnOb7C3r24QZUOEjADSpd985i8ViOp0m31ehWTYYDJIkmU6n9UsR1wu/YjQalZ8Wv3E0GpX3hqfFb9xkb0Vy14LK0G0+AkCT+vadky4WizRNkx/bjSrufMGdQmNS+SeEpq/w7/LPT9O0aHDaZO9yRTapBbSdjwDQpL595/w9pio08Nwm7N1kcNVsNgstXoXxeFzuUqzsLboIN98LANCAgxCV6m/0izL6O4wiD2PhsyyrRLT6AlT23uspAEAD1pqnasMbAMPbJ5NJGKg+Go1ms9lwOAzbV/7w8mD5TX71svSh4hYDAOiYRyG+hMkObnvROq1Z6yg6VsfjcZi7YbFY3KuNanO96twFABrzd0tVaEO6Tf3eO4VgVBn5VD+Kq76BapO9AADbcJAkSZjg4LZ5PkMkCq/ZRH2bkxQFALRaQ5N/Lk+pUN6yvDdN0+KXrtxbTJpQv7da257d2wkVPgJAk3r3nVOetKrSQ5ckyWAw2HDaz6B+ts/K3lCM295bmYy0fm9F0rNZyKDCRwBoUt++c1ZHyPpx6w9TNIYVca3cmJTn+XA4LJ5Op9NyATbZW9a7yAw/8hEAmtS375y1apvneZ7nUWarqr+RcHt7g74dXajwEQCa1LfvnLVqGxqZOvB36dvRhQofAaBJffvOWWvyTwAA6glVAAARCFUAABEIVQAAEQhVAAARPEqSJE3TXRcDAKDdHiVLSx0DAHBf/ZpAom8TZkCFjwDQpL595xhTBQAQgVAFABDBo10XoGn1o/J71UoJAETUu1AlNgEA26D7DwAgAqEKACACoQoAIILejakCgA1V7nkyWpdAqAKA+wkpqm8zW3In3X8AABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEfRu7b/KKpgVVnGieyxPBtCM3oUqVxf6ZrFYVHKVmAWwDbr/oPtCrgqPJSqALRGqoBeKXCVRAWxJ77r/oFeWBxGWtwhYABEJVdBlRWwKvX4aqwC2R/cfdF95HFV5fBUAEWmpgqpK5mh7u87yyPS21whgPwlVUBUyR2fukutGLQD2n+4/AIAItFRF1rGeIwBgTUJVZB3rOQIA1qT7DwAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIILeTalQv+qZeRAAgIfpXagSmwCAbdD9BwAQgVAFABCBUAUAEIFQBf+ov48BAGoIVfCPxWJRyVViFgBrEqrgB+Vclaap20UBWJNQBVVFrpKoAFhf7+apgtss9/SVtwhYANQTquBvRWwKvX4aqwC4F91/8IPyOKrlcesAcBuhCv6xPDJdSxUAaxKq4B8iFAAPJlQBAEQgVAHAuoyzpIZQBQDrsu4CNYQqALgH6y5wm97NU1X/XwqfDQDuZCo7VupdqPIBAOBhrLtAvd6FKgB4GOsuUM+YKgC4B+sucBuhCgDWZd0FaghVALAuEYoaQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEZlSHqvJSqYk7qAFYT3MtVemS8Xi85t4kScbjcZZlWZZVtq+zF+5l8aNdFweAdmiopSrP8wfvTZIky7LZbBYez2azPM/Lb6nfCwDQgEbHVFUaACqtStPpdOXe8Xg8m81Go1HYPhqNQnJaZy8AQDOqaxhtyXg8nkwmt/2u+r3L41rSNB0MBiE51e9d/lHN1LexXwT34syEiHyg7tS3P1FDLVVFxFnZN1e/N0mSwWBQeVr09925FwCgAY3e/VfcVJUkyXJjUs3eLMvKrywPorpzLwBAAxpqqQoppxj5FBqTilFTNXujj45avs1wTXGLAQB0TEOhKgxCL1JUnueDwWAymdy5t9IKtbnFQ8UtBgDQMQ2FquVsFLaEhqj6vcs2nKABACC6HS9TU98QVeyVogCAPddEqMrzfHmG9PIdfzV7k1V3881ms+KOv/q9AADNaCJUhQanyWSyPGNn2BtGUK3cG54mpVar8KA8NWjNXgCAZjQ3K1flBrrRaFRZ+69mb57nw+GweDqdTsv9hvV7K2Uw+Sd95syEiHyg7tS3P1GjtQ1ze4aVj++7N7l9VPs6ewOhip5zZkJEPlB36tufqGe1FaroN2cmROQDdae+/YkanVGdlqp0zvbqEwIAaxKquFtIUX37DwcA3MuO56kCAOgGoQoAIAKhCgAgAqEKACACA9WhF4pbOMMD9xzsJ3faskNOv80JVdALvh9bwZ227JDTb3O6/wAAIuhdS1WlebNCPAcAHqZ3oUpsAgC2QfcfAEAEQhUAQARCFQD0Wv1oY9YnVAFAry0Wi0quErMepncD1YHOMFchxBJylamqNiRUAW3lAgD3VfNfkaK9ygfqwYQqAOiL5f+KLPf0lbcIWPciVAFAf5XTlcaqDRmoDgB9V267Wh63zpqEKgDoteWBiVqqHkaoAoBeE6FiEaoAACIQqgAAInD3HwDcTzGO271ylAlVAHA/UhQr6f4D2sf93sAe6l1LVf13sf98QCuU1ykLLFYD7FzvWqoWtXZdOmBd5fkJJSpYx3w+f//+/fHxcZIkx8fH7969m8/nxd40TcNnqnjAffUuVAGdYUkNWN98Pn/9+vXV1dXl5eXV1dV0Or2+vj49PS1ylVaGzQlV3KH+fzbQvLTkti1AxZ9//vnLL7+cn5+/fPny6OjoxYsX5+fnz58//+OPP3ZdtO7oV5t5Y30EnemMCP+z+eWXX05PT09OTr5+/XpxcRH+ffr06a5LR991df3XznyBsFeOj48vLy9fvnxZ3vjp06csy759+7alX9q3k7lntRWq7un9+/dXV1fn5+ePHz8OW25ubt6+fXt4ePjbb7/ttmz0XPEpK6JVNz50SYe+QNgf19fXR0dHV1dXR0dH62yPpW8nc89qK1Td007+ZwN3Kn/EOvNxK3SvRuyD4+Pj6XT64sWL8kYtVXEZU8Wtrq+vv337dnJyUtl+cnISdu2kVJB0q7MPmnF2dnZxcXFzc1Nsubm5+fjx49nZ2Q5L1TFCFbc6PDw8Ojr6+vVrZfuXL1/Crl0UCoCHePPmzdevX9++ffvp06fr6+tPnz79+uuvnz9/fvPmza6L1h1CFXX8zwagG54+fXpxcXF4eJhl2dHRUZZlT548ubi4ePbs2a6L1h396uw0puq+5vP56enp8+fPw91/X758+fjx4+fPn30O2R+d+bgVulcj9o2r4ZZoqaKO/9kAwJr6FSFl8010slJ0QPfOzO7ViH3jarglWqoAACIQqgAAIni06wIA0GWVBRl71RlE3whVAGxReTWhXZcFtkv3HwBABFqqgLYq+pXCAw0hwG71LlRVevcrfClDi/jAAnuld6HKtzAAsA3GVMU0n8/fv39/fHycJMnx8fG7d+/m8/muCwUAf0vTNPTYFA+ISKiKZj6fv379+urq6vLy8urqajqdXl9fn56eylUA7InFj3ZdnK7p1z2uW72n9/3791dXV+fn548fPw5bbm5u3r59e3h4+Ntvv23plzbJHdHQjE5+1jpZKe7Ut+Pes9pu8+geHx9fXl6+fPmyvPHTp09Zln379m1Lv7RJfftswK508rPWyUpxp74dd91/cVxfX3/79u3k5KSy/eTkJOzaSakAgMYIVXEcHh4eHR19/fq1sv3Lly9h1y4KBQA0R6iK5uzs7OLi4ubmpthyc3Pz8ePHs7OzHZYKAGhGvzo7t9q5O5/PT09Pnz9/fnp6enJy8uXLl48fP37+/Pni4uLZs2db+qVN6lvXOOxKJz9rnawUd+rbcddSFc3Tp08vLi4ODw+zLDs6Osqy7MmTJ51JVABAvX5FyMYicyezeScrBXuok5+1TlaKO/XtuGupAgCIQKgCAIhAqAIAiECoAgCI4NGuC9AvlSXBezV8DwC6TahqVEhRfbsbAgD6oHehqtJWVCHrAAAP07tQJTYBNEOrPH1joDoAW7FYLCqdA/V9BdB2QhUA21LOVRqu6Lzedf/xAOXvxEQXKnAfRa7y1UHnCVXczVchcF/LPX3lLb5V6CShCoD4itgUev00VtEHxlQBsC3lcVTL49ahY4QqALZieWS6liq6TagCYCtEKPpGqAIAiMBAdTZiiWgACJprqUqXjMfj8gvG43GWZVmWVbZvvpftWSwWIUgVDwCgnxpqqcrzvP4FWZbNZrPweDab5XlefssmewEAGtDomKrFj4pWpfF4PJvNRqNR2D4ajUI22nwvAEAzGlqJaTweTyaT237X8qRwaZoOBoOQjTbZu/yLmqlv/S/q3gJY3asR+6OH4/Y6+YHqZKW4U9+Oe9Pdf+FBlmWVFwwGg8rTokdvw71Aq4Vv5L59NQNt1Gj3X5qmw+FwOBymaVrJVcsxq2bvvZ4CADSgoVAVmo6KkU+hMSmMqVrZT1cEo+ijo5ZvQlxT3GIAAB3TUKiaTqflkel5ng8Gg8lkktyzjWpzi4eKWwwAoGMaClXL2Shsua0hqr6BapO9ALAO/RXc145nVF+zm0+KAqBhbpLgvppoqcrzfHn+9HISWr5fL0w9VbO3uOOvfu/+8L8cAOi2JkJVaI6aTCbLM3YWT5NSq1V4UJ4a9MF798disajkKjEL6Aw9ZZA0NvlnspQhRqNROfrkeT4cDoun0+m0PAxrk72VMux28s+wfXnC0rbTPM629eoca29la0reyUpxp7799RqtbViVL6x8fNsLktvv+Ntkb7DzUJWsmgK+A/r2sdlbHZ58vFfnWHsrK1RR0be/Xs9qu6NQVd8Y3oFD0LePzZ7r5OHoZKVu097KClVU9O2vt+O7/3qiOKW62v0HADS6TE3PlQP78rh1AKDVhKqGLDeBaqkCgC4RqhoiQgFAtxlTBcAWFUMdDCel84QqALZIiqI/dP8BAEQgVAEARCBUAQBEIFQBAETQu4HqnV8xBgDYid6FKrEJ2EN9WyINOkn3H8DuLS9dZSUraB2hCmAvlHOVhitoI6EKYF8UuUqi2rn5fP7+/fvj4+MkSY6Pj9+9ezefz3ddKPadUAXsr55c2NKS27bQpPl8/vr166urq8vLy6urq+l0en19fXp62snTj4iEKnoq/dGui8MK/bmwLb5LSm1U5Y007M8///zll1/Oz89fvnx5dHT04sWL8/Pz58+f//HHH7suGnutX932jQ1TqP9F3Rst0d4atbfkNTpTqffv319dXZ2fnz9+/Dhsubm5efv27eHh4W+//bbbsm1DceDCgzYexzaWeaXj4+PLy8uXL1+WN3769CnLsm/fvu2qVG3UmVNiTT2rrVC1He2tUXtLXqMzlerVha181Np7BNtb8rLr6+ujo6Orq6ujo6N1tlOjG6fE+nT/Afvo+vr627dvJycnle0nJydh105KtT29uvDsucPDw6Ojo69fv1a2f/nyJezaRaFoB6EK2EcubOzQ2dnZxcXFzc1NseXm5ubjx49nZ2c7LBX7T6gC9pQLG7vy5s2br1+/vn379tOnT9fX158+ffr1118/f/785s2bXReNvSZUAXvKhY1defr06cXFxeHhYZZlR0dHWZY9efLk4uLi2bNnuy4ae61fI8gMVN+S9taovSWv0aVKzefzP/7448OHD9fX14eHh2dnZ2/evOn8ha29R7C9Ja/RyUo1pm9/vZ7VVqjajvbWqL0lr6FSbdfeyra35DU6WanG9O2vp/sPACACoYoHMgs5AJQJVY0qVkTpwNIoxcqvhbbXCAA28WjXBWha/YV/212/HetaLq+k0beOcwCo6F2ocuGPq2iv8ocFoOd6F6qIYrnBr7xFwAKgh4QqHqKy8qvGKgAwUJ2HK4+jWh63DgC9IlTxQMsj07VUAdBnQhUPJEIBQJlQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQRR/N5/P3798fHx8nSXJ8fPzu3bv5fL7rQgHQbkIVvTOfz1+/fn11dXV5eXl1dTWdTq+vr09PT+UqADZRnRS725YnAW/7L9oT7arv+/fvr66uzs/PHz9+HLbc3Ny8ffv28PDwt99+223ZomjX4VhTJyt1m/ZWtr0lr9HJSjWmb389LVX0zocPH05PT4tElSTJTz/9dHp6+uHDhx2WanP6NCGuNE3DkqbFA6gnVNEv173FoXoAAA3HSURBVNfX3759Ozk5qWw/OTkJu3ZSqs3p04ToFj/adXFogd6FqrTWrkvH1h0eHh4dHX39+rWy/cuXL2HXLgoVwZ9//vnLL7+cn5+/fPny6OjoxYsX5+fnz58//+OPP3ZdNIC+6FdnpzFVW9Ku+r579+76+vr8/Pynn34KW25ubn799dcnT57893//927L9mDHx8eXl5cvX74sb/z06VOWZe1tfitr1zm2ofZWtr0lZ0v6dkr0rLZC1Xa0q77z+fz09PT58+enp6cnJydfvnz5+PHj58+fLy4unj17tuvSPcT19fXR0dHV1VWlpe227W3UrnNsQ+2tbHtLzpb07ZToXfcfPH369OLi4vDwMMuyo6OjLMuePHnS3kSVdLdPk7ZwkwQE/YqQWqq2pL31bW/JKzrZpxlUBjt243jVa9dpGW6S+OWXX0LT79evXy8uLsK/T58+3XXp2LF2ncyb61lthartaG9921vyiu71afZZu07Lzk/8xibadTJvrme1Faq2o731bW/Jl83n8z/++OPDhw/X19eHh4dnZ2dv3ryRqNqoXadl52+SYBPtOpk317PaClXb0d76trfkNTpZqV5p0RHsw00SbKJFJ3MUBqoD8EBukoAyoQqAhzs7O7u4uLi5uSm23NzcfPz48ezsbIelgp3oV7uc7r8taW9921vyGp2sVK+06wi6SYIa7TqZN6elCoCH697Eb/Bg/YqQWqq2pL31bW/Ja3SyUr3S3iPY3pKzJX07JbRUAQBEIFQBAEQgVAEARCBUAQBEIFQBAETwaNcFaFplxfuKXt2kAABE1LtQJTYBANug+w8AIILetVRtW9G9GB5oGAOAnhCqIpOiAKCfdP8BAEQgVAEARCBUAQBEIFQBAERgoDo95T5NAOISqugpKQqAuHT/AQBEoKWKjehEA4BgNy1VWZbleV7eki4Zj8flF4zH4yzLsiyrbF9nL9uz+NGuiwMAO7ODlqosy2azWZ7nWZaFLZWAddtbwuPw3vJb6vcCADSg6ZaqPM+LAFQxnU7LbR5Fm9N4PJ7NZqPRKGwfjUYhOa2zFwCgGWnDXTZpmg4GgxCDyrFpMpncVpLlwTrhh4TkVL93+UfpoqLznOdt194j2N6SsyV9OyUabam6Le4UW27ruRsMBpWn5eau+r0AAA1oLlSFdqmajrk0TYfD4XA4TNO0GG4VbPIUAKABDYWqPM8nk8l0Ol25NzQsFeOiQlPTnSHsYZZvM1xT3GIAAB3TUKgaDoeDweC2NqQwRL0YYpXn+WAwmEwmyRaanRYPFbcYAEDHNDGlQkhLlUmk8jwvJpdaTk7L0y6U1TdfufUPAGhec/NUhZanwmw2m81m9Q1Ra05kJUUBADvXRPffeDxe7koLI6jC1OrL86eXc9Ly3Xyz2ay4469+LwBAM3a/oHKWZWEE1fJ8nsXTpNRqFR6U57iq2QsA0Iy9WFA5NFYNh8NiS3lq0CzLptNpmGohbJlOp+UUVbMXAKAZezTVaZj5c+W49eIFye33A9bvDfo2tSv95Dxvu/YewfaWnC3p2ynRs9r27OjST87ztmvvEWxvydmSvp0Sux9TBQDQAUIVAEAEQhUAQARCFQBABEIVAEAEQhUAQAR7MfknEEUxBW540Ks7mQF2TqiC7pCiAHZI9x8AQAS9a6kq+kdW8h99AOBheheqxCYAYBt0/wEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARNC7eaoAiM66k5AIVQBsToqCRPcfAEAUQhUAQAS6/wAaVVnWXccZdIZQBdCokKLSNBWnoGOEKoB94R46aDWhCmBfSFHQagaqAwBE0LuWqsoQ0Qr/TQQAHqZ3oUpsAgC2QfcfAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhVAc+bz+fv374+Pj5MkOT4+fvfu3Xw+33WhgDiEKoCGzOfz169fX11dXV5eXl1dTafT6+vr09NTuQq6Ie3VAsNp2q/6Anvl/fv3V1dX5+fnjx8/Dltubm7evn17eHj422+/7bZssA19u+z2rLY9O7rAXjk+Pr68vHz58mV546dPn7Is+/bt265KBdvTt8uu7j+AJlxfX3/79u3k5KSy/eTkJOzaSamAiIQqgCYcHh4eHR19/fq1sv3Lly9h1y4KBcQkVAE05Ozs7OLi4ubmpthyc3Pz8ePHs7OzHZYKiKVfnZ1pmta/oFd/DaBh8/n89PT0+fPnp6enJycnX758+fjx4+fPny8uLp49e7br0kF8xlR13KLWrksHdNnTp08vLi4ODw+zLDs6Osqy7MmTJxIVdEa/ImTfIjOwt3wd0Qd9O89711IFALANQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQASPdl0AgH5J07T8YLFY7LQ4QDS9C1XF19lKvt2AbfM9A13Vu1Dl6wwA2AZjqgAAIhCqAAAiEKoAACIQqlqvfuh9S3WvUt2rUaJS7dG9SnWvRklHK9U3QhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEuwlVWZbleV7ZOB6PsyzLsmw8Hi+/ZZO9AADblja/Fl6WZbPZbDQalQNQ2Fg8HQwG5dS1yd6yNN1BfbdNpVqhezVKVKo9ulep7tUoUalOaLqlKs/zcgAKxuNxiFmLxWKxWIxGo9lsVgSjTfYCADSj6QiZpulgMKi0VIVpZMslCS8L2WiTvcu/vXuRWaVaoXs1SlSqPbpXqe7VKFGpTmi0paom7gwGg8rTSo/eg/cCADSguVAV2qVu65jLsqzmvZW993oKANCAR838mjzPJ5PJdDpduWt5YzH2PProqE6uWKlSrdC9GiUq1R7dq1T3apR0tFK90lCoGg6Hg8FgZRvSvdqoNtSrnl0AoElNhKrQ8VeZRCrP82JyqeW31DdQbbIXAGAbGmqpSpJkMpmUn85ms9lsViQqKQoAaLUmBqqPx+PFj5IkCTNLhVC1fL9emHMhPF65t7jjr34vAEAz9mLtv6J/MDwND4q+wk32AgA0o7nuvxpZlk2n0+FwWNz4UL5PcOXecoqq2QsA0Iz9muo0jI66LRJtsncfVAbm53m+z6VdR7lGRZOhSu2h5ZtCxuNx2xt0b1tGvb31Kh+mPM/DV0SXzr1OfqDaXqk7L0yhgu39WDVsL7r/CvXn5SZ790FlQP1wOGz7EPui/FmWhcfD4XDPj8KdikqlaVpUqu1fKJUzbTweV24caaOiUpPJpO0fpaC8pOlwOEy6de518gNVrlRLT8L6C1OoYJ7nJtBa14K1TafTxWIxGo3Cg2LLdDottoxGo2J15+I14S2DwSC8rBhHX7xrV6LUKPTVFnt3flJFr9RoNBoMBk0U/Xaxzr2wcU8++7EqtQ91CaLXKGzcfsHrRP9AVb4xdiJKpcLdVMWuXVUq1lm3fGEqf/UNBoPKT2Clffky2n/FpSh8kMLpNRgMwolYfBUWW8K7ihcvp6idh6qINarUq+malEQ/TIs9CFVxK5V8H7O4k7oUIlYqbCkHx52IVaMiSO22OkHcw1RcvLvxgap8M+zkM7XVC1M5SO38a7AthKp1TVe1x5S/HVaef5W37GGoKp5uXqNFabKMXYlbqel+NOpErNRgb1p3YlUq/JzwH/Fw8Wim/Mti1WhUmk2mSx+ool4dO/eKl+0wVJULtoh3Yap8GQpV69ivMVXtUhkaX57LNMuy2WxW6ase7P3sWZvUaDwep2k6Go32bbTEJpUK95YObllhaYceVqnK/CP75mGVyrJssViEcfeLxaIya91ubXLuLRaLPM+n0+m+jX57WKWK5V8Xi8VoNNq3k/DB595oNErTdK/GG0W8MO3/NWsPCVXR3Hn+7dXX/TrWr1EYqB6ubVsv1mbWrFQYm5l8r9qeH7v1KzWbzYprQDHGdj/19gNVDhz7Fj6WrX/uFf85GY/He37s1j/3iqmt9/arb8PPUfEVsc/fFXtFqLqfcGKtnL6hfGv3eDwu/je25+fi5jUKl+q9qmaUSoXbr5Z37UqUShVt1EmSLL4vabBDm1cqNJGWf9pubV6j8rv2oUZJpEoVF+/OVCpJkuLcG4/HO2zX2dKFaTweFy+bTCZ7Gxz3S9P9ja1VDK8pD3QY/HhDxPIfthjnEf7dwzFVm9eoPFRi5ydVxMNUflk3jlT5Z+78sx+xUsXT3R6piDUqb+zMubcnh2mxtUo1X5HF9i9MYXhWstPRiu0iVK1rzVtywujF5Y1bKtUmulejRexKrXxZ8xyp5Y0rXxmnZA/VvRotVGrvvyUa+HLYh0PWIkLVuqadu/ehezVaqFR7dK9S3avRQqX2Xpfq0g3GVN3DzgegRNe9GiUq1R7dq1T3apSo1N7rUl06YL/W/gMAaCktVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAAR/H8qsVoc2PiquAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "l.Clear()\n",
    "g['lyAbs_MS3'].Draw(\"AP*\")\n",
    "l.AddEntry(g['lyAbs_MS3'],\"MS3\",\"PL\")\n",
    "g['lyAbs_MS3'].GetHistogram().GetYaxis().SetLimits(4000,6500)\n",
    "g['lyAbs_MS3'].GetHistogram().GetYaxis().SetRangeUser(4000,6500)\n",
    "g['lyAbs_MS3'].GetHistogram().GetYaxis().SetTitle(\"LO (pe/MeV)\")\n",
    "g['lyAbs_MS2'].SetMarkerStyle(24)\n",
    "g['lyAbs_MS2'].SetMarkerSize(1.2)\n",
    "l.AddEntry(g['lyAbs_MS2'],\"MS2\",\"PL\")\n",
    "g['lyAbs_MS2'].Draw(\"PSAME\")\n",
    "l.Draw()\n",
    "c1.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "id": "3eb40df9-e19a-4cc1-84cc-5a70aa0f54ce",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dv3LbSII/cEC/ScZVLlK5vGW9wfgJCDnbe4+ldbOh9wIrIZmsJ1iFO6qVLrjHuIgEH+Cm/ASWa4vxUYqYmb+gb7A0SVG01CQI4PMJXDT4r5toEl91NxrpfD5PAAB4nqOyCwAAUAdCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABD+UXYC9StO07CIAQIPM5/Oyi7A/zQpVScP2LgCUqGl9GfFDVZ7neZ6HG1mWJUnS7/ejvwsAwEFJY/Xc5Hne7/fH4/FDD+j1eqWnqzSNVl8AYLOmHXbjTFTPsuzs7CxJkl6vN18xGo16vd5gMEjTtPRcBQCwCxEiZJqm2/dChQHBMD64f02LzABQoqYddhtW24btXQAoUdMOuztfp6qsTikAgH2KHKr6/X6apiFIhdtnZ2fFFgCAuorZL5fneZiuHl4zrE4xGo3CWYGH0AHYtH5IAChR0w67MXuqwlz18PGFrqnRaJRlWdiuswoAqLHIi392Op1wI0SocK5fccZfuAEAzdSobpsGihmqsiwbDAbhdnEj+TZgAQBP9ughtRgXWn1MnufFGt2dTsdxOb7VtTqfbDQaJUnS6XRCf1VYCDRsjPtGT3YgxQCgmb4+T6/XC4fUXq/30GOK4/vixuFwuDYDDIfDpacPh8Ni0KnT6aw+4Ls07bAbc05VlmW9Xm88Ho/H406nE6ZSFSutR3wjAGiyxeGgRWunL+d5/vbt2yRJer1eCEnD4TAcl9++fbv4lPDIojdrPB6/ffs2HMfZxs6n5R/UVKqmnYYAwEF55jFoMBgUcWo4HK4eXs/OzopIVPRahY2rjw+v1ul0ijGlo6OjpVde3fJdjo6OGnXY3fnin4eTqPYj/VbZxQGgbkInUxGeFoWRotWNybrD8dLrhC6rpblWYdxw7XuxKnKoCv1S6Tpx3+hghVHVxRsAEN3qCGDY8l3zbRanW4XktPbpFkXaUsxQFRb/DDF5VcQ3AoD6CWNt2wjRZynrPHRiYDgEHx0drZ2JVTw+TH5fHSJc+5qsFXOOUZZlB7Jy+kMe7TCLVXiTtwBYtXpoODo6KiY/Ld5eK8yC6vV6vV7v6OhocTpUeHrYEsLZ4kstxbVer7d5SYWw+EL4d+ldvkvT5lTtavHPg9WovQvAgfv69etSlno0WhUWpzptHvv7+vVrnueDwSA8peiyCvls9fHhbMHiMdsUhiTu8F/oqYr4ggBQeyFXhdvbJ6qlEcBHFwXNsmw0Gi2up5AkyWAwWLtiwtKyCw+t4MCSyKNUaZr2er2wQtUB2tuonOE/AFZtODSsjtmtWhz+SxbG+1ZvP/pSycL6VZtXTNjy1R56bqOOhpEnqnc6ncFg0OSz/wDgUUcrHtq+efb64nDe2nG6wWBwdHS09vS9sGR38SJ5njvL75kiz6lKqjCtCgDKtdTxE0b9vrdPqNfrDQaDIgyt7W0KB+XBYPDoGXxhxtXTeqQImjVKZfgPgBKtPTQU86iKaPVQslka/kt+H/VbOkdvKZxtWBV98ZGrL54sDBEa/ttG/J4qAGBLqxHqe+PLhkU7g+Fw+Pbt27dv34ZlFELf1Xg8Xho0DP1eYePSYx66HjNLIoSq75qc3u/3jdoCQPDM4baQhJLHzvsLuWo8Hi+dpL/ULxUetnSu35Mv/NdAEUapwkLqSZJ0Op1+v7/2ow9ZKvRPlpioDP8BUKJyDw3hQBxub1j8s3jY5gVCt9G04b9ox/48z/v9/oZ1qsqNU4FQBUCJmnZoEKqeK/9d+G/2u7jv8jRCFQAlatqhQaiqM6EKgBI17dDQtFAVc/HPvcmybHUkMc/z0CV2sOu5AwA1Vr1QFa4wuBSq+v3+2dlZsbCsBdwBgD2rWKhaPHNhcWNYr2w+n+d5HlY/O5BZXABAQ1Rs8c+zs7OwdOzixjDeV4z6FRczAoCa2XBFmsUHrH3Mlksq8GRVClVpmoZ1GZZG98LyV4tbTKsCoH6KFc+XFu1cFK4qk3y7rGhxtZklaxf2XFxXXT/Fd4k8/Nfv93cUaMLLPrTSVZifnmVZmqZrp7EDQG0sLXpeWHv4KxJVr9cbDodfv34dDochKr19+3bxKXmeHx0dFS8eLlMTrg/IVuZRJUnS6XQWt4TdFiY8PVmYJjUajYp3KV6wuH5keOsiUBcPXirefj6o6B8sADXw9XmWLimz+oDFcZuljauPD6/W6XSKLeGJvV5v6bmLj/kuTTsa7jxUzefzsEsivuzaULX0+LXvuLe927RmBMA2nhZNlmJQ8e/aEFPkqsWNi/9denxxV3HV5A2P+V5NOxrGn1M1Ho/DUF0Ypl17vt53KV5tcWAxXBWnWKt9aU5VcYFJAKiK1SnCG4TT3pe2JEnS6/XWTp9a6+vXr8XwXzhYr86gCueHhcUgt3zZxtrVRPXFK9WECy0/8wWXQlK41PbmHawFAHD4JpPJ1dXV1dXV3d1du90+Pz//+eefT05ONjwl9B0sHeYeOjEwpKKjo6O109uLx/d6vbWnBIaw5Xi6lbgdX8m64b/okm8naa0OLz404Bi9vg/Z2xsBUCGrA2T//Oc//+3f/u38/Px//ud//vd///e3337793//9z/+8Y///Oc/Vx+8OPCXrEx1KraE4/vqEF6hmLH+KHOqvkvkKf2j0eiZg31PsDjaGP47Ho+dAgrA4bu6uvrDH/7wt7/97c2bN+12+6effvrb3/72+vXrX3/99dHnLh5wi7G/tY/8+vXrcDgsxhYHg8Hbt28XT/Rbq7hUyeI5YWwQOVRlWTafz/v9/j4XNQiraIzH4zRN0zQdDAZRBhwBYNeurq7+9Kc/vXjxotjy448//ulPf7q6utr8xJCfiqPto4uCZlk2Go2+LqynkCTJYDA4OztbfXBYSSHM8Vrt6OIh6bxGl49+tEml6Z7qu7c3AqBCikPDdy3+VMSasCxnMTXq6Oio0+mEbqSl28m3i3+uVaxftbgE6OIyoWuXBv0uR0dHjToa1mpFr+JkQAA4ZMWso3a7/dtvvy1NRfrtt99ardZDM6IWheG5DWN/oc9p7fBRsVp6MYy4uEzo169fHVK/V61CFQBUy/n5+c3NzWw2K7bMZrP//M//PD8/f/S5xQjghoGaMI9qy2WGio4r85KfRqgCgNL8/PPPX758+ctf/vLp06f7+/tPnz79x3/8x+fPn3/++edHnxuiz2Aw2LDAVUhaYaGp1XsXu7jCA3q9ng6qJxOqAKA0Jycn//jHP1qt1tnZ2fHx8dnZ2cuXL//xj3+8evVqy1d4aNHOQlgq/e3bt2dnZ2F1qzzPi4v6LU5aT34fLlz1zGo2xK4W/wQAtnFycvLXv/71r3/96xPWrC6uILL5vL/hcPj27duwbvbS0430RdSsk9Sc/QdAico9NCxeOG7t4unRNe3sv2Yd+4UqAErUtEND00KVUVIAgAiEKgCACIQqAIAInP3H49I0XfxvowbIAWBLjQtVS/lgibiwVvhYzL4HeKbNxyCqrnGhSiwAAHbBnCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACISqmCaTycXFxfHxcZIkx8fHHz58mEwmZRcKANgHoSqayWTy7t276XQ6HA6n0+loNLq/v+92u3IVADRBs67mttOr111cXEyn08vLyxcvXoQts9ns/fv3rVbr48ePO3rTfXLtPwC+S9MOHA2r7S737vHx8XA4fPPmzeLGT58+ZVl2d3e3ozfdp6Z9NwB4pqYdOAz/xXF/f393d3d6erq0/fT0NNxVSqmiMFEMALYhVMXRarXa7faXL1+Wtt/e3oa7yihUBCaKAcCWmtUvt9N+yA8fPtzf319eXv74449hy2w2+8tf/vLy5ctffvllR2+6a7WfKAbA7jRt+K9htd3l3p1MJt1u9/Xr191u9/T09Pb29ubm5vPnz9fX169evdrRm+5a7SeKAbA7TQtVhv+iOTk5ub6+brVaWZa12+0sy16+fFnpRFXjiWIAEF2zImSappsfEOvTqE02Pz4+Ho1GP/300+JGPVUAbKM2R8MtNa6nar5R2aU7OOfn59fX17PZrNgym81ubm7Oz89LLBUAHKBmRci9RebaZPNaThQDYD9qczTcUuN6qvgu9ZsoBgA70qwIqafqOWpZKQB2p2kHDj1VAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAAR/FB2AaiANE0XbzTqkuMAsKXGhaoiH6wlLqzlYwGARzUuVMkHAMAumFMFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQQSVDVZZleZ5vuLff7++vNAAAVQxVWZaNx+OHQlW4d78lAgCoWqjK83xDZtp8LwDA7lQsVJ2dnXU6nafdCwCwO1UKVWmadjqdDQN/SZJsmGsFALA7P5RdgG2FuecPZaZ+vz8ej+fz+aOvk6bp0wqwzYsDAI1VjVCV5/lgMBiNRk+4d4lsBADsQjWG/8JkqTDA9733AgDsQQV6qsLA39LqU3me9/v9YsGqh+6VtACA/ahAqAoGg8Hif8fj8Xg8LjLT5nsBAHYtreIcozRNe73eQ8umb7g3TfdU3729EQAcrKYdDasxpwoA4MAJVQAAEVRmTtWizX2JjeppBAAOhJ4qAIAIhCoAgAgqOfx3yIrL4IQbxiIBoCGEqsikKABoJsN/AAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAETRuSYViHam1LIgAADxN40KV2AQA7ILhPwCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIfii7APuWpumGe+fz+d5KAgDUSeNCldgEAOyC4T8AgAiEKgCACIQqAIAIhCoAgAiEKgCACBp39h8ANNbSukLOiI9LqAKApggpKk1TcWoXDP8BAESgp4pn0ZMMAIFQxbPoSQaAwPAfAEAEQhUAQARCFQBABEIVADTCZDK5uLg4Pj5OkuT4+PjDhw+TyaTsQtWKUAUA9TeZTN69ezedTofD4XQ6HY1G9/f33W5XroqoWSdtOUltR3ywAAfu4uJiOp1eXl6+ePEibJnNZu/fv2+1Wh8/ftzRmzbt6NCw2jZs7+6NDxbgwB0fHw+Hwzdv3ixu/PTpU5Zld3d3O3rTph0dGlbbb1eqXNWoTyOipn1tAKrl/v6+3W5Pp9N2u73N9liadnRo3Jyq+UZllw4A4mu1Wu12+8uXL0vbb29vw11lFKqGGheqAKCBzs/Pr6+vZ7NZsWU2m93c3Jyfn5dYqpppVr9c0/oh98YHC3DgJpNJt9t9/fp1t9s9PT29vb29ubn5/Pnz9fX1q1evdvSmTTs66KkCgPo7OTm5vr5utVpZlrXb7SzLXr58udNE1UDNipBNi8x744MFqJC9/Wg37eigpwoAIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACCCH8ouAACwJ8U1cMONRq13sAdCFU/UtNVHAGrA7/ZOVXL4L8uyPM+XNvb7/SzLsizr9/sllKl55vN58RdPsPRfAGiU6oWqLMvG4/FSqErTdDAYhNuDwSBN09XURXSLuUrH1SFIv1V2cQCapWLDf3mej8fjpY1ZliVJMhqNwo0kSdI0PTs7c4zfgyJX+bQPQdgLAi5AKSrWU3V2dtbpdJY2jsfjTqdTJKokSXq93l6L1Tyr3SE6SABouCr1VKVp2ul08jxfOmwvJSr2oOgICZ0iOqvYnaXvu2YGHKzKhKow/XztTKnVjcX8qlVP7kfxU75qcZgpRCufEtEZ0wSqohrDf3meDwaD0Wi0zSNDbHrowfOnilKROs0jXj3IOeYB0GTV6KkKU6keHeMLJwYm305aPyh1+pu7BlUAgIgqEKrCwN/SAlR5nhcLU4X/np2dJUnS6/WsUwUA7F8FQlWwNE1qPB6Px+PFRBXmsJdSNgCASo5DpWm62CO1/alnBzLudiDFiKh+Naq0Wu6OWlaK6nJS6paa9s2tTE/VQ4reqdVRP+OAAOxCnSbIElHlQ1VhdRkFoQoA2JtmpewD+aviQIoRUf1qVGlrd0fVRyu0MQ6QZvmopn1E9empAjYwWgGwa9VY/BMA4MAJVQAAEQhVUBOTyeTi4uL4+DhJkuPj4w8fPkwmk7ILBdAgQhXUwWQyeffu3XQ6HQ6H0+l0NBrd3993u125CmBvmjVr9UBm6R5IMSKqX40q5+LiYjqdXl5evnjxImyZzWbv379vtVofP34sHlbdPVXdklNjmuWjmvYRNay25e3dxbeuXyOrX40q5/j4eDgcvnnzZnHjp0+fsiy7v7+vQdurbsmpMc3yUU37iAz/7cl8Pl9aKGjpv/Bk9/f3d3d3p6enS9tPT0/v7++n06m2B7AHQtX+LOaqpoV3dqrVarXb7S9fvixtv729DXdpewB7IFTtVXFsc1QjrvPz8+vr69lsVmyZzWY3Nzfn5+fhv9oeRKGjlw2a9Tfro1+GHX0am9+3BrtA50fpJpNJt9t9/fp1t9s9PT29vb29ubn5/Pnzf//3f294VoX2mjbG4ShaY7ihcW7QtA+ncT1V8412/abJwmFs129Ko5ycnFxfX7darSzL2u12lmUvX768vr7W9iA64+k8pFmtodzWX+M/bupUlxpY3R01aHtVLDP1Zjx9G0375rqg8p6sNqxGtTNKpO1BLKtzORa3+GYhVO2JLxtl0fYglqUl33RWsaRxc6oA4DkWe39X1yCkyfRU8SyLszUTf7EBdWc8nQ2EKp7FrwnQKH702MDwH3C4JpPJxcXF8fFxkiTHx8cfPnyYTCZlFwpgPaEKOFCTyeTdu3fT6XQ4HE6n09FodH9/3+125SrYhfRbZRenkpq1gMSBLJhxIMWgrjY0sGq1vYuLi+l0enl5+eLFi7BlNpu9f/++1Wp9/Pix3LJBUrUv1JbiVqqWH9EGDavtYezdAykGdbW2gS393VmJFnh8fDwcDt+8ebO48dOnT1mW3d3dlVUqKNTyx1yoeg7Df9AI+7kiU0T39/d3d3enp6dL209PT8NdpZQKYAOhCjhErVar3W5/+fJlafvt7W24q4xCAWwiVAEH6vz8/Pr6ejabFVtms9nNzc35+XmJpQJ4SLMGOw9kcPdAikFd1aaBTSaTbrf7+vXrbrd7enp6e3t7c3Pz+fPn6+vrV69elV06qM93bZE5Vc+hpwo4UCcnJ9fX161WK8uydrudZdnLly8lKuBgNStCHkhkPpBiUD9VPMVvS741HKBaNks9Vc/hMjVQH4368QI4NELVXrn8MADUVeNC1eal93edcqQogBrwFzJrNS5UafoAPJNDCWs5+w8AIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAGm3zBdzYnlAFAI02n8+XcpWY9TRCFQA03WKuStPUxQ2fpnEXVIZg6e8wvyBAwxW5yu/hkwlVNFT41fAHGdBwqyN9i1v8Qn4XoQoAmquITeGPTJ1Vz2FOFQA03WK3/eq8dbYkVAFAo61OhNBT9TSNG/7bnL41IwCC5pzOUuOq7VnjQpWmA8A2nM7C9zL8BwAQgVAFAMsmk8nFxcXx8XGSJMfHxx8+fJhMJmUXikMnVAEHLU3TMLWluAG7NplM3r17N51Oh8PhdDodjUb39/fdbleuYrNmDRUbGmeJJgGsuri4mE6nl5eXL168CFtms9n79+9brdbHjx/LLduuxf1VbNpvbMNq27C9y6M0CWDV8fHxcDh88+bN4sZPnz5lWXZ3d1dWqfZDqHoOw38A8C/39/d3d3enp6dL209PT8NdpZSKShCqAOBfWq1Wu93+8uXL0vbb29twVxmFohqEKgD4xvn5+fX19Ww2K7bMZrObm5vz8/MSS8Xha9ZgZ9MGd3mUJgGsmkwm3W739evX3W739PT09vb25ubm8+fP19fXr169Krt0u2VO1XPoqQKAb5ycnFxfX7darSzL2u12lmUvX75sQqLimZoVIZsWmXmUJgFs1rRfCT1Vz1HJnqosy/I8X9rY7/ezLMuyrN/vl1AmAKDZqndB5SzLxuNxnudZli1tDLfDvaupCwC2V6zgH240qseFp6lYT1We50V4KvT7/fF43Ov15vP5fD7v9XohV5VRQKrBVb2AR82/VXZxqICKhaqzs7NOp7O0cTAYJElSjPqFGwYBeYirekF06bfKLg6Uo0qhKk3TTqeztgtqKWl1Op3VDi0Ifv311z/84Q+Xl5dv3rxpt9s//fTT5eXl69ev//73v5ddNKiqoi9Hpw5NVplQFXqeHhrUW5xfBZtdXV11u93iOqlJkvz444/dbvfq6qrEUgFQddUIVXmeDwaD0Wi09q7VjRsyVvpU8WpDmVzVC4AdqUaoClOp1kal7+2jmj9VlIpQOlf1AmBHKrCkQhj4W1qAKs/zYmGq1ac49Y8NwlW9Li8vf/zxx7DFVb0AeL4KhKognOJXGI/H4/G4SFRSFNv785//3O12379/v3pVr7KLBlAa63I9XwWG//r9/upIXFiVKoSq1XP9wrJVpZSWw+eqXgCrTHp5vgqEqkcV44Phv+GGdarY4OTk5OPHj2Fa+t3d3S+//CJRwZNZTReCOoSqLMtGo9F4PA6n6Y3H47XnCQIQndV0oVCry0eHmVWb11OoU315Pk0Cnuni4mI6nV5eXhZrv81ms/fv37darY8fP5ZbNkrXtN/YhtW2YXuXR2kS8EzHx8fD4fDNmzeLGz99+pRlmYXfaNpvbB2G/yAu677ClqymC4uEKljmKmawJavpwiKhCoCnC6vpzmazYovVdGmsZg12Nm1wl0dtaBJaC2xjMpl0u93Xr1+vrqZrpRKa9kOqpwqAp7OaLhSaFSGbFpl5lJ4qiMi3hiVNaxJ6qgAAIhCqAAAiEKoAACIQqgAAIvih7ALs2+YFshs1nw4AiKhxoUpsAoiu+Hs13PBLSzM1LlQBEJ0UBYk5VQAAUQhVAAARGP6joUwBASAuoYqGkqIAiMvwHwBABEIVAEAEQhUAQARCFQBABCaqA1W1dNUpJx8A5RKqgKoKKSpNU3EKOASG/wAAIhCqAAAiEKoAACIQqgAAIjBRHaieuk5Odz4jVJqeKqB65vP5Uv5Y+m9FzefzEKSKG0CFNK6navMvr18xqIqQq6yqAByOxoUqv7xQG0V/le81cAgaF6qAqlvtb17cImABZRGqgIopYlMY9dNZBRwIE9WBSlqcR7U6bx1g/4QqoHpWZ6brqQJKJ1QB1SNCAQdIqAIAiECoAgCIQKiCfzHZGYAnE6rgX+p68RMA9kCogm8s5ioXPwFge0IVLLOeJABPYEV1+D8uflI5i32KiX0ElE2ogv/j4ieVY+8AB8XwH3zDxU8AeBqhCv7FxU8AeDKhCv5FhALgyYQqAIAIGjdRffMUGR0VAMDTNC5UiU0AwC4Y/gMAiKBxPVUA5VqahKD7HGpDqALYq5CiXFkS6sfwHwBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEVQpV/X4/y7Isy/r9/vfeCwCwU9VY0jfP87OzsyRJOp1OkiTj8Tj59toO4bIPD927+LBK1JdDoLWwC5PJ5Ndff726urq7u2u32+fn53/+859PTk4WH1OztueyPE1Ws8b8qGr0VIVENZ/P8zzP83w0GiVJUvRIZVmWJMloNFp7L8CBmEwm7969m06nw+FwOp2ORqP7+/tutzuZTMou2g7N5/NwWC1uQF1VI0Kmadrr9RZzUpqmnU4nz/Pk9z+DljquinuXXqcS9eUQaC1Ed3FxMZ1OLy8vX7x4EbbMZrP379+3Wq2PHz8WD6tl26tlpXhU0/Z7NXqqRqPRas9T6KBKfh/1AzhwV1dX3W63SFRJkvz444/dbvfq6qrEUgGxVCxChs6nYjSw2Bi29Hq9JEkGg0FiThXPprUQ1/39fbvdnk6n7XZ78/Zatr1aVopHNW2//1B2Ab5PCE/J7/kpyLKs0+mMx+MQp5buXbI0ZXJ7jWoWQHStVqvdbn/58uWnn35a3H57exvuKqtgQCzVGP4rzOfz0WjU6/UGg8HiRPXxeNzr9cIsyHBvMTi4+gpPs7c6AnV1fn5+fX09m82KLbPZ7Obm5vz8vMRSAbFUtV8uBKlQ+NVp7Iv3LmpaPyTPobUQ3WQy6Xa7r1+/7na7p6ent7e3Nzc3nz9/vr6+fvXqVfGwWra9WlaKRzVtv1egpyrP89Vup2JLmGW19IDw39Wz/wBKdHJycn193Wq1sixrt9tZlr18+XIpUQHVVYFQlSTJeDxeOvuvCExr89PapAVQupOTk48fP97d3SVJcnd398svv9Q4UT15DitUVAVCVchGg8GgSE79fj9Mogr/7XQ6i1Oswr3WWQAo13w+X8pVYhb1VpnBzqWv4tLanmES1UP3Lr5IVepL6bQWdmpDA6tZ2wvVWV2lmSaoWWN+VJVquzTk99ADNoz6NW3v8hxaCzUmMhcAAAlOSURBVDvVnFCVrLvuBQ1Rv8a8WcNq27C9y3NoLexU7UPV5pG+GlSQbdSjMW+vYot/AlAJxaHU8B/NUYGJ6gBU1GJHxeq8dagZoQqAnVgd+tFTRb0JVQDshAhF0whVAAARCFUAABE4+w9gr4rJ2k6Ig5oRqgD2SoqCujL8BwAQgVAFy9I0DeMyxQ3YNS0NaqBxw3+unMCjNAP2LyyMudj2mnZ9D6iBxoUqP1LAYVrMVRIVVFHjQhXAwarlNfKc7UhzCFUAJVudlrC4peoppOrlh+0JVQAlK2JHGPXTowMV5ew/gIOwOI+qiFZAhQhVAOVbnZmupwoqR6gCKJ8IBTUgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBE8EPZBdi3zRd+d/ktAOBpGheqxCYAYBcM/wEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBE8EPZBdi3NE033Dufz/dWEgCgThoXqsQmAGAXDP8BAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARFClUNXv97Msy7Ks3++v3pvn+YZ7AQB2Kq3EBYbzPD87O0uSpNPpJEkyHo+Tby+N3O/3B4NBeMDqvYU0rUZ9gSbzS0VtNK0xV6OnKiSq+Xye53me56PRKEmSokcqz/PBYNDr9cIDwr1ZlpVWXACgeX4ouwDb6vV6xe0QmPI8D/8N6arIWFmWLT4YAGAPqhGqRqPRas9TsWU8HodhwYJpVQDAnlVssDP0ThWjgWFjmqahayrP8xCwwpT21ac3bXAXqCK/VNRG0xpzxWqbpmm40ev1QndUMYc9SZJOp5NlWZixvrZzq3j6E1TrgwKqq2nHIWqsaY25GsN/hWKuekhOi8N8xW7r9/tpmp6dna3dkY3auwDA3lQ1QmZZNh6PQ+HTNO10OsW89eT3FRZWq9a0yAxUkV8qaqNpjbkCSyqEVT2XNm6zYsJizAIA2KkKhKokScbj8dIJfYuBqVjwc+leS1UBAHtTgVAVstFgMFhcmGo8HheLUYW8VUSopXsBAPagMoOdSyfurZ1E9dC9iy9SlfoCjeWXitpoWmOuUm2LnPTQuN6jo35N27tAFfmlojaa1pgbVtuG7V2gWpa65P1eUXVNO+xWbJ0qgBpr1OEH6qcCE9UBAA6fUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQQeNWVF+6CsQSyxkDAE/TuFAlNgEAu2D4DwAgAqEKACACoQoAIAKhCgAgAqGq8jafz1hR9atU/WqUqFR11K9S9atRUtNKNY1QBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEEE6n8/LLsP+uFwlAOxTs2JGo2oLALAjhv8AACIQqgAAIhCqAAAiEKoAACIQqgAAIvh//X6/7DI0RfioX79+Hf6b53lxu6IWa9Tv9/M8TxYqWFG1r1SxJcuykooTR1GpsJsK1a3X4m7K8/y//uu/knq1vVp+oapeqUcPTKGC1f1a7Zmeqv0JX7zC2dnZ0pbKKcqfZVm4fXZ2VvXvXlGpNE2LSlX9b4+lltbv9weDQUlliaao1GAwqPpXKShq0e/3z87Oknq1vVp+oRYrVdFGuPnAFCqY57lVHrc1Z2uj0Wg+n/d6vXCj2DIajYotvV6v1+stPSs8pdPphId1Op3w4RfPKkuUGo1Go8WGVHqjil6pXq/X6XT2UfSHxWp7YeOBfPdjVeoQ6hJEr1HYuPuCbxL9C7X0i1GKKJXq9XqHUKlYrW71wLT409fpdJZegbUO5cfo8BWHovBFCs2r0+mEhlj8FBZbwrOKB6+mqNJDVcQaLdVr3zVZEH03zQ8gVMWtVLh9CEe1WJUKWxaDYyli1agIUuVWJ4i7m4qDdz2+UEu/DKV8p3Z6YFoMUqX/DFaFULWt0br+mMVfh7Xtb+kpBxiqiv8+v0bz3w9vOyz0Y+JWanQYnToRK9U5mN6dWJUKrxP+EA8Hj/2Uf1WsGhWHunAUrM0XqqhXzdpe8bASQ9ViwebxDkxLP4ZC1TbMqXq6MPBcTCEaj8fF7SzLxuPx0lh15/fO1YP1nBr1+/00TXu93qHNlnhOpbIsCz8lhzZR7GmVCrvm0OpSeFqlsiybz+f9fr/f78/n8/F4vL8SP+Y5bW8+n+d5PhqNDm3229Mqlef5YDAo0tWhNcInt71er5em6UHNN4p4YDr8Y9YBEqqiebT9HdTP/Ta2r1GYqB6ObTsv1vNsWakwNzP5vWoHvu+2r9R4PC6OAcUc28PU2C/UYuA4tPCxavu2V/xx0u/3D3zfbd/2QpQ/5J++Z36Pip+IQ/6tOChC1fcJDWvpT4Egy7Lie9Xv94u/xg68LT6/RuFQfVDVjFKpcPrV6l1liVKpoo86SZL5fF76Mfv5lQpdpIuvVq7n12jxWYdQoyRSpYqDd20qlSRJ0fb6/X6J/To7OjAVq0UkSTIYDA42OB6WfY83VlYxvWZxokPn2xMiVj/YYp5H+PcA51Q9v0aLUyVKb1QRd9Piw+qxpxZfs/TvfsRKFf8td09FrNHixtq0vQPZTfOdVWr/FZnv/sAUpmclpc5WrBahaltbnpITZi+ubtxRqZ6jfjWax67U2oftnz21unHtI+OU7KnqV6O5Sh38r8QefhwOYZdViFC1rVHtzn2oX43mKlUd9atU/Wo0V6mDV6e61IM5Vd+h9Ako0dWvRolKVUf9KlW/GiUqdfDqVJcaSOffDrgCAPAEeqoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIvj/HMQm89ErhtsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "g['tau_MS3'].Draw(\"AP*\")\n",
    "g['tau_MS3'].GetHistogram().GetYaxis().SetLimits(35,47)\n",
    "g['tau_MS3'].GetHistogram().GetYaxis().SetRangeUser(35,47)\n",
    "g['tau_MS3'].GetHistogram().GetYaxis().SetTitle(\"#tau (ns)\")\n",
    "g['tau_MS2'].SetMarkerStyle(24)\n",
    "g['tau_MS2'].SetMarkerSize(1.2)\n",
    "g['tau_MS2'].Draw(\"PSAME\")\n",
    "l.Draw()\n",
    "c1.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "id": "8fb2f65a-e5b3-45f4-90e1-925b0be0e7f1",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dv5LbRqI3bEDfSewq18xsLLukO7CuAKCuZJOVVevQPoF1qpagE2+wCteuHZ1gb8QD4Aa2dAWWgomPNNFk4ht0GR8McjjksEniz/MEKgogOd0DkPhNd6M7XS6XCQAA+3l06gIAAIyBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEEEfQ1We51VVrW6vqirP8zzPi6JY3VsUxYa9AAAHlS6Xy1OX4Q+qqprNZmVZ5nne3l4UxWKxSJIky7K6rpMkaZc8z/OwMciybG0sAwA4kB61VFVVVRTFbDZbu2uxWMzn8+VyWVVVWZZJkjSpqyiKuq7D3uVyOZ/P67oWqgCAY+pRS1Waps3jTktVaIhqFzX08YV/wwvbe9M01VgFABzTf526AP+/kIpC919nV13XWZa1t3QGTnX2Nl2EAADH0aPuv83CCPQ8z9M0XR3J3hmA1fkvAMCh9ail6i4hPzWj1Ofz+WKxaAaz79TH1+5hBAAOrT+jjI5gAKGq0RyYoijSNJ3NZsvlctdGqUkdXQA4oam1ZQyg+y/Eps6oqfl8vuElhqgDAEc2gFC1QROepCgA4LSGEapW7+YLKappxOrsXb1bEADgoIYRqsIECquzfa7dGx5YrAYAOKaHDFSvfpf8nmDConsxy/VHeZ6Hm/6aIW9ZljWxKc/zsixns1mzd3WVGwCAg9phRvWwjMyGSTXn8/mh24faSW7XvUmSpGmPZpAHgHGb2mV32+6/PM/DROdlWS5XlGU5n8+rqkrT9KC5anOT2KEbzAAA7rJVhEzTdPtWqBBr+nk73tQiMwCc0NQuuxOr7cSOLgCc0NQuu/vOqH7vMCYAIJhUwpignUNVURSLxSLcXhceh+1HGKUOABN3b1tGM/xm9TlVVTV3m2VZpkEkvtVR5xuUZdl+VXhclmWYaXOntzqJuL8NANjJp/00EzTO5/O7ntNc0dobr66u1l71rq6uOi+/urpqZs/Osmz1CTuZ2oV1t8k/Q1vUcrlMfs/CTZNV0tfB6R2bfx2nLh0A3K/pJupYeyGuqur58+dJkszn8xCSrq6uQj57/vx5+yXhmU1rVl3Xz58/D/f+s42dZ1RvAmxn8s9kIKEKAEZg7TV3bdgKG0OQaq7a8/k85Kr2S0L2ardOJUlS17Xr+5Z2C1V5njcBdrFYrA1YAMDhhDC0di7utUvfhmeuXqM77xMu5Z2xVqHfcMO837TtHKqS1hybzSp7oW1QqAKA41htlApbmnFX22gPtwrJae3LtVRtaedQNZ/P67oOWTgMpQpHsRnDDgA8wKNH216UQ/TpZJ27eo1C29WjR4/Wdg42zw+D3zsvDy/RaLKlac3KNbVZyADoldVr0KNHj5pb9tqP11osFovFIgyHevToUZZl7RaNZksIZ+236sS1+Xy+eUqFMPlC+LfzU3by6NGjSV129538EwB4sE+fPnWy1L3RqtEe6rS57+/Tp09VVS0Wi/CS9hyTa18SRqw3z9mmMCQPuPuvqqo8z9N1DlE+ABi3kKvC4+0TVacH8N47xvI8L8uyPZ9CkiSLxWLtjAmdaRfumsGBjt26w6qqCr/91ZsLkiEMZNP9B8AJbbgGrfbZrWp3/yWt/r7Vx/e+VdKav+rq6mpDGtvy3e567aQuu7t1/7Un/wQAHuauMelrt28INO3uvLX9dCGHrY1N4eaz0CeY57nZkfb38Mk/AYCHWV3OpUlOaxd7WavpAdyQh8JVe5v+u8Vi0R5KxQM8fPJPAGB/7XFU7fFV92qmRF8752cQktZds6K3m7jCMzvxq/8De3rlIWv/DbptcO0Qe2PtATiJ1ZHpu45e2jBpZxCm9wyr+C0Wi9CytVgsQnprL9Kc/N5d2DynGXS1W62marcxVVVVZVlW1/Xa/DGIsVaDKCQAE/GwAeCNMCgque++v6urq7BScqe7qTOlQnhap7Fq8zB22na++y80Vt21d/8CHZS7/wA4odNeg8JknuHxhsk/m6dtniB0G1O7+29aIUOoAuCEpnYNmlqo2mpMVZqmofHw0KUBABiorULVcrn829/+9uOPP4Z0ZcZ6AICObe/+e/78eVVVIV2FgerSFQBAY+fJP9vpKvl9hgLpCgCYuJ1DVSPcdblcLn/99dckSdI0NegKAJisre6GK4piw0wKA+LuPwBOaP9r0L0r9DXTG60+Z8spFSKa2t1/W4WMMNVnlmVFUQx6BjChCoAT2vMaFKY7T1Ym7WxrVrlpTytaVdXadf3WTuzZ/JQQvPYZ4TO1ULVV9998Pg8Tqc9mszBEvf/zfALAWN21QPLaq3OTqObz+dXV1adPn66urkJUCuOk28989OhR8+Z1XTer2bCV5S5Cumpem2VZWZY7vcNp7VpfAIjo0346S8qsPqF9je5sXH1+eLcsy5ot4YXz+bzz2vZzdjK1y+4Du8OKomh3zc7n8zzP+98zeO+SyQ/7bQDANva8yoSOubDe39oewEePHoWepaTV/Reamj6tW2Swvatp0Oo8c8PL76X7byshVC2Xy7IssyxbLBaz2SxuyQ5kc8Y8dekAmK7OasebrfYANsOttn+TT58+NXfuh5+++vLQWGXYzzb26igN6yvvdBIAAG3X19f/8z//86c//Wk2m/3pT3969erV9fX15peE6NMJOnfdGBhSUXuwVFvz/DDiajVUhat8/zuj+uAhoaqqqjzP0zSdzWZ1XTcjq6IXDgDG7fr6+uXLlx8+fPj111//7//+7+rq6ubm5ptvvtmcq0L06YSkcEVefXJZluFBGHUe0tXalqfV5BS6oda+Lat2CFWrWWo+ny+Xy7D9YCUEgNH65Zdfvvrqq3/84x/Pnj07Pz//+uuv//GPfzx58uTnn3++97XtnqLNfX+hm6/JRovF4vnz53e1XTXC5T5pxTI22ypUFUXRZKkkSZosNY4ZQQHgVH755Ze//OUvn3/+ebPls88++8tf/vLLL79sfmGnB/DeSUHzPC/L8lNrPoUkSe4aEh3atEIDysOGqE/TDpN/DuUWvw1M/gnACTXXoJ0mf2piTXP3X0hF4V6/0IzUeZxscb9ec7tfewrQ9jSha6cG3Ym7/9YIQ6aGPp06APREM5PT+fn5f/7zn870Tv/5z3/Ozs4600etFXqQNvT9hTanu0ZQhZc03YjtaUI/ffrkor+rrULVXb9WN1gCwD5evnz55s2b29vbZsvt7e3//u//vnz58t7XNj2AG/r+wjiqzWOnGk3D1T5L00zZQ+7+C8PVwyirJEnCwjWRywUAE/DXv/71/fv333///du3b29ubt6+ffvf//3fv/3221//+td7X9vcA3jXfX/J70mrruu17SDtJq7whDDU54GVmbydQ1WapnVdt9erCZO3OgYAsKvHjx//61//Ojs7m81mFxcXs9nsiy+++Ne//vXll19u+Q53TdrZCNN7Pn/+fDabhckUqqpqFvVrD1pPWtMudOxZzanYaVGb8KsPQ6zC4/b2nd7qJAZRSADGavNKeWuX82sLV9v22nxNJOqsuNfZ0kyb3rG6zN9dNhfsLlO77O52N1ye53Vdh5cURbFYLJqXp2lalmXP26vc/QfACZ32GtRetDfLsiNcsqd2999/nboAx7Z5TeVJHXsAJmXo8yL13269pGG2z9XBbuEgDeJQbW64O3XpAICh2q2lKs/zLMtms1nT8xo6AZMdl8UGABiZh4wxaoJUYz6fD2LJGmOqADihqV2Dpjamaq+QMbillIUqAE5oatcgoWrMhCoAOJqpXXa3GlO1ZXOUVWsAgMnaKkJunoagESuN5nm+unjzahk6A7mKomgWP7prgNfUIjMAnNDULru73f2XZdlq3ImrPTVZe+PmV4VZScPjsMKRZjMA4Ji2mqcqzOE0n8/rup7NZmEF5eippaqqoijCIs1rheVxGk1zVFEUYTnCdjmFKgDgmPadUiFi21W7g6+z4k1nSZy1L2zvTdM0y7LVXDW1dkgAOKGpXXYfsu50URTL5bIsy07b1Z5FCe1MZVmu7mri0V39ep1lILMsW+1DBAA4nIeEqiCMBw/dbUmSHCHEpGk6m83WZrjN/wUAOLSHh6ow82eapkdYpiYktmbUVGiIumshws3ShzpAtQCA8djt7r/k9+HkTbvUfD4/wqrXnSFWIc8tFosHDOeaVOcuAHA024aqk2SpxuoPCnMo3LVOjlv/AIAj2ypUNZ1fR85S92pKIkUBAKe125iqxWIRxokfc9RRVVVpmnYmSW+nqNV7/eq67twPCABwUFu1VJ02oOR5nmXZYrFoGsma2T7DE8KUoc18pM1zTlReAGCKtgpVJ+9cC41V7cnW2wv/5XlelmVoQgtbOgPbAQAObdsFlbMsm8/nG9aQOYIw8+eGQV3tlqq1pja1Kw/T6ct2zgA8zNQuu9vWtizLxWIRxiqdPF092NSOLvtwtgDsaWpfpDvXtp2uwqTqhynYQUzt6LIPZwvAnqb2RbrzjOqz2ayqqrA6TRjqtHprHgDA1Dx8mZomXV1dXSVJkqbp2rWQAQCmYK92ubsmNO+tqbVDsg9nC8CepvZF+pCWqrCOcjPHwbC6/6yaDAAcws4RMiSPMKCqruvlchmW4cuy7OTTWd1rapGZfThbAPY0tS/S3VqqQotUWZZFUbTX3ZvP552FYgAAJmW3UHXX7JohbPW/pQoA4EC2WqYGGDrTxAMc2kO6/1ZbpML2Yd0JCJOyXC5DkGoeABDXbi1VeZ5nWTabzbIsC1uKoggj1ufz+QGKBwAwDA8Zll8UxWKxaG+Zz+eDmFVharchsKv2GTLKs2WUlQJ6a2rfOSb/hD9oTpLwYGTnzMiqA/Tc1L5zHjhQvaqqZmTVsHIVbNYOUlP7OgBgHztfM6qqChOptw1i5s/ENZKthXvlxne2+AgAxzS175ydl6kJiWo+n5dluVwuy7LMsqyua+1VDN3qgkWWMAJge7t1/zUzqjcRKs/zqqrC0PXBDbGCts4Q9bE2VgFwIDFnVB8ECyqzWbuxuolWAHCv3ULVXQ1RgxhQFSw3OnXpOLHV7n9nBQBb2nkEWZqmq7NSpWk6iLHqUxsxxz5GebaMslJAb03tO2e32obhU3VdJ0nSTKre/LfdjtXPDsGpHV32McqzZZSVAnprat85u9V2dS71u/Tzlzi1o8s+Rnm2jLJSQG9N7TtnYrWd2NFlH6M8W0ZZKaC3pvads/M8VQAArNoqVKVpuv0YqaIozFYFAEzNVqGqLMvFYpGmaZjqc+1zQpZK07S9LCAAwERsNaN6nufL5TLc+re68F8jy7JJdZ0CADQeMoKs+l34b/67uCU7hKmNmGMfozxbRlkpoLem9p0zsdpO7Oiyj1GeLaOsFNBbU/vOcfcfAEAEW42pGpPN6+NOKlADABFNLlSJTQDAIej+AwCIYHItVXCvpo84PNC6CcA2dg5VYWr17SdYh8GRogB4gJ27/xaLRWfC9KIodlrHBgBgfHaeQCJN0yzLOrkqz/O6rvv/9/3UJsyADh8B4Jim9p3zkDFVdV2Hdqkwi3pVVXVdRy0VAMDAPHygenulmizLdP8BAFMWp/tvKKbWDgkdPgLAMU3tO2fngeplWersAwDo2DlU5Xm+XC6LohhoYxUAwCFMq11uau2Q0OEjABzT1L5zJjejugWVAYBD2Kr7L03TPM+vrq4OXZojWG506tIBAEO1VahaLpdZlv34449jSlcAABHt3Nl5dXX1448/1nWdZVmWZYvF4kAlO4Spde5Ch48AcExT+87Z+e6/58+fV1W1XC7/9re/1XWdpmmapvP5/BCFAwAYip1DVaNJV7/++muSJGma6hYEACbrIe1yzQxVYe2/AZlaOyR0+AgAxzS175wdWqqqqsrzPE3T2e9C31/0WUDzPN/8nnmery41WBRFnudrdwEAHNq281QVRRHGpGdZ1jRQVVVV1/VsNpvP57GiTHjPDU/I87yu604jWdgYHtd13V7sGQDgCLYKVVVVLRaLu9ZRzvN8sViEVqJ9ihKS0ObbCddGrqIo6rpugl3If6FdbZ/yAMBanXmkJ9XDxQZbdXaGs2fDM+99wlZFaZ2jZVmujURpmmZZ1o5Qa396eNpqBJxa5y50+AhARD5Q95rar2jbMVWbJ00Ie/fscQtzmpdledcTQsxa+1OyLOv8d3MfIgBAXPeHqhBiNnelHWFseOjjuyvwrg6xOnR5AADa4iyofOhR4WGs1dpGrF1/9OYFlTeYVAMmALCr+1uqNnS6NbZpzdrHbDZr33XYtusP3bygsrWWAYCH2balarFYbOjjO+gKgOHndiagqqqqmZhq9SXmU4DG1AaKApzKVqGqLMsw1efam/LClg0DzKPo5La6rtuzVUlRcJflctnJVWIWwCFs+926YfLPJEniTv45m83umlIhCEs4Nz8xzPxpSgXYIJz8UWY/AQLXlHtN7Ve0w4zqobstNBE127MsC9sPUrrtFEUxm82axW1CYSxWA20SFcCh7Xb3X9P206v5yvM8bzoow5bNDV0wHau3u7a3CFgAEcVplwsrzJy8cejemxCn1g4JDd1/g2Dxk2FxTbnX1H5FcWobRlz1/xc3taMLQXPmN9HKB6HPHKChcKTuNbVf0bbL1AADtfqlNqnvOICjEapg5EQogOMQqgAAIoiz9h/A8RnWDfSKUAUMVXv0/anLArBdqFqd6ma4NtfFVzMA8DBbhaosyw5djqMRmwCAQ5hWs7luAiZulB+B8VVqZDUa8dC3kR2pQ5jar8iYKgAOyNA3psOUCgAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVcDwjGmZB2A0hCpgeJbLZSdXiVnAyQlVwCC1c5WJJYE+mNyM6hZUhtFocpVPLtAHkwtVvnxh6Fb/NGpv8RkHTmVyoQoYuiY2hV4/jVUck75mNjCmChik9rVtddw6HIibJNhAqAKGZ7W1QOMBR+MmCe4iVAHD4zLGael3Zi1jqgBgK26SYDOhCoCDGF/XmJsk2Ez3HwAHMdYx3W6S4C5CFQCHMr4x3W6SYAOhCoADGlk32ThqwYEYUwVAfMZ0M0FCFQDxGdPNBE0uVFlQGeBoVsd0+5rtrc710ZF6gMmFKmcJwHEY0z0s4egIvvswUB2Ag3BtZmqEKgCACCbX/QeMRnsCpES7CHBqQhUwVFIU0Cu6/wAAIhCqAAAiEKoAACIQqgBg0jZPi832hCoAmLRmHaGGmPUwQhUATF07V5lU/cGEKgAgse71/iY3T5UFlQGgsXpZbG9xWdzJ5EKV8wMAGs1lMfT6aazah+4/AJi69jiq1XHrbEmoAoBJWx2ZrqXqYYQqAJg0ESoWoQoAIAKhCgAggj6GqjzPq6pa3V4URZ7neZ4XRbHrXpi4NE3DyNPmAQBx9W7W1KqqZrNZWZZ5nnc2JkmSZVmSJHVdJ3/sA87zPGwMsixbG8vMEgv038i+qTohfmRVG1N1griVGuWvaIMetVRVVVUURQhPHWHjcrmsqqqqqrIskyRpWqSKoqjrej6fL5fL5XI5n8/rul4bqgA4suUfnbo4cEA9ipDtv2Y6LVVpms7n83a/XpqmTXPU6kxl7b2dH9Gf+gKs5ZtqKEZ5pLRU7aNHLVXhj5jQCtXRyVhBe0voFmz/t90bCNBzBrrBCPQoVG0QRqAnSRK6/8LjdqjqRK7VBAbQZ6tzWItZMDgDW/uvGXE1n8+bmLXTOzz4e2pSDZjA8YVcFb5qptZpAuMwsFDVjFVfLBbJ79Mo7PoOBykZwN4sZ8sJNY0OTsIHG0b3X1uYiSrLspCr1nLrHzAgactdW+DQ3Ke5vwGEqjCIqpOTOg1UUhQwXO3LWHMxc22DwRlAqAoTe3ZiU/u/q/f61XXduR8QoOfa46hWx60D/TeAUBUsFosmSDWzfTb/TVptV+HBqRarSf/oJGUABmd1ZLo2qj6z7hNrDWOgevijrT3ZepZlTWzK87wsy9ls1pzZa+e1Og537gAP4BtjWBwv1hrStb9pqborMIUnbIhTR8s6QhXwYL5AGI2pncwTq61QBfSeLxBGY2on82DGVAEA9JlQBQAQgVAFABCBUAUAEMEwplSIaPOEIpMaTwcARDS5UCU2AQCHoPsPACACoQoAIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIILJrf1nQWUAJqtzEXTVi2tyocoJBMBkhYtgmqauhoeg+w8AIAKhCgAgAqEKACACoQoAIAKhCgAggsnd/QdAdG7Uh0RLFQD7Wy6XIUg1D+ih6+vrV69eXVxcJElycXHxww8/XF9fn7pQoyJUAcD4XV9ff/PNNx8+fLi6uvrw4UNZljc3Ny9evJCrIprW9F9Hm+7MvGrAgw33C2S4JZ+CV69effjw4fXr159//nnYcnt7+913352dnf30008H+qFTOyUmVluhCui94X6BDLfkU3BxcXF1dfXs2bP2xrdv3+Z5/vHjxwP90KmdErr/AGDkbm5uPn78+PTp0872p0+fhl0nKdX4TC5UpRudunQAEN/Z2dn5+fn79+8729+9exd2naJQIzS5ULXc6NSlA6AvRvZX98uXLy8vL29vb5stt7e3b968efny5QlLNTLT6uw0pgrov+F+gQy35BuMplLX19cvXrx48uTJixcvnj59+u7duzdv3vz222+Xl5dffvnlgX7oaH57W5pcSxUATNDjx48vLy/Pzs7yPD8/P8/z/IsvvjhoopqgaUVILVVA/w33C2S4Jd9ApQbxg3pCSxUAQARCFQBABEIVAEAEQhUAQAT/deoCAGzSmR9oUoNegWERqoBeCylqavcQAUOk+w8AIAItVUelIwOAE2ouQ+GBy1BckwtVm9dvOvTppSMDgBNy9TmoyYUq5xMAcAjGVAEARDC5lioIjG8DIC6hiokyvg2AuIQq7qdRBwDuJVRxP406AHCvPg5Uz/O8qqrV7UVR5Hme53lRFLvu5UDSPzp1cQDiuL6+fvXq1cXFRZIkFxcXP/zww/X19akLRd/1ru2hqqrZbFaWZZ7n7e3hgp1lWZIkdV0nSdJ+Tp7nYWOQZdnaWHa0tpbNP2igTT4bij3QGiVDLvnUTOpIDbeywy15x/X19TfffPPVV1+9ePHi6dOn79+/v7y8DP8+fvz41KUbktGcElvqUUtVVVVFUcxms9VdITyVZVlVVVVV4Qg1zyyKoq7r+Xy+XC6Xy+V8Pq/rem2oAoB7/fzzz1999dXr16+fPXt2fn7+9ddfv379+smTJ//85z9PXTR6rUcRst151GmpStO00/hUFMVisWjG+iR/HD29+vxmu5aqB9NSxQlN6kgNt7LDLXnHxcXF1dXVs2fP2hvfvn2b5/nHjx9PVaohGs0psaUetVSFdqayLFd3ZVnW6Q1cfULnv+3eQIBBaMYmGqR4Qjc3Nx8/fnz69Gln+9OnT8Ouk5SKQRjG3X+rbU6LxaL9307k6gyxAhiESf1N31tnZ2fn5+fv37//+uuv29vfvXsXdp2qYPRfj1qqtlRVVfgDLrRp7Tp2Kn2oPYvtj06AoXj58uXl5eXt7W2z5fb29s2bNy9fvjxhqei/gYWqPM/D+PRm0NXmbsFVy4fas+TL5bKTq8QsgH769ttv379//9133719+/bm5ubt27fff//9b7/99u233566aPTaYEJVaKBq7vLbnKX6eetfO1dNbewewIA8fvz48vLy7Owsz/Pz8/M8z7/44ovLy8svv/zy1EWj1wYzpmo2m901+1TS1xS1qslVEhWHYEEhiOXx48c//fTTTz/9lKapwelsaRihKnT53ZWcVu/1q+u6cz/gaa329LW3uPIRiwWFAE5oAKGqyVKr68+ELWHK0GZxm9Az2KvFaporXLjaaawCgPEZQKhqdKZRSH5PTnmel2U5m82a5p/VVW76oN1+EKKVXAUAozGq63q7pWqtE86o3t4y0DhlRvWhUClOZZSHaZSVOpqp/fYmVlvL1OxBqBoKleJURnmYRlmpo5nab28wUyoAAPSZUAUAEIFQxRRdX1+/evXq4uIiSZKLi4sffvjh+vr61IViDUcKGBChism5vr7+5ptvPnz4cHV19eHDh7Isb25uXrx44WrdN44UMCzTGkF273J7sX4bBqr32atXrz58+PD69evPP/88bLm9vf3uu+/Ozs5++umn05YtimEdjg1Gf6TGZzTnXtsoK3U0U/vtTay27v7bw2hC1cXFxdXV1bNnz9ob3759m+f5OBajGNbh2GD0R2o0rq+vf/75519++eXjx4/n5+cvX7789ttvHz9+fOpyxTGaD9RJTO23p/uPabm5ufn48ePTp087258+fRp2naRUrHKkhkIvLTSEKqbl7Ozs/Pz8/fv3ne3v3r0Lu05RKNZwpIbi559//uqrr16/fv3s2bPz8/Ovv/769evXT548+ec//3nqosGxCVXcY3y3X718+fLy8vL29rbZcnt7++bNm5cvX56wVKxypAbhl19+efHiRTPuLUmSzz777MWLF7/88ssJSwUnIVSxySgb9r/99tv3799/9913b9++vbm5efv27ffff//bb799++23py4af+BI9Z9eWmgTqthklA37jx8/vry8PDs7y/P8/Pw8z/Mvvvji8vLyyy+/PHXR+ANHqv/00rgkTUcAAA9ISURBVELbtIblu/tvV1vefjWgGnUMt+QbqBTH9MMPP9zc3Lx+/fqzzz4LW25vb7///vsvvvji73//+2nLFoVzbx9T++1NrLZC1S5ubm7Oz88/fPjQ+XNzdftQarRquCXfQKU4puvr6xcvXjx58uTFixdPnz599+7dmzdvfvvttxG0KXamNnQGPsDUPrm6/44qTdPwKW0e9JmGfeBeI+6lXf7RqYvDAEwrQvakpWpAtmzYH259h1vyDVSKU3GY6JjaKTGx2gpVO9qyYX+49R1uyTdQKU7FYaJjaqeE7j82GXHDPgDENa0I2ZMFlQdqNGv/tQ235BuoFKfiMNExtVPiv05dgGOb1NEFAI5G9x8AQARCFQBABEIVAEAEQhUAQARCVUzX19evXr26uLhIkuTi4uKHH364vr4+daEAgGMQqqK5vr7+5ptvPnz4cHV19eHDh7Isb25uXrx4IVcBwBRMawKJg06Y8erVqw8fPrx+/frzzz8PW25vb7/77ruzs7OffvrpQD/0mMxTNRQqxak4THRM7ZSYWG0PeXQvLi6urq6ePXvW3vj27ds8zz9+/HigH3pMQtVQqBSn4jDRMbVTQvdfHDc3Nx8/fnz69Gln+9OnT8Ouk5TqoO6dnh4AJkWoiuPs7Oz8/Pz9+/ed7e/evQu7TlGow1oul51cJWbBNtI/OnVxgGiEqmhevnx5eXl5e3vbbLm9vX3z5s3Lly9PWKqDaueqqbXxwoMtl8vwYWkeAOMwrQvhQRdUvr6+fvHixZMnT168ePH06dN37969efPmt99+u7y8/PLLLx/8tv1xV2wKv9WBnkijzIIqNQjjq1Ey0kqxj6mdEpNrqVputM87P378+PLy8uzsLM/z8/PzPM+/+OKL0SSqjtXOi8F1ZzRFHVCZAeizaUXIo0XmUWbztZUKGwfdWDU+0zn9Bm18NUpGWin2MbVTYnItVUTU/rSsjlsHgEkRqnig1b8/JvXnCAB0CFU8kAjVN5aeBDgtoQrGwNKTACc3rRFkBqrvY5TL1IzGlJeeHKjx1SgZaaXYx9ROiYnVVqjag1DVZyNeerJzA8RozrRRfmpGWSn2MbVTQvcfDN64l56MOJkcwEEJVdzPPJk9N8GlJwF6SKjifpoK+m+CS08C9M20OjuNqTqQqdW3h0a/9OT4jPJTM8pKsY+pnRKTa6lKNzp16eCBJrX0JEA/TStCaqk6kKnVt+ccjkEY5WEaZaXYx9ROicm1VAEAHIJQBQAQgVAFcDyWaIQRE6oAjsQSjTBu0xpBZqD6gUytvj3ncPSWJRqZmqmdEhOrrVB1GFOrb885HL014iUaA+ceHVM7JfrY/ZfneVVVu+4tiiLP8zzPi6I4WNEAHmjcSzQCSZL816kL0FVVVV3Xu+7N87zZXtd1VVUbYhnA8TVLNH799dft7ZZohNHoUUtVVVVFUcxms133FkVR1/V8Pg8r083n85CrDltcgB1ZohHGrUedne1VYsqyzPN8y71hV7siaZpmWbaaq4ypOpCp1bfnHI7eGv0Sjc49OqZ2SvSopSq0M5Vl+YC9WZZ1/ruhDxHgJCzRCOPWo1C1j06zVue/AD3x+PHjn376KQxL//jx49///neJCkajdwPVd7Xr2Kl2N+JOJtWACQDsavChatdGKdkIADiEwYeqtdz6B3BMTSfA6p1DMB0jCVVSFMAJSVGQjGOg+uq9fnVdd+4HBAA4qDGEqrAuTTO4KjywWA0AcExj6P7L87wsy9ls1nTqr84dCgBwUKOa6jSMrNoQp8yofiBTq2/PORyD4DAxBVM7zydWW6HqMKZW355zOAbBYWIKpnaej2FMFQDAyQlVAAARCFUAABEIVQAAEYxhSoWdbF5QeVLj6QCAiCYXqsQmAOAQdP8BAEQwuZYqgNNqBiGEB5rPYTSEKoCjkqJgrHT/AQBEIFQBAESg+4+9GB0CAIFQxV6kKAAIdP8BAEQgVAEARCBUAQBEIFQBAEQwuYHqFlRmxNyMCXBCkwtVLjOMmNMb4IR0/wEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEMLm1/yyoDAAcwuRCldgEAByC7j8AgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACCa39p8FlQGAQ5hcqBKbAIBD0P0HABCBUAUAEIFQBQAQgVAFABCBUAUAEEEfQ1We51VVrW4viiLP8zzPi6LYdS8AwEH1bkqFqqrqul7dnud5s72u66qq2sFr814AgEPrUUtVVVVFUcxms9VdRVHUdT2fz5fL5XK5nM/nITltsxcA4AjS/kyG2Z7rvCzLPM87u9pFTdM0y7KQnDbv7fyI49T3aD8IAHpralfDHrVUhXamsizX7s2yrPPfdi/h5r0AAIfWo1C1Wbvh6t69m58MABDdAELV2tFRTWzadexU+lA7vX/7AQAwBb27+2/VTm1U9zp05+6kOo8BgMYAWqrW2txA5dY/AODIBhOqpCgAoM+GEapW7+YLE1Nt2Nu5HxAA4KCGEarCyjPN8KnwoFmOZvNeAIAjGMBA9SRJ8jwvy3I2mzX307Wns1q716wKAMAxDWyq0zB26q7AtHlvMr2pXQHghKZ22Z1YbSd2dAHghKZ22R3GmCoAgJ4TqgAAIhCqAAAiEKoAACIYxpQKEW1e5HhS4+kAgIgmF6rEJgDgEHT/AQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQwubX/LKgMABzC5EKV2AQAHILuPwCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIJrf2nwWVAYBDmFyoEpsAgEPQ/QcAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQweRCVbrRqUv3EAMt9mbjq9T4apSo1HCMr1Ljq1Ey0kpNjQWVAQAimFxLFQDAIQhVAAARCFUAABEIVQAAEQhVAAARDClUFUWR53me50VR7LoXAOCg0kFMMVBV1Ww2S5Iky7K6rsODqqqaJ+R5HrYHnb2NNB1GfXeiUoMwvholKjUc46vU+GqUqNQoDKOlKiSq5XJZVdVyuSzLsq7rpkWqKIq6rufz+XK5XC6X8/m8ruu1oQoA4ECGESHTNJ3P5+1+vTRNm+aoMAttuyLtvZ33GUR9d6JSgzC+GiUqNRzjq9T4apSo1CgMoKUqZKM8z9sbm37A5r8b9gIAHNoAQlWIU51mp05m6kQuAIAjG8zaf4vFItzclyRJux9w7dipzrj1tlGuWKlSgzC+GiUqNRzjq9T4apSMtFKTMoxQVZblbDYLw9WDpoNvpzaqSfXsAgDHNIxQlef5crls3+63OUu59Q8AOLJhhKqg3evX6d2TogCA0xrAQPUkSdI0XW2ams/n4cHqvX5h2qrjlA0AIBlKqAqxKbRUVVUVhvK1ewOT1uCq1cHsAACHNphZuTo39JVl2W67ataxWbsXAODQ/r+htOj8+c9/zvP8yZMnRVH8+9//fvLkSXtv2J7n+Z///OfVvT0RftVN2aqq6mc5t9euUVEUYWSbSvVQ59xLtrjbo/+aSoXD1BhuvdqHqaqqf//738m4zr1RfqCGXql7L0yhgsP9WB3ZMLr/gjzPN18Jmoms+qkzmn42mw19fH1T/jzPw+PZbNbnQ7CNplJpmjaVGsrfHnfpnGlFUSwWixOVJZqmUovFYugfpaCpRVEUoel9TOfeKD9Q7UoN9CTcfGEKFWxG3XC/JVsry3K5XM7n8/Cg2VKWZbNlPp83Szs3zwkvybIsPK1ZVKd51alEqVFZlu0T6eQnVfRKzefzLMuOUfS7xTr3wsaefPZjVaoPdQmi1yhsPHzBN4n+gep8Y5xElEqF26GaXaeqVKyzbvXC1P7qy7Ks8w6s1Zcvo/5rLkXhgxROryzLwonYfBU2W8KrmievpqiTh6qINerU69g1aYl+mJY9CFVxKxUe9+GqFqtSYUs7OJ5ErBo1Qeq01QniHqbm4j2OD1Tnm+Ekn6mDXpjaQerkX4NDIVRtq1zXHtP+dlh7/nVe0sNQ1fx3/xotf7+8HbDQ94lbqbIfjToRK5X1pnUnVqXC+4Q/xMPF4zjlXxWrRvPWZDFj+kDNW9PcHKXsd4p77jVPO2GoahdsGe/C1PkyFKq2MaQxVX0TOp6bIUR1XbendajrutNXnf3euNpb+9SoKIo0Tefzed9GS+xTqTzPw1dJ3waKPaxSnclH+uZhlWqWWyiKYrlc3rXo50nsc+4tl8uqqsqy7Nvot4dVqqqqxWLRpKu+nYQPPvfm83mapr0abxTxwtT/a1YPCVXR3Hv+9errfhvb1ygMVF+2lhLqrS0rFcZmJr9XrefHbvtK1XXdXAOaMbb9NNkPVDtw9C18rNr+3Gv+OCmKoufHbvtzL0T5Pn/17fk5ar4i+vxd0StC1W7CidX5UyAINyeGx0VRNH+N9fxc3L9G4VLdq2pGqVQz81lPqhalUk0bdZIky+Xy5Nfs/SsVmkjb73Za+9eoMwPfoQq6iyiVai7eo6lUkiTNuVcUxQnbdQ50YWpmi0iSZLFY9DY49sux+xsHqxle0x7okP3xhojVX2wzziP828MxVfvXqD1U4uQnVcTD1H7aOI5U+z1P/tmPWKnmv6c9UhFr1N44mnOvJ4dpebBKHb8iy8NfmMLwrOSkoxWHRaja1pa35ITRi6sbD1SqfYyvRsvYlVr7tONzpFY3rn1mnJI91PhqtFSp3n9LHOHLoQ+HbECEqm2Vo7v3YXw1WqrUcIyvUuOr0VKlem9MdRkHY6p2cPIBKNGNr0aJSg3H+Co1vholKtV7Y6rLCAxmQWUAgD7TUgUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQwf8DQjUT7KCna+YAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "g['lyDT_MS3'].Draw(\"AP*\")\n",
    "g['lyDT_MS3'].GetHistogram().GetYaxis().SetLimits(90,160)\n",
    "g['lyDT_MS3'].GetHistogram().GetYaxis().SetRangeUser(90,160)\n",
    "g['lyDT_MS3'].GetHistogram().GetYaxis().SetTitle(\"LO/#tau (pe/MeV/ns)\")\n",
    "g['lyDT_MS2'].SetMarkerStyle(24)\n",
    "g['lyDT_MS2'].SetMarkerSize(1.2)\n",
    "g['lyDT_MS2'].Draw(\"PSAME\")\n",
    "l.Draw()\n",
    "c1.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "id": "6c70fee4-b956-4d0e-b797-03f7b50c1945",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxwAAAI8CAIAAAD0vjrdAAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dv3LbSII/cNC/TcZVU5I2tl3WG9h3D0DQT7LJyqpV6LlgfAHBZCdZh2PXyhfMi5wIvsCWnsBWwHglRcqGv6DLOAwpUbTUJP705xO4JAiUuwUQ+LK70T1YLBYZAACP86TpAgAA9IFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQQQdCVZ7nZVneu1tRFHme53leFMXWywQA8EdtD1VlWc5ms3v3GQwGk8kkfDuZTAaDwfaLBgDwf9obqsqyLIpiNBrdu2fYZ7FYlGVZluV0Os2yTHsVALBLf2q6AHfaJE5VxuNx9XWe51mWbdJjCAAQS3tD1WKxyLKsLMt701VomloSohUAwG60N1RtrspPoXUqdPwJVQDALvUhVFWqNq3xeHxrqDKAHQB2KfQ7JSJmqKqGMTXVSlSNVQ9PAt46Vj2powsADUqtLSPC039lWeZ5PhgMRt8MBoPBYNDIUPEwT9VwOKxmWAAA2IHHhqow68FsNhsOh+NvhsNhlmWj0WgH8xqESLcU4AyoAgB27FHdf6GjbTgc3toolef5ZDIJs5w/5n9ZL8/z2WwWolW9YNv7HwEAVj2qpSoMDL8rwYTt3zXd1IbCFOr1ZrDJZFIVoyiK2WxWn7kKAGDbHtv9tz67hJ9uqd2o+rVh7Hk1lis0nplRHQDYpcGDn4YL03JOp9P1vXuDweDefaLY5NnDweDh9QUAvktqt93tzlO1y7FNBqcDAA16ePffJkvshZ+KOwBA7z2qXS5M6rXmN9y7w46l1g4JAA1K7bb7qIHqYSXju+b5DA1Ut6523KDBWk2XDoA+WySm6b/3rj02QhZFEeYuHw6H9YWNZ7NZlmXj8bhVT+GlFpkBaJXH34PuHVez5rGt6u6c/fGuvT1PnjxJ67YbJYqGKdTrhsPhdDqN8ssjilVfAHiA3x+nmsZoPB7ftU91I65vPDs7uzUDnJ2dLb387OysuqcPh8PVHb5LarfdCGv/ZVlWlmX4dVWQWpriHACI5a71bW8djVOW5Zs3b7IsG4/HISSdnZ2FfPbmzZv6S8KeVWvWbDZ78+bNNibx7qs4oaqyFKTKsmxV9x8A9MOt+enWsBU2hiAVbtN5noe1epdeErJXvXUqy7KwFtwWatBDkUPVkrA44Fb/CwBISghDVXtS3Ww2Wx2QE/Zc7T5a+j0hOS2NtQr9hrf+X6zabqgCALZhtc0ibPmupW/rw62qJ8xWd9NStSGhCgBa4cmTTW/Kty6te9eDgaHt6smTJ7f2HVX7h8HvSy8PLzFIekNCFQA0ZvMgVbc6HCq7o+8vq00YOZlMnjx5EtLVvQuiTCaT0WgUfud3tX6lbLtr/wEAa/z+++9PnjypT4Ww9O0a9aFO6/v+fv/995CTwkuqNFYNV18SRqxX+2xSGDItVQDQrJCrwtebJ6qlHsB7JwXN83w6ndbnU8iyLDRH3Vqk+rQLnjnbUIS1/+71mP8iLjOqA9CgNfegkKvWJ6rJZDKZTKrmpSdPnoSptle/vvdXZbX5q87OztaksQ1/212vTeq2+6juv1v7bltufRBM6tgD0JS7hlLdun1NoKl3593aTxdy2K2xKcxWFfoE8zy/t62Lez0qVHXxGUuxCYDGLeWk0Ov3vW1CIRKVZbkmD4Xmj8lkcm9aCunqYS1SBMZUAUCT6uOo6uOr7lU9A3jXc3/Zt6R116zo9SausOfS8Kkutp40yNN/ANCY1ZHp39tWtGbSzuDs7OzNmzdv3rwJs6WH+DWbzZY6DUO7V9i4tM9d6zGzRKgCgMY8srstJKHsvuf+Qq6azWZLC84sTakQdltqrFo/jJ26tJ6G8/QfAA1q9h5UlmUVqpbW+Lt1tzX7bCi1p//SChlCFQANSu0elFqoMlAdACACoQoAIAKhCgAgAqEKACCC5KZUsEwNALANyYUqsQmA7rp3hb5qDvTVfTacUoEHS2uKAVMqANCgR96DqhnPlybtrKtWualPK1qW5Zs3b1Z3vnViz/q86mHR5QcX2JQKAECrLU16Xrl1qb4qUY3H47Ozs99///3s7CxEpTdv3tRfUpblkydPql8elqnZfC1C/KUAoHvWLJB868YQpEK7VGh/qtZjrvasstfv34RFAEej0Taq0D9CFQB0SQhDS6v4BbPZLMSgpY3ZbUOsln5PldLq/X3T6fSu/4tVQhUAtMJ3ZZfVRqlquNXmvyR0Bdb/99WXh5R2a8MYS4QqAGjSfD7/7//+7z//+c+j0ejPf/7z+/fv5/P5+peE6LMUdO56MDCkovpgqbpq/zDiajVU3dXQxSqhCgAaM5/Pj4+PLy8v//d///ff//732dnZ9fX127dv1+eq1eFQ2R19f9m3Lryw/5MnT0K6urXlaTU5hdFUt/5aVglVANCYT58+vXjx4h//+Mfr16/39/dfvXr1j3/84+XLlx8/frz3tfXuwvV9f6Gbr8pGk8nkzZs3d7VdVUajUfgvqljGekIVADTm06dPf/3rX58+fVpt+eGHH/76179++vRp/QuXegDvnRQ0z/PpdFqfTyHLsslkcuuTfaFNKzR91ee7Yr3kQtVgraZLB0ASnnxzdXX1n//5n0/+6D/+4z+ur6+rb2/9DUs9gHf1/a0K8ylUQ9Rns9ldU1WdnZ1po/ouyYWqxVpNlw6AJFQTQe3v7//rX//6/Y/+9a9/7e3tVd+u+T2he25N319oc7prBNXqrAr1qaoMTv9eyYUqAGiP4+Pjz58/39zcVFtubm7+53/+5/j4+N7XVj2Aa/r+QvPV+rFTlZCobn0GkE0IVQDQmL/97W8XFxc//fTT+fn59fX1+fn5f/3Xf3358uVvf/vbva+tegDX9P2FpLXUx1epN3GFHapZ13kAoQoAGvPs2bN//vOfe3t7o9Ho4OBgNBr9+OOP//znP58/f77hb7hr0s5KGDv15s2b0WgUJlMoy7Ja1K8+aD2rTbuw5JHVTMSfmi4AACTt2bNnf//73//+97+XZfm9rUTj8TiEofXP/Z2dnb1582Y2my1N2l6tAEgUg6RGZw8GadUXgFZp9h5UlmUVqobD4Q66+Z48eZLUbbcDISPP86Io7j32RVFUI/WKorh1H6EKgAaldg9KLVS1vZe0HqvXGAwG1aMNk8lkMBhY+hEA2KX2jqmqRtLdu2doxJpOp1Vr1mAwGI1GSaVjAKBZ7W2pCg8pbLJneJS03j9o2B0AsGPtDVVhivNNJsjfzWg7AIA12tv9t7nV4VNrmrgevMCfzkQAHskis/3W3paqhynLMpyydzVxrV/7z7KAAMDD9KGlqpLneXhUsD5oHQBgB3rSUhUaqGaz2Xg8XiwWEhUAsGN9aKkqy3I0Gg2HQ3NTAQBN6UNL1Wg0ym4brg4AsDOdbKkKTVPj8bhamibLstWlae5arAYAILpOhqpgqWlqdRoFoQoA2Jm0Fhi2oDIA7Exqt90+jKkCAGicUAUAEEGHx1Q9zPolApJqpQQAIkouVIlNAMA26P4DAIhAqAIAiECoAgCIQKgCAIgguYHq27b0dKFx8QCQCKEqspCiUptDFgDQ/QcAEIFQBQAQgVAFABCBUAUAEEFyA9Wt/QcAbENyoUpsAgC2QfcfAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABMnNqG6ZGgBgG5ILVWITALANuv8AACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACJIbpkaa/8BANuQXKgSmwCAbdD9BwAQgVAFABCBUAUAEEFyY6p4gKXR/calAcCqLrVU5XleluU2dma9xWIRglT1BQCwpDOhqizL2Wy2jZ0BAB6vA91/ZVmWZTmZTKLvDAAQSwdC1Wg02tLOAACxdKD7L4zjmU6n0XcGAIilAy1Vca1fpmYNA7QBgDWSC1WyEQCwDR3o/gMAaD+hCgAgAqEKACACoQoAIIJuh6qyLAeDQVEUTRcEAEhdt0NVYI0/AKBxg6SmGBgMdlTfnf1Hu9TLSgGwPandOPrQUgUA0DihCgAgAqEKACCC5JapWb/2X1JdvwBARMmFKrGJNC19nPBGAIguuVAFaQopKrUncQB2yZgq7jGfz9+/f39wcJBl2cHBwc8//zyfz5suFAC0jlDFOvP5/O3bt5eXl2dnZ5eXl9Pp9Pr6+ujoSK4CgCVp9QWY/PN7vX///vLy8sOHD0+fPg1bbm5u3r17t7e398svvzRbNh6gN2cm0AmpXXMSq61Q9Z0ODg7Ozs5ev35d33h+fp7n+dXVVVOl4sF6c2YCnZDaNUf3H3e6vr6+uro6PDxc2n54eBh+1EipAKCdhCrutLe3t7+/f3FxsbT969ev4UdNFAqgeYM/aro4tIVQxTrHx8enp6c3NzfVlpubm8+fPx8fHzdYKoBmLRaL0KtVfQGZUBVX/2YfODk5ubi4ePfu3fn5+fX19fn5+U8//fTly5eTk5OmiwYA7ZJcqBqs9Zjf3MvZB549e3Z6erq3t5fn+f7+fp7nP/744+np6fPnz5suGgC0S1rD8rf6GELvZx9I7SGOXnIQISJvqHul9idKrLbbPLq9n30gtfdGLzmIEJE31L1S+xMl1/23JWYfAIDECVVxmH0AABInVEVj9gEASFlanZ1b7dydz+dHR0cvX748Ojo6PDz8+vXr58+fv3z50ptn5VLrGu8lBxEi8oa6V2p/Ii1V0Zh9AABSllaEtKDyY/SyUqlxECEib6h7pfYn0lIFABCBUAUAEMGfmi7Arq1fiyapVkoAIKLkQpXYBABsg+4/AIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIkpunCu61NEOsuc0A2IRQxf2qkBG+6H3ICBVMbR1QYBOuDKwhVHE/V5BOcw+AiBaLxdJ7yluMSnKhytp/pMY9AOKqv6c6924yvGGrkgtVTiAS1Ol7ALRQeE9lHbynuA5sVXKhCtLU3XsAtMdqX0d9izcXQhX0mXsARFS9ZUJLjw8qLDFPFfyf9UPuumjxTVa79Nc3At+r3ndWRSvIhCqoW70+9uNy6R4AsayORvL5hEqfQ1We52VZNl0KOqaeOfoxltM9ACLq9NtnPp+/f//+4OAgy7KDg4Off/55Pp83Xahe6W2oKstyNps1XQo6qWdDJfpRC+CR5vP527dvLy8vz87OLi8vp9Pp9fX10dGRXBVRD0NVWZZFUYxGo6YLQscMau7aAtBRHz9+fPHixYcPH16/fr2/v//q1asPHz68fPny119/bbpo/dGH3o0l9fvfdDrN87z+o93Utx/dRsnq8XM9zkyIqFtvqIODg7Ozs9evX9c3np+f53l+dXW1pf+0W3+ix+thS1V4rGk6nTZdEDrJmG6gf66vr6+urg4PD5e2Hx4ehh81Uqr+6WGoWm/wUE0XnF0wphsexgWz5fb29vb39y8uLpa2f/36NfyoiUL1UHKhavFQTRecXXCg4WHq06F5H7XT8fHx6enpzc1NteXm5ubz58/Hx8cNlqpnkgtVAJCgk5OTi4uLd+/enZ+fX19fn5+f//TTT1++fDk5OWm6aP0hVAFA/z179uz09HRvby/P8/39/TzPf/zxx9PT0+fPnzddtP7o7bD8sixHo5Gn/3iwXh7EXlaK9kjtBOtufd0Nt0RLFQBABEIVAEAEQhUAQARpdXbqRWZzvTyIvawU7dGzE2xpwq3VqnWxvvdWKvp/17k/0WP8qekCQDPWXFnqa/9lJq+CVIX3fs9iQZ/q0kJCFYlac7l00QHgAYypAgCIILmWqvWLUmmigA7Z8egQgPWSC1Uuu9AbvRzyAnSX7j8AgAiEKgCACJLr/gO6xcApoCuEKqDVDJwCukKoAgC0CkcgVAEAWoUjMFAdAL7PYDAI7TrVF53Wgyq0hJYqAPg+PWvLWSwWSw1U2qseRksV0D0+WENcIVeFryWqB0suVA3Warp0wEbqN4DA+5e45vP5+/fvDw4Osiw7ODj4+eef5/N504XaruptJVE9WHKharFW06UDNtXLD9Y+5rXEfD5/+/bt5eXl2dnZ5eXldDq9vr4+OjrqZa5aPeWchA+WXKiCLMnPoL3Uvw/W1ac7H/Oa9fHjxxcvXnz48OH169f7+/uvXr368OHDy5cvf/3116aLFl+9WaE667Q1PExPPt5taGcfZ3vzubmXwmfQFy9eHB0dHR4eXlxcnJ6ehn+fPXvWdOm2q7tnZr3k6z89d7SCdf04TPWN9W87UbWDg4Ozs7PXr1/XN56fn+d5fnV11VSptqo6duGLWCdhd0/mh0mstkIVWfb+/fvLy8sPHz48ffo0bLm5uXn37t3e3t4vv/zSbNm2rbtn5l136/41VmW9O0z3/qhtrq+v9/f3Ly8v9/f3N9neA0sfWiIeqQ4d9ygSq61QRZKfQSvdPTNXS76lD9Zt0N269CNUZVl2cHAwnU5fvXpV3+gq0fhvaz9jqkjL9fX11dXV4eHh0vbDw8Pwo0ZKtQP9m6tw6Uqd1IWbbTs+Pj49Pb25uam23NzcfP78+fj4uMFS0X5CFWnZ29vb39+/uLhY2v7169fwoyYKFc2aZ8d69qBrD6pAm52cnFxcXLx79+78/Pz6+vr8/Pynn3768uXLyclJ00Xbop599GqEUEVyevwZ1LNjEMWzZ89OT0/39vbyPN/f38/z/Mcffzw9PX3+/HnTRduinn30akRanZ3GVJFl2Xw+Pzo6evnyZXj67+vXr58/f/7y5Utvrpi9PP16M1hnE92tUS8PU3dL3gap/fW0VJGcND+DArBtaUVILVUs6eWRSq1S/atv52o0n88/fvz46dOnq6ur/f394+Pjk5OTpVnfOlepSndL3gap/fWSa6karBXr92cG+gFpSGpFF1gvuVC1WCv673/8L4SUWVCo/ZJa0QXWS6tdLrV2SO7Vy1OiN5XacEGh3tS30q0abTibbrcqVdfdkrdBan+9xGqb2NHlXr08JXpTqQ0XFOpNfSsdqtH6FV2y2oxiHarUku6WvA1S++slVtvEji736uUp0ZtK9b4J5C7dqtH6FV16sJpQF8vcHqn99ZIbUwV0wiYLCnkupA3Wz6a7+LbidZbe/ZUECVVAG22yoJDnQtrg3hVdqlzlGNF7QhXQUj1eUKhP7ppN98WLF0uz1cSdvwZaKK3GWI3PLOnlKdGbSvV+QaG7dPcI3lryajRV1s3Gqu4ejjZI7a+npQpoKQsK9UD9nlofXwW9lFaETC0yc69enhIq1XXdrexSyevfdrFSSxGwc+Vvgy4e98f4U9MF2LX1n5OSOvYAW9X1K2rXy8/uJReqvEkAgG0wpgoAIAKhCgAgguS6/yCoz52T6RcG4NGEKhIlRQEQV6u7/4qiyPM8z/OiKOLuDAAQV3snkMjzfDabVd8Oh8OyLG/dsyzL0WgU9gkvuWvn1CbMIE29PM97Wam7dLeya0re3UrxGKkd95a2VBVFMZvNxuNxWCd1PB7PZrO7QlVIVIvFoizLxWIxnU5ns5n2KgBgl1oaIVfHDg8GgzXtT+PxuJ6iQivXXUtQbaPA0B69PM97Wam7dLeyWqpYktpxb2lLVZZlw+Fw6dt6b2AlxKw8z3dSKACA27U3VG2Yk8JuSy1Yt8Yv6Lf5fP7+/fuDg4Msyw4ODn7++ef5fN50oSAbDAah86H6AvqqjaHq1j6+9RlrMplUr1o/mmrwUA+tDezCfD5/+/bt5eXl2dnZ5eXldDq9vr4+OjqSq2jc4o+aLg5sURtD1ff25U2n0yzLRqNRSD+TyWSp67Bu8VCPrBRs1cePH1+8ePHhw4fXr1/v7++/evXqw4cPL1++/PXXX5suGkAqOjP5512P/mVZluf5YrGoGqjChFU7KRS0xadPn87Ozp4+fVpt+eGHH46OjvI8/+WXXxosGEA62huq1qSoW5lDgWRdX19fXV0dHh4ubT88PAw/2t/fb6RgAElpY/dfdtuzfmHaqlt3HgwGS01Ta3aG/tnb29vf37+4uFja/vXr1/CjJgoFkJyWhqrQ7FRFpfBF1RZVluVgMKi+DVODhm/DjzINVyTm+Pj49PT05uam2nJzc/P58+fj4+MGSwWQlPbOylUUxWQyqb6dTqdVxgrr0tTnAl1a06a+c11qs5CRjvl8fnR09PLly6Ojo8PDw69fv37+/PnLly+np6fPnz9vunQRJPXm7W5lu1tytiS1U6Lttd18bs+yLMuyDGsq37VPakeXpMzn819//fXTp0/X19d7e3vHx8cnJyc9SFRLE5qk8Bbu7pWquyVnS1I7JRKrbWJHlzQ5z7uuu0ewuyVnS1I7JVo6pgoAoFuEKgCACNo7T9WWrF9wJqlWSgAgouRCldgUV4KDiAHgVsmFqmb1L4KEKqQ2FBEAVglVOyWCAEBfGagOABCBUAUAEIHuPwAeqxowGr4wwoE0CVUAPJYUBZnuPwCAKIQqAIAIhCoAgAiSG1NlmRoAYBuSC1ViE9BC5gSGHtD9tyPrW8iAxC0Wi6WrhIsGdI5QtSOumMB69auEhivoIqFqd1wxgfWqq4TrA3RRcmOqmuWKCaxabbeub3G5gK4QqnbBFRNYo7oIhDZsH72go4SqXXDFBO5VHxUQLhSuEtAtxlTtzuoVs9nyAO2xGqEkKugcoWpHXDGBNVwQoAeEqh1xxQSAfjOmCvqjPmdHJsoD7FZyocraf7EYRdtCjghAg5ILVe46saw+nSRmAZAyY6p4OHPEA0BFqOJRTLsFAEFy3X9EYY54AFgiVPEQ5ogHgCW6/3g4c8QDQEWo4oHMEQ8AdUIVDyRCAUCdMVU7ZcJrAOgroWqnpChgafShywL0RnKhyjI1QLPCdcZ8udA/yYUqVzEAYBsMVAcAiECoAgCIQKgCAIhAqAIAiKDVoaooijzP8zwviiLuzgAAcbX3md48z2ezWfXtcDgsy/KuncNECcPhMMuy8Kpb6+UZ5i3xh4UNzefzjx8/fvr06erqan9///j4+OTk5NmzZ/V9vKHojdRO5pa2VBVFMZvNxuPxYrFYLBbj8Xg2m90VqvI8z7JsOp2WZVmW5XQ6Db9hd8UF2MB8Pn/79u3l5eXZ2dnl5eV0Or2+vj46OprP500XDYigpRFydRWXwWBwV2PV5junFpl3xh8WNvH+/fvLy8sPHz48ffo0bLm5uXn37t3e3t4vv/xS7eYNRW+kdjK3tKUq+9aXV/+23hu4Zk+Advr06dPR0VGVqLIs++GHH46Ojj59+tRgqYBY2huqQqfeJkJP32AwKIqiKIrQcLVmABbA7l1fX19dXR0eHi5tPzw8DD9qpFRARG0MVbfmoTUZK8/z0Fg1mUwmk0mWZePx+K6dBw/1+HoBKdvb29vf37+4uFja/vXr1/CjJgoFxNTGULV5G1W1/9Ko9slkctcvWTzU4+sFJO74+Pj09PTm5qbacnNz8/nz5+Pj4wZLBcTSxlB1qzXdeSFRVY/7FUWxZgAWQFNOTk4uLi7evXt3fn5+fX19fn7+008/ffny5eTkpOmiARG0N1RtOCgq7LbULhW+NawKaJVnz56dnp7u7e3leb6/v5/n+Y8//nh6evr8+fOmiwZE0NJnHUOP3tIsCfXmqLrVH62+vNqznfXtOn9Y+F5r3jXeUPRGaidzS1uqQkKq2p/CF1VsKssyPOsXvh0Oh5PJpN79N5vNzLMAAOxSeyNkURThUb5gOp1WGassy9FoVJ/ec8M1bVKLzDvjDwvfS0sVKUjtZG57bW8dMvXgnVM7ujuwNNmEPy9sSKgiBamdzInVNrGjC7SWUEUKUjuZWzqmCgCgW4QqAIAI/tR0AXZt/YIzSbVSAo2orkLhC5cd6I3kQpXrF9AsVyHoK91/AAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARCFUAABEIVQAAEQhVAAARJLdMjbX/AIBtSC5UiU0AwDbo/gMAiECoAgCIQKgCAIhAqAIAiECoAgCIQKgCAIhAqAIAiCC5eaoAWquanTh8YV496BahCqAtpCjotORClWVqAIBtSC5UiU0AwDYYqA4AEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABBBcjOqW6YGANiG5EKV2AQAbIPuPwCACIQqAIAIhCoAgAiEKgCACFodqoqiyPM8z/OiKNbsNrhDWZY7KigAkLz2Pv2X5/lsNgtfz2azsiyFJACgtVraUlUUxWw2G4/Hi8VisViMx+OQq27debFiOBwOh8M8z3daaAAgYYN2ztsUpuisl20wGAyHw00aq4qimEwmt9ZrMGhpfQGgf1K77ba0pSrLsuFwuPRt1Ru43mQyGY/H2ykUAMDt2huqHtZ5F161fmA7AEB0bQxVt/bxbZixZrPZdDpds8Ndjwre60FVAQBS0can/x48wDy8cP3Lk+rcBQB2po0tVbfaZIh6eGBw+2UBAFjW3lD1vbNShXFURlMBAI1o6bOOYebPpSkVxuPxmsy0OgvDrfu0s74A0D+p3XZb2lIVwlM1Omrpmb6yLAeDwWrA0vcHADSlpaEqz/Mwi3p48u7WZ/rq/YPha1OoAwBNaXu7XNy0lFo7JAA0KLXbbmK1TezoAkCDUrvttrT7DwCgW4QqAIAI2jij+latX3AmqcRVwB4AAAm7SURBVFZKACCi5EKV2AQAbIPuPwCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACIQqAIAIhCoAgAiEKgCACJJbpsbafwDANiQXqsQmAGAbdP8BAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEIFQBAEQgVAEARCBUAQBEkNyM6papAQC2IblQJTYBANug+w8AIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIAKhCgAgAqEKACACoQoAIILklqmx9h8AsA3JhSqxCQDYBt1/AAARCFUAABEIVQAAEbQ6VBVFked5nudFUdy7c1mWm+8MABDXoLUDt/M8n81m1bfD4bAsy7t2LopiMpmE3cKrbq3XYNDe+gJAz6R2221pS1VRFLPZbDweLxaLxWIxHo9ns9ldoaosy8lkEnYuy3I6nWZZluf5LgsMACSupREyzCZVL9tgMLirsSq0adV3Dj2Aq/2AqUVmAGhQarfdltZ2NUKtJqc1O6/5te2sLwD0T2q33ZZ2/2Xf2X8XxqfneT4YDPI83yRgAQBE1MZQdVcf35qdJ5NJGKgeRl+NRqO7ctXgoeLUDQDoqTYuU/OwMeZVA2NRFIPBYDQa3drkmFQ7JACwM21sqbrVXS1PIYENh8P6xvF4vP0SAQD8n/aGqsePizKyCgDYmZaGqmoOz0qYtmrDnUOcMlUVALAzLQ1VYYqpKhWFL6p5p8qyHAwG1bdLO1cTh+6qsAAArRyonmVZnufj8XgymVSP3YV50uuq3r3VnYfDoRUAAYBdavusXN/VkXfvzqnNQgYADUrttptYbRM7ugDQoNRuuy0dUwUA0C1CFQBABEIVAEAELX36b3vWr+KXVNcvABBRcqFKbAIAtkH3HwBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQARCFQBABEIVAEAEQhUAQATJzahumRoAYBuSC1ViEwCwDbr/AAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIkhuRnXL1AAA25BcqBKbAIBt0P0HABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQBQAQgVAFABCBUAUAEIFQ1Xnr54jvqP5Vqn81ylSqO/pXqf7VKOtppVIjVAEARJDcMjXW/gMAtiG5UCU2AQDboPsPACACoQoAIAKhCgAgglaHqqIo8jzP87woivV7Dlbc+xIAgIjaO1A9z/PZbBa+ns1mZVmWZXnrnndtBwDYmZa2VBVFMZvNxuPxYrFYLBbj8TjkqjUvmU6nixotVQDALg3aOcVAmE2qXrbBYDAcDm/NVUVRTCaTTSoyGLS0vo+hUp3QvxplKtUd/atU/2qUqVQvtLSlKsuy4XC49G3VG7ikSlpruggBALaq1WOqvmv/+lTpd7VpAQBsy6J9ptNplmXVgKpgPB7fVdpQkWr/0MS19PL6ngDAbsRNCC3Xxpaq722jmk6n9ZeUZZnn+WQyWR2rvpCrAIDtaO+YqiVruvNWQ1jYogcQANiZ9oaqx0ei723xAgB4sJaGqtVn/cK0Vat7lmW5On+6NioAYMdaGqpCSKqamsIXVXKqB6k8z4fD4WQyqYJUNXHoTksMAKStjQPVsyzL83w8Hk8mk2qihPBIYF19eqrBYDAajaofjcdjM6oDALv0/1obPsI6ynme/+Uvf/ntt99evnxZ/ejly5dFUfzlL3+ptoQ9w/bffvutnaOpwp+6qkhZlvVKdVG9RkVRhJirUi20dO5l394yDRUnjqpS4TBVuluv+mEqy/K3337L+nXu9fIN1fVK3XtjChXs7ttqx1ra/VfJ83zDY1mFsO0W6BGWRnqNRqOuj/2qyp/nefh6NBq1+RBsoqrUYDCoKtXazx4bWjrTwspODZUlmqpS9d7/TquPYQhN730693r5hqpXqqMn4fobU6hg6A7accG6qumJsrokrNk8Ho+rxZvDF9PptNoyHo+X5h2dTqfhJcPhMOxWrcCztAj07kWpUeiZrX7a+EkVvVLj8Xg4HO6i6HeLde4tvk2u2/hhWsSrVBvqEkSvUdi4/YKvE/0NtXTFaESUStWnpG6wUrHOutUbU/3SNxwOb51SmyVtuRi1X3UrCm+kcHoNh8NwIlaXwmpLeFW182qKajxURazRUr12XZOa6Idp0YJQFbdS2bcRio3UpRKxUmFLPTg2IlaNqiDVbHWCuIepunn34w21dGVo5D211RtTPUg1fhnsCqFqU9Pb2mPqV4dbz7+ll7QwVFXfPr5Gi2+3ty0W+j5xKzVtR6NOxEoNW9O6E6tS02+rWoU7RIPX/Vg1qm51w7tX3NqZiOfeuPZE9k7Kfqe45161W4Ohql6wRbwb09LFUKjaRNvHVLVZ6HiuhhDNZrP6HBCz2Wypr3r4rXG1tR5To6IoBoPBuH3PXT6mUnmeh0tJ2waKPaxSSzOVtM3DKpXn+WKxKIqiKIrFYrE0v12zHnPuLRaLsiyn02nbRr89rFJlWU4mkypdte0kfPC5Nx6PB4NBq8YbRbwxtf+e1UJCVTT3nn+tutxvYvMahYHq4d629WI9zoaVCmMzs29Va/mx27xSs9msugdUY2zbKdk3VD1wtC18rNr83Ks+nISpBHdQtgfb/NwLUb7Nl75Hvo/qUxdFKlHPCVXfJ5xYSx8FgvD4Yfi6KIrq01jLz8XH1yjcqltVzSiVqmY+a0nVolSqaqPOsmyxWDR+z358pUITaf23NevxNaq/qg01yiJVqrp596ZSWZZV515RFA2262zpxlTNFpFl2WQyaW1wbJdd9zd2VjW8pj7QYfjHByJW/7DVOI/wbwvHVD2+RvWhEo2fVBEPU323fhyp+u9s/L0fsVLVt80eqYg1qm/szbnXksO02Fqldl+RxfZvTGF4VtboaMVuEao2teEjOWH04urGLZXqMfpXo0XsSt262+45Uqsbb90zTskeqn81WqhU668SO7g4tOGQdYhQtalp75596F+NFirVHf2rVP9qtFCp1utTXfrBmKrv0PgAlOj6V6NMpbqjf5XqX40ylWq9PtWlBwaLP3a4AgDwAFqqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACIQqgAAIhCqAAAiEKoAACL4/+qvr5PQT+p7AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "g['lyNorm_MS3'].Draw(\"AP*\")\n",
    "g['lyNorm_MS3'].GetHistogram().GetYaxis().SetLimits(0.5,1.3)\n",
    "g['lyNorm_MS3'].GetHistogram().GetYaxis().SetRangeUser(0.5,1.3)\n",
    "g['lyNorm_MS3'].GetHistogram().GetYaxis().SetTitle(\"LO\")\n",
    "g['lyNorm_MS2'].SetMarkerStyle(24)\n",
    "g['lyNorm_MS2'].SetMarkerSize(1.2)\n",
    "g['lyNorm_MS2'].Draw(\"PSAME\")\n",
    "l.Draw()\n",
    "c1.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bb240a08-4871-494b-bd90-f32eda7a7ad6",
   "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.9.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}