SiStripTrackerMapCreator Class Reference

#include <SiStripTrackerMapCreator.h>

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, std::vector< uint32_t > const &detidList)
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_ {nullptr}
bool	useSSQuality_

Detailed Description

Definition at line 24 of file SiStripTrackerMapCreator.h.

Constructor & Destructor Documentation

SiStripTrackerMapCreator::SiStripTrackerMapCreator

(

edm::EventSetup const &

eSetup

)

Definition at line 22 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 36 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 73 of file SiStripTrackerMapCreator.cc.

References DQMStore::cd(), gather_cfg::cout, detCabling_, SiStripDetCabling::fedCabling(), edm::EventSetup::get(), SiStripUtility::goToDir(), meanToMaxFactor_, numTopModules_, edm::ESHandle< T >::product(), DQMStore::pwd(), ResidualsRMS_, sistrip::runNumber_, setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRangeOffline(), sRunNumber, 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,
		std::vector< uint32_t > const &	detidList
	)

Definition at line 450 of file SiStripTrackerMapCreator.cc.

References DQMStore::cd(), compare_using_db::dirname, getDetectorFlag(), SiStripUtility::goToDir(), triggerObjects_cff::id, 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);
      }
    }

    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 64 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 377 of file SiStripTrackerMapCreator.cc.

References DQMStore::cd(), constexpr, detCabling_, DQMStore::dirExists(), MonitorElement::DQM_KIND_INT, RemoveAddSevLevel::flag, DQMStore::get(), SiStripUtility::getBadModuleStatus(), SiStripDetCabling::getConnections(), SiStripFolderOrganizer::getSubDetFolder(), 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::DQM_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 242 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 333 of file SiStripTrackerMapCreator.cc.

References funct::abs(), cachedDetId_, cachedLayer_, cachedXYbin_, MonitorElement::DQM_KIND_TPROFILE2D, MonitorElement::getBinContent(), TkDetMap::getDetsForLayer(), TkDetMap::getLayerNum(), MonitorElement::getName(), MonitorElement::getTProfile2D(), TkDetMap::getXY(), TkLayerMap::XYbin::ix, TkLayerMap::XYbin::iy, MonitorElement::kind(), dataset::name, nDet_, ResidualsRMS_, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, tkDetMap_, tkMapMax_, 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::DQM_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 523 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 616 of file SiStripTrackerMapCreator.cc.

References comment, RemoveAddSevLevel::flag, edm::EventSetup::get(), getDetectorFlagAndComment(), numTopModules_, edm::ESHandle< T >::product(), jetUpdater_cfi::sort, 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 199 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 276 of file SiStripTrackerMapCreator.cc.

References DQMStore::cd(), DQMStore::dirExists(), RefreshWebPage::dname, DQMStore::getContents(), DQMStore::getSubdirs(), SiStripUtility::goToDir(), dataset::name, nDet_, paintTkMapFromHistogram(), printTopModules(), 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 662 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 686 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_;
  }
}

Member Data Documentation

DetId SiStripTrackerMapCreator::cachedDetId_ {}

private

Definition at line 88 of file SiStripTrackerMapCreator.h.

Referenced by paintTkMapFromHistogram().

int16_t SiStripTrackerMapCreator::cachedLayer_ {}

private

Definition at line 89 of file SiStripTrackerMapCreator.h.

Referenced by paintTkMapFromHistogram().

TkLayerMap::XYbin SiStripTrackerMapCreator::cachedXYbin_

private

Definition at line 91 of file SiStripTrackerMapCreator.h.

Referenced by paintTkMapFromHistogram().

edm::ESHandle<SiStripDetCabling> SiStripTrackerMapCreator::detCabling_

private

Definition at line 87 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), getDetectorFlagAndComment(), setTkMapFromAlarm(), and SiStripTrackerMapCreator().

std::map<uint32_t, uint16_t> SiStripTrackerMapCreator::detFlag_

private

Definition at line 90 of file SiStripTrackerMapCreator.h.

Referenced by paintTkMapFromAlarm().

edm::EventSetup const& SiStripTrackerMapCreator::eSetup_

private

Definition at line 86 of file SiStripTrackerMapCreator.h.

Referenced by setTkMapFromAlarm(), and SiStripTrackerMapCreator().

float SiStripTrackerMapCreator::meanToMaxFactor_ {2.5}

private

Definition at line 80 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), setTkMapRange(), and setTkMapRangeOffline().

int SiStripTrackerMapCreator::nDet_

private

Definition at line 84 of file SiStripTrackerMapCreator.h.

Referenced by create(), paintTkMapFromHistogram(), setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRange(), and setTkMapRangeOffline().

uint32_t SiStripTrackerMapCreator::numTopModules_

private

Definition at line 93 of file SiStripTrackerMapCreator.h.

Referenced by createForOffline(), and printTopModules().

bool SiStripTrackerMapCreator::ResidualsRMS_

private

Definition at line 82 of file SiStripTrackerMapCreator.h.

Referenced by createForOffline(), and paintTkMapFromHistogram().

std::string SiStripTrackerMapCreator::sRunNumber

private

Definition at line 74 of file SiStripTrackerMapCreator.h.

Referenced by createForOffline(), and printBadModuleList().

std::string SiStripTrackerMapCreator::ssqLabel_

private

Definition at line 83 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), and setTkMapFromAlarm().

std::string SiStripTrackerMapCreator::stripTopLevelDir_ {}

private

Definition at line 76 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), and getDetectorFlagAndComment().

TkDetMap const* SiStripTrackerMapCreator::tkDetMap_

private

Definition at line 85 of file SiStripTrackerMapCreator.h.

Referenced by createInfoFile(), paintTkMapFromHistogram(), and SiStripTrackerMapCreator().

float SiStripTrackerMapCreator::tkMapMax_

private

Definition at line 78 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), paintTkMapFromHistogram(), setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRange(), and setTkMapRangeOffline().

float SiStripTrackerMapCreator::tkMapMin_

private

Definition at line 79 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), setTkMapFromAlarm(), setTkMapFromHistogram(), setTkMapRange(), and setTkMapRangeOffline().

std::string SiStripTrackerMapCreator::tkMapName_

private

Definition at line 75 of file SiStripTrackerMapCreator.h.

std::string SiStripTrackerMapCreator::topModLabel_

private

Definition at line 94 of file SiStripTrackerMapCreator.h.

Referenced by createForOffline(), and printTopModules().

bool SiStripTrackerMapCreator::topModules_

private

Definition at line 92 of file SiStripTrackerMapCreator.h.

Referenced by createForOffline(), and setTkMapFromHistogram().

std::unique_ptr<TrackerMap> SiStripTrackerMapCreator::trackerMap_ {nullptr}

private

Definition at line 73 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), paintTkMapFromAlarm(), paintTkMapFromHistogram(), setTkMapFromAlarm(), and setTkMapRange().

bool SiStripTrackerMapCreator::useSSQuality_

private

Definition at line 81 of file SiStripTrackerMapCreator.h.

Referenced by create(), createForOffline(), and setTkMapFromAlarm().