InPhi0CurvTanLBox Class Reference

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

#include <InPhi0CurvTanLBox.h>

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

Public Member Functions
	InPhi0CurvTanLBox (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 iCLayer=-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 32 of file InPhi0CurvTanLBox.h.

Member Function Documentation

getDistanceSign()

ESign getDistanceSign ( const HoughBox &  houghBox, float  x, float  y, float  l, float  dxdz, float  dydz, ILayer  iCLayer = -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 46 of file InPhi0CurvTanLBox.h.

                             {
          // 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;
        }
 
        std::array<float, 2> s;
        s[0] = curv[0]->getArcLength2D(iCLayer, (xRot[0] < 0) and (xRot[1] < 0));
        s[1] = curv[1]->getArcLength2D(iCLayer, (xRot[0] < 0) and (xRot[1] < 0));
 
        std::array<float, 2> r2TimesHalfCurv;
        r2 -= l * l;
        r2TimesHalfCurv[0] = r2 * (static_cast<float>(*(curv[0])) / 2.0);
        r2TimesHalfCurv[1] = r2 * (static_cast<float>(*(curv[1])) / 2.0);
 
        // 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> z;
        z[0b00] = s[0] * static_cast<float>(tanL[0]);
        z[0b01] = s[1] * static_cast<float>(tanL[0]);
        z[0b10] = s[0] * static_cast<float>(tanL[1]);
        z[0b11] = s[1] * static_cast<float>(tanL[1]);
 
        // Using binary notation encoding lower and upper box bounds to fill the flat array.
        std::array<float, 8> dist;
        dist[0b000] = r2TimesHalfCurv[0] - yRotPlusL[0] - dydzRot[0] * z[0b00];
        dist[0b001] = r2TimesHalfCurv[0] - yRotPlusL[1] - dydzRot[1] * z[0b00];
        dist[0b010] = r2TimesHalfCurv[1] - yRotPlusL[0] - dydzRot[0] * z[0b01];
        dist[0b011] = r2TimesHalfCurv[1] - yRotPlusL[1] - dydzRot[1] * z[0b01];
 
        dist[0b100] = r2TimesHalfCurv[0] - yRotPlusL[0] - dydzRot[0] * z[0b10];
        dist[0b101] = r2TimesHalfCurv[0] - yRotPlusL[1] - dydzRot[1] * z[0b10];
        dist[0b110] = r2TimesHalfCurv[1] - yRotPlusL[0] - dydzRot[0] * z[0b11];
        dist[0b111] = r2TimesHalfCurv[1] - yRotPlusL[1] - dydzRot[1] * z[0b11];
 
        return ESignUtil::common(dist);
      }
 
    private:
      float m_curlCurv;
    };
  }
}

---

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

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