CDCTrajectory2D Class Reference

Particle trajectory as it is seen in xy projection represented as a circle. More...

#include <CDCTrajectory2D.h>

Public Member Functions
	CDCTrajectory2D ()
	Default constructor for ROOT compatibility.
	CDCTrajectory2D (const UncertainPerigeeCircle &perigeeCircle)
	Constructs a trajectory from a generalized circle.
	CDCTrajectory2D (const Vector2D &localOrigin, const UncertainPerigeeCircle &localPerigeeCircle, double flightTime=NAN)
	Constructs a trajectory from a generalized circle and a start point.
	CDCTrajectory2D (const Vector2D &pos2D, double time, const Vector2D &mom2D, double charge, double bZ)
	Construct a trajectory with given start point, transverse momentum at the start point, the given charge and the magnetic field value in z direction.
	CDCTrajectory2D (const Vector2D &pos2D, double time, const Vector2D &mom2D, double charge)
	Construct a trajectory with given start point, transverse momentum at the start point and given charge.
bool	isFitted () const
	Checks if the circle is already set to a valid value.
void	clear ()
	Clears all information from this trajectory.
void	reverse ()
	Reverses the trajectory in place.
CDCTrajectory2D	reversed () const
	Returns the reverse trajectory as a copy.
std::array< double, 2 >	reconstructBothZ (const CDC::WireLine &wireLine, double distance=0.0, double z=0) const
	Gives the two z positions where the given drift circle on the wire line touches the trajectory.
double	reconstructZ (const CDC::WireLine &wireLine, double distance=0.0, double z=0) const
	Gives the one z positions within the CDC closest to the given z where the given drift circle on the wire line touches the trajectory.
std::array< Vector3D, 2 >	reconstructBoth3D (const CDC::WireLine &wireLine, double distance=0.0, double z=0) const
	Gives the two three dimensional points where the drift circle touches the wire line.
Vector3D	reconstruct3D (const CDC::WireLine &wireLine, double distance=0.0, double z=0) const
	Gives the one three dimensional positions within the CDC closest to the given z where the given drift circle on the wire line touches the trajectory.
Vector2D	getClosest (const Vector2D &point) const
	Calculates the closest approach on the trajectory to the given point.
CDC::ISuperLayer	getNextISuperLayer () const
	Indicates which superlayer the trajectory traverses after the one, where the start point of the trajectory is located.
CDC::ISuperLayer	getPreviousISuperLayer () const
	Indicates which superlayer the trajectory traverses before the one, where the start point of the trajectory is located.
CDC::ISuperLayer	getNextAxialISuperLayer () const
	Indicates which axial superlayer the trajectory traverses after the one, where the start point of the trajectory is located.
CDC::ISuperLayer	getPreviousAxialISuperLayer () const
	Indicates which axial superlayer the trajectory traverses before the one, where the start point of the trajectory is located.
CDC::ISuperLayer	getMaximalISuperLayer () const
	Indicates the maximal superlayer the trajectory traverses.
CDC::ISuperLayer	getStartISuperLayer () const
	Indicates the superlayer the trajectory starts in.
CDC::ISuperLayer	getMinimalISuperLayer () const
	Indicates the minimal superlayer the trajectory traverses.
template<class AHits>
EForwardBackward	isForwardOrBackwardTo (const AHits &hits) const
	Calculates if this trajectory and the hits are coaligned Returns:
template<class AFromHits, class AToHits>
double	getArcLength2DGap (const AFromHits &fromHits, const AToHits &toHits) const
	Calculates the perpendicular travel distance from the last position of the fromHits to the first position of the toHits.
template<class AFromHits, class AToHits>
double	getArcLength2DFrontOffset (const AFromHits &fromHits, const AToHits &toHits) const
	Calculates the perpendicular travel distance from the first position of the fromHits to the first position of the toHits.
template<class AFromHits, class AToHits>
double	getArcLength2DBackOffset (const AFromHits &fromHits, const AToHits &toHits) const
	Calculates the perpendicular travel distance from the last position of the fromHits to the last position of the toHits.
template<class AHits>
double	getTotalArcLength2D (const AHits &hits) const
	Calculates the perpendicular travel distance from the first position of the hits to the last position of the hits.
double	calcArcLength2D (const Vector2D &point) const
	Calculate the travel distance from the start position of the trajectory.
double	calcArcLength2DBetween (const Vector2D &fromPoint, const Vector2D &toPoint) const
	Calculate the travel distance between the two given positions Returns the travel distance on the trajectory from the first given point to the second given point.
double	getArcLength2DPeriod () const
	Getter for the arc length for one round trip around the trajectory.
void	setPosMom2D (const Vector2D &pos2D, const Vector2D &mom2D, double charge)
	Setter for start point and momentum at the start point subjected to the charge sign.
ESign	getChargeSign () const
	Gets the charge sign of the trajectory.
double	getAbsMom2D (double bZ) const
	Get the estimation for the absolute value of the transvers momentum.
double	getAbsMom2D () const
	Get the estimation for the absolute value of the transvers momentum.
Vector2D	getMom2DAtSupport (const double bZ) const
	Get the momentum at the support point of the trajectory.
Vector2D	getMom2DAtSupport () const
	Get the momentum at the support point of the trajectory.
Vector2D	getFlightDirection2D (const Vector2D &point) const
	Get the unit direction of flight at the given point, where arcLength2D = 0.
Vector2D	getFlightDirection2DAtSupport () const
	Get the unit direction of flight at the support point, where arcLength2D = 0.
