#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	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

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

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	htopoToken_

std::vector< edm::InputTag >	interestingDetIdCollections

edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd >	qualToken_

edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	sevToken_

const HcalChannelQuality *	theHcalChStatus

const 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::HcalHitSelection ( const edm::ParameterSet & iConfig )

explicit

Definition at line 121 of file HcalHitSelection.cc.

     : hbheTag(iConfig.getParameter<edm::InputTag>("hbheTag")),
       hoTag(iConfig.getParameter<edm::InputTag>("hoTag")),
       hfTag(iConfig.getParameter<edm::InputTag>("hfTag")),
       theHcalTopology_(nullptr),
       theHcalChStatus(nullptr),
       theHcalSevLvlComputer(nullptr) {
   // 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());
  
   // ES tokens
   htopoToken_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   qualToken_ = esConsumes<HcalChannelQuality, HcalChannelQualityRcd>(edm::ESInputTag("", "withTopo"));
   sevToken_ = esConsumes<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd>();
 }

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

◆ ~HcalHitSelection()

HcalHitSelection::~HcalHitSelection ( )

override

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

◆ produce()

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

overrideprivate

Definition at line 161 of file HcalHitSelection.cc.

                                                                             {
   theHcalChStatus = &iSetup.getData(qualToken_);
   theHcalSevLvlComputer = &iSetup.getData(sevToken_);
   theHcalTopology_ = &iSetup.getData(htopoToken_);
  
   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());
 }

References edm::EDCollection< T >::begin(), edm::EDCollection< T >::end(), edm::EventSetup::getData(), photonIsolationHIProducer_cfi::hbhe, hbheTag, photonIsolationHIProducer_cfi::hf, hfTag, photonIsolationHIProducer_cfi::ho, hoSeverityLevel, hoTag, htopoToken_, iEvent, interestingDetIdCollections, edm::HandleBase::isValid(), edm::InputTag::label(), eostools::move(), qualToken_, sevToken_, skim(), submitPVValidationJobs::t, theHcalChStatus, theHcalSevLvlComputer, theHcalTopology_, toBeKept, tok_hbhe_, tok_hf_, tok_ho_, and toks_did_.

◆ skim()

template<class CollectionType >

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

private

Definition at line 80 of file HcalHitSelection.cc.

                                                          {
   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);
     }
   }
 }