Belle II Software development
CylinderDistancePXDPairFilter Class Reference

Base filter for CKF PXD states. More...

#include <CylinderDistancePXDPairFilter.h>

Inheritance diagram for CylinderDistancePXDPairFilter:
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

TrackFindingCDC::Weight operator() (const std::pair< const CKFToPXDState *, const CKFToPXDState * > &relation) override
 Return the weight based on azimuthal-angle separation.
 
void exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) override
 Expose the parameters.
 
virtual bool needsTruthInformation ()
 Indicates if the filter requires Monte Carlo information.
 
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.
 
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

bool m_param_isBackwardCKF = true
 Is backward or forward ToPXDCKF?
 
double m_param_CylinderExtrapolationToHitPhiCut = 0.1
 Filter potential relations in phi between seed states (based on MeasuredStateOnPlane) and hit states.
 
double m_param_CylinderExtrapolationToHitThetaCut = 0.1
 Filter potential relations in theta between seed states (based on MeasuredStateOnPlane) and hit states.
 
double m_param_RecoTrackToHitPhiCut = 0.5
 Filter potential relations in phi between seed states (based on RecoTracks) and hit states.
 
double m_param_RecoTrackToHitThetaCut = 0.5
 Filter potential relations in theta between seed states (based on RecoTracks) and hit states.
 
double m_param_HitHitPhiCut = 0.2
 Filter potential relations in phi between hit states.
 
double m_param_HitHitThetaCut = 0.2
 Filter potential relations in theta between hit states.
 
const std::array< double, 2 > c_PXDLayerRadii = {1.42854, 2.21218}
 PXD layer radii.
 
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

Base filter for CKF PXD states.

Definition at line 19 of file CylinderDistancePXDPairFilter.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.

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 
)
overridevirtual

Expose the parameters.

Reimplemented from Filter< AObject >.

Definition at line 99 of file CylinderDistancePXDPairFilter.cc.

100{
101 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "isBackwardCKF"), m_param_isBackwardCKF,
102 "Does this CKF extrapolate the tracks backward or forward?",
104 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "cylinderExtrapolationToHitPhiCut"),
106 "Cut in phi for the difference between extrapolated position of the Seed to a cylinder with the radius of the ToState layer and the ToState.",
108 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "cylinderExtrapolationToHitThetaCut"),
110 "Cut in theta for the difference between extrapolated position of the Seed to a cylinder with the radius of the ToState layer and the ToState.",
112 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "recoTrackToHitPhiCut"), m_param_RecoTrackToHitPhiCut,
113 "Cut in phi for the difference between RecoTrack information and current hit-based state.",
115 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "recoTrackToHitThetaCut"), m_param_RecoTrackToHitThetaCut,
116 "Cut in theta for the difference between RecoTrack information and current hit-based state.",
118 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "hitHitPhiCut"), m_param_HitHitPhiCut,
119 "Cut in phi between two hit-based states.", m_param_HitHitPhiCut);
120 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "hitHitThetaCut"), m_param_HitHitThetaCut,
121 "Cut in theta between two hit-based states.", m_param_HitHitThetaCut);
122}
double m_param_RecoTrackToHitPhiCut
Filter potential relations in phi between seed states (based on RecoTracks) and hit states.
bool m_param_isBackwardCKF
Is backward or forward ToPXDCKF?
double m_param_CylinderExtrapolationToHitPhiCut
Filter potential relations in phi between seed states (based on MeasuredStateOnPlane) and hit states.
double m_param_CylinderExtrapolationToHitThetaCut
Filter potential relations in theta between seed states (based on MeasuredStateOnPlane) and hit state...
double m_param_HitHitThetaCut
Filter potential relations in theta between hit states.
double m_param_HitHitPhiCut
Filter potential relations in phi between hit states.
double m_param_RecoTrackToHitThetaCut
Filter potential relations in theta between seed states (based on RecoTracks) and hit states.
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}

◆ 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/4]

◆ operator()() [2/4]

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()() [3/4]

TrackFindingCDC::Weight operator() ( const std::pair< const CKFToPXDState *, const CKFToPXDState * > &  relation)
override

Return the weight based on azimuthal-angle separation.

Definition at line 24 of file CylinderDistancePXDPairFilter.cc.

