IntegrityTask.cc

Go to the documentation of this file.

#include "DQM/EcalMonitorTasks/interface/IntegrityTask.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"
 
namespace ecaldqm {
  IntegrityTask::IntegrityTask() : DQWorkerTask() {}
 
  void IntegrityTask::beginEvent(edm::Event const& _evt,
                                 edm::EventSetup const& _es,
                                 bool const& ByLumiResetSwitch,
                                 bool&) {
    if (ByLumiResetSwitch) {
      MEs_.at("MapByLumi").reset(GetElectronicsMap());
      MEs_.at("ByLumi").reset(GetElectronicsMap());
    }
  }
 
  template <typename IDCollection>
  void IntegrityTask::runOnDetIdCollection(IDCollection const& _ids, Collections _collection) {
    // Collection is empty if there are no errors
    if (_ids.empty())
      return;
 
    MESet* set(nullptr);
    switch (_collection) {
      case kEBGainErrors:
      case kEEGainErrors:
        set = &MEs_.at("Gain");
        break;
      case kEBChIdErrors:
      case kEEChIdErrors:
        set = &MEs_.at("ChId");
        break;
      case kEBGainSwitchErrors:
      case kEEGainSwitchErrors:
        set = &MEs_.at("GainSwitch");
        break;
      default:
        return;
    }
 
    MESet& meMapByLumi(MEs_.at("MapByLumi"));
    MESet& meByLumi(MEs_.at("ByLumi"));
    MESet& meTotal(MEs_.at("Total"));
    MESet& meTrendNErrors(MEs_.at("TrendNErrors"));
 
    std::for_each(_ids.begin(), _ids.end(), [&](typename IDCollection::value_type const& id) {
      set->fill(getEcalDQMSetupObjects(), id);
      int dccid(dccId(id, GetElectronicsMap()));
      meByLumi.fill(getEcalDQMSetupObjects(), dccid);
      meTotal.fill(getEcalDQMSetupObjects(), dccid);
      // Fill Integrity Errors Map with channel errors for this lumi
      meMapByLumi.fill(getEcalDQMSetupObjects(), id);
 
      meTrendNErrors.fill(getEcalDQMSetupObjects(), double(timestamp_.iLumi), 1.);
    });
  }
 
  void IntegrityTask::runOnElectronicsIdCollection(EcalElectronicsIdCollection const& _ids, Collections _collection) {
    // Collection is empty if there are no errors
    if (_ids.empty())
      return;
 
    MESet* set(nullptr);
    switch (_collection) {
      case kTowerIdErrors:
        set = &MEs_.at("TowerId");
        break;
      case kBlockSizeErrors:
        set = &MEs_.at("BlockSize");
        break;
      default:
        return;
    }
 
    MESet& meMapByLumi(MEs_.at("MapByLumi"));
    MESet& meByLumi(MEs_.at("ByLumi"));
    MESet& meTotal(MEs_.at("Total"));
    MESet& meTrendNErrors(MEs_.at("TrendNErrors"));
 
    std::for_each(_ids.begin(), _ids.end(), [&](EcalElectronicsIdCollection::value_type const& id) {
      set->fill(getEcalDQMSetupObjects(), id);
      int dccid(id.dccId());
      double nCrystals(0.);
      std::vector<DetId> chIds(GetElectronicsMap()->dccTowerConstituents(dccid, id.towerId()));
      if (dccid <= kEEmHigh + 1 || dccid >= kEEpLow + 1)
        nCrystals = chIds.size();
      else
        nCrystals = 25.;
      meByLumi.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
      meTotal.fill(getEcalDQMSetupObjects(), dccid, nCrystals);
      // Fill Integrity Errors Map with tower errors for this lumi
      // Since binned by crystal for compatibility with channel errors,
      // fill with constituent channels of tower
      for (std::vector<DetId>::iterator chItr(chIds.begin()); chItr != chIds.end(); ++chItr)
        meMapByLumi.fill(getEcalDQMSetupObjects(), *chItr);
 
      meTrendNErrors.fill(getEcalDQMSetupObjects(), double(timestamp_.iLumi), nCrystals);
    });
  }
 
  DEFINE_ECALDQM_WORKER(IntegrityTask);
}  // namespace ecaldqm