This class is used to transfer SVD information to the track fit. More...

#include <AlignableSVDRecoHit.h>

Inheritance diagram for AlignableSVDRecoHit:

Public Member Functions
virtual	~AlignableSVDRecoHit ()
	Destructor.

genfit::AbsMeasurement *	clone () const override
	Creating a deep copy of this hit.

virtual std::pair< std::vector< int >, TMatrixD >	globalDerivatives (const genfit::StateOnPlane *sop) override
	Labels and derivatives of residuals (local measurement coordinates) w.r.t.

	SVDRecoHit (const SVDTrueHit *hit, bool uDirection, float sigma=-1)
	Inherit constructors.

	SVDRecoHit (const SVDCluster hit, const genfit::TrackCandHit trackCandHit=nullptr)
	Inherit constructors.

VxdID	getSensorID () const
	Get the compact ID.

const SVDTrueHit *	getTrueHit () const
	Get pointer to the TrueHit used when creating this RecoHit, can be nullptr if created from something else.

const SVDCluster *	getCluster () const
	Get pointer to the Cluster used when creating this RecoHit, can be nullptr if created from something else.

bool	isU () const
	Is the coordinate u or v?

float	getPosition () const
	Get coordinate.

float	getPositionVariance () const
	Get coordinate variance.

float	getEnergyDep () const
	Get deposited energy.

float	getRotation () const
	Get rotation angle.

virtual std::vector< genfit::MeasurementOnPlane * >	constructMeasurementsOnPlane (const genfit::StateOnPlane &state) const override
	Methods that actually interface to Genfit.

virtual const genfit::AbsHMatrix *	constructHMatrix (const genfit::AbsTrackRep *) const override
	Methods that actually interface to Genfit.

Static Public Attributes
static bool	s_enableLorentzGlobalDerivatives = false
	Static enabling(true) or disabling(false) addition of global derivatives for Lorentz shift.

Private Types
enum	{ HIT_DIMENSIONS = 1 }

Private Member Functions
	ClassDefOverride (AlignableSVDRecoHit, 3)
	ROOT Macro.

void	setDetectorPlane ()
	Set up Detector plane information.

TVectorD	applyPlanarDeformation (TVectorD rawHit, std::vector< double > planarParameters, const genfit::StateOnPlane &state) const
	Apply planar deformation of sensors.

Private Attributes
const SVDTrueHit *	m_trueHit
	Pointer to the Truehit used to generate this hit.

const SVDCluster *	m_cluster
	transient member (not written out during streaming)

float	m_energyDep
	transient member (not written out during streaming)

float	m_rotationPhi
	angle of the plane rotation, for u in wedge sensors.

unsigned short	m_sensorID
	Unique sensor identifier.

bool	m_isU
	True if the hit has u-coordinate, false if v.

Friends
class	SVDRecoHit

Detailed Description

This class is used to transfer SVD information to the track fit.

Definition at line 28 of file AlignableSVDRecoHit.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

privateinherited

Enumerator
HIT_DIMENSIONS	sensitive Dimensions of the Hit

Definition at line 115 of file SVDRecoHit.h.

115{ HIT_DIMENSIONS = 1 };

Constructor & Destructor Documentation

◆ ~AlignableSVDRecoHit()

virtual ~AlignableSVDRecoHit ( )

inlinevirtual

Destructor.

Definition at line 38 of file AlignableSVDRecoHit.h.

38{}

Member Function Documentation

◆ applyPlanarDeformation()

TVectorD applyPlanarDeformation	(	TVectorD	rawHit,
		std::vector< double >	planarParameters,
		const genfit::StateOnPlane &	state ) const

privateinherited

Apply planar deformation of sensors.

Definition at line 115 of file SVDRecoHit.cc.

{
  // Legendre parametrization of deformation
  auto L1 = [](double x) {return x;};
  auto L2 = [](double x) {return (3 * x * x - 1) / 2;};
  auto L3 = [](double x) {return (5 * x * x * x - 3 * x) / 2;};
  auto L4 = [](double x) {return (35 * x * x * x * x - 30 * x * x + 3) / 8;};
 
  const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
 
  double u = 0;
  double v = 0;
  double length = 0;
  double width = 0;
 
  if (m_isU) {
    u = rawHit[0];                 // U coordinate of hit
    v = state.getState()(4);       // V coordinate of hit
    width = geometry.getWidth(v);  // Width of sensor (U side) is function of V (slanted)
    length = geometry.getLength(); // Length of sensor (V side)
 
    u = u * 2 / width;             // Legendre parametrization required U in (-1, 1)
    v = v * 2 / length;            // Legendre parametrization required V in (-1, 1)
 
  } else {
    v = rawHit[0];                 // V coordinate of hit
    u = state.getState()(3);       // U coordinate of hit
    length = geometry.getLength(); // Length of sensor (V side) is function of V (slanted)
    width = geometry.getWidth(v);  // Width of sensor (U side)
 
    v = v * 2 / length;            // Legendre parametrization required V in (-1, 1)
    u = u * 2 / width;             // Legendre parametrization required U in (-1, 1)
 
  }
 
  /* Planar deformation using Legendre parametrization
     w(u, v) = L_{31} * L2(u) + L_{32} * L1(u) * L1(v) + L_{33} * L2(v) +
               L_{41} * L3(u) + L_{42} * L2(u) * L1(v) + L_{43} * L1(u) * L2(v) + L_{44} * L3(v) +
               L_{51} * L4(u) + L_{52} * L3(u) * L1(v) + L_{53} * L2(u) * L2(v) + L_{54} * L1(u) * L3(v) + L_{55} * L4(v); */
  double dw =
    planarParameters[0] * L2(u) + planarParameters[1] * L1(u) * L1(v) + planarParameters[2] * L2(v) +
    planarParameters[3] * L3(u) + planarParameters[4] * L2(u) * L1(v) + planarParameters[5] * L1(u) * L2(v) + planarParameters[6] * L3(
      v) +
    planarParameters[7] * L4(u) + planarParameters[8] * L3(u) * L1(v) + planarParameters[9] * L2(u) * L2(v) + planarParameters[10] * L1(
      u) * L3(v) + planarParameters[11] * L4(v);
 
  double du_dw = state.getState()[1]; // slope in U direction
  double dv_dw = state.getState()[2]; // slope in V direction
 
  u = u * width / 2;  // from Legendre to Local parametrization
  v = v * length / 2; // from Legendre to Local parametrization
 
  TVectorD pos(1);
 
  if (m_isU) {
    pos[0] = u + dw * du_dw;
  } else {
    pos[0] = v + dw * dv_dw;
  }
 
  return pos;
}

◆ clone()

genfit::AbsMeasurement * clone ( ) const

inlineoverride

Creating a deep copy of this hit.

Definition at line 41 of file AlignableSVDRecoHit.h.

    {
      return new AlignableSVDRecoHit(*this);
    }

◆ constructHMatrix()

virtual const genfit::AbsHMatrix * constructHMatrix ( const genfit::AbsTrackRep * ) const

inlineoverridevirtualinherited

Methods that actually interface to Genfit.

Definition at line 111 of file SVDRecoHit.h.

111{ if (m_isU) return new genfit::HMatrixU(); else return new genfit::HMatrixV(); }

◆ constructMeasurementsOnPlane()

std::vector< genfit::MeasurementOnPlane * > constructMeasurementsOnPlane ( const genfit::StateOnPlane & state ) const

overridevirtualinherited

Methods that actually interface to Genfit.

Definition at line 179 of file SVDRecoHit.cc.

{
  if (!m_isU || m_rotationPhi == 0.0) {
 
    // Apply planar deformation to rectangular sensor or V coordinate of slanted sensor
    TVectorD pos = applyPlanarDeformation(rawHitCoords_, VXD::GeoCache::getInstance().getSensorInfo(m_sensorID).getSurfaceParameters(),
                                          state);
 
    return std::vector<genfit::MeasurementOnPlane*>(1, new genfit::MeasurementOnPlane(pos, rawHitCov_, state.getPlane(),
                                                    state.getRep(), this->constructHMatrix(state.getRep())));
  }
 
  // Wedged sensor: the measured coordinate in U depends on V and the
  // rotation angle.  Namely, it needs to be scaled.
  double u = rawHitCoords_(0);
  double v = state.getState()(4);
  double uPrime = u - v * tan(m_rotationPhi);
  double scale = uPrime / u;
 
  TVectorD coords(1);
  coords(0) = uPrime;
 
  // Apply planar deformation to U coordinate of slanted sensor
  TVectorD pos = applyPlanarDeformation(coords, VXD::GeoCache::getInstance().getSensorInfo(m_sensorID).getSurfaceParameters(), state);
 
  TMatrixDSym cov(scale * scale * rawHitCov_);
 
  return std::vector<genfit::MeasurementOnPlane*>(1, new genfit::MeasurementOnPlane(pos, cov, state.getPlane(), state.getRep(),
                                                  this->constructHMatrix(state.getRep())));
}

◆ getCluster()

const SVDCluster * getCluster ( ) const

inlineinherited

Get pointer to the Cluster used when creating this RecoHit, can be nullptr if created from something else.

Definition at line 88 of file SVDRecoHit.h.

88{ return m_cluster; }

◆ getEnergyDep()

float getEnergyDep ( ) const

inlineinherited

Get deposited energy.

Definition at line 100 of file SVDRecoHit.h.

100{ return m_energyDep; }

◆ getPosition()

float getPosition ( ) const

inlineinherited

Get coordinate.

Definition at line 94 of file SVDRecoHit.h.

94{ return rawHitCoords_(0); }

◆ getPositionVariance()

float getPositionVariance ( ) const

inlineinherited

Get coordinate variance.

Definition at line 97 of file SVDRecoHit.h.

97{ return rawHitCov_(0, 0); }

◆ getRotation()

float getRotation ( ) const

inlineinherited

Get rotation angle.

Definition at line 103 of file SVDRecoHit.h.

103{ return m_rotationPhi; }

◆ getSensorID()

VxdID getSensorID ( ) const

inlineinherited

Get the compact ID.

Definition at line 82 of file SVDRecoHit.h.

82{ return m_sensorID; }

◆ getTrueHit()

const SVDTrueHit * getTrueHit ( ) const

inlineinherited

Get pointer to the TrueHit used when creating this RecoHit, can be nullptr if created from something else.

Definition at line 85 of file SVDRecoHit.h.

85{ return m_trueHit; }

◆ globalDerivatives()

std::pair< std::vector< int >, TMatrixD > globalDerivatives ( const genfit::StateOnPlane * sop )

overridevirtual

Labels and derivatives of residuals (local measurement coordinates) w.r.t.

alignment/calibration parameters Matrix "G" of derivatives valid for given prediction of track state:

G(i, j) = d_residual_i/d_parameter_j

For 2D measurement (u,v):

G = ( du/da du/db du/dc ... ) ( dv/da dv/db dv/dc ... )

for calibration parameters a, b, c.

For 1D measurement:

G = ( 0 0 0 ... ) ( dv/da dv/db dv/dc ... ) for V-strip,

G = ( du/da du/db du/dc ... ) ( 0 0 0 ... ) for U-strip,

Measurements with more dimensions (slopes, curvature) should provide full 4-5Dx(n params) matrix (state as (q/p, u', v', u, v) or (u', v', u, v))

Parameters

sop	Predicted state of the track as linearization point around which derivatives of alignment/calibration parameters shall be computed

Returns: pair<vector<int>, TMatrixD> With matrix with number of rows = dimension of residual, number of columns = number of parameters. number of columns must match vector<int>.size().

Definition at line 25 of file AlignableSVDRecoHit.cc.

{
  auto alignment = GlobalCalibrationManager::getInstance().getAlignmentHierarchy().getGlobalDerivatives<VXDAlignment>(getPlaneId(),
                   sop);
 
  auto globals = GlobalDerivatives(alignment);
 
  if (s_enableLorentzGlobalDerivatives) {
    auto lorentz = GlobalCalibrationManager::getInstance().getLorentzShiftHierarchy().getGlobalDerivatives<VXDAlignment>(getPlaneId(),
                   sop, BFieldManager::getInstance().getField(ROOT::Math::XYZVector(sop->getPos())));
    globals.add(lorentz);
  }
 
  const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(getSensorID()));
 
  // Legendre parametrization of deformation
  auto L1 = [](double x) {return x;};
  auto L2 = [](double x) {return (3 * x * x - 1) / 2;};
  auto L3 = [](double x) {return (5 * x * x * x - 3 * x) / 2;};
  auto L4 = [](double x) {return (35 * x * x * x * x - 30 * x * x + 3) / 8;};
 
  double du_dw = sop->getState()[1]; // slope in U direction
  double dv_dw = sop->getState()[2]; // slope in V direction
  double u = 0.0;
  double v = 0.0;
 
  if (isU()) {
    u = getPosition();              // U coordinate of hit
    v = sop->getState()[4];         // V coordinate of hit
  } else {
    v = getPosition();              // V coordinate of hit
    u = sop->getState()[3];         // U coordinate of hit
  }
 
  double width = geometry.getWidth(v);     // Width of sensor (U side)
  double length = geometry.getLength();    // Length of sensor (V side)
  u = u * 2 / width;                       // Legendre parametrization required U in (-1, 1)
  v = v * 2 / length;                      // Legendre parametrization required V in (-1, 1)
 
  // Add parameters of surface deformation to alignment
  // Numbering of VXD alignment parameters:
  //  -> 0-6:   Rigid body alignment
  //  -> 31-33: First level of surface deformation
  //  -> 41-44: Second level of surface deformation
  //  -> 51-55: Third level of surface deformation
 
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 31), std::vector<double> {L2(u)*du_dw,       L2(u)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 32), std::vector<double> {L1(u)*L1(v)*du_dw, L1(u)*L1(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 33), std::vector<double> {L2(v)*du_dw,       L2(v)*dv_dw});
 
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 41), std::vector<double> {L3(u)*du_dw,       L3(u)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 42), std::vector<double> {L2(u)*L1(v)*du_dw, L2(u)*L1(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 43), std::vector<double> {L1(u)*L2(v)*du_dw, L1(u)*L2(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 44), std::vector<double> {L3(v)*du_dw,       L3(v)*dv_dw});
 
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 51), std::vector<double> {L4(u)*du_dw,       L4(u)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 52), std::vector<double> {L3(u)*L1(v)*du_dw, L3(u)*L1(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 53), std::vector<double> {L2(u)*L2(v)*du_dw, L2(u)*L2(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 54), std::vector<double> {L1(u)*L3(v)*du_dw, L1(u)*L3(v)*dv_dw});
  globals.add(GlobalLabel::construct<VXDAlignment>(getSensorID(), 55), std::vector<double> {L4(v)*du_dw,       L4(v)*dv_dw});
 
  return globals;
}

◆ isU()

bool isU ( ) const

inlineinherited

Is the coordinate u or v?

Definition at line 91 of file SVDRecoHit.h.

91{ return m_isU; }

◆ setDetectorPlane()

void setDetectorPlane ( )

privateinherited

Set up Detector plane information.

Definition at line 89 of file SVDRecoHit.cc.

{
  // Construct a finite detector plane and set it.
  const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
  bool isWedgeU = m_isU && (geometry.getBackwardWidth() > geometry.getForwardWidth());
 
  // Construct vectors o, u, v
  ROOT::Math::XYZVector uLocal(1, 0, 0);
  ROOT::Math::XYZVector vLocal(0, 1, 0);
  ROOT::Math::XYZVector origin  = geometry.pointToGlobal(ROOT::Math::XYZVector(0, 0, 0), true);
  ROOT::Math::XYZVector uGlobal = geometry.vectorToGlobal(uLocal, true);
  ROOT::Math::XYZVector vGlobal = geometry.vectorToGlobal(vLocal, true);
 
  //Construct the detector plane
  VXD::SensorPlane* finitePlane = new VXD::SensorPlane(m_sensorID, 20.0, 20.0);
  if (isWedgeU) finitePlane->setRotation(m_rotationPhi);
  genfit::SharedPlanePtr detPlane(new genfit::DetPlane(XYZToTVector(origin), XYZToTVector(uGlobal), XYZToTVector(vGlobal),
                                                       finitePlane));
  setPlane(detPlane, m_sensorID);
}

◆ SVDRecoHit() [1/2]

SVDRecoHit	(	const SVDCluster *	hit,
		const genfit::TrackCandHit *	trackCandHit = nullptr )

explicit

Inherit constructors.

Definition at line 71 of file SVDRecoHit.cc.

                                                                      :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(0),
  m_cluster(hit), m_energyDep(0), m_rotationPhi(0), m_sensorID(0)
{
  // Set the sensor UID
  m_sensorID = hit->getSensorID();
  m_isU = hit->isUCluster();
 
  setStripV(!m_isU);
 
  // Determine if we have a wedge sensor.
  const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
 
  bool isWedgeU = m_isU && (geometry.getBackwardWidth() > geometry.getForwardWidth());
 
  // Set positions
  rawHitCoords_(0) = hit->getPosition();
  if (isWedgeU) {
    // For u coordinate in a wedge sensor, the position line is not u = const.
    // We have to rotate the coordinate system to achieve this.
    m_rotationPhi = atan2((geometry.getBackwardWidth() - geometry.getForwardWidth()) / geometry.getWidth(0) * hit->getPosition(),
                          geometry.getLength());
  }
  // Set the error covariance matrix (this does not scale with position)
  rawHitCov_(0, 0) = hit->getPositionSigma() * hit->getPositionSigma();
  // Set physical parameters
  m_energyDep = hit->getCharge();
  // Setup geometry information
  setDetectorPlane();
}

◆ SVDRecoHit() [2/2]

SVDRecoHit	(	const SVDTrueHit *	hit,
		bool	uDirection,
		float	sigma = -1 )

Inherit constructors.

Definition at line 63 of file SVDRecoHit.cc.

                                                                         :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(hit),
  m_cluster(0), m_energyDep(0), m_rotationPhi(0), m_sensorID(0), m_isU(uDirection)
{
  setStripV(!m_isU);
 
  // Smear the coordinate when constructing from a TrueHit.
  if (!gRandom) B2FATAL("gRandom not initialized, please set up gRandom first");
 
  // Set the sensor UID
  m_sensorID = hit->getSensorID();
 
  //If no error is given, estimate the error to pitch/sqrt(12)
  if (sigma < 0) {
    const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
    sigma = (m_isU) ? geometry.getUPitch(hit->getV()) / sqrt(12) : geometry.getVPitch() / sqrt(12);
  }
 
  // Set positions
  rawHitCoords_(0) = (m_isU) ? gRandom->Gaus(hit->getU(), sigma) : gRandom->Gaus(hit->getV(), sigma);
  // Set the error covariance matrix
  rawHitCov_(0, 0) = sigma * sigma;
  // Set physical parameters
  m_energyDep = hit->getEnergyDep();
  // Setup geometry information
  setDetectorPlane();
}

Friends And Related Symbol Documentation

◆ SVDRecoHit

friend class SVDRecoHit

friend

Definition at line 29 of file AlignableSVDRecoHit.h.

Member Data Documentation

◆ m_cluster

const SVDCluster* m_cluster

privateinherited

transient member (not written out during streaming)

Pointer to the Cluster used to generate this hit

Definition at line 120 of file SVDRecoHit.h.

◆ m_energyDep

float m_energyDep

privateinherited

transient member (not written out during streaming)

deposited energy.

Definition at line 121 of file SVDRecoHit.h.

◆ m_isU

bool m_isU

privateinherited

True if the hit has u-coordinate, false if v.

Definition at line 125 of file SVDRecoHit.h.

◆ m_rotationPhi

float m_rotationPhi

privateinherited

angle of the plane rotation, for u in wedge sensors.

Definition at line 123 of file SVDRecoHit.h.

◆ m_sensorID

unsigned short m_sensorID

privateinherited

Unique sensor identifier.

Definition at line 124 of file SVDRecoHit.h.

◆ m_trueHit

const SVDTrueHit* m_trueHit

privateinherited

Pointer to the Truehit used to generate this hit.

Definition at line 118 of file SVDRecoHit.h.

◆ s_enableLorentzGlobalDerivatives

bool s_enableLorentzGlobalDerivatives = false

static

Static enabling(true) or disabling(false) addition of global derivatives for Lorentz shift.

Definition at line 32 of file AlignableSVDRecoHit.h.

The documentation for this class was generated from the following files:

alignment/reconstruction/include/AlignableSVDRecoHit.h
alignment/reconstruction/src/AlignableSVDRecoHit.cc

Public Member Functions

Static Public Attributes

Private Types

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ ~AlignableSVDRecoHit()

Member Function Documentation

◆ applyPlanarDeformation()

◆ clone()

◆ constructHMatrix()

◆ constructMeasurementsOnPlane()

◆ getCluster()

◆ getEnergyDep()

◆ getPosition()

◆ getPositionVariance()

◆ getRotation()

◆ getSensorID()

◆ getTrueHit()

◆ globalDerivatives()

◆ isU()

◆ setDetectorPlane()

◆ SVDRecoHit() [1/2]

◆ SVDRecoHit() [2/2]

Friends And Related Symbol Documentation

◆ SVDRecoHit

Member Data Documentation

◆ m_cluster

◆ m_energyDep

◆ m_isU

◆ m_rotationPhi

◆ m_sensorID

◆ m_trueHit

◆ s_enableLorentzGlobalDerivatives