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]                          *
 *                                                                              *
 *    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.          *
 ********************************************************************************/
 
#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 = (*gdatabase)["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