#include <HLTEgammaDoubleEtPhiFilter.h>

Inheritance diagram for HLTEgammaDoubleEtPhiFilter:

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

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

	~HLTEgammaDoubleEtPhiFilter ()

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	candTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	candToken_

double	etcut1_

double	etcut2_

double	max_Acop_

double	max_EtBalance_

double	min_Acop_

double	min_EtBalance_

int	npaircut_

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

Author: Jonathan Hollar (LLNL)

Definition at line 26 of file HLTEgammaDoubleEtPhiFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 28 of file HLTEgammaDoubleEtPhiFilter.cc.

References candTag_, candToken_, etcut1_, etcut2_, edm::ParameterSet::getParameter(), max_Acop_, max_EtBalance_, min_Acop_, min_EtBalance_, and npaircut_.

                                                                                      : HLTFilter(iConfig)
 {
    candTag_ = iConfig.getParameter< edm::InputTag > ("candTag");
    etcut1_  = iConfig.getParameter<double> ("etcut1");
    etcut2_  = iConfig.getParameter<double> ("etcut2");
    min_Acop_ =   iConfig.getParameter<double> ("MinAcop");
    max_Acop_ =   iConfig.getParameter<double> ("MaxAcop");
    min_EtBalance_ = iConfig.getParameter<double> ("MinEtBalance");
    max_EtBalance_ = iConfig.getParameter<double> ("MaxEtBalance");
    npaircut_  = iConfig.getParameter<int> ("npaircut");
    candToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(candTag_);
 }

HLTEgammaDoubleEtPhiFilter::~HLTEgammaDoubleEtPhiFilter ( )

Definition at line 41 of file HLTEgammaDoubleEtPhiFilter.cc.

41 {}

Member Function Documentation

void HLTEgammaDoubleEtPhiFilter::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 44 of file HLTEgammaDoubleEtPhiFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLT_25ns14e33_v1_cff::InputTag, and HLTFilter::makeHLTFilterDescription().

                                                                                        {
    edm::ParameterSetDescription desc;
    makeHLTFilterDescription(desc);
    desc.add<edm::InputTag>("candTag",edm::InputTag("hltDoubleL1MatchFilter"));
    desc.add<double>("etcut1", 6.0);
    desc.add<double>("etcut2", 6.0);
    desc.add<double>("MinAcop", -0.1);
    desc.add<double>("MaxAcop", 0.6);
    desc.add<double>("MinEtBalance", -1.0);
    desc.add<double>("MaxEtBalance", 10.0);
    desc.add<int>("npaircut", 1);
    descriptions.add("hltEgammaDoubleEtPhiFilter",desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 60 of file HLTEgammaDoubleEtPhiFilter.cc.

References accept(), trigger::TriggerRefsCollections::addObject(), candToken_, etcut1_, etcut2_, edm::Event::getByToken(), i, j, LogDebug, M_PI, max_Acop_, max_EtBalance_, min_Acop_, min_EtBalance_, gen::n, npaircut_, python.multivaluedict::sort(), trigger::TriggerCluster, and trigger::TriggerPhoton.

 {
   using namespace trigger;
 
   // Ref to Candidate object to be recorded in filter object
   edm::Ref<reco::RecoEcalCandidateCollection> ref;
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken (candToken_,PrevFilterOutput);
 
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> >  mysortedrecoecalcands;
   PrevFilterOutput->getObjects(TriggerCluster,  mysortedrecoecalcands);
   if(mysortedrecoecalcands.empty()) PrevFilterOutput->getObjects(TriggerPhoton,mysortedrecoecalcands);  //we dont know if its type trigger cluster or trigger photon
   // Sort the list
   std::sort(mysortedrecoecalcands.begin(), mysortedrecoecalcands.end(), EgammaHLTEtSortCriterium());
   edm::Ref<reco::RecoEcalCandidateCollection> ref1, ref2;
 
   int n(0);
   for (unsigned int i=0; i<mysortedrecoecalcands.size(); i++) {
     ref1 = mysortedrecoecalcands[i];
     if( ref1->et() >= etcut1_){
 
       for (unsigned int j=i+1; j<mysortedrecoecalcands.size(); j++) {
     ref2 = mysortedrecoecalcands[j];
     if( ref2->et() >= etcut2_ ){
     
       // Acoplanarity
       double acop = fabs(ref1->phi()-ref2->phi());
       if (acop>M_PI) acop = 2*M_PI - acop;
       acop = M_PI - acop;
       LogDebug("HLTEgammaDoubleEtPhiFilter") << " ... 1-2 acop= " << acop;
 
       if ((acop>=min_Acop_) && (acop<=max_Acop_))
       {
             // Et balance
             double etbalance = fabs(ref1->et()-ref2->et());
             if ((etbalance>=min_EtBalance_) && (etbalance<=max_EtBalance_))
         {
           filterproduct.addObject(TriggerCluster, ref1);
           filterproduct.addObject(TriggerCluster, ref2);
           n++;
         }
       }
     }
       }
     }
   }
 
 
   // filter decision
   bool accept(n>=npaircut_);
 
   return accept;
 }