tss::SimpleClustering Class Reference

#include "SimpleClustering.h"

Public Member Functions
std::vector< tss::Cluster2D >	run (const std::vector< tss::Hit2D > &inp) const
std::vector< tss::Cluster2D >	run (const tss::Cluster2D &inp) const
bool	hitsTouching (const tss::Hit2D &h1, const tss::Hit2D &h2) const
bool	hitsTouching (const tss::Cluster2D &c1, const tss::Hit2D &h2) const
bool	hitsTouching (const tss::Cluster2D &c1, const tss::Cluster2D &c2) const

Private Member Functions
void	merge (std::vector< tss::Cluster2D > &clusters) const

Detailed Description

Definition at line 122 of file SimpleClustering.h.

Member Function Documentation

bool tss::SimpleClustering::hitsTouching	(	const tss::Hit2D &	h1,
		const tss::Hit2D &	h2
	)		const

Definition at line 180 of file SimpleClustering.cxx.

References tss::Hit2D::EndTick(), tss::Hit2D::PeakTime(), tss::Hit2D::StartTick(), and tss::Hit2D::Wire().

Referenced by tss::Segmentation2D::buildSegment().

{
  if ((h1.Wire() == h2.Wire()) && (h1.PeakTime() == h2.PeakTime())) return false;

  bool touches = false;
  if ((h1.Wire() >= h2.Wire() - 1) && (h1.Wire() <= h2.Wire() + 1)) {
    if (((h2.StartTick() <= h1.StartTick()) && (h1.StartTick() <= h2.EndTick() + 1)) ||
        ((h2.StartTick() <= h1.EndTick() + 1) && (h1.EndTick() <= h2.EndTick())) ||
        ((h2.StartTick() >= h1.StartTick()) && (h1.EndTick() >= h2.EndTick()))) {
      touches = true;
    }
  }
  return touches;
}

bool tss::SimpleClustering::hitsTouching	(	const tss::Cluster2D &	c1,
		const tss::Hit2D &	h2
	)		const

Definition at line 196 of file SimpleClustering.cxx.

References tss::Cluster2D::size().

{
  for (size_t i = 0; i < c1.size(); i++) {
    if (hitsTouching(c1[i], h2)) return true;
  }
  return false;
}

bool tss::SimpleClustering::hitsTouching	(	const tss::Cluster2D &	c1,
		const tss::Cluster2D &	c2
	)		const

Definition at line 205 of file SimpleClustering.cxx.

References tss::Cluster2D::size().

{
  for (unsigned int i = 0; i < c1.size(); i++) {
    if (hitsTouching(c2, c1[i])) return true;
  }
  return false;
}

void tss::SimpleClustering::merge ( std::vector< tss::Cluster2D > &  clusters ) const

private

Definition at line 214 of file SimpleClustering.cxx.

References tss::Cluster2D::hits(), and tss::Cluster2D::size().

{
  bool merged = true;
  while (merged) {
    merged = false;

    size_t i = 0;
    while (i < clusters.size() - 1) {
      size_t j = i + 1;
      while (j < clusters.size()) {
        if (hitsTouching(clusters[i], clusters[j])) {
          clusters[i].hits().reserve(clusters[i].size() + clusters[i].size());
          for (size_t h = 0; h < clusters[j].size(); ++h)
            clusters[i].hits().push_back(clusters[j].hits()[h]);
          clusters.erase(clusters.begin() + j);
          merged = true;
        }
        else
          ++j;
      }
      ++i;
    }
  }
}

std::vector< tss::Cluster2D > tss::SimpleClustering::run

(

const std::vector< tss::Hit2D > &

inp

)

const

Definition at line 240 of file SimpleClustering.cxx.

References tss::Cluster2D::merge(), and r.

Referenced by tss::TrackShowerHits::produce(), and tss::Segmentation2D::run().

{
  std::vector<tss::Cluster2D> result;
  for (size_t h = 0; h < inp.size(); ++h) {
    bool found = false;
    for (size_t r = 0; r < result.size(); ++r)
      if (hitsTouching(result[r], inp[h])) {
        result[r].hits().push_back(&(inp[h]));
        found = true;
        break;
      }
    if (!found) {
      result.push_back(tss::Cluster2D());
      result.back().hits().push_back(&(inp[h]));
    }
  }
  merge(result);

  return result;
}

std::vector< tss::Cluster2D > tss::SimpleClustering::run

(

const tss::Cluster2D &

inp

)

const

Definition at line 262 of file SimpleClustering.cxx.

References tss::Cluster2D::hits(), tss::Cluster2D::merge(), r, and tss::Cluster2D::size().

{
  std::vector<tss::Cluster2D> result;
  for (size_t h = 0; h < inp.size(); ++h) {
    bool found = false;
    for (size_t r = 0; r < result.size(); ++r)
      if (hitsTouching(result[r], inp[h])) {
        result[r].hits().push_back(inp.hits()[h]);
        found = true;
        break;
      }
    if (!found) {
      result.push_back(tss::Cluster2D());
      result.back().hits().push_back(inp.hits()[h]);
    }
  }
  merge(result);

  return result;
}

---

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

• larreco/v09_25_00/source/larreco/ClusterFinder/TrackShowerSplitter/Segmentation2D/SimpleClustering.h
• larreco/v09_25_00/source/larreco/ClusterFinder/TrackShowerSplitter/Segmentation2D/SimpleClustering.cxx