Weight.h

Go to the documentation of this file.

//-*-c++-*-
//-*-Weight.h-*-
//   Written by James Monk and Andrew Pilkington
#ifndef WEIGHT_HH
#define WEIGHT_HH

#include <iostream>
#include <map>
#include <vector>

namespace Exhume {

  class Weight {
  public:
    Weight() { NPoints = 1000; };
    virtual ~Weight(){};
    inline std::map<double, double> GetFuncMap() { return (FuncMap); };
    inline double GetTotalIntegral() { return (TotalIntegral); };

    inline std::map<double, double> GetLineShape() { return (LineShape); };

  protected:
    virtual double WeightFunc(const double &) = 0;

    void AddPoint(const double &, const double &);
    inline double GetFunc(const double &xx_) {
      if (xx_ > Max_) {
        return (WeightFunc(xx_));
      }

      std::map<double, double>::iterator high_, low_;
      high_ = FuncMap.upper_bound(xx_);
      low_ = high_;
      low_--;

      return (low_->second + (high_->second - low_->second) * (xx_ - low_->first) / (high_->first - low_->first));
    };

    inline double GetValue(const double &xx_) {
      std::map<double, double>::iterator high_, low_;
      high_ = LineShape.upper_bound(xx_);

      if (high_ == LineShape.end())
        high_--;

      low_ = high_;
      low_--;

      return (low_->second + (high_->second - low_->second) * (xx_ - low_->first) / (high_->first - low_->first));
    };
    void WeightInit(const double &, const double &);

    double Max_;
    double TotalIntegral;

  private:
    unsigned int NPoints;
    std::map<double, double> FuncMap;
    std::map<double, double> LineShape;
  };
}  // namespace Exhume

#endif