#include <EEDcsInfoTask.h>

Inheritance diagram for EEDcsInfoTask:

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

virtual	~EEDcsInfoTask ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

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

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 *	meEEDcsActive_ [18]

MonitorElement *	meEEDcsActiveMap_

MonitorElement *	meEEDcsFraction_

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

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 17 of file EEDcsInfoTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 29 of file EEDcsInfoTask.cc.

References dqmStore_, enableCleanup_, edm::ParameterSet::getUntrackedParameter(), i, meEEDcsActive_, meEEDcsActiveMap_, meEEDcsFraction_, mergeRuns_, cppFunctionSkipper::operator, prefixME_, and AlCaHLTBitMon_QueryRunRegistry::string.

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

EEDcsInfoTask::~EEDcsInfoTask ( )

virtual

Destructor.

Definition at line 47 of file EEDcsInfoTask.cc.

                               {
 
 }

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 291 of file EEDcsInfoTask.cc.

                                                                     {
 
 }

void EEDcsInfoTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 51 of file EEDcsInfoTask.cc.

References DQMStore::book2D(), DQMStore::bookFloat(), dqmStore_, MonitorElement::Fill(), i, meEEDcsActive_, meEEDcsActiveMap_, meEEDcsFraction_, mergeVDriftHistosByStation::name, prefixME_, Numbers::sEE(), MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                 {
 
   if ( dqmStore_ ) {
 
     std::string name;
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
 
     meEEDcsFraction_ = dqmStore_->bookFloat( "DCSSummary" );
     meEEDcsFraction_->Fill(0.0);
 
     name = "DCSSummaryMap";
     meEEDcsActiveMap_ = dqmStore_->book2D(name, name, 40, 0., 200., 20, 0., 100.);
     meEEDcsActiveMap_->setAxisTitle("ix / ix+100", 1);
     meEEDcsActiveMap_->setAxisTitle("iy", 2);
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/DCSContents");
 
     for (int i = 0; i < 18; i++) {
       meEEDcsActive_[i] = dqmStore_->bookFloat( "EcalEndcap_" + Numbers::sEE(i+1) );
       meEEDcsActive_[i]->Fill(-1.0);
     }
 
   }
 
 }

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

protectedvirtual

BeginLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 84 of file EEDcsInfoTask.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 lumiQTWidget::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<EcalDCSTowerStatusRcd>() ) ) {
     edm::LogWarning("EEDcsInfoTask") << "EcalDCSTowerStatus record not found";
     return;
   }
 
   edm::ESHandle<EcalDCSTowerStatus> pDCSStatus;
   iSetup.get<EcalDCSTowerStatusRcd>().get(pDCSStatus);
   if ( !pDCSStatus.isValid() ) {
     edm::LogWarning("EEDcsInfoTask") << "EcalDCSTowerStatus record not valid";
     return;
   }
   const EcalDCSTowerStatus* dcsStatus = pDCSStatus.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
     EcalDCSTowerStatus::const_iterator dcsStatusIt = dcsStatus->find( scs[u].rawId() );
     if ( dcsStatusIt != dcsStatus->end() ) dbStatus = dcsStatusIt->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;
     }
       }
     }
   }
 
 //   for(int iz=-1; iz<=1; iz+=2) {
 //     for(int itx=0 ; itx<20; itx++) {
 //       for(int ity=0 ; ity<20; ity++) {
 //         if (EcalScDetId::validDetId(itx+1,ity+1,iz )){
 
 //           EcalScDetId eeid(itx+1,ity+1,iz);
 
 //           uint16_t dbStatus = 0; // 0 = good
 //           EcalDCSTowerStatus::const_iterator dcsStatusIt = dcsStatus->find( eeid.rawId() );
 //           if ( dcsStatusIt != dcsStatus->end() ) dbStatus = dcsStatusIt->getStatusCode();
 
 //           if ( dbStatus > 0 ) {
 //             int offsetSC = (iz > 0) ? 0 : 20;
 //             readyRun[offsetSC+itx][ity] = 0;
 //             readyLumi[offsetSC+itx][ity] = 0;
 //           }
 
 //         }
 //       }
 //     }
 //   }
 
 }

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

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 175 of file EEDcsInfoTask.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 EEDcsInfoTask::cleanup ( void )

protected

Cleanup.

Definition at line 210 of file EEDcsInfoTask.cc.

References dqmStore_, MonitorElement::getName(), i, meEEDcsActive_, meEEDcsActiveMap_, meEEDcsFraction_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().

Referenced by endJob().

                                {
 
   if ( dqmStore_ ) {
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
 
     if ( meEEDcsFraction_ ) dqmStore_->removeElement( meEEDcsFraction_->getName() );
 
     if ( meEEDcsActiveMap_ ) dqmStore_->removeElement( meEEDcsActiveMap_->getName() );
 
     dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/DCSContents");
 
     for (int i = 0; i < 18; i++) {
       if ( meEEDcsActive_[i] ) dqmStore_->removeElement( meEEDcsActive_[i]->getName() );
     }
 
   }
 
 }

void EEDcsInfoTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 78 of file EEDcsInfoTask.cc.

References cleanup(), and enableCleanup_.

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

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

protectedvirtual

EndLuminosityBlock.

Reimplemented from edm::EDAnalyzer.

Definition at line 165 of file EEDcsInfoTask.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 EEDcsInfoTask::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 187 of file EEDcsInfoTask.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 EEDcsInfoTask::fillMonitorElements	(	int	ready[40][20],
		const EcalElectronicsMapping *	map
	)

private

Definition at line 230 of file EEDcsInfoTask.cc.

References DccId_, EcalElectronicsMapping::dccTowerConstituents(), MonitorElement::Fill(), EcalElectronicsMapping::getEcalScDetId(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), ecalpyutils::ism(), meEEDcsActive_, meEEDcsActiveMap_, meEEDcsFraction_, nTowerMax_, MonitorElement::setBinContent(), and lumiQTWidget::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(meEEDcsActiveMap_){
     for(int ix=1 ; ix<=meEEDcsActiveMap_->getNbinsX() ; ix++){
       for(int iy=1 ; iy<=meEEDcsActiveMap_->getNbinsY() ; iy++){
     meEEDcsActiveMap_->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(meEEDcsActiveMap_) meEEDcsActiveMap_->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( meEEDcsActive_[iDcc] ) meEEDcsActive_[iDcc]->Fill( readySum[iDcc]/float(nValidChannels[iDcc]) );
   }
 
   if( meEEDcsFraction_ ) meEEDcsFraction_->Fill( readySumTot/float(nValidChannelsTot) );
 
 }

void EEDcsInfoTask::reset ( void )

protected

Reset.

Definition at line 197 of file EEDcsInfoTask.cc.

References i, meEEDcsActive_, meEEDcsActiveMap_, meEEDcsFraction_, and MonitorElement::Reset().

Referenced by beginRun().

                               {
 
   if ( meEEDcsFraction_ ) meEEDcsFraction_->Reset();
 
   for (int i = 0; i < 18; i++) {
     if ( meEEDcsActive_[i] ) meEEDcsActive_[i]->Reset();
   }
 
   if ( meEEDcsActiveMap_ ) meEEDcsActiveMap_->Reset();
 
 }