Public Member Functions
def	__init__ (self, name, kwargs)

def	addToLegend (self, legend, legendLabels, denomUncertainty)

def	clone (self, kwargs)

def	create (self, tdirNEvents, requireAllHistograms=False)

def	draw (self, pad, ratio, ratioFactor, nrows)

def	drawRatioUncertainty (self)

def	getName (self)

def	getNumberOfHistograms (self)

def	isEmpty (self)

def	isRatio (self, ratio)

def	isTGraph2D (self)

def	setProperties (self, kwargs)

Private Member Functions
def	_createOne (self, name, index, tdir, nevents)

def	_normalize (self)

def	_setStats (self, histos, startingX, startingY)

Private Attributes
	_forLegend

	_frame

	_histograms

	_mainAdditional

	_name

	_ratioAdditional

	_ratios

	_tmp_histos

Detailed Description

Represents one plot, comparing one or more histograms.

Definition at line 1494 of file plotting.py.

Constructor & Destructor Documentation

def plotting.Plot.__init__	(	self,
		name,
		kwargs
	)

Constructor.

Arguments:
name -- String for name of the plot, or Efficiency object

Keyword arguments:
fallback     -- Dictionary for specifying fallback (default None)
outname      -- String for an output name of the plot (default None for the same as 'name')
title        -- String for a title of the plot (default None)
xtitle       -- String for x axis title (default None)
xtitlesize   -- Float for x axis title size (default None)
xtitleoffset -- Float for x axis title offset (default None)
xlabelsize   -- Float for x axis label size (default None)
ytitle       -- String for y axis title (default None)
ytitlesize   -- Float for y axis title size (default None)
ytitleoffset -- Float for y axis title offset (default None)
ztitle       -- String for z axis title (default None)
ztitleoffset -- Float for z axis title offset (default None)
xmin         -- Float for x axis minimum (default None, i.e. automatic)
xmax         -- Float for x axis maximum (default None, i.e. automatic)
ymin         -- Float for y axis minimum (default 0)
ymax         -- Float for y axis maximum (default None, i.e. automatic)
xlog         -- Bool for x axis log status (default False)
ylog         -- Bool for y axis log status (default False)
xgrid        -- Bool for x axis grid status (default True)
ygrid        -- Bool for y axis grid status (default True)
stat         -- Draw stat box? (default False)
fit          -- Do gaussian fit? (default False)
statx        -- Stat box x coordinate (default 0.65)
staty        -- Stat box y coordinate (default 0.8)
statyadjust  -- List of floats for stat box y coordinate adjustments (default None)
normalizeToUnitArea -- Normalize histograms to unit area? (default False)
normalizeToNumberOfEvents -- Normalize histograms to number of events? If yes, the PlotFolder needs 'numberOfEventsHistogram' set to a histogram filled once per event (default False)
profileX     -- Take histograms via ProfileX()? (default False)
fitSlicesY   -- Take histograms via FitSlicesY() (default False)
rebinX       -- rebin x axis (default None)
scale        -- Scale histograms by a number (default None)
xbinlabels   -- List of x axis bin labels (if given, default None)
xbinlabelsize -- Size of x axis bin labels (default None)
xbinlabeloption -- Option string for x axis bin labels (default None)
removeEmptyBins -- Bool for removing empty bins, but only if histogram has bin labels (default False)
printBins    -- Bool for printing bin values, but only if histogram has bin labels (default False)
drawStyle    -- If "hist", draw as line instead of points (default None)
drawCommand  -- Deliver this to Draw() (default: None for same as drawStyle)
lineWidth    -- If drawStyle=="hist", the width of line (default 2)
legendDx     -- Float for moving TLegend in x direction for separate=True (default None)
legendDy     -- Float for moving TLegend in y direction for separate=True (default None)
legendDw     -- Float for changing TLegend width for separate=True (default None)
legendDh     -- Float for changing TLegend height for separate=True (default None)
legend       -- Bool to enable/disable legend (default True)
adjustMarginRight  -- Float for adjusting right margin (default None)
ratio        -- Possibility to disable ratio for this particular plot (default None)
ratioYmin    -- Float for y axis minimum in ratio pad (default: list of values)
ratioYmax    -- Float for y axis maximum in ratio pad (default: list of values)
ratioFit     -- Fit straight line in ratio? (default None)
ratioUncertainty -- Plot uncertainties on ratio? (default True)
ratioCoverageXrange -- Range of x axis values (xmin,xmax) to limit the automatic ratio y axis range calculation to (default None for disabled)
histogramModifier -- Function to be called in create() to modify the histograms (default None)

