8#include <tracking/trackFindingCDC/filters/axialSegmentPair/MCAxialSegmentPairFilter.h>
10#include <tracking/trackFindingCDC/mclookup/CDCMCSegment2DLookUp.h>
12#include <tracking/trackFindingCDC/fitting/CDCRiemannFitter.h>
14#include <tracking/trackFindingCDC/eventdata/tracks/CDCAxialSegmentPair.h>
15#include <tracking/trackFindingCDC/eventdata/segments/CDCSegment2D.h>
17#include <tracking/trackFindingCDC/filters/base/MCSymmetricFilter.icc.h>
19#include <framework/core/ModuleParamList.templateDetails.h>
22using namespace TrackFindingCDC;
32 const std::string& prefix)
35 moduleParamList->
addParameter(prefixed(prefix,
"requireRLPure"),
37 "Switch to require the segment combination contain mostly correct rl information",
43 const CDCAxialSegment2D* ptrFromSegment = axialSegmentPair.
getStartSegment();
44 const CDCAxialSegment2D* ptrToSegment = axialSegmentPair.
getEndSegment();
46 assert(ptrFromSegment);
49 const CDCAxialSegment2D& fromSegment = *ptrFromSegment;
50 const CDCAxialSegment2D& toSegment = *ptrToSegment;
55 EForwardBackward pairFBInfo =
60 if (pairFBInfo == EForwardBackward::c_Invalid)
return NAN;
62 if (pairFBInfo == EForwardBackward::c_Forward or (
getAllowReverse() and pairFBInfo == EForwardBackward::c_Backward)) {
65 if (fromNPassedSuperLayers == c_InvalidIndex)
return NAN;
68 if (toNPassedSuperLayers == c_InvalidIndex)
return NAN;
71 if (abs(fromNPassedSuperLayers - toNPassedSuperLayers) == 0)
return NAN;
74 if (abs(fromNPassedSuperLayers - toNPassedSuperLayers) > 2)
return NAN;
81 return fromSegment.size() + toSegment.size();
The Module parameter list class.
Class representing a pair of reconstructed axial segments in adjacent superlayer.
void setTrajectory2D(const CDCTrajectory2D &trajectory2D) const
Setter for the trajectory of the two dimensional trajectory.
const CDCAxialSegment2D * getEndSegment() const
Getter for the end segment.
CDCTrajectory2D & getTrajectory2D() const
Getter for the trajectory of the two dimensional trajectory.
const CDCAxialSegment2D * getStartSegment() const
Getter for the start segment.
CDCTrajectory2D fit(const CDCObservations2D &observations2D) const
Fits a collection of observation drift circles.
Index getLastNPassedSuperLayers(const ACDCHitCollection *ptrHits) const
Getter for number of passed superlayer till the last hit the collection of hits which has the Monte C...
EForwardBackward areAlignedInMCTrackWithRLCheck(const ACDCHitCollection *ptrFromHits, const ACDCHitCollection *ptrToHits) const
Returns if the second collection of hits follows the first collection of hits in their common Monte C...
EForwardBackward areAlignedInMCTrack(const ACDCHitCollection *ptrFromHits, const ACDCHitCollection *ptrToHits) const
Returns if the second collection of hits follows the first collection of hits in their common Monte C...
Index getFirstNPassedSuperLayers(const ACDCHitCollection *ptrHits) const
Getter for number of passed superlayer till the first hit the collection of hits which has the Monte ...
Specialisation of the lookup for the truth values of two dimensional segments.
static const CDCMCSegment2DLookUp & getInstance()
Getter for the singletone instance.
static const CDCRiemannFitter & getFitter()
Static getter for a general Riemann fitter.
Particle trajectory as it is seen in xy projection represented as a circle.
bool isFitted() const
Checks if the circle is already set to a valid value.
bool m_param_requireRLPure
Parameter : Switch to require the segment combination contain mostly correct rl information.
Weight operator()(const CDCAxialSegmentPair &axialSegmentPair) final
Checks if a pair of axial segments is a good combination.
void exposeParameters(ModuleParamList *moduleParamList, const std::string &prefix) final
Expose the set of parameters of the filter to the module parameter list.
MCSymmetric< BaseAxialSegmentPairFilter > Super
Type of the super class.
MCAxialSegmentPairFilter(bool allowReverse=true)
Constructor.
Mixin for filters that use Monte Carlo information.
void exposeParameters(ModuleParamList *moduleParamList, const std::string &prefix) override
bool getAllowReverse() const
void addParameter(const std::string &name, T ¶mVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module list.
Abstract base class for different kinds of events.