Python Plotter.plot_mc_vs_dataの例

コード例 #1

ファイル: PlotControlZMM.py プロジェクト: sshalhou/FinalStateAnalysis

    'WW*',
    'TT*',
    'WplusJets*',
    "data_DoubleMu*",
]

files = []
lumifiles = []

for x in samples:
    files.extend(glob.glob('results/%s/ControlZMM/%s.root' % (jobid, x)))
    lumifiles.extend(glob.glob('inputs/%s/%s.lumicalc.sum' % (jobid, x)))

plotter = Plotter(files, lumifiles, output_dir)

plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=4)
plotter.save('mass')

plotter.plot_mc_vs_data('zmm', 'm1Pt')
plotter.save('m1Pt')
plotter.plot_mc_vs_data('zmm', 'm2Pt')
plotter.save('m2Pt')

plotter.plot_mc_vs_data('zmm', 'm1AbsEta')
plotter.save('m1AbsEta')
plotter.plot_mc_vs_data('zmm', 'm2AbsEta')
plotter.save('m2AbsEta')

plotter.plot_mc_vs_data('zmm', 'nvtx')
plotter.save('nvtx')

コード例 #2

ファイル: plotTauFakeRateMVA.py プロジェクト: fanbomeng/UWHiggs

print plotter.mc_samples

for i in tauiso:
    for s in sign:
        #for histo, axis in zip(histoname, axistitle):
        #        for n,h in enumerate(histoname) :
        foldername = s + '/' + i
        if not os.path.exists(outputdir + foldername):
            os.makedirs(outputdir + foldername)

        if i == 'tNoCuts':
            plotter.plot_mc_vs_data(foldername,
                                    'CUT_FLOW',
                                    rebin=1,
                                    xaxis='CUT_FLOW',
                                    leftside=False,
                                    show_ratio=False,
                                    ratio_range=1.5,
                                    sorted=True)
            plotter.save(foldername + '/CUT_FLOW')

        foldername = s + '/' + i
        for h in histoname:
            #plotter.canvas.SetLogy(True)
            #plotter.plot_data(foldername, h, rebin=rebins[n],  xaxis= axistitle[n], leftside=False)
            #plotter.plot_mc_vs_data(foldername,h, rebin=rebins[n], xaxis= axistitle[n], leftside=False, show_ratio=True, ratio_range=1.5, sort=True)

            plotter.plot_mc_vs_data(foldername,
                                    h[0],
                                    rebin=h[2],
                                    xaxis=h[1],

コード例 #3

    "data_DoubleMu*",
]

files = []
lumifiles = []

for x in samples:
    files.extend(glob.glob('results/%s/ControlZMM/%s.root' % (jobid, x)))
    lumifiles.extend(glob.glob('inputs/%s/%s.lumicalc.sum' % (jobid, x)))

plotter = Plotter(files, lumifiles, output_dir)
plotter.mc_samples = ['Zjets_M50']

sqrts = 7 if '7TeV' in jobid else 8

plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=2, xaxis='m_{#mu#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')

plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=6, xaxis='m_{#mu#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass_rebin')

plotter.plot_mc_vs_data('zmm', 'm1Pt')
plotter.save('m1Pt')
plotter.plot_mc_vs_data('zmm', 'm1Pt', 5)
plotter.save('m1Pt_rebin')
plotter.plot_mc_vs_data('zmm', 'm2Pt')
plotter.save('m2Pt')
plotter.plot_mc_vs_data('zmm', 'm2Pt', 5)
plotter.save('m2Pt_rebin')

コード例 #4

ファイル: PlotControlEM.py プロジェクト: lgray/UWHiggs

plotter.views['QCD'] = { 'view': qcd_view }

# Override ordering
plotter.mc_samples = [
    #'WZJetsTo3LNu*',
    #'ZZJetsTo4L*',
    'QCD',
    'WW*',
    'WplusJets_madgraph',
    'TTplusJets_madgraph',
    'Zjets_M50',
]


plotter.plot_mc_vs_data('em', 'emMass', rebin=5, leftside=False,
                        xaxis='m_{e#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')

plotter.plot_mc_vs_data('em', 'mPt')
plotter.save('mPt')
plotter.plot_mc_vs_data('em', 'ePt', rebin=10)
plotter.save('ePt')

plotter.plot_mc_vs_data('em', 'mAbsEta')
plotter.save('mAbsEta')
plotter.plot_mc_vs_data('em', 'eAbsEta', rebin=5)
plotter.save('eAbsEta')

plotter.plot_mc_vs_data('em', 'nvtx')
plotter.save('nvtx')

コード例 #5

ファイル: PlotEM.py プロジェクト: mverzett/UWHiggs

plotter = Plotter(files, lumifiles, output_dir)

# Override ordering
plotter.mc_samples = [
    'TTplusJets_madgraph',
    'WplusJets_madgraph',
    'Zjets_M50',
    'WZJetsTo3LNu*',
    'WW*',
    'ZZJetsTo4L*',
]

sqrts = 7 if '7TeV' in jobid else 8

plotter.plot_mc_vs_data('em', 'emMass', rebin=10, leftside=False,
                        xaxis='m_{e#mu} (GeV)')
plotter.add_cms_blurb(sqrts)
plotter.save('mass')

plotter.plot_mc_vs_data('em', 'mPt')
plotter.save('mPt')
plotter.plot_mc_vs_data('em', 'ePt')
plotter.save('ePt')

plotter.plot_mc_vs_data('em', 'mAbsEta')
plotter.save('mAbsEta')
plotter.plot_mc_vs_data('em', 'eAbsEta')
plotter.save('eAbsEta')

