InPhi0TanLBox Class Reference

Checker if a position is contained in a family of curves over phi0 and tan lambda. More...

#include <InPhi0TanLBox.h>

Public Types
using	HoughBox = Box< DiscretePhi0, ContinuousTanL >
	The box to which this object correspondes.

Public Member Functions
	InPhi0TanLBox (float curlCurv)
	Create a new box with the given curler curvature.
ESign	getDistanceSign (const HoughBox &houghBox, float x, float y, float l, float dxdz, float dydz, ILayer=-1) const
	Function that gives the sign of the distance from an observed drift circle to the sweeped object. More...

Private Attributes
float	m_curlCurv
	Curler curvature - set to value greater zero to activate one arm exclusive finding.

Detailed Description

Checker if a position is contained in a family of curves over phi0 and tan lambda.

Definition at line 22 of file InPhi0TanLBox.h.

Member Function Documentation

getDistanceSign()

ESign getDistanceSign ( const HoughBox &  houghBox, float  x, float  y, float  l, float  dxdz, float  dydz, ILayer  = -1  ) const

inline

Function that gives the sign of the distance from an observed drift circle to the sweeped object.

Two dimensional arc length

Definition at line 36 of file InPhi0TanLBox.h.

      {
        const std::array<DiscretePhi0, 2>& phi0Vec = houghBox.getBounds<DiscretePhi0>();
        const std::array<ContinuousTanL, 2>& tanL = houghBox.getBounds<ContinuousTanL>();
 
        std::array<float, 2> xRot;
        xRot[0] = x * phi0Vec[0]->x() + y * phi0Vec[0]->y();
        xRot[1] = x * phi0Vec[1]->x() + y * phi0Vec[1]->y();
 
        std::array<float, 2> yRot;
        yRot[0] = -x * phi0Vec[0]->y() + y * phi0Vec[0]->x();
        yRot[1] = -x * phi0Vec[1]->y() + y * phi0Vec[1]->x();
 
        std::array<float, 2> dxdzRot;
        dxdzRot[0] = dxdz * phi0Vec[0]->x() + dydz * phi0Vec[0]->y();
        dxdzRot[1] = dxdz * phi0Vec[1]->x() + dydz * phi0Vec[1]->y();
 
        std::array<float, 2> dydzRot;
        dydzRot[0] = -dxdz * phi0Vec[0]->y() + dydz * phi0Vec[0]->x();
        dydzRot[1] = -dxdz * phi0Vec[1]->y() + dydz * phi0Vec[1]->x();
 
        const std::array<float, 2>& s = xRot;
 
        const bool onlyPositiveArm = 0 < m_curlCurv;
        if (onlyPositiveArm) {
          // Reject hit if it is on the inward going branch but the curvature suggest it is no curler
          if (xRot[0] < 0 and xRot[1] < 0) return ESign::c_Invalid;
        }
 
        // cppcheck-suppress unreadVariable
        std::array<float, 4> xRotCor;
        xRotCor[0b00] = xRot[0] + dxdzRot[0] * s[0] * static_cast<float>(tanL[0]);
        xRotCor[0b01] = xRot[0] + dxdzRot[0] * s[0] * static_cast<float>(tanL[1]);
        xRotCor[0b10] = xRot[1] + dxdzRot[1] * s[1] * static_cast<float>(tanL[0]);
        xRotCor[0b11] = xRot[1] + dxdzRot[1] * s[1] * static_cast<float>(tanL[1]);
 
        const std::array<float, 4>& sCor = xRotCor;
 
        std::array<float, 4> dist;
        dist[0b00] = -(yRot[0] + dydzRot[0] * sCor[0b00] * static_cast<float>(tanL[0]) + l);
        dist[0b01] = -(yRot[0] + dydzRot[0] * sCor[0b01] * static_cast<float>(tanL[1]) + l);
        dist[0b10] = -(yRot[1] + dydzRot[1] * sCor[0b10] * static_cast<float>(tanL[0]) + l);
        dist[0b11] = -(yRot[1] + dydzRot[1] * sCor[0b11] * static_cast<float>(tanL[1]) + l);
 
        return ESignUtil::common(dist);
      }

---

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

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