DQM/HcalCommon/src/DQClient.cc - master - githubmirror/cmssw

DQM/HcalCommon/src/DQClient.cc
139 строк6 KB
Long
code modification upon request
03 апр 2026, 06:28
03 апр 2026, 06:286348b62
#include "DQM/HcalCommon/interface/DQClient.h"

namespace hcaldqm {
  using namespace constants;
  DQClient::DQClient(std::string const &name,
                     std::string const &taskname,
                     edm::ParameterSet const &ps,
                     edm::ConsumesCollector &iC)
      : DQModule(ps),
        _taskname(taskname),
        _totalLS(0),
        _maxProcessedLS(0),
        hcalDbServiceToken_(iC.esConsumes<HcalDbService, HcalDbRecord, edm::Transition::BeginRun>()),
        runInfoToken_(iC.esConsumes<RunInfo, RunInfoRcd, edm::Transition::BeginRun>()),
        hcalChannelQualityToken_(
            iC.esConsumes<HcalChannelQuality, HcalChannelQualityRcd, edm::Transition::BeginLuminosityBlock>(
                edm::ESInputTag("", "withTopo"))) {
    //	- SET THE TASK NAME YOU REFER TO
    //	- SET THE CLIENT'S NAME AS WELL - RUN SUMMARY PLOTS
    //	WILL BE GENERATED UNDER THAT FOLDER
    _name = name;
  }

  void DQClient::beginRun(edm::Run const &r, edm::EventSetup const &es) {
    //	TEMPORARY
    _vhashFEDs.clear();
    _vFEDs.clear();
    _vcdaqEids.clear();
    _vhashCrates.clear();

    //	get various FED lists
    edm::ESHandle<HcalDbService> dbs = es.getHandle(hcalDbServiceToken_);
    _emap = dbs->getHcalMapping();

    if (_ptype != fOffline) {  // hidefed2crate
      _vFEDs = utilities::getFEDList(_emap);
      for (std::vector<int>::const_iterator it = _vFEDs.begin(); it != _vFEDs.end(); ++it) {
        //
        //	FIXME
        //	until there exists a map of FED2Crate and Crate2FED,
        //	all the unknown Crates will be mapped to 0...
        //
        if (*it == 0) {
          _vhashFEDs.push_back(HcalElectronicsId(0, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
          continue;
        }

        if (*it > FED_VME_MAX) {
          std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
          _vhashFEDs.push_back(
              HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        } else
          _vhashFEDs.push_back(HcalElectronicsId(FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
      }

      //	get FEDs registered @cDAQ
      if (auto runInfoRec = es.tryToGet<RunInfoRcd>()) {
        const RunInfo &runInfo = es.getData(runInfoToken_);
        std::vector<int> vfeds = runInfo.m_fed_in;
        for (std::vector<int>::const_iterator it = vfeds.begin(); it != vfeds.end(); ++it) {
          if (*it >= constants::FED_VME_MIN && *it <= FED_VME_MAX)
            _vcdaqEids.push_back(
                HcalElectronicsId(constants::FIBERCH_MIN, constants::FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN)
                    .rawId());
          else if (*it >= constants::FED_uTCA_MIN && *it <= FEDNumbering::MAXHCALuTCAFEDID) {
            std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
            _vcdaqEids.push_back(
                HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
          }
        }
      }
    } else {
      _vCrates = utilities::getCrateList(_emap);
      std::map<int, uint32_t> crateHashMap = utilities::getCrateHashMap(_emap);
      for (auto &it_crate : _vCrates) {
        _vhashCrates.push_back(crateHashMap[it_crate]);
      }

      std::vector<int> vFEDsVME = utilities::getFEDVMEList(_emap);
      std::vector<int> vFEDsuTCA = utilities::getFEDuTCAList(_emap);
      for (std::vector<int>::const_iterator it = vFEDsVME.begin(); it != vFEDsVME.end(); ++it) {
        if (*it <= 0)
          continue;
        _vFEDs.push_back(*it);
        _vhashFEDs.push_back(
            HcalElectronicsId(constants::FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
      }
      for (std::vector<int>::const_iterator it = vFEDsuTCA.begin(); it != vFEDsuTCA.end(); ++it) {
        if (*it <= 0)
          continue;
        std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
        if (cspair.first == 0)
          continue;
        _vFEDs.push_back(*it);
        _vhashFEDs.push_back(
            HcalElectronicsId(cspair.first, cspair.second, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
      }
    }

    // Initialize channel quality masks, but do not load (changed for 10_4_X,
    // moving to LS granularity)
    _xQuality.initialize(hashfunctions::fDChannel);
  }

  void DQClient::beginLuminosityBlock(DQMStore::IBooker &,
                                      DQMStore::IGetter &,
                                      edm::LuminosityBlock const &lb,
                                      edm::EventSetup const &es) {
    //	get the Channel Quality masks
    _xQuality.reset();
    const HcalChannelQuality *cq = &es.getData(hcalChannelQualityToken_);
    std::vector<DetId> detids = cq->getAllChannels();
    for (std::vector<DetId>::const_iterator it = detids.begin(); it != detids.end(); ++it) {
      if (HcalGenericDetId(*it).genericSubdet() == HcalGenericDetId::HcalGenUnknown)
        continue;

      if (HcalGenericDetId(*it).isHcalDetId()) {
        HcalDetId did(*it);
        uint32_t mask = (cq->getValues(did))->getValue();
        if (mask != 0)
          _xQuality.push(did, mask);
      }
    }
  }

  void DQClient::endLuminosityBlock(DQMStore::IBooker &,
                                    DQMStore::IGetter &,
                                    edm::LuminosityBlock const &lb,
                                    edm::EventSetup const &) {
    _currentLS = lb.luminosityBlock();
    _totalLS++;
    if (_maxProcessedLS < _currentLS)
      _maxProcessedLS = _currentLS;
  }

  std::vector<flag::Flag> DQClient::endJob(DQMStore::IBooker &, DQMStore::IGetter &) {
    return std::vector<flag::Flag>();
  }
}  // namespace hcaldqm