Belle II Software development
WholeWireHitRelationFilter Class Reference

Class mapping the neighborhood of wires to the neighborhood of wire hits. More...

#include <WholeWireHitRelationFilter.h>

Inheritance diagram for WholeWireHitRelationFilter:
RelationFilter< CDCWireHit > Filter< AObject > CompositeProcessingSignalListener ProcessingSignalListener

Public Types

using Object = AObject
 Type of the object to be analysed.
 
using Interface = Filter< AObject >
 Mark this class as the basic interface.
 

Public Member Functions

 WholeWireHitRelationFilter (int neighborhoodDegree=2)
 Constructor form the default neighborhood degree.
 
 ~WholeWireHitRelationFilter ()
 Default destructor.
 
void exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) final
 Expose the parameters to a module.
 
std::vector< CDCWireHit * > getPossibleTos (CDCWireHit *from, const std::vector< CDCWireHit * > &wireHits) const final
 Returns a vector containing the neighboring wire hits of the given wire hit out of the sorted range given by the two iterator other arguments.
 
virtual std::vector< CDCWireHit * > getPossibleTos (CDCWireHit *from, const std::vector< CDCWireHit * > &objects) const
 Selects the objects possibly related to the given one from the given pool of objects.
 
virtual Weight operator() (const CDCWireHit &from, const CDCWireHit &to)
 Main filter method returning the weight of the neighborhood relation.
 
virtual std::vector< float > operator() (const std::vector< Relation< CDCWireHit > * > &objs) override
 Filter over a vector of relations.
 
Weight operator() (const Relation< CDCWireHit > &relation) override
 Main filter method overriding the filter interface method.
 
virtual Weight operator() (const Object &obj)
 Function to evaluate the object.
 
Weight operator() (const Object *obj)
 Function to evaluate the object.
 
virtual std::vector< float > operator() (const std::vector< Object * > &objs)
 Function to evaluate a vector of objects Base implementation applies the function to each object.
 
virtual bool needsTruthInformation ()
 Indicates if the filter requires Monte Carlo information.
 
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 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 = ProcessingSignalListener
 Type of the base class.
 

Private Attributes

int m_param_degree = 2
 Degree of the neighbor extend.
 
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 mapping the neighborhood of wires to the neighborhood of wire hits.

Definition at line 26 of file WholeWireHitRelationFilter.h.

Member Typedef Documentation

◆ Interface

using Interface = Filter<AObject>
inherited

Mark this class as the basic interface.

Definition at line 38 of file Filter.dcl.h.

◆ Object

using Object = AObject
inherited

Type of the object to be analysed.

Definition at line 35 of file Filter.dcl.h.

◆ Super

using Super = ProcessingSignalListener
privateinherited

Type of the base class.

Definition at line 25 of file CompositeProcessingSignalListener.h.

Constructor & Destructor Documentation

◆ WholeWireHitRelationFilter()

WholeWireHitRelationFilter ( int  neighborhoodDegree = 2)
explicit

Constructor form the default neighborhood degree.

Definition at line 31 of file WholeWireHitRelationFilter.cc.

32 : m_param_degree(neighborhoodDegree)
33{
34}

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.

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

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

32{
35 psl->beginEvent();
36 }
37}
Interface for an algorithm part that needs to receive the module processing signals.
virtual void beginEvent()
Receive signal for the start of a new event.

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

24{
27 psl->beginRun();
28 }
29}
virtual void beginRun()
Receive signal for the beginning of a new run.

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

40{
42 psl->endRun();
43 }
45}
virtual void endRun()
Receive signal for the end of the run.

◆ exposeParameters()

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

Expose the parameters to a module.

Reimplemented from Filter< AObject >.

Definition at line 38 of file WholeWireHitRelationFilter.cc.