plotter.plot_mc_vs_data('em', 'nvtx')
plotter.save('nvtx')

コード例 #6

    'm1 p_{T} (GeV)', 'm2 p_{T} (GeV)', 'm1-m2 Inv Mass (GeV)',
    'm1 Gen p_{T} (GeV)', 'm2 Gen p_{T} (GeV)'
]

#rebins = [5, 5, 2, 5, 5, 2, 1, 1, 1,  1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 5, 1, 5, 1, 5, 1, 5]
rebins = []
for n in histoname:
    rebins.append(1)

plotter.mc_samples = new_mc_samples
#plotter.mc_samples = mc_samples
for i in sign:
    for n, h in enumerate(histoname):
        foldername = i
        #plotter.pad.SetLogy(True)

        #        plotter.plot_mc(foldername, ['ggHiggsToETau','vbfHiggsToETau'],h, rebin=rebins[n], xaxis= axistitle[n], leftside=False, show_ratio=False, ratio_range=1.5,  rescale=10)
        plotter.plot_mc_vs_data(foldername,
                                h,
                                rebin=rebins[n],
                                xaxis=axistitle[n],
                                leftside=False,
                                show_ratio=True,
                                ratio_range=1.5,
                                sort=True)

        if not os.path.exists(outputdir + foldername):
            os.makedirs(outputdir + foldername)

        plotter.save(foldername + '/' + h)

コード例 #7

rebins = []
for n in histoname:
    rebins.append(1)

for n, h in enumerate(histoname):
    plotter.pad.SetLogy(True)

    #    plotter.plot_mc(foldername, h, 1, xaxis= axistitle[n], leftside=False, show_ratio=False, ratio_range=1.5,  sort=False)
    print "histo", h
    if h == 'm1Pt' or h == 'm2Pt':
        plotter.plot_mc_vs_data(foldername,
                                h,
                                1,
                                xaxis=axistitle[n],
                                legend_label=ll,
                                leftside=False,
                                xrange=[0, 200],
                                show_ratio=True,
                                ratio_range=1.5,
                                sort=True)
    elif h == 'm1_m2_Mass':
        plotter.plot_mc_vs_data(foldername,
                                h,
                                1,
                                xaxis=axistitle[n],
                                legend_label=ll,
                                leftside=False,
                                xrange=[20, 200],
                                show_ratio=True,
                                ratio_range=1.5,
                                sort=True)

コード例 #8

ファイル: PlotControlZMM.py プロジェクト: mverzett/UWHiggs

    "data_DoubleMu*",
]

files = []
lumifiles = []

for x in samples:
    files.extend(glob.glob('results/%s/ControlZMM/%s.root' % (jobid, x)))
    lumifiles.extend(glob.glob('inputs/%s/%s.lumicalc.sum' % (jobid, x)))

plotter = Plotter(files, lumifiles, output_dir)
plotter.mc_samples = filter(lambda x: 'data' not in x.lower(), samples) #['Zjets_M50']

sqrts = 7 if '7TeV' in jobid else 8

plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=2, xaxis='m_{#mu#mu} (GeV)')#, show_ratio=True)
plotter.add_cms_blurb(sqrts)
plotter.save('mass')

plotter.plot_mc_vs_data('zmm', 'm1m2Mass', rebin=6, xaxis='m_{#mu#mu} (GeV)')#, show_ratio=True)
plotter.add_cms_blurb(sqrts)
plotter.save('mass_rebin')

plotter.plot_mc_vs_data('zmm', 'm1Pt', show_ratio=True)
plotter.save('m1Pt')
plotter.plot_mc_vs_data('zmm', 'm1Pt', 5, show_ratio=True)
plotter.save('m1Pt_rebin')
plotter.plot_mc_vs_data('zmm', 'm2Pt', show_ratio=True)
plotter.save('m2Pt')
plotter.plot_mc_vs_data('zmm', 'm2Pt', 5, show_ratio=True)
plotter.save('m2Pt_rebin')

コード例 #9

ファイル: PlotControlEM.py プロジェクト: mverzett/UWHiggs

    'Fakes',
    'WWJetsTo2L2Nu-no-fakes',
    #'WW*',
    'TTplusJets_madgraph-no-fakes',
    #'TTplusJets_madgraph',
    #'WZJetsTo3LNu*',
    #'ZZJetsTo4L*',
    'Zjets_M50-no-fakes',
    #'Zjets_M50',
]

for suffix, samples in [('', os_ss_samples), ('-fakes', fakes_samples)]:
    plotter.mc_samples = samples

    #from pdb import set_trace; set_trace()
    plotter.plot_mc_vs_data('em', 'emMass', rebin=5, leftside=False,
                            xaxis='m_{e#mu} (GeV)')#, show_ratio=True)
    plotter.add_cms_blurb(sqrts)
    plotter.save('mass' + suffix)

    plotter.plot_mc_vs_data('em', 'emMass', rebin=10, leftside=False,
                            xaxis='m_{e#mu} (GeV)')#, show_ratio=True)
    plotter.add_cms_blurb(sqrts)
    plotter.save('mass_rebin' + suffix)

    plotter.plot_mc_vs_data('em', 'mPt', rebin=10, show_ratio=True)
    plotter.save('mPt' + suffix)
    plotter.plot_mc_vs_data('em', 'ePt', rebin=10, show_ratio=True)
    plotter.save('ePt' + suffix)

    plotter.plot_mc_vs_data('em', 'mAbsEta', rebin=5, show_ratio=True)
    plotter.save('mAbsEta' + suffix)