bool	isMovingOutward () const
	Indicates if the trajectory is moving outwards or inwards (to or away from the origin) from the start point on.
Vector2D	getPos2DAtArcLength2D (double arcLength2D)
	Getter for the position at a given two dimensional arc length.
Vector2D	getSupport () const
	Get the support point of the trajectory in global coordinates.
Vector2D	getGlobalPerigee () const
	Getter for the closest approach on the trajectory to the global origin.
Vector2D	getGlobalCenter () const
	Getter for the center of the trajectory in global coordinates.
Vector2D	getOuterExit (double factor=1) const
	Calculates the point where the trajectory meets the outer wall of the CDC.
Vector2D	getInnerExit () const
	Calculates the point where the trajectory meets the inner wall of the CDC.
Vector2D	getExit () const
	Calculates the point where the trajectory leaves the CDC.
bool	isCurler (double factor=1) const
	Checks if the trajectory leaves the outer radius of the CDC times the given tolerance factor.
bool	isOriginer (double factor=1) const
	Checks if the trajectory intersects with the inner radius of the CDC time the given tolerance factor.
double	getMaximalCylindricalR () const
	Getter for the maximal distance from the origin.
double	getMinimalCylindricalR () const
	Getter for the minimal distance from the origin - same as absolute value of the impact parameter.
double	getGlobalImpact () const
	Getter for the signed impact parameter of the trajectory.
double	getDist2D (const Vector2D &point) const
	Calculates the distance from the point to the trajectory as seen from the xy projection.
ERightLeft	isRightOrLeft (const Vector2D &point) const
	Checks if the given point is to the right or to the left of the trajectory.
double	getCurvature () const
	Getter for the curvature as seen from the xy projection.
double	getLocalCovariance (EPerigeeParameter iRow, EPerigeeParameter iCol) const
	Getter for an individual element of the covariance matrix of the local helix parameters.
double	getLocalVariance (EPerigeeParameter i) const
	Getter for an individual diagonal element of the covariance matrix of the local helix parameters.
PerigeeCircle	getGlobalCircle () const
	Getter for the circle in global coordinates.
void	setGlobalCircle (const UncertainPerigeeCircle &perigeeCircle)
	Setter for the generalized circle that describes the trajectory.
const UncertainPerigeeCircle &	getLocalCircle () const
	Getter for the circle in local coordinates.
double	getPValue () const
	Getter for p-value.
double	getChi2 () const
	Getter for the chi2 value of the circle fit.
void	setChi2 (const double chi2)
	Setter for the chi square value of the circle fit.
size_t	getNDF () const
	Getter for the number of degrees of freedom of the circle fit.
void	setNDF (std::size_t ndf)
	Setter for the number of degrees of freedom of the circle fit.
void	setLocalCircle (const UncertainPerigeeCircle &localPerigeeCircle)
	Setter for the generalized circle that describes the trajectory.
const Vector2D &	getLocalOrigin () const
	Getter for the origin of the local coordinate system.
double	setLocalOrigin (const Vector2D &localOrigin)
	Setter for the origin of the local coordinate system.
double	getFlightTime () const
	Getter for the time when the particle reached the support point position.
void	setFlightTime (double flightTime)
	Setter for the time when the particle reached the support point position.

Private Member Functions
CDC::ISuperLayer	getISuperLayerAfter (CDC::ISuperLayer iSuperLayer, bool movingOutward) const
	Returns which superlayer is traversed after the current one following the trajectory outward or inward as indicated by the boolean input.
CDC::ISuperLayer	getISuperLayerAfterStart (bool movingOutward) const
	Returns which superlayer is traversed after the current one following the trajectory outward or inward as indicated by the boolean input.
CDC::ISuperLayer	getISuperLayerAfterStart (EForwardBackward forwardBackwardInfo) const
	Indicates which superlayer is traversed after the current one following the trajectory forward or backward as indicated by the input.
CDC::ISuperLayer	getAxialISuperLayerAfterStart (EForwardBackward forwardBackwardInfo) const
	Indicates which axial superlayer is traversed after the one, where the start point of the trajectory is located considering if you want to follow the trajectory in the forward or backward direction.

Private Attributes
Vector2D	m_localOrigin
	Memory for local coordinate origin of the circle representing the trajectory in global coordinates.
UncertainPerigeeCircle	m_localPerigeeCircle
	Memory for the generalized circle describing the trajectory in coordinates from the local origin.
double	m_flightTime = NAN
	Memory for the estimation of the time at which the particle arrived at the support point.

Detailed Description

Particle trajectory as it is seen in xy projection represented as a circle.

Definition at line 38 of file CDCTrajectory2D.h.

Constructor & Destructor Documentation

CDCTrajectory2D() [1/5]

CDCTrajectory2D

(

)

Default constructor for ROOT compatibility.

Definition at line 38 of file CDCTrajectory2D.cc.

  : m_localOrigin()
  , m_localPerigeeCircle()
{
}

CDCTrajectory2D() [2/5]

CDCTrajectory2D ( const UncertainPerigeeCircle & perigeeCircle )

explicit

Constructs a trajectory from a generalized circle.

Constructs a trajectory which is described by the given line or circle. The start point is set to the closest approach to the origin.

Definition at line 44 of file CDCTrajectory2D.cc.

  : m_localOrigin(0.0, 0.0)
  , m_localPerigeeCircle(perigeeCircle)
{
}

CDCTrajectory2D() [3/5]

CDCTrajectory2D	(	const Vector2D &	localOrigin,
		const UncertainPerigeeCircle &	localPerigeeCircle,
		double	flightTime = NAN )

Constructs a trajectory from a generalized circle and a start point.

Constructs a trajectory which is described by the given line or circle and starts in the given point. The point is taken to be the closest approach to the circle.

Definition at line 50 of file CDCTrajectory2D.cc.

  : m_localOrigin(localOrigin)
  , m_localPerigeeCircle(localPerigeeCircle)
  , m_flightTime(flightTime)
{
}

CDCTrajectory2D() [4/5]

CDCTrajectory2D	(	const Vector2D &	pos2D,
		double	time,
		const Vector2D &	mom2D,
		double	charge,
		double	bZ )

Construct a trajectory with given start point, transverse momentum at the start point, the given charge and the magnetic field value in z direction.

Definition at line 59 of file CDCTrajectory2D.cc.

  : m_localOrigin(pos2D)
  , m_localPerigeeCircle(CDCBFieldUtil::absMom2DToCurvature(mom2D.norm(), charge, bZ),
                         mom2D.unit(),
                         0.0)
  , m_flightTime(time)
{
}

CDCTrajectory2D() [5/5]

CDCTrajectory2D	(	const Vector2D &	pos2D,
		double	time,
		const Vector2D &	mom2D,
		double	charge )

Construct a trajectory with given start point, transverse momentum at the start point and given charge.

Definition at line 72 of file CDCTrajectory2D.cc.

  : m_localOrigin(pos2D)
  , m_localPerigeeCircle(CDCBFieldUtil::absMom2DToCurvature(mom2D.norm(), charge, pos2D),
                         mom2D.unit(),
                         0.0)
  , m_flightTime(time)
{
}

Member Function Documentation

calcArcLength2D()

double calcArcLength2D ( const Vector2D & point ) const

inline

Calculate the travel distance from the start position of the trajectory.

Returns the travel distance on the trajectory from the start point to
the given point. This is subjected to a discontinuity at the far point
of the circle. Hence the value return is in the range from -pi*radius to pi*radius
If you have a heavily curling track you have care about the feasibility of this
calculation.

Definition at line 262 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->arcLengthBetween(Vector2D(0.0, 0.0), point - getLocalOrigin());
      }

calcArcLength2DBetween()

double calcArcLength2DBetween ( const Vector2D & fromPoint, const Vector2D & toPoint ) const

inline

Calculate the travel distance between the two given positions Returns the travel distance on the trajectory from the first given point to
the second given point.

This is subjected to a discontinuity at the far point
of the circle. Hence the value return is in the range from -pi*radius to pi*radius
If you have a heavily curling track you have care about the feasibility of this
calculation.

Definition at line 275 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->arcLengthBetween(fromPoint - getLocalOrigin(),
                                                  toPoint - getLocalOrigin());
      }

clear()

void clear

(

)

Clears all information from this trajectory.

Definition at line 89 of file CDCTrajectory2D.cc.

{
  m_localOrigin.set(0.0, 0.0);
  m_localPerigeeCircle.invalidate();
  m_flightTime = NAN;
}

getAbsMom2D() [1/2]

double getAbsMom2D

(

)

const

Get the estimation for the absolute value of the transvers momentum.

Definition at line 289 of file CDCTrajectory2D.cc.

{
  Vector2D position = getSupport();
  return CDCBFieldUtil::curvatureToAbsMom2D(getLocalCircle()->curvature(), position);
}

getAbsMom2D() [2/2]

double getAbsMom2D

(

double

bZ

)

const

Get the estimation for the absolute value of the transvers momentum.

Definition at line 284 of file CDCTrajectory2D.cc.

{
  return CDCBFieldUtil::curvatureToAbsMom2D(getLocalCircle()->curvature(), bZ);
}

getArcLength2DBackOffset()

template<class AFromHits, class AToHits>

double getArcLength2DBackOffset ( const AFromHits & fromHits, const AToHits & toHits ) const

inline

Calculates the perpendicular travel distance from the last position of the fromHits to the last position of the toHits.

Definition at line 236 of file CDCTrajectory2D.h.

      {
        const Vector2D& fromRecoPos2D = fromHits.back().getRecoPos2D();
        const Vector2D& toRecoPos2D = toHits.back().getRecoPos2D();
        return calcArcLength2DBetween(fromRecoPos2D, toRecoPos2D);
      }

getArcLength2DFrontOffset()

template<class AFromHits, class AToHits>

double getArcLength2DFrontOffset ( const AFromHits & fromHits, const AToHits & toHits ) const

inline

Calculates the perpendicular travel distance from the first position of the fromHits to the first position of the toHits.

Definition at line 224 of file CDCTrajectory2D.h.

      {
        const Vector2D& fromRecoPos2D = fromHits.front().getRecoPos2D();
        const Vector2D& toRecoPos2D = toHits.front().getRecoPos2D();
        return calcArcLength2DBetween(fromRecoPos2D, toRecoPos2D);
      }

getArcLength2DGap()

template<class AFromHits, class AToHits>

double getArcLength2DGap ( const AFromHits & fromHits, const AToHits & toHits ) const

inline

Calculates the perpendicular travel distance from the last position of the fromHits to the first position of the toHits.

Definition at line 212 of file CDCTrajectory2D.h.

      {
        const Vector2D& fromRecoPos2D = fromHits.back().getRecoPos2D();
        const Vector2D& toRecoPos2D = toHits.front().getRecoPos2D();
        return calcArcLength2DBetween(fromRecoPos2D, toRecoPos2D);
      }

getArcLength2DPeriod()

double getArcLength2DPeriod ( ) const

inline

Getter for the arc length for one round trip around the trajectory.

Definition at line 281 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->arcLengthPeriod();
      }

getAxialISuperLayerAfterStart()

ISuperLayer getAxialISuperLayerAfterStart ( EForwardBackward forwardBackwardInfo ) const

private

Indicates which axial superlayer is traversed after the one, where the start point of the trajectory is located considering if you want to follow the trajectory in the forward or backward direction.

Definition at line 219 of file CDCTrajectory2D.cc.

{
  bool movingOutward = isMovingOutward();
  if (forwardBackwardInfo == EForwardBackward::c_Backward) {
    movingOutward = not movingOutward;
  }
  ISuperLayer startISuperLayer = getStartISuperLayer();
  if (ISuperLayerUtil::isInvalid(startISuperLayer)) return ISuperLayerUtil::c_Invalid;
 
  ISuperLayer nextISuperLayer = getISuperLayerAfter(startISuperLayer, movingOutward);
  if (ISuperLayerUtil::isInvalid(nextISuperLayer)) return ISuperLayerUtil::c_Invalid;
  if (ISuperLayerUtil::isAxial(nextISuperLayer)) return nextISuperLayer;
 
  ISuperLayer iSuperLayerStep = nextISuperLayer - startISuperLayer;
  assert(std::abs(iSuperLayerStep) == 1);
  bool nextMovingOutward = iSuperLayerStep > 0;
  return getISuperLayerAfter(nextISuperLayer, nextMovingOutward);
}

getChargeSign()

ESign getChargeSign

(

)

const

Gets the charge sign of the trajectory.

Definition at line 279 of file CDCTrajectory2D.cc.

{
  return CDCBFieldUtil::ccwInfoToChargeSign(getLocalCircle()->orientation());
}

getChi2()

double getChi2 ( ) const

inline

Getter for the chi2 value of the circle fit.

Definition at line 470 of file CDCTrajectory2D.h.

      {
        return getLocalCircle().chi2();
      }

getClosest()

Vector2D getClosest

(

const Vector2D &

point

)

const

Calculates the closest approach on the trajectory to the given point.

Definition at line 172 of file CDCTrajectory2D.cc.

{
  return getLocalCircle()->closest(point - getLocalOrigin()) + getLocalOrigin();
}

getCurvature()

double getCurvature ( ) const

inline

Getter for the curvature as seen from the xy projection.

Definition at line 424 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->curvature();
      }

getDist2D()

double getDist2D ( const Vector2D & point ) const

inline

Calculates the distance from the point to the trajectory as seen from the xy projection.

Definition at line 411 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->distance(point - getLocalOrigin());
      }

getExit()

Vector2D getExit

(

)

const

Calculates the point where the trajectory leaves the CDC.

This method returns the first point in forward flight direction from the start point of the trajectory where it meets either the radius of the inner most layer or the outer radius of the outer most wall. If the trajectory does not leave the CDC by the inner or outer wall this will return Vector2D(nan,nan).

Definition at line 335 of file CDCTrajectory2D.cc.

{
  const Vector2D outerExit = getOuterExit();
  const Vector2D innerExit = getInnerExit();
  const Vector2D localExit =  getLocalCircle()->chooseNextForwardOf(Vector2D(0, 0),
                              outerExit - getLocalOrigin(),
                              innerExit - getLocalOrigin());
  return localExit + getLocalOrigin();
}

getFlightDirection2D()

Vector2D getFlightDirection2D ( const Vector2D & point ) const

inline

Get the unit direction of flight at the given point, where arcLength2D = 0.

Definition at line 312 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->tangential(point - getLocalOrigin());
      }

getFlightDirection2DAtSupport()

Vector2D getFlightDirection2DAtSupport ( ) const

inline

Get the unit direction of flight at the support point, where arcLength2D = 0.

Definition at line 318 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->tangential();
      }

getFlightTime()

double getFlightTime ( ) const

inline

Getter for the time when the particle reached the support point position.

Definition at line 524 of file CDCTrajectory2D.h.

      {
        return m_flightTime;
      }

getGlobalCenter()

Vector2D getGlobalCenter ( ) const

inline

Getter for the center of the trajectory in global coordinates.

Definition at line 349 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->center() + m_localOrigin;
      }

getGlobalCircle()

PerigeeCircle getGlobalCircle ( ) const

inline

Getter for the circle in global coordinates.

Definition at line 443 of file CDCTrajectory2D.h.

      {
        PerigeeCircle result = getLocalCircle();
        result.passiveMoveBy(-getLocalOrigin());
        return result;
      }

getGlobalImpact()

double getGlobalImpact ( ) const

inline

Getter for the signed impact parameter of the trajectory.

Definition at line 405 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->distance(-m_localOrigin);
      }

getGlobalPerigee()

Vector2D getGlobalPerigee ( ) const

inline

Getter for the closest approach on the trajectory to the global origin.

Definition at line 343 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->closest(-m_localOrigin) + m_localOrigin;
      }

getInnerExit()

Vector2D getInnerExit

(

)

const

Calculates the point where the trajectory meets the inner wall of the CDC.

This method returns the first point in forward flight direction from the start point of the trajectory where it meets the inner radius of the inner most layer. If the trajectory does not meet the CDC by the inner wall this will return Vector2D(nan,nan)