40{
41 moduleParamList
42 ->addParameter(prefixed(prefix, "degree"),
44 "Neighbor degree which are included. 1 for primary, 2 for secondary, 3 ...",
46}
void addParameter(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module list.

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )
protectedinherited

Get the number of currently registered listeners.

Definition at line 60 of file CompositeProcessingSignalListener.cc.

61{
63}

◆ getPossibleTos() [1/2]

std::vector< CDCWireHit * > getPossibleTos ( CDCWireHit from,
const std::vector< CDCWireHit * > &  objects 
) const
virtualinherited

Selects the objects possibly related to the given one from the given pool of objects.

Definition at line 39 of file RelationFilter.icc.h.

30 {
31 return objects;
32 }

◆ getPossibleTos() [2/2]

std::vector< CDCWireHit * > getPossibleTos ( CDCWireHit from,
const std::vector< CDCWireHit * > &  wireHits 
) const
final

Returns a vector containing the neighboring wire hits of the given wire hit out of the sorted range given by the two iterator other arguments.

Definition at line 48 of file WholeWireHitRelationFilter.cc.

51{
52 assert(std::is_sorted(wireHits.begin(), wireHits.end(), LessOf<Deref>()) &&
53 "Expected wire hits to be sorted");
54
55 const int nWireNeighbors = 8 + 10 * (m_param_degree - 1);
56 std::vector<const CDCWire*> m_wireNeighbors;
57 m_wireNeighbors.reserve(nWireNeighbors);
58
59 std::vector<CDCWireHit*> m_wireHitNeighbors;
60 m_wireHitNeighbors.reserve(2 * nWireNeighbors);
61
62 const CDCWireTopology& wireTopology = CDCWireTopology::getInstance();
63 const CDCWire& wire = from->getWire();
64
65 const CDCWire* ccwSixthSecondWireNeighbor = wireTopology.getSecondNeighborSixOClock(wire);
66 const CDCWire* ccwInWireNeighbor = wireTopology.getNeighborCCWInwards(wire);
67 const CDCWire* ccwWireNeighbor = wireTopology.getNeighborCCW(wire);
68 const CDCWire* ccwOutWireNeighbor = wireTopology.getNeighborCCWOutwards(wire);
69 const CDCWire* ccwTwelvethSecondWireNeighbor = wireTopology.getSecondNeighborTwelveOClock(wire);
70
71 const CDCWire* cwSixthSecondWireNeighbor = ccwSixthSecondWireNeighbor;
72 const CDCWire* cwInWireNeighbor = wireTopology.getNeighborCWInwards(wire);
73 const CDCWire* cwWireNeighbor = wireTopology.getNeighborCW(wire);
74 const CDCWire* cwOutWireNeighbor = wireTopology.getNeighborCWOutwards(wire);
75 const CDCWire* cwTwelvethSecondWireNeighbor = ccwTwelvethSecondWireNeighbor;
76
77 // Insert the neighbors such that they are most likely sorted.
78
79 // Degree 1 neighborhood - only add the six o'clock and the twelve o'clock neighbor once
80 if (m_param_degree > 1 and ccwSixthSecondWireNeighbor)
81 m_wireNeighbors.push_back(ccwSixthSecondWireNeighbor);
82
83 if (cwInWireNeighbor) m_wireNeighbors.push_back(cwInWireNeighbor);
84 if (ccwInWireNeighbor) m_wireNeighbors.push_back(ccwInWireNeighbor);
85
86 if (cwWireNeighbor) m_wireNeighbors.push_back(cwWireNeighbor);
87 if (ccwWireNeighbor) m_wireNeighbors.push_back(ccwWireNeighbor);
88
89 if (cwOutWireNeighbor) m_wireNeighbors.push_back(cwOutWireNeighbor);
90 if (ccwOutWireNeighbor) m_wireNeighbors.push_back(ccwOutWireNeighbor);
91
92 if (m_param_degree > 1 and ccwTwelvethSecondWireNeighbor)
93 m_wireNeighbors.push_back(ccwTwelvethSecondWireNeighbor);
94
95 for (int degree = 1; degree < m_param_degree; ++degree) {
96 if (cwSixthSecondWireNeighbor) {
97 cwSixthSecondWireNeighbor = cwSixthSecondWireNeighbor->getNeighborCW();
98 m_wireNeighbors.push_back(cwSixthSecondWireNeighbor);
99 }
100 if (ccwSixthSecondWireNeighbor) {
101 ccwSixthSecondWireNeighbor = ccwSixthSecondWireNeighbor->getNeighborCCW();
102 m_wireNeighbors.push_back(ccwSixthSecondWireNeighbor);
103 }
104
105 if (cwInWireNeighbor) {
106 cwInWireNeighbor = cwInWireNeighbor->getNeighborCW();
107 m_wireNeighbors.push_back(cwInWireNeighbor);
108 }
109 if (ccwInWireNeighbor) {
110 ccwInWireNeighbor = ccwInWireNeighbor->getNeighborCCW();
111 m_wireNeighbors.push_back(ccwInWireNeighbor);
112 }
113
114 if (cwWireNeighbor) {
115 cwWireNeighbor = cwWireNeighbor->getNeighborCW();
116 m_wireNeighbors.push_back(cwWireNeighbor);
117 }
118 if (ccwWireNeighbor) {
119 ccwWireNeighbor = ccwWireNeighbor->getNeighborCCW();
120 m_wireNeighbors.push_back(ccwWireNeighbor);
121 }
122
123 if (cwOutWireNeighbor) {
124 cwOutWireNeighbor = cwOutWireNeighbor->getNeighborCW();
125 m_wireNeighbors.push_back(cwOutWireNeighbor);
126 }
127 if (ccwOutWireNeighbor) {
128 ccwOutWireNeighbor = ccwOutWireNeighbor->getNeighborCCW();
129 m_wireNeighbors.push_back(ccwOutWireNeighbor);
130 }
131
132 if (cwTwelvethSecondWireNeighbor) {
133 cwTwelvethSecondWireNeighbor = cwTwelvethSecondWireNeighbor->getNeighborCW();
134 m_wireNeighbors.push_back(cwTwelvethSecondWireNeighbor);
135 }
136 if (ccwTwelvethSecondWireNeighbor) {
137 ccwTwelvethSecondWireNeighbor = ccwTwelvethSecondWireNeighbor->getNeighborCCW();
138 m_wireNeighbors.push_back(ccwTwelvethSecondWireNeighbor);
139 }
140 }
141
142 std::sort(std::begin(m_wireNeighbors), std::end(m_wireNeighbors));
143
144 for (const CDCWire* neighborWire : m_wireNeighbors) {
145 ConstVectorRange<CDCWireHit*> neighborWireHits{
146 std::equal_range(wireHits.begin(), wireHits.end(), neighborWire, LessOf<Deref>())};
147
148 m_wireHitNeighbors.insert(m_wireHitNeighbors.end(),
149 neighborWireHits.begin(),
150 neighborWireHits.end());
151 }
152
153 return m_wireHitNeighbors;
154}
Class representing the sense wire arrangement in the whole of the central drift chamber.
MayBePtr< const CDCWire > getNeighborCCWInwards(const WireID &wireID) const
Getter for the nearest counterclockwise neighbor in the next layer outwards.
MayBePtr< const CDCWire > getSecondNeighborSixOClock(const WireID &wireID) const
Getter for secondary neighbor at the six o'clock position.
MayBePtr< const CDCWire > getNeighborCCWOutwards(const WireID &wireID) const
Getter for the nearest counterclockwise neighbor in the next layer outwards.
MayBePtr< const CDCWire > getSecondNeighborTwelveOClock(const WireID &wireID) const
Getter for secondary neighbor at the twelve o'clock position.
MayBePtr< const CDCWire > getNeighborCWInwards(const WireID &wireID) const
Getter for the nearest clockwise neighbor in the next layer outwards.
MayBePtr< const CDCWire > getNeighborCCW(const WireID &wireID) const
Getter for the nearest counterclockwise neighbor.
MayBePtr< const CDCWire > getNeighborCWOutwards(const WireID &wireID) const
Getter for the nearest clockwise neighbor in the next layer outwards.
static CDCWireTopology & getInstance()
Getter for the singleton instance of the wire topology.
MayBePtr< const CDCWire > getNeighborCW(const WireID &wireID) const
Getter for the nearest clockwise neighbor.
Class representing a sense wire in the central drift chamber.
Definition: CDCWire.h:58
MayBePtr< const CDCWire > getNeighborCCW() const
Gives the closest neighbor in the counterclockwise direction - always exists.
Definition: CDCWire.cc:159
MayBePtr< const CDCWire > getNeighborCW() const
Gives the closest neighbor in the clockwise direction - always exists.
Definition: CDCWire.cc:164
A pair of iterators usable with the range base for loop.
Definition: Range.h:25
Functor factory turning a binary functor and two functors into a new functor which executes the binar...
Definition: Functional.h:127

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

16{
19 psl->initialize();
20 }
21}
virtual void initialize()
Receive signal before the start of the event processing.

◆ needsTruthInformation()

◆ operator()() [1/6]

Weight operator() ( const CDCWireHit from,
const CDCWireHit to 
)
virtualinherited

Main filter method returning the weight of the neighborhood relation.

Return always returns NAN to reject all segment neighbors.

Definition at line 46 of file RelationFilter.icc.h.

37 {
38 return 1;
39 }

◆ operator()() [2/6]

◆ operator()() [3/6]

Weight operator() ( const Object obj)
inherited

Function to evaluate the object.

Base implementation accepts all objects, except nullptr.

Parameters
objThe object to be accepted or rejected.
Returns
A finit float value if the object is accepted. NAN if the object is rejected. Nullptr is always rejected.

Definition at line 58 of file Filter.icc.h.

59 {
60 return obj ? operator()(*obj) : NAN;
61 }
virtual Weight operator()(const Object &obj)
Function to evaluate the object.
Definition: Filter.icc.h:52

◆ operator()() [4/6]

Weight operator() ( const Relation< CDCWireHit > &  relation)
overrideinherited

Main filter method overriding the filter interface method.

Checks the validity of the pointers in the relation and unpacks the relation to the method implementing the rejection.

Definition at line 69 of file RelationFilter.icc.h.

43 {
44 const AObject* from = relation.getFrom();
45 const AObject* to = relation.getTo();
46
47 if (from == to) return NAN; // Prevent relation to same.
48 if ((from == nullptr) or (to == nullptr)) return NAN;
49 return operator()(*from, *to);
50 }
virtual Weight operator()(const CDCWireHit &from, const CDCWireHit &to)
Main filter method returning the weight of the neighborhood relation.

