SVDRecoHit - an extended form of SVDHit containing geometry information. More...

#include <SVDRecoHit.h>

Inheritance diagram for SVDRecoHit:

[legend]

Collaboration diagram for SVDRecoHit:

[legend]

Public Member Functions
	SVDRecoHit ()
	Default constructor for ROOT IO.

	SVDRecoHit (const SVDTrueHit *hit, bool uDirection, float sigma=-1)
	Construct SVDRecoHit from a SVDTrueHit for Monte Carlo based tracking.

	SVDRecoHit (const SVDCluster hit, const genfit::TrackCandHit trackCandHit=nullptr)
	Construct SVDRecoHit from a SVD cluster.

virtual	~SVDRecoHit () override
	Destructor.

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

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.

Private Types
enum	{ HIT_DIMENSIONS = 1 }

Private Member Functions
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.

Detailed Description

SVDRecoHit - an extended form of SVDHit containing geometry information.

To create a list of SVDRecoHits for all SVDTrueHits belonging to one MCParticle do something like:

//Get the MCParticle in question
MCParticle* mcParticle = ...
//Iterate over the relation and create a list of hits
vector<SVDRecoHit*> hits;
RelationIndex<MCParticle,SVDTrueHit> relMCTrueHit;
RelationIndex<MCParticle,SVDTrueHit>::range_from it = relMCTrueHit.getElementsFrom(mcParticle);
for(; it.first!=it.second; ++it.first){
  hits.push_back(new SVDRecoHit(it.first->to));
}

Definition at line 47 of file SVDRecoHit.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

private

Enumerator
HIT_DIMENSIONS	sensitive Dimensions of the Hit

Definition at line 115 of file SVDRecoHit.h.

115{ HIT_DIMENSIONS = 1 };

Constructor & Destructor Documentation

◆ SVDRecoHit() [1/3]

SVDRecoHit ( )

Default constructor for ROOT IO.

Definition at line 23 of file SVDRecoHit.cc.

                      :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(0),
  m_cluster(0), m_energyDep(0), m_rotationPhi(0), m_sensorID(0), m_isU(0)
{
  setStripV(!m_isU);
}

◆ SVDRecoHit() [2/3]

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

Construct SVDRecoHit from a SVDTrueHit for Monte Carlo based tracking.

This requires a valid random number generator to be initialized at gRandom. The Hit position will be smeared using a gaussian smearing with sigma along u or v, depending on uDirection.

If the error is set to <0, the resolution will be set to pitch/sqrt(12).

Parameters

hit	SVDTrueHit to use as base
uDirection	If true, the hit has u coordinate, otherwise v
sigma	Error of the hit coordinate

Definition at line 30 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();
}

◆ SVDRecoHit() [3/3]

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

explicit

Construct SVDRecoHit from a SVD cluster.

Parameters

hit	SVDCluster to use as base.
trackCandHit	FIXME: Parameter sigma is no longer used and will be removed.

Definition at line 58 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()

virtual ~SVDRecoHit ( )

inlineoverridevirtual

Destructor.

Definition at line 74 of file SVDRecoHit.h.

74{}

Member Function Documentation

◆ applyPlanarDeformation()

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

private

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

override

Creating a deep copy of this hit.

Overrides the method inherited from GFRecoHit.

Definition at line 110 of file SVDRecoHit.cc.

{
  return new SVDRecoHit(*this);
}

◆ constructHMatrix()

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

inlineoverridevirtual

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

overridevirtual

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

inline

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

inline

Get deposited energy.

Definition at line 100 of file SVDRecoHit.h.

100{ return m_energyDep; }

◆ getPosition()

float getPosition ( ) const

inline

Get coordinate.

Definition at line 94 of file SVDRecoHit.h.

94{ return rawHitCoords_(0); }

◆ getPositionVariance()

float getPositionVariance ( ) const

inline

Get coordinate variance.

Definition at line 97 of file SVDRecoHit.h.

97{ return rawHitCov_(0, 0); }

◆ getRotation()

float getRotation ( ) const

inline

Get rotation angle.

Definition at line 103 of file SVDRecoHit.h.

103{ return m_rotationPhi; }

◆ getSensorID()

VxdID getSensorID ( ) const

inline

Get the compact ID.

Definition at line 82 of file SVDRecoHit.h.

82{ return m_sensorID; }

◆ getTrueHit()

const SVDTrueHit * getTrueHit ( ) const

inline

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

◆ isU()

bool isU ( ) const

inline

Is the coordinate u or v?

Definition at line 91 of file SVDRecoHit.h.

91{ return m_isU; }

◆ setDetectorPlane()

void setDetectorPlane ( )

private

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

Member Data Documentation

◆ m_cluster

const SVDCluster* m_cluster

private

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

private

transient member (not written out during streaming)

deposited energy.

Definition at line 121 of file SVDRecoHit.h.

◆ m_isU

bool m_isU

private

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

Definition at line 125 of file SVDRecoHit.h.

◆ m_rotationPhi

float m_rotationPhi

private

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

Definition at line 123 of file SVDRecoHit.h.

◆ m_sensorID

unsigned short m_sensorID

private

Unique sensor identifier.

Definition at line 124 of file SVDRecoHit.h.

◆ m_trueHit

const SVDTrueHit* m_trueHit

private

Pointer to the Truehit used to generate this hit.

Definition at line 118 of file SVDRecoHit.h.

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

svd/reconstruction/include/SVDRecoHit.h
svd/reconstruction/src/SVDRecoHit.cc

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ SVDRecoHit() [1/3]

◆ SVDRecoHit() [2/3]

◆ SVDRecoHit() [3/3]

◆ ~SVDRecoHit()

Member Function Documentation

◆ applyPlanarDeformation()

◆ clone()

◆ constructHMatrix()

◆ constructMeasurementsOnPlane()

◆ getCluster()

◆ getEnergyDep()

◆ getPosition()

◆ getPositionVariance()

◆ getRotation()

◆ getSensorID()

◆ getTrueHit()

◆ isU()

◆ setDetectorPlane()

Member Data Documentation

◆ m_cluster

◆ m_energyDep

◆ m_isU

◆ m_rotationPhi

◆ m_sensorID

◆ m_trueHit