draw_FY.C

#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 * 100);
   }
 
   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 * 100);
   }
   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 * 100);
   }
   nn->draw();
   nn->save_as("draw_cFY238U.png");
}