import ROOT,math

ROOT.gROOT.ProcessLine(".x rootlogon.C") 
ROOT.gROOT.ProcessLine('gROOT->SetStyle("Plain")')

cov=ROOT.TMatrixDSym(3)
cov[0][0]=1
cov[1][1]=0
cov[2][2]=3

cov_i=ROOT.TMatrixDSym(cov)
cov_i.Invert()
print "the cov"
cov.Print()


eigenproblem=ROOT.TMatrixDSymEigen(cov)
eigenwert=eigenproblem.GetEigenValues()
eigenvect=eigenproblem.GetEigenVectors()

print "eigenwert:"
eigenwert.Print()

print "eigenvectors:"
eigenvect.Print()


Axis=ROOT.TVectorD(3)
Axis[0]=0
Axis[1]=0
Axis[2]=1

cut=cov_i.Similarity(Axis)

print 'cut:',cut,'1/cut',1/cut,'sqrt(cut)',math.sqrt(cut),'cut^2',cut**2