#include <OccupancyClient.h>

Inheritance diagram for ecaldqm::OccupancyClient:

Public Member Functions
	OccupancyClient ()

void	producePlots (ProcessType) override

	~OccupancyClient () override

Public Member Functions inherited from ecaldqm::DQWorkerClient
void	bookMEs (DQMStore::IBooker &) override

	DQWorkerClient ()

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

void	releaseMEs () override

void	releaseSource ()

virtual void	resetMEs ()

void	resetPerLumi ()

bool	retrieveSource (DQMStore::IGetter &, ProcessType)

bool	runsOn (ProcessType _type) const

void	setStatusManager (StatusManager const &_manager)

	~DQWorkerClient () override

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

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

	DQWorker ()

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

const EcalDQMSetupObjects	getEcalDQMSetupObjects ()

const EcalElectronicsMapping *	GetElectronicsMap ()

const CaloGeometry *	GetGeometry ()

const std::string &	getName () const

const CaloTopology *	GetTopology ()

const EcalTrigTowerConstituentsMap *	GetTrigTowerMap ()

bool	onlineMode () const

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setSetupObjects (edm::EventSetup const &)

void	setTime (time_t _t)

virtual	~DQWorker () noexcept(false)

Private Member Functions
void	setParams (edm::ParameterSet const &) override

Private Attributes
float	deviationThreshold_

int	minHits_

Additional Inherited Members
Public Types inherited from ecaldqm::DQWorkerClient
enum	ProcessType { kLumi, kJob, nProcessType }

enum	Quality { kBad = 0, kGood = 1, kUnknown = 2, kMBad = 3, kMGood = 4, kMUnknown = 5 }

Static Public Member Functions inherited from ecaldqm::DQWorkerClient
static void	fillDescriptions (edm::ParameterSetDescription &)

Static Public Member Functions inherited from ecaldqm::DQWorker
static void	fillDescriptions (edm::ParameterSetDescription &_desc)

Protected Types inherited from ecaldqm::DQWorker
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Protected Member Functions inherited from ecaldqm::DQWorkerClient
void	setME (edm::ParameterSet const &_ps) final

void	setSource (edm::ParameterSet const &) override

void	towerAverage_ (MESet &, MESet const &, float)

bool	using_ (std::string const &_name, ProcessType _type=kJob) const

Protected Member Functions inherited from ecaldqm::DQWorker
void	initialize (std::string const &_name, edm::ParameterSet const &)

void	print_ (std::string const &, int=0) const

void	setVerbosity (int _verbosity)

Protected Attributes inherited from ecaldqm::DQWorkerClient
bool	hasLumiPlots_

std::set< std::string >	qualitySummaries_

MESetCollection	sources_

StatusManager const *	statusManager_

Protected Attributes inherited from ecaldqm::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 8 of file OccupancyClient.h.

Constructor & Destructor Documentation

◆ OccupancyClient()

ecaldqm::OccupancyClient::OccupancyClient ( )

Definition at line 12 of file OccupancyClient.cc.

                                    : DQWorkerClient(), minHits_(0), deviationThreshold_(0.) {
     qualitySummaries_.insert("QualitySummary");
   }

References ecaldqm::DQWorkerClient::qualitySummaries_.

◆ ~OccupancyClient()

ecaldqm::OccupancyClient::~OccupancyClient ( )

inlineoverride

Definition at line 11 of file OccupancyClient.h.

11 {}

Member Function Documentation

◆ producePlots()

