#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 &)

std::string const &	getName () const

bool	onlineMode () const

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setTime (time_t _t)

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

ecaldqm::PresampleClient::PresampleClient ( )

Definition at line 13 of file PresampleClient.cc.

References ecaldqm::DQWorkerClient::qualitySummaries_.

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

ecaldqm::PresampleClient::~PresampleClient ( )

inlineoverride

Definition at line 10 of file PresampleClient.h.

References producePlots(), and setParams().

10 {}

Member Function Documentation

void ecaldqm::PresampleClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 34 of file PresampleClient.cc.

Referenced by ~PresampleClient().

                                                 {
     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());
 
     MESet::const_iterator pItr(sPedestal);
     MESet::const_iterator pLSItr(sPedestalByLS);
     double maxEB(0.), minEB(0.), maxEE(0.), minEE(0.);
     double rmsMaxEB(0.), rmsMaxEE(0.);
     for (MESet::iterator qItr(meQuality.beginChannel()); qItr != qEnd; qItr.toNextChannel()) {
       pItr = qItr;
       pLSItr = qItr;
 
       DetId id(qItr->getId());
 
       bool doMask(meQuality.maskMatches(id, mask, statusManager_));
 
       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(id, doMask ? kMUnknown : kUnknown);
         meRMSMap.setBinContent(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));
 
       meMean.fill(dccid, mean);
       meRMS.fill(dccid, rms);
       meRMSMap.setBinContent(id, rms);
       meRMSMapAllByLumi.setBinContent(id, rmsLS);
 
       if (((mean > expectedMean_ + toleranceHigh_) || (mean < expectedMean_ - toleranceLow_)) || rms > rmsThresh) {
         qItr->setBinContent(doMask ? kMBad : kBad);
         meQualitySummary.setBinContent(id, doMask ? kMBad : kBad);
         if (!doMask)
           meErrorsSummary.fill(id);
       } else {
         qItr->setBinContent(doMask ? kMGood : kGood);
         meQualitySummary.setBinContent(id, doMask ? kMGood : kGood);
       }
 
       // Fill Presample Trend plots:
       // Use PedestalByLS which only contains digis from "current" LS
       float chStatus(sChStatus.getBinContent(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(EcalBarrel, double(timestamp_.iLumi), maxEB - minEB);
     meTrendMean.fill(EcalEndcap, double(timestamp_.iLumi), maxEE - minEE);
     meTrendRMS.fill(EcalBarrel, double(timestamp_.iLumi), rmsMaxEB);
     meTrendRMS.fill(EcalEndcap, double(timestamp_.iLumi), rmsMaxEE);
   }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 25 of file PresampleClient.cc.

References expectedMean_, edm::ParameterSet::getUntrackedParameter(), minChannelEntries_, toleranceHigh_, toleranceLow_, toleranceRMS_, and toleranceRMSFwd_.

Referenced by ~PresampleClient().

                                                                 {
     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");
   }