Bootstrapper.h

Go to the documentation of this file.

#ifndef _FdLightCollectionEfficiency_Bootstrapper_h_
#define _FdLightCollectionEfficiency_Bootstrapper_h_

#include <utl/TabulatedFunctionErrors.h>
#include <list>
#include <vector>


namespace FdLightCollectionEfficiencyKG {

  class Bootstrapper {

  public:
    Bootstrapper(const std::list<utl::TabulatedFunctionErrors>& data, const int verbosity = 0);
    ~Bootstrapper() { }

    const utl::TabulatedFunctionErrors& GetMean() { return fMean; }
    bool MaxRelUncertaintyBelowThreshold(const double relUncertaintyThreshold, const double absValueThreshold = 0.);

    void SetVerbosity(const int verbosity) { fVerbosity = verbosity; }
    int GetVerbosity() { return fVerbosity; }
    
  private:
    void Bootstrap();

    enum VarianceEstimator {
      eSampleVariance = 0,
      eSampleVarianceWeighted,
      eMedianAbsoluteDeviation
    };

    std::pair<double, double> MeanWithUncertainty(const std::vector<double>& samples,
                                                  const std::vector<double>& uncertainties,
                                                  const VarianceEstimator est = eSampleVariance);

    double Mean(const std::vector<double>& samples);
    double WeightedMean(const std::vector<double>& samples, const std::vector<double>& uncertainties);
    double Median(const std::vector<double>& sortedSamples);

    double SampleVariance(const std::vector<double>& samples, const double mean);
    double WeightedSampleVariance(const std::vector<double>& samples, const std::vector<double>& uncertainties, const double mean);

    double MedianAbsoluteDeviation(const std::vector<double>& samples, const double mean);

    const std::list<utl::TabulatedFunctionErrors>& fData;
    int fVerbosity;

    utl::TabulatedFunctionErrors fMean;
  };

 }

#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: