#include <HLTAcoFilter.h>

Inheritance diagram for HLTAcoFilter:

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

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

	~HLTAcoFilter ()

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
std::string	AcoString_

edm::InputTag	inputJetTag_

edm::InputTag	inputMETTag_

double	maxDPhi_

double	minDPhi_

double	minEtjet1_

double	minEtjet2_

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: Dominique J. Mangeol

Definition at line 18 of file HLTAcoFilter.h.

Constructor & Destructor Documentation

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

explicit

\ class HLTAcoFilter

Author: Dominique J. Mangeol

acoplanar dphi(jet1,jet2),dphi(jet2,met),dphi(jet1,met)

Definition at line 29 of file HLTAcoFilter.cc.

References AcoString_, edm::ParameterSet::getParameter(), inputJetTag_, inputMETTag_, maxDPhi_, minDPhi_, minEtjet1_, and minEtjet2_.

                                                          : HLTFilter(iConfig) 
 {
    inputJetTag_ = iConfig.getParameter< edm::InputTag > ("inputJetTag");
    inputMETTag_ = iConfig.getParameter< edm::InputTag > ("inputMETTag");
    minDPhi_     = iConfig.getParameter<double> ("minDeltaPhi");
    maxDPhi_     = iConfig.getParameter<double> ("maxDeltaPhi");
    minEtjet1_   = iConfig.getParameter<double> ("minEtJet1"); 
    minEtjet2_   = iConfig.getParameter<double> ("minEtJet2"); 
    AcoString_   = iConfig.getParameter<std::string> ("Acoplanar");
 }

HLTAcoFilter::~HLTAcoFilter ( )

Definition at line 40 of file HLTAcoFilter.cc.

40 {}

Member Function Documentation

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

virtual

Implements HLTFilter.

Definition at line 45 of file HLTAcoFilter.cc.

References abs, accept(), AcoString_, trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), edm::Event::getByLabel(), inputJetTag_, inputMETTag_, M_PI, maxDPhi_, minDPhi_, minEtjet1_, minEtjet2_, n, dt_dqm_sourceclient_common_cff::reco, HLTFilter::saveTags(), trigger::TriggerJet, and trigger::TriggerMET.

 {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   // The filter object
   if (saveTags()) {
     filterproduct.addCollectionTag(inputJetTag_);
     filterproduct.addCollectionTag(inputMETTag_);
   }
 
   Handle<CaloJetCollection> recocalojets;
   iEvent.getByLabel(inputJetTag_,recocalojets);
   Handle<trigger::TriggerFilterObjectWithRefs> metcal;
   iEvent.getByLabel(inputMETTag_,metcal);
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
   int JetNum = recocalojets->size();
 
   // events with two or more jets
 
   double etjet1=0.;
   double etjet2=0.;
   double phijet1=0.;
   double phijet2=0.;
   //double etmiss=0.;
   double phimiss=0.;
    
   VRcalomet vrefMET; 
   metcal->getObjects(TriggerMET,vrefMET);
   CaloMETRef metRef=vrefMET.at(0);
   //etmiss  = vrefMET.at(0)->et();
   phimiss = vrefMET.at(0)->phi();
 
   CaloJetRef ref1,ref2;
 
   if (JetNum>0) {
     CaloJetCollection::const_iterator recocalojet = recocalojets->begin(); 
     
     etjet1 = recocalojet->et();
     phijet1 = recocalojet->phi();
     ref1  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));
     
     if(JetNum>1) {
       recocalojet++;
       etjet2 = recocalojet->et();
       phijet2 = recocalojet->phi();
       ref2  = CaloJetRef(recocalojets,distance(recocalojets->begin(),recocalojet));
     }
     double Dphi= -1.;
     int JetSel = 0;
     
     if (AcoString_ == "Jet2Met") {
       Dphi = std::abs(phimiss-phijet2);
       if (JetNum>=2 && etjet1>minEtjet1_  && etjet2>minEtjet2_) {JetSel=1;} 
     }
     if (AcoString_ == "Jet1Jet2") {
       Dphi = std::abs(phijet1-phijet2);
       if (JetNum>=2 && etjet1>minEtjet1_  && etjet2>minEtjet2_) {JetSel=1;} 
     }
     if (AcoString_ == "Jet1Met") { 
       Dphi = std::abs(phimiss-phijet1);
       if (JetNum>=1 && etjet1>minEtjet1_ ) {JetSel=1;} 
     }
     
       
     if (Dphi>M_PI) {Dphi=2.0*M_PI-Dphi;}
     if(JetSel>0 && Dphi>=minDPhi_ && Dphi<=maxDPhi_){
       
       if (AcoString_=="Jet2Met" || AcoString_=="Jet1Met")  {filterproduct.addObject(TriggerMET,metRef);}
       if (AcoString_=="Jet1Met" || AcoString_=="Jet1Jet2") {filterproduct.addObject(TriggerJet,ref1);}
       if (AcoString_=="Jet2Met" || AcoString_=="Jet1Jet2") {filterproduct.addObject(TriggerJet,ref2);}
       n++;
     }
     
 
   } // at least one jet
   
     
   // filter decision
   bool accept(n>=1);
     
   return accept;
 }