#include <HLTDiJetAveFilter.h>

Inheritance diagram for HLTDiJetAveFilter< T >:

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

virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	~HLTDiJetAveFilter ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module (edm::Event const &) const

const std::string *	moduleLabel () const

int	path (edm::Event const &) const

const std::string *	pathName (edm::Event const &) const

std::pair< int, int >	pmid (edm::Event const &) const

bool	saveTags () const

virtual	~HLTFilter ()

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

Public Member Functions inherited from edm::global::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 HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
edm::InputTag	inputJetTag_

edm::EDGetTokenT< std::vector < T > >	m_theJetToken

double	minDphi_

double	minPtAve_

double	minPtJet3_

int	triggerType_

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

template<typename T>
class HLTDiJetAveFilter< T >

Author: Dominique J. Mangeol; Leonard Apanasevich

Definition at line 23 of file HLTDiJetAveFilter.h.

Constructor & Destructor Documentation

template<typename T >

HLTDiJetAveFilter< T >::HLTDiJetAveFilter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 25 of file HLTDiJetAveFilter.cc.

References edm::InputTag::encode(), HLTDiJetAveFilter< T >::inputJetTag_, LogDebug, HLTDiJetAveFilter< T >::m_theJetToken, HLTDiJetAveFilter< T >::minDphi_, HLTDiJetAveFilter< T >::minPtAve_, HLTDiJetAveFilter< T >::minPtJet3_, and HLTDiJetAveFilter< T >::triggerType_.

                                                                       : HLTFilter(iConfig),
   inputJetTag_ (iConfig.template getParameter< edm::InputTag > ("inputJetTag")),
   minPtAve_    (iConfig.template getParameter<double> ("minPtAve")),
   minPtJet3_   (iConfig.template getParameter<double> ("minPtJet3")),
   minDphi_     (iConfig.template getParameter<double> ("minDphi")),
   triggerType_ (iConfig.template getParameter<int> ("triggerType"))
 {
   m_theJetToken = consumes<std::vector<T>>(inputJetTag_);
   LogDebug("") << "HLTDiJetAveFilter: Input/minPtAve/minPtJet3/minDphi/triggerType : "
            << inputJetTag_.encode() << " "
            << minPtAve_ << " "
            << minPtJet3_ << " "
            << minDphi_ << " "
            << triggerType_;
 }

template<typename T >

HLTDiJetAveFilter< T >::~HLTDiJetAveFilter ( )

Definition at line 42 of file HLTDiJetAveFilter.cc.

42 {}

Member Function Documentation

template<typename T >

void HLTDiJetAveFilter< T >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 46 of file HLTDiJetAveFilter.cc.

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

                                                                                  {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("inputJetTag",edm::InputTag("hltIterativeCone5CaloJets"));
   desc.add<double>("minPtAve",100.0);
   desc.add<double>("minPtJet3",99999.0);
   desc.add<double>("minDphi",-1.0);
   desc.add<int>("triggerType",trigger::TriggerJet);
   descriptions.add(defaultModuleLabel<HLTDiJetAveFilter<T>>(), desc);
 }

template<typename T >

bool HLTDiJetAveFilter< T >::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

Implements HLTFilter.

Definition at line 60 of file HLTDiJetAveFilter.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), SiPixelRawToDigiRegional_cfi::deltaPhi, HLT_25ns14e33_v1_cff::distance, edm::Event::getByToken(), i, gen::n, dt_dqm_sourceclient_common_cff::reco, and CommPDSkim_cfg::saveTags.

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   typedef vector<T> TCollection;
   typedef Ref<TCollection> TRef;
 
   // The filter object
   if (saveTags()) filterproduct.addCollectionTag(inputJetTag_);
 
   // get hold of collection of objects
   Handle<TCollection> objects;
   iEvent.getByToken (m_theJetToken,objects);
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
   if(objects->size() > 1){
     // events with two or more jets
 
     double ptjet1=0., ptjet2=0.,ptjet3=0.;
     double phijet1=0.,phijet2=0;
     int countjets =0;
 
     int nmax=1;
     if (objects->size() > 2) nmax=2;
 
     TRef JetRef1,JetRef2;
 
     typename TCollection::const_iterator i ( objects->begin() );
     for (; i<=(objects->begin()+nmax); i++) {
       if(countjets==0) {
     ptjet1 = i->pt();
     phijet1 = i->phi();
     JetRef1 = TRef(objects,distance(objects->begin(),i));
       }
       if(countjets==1) {
     ptjet2 = i->pt();
     phijet2 = i->phi();
     JetRef2 = TRef(objects,distance(objects->begin(),i));
       }
       if(countjets==2) {
     ptjet3 = i->pt();
       }
       ++countjets;
     }
 
     double PtAve=(ptjet1 + ptjet2) / 2.;
     double Dphi = std::abs(deltaPhi(phijet1,phijet2));
 
     if( PtAve>minPtAve_ && ptjet3<minPtJet3_ && Dphi>minDphi_){
       filterproduct.addObject(triggerType_,JetRef1);
       filterproduct.addObject(triggerType_,JetRef2);
       ++n;
     }
 
   } // events with two or more jets
 
 
 
   // filter decision
   bool accept(n>=1);
 
   return accept;
 }