Definition at line 1496 of file plotting.py.

     def __init__(self, name, **kwargs):
         """ Constructor.
 
         Arguments:
         name -- String for name of the plot, or Efficiency object
 
         Keyword arguments:
         fallback     -- Dictionary for specifying fallback (default None)
         outname      -- String for an output name of the plot (default None for the same as 'name')
         title        -- String for a title of the plot (default None)
         xtitle       -- String for x axis title (default None)
         xtitlesize   -- Float for x axis title size (default None)
         xtitleoffset -- Float for x axis title offset (default None)
         xlabelsize   -- Float for x axis label size (default None)
         ytitle       -- String for y axis title (default None)
         ytitlesize   -- Float for y axis title size (default None)
         ytitleoffset -- Float for y axis title offset (default None)
         ztitle       -- String for z axis title (default None)
         ztitleoffset -- Float for z axis title offset (default None)
         xmin         -- Float for x axis minimum (default None, i.e. automatic)
         xmax         -- Float for x axis maximum (default None, i.e. automatic)
         ymin         -- Float for y axis minimum (default 0)
         ymax         -- Float for y axis maximum (default None, i.e. automatic)
         xlog         -- Bool for x axis log status (default False)
         ylog         -- Bool for y axis log status (default False)
         xgrid        -- Bool for x axis grid status (default True)
         ygrid        -- Bool for y axis grid status (default True)
         stat         -- Draw stat box? (default False)
         fit          -- Do gaussian fit? (default False)
         statx        -- Stat box x coordinate (default 0.65)
         staty        -- Stat box y coordinate (default 0.8)
         statyadjust  -- List of floats for stat box y coordinate adjustments (default None)
         normalizeToUnitArea -- Normalize histograms to unit area? (default False)
         normalizeToNumberOfEvents -- Normalize histograms to number of events? If yes, the PlotFolder needs 'numberOfEventsHistogram' set to a histogram filled once per event (default False)
         profileX     -- Take histograms via ProfileX()? (default False)
         fitSlicesY   -- Take histograms via FitSlicesY() (default False)
         rebinX       -- rebin x axis (default None)
         scale        -- Scale histograms by a number (default None)
         xbinlabels   -- List of x axis bin labels (if given, default None)
         xbinlabelsize -- Size of x axis bin labels (default None)
         xbinlabeloption -- Option string for x axis bin labels (default None)
         removeEmptyBins -- Bool for removing empty bins, but only if histogram has bin labels (default False)
         printBins    -- Bool for printing bin values, but only if histogram has bin labels (default False)
         drawStyle    -- If "hist", draw as line instead of points (default None)
         drawCommand  -- Deliver this to Draw() (default: None for same as drawStyle)
         lineWidth    -- If drawStyle=="hist", the width of line (default 2)
         legendDx     -- Float for moving TLegend in x direction for separate=True (default None)
         legendDy     -- Float for moving TLegend in y direction for separate=True (default None)
         legendDw     -- Float for changing TLegend width for separate=True (default None)
         legendDh     -- Float for changing TLegend height for separate=True (default None)
         legend       -- Bool to enable/disable legend (default True)
         adjustMarginRight  -- Float for adjusting right margin (default None)
         ratio        -- Possibility to disable ratio for this particular plot (default None)
         ratioYmin    -- Float for y axis minimum in ratio pad (default: list of values)
         ratioYmax    -- Float for y axis maximum in ratio pad (default: list of values)
         ratioFit     -- Fit straight line in ratio? (default None)
         ratioUncertainty -- Plot uncertainties on ratio? (default True)
         ratioCoverageXrange -- Range of x axis values (xmin,xmax) to limit the automatic ratio y axis range calculation to (default None for disabled)
         histogramModifier -- Function to be called in create() to modify the histograms (default None)
         """
         self._name = name
 
         def _set(attr, default):
             setattr(self, "_"+attr, kwargs.get(attr, default))
 
         _set("fallback", None)
         _set("outname", None)
 
         _set("title", None)
         _set("xtitle", None)
         _set("xtitlesize", None)
         _set("xtitleoffset", None)
         _set("xlabelsize", None)
         _set("ytitle", None)
         _set("ytitlesize", None)
         _set("ytitleoffset", None)
         _set("ztitle", None)
         _set("ztitleoffset", None)
 
         _set("xmin", None)
         _set("xmax", None)
         _set("ymin", 0.)
         _set("ymax", None)
 
         _set("xlog", False)
         _set("ylog", False)
         _set("xgrid", True)
         _set("ygrid", True)
 
         _set("stat", False)
         _set("fit", False)
 
         _set("statx", 0.65)
         _set("staty", 0.8)
         _set("statyadjust", None)
 
         _set("normalizeToUnitArea", False)
         _set("normalizeToNumberOfEvents", False)
         _set("profileX", False)
         _set("fitSlicesY", False)
         _set("rebinX", None)
 
         _set("scale", None)
         _set("xbinlabels", None)
         _set("xbinlabelsize", None)
         _set("xbinlabeloption", None)
         _set("removeEmptyBins", False)
         _set("printBins", False)
 
         _set("drawStyle", "EP")
         _set("drawCommand", None)
         _set("lineWidth", 2)
 
         _set("legendDx", None)
         _set("legendDy", None)
         _set("legendDw", None)
         _set("legendDh", None)
         _set("legend", True)
 
         _set("adjustMarginRight", None)
 
         _set("ratio", None)
         _set("ratioYmin", [0, 0.2, 0.5, 0.7, 0.8, 0.9, 0.95])
         _set("ratioYmax", [1.05, 1.1, 1.2, 1.3, 1.5, 1.8, 2, 2.5, 3, 4, 5])
         _set("ratioFit", None)
         _set("ratioUncertainty", True)
         _set("ratioCoverageXrange", None)
 
         _set("histogramModifier", None)
 
         self._histograms = []
 

Member Function Documentation

def plotting.Plot._createOne	(	self,
		name,
		index,
		tdir,
		nevents
	)

private

Create one histogram from a TDirectory.

Definition at line 1672 of file plotting.py.

References plotting._getOrCreateObject().

Referenced by plotting.Plot.create().

     def _createOne(self, name, index, tdir, nevents):
         """Create one histogram from a TDirectory."""
         if tdir == None:
             return None
 
         # If name is a list, pick the name by the index
         if isinstance(name, list):
             name = name[index]
 
         h = _getOrCreateObject(tdir, name)
         if h is not None and self._normalizeToNumberOfEvents and nevents is not None and nevents != 0:
             h.Scale(1.0/nevents)
         return h
 

def plotting.Plot._normalize ( self )

private

Normalise histograms to unit area

Definition at line 1777 of file plotting.py.

References RunHistogramManager._histograms, and plotting.Plot._histograms.

Referenced by plotting.Plot.draw().

     def _normalize(self):
         """Normalise histograms to unit area"""
 
         for h in self._histograms:
             if h is None:
                 continue
             i = h.Integral()
             if i == 0:
                 continue
             if h.GetSumw2().fN <= 0: # to suppress warning
                 h.Sumw2()
             h.Scale(1.0/i)
 

def plotting.Plot._setStats	(	self,
		histos,
		startingX,
		startingY
	)

private

Set stats box.

Definition at line 1736 of file plotting.py.

References DTSC._stat, and DTTrigGeom._stat.

     def _setStats(self, histos, startingX, startingY):
         """Set stats box."""
         if not self._stat:
             for h in histos:
                 if h is not None and hasattr(h, "SetStats"):
                     h.SetStats(0)
             return
 
         def _doStats(h, col, dy):
             if h is None:
                 return
             h.SetStats(True)
 
             if self._fit and h.GetEntries() > 0.5:
                 h.Fit("gaus", "Q")
                 f = h.GetListOfFunctions().FindObject("gaus")
                 if f == None:
                     h.SetStats(0)
                     return
                 f.SetLineColor(col)
                 f.SetLineWidth(1)
             h.Draw()
             ROOT.gPad.Update()
             st = h.GetListOfFunctions().FindObject("stats")
             if self._fit:
                 st.SetOptFit(0o010)
                 st.SetOptStat(1001)
             st.SetX1NDC(startingX)
             st.SetX2NDC(startingX+0.3)
             st.SetY1NDC(startingY+dy)
             st.SetY2NDC(startingY+dy+0.15)
             st.SetTextColor(col)
 
         dy = 0.0
         for i, h in enumerate(histos):
             if self._statyadjust is not None and i < len(self._statyadjust):
                 dy += self._statyadjust[i]
 
             _doStats(h, _plotStylesColor[i], dy)
             dy -= 0.19
 

def plotting.Plot.addToLegend	(	self,
		legend,
		legendLabels,
		denomUncertainty
	)

Add histograms to a legend.

Arguments:
legend       -- TLegend
legendLabels -- List of strings for the legend labels

Definition at line 2080 of file plotting.py.

References RunHistogramManager._histograms, plotting.Plot._histograms, and ComparisonHelper.zip().

     def addToLegend(self, legend, legendLabels, denomUncertainty):
         """Add histograms to a legend.
 
         Arguments:
         legend       -- TLegend
         legendLabels -- List of strings for the legend labels
         """
         first = denomUncertainty
         for h, label in zip(self._histograms, legendLabels):
             if h is None:
                 first = False
                 continue
             if first:
                 self._forLegend = h.Clone()
                 self._forLegend.SetFillStyle(1001)
                 self._forLegend.SetFillColor(ROOT.kGray)
                 entry = legend.AddEntry(self._forLegend, label, "lpf")
                 first = False
             else:
                 legend.AddEntry(h, label, "LP")
 

def plotting.Plot.clone	(	self,
		kwargs
	)

Definition at line 1634 of file plotting.py.

References GenericMVAComputerCache.isEmpty(), TkPhase2OTMeasurementDet.isEmpty(), TkPixelMeasurementDet.isEmpty(), MultipleKinematicConstraint.isEmpty(), KinematicTree.isEmpty(), reco::PFDisplacedVertexSeed.isEmpty(), L1MuBMTrackAssembler.isEmpty(), L1MuDTTrackAssembler.isEmpty(), BPHKinematicFit.isEmpty(), L1MuGMTInputEvent.isEmpty(), BXVector< T >.isEmpty(), TkStripMeasurementDet.isEmpty(), cond::persistency::IOVProxy.isEmpty(), HistoData.isEmpty, DDXMLElement.isEmpty(), html.Page.isEmpty(), and plotting.Plot.isEmpty().

     def clone(self, **kwargs):
         if not self.isEmpty():
             raise Exception("Plot can be cloned only before histograms have been created")
         cl = copy.copy(self)
         cl.setProperties(**kwargs)
         return cl
 

