#include <RecoLocalCalo/HcalHitSelection/src/HcalHitSelection.cc>

Inheritance diagram for HcalHitSelection:

Public Member Functions
	HcalHitSelection (const edm::ParameterSet &)

	~HcalHitSelection () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

template<typename CollectionType >
void	skim (const edm::Handle< CollectionType > &input, CollectionType &output, int severityThreshold=0) const

Private Attributes
edm::InputTag	hbheTag

edm::InputTag	hfTag

int	hoSeverityLevel

edm::InputTag	hoTag

std::vector< edm::InputTag >	interestingDetIdCollections

edm::ESHandle< HcalChannelQuality >	theHcalChStatus

edm::ESHandle< HcalSeverityLevelComputer >	theHcalSevLvlComputer

const HcalTopology *	theHcalTopology_

std::set< DetId >	toBeKept

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< HFRecHitCollection >	tok_hf_

edm::EDGetTokenT< HORecHitCollection >	tok_ho_

std::vector< edm::EDGetTokenT< DetIdCollection > >	toks_did_

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 49 of file HcalHitSelection.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 118 of file HcalHitSelection.cc.

References edm::ParameterSet::getParameter(), hbheTag, hfTag, hoSeverityLevel, hoTag, mps_fire::i, interestingDetIdCollections, edm::InputTag::label(), tok_hbhe_, tok_hf_, tok_ho_, and toks_did_.

                                                                  {
   hbheTag = iConfig.getParameter<edm::InputTag>("hbheTag");
   hfTag = iConfig.getParameter<edm::InputTag>("hfTag");
   hoTag = iConfig.getParameter<edm::InputTag>("hoTag");
 
   // register for data access
   tok_hbhe_ = consumes<HBHERecHitCollection>(hbheTag);
   tok_hf_ = consumes<HFRecHitCollection>(hfTag);
   tok_ho_ = consumes<HORecHitCollection>(hoTag);
 
   interestingDetIdCollections = iConfig.getParameter<std::vector<edm::InputTag> >("interestingDetIds");
 
   const unsigned nLabels = interestingDetIdCollections.size();
   for (unsigned i = 0; i != nLabels; i++)
     toks_did_.push_back(consumes<DetIdCollection>(interestingDetIdCollections[i]));
 
   hoSeverityLevel = iConfig.getParameter<int>("hoSeverityLevel");
 
   produces<HBHERecHitCollection>(hbheTag.label());
   produces<HFRecHitCollection>(hfTag.label());
   produces<HORecHitCollection>(hoTag.label());
 }

HcalHitSelection::~HcalHitSelection ( )

override

Definition at line 141 of file HcalHitSelection.cc.

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

Member Function Documentation

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

overrideprivate

Definition at line 151 of file HcalHitSelection.cc.

References edm::EDCollection< T >::begin(), DEFINE_FWK_MODULE, edm::EDCollection< T >::end(), edm::EventSetup::get(), edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, hbheTag, photonIsolationHIProducer_cfi::hf, hfTag, photonIsolationHIProducer_cfi::ho, hoSeverityLevel, hoTag, interestingDetIdCollections, edm::HandleBase::isValid(), edm::InputTag::label(), eostools::move(), edm::ESHandle< T >::product(), edm::Event::put(), skim(), OrderedSet::t, theHcalChStatus, theHcalSevLvlComputer, theHcalTopology_, toBeKept, tok_hbhe_, tok_hf_, tok_ho_, and toks_did_.

                                                                             {
   iSetup.get<HcalChannelQualityRcd>().get("withTopo", theHcalChStatus);
   iSetup.get<HcalSeverityLevelComputerRcd>().get(theHcalSevLvlComputer);
   edm::ESHandle<HcalTopology> topo;
   iSetup.get<HcalRecNumberingRecord>().get(topo);
   theHcalTopology_ = topo.product();
 
   edm::Handle<HBHERecHitCollection> hbhe;
   edm::Handle<HFRecHitCollection> hf;
   edm::Handle<HORecHitCollection> ho;
 
   iEvent.getByToken(tok_hbhe_, hbhe);
   iEvent.getByToken(tok_hf_, hf);
   iEvent.getByToken(tok_ho_, ho);
 
   toBeKept.clear();
   edm::Handle<DetIdCollection> detId;
   for (unsigned int t = 0; t < toks_did_.size(); ++t) {
     iEvent.getByToken(toks_did_[t], detId);
     if (!detId.isValid()) {
       edm::LogError("MissingInput") << "the collection of interesting detIds:" << interestingDetIdCollections[t]
                                     << " is not found.";
       continue;
     }
     toBeKept.insert(detId->begin(), detId->end());
   }
 
   auto hbhe_out = std::make_unique<HBHERecHitCollection>();
   skim(hbhe, *hbhe_out);
   iEvent.put(std::move(hbhe_out), hbheTag.label());
 
   auto hf_out = std::make_unique<HFRecHitCollection>();
   skim(hf, *hf_out);
   iEvent.put(std::move(hf_out), hfTag.label());
 
   auto ho_out = std::make_unique<HORecHitCollection>();
   skim(ho, *ho_out, hoSeverityLevel);
   iEvent.put(std::move(ho_out), hoTag.label());
 }

template<class CollectionType >

void HcalHitSelection::skim	(	const edm::Handle< CollectionType > &	input,
		CollectionType &	output,
		int	severityThreshold = `0`
	)		const

private

Definition at line 77 of file HcalHitSelection.cc.

References begin, end, HcalTopology::getMergePositionFlag(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), HcalEndcap, HcalTopology::idFront(), interestingDetIdCollectionProducer_cfi::severityLevel, theHcalChStatus, theHcalSevLvlComputer, theHcalTopology_, and toBeKept.

Referenced by produce().

                                                          {
   output.reserve(input->size());
   typename CollectionType::const_iterator begin = input->begin();
   typename CollectionType::const_iterator end = input->end();
   typename CollectionType::const_iterator hit = begin;
 
   for (; hit != end; ++hit) {
     //    edm::LogError("HcalHitSelection")<<"the hit pointer is"<<&(*hit);
     HcalDetId id = hit->detid();
     if (theHcalTopology_->getMergePositionFlag() && id.subdet() == HcalEndcap) {
       id = theHcalTopology_->idFront(id);
     }
     const uint32_t& recHitFlag = hit->flags();
     //    edm::LogError("HcalHitSelection")<<"the hit id and flag are "<<id.rawId()<<" "<<recHitFlag;
 
     const uint32_t& dbStatusFlag = theHcalChStatus->getValues(id)->getValue();
     int severityLevel = theHcalSevLvlComputer->getSeverityLevel(id, recHitFlag, dbStatusFlag);
     //anything that is not "good" goes in
     if (severityLevel > severityThreshold) {
       output.push_back(*hit);
     } else {
       //chek on the detid list
       if (toBeKept.find(id) != toBeKept.end())
         output.push_back(*hit);
     }
   }
 }