#include <QTest.h>

Inheritance diagram for CompareToMedian:

Public Member Functions
	CompareToMedian (const std::string &name)

float	runTest (const MonitorElement *me)

void	setEmptyBins (int eB)

void	setMax (float max)

void	setMaxMedian (float max)

void	setMin (float min)

void	setMinMedian (float min)

void	setStatCut (float cut)

	~CompareToMedian ()

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)

Private Member Functions
void	reset ()

Private Attributes
int	_emptyBins

float	_max

float	_maxMed

float	_min

float	_minMed

float	_statCut

std::vector< float >	binValues

int	nBins

int	nBinsX

int	nBinsY

Additional Inherited Members
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 593 of file QTest.h.

Constructor & Destructor Documentation

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

inline

Definition at line 597 of file QTest.h.

References _emptyBins, _max, _maxMed, _min, _minMed, _statCut, getAlgoName(), nBins, reset(), and QCriterion::setAlgoName().

                                          : SimpleTest(name,true){
     this->_min = 0.2;
     this->_max = 3.0;
     this->_emptyBins = 0;
     this->_maxMed = 10;
     this->_minMed = 0;
     this->nBins = 0;
     this->_statCut = 0;
     reset();
     setAlgoName( getAlgoName() );
   };

CompareToMedian::~CompareToMedian ( )

inline

Definition at line 609 of file QTest.h.

609 {};

Member Function Documentation

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

inlinestatic

Definition at line 611 of file QTest.h.

Referenced by CompareToMedian(), QTestConfigure::EnableCompareToMedianTest(), QTestConfigure::enableTests(), and QTestParameterNames::QTestParameterNames().

611 { return "CompareToMedian"; }

void CompareToMedian::reset ( void )

inlineprivate

Definition at line 642 of file QTest.h.

References binValues.

Referenced by CompareToMedian().

642 {binValues.clear();};

CompareToMedian::binValues

std::vector< float > binValues

Definition: QTest.h:640

float CompareToMedian::runTest ( const MonitorElement * me )

virtual

Reimplemented from QCriterion.

Definition at line 1264 of file QTest.cc.

