#include <SiStripTrackerMapCreator.h>

Public Types
typedef dqm::harvesting::DQMStore	DQMStore

typedef dqm::harvesting::MonitorElement	MonitorElement

Public Member Functions
void	create (edm::ParameterSet const &tkmapPset, DQMStore &dqm_store, std::string const &htype, edm::EventSetup const &eSetup)

void	createForOffline (edm::ParameterSet const &tkmapPset, DQMStore &dqm_store, std::string &htype, edm::EventSetup const &eSetup)

void	createInfoFile (std::vector< std::string > const &map_names, TTree *tkinfo_tree, DQMStore &dqm_store)

bool	readConfiguration ()

	SiStripTrackerMapCreator (edm::EventSetup const &eSetup)

Private Member Functions
uint16_t	getDetectorFlag (uint32_t const det_id)

uint16_t	getDetectorFlagAndComment (DQMStore dqm_store, uint32_t det_id, TrackerTopology const tTopo, std::ostringstream &comment)

void	paintTkMapFromAlarm (uint32_t det_id, const TrackerTopology *tTopo, DQMStore &dqm_store, bool isBad, std::map< unsigned int, std::string > &badmodmap)

void	paintTkMapFromHistogram (MonitorElement const me, std::string const &map_type, std::vector< std::pair< float, uint32_t >> topNmodVec)

void	printBadModuleList (std::map< unsigned int, std::string > const &badmodmap, edm::EventSetup const &eSetup)

void	printTopModules (std::vector< std::pair< float, uint32_t >> &topNmodVec, edm::EventSetup const &eSetup)

void	setTkMapFromAlarm (DQMStore &dqm_store, edm::EventSetup const &eSetup)

void	setTkMapFromHistogram (DQMStore &dqm_store, std::string const &htype, edm::EventSetup const &eSetup)

void	setTkMapRange (std::string const &map_type)

void	setTkMapRangeOffline ()

Private Attributes
DetId	cachedDetId_ {}

int16_t	cachedLayer_ {}

TkLayerMap::XYbin	cachedXYbin_

edm::ESHandle< SiStripDetCabling >	detCabling_

std::map< uint32_t, uint16_t >	detFlag_

edm::EventSetup const &	eSetup_

float	meanToMaxFactor_ {2.5}

int	nDet_

uint32_t	numTopModules_

bool	ResidualsRMS_

std::string	sRunNumber

std::string	ssqLabel_

std::string	stripTopLevelDir_ {}

TkDetMap const *	tkDetMap_

float	tkMapMax_

float	tkMapMin_

std::string	tkMapName_

std::string	topModLabel_

bool	topModules_

std::unique_ptr< TrackerMap >	trackerMap_ { 0 }

bool	useSSQuality_

Detailed Description

Definition at line 24 of file SiStripTrackerMapCreator.h.

Member Typedef Documentation

typedef dqm::harvesting::DQMStore SiStripTrackerMapCreator::DQMStore

Definition at line 27 of file SiStripTrackerMapCreator.h.

typedef dqm::harvesting::MonitorElement SiStripTrackerMapCreator::MonitorElement

Definition at line 26 of file SiStripTrackerMapCreator.h.

Constructor & Destructor Documentation

SiStripTrackerMapCreator::SiStripTrackerMapCreator ( edm::EventSetup const & eSetup )

Definition at line 23 of file SiStripTrackerMapCreator.cc.

References detCabling_, eSetup_, edm::EventSetup::get(), edm::ESHandle< T >::product(), and tkDetMap_.

                                                                               : eSetup_{eSetup} {
   eSetup_.get<SiStripDetCablingRcd>().get(detCabling_);
   edm::ESHandle<TkDetMap> tkDetMapHandle;
   eSetup_.get<TrackerTopologyRcd>().get(tkDetMapHandle);
   tkDetMap_ = tkDetMapHandle.product();
 }

Member Function Documentation

void SiStripTrackerMapCreator::create	(	edm::ParameterSet const &	tkmapPset,
		DQMStore &	dqm_store,
		std::string const &	htype,
		edm::EventSetup const &	eSetup
	)

Definition at line 33 of file SiStripTrackerMapCreator.cc.

References detCabling_, SiStripDetCabling::fedCabling(), edm::EventSetup::get(), meanToMaxFactor_, nDet_, edm::ESHandle< T >::product(), setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRange(), ssqLabel_, AlCaHLTBitMon_QueryRunRegistry::string, stripTopLevelDir_, tkMapMax_, tkMapMin_, trackerMap_, and useSSQuality_.

                                                                    {
   const SiStripFedCabling* fedcabling = detCabling_->fedCabling();
 
   if (!fedcabling)
     return;
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   trackerMap_ = std::make_unique<TrackerMap>(tkmapPset, fedcabling, tTopo);
   std::string tmap_title = " Tracker Map from  " + map_type;
   trackerMap_->setTitle(tmap_title);
 
   nDet_ = 0;
   tkMapMax_ = 0.0;
   tkMapMin_ = 0.0;
   meanToMaxFactor_ = 2.5;
   useSSQuality_ = false;
   ssqLabel_ = "";
   stripTopLevelDir_ = "SiStrip";
 
   if (map_type == "QTestAlarm") {
     setTkMapFromAlarm(dqm_store, eSetup);
   } else {
     trackerMap_->fill_all_blank();
     setTkMapFromHistogram(dqm_store, map_type, eSetup);
     setTkMapRange(map_type);
   }
   trackerMap_->printonline();
   trackerMap_.reset();
 }

void SiStripTrackerMapCreator::createForOffline	(	edm::ParameterSet const &	tkmapPset,
		DQMStore &	dqm_store,
		std::string &	htype,
		edm::EventSetup const &	eSetup
	)

Definition at line 69 of file SiStripTrackerMapCreator.cc.

References dqm::dqmstoreimpl::DQMStore::cd(), gather_cfg::cout, detCabling_, SiStripDetCabling::fedCabling(), edm::EventSetup::get(), SiStripUtility::goToDir(), meanToMaxFactor_, numTopModules_, edm::ESHandle< T >::product(), dqm::dqmstoreimpl::DQMStore::pwd(), ResidualsRMS_, sistrip::runNumber_, setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRangeOffline(), sRunNumber, contentValuesCheck::ss, ssqLabel_, AlCaHLTBitMon_QueryRunRegistry::string, stripTopLevelDir_, tkMapMax_, tkMapMin_, topModLabel_, topModules_, trackerMap_, and useSSQuality_.

                                                                              {
   // Determine the strip top level dirctory in the DQM file: it is the
   // path where MechanicalView is minus one directory
   std::string const mdir{"MechanicalView"};
   dqm_store.cd();
   if (!SiStripUtility::goToDir(dqm_store, mdir)) {
     edm::LogError("SiStripTopLevelDirNotFound") << "I cannot find the SiStrip top level directory in the DQM file";
   } else {
     const std::string& mechanicalview_dir = dqm_store.pwd();
     stripTopLevelDir_ = mechanicalview_dir.substr(0, mechanicalview_dir.find_last_of("/"));
     edm::LogInfo("SiStripTopLevelDirFound") << "SiStrip top level directory is " << stripTopLevelDir_;
   }
   dqm_store.cd();
 
   //
   const SiStripFedCabling* fedcabling = detCabling_->fedCabling();
 
   if (!fedcabling)
     return;
 
   edm::ESHandle<TrackerTopology> tTopoHandle;
   eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   trackerMap_ = std::make_unique<TrackerMap>(tkmapPset, fedcabling, tTopo);
 
   meanToMaxFactor_ = tkmapPset.getUntrackedParameter<double>("meanToMaxFact", 2.5);
   useSSQuality_ = tkmapPset.getUntrackedParameter<bool>("useSSQuality", false);
   ssqLabel_ = tkmapPset.getUntrackedParameter<std::string>("ssqLabel", "");
   bool tkMapPSU = tkmapPset.getUntrackedParameter<bool>("psuMap", false);
   bool tkMapFED = tkmapPset.getUntrackedParameter<bool>("fedMap", false);
   std::string namesuffix = tkmapPset.getUntrackedParameter<std::string>("mapSuffix", "");
 
   //  std::string tmap_title = " Tracker Map from  " + map_type;
   unsigned int runNumber_ = tkmapPset.getUntrackedParameter<unsigned int>("RunNumber", 1);
   std::stringstream ss;
   ss << runNumber_;
   sRunNumber = ss.str();
   std::cout << sRunNumber << "\n\n\n\n\n\n";
   std::string tmap_title;
   if (runNumber_ > 0) {
     tmap_title = " Run: " + sRunNumber + ", Tracker Map from " + map_type;
   } else {
     tmap_title = " Tracker Map from " + map_type;
   }
   trackerMap_->setTitle(tmap_title);
 
   if (tkmapPset.exists("TopModules"))
     topModules_ = tkmapPset.getUntrackedParameter<bool>("TopModules");
   else
     topModules_ = false;
 
   if (tkmapPset.exists("numberTopModules"))
     numTopModules_ = tkmapPset.getUntrackedParameter<uint32_t>("numberTopModules");
   else
     numTopModules_ = 20;
 
   if (tkmapPset.exists("topModLabel"))
     topModLabel_ = tkmapPset.getUntrackedParameter<int32_t>("topModLabel");
   else
     topModLabel_ = " top Modules " + map_type;
 
   if (map_type == "QTestAlarm") {
     setTkMapFromAlarm(dqm_store, eSetup);
   } else {
     setTkMapFromHistogram(dqm_store, map_type, eSetup);
   }
   // if not overwitten by manual configuration min=0 and max= mean value * meanToMaxFactor_
   setTkMapRangeOffline();
 
   // check manual setting
 
   if (tkmapPset.exists("mapMax"))
     tkMapMax_ = tkmapPset.getUntrackedParameter<double>("mapMax");
   if (tkmapPset.exists("mapMin"))
     tkMapMin_ = tkmapPset.getUntrackedParameter<double>("mapMin");
 
   if (map_type == "ResidualsMean") {
     ResidualsRMS_ = false;
     setTkMapFromHistogram(dqm_store, map_type, eSetup);
     edm::LogInfo("TkMapToBeSaved") << "Ready to save TkMap " << map_type << namesuffix << " with range set to "
                                    << tkMapMin_ << " - " << tkMapMax_;
     trackerMap_->save(true, 0.0000000001, 0.005, map_type + namesuffix + ".svg");
     trackerMap_->save(true, 0.0000000001, 0.005, map_type + namesuffix + ".png", 4500, 2400);
     ResidualsRMS_ = true;
     map_type = "ResidualsRMS";
     if (runNumber_ > 0) {
       tmap_title = " Run: " + sRunNumber + ", Tracker Map from " + map_type;
     } else {
       tmap_title = " Tracker Map from " + map_type;
     }
     trackerMap_->setTitle(tmap_title);
     setTkMapFromHistogram(dqm_store, map_type, eSetup);
     edm::LogInfo("TkMapToBeSaved") << "Ready to save TkMap " << map_type << namesuffix
                                    << " with range set to 0.0 - 1.0";
     trackerMap_->save(true, 0.0, 1.0, map_type + namesuffix + ".svg");
     trackerMap_->save(true, 0.0, 1.0, map_type + namesuffix + ".png", 4500, 2400);
   } else {
     edm::LogInfo("TkMapToBeSaved") << "Ready to save TkMap " << map_type << namesuffix << " with range set to "
                                    << tkMapMin_ << " - " << tkMapMax_;
     trackerMap_->save(true, tkMapMin_, tkMapMax_, map_type + namesuffix + ".svg");
     trackerMap_->save(true, tkMapMin_, tkMapMax_, map_type + namesuffix + ".png", 4500, 2400);
   }
 
   if (tkMapPSU) {
     edm::LogInfo("PSUMapToBeSaved") << "Ready to save PSU TkMap " << map_type << namesuffix << " with range set to "
                                     << tkMapMin_ << " - " << tkMapMax_;
     trackerMap_->save_as_psutrackermap(true, tkMapMin_, tkMapMax_, map_type + namesuffix + "_psu.png");
   }
 
   if (tkMapFED) {
     edm::LogInfo("FEDMapToBeSaved") << "Ready to save FED TkMap " << map_type << namesuffix << " with range set to "
                                     << tkMapMin_ << " - " << tkMapMax_;
     trackerMap_->save_as_fedtrackermap(true, tkMapMin_, tkMapMax_, map_type + namesuffix + "_fed.png");
   }
   trackerMap_.reset();
 }

