HiSpikeCleaner Class Reference

#include <RecoHI/HiSpikeCleaner/src/HiSpikeCleaner.cc>

Inheritance diagram for HiSpikeCleaner:

Public Member Functions
	HiSpikeCleaner (const edm::ParameterSet &)
	~HiSpikeCleaner () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_
edm::ESGetToken < EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	ecalSevLvlAlgoToken_
double	etCut_
std::string	outputCollection_
edm::EDGetTokenT < EcalRecHitCollection >	rHInputProducerBToken_
edm::EDGetTokenT < EcalRecHitCollection >	rHInputProducerEToken_
edm::EDGetTokenT < reco::SuperClusterCollection >	sCInputProducerToken_
double	swissCutThr_
double	TimingCut_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 45 of file HiSpikeCleaner.cc.

Constructor & Destructor Documentation

HiSpikeCleaner::HiSpikeCleaner ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 67 of file HiSpikeCleaner.cc.

    : ecalSevLvlAlgoToken_(esConsumes()), ecalClusterToolsESGetTokens_{consumesCollector()} {
  //register your products
  /* Examples
   produces<ExampleData2>();

   //if do put with a label
   produces<ExampleData2>("label");
*/
  //now do what ever other initialization is needed

  rHInputProducerBToken_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitProducerBarrel"));
  rHInputProducerEToken_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHitProducerEndcap"));

  sCInputProducerToken_ =
      consumes<reco::SuperClusterCollection>(iConfig.getParameter<edm::InputTag>("originalSuperClusterProducer"));
  TimingCut_ = iConfig.getUntrackedParameter<double>("TimingCut", 4.0);
  swissCutThr_ = iConfig.getUntrackedParameter<double>("swissCutThr", 0.95);
  etCut_ = iConfig.getParameter<double>("etCut");

  outputCollection_ = iConfig.getParameter<std::string>("outputColl");
  produces<reco::SuperClusterCollection>(outputCollection_);
}

HiSpikeCleaner::~HiSpikeCleaner ( )

override

Definition at line 91 of file HiSpikeCleaner.cc.

                                {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

Member Function Documentation

void HiSpikeCleaner::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 101 of file HiSpikeCleaner.cc.

References ecalClusterToolsESGetTokens_, ecalSevLvlAlgoToken_, etCut_, EcalClusterLazyToolsBase::ESGetTokens::get(), edm::Event::getByToken(), edm::EventSetup::getData(), eostools::move(), outputCollection_, edm::Event::put(), HI_PhotonSkim_cff::rechits, rHInputProducerBToken_, rHInputProducerEToken_, sCInputProducerToken_, fileCollector::seed, EcalTools::swissCross(), swissCutThr_, and TimingCut_.

                                                                          {
  using namespace edm;

  // Get raw SuperClusters from the event
  Handle<reco::SuperClusterCollection> pRawSuperClusters;
  iEvent.getByToken(sCInputProducerToken_, pRawSuperClusters);

  // Get the RecHits from the event
  Handle<EcalRecHitCollection> pRecHitsB;
  iEvent.getByToken(rHInputProducerBToken_, pRecHitsB);

  // Get the RecHits from the event
  Handle<EcalRecHitCollection> pRecHitsE;
  iEvent.getByToken(rHInputProducerEToken_, pRecHitsE);

  // get the channel status from the DB
  //   edm::ESHandle<EcalChannelStatus> chStatus;
  //   iSetup.get<EcalChannelStatusRcd>().get(chStatus);

  auto const& ecalSevLvlAlgo = iSetup.getData(ecalSevLvlAlgoToken_);

  // Create a pointer to the RecHits and raw SuperClusters
  const reco::SuperClusterCollection* rawClusters = pRawSuperClusters.product();

  EcalClusterLazyTools lazyTool(
      iEvent, ecalClusterToolsESGetTokens_.get(iSetup), rHInputProducerBToken_, rHInputProducerEToken_);

  // Define a collection of corrected SuperClusters to put back into the event
  auto corrClusters = std::make_unique<reco::SuperClusterCollection>();

  //  Loop over raw clusters and make corrected ones
  reco::SuperClusterCollection::const_iterator aClus;
  for (aClus = rawClusters->begin(); aClus != rawClusters->end(); aClus++) {
    double theEt = aClus->energy() / cosh(aClus->eta());
    //   std::cout << " et of SC = " << theEt << std::endl;

    if (theEt < etCut_)
      continue;  // cut off low pT superclusters

    bool flagS = true;
    float swissCrx(0);

    const reco::CaloClusterPtr seed = aClus->seed();
    DetId id = lazyTool.getMaximum(*seed).first;
    const EcalRecHitCollection& rechits = *pRecHitsB;
    EcalRecHitCollection::const_iterator it = rechits.find(id);

    if (it != rechits.end()) {
      ecalSevLvlAlgo.severityLevel(id, rechits);
      swissCrx = EcalTools::swissCross(id, rechits, 0., true);
      //        std::cout << "swissCross = " << swissCrx <<std::endl;
      // std::cout << " timing = " << it->time() << std::endl;
    }

    if (fabs(it->time()) > TimingCut_) {
      flagS = false;
      //        std::cout << " timing = " << it->time() << std::endl;
      //   std::cout << " timing is bad........" << std::endl;
    }
    if (swissCrx > (float)swissCutThr_) {
      flagS = false;  // swissCross cut
                      //        std::cout << "swissCross = " << swissCrx <<std::endl;
                      //   std::cout << " removed by swiss cross cut" << std::endl;
    }
    // - kGood        --> good channel
    // - kProblematic --> problematic (e.g. noisy)
    // - kRecovered   --> recovered (e.g. an originally dead or saturated)
    // - kTime        --> the channel is out of time (e.g. spike)
    // - kWeird       --> weird (e.g. spike)
    // - kBad         --> bad, not suitable to be used in the reconstruction
    //   enum EcalSeverityLevel { kGood=0, kProblematic, kRecovered, kTime, kWeird, kBad };

    reco::SuperCluster newClus;
    if (flagS == true)
      newClus = *aClus;
    else
      continue;
    corrClusters->push_back(newClus);
  }

  // Put collection of corrected SuperClusters into the event
  iEvent.put(std::move(corrClusters), outputCollection_);
}

Member Data Documentation

const EcalClusterLazyTools::ESGetTokens HiSpikeCleaner::ecalClusterToolsESGetTokens_

private

Definition at line 59 of file HiSpikeCleaner.cc.

Referenced by produce().

edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> HiSpikeCleaner::ecalSevLvlAlgoToken_

private

Definition at line 58 of file HiSpikeCleaner.cc.

Referenced by produce().

double HiSpikeCleaner::etCut_

private

Definition at line 64 of file HiSpikeCleaner.cc.

Referenced by produce().

std::string HiSpikeCleaner::outputCollection_

private

Definition at line 61 of file HiSpikeCleaner.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> HiSpikeCleaner::rHInputProducerBToken_

private

Definition at line 56 of file HiSpikeCleaner.cc.

Referenced by produce().

edm::EDGetTokenT<EcalRecHitCollection> HiSpikeCleaner::rHInputProducerEToken_

private

Definition at line 57 of file HiSpikeCleaner.cc.

Referenced by produce().

edm::EDGetTokenT<reco::SuperClusterCollection> HiSpikeCleaner::sCInputProducerToken_

private

Definition at line 55 of file HiSpikeCleaner.cc.

Referenced by produce().

double HiSpikeCleaner::swissCutThr_

private

Definition at line 63 of file HiSpikeCleaner.cc.

Referenced by produce().

double HiSpikeCleaner::TimingCut_

private

Definition at line 62 of file HiSpikeCleaner.cc.

Referenced by produce().