#plotter.save('mcdata-os-p1p2f3-nvtx')
# Lower stat but closer to signal region
#plotter.plot_mc_vs_data('os/p1p2p3', 'rho')
#plotter.save('mcdata-os-p1p2p3-rho')
#plotter.plot_mc_vs_data('os/p1p2p3', 'nvtx')
#plotter.save('mcdata-os-p1p2p3-nvtx')
# Make Z->mumu + tau jet control
weighted = plotter.plot('data',
'os/p1p2f3/w3/e_m_Mass',
'hist',
rebin=20,
styler=make_styler(2, 'hist'),
xaxis='m_{e#mu} (GeV)')
unweighted = plotter.plot('data',
'os/p1p2p3/e_m_Mass',
'same',
rebin=20,
styler=make_styler(1),
xaxis='m_{e#mu} (GeV)')
weighted.SetTitle('e^{+}#mu^{-} + fake #tau_{h} est.')
weighted.legendstyle = 'l'
unweighted.SetTitle('e^{+}#mu^{-} + fake #tau_{h} obs.')
unweighted.legendstyle = 'pe'
plotter.add_legend([weighted, unweighted])
plotter.add_cms_blurb(sqrts)
plotter.save('ztt-os-fr-control')
plotter.plot_mc_vs_data('os/p1p2f3', 'rho')
plotter.save('mcdata-os-p1p2f3-rho')
plotter.plot_mc_vs_data('os/p1p2f3', 'nvtx')
plotter.save('mcdata-os-p1p2f3-nvtx')
# Lower stat but closer to signal region
plotter.plot_mc_vs_data('os/p1p2p3', 'rho')
plotter.save('mcdata-os-p1p2p3-rho')
plotter.plot_mc_vs_data('os/p1p2p3', 'nvtx')
plotter.save('mcdata-os-p1p2p3-nvtx')
# Make Z->mumu + tau jet control
antiiso_m2JetPt = plotter.plot('data', 'ss/p1f2p3/m2JetPt', 'hist', styler=make_styler(2, 'hist'), xrange=(0, 120), rebin=10)
antiiso_m2JetPt.SetTitle("Anti-iso CR yield")
antiiso_m2JetPt.legendstyle='l'
antiiso_m2JetPt.GetXaxis().SetTitle("#mu_{2} Jet Pt")
plotter.save('data-p1f2p3-m2JetPt')
antiiso_m1JetPt = plotter.plot('data', 'ss/f1p2p3/m1JetPt', 'hist', styler=make_styler(2, 'hist'), xrange=(0, 120), rebin=10)
antiiso_m1JetPt.SetTitle("Anti-iso CR yield")
antiiso_m1JetPt.legendstyle='l'
antiiso_m1JetPt.GetXaxis().SetTitle("#mu_{1} Jet Pt")
plotter.save('data-f1p2p3-m1JetPt')
zmm_weighted = plotter.plot('data', 'os/p1p2f3/w3/m1_m2_Mass', 'hist', styler=make_styler(2, 'hist'), xrange=(60, 120))
zmm_weighted.SetTitle("Z#mu#mu + fake #tau_{h} est.")
zmm_weighted.legendstyle='l'
zmm_weighted.GetXaxis().SetTitle("m_{#mu#mu} (GeV)")
#plotter.plot_mc_vs_data('os/p1p2f3', 'rho')
#plotter.save('mcdata-os-p1p2f3-rho')
#plotter.plot_mc_vs_data('os/p1p2f3', 'nvtx')
#plotter.save('mcdata-os-p1p2f3-nvtx')
# Lower stat but closer to signal region
#plotter.plot_mc_vs_data('os/p1p2p3', 'rho')
#plotter.save('mcdata-os-p1p2p3-rho')
#plotter.plot_mc_vs_data('os/p1p2p3', 'nvtx')
#plotter.save('mcdata-os-p1p2p3-nvtx')
# Make Z->mumu + tau jet control
weighted = plotter.plot('data', 'os/p1p2f3/w3/e_m_Mass', 'hist', rebin=20, styler=make_styler(2, 'hist'), xaxis='m_{e#mu} (GeV)')
unweighted = plotter.plot('data', 'os/p1p2p3/e_m_Mass', 'same', rebin=20, styler=make_styler(1), xaxis='m_{e#mu} (GeV)')
weighted.SetTitle('e^{+}#mu^{-} + fake #tau_{h} est.')
weighted.legendstyle = 'l'
unweighted.SetTitle('e^{+}#mu^{-} + fake #tau_{h} obs.')
unweighted.legendstyle = 'pe'
plotter.add_legend([weighted, unweighted])
plotter.add_cms_blurb(sqrts)
plotter.save('ztt-os-fr-control')
#plotter.plot('data', 'os/p1p2p3/prescale', styler=make_styler(1))
#plotter.save('ztt-os-prescale-check')
#plotter.plot('Zjets_M50', 'os/p1p2f3/weight')
#plotter.save('ztt-mc-event-weights')
## Check MC weights
