lcsim/src/org/lcsim/fit/helicaltrack

diff -u -r1.9 -r1.10
--- HelixUtils.java	17 Dec 2009 20:32:29 -0000	1.9
+++ HelixUtils.java	12 May 2010 21:19:59 -0000	1.10
@@ -256,15 +256,15 @@

         double dca = helix.dca();
         double omega = helix.curvature();
         //  Calculate the non-zero derivatives of the direction with respect to the helix parameters

-        deriv.setElement(0, helix.curvatureIndex, (u.x() - cphi0 * sth) / omega);  // du_x / domega
-        deriv.setElement(1, helix.curvatureIndex, (u.y() - sphi0 * sth) / omega);  // du_y / domega
-        deriv.setElement(0, helix.dcaIndex, -omega * cphi0 * sth);  // du_x / ddca
-        deriv.setElement(1, helix.dcaIndex, -omega * sphi0 * sth);  // du_y / ddca
-        deriv.setElement(0, helix.phi0Index, -(1 - dca * omega) * sphi0 * sth);  // du_x / dphi0
-        deriv.setElement(1, helix.phi0Index,  (1 - dca * omega) * cphi0 * sth);  // du_y / dphi0
-        deriv.setElement(0, helix.slopeIndex, -u.x() * sth * cth);  // du_x / dslope
-        deriv.setElement(1, helix.slopeIndex, -u.y() * sth * cth);  // du_y / dslope
-        deriv.setElement(2, helix.slopeIndex, sth*sth*sth);  // du_z / dslope

+        deriv.setElement(0, HelicalTrackFit.curvatureIndex, (u.x() - cphi0 * sth) / omega);  // du_x / domega
+        deriv.setElement(1, HelicalTrackFit.curvatureIndex, (u.y() - sphi0 * sth) / omega);  // du_y / domega
+        deriv.setElement(0, HelicalTrackFit.dcaIndex, -omega * cphi0 * sth);  // du_x / ddca
+        deriv.setElement(1, HelicalTrackFit.dcaIndex, -omega * sphi0 * sth);  // du_y / ddca
+        deriv.setElement(0, HelicalTrackFit.phi0Index, -(1 - dca * omega) * sphi0 * sth);  // du_x / dphi0
+        deriv.setElement(1, HelicalTrackFit.phi0Index,  (1 - dca * omega) * cphi0 * sth);  // du_y / dphi0
+        deriv.setElement(0, HelicalTrackFit.slopeIndex, -u.x() * sth * cth);  // du_x / dslope
+        deriv.setElement(1, HelicalTrackFit.slopeIndex, -u.y() * sth * cth);  // du_y / dslope
+        deriv.setElement(2, HelicalTrackFit.slopeIndex, sth*sth*sth);  // du_z / dslope

         //  Return the derivative matrix
         return deriv;
     }