CompositeAlignmentDerivativesExtractor Class Reference

#include <CompositeAlignmentDerivativesExtractor.h>

Public Member Functions
	CompositeAlignmentDerivativesExtractor (const align::Alignables &alignables, const std::vector< AlignableDet * > &alignableDets, const std::vector< TrajectoryStateOnSurface > &tsos)
	CompositeAlignmentDerivativesExtractor (const align::Alignables &alignables, const std::vector< AlignableDetOrUnitPtr > &alignableDets, const std::vector< TrajectoryStateOnSurface > &tsos)
	constructor More...
const AlgebraicVector &	correctionTerm (void) const
const AlgebraicMatrix &	derivatives (void) const
	~CompositeAlignmentDerivativesExtractor (void)
	destructor More...

Private Member Functions
void	extractCurrentAlignment (const align::Alignables &alignables, const std::vector< AlignableDetOrUnitPtr > &alignableDets, const std::vector< TrajectoryStateOnSurface > &tsos)
void	extractWithMultipleHits (const std::vector< AlgebraicVector > &subCorrectionTerm, const std::vector< AlgebraicMatrix > &subDerivatives, const align::Alignables &alignables)
void	extractWithoutMultipleHits (const std::vector< AlgebraicVector > &subCorrectionTerm, const std::vector< AlgebraicMatrix > &subDerivatives)

Private Attributes
AlgebraicVector	theCorrectionTerm
AlgebraicMatrix	theDerivatives

Detailed Description

A helper class to extract derivatives from composite alignable objects

Date

2007/05/02 21:01:52

Revision

1.7

(last update by

Author

fronga

)

Definition at line 20 of file CompositeAlignmentDerivativesExtractor.h.

Constructor & Destructor Documentation

CompositeAlignmentDerivativesExtractor() [1/2]

CompositeAlignmentDerivativesExtractor::CompositeAlignmentDerivativesExtractor	(	const align::Alignables &	alignables,
		const std::vector< AlignableDet * > &	alignableDets,
		const std::vector< TrajectoryStateOnSurface > &	tsos
	)

deprecated constructor for backward compatibility (use mor general AlignableDetOrUnitPtr)

Definition at line 11 of file CompositeAlignmentDerivativesExtractor.cc.

                                                     {
  std::vector<AlignableDetOrUnitPtr> detOrUnits;
  detOrUnits.reserve(alignableDets.size());
 
  std::vector<AlignableDet *>::const_iterator it, itEnd;
  for (it = alignableDets.begin(), itEnd = alignableDets.end(); it != itEnd; ++it)
    detOrUnits.push_back(AlignableDetOrUnitPtr(*it));
 
  extractCurrentAlignment(alignables, detOrUnits, tsos);
}

References extractCurrentAlignment().

CompositeAlignmentDerivativesExtractor() [2/2]

CompositeAlignmentDerivativesExtractor::CompositeAlignmentDerivativesExtractor	(	const align::Alignables &	alignables,
		const std::vector< AlignableDetOrUnitPtr > &	alignableDets,
		const std::vector< TrajectoryStateOnSurface > &	tsos
	)

constructor

Definition at line 26 of file CompositeAlignmentDerivativesExtractor.cc.

                                                     {
  extractCurrentAlignment(alignables, alignableDets, tsos);
}

References extractCurrentAlignment().

~CompositeAlignmentDerivativesExtractor()

CompositeAlignmentDerivativesExtractor::~CompositeAlignmentDerivativesExtractor ( void  )

inline

destructor

Definition at line 33 of file CompositeAlignmentDerivativesExtractor.h.

{};

Member Function Documentation

correctionTerm()

const AlgebraicVector& CompositeAlignmentDerivativesExtractor::correctionTerm ( void  ) const

inline

Definition at line 36 of file CompositeAlignmentDerivativesExtractor.h.

{ return theCorrectionTerm; }

References theCorrectionTerm.

Referenced by CompositeAlignmentParameters::correctionTerm().

derivatives()

const AlgebraicMatrix& CompositeAlignmentDerivativesExtractor::derivatives ( void  ) const

inline

Definition at line 35 of file CompositeAlignmentDerivativesExtractor.h.

{ return theDerivatives; }

References theDerivatives.

Referenced by CompositeAlignmentParameters::derivatives().

extractCurrentAlignment()

void CompositeAlignmentDerivativesExtractor::extractCurrentAlignment ( const align::Alignables &  alignables, const std::vector< AlignableDetOrUnitPtr > &  alignableDets, const std::vector< TrajectoryStateOnSurface > &  tsos  )

private

Definition at line 35 of file CompositeAlignmentDerivativesExtractor.cc.

                                                     {
  // sanity check
  if (alignables.size() != alignableDets.size()) {
    edm::LogError("CompositeAlignmentDerivativesExtractor")
        << "Inconsistent length of arguments: alignables=" << alignables.size()
        << ", alignableDets=" << alignableDets.size();
    return;
  }
 
  if (alignables.size() != tsos.size()) {
    edm::LogError("CompositeAlignmentDerivativesExtractor")
        << "Inconsistent length of arguments: alignables=" << alignables.size() << ", tsos=" << tsos.size();
    return;
  }
 
  align::Alignables::const_iterator itAlignable = alignables.begin();
  std::vector<AlignableDetOrUnitPtr>::const_iterator itAlignableDet = alignableDets.begin();
  std::vector<TrajectoryStateOnSurface>::const_iterator itTsos = tsos.begin();
 
  int nRow = 0;
  int nCollumn = 0;
  unsigned int nAlignables = 0;
 
  std::vector<AlgebraicMatrix> subDerivatives;
  std::vector<AlgebraicVector> subCorrectionTerm;
 
  // get the individual derivatives and correction term and determine the
  // dimension
  while (itAlignable != alignables.end()) {
    // Get the current estimate on the alignment parameters
    AlgebraicVector subAlignmentParameters = (*itAlignable)->alignmentParameters()->selectedParameters();
 
    // Get the derivatives or the local coordinates w.r.t. the corresponding
    // alignment parameters
    AlgebraicMatrix subAlignmentDerivatives =
        (*itAlignable)->alignmentParameters()->selectedDerivatives(*itTsos, *itAlignableDet);
 
    subDerivatives.push_back(subAlignmentDerivatives.T());
    subCorrectionTerm.push_back(subAlignmentDerivatives.T() * subAlignmentParameters);
 
    nRow += 2;
    // check if it is the first occurrence of this Alignable
    if (count(alignables.begin(), itAlignable, *itAlignable) == 0) {
      // matrix is transposed -> num_row() instead of num_col()
      nCollumn += subAlignmentDerivatives.num_row();
      nAlignables++;
    }
 
    ++itAlignable;
    ++itAlignableDet;
    ++itTsos;
  }
 
  // construct derivatives and correction term with the right dimension
  theDerivatives = AlgebraicMatrix(nRow, nCollumn, 0);
  theCorrectionTerm = AlgebraicVector(nRow, 0);
 
  if (alignables.size() == nAlignables)
    // One hit per alignable
    extractWithoutMultipleHits(subCorrectionTerm, subDerivatives);
  else
    // At least one alignable has two hits
    extractWithMultipleHits(subCorrectionTerm, subDerivatives, alignables);
 
  return;
}

References KineDebug3::count(), extractWithMultipleHits(), extractWithoutMultipleHits(), theCorrectionTerm, and theDerivatives.

Referenced by CompositeAlignmentDerivativesExtractor().

extractWithMultipleHits()

void CompositeAlignmentDerivativesExtractor::extractWithMultipleHits ( const std::vector< AlgebraicVector > &  subCorrectionTerm, const std::vector< AlgebraicMatrix > &  subDerivatives, const align::Alignables &  alignables  )

private

Definition at line 132 of file CompositeAlignmentDerivativesExtractor.cc.

                                       {
  std::vector<AlgebraicVector>::const_iterator itSubCorrectionTerm = subCorrectionTerm.begin();
  std::vector<AlgebraicMatrix>::const_iterator itSubDerivatives = subDerivatives.begin();
  align::Alignables::const_iterator itAlignables = alignables.begin();
  align::Alignables::const_iterator itPosition;
  align::Alignables::const_iterator itLastPosition;
 
  int iRow = 1;
 
  // Fill in the individual terms
  while (itAlignables != alignables.end()) {
    theCorrectionTerm.sub(iRow, *itSubCorrectionTerm);
 
    int iCollumn = 1;
    int iAlignable = 0;
 
    itLastPosition = find(alignables.begin(), itAlignables, *itAlignables);
 
    for (itPosition = alignables.begin(); itPosition != itLastPosition; ++itPosition) {
      if (count(alignables.begin(), itPosition, *itPosition) == 0)
        iCollumn += subDerivatives[iAlignable].num_col();
      iAlignable++;
    }
 
    theDerivatives.sub(iRow, iCollumn, *itSubDerivatives);
 
    iRow += 2;
 
    ++itAlignables;
    ++itSubCorrectionTerm;
    ++itSubDerivatives;
  }
 
  return;
}

References KineDebug3::count(), spr::find(), theCorrectionTerm, and theDerivatives.

Referenced by extractCurrentAlignment().

extractWithoutMultipleHits()

void CompositeAlignmentDerivativesExtractor::extractWithoutMultipleHits ( const std::vector< AlgebraicVector > &  subCorrectionTerm, const std::vector< AlgebraicMatrix > &  subDerivatives  )

private

Definition at line 107 of file CompositeAlignmentDerivativesExtractor.cc.

                                                                                                           {
  std::vector<AlgebraicVector>::const_iterator itSubCorrectionTerm = subCorrectionTerm.begin();
  std::vector<AlgebraicMatrix>::const_iterator itSubDerivatives = subDerivatives.begin();
 
  int iRow = 1;
  int iCollumn = 1;
 
  // Fill in the individual terms
  while (itSubCorrectionTerm != subCorrectionTerm.end()) {
    theCorrectionTerm.sub(iRow, *itSubCorrectionTerm);
    theDerivatives.sub(iRow, iCollumn, *itSubDerivatives);
 
    iRow += 2;
    iCollumn += (*itSubDerivatives).num_col();
 
    ++itSubCorrectionTerm;
    ++itSubDerivatives;
  }
 
  return;
}

References theCorrectionTerm, and theDerivatives.

Referenced by extractCurrentAlignment().

Member Data Documentation

theCorrectionTerm

AlgebraicVector CompositeAlignmentDerivativesExtractor::theCorrectionTerm

private

Definition at line 51 of file CompositeAlignmentDerivativesExtractor.h.

Referenced by correctionTerm(), extractCurrentAlignment(), extractWithMultipleHits(), and extractWithoutMultipleHits().

theDerivatives

AlgebraicMatrix CompositeAlignmentDerivativesExtractor::theDerivatives

private

Definition at line 50 of file CompositeAlignmentDerivativesExtractor.h.

Referenced by derivatives(), extractCurrentAlignment(), extractWithMultipleHits(), and extractWithoutMultipleHits().