Skip to content
Snippets Groups Projects
Commit 4fbb7075 authored by Dave Casper's avatar Dave Casper
Browse files

Revert "Merge branch 'master' of https://gitlab.cern.ch:8443/faser/calypso"

This reverts commit 4c9e84aa
parent d89ba2c7
No related branches found
No related tags found
No related merge requests found
...@@ -112,30 +112,26 @@ Trk::SpacePoint* TrackerSpacePointMakerTool::makeSCT_SpacePoint(const Tracker::T ...@@ -112,30 +112,26 @@ Trk::SpacePoint* TrackerSpacePointMakerTool::makeSCT_SpacePoint(const Tracker::T
Amg::Vector3D s(a+b); Amg::Vector3D s(a+b);
Amg::Vector3D t(c+d); Amg::Vector3D t(c+d);
double lambda0 = r.cross(Amg::Vector3D(s-t)).z()/det; double lambda0 = r.cross(Amg::Vector3D(s-t)).z()/det;
double lambda1 = q.cross(Amg::Vector3D(s-t)).z()/det;
point = (s + lambda0 * q)/2; point = (s + lambda0 * q)/2;
ATH_MSG_VERBOSE( "Endpoints 1 : ( " << a.x() << " , " << a.y() << " , " << a.z() << " ) to (" << b.x() << " , " << b.y() << " , " << b.z() << " ) " ); ATH_MSG_VERBOSE( "Endpoints 1 : ( " << a.x() << " , " << a.y() << " , " << a.z() << " ) to (" << b.x() << " , " << b.y() << " , " << b.z() << " ) " );
ATH_MSG_VERBOSE( "Endpoints 2 : ( " << c.x() << " , " << c.y() << " , " << c.z() << " ) to (" << d.x() << " , " << d.y() << " , " << d.z() << " ) " ); ATH_MSG_VERBOSE( "Endpoints 2 : ( " << c.x() << " , " << c.y() << " , " << c.z() << " ) to (" << d.x() << " , " << d.y() << " , " << d.z() << " ) " );
ATH_MSG_VERBOSE( "Intersection: ( " << point.x() << " , " << point.y() << " , " << point.z() << " ) " ); ATH_MSG_VERBOSE( "Intersection: ( " << point.x() << " , " << point.y() << " , " << point.z() << " ) " );
// require that the point is within the bounds of at least one of the two strips
if (fabs(lambda0) > 1 + m_stripLengthTolerance)
{
if (fabs(lambda1) > 1 + m_stripLengthTolerance)
{
ATH_MSG_WARNING("Intersection lies outside the bounds of both strips");
ok = false;
}
}
} }
else else
{ {
ATH_MSG_WARNING("Alt Intersection failed"); ATH_MSG_WARNING("Intersection failed");
ok = false; ok = false;
} }
// Amg::Vector3D mab(c - a);
// double eaTeb = q.dot(r);
// double denom = 1. - eaTeb*eaTeb;
// if (fabs(denom)>10e-7){
// double lambda0 = (mab.dot(q) - mab.dot(r)*eaTeb)/denom;
// point = a+lambda0*q;
// ATH_MSG_VERBOSE( "Endpoints 1 : ( " << a.x() << " , " << a.y() << " , " << a.z() << " ) to (" << b.x() << " , " << b.y() << " , " << b.z() << " ) " );
// ATH_MSG_VERBOSE( "Endpoints 2 : ( " << c.x() << " , " << c.y() << " , " << c.z() << " ) to (" << d.x() << " , " << d.y() << " , " << d.z() << " ) " );
// ATH_MSG_VERBOSE( "Intersection: ( " << point.x() << " , " << point.y() << " , " << point.z() << " ) " );
// } else {
// ATH_MSG_WARNING("Intersection failed");
// ok = false;
// }
} else { } else {
Amg::Vector3D s(a+b-2.*vertexVec); // twice the vector from vertex to midpoint Amg::Vector3D s(a+b-2.*vertexVec); // twice the vector from vertex to midpoint
Amg::Vector3D t(c+d-2.*vertexVec); // twice the vector from vertex to midpoint Amg::Vector3D t(c+d-2.*vertexVec); // twice the vector from vertex to midpoint
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment