CMS 3D CMS Logo

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

Public Member Functions

def __init__
 
def absval
 
def check_filled_bins
 
def check_histograms
 
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

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

Constructor & Destructor Documentation

def utils.Chi2.__init__ (   self,
  threshold 
)

Definition at line 236 of file utils.py.

237  def __init__(self, threshold):
238  StatisticalTest.__init__(self,threshold)
239  self.name="Chi2"
def __init__
Definition: utils.py:236

Member Function Documentation

def utils.Chi2.absval (   self)

Definition at line 251 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().

252  def absval(self):
253  nbins=getNbins(self.h1)
254  binc=0
255  for i in xrange(1,nbins):
256  for h in self.h1,self.h2:
257  binc=h.GetBinContent(i)
258  if binc<0:
259  h.SetBinContent(i,-1*binc)
260  if h.GetBinError(i)==0 and binc!=0:
261  #print "Histo ",h.GetName()," Bin:",i,"-Content:",h.GetBinContent(i)," had zero error"
262  h.SetBinContent(i,0)
def absval
Definition: utils.py:251
def getNbins
Definition: utils.py:89
def utils.Chi2.check_filled_bins (   self,
  min_filled 
)

Definition at line 240 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().

241  def check_filled_bins(self,min_filled):
242  nbins=self.h1.GetNbinsX()
243  n_filled_l=[]
244  for h in self.h1,self.h2:
245  nfilled=0.
246  for ibin in xrange(1,nbins+1):
247  if h.GetBinContent(ibin)>0:
248  nfilled+=1
249  n_filled_l.append(nfilled)
250  return len(filter (lambda x:x>=min_filled,n_filled_l) )>0
def check_filled_bins
Definition: utils.py:240
def utils.Chi2.check_histograms (   self,
  histogram 
)

Definition at line 263 of file utils.py.

Referenced by utils.Chi2.do_test().

264  def check_histograms(self, histogram):
265  if histogram.InheritsFrom("TProfile") or (histogram.GetEntries()!=histogram.GetSumOfWeights()):
266  return 'W'
267  else:
268  return 'U'
def check_histograms
Definition: utils.py:263
def utils.Chi2.do_test (   self)

Definition at line 269 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.

270  def do_test(self):
271  self.absval()
272  if self.check_filled_bins(3):
273  #if self.h1.InheritsFrom("TProfile") or (self.h1.GetEntries()!=self.h1.GetSumOfWeights()):
274  # chi2=self.h1.Chi2Test(self.h2,'WW')
275  # #if chi2==0: print "DEBUG",self.h1.GetName(),"Chi2 is:", chi2
276  # return chi2
277  #else:
278  # return self.h1.Chi2Test(self.h2,'UU')
279  hist1 = self.check_histograms(self.h1)
280  hist2 = self.check_histograms(self.h2)
281  if hist1 =='W' and hist2 =='W': ##in case
282  chi2 = self.h1.Chi2Test(self.h2,'WW') ## the both histograms are weighted
283  return chi2
284  elif hist1 == 'U' and hist2 == 'U':
285  chi2 = self.h1.Chi2Test(self.h2,'UU') ##the both histograms are unweighted
286  return chi2
287  elif hist1 == 'U' and hist2 == 'W':
288  chi2 = self.h1.Chi2Test(self.h2,'UW') ## 1st histogram is unweighted, 2nd weighted
289  return chi2
290  elif hist1 == 'W' and hist2 == 'U':
291  chi2 = self.h2.Chi2Test(self.h1,'UW') ## 1 is wieghted, 2nd unweigthed. so flip order to make a UW comparison
292  return chi2
293  else:
294  return 1
295  #return test_codes["FEW_BINS"]
296 
297 #-------------------------------------------------------------------------------
def do_test
Definition: utils.py:269
def absval
Definition: utils.py:251
def check_filled_bins
Definition: utils.py:240
def check_histograms
Definition: utils.py:263

Member Data Documentation

utils.Chi2.name

Definition at line 238 of file utils.py.

Referenced by dataset.Dataset.__getFileInfoList(), cuy.divideElement.__init__(), cuy.plotElement.__init__(), cuy.additionElement.__init__(), cuy.superimposeElement.__init__(), cuy.graphElement.__init__(), and Vispa.Views.PropertyView.Property.valueChanged().