ecaldqm::PedestalWriter Class Reference

#include <DBWriterWorkers.h>

Inheritance diagram for ecaldqm::PedestalWriter:

Public Member Functions
	PedestalWriter (edm::ParameterSet const &)
bool	run (EcalCondDBInterface *, MonRunIOV &) override
	~PedestalWriter () override
Public Member Functions inherited from ecaldqm::DBWriterWorker
	DBWriterWorker (std::string const &, edm::ParameterSet const &)
std::string const &	getName () const
bool	isActive () const
void	retrieveSource (DQMStore::IGetter &)
bool	runsOn (std::string const &_runType) const
void	setVerbosity (int _v)
virtual	~DBWriterWorker ()

Private Attributes
std::map< int, unsigned >	gainToME_
std::map< int, unsigned >	pnGainToME_

Additional Inherited Members
Protected Attributes inherited from ecaldqm::DBWriterWorker
bool	active_
std::string const	name_
std::set< std::string >	runTypes_
MESetCollection	source_
int	verbosity_

Detailed Description

Definition at line 60 of file DBWriterWorkers.h.

Constructor & Destructor Documentation

ecaldqm::PedestalWriter::PedestalWriter

(

edm::ParameterSet const &

_ps

)

Definition at line 572 of file DBWriterWorkers.cc.

