HLTHPDFilter Class Reference

#include <HLTHPDFilter.h>

Inheritance diagram for HLTHPDFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	HLTHPDFilter (const edm::ParameterSet &)
	~HLTHPDFilter ()
Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default
Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDFilterBase ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	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 ()

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
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 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, and mInputTag.

                                                         :
     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

(

)

Definition at line 87 of file HLTHPDFilter.cc.

{}

Member Function Documentation

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

static

Definition at line 90 of file HLTHPDFilter.cc.

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

                                                                         {
  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  )

virtual

Implements edm::stream::EDFilterBase.

Definition at line 101 of file HLTHPDFilter.cc.

References edm::Event::getByToken(), i, j, m_theRecHitCollectionToken, mEnergyThreshold, mHPDSpikeEnergyThreshold, mHPDSpikeIsolationEnergyThreshold, mRBXSpikeEnergyThreshold, and mRBXSpikeUnbalanceThreshold.

{
  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 (size_t i = 0; i < 4; ++i) for (size_t j = 0; j < 73; ++j) hpdEnergy[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 (size_t partition = 0; partition < 4; ++partition) {
      for (size_t i = 1; i < 73; ++i) {
    if (hpdEnergy [partition][i] > mHPDSpikeEnergyThreshold) {
      int hpdPlus = i + 1;
      if (hpdPlus == 73) hpdPlus = 1;
      int hpdMinus = i - 1;
      if (hpdMinus == 0) hpdMinus = 72;
      double maxNeighborEnergy = fmax (hpdEnergy[partition][hpdPlus], hpdEnergy[partition][hpdMinus]);
      if (maxNeighborEnergy < mHPDSpikeIsolationEnergyThreshold)  return false; // HPD spike found
    }
      }
    }  
  }

  // not RBX flash
  if (mRBXSpikeEnergyThreshold > 0) {
    for (size_t partition = 0; partition < 4; ++partition) {
      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 += hpdEnergy[partition][hpd];
      if (minEnergy > maxEnergy) {
        minEnergy = maxEnergy = hpdEnergy[partition][hpd];
      }
      else {
        if (hpdEnergy[partition][hpd] < minEnergy) minEnergy = hpdEnergy[partition][hpd];
        if (hpdEnergy[partition][hpd] > maxEnergy) maxEnergy = hpdEnergy[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().