CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
pat::MatcherUsingTracks Class Reference

Matcher of reconstructed objects to other reconstructed objects using the tracks inside them. More...

Inheritance diagram for pat::MatcherUsingTracks:
edm::stream::EDProducer<>

Public Member Functions

 MatcherUsingTracks (const edm::ParameterSet &iConfig)
 
void produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
 
 ~MatcherUsingTracks () override
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
 EDProducer (const EDProducer &)=delete
 
bool hasAbilityToProduceInBeginLumis () const final
 
bool hasAbilityToProduceInBeginProcessBlocks () const final
 
bool hasAbilityToProduceInBeginRuns () const final
 
bool hasAbilityToProduceInEndLumis () const final
 
bool hasAbilityToProduceInEndProcessBlocks () const final
 
bool hasAbilityToProduceInEndRuns () const final
 
const EDProduceroperator= (const EDProducer &)=delete
 

Private Member Functions

template<typename T >
void storeValueMap (edm::Event &iEvent, const edm::Handle< edm::View< reco::Candidate > > &handle, const std::vector< T > &values, const std::string &label) const
 Store extra information in a ValueMap. More...
 

Private Attributes

MatcherUsingTracksAlgorithm algo_
 The real workhorse. More...
 
bool dontFailOnMissingInput_
 Some extra configurables. More...
 
edm::EDGetTokenT< edm::View
< reco::Candidate > > 
matchedToken_
 
edm::EDGetTokenT< edm::View
< reco::Candidate > > 
srcToken_
 Labels for input collections. More...
 
bool writeExtraPATOutput_
 

Additional Inherited Members

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

Detailed Description

Matcher of reconstructed objects to other reconstructed objects using the tracks inside them.

Author
Giovanni Petrucciani
Version
Id:
MatcherUsingTracks.cc,v 1.5 2010/07/12 20:56:11 gpetrucc Exp

Definition at line 32 of file MatcherUsingTracks.cc.

Constructor & Destructor Documentation

pat::MatcherUsingTracks::MatcherUsingTracks ( const edm::ParameterSet iConfig)
explicit

Definition at line 61 of file MatcherUsingTracks.cc.

References algo_, and writeExtraPATOutput_.

62  : srcToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("src"))),
63  matchedToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("matched"))),
64  algo_(iConfig, consumesCollector()),
65  dontFailOnMissingInput_(iConfig.existsAs<bool>("dontFailOnMissingInput")
66  ? iConfig.getParameter<bool>("dontFailOnMissingInput")
67  : false),
69  iConfig.existsAs<bool>("writeExtraPATOutput") ? iConfig.getParameter<bool>("writeExtraPATOutput") : false) {
70  // this is the basic output (edm::Association is not generic)
71  produces<edm::ValueMap<reco::CandidatePtr> >();
73  // this is the crazy stuff to get the same with UserData
74  produces<edm::OwnVector<pat::UserData> >();
75  produces<edm::ValueMap<edm::Ptr<pat::UserData> > >();
76  }
77  // this is the extra stuff
78  if (algo_.hasMetrics()) {
79  produces<edm::ValueMap<float> >("deltaR");
80  produces<edm::ValueMap<float> >("deltaEta");
81  produces<edm::ValueMap<float> >("deltaPhi");
82  produces<edm::ValueMap<float> >("deltaLocalPos");
83  produces<edm::ValueMap<float> >("deltaPtRel");
84  if (algo_.hasChi2()) {
85  produces<edm::ValueMap<float> >("chi2");
86  }
87  } else {
88  produces<edm::ValueMap<int> >("matched");
89  }
90 }
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:171
edm::EDGetTokenT< edm::View< reco::Candidate > > matchedToken_
edm::EDGetTokenT< edm::View< reco::Candidate > > srcToken_
Labels for input collections.
MatcherUsingTracksAlgorithm algo_
The real workhorse.
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
bool dontFailOnMissingInput_
Some extra configurables.
pat::MatcherUsingTracks::~MatcherUsingTracks ( )
inlineoverride

Definition at line 35 of file MatcherUsingTracks.cc.

35 {}

Member Function Documentation

void pat::MatcherUsingTracks::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 92 of file MatcherUsingTracks.cc.

References algo_, HLT_FULL_cff::chi2, HLT_FULL_cff::deltaEtas, HLT_FULL_cff::deltaPhis, edm::Event::getByToken(), iEvent, pat::UserData::make(), match(), eostools::move(), edm::Event::put(), and alcazmumu_cfi::src.

92  {
93  using namespace edm;
94  using namespace std;
95 
96  algo_.init(iSetup);
97 
99 
100  iEvent.getByToken(srcToken_, src);
101  iEvent.getByToken(matchedToken_, matched);
102 
103  // declare loop variables and some intermediate stuff
105  int isrc, nsrc = src->size();
106 
107  // working and output variables
108  vector<int> match(nsrc, -1);
109  vector<float> deltaRs(nsrc, 999);
110  vector<float> deltaEtas(nsrc, 999);
111  vector<float> deltaPhis(nsrc, 999);
112  vector<float> deltaPtRel(nsrc, 999);
113  vector<float> deltaLocalPos(nsrc, 999);
114  vector<float> chi2(nsrc, 999999);
115 
116  // don't try matching if the input collection is missing and the module is configured to fail silently
117  if (!(matched.failedToGet() && dontFailOnMissingInput_)) {
118  // loop on the source collection, and request for the match
119  for (itsrc = src->begin(), edsrc = src->end(), isrc = 0; itsrc != edsrc; ++itsrc, ++isrc) {
120  match[isrc] = algo_.match(*itsrc,
121  *matched,
122  deltaRs[isrc],
123  deltaEtas[isrc],
124  deltaPhis[isrc],
125  deltaLocalPos[isrc],
126  deltaPtRel[isrc],
127  chi2[isrc]);
128  }
129  }
130 
131  std::vector<reco::CandidatePtr> ptrs(nsrc);
132  for (isrc = 0; isrc < nsrc; ++isrc) {
133  if (match[isrc] != -1) {
134  ptrs[isrc] = matched->ptrAt(match[isrc]);
135  }
136  }
137  storeValueMap<reco::CandidatePtr>(iEvent, src, ptrs, "");
138 
139  if (writeExtraPATOutput_) {
140  auto outUDVect = std::make_unique<edm::OwnVector<pat::UserData> >();
141  std::vector<int> idxUD(nsrc, -1);
142  for (isrc = 0; isrc < nsrc; ++isrc) {
143  if (match[isrc] != -1) {
144  outUDVect->push_back(pat::UserData::make(ptrs[isrc]));
145  idxUD[isrc] = outUDVect->size() - 1;
146  }
147  }
148  edm::OrphanHandle<edm::OwnVector<pat::UserData> > doneUDVect = iEvent.put(std::move(outUDVect));
149  std::vector<edm::Ptr<pat::UserData> > ptrUD(nsrc);
150  for (isrc = 0; isrc < nsrc; ++isrc) {
151  if (idxUD[isrc] != -1)
152  ptrUD[isrc] = edm::Ptr<pat::UserData>(doneUDVect, idxUD[isrc]);
153  }
154  storeValueMap<edm::Ptr<pat::UserData> >(iEvent, src, ptrUD, "");
155  }
156 
157  if (algo_.hasMetrics()) {
158  storeValueMap<float>(iEvent, src, deltaRs, "deltaR");
159  storeValueMap<float>(iEvent, src, deltaEtas, "deltaEta");
160  storeValueMap<float>(iEvent, src, deltaPhis, "deltaPhi");
161  storeValueMap<float>(iEvent, src, deltaLocalPos, "deltaLocalPos");
162  storeValueMap<float>(iEvent, src, deltaPtRel, "deltaPtRel");
163  if (algo_.hasChi2()) {
164  storeValueMap<float>(iEvent, src, chi2, "chi2");
165  }
166  } else {
167  std::vector<int> ismatched(nsrc, 0);
168  for (isrc = 0; isrc < nsrc; ++isrc) {
169  ismatched[isrc] = (match[isrc] != -1);
170  }
171  storeValueMap<int>(iEvent, src, ismatched, "matched");
172  }
173 }
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
edm::EDGetTokenT< edm::View< reco::Candidate > > matchedToken_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:539
edm::EDGetTokenT< edm::View< reco::Candidate > > srcToken_
Labels for input collections.
static std::unique_ptr< UserData > make(const T &value, bool transientOnly=false)
MatcherUsingTracksAlgorithm algo_
The real workhorse.
int iEvent
Definition: GenABIO.cc:224
def move
Definition: eostools.py:511
bool dontFailOnMissingInput_
Some extra configurables.
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:86
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
template<typename T >
void pat::MatcherUsingTracks::storeValueMap ( edm::Event iEvent,
const edm::Handle< edm::View< reco::Candidate > > &  handle,
const std::vector< T > &  values,
const std::string &  label 
) const
private

Store extra information in a ValueMap.

Definition at line 176 of file MatcherUsingTracks.cc.

References edm::helper::Filler< Map >::fill(), patZpeak::handle, edm::helper::Filler< Map >::insert(), eostools::move(), and edm::Event::put().

179  {
180  using namespace edm;
181  using namespace std;
182  unique_ptr<ValueMap<T> > valMap(new ValueMap<T>());
183  typename edm::ValueMap<T>::Filler filler(*valMap);
184  filler.insert(handle, values.begin(), values.end());
185  filler.fill();
186  iEvent.put(std::move(valMap), label);
187 }
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:133
char const * label
def move
Definition: eostools.py:511

Member Data Documentation

MatcherUsingTracksAlgorithm pat::MatcherUsingTracks::algo_
private

The real workhorse.

Definition at line 45 of file MatcherUsingTracks.cc.

Referenced by MatcherUsingTracks().

bool pat::MatcherUsingTracks::dontFailOnMissingInput_
private

Some extra configurables.

Definition at line 48 of file MatcherUsingTracks.cc.

edm::EDGetTokenT<edm::View<reco::Candidate> > pat::MatcherUsingTracks::matchedToken_
private

Definition at line 42 of file MatcherUsingTracks.cc.

edm::EDGetTokenT<edm::View<reco::Candidate> > pat::MatcherUsingTracks::srcToken_
private

Labels for input collections.

Definition at line 41 of file MatcherUsingTracks.cc.

bool pat::MatcherUsingTracks::writeExtraPATOutput_
private

Definition at line 49 of file MatcherUsingTracks.cc.

Referenced by MatcherUsingTracks().