def plotting.Plot.create	(	self,
		tdirNEvents,
		requireAllHistograms = `False`
	)

Create histograms from list of TDirectories

Definition at line 1686 of file plotting.py.

References plotting.Plot._createOne(), RunHistogramManager._histograms, plotting.Plot._histograms, FP420HitsObject._name, TrackerHitsObject._name, PGeometricDet::Item._name, TrackingRecHitAlgorithm._name, LikelihoodSpecies._name, L1TMuon::PtAssignmentUnit._name, L1TMuon::PtRefinementUnit._name, LikelihoodPdfProduct._name, LikelihoodPdf._name, citk::IsolationConeDefinitionBase._name, Logger._name, hcaldqm::DQModule._name, DrellYanValidation._name, WValidation._name, HistoParams< T >._name, hcaldqm::flag::Flag._name, hcaldqm::quantity::Quantity._name, CutApplicatorBase._name, ElectronMVAEstimatorRun2Phys14NonTrig._name, PhotonMVAEstimatorRun2Spring15NonTrig._name, PhotonMVAEstimatorRun2Phys14NonTrig._name, ElectronMVAEstimatorRun2Spring15Trig._name, GeometricDetExtra._name, ElectronMVAEstimatorRun2Spring15NonTrig._name, HistoParams< TH2F >._name, HistoParams< TProfile2D >._name, Vispa.Views.PropertyView.Property._name, SequenceTypes.SequencePlaceholder._name, plotting.Subtract._name, plotting.Transform._name, plotting.FakeDuplicate._name, plotting.CutEfficiency._name, plotting.AggregateBins._name, plotting.AggregateHistos._name, plotting.ROC._name, plotting.Plot._name, Vispa.Gui.VispaWidget.VispaWidget._scale, plotting.AggregateBins._scale, HistoParams< TH1F >._title, HistoParams< T >._title, HistoParams< TProfile >._title, HistoParams< TH2F >._title, SiStripHitEffFromCalibTree._title, HistoParams< TProfile2D >._title, html.PageSet._title, plotting.Subtract._title, plotting.Transform._title, html.HtmlReport._title, plotting.FakeDuplicate._title, plotting.CutEfficiency._title, and genParticles_cff.map.

     def create(self, tdirNEvents, requireAllHistograms=False):
         """Create histograms from list of TDirectories"""
         self._histograms = [self._createOne(self._name, i, tdirNEvent[0], tdirNEvent[1]) for i, tdirNEvent in enumerate(tdirNEvents)]
 
         if self._fallback is not None:
             profileX = [self._profileX]*len(self._histograms)
             for i in xrange(0, len(self._histograms)):
                 if self._histograms[i] is None:
                     self._histograms[i] = self._createOne(self._fallback["name"], i, tdirNEvents[i][0], tdirNEvents[i][1])
                     profileX[i] = self._fallback.get("profileX", self._profileX)
 
         if self._histogramModifier is not None:
             self._histograms = self._histogramModifier(self._histograms)
 
         if len(self._histograms) > len(_plotStylesColor):
             raise Exception("More histograms (%d) than there are plot styles (%d) defined. Please define more plot styles in this file" % (len(self._histograms), len(_plotStylesColor)))
 
         # Modify histograms here in case self._name returns numbers
         # and self._histogramModifier creates the histograms from
         # these numbers
         def _modifyHisto(th1, profileX):
             if th1 is None:
                 return None
 
             if profileX:
                 th1 = th1.ProfileX()
 
             if self._fitSlicesY:
                 ROOT.TH1.AddDirectory(True)
                 th1.FitSlicesY()
                 th1 = ROOT.gDirectory.Get(th1.GetName()+"_2")
                 th1.SetDirectory(None)
                 #th1.SetName(th1.GetName()+"_ref")
                 ROOT.TH1.AddDirectory(False)
 
             if self._title is not None:
                 th1.SetTitle(self._title)
 
             if self._scale is not None:
                 th1.Scale(self._scale)
 
             return th1
 
         if self._fallback is not None:
             self._histograms = map(_modifyHisto, self._histograms, profileX)
         else:
             self._histograms = map(lambda h: _modifyHisto(h, self._profileX), self._histograms)
         if requireAllHistograms and None in self._histograms:
             self._histograms = [None]*len(self._histograms)
 

def plotting.Plot.draw	(	self,
		pad,
		ratio,
		ratioFactor,
		nrows
	)

Draw the histograms using values for a given algorithm.

Definition at line 1790 of file plotting.py.

References RunHistogramManager._histograms, plotting.Plot._histograms, plotting.Plot._normalize(), HistoData.getName(), emtf::Node.getName(), PhysicsTools::Source.getName(), PassiveHit.getName(), CocoaMaterialElementary.getName(), ecaldqm::DBReaderWorker.getName(), LikelihoodSpecies.getName(), PCaloHit.getName(), ParticleFlux.getName(), EcalLogicID.getName(), pftools::SpaceManager.getName(), ecaldqm::DBWriterWorker.getName(), LikelihoodPdf.getName(), TkOfflineVariables.getName(), hcaldqm::mapper::Mapper.getName(), FieldHandlerBase< TOutput >.getName(), AnyMVAEstimatorRun2Base.getName(), FittedEntry.getName(), TrackingRecHitAlgorithm.getName(), hcaldqm::mapper::HashMapper.getName(), jsoncollector::JsonMonitorable.getName(), CSGAction.getName(), L1TUtmAlgorithm.getName(), L1TUtmScale.getName(), DTTPGLutFile.getName(), PhysicsTools::MVATrainer.getName(), L1TriggerLutFile.getName(), SeedMergerPixelLayer.getName(), L1TUtmCondition.getName(), QCriterion.getName(), pos::PixelPortcardMap.getName(), L1MuGMTReg.getName(), GenericMVAComputer::TaggingVariableIterator< Iter_t >::Value.getName(), GEMStation.getName(), PhysicsTools::Variable::Value.getName(), pftools::SpaceVoxel.getName(), ecaldqm::DQWorker.getName(), RPCDCCLink.getName(), RPCAMCLink.getName(), PhotonMVAEstimatorRun2Spring16NonTrig.getName(), L1TUtmTriggerMenu.getName(), PhotonMVAEstimatorRun2Spring15NonTrig.getName(), L1TUtmObject.getName(), PhotonMVAEstimatorRun2Phys14NonTrig.getName(), L1TUtmCut.getName(), ElectronMVAEstimatorRun2Phys14NonTrig.getName(), PFTrackAlgoTools::edm::ELslProxy< ELgen >.getName(), edm::ELslProxy< ELgen >.getName(), L1MuGMTRegMMConfig.getName(), PhysicsTools::MVAModuleHelper< Record, Object, Filler >::Value.getName(), ElectronMVAEstimatorRun2Spring15Trig.getName(), ElectronMVAEstimatorRun2Spring15NonTrig.getName(), cscdqm::HistoDef.getName(), ElectronMVAEstimatorRun2Spring16GeneralPurpose.getName(), ElectronMVAEstimatorRun2Spring16HZZ.getName(), FourVectorHLT::PathInfo.getName(), RPCLBLink.getName(), jsoncollector::DataPoint.getName(), MonitorElement.getName(), PFTrackAlgoTools::edm::ELseverityLevel.getName(), edm::ELseverityLevel.getName(), PhysicsTools::Variable.getName(), PhysicsTools::TrainerMonitoring::Object.getName(), L1MuGMTRegSortRankOffset.getName(), L1MuGMTRegCDLConfig.getName(), cscdqm::CSCHistoDef.getName(), plotting.Plot.getName(), plotting.Plot.isTGraph2D(), and hpstanc_transforms.max.

     def draw(self, pad, ratio, ratioFactor, nrows):
         """Draw the histograms using values for a given algorithm."""
 #        if len(self._histograms) == 0:
 #            print "No histograms for plot {name}".format(name=self._name)
 #            return
 
         isTGraph2D = self.isTGraph2D()
         if isTGraph2D:
             # Ratios for the TGraph2Ds is not that interesting
             ratioOrig = ratio
             ratio = False
 
         if self._normalizeToUnitArea:
             self._normalize()
 
         if self._rebinX is not None:
             for h in self._histograms:
                 h.Rebin(self._rebinX)
 
         def _styleMarker(h, msty, col):
             h.SetMarkerStyle(msty)
             h.SetMarkerColor(col)
             h.SetMarkerSize(0.7)
             h.SetLineColor(1)
             h.SetLineWidth(1)
 
         def _styleHist(h, msty, col):
             _styleMarker(h, msty, col)
             h.SetLineColor(col)
             h.SetLineWidth(self._lineWidth)
 
         # Use marker or hist style
         style = _styleMarker
         if "hist" in self._drawStyle.lower():
             style = _styleHist
         if len(self._histograms) > 0 and isinstance(self._histograms[0], ROOT.TGraph):
             if "l" in self._drawStyle.lower():
                 style = _styleHist
 
         # Apply style to histograms, filter out Nones
         histos = []
         for i, h in enumerate(self._histograms):
             if h is None:
                 continue
             style(h, _plotStylesMarker[i], _plotStylesColor[i])
             histos.append(h)
         if len(histos) == 0:
             if verbose:
                 print "No histograms for plot {name}".format(name=self.getName())
             return
 
         # Extract x bin labels, make sure that only bins with same
         # label are compared with each other
         histosHaveBinLabels = len(histos[0].GetXaxis().GetBinLabel(1)) > 0
         xbinlabels = self._xbinlabels
         if xbinlabels is None:
             if histosHaveBinLabels:
                 xbinlabels = [histos[0].GetXaxis().GetBinLabel(i) for i in xrange(1, histos[0].GetNbinsX()+1)]
                 # Merge bin labels with difflib
                 for h in histos[1:]:
                     labels = [h.GetXaxis().GetBinLabel(i) for i in xrange(1, h.GetNbinsX()+1)]
                     diff = difflib.unified_diff(xbinlabels, labels, n=max(len(xbinlabels), len(labels)))
                     xbinlabels = []
                     operation = []
                     for item in diff: # skip the "header" lines
                         if item[:2] == "@@":
                             break
                     for item in diff:
                         operation.append(item[0])
                         lab = item[1:]
                         if lab in xbinlabels:
                             # pick the last addition of the bin
                             ind = xbinlabels.index(lab)
                             if operation[ind] == "-" and operation[-1] == "+":
                                 xbinlabels.remove(lab)
                                 del operation[ind] # to keep xbinlabels and operation indices in sync
                             elif operation[ind] == "+" and operation[-1] == "-":
                                 del operation[-1] # to keep xbinlabels and operation indices in sync
                                 continue
                             else:
                                 raise Exception("This should never happen")
                         xbinlabels.append(lab)
                     # unified_diff returns empty diff if xbinlabels and labels are equal
                     # so if xbinlabels is empty here, it can be just set to labels
                     if len(xbinlabels) == 0:
                         xbinlabels = labels
 
                 histos_new = []
                 for h in histos:
                     h_new = h.Clone(h.GetName()+"_xbinlabels")
                     h_new.SetBins(len(xbinlabels), h.GetBinLowEdge(1), h.GetBinLowEdge(1)+len(xbinlabels))
                     for i, label in enumerate(xbinlabels):
                         bin = h.GetXaxis().FindFixBin(label)
                         if bin >= 0:
                             h_new.SetBinContent(i+1, h.GetBinContent(bin))
                             h_new.SetBinError(i+1, h.GetBinError(bin))
                         else:
                             h_new.SetBinContent(i+1, 0)
                             h_new.SetBinError(i+1, 0)
                     histos_new.append(h_new)
                 self._tmp_histos = histos_new # need to keep these in memory too ...
                 histos = histos_new
 
         # Remove empty bins, but only if histograms have bin labels
         if self._removeEmptyBins and histosHaveBinLabels:
             # at this point, all histograms have been "equalized" by their x binning and labels
             # therefore remove bins which are empty in all histograms
             binsToRemove = set()
             for b in xrange(1, histos[0].GetNbinsX()+1):
                 binEmpty = True
                 for h in histos:
                     if h.GetBinContent(b) > 0:
                         binEmpty = False
                         break
                 if binEmpty:
                     binsToRemove.add(b)
 
             if len(binsToRemove) > 0:
                 # filter xbinlabels
                 xbinlab_new = []
                 for i in xrange(len(xbinlabels)):
                     if (i+1) not in binsToRemove:
                         xbinlab_new.append(xbinlabels[i])
                 xbinlabels = xbinlab_new
 
                 # filter histogram bins
                 histos_new = []
                 for h in histos:
                     values = []
                     for b in xrange(1, h.GetNbinsX()+1):
                         if b not in binsToRemove:
                             values.append( (h.GetXaxis().GetBinLabel(b), h.GetBinContent(b), h.GetBinError(b)) )
 
                     if len(values) > 0:
                         h_new = h.Clone(h.GetName()+"_empty")
                         h_new.SetBins(len(values), h.GetBinLowEdge(1), h.GetBinLowEdge(1)+len(values))
                         for b, (l, v, e) in enumerate(values):
                             h_new.GetXaxis().SetBinLabel(b+1, l)
                             h_new.SetBinContent(b+1, v)
                             h_new.SetBinError(b+1, e)
 
                         histos_new.append(h_new)
 
                 self._tmp_histos = histos_new # need to keep these in memory too ...
                 histos = histos_new
                 if len(histos) == 0:
                     if verbose:
                         print "No histograms with non-empty bins for plot {name}".format(name=self.getName())
                     return
 
         if self._printBins and histosHaveBinLabels:
             print "####################"
             print self._name
             width = max([len(l) for l in xbinlabels])
             tmp = "%%-%ds " % width
             for b in xrange(1, histos[0].GetNbinsX()+1):
                 s = tmp % xbinlabels[b-1]
                 for h in histos:
                     s += "%.3f " % h.GetBinContent(b)
                 print s
             print
 
         bounds = _findBounds(histos, self._ylog,
                              xmin=self._xmin, xmax=self._xmax,
                              ymin=self._ymin, ymax=self._ymax)
 
         # need to keep these in memory
         self._mainAdditional = []
         self._ratioAdditional = []
 
         if ratio:
             self._ratios = _calculateRatios(histos, self._ratioUncertainty) # need to keep these in memory too ...
             ratioHistos = filter(lambda h: h is not None, [r.getRatio() for r in self._ratios[1:]])
 
             if len(ratioHistos) > 0:
                 ratioBoundsY = _findBoundsY(ratioHistos, ylog=False, ymin=self._ratioYmin, ymax=self._ratioYmax, coverage=0.68, coverageRange=self._ratioCoverageXrange)
             else:
                 ratioBoundsY = (0.9, 1,1) # hardcoded default in absence of valid ratio calculations
 
             if self._ratioFit is not None:
                 for i, rh in enumerate(ratioHistos):
                     tf_line = ROOT.TF1("line%d"%i, "[0]+x*[1]")
                     tf_line.SetRange(self._ratioFit["rangemin"], self._ratioFit["rangemax"])
                     fitres = rh.Fit(tf_line, "RINSQ")
                     tf_line.SetLineColor(rh.GetMarkerColor())
                     tf_line.SetLineWidth(2)
                     self._ratioAdditional.append(tf_line)
                     box = PlotTextBox(xmin=self._ratioFit.get("boxXmin", 0.14), ymin=None, # None for automatix
                                       xmax=self._ratioFit.get("boxXmax", 0.35), ymax=self._ratioFit.get("boxYmax", 0.09),
                                       color=rh.GetMarkerColor(), font=43, size=11, lineheight=0.02)
                     box.move(dx=(box.width()+0.01)*i)
                     #box.addText("Const: %.4f" % fitres.Parameter(0))
                     #box.addText("Slope: %.4f" % fitres.Parameter(1))
                     box.addText("Const: %.4f#pm%.4f" % (fitres.Parameter(0), fitres.ParError(0)))
                     box.addText("Slope: %.4f#pm%.4f" % (fitres.Parameter(1), fitres.ParError(1)))
                     self._mainAdditional.append(box)
 
 
         # Create bounds before stats in order to have the
         # SetRangeUser() calls made before the fit
         #
         # stats is better to be called before frame, otherwise get
         # mess in the plot (that frame creation cleans up)
         if ratio:
             pad.cd(1)
         self._setStats(histos, self._statx, self._staty)
 
         # Create frame
         if isTGraph2D:
             frame = FrameTGraph2D(pad, bounds, histos, ratioOrig, ratioFactor)
         else:
             if ratio:
                 ratioBounds = (bounds[0], ratioBoundsY[0], bounds[2], ratioBoundsY[1])
                 frame = FrameRatio(pad, bounds, ratioBounds, ratioFactor, nrows, xbinlabels, self._xbinlabelsize, self._xbinlabeloption)
             else:
                 frame = Frame(pad, bounds, nrows, xbinlabels, self._xbinlabelsize, self._xbinlabeloption)
 
         # Set log and grid
         frame.setLogx(self._xlog)
         frame.setLogy(self._ylog)
         frame.setGridx(self._xgrid)
         frame.setGridy(self._ygrid)
 
         # Construct draw option string
         opt = "sames" # s for statbox or something?
         ds = ""
         if self._drawStyle is not None:
             ds = self._drawStyle
         if self._drawCommand is not None:
             ds = self._drawCommand
         if len(ds) > 0:
             opt += " "+ds
 
         # Set properties of frame
         frame.setTitle(histos[0].GetTitle())
         if self._xtitle is not None:
             frame.setXTitle(self._xtitle)
         if self._xtitlesize is not None:
             frame.setXTitleSize(self._xtitlesize)
         if self._xtitleoffset is not None:
             frame.setXTitleOffset(self._xtitleoffset)
         if self._xlabelsize is not None:
             frame.setXLabelSize(self._xlabelsize)
         if self._ytitle is not None:
             frame.setYTitle(self._ytitle)
         if self._ytitlesize is not None:
             frame.setYTitleSize(self._ytitlesize)
         if self._ytitleoffset is not None:
             frame.setYTitleOffset(self._ytitleoffset)
         if self._ztitle is not None:
             frame.setZTitle(self._ztitle)
         if self._ztitleoffset is not None:
             frame.setZTitleOffset(self._ztitleoffset)
         if self._adjustMarginRight is not None:
             frame.adjustMarginRight(self._adjustMarginRight)
         elif "z" in opt:
             frame.adjustMarginLeft(0.03)
             frame.adjustMarginRight(0.08)
 
         # Draw histograms
         if ratio:
             frame._pad.cd()
 
         for h in histos:
             h.Draw(opt)
 
         for addl in self._mainAdditional:
             addl.Draw("same")
 
         # Draw ratios
         if ratio and len(histos) > 0:
             frame._padRatio.cd()
             firstRatio = self._ratios[0].getRatio()
             if self._ratioUncertainty and firstRatio is not None:
                 firstRatio.SetFillStyle(1001)
                 firstRatio.SetFillColor(ROOT.kGray)
                 firstRatio.SetLineColor(ROOT.kGray)
                 firstRatio.SetMarkerColor(ROOT.kGray)
                 firstRatio.SetMarkerSize(0)
                 self._ratios[0].draw("E2")
                 frame._padRatio.RedrawAxis("G") # redraw grid on top of the uncertainty of denominator
             for r in self._ratios[1:]:
                 r.draw()
 
             for addl in self._ratioAdditional:
                 addl.Draw("same")
 
         frame.redrawAxis()
         self._frame = frame # keep the frame in memory for sure
 

