draw_FY.C

Explore Fission Yields

This example macro fills the nuclear chart with fission yields.

Start tkn-root, compile and execute this macro :

tkn-root
 _____    _   _   |  Documentation: https://tkn.in2p3.fr/
(_   _)  | \ | |  |  Source: https://gitlab.in2p3.fr/tkn/tkn-lib
  | |_  _|  \| |  |
  | | |/ /     |  |  Version 2.1
  | |   <| |\  |  |
  |_|_|\_\_| \_|  |  Database: TkN_ensdf_250201_xundl_240130_v2.1.db
 
tkn [0]
tkn [0] .L draw_FY.C+
tkn [1] draw_FY252Cf()
tkn [2] draw_FY239Pu()
tkn [3] draw_cumulative_238U()

It produces the following plot for ²⁵²Cf:

and the following plot for ²³⁹Pu:

and the following plot for ²³⁸U:

Source code :

#include "Rtypes.h"
#include "Riostream.h"
#include "tkmanager.h"
#include "tknuclear_chart.h"
 
using namespace  tkn;
 
void draw_FY239Pu(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("^{239}Pu fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
    for(const auto &nuc : dtm.get_nuclei()) {
        double yy = nuc->get_fission_yield("239Pu",false);
        if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
    }
 
    nn->draw();
    nn->save_as("draw_FY239Pu.png");
}
 
void draw_FY252Cf(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("^{252}Cf fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
    for(const auto &nuc : dtm.get_nuclei()) {
        double yy = nuc->get_fission_yield("252Cf",false);
        if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
    }
    nn->draw();
    nn->save_as("draw_FY252Cf.png");
}
 
void draw_cumulative_238U(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("^{238}U cumulative fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
    for(const auto &nuc : dtm.get_nuclei()) {
        double yy = nuc->get_fission_yield("238U",true);
        if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
    }
    nn->draw();
    nn->save_as("draw_cFY238U.png");
}