dqm_algorithms::MDTTubeCheck Class Reference

#include <MDTTubeCheck.h>

Inheritance diagram for dqm_algorithms::MDTTubeCheck:

Collaboration diagram for dqm_algorithms::MDTTubeCheck:

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

Private Attributes
std::string	m_name

Detailed Description

Definition at line 20 of file MDTTubeCheck.h.

Constructor & Destructor Documentation

MDTTubeCheck()

dqm_algorithms::MDTTubeCheck::MDTTubeCheck

(

)

Definition at line 37 of file MDTTubeCheck.cxx.

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

~MDTTubeCheck()

dqm_algorithms::MDTTubeCheck::~MDTTubeCheck ( )

virtual

Definition at line 44 of file MDTTubeCheck.cxx.

{
}

Member Function Documentation

clone()

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

virtual

Definition at line 50 of file MDTTubeCheck.cxx.

{
  return new MDTTubeCheck(*this);
}

execute()

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

virtual

Definition at line 57 of file MDTTubeCheck.cxx.

{
  const TH1 * histogram(0);
  const TH1 * refhist(0);
  
  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" );
  }   
  
  //Get Parameters and Thresholds
  double minstat; 
  double bin_threshold;
  double refcheck;
  double greenTh;
  double redTh;
  double LowStatLevel;
  try {
    bin_threshold = dqm_algorithms::tools::GetFirstFromMap( "BinThreshold", config.getParameters());
    refcheck      = dqm_algorithms::tools::GetFirstFromMap( "ReferenceCheck", config.getParameters());
    LowStatLevel  = dqm_algorithms::tools::GetFirstFromMap( "LowStatLevel", config.getParameters(), 10);
    minstat       = dqm_algorithms::tools::GetFirstFromMap( "MinStat", config.getParameters(), -1);
    redTh         = dqm_algorithms::tools::GetFromMap( "NBins", config.getRedThresholds());
    greenTh       = dqm_algorithms::tools::GetFromMap( "NBins", config.getGreenThresholds() );
  }
  catch ( dqm_core::Exception & ex ) {
    throw dqm_core::BadConfig( ERS_HERE, name, ex.what(), ex );
  }
 
  //Get Reference Histo
  if (refcheck>0) {
    try {
      refhist = static_cast<const TH1*>( config.getReference() );
    }
    catch ( dqm_core::Exception & ex ) {
      throw dqm_core::BadRefHist(ERS_HERE,name," Could not retreive reference");
    }
    if (histogram->GetDimension() != refhist->GetDimension() ) {
      throw dqm_core::BadRefHist( ERS_HERE, name, "Reference VS histo: Different dimension!" );
    } 
    if (histogram->GetNbinsX() != refhist->GetNbinsX() ) {
      throw dqm_core::BadRefHist( ERS_HERE, name, "Reference VS histo: Different number of bins!" );
    }
  }
    
 //Check of statistics
  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;
  }
  ERS_DEBUG(1,"Statistics: "<< histogram->GetEntries()<< " entries ");
 
  int NewDead_count = 0;
  int StillDead_count = 0;
  int Revived_count = 0;
  int LowStat_count = 0;
  std::vector<int> range;
  if (histogram->GetDimension() == 1){
    range = dqm_algorithms::tools::GetBinRange(histogram, config.getParameters()); 
  } else {
    ERS_INFO("Histogram has not dimension 1 : " <<histogram->GetName());
    throw dqm_core::Exception( ERS_HERE, histogram->GetName() );
  }
 
  std::vector<int> Tubes; 
  for ( int i = range[0]; i <= range[1]; ++i ) {
    double content = histogram->GetBinContent(i);
    if (content == LowStatLevel) {
      ++LowStat_count;
    } else if (content != 0.) {
      if (refcheck>0) {
    double RefCont = refhist->GetBinContent(i);
    if ((content < bin_threshold) && ((RefCont >= bin_threshold) || (RefCont == LowStatLevel))) {
      ++NewDead_count;
      Tubes.push_back(i);
    } else if ((content < bin_threshold) && ((RefCont < bin_threshold) || (RefCont == LowStatLevel))) {
      ++StillDead_count;
    } else if ((content >= bin_threshold) && (RefCont < bin_threshold)) {
      ++Revived_count;
    }
      } else {
    if (((content < bin_threshold) ) ) {
      ++NewDead_count;
      Tubes.push_back(i);
    }
      }
    }
  }
  ERS_DEBUG(1,"Number of bins " << name << " different from a treshold of " << bin_threshold << " is " << NewDead_count );
 
  dqm_core::Result* result = new dqm_core::Result(); 
  if (LowStat_count > 0) result->tags_["NTubes_withLowStat"] = LowStat_count;
  if (Revived_count > 0) result->tags_["NTubes_Revived"] = Revived_count;
  result->tags_["NewDead_Tubes"] = NewDead_count;
  if (Tubes.size()>0) {
    for (int k=0; k<(int)Tubes.size(); k++) {
      std::string ToDB="ChangedStatusTube_";
      ToDB += std::to_string(k+1);
      result->tags_[ToDB] = Tubes[k];
      ERS_DEBUG(1,"MDT Tube which changed status: "<<ToDB<<" = "<<Tubes[k] );
    }
  }
  if (NewDead_count >= redTh) {
    result->status_ = dqm_core::Result::Red; 
    ERS_DEBUG(1,"Red");
  } else if ( (NewDead_count > greenTh) || (Revived_count > 0)) {
    result->status_ = dqm_core::Result::Yellow;
    ERS_DEBUG(1,"Yellow");
  } else if (LowStat_count > 0) {
    result->status_ = dqm_core::Result::Undefined;
    ERS_DEBUG(1,"Undefined");
  } else if (StillDead_count == range[1]) {
    //All tubes are still dead
    result->status_ = dqm_core::Result::Undefined;
  } else {
    result->status_ = dqm_core::Result::Green;
    ERS_DEBUG(1,"Green");
  }
  return result;
}

printDescription()

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

virtual

Definition at line 186 of file MDTTubeCheck.cxx.

{
  std::string message;
  message += "\n";
  message += "Algorithm: \"" + m_name + "\"\n";
  message += "Description: Check if the number of entries in bins is less than BinThreshold and compare with Reference Histo \n";
  message += "Mandatory Parameter: BinThreshold = Look for bins less than BinTreshold; Count number of bins satifying requirement \n";
  message += "Mandatory Parameter: ReferenceCheck = 0 if no check on reference is requested \n";
  message += "Mandatory Green/Red Threshold: NBins = number of bins satifying requirement\n";
  message += "Optional Parameters: LowStatLevel = Level (in y) that means too few statistics for analysis\n";
  message += "Optional Parameters: MinStat = Minimum histogram statistics needed to perform Algorithm\n";
  message += "                     xmin: minimum x range\n";
  message += "                     xmax: maximum x range\n";
  message += "\n";
  
  out << message;
}

Member Data Documentation

m_name

std::string dqm_algorithms::MDTTubeCheck::m_name

private

Definition at line 33 of file MDTTubeCheck.h.

---

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

• MDTTubeCheck.h
• MDTTubeCheck.cxx