compare Namespace Reference

Functions
def	ensureDir (directory)
def	hoverlay (hists, xtitle, ytitle, name)
def	LegendSettings (leg, ncolumn)
def	makeCompareVars (tree, var, sel, leglist, nbin, xmin, xmax, xtitle, ytitle, scale, name)
def	makeEffPlotsVars (tree, varx, vary, sel, nbinx, xmin, xmax, nbiny, ymin, ymax, xtitle, ytitle, leglabel=None, header='', addon='', option='pt', marker=20)
def	overlay (hists, ytitle, header, addon)
def	save (canvas, name)

Variables
	argc
	argvs
	den_sel
	hist
	hists
	hvardict
	loose_id
	num_sel
	reco_cut
	runtype
	sampledict
	tfile
	tlabel
	tree
	vardict
	xlabel

Function Documentation

def compare.ensureDir

(

directory

)

Definition at line 42 of file compare.py.

Referenced by save().

def ensureDir(directory):
    if not os.path.exists(directory):
        os.makedirs(directory)

def compare.hoverlay	(		hists,
			xtitle,
			ytitle,
			name
	)

Definition at line 194 of file compare.py.

References LegendSettings(), and save().

def hoverlay(hists, xtitle, ytitle, name):
   
    c = TCanvas()

    ymax = 0
    for hist in hists:
        if ymax < hist.GetMaximum():
            ymax = hist.GetMaximum()


    leg = TLegend(0.6,0.65,0.91,0.9)
    LegendSettings(leg,1)

    for ii, ihist in enumerate(hists):
        ihist.SetMaximum(ymax*1.2)
        ihist.SetMinimum(0.)
        ihist.SetMarkerSize(0.)
        ihist.GetXaxis().SetTitle(xtitle)
        ihist.GetYaxis().SetTitle(ytitle)

        if ii ==0:
            ihist.Draw('h')
        else:
            ihist.Draw('hsame')

        leg.AddEntry(ihist, ihist.GetName(), "l")


    leg.Draw()


    save(c, 'compare_' + runtype + '/hist_' + name)

def compare.LegendSettings	(		leg,
			ncolumn
	)

Definition at line 52 of file compare.py.

Referenced by hoverlay(), makeCompareVars(), and overlay().

def LegendSettings(leg, ncolumn):
    leg.SetNColumns(ncolumn)
    leg.SetBorderSize(0)
    leg.SetFillColor(10)
    leg.SetLineColor(0)
    leg.SetFillStyle(0)
    leg.SetTextSize(0.03)
    leg.SetTextFont(42)

def compare.makeCompareVars	(		tree,
			var,
			sel,
			leglist,
			nbin,
			xmin,
			xmax,
			xtitle,
			ytitle,
			scale,
			name
	)

Definition at line 63 of file compare.py.

References LegendSettings(), and str.

def makeCompareVars(tree, var, sel, leglist, nbin, xmin, xmax, xtitle, ytitle, scale, name):
   
#    print leglist

    c = TCanvas()

    hists = []
    col = [1,2,4,8,6]
    ymax = 0

    for ii, isel in enumerate(sel):
        hist = TH1F('h_' + str(ii), 'h_' + str(ii), nbin, xmin, xmax)
        hist.GetXaxis().SetTitle(xtitle)
        hist.GetYaxis().SetTitle(ytitle)
        hist.GetYaxis().SetNdivisions(507)
        hist.SetLineColor(col[ii])
        hist.SetLineWidth(len(sel)+1-ii)
        hist.SetLineStyle(1)
        hist.SetMarkerSize(0)
        hist.SetMinimum(0)
        hist.Sumw2()

#        print hist.GetName(), var, isel
        tree.Project(hist.GetName(), var, isel)
        hist.Scale(1./hist.GetEntries())

        if ymax < hist.GetMaximum():
            ymax = hist.GetMaximum()

        hists.append(hist)


    leg = TLegend(0.6,0.65,0.91,0.9)
    LegendSettings(leg,1)

    for ii, ihist in enumerate(hists):
        ihist.SetMaximum(ymax*1.2)
        ihist.SetMinimum(0.)

        if ii ==0:
            ihist.Draw('h')
        else:
            ihist.Draw('hsame')

        if leglist[ii] != 'None':
            leg.AddEntry(ihist, leglist[ii], "l")


    if leglist[0]!='None':
        leg.Draw()


