#include <HLTJetPairDzMatchFilter.h>

Inheritance diagram for HLTJetPairDzMatchFilter< T >:

Public Member Functions
bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override

	HLTJetPairDzMatchFilter (const edm::ParameterSet &)

	~HLTJetPairDzMatchFilter () override

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

	~HLTFilter () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

Public Member Functions inherited from edm::global::EDFilterBase
	EDFilterBase ()

ModuleDescription const &	moduleDescription () const

	~EDFilterBase () 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 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
double	m_jetMaxDZ

double	m_jetMaxEta

double	m_jetMinDR

double	m_jetMinPt

edm::InputTag	m_jetTag

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

int	m_triggerType

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

template<typename T>
class HLTJetPairDzMatchFilter< T >

Definition at line 19 of file HLTJetPairDzMatchFilter.h.

Constructor & Destructor Documentation

template<typename T >

HLTJetPairDzMatchFilter< T >::HLTJetPairDzMatchFilter ( const edm::ParameterSet & conf )

explicit

Definition at line 17 of file HLTJetPairDzMatchFilter.cc.

References edm::ParameterSet::getParameter(), HLTJetPairDzMatchFilter< T >::m_jetMaxDZ, HLTJetPairDzMatchFilter< T >::m_jetMaxEta, HLTJetPairDzMatchFilter< T >::m_jetMinDR, HLTJetPairDzMatchFilter< T >::m_jetMinPt, HLTJetPairDzMatchFilter< T >::m_jetTag, HLTJetPairDzMatchFilter< T >::m_jetToken, and HLTJetPairDzMatchFilter< T >::m_triggerType.

                                                                                : HLTFilter(conf)
 {
   m_jetTag  = conf.getParameter<edm::InputTag>("JetSrc");
   m_jetToken    = consumes<std::vector<T> >(m_jetTag);
   m_jetMinPt    = conf.getParameter<double>("JetMinPt");
   m_jetMaxEta   = conf.getParameter<double>("JetMaxEta");
   m_jetMinDR    = conf.getParameter<double>("JetMinDR");
   m_jetMaxDZ    = conf.getParameter<double>("JetMaxDZ");
   m_triggerType = conf.getParameter<int>("TriggerType");
 
   // set the minimum DR between jets, so that one never has a chance
   // to create a pair out of the same Jet replicated with two
   // different vertices
   if (m_jetMinDR < 0.1) m_jetMinDR = 0.1;
 
 }

template<typename T >

HLTJetPairDzMatchFilter< T >::~HLTJetPairDzMatchFilter ( )

overridedefault

Referenced by HLTJetPairDzMatchFilter< T >::fillDescriptions().

Member Function Documentation

template<typename T >

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

static

Definition at line 36 of file HLTJetPairDzMatchFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), HLTFilter::makeHLTFilterDescription(), and HLTJetPairDzMatchFilter< T >::~HLTJetPairDzMatchFilter().

                                                                                        {
   edm::ParameterSetDescription desc;
   makeHLTFilterDescription(desc);
   desc.add<edm::InputTag>("JetSrc",edm::InputTag("hltMatchL2Tau30ToPixelTrk5"));
   desc.add<double>("JetMinPt",25.0);
   desc.add<double>("JetMaxEta",2.4);
   desc.add<double>("JetMinDR",0.2);
   desc.add<double>("JetMaxDZ",0.2);
   desc.add<int>("TriggerType",84);
   descriptions.add(defaultModuleLabel<HLTJetPairDzMatchFilter<T>>(), desc);
 }

template<typename T >

bool HLTJetPairDzMatchFilter< T >::hltFilter	(	edm::Event &	ev,
		const edm::EventSetup &	es,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)		const

overridevirtual

Implements HLTFilter.

Definition at line 52 of file HLTJetPairDzMatchFilter.cc.

References funct::abs(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), reco::deltaPhi(), PVValHelper::dz, PVValHelper::eta, edm::Event::getByToken(), fwrapper::jets, HLTJetPairDzMatchFilter< T >::m_jetMaxDZ, HLTJetPairDzMatchFilter< T >::m_jetMaxEta, HLTJetPairDzMatchFilter< T >::m_jetMinDR, HLTJetPairDzMatchFilter< T >::m_jetMinPt, HLTJetPairDzMatchFilter< T >::m_jetTag, HLTJetPairDzMatchFilter< T >::m_jetToken, HLTJetPairDzMatchFilter< T >::m_triggerType, phi, EnergyCorrector::pt, and HLTFilter::saveTags().

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
 
   typedef vector<T> TCollection;
   typedef Ref<TCollection> TRef;
 
   // The resuilting filter object to store in the Event
 
   if (saveTags()) filterproduct.addCollectionTag(m_jetTag);
 
   // Ref to Candidate object to be recorded in the filter object
   TRef ref;
 
   // *** Pick up L2 tau jets which have been equipped with some meaningful vertices before that ***
 
   edm::Handle<TCollection> jetsHandle;
   ev.getByToken(m_jetToken, jetsHandle );
   const TCollection & jets = *jetsHandle;
   const size_t n_jets = jets.size();
 
   // *** Combine jets into pairs and check the dz matching ***
 
   size_t npairs = 0, nfail_dz = 0;
   if (n_jets > 1) for(size_t j1 = 0; j1 < n_jets; ++j1)
   {
     if ( jets[j1].pt() < m_jetMinPt || std::abs(jets[j1].eta()) > m_jetMaxEta ) continue;
 
     float mindz = 99.f;
     for(size_t j2 = j1+1; j2 < n_jets; ++j2)
     {
       if ( jets[j2].pt() < m_jetMinPt || std::abs(jets[j2].eta()) > m_jetMaxEta ) continue;
 
       float deta = jets[j1].eta() - jets[j2].eta();
       float dphi = reco::deltaPhi(jets[j1].phi(), jets[j2].phi());
       float dr2 = dphi*dphi + deta*deta;
       float dz = jets[j1].vz() - jets[j2].vz();
 
       // skip pairs of jets that are close
       if ( dr2 < m_jetMinDR*m_jetMinDR ) {
     continue;
       }
 
       if (std::abs(dz) < std::abs(mindz)) mindz = dz;
 
       // do not form a pair if dz is too large
       if ( std::abs(dz) > m_jetMaxDZ ) {
     ++nfail_dz;
     continue;
       }
 
       // add references to both jets
       ref = TRef(jetsHandle, j1);
       filterproduct.addObject(m_triggerType, ref);
 
       ref = TRef(jetsHandle, j2);
       filterproduct.addObject(m_triggerType, ref);
 
       ++npairs;
 
     }
 
   }
 
   return npairs>0;
 
 }