---

FunctionStats.cxx

Go to the documentation of this file.

#include "FastShowerUtils/FunctionStats.h"
#include "FastShowerUtils/ParticleParameters.h"
#include <iostream>
#include <iomanip>
#include <stream.h>
#include <cmath>

using std::ostream;
namespace FastShower{

  //
  FunctionStats::FunctionStats():
    m_lastIsValid(false), m_lastX(0.), m_n(0), m_sumX(0.), m_sumXsq(0.){
  }
  //
  FunctionStats::FunctionStats(const FunctionStats& rhs){
    m_lastIsValid  = rhs.m_lastIsValid; 
    m_lastX        = rhs.m_lastX; 
    m_n            = rhs.m_n; 
    m_sumX         = rhs.m_sumX; 
    m_sumXsq       = rhs.m_sumXsq;
  }
  //
  FunctionStats& FunctionStats::operator=(const FunctionStats& rhs){
    if (&rhs == this) return *this;
    m_lastIsValid  = rhs.m_lastIsValid; 
    m_lastX        = rhs.m_lastX; 
    m_n            = rhs.m_n; 
    m_sumX         = rhs.m_sumX; 
    m_sumXsq       = rhs.m_sumXsq;
    return *this;
  }
  //
  void FunctionStats::bump(double x){
    m_lastIsValid=true;
    m_lastX=x;
    ++m_n;
    m_sumX+=x;
    m_sumXsq+=(x*x);
  }
  //
  void FunctionStats::bump(double x, const ParticleParameters*){
    this->bump(x);
  }
  //
  int FunctionStats::nCalls() const {return m_n;}
  //
  double FunctionStats::mean() const {
    return (m_n>0)? m_sumX/(static_cast<double>(m_n)) :0.;
  }
  //
  double FunctionStats::variance() const {
    if(m_n==0) return 0.;

    double calls = static_cast<double>(m_n);
    double xSqAv = m_sumXsq/calls;
    double xAvSq = (m_sumX*m_sumX)/(calls*calls);
    return sqrt(abs(xSqAv-xAvSq));
    
  }
  //
  double FunctionStats::lastX() const {return m_lastX;}
  //
  void FunctionStats::report(std::ostream& ostr) const { 
    //    this->header(ostr);
    ostr<<*this<<endl;
  }
  void FunctionStats::header(std::ostream& ostr) const {
    ostr.width(40);
    ostr<<"Calls";    
    ostr.width(10);
    ostr<<"Last x";
    ostr.width(10);
    ostr<<"Mean";
    ostr.width(10);
    ostr<<"Variance";
    ostr<<endl;
}
  std::ostream& operator<<(std::ostream& ostr, const FunctionStats& f){
    ostr.width(40);
    ostr<<f.nCalls()<<" ";
    
    ostr.setf(ios::scientific, ios::floatfield);
    ostr.precision(3);

    if(f.lastIsValid()){
      ostr<<f.lastX()<<" ";
    }else{
      ostr<<"---------"<<" ";
    }
    ostr<<f.mean()<<" "
        <<f.variance()<<" "
      ;
    return ostr;
  }
  void FunctionStats::prime(){m_lastIsValid=false;}
  bool FunctionStats::lastIsValid() const {return m_lastIsValid;}
  std::ostream& operator<<(std::ostream& ostr, const FunctionStats* f){
    ostr << *f;
    return ostr;
  }
}//namesapce

---