SiStripSummaryCreator Class Reference

#include <SiStripSummaryCreator.h>

Public Member Functions
void	createSummary (DQMStore &dqm_store)
void	fillLayout (DQMStore &dqm_store)
int	getFrequency () const
bool	readConfiguration (std::string const &file_path)
void	setSummaryMENames (std::map< std::string, std::string > &me_names)
	SiStripSummaryCreator ()
virtual	~SiStripSummaryCreator ()

Private Member Functions
void	fillGrandSummaryHistos (DQMStore &dqm_store)
void	fillHistos (int ival, int istep, std::string htype, MonitorElement *me_src, MonitorElement *me)
void	fillSummaryHistos (DQMStore &dqm_store)
MonitorElement *	getSummaryME (DQMStore &dqm_store, std::string &name, std::string htype)

Private Attributes
int	summaryFrequency_ {-1}
std::map< std::string, std::string >	summaryMEs_

Detailed Description

Definition at line 14 of file SiStripSummaryCreator.h.

Constructor & Destructor Documentation

SiStripSummaryCreator::SiStripSummaryCreator

(

)

Definition at line 9 of file SiStripSummaryCreator.cc.

{
  edm::LogInfo("SiStripSummaryCreator") << " Creating SiStripSummaryCreator "
                                        << "\n";
}

SiStripSummaryCreator::~SiStripSummaryCreator ( )

virtual

Definition at line 15 of file SiStripSummaryCreator.cc.

                                              {
  edm::LogInfo("SiStripSummaryCreator") << " Deleting SiStripSummaryCreator "
                                        << "\n";
}

Member Function Documentation

void SiStripSummaryCreator::createSummary

(

DQMStore &

dqm_store

)

Definition at line 56 of file SiStripSummaryCreator.cc.

References DQMStore::cd(), fillGrandSummaryHistos(), fillSummaryHistos(), DQMStore::getSubdirs(), DQMStore::goUp(), DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, and summaryMEs_.

{
  if (summaryMEs_.empty()) return;
  std::string currDir = dqm_store.pwd();
  std::vector<std::string> subdirs = dqm_store.getSubdirs();
  int nmod = 0;
  for (std::vector<std::string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
    if ( (*it).find("module_") == std::string::npos) continue;
    nmod++;
  }
  if (nmod > 0) {
    fillSummaryHistos(dqm_store);
  } else {
    for (std::vector<std::string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
      dqm_store.cd(*it);
      createSummary(dqm_store);
      dqm_store.goUp();
    }
    fillGrandSummaryHistos(dqm_store);
  }
}

void SiStripSummaryCreator::fillGrandSummaryHistos ( DQMStore &  dqm_store )

private

Definition at line 135 of file SiStripSummaryCreator.cc.

References DQMStore::cd(), relmon_rootfiles_spy::contents, compare_using_db::dir_name, fillHistos(), DQMStore::getContents(), MonitorElement::getName(), MonitorElement::getNbinsX(), DQMStore::getSubdirs(), getSummaryME(), DQMStore::goUp(), isum, dataset::name, DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, and summaryMEs_.

Referenced by createSummary(), and getFrequency().

{
  std::map<std::string, MonitorElement*> MEMap;
  std::string currDir = dqm_store.pwd();
  std::string dir_name =  currDir.substr(currDir.find_last_of("/")+1);
  if ((dir_name.find("SiStrip") == 0) ||
      (dir_name.find("Collector") == 0) ||
      (dir_name.find("MechanicalView") == 0) ||
      (dir_name.find("FU") == 0) ) return;
  std::vector<std::string> subdirs = dqm_store.getSubdirs();
  if (subdirs.empty()) return;;
  for (std::map<std::string, std::string>::const_iterator isum =
         summaryMEs_.begin();
       isum != summaryMEs_.end();
       isum++) {
    std::string name, summary_name;
    name = isum->first;
    if (isum->second == "sum" || isum->second == "sum")
      summary_name = "Summary_" + isum->first;
    else
      summary_name = "Summary_Mean" + isum->first;
    std::string htype = isum->second;
    int ibinStep =0;
    for (std::vector<std::string>::const_iterator it = subdirs.begin();
         it != subdirs.end();
         it++) {
      dqm_store.cd(*it);
      std::vector<MonitorElement*> contents =
        dqm_store.getContents(dqm_store.pwd());
      dqm_store.goUp();
      for (std::vector<MonitorElement *>::const_iterator im = contents.begin();
                im != contents.end(); im++) {
        MonitorElement * me_i = (*im);
        if (!me_i) continue;
        std::string name_i = me_i->getName();
        if (name_i.find((summary_name)) != std::string::npos) {

          std::map<std::string, MonitorElement*>::iterator iPos =
            MEMap.find(name);
          MonitorElement* me;
          if (iPos == MEMap.end()) {
            if (htype == "sum" || htype == "Sum") {
              me = getSummaryME(dqm_store, name, htype);
            } else {
              me = getSummaryME(dqm_store, name, "bin-by-bin");
            }
            MEMap.insert(std::pair<std::string, MonitorElement*>(name, me));
          } else  me =  iPos->second;
          if (htype == "sum" || htype == "Sum") {
            fillHistos(0, ibinStep, htype, me_i, me);
          } else {
            fillHistos(0, ibinStep, "bin-by-bin", me_i, me);
          }
          ibinStep += me_i->getNbinsX();
          break;
        }
      }
    }
  }
}

void SiStripSummaryCreator::fillHistos ( int  ival, int  istep, std::string  htype, MonitorElement *  me_src, MonitorElement *  me  )

private

Definition at line 308 of file SiStripSummaryCreator.cc.

References MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2F, MonitorElement::Fill(), MonitorElement::getBinContent(), MonitorElement::getEntries(), MonitorElement::getMean(), MonitorElement::getName(), MonitorElement::getNbinsX(), MonitorElement::getTH1(), MonitorElement::getTH1F(), MonitorElement::getTH2F(), gen::k, MonitorElement::kind(), dataset::name, pileupCalc::nbins, MonitorElement::setBinContent(), AlCaHLTBitMon_QueryRunRegistry::string, and heppy_batch::val.

Referenced by fillGrandSummaryHistos(), fillSummaryHistos(), and getFrequency().

{
  if (me->getTH1()) {
    TH1F* hist1 = nullptr;
    TH2F* hist2 = nullptr;
    if (me->kind() == MonitorElement::DQM_KIND_TH1F) hist1 = me->getTH1F();
    if (me->kind() == MonitorElement::DQM_KIND_TH2F) hist2 = me->getTH2F();

    int nbins = me_src->getNbinsX();
    const std::string& name = me_src->getName();
    if (htype == "mean" || htype == "Mean" ) {
      if (hist2 &&  name.find("NoisyStrips") != std::string::npos) {
        float bad = 0.0;
        float entries = me_src->getEntries();
        if (entries > 0.0) {
          float binEntry = entries / nbins;
          for (int k = 1; k < nbins + 1; k++) {
            float noisy =
              me_src->getBinContent(k, 3) + me_src->getBinContent(k, 5);
            float dead =
              me_src->getBinContent(k, 2) + me_src->getBinContent(k, 4);
            if (noisy >= binEntry * 0.5 || dead >= binEntry * 0.5) bad++;
          }
          bad = bad * 100.0 / nbins;
          me->Fill(ival, bad);
        }
      } else me->Fill(ival, me_src->getMean());
    } else if (htype == "bin-by-bin" || htype == "Bin-by-Bin") {
      for (int k=1; k<nbins+1; k++) {
        me->setBinContent(istep + k, me_src->getBinContent(k));
      }
    } else if (htype == "sum" || htype == "Sum") {
      if ( hist1) {
        for (int k = 1; k < nbins + 1; k++) {
          float val = me_src->getBinContent(k) + me->getBinContent(k);
          me->setBinContent(k, val);
        }
      }
    }
  }
}

void SiStripSummaryCreator::fillLayout

(

DQMStore &

dqm_store

)

void SiStripSummaryCreator::fillSummaryHistos ( DQMStore &  dqm_store )

private

Definition at line 83 of file SiStripSummaryCreator.cc.

References DQMStore::cd(), relmon_rootfiles_spy::contents, fillHistos(), DQMStore::getContents(), MonitorElement::getName(), MonitorElement::getNbinsX(), DQMStore::getSubdirs(), getSummaryME(), DQMStore::goUp(), isum, dataset::name, DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, and summaryMEs_.

Referenced by createSummary(), and getFrequency().

{
  std::string currDir = dqm_store.pwd();
  std::map<std::string, MonitorElement*> MEMap;
  std::vector<std::string> subdirs = dqm_store.getSubdirs();
  if (subdirs.empty()) return;

  for (std::map<std::string, std::string>::const_iterator isum =
         summaryMEs_.begin();
       isum != summaryMEs_.end();
       isum++) {
    std::string name = isum->first;
    int iBinStep = 0;
    int ndet = 0;
    std::string htype = isum->second;
    for (std::vector<std::string>::const_iterator it = subdirs.begin();
         it != subdirs.end();
         it++) {
      if ( (*it).find("module_") == std::string::npos) continue;
      dqm_store.cd(*it);
      ndet++;
      std::vector<MonitorElement*> contents =
        dqm_store.getContents(dqm_store.pwd());
      dqm_store.goUp();
      for (std::vector<MonitorElement *>::const_iterator im = contents.begin();
                im != contents.end(); im++) {
        MonitorElement * me_i = (*im);
        if (!me_i) continue;
        std::string name_i = me_i->getName();
        if (name_i.find(name) == std::string::npos) continue;
        std::map<std::string, MonitorElement*>::iterator iPos =
          MEMap.find(name);
        MonitorElement* me;
        // Get the Summary ME
        if (iPos == MEMap.end()) {
          me = getSummaryME(dqm_store, name, htype);
          MEMap.insert(std::pair<std::string, MonitorElement*>(name, me));
        } else
          me = iPos->second;
        // Fill it now
        fillHistos(ndet, iBinStep, htype, me_i, me);
        iBinStep += me_i->getNbinsX();
        break;
      }
    }
  }
}

int SiStripSummaryCreator::getFrequency ( ) const

inline

Definition at line 25 of file SiStripSummaryCreator.h.

References fillGrandSummaryHistos(), fillHistos(), fillSummaryHistos(), getSummaryME(), dataset::name, AlCaHLTBitMon_QueryRunRegistry::string, and summaryFrequency_.

  {
    return summaryFrequency_;
  }

MonitorElement * SiStripSummaryCreator::getSummaryME ( DQMStore &  dqm_store, std::string &  name, std::string  htype  )

private

Definition at line 197 of file SiStripSummaryCreator.cc.

References stringResolutionProvider_cfi::bin, DQMStore::book1D(), DQMStore::cd(), relmon_rootfiles_spy::contents, RefreshWebPage::dname, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TPROFILE, DQMStore::getContents(), MonitorElement::getName(), DQMStore::getSubdirs(), MonitorElement::getTH1(), MonitorElement::getTH1F(), DQMStore::goUp(), create_public_lumi_plots::hist, MonitorElement::kind(), label, dataset::name, DQMStore::pwd(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fillGrandSummaryHistos(), fillSummaryHistos(), and getFrequency().

{
  MonitorElement* me = nullptr;
  std::string currDir = dqm_store.pwd();
  std::string sum_name, tag_name;

  std::string dname = currDir.substr(currDir.find_last_of("/")+1);
  if (dname.find("_") != std::string::npos) dname.insert(dname.find("_"),"_");
  if (htype == "sum" && htype == "Sum") {
    sum_name = "Summary" + name + "__" + dname;
    tag_name = "Summary" + name;
  } else {
    sum_name = "Summary_Mean" + name + "__" + dname;
    tag_name = "Summary_Mean" + name;
  }
  // If already booked
  std::vector<MonitorElement*> contents = dqm_store.getContents(currDir);
  for (std::vector<MonitorElement *>::const_iterator im = contents.begin();
                im != contents.end(); im++) {
    MonitorElement * me = (*im);
    if (!me)  continue;
    std::string me_name = me->getName();
    if (me_name.find(sum_name) == 0) {
      if (me->kind() == MonitorElement::DQM_KIND_TH1F ||
          me->kind() == MonitorElement::DQM_KIND_TH2F ||
          me->kind() == MonitorElement::DQM_KIND_TPROFILE) {
        TH1* hist1 = me->getTH1();
        if (hist1) {
          hist1->Reset();
          return me;
        }
      }
    }
  }
  std::map<int, std::string> tags;
  if (!me) {
    int nBins = 0;
    std::vector<std::string> subdirs = dqm_store.getSubdirs();
    // set # of bins of the histogram
    if (htype == "mean" || htype == "Mean" ) {
       nBins = subdirs.size();
       me = dqm_store.book1D(sum_name, sum_name, nBins, 0.5, nBins + 0.5);
       int ibin = 0;
       for (auto const& subdir : subdirs) {
         std::string subdir_name = subdir.substr(subdir.find_last_of("/") + 1);
         ++ibin;
         tags.emplace(ibin, subdir_name);
       }
    } else if (htype == "bin-by-bin" || htype == "Bin-by-Bin") {
      for (auto const& subdir : subdirs) {
        dqm_store.cd(subdir);
        std::string subdir_name = subdir.substr(subdir.find_last_of("/") + 1);
        auto const& s_contents = dqm_store.getContents(dqm_store.pwd());
        for (auto const* s_me : s_contents) {
          if (!s_me) continue;
          std::string s_me_name = s_me->getName();
          if (s_me_name.find(name) == 0 || s_me_name.find(tag_name) == 0) {
            int ibin = s_me->getNbinsX();
            nBins += ibin;
            tags.emplace(nBins - ibin / 2, subdir_name);
            break;
          }
        }
        dqm_store.goUp();
      }
      me = dqm_store.book1D(sum_name, sum_name, nBins, 0.5, nBins + 0.5);
    } else if (htype == "sum" || htype == "Sum") {
      for (auto const& subdir : subdirs) {
        dqm_store.cd(subdir);
        auto const& s_contents = dqm_store.getContents(dqm_store.pwd());
        dqm_store.goUp();
        for (auto* s_me : s_contents) {
          if (!s_me) continue;
          std::string s_me_name = s_me->getName();
          if (s_me_name.find(name) == std::string::npos) continue;
          if (s_me->kind() == MonitorElement::DQM_KIND_TH1F) {
            TH1F* hist1 = s_me->getTH1F();
            if (hist1) {
              nBins = s_me->getNbinsX();
              me = dqm_store.book1D(sum_name,
                                    sum_name,
                                    nBins,
                                    hist1->GetXaxis()->GetXmin(),
                                    hist1->GetXaxis()->GetXmax());
              break;
            }
          }
        }
      }
    }
  }
  // Set the axis title
  if (me && me->kind() == MonitorElement::DQM_KIND_TH1F &&
      (htype != "sum" || htype != "Sum")) {
    TH1F* hist = me->getTH1F();
    if (hist) {
      if (name.find("NoisyStrips") != std::string::npos) hist->GetYaxis()->SetTitle("Noisy Strips (%)");
      else hist->GetYaxis()->SetTitle(name.c_str());

      for (auto const& [bin, label] : tags) {
        hist->GetXaxis()->SetBinLabel(bin, label.c_str());
      }
      hist->LabelsOption("uv");
    }
  }
  return me;
}

bool SiStripSummaryCreator::readConfiguration

(

std::string const &

file_path

)

Definition at line 21 of file SiStripSummaryCreator.cc.

References gather_cfg::cout, DQMParserBase::getDocument(), SiStripConfigParser::getFrequencyForSummary(), SiStripConfigParser::getMENamesForSummary(), summaryFrequency_, and summaryMEs_.

{
  summaryMEs_.clear();
  SiStripConfigParser config_parser;
  config_parser.getDocument(edm::FileInPath(file_path).fullPath());
  if (!config_parser.getFrequencyForSummary(summaryFrequency_)){
    std::cout << "SiStripSummaryCreator::readConfiguration: Failed to read Summary configuration parameters!! ";
    summaryFrequency_ = -1;
    return false;
  }
  if (!config_parser.getMENamesForSummary(summaryMEs_)) {
    std::cout << "SiStripSummaryCreator::readConfiguration: Failed to read Summary configuration parameters!! ";
    return false;
  }
  return true;
}

void SiStripSummaryCreator::setSummaryMENames

(

std::map< std::string, std::string > &

me_names

)

Definition at line 39 of file SiStripSummaryCreator.cc.

References isum, and summaryMEs_.

{
  summaryMEs_.clear();
  for (std::map<std::string, std::string>::const_iterator isum =
         me_names.begin();
       isum != me_names.end();
       isum++) {
    summaryMEs_.insert(
      std::pair<std::string, std::string>(isum->first, isum->second));
  }
}

Member Data Documentation

int SiStripSummaryCreator::summaryFrequency_ {-1}

private

Definition at line 43 of file SiStripSummaryCreator.h.

Referenced by getFrequency(), and readConfiguration().

std::map<std::string, std::string> SiStripSummaryCreator::summaryMEs_

private

Definition at line 42 of file SiStripSummaryCreator.h.

Referenced by createSummary(), fillGrandSummaryHistos(), fillSummaryHistos(), readConfiguration(), and setSummaryMENames().