#include <HLTExclDiJetFilter.h>

Inheritance diagram for HLTExclDiJetFilter< T >:

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

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

	~HLTExclDiJetFilter ()

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	caloTowerTag_

bool	HF_OR_

edm::InputTag	inputJetTag_

edm::EDGetTokenT < CaloTowerCollection >	m_theCaloTowerCollectionToken

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

double	minHFe_

double	minPtJet_

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

Author: Dominique J. Mangeol; Leonard Apanasevich

Definition at line 24 of file HLTExclDiJetFilter.h.

Constructor & Destructor Documentation

template<typename T >

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

explicit

Definition at line 27 of file HLTExclDiJetFilter.cc.

References HLTExclDiJetFilter< T >::caloTowerTag_, edm::InputTag::encode(), HLTExclDiJetFilter< T >::HF_OR_, HLTExclDiJetFilter< T >::inputJetTag_, LogDebug, HLTExclDiJetFilter< T >::m_theCaloTowerCollectionToken, HLTExclDiJetFilter< T >::m_theJetToken, HLTExclDiJetFilter< T >::minHFe_, HLTExclDiJetFilter< T >::minPtJet_, and HLTExclDiJetFilter< T >::triggerType_.

                                                                         :
   HLTFilter(iConfig),
   inputJetTag_ (iConfig.template getParameter<edm::InputTag> ("inputJetTag")),
   caloTowerTag_(iConfig.template getParameter<edm::InputTag> ("caloTowerTag")),
   minPtJet_    (iConfig.template getParameter<double> ("minPtJet")),
   minHFe_      (iConfig.template getParameter<double> ("minHFe")),
   HF_OR_       (iConfig.template getParameter<bool> ("HF_OR")),
   triggerType_ (iConfig.template getParameter<int> ("triggerType"))
 {
   m_theJetToken = consumes<std::vector<T>>(inputJetTag_);
   m_theCaloTowerCollectionToken = consumes<CaloTowerCollection>(caloTowerTag_);
   LogDebug("") << "HLTExclDiJetFilter: Input/minPtJet/minHFe/HF_OR/triggerType : "
            << inputJetTag_.encode() << " "
            << caloTowerTag_.encode() << " "
            << minPtJet_ << " "
            << minHFe_ << " "
            << HF_OR_ << " "
            << triggerType_;
 }

template<typename T >

HLTExclDiJetFilter< T >::~HLTExclDiJetFilter ( )

Definition at line 48 of file HLTExclDiJetFilter.cc.

48 {}

Member Function Documentation

template<typename T >

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

static

Definition at line 52 of file HLTExclDiJetFilter.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("hltMCJetCorJetIcone5HF07"));
   desc.add<edm::InputTag>("caloTowerTag",edm::InputTag("hltTowerMakerForAll"));
   desc.add<double>("minPtJet",30.0);
   desc.add<double>("minHFe",50.0);
   desc.add<bool>("HF_OR",false);
   desc.add<int>("triggerType",trigger::TriggerJet);
   descriptions.add(defaultModuleLabel<HLTExclDiJetFilter<T>>(), desc);
 }

template<typename T >

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

overridevirtual

Implements HLTFilter.

Definition at line 67 of file HLTExclDiJetFilter.cc.

References funct::abs(), accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), HLT_25ns14e33_v1_cff::distance, edm::Event::getByToken(), M_PI, gen::n, python.connectstrParser::o, 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_);
 
   Handle<TCollection> recojets; //recojets can be any jet collections
   iEvent.getByToken(m_theJetToken,recojets);
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
   double ptjet1=0., ptjet2=0.;
   double phijet1=0., phijet2=0.;
 
   if(recojets->size() > 1){
     // events with two or more jets
 
     int countjets =0;
 
     TRef JetRef1,JetRef2;
 
     typename TCollection::const_iterator recojet ( recojets->begin() );
     for (;recojet<=(recojets->begin()+1); ++recojet) {
       //
       if(countjets==0) {
     ptjet1 = recojet->pt();
         phijet1 = recojet->phi();
     
     JetRef1 = TRef(recojets,distance(recojets->begin(),recojet));
       }
       //
       if(countjets==1) {
     ptjet2 = recojet->pt();
         phijet2 = recojet->phi();
 
     JetRef2 = TRef(recojets,distance(recojets->begin(),recojet));
       }
       //
       ++countjets;
     }
 
     if(ptjet1>minPtJet_ && ptjet2>minPtJet_ ){
       double Dphi=std::abs(phijet1-phijet2);
       if(Dphi>M_PI) Dphi=2.0*M_PI-Dphi;
       if(Dphi>0.5*M_PI) {
     filterproduct.addObject(triggerType_,JetRef1);
     filterproduct.addObject(triggerType_,JetRef2);
     ++n;
       }
     } // pt(jet1,jet2) > minPtJet_
 
   } // events with two or more jets
 
   // calotowers
   bool hf_accept=false;
 
   if(n>0) {
      double ehfp(0.);
      double ehfm(0.);
 
      Handle<CaloTowerCollection> o;
      iEvent.getByToken(m_theCaloTowerCollectionToken,o);
 //     if( o.isValid()) {
       for( CaloTowerCollection::const_iterator cc = o->begin(); cc != o->end(); ++cc ) {
        if(std::abs(cc->ieta())>28 && cc->energy()<4.0) continue;
         if(cc->ieta()>28)  ehfp+=cc->energy();  // HF+ energy
         if(cc->ieta()<-28) ehfm+=cc->energy();  // HF- energy
       }
  //    }
 
      bool hf_accept_and  = (ehfp<minHFe_) && (ehfm<minHFe_);
      bool hf_accept_or  = (ehfp<minHFe_) || (ehfm<minHFe_);
 
      hf_accept = HF_OR_ ? hf_accept_or : hf_accept_and;
 
   } // n>0
 
 
 
 // filter decision
   bool accept(n>0 && hf_accept);
 
   return accept;
 }