#    save(c, 'compare_' + runtype + '/compare_' + name)




    

def compare.makeEffPlotsVars	(		tree,
			varx,
			vary,
			sel,
			nbinx,
			xmin,
			xmax,
			nbiny,
			ymin,
			ymax,
			xtitle,
			ytitle,
			leglabel = None,
			header = '',
			addon = '',
			option = 'pt',
			marker = 20
	)

Definition at line 230 of file compare.py.

References mps_check.array.

def makeEffPlotsVars(tree, varx, vary, sel, nbinx, xmin, xmax, nbiny, ymin, ymax, xtitle, ytitle, leglabel = None, header='', addon='', option='pt', marker=20):
   
    binning = [20,30,40,50,60,70,80,100,150,200]

    c = TCanvas()

    if option=='pt':
        _hist_ = TH1F('h_effp_' + addon, 'h_effp' + addon, len(binning)-1, array('d',binning))
        _ahist_ = TH1F('ah_effp_' + addon, 'ah_effp' + addon, len(binning)-1, array('d',binning))
    elif option=='eta':
        _hist_ = TH1F('h_effp_' + addon, 'h_effp' + addon, nbinx, xmin, xmax)
        _ahist_ = TH1F('ah_effp_' + addon, 'ah_effp' + addon, nbinx, xmin, xmax)
    elif option=='nvtx':
        _hist_ = TH1F('h_effp_' + addon, 'h_effp' + addon, len(vbinning)-1, array('d',vbinning))
        _ahist_ = TH1F('ah_effp_' + addon, 'ah_effp' + addon, len(vbinning)-1, array('d',vbinning))


    tree.Draw(varx + ' >> ' + _hist_.GetName(), sel)
    tree.Draw(varx + ' >> ' + _ahist_.GetName(), sel + ' && ' + vary)
    
    g_efficiency = TGraphAsymmErrors()
    g_efficiency.BayesDivide(_ahist_, _hist_)
    g_efficiency.GetXaxis().SetTitle(xtitle)
    g_efficiency.GetYaxis().SetTitle('efficiency')
    g_efficiency.GetYaxis().SetNdivisions(507)
    g_efficiency.SetLineWidth(3)
    g_efficiency.SetName(header)    
    g_efficiency.SetMinimum(0.)
    g_efficiency.GetYaxis().SetTitleOffset(1.3)
    g_efficiency.SetMarkerStyle(marker)
    g_efficiency.SetMarkerSize(1)
    g_efficiency.Draw('ap')

#    save(c, 'plots/' + addon)

    return copy.deepcopy(g_efficiency)

def compare.overlay	(		hists,
			ytitle,
			header,
			addon
	)

Definition at line 121 of file compare.py.

References LegendSettings(), edm.print(), and save().

Referenced by FWCaloTowerDetailView.build(), FWMuonDetailView.build(), FWPhotonDetailView.build(), FWElectronDetailView.build(), and SiStripHitEffFromCalibTree.makeSummary().

def overlay(hists, ytitle, header, addon):

    print('number of histograms = ', len(hists))

    canvas = TCanvas()
    leg = TLegend(0.2,0.7,0.5,0.9)
    LegendSettings(leg, 1)

    col = [1,2,4,6,8,9,12]

    ymax = -1
    ymin = 100

    for ii, hist in enumerate(hists):
        hist.GetYaxis().SetTitle('efficiency')
        hist.SetLineColor(col[ii])
        hist.SetMarkerColor(col[ii])
        hist.SetLineWidth(2)
        hist.SetMarkerSize(1)


        for ip in range(hist.GetN()):
            x = Double(-1)
            y = Double(-1)
            hist.GetPoint(ip, x, y)

            if ymin > y:
                ymin = y
            if ymax < y:
                ymax = y

#
#        if ymax < hist.GetMaximum():
#            ymax = hist.GetMaximum()
#        if ymin > hist.GetMinimum():
#            ymin = hist.GetMinimum()

        if ii==0:
            hist.Draw("Azp")
        else:
            
            hist.Draw("pzsame")
 
