#include <QTest.h>

Inheritance diagram for CompareToMedian:

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

float	runTest (const MonitorElement *me) override

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 () override=default

Public Member Functions inherited from SimpleTest
std::vector< DQMChannel >	getBadChannels () const override
	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 () const
	get algorithm name More...

std::string	getMessage () const
	get message attached to test More...

std::string	getName () const
	get name of quality test More...

int	getStatus () const
	get test status More...

void	init ()
	initialize values More...

	QCriterion (std::string qtname)

float	runTest (const MonitorElement *me, QReport &qr, DQMNet::QValue &qv)

void	setErrorProb (float prob)

void	setWarningProb (float prob)
	set probability limit for warning and error (default: 90% and 50%) More...

virtual	~QCriterion ()=default

Static Public Member Functions
static std::string	getAlgoName ()

Protected Member Functions
void	setMessage () override
	set status & message after test has run More...

Protected Member Functions inherited from QCriterion
void	setAlgoName (std::string name)
	set algorithm name More...

void	setVerbose (int verbose)
	probability limits for warnings, errors More...

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
Public Types inherited from QCriterion
typedef dqm::legacy::MonitorElement	MonitorElement
	(class should be created by DQMStore class) More...

Static Public Attributes inherited from QCriterion
static const float	ERROR_PROB_THRESHOLD = 0.50

static const float	WARNING_PROB_THRESHOLD = 0.90
	default "probability" values for setting warnings & errors when running tests More...

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 589 of file QTest.h.

Constructor & Destructor Documentation

◆ CompareToMedian()

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

inline

Definition at line 592 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::~CompareToMedian ( )

overridedefault

Member Function Documentation

◆ getAlgoName()

static std::string CompareToMedian::getAlgoName ( )

inlinestatic

Definition at line 607 of file QTest.h.

Referenced by CompareToMedian(), and QualityTester::makeQCriterion().

607 { return "CompareToMedian"; }

◆ reset()

void CompareToMedian::reset ( void )

inlineprivate

Definition at line 637 of file QTest.h.

References binValues.

Referenced by CompareToMedian().

637 { binValues.clear(); };

CompareToMedian::binValues

std::vector< float > binValues

Definition: QTest.h:635

◆ runTest()

float CompareToMedian::runTest ( const MonitorElement * me )

overridevirtual

Reimplemented from QCriterion.

Definition at line 1177 of file QTest.cc.

References newFWLiteAna::bin, officialStyle::chan, submitPVResolutionJobs::content, gather_cfg::cout, runEdmFileComparison::failed, ALCARECOEcalPhiSym_cff::float, h, createfilelist::int, hlt_dqm_clientPB-live_cfg::me, median(), LaserClient_cfi::nbins, reset(), and MonitorElementData::TPROFILE2D.

                                                        {
   int32_t nbins = 0, failed = 0;
   badChannels_.clear();
 
   if (!me)
     return -1;
   if (!me->getRootObject())
     return -1;
   TH1* h = nullptr;
 
   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::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);
       auto 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.empty()) {
       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);
         auto 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);
         auto 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);
       auto 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;
 }

◆ setEmptyBins()

void CompareToMedian::setEmptyBins ( int eB )

inline

Definition at line 612 of file QTest.h.

References _emptyBins.

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

CompareToMedian::_emptyBins

int _emptyBins

Definition: QTest.h:626

◆ setMax()

void CompareToMedian::setMax ( float max )

inline

Definition at line 611 of file QTest.h.

References _max, and WZElectronSkims53X_cff::max.

611 { _max = max; };

WZElectronSkims53X_cff.max

max

Definition: WZElectronSkims53X_cff.py:197

CompareToMedian::_max

float _max

Definition: QTest.h:625

◆ setMaxMedian()

void CompareToMedian::setMaxMedian ( float max )

inline

Definition at line 613 of file QTest.h.

References _maxMed, and WZElectronSkims53X_cff::max.

613 { _maxMed = max; };

WZElectronSkims53X_cff.max

max

Definition: WZElectronSkims53X_cff.py:197

CompareToMedian::_maxMed

float _maxMed

Definition: QTest.h:627

◆ setMessage()

void CompareToMedian::setMessage ( )

inlineoverrideprotectedvirtual

set status & message after test has run

Reimplemented from SimpleTest.

Definition at line 618 of file QTest.h.

References QCriterion::algoName_, QCriterion::message_, QCriterion::prob_, and QCriterion::qtname_.

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