def plotting.Plot.drawRatioUncertainty ( self )

Return true if the ratio uncertainty should be drawn

Definition at line 1668 of file plotting.py.

     def drawRatioUncertainty(self):
         """Return true if the ratio uncertainty should be drawn"""
         return self._ratioUncertainty
 

def plotting.Plot.getName ( self )

Definition at line 1660 of file plotting.py.

References FP420HitsObject._name, TrackerHitsObject._name, PGeometricDet::Item._name, TrackingRecHitAlgorithm._name, LikelihoodSpecies._name, L1TMuon::PtAssignmentUnit._name, L1TMuon::PtRefinementUnit._name, LikelihoodPdfProduct._name, LikelihoodPdf._name, citk::IsolationConeDefinitionBase._name, Logger._name, hcaldqm::DQModule._name, DrellYanValidation._name, WValidation._name, HistoParams< T >._name, hcaldqm::flag::Flag._name, hcaldqm::quantity::Quantity._name, CutApplicatorBase._name, ElectronMVAEstimatorRun2Phys14NonTrig._name, PhotonMVAEstimatorRun2Spring15NonTrig._name, PhotonMVAEstimatorRun2Phys14NonTrig._name, ElectronMVAEstimatorRun2Spring15Trig._name, GeometricDetExtra._name, ElectronMVAEstimatorRun2Spring15NonTrig._name, HistoParams< TH2F >._name, HistoParams< TProfile2D >._name, Vispa.Views.PropertyView.Property._name, SequenceTypes.SequencePlaceholder._name, plotting.Subtract._name, plotting.Transform._name, plotting.FakeDuplicate._name, plotting.CutEfficiency._name, plotting.AggregateBins._name, plotting.AggregateHistos._name, plotting.ROC._name, plotting.Plot._name, and harvestTrackValidationPlots.str.

