tklevel_builder.cpp

/********************************************************************************
 *   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                                               *
 ********************************************************************************/
 
#include "tklevel_builder.h"
 
#include <unistd.h>
 
#include "tklog.h"
#include "tkensdf_level_rec.h"
 
using namespace tkn;
using namespace std;
 
namespace tkn {
}
 
tklevel_builder::~tklevel_builder() = default;
 
tklevel_builder::tklevel_builder(tkdatabase *_database, const char *_table_name) : tkdb_builder(_database, _table_name)
{
   fTable = &fDataBase->new_table(fTableName);
 
   fTable->add_column("level_id", "INT NOT NULL PRIMARY KEY");
   fTable->add_column("isotope_id", "INT NOT NULL");
   fTable->add_column("dataset_id", "INT NOT NULL");
 
   add_measure("energy");
 
   fTable->add_column("level_spin", "REAL");
   fTable->add_column("level_parity", "INT");
   fTable->add_column("level_spin_parity", "TEXT");
 
   add_measure("lifetime");
 
   fTable->add_column("level_comment", "TEXT");
 
   fTable->add_column("level_uncertain", "TEXT");
 
   fTable->add_constraint("FOREIGN KEY (isotope_id) REFERENCES ISOTOPE (isotope_id),"); // hold this link with a constraints
   fTable->add_constraint("FOREIGN KEY (dataset_id) REFERENCES DATASET (dataset_id)");  // hold this link with a constraints
 
   fTable->write_to_database();
 
   glog << info << "Creating '" << _table_name << "' table" << do_endl;
}
 
void tklevel_builder::fill_level(int _dset_idx, int _isotope_idx, const tkensdf_level_rec &_lev_record, bool _adopted)
{
   fLevelIdx++;
   (*fTable)["level_id"].set_value(fLevelIdx);      // level id
   (*fTable)["isotope_id"].set_value(_isotope_idx); // link to isotope -> to be set for each level
   (*fTable)["dataset_id"].set_value(_dset_idx);    // link to dataset -> to be set for each level
 
   fill_measure("energy", _lev_record.get_energy());
 
   if (!_lev_record.get_Jpi_str().is_empty()) {
      if (_lev_record.get_J()->is_known()) (*fTable)["level_spin"].set_value(_lev_record.get_j());
      if (_lev_record.get_Pi()->is_known()) (*fTable)["level_parity"].set_value(_lev_record.get_pi());
      (*fTable)["level_spin_parity"].set_value(_lev_record.get_Jpi_str());
   }
 
   fill_measure("lifetime", _lev_record.get_lifetime());
 
   if (!_lev_record.get_comment_record().is_empty()) (*fTable)["level_comment"].set_value(_lev_record.get_comment_record());
 
   if (_lev_record.is_uncertain()) (*fTable)["level_uncertain"].set_value(_lev_record.get_uncertain_record().data());
 
   //    // fill the isotope ground state lifetime and spin
   if (_lev_record.get_energy().value == 0. && _adopted && !_lev_record.get_energy().info.contains("OFF")) {
      tkdb_table &isotope = (*fDataBase)["ISOTOPE"];
 
      if (!_lev_record.get_Jpi_str().is_empty()) {
         isotope["spin"].set_value(_lev_record.get_j());
         isotope["parity"].set_value(_lev_record.get_pi());
         isotope["spin_parity"].set_value(_lev_record.get_Jpi_str());
      }
 
      if (_lev_record.get_lifetime().filled) {
         isotope["lifetime"].set_value(_lev_record.get_lifetime().value);
         isotope["lifetime_unit"].set_value(_lev_record.get_lifetime().unit);
         if (_lev_record.get_lifetime().value == 0.) isotope["lifetime_unc"].set_value(0.);
         if (_lev_record.get_lifetime().err > 0.) isotope["lifetime_unc"].set_value(_lev_record.get_lifetime().err);
         if (_lev_record.get_lifetime().err_low > 0.) isotope["lifetime_unc_low"].set_value(_lev_record.get_lifetime().err_low);
         if (_lev_record.get_lifetime().err_high > 0.) isotope["lifetime_unc_high"].set_value(_lev_record.get_lifetime().err_high);
         if (!_lev_record.get_lifetime().info.is_empty()) isotope["lifetime_info"].set_value(_lev_record.get_lifetime().info);
      }
 
      isotope.update_row(tkstring::form("isotope_id=%d", _isotope_idx));
   }
 
   fTable->push_row();
}
 
#ifdef HAS_ROOT
ClassImp(tklevel_builder);
#endif