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

ensureDir()

def compare.ensureDir

(

directory

)

Definition at line 44 of file compare.py.

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

Referenced by save().

hoverlay()

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

Definition at line 196 of file compare.py.

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)
 
 
 
 

References LegendSettings(), and save().

LegendSettings()

def compare.LegendSettings	(		leg,
			ncolumn
	)

Definition at line 54 of file compare.py.

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)
 
 
 

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

makeCompareVars()

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

Definition at line 65 of file compare.py.

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)
 
 
 
 
    

References LegendSettings(), and str.

makeEffPlotsVars()

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 232 of file compare.py.

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)
 
 
 
 
 

References mps_check.array.

overlay()

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

Definition at line 123 of file compare.py.

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)
 
 
 

References LegendSettings(), print(), FastTimerService_cff.range, and save().

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

save()

def compare.save	(		canvas,
			name
	)

Definition at line 48 of file compare.py.

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

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

Referenced by hoverlay(), and overlay().

Variable Documentation

argc

compare.argc

Definition at line 18 of file compare.py.

argvs

compare.argvs

Definition at line 17 of file compare.py.

den_sel

compare.den_sel

Definition at line 329 of file compare.py.

hist

compare.hist

Definition at line 376 of file compare.py.

hists

compare.hists

Definition at line 319 of file compare.py.

Referenced by trklet::TrackletCalculatorBase.barrelSeeding(), AlignmentMonitorGeneric.book(), TrigObjTnPHistColl::HistDefs.bookHists(), ColourStatsBoxes(), CosmicGenFilterHelix.createHistsEnd(), CosmicGenFilterHelix.createHistsStart(), EcalSelectiveReadoutValidation.EcalSelectiveReadoutValidation(), trklet::MatchCalculator.execute(), PatZToMuMuAnalyzer.fill(), MillePedeMonitor.fillResidualHists(), MillePedeMonitor.fillResidualHitHists(), getHists(), PedsFullNoiseHistograms.histoAnalysis(), main(), trklet::MatchProcessor.matchCalculator(), CosmicGenFilterHelix.monitorEnd(), CosmicGenFilterHelix.monitorStart(), SetMinMaxRange(), and setStack().

hvardict

compare.hvardict

Definition at line 351 of file compare.py.

loose_id

compare.loose_id

Definition at line 308 of file compare.py.

num_sel

compare.num_sel

Definition at line 328 of file compare.py.

reco_cut

compare.reco_cut

Definition at line 307 of file compare.py.

runtype

compare.runtype

Definition at line 24 of file compare.py.

sampledict

compare.sampledict

Definition at line 310 of file compare.py.

tfile

compare.tfile

Definition at line 325 of file compare.py.

Referenced by HitParentTest.beginJob(), HcalSimHitAnalysis.beginJob(), HOSimHitStudy.beginJob(), EcalSimHitStudy.beginJob(), MaterialBudgetCastorHistos.book(), MaterialBudgetHcalHistos.book(), MaterialBudget.book(), MaterialBudgetForward.book(), CaloSimHitAnalysis.CaloSimHitAnalysis(), CaloSimHitStudy.CaloSimHitStudy(), CherenkovAnalysis.CherenkovAnalysis(), ECalSD.ECalSD(), ElectronStudy.ElectronStudy(), HCalSD.HCalSD(), HcalTB02Histo.HcalTB02Histo(), HcalTB04Histo.HcalTB04Histo(), HcalTB06Histo.HcalTB06Histo(), HFGflash.HFGflash(), HFShowerParam.HFShowerParam(), HGCalTBMB.HGCalTBMB(), HLTCSCOverlapFilter.HLTCSCOverlapFilter(), HGCSD.initRun(), and XtalDedxAnalysis.XtalDedxAnalysis().

tlabel

compare.tlabel

Definition at line 28 of file compare.py.

Referenced by cscdqm::EventProcessor.getCSCFromMap().

tree

compare.tree

Definition at line 326 of file compare.py.

vardict

compare.vardict

Definition at line 276 of file compare.py.

xlabel

compare.xlabel

Definition at line 29 of file compare.py.