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
};
 
} // namespace tkn
 
#endif