dqm_algorithms::HLTMETComponents Struct Reference

#include <HLTMETComponents.h>

Inheritance diagram for dqm_algorithms::HLTMETComponents:

Collaboration diagram for dqm_algorithms::HLTMETComponents:

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

Detailed Description

Definition at line 20 of file HLTMETComponents.h.

Constructor & Destructor Documentation

HLTMETComponents()

dqm_algorithms::HLTMETComponents::HLTMETComponents

(

)

Definition at line 27 of file HLTMETComponents.cxx.

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

Member Function Documentation

clone()

dqm_algorithms::HLTMETComponents * dqm_algorithms::HLTMETComponents::clone

(

)

Definition at line 34 of file HLTMETComponents.cxx.

{
  
  return new HLTMETComponents();
}

execute()

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

Definition at line 42 of file HLTMETComponents.cxx.

{  
  const TH2 *histogram;
  
  if( object.IsA()->InheritsFrom( "TH2" ) ) {
    histogram = static_cast<const TH2*>(&object);
    int dimension = histogram->GetDimension();
    if (dimension != 2 ){ 
      throw dqm_core::BadConfig( ERS_HERE, name, "dimension != 2 " );
    }
  } else {
    throw dqm_core::BadConfig( ERS_HERE, name, "does not inherit from TH2" );
  }
  
  int nXBins = histogram -> GetNbinsX(); // components
  if(nXBins > 25) {
    nXBins = 25;
  }
  const int nYBins = histogram -> GetNbinsY(); // SumEt per sampling
 
 
  // read thresholds
  int colflags[4];
  colflags[0] = static_cast<int>(dqm_algorithms::tools::GetFirstFromMap("METYellowMin", config.getParameters(), 1));
  colflags[1] = static_cast<int>(dqm_algorithms::tools::GetFirstFromMap("METYellowMax", config.getParameters(), nXBins));
  colflags[2] = static_cast<int>(dqm_algorithms::tools::GetFirstFromMap("METRedMin",    config.getParameters(), 1));
  colflags[3] = static_cast<int>(dqm_algorithms::tools::GetFirstFromMap("METRedMax",    config.getParameters(), nXBins));
  // read options
  bool doflags[2];
  doflags[0] = static_cast<bool>(dqm_algorithms::tools::GetFirstFromMap("DoYellowFlag", config.getParameters(), 1));
  doflags[1] = static_cast<bool>(dqm_algorithms::tools::GetFirstFromMap("DoRedFlag",    config.getParameters(), 1));
 
  dqm_core::Result* result = new dqm_core::Result();
  colflags[0] = ((colflags[0] > 0) && (colflags[0] <= nXBins)) ? colflags[0] : 1;
  colflags[1] = ((colflags[1] >= colflags[0]) && (colflags[1] <= nXBins)) ? colflags[1] : nXBins;
  colflags[2] = ((colflags[2] > 0) && (colflags[2] <= nXBins)) ? colflags[2] : 1;
  colflags[3] = ((colflags[3] >= colflags[2]) && (colflags[3] <= nXBins)) ? colflags[3] : nXBins;
  bool yellowFlagComp = false, redFlagComp = false;
  int ycnt = 0, rcnt = 0;
  
  // yellow flag
  float epsilon = 1.e-3;
  if(doflags[0]) {
    for( int j = colflags[0]; j <= colflags[1]; j++ ) {
      TString binLabel = histogram->GetXaxis()->GetBinLabel(j);
      TH1D *hprojy = histogram->ProjectionY(Form("hpy1yf%02d_%s",j,binLabel.Data()),j,j+1,"egoff");
      if(hprojy) {
        float uflow = std::abs(hprojy->GetBinContent(0));
        float oflow = std::abs(hprojy->GetBinContent(nYBins+1));
        float btotal = hprojy->Integral();
 
        // flag true if overflow/underflow has entries or if histo is empty
        // denoting missing component
        if(((uflow > epsilon) || (oflow > epsilon)) && (btotal < epsilon)) {
          yellowFlagComp = true;
          // count bins flagged as yellow
          ycnt++;
        } // if uflow
        delete hprojy; hprojy = 0;
      } // if hprojy
    } // for
  } // if doflags[0]
 
  // red flag
  if(doflags[1]) {
    for( int j = colflags[2]; j <= colflags[3]; j++ ) {
      TString binLabel = histogram->GetXaxis()->GetBinLabel(j);
      TH1D *hprojy = histogram->ProjectionY(Form("hpy1rf%02d_%s",j,binLabel.Data()),j,j+1,"egoff");
      if(hprojy) {
        float uflow = std::abs(hprojy->GetBinContent(0));
        float oflow = std::abs(hprojy->GetBinContent(nYBins+1));
        float btotal = hprojy->Integral();
        // flag true if overflow/underflow has entries or if histo is empty
        // denoting missing component
        if(((uflow > epsilon) || (oflow > epsilon)) && (btotal < epsilon)) {
          redFlagComp = true;
          rcnt++;
        } // if uflow
        delete hprojy; hprojy = 0;
      } // if hprojy
    } // for
  } // if doflags[1]
 
  // set results (default green, 
  result->status_ = dqm_core::Result::Green;
  if(doflags[0] && yellowFlagComp) {
    result->status_ = dqm_core::Result::Yellow;
  }
  if(doflags[1] && redFlagComp) {
    result->status_ = dqm_core::Result::Red;
  }
  result->tags_["YellowMETComponents"] = ycnt; 
  result->tags_["RedMETComponents"] = rcnt;
 
  
  return result;
}

printDescription()

void dqm_algorithms::HLTMETComponents::printDescription

(

std::ostream &

out

)

Definition at line 145 of file HLTMETComponents.cxx.

{
  
  out << "HLT MET Components: Of the 25 components (24 calo samplings and 1 muon component of MET) ";
  out << "Check if any sub-detector component is empty. Empty means any of the following:";
  out << " overflow/underflow bins are filled, or integral of histogram is zero. " << std::endl;
  out<<  "Each component (25 of them) are bins on abscissa and the 1d histograms represent the quantity (SumEt)" << std::endl;
  
}

---

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

• HLTMETComponents.h
• HLTMETComponents.cxx