tkensdf_gamma_rec.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 tkensdf_gamma_rec_H
#define tkensdf_gamma_rec_H
 
#include "tkensdf_record.h"
#include "tkstring.h"
 
namespace tkn{
 
class tkensdf_gamma_rec : public tkensdf_record {
 
protected:
 
    bool fIs_final_level_set = false;
    double fFinal_level_energy = -1.;
 
    tkdb_table::measure_data_struct frelative_intensity;
    tkstring fmultipolarity="";
    tkdb_table::measure_data_struct fmixing_ratio;
    tkdb_table::measure_data_struct fconv_coeff;
 
    // transition probabilities
    tkdb_table::measure_data_struct fBE;
    tkdb_table::measure_data_struct fBEW;
    tkdb_table::measure_data_struct fBM;
    tkdb_table::measure_data_struct fBMW;
 
public:
 
    tkensdf_gamma_rec();
    tkensdf_gamma_rec(const tkensdf_gamma_rec &_gamma_record)  ;
    ~tkensdf_gamma_rec() {}

    virtual bool set_record(const tkstring &_record) override;
    virtual void analyse_record() override;

    void analyse_continuation_record() override;

    bool is_final_level_set() { return fIs_final_level_set;}

    void set_final_level_set(bool _status) { fIs_final_level_set = _status;}

    double get_final_level_energy() { return fFinal_level_energy;}
 
    void clear() override;
    void print(std::ostream &) const override;
 
private:
    void set_relative_intensity(const tkstring &_record);
    void set_multipolarity(const tkstring &_record);
    void set_mixing_ratio(const tkstring &_record);
    void set_conv_coeff(const tkstring &_record);
    void set_red_trans_prob(const tkstring &_record);
    void check_red_trans_prob();
    double weisskopf_conversion(bool electric, double mass, double mult);
 
public:
    const tkdb_table::measure_data_struct  &get_relative_intensity() const {return frelative_intensity;}
    const tkstring &get_multipolarity() const {return fmultipolarity;}
    const tkdb_table::measure_data_struct  &get_mixing_ratio() const {return fmixing_ratio;}
    const tkdb_table::measure_data_struct  &get_conv_coeff() const {return fconv_coeff;}
    const tkdb_table::measure_data_struct  &get_BE() const {return fBE;}
    const tkdb_table::measure_data_struct  &get_BEW() const {return fBEW;}
    const tkdb_table::measure_data_struct  &get_BM() const {return fBM;}
    const tkdb_table::measure_data_struct  &get_BMW() const {return fBMW;}
 
};
}
 
#endif