tkdecay.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 tkdecay_H
#define tkdecay_H
 
#include "tkn_config.h"
 
#ifdef HAS_ROOT
#include "TClass.h"
#endif
 
#include <array>
 
#include "tkdb_table.h"
#include "tkproperty_list.h"
#include "tkmeasure.h"
 
namespace tkn {
 
class tklevel;
 
class tkdecay: public tkproperty_list {
    friend class tkdataset;
    friend class tklevel_scheme;
 
public:
 
    enum decay_type{ kundef, kbeta, kec, kalpha, kparticle, kgamma};
    static array<tkstring,6> decay_type_str;
 
protected:
 
    int fdecay_id;
    int fdecay_level_from_id;
    int fdecay_level_to_id;
 
    decay_type fdecay_type = kundef;
 
    tkstring funcertain_decay{};
 
    tkstring fcomment;
    bool fhas_comment=false;
 
    shared_ptr<tklevel> flevelfrom = nullptr;
    shared_ptr<tklevel> flevelto = nullptr;
 
public:
    tkdecay(int _id, int _from_id, int _to_id, const tkdb_table::measure_data_struct &_struc, decay_type _type);
    virtual ~tkdecay() = default;
 
    shared_ptr<tkmeasure> get_energy_meas() const {return get("energy");}
    double get_energy(const tkunit_manager::units_keys _unit=tkunit_manager::units_keys::keV);
 
    decay_type get_decay_type() {return fdecay_type;}
 
    bool has_comment() const {return fhas_comment;}
    const tkstring &get_comment() const {return fcomment;}
    bool is_uncertain() const {return !funcertain_decay.is_empty();}
 
    shared_ptr<tklevel> get_level_from() {return flevelfrom;}
    shared_ptr<tklevel> get_level_to() {return flevelto;}
 
    void print(const tkstring &_option="") const;
 
protected:
 
    int get_id() const {return fdecay_id;}
    int get_level_from_id() {return fdecay_level_from_id;}
    int get_level_to_id() {return fdecay_level_to_id;}
 
    void set_levels(shared_ptr<tklevel> _lvlfrom, shared_ptr<tklevel> _lvlto) {flevelfrom=_lvlfrom; flevelto=_lvlto;}
 
    void set_comment(const tkstring &_comment) {fcomment = _comment; fhas_comment=true;}
    void set_uncertain_decay(const tkstring &_uncertain) {funcertain_decay = _uncertain;}
 
#ifdef HAS_ROOT
    ClassDef(tkdecay,0);
#endif
};
 
class tkgammadecay: public tkdecay {
 
    friend class tkdataset;
    friend class tklevel_scheme;
 
public:
    tkgammadecay(int _id, int _from_id, int _to_id, const tkdb_table::measure_data_struct &_struc) : tkdecay(_id,_from_id,_to_id,_struc, kgamma) {;}
    virtual ~tkgammadecay() = default;
 
    shared_ptr<tkmeasure> get_relative_intensity_measure() {return get("RI");}
    double get_relative_intensity();
 
    shared_ptr<tkmeasure> get_mixing_ratio_measure() {return get("mix ratio");}
    double get_mixing_ratio();
 
    shared_ptr<tkmeasure> get_conv_coeff_measure() {return get("conv");}
    double get_conv_coeff();
 
    shared_ptr<tkmeasure> get_trans_prob_measure(bool _elec, int _L, bool _WU);
    double get_trans_prob(bool _elec, int _L, bool _WU);
 
    tkstring get_multipolarity() {return get_property("mult");}
 
    void set_relative_intensity(tkdb_table::measure_data_struct &_struc);
    void set_relative_intensity(double _val, double _unc);
 
    void set_mixing_ratio(tkdb_table::measure_data_struct &_struc);
 
    void set_multipolarity(const tkstring &_mul);
 
    void set_conv_coeff(tkdb_table::measure_data_struct &_struc);
 
    void set_transition_probability(tkdb_table::measure_data_struct &_struc);
 
#ifdef HAS_ROOT
    ClassDef(tkgammadecay,0);
#endif
};
 
}
 
#endif