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]                          *
 *                                                                              *
 *    This software is governed by the CeCILL-B license under French law and    *
 *    abiding by the  rules of distribution of free  software.  You can use,    *
 *    modify  and/ or  redistribute  the  software under  the  terms of  the    *
 *    CeCILL-B license as circulated by CEA, CNRS and INRIA at the following    *
 *    URL \"http://www.cecill.info\".                                           *
 *                                                                              *
 *    As a counterpart to the access  to the source code and rights to copy,    *
 *    modify  and redistribute granted  by the  license, users  are provided    *
 *    only with a limited warranty  and the software's author, the holder of    *
 *    the economic  rights, and the  successive licensors have  only limited    *
 *    liability.                                                                *
 *                                                                              *
 *    In this respect, the user's attention is drawn to the risks associated    *
 *    with loading,  using, modifying  and/or developing or  reproducing the    *
 *    software by the user in light of its specific status of free software,    *
 *    that  may mean that  it is  complicated to  manipulate, and  that also    *
 *    therefore  means that it  is reserved  for developers  and experienced    *
 *    professionals having in-depth  computer knowledge. Users are therefore    *
 *    encouraged  to load  and test  the software's  suitability  as regards    *
 *    their  requirements  in  conditions  enabling the  security  of  their    *
 *    systems  and/or data to  be ensured  and, more  generally, to  use and    *
 *    operate it in the same conditions as regards security.                    *
 *                                                                              *
 *    The fact that  you are presently reading this means  that you have had    *
 *    knowledge of the CeCILL-B license and that you accept its terms.          *
 ********************************************************************************/
 
#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() = 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
    ClassDef(tknucleus,0);
#endif
};
}
 
#endif