SiStripSummaryCreator Class Reference

#include <SiStripSummaryCreator.h>

Public Member Functions
void	createSummary (DQMStore *dqm_store)
void	fillLayout (DQMStore *dqm_store)
int	getFrequency ()
bool	readConfiguration (std::string &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_
std::map< std::string, std::string >	summaryMEMap

Detailed Description

Definition at line 14 of file SiStripSummaryCreator.h.

Constructor & Destructor Documentation

SiStripSummaryCreator::SiStripSummaryCreator

(

)

Definition at line 11 of file SiStripSummaryCreator.cc.

References summaryFrequency_, and summaryMEMap.

                                             {
  edm::LogInfo("SiStripSummaryCreator") << 
    " Creating SiStripSummaryCreator " << "\n" ;
  summaryMEMap.clear();
  summaryFrequency_ = -1;
}

SiStripSummaryCreator::~SiStripSummaryCreator ( )

virtual

Definition at line 20 of file SiStripSummaryCreator.cc.

References summaryMEMap.

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

Member Function Documentation

void SiStripSummaryCreator::createSummary

(

DQMStore *

dqm_store

)

Definition at line 57 of file SiStripSummaryCreator.cc.

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

Referenced by SiStripActionExecutor::createSummary(), and SiStripActionExecutor::createSummaryOffline().

                                                             {
  if (summaryMEMap.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 126 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 summaryMEMap.

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 = summaryMEMap.begin();
       isum != summaryMEMap.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 298 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()) {
    //    TProfile* prof = 0;
    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();
    //    if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) prof = me->getTProfile();
    
    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);
//      float good = me_src->getBinContent(k,1);
        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 82 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 summaryMEMap.

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 = summaryMEMap.begin();
       isum != summaryMEMap.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 ( )

inline

Definition at line 26 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 183 of file SiStripSummaryCreator.cc.

References 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(), MonitorElement::getNbinsX(), DQMStore::getSubdirs(), MonitorElement::getTH1(), MonitorElement::getTH1F(), DQMStore::goUp(), create_public_lumi_plots::hist, MonitorElement::kind(), 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 (std::vector<std::string>::const_iterator it = subdirs.begin();
          it != subdirs.end(); it++) {
     std::string subdir_name = (*it).substr((*it).find_last_of("/")+1);
     ibin++;
     tags.insert(std::pair<int,std::string>(ibin, (subdir_name)));        
       }
    } else if (htype == "bin-by-bin" || htype == "Bin-by-Bin") {
      for (std::vector<std::string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
    dqm_store->cd(*it);
        std::string subdir_name = (*it).substr((*it).find_last_of("/")+1);
    std::vector<MonitorElement*> s_contents = dqm_store->getContents(dqm_store->pwd());
    for (std::vector<MonitorElement *>::const_iterator iv = s_contents.begin();
                iv != s_contents.end(); iv++) {
          MonitorElement* s_me = (*iv);
          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.insert(std::pair<int,std::string>(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 (std::vector<std::string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
    dqm_store->cd(*it);
    std::vector<MonitorElement*> s_contents = dqm_store->getContents(dqm_store->pwd());
    dqm_store->goUp();        
    for (std::vector<MonitorElement *>::const_iterator iv = s_contents.begin();
                iv != s_contents.end(); iv++) {
          MonitorElement* s_me = (*iv);
          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 (std::map<int,std::string>::const_iterator ic = tags.begin();
       ic != tags.end(); ic++) {
    hist->GetXaxis()->SetBinLabel(ic->first, (ic->second).c_str());
      }
      hist->LabelsOption("uv");
    }
  }
  return me;
}

bool SiStripSummaryCreator::readConfiguration

(

std::string &

file_path

)

Definition at line 28 of file SiStripSummaryCreator.cc.

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

Referenced by SiStripActionExecutor::readConfiguration().

                                                                  {
    summaryMEMap.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(summaryMEMap)) {
    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 46 of file SiStripSummaryCreator.cc.

References isum, and summaryMEMap.

                                                                                    {

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

Member Data Documentation

int SiStripSummaryCreator::summaryFrequency_

private

Definition at line 41 of file SiStripSummaryCreator.h.

Referenced by getFrequency(), readConfiguration(), and SiStripSummaryCreator().

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

private

Definition at line 39 of file SiStripSummaryCreator.h.

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