ElementXsecInterpolation.h

#ifndef G4DarkBREM_ELEMENTXSECINTERPOLATION_H
#define G4DarkBREM_ELEMENTXSECINTERPOLATION_H
 
#include <memory>
 
#include "G4DarkBreM/PrototypeModel.h"
 
namespace g4db {
 
class ElementXsecInterpolation {
 public:
  ElementXsecInterpolation() = default;
 
  ElementXsecInterpolation(std::shared_ptr<PrototypeModel> model)
      : model_{model} {}
 
  G4double get(G4double energy, G4double A, G4double Z);
 
  void stream(std::ostream& o) const;
 
  friend std::ostream& operator<<(std::ostream& o,
                                  const ElementXsecInterpolation c) {
    c.stream(o);
    return o;
  }
 
 private:
  class SampleSet {
   public:
    SampleSet(const std::vector<double>& init_x,
              const std::vector<double>& init_y);
 
    const double& min() const;
 
    const double& max() const;
 
    void prepend(double x, double y);
 
    void append(double x, double y);
 
    double interpolate(double x) const;
 
    void dump(std::ostream& o, int z, double kinematic_min) const;
 
   private:
    std::vector<double> x_;
    std::vector<double> y_;
  };
 
 private:
  std::map<int, SampleSet> the_samples_;
 
  std::shared_ptr<PrototypeModel> model_;
 
};  // ElementXsecInterpolation
 
}  // namespace g4db
 
#endif