dqm_algorithms::TRTCheckPeakSimple Class Reference

#include <TRTCheckPeakSimple.h>

Inheritance diagram for dqm_algorithms::TRTCheckPeakSimple:

Collaboration diagram for dqm_algorithms::TRTCheckPeakSimple:

Public Member Functions
	TRTCheckPeakSimple ()
virtual	~TRTCheckPeakSimple ()
virtual dqm_core::Algorithm *	clone ()
virtual dqm_core::Result *	execute (const std::string &, const TObject &, const dqm_core::AlgorithmConfig &)
virtual void	printDescription (std::ostream &out)

Private Attributes
std::string	m_name

Detailed Description

Definition at line 20 of file TRTCheckPeakSimple.h.

Constructor & Destructor Documentation

TRTCheckPeakSimple()

dqm_algorithms::TRTCheckPeakSimple::TRTCheckPeakSimple

(

)

Definition at line 30 of file TRTCheckPeakSimple.cxx.

                                      : m_name("TRTCheckPeakSimple")
{
    dqm_core::AlgorithmManager::instance().registerAlgorithm(m_name, this);
}

~TRTCheckPeakSimple()

dqm_algorithms::TRTCheckPeakSimple::~TRTCheckPeakSimple ( )

virtual

Definition at line 35 of file TRTCheckPeakSimple.cxx.

{
}

Member Function Documentation

clone()

dqm_core::Algorithm * dqm_algorithms::TRTCheckPeakSimple::clone ( )

virtual

Definition at line 39 of file TRTCheckPeakSimple.cxx.

{
    return new TRTCheckPeakSimple(*this);
}

execute()

dqm_core::Result * dqm_algorithms::TRTCheckPeakSimple::execute ( const std::string &  name, const TObject &  object, const dqm_core::AlgorithmConfig &  config  )

virtual

Definition at line 45 of file TRTCheckPeakSimple.cxx.

{
    const TH1 *histogram;
 
    if (object.IsA()->InheritsFrom("TH1")) {
        histogram = static_cast<const TH1*>(&object);
        if (histogram->GetDimension() > 2) {
            throw dqm_core::BadConfig(ERS_HERE, name, "dimension > 2 ");
        }
    } else {
        throw dqm_core::BadConfig(ERS_HERE, name, "does not inherit from TH1");
    }
    const double minstat = dqm_algorithms::tools::GetFirstFromMap("MinStat", config.getParameters(), -1);
 
    if (histogram->GetEntries() < minstat) {
        // ERS_INFO("Histogram does not satisfy MinStat requirement " << histogram->GetName());
        dqm_core::Result *result = new dqm_core::Result(dqm_core::Result::Undefined);
        result->tags_["InsufficientEntries"] = histogram->GetEntries();
        return result;
    }
 
    dqm_core::Result *result = new dqm_core::Result();
 
    if (histogram->Integral() == 0 && histogram->GetEntries() > 0) {
        ERS_DEBUG(1, "Histogram " << histogram->GetName() << " is filled with zeroes!");
        result->status_ = dqm_core::Result::Red;
        result->tags_["Integral"] = histogram->Integral();
        return result;
    }
 
 
    const double gthreshold = dqm_algorithms::tools::GetFromMap("PeakPosition", config.getGreenThresholds());
    const double rthreshold = dqm_algorithms::tools::GetFromMap("PeakPosition", config.getRedThresholds());
    const std::vector<int> range = dqm_algorithms::tools::GetBinRange(histogram, config.getParameters());
    const int xmin = range[0];
    const int xmax = range[1];
 
    //compute the weighted mean
    float wmean = 0;
    float sum = 0;
    float integral = 0;
    for (int i = xmin; i <= xmax; i++) {
        integral += histogram->GetBinContent(i);
        sum += i * (histogram->GetBinContent(i));
    }
    wmean = sum / integral; // make sure that integer division does not truncate decimal of mean
 
    //find the most probable value (peak)
    float maxbinsum = 0;
    float binsum = 0;
    float peakbin = 0;
    float peakpos = 0;
 
    for (int i = xmin + 1; i <= xmax - 1; i++) {
        binsum = histogram->GetBinContent(i - 1) + histogram->GetBinContent(i) + histogram->GetBinContent(i + 1);
        if (binsum > maxbinsum) {
            maxbinsum = binsum;
            if (maxbinsum != 0) {
                peakbin = ((i - 1) * histogram->GetBinContent(i - 1) + i * histogram->GetBinContent(i) + (i + 1) * histogram->GetBinContent(i + 1)) /  maxbinsum;
                peakpos = (histogram->GetBinCenter(i - 1) * histogram->GetBinContent(i - 1) + histogram->GetBinCenter(i) * histogram->GetBinContent(i) + histogram->GetBinCenter(i + 1) * histogram->GetBinContent(i + 1)) / maxbinsum;
            } else {
                peakbin = 0;
                peakpos = 0;
            }
        }
    }
 
    result->tags_["Weighted_mean"] = wmean;
    result->tags_["PeakBin"] = peakbin;
    result->tags_["PeakPosition"] = peakpos;
 
    if      (peakpos >= gthreshold) result->status_ = dqm_core::Result::Green;
    else if (peakpos >  rthreshold) result->status_ = dqm_core::Result::Yellow;
    else                            result->status_ = dqm_core::Result::Red;
    return result;
}

printDescription()

void dqm_algorithms::TRTCheckPeakSimple::printDescription ( std::ostream &  out )

virtual

Definition at line 122 of file TRTCheckPeakSimple.cxx.

{
    out << m_name << ": Checks on the most probable value (peak) in the histogram." << std::endl;
}

Member Data Documentation

m_name

std::string dqm_algorithms::TRTCheckPeakSimple::m_name

private

Definition at line 32 of file TRTCheckPeakSimple.h.

---

The documentation for this class was generated from the following files:

• TRTCheckPeakSimple.h
• TRTCheckPeakSimple.cxx