ATLAS Offline Software
Loading...
Searching...
No Matches
VertexNewMatcher Class Referenceabstract

#include <VertexNewMatcher.h>

Inheritance diagram for VertexNewMatcher:
Collaboration diagram for VertexNewMatcher:

Public Types

typedef std::map< TIDA::Vertex *, TIDA::Vertex * > map_type
typedef std::map< TIDA::Vertex *, TIDA::Vertex * > rmap_type

Public Member Functions

 VertexNewMatcher (const std::string &s, double d=0.5)
 ~VertexNewMatcher ()
BestMatcher< TIDA::Vertex > * clone () override
virtual double distance (const TIDA::Vertex *v0, const TIDA::Vertex *v1) const override
virtual void match (const std::vector< TIDA::Vertex * > &ref, const std::vector< TIDA::Vertex * > &test)
virtual void match (const std::vector< TIDA::Vertex * > &s1, const std::vector< TIDA::Vertex * > &s2)=0
virtual const TIDA::Vertexmatched (TIDA::Vertex *t)
virtual const TIDA::Vertexrevmatched (TIDA::Vertex *t)
void clear ()
unsigned size () const

Protected Member Functions

std::map< int, int > matcher (const std::vector< Tp * > &ref, const std::vector< Tq * > &test)

Protected Attributes

double m_d
std::string mname
map_type mmatched
rmap_type mrevmatched

Detailed Description

Definition at line 22 of file VertexNewMatcher.h.

Member Typedef Documentation

◆ map_type

typedef std::map<TIDA::Vertex*, TIDA::Vertex*> TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::map_type
inherited

Definition at line 28 of file TIDAAssociator.h.

◆ rmap_type

typedef std::map<TIDA::Vertex*, TIDA::Vertex*> TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::rmap_type
inherited

Definition at line 29 of file TIDAAssociator.h.

Constructor & Destructor Documentation

◆ VertexNewMatcher()

VertexNewMatcher::VertexNewMatcher ( const std::string & s,
double d = 0.5 )
inline

Definition at line 26 of file VertexNewMatcher.h.

26: BestMatcher<TIDA::Vertex>( s, d*d ) { }
BestMatcher(const std::string &name, double d)
Definition BestMatcher.h:61

◆ ~VertexNewMatcher()

VertexNewMatcher::~VertexNewMatcher ( )
inline

Definition at line 28 of file VertexNewMatcher.h.

28{ }

Member Function Documentation

◆ clear()

void TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::clear ( )
inlineinherited

Definition at line 63 of file TIDAAssociator.h.

◆ clone()

BestMatcher< TIDA::Vertex > * VertexNewMatcher::clone ( )
inlineoverridevirtual

Implements TIDA::Associator< TIDA::Vertex, TIDA::Vertex >.

Definition at line 30 of file VertexNewMatcher.h.

30{ return new VertexNewMatcher(*this); }
VertexNewMatcher(const std::string &s, double d=0.5)

◆ distance()

virtual double VertexNewMatcher::distance ( const TIDA::Vertex * v0,
const TIDA::Vertex * v1 ) const
inlineoverridevirtual

Implements BestMatcher< TIDA::Vertex >.

Definition at line 32 of file VertexNewMatcher.h.

32 {
33
34 // vectors of pointers to tracks belonging to ref and test vertices
35 const std::vector<TIDA::Track*>& refTracks= v0->tracks();
36 const std::vector<TIDA::Track*>& testTracks = v1->tracks();
37
38 // match offline to trigger tracks using delta R
39 Associator_DeltaRMatcher track_matcher( "track_matcher", 0.05 );
40 track_matcher.match( refTracks, testTracks );
41
42 double n_matches = (double)track_matcher.size();
43 double eff = n_matches/refTracks.size();
44
45 // 1-eff as BestMatcher uses lower than threshold, squared as BestMatcher takes d*d
46 return (1-eff)*(1-eff);
47 };
const std::vector< TIDA::Track * > & tracks() const
Definition TIDAVertex.h:67

◆ match() [1/2]

virtual void BestMatcher< TIDA::Vertex, TIDA::Vertex >::match ( const std::vector< TIDA::Vertex * > & ref,
const std::vector< TIDA::Vertex * > & test )
inlinevirtualinherited

Implements TIDA::Associator< TIDA::Vertex, TIDA::Vertex >.

Definition at line 68 of file BestMatcher.h.

