Inheritance diagram for SiStripDetVOffTrendPlotter:

Public Member Functions
void	analyze (const edm::Event &evt, const edm::EventSetup &evtSetup) override

	SiStripDetVOffTrendPlotter (const edm::ParameterSet &iConfig)

	~SiStripDetVOffTrendPlotter () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Public Attributes
const std::vector< Color_t >	PLOT_COLORS {kRed, kBlue, kBlack, kOrange, kMagenta}

Private Member Functions
void	dumpCSV (bool isHV, std::size_t nModules)

std::string	formatIOV (cond::Time_t iov, std::string format="%Y-%m-%d__%H_%M_%S")

void	prepGraph (TGraph *gr, TString name, TString title, Color_t color)

Private Attributes
TFile *	fout

edm::ESGetToken< GeometricDet, IdealGeometryRecord >	geomDetToken_

std::map< cond::Time_t, std::map< std::string, std::pair< int, int > > >	iovMap

std::string	m_condDb

cond::persistency::ConnectionPool	m_connectionPool

std::string	m_endTime

int	m_interval

std::string	m_outputCSV

std::string	m_outputPlot

std::string	m_outputRootFile

std::vector< std::string >	m_plotTags

std::string	m_startTime

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 29 of file SiStripDetVOffTrendPlotter.cc.

Constructor & Destructor Documentation

◆ SiStripDetVOffTrendPlotter()

SiStripDetVOffTrendPlotter::SiStripDetVOffTrendPlotter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 65 of file SiStripDetVOffTrendPlotter.cc.

References cond::persistency::ConnectionPool::configure(), fout, edm::ParameterSet::getParameter(), m_connectionPool, m_outputRootFile, and cond::persistency::ConnectionPool::setParameters().

     : m_connectionPool(),
       m_condDb(iConfig.getParameter<std::string>("conditionDatabase")),
       m_plotTags(iConfig.getParameter<std::vector<std::string>>("plotTags")),
       m_interval(iConfig.getParameter<int>("timeInterval")),
       m_startTime(iConfig.getUntrackedParameter<std::string>("startTime", "")),
       m_endTime(iConfig.getUntrackedParameter<std::string>("endTime", "")),
       m_outputPlot(iConfig.getUntrackedParameter<std::string>("outputPlot", "")),
       m_outputRootFile(iConfig.getUntrackedParameter<std::string>("outputRootFile", "")),
       m_outputCSV(iConfig.getUntrackedParameter<std::string>("outputCSV", "")),
       fout(nullptr),
       geomDetToken_(esConsumes()) {
   m_connectionPool.setParameters(iConfig.getParameter<edm::ParameterSet>("DBParameters"));
   m_connectionPool.configure();
   if (!m_outputRootFile.empty())
     fout = new TFile(m_outputRootFile.data(), "RECREATE");
 }

◆ ~SiStripDetVOffTrendPlotter()

SiStripDetVOffTrendPlotter::~SiStripDetVOffTrendPlotter ( )

override

Definition at line 83 of file SiStripDetVOffTrendPlotter.cc.

References fout.

                                                         {
   if (fout)
     fout->Close();
 }

Member Function Documentation

◆ analyze()

void SiStripDetVOffTrendPlotter::analyze	(	const edm::Event &	evt,
		const edm::EventSetup &	evtSetup
	)

overridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 88 of file SiStripDetVOffTrendPlotter.cc.

References DummyCfis::c, cond::persistency::Transaction::commit(), gather_cfg::cout, cond::persistency::ConnectionPool::createSession(), data, dumpCSV(), Exception, cond::persistency::Session::fetchPayload(), cond::persistency::IOVArray::find(), dqmdumpme::first, formatIOV(), fout, cond::time::from_boost(), geomDetToken_, edm::EventSetup::getData(), TrackerGeometryUtils::getSiStripDetIds(), getRunInfo::hours, iovMap, m_condDb, m_connectionPool, m_endTime, m_interval, m_outputCSV, m_outputPlot, m_plotTags, m_startTime, jetsAK4_Puppi_cff::payload, PLOT_COLORS, prepGraph(), cond::persistency::Session::readIov(), cond::persistency::IOVProxy::selectAll(), findQualityFiles::size, cond::persistency::Transaction::start(), AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, cond::time::to_boost(), cond::persistency::Session::transaction(), and cond::time::unpack().

                                                                                            {
   // get total number of modules
   const auto num_modules = TrackerGeometryUtils::getSiStripDetIds(evtSetup.getData(geomDetToken_)).size();
 
   // get start and end time for DB query
   boost::posix_time::ptime p_start, p_end;
   if (m_interval > 0) {
     // from m_interval hours ago to now
     p_end = boost::posix_time::second_clock::universal_time();
     p_start = p_end - boost::posix_time::hours(m_interval);
   } else {
     // use start and end time from config file
     p_start = boost::posix_time::time_from_string(m_startTime);
     p_end = boost::posix_time::time_from_string(m_endTime);
   }
   cond::Time_t startIov = cond::time::from_boost(p_start);
   cond::Time_t endIov = cond::time::from_boost(p_end);
   if (startIov > endIov)
     throw cms::Exception("endTime must be greater than startTime!");
   edm::LogInfo("SiStripDetVOffTrendPlotter")
       << "[SiStripDetVOffTrendPlotter::" << __func__ << "] "
       << "Set start IOV " << startIov << " (" << boost::posix_time::to_simple_string(p_start) << ")"
       << "\n ... Set end IOV " << endIov << " (" << boost::posix_time::to_simple_string(p_end) << ")";
 
   // open db session
   edm::LogInfo("SiStripDetVOffTrendPlotter") << "[SiStripDetVOffTrendPlotter::" << __func__ << "] "
                                              << "Query the condition database " << m_condDb;
   cond::persistency::Session condDbSession = m_connectionPool.createSession(m_condDb);
   condDbSession.transaction().start(true);
 
   // loop over all tags to plot
   std::vector<TGraph *> hvgraphs, lvgraphs;
   TLegend *leg_hv = new TLegend(0.6, 0.87, 0.99, 0.99);
   TLegend *leg_lv = new TLegend(0.6, 0.87, 0.99, 0.99);
   for (unsigned itag = 0; itag < m_plotTags.size(); ++itag) {
     auto tag = m_plotTags.at(itag);
     auto color = itag < PLOT_COLORS.size() ? PLOT_COLORS.at(itag) : kGreen + itag;
 
     std::vector<double> vTime;
     std::vector<double> vHVOffPercent, vLVOffPercent;
 
     // query the database
     edm::LogInfo("SiStripDetVOffTrendPlotter") << "[SiStripDetVOffTrendPlotter::" << __func__ << "] "
                                                << "Reading IOVs from tag " << tag;
     cond::persistency::IOVProxy iovProxy = condDbSession.readIov(tag);
     auto iovs = iovProxy.selectAll();
     auto iiov = iovs.find(startIov);
     auto eiov = iovs.find(endIov);
     int niov = 0;
     while (iiov != iovs.end() && (*iiov).since <= (*eiov).since) {
       // convert cond::Time_t to seconds since epoch
       if ((*iiov).since < startIov)
         vTime.push_back(cond::time::unpack(startIov).first);
       else
         vTime.push_back(cond::time::unpack((*iiov).since).first);
       auto payload = condDbSession.fetchPayload<SiStripDetVOff>((*iiov).payloadId);
       vHVOffPercent.push_back(1.0 * payload->getHVoffCounts() / num_modules);
       vLVOffPercent.push_back(1.0 * payload->getLVoffCounts() / num_modules);
       iovMap[(*iiov).since][tag] = {payload->getHVoffCounts(), payload->getLVoffCounts()};
       // debug
       std::cout << boost::posix_time::to_simple_string(cond::time::to_boost((*iiov).since)) << " (" << (*iiov).since
                 << ")"
                 << ", # HV Off=" << std::setw(6) << payload->getHVoffCounts() << ", # LV Off=" << std::setw(6)
                 << payload->getLVoffCounts() << std::endl;
       ++iiov;
       ++niov;
     }
     edm::LogInfo("SiStripDetVOffTrendPlotter")
         << "[SiStripDetVOffTrendPlotter::" << __func__ << "] "
         << "Read " << niov << " IOVs from tag " << tag << " in the specified interval.";
 
     TGraph *hv = new TGraph(vTime.size(), vTime.data(), vHVOffPercent.data());
     prepGraph(hv, TString("HVOff_") + tag, ";UTC;Fraction of HV off", color);
     leg_hv->AddEntry(hv, tag.data(), "LP");
 
     TGraph *lv = new TGraph(vTime.size(), vTime.data(), vLVOffPercent.data());
     prepGraph(lv, TString("LVOff_") + tag, ";UTC;Fraction of LV off", color);
     leg_lv->AddEntry(lv, tag.data(), "LP");
 
     hvgraphs.push_back(hv);
     lvgraphs.push_back(lv);
   }
 
   condDbSession.transaction().commit();
 
   // Make plots
   TCanvas c("c", "c", 1800, 1200);
   c.SetTopMargin(0.12);
   c.SetBottomMargin(0.08);
   c.SetGridx();
   c.SetGridy();
   for (const auto hv : hvgraphs) {
     if (hv == hvgraphs.front())
       hv->Draw("ALP");
     else
       hv->Draw("LPsame");
     if (fout) {
       fout->cd();
       hv->Write();
     }
   }
   leg_hv->Draw();
   std::string plot_postfix =
       !m_outputPlot.empty() ? m_outputPlot : "from_" + formatIOV(startIov) + "_to_" + formatIOV(endIov) + ".png";
   c.Print(("HVOff_" + plot_postfix).data());
 
   c.Clear();
   for (const auto lv : lvgraphs) {
     if (lv == lvgraphs.front())
       lv->Draw("ALP");
     else
       lv->Draw("LPsame");
     if (fout) {
       fout->cd();
       lv->Write();
     }
   }
   leg_lv->Draw();
   c.Print(("LVOff_" + plot_postfix).data());
 
   if (!m_outputCSV.empty()) {
     dumpCSV(true, num_modules);
     dumpCSV(false, num_modules);
   }
 }

◆ dumpCSV()

void SiStripDetVOffTrendPlotter::dumpCSV	(	bool	isHV,
		std::size_t	nModules
	)

private

Definition at line 240 of file SiStripDetVOffTrendPlotter.cc.

References submitPVResolutionJobs::count, alignBH_cfg::fixed, iovMap, m_outputCSV, m_plotTags, AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, cond::timestamp, cond::time::to_boost(), and findQualityFiles::v.

Referenced by analyze().

                                                                       {
   std::string outCSV = isHV ? "HVOff_table_" + m_outputCSV : "LVOff_table_" + m_outputCSV;
   std::ofstream csv;
   csv.open(outCSV);
   csv << "IOV,Timestamp(UTC),";
   for (const auto &tag : m_plotTags)
     csv << tag << ",";
   csv << std::endl;
 
   csv << std::fixed << std::setprecision(1);
 
   for (const auto &v : iovMap) {
     auto iov = v.first;
     auto timestamp = boost::posix_time::to_simple_string(cond::time::to_boost(iov));
     csv << iov << "," << timestamp << ",";
     for (const auto &tag : m_plotTags) {
       if (v.second.find(tag) != v.second.end()) {
         int count = isHV ? v.second.at(tag).first : v.second.at(tag).second;
         csv << count << " (" << 100. * count / nModules << "%)";
       }
       csv << ",";
     }
     csv << std::endl;
   }
 
   csv.close();
 }

◆ formatIOV()

std::string SiStripDetVOffTrendPlotter::formatIOV	(	cond::Time_t	iov,
		std::string	format = `"%Y-%m-%d__%H_%M_%S"`
	)

private

Definition at line 214 of file SiStripDetVOffTrendPlotter.cc.

References cms::cuda::stream, and cond::time::to_boost().

Referenced by analyze().

                                                                                 {
   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();
 }

◆ prepGraph()

void SiStripDetVOffTrendPlotter::prepGraph	(	TGraph *	gr,
		TString	name,
		TString	title,
		Color_t	color
	)

private

Definition at line 222 of file SiStripDetVOffTrendPlotter.cc.

References Skims_PA_cff::name, and runGCPTkAlMap::title.

Referenced by analyze().

                                                                                                  {
   gr->SetName(name);
   gr->SetTitle(title);
   gr->SetLineColor(color);
   gr->SetLineWidth(2);
   gr->SetMarkerStyle(20);
   gr->SetMarkerSize(1.5);
   gr->SetMarkerColor(color);
   gr->GetXaxis()->SetTimeDisplay(1);
   gr->GetXaxis()->SetLabelOffset(0.02);
   gr->GetXaxis()->SetTimeFormat("#splitline{%b %d}{%H:%M}");
   gr->GetXaxis()->SetTimeOffset(0, "gmt");
   gr->GetXaxis()->SetLabelSize(0.025);
   gr->GetXaxis()->SetTitleSize(0.025);
   gr->GetXaxis()->SetTitleOffset(1.6);
   gr->GetYaxis()->SetRangeUser(0, 1.05);
 }