Generates axial tracks from hit using special leaf postprocessing. More...

#include <AxialTrackCreatorHitLegendre.h>

Inheritance diagram for AxialTrackCreatorHitLegendre:

Public Types
enum class	EPass { NonCurlers , NonCurlersWithIncreasingThreshold , FullRange , Straight }
	Pass keys for the different sets of predefined parameters for a pass if legendre search Note: Naming copied from original location. More...

using	IOTypes = std::tuple< AIOTypes... >
	Types that should be served to apply on invokation.

using	IOVectors = std::tuple< std::vector< AIOTypes >... >
	Vector types that should be served to apply on invokation.

Public Member Functions
	AxialTrackCreatorHitLegendre ()
	Constructor.

	AxialTrackCreatorHitLegendre (EPass pass)
	Constructor from a pass key.

std::string	getDescription () final
	Short description of the findlet.

void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) final
	Expose the parameters to a module.

void	initialize () final
	Initialisation before the event processing starts.

void	apply (const std::vector< const CDCWireHit * > &axialWireHits, std::vector< CDCTrack > &tracks) final
	Execute one pass over a quad tree.

virtual void	apply (ToVector< AIOTypes > &... ioVectors)=0
	Main function executing the algorithm.

void	beginRun () override
	Receive and dispatch signal for the beginning of a new run.

void	beginEvent () override
	Receive and dispatch signal for the start of a new event.

void	endRun () override
	Receive and dispatch signal for the end of the run.

void	terminate () override
	Receive and dispatch Signal for termination of the event processing.

Protected Types
using	ToVector = typename ToVectorImpl< T >::Type
	Short hand for ToRangeImpl.

Protected Member Functions
void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.

int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Private Types
using	Super = Findlet< const CDCWireHit *const, CDCTrack >
	Type of the base class.

using	CandidateReceiver = AxialHitQuadTreeProcessor::CandidateReceiver
	lambda function used for postprocessing

Private Member Functions
std::unique_ptr< AxialHitQuadTreeProcessor >	constructQTProcessor (EPass pass)
	Method to create QTProcessor that performs the search.

void	executeRelaxation (const CandidateReceiver &candidateReceiver, AxialHitQuadTreeProcessor &qtProcessor)
	Performs quadtree search.

Private Attributes
EPass	m_pass = EPass::NonCurlers
	The pass key for lookup of the parameters for this pass.

const double	m_param_stepScale = 0.75
	Parameter to define multiplier for hits threshold for the next quadtree iteration.

int	m_param_minNHits = 10
	Parameter to define minimal threshold of hit.

double	m_param_precision = 0.00000001
	Parameter to define precision of quadtree search in case of straight pass.

std::vector< ProcessingSignalListener * >	m_subordinaryProcessingSignalListeners
	References to subordinary signal processing listener contained in this findlet.

bool	m_initialized = false
	Flag to keep track whether initialization happend before.

bool	m_terminated = false
	Flag to keep track whether termination happend before.

std::string	m_initializedAs
	Name of the type during initialisation.

Detailed Description

Generates axial tracks from hit using special leaf postprocessing.

Implements Legendre transformation of the drift time circles. This is a module, performing tracking in the CDC. It is based on the paper "Implementation of the Legendre Transform for track segment reconstruction in drift tube chambers" by T. Alexopoulus, et al. NIM A592 456-462 (2008).

Definition at line 32 of file AxialTrackCreatorHitLegendre.h.

Member Typedef Documentation

◆ CandidateReceiver

using CandidateReceiver = AxialHitQuadTreeProcessor::CandidateReceiver

private

lambda function used for postprocessing

Definition at line 82 of file AxialTrackCreatorHitLegendre.h.

◆ IOTypes

using IOTypes = std::tuple<AIOTypes...>

inherited

Types that should be served to apply on invokation.

Definition at line 30 of file Findlet.h.

◆ IOVectors

using IOVectors = std::tuple< std::vector<AIOTypes>... >

inherited

Vector types that should be served to apply on invokation.

Definition at line 53 of file Findlet.h.

◆ Super

using Super = Findlet<const CDCWireHit* const, CDCTrack>

private

Type of the base class.

Definition at line 36 of file AxialTrackCreatorHitLegendre.h.

◆ ToVector

using ToVector = typename ToVectorImpl<T>::Type

protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Member Enumeration Documentation

◆ EPass

enum class EPass

strong

Pass keys for the different sets of predefined parameters for a pass if legendre search Note: Naming copied from original location.

Does not actually match with the associated parameters.

Enumerator
NonCurlers	Pass corresponds to High-pt track finding and more deeper quadtree.
NonCurlersWithIncreasingThreshold	Pass corresponds to High-pt track finding and more rough quadtree.
FullRange	Pass corresponds to full pt range and even more rough quadtree (non-ip tracks, tracks with energy losses etc)
Straight	Pass corresponds to very High-pt track and very rough quadtree (monopole tracks, including those with low hit count)

Definition at line 44 of file AxialTrackCreatorHitLegendre.h.

                       {
        NonCurlers,
        NonCurlersWithIncreasingThreshold,
        FullRange,
        Straight
      };

Constructor & Destructor Documentation

◆ AxialTrackCreatorHitLegendre()

AxialTrackCreatorHitLegendre ( EPass pass )

explicit

Constructor from a pass key.

Definition at line 69 of file AxialTrackCreatorHitLegendre.cc.

  : m_pass(pass)
{
}

Member Function Documentation

◆ addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener * psl )

protectedinherited

Register a processing signal listener to be notified.

Definition at line 55 of file CompositeProcessingSignalListener.cc.

{
  m_subordinaryProcessingSignalListeners.push_back(psl);
}

◆ apply()

void apply	(	const std::vector< const CDCWireHit * > &	axialWireHits,
		std::vector< CDCTrack > &	tracks
	)

final

Execute one pass over a quad tree.

Definition at line 99 of file AxialTrackCreatorHitLegendre.cc.

{
  // Prepare vector of unused hits to provide to the qt processor
  // Also reset the mask flag and select only the untaken hits
  std::vector<const CDCWireHit*> unusedAxialWireHits;
  for (const CDCWireHit* wireHit : axialWireHits) {
    (*wireHit)->setMaskedFlag(false);
    if ((*wireHit)->hasTakenFlag()) continue;
    unusedAxialWireHits.push_back(wireHit);
  }
 
  // Create quadtree processor
  std::unique_ptr<AxialHitQuadTreeProcessor> qtProcessor = constructQTProcessor(m_pass);
  qtProcessor->seed(unusedAxialWireHits);
//   qtProcessor->drawHits(unusedAxialWireHits, 9);
 
  // Create object which contains interface between quadtree processor and track processor (module)
  std::unique_ptr<BaseCandidateReceiver> receiver;
  if (not(m_pass == EPass::Straight)) {
    receiver = std::make_unique<OffOriginExtension>(unusedAxialWireHits);
  } else {
    receiver = std::make_unique<BaseCandidateReceiver>(unusedAxialWireHits);
  }
 
  // Start candidate finding
  this->executeRelaxation(std::ref(*receiver), *qtProcessor);
 
  const std::vector<CDCTrack>& newTracks = receiver->getTracks();
  tracks.insert(tracks.end(), newTracks.begin(), newTracks.end());
}

◆ beginEvent()

void beginEvent ( )

overridevirtualinherited

Receive and dispatch signal for the start of a new event.

Reimplemented from ProcessingSignalListener.

Definition at line 31 of file CompositeProcessingSignalListener.cc.

{
  Super::beginEvent();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginEvent();
  }
}

◆ beginRun()

void beginRun ( )

overridevirtualinherited

Receive and dispatch signal for the beginning of a new run.

Reimplemented from ProcessingSignalListener.

Reimplemented in LayerRelationFilter< AFilter >, FourHitFilter, QualityIndicatorFilter, ThreeHitFilter, TwoHitVirtualIPFilter, TwoHitVirtualIPQIFilter, RecoTrackStorer, ROIFinder, SpacePointLoaderAndPreparer, and TrackCandidateResultRefiner.

Definition at line 23 of file CompositeProcessingSignalListener.cc.

{
  Super::beginRun();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->beginRun();
  }
}

◆ constructQTProcessor()

std::unique_ptr< AxialHitQuadTreeProcessor > constructQTProcessor ( EPass pass )

private