69 {
70 this->clear();
71
73
75 while ( mitr!=matched.end() ) {
76 this->mmatched.insert( typename TIDA::Associator<T,S>::map_type::value_type( ref[mitr->first], test[mitr->second]) );
77 this->mrevmatched.insert( typename TIDA::Associator<S,T>::map_type::value_type( test[mitr->second], ref[mitr->first] ) );
78 ++mitr;
79 }
80
81 }
std::map< int, int > matcher(const std::vector< Tp * > &ref, const std::vector< Tq * > &test)
Definition BestMatcher.h:91
virtual const TIDA::Vertex * matched(TIDA::Vertex *t)

◆ match() [2/2]

virtual void TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::match ( const std::vector< TIDA::Vertex * > & s1,
const std::vector< TIDA::Vertex * > & s2 )
pure virtualinherited

◆ matched()

virtual const TIDA::Vertex * TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::matched ( TIDA::Vertex * t)
inlinevirtualinherited

Definition at line 45 of file TIDAAssociator.h.

45 {
46 typename map_type::const_iterator titr = mmatched.find(t);
47 if ( titr != mmatched.end() ) return titr->second;
48 else return 0;
49 }

◆ matcher()

std::map< int, int > BestMatcher< TIDA::Vertex, TIDA::Vertex >::matcher ( const std::vector< Tp * > & ref,
const std::vector< Tq * > & test )
inlineprotectedinherited

inserting into a multimatch sorts them into order as they are added so you get the list of best matches imediately NB: probably quicker to add them to a vector and then sort them afterward

now go through from best to worst, adding to another sorted set only the matches that do not use tracks already used by a better matching pair already found

set of unique track pairings NB: for some reason, using the multiset is faster than using a vector, even though we are only iterating through it later

hooray!! now print out the answer, make a map etc this isn't technically needed, could just use the set of "matched" objects, but hey ho

Definition at line 91 of file BestMatcher.h.

91 {
92
98
100
101 for (unsigned int i=0 ; i<ref.size() ; i++ ) {
102
103 for (unsigned int j=0 ; j<test.size() ; j++ ) {
104 double d = distance(ref[i], test[j]);
105 if ( d<m_d ){
106 m.insert( matched_(d, i, j) );
107 }
108 }
109 }
110
114
115 std::vector<bool> refused( ref.size(), false );
116 std::vector<bool> testused( test.size(), false );
117
123
124 typename std::multiset<matched_>::iterator mitr = m.begin();
125
126 double chi2 = 0;
127
128 for ( ; mitr!=m.end() ; ++mitr ) {
129
130 int rind = mitr->first();
131 int tind = mitr->second();
132
133 if ( refused[rind] ) continue;
134 if ( testused[tind] ) continue;
135
136 refused[rind] = true;
137 testused[tind] = true;
138
139 unique.insert( *mitr );
140
141 chi2 += (mitr->d()*mitr->d());
142
143 }
144
145
146 // std::cout << "chi2 of matches " << chi2 << std::endl;
147
151
152 // std::cout << "\nmatched" << std::endl;
153
155
156
157 mitr = unique.begin();
158 while ( mitr!=unique.end() ) {
159 matches.insert( std::map<int, int>::value_type( mitr->first(), mitr->second() ) );
160 // std::cout << "\tbest match " << *mitr << "\t" << ref[mitr->first()] << "\t" << test[mitr->second()] << std::endl;
161 ++mitr;
162 }
163
164 return matches;
165
166 }
virtual double distance(const TIDA::Vertex *t0, const TIDA::Vertex *t1) const=0
unsigned size() const

◆ revmatched()

virtual const TIDA::Vertex * TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::revmatched ( TIDA::Vertex * t)
inlinevirtualinherited

Definition at line 52 of file TIDAAssociator.h.

52 {
53 typename rmap_type::const_iterator titr = mrevmatched.find(t);
54 if ( titr != mrevmatched.end() ) return titr->second;
55 else return 0;
56 }

◆ size()

unsigned TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::size ( ) const
inlineinherited

Definition at line 66 of file TIDAAssociator.h.

66{ return mmatched.size(); }

Member Data Documentation

◆ m_d

double BestMatcher< TIDA::Vertex, TIDA::Vertex >::m_d
protectedinherited

Definition at line 171 of file BestMatcher.h.

◆ mmatched

map_type TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::mmatched
protectedinherited

Definition at line 73 of file TIDAAssociator.h.

◆ mname

std::string TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::mname
protectedinherited

Definition at line 71 of file TIDAAssociator.h.

◆ mrevmatched

rmap_type TIDA::Associator< TIDA::Vertex, TIDA::Vertex >::mrevmatched
protectedinherited

Definition at line 74 of file TIDAAssociator.h.


The documentation for this class was generated from the following file: