#include <HLTHPDFilter.h>

Inheritance diagram for HLTHPDFilter:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override

	HLTHPDFilter (const edm::ParameterSet &)

	~HLTHPDFilter () override

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< HBHERecHitCollection >	m_theRecHitCollectionToken

double	mEnergyThreshold

double	mHPDSpikeEnergyThreshold

double	mHPDSpikeIsolationEnergyThreshold

edm::InputTag	mInputTag

double	mRBXSpikeEnergyThreshold

double	mRBXSpikeUnbalanceThreshold

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
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

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Author: Fedor Ratnikov (UMd)

Definition at line 21 of file HLTHPDFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 76 of file HLTHPDFilter.cc.

References m_theRecHitCollectionToken, mInputTag, and ~HLTHPDFilter().

                                                          :
      mInputTag (iConfig.getParameter <edm::InputTag> ("inputTag")),
      mEnergyThreshold (iConfig.getParameter <double> ("energy")),
      mHPDSpikeEnergyThreshold (iConfig.getParameter <double> ("hpdSpikeEnergy")),
      mHPDSpikeIsolationEnergyThreshold (iConfig.getParameter <double> ("hpdSpikeIsolationEnergy")),
      mRBXSpikeEnergyThreshold (iConfig.getParameter <double> ("rbxSpikeEnergy")),
      mRBXSpikeUnbalanceThreshold (iConfig.getParameter <double> ("rbxSpikeUnbalance"))
 {
   m_theRecHitCollectionToken = consumes<HBHERecHitCollection>(mInputTag);
 }

HLTHPDFilter::~HLTHPDFilter ( )

overridedefault

Referenced by HLTHPDFilter().

Member Function Documentation

void HLTHPDFilter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 90 of file HLTHPDFilter.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                          {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("inputTag",edm::InputTag("hltHbhereco"));
   desc.add<double>("energy",-99.0);
   desc.add<double>("hpdSpikeEnergy",10.0);
   desc.add<double>("hpdSpikeIsolationEnergy",1.0);
   desc.add<double>("rbxSpikeEnergy",50.0);
   desc.add<double>("rbxSpikeUnbalance",0.2);
   descriptions.add("hltHPDFilter",desc);
 }

bool HLTHPDFilter::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 101 of file HLTHPDFilter.cc.

References edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, mps_fire::i, createfilelist::int, m_theRecHitCollectionToken, particleFlowClusterECALTimeSelected_cfi::maxEnergy, mEnergyThreshold, mHPDSpikeEnergyThreshold, mHPDSpikeIsolationEnergyThreshold, particleFlowClusterECALTimeSelected_cfi::minEnergy, mRBXSpikeEnergyThreshold, mRBXSpikeUnbalanceThreshold, and edm::SortedCollection< T, SORT >::size().

 {
   if (mHPDSpikeEnergyThreshold <= 0 && mRBXSpikeEnergyThreshold <= 0) return true; // nothing to filter
   // get hits
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(m_theRecHitCollectionToken,hbhe);
   
   // collect energies
   float hpdEnergy[4][73];
   for (auto & i : hpdEnergy) for (size_t j = 0; j < 73; ++j) i[j] = 0;
   
   // select hist above threshold
   for (unsigned i = 0; i < hbhe->size(); ++i) {
     if ((*hbhe)[i].energy() > mEnergyThreshold) {
       std::pair<Partition,int> hpd = hpdId ((*hbhe)[i].id());
       hpdEnergy[int (hpd.first)][hpd.second] += (*hbhe)[i].energy ();
     }
   }
   
   // not single HPD spike
   if (mHPDSpikeEnergyThreshold > 0) {
     for (auto & partition : hpdEnergy) {
       for (size_t i = 1; i < 73; ++i) {
     if (partition[i] > mHPDSpikeEnergyThreshold) {
       int hpdPlus = i + 1;
       if (hpdPlus == 73) hpdPlus = 1;
       int hpdMinus = i - 1;
       if (hpdMinus == 0) hpdMinus = 72;
       double maxNeighborEnergy = fmax (partition[hpdPlus], partition[hpdMinus]);
       if (maxNeighborEnergy < mHPDSpikeIsolationEnergyThreshold)  return false; // HPD spike found
     }
       }
     }  
   }
 
   // not RBX flash
   if (mRBXSpikeEnergyThreshold > 0) {
     for (auto & partition : hpdEnergy) {
       for (size_t rbx = 1; rbx < 19; ++rbx) {
     int ifirst = (rbx-1)*4-1;
     int iend = (rbx-1)*4+3;
     double minEnergy = 0;
     double maxEnergy = -1;
     double totalEnergy = 0;
     for (int irm = ifirst; irm < iend; ++irm) {
       int hpd = irm;
       if (hpd <= 0) hpd = 72 + hpd;
       totalEnergy += partition[hpd];
       if (minEnergy > maxEnergy) {
         minEnergy = maxEnergy = partition[hpd];
       }
       else {
         if (partition[hpd] < minEnergy) minEnergy = partition[hpd];
         if (partition[hpd] > maxEnergy) maxEnergy = partition[hpd];
       }
     }
     if (totalEnergy > mRBXSpikeEnergyThreshold) {
       if (minEnergy / maxEnergy > mRBXSpikeUnbalanceThreshold) return false; // likely HPF flash
     }
       }
     }
   }
   return true;
 }

Member Data Documentation

edm::EDGetTokenT<HBHERecHitCollection> HLTHPDFilter::m_theRecHitCollectionToken

private

Definition at line 30 of file HLTHPDFilter.h.

Referenced by filter(), and HLTHPDFilter().

double HLTHPDFilter::mEnergyThreshold

private

Definition at line 32 of file HLTHPDFilter.h.

Referenced by filter().

double HLTHPDFilter::mHPDSpikeEnergyThreshold

private

Definition at line 33 of file HLTHPDFilter.h.

Referenced by filter().

double HLTHPDFilter::mHPDSpikeIsolationEnergyThreshold

private

Definition at line 34 of file HLTHPDFilter.h.

Referenced by filter().

edm::InputTag HLTHPDFilter::mInputTag

private

Definition at line 31 of file HLTHPDFilter.h.

Referenced by HLTHPDFilter().

double HLTHPDFilter::mRBXSpikeEnergyThreshold

private

Definition at line 35 of file HLTHPDFilter.h.

Referenced by filter().

double HLTHPDFilter::mRBXSpikeUnbalanceThreshold

private

Definition at line 36 of file HLTHPDFilter.h.

Referenced by filter().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation