BeamSpotAlignmentParameters Class Reference

#include <BeamSpotAlignmentParameters.h>

Inheritance diagram for BeamSpotAlignmentParameters:

Public Types
enum	AlignmentParameterName {   dx = 0, dy, dxslope, dyslope,   N_PARAM }
	Give parameters a name. More...
Public Types inherited from AlignmentParameters
typedef AlignmentParametersData::DataContainer	DataContainer

Public Member Functions
virtual void	apply ()
	apply parameters to alignable More...
	BeamSpotAlignmentParameters (Alignable *alignable, bool calcMis)
	BeamSpotAlignmentParameters (Alignable *alignable, const AlgebraicVector &parameters, const AlgebraicSymMatrix &covMatrix)
	Constructor for full set of parameters. More...
	BeamSpotAlignmentParameters (Alignable *alignable, const AlgebraicVector &parameters, const AlgebraicSymMatrix &covMatrix, const std::vector< bool > &selection)
	Constructor for selection. More...
virtual BeamSpotAlignmentParameters *	clone (const AlgebraicVector &parameters, const AlgebraicSymMatrix &covMatrix) const
	Clone all parameters (for update of parameters) More...
virtual BeamSpotAlignmentParameters *	cloneFromSelected (const AlgebraicVector &parameters, const AlgebraicSymMatrix &covMatrix) const
	Clone selected parameters (for update of parameters) More...
virtual AlgebraicMatrix	derivatives (const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &) const
	Get all derivatives. More...
AlgebraicVector	globalParameters (void) const
	calculate and return parameters in global frame More...
void	print (void) const
	print parameters to screen More...
AlgebraicVector	rotation (void) const
	Get rotation parameters. More...
virtual AlgebraicMatrix	selectedDerivatives (const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &) const
	Get selected derivatives. More...
AlgebraicVector	translation (void) const
	Get translation parameters. More...
virtual int	type () const
	tell type (AlignmentParametersFactory::ParametersType - but no circular dependency) More...
virtual	~BeamSpotAlignmentParameters ()
	Destructor. More...
Public Member Functions inherited from AlignmentParameters
Alignable *	alignable (void) const
	Get pointer to corresponding alignable. More...
	AlignmentParameters ()
	Default constructor. More...
	AlignmentParameters (Alignable *object, const AlgebraicVector &par, const AlgebraicSymMatrix &cov)
	Constructor from given input. More...
	AlignmentParameters (Alignable *object, const AlgebraicVector &par, const AlgebraicSymMatrix &cov, const std::vector< bool > &sel)
	Constructor including selection of active parameters. More...
	AlignmentParameters (Alignable *object, const DataContainer &data)
	Constructor. More...
const AlgebraicSymMatrix &	covariance (void) const
	Get parameter covariance matrix. More...
virtual unsigned int	hierarchyLevel () const
bool	isValid (void) const
	Get validity flag. More...
int	numSelected (void) const
	Get number of selected parameters. More...
const AlgebraicVector &	parameters (void) const
	Get alignment parameters. More...
AlgebraicSymMatrix	selectedCovariance (void) const
	Get covariance matrix of selected parameters. More...
AlgebraicVector	selectedParameters (void) const
	Get selected parameters. More...
const std::vector< bool > &	selector (void) const
	Get alignment parameter selector vector. More...
void	setUserVariables (AlignmentUserVariables *auv)
	Set pointer to user variables. More...
void	setValid (bool v)
	Set validity flag. More...
int	size (void) const
	Get number of parameters. More...
AlignmentUserVariables *	userVariables (void) const
	Get pointer to user variables. More...
virtual	~AlignmentParameters ()
	Destructor. More...

Static Public Member Functions
static AlgebraicVector	displacementFromAlignable (const Alignable *ali)

Additional Inherited Members
Protected Member Functions inherited from AlignmentParameters
AlgebraicSymMatrix	collapseSymMatrix (const AlgebraicSymMatrix &m, const std::vector< bool > &sel) const
AlgebraicVector	collapseVector (const AlgebraicVector &m, const std::vector< bool > &sel) const
AlgebraicSymMatrix	expandSymMatrix (const AlgebraicSymMatrix &m, const std::vector< bool > &sel) const
AlgebraicVector	expandVector (const AlgebraicVector &m, const std::vector< bool > &sel) const
Protected Attributes inherited from AlignmentParameters
bool	bValid
	True if parameters are valid. More...
Alignable *	theAlignable
DataContainer	theData
AlignmentUserVariables *	theUserVariables

Detailed Description

Concrete class for alignment parameters and associated quantities [derived from AlignmentParameters]. The number of parameters N_PARAM is fixed to 4 (2 translations in x & y, and 2 slopes)

Date:

2008/09/02 15:08:12

Revision:

1.13

(last update by

Author:

flucke

)

Definition at line 21 of file BeamSpotAlignmentParameters.h.

Member Enumeration Documentation

enum BeamSpotAlignmentParameters::AlignmentParameterName

Give parameters a name.

Enumerator
dx
dy
dxslope
dyslope
N_PARAM

Definition at line 27 of file BeamSpotAlignmentParameters.h.

    {
      dx = 0, dy,
      dxslope, dyslope,
      N_PARAM
    };

Constructor & Destructor Documentation

BeamSpotAlignmentParameters::BeamSpotAlignmentParameters	(	Alignable *	alignable,
		bool	calcMis
	)

Constructor with empty parameters/covariance (if calcMis = false) or with parameters (no covariance) created from current (mis-)placement of alignable (if calcMis = true).

Definition at line 21 of file BeamSpotAlignmentParameters.cc.

Referenced by clone(), and cloneFromSelected().

                                                                                     :
  AlignmentParameters(ali, displacementFromAlignable(calcMis ? ali : 0),
              AlgebraicSymMatrix(N_PARAM, 0))
{

}

BeamSpotAlignmentParameters::BeamSpotAlignmentParameters	(	Alignable *	alignable,
		const AlgebraicVector &	parameters,
		const AlgebraicSymMatrix &	covMatrix
	)

Constructor for full set of parameters.

Definition at line 29 of file BeamSpotAlignmentParameters.cc.

References edm::hlt::Exception, and N_PARAM.

                                                                  :
  AlignmentParameters( alignable, parameters, covMatrix )
{  
  if (parameters.num_row() != N_PARAM) {
    throw cms::Exception("BadParameters") << "in BeamSpotAlignmentParameters(): "
                                          << parameters.num_row() << " instead of " << N_PARAM 
                                          << " parameters.";
  }
}

BeamSpotAlignmentParameters::BeamSpotAlignmentParameters	(	Alignable *	alignable,
		const AlgebraicVector &	parameters,
		const AlgebraicSymMatrix &	covMatrix,
		const std::vector< bool > &	selection
	)

Constructor for selection.

Definition at line 42 of file BeamSpotAlignmentParameters.cc.

References edm::hlt::Exception, and N_PARAM.

                                                                :
  AlignmentParameters( alignable, parameters, covMatrix, selection )
{  
  if (parameters.num_row() != N_PARAM) {
    throw cms::Exception("BadParameters") << "in BeamSpotAlignmentParameters(): "
                                          << parameters.num_row() << " instead of " << N_PARAM 
                                          << " parameters.";
  }
}

BeamSpotAlignmentParameters::~BeamSpotAlignmentParameters ( )

virtual

Destructor.

Definition at line 56 of file BeamSpotAlignmentParameters.cc.

{

} 

Member Function Documentation

void BeamSpotAlignmentParameters::apply ( )

virtual

apply parameters to alignable

Implements AlignmentParameters.

Definition at line 169 of file BeamSpotAlignmentParameters.cc.

References AlignmentParameters::alignable(), edm::hlt::Exception, Alignable::move(), align::rectify(), makeMuonMisalignmentScenario::rot, Alignable::rotateInGlobalFrame(), rotation(), edm::shift, Alignable::surface(), AlignableSurface::toGlobal(), align::toMatrix(), and translation().

Referenced by heavyIonTools.ConfigureHeavyIons::__call__(), editorTools.UserCodeTool::__call__(), HiCoreTools.RestrictInputToAOD::__call__(), coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), heavyIonTools.ProductionDefaults::__call__(), editorTools.ChangeSource::__call__(), HiCoreTools.RemoveMCMatching::__call__(), cmsswVersionTools.PickRelValInputFiles::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), jetTools.AddJetCollection::__call__(), trigTools.SwitchOnTrigger::__call__(), heavyIonTools.SelectionDefaults::__call__(), HiCoreTools.RemoveAllPATObjectsBut::__call__(), heavyIonTools.DisbaleMonteCarloDeps::__call__(), runType1CaloMEtUncertainties.RunType1CaloMEtUncertainties::__call__(), HiCoreTools.RemoveSpecificPATObjects::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), trackTools.MakeTrackCandidates::__call__(), runMVAMEtUncertainties.RunMVAMEtUncertainties::__call__(), tauTools.AddTauCollection::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), HiCoreTools.RemoveCleaning::__call__(), runType1PFMEtUncertainties.RunType1PFMEtUncertainties::__call__(), HiCoreTools.AddCleaning::__call__(), runNoPileUpMEtUncertainties.RunNoPileUpMEtUncertainties::__call__(), trigTools.SwitchOnTriggerMatchingStandAlone::__call__(), trigTools.SwitchOnTriggerMatchEmbedding::__call__(), jetTools.SwitchJetCollection::__call__(), jetTools.AddJetID::__call__(), and jetTools.SetTagInfos::__call__().