void SiStripTrackerMapCreator::createInfoFile	(	std::vector< std::string > const &	map_names,
		TTree *	tkinfo_tree,
		DQMStore &	dqm_store
	)

Definition at line 437 of file SiStripTrackerMapCreator.cc.

References dqm::dqmstoreimpl::DQMStore::cd(), BTVHLTOfflineSource_cfi::dirname, eSetup_, edm::EventSetup::get(), getDetectorFlag(), TrackerGeometryUtils::getSiStripDetIds(), SiStripUtility::goToDir(), triggerObjects_cff::id, dqm::dqmstoreimpl::DQMStore::pwd(), AlCaHLTBitMon_QueryRunRegistry::string, and tkDetMap_.

                                                                    {
   std::map<std::string, float> tkhmap_value;
   int qtalarm_flag = 0;
   uint32_t det_id = 0;
 
   if (!tkinfo_tree) {
     edm::LogError("SiStripTrackerMapCreator::createInfoFile") << "Tree not found!";
   } else {
     tkinfo_tree->Branch("DetId", &det_id, "DetId/i");
     for (auto const& mapname : map_names) {
       if (mapname == "QTestAlarm") {
         qtalarm_flag = 0;
         tkinfo_tree->Branch(mapname.c_str(), &qtalarm_flag, std::string(mapname + "/I").c_str());
       } else {
         tkhmap_value[mapname] = -1.0;
         tkinfo_tree->Branch(mapname.c_str(), &tkhmap_value[mapname], std::string(mapname + "/F").c_str());
       }
     }
 
     std::string dirname = "";
 
     std::string mdir = "MechanicalView";
     dqm_store.cd();
     if (!SiStripUtility::goToDir(dqm_store, mdir)) {
       edm::LogError("SiStripTrackerMapCreator::createInfoFile")
           << "I cannot find the SiStrip top level directory in the DQM file";
     } else {
       const std::string& mechanicalview_dir = dqm_store.pwd();
       dirname = mechanicalview_dir.substr(0, mechanicalview_dir.find_last_of("/"));
       edm::LogInfo("SiStripTrackerMapCreator::createInfoFile") << "SiStrip top level directory is " << dirname;
     }
     dqm_store.cd();
 
     std::vector<TkHistoMap> tkHMaps;
 
     uint32_t nHists = map_names.size();
 
     for (uint32_t ih = 0; ih < nHists; ++ih) {
       tkHMaps.emplace_back(tkDetMap_);
       if (map_names.at(ih) != "QTestAlarm") {
         std::string tkhmap_name = "TkHMap_" + map_names.at(ih);
         tkHMaps.at(ih).loadTkHistoMap(dirname, tkhmap_name, true);
       }
     }
 
     edm::ESHandle<GeometricDet> geomDetHandle;
     eSetup_.get<IdealGeometryRecord>().get(geomDetHandle);
     const auto detidList = TrackerGeometryUtils::getSiStripDetIds(*geomDetHandle);
     for (auto const id : detidList) {
       det_id = id;
       for (uint32_t ih = 0; ih < nHists; ++ih) {
         if (map_names.at(ih) == "QTestAlarm") {
           qtalarm_flag = getDetectorFlag(id);
         } else {
           tkhmap_value[map_names.at(ih)] = tkHMaps.at(ih).getValue(id);
         }
       }
       if (!tkinfo_tree) {
         edm::LogError("SiStripTrackerMapCreator::createInfoFile") << "Tree not found!";
       } else {
         tkinfo_tree->Fill();
       }
     }
   }
 }

