StatisticsFinder.h

Go to the documentation of this file.

#ifndef STATISTICS_FINDER_HH
#define STATISTICS_FINDER_HH

#include "Mint/MessageService.h"

class StatisticsFinder {

  double _min;     
  double _max;     
  double _nEvents; 
  double _wSum;    
  double _wSum2;   
  double _wSum3;   
  double _wSum4;   
  double _sumW;    
  //Can specify to keep track of the ordered events
  //allowing things like the median to be calculated

  int _keepOrderedEvents; 
  mutable std::vector<double> _orderedEvents; 
  bool needOrderedEvents() const;
  void warnIfWeightedEvents() const;
  bool notEnoughInformation(const double& val) const;

  public:

  StatisticsFinder(bool mean = 1, bool width = 1, bool widthError = 1, bool keepOrderedEvents = 0);

  double median() const;
  
  double numEvents() const{return _nEvents;}  
  void add(const double& x, const double& weight = 1.0);

  double mean() const;
  double meanError() const;
  double varience() const;
  double width() const;
  double widthError() const;
 
  double expX() const; 
  double expX2() const;
  double expX3() const;
  double expX4() const;
  
  double secondCentralMom() const;
  double fourthCentralMom() const;

  const double& getMin() const{return _min;}  
  const double& getMax() const{return _max;}  
  double range() const{return _max - _min;} 
  virtual ~StatisticsFinder();

};

class MinMaxFinder : public StatisticsFinder{
  public:
  MinMaxFinder() : StatisticsFinder(0,0,0,0){} 
  ~MinMaxFinder(){} 
};

class MeanFinder : public StatisticsFinder{
  public:
  MeanFinder() : StatisticsFinder(1,0,0,0){} 
  ~MeanFinder(){} 
};

class WidthFinder : public StatisticsFinder{
  public:
  WidthFinder() : StatisticsFinder(1,1,0,0){} 
  ~WidthFinder(){} 
};

class WidthErrorFinder : public StatisticsFinder{
  public:
  WidthErrorFinder() : StatisticsFinder(1,1,1,0){} 
  ~WidthErrorFinder(){} 
};

class MedianFinder : public StatisticsFinder{
  public:
  MedianFinder() : StatisticsFinder(1,1,1,1){} 
  ~MedianFinder(){} 
};

#endif