|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
|
#ifndef JetReco_JetTrackMatch_h
#define JetReco_JetTrackMatch_h
/** \class reco::JetTrackMatch
*
* \short Association between Jets from jet collection and tracks from track collection
*
* Every jet may have several tracks associated with it
*
* class definition:
*
* template <typename JetC>
* class JetTrackMatch {
* public:
* typedef edm::Ref<JetC> JetRef;
* typedef edm::Ref<reco::TrackCollection> TrackRef;
* JetTrackMatch ();
* ~JetTrackMatch ();
* // insert orphan jet
* void insert (const JetRef& fJet);
* // assign track to jet.
* void insert (const JetRef& fJet, const TrackRef& fTrack);
* // get list of all tracks in the map
* std::vector <JetRef> allJets () const;
* // get all tracks associated with jet
* std::vector <TrackRef> getTracks (const JetRef& mJet) const;
* };
*
* \author Fedor Ratnikov, UMd
*
************************************************************/
#include "DataFormats/Common/interface/AssociationMap.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
namespace reco {
template <typename JetC>
class JetTrackMatch {
public:
typedef edm::Ref<JetC> JetRef;
typedef edm::Ref<reco::TrackCollection> TrackRef;
typedef edm::AssociationMap<edm::OneToMany<JetC, reco::TrackCollection> > Map;
private:
Map mMap;
public:
JetTrackMatch() {}
~JetTrackMatch() {}
/// insert orphan jet
void insert(const JetRef& fJet) { mMap.insert(fJet, TrackRef()); }
/// assign track to jet.
void insert(const JetRef& fJet, const TrackRef& fTrack) { mMap.insert(fJet, fTrack); }
/// get list of all jats in the map
std::vector<JetRef> allJets() const {
std::vector<JetRef> result;
typename Map::const_iterator it = mMap.begin();
for (; it != mMap.end(); ++it) {
result.push_back(it->key);
}
return result;
}
/// get all tracks associated with jet
std::vector<TrackRef> getTracks(const JetRef& mJet) const {
std::vector<TrackRef> result;
reco::TrackRefVector tracks = mMap[mJet];
int i = tracks.size();
while (--i >= 0) {
if (!tracks[i].isNull())
result.push_back(tracks[i]);
}
return result;
}
};
} // namespace reco
#endif
|