#include <AlignmentCorrelationsStore.h>

Inheritance diagram for AlignmentCorrelationsStore:

Public Types
typedef std::map< Alignable , CorrelationsTable >	Correlations

typedef std::map< Alignable *, AlgebraicMatrix >	CorrelationsTable

Public Member Functions
	AlignmentCorrelationsStore (void)

virtual void	correlations (Alignable ap1, Alignable ap2, AlgebraicSymMatrix &cov, int row, int col) const

virtual bool	correlationsAvailable (Alignable ap1, Alignable ap2) const
	Check whether correlations are stored for a given pair of alignables. More...

virtual void	resetCorrelations (void)
	Reset correlations. More...

virtual void	setCorrelations (Alignable ap1, Alignable ap2, const AlgebraicSymMatrix &cov, int row, int col)

virtual void	setCorrelations (Alignable ap1, Alignable ap2, AlgebraicMatrix &mat)
	Set correlations. More...

virtual unsigned int	size (void) const
	Get number of stored correlations. More...

virtual	~AlignmentCorrelationsStore (void)

Private Member Functions
void	fillCorrelationsTable (Alignable ap1, Alignable ap2, CorrelationsTable *table, const AlgebraicSymMatrix &cov, int row, int col, bool transpose)

void	fillCovariance (Alignable ap1, Alignable ap2, const AlgebraicMatrix &entry, AlgebraicSymMatrix &cov, int row, int col) const

void	fillCovarianceT (Alignable ap1, Alignable ap2, const AlgebraicMatrix &entry, AlgebraicSymMatrix &cov, int row, int col) const

void	readFromCovariance (Alignable ap1, Alignable ap2, AlgebraicMatrix &entry, const AlgebraicSymMatrix &cov, int row, int col)

void	readFromCovarianceT (Alignable ap1, Alignable ap2, AlgebraicMatrix &entry, const AlgebraicSymMatrix &cov, int row, int col)

Private Attributes
Correlations	theCorrelations

Detailed Description

Definition at line 16 of file AlignmentCorrelationsStore.h.

Member Typedef Documentation

typedef std::map< Alignable*, CorrelationsTable* > AlignmentCorrelationsStore::Correlations

Definition at line 22 of file AlignmentCorrelationsStore.h.

typedef std::map< Alignable*, AlgebraicMatrix > AlignmentCorrelationsStore::CorrelationsTable

Definition at line 21 of file AlignmentCorrelationsStore.h.

Constructor & Destructor Documentation

AlignmentCorrelationsStore::AlignmentCorrelationsStore ( void )

Definition at line 8 of file AlignmentCorrelationsStore.cc.

 {
   edm::LogInfo("Alignment") << "@SUB=AlignmentCorrelationsStore::AlignmentCorrelationsStore "
                             << "\nCreated.";
 }

virtual AlignmentCorrelationsStore::~AlignmentCorrelationsStore ( void )

inlinevirtual

Definition at line 26 of file AlignmentCorrelationsStore.h.

References cuy::col, correlations(), correlationsAvailable(), mps_splice::entry, fillCorrelationsTable(), fillCovariance(), fillCovarianceT(), readFromCovariance(), readFromCovarianceT(), resetCorrelations(), setCorrelations(), size(), and geometryDiff::transpose().

26 {}

Member Function Documentation

void AlignmentCorrelationsStore::correlations	(	Alignable *	ap1,
		Alignable *	ap2,
		AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)		const

virtual

Write correlations directly to the covariance matrix starting at the given position. Indices are assumed to start from 0.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 15 of file AlignmentCorrelationsStore.cc.

References fillCovariance(), fillCovarianceT(), std::swap(), theCorrelations, and geometryDiff::transpose().

Referenced by resetCorrelations(), AlignmentParameterStore::selectParameters(), and ~AlignmentCorrelationsStore().

 {
   static Alignable* previousAlignable = nullptr;
   static CorrelationsTable* previousCorrelations;
 
   // Needed by 'resetCorrelations()' to reset the static pointer:
   if ( ap1 == nullptr ) { previousAlignable = nullptr; return; }
 
   bool transpose = ( ap2 > ap1 );
   if ( transpose ) std::swap( ap1, ap2 ); 
 
   if ( ap1 == previousAlignable )
   {
     CorrelationsTable::const_iterator itC2 = previousCorrelations->find( ap2 );
     if ( itC2 != previousCorrelations->end() )
     {
       transpose ?
     fillCovarianceT( ap1, ap2, (*itC2).second, cov, row, col ) :
     fillCovariance( ap1, ap2, (*itC2).second, cov, row, col );
     }
   }
   else
   {
     Correlations::const_iterator itC1 = theCorrelations.find( ap1 );
     if ( itC1 != theCorrelations.end() )
     {
       previousAlignable = ap1;
       previousCorrelations = (*itC1).second;
 
       CorrelationsTable::const_iterator itC2 = (*itC1).second->find( ap2 );
       if ( itC2 != (*itC1).second->end() )
       {
     transpose ?
       fillCovarianceT( ap1, ap2, (*itC2).second, cov, row, col ) :
       fillCovariance( ap1, ap2, (*itC2).second, cov, row, col );
       }
     }
   }
 
   // don't fill anything into the covariance if there's no entry
   return;
 }

