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

Referenced by save().

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

hoverlay()

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

Definition at line 195 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)

LegendSettings()

def compare.LegendSettings	(		leg,
			ncolumn
	)

Definition at line 53 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)

makeCompareVars()

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

Definition at line 64 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)




    

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 231 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)

overlay()

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

Definition at line 122 of file compare.py.

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

Referenced by FWCaloTowerDetailView.build(), FWMuonDetailView.build(), FWPhotonDetailView.build(), FWElectronDetailView.build(), SiStripHitResolFromCalibTree.makeSummary(), SiStripHitEfficiencyHarvester.makeSummary(), 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)

save()

def compare.save	(		canvas,
			name
	)

Definition at line 47 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

argc

compare.argc

Definition at line 17 of file compare.py.

argvs

compare.argvs

Definition at line 16 of file compare.py.

den_sel

compare.den_sel

Definition at line 328 of file compare.py.

hist

compare.hist

Definition at line 375 of file compare.py.

hists

compare.hists

Definition at line 318 of file compare.py.

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

hvardict

compare.hvardict

Definition at line 350 of file compare.py.

loose_id

compare.loose_id

Definition at line 307 of file compare.py.

num_sel

compare.num_sel

Definition at line 327 of file compare.py.

reco_cut

compare.reco_cut

Definition at line 306 of file compare.py.

runtype

compare.runtype

Definition at line 23 of file compare.py.

sampledict

compare.sampledict

Definition at line 309 of file compare.py.

tfile

compare.tfile

Definition at line 324 of file compare.py.

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

tlabel

compare.tlabel

Definition at line 27 of file compare.py.

Referenced by cscdqm::EventProcessor.getCSCFromMap().

tree

compare.tree

Definition at line 325 of file compare.py.

Referenced by SiPixelFedCablingMap.cablingTree(), fwlite::DataGetterHelper.DataGetterHelper(), TrackerOfflineValidationSummary.endJob(), edm::RootOutputFile.finishEndFile(), emtf::Forest.Forest(), hcalCalib.Init(), lutNN::LutNetworkFixedPointRegression2Outputs< input_I, input_F, inputSize, layer1_lut_I, layer1_lut_F, layer1_neurons, layer1_output_I, layer2_input_I, layer2_lut_I, layer2_lut_F, layer2_neurons, layer3_input_I, layer3_0_inputCnt, layer3_0_lut_I, layer3_0_lut_F, output0_I, output0_F, layer3_1_inputCnt, layer3_1_lut_I, layer3_1_lut_F, output1_I, output1_F >.load(), OMTFProcessor< GoldenPatternType >.loadExtrapolFactors(), emtf::Forest.operator=(), lutNN::LutNetworkFixedPointRegression2Outputs< input_I, input_F, inputSize, layer1_lut_I, layer1_lut_F, layer1_neurons, layer1_output_I, layer2_input_I, layer2_lut_I, layer2_lut_F, layer2_neurons, layer3_input_I, layer3_0_inputCnt, layer3_0_lut_I, layer3_0_lut_F, output0_I, output0_F, layer3_1_inputCnt, layer3_1_lut_I, layer3_1_lut_F, output1_I, output1_F >.save(), OMTFProcessor< GoldenPatternType >.saveExtrapolFactors(), PVFitter.setTree(), LagrangeChildUpdator.update(), and SummaryTableOutputBranches.updateBranches().

vardict

compare.vardict

Definition at line 275 of file compare.py.

xlabel

compare.xlabel

Definition at line 28 of file compare.py.