◆ operator()() [5/6]

std::vector< float > operator() ( const std::vector< Object * > &  objs)
virtualinherited

Function to evaluate a vector of objects Base implementation applies the function to each object.

Can be optimized for MVA filters

Parameters
objsA vector of pointers to objects
Returns
A vector of float or NAN values. See above

Definition at line 64 of file Filter.icc.h.

65 {
66 std::vector<float> out;
67 for (const auto& obj : objs) {
68 out.push_back(operator()(obj));
69 }
70 return out;
71 }

◆ operator()() [6/6]

virtual std::vector< float > operator() ( const std::vector< Relation< CDCWireHit > * > &  objs)
inlineoverridevirtualinherited

Filter over a vector of relations.

We need to override it since derived relation classes do not implement a call with a pointer to a relation.

Definition at line 51 of file RelationFilter.dcl.h.

52 {
53 std::vector<float> out(objs.size());
54 for (size_t iObj = 0; iObj < objs.size(); iObj += 1) {
55 if (objs[iObj]) {
56 out[iObj] = operator()(*objs[iObj]);
57 } else {
58 out[iObj] = NAN;
59 }
60 }
61 return out;
62 }

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

48{
50 psl->terminate();
51 }
53}
virtual void terminate()
Receive Signal for termination of the event processing.

Member Data Documentation

◆ 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_degree

int m_param_degree = 2
private

Degree of the neighbor extend.

Definition at line 47 of file WholeWireHitRelationFilter.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: