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

#include <AlignableSVDRecoHit.h>

Inheritance diagram for AlignableSVDRecoHit:

[legend]

Collaboration diagram for AlignableSVDRecoHit:

[legend]

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. More...

	SVDRecoHit ()
	Inherit constructors.

	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. More...

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

int	getPlaneId () const

virtual SharedPlanePtr	constructPlane (const StateOnPlane &state) const override
	Construct (virtual) detector plane (use state's AbsTrackRep). More...

virtual void	setPlane (const SharedPlanePtr &physicalPlane, int planeId=-1)

void	setStripV (bool v=true)
	Use if the coordinate for 1D hits measured in V direction. More...

TrackPoint *	getTrackPoint () const

void	setTrackPoint (TrackPoint *tp)

const TVectorD &	getRawHitCoords () const

TVectorD &	getRawHitCoords ()

const TMatrixDSym &	getRawHitCov () const

TMatrixDSym &	getRawHitCov ()

int	getDetId () const

int	getHitId () const

virtual bool	isLeftRightMeasurement () const
	If the AbsMeasurement is a wire hit, the left/right resolution will be used.

virtual int	getLeftRightResolution () const

unsigned int	getDim () const

void	setRawHitCoords (const TVectorD &coords)

void	setRawHitCov (const TMatrixDSym &cov)

void	setDetId (int detId)

void	setHitId (int hitId)

virtual void	Print (const Option_t *="") const

virtual std::vector< int >	labels ()
	Vector of integer labels for calibration/alignment parameters available (must match #columns of derivatives(...)) More...

virtual TMatrixD	derivatives (const genfit::StateOnPlane *)
	Derivatives of residuals (local measurement coordinates) w.r.t. More...

virtual TMatrixD	localDerivatives (const genfit::StateOnPlane *)
	Derivatives for additional local parameters to be fitted in global calibration algorithms together with with global parameters. More...

virtual std::vector< int >	localLabels ()
	Vector of integer labels for local calibration parameters available (must match #columns of localDerivatives(...)) More...

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

Protected Attributes
SharedPlanePtr	physicalPlane_

int	planeId_
	This is persistent, but '!' makes ROOT shut up.

bool	stripV_

TVectorD	rawHitCoords_

TMatrixDSym	rawHitCov_

int	detId_

int	hitId_

TrackPoint *	trackPoint_
	Pointer to TrackPoint where the measurement belongs to.

Private Types
enum	{ HIT_DIMENSIONS = 1 }

Private Member Functions
	ClassDefOverride (AlignableSVDRecoHit, 2)
	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
unsigned short	m_sensorID
	Unique sensor identifier.

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

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

bool	m_isU
	transient member (not written out during streaming) More...

float	m_energyDep
	deposited energy.

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

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 };

Belle2::SVDRecoHit::HIT_DIMENSIONS

@ HIT_DIMENSIONS

sensitive Dimensions of the Hit

Definition: SVDRecoHit.h:115

Member Function Documentation

◆ constructHMatrix()

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

inlineoverridevirtualinherited

Methods that actually interface to Genfit.

Reimplemented from PlanarMeasurement.

Definition at line 111 of file SVDRecoHit.h.

◆ constructMeasurementsOnPlane()

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

overridevirtualinherited

Methods that actually interface to Genfit.

Reimplemented from PlanarMeasurement.

Definition at line 180 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::get(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::get(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())));
 }

◆ constructPlane()

SharedPlanePtr constructPlane ( const StateOnPlane & state ) const

overridevirtualinherited

Construct (virtual) detector plane (use state's AbsTrackRep).

It's possible to make corrections to the plane here. The state should be defined somewhere near the measurement. For virtual planes, the state will be extrapolated to the POCA to point (SpacepointMeasurement) or line (WireMeasurement), and from this info the plane will be constructed.

Implements AbsMeasurement.

Definition at line 46 of file PlanarMeasurement.cc.

◆ derivatives()

virtual TMatrixD derivatives ( const genfit::StateOnPlane * )

inlinevirtualinherited

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 both forms are allowed:

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,

or :

G = ( d_sensitive/da d_sensitive/db d_sensitive/dc ... ) as matrix with one row.

A possible algorithm using these derivatives should be able to resolve this based on the measurement HMatrix. Measurements with more dimesions (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: TMatrixD Matrix with #rows = dimension of residual, #columns = number of parameters. #columns must match labels().size().

Definition at line 129 of file ICalibrationParametersDerivatives.h.

◆ 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 dimesions (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().

Reimplemented from ICalibrationParametersDerivatives.

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::get(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;
 }

◆ labels()

virtual std::vector<int> labels ( )

inlinevirtualinherited

Vector of integer labels for calibration/alignment parameters available (must match #columns of derivatives(...))

unique across all sub-detectors in calibration

Returns: std::vector< int > Vector of integer labels

Definition at line 90 of file ICalibrationParametersDerivatives.h.

◆ localDerivatives()

virtual TMatrixD localDerivatives ( const genfit::StateOnPlane * )

inlinevirtualinherited

Derivatives for additional local parameters to be fitted in global calibration algorithms together with with global parameters.

Local parameters are not neccesarily identified by label because their number is proportional to number of measurements included in calibration (possibly very huge number!)

Returns: TMatrixD Matrix in form d_residual_i/d_parameter_j

Reimplemented in AlignableCDCRecoHit.

Definition at line 141 of file ICalibrationParametersDerivatives.h.

◆ localLabels()

virtual std::vector<int> localLabels ( )

inlinevirtualinherited

Vector of integer labels for local calibration parameters available (must match #columns of localDerivatives(...))

This will be usually ignored (e.g. does not have to match localDerivatives), but it is a good practice to return vector of zeros of correct size

Returns: std::vector< int > Vector of integer labels

Definition at line 151 of file ICalibrationParametersDerivatives.h.

◆ setStripV()

void setStripV ( bool v = true )

inlineinherited

Use if the coordinate for 1D hits measured in V direction.

Per default for 1D planar hits, the coordinate is measured in U direction. With this function you can set it to be measured in V direction. This affects the outcoe of constructHMatrix().

Definition at line 70 of file PlanarMeasurement.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 121 of file SVDRecoHit.h.

◆ m_isU

bool m_isU

privateinherited

transient member (not written out during streaming)

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

Definition at line 122 of file SVDRecoHit.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

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Friends

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Member Function Documentation

◆ constructHMatrix()

◆ constructMeasurementsOnPlane()

◆ constructPlane()

◆ derivatives()

◆ globalDerivatives()

◆ labels()

◆ localDerivatives()

◆ localLabels()

◆ setStripV()

Member Data Documentation

◆ m_cluster

◆ m_isU