#include <vector>
#include <iostream>
#include <math.h>
//----------------------------------------------------------------------
// temporary storage of theta values
//----------------------------------------------------------------------
class Theta {
  private:
     std::vector<double>  m_array;  
  public:
     Theta(){}; 
     ~Theta(){};    // destructor
     void initTheta() { m_array.clear();}
     void addTheta(double theta){ m_array.push_back( theta ); }
     double getTheta();
};

int main()
{
	double x;
        Theta s;
        s.initTheta();
        for(int i = 0; i < 15 ; ++i ) {
          x=0.5982;
          if(i>10) x = x*2.;
	  s.addTheta(x);
	}
	std::cout<<" getTheta="<<s.getTheta()<< std::endl;
 	return 0;
} 

double Theta::getTheta() {
      double xmax, countmax;
      if ( m_array.size()< 5) return -999.;
      double count = 1.;
      double x = m_array[0];
      double xsum = x;
      xmax = x;
      countmax = 1.;
      for (unsigned int i=1; i<m_array.size(); i++) {
          if ( fabs (x-m_array[i]) < 0.2 ) {
             count += 1.;
             xsum += m_array[i];
             x = xsum / count;
          } else
          { if ( count > countmax ) {
                xmax = x;
                countmax = count;
             }
            count = 1.;
            x = m_array[i];
            xsum = x;
          }
          if (i == m_array.size()-1 && count > countmax) {
            xmax = x;
          }
      }
      if( countmax < 5)  return -999.;
      return xmax;
}