2017-10-13-lpc-testdrive/nested-structures-evaluated.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "# Calculations on nested structures"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [],
   "source": [
    "import numpy\n",
    "MHT_pt = numpy.empty(55432, dtype=numpy.float32)\n",
    "MHT_phi = numpy.empty(55432, dtype=numpy.float32)\n",
    "\n",
    "def calculateMHT(events, MHT_pt, MHT_phi):\n",
    "    i = 0\n",
    "    for event in events:\n",
    "        mht_px = 0.0\n",
    "        mht_py = 0.0\n",
    "        mht_pz = 0.0\n",
    "        for electron in event.Electron:\n",
    "            mht_px += electron.pt * cosh(electron.eta) * sin(electron.phi)\n",
    "            mht_py += electron.pt * cosh(electron.eta) * cos(electron.phi)\n",
    "            mht_pz += electron.pt * sinh(electron.eta)\n",
    "        for muon in event.Muon:\n",
    "            mht_px += muon.pt * cosh(muon.eta) * sin(muon.phi)\n",
    "            mht_py += muon.pt * cosh(muon.eta) * cos(muon.phi)\n",
    "            mht_pz += muon.pt * sinh(muon.eta)\n",
    "        for jet in event.Jet:\n",
    "            if jet.pt > 150:\n",
    "                mht_px += jet.pt * cosh(jet.eta) * sin(jet.phi)\n",
    "                mht_py += jet.pt * cosh(jet.eta) * cos(jet.phi)\n",
    "                mht_pz += jet.pt * sinh(jet.eta)\n",
    "\n",
    "        MHT_pt[i] = sqrt(mht_px**2 + mht_py**2)\n",
    "        MHT_phi[i] = -atan2(MHT_pt[i], mht_pz)\n",
    "        i += 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Time for one run: 1.84681391716 seconds --> 30.01 kHz\n"
     ]
    }
   ],
   "source": [
    "import time\n",
    "import math\n",
    "import uproot\n",
    "tree = uproot.open(\"~/storage/data/nano-TTLHE-2017-09-04-uncompressed.root\")[\"Events\"]\n",
    "\n",
    "startTime = time.time()\n",
    "tree.arrowed.run(calculateMHT, (MHT_pt, MHT_phi), env=math.__dict__)\n",
    "endTime = time.time()\n",
    "\n",
    "print(\"Time for one run: {} seconds --> {:.2f} kHz\".format(endTime - startTime, len(tree)/(endTime - startTime)*1e-3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>MHT_phi</th>\n",
       "      <th>MHT_pt</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-1.016380</td>\n",
       "      <td>124.105011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1.786605</td>\n",
       "      <td>5.250953</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-1.255930</td>\n",
       "      <td>14.964581</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-0.636966</td>\n",
       "      <td>213.986099</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-2.210669</td>\n",
       "      <td>41.393135</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-0.397885</td>\n",
       "      <td>161.648438</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-0.786898</td>\n",
       "      <td>1370.723755</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-2.393102</td>\n",
       "      <td>2427.410400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-1.258100</td>\n",
       "      <td>371.518036</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-1.906085</td>\n",
       "      <td>412.545441</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>-2.242043</td>\n",
       "      <td>87.483330</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>-2.340775</td>\n",
       "      <td>478.170929</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>-1.934738</td>\n",
       "      <td>379.618561</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>-0.980049</td>\n",
       "      <td>68.327057</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>-0.928512</td>\n",
       "      <td>12.108970</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>-0.228699</td>\n",
       "      <td>57.812366</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>-2.305995</td>\n",
       "      <td>81.671326</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>-0.866310</td>\n",
       "      <td>36.737701</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>-0.812940</td>\n",
       "      <td>776.566467</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>-0.781964</td>\n",
       "      <td>431.855164</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>-0.994510</td>\n",
       "      <td>159.442017</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>-0.873884</td>\n",
       "      <td>502.132050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>-1.294286</td>\n",
       "      <td>71.878319</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>-2.231563</td>\n",
       "      <td>83.048508</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>-2.214998</td>\n",
       "      <td>35.620182</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>-2.213963</td>\n",
       "      <td>753.088745</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>-0.872133</td>\n",
       "      <td>29.828629</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>-0.939470</td>\n",
       "      <td>29.672432</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55402</th>\n",
       "      <td>-2.288311</td>\n",
       "      <td>77.511719</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55403</th>\n",
       "      <td>-0.894256</td>\n",
       "      <td>391.617401</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55404</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55405</th>\n",
       "      <td>-2.208956</td>\n",
       "      <td>105.504913</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55406</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55407</th>\n",
       "      <td>-1.587248</td>\n",
       "      <td>90.397926</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55408</th>\n",
       "      <td>-0.449359</td>\n",
       "      <td>36.007977</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55409</th>\n",
       "      <td>-2.206260</td>\n",
       "      <td>120.284187</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55410</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55411</th>\n",
       "      <td>-1.594901</td>\n",
       "      <td>117.565231</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55412</th>\n",
       "      <td>-2.494349</td>\n",
       "      <td>127.289001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55413</th>\n",
       "      <td>-2.317117</td>\n",
       "      <td>23.957169</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55414</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55415</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55416</th>\n",
       "      <td>-1.218979</td>\n",
       "      <td>167.833115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55417</th>\n",
       "      <td>-0.790888</td>\n",
       "      <td>1033.098755</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55418</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55419</th>\n",
       "      <td>-2.708817</td>\n",
       "      <td>10.382696</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55420</th>\n",
       "      <td>-0.621232</td>\n",
       "      <td>122.020813</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55421</th>\n",
       "      <td>-2.138939</td>\n",
       "      <td>33.520893</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55422</th>\n",
       "      <td>-2.320051</td>\n",
       "      <td>31.969828</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55423</th>\n",
       "      <td>-0.803855</td>\n",
       "      <td>27.767122</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55424</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55425</th>\n",
       "      <td>-2.343607</td>\n",
       "      <td>227.071152</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55426</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55427</th>\n",
       "      <td>-0.882814</td>\n",
       "      <td>147.067337</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55428</th>\n",
       "      <td>-0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55429</th>\n",
       "      <td>-0.950985</td>\n",
       "      <td>14.167120</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55430</th>\n",
       "      <td>-0.423686</td>\n",
       "      <td>5.102069</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>55431</th>\n",
       "      <td>-2.159389</td>\n",
       "      <td>594.546265</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>55432 rows × 2 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        MHT_phi       MHT_pt\n",
       "0     -1.016380   124.105011\n",
       "1     -1.786605     5.250953\n",
       "2     -1.255930    14.964581\n",
       "3     -0.636966   213.986099\n",
       "4     -2.210669    41.393135\n",
       "5     -0.397885   161.648438\n",
       "6     -0.786898  1370.723755\n",
       "7     -2.393102  2427.410400\n",
       "8     -1.258100   371.518036\n",
       "9     -1.906085   412.545441\n",
       "10    -2.242043    87.483330\n",
       "11    -2.340775   478.170929\n",
       "12    -0.000000     0.000000\n",
       "13    -1.934738   379.618561\n",
       "14    -0.980049    68.327057\n",
       "15    -0.928512    12.108970\n",
       "16    -0.228699    57.812366\n",
       "17    -0.000000     0.000000\n",
       "18    -2.305995    81.671326\n",
       "19    -0.866310    36.737701\n",
       "20    -0.812940   776.566467\n",
       "21    -0.781964   431.855164\n",
       "22    -0.994510   159.442017\n",
       "23    -0.873884   502.132050\n",
       "24    -1.294286    71.878319\n",
       "25    -2.231563    83.048508\n",
       "26    -2.214998    35.620182\n",
       "27    -2.213963   753.088745\n",
       "28    -0.872133    29.828629\n",
       "29    -0.939470    29.672432\n",
       "...         ...          ...\n",
       "55402 -2.288311    77.511719\n",
       "55403 -0.894256   391.617401\n",
       "55404 -0.000000     0.000000\n",
       "55405 -2.208956   105.504913\n",
       "55406 -0.000000     0.000000\n",
       "55407 -1.587248    90.397926\n",
       "55408 -0.449359    36.007977\n",
       "55409 -2.206260   120.284187\n",
       "55410 -0.000000     0.000000\n",
       "55411 -1.594901   117.565231\n",
       "55412 -2.494349   127.289001\n",
       "55413 -2.317117    23.957169\n",
       "55414 -0.000000     0.000000\n",
       "55415 -0.000000     0.000000\n",
       "55416 -1.218979   167.833115\n",
       "55417 -0.790888  1033.098755\n",
       "55418 -0.000000     0.000000\n",
       "55419 -2.708817    10.382696\n",
       "55420 -0.621232   122.020813\n",
       "55421 -2.138939    33.520893\n",
       "55422 -2.320051    31.969828\n",
       "55423 -0.803855    27.767122\n",
       "55424 -0.000000     0.000000\n",
       "55425 -2.343607   227.071152\n",
       "55426 -0.000000     0.000000\n",
       "55427 -0.882814   147.067337\n",
       "55428 -0.000000     0.000000\n",
       "55429 -0.950985    14.167120\n",
       "55430 -0.423686     5.102069\n",
       "55431 -2.159389   594.546265\n",
       "\n",
       "[55432 rows x 2 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas\n",
    "pandas.DataFrame({\"MHT_pt\": MHT_pt, \"MHT_phi\": MHT_phi})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "\n",
       "require(['notebook'],\n",
       "  function() {\n",
       "    IPython.CodeCell.config_defaults.highlight_modes['magic_text/x-c++src'] = {'reg':[/^%%cpp/]};\n",
       "    console.log(\"JupyROOT - %%cpp magic configured\");\n",
       "  }\n",
       ");\n"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to JupyROOT 6.08/04\n",
      "Will write selected trees to .\n",
      "Processed    10000/   55432 entries (elapsed time     9.9s, curr speed    1.006 kHz, avg speed    1.006 kHz), accepted    10001/   10001 events (100.00%)\n",
      "Processed    20000/   55432 entries (elapsed time    19.5s, curr speed    1.044 kHz, avg speed    1.025 kHz), accepted    20001/   20001 events (100.00%)\n",
      "Processed    30000/   55432 entries (elapsed time    29.0s, curr speed    1.052 kHz, avg speed    1.034 kHz), accepted    30001/   30001 events (100.00%)\n",
      "Processed    40000/   55432 entries (elapsed time    38.5s, curr speed    1.053 kHz, avg speed    1.038 kHz), accepted    40001/   40001 events (100.00%)\n",
      "Processed    50000/   55432 entries (elapsed time    48.4s, curr speed    1.011 kHz, avg speed    1.033 kHz), accepted    50001/   50001 events (100.00%)\n",
      "Processed 55432 entries from /home/pivarski/storage/data/nano-TTLHE-2017-09-04-uncompressed.root, selected 55432 entries\n",
      "Done ./nano-TTLHE-2017-09-04-uncompressed_Friend.root\n",
      "Time for one run: 54.5998768806 seconds\n"
     ]
    }
   ],
   "source": [
    "import time\n",
    "from PhysicsTools.NanoAODTools.postprocessing.framework.postprocessor import PostProcessor\n",
    "from PhysicsTools.NanoAODTools.postprocessing.examples.mhtProducer import *\n",
    "\n",
    "p = PostProcessor(\".\", [\"/home/pivarski/storage/data/nano-TTLHE-2017-09-04-uncompressed.root\"], \"Jet_pt > 150\", \"keep_and_drop.txt\", [mht()], compression=\"ZLIB:0\", friend=True)\n",
    "\n",
    "startTime = time.time()\n",
    "p.run()\n",
    "endTime = time.time()\n",
    "\n",
    "print(\"Time for one run: {0} seconds\".format(endTime - startTime))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Just as `tree` has a `.pandas` connector, it has a `.arrowed` connector to present itself as an Arrow dataset, which is used to compile the function.\n",
    "\n",
    "This includes a ROOT branch ⟷ nested schema mapping."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "List [\n",
      "  countarray = None\n",
      "  Record {\n",
      "    name = u'Events'\n",
      "    run: 'run'\n",
      "    luminosityBlock: 'luminosityBlock'\n",
      "    event: 'event'\n",
      "    nElectron: 'nElectron'\n",
      "    Electron: List [\n",
      "      countarray = 'nElectron'\n",
      "      Record {\n",
      "        name = u'Electron'\n",
      "        deltaEtaSC: 'Electron_deltaEtaSC'\n",
      "        dxy: 'Electron_dxy'\n",
      "        eta: 'Electron_eta'\n",
      "        mass: 'Electron_mass'\n",
      "        mvaSpring16GP: 'Electron_mvaSpring16GP'\n",
      "        mvaSpring16HZZ: 'Electron_mvaSpring16HZZ'\n",
      "        phi: 'Electron_phi'\n",
      "        pt: 'Electron_pt'\n",
      "        pterr: 'Electron_pterr'\n",
      "        r9: 'Electron_r9'\n",
      "        sieie: 'Electron_sieie'\n",
      "        sip3d: 'Electron_sip3d'\n",
      "        charge: 'Electron_charge'\n",
      "        cutBased: 'Electron_cutBased'\n",
      "        cutBased_HLTPreSel: 'Electron_cutBased_HLTPreSel'\n",
      "        jet: 'Electron_jet'\n",
      "        mvaSpring16GP_WP80: 'Electron_mvaSpring16GP_WP80'\n",
      "        mvaSpring16GP_WP90: 'Electron_mvaSpring16GP_WP90'\n",
      "        mvaSpring16HZZ_WPL: 'Electron_mvaSpring16HZZ_WPL'\n",
      "        pdgId: 'Electron_pdgId'\n",
      "        cleanmask: 'Electron_cleanmask'\n",
      "        }\n",
      "      ]\n",
      "    nFatJet: 'nFatJet'\n",
      "    FatJet: List [\n",
      "      countarray = 'nFatJet'\n",
      "      Record {\n",
      "        name = u'FatJet'\n",
      "        area: 'FatJet_area'\n",
      "        btagCMVA: 'FatJet_btagCMVA'\n",
      "        btagDeepB: 'FatJet_btagDeepB'\n",
      "        btagDeepBB: 'FatJet_btagDeepBB'\n",
      "        eta: 'FatJet_eta'\n",
      "        mass: 'FatJet_mass'\n",
      "        mpruned: 'FatJet_mpruned'\n",
      "        msoftdrop: 'FatJet_msoftdrop'\n",
      "        phi: 'FatJet_phi'\n",
      "        pt: 'FatJet_pt'\n",
      "        tau1: 'FatJet_tau1'\n",
      "        tau2: 'FatJet_tau2'\n",
      "        tau3: 'FatJet_tau3'\n",
      "        subJet1: 'FatJet_subJet1'\n",
      "        subJet2: 'FatJet_subJet2'\n",
      "        subJet3: 'FatJet_subJet3'\n",
      "        }\n",
      "      ]\n",
      "    nGenJet: 'nGenJet'\n",
      "    GenJet: List [\n",
      "      countarray = 'nGenJet'\n",
      "      Record {\n",
      "        name = u'GenJet'\n",
      "        eta: 'GenJet_eta'\n",
      "        mass: 'GenJet_mass'\n",
      "        phi: 'GenJet_phi'\n",
      "        pt: 'GenJet_pt'\n",
      "        }\n",
      "      ]\n",
      "    nGenPart: 'nGenPart'\n",
      "    GenPart: List [\n",
      "      countarray = 'nGenPart'\n",
      "      Record {\n",
      "        name = u'GenPart'\n",
      "        eta: 'GenPart_eta'\n",
      "        phi: 'GenPart_phi'\n",
      "        pt: 'GenPart_pt'\n",
      "        mother: 'GenPart_mother'\n",
      "        pdgId: 'GenPart_pdgId'\n",
      "        status: 'GenPart_status'\n",
      "        }\n",
      "      ]\n",
      "    genWeight: 'genWeight'\n",
      "    LHEWeight: Record {\n",
      "      name = u'LHEWeight'\n",
      "      originalXWGTUP: 'LHEWeight_originalXWGTUP'\n",
      "      }\n",
      "    nLHEPdfWeight: 'nLHEPdfWeight'\n",
      "    LHEPdfWeight: List [\n",
      "      countarray = 'nLHEPdfWeight'\n",
      "      'LHEPdfWeight'\n",
      "      ]\n",
      "    nLHEScaleWeight: 'nLHEScaleWeight'\n",
      "    LHEScaleWeight: List [\n",
      "      countarray = 'nLHEScaleWeight'\n",
      "      'LHEScaleWeight'\n",
      "      ]\n",
      "    nJet: 'nJet'\n",
      "    Jet: List [\n",
      "      countarray = 'nJet'\n",
      "      Record {\n",
      "        name = u'Jet'\n",
      "        area: 'Jet_area'\n",
      "        btagCMVA: 'Jet_btagCMVA'\n",
      "        btagDeepB: 'Jet_btagDeepB'\n",
      "        btagDeepBB: 'Jet_btagDeepBB'\n",
      "        btagDeepC: 'Jet_btagDeepC'\n",
      "        eta: 'Jet_eta'\n",
      "        mass: 'Jet_mass'\n",
      "        phi: 'Jet_phi'\n",
      "        pt: 'Jet_pt'\n",
      "        qgl: 'Jet_qgl'\n",
      "        rawFactor: 'Jet_rawFactor'\n",
      "        electron1: 'Jet_electron1'\n",
      "        electron2: 'Jet_electron2'\n",
      "        muon1: 'Jet_muon1'\n",
      "        muon2: 'Jet_muon2'\n",
      "        nConstituents: 'Jet_nConstituents'\n",
      "        nElectrons: 'Jet_nElectrons'\n",
      "        nMuons: 'Jet_nMuons'\n",
      "        puId: 'Jet_puId'\n",
      "        genJetIdx: 'Jet_genJetIdx'\n",
      "        hadroFlavour: 'Jet_hadroFlavour'\n",
      "        partonFlavour: 'Jet_partonFlavour'\n",
      "        cleanmask: 'Jet_cleanmask'\n",
      "        }\n",
      "      ]\n",
      "    LHE: Record {\n",
      "      name = u'LHE'\n",
      "      HT: 'LHE_HT'\n",
      "      HTIncoming: 'LHE_HTIncoming'\n",
      "      Vpt: 'LHE_Vpt'\n",
      "      Njets: 'LHE_Njets'\n",
      "      Nb: 'LHE_Nb'\n",
      "      Nc: 'LHE_Nc'\n",
      "      Nuds: 'LHE_Nuds'\n",
      "      Nglu: 'LHE_Nglu'\n",
      "      }\n",
      "    MET: Record {\n",
      "      name = u'MET'\n",
      "      phi: 'MET_phi'\n",
      "      pt: 'MET_pt'\n",
      "      rawEt: 'MET_rawEt'\n",
      "      significance: 'MET_significance'\n",
      "      sumEt: 'MET_sumEt'\n",
      "      }\n",
      "    nMuon: 'nMuon'\n",
      "    Muon: List [\n",
      "      countarray = 'nMuon'\n",
      "      Record {\n",
      "        name = u'Muon'\n",
      "        dxy: 'Muon_dxy'\n",
      "        eta: 'Muon_eta'\n",
      "        mass: 'Muon_mass'\n",
      "        phi: 'Muon_phi'\n",
      "        pt: 'Muon_pt'\n",
      "        ptErr: 'Muon_ptErr'\n",
      "        segmentComp: 'Muon_segmentComp'\n",
      "        charge: 'Muon_charge'\n",
      "        jet: 'Muon_jet'\n",
      "        nStations: 'Muon_nStations'\n",
      "        pdgId: 'Muon_pdgId'\n",
      "        mediumId: 'Muon_mediumId'\n",
      "        cleanmask: 'Muon_cleanmask'\n",
      "        }\n",
      "      ]\n",
      "    nPhoton: 'nPhoton'\n",
      "    Photon: List [\n",
      "      countarray = 'nPhoton'\n",
      "      Record {\n",
      "        name = u'Photon'\n",
      "        eta: 'Photon_eta'\n",
      "        mass: 'Photon_mass'\n",
      "        phi: 'Photon_phi'\n",
      "        pt: 'Photon_pt'\n",
      "        charge: 'Photon_charge'\n",
      "        jet: 'Photon_jet'\n",
      "        pdgId: 'Photon_pdgId'\n",
      "        cleanmask: 'Photon_cleanmask'\n",
      "        }\n",
      "      ]\n",
      "    METpuppi: Record {\n",
      "      name = u'METpuppi'\n",
      "      phi: 'METpuppi_phi'\n",
      "      pt: 'METpuppi_pt'\n",
      "      sumEt: 'METpuppi_sumEt'\n",
      "      }\n",
      "    fixedGridRhoFastjetAll: 'fixedGridRhoFastjetAll'\n",
      "    fixedGridRhoFastjetCentralCalo: 'fixedGridRhoFastjetCentralCalo'\n",
      "    fixedGridRhoFastjetCentralNeutral: 'fixedGridRhoFastjetCentralNeutral'\n",
      "    nSubJet: 'nSubJet'\n",
      "    SubJet: List [\n",
      "      countarray = 'nSubJet'\n",
      "      Record {\n",
      "        name = u'SubJet'\n",
      "        btagCMVA: 'SubJet_btagCMVA'\n",
      "        btagDeepB: 'SubJet_btagDeepB'\n",
      "        btagDeepBB: 'SubJet_btagDeepBB'\n",
      "        eta: 'SubJet_eta'\n",
      "        mass: 'SubJet_mass'\n",
      "        phi: 'SubJet_phi'\n",
      "        pt: 'SubJet_pt'\n",
      "        }\n",
      "      ]\n",
      "    nTau: 'nTau'\n",
      "    Tau: List [\n",
      "      countarray = 'nTau'\n",
      "      Record {\n",
      "        name = u'Tau'\n",
      "        dxy: 'Tau_dxy'\n",
      "        dz: 'Tau_dz'\n",
      "        eta: 'Tau_eta'\n",
      "        mass: 'Tau_mass'\n",
      "        phi: 'Tau_phi'\n",
      "        pt: 'Tau_pt'\n",
      "        rawMVArun2: 'Tau_rawMVArun2'\n",
      "        decayMode: 'Tau_decayMode'\n",
      "        idDecayMode: 'Tau_idDecayMode'\n",
      "        idDecayModeNewDMs: 'Tau_idDecayModeNewDMs'\n",
      "        jet: 'Tau_jet'\n",
      "        cleanmask: 'Tau_cleanmask'\n",
      "        }\n",
      "      ]\n",
      "    PV: Record {\n",
      "      name = u'PV'\n",
      "      ndof: 'PV_ndof'\n",
      "      x: 'PV_x'\n",
      "      y: 'PV_y'\n",
      "      z: 'PV_z'\n",
      "      chi2: 'PV_chi2'\n",
      "      }\n",
      "    nSV: 'nSV'\n",
      "    SV: List [\n",
      "      countarray = 'nSV'\n",
      "      Record {\n",
      "        name = u'SV'\n",
      "        dlen: 'SV_dlen'\n",
      "        dlenSig: 'SV_dlenSig'\n",
      "        chi2: 'SV_chi2'\n",
      "        eta: 'SV_eta'\n",
      "        mass: 'SV_mass'\n",
      "        ndof: 'SV_ndof'\n",
      "        phi: 'SV_phi'\n",
      "        pt: 'SV_pt'\n",
      "        x: 'SV_x'\n",
      "        y: 'SV_y'\n",
      "        z: 'SV_z'\n",
      "        }\n",
      "      ]\n",
      "    HLTriggerFirstPath: 'HLTriggerFirstPath'\n",
      "    HLT: Record {\n",
      "      name = u'HLT'\n",
      "      AK8PFJet360_TrimMass30: 'HLT_AK8PFJet360_TrimMass30'\n",
      "      AK8PFJet380_TrimMass30: 'HLT_AK8PFJet380_TrimMass30'\n",
      "      AK8PFJet400_TrimMass30: 'HLT_AK8PFJet400_TrimMass30'\n",
      "      AK8PFJet420_TrimMass30: 'HLT_AK8PFJet420_TrimMass30'\n",
      "      AK8PFHT750_TrimMass50: 'HLT_AK8PFHT750_TrimMass50'\n",
      "      AK8PFHT800_TrimMass50: 'HLT_AK8PFHT800_TrimMass50'\n",
      "      AK8PFHT850_TrimMass50: 'HLT_AK8PFHT850_TrimMass50'\n",
      "      AK8PFHT900_TrimMass50: 'HLT_AK8PFHT900_TrimMass50'\n",
      "      CaloJet500_NoJetID: 'HLT_CaloJet500_NoJetID'\n",
      "      CaloJet550_NoJetID: 'HLT_CaloJet550_NoJetID'\n",
      "      Trimuon5_3p5_2_Upsilon_Muon: 'HLT_Trimuon5_3p5_2_Upsilon_Muon'\n",
      "      DoubleEle33_CaloIdL_MW: 'HLT_DoubleEle33_CaloIdL_MW'\n",
      "      DoubleEle24_eta2p1_WPTight_Gsf: 'HLT_DoubleEle24_eta2p1_WPTight_Gsf'\n",
      "      DoubleEle8_CaloIdM_TrackIdM_Mass8_DZ_PFHT350: 'HLT_DoubleEle8_CaloIdM_...'\n",
      "      DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT350: 'HLT_DoubleEle8_CaloIdM_Tra...'\n",
      "      DoubleMu4_3_Bs: 'HLT_DoubleMu4_3_Bs'\n",
      "      DoubleMu4_3_Jpsi_Displaced: 'HLT_DoubleMu4_3_Jpsi_Displaced'\n",
      "      DoubleMu4_JpsiTrk_Displaced: 'HLT_DoubleMu4_JpsiTrk_Displaced'\n",
      "      DoubleMu4_LowMassNonResonantTrk_Displaced: 'HLT_DoubleMu4_LowMassNonRe...'\n",
      "      DoubleMu3_Trk_Tau3mu: 'HLT_DoubleMu3_Trk_Tau3mu'\n",
      "      DoubleMu4_PsiPrimeTrk_Displaced: 'HLT_DoubleMu4_PsiPrimeTrk_Displaced'\n",
      "      DoubleMu4_Mass8_DZ_PFHT350: 'HLT_DoubleMu4_Mass8_DZ_PFHT350'\n",
      "      DoubleMu8_Mass8_PFHT350: 'HLT_DoubleMu8_Mass8_PFHT350'\n",
      "      Mu3_PFJet40: 'HLT_Mu3_PFJet40'\n",
      "      Mu7p5_L2Mu2_Jpsi: 'HLT_Mu7p5_L2Mu2_Jpsi'\n",
      "      Mu7p5_L2Mu2_Upsilon: 'HLT_Mu7p5_L2Mu2_Upsilon'\n",
      "      Mu7p5_Track2_Jpsi: 'HLT_Mu7p5_Track2_Jpsi'\n",
      "      Mu7p5_Track3p5_Jpsi: 'HLT_Mu7p5_Track3p5_Jpsi'\n",
      "      Mu7p5_Track7_Jpsi: 'HLT_Mu7p5_Track7_Jpsi'\n",
      "      Mu7p5_Track2_Upsilon: 'HLT_Mu7p5_Track2_Upsilon'\n",
      "      Mu7p5_Track3p5_Upsilon: 'HLT_Mu7p5_Track3p5_Upsilon'\n",
      "      Mu7p5_Track7_Upsilon: 'HLT_Mu7p5_Track7_Upsilon'\n",
      "      DoublePhoton33_CaloIdL: 'HLT_DoublePhoton33_CaloIdL'\n",
      "      DoublePhoton70: 'HLT_DoublePhoton70'\n",
      "      DoublePhoton85: 'HLT_DoublePhoton85'\n",
      "      Ele20_WPTight_Gsf: 'HLT_Ele20_WPTight_Gsf'\n",
      "      Ele20_WPLoose_Gsf: 'HLT_Ele20_WPLoose_Gsf'\n",
      "      Ele20_eta2p1_WPLoose_Gsf: 'HLT_Ele20_eta2p1_WPLoose_Gsf'\n",
      "      DiEle27_WPTightCaloOnly_L1DoubleEG: 'HLT_DiEle27_WPTightCaloOnly_L1Dou...'\n",
      "      Ele27_WPTight_Gsf: 'HLT_Ele27_WPTight_Gsf'\n",
      "      Ele35_WPTight_Gsf: 'HLT_Ele35_WPTight_Gsf'\n",
      "      Ele35_WPTight_Gsf_L1EGMT: 'HLT_Ele35_WPTight_Gsf_L1EGMT'\n",
      "      Ele38_WPTight_Gsf: 'HLT_Ele38_WPTight_Gsf'\n",
      "      Ele40_WPTight_Gsf: 'HLT_Ele40_WPTight_Gsf'\n",
      "      Ele32_WPTight_Gsf_L1DoubleEG: 'HLT_Ele32_WPTight_Gsf_L1DoubleEG'\n",
      "      HT450_Beamspot: 'HLT_HT450_Beamspot'\n",
      "      HT300_Beamspot: 'HLT_HT300_Beamspot'\n",
      "      IsoMu20_eta2p1_LooseChargedIsoPFTau27_eta2p1_CrossL1: 'HLT_IsoMu20_eta...'\n",
      "      IsoMu20_eta2p1_MediumChargedIsoPFTau27_eta2p1_CrossL1: 'HLT_IsoMu20_et...'\n",
      "      IsoMu20_eta2p1_TightChargedIsoPFTau27_eta2p1_CrossL1: 'HLT_IsoMu20_eta...'\n",
      "      IsoMu20_eta2p1_LooseChargedIsoPFTau27_eta2p1_TightID_CrossL1: 'HLT_Iso...'\n",
      "      IsoMu20_eta2p1_MediumChargedIsoPFTau27_eta2p1_TightID_CrossL1: 'HLT_Is...'\n",
      "      IsoMu20_eta2p1_TightChargedIsoPFTau27_eta2p1_TightID_CrossL1: 'HLT_Iso...'\n",
      "      IsoMu24_eta2p1_LooseChargedIsoPFTau20_SingleL1: 'HLT_IsoMu24_eta2p1_Lo...'\n",
      "      IsoMu24_eta2p1_MediumChargedIsoPFTau20_SingleL1: 'HLT_IsoMu24_eta2p1_M...'\n",
      "      IsoMu24_eta2p1_TightChargedIsoPFTau20_SingleL1: 'HLT_IsoMu24_eta2p1_Ti...'\n",
      "      IsoMu24_eta2p1_LooseChargedIsoPFTau20_TightID_SingleL1: 'HLT_IsoMu24_e...'\n",
      "      IsoMu24_eta2p1_MediumChargedIsoPFTau20_TightID_SingleL1: 'HLT_IsoMu24_...'\n",
      "      IsoMu24_eta2p1_TightChargedIsoPFTau20_TightID_SingleL1: 'HLT_IsoMu24_e...'\n",
      "      IsoMu20: 'HLT_IsoMu20'\n",
      "      IsoMu24: 'HLT_IsoMu24'\n",
      "      IsoMu24_eta2p1: 'HLT_IsoMu24_eta2p1'\n",
      "      IsoMu27: 'HLT_IsoMu27'\n",
      "      UncorrectedJetE30_NoBPTX: 'HLT_UncorrectedJetE30_NoBPTX'\n",
      "      UncorrectedJetE30_NoBPTX3BX: 'HLT_UncorrectedJetE30_NoBPTX3BX'\n",
      "      UncorrectedJetE60_NoBPTX3BX: 'HLT_UncorrectedJetE60_NoBPTX3BX'\n",
      "      UncorrectedJetE70_NoBPTX3BX: 'HLT_UncorrectedJetE70_NoBPTX3BX'\n",
      "      L1SingleMu18: 'HLT_L1SingleMu18'\n",
      "      L1SingleMu25: 'HLT_L1SingleMu25'\n",
      "      L2Mu10: 'HLT_L2Mu10'\n",
      "      L2Mu10_NoVertex_NoBPTX3BX: 'HLT_L2Mu10_NoVertex_NoBPTX3BX'\n",
      "      L2Mu10_NoVertex_NoBPTX: 'HLT_L2Mu10_NoVertex_NoBPTX'\n",
      "      L2Mu45_NoVertex_3Sta_NoBPTX3BX: 'HLT_L2Mu45_NoVertex_3Sta_NoBPTX3BX'\n",
      "      L2Mu40_NoVertex_3Sta_NoBPTX3BX: 'HLT_L2Mu40_NoVertex_3Sta_NoBPTX3BX'\n",
      "      L2Mu50: 'HLT_L2Mu50'\n",
      "      DoubleL2Mu50: 'HLT_DoubleL2Mu50'\n",
      "      Mu17_TrkIsoVVL_Mu8_TrkIsoVVL: 'HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL'\n",
      "      Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ: 'HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ'\n",
      "      Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ_Mass8: 'HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoV...'\n",
      "      Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ_Mass3p8: 'HLT_Mu17_TrkIsoVVL_Mu8_TrkIs...'\n",
      "      Mu25_TkMu0_Onia: 'HLT_Mu25_TkMu0_Onia'\n",
      "      Mu30_TkMu0_Onia: 'HLT_Mu30_TkMu0_Onia'\n",
      "      Mu20_TkMu0_Phi: 'HLT_Mu20_TkMu0_Phi'\n",
      "      Mu25_TkMu0_Phi: 'HLT_Mu25_TkMu0_Phi'\n",
      "      Mu20: 'HLT_Mu20'\n",
      "      Mu27: 'HLT_Mu27'\n",
      "      Mu50: 'HLT_Mu50'\n",
      "      OldMu100: 'HLT_OldMu100'\n",
      "      TkMu100: 'HLT_TkMu100'\n",
      "      DiPFJet15_NoCaloMatched: 'HLT_DiPFJet15_NoCaloMatched'\n",
      "      DiPFJet25_NoCaloMatched: 'HLT_DiPFJet25_NoCaloMatched'\n",
      "      DiPFJet15_FBEta3_NoCaloMatched: 'HLT_DiPFJet15_FBEta3_NoCaloMatched'\n",
      "      DiPFJet25_FBEta3_NoCaloMatched: 'HLT_DiPFJet25_FBEta3_NoCaloMatched'\n",
      "      DiPFJetAve40: 'HLT_DiPFJetAve40'\n",
      "      DiPFJetAve60: 'HLT_DiPFJetAve60'\n",
      "      DiPFJetAve80: 'HLT_DiPFJetAve80'\n",
      "      DiPFJetAve140: 'HLT_DiPFJetAve140'\n",
      "      DiPFJetAve200: 'HLT_DiPFJetAve200'\n",
      "      DiPFJetAve260: 'HLT_DiPFJetAve260'\n",
      "      DiPFJetAve320: 'HLT_DiPFJetAve320'\n",
      "      DiPFJetAve400: 'HLT_DiPFJetAve400'\n",
      "      DiPFJetAve500: 'HLT_DiPFJetAve500'\n",
      "      DiPFJetAve15_HFJEC: 'HLT_DiPFJetAve15_HFJEC'\n",
      "      DiPFJetAve25_HFJEC: 'HLT_DiPFJetAve25_HFJEC'\n",
      "      DiPFJetAve35_HFJEC: 'HLT_DiPFJetAve35_HFJEC'\n",
      "      DiPFJetAve60_HFJEC: 'HLT_DiPFJetAve60_HFJEC'\n",
      "      DiPFJetAve80_HFJEC: 'HLT_DiPFJetAve80_HFJEC'\n",
      "      DiPFJetAve100_HFJEC: 'HLT_DiPFJetAve100_HFJEC'\n",
      "      DiPFJetAve160_HFJEC: 'HLT_DiPFJetAve160_HFJEC'\n",
      "      DiPFJetAve220_HFJEC: 'HLT_DiPFJetAve220_HFJEC'\n",
      "      DiPFJetAve300_HFJEC: 'HLT_DiPFJetAve300_HFJEC'\n",
      "      AK8PFJet40: 'HLT_AK8PFJet40'\n",
      "      AK8PFJet60: 'HLT_AK8PFJet60'\n",
      "      AK8PFJet80: 'HLT_AK8PFJet80'\n",
      "      AK8PFJet140: 'HLT_AK8PFJet140'\n",
      "      AK8PFJet200: 'HLT_AK8PFJet200'\n",
      "      AK8PFJet260: 'HLT_AK8PFJet260'\n",
      "      AK8PFJet320: 'HLT_AK8PFJet320'\n",
      "      AK8PFJet400: 'HLT_AK8PFJet400'\n",
      "      AK8PFJet450: 'HLT_AK8PFJet450'\n",
      "      AK8PFJet500: 'HLT_AK8PFJet500'\n",
      "      AK8PFJet550: 'HLT_AK8PFJet550'\n",
      "      PFJet40: 'HLT_PFJet40'\n",
      "      PFJet60: 'HLT_PFJet60'\n",
      "      PFJet80: 'HLT_PFJet80'\n",
      "      PFJet140: 'HLT_PFJet140'\n",
      "      PFJet200: 'HLT_PFJet200'\n",
      "      PFJet260: 'HLT_PFJet260'\n",
      "      PFJet320: 'HLT_PFJet320'\n",
      "      PFJet400: 'HLT_PFJet400'\n",
      "      PFJet450: 'HLT_PFJet450'\n",
      "      PFJet500: 'HLT_PFJet500'\n",
      "      PFJet550: 'HLT_PFJet550'\n",
      "      PFJetFwd40: 'HLT_PFJetFwd40'\n",
      "      PFJetFwd60: 'HLT_PFJetFwd60'\n",
      "      PFJetFwd80: 'HLT_PFJetFwd80'\n",
      "      PFJetFwd140: 'HLT_PFJetFwd140'\n",
      "      PFJetFwd200: 'HLT_PFJetFwd200'\n",
      "      PFJetFwd260: 'HLT_PFJetFwd260'\n",
      "      PFJetFwd320: 'HLT_PFJetFwd320'\n",
      "      PFJetFwd400: 'HLT_PFJetFwd400'\n",
      "      PFJetFwd450: 'HLT_PFJetFwd450'\n",
      "      PFJetFwd500: 'HLT_PFJetFwd500'\n",
      "      AK8PFJetFwd40: 'HLT_AK8PFJetFwd40'\n",
      "      AK8PFJetFwd60: 'HLT_AK8PFJetFwd60'\n",
      "      AK8PFJetFwd80: 'HLT_AK8PFJetFwd80'\n",
      "      AK8PFJetFwd140: 'HLT_AK8PFJetFwd140'\n",
      "      AK8PFJetFwd200: 'HLT_AK8PFJetFwd200'\n",
      "      AK8PFJetFwd260: 'HLT_AK8PFJetFwd260'\n",
      "      AK8PFJetFwd320: 'HLT_AK8PFJetFwd320'\n",
      "      AK8PFJetFwd400: 'HLT_AK8PFJetFwd400'\n",
      "      AK8PFJetFwd450: 'HLT_AK8PFJetFwd450'\n",
      "      AK8PFJetFwd500: 'HLT_AK8PFJetFwd500'\n",
      "      PFHT180: 'HLT_PFHT180'\n",
      "      PFHT250: 'HLT_PFHT250'\n",
      "      PFHT370: 'HLT_PFHT370'\n",
      "      PFHT430: 'HLT_PFHT430'\n",
      "      PFHT510: 'HLT_PFHT510'\n",
      "      PFHT590: 'HLT_PFHT590'\n",
      "      PFHT680: 'HLT_PFHT680'\n",
      "      PFHT780: 'HLT_PFHT780'\n",
      "      PFHT890: 'HLT_PFHT890'\n",
      "      PFHT1050: 'HLT_PFHT1050'\n",
      "      PFHT500_PFMET100_PFMHT100_IDTight: 'HLT_PFHT500_PFMET100_PFMHT100_IDTight'\n",
      "      PFHT500_PFMET110_PFMHT110_IDTight: 'HLT_PFHT500_PFMET110_PFMHT110_IDTight'\n",
      "      PFHT700_PFMET85_PFMHT85_IDTight: 'HLT_PFHT700_PFMET85_PFMHT85_IDTight'\n",
      "      PFHT700_PFMET95_PFMHT95_IDTight: 'HLT_PFHT700_PFMET95_PFMHT95_IDTight'\n",
      "      PFHT800_PFMET75_PFMHT75_IDTight: 'HLT_PFHT800_PFMET75_PFMHT75_IDTight'\n",
      "      PFHT800_PFMET85_PFMHT85_IDTight: 'HLT_PFHT800_PFMET85_PFMHT85_IDTight'\n",
      "      PFMET110_PFMHT110_IDTight: 'HLT_PFMET110_PFMHT110_IDTight'\n",
      "      PFMET120_PFMHT120_IDTight: 'HLT_PFMET120_PFMHT120_IDTight'\n",
      "      PFMET130_PFMHT130_IDTight: 'HLT_PFMET130_PFMHT130_IDTight'\n",
      "      PFMET140_PFMHT140_IDTight: 'HLT_PFMET140_PFMHT140_IDTight'\n",
      "      PFMET100_PFMHT100_IDTight_CaloBTagCSV_3p1: 'HLT_PFMET100_PFMHT100_IDTi...'\n",
      "      PFMET110_PFMHT110_IDTight_CaloBTagCSV_3p1: 'HLT_PFMET110_PFMHT110_IDTi...'\n",
      "      PFMET120_PFMHT120_IDTight_CaloBTagCSV_3p1: 'HLT_PFMET120_PFMHT120_IDTi...'\n",
      "      PFMET130_PFMHT130_IDTight_CaloBTagCSV_3p1: 'HLT_PFMET130_PFMHT130_IDTi...'\n",
      "      PFMET140_PFMHT140_IDTight_CaloBTagCSV_3p1: 'HLT_PFMET140_PFMHT140_IDTi...'\n",
      "      PFMET120_PFMHT120_IDTight_PFHT60: 'HLT_PFMET120_PFMHT120_IDTight_PFHT60'\n",
      "      PFMETNoMu120_PFMHTNoMu120_IDTight_PFHT60: 'HLT_PFMETNoMu120_PFMHTNoMu1...'\n",
      "      PFMETTypeOne120_PFMHT120_IDTight_PFHT60: 'HLT_PFMETTypeOne120_PFMHT120...'\n",
      "      PFMETTypeOne110_PFMHT110_IDTight: 'HLT_PFMETTypeOne110_PFMHT110_IDTight'\n",
      "      PFMETTypeOne120_PFMHT120_IDTight: 'HLT_PFMETTypeOne120_PFMHT120_IDTight'\n",
      "      PFMETTypeOne130_PFMHT130_IDTight: 'HLT_PFMETTypeOne130_PFMHT130_IDTight'\n",
      "      PFMETTypeOne140_PFMHT140_IDTight: 'HLT_PFMETTypeOne140_PFMHT140_IDTight'\n",
      "      PFMETNoMu110_PFMHTNoMu110_IDTight: 'HLT_PFMETNoMu110_PFMHTNoMu110_IDTight'\n",
      "      PFMETNoMu120_PFMHTNoMu120_IDTight: 'HLT_PFMETNoMu120_PFMHTNoMu120_IDTight'\n",
      "      PFMETNoMu130_PFMHTNoMu130_IDTight: 'HLT_PFMETNoMu130_PFMHTNoMu130_IDTight'\n",
      "      PFMETNoMu140_PFMHTNoMu140_IDTight: 'HLT_PFMETNoMu140_PFMHTNoMu140_IDTight'\n",
      "      PFMET120_PFMHT120_IDTight_HFCleaned: 'HLT_PFMET120_PFMHT120_IDTight_HF...'\n",
      "      PFMETNoMu120_PFMHTNoMu120_IDTight_HFCleaned: 'HLT_PFMETNoMu120_PFMHTNo...'\n",
      "      PFMETTypeOne120_PFMHT120_IDTight_HFCleaned: 'HLT_PFMETTypeOne120_PFMHT...'\n",
      "      PFMET120_PFMHT120_IDTight_PFHT60_HFCleaned: 'HLT_PFMET120_PFMHT120_IDT...'\n",
      "      MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight: 'HLT_MonoCentral...'\n",
      "      MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight: 'HLT_MonoCentral...'\n",
      "      MonoCentralPFJet80_PFMETNoMu130_PFMHTNoMu130_IDTight: 'HLT_MonoCentral...'\n",
      "      MonoCentralPFJet80_PFMETNoMu140_PFMHTNoMu140_IDTight: 'HLT_MonoCentral...'\n",
      "      L1ETMHadSeeds: 'HLT_L1ETMHadSeeds'\n",
      "      CaloMHT90: 'HLT_CaloMHT90'\n",
      "      CaloMET80_NotCleaned: 'HLT_CaloMET80_NotCleaned'\n",
      "      CaloMET90_NotCleaned: 'HLT_CaloMET90_NotCleaned'\n",
      "      CaloMET100_NotCleaned: 'HLT_CaloMET100_NotCleaned'\n",
      "      CaloMET110_NotCleaned: 'HLT_CaloMET110_NotCleaned'\n",
      "      CaloMET250_NotCleaned: 'HLT_CaloMET250_NotCleaned'\n",
      "      CaloMET70_HBHECleaned: 'HLT_CaloMET70_HBHECleaned'\n",
      "      CaloMET80_HBHECleaned: 'HLT_CaloMET80_HBHECleaned'\n",
      "      CaloMET90_HBHECleaned: 'HLT_CaloMET90_HBHECleaned'\n",
      "      CaloMET100_HBHECleaned: 'HLT_CaloMET100_HBHECleaned'\n",
      "      CaloMET250_HBHECleaned: 'HLT_CaloMET250_HBHECleaned'\n",
      "      CaloMET300_HBHECleaned: 'HLT_CaloMET300_HBHECleaned'\n",
      "      CaloMET350_HBHECleaned: 'HLT_CaloMET350_HBHECleaned'\n",
      "      PFMET200_NotCleaned: 'HLT_PFMET200_NotCleaned'\n",
      "      PFMET200_HBHECleaned: 'HLT_PFMET200_HBHECleaned'\n",
      "      PFMET250_HBHECleaned: 'HLT_PFMET250_HBHECleaned'\n",
      "      PFMET300_HBHECleaned: 'HLT_PFMET300_HBHECleaned'\n",
      "      PFMET200_HBHE_BeamHaloCleaned: 'HLT_PFMET200_HBHE_BeamHaloCleaned'\n",
      "      PFMETTypeOne200_HBHE_BeamHaloCleaned: 'HLT_PFMETTypeOne200_HBHE_BeamHa...'\n",
      "      MET105_IsoTrk50: 'HLT_MET105_IsoTrk50'\n",
      "      MET120_IsoTrk50: 'HLT_MET120_IsoTrk50'\n",
      "      SingleJet30_Mu12_SinglePFJet40: 'HLT_SingleJet30_Mu12_SinglePFJet40'\n",
      "      Mu12_DoublePFJets40_CaloBTagCSV_p33: 'HLT_Mu12_DoublePFJets40_CaloBTag...'\n",
      "      Mu12_DoublePFJets100_CaloBTagCSV_p33: 'HLT_Mu12_DoublePFJets100_CaloBT...'\n",
      "      Mu12_DoublePFJets200_CaloBTagCSV_p33: 'HLT_Mu12_DoublePFJets200_CaloBT...'\n",
      "      Mu12_DoublePFJets350_CaloBTagCSV_p33: 'HLT_Mu12_DoublePFJets350_CaloBT...'\n",
      "      Mu12_DoublePFJets40MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_Mu12_DoubleP...'\n",
      "      Mu12_DoublePFJets54MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_Mu12_DoubleP...'\n",
      "      Mu12_DoublePFJets62MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_Mu12_DoubleP...'\n",
      "      DoublePFJets40_CaloBTagCSV_p33: 'HLT_DoublePFJets40_CaloBTagCSV_p33'\n",
      "      DoublePFJets100_CaloBTagCSV_p33: 'HLT_DoublePFJets100_CaloBTagCSV_p33'\n",
      "      DoublePFJets200_CaloBTagCSV_p33: 'HLT_DoublePFJets200_CaloBTagCSV_p33'\n",
      "      DoublePFJets350_CaloBTagCSV_p33: 'HLT_DoublePFJets350_CaloBTagCSV_p33'\n",
      "      DoublePFJets100MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_DoublePFJets100M...'\n",
      "      DoublePFJets116MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_DoublePFJets116M...'\n",
      "      DoublePFJets128MaxDeta1p6_DoubleCaloBTagCSV_p33: 'HLT_DoublePFJets128M...'\n",
      "      Photon300_NoHE: 'HLT_Photon300_NoHE'\n",
      "      Mu8_TrkIsoVVL: 'HLT_Mu8_TrkIsoVVL'\n",
      "      Mu8_DiEle12_CaloIdL_TrackIdL_DZ: 'HLT_Mu8_DiEle12_CaloIdL_TrackIdL_DZ'\n",
      "      Mu8_DiEle12_CaloIdL_TrackIdL: 'HLT_Mu8_DiEle12_CaloIdL_TrackIdL'\n",
      "      Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350_DZ: 'HLT_Mu8_Ele8_CaloIdM_Trac...'\n",
      "      Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT350: 'HLT_Mu8_Ele8_CaloIdM_TrackId...'\n",
      "      Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ: 'HLT_Mu8_TrkIsoVVL_Ele2...'\n",
      "      Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL: 'HLT_Mu8_TrkIsoVVL_Ele23_C...'\n",
      "      Mu17_TrkIsoVVL: 'HLT_Mu17_TrkIsoVVL'\n",
      "      BTagMu_AK4DiJet20_Mu5: 'HLT_BTagMu_AK4DiJet20_Mu5'\n",
      "      BTagMu_AK4DiJet40_Mu5: 'HLT_BTagMu_AK4DiJet40_Mu5'\n",
      "      BTagMu_AK4DiJet70_Mu5: 'HLT_BTagMu_AK4DiJet70_Mu5'\n",
      "      BTagMu_AK4DiJet110_Mu5: 'HLT_BTagMu_AK4DiJet110_Mu5'\n",
      "      BTagMu_AK4DiJet170_Mu5: 'HLT_BTagMu_AK4DiJet170_Mu5'\n",
      "      BTagMu_AK4Jet300_Mu5: 'HLT_BTagMu_AK4Jet300_Mu5'\n",
      "      BTagMu_AK8DiJet170_Mu5: 'HLT_BTagMu_AK8DiJet170_Mu5'\n",
      "      BTagMu_AK8Jet300_Mu5: 'HLT_BTagMu_AK8Jet300_Mu5'\n",
      "      Ele23_Ele12_CaloIdL_TrackIdL_IsoVL_DZ: 'HLT_Ele23_Ele12_CaloIdL_TrackI...'\n",
      "      Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_DZ: 'HLT_Mu23_TrkIsoVVL_El...'\n",
      "      Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL: 'HLT_Mu23_TrkIsoVVL_Ele12...'\n",
      "      Mu12_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL: 'HLT_Mu12_TrkIsoVVL_Ele23...'\n",
      "      Mu12_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL_DZ: 'HLT_Mu12_TrkIsoVVL_El...'\n",
      "      Ele23_Ele12_CaloIdL_TrackIdL_IsoVL: 'HLT_Ele23_Ele12_CaloIdL_TrackIdL_...'\n",
      "      Photon33: 'HLT_Photon33'\n",
      "      Photon50: 'HLT_Photon50'\n",
      "      Photon75: 'HLT_Photon75'\n",
      "      Photon90: 'HLT_Photon90'\n",
      "      Photon120: 'HLT_Photon120'\n",
      "      Photon150: 'HLT_Photon150'\n",
      "      Photon175: 'HLT_Photon175'\n",
      "      Photon200: 'HLT_Photon200'\n",
      "      Photon50_R9Id90_HE10_IsoM: 'HLT_Photon50_R9Id90_HE10_IsoM'\n",
      "      Photon75_R9Id90_HE10_IsoM: 'HLT_Photon75_R9Id90_HE10_IsoM'\n",
      "      Photon90_R9Id90_HE10_IsoM: 'HLT_Photon90_R9Id90_HE10_IsoM'\n",
      "      Photon120_R9Id90_HE10_IsoM: 'HLT_Photon120_R9Id90_HE10_IsoM'\n",
      "      Photon165_R9Id90_HE10_IsoM: 'HLT_Photon165_R9Id90_HE10_IsoM'\n",
      "      Photon90_CaloIdL_PFHT700: 'HLT_Photon90_CaloIdL_PFHT700'\n",
      "      Diphoton30_22_R9Id_OR_IsoCaloId_AND_HE_R9Id_Mass90: 'HLT_Diphoton30_22...'\n",
      "      Diphoton30_22_R9Id_OR_IsoCaloId_AND_HE_R9Id_Mass95: 'HLT_Diphoton30_22...'\n",
      "      Diphoton30PV_18PV_R9Id_AND_IsoCaloId_AND_HE_R9Id_PixelVeto_Mass55: 'HL...'\n",
      "      Diphoton30PV_18PV_R9Id_AND_IsoCaloId_AND_HE_R9Id_NoPixelVeto_Mass55: '...'\n",
      "      Diphoton30EB_18EB_R9Id_OR_IsoCaloId_AND_HE_R9Id_NoPixelVeto_Mass55: 'H...'\n",
      "      Diphoton30EB_18EB_R9Id_OR_IsoCaloId_AND_HE_R9Id_PixelVeto_Mass55: 'HLT...'\n",
      "      Dimuon0_Jpsi_L1_NoOS: 'HLT_Dimuon0_Jpsi_L1_NoOS'\n",
      "      Dimuon0_Jpsi_NoVertexing_NoOS: 'HLT_Dimuon0_Jpsi_NoVertexing_NoOS'\n",
      "      Dimuon0_Jpsi: 'HLT_Dimuon0_Jpsi'\n",
      "      Dimuon0_Jpsi_NoVertexing: 'HLT_Dimuon0_Jpsi_NoVertexing'\n",
      "      Dimuon0_Jpsi_L1_4R_0er1p5R: 'HLT_Dimuon0_Jpsi_L1_4R_0er1p5R'\n",
      "      Dimuon0_Jpsi_NoVertexing_L1_4R_0er1p5R: 'HLT_Dimuon0_Jpsi_NoVertexing_...'\n",
      "      Dimuon0_Jpsi3p5_Muon2: 'HLT_Dimuon0_Jpsi3p5_Muon2'\n",
      "      Dimuon0_Upsilon_L1_4p5: 'HLT_Dimuon0_Upsilon_L1_4p5'\n",
      "      Dimuon0_Upsilon_L1_5: 'HLT_Dimuon0_Upsilon_L1_5'\n",
      "      Dimuon0_Upsilon_L1_4p5NoOS: 'HLT_Dimuon0_Upsilon_L1_4p5NoOS'\n",
      "      Dimuon0_Upsilon_L1_4p5er2p0: 'HLT_Dimuon0_Upsilon_L1_4p5er2p0'\n",
      "      Dimuon0_Upsilon_L1_4p5er2p0M: 'HLT_Dimuon0_Upsilon_L1_4p5er2p0M'\n",
      "      Dimuon0_Upsilon_NoVertexing: 'HLT_Dimuon0_Upsilon_NoVertexing'\n",
      "      Dimuon0_Upsilon_L1_5M: 'HLT_Dimuon0_Upsilon_L1_5M'\n",
      "      Dimuon0_LowMass_L1_0er1p5R: 'HLT_Dimuon0_LowMass_L1_0er1p5R'\n",
      "      Dimuon0_LowMass_L1_0er1p5: 'HLT_Dimuon0_LowMass_L1_0er1p5'\n",
      "      Dimuon0_LowMass: 'HLT_Dimuon0_LowMass'\n",
      "      Dimuon0_LowMass_L1_4: 'HLT_Dimuon0_LowMass_L1_4'\n",
      "      Dimuon0_LowMass_L1_4R: 'HLT_Dimuon0_LowMass_L1_4R'\n",
      "      Dimuon0_LowMass_L1_TM530: 'HLT_Dimuon0_LowMass_L1_TM530'\n",
      "      Dimuon0_Upsilon_Muon_L1_TM0: 'HLT_Dimuon0_Upsilon_Muon_L1_TM0'\n",
      "      Dimuon0_Upsilon_Muon_NoL1Mass: 'HLT_Dimuon0_Upsilon_Muon_NoL1Mass'\n",
      "      TripleMu_5_3_3_Mass3p8to60_DZ: 'HLT_TripleMu_5_3_3_Mass3p8to60_DZ'\n",
      "      TripleMu_10_5_5_DZ: 'HLT_TripleMu_10_5_5_DZ'\n",
      "      TripleMu_12_10_5: 'HLT_TripleMu_12_10_5'\n",
      "      DoubleMu3_DZ_PFMET50_PFMHT60: 'HLT_DoubleMu3_DZ_PFMET50_PFMHT60'\n",
      "      DoubleMu3_DZ_PFMET70_PFMHT70: 'HLT_DoubleMu3_DZ_PFMET70_PFMHT70'\n",
      "      DoubleMu3_DZ_PFMET90_PFMHT90: 'HLT_DoubleMu3_DZ_PFMET90_PFMHT90'\n",
      "      DoubleMu3_Trk_Tau3mu_NoL1Mass: 'HLT_DoubleMu3_Trk_Tau3mu_NoL1Mass'\n",
      "      DoubleMu4_Jpsi_Displaced: 'HLT_DoubleMu4_Jpsi_Displaced'\n",
      "      DoubleMu4_Jpsi_NoVertexing: 'HLT_DoubleMu4_Jpsi_NoVertexing'\n",
      "      DoubleMu4_JpsiTrkTrk_Displaced: 'HLT_DoubleMu4_JpsiTrkTrk_Displaced'\n",
      "      DoubleMu43NoFiltersNoVtx: 'HLT_DoubleMu43NoFiltersNoVtx'\n",
      "      DoubleMu48NoFiltersNoVtx: 'HLT_DoubleMu48NoFiltersNoVtx'\n",
      "      Mu43NoFiltersNoVtx_Photon43_CaloIdL: 'HLT_Mu43NoFiltersNoVtx_Photon43_...'\n",
      "      Mu48NoFiltersNoVtx_Photon48_CaloIdL: 'HLT_Mu48NoFiltersNoVtx_Photon48_...'\n",
      "      DoubleMu20_7_Mass0to30_L1_DM4: 'HLT_DoubleMu20_7_Mass0to30_L1_DM4'\n",
      "      DoubleMu20_7_Mass0to30_L1_DM4EG: 'HLT_DoubleMu20_7_Mass0to30_L1_DM4EG'\n",
      "      HT425: 'HLT_HT425'\n",
      "      HT430_DisplacedDijet40_DisplacedTrack: 'HLT_HT430_DisplacedDijet40_Dis...'\n",
      "      HT430_DisplacedDijet60_DisplacedTrack: 'HLT_HT430_DisplacedDijet60_Dis...'\n",
      "      HT430_DisplacedDijet80_DisplacedTrack: 'HLT_HT430_DisplacedDijet80_Dis...'\n",
      "      HT400_DisplacedDijet40_DisplacedTrack: 'HLT_HT400_DisplacedDijet40_Dis...'\n",
      "      HT650_DisplacedDijet60_Inclusive: 'HLT_HT650_DisplacedDijet60_Inclusive'\n",
      "      HT550_DisplacedDijet80_Inclusive: 'HLT_HT550_DisplacedDijet80_Inclusive'\n",
      "      HT550_DisplacedDijet60_Inclusive: 'HLT_HT550_DisplacedDijet60_Inclusive'\n",
      "      HT650_DisplacedDijet80_Inclusive: 'HLT_HT650_DisplacedDijet80_Inclusive'\n",
      "      HT750_DisplacedDijet80_Inclusive: 'HLT_HT750_DisplacedDijet80_Inclusive'\n",
      "      DiJet110_35_Mjj650_PFMET110: 'HLT_DiJet110_35_Mjj650_PFMET110'\n",
      "      DiJet110_35_Mjj650_PFMET120: 'HLT_DiJet110_35_Mjj650_PFMET120'\n",
      "      DiJet110_35_Mjj650_PFMET130: 'HLT_DiJet110_35_Mjj650_PFMET130'\n",
      "      TripleJet110_35_35_Mjj650_PFMET110: 'HLT_TripleJet110_35_35_Mjj650_PFM...'\n",
      "      TripleJet110_35_35_Mjj650_PFMET120: 'HLT_TripleJet110_35_35_Mjj650_PFM...'\n",
      "      TripleJet110_35_35_Mjj650_PFMET130: 'HLT_TripleJet110_35_35_Mjj650_PFM...'\n",
      "      VBF_DoubleLooseChargedIsoPFTau20_Trk1_eta2p1_Reg: 'HLT_VBF_DoubleLoose...'\n",
      "      VBF_DoubleMediumChargedIsoPFTau20_Trk1_eta2p1_Reg: 'HLT_VBF_DoubleMedi...'\n",
      "      VBF_DoubleTightChargedIsoPFTau20_Trk1_eta2p1_Reg: 'HLT_VBF_DoubleTight...'\n",
      "      Ele30_eta2p1_WPTight_Gsf_CentralPFJet35_EleCleaned: 'HLT_Ele30_eta2p1_...'\n",
      "      Ele28_eta2p1_WPTight_Gsf_HT150: 'HLT_Ele28_eta2p1_WPTight_Gsf_HT150'\n",
      "      Ele28_HighEta_SC20_Mass55: 'HLT_Ele28_HighEta_SC20_Mass55'\n",
      "      DoubleMu20_7_Mass0to30_Photon23: 'HLT_DoubleMu20_7_Mass0to30_Photon23'\n",
      "      Ele15_IsoVVVL_PFHT450_CaloBTagCSV_4p5: 'HLT_Ele15_IsoVVVL_PFHT450_Calo...'\n",
      "      Ele15_IsoVVVL_PFHT450_PFMET50: 'HLT_Ele15_IsoVVVL_PFHT450_PFMET50'\n",
      "      Ele15_IsoVVVL_PFHT450: 'HLT_Ele15_IsoVVVL_PFHT450'\n",
      "      Ele50_IsoVVVL_PFHT450: 'HLT_Ele50_IsoVVVL_PFHT450'\n",
      "      Ele15_IsoVVVL_PFHT600: 'HLT_Ele15_IsoVVVL_PFHT600'\n",
      "      Mu8_TrkIsoVVL_DiPFJet40_DEta3p5_MJJ750_HTT300_PFMETNoMu60: 'HLT_Mu8_Tr...'\n",
      "      Mu10_TrkIsoVVL_DiPFJet40_DEta3p5_MJJ750_HTT350_PFMETNoMu60: 'HLT_Mu10_...'\n",
      "      Mu15_IsoVVVL_PFHT450_CaloBTagCSV_4p5: 'HLT_Mu15_IsoVVVL_PFHT450_CaloBT...'\n",
      "      Mu15_IsoVVVL_PFHT450_PFMET50: 'HLT_Mu15_IsoVVVL_PFHT450_PFMET50'\n",
      "      Mu15_IsoVVVL_PFHT450: 'HLT_Mu15_IsoVVVL_PFHT450'\n",
      "      Mu50_IsoVVVL_PFHT450: 'HLT_Mu50_IsoVVVL_PFHT450'\n",
      "      Mu15_IsoVVVL_PFHT600: 'HLT_Mu15_IsoVVVL_PFHT600'\n",
      "      Dimuon10_PsiPrime_Barrel_Seagulls: 'HLT_Dimuon10_PsiPrime_Barrel_Seagulls'\n",
      "      Dimuon20_Jpsi_Barrel_Seagulls: 'HLT_Dimuon20_Jpsi_Barrel_Seagulls'\n",
      "      Dimuon10_Upsilon_Barrel_Seagulls: 'HLT_Dimuon10_Upsilon_Barrel_Seagulls'\n",
      "      Dimuon12_Upsilon_eta1p5: 'HLT_Dimuon12_Upsilon_eta1p5'\n",
      "      Dimuon14_Phi_Barrel_Seagulls: 'HLT_Dimuon14_Phi_Barrel_Seagulls'\n",
      "      Dimuon18_PsiPrime: 'HLT_Dimuon18_PsiPrime'\n",
      "      Dimuon25_Jpsi: 'HLT_Dimuon25_Jpsi'\n",
      "      Dimuon18_PsiPrime_noCorrL1: 'HLT_Dimuon18_PsiPrime_noCorrL1'\n",
      "      Dimuon24_Upsilon_noCorrL1: 'HLT_Dimuon24_Upsilon_noCorrL1'\n",
      "      Dimuon24_Phi_noCorrL1: 'HLT_Dimuon24_Phi_noCorrL1'\n",
      "      Dimuon25_Jpsi_noCorrL1: 'HLT_Dimuon25_Jpsi_noCorrL1'\n",
      "      DiMu9_Ele9_CaloIdL_TrackIdL_DZ: 'HLT_DiMu9_Ele9_CaloIdL_TrackIdL_DZ'\n",
      "      DiMu9_Ele9_CaloIdL_TrackIdL: 'HLT_DiMu9_Ele9_CaloIdL_TrackIdL'\n",
      "      DoubleIsoMu20_eta2p1: 'HLT_DoubleIsoMu20_eta2p1'\n",
      "      DoubleIsoMu24_eta2p1: 'HLT_DoubleIsoMu24_eta2p1'\n",
      "      TrkMu12_DoubleTrkMu5NoFiltersNoVtx: 'HLT_TrkMu12_DoubleTrkMu5NoFilters...'\n",
      "      TrkMu16_DoubleTrkMu6NoFiltersNoVtx: 'HLT_TrkMu16_DoubleTrkMu6NoFilters...'\n",
      "      TrkMu17_DoubleTrkMu8NoFiltersNoVtx: 'HLT_TrkMu17_DoubleTrkMu8NoFilters...'\n",
      "      Mu8: 'HLT_Mu8'\n",
      "      Mu17: 'HLT_Mu17'\n",
      "      Mu17_Photon30_IsoCaloId: 'HLT_Mu17_Photon30_IsoCaloId'\n",
      "      Ele8_CaloIdL_TrackIdL_IsoVL_PFJet30: 'HLT_Ele8_CaloIdL_TrackIdL_IsoVL_...'\n",
      "      Ele12_CaloIdL_TrackIdL_IsoVL_PFJet30: 'HLT_Ele12_CaloIdL_TrackIdL_IsoV...'\n",
      "      Ele23_CaloIdL_TrackIdL_IsoVL_PFJet30: 'HLT_Ele23_CaloIdL_TrackIdL_IsoV...'\n",
      "      Ele8_CaloIdM_TrackIdM_PFJet30: 'HLT_Ele8_CaloIdM_TrackIdM_PFJet30'\n",
      "      Ele17_CaloIdM_TrackIdM_PFJet30: 'HLT_Ele17_CaloIdM_TrackIdM_PFJet30'\n",
      "      Ele23_CaloIdM_TrackIdM_PFJet30: 'HLT_Ele23_CaloIdM_TrackIdM_PFJet30'\n",
      "      Ele50_CaloIdVT_GsfTrkIdT_PFJet165: 'HLT_Ele50_CaloIdVT_GsfTrkIdT_PFJet165'\n",
      "      Ele115_CaloIdVT_GsfTrkIdT: 'HLT_Ele115_CaloIdVT_GsfTrkIdT'\n",
      "      Ele135_CaloIdVT_GsfTrkIdT: 'HLT_Ele135_CaloIdVT_GsfTrkIdT'\n",
      "      Ele145_CaloIdVT_GsfTrkIdT: 'HLT_Ele145_CaloIdVT_GsfTrkIdT'\n",
      "      Ele200_CaloIdVT_GsfTrkIdT: 'HLT_Ele200_CaloIdVT_GsfTrkIdT'\n",
      "      Ele250_CaloIdVT_GsfTrkIdT: 'HLT_Ele250_CaloIdVT_GsfTrkIdT'\n",
      "      Ele300_CaloIdVT_GsfTrkIdT: 'HLT_Ele300_CaloIdVT_GsfTrkIdT'\n",
      "      PFHT300PT30_QuadPFJet_75_60_45_40: 'HLT_PFHT300PT30_QuadPFJet_75_60_45_40'\n",
      "      PFHT300PT30_QuadPFJet_75_60_45_40_TriplePFBTagCSV_3p0: 'HLT_PFHT300PT3...'\n",
      "      PFHT380_SixPFJet32_DoublePFBTagCSV_2p2: 'HLT_PFHT380_SixPFJet32_Double...'\n",
      "      PFHT380_SixPFJet32_DoublePFBTagDeepCSV_2p2: 'HLT_PFHT380_SixPFJet32_Do...'\n",
      "      PFHT380_SixPFJet32: 'HLT_PFHT380_SixPFJet32'\n",
      "      PFHT430_SixPFJet40_PFBTagCSV_1p5: 'HLT_PFHT430_SixPFJet40_PFBTagCSV_1p5'\n",
      "      PFHT430_SixPFJet40: 'HLT_PFHT430_SixPFJet40'\n",
      "      PFHT350: 'HLT_PFHT350'\n",
      "      PFHT350MinPFJet15: 'HLT_PFHT350MinPFJet15'\n",
      "      Photon60_R9Id90_CaloIdL_IsoL: 'HLT_Photon60_R9Id90_CaloIdL_IsoL'\n",
      "      Photon60_R9Id90_CaloIdL_IsoL_DisplacedIdL: 'HLT_Photon60_R9Id90_CaloId...'\n",
      "      Photon60_R9Id90_CaloIdL_IsoL_DisplacedIdL_PFHT350MinPFJet15: 'HLT_Phot...'\n",
      "      FullTrack_Multiplicity85: 'HLT_FullTrack_Multiplicity85'\n",
      "      FullTrack_Multiplicity100: 'HLT_FullTrack_Multiplicity100'\n",
      "      FullTrack_Multiplicity130: 'HLT_FullTrack_Multiplicity130'\n",
      "      FullTrack_Multiplicity155: 'HLT_FullTrack_Multiplicity155'\n",
      "      ECALHT800: 'HLT_ECALHT800'\n",
      "      DiSC30_18_EIso_AND_HE_Mass70: 'HLT_DiSC30_18_EIso_AND_HE_Mass70'\n",
      "      Physics: 'HLT_Physics'\n",
      "      Random: 'HLT_Random'\n",
      "      ZeroBias: 'HLT_ZeroBias'\n",
      "      AK4CaloJet30: 'HLT_AK4CaloJet30'\n",
      "      AK4CaloJet40: 'HLT_AK4CaloJet40'\n",
      "      AK4CaloJet50: 'HLT_AK4CaloJet50'\n",
      "      AK4CaloJet80: 'HLT_AK4CaloJet80'\n",
      "      AK4CaloJet100: 'HLT_AK4CaloJet100'\n",
      "      AK4CaloJet120: 'HLT_AK4CaloJet120'\n",
      "      AK4PFJet30: 'HLT_AK4PFJet30'\n",
      "      AK4PFJet50: 'HLT_AK4PFJet50'\n",
      "      AK4PFJet80: 'HLT_AK4PFJet80'\n",
      "      AK4PFJet100: 'HLT_AK4PFJet100'\n",
      "      AK4PFJet120: 'HLT_AK4PFJet120'\n",
      "      HISinglePhoton10_Eta3p1ForPPRef: 'HLT_HISinglePhoton10_Eta3p1ForPPRef'\n",
      "      HISinglePhoton20_Eta3p1ForPPRef: 'HLT_HISinglePhoton20_Eta3p1ForPPRef'\n",
      "      HISinglePhoton30_Eta3p1ForPPRef: 'HLT_HISinglePhoton30_Eta3p1ForPPRef'\n",
      "      HISinglePhoton40_Eta3p1ForPPRef: 'HLT_HISinglePhoton40_Eta3p1ForPPRef'\n",
      "      HISinglePhoton50_Eta3p1ForPPRef: 'HLT_HISinglePhoton50_Eta3p1ForPPRef'\n",
      "      HISinglePhoton60_Eta3p1ForPPRef: 'HLT_HISinglePhoton60_Eta3p1ForPPRef'\n",
      "      Photon20_HoverELoose: 'HLT_Photon20_HoverELoose'\n",
      "      Photon30_HoverELoose: 'HLT_Photon30_HoverELoose'\n",
      "      Photon40_HoverELoose: 'HLT_Photon40_HoverELoose'\n",
      "      Photon50_HoverELoose: 'HLT_Photon50_HoverELoose'\n",
      "      Photon60_HoverELoose: 'HLT_Photon60_HoverELoose'\n",
      "      EcalCalibration: 'HLT_EcalCalibration'\n",
      "      HcalCalibration: 'HLT_HcalCalibration'\n",
      "      L1UnpairedBunchBptxMinus: 'HLT_L1UnpairedBunchBptxMinus'\n",
      "      L1UnpairedBunchBptxPlus: 'HLT_L1UnpairedBunchBptxPlus'\n",
      "      L1NotBptxOR: 'HLT_L1NotBptxOR'\n",
      "      L1MinimumBiasHF_OR: 'HLT_L1MinimumBiasHF_OR'\n",
      "      L1MinimumBiasHF0OR: 'HLT_L1MinimumBiasHF0OR'\n",
      "      HcalNZS: 'HLT_HcalNZS'\n",
      "      HcalPhiSym: 'HLT_HcalPhiSym'\n",
      "      HcalIsolatedbunch: 'HLT_HcalIsolatedbunch'\n",
      "      IsoTrackHB: 'HLT_IsoTrackHB'\n",
      "      IsoTrackHE: 'HLT_IsoTrackHE'\n",
      "      ZeroBias_FirstCollisionAfterAbortGap: 'HLT_ZeroBias_FirstCollisionAfte...'\n",
      "      ZeroBias_IsolatedBunches: 'HLT_ZeroBias_IsolatedBunches'\n",
      "      ZeroBias_FirstCollisionInTrain: 'HLT_ZeroBias_FirstCollisionInTrain'\n",
      "      ZeroBias_LastCollisionInTrain: 'HLT_ZeroBias_LastCollisionInTrain'\n",
      "      ZeroBias_FirstBXAfterTrain: 'HLT_ZeroBias_FirstBXAfterTrain'\n",
      "      Ele24_eta2p1_WPTight_Gsf_LooseChargedIsoPFTau30_eta2p1_CrossL1: 'HLT_E...'\n",
      "      Ele24_eta2p1_WPTight_Gsf_MediumChargedIsoPFTau30_eta2p1_CrossL1: 'HLT_...'\n",
      "      Ele24_eta2p1_WPTight_Gsf_TightChargedIsoPFTau30_eta2p1_CrossL1: 'HLT_E...'\n",
      "      Ele24_eta2p1_WPTight_Gsf_LooseChargedIsoPFTau30_eta2p1_TightID_CrossL1: 'HLT_Ele24_eta2p1_WPTight_Gsf_LooseChargedIsoPFTau30_eta2p1_TightID_CrossL...'\n",
      "      Ele24_eta2p1_WPTight_Gsf_MediumChargedIsoPFTau30_eta2p1_TightID_CrossL1: 'HLT_Ele24_eta2p1_WPTight_Gsf_MediumChargedIsoPFTau30_eta2p1_TightID_Cross...'\n",
      "      Ele24_eta2p1_WPTight_Gsf_TightChargedIsoPFTau30_eta2p1_TightID_CrossL1: 'HLT_Ele24_eta2p1_WPTight_Gsf_TightChargedIsoPFTau30_eta2p1_TightID_CrossL...'\n",
      "      DoubleLooseChargedIsoPFTau35_Trk1_eta2p1_Reg: 'HLT_DoubleLooseChargedI...'\n",
      "      DoubleLooseChargedIsoPFTau40_Trk1_eta2p1_Reg: 'HLT_DoubleLooseChargedI...'\n",
      "      DoubleMediumChargedIsoPFTau35_Trk1_eta2p1_Reg: 'HLT_DoubleMediumCharge...'\n",
      "      DoubleMediumChargedIsoPFTau40_Trk1_eta2p1_Reg: 'HLT_DoubleMediumCharge...'\n",
      "      DoubleTightChargedIsoPFTau35_Trk1_eta2p1_Reg: 'HLT_DoubleTightChargedI...'\n",
      "      DoubleTightChargedIsoPFTau40_Trk1_eta2p1_Reg: 'HLT_DoubleTightChargedI...'\n",
      "      DoubleLooseChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleLoose...'\n",
      "      DoubleLooseChargedIsoPFTau40_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleLoose...'\n",
      "      DoubleMediumChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleMedi...'\n",
      "      DoubleMediumChargedIsoPFTau40_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleMedi...'\n",
      "      DoubleTightChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleTight...'\n",
      "      DoubleTightChargedIsoPFTau40_Trk1_TightID_eta2p1_Reg: 'HLT_DoubleTight...'\n",
      "      IsoMu24_eta2p1_MediumChargedIsoPFTau50_Trk30_eta2p1_1pr: 'HLT_IsoMu24_...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr_MET90: 'HLT_MediumChargedIsoP...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr_MET100: 'HLT_MediumChargedIso...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr_MET110: 'HLT_MediumChargedIso...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr_MET120: 'HLT_MediumChargedIso...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr_MET130: 'HLT_MediumChargedIso...'\n",
      "      MediumChargedIsoPFTau50_Trk30_eta2p1_1pr: 'HLT_MediumChargedIsoPFTau50...'\n",
      "      MediumChargedIsoPFTau180HighPtRelaxedIso_Trk50_eta2p1_1pr: 'HLT_Medium...'\n",
      "      MediumChargedIsoPFTau180HighPtRelaxedIso_Trk50_eta2p1: 'HLT_MediumChar...'\n",
      "      IsoMu24_eta2p1_LooseChargedIsoPFTau35_Trk1_eta2p1_Reg_CrossL1: 'HLT_Is...'\n",
      "      IsoMu24_eta2p1_LooseChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL1: 'HLT_IsoMu24_eta2p1_LooseChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL1...'\n",
      "      IsoMu24_eta2p1_MediumChargedIsoPFTau35_Trk1_eta2p1_Reg_CrossL1: 'HLT_I...'\n",
      "      IsoMu24_eta2p1_MediumChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL1: 'HLT_IsoMu24_eta2p1_MediumChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL...'\n",
      "      IsoMu24_eta2p1_TightChargedIsoPFTau35_Trk1_eta2p1_Reg_CrossL1: 'HLT_Is...'\n",
      "      IsoMu24_eta2p1_TightChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL1: 'HLT_IsoMu24_eta2p1_TightChargedIsoPFTau35_Trk1_TightID_eta2p1_Reg_CrossL1...'\n",
      "      Ele16_Ele12_Ele8_CaloIdL_TrackIdL: 'HLT_Ele16_Ele12_Ele8_CaloIdL_TrackIdL'\n",
      "      Rsq0p35: 'HLT_Rsq0p35'\n",
      "      Rsq0p40: 'HLT_Rsq0p40'\n",
      "      RsqMR300_Rsq0p09_MR200: 'HLT_RsqMR300_Rsq0p09_MR200'\n",
      "      RsqMR320_Rsq0p09_MR200: 'HLT_RsqMR320_Rsq0p09_MR200'\n",
      "      RsqMR300_Rsq0p09_MR200_4jet: 'HLT_RsqMR300_Rsq0p09_MR200_4jet'\n",
      "      RsqMR320_Rsq0p09_MR200_4jet: 'HLT_RsqMR320_Rsq0p09_MR200_4jet'\n",
      "      }\n",
      "    HLTriggerFinalPath: 'HLTriggerFinalPath'\n",
      "    Flag: Record {\n",
      "      name = u'Flag'\n",
      "      HBHENoiseFilter: 'Flag_HBHENoiseFilter'\n",
      "      HBHENoiseIsoFilter: 'Flag_HBHENoiseIsoFilter'\n",
      "      CSCTightHaloFilter: 'Flag_CSCTightHaloFilter'\n",
      "      CSCTightHaloTrkMuUnvetoFilter: 'Flag_CSCTightHaloTrkMuUnvetoFilter'\n",
      "      CSCTightHalo2015Filter: 'Flag_CSCTightHalo2015Filter'\n",
      "      globalTightHalo2016Filter: 'Flag_globalTightHalo2016Filter'\n",
      "      globalSuperTightHalo2016Filter: 'Flag_globalSuperTightHalo2016Filter'\n",
      "      HcalStripHaloFilter: 'Flag_HcalStripHaloFilter'\n",
      "      hcalLaserEventFilter: 'Flag_hcalLaserEventFilter'\n",
      "      EcalDeadCellTriggerPrimitiveFilter: 'Flag_EcalDeadCellTriggerPrimitive...'\n",
      "      EcalDeadCellBoundaryEnergyFilter: 'Flag_EcalDeadCellBoundaryEnergyFilter'\n",
      "      goodVertices: 'Flag_goodVertices'\n",
      "      eeBadScFilter: 'Flag_eeBadScFilter'\n",
      "      ecalLaserCorrFilter: 'Flag_ecalLaserCorrFilter'\n",
      "      trkPOGFilters: 'Flag_trkPOGFilters'\n",
      "      chargedHadronTrackResolutionFilter: 'Flag_chargedHadronTrackResolution...'\n",
      "      muonBadTrackFilter: 'Flag_muonBadTrackFilter'\n",
      "      BadChargedCandidateFilter: 'Flag_BadChargedCandidateFilter'\n",
      "      BadPFMuonFilter: 'Flag_BadPFMuonFilter'\n",
      "      BadChargedCandidateSummer16Filter: 'Flag_BadChargedCandidateSummer16Fi...'\n",
      "      BadPFMuonSummer16Filter: 'Flag_BadPFMuonSummer16Filter'\n",
      "      trkPOG_manystripclus53X: 'Flag_trkPOG_manystripclus53X'\n",
      "      trkPOG_toomanystripclus53X: 'Flag_trkPOG_toomanystripclus53X'\n",
      "      trkPOG_logErrorTooManyClusters: 'Flag_trkPOG_logErrorTooManyClusters'\n",
      "      METFilters: 'Flag_METFilters'\n",
      "      }\n",
      "    }\n",
      "  ]\n"
     ]
    }
   ],
   "source": [
    "print(tree.arrowed.schema().format())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Using this mapping, we walk through the structure of the code and replace object references with array lookups."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": true,
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Before transformation:\n",
      "----------------------\n",
      "def calculateMHT(events, MHT_pt, MHT_phi):\n",
      "    i = 0\n",
      "    for event in events:    \n",
      "        mht_px = 0.0\n",
      "        mht_py = 0.0\n",
      "        mht_pz = 0.0\n",
      "        for electron in event.Electron:        \n",
      "            mht_px += ((electron.pt * cosh(electron.eta)) * sin(electron.phi))\n",
      "            mht_py += ((electron.pt * cosh(electron.eta)) * cos(electron.phi))\n",
      "            mht_pz += (electron.pt * sinh(electron.eta))\n",
      "        for muon in event.Muon:        \n",
      "            mht_px += ((muon.pt * cosh(muon.eta)) * sin(muon.phi))\n",
      "            mht_py += ((muon.pt * cosh(muon.eta)) * cos(muon.phi))\n",
      "            mht_pz += (muon.pt * sinh(muon.eta))\n",
      "        for jet in event.Jet:        \n",
      "            if (jet.pt > 150):            \n",
      "                mht_px += ((jet.pt * cosh(jet.eta)) * sin(jet.phi))\n",
      "                mht_py += ((jet.pt * cosh(jet.eta)) * cos(jet.phi))\n",
      "                mht_pz += (jet.pt * sinh(jet.eta))\n",
      "            \n",
      "        MHT_pt[i] = sqrt(((mht_px ** 2) + (mht_py ** 2)))\n",
      "        MHT_phi[i] = (-atan2(MHT_pt[i], mht_pz))\n",
      "        i += 1\n",
      "    return None\n",
      "\n",
      "After transformation:\n",
      "---------------------\n",
      "def calculateMHT(par0_mem0_None_beginarray, par0_mem0_None_endarray, par0_mem6_nElectron_beginarray, par0_mem6_nElectron_endarray, par0_mem15_Electron_pt_array, par0_mem10_Electron_eta_array, par0_mem14_Electron_phi_array, par0_mem115_nMuon_beginarray, par0_mem115_nMuon_endarray, par0_mem121_Muon_pt_array, par0_mem118_Muon_eta_array, par0_mem120_Muon_phi_array, par0_mem74_nJet_beginarray, par0_mem74_nJet_endarray, par0_mem84_Jet_pt_array, par0_mem81_Jet_eta_array, par0_mem83_Jet_phi_array, MHT_pt, MHT_phi):\n",
      "    i = 0\n",
      "    for event in range(par0_mem0_None_beginarray[0], par0_mem0_None_endarray[0]):    \n",
      "        mht_px = 0.0\n",
      "        mht_py = 0.0\n",
      "        mht_pz = 0.0\n",
      "        for electron in range(par0_mem6_nElectron_beginarray[event], par0_mem6_nElectron_endarray[event]):        \n",
      "            mht_px += ((par0_mem15_Electron_pt_array[electron] * cosh(par0_mem10_Electron_eta_array[electron])) * sin(par0_mem14_Electron_phi_array[electron]))\n",
      "            mht_py += ((par0_mem15_Electron_pt_array[electron] * cosh(par0_mem10_Electron_eta_array[electron])) * cos(par0_mem14_Electron_phi_array[electron]))\n",
      "            mht_pz += (par0_mem15_Electron_pt_array[electron] * sinh(par0_mem10_Electron_eta_array[electron]))\n",
      "        for muon in range(par0_mem115_nMuon_beginarray[event], par0_mem115_nMuon_endarray[event]):        \n",
      "            mht_px += ((par0_mem121_Muon_pt_array[muon] * cosh(par0_mem118_Muon_eta_array[muon])) * sin(par0_mem120_Muon_phi_array[muon]))\n",
      "            mht_py += ((par0_mem121_Muon_pt_array[muon] * cosh(par0_mem118_Muon_eta_array[muon])) * cos(par0_mem120_Muon_phi_array[muon]))\n",
      "            mht_pz += (par0_mem121_Muon_pt_array[muon] * sinh(par0_mem118_Muon_eta_array[muon]))\n",
      "        for jet in range(par0_mem74_nJet_beginarray[event], par0_mem74_nJet_endarray[event]):        \n",
      "            if (par0_mem84_Jet_pt_array[jet] > 150):            \n",
      "                mht_px += ((par0_mem84_Jet_pt_array[jet] * cosh(par0_mem81_Jet_eta_array[jet])) * sin(par0_mem83_Jet_phi_array[jet]))\n",
      "                mht_py += ((par0_mem84_Jet_pt_array[jet] * cosh(par0_mem81_Jet_eta_array[jet])) * cos(par0_mem83_Jet_phi_array[jet]))\n",
      "                mht_pz += (par0_mem84_Jet_pt_array[jet] * sinh(par0_mem81_Jet_eta_array[jet]))\n",
      "            \n",
      "        MHT_pt[i] = sqrt(((mht_px ** 2) + (mht_py ** 2)))\n",
      "        MHT_phi[i] = (-atan2(MHT_pt[i], mht_pz))\n",
      "        i += 1\n",
      "    return None\n",
      "\n"
     ]
    }
   ],
   "source": [
    "tree.arrowed.run(calculateMHT, (MHT_pt, MHT_phi), env=math.__dict__, debug=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "This is a formal process, like a compilation pass. To be rigorous, we need to do type inference through the code tree."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "Finally, the transformed code, which consists of nothing but numbers and arrays of numbers, is passed to Numba, a numerical Python compiler.\n",
    "\n",
    "<table>\n",
    "<tr style=\"background: white\"><td><center><img src=\"http://numba.pydata.org/_static/numba_blue_icon_rgb.png\" width=\"200\"></center></td></tr>\n",
    "<tr><td style=\"font-size: 18pt\">Numba <a href=\"http://numba.pydata.org/\">numba.pydata.org</a></td></tr>\n",
    "</table>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Even without the full compilation, the transformation by itself speeds up access considerably."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Time for one run: 11.898900032 seconds\n"
     ]
    }
   ],
   "source": [
    "startTime = time.time()\n",
    "tree.arrowed.run(calculateMHT, (MHT_pt, MHT_phi), env=math.__dict__, numba=None)\n",
    "endTime = time.time()\n",
    "\n",
    "print(\"Time for one run: {0} seconds\".format(endTime - startTime))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Compilation results in fast code, but it's hard to debug. The thing we like about Python is the ease of *exploring* data. For interacting with individual objects, use proxies:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [],
   "source": [
    "events = tree.arrowed.proxy()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<Record at index 0>, <Record at index 1>, <Record at index 2>, <Record at index 3>, ...]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "events"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "fragment"
    }
   },
   "source": [
    "A proxy pretends to be a Python object, fetching the data you ask for on demand. Every time you look at a particle attribute, it reads one TBasket from the ROOT file or reuses a previously read one.\n",
    "\n",
    "Compilation is for high throughput; proxies are for low latency."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Compiled functions can return proxies."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [],
   "source": [
    "def event_with_bigjet(events):\n",
    "    for event in events:\n",
    "        for jet in event.Jet:\n",
    "            if jet.pt > 1000:\n",
    "                return event\n",
    "\n",
    "event = tree.arrowed.run(event_with_bigjet)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[1030.0,\n",
       " 628.0,\n",
       " 277.5,\n",
       " 140.25,\n",
       " 92.125,\n",
       " 34.84375,\n",
       " 32.21875,\n",
       " 23.6875,\n",
       " 21.9375,\n",
       " 21.84375,\n",
       " 20.78125]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "[jet.pt for jet in event.Jet]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "You *can* use existing code in these compiled functions. If the code references other Python functions, it will (someday) extend the compilation process to include those other functions. It can also call functions defined in C.\n",
    "\n",
    "Here's the catch: C, not C++.\n",
    "\n",
    "C functions can be loaded from `.so` files without name mangling and executed as pointers, needing only the argument signature. C++ files require more vital information from the header files."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "**Solution:** use ROOT's Cling to make C functions on demand."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [],
   "source": [
    "import ROOT\n",
    "import ctypes\n",
    "\n",
    "# we want a ROOT histogram\n",
    "h = ROOT.TH1I(\"h\", \"\", 150, 0, 150)\n",
    "\n",
    "# make a function (to get a function pointer) that fills h\n",
    "ROOT.gInterpreter.Declare(\"void h_Fill(double x) { h->Fill(x); }\")\n",
    "\n",
    "# turn it into an integer and pass it from C++ to Python\n",
    "ROOT.gInterpreter.Declare(\"size_t ptr_h_Fill = reinterpret_cast<size_t>(h_Fill);\")\n",
    "\n",
    "# use Python's ctypes to interpret that integer as a callable function\n",
    "h_Fill = ctypes.CFUNCTYPE(None, ctypes.c_double)(ROOT.ptr_h_Fill)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [],
   "source": [
    "def process(events):\n",
    "    for event in events:\n",
    "        for jet in event.Jet:\n",
    "            h_Fill(jet.pt)\n",
    "\n",
    "tree.arrowed.run(process, env={\"h_Fill\": h_Fill})"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3da3KruKIGUHHrzAs8mT6zaONh9J5Ll2Fk3B+60VUAOSTBD+S1qqvLwYBRnG191otmmqYA\nALDmf559AQDA6xIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIU\nAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAi\nQQEAKBIUAIAiQQEAKBIUAIAiQQEAKPrPsy8AgNo0TfPsS+Br0zRt2U1QAGB/GyshnmV7mNP1AAAU\nCQoAQJGgAAAUCQoAsK5pmr7vn30VTyYoAEDRMAzPvoQnExQAgCJBAQBu6fu+67q37YZoTHUFYF9N\nU0nlkhYbaNt2HMf4oI7OiO3vUSXvJQCvo7KgkMoy+/HQtr9Huh4AoKht29XH70NQAACKBAUAoEhQ\nAACKBAUAoEhQAACKKpnBkmy/wTYA91NZ5VKf7dMj/3PvS3k8f50Az+U7W00qDAoAPJ2sUA1BAYD9\nadx9cduTnMGMAECRoAAAFFXY9VBqTtEOBvB0q7de7LruNyf8zeF8qcIWhang2dcF8O66rjut+fLA\n0p2d4wl3vsr9NE3T9336MWaapmm6rksl6vu+WUjRJz2bHzJ7Kn+JG9t/rMIWBQBe2fV6/db+fd9f\nLpfV73uv3JYwu7ZhGGKmOZ/P8fH5fI51+eymlOM4xgex4PHZcRxPp1P6JXRdN45jfOpyucSdw0eb\neto+DEMpY31D6fv3QdVXIoDDKX0Uxwrsy8Ov1+v1ek0/ns/nEEK+Je6z5dgvt99JvOYYC+KWWdlL\nVXA8MO3Ttm18HNNVPFt8nIqTzjz7ReWnWtpeXdZWrQoKAE/3s6Cw/G49ZTVurALbtk37xB/TCWcN\nFbP6crn9rlJZ8qCQav2pHBTyK5xdbTpb6de43J5fwPKFNpZF1wMAD7XsO09bxnG8Xq9d18VW92EY\n4lOXyyVuj/u0bdv3/axtPzbsx7jQ931qqI+t9/E8cUjEdOdRa7H9fxiGfHB9PighPpglmPDRW5HK\nFa9z+JA/FUsUOx1SF0bsj5i9oq6HufpKBHA4pY/iWYPBrCYKWUv7lH0bzpvQZ1+aZ63us2uIh4fP\nDfj3blGY9R3MvtDnzR6zA2cdCkn6paVSpDOk9pX4VDxD27axmLOjZrZXlxXOegDgla1WWtGPByfG\n7835xIG0sW3bcRzTxIG7jn8chuFyuSybCpKu61KNPpvMH1tElpc3DEPMH+M45s9O0xRbGuJTsWjx\ncWw4iRni9+UVFACoxPmz2CA/DEMczTCO4+VyuetNKGI0ia8SXyi9Yj4BIdboIesXiA9mCSPvNYhl\nyXsWkjwK9H0/fYz03KHTIYRgeiQAFYjd8/noh67rYhNC3/d938daMw19uFO7wuy0cXhE3Bi/4s8q\n77R/urz0VJxCmUZm5GLoyfdMD2Jh4yFp4MJvShTC3j36KQ3FPpK4cTllNu8ySYcsp74sd45Sl8yy\nL2f3EgHwXaWP4jQRYCZ+mIfPPfrhO2MUpo+xCLF7frl9mqbZ9gcIi+mRsaSpgsv3XF5YXjnGX0I+\nPTIfl7D6S1itQPOTby3Fxv22SPV3PndlyoqUzH5x6ZDZbJbZMI3/u+KP9FD6te5YIgB+4HZQWEqj\nDleDQqoLlzX9D6ZHlqYL3sOyRLnZnstKfVai1e/Ys1Plh9xICdOzgsKyRo8/ltZ8yJePmNaCYXyc\nryCxXLBi9pYfJSiE8E/679nXArCze3wUb5yt8CILLpVcP3z3kG89teUlnhAUlvM6UsxJtfvsumcB\nIgWCWYCYCrNcpsLiEnuV6K5SPhAUgPoc5aP4nW1/j3ab9RCnfMzGZeRtI03TnE6n/F4XpTEdqz/u\nNXoTANjuLtMj45JY4fN0jjhcJZ8JujrN48s7WOQBYnUw5/I2XF/6cUkBoG47T49Md8fKJ4FMnxfT\nGIYhRoTZBI98h+0vt9w4LQaMAAA/s2eLQkoJ1+v1RmX/5ZzO2zvogwCAh9kzKMTVJPKRCuHjrhiz\nRSTig9VRCHGJjLBYrypunDVCpLtxA/D6bnQu+xL4uvYaP5mvfDBbLyG+0HLJiOnzOhghm9GQL6uw\nnCabL8Exm0mxY4nuyqwHoGKlj+JYHSxXMihtfwexusy3pHkAN34hq8sbLne48brbq8vdqtXVb/bL\nW12Fz/MbZ6tJfLqyTP5byGdSrP617VWiuxIUgIrdDgqzejFVBG8YFNJ37LRltg7h6qJJeUpY/c6/\nXPlx6QlBYYsdl4wobRcUAJ7udlCYPVv6Av2zquEVVlXa6Lp2p4JQWIcwN/sdLn91q7/n5Uk2Xucx\nqtXtBAWAp7sRFNLtHm5sLK1PPNueMkH+3frLCvJ1pMuerUS88cD8x/z3mc65V1Bwm2kAHm02XH02\n4P1yueR3eUgj3OPNGNMCvnGeXTSOYz58LT/ha5oN24/Sb6Prutk8gFxakajv+7gUUJoc0Pf9OI7X\nxb0Yf2VjoDiKo5RIiwJQsdJHcQghjnNPO6THIVuqPyxaEaZpikv2pe2z2wOFz03xt++H9HT5YPyw\nGMgfPt9ecfUMeT2eDs/H++/YonCManU7QQHg6W4HhWntDj6zoDCT3woy72WYnXb1xxcU1mb/TZ8n\nCU6FyX3T52CUzxlcZo7b17DxanU9APBocfXeZb9Dcv4sfKzKc7lcuq5L95s+oljey+WS7iEwjmPT\nNMMwxELl/Qih0D2ROhe6rosrDMXt8VRN08Q1h270X2y38xLOAPClrusul0usw1ZX403V2zAMfd/H\nbvtQxSL9s/JeLpcQwvl8jusN5osK5usNbjltPtgzP+1vr3hjy8NRHKVEuh6AipU+isOimT3/MR/A\nmD+OZ8sX4putExCO1vWQW734NGAzlXE5JiNfxnBZ3h27HipsUSjdDXI6fg4FqEZsMF82jMdvxpfL\nJX4nDh+f3rHrIc10iIc3TVPZZ/v1ej2dTqmYq/MXZvvkd2G8h9p+xUf5o2maP9P0V/4AoBq7fBSv\nNrznG7/VMn8sW4r2y+Jvf4+OUa1uJygAPN1RPorf2fb3yKwHAKBIUAAAigQFAKBIUAAAigQFAKCo\nwnUUAHhxe81szNcPqHKe5CvQogDA48QbKMf1gpqmyWv3dPeH0lHL7adMvMfB699gOkr3ku66blnq\nuLrUjcPjDqXD4/lvn+EbNq7geBRHKZElnIGKlT6K47rC6caJs+WHQ3np5dKCxCFb8DjdWPL1129O\nBU93Z5jdIvL2AszpqHRnrNnhs+WfV22vLo9RrW4nKAA8XemjeFmLz25hMHs2JoDpq6CQb1m9NXM6\nz4uY1eJhcQuM29X88vDZzlvaArZXl7oeAHic/O6IIYRhGOLX3/yGyyG7s0PsU5gddUPsekgdELEF\nPp0nbpx1eZzP591a6TfL7/Q4256aRjYePts5Fm3P23BvDBRHcZQSaVEAKlb6KM7bzFMHRH5UWLtB\nYqrzVl9o2dGQNi5vRBm3z9onVk/yGCkTfKvrYSa//tSgsuPdI49RrW53lEgkKAAVu/GRm+qwJNWR\nswo+rzu/GxRCdmfqtD1WotM0/fvvvyGEv//+Oz3+999/f1TQX8lbBWZPbQwKaSxCfvftGIwEhaJX\nCwQlggJQsY2V3KypINX6qUbPT/iDFoXVL42xTi21Ljxe3tSRbLmqZUNLfp4dg4IxCgA8wjAM+ciA\nOK8vZoLd5zTOhiDk4sa///47jnsYx3HP7vxtZss/tG27fRBG+BjAMY5jCj3pnHGQRxrVsc980Y2B\n4iiOUiItCkDFSh/FIZsbGeW9DGGProfZNIr8qNiGER/HHoe48+P7HWa/h+W3/9vtAaXfxjIPnc/n\n0nSP7dXlMarV7QQFgKcrfRSndRTSygd5nbd8HHe7PZgxLGYKpDp4OZgxfG6of9YX5lS6abGYRHQj\nOqSCtJnlyFBjFIoEBYCnu/FRvGznT0/luSFVh6lSvBEU8t1mX6BnsxDzp/7+++/wMaTxwZalm+1w\nIyiU5lV+eYaZ7dVlM92cKXA4TXOMEjXNn2n6K38AUI0vP4pL93qYbd/3lhCvdjOI517V9uryGNXq\ndoICwNMd5aP4nW1/j8x6AACKBAUAoOg/z74AACr0+LsncCeCAgA7e+cBCvWNz9D1AAAUCQoAQJGg\nAAAUVThGoTSCprJOIwB4gAqDgkAAAHvR9QAAFAkKAECRoAAAFAkKAECRoAAAFAkKAECRoAAAFAkK\nAECRoAAAFAkKAECRoAAAFAkKAEBRhTeFcvdIANhLhUFBIACAveh6AACKBAUAoEhQAACKKhyjcDhN\n8yc9nqa/nnglADAjKDxZngzyxAAAr0DXAwBQJCgAAEU7B4W+75umaZqm67phGNL2YRi6rovbS4fk\n+8dD0qlmh6RTzQ4BAPbV7Lg8Udd14zi2bRtCGMcxhHC9XmM1H1dLbNs2bk8vOjsk7d/3/eVySdvb\ntk2ZoHSq9OwhFlxqmj/LcYurGwE4kKNUQ98w7SeE0Lbt8sdY3+fbz+fzNE3X6zU9nu2Wn+p8PocQ\nrtfrjVPlW3Ys0f2E8M/GjQAcyFGqoe1263qI3/j7vs83xi/9qc0gatv2crksD4ltCcMwxO2pxyHu\nEDeWTgUA3MNuQaHrummaUu0e6/XYGBCyWj83G2Ew22f2o+EIAPB4d5n10Pf96XQKiwaGmdjeMJNa\nFEryALGaP5rv21IoAHhDOweFOFUhjkOcvhrNkXciJKt1/42XW278QQfM9lcEgLeyZ1AYhiE2JFyv\n1+Vcx+3nuZ0V9EEAwMPsGRROp1NsSFjW9HkvQxqQuDoKoeu6NKpxtj1ksyJnpwIA7mG3oJBPXkji\nxjgNMj5e/j/FgrRwQviYzpBPi7j9fwDgLvaaZ7n6zX62FkKUr3wQM8TqxeTb4yIKt0+VjtqrRHdl\nHQWAKh2lGtruoQtIxYWcV7eHtaEJpe03TnWUJbGszAhQpaNUQ9tVV56DvEOCAkCVjlINbefukQBA\nkaAAABQJCgBAkaAAABQJCgBAkaAAABT959kXsL/S3SArm68CAA9QYVAQCABgL7oeAIAiQQEAKBIU\nAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAIAiQQEAKBIUAICiCu8e6TbTALCX\nCoOCQAAAe9H1AAAUCQoAQFGFXQ+H1jR/0uNp+uuJVwIAQVB4KXkyyBMDADyLrgcAoEhQAACKBAUA\noEhQAACKBAUAoEhQAACKBAUAoEhQAACKKlxwyd0jAWAvFQYFgQAA9qLrAQAoEhQAgCJBAQAoEhQA\ngCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgKIKbwrl7pEAsJcKg4JAAAB70fUA\nABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABTdJSg0TTMMQ/pxGIbms67r0rN938eN+SH5UfnO\nUdd1cfvsEABgX/svuNT3/er2tm3T41T3d103jmN86nQ6Xa/X+FTf95fLJW4fxzHPBHHhxbZtx3E8\nnU6WVwKAO5r2cz6f02mv1+ts+3L/6/UaQjifz/HHGAvi4xBC27b54fGE+T5xt3R42rJXce4qhH9+\nuQMAL+go1dB2e3Y9dF13Pp/zloMoNQYsOxdC1gIR2xKGYYjbU6tD3CFuTM0PUdu2l8tlvxIAAJ/s\nHBT6vi91PTRNczqd8jEHs9wwG4sw+9FwBAB4vEfMehjHMYRwPp+v1+v5fI5jDtL2mdSiUJIHiOU4\nxxBC832/KR0AVOwRd4+csvGGcVhijAhxQOJs52/NZVjdczK8EQB28oR1FFabAbbv8D59EE3zJ/33\n7GsB4E3dPSjE5RDygQuppl8dhdB1XRrVONseFo0Qs7GNNZmmv9J/z74WAN7Y7vMo4qTHfHpkfKG4\nJc51THMaw8c0yHhUmhIZq//llMh8RmU+bTJ/rd1LdA/fmv1oqiTAURylGtruQWMU4pSH+GPbtqmB\n4Xq9xqkQ8ce8FSE/JOaD8DED83K5xFmR5/P5y44MAODHmumBQ/+GYVit12cLJ3y5/capmuahJfqx\npvmzvU/hWzsD8ERHqYa2q648B3mHBAWAKh2lGtrO3SMBgCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAo\nEhQAgKJHrMz4YKXbRlc2sRUAHqDCoCAQAMBedD0AAEWCAgBQJCgAAEWCAgBQJCgAAEWCAgBQJCgA\nAEWCAgBQJCgAAEWCAgBQJCgAAEWCAgBQVOFNoaq8e2TT/IkPpumv514JAG+lwqBw6ECwKoWDFBcA\n4DF0PQAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFBU\n4U2hqrx7JAA8RYVBQSAAgL3oegAAigQFAKCowq6HujXNn/R4mv564pUA8A4EhSPJk0GeGADgTnQ9\nAABFggIAUCQoAABFggIAUCQoAABFggIAUCQoAABFggIAUFThgkvuHgkAe6kwKAgEALAXXQ8AQJGg\nAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAU3SUoNE0zDEO+ZRiGruuapum6brZz3/dN06we\nErcvD0mnmh0CAOxr/5UZ+75fbjydTiGEtm3HcWyaJi2e2HXdOI5t28Z9rtdrjAV9318ul7h9HMc8\nE8QVmuOpTqfTO6/D2DR/0uNp+uuJVwJAtab9nM/ndNrr9Zq2x/o+/RhCOJ/P0zRdr9f0eLZbCKFt\n2/y08YSlU+VbdizR/YTwz8ueDYAfO0o1tN2eXQ9d153P51iX51KbQdS27eVyCSHERoLUAhHbEoZh\niNtTj0PcIW4snQoAuIedg0Lf96tdD8txBuGj7i/tM/vRcAQAeLxn3j1yHMflxi8DQR4g4hCH2Q6l\n20zfML3xQAcAuOGZ0yOXnRSh0PZQspoqftAB89MSAEDlHhQUvtVxcDsr6IMAgId5UFDIOwjSgMTV\nUQhd16VRjbPt4WNW5PJUAMA9PCIoxGmQcZDj8v8pFqSFE8LHdIZ8WsTt/wMAd7H3fMv/Wx0hX0dh\n+rzEQr7yQdx59WLy7fnZSqdKR+1cnvuwjgJAlY5SDW33/4skPkBcyHl1e1gbmlDafuNU+bKPr6xp\n/uy4luK+ZwPgx45SDW1XXXkO8g4JCgBVOko1tJ27RwIARc9ccIkdpRtEaVoAYEeCQg1SOMjvJwkA\nv6frAQAoEhQAgCJBAQAoqnCMQunukZXNVwGAB6gwKAgEALAXXQ8AQFGFLQpvLp8haU0FAH5JUKhK\nngysqQDA7+l6AACKBAUAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKKlyZ0d0jAWAv\nFQYFgQAA9lJhUCBxgygAfklQqJYbRAHwewYzAgBFggIAUCQoAABFggIAUGQw47swAwKAHxAU3oIZ\nEAD8jK4HAKBIUAAAigQFAKBIUAAAiioczOjukQCwlwqDgkAAAHupMCjwpTRD0oIKANwmKLydFA4s\nqADAlwSFt2a5RgBuExTel+UaAfiS6ZEAQJGgAAAUCQoAQJExCvwfAxsBWBIUCMHARgAKdD0AAEVa\nFFhh6UYAIkGBOUs3ApDoegAAiipsUXCbaQDYS4VBQSDYkTmTAG+uwqDAXmZzJoUGgDckKLCJhRYA\n3pPBjABAkRYFfkI3BMCbEBT4Nt0QAO9DUOC3LOMIUDFBgV+xjCNA3R40mLHv++azYRjiU8MwdF3X\nNE3XdaWj0s7pkLh9eQgAsKMHtSjEmr5t2+VTp9MpPjWOY9M0abmkruvGcYyHnE6n6/UaY0Hf95fL\nJW4fx7HrulmM4FmMcASo0PQQIYS2bZfbY32f73Y+n6dpul6v6fFst/xU5/M5hHC9XvMz7H/1dxDC\nP8++hDsK4Z/8v2dfDsDjHKUa2u7R6yjMvv2nNoOobdvL5ZJ26/s+bo9tCcMwxO2pxyHuoEXh1UzT\nX+m/Z18LAL/yuKAQexZOp1PTNCkBhKzWz83q/tk+sx8FhRcXl3822hHgiB4xRiFV5LFDIQ4y6Lru\nxlDEcRxvnOe20t0jb5jcR+puTIsAOLRHBIWu6/KaOM5Z6Pv+RsUfxzYuz7MlK6j1X5bRjgCH87R1\nFFIO+FbHQdd1cRADh+NelABH9IgxCnE5hNnGNIYxbzlIYxtXRyGk3op8DYblzrw+ox0BjuIRQSFW\n5LHjIC6vFD4mLKRRC6v/T7EgLZwQPmZG5NMi8qGRAMCOmsf06A/DEBdWis7nc6rd49jG5fbZIfl1\n5u0TaSGm9NQhxig0zR9fpqPZIEe/FuDQjlINbffQ8tzoKUgtDRsPKW0/yjskKKwqzYzwuwKO4ijV\n0HbVlecg75CgsJ0mB+BAjlINbefukby60nQJiQHgAQQFjsTyTQAPJihwVAY0ADyAoMAhldKAlgaA\nfT367pEAwIFoUaA2FocG2FGFQaF098jK5quw6sYdJVb3AeC2CoOCQEC0GggMYgD4FmMUAICiClsU\n4Db9EQDbCQq8F/0RAN8iKEAI5koAFAgKsGmuRBAggLckKMAn1nwEyAkKsJUbVwJvSFCATfIbVxrQ\nALwPQQG+Zzag4YlXAvAAggL8yreyguYH4HAEBfi5b1X8mh+AIxIU4HEMbgAOR1CABymt1iAxAK+s\nwqDgNtO8PnMogKOoMCgIBByIORTAi6swKMBxubMl8GoEBXgVpTtb6psAnkhQgJd244ZVQgPwAIIC\nHMaWu1xKD8C+BAU4JP0UwGMIClAP/RTA7gQFqNPqxMtZb4X0AHxJUIC3sFwIMg8NAgRQIihA/UoV\n/2rfhHWfgJygAG/qRrOB+RRAIigAn2yZT3F7Z6AmggLwtVJ6ePyVAA9WYVBw90h4mFJW0NIA1agw\nKAgE8Bg3xkguJ1kAB1VhUACeK59AockBjk5QAO5lS5PD6m7WdYDX0VTWUN80xyhR0/zx2QdhwygH\nC1FzLEephrbTogA805d1v5Wg4LkEBeAwjHiAxxMUgGPYOOLhy/2BbxEUgGP71lKSpf2BEkEBqJDm\nB9iLoAC8Ec0P8F2CAvDuvtv8sOVYqIagALDuyxCgI4N3ICgA/NB3OzJuHwivqcKg4O6RwBP9uB1i\ny7HweBUGBYEAeGWGRHAsFQYFgCP6TVPEz04IWwgKR1LfvUZm6i6g0h3Xi5TuWxV/ihTbssV/X6GA\nd/Iib99xCQoAdUoRYUNbxX/1elBSW846SnL82W2mj1K6H6u7gEp3XHWXLmwo4LfuBl7yrKjx4Lev\nvr+W6spzkHdIUFhVdwGV7rjqLl34aQFn4eD2Z9pvokbc51sv9/mlBYVfqa48B3mHBIVVdRdQ6Y6r\n7tKF5xVwS92f7/NlqiicQVD4lerKc5B3SFBYVXcBle646i5dqKWA5TaJ/4bwT3y0ZWrJ6vaXbcB4\ngEOWp+/7y+USQmjbdhiG/Kkfv0M/bXnzcs8/0Mu9woFe7rgv9+MDD/dyGyeXfvR0/P/LfWtWan1T\nSI4XFGJKaNs2hDCO4ywrHO4P18v9/kAv9woHernjvtyPD/Ry+x74so5XnqZpUjiIoSEvQt1/SXW/\n3I8P9HKvcKCXO+7L/fhAL7fvgS/rf559Ad+T8kH8MT5IPwIA+zpYUFg1G6YAAOzlYCszxkzQdd2N\nfZrC3SO/9LMDvdwrHOjlXuFAL3fcl/vxgV7uHRwsKETDMJSyQmU9QwDwXAfreljNB7cbGACAHztk\nUEiDErb0RAAAP3a8WRxd143jeL1eu66LHUiHKwIAHMXBWhTCRyvC6XSKKeF6vcbtfd83TdM0TR0N\nDHlx8mkdcXxGNcUMITRNU18BS6Woo3Slf2tHL93sTzHcLFH6JRxo1tXtAs6eOlwBb1xq13Wzd/Bw\npXuy6Ziu1+v1ek0/ns/nEELbtnHFxrZtn3Zle0iliA9CCKmw8ce0/amXuYNZ6aYqCpjC67IUFZRu\n+ceZnjp06eJnSP6nOJVLNPslzI56TcsCLv9Qz+dzfOpwBVx9+/Kn8krhcKV7uoP9Yy7J/w7in8VT\nL+e3Zn/W6cf4Z51vT/+wjyi+U/k/1DoKmL998bM4lqKa0qXLjqWLb99xS5f+Dmd1RqlE+Xu63O0F\nbSxg+vFYBSyVLsrDUL7lKKV7ETX8gvJPq+hAH1JLq8UJH0Mxlrn44Re4j1jMWaKvoIDLty+poHRT\nOSgct3TX6/V8Pi+/XJZKNPsqcuO77IvYWMBUrmMVsFS6KH3IlL5JvnjpXsTxxihsdNyep67rpmlK\nPWqxICk1H7T3d+l0OrVtu1x+u44Clvq2Kyhd27aXy6Xrur7vT6dTyAp10NLFsqyuBL9aotlny+uX\nulTAaZryssRb8oajFfDLt285LGO5D7fVEBQqniSZPosru5/F6j/gOsRCpdugj+NY2eJu8U9xHMdY\nxrzh902M47jceOg/5mEY8rHhdRSw7/s4P262vY7SPVgNQSGq7J2O/3TjDbWnuuZ/lv4B1+R8Pg/D\nMAxDfO9qSrGxKSg2SJ7P58vlUlmK/VIa+pc77lvcdV38NhLnnIdaCni5XM7n8/Ky6yjdg9UQFOpb\nrnEYhvRP93a72RHlE1xjMdNk13D8AsY/vNmfX/oSc/TSxUyQShF/LDVZV6C+Es00TTOOYxx0cujP\nzJn0lxnnQIYQYtte9W/onRzyXg8zqR07b9A+9B99/NK2+jedt5uN47iajl9c3/f52pqxFOn9OnoB\nZ3+NUSrF0Uu3quLSrZYorvmWth/3Ayde87LBsoICzq42dZPFBRWOXroneMoQyt3lQ16PXq58TYhk\nOS+rjsG6szkCdRQw/2vMH9dRuvjHmZeujj/O5XSVGyUKH6Po0+SdR1/u95WmU80+atJTxyrgjdlG\nU3nC+VFK93QHrlBn8vRzrE+omdXvYbO5PdFxp4Amy3/edRQwf+/yUlRQuuGGu88AAADaSURBVNng\nkmpKt1rTlEo0+yU8+FJ/ZjWUL60++7SL3uxbQeFwpXu6493r4Yb3aUS6caPtOlRQwBt/jUp3LKUS\nVf+BU3cB6y7dvqoKCgDAvmqY9QAA3ImgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoA\nQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGg\nAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAU/S8e3MOfukra9wAAAABJRU5E\nrkJggg==\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "c = ROOT.TCanvas()\n",
    "h.Draw()\n",
    "c.Draw()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [],
   "source": [
    "ROOT.TLorentzVector   # cause ROOT to dyload the library\n",
    "ROOT.gInterpreter.Declare(\"\"\"\n",
    "double mass(double pt1, double eta1, double phi1, double pt2, double eta2, double phi2) {\n",
    "    TLorentzVector v1, v2;\n",
    "    v1.SetPtEtaPhiM(pt1, eta1, phi1, 0);\n",
    "    v2.SetPtEtaPhiM(pt2, eta2, phi2, 0);\n",
    "    return (v1 + v2).M();\n",
    "}\n",
    "size_t ptr_mass = reinterpret_cast<size_t>(mass);\n",
    "\"\"\")\n",
    "\n",
    "mass = ctypes.CFUNCTYPE(ctypes.c_double,                                    # return type\n",
    "                        ctypes.c_double, ctypes.c_double, ctypes.c_double,  # parameter types\n",
    "                        ctypes.c_double, ctypes.c_double, ctypes.c_double\n",
    "                       )(ROOT.ptr_mass)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [],
   "source": [
    "def calc_masses(events):\n",
    "    i = 0\n",
    "    for event in events:\n",
    "        e1 = None\n",
    "        e2 = None\n",
    "        for electron in event.Electron:\n",
    "            if e1 is None or electron.pt > e1.pt:\n",
    "                e2 = e1\n",
    "                e1 = electron\n",
    "                continue\n",
    "            if e2 is None or electron.pt > e2.pt:\n",
    "                e2 = electron\n",
    "        if e1 is not None and e2 is not None:\n",
    "            h_Fill(mass(e1.pt, e1.eta, e1.phi, e2.pt, e2.eta, e2.phi))\n",
    "        i += 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "scrolled": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "10000 events / 1.80 seconds = 5.5 kHz\n",
      "20000 events / 1.84 seconds = 10.9 kHz\n",
      "30000 events / 2.79 seconds = 10.8 kHz\n",
      "40000 events / 2.82 seconds = 14.2 kHz\n",
      "50000 events / 3.76 seconds = 13.3 kHz\n",
      "60000 events / 3.80 seconds = 15.8 kHz\n",
      "70000 events / 4.78 seconds = 14.6 kHz\n",
      "80000 events / 6.05 seconds = 13.2 kHz\n",
      "90000 events / 6.09 seconds = 14.8 kHz\n",
      "100000 events / 7.07 seconds = 14.1 kHz\n",
      "110000 events / 8.02 seconds = 13.7 kHz\n",
      "120000 events / 8.08 seconds = 14.9 kHz\n",
      "130000 events / 8.11 seconds = 16.0 kHz\n",
      "140000 events / 8.97 seconds = 15.6 kHz\n",
      "150000 events / 9.85 seconds = 15.2 kHz\n",
      "160000 events / 9.89 seconds = 16.2 kHz\n",
      "170000 events / 11.02 seconds = 15.4 kHz\n",
      "180000 events / 11.06 seconds = 16.3 kHz\n",
      "190000 events / 11.10 seconds = 17.1 kHz\n",
      "200000 events / 11.98 seconds = 16.7 kHz\n",
      "210000 events / 12.01 seconds = 17.5 kHz\n",
      "220000 events / 13.70 seconds = 16.1 kHz\n",
      "230000 events / 13.73 seconds = 16.8 kHz\n",
      "234805 events / 13.76 seconds = 17.1 kHz\n"
     ]
    }
   ],
   "source": [
    "tree = uproot.open(\"~/storage/data/SingleElectronRun2015D_16Dec2015_v1_0000_1.root\")[\"Events\"]\n",
    "\n",
    "function = tree.arrowed.compile(calc_masses, env={\"mass\": mass, \"h_Fill\": h_Fill})\n",
    "\n",
    "h.Reset()\n",
    "startTime = time.time()\n",
    "\n",
    "for entrystart, entryend, dummy in function.run(10000, \"~/storage/data/SingleElectronRun2015D_16Dec2015_v1_0000_*.root\", \"Events\", reportentries=True):\n",
    "    elapsed = time.time() - startTime\n",
    "    print(\"{} events / {:.2f} seconds = {:.1f} kHz\".format(entryend, elapsed, entryend / elapsed * 1e-3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAf3ElE\nQVR4nO3dbZKiSqIGYJiYfYGbOXcXB1zG9F5OCCvj/sjpnGwgLatEgfR5oqPDQrFItcjX/CDraZoq\nAIA1/9r7AACA4xIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIU\nAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAsQQEAyBIUAIAs\nQQEAyBIUAIAsQQEAyBIUAIAsQQEAyPr33gcAQGnqut77EPjaNE2PPExQAGB7D1ZC7OXxMKfrAQDI\nEhQAgCxBAQDIEhQAYF1d133f730UOxMUACBrGIa9D2FnggIAkCUoAMA9fd+3bfux3RC1qa4AbKuu\nC6lc4sUGmqYZxzHcKKMz4vH3qJD3EoDjKCwoxLLMfjy1x98jXQ8AkNU0zertzyEoAABZggIAkCUo\nAABZggIAkCUoAABZhcxgiR5fYBuA1ymscinP49Mj//3qQzkOn1qA9/CdrSQFBgWBAGB3skIxCgwK\nAOzOd7aDezzJGcwIAGQJCgBAlq4HAN5ndenFtm2fecJndudLWhQAeJO2bS9rvtwxt7JzeMKNj/Jl\n+r6v67qu67Ztc4FpdejAMAx1Xa/uUtd13/fpI8OTzLY/Q4sCAG91u92+9fi+76/X6+royBO1JYRS\nhPUnx3G8XC632y09/mEYxnFc3TcXhsLuaYAIj2yaZhzH6/Uafu+TRy4oAPBWX9buoeabPWzWxRB+\nXNaCq/ve2f42ISXESr2u68vlkqaf+2lgaRkswqsRn7Ou6+v1+nxQ0PUAwCGENvlQg14ul9AIH76I\nV1V1uVxCOIiPic3sYffQPh/3jVVyaPBfbn+/tMoPTQuzu2Ybq6rq+34cx+X2qqoul8vq9tyv+Lmp\nLOWVCOB0cqfiUHV1C3Gvqqput9s0TV3Xrd4Oz9A0TfrjbPfb7TbbHh4fuzxeV/DHhaMKt2MB08Oe\npikccHiJ4isQd69+Nx7E54kF7LouvtR3DuDRQ/1OuU7gIJ8AgE92PyjkvrKmdV74MdRzoZpMn2H2\nhLPHzHafVaVpdbuLWJ2HI4lpYFqULr44s6CQFnb2os1e4TuH8Xh1qesBgLdarbSCH48hCB0KdSJu\nDCP74lyDfcc/hk6QULWHIwk/LkcShHunxRDOYRiu12vICjOhnyLmiWqjC2kXOJgx97osX24ASjKr\nPuOkgGEYQiUa5gLsUh0MwxCGK6aTHUKUCTkmPrKu667rwkDFdHvYPZTxer2GoRtx92mawnjJ+OS3\n2y2O5HjmyAsMCgIBwKdp23Ycx/R7eRj2GCJC3/ehSg5DI3dpVwjV/KyGats2DTeh7u+6bjnjIwxp\nDIVK74q7vOq4q+J69MsrEcDp5E7FucGMobW8+nPwXfWdMQrTn4MWl9unaVqOFnybOC6hSSxHGuYO\nbzmYMaqSMQrpAMYvR24+/joU2KIAwJHFNvPUna/4bdter9e6ru9fqSm0tMerEaRTAMLu8cefHfYz\n4pzM9MoHsyaQTX5LeK3iK/zda1utqqeyGupDP83eRwG8Vl3/Sn+cpr/2OhJWveJU/GB/wWEvuPRO\njxT28feotGpVUIBPUNe/YjhIb3MQTsXH9/h7ZHokAJAlKAAAWYICAJAlKAAAWYICAG8SLpKYu+ut\nh8LDBAUA3iRc52B58YCwAsK2FxU4vrBG9nLl67Audl3Xd16QuMT2cg7kI7t/iwsuAfBWs3rxM9sS\nwgWgwoJVl8slLgARrkUdrrEYrja9fH3ishFxvas40XG2e1VVz8eFAlsU6oy9jwuAqvrz6oRVPijk\n+inu91+cInaEynuapmEYQh0fLyg5jmPXdaEgIQfc3z1cezFsCUtC3G63uPvqRTC/7cFLPZ9FeSUC\nlqrqP6u3OYjcqbiqqnQ9gtzG2VWW4yNn29P1lMMznKVqq5IFGlLpqhb3d1++gFN+qYjckzz4yAJb\nFAA4uPi9Py7qmN51vV7TlY3iatFhGeWw7FOVfAuvfn8Rj3XtKYY7xEEGy1cj3PV4KWLDQ9M0m49R\nOHrs+q7ySgQsaVE4uNypONRk6ffmeLv6/S15VjfFB9xut/Rr9Gx9yHSXKvN9/ThiFRwbQkLR4o9h\necnVajqEpHQlyfhSpE8bNt55HR6vLkurVgUF+ASCwsHdDwrTn7FguWUp1otd16W9DLOnXf3xgGYJ\nYFbTp4WtFqtLh6CwmplWM9bq4tSTrgcAjqxpmjjwcLWFvPtTVVXDMNR1fb1e27a93W6zQQmnkx5/\nnOIYNsYf4xDFdMflfMg4zWFmq6UyTY8E4N3atr1er6EiXK3PYnoYhqHv+9BtX/3ZwF6MGAXC5MbZ\n9tXXZzU9zGZJbDYB5MGWh7Mor0TAkq6Hg8udiqukUyDUQemP6QDG9HZ4tvClOTSkp9unE3Y9hDaS\ntIzp65AO2IxlTMdkpH0K6fb0qWa7Lz1eXZZWrQoK8AkEhYN7JCikFf/055S/3PTIdGM61u90QWFK\nxi3OjjZW8EF8fWZTH9Pd01EIuZdu6fHq8n+XcypDeoEqoFR1/Wua/lre5iA2ORWvNrynG++0zJ9F\nuBLz6vbqgaI9s/vj71Fp1aqgAJ9AUDg4p+Lje/w9MusBAMgSFACArAKnR+bWf9IOBgDfVWBQEAgA\nYCu6HgCArAJbFAA4uK1mNqYXHzz1PMkj06IAwPuEBZQvl8vlcgmLLMe74uoPub2W2y+JjddWfrG4\nlnS6zHQQtt8vS1xLerZ7/adtwtODF2Y6i/JKBCy5MuPB5U7F4XqCsyswzq5efGfH1V8UL3gcF5Y8\n/mUZY8HjhRTTK1RWd5eZnpLLLy4fk+4enj93DI9Xl6VVq4ICfIJZUEj/7XhURLlT8bIWTxPA8t6Q\nAKavgkK6Jbc0c2615V0sa/dQitnB59aJXu6eLhvxYEkfry51PQDnNk1/xX97HwtfS5c3rKpqGIZQ\nvYWehXEcw42wqHTsU5jtdUdoro+N9qGJPj5P2Dhrk++6Ljev/nVmizIEs6Eb6ZWqZ1bXlY6P3Gzd\nyODBQPEiy1cqRqG43PiyEWnZVhPtXiLgDXItB1oUDiJ3Kk7bzJet4mmtlJ7k7zTCr9YR1Z8rKC4X\naZy1T6w+yXvE7pLVZoDcXbmei2WVeudXP15d7lytxigQLftpZqVNd1m+goICfAJB4eDunIrjl8Dl\n98NZBZ+e3r8bFKpkZeq4PdSj0zT9888/VVX9/fff8fY///zzo4I+Ja3Xl/feH2+RvoCzBTm7rgsj\nNu5nhdMEhdyrsIx7y3XKlw+bBAX4DILCwT1yKk4TQ9wr7aqfPeEPWhSqNSF/5FoX3i9t6ghiHZ8b\nipi+GsvdZ8+TG7LweKmPMkZh1qEyjuNsre7r9RofFjuf7vTfAHAoszWRw7y+UJltPqdxNgQhFTb+\n/fffYdzDrLp5j9nlH5qmiYMw2ra9Xq+h4l99WcK+IR8sd09tdWGJ/YNCGLoShpmkL8qdNbbvPwaA\nA2rbdhzHWeUXTuOzk/kzXwLDvulv6X8L0SQ8IISDcPvvv//+wS96xuVyWQ0Bfd+P49h13c/Kvrwu\nxc8Ob+7BlodXiIEoTFxJxxxUfza5xKah2THPeiKmTEPTkV8E4Ad0PRxc7ryadqIv+9GXt9PaYfU5\nw/bQVBAfFiuF5WDG6s9O7b2qgFi66c+LSSzH7cWhe7N5pHd2Xw7ezB3Do0f7s0K+SJV0LK0GhVln\n0uqAl3cdLLAbQeHg7pyKl+388a40N8R6PbhzHYXZw2Zd8rOJAOldoSEhDGl8s2XpwvZqTagN01cg\nt/u0eHnvHMPj1WWdO7JdhJms0zTVdd00TWw2CQ1W0zTFG2H7MAyXy+V2u8XGlro+VomAV6jrX6tX\nTcht582+PBXn1nqYbd92SYij9VY/eVR3dn/kmR+vLvesVvu+v16v6QHEfBATw2z7bJfVZxAUoHiC\nwsE5FR/f4+/RnoMZ4xiWsBBIOgIltKukF9hK/4/jXMLQ0LcfOAB8ip1DX+g7iD92XZdedzNMiZxt\nn+0yO34xFj6BFoWDcyo+vnN0PUT3O1q+1QHj0wmfQFA4uPcvncAPnCkobEhQgE8gKHBY5VVD+19w\nCQA4LEEBAMgSFACALEEBAMgSFACArH/vfQDby03LKWwYKgC8QYFBQSAAgK3oegAAsgQFACBLUAAA\nsgQFACBLUAAAsgQFACBLUAAAsgQFACBLUAAAsgQFACBLUAAAsgpc68GiUACwlQKDgkAAAFvR9QAA\nZAkKAECWoAAAZAkKAECWoAAAZBU46wEoUl3/2vsQ4BMJCsBpTNNfex8CfBxdDwBAlqAAAGQJCgBA\nlqAAAGQVOJjRolAAsJUCg4JAAABb0fUAAGQJCgBAlqAAAGQJCgBAlqAAAGQJCgBAlqAAAGQJCgBA\nlqAAAGQVeGVGoCR1/WvvQ4CPJigARzdNf+19CPC5CgwKFoWCTxZbIMQL2ESBQUEggI8Vw4EOC9iK\nwYwAQJagAABkCQoAQJagAABkCQoAQJagAABkCQoAQNaBgkJd18MwxB+HYWjbtq7rtm1nj+z7vq7r\n2eMBgM0d5YJLy8spXi6XqqqaphnHsa7reBmltm3HcWyaJjzmdrstkwQAsIlDtCj0fT/bEur+aZqG\nYQgRITxmGIZxHLuuG4ZhGIamaUKeAABeYf+gMAzD9XoNLQRRbDMImqa5Xq/hwVUSLEKe0AEBAC+y\nf1C4XC5N0+QaFWZmmUCnAwC81M5B4btNAuM4LjfOdq+/78lSAECp9gwKfd+P43i73R7fZdZDEcza\nFabve7IgAFCqPWc9hJaAdDRiuB1qbiMPAGB3ewaFvu9jGgjTGZqmic0DaS9DHNsY5kbG7WF3IxUA\n4FV+0FD/CqED4na7pT92XTdNU9d16V1VVTVNEx8TbkfHKRGwiar6z9v2gueVVw0d5YJLM23bdl13\nvV7DrMiu62Kzwe12u1wucQSiHgoAeJ3/XfHwmMKFnFe3V2udDuk1HIEC1PWvafrrPXvB88qrhg7a\nohDlxh8YlwAAb7D/BZcAgMMSFACALEEBAMgSFACALEEBAMg6+qyHH8gt8lTYfBUAeIMCg4JAAABb\n0fUAAGQJCgBAlqAAAGQJCgBAlqAAAGQJCgBAlqAAAGQJCgBAlqAAAGQJCgBAlqAAAGQVuNaDRaEA\nYCsFBgWBAAC2ousBAMgSFACALEEBAMgSFACALEEBAMgSFACALEEBAMgSFACALEEBAMgSFACALEEB\nAMgqcK0Hi0IBwFYKDAoCAQBsRdcDAJAlKAAAWQV2PQBnV9e/9j4E4L8EBeCIpumvvQ8BqCpdDwDA\nHYICAJAlKAAAWYICAJAlKAAAWYICAJAlKAAAWQVeR8GiUACwlQKDgkAAAFvR9QAAZAkKAECWoAAA\nZAkKAECWoAAAZAkKAEDW/kGh7/u6ruu6bts23T4MQ9u2y+3pLsMwvOswgZOp61/x397HAie283UU\n2rYdx7FpmqqqxnGs6zpeBeFyuVRV1TTNbHu6y+Vyud1uyyQBfLhp+iveFhTgGXu2KAzDMI5j13XD\nMAzD0HVd2FhVVaj7p2kahiFEhL7vl7s0TRPyBADwCvt3PcT2gLRhILYZBE3TXK/X6neMCKEh7qID\nAgBeZM+g0LbtNE1t2w7D0Pd9aBtYzQ3RLBPodACAlzrEWg+x+yD0PtwxjuNyYxj2GH/MLQp1h+Uh\nAGDVIYJCGIswDEPoX4g9C0thbONs46xdQa0PAFvZeTBj7Epo2zYOV4z37nJUAEC0c1C4M2chbTmI\nYxuX11pYbuS8THwHOJqdBzOG/9P6PrQr3G63eHv5f5zscL1e08kRFGCa/kpnwFM86RCObtrVbPRi\n13Wrd6XbQ4bIHf/uJeIZVfWf2Q2Kt/qmb/sB8HHincqrhv53xcMd3elBmM1o+HKX9BqOnE5d/wrN\nCfEGxVt907f9APg48U7lVUOHmPVwZ5BB7i7jEj7HrEXaGR/gnQ4RFOC+9IvmvkcC8Gn2v4QzAHBY\nWhTYmUYCgCMTFNifYQcAhyUosA1DDgGKVGBQyC0KVdh8lQMy5BCgPAUGBYHgCNKsoHUB4LwKDArs\nLk0Gm7cuxCeUPwDeQFDgTPRuALyZ6ygAAFlaFDg9Ey6KoaEIDkhQoAS6JAog4cExCQqUT5MDwI8J\nCnwETQ4APyMocFZfVvkyAcDzBAVO6cHuA70MAE8SFDioN7QHuHwkwJcEBY7obdV2+EU6KQByCgwK\nFoUCgK0UGBQEAgDYSoFBgVdzWQKAzyEo8BPpZQken6YoUgCcjqDAU76s+13pCODUBAXex3REgNMR\nFHiTNBkcsHVBiAFYJSjwkM2r9kNlhYOHGIAdCQo8asPv2b6yA5yFoMDH0WYA8DhBgc/yfGOGy0gA\nH0VQoDRvaDAw5/N0DFaFHxMUKIo6gCWDVeEZBQYFi0IBwFYKDAoCAQBspcCgAFsxbhFAUIB7jFsE\nPpygQJaqEQBBgXs0tgN8uH/tfQAAwHEJCgBAlqAAAGQZowBzuVGcRncCH0hQgD/kxm8a1wl8JkEB\ntmHZIaBIxijAZqbpLxEBKEyBLQoWhQKArRQYFAQCANiKrgcAIEtQAACyBAUAIKvAMQrAwbl0FZyI\noADswDxSOAtdDwBAlqAAAGTtHxT6vq/ruq7rtm2HYYjbh2Fo2zZsz+2SPh4A2NzOYxTath3HsWma\nqqrGcbxcLrfbLSSDy+VSVVXTNOM41nUdL6OU7pI+HgDY3M4tCqHKH4ZhGIYQBfq+r6oq1P3TNM22\nD8MwjmPXdWGXpmlCngAAXmHPoBA6DkICiMZxrH4HiLixaZrr9brcJeQJHRCcTl3/iv/2PhaAe/YM\nCm3bTtMUOw5Cfd91Xbx3ucssE+h04LwsNQmcwlGuo9D3fWgzmDUwzIT2hpkw7DH+mFs98g7rSKV8\nx/2W1Zcr3SgNAKe2f1AYhiGOW/yyEyGMbZxtnLUrqPWfp2570OoLlW6UuoCz23kwY0wJt9ttlhKM\nPACA3e0cFC6XS9M06UiFKG05iGMbZw8LYcJIBQB4kT2DQjp5IQobb7dbfMDy/zjZ4Xq9ppMjgENJ\nJ3fohYGT2nOMQmgPWB2f2LZt13XX6zWMcOy6LjYb3G63y+USRyzqoeAsPrOmjCM2PrP4UID9g0JO\n3/d9389mNFS/J1XqdOCMjBIFTmf/WQ/35aKAiAAAb3D0oABnp8kdODVBAV7oPX0Nru8EvI6gACUI\n+eDgrRcHPzxglaAAvIOmDjipnS+4BAAcWYEtCrlFoawBwRnNmut9LwferMCgIBBQGNcsAnak6wEA\nyBIUAICsArse4FxWr4KglwE4CEEB9pQOTpwtsfizcYt7XXzJoEsolaAAR/F85TqLHU8+249/u+YQ\nKImgAGfyZINB7nv/t9oDXDEaPoqgAKexSYPBl+HgkWc+xRWjgU0ICnBWX9bTvvqvii+L1wQeISh8\nLqPPTu3L92vf8QqHZSAFfJeg8NGcNAG4zwWXAICsAlsULAoFb6MtCopXYFCIgcApDF7KuBb4BAUG\nhZQ+eNiFvzgoRuFBAXg/LQ1QEkGB//IVkBkTaIFKUCBQB7BK5x0gKABzq7FAVoDPJCgAf1htXtLm\nBB9LUIDP9eGNBNbCgEcICvChPrxqtBYGPEhQgI+gLgR+RlCA8n144wHwDItCAQBZBbYopItC1fX/\nJff4UgUA31Ngi8L0W1X9J72993EBwPkUGBQAgK0ICgBAlqAAAGQJCgBAlqAAAGQVOD2S+1yhD4DH\nCQqfyHX6AHiQrgcAIEuLAsD/zPrmNL+BoADwhxgODOiBStcDAHBHgS0KFoUCgK0U2KJgUSgA2EqB\nQQEA2EqBXQ8A32XcIuQICsCnMwcS7hAUAL7g4gp8sgKDgiZEYHMursDHKjAoCPsAsJWjzHqo63oY\nhnTLMAxt29Z13bbt7MF939d1vdyFnLr+Ff/tfSwAnMkhWhT6vl9uvFwuVVU1TTOOY13X0zSF7W3b\njuPYNE14zO12WyYJljS0APADO7cohLaB6/U62x7q/mmahmEIESGEiWEYxnHsum4YhmEYmqYJeQIA\neIWdg0Lbtl3XheaBVGwzCJqmCWEi9DXEFoiQJ3RAAMCL7B8U+r5f7XpY7VCYZQKdDgDwUkcZzPig\ncRyXG2fpoV5TVf+3uv33vQDAikMMZnxcGNs42zhrV4jDHlPpaH/D+gDgQccNCtuOPHCxFAD4geN2\nPaQtB3Fs46zxIIQJIxUA4EUOGhRut1v1e3bD8v842eF6vS5nTAAAWzlo10OYNnm9XsOsyK7rYrPB\n7Xa7XC5xBKK5kQDwOvXq0L/jCBdyXt1erXU6pNdwXFXXvz5wMONnlhqetDoI2l8T931ZDZ1OceU5\nW1B4z/K1Rys1nNe2f7MWsC5PeUHhoF0PH8WMDDiRtC7f5G/WGYCDExR+Ivf37NsA8DqaH9iFoPBD\nyz/Rbb8NPH9G8O0EyqP5gfcTFI7rB2cEXzgA2JagUBrhAIANFRgUcos8FTYMFQDeoMCgcNhAoF8A\ngNMpMCgcSi4cvHTkI3BAb17D1pK5bEVQeLlN/kS/jAJOBHB84e+0rn+9J9zHX/eG30XBBIUtvSjC\nCwFQEn/RnIugsJnNr9cGHJ/WPoonKAD80Jch4PleBt862J2g8Cr+vIFNmhPuP4npVLyaoPBf2w4v\n8LcK/Nh3v2Z8OZ3KDAieISj8jxHCwL7i+efH1fnyDPbj4VPaKgg+MSgI18ABHbAt0xpUVB8YFGbh\n+mif/qMdDwAf7uOCQupozQlHOx4AKDAobLgolC46gEc4WxaswKCw7aJQr+ui08sAPGPDc8gmT2VA\nQ6kKDAqnIG4Dz3jduEWYERQA+NpqO4F48QkEBYAP8rP54auPvNPF8PwFIR5nxvurCQoAn+I988Of\nH6ywOjTyznhJl8t7KUEB4BMd/Mv3atowXnIXgsI3+GgCbOjIkyqPfGxvJig86pM/JQBPytW7aSPB\nI4tavZMGjEBQWPHhnwmAV7hf7+a+jOW2O1G/jaAwp+UA4EG7f9fP0XGwIUEBgJ84eO27bcfBJ0/C\nFBQA+DjfanKYTSt91TEdVYFBYcNFoQB4hSNUt7lxlD9oMPizOP959sgOpsCgIBAAHNn7m+7vXyly\nk+tQJbGjrqqi+iYKDAoAHMSRWw7uP5hIUADgJY5W6R7teM5CUACgEEdowCiPoABACTQYvIigAADf\n8GntFnVhcwTqurQSAXAi5VVD/9r7AACA4xIUAIAsQQEAyBIUAIAsQeFMcstYFKPsAirdeZVduqr0\nApZdujcocHqkRaEAYCsFBgWBAAC2ousBAMgSFACALEEBAMgSFP7rZ8NifzyY9s2jcJVuwx3LLt0z\nO77z1yndtju++deV/adXHkEBAMgSFACALEEBAMgSFACArNKWzTb8BIB9lVaxFlYeAGBDuh4AgCxB\nAQDIEhQAgKxygkLf93Vd13Xdtu3ex7KBtDjDMMTtwzC0bVtMMauqquu6vALmSlFG6XJ/a2cv3eyj\nWN0tUXwRZrsc2f0Czu46XQHvHGrbtrN38HSl29lUhK7rqqpqmqZpmnBj7yN6SixFuFFV1e12C3eF\nH+P2XQ9zA7PSTUUU8Ha75UpRQOmWH85416lLF84h6Udxypdo9iLM9jqmZQGXH9Su68Jdpyvg6tuX\n3pVWCqcr3e5O9seck34Owsdi18N51uxjHX8MH+t0e/zDPqPwTqV/qGUUMH37wrk4lKKY0sXDDqUL\nb995Sxc/h7M6I1ei9D1dPuyAHixg/PFcBcyVLkjDULrlLKU7iBJeoPRsFZzoJLW0Wpzq98TcZS5+\n+wFuIxRzlugLKODy7YsKKN2UDwrnLd3tduu6bvnlMlei2VeRO99lD+LBAsZynauAudIF8SST+yZ5\n8NIdRDljFGbO2/PUtu00TbFHLRQkpuaT9v4uXS6Xpmn6vp9tL6OAub7tAkrXNM31em3btu/7y+VS\nJYU6aelCWZYfxSpTotm55filzhVwmqa0LNfrNdw4VwG/fPuWwzKWj+G+EoJCeOOLfL/juXj1z+C8\nVv+AyxAKFc65TdOM41jY1ULDR3Ecx1DGtOH3Q4zjuNx46g/zMAzhUxqaiMooYN/34zjGroeojNK9\nWQlBISjsnQ5/utfrNbSY7X04W8r9AZek67phGIZhCO9dSSk2NAWFBsmu667Xa2Ep9ktx6F/qvG9x\n27bh28jtdgulKKOA1+u167rlYZdRujcrISisvsenfuOHYYh/uvfbzc4oFOFyudR1HYoZbqf3nlf4\n4M0+fvFLzNlLFzJBLEX4MddkXYDySjRT1/U4jmHQyanPmTPxkxnmQFZVFdr2in9DX+Tfex/ABmI7\ndtqgfeoPffjStvqZTtvNxnFcTccH1/d9LNowDKEU8f06ewFnn8YgluLspVtVcOlWS9S2bbr9vCec\ncMzLBssCCjg72thNFi6ocPbS7WCXIZSbS4e8nr1c6TUhouW8rDIG687mCJRRwPTTmN4uo3Thw5mW\nrowP53K6yp0SVb9H0cfJO+8+3O/LTaeanWriXecq4J3ZRlN+wvlZSre7E1eoM2n6OdcZamb1e9hs\nbk9w3img0fLPu4wCpu9dWooCSjcbXFJM6VZrmlyJZi/Cmw/1Z1ZD+dLqvbsd9MO+FRROV7rdFbXM\n9Oc0Is1atstTQAHvfBqV7lxyJSr+hFN2Acsu3baKCgoAwLZKmPUAALyIoAAAZAkKAECWoAAAZAkK\nAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECW\noAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAAZAkKAECWoAAA\nZP0/mV6f7d3+ivIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "h.Draw()\n",
    "c.Draw()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "The interface I've shown is not how this package will ultimately look.\n",
    "\n",
    "   * This `run` function is powerful, but side-effect changing; it encourages non-parallelizable code. I'll leave `run` as a low-level interface, but build the standard `map`, `filter`, `reduce`, etc. on top of it.\n",
    "   * I should probably generalize and wrap up this method of exposing ROOT functions to Arrowed.\n",
    "   * Histogrammar, a small suite of physics functions, ...\n",
    "\n",
    "Interested in trying it out? Giving me feedback?"
   ]
  }
 ],
 "metadata": {
  "celltoolbar": "Slideshow",
  "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.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}