lcsim/src/org/lcsim/contrib/seedtracker
diff -u -r1.1 -r1.2
--- ConstrainHelix.java 4 Apr 2008 23:41:33 -0000 1.1
+++ ConstrainHelix.java 9 May 2008 01:14:35 -0000 1.2
@@ -35,21 +35,21 @@
double curvmax = 0.0003 * _bfield / strategy.getMinPT();
double curv = Math.abs(params[2]);
if (curv > curvmax) {
- nhchisq += Math.pow(curv - curvmax, 2) / helix.covariance().e(2, 2);
+ nhchisq += Math.pow(curv - curvmax, 2) / Math.abs(helix.covariance().e(2, 2));
}
// Inflate chisq if |DCA| is too large
double dcamax = strategy.getMaxDCA();
double dca = Math.abs(params[0]);
if (dca > dcamax) {
- nhchisq += Math.pow(dca - dcamax, 2) / helix.covariance().e(0, 0);
+ nhchisq += Math.pow(dca - dcamax, 2) / Math.abs(helix.covariance().e(0, 0));
}
// Inflate chisq if |z0| is too large
double z0max = strategy.getMaxZ0();
double z0 = Math.abs(params[3]);
if (z0 > z0max) {
- nhchisq += Math.pow(z0 - z0max, 2) / helix.covariance().e(3, 3);
+ nhchisq += Math.pow(z0 - z0max, 2) / Math.abs(helix.covariance().e(3, 3));
}
// Set the non-holenomic chi squared term in the helix