uint16_t SiStripTrackerMapCreator::getDetectorFlag ( uint32_t const det_id )

inlineprivate

Definition at line 60 of file SiStripTrackerMapCreator.h.

Referenced by createInfoFile().

                                                   {
     return detFlag_.find(det_id) != detFlag_.end() ? detFlag_[det_id] : 0;
   }

uint16_t SiStripTrackerMapCreator::getDetectorFlagAndComment	(	DQMStore *	dqm_store,
		uint32_t	det_id,
		TrackerTopology const *	tTopo,
		std::ostringstream &	comment
	)

private

Definition at line 369 of file SiStripTrackerMapCreator.cc.

References dqm::dqmstoreimpl::DQMStore::cd(), SiStripCommissioningClient_cfg::conn, constexpr, detCabling_, dqm::dqmstoreimpl::DQMStore::dirExists(), RemoveAddSevLevel::flag, dqm::dqmstoreimpl::DQMStore::get(), SiStripUtility::getBadModuleStatus(), SiStripDetCabling::getConnections(), SiStripFolderOrganizer::getSubDetFolder(), MonitorElementData::INT, LogDebug, SiStripFolderOrganizer::setSiStripFolderName(), AlCaHLTBitMon_QueryRunRegistry::string, and stripTopLevelDir_.

Referenced by paintTkMapFromAlarm(), and printTopModules().

                                                                                         {
   comment << "Module " << det_id;
   uint16_t flag = 0;
 
   // get FED channels corresponding to the det_id
 
   comment << " FEDCHs ";
   constexpr std::size_t indent_width{7ull};
   std::vector<const FedChannelConnection*> conns = detCabling_->getConnections(det_id);
   for (auto const conn : conns) {
     if (conn) {
       comment << std::setw(3) << conn->fedId() << "/" << std::setw(2) << conn->fedCh() << " ";
     } else {
       comment << std::string(indent_width, ' ');
     }
   }
   if (conns.empty()) {
     comment << std::string(indent_width * 3, ' ');
   } else if (conns.size() == 1) {
     comment << std::string(indent_width * 2, ' ');
   } else if (conns.size() == 2) {
     comment << std::string(indent_width, ' ');
   }
 
   if (!dqm_store)
     return flag;
 
   SiStripFolderOrganizer folder_organizer;
   std::string subdet_folder, badmodule_folder;
 
   dqm_store->cd();
 
   folder_organizer.setSiStripFolderName(stripTopLevelDir_);
   folder_organizer.getSubDetFolder(det_id, tTopo, subdet_folder);
 
   LogDebug("SearchBadModule") << det_id << " " << subdet_folder << " " << stripTopLevelDir_;
 
   if (dqm_store->dirExists(subdet_folder)) {
     badmodule_folder = subdet_folder + "/BadModuleList";
     LogDebug("SearchBadModule") << subdet_folder << " exists: " << badmodule_folder;
   } else {
     edm::LogError("SubDetFolderNotFound") << subdet_folder << " does not exist for detid " << det_id;
     return flag;
   }
   if (!dqm_store->dirExists(badmodule_folder)) {
     LogDebug("BadModuleFolderNotFound") << badmodule_folder << " does not exist for detid " << det_id;
     return flag;
   }
   std::ostringstream badmodule_path;
   badmodule_path << badmodule_folder << "/" << det_id;
   LogDebug("SearchBadModule") << badmodule_folder << " exists: " << badmodule_path;
 
   auto const* bad_module_me = dqm_store->get(badmodule_path.str());
   if (bad_module_me && bad_module_me->kind() == MonitorElement::Kind::INT) {
     LogDebug("SearchBadModule") << "Monitor Element found";
     flag = bad_module_me->getIntValue();
     std::string message;
     SiStripUtility::getBadModuleStatus(flag, message);
     comment << message.c_str();
   }
   return flag;
 }

void SiStripTrackerMapCreator::paintTkMapFromAlarm	(	uint32_t	det_id,
		const TrackerTopology *	tTopo,
		DQMStore &	dqm_store,
		bool	isBad,
		std::map< unsigned int, std::string > &	badmodmap
	)

private

Definition at line 235 of file SiStripTrackerMapCreator.cc.

References comment, detFlag_, RemoveAddSevLevel::flag, getDetectorFlagAndComment(), SiStripUtility::getDetectorStatusColor(), TrackerTopology::layer(), AlCaHLTBitMon_QueryRunRegistry::string, trackerMap_, and parallelization::uint.

Referenced by setTkMapFromAlarm().

                                                                                                {
   std::ostringstream comment;
   uint16_t flag = 0;
   flag = getDetectorFlagAndComment(&dqm_store, det_id, tTopo, comment);
 
   int rval, gval, bval;
   SiStripUtility::getDetectorStatusColor(flag, rval, gval, bval);
   if (isBad) {
     rval = 255;
     gval = 255;
     bval = 0;
     comment << " PCLBadModule ";
   }
   trackerMap_->setText(det_id, comment.str());
   trackerMap_->fillc(det_id, rval, gval, bval);
 
   if (flag != 0 || isBad) {
     uint lay = tTopo->layer(SiStripDetId(det_id));
     std::string layer = " Layer " + std::to_string(lay);
     badmodmap[det_id] = comment.str() + layer;
   }
 
   detFlag_[det_id] = flag;
 }

void SiStripTrackerMapCreator::paintTkMapFromHistogram	(	MonitorElement const *	me,
		std::string const &	map_type,
		std::vector< std::pair< float, uint32_t >> *	topNmodVec
	)

private

Definition at line 325 of file SiStripTrackerMapCreator.cc.

References funct::abs(), cachedDetId_, cachedLayer_, cachedXYbin_, dqm::impl::MonitorElement::getBinContent(), TkDetMap::getDetsForLayer(), TkDetMap::getLayerNum(), dqm::impl::MonitorElement::getName(), dqm::legacy::MonitorElement::getTProfile2D(), TkDetMap::getXY(), TkLayerMap::XYbin::ix, TkLayerMap::XYbin::iy, dqm::impl::MonitorElement::kind(), Skims_PA_cff::name, nDet_, ResidualsRMS_, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, tkDetMap_, tkMapMax_, MonitorElementData::TPROFILE2D, and trackerMap_.

Referenced by setTkMapFromHistogram().

                                                                                                         {
   const std::string& name = me->getName();
   std::string lname = name.substr(name.find("TkHMap_") + 7);
   lname = lname.substr(lname.find("_T") + 1);
   for (DetId det_id : tkDetMap_->getDetsForLayer(TkDetMap::getLayerNum(lname))) {
     if (det_id.rawId() <= 0)
       continue;
     nDet_++;
     const TkLayerMap::XYbin& xyval = tkDetMap_->getXY(det_id, cachedDetId_, cachedLayer_, cachedXYbin_);
     float fval = 0.0;
     if (name.find("Residuals") != std::string::npos) {
       if (ResidualsRMS_ == true) {
         if (me->kind() == MonitorElement::Kind::TPROFILE2D) {
           TProfile2D* tp = me->getTProfile2D();
           float fval_prov =
               tp->GetBinError(xyval.ix, xyval.iy) * sqrt(tp->GetBinEntries(tp->GetBin(xyval.ix, xyval.iy)));
           fval = fval_prov;
         }
       } else {
         float fval_prov = me->getBinContent(xyval.ix, xyval.iy);
         fval = std::abs(fval_prov);
       }
     } else
       fval = me->getBinContent(xyval.ix, xyval.iy);
     if (htype == "QTestAlarm") {
       edm::LogError("ItShouldNotBeHere") << "QTestAlarm map: you should not be here!";
     } else {
       if (fval == 0.0)
         trackerMap_->fillc(det_id, 255, 255, 255);
       else {
         trackerMap_->fill_current_val(det_id, fval);
         if (topNmodVec) {
           topNmodVec->emplace_back(fval, det_id);
         }
       }
       tkMapMax_ += fval;
     }
   }
 }

