CompositeAlignmentDerivativesExtractor Class Reference

#include <CompositeAlignmentDerivativesExtractor.h>

Public Member Functions
	CompositeAlignmentDerivativesExtractor (const align::Alignables &alignables, const std::vector< AlignableDet * > &alignableDets, const std::vector< TrajectoryStateOnSurface > &tsos)
	deprecated constructor for backward compatibility (use mor general AlignableDetOrUnitPtr) More...
	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::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 12 of file CompositeAlignmentDerivativesExtractor.cc.

References extractCurrentAlignment().

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

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

constructor

Definition at line 28 of file CompositeAlignmentDerivativesExtractor.cc.

References extractCurrentAlignment().

{
  extractCurrentAlignment( alignables, alignableDets, tsos );
}

CompositeAlignmentDerivativesExtractor::~CompositeAlignmentDerivativesExtractor ( void  )

inline

destructor

Definition at line 35 of file CompositeAlignmentDerivativesExtractor.h.

{};

Member Function Documentation

const AlgebraicVector& CompositeAlignmentDerivativesExtractor::correctionTerm ( void  ) const

inline

Definition at line 38 of file CompositeAlignmentDerivativesExtractor.h.

References extractCurrentAlignment(), extractWithMultipleHits(), extractWithoutMultipleHits(), and theCorrectionTerm.

Referenced by CompositeAlignmentParameters::correctionTerm().

{ return theCorrectionTerm; }

const AlgebraicMatrix& CompositeAlignmentDerivativesExtractor::derivatives ( void  ) const

inline

Definition at line 37 of file CompositeAlignmentDerivativesExtractor.h.

References theDerivatives.

Referenced by CompositeAlignmentParameters::derivatives().

{ return theDerivatives; }

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

private

Definition at line 38 of file CompositeAlignmentDerivativesExtractor.cc.

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

Referenced by CompositeAlignmentDerivativesExtractor(), and correctionTerm().

{

  // 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;

}

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

private

Definition at line 146 of file CompositeAlignmentDerivativesExtractor.cc.

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

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

{

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

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

private

Definition at line 118 of file CompositeAlignmentDerivativesExtractor.cc.

References extractWithMultipleHits(), theCorrectionTerm, and theDerivatives.

Referenced by correctionTerm(), and extractCurrentAlignment().

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

Member Data Documentation

AlgebraicVector CompositeAlignmentDerivativesExtractor::theCorrectionTerm

private

Definition at line 54 of file CompositeAlignmentDerivativesExtractor.h.

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

AlgebraicMatrix CompositeAlignmentDerivativesExtractor::theDerivatives

private

Definition at line 53 of file CompositeAlignmentDerivativesExtractor.h.

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