list_BE2_WU.C

#include "TH2F.h"
#include "TStyle.h"
 
#include "tkmanager.h"
#include "tknuclear_chart.h"
 
using namespace  tkn;
 
void list_BE2_WU()
{
    // to remove warning of emty datasets
    glog.set_warnings(false);
 
    // creation of a 2D histogram to plot the BE2 values on a nuclear chart view
    tknuclear_chart* nn = nullptr;
    nn = new tknuclear_chart("B(E2) (in Weisskopf unit)",tknuclear_chart::kEven,0,true);
 
    std::cout << "*-----------------------------------------*" << std::endl;
    std::cout << "* List of known B(E2) in Weisskopf unit   *" << std::endl;
    std::cout << "*-----------------------------------------*" << std::endl;
 
    cout << left << setw(7) << "Nuc" << "Data    " << setw(15) << "E gamma (keV)" << setw(6) << "mult" << setw(28) << "E from -> E to (keV)" << "Transition probability" << endl;
 
    // loop on all the even-even nuclei using a lambda function as parameter of the get_nuclei method
    for(const auto &nuc : gmanager->get_nuclei([](auto nuc) {return (nuc->get_z()%2==0 && nuc->get_n()%2==0);})) {
        // loop on all the datasets available on the ENSDF database
        for(const auto &dataset: nuc->get_level_scheme()->get_datasets()) {
 
            // Here we ask only data comming from the aggregated adopted datasets from ensdf, and the non evaluated XUNDL datasets
            bool xundl = dataset.second->get_name().contains("XUNDL");
            bool adopted = dataset.second->get_name().contains("ADOPTED");
            if(!(adopted || xundl)) continue;
 
            // we extract the pointer to the level scheme
            auto lev_scheme = nuc->get_level_scheme();
            // and we select the selected dataset
            lev_scheme->select_dataset(dataset.first);
 
            // here we ask for (if existing) the gamma transition from the first 2+ state to the ground state
            auto gamma = lev_scheme->get_decay<tkgammadecay>("2+1->0+1",false);
            if(!gamma) continue;
 
            // here we ask, if known, the BE2 value (1: electric, 2:L, 1: Weisskopf units)
            auto BE2W = gamma->get_trans_prob(true,2,true);
            if(BE2W<=0||std::isnan(BE2W)) continue;
 
            // the following lines are for printing in the terminal some informations on the transition
            cout <<  setw(6) << nuc->get_symbol();
            if(xundl) cout << " [XUNDL] ";
            else cout << " [ENSDF] ";
            cout << setw(15) << setprecision(2) << gamma->get_energy();
            cout << setw(6) <<gamma->get_multipolarity();
            tkstring transname = tkstring::Form("%-4s (%-6.1f) --> %s (%2.1f)",
                                                gamma->get_level_from()->get_spin_parity()->get_jpi_str().data(),
                                                gamma->get_level_from()->get_energy(),
                                                gamma->get_level_to()->get_spin_parity_str().data(),
                                                gamma->get_level_to()->get_energy());
            cout << setw(28) <<  transname  << " :" << BE2W;
            cout<<endl;
 
            // and finally we fill the histogram with the given BE2 value in Weisskopf unit
            if(adopted) nn->set_value(nuc->get_z(),nuc->get_n(),BE2W);
        }
    }
    // we plot the 2D histogram
    nn->draw();
    nn->save_as("list_BE2_WU.png");
}