References Exception, muonCSCDigis_cfi::gain, gainToME_, ecaldqm::MESetMulti::getIndex(), edm::ParameterSet::getUntrackedParameter(), muonCSCDigis_cfi::pedestal, pnGainToME_, and ecaldqm::DBWriterWorker::source_.

      : DBWriterWorker("Pedestal", _ps), gainToME_(), pnGainToME_() {
    std::vector<int> MGPAGains(_ps.getUntrackedParameter<std::vector<int>>("MGPAGains"));
    std::vector<int> MGPAGainsPN(_ps.getUntrackedParameter<std::vector<int>>("MGPAGainsPN"));

    MESet::PathReplacements repl;

    MESetMulti const &pedestal(static_cast<MESetMulti const &>(source_.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 &>(source_.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);
    }
  }

ecaldqm::PedestalWriter::~PedestalWriter ( )

inlineoverride

Definition at line 63 of file DBWriterWorkers.h.

References ecaldqm::DBWriterWorker::run().

{}

Member Function Documentation

bool ecaldqm::PedestalWriter::run ( EcalCondDBInterface *  _db, MonRunIOV &  _iov  )

overridevirtual

Implements ecaldqm::DBWriterWorker.

Definition at line 602 of file DBWriterWorkers.cc.

References ecaldqm::crystalID(), data, MillePedeFileConverter_cfg::e, EcalBarrel, EcalEndcap, Exception, spr::find(), muonCSCDigis_cfi::gain, gainToME_, EcalCondDBInterface::insertDataArraySet(), ecaldqm::kBad, ecaldqm::kEEpLow, ecaldqm::kMBad, ecaldqm::lmPNID(), SiStripPI::mean, ecaldqm::memDCC, pnGainToME_, ecaldqm::qualityOK(), mps_fire::result, SiStripPI::rms, MonPedestalsDat::setPedMeanG1(), MonPNPedDat::setPedMeanG1(), MonPedestalsDat::setPedMeanG12(), MonPNPedDat::setPedMeanG16(), MonPedestalsDat::setPedMeanG6(), MonPNPedDat::setPedRMSG1(), MonPedestalsDat::setPedRMSG1(), MonPedestalsDat::setPedRMSG12(), MonPNPedDat::setPedRMSG16(), MonPedestalsDat::setPedRMSG6(), MonPNPedDat::setTaskStatus(), MonPedestalsDat::setTaskStatus(), ecaldqm::DBWriterWorker::source_, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, and ecaldqm::MESet::const_iterator::toNextChannel().

Referenced by Types.EventID::cppID(), and Types.LuminosityBlockID::cppID().

                                                                    {
    /*
    uses
    PedestalTask.Pedestal (h01, h02, h03)
    PedestalTask.PNPedestal (i01, i02)
    PedestalClient.Quality (meg01, meg02, meg03)
    PedestalClient.PNQualitySummary (meg04, meg05)
  */

    bool result(true);

    std::map<EcalLogicID, MonPedestalsDat> pedestals;
    std::map<EcalLogicID, MonPNPedDat> pnPedestals;

    MESet const &pedestalME(source_.at("Pedestal"));
    MESet const &qualityME(source_.at("Quality"));

    MESet const &pnPedestalME(source_.at("PNPedestal"));
    MESet const &pnQualityME(source_.at("PNQuality"));

    for (std::map<int, unsigned>::iterator gainItr(gainToME_.begin()); gainItr != gainToME_.end(); ++gainItr) {
      int gain(gainItr->first);
      int iM(gainItr->second);

      static_cast<MESetMulti const &>(pedestalME).use(iM);
      static_cast<MESetMulti const &>(qualityME).use(iM);

      MESet::const_iterator pEnd(pedestalME.end());
      MESet::const_iterator qItr(qualityME);
      for (MESet::const_iterator pItr(pedestalME.beginChannel()); pItr != pEnd; pItr.toNextChannel()) {
        float entries(pItr->getBinEntries());
        if (entries < 1.)
          continue;

        qItr = pItr;

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

        EcalLogicID logicID(crystalID(pItr->getId()));
        if (pedestals.find(logicID) == pedestals.end()) {
          MonPedestalsDat &insertion(pedestals[logicID]);
          insertion.setPedMeanG1(-1.);
          insertion.setPedRMSG1(-1.);
          insertion.setPedMeanG6(-1.);
          insertion.setPedRMSG6(-1.);
          insertion.setPedMeanG12(-1.);
          insertion.setPedRMSG12(-1.);
          insertion.setTaskStatus(false);
        }

        MonPedestalsDat &data(pedestals[logicID]);
        switch (gain) {
          case 1:
            data.setPedMeanG1(mean);
            data.setPedRMSG1(rms);
            break;
          case 6:
            data.setPedMeanG6(mean);
            data.setPedRMSG6(rms);
            break;
          case 12:
            data.setPedMeanG12(mean);
            data.setPedRMSG12(rms);
            break;
        }

        int channelStatus(qItr->getBinContent());
        bool channelBad(channelStatus == kBad || channelStatus == kMBad);
        if (channelBad)
          data.setTaskStatus(true);

        result &= qualityOK(channelStatus);
      }
    }

    for (std::map<int, unsigned>::iterator gainItr(pnGainToME_.begin()); gainItr != pnGainToME_.end(); ++gainItr) {
      int gain(gainItr->first);
      int iM(gainItr->second);

      static_cast<MESetMulti const &>(pnPedestalME).use(iM);
      static_cast<MESetMulti const &>(pnQualityME).use(iM);

      for (unsigned iMD(0); iMD < memDCC.size(); ++iMD) {
        unsigned iDCC(memDCC[iMD]);

        int subdet(iDCC <= kEEmHigh || iDCC >= kEEpLow ? EcalEndcap : EcalBarrel);

        for (unsigned iPN(1); iPN <= 10; ++iPN) {
          EcalPnDiodeDetId pnid(subdet, iDCC + 1, iPN);

          float entries(pnPedestalME.getBinEntries(pnid));
          if (entries < 1.)
            continue;

          float mean(pnPedestalME.getBinContent(pnid));
          float rms(pnPedestalME.getBinError(pnid) * std::sqrt(entries));

          EcalLogicID logicID(lmPNID(pnid));
          if (pnPedestals.find(logicID) == pnPedestals.end()) {
            MonPNPedDat &insertion(pnPedestals[logicID]);
            insertion.setPedMeanG1(-1.);
            insertion.setPedRMSG1(-1.);
            insertion.setPedMeanG16(-1.);
            insertion.setPedRMSG16(-1.);
            insertion.setTaskStatus(false);
          }

          MonPNPedDat &data(pnPedestals[lmPNID(pnid)]);
          switch (gain) {
            case 1:
              data.setPedMeanG1(mean);
              data.setPedRMSG1(rms);
              break;
            case 16:
              data.setPedMeanG16(mean);
              data.setPedRMSG16(rms);
              break;
          }

          int channelStatus(pnQualityME.getBinContent(pnid));
          bool channelBad(channelStatus == kBad || channelStatus == kMBad);
          if (channelBad)
            data.setTaskStatus(true);

          result &= qualityOK(channelStatus);
        }
      }
    }

    try {
      if (!pedestals.empty())
        _db->insertDataArraySet(&pedestals, &_iov);
      if (!pnPedestals.empty())
        _db->insertDataArraySet(&pnPedestals, &_iov);
    } catch (std::runtime_error &e) {
      if (std::string(e.what()).find("unique constraint") != std::string::npos)
        edm::LogWarning("EcalDQM") << e.what();
      else
        throw cms::Exception("DBError") << e.what();
    }

    return result;
  }

Member Data Documentation

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

private

Definition at line 68 of file DBWriterWorkers.h.

Referenced by PedestalWriter(), and run().

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

private

Definition at line 69 of file DBWriterWorkers.h.

Referenced by PedestalWriter(), and run().