Definition at line 295 of file CDCTrajectory2D.cc.

{
  const CDCWireTopology& topology = CDCWireTopology::getInstance();
  const CDCWireLayer& innerMostLayer = topology.getWireLayers().front();
  double innerCylindricalR = innerMostLayer.getInnerCylindricalR();
 
  const Vector2D support = getSupport();
  const PerigeeCircle globalCircle = getGlobalCircle();
  if (support.cylindricalR() < innerCylindricalR) {
    // If we start within the inner volume of the CDC we want the trajectory to enter the CDC
    // and not stop at first intersection with the inner wall.
    // Therefore we take the inner exit that comes after the apogee (far point of the circle).
    const Vector2D apogee = globalCircle.apogee();
    return globalCircle.atCylindricalRForwardOf(apogee, innerCylindricalR);
 
  } else {
    return globalCircle.atCylindricalRForwardOf(support, innerCylindricalR);
  }
}

getISuperLayerAfter()

ISuperLayer getISuperLayerAfter ( CDC::ISuperLayer iSuperLayer, bool movingOutward ) const

private

Returns which superlayer is traversed after the current one following the trajectory outward or inward as indicated by the boolean input.

Definition at line 177 of file CDCTrajectory2D.cc.

{
  if (ISuperLayerUtil::isInvalid(iSuperLayer)) return ISuperLayerUtil::c_Invalid;
 
  ISuperLayer minimalISuperLayer = getMinimalISuperLayer();
  ISuperLayer maximalISuperLayer = getMaximalISuperLayer();
  if (minimalISuperLayer == maximalISuperLayer) return ISuperLayerUtil::c_Invalid; // No next super layer to go to
  if (iSuperLayer == minimalISuperLayer) return ISuperLayerUtil::getNextOutwards(iSuperLayer);
  if (iSuperLayer == maximalISuperLayer) return ISuperLayerUtil::getNextInwards(iSuperLayer);
 
  if (movingOutward) {
    return ISuperLayerUtil::getNextOutwards(iSuperLayer);
  } else {
    return ISuperLayerUtil::getNextInwards(iSuperLayer);
  }
}

getISuperLayerAfterStart() [1/2]

ISuperLayer getISuperLayerAfterStart ( bool movingOutward ) const

private

Returns which superlayer is traversed after the current one following the trajectory outward or inward as indicated by the boolean input.

Definition at line 194 of file CDCTrajectory2D.cc.

{
  ISuperLayer iSuperLayer = getStartISuperLayer();
  return getISuperLayerAfter(iSuperLayer, movingOutward);
}

getISuperLayerAfterStart() [2/2]

ISuperLayer getISuperLayerAfterStart ( EForwardBackward forwardBackwardInfo ) const

private

Indicates which superlayer is traversed after the current one following the trajectory forward or backward as indicated by the input.

Definition at line 200 of file CDCTrajectory2D.cc.

{
  bool movingOutward = isMovingOutward();
  if (forwardBackwardInfo == EForwardBackward::c_Backward) {
    movingOutward = not movingOutward;
  }
  return getISuperLayerAfterStart(movingOutward);
}

getLocalCircle()

const UncertainPerigeeCircle & getLocalCircle ( ) const

inline

Getter for the circle in local coordinates.

Definition at line 458 of file CDCTrajectory2D.h.

      {
        return m_localPerigeeCircle;
      }

getLocalCovariance()

double getLocalCovariance ( EPerigeeParameter iRow, EPerigeeParameter iCol ) const

inline

Getter for an individual element of the covariance matrix of the local helix parameters.

Definition at line 430 of file CDCTrajectory2D.h.

      {
        return getLocalCircle().covariance(iRow, iCol);
      }

getLocalOrigin()

const Vector2D & getLocalOrigin ( ) const

inline

Getter for the origin of the local coordinate system.

Definition at line 500 of file CDCTrajectory2D.h.

      {
        return m_localOrigin;
      }

getLocalVariance()

double getLocalVariance ( EPerigeeParameter i ) const

inline

Getter for an individual diagonal element of the covariance matrix of the local helix parameters.

Definition at line 437 of file CDCTrajectory2D.h.

      {
        return getLocalCircle().variance(i);
      }

getMaximalCylindricalR()

double getMaximalCylindricalR ( ) const

inline

Getter for the maximal distance from the origin.

Definition at line 392 of file CDCTrajectory2D.h.

      {
        return std::fabs(getGlobalImpact() + 2 * getLocalCircle()->radius());
      }

getMaximalISuperLayer()

ISuperLayer getMaximalISuperLayer

(

)

const

Indicates the maximal superlayer the trajectory traverses.

Definition at line 248 of file CDCTrajectory2D.cc.

{
  double maximalCylindricalR = getMaximalCylindricalR();
  return CDCWireTopology::getInstance().getISuperLayerAtCylindricalR(maximalCylindricalR);
}

getMinimalCylindricalR()

double getMinimalCylindricalR ( ) const

inline

Getter for the minimal distance from the origin - same as absolute value of the impact parameter.

Definition at line 399 of file CDCTrajectory2D.h.

      {
        return std::fabs(getGlobalImpact());
      }

getMinimalISuperLayer()

ISuperLayer getMinimalISuperLayer

(

)

const

Indicates the minimal superlayer the trajectory traverses.

Definition at line 260 of file CDCTrajectory2D.cc.

{
  double minimalCylindricalR = getMinimalCylindricalR();
  return CDCWireTopology::getInstance().getISuperLayerAtCylindricalR(minimalCylindricalR);
}

getMom2DAtSupport() [1/2]

Vector2D getMom2DAtSupport ( ) const

inline

Get the momentum at the support point of the trajectory.

Definition at line 306 of file CDCTrajectory2D.h.

      {
        return getFlightDirection2DAtSupport() *= getAbsMom2D();
      }

getMom2DAtSupport() [2/2]

Vector2D getMom2DAtSupport ( const double bZ ) const

inline

Get the momentum at the support point of the trajectory.

Definition at line 300 of file CDCTrajectory2D.h.

      {
        return getFlightDirection2DAtSupport() *= getAbsMom2D(bZ);
      }

getNDF()

size_t getNDF ( ) const

inline

Getter for the number of degrees of freedom of the circle fit.

Definition at line 482 of file CDCTrajectory2D.h.

      {
        return getLocalCircle().ndf();
      }

getNextAxialISuperLayer()

ISuperLayer getNextAxialISuperLayer

(

)

const

Indicates which axial superlayer the trajectory traverses after the one, where the start point of the trajectory is located.

Definition at line 238 of file CDCTrajectory2D.cc.

{
  return getAxialISuperLayerAfterStart(EForwardBackward::c_Forward);
}

getNextISuperLayer()

ISuperLayer getNextISuperLayer

(

)

const

Indicates which superlayer the trajectory traverses after the one, where the start point of the trajectory is located.

Definition at line 209 of file CDCTrajectory2D.cc.

{
  return getISuperLayerAfterStart(EForwardBackward::c_Forward);
}

getOuterExit()

Vector2D getOuterExit

(

double

factor = 1

)

const

Calculates the point where the trajectory meets the outer wall of the CDC.

This method returns the first point in forward flight direction from the start point of the trajectory where it meets the outer radius of the outer most layer. If the trajectory does not meet the CDC by the outer wall this will return Vector2D(nan,nan) The factor can be used to virtually resize the CDC.

Definition at line 315 of file CDCTrajectory2D.cc.

{
  const CDCWireTopology& topology = CDCWireTopology::getInstance();
  const CDCWireLayer& outerMostLayer = topology.getWireLayers().back();
  double outerCylindricalR = outerMostLayer.getOuterCylindricalR() * factor;
 
  const Vector2D support = getSupport();
  const PerigeeCircle globalCircle = getGlobalCircle();
  if (support.cylindricalR() > outerCylindricalR) {
    // If we start outside of the volume of the CDC we want the trajectory to enter the CDC
    // and not stop at first intersection with the outer wall.
    // Therefore we take the outer exit that comes after the perigee.
    const Vector2D perigee = globalCircle.perigee();
    return globalCircle.atCylindricalRForwardOf(perigee, outerCylindricalR);
 
  } else {
    return getGlobalCircle().atCylindricalRForwardOf(support, outerCylindricalR);
  }
}

getPos2DAtArcLength2D()

Vector2D getPos2DAtArcLength2D ( double arcLength2D )

inline

Getter for the position at a given two dimensional arc length.

Definition at line 331 of file CDCTrajectory2D.h.

      {
        return getLocalOrigin() + getLocalCircle()->atArcLength(arcLength2D);
      }

getPreviousAxialISuperLayer()

ISuperLayer getPreviousAxialISuperLayer

(

)

const

Indicates which axial superlayer the trajectory traverses before the one, where the start point of the trajectory is located.

Definition at line 243 of file CDCTrajectory2D.cc.

{
  return getAxialISuperLayerAfterStart(EForwardBackward::c_Backward);
}

getPreviousISuperLayer()

ISuperLayer getPreviousISuperLayer

(

)

const

Indicates which superlayer the trajectory traverses before the one, where the start point of the trajectory is located.

Definition at line 214 of file CDCTrajectory2D.cc.

{
  return getISuperLayerAfterStart(EForwardBackward::c_Backward);
}

getPValue()

double getPValue ( ) const

inline

Getter for p-value.

Definition at line 464 of file CDCTrajectory2D.h.

      {
        return TMath::Prob(getChi2(), getNDF());
      }

getStartISuperLayer()

ISuperLayer getStartISuperLayer

(

)

const

Indicates the superlayer the trajectory starts in.

Definition at line 254 of file CDCTrajectory2D.cc.

{
  double startCylindricalR = getLocalOrigin().cylindricalR();
  return CDCWireTopology::getInstance().getISuperLayerAtCylindricalR(startCylindricalR);
}

getSupport()

Vector2D getSupport ( ) const

inline

Get the support point of the trajectory in global coordinates.

Definition at line 337 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->perigee() + getLocalOrigin();
      }

getTotalArcLength2D()

template<class AHits>

double getTotalArcLength2D ( const AHits & hits ) const

inline

Calculates the perpendicular travel distance from the first position of the hits to the last position of the hits.

Definition at line 246 of file CDCTrajectory2D.h.

      {
        Vector2D frontRecoPos2D = hits.front().getRecoPos2D();
        Vector2D backRecoPos2D = hits.back().getRecoPos2D();
        return calcArcLength2DBetween(frontRecoPos2D, backRecoPos2D);
      }

isCurler()

bool isCurler

(

double

factor = 1

)

const

Checks if the trajectory leaves the outer radius of the CDC times the given tolerance factor.

Definition at line 266 of file CDCTrajectory2D.cc.

{
  const CDCWireTopology& topology = CDCWireTopology::getInstance();
  return getMaximalCylindricalR() < factor * topology.getOuterCylindricalR();
}

isFitted()

bool isFitted

(

)

const

Checks if the circle is already set to a valid value.

Definition at line 84 of file CDCTrajectory2D.cc.

{
  return not getLocalCircle()->isInvalid();
}

isForwardOrBackwardTo()

template<class AHits>

EForwardBackward isForwardOrBackwardTo ( const AHits & hits ) const

inline

Calculates if this trajectory and the hits are coaligned Returns:

• EForwardBackward::c_Forward if the last entity lies behind the first.
• EForwardBackward::c_Backward if the last entity lies before the first.

Definition at line 202 of file CDCTrajectory2D.h.

      {
        return static_cast<EForwardBackward>(sign(getTotalArcLength2D(hits)));
      }

isMovingOutward()

bool isMovingOutward ( ) const

inline

Indicates if the trajectory is moving outwards or inwards (to or away from the origin) from the start point on.

Definition at line 325 of file CDCTrajectory2D.h.

      {
        return getFlightDirection2DAtSupport().dot(getSupport()) > 0;
      }

isOriginer()

bool isOriginer

(

double

factor = 1

)

const

Checks if the trajectory intersects with the inner radius of the CDC time the given tolerance factor.

Definition at line 272 of file CDCTrajectory2D.cc.

{
  const CDCWireTopology& topology = CDCWireTopology::getInstance();
  return getMinimalCylindricalR() < factor * topology.getInnerCylindricalR();
}

isRightOrLeft()

ERightLeft isRightOrLeft ( const Vector2D & point ) const

inline

Checks if the given point is to the right or to the left of the trajectory.

Definition at line 417 of file CDCTrajectory2D.h.

      {
        return getLocalCircle()->isRightOrLeft(point - getLocalOrigin());
      }

reconstruct3D()

Vector3D reconstruct3D	(	const CDC::WireLine &	wireLine,
		double	distance = 0.0,
		double	z = 0 ) const

Gives the one three dimensional positions within the CDC closest to the given z where the given drift circle on the wire line touches the trajectory.

This method makes the reconstruction of the z coordinate possible by using the skewness
stereo layer of the stereo wires. The point is determined such that it is at the (signed) distance to the wire line.

Parameters

wireLine	The geometrical wire line on which the hit is located-
distance	The desired distance from the wire line a.k.a. drift length
z	The expected value of z to which to closest solution should be selected.

Definition at line 163 of file CDCTrajectory2D.cc.

{
  const double recoZ = reconstructZ(wireLine, distance, z);
  const Vector3D recoWirePos2D = wireLine.sagPos3DAtZ(recoZ);
  return Vector3D(getClosest(recoWirePos2D.xy()), recoZ);
}

reconstructBoth3D()

std::array< Vector3D, 2 > reconstructBoth3D	(	const CDC::WireLine &	wireLine,
		double	distance = 0.0,
		double	z = 0 ) const

Gives the two three dimensional points where the drift circle touches the wire line.

Only works for the skew stereo wires.

Parameters

wireLine	The geometrical wire line on which the hit is located-
distance	The desired distance from the wire line a.k.a. drift length
z	The expected value of z to which to closest solution should be selected.

Definition at line 150 of file CDCTrajectory2D.cc.

{
  const std::array<double, 2> solutionsZ = reconstructBothZ(wireLine, distance, z);
 
  const Vector3D firstRecoWirePos3D = wireLine.sagPos3DAtZ(solutionsZ[0]);
  const Vector3D secondRecoWirePos3D = wireLine.sagPos3DAtZ(solutionsZ[1]);
  return {{{getClosest(firstRecoWirePos3D.xy()), firstRecoWirePos3D.z()},
      {getClosest(secondRecoWirePos3D.xy()), secondRecoWirePos3D.z()}
    }};
}

reconstructBothZ()

std::array< double, 2 > reconstructBothZ	(	const CDC::WireLine &	wireLine,
		double	distance = 0.0,
		double	z = 0 ) const

Gives the two z positions where the given drift circle on the wire line touches the trajectory.

Only works for the skew stereo wires

Parameters

wireLine	The geometrical wire line on which the hit is located-
distance	The desired distance from the wire line a.k.a. drift length
z	The expected value of z to which to closest solution should be selected.

Definition at line 110 of file CDCTrajectory2D.cc.

{
  Vector2D globalPos2D = wireLine.sagPos2DAtZ(z);
  Vector2D movePerZ = wireLine.sagMovePerZ(z);
 
  Vector2D localPos2D = globalPos2D - getLocalOrigin();
  const PerigeeCircle& localCircle = getLocalCircle();
 
  double fastDistance = distance != 0.0 ? localCircle.fastDistance(distance) : 0.0;
 
  double c = localCircle.fastDistance(localPos2D) - fastDistance;
  double b = localCircle.gradient(localPos2D).dot(movePerZ);
  double a = localCircle.n3() * movePerZ.normSquared();
 
  const std::pair<double, double> solutionsDeltaZ = solveQuadraticABC(a, b, c);
 
  // Put the solution of smaller deviation first
  const std::array<double, 2> solutionsZ{solutionsDeltaZ.second + z, solutionsDeltaZ.first + z};
  return solutionsZ;
}

reconstructZ()

double reconstructZ	(	const CDC::WireLine &	wireLine,
		double	distance = 0.0,
		double	z = 0 ) const