Method to create QTProcessor that performs the search.

Definition at line 23 of file AxialTrackCreatorHitLegendre.cc.

{
  using XYSpans = AxialHitQuadTreeProcessor::XYSpans;
  using PrecisionFunction = PrecisionUtil::PrecisionFunction;
  const int maxTheta = std::pow(2, PrecisionUtil::getLookupGridLevel());
 
  if (pass == EPass::NonCurlers) {
    int maxLevel = 12;
    int seedLevel = 4;
    XYSpans xySpans({{0, maxTheta}, { -0.02, 0.14}});
    PrecisionFunction precisionFunction = &PrecisionUtil::getOriginCurvPrecision;
 
    return std::make_unique<AxialHitQuadTreeProcessor>(maxLevel, seedLevel, xySpans, precisionFunction);
 
  } else if (pass == EPass::NonCurlersWithIncreasingThreshold) {
    int maxLevel = 10;
    int seedLevel = 4;
    XYSpans xySpans({{0, maxTheta}, { -0.02, 0.14}});
    PrecisionFunction precisionFunction = &PrecisionUtil::getNonOriginCurvPrecision;
 
    return std::make_unique<AxialHitQuadTreeProcessor>(maxLevel, seedLevel, xySpans, precisionFunction);
 
  } else if (pass == EPass::FullRange) {
    int maxLevel = 10;
    int seedLevel = 1;
    XYSpans xySpans({{0, maxTheta}, {0.00, 0.30}});
    PrecisionFunction precisionFunction = &PrecisionUtil::getNonOriginCurvPrecision;
 
    return std::make_unique<AxialHitQuadTreeProcessor>(maxLevel, seedLevel, xySpans, precisionFunction);
 
  } else if (pass == EPass::Straight) {
    int maxLevel = 10;
    int seedLevel = 4;
    XYSpans xySpans({{0, maxTheta}, { -0.02, 0.02}});
//       PrecisionFunction precisionFunction = &PrecisionUtil::getBasicCurvPrecision; //That is 0.3 / pow(2, 16)
    PrecisionFunction precisionFunction = [this](double curv __attribute__((unused))) {return m_param_precision ;};
 
    return std::make_unique<AxialHitQuadTreeProcessor>(maxLevel, seedLevel, xySpans, precisionFunction);
 
  }
  B2FATAL("Invalid pass");
}

◆ endRun()

void endRun ( )

overridevirtualinherited

Receive and dispatch signal for the end of the run.

Reimplemented from ProcessingSignalListener.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->endRun();
  }
  Super::endRun();
}

◆ executeRelaxation()

void executeRelaxation	(	const CandidateReceiver &	candidateReceiver,
		AxialHitQuadTreeProcessor &	qtProcessor
	)

private

Performs quadtree search.

Parameters

candidateReceiver	lambda interface process found leaves in the legendre search
qtProcessor	reference to the AxialHitQuadTreeProcessor instance

Definition at line 131 of file AxialTrackCreatorHitLegendre.cc.

{
  // radius of the CDC
  const double rCDC = 113.;
 
  // Curvature for high pt particles that leave the CDC
  const double curlCurv = 2. / rCDC;
 
  // find leavers
  qtProcessor.fill(candidateReceiver, 50, curlCurv);
 
  // find curlers with diameter higher than half of radius of CDC
  qtProcessor.fill(candidateReceiver, 70, 2 * curlCurv);
 
  // Start relaxation loop
  int minNHits = m_pass == EPass::FullRange ? 30 : 50;
  double maxCurv = m_pass == EPass::FullRange ? 0.15 : 0.07;
  do {
    qtProcessor.fill(candidateReceiver, minNHits, maxCurv);
 
    minNHits = minNHits * m_param_stepScale;
 
    if (m_pass != EPass::NonCurlers) {
      maxCurv *= 2.;
      if (maxCurv > 0.15) maxCurv = 0.15;
    }
 
  } while (minNHits >= m_param_minNHits);
}

◆ exposeParameters()

void exposeParameters	(	ModuleParamList *	moduleParamList,
		const std::string &	prefix
	)

finalvirtual

Expose the parameters to a module.

Reimplemented from Findlet< const CDCWireHit *const, CDCTrack >.