void SiStripTrackerMapCreator::printBadModuleList	(	std::map< unsigned int, std::string > const &	badmodmap,
		edm::EventSetup const &	eSetup
	)

private

Definition at line 507 of file SiStripTrackerMapCreator.cc.

References constexpr, edm::EventSetup::get(), edm::ESHandle< T >::product(), sRunNumber, SiStripDetId::subDetector(), SiStripDetId::TEC, TrackerTopology::tecSide(), SiStripDetId::TIB, SiStripDetId::TID, TrackerTopology::tidSide(), and SiStripDetId::TOB.

Referenced by setTkMapFromAlarm().

                                                                                {
   // Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   bool tibDone = false, tidSide1Done = false, tidSide2Done = false, tobDone = false, tecSide1Done = false,
        tecSide2Done = false;
   constexpr unsigned int tibFirst = 369120277 - 1, tidSide1First = 402664197 - 1, tidSide2First = 402672389 - 1,
                          tobFirst = 436228134 - 1, tecSide1First = 470044965 - 1, tecSide2First = 470307109 - 1;
 
   int ntib = 0, ntids1 = 0, ntids2 = 0, ntob = 0, ntecs1 = 0, ntecs2 = 0;
 
   for (auto const& badmod : badmodmap) {
     SiStripDetId ssdetid(badmod.first);
     if (ssdetid.subDetector() == SiStripDetId::TIB)
       ntib++;
     if (ssdetid.subDetector() == SiStripDetId::TID) {
       if (tTopo->tidSide(ssdetid) == 1)
         ntids1++;
       if (tTopo->tidSide(ssdetid) == 2)
         ntids2++;
     }
     if (ssdetid.subDetector() == SiStripDetId::TOB)
       ntob++;
     if (ssdetid.subDetector() == SiStripDetId::TEC) {
       if (tTopo->tecSide(ssdetid) == 1)
         ntecs1++;
       if (tTopo->tecSide(ssdetid) == 2)
         ntecs2++;
     }
   }
 
   edm::LogVerbatim("BadModuleList") << "Run: " << sRunNumber << ", Number of bad modules in total:";
   edm::LogVerbatim("BadModuleList") << "--------------------------------------------------------------";
   edm::LogVerbatim("BadModuleList") << "TIB: " << ntib;
   edm::LogVerbatim("BadModuleList") << "TID/MINUS: " << ntids1;
   edm::LogVerbatim("BadModuleList") << "TID/PLUS: " << ntids2;
   edm::LogVerbatim("BadModuleList") << "TOB: " << ntob;
   edm::LogVerbatim("BadModuleList") << "TEC/MINUS: " << ntecs1;
   edm::LogVerbatim("BadModuleList") << "TEC/PLUS: " << ntecs2;
   edm::LogVerbatim("BadModuleList") << "-------------------------------";
   edm::LogVerbatim("BadModuleList");
   edm::LogVerbatim("BadModuleList") << "List of bad modules per partition:";
   edm::LogVerbatim("BadModuleList") << "----------------------------------";
 
   for (auto const& badmod : badmodmap) {
     if (!tibDone && badmod.first >= tibFirst) {
       tibDone = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TIB";
       edm::LogVerbatim("BadModuleList");
     }
     if (!tidSide1Done && badmod.first >= tidSide1First) {
       tidSide1Done = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TID/MINUS";
       edm::LogVerbatim("BadModuleList");
     }
     if (!tidSide2Done && badmod.first >= tidSide2First) {
       tidSide2Done = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TID/PLUS";
       edm::LogVerbatim("BadModuleList");
     }
     if (!tobDone && badmod.first >= tobFirst) {
       tobDone = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TOB";
       edm::LogVerbatim("BadModuleList");
     }
     if (!tecSide1Done && badmod.first >= tecSide1First) {
       tecSide1Done = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TEC/MINUS";
       edm::LogVerbatim("BadModuleList");
     }
     if (!tecSide2Done && badmod.first >= tecSide2First) {
       tecSide2Done = true;
       edm::LogVerbatim("BadModuleList");
       edm::LogVerbatim("BadModuleList") << "SubDetector TEC/PLUS";
       edm::LogVerbatim("BadModuleList");
     }
     edm::LogVerbatim("BadModuleList") << badmod.second;
   }
 }

