CurrentAlignmentKFUpdator Class Reference

#include <CurrentAlignmentKFUpdator.h>

Inheritance diagram for CurrentAlignmentKFUpdator:

Public Member Functions
virtual CurrentAlignmentKFUpdator *	clone (void) const
	CurrentAlignmentKFUpdator (void)
	CurrentAlignmentKFUpdator (AlignableNavigator *navigator)
template<unsigned int D>
TrajectoryStateOnSurface	update (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const
TrajectoryStateOnSurface	update (const TrajectoryStateOnSurface &, const TransientTrackingRecHit &) const
	~CurrentAlignmentKFUpdator (void)
Public Member Functions inherited from TrajectoryStateUpdator
	TrajectoryStateUpdator ()
virtual	~TrajectoryStateUpdator ()

Private Member Functions
AlignmentParameters *	getAlignmentParameters (const AlignableDetOrUnitPtr alignableDet) const
AlignmentParameters *	getHigherLevelParameters (const Alignable *aAlignable) const
template<unsigned int D>
void	includeCurrentAlignmentEstimate (const TransientTrackingRecHit &aRecHit, const TrajectoryStateOnSurface &tsos, typename AlgebraicROOTObject< D >::Vector &vecR, typename AlgebraicROOTObject< D >::SymMatrix &matV) const

Private Attributes
AlignableNavigator *	theAlignableNavigator

Detailed Description

Definition at line 14 of file CurrentAlignmentKFUpdator.h.

Constructor & Destructor Documentation

CurrentAlignmentKFUpdator::CurrentAlignmentKFUpdator ( void  )

inline

Definition at line 19 of file CurrentAlignmentKFUpdator.h.

Referenced by clone().

: theAlignableNavigator( 0 ) {}

CurrentAlignmentKFUpdator::CurrentAlignmentKFUpdator ( AlignableNavigator *  navigator )

inline

Definition at line 20 of file CurrentAlignmentKFUpdator.h.

: theAlignableNavigator( navigator ) {}

CurrentAlignmentKFUpdator::~CurrentAlignmentKFUpdator ( void  )

inline

Definition at line 21 of file CurrentAlignmentKFUpdator.h.

{}

Member Function Documentation

virtual CurrentAlignmentKFUpdator* CurrentAlignmentKFUpdator::clone ( void  ) const

inlinevirtual

Implements TrajectoryStateUpdator.

Definition at line 28 of file CurrentAlignmentKFUpdator.h.

References CurrentAlignmentKFUpdator().

{ return new CurrentAlignmentKFUpdator( *this ); }

AlignmentParameters * CurrentAlignmentKFUpdator::getAlignmentParameters ( const AlignableDetOrUnitPtr  alignableDet ) const

private

Definition at line 122 of file CurrentAlignmentKFUpdator.cc.

References Alignable::alignmentParameters(), and getHigherLevelParameters().

Referenced by includeCurrentAlignmentEstimate().

{
  // Get alignment parameters from AlignableDet ...
  AlignmentParameters* alignmentParameters = alignableDet->alignmentParameters();
  // ... or any higher level alignable.
  if ( !alignmentParameters ) alignmentParameters = getHigherLevelParameters( alignableDet );
  return alignmentParameters;
}

AlignmentParameters * CurrentAlignmentKFUpdator::getHigherLevelParameters ( const Alignable *  aAlignable ) const

private

Definition at line 132 of file CurrentAlignmentKFUpdator.cc.

References Alignable::alignmentParameters(), and Alignable::mother().

Referenced by getAlignmentParameters().

{
  Alignable* higherLevelAlignable = aAlignable->mother();
  // Alignable has no mother ... most probably the alignable is already the full tracker.
  if ( !higherLevelAlignable ) return 0;
  AlignmentParameters* higherLevelParameters = higherLevelAlignable->alignmentParameters();
  // Found alignment parameters? If not, go one level higher in the hierarchy.
  return higherLevelParameters ? higherLevelParameters : getHigherLevelParameters( higherLevelAlignable );
}

template<unsigned int D>

void CurrentAlignmentKFUpdator::includeCurrentAlignmentEstimate ( const TransientTrackingRecHit &  aRecHit, const TrajectoryStateOnSurface &  tsos, typename AlgebraicROOTObject< D >::Vector &  vecR, typename AlgebraicROOTObject< D >::SymMatrix &  matV  ) const

private

Definition at line 80 of file CurrentAlignmentKFUpdator.cc.

References AlignableNavigator::alignableFromGeomDet(), gather_cfg::cout, deltaR(), TransientTrackingRecHit::det(), getAlignmentParameters(), AlignableDetOrUnitPtr::isNull(), AlignmentParameters::selectedCovariance(), AlignmentParameters::selectedDerivatives(), AlignmentParameters::selectedParameters(), and theAlignableNavigator.

{
  const GeomDet* det = aRecHit.det();
  if ( !det ) return;

  AlignableDetOrUnitPtr alignableDet = theAlignableNavigator->alignableFromGeomDet( det );
  if ( alignableDet.isNull() )
  {
    //std::cout << "[CurrentAlignmentKFUpdator::includeCurrentAlignmentEstimate] No AlignableDet associated with RecHit." << std::endl;
    return;
  }

  AlignmentParameters* alignmentParameters = getAlignmentParameters( alignableDet );

  if ( alignmentParameters )
  {
    AlgebraicMatrix selectedDerivatives = alignmentParameters->selectedDerivatives( tsos, alignableDet );
    AlgebraicVector selectedParameters = alignmentParameters->selectedParameters();
    AlgebraicSymMatrix selectedCovariance = alignmentParameters->selectedCovariance();

    AlgebraicSymMatrix deltaV = selectedCovariance.similarityT( selectedDerivatives );
    AlgebraicVector deltaR = selectedDerivatives.T()*selectedParameters;

    //AlignmentUserVariables* auv = alignmentParameters->userVariables();
    //if ( !auv ) std::cout << "[CurrentAlignmentKFUpdator::includeCurrentAlignmentEstimate] No AlignmentUserVariables associated with AlignableDet." << std::endl;
    //if ( theAnnealing ) matV *= (*theAnnealing)( auv );

    if ( deltaR.num_row() == D )
    {
      vecR += asSVector<D>(deltaR);
      matV += asSMatrix<D>(deltaV);
    }
    else std::cout << "[CurrentAlignmentKFUpdator::includeCurrentAlignmentEstimate] Predicted state and misalignment correction not compatible." << std::endl;
  } else std::cout << "[CurrentAlignmentKFUpdator::includeCurrentAlignmentEstimate] No AlignmentParameters associated with AlignableDet." << std::endl;

  return;
}

template<unsigned int D>

TrajectoryStateOnSurface CurrentAlignmentKFUpdator::update ( const TrajectoryStateOnSurface &  tsos, const TransientTrackingRecHit &  aRecHit  ) const

virtual

Implements TrajectoryStateUpdator.

Definition at line 11 of file CurrentAlignmentKFUpdator.cc.

References TrackingRecHit::dimension(), and edm::hlt::Exception.

Referenced by python.seqvaluedict.seqdict::__setslice__(), python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.seqvaluedict.seqdict::extend(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
    switch (aRecHit.dimension()) {
        case 1: return update<1>(tsos,aRecHit);
        case 2: return update<2>(tsos,aRecHit);
        case 3: return update<3>(tsos,aRecHit);
        case 4: return update<4>(tsos,aRecHit);
        case 5: return update<5>(tsos,aRecHit);
    }
    throw cms::Exception("Rec hit of invalid dimension (not 1,2,3,4,5)");
}

TrajectoryStateOnSurface CurrentAlignmentKFUpdator::update ( const TrajectoryStateOnSurface &  , const TransientTrackingRecHit &    ) const

virtual

Implements TrajectoryStateUpdator.

Referenced by python.seqvaluedict.seqdict::__setslice__(), python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.seqvaluedict.seqdict::extend(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

Member Data Documentation

AlignableNavigator* CurrentAlignmentKFUpdator::theAlignableNavigator

private

Definition at line 41 of file CurrentAlignmentKFUpdator.h.

Referenced by includeCurrentAlignmentEstimate().