#include <EBDataCertificationTask.h>

Inheritance diagram for EBDataCertificationTask:

Public Member Functions
	EBDataCertificationTask (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~EBDataCertificationTask ()
	Destructor. More...

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)
	Analyze. More...

void	beginJob (void)
	BeginJob. More...

void	beginLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
	BeginLuminosityBlock. More...

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

void	cleanup (void)
	Cleanup. More...

void	endJob (void)
	EndJob. More...

void	endLuminosityBlock (const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
	EndLuminosityBlock. More...

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

void	reset (void)
	Reset. More...

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
bool	cloneME_

DQMStore *	dqmStore_

bool	enableCleanup_

TH2F *	hDAQ_

TH2F *	hDCS_

TH2F *	hDQM_

TH1F *	hFrontendByLumi_

TH1F *	hIntegrityByLumi_

TH1F *	hSynchronizationByLumi_

MonitorElement *	meEBDataCertification_ [36]

MonitorElement *	meEBDataCertificationSummary_

MonitorElement *	meEBDataCertificationSummaryMap_

MonitorElement *	meEBReportSummary_

MonitorElement *	meEBReportSummaryContents_ [36]

bool	mergeRuns_

std::string	prefixME_

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)

Detailed Description

Definition at line 19 of file EBDataCertificationTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 28 of file EBDataCertificationTask.cc.

References cloneME_, dqmStore_, enableCleanup_, edm::ParameterSet::getUntrackedParameter(), hDAQ_, hDCS_, hDQM_, hFrontendByLumi_, hIntegrityByLumi_, hSynchronizationByLumi_, i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, mergeRuns_, cmsCodeRules.cppFunctionSkipper::operator, and prefixME_.

                                                                           {
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   // cloneME switch
   cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
 
   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
 
   meEBDataCertificationSummary_ = 0;
   meEBDataCertificationSummaryMap_ = 0;
   for (int i = 0; i < 36; i++) {
     meEBDataCertification_[i] = 0;
   }
 
   hDQM_ = 0;
   hDAQ_ = 0;
   hDCS_ = 0;
   hIntegrityByLumi_ = 0;
   hFrontendByLumi_ = 0;
   hSynchronizationByLumi_ = 0;
 
 }

EBDataCertificationTask::~EBDataCertificationTask ( )

virtual

Destructor.

Definition at line 56 of file EBDataCertificationTask.cc.

                                                   {
 
 }

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 383 of file EBDataCertificationTask.cc.

                                                                               {
 
 }

void EBDataCertificationTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 60 of file EBDataCertificationTask.cc.

References DQMStore::book2D(), DQMStore::bookFloat(), dqmStore_, MonitorElement::Fill(), trackerHits::histo, i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, prefixME_, Numbers::sEB(), MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

                                           {
 
   char histo[200];
 
   if ( dqmStore_ ) {
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
 
     sprintf(histo, "CertificationSummary");
     meEBDataCertificationSummary_ = dqmStore_->bookFloat(histo);
     meEBDataCertificationSummary_->Fill(-1.0);
 
     sprintf(histo, "CertificationSummaryMap");
     meEBDataCertificationSummaryMap_ = dqmStore_->book2D(histo,histo, 72, 0., 72., 34, 0., 34.);
     meEBDataCertificationSummaryMap_->setAxisTitle("jphi", 1);
     meEBDataCertificationSummaryMap_->setAxisTitle("jeta", 2);
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
 
     for (int i = 0; i < 36; i++) {
       sprintf(histo, "EcalBarrel_%s", Numbers::sEB(i+1).c_str());
       meEBDataCertification_[i] = dqmStore_->bookFloat(histo);
       meEBDataCertification_[i]->Fill(-1.0);
     }
 
   }
 
 }

void EBDataCertificationTask::beginLuminosityBlock	(	const edm::LuminosityBlock &	lumiBlock,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

BeginLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 95 of file EBDataCertificationTask.cc.

                                                                                                                    {
 
 }

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

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 252 of file EBDataCertificationTask.cc.

References mergeRuns_, and reset().

                                                                               {
 
   if ( ! mergeRuns_ ) this->reset();
 
 }

void EBDataCertificationTask::cleanup ( void )

protected

Cleanup.

Definition at line 353 of file EBDataCertificationTask.cc.

References cloneME_, dqmStore_, edm::getName(), MonitorElement::getName(), hDAQ_, hDCS_, hDQM_, hFrontendByLumi_, hIntegrityByLumi_, hSynchronizationByLumi_, i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().

Referenced by endJob().

                                          {
 
   if ( cloneME_ ) {
     if( hDQM_ ) delete hDQM_;
     if( hDAQ_ ) delete hDAQ_;
     if( hDCS_ ) delete hDCS_;
     if( hIntegrityByLumi_ ) delete hIntegrityByLumi_;
     if( hFrontendByLumi_ ) delete hFrontendByLumi_;
     if( hSynchronizationByLumi_ ) delete hSynchronizationByLumi_;
   }
   hDQM_ = 0;
   hDAQ_ = 0;
   hDCS_ = 0;
   hIntegrityByLumi_ = 0;
   hFrontendByLumi_ = 0;
   hSynchronizationByLumi_ = 0;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
     if ( meEBDataCertificationSummary_ ) dqmStore_->removeElement( meEBDataCertificationSummary_->getName() );
     if ( meEBDataCertificationSummaryMap_ ) dqmStore_->removeElement( meEBDataCertificationSummaryMap_->getName() );
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
     for (int i = 0; i < 36; i++) {
       if ( meEBDataCertification_[i] ) dqmStore_->removeElement( meEBDataCertification_[i]->getName() );
     }
   }
 
 }

void EBDataCertificationTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file EBDataCertificationTask.cc.

References cleanup(), and enableCleanup_.

                                          {
 
   if ( enableCleanup_ ) this->cleanup();
 
 }

void EBDataCertificationTask::endLuminosityBlock	(	const edm::LuminosityBlock &	lumiBlock,
		const edm::EventSetup &	iSetup
	)

protectedvirtual

EndLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 99 of file EBDataCertificationTask.cc.

References abs, cloneME_, dqmStore_, ExpressReco_HICollisions_FallBack::errors, MonitorElement::Fill(), DQMStore::get(), hDAQ_, hDCS_, hFrontendByLumi_, hIntegrityByLumi_, trackerHits::histo, hSynchronizationByLumi_, i, ecalpyutils::ism(), meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, min, prefixME_, reset(), Numbers::sEB(), and MonitorElement::setBinContent().

                                                                                                                    {
 
   this->reset();
 
   char histo[200];
 
   MonitorElement* me;
 
   // evaluate the DQM quality of observables checked by lumi
   float DQMVal[36];
   for (int i = 0; i < 36; i++) {
     DQMVal[i] = -1.;
   }
 
   sprintf(histo, (prefixME_ + "/EBIntegrityTask/EBIT weighted integrity errors by lumi").c_str());
   me = dqmStore_->get(histo);
   hIntegrityByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hIntegrityByLumi_ );
 
   sprintf(histo, (prefixME_ + "/EBStatusFlagsTask/FEStatus/EBSFT weighted frontend errors by lumi").c_str());
   me = dqmStore_->get(histo);
   hFrontendByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hFrontendByLumi_ );
 
   sprintf(histo, (prefixME_ + "/EBRawDataTask/EBRDT FE synchronization errors by lumi").c_str());
   me = dqmStore_->get(histo);
   hSynchronizationByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hSynchronizationByLumi_ );
 
   if( hIntegrityByLumi_ && hFrontendByLumi_ && hSynchronizationByLumi_ ) {
 
     float integrityErrSum = 0.;
     float integrityQual = 1.0;
     float frontendErrSum = 0.;
     float frontendQual = 1.0;
     float synchronizationErrSum = 0.;
     float synchronizationQual = 1.0;
 
     for ( int i=0; i<36; i++) {
       float ismIntegrityQual = 1.0;
       if( hIntegrityByLumi_->GetBinContent(0) > 0 ) {
         float errors = hIntegrityByLumi_->GetBinContent(i+1);
         ismIntegrityQual = 1.0 - errors/hIntegrityByLumi_->GetBinContent(0);
         integrityErrSum += errors;
       }
       float ismFrontendQual = 1.0;
       if( hFrontendByLumi_->GetBinContent(0) > 0 ) {
         float errors = hFrontendByLumi_->GetBinContent(i+1);
         ismFrontendQual = 1.0 - errors/hFrontendByLumi_->GetBinContent(0);
         frontendErrSum += errors;
       }
       float ismSynchronizationQual = 1.0;
       if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) {
         float errors = hSynchronizationByLumi_->GetBinContent(i+1);
         ismSynchronizationQual = 1.0 - errors/hSynchronizationByLumi_->GetBinContent(0);
         synchronizationErrSum += errors;
       }
       float minVal= std::min(ismIntegrityQual,ismFrontendQual);
       DQMVal[i] = std::min(minVal,ismSynchronizationQual);
     }
 
     if( hIntegrityByLumi_->GetBinContent(0) > 0 ) integrityQual = 1.0 - integrityErrSum/hIntegrityByLumi_->GetBinContent(0)/36.;
     if( hFrontendByLumi_->GetBinContent(0) > 0 ) frontendQual = 1.0 - frontendErrSum/hFrontendByLumi_->GetBinContent(0)/36.;
     if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) synchronizationQual = 1.0 - synchronizationErrSum/hSynchronizationByLumi_->GetBinContent(0)/36.;
     float minVal = std::min(integrityQual,frontendQual);
     float totDQMVal = std::min(minVal,synchronizationQual);
 
     sprintf(histo, (prefixME_ + "/EventInfo/reportSummary").c_str());
     me = dqmStore_->get(histo);
     if( me ) me->Fill(totDQMVal);
 
     for ( int i=0; i<36; i++) {
       sprintf(histo, "EcalBarrel_%s", Numbers::sEB(i+1).c_str());
       me = dqmStore_->get(prefixME_ + "/EventInfo/reportSummaryContents/" + histo);
       if( me ) me->Fill(DQMVal[i]);
 
       sprintf(histo, "reportSummaryMap");
       me = dqmStore_->get(prefixME_ + "/EventInfo/" + histo );
       if( me ) {
         for ( int iett = 0; iett < 34; iett++ ) {
           for ( int iptt = 0; iptt < 72; iptt++ ) {
             int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
             if( i == (ism-1) ) me->setBinContent(iptt+1, iett+1, DQMVal[ism-1]);
           }
         }
       }
     }
 
   }
 
   // now combine reduced DQM with DCS and DAQ
   sprintf(histo, (prefixME_ + "/EventInfo/DAQSummaryMap").c_str());
   me = dqmStore_->get(histo);
   hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
 
   sprintf(histo, (prefixME_ + "/EventInfo/DCSSummaryMap").c_str());
   me = dqmStore_->get(histo);
   hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
 
   float sumCert = 0.;
   float sumCertEB[36];
   int nValidChannels = 0;
   int nValidChannelsEB[36];
 
   for (int i = 0; i < 36; i++) {
     sumCertEB[i] = 0.;
     nValidChannelsEB[i] = 0;
   }
 
   for ( int iett = 0; iett < 34; iett++ ) {
     for ( int iptt = 0; iptt < 72; iptt++ ) {
 
       int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
 
       float xvalDQM = DQMVal[ism-1];
 
       float xvalDAQ, xvalDCS;
       xvalDAQ = xvalDCS = -1.;
       float xcert = -1.;
 
       if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( iptt+1, iett+1 );
       if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( iptt+1, iett+1 );
 
       if ( xvalDQM == -1  || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
         // problems: DQM empty or DAQ and DCS not available
         xcert = 0.0;
       } else {
         // do not consider the white value of DAQ and DCS (problems with DB)
         xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
       }
 
       if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->setBinContent( iptt+1, iett+1, xcert );
 
       sumCertEB[ism-1] += xcert;
       nValidChannelsEB[ism-1]++;
 
       sumCert += xcert;
       nValidChannels++;
 
     }
   }
 
   if( meEBDataCertificationSummary_ ) {
     if( nValidChannels>0 ) meEBDataCertificationSummary_->Fill( sumCert/nValidChannels );
     else meEBDataCertificationSummary_->Fill( 0.0 );
   }
 
   for (int i = 0; i < 36; i++) {
     if( meEBDataCertification_[i] ) {
       if( nValidChannelsEB[i]>0 ) meEBDataCertification_[i]->Fill( sumCertEB[i]/nValidChannelsEB[i] );
       else meEBDataCertification_[i]->Fill( 0.0 );
     }
   }
 
 }

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

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 258 of file EBDataCertificationTask.cc.

References abs, cloneME_, dqmStore_, MonitorElement::Fill(), DQMStore::get(), hDAQ_, hDCS_, hDQM_, trackerHits::histo, i, ecalpyutils::ism(), meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, prefixME_, reset(), and MonitorElement::setBinContent().

                                                                             {
 
   this->reset();
 
   char histo[200];
 
   MonitorElement* me;
 
   sprintf(histo, (prefixME_ + "/EventInfo/reportSummaryMap").c_str());
   me = dqmStore_->get(histo);
   hDQM_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDQM_ );
 
   sprintf(histo, (prefixME_ + "/EventInfo/DAQSummaryMap").c_str());
   me = dqmStore_->get(histo);
   hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
 
   sprintf(histo, (prefixME_ + "/EventInfo/DCSSummaryMap").c_str());
   me = dqmStore_->get(histo);
   hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
 
   float sumCert = 0.;
   float sumCertEB[36];
   int nValidChannels = 0;
   int nValidChannelsEB[36];
 
   for (int i = 0; i < 36; i++) {
     sumCertEB[i] = 0.;
     nValidChannelsEB[i] = 0;
   }
 
   for ( int iett = 0; iett < 34; iett++ ) {
     for ( int iptt = 0; iptt < 72; iptt++ ) {
 
       float xvalDQM, xvalDAQ, xvalDCS;
       xvalDQM = xvalDAQ = xvalDCS = -1.;
       float xcert = -1;
 
       if ( hDQM_ ) xvalDQM = hDQM_->GetBinContent( iptt+1, iett+1 );
       if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( iptt+1, iett+1 );
       if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( iptt+1, iett+1 );
 
       if ( xvalDQM == -1 || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
         // problems: DQM empty or DAQ and DCS not available
         xcert = 0.0;
       } else {
         // do not consider the white value of DAQ and DCS (problems with DB)
         xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
       }
 
       if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->setBinContent( iptt+1, iett+1, xcert );
 
       int ism = ( iett<17 ) ? iptt/4+1 : 18+iptt/4+1;
 
       sumCertEB[ism-1] += xcert;
       nValidChannelsEB[ism-1]++;
 
       sumCert += xcert;
       nValidChannels++;
 
     }
   }
 
   if( meEBDataCertificationSummary_ ) {
     if( nValidChannels>0 ) {
       meEBDataCertificationSummary_->Fill( sumCert/nValidChannels );
     } else {
       meEBDataCertificationSummary_->Fill( 0.0 );
     }
   }
 
   for (int i = 0; i < 36; i++) {
     if( meEBDataCertification_[i] ) {
       if( nValidChannelsEB[i]>0 ) {
         meEBDataCertification_[i]->Fill( sumCertEB[i]/nValidChannelsEB[i] );
       } else {
         meEBDataCertification_[i]->Fill( 0.0 );
       }
     }
   }
 
 }

void EBDataCertificationTask::reset ( void )

protected

Reset.

Definition at line 340 of file EBDataCertificationTask.cc.

References i, meEBDataCertification_, meEBDataCertificationSummary_, meEBDataCertificationSummaryMap_, and MonitorElement::Reset().

Referenced by beginRun(), endLuminosityBlock(), and endRun().

                                         {
 
   if ( meEBDataCertificationSummary_ ) meEBDataCertificationSummary_->Reset();
 
   for (int i = 0; i < 36; i++) {
     if ( meEBDataCertification_[i] ) meEBDataCertification_[i]->Reset();
   }
 
   if ( meEBDataCertificationSummaryMap_ ) meEBDataCertificationSummaryMap_->Reset();
 
 }