#include <HBHENoiseFilter.h>

Inheritance diagram for HBHENoiseFilter:

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

	~HBHENoiseFilter ()

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

virtual	~EDFilter ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)

Private Attributes
bool	IgnoreTS4TS5ifJetInLowBVRegion_

edm::InputTag	jetlabel_

double	maxHighEHitTime_

int	maxjetindex_

double	maxNHF_

double	maxRatio_

double	maxRBXEMF_

double	minHighEHitTime_

int	minHPDHits_

int	minHPDNoOtherHits_

double	minIsolatedNoiseSumE_

double	minIsolatedNoiseSumEt_

int	minNumIsolatedNoiseChannels_

double	minRatio_

int	minRBXHits_

int	minZeros_

edm::InputTag	noiselabel_

bool	useTS4TS5_

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

typedef WorkerT< EDFilter >	WorkerType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Description: Filter that identifies events containing an RBX with bad pulse-shape, timing, hit multiplicity, and ADC 0 counts Designed to reduce noise rate by factor of 100

Implementation: Use the HcalNoiseSummary to make cuts on an event-by-event basis

Definition at line 21 of file HBHENoiseFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 38 of file HBHENoiseFilter.cc.

References edm::ParameterSet::getParameter(), IgnoreTS4TS5ifJetInLowBVRegion_, jetlabel_, maxHighEHitTime_, maxjetindex_, maxNHF_, maxRatio_, maxRBXEMF_, minHighEHitTime_, minHPDHits_, minHPDNoOtherHits_, minIsolatedNoiseSumE_, minIsolatedNoiseSumEt_, minNumIsolatedNoiseChannels_, minRatio_, minRBXHits_, minZeros_, noiselabel_, and useTS4TS5_.

 {
   //now do what ever initialization is needed
   noiselabel_ = iConfig.getParameter<edm::InputTag>("noiselabel");
   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");
   minNumIsolatedNoiseChannels_ = iConfig.getParameter<int>("minNumIsolatedNoiseChannels");
   minIsolatedNoiseSumE_ = iConfig.getParameter<double>("minIsolatedNoiseSumE");
   minIsolatedNoiseSumEt_ = iConfig.getParameter<double>("minIsolatedNoiseSumEt");  
   useTS4TS5_ = iConfig.getParameter<bool>("useTS4TS5");
 
   IgnoreTS4TS5ifJetInLowBVRegion_ = iConfig.getParameter<bool>("IgnoreTS4TS5ifJetInLowBVRegion");
   jetlabel_ =  iConfig.getParameter<edm::InputTag>("jetlabel");
   maxjetindex_ = iConfig.getParameter<int>("maxjetindex");
   maxNHF_ = iConfig.getParameter<double>("maxNHF");
 
 }

HBHENoiseFilter::~HBHENoiseFilter ( )

Definition at line 64 of file HBHENoiseFilter.cc.

 {
 
 }

Member Function Documentation

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

privatevirtual

Implements edm::EDFilter.

Definition at line 76 of file HBHENoiseFilter.cc.

References edm::hlt::Exception, edm::Event::getByLabel(), IgnoreTS4TS5ifJetInLowBVRegion_, edm::HandleBase::isValid(), metsig::jet, jetlabel_, maxHighEHitTime_, maxjetindex_, maxNHF_, maxRatio_, maxRBXEMF_, minHighEHitTime_, minHPDHits_, minHPDNoOtherHits_, minIsolatedNoiseSumE_, minIsolatedNoiseSumEt_, minNumIsolatedNoiseChannels_, minRatio_, minRBXHits_, minZeros_, noiselabel_, edm::errors::ProductNotFound, edmLumisInFiles::summary, and useTS4TS5_.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

 {
    using namespace edm;
 
   // get the Noise summary object
   edm::Handle<HcalNoiseSummary> summary_h;
   iEvent.getByLabel(noiselabel_, summary_h);
   if(!summary_h.isValid()) {
     throw edm::Exception(edm::errors::ProductNotFound) << " could not find HcalNoiseSummary.\n";
     return true;
   }
   const HcalNoiseSummary summary = *summary_h;
 
   //  if(summary.HasBadRBXTS4TS5() == true) std::cout << "TS4TS5 rejection!" << std::endl;
   //  else                                  std::cout << "TS4TS5 passing!" << std::endl;
   
 
   // One HBHE region has HPD with low BV, which increases hits in HPD, and
   // often causes otherwise good events to be flagged as noise.
 
   bool goodJetFoundInLowBVRegion=false; // checks whether a jet is in a low BV region, where false noise flagging rate is higher.
   if ( IgnoreTS4TS5ifJetInLowBVRegion_==true)
     {
       edm::Handle<reco::PFJetCollection> pfjet_h;
       iEvent.getByLabel(jetlabel_, pfjet_h);
       if (pfjet_h.isValid()) 
     {
       // Loop over all jets up to (and including) maxjetindex.
       // If jet found in low-BV region, set goodJetFoundInLowBVRegion = true
       int jetindex=0;
       for( reco::PFJetCollection::const_iterator jet = pfjet_h->begin(); jet != pfjet_h->end(); ++jet) 
         {
           if (jetindex>maxjetindex_) break; // only look at first N jets (N specified by user)
           // Check whether jet is in low-BV region (0<eta<1.4, -1.8<phi<-1.4)
           if (jet->eta()>0 && jet->eta()<1.4 && 
           jet->phi()>-1.8 && jet->phi()<-1.4) 
         {
           if (maxNHF_<0 || (maxNHF_>0 &&  jet->neutralHadronEnergyFraction()<maxNHF_))
             {
               //std::cout <<"JET eta = "<<jet->eta()<<"  phi = "<<jet->phi()<<"  Energy = "<<jet->energy()<<"  neutral had fraction = "<<jet->neutralHadronEnergy()/jet->energy() <<std::endl;
               goodJetFoundInLowBVRegion=true;
               break;
             }
         }
           ++jetindex;
         }
     } //if (pfjet_h.isValid())
       else // no valid jet collection found
     {
       // If no jet collection found, do we want to throw a fatal exception?  Or just proceed normally, not treating the lowBV region as special?
       //throw edm::Exception(edm::errors::ProductNotFound) << " could not find PFJetCollection with label "<<jetlabel_<<".\n";
     }
     } // if (IgnoreTS4TS5ifJetInLowBVRegion_==true)
 
   if(summary.minE2Over10TS()<minRatio_) return false;
   if(summary.maxE2Over10TS()>maxRatio_) return false;
   if(summary.maxHPDHits()>=minHPDHits_) return false;
   if(summary.maxRBXHits()>=minRBXHits_) return false;
   if(summary.maxHPDNoOtherHits()>=minHPDNoOtherHits_) return false;
   if(summary.maxZeros()>=minZeros_) return false;
   if(summary.min25GeVHitTime()<minHighEHitTime_) return false;
   if(summary.max25GeVHitTime()>maxHighEHitTime_) return false;
   if(summary.minRBXEMF()<maxRBXEMF_) return false;
   if(summary.numIsolatedNoiseChannels()>=minNumIsolatedNoiseChannels_) return false;
   if(summary.isolatedNoiseSumE()>=minIsolatedNoiseSumE_) return false;
   if(summary.isolatedNoiseSumEt()>=minIsolatedNoiseSumEt_) return false;
   // Only use TS4TS5 test if jet is not in low BV region
   if(useTS4TS5_ == true && summary.HasBadRBXTS4TS5() == true && goodJetFoundInLowBVRegion==false) return false;
 
   return true;
 }