#include <HLTEgammaDoubleEtFilter.h>

Inheritance diagram for HLTEgammaDoubleEtFilter:

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

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

	~HLTEgammaDoubleEtFilter ()

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_

edm::InputTag	L1IsoCollTag_

edm::InputTag	L1NonIsoCollTag_

int	npaircut_

bool	relaxed_

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: Monica Vazquez Acosta (CERN)

Definition at line 26 of file HLTEgammaDoubleEtFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 of file HLTEgammaDoubleEtFilter.cc.

References candTag_, candToken_, etcut1_, etcut2_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), L1IsoCollTag_, L1NonIsoCollTag_, npaircut_, and relaxed_.

                                                                                : HLTFilter(iConfig)
 {
   candTag_ = iConfig.getParameter< edm::InputTag > ("candTag");
   etcut1_  = iConfig.getParameter<double> ("etcut1");
   etcut2_  = iConfig.getParameter<double> ("etcut2");
   npaircut_  = iConfig.getParameter<int> ("npaircut");
   relaxed_ = iConfig.getUntrackedParameter<bool> ("relaxed",true) ;
   L1IsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1IsoCand");
   L1NonIsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1NonIsoCand");
   candToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(candTag_);
 }

HLTEgammaDoubleEtFilter::~HLTEgammaDoubleEtFilter ( )

Definition at line 42 of file HLTEgammaDoubleEtFilter.cc.

42 {}

Member Function Documentation

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

static

Definition at line 45 of file HLTEgammaDoubleEtFilter.cc.

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

                                                                                     {
    edm::ParameterSetDescription desc;
    makeHLTFilterDescription(desc);
    desc.add<edm::InputTag>("candTag",edm::InputTag("hltTrackIsolFilter"));
    desc.add<edm::InputTag>("L1IsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
    desc.add<edm::InputTag>("L1NonIsoCand",edm::InputTag("hltL1NonIsoRecoEcalCandidate"));
    desc.addUntracked<bool>("relaxed",true);
    desc.add<double>("etcut1", 30.0);
    desc.add<double>("etcut2", 20.0);
    desc.add<int>("npaircut", 1);
    descriptions.add("hltEgammaDoubleEtFilter",desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 60 of file HLTEgammaDoubleEtFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), candToken_, etcut1_, etcut2_, edm::Event::getByToken(), i, j, L1IsoCollTag_, L1NonIsoCollTag_, gen::n, npaircut_, relaxed_, HLTFilter::saveTags(), python.multivaluedict::sort(), trigger::TriggerCluster, and trigger::TriggerPhoton.

 {
   using namespace trigger;
 
   // The filter object
   if (saveTags()) {
     filterproduct.addCollectionTag(L1IsoCollTag_);
     if (relaxed_) filterproduct.addCollectionTag(L1NonIsoCollTag_);
   }
   // Ref to Candidate object to be recorded in filter object
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken (candToken_,PrevFilterOutput);
 
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> >  mysortedrecoecalcands;
   PrevFilterOutput->getObjects(TriggerPhoton,  mysortedrecoecalcands);
   if(mysortedrecoecalcands.empty()) PrevFilterOutput->getObjects(TriggerCluster,mysortedrecoecalcands);  //we dont know if its type trigger cluster or trigger photon
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
   // Sort the list
   std::sort(mysortedrecoecalcands.begin(), mysortedrecoecalcands.end(), EgammaHLTEtSortCriterium());
   edm::Ref<reco::RecoEcalCandidateCollection> ref1, ref2;
   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_ ){
       filterproduct.addObject(TriggerPhoton, ref1);
       filterproduct.addObject(TriggerPhoton, ref2);
       n++;
     }
       }
     }
   }
 
 
   // filter decision
   bool accept(n>=npaircut_);
 
   return accept;
 }