ecaldqm::PresampleClient Class Reference

#include <PresampleClient.h>

Inheritance diagram for ecaldqm::PresampleClient:

Public Member Functions
	PresampleClient ()
void	producePlots (ProcessType) override
	~PresampleClient () 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	expectedMean_
int	minChannelEntries_
float	toleranceHigh_
float	toleranceLow_
float	toleranceRMS_
float	toleranceRMSFwd_

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

Constructor & Destructor Documentation

PresampleClient()

ecaldqm::PresampleClient::PresampleClient

(

)

Definition at line 13 of file PresampleClient.cc.

      : DQWorkerClient(),
        minChannelEntries_(0),
        expectedMean_(0.),
        toleranceLow_(0.),
        toleranceHigh_(0.),
        toleranceRMS_(0.),
        toleranceRMSFwd_(0.) {
    qualitySummaries_.insert("Quality");
    qualitySummaries_.insert("QualitySummary");
  }

References ecaldqm::DQWorkerClient::qualitySummaries_.

~PresampleClient()

ecaldqm::PresampleClient::~PresampleClient ( )

inlineoverride

Definition at line 10 of file PresampleClient.h.

{}

Member Function Documentation

producePlots()

void ecaldqm::PresampleClient::producePlots ( ProcessType  )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 34 of file PresampleClient.cc.

                                                {
    MESet& meQualitySummary(MEs_.at("QualitySummary"));
    MESet& meQuality(MEs_.at("Quality"));
    MESet& meErrorsSummary(MEs_.at("ErrorsSummary"));
    MESet& meMean(MEs_.at("Mean"));
    MESet& meRMS(MEs_.at("RMS"));
    MESet& meRMSMap(MEs_.at("RMSMap"));
    MESet& meRMSMapAll(MEs_.at("RMSMapAll"));
    MESet& meRMSMapAllByLumi(MEs_.at("RMSMapAllByLumi"));
 
    MESet const& sPedestal(sources_.at("Pedestal"));
    MESet const& sPedestalByLS(sources_.at("PedestalByLS"));
    MESet const& sChStatus(sources_.at("ChStatus"));
 
    uint32_t mask(1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR |
                  1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR);
 
    MESet::iterator qEnd(meQuality.end(GetElectronicsMap()));
 
    MESet::const_iterator pItr(GetElectronicsMap(), sPedestal);
    MESet::const_iterator pLSItr(GetElectronicsMap(), sPedestalByLS);
    double maxEB(0.), minEB(0.), maxEE(0.), minEE(0.);
    double rmsMaxEB(0.), rmsMaxEE(0.);
    for (MESet::iterator qItr(meQuality.beginChannel(GetElectronicsMap())); qItr != qEnd;
         qItr.toNextChannel(GetElectronicsMap())) {
      pItr = qItr;
      pLSItr = qItr;
 
      DetId id(qItr->getId());
 
      bool doMask(meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));
 
      double rmsThresh(toleranceRMS_);
 
      if (isForward(id))
        rmsThresh = toleranceRMSFwd_;
 
      double entries(pItr->getBinEntries());
      double entriesLS(pLSItr->getBinEntries());
 
      if (entries < minChannelEntries_) {
        qItr->setBinContent(doMask ? kMUnknown : kUnknown);
        meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMUnknown : kUnknown);
        meRMSMap.setBinContent(getEcalDQMSetupObjects(), id, -1.);
        continue;
      }
 
      double mean(pItr->getBinContent());
      double meanLS(pLSItr->getBinContent());
      double rms(pItr->getBinError() * std::sqrt(entries));
      double rmsLS(pLSItr->getBinError() * std::sqrt(entriesLS));
 
      int dccid(dccId(id, GetElectronicsMap()));
 
      meMean.fill(getEcalDQMSetupObjects(), dccid, mean);
      meRMS.fill(getEcalDQMSetupObjects(), dccid, rms);
      meRMSMap.setBinContent(getEcalDQMSetupObjects(), id, rms);
      meRMSMapAllByLumi.setBinContent(getEcalDQMSetupObjects(), id, rmsLS);
 
      if (((mean > expectedMean_ + toleranceHigh_) || (mean < expectedMean_ - toleranceLow_)) || rms > rmsThresh) {
        qItr->setBinContent(doMask ? kMBad : kBad);
        meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMBad : kBad);
        if (!doMask)
          meErrorsSummary.fill(getEcalDQMSetupObjects(), id);
      } else {
        qItr->setBinContent(doMask ? kMGood : kGood);
        meQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMGood : kGood);
      }
 
      // Fill Presample Trend plots:
      // Use PedestalByLS which only contains digis from "current" LS
      float chStatus(sChStatus.getBinContent(getEcalDQMSetupObjects(), id));
      if (entriesLS < minChannelEntries_)
        continue;
      if (chStatus != EcalChannelStatusCode::kOk)
        continue;  // exclude problematic channels
 
      // Get max/min
      // Min is effectively just 0
      if (id.subdetId() == EcalBarrel) {
        if (meanLS > maxEB)
          maxEB = meanLS;
        if (meanLS < minEB)
          minEB = meanLS;
        if (rmsLS > rmsMaxEB)
          rmsMaxEB = rmsLS;
      } else {
        if (meanLS > maxEE)
          maxEE = meanLS;
        if (meanLS < minEE)
          minEE = meanLS;
        if (rmsLS > rmsMaxEE)
          rmsMaxEE = rmsLS;
      }
 
    }  // qItr
 
    towerAverage_(meRMSMapAll, meRMSMap, -1.);
 
    MESet& meTrendMean(MEs_.at("TrendMean"));
    MESet& meTrendRMS(MEs_.at("TrendRMS"));
    meTrendMean.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), maxEB - minEB);
    meTrendMean.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), maxEE - minEE);
    meTrendRMS.fill(getEcalDQMSetupObjects(), EcalBarrel, double(timestamp_.iLumi), rmsMaxEB);
    meTrendRMS.fill(getEcalDQMSetupObjects(), EcalEndcap, double(timestamp_.iLumi), rmsMaxEE);
  }

References ecaldqm::MESetCollection::at(), ecaldqm::dccId(), EcalBarrel, EcalEndcap, expectedMean_, ecaldqm::MESet::fill(), ecaldqm::MESet::ConstBin::getBinContent(), ecaldqm::MESet::ConstBin::getBinEntries(), ecaldqm::MESet::ConstBin::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), triggerObjects_cff::id, ecaldqm::DQWorker::Timestamp::iLumi, ecaldqm::isForward(), ecaldqm::DQWorkerClient::kBad, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, EcalChannelStatusCode::kOk, ecaldqm::DQWorkerClient::kUnknown, SiStripPI::mean, ecaldqm::DQWorker::MEs_, minChannelEntries_, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR, SiStripPI::rms, ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, ecaldqm::DQWorker::timestamp_, toleranceHigh_, toleranceLow_, toleranceRMS_, toleranceRMSFwd_, ecaldqm::MESet::iterator::toNextChannel(), and ecaldqm::DQWorkerClient::towerAverage_().

setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 25 of file PresampleClient.cc.

                                                                {
    minChannelEntries_ = _params.getUntrackedParameter<int>("minChannelEntries");
    expectedMean_ = _params.getUntrackedParameter<double>("expectedMean");
    toleranceLow_ = _params.getUntrackedParameter<double>("toleranceLow");
    toleranceHigh_ = _params.getUntrackedParameter<double>("toleranceHigh");
    toleranceRMS_ = _params.getUntrackedParameter<double>("toleranceRMS");
    toleranceRMSFwd_ = _params.getUntrackedParameter<double>("toleranceRMSFwd");
  }

References l1extraParticles_cfi::_params, expectedMean_, minChannelEntries_, toleranceHigh_, toleranceLow_, toleranceRMS_, and toleranceRMSFwd_.

Member Data Documentation

expectedMean_

float ecaldqm::PresampleClient::expectedMean_

private

Definition at line 18 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

minChannelEntries_

int ecaldqm::PresampleClient::minChannelEntries_

private

Definition at line 17 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

toleranceHigh_

float ecaldqm::PresampleClient::toleranceHigh_

private

Definition at line 20 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

toleranceLow_

float ecaldqm::PresampleClient::toleranceLow_

private

Definition at line 19 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

toleranceRMS_

float ecaldqm::PresampleClient::toleranceRMS_

private

Definition at line 21 of file PresampleClient.h.

Referenced by producePlots(), and setParams().

toleranceRMSFwd_

float ecaldqm::PresampleClient::toleranceRMSFwd_

private

Definition at line 22 of file PresampleClient.h.

Referenced by producePlots(), and setParams().