Integrator.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include <functional>
#include <vector>
 
#include "marley/marley_utils.hh"
 
namespace marley {

  class Integrator {
    public:

      Integrator(size_t num = N_DEFAULT_);

      double num_integrate(const std::function<double(double)>& f,
        double a, double b) const;
 
    private:

      size_t N_;

      std::vector<double> weights_;

      std::vector<double> offsets_;
 
      static constexpr size_t N_DEFAULT_ = 20; 
  };
 
}