hzz4l_analysis.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# HZZ4l mini analysis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First a few useful Jupyter shortcuts:\n",
    "\n",
    "* Enter: edit a cell\n",
    "* Esc: escape from edit mode\n",
    "* b: Add a cell below\n",
    "* a: Add a cell above\n",
    "* m: switch to Markdown cell\n",
    "* Ctrl-Enter: Run cell\n",
    "* Alt-Enter: Run cell and add cell below\n",
    "* Shift-Enter: Run cell and move below"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "require(['codemirror/mode/clike/clike'], function(Clike) { console.log('ROOTaaS - C++ CodeMirror module loaded'); });"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": [
       "IPython.CodeCell.config_defaults.highlight_modes['magic_text/x-c++src'] = {'reg':[/^%%cpp/]};"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Welcome to ROOTaaS 6.06/00\n"
     ]
    }
   ],
   "source": [
    "import ROOT"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are signal and background simulated data we are interested in"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ZZ4lAnalysis.root\n",
      "ZZ4lAnalysis.root\n"
     ]
    }
   ],
   "source": [
    "# Signal tree\n",
    "! ls /data_CMS/cms/ochando/CJLSTReducedTree/160726/ZZTo4l\n",
    "# Background tree\n",
    "! ls /data_CMS/cms/ochando/CJLSTReducedTree/160726/ggH125"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's look at what is inside the files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "TFile**\t\t/data_CMS/cms/ochando/CJLSTReducedTree/160726/ggH125/ZZ4lAnalysis.root\t\r\n",
      " TFile*\t\t/data_CMS/cms/ochando/CJLSTReducedTree/160726/ggH125/ZZ4lAnalysis.root\t\r\n",
      "  TDirectoryFile*\t\tZZTree\tZZTree\r\n",
      "   KEY: TTree\tcandTree;1\tEvent Summary\r\n",
      "   KEY: TH1F\tCounters;1\tCounters\r\n",
      "   KEY: TH2F\tCounters_reweighted;1\tCounters_reweighted\r\n",
      "   KEY: TTree\tcouplings;1\treweighting couplings\r\n",
      "  KEY: TDirectoryFile\tZZTree;1\tZZTree\r\n",
      "  KEY: TDirectoryFile\tPlotsZZ;1\tPlotsZZ\r\n"
     ]
    }
   ],
   "source": [
    "signal_file = ROOT.TFile.Open('/data_CMS/cms/ochando/CJLSTReducedTree/160726/ggH125/ZZ4lAnalysis.root')\n",
    "signal_file.cd('ZZTree')\n",
    "signal_file.ls()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We retrieve the signal ROOT tree where the events are stored. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "signal_tree = signal_file.Get('ZZTree/candTree')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We load the first event and look at the variables values. Type Esc-o to toggle the output cell."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "======> EVENT:0\r\n",
      " RunNumber       = 1\r\n",
      " EventNumber     = 254049\r\n",
      " LumiNumber      = 1278\r\n",
      " NRecoMu         = 2\r\n",
      " NRecoEle        = 2\r\n",
      " Nvtx            = 12\r\n",
      " NObsInt         = 15\r\n",
      " NTrueInt        = 20.1439\r\n",
      " PFMET           = 23.4662\r\n",
      " PFMETPhi        = -2.08841\r\n",
      " PFMETNoHF       = 43.1148\r\n",
      " PFMETNoHFPhi    = -1.85145\r\n",
      " nCleanedJets    = 2\r\n",
      " nCleanedJetsPt30 = 1\r\n",
      " nCleanedJetsPt30_jecUp = 1\r\n",
      " nCleanedJetsPt30_jecDn = 1\r\n",
      " nCleanedJetsPt30BTagged = 0\r\n",
      " nCleanedJetsPt30BTagged_bTagSF = 0\r\n",
      " nCleanedJetsPt30BTagged_bTagSFUp = 0\r\n",
      " nCleanedJetsPt30BTagged_bTagSFDn = 0\r\n",
      " trigWord        = -32565\r\n",
      " ZZMass          = 119.048\r\n",
      " ZZMassErr       = 1.17973\r\n",
      " ZZMassErrCorr   = 1.17973\r\n",
      " ZZMassPreFSR    = 119.048\r\n",
      " ZZsel           = 100\r\n",
      " ZZPt            = 105.919\r\n",
      " ZZEta           = 0.323407\r\n",
      " ZZPhi           = -2.33697\r\n",
      " CRflag          = 805306368\r\n",
      " Z1Mass          = 83.7659\r\n",
      " Z1Pt            = 71.7932\r\n",
      " Z1Flav          = -169\r\n",
      " ZZMassRefit     = 119.327\r\n",
      " ZZMassRefitErr  = 1.44506\r\n",
      " ZZMassUnrefitErr = 1.46347\r\n",
      " Z2Mass          = 26.9499\r\n",
      " Z2Pt            = 34.1295\r\n",
      " Z2Flav          = -121\r\n",
      " costhetastar    = 0.964909\r\n",
      " helphi          = 0.106953\r\n",
      " helcosthetaZ1   = -0.871139\r\n",
      " helcosthetaZ2   = -0.290205\r\n",
      " phistarZ1       = -0.773818\r\n",
      " phistarZ2       = 0\r\n",
      " xi              = -2.33697\r\n",
      " xistar          = 0.724765\r\n",
      " LepPt           = (vector<float>*)0x5814420\r\n",
      " LepEta          = (vector<float>*)0x578cdb0\r\n",
      " LepPhi          = (vector<float>*)0x5792450\r\n",
      " LepLepId        = (vector<short>*)0x57821f0\r\n",
      " LepSIP          = (vector<float>*)0x572c9d0\r\n",
      " LepTime         = (vector<float>*)0x579bdf0\r\n",
      " LepisID         = (vector<bool>*)0x5afe4c0\r\n",
      " LepisLoose      = (vector<short>*)0x575fd00\r\n",
      " LepBDT          = (vector<float>*)0x5814e90\r\n",
      " LepMissingHit   = (vector<char>*)0x5b40e20\r\n",
      " LepCombRelIsoPF = (vector<float>*)0x5b3f800\r\n",
      " fsrPt           = (vector<float>*)0x5b40430\r\n",
      " fsrEta          = (vector<float>*)0x5b21c30\r\n",
      " fsrPhi          = (vector<float>*)0x5ac8510\r\n",
      " fsrLept         = (vector<short>*)0x5ac78b0\r\n",
      " passIsoPreFSR   = 1\r\n",
      " p0plus_VAJHU    = 0.180621\r\n",
      " p0_g1prime2_VAJHU = 1.0829e-09\r\n",
      " p0hplus_VAJHU   = 0.0936451\r\n",
      " p0minus_VAJHU   = 0.0406323\r\n",
      " p0_g1prime2_zgs_VAJHU = 2.36812e-09\r\n",
      " p0hplus_zgs_VAJHU = 28.8962\r\n",
      " p0minus_zgs_VAJHU = 12.6586\r\n",
      " p0hplus_gsgs_VAJHU = 2.89802\r\n",
      " p0minus_gsgs_VAJHU = 1.23666\r\n",
      " pg1g1prime2_VAJHU = 2.79695e-05\r\n",
      " pg1g2_VAJHU     = 0.246853\r\n",
      " pg1g2_pi2_VAJHU = 0.00152239\r\n",
      " pg1g4_VAJHU     = 0.000651017\r\n",
      " pg1g4_pi2_VAJHU = 0.0175876\r\n",
      " p0plus_zz_g1prime2_zgs_VAJHU = -1.62423e-06\r\n",
      " p0plus_zz_g1prime2_zgs_pi2_VAJHU = -7.45058e-08\r\n",
      " p0plus_zz_0hplus_zgs_VAJHU = 0.0156174\r\n",
      " p0plus_zz_0minus_zgs_VAJHU = -0.0331306\r\n",
      " p0plus_zz_0hplus_gsgs_VAJHU = -0.63163\r\n",
      " p0plus_zz_0minus_gsgs_VAJHU = -0.076394\r\n",
      " p1_VAJHU        = 1.33771\r\n",
      " p1_prodIndep_VAJHU = 3.17499e-06\r\n",
      " p1plus_VAJHU    = 1.79188\r\n",
      " p1plus_prodIndep_VAJHU = 4.30847e-06\r\n",
      " p2plus_gg_VAJHU = 12.1221\r\n",
      " p2plus_prodIndep_VAJHU = 3.03784e-05\r\n",
      " p2plus_qqb_VAJHU = 0.848274\r\n",
      " p2h2plus_gg_VAJHU = 8.04719e-09\r\n",
      " p2h2plus_qqb_VAJHU = 7.47192e-06\r\n",
      " p2h2plus_prodIndep_VAJHU = 2.95759e-10\r\n",
      " p2h3plus_gg_VAJHU = 1.26054e-11\r\n",
      " p2h3plus_qqb_VAJHU = 1.98745e-08\r\n",
      " p2h3plus_prodIndep_VAJHU = 4.43768e-13\r\n",
      " p2h4plus_gg_VAJHU = 2.46815e-10\r\n",
      " p2h4plus_qqb_VAJHU = 4.59524e-08\r\n",
      " p2h4plus_prodIndep_VAJHU = 1.47867e-12\r\n",
      " p2bplus_gg_VAJHU = 6.56845\r\n",
      " p2bplus_qqb_VAJHU = 0.532777\r\n",
      " p2bplus_prodIndep_VAJHU = 1.81795e-05\r\n",
      " p2h6plus_gg_VAJHU = 5.8327e-06\r\n",
      " p2h6plus_qqb_VAJHU = 3.06779e-06\r\n",
      " p2h6plus_prodIndep_VAJHU = 7.82215e-11\r\n",
      " p2h7plus_gg_VAJHU = 3.57955e-08\r\n",
      " p2h7plus_qqb_VAJHU = 8.86321e-08\r\n",
      " p2h7plus_prodIndep_VAJHU = 2.85202e-12\r\n",
      " p2hminus_gg_VAJHU = 1.07092e-10\r\n",
      " p2hminus_qqb_VAJHU = 1.99386e-08\r\n",
      " p2hminus_prodIndep_VAJHU = 6.41588e-13\r\n",
      " p2h9minus_gg_VAJHU = 5.44114e-09\r\n",
      " p2h9minus_qqb_VAJHU = 4.11016e-06\r\n",
      " p2h9minus_prodIndep_VAJHU = 7.31853e-11\r\n",
      " p2h10minus_gg_VAJHU = 1.44112e-14\r\n",
      " p2h10minus_qqb_VAJHU = 2.06738e-11\r\n",
      " p2h10minus_prodIndep_VAJHU = 3.2332e-16\r\n",
      " p0plus_VAMCFM   = 0.182391\r\n",
      " ggzz_VAMCFM     = 0.155727\r\n",
      " ggzz_p0plus_VAMCFM = 0.155218\r\n",
      " bkg_VAMCFM      = 0.0732347\r\n",
      " bkg_prodIndep_VAMCFM = 2.7394e-09\r\n",
      " pZJJ_VAMCFM     = 233.498\r\n",
      " Dgg10_VAMCFM    = 2.94708e-11\r\n",
      " p0plus_m4l      = 0.0188059\r\n",
      " p0plus_m4l_ScaleUp = 0.0214892\r\n",
      " p0plus_m4l_ScaleDown = 0.0165231\r\n",
      " p0plus_m4l_ResUp = 0.0183345\r\n",
      " p0plus_m4l_ResDown = 0.0163023\r\n",
      " bkg_m4l         = 0.000717506\r\n",
      " bkg_m4l_ScaleUp = 0.000723788\r\n",
      " bkg_m4l_ScaleDown = 0.0007109\r\n",
      " bkg_m4l_ResUp   = 0.000724485\r\n",
      " bkg_m4l_ResDown = 0.000714229\r\n",
      " pwh_leptonic_VAJHU = -1\r\n",
      " pzh_leptonic_VAJHU = -1\r\n",
      " phjj_VAJHU_highestPTJets = -1\r\n",
      " pvbf_VAJHU_highestPTJets = 0.00534728\r\n",
      " phjj_VAJHU_highestPTJets_up = -1\r\n",
      " pvbf_VAJHU_highestPTJets_up = 0.0053556\r\n",
      " phjj_VAJHU_highestPTJets_dn = -1\r\n",
      " pvbf_VAJHU_highestPTJets_dn = 0.00531746\r\n",
      " phjj_VAJHU_bestDjet = -1\r\n",
      " pvbf_VAJHU_bestDjet = 0.00534728\r\n",
      " phjj_VAJHU_bestDjet_up = -1\r\n",
      " pvbf_VAJHU_bestDjet_up = 0.0053556\r\n",
      " phjj_VAJHU_bestDjet_dn = -1\r\n",
      " pvbf_VAJHU_bestDjet_dn = 0.00531746\r\n",
      " pAux_vbf_VAJHU  = 15.0424\r\n",
      " pAux_vbf_VAJHU_up = 14.6885\r\n",
      " pAux_vbf_VAJHU_dn = 15.4201\r\n",
      " phj_VAJHU       = 0.0600557\r\n",
      " phj_VAJHU_up    = 0.0594193\r\n",
      " phj_VAJHU_dn    = 0.0607074\r\n",
      " pwh_hadronic_VAJHU = -1\r\n",
      " pwh_hadronic_VAJHU_up = -1\r\n",
      " pwh_hadronic_VAJHU_dn = -1\r\n",
      " pzh_hadronic_VAJHU = -1\r\n",
      " pzh_hadronic_VAJHU_up = -1\r\n",
      " pzh_hadronic_VAJHU_dn = -1\r\n",
      " ptth_VAJHU      = -1\r\n",
      " ptth_VAJHU_up   = -1\r\n",
      " ptth_VAJHU_dn   = -1\r\n",
      " pbbh_VAJHU      = -1\r\n",
      " pbbh_VAJHU_up   = -1\r\n",
      " pbbh_VAJHU_dn   = -1\r\n",
      " JetPt           = (vector<float>*)0x5eb5ef0\r\n",
      " JetEta          = (vector<float>*)0x5ebf530\r\n",
      " JetPhi          = (vector<float>*)0x5ebf930\r\n",
      " JetMass         = (vector<float>*)0x5ebfd00\r\n",
      " JetBTagger      = (vector<float>*)0x5edb910\r\n",
      " JetIsBtagged    = (vector<float>*)0x5edbd10\r\n",
      " JetIsBtaggedWithSF = (vector<float>*)0x5edc0e0\r\n",
      " JetIsBtaggedWithSFUp = (vector<float>*)0x5ef7d20\r\n",
      " JetIsBtaggedWithSFDn = (vector<float>*)0x5ef8180\r\n",
      " JetQGLikelihood = (vector<float>*)0x5f0ac00\r\n",
      " JetAxis2        = (vector<float>*)0x5f14250\r\n",
      " JetMult         = (vector<float>*)0x5f14650\r\n",
      " JetPtD          = (vector<float>*)0x5f27020\r\n",
      " JetSigma        = (vector<float>*)0x5f30660\r\n",
      " JetHadronFlavour = (vector<short>*)0x5f30a60\r\n",
      " DiJetMass       = -99\r\n",
      " DiJetDEta       = -99\r\n",
      " DiJetFisher     = -99\r\n",
      " nExtraLep       = 0\r\n",
      " nExtraZ         = 0\r\n",
      " ExtraLepPt      = (vector<float>*)0x5f6b290\r\n",
      " ExtraLepEta     = (vector<float>*)0x5f74ee0\r\n",
      " ExtraLepPhi     = (vector<float>*)0x5f7eba0\r\n",
      " ExtraLepLepId   = (vector<short>*)0x5f88860\r\n",
      " ZXFakeweight    = 1\r\n",
      " KFactor_QCD_ggZZ_Nominal = 1.25693\r\n",
      " KFactor_QCD_ggZZ_PDFScaleDn = 1.24773\r\n",
      " KFactor_QCD_ggZZ_PDFScaleUp = 1.26709\r\n",
      " KFactor_QCD_ggZZ_QCDScaleDn = 1.36823\r\n",
      " KFactor_QCD_ggZZ_QCDScaleUp = 1.12941\r\n",
      " KFactor_QCD_ggZZ_AsDn = 1.225\r\n",
      " KFactor_QCD_ggZZ_AsUp = 1.28233\r\n",
      " KFactor_QCD_ggZZ_PDFReplicaDn = 1.24056\r\n",
      " KFactor_QCD_ggZZ_PDFReplicaUp = 1.27715\r\n",
      " genFinalState   = 2\r\n",
      " genProcessId    = 9999\r\n",
      " genHEPMCweight  = 1\r\n",
      " PUWeight        = 1.11485\r\n",
      " dataMCWeight    = 0.845026\r\n",
      " trigEffWeight   = 1\r\n",
      " overallEventWeight = 0.942077\r\n",
      " HqTMCweight     = 1\r\n",
      " xsec            = 0.0134081\r\n",
      " genExtInfo      = -1\r\n",
      " GenHMass        = 125.002\r\n",
      " GenHPt          = 101.222\r\n",
      " GenHRapidity    = 0.295107\r\n",
      " GenZ1Mass       = 90.1455\r\n",
      " GenZ1Pt         = 79.1216\r\n",
      " GenZ1Phi        = -2.37552\r\n",
      " GenZ1Flav       = -169\r\n",
      " GenZ2Mass       = 27.0515\r\n",
      " GenZ2Pt         = 34.2134\r\n",
      " GenZ2Phi        = -2.3526\r\n",
      " GenZ2Flav       = -121\r\n",
      " GenLep1Pt       = 65.74\r\n",
      " GenLep1Eta      = -0.204061\r\n",
      " GenLep1Phi      = -2.18823\r\n",
      " GenLep1Id       = 13\r\n",
      " GenLep2Pt       = 18.9999\r\n",
      " GenLep2Eta      = 1.82906\r\n",
      " GenLep2Phi      = -3.07557\r\n",
      " GenLep2Id       = -13\r\n",
      " GenLep3Pt       = 28.4689\r\n",
      " GenLep3Eta      = 0.115819\r\n",
      " GenLep3Phi      = -2.08893\r\n",
      " GenLep3Id       = 11\r\n",
      " GenLep4Pt       = 10.016\r\n",
      " GenLep4Eta      = -1.03506\r\n",
      " GenLep4Phi      = 3.09637\r\n",
      " GenLep4Id       = -11\r\n",
      " GenAssocLep1Pt  = 0\r\n",
      " GenAssocLep1Eta = 0\r\n",
      " GenAssocLep1Phi = 0\r\n",
      " GenAssocLep1Id  = 0\r\n",
      " GenAssocLep2Pt  = 0\r\n",
      " GenAssocLep2Eta = 0\r\n",
      " GenAssocLep2Phi = 0\r\n",
      " GenAssocLep2Id  = 0\r\n",
      " reweightingweights = (vector<float>*)0x6147a20\r\n",
      " LHEPDFScale     = 125.002\r\n",
      " LHEweight_QCDscale_muR1_muF1 = 1\r\n",
      " LHEweight_QCDscale_muR1_muF2 = 0.990736\r\n",
      " LHEweight_QCDscale_muR1_muF0p5 = 0.997017\r\n",
      " LHEweight_QCDscale_muR2_muF1 = 0.752518\r\n",
      " LHEweight_QCDscale_muR2_muF2 = 0.745547\r\n",
      " LHEweight_QCDscale_muR2_muF0p5 = 0.750273\r\n",
      " LHEweight_QCDscale_muR0p5_muF1 = 1.37145\r\n",
      " LHEweight_QCDscale_muR0p5_muF2 = 1.35874\r\n",
      " LHEweight_QCDscale_muR0p5_muF0p5 = 1.36736\r\n",
      " LHEweight_PDFVariation_Up = 0.984894\r\n",
      " LHEweight_PDFVariation_Dn = 1.02047\r\n",
      " LHEweight_AsMZ_Up = 1\r\n",
      " LHEweight_AsMZ_Dn = 1\r\n"
     ]
    }
   ],
   "source": [
    "signal_tree.GetEntry(0)\n",
    "signal_tree.Show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We choose to read only a subset of the variables, that we are going to use"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "signal_tree.SetBranchStatus('*', 0)\n",
    "signal_tree.SetBranchStatus('GenLep*',1)\n",
    "signal_tree.SetBranchStatus('GenZ*',1)\n",
    "signal_tree.SetBranchStatus('GenH*',1)\n",
    "signal_tree.SetBranchStatus('Lep*',1)\n",
    "signal_tree.SetBranchStatus('Z*',1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "======> EVENT:0\r\n",
      " ZZMass          = 119.048\r\n",
      " ZZMassErr       = 1.17973\r\n",
      " ZZMassErrCorr   = 1.17973\r\n",
      " ZZMassPreFSR    = 119.048\r\n",
      " ZZsel           = 100\r\n",
      " ZZPt            = 105.919\r\n",
      " ZZEta           = 0.323407\r\n",
      " ZZPhi           = -2.33697\r\n",
      " Z1Mass          = 83.7659\r\n",
      " Z1Pt            = 71.7932\r\n",
      " Z1Flav          = -169\r\n",
      " ZZMassRefit     = 119.327\r\n",
      " ZZMassRefitErr  = 1.44506\r\n",
      " ZZMassUnrefitErr = 1.46347\r\n",
      " Z2Mass          = 26.9499\r\n",
      " Z2Pt            = 34.1295\r\n",
      " Z2Flav          = -121\r\n",
      " LepPt           = (vector<float>*)0x5814420\r\n",
      " LepEta          = (vector<float>*)0x578cdb0\r\n",
      " LepPhi          = (vector<float>*)0x5792450\r\n",
      " LepLepId        = (vector<short>*)0x57821f0\r\n",
      " LepSIP          = (vector<float>*)0x572c9d0\r\n",
      " LepTime         = (vector<float>*)0x579bdf0\r\n",
      " LepisID         = (vector<bool>*)0x5afe4c0\r\n",
      " LepisLoose      = (vector<short>*)0x575fd00\r\n",
      " LepBDT          = (vector<float>*)0x5814e90\r\n",
      " LepMissingHit   = (vector<char>*)0x5b40e20\r\n",
      " LepCombRelIsoPF = (vector<float>*)0x5b3f800\r\n",
      " ZXFakeweight    = 1\r\n",
      " GenHMass        = 125.002\r\n",
      " GenHPt          = 101.222\r\n",
      " GenHRapidity    = 0.295107\r\n",
      " GenZ1Mass       = 90.1455\r\n",
      " GenZ1Pt         = 79.1216\r\n",
      " GenZ1Phi        = -2.37552\r\n",
      " GenZ1Flav       = -169\r\n",
      " GenZ2Mass       = 27.0515\r\n",
      " GenZ2Pt         = 34.2134\r\n",
      " GenZ2Phi        = -2.3526\r\n",
      " GenZ2Flav       = -121\r\n",
      " GenLep1Pt       = 65.74\r\n",
      " GenLep1Eta      = -0.204061\r\n",
      " GenLep1Phi      = -2.18823\r\n",
      " GenLep1Id       = 13\r\n",
      " GenLep2Pt       = 18.9999\r\n",
      " GenLep2Eta      = 1.82906\r\n",
      " GenLep2Phi      = -3.07557\r\n",
      " GenLep2Id       = -13\r\n",
      " GenLep3Pt       = 28.4689\r\n",
      " GenLep3Eta      = 0.115819\r\n",
      " GenLep3Phi      = -2.08893\r\n",
      " GenLep3Id       = 11\r\n",
      " GenLep4Pt       = 10.016\r\n",
      " GenLep4Eta      = -1.03506\r\n",
      " GenLep4Phi      = 3.09637\r\n",
      " GenLep4Id       = -11\r\n"
     ]
    }
   ],
   "source": [
    "signal_tree.Show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Generated objects"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First we look at generated objects. Below we fill an histogram with the transverse momentum (pT) of the first lepton."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "127817L"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "histo_pt = ROOT.TH1F('histo_pt', 'histo_pt', 50, 0, 100)\n",
    "signal_tree.Project('histo_pt', 'GenLep1Pt')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElE\nQVR4nO3dza7rxoEu0OJFv4Xnto/jeQ9JzuO07cCN9Au45zZsxIgBkUCMbiRw5u0XiNFBp3/iuahh\n5onj43k/h+6g7qlbh1Jpa+9NURS1FgxDokidoqSt+lR/rPb7fQAAOOb/XLsAAMByCQoAQJGgAAAU\nCQoAQJGgAAv1xRdfVFX1wQcfnN7t7//+76uq+vDDD+cpFXBvBAVYqB9++OGc3f73f//3sc/c933f\n948v0cS++uqrr7766tqlAB4gKMBte+ONN954440//vGP5x/y17/+teu6f/zHf7xcqc7x7bfffvnl\nl9ctA/Cgv7t2AYBn+fOf/3ztIgBrpkUBbsMXX3zxwQcffPHFF2fu/+mnn7733nuffvrpo/6Vzz//\n/PPPP/+Hf/iHzz///PFlPOLLL7/8+c9/frTl4C9/+UsIQe8DLN0eWKT3338/hPD+++/HG7lf/vKX\nabc33ngjhPDBBx+kLT/96U9H+7/55pvp0dFDXdelh956663Ro5999tlji/2rX/0qhPDhhx8ejq/8\n1a9+Fff5xS9+4YsIboWuB1i0//qv/wohvP322++///4PP/wQ7/7rv/7rv/zLvxzd/7333vvuu+9C\nCCkufPfddz/++ONbb7318uXLEMJHH330hz/8IYTwzjvv/OQnP0kHVlUVb/zsZz8LIbx8+fLly5e/\n/e1vQwi/+c1vHlvsNGYixoV4NzYe/PrXv3733Xd/8YtffPvttyGEd9999913333s8wPzuXZSAY5L\nDQnvv/9+vn30lztqUYiPfvLJJ2mHTz75ZHTIRx99FEL46KOP0pYYDkZfCEc3Pii2KIQQXrx4kW9/\n8eLF6NliPnjUkwPzM0YBlu4///M/87tvv/12COH0YIWvv/46v/3JJ5+kuHDU//zP/4QQPvvss3zj\nf//3f8cbTxuv8P3335+4C9wKQQEWLcaC87355pshhKqq8mGMX3/9dR4dSg67GFI3xKPKUBIbFX7+\n859P8mzAPAQFWLTHBoVUqf/ud7+rqqqqqgfnPjzYYPCEoJA6Gh7cCCycoABrs9/vP/nkk9i0EEL4\n7rvvfve73x3OaAA4h6AAK/T111+/fPkyDkSK0x9+/PHHx66pkDwhZPztb387cyOwcIICrMpbb731\n3nvv5Vv+9Kc/paxw9JATsx/jIMepWiNiUNABAbdFUIBV+fHHH+M6CqONp4+KUWA0HiKNXXjCOgoh\nhHfeeSe/m8Yw/vrXv37CswHXYsElWJU333zzxx9/rKrqpz/9aRymEBdcCiH86U9/yvf861//mm7/\n8MMPVVW9fPny7bffTjMdYnNCWk3hsf72t7+98847acGl2JyQVllI/umf/un3v//90/4JYAaCAqzK\ny5cv4xqLo3aFfB2Fn/zkJ3/4wx++//77qqq6rttsNiGEzz777Le//W1ajTH62c9+llZTeJQXL168\nePHij3/8Y34phw8//DBvTnj33Xf/8pe/fPvtt99+++3+YG1pYCEqf5+wPp9++mnqbnjzzTcPF1Ho\n+z62KPz7v/97vv3zzz9PkyGfFhG+/PLLr7766sWLF99///2XX36ZBjD+x3/8x+HOX331Vbw0lEYF\nWCxBAZhSHhSuXRZgAgYzAgBFggLwsOps1y4pMDFBAXjY+Reau3ZJgYkZowAAFGlRAACKBAUAoMiC\nSwBMzLDWm3Dm2ANBAYDpGQC3cOeHOV0PAECRoADA2jRNMwzDYx/iKEEBgLXZ7XalNLDb7R48/KJh\n4uaSiqAAwB3Z7/dN05ze55ww8WQnQswyCQoArNMwDE3TjH7Bd12X79N1XdM0+cZ4exiGdNQwDIe7\nPfhPn/nkN+D8lVkB4BxXr1xCCHVdx//HG5vNJj203W73+/12u431YNwhlTndjYeUdjshHfLgk1/R\n+e+RoADAxJYQFPIyxOo5PRSDQr4xbj8ME/tXmSPfLb97VAwKpaOWkBL2j3mPdD0AsEKbzSbdPjEo\nIXUB7Pf7w56F1IOQP+2ZIxiedtQCCQoA3KMYAtq2raqqNBMhbsxzxoMDIY/uef5RCyQoAHCn9vv9\ndruNP/fbtj1/rOJdERQAuEdd18VpEV3X7ff7uq77vh/tE1sC8saG82crPO2oBRIUALhHwzC0bZvu\n7na7NCUhiUEh7TYMQ9/3h7sd9bSjFkhQAOAexY6G6pXw+vDD1BMRpzDEfdq2rev6zOaBuq7TUeH1\nRoW+729o1EK1d4EvACZVVTdTucS1j+K6TPn2uFZS2ljarfScbdvu9/vSUYdjJOd3/nt0M+8lALfi\nhoLCJaSgcO2CnHL+e/R3ly7KzM6/wDYAPE3XdYcjH5N8ZcYVWFvou/MYC7AEvoq7rlv4ZMv77Xrw\n6QS4Oo27N+FOux4AWAK/2Rbu/DBneiQAUCQoAABFggIAULTCoFAVXLtcAPeu67ppv5/jskjTFXBK\n8VoShxvj+kvpobQiUy7NmBiGIR5SmkORbz/9VE+2wsGMRtAALNljlxk4sX7RYlNCCCEutJCXMEai\nzWYTz2iz2RytxXe7XbwRd4sXiej7vu/70Yvw4Omnp3qW/bqs74wAbk7pq3iz2TzhWzoGi2cXaj7p\n+k+bzSZtHJ176aXIT7au67qu00P5E8bDQwj5Diee6tD5L+kKux4AuEVVVQ3DkPoj4q/tdI3H+Ghc\nyKhpmnQ3/1Udt0d52/7M3dBN06SKPBmGIb+A5OEFrKO2bVOLy263O9FmsNlsTl+RMn+qZzkzUNyK\n9Z0RwM0pfRXH6nN7IB0VH8333L/+4ztWjXVdHz6Ub4+75dv3r35hn/gJPq3weovCyNEWhVj3j/Y5\nfEGSUZPDiac6WrwTj76255n73QpBAeDqTgeF0k/WUc2a6sVRa3w41oB/2MyeDs/r1zyXXNqJoBCL\nffjo4euWv0qH+58ICoep4sF/q2SFgxkBWLJ9ecj5OYMTj7a3xzb80djAOAugrus4JHCSKQDPlHpS\nttvt6GRjUUdbwquXKx14zil0XRfPd5IyG6MAwAptNps0DiC2N/R9X1XVFSdKdF0XJzvs9/vDYux2\nu9Fcx91ulwYZxHEPJy5Zmev7fsJIpEUBgJXIa8c0zjG2K6QmhzPr2kuUre/7w4aE9GiYaLZnPNMJ\n85AWBQBmNRx4/nPGijYFhTwQtG179R6HEELf97Fn4eiJj+ZEhFc1fV7y9AynxX6HaQodQlhli0Jp\n9suJXjF4UFV98+A++/3HM5QEbl3sa8+Vfmfnqqo6Pdlvu922bZvyQb7SQFytaLR9frvdbnTuqWLa\n7XaHBYtnlCq1uq7PCVVHn+o5qpVVn1W1tjNiCc5JCZGsAOECX8XnN6eX9py8QX42KRxMW/jz36O1\nVauCApdwflA4TYzgTvgqXr7z36MVdj3A5ZRq+jOTRFV9IysAt8VgRpiA6h9YKy0KMI3TWWGqzguA\nmWlRAACKtCjArE43LejCAJZGiwIsiB4K1m10AegkXh56CcsizSZdI/vwrNP2E6smHD08bhl5/npW\nggLMQVMBJIdV1263u0ZBriZe4SmuLR2vQJEeqqqq7/tYwcd9ThweXl3AovQPTfLCrm2q64nXa2Vn\nypzSD/3L1fcz/BMwm9Ic/fQVnT+arou42WzupFEhXYQi3q2qKp57XHk6vTijuw8efnq3kfPXUVhh\ni8KJa28DcF1xdeG89uq6brTk8DAM1Sv5coT59rwLI7bAx/6L614f8nyHV8QOB1d8KJ3Ibrd78Bzj\nxScnuY7GCoMCAEs2umbBqNqLDQzxWszb7TZ/NN9e13Wqa3e7XWyuT4csvGUidjqk2+EgN4x2Hm2J\n15vODzzMDW3bnr40xiOUfn/fqPWdETMI4d/O/G+GMlzun4DZlL6KYyUXK7C4Jd2OD+33+7qu67pO\nh+Q7xx3S7bQ9hJAfkp5q4VJFngofT2q73ca7sXUh3c2NcsPo0fzFKTm/upxjemQcnDnamI/VjGM7\nDwd/plGd0eETpsuNw5OZaAAzS33n8Wt/1O8Qx9+l6mDUSZEu0DwappfXBdNeZPlyYhNIbEHpXkkj\nNkL5ROLJxso+7T+6IPVkzQnX6nrIO07iCM9wMHSzaZp48qORn/mLuJCrjHMnDDOEqaSOg6Pd7aN+\n+pQkqqpKkwWueLXoacVzSVfBHoZh+8rRCi6GpJQDRoeHC1wnc44WhdE8zrzij69CakWJ82jTpydd\noTyfUdq2berfSkFshrNg9eQAmE3Xdenbe1SlxZQwWh4gHNQXN/3NX1VVquBysbEkb24Pj6/yu66b\nuE3lzC6KCYWsQyW83rGUuqZG/St5H1V4vcMmvN49c5Uz4qYtZHDAQooBkyh9FY++/0d34+2j/fT7\ng3ohr8JGFUFd1wsfoxAOBlXEu7GyS+een1ccv3n68HT36LCGwzKcWdq5l3DOg2GUZ6W49EQ4GOQ5\nmgxqXAIrZo1n7kRd10f7HUb99OFVK0JcVCD1UKfDJ5kBOLPtdtu2bX4uqb+grut8jELeupCGZZQO\nDxfodwjzX+vhcIRF6XzObDkZTbMJJ9dcKtlbZYEbUVXfyArcrvzLdvTVffjQYZ233+/zjekZRt/h\ny48OaRhjupseir0PR7efc3h8aNrSzhoUntbd8lhqfW7Rfv+x+ReQO1pZjCrFucpyEaXyn3les53+\nrEGh7/vDcar5uhO5fOrLiXi42+1WM/aVO3e6qUCMAK5ivumRRxefGnUcpNUrR9EhtTUdXbn61kMl\nACzW3EFhJB+dEeeG5qMdU6ro+z4NWcgHeszTlwEAd2u+q0dWVXU48DCE0DRNygqjoZuHo17TU6Xb\no6mo518OC6JbuWzjrZQTgq/iW3D+ezTfGIVSgUojPPNRnaM2g9J2AGBaawt9YiyHzhwGuPBf6loU\nuCG+ipfv/PfIZaZZOZMFYFG6rqsyZ67EPLo0YFIdWP4iCuc7vFBi87qjR6XrKU61yrWgACH4mQ6z\niKsrpmWJ49WMUn12ovI7oa7rdBWlONr9pi8DkRwuOpkP/y9JV1OKKx0/YQXCI85c6vlW3M+ZcqbV\nXENhNSfCPSh95YbXL8qwf/0KDulyP4dKD51+whuVFgcanfLhyR4avVAnDjn/VVphi0LpVK9dLgDG\nU+W7rovr+scZcPnVH5qmib0Jj2pmGM2uH4Yh9Urkz5y3OsR9nnIyF1NaSPCclyI/taOTDR/tzEBx\nK9Z3RjzTan6Ir+ZEuAelr+JY/5Wu7pj/Go41VLwKYlxH58wWhfx5YgSJO8TbRy9QfKIl44oOSxWy\nqyDFDpfTzzC6FuXI+dXlClsUAFimrutiJR27z0uDGeOP4LRMzpN/E3ddly7A2DTNdruNffxxS3ra\ntNbfkqXSxtEYu90uX2rocOeqqtq2rev6+esIzH31SADuWd4pEMc29n2/P3b5x7yGO/N6wiN5LAjZ\nRSlT/kh9EMtflWd0Wcj9fh9j1tGIkxYiikM7nxmDtCgAMJO8xoqj+mPz+LRzGvOBDnnCaJpms9nE\nh+KEi/hPr/XKgvF842k+h6AAwByGYej7/pyr+h1ufHBaYJK3EMSU0GVGu+XXGFq42ADz4G6XWElC\nUABgDrHyjrP8RxtHyWC08XQFOWTSOg3x2LjwQH6VgPzndVx04WmdGleRx6x8ckfqTMkfSoc8/wSN\nUQBgJrFnfTQKL/Y+hFcVeeqSaNs2zVo8UduNhvVtNpt8Baf8gsPh9asOxQmZyx+dEMV+hPxc0usW\n20ViXBi9bmGKbp21LcdtgXFGVnOJhNWcCPfgwa/io9cCDAfDGKe6/t+ariN4zrmUXt7c+dXl2qrV\nE4tmrOxMOdNq6tfVnAj3wG+25VviZaZn49MJAFMxmBEAKFphiwIAV7e0qyfwZIICABO75y7g9Y3P\n0PUAABQJCgBAkaAAABQJCgBAkaAAABStcNZDaU7OyoahcrfSEo1HWbcRmNYKWxT2BdcuF8zhdIwA\neKwVBgVYJU0FwFWssOsB1up0VtCWAFyCFgUAoEhQAACKBAUAoEhQAACKBAUAoGiFsx4suAQAU1lh\nUBAIAGAquh4AgCJBAQAoEhQAgKJZg0LXdU3TdF135vZhGOL2YRgOtzdNM9rOfaqqb078d+3SAdy2\n+YJCVVV934cQ+r7PJyaUtjdN07ZtCGEYhrZtUyaId+Pttm0P4wV3RRQAuKhqnjkCTdOEEFJlX1XV\ndruNrQV936cyVFW12Wxi3Z/2GR1eVVVd1/F2DA35KVTVTGfEQpwZFO7h0ovppbiHk4XFWl81NNP5\n5LX+aHuq9UMWCEYBIg8Eo6fKs0VY4zvEaWrHxEsBS7C+amiOroeYA2L7weFYhDw9NE2z2+1C1vaQ\n7zMMQ3qqS5YXAPh/Zh2jEKv5w7EIR/ev6/qcp80bJNI/9FhPPScAWLn5gsJms4lNArFN5nKDEPeP\nd6GSAMCtmy8o5MkgbwYoTXGMfRCn94m76YkAgAuZIyiMpjzkRh0HwzDEHodR3Z+GJhx9KkEBAC5k\nphaFuq5Ti8IwDLvdLt7tui61HIy2h6wRou/7NGShruu0jkLcQVAAgAuZbxZHPmYwn9CYZjqEgy6J\nFAjC69eEzJ9qNOtyffNSOM2cwMRLAUuwvmpovstM7/f7o5Mb06THcDBVsnRIaTsAMK21BZ/1RTlO\n8zM68VLAEqyvGnL1SACgSFAAAIrmG6Mwm9JKiytrC4KS0xfK0jEBPMoKWxQsvwgnuDA38CgrDApw\nnzQVAJewwq4HuFuns4K2BOAJtCgAAEWCAgBQJCgAAEWCAgBQZDAjS2cIHsAVaVFg0aQEgOtaYYuC\nlRnvkCUEAC5khUFBIFglUQDgKnQ9AABFggIAUCQoAABFggIAUCQoAABFggIAUCQoAABFK1xHwYJL\nADCVFQYFgQAApqLrAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgKIVLrhkZUYAmMoK\ng4JAAABTWWFQ4LZU1TfXLgIARcYocE1SAsDCCQrcgP3+42sXAeBOzdT1MAxD13WjLel213XDMDRN\nc7hP13XNK4fPFh+9WKmZjygAsEwztSgMw7Db7Y4+VFVV3/chhL7v8wkLTdO0bRuPbds2BYt4N95u\n23aULQCACVXzzBGIv/vzVoSo67q+71MZqqrabDax7q+qarvdxgPzw6uqqus63o6hIT+FqprpjJhE\nGqOgRWEGXm2YwfqqoZlaFHa73dE+gr7v67pOd1MCiFkhHdJ1Xd4gkVoR4g4aFQDgQuabHjkMQ+pZ\nSE0FIUsD8Xbshhi1PYwaJIxLgOd4cLKJJgcgmXXWw3a73W63dV2nQQahXOvnLQ0npEaIpHq8p58S\nrJFpq0AyU4tC3mETmxaapjkcsjD5vwUk+/3HEgDwWNdZmTFvLYgTIw/3yQclnIgUu91us9lMWThY\nrwf7FCQJYGSOrofDKJBCwKjjYBiGmCFG+8d90moKR0cwAACTmyMoNE2z2+3S3IR8XkM+nSGutZC2\nh2w6Qz45Ih/iMJocAQBMa6bpnnG9hHQ3LZYQXsWIeDtvYMgXVgqvjzzIhx/mEyjCGiewrpuZ/Uvj\nHYFnWl81NOv5pB6Eo9tPPHT+9vW9Q+umWloa7wg80/qqodWdz+reoXVTLS2NdwSeaX3VkKtHAgBF\nggIAUHSddRQuqrTS4sraggBgBisMCgIBAExF1wMAUCQoAABFggIAUCQoAABFggIAUCQoAABFK5we\nyaKkJYEBuEUrDAoWXFoOKQHg1q0wKAgEN8f1hwAWa4VBgQUSBQBulMGMAECRoAAAFAkKAECRoAAA\nFAkKAECRoAAAFAkKAEDRCtdRsDIjAExlhUFBIACAqeh6AACKBAUAoEhQAACKBAUAoEhQAACKVjjr\nAXi+qvrmxKOuGw73Q4sC8GinYwSwJoIC8P9pKgBGVtj1YGVGeI7TWUFbAtybFQYFgQAApqLrAQAo\nEhQAgCJBAQAoukJQaJpmtKXruqZpuq4bbR+GIW4fhuFwe9M0o+0AwLTmDgpN0+x2u7yCr6qq7/sQ\nQt/3+YSFpmnatg0hDMPQtm06JN6Nt9u2PYwXAMBUZg0KwzDsdrt8S6zm9/v9MAxxtkKq+He73Xa7\nHYZhGIa6rtP2tm3ruo7bt9ttDBkAwCXMGhTatt1sNvmWvu/ruk53YwIIr+JC6qToui5PGCk0xB00\nKgDAhcwXFJqm2Ww2h5V6PmQhdkyEEEaDD+I+sRUhHBvlAABcwkwLLsUmgaNjD0u1ft7ScEJqhEhK\nKzOeYI2mJ7NOH8C6zREUhmHo+3673c7wbwW1/oykBIDVm6PrIXY3xDmQ8ed+27apIaE0xTEflHBi\nGuRut9MTsXCuMwRwu+ZoUcgXQmiaJg5gjLX7qOMgTnAI2WCFtD1knRRxHYX0qKBwdaIAwFpV8zfU\nV1W13W7T+MS2bWMZ4u30UFVVafBjVVUpUsQMkeZS9n2fn0JVXeGM7lbqehAU7oc3HU5bXzV05atH\nNk1T13UafphaGkII2+22bdu0TEK+4FJVVemQ2YY+AMAdukJQGEWtNOkxHEyVjAsxhYPOhdJ2AGBa\na2shWV+bz5Jphb5D3nQ4bX3VkKtHAgBFVx6jcAmlBZdWFvEAYAYrDAoCAQBMRdcDAFAkKAAARYIC\nAFAkKAAARYICAFAkKAAARSucHgnMIC3ReJR1G2E1tCgA0zsdI4AbssKgUBVcu1ywBpoK4N6ssOvB\nyoxwUaezgrYEWJkVtigAAFMRFACAIkEBACgSFACAIkEBACgSFACAIkEBACgSFACAohUuuFRahNFC\nTADwWCsMCgLBVCyxB4CuB46TEgAIggLP5ypBACu2wq4HpiUHANwzLQoAQJGgAAAUCQoAQJGgAAAU\nrXAwowWXAGAqKwwKAgEATEXXAwBQJCgAAEWCAgBQNF9QGIahaZqmabquGz3Udd3R7fGQruuGYTj6\nVKPtAMC0ZgoKTdO0bRtv933fNE16qKqqvu/j9nzCQjpkGIa2bVMmiHfj7bZtD+MFADCVap45AlVV\nbbfbmA9iTR/vdl3X930qQ1VVm80m1v35IenAuL2u63g7PlV+ClU10xmtXrp6pGs98Cg+Ody59VVD\nc7QoxEo9tSLkzQl939d1ne6mBBCzQtqz67rdbpd2S60IcQeNCgBwIXMEhaZp8niVNxKE13ND0zQx\nEIwGH6QWhVHmAAAuatYFl1IO2G63+cajO+ctDSekRoiktDLjCStrJgKAqcw6PXIYhv1+v9ls8sGJ\nk9s/3oVKAgC37grrKMQhBfkshqO75YMSTqSK3W6nJwIALmSOoNB1Xak7YNRxMAxD7HEY1f1paEI+\n/SERFADgQmYazBiyuQn5gMR8OsMwDLvdLu6W/z+8Pjmiruu0jsJocgQAMK2ZpnvG9RLS3bRYQshG\nOIbXGxjyhZXC6+MN8/aJtNZCesiYg0mYDc/T+ORw59ZXDc16PqXJjSkclB46f/v63qFr8XXP06RP\nzmk+V6zV+qqh1Z3P6t6haxEUeJozg0Lw0WKl1lcNuXokMCXVP6zMrAsusSjn//KDRzmdFXzw4Las\nMCiUpmKurC3omXxZA3COFQYFgWBCmpEB7twKgwKPIgoAcILBjABAkaAAABQJCgBAkaAAABQJCgBA\nkaAAABQJCgBA0QrXUbAyIwBMZYVBQSAAgKnoegAAigQFAKBIUAAAigQFAKBIUAAAigQFAKBIUAAA\nila4joIFlwBgKisMCgIBAExlhUEBuAlV9c2JR/f7j2crCXCCMQrAEp2OEcBsBAVgVpoK4LboegDm\ndjoraEuARdGiAAAUaVFYLT/LAHg+LQrrJCUAMAlB4a4ZVgbAaSvserAyY04UAOA5VhgU7jMQAMAl\n6HoAAIoEBQCgSFAAAIrmCwrDMDRN0zRN13Wjh7quO7o9HtJ13TAMR59qtB0AmNZMQaHrurZt4+2+\n7/OJCVVV9X1/uL1pmnjIMAxt26ZMEO/G223bHsYLAGAq1TxzBKqq2mw2qVJPd7uu6/s+lSHfraqq\n7XbbNE0IIf4/ZoWqquq6jrdjaMhPoapmOqOFSwsumR7JzfHp5aatrxqao0UhVur5T/9U0/d9X9f1\n4fa4c8wH8e5ut0u7paeKO2hUAIALmSMoNE0zile73S6FgHQj3o6BYDT4ILUoxO35IQDA5cw962EY\nhjgQYdQqcChvaTghNUIk1eM944QAYM1mDQpxfGJd1xftv9k/3uUKAwA3baagkBoSttvt4VzHo4fk\ngxJOTIPMezEAgGnNdK2H2JBwWN+PNg7DEHsc0mCFtD1knRRxHYX0qKAAq3T6aunmRMA85pjFEScx\nbrfb0fa4YlKa35h2ixX/aKpkihQxQ8RDRrMrwxrnpTyNCWbcrtP5IOfjzQKtrxqao0UhVvBplaQo\nVvxN09R1nYYT1nWdmge2223btnEtppD1PsRejHTIYf4Abtp+//H5WQG4tEUEnxQCDjsRSvMhS9vX\nF+WeRosCK+bjzZKtrxpa3fms7h16Gt+krJiPN0u2vmrI1SMBgCJBAQAomml65JxKKy2urC0IAGaw\nwqAgEADAVHQ9AABFggIAUCQoAABFKxyjcCcsXQfADASFmyQlQHDVKJiFroc180XJPZOnYRJaFG6b\nKMAdctUomNMKg4IFl2D1TkdkMQImtMKgIBAAwFSMUQAAigQFAKBIUAAAigQFAKBIUAAAigQFAKBo\nhdMjASJrPMPzaVEA7pR1meAcKwwKVcG1ywXMRFMBTGiFXQ9WZgSs8QxTWWGLAgAwFUEBACgSFACA\nIkEBACha4WDGdTDYCoAl0KKwRFICAAshKNwwk8UBuDRdD4smCgBwXSsMCqVFGC3EBACPtcKgIBAA\nwFRWGBQAzvTgwGHdf2AwI0CRKUggKAB3RzsBnG/uoNB13dGNTdMcPjQMQ9w+DMPh9qZpRtsBzrTf\nf3z6v2sXEJZi7qDQ9/2odq+qqu/7+FA+YaFpmrZtQwjDMLRtm46Kd+PtthijbXQAAAkkSURBVG2P\nJg8AYBIzBYXYBnA4cTFW8/v9fhiGOFshVfy73W673Q7DMAxDXddpe9u2dV3H7dvtNoYMAOAS5mtR\naJpms9mMNvZ9X9d1uhsTQHgVF5qmidu7rtvtdmm3FBriDhoVAOBCZpoeGYcUhBAOGwBSGoi34w6j\n7om4T9qYHwIAXM7111Eo1fp5S8MJqREiKa3MeII1mgDgqOsHhcndRK1vcjbcitN/reZHsHrXX0eh\nNMUxH5RwYhrkbre7uZ4IKQFWw58zq3floDDqOIgTHMJBf0TcJw10ODqCYX38UoEr8gcIUTVzQ31V\nVdvtNlXtcVGEWIZ4Oz1aVdVms4kzGqqqSpGiaZrdbpfmUvZ9n59CVc19Rk+QfoL4JoLb5Q+Zo26i\nGnqUK49RaJqmrus0/LCu65Qhtttt27ZplkS+4FJVVemQ7XY7Z4EB4K7MHRQOc1ZcOineHk2VjAsx\nhYPOhdJ2AGBaa2shuYk2Hy2WsAJnDmP0Z35vbqIaepTrz3oAWDHTIrh1KwwKVcG1ywWsiqYC7oQF\nlwCe6HRW0JbAOqywRQEAmIqgAAAUrbDrAWBRXC2Cm6ZFAeCaDGVg4QQFgIvQVMA6rG1diIWsdGEl\nFuBBvihWaSHV0IS0KExPQyIwIV8pXJegcDV+JcCd8yXATVjhrIfSIozztwX5FgBOs2QTy7fCoLCy\nziEAuKIVBgWAlbESA1dkjALAbdNDwUUJCgALpamAJdD1ALBcRjtydYICwM0ziIHLERQAVu6chgdh\nghJjFABu1YS1u14MSrQoANywB7PC+QlA/wVHrTAoLGdlRoCrmypJVNU3ssJ9WmFQEAgAzrfff3x+\nVjj9PBOViGVZYVAA4FEmmYSpyWGtBAUATpmqyeGcf+g5h3MhggIAD5hn3SdtEsskKDyFeUQAyflN\nDtyiamVD/6rq4md0/t+DaAxwjqlyxhK+dWeohmZmwaVLWcLnFeCuaNi4BF0PT/ecKHATkfMmChmU\nc2o3Uc6bKGRQzrOd13nxzyH82xyl4XW38SE+35xdD4LCQijntG6inDdRyKCckzpdyDPbEmZo672J\nF/NRVtiiYGVGAI7SN/EEKxyjsC+4drkAuBTDwi5nhS0KANyhCa+PRU5QAOAuzNPqUFX/PMO/Mqeb\n7HoYhqFpmqZphmF47LGlEQzn7xBCCOGBz8F5T/KsZ5hhh3PcRDkn+Sduopze9Gn/iZsopzd92n9i\nntfzttxei8IwDG3b1nUdQmjbdrPZdF037T+heQoAottrUYgpYRiGYRg2m03f9xM+uYgAALnbCwoh\nhNSEEG9M3qJwDiNsAbgHN9b1EAclNE1zYp8HWwXOaTZ46FJpaxurAgBH3VhQOOpgSOODtfjDOzwY\nBZYwtGcJOyykGM50wh0WUowl7LCQYjjTOXeY6knWZA1B4XUTrASuWwEAojUEhbwnQh0PABO6scGM\nMROM+hpOD1kAAJ7sxoJCCKGu69GsB0EBAC7k9oLCMAy73a6qqqqq+r7fbrdp+5OXa5zH0WmcXdc1\nTXOVGZ4j6QU8LMxyChlup5zJYZBdTjnTi5nkjy6nnKFcmIUU8vCVHJVqIeWMbuKP6MRX+hIK+ajv\n83guXdcttnp6QOlaiwu33W63221+N4RQ13VcsXGz2VytZGUhhLzM+1cXtIxlvu57sdls8hcwL8xy\nCrnf72MxUjnruk4PLaqcSSxP/r4vqpz5+x6lhxZVzlJhllPI7XZbvy7/fC6nnPuDP/ZlvunLL+T5\n3+ejb63RUTdhKd+nz5R/kuIn7KrFeU38BokfnfwjMirndfPN6F9PdxdVyP3rf5wxHca7SytnlJq7\nUpmXVs5ROEgWVc7DBLPkNz2Kb328vbRy5gVY7B/R4V96vH31Qj7h+zzfs/QXt3ALqlCf4/BH26K+\nMjabTfwYjQqZf2Ku+AHK/w5HhVlOIffHyple0kWVM4mfw/x9X1o5S38piypn6UfYogo5MqonFlXO\nwy+ipf0RHf7YS6/n1Qv52O/z0bkcfondhNsbo3DonOUaryh2TR3t0MrL3DTNbrebrVSjYuxftZtF\nu90ulW0hhQwH5YwFW2A5UxmOXrFsaeUchqF6Je9AXUg501/30d7fhRRy5HBp+UWVs67rtm3j5XIW\n/keUSx/O6xbysd/nR+fo3dxIhTUEhaNu4p1YYLiJ1UbIvuYWWMimaaqq2u12qW0/LKycXdftdrsH\nv00WIo74ifVH2riocqYQ0/d9vijeogqZ5IOso0WVM76Sbdu2bbvMP6JR0jqRDhelVLDUVXG7VhsU\neIKmaeLFOUcNDEszDMN+v99sNvGH0bWLMzYMw2FVsUz7/X40tnyZ38KbzSb+Ao6fzCWMyS9Z/pzt\nqqryEV3L/COKVwZOs9uuXZx7t9qgsOQ/1GQ5f5+pIWG73Y5KtZxCjsRv5FS85ZQz/R6KLR8hhLZt\n0wdyOeU8lP/0WVQ582RQ1/UC3/Sk7/vYLZ1bTjlHfzWj3+6LKuc+m90Wsq/05RRypFSwvHNksYU/\nbQ1B4Wivz/KDQv59F0IYhuGKLVSxISH1WSaLKmTXdaVrsSytnJvNJv5MT7O84gu7qHIevt3pG205\n5TzRp7ucQuZlCAcNHgssZy4VZlHlTA0zo/ETiypkrlSw0Z/YktvtTrnCAMoLyMe+Lm16ZBJeHyWb\nD3/NJynNL/3rI4sqZCrA0ZldiyrnSChM9FpCOW/i9cz/uhdbyOjol8/SynkTb3peyFjj5gXLb1+r\nkOd/n4fyFJhbscQK9Wny9HP174ujDguWZ+ErfnoOW0rz8iykkNGoqPnUvkWVMzd63xdVzlt5PW+i\nkPtyNbCoco5Gzyzz9RwVcoF/QeHs7/PRucxczklU+2UPW3uUW2zVOTrnZ2mWVsjSG720cpYsrZw3\n8XreRCFPWFo5b+L1vIlC5k4U7Barp2RVQQEAmNYaBjMCABciKAAARYICAFAkKAAARYICAFAkKAAA\nRYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYIC\nAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARYICAFAkKAAARf8X4MdU\nEId0tpgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "canvas = ROOT.TCanvas()\n",
    "histo_pt.Draw()\n",
    "canvas.Draw()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* Look at the ID and kinematic variables of the four leptons, for signal and background\n",
    "* Compute the two Z 4-vectors and compare them to the variables stored in the tree\n",
    "* Compute the ZZ (aka H) 4-vector and compare it to the variables stored in the tree"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Reconstructed objects"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* Compare generation-level with reconstruction-level lepton variables\n",
    "* Compare the H mass at generation- and reconstruction-levels"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Comparison between data and simulation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "* Apply the correct normalization to the simulation distributions\n",
    "* Stack the background and signal distributions\n",
    "* Compare real data and simulation"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Measure the Higgs mass"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* Measure the reconstructed masses for different generated masses\n",
    "* Measure the mass in data"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}