Referenced by plotting.Plot.draw().

     def getName(self):
         if self._outname is not None:
             return self._outname
         if isinstance(self._name, list):
             return str(self._name[0])
         else:
             return str(self._name)
 

def plotting.Plot.getNumberOfHistograms ( self )

Return number of existing histograms.

Definition at line 1641 of file plotting.py.

References RunHistogramManager._histograms, plotting.Plot._histograms, and ALCARECOTkAlBeamHalo_cff.filter.

Referenced by plotting.Plot.isEmpty().

     def getNumberOfHistograms(self):
         """Return number of existing histograms."""
         return len(filter(lambda h: h is not None, self._histograms))
 

def plotting.Plot.isEmpty ( self )

Return true if there are no histograms created for the plot

Definition at line 1645 of file plotting.py.

References plotting.Plot.getNumberOfHistograms().

Referenced by plotting.Plot.clone().

     def isEmpty(self):
         """Return true if there are no histograms created for the plot"""
         return self.getNumberOfHistograms() == 0
 

def plotting.Plot.isRatio	(	self,
		ratio
	)

Definition at line 1655 of file plotting.py.

     def isRatio(self, ratio):
         if self._ratio is None:
             return ratio
         return ratio and self._ratio
 

def plotting.Plot.isTGraph2D ( self )

