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

Constructor & Destructor Documentation

def utils.Chi2.__init__ (   self,
  threshold 
)

Definition at line 235 of file utils.py.

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

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().

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

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

Definition at line 262 of file utils.py.

Referenced by utils.Chi2.do_test().

262 
263  def check_histograms(self, histogram):
264  if histogram.InheritsFrom("TProfile") or (histogram.GetEntries()!=histogram.GetSumOfWeights()):
265  return 'W'
266  else:
267  return 'U'
utils.Chi2.check_histograms
def check_histograms
Definition: utils.py:262
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.

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

Member Data Documentation

utils.Chi2.name

Definition at line 237 of file utils.py.

Referenced by cuy.divideElement.__init__(), cuy.plotElement.__init__(), cuy.additionElement.__init__(), cuy.superimposeElement.__init__(), cuy.graphElement.__init__(), config.CFG.__str__(), validation.Sample.digest(), VIDSelectorBase.VIDSelectorBase.initialize(), and Vispa.Views.PropertyView.Property.valueChanged().

Generated for CMSSW Reference Manual by doxygen 1.8.5