#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 50 of file HcalHitSelection.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 116 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 142 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 157 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(), protons_cff::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 76 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);
     }
   }
 }