#include <HLT2jetGapFilter.h>

Inheritance diagram for HLT2jetGapFilter:

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

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

	~HLT2jetGapFilter ()

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

Private Attributes
edm::InputTag	inputTag_

double	minEt_

double	minEta_

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

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

Author: Monica Vazquez Acosta (CERN)

Definition at line 15 of file HLT2jetGapFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file HLT2jetGapFilter.cc.

References edm::ParameterSet::getParameter(), inputTag_, minEt_, and minEta_.

                                                                  : HLTFilter(iConfig) 
 {
    inputTag_ = iConfig.getParameter<edm::InputTag> ("inputTag");
    minEt_    = iConfig.getParameter<double> ("minEt");
    minEta_   = iConfig.getParameter<double> ("minEta"); 
 }

HLT2jetGapFilter::~HLT2jetGapFilter ( )

Definition at line 31 of file HLT2jetGapFilter.cc.

31 {}

Member Function Documentation

bool HLT2jetGapFilter::hltFilter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)

virtual

Implements HLTFilter.

Definition at line 36 of file HLT2jetGapFilter.cc.

References abs, accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), edm::Event::getByLabel(), inputTag_, minEt_, minEta_, n, HLTFilter::saveTags(), and trigger::TriggerJet.

 {
   using namespace trigger;
 
   // The filter object
   if (saveTags()) filterproduct.addCollectionTag(inputTag_);
 
   edm::Handle<reco::CaloJetCollection> recocalojets;
   iEvent.getByLabel(inputTag_,recocalojets);
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
 //  std::cout << "HLT2jetGapFilter " << recocalojets->size() << " jets in this event" << std::endl;
 
 
   if(recocalojets->size() > 1){
     // events with two or more jets
 
     double etjet1=0.;
     double etjet2=0.;
     double etajet1=0.;
     double etajet2=0.;
     int countjets =0;
 
     for (reco::CaloJetCollection::const_iterator recocalojet = recocalojets->begin(); 
      recocalojet<=(recocalojets->begin()+1); recocalojet++) {
       
       if(countjets==0) {
     etjet1 = recocalojet->et();
     etajet1 = recocalojet->eta();
       }
       if(countjets==1) {
     etjet2 = recocalojet->et();
     etajet2 = recocalojet->eta();
       }
       countjets++;
     }
 
   if(etjet1>minEt_ && etjet2>minEt_ && (etajet1*etajet2)<0 && std::abs(etajet1)>minEta_ && std::abs(etajet2)>minEta_){
       for (reco::CaloJetCollection::const_iterator recocalojet = recocalojets->begin(); 
        recocalojet<=(recocalojets->begin()+1); recocalojet++) {
     reco::CaloJetRef ref(reco::CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet)));
     filterproduct.addObject(TriggerJet,ref);
     n++;
       }
     }
     
   } // events with two or more jets
   
   
   
   // filter decision
   bool accept(n>=2);
   
   return accept;
 }