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 41 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 193 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 51 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 62 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 229 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 120 of file compare.py.

References LegendSettings(), 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 45 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 15 of file compare.py.

compare.argvs

Definition at line 14 of file compare.py.

compare.den_sel

Definition at line 326 of file compare.py.

compare.hist

Definition at line 373 of file compare.py.

compare.hists

Definition at line 316 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 348 of file compare.py.

compare.loose_id

Definition at line 305 of file compare.py.

compare.num_sel

Definition at line 325 of file compare.py.

compare.reco_cut

Definition at line 304 of file compare.py.

compare.runtype

Definition at line 21 of file compare.py.

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

compare.sampledict

Definition at line 307 of file compare.py.

compare.tfile

Definition at line 322 of file compare.py.

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

compare.tlabel

Definition at line 25 of file compare.py.

Referenced by cscdqm::EventProcessor.getCSCFromMap().

compare.tree

Definition at line 323 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 273 of file compare.py.

compare.xlabel

Definition at line 26 of file compare.py.