#include <ESTrendTask.h>

Inheritance diagram for ESTrendTask:

Public Member Functions
	ESTrendTask (const edm::ParameterSet &ps)

	~ESTrendTask () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override

void	shift2Left (TProfile *p, int bins=1)

void	shift2Right (TProfile *p, int bins=1)

void	updateTime (const edm::Event &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
long int	current_time_

edm::EDGetTokenT< ESRawDataCollection >	dccCollections_

MonitorElement *	hESFiberErrTrend_

MonitorElement *	hESFiberErrTrendHr_

MonitorElement *	hESRecHitTrend_ [2][2]

MonitorElement *	hESRecHitTrendHr_ [2][2]

MonitorElement *	hESSLinkErrTrend_

MonitorElement *	hESSLinkErrTrendHr_

int	ievt_

long int	last_time_

std::string	prefixME_

edm::EDGetTokenT< ESRecHitCollection >	rechittoken_

long int	start_time_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 16 of file ESTrendTask.h.

Constructor & Destructor Documentation

◆ ESTrendTask()

ESTrendTask::ESTrendTask ( const edm::ParameterSet & ps )

Definition at line 22 of file ESTrendTask.cc.

                                                {
   prefixME_ = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower");
   rechittoken_ = consumes<ESRecHitCollection>(ps.getParameter<InputTag>("RecHitLabel"));
   dccCollections_ = consumes<ESRawDataCollection>(ps.getParameter<InputTag>("ESDCCCollections"));
  
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       hESRecHitTrend_[i][j] = nullptr;
       hESRecHitTrendHr_[i][j] = nullptr;
     }
  
   hESSLinkErrTrend_ = nullptr;
   hESFiberErrTrend_ = nullptr;
   hESSLinkErrTrendHr_ = nullptr;
   hESFiberErrTrendHr_ = nullptr;
  
   start_time_ = 0;
   current_time_ = 0;
   last_time_ = 0;
  
   ievt_ = 0;
 }

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, and dqmiolumiharvest::j.

◆ ~ESTrendTask()

ESTrendTask::~ESTrendTask ( )

inlineoverride

Definition at line 22 of file ESTrendTask.h.

22 {}

Member Function Documentation

◆ analyze()

