Skip to content
Snippets Groups Projects
Commit b779588d authored by Whitney Armstrong's avatar Whitney Armstrong
Browse files

Updated SimpleClustering algo: use minimum energy when seeding clusters.

	modified:   SimpleClustering.cpp
parent 1b995670
No related branches found
No related tags found
No related merge requests found
...@@ -71,6 +71,7 @@ namespace Jug::Reco { ...@@ -71,6 +71,7 @@ namespace Jug::Reco {
std::vector<eic::CalorimeterHit>& remaining_hits = hits_B; std::vector<eic::CalorimeterHit>& remaining_hits = hits_B;
double max_dist = m_maxDistance.value() / mm; double max_dist = m_maxDistance.value() / mm;
double min_energy = m_minModuleEdep.value() / GeV;
eic::CalorimeterHit ref_hit; eic::CalorimeterHit ref_hit;
ref_hit.energy(0.0); ref_hit.energy(0.0);
...@@ -82,7 +83,7 @@ namespace Jug::Reco { ...@@ -82,7 +83,7 @@ namespace Jug::Reco {
} }
debug() << " max_dist = " << max_dist << endmsg; debug() << " max_dist = " << max_dist << endmsg;
bool continue_clustering = true; bool continue_clustering = (ref_hit.energy() > min_energy);
while (continue_clustering) { while (continue_clustering) {
...@@ -136,6 +137,8 @@ namespace Jug::Reco { ...@@ -136,6 +137,8 @@ namespace Jug::Reco {
std::swap( remaining_hits, the_hits); std::swap( remaining_hits, the_hits);
remaining_hits.clear(); remaining_hits.clear();
continue_clustering = (ref_hit.energy() > min_energy);
} else { } else {
continue_clustering = false; continue_clustering = false;
break; break;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment