UncertainParameters.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Oliver Frost                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <tracking/trackFindingCDC/numerics/CovarianceMatrix.h>
#include <tracking/trackFindingCDC/numerics/PrecisionMatrix.h>
#include <tracking/trackFindingCDC/numerics/JacobianMatrix.h>
#include <tracking/trackFindingCDC/numerics/ParameterVector.h>
 
namespace Belle2 {
  namespace TrackFindingCDC {
 
    template <class T, class AEParameters>
    struct UncertainParametersUtil {
 
      using EParamaters = AEParameters;
 
      static const size_t c_N = EParamaters::c_N;
 
      using ParameterVector = TrackFindingCDC::ParameterVector<c_N>;
 
      using CovarianceMatrix = TrackFindingCDC::CovarianceMatrix<c_N>;
 
      using PrecisionMatrix = TrackFindingCDC::PrecisionMatrix<c_N>;
 
      using JacobianMatrix = TrackFindingCDC::JacobianMatrix<c_N>;
 
      static void transport(const JacobianMatrix& jacobian, CovarianceMatrix& cov);
 
      static CovarianceMatrix transported(const JacobianMatrix& jacobian,
                                          const CovarianceMatrix& cov);
 
      static JacobianMatrix reversalJacobian();
 
      static void reverse(CovarianceMatrix& cov);
 
      static CovarianceMatrix reversed(const CovarianceMatrix& cov);
 
      static CovarianceMatrix identity();
 
      static CovarianceMatrix covarianceFromPrecision(const PrecisionMatrix& prec);
 
      static PrecisionMatrix precisionFromCovariance(const CovarianceMatrix& cov);
 
      static CovarianceMatrix covarianceFromFullPrecision(const PrecisionMatrix& prec);
 
      static PrecisionMatrix precisionFromFullCovariance(const CovarianceMatrix& cov);
 
      template <class AParameterVector, int I = 0>
      static AParameterVector getSubParameterVector(const ParameterVector& par);
 
      template <class ACovarianceMatrix, int I = 0>
      static ACovarianceMatrix getSubCovarianceMatrix(const CovarianceMatrix& cov);
    };
  }
}