Skip to content

Commit

Permalink
fix: Assertion failure in proapgation example
Browse files Browse the repository at this point in the history
This was for some reason using matrix access with two arguments to set
elements of a vector (the diagonal of a matrix). Switched to use single
commponent element access.
  • Loading branch information
paulgessinger committed Mar 2, 2022
1 parent 3ee4e2e commit ea1b7ca
Showing 1 changed file with 6 additions and 12 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -134,23 +134,17 @@ inline ActsExamples::PropagationAlgorithm::Config readPropagationConfig(
/// Set the covariance transport to true
pAlgConfig.covarianceTransport = true;
/// Set the covariance matrix
pAlgConfig.covariances(Acts::BoundIndices::eBoundLoc0,
Acts::BoundIndices::eBoundLoc0) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundLoc0) =
pAlgConfig.d0Sigma * pAlgConfig.d0Sigma;
pAlgConfig.covariances(Acts::BoundIndices::eBoundLoc1,
Acts::BoundIndices::eBoundLoc1) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundLoc1) =
pAlgConfig.z0Sigma * pAlgConfig.z0Sigma;
pAlgConfig.covariances(Acts::BoundIndices::eBoundPhi,
Acts::BoundIndices::eBoundPhi) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundPhi) =
pAlgConfig.phiSigma * pAlgConfig.phiSigma;
pAlgConfig.covariances(Acts::BoundIndices::eBoundTheta,
Acts::BoundIndices::eBoundTheta) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundTheta) =
pAlgConfig.thetaSigma * pAlgConfig.thetaSigma;
pAlgConfig.covariances(Acts::BoundIndices::eBoundQOverP,
Acts::BoundIndices::eBoundQOverP) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundQOverP) =
pAlgConfig.qpSigma * pAlgConfig.qpSigma;
pAlgConfig.covariances(Acts::BoundIndices::eBoundTime,
Acts::BoundIndices::eBoundTime) =
pAlgConfig.covariances(Acts::BoundIndices::eBoundTime) =
pAlgConfig.tSigma * pAlgConfig.tSigma;

// Only if they are properly defined, assign off-diagonals
Expand Down

0 comments on commit ea1b7ca

Please sign in to comment.