#include <ESTrendTask.h>

Inheritance diagram for ESTrendTask:

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

virtual	~ESTrendTask ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

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

void	beginJob (void)

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

void	cleanup (void)

void	endJob (void)

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

void	reset (void)

void	setup (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::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
long int	current_time_

edm::InputTag	dccCollections_

DQMStore *	dqmStore_

bool	enableCleanup_

MonitorElement *	hESFiberErrTrend_

MonitorElement *	hESFiberErrTrendHr_

MonitorElement *	hESRecHitTrend_ [2][2]

MonitorElement *	hESRecHitTrendHr_ [2][2]

MonitorElement *	hESSLinkErrTrend_

MonitorElement *	hESSLinkErrTrendHr_

int	ievt_

bool	init_

long int	last_time_

bool	mergeRuns_

std::string	prefixME_

edm::InputTag	rechitlabel_

long int	start_time_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Detailed Description

Definition at line 13 of file ESTrendTask.h.

Constructor & Destructor Documentation

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

Definition at line 25 of file ESTrendTask.cc.

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

                                                {
 
   init_ = false;
 
   dqmStore_ = Service<DQMStore>().operator->();
 
   prefixME_       = ps.getUntrackedParameter<string>("prefixME", "EcalPreshower"); 
   enableCleanup_  = ps.getUntrackedParameter<bool>("enableCleanup", false);
   mergeRuns_      = ps.getUntrackedParameter<bool>("mergeRuns", false);
   rechitlabel_    = ps.getParameter<InputTag>("RecHitLabel");
   dccCollections_ = 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;
 
 }

ESTrendTask::~ESTrendTask ( )

virtual

Definition at line 54 of file ESTrendTask.cc.

54 {

55 }

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 180 of file ESTrendTask.cc.

References edm::Event::getByLabel(), ESDCCHeaderBlock::getDCCErrors(), ESDCCHeaderBlock::getFEChannelStatus(), getTProfile(), i, j, HcalObjRepresent::setup(), ecaldqm::shift2Right(), and ESDetId::zside().

                                                              {
 
   if ( ! init_ ) this->setup();
 
   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.getByLabel(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<EcalRecHitCollection> ESRecHit;
   if ( e.getByLabel(rechitlabel_, 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") << rechitlabel_ << " 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::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 57 of file ESTrendTask.cc.

                                {
 
   ievt_ = 0;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/ESTrendTask");
     dqmStore_->rmdir(prefixME_ + "/ESTrendTask");
   }
 
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 68 of file ESTrendTask.cc.

References edm::RunBase::beginTime(), and reset().

                                                             {
 
   if ( ! mergeRuns_ ) this->reset();
 
   start_time_ = (r.beginTime()).unixTime();
 
   //std::cout << "start time : " << start_time_ << std::endl;
 
 }

void ESTrendTask::cleanup ( void )

protected

Definition at line 143 of file ESTrendTask.cc.

References edm::getName(), i, and j.

                               {
 
   if ( ! init_ ) return;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/ESTrendTask");
 
     for (int i=0 ; i<2; ++i)  
       for (int j=0 ; j<2; ++j) {
     if (hESRecHitTrend_[i][j]) dqmStore_->removeElement(hESRecHitTrend_[i][j]->getName());
     hESRecHitTrend_[i][j] = 0;
     if (hESRecHitTrendHr_[i][j]) dqmStore_->removeElement(hESRecHitTrendHr_[i][j]->getName());
     hESRecHitTrendHr_[i][j] = 0;
       }
 
     if (hESSLinkErrTrend_) dqmStore_->removeElement(hESSLinkErrTrend_->getName());
     hESSLinkErrTrend_ = 0;
     if (hESFiberErrTrend_) dqmStore_->removeElement(hESFiberErrTrend_->getName());
     hESFiberErrTrend_ = 0;
     if (hESSLinkErrTrendHr_) dqmStore_->removeElement(hESSLinkErrTrendHr_->getName());
     hESSLinkErrTrendHr_ = 0;
     if (hESFiberErrTrendHr_) dqmStore_->removeElement(hESFiberErrTrendHr_->getName());
     hESFiberErrTrendHr_ = 0;
   }
 
   init_ = false;
 
 }

void ESTrendTask::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 172 of file ESTrendTask.cc.

References edm::cleanup().

                              {
 
   LogInfo("ESTrendTask") << "analyzed " << ievt_ << " events";
 
   if ( enableCleanup_ ) this->cleanup();
 
 }

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

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file ESTrendTask.cc.

                                                           {
 
 }

void ESTrendTask::reset ( void )

protected

Definition at line 82 of file ESTrendTask.cc.

References i, and j.

                             {
 
   for (int i=0 ; i<2; ++i)  
     for (int j=0 ; j<2; ++j) {
       if (hESRecHitTrend_[i][j]) hESRecHitTrend_[i][j]->Reset();
       if (hESRecHitTrendHr_[i][j]) hESRecHitTrendHr_[i][j]->Reset();
     }
 
   if (hESSLinkErrTrend_) hESSLinkErrTrend_->Reset();
   if (hESFiberErrTrend_) hESFiberErrTrend_->Reset();
   if (hESSLinkErrTrendHr_) hESSLinkErrTrendHr_->Reset();
   if (hESFiberErrTrendHr_) hESFiberErrTrendHr_->Reset();
 
 }

void ESTrendTask::setup ( void )

protected

Definition at line 97 of file ESTrendTask.cc.

References interpolateCardsSimple::histo, i, and j.

                             {
 
   init_ = true;
 
   char histo[200];
 
   if ( dqmStore_ ) {
     dqmStore_->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] = dqmStore_->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] = dqmStore_->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_ = dqmStore_->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_ = dqmStore_->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_ = dqmStore_->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_ = dqmStore_->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::shift2Left	(	TProfile *	p,
		int	bins = `1`
	)

protected

Definition at line 319 of file ESTrendTask.cc.

References combineCards::bins, and 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 292 of file ESTrendTask.cc.

References combineCards::bins, and 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 285 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 84 of file ESTrendTask.h.

edm::InputTag ESTrendTask::dccCollections_

private

Definition at line 71 of file ESTrendTask.h.

DQMStore* ESTrendTask::dqmStore_

private

Definition at line 62 of file ESTrendTask.h.

bool ESTrendTask::enableCleanup_

private

Definition at line 66 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESFiberErrTrend_

private

Definition at line 75 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESFiberErrTrendHr_

private

Definition at line 79 of file ESTrendTask.h.

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

private

Definition at line 73 of file ESTrendTask.h.

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

private

Definition at line 77 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESSLinkErrTrend_

private

Definition at line 74 of file ESTrendTask.h.

MonitorElement* ESTrendTask::hESSLinkErrTrendHr_

private

Definition at line 78 of file ESTrendTask.h.

int ESTrendTask::ievt_

private

Definition at line 60 of file ESTrendTask.h.

bool ESTrendTask::init_

private

Definition at line 81 of file ESTrendTask.h.

long int ESTrendTask::last_time_

private

Definition at line 85 of file ESTrendTask.h.

bool ESTrendTask::mergeRuns_

private

Definition at line 68 of file ESTrendTask.h.

std::string ESTrendTask::prefixME_

private

Definition at line 64 of file ESTrendTask.h.

edm::InputTag ESTrendTask::rechitlabel_

private

Definition at line 70 of file ESTrendTask.h.

long int ESTrendTask::start_time_

private

Definition at line 83 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