lms_1d.h File Reference

#include "CommonTools/Statistics/interface/StatisticsException.h"
#include <vector>
#include <algorithm>
#include <cmath>

Go to the source code of this file.

Functions
template<class T >
T	lms_1d (std::vector< T > values)

Function Documentation

lms_1d()

template<class T >

T lms_1d

(

std::vector< T >

values

)

Least Median Sum of Squares in one dimension. T must be sortable, 'addable', and dividable by 2.

Definition at line 15 of file lms_1d.h.

References mps_fire::i, findQualityFiles::size, jetsAK4_CHS_cff::sort, and contentValuesCheck::values.

Referenced by lms_3d().

                              {
  switch (values.size()) {
    case 0:
      throw StatisticsException("Lms of empty vector undefined");
    case 1:
      return (T) * (values.begin());
    case 2:
      return (T)((*(values.begin()) + *(values.end() - 1)) / 2);
  };
  const int size = values.size();
  const int half = size / 2;
  const int n_steps = (size + 1) / 2;

  sort(values.begin(), values.end());

  T i_begin = *(values.begin());
  T i_end = *(values.begin() + half);
  T div = (T)fabs(i_begin - i_end);
  for (typename std::vector<T>::const_iterator i = values.begin(); i != (values.begin() + n_steps); ++i) {
    if (fabs((*i) - (*(i + half))) < div) {
      i_begin = (*i);
      i_end = *(i + half);
      div = (T)fabs(i_begin - i_end);
    }
  }

  return (T)((i_begin + i_end) / 2);
}