MtvClusterizer1D< T > Class Template Reference

#include <MtvClusterizer1D.h>

Inheritance diagram for MtvClusterizer1D< T >:

Public Member Functions
virtual MtvClusterizer1D *	clone () const
	MtvClusterizer1D (const MtvClusterizer1D &)
	MtvClusterizer1D (const WeightEstimator< T > &est=TrivialWeightEstimator< T >())
std::pair< std::vector< Cluster1D< T > >, std::vector< const T * > >	operator() (const std::vector< Cluster1D< T > > &) const
	~MtvClusterizer1D ()
Public Member Functions inherited from Clusterizer1D< T >
virtual	~Clusterizer1D ()

Private Attributes
float	theErrorStretchFactor
WeightEstimator< T > *	theEstimator

Detailed Description

template<class T>
class MtvClusterizer1D< T >

Hsm clusterizer in one dimension, originally designed for ApexPoint Finding

Definition at line 17 of file MtvClusterizer1D.h.

Constructor & Destructor Documentation

MtvClusterizer1D() [1/2]

template<class T >

MtvClusterizer1D< T >::MtvClusterizer1D

(

const WeightEstimator< T > &

est = TrivialWeightEstimator<T>()

)

Definition at line 41 of file MtvClusterizer1D.h.

: theEstimator(est.clone()) {}

MtvClusterizer1D() [2/2]

template<class T>

MtvClusterizer1D< T >::MtvClusterizer1D

(

const MtvClusterizer1D< T > &

)

~MtvClusterizer1D()

template<class T >

MtvClusterizer1D< T >::~MtvClusterizer1D

(

)

Definition at line 44 of file MtvClusterizer1D.h.

                                       {
  delete theEstimator;
}

Member Function Documentation

clone()

template<class T >

MtvClusterizer1D< T > * MtvClusterizer1D< T >::clone ( void  ) const

virtual

Implements Clusterizer1D< T >.

Definition at line 49 of file MtvClusterizer1D.h.

                                                      {
  return new MtvClusterizer1D<T>(*this);
}

operator()()

template<class T >

std::pair< std::vector< Cluster1D< T > >, std::vector< const T * > > MtvClusterizer1D< T >::operator() ( const std::vector< Cluster1D< T > > &  ov ) const

virtual

Implements Clusterizer1D< T >.

Definition at line 54 of file MtvClusterizer1D.h.

                                              {
  typedef Cluster1D<T> Cluster1D;
  using namespace Clusterizer1DCommons;
  std::vector<const T*> unusedtracks;
  switch (ov.size()) {
    case 0:
      throw Clustering1DException("[MtvClusterizer1D] no values given");
    case 1:
      std::pair<std::vector<Cluster1D>, std::vector<const T*> > ret(ov, unusedtracks);
      return ret;
  };
  std::vector<Cluster1D> v = ov;
  sort(v.begin(), v.end(), ComparePairs<T>());
  std::vector<Cluster1D> sols;
  std::vector<const T*> trks;
 
  typename std::vector<Cluster1D>::iterator cur = v.begin();
  typename std::vector<Cluster1D>::iterator end = (v.end() - 1);
  double cur_min = cur->weight() + (cur + 1)->weight();
 
  for (typename std::vector<Cluster1D>::iterator i = v.begin(); i != end; ++i) {
    double cur_val = i->weight() + (i + 1)->weight();
    if (cur_val > cur_min) {
      cur_min = cur_val;
      cur = i;
    };
  };
 
  double weight = (cur->weight() + (cur + 1)->weight());
  double est = (cur->weight() * cur->position().value() + (cur + 1)->weight() * (cur + 1)->position().value()) / weight;
  double sigma = sqrt(square(cur->position().value() - est) + square((cur + 1)->position().value() - est));
  double err = 0.;
  int inliers = 0;
 
  for (typename std::vector<Cluster1D>::iterator i = v.begin(); i != v.end(); ++i) {
    if (fabs(i->position().value() - est) < 3 * sigma) {
      // all within 3 sigma are 'in'
      add(i->tracks(), trks);
      err += square(i->position().value() - est);
      inliers++;
    } else {
      add(i->tracks(), unusedtracks);
    };
  };
  err /= (inliers - 1);  // the algo definitely produces 2 or more inliers
  err = sqrt(err);
 
  sols.push_back(Cluster1D(Measurement1D(est, err), trks, weight));
  std::pair<std::vector<Cluster1D>, std::vector<const T*> > ret(sols, unusedtracks);
  return ret;
}

References Clusterizer1DCommons::add(), mps_fire::end, submitPVResolutionJobs::err, mps_fire::i, position, runTheMatrix::ret, jetUpdater_cfi::sort, mathSSE::sqrt(), Clusterizer1DCommons::square(), hltEgammaHLTExtra_cfi::trks, findQualityFiles::v, and mps_merge::weight.

Member Data Documentation

theErrorStretchFactor

template<class T>

float MtvClusterizer1D< T >::theErrorStretchFactor

private

Definition at line 29 of file MtvClusterizer1D.h.

theEstimator

template<class T>

WeightEstimator<T>* MtvClusterizer1D< T >::theEstimator

private

Definition at line 28 of file MtvClusterizer1D.h.