DimuonMass_NanoAOD_Demo.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Set python for pyspark\n",
    "_this is fixed in the next lcg release_"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import os\n",
    "os.environ[\"PYSPARK_PYTHON\"]=\"/cvmfs/sft.cern.ch/lcg/views/LCG_93/x86_64-centos7-gcc62-opt/bin/python\""
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Setting up Spark\n",
    "\n",
    "Start a pySpark session including Diana-Hep spark-root and histogrammar APIs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Once deleted, variables cannot be recovered. Proceed (y/[n])? n\n",
      "Nothing done.\n",
      "SparkSQL sesssion created\n"
     ]
    }
   ],
   "source": [
    "%reset\n",
    "import pyspark.sql\n",
    "session = pyspark.sql.SparkSession.builder \\\n",
    "    .master('spark://10.64.22.198:7077') \\\n",
    "    .appName('Demo') \\\n",
    "    .config('spark.jars.packages','org.diana-hep:spark-root_2.11:0.1.16,org.diana-hep:histogrammar-sparksql_2.11:1.0.4') \\\n",
    "    .config('spark.driver.extraClassPath','/opt/hadoop/share/hadoop/common/lib/EOSfs.jar') \\\n",
    "    .config('spark.executor.extraClassPath','/opt/hadoop/share/hadoop/common/lib/EOSfs.jar') \\\n",
    "    .config('py-files','helper.py') \\\n",
    "    .getOrCreate()\n",
    "    \n",
    "sqlContext = session\n",
    "\n",
    "print 'SparkSQL sesssion created'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Loading root files stored remotely on EOS via xrootd\n",
    "\n",
    "Loading root files (NanoAOD CMS format) from CERN public EOS area via xrootd.\n",
    "Trees are read using the spark-root"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Loading SingleMuon sample from EOS file\n",
      "Loading DYJetsToLL sample from EOS file\n",
      "Loading WZ sample from EOS file\n"
     ]
    }
   ],
   "source": [
    "from pyspark.sql.functions import lit\n",
    "from samples import *\n",
    "\n",
    "DFList = [] \n",
    "\n",
    "for s in samples:\n",
    "    print 'Loading {0} sample from EOS file'.format(s) \n",
    "    dsPath = \"root://eospublic.cern.ch//eos/opstest/cmspd-bigdata/\"+samples[s]['filename']    \n",
    "    tempDF = sqlContext.read \\\n",
    "                .format(\"org.dianahep.sparkroot\") \\\n",
    "                .option(\"tree\", \"Events\") \\\n",
    "                .load(dsPath)\\\n",
    "                .withColumn(\"sample\", lit(s))                \n",
    "    DFList.append(tempDF)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Access DataFrame content\n",
    "\n",
    "Get list of columns in the DataFrame (\"branches\" of the equivalent ROOT TTree)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "root\n",
      " |-- run: integer (nullable = true)\n",
      " |-- luminosityBlock: integer (nullable = true)\n",
      " |-- event: long (nullable = true)\n",
      " |-- CaloMET_phi: float (nullable = true)\n",
      " |-- CaloMET_pt: float (nullable = true)\n",
      " |-- CaloMET_sumEt: float (nullable = true)\n",
      " |-- nElectron: integer (nullable = true)\n",
      " |-- Electron_deltaEtaSC: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dr03EcalRecHitSumEt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dr03HcalDepth1TowerSumEt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dr03TkSumPt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dxy: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dxyErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dxyErr_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dxy_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dz: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_dzErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_eCorr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_eInvMinusPInv: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_energyErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_hoe: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_ip3d: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_ip3d_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_miniPFRelIso_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_miniPFRelIso_chg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_mvaSpring16GP: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_mvaSpring16HZZ: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_pfRelIso03_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_pfRelIso03_chg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_r9: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_sieie: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_sip3d: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_sip3d_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_mvaTTH: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_cutBased: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_cutBased_HLTPreSel: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_jetIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_pdgId: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_photonIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_tightCharge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_vidNestedWPBitmap: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Electron_convVeto: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Electron_cutBased_HEEP: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Electron_isPFcand: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Electron_lostHits: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Electron_mvaSpring16GP_WP80: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Electron_mvaSpring16GP_WP90: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Electron_mvaSpring16HZZ_WPL: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- nFatJet: integer (nullable = true)\n",
      " |-- FatJet_area: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_btagCMVA: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_btagCSVV2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_btagDeepB: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_btagHbb: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_msoftdrop: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_n2b1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_n3b1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_tau1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_tau2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_tau3: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_tau4: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- FatJet_subJetIdx1: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- FatJet_subJetIdx2: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- nJet: integer (nullable = true)\n",
      " |-- Jet_area: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_btagCMVA: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_btagCSVV2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_btagDeepB: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_btagDeepC: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_chEmEF: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_chHEF: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_neEmEF: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_neHEF: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_qgl: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_rawFactor: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_bReg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Jet_electronIdx1: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_electronIdx2: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_jetId: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_muonIdx1: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_muonIdx2: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_nConstituents: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_nElectrons: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_nMuons: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Jet_puId: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- MET_MetUnclustEnUpDeltaX: float (nullable = true)\n",
      " |-- MET_MetUnclustEnUpDeltaY: float (nullable = true)\n",
      " |-- MET_covXX: float (nullable = true)\n",
      " |-- MET_covXY: float (nullable = true)\n",
      " |-- MET_covYY: float (nullable = true)\n",
      " |-- MET_phi: float (nullable = true)\n",
      " |-- MET_pt: float (nullable = true)\n",
      " |-- MET_significance: float (nullable = true)\n",
      " |-- MET_sumEt: float (nullable = true)\n",
      " |-- nMuon: integer (nullable = true)\n",
      " |-- Muon_dxy: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_dxyErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_dxyErr_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_dxy_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_dz: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_dzErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_ip3d: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_ip3d_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_miniPFRelIso_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_miniPFRelIso_chg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_pfRelIso03_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_pfRelIso03_chg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_pfRelIso04_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_ptErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_segmentComp: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_sip3d: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_sip3d_ugo: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_mvaTTH: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_jetIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_nStations: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_nTrackerLayers: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_pdgId: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_tightCharge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_highPtId: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Muon_isPFcand: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Muon_mediumId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Muon_softId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Muon_tightId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- nPhoton: integer (nullable = true)\n",
      " |-- Photon_eCorr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_energyErr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_hoe: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_mvaID: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_pfRelIso03_all: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_pfRelIso03_chg: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_r9: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_sieie: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Photon_charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_cutBased: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_electronIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_jetIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_pdgId: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_vidNestedWPBitmap: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Photon_electronVeto: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Photon_mvaID_WP80: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Photon_mvaID_WP90: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Photon_pixelSeed: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- PuppiMET_phi: float (nullable = true)\n",
      " |-- PuppiMET_pt: float (nullable = true)\n",
      " |-- PuppiMET_sumEt: float (nullable = true)\n",
      " |-- RawMET_phi: float (nullable = true)\n",
      " |-- RawMET_pt: float (nullable = true)\n",
      " |-- RawMET_sumEt: float (nullable = true)\n",
      " |-- fixedGridRhoFastjetAll: float (nullable = true)\n",
      " |-- fixedGridRhoFastjetCentralCalo: float (nullable = true)\n",
      " |-- fixedGridRhoFastjetCentralNeutral: float (nullable = true)\n",
      " |-- nSoftActivityJet: integer (nullable = true)\n",
      " |-- SoftActivityJet_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SoftActivityJet_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SoftActivityJet_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SoftActivityJetHT: float (nullable = true)\n",
      " |-- SoftActivityJetHT10: float (nullable = true)\n",
      " |-- SoftActivityJetHT2: float (nullable = true)\n",
      " |-- SoftActivityJetHT5: float (nullable = true)\n",
      " |-- SoftActivityJetNjets10: integer (nullable = true)\n",
      " |-- SoftActivityJetNjets2: integer (nullable = true)\n",
      " |-- SoftActivityJetNjets5: integer (nullable = true)\n",
      " |-- nSubJet: integer (nullable = true)\n",
      " |-- SubJet_btagCMVA: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_btagCSVV2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_btagDeepB: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_n2b1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_n3b1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_tau1: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_tau2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_tau3: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SubJet_tau4: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- nTau: integer (nullable = true)\n",
      " |-- Tau_chargedIso: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_dxy: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_dz: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_footprintCorr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_leadTkDeltaEta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_leadTkDeltaPhi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_leadTkPtOverTauPt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_neutralIso: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_photonsOutsideSignalCone: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_puCorr: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_rawAntiEle: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_rawIso: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_rawMVAnewDM: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_rawMVAoldDM: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_rawMVAoldDMdR03: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Tau_charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Tau_decayMode: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Tau_jetIdx: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Tau_rawAntiEleCat: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Tau_idAntiEle: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Tau_idAntiMu: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Tau_idDecayMode: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Tau_idDecayModeNewDMs: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Tau_idMVAnewDM: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Tau_idMVAoldDM: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Tau_idMVAoldDMdR03: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- TkMET_phi: float (nullable = true)\n",
      " |-- TkMET_pt: float (nullable = true)\n",
      " |-- TkMET_sumEt: float (nullable = true)\n",
      " |-- nTrigObj: integer (nullable = true)\n",
      " |-- TrigObj_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_l1pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_l1pt_2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_l2pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- TrigObj_id: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- TrigObj_l1iso: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- TrigObj_l1charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- TrigObj_filterBits: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- nOtherPV: integer (nullable = true)\n",
      " |-- OtherPV_z: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- PV_ndof: float (nullable = true)\n",
      " |-- PV_x: float (nullable = true)\n",
      " |-- PV_y: float (nullable = true)\n",
      " |-- PV_z: float (nullable = true)\n",
      " |-- PV_chi2: float (nullable = true)\n",
      " |-- PV_score: float (nullable = true)\n",
      " |-- PV_npvs: integer (nullable = true)\n",
      " |-- nSV: integer (nullable = true)\n",
      " |-- SV_dlen: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_dlenSig: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_pAngle: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Electron_cleanmask: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Jet_cleanmask: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Muon_cleanmask: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Photon_cleanmask: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- Tau_cleanmask: array (nullable = true)\n",
      " |    |-- element: byte (containsNull = true)\n",
      " |-- SV_chi2: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_ndof: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_x: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_y: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- SV_z: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- HLTriggerFirstPath: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet360_TrimMass30: boolean (nullable = true)\n",
      " |-- HLT_AK8PFHT700_TrimR0p1PT0p03Mass50: boolean (nullable = true)\n",
      " |-- HLT_AK8PFHT650_TrimR0p1PT0p03Mass50: boolean (nullable = true)\n",
      " |-- HLT_AK8PFHT600_TrimR0p1PT0p03Mass50_BTagCSV_p20: boolean (nullable = true)\n",
      " |-- HLT_CaloJet500_NoJetID: boolean (nullable = true)\n",
      " |-- HLT_Dimuon13_PsiPrime: boolean (nullable = true)\n",
      " |-- HLT_Dimuon13_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Dimuon20_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle24_22_eta2p1_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle33_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle33_CaloIdL_MW: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle33_CaloIdL_GsfTrkIdVL_MW: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle33_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle37_Ele27_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_DoubleMediumIsoPFTau32_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_DoubleMediumIsoPFTau35_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_DoubleMediumIsoPFTau40_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu33NoFiltersNoVtx: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu38NoFiltersNoVtx: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu23NoFiltersNoVtxDisplaced: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu28NoFiltersNoVtxDisplaced: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu4_3_Bs: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu4_3_Jpsi_Displaced: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu4_JpsiTrk_Displaced: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu4_LowMassNonResonantTrk_Displaced: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu3_Trk_Tau3mu: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu4_PsiPrimeTrk_Displaced: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_L2Mu2_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_L2Mu2_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track2_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track3p5_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track7_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track2_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track3p5_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Mu7p5_Track7_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Dimuon0er16_Jpsi_NoOS_NoVertexing: boolean (nullable = true)\n",
      " |-- HLT_Dimuon0er16_Jpsi_NoVertexing: boolean (nullable = true)\n",
      " |-- HLT_Dimuon6_Jpsi_NoVertexing: boolean (nullable = true)\n",
      " |-- HLT_DoublePhoton60: boolean (nullable = true)\n",
      " |-- HLT_DoublePhoton85: boolean (nullable = true)\n",
      " |-- HLT_Ele22_eta2p1_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele22_eta2p1_WPLoose_Gsf_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_Ele23_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele23_WPLoose_Gsf_WHbbBoost: boolean (nullable = true)\n",
      " |-- HLT_Ele24_eta2p1_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele24_eta2p1_WPLoose_Gsf_LooseIsoPFTau20: boolean (nullable = true)\n",
      " |-- HLT_Ele24_eta2p1_WPLoose_Gsf_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_Ele25_WPTight_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele25_eta2p1_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele25_eta2p1_WPTight_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele27_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele27_WPLoose_Gsf_WHbbBoost: boolean (nullable = true)\n",
      " |-- HLT_Ele27_WPTight_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf_DoubleMediumIsoPFTau32_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf_DoubleMediumIsoPFTau35_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf_DoubleMediumIsoPFTau40_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPTight_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele32_eta2p1_WPLoose_Gsf_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_Ele32_eta2p1_WPTight_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele35_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele35_CaloIdVT_GsfTrkIdT_PFJet150_PFJet50: boolean (nullable = true)\n",
      " |-- HLT_Ele45_WPLoose_Gsf: boolean (nullable = true)\n",
      " |-- HLT_Ele45_CaloIdVT_GsfTrkIdT_PFJet200_PFJet50: boolean (nullable = true)\n",
      " |-- HLT_Ele105_CaloIdVT_GsfTrkIdT: boolean (nullable = true)\n",
      " |-- HLT_Ele30WP60_SC4_Mass55: boolean (nullable = true)\n",
      " |-- HLT_Ele30WP60_Ele8_Mass55: boolean (nullable = true)\n",
      " |-- HLT_HT200: boolean (nullable = true)\n",
      " |-- HLT_HT275: boolean (nullable = true)\n",
      " |-- HLT_HT325: boolean (nullable = true)\n",
      " |-- HLT_HT425: boolean (nullable = true)\n",
      " |-- HLT_HT575: boolean (nullable = true)\n",
      " |-- HLT_HT410to430: boolean (nullable = true)\n",
      " |-- HLT_HT430to450: boolean (nullable = true)\n",
      " |-- HLT_HT450to470: boolean (nullable = true)\n",
      " |-- HLT_HT470to500: boolean (nullable = true)\n",
      " |-- HLT_HT500to550: boolean (nullable = true)\n",
      " |-- HLT_HT550to650: boolean (nullable = true)\n",
      " |-- HLT_HT650: boolean (nullable = true)\n",
      " |-- HLT_Mu16_eta2p1_MET30: boolean (nullable = true)\n",
      " |-- HLT_IsoMu16_eta2p1_MET30: boolean (nullable = true)\n",
      " |-- HLT_IsoMu16_eta2p1_MET30_LooseIsoPFTau50_Trk30_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_IsoMu17_eta2p1_LooseIsoPFTau20: boolean (nullable = true)\n",
      " |-- HLT_IsoMu17_eta2p1_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_DoubleIsoMu17_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_DoubleIsoMu17_eta2p1_noDzCut: boolean (nullable = true)\n",
      " |-- HLT_IsoMu18: boolean (nullable = true)\n",
      " |-- HLT_IsoMu19_eta2p1_LooseIsoPFTau20: boolean (nullable = true)\n",
      " |-- HLT_IsoMu19_eta2p1_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_IsoMu19_eta2p1_MediumIsoPFTau32_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_IsoMu20: boolean (nullable = true)\n",
      " |-- HLT_IsoMu21_eta2p1_LooseIsoPFTau20_SingleL1: boolean (nullable = true)\n",
      " |-- HLT_IsoMu21_eta2p1_MediumIsoPFTau32_Trk1_eta2p1_Reg: boolean (nullable = true)\n",
      " |-- HLT_IsoMu22: boolean (nullable = true)\n",
      " |-- HLT_IsoMu22_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_IsoMu24: boolean (nullable = true)\n",
      " |-- HLT_IsoMu27: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu18: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu20: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu22: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu22_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu24: boolean (nullable = true)\n",
      " |-- HLT_IsoTkMu27: boolean (nullable = true)\n",
      " |-- HLT_JetE30_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_JetE30_NoBPTX: boolean (nullable = true)\n",
      " |-- HLT_JetE50_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_JetE70_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_L1SingleMu18: boolean (nullable = true)\n",
      " |-- HLT_L2Mu10: boolean (nullable = true)\n",
      " |-- HLT_L1SingleMuOpen: boolean (nullable = true)\n",
      " |-- HLT_L1SingleMuOpen_DT: boolean (nullable = true)\n",
      " |-- HLT_L2DoubleMu23_NoVertex: boolean (nullable = true)\n",
      " |-- HLT_L2DoubleMu28_NoVertex_2Cha_Angle2p5_Mass10: boolean (nullable = true)\n",
      " |-- HLT_L2DoubleMu38_NoVertex_2Cha_Angle2p5_Mass10: boolean (nullable = true)\n",
      " |-- HLT_L2Mu10_NoVertex_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_L2Mu10_NoVertex_NoBPTX: boolean (nullable = true)\n",
      " |-- HLT_L2Mu35_NoVertex_3Sta_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_L2Mu40_NoVertex_3Sta_NoBPTX3BX: boolean (nullable = true)\n",
      " |-- HLT_LooseIsoPFTau50_Trk30_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_LooseIsoPFTau50_Trk30_eta2p1_MET80: boolean (nullable = true)\n",
      " |-- HLT_LooseIsoPFTau50_Trk30_eta2p1_MET90: boolean (nullable = true)\n",
      " |-- HLT_LooseIsoPFTau50_Trk30_eta2p1_MET110: boolean (nullable = true)\n",
      " |-- HLT_LooseIsoPFTau50_Trk30_eta2p1_MET120: boolean (nullable = true)\n",
      " |-- HLT_PFTau120_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_VLooseIsoPFTau120_Trk50_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_VLooseIsoPFTau140_Trk50_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Mu8: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Mu8_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Mu8_SameSign: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Mu8_SameSign_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu20_Mu10: boolean (nullable = true)\n",
      " |-- HLT_Mu20_Mu10_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu20_Mu10_SameSign: boolean (nullable = true)\n",
      " |-- HLT_Mu20_Mu10_SameSign_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TkMu8_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL_Mu8_TrkIsoVVL_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL_TkMu8_TrkIsoVVL: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL_TkMu8_TrkIsoVVL_DZ: boolean (nullable = true)\n",
      " |-- HLT_Mu25_TkMu0_dEta18_Onia: boolean (nullable = true)\n",
      " |-- HLT_Mu27_TkMu8: boolean (nullable = true)\n",
      " |-- HLT_Mu30_TkMu11: boolean (nullable = true)\n",
      " |-- HLT_Mu30_eta2p1_PFJet150_PFJet50: boolean (nullable = true)\n",
      " |-- HLT_Mu40_TkMu11: boolean (nullable = true)\n",
      " |-- HLT_Mu40_eta2p1_PFJet200_PFJet50: boolean (nullable = true)\n",
      " |-- HLT_Mu20: boolean (nullable = true)\n",
      " |-- HLT_TkMu20: boolean (nullable = true)\n",
      " |-- HLT_Mu24_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_TkMu24_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_Mu27: boolean (nullable = true)\n",
      " |-- HLT_TkMu27: boolean (nullable = true)\n",
      " |-- HLT_Mu45_eta2p1: boolean (nullable = true)\n",
      " |-- HLT_Mu50: boolean (nullable = true)\n",
      " |-- HLT_TkMu50: boolean (nullable = true)\n",
      " |-- HLT_Mu38NoFiltersNoVtx_Photon38_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu42NoFiltersNoVtx_Photon42_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu28NoFiltersNoVtxDisplaced_Photon28_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu33NoFiltersNoVtxDisplaced_Photon33_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu23NoFiltersNoVtx_Photon23_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu18NoFiltersNoVtx: boolean (nullable = true)\n",
      " |-- HLT_Mu33NoFiltersNoVtxDisplaced_DisplacedJet50_Tight: boolean (nullable = true)\n",
      " |-- HLT_Mu33NoFiltersNoVtxDisplaced_DisplacedJet50_Loose: boolean (nullable = true)\n",
      " |-- HLT_Mu28NoFiltersNoVtx_DisplacedJet40_Loose: boolean (nullable = true)\n",
      " |-- HLT_Mu38NoFiltersNoVtxDisplaced_DisplacedJet60_Tight: boolean (nullable = true)\n",
      " |-- HLT_Mu38NoFiltersNoVtxDisplaced_DisplacedJet60_Loose: boolean (nullable = true)\n",
      " |-- HLT_Mu38NoFiltersNoVtx_DisplacedJet60_Loose: boolean (nullable = true)\n",
      " |-- HLT_Mu28NoFiltersNoVtx_CentralCaloJet40: boolean (nullable = true)\n",
      " |-- HLT_PFHT300_PFMET100: boolean (nullable = true)\n",
      " |-- HLT_PFHT300_PFMET110: boolean (nullable = true)\n",
      " |-- HLT_PFHT550_4JetPt50: boolean (nullable = true)\n",
      " |-- HLT_PFHT650_4JetPt50: boolean (nullable = true)\n",
      " |-- HLT_PFHT750_4JetPt50: boolean (nullable = true)\n",
      " |-- HLT_PFJet15_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_PFJet25_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet15_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet25_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet15_FBEta3_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet25_FBEta3_NoCaloMatched: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve15_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve25_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve35_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet40: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet60: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet80: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet140: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet200: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet260: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet320: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet400: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet450: boolean (nullable = true)\n",
      " |-- HLT_AK8PFJet500: boolean (nullable = true)\n",
      " |-- HLT_PFJet40: boolean (nullable = true)\n",
      " |-- HLT_PFJet60: boolean (nullable = true)\n",
      " |-- HLT_PFJet80: boolean (nullable = true)\n",
      " |-- HLT_PFJet140: boolean (nullable = true)\n",
      " |-- HLT_PFJet200: boolean (nullable = true)\n",
      " |-- HLT_PFJet260: boolean (nullable = true)\n",
      " |-- HLT_PFJet320: boolean (nullable = true)\n",
      " |-- HLT_PFJet400: boolean (nullable = true)\n",
      " |-- HLT_PFJet450: boolean (nullable = true)\n",
      " |-- HLT_PFJet500: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve40: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve60: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve80: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve140: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve200: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve260: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve320: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve400: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve500: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve60_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve80_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve100_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve160_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve220_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJetAve300_HFJEC: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet40_DEta3p5_MJJ600_PFMETNoMu140: boolean (nullable = true)\n",
      " |-- HLT_DiPFJet40_DEta3p5_MJJ600_PFMETNoMu80: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet55_PFMET110: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet170: boolean (nullable = true)\n",
      " |-- HLT_SingleCentralPFJet170_CFMax0p1: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet170_CFMax0p1: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet220_CFMax0p3: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet330_CFMax0p5: boolean (nullable = true)\n",
      " |-- HLT_DiCentralPFJet430: boolean (nullable = true)\n",
      " |-- HLT_PFHT125: boolean (nullable = true)\n",
      " |-- HLT_PFHT200: boolean (nullable = true)\n",
      " |-- HLT_PFHT250: boolean (nullable = true)\n",
      " |-- HLT_PFHT300: boolean (nullable = true)\n",
      " |-- HLT_PFHT350: boolean (nullable = true)\n",
      " |-- HLT_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_PFHT475: boolean (nullable = true)\n",
      " |-- HLT_PFHT600: boolean (nullable = true)\n",
      " |-- HLT_PFHT650: boolean (nullable = true)\n",
      " |-- HLT_PFHT800: boolean (nullable = true)\n",
      " |-- HLT_PFHT900: boolean (nullable = true)\n",
      " |-- HLT_PFHT200_PFAlphaT0p51: boolean (nullable = true)\n",
      " |-- HLT_PFHT200_DiPFJetAve90_PFAlphaT0p57: boolean (nullable = true)\n",
      " |-- HLT_PFHT200_DiPFJetAve90_PFAlphaT0p63: boolean (nullable = true)\n",
      " |-- HLT_PFHT250_DiPFJetAve90_PFAlphaT0p55: boolean (nullable = true)\n",
      " |-- HLT_PFHT250_DiPFJetAve90_PFAlphaT0p58: boolean (nullable = true)\n",
      " |-- HLT_PFHT300_DiPFJetAve90_PFAlphaT0p53: boolean (nullable = true)\n",
      " |-- HLT_PFHT300_DiPFJetAve90_PFAlphaT0p54: boolean (nullable = true)\n",
      " |-- HLT_PFHT350_DiPFJetAve90_PFAlphaT0p52: boolean (nullable = true)\n",
      " |-- HLT_PFHT350_DiPFJetAve90_PFAlphaT0p53: boolean (nullable = true)\n",
      " |-- HLT_PFHT400_DiPFJetAve90_PFAlphaT0p51: boolean (nullable = true)\n",
      " |-- HLT_PFHT400_DiPFJetAve90_PFAlphaT0p52: boolean (nullable = true)\n",
      " |-- HLT_MET60_IsoTrk35_Loose: boolean (nullable = true)\n",
      " |-- HLT_MET75_IsoTrk50: boolean (nullable = true)\n",
      " |-- HLT_MET90_IsoTrk50: boolean (nullable = true)\n",
      " |-- HLT_PFMET120_BTagCSV_p067: boolean (nullable = true)\n",
      " |-- HLT_PFMET120_Mu5: boolean (nullable = true)\n",
      " |-- HLT_PFMET170_NoiseCleaned: boolean (nullable = true)\n",
      " |-- HLT_PFMET170_HBHECleaned: boolean (nullable = true)\n",
      " |-- HLT_PFMET170_JetIdCleaned: boolean (nullable = true)\n",
      " |-- HLT_PFMET170_NotCleaned: boolean (nullable = true)\n",
      " |-- HLT_PFMET170_BeamHaloCleaned: boolean (nullable = true)\n",
      " |-- HLT_PFMET90_PFMHT90_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMET100_PFMHT100_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMET110_PFMHT110_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMET120_PFMHT120_IDTight: boolean (nullable = true)\n",
      " |-- HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDTight_BTagCSV_p067: boolean (nullable = true)\n",
      " |-- HLT_CaloMHTNoPU90_PFMET90_PFMHT90_IDTight: boolean (nullable = true)\n",
      " |-- HLT_QuadPFJet_BTagCSV_p016_p11_VBF_Mqq200: boolean (nullable = true)\n",
      " |-- HLT_QuadPFJet_BTagCSV_p016_VBF_Mqq460: boolean (nullable = true)\n",
      " |-- HLT_QuadPFJet_BTagCSV_p016_p11_VBF_Mqq240: boolean (nullable = true)\n",
      " |-- HLT_QuadPFJet_BTagCSV_p016_VBF_Mqq500: boolean (nullable = true)\n",
      " |-- HLT_QuadPFJet_VBF: boolean (nullable = true)\n",
      " |-- HLT_L1_TripleJet_VBF: boolean (nullable = true)\n",
      " |-- HLT_QuadJet45_TripleBTagCSV_p087: boolean (nullable = true)\n",
      " |-- HLT_QuadJet45_DoubleBTagCSV_p087: boolean (nullable = true)\n",
      " |-- HLT_DoubleJet90_Double30_TripleBTagCSV_p087: boolean (nullable = true)\n",
      " |-- HLT_DoubleJet90_Double30_DoubleBTagCSV_p087: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_DoubleBTagCSV_p026_DoublePFJetsC160: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_DoubleBTagCSV_p014_DoublePFJetsC100MaxDeta1p6: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC112_DoubleBTagCSV_p026_DoublePFJetsC172: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC112_DoubleBTagCSV_p014_DoublePFJetsC112MaxDeta1p6: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_SingleBTagCSV_p026: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_SingleBTagCSV_p014: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_SingleBTagCSV_p026_SinglePFJetC350: boolean (nullable = true)\n",
      " |-- HLT_DoubleJetsC100_SingleBTagCSV_p014_SinglePFJetC350: boolean (nullable = true)\n",
      " |-- HLT_Photon135_PFMET100: boolean (nullable = true)\n",
      " |-- HLT_Photon22_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon22_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Photon250_NoHE: boolean (nullable = true)\n",
      " |-- HLT_Photon300_NoHE: boolean (nullable = true)\n",
      " |-- HLT_Photon26_R9Id85_OR_CaloId24b40e_Iso50T80L_Photon16_AND_HE10_R9Id65_Eta2_Mass60: boolean (nullable = true)\n",
      " |-- HLT_Photon36_R9Id85_OR_CaloId24b40e_Iso50T80L_Photon22_AND_HE10_R9Id65_Eta2_Mass15: boolean (nullable = true)\n",
      " |-- HLT_Photon36_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon36_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Photon50_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon50_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Photon75_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon75_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Photon90_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon90_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Photon120_R9Id90_HE10_Iso40_EBOnly_PFMET40: boolean (nullable = true)\n",
      " |-- HLT_Photon120_R9Id90_HE10_Iso40_EBOnly_VBF: boolean (nullable = true)\n",
      " |-- HLT_Mu8_TrkIsoVVL: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL: boolean (nullable = true)\n",
      " |-- HLT_Ele8_CaloIdL_TrackIdL_IsoVL_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele12_CaloIdL_TrackIdL_IsoVL_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele17_CaloIdL_TrackIdL_IsoVL_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele23_CaloIdL_TrackIdL_IsoVL_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_DiJet20_Mu5: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_DiJet40_Mu5: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_DiJet70_Mu5: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_DiJet110_Mu5: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_Jet300_Mu5: boolean (nullable = true)\n",
      " |-- HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL_DZ: boolean (nullable = true)\n",
      " |-- HLT_Ele17_Ele12_CaloIdL_TrackIdL_IsoVL_DZ: boolean (nullable = true)\n",
      " |-- HLT_Ele16_Ele12_Ele8_CaloIdL_TrackIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu8_TrkIsoVVL_Ele23_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Mu8_TrkIsoVVL_Ele17_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Mu23_TrkIsoVVL_Ele8_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Mu23_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Mu17_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Mu30_Ele30_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_Mu37_Ele27_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_Mu27_Ele37_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_Mu8_DiEle12_CaloIdL_TrackIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu12_Photon25_CaloIdL: boolean (nullable = true)\n",
      " |-- HLT_Mu12_Photon25_CaloIdL_L1ISO: boolean (nullable = true)\n",
      " |-- HLT_Mu12_Photon25_CaloIdL_L1OR: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Photon22_CaloIdL_L1ISO: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Photon30_CaloIdL_L1ISO: boolean (nullable = true)\n",
      " |-- HLT_Mu17_Photon35_CaloIdL_L1ISO: boolean (nullable = true)\n",
      " |-- HLT_DiMu9_Ele9_CaloIdL_TrackIdL: boolean (nullable = true)\n",
      " |-- HLT_TripleMu_5_3_3: boolean (nullable = true)\n",
      " |-- HLT_TripleMu_12_10_5: boolean (nullable = true)\n",
      " |-- HLT_Mu3er_PFHT140_PFMET125: boolean (nullable = true)\n",
      " |-- HLT_Mu6_PFHT200_PFMET80_BTagCSV_p067: boolean (nullable = true)\n",
      " |-- HLT_Mu6_PFHT200_PFMET100: boolean (nullable = true)\n",
      " |-- HLT_Mu14er_PFMET100: boolean (nullable = true)\n",
      " |-- HLT_Ele17_Ele12_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Ele12_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Ele17_CaloIdL_GsfTrkIdVL: boolean (nullable = true)\n",
      " |-- HLT_Ele17_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_Ele23_CaloIdL_TrackIdL_IsoVL: boolean (nullable = true)\n",
      " |-- HLT_AK8DiPFJet280_200_TrimMass30: boolean (nullable = true)\n",
      " |-- HLT_AK8DiPFJet250_200_TrimMass30: boolean (nullable = true)\n",
      " |-- HLT_AK8DiPFJet280_200_TrimMass30_BTagCSV_p20: boolean (nullable = true)\n",
      " |-- HLT_AK8DiPFJet250_200_TrimMass30_BTagCSV_p20: boolean (nullable = true)\n",
      " |-- HLT_PFHT650_WideJetMJJ900DEtaJJ1p5: boolean (nullable = true)\n",
      " |-- HLT_PFHT650_WideJetMJJ950DEtaJJ1p5: boolean (nullable = true)\n",
      " |-- HLT_Photon22: boolean (nullable = true)\n",
      " |-- HLT_Photon30: boolean (nullable = true)\n",
      " |-- HLT_Photon36: boolean (nullable = true)\n",
      " |-- HLT_Photon50: boolean (nullable = true)\n",
      " |-- HLT_Photon75: boolean (nullable = true)\n",
      " |-- HLT_Photon90: boolean (nullable = true)\n",
      " |-- HLT_Photon120: boolean (nullable = true)\n",
      " |-- HLT_Photon175: boolean (nullable = true)\n",
      " |-- HLT_Photon165_HE10: boolean (nullable = true)\n",
      " |-- HLT_Photon22_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon30_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon36_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon50_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon75_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon90_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon120_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Photon165_R9Id90_HE10_IsoM: boolean (nullable = true)\n",
      " |-- HLT_Diphoton30_18_R9Id_OR_IsoCaloId_AND_HE_R9Id_Mass90: boolean (nullable = true)\n",
      " |-- HLT_Diphoton30_18_R9Id_OR_IsoCaloId_AND_HE_R9Id_DoublePixelSeedMatch_Mass70: boolean (nullable = true)\n",
      " |-- HLT_Diphoton30PV_18PV_R9Id_AND_IsoCaloId_AND_HE_R9Id_DoublePixelVeto_Mass55: boolean (nullable = true)\n",
      " |-- HLT_Diphoton30_18_Solid_R9Id_AND_IsoCaloId_AND_HE_R9Id_Mass55: boolean (nullable = true)\n",
      " |-- HLT_Diphoton30EB_18EB_R9Id_OR_IsoCaloId_AND_HE_R9Id_DoublePixelVeto_Mass55: boolean (nullable = true)\n",
      " |-- HLT_Dimuon0_Jpsi_Muon: boolean (nullable = true)\n",
      " |-- HLT_Dimuon0_Upsilon_Muon: boolean (nullable = true)\n",
      " |-- HLT_QuadMuon0_Dimuon0_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_QuadMuon0_Dimuon0_Upsilon: boolean (nullable = true)\n",
      " |-- HLT_Rsq0p25: boolean (nullable = true)\n",
      " |-- HLT_Rsq0p30: boolean (nullable = true)\n",
      " |-- HLT_RsqMR240_Rsq0p09_MR200: boolean (nullable = true)\n",
      " |-- HLT_RsqMR240_Rsq0p09_MR200_4jet: boolean (nullable = true)\n",
      " |-- HLT_RsqMR270_Rsq0p09_MR200: boolean (nullable = true)\n",
      " |-- HLT_RsqMR270_Rsq0p09_MR200_4jet: boolean (nullable = true)\n",
      " |-- HLT_Rsq0p02_MR300_TriPFJet80_60_40_BTagCSV_p063_p20_Mbb60_200: boolean (nullable = true)\n",
      " |-- HLT_Rsq0p02_MR300_TriPFJet80_60_40_DoubleBTagCSV_p063_Mbb60_200: boolean (nullable = true)\n",
      " |-- HLT_HT200_DisplacedDijet40_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_HT250_DisplacedDijet40_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_HT350_DisplacedDijet40_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_HT350_DisplacedDijet80_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_HT350_DisplacedDijet80_Tight_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_HT350_DisplacedDijet40_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_HT400_DisplacedDijet40_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_HT500_DisplacedDijet40_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_HT550_DisplacedDijet40_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_HT650_DisplacedDijet80_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_HT750_DisplacedDijet80_Inclusive: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_DisplacedTrack_2TrackIP2DSig5: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_TightID_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_Hadronic: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_Hadronic_2PromptTrack: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_TightID_Hadronic: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_VTightID_Hadronic: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_VVTightID_Hadronic: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_VTightID_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_VBF_DisplacedJet40_VVTightID_DisplacedTrack: boolean (nullable = true)\n",
      " |-- HLT_PFMETNoMu90_PFMHTNoMu90_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMETNoMu100_PFMHTNoMu100_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMETNoMu110_PFMHTNoMu110_IDTight: boolean (nullable = true)\n",
      " |-- HLT_PFMETNoMu120_PFMHTNoMu120_IDTight: boolean (nullable = true)\n",
      " |-- HLT_MonoCentralPFJet80_PFMETNoMu90_PFMHTNoMu90_IDTight: boolean (nullable = true)\n",
      " |-- HLT_MonoCentralPFJet80_PFMETNoMu100_PFMHTNoMu100_IDTight: boolean (nullable = true)\n",
      " |-- HLT_MonoCentralPFJet80_PFMETNoMu110_PFMHTNoMu110_IDTight: boolean (nullable = true)\n",
      " |-- HLT_MonoCentralPFJet80_PFMETNoMu120_PFMHTNoMu120_IDTight: boolean (nullable = true)\n",
      " |-- HLT_Ele27_eta2p1_WPLoose_Gsf_HT200: boolean (nullable = true)\n",
      " |-- HLT_Photon90_CaloIdL_PFHT500: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu8_Mass8_PFHT250: boolean (nullable = true)\n",
      " |-- HLT_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT250: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT250: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu8_Mass8_PFHT300: boolean (nullable = true)\n",
      " |-- HLT_Mu8_Ele8_CaloIdM_TrackIdM_Mass8_PFHT300: boolean (nullable = true)\n",
      " |-- HLT_DoubleEle8_CaloIdM_TrackIdM_Mass8_PFHT300: boolean (nullable = true)\n",
      " |-- HLT_Mu10_CentralPFJet30_BTagCSV_p13: boolean (nullable = true)\n",
      " |-- HLT_DoubleMu3_PFMET50: boolean (nullable = true)\n",
      " |-- HLT_Ele10_CaloIdM_TrackIdM_CentralPFJet30_BTagCSV_p13: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_BTagCSV_p067_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_PFHT350_PFMET50: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_PFHT600: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_PFHT350: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_PFHT400_PFMET50: boolean (nullable = true)\n",
      " |-- HLT_Ele15_IsoVVVL_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Ele50_IsoVVVL_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Mu8_TrkIsoVVL_DiPFJet40_DEta3p5_MJJ750_HTT300_PFMETNoMu60: boolean (nullable = true)\n",
      " |-- HLT_Mu10_TrkIsoVVL_DiPFJet40_DEta3p5_MJJ750_HTT350_PFMETNoMu60: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_BTagCSV_p067_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_PFHT350_PFMET50: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_PFHT600: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_PFHT350: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_PFHT400_PFMET50: boolean (nullable = true)\n",
      " |-- HLT_Mu15_IsoVVVL_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Mu50_IsoVVVL_PFHT400: boolean (nullable = true)\n",
      " |-- HLT_Dimuon16_Jpsi: boolean (nullable = true)\n",
      " |-- HLT_Dimuon10_Jpsi_Barrel: boolean (nullable = true)\n",
      " |-- HLT_Dimuon8_PsiPrime_Barrel: boolean (nullable = true)\n",
      " |-- HLT_Dimuon8_Upsilon_Barrel: boolean (nullable = true)\n",
      " |-- HLT_Dimuon0_Phi_Barrel: boolean (nullable = true)\n",
      " |-- HLT_Mu16_TkMu0_dEta18_Onia: boolean (nullable = true)\n",
      " |-- HLT_Mu16_TkMu0_dEta18_Phi: boolean (nullable = true)\n",
      " |-- HLT_TrkMu15_DoubleTrkMu5NoFiltersNoVtx: boolean (nullable = true)\n",
      " |-- HLT_TrkMu17_DoubleTrkMu8NoFiltersNoVtx: boolean (nullable = true)\n",
      " |-- HLT_Mu8: boolean (nullable = true)\n",
      " |-- HLT_Mu17: boolean (nullable = true)\n",
      " |-- HLT_Mu3_PFJet40: boolean (nullable = true)\n",
      " |-- HLT_Ele8_CaloIdM_TrackIdM_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele12_CaloIdM_TrackIdM_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele17_CaloIdM_TrackIdM_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele23_CaloIdM_TrackIdM_PFJet30: boolean (nullable = true)\n",
      " |-- HLT_Ele50_CaloIdVT_GsfTrkIdT_PFJet140: boolean (nullable = true)\n",
      " |-- HLT_Ele50_CaloIdVT_GsfTrkIdT_PFJet165: boolean (nullable = true)\n",
      " |-- HLT_PFHT400_SixJet30_DoubleBTagCSV_p056: boolean (nullable = true)\n",
      " |-- HLT_PFHT450_SixJet40_BTagCSV_p056: boolean (nullable = true)\n",
      " |-- HLT_PFHT400_SixJet30: boolean (nullable = true)\n",
      " |-- HLT_PFHT450_SixJet40: boolean (nullable = true)\n",
      " |-- HLT_Ele115_CaloIdVT_GsfTrkIdT: boolean (nullable = true)\n",
      " |-- HLT_Mu55: boolean (nullable = true)\n",
      " |-- HLT_Photon42_R9Id85_OR_CaloId24b40e_Iso50T80L_Photon25_AND_HE10_R9Id65_Eta2_Mass15: boolean (nullable = true)\n",
      " |-- HLT_Photon90_CaloIdL_PFHT600: boolean (nullable = true)\n",
      " |-- HLT_PixelTracks_Multiplicity60ForEndOfFill: boolean (nullable = true)\n",
      " |-- HLT_PixelTracks_Multiplicity85ForEndOfFill: boolean (nullable = true)\n",
      " |-- HLT_PixelTracks_Multiplicity110ForEndOfFill: boolean (nullable = true)\n",
      " |-- HLT_PixelTracks_Multiplicity135ForEndOfFill: boolean (nullable = true)\n",
      " |-- HLT_PixelTracks_Multiplicity160ForEndOfFill: boolean (nullable = true)\n",
      " |-- HLT_FullTracks_Multiplicity80: boolean (nullable = true)\n",
      " |-- HLT_FullTracks_Multiplicity100: boolean (nullable = true)\n",
      " |-- HLT_FullTracks_Multiplicity130: boolean (nullable = true)\n",
      " |-- HLT_FullTracks_Multiplicity150: boolean (nullable = true)\n",
      " |-- HLT_ECALHT800: boolean (nullable = true)\n",
      " |-- HLT_DiSC30_18_EIso_AND_HE_Mass70: boolean (nullable = true)\n",
      " |-- HLT_MET200: boolean (nullable = true)\n",
      " |-- HLT_Ele27_HighEta_Ele20_Mass55: boolean (nullable = true)\n",
      " |-- HLT_L1FatEvents: boolean (nullable = true)\n",
      " |-- HLT_Physics: boolean (nullable = true)\n",
      " |-- HLT_Physics_part0: boolean (nullable = true)\n",
      " |-- HLT_Physics_part1: boolean (nullable = true)\n",
      " |-- HLT_Physics_part2: boolean (nullable = true)\n",
      " |-- HLT_Physics_part3: boolean (nullable = true)\n",
      " |-- HLT_Random: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_part0: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_part1: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_part2: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_part3: boolean (nullable = true)\n",
      " |-- HLT_AK4CaloJet30: boolean (nullable = true)\n",
      " |-- HLT_AK4CaloJet40: boolean (nullable = true)\n",
      " |-- HLT_AK4CaloJet50: boolean (nullable = true)\n",
      " |-- HLT_AK4CaloJet80: boolean (nullable = true)\n",
      " |-- HLT_AK4CaloJet100: boolean (nullable = true)\n",
      " |-- HLT_AK4PFJet30: boolean (nullable = true)\n",
      " |-- HLT_AK4PFJet50: boolean (nullable = true)\n",
      " |-- HLT_AK4PFJet80: boolean (nullable = true)\n",
      " |-- HLT_AK4PFJet100: boolean (nullable = true)\n",
      " |-- HLT_HISinglePhoton10: boolean (nullable = true)\n",
      " |-- HLT_HISinglePhoton15: boolean (nullable = true)\n",
      " |-- HLT_HISinglePhoton20: boolean (nullable = true)\n",
      " |-- HLT_HISinglePhoton40: boolean (nullable = true)\n",
      " |-- HLT_HISinglePhoton60: boolean (nullable = true)\n",
      " |-- HLT_EcalCalibration: boolean (nullable = true)\n",
      " |-- HLT_HcalCalibration: boolean (nullable = true)\n",
      " |-- HLT_GlobalRunHPDNoise: boolean (nullable = true)\n",
      " |-- HLT_L1BptxMinus: boolean (nullable = true)\n",
      " |-- HLT_L1BptxPlus: boolean (nullable = true)\n",
      " |-- HLT_L1NotBptxOR: boolean (nullable = true)\n",
      " |-- HLT_L1BeamGasMinus: boolean (nullable = true)\n",
      " |-- HLT_L1BeamGasPlus: boolean (nullable = true)\n",
      " |-- HLT_L1BptxXOR: boolean (nullable = true)\n",
      " |-- HLT_L1MinimumBiasHF_OR: boolean (nullable = true)\n",
      " |-- HLT_L1MinimumBiasHF_AND: boolean (nullable = true)\n",
      " |-- HLT_HcalNZS: boolean (nullable = true)\n",
      " |-- HLT_HcalPhiSym: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_FirstCollisionAfterAbortGap: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_FirstCollisionAfterAbortGap_TCDS: boolean (nullable = true)\n",
      " |-- HLT_ZeroBias_IsolatedBunches: boolean (nullable = true)\n",
      " |-- HLT_Photon500: boolean (nullable = true)\n",
      " |-- HLT_Photon600: boolean (nullable = true)\n",
      " |-- HLT_Mu300: boolean (nullable = true)\n",
      " |-- HLT_Mu350: boolean (nullable = true)\n",
      " |-- HLT_MET250: boolean (nullable = true)\n",
      " |-- HLT_MET300: boolean (nullable = true)\n",
      " |-- HLT_MET600: boolean (nullable = true)\n",
      " |-- HLT_MET700: boolean (nullable = true)\n",
      " |-- HLT_PFMET300: boolean (nullable = true)\n",
      " |-- HLT_PFMET400: boolean (nullable = true)\n",
      " |-- HLT_PFMET500: boolean (nullable = true)\n",
      " |-- HLT_PFMET600: boolean (nullable = true)\n",
      " |-- HLT_HT2000: boolean (nullable = true)\n",
      " |-- HLT_HT2500: boolean (nullable = true)\n",
      " |-- HLT_IsoTrackHE: boolean (nullable = true)\n",
      " |-- HLT_IsoTrackHB: boolean (nullable = true)\n",
      " |-- HLT_Ele27_WPTight_Gsf_L1JetTauSeeded: boolean (nullable = true)\n",
      " |-- HLT_Ele45_WPLoose_Gsf_L1JetTauSeeded: boolean (nullable = true)\n",
      " |-- HLT_Ele23_Ele12_CaloIdL_TrackIdL_IsoVL_DZ_L1JetTauSeeded: boolean (nullable = true)\n",
      " |-- HLTriggerFinalPath: boolean (nullable = true)\n",
      " |-- Flag_duplicateMuons: boolean (nullable = true)\n",
      " |-- Flag_badMuons: boolean (nullable = true)\n",
      " |-- Flag_noBadMuons: boolean (nullable = true)\n",
      " |-- Flag_HBHENoiseFilter: boolean (nullable = true)\n",
      " |-- Flag_HBHENoiseIsoFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHaloFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHaloTrkMuUnvetoFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHalo2015Filter: boolean (nullable = true)\n",
      " |-- Flag_globalTightHalo2016Filter: boolean (nullable = true)\n",
      " |-- Flag_globalSuperTightHalo2016Filter: boolean (nullable = true)\n",
      " |-- Flag_HcalStripHaloFilter: boolean (nullable = true)\n",
      " |-- Flag_hcalLaserEventFilter: boolean (nullable = true)\n",
      " |-- Flag_EcalDeadCellTriggerPrimitiveFilter: boolean (nullable = true)\n",
      " |-- Flag_EcalDeadCellBoundaryEnergyFilter: boolean (nullable = true)\n",
      " |-- Flag_goodVertices: boolean (nullable = true)\n",
      " |-- Flag_eeBadScFilter: boolean (nullable = true)\n",
      " |-- Flag_ecalLaserCorrFilter: boolean (nullable = true)\n",
      " |-- Flag_trkPOGFilters: boolean (nullable = true)\n",
      " |-- Flag_chargedHadronTrackResolutionFilter: boolean (nullable = true)\n",
      " |-- Flag_muonBadTrackFilter: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_manystripclus53X: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_toomanystripclus53X: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_logErrorTooManyClusters: boolean (nullable = true)\n",
      " |-- Flag_METFilters: boolean (nullable = true)\n",
      " |-- L1Reco_step: boolean (nullable = true)\n",
      " |-- Flag_HBHENoiseFilter: boolean (nullable = true)\n",
      " |-- Flag_HBHENoiseIsoFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHaloFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHaloTrkMuUnvetoFilter: boolean (nullable = true)\n",
      " |-- Flag_CSCTightHalo2015Filter: boolean (nullable = true)\n",
      " |-- Flag_globalTightHalo2016Filter: boolean (nullable = true)\n",
      " |-- Flag_globalSuperTightHalo2016Filter: boolean (nullable = true)\n",
      " |-- Flag_HcalStripHaloFilter: boolean (nullable = true)\n",
      " |-- Flag_hcalLaserEventFilter: boolean (nullable = true)\n",
      " |-- Flag_EcalDeadCellTriggerPrimitiveFilter: boolean (nullable = true)\n",
      " |-- Flag_EcalDeadCellBoundaryEnergyFilter: boolean (nullable = true)\n",
      " |-- Flag_goodVertices: boolean (nullable = true)\n",
      " |-- Flag_eeBadScFilter: boolean (nullable = true)\n",
      " |-- Flag_ecalLaserCorrFilter: boolean (nullable = true)\n",
      " |-- Flag_trkPOGFilters: boolean (nullable = true)\n",
      " |-- Flag_chargedHadronTrackResolutionFilter: boolean (nullable = true)\n",
      " |-- Flag_muonBadTrackFilter: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_manystripclus53X: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_toomanystripclus53X: boolean (nullable = true)\n",
      " |-- Flag_trkPOG_logErrorTooManyClusters: boolean (nullable = true)\n",
      " |-- Flag_METFilters: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_DiJet170_Mu5: boolean (nullable = true)\n",
      " |-- HLT_BTagMu_AK8Jet300_Mu5: boolean (nullable = true)\n",
      " |-- HLT_L1FatEvents_part0: boolean (nullable = true)\n",
      " |-- HLT_L1FatEvents_part1: boolean (nullable = true)\n",
      " |-- HLT_L1FatEvents_part2: boolean (nullable = true)\n",
      " |-- HLT_L1FatEvents_part3: boolean (nullable = true)\n",
      " |-- HLT_Ele250_CaloIdVT_GsfTrkIdT: boolean (nullable = true)\n",
      " |-- HLT_Ele300_CaloIdVT_GsfTrkIdT: boolean (nullable = true)\n",
      " |-- sample: string (nullable = false)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "DFList[0].printSchema()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Data reduction\n",
    "\n",
    "Subsets of interesting attributes can be selected via 'select' operations on the DataFrames (equivalent to \"pruning\" steps in ROOT-based frameworks).\n",
    "\n",
    "All datasets can be joint into a single DataFrame (e.g. collecting data from various samples)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "root\n",
      " |-- Muon_pt: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_eta: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_phi: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_mass: array (nullable = true)\n",
      " |    |-- element: float (containsNull = true)\n",
      " |-- Muon_charge: array (nullable = true)\n",
      " |    |-- element: integer (containsNull = true)\n",
      " |-- Muon_mediumId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Muon_softId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- Muon_tightId: array (nullable = true)\n",
      " |    |-- element: boolean (containsNull = true)\n",
      " |-- nMuon: integer (nullable = true)\n",
      " |-- sample: string (nullable = false)\n",
      "\n",
      "Partitions: 3\n"
     ]
    }
   ],
   "source": [
    "columns = [\n",
    "    ### MUON\n",
    "    'Muon_pt',\n",
    "    'Muon_eta',\n",
    "    'Muon_phi',\n",
    "    'Muon_mass',\n",
    "    'Muon_charge',\n",
    "    'Muon_mediumId',\n",
    "    'Muon_softId',\n",
    "    'Muon_tightId',\n",
    "    'nMuon',\n",
    "    ### SAMPLE\n",
    "    'sample',\n",
    "]\n",
    "\n",
    "# Select columns from dataframe\n",
    "DF = DFList[0].select(columns)\n",
    "DF.printSchema()\n",
    "\n",
    "# Merge all dataset into a single dataframe\n",
    "for df_ in DFList[1:]:\n",
    "    DF = DF.union(df_.select(columns))\n",
    "\n",
    "print 'Partitions: {}'.format(DF.rdd.getNumPartitions())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Data selection\n",
    "\n",
    "Selection of events based on features is obtained via a 'filter' operation.\n",
    "\n",
    "Number of entries is obtained by 'count'."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "total number of events in the DataFrame  =  1451024\n",
      "events in the DataFrame with \"nMuon > 0\" =  760701\n"
     ]
    }
   ],
   "source": [
    "print 'total number of events in the DataFrame  = ', DF.count()\n",
    "print 'events in the DataFrame with \\\"nMuon > 0\\\" = ', DF.filter('nMuon > 0').count()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Caching - 1\n",
    "\n",
    "Dataframes can be cached into memory, shared across the Spark cluster nodes, for a faster access.\n",
    "\n",
    "Caching takes some time..."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "DF = DF.cache()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Caching - 2\n",
    "\n",
    "... but ensures fast data-handling operations afterwards."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "total number of events in the DataFrame  =  1451024\n",
      "events in the DataFrame with \"nMuon > 0\" =  760701\n"
     ]
    }
   ],
   "source": [
    "print 'total number of events in the DataFrame  = ', DF.count()\n",
    "print 'events in the DataFrame with \\\"nMuon > 0\\\" = ', DF.filter('nMuon > 0').count()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Data inspection\n",
    "\n",
    "\n",
    "Events can be inspected with 'show' (as in TTree.Show() ), also concatenating 'select' and 'filter'."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "+----------+-----+-----------+------------+------------+-----------+\n",
      "|    sample|nMuon|    Muon_pt|    Muon_eta|    Muon_phi|Muon_charge|\n",
      "+----------+-----+-----------+------------+------------+-----------+\n",
      "|DYJetsToLL|    1| [34.75507]|[-1.3212891]|[-1.0375977]|        [1]|\n",
      "|DYJetsToLL|    0|         []|          []|          []|         []|\n",
      "|DYJetsToLL|    1|[3.3921983]|[-1.5285645]|[-0.2514038]|       [-1]|\n",
      "|DYJetsToLL|    0|         []|          []|          []|         []|\n",
      "|DYJetsToLL|    0|         []|          []|          []|         []|\n",
      "+----------+-----+-----------+------------+------------+-----------+\n",
      "only showing top 5 rows\n",
      "\n"
     ]
    }
   ],
   "source": [
    "DF.filter(DF['sample'] == 'DYJetsToLL')\\\n",
    "  .select('sample','nMuon','Muon_pt','Muon_eta','Muon_phi','Muon_charge')\\\n",
    "  .show(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create derivate quantities and structures - 1\n",
    "\n",
    "User defined functions can be used for transformations evalueted row by row to compute derived quantity, such as invaraint mass of two physics objects involving multiple column.\n",
    "The return value is added as a new column in the output DataFrame.\n",
    "\n",
    "Dimuon candidate structure is created as an example."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from pyspark.sql.types import *\n",
    "\n",
    "dimuonSchema = StructType([\n",
    "    StructField(\"pass\", BooleanType(), False),   # True if filled / False if default(empty) \n",
    "    #\n",
    "    StructField(\"mass\", FloatType(), False),     # Dimuon mass\n",
    "    StructField(\"pt\", FloatType(), False),       # Dimuon pt\n",
    "    StructField(\"eta\", FloatType(), False),      # Dimuon eta\n",
    "    StructField(\"phi\", FloatType(), False),      # Dimuon phi\n",
    "    StructField(\"dPhi\", FloatType(), False),     # DeltaPhi(mu1,mu2)\n",
    "    StructField(\"dR\", FloatType(), False),       # DeltaR(mu1,mu2)\n",
    "    StructField(\"dEta\", FloatType(), False),     # DeltaEta(mu1,mu2)\n",
    "    #\n",
    "    StructField(\"mu1_pt\", FloatType(), False),   # leading mu pT \n",
    "    StructField(\"mu2_pt\", FloatType(), False),   # sub-leading mu pT \n",
    "    StructField(\"mu1_eta\", FloatType(), False),  # leading mu eta\n",
    "    StructField(\"mu2_eta\", FloatType(), False),  # sub-leading mu eta\n",
    "    StructField(\"mu1_phi\", FloatType(), False),  # leading mu phi\n",
    "    StructField(\"mu2_phi\", FloatType(), False),  # sub-leading mu phi\n",
    "])\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create derivate quantities and structures - 2\n",
    "\n",
    "And a generic function filling the candidate structure can be defined."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from math import *\n",
    "\n",
    "def deltaPhi(phi1,phi2):\n",
    "    ## Catch if being called with two objects\n",
    "    if type(phi1) != float and type(phi1) != int:\n",
    "        phi1 = phi1.phi\n",
    "    if type(phi2) != float and type(phi2) != int:\n",
    "        phi2 = phi2.phi\n",
    "    ## Otherwise\n",
    "    dphi = (phi1-phi2)\n",
    "    while dphi >  pi: dphi -= 2*pi\n",
    "    while dphi < -pi: dphi += 2*pi\n",
    "    return dphi\n",
    "\n",
    "def deltaR(eta1,phi1,eta2=None,phi2=None):\n",
    "    ## catch if called with objects\n",
    "    if eta2 == None:\n",
    "        return deltaR(eta1.eta,eta1.phi,phi1.eta,phi1.phi)\n",
    "    ## otherwise\n",
    "    return hypot(eta1-eta2, deltaPhi(phi1,phi2))\n",
    "\n",
    "def invMass(pt1, pt2, eta1, eta2, phi1, phi2, mass1, mass2):\n",
    "    #\n",
    "    theta1 = 2.0*atan(exp(-eta1))\n",
    "    px1    = pt1 * cos(phi1)\n",
    "    py1    = pt1 * sin(phi1)\n",
    "    pz1    = pt1 / tan(theta1)\n",
    "    E1     = sqrt(px1**2 + py1**2 + pz1**2 + mass1**2)\n",
    "\n",
    "    theta2 = 2.0*atan(exp(-eta2))\n",
    "    px2    = pt2 * cos(phi2)\n",
    "    py2    = pt2 * sin(phi2)\n",
    "    pz2    = pt2 / tan(theta2)\n",
    "    E2     = sqrt(px2**2 + py2**2 + pz2**2 + mass2**2)\n",
    "\n",
    "    themass  = sqrt((E1 + E2)**2 - (px1 + px2)**2 - (py1 + py2)**2 - (pz1 + pz2)**2)\n",
    "    thept    = sqrt((px1 + px2)**2 + (py1 + py2)**2)\n",
    "    thetheta = atan( thept / (pz1 + pz2) )        \n",
    "    theeta   = 0.5*log( (sqrt((px1 + px2)**2 + (py1 + py2)**2 + (pz1 + pz2)**2)+(pz1 + pz2))/(sqrt((px1 + px2)**2 + (py1 + py2)**2 + (pz1 + pz2)**2)-(pz1 + pz2)) )\n",
    "    thephi   = asin((py1 + py2)/thept)\n",
    "\n",
    "    delPhi = deltaPhi(phi1,phi2)\n",
    "    delR   = deltaR(eta1,phi1,eta2,phi2)\n",
    "    delEta = eta1-eta2\n",
    "\n",
    "    return (themass, thept, theeta, thephi, delPhi, delR, delEta)\n",
    "\n",
    "def dimuonCandidate(pt, eta, phi, mass, charge, mediumid):\n",
    "    # default class implementation   \n",
    "    default_ = (False, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)\n",
    "    \n",
    "    \"\"\"\n",
    "    Z->mm candidate from arbitrary muon selection:\n",
    "      N(mu) >= 2\n",
    "      pT > 30, 10\n",
    "      abs(eta) < 2.4, 2.4\n",
    "      mediumId muon\n",
    "      opposite charge\n",
    "    \"\"\"\n",
    "    \n",
    "    if len(pt) < 2:\n",
    "        return default_\n",
    "    \n",
    "    #Identify muon candidate\n",
    "    leadingIdx = None\n",
    "    trailingIdx = None\n",
    " \n",
    "    for idx in range(len(pt)):\n",
    "        if leadingIdx == None:\n",
    "            if pt[idx] > 30 and abs(eta[idx]) < 2.4 and mediumid[idx]:\n",
    "                leadingIdx = idx\n",
    "        elif trailingIdx == None:\n",
    "            if pt[idx] > 10 and abs(eta[idx]) < 2.4 and mediumid[idx]:\n",
    "                trailingIdx = idx\n",
    "        else:\n",
    "            if pt[idx] > 10 and abs(eta[idx]) < 2.4 and mediumid[idx]:\n",
    "                return default_\n",
    "\n",
    "    if leadingIdx != None and trailingIdx != None and charge[leadingIdx] != charge[trailingIdx]:            \n",
    "        # Candidate found\n",
    "        dimuon_   = (True,) + \\\n",
    "                    invMass(pt[leadingIdx], pt[trailingIdx],\n",
    "                            eta[leadingIdx], eta[trailingIdx],\n",
    "                            phi[leadingIdx], phi[trailingIdx],\n",
    "                            mass[leadingIdx], mass[trailingIdx]) + \\\n",
    "                    (pt[leadingIdx], pt[trailingIdx],\n",
    "                     eta[leadingIdx], eta[trailingIdx],\n",
    "                     phi[leadingIdx], phi[trailingIdx])\n",
    "                \n",
    "        return dimuon_\n",
    "    else:\n",
    "        return default_    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Create derivate quantities and structures - 3\n",
    "\n",
    "Finally, a dimuon candidate structure can be appended to the DataFrame as an additional column"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "+--------------------+\n",
      "|              Dimuon|\n",
      "+--------------------+\n",
      "|[true,89.35462,18...|\n",
      "|[true,88.32019,21...|\n",
      "|[true,92.423805,5...|\n",
      "+--------------------+\n",
      "only showing top 3 rows\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from pyspark.sql.functions import udf\n",
    "dimuonUDF = udf(dimuonCandidate, dimuonSchema)\n",
    "\n",
    "DF = DF.withColumn('Dimuon', dimuonUDF (\"Muon_pt\",\n",
    "                                        \"Muon_eta\",\n",
    "                                        \"Muon_phi\",\n",
    "                                        \"Muon_mass\",\n",
    "                                        \"Muon_charge\",\n",
    "                                        \"Muon_mediumId\")\n",
    "                  )\n",
    "\n",
    "DF.where('Dimuon.pass == True').select('Dimuon').show(3)\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get statistics info about the data\n",
    "\n",
    "Exploit pySparkSql functions to get statistical insights on the DataFrame"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Number of events, pre-selection level\n",
      "+----------+------+\n",
      "|    sample| count|\n",
      "+----------+------+\n",
      "|DYJetsToLL|549790|\n",
      "|SingleMuon|477186|\n",
      "|        WZ|424048|\n",
      "+----------+------+\n",
      "\n",
      "Number of events, Dimuon invariant mass in [70-110] GeV\n",
      "+----------+-----+\n",
      "|    sample|count|\n",
      "+----------+-----+\n",
      "|DYJetsToLL|67831|\n",
      "|SingleMuon|13934|\n",
      "|        WZ| 7209|\n",
      "+----------+-----+\n",
      "\n",
      "Mean of Dimuon mass, evaluated in [70-110] GeV range\n",
      "+----------+--------------------------+\n",
      "|    sample|avg(Dimuon.mass AS `mass`)|\n",
      "+----------+--------------------------+\n",
      "|DYJetsToLL|         90.19341036229169|\n",
      "|SingleMuon|         89.91447557800039|\n",
      "|        WZ|         90.19251177561962|\n",
      "+----------+--------------------------+\n",
      "\n",
      "Description of Dimuon mass features for SingleMuon dataset only, evaluated in [70-110] GeV range\n",
      "+-------+-----------------+\n",
      "|summary|      Dimuon.mass|\n",
      "+-------+-----------------+\n",
      "|  count|            13934|\n",
      "|   mean|89.91447557800039|\n",
      "| stddev|5.392174086677417|\n",
      "|    min|         70.00491|\n",
      "|    max|         109.9901|\n",
      "+-------+-----------------+\n",
      "\n"
     ]
    }
   ],
   "source": [
    "from pyspark.sql.functions import *\n",
    "\n",
    "print 'Number of events, pre-selection level'\n",
    "\n",
    "DF.groupBy(\"sample\").count().show()\n",
    "\n",
    "print 'Number of events, Dimuon invariant mass in [70-110] GeV'\n",
    "\n",
    "DF.where( (col(\"Dimuon.mass\") > 70) & (col(\"Dimuon.mass\") < 110) ).groupBy(\"sample\").count().show()\n",
    "\n",
    "print 'Mean of Dimuon mass, evaluated in [70-110] GeV range'\n",
    "\n",
    "DF.where( (col(\"Dimuon.mass\") > 70) & (col(\"Dimuon.mass\") < 110) ).groupBy('sample').mean('Dimuon.mass').show()\n",
    "\n",
    "print 'Description of Dimuon mass features for SingleMuon dataset only, evaluated in [70-110] GeV range'\n",
    "\n",
    "DF.where( (col(\"Dimuon.mass\") > 70) & (col(\"Dimuon.mass\") < 110) & (DF[\"sample\"] == \"SingleMuon\") ).describe('Dimuon.mass').show()\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plotting the Zpeak mass"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# add histogrammar to driver python path\n",
    "sc.addPyFile(\"/eos/project/s/swan/public/PhysicsAnalysis/histogrammar.zip\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Load libraries, and append histogrammar functionalities to dataframe\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "import histogrammar as hg\n",
    "import histogrammar.sparksql\n",
    "import numpy as np\n",
    "\n",
    "DF = DF.where( (col(\"Dimuon.mass\") > 70) & (col(\"Dimuon.mass\") < 110) )\n",
    "\n",
    "hg.sparksql.addMethods(DF)\n",
    "\n",
    "\n",
    "plots = hg.UntypedLabel(\n",
    "    # 1d histograms\n",
    "    LeadPt       = hg.Bin(50, 30, 180,   DF['Dimuon.mu1_pt']),\n",
    "    LeadPtEta    = hg.Bin(48, -2.4, 2.4, DF['Dimuon.mu1_eta']),\n",
    "    SubLeadPt    = hg.Bin(100, 0, 200,   DF['Dimuon.mu2_pt']),\n",
    "    SubLeadPtEta = hg.Bin(48, -2.4, 2.4, DF['Dimuon.mu2_eta']),\n",
    "    InvMass      = hg.Bin(80, 70, 110,   DF['Dimuon.mass']),\n",
    "    DeltaR       = hg.Bin(50, 0, 5,      DF['Dimuon.dPhi']),\n",
    "    DeltaPhi     = hg.Bin(64, -3.2, 3.2, DF['Dimuon.dR']),\n",
    ")\n",
    "\n",
    "# Make a set of histograms, categorized per-sample\n",
    "bulkHisto = hg.Categorize(quantity = DF['sample'], value = plots)\n",
    "\n",
    "# Fill from spark\n",
    "bulkHisto.fillsparksql(df=DF)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x7fd0a7f308d0>"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGiCAYAAAAfnjf+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAMTQAADE0B0s6tTgAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xt40/Xd//Hnt01CLafSCgjUhrYp\nWouKIBuaWqbFtXNszCmbJxQPk02Z98aube52Cv7cPM2bzdMudVOZY1PAwxzqqKeN2g7vG1ScVkZp\ngBSwDqWgYqFp2u/vj9DYQ5ImbXPs63FduUq+3yR9f2mbvvo5GqZpmoiIiIgkubR4FyAiIiIyGBRq\nREREJCUo1IiIiEhKUKgRERGRlKBQIyIiIilBoUZERERSgkKNiIiIpASFGhEREUkJCjUiIiKSEhRq\nREREJCVY4l1ArAwbNoyxY8fGuwwRERGJwIcffkhra2tYjx0yoWbs2LHs3r073mWIiIhIBHJzc8N+\nrLqfREREJCUo1IiIiEhKUKgRERGRlKBQIyIiIilBoUZERERSwpCZ/SQiItFhmqb/JhIOwzD8t8Gk\nUCMiIv3W3NzMvn378Hq98S5FkozFYiEnJ4fs7OzBe81BeyURERlSWlpa+PDDD5k4cSKZmZnxLkeS\nTEtLC++//z4ZGRmD9v2jUCMiIv2yd+9ecnJyGDlyZLxLkSQ0cuRIcnJy2Lt3L5MnTx6U19RAYRER\niZhpmrS2tirQyICMHDmS1tbWQRuPpVAjIiIRM02Tjo4OLBY1+Ev/WSwWOjo6FGpERCR+NNNJBpNC\njYiIiEgXCjUiIjLoDnsP80nrJ3G7HfYe7lfdK1as6LaGyvDhw5k8eTLnnnsuq1evpqOjA4DnnnsO\nwzB46KGHer1Ga2srxcXFzJw5k/b2dv9rNjQ0RFzPsmXLePXVVyN6zuTJk7tdQ7BbJC644AIcDkfI\nx8yaNYs5c+ZE9LqDTZ2hIiIyqA57DzNp+SSaDzXHrYbso7LZs2QPGZaMfj1/zZo15Obm0traSmNj\nI88//zwXXnghDz30EGvXrmXu3LlcfPHF/OQnP2Hu3LlMnDjR/9xbbrkFl8vFG2+8QXp6+oCu4+ab\nb+aGG27grLPOCvs5zzzzDK2trf7711xzDe3t7Tz44IMDqiUZKNSIiMig8rR74hpoAJoPNeNp9/Q7\n1EybNq1by8SCBQuYP38+8+fP5yc/+Qn33nsvd999Ny+//DLXXHMNf/nLXwD417/+xZ133skNN9zA\niSeeOCjXEqlTTjml2/1Ro0bh9XqZNWtWXOqJJXU/iciQUFJSQklJSbzLkCR23nnnMW/ePH73u9/R\n0tJCTk4O9913H88++yyrV6+mvb2dK6+8kuLiYv77v/+7z9d7+umnmTVrFpmZmWRlZTF//nwaGxv9\n5zu7iH75y1/6u4yWLVsGwMaNGzn77LPJyckhMzOTgoICrrnmmn5d19q1a/nCF77AUUcdxZgxYzjv\nvPNwuVz9eq14U6gREREJ0znnnENrayubNm0C4Pzzz+e8887j+9//PjfeeCNvvfUWjzzyCFarNeTr\nPPDAA5x33nmccMIJPPnkkzz44IO8++67zJ49m08//RSADRs2ALBw4UI2bNjAhg0buOqqqzh48CAV\nFRWkp6ezYsUKXnjhBW666aZ+bVXx7LPPMm/ePI4++mhWr17NvffeyxtvvEFpaSl79+6N+PXiTd1P\nIiIiYcrLywOgqanJf+z++++npKSE2267jeuvv54ZM2aEfI2DBw/y05/+lMsvv5xHHnnEf/yLX/wi\nU6ZM4eGHH+YHP/iBv7to0qRJ3bqONm3axP79+7nzzjs56aST/McXLlwY8fXccMMNFBcX89xzz5GW\n5mvnmDlzJieccAK/+c1vuPXWWyN+zXhSS42IiEiYOtdT6Tp7aPz48Vx99dUA3HjjjX2+xoYNG/jk\nk0+4+OKL8Xq9/ltubi7HH3881dXVIZ9fVFREVlYWixYtYuXKlezatatf19Lc3ExdXR0XXnihP9AA\nHHfcccycOZP169f363XjSaFGREQkTJ0BYsKECd2O22y2bh9D6ezWmTNnDlartdvtnXfeYd++fSGf\nP3r0aP7+978zceJErrnmGvLy8pg6dSpPPfVURNfS3Nwc8FoAjjnmGP/5ZKLuJxERkTA9//zzZGRk\n9NnFFEpOTg7gWxMn0OD1cPbTmjZtGk899RRer5dNmzZx22238a1vfYu3336bqVOnhlVHdnY2AB98\n8EGvcx988IG/zmSilhoRSXmmadLS0sKBAweoqanREv/SL08//TR//etf+e53v0tmZma/X+f0009n\n5MiRNDQ0cOqpp/a6HXfccf7H2mw2Dh06FPS1LBYLs2bN4pZbbqGjo4MtW7aEXUd2djZTp05l9erV\n3X4mtm3bxqZNm5g9e3b/LjCO1FIjIinN7XZTUVGB2+0GoLy8nPz8fKqqqrDb7XGuThLV5s2b+eij\nj/B4PDQ2NvLcc8+xZs0azj77bG677bZ+vWbnOJxRo0bxq1/9imuvvZYPP/yQr3zlK4wePZo9e/aw\nfv16vvSlL3HRRRcBcMIJJ/D8889TWVnJmDFjmDhxIm+++SYPPfQQ3/jGN8jPz+ezzz7jnnvuYeTI\nkZx22mkR1fSLX/yCc889l3nz5rFo0SIOHDjAjTfeyNixY/mv//qvbo/97LPPePLJJ3u9xnHHHedf\nk2fv3r0BH3PKKadQWFgYUW39oVAjIinLNE0qKipwuVz+v0Q9Hg8ul4vKykree++9iJeLl6Fh/vz5\nAGRkZDBu3DimT5/OE088wfnnnx/x90xnS8uwYcP8xxYtWsSxxx7Lr371K/785z/T1tbGpEmTKCsr\nY9q0af7H3XfffVx33XV87Wtfo7W1laVLl3LhhRdy1FFHccstt9DU1MTIkSOZOXMmL730Erm5uRHV\nNm/ePJ599lluueUWzjvvPDIyMigvL+fOO+9k3Lhx3R77wQcf+P9fuvrRj37EXXfdBcA777wT8DH3\n3nsvixcvjqi2/jDMIdIOm5uby+7du+NdhojEUE1NDWVlZQG7m2w2G6+88gqlpaVxqCz5tbe3U19f\nz5QpU3ptBZAK2yQMpiVLlvDb3/6WgwcPYrGoLaGrUN9HnSL5/a3/XRFJWQ0NDUHHz9hsNhoaGhRq\noiDDksGeJXvwtHviVoMt3Rb3QFNXV8frr7/OI488wrnnnqtAEwP6HxaRlOVwOCAdaO99zuPx9Lnr\nsPRfhiUj7qEi3hYvXkx9fT3nn38+y5cvj3c5Q4JCjYikLKfTCVlAM9C1wcaAgoIC33mRKPn73/8e\n7xKGHE3pFpGU0nXjSsMwYAGQ3eUB6UAOVFVVYRiGNroUSSFqqRGR1JYFLAbuxtcNdT6Q9/kePiKS\nOhRqRCT1GYD1yE1L04ikrKh3P1133XVMnjwZwzDYvHmz//i2bds4/fTTmTJlCjNnzqSurm7A50RE\nRGToinqoOf/886mpqem1cueiRYu4+uqrqa+v56c//Wm3LdP7e05ERESGrqiHmrKysl4rHO7du5dN\nmzZxySWXAHDeeeexa9cuGhoa+n1OREREhra4jKnZtWsXEyZM8C9EZBgGeXl5NDY2Mnr06H6d03oT\nImLcbMCHXf7dh86NLj0eDzU1NTidTm2bIJLEEmZKd6jdGvpzbvny5eTm5vpvBw8eHHCNIpLgTKAN\nOAy46b42zbVHbke43W6Ki4txu900NTVRXl7uvy8Dd/hwB5980h632+HDHRHX/Pjjj2MYBtXV1d2O\n/+c//8EwDMaPH9/rOffffz+GYfDuu+9iGEafN4muuLTUHHvssTQ1NeH1erFYLJimya5du8jLy2PU\nqFH9OtfTkiVLWLJkif9+pJt8iUhycbvdcB9w4MiBx/BN515w5GNXJtroMooOH+5g0qR3aG4OsJRz\njGRnp7Nnz4lkZIT/t/vs2bMBqK6upqyszH+8urqazMxM9u7dy7///W+OP/74budycnIoKipiw4YN\nAV/3jTfeYPHixXzta1/r59VIuOLSUtO54+nKlSsBeOqpp8jNzcXhcPT7nIgMXZ27cdN1/8R2fPdX\n0r3FBqAR6uvr8Xq93Q57vV62b99ObW1tdAtOcR6PGddAA9Dc3I7HE9l+zRMnTqSgoKBXS011dTVn\nnXUWhYWFvc699tprnHHGGQwbNoxZs2b1uh1//PH8z//8D3a7nT/84Q8Dvi4JLeqhZtGiRf4dNisq\nKvwB5MEHH+TBBx9kypQp3H777Tz66KP+5/T3nIgMTbW1tezcubN3eDGB/UBjj+MhNo/u3OhShqbZ\ns2ezYcOGboG3urqaM844g9LS0m6hZtu2bTQ1NXVr1enKNE0uvfRS9uzZw5o1axgzZkzU6x/qot79\n9OCDDwY8ftxxxwVtquvvOREZmhoaGrBarbS2tvY+mYYvxHRdVSIbrFYrHk/vXaS10eXQVlZWxqOP\nPsqbb77JF77wBQ4cOMC7777LGWecQU5ODjfffLP/sZ0BJ1iouf3221m7di333nsvM2fOjEn9Q13C\nDBQWEekvh8MRMKAA0EH3vZ8A8iA/P98/k7KTxWLRRpdDXGdA6Qwsr732GsOGDWPGjBmcccYZ7Nq1\ny9cqeOQxo0aNYtq0ab1e5x//+Ac33ngj8+fPZ/HixTGrf6hTqBGRpOd0OsnPz/dth9CVAYwBes4l\nMHwbWhYWFvpnpdhsNhwOh3+jSxmaCgoKyM3N9Yea6upqvvjFL2Kz2ZgyZQrjxo3rds7pdJKent7t\nNZqamrjgggsoKCjg4Ycfjvk1DGUKNSKS9AzDoKqqKuBu3FxC77AD2O12tmzZgt1uZ8KECbzyyiu8\n99572uhSKCsro6amBtM0/eNpOnWOq9m9ezc7d+7s1fXk9Xq54IIL+OSTT3jyyScZOXJkrMsf0hRq\nRCQl2O12327cWcBI4FJ869L0nM7dhWEYZGZmkpWVRWlpqVpoBPCFmv379/P666/z5ptvdgs1Z5xx\nBtXV1axfvx74fBp4p5/97GdUV1dz7733ctJJJ8W0blGoEZFU0rkbdwa+gcH9zCglJSWUlJQMXl2S\nVDqDyu23345pmpx22mn+c6WlpWzbto3Vq1eTmZnJqaee6j/37LPPctddd3HZZZdx5ZVXxrxuidPi\neyIiIonq+OOPZ9y4caxdu5YZM2YwYsQI/7lTTjmFESNGsHbtWs4880ysVisAO3bs4LLLLiMnJ4dL\nL72U119/PeBrn3LKKQwbNiwm1zEUKdSIiIj0UFZWxpNPPtmt6wkgPT2d0047jZdeeqnbeJr169fz\n8ccfA1BeXh70dXfs2MHkyZOjUrMo1IiIyCCz2Qyys9Pjvk2Czdb/MVJr1qwJeu7FF1/sdWzhwoUs\nXLiw359PBodCjYgMaXV1dfEuIeVkZKSxZ8+JEW9TMJhsNiOifZ8kNSjUiEhqubbvh0j0ZWSkkZER\n7ypkqFGMFRHpwjRNWlpaOHDggH+tEhFJDmqpEZGkY9wcnfVk3G43FRUVuN1uwDfgMz8/n6qqKt86\nOCKS0NRSIyKCr4WmoqICl8uFaZqYponH48HlclFZWakWG5EkoJYaERmSerX2uMGoN3qFF6/Xy/bt\n26mtraW0tDSGFYpIpNRSIyIC0Bz8lM1mo6GhIXa1iEi/KNSIiABk418dtiePx4PD4YhxQSISKYUa\nERGAPMjPz8di6d4rb7FYKCgowOl0xqmwoUd7b0l/KdSIiAAYUFVVRWFhIYZhYBgGNpsNh8NBVVWV\ndvAWSQIKNSIiR9jtdrZs2YLdbmfChAm88sorvPfee+Tl5cW7NImRFStW+EOtYRgMHz6cyZMnc+65\n57J69Wo6Ojoifs3NmzezbNkymptDDNySQaFQIyLShWEYZGZmkpWVRWlpqVpohqg1a9awYcMGXnjh\nBW655RaGDRvGhRdeyJe//GUOHToU0Wtt3ryZm2++WaEmBhRqREQkYSTKis7Tpk1j1qxZzJ49mwUL\nFvDEE0+wevVqXn31VX7yk5/EpSbpm0KNiIgkBLfbTXFxMW63m6amJsrLy/33E8F5553HvHnz+N3v\nfkdLSwsAS5cuZfr06YwePZqjjz6as846i9dff93/nBUrVnD55ZcDUFRU5O/W2rlzJwD33Xcfp512\nGtnZ2WRlZTFr1iyef/75mF9bqlCoERGRuEuWFZ3POeccWltb2bRpEwB79uzhhz/8IX/5y19YsWIF\n48aNo6ysjH/9618AfPWrX+XnP/858HmX1oYNG5gwYQIAO3fu5KqrrmLNmjWsWrWKU089lblz5/K3\nv/0tPheY5LSisIiIxFzPKdstLS243e6AKzpv3bqVgoICMjMz/cfr6upiUmdPnYPGm5qaAPj973/v\nP9fe3k5lZSUlJSU8/PDD3H333YwdO5bCwkLA16XVc72ju+66y//vjo4OysvLqa+v54EHHuArX/lK\ntC8n5ailRkSkh7q6urj90hyqPB7PgM7HSmfo6hxA/vLLL3PmmWeSk5ODxWLBarVSX1/P1q1bw3q9\nN954g7lz5zJ+/Hj/81966aWwny/dqaVGRERirmdorKmpoby8PGB4sVqtrFq1KiH23tq1axcAEyZM\n4M033+Scc86hoqKChx9+mAkTJpCens5VV13F4cOHw3qt8vJyTjjhBO69917y8vKwWCzceOONbNmy\nJdqXkpIUakREJO6cTif5+fm4XC68Xq//eKKt6Pz888+TkZHBjBkz+OUvf4nFYuHpp5/utsXG/v37\nycrK6vO11q1bx8cff8zq1avJzc31H+8chCyRU/eTiIjEnWEYCb+i89NPP81f//pXvvvd75KZmUlL\nSwvp6endanv11VdpbGzs9rxhw4YB9FrfpjO8dA1E9fX11NbWRusSUp5aakREJCF0ruhcUFCAx+Nh\n1apVOJ3OuASazZs389FHH+HxeGhsbOS5555jzZo1nH322dx2220AVFZW8pvf/IaFCxdy+eWXU19f\nzy233MKkSZO6vdYJJ5wAwP33389ll12G1WrlpJNOYs6cOVgsFi699FJ+9KMf0dTUxNKlS8nLy+vX\nysWilhoREUkgibKi8/z58znttNOoqKjghhtuoLW1lSeeeIJ169aRkZEBQEVFBffccw+1tbXMnTuX\nRx55hMcee6zXDKeTTz6ZZcuWsXbtWkpLS5k5cybvv/8+JSUl/OlPf8LtdvP1r3+dO++8k9tvv52y\nsrJ4XHJKMMxEmfwfZbm5uezevTveZYjIIDBujs4vOnPpkHg7HBTt7e3U19czZcoU0tPTB/W1O6d7\nawZa6gvn+yiS399qqREREZGUoDE1IiKSUNRCI/2llhoRERFJCQo1IiIikhIUakRERCQlKNSIiEjE\nOqdaD5EJtBIlPffSGiiFGhERiVhaWhrp6elh7XEkEszhw4dJT08nLW1w4ohmP4mISL+MHTuWPXv2\nMGnSJDIyMhJiKwNJDqZpcvjwYfbs2cO4ceMG7XUVakQked1/5OO1ca1iyBozZgwA77//Pu3t7XGu\nRpJNeno648aN838fDQaFGhER6bcxY8YwZswYOjo6NL5GwmYYxqB1OXWlUCMiIgMWjV9QIpHSd6GI\niIikBIUaERERSQkKNSKSnEygDTgMuI/cF5EhTWNqRCT5HAD+eOQjwGNAFrDgyEcRGZLUUiMiScU0\nTV+gae5ysP3I/ZWoxUZkCFOoEZGkUltbC/voHV5MYD/QGPuaRCQxKNSISFJpaGgIfjKN7i04IjKk\nKNSISFJxOByQHuRkB5Ady2pEJJEo1IhIUnE6nb7BwD23GTKAMUBe7GsSkcSgUCMiScUwDN8sp64t\nMulADnAJvcOOiAwZmtItIsknC1gM3I1v5tP5+FpoFGhEhjSFGhFJTgZgPXKzx7kWEUkI6n4SERGR\nlKBQIyIiIikhrqHmhRdeYPr06UybNo2pU6fyhz/8AYC9e/dSWVlJUVERU6dOpbq62v+cUOdERERk\n6IrbmBrTNLnkkkv4xz/+wUknncTOnTs5/vjj+eY3v8n111/PrFmzWLduHRs3buTcc89lx44dWK3W\nkOdERERk6IrrQGHDMDhwwLcj3SeffEJOTg7Dhg1j9erV/lVDZ86cycSJE1m/fj1z5swJeU5Ehphr\n412AiCSSuIUawzBYtWoV3/zmNxk+fDj79+/n6aef5tNPP6WtrY1jjjnG/9jJkyfT2NjIvn37gp4T\nERGRoS1uY2q8Xi+/+MUvePrpp3G73bzyyissWLAAr9frW1yrC9P8fOe6UOe6Wr58Obm5uf7bwYMH\nB/8iRGRIKSkpoaSkJN5liEgQcWup2bx5M++//z5lZWWAryspNzeXf/3rX6Snp/PBBx/4W2Tcbjd5\neXnk5OQEPdfTkiVLWLJkif9+bm5uDK5KRJKZcXPg1fvMpYH/eBKRxBK3lppjjz2WpqYmtmzZAvh2\n3nW5XBx33HHMnz+fBx54AICNGzeyZ88eZs+eDRDynIiIiAxdcWupGT9+PA899BDf+ta3SEtLo6Oj\ng/vuu4+8vDzuuOMOFixYQFFRETabjZUrV/pnN4U6JyIiIkOXYQYblJJicnNz2b17d7zLEJFBEKyb\nKFo6u586x9PU1dXF9POLDGWR/P7W3k8iktBiHWCCMU2TlpYWPB4PNTU1OJ3OXhMXRCS+tE2CiEgf\n3G43xcXFuN1umpqaKC8v998XkcShUCMiEooJFRUVuFwuTNPENE08Hg8ul4vKysqgy0qISOwp1IiI\nhNII9fX1eL3eboe9Xi/bt2+ntrY2ToWJSE8KNSIioTQHP2Wz2fzbtohI/CnUiIiEkk3QZSM8Hg8O\nhyPGBYlIMAo1IpL47j9yi4c8yM/Px2LpPlnUYrFQUFCA0+mMU2Ei0pNCjYhIKAZUVVVRWFiIYRgY\nhoHNZsPhcFBVVaVp3SIJROvUiIj0wW63s2XLFgoKCvB4PKxatUrr1IgkIIUaEZEwGIZBZmYmmZmZ\nlJaWxrscEQlA3U8iIiKSEhRqRCSxmUAbcBhwH7kvIhKAup9EJGG53W64Dzhw5MBjQBaw4MhHEZEu\nFGpEJCGZpklFRUX3xe/a8d1fCVwLxHicrnbnFkls6n4SkYRUW1tLfX197+4mE9gPNMahKBFJaAo1\nIpKQQm4/kEbI7QtEZGhSqBGRhORwOIJuT0AHkB3TckQkCSjUiEhCcjqd5Ofn9x43YwBjgLw4FCUi\nCU2hRkQSkmEYVFVVdW+RSQdygEuI+SBhEUl8mv0kIgnLbrfDYuBufDOfzsfXQqNAIyIBKNSISGIz\nAOuRmz3OtYhIQlP3k4iIiKQEhRoRERFJCep+EpHEd228CxCRZKCWGhEREUkJCjUiIiKSEhRqRERE\nJCUo1IiIiEhKCBlq7r77bpqbtWuciIiIJL6QoebVV1/Fbrfz7W9/mxdffDFWNYmIiIhELGSoefbZ\nZ3G5XJx66qn88Ic/xG63s2zZMtxud6zqExEREQlLn2Nqxo0bx49//GPq6up4/PHH2bNnD6eccgpf\n/vKXY1GfiIiISFgiGihcVFREcXEx48ePx+VyRasmERERkYj1GWo6OjpYu3Yt5557Lg6Hg02bNnHv\nvfcq1IiIiEhCCblNwo9//GP+9Kc/MX78eK644goeeeQRxowZE6vaRERERMIWMtR89tlnrF27lhkz\nZsSqHhEREZF+Cdn99Nvf/tYfaJ566iluvfVWAPbs2cM777wT/epEREREwhTWQOGbbrqJ3//+96xY\nscL3pLQ0Fi1aFM26RERERCISVqh59tlnee655xg+fDgAEyZM4ODBg1EtTERERCQSYYWao446ivT0\n9G7HTNOMSkEiIiIi/RFyoHAnu93Oa6+9hmEYtLW1ceuttzJt2rRo1yYikvBKSkoAqKuri3MlIhJW\nqLnnnnu47LLLePfddxk+fDhnnnkmf/rTn6Jdm4iIiEjYwgo148ePZ926dbS0tGCapn9sjYiIiEii\nCBlq2tra+PTTT8nOzgagpqaGw4cPA3DGGWdoIT4RERFJGCFDzc9//nNGjRrFDTfcAMB3vvMdCgoK\nOHz4MBs2bOC2226LSZEiIonINE1aWlrweDzU1NTgdDoxDCPeZYkMWSFDzYsvvsg///lP//3s7Gz+\n/ve/09HRQVlZWdSLExFJVG63m4qKCtxuNwDl5eXk5+dTVVWF3W6Pc3UiQ1PIKd1paWkcddRR/vuX\nXXaZ/3hra2t0KxMRSVCmaVJRUYHL5cI0TUzTxOPx4HK5qKys1JIXInESsqXmwIED3e7/4Ac/8P97\n37590alIRCTBGDf36FJyw7Cdw/B6vd0Oe71etm/fTm1tLaWlpbErUESAPlpqZsyY4d8aoavHHnuM\n6dOnR6smEZHE1gxWqzXgKZvNRkNDQ4wLEhHoo6Xmtttuo7S0lBdffJEvfvGLGIbB66+/zquvvspr\nr70WqxpFZAjo1RqSyLLB4/EEPOXxeHA4HDEuSESgj5aawsJC3nrrLQoLC3nxxRepqqqioKCAt956\ni6KioljVKCKSWPIgPz8fi6X734UWi4WCggKcTmecChMZ2vpcfO+YY47hlltuiUUtIiLJwYCqqioq\nKiqor68HfN1RBQUFVFVVaVq3SJyEtaKwiIh0Z7fb2bJlCwUFBXg8HlatWqV1akTiTKFGRKSfDMMg\nMzOTzMxMzXYSSQAhx9SIiIiIJIuIQ83+/fujUYeIiIjIgIQMNXfffbf/3zt27KCkpISJEyeSn5/P\nO++8M+BP3trayuLFiykqKuLEE0/kkksuAWDbtm2cfvrpTJkyhZkzZ1JXV+d/TqhzIiKxVldXp/ch\nkQQRMtT84Q9/8P/7v//7v7nmmms4dOgQd911F0uWLBnwJ7/++usxDIP6+nreeecd7rrrLgAWLVrE\n1VdfTX19PT/96U9ZuHCh/zmhzolICrj/yE1EJEKGGWKTkunTp/Pmm28CcPLJJ/P222/7z02bNo3N\nmzf3+xN/9tlnTJgwgd27dzNq1Cj/8b179+JwOGhubsZisWCaJhMmTKCmpoZRo0YFPdfXYle5ubns\n3r273/WKSHT5F9/rDDTXxq2UsJhLtb+TSCxE8vu7z72f1q5di2matLW1dTs30A3bXC4X2dnZ3Hrr\nrbz88sscddRRLFu2jKysLCYl9ugaAAAgAElEQVRMmOBf1MowDPLy8mhsbGT06NFBz2kFTxERkaEt\nZPdTXl4ey5cv59e//jXjx49nz549gK81xWazDegTt7W14Xa7OeGEE9i0aRP33HMP3/72t3ttEAeh\nA1Swc8uXLyc3N9d/O3jw4IDqFRERkcQWsqXmH//4R8DjOTk5rF+/fkCf2G63k5aWxsUXXwzAKaec\nQn5+Pm63m6amJrxer7+LadeuXeTl5TFq1Kig53pasmRJt3E/ubm5A6pXRGLABNqAdsAN5AFay05E\nwtSvdWrS09PJzMwc0Cc++uijKS8vp6qqCvDNrtqxYwdOp5Pp06ezcuVKAJ566ilyc3NxOByMGzcu\n6DkRSXIHgPuOfPwUeKzLfRGRMIQcKBzK3Llzee655wb0ybdv386VV17JRx99RFpaGjfddBPnnXce\nW7duZeHChezbt49Ro0bx6KOPcuKJJwKEPBeKBgqLJC7TNEkbmwbN+FprOhlADr5BwwnWYqOBwiKx\nEcnv736HmjfeeIMZM2b056lxoVAjkrhqamo444wzAp9MBy4F7LGsqG8KNSKxEcnv735vk5BMgUZE\nEltDQ0Pwk2n4WnBERPrQ54aWf/vb33j88cdpbGwEfDOiLrjgAs4555yoFyciQ4PD4fC1yLQHONkB\nZMe4IBFJSiFDzc9//nOqqqq44ooruPDCCzFNE7fbzdKlS/nnP//JL37xi1jVKSIpzOl0QhaBx9SM\nwTcLSkSkDyHH1BQVFVFXV9drTZrW1lZKSkpCNxknGI2pEUlsxg8N+COw78iBdHyB5hJ8gSfBaEyN\nSGwM2orCpmnS0dHR63hHR8eAVxQWEekmC1gM3I2vG+p8tE6NiEQkZKhZuHAhM2fOZMGCBdjtvqkH\nbrebP/7xj1x++eUxKVBEhhADsB65JdhsJxFJfH2OqSkrK2P16tXU1NQAvoHC9913H7Nnz45JgSIi\nIiLh6HP2U1lZGWVlZbGoRURERKTf+gw1IiIxdW28CxCRZNXvxfdEREREEolCjYiIiKSEkKHmxBNP\n5IYbbmDTpk2xqkdERESkX0KGmldffZX8/HyWLVtGUVER3//+93nllVdobw+0lrmIiIhI/IQMNWPH\njuWqq67iueee46233qKsrIyHH36YoqIiLr30Up5++ulY1SkiIiISUthjakaMGMH8+fP585//zNat\nW7n44ot56aWXolmbiIiISNj6NaXbarVSUVFBRUXFYNcjIiIi0i+a/SQiIiIpQaFGREREUkLEoebj\njz/m3XffjUYtIiIiIv0WVqiprKzkwIEDHDx4kJNPPpm5c+dy0003Rbs2ERERkbCFFWr+85//kJWV\nxQsvvMC8efPYtm0bzzzzTLRrExEREQlbWKGmra0NgOrqas4++2ysVisWi/bCFBERkcQRVqiZOnUq\nX/nKV3juuec466yzaGlpwTTNaNcmIiIiErawmltWrFjBunXrOPnkk8nMzGTPnj3cfvvt0a5NRFKQ\ncbMR7xJEJEWF1VKzZMkSvvGNb5Cfnw/ApEmT+Otf/xrVwkREREQiEVaoef3113sd++c//znoxYiI\niIj0V8jup1WrVvHEE0+wY8cOvvnNb/qPf/zxxwwfPjzqxYmIiIiEK2SoOf7445k3bx5vvvkm8+bN\n8x8fNWoU5eXlUS9OREREJFwhQ83JJ5/MySefzFe/+lXGjh0bq5pERJJeSUkJAHV1dXGuRGToCGv2\nk9Vq5Y477sDlcuH1ev3HH3nkkagVJiIiIhKJsELN+eefz9ixYznttNNIT0+Pdk0iIgkv2NR0c6nW\n8BKJl7BCTVNTEy+//HK0axERSQmmadLS0oLH46Gmpgan04lhaH0ekWgLa0p3YWEhBw4ciHYtIjKU\n3H/klmLcbjfFxcW43W6ampooLy/33xeR6AqrpSYzM5Pp06dTWVlJRkaG//jy5cujVpiISNIxoaKi\nApfL5d9KxuPx4HK5qKys5L333lOLjUgUhRVqiouLKS4ujnYtIiLJrRHq6+t77Y3n9XrZvn07tbW1\nlJaWxqk4kdQXVqhZunRptOsQkaHEBNqAdsAN5AGp0IDRHPyUzWajoaFBoUYkisIaU7Nr1y7mzp3L\ntGnTANi8eTO//vWvo1qYiKSoA8B9Rz5+CjzW5X6yy/YtgRGIx+PB4XDEuCCRoSWsULNo0SIuuOAC\nf5Pq1KlTefjhh6NamIikHtM04Y90b9FoP3J/Jb4WnGSWB/n5+Vgs3RvBLRYLBQUFOJ3OOBUmMjSE\nFWr27t3LJZdcQlqa7+EWi6XXD62ISF9qa2thH73DiwnsBxpjX9OgMqCqqorCwkIMw8AwDGw2Gw6H\ng6qqKg0SFomysEKNxWLpNvBt//79UStIRFJXQ0ND8JNphByTkizsdjtbtmzBbrczYcIEXnnlFd57\n7z3y8vLiXZpIygsr1MyfP59FixbxySef8Pvf/56zzz6bq666Ktq1iUiKcTgcEGxR8g4gO5bVRI9h\nGGRmZpKVlUVpaalaaERiJKw+pB/96Ec8/vjjfPzxx7z44ossWbKEiy66KNq1iUiKcTqdkIWvRaZr\nF5QBjME3C0pEpJ/CCjXNzc1ceOGFXHjhhdGuR0RSmGEYsADfYOF9Rw6m4ws0l5Aa07pFJG7C6n4q\nKipi/vz5/O1vf+u1qJSISESygMVHPo4ELgWuPXI/hdTV1VFXVxfvMkSGlLBCTWNjI1/96le54447\nyMvL42c/+xn19fXRrk1EUpUBWIEMwI5aaERkUIQVaoYPH87ChQv5xz/+QXV1NR999JG2TRAREZGE\nEvZiM21tbfz1r3/lkUceYePGjXzve9+LZl0ikuqujXcBIpJqwgo11113HatWrWL69OksXLiQZ555\nBpvNFu3aRERERMIWVqgZP348b775JpMmTYp2PSIiIiL9EnJMTefI/RtuuIGjjz6627n169dHryoR\nERGRCIUMNQsWLPD/+7TTTut27oc//GF0KhIRERHph5ChpuuaND3Xp9F6NSIiIpJIQoaarvuV9Ny7\nRHuZiIiISCIJOVD40KFDvPPOO5im2e3fnedEREREEkWfoebrX/+6/37Xf6ulRkRERBJJyFCzc+fO\nGJUhIjI0lJSUAGhfKJEoCGubhGh79NFHMQyDv/zlLwDs3buXyspKioqKmDp1KtXV1f7HhjonIiIi\nQ1fY2yREi9vt5ne/+x2zZs3yH7v++uuZNWsW69atY+PGjZx77rns2LEDq9Ua8pyIJBbjZnVTd2Wa\nJi0tLXg8HmpqanA6nerKFxlEcW2p6ejo4Morr+Tee+9l2LBh/uOrV6/mu9/9LgAzZ85k4sSJ/sX+\nQp0TEUlUbreb4uJi3G43TU1NlJeX+++LyOCIa6hZvnw5TqeTGTNm+I/t27ePtrY2jjnmGP+xyZMn\n09jYGPKciEiiMk2TiooKXC4XpmlimiYejweXy0VlZaXW/RIZJHELNXV1dTz55JP8/Oc/73WuZ3Ns\n1x/4UOe6Wr58Obm5uf7bwYMHB6FqEZHI1dbWUl9fj9fr7Xbc6/Wyfft2amtr41SZSGqJW6iprq7G\n7XZTVFTE5MmTef3117n66qtZvXo16enpfPDBB/7Hut1u8vLyyMnJCXqupyVLlrB7927/bcSIETG5\nLhGRnhoaGoKes9lsIc+LSPjiFmq+973v0dTUxM6dO9m5cyezZs3ioYce4nvf+x7z58/ngQceAGDj\nxo3s2bOH2bNnA4Q8JyKSiBwOR9DJDB6PB4fDEeOKRFJT3Gc/BXLHHXewYMECioqKsNlsrFy50v+G\nEOqciEi8BZrx1XFTB/n5+bhcrm5dUBaLhYKCApxOZyxLFElZhjlERqjl5uaye/fueJchMqRoSreP\nudTE7XZTUVFBfX09AFarlYKCAqqqqgJ2oYuITyS/vxOypUZEJNXY7Xa2bNlCQUEBHo+HVatWaZ0a\nkUGmUCMiEiOGYZCZmUlmZialpaXxLkck5STENgkiIiIiA6WWGhGRGNJGliLRo5YaERERSQkKNSIS\nXfcfuYmIRJlCjYiIiKQEhRoRiR4TaAMOA+4j90VEokQDhUUkKtxuN9wHHDhy4DEgC1hw5KOIyCBT\nS42IDDrTNKmoqIDmLgfb8d1fiVpsRCQqFGpEZNDV1tayc+fO3uHFBPYDjbGvSURSn0KNiAy6hoaG\n4BvNptG9BUdEZJAo1IjIoHM4HHg8nsAnO4DsmJYjIkOEQo2IDDqn00l+fj703KvRAMYA2pRaRKJA\noUZEBp1hGFRVVXVvkUkHcoBL6B12REQGgaZ0i8igKSkpAXz7G9ntdlgM3I1v5tP5+FpoFGhEJEoU\nakQkegzAeuRmj3MtIpLyFGpEZMCMmw3fdO33gXYwrjA+b5W5Nr61icjQoVAjIgN3APgjWj1YROJK\nA4VFZEBM0/QFGq0eLCJxplAjIgNSW1sL+9DqwSISdwo1IjIgDQ0NwU9q9WARiSGFGhEZEIfD4VuD\nJhCtHiwiMaRQIyID4nQ6fYOBtXpwv5WUlPjX+BGR/lOoEZEBMQzDN8tJqweLSJwp1IjIwGXhWz04\nCxgJXIpvfRpN5+6TaZq0tLRw4MABampqfLPJRKRfFGpEZHB0rh6cgW/1YLXQ9MntdlNcXIzb7aap\nqYny8nL/fRGJnEKNiAyea9EKwmEyTZOKigpcLhemaWKaJh6PB5fLRWVlpVpsRPpBKwqLiESZcXOA\nZis3GPVGr/Di9XrZvn07tbW1lJaWxqhCkdSglhoRkXgIsX6PzWYLvf6PiASkUCMiEg/ZYLVaA57y\neDy+9X9EJCIKNSIi8ZAH+fn5WCzdRwFYLBYKCgp86/+ISEQUakRE4sGAqqoqCgsLMQwDwzCw2Ww4\nHA6qqqp86/+ISEQ0UFhEJE7sdjtbtmyhoKAAj8fDqlWrcDqdCjQi/aRQIyISR4ZhsGPHjniXIZIS\n1P0kIiIiKUGhRkQipg0YRSQRKdSISES0V5GIJCqFGhEJW7C9ijgQ78pERBRqRCRMnXsVba3f2m2v\noq31W2EloAYbEYkzhRoRCUttbS319fW9w4sJ7Aca41BUitPYJZHIKNSISFhC7kWURsi9jCRyGrsk\nEjmFGhEJi8PhCLpXER1AdkzLSWnBxi653e54lyaS0BRqRCQsTqeT/Px86LnYrQGMAfLiUFQK6hy7\n5HK5uo1dcrlcVFZWqsVGJASFGhEJi2EYVFVVdW+RSQdygEvoHXakXzrHLnm93m7HvV4v27dvp7a2\nNk6ViSQ+hRoRCZvdbofFQBYwErgUuPbIfRkUocYu2Wy20GObRIY47f0kIpExgB/Eu4jU1Tl2yePx\n9Drn8XhwOBxxqEokOailRkQkgXSOXbJYuv/NabFYKCgowOl0xqkykcSnUCMikkA6xy4VFhZiGAaG\nYWCz2XA4HFRVVWEYGrwkEoy6n0REEozdbmfLli0UFBTg8XhYtWoVTqdTgUakDwo1IiIJyDAMMjMz\nyczMpLS0NN7liCQFhRoRkQRVV1cX7xJEkorG1IiIiEhKUEuNiARk3KzxGyKSXBRqRCQ4E7gbaAfO\nx7cVgrLOoAkUHM2l2gZBpL8UakQksAPAH498BHgM38rBC9AKwiKSkOI2pubw4cN84xvfYMqUKZx8\n8smcffbZ/uW/9+7dS2VlJUVFRUydOpXq6mr/80KdE5HBYZqmL9A0dznYfuT+SnwtOCIiCSauA4Wv\nvvpqtm7dyttvv828efO46qqrALj++uuZNWsW27Zt49FHH+Wiiy6ira2tz3MiMjhqa2thH73Diwns\nBxpjX5OISF/iFmoyMjI455xz/ItJzZo1i507dwKwevVqvvvd7wIwc+ZMJk6cyPr16/s8JyKDI+Sm\niWl0b8EREUkQCTOl+5577mHevHns27ePtrY2jjnmGP+5yZMn09jYGPKciPSfaZrk5+czadIkampq\nKCwshPQgD+4AsmNZnXTq+XUyTfUDinSVEAOFb731VrZt28Yrr7zCoUOHei0F3vUHN9S5rpYvX87y\n5cv99w8ePDiIFYukDrfbTUVFBW63G4Dy8nImT54Mo/ANEu76I2YAY/DNgpKYCvR1ys/Pp6qqCrvd\nHufqRBJD3Ftq7rrrLp5++mn+9re/kZmZSU5ODunp6XzwwQf+x7jdbvLy8kKe62nJkiXs3r3bfxsx\nYkRMrkckmZimSUVFBS6XC9M0MU0Tj8fD9u3bfQ/o2iKTDuQAl6Bp3TEW7OvkcrmorKxUi43IEXEN\nNcuXL+fxxx/npZdeIivr8zmi8+fP54EHHgBg48aN7Nmzh9mzZ/d5TkQiU1tbS319PV6vt9txr9cL\nnwBfwzd9eyRwKXAtms4dB7W1tezcuTPg12n79u2+gd0iEr/up927d/OjH/2IgoICzjzzTACGDRvG\n//7v/3LHHXewYMECioqKsNlsrFy5EqvVChDynIhEps8BwfuBH8SqGgmmoaEBq9VKa2trr3M2m42G\nhgZteilCHENNbm5u0CbT8ePH8+KLL0Z8TkQi43A4sFqteDye3ic1IDhhOByOwF8jwOPx4HA4YlyR\nSGKK+5gaEYkfp9NJfn5+7zEyGhCcUDq/ThZL979DLRYLBQUFOJ3OOFUmklgUakSGMMMwqKqq0oDg\nBNf5dSosLMQwDAzDwGaz4XA4qKqq6jUrVGSoSogp3SISP3a7HRajjSsTnN1uZ8uWLRQUFODxeFi1\nahVOp1OBRqQLhRoR8QUYDQhOeIZhsGPHjniXIZKw1P0kIiIiKUEtNSIiCcS4OXB3krlUC+yJ9EUt\nNSJDiPYOEpFUplAjMkS43W6Ki4txu900NTVRXl7uvy8ikgoUakSGgL72DkINNiKSAhRqRIaAvvYO\nojFOhYmIDCKFGpEhoHPvoEBsNhs0x7ggEZEoUKgRGQJC7R108NBB7fEkIilBoUYkRXWd6dT5b+3x\nlJo0q03ER6FGJMkF+oXWc6bTnDlzaGtrg6wuT9QeTylBs9pEPmeYQyTS5+bmsnv37niXITKo3G43\nFRUV1NfXA2C1Wpk8eTJer5fGxsZuA4MtFgve0V74Or4xNNloj6ckEmjxPdM0KS4uxuVy9fpaOxwO\n3nvvPe0NJUkvkt/faqkRSVKhpmlv37494EwnDhy5cwpgR4EmydXW1lJfXx90VltNTY26pWRIUagR\nSVLBfqG1t7cHf1IamumUQhoaGoKes1gsfOtb31K3lAwp2vtJJImYpkltbS0NDQ3+LqeIdKCZTkkq\n0J5Qr5W/htVqDTizraWlhdbWVn/rTNfFFtUtJalKoUYkSfQcP5Oenh66O8Gg+0rBmumUcpxOJ/n5\n+b3G1KSnp9Pe3t6r1a6zW6q2tpbS0tJYlysSdep+EkkCgcbP9Ox26mSxWCgsLOzeIqOZTinJMAyq\nqqooLCzEZrMxYsQIbDYbY8eODdoSY7PZQnZbiSQztdSIJIHO8TOhWmYMw8BqtVJQUEBVVRX2R+y+\n7Q800yml2e12tmzZ4u+WdDgcmKbJnDlzAnZLeTweHA5HHCoViT6FGpEkEOov6+HDh3PdddcxZcoU\nHA4HTqfT91e6gW+Gkz1mZUqcGIZBaWmpv0upc+2iQFO9CwoKcDqd8SpVJKoUakSSgMPhCDogtK2t\njXPOOUdjJMSvs1uqoqKCHTt2YLPZ8Hg8FBQUsG7dum6tOv4QLJICtPieSBIItciaN8sL16KupSEo\n0IJ83c53mS3ncDjIzc2lsrKy22KN+fn5vu5Ku5r0JDFF8vtboUYkSQRaPbigoIB/V/67+/YH0t2y\nN4IcnxHbOqKgr1DT7bFafViSVCS/v9X9JJKATNOkoKAAj8fDqlWrcDqdAQeEOp1O0v6fJjH6BQsw\nKSrQ2jUQOOwEG2yuad6SShRqRBJMZ4tM58qv5eXl3boIug4IlQEIFIBSoPUmmFCDzTuneev7SpKd\nQo1IAum5Hg10Xwn239/+t8bOSL+EGmyuad6SKhRqRAZZz8GZoWaX9HysaZohuwhoRFO0pV+CrT6s\nad6SShRqRAZRoMG8wWaXBHpsdnbwjZlsNhueZo9CjfRLqGneVVVVGiQsKUGhRmSQ9NV1VFdXxz//\n+U8aGhooLCzkO9/5Tq/Hfvjhh0Ff/+Chg9qMUgYk2GBzwzAiamEUSVQKNSKDJNTsEpfLRX5+Prt2\n7QI+33Cw52P9GxBqM0oZgFCzonquPgyRtTBCZF2sIrGkUCPSxUDerEPNLvF6vezevdsfYoJtRgmQ\nmZlJS1oLHMK35WwHvkCjzSg/N8SmbkdTqBbGiooKfve73+Fyufw/D42Njb26sLSAnyQKhRqRIyL9\na7WnULNLOnfWDofX64WL8AUYbUYZ2wCTwgv1BROqhXHr1q2UlZX5N0udPHkyXq+XxsZGvF6v/3u9\ns4tVC/hJvGnVLhF6/7Vqmma38TDBAolpmtTU1LBixQr/JoIWS/e/FdLSwv8x65yJ4t+I8pQjH/V7\nQqIkVAtjp64/D9u3b+/V0th1AT+ReFKokZTSNWTU1NSE3TpSW1vLzp07I3qzdrvdFBcXU1ZWxhVX\nXMGcOXNoa2sjLy8Pm83GiBEjsNls5ObmYrVag3/yNMAKpIM3y+vb9kAhRmKks4UxHP4xXwF0LuAn\nEk/qfpKUMZDuo4aGBqxWK62trb3OBVptNdg4hMbGRgoLC/Fc7PFNv86GxmMb4X58XUk9B/9mA19H\n3UwSE4EGEHfc1BFw/ZpIaQE/SQQKNZKUeg7oPf300yMa7Niz39/hcAQcC9P5OoWFhdTU1IS1SN6O\nHTvgdHxdR50WAH8EDtB78G8WWnsGNPg3TtL+Xxp8he7fn+1gmEbQls7O2XudtICfJAqFGkkY4c48\nCtQiM378eP906a46BzuWl5czbNiwbjM18vLy/J+vsLAw6Gqrxx57LN/5znfCXiTPY3p8LS9dg0oW\nsBjfisBqlZFE0/P7cwxM2TAl4M+D3W7HYrFoAT9JSIYZ7qCDJBfJ1uUSe51Bpec00XXr1rF79+5u\nLTInnHBCrzfbtLQ0Ojo6wvpcXd+YuwaV3NxcAHbv3t2tBq/Xi9vt7vb50tPT6ejoCPyXbDpwKWp9\niVSytdSk8IwogJ0LdwZdffjYY4/VAn4SM5H8/laokbgzTZPi4uJeQSU9PZ309HQA/5tqsBaZ/gjU\nhF5YWNitq8o0TWbPnh18wHGgRfJygGtRK0ykFGoSirnUjCikDHRJBJFgIvn9re4niapAb4pAWJs4\ntre3+0NH53iXPXv2BP1chmFgGEbYLTY9Z3J0jocxDIOFCxcCsGLFiuAvYAEy0CJ5kpJ6DSp2g1ka\nfGmDUFuEaP0aiRWFGomaQF1Kgbp4Qo1P6SlUYLFYLEyYMIEPPvgAm83G4cOHA25FEErPmU6hFtTD\nBM5Hi+RFKtlaZKRPoRbw61wSoevsQZFoUahJIrHur47k84WajdR15dHt27f7n9N5LNQmjoEEapHp\n7DoKtmlkuONhPB4Pl1dfzuXuy49cGDCSwNOxx/D5wnhqXR96huDqw8GEWp/GarXywgsvaPyNxIRC\nTZxFOuOnr4G0/XnDCPTYSPZ3CVRbJGNfOruBeo5xCaZni0zXAYxpaWndNuurqqoKeB2BBv92Tkv9\nd96/P/9kBqGnY+s9WCRki+Znn33GXXfd1W324cMPP8yVV14Z9v5RCkASLg0UHmTBfvgGEhwiGUgb\n6RtGoEDSc3+XThaLBYfD0a1/PFhtkcxGAt8mjqNGjaK5uRmbzUZrayuGYXQbV9O1hq4tMl3/nwPu\nTmzSeyr1x4ReNyac19B7av+kevdTCrfUmEuDj6kJ9D4QSOf7VkdHR6/3l54D9ft6n+y6LIPCTurS\n7KcAYhFqgrWmBAoZkQSHmpoaysrKwhobEuoNI9xAEqrFxGq18pvf/IbMzMzwZgeFyWaz8fLLL2MY\nhv8NKjc3l8rKyoBTSvPy8gIHmEgoqERfqgeYQFI41ARjLjV7vf8dPny4XysUW61Wf6tOqPfJQOvl\naLZValKoCSDaoaav1pSeISNUcLDZbLzyyiv+LpQVK1ZwxRVXDHpwaGlpYfHixRG9rmEYpKenk5GR\n4R/k+5///Cfga4Q9G6nLNGhzWffXMU2TtCvTFDyS1VAMNYEMlaDT5Q+Fn53wM26//fYBvW/11SUd\nzrIMkXa9q6Un8WhKdxwE2xCxZ/dJ1+PBRDQDJwIWi4Vvfetb/i6eQ4cORfyGY5omXq+XgwcPAqEH\n+XaOfWnc0/h5F8+oIyc/IeD4lICtL507VkviUngR+HzQvB1uc9/m33Khv/oaYxdoWYZwVhAP1hLc\nV0tPIk/WEB+FmkESakPESLW2ttLS0sKKFSv8M4kGY8O5lpYWWltbaW9vDysghTNwN9gg366zkdKv\nSu/e0gLq9hFJdXn4xqf1nDkIvlW3OwIcHyRtbW20tbUB4HK5KC8v79ZV1fk+3dHR0e39MNS6OpFO\n1hioYJ9P3WuhqftpkNTU1FBeXh5xa0rPMBDu4N9gA2l9L0KvN4xwZxZB4P7qQ4cOBZ0G3XOQ76CP\nfZHEp5aa0IZK91NPBwg8IP/rwLN020CTEL3Ugd4ng25TEuZrBBNofF+w7VmCvV8PNHgEG84QaGxk\n1+ekaquOxtQEkGhjagIFh3Bm/HRr9RgFrKTvN4wO4CjgYIgLSOPzMNTZHTSaz1tU2oAqAjclp+Ob\n9qxF6IYGBZjoSrUAFGxAftfjY4C19GrVCTUgONCyDIMh0B9pkSxR0dcMzXCEmhzSc7JGpEtwQPIF\nIIWaAKIZag57D+Np99DobuSbX/sm7p1urFYrbW1tFBQU8PDDD3PFFVewY8cO3w7OPYLDa3Ne63vg\nbrBNEsN5w+hcsPcxAoeSNKASsBI8kJjAfQRehE57HaUmhZf4SLVQE65grTo9/8AKsiyDxbREvIJ4\nMD1bdSJdosJqtfZaSyuSFpxQk0N6TtYINUMsmlPkYxmMFGoCiFaoOew9zKTlk2g+1Ow7EG7ICBQc\n3sLXwhKIFTgHOKWfheZjyBwAABFCSURBVA5GKAn1phNobRdJbgo1iWUohJ1Il1kIo7Unkq6qSLrp\nQwm26nm4XUemaTJnzpxBGfs40CnygcYMRdoyNFAKNQFEK9R80voJo28fPTgv5iZ4a0qwlppIDEYo\n0douyUOhJLUMhVAzUIHe47KOfDxA4AHLEFY3faQb5gYSbtdRsPAxUJFOkQ80ZiiSNdYGi0JNAEkR\namLRxaNQknoUXoYGhZrwRLKC+MX4lpcIo5u+Z5fSwZaDvvfOHpMyQnVVRdJ1FMlkjViKZI21wZLy\n69Rs27aNyy67jI8++ojRo0ezYsUKSkpK4l3WwMVin6Eu60hIElKAGboG42s/FIJRoPe4LGAxgf+g\nG9PlsSaBp6Eb0Da6jcaFjbALPM2eoJM1OkZ0wKcEnM3Vc50vl8sVMCB4vV5cO12+3wenH/l8oSZr\nMHhdZ32JZI21eEjKULNo0SKuvvpqFi5cyJNPPsnChQvZuHFjvMsaHKF++OJBOxGHppAhySSS79dU\n+xkP5w+6vv6wTAvwGj3fr48F7ifw+jw9hAwhacB+fOMo7Ude638DvK4BZEF7WnvYU+T94zwHmcfj\nweFwDP4LRyDpup/27t2Lw+GgubkZi8WCaZpMmDCBmpqakP+ZUe1+umUstFu7n0hvA+vAVgCOuYG+\n4UUSgGIdlhQ+RGIr2M9yoJ/FwXjsYBpoN33PsT1eIg8RgcZRhjtDLMig6c4ARBrhjzsKVl6ABVc1\npqYf3njjDS666CK2bt3qP/aFL3yB22+/nbPOOivo86IVavZ+fIDxE96CQ4M0rkZERFKACWwGdsE5\nN4fsOvo8XXRKx9fk8ySB19c48rqX3xw8cA1kinz7Ufj6uzo/X/fwEmimVNcFVwdbyo+p6SlQLlu+\nfDnLly/33//ggw/Izc2NyueflP35vw8ePMiIESOi8nniTdeWvFL5+nRtySmVrw16XN/bk/iP8R+8\n9J7JlJ6ejmEYeL0dGIaBaZpYLAZHH32Y9PR5oT/Ji5NCn88AT7YHr9eLxWLBZrPBmnAfP4L09HQ+\n+ugjvF78tYGvx2TUqFGMHj0ar9fLUUcdxaeffsrpp58eup5+CrXHYE9J11LT3+6nWIn2ysXxpGtL\nXql8fbq25JTK1wapfX2JfG1p8S4gUuPGjWP69OmsXLkSgKeeeorc3NyECDQiIiISP0nZ/fTggw+y\ncOFCbr31VkaNGsWjjz4a75JEREQkztKXLVu2LN5FROroo4/mqquu4vvf/z5XX30148ePj3dJ3Zx2\n2mnxLiFqdG3JK5WvT9eWnFL52iC1ry9Rry3pxtSIiIiIBJJ0Y2pEREREAlGoERERkZSgUNMP+/bt\nY9q0af7blClTsFgsNDc3s3fvXiorKykqKmLq1KlUV1fHu9yIhbq+L33pS+Tn5/vP/frXv453uRF7\n4YUXmD59OtOmTWPq1Kn84Q9/AEiJr12wa0uFrxvAunXrOPXUUznppJOYNWsWb7/9NpAaX7tg15aM\nX7vrrruOyZMnYxgGmzdv9h/ftm0bp59+OlOmTGHmzJnU1dWFdS6R9OfaJk+ezHHHHef/Gq5atSoe\npYcl2PUFOw4J9rUzZcB+9atfmXPnzjVN0zQvv/xyc+nSpaZpmub//d//mZMmTTI9Hk8cqxu4rtc3\ne/Zs85lnnolzRf3X0dFhjhkzxnz77bdN0zTNHTt2mMOGDTM/+eSTpP/ahbq2ZP+6maZpNjc3m9nZ\n2ea7775rmqZpVldXmyUlJaZpJv/PXahrS8av3fr1681du3aZdrvdfOutt/zHzzzzTPPRRx81TdM0\n16xZY5566qlhnUsk/bm2no9NZMGuL9hx00ysr51CzSAoLi72v+kMHz7cbGpq8p+bOXOm+dJLL8Wr\ntEHR9fqS8Q22q46ODjM7O9tcv369aZqm+fbbb5sTJ040W1tbk/5rF+rakv3rZpqmuXHjRrOoqKjb\nsZEjR/7/9u4/Jur6jwP4k2FaeIQEasVI8o4Q7icepVLoZMNla1KOsKaQSsu5tbl+bMI/LSfRP65Y\na6NaYC0tkICl06Vhba1MvcoLSCTAIzAMEY8Ckl/y/P7h+OzgfvDjiOPu+3psbH7e73t/Xq/3vdju\n5eeO+/Dnn3/2+9p52ps/187xBbCjo4OhoaEcGhoiefv3denSpWxsbPQ4N1dNdm/jH+sv3OU8fnyu\n1U7efvLSmTNncOPGDTz55JPo6urC0NAQ7r33XmU+JiYGra2tPszQO477G5Wbmwu9Xo8tW7bg8uXL\nPsxu6oKCglBWVobNmzdj2bJleOyxx/DJJ5+gp6fH72vnbm/z588H4N91A4DY2Fh0dXXhzJkzAICj\nR4+ip6cHNpvN72vnbm8tLS0A/L92ANDW1ob77rsP8+bd/nq0oKAgPPDAA2htbfU45w8mk392djb0\nej1ycnKm9LX/c91cq500NV4qLi5Gdnb2mII6op//xfz4/X366ae4dOkSampqkJKSMqbZ8QfDw8PI\nz89HZWUl/vjjD5w+fRpZWVkYHh72+9q529v169f9vm4AEBYWhi+++AJ5eXkwm804deoUEhIS0Nvb\n6/e1c7e3O+64IyBq546nOvlbDcdzzP+7775DTU0NfvnlF0RGRuL555/3YWb/PZ/WzifXhwJET08P\nVSoV6+vrlbGQkBC/vgzuyNX+xluwYAGvX78+i1l5x9Vl/qSkJJ46dcrva+dpb+P5W91c6e/v56JF\ni9jY2Oj3tRvPcW/j+VPt5O0n5/zb29upUqlmNdfpkLef/g+VlZXBaDRixYoVytgzzzyD999/HwBg\nsVjw559/Yt26db5K0Svj9zc8PIyOjg5lvqKiAkuXLkVERISvUpyy6OhoXL16FfX19QCApqYmNDc3\nIy4uzu9r52lv/l63UVevXlX+vX//fqSmpkKj0fh97QDXe4uJiQmY2nm6b5+/39PPU/59fX3o7u5W\nHvv5558jMTHRV6nOuDlXO5+0UgFizZo1LCkpGTP2119/MS0tjRqNhgkJCfzmm298lJ33xu+vt7eX\nZrOZOp2OBoOBqamptFqtPsxwej777DNlDzqdjocPHyYZGLVztbdAqRtJ5uTkMC4ujmq1mtu2baPd\nbicZGLVztTd/rd2LL77IqKgoBgcHc8mSJVSr1STJS5cucfXq1YyNjaXZbGZNTY2yxtPcXDLVvTU3\nN9NkMlGv11On03HTpk202Ww+3IFn7vbnbpycW7WT2yQIIYQQIiDI209CCCGECAjS1AghhBAiIEhT\nI4QQQoiAIE2NEEIIIQKCNDVCCCGECAjS1AghhBAiIEhTI4QQQoiAIE2NENMUExODuLg4GI1GaDQa\npKenKzckBACTyYSenh4fZuiaL/KaiZhvvPEG+vv7Zygj38vIyMDZs2e9OgdJpKSkwGazTTrO0NAQ\n9u3bhxUrVkCr1SIxMRFPPfUUrFarx1gbN27Ee++95zRuNBpRVVU1qVyE+M/57Gv/hPBz4++BUlFR\nwbCwMJ49e9aHWc0to/eDmQkAlG8Q9nfnzp1jamrqjJyrqqqKWVlZk46zdetWpqen88aNG8rY119/\nzdLSUo9xysvLuXLlyjFjFouFS5Ys4eDg4IS5CDEbpKkRYppc3fBt7969zMjIIDn2RRgA8/Pz+cgj\nj3DZsmWsqqpiQUEBzWYzNRoNv/32W+Uc41+8IyIilK9V/+qrr5iYmEi9Xs+1a9fyt99+G7PuzTff\n5MMPP8yYmBinW3i4O7+7dfv37+dLL72kPK6np4fh4eHs7Ozk1q1baTabqdfr+cQTT4y5mSQAvv76\n60xKSuJrr73mFNPdWk/579q1iwCo0+loNBrZ0dHhtKfJPr/u4v/777/MzMxkfHw8DQYD09LSPI67\nel4nm8POnTv50UcfKcfbt2/nhx9+qBzn5eWxoKBgwjmSHBwcZGRkJLu7u51yGh/n999/Z0hICLu6\nulzugSTPnz/P9evX02w202Qy8ciRIyTJgYEBLl68eMxtGnbv3q3UeKJchJgN0tQIMU2umprKykrG\nx8eTdG5qCgsLSZLV1dVcuHAhDx48SJI8cuQIk5KSlHO4a2o6Ojp4zz33KPdVOXToEOPj4zkyMqKs\nO3DgAEmyvr6eKpXK5ZUSV02Nq3VtbW2MjIxkf38/SbKkpISbN28mSXZ2dirr33rrLe7atWvM+fbt\n2+c2pru1E+Xv6UrNVJ5fd/ErKyu5YcMGZW70hd/duDc5LF++nLW1tcqxTqcb87uUlpbGkydPTjg3\nav369Tx27JhTTuPjlJWV0WAwuMyfJO12O00mE9vb20nefq6io6N55coVkuTLL7/MPXv2kCRv3rzJ\n8PBwXrx4cVK5CDEbpKkRYppcNTUVFRVum5rRKwLd3d0EwJs3b5IkW1paGBYWppzDXVNz9OhRrlu3\nbky8sLAwtrW1OcUgyUWLFilzjlw1Ne7WbdiwQfmf+tq1a5UXq8LCQprNZmq1WqrVaq5atWrM+cbH\ndYzpbu1E+U/U1Ez2+XUXv7m5mdHR0dy9ezdLS0v5zz//eBz3Jof58+fz2rVrJMm+vj6nBi4iIoKd\nnZ0e5xw999xzLCoqcsrJMQ55u6kxGo3KcVNTE41GIx966CFu376dx48f5913302j0aj8REdH8/Tp\n0yTJuro6RkZGcmBggIcPH2ZycrJTTHe5CDEb5IPCQswgi8UCnU7ncu7OO+8EAAQHBzsdDw8PK48L\nDg7GrVu3lOPRD8eSRFBQkMf4o+d0dd7prNuxYwcOHjyIy5cvo6mpCRs3bsT333+Pd999FydOnEBd\nXR3efvttpw/wqlQql3EmWjvd/B3Xenp+PcVfvnw5Ll68iMcffxw//PADdDod7Ha72/Hp5gAAISEh\nStwLFy5Aq9Vi3rx5AICGhgaoVCpERkZ6nHPU39+Pu+66yykfxzgAkJiYiMbGRiV/tVoNq9WKvLw8\n2O12kIRWq4XValV+WltbkZqaCgDQarVQq9U4duwYSkpKkJOT4xTTXS5CzAZpaoSYIV9++SWKiorw\nyiuveHUetVqNc+fOAQAqKyvR19cHAFizZg1qa2tRV1cHACgtLUVUVBSioqK8S9yDp59+GhaLBQUF\nBdi2bRuCg4Nht9sRGhqKiIgIDA4O4oMPPpj0+bxZGxoair///ns625hU/CtXriAoKAibNm3CgQMH\nQBJtbW1ux71hMBjQ0NAAAPjpp58wMDCAkZERDA0NIT8/HyaTacI5R/X19TAajR7jAEBsbCzS09OR\nk5OD7u5uZXz0dyw5ORk2mw3V1dXKnNVqxeDgoHKck5ODgoICWCwWZGZmTjoXIWaDNDVCeGHLli3K\nn3QXFxfjxIkTWL16tVfnLCwsxJ49e7By5UpcuHABERERAIDFixfj0KFDyM7OhsFgQFFREcrLyye8\neuONBQsWIDMzE8XFxdixYweA23/aGxcXh7i4OKSkpLh8kXXHm7Wvvvoq0tLSYDKZcO3atSnvZaL4\ntbW1SE5OhtFoRGJiIrKysmAwGNyOeyMjIwMnT54EcLtxefDBB2EymZCSkoLo6GicP38ev/76q8e5\nUS0tLbh165bLRsIxzqiPP/4Yer0eq1atQkJCAh599FFUV1dj7969CA8Px/Hjx1FQUACj0YiEhATk\n5uZiZGREWf/ss8+ioaEBGRkZTlfkPOUixGwIIklfJyGEEP9Pent7kZycjB9//BFJSUkoLy93+bZl\nfHy827lRubm50Gg0eOGFFzzGWbhw4YzuYaq5CDEb5EqNEELMMpVKhXfeeQc2mw3t7e1ISEhwekxv\nb6/bOUf3338/du7cOWGc2eApFyFmg1ypEUIIIURAkCs1QgghhAgI0tQIIYQQIiBIUyOEEEKIgCBN\njRBCCCECgjQ1QgghhAgI0tQIIYQQIiBIUyOEEEKIgCBNjRBCCCECgjQ1QgghhAgI/wPOCbxR4ZSt\nJAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fd0a7f30890>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "VARIABLE = 'InvMass'\n",
    "\n",
    "fig = plt.figure(num=None, figsize=(8, 6), dpi=80, facecolor='w', edgecolor='k')\n",
    "\n",
    "aHisto   = bulkHisto(\"SingleMuon\")(VARIABLE)\n",
    "nBins    = len(aHisto.values)\n",
    "edges    = np.linspace(aHisto.low, aHisto.high, nBins + 1)\n",
    "width    = (aHisto.high - aHisto.low) / nBins\n",
    "\n",
    "plotVals = {}\n",
    "for k in bulkHisto.bins:\n",
    "    if k == 'SingleMuon':\n",
    "        continue\n",
    "    plotVals[k] = [x.toJson()['data']*0.19 for x in bulkHisto(k)(VARIABLE).values]\n",
    "    plt.bar(edges[:-1], plotVals[k], width=width, label=k, color=samples[k]['color'], edgecolor=samples[k]['color'], fill=True)\n",
    "    \n",
    "xdata   = np.linspace(aHisto.low+0.5*width, aHisto.high+0.5*width, nBins)    \n",
    "ydata   = [x.toJson()['data'] for x in bulkHisto('SingleMuon')(VARIABLE).values]\n",
    "yerror  = [x**0.5 for x in ydata]\n",
    "\n",
    "plt.errorbar(xdata, ydata, fmt='ko', label=\"Data\", xerr=width/2, yerr=yerror, ecolor='black')\n",
    "\n",
    "plt.xlabel('Dimuon invariant mass m($\\mu\\mu$) (GeV)')\n",
    "plt.ylabel('Events / 0.5 GeV')\n",
    "plt.legend(loc='upper right', fontsize='x-large', )\n"
   ]
  }
 ],
 "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.13"
  },
  "sparkconnect": {
   "bundled_options": [],
   "list_of_options": [
    {
     "name": "spark.executor.extraJavaOptions",
     "value": "-Djava.library.path=/usr/lib/hadoop/lib/native/lib"
    },
    {
     "name": "spark.jars.packages",
     "value": "org.diana-hep:spark-root_2.11:0.1.15,org.diana-hep:histogrammar-sparksql_2.11:1.0.4"
    },
    {
     "name": "spark.yarn.dist.files",
     "value": "{KRB5CCNAME}#krbcache,/eos/project/s/swan/public/PhysicsAnalysis/histogrammar.zip"
    },
    {
     "name": "spark.executorEnv.KRB5CCNAME",
     "value": "./krbcache"
    },
    {
     "name": "spark.driver.extraLibraryPath",
     "value": "/eos/project/s/swan/public/PhysicsAnalysis/lib"
    },
    {
     "name": "spark.driver.extraClassPath",
     "value": "/eos/project/s/swan/public/PhysicsAnalysis/EOSfs.jar"
    }
   ]
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}