AlignmentCorrelationsStore Class Reference

#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

Correlations

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

Definition at line 19 of file AlignmentCorrelationsStore.h.

CorrelationsTable

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

Definition at line 18 of file AlignmentCorrelationsStore.h.

Constructor & Destructor Documentation

AlignmentCorrelationsStore()

AlignmentCorrelationsStore::AlignmentCorrelationsStore

(

void

)

Definition at line 7 of file AlignmentCorrelationsStore.cc.

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

~AlignmentCorrelationsStore()

virtual AlignmentCorrelationsStore::~AlignmentCorrelationsStore ( void  )

inlinevirtual

Definition at line 23 of file AlignmentCorrelationsStore.h.

{}

Member Function Documentation

correlations()

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 12 of file AlignmentCorrelationsStore.cc.

References cuy::col, fillCovariance(), fillCovarianceT(), edm::swap(), theCorrelations, and geometryDiff::transpose().

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

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

correlationsAvailable()

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 100 of file AlignmentCorrelationsStore.cc.

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

Referenced by AlignmentParameterStore::attachCorrelations().

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

fillCorrelationsTable()

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

private

Definition at line 135 of file AlignmentCorrelationsStore.cc.

References Alignable::alignmentParameters(), cuy::col, AlignmentParameters::numSelected(), readFromCovariance(), readFromCovarianceT(), TableParser::table, and geometryDiff::transpose().

Referenced by setCorrelations().

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

fillCovariance()

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

private

Definition at line 159 of file AlignmentCorrelationsStore.cc.

References cuy::col, and mps_splice::entry.

Referenced by correlations().

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

fillCovarianceT()

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

private

Definition at line 169 of file AlignmentCorrelationsStore.cc.

References cuy::col, and mps_splice::entry.

Referenced by correlations().

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

readFromCovariance()

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

private

Definition at line 179 of file AlignmentCorrelationsStore.cc.

References cuy::col, and mps_splice::entry.

Referenced by fillCorrelationsTable().

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

readFromCovarianceT()

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

private

Definition at line 189 of file AlignmentCorrelationsStore.cc.

References cuy::col, and mps_splice::entry.

Referenced by fillCorrelationsTable().

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

resetCorrelations()

void AlignmentCorrelationsStore::resetCorrelations ( void  )

virtual

Reset correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 114 of file AlignmentCorrelationsStore.cc.

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

Referenced by AlignmentParameterStore::resetParameters().

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

setCorrelations() [1/2]

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 51 of file AlignmentCorrelationsStore.cc.

References cuy::col, fillCorrelationsTable(), edm::swap(), theCorrelations, and geometryDiff::transpose().

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

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

setCorrelations() [2/2]

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

virtual

Set correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 85 of file AlignmentCorrelationsStore.cc.

References edm::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;
  }
}

size()

unsigned int AlignmentCorrelationsStore::size ( void  ) const

virtual

Get number of stored correlations.

Reimplemented in AlignmentExtendedCorrelationsStore.

Definition at line 126 of file AlignmentCorrelationsStore.cc.

References theCorrelations.

Referenced by ntupleDataFormat._Collection::__iter__(), ntupleDataFormat._Collection::__len__(), and AlignmentParameterStore::numCorrelations().

                                                        {
  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

theCorrelations

Correlations AlignmentCorrelationsStore::theCorrelations

private

Definition at line 66 of file AlignmentCorrelationsStore.h.

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