Trk::iPatGlobalFitter Class Reference

GlobalTrackFitter tool providing methods used during alignment. More...

#include <iPatGlobalFitter.h>

Inheritance diagram for Trk::iPatGlobalFitter:

Collaboration diagram for Trk::iPatGlobalFitter:

Public Member Functions
Track *	alignmentFit (const EventContext &ctx, AlignmentCache &, const Track &, const RunOutlierRemoval runOutlier=false, const ParticleHypothesis matEffects=Trk::nonInteracting) const

Private Member Functions
std::unique_ptr< Amg::MatrixX >	derivMatrix (const FitState &fitState) const
	GlobalTrackFitter methods: access to the matrix of derivatives used during the latest track fit.
std::unique_ptr< Amg::MatrixX >	fullCovarianceMatrix (const FitState &fitState) const
	access to the global fitter's full covariance matrix
void	setMinIterations (int minIterations)
	set method for the minimum number of iterations for (alignment) friend

Static Private Member Functions
static int	iterationsOfLastFit (const FitState &fitState)
	access to the number of iterations taken by the latest track fit

Private Attributes
Gaudi::Property< bool >	m_allParameters

Detailed Description

GlobalTrackFitter tool providing methods used during alignment.

Definition at line 26 of file iPatGlobalFitter.h.

Member Function Documentation

alignmentFit()

Track * Trk::iPatGlobalFitter::alignmentFit	(	const EventContext &	ctx,
		AlignmentCache &	alignCache,
		const Track &	trk,
		const RunOutlierRemoval	runOutlier = false,
		const ParticleHypothesis	matEffects = Trk::nonInteracting ) const

Definition at line 22 of file iPatGlobalFitter.cxx.

                                               {
  //  @TODO ensure the number of iterations is passed through to the fitter
  // setMinIterations (alignCache.m_minIterations);
  alignCache.m_derivMatrix.reset();
  alignCache.m_fullCovarianceMatrix.reset();
  alignCache.m_iterationsOfLastFit = 0;
 
  auto [refittedTrack, fitState] =
      fitWithState(ctx, trk, runOutlier, matEffects);
 
  if (refittedTrack) {
    alignCache.m_derivMatrix = derivMatrix(*fitState);
    alignCache.m_fullCovarianceMatrix = fullCovarianceMatrix(*fitState);
    alignCache.m_iterationsOfLastFit = iterationsOfLastFit(*fitState);
  }
  return refittedTrack.release();
}

derivMatrix()

std::unique_ptr< Amg::MatrixX > Trk::iPatGlobalFitter::derivMatrix ( const FitState & fitState ) const

private

GlobalTrackFitter methods: access to the matrix of derivatives used during the latest track fit.

Definition at line 44 of file iPatGlobalFitter.cxx.

                                    {
  // copy derivatives to a new HepMatrix
  if (!fitState.hasMeasurements() || !fitState.parameters) {
    return nullptr;
  }
 
  int numberParameters = 5;
  if (m_allParameters) {
    numberParameters = fitState.parameters->numberParameters();
  }
  int rows = 0;
 
  for (const FitMeasurement* m : fitState.getMeasurements()) {
    if (!m->isPositionMeasurement()) {
      continue;
    }
    rows += m->numberDoF();
  }
 
  if (!numberParameters || !rows) {
    return nullptr;
  }
 
  ATH_MSG_VERBOSE(" DerivMatrix : " << fitState.getMeasurements().size()
                                    << " measurement objects giving " << rows
                                    << " rows and " << numberParameters
                                    << " columns (parameters)");
 
  auto derivativeMatrix =
      std::make_unique<Amg::MatrixX>(rows, numberParameters);
  int row = 0;
  for (const FitMeasurement* m : fitState.getMeasurements()) {
    if (!m->numberDoF() || !m->isPositionMeasurement()) {
      continue;
    }
    double norm = 0.;
    if (m->weight() > 0.) {
      norm = 1. / m->weight();
    }
 
    for (int col = 0; col < numberParameters; ++col) {
      (*derivativeMatrix)(row, col) = norm * m->derivative(col);
    }
 
    // take care of units for momentum derivs
    (*derivativeMatrix)(row, 4) *= Gaudi::Units::TeV;
    if (fitState.parameters->fitEnergyDeposit()) {
      (*derivativeMatrix)(row, 5) *= Gaudi::Units::TeV;
    }
    ++row;
    if (m->numberDoF() < 2) {
      continue;
    }
 
    // pixel measurements
    norm = 0.;
    if (m->weight2() > 0.)
      norm = 1. / m->weight2();
    for (int col = 0; col < numberParameters; ++col) {
      (*derivativeMatrix)(row, col) = norm * m->derivative2(col);
    }
    (*derivativeMatrix)(row, 4) *= Gaudi::Units::TeV;
    if (fitState.parameters->fitEnergyDeposit()) {
      (*derivativeMatrix)(row, 5) *= Gaudi::Units::TeV;
    }
    ++row;
  }
 
  if (row != rows) {
    ATH_MSG_WARNING("iPatGlobalFitter: inconsistent #rows in deriv matrix ");
  }
 
  return derivativeMatrix;
}

fullCovarianceMatrix()

std::unique_ptr< Amg::MatrixX > Trk::iPatGlobalFitter::fullCovarianceMatrix ( const FitState & fitState ) const

private

access to the global fitter's full covariance matrix

Definition at line 120 of file iPatGlobalFitter.cxx.

                                    {
  int numberParameters = 5;
 
  if (m_allParameters) {
    numberParameters = fitState.parameters->numberParameters();
  }
  ATH_MSG_VERBOSE(" FullCovarianceMatrix for " << numberParameters
                                               << " parameters");
 
  return std::make_unique<Amg::MatrixX>(Trk::FitProcedure::fullCovariance()->block(
      0, 0, numberParameters, numberParameters));
}

iterationsOfLastFit()

int Trk::iPatGlobalFitter::iterationsOfLastFit ( const FitState & fitState )

staticprivate

access to the number of iterations taken by the latest track fit

Definition at line 134 of file iPatGlobalFitter.cxx.

                                                                  {
  return fitState.iterations;
}

setMinIterations()

void Trk::iPatGlobalFitter::setMinIterations ( int minIterations )

private

set method for the minimum number of iterations for (alignment) friend

Definition at line 138 of file iPatGlobalFitter.cxx.

                                                         {
  m_fitProcedure->setMinIterations(minIterations);
}

Member Data Documentation

m_allParameters

Gaudi::Property<bool> Trk::iPatGlobalFitter::m_allParameters

private

Initial value:

{
    this, "AllParameters", false, "all or 5 parameters for matrix methods"}

Definition at line 51 of file iPatGlobalFitter.h.

                                     {
    this, "AllParameters", false, "all or 5 parameters for matrix methods"};

---

The documentation for this class was generated from the following files:

• iPatGlobalFitter.h
• iPatGlobalFitter.cxx