UncertainParameters.icc.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#pragma once
 
#include <tracking/trackingUtilities/geometry/UncertainParameters.h>
 
#include <tracking/trackingUtilities/numerics/CovarianceMatrixUtil.h>
#include <tracking/trackingUtilities/numerics/JacobianMatrixUtil.h>
#include <tracking/trackingUtilities/numerics/ParameterVectorUtil.h>
 
namespace Belle2 {
  namespace TrackingUtilities {
 
    template <class T, class AEParameters>
    void UncertainParametersUtil<T, AEParameters>::transport(const JacobianMatrix& jacobian,
                                                             CovarianceMatrix& cov)
    {
      CovarianceMatrixUtil::transport(jacobian, cov);
    }
 
    template <class T, class AEParameters>
    auto UncertainParametersUtil<T, AEParameters>::transported(const JacobianMatrix& jacobian,
                                                               const CovarianceMatrix& cov)
    -> CovarianceMatrix
    {
      return CovarianceMatrixUtil::transported(jacobian, cov);
    }
 
    template <class T, class AEParameters>
    auto UncertainParametersUtil<T, AEParameters>::reversalJacobian() -> JacobianMatrix
    {
      return JacobianMatrixUtil::scale(T::reversalSigns());
    }
 
    template <class T, class AEParameters>
    void UncertainParametersUtil<T, AEParameters>::reverse(CovarianceMatrix& cov)
    {
      return transport(reversalJacobian(), cov);
    }
 
    template <class T, class AEParameters>
    auto UncertainParametersUtil<T, AEParameters>::reversed(const CovarianceMatrix& cov)
    -> CovarianceMatrix
    {
      return transported(reversalJacobian(), cov);
    }
 
    template <class T, class AEParameters>
    auto UncertainParametersUtil<T, AEParameters>::identity() -> CovarianceMatrix
    {
      return CovarianceMatrixUtil::identity<c_N>();
    }
 
    template <class T, class AEParameters>
    template <class AParameterVector, int I>
    AParameterVector
    UncertainParametersUtil<T, AEParameters>::getSubParameterVector(const ParameterVector& par)
    {
      return ParameterVectorUtil::getSub<AParameterVector, I>(par);
    }
 
    template <class T, class AEParameters>
    template <class ACovarianceMatrix, int I>
    ACovarianceMatrix
    UncertainParametersUtil<T, AEParameters>::getSubCovarianceMatrix(const CovarianceMatrix& cov)
    {
      return CovarianceMatrixUtil::getSub<ACovarianceMatrix, I>(cov);
    }
 
    template <class T, class AEParameters>
    auto
    UncertainParametersUtil<T, AEParameters>::covarianceFromPrecision(const PrecisionMatrix& prec)
    -> CovarianceMatrix
    {
      return CovarianceMatrixUtil::fromPrecision(prec);
    }
 
    template <class T, class AEParameters>
    auto
    UncertainParametersUtil<T, AEParameters>::precisionFromCovariance(const CovarianceMatrix& cov)
    -> PrecisionMatrix
    {
      return CovarianceMatrixUtil::toPrecision(cov);
    }
 
    template <class T, class AEParameters>
    auto
    UncertainParametersUtil<T, AEParameters>::covarianceFromFullPrecision(const PrecisionMatrix& prec)
    -> CovarianceMatrix
    {
      return CovarianceMatrixUtil::fromFullPrecision(prec);
    }
 
    template <class T, class AEParameters>
    auto
    UncertainParametersUtil<T, AEParameters>::precisionFromFullCovariance(const CovarianceMatrix& cov)
    -> PrecisionMatrix
    {
      return CovarianceMatrixUtil::fullToPrecision(cov);
    }
  }
}