HcalHitSelection Class Reference

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

explicit

Definition at line 118 of file HcalHitSelection.cc.

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

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

~HcalHitSelection()

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

produce()

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

overrideprivate

Definition at line 151 of file HcalHitSelection.cc.

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

References edm::EDCollection< T >::begin(), edm::EDCollection< T >::end(), edm::EventSetup::get(), get, photonIsolationHIProducer_cfi::hbhe, hbheTag, photonIsolationHIProducer_cfi::hf, hfTag, photonIsolationHIProducer_cfi::ho, hoSeverityLevel, hoTag, iEvent, interestingDetIdCollections, edm::HandleBase::isValid(), edm::InputTag::label(), eostools::move(), edm::ESHandle< T >::product(), skim(), OrderedSet::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 77 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);
    }
  }
}

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

Referenced by produce().

Member Data Documentation

hbheTag

edm::InputTag HcalHitSelection::hbheTag

private

Definition at line 57 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

hfTag

edm::InputTag HcalHitSelection::hfTag

private

Definition at line 57 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

hoSeverityLevel

int HcalHitSelection::hoSeverityLevel

private

Definition at line 62 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

hoTag

edm::InputTag HcalHitSelection::hoTag

private

Definition at line 57 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

interestingDetIdCollections

std::vector<edm::InputTag> HcalHitSelection::interestingDetIdCollections

private

Definition at line 63 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

theHcalChStatus

edm::ESHandle<HcalChannelQuality> HcalHitSelection::theHcalChStatus

private

Definition at line 67 of file HcalHitSelection.cc.

Referenced by produce(), and skim().

theHcalSevLvlComputer

edm::ESHandle<HcalSeverityLevelComputer> HcalHitSelection::theHcalSevLvlComputer

private

Definition at line 68 of file HcalHitSelection.cc.

Referenced by produce(), and skim().

theHcalTopology_

const HcalTopology* HcalHitSelection::theHcalTopology_

private

Definition at line 64 of file HcalHitSelection.cc.

Referenced by produce(), and skim().

toBeKept

std::set<DetId> HcalHitSelection::toBeKept

private

Definition at line 69 of file HcalHitSelection.cc.

Referenced by produce(), and skim().

tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> HcalHitSelection::tok_hbhe_

private

Definition at line 58 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

tok_hf_

edm::EDGetTokenT<HFRecHitCollection> HcalHitSelection::tok_hf_

private

Definition at line 60 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

tok_ho_

edm::EDGetTokenT<HORecHitCollection> HcalHitSelection::tok_ho_

private

Definition at line 59 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().

toks_did_

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

private

Definition at line 61 of file HcalHitSelection.cc.

Referenced by HcalHitSelection(), and produce().