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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () const

	~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.

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