25{
26 const CKFToPXDState& fromState = *(relation.first);
27 const CKFToPXDState& toState = *(relation.second);
28
29 const CKFToPXDState::stateCache& fromStateCache = fromState.getStateCache();
30 const CKFToPXDState::stateCache& toStateCache = toState.getStateCache();
31
32 B2ASSERT("You have filled the wrong states into this!", toStateCache.isHitState);
33
34 if (not fromStateCache.isHitState) {
35 // We are coming from an SVD / CDC-SVD track, so we can use its position to only look for matching ladders
36 // RecoTracks were fitted before CKF, so the measured states on planes should exist
37 const RecoTrack* seedRecoTrack = fromState.getSeed();
38 genfit::MeasuredStateOnPlane measuredStateOnPlane =
40
41 bool extrapolationSuccessful = true;
42 try {
43 measuredStateOnPlane.extrapolateToCylinder(c_PXDLayerRadii[toStateCache.geoLayer - 1]);
44 } catch (...) {
45 B2DEBUG(20, "ToPXDCKF CylinderDistancePXDPairFilter: extrapolation to cylinder of layer " << toStateCache.geoLayer <<
46 " with radius of " << c_PXDLayerRadii[toStateCache.geoLayer - 1] << " failed");
47 extrapolationSuccessful = false;
48 }
49
50 if (extrapolationSuccessful) {
51 const ROOT::Math::XYZVector& positionOnCylinder = ROOT::Math::XYZVector(measuredStateOnPlane.getPos());
52
53 float phiDiff = positionOnCylinder.Phi() - toStateCache.phi;
54 while (phiDiff > M_PI) phiDiff -= 2. * M_PI;
55 while (phiDiff < -M_PI) phiDiff += 2. * M_PI;
56 const float thetaDiff = positionOnCylinder.Theta() - toStateCache.theta;
57 if (abs(phiDiff) < static_cast<float>(m_param_CylinderExtrapolationToHitPhiCut) and
58 abs(thetaDiff) < static_cast<float>(m_param_CylinderExtrapolationToHitThetaCut)) {
59 return 1.0;
60 } else {
61 return NAN;
62 }
63 } else {
64 // We don't have have an extrapolated position on the layer of the ToState based on this Seed (RecoTrack), so use simple angular filters.
65 float phiDiff = fromStateCache.phi - toStateCache.phi;
66 while (phiDiff > M_PI) phiDiff -= 2. * M_PI;
67 while (phiDiff < -M_PI) phiDiff += 2. * M_PI;
68 const float thetaDiff = fromStateCache.theta - toStateCache.theta;
69 if (abs(phiDiff) < static_cast<float>(m_param_RecoTrackToHitPhiCut) and
70 abs(thetaDiff) < static_cast<float>(m_param_RecoTrackToHitThetaCut)) {
71 return 1.0;
72 }
73 return NAN;
74 }
75 }
76
77 // On same layer we already know from LayerPXDRelationFilter, that we only deal with overlaps in r-phi.
78 // So it's sufficient here to check for same layer number to accept states in the overlap region.
79 if (fromStateCache.geoLayer == toStateCache.geoLayer and
80 fromStateCache.sensorID.getSensorNumber() == toStateCache.sensorID.getSensorNumber()) {
81 // TODO: Checking for equality of sensor numbers seems not to harm the hit efficiency,
82 // but maybe it's safer to allow for a sensor number difference of 1?
83 return 1.0;
84 }
85
86 float phiDiff = fromStateCache.phi - toStateCache.phi;
87 while (phiDiff > M_PI) phiDiff -= 2. * M_PI;
88 while (phiDiff < -M_PI) phiDiff += 2. * M_PI;
89 const float thetaDiff = fromStateCache.theta - toStateCache.theta;
90
91 if (abs(phiDiff) < static_cast<float>(m_param_HitHitPhiCut) and
92 abs(thetaDiff) < static_cast<float>(m_param_HitHitThetaCut)) {
93 return 1.0;
94 }
95
96 return NAN;
97}
const Seed * getSeed() const
Return the track this state is related to.
Definition: CKFState.h:60
Specialized CKF State for extrapolating into the PXD.
Definition: CKFToPXDState.h:27
const struct stateCache & getStateCache() const
Get the cached data of this state.
Definition: CKFToPXDState.h:48
const std::array< double, 2 > c_PXDLayerRadii
PXD layer radii.
This is the Reconstruction Event-Data Model Track.
Definition: RecoTrack.h:79
const genfit::MeasuredStateOnPlane & getMeasuredStateOnPlaneFromLastHit(const genfit::AbsTrackRep *representation=nullptr) const
Return genfit's MeasuredStateOnPlane for the last hit in a fit useful for extrapolation of measuremen...
Definition: RecoTrack.cc:619
const genfit::MeasuredStateOnPlane & getMeasuredStateOnPlaneFromFirstHit(const genfit::AbsTrackRep *representation=nullptr) const
Return genfit's MeasuredStateOnPlane for the first hit in a fit useful for extrapolation of measureme...
Definition: RecoTrack.cc:605

◆ operator()() [4/4]

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 }

◆ 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

◆ c_PXDLayerRadii

const std::array<double, 2> c_PXDLayerRadii = {1.42854, 2.21218}
private

PXD layer radii.

Definition at line 42 of file CylinderDistancePXDPairFilter.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_CylinderExtrapolationToHitPhiCut

double m_param_CylinderExtrapolationToHitPhiCut = 0.1
private

Filter potential relations in phi between seed states (based on MeasuredStateOnPlane) and hit states.

Definition at line 30 of file CylinderDistancePXDPairFilter.h.

◆ m_param_CylinderExtrapolationToHitThetaCut

double m_param_CylinderExtrapolationToHitThetaCut = 0.1
private

Filter potential relations in theta between seed states (based on MeasuredStateOnPlane) and hit states.

Definition at line 32 of file CylinderDistancePXDPairFilter.h.

◆ m_param_HitHitPhiCut

double m_param_HitHitPhiCut = 0.2
private

Filter potential relations in phi between hit states.

Definition at line 38 of file CylinderDistancePXDPairFilter.h.

◆ m_param_HitHitThetaCut

double m_param_HitHitThetaCut = 0.2
private

Filter potential relations in theta between hit states.

Definition at line 40 of file CylinderDistancePXDPairFilter.h.

◆ m_param_isBackwardCKF

bool m_param_isBackwardCKF = true
private

Is backward or forward ToPXDCKF?

Definition at line 28 of file CylinderDistancePXDPairFilter.h.

◆ m_param_RecoTrackToHitPhiCut

double m_param_RecoTrackToHitPhiCut = 0.5
private

Filter potential relations in phi between seed states (based on RecoTracks) and hit states.

Definition at line 34 of file CylinderDistancePXDPairFilter.h.

◆ m_param_RecoTrackToHitThetaCut

double m_param_RecoTrackToHitThetaCut = 0.5
private

Filter potential relations in theta between seed states (based on RecoTracks) and hit states.

Definition at line 36 of file CylinderDistancePXDPairFilter.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: