Convenience class for the typical usage-case: A box divisioned hough tree with maximum and minimum values in both directions. More...

#include <SimpleBoxDivisionHoughTree.h>

Inheritance diagram for SimpleBoxDivisionHoughTree< AHitPtr, AInBoxAlgorithm, divisionX, divisionY >:

Public Types
using	BoxDivision = SectoredLinearDivision< HoughBox, divisions... >
	Type of the box division strategy.

using	HoughTree = WeightedFastHoughTree< AHitPtr, HoughBox, BoxDivision >
	Type of the fast hough tree structure.

using	Type = typename HoughBox::template Type< I >
	Type of the coordinate I.

using	HasType = typename HoughBox::template HasType< T >
	Predicate that the given type is indeed a coordinate of the hough space.

using	TypeIndex = typename HoughBox::template TypeIndex< T >
	Function to get the coordinate index from its type.

using	Node = typename HoughTree::Node
	Type of the nodes used in the tree for the search.

Public Member Functions
	SimpleBoxDivisionHoughTree (float maximumX, float maximumY, Width< 0 > overlapX=0, Width< 1 > overlapY=0)
	Constructor using the given maximal level.

void	initialize ()
	Initialize the tree with the given values.

std::vector< std::pair< HoughBox, std::vector< AHitPtr > > >	findSingleBest (const Weight &minWeight)
	Find only the leave with the highest weight = number of items.

float	getMaximumX () const
	Return the maximum value in x direction.

float	getMaximumY () const
	Return the maximum value in y direction.

Width< 0 >	getOverlapX () const
	Return the overlap in x direction.

Width< 1 >	getOverlapY () const
	Return the overlap in y direction.

size_t	getDivision (size_t i) const
	Getter the number of divisions at each level for coordinate index I.

void	constructArray (double lowerBound, double upperBound, Width< I > nBinOverlap=0, Width< I > nBinWidth=0)
	Construct the discrete value array at coordinate index I.

void	assignArray (Array< I > array, Width< I > overlap=0)
	Provide an externally constructed array by coordinate index.

std::enable_if_t< HasType< T >::value, void >	assignArray (Array< TypeIndex< T >::value > array, Width< TypeIndex< T >::value > overlap=0)
	Provide an externally constructed array by coordinate type.

void	seed (const AItemPtrs &items)
	Prepare the leave finding by filling the top node with given hits.

void	fell ()
	Terminates the processing by striping all hit information from the tree.

void	raze ()
	Release all memory that the tree aquired during the runs.

HoughTree *	getTree () const
	Getter for the tree used in the search in the hough plane.

int	getMaxLevel () const
	Getter for the currently set maximal level.

void	setMaxLevel (int maxLevel)
	Setter maximal level of the hough tree.

int	getSectorLevelSkip () const
	Getter for number of levels to skip in first level to form a finer sectored hough space.

void	setSectorLevelSkip (int sectorLevelSkip)
	Setter for number of levels to skip in first level to form a finer sectored hough space.

const Array< I > &	getArray () const
	Getter for the array of discrete value for coordinate I.

Private Types
using	Super = BoxDivisionHoughTree< AHitPtr, typename AInBoxAlgorithm::HoughBox, divisionX, divisionY >
	The Super class.

using	HoughBox = typename AInBoxAlgorithm::HoughBox
	The HoughBox we use.

template<size_t I>
using	Width = typename HoughBox::template Width< I >
	Type of the width in coordinate I.

using	Array = typename Type< I >::Array
	Type of the discrete value array to coordinate index I.

using	Arrays = TupleGenerateN< Array, sizeof...(divisions)>
	Tuple type of the discrete value arrays.

Private Member Functions
HoughBox	constructHoughPlaneImpl (const std::index_sequence< Is... > &is)
	Construct the box of the top node of the tree. Implementation unroling the indices.

HoughBox	constructHoughPlane ()
	Construct the box of the top node of the tree.

Private Attributes
float	m_maximumX = 0
	The maximum value in X direction.

float	m_maximumY = 0
	The maximum value in y direction.

Width< 0 >	m_overlapX = 0
	The overlap in X direction.

Width< 1 >	m_overlapY = 0
	The overlap in Y direction.

int	m_maxLevel
	Number of the maximum tree level.

int	m_sectorLevelSkip
	Number of levels to skip in first level to form a finer sectored hough space.

const std::array< size_t, sizeof ...(divisions)>	m_divisions
	Array of the number of divisions at each level.

HoughBox::Delta	m_overlaps
	An tuple of division overlaps in each coordinate.

Arrays	m_arrays
	A tuple of value arrays providing the memory for the discrete bin bounds.

std::unique_ptr< HoughTree >	m_houghTree
	Dynamic hough tree structure traversed in the leaf search.

Detailed Description

template<class AHitPtr, class AInBoxAlgorithm, size_t divisionX, size_t divisionY>
class Belle2::TrackFindingCDC::SimpleBoxDivisionHoughTree< AHitPtr, AInBoxAlgorithm, divisionX, divisionY >

Convenience class for the typical usage-case: A box divisioned hough tree with maximum and minimum values in both directions.

Definition at line 20 of file SimpleBoxDivisionHoughTree.h.

Member Typedef Documentation

◆ Array

using Array = typename Type<I>::Array

privateinherited

Type of the discrete value array to coordinate index I.

Definition at line 77 of file BoxDivisionHoughTree.h.

◆ Arrays

using Arrays = TupleGenerateN<Array, sizeof...(divisions)>

privateinherited

Tuple type of the discrete value arrays.

Definition at line 80 of file BoxDivisionHoughTree.h.

◆ BoxDivision

using BoxDivision = SectoredLinearDivision<HoughBox, divisions...>

inherited

Type of the box division strategy.

Definition at line 41 of file BoxDivisionHoughTree.h.

◆ HasType

using HasType = typename HoughBox::template HasType<T>

inherited

Predicate that the given type is indeed a coordinate of the hough space.

Definition at line 52 of file BoxDivisionHoughTree.h.

◆ HoughBox

using HoughBox = typename AInBoxAlgorithm::HoughBox

private

The HoughBox we use.

Definition at line 27 of file SimpleBoxDivisionHoughTree.h.

◆ HoughTree

using HoughTree = WeightedFastHoughTree<AHitPtr , HoughBox, BoxDivision>

inherited

Type of the fast hough tree structure.

Definition at line 44 of file BoxDivisionHoughTree.h.

◆ Node

using Node = typename HoughTree::Node

inherited

Type of the nodes used in the tree for the search.

Definition at line 63 of file BoxDivisionHoughTree.h.

◆ Super

using Super = BoxDivisionHoughTree<AHitPtr, typename AInBoxAlgorithm::HoughBox, divisionX, divisionY>

private

The Super class.

Definition at line 24 of file SimpleBoxDivisionHoughTree.h.

◆ Type

using Type = typename HoughBox::template Type<I>

inherited

Type of the coordinate I.

Definition at line 48 of file BoxDivisionHoughTree.h.

◆ TypeIndex

using TypeIndex = typename HoughBox::template TypeIndex<T>

inherited

Function to get the coordinate index from its type.

Definition at line 56 of file BoxDivisionHoughTree.h.

◆ Width

using Width = typename HoughBox::template Width<I>

private

Type of the width in coordinate I.

Definition at line 31 of file SimpleBoxDivisionHoughTree.h.

Constructor & Destructor Documentation

◆ SimpleBoxDivisionHoughTree()

SimpleBoxDivisionHoughTree	(	float	maximumX,
		float	maximumY,
		Width< 0 >	overlapX = `0`,
		Width< 1 >	overlapY = `0`
	)

inline

Constructor using the given maximal level.

Definition at line 35 of file SimpleBoxDivisionHoughTree.h.

        : Super(0)
        , m_maximumX(maximumX)
        , m_maximumY(maximumY)
        , m_overlapX(overlapX)
        , m_overlapY(overlapY)
      {
      }

Member Function Documentation

◆ assignArray() [1/2]

void assignArray	(	Array< I >	array,
		Width< I >	overlap = `0`
	)

inlineinherited

Provide an externally constructed array by coordinate index.

Definition at line 128 of file BoxDivisionHoughTree.h.

      {
        std::get<I>(m_arrays) = std::move(array);
        std::get<I>(m_overlaps) = overlap;
 
        // In case of a discrete axes, check whether the size of the array is sensible
        // such that the bin width at the highest granularity level is a whole number given the overlap.
        if (std::is_integral<Width<I> >::value) {
          const int division = getDivision(I);
          const int granularityLevel = m_maxLevel + m_sectorLevelSkip;
          const long int nBins = std::pow(division, granularityLevel);
          const long int nPositions = std::get<I>(m_arrays).size();
          const long int nWidthTimeNBins = nPositions - 1 + (nBins - 1) * overlap;
 
          B2ASSERT("Not enough positions in array to cover all bins.\n"
                   "Expected: positions = " << nBins - (nBins - 1) * overlap + 1 << " at least.\n"
                   "Actual:   positions = " << nPositions << " (overlap = " << overlap << ", bins = " << nBins << ")\n",
                   nWidthTimeNBins >= nBins);
 
          B2ASSERT("Number of positions in array introduces inhomogeneous width at the highest granularity level.\n"
                   "Expected: positions = 'width * bins - (bins - 1) * overlap + 1'\n"
                   "Actual:   positions = " << nPositions << " (overlap = " << overlap << ", bins = " << nBins << ")\n",
                   nWidthTimeNBins % nBins == 0);
        }
      }

◆ assignArray() [2/2]

std::enable_if_t< HasType< T >::value, void > assignArray	(	Array< TypeIndex< T >::value >	array,
		Width< TypeIndex< T >::value >	overlap = `0`
	)

inlineinherited

Provide an externally constructed array by coordinate type.

Definition at line 157 of file BoxDivisionHoughTree.h.

      {
        assignArray<TypeIndex<T>::value>(std::move(array), overlap);
      }

◆ constructArray()

void constructArray	(	double	lowerBound,
		double	upperBound,
		Width< I >	nBinOverlap = `0`,
		Width< I >	nBinWidth = `0`
	)

inlineinherited

Construct the discrete value array at coordinate index I.

This function is only applicable for discrete axes. For continuous axes assignArray should be call with an array containing only the lower and upper bound of the axes range and an optional overlap.

Parameters

lowerBound	Lower bound of the value range
upperBound	Upper bound of the value range
nBinOverlap	Overlap of neighboring bins. Default is no overlap. Usuallly this is counted in number of discrete values
nBinWidth	Width of the bins at lowest level. Default is width of 1. Usually this is counted in numbers of discrete values

Definition at line 104 of file BoxDivisionHoughTree.h.

      {
        static_assert(std::is_integral<Width<I> >::value, "Method only applicable to discrete axes");
        const size_t division = getDivision(I);
        const int granularityLevel = m_maxLevel + m_sectorLevelSkip;
        const size_t nBins = std::pow(division, granularityLevel);
 
        if (nBinWidth == 0) {
          nBinWidth = nBinOverlap + 1;
        }
 
        B2ASSERT("Width " << nBinWidth << "is not bigger than overlap " << nBinOverlap,
                 nBinOverlap < nBinWidth);
 
        const auto nPositions = (nBinWidth - nBinOverlap) * nBins + nBinOverlap + 1;
        std::get<I>(m_arrays) = linspace<float>(lowerBound, upperBound, nPositions);
        std::get<I>(m_overlaps) = nBinOverlap;
      }

◆ constructHoughPlane()

HoughBox constructHoughPlane ( )

inlineprivateinherited

Construct the box of the top node of the tree.

Definition at line 239 of file BoxDivisionHoughTree.h.

      {
        return constructHoughPlaneImpl(std::make_index_sequence<sizeof...(divisions)>());
      }

◆ constructHoughPlaneImpl()

HoughBox constructHoughPlaneImpl ( const std::index_sequence< Is... > & is )

inlineprivateinherited

Construct the box of the top node of the tree. Implementation unroling the indices.

Definition at line 233 of file BoxDivisionHoughTree.h.

      {
        return HoughBox(Type<Is>::getRange(std::get<Is>(m_arrays))...);
      }

◆ fell()

void fell ( )

inlineinherited

Terminates the processing by striping all hit information from the tree.

Definition at line 181 of file BoxDivisionHoughTree.h.

      {
        m_houghTree->fell();
      }

◆ findSingleBest()

std::vector< std::pair< HoughBox, std::vector< AHitPtr > > > findSingleBest ( const Weight & minWeight )

inline

Find only the leave with the highest weight = number of items.

Definition at line 58 of file SimpleBoxDivisionHoughTree.h.

      {
        AInBoxAlgorithm inBoxAlgorithm;
        auto skipLowWeightNode = [minWeight](const typename Super::Node * node) {
          return not(node->getWeight() >= minWeight);
        };
        auto found = this->getTree()->findHeaviestLeafSingle(inBoxAlgorithm, this->getMaxLevel(), skipLowWeightNode);
 
        std::vector<std::pair<HoughBox, std::vector<AHitPtr>>> result;
        if (found) {
          // Move the found content over. unique_ptr still destroys the left overs.
          result.push_back(std::move(*found));
        }
        return result;
      }

◆ getArray()

const Array< I > & getArray ( ) const

inlineinherited

Getter for the array of discrete value for coordinate I.

Definition at line 225 of file BoxDivisionHoughTree.h.

      {
        return std::get<I>(m_arrays);
      }

◆ getDivision()

size_t getDivision ( size_t i ) const

inlineinherited

Getter the number of divisions at each level for coordinate index I.

Definition at line 84 of file BoxDivisionHoughTree.h.

      {
        return m_divisions[i];
      }

◆ getMaximumX()

float getMaximumX ( ) const

inline

Return the maximum value in x direction.

Definition at line 75 of file SimpleBoxDivisionHoughTree.h.

      {
        return m_maximumX;
      }

◆ getMaximumY()

float getMaximumY ( ) const

inline

Return the maximum value in y direction.

Definition at line 81 of file SimpleBoxDivisionHoughTree.h.

      {
        return m_maximumY;
      }

◆ getMaxLevel()

int getMaxLevel ( ) const

inlineinherited

Getter for the currently set maximal level.

Definition at line 200 of file BoxDivisionHoughTree.h.

      {
        return m_maxLevel;
      }

◆ getOverlapX()

Width< 0 > getOverlapX ( ) const

inline

Return the overlap in x direction.

Definition at line 87 of file SimpleBoxDivisionHoughTree.h.

      {
        return m_overlapX;
      }

◆ getOverlapY()

Width< 1 > getOverlapY ( ) const

inline

Return the overlap in y direction.

Definition at line 93 of file SimpleBoxDivisionHoughTree.h.

      {
        return m_overlapY;
      }

◆ getSectorLevelSkip()

int getSectorLevelSkip ( ) const

inlineinherited

Getter for number of levels to skip in first level to form a finer sectored hough space.

Definition at line 212 of file BoxDivisionHoughTree.h.

      {
        return m_sectorLevelSkip;
      }

◆ getTree()

HoughTree * getTree ( ) const

inlineinherited

Getter for the tree used in the search in the hough plane.

Definition at line 194 of file BoxDivisionHoughTree.h.

      {
        return m_houghTree.get();
      }

◆ initialize()

void initialize ( )

inline

Initialize the tree with the given values.

Definition at line 48 of file SimpleBoxDivisionHoughTree.h.

      {
        Super::template constructArray<0>(-getMaximumX(), getMaximumX(), getOverlapX());
        Super::template constructArray<1>(-getMaximumY(), getMaximumY(), getOverlapY());
 
        Super::initialize();
      }

◆ raze()

void raze ( )

inlineinherited

Release all memory that the tree aquired during the runs.

Definition at line 187 of file BoxDivisionHoughTree.h.

      {
        m_houghTree->raze();
      }