void ESTrendTask::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 91 of file ESTrendTask.cc.

                                                              {
   ievt_++;
  
   // Collect time information
   updateTime(e);
  
   long int diff_current_start = current_time_ - start_time_;
   long int diff_last_start = last_time_ - start_time_;
   //LogInfo("ESTrendTask") << "time difference is negative in " << ievt_ << " events\n"
   //<< "\tcurrent - start time = " << diff_current_start
   //<< ", \tlast - start time = " << diff_last_start << endl;
  
   //  std::cout << "current_time : " << current_time_ << ", diff : " << diff_current_start << std::endl;
  
   // Calculate time interval and bin width
   //  int minuteBinWidth = int(nBasicClusterMinutely_->getTProfile()->GetXaxis()->GetBinWidth(1));
   int minuteBinWidth = 5;
   long int minuteBinDiff = diff_current_start / 60 / minuteBinWidth - diff_last_start / 60 / minuteBinWidth;
   long int minuteDiff = (current_time_ - last_time_) / 60;
  
   //  int hourBinWidth = int(nBasicClusterHourly_->getTProfile()->GetXaxis()->GetBinWidth(1));
   int hourBinWidth = 1;
   long int hourBinDiff = diff_current_start / 3600 / hourBinWidth - diff_last_start / 3600 / hourBinWidth;
   long int hourDiff = (current_time_ - last_time_) / 3600;
  
   if (minuteDiff >= minuteBinWidth) {
     while (minuteDiff >= minuteBinWidth)
       minuteDiff -= minuteBinWidth;
   }
   if (hourDiff >= hourBinWidth) {
     while (hourDiff >= hourBinWidth)
       hourDiff -= hourBinWidth;
   }
  
   // ES DCC
   Int_t slinkCRCErr = 0;
   Int_t fiberErr = 0;
   vector<int> fiberStatus;
   Handle<ESRawDataCollection> dccs;
   if (e.getByToken(dccCollections_, dccs)) {
     for (ESRawDataCollection::const_iterator dccItr = dccs->begin(); dccItr != dccs->end(); ++dccItr) {
       ESDCCHeaderBlock dcc = (*dccItr);
  
       if (dcc.getDCCErrors() == 101)
         slinkCRCErr++;
  
       fiberStatus = dcc.getFEChannelStatus();
       for (unsigned int i = 0; i < fiberStatus.size(); ++i) {
         if (fiberStatus[i] == 4 || fiberStatus[i] == 8 || fiberStatus[i] == 10 || fiberStatus[i] == 11 ||
             fiberStatus[i] == 12)
           fiberErr++;
       }
     }
   }
  
   shift2Right(hESSLinkErrTrend_->getTProfile(), minuteBinDiff);
   hESSLinkErrTrend_->Fill(minuteDiff, slinkCRCErr);
  
   shift2Right(hESFiberErrTrend_->getTProfile(), minuteBinDiff);
   hESFiberErrTrend_->Fill(minuteDiff, fiberErr);
  
   shift2Right(hESSLinkErrTrendHr_->getTProfile(), hourBinDiff);
   hESSLinkErrTrendHr_->Fill(hourDiff, slinkCRCErr);
  
   shift2Right(hESFiberErrTrendHr_->getTProfile(), hourBinDiff);
   hESFiberErrTrendHr_->Fill(hourDiff, fiberErr);
  
   // ES RecHits
   int zside, plane;
   int nrh[2][2];
   for (int i = 0; i < 2; i++)
     for (int j = 0; j < 2; j++) {
       nrh[i][j] = 0;
     }
  
   Handle<ESRecHitCollection> ESRecHit;
   if (e.getByToken(rechittoken_, ESRecHit)) {
     for (ESRecHitCollection::const_iterator hitItr = ESRecHit->begin(); hitItr != ESRecHit->end(); ++hitItr) {
       ESDetId id = ESDetId(hitItr->id());
  
       zside = id.zside();
       plane = id.plane();
  
       int i = (zside == 1) ? 0 : 1;
       int j = plane - 1;
  
       nrh[i][j]++;
     }
   } else {
     LogWarning("ESTrendTask") << "RecHitCollection not available";
   }
  
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       shift2Right(hESRecHitTrend_[i][j]->getTProfile(), minuteBinDiff);
       hESRecHitTrend_[i][j]->Fill(minuteDiff, nrh[i][j] / (1072 * 32.));
  
       shift2Right(hESRecHitTrendHr_[i][j]->getTProfile(), hourBinDiff);
       hESRecHitTrendHr_[i][j]->Fill(hourDiff, nrh[i][j] / (1072 * 32.));
     }
 }

References edm::SortedCollection< T, SORT >::begin(), MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), ESDCCHeaderBlock::getDCCErrors(), ESDCCHeaderBlock::getFEChannelStatus(), mps_fire::i, dqmiolumiharvest::j, and ecaldqm::zside().

◆ bookHistograms()

void ESTrendTask::bookHistograms	(	DQMStore::IBooker &	iBooker,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 51 of file ESTrendTask.cc.

                                                                                         {
   char histo[200];
  
   iBooker.setCurrentFolder(prefixME_ + "/ESTrendTask");
  
   for (int i = 0; i < 2; ++i)
     for (int j = 0; j < 2; ++j) {
       int iz = (i == 0) ? 1 : -1;
       sprintf(histo, "ES Trending RH Occ per 5 mins Z %d P %d", iz, j + 1);
       hESRecHitTrend_[i][j] = iBooker.bookProfile(histo, histo, 36, 0.0, 180.0, 100, 0.0, 1.0e6, "s");
       hESRecHitTrend_[i][j]->setAxisTitle("Elapse time (Minutes)", 1);
       hESRecHitTrend_[i][j]->setAxisTitle("ES RecHit Occupancy / 5 minutes", 2);
  
       sprintf(histo, "ES Trending RH Occ per hour Z %d P %d", iz, j + 1);
       hESRecHitTrendHr_[i][j] = iBooker.bookProfile(histo, histo, 24, 0.0, 24.0, 100, 0.0, 1.0e6, "s");
       hESRecHitTrendHr_[i][j]->setAxisTitle("Elapse time (Hours)", 1);
       hESRecHitTrendHr_[i][j]->setAxisTitle("ES RecHit Occupancy / hour", 2);
     }
  
   sprintf(histo, "ES Trending SLink CRC Error per 5 mins");
   hESSLinkErrTrend_ = iBooker.bookProfile(histo, histo, 36, 0.0, 180.0, 100, 0.0, 1.0e6, "s");
   hESSLinkErrTrend_->setAxisTitle("Elapse time (Minutes)", 1);
   hESSLinkErrTrend_->setAxisTitle("ES SLink CRC Err / 5 minutes", 2);
  
   sprintf(histo, "ES Trending Fiber Error per 5 mins");
   hESFiberErrTrend_ = iBooker.bookProfile(histo, histo, 36, 0.0, 180.0, 100, 0.0, 1.0e6, "s");
   hESFiberErrTrend_->setAxisTitle("Elapse time (Minutes)", 1);
   hESFiberErrTrend_->setAxisTitle("ES Fiber Err / 5 minutes", 2);
  
   sprintf(histo, "ES Trending SLink CRC Error per hour");
   hESSLinkErrTrendHr_ = iBooker.bookProfile(histo, histo, 24, 0.0, 24.0, 100, 0.0, 1.0e6, "s");
   hESSLinkErrTrendHr_->setAxisTitle("Elapse time (Hours)", 1);
   hESSLinkErrTrendHr_->setAxisTitle("ES SLink CRC Err / hour", 2);
  
   sprintf(histo, "ES Trending Fiber Error per hour");
   hESFiberErrTrendHr_ = iBooker.bookProfile(histo, histo, 24, 0.0, 24.0, 100, 0.0, 1.0e6, "s");
   hESFiberErrTrendHr_->setAxisTitle("Elapse time (Hours)", 1);
   hESFiberErrTrendHr_->setAxisTitle("ES Fiber Err / hour", 2);
 }

References dqm::implementation::IBooker::bookProfile(), timingPdfMaker::histo, mps_fire::i, dqmiolumiharvest::j, dqm::impl::MonitorElement::setAxisTitle(), and dqm::implementation::NavigatorBase::setCurrentFolder().

◆ dqmBeginRun()

void ESTrendTask::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 45 of file ESTrendTask.cc.

                                                                {
   start_time_ = (r.beginTime()).unixTime();
  
   //std::cout << "start time : " << start_time_ << std::endl;
 }

References alignCSCRings::r.

◆ shift2Left()

void ESTrendTask::shift2Left	(	TProfile *	p,
		int	bins = `1`
	)

protected

Definition at line 226 of file ESTrendTask.cc.

                                                   {
   if (bins <= 0)
     return;
  
   if (!p->GetSumw2())
     p->Sumw2();
   int nBins = p->GetXaxis()->GetNbins();
  
   // by shifting n bin to the left, the number of entries are
   // reduced by the number in n bins including the underflow bin.
   double nentries = p->GetEntries();
   for (int i = 0; i < bins; i++)
     nentries -= p->GetBinEntries(i);
   p->SetEntries(nentries);
  
   // the first bin goes to underflow
   // each bin moves to the right
  
   TArrayD* sumw2 = p->GetSumw2();
  
   for (int i = 0; i <= nBins + 1 - bins; i++) {
     // GetBinContent return binContent/binEntries
     p->SetBinContent(i, p->GetBinContent(i + bins) * p->GetBinEntries(i + bins));
     p->SetBinEntries(i, p->GetBinEntries(i + bins));
     sumw2->SetAt(sumw2->GetAt(i + bins), i);
   }
 }

References trigObjTnPSource_cfi::bins, mps_fire::i, seedmultiplicitymonitor_newtracking_cfi::nBins, and AlCaHLTBitMon_ParallelJobs::p.

◆ shift2Right()

void ESTrendTask::shift2Right	(	TProfile *	p,
		int	bins = `1`
	)

protected

Definition at line 198 of file ESTrendTask.cc.

                                                    {
   if (bins <= 0)
     return;
  
   if (!p->GetSumw2())
     p->Sumw2();
   int nBins = p->GetXaxis()->GetNbins();
  
   // by shifting n bin to the right, the number of entries are
   // reduced by the number in n bins including the overflow bin.
   double nentries = p->GetEntries();
   for (int i = 0; i < bins; i++)
     nentries -= p->GetBinEntries(nBins + 1 - bins);
   p->SetEntries(nentries);
  
   // the last bin goes to overflow
   // each bin moves to the right
  
   TArrayD* sumw2 = p->GetSumw2();
  
   for (int i = nBins + 1; i > bins; i--) {
     // GetBinContent return binContent/binEntries
     p->SetBinContent(i, p->GetBinContent(i - bins) * p->GetBinEntries(i - bins));
     p->SetBinEntries(i, p->GetBinEntries(i - bins));
     sumw2->SetAt(sumw2->GetAt(i - bins), i);
   }
 }

References trigObjTnPSource_cfi::bins, mps_fire::i, seedmultiplicitymonitor_newtracking_cfi::nBins, and AlCaHLTBitMon_ParallelJobs::p.

◆ updateTime()

void ESTrendTask::updateTime ( const edm::Event & e )

protected

Definition at line 193 of file ESTrendTask.cc.

                                               {
   last_time_ = current_time_;
   current_time_ = e.time().unixTime();
 }

References MillePedeFileConverter_cfg::e.

Member Data Documentation

◆ current_time_

long int ESTrendTask::current_time_

private

Definition at line 59 of file ESTrendTask.h.

◆ dccCollections_

edm::EDGetTokenT<ESRawDataCollection> ESTrendTask::dccCollections_

private

Definition at line 48 of file ESTrendTask.h.

◆ hESFiberErrTrend_

MonitorElement* ESTrendTask::hESFiberErrTrend_

private

Definition at line 52 of file ESTrendTask.h.

◆ hESFiberErrTrendHr_

MonitorElement* ESTrendTask::hESFiberErrTrendHr_

private

Definition at line 56 of file ESTrendTask.h.

◆ hESRecHitTrend_

MonitorElement* ESTrendTask::hESRecHitTrend_[2][2]

private

Definition at line 50 of file ESTrendTask.h.

◆ hESRecHitTrendHr_

MonitorElement* ESTrendTask::hESRecHitTrendHr_[2][2]

private

Definition at line 54 of file ESTrendTask.h.

◆ hESSLinkErrTrend_

MonitorElement* ESTrendTask::hESSLinkErrTrend_

private

Definition at line 51 of file ESTrendTask.h.

◆ hESSLinkErrTrendHr_

MonitorElement* ESTrendTask::hESSLinkErrTrendHr_

private

Definition at line 55 of file ESTrendTask.h.

◆ ievt_

int ESTrendTask::ievt_

private

Definition at line 43 of file ESTrendTask.h.

◆ last_time_

long int ESTrendTask::last_time_

private

Definition at line 60 of file ESTrendTask.h.

◆ prefixME_

std::string ESTrendTask::prefixME_

private

Definition at line 45 of file ESTrendTask.h.

◆ rechittoken_

edm::EDGetTokenT<ESRecHitCollection> ESTrendTask::rechittoken_

private

Definition at line 47 of file ESTrendTask.h.

◆ start_time_

long int ESTrendTask::start_time_

private

Definition at line 58 of file ESTrendTask.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ ESTrendTask()

◆ ~ESTrendTask()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dqmBeginRun()

◆ shift2Left()

◆ shift2Right()

◆ updateTime()

Member Data Documentation

◆ current_time_

◆ dccCollections_

◆ hESFiberErrTrend_

◆ hESFiberErrTrendHr_

◆ hESRecHitTrend_

◆ hESRecHitTrendHr_

◆ hESSLinkErrTrend_

◆ hESSLinkErrTrendHr_

◆ ievt_

◆ last_time_

◆ prefixME_

◆ rechittoken_

◆ start_time_