Changeset 215 for trunk/src/Histogram.cc

Timestamp:

Nov 7, 2004, 12:07:03 AM (18 years ago)

Author:

Jari Häkkinen

Message:

Continued work on new Histogram.

File:

• trunk/src/Histogram.cc (modified) (5 diffs)

trunk/src/Histogram.cc

@@ -25 +25 @@
 
 
-Histogram::Histogram(const double min, const double max, const size_t)
-  : xmax_(max), xmin_(min), sum_all_(), sum_histogram_()
+Histogram::Histogram(const double min, const double max, const size_t n)
+  : histogram_(std::vector<double>(n,0.0)),
+    xmax_(max), xmin_(min),
+    sum_all_(), sum_histogram_()
 {
-  // Peter, fix me
 }
-
-
-/*
-Histogram::Histogram(const u_int nof_bins, const double min, const double max)
-  : spacing_((max-min)/nof_bins), xmax_(max), xmin_(min), sum_all_(),
-    sum_histogram_()
-{
-  histogram_=vector<statistics::WeightedAverager>(nof_bins);
-}
-*/
 
 
@@ -51 +42 @@
 int Histogram::add(const double x, const double w)
 {
+  std::cout << "Histogram::add Jari, check me!" << std::endl;
+
   sum_all_.add(x,w);
   if (x<xmin_)
@@ -64 +57 @@
 
 
-void Histogram::normalize(double factor)
-{
-  for (u_int i=0; i<histogram_.size(); i++)
-    histogram_[i]/=factor;
-}
-
-
-
-/*
 void Histogram::normalize(double scale_factor)
 {
-  for (u_int i=0; i<histogram_.size(); i++)
-    histogram_[i]/=_sum_histogram.nof_entries()*spacing_*scale_factor;
+  for (size_t i=0; i<histogram_.size(); i++)
+    histogram_[i]/=scale_factor;
 }
-*/
 
 
@@ -89 +72 @@
   sum_histogram_.reset();
 }
-
-
-
-/*
-const statistics::WeightedAverager& Histogram::operator[](double d) const
-{
-  if ((d<xmin_) || (d>xmax_))
-    return statistics::WeightedAverager();// requested bin not within boundaries
-  return histogram_[static_cast<u_int>(floor((d-xmin_) /
-               spacing_))];
-}
-*/
 
 
@@ -118 +89 @@
 
 
-std::ostream& operator<<(std::ostream& s,const Histogram& b)
+std::ostream& operator<<(std::ostream& s,const Histogram& histogram)
 {
-  s << "# min           : " << b.xmin() << '\n';
-  s << "# max           : " << b.xmax() << '\n';
-  s << "# spacing       : " << b.spacing() << '\n';
-  s << "# number of bins: " << b.nof_bins() << '\n';
-
-//   s << "# histogram averager: " << b.averager_histogram() << '\n';
-  s << "#            entries: " << b.averager_histogram().sum_w() << '\n';
-//   s << "# all averager      : " << b.averager_all() << '\n';
-  s << "#            entries: " << b.averager_all().sum_w() << '\n';
+  s << "# histogram min : " << histogram.xmin() << '\n';
+  s << "# histogram max : " << histogram.xmax() << '\n';
+  s << "# number of bins: " << histogram.nof_bins() << '\n';
+  s << "# histogram averager: " << histogram.averager_histogram() << '\n';
+  s << "#       all averager: " << histogram.averager_all() << '\n';
 
   s << "# column 1: center of observation bin\n"
-    << "# column 2: frequency\n"
-    << "# column 3: relative frequency\n"
-    << "# column 4: average of observations in bin\n"
-    << "# column 5: stddev of frequency\n";
+    << "# column 2: frequency\n";
 
-  double observation=b.xmin()+b.spacing()/2;
-  for (u_int i=0; i<b.nof_bins(); i++) {
+  for (u_int i=0; i<histogram.nof_bins(); i++) {
     s.width(12);
-    s << observation;
+    s << histogram.xmin()+(i+0.5)*histogram.spacing();
     s.width(12);
-    s << b[i];
-    s.width(12);
-// Peter, fix me
-//     if (b.nof_entries_all())
-//       s << b[i].sum_w()/b.nof_entries_all();
-//     else
-//       s << 0;
-    s.width(12);
-//     s << b[i] << "\n";
-    observation+=b.spacing();
+    s << histogram[i] << '\n';
   }