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 44 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 196 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 54 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 65 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 232 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 123 of file compare.py.

References LegendSettings(), edm.print(), FastTimerService_cff.range, 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 48 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 18 of file compare.py.

compare.argvs

Definition at line 17 of file compare.py.

compare.den_sel

Definition at line 329 of file compare.py.

compare.hist

Definition at line 376 of file compare.py.

Referenced by DDSpecificsHasNamedValueFilter.accept(), DDSpecificsMatchesValueFilter.accept(), DDSpecificsFilter.accept_impl(), PVValHelper.add(), DTT0Calibration.analyze(), APVCyclePhaseMonitor.beginRun(), Book.book(), HLTDQMFilterEffHists< ObjType >.book1D(), HLTDQMFilterEffHists< ObjType >.book2D(), L1TStage2EMTF.bookHistograms(), TrackerG4SimHitNumberingScheme.buildAll(), PlotMillePede.ClearAdditionalSel(), ClusMultInvestPlots(), Combined2DHisto(), CombinedHisto(), examples::TrackAnalysisAlgorithm.draw(), fwlite::Scanner< Collection >.draw(), fwlite::Scanner< Collection >.draw2D(), fwlite::Scanner< Collection >.drawProf(), L1TEMTFEventInfoClient.dumpContentMonitorElements(), L1TEventInfoClient.dumpContentMonitorElements(), EcalPedHists.endJob(), EcalURecHitHists.endJob(), EcalCosmicsHists.endJob(), DQMHistNormalizer.endRun(), TrigObjTnPHistColl::HistColl.fill(), HLTDQMFilterEffHists< ObjType >.fillHists(), edm::Principal.fillPrincipal(), Rivet::HiggsTemplateCrossSections.finalize(), EgHLTOfflineSummaryClient.getEgHLTSumHist_(), fwlite::Scanner< Collection >.getSameH1(), fwlite::Scanner< Collection >.getSameH2(), fwlite::Scanner< Collection >.getSameProf(), SiStripSummaryCreator.getSummaryME(), HGCalTBMB.HGCalTBMB(), FWEveViewManager.highlightAdded(), HistoAnalyzer< C >.HistoAnalyzer(), MuonResidualsFitter.histogramChi2GaussianFit(), CTPPSCommonDQMSource::GlobalPlots.Init(), EcalDisplaysByEvent.init2DEcalHist(), EcalDisplaysByEvent.init3DEcalHist(), EcalURecHitHists.initHists(), EcalCosmicsHists.initHists(), TopDiLeptonOffline::MonitorEnsemble.loggerBinLabels(), main(), LA_Filler_Fitter.make_and_fit_symmchi2(), SiStripSpyDisplayModule.MakeDigiHist_(), DTGeometryValidate.makeHistogram(), RPCGeometryValidate.makeHistogram(), ValidateGeometry.makeHistogram(), SiStripSpyDisplayModule.MakeProcessedRawDigiHist_(), SiStripSpyDisplayModule.MakeRawDigiHist_(), DDExpandedView.mergedSpecificsV(), SymmetryFit.name(), DQMGenericClient.normalizeToEntries(), DQMRivetClient.normalizeToIntegral(), DQMRivetClient.normalizeToLumi(), helper::ScannerBase.numberOfExtraCuts(), CSCAlignmentCorrections.plot(), MuonResidualsFitter.plotsimple(), MuonResidualsFitter.plotweighted(), MuonGEMRecHitsHarvestor.ProcessBooking(), MuonGEMHitsHarvestor.ProcessBooking(), MuonGEMDigisHarvestor.ProcessBooking(), EcalPedHists.readEBdigis(), EcalPedHists.readEEdigis(), edm::OneLumiPoolSource.readLuminosityBlockAuxiliary_(), DQMRivetClient.scaleByFactor(), TopDiLeptonOffline::MonitorEnsemble.selectionPath(), DDExpandedView.specificsV(), tfwliteselectortest::ThingsTSelector.terminate(), tfwliteselectortest::ThingsTSelector2.terminate(), __class__< T >.terminate(), sistrip::EnsembleCalibrationLA.write_ensembles_plots(), sistrip::MeasureLA.write_report_plots(), and sistrip::EnsembleCalibrationLA.write_samples_plots().

compare.hists

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

compare.loose_id

Definition at line 308 of file compare.py.

compare.num_sel

Definition at line 328 of file compare.py.

compare.reco_cut

Definition at line 307 of file compare.py.

compare.runtype

Definition at line 24 of file compare.py.

compare.sampledict

Definition at line 310 of file compare.py.

compare.tfile

Definition at line 325 of file compare.py.

Referenced by HitParentTest.beginJob(), HOSimHitStudy.beginJob(), EcalSimHitStudy.beginJob(), MaterialBudgetCastorHistos.book(), MaterialBudgetHcalHistos.book(), MaterialBudget.book(), MaterialBudgetForward.book(), CaloSimHitAnalysis.CaloSimHitAnalysis(), CaloSimHitStudy.CaloSimHitStudy(), 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 28 of file compare.py.

Referenced by cscdqm::EventProcessor.getCSCFromMap().

compare.tree

Definition at line 326 of file compare.py.

compare.vardict

Definition at line 276 of file compare.py.

compare.xlabel

Definition at line 29 of file compare.py.