Comp2Ref2DChi2 Class Reference

#include <QTest.h>

Inheritance diagram for Comp2Ref2DChi2:

Public Member Functions
	Comp2Ref2DChi2 (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)

Protected Member Functions
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
double	chi2_
int	Ndof_
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 212 of file QTest.h.

Constructor & Destructor Documentation

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

inline

Definition at line 215 of file QTest.h.

References QCriterion::setAlgoName().

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

Member Function Documentation

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

inlinestatic

Definition at line 219 of file QTest.h.

References QCriterion::runTest().

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

{ return "Comp2Ref2DChi2"; }

float Comp2Ref2DChi2::runTest ( const MonitorElement *  me )

virtual

Reimplemented from QCriterion.

Definition at line 292 of file QTest.cc.

References gather_cfg::cout, MonitorElement::DQM_KIND_TH2D, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TH2S, MonitorElement::DQM_KIND_TPROFILE2D, MonitorElement::getEntries(), MonitorElement::getRefRootObject(), MonitorElement::getRefTH2D(), MonitorElement::getRefTH2F(), MonitorElement::getRefTH2S(), MonitorElement::getRefTProfile2D(), MonitorElement::getRootObject(), MonitorElement::getTH2D(), MonitorElement::getTH2F(), MonitorElement::getTH2S(), MonitorElement::getTProfile2D(), and MonitorElement::kind().

{
  if (!me) 
    return -1;
  if (!me->getRootObject() || !me->getRefRootObject()) 
    return -1;
  if (minEntries_ != 0 && me->getEntries() < minEntries_)
    return -1;

  TH2* h=0;
  TH2* ref_=0;
 
  if (verbose_>1) 
    std::cout << "QTest:" << getAlgoName() << "::runTest called on " 
              << me-> getFullname() << "\n";
  //-- TH2F
  if (me->kind()==MonitorElement::DQM_KIND_TH2F)
  { 
    h = me->getTH2F(); // access Test histo
    ref_ = me->getRefTH2F(); //access Ref histo
  } 
  //-- TH2S
  else if (me->kind()==MonitorElement::DQM_KIND_TH2S)
  { 
    h = me->getTH2S(); // access Test histo
    ref_ = me->getRefTH2S(); //access Ref histo
  } 
  //-- TH2D
  else if (me->kind()==MonitorElement::DQM_KIND_TH2D)
  { 
    h = me->getTH2D(); // access Test histo
    ref_ = me->getRefTH2D(); //access Ref histo
  } 
  //-- TProfile
  else if (me->kind()==MonitorElement::DQM_KIND_TPROFILE2D)
  {
    h = me->getTProfile2D(); // access Test histo
    ref_ = me->getRefTProfile2D(); //access Ref histo
  } 
  else
  { 
    if (verbose_>0) 
      std::cout << "QTest::Comp2Ref2DChi2"
                << " ME does not contain TH2F/TH2S/TH2D/TProfile2D, exiting\n"; 
    return -1;
  } 

  //-- isInvalid ? - Check consistency in number of channels
  int ncx1  = h->GetXaxis()->GetNbins(); 
  int ncx2  = ref_->GetXaxis()->GetNbins();
  int ncy1  = h->GetYaxis()->GetNbins(); 
  int ncy2  = ref_->GetYaxis()->GetNbins();
  if ( ( ncx1 !=  ncx2) || ( ncy1 !=  ncy2) )
  {
    if (verbose_>0) 
      std::cout << "QTest:Comp2Ref2DChi2"
                << " different number of channels! ("
                << ncx1 << ", " << ncx2 << "), ("
                << ncy1 << ", " << ncy2 << "), exiting\n";
    return -1;
  } 

  //--  QUALITY TEST itself 
  //reset Results
  Ndof_ = 0; chi2_ = -1.;

  //check that the histograms are not empty
  int i_start = h->GetXaxis()->GetFirst();
  int i_end   = h->GetXaxis()->GetLast();
  int j_start = h->GetYaxis()->GetFirst();
  int j_end   = h->GetYaxis()->GetLast();
  if (h->Integral(i_start, i_end, j_start, j_end) == 0)
  {
    if (verbose_>0) 
      std::cout << "QTest:Comp2Ref2DChi2"
                << " Test Histogram " << h->GetName() 
        << " is empty, exiting\n";
    return -1;
  }
  if (ref_->Integral(i_start, i_end, j_start, j_end) == 0)
  {
    if (verbose_>0) 
      std::cout << "QTest:Comp2Ref2DChi2"
                << " Ref Histogram " << ref_->GetName() 
                << " is empty, exiting\n";
    return -1;
  }

  //use the chi2 test for 2D histograms defined in ROOT
  int igood = 0;
  double pValue = h->Chi2TestX(ref_, chi2_, Ndof_, igood, "");

  if (chi2_==-1. && Ndof_==0)
    return -1;

  return pValue;
}

void Comp2Ref2DChi2::setMessage ( void  )

inlineprotectedvirtual

set status & message after test has run

Reimplemented from SimpleTest.

Definition at line 224 of file QTest.h.

References QCriterion::errorProb_, python.rootplot.argparse::message, QCriterion::message_, and QCriterion::warningProb_.

  {
    std::ostringstream message;
    message << "chi2/Ndof = " << chi2_ << "/" << Ndof_
        << ", minimum needed statistics = " << minEntries_
        << " warning threshold = " << this->warningProb_
        << " error threshold = " << this->errorProb_;
    message_ = message.str();
  }

Member Data Documentation

double Comp2Ref2DChi2::chi2_

protected

Definition at line 235 of file QTest.h.

int Comp2Ref2DChi2::Ndof_

protected

Definition at line 235 of file QTest.h.