#include <HLTEgammaEtFilterPairs.h>

Inheritance diagram for HLTEgammaEtFilterPairs:

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

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

	~HLTEgammaEtFilterPairs ()

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
double	etcutEB1_

double	etcutEB2_

double	etcutEE1_

double	etcutEE2_

edm::InputTag	inputTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	inputToken_

edm::InputTag	L1IsoCollTag_

edm::InputTag	L1NonIsoCollTag_

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: Alessio Ghezzi

Definition at line 22 of file HLTEgammaEtFilterPairs.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 21 of file HLTEgammaEtFilterPairs.cc.

References etcutEB1_, etcutEB2_, etcutEE1_, etcutEE2_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inputTag_, inputToken_, L1IsoCollTag_, L1NonIsoCollTag_, and relaxed_.

                                                                              : HLTFilter(iConfig)
 {
    inputTag_ = iConfig.getParameter< edm::InputTag > ("inputTag");
    etcutEB1_  = iConfig.getParameter<double> ("etcut1EB");
    etcutEE1_  = iConfig.getParameter<double> ("etcut1EE");
    etcutEB2_  = iConfig.getParameter<double> ("etcut2EB");
    etcutEE2_  = iConfig.getParameter<double> ("etcut2EE");
    relaxed_ = iConfig.getUntrackedParameter<bool> ("relaxed",true) ;
    L1IsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1IsoCand");
    L1NonIsoCollTag_= iConfig.getParameter< edm::InputTag > ("L1NonIsoCand");
    inputToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(inputTag_);
 }

HLTEgammaEtFilterPairs::~HLTEgammaEtFilterPairs ( )

Definition at line 49 of file HLTEgammaEtFilterPairs.cc.

49 {}

Member Function Documentation

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

static

Definition at line 35 of file HLTEgammaEtFilterPairs.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>("inputTag",edm::InputTag("HLTEgammaL1MatchFilter"));
    desc.add<edm::InputTag>("L1IsoCand",edm::InputTag("hltL1IsoRecoEcalCandidate"));
    desc.add<edm::InputTag>("L1NonIsoCand",edm::InputTag("hltL1NonIsoRecoEcalCandidate"));
    desc.addUntracked<bool>("relaxed",true);
    desc.add<double>("etcut1EB", 1.0);
    desc.add<double>("etcut1EE", 1.0);
    desc.add<double>("etcut2EB", 1.0);
    desc.add<double>("etcut2EE", 1.0);
    descriptions.add("hltEgammaEtFilterPairs",desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 54 of file HLTEgammaEtFilterPairs.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), trigger::TriggerRefsCollections::addObject(), etcutEB1_, etcutEB2_, etcutEE1_, etcutEE2_, plotBeamSpotDB::first, edm::Event::getByToken(), i, inputToken_, L1IsoCollTag_, L1NonIsoCollTag_, gen::n, diffTwoXMLs::r1, diffTwoXMLs::r2, relaxed_, HLTFilter::saveTags(), edm::second(), trigger::TriggerCluster, and trigger::TriggerPhoton.

 {
   using namespace trigger;
   // The filter object
   if (saveTags()) {
     filterproduct.addCollectionTag(L1IsoCollTag_);
     if (relaxed_) filterproduct.addCollectionTag(L1NonIsoCollTag_);
   }
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> PrevFilterOutput;
   iEvent.getByToken (inputToken_,PrevFilterOutput);
 
   std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > recoecalcands;                // vref with your specific C++ collection type
   PrevFilterOutput->getObjects(TriggerCluster, recoecalcands);
   if(recoecalcands.empty()) PrevFilterOutput->getObjects(TriggerPhoton, recoecalcands);
   // they list should be interpreted as pairs:
   // <recoecalcands[0],recoecalcands[1]>
   // <recoecalcands[2],recoecalcands[3]>
   // <recoecalcands[4],recoecalcands[5]>
   // .......
 
   // Should I check that the size of recoecalcands is even ?
   int n(0);
 
   for (unsigned int i=0; i<recoecalcands.size(); i=i+2) {
 
      edm::Ref<reco::RecoEcalCandidateCollection> r1 = recoecalcands[i];
      edm::Ref<reco::RecoEcalCandidateCollection> r2 = recoecalcands[i+1];
      //  std::cout<<"EtFilter: 1) Et Eta phi: "<<r1->et()<<" "<<r1->eta()<<" "<<r1->phi()<<" 2) Et eta phi: "<<r2->et()<<" "<<r2->eta()<<" "<<r2->phi()<<std::endl;
      bool first  = (fabs(r1->eta()) < 1.479 &&  r1->et()  >= etcutEB1_) || (fabs(r1->eta()) >= 1.479 &&  r1->et()  >= etcutEE1_);
      bool second = (fabs(r2->eta()) < 1.479 &&  r2->et()  >= etcutEB2_) || (fabs(r2->eta()) >= 1.479 &&  r2->et()  >= etcutEE2_);
 
     if ( first && second ) {
       n++;
       filterproduct.addObject(TriggerCluster,r1 );
       filterproduct.addObject(TriggerCluster,r2 );
     }
   }
 
 
   // filter decision
   bool accept(n>=1);
 
   return accept;
 }