ratio_types_base = [
["Base: Raw $H_T$ ratio" , "kinem_base_ht_raw__comparison_mc.pdf"],
["Base: Raw $\\met$ ratio" , "kinem_base_rawmet__comparison_mc.pdf"],
["Base: Raw Njets ratio" , "kinem_base_njets_raw__comparison_mc.pdf"],
["Base: Raw Nbtags ratio" , "kinem_base_nbtags_raw_comparison_mc.pdf"],
]
sm.addObject("left7","text",(0.27,0.91),width=0.3, text="74X", color="red", size=1, bold=True)
sm.addObject("right7","text",(0.74,0.91),width=0.3, text="76X", color="red", size=1, bold=True)
sm.addObject("leftDMC","text",(0.27,0.91),width=0.3, text="MC", color="red", size=1, bold=True)
sm.addObject("rightDMC","text",(0.74,0.91),width=0.3, text="Data", color="red", size=1, bold=True)
sm.addSlide(title="Overview", text="""
- Kinematic plots are of events falling into HH, HL, or LL SRs
- 74/76 ratio plots at end
""")
for label, plot74, plot76 in plot_types:
sm.addSlide(title=label, p1=plot74, p2=plot76, objects=["left7","right7"])
for label, ratioplot in ratio_types:
sm.addSlide(title=label, p1=ratioplot, p2=ratioplot.replace("_mc.pdf","_data.pdf"), objects=["leftDMC","rightDMC"])
sm.addSlide(title="base", text="- base = all hyp\\_class==3 events (no requirement that it falls into a SR)")
for label, ratioplot in ratio_types_base:
sm.addSlide(title=label, p1=ratioplot, p2=ratioplot.replace("_mc.pdf","_data.pdf"), objects=["leftDMC","rightDMC"])
sm.writeSlides("SScomparison7476_Apr25.tex", opts="--compile --copy --dump")
import slideMaker as sm
import fnmatch
import os
import utils
#PATH=/nfs-7/texlive/2015/bin/x86_64-linux:$PATH
#export PYTHONPATH=$PYTHONPATH:/home/users/cgeorge/software/slideMaker/
t = "madrid"
exprs = {'_mm_OS', '_ee_OS', '_em_OS', '_all_SF', '_all_EZ', '_all_SS'}
#exprs = {'_all_SF','_all_EZ','_all_SS'}
for expr in exprs:
sm.initSlides(me="GC", themeName=t, opts="--themecolor 51,51,179")
for file in os.listdir('.'):
if fnmatch.fnmatch(file, '*' + expr + '.pdf'):
print file
sm.addSlide(p1=file)
sm.writeSlides("kinem_plot_dump" + expr + ".tex", opts="--compile --copy")
sm.addObject("t1","text",(0.25,0.15),width=0.3, text="testlabel", color="red", size=0, bold=False,opts="--rotate -45")
sm.addObject("t2","text",(0.75,0.15),width=0.3, text="testlabel", color="coolblue", size=0, bold=False)
sm.addObject("a1","arrow", (0.31,0.15), (0.69,0.15))
sm.addObject("a2","arrow", (0.31,0.15), (0.65,0.46), opts="--crayon --noarrowhead")
sm.addObject("a0","arrow")
sm.addObject("c0","circle")
# supported themes are "nick", "alex", and "madrid"
# test2 and test3 can be folders with your plots, so you can execute this script anywhere really
# themecolor takes RGB. could put in explicit color names, but RGB more robust
# also try the --modernfont option
sm.initSlides(me="Nick",themeName="nick",opts="--graphicspaths ./test2/,./test3/ --themecolor 51,51,179")
sm.addSlide(title="Perturbation Theory on $H_m(dS_n,\\mathbb{R})$ Orbifolds", opts="--shorttitle snarxiv hep-th")
# pass in a list of textobjects (which are just dicts, so they can be modified too)
sm.addSlide(p1="yields.pdf",p2="yields.pdf", objects=["t1","t2","a1","a2"])
sm.addSlide(p1="yields.pdf",p2="yields.pdf")
# when I specify an empty object, a helper grid gets printed on the slide if I've specified
# the global option "--makegrid". A gui is produced if I specify "--makegui" instead
sm.addSlide(p1="zmass.pdf", objects=["a0","c0"])
sm.addSlide(text=content+content)
# slides reset their numbering when you start the backup section
sm.startBackup()
sm.addSlide(text=content, p1="filt.pdf", opts="--sidebyside")
sm.addSlide(text=content, p1="zmass.pdf", p2="zmass.pdf")
sm.writeSlides("test_1.tex", opts="--compile --copy")
-- Pair of leptons has $|m_{\\ell\\ell}-m_Z|\\leq15$
-- Njets$<5$
-- 0 btags
-- in the HighHigh region (each lepton has $\\pt>25$)
""")
sm.addSlide(p1="h1D_njets.pdf", p2="h1D_nbtags.pdf", title="N-1 plots: njets, nbtags")
sm.addSlide(p1="h1D_lep1pt.pdf", p2="h1D_lep2pt.pdf", title="N-1 plots: lepton $\\pt$")
sm.addSlide(p1="h1D_met.pdf", p2="h1D_mtmin.pdf", title="CR plots: $\\met$, $\\mtmin$")
sm.addSlide(p1="h1D_ht.pdf", p2="h1D_zmass.pdf", title="CR plots: $\\Ht$, $m_{\\ell\\ell}$")
# DATA
# data = 14.00 pm 3.74
# WZ
# wz = 8.44 pm 0.13
# NONWZ
# non-wz = 1.91 pm 2.03
# SCALE FACTOR
# data/MC = 1.433 pm 0.505
sm.addSlide(title="CR scale factor",text="""
-- Put 100%% systematic uncertainty on nonWZ contributions
-- From this CR, after subtracting out nonWZ contributions, I find a WZ data/MC scale factor of $\\mathbf{1.45\\pm 0.50}$.
""")
sm.writeSlides("WZ101915.tex", opts="--compile --copy --dump")
- Using latest json with 1.264 fb$^{-1}$
- Construct \\coolblue{WZ control region} with baseline requirement of $\\met > 30$ and Njets$>2$ and
-- 3 tight leptons
-- Pair of leptons has $|m_{\\ell\\ell}-m_Z|\\leq15$
-- $2<$Njets$<5$
-- 0 btags
-- in the HighHigh region (each lepton has $\\pt>25$)
""")
sm.addSlide(p1="h1D_njets.pdf",
p2="h1D_nbtags.pdf",
title="N-1 plots: njets, nbtags")
sm.addSlide(p1="h1D_lep1pt.pdf",
p2="h1D_lep2pt.pdf",
title="N-1 plots: lepton $\\pt$")
sm.addSlide(p1="h1D_met.pdf",
p2="h1D_mtmin.pdf",
title="CR plots: $\\met$, $\\mtmin$")
sm.addSlide(p1="h1D_ht.pdf",
p2="h1D_zmass.pdf",
title="CR plots: $\\Ht$, $m_{\\ell\\ell}$")
sm.addSlide(title="CR scale factor",
text="""
-- Put 100%% systematic uncertainty on nonWZ contributions
-- From this CR, after subtracting out nonWZ contributions, I find a WZ data/MC scale factor of $\\mathbf{1.340 \\pm 0.435}$.
""")
sm.writeSlides("WZ062316.tex", opts="--compile --copy --dump")
sm.startBackup()
sm.addSlide(title="L1 Parameters",
p1="../../l1/parameters/Parameter_%s.pdf" % alg,
opts="--fithorizontal",
objects=["t" + alg],
text="""
- Ignoring binning issues, parameters look reasonable
""")
sm.addSlide(title="L1 $\\chi^2$/ndof",
p1=alg + "_L1chisq.pdf",
objects=["t" + alg],
text="""
- Find reduced $\\chi^2$ for each $\\eta$ slice fit
- The threshold for a ``bad'' fit is 9, so these are acceptable
""")
sm.addSlide(title="L2 $\\chi^2$/ndof",
p1=alg + "_L2chisq.pdf",
objects=["t" + alg],
text="""
""")
# sm.addSlide(title="L2L3 Fit Results", p1="L2AbsCorGraphs_%s.pdf" % alg, objects=["t"+alg], text="""
# - Fit correction as function of $\\pt$ for each of the 82 $\\eta$ bins
# """)
sm.writeSlides("jecValidation2016June17_v2.tex",
opts="--compile --copy --dump")
- Recommend using AK8 jets above 60 GeV
""",
objects=["tbottomak4pfchs","tbottomak4pf", "tbottomak8pfchs", "tbottomak8pf"]
)
sm.addSlide(title="Conclusions", text="""
- Distributions look acceptable after applying L1L2L3 corrections
- Text files in \\url{http://uaf-6.t2.ucsd.edu/~namin/dump/jecs_FastSim_80X_V1/}
""")
sm.startBackup()
sm.addSlide(title="L1 Parameters", p1="../../l1/parameters/Parameter_%s.pdf" % alg, opts="--fithorizontal", objects=["t"+alg], text="""
- Ignoring binning issues, parameters look reasonable
""")
sm.addSlide(title="L1 $\\chi^2$/ndof", p1=alg+"_L1chisq.pdf", objects=["t"+alg], text="""
- Find reduced $\\chi^2$ for each $\\eta$ slice fit
- The threshold for a ``bad'' fit is 9, so these are acceptable
""")
sm.addSlide(title="L2 $\\chi^2$/ndof", p1=alg+"_L2chisq.pdf", objects=["t"+alg], text="""
""")
# sm.addSlide(title="L2L3 Fit Results", p1="L2AbsCorGraphs_%s.pdf" % alg, objects=["t"+alg], text="""
# - Fit correction as function of $\\pt$ for each of the 82 $\\eta$ bins
# """)
sm.writeSlides("jecValidation2016June17_v2.tex", opts="--compile --copy --dump")
sm.addSlide(p1="tables/table1.pdf",
p2="zmass.pdf",
opts="--vertical --sizeratio 0.2")
### END TABLE EXAMPLE ###
# four plots with text
sm.addSlide(p1="tables/table1.pdf",
p2="tables/table1.pdf",
p3="tables/table1.pdf",
p4="tables/table1.pdf",
text="""
--four plots here...
--another bullet
""")
sm.addSlide(text1="""
- this is text that goes in the left column
""",
text2="""
- this is text that goes in the right column
""",
opts="--texttop")
# slides reset their numbering when you start the backup section if
# you specify opts="--resetnumbering"
sm.startBackup()
sm.addSlide(text=content, p1="filt.pdf", opts="--sidebyside")
sm.addSlide(text=content, p1="zmass.pdf", p2="zmass.pdf")
sm.writeSlides("test_1.tex", opts="--compile --copy --dump")
p2="../met3p0_OS.pdf",
text="")
sm.addSlide(title="$M_T$ min",
p1="../mtmin.pdf",
p2="../mtmin_OS.pdf",
text="Hi")
sm.addSlide(title="Conclusions", text=content5)
sm.startBackup()
sm.addSlide(title="HT uncorrected",
p1="../htuncorr.pdf",
p2="../htuncorr_OS.pdf",
text="Hi")
sm.addSlide(title="HT corrected", p1="../ht.pdf", p2="../ht_OS.pdf", text="Hi")
sm.addSlide(title="nJets uncorrected",
p1="../njets_uncorr.pdf",
p2="../njets_uncorr_OS.pdf",
text="Hi")
sm.addSlide(title="nJets corrected",
p1="../njets.pdf",
p2="../njets_OS.pdf",
text="Hi")
sm.addSlide(title="nBjets uncorrected",
p1="../nbtags_uncorr.pdf",
p2="../nbtags_uncorr_OS.pdf",
text="Hi")
sm.addSlide(title="nBjets corrected",
p1="../nbtags.pdf",
p2="../nbtags_OS.pdf",
text="Hi")
sm.writeSlides("Phys14_vs_Run2-25ns_v2.tex", opts="--compile --copy")
