CKFToSVDFindlet Class Reference

Combinatorial Kalman Filter to extrapolate CDC Reco Tracks into the VXD (SVD) and collect space points. More...

#include <CKFToSVDFindlet.h>

Inheritance diagram for CKFToSVDFindlet:

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

Public Member Functions
	CKFToSVDFindlet ()
	Constructor, for setting module description and parameters.
	~CKFToSVDFindlet () override
	Default destructor.
void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) override
	Expose the parameters of the sub findlets.
void	apply () override
	Do the track/hit finding/merging.
void	beginEvent () override
	Clear the object pools.
virtual std::string	getDescription ()
	Brief description of the purpose of the concrete findlet.
virtual void	apply (ToVector< AIOTypes > &... ioVectors)=0
	Main function executing the algorithm.
void	initialize () override
	Receive and dispatch signal before the start of the event processing.
void	beginRun () override
	Receive and dispatch signal for the beginning of a new run.
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 = TrackFindingCDC::Findlet<>
	Parent class.

Private Attributes
unsigned int	m_param_minimalHitRequirement = 1
	Minimal hit requirement for the results (counted in number of space points)
TrackLoader	m_dataHandler
	Findlet for retrieving the cdc tracks and writing the result out.
SpacePointLoader	m_hitsLoader
	Findlet for loading the space points.
StateCreatorWithReversal< CKFToSVDState >	m_stateCreatorFromTracks
	Findlet for creating states out of tracks.
StateCreator< const SpacePoint, CKFToSVDState >	m_stateCreatorFromHits
	Findlet for creating states out of hits.
CKFRelationCreator< CKFToSVDState, ChooseableSVDRelationFilter >	m_relationCreator
	Findlet for creating relations between states.
TreeSearcher< CKFToSVDState, SVDStateRejecter, CKFToSVDResult >	m_treeSearchFindlet
	Findlet doing the main work: the tree finding.
OverlapResolver< ChooseableSVDResultFilter >	m_overlapResolver
	Findlet for resolving overlaps.
SpacePointTagger< CKFToSVDResult, SVDCluster >	m_spacePointTagger
	Findlet for tagging the used space points.
ResultStorer< CKFToSVDResult >	m_resultStorer
	Findlet for storing the results.
std::vector< RecoTrack * >	m_cdcRecoTrackVector
	Pointers to the CDC Reco tracks as a vector.
std::vector< const SpacePoint * >	m_spacePointVector
	Pointers to the (const) SpacePoints as a vector.
std::vector< CKFToSVDState >	m_seedStates
	States for the tracks.
std::vector< CKFToSVDState >	m_states
	States for the hits.
std::vector< TrackFindingCDC::WeightedRelation< CKFToSVDState > >	m_relations
	Relations between states.
std::vector< CKFToSVDResult >	m_results
	Vector for storing the results.
std::vector< CKFToSVDResult >	m_filteredResults
	Vector for storing the filtered results.
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 happened before.
bool	m_terminated = false
	Flag to keep track whether termination happened before.
std::string	m_initializedAs
	Name of the type during initialisation.

Detailed Description

Combinatorial Kalman Filter to extrapolate CDC Reco Tracks into the VXD (SVD) and collect space points.

The implementation is split up in four parts and factored out into three sub findlets.

• Fetch the SpacePoints and the reco tracks from the data store (CDCTrackSpacePointStoreArrayHandler)
• Construct all possible candidates starting from a RecoTrack and going through the layers of the VXD collecting space points (this is handles by the TreeSearchFindlet, which works on States. The selection of space points is handled by the CDCToSpacePointHitSelector)
• Find a non-overlapping set of results (only one candidate per space point and seed) (OverlapResolverFindlet, quality is determined by a filter)
• Write the results out to the data store (again the CDCTrackSpacePointStoreArrayHandler)

If you want to reimplement this algorithm for a different purpose, you probably only have to implement a new state class and a new hit selector (and maybe the store array handling). The rest should be taken care by the framework.

Definition at line 56 of file CKFToSVDFindlet.h.

Member Typedef Documentation

IOTypes

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

inherited

Types that should be served to apply on invocation.

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 invocation.

Definition at line 53 of file Findlet.h.

Super

using Super = TrackFindingCDC::Findlet<>

private

Parent class.

Definition at line 58 of file CKFToSVDFindlet.h.

ToVector

using ToVector = typename ToVectorImpl<T>::Type

protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Constructor & Destructor Documentation

CKFToSVDFindlet()

CKFToSVDFindlet

(

)

Constructor, for setting module description and parameters.

Definition at line 31 of file CKFToSVDFindlet.cc.

{
  addProcessingSignalListener(&m_dataHandler);
  addProcessingSignalListener(&m_hitsLoader);
  addProcessingSignalListener(&m_stateCreatorFromTracks);
  addProcessingSignalListener(&m_stateCreatorFromHits);
  addProcessingSignalListener(&m_relationCreator);
  addProcessingSignalListener(&m_treeSearchFindlet);
  addProcessingSignalListener(&m_overlapResolver);
  addProcessingSignalListener(&m_spacePointTagger);
  addProcessingSignalListener(&m_resultStorer);
}

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 ( )

override

Do the track/hit finding/merging.

Definition at line 100 of file CKFToSVDFindlet.cc.

{
  m_dataHandler.apply(m_cdcRecoTrackVector);
  m_hitsLoader.apply(m_spacePointVector);
 
  B2DEBUG(29, "Now have " << m_spacePointVector.size() << " hits.");
 
  if (m_spacePointVector.empty() or m_cdcRecoTrackVector.empty()) {
    return;
  }
 
  m_stateCreatorFromTracks.apply(m_cdcRecoTrackVector, m_seedStates);
  m_stateCreatorFromHits.apply(m_spacePointVector, m_states);
  m_relationCreator.apply(m_seedStates, m_states, m_relations);
 
  B2DEBUG(29, "Created " << m_relations.size() << " relations.");
  m_treeSearchFindlet.apply(m_seedStates, m_states, m_relations, m_results);
 
  B2DEBUG(29, "Having found " << m_results.size() << " results before overlap check");
 
  const auto hasLowHitNumber = [this](const CKFResult<RecoTrack, SpacePoint>& result) {
    return result.getHits().size() < m_param_minimalHitRequirement;
  };
  TrackFindingCDC::erase_remove_if(m_results, hasLowHitNumber);
 
  m_overlapResolver.apply(m_results, m_filteredResults);
 
  B2DEBUG(29, "Having found " << m_filteredResults.size() << " results");
 
  m_resultStorer.apply(m_filteredResults);
  m_spacePointTagger.apply(m_filteredResults, m_spacePointVector);
}

beginEvent()

void beginEvent ( )

overridevirtual

Clear the object pools.

Reimplemented from ProcessingSignalListener.

Definition at line 80 of file CKFToSVDFindlet.cc.

{
  Super::beginEvent();
 
  // If the capacity of a std::vector is very large without being used, it just allocates RAM for no reason, increasing the RAM
  // usage unnecessarily. In this case, start with a fresh one.
  // Since std::vector.shrink() or std::vector.shrink_to_fit() not necessarily reduce the capacity in the desired way, create a
  // temporary vector of the same type, swap them to use the vector at the new location afterwards, and clear the tempoary vector.
  m_cdcRecoTrackVector.clear();
  checkResizeClear<const SpacePoint*>(m_spacePointVector, 40000);
 
  m_seedStates.clear();
  checkResizeClear<CKFToSVDState>(m_states, 40000);
 
  checkResizeClear<TrackFindingCDC::WeightedRelation<CKFToSVDState>>(m_relations, 100000);
 
  m_results.clear();
  m_filteredResults.clear();
}

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();
  }
}

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();
}

exposeParameters()

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

overridevirtual

Expose the parameters of the sub findlets.

Reimplemented from Findlet<>.

Definition at line 44 of file CKFToSVDFindlet.cc.

{
  Super::exposeParameters(moduleParamList, prefix);
 
  m_dataHandler.exposeParameters(moduleParamList, prefix);
  m_hitsLoader.exposeParameters(moduleParamList, prefix);
  m_stateCreatorFromTracks.exposeParameters(moduleParamList, prefix);
  m_stateCreatorFromHits.exposeParameters(moduleParamList, prefix);
  m_relationCreator.exposeParameters(moduleParamList, prefix);
  m_treeSearchFindlet.exposeParameters(moduleParamList, prefix);
  m_overlapResolver.exposeParameters(moduleParamList, prefix);
  m_spacePointTagger.exposeParameters(moduleParamList, prefix);
  m_resultStorer.exposeParameters(moduleParamList, prefix);
 
  moduleParamList->addParameter("minimalHitRequirement", m_param_minimalHitRequirement,
                                "Minimal Hit requirement for the results (counted in space points)",
                                m_param_minimalHitRequirement);
 
  moduleParamList->getParameter<std::string>("firstHighFilter").setDefaultValue("mva_with_direction_check");
  moduleParamList->getParameter<std::string>("advanceHighFilter").setDefaultValue("advance");
  moduleParamList->getParameter<std::string>("secondHighFilter").setDefaultValue("mva");
  moduleParamList->getParameter<std::string>("updateHighFilter").setDefaultValue("fit");
  moduleParamList->getParameter<std::string>("thirdHighFilter").setDefaultValue("mva");
 
  moduleParamList->getParameter<std::string>("filter").setDefaultValue("mva");
 
  moduleParamList->getParameter<std::string>("hitFilter").setDefaultValue("sensor");
  moduleParamList->getParameter<std::string>("seedFilter").setDefaultValue("all");
  moduleParamList->getParameter<std::string>("preSeedFilter").setDefaultValue("loose");
  moduleParamList->getParameter<std::string>("preHitFilter").setDefaultValue("loose");
 
  moduleParamList->getParameter<std::string>("hitsSpacePointsStoreArrayName").setDefaultValue("SVDSpacePoints");
 
  moduleParamList->getParameter<bool>("useAssignedHits").setDefaultValue(false);
}

getDescription()

virtual std::string getDescription ( )

inlinevirtualinherited

Brief description of the purpose of the concrete findlet.

Reimplemented in CosmicsTrackMergerFindlet, TrackFinderCosmics, TrackFinder, TrackFinderAutomaton, and FindletStoreArrayInput< TWrappedFindlet >.

Definition at line 60 of file Findlet.h.

      {
        return "(no description)";
      }

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 ( )

overridevirtualinherited

Receive and dispatch signal before the start of the event processing.

Reimplemented from ProcessingSignalListener.

Reimplemented in UnionVarSet< AObject >, UnionVarSet< Object >, VariadicUnionVarSet< AVarSets >, ResultStorer< Belle2::CKFToPXDResult >, ResultStorer< Belle2::CKFToSVDResult >, BaseEventTimeExtractor< RecoTrack * >, BaseEventTimeExtractor< TrackFindingCDC::CDCWireHit & >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::HyperHough >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::QuadraticLegendre >, StereoHitTrackQuadTreeMatcher< Belle2::TrackFindingCDC::Z0TanLambdaLegendre >, OnVarSet< Filter< ATruthVarSet::Object > >, OnVarSet< Filter< AVarSet::Object > >, OnVarSet< BaseFacetFilter >, OnVarSet< BaseFacetRelationFilter >, OnVarSet< BaseAxialSegmentPairFilter >, OnVarSet< BaseSegmentRelationFilter >, OnVarSet< BaseTrackRelationFilter >, OnVarSet< BaseSegmentPairRelationFilter >, MCSymmetric< BaseAxialSegmentPairFilter >, MCSymmetric< BaseFacetFilter >, MCSymmetric< BaseFacetRelationFilter >, MCSymmetric< BaseSegmentPairFilter >, MCSymmetric< BaseSegmentPairRelationFilter >, MCSymmetric< BaseSegmentRelationFilter >, MCSymmetric< BaseSegmentTripleFilter >, MCSymmetric< BaseSegmentTripleRelationFilter >, MCSymmetric< BaseTrackRelationFilter >, StoreArrayLoader< const Belle2::SpacePoint >, StoreArrayLoader< DataStoreInputTypeRefType >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCFacet >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCWireHit, true >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegment2D >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCTrack >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentPair >, StoreVectorSwapper< Belle2::TrackFindingCDC::CDCSegmentTriple >, RelationVarSet< ABaseVarSet >, QualityIndicatorFilter, TwoHitVirtualIPQIFilter, MultiHoughSpaceFastInterceptFinder, RawTrackCandCleaner< AHit >, RawTrackCandCleaner< Belle2::vxdHoughTracking::VXDHoughState >, RecoTrackStorer, ROIFinder, SingleHoughSpaceFastInterceptFinder, SpacePointLoaderAndPreparer, TrackCandidateOverlapResolver, and TrackCandidateResultRefiner.

Definition at line 15 of file CompositeProcessingSignalListener.cc.

{
  Super::initialize();
  for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
    psl->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_cdcRecoTrackVector

std::vector<RecoTrack*> m_cdcRecoTrackVector

private

Pointers to the CDC Reco tracks as a vector.

Definition at line 103 of file CKFToSVDFindlet.h.

m_dataHandler

TrackLoader m_dataHandler

private

Findlet for retrieving the cdc tracks and writing the result out.

Definition at line 83 of file CKFToSVDFindlet.h.

m_filteredResults

std::vector<CKFToSVDResult> m_filteredResults

private

Vector for storing the filtered results.

Definition at line 115 of file CKFToSVDFindlet.h.

m_hitsLoader

SpacePointLoader m_hitsLoader

private

Findlet for loading the space points.

Definition at line 85 of file CKFToSVDFindlet.h.

m_initialized

bool m_initialized = false

privateinherited

Flag to keep track whether initialization happened 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_overlapResolver

OverlapResolver<ChooseableSVDResultFilter> m_overlapResolver

private

Findlet for resolving overlaps.

Definition at line 95 of file CKFToSVDFindlet.h.

m_param_minimalHitRequirement

unsigned int m_param_minimalHitRequirement = 1

private

Minimal hit requirement for the results (counted in number of space points)

Definition at line 79 of file CKFToSVDFindlet.h.

m_relationCreator

CKFRelationCreator<CKFToSVDState, ChooseableSVDRelationFilter> m_relationCreator

private

Findlet for creating relations between states.

Definition at line 91 of file CKFToSVDFindlet.h.

m_relations

std::vector<TrackFindingCDC::WeightedRelation<CKFToSVDState> > m_relations

private

Relations between states.

Definition at line 111 of file CKFToSVDFindlet.h.

m_results

std::vector<CKFToSVDResult> m_results

private

Vector for storing the results.

Definition at line 113 of file CKFToSVDFindlet.h.

m_resultStorer

ResultStorer<CKFToSVDResult> m_resultStorer

private

Findlet for storing the results.

Definition at line 99 of file CKFToSVDFindlet.h.

m_seedStates

std::vector<CKFToSVDState> m_seedStates

private

States for the tracks.

Definition at line 107 of file CKFToSVDFindlet.h.

m_spacePointTagger

SpacePointTagger<CKFToSVDResult, SVDCluster> m_spacePointTagger

private

Findlet for tagging the used space points.

Definition at line 97 of file CKFToSVDFindlet.h.

m_spacePointVector

std::vector<const SpacePoint*> m_spacePointVector

private

Pointers to the (const) SpacePoints as a vector.

Definition at line 105 of file CKFToSVDFindlet.h.

m_stateCreatorFromHits

StateCreator<const SpacePoint, CKFToSVDState> m_stateCreatorFromHits

private

Findlet for creating states out of hits.

Definition at line 89 of file CKFToSVDFindlet.h.

m_stateCreatorFromTracks

StateCreatorWithReversal<CKFToSVDState> m_stateCreatorFromTracks

private

Findlet for creating states out of tracks.

Definition at line 87 of file CKFToSVDFindlet.h.

m_states

std::vector<CKFToSVDState> m_states

private

States for the hits.

Definition at line 109 of file CKFToSVDFindlet.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 happened before.

Definition at line 55 of file ProcessingSignalListener.h.

m_treeSearchFindlet

TreeSearcher<CKFToSVDState, SVDStateRejecter, CKFToSVDResult> m_treeSearchFindlet

private

Findlet doing the main work: the tree finding.

Definition at line 93 of file CKFToSVDFindlet.h.

---

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

• tracking/ckf/svd/findlets/include/CKFToSVDFindlet.h
• tracking/ckf/svd/findlets/src/CKFToSVDFindlet.cc