Class performs extension (adding new hits) of given candidate using conformal transformation w.r.t point on the trajectory. More...

#include <OffOriginExtension.h>

Inheritance diagram for OffOriginExtension:

Public Member Functions
	OffOriginExtension (std::vector< const TrackingUtilities::CDCWireHit * > allAxialWireHits, double levelPrecision=9)
	Constructor.

void	operator() (const std::vector< const TrackingUtilities::CDCWireHit * > &inputWireHits, void *qt) final
	Main entry point for the post processing call from the QuadTreeProcessor.

std::vector< const TrackingUtilities::CDCWireHit * >	roadSearch (const std::vector< const TrackingUtilities::CDCWireHit * > &wireHits)
	Perform transformation for set of given hits; reference position taken as POCA of the fitted trajectory.

std::vector< const TrackingUtilities::CDCWireHit * >	getHitsWRTtoRefPos (const TrackingUtilities::Vector2D &refPos, float curv, float theta)
	Get hits which are compatible with given trajectory.

const std::vector< TrackingUtilities::CDCTrack > &	getTracks () const
	Get the collected tracks.

Protected Attributes
std::vector< const TrackingUtilities::CDCWireHit * >	m_allAxialWireHits
	Pool of all axial hits from which the road search may select additional hits.

std::vector< TrackingUtilities::CDCTrack >	m_tracks
	Collected tracks.

Private Attributes
double	m_levelPrecision
	Precision level for the width of the off origin hough search.

Detailed Description

Class performs extension (adding new hits) of given candidate using conformal transformation w.r.t point on the trajectory.

Definition at line 28 of file OffOriginExtension.h.

Constructor & Destructor Documentation

◆ OffOriginExtension()

OffOriginExtension	(	std::vector< const TrackingUtilities::CDCWireHit * >	allAxialWireHits,
		double	levelPrecision = 9 )

explicit

Constructor.

Definition at line 27 of file OffOriginExtension.cc.

  : BaseCandidateReceiver(std::move(allAxialWireHits))
  , m_levelPrecision(levelPrecision)
{
}

Member Function Documentation

◆ getHitsWRTtoRefPos()

std::vector< const CDCWireHit * > getHitsWRTtoRefPos	(	const TrackingUtilities::Vector2D &	refPos,
		float	curv,
		float	theta )

Get hits which are compatible with given trajectory.

Parameters

refPos	defines 2D reference position with respect to which transformation will be performed
curv	curvarute of the track trajectory
theta	angle between x-axis and vector to the center of the circle which represents trajectory

Returns: vector of CDCWireHit objects which satisfy legendre transformation with respect to the given parameters

Definition at line 99 of file OffOriginExtension.cc.

{
  float thetaPrecision = M_PI / (pow(2., m_levelPrecision + 1));
  float curvPrecision = 0.15 / (pow(2., m_levelPrecision));
 
  using YSpan = AxialHitQuadTreeProcessor::YSpan;
  YSpan curvSpan{curv - curvPrecision, curv + curvPrecision};
  LookupTable<Vector2D> thetaSpan(&Vector2D::Phi, 1, theta - thetaPrecision, theta + thetaPrecision);
 
  AxialHitQuadTreeProcessor qtProcessor(refPos, curvSpan, &thetaSpan);
  qtProcessor.seed(m_allAxialWireHits);
 
  std::vector<const CDCWireHit*> newWireHits = qtProcessor.getAssignedItems();
  return newWireHits;
}

◆ getTracks()

const std::vector< CDCTrack > & getTracks ( ) const

inherited

Get the collected tracks.

Definition at line 41 of file BaseCandidateReceiver.cc.

{
  return m_tracks;
}

◆ operator()()

void operator()	(	const std::vector< const TrackingUtilities::CDCWireHit * > &	inputWireHits,
		void *	qt )

finalvirtual

Main entry point for the post processing call from the QuadTreeProcessor.

Reimplemented from BaseCandidateReceiver.

Definition at line 34 of file OffOriginExtension.cc.

{
  // Unset the taken flag and let the postprocessing decide
  for (const CDCWireHit* wireHit : inputWireHits) {
    (*wireHit)->setTakenFlag(false);
  }
 
  std::vector<const CDCWireHit*> candidateHits = roadSearch(inputWireHits);
  AxialTrackUtil::addCandidateFromHits(candidateHits,
                                       m_allAxialWireHits,
                                       m_tracks,
                                       true);
}

◆ roadSearch()

std::vector< const CDCWireHit * > roadSearch ( const std::vector< const TrackingUtilities::CDCWireHit * > & wireHits )

Perform transformation for set of given hits; reference position taken as POCA of the fitted trajectory.

Definition at line 50 of file OffOriginExtension.cc.

{
  const CDCKarimakiFitter& fitter = CDCKarimakiFitter::getNoDriftVarianceFitter();
  CDCTrajectory2D trackTrajectory2D = fitter.fit(wireHits);
 
  double chi2 = trackTrajectory2D.getChi2();
  Vector2D refPos = trackTrajectory2D.getGlobalPerigee();
 
  // change sign of the curvature; should be the same as the charge of the candidate
  double curv = trackTrajectory2D.getCurvature();
 
  // theta is the vector normal to the trajectory at the perigee
  double theta = trackTrajectory2D.getGlobalCircle().phi0() + M_PI_2;
 
  // Hide the current hits from the road search
  for (const CDCWireHit* hit : wireHits) {
    hit->getAutomatonCell().setTakenFlag(true);
  }
 
  std::vector<const CDCWireHit*> newWireHits = getHitsWRTtoRefPos(refPos, curv, theta);
 
  // Unhide the current hits from the road search
  for (const CDCWireHit* hit : wireHits) {
    hit->getAutomatonCell().setTakenFlag(false);
  }
 
  if (newWireHits.size() == 0) return wireHits;
 
  std::vector<const CDCWireHit*> combinedWireHits;
 
  for (const CDCWireHit* hit : wireHits) {
    combinedWireHits.push_back(hit);
  }
 
  for (const CDCWireHit* hit : newWireHits) {
    combinedWireHits.push_back(hit);
  }
 
  CDCTrajectory2D combinedTrajectory2D = fitter.fit(wireHits);
  double combinedChi2 = combinedTrajectory2D.getChi2();
 
  if (combinedChi2 < chi2 * 2.) {
    return combinedWireHits;
  }
 
  return wireHits;
}

Member Data Documentation

◆ m_allAxialWireHits

std::vector<const TrackingUtilities::CDCWireHit*> m_allAxialWireHits

protectedinherited

Pool of all axial hits from which the road search may select additional hits.

Definition at line 39 of file BaseCandidateReceiver.h.

◆ m_levelPrecision

double m_levelPrecision

private

Precision level for the width of the off origin hough search.

Definition at line 53 of file OffOriginExtension.h.

◆ m_tracks

std::vector<TrackingUtilities::CDCTrack> m_tracks

protectedinherited

Collected tracks.

Definition at line 42 of file BaseCandidateReceiver.h.

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

tracking/trackFindingCDC/legendre/quadtree/include/OffOriginExtension.h
tracking/trackFindingCDC/legendre/quadtree/src/OffOriginExtension.cc

Public Member Functions

Protected Attributes

Private Attributes