TkN 2.1
Toolkit for Nuclei
draw_r42.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_r42(){
9
11 auto* nn = new tknuclear_chart("E(4^{+})/E(2^{+})",tknuclear_chart::kEven,0,true);
12 for(const auto &nuc : dtm.get_nuclei([](auto nuc) {return (nuc->get_z()%2==0 && nuc->get_n()%2==0);})) {
13
14 auto lev_scheme = nuc->get_level_scheme();
15
16 auto lvl2 = lev_scheme->get_level("2+1",true);
17 auto lvl4 = lev_scheme->get_level("4+1",true);
18 if(!lvl2||!lvl4) continue;
19
20 if(!lvl2->has_property("energy")||!lvl4->has_property("energy")) continue;
21 double e2 = lvl2->get_energy();
22 double e4 = lvl4->get_energy();
23 double r42 = e4/e2;
24
25 nn->set_value(nuc->get_z(), nuc->get_n(), r42);
26 }
27
28 nn->draw();
29 nn->save_as("draw_r42.png");
30}
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