TkN 2.1
Toolkit for Nuclei
draw_FY.C
1#include "Rtypes.h"
2#include "Riostream.h"
3#include "tkmanager.h"
4#include "tknuclear_chart.h"
5
6using namespace tkn;
7
8void draw_FY239Pu(){
9
11 auto* nn = new tknuclear_chart("^{239}Pu fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
12 for(const auto &nuc : dtm.get_nuclei()) {
13 double yy = nuc->get_fission_yield("239Pu",false);
14 if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
15 }
16
17 nn->draw();
18 nn->save_as("draw_FY239Pu.png");
19}
20
21void draw_FY252Cf(){
22
24 auto* nn = new tknuclear_chart("^{252}Cf fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
25 for(const auto &nuc : dtm.get_nuclei()) {
26 double yy = nuc->get_fission_yield("252Cf",false);
27 if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
28 }
29 nn->draw();
30 nn->save_as("draw_FY252Cf.png");
31}
32
33void draw_cumulative_238U(){
34
36 auto* nn = new tknuclear_chart("^{238}U cumulative fission yield (in %)",tknuclear_chart::kAll,1e-6,true);
37 for(const auto &nuc : dtm.get_nuclei()) {
38 double yy = nuc->get_fission_yield("238U",true);
39 if(yy>0) nn->set_value(nuc->get_z(), nuc->get_n(), yy*50);
40 }
41 nn->draw();
42 nn->save_as("draw_cFY238U.png");
43}
Manages the database loading and provides access to the physics properties.
Definition: tkmanager.h:55
vector< shared_ptr< tknucleus > > get_nuclei(std::function< bool(shared_ptr< tknucleus >)>const &_selection)
return a vector containing all the known nuclei filtered by the lambda function
Definition: tkmanager.cpp:367
nuclear chart plot with ROOT
Definition: tklog.cpp:39