reprocessing_compare.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import uproot"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Samples\n",
    "These ntuples were created using DHNLAlgorithm using all 10K events of\n",
    "new: mc16_13TeV.311625.Pythia8EvtGen_A14NNPDF23LO_WmuHNL50_5G_lt100dd_el.recon.DAOD_RPVLL.e7660_e5984_a875_r11915\n",
    "old: mc16_13TeV.311625.Pythia8EvtGen_A14NNPDF23LO_WmuHNL50_5G_lt100dd_el.merge.DAOD_RPVLL.e7660_e5984_a875_r10740_r10706\n",
    "\n",
    "#### Sample Updated\n",
    "I'm now trying this with the same sample but updated MC campaign. mc16e. (previously mc16a)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "old_tree = uproot.open(\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_old/data-tree/files_old.root\")['outTree']\n",
    "new_tree = uproot.open(\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new/data-tree/files_new.root\")['outTree']\n",
    "# old_tree = uproot.open(\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_old_mc16e/data-tree/files_old.root\")['outTree']\n",
    "# new_tree = uproot.open(\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new_mc16e/data-tree/files_new.root\")['outTree']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "event_columns = [\n",
    "'runNumber',\n",
    "'eventNumber',\n",
    "'lumiBlock',\n",
    "'coreFlags',\n",
    "'bcid',\n",
    "'mcEventNumber',\n",
    "'mcChannelNumber',\n",
    "'mcEventWeight',\n",
    "'NPV',\n",
    "'actualInteractionsPerCrossing',\n",
    "'averageInteractionsPerCrossing',\n",
    "'passedTriggers',\n",
    "'weight_pileup',\n",
    "'correctedAverageMu',\n",
    "'correctedAndScaledAverageMu',\n",
    "'correctedActualMu',\n",
    "'correctedAndScaledActualMu',\n",
    "'rand_run_nr',\n",
    "'rand_lumiblock_nr',]\n",
    "\n",
    "old_df = old_tree.pandas.df(event_columns).set_index(['lumiBlock'])\n",
    "new_df = new_tree.pandas.df(event_columns).set_index(['lumiBlock'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>runNumber</th>\n",
       "      <th>eventNumber</th>\n",
       "      <th>coreFlags</th>\n",
       "      <th>bcid</th>\n",
       "      <th>mcEventNumber</th>\n",
       "      <th>mcChannelNumber</th>\n",
       "      <th>mcEventWeight</th>\n",
       "      <th>NPV</th>\n",
       "      <th>actualInteractionsPerCrossing</th>\n",
       "      <th>averageInteractionsPerCrossing</th>\n",
       "      <th>passedTriggers</th>\n",
       "      <th>weight_pileup</th>\n",
       "      <th>correctedAverageMu</th>\n",
       "      <th>correctedAndScaledAverageMu</th>\n",
       "      <th>correctedActualMu</th>\n",
       "      <th>correctedAndScaledActualMu</th>\n",
       "      <th>rand_run_nr</th>\n",
       "      <th>rand_lumiblock_nr</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lumiBlock</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>5006</td>\n",
       "      <td>0</td>\n",
       "      <td>3120</td>\n",
       "      <td>5006</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>6</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279169</td>\n",
       "      <td>642</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>3000</td>\n",
       "      <td>0</td>\n",
       "      <td>2103</td>\n",
       "      <td>3000</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279259</td>\n",
       "      <td>470</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>4998</td>\n",
       "      <td>0</td>\n",
       "      <td>555</td>\n",
       "      <td>4998</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>4</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279598</td>\n",
       "      <td>853</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>3008</td>\n",
       "      <td>0</td>\n",
       "      <td>685</td>\n",
       "      <td>3008</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>2</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279984</td>\n",
       "      <td>287</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>8003</td>\n",
       "      <td>0</td>\n",
       "      <td>907</td>\n",
       "      <td>8003</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>280464</td>\n",
       "      <td>611</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>6002</td>\n",
       "      <td>0</td>\n",
       "      <td>2550</td>\n",
       "      <td>6002</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>28</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>305380</td>\n",
       "      <td>350</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>7008</td>\n",
       "      <td>0</td>\n",
       "      <td>2942</td>\n",
       "      <td>7008</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>34</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>311321</td>\n",
       "      <td>191</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>5997</td>\n",
       "      <td>0</td>\n",
       "      <td>1917</td>\n",
       "      <td>5997</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>22</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>305920</td>\n",
       "      <td>227</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>2340</td>\n",
       "      <td>1</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>21</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>302393</td>\n",
       "      <td>292</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>4005</td>\n",
       "      <td>0</td>\n",
       "      <td>2768</td>\n",
       "      <td>4005</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>23</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>305723</td>\n",
       "      <td>374</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>10000 rows × 18 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           runNumber  eventNumber  coreFlags  bcid  mcEventNumber  \\\n",
       "lumiBlock                                                           \n",
       "5             284500         5006          0  3120           5006   \n",
       "5             284500         3000          0  2103           3000   \n",
       "5             284500         4998          0   555           4998   \n",
       "5             284500         3008          0   685           3008   \n",
       "5             284500         8003          0   907           8003   \n",
       "...              ...          ...        ...   ...            ...   \n",
       "45            284500         6002          0  2550           6002   \n",
       "45            284500         7008          0  2942           7008   \n",
       "45            284500         5997          0  1917           5997   \n",
       "45            284500            1          0  2340              1   \n",
       "45            284500         4005          0  2768           4005   \n",
       "\n",
       "           mcChannelNumber  mcEventWeight  NPV  actualInteractionsPerCrossing  \\\n",
       "lumiBlock                                                                       \n",
       "5                   311625            1.0    6                            4.5   \n",
       "5                   311625            1.0    5                            4.5   \n",
       "5                   311625            1.0    4                            4.5   \n",
       "5                   311625            1.0    2                            4.5   \n",
       "5                   311625            1.0    5                            4.5   \n",
       "...                    ...            ...  ...                            ...   \n",
       "45                  311625            1.0   28                           44.5   \n",
       "45                  311625            1.0   34                           44.5   \n",
       "45                  311625            1.0   22                           44.5   \n",
       "45                  311625            1.0   21                           44.5   \n",
       "45                  311625            1.0   23                           44.5   \n",
       "\n",
       "           averageInteractionsPerCrossing  \\\n",
       "lumiBlock                                   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "...                                   ...   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "\n",
       "                                              passedTriggers  weight_pileup  \\\n",
       "lumiBlock                                                                     \n",
       "5                                                         []       0.032687   \n",
       "5                                                         []       0.032687   \n",
       "5          [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.032687   \n",
       "5          [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.032687   \n",
       "5                                                         []       0.032687   \n",
       "...                                                      ...            ...   \n",
       "45                                                        []       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "45         [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "\n",
       "           correctedAverageMu  correctedAndScaledAverageMu  correctedActualMu  \\\n",
       "lumiBlock                                                                       \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "...                       ...                          ...                ...   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "\n",
       "           correctedAndScaledActualMu  rand_run_nr  rand_lumiblock_nr  \n",
       "lumiBlock                                                              \n",
       "5                                 4.5       279169                642  \n",
       "5                                 4.5       279259                470  \n",
       "5                                 4.5       279598                853  \n",
       "5                                 4.5       279984                287  \n",
       "5                                 4.5       280464                611  \n",
       "...                               ...          ...                ...  \n",
       "45                               44.5       305380                350  \n",
       "45                               44.5       311321                191  \n",
       "45                               44.5       305920                227  \n",
       "45                               44.5       302393                292  \n",
       "45                               44.5       305723                374  \n",
       "\n",
       "[10000 rows x 18 columns]"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "old_df.sort_index(inplace=True)\n",
    "old_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>runNumber</th>\n",
       "      <th>eventNumber</th>\n",
       "      <th>coreFlags</th>\n",
       "      <th>bcid</th>\n",
       "      <th>mcEventNumber</th>\n",
       "      <th>mcChannelNumber</th>\n",
       "      <th>mcEventWeight</th>\n",
       "      <th>NPV</th>\n",
       "      <th>actualInteractionsPerCrossing</th>\n",
       "      <th>averageInteractionsPerCrossing</th>\n",
       "      <th>passedTriggers</th>\n",
       "      <th>weight_pileup</th>\n",
       "      <th>correctedAverageMu</th>\n",
       "      <th>correctedAndScaledAverageMu</th>\n",
       "      <th>correctedActualMu</th>\n",
       "      <th>correctedAndScaledActualMu</th>\n",
       "      <th>rand_run_nr</th>\n",
       "      <th>rand_lumiblock_nr</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>lumiBlock</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>4013</td>\n",
       "      <td>0</td>\n",
       "      <td>1467</td>\n",
       "      <td>4013</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>3</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>280319</td>\n",
       "      <td>714</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>3052</td>\n",
       "      <td>0</td>\n",
       "      <td>2047</td>\n",
       "      <td>3052</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>2</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279598</td>\n",
       "      <td>621</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>49</td>\n",
       "      <td>0</td>\n",
       "      <td>2690</td>\n",
       "      <td>49</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>3</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>279598</td>\n",
       "      <td>793</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>7916</td>\n",
       "      <td>0</td>\n",
       "      <td>2853</td>\n",
       "      <td>7916</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>6</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>280464</td>\n",
       "      <td>550</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>284500</td>\n",
       "      <td>6985</td>\n",
       "      <td>0</td>\n",
       "      <td>3028</td>\n",
       "      <td>6985</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>2</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.032687</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>4.5</td>\n",
       "      <td>280319</td>\n",
       "      <td>637</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>8997</td>\n",
       "      <td>0</td>\n",
       "      <td>3261</td>\n",
       "      <td>8997</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>24</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>304198</td>\n",
       "      <td>285</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>1</td>\n",
       "      <td>0</td>\n",
       "      <td>2340</td>\n",
       "      <td>1</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>21</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>302393</td>\n",
       "      <td>292</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>4005</td>\n",
       "      <td>0</td>\n",
       "      <td>2550</td>\n",
       "      <td>4005</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>24</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>305723</td>\n",
       "      <td>374</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>1999</td>\n",
       "      <td>0</td>\n",
       "      <td>3290</td>\n",
       "      <td>1999</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>28</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>279284</td>\n",
       "      <td>852</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>284500</td>\n",
       "      <td>7008</td>\n",
       "      <td>0</td>\n",
       "      <td>2284</td>\n",
       "      <td>7008</td>\n",
       "      <td>311625</td>\n",
       "      <td>1.0</td>\n",
       "      <td>25</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>[]</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>44.5</td>\n",
       "      <td>311321</td>\n",
       "      <td>191</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>10000 rows × 18 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           runNumber  eventNumber  coreFlags  bcid  mcEventNumber  \\\n",
       "lumiBlock                                                           \n",
       "5             284500         4013          0  1467           4013   \n",
       "5             284500         3052          0  2047           3052   \n",
       "5             284500           49          0  2690             49   \n",
       "5             284500         7916          0  2853           7916   \n",
       "5             284500         6985          0  3028           6985   \n",
       "...              ...          ...        ...   ...            ...   \n",
       "45            284500         8997          0  3261           8997   \n",
       "45            284500            1          0  2340              1   \n",
       "45            284500         4005          0  2550           4005   \n",
       "45            284500         1999          0  3290           1999   \n",
       "45            284500         7008          0  2284           7008   \n",
       "\n",
       "           mcChannelNumber  mcEventWeight  NPV  actualInteractionsPerCrossing  \\\n",
       "lumiBlock                                                                       \n",
       "5                   311625            1.0    3                            4.5   \n",
       "5                   311625            1.0    2                            4.5   \n",
       "5                   311625            1.0    3                            4.5   \n",
       "5                   311625            1.0    6                            4.5   \n",
       "5                   311625            1.0    2                            4.5   \n",
       "...                    ...            ...  ...                            ...   \n",
       "45                  311625            1.0   24                           44.5   \n",
       "45                  311625            1.0   21                           44.5   \n",
       "45                  311625            1.0   24                           44.5   \n",
       "45                  311625            1.0   28                           44.5   \n",
       "45                  311625            1.0   25                           44.5   \n",
       "\n",
       "           averageInteractionsPerCrossing  \\\n",
       "lumiBlock                                   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "5                                     4.5   \n",
       "...                                   ...   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "45                                   44.5   \n",
       "\n",
       "                                              passedTriggers  weight_pileup  \\\n",
       "lumiBlock                                                                     \n",
       "5                                                         []       0.032687   \n",
       "5                                                         []       0.032687   \n",
       "5          [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.032687   \n",
       "5                                                         []       0.032687   \n",
       "5          [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.032687   \n",
       "...                                                      ...            ...   \n",
       "45                                                        []       0.000000   \n",
       "45         [b'HLT_mu24_ivarloose', b'HLT_mu24_iloose', b'...       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "45                                                        []       0.000000   \n",
       "\n",
       "           correctedAverageMu  correctedAndScaledAverageMu  correctedActualMu  \\\n",
       "lumiBlock                                                                       \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "5                         4.5                          4.5                4.5   \n",
       "...                       ...                          ...                ...   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "45                       44.5                         44.5               44.5   \n",
       "\n",
       "           correctedAndScaledActualMu  rand_run_nr  rand_lumiblock_nr  \n",
       "lumiBlock                                                              \n",
       "5                                 4.5       280319                714  \n",
       "5                                 4.5       279598                621  \n",
       "5                                 4.5       279598                793  \n",
       "5                                 4.5       280464                550  \n",
       "5                                 4.5       280319                637  \n",
       "...                               ...          ...                ...  \n",
       "45                               44.5       304198                285  \n",
       "45                               44.5       302393                292  \n",
       "45                               44.5       305723                374  \n",
       "45                               44.5       279284                852  \n",
       "45                               44.5       311321                191  \n",
       "\n",
       "[10000 rows x 18 columns]"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "new_df.sort_index(inplace=True)\n",
    "new_df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Event number mismatch\n",
    "It appears that the events used in the reprocessing are different at least for some lumiblocks.\n",
    "In lumiblock 5-24 the mc event numbers do not match, but in lumiblock 25-45, there is a perfect match.\n",
    "In each lumiblock, the number of events matches."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "lumiblock 5 has 10 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 6 has 10 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 7 has 10 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 8 has 10 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 9 has 10 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 10 has 30 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 11 has 110 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 12 has 170 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 13 has 250 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 14 has 300 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 15 has 340 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 16 has 350 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 17 has 360 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 18 has 400 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 19 has 460 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 20 has 470 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 21 has 480 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 22 has 490 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 23 has 460 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 24 has 490 events\n",
      "event numbers match: False\n",
      "\n",
      "lumiblock 25 has 480 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 26 has 550 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 27 has 450 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 28 has 440 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 29 has 440 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 30 has 400 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 31 has 380 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 32 has 340 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 33 has 300 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 34 has 270 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 35 has 200 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 36 has 130 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 37 has 110 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 38 has 80 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 39 has 70 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 40 has 60 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 41 has 40 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 42 has 20 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 43 has 10 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 44 has 10 events\n",
      "event numbers match: True\n",
      "\n",
      "lumiblock 45 has 10 events\n",
      "event numbers match: True\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for i in set(new_df.index):\n",
    "    print(\"lumiblock\", i, \"has\", len(new_df.loc[i]), \"events\")\n",
    "    if not len(new_df.loc[i]) == len(old_df.loc[i]): print(\"event count mismatch\")\n",
    "    print(\"event numbers match:\", set(new_df.loc[i].eventNumber) == set(old_df.loc[i].eventNumber))\n",
    "    print()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### No overlap in mismatched blocks\n",
    "There is no overlap in event number for the unmatched lumiblocks.\n",
    "This pattern is the same for all mismatched blocks."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "lumiBlock\n",
      "5      49\n",
      "5    1060\n",
      "5    2075\n",
      "5    3052\n",
      "5    4013\n",
      "5    5011\n",
      "5    6985\n",
      "5    7916\n",
      "5    8065\n",
      "5    9055\n",
      "Name: eventNumber, dtype: int64 lumiBlock\n",
      "5     999\n",
      "5    1002\n",
      "5    3000\n",
      "5    3008\n",
      "5    4998\n",
      "5    5006\n",
      "5    6997\n",
      "5    7993\n",
      "5    8003\n",
      "5    9007\n",
      "Name: eventNumber, dtype: int64\n"
     ]
    }
   ],
   "source": [
    "# Mismatched block 5\n",
    "print(new_df.loc[5].eventNumber.sort_values() , old_df.loc[5].eventNumber.sort_values())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "lumiBlock\n",
      "45       1\n",
      "45    1999\n",
      "45    2003\n",
      "45    3996\n",
      "45    4005\n",
      "45    5997\n",
      "45    6002\n",
      "45    7008\n",
      "45    8997\n",
      "45    9995\n",
      "Name: eventNumber, dtype: int64 lumiBlock\n",
      "45       1\n",
      "45    1999\n",
      "45    2003\n",
      "45    3996\n",
      "45    4005\n",
      "45    5997\n",
      "45    6002\n",
      "45    7008\n",
      "45    8997\n",
      "45    9995\n",
      "Name: eventNumber, dtype: int64\n"
     ]
    }
   ],
   "source": [
    "# Matching block 45\n",
    "print(new_df.loc[45].eventNumber.sort_values() , old_df.loc[45].eventNumber.sort_values())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [],
   "source": [
    "old_df = old_tree.pandas.df(flatten=False).set_index(['eventNumber']).sort_index()\n",
    "new_df = new_tree.pandas.df(flatten=False).set_index(['eventNumber']).sort_index()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "el_pt    [10.549664, 3.395912, 2.9742491]\n",
      "Name: 5997, dtype: object\n",
      "\n",
      "el_pt    [12.493513, 5.387406, 4.7612166, 3.897623, 3.3...\n",
      "Name: 5997, dtype: object\n"
     ]
    }
   ],
   "source": [
    "idx = 5997\n",
    "print(old_df[['el_pt']].loc[idx])\n",
    "print()\n",
    "print(new_df[['el_pt']].loc[idx])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "el_LHLoose     [0, 0, 0]\n",
      "el_LHMedium    [0, 0, 0]\n",
      "el_LHTight     [0, 0, 0]\n",
      "Name: 5997, dtype: object\n",
      "\n",
      "el_LHLoose     [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
      "el_LHMedium    [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
      "el_LHTight     [0, 0, 0, 0, 0, 0, 0, 0, 0, 0]\n",
      "Name: 5997, dtype: object\n"
     ]
    }
   ],
   "source": [
    "print(old_df[[ 'el_LHLoose', 'el_LHMedium', 'el_LHTight',]].loc[idx])\n",
    "print()\n",
    "print(new_df[[ 'el_LHLoose', 'el_LHMedium', 'el_LHTight',]].loc[idx])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "el_LHLoose\n",
      "old: 620\n",
      "new: 585\n",
      "\n",
      "el_LHMedium\n",
      "old: 185\n",
      "new: 183\n",
      "\n",
      "el_LHTight\n",
      "old: 72\n",
      "new: 82\n"
     ]
    }
   ],
   "source": [
    "print('el_LHLoose')\n",
    "print('old:', old_df.el_LHLoose.explode().sum())\n",
    "print('new:', new_df.el_LHLoose.explode().sum())\n",
    "print()\n",
    "print('el_LHMedium')\n",
    "print('old:', old_df.el_LHMedium.explode().sum())\n",
    "print('new:', new_df.el_LHMedium.explode().sum())\n",
    "print()\n",
    "print('el_LHTight')\n",
    "print('old:', old_df.el_LHTight.explode().sum())\n",
    "print('new:', new_df.el_LHTight.explode().sum())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['runNumber',\n",
       " 'lumiBlock',\n",
       " 'coreFlags',\n",
       " 'bcid',\n",
       " 'mcEventNumber',\n",
       " 'mcChannelNumber',\n",
       " 'mcEventWeight',\n",
       " 'NPV',\n",
       " 'actualInteractionsPerCrossing',\n",
       " 'averageInteractionsPerCrossing',\n",
       " 'weight_pileup',\n",
       " 'correctedAverageMu',\n",
       " 'correctedAndScaledAverageMu',\n",
       " 'correctedActualMu',\n",
       " 'correctedAndScaledActualMu',\n",
       " 'rand_run_nr',\n",
       " 'rand_lumiblock_nr',\n",
       " 'pdgId1',\n",
       " 'pdgId2',\n",
       " 'pdfId1',\n",
       " 'pdfId2',\n",
       " 'x1',\n",
       " 'x2',\n",
       " 'q',\n",
       " 'xf1',\n",
       " 'xf2',\n",
       " 'weight',\n",
       " 'passesHnlMuMuFilter',\n",
       " 'passesHnlElMuFilter',\n",
       " 'passesHnlElElFilter',\n",
       " 'passesHnlMuElFilter',\n",
       " 'passesVH4bFilter',\n",
       " 'passL1',\n",
       " 'passHLT',\n",
       " 'passedTriggers',\n",
       " 'disabledTriggers',\n",
       " 'nmuon',\n",
       " 'muon_m',\n",
       " 'muon_pt',\n",
       " 'muon_phi',\n",
       " 'muon_eta',\n",
       " 'muon_charge',\n",
       " 'muon_isTrigMatched',\n",
       " 'muon_isTrigMatchedToChain',\n",
       " 'muon_listTrigChains',\n",
       " 'muon_isTight',\n",
       " 'muon_isMedium',\n",
       " 'muon_isLoose',\n",
       " 'muon_trkd0',\n",
       " 'muon_trkd0sig',\n",
       " 'muon_trkz0',\n",
       " 'muon_trkz0sintheta',\n",
       " 'muon_trkphi0',\n",
       " 'muon_trktheta',\n",
       " 'muon_trkcharge',\n",
       " 'muon_trkqOverP',\n",
       " 'muon_EnergyLoss',\n",
       " 'muon_EnergyLossSigma',\n",
       " 'muon_energyLossType',\n",
       " 'muon_MeasEnergyLoss',\n",
       " 'muon_MeasEnergyLossSigma',\n",
       " 'muon_ParamEnergyLoss',\n",
       " 'muon_ParamEnergyLossSigmaMinus',\n",
       " 'muon_ParamEnergyLossSigmaPlus',\n",
       " 'muon_index',\n",
       " 'muon_type',\n",
       " 'muon_passesPromptCuts',\n",
       " 'muon_passesDisplacedCuts',\n",
       " 'muon_px',\n",
       " 'muon_py',\n",
       " 'muon_pz',\n",
       " 'muon_chi2',\n",
       " 'muon_msInnerMatchDOF',\n",
       " 'muon_isLRT',\n",
       " 'muon_myptcone30',\n",
       " 'muon_myptcone30noLRT',\n",
       " 'nel',\n",
       " 'el_m',\n",
       " 'el_pt',\n",
       " 'el_phi',\n",
       " 'el_eta',\n",
       " 'el_caloCluster_eta',\n",
       " 'el_charge',\n",
       " 'el_isTrigMatched',\n",
       " 'el_isTrigMatchedToChain',\n",
       " 'el_listTrigChains',\n",
       " 'el_Loose',\n",
       " 'el_Medium',\n",
       " 'el_Tight',\n",
       " 'el_LHLoose',\n",
       " 'el_LHMedium',\n",
       " 'el_LHTight',\n",
       " 'el_MultiLepton',\n",
       " 'el_trkd0',\n",
       " 'el_trkd0sig',\n",
       " 'el_trkz0',\n",
       " 'el_trkz0sintheta',\n",
       " 'el_trkphi0',\n",
       " 'el_trktheta',\n",
       " 'el_trkcharge',\n",
       " 'el_trkqOverP',\n",
       " 'el_index',\n",
       " 'el_passesPromptCuts',\n",
       " 'el_passesDisplacedCuts',\n",
       " 'el_px',\n",
       " 'el_py',\n",
       " 'el_pz',\n",
       " 'vertex_x',\n",
       " 'vertex_y',\n",
       " 'vertex_z',\n",
       " 'nsecVtx_VSI',\n",
       " 'secVtx_VSI_ID',\n",
       " 'secVtx_VSI_x',\n",
       " 'secVtx_VSI_y',\n",
       " 'secVtx_VSI_z',\n",
       " 'secVtx_VSI_r',\n",
       " 'secVtx_VSI_distFromPV',\n",
       " 'secVtx_VSI_pt',\n",
       " 'secVtx_VSI_eta',\n",
       " 'secVtx_VSI_phi',\n",
       " 'secVtx_VSI_mass',\n",
       " 'secVtx_VSI_mass_nonAssoc',\n",
       " 'secVtx_VSI_covariance',\n",
       " 'secVtx_VSI_chi2',\n",
       " 'secVtx_VSI_direction',\n",
       " 'secVtx_VSI_charge',\n",
       " 'secVtx_VSI_H',\n",
       " 'secVtx_VSI_Ht',\n",
       " 'secVtx_VSI_minOpAng',\n",
       " 'secVtx_VSI_maxOpAng',\n",
       " 'secVtx_VSI_maxDR_wrtSV',\n",
       " 'secVtx_VSI_maxDR',\n",
       " 'secVtx_VSI_mind0',\n",
       " 'secVtx_VSI_maxd0',\n",
       " 'secVtx_VSI_minOneTrackRemovedMass',\n",
       " 'secVtx_VSI_twoTracksMass',\n",
       " 'secVtx_VSI_twoTracksMassRest',\n",
       " 'secVtx_VSI_twoTracksCharge',\n",
       " 'secVtx_VSI_ntrk',\n",
       " 'secVtx_VSI_ntrk_lrt',\n",
       " 'secVtx_VSI_ntrk_sel',\n",
       " 'secVtx_VSI_ntrk_assoc',\n",
       " 'secVtx_VSI_pass_mat',\n",
       " 'secVtx_VSI_trk_ID',\n",
       " 'secVtx_VSI_trk_qOverP',\n",
       " 'secVtx_VSI_trk_theta',\n",
       " 'secVtx_VSI_trk_vz',\n",
       " 'secVtx_VSI_trk_E',\n",
       " 'secVtx_VSI_trk_M',\n",
       " 'secVtx_VSI_trk_pt',\n",
       " 'secVtx_VSI_trk_px',\n",
       " 'secVtx_VSI_trk_py',\n",
       " 'secVtx_VSI_trk_pz',\n",
       " 'secVtx_VSI_trk_eta',\n",
       " 'secVtx_VSI_trk_phi',\n",
       " 'secVtx_VSI_trk_d0',\n",
       " 'secVtx_VSI_trk_z0',\n",
       " 'secVtx_VSI_trk_errd0',\n",
       " 'secVtx_VSI_trk_errz0',\n",
       " 'secVtx_VSI_trk_chi2',\n",
       " 'secVtx_VSI_trk_nDoF',\n",
       " 'secVtx_VSI_trk_pt_wrtSV',\n",
       " 'secVtx_VSI_trk_eta_wrtSV',\n",
       " 'secVtx_VSI_trk_phi_wrtSV',\n",
       " 'secVtx_VSI_trk_d0_wrtSV',\n",
       " 'secVtx_VSI_trk_z0_wrtSV',\n",
       " 'secVtx_VSI_trk_errd0_wrtSV',\n",
       " 'secVtx_VSI_trk_errz0_wrtSV',\n",
       " 'secVtx_VSI_trk_chi2_toSV',\n",
       " 'secVtx_VSI_trk_charge',\n",
       " 'secVtx_VSI_trk_isFinal',\n",
       " 'secVtx_VSI_trk_isSelected',\n",
       " 'secVtx_VSI_trk_isAssociated',\n",
       " 'secVtx_VSI_trk_isLRT',\n",
       " 'secVtx_VSI_trk_electronIndex',\n",
       " 'secVtx_VSI_trk_muonIndex',\n",
       " 'secVtx_VSI_trk_nPixelHits',\n",
       " 'secVtx_VSI_trk_nSCTHits',\n",
       " 'secVtx_VSI_trk_nTRTHits',\n",
       " 'secVtx_VSI_trk_nPixelBarrelLayers',\n",
       " 'secVtx_VSI_trk_nPixelEndCapLayers',\n",
       " 'secVtx_VSI_trk_nSCTBarrelLayers',\n",
       " 'secVtx_VSI_trk_nSCTEndCapLayers',\n",
       " 'secVtx_VSI_trk_hitPattern',\n",
       " 'secVtx_VSI_trk_dEdx',\n",
       " 'secVtx_VSI_trk_radiusOfFirstHit',\n",
       " 'secVtx_VSI_trk_truthProb',\n",
       " 'secVtx_VSI_trk_truthID',\n",
       " 'secVtx_VSI_trk_truthBarcode',\n",
       " 'secVtx_VSI_trk_truthPdgId',\n",
       " 'secVtx_VSI_trk_truthPointsToClosestTV',\n",
       " 'secVtx_VSI_trk_truthPointsToMaxlinkTV',\n",
       " 'secVtx_VSI_trk_truthParentPointsToMaxlinkPTV',\n",
       " 'secVtx_VSI_closeTruth_distance',\n",
       " 'secVtx_VSI_closeTruth_ID',\n",
       " 'secVtx_VSI_closeTruth_barcode',\n",
       " 'secVtx_VSI_closestTruth_distance',\n",
       " 'secVtx_VSI_closestTruth_ID',\n",
       " 'secVtx_VSI_closestTruth_x',\n",
       " 'secVtx_VSI_closestTruth_y',\n",
       " 'secVtx_VSI_closestTruth_z',\n",
       " 'secVtx_VSI_closestTruth_r',\n",
       " 'secVtx_VSI_closestTruth_eta',\n",
       " 'secVtx_VSI_closestTruth_phi',\n",
       " 'secVtx_VSI_closestTruth_barcode',\n",
       " 'secVtx_VSI_linkTruth_score',\n",
       " 'secVtx_VSI_linkTruth_ID',\n",
       " 'secVtx_VSI_linkTruth_barcode',\n",
       " 'secVtx_VSI_maxlinkTruth_score',\n",
       " 'secVtx_VSI_maxlinkTruth_ID',\n",
       " 'secVtx_VSI_maxlinkTruth_x',\n",
       " 'secVtx_VSI_maxlinkTruth_y',\n",
       " 'secVtx_VSI_maxlinkTruth_z',\n",
       " 'secVtx_VSI_maxlinkTruth_r',\n",
       " 'secVtx_VSI_maxlinkTruth_eta',\n",
       " 'secVtx_VSI_maxlinkTruth_phi',\n",
       " 'secVtx_VSI_maxlinkTruth_barcode',\n",
       " 'secVtx_VSI_maxlinkTruth_parent_pdgId',\n",
       " 'nsecVtx_VSI_Leptons',\n",
       " 'secVtx_VSI_Leptons_ID',\n",
       " 'secVtx_VSI_Leptons_x',\n",
       " 'secVtx_VSI_Leptons_y',\n",
       " 'secVtx_VSI_Leptons_z',\n",
       " 'secVtx_VSI_Leptons_r',\n",
       " 'secVtx_VSI_Leptons_distFromPV',\n",
       " 'secVtx_VSI_Leptons_pt',\n",
       " 'secVtx_VSI_Leptons_eta',\n",
       " 'secVtx_VSI_Leptons_phi',\n",
       " 'secVtx_VSI_Leptons_mass',\n",
       " 'secVtx_VSI_Leptons_mass_nonAssoc',\n",
       " 'secVtx_VSI_Leptons_covariance',\n",
       " 'secVtx_VSI_Leptons_chi2',\n",
       " 'secVtx_VSI_Leptons_direction',\n",
       " 'secVtx_VSI_Leptons_charge',\n",
       " 'secVtx_VSI_Leptons_H',\n",
       " 'secVtx_VSI_Leptons_Ht',\n",
       " 'secVtx_VSI_Leptons_minOpAng',\n",
       " 'secVtx_VSI_Leptons_maxOpAng',\n",
       " 'secVtx_VSI_Leptons_maxDR_wrtSV',\n",
       " 'secVtx_VSI_Leptons_maxDR',\n",
       " 'secVtx_VSI_Leptons_mind0',\n",
       " 'secVtx_VSI_Leptons_maxd0',\n",
       " 'secVtx_VSI_Leptons_minOneTrackRemovedMass',\n",
       " 'secVtx_VSI_Leptons_twoTracksMass',\n",
       " 'secVtx_VSI_Leptons_twoTracksMassRest',\n",
       " 'secVtx_VSI_Leptons_twoTracksCharge',\n",
       " 'secVtx_VSI_Leptons_ntrk',\n",
       " 'secVtx_VSI_Leptons_ntrk_lrt',\n",
       " 'secVtx_VSI_Leptons_ntrk_sel',\n",
       " 'secVtx_VSI_Leptons_ntrk_assoc',\n",
       " 'secVtx_VSI_Leptons_pass_mat',\n",
       " 'secVtx_VSI_Leptons_trk_ID',\n",
       " 'secVtx_VSI_Leptons_trk_qOverP',\n",
       " 'secVtx_VSI_Leptons_trk_theta',\n",
       " 'secVtx_VSI_Leptons_trk_vz',\n",
       " 'secVtx_VSI_Leptons_trk_E',\n",
       " 'secVtx_VSI_Leptons_trk_M',\n",
       " 'secVtx_VSI_Leptons_trk_pt',\n",
       " 'secVtx_VSI_Leptons_trk_px',\n",
       " 'secVtx_VSI_Leptons_trk_py',\n",
       " 'secVtx_VSI_Leptons_trk_pz',\n",
       " 'secVtx_VSI_Leptons_trk_eta',\n",
       " 'secVtx_VSI_Leptons_trk_phi',\n",
       " 'secVtx_VSI_Leptons_trk_d0',\n",
       " 'secVtx_VSI_Leptons_trk_z0',\n",
       " 'secVtx_VSI_Leptons_trk_errd0',\n",
       " 'secVtx_VSI_Leptons_trk_errz0',\n",
       " 'secVtx_VSI_Leptons_trk_chi2',\n",
       " 'secVtx_VSI_Leptons_trk_nDoF',\n",
       " 'secVtx_VSI_Leptons_trk_pt_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_eta_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_phi_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_d0_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_z0_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_errd0_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_errz0_wrtSV',\n",
       " 'secVtx_VSI_Leptons_trk_chi2_toSV',\n",
       " 'secVtx_VSI_Leptons_trk_charge',\n",
       " 'secVtx_VSI_Leptons_trk_isFinal',\n",
       " 'secVtx_VSI_Leptons_trk_isSelected',\n",
       " 'secVtx_VSI_Leptons_trk_isAssociated',\n",
       " 'secVtx_VSI_Leptons_trk_isLRT',\n",
       " 'secVtx_VSI_Leptons_trk_electronIndex',\n",
       " 'secVtx_VSI_Leptons_trk_muonIndex',\n",
       " 'secVtx_VSI_Leptons_trk_nPixelHits',\n",
       " 'secVtx_VSI_Leptons_trk_nSCTHits',\n",
       " 'secVtx_VSI_Leptons_trk_nTRTHits',\n",
       " 'secVtx_VSI_Leptons_trk_nPixelBarrelLayers',\n",
       " 'secVtx_VSI_Leptons_trk_nPixelEndCapLayers',\n",
       " 'secVtx_VSI_Leptons_trk_nSCTBarrelLayers',\n",
       " 'secVtx_VSI_Leptons_trk_nSCTEndCapLayers',\n",
       " 'secVtx_VSI_Leptons_trk_hitPattern',\n",
       " 'secVtx_VSI_Leptons_trk_dEdx',\n",
       " 'secVtx_VSI_Leptons_trk_radiusOfFirstHit',\n",
       " 'secVtx_VSI_Leptons_trk_truthProb',\n",
       " 'secVtx_VSI_Leptons_trk_truthID',\n",
       " 'secVtx_VSI_Leptons_trk_truthBarcode',\n",
       " 'secVtx_VSI_Leptons_trk_truthPdgId',\n",
       " 'secVtx_VSI_Leptons_trk_truthPointsToClosestTV',\n",
       " 'secVtx_VSI_Leptons_trk_truthPointsToMaxlinkTV',\n",
       " 'secVtx_VSI_Leptons_trk_truthParentPointsToMaxlinkPTV',\n",
       " 'secVtx_VSI_Leptons_closeTruth_distance',\n",
       " 'secVtx_VSI_Leptons_closeTruth_ID',\n",
       " 'secVtx_VSI_Leptons_closeTruth_barcode',\n",
       " 'secVtx_VSI_Leptons_closestTruth_distance',\n",
       " 'secVtx_VSI_Leptons_closestTruth_ID',\n",
       " 'secVtx_VSI_Leptons_closestTruth_x',\n",
       " 'secVtx_VSI_Leptons_closestTruth_y',\n",
       " 'secVtx_VSI_Leptons_closestTruth_z',\n",
       " 'secVtx_VSI_Leptons_closestTruth_r',\n",
       " 'secVtx_VSI_Leptons_closestTruth_eta',\n",
       " 'secVtx_VSI_Leptons_closestTruth_phi',\n",
       " 'secVtx_VSI_Leptons_closestTruth_barcode',\n",
       " 'secVtx_VSI_Leptons_linkTruth_score',\n",
       " 'secVtx_VSI_Leptons_linkTruth_ID',\n",
       " 'secVtx_VSI_Leptons_linkTruth_barcode',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_score',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_ID',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_x',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_y',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_z',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_r',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_eta',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_phi',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_barcode',\n",
       " 'secVtx_VSI_Leptons_maxlinkTruth_parent_pdgId',\n",
       " 'nxAH_truth',\n",
       " 'xAH_truth_pdgId',\n",
       " 'xAH_truth_status',\n",
       " 'xAH_truth_barcode',\n",
       " 'ntruthVtx',\n",
       " 'truthVtx_ID',\n",
       " 'truthVtx_x',\n",
       " 'truthVtx_y',\n",
       " 'truthVtx_z',\n",
       " 'truthVtx_r',\n",
       " 'truthVtx_pt',\n",
       " 'truthVtx_eta',\n",
       " 'truthVtx_phi',\n",
       " 'truthVtx_mass',\n",
       " 'truthVtx_nOutP',\n",
       " 'truthVtx_barcode',\n",
       " 'truthVtx_parent_ID',\n",
       " 'truthVtx_parent_pt',\n",
       " 'truthVtx_parent_eta',\n",
       " 'truthVtx_parent_phi',\n",
       " 'truthVtx_parent_E',\n",
       " 'truthVtx_parent_M',\n",
       " 'truthVtx_parent_charge',\n",
       " 'truthVtx_parent_pdgId',\n",
       " 'truthVtx_parent_status',\n",
       " 'truthVtx_parent_barcode',\n",
       " 'truthVtx_parent_prod_ID',\n",
       " 'truthVtx_parent_prod_x',\n",
       " 'truthVtx_parent_prod_y',\n",
       " 'truthVtx_parent_prod_z',\n",
       " 'truthVtx_parent_prod_r',\n",
       " 'truthVtx_parent_prod_eta',\n",
       " 'truthVtx_parent_prod_phi',\n",
       " 'truthVtx_parent_prod_barcode',\n",
       " 'truthVtx_outP_ID',\n",
       " 'truthVtx_outP_pt',\n",
       " 'truthVtx_outP_eta',\n",
       " 'truthVtx_outP_phi',\n",
       " 'truthVtx_outP_E',\n",
       " 'truthVtx_outP_M',\n",
       " 'truthVtx_outP_charge',\n",
       " 'truthVtx_outP_pdgId',\n",
       " 'truthVtx_outP_status',\n",
       " 'truthVtx_outP_barcode',\n",
       " 'truthVtx_outP_isReco',\n",
       " 'truthVtx_outP_recoProb',\n",
       " 'truthVtx_outP_recoID',\n",
       " 'truthVtx_outP_recoIsSelected',\n",
       " 'truthVtx_outP_recoIsAssociated',\n",
       " 'truthVtx_outP_isStable',\n",
       " 'truthVtx_outP_isInteracting',\n",
       " 'truthVtx_outP_isReconstructible',\n",
       " 'truthVtx_linkedRecoVtx_ID',\n",
       " 'truthVtx_linkedRecoVtx_score',\n",
       " 'truthVtx_maxlinkedRecoVtx_ID',\n",
       " 'truthVtx_maxlinkedRecoVtx_score']"
      ]
     },
     "execution_count": 70,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(old_df.columns)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [
    {
     "ename": "ValueError",
     "evalue": "Lengths must match to compare",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mValueError\u001b[0m                                Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-72-e2bc2308b2ab>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mold_df\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msecVtx_VSI_maxlinkTruth_barcode\u001b[0m \u001b[0;34m!=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/pandas/core/ops/common.py\u001b[0m in \u001b[0;36mnew_method\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m     62\u001b[0m         \u001b[0mother\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mitem_from_zerodim\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mother\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     63\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 64\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mmethod\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mother\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     65\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     66\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0mnew_method\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/pandas/core/ops/__init__.py\u001b[0m in \u001b[0;36mwrapper\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m    527\u001b[0m         \u001b[0mrvalues\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mextract_array\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mother\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mextract_numpy\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    528\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 529\u001b[0;31m         \u001b[0mres_values\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcomparison_op\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlvalues\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrvalues\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mop\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    530\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    531\u001b[0m         \u001b[0;32mreturn\u001b[0m \u001b[0m_construct_result\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mres_values\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mindex\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mindex\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mname\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mres_name\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/pandas/core/ops/array_ops.py\u001b[0m in \u001b[0;36mcomparison_op\u001b[0;34m(left, right, op)\u001b[0m\n\u001b[1;32m    232\u001b[0m         \u001b[0;31m#  The ambiguous case is object-dtype.  See GH#27803\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    233\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlvalues\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m!=\u001b[0m \u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 234\u001b[0;31m             \u001b[0;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Lengths must match to compare\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    235\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    236\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mshould_extension_dispatch\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlvalues\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrvalues\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mValueError\u001b[0m: Lengths must match to compare"
     ]
    }
   ],
   "source": [
    "old_df.secVtx_VSI_maxlinkTruth_barcode"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {},
   "outputs": [],
   "source": [
    "nicky_old = uproot.open('/data/hnl/forRobin/mc16_13TeV.311636.Pythia8EvtGen_A14NNPDF23LO_WmuHNL50_10G_lt10dd_el.merge.DAOD_RPVLL.e7422_e5984_a875_r10790_r10726/data-tree/mc16eONLY.311636.promptmuHNL_10G_10mm_el_DAOD_RPVLL.root')['outTree']\n",
    "nicky_old_df = nicky_old.pandas.df('el_*')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>el_LHLoose</th>\n",
       "      <th>el_LHMedium</th>\n",
       "      <th>el_LHTight</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>entry</th>\n",
       "      <th>subentry</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"5\" valign=\"top\">0</th>\n",
       "      <th>0</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19998</th>\n",
       "      <th>0</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"4\" valign=\"top\">19999</th>\n",
       "      <th>0</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>101676 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                el_LHLoose  el_LHMedium  el_LHTight\n",
       "entry subentry                                     \n",
       "0     0                 -1           -1          -1\n",
       "      1                 -1           -1          -1\n",
       "      2                 -1           -1          -1\n",
       "      3                 -1           -1          -1\n",
       "      4                 -1           -1          -1\n",
       "...                    ...          ...         ...\n",
       "19998 0                 -1           -1          -1\n",
       "19999 0                 -1           -1          -1\n",
       "      1                 -1           -1          -1\n",
       "      2                 -1           -1          -1\n",
       "      3                 -1           -1          -1\n",
       "\n",
       "[101676 rows x 3 columns]"
      ]
     },
     "execution_count": 84,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nicky_old_df[['el_LHLoose', 'el_LHMedium', 'el_LHTight']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 143,
   "metadata": {},
   "outputs": [
    {
     "ename": "FileNotFoundError",
     "evalue": "[Errno 2] No such file or directory: '/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new/data-tree/files_new.root'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mFileNotFoundError\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-143-00ba258359ef>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mnicky_new\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0muproot\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new/data-tree/files_new.root\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'outTree'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      2\u001b[0m \u001b[0mnicky_new_df\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnicky_new\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpandas\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdf\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'el_*'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/uproot/rootio.py\u001b[0m in \u001b[0;36mopen\u001b[0;34m(path, localsource, xrootdsource, httpsource, **options)\u001b[0m\n\u001b[1;32m     52\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     53\u001b[0m             \u001b[0mopenfcn\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlocalsource\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 54\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mROOTDirectory\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mread\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mopenfcn\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0moptions\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     55\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     56\u001b[0m     \u001b[0;32melif\u001b[0m \u001b[0m_bytesid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mparsed\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mscheme\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34mb\"root\"\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/uproot/rootio.py\u001b[0m in \u001b[0;36m<lambda>\u001b[0;34m(path)\u001b[0m\n\u001b[1;32m     49\u001b[0m                 \u001b[0;32mif\u001b[0m \u001b[0mn\u001b[0m \u001b[0;32min\u001b[0m \u001b[0moptions\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     50\u001b[0m                     \u001b[0mkwargs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mn\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0moptions\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpop\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mn\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 51\u001b[0;31m             \u001b[0mopenfcn\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mlambda\u001b[0m \u001b[0mpath\u001b[0m\u001b[0;34m:\u001b[0m \u001b[0mMemmapSource\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     52\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     53\u001b[0m             \u001b[0mopenfcn\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlocalsource\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/uproot/source/memmap.py\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, path)\u001b[0m\n\u001b[1;32m     19\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__init__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     20\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mexpanduser\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 21\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_source\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmemmap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdtype\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mnumpy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0muint8\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"r\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     22\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclosed\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mFalse\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     23\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/miniconda3/lib/python3.7/site-packages/numpy/core/memmap.py\u001b[0m in \u001b[0;36m__new__\u001b[0;34m(subtype, filename, dtype, mode, offset, shape, order)\u001b[0m\n\u001b[1;32m    223\u001b[0m             \u001b[0mf_ctx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcontextlib_nullcontext\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    224\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 225\u001b[0;31m             \u001b[0mf_ctx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mos_fspath\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfilename\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;34m'r'\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mmode\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;34m'c'\u001b[0m \u001b[0;32melse\u001b[0m \u001b[0mmode\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0;34m'b'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    226\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    227\u001b[0m         \u001b[0;32mwith\u001b[0m \u001b[0mf_ctx\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mfid\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mFileNotFoundError\u001b[0m: [Errno 2] No such file or directory: '/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new/data-tree/files_new.root'"
     ]
    }
   ],
   "source": [
    "nicky_new = uproot.open(\"/home/newhouse/public/Analysis/HNL/DHNLAlgorithm/run/testRun_new/data-tree/files_new.root\")['outTree']\n",
    "nicky_new_df = nicky_new.pandas.df('el_*')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 142,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>el_LHLoose</th>\n",
       "      <th>el_LHMedium</th>\n",
       "      <th>el_LHTight</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>entry</th>\n",
       "      <th>subentry</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"5\" valign=\"top\">0</th>\n",
       "      <th>0</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th rowspan=\"5\" valign=\"top\">99</th>\n",
       "      <th>9</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "      <td>-1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>737 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                el_LHLoose  el_LHMedium  el_LHTight\n",
       "entry subentry                                     \n",
       "0     0                 -1           -1          -1\n",
       "      1                 -1           -1          -1\n",
       "      2                 -1           -1          -1\n",
       "      3                 -1           -1          -1\n",
       "      4                 -1           -1          -1\n",
       "...                    ...          ...         ...\n",
       "99    9                 -1           -1          -1\n",
       "      10                -1           -1          -1\n",
       "      11                -1           -1          -1\n",
       "      12                -1           -1          -1\n",
       "      13                -1           -1          -1\n",
       "\n",
       "[737 rows x 3 columns]"
      ]
     },
     "execution_count": 142,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "pd.set_option('display.max_rows', 500)\n",
    "nicky_new_df[['el_LHLoose', 'el_LHMedium', 'el_LHTight']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 113,
   "metadata": {},
   "outputs": [],
   "source": [
    "nicky_old = uproot.open('/data/hnl/forRobin/mc16_13TeV.311636.Pythia8EvtGen_A14NNPDF23LO_WmuHNL50_10G_lt10dd_el.merge.DAOD_RPVLL.e7422_e5984_a875_r10790_r10726/data-tree/mc16eONLY.311636.promptmuHNL_10G_10mm_el_DAOD_RPVLL.root')['outTree']\n",
    "nicky_old_df = nicky_old.pandas.df('el_*', entrystop=30, flatten=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 135,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "          el_LHLoose  el_LHMedium  el_LHTight\n",
      "subentry                                     \n",
      "0                 -1            1           0\n",
      "1                 -1           -1          -1\n",
      "2                 -1           -1          -1\n",
      "3                 -1           -1          -1\n",
      "el_LHLoose[0] -1: accepted\n",
      "el_LHLoose[1] -1: accepted\n",
      "el_LHLoose[2] -1: accepted\n",
      "el_LHLoose[3] -1: accepted\n",
      "el_LHMedium[0] 1: accepted\n",
      "el_LHMedium[1] -1: accepted\n",
      "el_LHMedium[2] -1: accepted\n",
      "el_LHMedium[3] -1: accepted\n",
      "el_LHTight[0] 0: rejected\n",
      "el_LHTight[1] -1: accepted\n",
      "el_LHTight[2] -1: accepted\n",
      "el_LHTight[3] -1: accepted\n"
     ]
    }
   ],
   "source": [
    "electron_quality = nicky_old_df[['el_LHLoose', 'el_LHMedium', 'el_LHTight']].loc[20]\n",
    "print(electron_quality)\n",
    "for wp in electron_quality:\n",
    "    for i, val in enumerate(electron_quality[wp]):\n",
    "        print(f\"{wp}[{i}]\", val, end=': ')\n",
    "        if val: print('accepted') \n",
    "        else: print(\"rejected\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 154,
   "metadata": {},
   "outputs": [],
   "source": [
    "data_ntuple = uproot.open(\"/data/hnl/data_rep_ntuples/alignment_off/user.rnewhous.data17_13TeV.00340368.physics_Main.merge.DAOD_RPVLL.r11761_r11764_p4054_HNLNtuple_01_tree.root/user.rnewhous.21710573._000001.tree.root\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 165,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "count    42857\n",
       "freq     38375\n",
       "Name: passesHnlMuMuFilter, dtype: object"
      ]
     },
     "execution_count": 165,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data_ntuple['outTree'].pandas.df(\"*Filter*\").describe()['passesHnlMuMuFilter'].loc[['count', 'freq']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}