#include <OccupancyClient.h>

Inheritance diagram for ecaldqm::OccupancyClient:

Public Types
enum	MESets { kHotDigi, kHotRecHitThr, kHotTPDigiThr, kQualitySummary, nMESets }

enum	Sources { sDigi, sRecHitThr, sTPDigiThr, nSources }

Public Types inherited from ecaldqm::DQWorkerClient
enum	Sources { nSources }

Public Types inherited from ecaldqm::DQWorker
enum	MESets { nMESets }

Public Member Functions
void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	bookMEs () override

	OccupancyClient (const edm::ParameterSet &, const edm::ParameterSet &)

void	producePlots () override

	~OccupancyClient ()

Public Member Functions inherited from ecaldqm::DQWorkerClient
	DQWorkerClient (const edm::ParameterSet &, const edm::ParameterSet &, std::string const &)

virtual void	endLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)

virtual void	initialize ()

void	reset ()

virtual	~DQWorkerClient ()

Public Member Functions inherited from ecaldqm::DQWorker
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &)

	DQWorker (const edm::ParameterSet &, const edm::ParameterSet &, std::string const &)

virtual void	endRun (const edm::Run &, const edm::EventSetup &)

const std::vector< MESet * > &	getMEs ()

virtual std::string const &	getName ()

virtual bool	isInitialized ()

virtual void	setInitialized (bool _init)

virtual void	setVerbosity (int _verbosity)

virtual	~DQWorker ()

Static Public Member Functions
static void	setMEData (std::vector< MEData > &)

Static Public Member Functions inherited from ecaldqm::DQWorker
static void	setMEData (std::vector< MEData > &)

Private Attributes
float	deviationThreshold_

const CaloGeometry *	geometry_

int	minHits_

Additional Inherited Members
Static Public Attributes inherited from ecaldqm::DQWorkerClient
static EcalDQMChannelStatus const *	channelStatus

static EcalDQMTowerStatus const *	towerStatus

Static Public Attributes inherited from ecaldqm::DQWorker
static std::map< std::string, std::vector< MEData > >	meData

Protected Member Functions inherited from ecaldqm::DQWorkerClient
void	fillQuality_ (unsigned, DetId const &, uint32_t, float)

void	source_ (unsigned, std::string const &, unsigned, edm::ParameterSet const &)

Protected Member Functions inherited from ecaldqm::DQWorker
MESet *	createMESet_ (std::string const &, MEData const &, bool _readOnly=false) const

void	meSet_ (unsigned, edm::ParameterSet const &)

Protected Attributes inherited from ecaldqm::DQWorkerClient
std::vector< MESet const * >	sources_

Protected Attributes inherited from ecaldqm::DQWorker
bool	initialized_

std::vector< MESet * >	MEs_

std::string	name_

int	verbosity_

Detailed Description

Definition at line 10 of file OccupancyClient.h.

Member Enumeration Documentation

enum ecaldqm::OccupancyClient::MESets

Enumerator
kHotDigi
kHotRecHitThr
kHotTPDigiThr
kQualitySummary
nMESets

Definition at line 21 of file OccupancyClient.h.

                 {
       kHotDigi,
       kHotRecHitThr,
       kHotTPDigiThr,
       kQualitySummary,
       nMESets
     };

enum ecaldqm::OccupancyClient::Sources

Enumerator
sDigi
sRecHitThr
sTPDigiThr
nSources

Definition at line 31 of file OccupancyClient.h.

                  {
       sDigi,
       sRecHitThr,
       sTPDigiThr,
       nSources
     };

Constructor & Destructor Documentation

ecaldqm::OccupancyClient::OccupancyClient	(	const edm::ParameterSet &	_params,
		const edm::ParameterSet &	_paths
	)

Definition at line 15 of file OccupancyClient.cc.

References deviationThreshold_, edm::ParameterSet::getUntrackedParameterSet(), ecaldqm::OccupancyTask::kDigi, ecaldqm::OccupancyTask::kRecHitThr, ecaldqm::OccupancyTask::kTPDigiThr, minHits_, ecaldqm::DQWorker::name_, sDigi, ecaldqm::DQWorkerClient::source_(), sRecHitThr, and sTPDigiThr.

                                                                                                 :
     DQWorkerClient(_params, _paths, "OccupancyClient"),
     geometry_(0),
     minHits_(0),
     deviationThreshold_(0.)
   {
     edm::ParameterSet const& taskParams(_params.getUntrackedParameterSet(name_));
     minHits_ = taskParams.getUntrackedParameter<int>("minHits");
     deviationThreshold_ = taskParams.getUntrackedParameter<double>("deviationThreshold");
 
     edm::ParameterSet const& sources(_params.getUntrackedParameterSet("sources"));
     source_(sDigi, "OccupancyTask", OccupancyTask::kDigi, sources);
     source_(sRecHitThr, "OccupancyTask", OccupancyTask::kRecHitThr, sources);
     source_(sTPDigiThr, "OccupancyTask", OccupancyTask::kTPDigiThr, sources);
   }

ecaldqm::OccupancyClient::~OccupancyClient ( )

inline

Definition at line 13 of file OccupancyClient.h.

13 {}

Member Function Documentation

void ecaldqm::OccupancyClient::beginRun	(	const edm::Run &	,
		const edm::EventSetup &	_es
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 32 of file OccupancyClient.cc.

References edm::hlt::Exception, geometry_, edm::EventSetup::get(), and edm::ESHandle< class >::product().

   {
     edm::ESHandle<CaloGeometry> geomHndl;
     _es.get<CaloGeometryRecord>().get(geomHndl);
     geometry_ = geomHndl.product();
     if(!geometry_)
       throw cms::Exception("EventSetup") << "CaloGeometry invalid";
   }

void ecaldqm::OccupancyClient::bookMEs ( )

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 42 of file OccupancyClient.cc.

References ecaldqm::DQWorker::bookMEs(), kQualitySummary, and ecaldqm::DQWorker::MEs_.

   {
     DQWorker::bookMEs();
 
     MEs_[kQualitySummary]->resetAll(-1.);
   }

void ecaldqm::OccupancyClient::producePlots ( )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 50 of file OccupancyClient.cc.

References EcalTrigTowerDetId::detIdFromDenseIndex(), deviationThreshold_, python.tagInventory::entries, ecaldqm::DQWorkerClient::fillQuality_(), ecaldqm::getElectronicsMap(), ecaldqm::getNSuperCrystals(), ecaldqm::getTrigTowerMap(), EcalTrigTowerDetId::ietaAbs(), kHotDigi, kHotRecHitThr, kHotTPDigiThr, kQualitySummary, EcalTrigTowerDetId::kSizeForDenseIndexing, ecaldqm::DQWorker::MEs_, minHits_, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, sDigi, ecaldqm::DQWorkerClient::sources_, sRecHitThr, and sTPDigiThr.

   {
     using namespace std;
 
     MEs_[kHotDigi]->reset();
     MEs_[kHotRecHitThr]->reset();
     MEs_[kHotTPDigiThr]->reset();
 
     uint32_t mask(1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR);
 
     vector<double> digiPhiRingMean(28, 0.);
     vector<double> rechitPhiRingMean(28, 0.);
     vector<int> numCrystals(28, 0); // this is static, but is easier to count now
 
     for(unsigned dccid(1); dccid <= 54; dccid++){
       for(unsigned tower(1); tower <= getNSuperCrystals(dccid); tower++){
     vector<DetId> ids(getElectronicsMap()->dccTowerConstituents(dccid, tower));
 
     if(ids.size() == 0) continue;
 
     for(vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
       float entries(sources_[sDigi]->getBinContent(*idItr));
       float rhentries(sources_[sRecHitThr]->getBinContent(*idItr));
 
       int ieta(getTrigTowerMap()->towerOf(*idItr).ietaAbs());
       digiPhiRingMean.at(ieta - 1) += entries;
       rechitPhiRingMean.at(ieta - 1) += rhentries;
 
       numCrystals.at(ieta - 1) += 1;
     }
       }
     }
 
     for(int ie(0); ie < 28; ie++){
       digiPhiRingMean[ie] /= numCrystals[ie];
       rechitPhiRingMean[ie] /= numCrystals[ie];
     }
 
     // second round to find hot towers
     for(unsigned dccid(1); dccid <= 54; dccid++){
       for(unsigned tower(1); tower <= getNSuperCrystals(dccid); tower++){
     vector<DetId> ids(getElectronicsMap()->dccTowerConstituents(dccid, tower));
 
     if(ids.size() == 0) continue;
 
     float quality(1.);
     for(vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
       float entries(sources_[sDigi]->getBinContent(*idItr));
       float rhentries(sources_[sRecHitThr]->getBinContent(*idItr));
 
       int ieta(getTrigTowerMap()->towerOf(*idItr).ietaAbs());
 
       if(entries > minHits_ && entries > digiPhiRingMean.at(ieta - 1) * deviationThreshold_){
         MEs_[kHotDigi]->fill(*idItr);
         quality = 0.;
       }
       if(rhentries > minHits_ && rhentries > rechitPhiRingMean.at(ieta - 1) * deviationThreshold_){
         MEs_[kHotRecHitThr]->fill(*idItr);
         quality = 0.;
       }
     }
     if(dccid <= 9 || dccid >= 46){
       vector<EcalScDetId> scs(getElectronicsMap()->getEcalScDetId(dccid, tower));
       for(vector<EcalScDetId>::iterator scItr(scs.begin()); scItr != scs.end(); ++scItr)
         fillQuality_(kQualitySummary, *scItr, mask, quality);
     }
     else
       fillQuality_(kQualitySummary, ids[0], mask, quality);
       }
     }
 
     vector<double> tpdigiPhiRingMean(28, 0.);
 
     for(unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++){
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
       float entries(sources_[sTPDigiThr]->getBinContent(ttid));
 
       tpdigiPhiRingMean.at(ttid.ietaAbs() - 1) += entries;
     }
 
     for(int ie(0); ie < 28; ie++){
       float denom(-1.);
       if(ie < 27) denom = 72.;
       else denom = 36.;
       tpdigiPhiRingMean[ie] /= denom;
     }
 
     for(unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; iTT++){
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
       float entries(sources_[sTPDigiThr]->getBinContent(ttid));
 
       if(entries > minHits_ && entries > tpdigiPhiRingMean.at(ttid.ietaAbs() - 1) * deviationThreshold_){
     MEs_[kHotTPDigiThr]->fill(ttid);
     vector<DetId> ids(getTrigTowerMap()->constituentsOf(ttid));
     for(vector<DetId>::iterator idItr(ids.begin()); idItr != ids.end(); ++idItr){
       if(MEs_[kQualitySummary]->getBinContent(*idItr) > 0.)
         fillQuality_(kQualitySummary, *idItr, mask, 0.);
     }
       }   
     }
 
   }

void ecaldqm::OccupancyClient::setMEData ( std::vector< MEData > & _data )

static

Definition at line 155 of file OccupancyClient.cc.

References MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2F, EcalDQMBinningService::kChannel, EcalDQMBinningService::kCrystal, EcalDQMBinningService::kEcal2P, kHotDigi, kHotRecHitThr, kHotTPDigiThr, kQualitySummary, EcalDQMBinningService::kSuperCrystal, and EcalDQMBinningService::kTriggerTower.

   {
     _data[kHotDigi] = MEData("HotDigi", BinService::kChannel, BinService::kCrystal, MonitorElement::DQM_KIND_TH1F);
     _data[kHotRecHitThr] = MEData("HotRecHitThr", BinService::kChannel, BinService::kCrystal, MonitorElement::DQM_KIND_TH1F);
     _data[kHotTPDigiThr] = MEData("HotTPDigiThr", BinService::kChannel, BinService::kTriggerTower, MonitorElement::DQM_KIND_TH1F);
     _data[kQualitySummary] = MEData("QualitySummary", BinService::kEcal2P, BinService::kSuperCrystal, MonitorElement::DQM_KIND_TH2F);
   }

Member Data Documentation

float ecaldqm::OccupancyClient::deviationThreshold_

private

Definition at line 42 of file OccupancyClient.h.

Referenced by OccupancyClient(), and producePlots().

const CaloGeometry* ecaldqm::OccupancyClient::geometry_

private

Definition at line 39 of file OccupancyClient.h.

Referenced by beginRun().

int ecaldqm::OccupancyClient::minHits_

private

Definition at line 41 of file OccupancyClient.h.

Referenced by OccupancyClient(), and producePlots().

Public Types

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation