#include <HLTHcalTowerNoiseCleaner.h>

Inheritance diagram for HLTHcalTowerNoiseCleaner:

Classes
struct	noisedatacomp

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

virtual void	produce (edm::Event &, const edm::EventSetup &)

	~HLTHcalTowerNoiseCleaner ()

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Private Attributes
edm::InputTag	HcalNoiseRBXCollectionTag_

edm::EDGetTokenT < CaloTowerCollection >	m_theCaloTowerCollectionToken

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_

edm::InputTag	TowerCollectionTag_

double	TS4TS5EnergyThreshold_

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

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

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

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

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: HLT filter module for cleaning HCal Noise from MET or MHT

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

Definition at line 23 of file HLTHcalTowerNoiseCleaner.h.

Member Typedef Documentation

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

private

Definition at line 65 of file HLTHcalTowerNoiseCleaner.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 50 of file HLTHcalTowerNoiseCleaner.cc.

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), HcalNoiseRBXCollectionTag_, i, m_theCaloTowerCollectionToken, m_theHcalNoiseToken, minR45HitE_, python.multivaluedict::sort(), TowerCollectionTag_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

   : HcalNoiseRBXCollectionTag_(iConfig.getParameter<edm::InputTag>("HcalNoiseRBXCollection")),
     TowerCollectionTag_(iConfig.getParameter<edm::InputTag>("CaloTowerCollection")),
     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());
 
   m_theHcalNoiseToken = consumes<reco::HcalNoiseRBXCollection>(HcalNoiseRBXCollectionTag_);
   m_theCaloTowerCollectionToken = consumes<CaloTowerCollection>(TowerCollectionTag_);
 
   if(iConfig.existsAs<double>("minR45HitE"))
      minR45HitE_ = iConfig.getParameter<double>("minR45HitE");
 
   produces<CaloTowerCollection>();
 }

HLTHcalTowerNoiseCleaner::~HLTHcalTowerNoiseCleaner ( )

Definition at line 97 of file HLTHcalTowerNoiseCleaner.cc.

97 {}

Member Function Documentation

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

static

Definition at line 100 of file HLTHcalTowerNoiseCleaner.cc.

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

                                                                                      {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("HcalNoiseRBXCollection",edm::InputTag("hltHcalNoiseInfoProducer"));
   desc.add<edm::InputTag>("CaloTowerCollection",edm::InputTag("hltTowerMakerForAll"));
   desc.add<double>("maxTowerNoiseEnergyFraction",0.5);
   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("hltHcalTowerNoiseCleaner",desc);
 }

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

virtual

Implements edm::EDProducer.

Definition at line 145 of file HLTHcalTowerNoiseCleaner.cc.

References edm::RefVector< C, T, F >::begin(), data, edm::RefVector< C, T, F >::end(), edm::Event::getByToken(), HcalNoiseRBXCollectionTag_, CaloTower::id(), edm::Ref< C, T, F >::id(), LogDebug, m_theCaloTowerCollectionToken, m_theHcalNoiseToken, maxHighEHitTime_, maxRatio_, maxRBXEMF_, minHighEHitTime_, minHighHitE_, minHPDHits_, minHPDNoOtherHits_, minLowHitE_, minR45HitE_, minRatio_, minRBXEnergy_, minRBXHits_, minRecHitE_, minZeros_, needEMFCoincidence_, edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, severity_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

 {
   using namespace reco;
 
   //get the calo MET / MHT
   edm::Handle<CaloTowerCollection> tower_h;
   iEvent.getByToken(m_theCaloTowerCollectionToken,tower_h);
   
   std::set<unsigned int> noisyTowers;
 
   if(not tower_h.isValid()){ //No towers MET, don't do anything and accept the event
     edm::LogError("HLTHcalTowerNoiseCleaner") << "Input Tower Collection is not Valid";
     return;
   }
 
   
   // get the RBXs produced by RecoMET/METProducers/HcalNoiseInfoProducer
   edm::Handle<HcalNoiseRBXCollection> rbxs_h;
   iEvent.getByToken(m_theHcalNoiseToken,rbxs_h);
   if(!rbxs_h.isValid()) {
     edm::LogWarning("HLTHcalTowerNoiseCleaner") << "Could not find HcalNoiseRBXCollection product named "
                           << HcalNoiseRBXCollectionTag_ << "." << std::endl;
     severity_=0;
   }
         
   // create a sorted set of the RBXs, ordered by energy
   noisedataset_t data;
   for(HcalNoiseRBXCollection::const_iterator it=rbxs_h->begin(); it!=rbxs_h->end(); ++it) {
     const HcalNoiseRBX &rbx=(*it);
     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
   if(severity_>0){
     for(noisedataset_t::const_iterator it=data.begin();
     it!=data.end();
     it++) {
       
       
       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)) { // check for noise
     LogDebug("") << "HLTHcalTowerNoiseCleaner 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;
     // add calotowers associated with this RBX to the noise list
     edm::RefVector<CaloTowerCollection> noiseTowers = it->rbxTowers();
     edm::RefVector<CaloTowerCollection>::const_iterator noiseTowersIt;
     //add these calotowers to the noisy list
     for( noiseTowersIt = noiseTowers.begin(); noiseTowersIt != noiseTowers.end(); noiseTowersIt++){
       edm::Ref<edm::SortedCollection<CaloTower> > tower_ref = *noiseTowersIt;
       CaloTowerDetId id = tower_ref->id();
       noisyTowers.insert( id.denseIndex() );
     }}
     } // done with noise loop
   }//if(severity_>0)
   
   //output collection
   std::auto_ptr<CaloTowerCollection> OutputTowers(new CaloTowerCollection() );
 
   CaloTowerCollection::const_iterator inTowersIt;
   
   for(inTowersIt = tower_h->begin(); inTowersIt != tower_h->end(); inTowersIt++){
     const CaloTower & tower = (*inTowersIt);
     CaloTowerDetId id = tower.id();
     if(noisyTowers.find( id.denseIndex() ) == noisyTowers.end()){ // the tower is not noisy
       OutputTowers->push_back(*inTowersIt);
     }
   }
   iEvent.put(OutputTowers);
 
 }

Member Data Documentation

edm::InputTag HLTHcalTowerNoiseCleaner::HcalNoiseRBXCollectionTag_

private

Definition at line 35 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

edm::EDGetTokenT<CaloTowerCollection> HLTHcalTowerNoiseCleaner::m_theCaloTowerCollectionToken

private

Definition at line 33 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

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

private

Definition at line 32 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

double HLTHcalTowerNoiseCleaner::maxHighEHitTime_

private

Definition at line 49 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::maxNumRBXs_

private

Definition at line 38 of file HLTHcalTowerNoiseCleaner.h.

double HLTHcalTowerNoiseCleaner::maxRatio_

private

Definition at line 43 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::maxRBXEMF_

private

Definition at line 50 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minHighEHitTime_

private

Definition at line 48 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minHighHitE_

private

Definition at line 53 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::minHPDHits_

private

Definition at line 44 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::minHPDNoOtherHits_

private

Definition at line 46 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minLowHitE_

private

Definition at line 53 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minR45HitE_

private

Definition at line 53 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

double HLTHcalTowerNoiseCleaner::minRatio_

private

Definition at line 42 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minRBXEnergy_

private

Definition at line 41 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::minRBXHits_

private

Definition at line 45 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

double HLTHcalTowerNoiseCleaner::minRecHitE_

private

Definition at line 53 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::minZeros_

private

Definition at line 47 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

bool HLTHcalTowerNoiseCleaner::needEMFCoincidence_

private

Definition at line 40 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

int HLTHcalTowerNoiseCleaner::numRBXsToConsider_

private

Definition at line 39 of file HLTHcalTowerNoiseCleaner.h.

int HLTHcalTowerNoiseCleaner::severity_

private

Definition at line 37 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

edm::InputTag HLTHcalTowerNoiseCleaner::TowerCollectionTag_

private

Definition at line 36 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner().

double HLTHcalTowerNoiseCleaner::TS4TS5EnergyThreshold_

private

Definition at line 55 of file HLTHcalTowerNoiseCleaner.h.

Referenced by produce().

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

private

Definition at line 57 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

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

private

Definition at line 56 of file HLTHcalTowerNoiseCleaner.h.

Referenced by HLTHcalTowerNoiseCleaner(), and produce().

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