SubsetHsmModeFinder3d Class Reference

#include <SubsetHsmModeFinder3d.h>

Inheritance diagram for SubsetHsmModeFinder3d:

Public Member Functions
SubsetHsmModeFinder3d *	clone () const override
GlobalPoint	operator() (const std::vector< PointAndDistance > &values) const override
Public Member Functions inherited from ModeFinder3d
virtual	~ModeFinder3d ()

Additional Inherited Members
Public Types inherited from ModeFinder3d
typedef std::pair< GlobalPoint, float >	PointAndDistance

Detailed Description

Definition at line 14 of file SubsetHsmModeFinder3d.h.

Member Function Documentation

SubsetHsmModeFinder3d * SubsetHsmModeFinder3d::clone ( void  ) const

overridevirtual

Implements ModeFinder3d.

Definition at line 52 of file SubsetHsmModeFinder3d.cc.

{ return new SubsetHsmModeFinder3d(*this); }

GlobalPoint SubsetHsmModeFinder3d::operator() ( const std::vector< PointAndDistance > &  values ) const

overridevirtual

Implements ModeFinder3d.

Definition at line 16 of file SubsetHsmModeFinder3d.cc.

References SplitLinear::begin, dataset::end, hsm_3d(), mps_fire::i, RPCpg::pts, and runTheMatrix::ret.

                                                                                             {
  if (values.empty()) {
    throw VertexException("SubsetHsmModeFinder3d: no value given.");
  };

  std::vector<GlobalPoint> pts;
  pts.reserve(values.size() - 1);
  std::vector<PointAndDistance> sorted_values(values.size());
  partial_sort_copy(values.begin(), values.end(), sorted_values.begin(), sorted_values.end(), compareByDistance());

  std::vector<PointAndDistance>::iterator end = sorted_values.end();
  std::vector<PointAndDistance>::iterator begin = sorted_values.begin();

  float dmax = 0.004;  // 40 microns, as a first try.

  // we want at least 30 values
  unsigned int min_num = values.size() < 30 ? values.size() : 30;

  // we also want at least 50 % of all values
  if (values.size() > 2 * min_num)
    min_num = (int)values.size() / 2;

  while (pts.size() < min_num) {
    // we cut at a dmax
    std::vector<PointAndDistance>::iterator i;
    for (i = begin; i != end && (i->second < dmax); ++i) {
      pts.push_back(i->first);
    };
    dmax += 0.003;  // add 30 microns with every iteration
    begin = i;
  };

  GlobalPoint ret = hsm_3d(pts);
  return ret;
}