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 << "unknown";
      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;
   }
}