Skip to content

Commit 5f9b6ad

Browse files
rovereJanGerritSchulz
authored andcommitted
Fix Jacobian for linefit in Brokenline.
Co-authored-by: Jan Schulz (LXPLUS) <[email protected]>
1 parent 97e772d commit 5f9b6ad

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

RecoTracker/PixelTrackFitting/interface/alpaka/BrokenLine.h

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -495,8 +495,8 @@ namespace ALPAKA_ACCELERATOR_NAMESPACE::brokenline {
495495
vMat(2, 1) = vMat(1, 2) = hits_ge.col(i)[4]; // cov_yz
496496
vMat(2, 2) = hits_ge.col(i)[5]; // z errors
497497
auto tmp = 1. / radii.block(0, i, 2, 1).norm();
498-
jacobXYZtosZ(0, 0) = radii(1, i) * tmp;
499-
jacobXYZtosZ(0, 1) = -radii(0, i) * tmp;
498+
jacobXYZtosZ(0, 0) = data.qCharge * radii(1, i) * tmp;
499+
jacobXYZtosZ(0, 1) = -data.qCharge * radii(0, i) * tmp;
500500
jacobXYZtosZ(1, 2) = 1.;
501501
weights(i) = 1. / ((rotMat * jacobXYZtosZ * vMat * jacobXYZtosZ.transpose() * rotMat.transpose())(
502502
1, 1)); // compute the orthogonal weight point by point

0 commit comments

Comments
 (0)