EcalNextToDeadChannelESProducer.cc

Go to the documentation of this file.

#include <memory>
#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/ESProductHost.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/ReusableObjectHolder.h"
#include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalNextToDeadChannel.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalNextToDeadChannelRcd.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalTools.h"

class EcalNextToDeadChannelESProducer : public edm::ESProducer {
public:
  EcalNextToDeadChannelESProducer(const edm::ParameterSet& iConfig);

  typedef std::shared_ptr<EcalNextToDeadChannel> ReturnType;

  ReturnType produce(const EcalNextToDeadChannelRcd& iRecord);

private:
  void setupNextToDeadChannels(const EcalChannelStatusRcd&, EcalNextToDeadChannel*);

  using HostType = edm::ESProductHost<EcalNextToDeadChannel, EcalChannelStatusRcd>;

  edm::ReusableObjectHolder<HostType> holder_;

  edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> const channelToken_;
  // threshold above which a channel will be considered "dead"
  int statusThreshold_;
};

EcalNextToDeadChannelESProducer::EcalNextToDeadChannelESProducer(const edm::ParameterSet& iConfig)
    : channelToken_(setWhatProduced(this).consumesFrom<EcalChannelStatus, EcalChannelStatusRcd>()) {
  statusThreshold_ = iConfig.getParameter<int>("channelStatusThresholdForDead");
}

EcalNextToDeadChannelESProducer::ReturnType EcalNextToDeadChannelESProducer::produce(
    const EcalNextToDeadChannelRcd& iRecord) {
  auto host = holder_.makeOrGet([]() { return new HostType; });

  host->ifRecordChanges<EcalChannelStatusRcd>(
      iRecord, [this, h = host.get()](auto const& rec) { setupNextToDeadChannels(rec, h); });

  return host;
}

void EcalNextToDeadChannelESProducer::setupNextToDeadChannels(const EcalChannelStatusRcd& chs,
                                                              EcalNextToDeadChannel* rcd) {
  rcd->clear();

  // Find channels next to dead ones and fill corresponding record

  EcalChannelStatus const& h = chs.get(channelToken_);

  for (int ieta = -EBDetId::MAX_IETA; ieta <= EBDetId::MAX_IETA; ++ieta) {
    if (ieta == 0)
      continue;
    for (int iphi = EBDetId::MIN_IPHI; iphi <= EBDetId::MAX_IPHI; ++iphi) {
      if (EBDetId::validDetId(ieta, iphi)) {
        EBDetId detid(ieta, iphi);

        if (EcalTools::isNextToDeadFromNeighbours(detid, h, statusThreshold_)) {
          rcd->setValue(detid, 1);
        };
      }
    }  // for phi
  }    // for eta

  // endcap

  for (int iX = EEDetId::IX_MIN; iX <= EEDetId::IX_MAX; ++iX) {
    for (int iY = EEDetId::IY_MIN; iY <= EEDetId::IY_MAX; ++iY) {
      if (EEDetId::validDetId(iX, iY, 1)) {
        EEDetId detid(iX, iY, 1);

        if (EcalTools::isNextToDeadFromNeighbours(detid, h, statusThreshold_)) {
          rcd->setValue(detid, 1);
        };
      }

      if (EEDetId::validDetId(iX, iY, -1)) {
        EEDetId detid(iX, iY, -1);

        if (EcalTools::isNextToDeadFromNeighbours(detid, h, statusThreshold_)) {
          rcd->setValue(detid, 1);
        };
      }
    }  // for iy
  }    // for ix
}

//define this as a plug-in
DEFINE_FWK_EVENTSETUP_MODULE(EcalNextToDeadChannelESProducer);