RhoError Class Reference

#include <RhoError.h>

Inheritance diagram for RhoError:

Public Member Functions
	RhoError (Int_t n, Double_t init=0.0)
	RhoError (const TMatrixD &p)
	RhoError (const RhoError &v)
RhoError &	operator= (const RhoError &v)
RhoError &	operator= (const TMatrixD &v)
virtual	~RhoError ()
Double_t	DetermineChisq (TVectorD &diff)
RhoError &	operator*= (Double_t t)
RhoError &	operator- ()
RhoError	Similarity (const TMatrixD &m1) const
RhoError &	SimilarityWith (const RhoError &m, const TMatrixD &m1)
RhoError	Similarity (const TRotation &rot) const
RhoError	Similarity (const TLorentzRotation &rot) const
Double_t	Similarity (TVectorD &m1)
TMatrixD	SimilarityT (TMatrixD &m1)

Static Public Attributes
static const Double_t	chisqUndef = -1.

Friends
RhoError	operator* (Double_t t, const RhoError &m1)
RhoError	operator* (const RhoError &m1, Double_t t)
RhoError	operator+ (const RhoError &m1, const RhoError &m2)
RhoError	operator- (const RhoError &m1, const RhoError &m2)
std::ostream &	operator<< (std::ostream &out, const RhoError &mat)
std::istream &	operator>> (std::istream &in, RhoError &mat)

Detailed Description

Definition at line 34 of file RhoError.h.

Constructor & Destructor Documentation

RhoError::RhoError	(	Int_t	n,
		Double_t	init = 0.0
	)

Definition at line 25 of file RhoError.cxx.

References init, and n.

{
  TMatrixD::Allocate ( n,n );
  for ( int i1=0; i1<n; i1++ )
    for ( int i2=0; i2<n; i2++ ) {
      TMatrixD::operator() ( i1,i2 ) = init;
    }
}

RhoError::RhoError

(

const TMatrixD &

p

)

Definition at line 41 of file RhoError.cxx.

{
  TMatrixD::Allocate ( p.GetNrows(),p.GetNcols() );
  TMatrixD::operator= ( p );
}

RhoError::RhoError

(

const RhoError &

v

)

Definition at line 34 of file RhoError.cxx.

                                       :
TMatrixD(v)
{
  TMatrixD::Allocate ( v.GetNrows(),v.GetNcols() );
  TMatrixD::operator= ( v );
}

virtual RhoError::~RhoError ( )

inlinevirtual

Definition at line 56 of file RhoError.h.

{};

Member Function Documentation

Double_t RhoError::DetermineChisq

(

TVectorD &

diff

)

Definition at line 166 of file RhoError.cxx.

References Double_t, i, and SimilarityT().

Referenced by RhoVector3Err::DetermineChisq(), and RhoLorentzVectorErr::DetermineChisq().

{
  TMatrixD dMat ( diff.GetNrows(),diff.GetNrows() );
  dMat.UnitMatrix();

  for ( int i = 0; i < diff.GetNrows(); i++ ) { dMat ( i,0 ) = diff ( i ); }

  Invert();
  Double_t chisq = SimilarityT ( dMat ) ( 0,0 );

  return chisq;
}

RhoError & RhoError::operator*=

(

Double_t

t

)

Definition at line 61 of file RhoError.cxx.

References operator*=().

{
  return ( RhoError& ) TMatrixD::operator*= ( t );
}

RhoError & RhoError::operator-

(

)

Definition at line 81 of file RhoError.cxx.

{
  return *this;
}

RhoError & RhoError::operator=

(

const RhoError &

v

)

Definition at line 47 of file RhoError.cxx.

{
  TMatrixD::ResizeTo ( v );
  TMatrixD::operator= ( v );
  return *this;
}

RhoError & RhoError::operator=

(

const TMatrixD &

v

)

Definition at line 54 of file RhoError.cxx.

{
  TMatrixD::ResizeTo ( v );
  TMatrixD::operator= ( v );
  return *this;
}

RhoError RhoError::Similarity

(

const TMatrixD &

m1

)

const

Definition at line 104 of file RhoError.cxx.

References SimilarityWith().

Referenced by RhoCandidate::Set(), Similarity(), RhoVector3Err::Transform(), and RhoLorentzVectorErr::Transform().

{
  RhoError mret ( m1.GetNrows() );
  mret.ResizeTo ( m1.GetNrows(),m1.GetNrows() );
  mret.SimilarityWith ( *this, m1 );
  return mret;
}

RhoError RhoError::Similarity

(

const TRotation &

rot

)

const

Definition at line 88 of file RhoError.cxx.

References Similarity().

{
  TMatrixD mat ( 3,3 );
  mat ( 0,0 ) =rot.XX();
  mat ( 0,1 ) =rot.XY();
  mat ( 0,2 ) =rot.XZ();
  mat ( 1,0 ) =rot.YX();
  mat ( 1,1 ) =rot.YY();
  mat ( 1,2 ) =rot.YZ();
  mat ( 2,0 ) =rot.ZX();
  mat ( 2,1 ) =rot.ZY();
  mat ( 2,2 ) =rot.ZZ();

  return Similarity ( mat );
}

RhoError RhoError::Similarity

(

const TLorentzRotation &

rot

)

const

Definition at line 112 of file RhoError.cxx.

References Similarity().

{
  TMatrixD mat ( 4,4 );
  mat ( 0,0 ) =rot.XX();
  mat ( 0,1 ) =rot.XY();
  mat ( 0,2 ) =rot.XZ();
  mat ( 0,3 ) =rot.XT();
  mat ( 1,0 ) =rot.YX();
  mat ( 1,1 ) =rot.YY();
  mat ( 1,2 ) =rot.YZ();
  mat ( 1,3 ) =rot.YT();
  mat ( 2,0 ) =rot.ZX();
  mat ( 2,1 ) =rot.ZY();
  mat ( 2,2 ) =rot.ZZ();
  mat ( 2,3 ) =rot.ZT();
  mat ( 3,0 ) =rot.TX();
  mat ( 3,1 ) =rot.TY();
  mat ( 3,2 ) =rot.TZ();
  mat ( 3,3 ) =rot.TT();

  return Similarity ( mat );
}

Double_t RhoError::Similarity

(

TVectorD &

m1

)

Definition at line 264 of file RhoError.cxx.

References a, b, and Double_t.

{
  Double_t mret = 0.0;
  TVectorD temp ( m1 );
  temp*= ( *this );
  // If m1*(*this) has correct dimensions, then so will the m1.T multiplication.
  // So there is no need to check dimensions again.
  Double_t* a=&temp ( 0 );
  Double_t* b=&m1 ( 0 );
  Double_t* e=a+m1.GetNrows();
  for ( ; a<e; ) { mret += ( * ( a++ ) ) * ( * ( b++ ) ); }
  return mret;
}

TMatrixD RhoError::SimilarityT

(

TMatrixD &

m1

)

Definition at line 278 of file RhoError.cxx.

References c, Double_t, i, n, and r.

Referenced by DetermineChisq().

{
  TMatrixD mret;
  mret.ResizeTo ( m1 );
  TMatrixD::Mult ( *this,m1 );
  int n = m1.GetNcols();
  Double_t* mrc = &mret ( 0,0 );
  Double_t* temp1r = &TMatrixD::operator() ( 0,0 );
  for ( int r=1; r<=mret.GetNrows(); r++ ) {
    Double_t* m11c = &m1 ( 0,0 );
    for ( int c=1; c<=r; c++ ) {
      Double_t tmp = 0.0;
      Double_t* tempir = temp1r;
      Double_t* m1ic = m11c;
      for ( int i=1; i<=m1.GetNrows(); i++ ) {
        tmp+= ( * ( tempir ) ) * ( * ( m1ic ) );
        tempir += n;
        m1ic += n;
      }
      * ( mrc++ ) = tmp;
      m11c++;
    }
    temp1r++;
  }
  return mret;
}

RhoError & RhoError::SimilarityWith	(	const RhoError &	m,
		const TMatrixD &	m1
	)

Definition at line 140 of file RhoError.cxx.

References c, Double_t, and r.

Referenced by Similarity().

{
  assert ( GetNrows() == m1.GetNrows() );
  TMatrixD temp;
  temp.ResizeTo ( m1.GetNrows(),m1.GetNrows() );
  temp.Mult ( m1,mat );
  Double_t tmp;

  for ( int r = 0; r < 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 );
      }
      ( *this ) ( r,c ) = tmp;
    }
  }
  return *this;
}

Friends And Related Function Documentation

RhoError operator* ( Double_t  t, const RhoError &  m1  )

friend

Definition at line 213 of file RhoError.cxx.

{
  RhoError mret ( m1 );
  //mret.ResizeTo ( m1 );
  //mret = m1;
  mret *= t;
  return mret;
}

RhoError operator* ( const RhoError &  m1, Double_t  t  )

friend

Definition at line 222 of file RhoError.cxx.

{
  RhoError mret ( m1 );
  //mret.ResizeTo ( m1 );
  //mret = m1;
  mret *= t;
  return mret;
}

RhoError operator+ ( const RhoError &  m1, const RhoError &  m2  )

friend

Definition at line 245 of file RhoError.cxx.

{
  RhoError mret ( m1 );
  //mret.ResizeTo ( m1 );
  //mret = m1;
  mret += m2;
  //std::cout<<" -- Adding two matrices: -- m1:"<<&m1<<"  m2:"<<&m2<<std::endl;
  return mret;
}

RhoError operator- ( const RhoError &  m1, const RhoError &  m2  )

friend

Definition at line 255 of file RhoError.cxx.

{
  RhoError mret ( m1 );
  //mret.ResizeTo ( m1 );
  mret -= m2;
  return mret;
}

std::ostream& operator<< ( std::ostream &  out, const RhoError &  mat  )

friend

Definition at line 179 of file RhoError.cxx.

{
  out << "Bbr Covariance Matrix:";
  out << ( TMatrixD& ) mat;
  return out;
}

std::istream& operator>> ( std::istream &  in, RhoError &  mat  )

friend

Definition at line 186 of file RhoError.cxx.

{
  // Peek at the next non-space character:
  char nextChar = ' ';
  while ( isspace ( nextChar ) ) {
    nextChar = in.get();
  }
  in.putback ( nextChar );

  if ( EOF != nextChar ) {
    if ( !isdigit ( nextChar ) ) {
      // Remove the "Bbr Covariance Matrix:" line:
      const int DUMMY_SIZE = 1000;
      char dummy[DUMMY_SIZE];
      in.getline ( dummy, DUMMY_SIZE );
    }
    // Read in the matrix:
    for ( int row = 1; row <= mat.GetNrows(); ++row ) {
      for ( int col = 1; col <= mat.GetNcols(); ++col ) {
        in >> mat ( row, col );
      }
    }
  }
  return in;
}

Member Data Documentation

const Double_t RhoError::chisqUndef = -1.

static

Definition at line 39 of file RhoError.h.

---

The documentation for this class was generated from the following files:

• RhoError.h
• RhoError.cxx