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;
  }