CMS 3D CMS Logo
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
reco::PhysObjectMatcher< C1, C2, S, D, Q > Class Template Reference

#include <PhysObjectMatcher.h>

Inheritance diagram for reco::PhysObjectMatcher< C1, C2, S, D, Q >:
edm::stream::EDProducer<>

Public Member Functions

 PhysObjectMatcher (const edm::ParameterSet &cfg)
 
 ~PhysObjectMatcher () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 

Private Types

typedef std::pair< size_t, size_t > IndexPair
 
typedef std::vector< IndexPair > MatchContainer
 
typedef edm::Association< C2 > MatchMap
 
typedef C1::value_type T1
 
typedef C2::value_type T2
 

Private Member Functions

void produce (edm::Event &, const edm::EventSetup &) override
 
bool select (const T1 &c1, const T2 &c2) const
 

Private Attributes

edm::ParameterSet config_
 
D distance_
 
edm::EDGetTokenT< C2 > matchedToken_
 
bool resolveAmbiguities_
 
bool resolveByMatchQuality_
 
S select_
 
edm::EDGetTokenT< C1 > srcToken_
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

template<typename C1, typename C2, typename S, typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
class reco::PhysObjectMatcher< C1, C2, S, D, Q >

Definition at line 58 of file PhysObjectMatcher.h.

Member Typedef Documentation

◆ IndexPair

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
typedef std::pair<size_t, size_t> reco::PhysObjectMatcher< C1, C2, S, D, Q >::IndexPair
private

Definition at line 67 of file PhysObjectMatcher.h.

◆ MatchContainer

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
typedef std::vector<IndexPair> reco::PhysObjectMatcher< C1, C2, S, D, Q >::MatchContainer
private

Definition at line 68 of file PhysObjectMatcher.h.

◆ MatchMap

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
typedef edm::Association<C2> reco::PhysObjectMatcher< C1, C2, S, D, Q >::MatchMap
private

Definition at line 66 of file PhysObjectMatcher.h.

◆ T1

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
typedef C1::value_type reco::PhysObjectMatcher< C1, C2, S, D, Q >::T1
private

Definition at line 64 of file PhysObjectMatcher.h.

◆ T2

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
typedef C2::value_type reco::PhysObjectMatcher< C1, C2, S, D, Q >::T2
private

Definition at line 65 of file PhysObjectMatcher.h.

Constructor & Destructor Documentation

◆ PhysObjectMatcher()

template<typename C1 , typename C2 , typename S , typename D , typename Q >
reco::PhysObjectMatcher< C1, C2, S, D, Q >::PhysObjectMatcher ( const edm::ParameterSet &  cfg)

Definition at line 85 of file PhysObjectMatcher.h.

86  : config_(cfg),
87  srcToken_(consumes<C1>(cfg.template getParameter<edm::InputTag>("src"))),
88  matchedToken_(consumes<C2>(cfg.template getParameter<edm::InputTag>("matched"))),
89  resolveAmbiguities_(cfg.template getParameter<bool>("resolveAmbiguities")),
90  resolveByMatchQuality_(cfg.template getParameter<bool>("resolveByMatchQuality")),
91  select_(reco::modules::make<S>(cfg)),
92  distance_(reco::modules::make<D>(cfg)) {
93  // definition of the product
94  produces<MatchMap>();
95  // set resolveByMatchQuality only if ambiguities are to be resolved
96  resolveByMatchQuality_ = resolveByMatchQuality_ && resolveAmbiguities_;
97  }

References reco::PhysObjectMatcher< C1, C2, S, D, Q >::resolveAmbiguities_, and reco::PhysObjectMatcher< C1, C2, S, D, Q >::resolveByMatchQuality_.

◆ ~PhysObjectMatcher()

template<typename C1 , typename C2 , typename S , typename D , typename Q >
reco::PhysObjectMatcher< C1, C2, S, D, Q >::~PhysObjectMatcher ( )
override

Definition at line 100 of file PhysObjectMatcher.h.

100 {}

Member Function Documentation

◆ produce()

template<typename C1 , typename C2 , typename S , typename D , typename Q >
void reco::PhysObjectMatcher< C1, C2, S, D, Q >::produce ( edm::Event &  evt,
const edm::EventSetup &   
)
overrideprivate

Definition at line 103 of file PhysObjectMatcher.h.

103  {
104  using namespace edm;
105  using namespace std;
106  typedef std::pair<size_t, size_t> IndexPair;
107  typedef std::vector<IndexPair> MatchContainer;
108  // get collections from event
109  Handle<C2> matched;
110  evt.getByToken(matchedToken_, matched);
111  Handle<C1> cands;
112  evt.getByToken(srcToken_, cands);
113  // create product
114  unique_ptr<MatchMap> matchMap(new MatchMap(matched));
115  size_t size = cands->size();
116  if (size != 0) {
117  //
118  // create helpers
119  //
120  Q comparator(config_, *cands, *matched);
121  typename MatchMap::Filler filler(*matchMap);
122  ::helper::MasterCollection<C1> master(cands, evt);
123  vector<int> indices(master.size(), -1); // result: indices in target collection
124  vector<bool> mLock(matched->size(), false); // locks in target collection
125  MatchContainer matchPairs; // container of matched pairs
126  // loop over candidates
127  for (size_t c = 0; c != size; ++c) {
128  const T1& cand = (*cands)[c];
129  // no global comparison of match quality -> reset the container for each candidate
130  if (!resolveByMatchQuality_)
131  matchPairs.clear();
132  // loop over target collection
133  for (size_t m = 0; m != matched->size(); ++m) {
134  const T2& match = (*matched)[m];
135  // check lock and preselection
136  if (!mLock[m] && select(cand, match)) {
137  // double dist = testDR_(cand,match);
138  // cout << "dist between c = " << c << " and m = "
139  // << m << " is " << dist << " at pts of "
140  // << cand.pt() << " " << match.pt() << endl;
141  // matching requirement fulfilled -> store pair of indices
142  if (distance_(cand, match))
143  matchPairs.push_back(make_pair(c, m));
144  }
145  }
146  // if match(es) found and no global ambiguity resolution requested
147  if (!matchPairs.empty() && !resolveByMatchQuality_) {
148  // look for and store best match
149  size_t idx = master.index(c);
150  assert(idx < indices.size());
151  size_t index = min_element(matchPairs.begin(), matchPairs.end(), comparator)->second;
152  indices[idx] = index;
153  // if ambiguity resolution by order of (reco) candidates:
154  // lock element in target collection
155  if (resolveAmbiguities_)
156  mLock[index] = true;
157  // {
158  // MatchContainer::const_iterator i = min_element(matchPairs.begin(), matchPairs.end(), comparator);
159  // cout << "smallest distance for c = " << c << " is "
160  // << testDR_((*cands)[(*i).first],
161  // (*matched)[(*i).second]) << endl;
162  // }
163  }
164  }
165  // ambiguity resolution by global match quality (if requested)
166  if (resolveByMatchQuality_) {
167  // sort container of all matches by quality
168  sort(matchPairs.begin(), matchPairs.end(), comparator);
169  vector<bool> cLock(master.size(), false);
170  // loop over sorted container
171  for (MatchContainer::const_iterator i = matchPairs.begin(); i != matchPairs.end(); ++i) {
172  size_t c = (*i).first;
173  size_t m = (*i).second;
174  // cout << "rel dp = " << ((*cands)[c].pt()-(*matched)[m].pt())/(*matched)[m].pt() << endl;
175  // accept only pairs without any lock
176  if (mLock[m] || cLock[c])
177  continue;
178  // store index to target collection and lock the two items
179  size_t idx = master.index(c);
180  assert(idx < indices.size());
181  indices[idx] = m;
182  mLock[m] = true;
183  cLock[c] = true;
184  }
185  }
186  filler.insert(master.get(), indices.begin(), indices.end());
187  filler.fill();
188  }
189  evt.put(std::move(matchMap));
190  }

References cms::cuda::assert(), HltBtagPostValidation_cff::c, HLT_FULL_cff::cands, distance_, trigObjTnPSource_cfi::filler, edm::Event::getByToken(), mps_fire::i, heavyIonCSV_trainingSettings::idx, bTagCombinedSVVariables_cff::indices, visualization-live-secondInstance_cfg::m, volumeBasedMagneticField_160812_cfi::master, match(), muonTagProbeFilters_cff::matched, eostools::move(), edm::Event::put(), class-composition::Q, edm::second(), singleTopDQM_cfi::select, and findQualityFiles::size.

◆ select()

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
bool reco::PhysObjectMatcher< C1, C2, S, D, Q >::select ( const T1 &  c1,
const T2 &  c2 
) const
inlineprivate

Definition at line 78 of file PhysObjectMatcher.h.

78 { return select_(c1, c2); }

References alignmentValidation::c1, and reco::PhysObjectMatcher< C1, C2, S, D, Q >::select_.

Member Data Documentation

◆ config_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
edm::ParameterSet reco::PhysObjectMatcher< C1, C2, S, D, Q >::config_
private

Definition at line 70 of file PhysObjectMatcher.h.

◆ distance_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
D reco::PhysObjectMatcher< C1, C2, S, D, Q >::distance_
private

Definition at line 80 of file PhysObjectMatcher.h.

◆ matchedToken_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
edm::EDGetTokenT<C2> reco::PhysObjectMatcher< C1, C2, S, D, Q >::matchedToken_
private

Definition at line 72 of file PhysObjectMatcher.h.

◆ resolveAmbiguities_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
bool reco::PhysObjectMatcher< C1, C2, S, D, Q >::resolveAmbiguities_
private

Definition at line 73 of file PhysObjectMatcher.h.

Referenced by reco::PhysObjectMatcher< C1, C2, S, D, Q >::PhysObjectMatcher().

◆ resolveByMatchQuality_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
bool reco::PhysObjectMatcher< C1, C2, S, D, Q >::resolveByMatchQuality_
private

Definition at line 75 of file PhysObjectMatcher.h.

Referenced by reco::PhysObjectMatcher< C1, C2, S, D, Q >::PhysObjectMatcher().

◆ select_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
S reco::PhysObjectMatcher< C1, C2, S, D, Q >::select_
private

Definition at line 79 of file PhysObjectMatcher.h.

Referenced by reco::PhysObjectMatcher< C1, C2, S, D, Q >::select().

◆ srcToken_

template<typename C1 , typename C2 , typename S , typename D = reco::MatchByDR<typename C1::value_type, typename C2::value_type>, typename Q = helper::LessByMatchDistance<DeltaR<typename C1::value_type, typename C2::value_type>, C1, C2>>
edm::EDGetTokenT<C1> reco::PhysObjectMatcher< C1, C2, S, D, Q >::srcToken_
private

Definition at line 71 of file PhysObjectMatcher.h.

reco::PhysObjectMatcher::MatchMap
edm::Association< C2 > MatchMap
Definition: PhysObjectMatcher.h:66
muonTagProbeFilters_cff.matched
matched
Definition: muonTagProbeFilters_cff.py:62
bTagCombinedSVVariables_cff.indices
indices
Definition: bTagCombinedSVVariables_cff.py:67
reco::PhysObjectMatcher::select_
S select_
Definition: PhysObjectMatcher.h:79
mps_fire.i
i
Definition: mps_fire.py:428
edm
HLT enums.
Definition: AlignableModifier.h:19
reco::PhysObjectMatcher::T2
C2::value_type T2
Definition: PhysObjectMatcher.h:65
reco::PhysObjectMatcher::T1
C1::value_type T1
Definition: PhysObjectMatcher.h:64
reco::PhysObjectMatcher::matchedToken_
edm::EDGetTokenT< C2 > matchedToken_
Definition: PhysObjectMatcher.h:72
cms::cuda::assert
assert(be >=bs)
reco::PhysObjectMatcher::distance_
D distance_
Definition: PhysObjectMatcher.h:80
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
reco::PhysObjectMatcher::select
bool select(const T1 &c1, const T2 &c2) const
Definition: PhysObjectMatcher.h:78
edm::Handle
Definition: AssociativeIterator.h:50
class-composition.Q
Q
Definition: class-composition.py:82
heavyIonCSV_trainingSettings.idx
idx
Definition: heavyIonCSV_trainingSettings.py:5
visualization-live-secondInstance_cfg.m
m
Definition: visualization-live-secondInstance_cfg.py:72
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:531
reco::PhysObjectMatcher::MatchContainer
std::vector< IndexPair > MatchContainer
Definition: PhysObjectMatcher.h:68
reco::PhysObjectMatcher::srcToken_
edm::EDGetTokenT< C1 > srcToken_
Definition: PhysObjectMatcher.h:71
HLT_FULL_cff.cands
cands
Definition: HLT_FULL_cff.py:15142
trigObjTnPSource_cfi.filler
filler
Definition: trigObjTnPSource_cfi.py:21
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
cand
Definition: decayParser.h:32
alignmentValidation.c1
c1
do drawing
Definition: alignmentValidation.py:1025
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
reco::PhysObjectMatcher::resolveAmbiguities_
bool resolveAmbiguities_
Definition: PhysObjectMatcher.h:73
edm::ValueMap< int >::Filler
helper::Filler< ValueMap< int > > Filler
Definition: ValueMap.h:168
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
looper.cfg
cfg
Definition: looper.py:297
eostools.move
def move(src, dest)
Definition: eostools.py:511
std
Definition: JetResolutionObject.h:76
reco::PhysObjectMatcher::IndexPair
std::pair< size_t, size_t > IndexPair
Definition: PhysObjectMatcher.h:67
reco::PhysObjectMatcher::config_
edm::ParameterSet config_
Definition: PhysObjectMatcher.h:70
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
volumeBasedMagneticField_160812_cfi.master
master
Definition: volumeBasedMagneticField_160812_cfi.py:60
reco::PhysObjectMatcher::resolveByMatchQuality_
bool resolveByMatchQuality_
Definition: PhysObjectMatcher.h:75
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
Generated for CMSSW Reference Manual by doxygen 1.8.16