DQM/EcalMonitorClient/src/RawDataClient.cc - master - githubmirror/cmssw

DQM/EcalMonitorClient/src/RawDataClient.cc
76 строк3 KB
Abhirami Harilal
Change FEStatus quality threshold
11 июл 2024, 17:08
11 июл 2024, 17:085c43000
#include "DQM/EcalMonitorClient/interface/RawDataClient.h"

#include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"
#include "DQM/EcalCommon/interface/FEFlags.h"

#include "CondFormats/EcalObjects/interface/EcalDQMStatusHelper.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <cmath>

namespace ecaldqm {

  RawDataClient::RawDataClient() : DQWorkerClient(), minEvents_(0), synchErrThresholdFactor_(0.) {
    qualitySummaries_.insert("QualitySummary");
  }

  void RawDataClient::setParams(edm::ParameterSet const& _params) {
    minEvents_ = _params.getUntrackedParameter<int>("minEvents");
    synchErrThresholdFactor_ = _params.getUntrackedParameter<double>("synchErrThresholdFactor");
  }

  void RawDataClient::producePlots(ProcessType) {
    MESet& meQualitySummary(MEs_.at("QualitySummary"));
    MESet& meErrorsSummary(MEs_.at("ErrorsSummary"));

    MESet const& sEntries(sources_.at("Entries"));
    MESet const& sL1ADCC(sources_.at("L1ADCC"));
    MESet const& sFEStatus(sources_.at("FEStatus"));

    uint32_t mask(1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR);

    std::vector<int> dccStatus(nDCC, 1);

    for (unsigned iDCC(0); iDCC < nDCC; ++iDCC) {
      double entries(sEntries.getBinContent(getEcalDQMSetupObjects(), iDCC + 1));
      if (entries > 1. && sL1ADCC.getBinContent(getEcalDQMSetupObjects(), iDCC + 1) >
                              synchErrThresholdFactor_ * std::log(entries) / std::log(10.))
        dccStatus[iDCC] = 0;
    }

    MESet::iterator meEnd(meQualitySummary.end(GetElectronicsMap()));
    for (MESet::iterator meItr(meQualitySummary.beginChannel(GetElectronicsMap())); meItr != meEnd;
         meItr.toNextChannel(GetElectronicsMap())) {
      DetId id(meItr->getId());

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

      int dccid(dccId(id, GetElectronicsMap()));

      if (dccStatus[dccid - 1] == 0) {
        meItr->setBinContent(doMask ? kMUnknown : kUnknown);
        continue;
      }

      int towerStatus(doMask ? kMGood : kGood);
      float towerEntries(0.);
      for (unsigned iS(0); iS < nFEFlags; iS++) {
        float entries(sFEStatus.getBinContent(getEcalDQMSetupObjects(), id, iS + 1));
        towerEntries += entries;
        if (entries > minEvents_ && iS != Enabled && iS != Suppressed && iS != ForcedFullSupp && iS != FIFOFull &&
            iS != ForcedZS)
          towerStatus = doMask ? kMBad : kBad;
      }

      if (towerEntries < 1.)
        towerStatus = doMask ? kMUnknown : kUnknown;

      meItr->setBinContent(towerStatus);
      if (towerStatus == kBad)
        meErrorsSummary.fill(getEcalDQMSetupObjects(), dccid);
    }
  }

  DEFINE_ECALDQM_WORKER(RawDataClient);
}  // namespace ecaldqm