Gives the one z positions within the CDC closest to the given z where the given drift circle on the wire line touches the trajectory.

Only works for the skew stereo wires.

Parameters

wireLine	The geometrical wire line on which the hit is located-
distance	The desired distance from the wire line a.k.a. drift length
z	The expected value of z to which to closest solution should be selected.

Definition at line 133 of file CDCTrajectory2D.cc.

{
  const std::array<double, 2> solutionsZ = reconstructBothZ(wireLine, distance, z);
 
  bool firstIsInCDC = (wireLine.backwardZ() < solutionsZ[0] and
                       solutionsZ[0] < wireLine.forwardZ());
  bool secondIsInCDC = (wireLine.backwardZ() < solutionsZ[1] and
                        solutionsZ[1] < wireLine.forwardZ());
 
  // Prefer the solution with the smaller deviation from the given z position which is the first
  assert(not(std::fabs(solutionsZ[0] - z) > std::fabs(solutionsZ[1] - z)));
  const double recoZ = (firstIsInCDC or not secondIsInCDC) ? solutionsZ[0] : solutionsZ[1];
  return recoZ;
}

reverse()

void reverse

(

)

Reverses the trajectory in place.

Definition at line 96 of file CDCTrajectory2D.cc.

{
  m_localPerigeeCircle.reverse();
  m_flightTime = -m_flightTime;
}

reversed()

CDCTrajectory2D reversed

(

)

const

Returns the reverse trajectory as a copy.

Definition at line 103 of file CDCTrajectory2D.cc.

{
  CDCTrajectory2D result = *this;
  result.reverse();
  return result;
}

setChi2()

void setChi2 ( const double chi2 )

inline

Setter for the chi square value of the circle fit.

Definition at line 476 of file CDCTrajectory2D.h.

      {
        return m_localPerigeeCircle.setChi2(chi2);
      }

setFlightTime()

void setFlightTime ( double flightTime )

inline

Setter for the time when the particle reached the support point position.

Definition at line 530 of file CDCTrajectory2D.h.

      {
        m_flightTime = flightTime;
      }

setGlobalCircle()

void setGlobalCircle ( const UncertainPerigeeCircle & perigeeCircle )

inline

Setter for the generalized circle that describes the trajectory.

Definition at line 451 of file CDCTrajectory2D.h.

      {
        m_localPerigeeCircle = perigeeCircle;
        m_localPerigeeCircle.passiveMoveBy(getLocalOrigin());
      }

setLocalCircle()

void setLocalCircle ( const UncertainPerigeeCircle & localPerigeeCircle )

inline

Setter for the generalized circle that describes the trajectory.

Definition at line 494 of file CDCTrajectory2D.h.

      {
        m_localPerigeeCircle = localPerigeeCircle;
      }

setLocalOrigin()

double setLocalOrigin

(

const Vector2D &

localOrigin

)

Setter for the origin of the local coordinate system.

This sets the origin point the local circle representation is subjected. The local circle is also changed such that the set of points in global space is not changed on repositioning the local parameterisation. It also implicitly sets up a new reference position for all travel distances.
To be able to keep track how the travel distances have to be shifted by this change
the setter returns the value by which the coordinate s parameter was moved ( passively ).
The relation "old traveldistance - return value == new traveldistance" holds.
(if they are not to far away from the reference points, up to the discontinuity at the
far point on the circle)

Parameters

localOrigin

New local reference point in the global coordinate system

Returns

Travel distance from the old to the new origin point

Definition at line 356 of file CDCTrajectory2D.cc.

{
  double arcLength2D = calcArcLength2D(localOrigin);
  m_flightTime += arcLength2D / Const::speedOfLight;
  m_localPerigeeCircle.passiveMoveBy(localOrigin - m_localOrigin);
  m_localOrigin = localOrigin;
  return arcLength2D;
}

setNDF()

void setNDF ( std::size_t ndf )

inline

Setter for the number of degrees of freedom of the circle fit.

Definition at line 488 of file CDCTrajectory2D.h.

      {
        return m_localPerigeeCircle.setNDF(ndf);
      }

setPosMom2D()

void setPosMom2D	(	const Vector2D &	pos2D,
		const Vector2D &	mom2D,
		double	charge )

Setter for start point and momentum at the start point subjected to the charge sign.

Definition at line 345 of file CDCTrajectory2D.cc.

{
  m_localOrigin = pos2D;
  double curvature = CDCBFieldUtil::absMom2DToCurvature(mom2D.norm(), charge, pos2D);
  Vector2D phiVec = mom2D.unit();
  double impact = 0.0;
  m_localPerigeeCircle = UncertainPerigeeCircle(curvature, phiVec, impact);
}

Member Data Documentation

m_flightTime

double m_flightTime = NAN

private

Memory for the estimation of the time at which the particle arrived at the support point.

Definition at line 543 of file CDCTrajectory2D.h.

m_localOrigin

Vector2D m_localOrigin

private

Memory for local coordinate origin of the circle representing the trajectory in global coordinates.

Definition at line 537 of file CDCTrajectory2D.h.

m_localPerigeeCircle

UncertainPerigeeCircle m_localPerigeeCircle

private

Memory for the generalized circle describing the trajectory in coordinates from the local origin.

Definition at line 540 of file CDCTrajectory2D.h.

---

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

• tracking/trackingUtilities/eventdata/trajectories/include/CDCTrajectory2D.h
• tracking/trackingUtilities/eventdata/trajectories/src/CDCTrajectory2D.cc