TTTT Analysis

This Notebook is simply a playground to examine the resulting histograms from the main TTTT analysis executable.

In [8]:
import ROOT
from IPython.display import Image
%load_ext autoreload
%autoreload 2
from utils import HistCollection as HC
from utils import show_event

The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

First, we need to load(and optionally rebuild) the histogram datafiles. These will generally contain a set of histograms of various quantities calculated from data in the input MiniTrees.

In [9]:
rebuild_hists=False
hists_TTTT = HC("TTTT", "../data/TTTT_ext_treeProducerSusyMultilepton_tree.root", rebuild_hists=rebuild_hists)
hists_TTZ = HC("TTZ", "../data/TTZToLLNuNu_treeProducerSusyMultilepton_tree.root", rebuild_hists=rebuild_hists)
hists_TTW = HC("TTW", "../data/TTWToLNu_treeProducerSusyMultilepton_tree.root", rebuild_hists=rebuild_hists)
In [10]:
HC.canvas.Clear()
HC.stack_hist_array(*zip(('jet_count_os_dilepton','Jet Multiplicity for Opposite-Sign Dilepton Events'),
                         ('jet_count_ss_dilepton','Jet Multiplicity for Same-Sign Dilepton Events'),
                         ('jet_count_trilepton', 'Jet Multiplicity for Trilepton Events')
                        ),
                    normalize_to=1,
                    enable_fill=True,
                    shape=(3,1),
                   )
                    
HC.canvas.Draw()

We can use the show_event function to look at the Generator-Level particles for the event. They are color-coded based on their pt relative to the maximum pt of a particles in the event. Green means higher, red means lower

In [11]:
show_event(hists_TTZ, 3)
Out[11]:
In [12]:
HC.stack_hist("lepton_count", title="Lepton Multiplicity",
              enable_fill=True, normalize_to=1, make_legend=True, draw=True)
Out[12]:
<ROOT.THStack object ("lepton_count_stack") at 0x23712980>
In [16]:
HC.stack_hist("b_jet_count", title="B-Jet Multiplicity",
              enable_fill=True, normalize_to=1, make_legend=True, draw=True)
Out[16]:
<ROOT.THStack object ("b_jet_count_stack") at 0x233e2b50>

The draw method of the HistCollection class simply creates a grid of plots showing all of the drawable objects(i.e. objects with a Draw method) contained in the input file.

In [13]:
hists_TTTT.draw()

Warning in <TROOT::Append>: Replacing existing TH1: TTTT:nLepvsnJet (Potential memory leak).
In [14]:
hists_TTZ.draw()

Warning in <TROOT::Append>: Replacing existing TH1: TTZ:nLepvsnJet (Potential memory leak).
In [15]:
hists_TTW.draw()

Warning in <TROOT::Append>: Replacing existing TH1: TTW:nLepvsnJet (Potential memory leak).