#include <HLTEgammaL1MatchFilterPairs.h>

Inheritance diagram for HLTEgammaL1MatchFilterPairs:

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

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

	~HLTEgammaL1MatchFilterPairs ()

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 Member Functions
bool	CheckL1Matching (edm::Ref< reco::RecoEcalCandidateCollection >ref, std::vector< l1extra::L1EmParticleRef > &l1EGIso, std::vector< l1extra::L1EmParticleRef > &l1EGNonIso) const

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 HLTEgammaL1MatchFilterPairs.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 33 of file HLTEgammaL1MatchFilterPairs.cc.

References AlsoNonIsolatedFirst_, AlsoNonIsolatedSecond_, barrel_end_, candIsolatedTag_, candIsolatedToken_, candNonIsolatedTag_, candNonIsolatedToken_, endcap_end_, edm::ParameterSet::getParameter(), l1IsolatedTag_, l1NonIsolatedTag_, L1SeedFilterTag_, L1SeedFilterToken_, region_eta_size_, region_eta_size_ecap_, and region_phi_size_.

                                                                                        : HLTFilter(iConfig)
 {
    candIsolatedTag_ = iConfig.getParameter< edm::InputTag > ("candIsolatedTag");
    l1IsolatedTag_ = iConfig.getParameter< edm::InputTag > ("l1IsolatedTag");
    candNonIsolatedTag_ = iConfig.getParameter< edm::InputTag > ("candNonIsolatedTag");
    l1NonIsolatedTag_ = iConfig.getParameter< edm::InputTag > ("l1NonIsolatedTag");
    L1SeedFilterTag_ = iConfig.getParameter< edm::InputTag > ("L1SeedFilterTag");
 
    AlsoNonIsolatedFirst_ = iConfig.getParameter<bool>("AlsoNonIsolatedFirst");
    AlsoNonIsolatedSecond_  = iConfig.getParameter<bool>("AlsoNonIsolatedSecond");
 
    region_eta_size_      = iConfig.getParameter<double> ("region_eta_size");
    region_eta_size_ecap_ = iConfig.getParameter<double> ("region_eta_size_ecap");
    region_phi_size_      = iConfig.getParameter<double> ("region_phi_size");
    barrel_end_           = iConfig.getParameter<double> ("barrel_end");
    endcap_end_           = iConfig.getParameter<double> ("endcap_end");
 
    candIsolatedToken_ = consumes<reco::RecoEcalCandidateCollection>(candIsolatedTag_);
    candNonIsolatedToken_ = consumes<reco::RecoEcalCandidateCollection>(candNonIsolatedTag_);
    L1SeedFilterToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(L1SeedFilterTag_);
 }

HLTEgammaL1MatchFilterPairs::~HLTEgammaL1MatchFilterPairs ( )

Definition at line 55 of file HLTEgammaL1MatchFilterPairs.cc.

55 {}

Member Function Documentation

bool HLTEgammaL1MatchFilterPairs::CheckL1Matching	(	edm::Ref< reco::RecoEcalCandidateCollection >	ref,
		std::vector< l1extra::L1EmParticleRef > &	l1EGIso,
		std::vector< l1extra::L1EmParticleRef > &	l1EGNonIso
	)		const

private

Definition at line 169 of file HLTEgammaL1MatchFilterPairs.cc.

References barrel_end_, i, M_PI, region_eta_size_, region_eta_size_ecap_, region_phi_size_, and TWOPI.

