CheckVariance Class Reference

#include <QTest.h>

Inheritance diagram for CheckVariance:

Public Member Functions
	CheckVariance (const std::string &name)
float	runTest (const MonitorElement *me)
Public Member Functions inherited from SimpleTest
virtual std::vector< DQMChannel >	getBadChannels (void) const
	get vector of channels that failed test (not always relevant!) More...
void	setMinimumEntries (unsigned n)
	set minimum # of entries needed More...
	SimpleTest (const std::string &name, bool keepBadChannels=false)
Public Member Functions inherited from QCriterion
std::string	algoName (void) const
	get algorithm name More...
std::string	getMessage (void) const
	get message attached to test More...
std::string	getName (void) const
	get name of quality test More...
int	getStatus (void) const
	(class should be created by DQMStore class) More...
void	setErrorProb (float prob)
void	setWarningProb (float prob)
	set probability limit for warning and error (default: 90% and 50%) More...

Static Public Member Functions
static std::string	getAlgoName (void)
	get algorithm name More...

Additional Inherited Members
Protected Member Functions inherited from SimpleTest
virtual void	setMessage (void)
	set status & message after test has run More...
Protected Member Functions inherited from QCriterion
void	init (void)
	initialize values More...
	QCriterion (std::string qtname)
float	runTest (const MonitorElement *me, QReport &qr, DQMNet::QValue &qv)
void	setAlgoName (std::string name)
	set algorithm name More...
void	setVerbose (int verbose)
	probability limits for warnings, errors More...
virtual	~QCriterion (void)
Protected Attributes inherited from SimpleTest
std::vector< DQMChannel >	badChannels_
bool	keepBadChannels_
unsigned	minEntries_
Protected Attributes inherited from QCriterion
std::string	algoName_
	name of quality test More...
float	errorProb_
std::string	message_
	quality test status More...
float	prob_
	name of algorithm More...
std::string	qtname_
int	status_
int	verbose_
float	warningProb_
	message attached to test More...

Detailed Description

Definition at line 743 of file QTest.h.

Constructor & Destructor Documentation

CheckVariance::CheckVariance ( const std::string &  name )

inline

Definition at line 746 of file QTest.h.

References QCriterion::setAlgoName().

                                       : SimpleTest(name)
    {
      setAlgoName(getAlgoName());
    }

Member Function Documentation

static std::string CheckVariance::getAlgoName ( void  )

inlinestatic

get algorithm name

Definition at line 751 of file QTest.h.

References QCriterion::runTest().

Referenced by QTestConfigure::EnableCheckVarianceTest(), QTestConfigure::enableTests(), and QTestParameterNames::QTestParameterNames().

{ return "CheckVariance"; }

float CheckVariance::runTest ( const MonitorElement *  me )

virtual

Reimplemented from QCriterion.

Definition at line 1554 of file QTest.cc.

References stringResolutionProvider_cfi::bin, relmon_rootfiles_spy::contents, gather_cfg::cout, MonitorElement::DQM_KIND_TH1D, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TPROFILE, MonitorElement::getFullname(), MonitorElement::getRootObject(), MonitorElement::getTH1D(), MonitorElement::getTH1F(), MonitorElement::getTProfile(), MonitorElement::kind(), and combinedConstraintHelpers::sum2().

{
  badChannels_.clear();

  if (!me)
    return -1;
  if (!me->getRootObject())
    return -1;
  TH1* h=0;

  if (verbose_>1)
    std::cout << "QTest:" << getAlgoName() << "::runTest called on "
              << me-> getFullname() << "\n";
  // -- TH1F
  if (me->kind()==MonitorElement::DQM_KIND_TH1F )
    {
      h = me->getTH1F();
    }
  // -- TH1D
  else if ( me->kind()==MonitorElement::DQM_KIND_TH1D )
    {
      h = me->getTH1D();
    }
  else if (me->kind()==MonitorElement::DQM_KIND_TPROFILE)
    {
      h  = me->getTProfile(); // access Test histo
    }
  else
    {
      if (verbose_>0) std::cout << "QTest:CheckVariance"
                                << " ME " << me->getFullname()
                                << " does not contain TH1F/TH1D/TPROFILE, exiting\n";
      return -1;
    }

  int ncx = h->GetXaxis()->GetNbins();

  double sum = 0;
  double sum2 = 0;
  for (int bin = 1; bin <= ncx; ++bin)
    {
      double contents = h->GetBinContent(bin);
      sum += contents;
    }
  if (sum==0) return -1;
  double avg = sum/ncx;

  for (int bin = 1; bin <= ncx; ++bin)
    {
      double contents = h->GetBinContent(bin);
      sum2 += (contents - avg)*(contents - avg);
    }

  double Variance = TMath::Sqrt(sum2/ncx);
  return Variance;
}