SiStripDetVOffTkMapPlotter.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include <iostream>
#include <sstream>
 
#include "CondCore/CondDB/interface/ConnectionPool.h"
 
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "Geometry/TrackerNumberingBuilder/interface/utils.h"
 
#include "CondFormats/SiStripObjects/interface/SiStripDetVOff.h"
#include "CondFormats/Common/interface/Time.h"
#include "CondFormats/Common/interface/TimeConversions.h"
 
#include "DQM/SiStripCommon/interface/TkHistoMap.h"
#include "CommonTools/TrackerMap/interface/TrackerMap.h"
 
class SiStripDetVOffTkMapPlotter : public edm::EDAnalyzer {
public:
  explicit SiStripDetVOffTkMapPlotter(const edm::ParameterSet& iConfig);
  ~SiStripDetVOffTkMapPlotter() override;
  void analyze(const edm::Event& evt, const edm::EventSetup& evtSetup) override;
  void endJob() override;
 
private:
  std::string formatIOV(cond::Time_t iov, std::string format = "%Y-%m-%d__%H_%M_%S");
 
  cond::persistency::ConnectionPool m_connectionPool;
  std::string m_condDb;
  std::string m_plotTag;
 
  // IOV of plotting.
  cond::Time_t m_IOV;
  // Or use datatime string. Format: "2002-01-20 23:59:59.000". Set IOV to 0 to use this.
  std::string m_Time;
  // Set the plot format. Default: png.
  std::string m_plotFormat;
  // Specify output root file name. Leave empty if do not want to save plots in a root file.
  std::string m_outputFile;
};
 
SiStripDetVOffTkMapPlotter::SiStripDetVOffTkMapPlotter(const edm::ParameterSet& iConfig)
    : m_connectionPool(),
      m_condDb(iConfig.getParameter<std::string>("conditionDatabase")),
      m_plotTag(iConfig.getParameter<std::string>("Tag")),
      m_IOV(iConfig.getUntrackedParameter<cond::Time_t>("IOV", 0)),
      m_Time(iConfig.getUntrackedParameter<std::string>("Time", "")),
      m_plotFormat(iConfig.getUntrackedParameter<std::string>("plotFormat", "png")),
      m_outputFile(iConfig.getUntrackedParameter<std::string>("outputFile", "")) {
  m_connectionPool.setParameters(iConfig.getParameter<edm::ParameterSet>("DBParameters"));
  m_connectionPool.configure();
}
 
SiStripDetVOffTkMapPlotter::~SiStripDetVOffTkMapPlotter() {}
 
void SiStripDetVOffTkMapPlotter::analyze(const edm::Event& evt, const edm::EventSetup& evtSetup) {
  cond::Time_t theIov = 0;
  if (m_IOV != 0) {
    theIov = m_IOV;
  } else if (!m_Time.empty()) {
    theIov = cond::time::from_boost(boost::posix_time::time_from_string(m_Time));
  } else {
    // Use the current time if no input. Will get the last IOV.
    theIov = cond::time::from_boost(boost::posix_time::second_clock::universal_time());
  }
 
  // open db session
  edm::LogInfo("SiStripDetVOffMapPlotter") << "[SiStripDetVOffMapPlotter::" << __func__ << "] "
                                           << "Query the condition database " << m_condDb << " for tag " << m_plotTag;
  cond::persistency::Session condDbSession = m_connectionPool.createSession(m_condDb);
  condDbSession.transaction().start(true);
  cond::persistency::IOVProxy iovProxy = condDbSession.readIov(m_plotTag);
  auto iovs = iovProxy.selectAll();
  auto iiov = iovs.find(theIov);
  if (iiov == iovs.end())
    throw cms::Exception("Input IOV " + std::to_string(m_IOV) + "/" + m_Time + " is invalid!");
 
  theIov = (*iiov).since;
  edm::LogInfo("SiStripDetVOffMapPlotter") << "[SiStripDetVOffMapPlotter::" << __func__ << "] "
                                           << "Make tkMap for IOV " << theIov << " ("
                                           << boost::posix_time::to_simple_string(cond::time::to_boost(theIov)) << ")";
  auto payload = condDbSession.fetchPayload<SiStripDetVOff>((*iiov).payloadId);
 
  edm::ESHandle<TkDetMap> tkDetMapHandle;
  evtSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
  const TkDetMap* tkDetMap = tkDetMapHandle.product();
  TrackerMap lvmap, hvmap;
  TkHistoMap lvhisto(tkDetMap, "LV_Status", "LV_Status", -1);
  TkHistoMap hvhisto(tkDetMap, "HV_Status", "HV_Status", -1);
 
  edm::ESHandle<GeometricDet> geomDetHandle;
  evtSetup.get<IdealGeometryRecord>().get(geomDetHandle);
 
  const auto detids = TrackerGeometryUtils::getSiStripDetIds(*geomDetHandle);
  for (auto id : detids) {
    if (payload->IsModuleLVOff(id))
      lvhisto.fill(id, 1);  // RED
    else
      lvhisto.fill(id, 0.5);
 
    if (payload->IsModuleHVOff(id))
      hvhisto.fill(id, 1);  // RED
    else
      hvhisto.fill(id, 0.5);
  }
 
  lvhisto.dumpInTkMap(&lvmap);
  hvhisto.dumpInTkMap(&hvmap);
  lvmap.setPalette(1);
  hvmap.setPalette(1);
  lvmap.save(true, 0, 0, "LV_tkMap_" + formatIOV(theIov) + "." + m_plotFormat);
  hvmap.save(true, 0, 0, "HV_tkMap_" + formatIOV(theIov) + "." + m_plotFormat);
 
  if (!m_outputFile.empty()) {
    lvhisto.save(m_outputFile);
    hvhisto.save(m_outputFile);
  }
}
 
void SiStripDetVOffTkMapPlotter::endJob() {}
 
std::string SiStripDetVOffTkMapPlotter::formatIOV(cond::Time_t iov, std::string format) {
  auto facet = new boost::posix_time::time_facet(format.c_str());
  std::ostringstream stream;
  stream.imbue(std::locale(stream.getloc(), facet));
  stream << cond::time::to_boost(iov);
  return stream.str();
}
 
DEFINE_FWK_MODULE(SiStripDetVOffTkMapPlotter);