#        print hist.GetName(), hist.GetTitle()
        legname = hist.GetName()

        leg.AddEntry(hist, legname, 'lep')


    for hist in hists:
        hist.SetMaximum(ymax*2)
#        hist.SetMinimum(ymin*0.5)

    leg.Draw()

    tex = TLatex( hists[-1].GetXaxis().GetXmin() + 0.01*(hists[-1].GetXaxis().GetXmax() - hists[-1].GetXaxis().GetXmin()), ymax*2.1, addon.replace('tau_',''))

    tex.SetTextFont(42)
    tex.SetTextSize(0.03)
    tex.Draw()

    tex2 = TLatex( hists[-1].GetXaxis().GetXmin() + 0.87*(hists[-1].GetXaxis().GetXmax() - hists[-1].GetXaxis().GetXmin()), ymax*2.1, tlabel)

    tex2.SetTextFont(42)
    tex2.SetTextSize(0.03)
    tex2.Draw()



    save(canvas, 'compare_' + runtype + '/' + header)

def compare.save	(		canvas,
			name
	)

Definition at line 46 of file compare.py.

References ensureDir(), and python.rootplot.root2matplotlib.replace().

Referenced by hoverlay(), and overlay().

def save(canvas, name):
    ensureDir('compare_' + runtype)
    canvas.SaveAs(name.replace(' ','').replace('&&','')+'.pdf')
    canvas.SaveAs(name.replace(' ','').replace('&&','')+'.gif')

Variable Documentation

compare.argc

Definition at line 16 of file compare.py.

compare.argvs

Definition at line 15 of file compare.py.

compare.den_sel

Definition at line 327 of file compare.py.

compare.hist

Definition at line 374 of file compare.py.

compare.hists

Definition at line 317 of file compare.py.

Referenced by AlignmentMonitorGeneric.book(), TrigObjTnPHistColl::HistDefs.bookHists(), EcalSelectiveReadoutValidation.EcalSelectiveReadoutValidation(), PatZToMuMuAnalyzer.fill(), MillePedeMonitor.fillResidualHitHists(), GFOverlay.GetHistManager(), PedsFullNoiseHistograms.histoAnalysis(), and main().

compare.hvardict

Definition at line 349 of file compare.py.

compare.loose_id

Definition at line 306 of file compare.py.

compare.num_sel

Definition at line 326 of file compare.py.

compare.reco_cut

Definition at line 305 of file compare.py.

compare.runtype

Definition at line 22 of file compare.py.

Referenced by RPCTechTriggerConfig.RPCTechTriggerConfig(), and RunTypeDef.setRunType().

compare.sampledict

Definition at line 308 of file compare.py.

compare.tfile

Definition at line 323 of file compare.py.

Referenced by MaterialBudgetCastorHistos.book(), MaterialBudgetHcalHistos.book(), MaterialBudget.book(), MaterialBudgetForward.book(), DreamSD.DreamSD(), ECalSD.ECalSD(), ElectronStudy.ElectronStudy(), HCalSD.HCalSD(), HcalTB02Histo.HcalTB02Histo(), HcalTB04Histo.HcalTB04Histo(), HcalTB06Histo.HcalTB06Histo(), HFGflash.HFGflash(), HFShowerParam.HFShowerParam(), HGCalTBMB.HGCalTBMB(), HLTCSCOverlapFilter.HLTCSCOverlapFilter(), and HGCSD.initRun().

compare.tlabel

Definition at line 26 of file compare.py.

Referenced by cscdqm::EventProcessor.getCSCFromMap().

compare.tree

Definition at line 324 of file compare.py.

Referenced by SiPixelFedCablingMap.cablingTree(), pftools::Exercises3.calibrateCalibratables(), pftools::CalibCompare.calibrateCalibratables(), fwlite::DataGetterHelper.DataGetterHelper(), SummaryTableOutputBranches.defineBranchesFromFirstEvent(), TrackerOfflineValidationSummary.endJob(), edm::RootOutputFile.finishEndFile(), emtf::Forest.Forest(), emtf::Forest.operator=(), PVFitter.setTree(), and LagrangeChildUpdator.update().

compare.vardict

Definition at line 274 of file compare.py.

compare.xlabel

Definition at line 27 of file compare.py.