lcsim/src/org/lcsim/recon/pfa/identifier
diff -u -r1.6 -r1.7
--- TrackHelixPlusHitExtrapolator.java 1 Mar 2011 16:26:47 -0000 1.6
+++ TrackHelixPlusHitExtrapolator.java 1 Mar 2011 17:53:28 -0000 1.7
@@ -53,7 +53,6 @@
double alpha = Double.NaN;
Hep3Vector pointOfClosestApproachToTrackHit = null;
Hep3Vector offset = null;
-
Hep3Vector positionOfOutermostHit = new BasicHep3Vector(outermostHit.getPosition());
double rr = Math.sqrt(positionOfOutermostHit.x()*positionOfOutermostHit.x()+positionOfOutermostHit.y()*positionOfOutermostHit.y());
double z = positionOfOutermostHit.z();
@@ -75,8 +74,9 @@
{
offset = new BasicHep3Vector(0.,0.,0.);
pointOfClosestApproachToTrackHit = positionOfOutermostHit;
+ alpha = m_swimmer.getTrackLengthToPoint(positionOfOutermostHit);
}
- /*
+/*
if (outermostHit instanceof org.lcsim.event.base.BaseTrackerHitMC) {
// This cheated hit has exact 3D info
Hep3Vector positionOfOutermostHit = new BasicHep3Vector(outermostHit.getPosition());
@@ -117,8 +117,8 @@
double z = positionOfOutermostHit.z();
throw new AssertionError("ERROR: Unknown hit of type "+outermostHit.getClass().getName()+" at r="+r+", z="+z);
}
- *
- */
+ *
+ */
Hep3Vector momentumAtPOCA = m_swimmer.getMomentumAtLength(alpha);
Hep3Vector newPoint = VecOp.add(pointOfClosestApproachToTrackHit, offset);