PndEmcErrorMatrix.cxx File Reference

#include "PndEmcErrorMatrix.h"
#include "PndEmcErrorMatrixPar.h"
#include "TSystem.h"
#include <cmath>

Go to the source code of this file.

Functions
TMatrixD	similarityWith (const TMatrixD &mat, const TMatrixD &m1)

Function Documentation

TMatrixD similarityWith	(	const TMatrixD &	mat,
		const TMatrixD &	m1
	)

Definition at line 305 of file PndEmcErrorMatrix.cxx.

References c, and r.

Referenced by PndEmcErrorMatrix::Get4MomentumErrorMatrix(), PndEmcErrorMatrix::GetErrorMatrix(), and PndEmcErrorMatrix::GetErrorP7().

{
        TMatrixD result(m1.GetNrows(),m1.GetNrows());
        
        TMatrixD temp = m1*mat;
        double tmp;
        
        for (int r = 0; r < m1.GetNrows(); r++) {
                for (int c = 0; c <= r; c++) {
                        tmp = 0.;
                        for (int k = 0; k < m1.GetNcols(); k++) {
                                tmp += temp(r,k)*m1(c,k);
                        }
                        result(r,c) = tmp;
                        // Modification from original code, to make output matrix explicitly symmetric
                        // In original babar code symmetric matrix were stored as lower triangular.
                        if (r!=c)
                                result(c,r) = tmp;
                }
        }
        return result;
}