utils.BinToBin Class Reference

Inheritance diagram for utils.BinToBin:

Public Member Functions
def	__init__ (self, threshold=1)
def	checkBinningMatches (self)
def	do_test (self)
Public Member Functions inherited from utils.StatisticalTest
def	__init__ (self, threshold)
def	do_test (self)
def	get_rank (self)
def	get_status (self)
def	set_operands (self, h1, h2)

Public Attributes
	epsilon
	name
Public Attributes inherited from utils.StatisticalTest
	h1
	h2
	is_init
	name
	rank
	2D! return test_codes["2D"] More...
	threshold

Detailed Description

The bin to bin comparison builds a fake pvalue. It is 0 if the number of 
bins is different. It is % of corresponding bins otherwhise.
A threshold of 1 is needed to require a 1 to 1 correspondance between 
hisograms.

Definition at line 298 of file utils.py.

Constructor & Destructor Documentation

def utils.BinToBin.__init__	(		self,
			threshold = 1
	)

Definition at line 304 of file utils.py.

  def __init__(self, threshold=1):
    StatisticalTest.__init__(self, threshold)
    self.name='BinToBin'
    self.epsilon= 0.000001

Member Function Documentation

def utils.BinToBin.checkBinningMatches

(

self

)

Definition at line 309 of file utils.py.

References funct.abs(), learn_.epsilon, HcalAutoPedestalValidator.epsilon, HcalPedestalWidthsCheck.epsilon, HcalPedestalsCheck.epsilon, HcalGainsCheck.epsilon, and utils.BinToBin.epsilon.

Referenced by utils.BinToBin.do_test(), and utils.BinToBin1percent.do_test().

  def checkBinningMatches(self):
    if self.h1.GetNbinsX() != self.h2.GetNbinsX() \
           or self.h1.GetNbinsY() != self.h2.GetNbinsY() \
           or self.h1.GetNbinsZ() != self.h2.GetNbinsZ() \
           or abs(self.h1.GetXaxis().GetXmin() - self.h2.GetXaxis().GetXmin()) >self.epsilon \
           or abs(self.h1.GetYaxis().GetXmin() - self.h2.GetYaxis().GetXmin()) >self.epsilon \
           or abs(self.h1.GetZaxis().GetXmin() - self.h2.GetZaxis().GetXmin()) >self.epsilon \
           or abs(self.h1.GetXaxis().GetXmax() - self.h2.GetXaxis().GetXmax()) >self.epsilon \
           or abs(self.h1.GetYaxis().GetXmax() - self.h2.GetYaxis().GetXmax()) >self.epsilon \
           or abs(self.h1.GetZaxis().GetXmax() - self.h2.GetZaxis().GetXmax()) >self.epsilon:
      return False
    return True

def utils.BinToBin.do_test

(

self

)

Definition at line 322 of file utils.py.

References funct.abs(), utils.BinToBin.checkBinningMatches(), DQMStore.checkBinningMatches(), learn_.epsilon, HcalAutoPedestalValidator.epsilon, HcalPedestalWidthsCheck.epsilon, HcalPedestalsCheck.epsilon, HcalGainsCheck.epsilon, utils.BinToBin.epsilon, utils.getNbins(), FitSlicesYTool.h1, ConfigurableHisto.h1, DQMSourceExample.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, and edm.print().

  def do_test(self):
    # fist check that binning matches
    if not self.checkBinningMatches():
      return test_codes["DIFF_BIN"]
    # then do the real check
    equal = 1
    nbins = getNbins(self.h1)
    n_ok_bins=0.0
    for ibin in xrange(0, nbins+2):
      h1bin=self.h1.GetBinContent(ibin)
      h2bin=self.h2.GetBinContent(ibin)
      bindiff=h1bin-h2bin

      binavg=.5*(h1bin+h2bin)

      if binavg==0 or abs(bindiff) < self.epsilon:
        n_ok_bins+=1
        #print "Bin %ibin: bindiff %s" %(ibin,bindiff)
      else:
        print("Bin %ibin: bindiff %s" %(ibin,bindiff))

      #if abs(bindiff)!=0 :
        #print "Bin %ibin: bindiff %s" %(ibin,bindiff)
    
    rank=n_ok_bins/nbins
    
    if rank!=1:
      print("Histogram %s differs: nok: %s ntot: %s" %(self.h1.GetName(),n_ok_bins,nbins))
    
    return rank

#-------------------------------------------------------------------------------

Member Data Documentation

utils.BinToBin.epsilon

Definition at line 307 of file utils.py.

Referenced by utils.BinToBin.checkBinningMatches(), utils.BinToBin1percent.checkBinningMatches(), and utils.BinToBin.do_test().

utils.BinToBin.name

Definition at line 306 of file utils.py.

Referenced by config.CFG.__str__(), validation.Sample.digest(), VIDSelectorBase.VIDSelectorBase.initialize(), and Vispa.Views.PropertyView.Property.valueChanged().