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::beginLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) {
    // Reset by LS plots at beginning of every LS
    MEs_.at("MapByLumi").reset();
    MEs_.at("ByLumi").reset();
  }
 
  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(id);
      int dccid(dccId(id));
      meByLumi.fill(dccid);
      meTotal.fill(dccid);
      // Fill Integrity Errors Map with channel errors for this lumi
      meMapByLumi.fill(id);
 
      meTrendNErrors.fill(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(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(dccid, nCrystals);
      meTotal.fill(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(*chItr);
 
      meTrendNErrors.fill(double(timestamp_.iLumi), nCrystals);
    });
  }
 
  DEFINE_ECALDQM_WORKER(IntegrityTask);
}  // namespace ecaldqm