utils::BinToBin Class Reference

Inheritance diagram for utils::BinToBin:

List of all members.

Public Member Functions
def	__init__
def	checkBinningMatches
def	do_test
Public Attributes
	epsilon
	name

---

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

---

Constructor & Destructor Documentation

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

Reimplemented from utils::StatisticalTest.

Definition at line 288 of file utils.py.

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

---

Member Function Documentation

def utils::BinToBin::checkBinningMatches

(

self

)

Definition at line 293 of file utils.py.

                               :
    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

)

Reimplemented from utils::StatisticalTest.

Definition at line 306 of file utils.py.

                   :
    # 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):
      ibin+=1
      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 288 of file utils.py.

utils::BinToBin::name

Reimplemented from utils::StatisticalTest.

Definition at line 288 of file utils.py.