Definition at line 79 of file AxialTrackCreatorHitLegendre.cc.

{
  moduleParamList->addParameter(prefixed(prefix, "minNHits"),
                                m_param_minNHits,
                                "Parameter to define minimal threshold of number of hits.",
                                m_param_minNHits);
  if (m_pass == EPass::Straight) {
    moduleParamList->addParameter(prefixed(prefix, "precision"),
                                  m_param_precision,
                                  "Parameter to define precision of quadtree search.",
                                  m_param_precision);
  }
}

◆ getDescription()

std::string getDescription ( )

finalvirtual

Short description of the findlet.

Reimplemented from Findlet< const CDCWireHit *const, CDCTrack >.

Definition at line 74 of file AxialTrackCreatorHitLegendre.cc.

{
  return "Generates axial tracks from hits using several increasingly relaxed legendre space search over phi0 and curvature.";
}

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )

protectedinherited

Get the number of currently registered listeners.

Definition at line 60 of file CompositeProcessingSignalListener.cc.

{
  return m_subordinaryProcessingSignalListeners.size();
}

◆ initialize()

void initialize ( )

finalvirtual

Initialisation before the event processing starts.

Reimplemented from ProcessingSignalListener.

Definition at line 94 of file AxialTrackCreatorHitLegendre.cc.

{
  Super::initialize();
}

◆ terminate()

void terminate ( )

overridevirtualinherited

Receive and dispatch Signal for termination of the event processing.

Reimplemented from ProcessingSignalListener.

Reimplemented in StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::HyperHough >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::QuadraticLegendre >, and StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::Z0TanLambdaLegendre >.

Definition at line 47 of file CompositeProcessingSignalListener.cc.

{
  for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
    psl->terminate();
  }
  Super::terminate();
}

Member Data Documentation

◆ m_initialized

bool m_initialized = false

privateinherited

Flag to keep track whether initialization happend before.

Definition at line 52 of file ProcessingSignalListener.h.

◆ m_initializedAs

std::string m_initializedAs

privateinherited

Name of the type during initialisation.

Definition at line 58 of file ProcessingSignalListener.h.

◆ m_param_minNHits

int m_param_minNHits = 10

private

Parameter to define minimal threshold of hit.

Definition at line 99 of file AxialTrackCreatorHitLegendre.h.

◆ m_param_precision

double m_param_precision = 0.00000001

private

Parameter to define precision of quadtree search in case of straight pass.

Definition at line 102 of file AxialTrackCreatorHitLegendre.h.

◆ m_param_stepScale

const double m_param_stepScale = 0.75

private

Parameter to define multiplier for hits threshold for the next quadtree iteration.

Definition at line 96 of file AxialTrackCreatorHitLegendre.h.

◆ m_pass

EPass m_pass = EPass::NonCurlers

private

The pass key for lookup of the parameters for this pass.

Definition at line 93 of file AxialTrackCreatorHitLegendre.h.

◆ m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners

privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 52 of file CompositeProcessingSignalListener.h.

◆ m_terminated

bool m_terminated = false

privateinherited

Flag to keep track whether termination happend before.

Definition at line 55 of file ProcessingSignalListener.h.

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

tracking/trackFindingCDC/findlets/minimal/include/AxialTrackCreatorHitLegendre.h
tracking/trackFindingCDC/findlets/minimal/src/AxialTrackCreatorHitLegendre.cc

Public Types

Public Member Functions

Protected Types

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ CandidateReceiver

◆ IOTypes

◆ IOVectors

◆ Super

◆ ToVector

Member Enumeration Documentation

◆ EPass

Constructor & Destructor Documentation

◆ AxialTrackCreatorHitLegendre()

Member Function Documentation

◆ addProcessingSignalListener()

◆ apply()

◆ beginEvent()

◆ beginRun()

◆ constructQTProcessor()

◆ endRun()

◆ executeRelaxation()

◆ exposeParameters()

◆ getDescription()

◆ getNProcessingSignalListener()

◆ initialize()

◆ terminate()

Member Data Documentation

◆ m_initialized

◆ m_initializedAs

◆ m_param_minNHits

◆ m_param_precision

◆ m_param_stepScale

◆ m_pass

◆ m_subordinaryProcessingSignalListeners

◆ m_terminated