Belle II Software development
SegmentAliasResolver Class Referenceabstract

Resolves between the potential alias versions of the segments and contained hits. More...

#include <SegmentAliasResolver.h>

Inheritance diagram for SegmentAliasResolver:
Findlet< TrackingUtilities::CDCSegment2D & > CompositeProcessingSignalListener ProcessingSignalListener

Public Types

using IOTypes
 Types that should be served to apply on invocation.
 
using IOTypes
 Types that should be served to apply on invocation.
 
using IOVectors
 Vector types that should be served to apply on invocation.
 
using IOVectors
 Vector types that should be served to apply on invocation.
 

Public Member Functions

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
 Signals the beginning of the event processing.
 
void apply (std::vector< TrackingUtilities::CDCSegment2D > &outputSegments) final
 Main algorithm applying the fit to each segment.
 
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
 Short hand for ToRangeImpl.
 
using ToVector
 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 = TrackingUtilities::Findlet<TrackingUtilities::CDCSegment2D&>
 Type of the base class.
 

Private Member Functions

void refit (TrackingUtilities::CDCSegment2D &segment, bool reestimate)
 Fit the alias segment.
 

Private Attributes

std::vector< std::string > m_param_investigateAlias = { "cross", "borders", }
 Parameter : Which alias resolutions should be applied.
 
bool m_fullAlias = false
 Switch whether the complete segment should be aliased.
 
bool m_crossAliases = false
 Switch whether the segment should be aliased at the waist where an rl side switch was missed.
 
bool m_borderAliases = false
 Switch whether the border hits at the beginning and the end should be aliased.
 
bool m_middleAliases = false
 Switch whether the hits in the middle with lowest drift radius should be aliased.
 
bool m_param_reestimateDriftLength = true
 Parameter : Switch to reestimate the drift length before each fit.
 
bool m_param_reestimatePositions = false
 Parameter : Switch to restimate the reconstructed positions with a pass of facets.
 
CDCRiemannFitter m_riemannFitter
 Instance of the riemann fitter to be used.
 
DriftLengthEstimator m_driftLengthEstimator
 Instance of the drift length estimator to be used.
 
std::vector< ProcessingSignalListener * > m_subordinaryProcessingSignalListeners
 References to subordinary signal processing listener contained in this findlet.
 
bool m_initialized
 Flag to keep track whether initialization happened before.
 
bool m_terminated
 Flag to keep track whether termination happened before.
 
std::string m_initializedAs
 Name of the type during initialisation.
 

Detailed Description

Resolves between the potential alias versions of the segments and contained hits.

Definition at line 31 of file SegmentAliasResolver.h.

Member Typedef Documentation

◆ IOTypes [1/2]

using IOTypes
inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOTypes [2/2]

using IOTypes
inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOVectors [1/2]

using IOVectors
inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ IOVectors [2/2]

using IOVectors
inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ Super

Type of the base class.

Definition at line 35 of file SegmentAliasResolver.h.

◆ ToVector [1/2]

using ToVector
protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

◆ ToVector [2/2]

using ToVector
protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Member Function Documentation

◆ addProcessingSignalListener()

void addProcessingSignalListener ( ProcessingSignalListener * psl)
protectedinherited

Register a processing signal listener to be notified.

Definition at line 53 of file CompositeProcessingSignalListener.cc.

56{
58}
Interface for a minimal algorithm part that wants to expose some parameters to a module.
Definition Findlet.h:26

◆ apply()

void apply ( std::vector< TrackingUtilities::CDCSegment2D > & outputSegments)
final

Main algorithm applying the fit to each segment.

Definition at line 81 of file SegmentAliasResolver.cc.

82{
83 if (m_fullAlias) {
84 for (CDCSegment2D& segment : outputSegments) {
85 if (segment->hasAliasFlag()) continue; // A full alias has been found in the facet ca.
86 int nRLSwitches = segment.getNRLSwitches();
87 // Sufficiently right left constrained that the alias is already fixed.
88 bool aliasStable = nRLSwitches > 2;
89 if (aliasStable) continue;
90
91 CDCSegment2D aliasSegment = segment.getAlias();
92 refit(aliasSegment, true);
93 swapBetterChi2(segment, aliasSegment);
94 } // end alias loop
95 }
96
97 // Detect whether the segment has a waist and flip the rest of the segment to right
98 if (m_crossAliases) {
99 for (CDCSegment2D& segment : outputSegments) {
100 int nRLSwitches = segment.getNRLSwitches();
101 // Sufficiently right left constrained that the alias is already fixed.
102 bool aliasStable = nRLSwitches > 2;
103 if (aliasStable) continue;
104
105 if (nRLSwitches == 1) {
106 // One RL switch. Try the all left and all right aliases
107 CDCSegment2D rightSegment = segment;
108 CDCSegment2D leftSegment = segment;
109 for (CDCRecoHit2D& recoHit2D : rightSegment) {
110 if (recoHit2D.getRLInfo() == ERightLeft::c_Left) {
111 recoHit2D = recoHit2D.getAlias();
112 }
113 }
114
115 for (CDCRecoHit2D& recoHit2D : leftSegment) {
116 if (recoHit2D.getRLInfo() == ERightLeft::c_Right) {
117 recoHit2D = recoHit2D.getAlias();
118 }
119 }
120 refit(rightSegment, true);
121 refit(leftSegment, true);
122 swapBetterChi2(segment, rightSegment);
123 swapBetterChi2(segment, leftSegment);
124
125 } else if (nRLSwitches == 0) {
126 // No RL switch. Try to introduce one at the smallest drift length
127 CDCSegment2D frontCrossSegment = segment;
128 CDCSegment2D backCrossSegment = segment;
129
130 auto lessDriftLength = [](const CDCRecoHit2D & lhs, const CDCRecoHit2D & rhs) {
131 return lhs.getRefDriftLength() < rhs.getRefDriftLength();
132 };
133
134 // Alias before the minimal drift length
135 {
136 auto itMinLRecoHit2D =
137 std::min_element(frontCrossSegment.begin(), frontCrossSegment.end(), lessDriftLength);
138 for (CDCRecoHit2D& recoHit2D : asRange(frontCrossSegment.begin(), itMinLRecoHit2D)) {
139 recoHit2D = recoHit2D.getAlias();
140 }
141 }
142
143 // Alias after the minimal drift length
144 {
145 auto itMinLRecoHit2D =
146 std::min_element(backCrossSegment.begin(), backCrossSegment.end(), lessDriftLength);
147 for (CDCRecoHit2D& recoHit2D : asRange(itMinLRecoHit2D, backCrossSegment.end())) {
148 recoHit2D = recoHit2D.getAlias();
149 }
150 }
151 refit(frontCrossSegment, true);
152 refit(backCrossSegment, true);
153 swapBetterChi2(segment, frontCrossSegment);
154 swapBetterChi2(segment, backCrossSegment);
155 }
156 }
157 }
158
159 if (m_borderAliases) {
160 for (CDCSegment2D& segment : outputSegments) {
161 // Check aliasing last hit
162 CDCSegment2D aliasSegment = segment;
163 CDCRecoHit2D& aliasHit = aliasSegment.back();
164 aliasHit.reverse();
165 aliasHit.setRecoPos2D(Vector2D(NAN, NAN));
166 if (aliasHit.getRefDriftLength() < 0.2) {
167 refit(aliasSegment, false);
168 Vector2D recoPos2D = aliasSegment.getTrajectory2D().getClosest(aliasHit.getRefPos2D());
169 ERightLeft rlInfo = aliasSegment.getTrajectory2D().isRightOrLeft(aliasHit.getRefPos2D());
170 aliasHit.setRecoPos2D(recoPos2D);
171 aliasHit.setRLInfo(rlInfo);
172 swapBetterChi2(segment, aliasSegment);
173 }
174 }
175
176 for (CDCSegment2D& segment : outputSegments) {
177 // Check aliasing first hit
178 CDCSegment2D aliasSegment = segment;
179 CDCRecoHit2D& aliasHit = aliasSegment.front();
180 aliasHit.reverse();
181 aliasHit.setRecoPos2D(Vector2D(NAN, NAN));
182 if (aliasHit.getRefDriftLength() < 0.2) {
183 refit(aliasSegment, false);
184 Vector2D recoPos2D = aliasSegment.getTrajectory2D().getClosest(aliasHit.getRefPos2D());
185 ERightLeft rlInfo = aliasSegment.getTrajectory2D().isRightOrLeft(aliasHit.getRefPos2D());
186 aliasHit.setRecoPos2D(recoPos2D);
187 aliasHit.setRLInfo(rlInfo);
188 swapBetterChi2(segment, aliasSegment);
189 }
190 }
191 }// end if (m_borderAliases)
192
193 if (m_middleAliases) {
194 for (CDCSegment2D& segment : outputSegments) {
195 // Check aliasing hit with lowest drift length
196 CDCSegment2D aliasSegment = segment;
197
198 // Find hit with minimal drift length
199 double minimalDriftLength = 0.1;
200 CDCRecoHit2D* minimalHit = nullptr;
201 for (CDCRecoHit2D& aliasHit : aliasSegment) {
202 if (aliasHit.getRefDriftLength() < minimalDriftLength) {
203 minimalHit = &aliasHit;
204 minimalDriftLength = aliasHit.getRefDriftLength();
205 }
206 }
207 if (minimalHit) {
208 CDCRecoHit2D& aliasHit = *minimalHit;
209 aliasHit.reverse();
210 aliasHit.setRecoPos2D(Vector2D(NAN, NAN));
211 refit(aliasSegment, false);
212 Vector2D recoPos2D = aliasSegment.getTrajectory2D().getClosest(aliasHit.getRefPos2D());
213 ERightLeft rlInfo = aliasSegment.getTrajectory2D().isRightOrLeft(aliasHit.getRefPos2D());
214 aliasHit.setRecoPos2D(recoPos2D);
215 aliasHit.setRLInfo(rlInfo);
216 swapBetterChi2(segment, aliasSegment);
217 }
218 }
219 }
220
221 std::sort(outputSegments.begin(), outputSegments.end());
222}
bool m_middleAliases
Switch whether the hits in the middle with lowest drift radius should be aliased.
void refit(TrackingUtilities::CDCSegment2D &segment, bool reestimate)
Fit the alias segment.
bool m_fullAlias
Switch whether the complete segment should be aliased.
bool m_borderAliases
Switch whether the border hits at the beginning and the end should be aliased.
bool m_crossAliases
Switch whether the segment should be aliased at the waist where an rl side switch was missed.
void setRLInfo(ERightLeft &rlInfo)
Setter the right left passage information.
void reverse()
Turns the orientation in place.
double getRefDriftLength() const
Getter for the drift length at the wire reference position.
const ROOT::Math::XYVector & getRefPos2D() const
Getter for the reference position of the wire.
void setRecoPos2D(const Vector2D &recoPos2D)
Setter for the position in the reference plane.
CDCTrajectory2D & getTrajectory2D() const
Getter for the two dimensional trajectory fitted to the segment.
Definition CDCSegment.h:69
ERightLeft isRightOrLeft(const Vector2D &point) const
Checks if the given point is to the right or to the left of the trajectory.
Vector2D getClosest(const Vector2D &point) const
Calculates the closest approach on the trajectory to the given point.
ERightLeft
Enumeration to represent the distinct possibilities of the right left passage.
Definition ERightLeft.h:25

◆ beginEvent()

void beginEvent ( )
overrideinherited

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

Definition at line 36 of file CompositeProcessingSignalListener.cc.

32{
35 psl->beginEvent();
36 }
37}
void beginEvent() override
Receive and dispatch signal for the start of a new event.
virtual void beginEvent()
Receive signal for the start of a new event.

◆ beginRun()

void beginRun ( )
overrideinherited

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

Definition at line 33 of file CompositeProcessingSignalListener.cc.

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

◆ endRun()

void endRun ( )
overrideinherited

Receive and dispatch signal for the end of the run.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

40{
42 psl->endRun();
43 }
45}
void endRun() override
Receive and dispatch signal for the end of the run.
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 CompositeProcessingSignalListener.

Definition at line 44 of file SegmentAliasResolver.cc.

45{
46 moduleParamList->addParameter(prefixed(prefix, "investigateAlias"),
48 "Which alias resolutions should be applied. "
49 "Options are trailing_hit, last_hit, full.",
51
52 moduleParamList->addParameter(prefixed(prefix, "reestimateDriftLength"),
54 "Switch to reestimate the drift length in the alias resolver",
56
57 moduleParamList->addParameter(prefixed(prefix, "reestimatePositions"),
59 "Switch to reestimate the reconstructed positions in the alias resolver",
61}
bool m_param_reestimatePositions
Parameter : Switch to restimate the reconstructed positions with a pass of facets.
std::vector< std::string > m_param_investigateAlias
Parameter : Which alias resolutions should be applied.
bool m_param_reestimateDriftLength
Parameter : Switch to reestimate the drift length before each fit.
void addParameter(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module list.

◆ getDescription()

std::string getDescription ( )
finalvirtual

Short description of the findlet.

Reimplemented from Findlet< TrackingUtilities::CDCSegment2D & >.

Definition at line 39 of file SegmentAliasResolver.cc.

40{
41 return "Resolves the rl aliasing of segments in various scenarios";
42}

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )
protectedinherited

Get the number of currently registered listeners.

Definition at line 56 of file CompositeProcessingSignalListener.cc.

61{
63}

◆ initialize()

void initialize ( )
finalvirtual

Signals the beginning of the event processing.

Reimplemented from CompositeProcessingSignalListener.

Definition at line 63 of file SegmentAliasResolver.cc.

64{
66 for (const std::string& investigatedAlias : m_param_investigateAlias) {
67 if (investigatedAlias == "full") {
68 m_fullAlias = true;
69 } else if (investigatedAlias == "cross") {
70 m_crossAliases = true;
71 } else if (investigatedAlias == "borders") {
72 m_borderAliases = true;
73 } else if (investigatedAlias == "middle") {
74 m_middleAliases = true;
75 } else {
76 B2ERROR("Unknown alias to investigate " << investigatedAlias);
77 }
78 }
79}

◆ refit()

void refit ( TrackingUtilities::CDCSegment2D & segment,
bool reestimate )
private

Fit the alias segment.

Definition at line 224 of file SegmentAliasResolver.cc.

225{
226 if (reestimate) {
228 CDCRLWireHitSegment rlWireHitSegment = segment.getRLWireHitSegment();
229 // Reset the drift length to the default values.
230 for (CDCRLWireHit& rlWireHit : rlWireHitSegment) {
231 rlWireHit.setRefDriftLength(rlWireHit.getWireHit().getRefDriftLength());
232 }
233 CDCFacetSegment facetSegment = CDCFacetSegment::create(rlWireHitSegment);
234 for (CDCFacet& facet : facetSegment) {
235 facet.adjustFitLine();
237 m_driftLengthEstimator.updateDriftLength(facet);
238 }
239 }
240 CDCSegment2D replacement = CDCSegment2D::condense(facetSegment);
241 segment = replacement;
242 }
243
245 m_driftLengthEstimator.updateDriftLength(segment);
246 }
247 }
248
249 const EFitPos fitPos = EFitPos::c_RecoPos;
250 const EFitVariance fitVariance = EFitVariance::c_Proper;
251 CDCObservations2D observations2D(fitPos, fitVariance);
252 observations2D.appendRange(segment);
253 if (observations2D.size() < 4) {
254 segment.getTrajectory2D().clear();
255 } else {
256 CDCTrajectory2D trajectory2D = m_riemannFitter.fit(observations2D);
257 segment.setTrajectory2D(trajectory2D);
258 }
259
260}
DriftLengthEstimator m_driftLengthEstimator
Instance of the drift length estimator to be used.
CDCRiemannFitter m_riemannFitter
Instance of the riemann fitter to be used.
static CDCFacetSegment create(const CDCRLWireHitSegment &rlWireHitSegment)
Construct a train of facets from the given oriented wire hits.
static CDCSegment2D condense(const CDCTangentSegment &tangentSegment)
Averages the reconstructed positions from hits that overlap in adjacent tangents in the given tangent...

◆ terminate()

void terminate ( )
overrideinherited

Receive and dispatch Signal for termination of the event processing.

Definition at line 42 of file CompositeProcessingSignalListener.cc.

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

Member Data Documentation

◆ m_borderAliases

bool m_borderAliases = false
private

Switch whether the border hits at the beginning and the end should be aliased.

Definition at line 66 of file SegmentAliasResolver.h.

◆ m_crossAliases

bool m_crossAliases = false
private

Switch whether the segment should be aliased at the waist where an rl side switch was missed.

Definition at line 63 of file SegmentAliasResolver.h.

◆ m_driftLengthEstimator

DriftLengthEstimator m_driftLengthEstimator
private

Instance of the drift length estimator to be used.

Definition at line 81 of file SegmentAliasResolver.h.

◆ m_fullAlias

bool m_fullAlias = false
private

Switch whether the complete segment should be aliased.

Definition at line 60 of file SegmentAliasResolver.h.

◆ m_initialized

bool m_initialized
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_middleAliases

bool m_middleAliases = false
private

Switch whether the hits in the middle with lowest drift radius should be aliased.

Definition at line 69 of file SegmentAliasResolver.h.

◆ m_param_investigateAlias

std::vector<std::string> m_param_investigateAlias = { "cross", "borders", }
private

Parameter : Which alias resolutions should be applied.

Definition at line 57 of file SegmentAliasResolver.h.

57{ /*"full",*/ "cross", "borders", /*"middle"*/};

◆ m_param_reestimateDriftLength

bool m_param_reestimateDriftLength = true
private

Parameter : Switch to reestimate the drift length before each fit.

Definition at line 72 of file SegmentAliasResolver.h.

◆ m_param_reestimatePositions

bool m_param_reestimatePositions = false
private

Parameter : Switch to restimate the reconstructed positions with a pass of facets.

Definition at line 75 of file SegmentAliasResolver.h.

◆ m_riemannFitter

CDCRiemannFitter m_riemannFitter
private

Instance of the riemann fitter to be used.

Definition at line 78 of file SegmentAliasResolver.h.

◆ m_subordinaryProcessingSignalListeners

std::vector<ProcessingSignalListener*> m_subordinaryProcessingSignalListeners
privateinherited

References to subordinary signal processing listener contained in this findlet.

Definition at line 60 of file CompositeProcessingSignalListener.h.

◆ m_terminated

bool m_terminated
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: