BowedSurfaceAlignmentDerivatives Class Reference

#include <BowedSurfaceAlignmentDerivatives.h>

Public Types
enum	AlignmentParameterName {   dx = 0, dy, dz, dslopeX,   dslopeY, drotZ, dsagittaX, dsagittaXY,   dsagittaY, N_PARAM }

Public Member Functions
AlgebraicMatrix	operator() (const TrajectoryStateOnSurface &tsos, double uWidth, double vLength, bool doSplit=false, double ySplit=0.) const
	Returns 9x2 jacobian matrix. More...

Static Public Member Functions
static double	gammaScale (double width, double splitLength)

Detailed Description

Calculates alignment derivatives for a bowed surface using Legendre polynomials for the surface structure (as studied by Claus Kleinwort), i.e.

• rigid body part (partially from KarimakiAlignmentDerivatives)
• bow in local u, v and mixed term.

If a surface is split into two parts at a given ySplit value, rotation axes are re-centred to that part hit by the track (as predicted by TSOS) and the length of the surface is re-scaled.

by Gero Flucke, October 2010 $Date$ $Revision$ (last update by $Author$)

Definition at line 25 of file BowedSurfaceAlignmentDerivatives.h.

Member Enumeration Documentation

enum BowedSurfaceAlignmentDerivatives::AlignmentParameterName

Enumerator
dx
dy
dz
dslopeX
dslopeY
drotZ
dsagittaX
dsagittaXY
dsagittaY
N_PARAM

Definition at line 29 of file BowedSurfaceAlignmentDerivatives.h.

                              {
    dx = 0, dy, dz,
    dslopeX, // NOTE: slope(u) -> k*tan(beta), 
    dslopeY, //       slope(v) -> k*tan(alpha)
    drotZ,   // rotation around w axis, scaled by gammaScale
    dsagittaX, dsagittaXY, dsagittaY,
    N_PARAM
  };

Member Function Documentation

double BowedSurfaceAlignmentDerivatives::gammaScale ( double  width, double  splitLength  )

static

scale to apply to convert drotZ to karimaki-gamma, depending on module width and length (the latter after splitting!)

Definition at line 82 of file BowedSurfaceAlignmentDerivatives.cc.

Referenced by TwoBowedSurfacesAlignmentParameters::apply(), operator()(), and BowedSurfaceAlignmentParameters::rotation().

{
//   return 0.5 * std::sqrt(width*width + splitLength*splitLength);
//   return 0.5 * (std::fabs(width) + std::fabs(splitLength));
  return 0.5 * (width + splitLength);
}

AlgebraicMatrix BowedSurfaceAlignmentDerivatives::operator()	(	const TrajectoryStateOnSurface &	tsos,
		double	uWidth,
		double	vLength,
		bool	doSplit = false,
		double	ySplit = 0.
	)		const

Returns 9x2 jacobian matrix.

Definition at line 14 of file BowedSurfaceAlignmentDerivatives.cc.

References drotZ, dsagittaX, dsagittaXY, dsagittaY, dslopeX, dslopeY, dx, dy, dz, gammaScale(), TrajectoryStateOnSurface::localParameters(), LocalTrajectoryParameters::mixedFormatVector(), N_PARAM, query::result, and jetcorrextractor::sign().

{ 

  AlgebraicMatrix result(N_PARAM, 2);

  // track parameters on surface:
  const AlgebraicVector5 tsosPar(tsos.localParameters().mixedFormatVector());
  // [1] dxdz : direction tangent in local xz-plane
  // [2] dydz : direction tangent in local yz-plane
  // [3] x    : local x-coordinate
  // [4] y    : local y-coordinate
  double myY = tsosPar[4];
  double myLengthV = vLength;
  if (doSplit) { // re-'calibrate' y length and transform myY to be w.r.t. surface middle
    // Some signs depend on whether we are in surface part below or above ySplit:
    const double sign = (tsosPar[4] < ySplit ? +1. : -1.); 
    const double yMiddle = ySplit * 0.5 - sign * vLength * .25; // middle of surface
    myY = tsosPar[4] - yMiddle;
    myLengthV = vLength * 0.5 + sign * ySplit;
  }

  const AlgebraicMatrix karimaki(KarimakiAlignmentDerivatives()(tsos)); // it's just 6x2...
  // copy u, v, w from Karimaki - they are independent of splitting
  result[dx][0] = karimaki[0][0];
  result[dx][1] = karimaki[0][1];
  result[dy][0] = karimaki[1][0];
  result[dy][1] = karimaki[1][1];
  result[dz][0] = karimaki[2][0];
  result[dz][1] = karimaki[2][1];
  const double aScale = gammaScale(uWidth, myLengthV);
  result[drotZ][0] = myY / aScale; // Since karimaki[5][0] == vx;
  result[drotZ][1] = karimaki[5][1] / aScale;

  double uRel = 2. * tsosPar[3] / uWidth;  // relative u (-1 .. +1)
  double vRel = 2. * myY / myLengthV;       // relative v (-1 .. +1)
  // 'range check':
  const double cutOff = 1.5;
  if (uRel < -cutOff) { uRel = -cutOff; } else if (uRel > cutOff) { uRel = cutOff; }
  if (vRel < -cutOff) { vRel = -cutOff; } else if (vRel > cutOff) { vRel = cutOff; }

  // Legendre polynomials renormalized to LPn(1)-LPn(0)=1 (n>0)
  const double uLP0 = 1.0;
  const double uLP1 = uRel;
  const double uLP2 = uRel * uRel - 1./3.;
  const double vLP0 = 1.0;
  const double vLP1 = vRel;
  const double vLP2 = vRel * vRel - 1./3.;

  // 1st order (slopes, replacing angles beta, alpha)
  result[dslopeX][0] = tsosPar[1] * uLP1 * vLP0;
  result[dslopeX][1] = tsosPar[2] * uLP1 * vLP0;
  result[dslopeY][0] = tsosPar[1] * uLP0 * vLP1;
  result[dslopeY][1] = tsosPar[2] * uLP0 * vLP1;
  
  // 2nd order (sagitta)
  result[dsagittaX ][0] = tsosPar[1] * uLP2 * vLP0;
  result[dsagittaX ][1] = tsosPar[2] * uLP2 * vLP0;
  result[dsagittaXY][0] = tsosPar[1] * uLP1 * vLP1;
  result[dsagittaXY][1] = tsosPar[2] * uLP1 * vLP1;
  result[dsagittaY ][0] = tsosPar[1] * uLP0 * vLP2;
  result[dsagittaY ][1] = tsosPar[2] * uLP0 * vLP2;
   
  return result;
}