55 const std::vector<TIDA::Track*>& ,
59 std::vector<TIDA::Track*>::const_iterator
reference, referenceEnd=referenceTracks.end();
63 double referenceEta = (*reference)->eta();
64 double referencePhi =
phi((*reference)->phi());
65 double referenceZ0 = (*reference)->z0();
66 double referenceD0 = (*reference)->a0();
67 double referencePT = (*reference)->pT();
76 double testEta = -
test->eta();
78 double testZ0 =
test->z0();
79 double testD0 = -
test->a0();
80 double testPT = -
test->pT();
83 if(referencePT==0 || testPT==0)
continue;
86 double eeta = sqrt(
test->deta()*
test->deta() + (*reference)->deta()*(*reference)->deta() );
87 double ephi = sqrt(
test->dphi()*
test->dphi() + (*reference)->dphi()*(*reference)->dphi() );
88 double ez0 = sqrt(
test->dz0() *
test->dz0() + (*reference)->dz0() *(*reference)->dz0() );
89 double ed0 = sqrt(
test->da0() *
test->da0() + (*reference)->da0() *(*reference)->da0() );
90 double einvpT = sqrt(
test->dpT() *
test->dpT() + (*reference)->dpT() *(*reference)->dpT() );