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

#include <SVDRecoHit2D.h>

Inheritance diagram for SVDRecoHit2D:

[legend]

Collaboration diagram for SVDRecoHit2D:

[legend]

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

	SVDRecoHit2D (const SVDTrueHit hit, const genfit::TrackCandHit trackCandHit=nullptr, float sigmaU=-1, float sigmaV=-1)
	Construct SVDRecoHit from a SVDTrueHit for Monte Carlo based tracking.

	SVDRecoHit2D (VxdID::baseType vxdid, const double u, const double v, double sigmaU=-1, double sigmaV=-1)
	Construct SVDRecoHit directly from vxdid and coordinates without using a hit object from the storearray in the constructor.

	SVDRecoHit2D (const SVDCluster &uHit, const SVDCluster &vHit)
	Construct SVDRecoHit from a u- and v- SVDCluster.

	SVDRecoHit2D (const SVDRecoHit &uRecoHit, const SVDRecoHit &vRecoHit)
	Construct SVDRecoHit from a u- and v- (1D) SVDRecoHit.

virtual	~SVDRecoHit2D () 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 *	getUCluster () const
	Get pointer to the u cluster used to create this RecoHit.

const SVDCluster *	getVCluster () const
	Get pointer to the u cluster used to create this RecoHit.

float	getU () const
	Get u coordinate.

float	getV () const
	Get v coordinate.

float	getUVariance () const
	Get u coordinate variance.

float	getVVariance () const
	Get v coordinate variance.

float	getUVCov () const
	Get u-v error covariance.

float	getEnergyDep () const
	Get deposited energy.

virtual std::vector< genfit::MeasurementOnPlane * >	constructMeasurementsOnPlane (const genfit::StateOnPlane &state) const override
	Get deposited energy error.

Private Types
enum	{ HIT_DIMENSIONS = 2 }

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_uCluster
	Pointer to mother uCluster.

const SVDCluster *	m_vCluster
	Pointer to mother vCluster.

float	m_energyDep
	deposited energy.

unsigned short	m_sensorID
	Unique sensor identifier.

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 46 of file SVDRecoHit2D.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

private

Enumerator
HIT_DIMENSIONS	sensitive Dimensions of the Hit

Definition at line 133 of file SVDRecoHit2D.h.

133{ HIT_DIMENSIONS = 2 };

Constructor & Destructor Documentation

◆ SVDRecoHit2D() [1/5]

SVDRecoHit2D ( )

Default constructor for ROOT IO.

Definition at line 23 of file SVDRecoHit2D.cc.

                          :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(0), m_uCluster(0), m_vCluster(0),
  m_energyDep(0), m_sensorID(0)//, m_energyDepError(0)
{}

◆ SVDRecoHit2D() [2/5]

SVDRecoHit2D	(	const SVDTrueHit *	hit,
		const genfit::TrackCandHit *	trackCandHit = nullptr,
		float	sigmaU = -1,
		float	sigmaV = -1 )

explicit

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 sigmaU and sigmaV along u and v respectively

If one of the errors is set <0, a default resolution will be assumed for both values by dividing the pixel size by sqrt(12).

Parameters

hit	SVDTrueHit to use as base
trackCandHit
sigmaU	Error of the Hit along u
sigmaV	Error of the Hit along v

Definition at line 28 of file SVDRecoHit2D.cc.

                                                                                                      :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(hit), m_uCluster(0), m_vCluster(0),
  m_energyDep(0), m_sensorID(0)//, m_energyDepError(0)
{
  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 by dividing the pixel size by sqrt(12)
  if (sigmaU < 0 || sigmaV < 0) {
    const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
    sigmaU = geometry.getUPitch(hit->getV()) / sqrt(12);
    sigmaV = geometry.getVPitch(hit->getV()) / sqrt(12);
  }
 
  // Set positions
  rawHitCoords_(0) = gRandom->Gaus(hit->getU(), sigmaU);
  rawHitCoords_(1) = gRandom->Gaus(hit->getV(), sigmaV);
  // Set the error covariance matrix
  rawHitCov_(0, 0) = sigmaU * sigmaU;
  rawHitCov_(0, 1) = 0;
  rawHitCov_(1, 0) = 0;
  rawHitCov_(1, 1) = sigmaV * sigmaV;
  // Set physical parameters
  m_energyDep = hit->getEnergyDep();
  // Setup geometry information
  setDetectorPlane();
}

◆ SVDRecoHit2D() [3/5]

SVDRecoHit2D	(	VxdID::baseType	vxdid,
		const double	u,
		const double	v,
		double	sigmaU = -1,
		double	sigmaV = -1 )

Construct SVDRecoHit directly from vxdid and coordinates without using a hit object from the storearray in the constructor.

Parameters

vxdid	the vxd id so the reco hit knows on which sensor the measurement took place
u	coordinate of the Hit along u
v	coordinate of the Hit along v
sigmaU	Error of the Hit along u
sigmaV	Error of the Hit along v

Definition at line 58 of file SVDRecoHit2D.cc.

                                                                                                           :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(nullptr), m_uCluster(0), m_vCluster(0),
  m_energyDep(0), m_sensorID(vxdid)//, m_energyDepError(0)
{
  //If no error is given, estimate the error by dividing the pixel size by sqrt(12)
  if (sigmaU < 0 || sigmaV < 0) {
    const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
    sigmaU = geometry.getUPitch(v) / sqrt(12);
    sigmaV = geometry.getVPitch(v) / sqrt(12);
  }
  // Set positions
  rawHitCoords_(0) = u;
  rawHitCoords_(1) = v;
  // Set the error covariance matrix
  rawHitCov_(0, 0) = sigmaU * sigmaU;
  rawHitCov_(0, 1) = 0;
  rawHitCov_(1, 0) = 0;
  rawHitCov_(1, 1) = sigmaV * sigmaV;
  // Setup geometry information
  setDetectorPlane();
}

◆ SVDRecoHit2D() [4/5]

SVDRecoHit2D	(	const SVDCluster &	uHit,
		const SVDCluster &	vHit )

Construct SVDRecoHit from a u- and v- SVDCluster.

Parameters

uHit	the u-type SVDCluster
vHit	the v-type SVDCluster

Definition at line 80 of file SVDRecoHit2D.cc.

                                                                        :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(nullptr), m_uCluster(&uHit), m_vCluster(&vHit),
  m_energyDep(0)
{
  if ((uHit.getRawSensorID() != vHit.getRawSensorID()) || !uHit.isUCluster() || vHit.isUCluster())
    B2FATAL("Error in SVDRecoHit2D: Incorrect SVDCluster instances on input!");
 
  m_sensorID = uHit.getRawSensorID();
 
  // Now that we have a v coordinate, we can rescale u.
  const SVD::SensorInfo& info =
    dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
 
  double DeltaU =
    (info.getForwardWidth() - info.getBackwardWidth()) / info.getLength() / info.getWidth(0);
  double scaleFactorU = 1 + DeltaU * vHit.getPosition();
  double tan_phi = DeltaU * uHit.getPosition(); // need u at v=0!
  double one_over_cos_phi_sqr = 1 + tan_phi * tan_phi;
 
  rawHitCoords_[0] = uHit.getPosition() * scaleFactorU;
  rawHitCoords_[1] = vHit.getPosition();
 
  double sigmaU = uHit.getPositionSigma() * scaleFactorU;
  double sigmaU_sq = sigmaU * sigmaU;
  double sigmaV = vHit.getPositionSigma();
  double sigmaV_sq = sigmaV * sigmaV;
 
  m_energyDep = 0.5 * (uHit.getCharge() + vHit.getCharge());
 
  rawHitCov_(0, 0) = sigmaV_sq * tan_phi * tan_phi + sigmaU_sq * one_over_cos_phi_sqr;
  rawHitCov_(0, 1) = sigmaV_sq * tan_phi;
  rawHitCov_(1, 0) = sigmaV_sq * tan_phi;
  rawHitCov_(1, 1) = sigmaV_sq;
  // Setup geometry information
  setDetectorPlane();
}

◆ SVDRecoHit2D() [5/5]

SVDRecoHit2D	(	const SVDRecoHit &	uRecoHit,
		const SVDRecoHit &	vRecoHit )

Construct SVDRecoHit from a u- and v- (1D) SVDRecoHit.

Parameters

uRecoHit	the u-type SVDRecoHit
vRecoHit	the v-type SVDRecoHit

Definition at line 117 of file SVDRecoHit2D.cc.

                                                                                :
  genfit::PlanarMeasurement(HIT_DIMENSIONS), m_trueHit(nullptr), m_uCluster(0), m_vCluster(0), m_energyDep(0)
{
  const SVDCluster& uHit = *(uRecoHit.getCluster());
  const SVDCluster& vHit = *(vRecoHit.getCluster());
  if ((uHit.getRawSensorID() != vHit.getRawSensorID()) || !uHit.isUCluster() || vHit.isUCluster())
    B2FATAL("Error in SVDRecoHit2D: Incorrect SVDCluster instances on input!");
 
  m_sensorID = uHit.getRawSensorID();
  m_uCluster = &uHit;
  m_vCluster = &vHit;
 
  // Now that we have a v coordinate, we can rescale u.
  const SVD::SensorInfo& info =
    dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
 
  double DeltaU =
    (info.getForwardWidth() - info.getBackwardWidth()) / info.getLength() / info.getWidth(0);
  double scaleFactorU = 1 + DeltaU * vHit.getPosition();
  double tan_phi = DeltaU * uHit.getPosition(); // need u at v=0!
  double one_over_cos_phi_sqr = 1 + tan_phi * tan_phi;
 
  rawHitCoords_[0] = uHit.getPosition() * scaleFactorU;
  rawHitCoords_[1] = vHit.getPosition();
 
  double sigmaU = uHit.getPositionSigma() * scaleFactorU;
  double sigmaU_sq = sigmaU * sigmaU;
  double sigmaV = vHit.getPositionSigma();
  double sigmaV_sq = sigmaV * sigmaV;
 
  m_energyDep = 0.5 * (uHit.getCharge() + vHit.getCharge());
 
  rawHitCov_(0, 0) = sigmaV_sq * tan_phi * tan_phi + sigmaU_sq * one_over_cos_phi_sqr;
  rawHitCov_(0, 1) = sigmaV_sq * tan_phi;
  rawHitCov_(1, 0) = sigmaV_sq * tan_phi;
  rawHitCov_(1, 1) = sigmaV_sq;
  // Setup geometry information
  setDetectorPlane();
}

◆ ~SVDRecoHit2D()

virtual ~SVDRecoHit2D ( )

inlineoverridevirtual

Destructor.

Definition at line 92 of file SVDRecoHit2D.h.

92{}

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 173 of file SVDRecoHit2D.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 = rawHit[0];                 // U coordinate of hit
  double v = rawHit[1];                 // V 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)
 
  /* 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(2);
 
  pos[0] = u + dw * du_dw;
  pos[1] = 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 226 of file SVDRecoHit2D.cc.

{
  return new SVDRecoHit2D(*this);
}

◆ constructMeasurementsOnPlane()

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

overridevirtual

Get deposited energy error.

Methods that actually interface to Genfit.

Definition at line 216 of file SVDRecoHit2D.cc.

{
  // Apply planar deformation
  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())));
}

◆ getEnergyDep()

float getEnergyDep ( ) const

inline

Get deposited energy.

Definition at line 122 of file SVDRecoHit2D.h.

122{ return m_energyDep; }

◆ getSensorID()

VxdID getSensorID ( ) const

inline

Get the compact ID.

Definition at line 100 of file SVDRecoHit2D.h.

100{ 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 103 of file SVDRecoHit2D.h.

103{ return m_trueHit; }

◆ getU()

float getU ( ) const

inline

Get u coordinate.

Definition at line 110 of file SVDRecoHit2D.h.

110{ return rawHitCoords_(0); }

◆ getUCluster()

const SVDCluster * getUCluster ( ) const

inline

Get pointer to the u cluster used to create this RecoHit.

Definition at line 105 of file SVDRecoHit2D.h.

105{return m_uCluster; }

◆ getUVariance()

float getUVariance ( ) const

inline

Get u coordinate variance.

Definition at line 115 of file SVDRecoHit2D.h.

115{ return rawHitCov_(0, 0); }

◆ getUVCov()

float getUVCov ( ) const

inline

Get u-v error covariance.

Definition at line 119 of file SVDRecoHit2D.h.

119{ return rawHitCov_(0, 1); }

◆ getV()

float getV ( ) const

inline

Get v coordinate.

Definition at line 112 of file SVDRecoHit2D.h.

112{ return rawHitCoords_(1); }

◆ getVCluster()

const SVDCluster * getVCluster ( ) const

inline

Get pointer to the u cluster used to create this RecoHit.

Definition at line 107 of file SVDRecoHit2D.h.

107{return m_vCluster; }

◆ getVVariance()

float getVVariance ( ) const

inline

Get v coordinate variance.

Definition at line 117 of file SVDRecoHit2D.h.

117{ return rawHitCov_(1, 1); }

◆ setDetectorPlane()

void setDetectorPlane ( )

private

Set up Detector plane information.

Definition at line 157 of file SVDRecoHit2D.cc.

{
  // Construct a finite detector plane and set it.
  const SVD::SensorInfo& geometry = dynamic_cast<const SVD::SensorInfo&>(VXD::GeoCache::getInstance().getSensorInfo(m_sensorID));
 
  // Construct vectors o, u, v
  ROOT::Math::XYZVector origin  = geometry.pointToGlobal(ROOT::Math::XYZVector(0, 0, 0), true);
  ROOT::Math::XYZVector uGlobal = geometry.vectorToGlobal(ROOT::Math::XYZVector(1, 0, 0), true);
  ROOT::Math::XYZVector vGlobal = geometry.vectorToGlobal(ROOT::Math::XYZVector(0, 1, 0), true);
 
  //Construct the detector plane
  genfit::SharedPlanePtr detPlane(new genfit::DetPlane(XYZToTVector(origin), XYZToTVector(uGlobal), XYZToTVector(vGlobal),
                                                       new VXD::SensorPlane(m_sensorID, 20, 20)));
  setPlane(detPlane, m_sensorID);
}

Member Data Documentation

◆ m_energyDep

float m_energyDep

private

deposited energy.

Definition at line 139 of file SVDRecoHit2D.h.

◆ m_sensorID

unsigned short m_sensorID

private

Unique sensor identifier.

Definition at line 141 of file SVDRecoHit2D.h.

◆ m_trueHit

const SVDTrueHit* m_trueHit

private

Pointer to the Truehit used to generate this hit.

Definition at line 135 of file SVDRecoHit2D.h.

◆ m_uCluster

const SVDCluster* m_uCluster

private

Pointer to mother uCluster.

Definition at line 136 of file SVDRecoHit2D.h.

◆ m_vCluster

const SVDCluster* m_vCluster

private

Pointer to mother vCluster.

Definition at line 137 of file SVDRecoHit2D.h.

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

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

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ SVDRecoHit2D() [1/5]

◆ SVDRecoHit2D() [2/5]

◆ SVDRecoHit2D() [3/5]

◆ SVDRecoHit2D() [4/5]

◆ SVDRecoHit2D() [5/5]

◆ ~SVDRecoHit2D()

Member Function Documentation

◆ applyPlanarDeformation()

◆ clone()

◆ constructMeasurementsOnPlane()

◆ getEnergyDep()

◆ getSensorID()

◆ getTrueHit()

◆ getU()

◆ getUCluster()

◆ getUVariance()

◆ getUVCov()

◆ getV()

◆ getVCluster()

◆ getVVariance()

◆ setDetectorPlane()

Member Data Documentation

◆ m_energyDep

◆ m_sensorID

◆ m_trueHit

◆ m_uCluster

◆ m_vCluster