/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/RecoVertex/VertexTools/interface/lms_3d.h File Reference

#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include <vector>

Go to the source code of this file.

Functions
GlobalPoint	lms_3d (std::vector< GlobalPoint > values)

---

Function Documentation

GlobalPoint lms_3d

(

std::vector< GlobalPoint >

values

)

least median of squares in three dimensions, doing every dimension separately

Definition at line 8 of file lms_3d.cc.

References i.

Referenced by Lms3d::operator()(), and LmsModeFinder3d::operator()().

{
  const int sze = values.size();
  if ( sze == 0 ) {
      throw VertexException("lms_3d: no values given.");
  };
  std::vector <float> x_vals, y_vals, z_vals;
  x_vals.reserve(sze-1);
  y_vals.reserve(sze-1);
  z_vals.reserve(sze-1);
  for (std:: vector<GlobalPoint>::iterator i=values.begin();
      i!=values.end() ; i++ ) {
    x_vals.push_back( i->x() );
    y_vals.push_back( i->y() );
    z_vals.push_back( i->z() );
  };
  return GlobalPoint ( lms_1d(x_vals), lms_1d(y_vals), lms_1d(z_vals) );
}