parallel.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
       "    return false;\n",
       "}"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%%javascript\n",
    "IPython.OutputArea.prototype._should_scroll = function(lines) {\n",
    "    return false;\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Linux vinavx3.cern.ch 3.10.0-514.10.2.el7.x86_64 #1 SMP Fri Mar 3 00:04:05 UTC 2017 x86_64 x86_64 x86_64 GNU/Linux\n",
      "innocent\n",
      "root://eoscms.cern.ch///store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root\n",
      "ERR Could not open file root://eoscms.cern.ch///store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root\n",
      "root://cms-xrd-tzero.cern.ch///store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/061/00000/5C4984B0-B86A-E711-837F-02163E01A792.root\n",
      "root://cms-xrd-tzero.cern.ch///store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/061/00000/5C4984B0-B86A-E711-837F-02163E01A792.root (1 runs, 2 lumis, 9642 events, 3107653772 bytes)\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "ls: cannot access /eos/cms/tier0/store/data/Run2017: Permission denied\n",
      "Error in <TNetXNGFile::Open>: [ERROR] Server responded with an error: [3011] Unable to open file /eos/cms/store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root; No such file or directory\n",
      "\n"
     ]
    }
   ],
   "source": [
    "%%bash\n",
    "uname -a\n",
    "whoami\n",
    "ls /eos/cms/tier0/store/data/Run2017\n",
    "edmFileUtil root://eoscms.cern.ch///store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root\n",
    "edmFileUtil root://cms-xrd-tzero.cern.ch///store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/061/00000/5C4984B0-B86A-E711-837F-02163E01A792.root\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.10/01\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "import math\n",
    "import time\n",
    "from multiprocessing import Pool\n",
    "from DataFormats.FWLite import Handle, Events"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "xrd = 'root://cms-xrd-global.cern.ch//'\n",
    "tier0 = 'root://cms-xrd-tzero.cern.ch//' # 'file:/eos/cms/tier0'\n",
    "iter2 = 'file:/eos/cms'\n",
    "\n",
    "#LS50\n",
    "events2016HL = xrd+\"/store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root\"\n",
    "#LS1000\n",
    "events2016LL = xrd+\"/store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/88FB6BA5-7B98-E611-B25E-FA163EE1CC08.root\"\n",
    "\n",
    "events2017B_HL = tier0+\"/store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/061/00000/5C4984B0-B86A-E711-837F-02163E01A792.root\"\n",
    "events2017B_LL = tier0+\"/store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/067/00000/B4146406-F66A-E711-B9C9-02163E014389.root\"\n",
    "\n",
    "label = \"generalTracks\"\n",
    "quality = \"highPurity\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def nt(fevents) :\n",
    "    print fevents\n",
    "    events = Events(fevents)\n",
    "    tracks = Handle(\"std::vector<reco::Track>\")\n",
    "    lumi = Handle(\"std::vector<LumiScalers>\")\n",
    "    nt = []\n",
    "    for i in range(0,1000) : # events.size()):\n",
    "      a= events.to(i)\n",
    "      if (i%500==0) : print \"Event\", i \n",
    "      a=events.getByLabel(label, tracks)\n",
    "      nt.append(tracks.product().size())\n",
    "    print len(nt)\n",
    "    id = events.object().id()\n",
    "    a=events.getByLabel('scalersRawToDigi', lumi)\n",
    "    runid = '{:d}:{:d}:{:.2e}'.format(int(id.run()),int(id.luminosityBlock()),lumi.product()[0].instantLumi())\n",
    "    print runid\n",
    "    return (runid,nt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "root://cms-xrd-tzero.cern.ch///store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/061/00000/5C4984B0-B86A-E711-837F-02163E01A792.root\n",
      "root://cms-xrd-global.cern.ch///store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/AE49FDAB-8398-E611-9B98-02163E014388.root\n",
      "root://cms-xrd-tzero.cern.ch///store/data/Run2017B/SingleMuon/AOD/PromptReco-v2/000/299/067/00000/B4146406-F66A-E711-B9C9-02163E014389.root\n",
      "root://cms-xrd-global.cern.ch///store/data/Run2016H/JetHT/AOD/PromptReco-v2/000/283/408/00000/88FB6BA5-7B98-E611-B25E-FA163EE1CC08.root\n",
      "Event 0\n",
      "Event 0\n",
      "Event 0\n",
      "Event 500\n",
      "Event 0\n",
      "1000\n",
      "299061:71:1.20e+04\n",
      "Event 500\n",
      "1000\n",
      "299067:23:7.98e+03\n",
      "Event 500\n",
      "Event 500\n",
      "1000\n",
      "283408:999:8.97e+03\n",
      "1000\n",
      "283408:50:1.43e+04\n"
     ]
    }
   ],
   "source": [
    "p = Pool(4)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "78.6919720173\n",
      "283408:50:1.43e+04 1000\n",
      "283408:999:8.97e+03 1000\n",
      "299061:71:1.20e+04 1000\n",
      "299067:23:7.98e+03 1000\n"
     ]
    }
   ],
   "source": [
    "t = -time.time()\n",
    "ret = p.map(nt, [events2016HL,events2016LL,events2017B_HL,events2017B_LL])\n",
    "#ret = map(nt, [events2016HL,events2016LL,events2017B_HL,events2017B_LL])\n",
    "t +=time.time()\n",
    "print t\n",
    "for n in ret : print n[0],len(n[1])\n",
    "runids,nts = zip(*ret)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "p.terminate()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def plotNT(nt):\n",
    "    y,x = np.histogram(nt,np.linspace(0.,3000.,300))\n",
    "    x.resize(len(y))\n",
    "    return (y,x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "4\n",
      "4 299\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAEACAYAAAC9Gb03AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXt8FdW1x387BsnrJJkJr4BIAopUEYNKW0mUIFaiXoig\nBaUgQbRC8f24gFhILFrQWkGlt4pWwF70FkEhtEWLcMBApQqEl8hLCfJGZkIeBMxj3T/OzGTOMyc5\nJ6/N+n4+B87svWfvvfZMVibr/M7agojAMAzDyEtEc0+AYRiGaVzY0TMMw0gOO3qGYRjJYUfPMAwj\nOezoGYZhJIcdPcMwjOTU6eiFEG2FEJuEEFuFELuEEC8a5YoQ4lMhxB4hxCdCiIQmmTHDMAxTL0Qw\nOnohRAwRnRVCXARgA4CnAAwFcJqIXhJCTAagENGUJpk1wzAMEzRBhW6I6Kzxtq1xjg4gG8BCo3wh\ngDsbb5oMwzBMQwnK0QshIoQQWwEcB+Akoq8BdCSiE3D9IjgOoEOjz5ZhGIapN5HBNCKiGgB9hRDx\nAD4RQmQC8Iz5cC4FhmGYFkhQjt6EiEqEEP8AcD2AE0KIjkR0QgjRCcBJX+cIIfgXAMMwTAMgIhGO\nfoJR3bQzFTVCiGgAvwCwFcAKADlGs7EAlgeYrLSvGTNmNPsc2D62je2T7xVOgnmiTwawUAghjF8M\n7xHRZ0bM/m9CiPsBFAEYEdaZtRIOHjzY3FNoVGS2T2bbwPYxNup09ES0A8C1Pso1ALc02swYhmGY\nsMDfjA2RnJycIFq1XmS2T2bbwPYxNoL6wlRIAwhBjT0GwzCMbAgh0GQfxjKBcTqdzT2FRqWl2peS\nkgIhBL/41epfKSkpjf7zUi95JcO0FIqKisKuTGCY5kCIsDy0Bx6DQzdMa0S4/qxt7mkwTMj4u5cF\nh24YhmGYYGFHHyItNYYdLmS3j2EuBNjRMwzDSA47+hDJzMxs7ik0KrLb1xj8+OOPeOCBB5CSkoKE\nhARce+21WLVqlVX/8ccfo3fv3khISEDv3r2xfHlt9pA5c+agR48eiI+PR6dOnXD//fejrKzMa4x1\n69YhIiIC06dPdytfvHgxUlJS4HA4MHz4cBQXF1t1paWlGD16NNq3b48OHTpgzJgxPvs2iYiIgMPh\ngMPhQHx8PH7961+71b/66qtITk5GYmIiHnjgAVRWVvrta968eejXrx+ioqJw//33B7GKLgYNGoSI\niAjU1NRYZWPGjEFycjISEhLQo0cPvPDCC0H3Fyy6rmPYsGGIi4tDamoq3n//fZ/tnn/+eURERGDN\nmjVhn0NYaYJ8DcQw4aYl31fl5eWUl5dHhw4dIiKilStXksPhoKKiIjp58iTFxMTQJ598QkREf//7\n3ykmJoZOnTpFRETffvstaZpGRES6rtPNN99MU6ZMceu/srKS0tLS6IYbbqDf/va3VvnOnTvJ4XBQ\nQUEBlZeX06hRo+iee+6x6p944gkaPHgwlZWVUUlJCd1yyy301FNP+bUjIiKCvv32W591q1atok6d\nOtHu3bupuLiYMjMzaerUqX77+uijj2j58uX0m9/8hsaNGxfUOv7v//4v3XTTTRQREUHV1dVW+a5d\nu6iiooKIiPbs2UMdO3akVatWBdWnndzcXMrLy/NZd88999A999xDZ8+epYKCAkpISKCvv/7arc2B\nAwfo6quvpi5dutBnn31W7/FN/N3LRnlY/DA/0YdIOGPY6mwV6mw1bP2FA47R15+YmBhMnz4dXbt2\nBQDccccdSE1NxebNm7F//344HA7ceuutAIDbb78dsbGxOHDgAAAgNTUViqIAAKqrqxEREYHk5GS3\n/l955RUMHjwYvXr1citfvHgxhg4divT0dMTExOB3v/sdli1bhvLycgDArl27cOeddyI2NhYOhwPD\nhg3Drl27/NpBRG5P0nYWLVqE8ePHo1evXkhISMD06dPx7rvv+u3rzjvvxNChQ6Gqwd3fJSUleP75\n5/Hyyy971V155ZWIioqy5timTRu0b9/eql+5ciX69u0LRVGQkZGBHTt2BDWmydmzZ7Fs2TLMnDkT\n0dHRSE9PR3Z2Nt577z23dpMmTcJLL72ENm3a1Kv/5oAdfQtCP6dDP6c39zSYMHPixAns3bsXvXv3\nRp8+fRAZGYmVK1eipqYGH3/8MaKiotCnTx+r/fvvv4+EhAR06NABHTp0wKOPPmrVFRUV4d1338X0\n6dO9JHm7du3CNddcYx13794dbdu2xd69ewEAgwcPxrJly1BcXAxd17F06VLcfvvtVvtrrrkGH3zw\ngVufAwYMQOfOnXH33XejqKjI71jXXHMNTp48CV133b+TJk3Cww8/3OA1e/bZZ/Gb3/wGHTt29Fk/\nadIkxMbGonfv3pg2bRquvdaVjmvr1q0YP3485s+fD03T8NBDD2Ho0KEBw0qe7N27F23atEGPHj3c\n7LP/UlyyZAmioqKQlZXVYBubEnb0ISJ7DLu12idEeF6hUlVVhdGjR2PcuHG4/PLLERsbizfffBMj\nR45EVFQURo8ejTfffBPR0dHWOffeey/OnDmDvXv34uuvv8acOXOsusceewwzZ85ETEyM11hlZWVI\nSEhwK4uPj0dpaSkAWI43KSkJ7du3R2RkJCZOnGi13bZtG+655x7reP369Th48CC++eYbJCcn47/+\n67+sJ3zPseLj40FE1ljz5s3DG2+80aA1++qrr7Bx40Y88sgjftvMmzcPZWVlWL16NZ577jl8+eWX\nAID58+djwoQJuP766yGEwJgxY9C2bVt88cUXQY9fVlaG+Ph4tzL7OpaWlmLatGl47bXXGmRfc8CO\nnpESovC8QpsDYfTo0Wjbti1ef/11AMCWLVvw4IMPoqCgAD/++COcTifGjx+P7du3e53fo0cPTJky\nBYsWLQIA5Ofno7S0FHfffbfP8eLi4lBSUuJWdubMGTgcDgDAqFGj0LNnT5SXl6OkpATdu3fHr371\nK7/zz8jIQGRkJOLj4zF37lwcPHgQu3fv9jnWmTNnIISwxmooRIRJkyZh7ty5dX4pTgiBAQMG4Je/\n/KX1YWlRURFeeeUVqKoKVVWhKAoOHz6Mo0ePAgCGDBkCRVGgqipmzZqFWbNmWW2HDh0a1Drm5ubi\nvvvus0JzrQFOgRAiTqez1T71BoPs9jUm48ePxw8//IB//OMfuOiiiwAAa9aswQ033IC+ffsCAK6/\n/nr87Gc/w+rVq93CNyaVlZXW0/uaNWuwefNmK2Z/5swZREZGYseOHfjoo49w1VVXYdu2bda5Bw4c\nQGVlJa644goAwKpVq/Dvf//bim9PmDABN954Y1C2mA7X/N8cy/ylU1hYiI4dO1qfLzSUkpISbN68\nGSNHjgQRobq6GkSESy65BEuWLEF6errXOVVVVdYade3aFdOmTcPUqVN99p+fn2+9z8vLgxDCS7nU\ns2dPVFVV4cCBA1b4Ztu2bbjqqqsA4zocOXIE8+bNAwCcOnUKI0aMwOTJk/HMM8+EZH+jwaqb0Fi7\ndm3Y+kIuCLkta73CaV84aen31UMPPUQ33HADlZeXu5V/8skn1L59eyosLCQioi1btlBSUhKtXr2a\niIjefvttOnnyJJGhLrnqqqvoj3/8IxERlZWV0YkTJ6zXyJEj6cknnyRd1632CQkJVFBQQGVlZTRq\n1CgaNWqUNXb//v3p0UcfpYqKCjp79ixNnDiR0tPTfc5/165dVFhYSNXV1VRaWkqPPvoo9erVi6qq\nqogM1U1ycjJ9/fXXpGkaZWZm0rPPPut3PaqqqqiiooKmTp1KY8aMoXPnzll9eWK38csvvyQhBB07\ndowqKyvp5MmT9MEHH1BZWRlVV1fTqlWrKCEhgb788ksiIvrqq6/o0ksvpU2bNllr9ve//53Kysq8\nxgmkurn33ntp1KhRVF5eTp9//jklJiZaqhtN09zm2LVrV1q6dKnXtQ6WplDdsKNvQbRER99Sacn3\nVVFREQkhKDo6muLi4iguLo4cDgctXryYiIhefvll6t69OzkcDurRowe9+uqr1rnjxo2jjh07ksPh\noCuuuIJeeuklv+Pk5OS4ySuJiN5//3269NJLKS4ujoYNG2b9EiAi2rt3Lw0ePJhUVaWkpCS67bbb\naP/+/Vb9VVddZc1xzZo1dMUVV1BcXBx17NiRhg0b5taWiOjVV1+ljh07UkJCAo0fP55+/PFHq27C\nhAk0ceJE6zg3N5eEEBQREWG9TCd76NAhcjgc9P3333vZePDgQTd55alTp2jAgAGkKAolJiZSv379\naMWKFW7nfPLJJ9SvXz9SFIU6d+5MI0aMqLej1zSN7rzzToqNjaVu3brRBx984Pc6pKamtnh5JSc1\na0GIPNenfzSD16suOKkZIwuc1KwV0JQ68+bQ2bOOnmFaP/xhbCuCNfYMwzQEfqIPEdkVKbLbxzAX\nAuzoGYZhJIcdfYjIHsOW3T6GuRBgR88wDCM57OhDRPYYtuz2McyFAKtumghTFqlN1vw36r/cfx3D\nMEwD4Sf6EAk2hh1UCuI28a5XC4Jj9AzT+mFHzzBhRpatBPPz83H11VcjPj4eGRkZVuZK08YnnngC\nXbp0QVJSEh5++GFUV1db9d988w0GDRqExMRE9OzZEx9//HHANQu0Jp4E6nvx4sXW1ofx8fGIjY1F\nREQEtm7dGnD8+vLZZ5/hJz/5CeLi4jBo0CAcOnTIqgv2GjYpQeSquQTAGgC7AOwA8IhRPgPAYQBb\njFcW57rxTzB5bLB2LSFAEjHOhVNLS76vZNhKcN++fRQfH08bN26k6upq+v3vf0+XXXaZlXMmNzeX\nbrrpJiouLqYffviBfv7zn1Nubi6RkcCsZ8+eNGfOHKqpqaE1a9ZQbGws7du3z+dYda2Jnfr2vWDB\nArrssssCXC3/pKSkUFFRkVf5Dz/8QAkJCbR06VI6f/48PfPMM/Tzn//cqg/mGtppEUnNAHQCkGa8\njwOwB0Avw9E/yUnNgoMdfXhpbfdVnz59aNmyZbRx40bq2LGjW1379u3piy++8Drnhx9+oFtuuYXm\nzp3rVj5r1iyaPHkyjRs3zs3RP/vss/SrX/3KOj5w4ABdfPHFVkKvW2+9lf7nf/7Hqp83bx5lZWX5\nnO8bb7xBd9xxh3VcU1ND0dHRtGbNGiIiuv7662nJkiVW/eLFi+nSSy8lsv3CsXPrrbfS9OnTfY5V\nnzWpb98DBw6k559/3jo+f/48PfXUU3TppZdSp06daOLEiXTu3Dmf56ampvp09G+99ZZb1s/y8nKK\njo6mPXv2eLX1dw3ttIg9Y4noOBEVGu/LAOwG0MWoDkvCndaM7DFs2e1rClrzVoImNTU1ICLs3LnT\nb/3hw4etXZg88TxXURRs3LgRAIJak0D4m1dRURE+//xz3HfffVbZ5MmTsX//fmzfvh379+/HkSNH\n8Pzzz/vt1xee6xwTE4PLLrvMbavBQNewOahXjF4IkQIgDcAmo+hhIUShEOJtIURCHaczTJMh8kRY\nXqHSWrcSvOWWW7Bu3TqsX78elZWVePHFF1FZWYmzZ88CALKysjB37lz88MMPOH78uLWD1tmzZ3HF\nFVegQ4cO+MMf/oCqqip8+umnWLdunXUuAOi6jv79+wNAUGtiEkzfJosWLcKNN96Ibt26WWXz58/H\nq6++ioSEBMTGxmLKlCnW7lS+8OXs61pn1HENm4V65JWPA/AVgGzjuD1gpTmeCeAdDt34h0M34aU1\n3Fc1NTU0cuRIuuOOO6xNNjZv3kzJycm0ZcsWIiL68ssvKTk5mbZt2+azjw8++ID69u1LREQrVqyg\nm2++2arzzEefnZ1NL7/8stv5cXFx1lh33XUXTZo0iSoqKqi8vJwmTJhAI0aM8Dv/pUuXUu/evald\nu3b0+OOP09VXX01//etfiYiooqKCHnnkEerSpQv16NGDZs2aRW3btrXO3bFjBw0YMIDatWtHWVlZ\nNGbMGHrggQd8jlPfNQm278svv5wWLlxoHZ88eZKEEKQoivVKSEig+Ph4IiMvfmJiopXrPiIighIS\nEqyy999/n4iIHnvsMZo0aZLbWL1796Zly5b5nK/9GvqiKUI3QenohRCRAD4E8B4RLTe89ylbk/kA\n8v2dn5OTg5SUFABAYmIi0tLSrC/imKEB2Y9N6mqPwkI4A9V/5769X0uxr7nWsyXT2rcSHD58OIYP\nH26N9fbbb6Nfv34AgKioKLz22mvWBtlvvfUWrrvuOuvc3r17u12r9PR05OTk+BynvmsSTN8bNmzA\nsWPHcNddd1ll7dq1Q0xMDHbt2mWtoZ2uXbtC12sl0N27d8e6deu89oa96qqrsHDhQuu4vLwcBw4c\nsLYa9MR+DQPhdDqxYMECALD8ZdgI8ml+EYA/en5Ia3v/BIDFF+ITfbBb7bXWJ3reSrBhtPatBMl4\n0q6urqaTJ0/SiBEjaPTo0VbdkSNH6OjRo0RE9O9//5u6du1q2UBEtH37djp37hyVl5dbO2rZd6Cq\na03atWtH//rXv3y2D6bvBx98kMaOHet17uOPP04jRoyw1vjw4cOW2scTf6qbU6dOUWJiIi1btozO\nnTtHzzzzDN1www1WfaBr6IuWorpJB1ANoBDAVlNKaTj/7Ub5xwA6sqP3Dzv68NKS7ysZthIkIsrI\nyCCHw0FJSUk0ceJEOnv2rFW3fv16SklJodjYWOrVq5cV1jB55plnSFEUcjgcdPvtt9OBAwfc6uPi\n4qigoMA6DrQmL774It1+++1B933u3DlSFMXnvXv+/Hl69tlnqXv37pSQkEBXXnklvf766z7X15/q\nhojos88+o169elFMTAwNHDjQrV19riE1kaPnrQSbiGC2CRTGn6PkJ78MbzVYC28lyMgCbyXIMAzD\nhAw7+hBpTR8ONgTZ7WOYCwF29AzDMJLDaYpDJBz52s0UxvjZstAnFGY4Hz3DtH7Y0bcA6kxfzDAM\nEwIcugkR2WPYstvHMBcC7OgZhmEkhx19iMgew5bdPoa5EGBHzzAMIzns6ENE9hi27PY1BnVtJfj2\n22/j8ssvR3x8PG6//XYcO3bMqjtz5gxycnLQsWNHdOrUCXl5eW59FxUV4eabb0ZsbCyuvPJKfPbZ\nZ271P/zwA371q18hMTERSUlJGDNmjFW3ZMkSpKenIzY2FjfffHNAGwoKCty25HM4HIiIiMBHH30E\nGDnZs7Ky0L59eythWyDmzZuHfv36ISoqCvfff3/AtosWLcL111+PhIQEXHrppZg8eTJqamqsel3X\nMWzYMMTFxSE1NTVgmuGGEuwYzz//PCIiIrBmzZqwzyGcsKNnmDBTVVWFSy+9FJ9//jnOnDmD3/3u\ndxgxYgQOHToEp9OJadOmIT8/H5qmISUlBffee6917uOPP46KigocOnQImzZtwnvvveeWKfHee+/F\nddddB03TMHPmTNx99904ffq0VT98+HB07twZhw8fxsmTJ/H0009bdUlJSXjiiScwderUOm3IyMhA\naWkpSkpKUFJSgpUrV8LhcCArKwsA0KZNG4wcORJ/+ctfglqTLl264Le//S3Gjx9fZ9uKigrMnTsX\np0+fxqZNm/DZZ5/hD3/4g1X/m9/8BlFRUTh16hT++te/YuLEiW772QZLXl6e301Hghnj22+/xYcf\nfojOnTvXe+wmJ1xJcwIkRQuY0OdCIVBCMrOuJSY1a6m0tvvK3Erw6aefdstlfvToURJC0LfffktE\nRO3ataOvvvrKqn/xxRfppptuIiKiPXv2UFRUlLU1IBHRTTfdRG+++SaRkQUyNTWVampqAs7l7bff\npoEDB9Zr/jk5OXT//fd7le/fv58iIiKC7ue5556jcePG1WvsP/7xjzR06FAiY9u+iy++2C0Z2333\n3UdTp061jvPz8yktLY0SExMpPT2dtm/f7rPf3NxcysvL8yoPZgwioqysLPrnP/9JKSkp9Nlnn9XL\nJjstYitBhmFC48SJE9i3b5/PfOVmSMK+FZ49wVVNTY1V9/XXX6N79+6IjY216q+55hprC7tNmzah\nZ8+euO+++9CuXTv87Gc/w/r164Oep7+tBM+ePYulS5f6zSfvi0mTJlk7WoXK+vXrrbXbu3cv2rRp\ngx49erjN21yDrVu3Yvz48Zg/fz40TcNDDz2EoUOHorKyMujx6hoDRhgsKirK+gunpcOOPkRkj2G3\nWvuECM8rRMytBHNyctCzZ09kZWVhyZIl2LlzJyoqKqwYr32LvtmzZ6OsrAz79+/Hu+++a9XVtYXd\n4cOH8a9//QuDBg3CiRMn8OSTTyI7OxuapgU1V/tWgnaWLl2K9u3bB9ykxJN58+bhjTfeCLq9P/7y\nl79g8+bNVgiqrKwM8fHxbm3sazB//nxMmDAB119/PYQQGDNmDNq2bYsvvvgi6DHrGqO0tBTTpk2z\nNl1pDbCjZ+TEtdlC6K+QpkAYPXo02rZta+2pOmjQIOTm5mL48OHo3r07unfvDofDgUsuuQQA8Npr\nr6Ft27a4/PLLMWzYMIwaNcqqi4uLQ0lJidsYZ86cgcPhAABER0cjJSUFOTk5uOiiizBy5Eh07doV\nGzZsCMmORYsWuW2w3VR8/PHHmDZtGlatWgVVdaUJqWsNioqK8Morr0BVVaiqCkVRcPjwYRw9ehQA\nMGTIECiKAlVVMWvWLMyaNctqO3To0KDGyM3NxX333ee181RLhh19iMiuM5fdvsbE3Epw2bJlbsqU\niRMnYu/evTh27BiGDx+Oqqoq9O7dGwCgKAr++te/4tixY9ixYweqq6vx05/+FDC2sPv2229RXl5u\n9bVt2zYrrNGnTx8Ij79CPI/ry+HDh+F0Opvc0a9atQoPPfQQVq5ciSuvvNIq79mzJ6qqqnDgwAGr\nzL4GXbt2xbRp06BpGjRNg67rKCsrw8iRIwEA+fn50HUdmqZhypQpmDJlitV2xYoVQY2xZs0avPba\na0hOTkZycjK+//57jBgxAi+//HKTrU99YUfPMI3AhAkT8M0332DFihW4+OKLrfLz589bsd5Dhw7h\n17/+NR5//HErJPPtt99C0zTU1NTgn//8J+bPn4/nnnsOAHD55ZcjLS0NeXl5OH/+PJYtW4adO3da\n+6IOGzYMuq7jvffeQ01NDT788EMcOXIE6enpgBHvP3/+PCorK1FdXY3z58+jqqoqoB2LFi1Ceno6\nUlNTverOnz+P8+fPg4hw/vx5/Pjjj377qa6uxrlz51BdXY2qqiqcP38e1dXVPtuuWbMGo0ePxtKl\nS932oQWAmJgYDB8+HNOnT8fZs2dRUFCA/Px8S0b64IMP4s9//jP+85//AMZ+rv/4xz/cfjnWRV1j\nrFmzBjt37sS2bduwbds2dO7cGW+99RYmTZoU9BhNDqtuQiMcWwm2ZNUNbyVYfwJtJVhcXEx9+vSh\nuLg4Sk5OpmnTprmpZP72t79R586dKTY2lvr27eu1b2pRURFlZmZSdHQ09erVi9asWeNWX1BQQFdf\nfTU5HA7q168fbdiwwapbsGABCSEoIiLCetkVMJ5bCRIR/eQnP6F3333Xy8aDBw+69SWEoNTUVKt+\nwoQJNHHiROs4NzfXa2xT8XLo0CFyOBz0/fffExHRwIEDqU2bNuRwOKy1s28lqGka3XnnnRQbG0vd\nunWjDz74wG1un3zyCfXr148URaHOnTvTiBEj3JRK9jn5Ut0EM4ad1NTUFq+64a0EQ8TpdAYV3rBv\nA2imJdYma1adEqVAN9IUh7KVoGffQWHEP+HjQ7tg7WtqeCtBRhZ4K8FWQEOcoH5O90pNXC/HXM++\n6z5Jd7180BKdPMMw9YMdPcMwjOSwow+RVqszDxLZ7WOYCwF29AzDMJLDjj5EZI9hy24fw1wIsKNn\nGIaRHHb0ISJ7DFt2+xjmQiCyuScgO6auvaHnmbJLdbbq0trXVzoJQC0ocPWVkdGguTAM07phRx8i\ndcWwG+KYfZ2nn9NBM8j60lS9+qrja+6B4Bg9w7R+OHTDMGGmpW4l2Lt3b2trwPj4eLRp0wbZ2dk+\nbdi3bx/uvPNOdOjQAe3atcNtt92GvXv3WvX/93//h169eiEhIQHt27fHXXfdZWWI9IXntoSRkZF4\n7LHH/K7fE088gS5duiApKQkPP/ywW16cI0eOYOjQoUhKSkLnzp3xyCOPuG01GA4CbSW4e/du9OvX\nz8qOmZGRgQLjr+aWSp2OXghxiRBijRBilxBihxDiUaNcEUJ8KoTYI4T4RAiRUFdfMiJ7DFt2+xqD\nlrqV4M6dO62tAUtKStC1a1eMGDHCpw3FxcXIzs7G3r17ceLECfTr18/tl0J6ejrWr1+PM2fOoKio\nCNHR0Xjqqaf8rol9W8Ljx48jJibG79i///3vsWXLFnz99dfYu3cvNm/ejJkzZ1r1jz76KJKSknD8\n+HEUFhZi3bp1+NOf/hTUtbHT0K0EO3fujL/97W9W1suRI0fi7rvvrvf4TUkwT/RVAJ4koqsA3ABg\nkhCiF4ApAFYT0RUA1gCoeyNKhrkAiImJwfTp06185XfccQdSU1OxefNm/P3vf8cvf/lL9OrVC5GR\nkfjtb3+L9evX47vvvgMArFy5Ev/93/+Ntm3bolu3bhg/fry1L+vevXuxdetW5Obmom3bthg+fDj6\n9OmDpUuXAgA+/fRTHD58GC+99BLi4uJw0UUX4ZprrvE5x3Xr1uH06dMYPny4z/p+/fph3LhxSExM\nxEUXXYQnnngCe/bsgW6kyrjkkkvQoUMHwMiKedFFFyE5OTmo9fnwww/RoUMHK6umJytXrsQjjzyC\nhIQEJCUl4dFHH3Xbm3bnzp0YOXIk2rRpgw4dOiArK8tt96eVK1eib9++1tP2jh07gpqXydmzZ7Fs\n2TLMnDkT0dHRSE9PR3Z2Nt577z0AQEJCgpXNs7q6GhERES1+39g6HT0RHSeiQuN9GYDdAC4BkA3A\nfNRYCODORp9tC0T2GLbs9jUFLXErwUWLFuGuu+5CdHS0W1++thKE8YshOTkZiqJYZRs2bEBiYiIS\nEhLw/fffY/bs2VZdoK0E67uRSU1NDQ4fPmzt8JSVlYXFixejoqICR44cwT//+U/cdtttQBNuJQhj\n74CYmBifLJsVAAAgAElEQVS8/PLLWLJkSdD9Nwf1+jBWCJECIA3AFwA6EtEJGL8MhBAdGm2WDFNP\nRJhCTv4yiQaLr60ER40ahQkTJqBHjx5+txJ89913cfz48aC2EjRj4+ZWgu+88w4WLFiADz/8ENnZ\n2Thw4IC1QxMAVFRU4MMPP8TKlSvd+tq2bZtPGw4fPoyHH34Yr776qlt5eno6iouLcezYMYwdOxZP\nP/005s6dCxhbCfqiqKgI69evd3tC9yQrKwtz585FZmYmqqqqrN25zp49C4fDgdzcXAwaNAjx8fGo\nqanB2LFjrd2h7FsJAsCYMWPwwgsv4Isvvgh6K8S6thI00XUdFRUVyM3NxS9/+Uts2bIlqP6bg6Ad\nvRAiDsCHAB4jojIhhGdeTb85Y3NycpCSkgIASExMRFpamvWkaMaAW+vxnDlz/NqjFhQAsTOBnc8B\nxr4NTqfTVXbN1bXH37kkkEpkJPSvvoJj8RAoPRTo1y6AmDMHcbbdiRA7E47XX0fpI4/4nA++c197\np9MJFBYCaWmB7bG3D8K+4cNdx8uWNc/610WoDjoc1LWVYGlpKR5//HGvrQQfeeQRXH755WjXrh1G\njRplfRBYn60EAWDkyJF44YUXsGHDBgwZMsQ6Z+nSpUhKSgrK8Z06dQqDBw/Gww8/7DemnpycjN/9\n7ne47bbbLEfvj/feew8ZGRno1q2b3zbTpk3DmTNnkJaWhqioKDz44IMoLCxEx44dAQCDBw/GiBEj\n8J///AelpaUYN24cpkyZglmzZqGoqAiLFi2y1puIUFlZ6baVYEFBAYQQqKiogBACc+bMAQBkZGRg\nxYoVda6znejoaMyaNQvz5s3D9u3b0adPnzrX1B9OpxMLFiwAAMtfho0gNw+JBLDKcPJm2W7jqR4A\nOgHYzRuPuGNuJGJtLGJsGuK5wYi56YhVZ2vnq6y+m5PUdY61Q2o97AtwSpPQGu6rcePG0aBBg+j8\n+fN+2+zdu5fi4uKouLjYZ/2zzz5Lo0aNstpGR0e7baJx44030ptvvklERO+88w716NHD7fw+ffrQ\nihUr3Mp+8Ytf0IwZM+qcv67r1LdvX3r22WfrbPv5559Tly5d6mzXs2dPWrBgQZ3t7Lz55pvUv39/\nIiI6deoUCSGopKTEqv/444/p6quvJiKihx56iF588cWg+vW38Uh5eTm1bduW9u/fb5WNGTOGpk6d\n6rOfyspKio6Opn379tXLLpOm2HgkWHnlXwB8TUT2X9crAOQY78cCWB7eX0GtA9lj2LLb11i0xK0E\nYYRh1q5di7Fjxwacf2lpKW699VZkZGTghRde8KpfvHgxvv/+e8AIxzz33HPWPPyxceNGHD16tE6F\nytGjRy3J6RdffIGZM2da6ph27dqhc+fO+POf/4zq6moUFxdj4cKF1ofOjbmVoPm5wurVq1FYWIia\nmhqUlJTgySefxBVXXIHLLrss6DGanCCe5tMBVAMoBLAVwBYAWQBUAKsB7AHwKYDEC/GJPhAyPNGH\n8ZSw0pLvq5a6lSAR0e9//3saMGCAz3nbtxJcuHAhRUREWPM3bTC3+5s2bRpdcsklFBcXR6mpqTRl\nyhSqqKiw+vLcSpCMp+2xY8d6jeu5leD69espJSWFYmNjqVevXvT++++7td+0aRNlZGRQYmIitW/f\nnkaOHEknT5606ht7K8ElS5ZQr169yOFwUHJyMt1zzz106NAhn/0EA28l2AoItNWe9YHguoFWGc0g\nq9yMI4s8AQxYC8rMdNWtG2i181WGem43WNc5EMa3bX1cJ3/2BTilSeCtBBlZ4K0EGYZhmJBhRx8i\nssewZbePYS4EOKlZE6BEub5kElSCs8oSoL/759qiqhRIX+FW1pCMlCJPuOZi9KVlZLiyZC53jRee\n7ckZhmlp8BN9iASj69Yma1a64TrZmA20cf+yBm0YCop01/DqVVX1zkpJMwj6Od3tXP2cDj0+HrrH\nF0RMONcNw7R+2NEzDMNIDjv6EJE9hi27fQxzIcAxeqZV0q1bNwgRFuUZwzQrgdJBhAt+og8R2WPY\nLdW+gwcPhvwlkrVr14blyygt9cX2tY7XwYMHG/3nhR09wzCM5LCjDxHZY9gy2yezbWD7GBvs6JsI\ntaDAWx+fJ1w69nDTf7mlsw84p9mqpfFvdlTV9WIYJuywow+RYGPYelWVtz7e0LWHnTbxQWns9XN6\nnfr+JovR67rr1YS01M8fwgXbx5iwo2cYhpEcdvQhInucUGb7ZLYNbB9jgx09wzCM5LCjDxHZ44Qy\n2yezbWD7GBvs6BmGYSSHHX2I1CtOWFkSlOyxPjSGPNOudKyPfX4VkkaFVd9QKWWYJZiyx3jZPsaE\nc900JRuzoQ9YG9YuG0Oe2VCVo9/zjIra6nAPwDBMIPiJPkRkjxPKbJ/MtoHtY2ywo2cYhpEcdvQh\nInucUGb7ZLYNbB9jgx09wzCM5LCjDxHZ44Qy2yezbWD7GBvs6BmGYSSHHX2I1DdOqERGApUlQbcX\necJKJawWFEBUlQZ3XlUphNMJ4XRCKQl+PE887VNVQAhACZTdOFS9exOlLJY9xsv2MSaso29itIwM\niDwBBKmnpxkEABBOJ/SqKtAtQyCC+JOVNgytPV8IiLXh0e/rOkCuKcHvlq2h6t1ZL88wYYWf6ENE\n9jihzPbJbBvYPsZGnY5eCPGOEOKEEGK7rWyGEOKwEGKL8cpq9JkyDMMwDSKYJ/p3AQz2Uf5HIrrW\neK1qhLm1CmSPE8psn8y2ge1jbNTp6ImowE9yEn8RWoZhGKYFEUqM/mEhRKEQ4m0hREIY59SqkD1O\nKLN9MtsGto+x0VDVzZ8APE9EJISYCeCPAMb7a5yTk4OUlBQAQGJiItLS0qw/u8yL1VqPCwsLrWN1\ntorKA5XIH5Vf+2dlYSHst6PT6QS+AzDAKIid6b5Y3wFIKAQ8xnM7v7AQSEuz2judTtf4BQWu/goL\noUQp0M/pcPzagXwASkmJS61j1CMtDYnXLofjxReBmBirv9rxvO0zZgCn07YeDoervrS0tt4627Ze\ntRYErPd7bI7nKmwx15+P+Thcx06nEwsWLAAAy1+GC0GmVi5QIyG6Acgnoj71qTPqKZgxZEDkuaJZ\ndkkkAFBmpuv9uoGgGWTJK61yA8rM9Koj44bw6svE6NPexq0vAJRrFlJAaaY1lrCae9so3OWVhNrG\n1nnw6MCoECDf9Z6D+ptAoIkxjGQIIUBEYQmRBxu6EfaYvBCik61uOICd4ZgMwzAME36CkVcuBrAR\nQE8hxCEhxDgALwkhtgshCo0gxBNNM92Wh+xxQpntk9k2sH2MjTpj9EQ0ykfxu40zHYZhGCbc8Ddj\nQ0R2La/M9slsG9g+xgY7eoZhGMlhRx8isscJZbZPZtvA9jE2OHtluOm/3JUaODLw0iqRka40wlWl\noEiHV71aUFBnHwCgRCkuGWX/5S45lo++Ap5vpDDW4+P9tjEzBmuae+FpH21OQ4UGBZctXwDdWAcN\nABQFp3WjkaKENUOlz/kxDGPBjj5EvOKEbeK9te4+0DIyrPe+2upVVZauPWA/kzWrD1o3MOj0x9b5\n2dmu8/2kMc7MzPTtk3UdqvshAECF7tLLxzvd10HToJo6eI0C5DiuPw39nSF7jJftY0w4dMMwDCM5\n7OhDRPY4ocz2yWwb2D7GBjt6hmEYyWFHHyKyxwlltk9m28D2MTbY0TMMw0gOO/oQ8RUnFHkCoqrU\nkkeaUkklSoE6W/XRS90okZFeckulpMSVmtg8jlIa1LeJWlDgUsksL8BpqIAQyBcuuaai1MoYTTQo\nrgqD01ChC8VeFBAhal/NgewxXraPMWFH3wjQDAJtGGpJKPWqKmgZGdAma9DPNUwLqGVkuEkyYUgj\n9aqq2uPJoQnJLUlnfJUlk3SgzNW35i1jTILmJl5XoUOp0YLWsxPVvhiGaTzY0YeI7HFCma2T/tqx\nfYwBO3qGYRjJYUcfIrLHCWW2Tvprx/YxBuzoGYZhJIcdfYjIHieU2Trprx3bxxiwo2cYhpEcdvQh\n4nj9dTctu0l9NPO+NPLBYqY7ViIjvYTuSpQCca4OUbuiQIuyHRuidnX5cmB5ARwzZyJxeb7Vt6mv\n9xTLB9TOe1R66vGhqnV04NGX0YGq+uirHsge42X7GBNOUxwiZdXVgE3LbqJN1lx54oPAUx9fH9zO\n9RC6a5M1w28HmIemIck+TyLA6TTy01ehLDYWiHdYfavQa4Xvwq0b/8N4COu90grrRp/BfHNK06x2\nYUxpzzBSw0/0oZKW1twzaFwktk/2GC/bx5iwo2cYhpEcdvShUljY3DNoXCS2T/YYL9vHmLCjZxiG\nkRx29KEicQwbkNs+2WO8bB9jwo4+BOzySXW26pUm2FfaYCVKCZhOWIlS3FIc+2xTWorE0pLAc/Mh\nPRTCldoYlYHP9RqvpMQlt6xrwGAlkvCd91gXCjShWu/1iBC0kwzDWLCjDwH9nG7FsPVzuleaYF9p\ng7XJWsB0wtpkDTW3DAkoudSGDoU+ZGjguene8kMiV2pjbMwOeK4bhYWudMjx8YHb6bqXjDIgPvIe\nKzWaS74JQCUNCjWuflL2GC/bx5iwo2cYhpGcOh29EOIdIcQJIcR2W5kihPhUCLFHCPGJECKh0Wfa\nUpE4hg3IbZ/sMV62jzEJ5on+XQCDPcqmAFhNRFcAWANgaiPNj2EYhgmROh09ERUA8AyWZgNYaLxf\nCODOxpleK0BinTkgt32yx3jZPsakoTH6DkR0Aq5fBMcBdAjvtBiGYZhwEa6kZgG3d87JyUFKSgoA\nIDExEWlpaVZ8zfyt3FKOHb92AABK3yr1qldnq9B368DVL0D52U1A/+XA9p1AYWGtZLKwEE57/NB8\nIg73fI219ayP27ABiI9HEYYDAJKwzDUfsz0AfAdggOvY4XACd9UeixwBpL3qis1XliCuvNxtlymn\na1AM3zQcmAw4nct8zsfVmQBefdVtPZzWv37sgTtOhwMoK0OmId30bF97RmZtewCZpaW+j23rlZmZ\nGXi9VRXOykogP7/F3J/1Oa7TvlZ+LJt9TqcTCxYsAADLX4YLQRTQR7saCdENQD4R9TGOdwPIJKIT\nQohOANYS0U/8nEvBjNFSMDNO0gzvOYs8AZpBEDZnRpmZEE4nyLhw9vfmsdkufJMUtdke7WtrZn+0\nZYIUIBDcy0UugAFrXWUDM4FcYR2b9gAA1g0E5Rpdr10LDMx09UXkvk72cW12e66N25zIjx0B5u32\nnsgr2aXVp30uPuZWr3Vu6LkMEyJCCBBRcClw6yDY0I3wSEK7AkCO8X4sgDq+TSMvsscJZbZO+mvH\n9jEGwcgrFwPYCKCnEOKQEGIcgFkAfiGE2ANgkHHMMAzDtEDqjNET0Sg/VbeEfzqtj8zMTEDiJwuZ\nlcqy67DZPsaEvxnLMAwjOezoQ0T2OKHM1kl/7dg+xoAdPcMwjOSwow8SzzTESmQkUFmCzMxMKJGR\nUAsKoBYUBEwvHHgAH3mFg68OWBcIz8zCnpLFTEWxGp2GCnX5cpf8sn9wQitRGomI/ILg56co0IUC\nRTHSFkNxm5O/fhTFe+5WWx8pkQFg+PDMBq9ba0D2GLbs9oUTdvRB4pmGWMvIsNL9ahkZ0KuqoFdV\nBUwvHHgAH3mFg68OWBcIz8zCXpJxTbMaqdChx8e7dPFt6khbbFAzJAPkqAp+fpoGpUaDphlpi0lz\nm5O/fjTNe+5WWx8pkRHEmjKMLLCjDxHZ44Ry2yezbbJfO/ntCyfs6BmGYSSHHX2IyB4nlNs+mW2T\n/drJb184YUfPMAwjOezoQ0T2OKHc9slsm+zXTn77wkm40hRf8PiTVTZYbhkkpqRQUYztYTykhKeh\nQoNLHqkKAS1KQIlKdN9JpkIBKkuABibK06BAVYCI1fkAgJpbhni1sebnWSiEt8bTF2Ybo496yyIV\nBXqECpU0aEKFkggsRyUy0AZA4E3NVbV++54zTEuDHX2ImHFCf7LKBsstg8RNUigMKaFNUK5ChzC2\nCyACkgRAk+GWahmzNdAs4zQPiWIwcdAkaCANoABPWJrmkf8UPrSdgTDbGn3UWxapaVCEMFIZ64AO\nDA3y1NYqwZQ9hi27feGEQzcMwzCSw44+RGSPE8psn7yWuZD52uECsC+csKNnGIaRHHb0ISJ7nFBm\n++S1zIXM1w4XgH3hhB09wzCM5LCjDxHZ44Qy2yevZS5kvna4AOwLJyyvrAO1oADF53QoUQrUggLo\nVVVAZYlX2uKwjmloxH2pD1XVpfpWVeCAUKC4Celr6wFXml/AXTJpyewXRELxl/LX1P5HKRiyeAhK\nM0uhREZCrF2LxJJSt7YRU1VgMurUogOAumIFiAAlp86mAdGg4DRUJMGliTdLTXShABGqW5l9vXSh\nQKFazaQmVJecv4bF8oycsKOvA72qCoh0QJusQTidoMxMiDwBHQDN8MzpG6Yxg0hHrOuAQjbHZDh8\n+7mW47Lp1y2ZfXYGND/Tt7T/GRpEnnAv84CigheZ6w6H65wQ/WkSNJBhlN1hm6hUW2+NbVsvBa71\nMiO8KnSv7w/IgOwxbNntCyccumEYhpEcdvQhIn2c8LvmnkDjIfmVk/7elN2+cMKOnmEYRnLY0YeI\n9HHC1OaeQOMh+ZWT/t6U3b5wwo6eYRhGctjRG6izVYg84ZJM9l+OiNX5UGfXnQs37HFCVcVpqOZb\nCAHoEbVlJsFk9rW3rau9qvppEyBGL9JXAAPWAv2XQy0ocKU6DgK1oAARq/Mh8gREnghqnU3stmjC\n96SDWRtdKMhHnNuxKbMUwj0NskfmZ6iq65p4NlS9i5oV2WPYstsXTlheaaCf0y25pHA6QUZZ009E\nt1y6rqM2ra4HDcnwW8ew7imPg4AiHS65qdPpkqFuzHY5/rrGMiSr1nrnBZ8H37JFGNJKjdzSMsNf\nSmQPlBoNDqcTGDjQOoaZxhjuXXpkfnZJNWEsmK1C9y5imBYBP9GHiPRxQplj9JJfO7aPMQnpiV4I\ncRDAGQA1ACqJ6KfhmxrDMAwTDkJ9oq8BkElEfS9UJy99nFBmHb3k147tY0xCdfSCwz8MwzAtm1Cd\nNAH4lxDiSyHEg2GaU6tC+jghx+hbLWwfYxKq6iadiI4JIdobDn83ERWEaW4MwzBMGAjJ0RPRMeP/\nU0KIjwD8FICXo8/JyUFKSgoAIDExEWlpadZvYzPO1pzHQ3bsAPovr60vLATS0lxlX21AXNt4qEbq\nXqfT6YpbG0+6c+bMqbc9Q4YAbdpkQtOM+iFDkGmIv52VlUBZmaWbdzqXAQD6WKl5l8Fp+1an1b9x\nHBfnhLPMR71t/Lg41/j+ju3to76MguPXDpS+VepWD3v/5noBQGEhxIKBQM5a93rz6auw0PV/Wpor\n7fNXXwGxM937C+b6+bHf6SrENUKBBhVOch2bLez92W2xznc4jPVX3PrzPB+GBHMtADgcyGzTBoBm\nq6///aiqQGWlE/n54bm/3exrQT9v4TqWzT6n04kFCxYAgOUvw4Wg+oqnzROFiAEQQURlQohYAJ8C\nyCOiTz3aUUPHaCqEsdhkLL6w3UBYNxA0g9zamLpvmkFwOp31/hPS1FlbyyKEu4jdLsQ2hNkCBELt\n/16id3unXgM0HKfTiYHGGrgNZ18Pc73WDQQGrLXK7Ovp1s6GZ9ug8bTXLoD3FMMTeS2xaVumoaP3\nWjePa+A1hHkNbO0EyGs6DTUpHDTk3mxNyG6fEAJEFJZvZYTyRN8RwEdCCDL6+V9PJ38hIPONBtO+\ndc09i8bhgrh2EiO7feGkwY6eiL4DkBbe6TAMwzDhhqWRISK7lldm+2S2DWwfY4MdPcMwjOSwow8R\n2eOEMtsns21g+xgbnL0ySJTIBi6VmbO2Pukm/aALBUTCdx7e+uQtridKlOKVYVKkr0BilI8UwZGR\nrpTFdZShsgRKtP98vmbqYm2y5vPYzV4zj7DnGisK9GKgR6C0wX7WTYNiqWAUVQWgAaoKvRjQBIDE\nIPI+w3XdrbdQAV2HLhT0SKydq33avkzxvIU80yDX59aq63YM4+3KtCDY0QeJlpHhs7xOiZcevlTH\nSk2An75G+sl0Op21zjUItIwMLxmlrzJszIY2w7+O0DNFtFfKaLu9nnmEbeUKAN2PQM3pdCLTz7ol\nQQPITBNtdKDrsNy7VocG0nbda9+68hgrQvi9LXyZ4tk22FvK171Z17lhvF0bHdnlleGEQzcMwzCS\nw44+RGR/opDZPpltA9vH2GBHzzAMIzns6ENEdi2vzPbJbBvYPsYGO3qGYRjJueBUN2pBAfSqKoiq\nUhARhBBApMOt3pRS6hW+FRlKlIIDecXAXLVWtWHXvOl6rWzP/N9DM6coLnWFogD2UVTV/djePmjC\nKLUMJg6qREa6rZXnsTpbBX62zKorPqf7lGbCdn3QfzmULTlQZ6u1qp/+yxGxOh+Jm8d6K4HsC+o5\nPz/L4WVbCOumG9kyoQKnPepOQ3Upaaj2+tqHEgLQhGqcp9VKLI0MprX9uaSaihKcOqYlxrDDKd9s\nifa1VC44R69XVbllUvQUyZn1ALy04ybaZA2YIoAK20+b/SfPVwZJD82ceaPXJaXzbB8UTSyC1jIy\n3NbKOv6Fy3a7NNKs8yettF8fbbLmfg3axIN8SS0R2OaglyOEdVNqDF2kDnjK9lXo1o1mXl+voYRt\njSyJpW6cX4uueyfXbE20JvmmTHDoJkRkjxPKbJ/MtoHtY2ywo2cYhpEcdvQhInucUGb7ZLYNbB9j\ngx09wzCM5LCjDxHZ44Qy2yezbWD7GBvs6BmGYSRHenmlOluFfk6HEqWg+LqFsFRpRppcvaoKACCq\nSiGcTrd0xIo/rbdN6545fLiVelYhm3bMU5Nt6Lw1KK40t4ac7rSRD9HUFyuKpaprETQ0DmqmNvZc\nQzPdsCcRq/Ndb2zfaYCRDlk4nRDpK9yksOpsFcXXLbT0+OZ1hE2rDwA1twxxO8e8F7TJmpdtZipl\nLSPD6/JpQoUOBUqip6GK13tNNyT9pBsbifvBn6hcVV33iO1EDQpUIYz7xZbiWLjaXqZobmmQVdKg\nKC77zK91mMOYXzmwT9+UPQar0Q8HAcz3We6J/fqFos+/EFIzS+/o9XM6aAa59NiRDmDdQOAWstLk\nmnp6u0Mw8Zee1+0HwRA2qwIg2MTNnneNphk/9MZPvqGvV42fOLNPV1rcUK1ufsy18/wugk8NPADy\ncPD2cvv3Htz6iXTUOvh1A63rrA9Y6/ULA573gg/svyy8UtuTXvudCDdDNa/3SYaTJojA19OfR7Wn\nQzZIggai2vvFbV5WV7VpkD2/yiH83JqeX/eAj+92NBYBzA9bX419bmuBQzchInuUUOY4qMy2uZDb\nPvmvX/hgR88wDCM57OhDRHYlr8xaZZltcyG3ffJfv/DBjp5hGEZy2NGHiOxRQpnjoDLb5kJu++S/\nfuGjaVQ3ntvah9qdIYPDhqFucjl7efF1CwGbPA/9lwOG7M8u8RNVpd5zBaBC85KlqaotA7GvT+oV\nBZqRvdAu2VJ9KwoBI70tAHfZnmcO2zCmHW5qLHmlKVv1IVm1XwMiQxUzYK1bymiLyhJgwFrrvRKt\nBk57bKQ5Vmer1lyse6KwEKKqFImbxwLpK1zplW3Km1ojAq+/XZbpdZ5x/dSFrvtRGQvoxQoUu55W\nCJ/yXGXBQggQtOxsiCkqhNCgCde5miHf1KMFlKhEaMWqJcfUoEBVVcTFLcOhchUQOgiAumI5lOUA\n1BzrplYU4ECxWps9U1GgKBo0oUKFDg2KS7qpeesQ9QgVCrnaJBnplU3MFMu+fu7NNMwwUjhD9aNt\n9BjP/vNnb26mboaHYsiznWfXZl+nrfygmvuQ/gZshTSNow+zfsn6YfSQy9nL7fI6uzzPM/Wtl6zS\nmKvuQ5am2xV2RrlblFDTkGSoJ3U/GYw9UWpcN5DbbdSCbqpQ46BeEtUMb9u8rsEvXLJXe8poi43Z\nICPFsZnuOGDaY0POaUorYUo5TTmm0+lq45G+2t2IwNfD5y8Hj/N0My22BtfVNm8s0+ECXhrIYttc\nKEo37j3XuQoAkQsABJqhQbGdmwQNpAuUUqbLgRO5uo43+tPdUyKbbWDMQaPaMlWI2uYeN7IpOVU9\n5JxGY7/rpdnMV6H7b+oxnmeKZvPe1P0oXwPJRO19qbYJ+JJOt9qc0DY4dMMwDCM5ITl6IUSWEOIb\nIcReIcTk8E2r9SB7lFDmOKjMtoHtY2w02NELISIAvAFgMICrANwrhOgV3um1fAqbewKNTGGhvBbK\nbBvYPsZGKE/0PwWwj4iKiKgSwAcAssM4t1ZBcXNPoJEpLpbXQpltA9vH2AjF0XcB8L3t+LBRxjAM\nw7Qg+MPYEDnY3BNoZA4elNdCmW0D28fYEORLlxTMiUL8HEAuEWUZx1NcMmia7dGuYQMwDMNc4BBR\nWLSdoTj6iwDsATAIwDEA/wFwLxHtDsfEGIZhmPDQ4C9MEVG1EOJhAJ8aIaB32MkzDMO0PBr8RM8w\nDMO0Dhrtw1hZvkwlhDgohNgmhNgqhPiPUaYIIT4VQuwRQnwihEiwtZ8qhNgnhNgthLi1WSfvAyHE\nO0KIE0KI7bayetsjhLhWCLHduL5zmsMWX/ixb4YQ4rAQYovxyrLVtRr7hBCXCCHWCCF2CSF2CCEe\nhUTXz4d9j0Cu69dWCLHJ8CW7hBAvoqmuHxGF/WX8AtkPoBuANsb3ino1xliN/QLwLQDFo2w2gP82\n3k8GMMt4fyWArUZILMVYA9HcNnjMPQNAGoDtodgDYBOAfsb7fwAY3Ny2BbBvBoAnfbT9SWuyD0An\nAGnG+zjjM7Jesly/APZJcf2MucQY/18E4AsA6U1x/RrriV6mL1MJH3/5ZANYaLxfCOBO4/1QAB8Q\nURURHQSwz1iLFgMRFfjIOFUve4QQnQA4iOhLo90i2znNih/74Gfn1uzWZB8RHSeiQuN9GYDdAC6R\n5UpYMqsAAAJiSURBVPr5sc/8bk6rv35w2XXWeNvW8Ct6U1y/xnL0Mn2ZigD8SwjxpRDiAaOsIxGd\ngHFzAuhglHvafaSV2N2hnvZ0Ma6pSWu4vg8LIQqFEG/b/jRutfYJIVKMv1y+aMD92Jrs22QUSXH9\nhBARQoitAI4DcBLR101x/fgLU3WTTkTXArgdwCQhxI2wMn9byPaJtmz2/AlAdyJKM37AXmnuCYWC\nECIOwIcAHjOefKW6H33YJ831I6IaIupr/CV2oxAisymuX2M5+iMALrUdX2KUtTqI6Jjx/ykAHxuh\nmBNCiI5w3ZSdAJw0mh8B0NV2emuxu772tCo7iegUkSUvm28Lp7U6+4QQkYYTfI+IlhvF0lw/X/bJ\ndP1MiKjEiK1f3xTXr7Ec/ZcALhNCdBNCXAzgHgArGmmsRkMIEWM8XUAIEQvgVgA7DFtyjGZjAZg/\ncCsA3COEuFgIkQrgMuOLZC0N4RHzrJc9xp+XZ4QQPxVCCAD32c5pCbjZZ/zwmAwHsNN43xrt+wuA\nr4lorq1MpuvnZZ8s108I0c4MOwkhogH8wviwtfGvXyN+upxlfGq+D8CU5v60u4E2pBqKoa2Gg59i\nlKsAVhv2fQog0XbOVOPT8d0Abm1uG3zYtBjAUQDnARwCMM7Y4Khe9gC4zliTfQDmNrddddi3CMB2\n41p+bMREW519hkKj2nZPbjF+zup9P7Yy+2S5flcbNm0FsA3A09RAf1Jf+/gLUwzDMJLDH8YyDMNI\nDjt6hmEYyWFHzzAMIzns6BmGYSSHHT3DMIzksKNnGIaRHHb0DMMwksOOnmEYRnL+HwbiW8IrP3Uk\nAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fc525c65e10>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "print len(nts)\n",
    "#(y, x) = zip(*p.map(plotNT,nts))\n",
    "(y, x) = zip(*map(plotNT,nts))\n",
    "print len(x),len(y[0])\n",
    "for i in range(0,len(y)) :\n",
    "    plt.step(x[0],y[i],where='post',label=runids[i])\n",
    "plt.grid(True)\n",
    "plt.legend(loc='upper right')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}