#include <HLTJetPairDzMatchFilter.h>

Inheritance diagram for HLTJetPairDzMatchFilter< T >:

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

	HLTJetPairDzMatchFilter (const edm::ParameterSet &)

	~HLTJetPairDzMatchFilter ()

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

int	module () const

const std::string *	moduleLabel () const

int	path () const

const std::string *	pathName () const

std::pair< int, int >	pmid () const

bool	saveTags () const

virtual	~HLTFilter ()

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Private Attributes
double	m_jetMaxDZ

double	m_jetMaxEta

double	m_jetMinDR

double	m_jetMinPt

edm::InputTag	m_jetSrc

int	m_triggerType

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

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

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 21 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_jetSrc, and HLTJetPairDzMatchFilter< T >::m_triggerType.

                                                                                : HLTFilter(conf)
 { 
   m_jetSrc  = conf.getParameter<edm::InputTag>("JetSrc");
   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 ( )

Definition at line 52 of file HLTJetPairDzMatchFilter.cc.

52 {}

Member Function Documentation

template<typename T >

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

static

Definition at line 39 of file HLTJetPairDzMatchFilter.cc.

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

                                                                                        {
   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(std::string("hlt")+std::string(typeid(HLTJetPairDzMatchFilter<T>).name()),desc);
 }

template<typename T >

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

virtual

Implements HLTFilter.

Definition at line 55 of file HLTJetPairDzMatchFilter.cc.

References abs, trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), reco::deltaPhi(), eta(), edm::Event::getByLabel(), fwrapper::jets, phi, dt_dqm_sourceclient_common_cff::reco, and CommPDSkim_cfg::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_jetSrc);
 
   // 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.getByLabel( m_jetSrc, 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;
 
 }