◆ seed()

void seed ( const AItemPtrs & items )

inlineinherited

Prepare the leave finding by filling the top node with given hits.

Definition at line 174 of file BoxDivisionHoughTree.h.

      {
        if (not m_houghTree) initialize();
        m_houghTree->seed(items);
      }

◆ setMaxLevel()

void setMaxLevel ( int maxLevel )

inlineinherited

Setter maximal level of the hough tree.

Definition at line 206 of file BoxDivisionHoughTree.h.

      {
        m_maxLevel = maxLevel;
      }

◆ setSectorLevelSkip()

void setSectorLevelSkip ( int sectorLevelSkip )

inlineinherited

Setter for number of levels to skip in first level to form a finer sectored hough space.

Definition at line 218 of file BoxDivisionHoughTree.h.

      {
        m_sectorLevelSkip = sectorLevelSkip;
      }

Member Data Documentation

◆ m_arrays

Arrays m_arrays

privateinherited

A tuple of value arrays providing the memory for the discrete bin bounds.

Definition at line 258 of file BoxDivisionHoughTree.h.

◆ m_divisions

const std::array<size_t, sizeof ...(divisions)> m_divisions

privateinherited

Array of the number of divisions at each level.

Definition at line 252 of file BoxDivisionHoughTree.h.

◆ m_houghTree

std::unique_ptr<HoughTree> m_houghTree

privateinherited

Dynamic hough tree structure traversed in the leaf search.

Definition at line 261 of file BoxDivisionHoughTree.h.

◆ m_maximumX

float m_maximumX = 0

private

The maximum value in X direction.

Definition at line 100 of file SimpleBoxDivisionHoughTree.h.

◆ m_maximumY

float m_maximumY = 0

private

The maximum value in y direction.

Definition at line 103 of file SimpleBoxDivisionHoughTree.h.

◆ m_maxLevel

int m_maxLevel

privateinherited

Number of the maximum tree level.

Definition at line 246 of file BoxDivisionHoughTree.h.

◆ m_overlaps

HoughBox::Delta m_overlaps

privateinherited

An tuple of division overlaps in each coordinate.

Definition at line 255 of file BoxDivisionHoughTree.h.

◆ m_overlapX

Width<0> m_overlapX = 0

private

The overlap in X direction.

Definition at line 106 of file SimpleBoxDivisionHoughTree.h.

◆ m_overlapY

Width<1> m_overlapY = 0

private

The overlap in Y direction.

Definition at line 109 of file SimpleBoxDivisionHoughTree.h.

◆ m_sectorLevelSkip

int m_sectorLevelSkip

privateinherited

Number of levels to skip in first level to form a finer sectored hough space.

Definition at line 249 of file BoxDivisionHoughTree.h.

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

tracking/trackFindingCDC/hough/trees/include/SimpleBoxDivisionHoughTree.h

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ Array

◆ Arrays

◆ BoxDivision

◆ HasType

◆ HoughBox

◆ HoughTree

◆ Node

◆ Super

◆ Type

◆ TypeIndex

◆ Width

Constructor & Destructor Documentation

◆ SimpleBoxDivisionHoughTree()

Member Function Documentation

◆ assignArray() [1/2]

◆ assignArray() [2/2]

◆ constructArray()

◆ constructHoughPlane()

◆ constructHoughPlaneImpl()

◆ fell()

◆ findSingleBest()

◆ getArray()

◆ getDivision()

◆ getMaximumX()

◆ getMaximumY()

◆ getMaxLevel()

◆ getOverlapX()

◆ getOverlapY()

◆ getSectorLevelSkip()

◆ getTree()

◆ initialize()

◆ raze()

◆ seed()

◆ setMaxLevel()

◆ setSectorLevelSkip()

Member Data Documentation

◆ m_arrays

◆ m_divisions

◆ m_houghTree

◆ m_maximumX

◆ m_maximumY

◆ m_maxLevel

◆ m_overlaps

◆ m_overlapX

◆ m_overlapY

◆ m_sectorLevelSkip