#include <SiStripQualityChecker.h>

Classes
struct	SubDetMEs

struct	TrackingMEs

Public Member Functions
void	bookStatus (DQMStore *dqm_store)

void	fillDummyStatus ()

void	fillFaultyModuleStatus (DQMStore *dqm_store)

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

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	fillTrackingStatus (DQMStore *dqm_store)

void	fillTrackingStatusAtLumi (DQMStore *dqm_store)

void	getModuleStatus (DQMStore dqm_store, std::vector< MonitorElement > &layer_mes, int &errdet)

void	initialiseBadModuleList ()

Private Attributes
std::map< uint32_t, uint16_t >	badModuleList

bool	bookedStripStatus_

bool	bookedTrackingStatus_

float	cutoffChi2overDoF_

float	cutoffRecHits_

float	cutoffTrackRate_

MonitorElement *	DetFractionReportMap

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_

std::map< std::string, TrackingMEs >	TrackingMEsMap

MonitorElement *	TrackSummaryReportGlobal

MonitorElement *	TrackSummaryReportMap

bool	useGoodTracks_

Detailed Description

Definition at line 20 of file SiStripQualityChecker.h.

Constructor & Destructor Documentation

SiStripQualityChecker::SiStripQualityChecker ( edm::ParameterSet const & ps )

Definition at line 28 of file SiStripQualityChecker.cc.

References badModuleList, bookedStripStatus_, bookedTrackingStatus_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), SiStripQualityChecker::TrackingMEs::HistoName, SiStripQualityChecker::TrackingMEs::LowerCut, cppFunctionSkipper::operator, pSet_, SubDetFolderMap, tkDetMap_, SiStripQualityChecker::TrackingMEs::TrackingFlag, TrackingMEsMap, SiStripQualityChecker::TrackingMEs::UpperCut, and useGoodTracks_.

                                                                      :pSet_(ps) {
   edm::LogInfo("SiStripQualityChecker") << 
     " Creating SiStripQualityChecker " << "\n" ;
 
   bookedStripStatus_ = false;
   bookedTrackingStatus_ = 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/side_2"));
   SubDetFolderMap.insert(std::pair<std::string, std::string>("TECB", "TEC/side_1"));
   SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDF", "TID/side_2"));
   SubDetFolderMap.insert(std::pair<std::string, std::string>("TIDB", "TID/side_1"));
   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->();
 
   TrackingMEs tracking_mes;
   edm::ParameterSet TkPSet;
  
   TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackRatePSet"); 
   tracking_mes.TrackingFlag = 0;
   tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
   tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
   tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
   TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("Rate", tracking_mes));
 
   TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackChi2PSet"); 
   tracking_mes.TrackingFlag = 0;
   tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
   tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
   tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
   TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("Chi2", tracking_mes));
 
   TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackHitPSet"); 
   tracking_mes.TrackingFlag = 0;
   tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
   tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
   tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
   TrackingMEsMap.insert(std::pair<std::string, TrackingMEs>("RecHits", tracking_mes));
 
 //  // LS analysis
 //  TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackRateLSPSet"); 
 //  tracking_mes.TrackingFlag = 0;
 //  tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
 //  tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
 //  tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
 //  TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("Rate", tracking_mes));
 //
 //  TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackChi2LSPSet"); 
 //  tracking_mes.TrackingFlag = 0;
 //  tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
 //  tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
 //  tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
 //  TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("Chi2", tracking_mes));
 //
 //  TkPSet = pSet_.getParameter<edm::ParameterSet>("TrackHitLSPSet"); 
 //  tracking_mes.TrackingFlag = 0;
 //  tracking_mes.HistoName    = TkPSet.getParameter<std::string>("Name");
 //  tracking_mes.LowerCut = TkPSet.getParameter<double>("LowerCut");
 //  tracking_mes.UpperCut = TkPSet.getParameter<double>("UpperCut");
 //  TrackingMEsLSMap.insert(std::pair<std::string, TrackingMEs>("RecHits", tracking_mes));
 
   useGoodTracks_  = pSet_.getUntrackedParameter<bool>("UseGoodTracks", false);
   if (useGoodTracks_) edm::LogInfo("SiStripQualityChecker") <<  " use GoodTrack histograms for certification " << "\n" ;
 }

SiStripQualityChecker::~SiStripQualityChecker ( )

virtual

Definition at line 104 of file SiStripQualityChecker.cc.

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

Member Function Documentation

void SiStripQualityChecker::bookStatus ( DQMStore * dqm_store )

Definition at line 111 of file SiStripQualityChecker.cc.

References DQMStore::book2D(), bookedStripStatus_, bookedTrackingStatus_, DQMStore::bookFloat(), DQMStore::cd(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, SiStripUtility::getTopFolderPath(), mergeVDriftHistosByStation::name, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetFolderMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.

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.size() == 0) 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  = "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,det);
       SToNReportMap->setBinLabel(ibin,det);
       SummaryReportMap->setBinLabel(ibin,det);
       
       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;
   }  
   if (!bookedTrackingStatus_) {
     dqm_store->cd();     
     std::string tracking_dir = "";
     SiStripUtility::getTopFolderPath(dqm_store, "Tracking", tracking_dir);
     if (tracking_dir.size() ==  0) tracking_dir = "Tracking";
     dqm_store->setCurrentFolder(tracking_dir+"/EventInfo"); 
     
     TrackSummaryReportGlobal = dqm_store->bookFloat("reportSummary");
     
     std::string hname, htitle;
     hname  = "reportSummaryMap";
     htitle = "Tracking Report Summary Map";
     
     TrackSummaryReportMap    = dqm_store->book2D(hname, htitle, 3,0.5,3.5,1,0.5,1.5);
     TrackSummaryReportMap->setAxisTitle("Track Quality Type", 1);
     TrackSummaryReportMap->setAxisTitle("QTest Flag", 2);
     TrackSummaryReportMap->setBinLabel(1,"Rate");
     TrackSummaryReportMap->setBinLabel(2,"Chi2");
     TrackSummaryReportMap->setBinLabel(3,"RecHits");
     
     dqm_store->setCurrentFolder(tracking_dir+"/EventInfo/reportSummaryContents");  
     int ibin = 0;
     for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
          it != TrackingMEsMap.end(); it++) {
       ibin++;
       std::string name = it->first;
       it->second.TrackingFlag = dqm_store->bookFloat("Track"+name);
       TrackSummaryReportMap->setBinLabel(ibin,name);
     }
     bookedTrackingStatus_ = true;
     dqm_store->cd();
   }
 }

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

private

Definition at line 303 of file SiStripQualityChecker.cc.

References DQMStore::cd(), DQMStore::dirExists(), RefreshWebPage::dname, MonitorElement::Fill(), fillSubDetStatus(), SiStripUtility::goToDir(), initialiseBadModuleList(), DQMStore::pwd(), 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 623 of file SiStripQualityChecker.cc.

References DQMStore::cd(), MonitorElement::DQM_KIND_TH1F, MonitorElement::Fill(), DQMStore::get(), MonitorElement::getTH1F(), SiStripUtility::goToDir(), MonitorElement::kind(), diffTwoXMLs::label, DQMStore::pwd(), 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/"
                           + "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 213 of file SiStripQualityChecker.cc.

References bookedStripStatus_, bookedTrackingStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, MonitorElement::Fill(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), resetStatus(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.

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.DetFraction->Fill(-1.0);
       local_mes.SToNFlag->Fill(-1.0);
       local_mes.SummaryFlag->Fill(-1.0);
     }
     
     for (int xbin = 1; xbin < SummaryReportMap->getNbinsX()+1; xbin++) {
       for (int ybin = 1; ybin < SummaryReportMap->getNbinsY()+1; ybin++) {
     DetFractionReportMap->Fill(xbin, ybin, -1.0);
     SToNReportMap->Fill(xbin, ybin, -1.0);
     SummaryReportMap->Fill(xbin, ybin, -1.0);
       }
     }
     SummaryReportGlobal->Fill(-1.0);
   }
   if (bookedTrackingStatus_) {  
     TrackSummaryReportGlobal->Fill(-1.0);
     for (int xbin = 1; xbin < TrackSummaryReportMap->getNbinsX()+1; xbin++) {
       for (int ybin = 1; ybin < TrackSummaryReportMap->getNbinsY()+1; ybin++) {
         TrackSummaryReportMap->Fill(xbin, ybin, -1.0);
       }
     }
     for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
          it != TrackingMEsMap.end(); it++) {
       it->second.TrackingFlag->Fill(-1.0);
     }
     /*
     for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsLSMap.begin();
          it != TrackingMEsLSMap.end(); it++) {
       it->second.TrackingFlag->Fill(-1.0);
     }
     */
   }
 }

void SiStripQualityChecker::fillFaultyModuleStatus ( DQMStore * dqm_store )

Definition at line 574 of file SiStripQualityChecker.cc.

References badModuleList, DQMStore::bookInt(), DQMStore::cd(), DQMStore::dirExists(), MonitorElement::Fill(), DQMStore::get(), SiStripFolderOrganizer::getSubDetFolder(), SiStripUtility::goToDir(), DQMStore::pwd(), MonitorElement::Reset(), and DQMStore::setCurrentFolder().

Referenced by fillStatus().

                                                                       {
   if (badModuleList.size() == 0) return;
   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,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
	)

Definition at line 290 of file SiStripQualityChecker.cc.

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

Referenced by SiStripActionExecutor::fillStatus().

                                                                                                            {
   if (!bookedStripStatus_ || !bookedTrackingStatus_) bookStatus(dqm_store);
 
   fillDummyStatus();
   fillDetectorStatus(dqm_store, cabling);
   fillTrackingStatus(dqm_store);
 
   int faulty_moduleflag  = pSet_.getUntrackedParameter<bool>("PrintFaultyModuleList", false);
   if (faulty_moduleflag) fillFaultyModuleStatus(dqm_store);   
 }

void SiStripQualityChecker::fillStatusAtLumi ( DQMStore * dqm_store )

Definition at line 614 of file SiStripQualityChecker.cc.

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

Referenced by SiStripActionExecutor::fillStatusAtLumi().

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

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

private

Definition at line 565 of file SiStripQualityChecker.cc.

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

Referenced by fillSubDetStatus(), fillTrackingStatus(), and fillTrackingStatusAtLumi().

                                                                                                  {
    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 383 of file SiStripQualityChecker.cc.

References DQMStore::cd(), SiStripQualityChecker::SubDetMEs::detectorTag, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, RefreshWebPage::dname, dqm::qstatus::ERROR, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), SiStripUtility::getMEStatus(), getModuleStatus(), MonitorElement::getName(), MonitorElement::getQReports(), DQMStore::getSubdirs(), edm::ParameterSet::getUntrackedParameter(), min, mergeVDriftHistosByStation::name, pSet_, dqm::qstatus::STATUS_OK, SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, 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;
     std::vector<MonitorElement*> meVec;
     ybin++;
     dqm_store->cd((*ic));
     meVec = dqm_store->getContents((*ic));
     uint16_t ndet = 100;
     int errdet = 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);
 
     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.size() == 0) 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);
       fillStatusHistogram(SToNReportMap,        xbin, ybin, ston_stat);
       fillStatusHistogram(DetFractionReportMap, xbin, ybin, eff_fac);
       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::fillTrackingStatus ( DQMStore * dqm_store )

private

Definition at line 333 of file SiStripQualityChecker.cc.

References DQMStore::cd(), dir, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), MonitorElement::getName(), MonitorElement::getQReports(), SiStripUtility::goToDir(), relval_steps::merge(), mergeVDriftHistosByStation::name, DQMStore::pwd(), ntuplemaker::status, TrackingMEsMap, TrackSummaryReportGlobal, TrackSummaryReportMap, and useGoodTracks_.

Referenced by fillStatus().

                                                                   {
 
   dqm_store->cd();
   std::string dir = "Tracking"; 
   if (!SiStripUtility::goToDir(dqm_store, dir)) return;
   dir = "TrackParameters"; 
   if (!SiStripUtility::goToDir(dqm_store, dir)) return;
   
   std::vector<MonitorElement*> meVec1;
   std::vector<MonitorElement*> meVec2;
   if (useGoodTracks_){
     meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties/GoodTracks");
     meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties/GoodTracks");
   }else{
     meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties");
     meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties");
   }
   std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size()); 
   std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
 
   float gstatus = 1.0;
   for (std::vector<MonitorElement*>::const_iterator itME = meVec.begin(); itME != meVec.end(); itME++) {
     MonitorElement * me = (*itME);     
     if (!me) continue;     
     std::vector<QReport *> qt_reports = me->getQReports();          
     if (qt_reports.size() == 0) continue;
     std::string name = me->getName();
 
     float status = 1.0; 
 
     int ibin = 0;
     for (std::map<std::string, TrackingMEs>::const_iterator it = TrackingMEsMap.begin();
          it != TrackingMEsMap.end(); it++) {
       ibin++;
       std::string hname = it->second.HistoName;
       if (name.find(hname) != std::string::npos) {
     status = qt_reports[0]->getQTresult();
     it->second.TrackingFlag->Fill(status);
     fillStatusHistogram(TrackSummaryReportMap, ibin, 1, status);
         break;
       }
     }
     gstatus = gstatus * status; 
   }
   TrackSummaryReportGlobal->Fill(gstatus);
   dqm_store->cd();
 }

void SiStripQualityChecker::fillTrackingStatusAtLumi ( DQMStore * dqm_store )

private

Definition at line 655 of file SiStripQualityChecker.cc.

References DQMStore::cd(), dir, MonitorElement::Fill(), fillStatusHistogram(), DQMStore::getContents(), MonitorElement::getMean(), MonitorElement::getName(), SiStripUtility::goToDir(), relval_steps::merge(), mergeVDriftHistosByStation::name, DQMStore::pwd(), ntuplemaker::status, TrackingMEsMap, TrackSummaryReportGlobal, TrackSummaryReportMap, and useGoodTracks_.

Referenced by fillStatusAtLumi().

                                                                        {
   dqm_store->cd();
   std::string dir = "Tracking"; 
   if (!SiStripUtility::goToDir(dqm_store, dir)) return;
   dir = "TrackParameters"; 
   if (!SiStripUtility::goToDir(dqm_store, dir)) return;
 
   std::vector<MonitorElement*> meVec1;
   std::vector<MonitorElement*> meVec2;
   if (useGoodTracks_){
     meVec1 = dqm_store->getContents(dqm_store->pwd()+"/LSanalysis");
     //    meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties/LSanalysis");
   }else{
     meVec1 = dqm_store->getContents(dqm_store->pwd()+"/GeneralProperties");
     meVec2 = dqm_store->getContents(dqm_store->pwd()+"/HitProperties");
   }
   std::vector<MonitorElement*> meVec(meVec1.size() + meVec2.size()); 
   std::merge(meVec1.begin(), meVec1.end(), meVec2.begin(), meVec2.end(), meVec.begin());
 
   float gstatus = 1.0;
   for (std::vector<MonitorElement*>::const_iterator itME = meVec.begin(); itME != meVec.end(); itME++) {
     MonitorElement * me = (*itME);     
     if (!me) continue;     
     std::string name = me->getName();
 
     float status = -1.0; 
     int ibin = 0;
     for (std::map<std::string, TrackingMEs>::const_iterator it = TrackingMEsMap.begin();
          it != TrackingMEsMap.end(); it++) {
       ibin++;
       std::string hname = it->second.HistoName+"lumiFlag_";
       float lower_cut = it->second.LowerCut; 
       float upper_cut = it->second.UpperCut; 
       if (name.find(hname) != std::string::npos) {
         if (me->getMean() <= lower_cut || me->getMean() > upper_cut) status = 0.0;
         else status = 1.0; 
     it->second.TrackingFlag->Fill(status);
     fillStatusHistogram(TrackSummaryReportMap, ibin, 1, status);
         break;
       } else {
       }
     }
     if (status == -1.0) gstatus = -1.0;
     else gstatus = gstatus * status; 
   }
   TrackSummaryReportGlobal->Fill(gstatus);
   dqm_store->cd();
 }

void SiStripQualityChecker::getModuleStatus	(	DQMStore *	dqm_store,
		std::vector< MonitorElement * > &	layer_mes,
		int &	errdet
	)

private

Definition at line 493 of file SiStripQualityChecker.cc.

References badModuleList, lumiPlot::currentdir, MonitorElement::DQM_KIND_TPROFILE, MonitorElement::DQM_KIND_TPROFILE2D, DQMStore::get(), TkDetMap::getDetFromBin(), SiStripUtility::getMEValue(), MonitorElement::getName(), MonitorElement::getQReports(), MonitorElement::kind(), mergeVDriftHistosByStation::name, DQMStore::pwd(), SiStripUtility::setBadModuleFlag(), and tkDetMap_.

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.size() == 0) 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
     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;
     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 606 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 462 of file SiStripQualityChecker.cc.

References SiStripQualityChecker::SubDetMEs::DetFraction, SiStripUtility::getMEValue(), SiStripQualityChecker::SubDetMEs::SToNFlag, 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.size() > 0) fval1 = atof(sval.c_str());
     SiStripUtility::getMEValue(local_mes.SToNFlag, sval); 
     if (sval.size() > 0) fval2 = atof(sval.c_str());
     SiStripUtility::getMEValue(local_mes.SummaryFlag, sval); 
     if (sval.size() > 0) 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 256 of file SiStripQualityChecker.cc.

References bookedStripStatus_, bookedTrackingStatus_, SiStripQualityChecker::SubDetMEs::DetFraction, DetFractionReportMap, MonitorElement::Reset(), SiStripQualityChecker::SubDetMEs::SToNFlag, SToNReportMap, SubDetMEsMap, SiStripQualityChecker::SubDetMEs::SummaryFlag, SummaryReportGlobal, SummaryReportMap, TrackingMEsMap, TrackSummaryReportGlobal, and TrackSummaryReportMap.

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.SToNFlag->Reset();
       local_mes.SummaryFlag->Reset();
     }
 
     DetFractionReportMap->Reset();
     SToNReportMap->Reset();
     SummaryReportMap->Reset();
 
     SummaryReportGlobal->Reset();
   }
   if (bookedTrackingStatus_) {  
     TrackSummaryReportGlobal->Reset();
     TrackSummaryReportMap->Reset();
     for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsMap.begin();
          it != TrackingMEsMap.end(); it++) {
       it->second.TrackingFlag->Reset();
     }
     /*
     for (std::map<std::string, TrackingMEs>::iterator it = TrackingMEsLSMap.begin();
          it != TrackingMEsLSMap.end(); it++) {
       it->second.TrackingFlag->Reset();
     }
     */
   }
 }