Inheritance diagram for SiStripDetVOffTrendPlotter:

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

void	endJob () override

	SiStripDetVOffTrendPlotter (const edm::ParameterSet &iConfig)

	~SiStripDetVOffTrendPlotter () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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::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
edm::Service< SiStripDetInfoFileReader >	detidReader

TFile *	fout

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::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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<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)

Detailed Description

Definition at line 26 of file SiStripDetVOffTrendPlotter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 62 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) {
   m_connectionPool.setParameters(iConfig.getParameter<edm::ParameterSet>("DBParameters"));
   m_connectionPool.configure();
   if (!m_outputRootFile.empty())
     fout = new TFile(m_outputRootFile.data(), "RECREATE");
 }

SiStripDetVOffTrendPlotter::~SiStripDetVOffTrendPlotter ( )

override

Definition at line 79 of file SiStripDetVOffTrendPlotter.cc.

References fout.

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

Member Function Documentation

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

override

Definition at line 84 of file SiStripDetVOffTrendPlotter.cc.

References EnergyCorrector::c, create_public_lumi_plots::color, cond::persistency::Transaction::commit(), gather_cfg::cout, cond::persistency::ConnectionPool::createSession(), data, detidReader, dumpCSV(), cond::persistency::IOVProxy::end(), Exception, cond::persistency::Session::fetchPayload(), cond::persistency::IOVProxy::find(), plotBeamSpotDB::first, formatIOV(), fout, cond::time::from_boost(), SiStripDetInfoFileReader::getAllDetIds(), AnalysisDataFormats_SUSYBSMObjects::hv, iovMap, m_condDb, m_connectionPool, m_endTime, m_interval, m_outputCSV, m_outputPlot, m_plotTags, m_startTime, jets_cff::payload, PLOT_COLORS, prepGraph(), cond::persistency::Session::readIov(), cond::persistency::Transaction::start(), AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, cond::time::to_boost(), cond::persistency::Session::transaction(), and cond::time::unpack().

                                                                                            {
   // get total number of modules
   auto num_modules = detidReader->getAllDetIds().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, true);  // load all?
     auto iiov = iovProxy.find(startIov);
     auto eiov = iovProxy.find(endIov);
     int niov = 0;
     while (iiov != iovProxy.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);
     dumpCSV(false);
   }
 }

void SiStripDetVOffTrendPlotter::dumpCSV ( bool isHV )

private

Definition at line 237 of file SiStripDetVOffTrendPlotter.cc.

References KineDebug3::count(), DEFINE_FWK_MODULE, detidReader, alignBH_cfg::fixed, SiStripDetInfoFileReader::getAllDetIds(), iovMap, m_outputCSV, m_plotTags, AlCaHLTBitMon_QueryRunRegistry::string, GlobalPosition_Frontier_DevDB_cff::tag, cond::timestamp, cond::time::to_boost(), and findQualityFiles::v.

Referenced by analyze().

                                                   {
   // get total number of modules
   auto num_modules = detidReader->getAllDetIds().size();
 
   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 / num_modules << "%)";
       }
       csv << ",";
     }
     csv << std::endl;
   }
 
   csv.close();
 }

void SiStripDetVOffTrendPlotter::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 209 of file SiStripDetVOffTrendPlotter.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

209 {}

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

private

Definition at line 211 of file SiStripDetVOffTrendPlotter.cc.

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

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

private

Definition at line 219 of file SiStripDetVOffTrendPlotter.cc.

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