histogram.h

Go to the documentation of this file.

//Copyright 2009 Thomas A Caswell
//tcaswell@uchicago.edu
//http://jfi.uchicago.edu/~tcaswell
//
//This program is free software; you can redistribute it and/or modify
//it under the terms of the GNU General Public License as published by
//the Free Software Foundation; either version 3 of the License, or (at
//your option) any later version.
//
//This program is distributed in the hope that it will be useful, but
//WITHOUT ANY WARRANTY; without even the implied warranty of
//MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
//General Public License for more details.
//
//You should have received a copy of the GNU General Public License
//along with this program; if not, see <http://www.gnu.org/licenses>.
//
//Additional permission under GNU GPL version 3 section 7
//
//If you modify this Program, or any covered work, by linking or
//combining it with MATLAB (or a modified version of that library),
//containing parts covered by the terms of MATLAB User License, the
//licensors of this Program grant you additional permission to convey
//the resulting work.
#ifndef HIST
#define HIST

#include <vector>
#include <exception>
#include <iostream>

namespace gnuplot
{
class Gnuplot;
}

namespace utilities{
//forward declare
class Generic_wrapper;
class Md_store;

class Histogram{
public:
  template<class T>
  void add_data_point(T in);
  
  std::vector<int> get_bin_values() const;

  std::vector<float> get_bin_edges() const;

  void print() const;
  int get_over_count(){
    return over_count_;
  }

  int get_under_count(){
    return under_count_;
  }
  Histogram(int num_bins, float bottom , float top);

  void output_to_wrapper(Generic_wrapper * wrapper_out,
                         const std::string & g_name,
                         const std::string & count_name,
                         const std::string & edges_name,
                         const Md_store* g_md_store) const;
  
  void display()const;
  void display(gnuplot::Gnuplot & g)const;
  
  
  
  ~Histogram(){};

protected:
  std::vector<int> hist_array_;




  int under_count_;

  int over_count_;

  int number_bins_;

  float top_edge_;
  
  float bottom_edge_;

  float bin_width_;
private:


};
  
template<class T>
void Histogram::add_data_point(T in){
  using std::cout;
  using std::endl;
  float tmp_in = (float) in;
  if(tmp_in < bottom_edge_){
    ++under_count_;
    return;
  }
  if(tmp_in>= top_edge_){
    ++over_count_;
    return;
  }
  
  //  ++(*(hist_array_ptr_ + ((int)((tmp_in - bottom_edge_)/bin_width_))));
  //  cout<<in<<"\t"<<((int)((tmp_in - bottom_edge_)/bin_width_))<<endl;
  try
  {
    ++hist_array_.at(((unsigned)((tmp_in - bottom_edge_)/bin_width_)));
  }
  catch(std::exception & e)
  {
    cout<<e.what()<<endl;
    cout<<tmp_in<<"\t"
        <<bottom_edge_<<"\t"
        <<bin_width_<<"\t"<<endl;
    
  }
}




}

#endif