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;}

    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() {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
};
}
 
#endif