NuclearData.hpp

#pragma once
#include <cmath>
#include <optional>
#include <vector>
 
namespace Nuclear {
 
// skin-thickness. Always same?
constexpr double default_t = 2.30;
 
constexpr auto FourLn3 = 4.0 * 1.098612289;
// 4.0 * std::log(3.0); // std::log not constepr
constexpr auto Pi2 = M_PI * M_PI;
 
struct Isotope {
  int Z;
  int A;
  std::optional<double> r_rms{};
  std::optional<double> I_N{};
  std::optional<int> parity{};
  std::optional<double> mu{};
  std::optional<double> q{};
};
 
//==============================================================================
Isotope findIsotopeData(int z, int a);
std::vector<Isotope> findIsotopeList(int z);
double find_rrms(int z, int a);
double find_mu(int z, int a);
int find_parity(int z, int a);
double find_spin(int z, int a);
 
//==============================================================================
double c_hdr_formula_rrms_t(double rrms, double t = default_t);
double rrms_formula_c_t(double c, double t = default_t);
 
double approximate_r_rms(int a, int z);
 
double deformation_effective_t(double c, double t, double beta);
 
double approximate_t_skin(int a);
 
} // namespace Nuclear