Belle II Software development
Public Types | Public Member Functions | Protected Types | Protected Member Functions | Private Types | Private Attributes | List of all members
CDCCKFDuplicateRemover Class Referenceabstract
Packages » tracking

Remove duplicate paths created from ECLShowers These typically come from the seeding with two charge assumptions and Bremsstrahlung. More...

#include <CDCCKFDuplicateRemover.h>

Inheritance diagram for CDCCKFDuplicateRemover:
Findlet< CDCCKFResult > CompositeProcessingSignalListener ProcessingSignalListener

Public Types

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

Public Member Functions

void exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) override
 Expose the parameters of the sub findlets.
 
void apply (std::vector< CDCCKFResult > &results) override
 main method of the findlet, merges and filters paths
 
virtual std::string getDescription ()
 Brief description of the purpose of the concrete findlet.
 
virtual void apply (ToVector< AIOTypes > &... ioVectors)=0
 Main function executing the algorithm.
 
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 Types

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
 Type of the base class.
 

Private Attributes

TrackFindingCDC::ChooseableFilter< CDCPathFilterFactory > m_filter_badTracks
 Filter to remove badly reconstructed tracks (e.g. wrongly assigned charge)
 
TrackFindingCDC::ChooseableFilter< CDCPathPairFilterFactory > m_filter_duplicateTrack
 Filter to remove duplicates from helix extrapolation (2 charge assumptions)
 
TrackFindingCDC::ChooseableFilter< CDCPathPairFilterFactory > m_filter_duplicateSeed
 Merge duplicate paths (mostly seeds from Bremstrahlung)
 
double duplicateSeed_maxPhi = 2.
 Seeds within this dPhi can be considered as duplicates.
 
double duplicateSeed_maxTheta = 0.1
 Seeds within this dTheta can be considered as duplicates.
 
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

Remove duplicate paths created from ECLShowers These typically come from the seeding with two charge assumptions and Bremsstrahlung.

Definition at line 29 of file CDCCKFDuplicateRemover.h.

Member Typedef Documentation

◆ IOTypes

using IOTypes
inherited

Types that should be served to apply on invocation.

Definition at line 30 of file Findlet.h.

◆ IOVectors

using IOVectors
inherited

Vector types that should be served to apply on invocation.

Definition at line 53 of file Findlet.h.

◆ Super

using Super
privateinherited

Type of the base class.

Definition at line 26 of file CompositeProcessingSignalListener.h.

◆ ToVector

using ToVector
protectedinherited

Short hand for ToRangeImpl.

Definition at line 49 of file Findlet.h.

Constructor & Destructor Documentation

◆ CDCCKFDuplicateRemover()

CDCCKFDuplicateRemover ( )
inline

Definition at line 31 of file CDCCKFDuplicateRemover.h.

32 {
33 addProcessingSignalListener(&m_filter_badTracks);
34 addProcessingSignalListener(&m_filter_duplicateTrack);
35 addProcessingSignalListener(&m_filter_duplicateSeed);
36 }

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{
57 m_subordinaryProcessingSignalListeners.push_back(psl);
58}
Belle2::TrackFindingCDC::Findlet
Interface for a minimal algorithm part that wants to expose some parameters to a module.
Definition Findlet.h:26

◆ apply()

void apply ( std::vector< CDCCKFResult > & results)
inlineoverride

main method of the findlet, merges and filters paths

Definition at line 57 of file CDCCKFDuplicateRemover.h.

58 {
59 B2DEBUG(29, "CDCCKFDuplicateRemover: " << results.size() << " paths created (might be without any hits)");
60
61 std::vector<CDCCKFResult> goodResults;
62
63 // Additional filter (typically check if charge of reconstructed track is equal to charge of seed)
64 TrackFindingCDC::Weight weight;
65 std::unordered_map<double, CDCCKFResult> resultToWeightList;
66 for (const auto& result : results) {
67 weight = m_filter_badTracks(result);
68 if (not std::isnan(weight)) {
69 goodResults.push_back(result);
70 }
71 }
72
73 int n_goodresults = goodResults.size();
74
75 B2DEBUG(29, "CDCCKFDuplicateRemover: " << goodResults.size() << " paths created (after filtering)");
76
77 if (n_goodresults > 1) {
78 for (const auto& result : goodResults) {
79 B2DEBUG(29, "charge = " << result.front().getSeed()->getChargeSeed() << "; "
80 << "theta = " << result.front().getSeed()->getPositionSeed().Theta() * 180. / M_PI << "; "
81 << (result.size() - 1) << " hits (" << result.at(1).getWireHit()->getWire().getICLayer() << "->" <<
82 result.back().getWireHit()->getWire().getICLayer() << "); "
83 << "r/z = " << result.front().getSeed()->getPositionSeed().Rho() << "; " << result.front().getSeed()->getPositionSeed().Z());
84 }
85 }
86
87 // If both charge assumptions lead to a good track, only pick one to avoid duplicate tracks
88 auto iter = goodResults.begin();
89 while (iter < goodResults.end()) {
90 auto iter2 = iter + 1;
91 bool increaseIter = true;
92 while (iter2 < goodResults.end()) {
93 // find tracks from same seed
94 if (iter2->front().getSeed()->getRelated<ECLShower>() == iter2->front().getSeed()->getRelated<ECLShower>()) {
95 // let filter decide which one to keep
96 bool selectFirst = m_filter_duplicateTrack(std::pair(&*iter, &*iter2)) > 0;
97 if (selectFirst) {
98 iter2 = goodResults.erase(iter2);
99 } else {
100 iter = goodResults.erase(iter);
101 increaseIter = false;
102 break;
103 }
104 } else {
105 ++iter2;
106 }
107 }
108 if (increaseIter) {
109 ++iter;
110 }
111 }
112
113 B2DEBUG(29, "CDCCKFDuplicateRemover: " << goodResults.size() << " paths created (after duplicates)");
114
115 if (n_goodresults > 1) {
116 for (const auto& result : goodResults) {
117 B2DEBUG(29, "charge = " << result.front().getSeed()->getChargeSeed() << "; "
118 << "theta = " << result.front().getSeed()->getPositionSeed().Theta() * 180. / M_PI << "; "
119 << (result.size() - 1) << " hits (" << result.at(1).getWireHit()->getWire().getICLayer() << "->" <<
120 result.back().getWireHit()->getWire().getICLayer() << "); "
121 << "r/z = " << result.front().getSeed()->getPositionSeed().Rho() << "; " << result.front().getSeed()->getPositionSeed().Z());
122 }
123 }
124
125 // Remove duplicate tracks from Bremsstrahlung
126 // Be careful as this might also remove photon conversions (m_filter_duplicateSeed decides if both should be kept)
127 iter = goodResults.begin();
128 while (iter < goodResults.end()) {
129 double phiClus = iter->front().getSeed()->getPositionSeed().Phi();
130 double thetaClus = iter->front().getSeed()->getPositionSeed().Theta();
131
132 auto iter2 = iter + 1;
133 bool increaseIter = true;
134 while (iter2 < goodResults.end()) {
135 // find tracks from close-by seeds (use small strip in phi direction as expected from Bremsstrahlung)
136 // to disregard this filter set duplicateSeed_maxPhi and duplicateSeed_maxTheta to negative values
137 if (std::abs(ROOT::Math::VectorUtil::Phi_mpi_pi(iter2->front().getSeed()->getPositionSeed().Phi() - phiClus)) < duplicateSeed_maxPhi
138 && std::abs(iter2->front().getSeed()->getPositionSeed().Theta() - thetaClus) < duplicateSeed_maxTheta) {
139 // let filter decide which one to keep
140 bool isDuplicate = m_filter_duplicateSeed(std::pair(&*iter, &*iter2)) > 0;
141 if (! isDuplicate) {
142 B2DEBUG(29, "Keeping both tracks");
143 ++iter2;
144 } else {
145 B2DEBUG(29, "Duplicate hits found");
146 bool selectFirst = m_filter_duplicateTrack(std::pair(&*iter, &*iter2)) > 0;
147 if (selectFirst) {
148 iter2 = goodResults.erase(iter2);
149 } else {
150 iter = goodResults.erase(iter);
151 increaseIter = false;
152 break;
153 }
154 }
155 } else {
156 ++iter2;
157 }
158 }
159 if (increaseIter) {
160 ++iter;
161 }
162 }
163
164 results = goodResults;
165
166 B2DEBUG(29, "CDCCKFDuplicateRemover: " << results.size() << " paths created (after merging)");
167
168 if (n_goodresults > 1) {
169 for (const auto& result : results) {
170 B2DEBUG(29, "charge = " << result.front().getSeed()->getChargeSeed() << "; "
171 << "theta = " << result.front().getSeed()->getPositionSeed().Theta() * 180. / M_PI << "; "
172 << (result.size() - 1) << " hits (" << result.at(1).getWireHit()->getWire().getICLayer() << "->" <<
173 result.back().getWireHit()->getWire().getICLayer() << "); "
174 << "r/z = " << result.front().getSeed()->getPositionSeed().Rho() << "; " << result.front().getSeed()->getPositionSeed().Z());
175 }
176 }
177 }

◆ beginEvent()

void beginEvent ( )
overrideinherited

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

Definition at line 36 of file CompositeProcessingSignalListener.cc.

32{
33 Super::beginEvent();
34 for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
35 psl->beginEvent();
36 }
37}
Belle2::TrackFindingCDC::CompositeProcessingSignalListener::beginEvent
void beginEvent() override
Receive and dispatch signal for the start of a new event.
Definition CompositeProcessingSignalListener.cc:31
Belle2::TrackFindingCDC::ProcessingSignalListener::beginEvent
virtual void beginEvent()
Receive signal for the start of a new event.
Definition ProcessingSignalListener.cc:45

◆ beginRun()

void beginRun ( )
overrideinherited

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

Definition at line 33 of file CompositeProcessingSignalListener.cc.

24{
25 Super::beginRun();
26 for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
27 psl->beginRun();
28 }
29}
Belle2::TrackFindingCDC::CompositeProcessingSignalListener::beginRun
void beginRun() override
Receive and dispatch signal for the beginning of a new run.
Definition CompositeProcessingSignalListener.cc:23
Belle2::TrackFindingCDC::ProcessingSignalListener::beginRun
virtual void beginRun()
Receive signal for the beginning of a new run.
Definition ProcessingSignalListener.cc:38

◆ endRun()

void endRun ( )
overrideinherited

Receive and dispatch signal for the end of the run.

Definition at line 39 of file CompositeProcessingSignalListener.cc.

40{
41 for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
42 psl->endRun();
43 }
44 Super::endRun();
45}
Belle2::TrackFindingCDC::CompositeProcessingSignalListener::endRun
void endRun() override
Receive and dispatch signal for the end of the run.
Definition CompositeProcessingSignalListener.cc:39
Belle2::TrackFindingCDC::ProcessingSignalListener::endRun
virtual void endRun()
Receive signal for the end of the run.
Definition ProcessingSignalListener.cc:52

◆ exposeParameters()

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

Expose the parameters of the sub findlets.

Reimplemented from CompositeProcessingSignalListener.

Definition at line 39 of file CDCCKFDuplicateRemover.h.

40 {
41 m_filter_badTracks.exposeParameters(moduleParamList, TrackFindingCDC::prefixed("badTracks", prefix));
42 m_filter_duplicateTrack.exposeParameters(moduleParamList, TrackFindingCDC::prefixed("duplicateTrack", prefix));
43 m_filter_duplicateSeed.exposeParameters(moduleParamList, TrackFindingCDC::prefixed("duplicateSeed", prefix));
44
45 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "duplicateSeed_maxPhi"),
46 duplicateSeed_maxPhi,
47 "Seeds within this dPhi can be considered as duplicates (-1 to neglect)",
48 duplicateSeed_maxPhi);
49
50 moduleParamList->addParameter(TrackFindingCDC::prefixed(prefix, "duplicateSeed_maxTheta"),
51 duplicateSeed_maxTheta,
52 "Seeds within this dTheta can be considered as duplicates (-1 to neglect)",
53 duplicateSeed_maxTheta);
54 }

◆ getDescription()

virtual std::string getDescription ( )
inlinevirtualinherited

Brief description of the purpose of the concrete findlet.

Definition at line 60 of file Findlet.h.

61 {
62 return "(no description)";
63 }

◆ getNProcessingSignalListener()

int getNProcessingSignalListener ( )
protectedinherited

Get the number of currently registered listeners.

Definition at line 56 of file CompositeProcessingSignalListener.cc.

61{
62 return m_subordinaryProcessingSignalListeners.size();
63}

◆ initialize()

void initialize ( )
overrideinherited

Receive and dispatch signal before the start of the event processing.

Definition at line 30 of file CompositeProcessingSignalListener.cc.

16{
17 Super::initialize();
18 for (ProcessingSignalListener* psl : m_subordinaryProcessingSignalListeners) {
19 psl->initialize();
20 }
21}
Belle2::TrackFindingCDC::CompositeProcessingSignalListener::initialize
void initialize() override
Receive and dispatch signal before the start of the event processing.
Definition CompositeProcessingSignalListener.cc:15
Belle2::TrackFindingCDC::ProcessingSignalListener::initialize
virtual void initialize()
Receive signal before the start of the event processing.
Definition ProcessingSignalListener.cc:28

◆ terminate()

void terminate ( )
overrideinherited

Receive and dispatch Signal for termination of the event processing.

Definition at line 42 of file CompositeProcessingSignalListener.cc.

48{
49 for (ProcessingSignalListener* psl : reversedRange(m_subordinaryProcessingSignalListeners)) {
50 psl->terminate();
51 }
52 Super::terminate();
53}
Belle2::TrackFindingCDC::CompositeProcessingSignalListener::terminate
void terminate() override
Receive and dispatch Signal for termination of the event processing.
Definition CompositeProcessingSignalListener.cc:47
Belle2::TrackFindingCDC::ProcessingSignalListener::terminate
virtual void terminate()
Receive Signal for termination of the event processing.
Definition ProcessingSignalListener.cc:59

Member Data Documentation

◆ duplicateSeed_maxPhi

double duplicateSeed_maxPhi = 2.
private

Seeds within this dPhi can be considered as duplicates.

Definition at line 188 of file CDCCKFDuplicateRemover.h.

◆ duplicateSeed_maxTheta

double duplicateSeed_maxTheta = 0.1
private

Seeds within this dTheta can be considered as duplicates.

Definition at line 190 of file CDCCKFDuplicateRemover.h.

◆ m_filter_badTracks

TrackFindingCDC::ChooseableFilter<CDCPathFilterFactory> m_filter_badTracks
private

Filter to remove badly reconstructed tracks (e.g. wrongly assigned charge)

Definition at line 181 of file CDCCKFDuplicateRemover.h.

◆ m_filter_duplicateSeed

TrackFindingCDC::ChooseableFilter<CDCPathPairFilterFactory> m_filter_duplicateSeed
private

Merge duplicate paths (mostly seeds from Bremstrahlung)

Definition at line 185 of file CDCCKFDuplicateRemover.h.

◆ m_filter_duplicateTrack

TrackFindingCDC::ChooseableFilter<CDCPathPairFilterFactory> m_filter_duplicateTrack
private

Filter to remove duplicates from helix extrapolation (2 charge assumptions)

Definition at line 183 of file CDCCKFDuplicateRemover.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_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 file: