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

Inheritance diagram for HcalHitSelection:

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

	~HcalHitSelection ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

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

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::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 46 of file HcalHitSelection.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 108 of file HcalHitSelection.cc.

References edm::ParameterSet::getParameter(), hbheTag, hfTag, hoSeverityLevel, hoTag, 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 ( )

Definition at line 134 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
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 149 of file HcalHitSelection.cc.

References edm::EventSetup::get(), edm::Event::getByToken(), hbheTag, hfTag, hoSeverityLevel, hoTag, interestingDetIdCollections, edm::HandleBase::isValid(), edm::InputTag::label(), edm::Event::put(), skim(), edmStreamStallGrapher::t, theHcalChStatus, theHcalSevLvlComputer, toBeKept, tok_hbhe_, tok_hf_, tok_ho_, and toks_did_.

 {
   iSetup.get<HcalChannelQualityRcd>().get("withTopo", theHcalChStatus);
   iSetup.get<HcalSeverityLevelComputerRcd>().get(theHcalSevLvlComputer);
 
   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());
     }
 
   std::auto_ptr<HBHERecHitCollection> hbhe_out(new HBHERecHitCollection());
   skim(hbhe,hbhe_out);
   iEvent.put(hbhe_out,hbheTag.label());
 
   std::auto_ptr<HFRecHitCollection> hf_out(new HFRecHitCollection());
   skim(hf,hf_out);
   iEvent.put(hf_out,hfTag.label());
 
   std::auto_ptr<HORecHitCollection> ho_out(new HORecHitCollection());
   skim(ho,ho_out,hoSeverityLevel);
   iEvent.put(ho_out,hoTag.label());
   
 }

template<class CollectionType >

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

private

Definition at line 71 of file HcalHitSelection.cc.

References begin, end, theHcalChStatus, theHcalSevLvlComputer, 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);
     const DetId & id = hit->detid();
     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);
     }
   }
 }