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.

References extractCurrentAlignment(), and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

                                                     {
  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() [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.

References extractCurrentAlignment().

                                                     {
  extractCurrentAlignment(alignables, alignableDets, tsos);
}

~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.

References theCorrectionTerm.

Referenced by CompositeAlignmentParameters::correctionTerm().

{ return theCorrectionTerm; }

derivatives()

const AlgebraicMatrix& CompositeAlignmentDerivativesExtractor::derivatives ( void  ) const

inline

Definition at line 35 of file CompositeAlignmentDerivativesExtractor.h.

References theDerivatives.

Referenced by CompositeAlignmentParameters::derivatives().

{ return theDerivatives; }

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.

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

Referenced by CompositeAlignmentDerivativesExtractor().

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

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.

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

Referenced by extractCurrentAlignment().

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

extractWithoutMultipleHits()

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

private

Definition at line 107 of file CompositeAlignmentDerivativesExtractor.cc.

References theCorrectionTerm, and theDerivatives.

Referenced by 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

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().