#include <EEDaqInfoTask.h>

Inheritance diagram for EEDaqInfoTask:

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

virtual	~EEDaqInfoTask ()
	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

Static Protected Attributes
static const int	DccId_ [18]

static const int	nTowerMax_ = 44

Private Member Functions
void	fillMonitorElements (int ready[40][20], const EcalElectronicsMapping *)

Private Attributes
DQMStore *	dqmStore_

bool	enableCleanup_

MonitorElement *	meEEDaqActive_ [18]

MonitorElement *	meEEDaqActiveMap_

MonitorElement *	meEEDaqFraction_

bool	mergeRuns_

std::string	prefixME_

int	readyLumi [40][20]

int	readyRun [40][20]

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 EEDaqInfoTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 34 of file EEDaqInfoTask.cc.

References dqmStore_, enableCleanup_, edm::ParameterSet::getUntrackedParameter(), i, meEEDaqActive_, meEEDaqActiveMap_, meEEDaqFraction_, mergeRuns_, cmsCodeRules.cppFunctionSkipper::operator, and prefixME_.

                                                       {
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
 
   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
 
   meEEDaqFraction_ = 0;
   meEEDaqActiveMap_ = 0;
   for (int i = 0; i < 18; i++) {
     meEEDaqActive_[i] = 0;
   }
 
 }

EEDaqInfoTask::~EEDaqInfoTask ( )

virtual

Destructor.

Definition at line 52 of file EEDaqInfoTask.cc.

                               {
 
 }

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 274 of file EEDaqInfoTask.cc.

                                                                     {
 
 }

void EEDaqInfoTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 56 of file EEDaqInfoTask.cc.

References DQMStore::book2D(), DQMStore::bookFloat(), dqmStore_, MonitorElement::Fill(), i, meEEDaqActive_, meEEDaqActiveMap_, meEEDaqFraction_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEE(), MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().

                                 {
 
   if ( dqmStore_ ) {
 
     std::string name;
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
 
     meEEDaqFraction_ = dqmStore_->bookFloat( "DAQSummary" );
     meEEDaqFraction_->Fill(0.0);
 
     name = "DAQSummaryMap";
     meEEDaqActiveMap_ = dqmStore_->book2D(name, name, 40, 0., 200., 20, 0., 100.);
     meEEDaqActiveMap_->setAxisTitle("ix / ix+100", 1);
     meEEDaqActiveMap_->setAxisTitle("iy", 2);
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/DAQContents");
 
     for (int i = 0; i < 18; i++) {
       meEEDaqActive_[i] = dqmStore_->bookFloat( "EcalEndcap_" + Numbers::sEE(i+1) );
       meEEDaqActive_[i]->Fill(0.0);
     }
 
   }
 
 }

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

protectedvirtual

BeginLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 89 of file EEDaqInfoTask.cc.

References DccId_, EcalElectronicsMapping::dccTowerConstituents(), EcalCondTowerObjectContainer< T >::end(), EcalCondTowerObjectContainer< T >::find(), edm::EventSetup::find(), edm::EventSetup::get(), EcalElectronicsMapping::getEcalScDetId(), i, edm::ESHandleBase::isValid(), python.multivaluedict::map(), nTowerMax_, edm::ESHandle< class >::product(), readyLumi, readyRun, and matplotRender::t.

                                                                                                          {
 
   // information is by run, so fill the same for the run and for every lumi section
   for ( int itx = 0; itx < 40; itx++ ) {
     for ( int ity = 0; ity < 20; ity++ ) {
       readyLumi[itx][ity] = 1;
     }
   }
 
   if ( !iSetup.find( edm::eventsetup::EventSetupRecordKey::makeKey<EcalDAQTowerStatusRcd>() ) ) {
     edm::LogWarning("EEDaqInfoTask") << "EcalDAQTowerStatus record not found";
     return;
   }
 
   edm::ESHandle<EcalDAQTowerStatus> pDAQStatus;
   iSetup.get<EcalDAQTowerStatusRcd>().get(pDAQStatus);
   if ( !pDAQStatus.isValid() ) {
     edm::LogWarning("EEDaqInfoTask") << "EcalDAQTowerStatus record not valid";
     return;
   }
   const EcalDAQTowerStatus* daqStatus = pDAQStatus.product();
 
   edm::ESHandle< EcalElectronicsMapping > pElecMapping;
   iSetup.get< EcalMappingRcd >().get(pElecMapping);
   if( !pElecMapping.isValid() ) {
     edm::LogWarning("EEDaqInfoTask") << "EcalElectronicsMapping not available";
     return;
   }
   const EcalElectronicsMapping *map = pElecMapping.product();
 
   std::vector<DetId> crystals;
   std::vector<EcalScDetId> scs;
 
   for(unsigned i=0 ; i<sizeof(DccId_)/sizeof(int) ; i++){
     for(int t=1 ; t<=nTowerMax_ ; t++){
 
       crystals = map->dccTowerConstituents(DccId_[i], t);
       if(!crystals.size()) continue;
 
       scs = map->getEcalScDetId(DccId_[i], t, false);
 
       for(unsigned u=0 ; u<scs.size() ; u++){
 
     uint16_t dbStatus = 0; // 0 = good
     EcalDAQTowerStatus::const_iterator daqStatusIt = daqStatus->find( scs[u].rawId() );
     if ( daqStatusIt != daqStatus->end() ) dbStatus = daqStatusIt->getStatusCode();
     
     if ( dbStatus > 0 ) {
       int jx = scs[u].ix() - 1 + (scs[u].zside()<0 ? 0 : 20);
       int jy = scs[u].iy() - 1;
       readyRun[jx][jy] = 0;
       readyLumi[jx][jy] = 0;
     }
       }
     }
   }
 
 }

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

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 158 of file EEDaqInfoTask.cc.

References mergeRuns_, readyRun, and reset().

                                                                     {
 
   if ( ! mergeRuns_ ) this->reset();
 
   for ( int itx = 0; itx < 40; itx++ ) {
     for ( int ity = 0; ity < 20; ity++ ) {
       readyRun[itx][ity] = 1;
     }
   }
 
 }

void EEDaqInfoTask::cleanup ( void )

protected

Cleanup.

Definition at line 193 of file EEDaqInfoTask.cc.

References dqmStore_, edm::getName(), MonitorElement::getName(), i, meEEDaqActive_, meEEDaqActiveMap_, meEEDaqFraction_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().

Referenced by endJob().

                                {
 
   if ( dqmStore_ ) {
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
 
     if ( meEEDaqFraction_ ) dqmStore_->removeElement( meEEDaqFraction_->getName() );
 
     if ( meEEDaqActiveMap_ ) dqmStore_->removeElement( meEEDaqActiveMap_->getName() );
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/DAQContents");
 
     for (int i = 0; i < 18; i++) {
       if ( meEEDaqActive_[i] ) dqmStore_->removeElement( meEEDaqActive_[i]->getName() );
     }
 
   }
 
 }

void EEDaqInfoTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 83 of file EEDaqInfoTask.cc.

References cleanup(), and enableCleanup_.

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

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

protectedvirtual

EndLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 148 of file EEDaqInfoTask.cc.

References fillMonitorElements(), edm::EventSetup::get(), patZpeak::handle, python.multivaluedict::map(), edm::ESHandle< class >::product(), and readyLumi.

                                                                                                          {
 
   edm::ESHandle< EcalElectronicsMapping > handle;
   iSetup.get< EcalMappingRcd >().get(handle);
   const EcalElectronicsMapping *map = handle.product();
   if( ! map ) edm::LogWarning("EEDaqInfoTask") << "EcalElectronicsMapping not available";
   else this->fillMonitorElements(readyLumi, map);
 
 }

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

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 170 of file EEDaqInfoTask.cc.

References fillMonitorElements(), edm::EventSetup::get(), patZpeak::handle, python.multivaluedict::map(), edm::ESHandle< class >::product(), and readyRun.

                                                                   {
 
   edm::ESHandle< EcalElectronicsMapping > handle;
   c.get< EcalMappingRcd >().get(handle);
   const EcalElectronicsMapping *map = handle.product();
   if( ! map ) edm::LogWarning("EEDaqInfoTask") << "EcalElectronicsMapping not available";
   else this->fillMonitorElements(readyRun, map);
 
 }

void EEDaqInfoTask::fillMonitorElements	(	int	ready[40][20],
		const EcalElectronicsMapping *	map
	)

private

Definition at line 213 of file EEDaqInfoTask.cc.

References DccId_, EcalElectronicsMapping::dccTowerConstituents(), MonitorElement::Fill(), EcalElectronicsMapping::getEcalScDetId(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), ecalpyutils::ism(), meEEDaqActive_, meEEDaqActiveMap_, meEEDaqFraction_, nTowerMax_, MonitorElement::setBinContent(), and matplotRender::t.

Referenced by endLuminosityBlock(), and endRun().

                                                                                             {
 
   float readySum[18];
   int nValidChannels[18];
   for ( int ism = 0; ism < 18; ism++ ) {
     readySum[ism] = 0;
     nValidChannels[ism] = 0;
   }
   float readySumTot = 0.;
   int nValidChannelsTot = 0;
 
   if(meEEDaqActiveMap_){
     for(int ix=1 ; ix<=meEEDaqActiveMap_->getNbinsX() ; ix++){
       for(int iy=1 ; iy<=meEEDaqActiveMap_->getNbinsY() ; iy++){
     meEEDaqActiveMap_->setBinContent( ix, iy, -1.0 );
       }
     }
   }
 
   std::vector<DetId> crystals;
   std::vector<EcalScDetId> scs;
 
   for ( unsigned iDcc = 0; iDcc < sizeof(DccId_)/sizeof(int); iDcc++) {
     for ( int t = 1; t<=nTowerMax_; t++ ) {
 
       crystals = map->dccTowerConstituents(DccId_[iDcc], t);
       if(!crystals.size()) continue;
 
       scs = map->getEcalScDetId(DccId_[iDcc], t, false);
 
       for(unsigned u=0 ; u<scs.size() ; u++){ // most of the time one DCC tower = one SC
 
     int jx = scs[u].ix() - 1 + (scs[u].zside()<0 ? 0 : 20);
     int jy = scs[u].iy() - 1;
 
     if(meEEDaqActiveMap_) meEEDaqActiveMap_->setBinContent( jx+1, jy+1, ready[jx][jy] );
 
     int ncrystals = 0;
 
     for(std::vector<DetId>::const_iterator it=crystals.begin() ; it!=crystals.end() ; ++it){
       EEDetId id(*it);
       if( id.zside() == scs[u].zside() && (id.ix()-1)/5+1 == scs[u].ix() && (id.iy()-1)/5+1 == scs[u].iy() ) ncrystals++;
     }
 
     if(ready[jx][jy]) {
       readySum[iDcc] += ncrystals;
       readySumTot += ncrystals;
     }
 
     nValidChannels[iDcc] += ncrystals;
     nValidChannelsTot += ncrystals;
 
       }
     }
     if( meEEDaqActive_[iDcc] ) meEEDaqActive_[iDcc]->Fill( readySum[iDcc]/float(nValidChannels[iDcc]) );
   }
 
   if( meEEDaqFraction_ ) meEEDaqFraction_->Fill( readySumTot/float(nValidChannelsTot) );
 
 }

void EEDaqInfoTask::reset ( void )

protected

Reset.

Definition at line 180 of file EEDaqInfoTask.cc.

References i, meEEDaqActive_, meEEDaqActiveMap_, meEEDaqFraction_, and MonitorElement::Reset().

Referenced by beginRun().

                               {
 
   if ( meEEDaqFraction_ ) meEEDaqFraction_->Reset();
 
   for (int i = 0; i < 18; i++) {
     if ( meEEDaqActive_[i] ) meEEDaqActive_[i]->Reset();
   }
 
   if ( meEEDaqActiveMap_ ) meEEDaqActiveMap_->Reset();
 
 }