CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Attributes
HLTEgammaDoubleLegCombFilter Class Reference

#include <HLTEgammaDoubleLegCombFilter.h>

Inheritance diagram for HLTEgammaDoubleLegCombFilter:
HLTFilter edm::global::EDFilter<> edm::global::EDFilterBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 HLTEgammaDoubleLegCombFilter (const edm::ParameterSet &)
 
bool hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
 
void matchCands (const std::vector< math::XYZPoint > &firstLegP3s, const std::vector< math::XYZPoint > &secondLegP3s, std::vector< std::pair< int, int > > &matchedCands) const
 
 ~HLTEgammaDoubleLegCombFilter () override
 
- 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
 
 ~HLTFilter () override
 
- Public Member Functions inherited from edm::global::EDFilter<>
 EDFilter ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsStreamLuminosityBlocks () const final
 
bool wantsStreamRuns () const final
 
- Public Member Functions inherited from edm::global::EDFilterBase
 EDFilterBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDFilterBase () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
static void getP3OfLegCands (const edm::Event &iEvent, const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > &filterToken, std::vector< math::XYZPoint > &p3s)
 
- 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 firstLegLastFilterTag_
 
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefsfirstLegLastFilterToken_
 
double maxMatchDR_
 
int nrRequiredFirstLeg_
 
int nrRequiredSecondLeg_
 
int nrRequiredUniqueLeg_
 
edm::InputTag secondLegLastFilterTag_
 
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefssecondLegLastFilterToken_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
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

$Id: HLTEgammaDoubleLegCombFilter.cc,

Author
Sam Harper (RAL)

Definition at line 19 of file HLTEgammaDoubleLegCombFilter.h.

Constructor & Destructor Documentation

HLTEgammaDoubleLegCombFilter::HLTEgammaDoubleLegCombFilter ( const edm::ParameterSet iConfig)
explicit

Definition at line 26 of file HLTEgammaDoubleLegCombFilter.cc.

References firstLegLastFilterTag_, firstLegLastFilterToken_, edm::ParameterSet::getParameter(), maxMatchDR_, nrRequiredFirstLeg_, nrRequiredSecondLeg_, nrRequiredUniqueLeg_, secondLegLastFilterTag_, and secondLegLastFilterToken_.

26  : HLTFilter(iConfig)
27 {
28  firstLegLastFilterTag_ = iConfig.getParameter<edm::InputTag>("firstLegLastFilter");
29  secondLegLastFilterTag_= iConfig.getParameter<edm::InputTag>("secondLegLastFilter");
30  nrRequiredFirstLeg_ = iConfig.getParameter<int> ("nrRequiredFirstLeg");
31  nrRequiredSecondLeg_ = iConfig.getParameter<int> ("nrRequiredSecondLeg");
32  nrRequiredUniqueLeg_ = iConfig.getParameter<int> ("nrRequiredUniqueLeg");
33  maxMatchDR_ = iConfig.getParameter<double> ("maxMatchDR");
34  firstLegLastFilterToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(firstLegLastFilterTag_);
35  secondLegLastFilterToken_ = consumes<trigger::TriggerFilterObjectWithRefs>(secondLegLastFilterTag_);
36 }
T getParameter(std::string const &) const
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > firstLegLastFilterToken_
HLTFilter(const edm::ParameterSet &config)
Definition: HLTFilter.cc:20
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > secondLegLastFilterToken_
HLTEgammaDoubleLegCombFilter::~HLTEgammaDoubleLegCombFilter ( )
overridedefault

Referenced by fillDescriptions().

Member Function Documentation

void HLTEgammaDoubleLegCombFilter::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 39 of file HLTEgammaDoubleLegCombFilter.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), HLTFilter::makeHLTFilterDescription(), and ~HLTEgammaDoubleLegCombFilter().

39  {
42  desc.add<edm::InputTag>("firstLegLastFilter",edm::InputTag("firstFilter"));
43  desc.add<edm::InputTag>("secondLegLastFilter",edm::InputTag("secondFilter"));
44  desc.add<int>("nrRequiredFirstLeg",0);
45  desc.add<int>("nrRequiredSecondLeg",0);
46  desc.add<int>("nrRequiredUniqueLeg",0);
47  desc.add<double>("maxMatchDR",0.01);
48  descriptions.add("hltEgammaDoubleLegCombFilter",desc);
49 }
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void makeHLTFilterDescription(edm::ParameterSetDescription &desc)
Definition: HLTFilter.cc:29
void add(std::string const &label, ParameterSetDescription const &psetDescription)
void HLTEgammaDoubleLegCombFilter::getP3OfLegCands ( const edm::Event iEvent,
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > &  filterToken,
std::vector< math::XYZPoint > &  p3s 
)
static

Definition at line 120 of file HLTEgammaDoubleLegCombFilter.cc.

References DEFINE_FWK_MODULE, edm::Event::getByToken(), trigger::TriggerRefsCollections::getObjects(), p3, trigger::TriggerCluster, trigger::TriggerElectron, trigger::TriggerJet, and trigger::TriggerPhoton.

Referenced by hltFilter().

121 {
123  iEvent.getByToken(filterToken,filterOutput);
124 
125  //its easier on the if statement flow if I try everything at once, shouldnt add to timing
126  std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > phoCands;
127  filterOutput->getObjects(trigger::TriggerPhoton,phoCands);
128  std::vector<edm::Ref<reco::RecoEcalCandidateCollection> > clusCands;
129  filterOutput->getObjects(trigger::TriggerCluster,clusCands);
130  std::vector<edm::Ref<reco::ElectronCollection> > eleCands;
131  filterOutput->getObjects(trigger::TriggerElectron,eleCands);
132  std::vector<edm::Ref<reco::CaloJetCollection> > jetCands;
133  filterOutput->getObjects(trigger::TriggerJet,jetCands);
134 
135  if(!phoCands.empty()){ //its photons
136  for(auto & phoCand : phoCands){
137  p3s.push_back(phoCand->superCluster()->position());
138  }
139  }else if(!clusCands.empty()){ //try trigger cluster (should never be this, at the time of writing (17/1/11) this would indicate an error)
140  for(auto & clusCand : clusCands){
141  p3s.push_back(clusCand->superCluster()->position());
142  }
143  }else if(!eleCands.empty()){
144  for(auto & eleCand : eleCands){
145  p3s.push_back(eleCand->superCluster()->position());
146  }
147  }else if(!jetCands.empty()){
148  for(auto & jetCand : jetCands){
150  p3.SetX(jetCand->p4().x());
151  p3.SetY(jetCand->p4().y());
152  p3.SetZ(jetCand->p4().z());
153  p3s.push_back(p3);
154  }
155  }
156 }
void getObjects(Vids &ids, VRphoton &refs) const
various physics-level getters:
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
double p3[4]
Definition: TauolaWrapper.h:91
bool HLTEgammaDoubleLegCombFilter::hltFilter ( edm::Event iEvent,
const edm::EventSetup iSetup,
trigger::TriggerFilterObjectWithRefs filterproduct 
) const
overridevirtual

Implements HLTFilter.

Definition at line 55 of file HLTEgammaDoubleLegCombFilter.cc.

References accept(), firstLegLastFilterToken_, getP3OfLegCands(), matchCands(), SiStripPI::max, nrRequiredFirstLeg_, nrRequiredSecondLeg_, nrRequiredUniqueLeg_, and secondLegLastFilterToken_.

56 {
57  //right, issue 1, we dont know if this is a TriggerElectron, TriggerPhoton, TriggerCluster (should never be a TriggerCluster btw as that implies the 4-vectors are not stored in AOD)
58 
59  //trigger::TriggerObjectType firstLegTrigType;
60  std::vector<math::XYZPoint> firstLegP3s;
61 
62  //trigger::TriggerObjectType secondLegTrigType;
63  std::vector<math::XYZPoint> secondLegP3s;
64 
65  getP3OfLegCands(iEvent,firstLegLastFilterToken_,firstLegP3s);
66  getP3OfLegCands(iEvent,secondLegLastFilterToken_,secondLegP3s);
67 
68  std::vector<std::pair<int,int> > matchedCands;
69  matchCands(firstLegP3s,secondLegP3s,matchedCands);
70 
71 
72  int nr1stLegOnly=0;
73  int nr2ndLegOnly=0;
74  int nrBoth=0;;
75 
76  for(auto & matchedCand : matchedCands){
77  if(matchedCand.first>=0){ //we found a first leg cand
78  if(matchedCand.second>=0) nrBoth++;//we also found a second leg cand
79  else nr1stLegOnly++; //we didnt find a second leg cand
80  }else if(matchedCand.second>=0) nr2ndLegOnly++; //we found a second leg cand but we didnt find a first leg
81 
82  }
83 
84  bool accept=false;
85  if(nr1stLegOnly + nr2ndLegOnly + nrBoth >= nrRequiredUniqueLeg_) {
86  if(nr1stLegOnly>=nrRequiredFirstLeg_ && nr2ndLegOnly>=nrRequiredSecondLeg_) accept=true;
87  else{
88  int nrNeededFirstLeg = std::max(0,nrRequiredFirstLeg_ - nr1stLegOnly);
89  int nrNeededSecondLeg = std::max(0,nrRequiredSecondLeg_ - nr2ndLegOnly);
90 
91  if(nrBoth >= nrNeededFirstLeg + nrNeededSecondLeg) accept = true;
92  }
93  }
94 
95  return accept;
96 }
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > firstLegLastFilterToken_
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:30
static void getP3OfLegCands(const edm::Event &iEvent, const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > &filterToken, std::vector< math::XYZPoint > &p3s)
void matchCands(const std::vector< math::XYZPoint > &firstLegP3s, const std::vector< math::XYZPoint > &secondLegP3s, std::vector< std::pair< int, int > > &matchedCands) const
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > secondLegLastFilterToken_
void HLTEgammaDoubleLegCombFilter::matchCands ( const std::vector< math::XYZPoint > &  firstLegP3s,
const std::vector< math::XYZPoint > &  secondLegP3s,
std::vector< std::pair< int, int > > &  matchedCands 
) const

Definition at line 99 of file HLTEgammaDoubleLegCombFilter.cc.

References reco::deltaR2(), maxMatchDR_, and jetUpdater_cfi::sort.

Referenced by hltFilter().

100 {
101  std::vector<size_t> matched2ndLegs;
102  for(size_t firstLegNr=0;firstLegNr<firstLegP3s.size();firstLegNr++){
103  int matchedNr = -1;
104  for(size_t secondLegNr=0;secondLegNr<secondLegP3s.size() && matchedNr==-1;secondLegNr++){
105  if(reco::deltaR2(firstLegP3s[firstLegNr],secondLegP3s[secondLegNr])<maxMatchDR_*maxMatchDR_) matchedNr=secondLegNr;
106  }
107  matchedCands.push_back(std::make_pair(firstLegNr,matchedNr));
108  if(matchedNr>=0) matched2ndLegs.push_back(static_cast<size_t>(matchedNr));
109  }
110  std::sort(matched2ndLegs.begin(),matched2ndLegs.end());
111 
112  for(size_t secondLegNr=0;secondLegNr<secondLegP3s.size();secondLegNr++){
113  if(!std::binary_search(matched2ndLegs.begin(),matched2ndLegs.end(),secondLegNr)){ //wasnt matched already
114  matchedCands.push_back(std::make_pair<int,int>(-1, static_cast<int>(secondLegNr)));
115  }
116  }
117 }
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:16

Member Data Documentation

edm::InputTag HLTEgammaDoubleLegCombFilter::firstLegLastFilterTag_
private

Definition at line 31 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTEgammaDoubleLegCombFilter::firstLegLastFilterToken_
private

Definition at line 33 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and hltFilter().

double HLTEgammaDoubleLegCombFilter::maxMatchDR_
private

Definition at line 38 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and matchCands().

int HLTEgammaDoubleLegCombFilter::nrRequiredFirstLeg_
private

Definition at line 35 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and hltFilter().

int HLTEgammaDoubleLegCombFilter::nrRequiredSecondLeg_
private

Definition at line 36 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and hltFilter().

int HLTEgammaDoubleLegCombFilter::nrRequiredUniqueLeg_
private

Definition at line 37 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and hltFilter().

edm::InputTag HLTEgammaDoubleLegCombFilter::secondLegLastFilterTag_
private

Definition at line 32 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter().

edm::EDGetTokenT<trigger::TriggerFilterObjectWithRefs> HLTEgammaDoubleLegCombFilter::secondLegLastFilterToken_
private

Definition at line 34 of file HLTEgammaDoubleLegCombFilter.h.

Referenced by HLTEgammaDoubleLegCombFilter(), and hltFilter().