sigma_sum.hpp

/* sigma_sum.hpp */
/* $Id$ */
#ifndef __sigma_sum_hpp
#define __sigma_sum_hpp

#include <vector>

template<typename T> 
T sigma_sum(const size_t size,
        const T* values, 
        const float mult = 3)
{
  if(0 == size)
    return 0.;
  
  T M = 0.;
  for(size_t i=0; i<size; ++i)
    M += values[i];
  M /= size;
  
  T avg = 0.;
  size_t count = 0;
  for(size_t i=0; i<size; ++i)
    {
      if(fabs(values[i] - M) < M*mult)
    {
      avg += values[i]; 
      count++;
    }
    }
  
  if(count > 1)
    avg /= count;
  
  return avg;
};

template<typename T> 
T sigma_sum(const std::vector<T> &values, 
        const float mult = 3)
{
  return sigma_sum<T>(values.size(), &(values.front()), mult);
}

#endif /* sigma_sum.hpp */