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]                          *
 *                                                                              *
 *   Licensed under the MIT License <http://opensource.org/licenses/MIT>.       *   
 *   SPDX-License-Identifier: MIT                                               *
 ********************************************************************************/
 
#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() override = 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
    ClassDefOverride(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() override = 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
    ClassDefOverride(tkgammadecay,0);
#endif
};
 
}
 
#endif