Definition at line 1649 of file plotting.py.

References RunHistogramManager._histograms, and plotting.Plot._histograms.

Referenced by plotting.Plot.draw().

     def isTGraph2D(self):
         for h in self._histograms:
             if isinstance(h, ROOT.TGraph2D):
                 return True
         return False
 

def plotting.Plot.setProperties	(	self,
		kwargs
	)

Definition at line 1628 of file plotting.py.

     def setProperties(self, **kwargs):
         for name, value in kwargs.iteritems():
             if not hasattr(self, "_"+name):
                 raise Exception("No attribute '%s'" % name)
             setattr(self, "_"+name, value)
 

Member Data Documentation

plotting.Plot._forLegend

private

Definition at line 2093 of file plotting.py.

plotting.Plot._frame

private

Definition at line 2078 of file plotting.py.

plotting.Plot._histograms

private

Definition at line 1626 of file plotting.py.

Referenced by plotting.Plot._normalize(), plotting.Plot.addToLegend(), plotting.Plot.create(), plotting.Plot.draw(), plotting.Plot.getNumberOfHistograms(), and plotting.Plot.isTGraph2D().

plotting.Plot._mainAdditional

private

Definition at line 1957 of file plotting.py.

plotting.Plot._name

private

Definition at line 1556 of file plotting.py.

Referenced by trackingPlots.TimePerTrackPlot.__str__(), plotting.PlotGroup._drawSeparate(), plotting.PlotGroup._save(), trackingPlots.TimePerTrackPlot.create(), plotting.Plot.create(), plotting.PlotGroup.draw(), plotting.Plot.getName(), plotting.PlotGroup.getName(), plotting.PlotterFolder.getName(), plotting.PlotterItem.getName(), plotting.PlotterFolder.getSelectionNameIterator(), validation.SimpleSample.name(), trackingPlots.Iteration.name(), plotting.PlotterItem.readDirs(), and plotting.PlotGroup.remove().

plotting.Plot._ratioAdditional

private

Definition at line 1958 of file plotting.py.

plotting.Plot._ratios

private

Definition at line 1961 of file plotting.py.

plotting.Plot._tmp_histos

private

Definition at line 1890 of file plotting.py.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation