Class for measurements implementing a space point hit geometry. More...

#include <SpacepointMeasurement.h>

Inheritance diagram for genfit::SpacepointMeasurement:

[legend]

Collaboration diagram for genfit::SpacepointMeasurement:

[legend]

Public Member Functions
	SpacepointMeasurement (int nDim=3)

	SpacepointMeasurement (const TVectorD &rawHitCoords, const TMatrixDSym &rawHitCov, int detId, int hitId, TrackPoint *trackPoint, bool weightedPlaneContruction=true, bool cutCov=true)

virtual	~SpacepointMeasurement ()

virtual AbsMeasurement *	clone () const
	Deep copy ctor for polymorphic class.

virtual SharedPlanePtr	constructPlane (const StateOnPlane &state) const
	Contruct the virtual detector plane.

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

virtual const AbsHMatrix *	constructHMatrix (const AbsTrackRep *) const

Public Member Functions inherited from genfit::AbsMeasurement
	AbsMeasurement ()

	AbsMeasurement (int nDims)

	AbsMeasurement (const TVectorD &rawHitCoords, const TMatrixDSym &rawHitCov, int detId, int hitId, TrackPoint *trackPoint)

virtual	~AbsMeasurement ()

TrackPoint *	getTrackPoint () const

void	setTrackPoint (TrackPoint *tp)

const TVectorD &	getRawHitCoords () const

const TMatrixDSym &	getRawHitCov () const

int	getDetId () const

int	getHitId () const

unsigned int	getDim () const

void	setDetId (int detId)

void	setHitId (int hitId)

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

Private Member Functions
void	initG ()

Private Attributes
bool	weightedPlaneContruction_

TMatrixDSym	G_

bool	cutCov_
	inverse of 3x3 cov

Additional Inherited Members
Protected Member Functions inherited from genfit::AbsMeasurement
	AbsMeasurement (const AbsMeasurement &)
	protect from calling copy c'tor from outside the class. Use clone() if you want a copy!

Protected Attributes inherited from genfit::AbsMeasurement
TVectorD	rawHitCoords_

TMatrixDSym	rawHitCov_

int	detId_

int	hitId_

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

Detailed Description

Class for measurements implementing a space point hit geometry.

Author: Christian Höppner (Technische Universität München, original author); Sebastian Neubert (Technische Universität München, original author); Johannes Rauch (Technische Universität München, original author)

For a space point the detector plane has to be defined with respect to a track representation. SpacepointMeasurement implements a scheme where the detectorplane is chosen perpendicular to the track. In a track fit, only two of the three coordinates of a space point are independent (the track is a one-dimensional object). Therefore the 3D data of the hit is used to define a proper detector plane into which the hit coordinates are then projected.

Definition at line 46 of file SpacepointMeasurement.h.

Constructor & Destructor Documentation

◆ SpacepointMeasurement() [1/2]

genfit::SpacepointMeasurement::SpacepointMeasurement ( int nDim = 3 )

Definition at line 33 of file SpacepointMeasurement.cc.

  : AbsMeasurement(nDim), weightedPlaneContruction_(true), G_(3), cutCov_(true)
{
  assert(nDim >= 3);
 
  //this->initG();
}

◆ SpacepointMeasurement() [2/2]

genfit::SpacepointMeasurement::SpacepointMeasurement	(	const TVectorD &	rawHitCoords,
		const TMatrixDSym &	rawHitCov,
		int	detId,
		int	hitId,
		TrackPoint *	trackPoint,
		bool	weightedPlaneContruction = `true`,
		bool	cutCov = `true`
	)

Definition at line 41 of file SpacepointMeasurement.cc.

  : AbsMeasurement(rawHitCoords, rawHitCov, detId, hitId, trackPoint),
    weightedPlaneContruction_(weightedPlaneContruction), cutCov_(cutCov)
{
  assert(rawHitCoords_.GetNrows() >= 3);
 
  if (weightedPlaneContruction_)
    this->initG();
}

◆ ~SpacepointMeasurement()

virtual genfit::SpacepointMeasurement::~SpacepointMeasurement ( )

inlinevirtual

Definition at line 53 of file SpacepointMeasurement.h.

53{;}

Member Function Documentation

◆ clone()

virtual AbsMeasurement * genfit::SpacepointMeasurement::clone ( ) const

inlinevirtual

Deep copy ctor for polymorphic class.

Implements genfit::AbsMeasurement.

Reimplemented in genfit::ProlateSpacepointMeasurement, and genfit::mySpacepointMeasurement.

Definition at line 55 of file SpacepointMeasurement.h.

55{return new SpacepointMeasurement(*this);}

genfit::SpacepointMeasurement::SpacepointMeasurement

SpacepointMeasurement(int nDim=3)

Definition SpacepointMeasurement.cc:33

◆ constructHMatrix()

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

virtual

Returns a new AbsHMatrix object. Caller must take ownership.

Implements genfit::AbsMeasurement.

Definition at line 120 of file SpacepointMeasurement.cc.

                                                                                      {
  if (dynamic_cast<const RKTrackRep*>(rep) == NULL) {
    Exception exc("SpacepointMeasurement default implementation can only handle state vectors of type RKTrackRep!", __LINE__,__FILE__);
    throw exc;
  }
 
  return new HMatrixUV();
}

◆ constructMeasurementsOnPlane()

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

virtual

Construct MeasurementOnPlane on plane of the state and wrt the states TrackRep. The state will usually be the prediction or reference state, and has to be defined AT the measurement. The AbsMeasurement will be projected onto the plane. It's possible to make corrections to the coordinates here (e.g. by using the state coordinates). Usually the vector will contain only one element. But in the case of e.g. a WireMeasurement, it will be 2 (left and right).

Implements genfit::AbsMeasurement.

Definition at line 69 of file SpacepointMeasurement.cc.

{
  MeasurementOnPlane* mop = new MeasurementOnPlane(TVectorD(2),
       TMatrixDSym(3), // will be resized to 2x2 by Similarity later
       state.getPlane(), state.getRep(), constructHMatrix(state.getRep()));
 
  TVectorD& m = mop->getState();
  TMatrixDSym& V = mop->getCov();
 
  const TVector3& o(state.getPlane()->getO());
  const TVector3& u(state.getPlane()->getU());
  const TVector3& v(state.getPlane()->getV());
 
  // m
  m(0) = (rawHitCoords_(0)-o.X()) * u.X() +
         (rawHitCoords_(1)-o.Y()) * u.Y() +
         (rawHitCoords_(2)-o.Z()) * u.Z();
 
  m(1) = (rawHitCoords_(0)-o.X()) * v.X() +
         (rawHitCoords_(1)-o.Y()) * v.Y() +
         (rawHitCoords_(2)-o.Z()) * v.Z();
 
  //
  // V
  //
  V = rawHitCov_;
 
  // jac = dF_i/dx_j = s_unitvec * t_untivec, with s=u,v and t=x,y,z
  TMatrixD jac(3,2);
  jac(0,0) = u.X();
  jac(1,0) = u.Y();
  jac(2,0) = u.Z();
  jac(0,1) = v.X();
  jac(1,1) = v.Y();
  jac(2,1) = v.Z();
 
  if (cutCov_) { // cut (default)
    tools::invertMatrix(V);
    V.SimilarityT(jac);
    tools::invertMatrix(V);
  }
  else { // projection
    V.SimilarityT(jac);
  }
 
  std::vector<MeasurementOnPlane*> retVal;
  retVal.push_back(mop);
  return retVal;
}

◆ constructPlane()

SharedPlanePtr genfit::SpacepointMeasurement::constructPlane ( const StateOnPlane & state ) const

virtual

Contruct the virtual detector plane.

Per default, the plane will be constructed such that it contains the measurement and POCA to the measurement in cartesian space. The plane is perpendicular to the track (at the POCA).

If weightedPlaneContruction_ is set, the POCA will be calculated in a space weighted with the inverse of the 3D covariance. E.g. if the covariance is very oblate, the plane will be almost defined by the covariance shape. If the covariance is very prolate, the behaviour will be very similar to the ProlateSpacepointHit.

Implements genfit::AbsMeasurement.

Reimplemented in genfit::ProlateSpacepointMeasurement.

Definition at line 53 of file SpacepointMeasurement.cc.

                                                                                    {
 
  // copy state. Neglect covariance.
  StateOnPlane st(state);
 
  const TVector3 point(rawHitCoords_(0), rawHitCoords_(1), rawHitCoords_(2));
 
  if (weightedPlaneContruction_)
    st.extrapolateToPoint(point, G_);
  else
    st.extrapolateToPoint(point);
 
  return st.getPlane();
}

◆ initG()

void genfit::SpacepointMeasurement::initG ( )

private

Definition at line 130 of file SpacepointMeasurement.cc.

                                  {
  rawHitCov_.GetSub(0, 2, G_);
  tools::invertMatrix(G_);
}

Member Data Documentation

◆ cutCov_

bool genfit::SpacepointMeasurement::cutCov_

private

inverse of 3x3 cov

Definition at line 79 of file SpacepointMeasurement.h.

◆ G_

TMatrixDSym genfit::SpacepointMeasurement::G_

private

Definition at line 78 of file SpacepointMeasurement.h.

◆ weightedPlaneContruction_

bool genfit::SpacepointMeasurement::weightedPlaneContruction_

private

Definition at line 77 of file SpacepointMeasurement.h.

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

genfit/measurements/include/SpacepointMeasurement.h
genfit/measurements/src/SpacepointMeasurement.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ SpacepointMeasurement() [1/2]

◆ SpacepointMeasurement() [2/2]

◆ ~SpacepointMeasurement()

Member Function Documentation

◆ clone()

◆ constructHMatrix()

◆ constructMeasurementsOnPlane()

◆ constructPlane()

◆ initG()

Member Data Documentation

◆ cutCov_

◆ G_

◆ weightedPlaneContruction_