void ecaldqm::OccupancyClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 21 of file OccupancyClient.cc.

                                                 {
     using namespace std;
  
     // number of allowed ieta indices
     // EE-: -28 to -1 with -27, -25 empty
     // EE+: 1 to 28 with 26, 28 empty
     unsigned const nPhiRings(56);
  
     MESet& meQualitySummary(MEs_.at("QualitySummary"));
     //     MESet& meHotDigi(MEs_.at("HotDigi"));
     //     MESet& meHotRecHitThr(MEs_.at("HotRecHitThr"));
     //     MESet& meHotTPDigiThr(MEs_.at("HotTPDigiThr"));
  
     MESet const& sDigi(sources_.at("DigiAll"));
     MESet const& sRecHitThr(sources_.at("RecHitThrAll"));
     MESet const& sTPDigiThr(sources_.at("TPDigiThrAll"));
  
     uint32_t mask(1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR |
                   1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING |
                   1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_ERROR);
  
     double digiPhiRingMean[nPhiRings];
     std::fill_n(digiPhiRingMean, nPhiRings, 0.);
     double rechitPhiRingMean[nPhiRings];
     std::fill_n(rechitPhiRingMean, nPhiRings, 0.);
     int numCrystals[nPhiRings];  // this is static, but is easier to count now
     std::fill_n(numCrystals, nPhiRings, 0);
  
     MESet::const_iterator dEnd(sDigi.end(GetElectronicsMap()));
     MESet::const_iterator rItr(GetElectronicsMap(), sRecHitThr);
     for (MESet::const_iterator dItr(sDigi.beginChannel(GetElectronicsMap())); dItr != dEnd;
          dItr.toNextChannel(GetElectronicsMap())) {
       rItr = dItr;
  
       float entries(dItr->getBinContent());
       float rhentries(rItr->getBinContent());
  
       DetId id(dItr->getId());
       int ieta(0);
       if (id.subdetId() == EcalTriggerTower)  // barrel
         ieta = EcalTrigTowerDetId(id).ieta();
       else {
         std::vector<DetId> ids(scConstituents(EcalScDetId(id)));
         if (ids.empty())
           continue;
         ieta = GetTrigTowerMap()->towerOf(ids[0]).ieta();
       }
  
       unsigned index(ieta < 0 ? ieta + 28 : ieta + 27);
  
       digiPhiRingMean[index] += entries;
       rechitPhiRingMean[index] += rhentries;
       numCrystals[index] += 1;
     }
  
     for (unsigned ie(0); ie < nPhiRings; ie++) {
       digiPhiRingMean[ie] /= numCrystals[ie];
       rechitPhiRingMean[ie] /= numCrystals[ie];
     }
  
     // Store # of entries for Occupancy analysis
     std::vector<float> Nentries(nDCC, 0.);    // digis
     std::vector<float> Nrhentries(nDCC, 0.);  // (filtered) rechits
  
     // second round to find hot towers
     for (MESet::const_iterator dItr(sDigi.beginChannel(GetElectronicsMap())); dItr != dEnd;
          dItr.toNextChannel(GetElectronicsMap())) {
       DetId id(dItr->getId());
  
       bool doMask(meQualitySummary.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
  
       rItr = dItr;
  
       float entries(dItr->getBinContent());
       float rhentries(rItr->getBinContent());
  
       int ieta(0);
       if (id.subdetId() == EcalTriggerTower)  // barrel
         ieta = EcalTrigTowerDetId(id).ieta();
       else {
         std::vector<DetId> ids(scConstituents(EcalScDetId(id)));
         if (ids.empty())
           continue;
         ieta = GetTrigTowerMap()->towerOf(ids[0]).ieta();
       }
  
       unsigned index(ieta < 0 ? ieta + 28 : ieta + 27);
  
       int quality(doMask ? kMGood : kGood);
  
       if (entries > minHits_ && entries > digiPhiRingMean[index] * deviationThreshold_) {
         //        meHotDigi->fill(id);
         quality = doMask ? kMBad : kBad;
       }
       if (rhentries > minHits_ && rhentries > rechitPhiRingMean[index] * deviationThreshold_) {
         //        meHotRecHitThr->fill(id);
         quality = doMask ? kMBad : kBad;
       }
  
       meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, double(quality));
  
       // Keep count of digis & rechits for Occupancy analysis
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (entries > minHits_)
         Nentries[iDCC] += entries;
       if (rhentries > minHits_)
         Nrhentries[iDCC] += rhentries;
     }
  
     double tpdigiPhiRingMean[nPhiRings];
     std::fill_n(tpdigiPhiRingMean, nPhiRings, 0.);
  
     for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; ++iTT) {
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
       float entries(sTPDigiThr.getBinContent(getEcalDQMSetupObjects(), ttid));
  
       unsigned index(ttid.ieta() < 0 ? ttid.ieta() + 28 : ttid.ieta() + 27);
  
       tpdigiPhiRingMean[index] += entries;
     }
  
     for (int ie(0); ie < 28; ie++) {
       float denom(-1.);
       if (ie < 27)
         denom = 72.;
       else
         denom = 36.;
       tpdigiPhiRingMean[ie] /= denom;
       tpdigiPhiRingMean[55 - ie] /= denom;
     }
  
     for (unsigned iTT(0); iTT < EcalTrigTowerDetId::kSizeForDenseIndexing; ++iTT) {
       EcalTrigTowerDetId ttid(EcalTrigTowerDetId::detIdFromDenseIndex(iTT));
  
       float entries(sTPDigiThr.getBinContent(getEcalDQMSetupObjects(), ttid));
  
       unsigned index(ttid.ieta() < 0 ? ttid.ieta() + 28 : ttid.ieta() + 27);
  
       int quality(kGood);
  
       if (entries > minHits_ && entries > tpdigiPhiRingMean[index] * deviationThreshold_) {
         //        meHotTPDigiThr.fill(ttid);
         quality = kBad;
       }
  
       if (quality != kBad)
         continue;
  
       std::vector<DetId> ids(GetTrigTowerMap()->constituentsOf(ttid));
       for (unsigned iD(0); iD < ids.size(); ++iD) {
         DetId& id(ids[iD]);
  
         int quality(meQualitySummary.getBinContent(getEcalDQMSetupObjects(), id));
         if (quality == kMBad || quality == kBad)
           continue;
  
         meQualitySummary.setBinContent(
             getEcalDQMSetupObjects(),
             id,
             meQualitySummary.maskMatches(id, mask, statusManager_, GetTrigTowerMap()) ? kMBad : kBad);
       }
     }
  
     // Quality check: set entire FED to BAD if its occupancy begins to vanish
     // Fill FED statistics from (filtered) RecHit Occupancy
     float meanFEDEB(0), meanFEDEE(0), rmsFEDEB(0), rmsFEDEE(0);
     unsigned int nFEDEB(0), nFEDEE(0);
     for (unsigned iDCC(0); iDCC < nDCC; iDCC++) {
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh) {
         meanFEDEB += Nrhentries[iDCC];
         rmsFEDEB += Nrhentries[iDCC] * Nrhentries[iDCC];
         nFEDEB++;
       } else {
         meanFEDEE += Nrhentries[iDCC];
         rmsFEDEE += Nrhentries[iDCC] * Nrhentries[iDCC];
         nFEDEE++;
       }
     }
     meanFEDEB /= float(nFEDEB);
     rmsFEDEB /= float(nFEDEB);
     meanFEDEE /= float(nFEDEE);
     rmsFEDEE /= float(nFEDEE);
     rmsFEDEB = sqrt(abs(rmsFEDEB - meanFEDEB * meanFEDEB));
     rmsFEDEE = sqrt(abs(rmsFEDEE - meanFEDEE * meanFEDEE));
     // Analyze FED statistics
     float meanFED(0.), rmsFED(0.), nRMS(5.);
     for (MESet::iterator qsItr(meQualitySummary.beginChannel(GetElectronicsMap()));
          qsItr != meQualitySummary.end(GetElectronicsMap());
          qsItr.toNextChannel(GetElectronicsMap())) {
       DetId id(qsItr->getId());
       unsigned iDCC(dccId(id, GetElectronicsMap()) - 1);
       if (iDCC >= kEBmLow && iDCC <= kEBpHigh) {
         meanFED = meanFEDEB;
         rmsFED = rmsFEDEB;
       } else {
         meanFED = meanFEDEE;
         rmsFED = rmsFEDEE;
       }
       float threshold(meanFED < nRMS * rmsFED ? minHits_ : meanFED - nRMS * rmsFED);
       if (meanFED > 1000. && Nrhentries[iDCC] < threshold)
         meQualitySummary.setBinContent(
             getEcalDQMSetupObjects(),
             id,
             meQualitySummary.maskMatches(id, mask, statusManager_, GetTrigTowerMap()) ? kMBad : kBad);
     }
  
   }  // producePlots()

References funct::abs(), ecaldqm::MESetCollection::at(), ecaldqm::dccId(), makePileupJSON::denom, EcalTrigTowerDetId::detIdFromDenseIndex(), deviationThreshold_, EcalTriggerTower, dqmMemoryStats::float, ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), triggerObjects_cff::id, EcalTrigTowerDetId::ieta(), LEDCalibrationChannels::ieta, ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, EcalTrigTowerDetId::kSizeForDenseIndexing, ecaldqm::DQWorker::MEs_, minHits_, ecaldqm::nDCC, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_ERROR, EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING, quality, ecaldqm::scConstituents(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, remoteMonitoring_LED_IterMethod_cfg::threshold, ecaldqm::MESet::const_iterator::toNextChannel(), ecaldqm::MESet::iterator::toNextChannel(), and EcalTrigTowerConstituentsMap::towerOf().

◆ setParams()

void ecaldqm::OccupancyClient::setParams ( edm::ParameterSet const & _params )

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 16 of file OccupancyClient.cc.

                                                                 {
     minHits_ = _params.getUntrackedParameter<int>("minHits");
     deviationThreshold_ = _params.getUntrackedParameter<double>("deviationThreshold");
   }

References l1extraParticles_cfi::_params, deviationThreshold_, and minHits_.

Member Data Documentation

◆ deviationThreshold_

float ecaldqm::OccupancyClient::deviationThreshold_

private

Definition at line 19 of file OccupancyClient.h.

Referenced by producePlots(), and setParams().

◆ minHits_

int ecaldqm::OccupancyClient::minHits_