SiStripQualityChecker Class Reference

#include <SiStripQualityChecker.h>

Classes
struct	SubDetMEs

Public Member Functions
void	bookStatus (DQMStore *dqm_store)
void	fillDummyStatus ()
void	fillFaultyModuleStatus (DQMStore *dqm_store, const edm::EventSetup &eSetup)
void	fillStatus (DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling, const edm::EventSetup &eSetup)
void	fillStatusAtLumi (DQMStore *dqm_store)
void	printStatusReport ()
void	resetStatus ()
	SiStripQualityChecker (edm::ParameterSet const &ps)
virtual	~SiStripQualityChecker ()

Private Member Functions
void	fillDetectorStatus (DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling)
void	fillDetectorStatusAtLumi (DQMStore *dqm_store)
void	fillStatusHistogram (MonitorElement *, int xbin, int ybin, float val)
void	fillSubDetStatus (DQMStore *dqm_store, const edm::ESHandle< SiStripDetCabling > &cabling, SubDetMEs &mes, unsigned int xbin, float &gflag)
void	getModuleStatus (DQMStore *dqm_store, std::vector< MonitorElement * > &layer_mes, int &errdet, int &errdet_hasBadChan, int &errdet_hasTooManyDigis, int &errdet_hasTooManyClu, int &errdet_hasExclFed, int &errdet_hasDcsErr)
void	initialiseBadModuleList ()

Private Attributes
std::map< uint32_t, uint16_t >	badModuleList
bool	bookedStripStatus_
float	cutoffChi2overDoF_
float	cutoffRecHits_
float	cutoffTrackRate_
MonitorElement *	DetFractionReportMap
MonitorElement *	DetFractionReportMap_hasBadChan
MonitorElement *	DetFractionReportMap_hasDcsErr
MonitorElement *	DetFractionReportMap_hasExclFed
MonitorElement *	DetFractionReportMap_hasTooManyClu
MonitorElement *	DetFractionReportMap_hasTooManyDigis
int	globalStatusFilling_
edm::ParameterSet	pSet_
MonitorElement *	SToNReportMap
std::map< std::string, std::string >	SubDetFolderMap
std::map< std::string, SubDetMEs >	SubDetMEsMap
MonitorElement *	SummaryReportGlobal
MonitorElement *	SummaryReportMap
TkDetMap *	tkDetMap_
MonitorElement *	TrackSummaryReportGlobal
MonitorElement *	TrackSummaryReportMap
bool	useGoodTracks_

Detailed Description

Definition at line 21 of file SiStripQualityChecker.h.

Constructor & Destructor Documentation

SiStripQualityChecker::SiStripQualityChecker

(

edm::ParameterSet const &

ps

)

Definition at line 26 of file SiStripQualityChecker.cc.

References badModuleList, bookedStripStatus_, Utilities::operator, SubDetFolderMap, and tkDetMap_.

                                                                     :pSet_(ps) {
  edm::LogInfo("SiStripQualityChecker") << 
    " Creating SiStripQualityChecker " << "\n" ;

  bookedStripStatus_ = false;

  SubDetFolderMap.insert(std::pair<std::string, std::string>("TIB",  "TIB"));
  SubDetFolderMap.insert(std::pair<std::string, std::string>("TOB",  "TOB"));
  SubDetFolderMap.insert(std::pair<std::string, std::string>("TECF", "TEC/PLUS"));
  SubDetFolderMap.insert(std::pair<std::string, std::string>("TECB", "TEC/MINUS"));
  SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDF", "TID/PLUS"));
  SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDB", "TID/MINUS"));
  badModuleList.clear();

  if(!edm::Service<TkDetMap>().isAvailable()){
    edm::LogError("TkHistoMap") <<
      "\n------------------------------------------"
      "\nUnAvailable Service TkHistoMap: please insert in the configuration file an instance like"
      "\n\tprocess.TkDetMap = cms.Service(\"TkDetMap\")"
      "\n------------------------------------------";
  }
  tkDetMap_=edm::Service<TkDetMap>().operator->();

}

SiStripQualityChecker::~SiStripQualityChecker ( )

virtual

Definition at line 53 of file SiStripQualityChecker.cc.

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

Member Function Documentation

void SiStripQualityChecker::bookStatus

(

DQMStore *

dqm_store

)

Definition at line 60 of file SiStripQualityChecker.cc.

References DQMStore::book2D(), bookedStripStatus_, DQMStore::bookFloat(), DQMStore::cd(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, DetFractionReportMap_hasBadChan, DetFractionReportMap_hasDcsErr, DetFractionReportMap_hasExclFed, DetFractionReportMap_hasTooManyClu, DetFractionReportMap_hasTooManyDigis, SiStripUtility::getTopFolderPath(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, AlCaHLTBitMon_QueryRunRegistry::string, SubDetFolderMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, and SummaryReportMap.

Referenced by SiStripActionExecutor::createStatus(), fillStatus(), and fillStatusAtLumi().

                                                          {

  if (!bookedStripStatus_) {
    dqm_store->cd();
    std::string strip_dir = "";
    SiStripUtility::getTopFolderPath(dqm_store, "SiStrip", strip_dir); 
    if (strip_dir.empty()) strip_dir = "SiStrip";

    // Non Standard Plots and should be put outside EventInfo folder

    dqm_store->setCurrentFolder(strip_dir+"/MechanicalView"); 
      
    std::string hname, htitle;
    hname  = "detFractionReportMap";
    htitle = "SiStrip Report for Good Detector Fraction";
    DetFractionReportMap  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap->setAxisTitle("Layer/Disc Number", 2);

    hname  = "detFractionReportMap_hasBadChan";
    htitle = "SiStrip Report for Good Detector Fraction due to bad channels";
    DetFractionReportMap_hasBadChan  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap_hasBadChan->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap_hasBadChan->setAxisTitle("Layer/Disc Number", 2);
    hname  = "detFractionReportMap_hasTooManyDigis";
    htitle = "SiStrip Report for Good Detector Fraction due to too many digis";
    DetFractionReportMap_hasTooManyDigis  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap_hasTooManyDigis->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap_hasTooManyDigis->setAxisTitle("Layer/Disc Number", 2);
    hname  = "detFractionReportMap_hasTooManyClu";
    htitle = "SiStrip Report for Good Detector Fraction due to too many clusters";
    DetFractionReportMap_hasTooManyClu  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap_hasTooManyClu->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap_hasTooManyClu->setAxisTitle("Layer/Disc Number", 2);
    hname  = "detFractionReportMap_hasExclFed";
    htitle = "SiStrip Report for Good Detector Fraction due to excluded FEDs";
    DetFractionReportMap_hasExclFed  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap_hasExclFed->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap_hasExclFed->setAxisTitle("Layer/Disc Number", 2);
    hname  = "detFractionReportMap_hasDcsErr";
    htitle = "SiStrip Report for Good Detector Fraction due to DCS errors";
    DetFractionReportMap_hasDcsErr  = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    DetFractionReportMap_hasDcsErr->setAxisTitle("Sub Detector Type", 1);
    DetFractionReportMap_hasDcsErr->setAxisTitle("Layer/Disc Number", 2);

    hname  = "sToNReportMap";
    htitle = "SiStrip Report for Signal-to-Noise";
    SToNReportMap         = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    SToNReportMap->setAxisTitle("Sub Detector Type", 1);
    SToNReportMap->setAxisTitle("Layer/Disc Number", 2);

    // this is the main reportSummary 2D plot and should be in EventInfo    
    dqm_store->setCurrentFolder(strip_dir+"/EventInfo"); 

    hname  = "reportSummaryMap";
    htitle = "SiStrip Report Summary Map";
    SummaryReportMap      = dqm_store->book2D(hname, htitle, 6,0.5,6.5,9,0.5,9.5);
    SummaryReportMap->setAxisTitle("Sub Detector Type", 1);
    SummaryReportMap->setAxisTitle("Layer/Disc Number", 2);
    
    SummaryReportGlobal = dqm_store->bookFloat("reportSummary");
    int ibin = 0;
    
    dqm_store->setCurrentFolder(strip_dir+"/EventInfo/reportSummaryContents");      
    for (std::map<std::string, std::string>::const_iterator it = SubDetFolderMap.begin(); 
     it != SubDetFolderMap.end(); it++) {
      ibin++;
      std::string det = it->first;
      DetFractionReportMap->setBinLabel(ibin,it->second);
      DetFractionReportMap_hasBadChan     ->setBinLabel(ibin,it->second);
      DetFractionReportMap_hasTooManyDigis->setBinLabel(ibin,it->second);
      DetFractionReportMap_hasTooManyClu  ->setBinLabel(ibin,it->second);
      DetFractionReportMap_hasExclFed     ->setBinLabel(ibin,it->second);
      DetFractionReportMap_hasDcsErr      ->setBinLabel(ibin,it->second);
      SToNReportMap->setBinLabel(ibin,it->second);
      SummaryReportMap->setBinLabel(ibin,it->second);
      
      SubDetMEs local_mes;
      
      if (det == "TECF")      local_mes.detectorTag = "TEC+";
      else if (det == "TECB") local_mes.detectorTag = "TEC-";         
      else if (det == "TIDF") local_mes.detectorTag = "TID+";
      else if (det == "TIDB") local_mes.detectorTag = "TID-";
      else                    local_mes.detectorTag = det;
      
      std::string me_name;
      me_name = "SiStrip_" + det;
      local_mes.SummaryFlag = dqm_store->bookFloat(me_name);
      
      me_name = "SiStrip_DetFraction_" + det;
      local_mes.DetFraction = dqm_store->bookFloat(me_name);
      
      me_name = "SiStrip_SToNFlag_" + det;
      local_mes.SToNFlag    = dqm_store->bookFloat(me_name);
      SubDetMEsMap.insert(std::pair<std::string, SubDetMEs>(det, local_mes));
    }
    bookedStripStatus_ = true;
  }  
}

void SiStripQualityChecker::fillDetectorStatus ( DQMStore *  dqm_store, const edm::ESHandle< SiStripDetCabling > &  cabling  )

private

Definition at line 228 of file SiStripQualityChecker.cc.

References DQMStore::cd(), DQMStore::dirExists(), RefreshWebPage::dname, MonitorElement::Fill(), fillSubDetStatus(), RemoveAddSevLevel::flag, SiStripUtility::goToDir(), initialiseBadModuleList(), DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, SubDetFolderMap, SubDetMEsMap, and SummaryReportGlobal.

Referenced by fillStatus().

                                                                                                                   {
  unsigned int xbin = 0;
  float global_flag = 0;
  dqm_store->cd();
  std::string mdir = "MechanicalView"; 
  if (!SiStripUtility::goToDir(dqm_store, mdir)) return;
  std::string mechanicalview_dir = dqm_store->pwd();

  initialiseBadModuleList();
  for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin(); 
       it != SubDetMEsMap.end(); it++) {
    std::string det = it->first;
    std::map<std::string, std::string>::const_iterator cPos = SubDetFolderMap.find(det);
    if (cPos == SubDetFolderMap.end()) continue; 
    std::string dname = mechanicalview_dir + "/" + cPos->second;
    if (!dqm_store->dirExists(dname)) continue;
    dqm_store->cd(dname);
    SubDetMEs local_mes = it->second;
    xbin++;
    float flag;
    fillSubDetStatus(dqm_store, cabling, local_mes, xbin,flag);
    global_flag += flag; 
  }
  global_flag = global_flag/xbin*1.0;
  if (SummaryReportGlobal) SummaryReportGlobal->Fill(global_flag);
  dqm_store->cd();
}

void SiStripQualityChecker::fillDetectorStatusAtLumi ( DQMStore *  dqm_store )

private

Definition at line 525 of file SiStripQualityChecker.cc.

References DQMStore::cd(), MonitorElement::DQM_KIND_TH1F, MonitorElement::Fill(), HcalObjRepresent::Fill(), dedxEstimators_cff::fraction, DQMStore::get(), MonitorElement::getTH1F(), SiStripUtility::goToDir(), MonitorElement::kind(), diffTwoXMLs::label, DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, SubDetMEsMap, and SummaryReportGlobal.

Referenced by fillStatusAtLumi().

                                                                       {
  
  

  dqm_store->cd();
  std::string rdir = "ReadoutView"; 
  if (!SiStripUtility::goToDir(dqm_store, rdir)) return;
  std::string fullpath = dqm_store->pwd() 
    //                          + "/FedSummary/PerLumiSection/"
                          + "/PerLumiSection/"
                          + "lumiErrorFraction";  
  MonitorElement* me = dqm_store->get(fullpath);
  if (me && me->kind() == MonitorElement::DQM_KIND_TH1F) {
    TH1F* th1 = me->getTH1F(); 
    float global_fraction = 0.0;
    float dets = 0.0;
    for (int ibin = 1; ibin <= th1->GetNbinsX(); ibin++) {
      std::string label = th1->GetXaxis()->GetBinLabel(ibin);
      std::map<std::string, SubDetMEs>::iterator iPos = SubDetMEsMap.find(label);
      if (iPos != SubDetMEsMap.end()) {
        float fraction = 1.0 - th1->GetBinContent(ibin);
        global_fraction +=  fraction;
        dets++; 
        iPos->second.DetFraction -> Fill(fraction);
        iPos->second.SToNFlag    -> Fill(-1.0);
        iPos->second.SummaryFlag -> Fill(fraction);
      }
    }
    global_fraction = global_fraction/dets;
    if (SummaryReportGlobal) SummaryReportGlobal->Fill(global_fraction);    
  }
  dqm_store->cd();
}

void SiStripQualityChecker::fillDummyStatus

(

)

Definition at line 162 of file SiStripQualityChecker.cc.

References bookedStripStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, DetFractionReportMap_hasBadChan, DetFractionReportMap_hasDcsErr, DetFractionReportMap_hasExclFed, DetFractionReportMap_hasTooManyClu, DetFractionReportMap_hasTooManyDigis, MonitorElement::Fill(), HcalObjRepresent::Fill(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), resetStatus(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, and SummaryReportMap.

Referenced by SiStripActionExecutor::fillDummyStatus(), fillStatus(), and fillStatusAtLumi().

                                           {
 
  resetStatus();
  if (bookedStripStatus_) {
    for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin(); 
     it != SubDetMEsMap.end(); it++) {
      SubDetMEs local_mes = it->second;
      local_mes.SummaryFlag -> Fill(-1.0);
      local_mes.DetFraction -> Fill(-1.0);
      local_mes.SToNFlag    -> Fill(-1.0);
    }
    
    for (int xbin = 1; xbin < SummaryReportMap->getNbinsX()+1; xbin++) {
      for (int ybin = 1; ybin < SummaryReportMap->getNbinsY()+1; ybin++) {
    SummaryReportMap     -> Fill(xbin, ybin, -1.0);
    DetFractionReportMap -> Fill(xbin, ybin, -1.0);
  DetFractionReportMap_hasBadChan      -> Fill(xbin, ybin, -1.0);
    DetFractionReportMap_hasTooManyDigis -> Fill(xbin, ybin, -1.0);
    DetFractionReportMap_hasTooManyClu   -> Fill(xbin, ybin, -1.0);
    DetFractionReportMap_hasExclFed      -> Fill(xbin, ybin, -1.0);
    DetFractionReportMap_hasDcsErr       -> Fill(xbin, ybin, -1.0);
    SToNReportMap        -> Fill(xbin, ybin, -1.0);
      }
    }
    SummaryReportGlobal->Fill(-1.0);
  }
}

void SiStripQualityChecker::fillFaultyModuleStatus	(	DQMStore *	dqm_store,
		const edm::EventSetup &	eSetup
	)

Definition at line 471 of file SiStripQualityChecker.cc.

References badModuleList, DQMStore::bookInt(), DQMStore::cd(), DQMStore::dirExists(), MonitorElement::Fill(), edm::EventSetup::get(), DQMStore::get(), SiStripFolderOrganizer::getSubDetFolder(), SiStripUtility::goToDir(), edm::ESHandle< T >::product(), DQMStore::pwd(), MonitorElement::Reset(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by fillStatus().

                                                                                                   {
  if (badModuleList.empty()) return;

  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology* const tTopo = tTopoHandle.product();

  dqm_store->cd();
  std::string mdir = "MechanicalView";
  if (!SiStripUtility::goToDir(dqm_store, mdir)) return;
  std::string mechanical_dir = dqm_store->pwd();

  SiStripFolderOrganizer folder_organizer;
  for (std::map<uint32_t,uint16_t>::const_iterator it =  badModuleList.begin() ; it != badModuleList.end(); it++) {
    uint32_t detId =  it->first;
    std::string subdet_folder ;
    folder_organizer.getSubDetFolder(detId,tTopo,subdet_folder);
    if (!dqm_store->dirExists(subdet_folder)) {
      subdet_folder = mechanical_dir + subdet_folder.substr(subdet_folder.find("MechanicalView")+14);
      if (!dqm_store->dirExists(subdet_folder)) continue;
    }
    std::string bad_module_folder = subdet_folder + "/" + "BadModuleList";
    dqm_store->setCurrentFolder(bad_module_folder);

    std::ostringstream detid_str;
    detid_str << detId;
    std::string full_path = bad_module_folder + "/" + detid_str.str();
    MonitorElement* me = dqm_store->get(full_path);
    if (me) me->Reset();
    else me = dqm_store->bookInt(detid_str.str());
    me->Fill(it->second);
  }
  dqm_store->cd();
}

void SiStripQualityChecker::fillStatus	(	DQMStore *	dqm_store,
		const edm::ESHandle< SiStripDetCabling > &	cabling,
		const edm::EventSetup &	eSetup
	)

Definition at line 216 of file SiStripQualityChecker.cc.

References bookedStripStatus_, bookStatus(), fillDetectorStatus(), fillDummyStatus(), fillFaultyModuleStatus(), edm::ParameterSet::getUntrackedParameter(), and pSet_.

Referenced by SiStripActionExecutor::fillStatus().

                                                                                                                                        {
  if (!bookedStripStatus_) bookStatus(dqm_store);

  fillDummyStatus();
  fillDetectorStatus(dqm_store, cabling);

  int faulty_moduleflag  = pSet_.getUntrackedParameter<bool>("PrintFaultyModuleList", false);
  if (faulty_moduleflag) fillFaultyModuleStatus(dqm_store, eSetup);   
}

void SiStripQualityChecker::fillStatusAtLumi

(

DQMStore *

dqm_store

)

Definition at line 517 of file SiStripQualityChecker.cc.

References bookedStripStatus_, bookStatus(), fillDetectorStatusAtLumi(), and fillDummyStatus().

Referenced by SiStripActionExecutor::fillStatusAtLumi().

                                                               {
  if (!bookedStripStatus_) bookStatus(dqm_store);
  fillDummyStatus();
  fillDetectorStatusAtLumi(dqm_store);
}

void SiStripQualityChecker::fillStatusHistogram ( MonitorElement *  me, int  xbin, int  ybin, float  val  )

private

Definition at line 462 of file SiStripQualityChecker.cc.

References MonitorElement::DQM_KIND_TH2F, MonitorElement::getTH2F(), and MonitorElement::kind().

Referenced by fillSubDetStatus().

                                                                                                 {
   if (me &&  me->kind() == MonitorElement::DQM_KIND_TH2F) {
     TH2F*  th2d = me->getTH2F();
     th2d->SetBinContent(xbin, ybin, val);
   }
 }

void SiStripQualityChecker::fillSubDetStatus ( DQMStore *  dqm_store, const edm::ESHandle< SiStripDetCabling > &  cabling, SubDetMEs &  mes, unsigned int  xbin, float &  gflag  )

private

Definition at line 258 of file SiStripQualityChecker.cc.

References DQMStore::cd(), SiStripDetCabling::connectedNumber(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, DetFractionReportMap_hasBadChan, DetFractionReportMap_hasDcsErr, DetFractionReportMap_hasExclFed, DetFractionReportMap_hasTooManyClu, DetFractionReportMap_hasTooManyDigis, RefreshWebPage::dname, dqm::qstatus::ERROR, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), SiStripUtility::getMEStatus(), getModuleStatus(), MonitorElement::getName(), MonitorElement::getQReports(), DQMStore::getSubdirs(), edm::ParameterSet::getUntrackedParameter(), min(), dataset::name, pSet_, dqm::qstatus::STATUS_OK, SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, AlCaHLTBitMon_QueryRunRegistry::string, SiStripQualityChecker::SubDetMEs::SummaryFlag, and SummaryReportMap.

Referenced by fillDetectorStatus().

                                                                          {

  int status_flag  = pSet_.getUntrackedParameter<int>("GlobalStatusFilling", 1);
  if (status_flag < 1) return;

  std::vector<std::string> subDirVec = dqm_store->getSubdirs();

  unsigned int ybin   = 0;
  int tot_ndet        = 0;
  int tot_errdet      = 0;
  float tot_ston_stat = 0;

  for (std::vector<std::string>::const_iterator ic = subDirVec.begin();
       ic != subDirVec.end(); ic++) {
    std::string dname = (*ic);
    if (dname.find("BadModuleList") != std::string::npos) continue;
    if (dname.find("ring") !=std::string::npos) continue;
    std::vector<MonitorElement*> meVec;
    
    ybin++;
    dqm_store->cd((*ic));
    meVec = dqm_store->getContents((*ic));
    uint16_t ndet = 100;
    int errdet = 0;
    int errdet_hasBadChan = 0;
    int errdet_hasTooManyDigis = 0;
    int errdet_hasTooManyClu = 0;
    int errdet_hasExclFed = 0;
    int errdet_hasDcsErr = 0;      

    int ston_stat = 1;
    int lnum = atoi(dname.substr(dname.find_last_of("_")+1).c_str());
    ndet = cabling->connectedNumber(mes.detectorTag, lnum);
     
    getModuleStatus(dqm_store, meVec, errdet, errdet_hasBadChan, errdet_hasTooManyDigis, errdet_hasTooManyClu, errdet_hasExclFed, errdet_hasDcsErr);

    for (std::vector<MonitorElement*>::const_iterator it = meVec.begin();
     it != meVec.end(); it++) {
      MonitorElement * me = (*it);
      if (!me) continue;
      std::vector<QReport *> reports = me->getQReports();

      if (reports.empty()) continue;
      std::string name = me->getName();
      
      if( name.find("Summary_ClusterStoNCorr__OnTrack") != std::string::npos){
    int istat =  SiStripUtility::getMEStatus((*it)); 
        if (reports[0]->getQTresult() == -1) {
      ston_stat =-1;
        } else {
          if (istat == dqm::qstatus::ERROR) ston_stat = 0;
          else if (istat == dqm::qstatus::STATUS_OK) ston_stat = 1;
        }
      }
    }
    if (ndet > 0) {
      float eff_fac = 1 - (errdet*1.0/ndet);
      float eff_fac_hasBadChan      = 1 - (errdet_hasBadChan     *1.0/ndet);
      float eff_fac_hasTooManyDigis = 1 - (errdet_hasTooManyDigis*1.0/ndet);
      float eff_fac_hasTooManyClu   = 1 - (errdet_hasTooManyClu  *1.0/ndet);
      float eff_fac_hasExclFed      = 1 - (errdet_hasExclFed     *1.0/ndet);
      float eff_fac_hasDcsErr       = 1 - (errdet_hasDcsErr      *1.0/ndet);
      fillStatusHistogram(SToNReportMap,        xbin, ybin, ston_stat);
      fillStatusHistogram(DetFractionReportMap, xbin, ybin, eff_fac);
      fillStatusHistogram(DetFractionReportMap_hasBadChan     , xbin, ybin, eff_fac_hasBadChan     );
      fillStatusHistogram(DetFractionReportMap_hasTooManyDigis, xbin, ybin, eff_fac_hasTooManyDigis);
      fillStatusHistogram(DetFractionReportMap_hasTooManyClu  , xbin, ybin, eff_fac_hasTooManyClu  );
      fillStatusHistogram(DetFractionReportMap_hasExclFed     , xbin, ybin, eff_fac_hasExclFed     );
      fillStatusHistogram(DetFractionReportMap_hasDcsErr      , xbin, ybin, eff_fac_hasDcsErr      );
      if (ston_stat < 0) fillStatusHistogram(SummaryReportMap, xbin, ybin, eff_fac);
      else               fillStatusHistogram(SummaryReportMap, xbin, ybin, ston_stat*eff_fac);

      tot_ndet      += ndet;
      tot_errdet    += errdet;
      tot_ston_stat += ston_stat;  
    }
    dqm_store->cd((*ic));
  }
  if (tot_ndet > 0) { 
    float tot_eff_fac = 1 - (tot_errdet*1.0/tot_ndet);
    if (mes.DetFraction) mes.DetFraction->Fill(tot_eff_fac);
    float tot_ston_fac = tot_ston_stat/ybin;
    if (mes.SToNFlag) mes.SToNFlag->Fill(tot_ston_fac);
    if (tot_ston_fac < 0){
      gflag = tot_eff_fac;    
    }else{
      gflag = std::min(tot_eff_fac,tot_ston_fac);    
    }
    if (mes.SummaryFlag) mes.SummaryFlag->Fill(gflag);
  }
}    

void SiStripQualityChecker::getModuleStatus ( DQMStore *  dqm_store, std::vector< MonitorElement * > &  layer_mes, int &  errdet, int &  errdet_hasBadChan, int &  errdet_hasTooManyDigis, int &  errdet_hasTooManyClu, int &  errdet_hasExclFed, int &  errdet_hasDcsErr  )

private

Definition at line 385 of file SiStripQualityChecker.cc.

References badModuleList, lumiPlot::currentdir, MonitorElement::DQM_KIND_TPROFILE, MonitorElement::DQM_KIND_TPROFILE2D, RemoveAddSevLevel::flag, DQMStore::get(), TkDetMap::getDetFromBin(), SiStripUtility::getMEValue(), MonitorElement::getName(), MonitorElement::getQReports(), MonitorElement::kind(), dataset::name, DQMStore::pwd(), SiStripUtility::setBadModuleFlag(), AlCaHLTBitMon_QueryRunRegistry::string, tkDetMap_, and globals_cff::x1.

Referenced by fillSubDetStatus().

                                                                                                                                                                                                                                                 { 
  
  std::string lname;
  std::map<uint32_t,uint16_t> bad_modules;
  for (std::vector<MonitorElement*>::const_iterator it = layer_mes.begin();
       it != layer_mes.end(); it++) {
    MonitorElement * me = (*it);
    if (!me) continue;
    std::vector<QReport *> qreports = me->getQReports();
    if (qreports.empty()) continue;
    std::string name = me->getName();
    std::vector<DQMChannel>  bad_channels_me;
    if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) {
      bad_channels_me = qreports[0]->getBadChannels();
      lname = "";
    } else if (me->kind() == MonitorElement::DQM_KIND_TPROFILE2D && name.find("TkHMap") != std::string::npos) {
      bad_channels_me = qreports[0]->getBadChannels();
      lname = name.substr(name.find("TkHMap_")+7);
      lname = lname.substr(lname.find("_T")+1);

    }
    for (std::vector<DQMChannel>::iterator it = bad_channels_me.begin(); it != bad_channels_me.end(); it++){
      int xval = (*it).getBinX();
      int yval = (*it).getBinY();
      uint32_t detId = tkDetMap_->getDetFromBin(lname, xval, yval);       
      std::map<uint32_t,uint16_t>::iterator iPos = bad_modules.find(detId);
      uint16_t flag;
      if (iPos != bad_modules.end()){
    flag = iPos->second;
    SiStripUtility::setBadModuleFlag(name,flag);            
    iPos->second = flag;
      } else {
        //  
    //if not in the local bad module list, check the BadModuleList dir  
    //  
    std::ostringstream detid_str;  
    detid_str << detId;  
    //now in the layer/wheel dir  
    const std::string& currentdir = dqm_store->pwd();  
    std::string thisMEpath = currentdir.substr( 0 , currentdir.rfind( "/" ) ) + "/BadModuleList/" + detid_str.str() ;  
    
    MonitorElement *meBadModule = dqm_store->get ( thisMEpath );  
    if ( meBadModule )  
      {  
        std::string val_str;  
        SiStripUtility::getMEValue ( meBadModule , val_str );  
        flag = atoi ( val_str.c_str() );  
      }  
    else
      flag = 0;
    
    SiStripUtility::setBadModuleFlag(name,flag);              
    bad_modules.insert(std::pair<uint32_t,uint16_t>(detId,flag));
      }
    }
  }
  for(std::map<uint32_t,uint16_t>::const_iterator it = bad_modules.begin();
      it != bad_modules.end(); it++) {
    uint32_t detId = it->first;
    uint16_t flag  = it->second;
    if (((flag >> 0) & 0x1) > 0) errdet_hasBadChan++;
    if (((flag >> 1) & 0x1) > 0) errdet_hasTooManyDigis++;
    if (((flag >> 2) & 0x1) > 0) errdet_hasTooManyClu++;
    if (((flag >> 3) & 0x1) > 0) errdet_hasExclFed++;
    if (((flag >> 4) & 0x1) > 0) errdet_hasDcsErr++;
    std::map<uint32_t,uint16_t>::iterator iPos = badModuleList.find(detId);
    if (iPos != badModuleList.end()){
      iPos->second = flag;
    } else {
      badModuleList.insert(std::pair<uint32_t,uint16_t>(detId,flag));
    }
  }    
  errdet = bad_modules.size();
}

void SiStripQualityChecker::initialiseBadModuleList ( )

private

Definition at line 509 of file SiStripQualityChecker.cc.

References badModuleList.

Referenced by fillDetectorStatus().

                                                    {
  for (std::map<uint32_t,uint16_t>::iterator it=badModuleList.begin(); it!=badModuleList.end(); it++) {
    it->second = 0;
  }
}

void SiStripQualityChecker::printStatusReport

(

)

Definition at line 354 of file SiStripQualityChecker.cc.

References SiStripQualityChecker::SubDetMEs::DetFraction, alignBH_cfg::fixed, SiStripUtility::getMEValue(), SiStripQualityChecker::SubDetMEs::SToNFlag, AlCaHLTBitMon_QueryRunRegistry::string, SubDetMEsMap, and SiStripQualityChecker::SubDetMEs::SummaryFlag.

                                              {
  std::ostringstream det_summary_str;
  for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin(); 
       it != SubDetMEsMap.end(); it++) {
    std::string det = it->first;
    det_summary_str << std::setprecision(4);
    det_summary_str << std::setiosflags(std::ios::fixed);

    det_summary_str << " Printing Status for " <<   det << " : " << std::endl;
    SubDetMEs local_mes = it->second;

    std::string sval;
    float fval1, fval2, fval3;
    fval1 = fval2 = fval3 = -1.0;   
    

    SiStripUtility::getMEValue(local_mes.DetFraction, sval); 
    if (!sval.empty()) fval1 = atof(sval.c_str());
    SiStripUtility::getMEValue(local_mes.SToNFlag, sval); 
    if (!sval.empty()) fval2 = atof(sval.c_str());
    SiStripUtility::getMEValue(local_mes.SummaryFlag, sval); 
    if (!sval.empty()) fval3 = atof(sval.c_str());

    det_summary_str << std::setw(7) << " % of good detectors " << fval1
            << " SToN Flag           " << fval2
            << " Summary Flag        " << fval3 << std::endl;
  }
}

