#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

	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

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

Constructor & Destructor Documentation

◆ HcalHitSelection()

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

explicit

Definition at line 121 of file HcalHitSelection.cc.

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

     : 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>();
 }

◆ ~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.

References hcalRecHitTable_cff::detId, edm::EventSetup::getData(), photonIsolationHIProducer_cfi::hbhe, hbheTag, Hcal_Conditions_forGlobalTag_cff::hf, hfTag, Hcal_Conditions_forGlobalTag_cff::ho, hoSeverityLevel, hoTag, htopoToken_, iEvent, interestingDetIdCollections, edm::InputTag::label(), eostools::move(), qualToken_, sevToken_, skim(), submitPVValidationJobs::t, theHcalChStatus, theHcalSevLvlComputer, theHcalTopology_, toBeKept, tok_hbhe_, tok_hf_, tok_ho_, and toks_did_.

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

◆ 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.

References SplitLinear::begin, mps_fire::end, HcalTopology::getMergePositionFlag(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), HcalEndcap, HcalTopology::idFront(), input, 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);
     }
   }
 }