draw_r42.C

Explore R_4/2

This example macro fills the nuclear chart with R_4/2 values.

Start tkn-root, compile and execute this macro :

tkn-root
 _____    _   _   |  Documentation: https://tkn.in2p3.fr/
(_   _)  | \ | |  |  Source: https://gitlab.in2p3.fr/tkn/tkn-lib
  | |_  _|  \| |  |
  | | |/ /     |  |  Version 1.0
  | |   <| |\  |  |
  |_|_|\_\_| \_|  |  Database: TkN_ensdf_221101_xundl_210701_v1.0.db
 
 
tkn [0] .L draw_r42.C+
tkn [1] draw_r42()

It produces the following plot:

Source code :

#include "Rtypes.h"
#include "Riostream.h"
#include "tkmanager.h"
#include "tknuclear_chart.h"
 
using namespace  tkn;
 
void draw_r42(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("E(4^{+})/E(2^{+})",tknuclear_chart::kEven,0,true);
    for(const auto &nuc : dtm.get_nuclei([](auto nuc) {return (nuc->get_z()%2==0 && nuc->get_n()%2==0);})) {
 
        auto lev_scheme = nuc->get_level_scheme();
 
        auto lvl2 = lev_scheme->get_level("2+1",true);
        auto lvl4 = lev_scheme->get_level("4+1",true);
        if(!lvl2||!lvl4) continue;
 
        if(!lvl2->has_property("energy")||!lvl4->has_property("energy")) continue;
        double e2 = lvl2->get_energy();
        double e4 = lvl4->get_energy();
        double r42 = e4/e2;
 
        nn->set_value(nuc->get_z(), nuc->get_n(), r42);
    }
 
    nn->draw();
    nn->save_as("draw_r42.png");
}