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_;

  // threshold above which a channel will be considered "dead"
  int statusThreshold_;
};

EcalNextToDeadChannelESProducer::EcalNextToDeadChannelESProducer(const edm::ParameterSet& iConfig){
  //the following line is needed to tell the framework what
  // data is being produced
  setWhatProduced(this);

  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

  edm::ESHandle <EcalChannelStatus> h;
  chs.get (h);
  
  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);

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "cd .. ; scram b"
// End: