SegmentPairCreator Class Reference

Class providing construction combinatorics for the axial stereo segment pairs. More...

#include <SegmentPairCreator.h>

Inheritance diagram for SegmentPairCreator:

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
	SegmentPairCreator ()
	Constructor adding the filter as a subordinary processing signal listener.
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	apply (const std::vector< CDCSegment2D > &inputSegments, std::vector< CDCSegmentPair > &segmentPairs) final
	Main method constructing pairs in adjacent super layers.
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	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 CDCSegment2D, CDCSegmentPair >
	Type of the base class.

Private Member Functions
void	create (const std::vector< const CDCSegment2D * > &fromSegments, const std::vector< const CDCSegment2D * > &toSegments, std::vector< CDCSegmentPair > &segmentPairs)
	Creates segment pairs from a combination of from segments and to segments.

Private Attributes
bool	m_param_axialBridging = false
	Parameter : Switch to enable the search for axial to axial pairs to enable more stable reconstruction of the middle stereo.
AxialSegmentPairCreator	m_axialSegmentPairCreator
	Findlet responsible for the creation of axial axial segment pairs.
ChooseableSegmentPairFilter	m_segmentPairFilter
	The filter to be used for the segment pair generation.
std::array< std::vector< const CDCSegment2D * >, ISuperLayerUtil::c_N >	m_segmentsBySuperLayer
	Structure for the segments grouped by super layer id.
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

Class providing construction combinatorics for the axial stereo segment pairs.

Definition at line 34 of file SegmentPairCreator.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 = Findlet<const CDCSegment2D, CDCSegmentPair>

private

Type of the base class.

Definition at line 38 of file SegmentPairCreator.h.

ToVector

using ToVector = typename ToVectorImpl<T>::Type

protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Constructor & Destructor Documentation

SegmentPairCreator()

SegmentPairCreator

(

)

Constructor adding the filter as a subordinary processing signal listener.

Definition at line 26 of file SegmentPairCreator.cc.

{
  this->addProcessingSignalListener(&m_segmentPairFilter);
  this->addProcessingSignalListener(&m_axialSegmentPairCreator);
}

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< CDCSegment2D > &  inputSegments, std::vector< CDCSegmentPair > &  segmentPairs  )

final

Main method constructing pairs in adjacent super layers.

Definition at line 49 of file SegmentPairCreator.cc.

{
  segmentPairs.reserve(500);
 
  // Group the segments by their super layer id
  for (std::vector<const CDCSegment2D*>& segmentsInSuperLayer : m_segmentsBySuperLayer) {
    segmentsInSuperLayer.clear();
  }
 
  for (const CDCSegment2D& segment : inputSegments) {
    if (segment.empty()) continue;
    ISuperLayer iSuperLayer = segment.front().getISuperLayer();
    if (not ISuperLayerUtil::isInvalid(iSuperLayer)) {
      const CDCSegment2D* ptrSegment = &segment;
      m_segmentsBySuperLayer[iSuperLayer].push_back(ptrSegment);
    }
  }
 
  // Make pairs of closeby superlayers
  for (ISuperLayer iSuperLayer = 0; iSuperLayer < ISuperLayerUtil::c_N; ++iSuperLayer) {
    const std::vector<const CDCSegment2D*>& fromSegments = m_segmentsBySuperLayer[iSuperLayer];
 
    // Make pairs of this superlayer and the superlayer more to the inside
    ISuperLayer iSuperLayerIn = ISuperLayerUtil::getNextInwards(iSuperLayer);
    if (ISuperLayerUtil::isInCDC(iSuperLayerIn)) {
      const std::vector<const CDCSegment2D*>& toSegments = m_segmentsBySuperLayer[iSuperLayerIn];
      create(fromSegments, toSegments, segmentPairs);
    }
 
    // Make pairs of this superlayer and the superlayer more to the outside
    ISuperLayer iSuperLayerOut = ISuperLayerUtil::getNextOutwards(iSuperLayer);
    if (ISuperLayerUtil::isInCDC(iSuperLayerOut)) {
      const std::vector<const CDCSegment2D*>& toSegments = m_segmentsBySuperLayer[iSuperLayerOut];
      create(fromSegments, toSegments, segmentPairs);
    }
  }
 
  std::sort(segmentPairs.begin(), segmentPairs.end());
 
  if (not m_param_axialBridging) return;
 
  // Memory for the axial axial segment pairs.
  std::vector<CDCAxialSegmentPair> axialSegmentPairs;
 
  // Create the axial to axial bridges
  m_axialSegmentPairCreator.apply(inputSegments, axialSegmentPairs);
 
  // Make pairs by splitting the axial segment pairs and using the *common* fit as the 3d reconstruction bases
  for (const CDCAxialSegmentPair& axialSegmentPair : axialSegmentPairs) {
    const CDCSegment2D* startSegment = axialSegmentPair.getStartSegment();
    const CDCSegment2D* endSegment = axialSegmentPair.getEndSegment();
 
    CDCTrajectory2D startCommonTrajectory = axialSegmentPair.getTrajectory2D();
    if (not startCommonTrajectory.isFitted()) continue;
 
    CDCTrajectory2D startTrajectory2D = startSegment->getTrajectory2D();
    startCommonTrajectory.setLocalOrigin(startTrajectory2D.getLocalOrigin());
    startSegment->setTrajectory2D(startCommonTrajectory);
    std::vector<const CDCSegment2D*> startSegments{startSegment};
 
    CDCTrajectory2D endCommonTrajectory = axialSegmentPair.getTrajectory2D();
    CDCTrajectory2D endTrajectory2D = endSegment->getTrajectory2D();
    endCommonTrajectory.setLocalOrigin(endTrajectory2D.getLocalOrigin());
    endSegment->setTrajectory2D(endCommonTrajectory);
    std::vector<const CDCSegment2D*> endSegments{endSegment};
 
    if (startSegment->getISuperLayer() == endSegment->getISuperLayer()) {
      ISuperLayer iSuperLayerCommon = startSegment->getISuperLayer();
 
      ISuperLayer iSuperLayerIn = ISuperLayerUtil::getNextInwards(iSuperLayerCommon);
      if (ISuperLayerUtil::isInCDC(iSuperLayerIn)) {
        const std::vector<const CDCSegment2D*>& middleSegments = m_segmentsBySuperLayer[iSuperLayerIn];
        create(startSegments, middleSegments, segmentPairs);
        create(middleSegments, endSegments, segmentPairs);
      }
 
      ISuperLayer iSuperLayerOut = ISuperLayerUtil::getNextOutwards(iSuperLayerCommon);
      if (ISuperLayerUtil::isInCDC(iSuperLayerOut)) {
        const std::vector<const CDCSegment2D*>& middleSegments = m_segmentsBySuperLayer[iSuperLayerOut];
        create(startSegments, middleSegments, segmentPairs);
        create(middleSegments, endSegments, segmentPairs);
      }
 
    } else {
      ISuperLayer iSuperLayerMiddle =
        (startSegment->getISuperLayer() + endSegment->getISuperLayer()) / 2;
      const std::vector<const CDCSegment2D*>& middleSegments = m_segmentsBySuperLayer[iSuperLayerMiddle];
      create(startSegments, middleSegments, segmentPairs);
      create(middleSegments, endSegments, segmentPairs);
    }
    startSegment->setTrajectory2D(startTrajectory2D);
    endSegment->setTrajectory2D(endTrajectory2D);
  }
 
  auto lessPairAndgreaterWeight = [](const CDCSegmentPair & lhs, const CDCSegmentPair & rhs) {
    if (lhs < rhs) return true;
    if (rhs < lhs) return false;
    return lhs.getAutomatonCell().getCellWeight() > rhs.getAutomatonCell().getCellWeight();
  };
 
  std::sort(segmentPairs.begin(), segmentPairs.end(), lessPairAndgreaterWeight);
 
  // auto samePair = [](const CDCSegmentPair& lhs, const CDCSegmentPair& rhs) {
  //   if (lhs < rhs) return false;
  //   if (rhs < lhs) return false;
  //   return true;
  // };
  // erase_unique(segmentPairs, samePair);
}

beginEvent()

void beginEvent ( )

overridevirtualinherited

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

Reimplemented from ProcessingSignalListener.

Reimplemented in SpacePointTagger< Belle2::CKFToPXDResult, Belle2::PXDCluster >, SpacePointTagger< Belle2::CKFToSVDResult, Belle2::SVDCluster >, BaseEventTimeExtractor< RecoTrack * >, BaseEventTimeExtractor< TrackFindingCDC::CDCWireHit & >, SharingHitsMatcher< Belle2::TrackFindingCDC::CDCTrack, Belle2::TrackFindingCDC::CDCSegment2D >, MCSymmetric< BaseAxialSegmentPairFilter >, MCSymmetric< BaseFacetFilter >, MCSymmetric< BaseFacetRelationFilter >, MCSymmetric< BaseSegmentPairFilter >, MCSymmetric< BaseSegmentPairRelationFilter >, MCSymmetric< BaseSegmentRelationFilter >, MCSymmetric< BaseSegmentTripleFilter >, MCSymmetric< BaseSegmentTripleRelationFilter >, MCSymmetric< BaseTrackRelationFilter >, 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 >, RecoTrackStorer, ROIFinder, and SVDHoughTracking.

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

create()

void create ( const std::vector< const CDCSegment2D * > &  fromSegments, const std::vector< const CDCSegment2D * > &  toSegments, std::vector< CDCSegmentPair > &  segmentPairs  )

private

Creates segment pairs from a combination of from segments and to segments.

Definition at line 160 of file SegmentPairCreator.cc.

{
  CDCSegmentPair segmentPair;
  for (const CDCSegment2D* ptrFromSegment : fromSegments) {
    for (const CDCSegment2D* ptrToSegment : toSegments) {
 
      if (ptrFromSegment == ptrToSegment) continue;
      segmentPair.setSegments(ptrFromSegment, ptrToSegment);
      segmentPair.clearTrajectory3D();
 
      Weight pairWeight = m_segmentPairFilter(segmentPair);
      if (not std::isnan(pairWeight)) {
        segmentPair.getAutomatonCell().setCellWeight(pairWeight);
        segmentPairs.push_back(segmentPair);
      }
    }
  }
}

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  )

finalvirtual

Expose the parameters to a module.

Reimplemented from Findlet< const CDCSegment2D, CDCSegmentPair >.

Definition at line 37 of file SegmentPairCreator.cc.

{
  m_segmentPairFilter.exposeParameters(moduleParamList, prefix);
  m_axialSegmentPairCreator.exposeParameters(moduleParamList, prefixed("Axial", prefix));
 
  moduleParamList->addParameter(prefixed(prefix, "axialBridging"),
                                m_param_axialBridging,
                                "Switch to enable the search for axial to axial pairs "
                                "to enable more stable reconstruction of the middle stereo",
                                m_param_axialBridging);
}

getDescription()

std::string getDescription ( )

finalvirtual

Short description of the findlet.

Reimplemented from Findlet< const CDCSegment2D, CDCSegmentPair >.

Definition at line 32 of file SegmentPairCreator.cc.

{
  return "Creates axial stereo segment pairs from a set of segments filtered by some acceptance criterion";
}

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_axialSegmentPairCreator

AxialSegmentPairCreator m_axialSegmentPairCreator

private

Findlet responsible for the creation of axial axial segment pairs.

Definition at line 66 of file SegmentPairCreator.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_param_axialBridging

bool m_param_axialBridging = false

private

Parameter : Switch to enable the search for axial to axial pairs to enable more stable reconstruction of the middle stereo.

Definition at line 62 of file SegmentPairCreator.h.

m_segmentPairFilter

ChooseableSegmentPairFilter m_segmentPairFilter

private

The filter to be used for the segment pair generation.

Definition at line 69 of file SegmentPairCreator.h.

m_segmentsBySuperLayer

std::array<std::vector<const CDCSegment2D*>, ISuperLayerUtil::c_N> m_segmentsBySuperLayer

private

Structure for the segments grouped by super layer id.

Definition at line 73 of file SegmentPairCreator.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.

---

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

• tracking/trackFindingCDC/findlets/minimal/include/SegmentPairCreator.h
• tracking/trackFindingCDC/findlets/minimal/src/SegmentPairCreator.cc