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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
const EcalClusterLazyTools::ESGetTokens	ecalClusterToolsESGetTokens_
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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 46 of file HiSpikeCleaner.cc.

Constructor & Destructor Documentation

HiSpikeCleaner()

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

explicit

Definition at line 67 of file HiSpikeCleaner.cc.

                                                             : 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::~HiSpikeCleaner ( )

override

Definition at line 90 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

produce()

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

overrideprivate

Definition at line 100 of file HiSpikeCleaner.cc.

                                                                          {
  using namespace edm;
 
  // Get raw SuperClusters from the event
  Handle<reco::SuperClusterCollection> pRawSuperClusters;
  try {
    iEvent.getByToken(sCInputProducerToken_, pRawSuperClusters);
  } catch (cms::Exception& ex) {
    edm::LogError("EgammaSCCorrectionMakerError") << "Error! can't get the rawSuperClusters ";
  }
 
  // Get the RecHits from the event
  Handle<EcalRecHitCollection> pRecHitsB;
  try {
    iEvent.getByToken(rHInputProducerBToken_, pRecHitsB);
  } catch (cms::Exception& ex) {
    edm::LogError("EgammaSCCorrectionMakerError") << "Error! can't get the RecHits ";
  }
 
  // Get the RecHits from the event
  Handle<EcalRecHitCollection> pRecHitsE;
  try {
    iEvent.getByToken(rHInputProducerEToken_, pRecHitsE);
  } catch (cms::Exception& ex) {
    edm::LogError("EgammaSCCorrectionMakerError") << "Error! can't get the RecHits ";
  }
 
  // get the channel status from the DB
  //   edm::ESHandle<EcalChannelStatus> chStatus;
  //   iSetup.get<EcalChannelStatusRcd>().get(chStatus);
 
  edm::ESHandle<EcalSeverityLevelAlgo> ecalSevLvlAlgoHndl;
  iSetup.get<EcalSeverityLevelAlgoRcd>().get(ecalSevLvlAlgoHndl);
 
  // 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()) {
      ecalSevLvlAlgoHndl->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_);
}

References ecalClusterToolsESGetTokens_, etCut_, EcalClusterLazyToolsBase::ESGetTokens::get(), edm::EventSetup::get(), get, iEvent, eostools::move(), outputCollection_, edm::Handle< T >::product(), HI_PhotonSkim_cff::rechits, rHInputProducerBToken_, rHInputProducerEToken_, sCInputProducerToken_, fileCollector::seed, EcalSeverityLevelAlgo::severityLevel(), EcalTools::swissCross(), swissCutThr_, and TimingCut_.

Member Data Documentation

ecalClusterToolsESGetTokens_

const EcalClusterLazyTools::ESGetTokens HiSpikeCleaner::ecalClusterToolsESGetTokens_

private

Definition at line 59 of file HiSpikeCleaner.cc.

Referenced by produce().

etCut_

double HiSpikeCleaner::etCut_

private

Definition at line 64 of file HiSpikeCleaner.cc.

Referenced by produce().

outputCollection_

std::string HiSpikeCleaner::outputCollection_

private

Definition at line 61 of file HiSpikeCleaner.cc.

Referenced by produce().

rHInputProducerBToken_

edm::EDGetTokenT<EcalRecHitCollection> HiSpikeCleaner::rHInputProducerBToken_

private

Definition at line 57 of file HiSpikeCleaner.cc.

Referenced by produce().

rHInputProducerEToken_

edm::EDGetTokenT<EcalRecHitCollection> HiSpikeCleaner::rHInputProducerEToken_

private

Definition at line 58 of file HiSpikeCleaner.cc.

Referenced by produce().

sCInputProducerToken_

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

private

Definition at line 56 of file HiSpikeCleaner.cc.

Referenced by produce().

swissCutThr_

double HiSpikeCleaner::swissCutThr_

private

Definition at line 63 of file HiSpikeCleaner.cc.

Referenced by produce().

TimingCut_

double HiSpikeCleaner::TimingCut_

private

Definition at line 62 of file HiSpikeCleaner.cc.

Referenced by produce().