tknucleus.cpp

/********************************************************************************
 *   Copyright (c) : Université de Lyon 1, CNRS/IN2P3, UMR5822,                 *
 *                   IP2I, F-69622 Villeurbanne Cedex, France                   *
 *                   Normandie Université, ENSICAEN, UNICAEN, CNRS/IN2P3,       *
 *                   LPC Caen, F-14000 Caen, France                             *
 *   Contibutor(s) :                                                            *
 *      Jérémie Dudouet jeremie.dudouet@cnrs.fr [2020]                          *
 *      Diego Gruyer    diego.gruyer@cnrs.fr    [2020]                          *
 *                                                                              *
 *   Licensed under the MIT License <http://opensource.org/licenses/MIT>.       *
 *   SPDX-License-Identifier: MIT                                               *
 ********************************************************************************/
 
#include "tknucleus.h"
#include "tkstate_time.h"
#include "tkn_sqlite3.h"
#include "tklog.h"
#include "tkstring.h"
#include "tkmanager.h"
 
namespace tkn {
}
 
using namespace tkn;
 
vector<int> tknucleus::fmagic_number_list{2, 8, 20, 28, 50, 82, 126};

tknucleus::tknucleus(int _Z, int _A) : tkproperty_list("nucleus")
{
   fZ = _Z;
   fN = _A - _Z;
   fA = _A;
   tkstring symbol = gmanager->get_element_symbol(_Z);
   if (symbol.empty()) {
      glog << warning << "Unknown nucleus. No properties loaded " << tkstring::form("(%d,%d)", _Z, _A) << do_endl;
      return;
   }
   symbol.prepend(tkstring::form("%d", _A));
   if (!gmanager->known_nucleus(symbol)) {
      glog << warning << "Unknown nucleus. No properties loaded (" << symbol << ")" << do_endl;
      return;
   }
   const tknucleus &nuc = *gmanager->get_nucleus(symbol).get(); // vérifier que ca fasse une copie et pas une reference
   (*this) = nuc;
}

tknucleus::tknucleus(int _Z) : tkproperty_list("nucleus")
{
   if (!gmanager->known_element(_Z))
      glog << warning << "Unknown element. No properties loaded" << do_endl;
   else {
      shared_ptr<tknucleus> bestnuc;
      for (const auto &nuc : gmanager->get_nuclei_for_z(_Z)) {
         if (bestnuc == nullptr)
            bestnuc = nuc;
         else if (isnan(bestnuc->get_lifetime()) || (nuc->get_lifetime() > bestnuc->get_lifetime()))
            bestnuc = nuc;
         else if (nuc->get_lifetime() == bestnuc->get_lifetime() && nuc->get_abundance() > bestnuc->get_abundance())
            bestnuc = nuc;
      }
      const tknucleus &nuc = *bestnuc.get();
      (*this) = nuc;
   }
}

tknucleus::tknucleus(const char *_symbol) : tkproperty_list("nucleus")
{
   if (!gmanager->known_nucleus(tkstring(_symbol))) {
      int z;
      if (gmanager->known_element(tkstring(_symbol), z)) {
         if (!tkstring(_symbol).is_alpha()) {
            glog << warning << "Unknown nucleus. TkN will choose one for you ;)" << do_endl;
         }
         const tknucleus nuc(z);
         (*this) = nuc;
         return;
      }
      glog << warning << "Unknown nucleus. No properties loaded (" << _symbol << ")" << do_endl;
      return;
   }
   const tknucleus &nuc = *gmanager->get_nucleus(_symbol);
   (*this) = nuc;
}
 
void tknucleus::set_name(const char *_name)
{
   fSymbol = _name;
   felement_symbol = fSymbol.extract_alpha();
   fA = fSymbol.atoi();
}
 
bool tknucleus::is_z_magic()
{
   if (std::find(fmagic_number_list.begin(), fmagic_number_list.end(), fZ) != fmagic_number_list.end()) return true;
   return false;
}
 
bool tknucleus::is_n_magic()
{
   if (std::find(fmagic_number_list.begin(), fmagic_number_list.end(), fN) != fmagic_number_list.end()) return true;
   return false;
}
 
shared_ptr<tklevel_scheme> tknucleus::get_level_scheme()
{
   return gmanager->get_level_scheme(get_symbol(), get_z(), get_a());
}

double tknucleus::get_radius() const
{
   if (!has_property("radius")) return -1.;
   return get("radius")->get_value();
}

double tknucleus::get_magnetic_dipole()
{
   if (!has_property("magnetic_dipole")) return std::nan("1");
   return get("magnetic_dipole")->get_value();
}

double tknucleus::get_electric_quadrupole()
{
   if (!has_property("electric_quadrupole")) return std::nan("1");
   return get("electric_quadrupole")->get_value();
}

double tknucleus::get_abundance() const
{
   if (!has_property("abundance")) return -1.;
   return get("abundance")->get_value();
}

double tknucleus::get_mass_excess(const tkunit_manager::units_keys _unit)
{
   if (!has_property("mass_excess")) return std::nan("1");
   if (gunits->get_type(_unit) != tkunit_manager::units_type::kEnergy) {
      glog << error << gunits->get_name(_unit) << " is not of 'energy' type (returns nan(1))" << do_endl;
      return std::nan("1");
   }
   return get("mass_excess")->get_value(_unit);
}

double tknucleus::get_binding_energy_over_a(const tkunit_manager::units_keys _unit)
{
   if (!has_property("binding_energy_overA")) return std::nan("1");
   if (gunits->get_type(_unit) != tkunit_manager::units_type::kEnergy) {
      glog << error << gunits->get_name(_unit) << " is not of 'energy' type (returns nan(1))" << do_endl;
      return std::nan("1");
   }
   return get("binding_energy_overA")->get_value(_unit);
}

double tknucleus::get_half_life(tkunit_manager::units_keys unit)
{
   return detail::state_time_value(get("lifetime"), is_stable(), unit, false);
}
 
double tknucleus::get_mean_lifetime(tkunit_manager::units_keys unit)
{
   return get_half_life(unit) / std::log(2.);
}
 
double tknucleus::get_decay_width(tkunit_manager::units_keys unit)
{
   return detail::state_time_value(get("lifetime"), is_stable(), unit, true);
}
 
double tknucleus::get_lifetime(const tkunit_manager::units_keys _unit)
{
   if (!has_property("lifetime")) return std::nan("1");
   if (gunits->get_type(_unit) != tkunit_manager::units_type::kTime && gunits->get_type(_unit) != tkunit_manager::units_type::kEnergy) {
      glog << error << gunits->get_name(_unit) << " is not of 'time' or 'energy' type (returns nan(1))" << do_endl;
      return std::nan("1");
   }
   return get("lifetime")->get_value(_unit);
}

double tknucleus::get_fission_yield(tkstring _parent, bool _cumulative)
{
   if (!_parent.begins_with("FY")) _parent.prepend("FY");
   if (_cumulative) _parent.prepend("c");
   if (!has_property(_parent)) return 0.;
   return get(_parent)->get_value();
}

tkstring tknucleus::get_lifetime_str()
{
   if (!has_property("lifetime")) return "";
   tkstring lifetime_str;
 
   if (is_stable())
      lifetime_str = "STABLE";
   else {
      vector<tkunit_manager::units_keys> units{tkunit_manager::y, tkunit_manager::d, tkunit_manager::h, tkunit_manager::min,
                                               tkunit_manager::s, tkunit_manager::ms, tkunit_manager::us, tkunit_manager::ns,
                                               tkunit_manager::ps, tkunit_manager::fs, tkunit_manager::as,
                                               tkunit_manager::MeV, tkunit_manager::keV, tkunit_manager::eV};
 
      for (auto &unit : units) {
         double lifetime = get_lifetime(unit);
         if (lifetime > 1.) {
            lifetime_str = tkstring::Form("%.3g %s", lifetime, std::get<0>(tkunit_manager::funits_properties.at(unit)).data());
            break;
         }
      }
   }
 
   return lifetime_str;
}

shared_ptr<tklevel> tknucleus::get_ground_state()
{
   return get_level_scheme()->get_ground_state();
}

bool tknucleus::is_bound() const
{
   if (!has_property("binding_energy_overA")) return false;
   return get("binding_energy_overA")->get_value() > 0;
}

vector<pair<tkstring, double>> tknucleus::get_decay_modes()
{
   vector<pair<tkstring, double>> decay_list;
   auto decays = get_decay_mode_str();
   auto tokkens = decays.replace_all("[", "]").tokenize("]");
   for (auto &dec : tokkens) {
      auto tokk = dec.tokenize(";");
      decay_list.push_back({tokk.front(), tokk.back().atof()});
   }
   return decay_list;
}

void tknucleus::print() const
{
   glog << info << fSymbol << " (Z=" << fZ << ", N=" << fN << ") properties:" << do_endl;
   glog << info << "Ground state configuration: ";
   if (get_jpi().is_empty())
      glog << "Unknown" << do_endl;
   else
      glog << get_jpi() << do_endl;
   if (is_stable())
      glog << info << "Stable nucleus, abundance: " << tkstring::form("%3.2f", get_abundance()) << " %" << do_endl;
   else {
      glog << info << "Radioactive nucleus:" << do_endl;
      if (get_lifetime_measure()) glog << info << get_lifetime_measure() << do_endl;
      glog << info << "Decay: " << get_property("decay_modes") << do_endl;
   }
}
 
#ifdef HAS_ROOT
ClassImp(tknucleus);
#endif