Trk::MultiComponentStateHelpers Namespace Reference

Functions
MultiComponentState	clone (const MultiComponentState &in)
	Clone TrackParameters method.
MultiComponentState	WithScaledError (MultiComponentState &&in, double errorScaleLocX, double errorScaleLocY, double errorScalePhi, double errorScaleTheta, double errorScaleQoverP)
	Scale the covariance matrix components by individual factors.
bool	allHaveCovariance (const MultiComponentState &in)
	Check to see if all components have covariance Matrix.
void	renormaliseState (MultiComponentState &, double norm=1)
	Performing renormalisation of total state weighting to one.

Function Documentation

allHaveCovariance()

bool Trk::MultiComponentStateHelpers::allHaveCovariance

(

const MultiComponentState &

in

)

Check to see if all components have covariance Matrix.

Definition at line 66 of file ComponentParameters.cxx.

{
  bool allHaveCovariance = true;
  for (const ComponentParameters& component : in) {
    const AmgSymMatrix(5)* originalMatrix = component.params->covariance();
    if (!originalMatrix) {
      allHaveCovariance = false;
      break;
    }
  }
  return allHaveCovariance;
}

clone()

Trk::MultiComponentState Trk::MultiComponentStateHelpers::clone

(

const MultiComponentState &

in

)

Clone TrackParameters method.

Definition at line 15 of file ComponentParameters.cxx.

{
  auto clonedState = Trk::MultiComponentState();
  clonedState.reserve(in.size());
  for (const ComponentParameters& component : in) {
    clonedState.push_back({component.params->uniqueClone(), component.weight});
  }
  return clonedState;
}

renormaliseState()

void Trk::MultiComponentStateHelpers::renormaliseState	(	Trk::MultiComponentState &	in,
		double	norm = 1 )

Performing renormalisation of total state weighting to one.

Definition at line 80 of file ComponentParameters.cxx.

{
  // Determine total weighting of state
  double sumWeights = 0.;
  for (const ComponentParameters& component : in) {
    sumWeights += component.weight;
  }
  if (sumWeights == 0) {
    return;
  }
  double normalise = norm / sumWeights;
  for (ComponentParameters& component : in) {
    component.weight = component.weight * normalise;
  }
}

WithScaledError()

Trk::MultiComponentState Trk::MultiComponentStateHelpers::WithScaledError	(	Trk::MultiComponentState &&	in,
		double	errorScaleLocX,
		double	errorScaleLocY,
		double	errorScalePhi,
		double	errorScaleTheta,
		double	errorScaleQoverP )

Scale the covariance matrix components by individual factors.

Definition at line 26 of file ComponentParameters.cxx.

{
  AmgSymMatrix(5) coefficients;
  coefficients(0, 0) = (errorScaleLocX * errorScaleLocX);
  coefficients(1, 1) = (errorScaleLocY * errorScaleLocY);
  coefficients(2, 2) = (errorScalePhi * errorScalePhi);
  coefficients(3, 3) = (errorScaleTheta * errorScaleTheta);
  coefficients(4, 4) = (errorScaleQoverP * errorScaleQoverP);
  coefficients.fillSymmetric(0, 1, (errorScaleLocX * errorScaleLocY));
  coefficients.fillSymmetric(0, 2, (errorScaleLocX * errorScalePhi));
  coefficients.fillSymmetric(0, 3, (errorScaleLocX * errorScaleTheta));
  coefficients.fillSymmetric(0, 4, (errorScaleLocX * errorScaleQoverP));
  coefficients.fillSymmetric(1, 2, (errorScaleLocY * errorScalePhi));
  coefficients.fillSymmetric(1, 3, (errorScaleLocY * errorScaleTheta));
  coefficients.fillSymmetric(1, 4, (errorScaleLocY * errorScaleQoverP));
  coefficients.fillSymmetric(2, 3, (errorScalePhi * errorScaleTheta));
  coefficients.fillSymmetric(2, 4, (errorScalePhi * errorScaleQoverP));
  coefficients.fillSymmetric(3, 4, (errorScaleTheta * errorScaleQoverP));
 
  for (ComponentParameters& component : in) {
    Trk::TrackParameters* trackParameters = component.params.get();
    const AmgSymMatrix(5)* originalMatrix = trackParameters->covariance();
    // If no covariance skip
    if (!originalMatrix) {
      continue;
    }
    AmgSymMatrix(5) newCovarianceMatrix =
      (*originalMatrix).cwiseProduct(coefficients);
    trackParameters->updateParameters(trackParameters->parameters(),
                                      newCovarianceMatrix);
  }
  return { std::move(in) };
}