OpticalModel.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 <complex>
 
namespace marley {
 
  class OpticalModel {
 
    public:
 
      OpticalModel(int Z, int A) : Z_(Z), A_(A) {}
 
      virtual ~OpticalModel() = default;
 
      virtual std::complex<double> optical_model_potential(double r,
        double fragment_KE_lab, int fragment_pdg, int two_j, int l, int two_s,
        int target_charge = 0) = 0;
 
      virtual double transmission_coefficient(double total_KE_CM, int fragment_pdg,
        int two_j, int l, int two_s, int target_charge = 0) = 0;
 
      virtual double total_cross_section(double fragment_KE_lab,
        int fragment_pdg, int two_s, size_t l_max, int target_charge = 0) = 0;
 
      inline int Z() const;
 
      inline int A() const;
 
    protected:
 
      // Nuclear atomic and mass numbers
      int Z_, A_;
  };
 
  // Inline function definitions
  inline int OpticalModel::Z() const { return Z_; }
 
  inline int OpticalModel::A() const { return A_; }
}