void SiStripQualityChecker::resetStatus

(

)

Definition at line 192 of file SiStripQualityChecker.cc.

References bookedStripStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, DetFractionReportMap_hasBadChan, DetFractionReportMap_hasDcsErr, DetFractionReportMap_hasExclFed, DetFractionReportMap_hasTooManyClu, DetFractionReportMap_hasTooManyDigis, MonitorElement::Reset(), HcalObjRepresent::Reset(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, and SummaryReportMap.

Referenced by fillDummyStatus().

                                        {
  if (bookedStripStatus_) {
    for (std::map<std::string, SubDetMEs>::const_iterator it = SubDetMEsMap.begin(); 
     it != SubDetMEsMap.end(); it++) {
      SubDetMEs local_mes = it->second;
      local_mes.DetFraction -> Reset();
      local_mes.SummaryFlag -> Reset();
      local_mes.SToNFlag    -> Reset();
    }
    SummaryReportMap->Reset();
    DetFractionReportMap->Reset();
    DetFractionReportMap_hasBadChan     ->Reset();
    DetFractionReportMap_hasTooManyDigis->Reset();
    DetFractionReportMap_hasTooManyClu  ->Reset();
    DetFractionReportMap_hasExclFed     ->Reset();
    DetFractionReportMap_hasDcsErr      ->Reset();
    SToNReportMap->Reset();

    SummaryReportGlobal->Reset();
  }
}

Member Data Documentation

std::map<uint32_t,uint16_t> SiStripQualityChecker::badModuleList

private

Definition at line 74 of file SiStripQualityChecker.h.

Referenced by fillFaultyModuleStatus(), getModuleStatus(), initialiseBadModuleList(), and SiStripQualityChecker().

bool SiStripQualityChecker::bookedStripStatus_

private

Definition at line 78 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillStatus(), fillStatusAtLumi(), resetStatus(), and SiStripQualityChecker().

float SiStripQualityChecker::cutoffChi2overDoF_

private

Definition at line 85 of file SiStripQualityChecker.h.

float SiStripQualityChecker::cutoffRecHits_

private

Definition at line 86 of file SiStripQualityChecker.h.

float SiStripQualityChecker::cutoffTrackRate_

private

Definition at line 84 of file SiStripQualityChecker.h.

MonitorElement* SiStripQualityChecker::DetFractionReportMap

private

Definition at line 59 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::DetFractionReportMap_hasBadChan

private

Definition at line 60 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::DetFractionReportMap_hasDcsErr

private

Definition at line 64 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::DetFractionReportMap_hasExclFed

private

Definition at line 63 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::DetFractionReportMap_hasTooManyClu

private

Definition at line 62 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::DetFractionReportMap_hasTooManyDigis

private

Definition at line 61 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

int SiStripQualityChecker::globalStatusFilling_

private

Definition at line 79 of file SiStripQualityChecker.h.

edm::ParameterSet SiStripQualityChecker::pSet_

private

Definition at line 76 of file SiStripQualityChecker.h.

Referenced by fillStatus(), and fillSubDetStatus().

MonitorElement* SiStripQualityChecker::SToNReportMap

private

Definition at line 65 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

std::map<std::string, std::string> SiStripQualityChecker::SubDetFolderMap

private

Definition at line 57 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDetectorStatus(), and SiStripQualityChecker().

std::map<std::string, SubDetMEs> SiStripQualityChecker::SubDetMEsMap

private

Definition at line 56 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDetectorStatus(), fillDetectorStatusAtLumi(), fillDummyStatus(), printStatusReport(), and resetStatus().

MonitorElement* SiStripQualityChecker::SummaryReportGlobal

private

Definition at line 68 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDetectorStatus(), fillDetectorStatusAtLumi(), fillDummyStatus(), and resetStatus().

MonitorElement* SiStripQualityChecker::SummaryReportMap

private

Definition at line 66 of file SiStripQualityChecker.h.

Referenced by bookStatus(), fillDummyStatus(), fillSubDetStatus(), and resetStatus().

TkDetMap* SiStripQualityChecker::tkDetMap_

private

Definition at line 82 of file SiStripQualityChecker.h.

Referenced by getModuleStatus(), and SiStripQualityChecker().

MonitorElement* SiStripQualityChecker::TrackSummaryReportGlobal

private

Definition at line 72 of file SiStripQualityChecker.h.

MonitorElement* SiStripQualityChecker::TrackSummaryReportMap

private

Definition at line 70 of file SiStripQualityChecker.h.

bool SiStripQualityChecker::useGoodTracks_

private

Definition at line 80 of file SiStripQualityChecker.h.