macro_explore.C

#include "Rtypes.h"
#include "tkmanager.h"
 
using namespace  tkn;
 
void macro_explore()
{
    std::cout << "*******************************************" << std::endl;
    std::cout << "* Excited states on Carbon isotopic chain *" << std::endl;
    std::cout << "*******************************************" << std::endl;
 
    std::cout << "* Ground state spin parity                *" << std::endl;
    std::cout << "*-----------------------------------------*" << std::endl;
    for(const auto &nuc : gmanager->get_nuclei([](auto nuc) {return (nuc->get_z()==6);})) {
 
        tkn::tkstring test = nuc->get_symbol();
        auto lvls =nuc->get_level_scheme()->get_levels([](auto lvl) {return (lvl->get_energy()==0);});
 
        std::cout << std::setw(5) << std::left << nuc->get_symbol();
        if(!lvls.empty()) std::cout << lvls.at(0)->get_spin_parity()->get_jpi_str();
        std::cout << std::endl;
    }
 
    std::cout << "*-----------------------------------------*" << std::endl;
    std::cout << "* Energy of second 0+ state               *" << std::endl;
    std::cout << "*-----------------------------------------*" << std::endl;
    for(const auto &nuc : gmanager->get_nuclei([](auto nuc) {return (nuc->get_z()==6&&!(nuc->get_a()%2));})) {
        auto lvls = nuc->get_level_scheme()->get_levels([](auto lvl) {return (lvl->get_spin_parity()->get_jpi_str()=="0+");});
 
        std::cout << std::setw(5) << std::left << nuc->get_symbol();
        if(lvls.size()>1) lvls[1]->get_energy_measure()->print(false);
        else std::cout << "unkown";
        std::cout << std::endl;
    }
 
    std::cout << "*-----------------------------------------*" << std::endl;
    std::cout << "* All positive parity states              *" << std::endl;
    std::cout << "*-----------------------------------------*" << std::endl;
    for(const auto &nuc : gmanager->get_nuclei([](auto nuc) {return (nuc->get_z()==6);})) {
        std::cout << std::setw(5) << std::left << nuc->get_symbol();
        for(const auto &lvl : nuc->get_level_scheme()->get_levels([](auto lvl) {return (lvl->get_spin_parity()->get_parity().get_value()==tkparity::kParityPlus);})) {
            std::cout << std::setw(5) << std::left << lvl->get_spin_parity()->get_jpi_str();
        }
        std::cout << std::endl;
    }
}