SegmentInPhi0CurvBox Class Reference

Predicate class to check for the containment of segments in a phi0 curv hough space part. More...

#include <SegmentInPhi0CurvBox.h>

Private Member Functions
TrackingUtilities::Weight	operator() (const TrackingUtilities::CDCSegment2D *segment, const Phi0CurvBox *phi0CurvBox)
	Checks if the segment is contained in a phi0 curv hough space.
TrackingUtilities::Weight	operator() (const TrackingUtilities::CDCTrajectory2D *trajectory2D, const Phi0CurvBox *phi0CurvBox)
	Checks if the given trajectory is considered part of the hough box in phi0, curvature space.
TrackingUtilities::Weight	operator() (const TrackingUtilities::CDCFacet *facet, const Phi0CurvBox *phi0CurvBox)
	Checks if the given hit triplet is considered part of the hough box in phi0, curvature space.
bool	isPhi0CurvPointIn (const ROOT::Math::XYVector &phi0Vec, const double curvature, const Phi0CurvBox *phi0CurvBox)
	Predicate checking if the phi0 vector and curvature are contained in the given box.

Detailed Description

Predicate class to check for the containment of segments in a phi0 curv hough space part.

Unlike hits segments occupy only a point in the hough space.

Definition at line 28 of file SegmentInPhi0CurvBox.h.

Member Function Documentation

isPhi0CurvPointIn()

bool isPhi0CurvPointIn ( const ROOT::Math::XYVector & phi0Vec, const double curvature, const Phi0CurvBox * phi0CurvBox )

inlineprivate

Predicate checking if the phi0 vector and curvature are contained in the given box.

Definition at line 87 of file SegmentInPhi0CurvBox.h.

      {
 
        const ROOT::Math::XYVector& lowerPhi0Vec = phi0CurvBox->getLowerPhi0Vec();
        const ROOT::Math::XYVector& upperPhi0Vec = phi0CurvBox->getUpperPhi0Vec();
 
        // Allow containment to keep the reversal symmetry
        if (phi0CurvBox->isIn<1>(curvature)) {
          return VectorUtil::isBetween(phi0Vec, lowerPhi0Vec, upperPhi0Vec);
 
        } else if (phi0CurvBox->isIn<1>(-curvature)) {
          return VectorUtil::isBetween(-phi0Vec, lowerPhi0Vec, upperPhi0Vec);
 
        } else {
          return false;
 
        }
 
      }

operator()() [1/3]

TrackingUtilities::Weight operator() ( const TrackingUtilities::CDCFacet * facet, const Phi0CurvBox * phi0CurvBox )

inlineprivate

Checks if the given hit triplet is considered part of the hough box in phi0, curvature space.

Returns a finit weight if it is contained, NAN if not contained.

Definition at line 65 of file SegmentInPhi0CurvBox.h.

      {
        const ParameterLine2D& line = facet->getStartToEndLine();
 
        const ROOT::Math::XYVector pos2D = facet->getMiddleRecoPos2D();
        const ROOT::Math::XYVector phiVec = VectorUtil::unit(line.tangential());
 
        // Calculate the curvature and phi0 of the circle through the origin
        // that touches the position pos2D under the angle phiVec.
        double curvature = 2 * VectorUtil::Cross(pos2D, phiVec) / pos2D.Mag2();
        ROOT::Math::XYVector phi0Vec = VectorUtil::flippedOver(phiVec, pos2D);
 
        bool in = isPhi0CurvPointIn(phi0Vec, curvature, phi0CurvBox);
        return in ? 1 : NAN;
        // All facets are currently worth the same weight
        // Facets contain three hits but in general contribute only one
        // additional hit due to overlaps with neighboring facets.
        // Hence we stick with 1 as weight here.
      }

operator()() [2/3]

TrackingUtilities::Weight operator() ( const TrackingUtilities::CDCSegment2D * segment, const Phi0CurvBox * phi0CurvBox )

inlineprivate

Checks if the segment is contained in a phi0 curv hough space.

Returns the size of the segment if it is contained, returns NAN if it is not contained. The method simply assumes that the segment has been fitted with origin constraint and simply uses the values from the contained trajectory The calling code has to ensure this situation

Definition at line 37 of file SegmentInPhi0CurvBox.h.

      {
        if (segment->empty()) return NAN;
        const TrackingUtilities::CDCTrajectory2D& trajectory2D = segment->getTrajectory2D();
        TrackingUtilities::Weight weight = operator()(&trajectory2D, phi0CurvBox);
        return weight * segment->size();
      }

operator()() [3/3]

TrackingUtilities::Weight operator() ( const TrackingUtilities::CDCTrajectory2D * trajectory2D, const Phi0CurvBox * phi0CurvBox )

inlineprivate

Checks if the given trajectory is considered part of the hough box in phi0, curvature space.

Returns a finit weight if it is contained, NAN if not contained.

Definition at line 50 of file SegmentInPhi0CurvBox.h.

      {
        double curvature = trajectory2D->getCurvature();
        const ROOT::Math::XYVector& phi0Vec = trajectory2D->getStartUnitMom2D();
 
        bool in = isPhi0CurvPointIn(phi0Vec, curvature, phi0CurvBox);
        return in ? 1 : NAN;
      }

---

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

• tracking/trackFindingCDC/hough/algorithms/include/SegmentInPhi0CurvBox.h