CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Public Attributes
utils.BinToBin1percent Class Reference
Inheritance diagram for utils.BinToBin1percent:
utils.StatisticalTest

Public Member Functions

def __init__
 
def checkBinningMatches
 
def do_test
 
- Public Member Functions inherited from utils.StatisticalTest
def __init__
 
def do_test
 
def get_rank
 
def get_status
 
def set_operands
 

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

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

Definition at line 372 of file utils.py.

372 
373  def __init__(self, threshold=1):
374  StatisticalTest.__init__(self, threshold)
375  self.name='BinToBin1percent'
376  self.epsilon= 0.000001
377  self.tolerance= 0.01

Member Function Documentation

def utils.BinToBin1percent.checkBinningMatches (   self)

Definition at line 378 of file utils.py.

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

Referenced by utils.BinToBin1percent.do_test().

378 
379  def checkBinningMatches(self):
380  if self.h1.GetNbinsX() != self.h2.GetNbinsX() \
381  or self.h1.GetNbinsY() != self.h2.GetNbinsY() \
382  or self.h1.GetNbinsZ() != self.h2.GetNbinsZ() \
383  or abs(self.h1.GetXaxis().GetXmin() - self.h2.GetXaxis().GetXmin()) >self.epsilon \
384  or abs(self.h1.GetYaxis().GetXmin() - self.h2.GetYaxis().GetXmin()) >self.epsilon \
385  or abs(self.h1.GetZaxis().GetXmin() - self.h2.GetZaxis().GetXmin()) >self.epsilon \
386  or abs(self.h1.GetXaxis().GetXmax() - self.h2.GetXaxis().GetXmax()) >self.epsilon \
387  or abs(self.h1.GetYaxis().GetXmax() - self.h2.GetYaxis().GetXmax()) >self.epsilon \
388  or abs(self.h1.GetZaxis().GetXmax() - self.h2.GetZaxis().GetXmax()) >self.epsilon:
389  return False
390  return True
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
def utils.BinToBin1percent.do_test (   self)

Definition at line 391 of file utils.py.

References funct.abs(), utils.BinToBin.checkBinningMatches(), utils.BinToBin1percent.checkBinningMatches(), utils.getNbins(), ConfigurableHisto.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h1, print(), sistrip::SpyUtilities.range(), MuonGeometrySanityCheck.tolerance, and utils.BinToBin1percent.tolerance.

391 
392  def do_test(self):
393  # fist check that binning matches
394  if not self.checkBinningMatches():
395  return test_codes["DIFF_BIN"]
396  # then do the real check
397  equal = 1
398  nbins = getNbins(self.h1)
399  n_ok_bins=0.0
400  for ibin in range(0,nbins):
401  ibin+=1
402  h1bin=self.h1.GetBinContent(ibin)
403  h2bin=self.h2.GetBinContent(ibin)
404  bindiff=h1bin-h2bin
405 
406  binavg=.5*(h1bin+h2bin)
407 
408  if binavg==0 or 100*abs(bindiff)/binavg < self.tolerance:
409  n_ok_bins+=1
410  #print "Bin %i bin: bindiff %s" %(ibin,bindiff)
411  else:
412  print("-->Bin %i bin: bindiff %s (%s - %s )" %(ibin,bindiff,h1bin,h2bin))
413 
414  #if abs(bindiff)!=0 :
415  #print "Bin %ibin: bindiff %s" %(ibin,bindiff)
416 
417  rank=n_ok_bins/nbins
418 
419  if rank!=1:
420  print("%s nok: %s ntot: %s" %(self.h1.GetName(),n_ok_bins,nbins))
421 
422  return rank
423 #-------------------------------------------------------------------------------
424 Statistical_Tests={"KS":KS,
425  "Chi2":Chi2,
426  "BinToBin":BinToBin,
427  "BinToBin1percent":BinToBin1percent,
428  "Bin2Bin":BinToBin,
429  "b2b":BinToBin,}
430 #-------------------------------------------------------------------------------
sistrip::SpyUtilities::range
const uint16_t range(const Frame &aFrame)
Definition: SiStripSpyUtilities.cc:76
print
void print(TMatrixD &m, const char *label=nullptr, bool mathematicaFormat=false)
Definition: Utilities.cc:47
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
utils.getNbins
def getNbins
Definition: utils.py:94

Member Data Documentation

utils.BinToBin1percent.epsilon

Definition at line 375 of file utils.py.

Referenced by utils.BinToBin1percent.checkBinningMatches().

utils.BinToBin1percent.name

Definition at line 374 of file utils.py.

Referenced by config.CFG.__str__(), geometryComparison.GeometryComparison.createScript(), validation.Sample.digest(), primaryVertexResolution.PrimaryVertexResolution.getRepMap(), primaryVertexValidation.PrimaryVertexValidation.getRepMap(), zMuMuValidation.ZMuMuValidation.getRepMap(), and VIDSelectorBase.VIDSelectorBase.initialize().

utils.BinToBin1percent.tolerance

Definition at line 376 of file utils.py.

Referenced by utils.BinToBin1percent.do_test().

Generated for CMSSW Reference Manual by doxygen 1.8.5