TTTT Analysis¶
This Notebook is simply a playground to examine the resulting histograms from the main TTTT analysis executable.
In [1]:
import ROOT
from IPython.display import display, display_markdown
%load_ext autoreload
%autoreload 2
from utils import HistCollection as HC
from utils import show_event, CANVAS
First, we need to load the pre-processed datafiles. These will generally contain a set of histograms of various quantities calculated from data in the input MiniTrees. However, they can also contain things besides histograms. For example, C++ STL containers can be serialized to the ROOT file to save things such as counters or even "raw" event information.
In [2]:
hists_TTZ = HC("TTZ", "../data/TTZToLLNuNu_treeProducerSusyMultilepton_tree.root")
hists_TTW = HC("TTW", "../data/TTWToLNu_treeProducerSusyMultilepton_tree.root")
hists_TTTT = HC("TTTT", "../data/TTTT_ext_treeProducerSusyMultilepton_tree.root")
In [3]:
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),
)
CANVAS.Draw()
In [4]:
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),
draw_option='nostack',
)
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
- TTZ Generator-Level Particles
- TTW Generator-Level Particles
- TTTT Generator-Level Particles
In [5]:
HC.stack_hist("lepton_count", title="Lepton Multiplicity",
enable_fill=True, normalize_to=1, make_legend=True, draw=True)
Out[5]:
In [6]:
HC.stack_hist("b_jet_count", title="B-Jet Multiplicity",
enable_fill=True, normalize_to=1, make_legend=True, draw=True)
Out[6]:
