tknucleus.h

/********************************************************************************
 *   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                                               *
 ********************************************************************************/
 
#ifndef tknucleus_H
#define tknucleus_H
 
#include "tkn_config.h"
 
#ifdef HAS_ROOT
#include "TClass.h"
#endif
 
#include "tkproperty_list.h"
#include "tklevel_scheme.h"
#include "tkmeasure.h"
 
namespace tkn {
 
class tkmanager;

class tknucleus : public tkproperty_list {
 
   friend class tkmanager;
 
protected:
   static vector<int> fmagic_number_list;
 
   int fZ{};
   int fN{};
   int fA{};
 
   tkstring felement_symbol{};
   tkstring felement_name{};
   tkstring fSymbol{};
 
   bool fis_known = false;
   bool fis_stable = false;
 
   tkspin_parity fspin_parity{};
 
public:
   tknucleus() : tkproperty_list("nucleus") {}
   tknucleus(int _Z, int _A);
   tknucleus(int _Z);
   tknucleus(const char *_symbol);
   virtual ~tknucleus() override = default;

   const tkstring &get_element_symbol() { return felement_symbol; }
   const tkstring &get_element_name() { return felement_name; }
   tkstring get_element_state() const { return get_property("state"); }
   bool is_gas() const { return get_element_state() == "Gas"; }
   bool is_liquid() const { return get_element_state() == "Liquid"; }
   bool is_solid() const { return get_element_state() == "Solid"; }

   int get_z() { return fZ; }
   int get_n() { return fN; }
   int get_a() { return fA; }

   const tkstring &get_symbol() { return fSymbol; }

   bool is_z_magic();
   bool is_n_magic();
   bool is_doubly_magic() { return is_z_magic() && is_n_magic(); }

   shared_ptr<tklevel_scheme> get_level_scheme();

   shared_ptr<tkmeasure> get_radius_measure() { return get("radius"); }
   double get_radius() const;

   shared_ptr<tkmeasure> get_magnetic_dipole_measure() { return get("magnetic_dipole"); }
   double get_magnetic_dipole();

   shared_ptr<tkmeasure> get_electric_quadrupole_measure() { return get("electric_quadrupole"); }
   double get_electric_quadrupole();

   shared_ptr<tkmeasure> get_abundance_measure() { return get("abundance"); }
   double get_abundance() const;

   shared_ptr<tkmeasure> get_mass_excess_measure() { return get("mass_excess"); }
   double get_mass_excess(const tkunit_manager::units_keys _unit = tkunit_manager::units_keys::keV);

   shared_ptr<tkmeasure> get_binding_energy_over_a_measure() { return get("binding_energy_overA"); }
   double get_binding_energy_over_a(const tkunit_manager::units_keys _unit = tkunit_manager::units_keys::MeV);

   shared_ptr<tkmeasure> get_lifetime_measure() const { return get("lifetime"); }
   double get_lifetime(const tkunit_manager::units_keys _unit = tkunit_manager::units_keys::s);
   tkstring get_lifetime_str();

   double get_fission_yield(tkstring _parent, bool _cumulative = false);

   shared_ptr<tklevel> get_ground_state();

   tkstring get_decay_mode_str() { return get_property("decay_modes"); }

   int get_discovery_year() { return get_property("isotope_year_discovered").atoi(); }

   vector<pair<tkstring, double>> get_decay_modes();

   bool is_known() const { return fis_known; }
   bool is_stable() const { return fis_stable; }
   bool is_bound() const;

   const tkspin_parity &get_spin_parity() const { return fspin_parity; }
   const tkstring &get_jpi() const { return fspin_parity.get_jpi_str(); }

   const vector<shared_ptr<tkmeasure>> get_xrays() { return get_data_properties("XRay*"); }

   void print() const;
 
private:
   void set_name(const char *_name);
 
#ifdef HAS_ROOT
   ClassDefOverride(tknucleus, 0);
#endif
};
} // namespace tkn
 
#endif