#include <ESTrendTask.h>

Inheritance diagram for ESTrendTask:

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

virtual	~ESTrendTask ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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	endJob (void)

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

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

void	updateTime (const edm::Event &)

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

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

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 17 of file ESTrendTask.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file ESTrendTask.cc.

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

                                                {
 
   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] = 0;
       hESRecHitTrendHr_[i][j] = 0;
     }
 
   hESSLinkErrTrend_ = 0;
   hESFiberErrTrend_ = 0;
   hESSLinkErrTrendHr_ = 0;
   hESFiberErrTrendHr_ = 0;
 
   start_time_ = 0;
   current_time_ = 0;
   last_time_ = 0;
 
   ievt_ = 0;
 }

virtual ESTrendTask::~ESTrendTask ( )

inlinevirtual

Definition at line 25 of file ESTrendTask.h.

25 {}

Member Function Documentation

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

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 102 of file ESTrendTask.cc.

References edm::Event::getByToken(), ESDCCHeaderBlock::getDCCErrors(), ESDCCHeaderBlock::getFEChannelStatus(), getTProfile(), i, j, ESDetId::zside(), and ecaldqm::zside().

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

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 55 of file ESTrendTask.cc.

References DQMStore::IBooker::bookProfile(), timingPdfMaker::histo, i, j, MonitorElement::setAxisTitle(), and DQMStore::IBooker::setCurrentFolder().

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

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 47 of file ESTrendTask.cc.

References edm::RunBase::beginTime().

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

void ESTrendTask::endJob ( void )

protected

Definition at line 96 of file ESTrendTask.cc.

                              {
 
   LogInfo("ESTrendTask") << "analyzed " << ievt_ << " events";
 
 }

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

protected

Definition at line 239 of file ESTrendTask.cc.

References i.

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

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

protected

Definition at line 212 of file ESTrendTask.cc.

References i.

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

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

protected

Definition at line 205 of file ESTrendTask.cc.

References edm::EventBase::time(), and edm::Timestamp::unixTime().

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

Member Data Documentation

long int ESTrendTask::current_time_

private

Definition at line 67 of file ESTrendTask.h.

edm::EDGetTokenT<ESRawDataCollection> ESTrendTask::dccCollections_

private

Definition at line 56 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESFiberErrTrend_

private

Definition at line 60 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESFiberErrTrendHr_

private

Definition at line 64 of file ESTrendTask.h.

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

private

Definition at line 58 of file ESTrendTask.h.

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

private

Definition at line 62 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESSLinkErrTrend_

private

Definition at line 59 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESSLinkErrTrendHr_

private

Definition at line 63 of file ESTrendTask.h.

int ESTrendTask::ievt_

private

Definition at line 51 of file ESTrendTask.h.

long int ESTrendTask::last_time_

private

Definition at line 68 of file ESTrendTask.h.

std::string ESTrendTask::prefixME_

private

Definition at line 53 of file ESTrendTask.h.

edm::EDGetTokenT<ESRecHitCollection> ESTrendTask::rechittoken_

private

Definition at line 55 of file ESTrendTask.h.

long int ESTrendTask::start_time_

private

Definition at line 66 of file ESTrendTask.h.

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation