#include <PresampleClient.h>

Inheritance diagram for ecaldqm::PresampleClient:

Public Member Functions
	PresampleClient ()

void	producePlots (ProcessType) override

	~PresampleClient ()

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

Constructor & Destructor Documentation

ecaldqm::PresampleClient::PresampleClient ( )

Definition at line 14 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 ( )

inline

Definition at line 11 of file PresampleClient.h.

References producePlots(), and setParams().

11 {}

Member Function Documentation

void ecaldqm::PresampleClient::producePlots ( ProcessType )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 39 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 28 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");
   }