Belle II Software  release-08-01-10
SimpleSegmentTripleFilter Class Reference

Filter for the constuction of segment triples based on simple criterions. More...

#include <SimpleSegmentTripleFilter.h>

Inheritance diagram for SimpleSegmentTripleFilter:
Collaboration diagram for SimpleSegmentTripleFilter:

Public Types

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

Public Member Functions

Weight operator() (const CDCSegmentTriple &segmentTriple) final
 Checks if a triple of axial, stereo and axial segments is a good combination to be stored. More...
 
virtual void exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix)
 Expose the set of parameters of the filter to the module parameter list. More...
 
virtual bool needsTruthInformation ()
 Indicates if the filter requires Monte Carlo information.
 
virtual Weight operator() (const Object &obj)
 Function to evaluate the object. More...
 
Weight operator() (const Object *obj)
 Function to evaluate the object. More...
 
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 = Filter< CDCSegmentTriple >
 Type of the base class.
 

Private Member Functions

const CDCRiemannFittergetFitter2D () const
 Returns the xy fitter instance that is used by this filter.
 
const CDCSZFittergetSZFitter () const
 Returns the sz fitter instance that is used by this filter.
 

Private Attributes

CDCSZFitter m_szFitter
 Memory of the SZ fitter fitting sz lines to the stereo segments.
 
CDCRiemannFitter m_riemannFitter
 Memory of the Riemann fitter for the circle fits.
 
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

Filter for the constuction of segment triples based on simple criterions.

Definition at line 22 of file SimpleSegmentTripleFilter.h.

Member Function Documentation

◆ exposeParameters()

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

Expose the set of parameters of the filter to the module parameter list.

Note that not all filters have yet exposed their parameters in this way.

This method is deprecated as the exposeParams below uses a less compile heavy equivalent.

Reimplemented in MCSymmetric< BaseSegmentTripleFilter >, and MCSegmentTripleFilter.

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

◆ operator()() [1/3]

Weight operator() ( const CDCSegmentTriple segmentTriple)
final

Checks if a triple of axial, stereo and axial segments is a good combination to be stored.

Checks the proper alignement and the quality of connection between all three segments. Returns NAN if the connection shall not be made or a finit value be used as the cell weight of the cell to constructed.

Definition at line 23 of file SimpleSegmentTripleFilter.cc.

24 {
25  const CDCAxialSegment2D* ptrStartSegment = segmentTriple.getStartSegment();
26  const CDCStereoSegment2D* ptrMiddleSegment = segmentTriple.getMiddleSegment();
27  const CDCAxialSegment2D* ptrEndSegment = segmentTriple.getEndSegment();
28 
29  assert(ptrStartSegment);
30  assert(ptrMiddleSegment);
31  assert(ptrEndSegment);
32 
33  const CDCAxialSegment2D& startSegment = *ptrStartSegment;
34  const CDCStereoSegment2D& middleSegment = *ptrMiddleSegment;
35  const CDCAxialSegment2D& endSegment = *ptrEndSegment;
36 
37  const double toleranceFraction = 0.0;
38 
39  // Check if the middle segment lies within the acceptable bounds in angular deviation
40  {
41  //get the remembered fits
42  const CDCTrajectory2D& startFit = startSegment.getTrajectory2D();
43  const CDCTrajectory2D& endFit = endSegment.getTrajectory2D();
44 
45  //use only the first and last hit for this check
46  const CDCRecoHit2D& firstHit = middleSegment.front();
47  const CDCRecoHit2D& lastHit = middleSegment.back();
48 
49  Vector3D firstRecoPos = firstHit.reconstruct3D(startFit);
50  Vector3D lastRecoPos = lastHit.reconstruct3D(endFit);
51 
52  const CDCWire& firstWire = firstHit.getWire();
53  const CDCWire& lastWire = lastHit.getWire();
54 
55  const bool agrees =
56  firstWire.isInCellZBounds(firstRecoPos, toleranceFraction) and
57  lastWire.isInCellZBounds(lastRecoPos, toleranceFraction);
58 
59  if (not agrees) return NAN;
60  }
61 
62  // make more complex judgement on fitness
63 
64  // Get the combined fit of start and end axial segment
65  CDCTrajectory2D trajectory2D = getFitter2D().fit(*(segmentTriple.getStartSegment()),
66  *(segmentTriple.getEndSegment()));
67 
68  // Check if the middle segment is actually coaligned with the trajectory
69  EForwardBackward fbInfo = trajectory2D.isForwardOrBackwardTo(middleSegment);
70  if (fbInfo != EForwardBackward::c_Forward) return NAN;
71 
72  // Reconstruct the middle stereo segment
73  CDCSegment3D reconstructedMiddle;
74  for (const CDCRecoHit2D& recoHit2D : middleSegment) {
75  reconstructedMiddle.push_back(CDCRecoHit3D::reconstruct(recoHit2D, trajectory2D));
76  if (not reconstructedMiddle.back().isInCellZBounds(toleranceFraction)) {
77  B2DEBUG(25, " RecoHit out of CDC");
78  return NAN;
79  }
80  }
81 
82  // Fit the sz slope and intercept
83  CDCTrajectorySZ trajectorySZ;
84  getSZFitter().update(trajectorySZ, middleSegment, trajectory2D);
85  segmentTriple.setTrajectory3D(CDCTrajectory3D(trajectory2D, trajectorySZ));
86 
87  Weight result = startSegment.size() + middleSegment.size() + endSegment.size();
88 
89  return result;
90 }
CDCTrajectory2D fit(const CDCObservations2D &observations2D) const
Fits a collection of observation drift circles.
Class representing a two dimensional reconstructed hit in the central drift chamber.
Definition: CDCRecoHit2D.h:47
Vector3D reconstruct3D(const CDCTrajectory2D &trajectory2D, const double z=0) const
Reconstruct the three dimensional position (especially of stereo hits) by determinating the z coordin...
const CDCWire & getWire() const
Getter for the wire the reconstructed hit assoziated to.
Definition: CDCRecoHit2D.h:175
static CDCRecoHit3D reconstruct(const CDCRecoHit2D &recoHit2D, const CDCTrajectory2D &trajectory2D)
Reconstructs the three dimensional hit from the two dimensional and the two dimensional trajectory.
Definition: CDCRecoHit3D.cc:56
void update(const CDCSegmentPair &segmentPair) const
Updates the trajectory of the axial stereo segment pair inplace.
Definition: CDCSZFitter.cc:163
A reconstructed sequence of two dimensional hits in one super layer.
Definition: CDCSegment2D.h:39
A segment consisting of three dimensional reconstructed hits.
Definition: CDCSegment3D.h:26
const CDCAxialSegment2D * getStartSegment() const
Getter for the start axial segment.
void setTrajectory3D(const CDCTrajectory3D &trajectory3D) const
Setter for the three dimensional helix trajectory.
const CDCAxialSegment2D * getEndSegment() const
Getter for the end axial segment.
const CDCStereoSegment2D * getMiddleSegment() const
Getter for the middle stereo segment.
CDCTrajectory2D & getTrajectory2D() const
Getter for the two dimensional trajectory fitted to the segment.
Definition: CDCSegment.h:69
Particle trajectory as it is seen in xy projection represented as a circle.
EForwardBackward isForwardOrBackwardTo(const AHits &hits) const
Calculates if this trajectory and the hits are coaligned Returns:
Particle full three dimensional trajectory.
Linear trajectory in sz space.
Class representing a sense wire in the central drift chamber.
Definition: CDCWire.h:58
bool isInCellZBounds(const Vector3D &pos3D, const double factor=1) const
Checks whether the position is in the z bounds of the drift cell (scaled by the factor) surrounding t...
Definition: CDCWire.h:295
const CDCSZFitter & getSZFitter() const
Returns the sz fitter instance that is used by this filter.
const CDCRiemannFitter & getFitter2D() const
Returns the xy fitter instance that is used by this filter.
A three dimensional vector.
Definition: Vector3D.h:33
EForwardBackward
Enumeration to represent the distinct possibilities of the right left passage information.

◆ operator()() [2/3]

Weight operator() ( const Object obj)
virtualinherited

Function to evaluate the object.

Base implementation accepts all objects.

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

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

◆ operator()() [3/3]

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 76 of file Filter.icc.h.


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