Skip to content

Commit

Permalink
step 1
Browse files Browse the repository at this point in the history
  • Loading branch information
AJPfleger committed Oct 31, 2024
1 parent 21aae58 commit ff1ef92
Showing 1 changed file with 23 additions and 21 deletions.
44 changes: 23 additions & 21 deletions Core/include/Acts/TrackFitting/GlobalChiSquareFitter.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -372,30 +372,30 @@ void addMeasurementToGx2fSums(Eigen::MatrixXd& aMatrixExtended,

ACTS_VERBOSE(
"Contributions in addMeasurementToGx2fSums:\n"
<< "kMeasDim: " << kMeasDim << "\n"
<< "predicted" << predicted.transpose() << "\n"
<< "measurement: " << measurement.transpose() << "\n"
<< "covarianceMeasurement:\n"
<< " kMeasDim: " << kMeasDim << "\n"
<< " predicted" << predicted.transpose() << "\n"
<< " measurement: " << measurement.transpose() << "\n"
<< " covarianceMeasurement:\n"
<< covarianceMeasurement << "\n"
<< "projector:\n"
<< " projector:\n"
<< projector.eval() << "\n"
<< "projJacobian:\n"
<< " projJacobian:\n"
<< projJacobian.eval() << "\n"
<< "projPredicted: " << (projPredicted.transpose()).eval() << "\n"
<< "residual: " << (residual.transpose()).eval() << "\n"
<< "extendedJacobian:\n"
<< " projPredicted: " << (projPredicted.transpose()).eval() << "\n"
<< " residual: " << (residual.transpose()).eval() << "\n"
<< " extendedJacobian:\n"
<< extendedJacobian << "\n"
<< "aMatrixMeas:\n"
<< " aMatrix contribution:\n"
<< (projJacobian.transpose() * (*safeInvCovMeasurement) * projJacobian)
.eval()
<< "\n"
<< "bVectorMeas: "
<< " bVector contribution: "
<< (residual.transpose() * (*safeInvCovMeasurement) * projJacobian).eval()
<< "\n"
<< "chi2sumMeas: "
<< " chi2sumMeas: "
<< (residual.transpose() * (*safeInvCovMeasurement) * residual)(0, 0)
<< "\n"
<< "safeInvCovMeasurement:\n"
<< " safeInvCovMeasurement:\n"
<< (*safeInvCovMeasurement));

return;
Expand Down Expand Up @@ -727,13 +727,14 @@ class Gx2Fitter {
if (doMaterial) {
ACTS_DEBUG(" Update parameters with scattering angles.");
const auto scatteringMapId = scatteringMap->find(geoId);
ACTS_VERBOSE(" scatteringAngles:\n"
ACTS_VERBOSE(" scatteringAngles:\n"
<< scatteringMapId->second.scatteringAngles()
<< "\n boundParams before the update:\n"
<< boundParams);
<< "\n boundParams before the update: "
<< boundParams.parameters().transpose());
boundParams.parameters() +=
scatteringMapId->second.scatteringAngles();
ACTS_VERBOSE(" boundParams after the update:\n" << boundParams);
ACTS_VERBOSE(" boundParams after the update: "
<< boundParams.parameters().transpose());
}

// Fill the track state
Expand Down Expand Up @@ -829,13 +830,14 @@ class Gx2Fitter {
// multipleScattering and have material
ACTS_DEBUG(" Update parameters with scattering angles.");
const auto scatteringMapId = scatteringMap->find(geoId);
ACTS_VERBOSE(" scatteringAngles:\n"
ACTS_VERBOSE(" scatteringAngles:\n"
<< scatteringMapId->second.scatteringAngles()
<< "\n boundParams before the update:\n"
<< boundParams);
<< "\n boundParams before the update: "
<< boundParams.parameters().transpose());
boundParams.parameters() +=
scatteringMapId->second.scatteringAngles();
ACTS_VERBOSE(" boundParams after the update:\n" << boundParams);
ACTS_VERBOSE(" boundParams after the update: "
<< boundParams.parameters().transpose());

// Fill the track state
trackStateProxy.smoothed() = boundParams.parameters();
Expand Down

0 comments on commit ff1ef92

Please sign in to comment.