utils.BinToBin1percent Class Reference

Inheritance diagram for utils.BinToBin1percent:

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
	tolerance
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 366 of file utils.py.

Constructor & Destructor Documentation

__init__()

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

Definition at line 372 of file utils.py.

References HcalAutoPedestalValidator.epsilon, HcalPedestalWidthsCheck.epsilon, HcalPedestalsCheck.epsilon, HcalGainsCheck.epsilon, utils.BinToBin.epsilon, utils.BinToBin1percent.epsilon, MuonGeometrySanityCheck.tolerance, and utils.BinToBin1percent.tolerance.

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

Member Function Documentation

checkBinningMatches()

def utils.BinToBin1percent.checkBinningMatches

(

self

)

Definition at line 378 of file utils.py.

References funct.abs(), utils.BinToBin.checkBinningMatches(), utils.BinToBin1percent.checkBinningMatches(), utils.BinToBin1percent.do_test(), HcalAutoPedestalValidator.epsilon, HcalPedestalsCheck.epsilon, HcalPedestalWidthsCheck.epsilon, HcalGainsCheck.epsilon, utils.BinToBin.epsilon, utils.BinToBin1percent.epsilon, ConfigurableHisto.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h1, EgammaTowerIsolationNew< NC >::Sum.h2, ConfigurableHisto.h2, EgammaTowerIsolationNew< NC >.h2, utils.StatisticalTest.h2, PlotAlignmentValidation::DMRPlotInfo.h2, and FitSlicesYTool.h2.

Referenced by utils.BinToBin1percent.checkBinningMatches().

  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

do_test()

def utils.BinToBin1percent.do_test

(

self

)

Definition at line 391 of file utils.py.

References funct.abs(), utils.ask_ok(), utils.getNbins(), ConfigurableHisto.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h1, EgammaTowerIsolationNew< NC >::Sum.h2, ConfigurableHisto.h2, EgammaTowerIsolationNew< NC >.h2, utils.StatisticalTest.h2, PlotAlignmentValidation::DMRPlotInfo.h2, FitSlicesYTool.h2, print(), FastTimerService_cff.range, MuonGeometrySanityCheck.tolerance, and utils.BinToBin1percent.tolerance.

Referenced by utils.BinToBin1percent.checkBinningMatches().

  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 range(0,nbins):
      ibin+=1
      h1bin=self.h1.GetBinContent(ibin)
      h2bin=self.h2.GetBinContent(ibin)
      bindiff=h1bin-h2bin

      binavg=.5*(h1bin+h2bin)

      if binavg==0 or 100*abs(bindiff)/binavg < self.tolerance:
        n_ok_bins+=1
        #print "Bin %i bin: bindiff %s" %(ibin,bindiff)
      else:
        print("-->Bin %i bin: bindiff %s (%s - %s )" %(ibin,bindiff,h1bin,h2bin))

      #if abs(bindiff)!=0 :
        #print "Bin %ibin: bindiff %s" %(ibin,bindiff)
    
    rank=n_ok_bins/nbins
    
    if rank!=1:
      print("%s nok: %s ntot: %s" %(self.h1.GetName(),n_ok_bins,nbins))
    
    return rank
#-------------------------------------------------------------------------------
Statistical_Tests={"KS":KS,
                   "Chi2":Chi2,
                   "BinToBin":BinToBin,
                   "BinToBin1percent":BinToBin1percent,
                   "Bin2Bin":BinToBin,
                   "b2b":BinToBin,}
#-------------------------------------------------------------------------------  

Member Data Documentation

epsilon

utils.BinToBin1percent.epsilon

Definition at line 375 of file utils.py.

Referenced by utils.BinToBin1percent.__init__(), and utils.BinToBin1percent.checkBinningMatches().

name

utils.BinToBin1percent.name

Definition at line 374 of file utils.py.

Referenced by config.CFG.__str__(), validation.Sample.digest(), and VIDSelectorBase.VIDSelectorBase.initialize().

tolerance

utils.BinToBin1percent.tolerance

Definition at line 376 of file utils.py.

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