void SiStripTrackerMapCreator::printTopModules	(	std::vector< std::pair< float, uint32_t >> &	topNmodVec,
		edm::EventSetup const &	eSetup
	)

private

Definition at line 595 of file SiStripTrackerMapCreator.cc.

References comment, RemoveAddSevLevel::flag, edm::EventSetup::get(), getDetectorFlagAndComment(), numTopModules_, edm::ESHandle< T >::product(), AlCaHLTBitMon_QueryRunRegistry::string, SiStripDetId::subDetector(), subdetector, SiStripDetId::TEC, TrackerTopology::tecSide(), SiStripDetId::TIB, SiStripDetId::TID, TrackerTopology::tidSide(), SiStripDetId::TOB, and topModLabel_.

Referenced by setTkMapFromHistogram().

                                                                             {
   edm::ESHandle<TrackerTopology> tTopoHandle;
   eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   if (topNmodVec.empty())
     return;
 
   std::sort(topNmodVec.rbegin(), topNmodVec.rend());
   if (topNmodVec.size() > numTopModules_)
     topNmodVec.resize(numTopModules_);
 
   edm::LogVerbatim("TopModules") << topModLabel_;
   edm::LogVerbatim("TopModules") << "------------------------------------------------------";
 
   for (auto const& aPair : topNmodVec) {
     uint32_t det_id = aPair.second;
     std::ostringstream comment;
     std::string subdetector;
     SiStripDetId ssdetid(aPair.second);
     if (ssdetid.subDetector() == SiStripDetId::TIB)
       subdetector = "TIB       ";
     if (ssdetid.subDetector() == SiStripDetId::TID) {
       if (tTopo->tidSide(ssdetid) == 1)
         subdetector = "TID/MINUS ";
       if (tTopo->tidSide(ssdetid) == 2)
         subdetector = "TID/PLUS  ";
     }
     if (ssdetid.subDetector() == SiStripDetId::TOB)
       subdetector = "TOB       ";
     if (ssdetid.subDetector() == SiStripDetId::TEC) {
       if (tTopo->tecSide(ssdetid) == 1)
         subdetector = "TEC/MINUS ";
       if (tTopo->tecSide(ssdetid) == 2)
         subdetector = "TEC/PLUS  ";
     }
     uint16_t flag = getDetectorFlagAndComment(nullptr, det_id, tTopo, comment);
     if (flag == 0)
       edm::LogVerbatim("TopModules") << subdetector << comment.str() << " value: " << aPair.first;
   }
   edm::LogVerbatim("TopModules") << "------------------------------------------------------";
 }

bool SiStripTrackerMapCreator::readConfiguration ( )

void SiStripTrackerMapCreator::setTkMapFromAlarm	(	DQMStore &	dqm_store,
		edm::EventSetup const &	eSetup
	)

private

Definition at line 193 of file SiStripTrackerMapCreator.cc.

References detCabling_, eSetup_, edm::EventSetup::get(), SiStripDetCabling::getDetCabling(), SiStripQuality::IsModuleBad(), nDet_, paintTkMapFromAlarm(), printBadModuleList(), edm::ESHandle< T >::product(), ssqLabel_, tkMapMax_, tkMapMin_, trackerMap_, and useSSQuality_.

Referenced by create(), and createForOffline().

                                                                                                  {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   eSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   nDet_ = 0;
   tkMapMax_ = 0.0;
   tkMapMin_ = 0.0;
 
   edm::ESHandle<SiStripQuality> ssq;
 
   if (useSSQuality_) {
     eSetup_.get<SiStripQualityRcd>().get(ssqLabel_, ssq);
   }
 
   trackerMap_->fillc_all_blank();
 
   std::map<unsigned int, std::string> badmodmap;
 
   // used to avoid multiple checks on the same detid since the loop is done on the FED channels
   uint32_t detId_save = 0;
   // example of loop using SiStripDetCabling
   for (auto const& module : detCabling_->getDetCabling()) {
     uint32_t detId = module.first;
     if (detId == 0 || detId == 0xFFFFFFFF)
       continue;
     if (detId_save != detId) {
       detId_save = detId;
       bool isBad = useSSQuality_ && ssq->IsModuleBad(detId);
       paintTkMapFromAlarm(detId, tTopo, dqm_store, isBad, badmodmap);
     } else {
       edm::LogWarning("TwiceTheSameDetId")
           << "The detid " << detId << " was found already in the loop on SiStripDetCabling";
     }
   }
   printBadModuleList(badmodmap, eSetup);
 }

void SiStripTrackerMapCreator::setTkMapFromHistogram	(	DQMStore &	dqm_store,
		std::string const &	htype,
		edm::EventSetup const &	eSetup
	)

private

Definition at line 266 of file SiStripTrackerMapCreator.cc.

References dqm::dqmstoreimpl::DQMStore::cd(), dqm::dqmstoreimpl::DQMStore::dirExists(), RefreshWebPage::dname, dqm::dqmstoreimpl::DQMStore::getContents(), dqm::dqmstoreimpl::DQMStore::getSubdirs(), SiStripUtility::goToDir(), Skims_PA_cff::name, nDet_, paintTkMapFromHistogram(), printTopModules(), dqm::dqmstoreimpl::DQMStore::pwd(), sd, AlCaHLTBitMon_QueryRunRegistry::string, tkMapMax_, tkMapMin_, and topModules_.

Referenced by create(), and createForOffline().

                                                                                   {
   dqm_store.cd();
 
   std::string const mdir{"MechanicalView"};
   if (!SiStripUtility::goToDir(dqm_store, mdir))
     return;
   std::string mechanicalview_dir = dqm_store.pwd();
 
   std::vector<std::string> subdet_folder;
   subdet_folder.push_back("TIB");
   subdet_folder.push_back("TOB");
   subdet_folder.push_back("TEC/MINUS");
   subdet_folder.push_back("TEC/PLUS");
   subdet_folder.push_back("TID/MINUS");
   subdet_folder.push_back("TID/PLUS");
 
   nDet_ = 0;
   tkMapMax_ = 0.0;
   tkMapMin_ = 0.0;
 
   auto topNmodVec = topModules_ ? new std::vector<std::pair<float, uint32_t>>{} : nullptr;
   for (auto const& sd : subdet_folder) {
     std::string dname = mechanicalview_dir + "/" + sd;
     if (!dqm_store.dirExists(dname))
       continue;
     dqm_store.cd(dname);
     std::vector<std::string> layerVec = dqm_store.getSubdirs();
     for (auto const& layer : layerVec) {
       if (layer.find("BadModuleList") != std::string::npos)
         continue;
       std::vector<MonitorElement*> meVec = dqm_store.getContents(layer);
       MonitorElement* tkhmap_me{nullptr};
       std::string name;
       for (auto tkh : meVec) {
         name = tkh->getName();
         if (name.find("TkHMap") == std::string::npos)
           continue;
         if (htype == "QTestAlarm") {
           edm::LogError("ItShouldNotBeHere") << "QTestAlarm map: you should not be here!";
           tkhmap_me = tkh;
           break;
         } else if (name.find(htype) != std::string::npos) {
           tkhmap_me = tkh;
           break;
         }
       }
       if (tkhmap_me != nullptr) {
         paintTkMapFromHistogram(tkhmap_me, htype, topNmodVec);
       }
     }
     dqm_store.cd(mechanicalview_dir);
   }
   dqm_store.cd();
   if (topNmodVec)
     printTopModules(*topNmodVec, eSetup);
 }

void SiStripTrackerMapCreator::setTkMapRange ( std::string const & map_type )

private

Definition at line 643 of file SiStripTrackerMapCreator.cc.

References meanToMaxFactor_, nDet_, tkMapMax_, tkMapMin_, and trackerMap_.

Referenced by create().

                                                                       {
   tkMapMin_ = 0.0;
   if (tkMapMax_ == 0.0) {
     if (map_type.find("FractionOfBadChannels") != std::string::npos)
       tkMapMax_ = 1.0;
     else if (map_type.find("NumberOfCluster") != std::string::npos)
       tkMapMax_ = 0.01;
     else if (map_type.find("NumberOfDigi") != std::string::npos)
       tkMapMax_ = 0.6;
     else if (map_type.find("NumberOfOffTrackCluster") != std::string::npos)
       tkMapMax_ = 100.0;
     else if (map_type.find("NumberOfOnTrackCluster") != std::string::npos)
       tkMapMax_ = 50.0;
     else if (map_type.find("StoNCorrOnTrack") != std::string::npos)
       tkMapMax_ = 200.0;
   } else {
     tkMapMax_ = tkMapMax_ / nDet_ * 1.0;
     tkMapMax_ = tkMapMax_ * meanToMaxFactor_;
   }
   trackerMap_->setRange(tkMapMin_, tkMapMax_);
 }

void SiStripTrackerMapCreator::setTkMapRangeOffline ( )

private

Definition at line 665 of file SiStripTrackerMapCreator.cc.

References meanToMaxFactor_, nDet_, tkMapMax_, and tkMapMin_.

Referenced by createForOffline().

                                                     {
   tkMapMin_ = 0.0;
   if (tkMapMax_ != 0.0) {
     tkMapMax_ = tkMapMax_ / (nDet_ * 1.0);
     tkMapMax_ = tkMapMax_ * meanToMaxFactor_;
   }
 }