ecaldqm::PedestalClient Class Reference

#include <PedestalClient.h>

Inheritance diagram for ecaldqm::PedestalClient:

Public Member Functions
	PedestalClient ()
void	producePlots (ProcessType) override
	~PedestalClient () 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)
virtual void	setTokens (edm::ConsumesCollector &)
	~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 &)
bool	checkElectronicsMap (bool=true)
bool	checkGeometry (bool=true)
bool	checkTopology (bool=true)
bool	checkTrigTowerMap (bool=true)
	DQWorker ()
virtual void	endRun (edm::Run const &, edm::EventSetup const &)
EcalDQMSetupObjects const	getEcalDQMSetupObjects ()
EcalElectronicsMapping const *	GetElectronicsMap ()
CaloGeometry const *	GetGeometry ()
std::string const &	getName () const
CaloTopology const *	GetTopology ()
EcalTrigTowerConstituentsMap const *	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	setSetupObjectsEndLumi (edm::EventSetup const &)
void	setTime (time_t _t)
void	setTokens (edm::ConsumesCollector &)
virtual	~DQWorker () noexcept(false)

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

Private Attributes
float	expectedMean_
float	expectedPNMean_
std::map< int, unsigned >	gainToME_
int	minChannelEntries_
std::map< int, unsigned >	pnGainToME_
float	toleranceMean_
float	tolerancePNMean_
std::vector< float >	tolerancePNRMS_
std::vector< float >	toleranceRMSEB_
std::vector< float >	toleranceRMSEE_

Additional Inherited Members
Public Types inherited from ecaldqm::DQWorkerClient
enum	ProcessType { kLumi, kJob, nProcessType }
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)
Public Attributes inherited from ecaldqm::DQWorker
edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	elecMapHandle
edm::ESGetToken< EcalElectronicsMapping, EcalMappingRcd >	elecMapHandleEndLumi
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandle
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geomHandleEndLumi
edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandle
edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topoHandleEndLumi
edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandle
edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapHandleEndLumi
Static Public Attributes inherited from ecaldqm::DQWorkerClient
static constexpr int	kBad = 0
static constexpr int	kGood = 1
static constexpr int	kMBad = 3
static constexpr int	kMGood = 4
static constexpr int	kMUnknown = 5
static constexpr int	kUnknown = 2
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 PedestalClient.h.

Constructor & Destructor Documentation

PedestalClient()

ecaldqm::PedestalClient::PedestalClient

(

)

Definition at line 15 of file PedestalClient.cc.

      : DQWorkerClient(),
        gainToME_(),
        pnGainToME_(),
        minChannelEntries_(0),
        expectedMean_(0.),
        toleranceMean_(0.),
        toleranceRMSEB_(0),
        toleranceRMSEE_(0),
        expectedPNMean_(0.),
        tolerancePNMean_(0.),
        tolerancePNRMS_(0) {}

~PedestalClient()

ecaldqm::PedestalClient::~PedestalClient ( )

inlineoverride

Definition at line 10 of file PedestalClient.h.

{}

Member Function Documentation

producePlots()

void ecaldqm::PedestalClient::producePlots ( ProcessType  )

overridevirtual

Implements ecaldqm::DQWorkerClient.

Definition at line 111 of file PedestalClient.cc.

References funct::abs(), ecaldqm::MESetCollection::at(), ecaldqm::MESetMulti::beginChannel(), EcalBarrel, EcalEndcap, ecaldqm::MESetMulti::end(), expectedMean_, expectedPNMean_, ecaldqm::MESetMulti::fill(), gainToME_, ecaldqm::MESetMulti::getBinContent(), ecaldqm::MESetMulti::getBinEntries(), ecaldqm::MESetMulti::getBinError(), ecaldqm::DQWorker::getEcalDQMSetupObjects(), ecaldqm::DQWorker::GetElectronicsMap(), ecaldqm::DQWorker::GetTrigTowerMap(), l1ctLayer2EG_cff::id, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), ecaldqm::DQWorkerClient::kBad, ecaldqm::kEBmLow, ecaldqm::kEBpHigh, ecaldqm::DQWorkerClient::kGood, ecaldqm::DQWorkerClient::kMBad, ecaldqm::DQWorkerClient::kMGood, ecaldqm::DQWorkerClient::kMUnknown, ecaldqm::DQWorkerClient::kUnknown, ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::mask, ecaldqm::MESetMulti::maskMatches(), SiStripPI::mean, ecaldqm::memDCCIndex(), ecaldqm::DQWorker::MEs_, minChannelEntries_, ecaldqm::nDCC, EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR, EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR, EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR, pnGainToME_, SiStripPI::rms, ecaldqm::MESetMulti::setBinContent(), ecaldqm::DQWorkerClient::sources_, mathSSE::sqrt(), ecaldqm::DQWorkerClient::statusManager_, toleranceMean_, tolerancePNMean_, tolerancePNRMS_, toleranceRMSEB_, toleranceRMSEE_, ecaldqm::MESet::iterator::toNextChannel(), ecaldqm::DQWorkerClient::towerAverage_(), and ecaldqm::MESetMulti::use().

                                               {
    using namespace std;

    MESetMulti& meQuality(static_cast<MESetMulti&>(MEs_.at("Quality")));
    MESetMulti& meQualitySummary(static_cast<MESetMulti&>(MEs_.at("QualitySummary")));
    MESetMulti& meMean(static_cast<MESetMulti&>(MEs_.at("Mean")));
    MESetMulti& meRMS(static_cast<MESetMulti&>(MEs_.at("RMS")));
    MESetMulti& mePNQualitySummary(static_cast<MESetMulti&>(MEs_.at("PNQualitySummary")));
    MESetMulti& mePNRMS(static_cast<MESetMulti&>(MEs_.at("PNRMS")));

    MESetMulti const& sPedestal(static_cast<MESetMulti const&>(sources_.at("Pedestal")));
    MESetMulti const& sPNPedestal(static_cast<MESetMulti const&>(sources_.at("PNPedestal")));

    for (map<int, unsigned>::iterator gainItr(gainToME_.begin()); gainItr != gainToME_.end(); ++gainItr) {
      meQuality.use(gainItr->second);
      meQualitySummary.use(gainItr->second);
      meMean.use(gainItr->second);
      meRMS.use(gainItr->second);

      sPedestal.use(gainItr->second);

      uint32_t mask(0);
      switch (gainItr->first) {
        case 1:
          mask |= (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR);
          break;
        case 6:
          mask |= (1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR);
          break;
        case 12:
          mask |= (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR);
          break;
        default:
          break;
      }

      MESet::iterator qEnd(meQuality.end(GetElectronicsMap()));
      MESet::const_iterator pItr(GetElectronicsMap(), sPedestal);
      for (MESet::iterator qItr(meQuality.beginChannel(GetElectronicsMap())); qItr != qEnd;
           qItr.toNextChannel(GetElectronicsMap())) {
        DetId id(qItr->getId());

        bool doMask(meQuality.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));

        pItr = qItr;

        float entries(pItr->getBinEntries());

        if (entries < minChannelEntries_) {
          qItr->setBinContent(doMask ? kMUnknown : kUnknown);
          continue;
        }

        float mean(pItr->getBinContent());
        float rms(pItr->getBinError() * sqrt(entries));

        meMean.fill(getEcalDQMSetupObjects(), id, mean);
        meRMS.fill(getEcalDQMSetupObjects(), id, rms);

        float toleranceRMS_ =
            (id.subdetId() == EcalBarrel) ? toleranceRMSEB_[gainItr->second] : toleranceRMSEE_[gainItr->second];

        if (abs(mean - expectedMean_) > toleranceMean_ || rms > toleranceRMS_)
          qItr->setBinContent(doMask ? kMBad : kBad);
        else
          qItr->setBinContent(doMask ? kMGood : kGood);
      }

      towerAverage_(meQualitySummary, meQuality, 0.2);
    }

    for (map<int, unsigned>::iterator gainItr(pnGainToME_.begin()); gainItr != pnGainToME_.end(); ++gainItr) {
      mePNQualitySummary.use(gainItr->second);
      mePNRMS.use(gainItr->second);

      sPNPedestal.use(gainItr->second);

      uint32_t mask(0);
      switch (gainItr->first) {
        case 1:
          mask |= (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR);
          break;
        case 16:
          mask |= (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR |
                   1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR);
          break;
        default:
          break;
      }

      for (unsigned iDCC(0); iDCC < nDCC; ++iDCC) {
        if (memDCCIndex(iDCC + 1) == unsigned(-1))
          continue;

        for (unsigned iPN(0); iPN < 10; ++iPN) {
          int subdet(0);
          if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
            subdet = EcalBarrel;
          else
            subdet = EcalEndcap;

          EcalPnDiodeDetId id(subdet, iDCC + 1, iPN + 1);

          bool doMask(mePNQualitySummary.maskMatches(id, mask, statusManager_, GetTrigTowerMap()));

          float entries(sPNPedestal.getBinEntries(getEcalDQMSetupObjects(), id));

          if (entries < minChannelEntries_) {
            mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMUnknown : kUnknown);
            continue;
          }

          float mean(sPNPedestal.getBinContent(getEcalDQMSetupObjects(), id));
          float rms(sPNPedestal.getBinError(getEcalDQMSetupObjects(), id) * sqrt(entries));

          mePNRMS.fill(getEcalDQMSetupObjects(), id, rms);

          if (abs(mean - expectedPNMean_) > tolerancePNMean_ || rms > tolerancePNRMS_[gainItr->second])
            mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMBad : kBad);
          else
            mePNQualitySummary.setBinContent(getEcalDQMSetupObjects(), id, doMask ? kMGood : kGood);
        }
      }
    }
  }

setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 28 of file PedestalClient.cc.

References l1extraParticles_cfi::_params, ecaldqm::MESetCollection::at(), Exception, expectedMean_, expectedPNMean_, PedestalClient_cfi::gain, gainToME_, ecaldqm::MESetMulti::getIndex(), PedestalClient_cfi::MGPAGains, PedestalClient_cfi::MGPAGainsPN, minChannelEntries_, EcalCondDBWriter_cfi::pedestal, pnGainToME_, ecaldqm::DQWorkerClient::qualitySummaries_, ecaldqm::DQWorkerClient::sources_, to_string(), toleranceMean_, tolerancePNMean_, tolerancePNRMS_, toleranceRMSEB_, and toleranceRMSEE_.

                                                               {
    minChannelEntries_ = _params.getUntrackedParameter<int>("minChannelEntries");
    expectedMean_ = _params.getUntrackedParameter<double>("expectedMean");
    toleranceMean_ = _params.getUntrackedParameter<double>("toleranceMean");
    expectedPNMean_ = _params.getUntrackedParameter<double>("expectedPNMean");
    tolerancePNMean_ = _params.getUntrackedParameter<double>("tolerancePNMean");

    std::vector<int> MGPAGains(_params.getUntrackedParameter<std::vector<int> >("MGPAGains"));
    std::vector<int> MGPAGainsPN(_params.getUntrackedParameter<std::vector<int> >("MGPAGainsPN"));

    MESet::PathReplacements repl;

    MESetMulti const& pedestal(static_cast<MESetMulti const&>(sources_.at("Pedestal")));
    unsigned nG(MGPAGains.size());
    for (unsigned iG(0); iG != nG; ++iG) {
      int gain(MGPAGains[iG]);
      if (gain != 1 && gain != 6 && gain != 12)
        throw cms::Exception("InvalidConfiguration") << "MGPA gain";
      repl["gain"] = std::to_string(gain);
      gainToME_[gain] = pedestal.getIndex(repl);
    }

    repl.clear();

    MESetMulti const& pnPedestal(static_cast<MESetMulti const&>(sources_.at("PNPedestal")));
    unsigned nGPN(MGPAGainsPN.size());
    for (unsigned iG(0); iG != nGPN; ++iG) {
      int gain(MGPAGainsPN[iG]);
      if (gain != 1 && gain != 16)
        throw cms::Exception("InvalidConfiguration") << "PN MGPA gain";
      repl["pngain"] = std::to_string(gain);
      pnGainToME_[gain] = pnPedestal.getIndex(repl);
    }

    toleranceRMSEB_.resize(nG);
    toleranceRMSEE_.resize(nG);

    std::vector<double> inToleranceRMSEB(_params.getUntrackedParameter<std::vector<double> >("toleranceRMSEB"));
    std::vector<double> inToleranceRMSEE(_params.getUntrackedParameter<std::vector<double> >("toleranceRMSEE"));

    for (std::map<int, unsigned>::iterator gainItr(gainToME_.begin()); gainItr != gainToME_.end(); ++gainItr) {
      unsigned iME(gainItr->second);
      unsigned iGain(0);
      switch (gainItr->first) {
        case 1:
          iGain = 0;
          break;
        case 6:
          iGain = 1;
          break;
        case 12:
          iGain = 2;
          break;
      }

      toleranceRMSEB_[iME] = inToleranceRMSEB[iGain];
      toleranceRMSEE_[iME] = inToleranceRMSEE[iGain];
    }

    tolerancePNRMS_.resize(nGPN);

    std::vector<double> inTolerancePNRMS(_params.getUntrackedParameter<std::vector<double> >("tolerancePNRMS"));

    for (std::map<int, unsigned>::iterator gainItr(pnGainToME_.begin()); gainItr != pnGainToME_.end(); ++gainItr) {
      unsigned iME(gainItr->second);
      unsigned iGain(0);
      switch (gainItr->first) {
        case 1:
          iGain = 0;
          break;
        case 16:
          iGain = 1;
          break;
      }

      tolerancePNRMS_[iME] = inTolerancePNRMS[iGain];
    }

    qualitySummaries_.insert("Quality");
    qualitySummaries_.insert("QualitySummary");
    qualitySummaries_.insert("PNQualitySummary");
  }

Member Data Documentation

expectedMean_

float ecaldqm::PedestalClient::expectedMean_

private

Definition at line 21 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

expectedPNMean_

float ecaldqm::PedestalClient::expectedPNMean_

private

Definition at line 25 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

gainToME_

std::map<int, unsigned> ecaldqm::PedestalClient::gainToME_

private

Definition at line 17 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

minChannelEntries_

int ecaldqm::PedestalClient::minChannelEntries_

private

Definition at line 20 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

pnGainToME_

std::map<int, unsigned> ecaldqm::PedestalClient::pnGainToME_

private

Definition at line 18 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

toleranceMean_

float ecaldqm::PedestalClient::toleranceMean_

private

Definition at line 22 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

tolerancePNMean_

float ecaldqm::PedestalClient::tolerancePNMean_

private

Definition at line 26 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

tolerancePNRMS_

std::vector<float> ecaldqm::PedestalClient::tolerancePNRMS_

private

Definition at line 27 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

toleranceRMSEB_

std::vector<float> ecaldqm::PedestalClient::toleranceRMSEB_

private

Definition at line 23 of file PedestalClient.h.

Referenced by producePlots(), and setParams().

toleranceRMSEE_

std::vector<float> ecaldqm::PedestalClient::toleranceRMSEE_

private

Definition at line 24 of file PedestalClient.h.

Referenced by producePlots(), and setParams().