bool AlignmentCorrelationsStore::correlationsAvailable	(	Alignable *	ap1,
		Alignable *	ap2
	)		const

virtual

Check whether correlations are stored for a given pair of alignables.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 117 of file AlignmentCorrelationsStore.cc.

References std::swap(), theCorrelations, and geometryDiff::transpose().

Referenced by AlignmentParameterStore::attachCorrelations(), and ~AlignmentCorrelationsStore().

 {
   bool transpose = ( ap2 > ap1 );
   if ( transpose ) std::swap( ap1, ap2 );
 
   Correlations::const_iterator itC1 = theCorrelations.find( ap1 );
   if ( itC1 != theCorrelations.end() )
   {
     CorrelationsTable::const_iterator itC2 = itC1->second->find( ap2 );
     if ( itC2 != itC1->second->end() ) return true;
   }
   return false;
 }

void AlignmentCorrelationsStore::fillCorrelationsTable	(	Alignable *	ap1,
		Alignable *	ap2,
		CorrelationsTable *	table,
		const AlgebraicSymMatrix &	cov,
		int	row,
		int	col,
		bool	transpose
	)

private

Definition at line 157 of file AlignmentCorrelationsStore.cc.

References Alignable::alignmentParameters(), AlignmentParameters::numSelected(), readFromCovariance(), and readFromCovarianceT().

Referenced by setCorrelations(), and ~AlignmentCorrelationsStore().

 {
   CorrelationsTable::iterator itC = table->find( ap2 );
 
   if ( itC != table->end() )
   {
     transpose ?
       readFromCovarianceT( ap1, ap2, itC->second, cov, row, col ) :
       readFromCovariance( ap1, ap2, itC->second, cov, row, col );
   }
   else
   {
     int nRow = ap1->alignmentParameters()->numSelected();
     int nCol = ap2->alignmentParameters()->numSelected();
     AlgebraicMatrix newEntry( nRow, nCol );
 
     transpose ?
       readFromCovarianceT( ap1, ap2, newEntry, cov, row, col ) :
       readFromCovariance( ap1, ap2, newEntry, cov, row, col );
 
     (*table)[ap2] = newEntry;
   }
 }

void AlignmentCorrelationsStore::fillCovariance	(	Alignable *	ap1,
		Alignable *	ap2,
		const AlgebraicMatrix &	entry,
		AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)		const

private

Definition at line 184 of file AlignmentCorrelationsStore.cc.

Referenced by correlations(), and ~AlignmentCorrelationsStore().

 {
   int nRow = entry.num_row();
   int nCol = entry.num_col();
 
   for ( int iRow = 0; iRow < nRow; ++iRow )
     for ( int jCol = 0; jCol < nCol; ++jCol )
       cov[row+iRow][col+jCol] = entry[iRow][jCol];
 }

void AlignmentCorrelationsStore::fillCovarianceT	(	Alignable *	ap1,
		Alignable *	ap2,
		const AlgebraicMatrix &	entry,
		AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)		const

private

Definition at line 197 of file AlignmentCorrelationsStore.cc.

Referenced by correlations(), and ~AlignmentCorrelationsStore().

 {
   int nRow = entry.num_row();
   int nCol = entry.num_col();
 
   for ( int iRow = 0; iRow < nRow; ++iRow )
     for ( int jCol = 0; jCol < nCol; ++jCol )
       cov[row+jCol][col+iRow] = entry[iRow][jCol];
 }

void AlignmentCorrelationsStore::readFromCovariance	(	Alignable *	ap1,
		Alignable *	ap2,
		AlgebraicMatrix &	entry,
		const AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)

private

Definition at line 210 of file AlignmentCorrelationsStore.cc.

Referenced by fillCorrelationsTable(), and ~AlignmentCorrelationsStore().

 {
   int nRow = entry.num_row();
   int nCol = entry.num_col();
 
   for ( int iRow = 0; iRow < nRow; ++iRow )
     for ( int jCol = 0; jCol < nCol; ++jCol )
       entry[iRow][jCol] = cov[row+iRow][col+jCol];
 }

void AlignmentCorrelationsStore::readFromCovarianceT	(	Alignable *	ap1,
		Alignable *	ap2,
		AlgebraicMatrix &	entry,
		const AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)

private

Definition at line 223 of file AlignmentCorrelationsStore.cc.

Referenced by fillCorrelationsTable(), and ~AlignmentCorrelationsStore().

 {
   int nRow = entry.num_row();
   int nCol = entry.num_col();
 
   for ( int iRow = 0; iRow < nRow; ++iRow )
     for ( int jCol = 0; jCol < nCol; ++jCol )
       entry[iRow][jCol] = cov[row+jCol][col+iRow];
 }

void AlignmentCorrelationsStore::resetCorrelations ( void )

virtual

Reset correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 132 of file AlignmentCorrelationsStore.cc.

References correlations(), setCorrelations(), and theCorrelations.

Referenced by AlignmentParameterStore::resetParameters(), and ~AlignmentCorrelationsStore().

 {
   Correlations::iterator itC;
   for ( itC = theCorrelations.begin(); itC != theCorrelations.end(); ++itC ) delete (*itC).second;
   theCorrelations.erase( theCorrelations.begin(), theCorrelations.end() );
 
   // Reset the static pointers to the 'previous alignables'
   AlgebraicSymMatrix dummy;
   correlations( nullptr, nullptr, dummy, 0, 0 );
   setCorrelations( nullptr, nullptr, dummy, 0, 0 );
 }

void AlignmentCorrelationsStore::setCorrelations	(	Alignable *	ap1,
		Alignable *	ap2,
		const AlgebraicSymMatrix &	cov,
		int	row,
		int	col
	)

virtual

Get correlations directly from the given position of the covariance matrix and store them. Indices are assumed to start from 0.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 60 of file AlignmentCorrelationsStore.cc.

References fillCorrelationsTable(), std::swap(), theCorrelations, and geometryDiff::transpose().

Referenced by AlignmentParameterStore::attachCorrelations(), resetCorrelations(), AlignmentParameterStore::updateParameters(), and ~AlignmentCorrelationsStore().

 {
   static Alignable* previousAlignable = nullptr;
   static CorrelationsTable* previousCorrelations;
 
   // Needed by 'resetCorrelations()' to reset the static pointer:
   if ( ap1 == nullptr ) { previousAlignable = nullptr; return; }
 
   bool transpose = ( ap2 > ap1 );
   if ( transpose ) std::swap( ap1, ap2 );
 
   if ( ap1 == previousAlignable )
   {
     fillCorrelationsTable( ap1, ap2, previousCorrelations, cov, row, col, transpose );
   }
   else
   {
     Correlations::iterator itC = theCorrelations.find( ap1 );
     if ( itC != theCorrelations.end() )
     {
       fillCorrelationsTable( ap1, ap2, itC->second, cov, row, col, transpose );
       previousAlignable = ap1;
       previousCorrelations = itC->second;
     }
     else
     {
       CorrelationsTable* newTable = new CorrelationsTable;
       fillCorrelationsTable( ap1, ap2, newTable, cov, row, col, transpose );
 
       theCorrelations[ap1] = newTable;
       previousAlignable = ap1;
       previousCorrelations = newTable;
     }
   }
 }

void AlignmentCorrelationsStore::setCorrelations	(	Alignable *	ap1,
		Alignable *	ap2,
		AlgebraicMatrix &	mat
	)

virtual

Set correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 98 of file AlignmentCorrelationsStore.cc.

References std::swap(), theCorrelations, and geometryDiff::transpose().

 {
   bool transpose = ( ap2 > ap1 );
   if ( transpose ) std::swap( ap1, ap2 );
 
   Correlations::iterator itC1 = theCorrelations.find( ap1 );
   if ( itC1 != theCorrelations.end() )
   {
     (*itC1->second)[ap2] = transpose ? mat.T() : mat;
   }
   else
   {
     CorrelationsTable* newTable = new CorrelationsTable;
     (*newTable)[ap2] = transpose ? mat.T() : mat;
     theCorrelations[ap1] = newTable;
   }
 }

unsigned int AlignmentCorrelationsStore::size ( void ) const

virtual

Get number of stored correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 145 of file AlignmentCorrelationsStore.cc.

References theCorrelations.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and ~AlignmentCorrelationsStore().

 {
   unsigned int size = 0;
   Correlations::const_iterator itC;
   for ( itC = theCorrelations.begin(); itC != theCorrelations.end(); ++itC )
     size += itC->second->size();
 
   return size;
 }

Member Data Documentation

Correlations AlignmentCorrelationsStore::theCorrelations

private

Definition at line 67 of file AlignmentCorrelationsStore.h.

Referenced by correlations(), correlationsAvailable(), resetCorrelations(), setCorrelations(), and size().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation