Exemplo n.º 1
0
        #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')
Exemplo n.º 2
0
 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)")
Exemplo n.º 3
0
    #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