References newFWLiteAna::bin, ws_sso_content_reader::content, gather_cfg::cout, MonitorElement::DQM_KIND_TPROFILE2D, python.tagInventory::entries, MonitorElement::getRootObject(), MonitorElement::getTProfile2D(), h, MonitorElement::kind(), pileupCalc::nbins, and reset().

                                                       {
   int32_t nbins=0, failed=0;
   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";
     std::cout << "\tMin = " << _min << "; Max = " << _max << "\n";
     std::cout << "\tMinMedian = " << _minMed << "; MaxMedian = " << _maxMed << "\n";
     std::cout << "\tUseEmptyBins = " << (_emptyBins?"Yes":"No" )<< "\n";
   }
 
   if (me->kind()==MonitorElement::DQM_KIND_TPROFILE2D)
     {
       h  = me->getTProfile2D(); // access Test histo
     }
   else
     {
       if (verbose_>0) 
         std::cout << "QTest:ContentsWithinExpected" 
      << " ME does not contain TPROFILE2D, exiting\n"; 
       return -1;
     }
   
   nBinsX = h->GetNbinsX();
   nBinsY = h->GetNbinsY();
   int entries = 0;
   float median = 0.0;
 
   //Median calculated with partially sorted vector
   for (int binX = 1; binX <= nBinsX; binX++ ){
     reset();
     // Fill vector
     for (int binY = 1; binY <= nBinsY; binY++){
       int bin = h->GetBin(binX, binY);
       double content = (double)h->GetBinContent(bin);
       if ( content == 0 && !_emptyBins)
     continue;
       binValues.push_back(content);
       entries = me->getTProfile2D()->GetBinEntries(bin);
     }
     if (binValues.empty())
       continue;
     nbins+=binValues.size();
 
     //calculate median
     if(binValues.size()>0){
       int medPos = (int)binValues.size()/2;
       nth_element(binValues.begin(),binValues.begin()+medPos,binValues.end());
       median = *(binValues.begin()+medPos);
     }
 
     // if median == 0, use the absolute cut
     if(median == 0){
       if (verbose_ > 0){
         std::cout << "QTest: Median is 0; the fixed cuts: [" << _minMed << "; " << _maxMed << "]  are used\n";
       }
       for(int  binY = 1; binY <= nBinsY; binY++){
         int bin = h->GetBin(binX, binY);
         double content = (double)h->GetBinContent(bin);
     entries = me->getTProfile2D()->GetBinEntries(bin);
     if ( entries == 0 )
           continue;
     if (content > _maxMed || content < _minMed){ 
         DQMChannel chan(binX,binY, 0, content, h->GetBinError(bin));
         badChannels_.push_back(chan);
         failed++;
     }
       }
       continue;
     }
 
     //Cut on stat: will mask rings with no enought of statistics
     if (median*entries < _statCut )
       continue;
      
     // If median is off the absolute cuts, declare everything bad (if bin has non zero entries)
     if(median > _maxMed || median < _minMed){
       for(int  binY = 1; binY <= nBinsY; binY++){
         int bin = h->GetBin(binX, binY);
         double content = (double)h->GetBinContent(bin);
     entries = me->getTProfile2D()->GetBinEntries(bin);
          if ( entries == 0 )
           continue;
      DQMChannel chan(binX,binY, 0, content/median, h->GetBinError(bin));
      badChannels_.push_back(chan);
      failed++;
       }
       continue;
     }
 
     // Test itself  
     float minCut = median*_min;
     float maxCut = median*_max;
     for(int  binY = 1; binY <= nBinsY; binY++){
         int bin = h->GetBin(binX, binY);
         double content = (double)h->GetBinContent(bin);
     entries = me->getTProfile2D()->GetBinEntries(bin);
         if ( entries == 0 )
           continue;
         if (content > maxCut || content < minCut){
           DQMChannel chan(binX,binY, 0, content/median, h->GetBinError(bin));
           badChannels_.push_back(chan);
           failed++;
         }
     }
   }
 
   if (nbins==0){
     if (verbose_ > 0)
       std::cout << "QTest:CompareToMedian: Histogram is empty" << std::endl;
     return 1.;
     }
   return 1 - (float)failed/nbins;
 }

void CompareToMedian::setEmptyBins ( int eB )

inline

Definition at line 616 of file QTest.h.

References _emptyBins.

Referenced by QTestConfigure::EnableCompareToMedianTest().

616 {eB > 0 ? _emptyBins = 1 : _emptyBins = 0;};

CompareToMedian::_emptyBins

int _emptyBins

Definition: QTest.h:631

void CompareToMedian::setMax ( float max )

inline

Definition at line 615 of file QTest.h.

References _max, and max().

Referenced by QTestConfigure::EnableCompareToMedianTest().

615 {_max = max;};

max

const T & max(const T &a, const T &b)

Definition: MaterialBudgetTrackerHistos.cc:4

CompareToMedian::_max

float _max

Definition: QTest.h:630

void CompareToMedian::setMaxMedian ( float max )

inline

Definition at line 617 of file QTest.h.

References _maxMed, and max().

Referenced by QTestConfigure::EnableCompareToMedianTest().

617 {_maxMed = max;};

max

const T & max(const T &a, const T &b)

Definition: MaterialBudgetTrackerHistos.cc:4

CompareToMedian::_maxMed

float _maxMed

Definition: QTest.h:632

void CompareToMedian::setMessage ( void )

inlineprotectedvirtual

set status & message after test has run

Reimplemented from SimpleTest.

Definition at line 622 of file QTest.h.

References QCriterion::algoName_, python.rootplot.argparse::message, QCriterion::message_, QCriterion::prob_, and QCriterion::qtname_.

                        {
     std::ostringstream message;
     message << "Test " << qtname_ << " (" << algoName_
             << "): Entry fraction within range = " << prob_;
     message_ = message.str();
   }