Referenced by hltFilter().

                                                                                                                                                                                                 {
 
   for (unsigned int i=0; i<l1EGIso.size(); i++) {
     //ORCA matching method
     double etaBinLow  = 0.;
     double etaBinHigh = 0.; 
     if(fabs(ref->eta()) < barrel_end_){
       etaBinLow = l1EGIso[i]->eta() - region_eta_size_/2.;
       etaBinHigh = etaBinLow + region_eta_size_;
     }
     else{
       etaBinLow = l1EGIso[i]->eta() - region_eta_size_ecap_/2.;
       etaBinHigh = etaBinLow + region_eta_size_ecap_;
     }
 
     float deltaphi=fabs(ref->phi() -l1EGIso[i]->phi());
     if(deltaphi>TWOPI) deltaphi-=TWOPI;
     if(deltaphi>M_PI) deltaphi=TWOPI-deltaphi;
 
     if(ref->eta() < etaBinHigh && ref->eta() > etaBinLow && deltaphi <region_phi_size_/2. )  {return true;}
 
   }
 
   for (unsigned int i=0; i<l1EGNonIso.size(); i++) {
     //ORCA matching method
     double etaBinLow  = 0.;
     double etaBinHigh = 0.; 
     if(fabs(ref->eta()) < barrel_end_){
       etaBinLow = l1EGNonIso[i]->eta() - region_eta_size_/2.;
       etaBinHigh = etaBinLow + region_eta_size_;
     }
     else{
       etaBinLow = l1EGNonIso[i]->eta() - region_eta_size_ecap_/2.;
       etaBinHigh = etaBinLow + region_eta_size_ecap_;
     }
 
     float deltaphi=fabs(ref->phi() - l1EGNonIso[i]->phi());
     if(deltaphi>TWOPI) deltaphi-=TWOPI;
     if(deltaphi>M_PI) deltaphi=TWOPI-deltaphi;
 
     if(ref->eta() < etaBinHigh && ref->eta() > etaBinLow && deltaphi <region_phi_size_/2. )  {return true;}
     
   }
 
   return false;
 }

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

static

Definition at line 58 of file HLTEgammaL1MatchFilterPairs.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>("candIsolatedTag",edm::InputTag("hltRecoIsolatedEcalCandidate"));
   desc.add<edm::InputTag>("l1IsolatedTag",edm::InputTag("l1extraParticles","Isolated"));
   desc.add<edm::InputTag>("candNonIsolatedTag",edm::InputTag("hltRecoNonIsolatedEcalCandidate"));
   desc.add<edm::InputTag>("l1NonIsolatedTag",edm::InputTag("l1extraParticles","NonIsolated"));
   desc.add<edm::InputTag>("L1SeedFilterTag",edm::InputTag("theL1SeedFilter"));
   desc.add<bool>("AlsoNonIsolatedFirst",false);
   desc.add<bool>("AlsoNonIsolatedSecond",false);
   desc.add<double>("region_eta_size",0.522);
   desc.add<double>("region_eta_size_ecap",1.0);
   desc.add<double>("region_phi_size",1.044);
   desc.add<double>("barrel_end",1.4791);
   desc.add<double>("endcap_end",2.65);
   descriptions.add("hltEgammaL1MatchFilterPairs",desc);
 }

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

overridevirtual

Implements HLTFilter.

Definition at line 80 of file HLTEgammaL1MatchFilterPairs.cc.

References accept(), trigger::TriggerRefsCollections::addObject(), AlsoNonIsolatedFirst_, AlsoNonIsolatedSecond_, candIsolatedToken_, candNonIsolatedToken_, CheckL1Matching(), HLT_25ns14e33_v1_cff::distance, edm::EventSetup::get(), edm::Event::getByToken(), L1SeedFilterToken_, gen::n, trigger::TriggerCluster, trigger::TriggerL1IsoEG, and trigger::TriggerL1NoIsoEG.

 {
 
   using namespace trigger;
   using namespace l1extra;
   std::vector < std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > > thePairs;
 
   edm::Handle<reco::RecoEcalCandidateCollection> recoIsolecalcands;
   iEvent.getByToken(candIsolatedToken_,recoIsolecalcands);
   edm::Handle<reco::RecoEcalCandidateCollection> recoNonIsolecalcands;
   iEvent.getByToken(candNonIsolatedToken_,recoNonIsolecalcands);
 
   // create pairs <L1Iso,L1Iso> and optionally <L1Iso, L1NonIso>
    for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand1= recoIsolecalcands->begin(); recoecalcand1!=recoIsolecalcands->end(); recoecalcand1++) {
      edm::Ref<reco::RecoEcalCandidateCollection> ref1 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoIsolecalcands, distance(recoIsolecalcands->begin(),recoecalcand1) );
        for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand2= recoIsolecalcands->begin(); recoecalcand2!=recoIsolecalcands->end(); recoecalcand2++) {
      edm::Ref<reco::RecoEcalCandidateCollection> ref2 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoIsolecalcands, distance(recoIsolecalcands->begin(),recoecalcand2) );
      if( &(*ref1) != &(*ref2) ) {thePairs.push_back(std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > (ref1,ref2) );}
        }
        if (AlsoNonIsolatedSecond_){
      for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand2= recoNonIsolecalcands->begin(); recoecalcand2!=recoNonIsolecalcands->end(); recoecalcand2++) {
        edm::Ref<reco::RecoEcalCandidateCollection> ref2 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoNonIsolecalcands, distance(recoNonIsolecalcands->begin(),recoecalcand2) );
        if( &(*ref1) != &(*ref2) ) {thePairs.push_back(std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > (ref1,ref2) );}
      }
        }
    }
 
 
   // create pairs <L1NonIso,L1Iso> and optionally <L1NonIso, L1NonIso>
    if (AlsoNonIsolatedFirst_){
      for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand1= recoNonIsolecalcands->begin(); recoecalcand1!=recoNonIsolecalcands->end(); recoecalcand1++) {
        edm::Ref<reco::RecoEcalCandidateCollection> ref1 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoNonIsolecalcands, distance(recoNonIsolecalcands->begin(),recoecalcand1) );
        for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand2= recoIsolecalcands->begin(); recoecalcand2!=recoIsolecalcands->end(); recoecalcand2++) {
      edm::Ref<reco::RecoEcalCandidateCollection> ref2 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoIsolecalcands, distance(recoIsolecalcands->begin(),recoecalcand2) );
      if( &(*ref1) != &(*ref2) ) {thePairs.push_back(std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > (ref1,ref2) );}
        }
        if (AlsoNonIsolatedSecond_){
      for (reco::RecoEcalCandidateCollection::const_iterator recoecalcand2= recoNonIsolecalcands->begin(); recoecalcand2!=recoNonIsolecalcands->end(); recoecalcand2++) {
        edm::Ref<reco::RecoEcalCandidateCollection> ref2 =  edm::Ref<reco::RecoEcalCandidateCollection>(recoNonIsolecalcands, distance(recoNonIsolecalcands->begin(),recoecalcand2) );
        if( &(*ref1) != &(*ref2) ) {thePairs.push_back(std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > (ref1,ref2) );}
      }
        }
      }
    }
 
 
    // Get the CaloGeometry
   edm::ESHandle<L1CaloGeometry> l1CaloGeom ;
   iSetup.get<L1CaloGeometryRecord>().get(l1CaloGeom) ;
 
 
   // look at all candidates,  check cuts and add to filter object
   int n(0);
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> L1SeedOutput;
   iEvent.getByToken (L1SeedFilterToken_,L1SeedOutput);
 
   std::vector<l1extra::L1EmParticleRef > l1EGIso;
   L1SeedOutput->getObjects(TriggerL1IsoEG, l1EGIso);
   std::vector<l1extra::L1EmParticleRef > l1EGNonIso;
   L1SeedOutput->getObjects(TriggerL1NoIsoEG, l1EGNonIso);
 
 //   std::cout<<"L1EGIso size: "<<l1EGIso.size()<<std::endl;
 //   for (unsigned int i=0; i<l1EGIso.size(); i++){std::cout<<"L1EGIso Et Eta phi: "<<l1EGIso[i]->et()<<" "<<l1EGIso[i]->eta()<<" "<<l1EGIso[i]->phi()<<std::endl;}
 //   std::cout<<"L1EGNonIso size: "<<l1EGNonIso.size()<<std::endl;
 //   for (unsigned int i=0; i<l1EGNonIso.size(); i++){std::cout<<"L1EGNonIso Et Eta phi: "<<l1EGNonIso[i]->et()<<" "<<l1EGNonIso[i]->eta()<<" "<<l1EGNonIso[i]->phi()<<std::endl;}
 //   std::cout<<"Lpair vector size: "<<thePairs.size()<<std::endl;
 
   std::vector < std::pair< edm::Ref<reco::RecoEcalCandidateCollection>, edm::Ref<reco::RecoEcalCandidateCollection> > >::iterator  pairsIt;
   for (pairsIt = thePairs.begin(); pairsIt != thePairs.end(); pairsIt++){
 //      edm::Ref<reco::RecoEcalCandidateCollection> r1 = pairsIt->first;
 //      edm::Ref<reco::RecoEcalCandidateCollection> r2 = pairsIt->second;
 //      std::cout<<"1) Et Eta phi: "<<r1->et()<<" "<<r1->eta()<<" "<<r1->phi()<<" 2) Et eta phi: "<<r2->et()<<" "<<r2->eta()<<" "<<r2->phi()<<std::endl;
 
     if ( CheckL1Matching(pairsIt->first,l1EGIso,l1EGNonIso) && CheckL1Matching(pairsIt->second,l1EGIso,l1EGNonIso) ){
       filterproduct.addObject(TriggerCluster, pairsIt->first);
       filterproduct.addObject(TriggerCluster, pairsIt->second);
       n++;
     }
   }
 
 
   //  std::cout<<"#####################################################"<<std::endl;
   // filter decision
   bool accept(n>=1);
 
   return accept;
 }

Member Data Documentation

bool HLTEgammaL1MatchFilterPairs::AlsoNonIsolatedFirst_

private

Definition at line 46 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs(), and hltFilter().

bool HLTEgammaL1MatchFilterPairs::AlsoNonIsolatedSecond_

private

Definition at line 46 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs(), and hltFilter().

double HLTEgammaL1MatchFilterPairs::barrel_end_

private

Definition at line 52 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by CheckL1Matching(), and HLTEgammaL1MatchFilterPairs().

edm::InputTag HLTEgammaL1MatchFilterPairs::candIsolatedTag_

private

Definition at line 37 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::EDGetTokenT<reco::RecoEcalCandidateCollection> HLTEgammaL1MatchFilterPairs::candIsolatedToken_

private

Definition at line 41 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs(), and hltFilter().

edm::InputTag HLTEgammaL1MatchFilterPairs::candNonIsolatedTag_

private

Definition at line 39 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::EDGetTokenT<reco::RecoEcalCandidateCollection> HLTEgammaL1MatchFilterPairs::candNonIsolatedToken_

private

Definition at line 42 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs(), and hltFilter().

double HLTEgammaL1MatchFilterPairs::endcap_end_

private

Definition at line 53 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::InputTag HLTEgammaL1MatchFilterPairs::l1IsolatedTag_

private

Definition at line 38 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::InputTag HLTEgammaL1MatchFilterPairs::l1NonIsolatedTag_

private

Definition at line 40 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::InputTag HLTEgammaL1MatchFilterPairs::L1SeedFilterTag_

private

Definition at line 44 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTEgammaL1MatchFilterPairs::L1SeedFilterToken_

private

Definition at line 45 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by HLTEgammaL1MatchFilterPairs(), and hltFilter().

double HLTEgammaL1MatchFilterPairs::region_eta_size_

private

Definition at line 49 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by CheckL1Matching(), and HLTEgammaL1MatchFilterPairs().

double HLTEgammaL1MatchFilterPairs::region_eta_size_ecap_

private

Definition at line 50 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by CheckL1Matching(), and HLTEgammaL1MatchFilterPairs().

double HLTEgammaL1MatchFilterPairs::region_phi_size_

private

Definition at line 51 of file HLTEgammaL1MatchFilterPairs.h.

Referenced by CheckL1Matching(), and HLTEgammaL1MatchFilterPairs().

Private Attributes
bool	AlsoNonIsolatedFirst_

bool	AlsoNonIsolatedSecond_

double	barrel_end_

edm::InputTag	candIsolatedTag_

edm::EDGetTokenT < reco::RecoEcalCandidateCollection >	candIsolatedToken_

edm::InputTag	candNonIsolatedTag_

edm::EDGetTokenT < reco::RecoEcalCandidateCollection >	candNonIsolatedToken_

double	endcap_end_

edm::InputTag	l1IsolatedTag_

edm::InputTag	l1NonIsolatedTag_

edm::InputTag	L1SeedFilterTag_

edm::EDGetTokenT < trigger::TriggerFilterObjectWithRefs >	L1SeedFilterToken_

double	region_eta_size_

double	region_eta_size_ecap_

double	region_phi_size_

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation