|
|
|
#include <ParametersToParametersDerivatives.h>
Public Member Functions | |
| bool | isOK () const |
| Indicate whether able to provide the derivatives. | |
| double | operator() (unsigned int indParMother, unsigned int indParComp) const |
| ParametersToParametersDerivatives (const Alignable &component, const Alignable &mother) | |
Private Types | |
| typedef ROOT::Math::SMatrix < double, 6, 9, ROOT::Math::MatRepStd < double, 6, 9 > > | AlgebraicMatrix69 |
Private Member Functions | |
| AlgebraicMatrix69 | dBowed_dRigid (const AlgebraicMatrix66 &f2f, double halfWidth, double halfLength) const |
| bool | init (const Alignable &component, int typeComponent, const Alignable &mother, int typeMother) |
| init by choosing the correct detailed init method depending on parameter types | |
| bool | init2BowedRigid (const Alignable &component, const Alignable &mother) |
| init for component with TwoBowedSurfaces and mother with RigidBody parameters | |
| bool | initBowedRigid (const Alignable &component, const Alignable &mother) |
| init for component with BowedSurface and mother with RigidBody parameters | |
| bool | initRigidRigid (const Alignable &component, const Alignable &mother) |
| init for component and mother both with RigidBody parameters | |
Private Attributes | |
| TMatrixD | derivatives_ |
| can we provide the desired? | |
| bool | isOK_ |
| data members | |
Class for calculating derivatives for hierarchies between different kind of alignment parameters (note that not all combinations might be supported!), needed e.g. to formulate constraints to remove the additional degrees of freedom introduced if larger structure and their components are aligned simultaneously.
(last update by
)
Definition at line 21 of file ParametersToParametersDerivatives.h.
typedef ROOT::Math::SMatrix<double,6,9,ROOT::Math::MatRepStd<double,6,9> > ParametersToParametersDerivatives::AlgebraicMatrix69 [private] |
Definition at line 47 of file ParametersToParametersDerivatives.h.
| ParametersToParametersDerivatives::ParametersToParametersDerivatives | ( | const Alignable & | component, |
| const Alignable & | mother | ||
| ) |
Definition at line 24 of file ParametersToParametersDerivatives.cc.
References Alignable::alignmentParameters(), init, and AlignmentParameters::type().
: isOK_(component.alignmentParameters() && mother.alignmentParameters()) { if (isOK_) { isOK_ = this->init(component, component.alignmentParameters()->type(), mother, mother .alignmentParameters()->type()); } }
| ParametersToParametersDerivatives::AlgebraicMatrix69 ParametersToParametersDerivatives::dBowed_dRigid | ( | const AlgebraicMatrix66 & | f2f, |
| double | halfWidth, | ||
| double | halfLength | ||
| ) | const [private] |
Definition at line 146 of file ParametersToParametersDerivatives.cc.
Referenced by init2BowedRigid(), and initBowedRigid().
{
typedef BowedSurfaceAlignmentDerivatives BowedDerivs;
const double gammaScale = BowedDerivs::gammaScale(2.*halfWidth, 2.*halfLength);
// 1st index (column) is parameter of the mother (<6),
// 2nd index (row) that of component (<9):
AlgebraicMatrix69 derivs;
for (unsigned int iRow = 0; iRow < 6; ++iRow) { // 6 rigid body parameters of mother
// First copy the common rigid body part, e.g.:
// - (0,0): du_comp/du_moth
// - (0,1): dv_comp/du_moth
// - (1,2): dw_comp/dv_moth
for (unsigned int iCol = 0; iCol < 3; ++iCol) { // 3 movements of component
derivs(iRow, iCol) = f2f(iRow, iCol);
}
// Now we have to take care of order and scales for rotation-like parameters:
// slopeX -> halfWidth * beta
derivs(iRow, 3) = halfWidth * f2f(iRow, 4); // = dslopeX_c/dpar_m = hw * db_c/dpar_m
// slopeY -> halfLength * alpha
derivs(iRow, 4) = halfLength * f2f(iRow, 3); // = dslopeY_c/dpar_m = hl * da_c/dpar_m
// rotZ -> gammaScale * gamma
derivs(iRow, 5) = gammaScale * f2f(iRow, 5); // = drotZ_c/dpar_m = gscale * dg_c/dpar_m
// Finally, movements and rotations have no influence on surface internals:
for (unsigned int iCol = 6; iCol < 9; ++iCol) { // 3 sagittae of component
derivs(iRow, iCol) = 0.;
}
}
return derivs;
}
| bool ParametersToParametersDerivatives::init | ( | const Alignable & | component, |
| int | typeComponent, | ||
| const Alignable & | mother, | ||
| int | typeMother | ||
| ) | [private] |
init by choosing the correct detailed init method depending on parameter types
Definition at line 34 of file ParametersToParametersDerivatives.cc.
References init2BowedRigid(), initBowedRigid(), initRigidRigid(), AlignmentParametersFactory::kBowedSurface, AlignmentParametersFactory::kRigidBody, AlignmentParametersFactory::kRigidBody4D, AlignmentParametersFactory::kTwoBowedSurfaces, AlignmentParametersFactory::parametersType(), and AlignmentParametersFactory::parametersTypeName().
{
using namespace AlignmentParametersFactory; // for kRigidBody etc.
if ((typeMother == kRigidBody || typeMother == kRigidBody4D) &&
(typeComponent == kRigidBody || typeComponent == kRigidBody4D)) {
return this->initRigidRigid(component, mother);
} else if ((typeMother == kRigidBody || typeMother == kRigidBody4D) &&
typeComponent == kBowedSurface) {
return this->initBowedRigid(component, mother);
} else if ((typeMother == kRigidBody || typeMother == kRigidBody4D) &&
typeComponent == kTwoBowedSurfaces) {
return this->init2BowedRigid(component, mother);
} else {
// missing: mother with bows and component without, i.e. having 'common' bow parameters
edm::LogError("Alignment") << "@SUB=ParametersToParametersDerivatives::init"
<< "Mother " << parametersTypeName(parametersType(typeMother))
<< ", component " << parametersTypeName(parametersType(typeComponent))
<< ": not supported.";
return false;
}
}
| bool ParametersToParametersDerivatives::init2BowedRigid | ( | const Alignable & | component, |
| const Alignable & | mother | ||
| ) | [private] |
init for component with TwoBowedSurfaces and mother with RigidBody parameters
Definition at line 92 of file ParametersToParametersDerivatives.cc.
References Alignable::alignmentParameters(), dBowed_dRigid(), derivatives_, Alignable::globalPosition(), Alignable::globalRotation(), AlignableSurface::length(), Alignable::surface(), AlignableSurface::toGlobal(), AlignableSurface::width(), and TwoBowedSurfacesAlignmentParameters::ySplit().
Referenced by init().
{
// component is two bowed surfaces, mother rigid body
const TwoBowedSurfacesAlignmentParameters *aliPar =
dynamic_cast<TwoBowedSurfacesAlignmentParameters*>(component.alignmentParameters());
if (!aliPar) {
edm::LogError("Alignment") << "@SUB=ParametersToParametersDerivatives::init2BowedRigid"
<< "dynamic_cast to TwoBowedSurfacesAlignmentParameters failed.";
return false;
}
// We treat the two surfaces as independent objects, i.e.
// 1) get the global position of each surface, depending on the ySplit value,
const double ySplit = aliPar->ySplit();
const double halfWidth = 0.5 * component.surface().width();
const double halfLength = 0.5 * component.surface().length();
const double halfLength1 = 0.5 * (halfLength + ySplit);
const double halfLength2 = 0.5 * (halfLength - ySplit);
const double yM1 = 0.5 * (ySplit - halfLength); // y_mean of surface 1
const double yM2 = yM1 + halfLength; // y_mean of surface 2
// The sensor positions and orientations could be adjusted using
// TwoBowedSurfacesDeformation attached to the component,
// but that should be 2nd order effect.
const align::GlobalPoint posSurf1(component.surface().toGlobal(align::LocalPoint(0.,yM1,0.)));
const align::GlobalPoint posSurf2(component.surface().toGlobal(align::LocalPoint(0.,yM2,0.)));
// 2) get derivatives for both,
FrameToFrameDerivative f2fMaker;
const AlgebraicMatrix66 f2fSurf1(f2fMaker.getDerivative(component.globalRotation(),
mother.globalRotation(),
posSurf1, mother.globalPosition()));
const AlgebraicMatrix66 f2fSurf2(f2fMaker.getDerivative(component.globalRotation(),
mother.globalRotation(),
posSurf2, mother.globalPosition()));
const AlgebraicMatrix69 derivs1(this->dBowed_dRigid(f2fSurf1, halfWidth, halfLength1));
const AlgebraicMatrix69 derivs2(this->dBowed_dRigid(f2fSurf2, halfWidth, halfLength2));
// 3) fill the common matrix by merging the two.
typedef ROOT::Math::SMatrix<double,6,18,ROOT::Math::MatRepStd<double,6,18> > AlgebraicMatrix6_18;
AlgebraicMatrix6_18 derivs;
derivs.Place_at(derivs1, 0, 0); // left half
derivs.Place_at(derivs2, 0, 9); // right half
// copy to TMatrix
derivatives_.ResizeTo(6, 18);
derivatives_.SetMatrixArray(derivs.begin());
return true;
}
| bool ParametersToParametersDerivatives::initBowedRigid | ( | const Alignable & | component, |
| const Alignable & | mother | ||
| ) | [private] |
init for component with BowedSurface and mother with RigidBody parameters
Definition at line 74 of file ParametersToParametersDerivatives.cc.
References dBowed_dRigid(), derivatives_, FrameToFrameDerivative::frameToFrameDerivative(), AlignableSurface::length(), m, Alignable::surface(), and AlignableSurface::width().
Referenced by init().
{
// component is bowed surface, mother rigid body
FrameToFrameDerivative f2fMaker;
const AlgebraicMatrix66 f2f(asSMatrix<6,6>(f2fMaker.frameToFrameDerivative(&component,&mother)));
const double halfWidth = 0.5 * component.surface().width();
const double halfLength = 0.5 * component.surface().length();
const AlgebraicMatrix69 m(this->dBowed_dRigid(f2f, halfWidth, halfLength));
// copy to TMatrix
derivatives_.ResizeTo(6,9);
derivatives_.SetMatrixArray(m.begin());
return true;
}
| bool ParametersToParametersDerivatives::initRigidRigid | ( | const Alignable & | component, |
| const Alignable & | mother | ||
| ) | [private] |
init for component and mother both with RigidBody parameters
Definition at line 59 of file ParametersToParametersDerivatives.cc.
References derivatives_, FrameToFrameDerivative::frameToFrameDerivative(), and m.
Referenced by init().
{
// simply frame to frame!
FrameToFrameDerivative f2fDerivMaker;
AlgebraicMatrix66 m(asSMatrix<6,6>(f2fDerivMaker.frameToFrameDerivative(&component, &mother)));
// copy to TMatrix
derivatives_.ResizeTo(6,6);
derivatives_.SetMatrixArray(m.begin());
return true;
}
| bool ParametersToParametersDerivatives::isOK | ( | ) | const [inline] |
Indicate whether able to provide the derivatives.
Definition at line 27 of file ParametersToParametersDerivatives.h.
References isOK_.
Referenced by AlignmentParameterStore::hierarchyConstraints().
{ return isOK_;}
| double ParametersToParametersDerivatives::operator() | ( | unsigned int | indParMother, |
| unsigned int | indParComp | ||
| ) | const |
Return the derivative DeltaParam(object)/DeltaParam(composedobject), indices start with 0. But check isOK() first!
Definition at line 183 of file ParametersToParametersDerivatives.cc.
References derivatives_.
{
// Do range checks?
return derivatives_(indParMother, indParComp);
}
TMatrixD ParametersToParametersDerivatives::derivatives_ [private] |
can we provide the desired?
Definition at line 53 of file ParametersToParametersDerivatives.h.
Referenced by init2BowedRigid(), initBowedRigid(), initRigidRigid(), and operator()().
bool ParametersToParametersDerivatives::isOK_ [private] |
data members
Definition at line 52 of file ParametersToParametersDerivatives.h.
Referenced by isOK().
1.7.3