utils.Chi2 Class Reference

Inheritance diagram for utils.Chi2:

Public Member Functions
def	__init__ (self, threshold)
def	absval (self)
def	check_filled_bins (self, min_filled)
def	check_histograms (self, histogram)
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
	name
Public Attributes inherited from utils.StatisticalTest
	h1
	h2
	is_init
	name
	rank
	2D! return test_codes["2D"] More...
	threshold

Detailed Description

Definition at line 234 of file utils.py.

Constructor & Destructor Documentation

def utils.Chi2.__init__	(		self,
			threshold
	)

Definition at line 235 of file utils.py.

  def __init__(self, threshold):
    StatisticalTest.__init__(self,threshold)
    self.name="Chi2"     

Member Function Documentation

def utils.Chi2.absval

(

self

)

Definition at line 250 of file utils.py.

References utils.getNbins(), FitSlicesYTool.h1, ConfigurableHisto.h1, DQMSourceExample.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h2, EgammaTowerIsolationNew< NC >::Sum.h2, ConfigurableHisto.h2, EgammaTowerIsolationNew< NC >.h2, DQMSourceExample.h2, utils.StatisticalTest.h2, and PlotAlignmentValidation::DMRPlotInfo.h2.

Referenced by utils.Chi2.do_test().

  def absval(self):
    nbins=getNbins(self.h1)
    binc=0
    for i in xrange(1,nbins):
      for h in self.h1,self.h2:
        binc=h.GetBinContent(i)
        if binc<0:
          h.SetBinContent(i,-1*binc)
        if h.GetBinError(i)==0 and binc!=0:
          #print "Histo ",h.GetName()," Bin:",i,"-Content:",h.GetBinContent(i)," had zero error"
          h.SetBinContent(i,0)

def utils.Chi2.check_filled_bins	(		self,
			min_filled
	)

Definition at line 239 of file utils.py.

References FitSlicesYTool.h1, ConfigurableHisto.h1, DQMSourceExample.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h2, EgammaTowerIsolationNew< NC >::Sum.h2, ConfigurableHisto.h2, EgammaTowerIsolationNew< NC >.h2, DQMSourceExample.h2, utils.StatisticalTest.h2, and PlotAlignmentValidation::DMRPlotInfo.h2.

Referenced by utils.Chi2.do_test().

  def check_filled_bins(self,min_filled):
    nbins=self.h1.GetNbinsX()
    n_filled_l=[]
    for h in self.h1,self.h2:
      nfilled=0.
      for ibin in xrange(1,nbins+1):
        if h.GetBinContent(ibin)>0:
          nfilled+=1
      n_filled_l.append(nfilled)
    return len(filter (lambda x:x>=min_filled,n_filled_l) )>0

def utils.Chi2.check_histograms	(		self,
			histogram
	)

Definition at line 262 of file utils.py.

Referenced by utils.Chi2.do_test().

  def check_histograms(self, histogram):
      if histogram.InheritsFrom("TProfile") or  (histogram.GetEntries()!=histogram.GetSumOfWeights()):
          return 'W'
      else:
          return 'U'

def utils.Chi2.do_test

(

self

)

Definition at line 268 of file utils.py.

References utils.Chi2.absval(), utils.Chi2.check_filled_bins(), utils.Chi2.check_histograms(), FitSlicesYTool.h1, ConfigurableHisto.h1, DQMSourceExample.h1, utils.StatisticalTest.h1, PlotAlignmentValidation::DMRPlotInfo.h1, FitSlicesYTool.h2, EgammaTowerIsolationNew< NC >::Sum.h2, ConfigurableHisto.h2, EgammaTowerIsolationNew< NC >.h2, DQMSourceExample.h2, utils.StatisticalTest.h2, and PlotAlignmentValidation::DMRPlotInfo.h2.

  def do_test(self):
    self.absval()
    if self.check_filled_bins(3):
      #if self.h1.InheritsFrom("TProfile") or  (self.h1.GetEntries()!=self.h1.GetSumOfWeights()):
      #  chi2=self.h1.Chi2Test(self.h2,'WW')
      #  #if chi2==0: print "DEBUG",self.h1.GetName(),"Chi2 is:", chi2
      #  return chi2
      #else:
      #  return self.h1.Chi2Test(self.h2,'UU')
      hist1 = self.check_histograms(self.h1)
      hist2 = self.check_histograms(self.h2)
      if hist1 =='W' and hist2 =='W': ##in case 
          chi2 = self.h1.Chi2Test(self.h2,'WW')     ## the both histograms are weighted
          return chi2
      elif hist1 == 'U' and hist2 == 'U':
          chi2 = self.h1.Chi2Test(self.h2,'UU')    ##the both histograms are unweighted
          return chi2
      elif hist1 == 'U' and hist2 == 'W':
          chi2 = self.h1.Chi2Test(self.h2,'UW')   ## 1st histogram is unweighted, 2nd weighted
          return chi2
      elif hist1 == 'W' and hist2 == 'U':
          chi2 = self.h2.Chi2Test(self.h1,'UW')   ## 1 is wieghted, 2nd unweigthed. so flip order to make a UW comparison
          return chi2
    else:
      return 1
      #return test_codes["FEW_BINS"]

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

Member Data Documentation

utils.Chi2.name

Definition at line 237 of file utils.py.

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