#include <HLTHcalMETNoiseFilter.h>

Inheritance diagram for HLTHcalMETNoiseFilter:

Classes
struct	noisedatacomp

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

	HLTHcalMETNoiseFilter (const edm::ParameterSet &)

	~HLTHcalMETNoiseFilter () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () 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::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Private Types
typedef std::set< CommonHcalNoiseRBXData, noisedatacomp >	noisedataset_t

Private Attributes
edm::InputTag	HcalNoiseRBXCollectionTag_

edm::EDGetTokenT< reco::HcalNoiseRBXCollection >	m_theHcalNoiseToken

double	maxHighEHitTime_

int	maxNumRBXs_

double	maxRatio_

double	maxRBXEMF_

double	minHighEHitTime_

double	minHighHitE_

int	minHPDHits_

int	minHPDNoOtherHits_

double	minLowHitE_

double	minR45HitE_

double	minRatio_

double	minRBXEnergy_

int	minRBXHits_

double	minRecHitE_

int	minZeros_

bool	needEMFCoincidence_

int	numRBXsToConsider_

int	severity_

double	TS4TS5EnergyThreshold_

std::vector< std::pair< double, double > >	TS4TS5LowerCut_

std::vector< std::pair< double, double > >	TS4TS5UpperCut_

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

Description: HLT filter module for rejecting MET events due to noise in the HCAL

Implementation: <Notes on="" implementation>="">

Definition at line 23 of file HLTHcalMETNoiseFilter.h.

Member Typedef Documentation

typedef std::set<CommonHcalNoiseRBXData, noisedatacomp> HLTHcalMETNoiseFilter::noisedataset_t

private

Definition at line 63 of file HLTHcalMETNoiseFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 34 of file HLTHcalMETNoiseFilter.cc.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), HcalNoiseRBXCollectionTag_, mps_fire::i, createfilelist::int, m_theHcalNoiseToken, minR45HitE_, jetUpdater_cfi::sort, TS4TS5LowerCut_, TS4TS5UpperCut_, and ~HLTHcalMETNoiseFilter().

                                                                            :
     HcalNoiseRBXCollectionTag_(iConfig.getParameter<edm::InputTag>("HcalNoiseRBXCollection")),
     severity_(iConfig.getParameter<int> ("severity")),
     maxNumRBXs_(iConfig.getParameter<int>("maxNumRBXs")),
     numRBXsToConsider_(iConfig.getParameter<int>("numRBXsToConsider")),
     needEMFCoincidence_(iConfig.getParameter<bool>("needEMFCoincidence")),
     minRBXEnergy_(iConfig.getParameter<double>("minRBXEnergy")),
     minRatio_(iConfig.getParameter<double>("minRatio")),
     maxRatio_(iConfig.getParameter<double>("maxRatio")),
     minHPDHits_(iConfig.getParameter<int>("minHPDHits")),
     minRBXHits_(iConfig.getParameter<int>("minRBXHits")),
     minHPDNoOtherHits_(iConfig.getParameter<int>("minHPDNoOtherHits")),
     minZeros_(iConfig.getParameter<int>("minZeros")),
     minHighEHitTime_(iConfig.getParameter<double>("minHighEHitTime")),
     maxHighEHitTime_(iConfig.getParameter<double>("maxHighEHitTime")),
     maxRBXEMF_(iConfig.getParameter<double>("maxRBXEMF")),
     minRecHitE_(iConfig.getParameter<double>("minRecHitE")),
     minLowHitE_(iConfig.getParameter<double>("minLowHitE")),
     minHighHitE_(iConfig.getParameter<double>("minHighHitE")),
     minR45HitE_(5.0),
     TS4TS5EnergyThreshold_(iConfig.getParameter<double>("TS4TS5EnergyThreshold"))
 {
 
   std::vector<double> TS4TS5UpperThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperThreshold");
   std::vector<double> TS4TS5UpperCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperCut");
   std::vector<double> TS4TS5LowerThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerThreshold");
   std::vector<double> TS4TS5LowerCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerCut");
 
   for(int i = 0; i < (int)TS4TS5UpperThresholdTemp.size() && i < (int)TS4TS5UpperCutTemp.size(); i++)
      TS4TS5UpperCut_.push_back(std::pair<double, double>(TS4TS5UpperThresholdTemp[i], TS4TS5UpperCutTemp[i]));
   sort(TS4TS5UpperCut_.begin(), TS4TS5UpperCut_.end());
 
   for(int i = 0; i < (int)TS4TS5LowerThresholdTemp.size() && i < (int)TS4TS5LowerCutTemp.size(); i++)
      TS4TS5LowerCut_.push_back(std::pair<double, double>(TS4TS5LowerThresholdTemp[i], TS4TS5LowerCutTemp[i]));
   sort(TS4TS5LowerCut_.begin(), TS4TS5LowerCut_.end());
 
   if(iConfig.existsAs<double>("minR45HitE"))
      minR45HitE_ = iConfig.getParameter<double>("minR45HitE");
 
   m_theHcalNoiseToken = consumes<reco::HcalNoiseRBXCollection>(HcalNoiseRBXCollectionTag_);
 }

HLTHcalMETNoiseFilter::~HLTHcalMETNoiseFilter ( )

overridedefault

Referenced by HLTHcalMETNoiseFilter().

Member Function Documentation

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

static

Definition at line 80 of file HLTHcalMETNoiseFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HBHEPulseShapeFlagSetter_cfi::TS4TS5LowerCut, HBHEPulseShapeFlagSetter_cfi::TS4TS5LowerThreshold, HBHEPulseShapeFlagSetter_cfi::TS4TS5UpperCut, and HBHEPulseShapeFlagSetter_cfi::TS4TS5UpperThreshold.

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("HcalNoiseRBXCollection",edm::InputTag("hltHcalNoiseInfoProducer"));
   desc.add<int>("severity",1);
   desc.add<int>("maxNumRBXs",2);
   desc.add<int>("numRBXsToConsider",2);
   desc.add<bool>("needEMFCoincidence",true);
   desc.add<double>("minRBXEnergy",50.0);
   desc.add<double>("minRatio",-999.);
   desc.add<double>("maxRatio",999.);
   desc.add<int>("minHPDHits",17);
   desc.add<int>("minRBXHits",999);
   desc.add<int>("minHPDNoOtherHits",10);
   desc.add<int>("minZeros",10);
   desc.add<double>("minHighEHitTime",-9999.0);
   desc.add<double>("maxHighEHitTime",9999.0);
   desc.add<double>("maxRBXEMF",0.02);
   desc.add<double>("minRecHitE",1.5);
   desc.add<double>("minLowHitE",10.0);
   desc.add<double>("minHighHitE",25.0);
   desc.add<double>("minR45HitE",5.0);
   desc.add<double>("TS4TS5EnergyThreshold",50.0);
 
   double TS4TS5UpperThresholdArray[5] = {70, 90, 100, 400, 4000 };
   double TS4TS5UpperCutArray[5] = {1, 0.8, 0.75, 0.72, 0.72};
   double TS4TS5LowerThresholdArray[7] = {100, 120, 150, 200, 300, 400, 500};
   double TS4TS5LowerCutArray[7] = {-1, -0.7, -0.4, -0.2, -0.08, 0, 0.1};
   std::vector<double> TS4TS5UpperThreshold(TS4TS5UpperThresholdArray, TS4TS5UpperThresholdArray+5);
   std::vector<double> TS4TS5UpperCut(TS4TS5UpperCutArray, TS4TS5UpperCutArray+5);
   std::vector<double> TS4TS5LowerThreshold(TS4TS5LowerThresholdArray, TS4TS5LowerThresholdArray+7);
   std::vector<double> TS4TS5LowerCut(TS4TS5LowerCutArray, TS4TS5LowerCutArray+7);
 
   desc.add<std::vector<double> >("TS4TS5UpperThreshold", TS4TS5UpperThreshold);
   desc.add<std::vector<double> >("TS4TS5UpperCut", TS4TS5UpperCut);
   desc.add<std::vector<double> >("TS4TS5LowerThreshold", TS4TS5LowerThreshold);
   desc.add<std::vector<double> >("TS4TS5LowerCut", TS4TS5LowerCut);
   descriptions.add("hltHcalMETNoiseFilter",desc);
 }

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

override

Definition at line 123 of file HLTHcalMETNoiseFilter.cc.

References edmIntegrityCheck::d, data, edm::Event::getByToken(), HcalNoiseRBXCollectionTag_, edm::HandleBase::isValid(), LogDebug, m_theHcalNoiseToken, maxHighEHitTime_, maxNumRBXs_, maxRatio_, maxRBXEMF_, minHighEHitTime_, minHighHitE_, minHPDHits_, minHPDNoOtherHits_, minLowHitE_, minR45HitE_, minRatio_, minRBXEnergy_, minRBXHits_, minRecHitE_, minZeros_, needEMFCoincidence_, numRBXsToConsider_, severity_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

 {
   using namespace reco;
 
   // in this case, do not filter anything
   if(severity_==0) return true;
 
   // get the RBXs produced by RecoMET/METProducers/HcalNoiseInfoProducer
   edm::Handle<HcalNoiseRBXCollection> rbxs_h;
   iEvent.getByToken(m_theHcalNoiseToken,rbxs_h);
   if(!rbxs_h.isValid()) {
     edm::LogError("DataNotFound") << "HLTHcalMETNoiseFilter: Could not find HcalNoiseRBXCollection product named "
                   << HcalNoiseRBXCollectionTag_ << "." << std::endl;
     return true;
   }
 
   // reject events with too many RBXs
   if(static_cast<int>(rbxs_h->size())>maxNumRBXs_) return true;
 
   // create a sorted set of the RBXs, ordered by energy
   noisedataset_t data;
   for(auto const & rbx : *rbxs_h) {
     CommonHcalNoiseRBXData d(rbx, minRecHitE_, minLowHitE_, minHighHitE_, TS4TS5EnergyThreshold_,
                  TS4TS5UpperCut_, TS4TS5LowerCut_, minR45HitE_);
     data.insert(d);
   }
 
   // data is now sorted by RBX energy
   // only consider top N=numRBXsToConsider_ energy RBXs
   int cntr=0;
   for(auto it=data.begin();
       it!=data.end() && cntr<numRBXsToConsider_;
       ++it, ++cntr) {
 
     bool passFilter=true;
     bool passEMF=true;
     if(it->energy()>minRBXEnergy_) {
       if(it->validRatio() && it->ratio()<minRatio_)        passFilter=false;
       else if(it->validRatio() && it->ratio()>maxRatio_)   passFilter=false;
       else if(it->numHPDHits()>=minHPDHits_)               passFilter=false;
       else if(it->numRBXHits()>=minRBXHits_)               passFilter=false;
       else if(it->numHPDNoOtherHits()>=minHPDNoOtherHits_) passFilter=false;
       else if(it->numZeros()>=minZeros_)                   passFilter=false;
       else if(it->minHighEHitTime()<minHighEHitTime_)      passFilter=false;
       else if(it->maxHighEHitTime()>maxHighEHitTime_)      passFilter=false;
       else if(!it->PassTS4TS5())                           passFilter=false;
 
       if(it->RBXEMF()<maxRBXEMF_) passEMF=false;
     }
 
     if((needEMFCoincidence_ && !passEMF && !passFilter) ||
        (!needEMFCoincidence_ && !passFilter)) {
       LogDebug("") << "HLTHcalMETNoiseFilter debug: Found a noisy RBX: "
            << "energy=" << it->energy() << "; "
            << "ratio=" << it->ratio() << "; "
            << "# RBX hits=" << it->numRBXHits() << "; "
            << "# HPD hits=" << it->numHPDHits() << "; "
            << "# Zeros=" << it->numZeros() << "; "
            << "min time=" << it->minHighEHitTime() << "; "
            << "max time=" << it->maxHighEHitTime() << "; "
            << "passTS4TS5=" << it->PassTS4TS5() << "; "
            << "RBX EMF=" << it->RBXEMF()
            << std::endl;
       return false;
     }
   }
 
   // no problems found
   return true;
 }

Member Data Documentation

edm::InputTag HLTHcalMETNoiseFilter::HcalNoiseRBXCollectionTag_

private

Definition at line 34 of file HLTHcalMETNoiseFilter.h.

Referenced by filter(), and HLTHcalMETNoiseFilter().

edm::EDGetTokenT<reco::HcalNoiseRBXCollection> HLTHcalMETNoiseFilter::m_theHcalNoiseToken

private

Definition at line 32 of file HLTHcalMETNoiseFilter.h.

Referenced by filter(), and HLTHcalMETNoiseFilter().

double HLTHcalMETNoiseFilter::maxHighEHitTime_

private

Definition at line 47 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::maxNumRBXs_

private

Definition at line 36 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::maxRatio_

private

Definition at line 41 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::maxRBXEMF_

private

Definition at line 48 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minHighEHitTime_

private

Definition at line 46 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minHighHitE_

private

Definition at line 51 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::minHPDHits_

private

Definition at line 42 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::minHPDNoOtherHits_

private

Definition at line 44 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minLowHitE_

private

Definition at line 51 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minR45HitE_

private

Definition at line 51 of file HLTHcalMETNoiseFilter.h.

Referenced by filter(), and HLTHcalMETNoiseFilter().

double HLTHcalMETNoiseFilter::minRatio_

private

Definition at line 40 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minRBXEnergy_

private

Definition at line 39 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::minRBXHits_

private

Definition at line 43 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::minRecHitE_

private

Definition at line 51 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::minZeros_

private

Definition at line 45 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

bool HLTHcalMETNoiseFilter::needEMFCoincidence_

private

Definition at line 38 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::numRBXsToConsider_

private

Definition at line 37 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

int HLTHcalMETNoiseFilter::severity_

private

Definition at line 35 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

double HLTHcalMETNoiseFilter::TS4TS5EnergyThreshold_

private

Definition at line 53 of file HLTHcalMETNoiseFilter.h.

Referenced by filter().

std::vector<std::pair<double, double> > HLTHcalMETNoiseFilter::TS4TS5LowerCut_

private

Definition at line 55 of file HLTHcalMETNoiseFilter.h.

Referenced by filter(), and HLTHcalMETNoiseFilter().

std::vector<std::pair<double, double> > HLTHcalMETNoiseFilter::TS4TS5UpperCut_

private

Definition at line 54 of file HLTHcalMETNoiseFilter.h.

Referenced by filter(), and HLTHcalMETNoiseFilter().

Classes

Public Member Functions

Static Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation