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

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

Protected Member Functions inherited from QCriterion
void	init ()
	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 ()=default

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

Constructor & Destructor Documentation

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

inline

Definition at line 628 of file QTest.h.

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

overridedefault

Member Function Documentation

static std::string CompareToMedian::getAlgoName ( )

inlinestatic

Definition at line 642 of file QTest.h.

References QCriterion::runTest().

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

642 { return "CompareToMedian"; }

void CompareToMedian::reset ( void )

inlineprivate

Definition at line 673 of file QTest.h.

673 {binValues.clear();};

CompareToMedian::binValues

std::vector< float > binValues

Definition: QTest.h:671

float CompareToMedian::runTest ( const MonitorElement * me )

overridevirtual

Reimplemented from QCriterion.

Definition at line 1375 of file QTest.cc.

References stringResolutionProvider_cfi::bin, officialStyle::chan, ws_sso_content_reader::content, gather_cfg::cout, MonitorElement::DQM_KIND_TPROFILE2D, runEdmFileComparison::failed, objects.autophobj::float, MonitorElement::getRootObject(), MonitorElement::getTProfile2D(), createfilelist::int, 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=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::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);
       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;
 }

void CompareToMedian::setEmptyBins ( int eB )

inline

Definition at line 647 of file QTest.h.

Referenced by QTestConfigure::EnableCompareToMedianTest().

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

CompareToMedian::_emptyBins

int _emptyBins

Definition: QTest.h:662

void CompareToMedian::setMax ( float max )

inline

Definition at line 646 of file QTest.h.

References SiStripPI::max.

Referenced by QTestConfigure::EnableCompareToMedianTest().

646 {_max = max;};

SiStripPI::max

Definition: SiStripPayloadInspectorHelper.h:19

CompareToMedian::_max

float _max

Definition: QTest.h:661

void CompareToMedian::setMaxMedian ( float max )

inline

Definition at line 648 of file QTest.h.

References SiStripPI::max.

Referenced by QTestConfigure::EnableCompareToMedianTest().

648 {_maxMed = max;};

CompareToMedian::_maxMed

float _maxMed

Definition: QTest.h:663

SiStripPI::max

Definition: SiStripPayloadInspectorHelper.h:19

void CompareToMedian::setMessage ( )

inlineoverrideprotectedvirtual

set status & message after test has run

Reimplemented from SimpleTest.

Definition at line 653 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();
   }