draw_binding.C

#include "Rtypes.h"
#include "Riostream.h"
#include "tkmanager.h"
#include "tknuclear_chart.h"
 
using namespace  tkn;
 
void draw_binding(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("Binding energy (in MeV/nuc)",tknuclear_chart::kAll,true);
    for(const auto &nuc : dtm.get_nuclei()) {
        nn->set_value(nuc->get_z(), nuc->get_n(), nuc->get_binding_energy_over_a());
    }
 
    nn->draw();
    nn->hide_empty_bins();
    nn->save_as("draw_binding.png");
}