GsfMatrixTools Namespace Reference

Functions
template<typename T , unsigned int N>
T	trace (ROOT::Math::SMatrix< T, N, N > const &matrix)
template<typename T , unsigned int N>
T	trace (ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &a, ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &b)

Function Documentation

template<typename T , unsigned int N>

T GsfMatrixTools::trace

(

ROOT::Math::SMatrix< T, N, N > const &

matrix

)

Definition at line 18 of file GsfMatrixTools.h.

References mps_fire::i, makeMuonMisalignmentScenario::matrix, N, and mps_fire::result.

                                                  {
    T result = 0;
    for (unsigned int i = 0; i < N; ++i)
      result += matrix(i, i);
    return result;
  }

template<typename T , unsigned int N>

T GsfMatrixTools::trace	(	ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &	a,
		ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &	b
	)

Definition at line 29 of file GsfMatrixTools.h.

References a, b, mps_fire::i, and N.

                                                                       {
    auto i1 = a.begin();
    auto e1 = a.end();
    auto i2 = b.begin();

    T res = 0;
    // sum of the lower triangle;
    for (; i1 != e1; i1++, i2++)
      res += (*i1) * (*i2);
    res *= T(2.);
    // remove the duplicated diagonal...
    for (unsigned int i = 0; i < N; ++i)
      res -= a(i, i) * b(i, i);
    return res;
  }