draw_mass_excess.C

#include "TStyle.h"
#include "TROOT.h"
#include "tkmanager.h"
#include "tknuclear_chart.h"
 
using namespace  tkn;
 
void draw_mass_excess(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("Mass Excess in MeV",tknuclear_chart::kAll,false);
    for(const auto &nuc : dtm.get_nuclei()) {
        auto mass_excess = nuc->get_mass_excess_measure();
        if(!mass_excess) continue;
        double yy = mass_excess->get_value(tkn::tkunit_manager::units_keys::MeV);
        nn->set_value(nuc->get_z(), nuc->get_n(), yy);
    }
 
    nn->draw();
    nn->hide_empty_bins();
 
    gROOT->SetBatch(true);
    gROOT->ForceStyle(true);
    gStyle->SetImageScaling(2.);
    nn->save_as("draw_mass_excess.png");
    gROOT->SetBatch(false);
}
 
void draw_mass_excess_no_systematic(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("Measured Mass Excess in MeV",tknuclear_chart::kAll,false);
    for(const auto &nuc : dtm.get_nuclei()) {
        auto mass_excess = nuc->get_mass_excess_measure();
        if(!mass_excess) continue;
        double yy = mass_excess->get_value(tkn::tkunit_manager::units_keys::MeV);
        if(!(mass_excess->get_info()==tkproperty::kSystematic)) nn->set_value(nuc->get_z(), nuc->get_n(), yy);
    }
 
    nn->draw();
    nn->hide_empty_bins();
 
    gROOT->SetBatch(true);
    gROOT->ForceStyle(true);
    gStyle->SetImageScaling(2.);
    nn->save_as("draw_mass_excess_no_systematic.png");
    gROOT->SetBatch(false);
}
 
void draw_mass_excess_from_systematic(){
 
    tkn::tkmanager dtm;
    auto* nn = new tknuclear_chart("Mass Excess in MeV from systematics",tknuclear_chart::kAll,false);
    for(const auto &nuc : dtm.get_nuclei()) {
        auto mass_excess = nuc->get_mass_excess_measure();
        if(!mass_excess) continue;
        double yy = mass_excess->get_value(tkn::tkunit_manager::units_keys::MeV);
        if((mass_excess->get_info()==tkproperty::kSystematic)) nn->set_value(nuc->get_z(), nuc->get_n(), yy);
    }
 
    nn->draw();
    nn->hide_empty_bins();
 
    gROOT->SetBatch(true);
    gROOT->ForceStyle(true);
    gStyle->SetImageScaling(2.);
    nn->save_as("draw_mass_excess_from_systematic.png");
    gROOT->SetBatch(false);
}