utils::Chi2 Class Reference

Inheritance diagram for utils::Chi2:

List of all members.

Public Member Functions
def	__init__
def	absval
def	check_filled_bins
def	check_histograms
def	do_test
Public Attributes
	name

---

Detailed Description

Definition at line 238 of file utils.py.

---

Constructor & Destructor Documentation

def utils::Chi2::__init__	(		self,
			threshold
	)

Reimplemented from utils::StatisticalTest.

Definition at line 239 of file utils.py.

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

---

Member Function Documentation

def utils::Chi2::absval

(

self

)

Definition at line 254 of file utils.py.

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

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

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

def utils::Chi2::do_test

(

self

)

Reimplemented from utils::StatisticalTest.

Definition at line 272 of file utils.py.

                   :
    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

Reimplemented from utils::StatisticalTest.

Definition at line 239 of file utils.py.