{
  Alignable *alignable = this->alignable();
  if (!alignable) {
    throw cms::Exception("BadParameters") 
      << "BeamSpotAlignmentParameters::apply: parameters without alignable";
  }
  
  // Translation in local frame
  AlgebraicVector shift = this->translation(); // fixme: should be LocalVector
  
  // Translation local->global
  align::GlobalVector gv(shift[0], shift[1], shift[2]);
  alignable->move( gv );

  // Rotation in local frame
  align::EulerAngles angles = this->rotation();
  // original code:
  //  alignable->rotateInLocalFrame( align::toMatrix(angles) );
  // correct for rounding errors:
  align::RotationType rot = alignable->surface().toGlobal( align::toMatrix(angles) );
  align::rectify(rot);
  alignable->rotateInGlobalFrame(rot);
}

BeamSpotAlignmentParameters * BeamSpotAlignmentParameters::clone ( const AlgebraicVector &  parameters, const AlgebraicSymMatrix &  covMatrix  ) const

virtual

Clone all parameters (for update of parameters)

Implements AlignmentParameters.

Definition at line 63 of file BeamSpotAlignmentParameters.cc.

References AlignmentParameters::alignable(), BeamSpotAlignmentParameters(), AlignmentParameters::isValid(), AlignmentParameters::selector(), AlignmentParameters::setUserVariables(), AlignmentParameters::setValid(), and AlignmentParameters::userVariables().

Referenced by cloneFromSelected().

{
  BeamSpotAlignmentParameters* rbap = 
    new BeamSpotAlignmentParameters( alignable(), parameters, covMatrix, selector());

  if (userVariables()) rbap->setUserVariables(userVariables()->clone());
  rbap->setValid(isValid());

  return rbap;
}

BeamSpotAlignmentParameters * BeamSpotAlignmentParameters::cloneFromSelected ( const AlgebraicVector &  parameters, const AlgebraicSymMatrix &  covMatrix  ) const

virtual

Clone selected parameters (for update of parameters)

Implements AlignmentParameters.

Definition at line 77 of file BeamSpotAlignmentParameters.cc.

References AlignmentParameters::alignable(), BeamSpotAlignmentParameters(), clone(), AlignmentParameters::expandSymMatrix(), AlignmentParameters::expandVector(), AlignmentParameters::isValid(), AlignmentParameters::selector(), AlignmentParameters::setUserVariables(), AlignmentParameters::setValid(), and AlignmentParameters::userVariables().

{
  BeamSpotAlignmentParameters* rbap = 
    new BeamSpotAlignmentParameters(alignable(), expandVector( parameters, selector()),
                    expandSymMatrix(covMatrix, selector()), selector());

  if ( userVariables() ) rbap->setUserVariables(userVariables()->clone());
  rbap->setValid(isValid());
  
  return rbap;
}

AlgebraicMatrix BeamSpotAlignmentParameters::derivatives ( const TrajectoryStateOnSurface &  tsos, const AlignableDetOrUnitPtr &  alidet  ) const

virtual

Get all derivatives.

Implements AlignmentParameters.

Definition at line 92 of file BeamSpotAlignmentParameters.cc.

References AlignmentParameters::alignable(), edm::hlt::Exception, and N_PARAM.

Referenced by selectedDerivatives().

{
  const Alignable *ali = this->alignable(); // Alignable of these parameters

  if (ali == alidet) { // same alignable => same frame
    return BeamSpotAlignmentDerivatives()(tsos);
  } else {
    throw cms::Exception("MisMatch")
      << "BeamSpotAlignmentParameters::derivatives: The hit alignable must match the "
      << "aligned one, i.e. these parameters make only sense for AlignableBeamSpot.\n";
    return AlgebraicMatrix(N_PARAM, 2); // please compiler
  }
}

AlgebraicVector BeamSpotAlignmentParameters::displacementFromAlignable ( const Alignable *  ali )

static

Calculate parameter vector of misplacements (shift+rotation) from alignable. (If ali=0, return empty AlgebraicVector of proper length.)

Definition at line 234 of file BeamSpotAlignmentParameters.cc.

References PV3DBase< T, PVType, FrameType >::basicVector(), Alignable::displacement(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, Alignable::globalRotation(), N_PARAM, Alignable::rotation(), TkRotation< T >::transposed(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

{
  AlgebraicVector displacement(N_PARAM);

  if (ali) {
    const align::RotationType& dR = ali->rotation();
    
    const align::LocalVector shifts( ali->globalRotation() * 
                     ( dR.transposed() * ali->displacement().basicVector() ) );
    
    align::GlobalVector gv(0.0, 0.0, 1.0);
    align::LocalVector lv(dR.transposed() * gv.basicVector());
    
    displacement[0] = shifts.x();
    displacement[1] = shifts.y();
    displacement[2] = lv.x()/lv.z();
    displacement[3] = lv.y()/lv.z();
  }

  return displacement;
}

AlgebraicVector BeamSpotAlignmentParameters::globalParameters

(

void

)

const

calculate and return parameters in global frame

Definition at line 201 of file BeamSpotAlignmentParameters.cc.

References N_PARAM, edm::shift, Alignable::surface(), AlignmentParameters::theAlignable, AlignableSurface::toGlobal(), translation(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{
  AlgebraicVector m_GlobalParameters(N_PARAM, 0);

  const AlgebraicVector shift = translation(); // fixme: should return LocalVector

  const align::GlobalVector dg(shift[0], shift[1], 0);
  
  m_GlobalParameters[0] = dg.x();
  m_GlobalParameters[1] = dg.y();
  
  align::LocalVector lv(0.0, 0.0, 1.0);
  align::GlobalVector gv = theAlignable->surface().toGlobal(lv);
  
  double dxdz = gv.x()/gv.z();
  double dydz = gv.x()/gv.z();
  
  m_GlobalParameters[2] = dxdz;
  m_GlobalParameters[3] = dydz;
  
  return m_GlobalParameters;
}

void BeamSpotAlignmentParameters::print

(

void

)

const

print parameters to screen

Definition at line 225 of file BeamSpotAlignmentParameters.cc.

References gather_cfg::cout, and AlignmentParameters::theData.

{

  std::cout << "Contents of BeamSpotAlignmentParameters:"
            << "\nParameters: " << theData->parameters()
            << "\nCovariance: " << theData->covariance() << std::endl;
}

AlgebraicVector BeamSpotAlignmentParameters::rotation

(

void

)

const

Get rotation parameters.

Definition at line 142 of file BeamSpotAlignmentParameters.cc.

References funct::cos(), makeMuonMisalignmentScenario::rot, funct::sin(), AlignmentParameters::theData, and align::toAngles().

Referenced by apply().

{
  AlgebraicVector rot(3);

  double dxdz = theData->parameters()[2];  
  double dydz = theData->parameters()[3];  
  double angleY = std::atan(dxdz);
  double angleX = -std::atan(dydz);

  align::RotationType rotY( std::cos(angleY),  0., -std::sin(angleY), 
                0.,                1.,  0.,
                std::sin(angleY),  0.,  std::cos(angleY) );

  align::RotationType rotX( 1.,  0.,                0.,
                0.,  std::cos(angleX),  std::sin(angleX),
                0., -std::sin(angleX),  std::cos(angleX) );

  align::EulerAngles angles = align::toAngles(rotY * rotX);

  rot[0] = angles(1);
  rot[1] = angles(2);
  rot[2] = angles(3);

  return rot;
}

AlgebraicMatrix BeamSpotAlignmentParameters::selectedDerivatives ( const TrajectoryStateOnSurface &  tsos, const AlignableDetOrUnitPtr &  alignableDet  ) const

virtual

Get selected derivatives.

Reimplemented from AlignmentParameters.

Definition at line 109 of file BeamSpotAlignmentParameters.cc.

References derivatives(), AlignmentParameters::numSelected(), and AlignmentParameters::selector().

{
  const AlgebraicMatrix dev = this->derivatives( tsos, alignableDet );

  int ncols  = dev.num_col();
  int nrows  = dev.num_row();
  int nsel   = numSelected();

  AlgebraicMatrix seldev( nsel, ncols );

  int ir2=0;
  for ( int irow=0; irow<nrows; ++irow ) {
    if (selector()[irow]) {
      for ( int icol=0; icol<ncols; ++icol ) seldev[ir2][icol] = dev[irow][icol];
      ++ir2;
    }
  }

  return seldev;
}

AlgebraicVector BeamSpotAlignmentParameters::translation

(

void

)

const

Get translation parameters.

Definition at line 132 of file BeamSpotAlignmentParameters.cc.

References i, edm::shift, and AlignmentParameters::theData.

Referenced by apply(), and globalParameters().

{ 
  AlgebraicVector shift(3);
  for ( int i=0;i<2;++i ) shift[i] = theData->parameters()[i];
  shift[2] = 0.0;

  return shift;
}

int BeamSpotAlignmentParameters::type ( ) const

virtual

tell type (AlignmentParametersFactory::ParametersType - but no circular dependency)

Implements AlignmentParameters.

Definition at line 195 of file BeamSpotAlignmentParameters.cc.

References AlignmentParametersFactory::kBeamSpot.

Referenced by cuy.ValElement::__init__(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputCommands(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::outputEventContent(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties().

{
  return AlignmentParametersFactory::kBeamSpot;
}