material.cc

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) BDSIM Collaboration, 2001 - 2026.
 
This file is part of BDSIM.
 
BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.
 
BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "material.h"
 
#include "array.h"
 
using namespace GMAD;
 
Material::Material()
{
  clear();
  PublishMembers();
}
 
void Material::clear()
{
  A = 0;
  Z = 0;
  density = 0;      //g*cm-3
  temper = 300;     //kelvin
  pressure = 0;     //atm
  state = "solid";  //allowed values: "solid", "liquid", "gas"
 
  components.clear();
  componentsFractions.clear();
  componentsWeights.clear();
}
 
void Material::PublishMembers()
{
  publish("name",   &Material::name);
  publish("A",      &Material::A);
  publish("Z",      &Material::Z);
  publish("density",&Material::density);
  publish("T"      ,&Material::temper);
  publish("P"      ,&Material::pressure);
  publish("state",  &Material::state);
  publish("components"         ,&Material::components);
  publish("componentsWeights"  ,&Material::componentsWeights);
  publish("componentsFractions",&Material::componentsFractions);
}
 
void Material::print()const
{
  std::cout << "material: " << name << " ";
 
  if (components.empty()) {
    std::cout << "A= " << A << "g/mole "
              << "Z= " << Z << " ";
  } else {
    std::cout << "ncomponents= " << components.size() << " ";
  }
  std::cout << "density= " << density  << "g/cm3 "
            << "state= "   << state    << " "
            << "T= "       << temper   << "K "
            << "P= "       << pressure << "atm "
            << std::endl;
}
 
// TODO template specialisation for Array pointers, to be merged into templated function
void Material::set_value(const std::string& property, Array* value, bool bExit)
{
#ifdef BDSDEBUG
  std::cout << "parser> Setting value " << std::setw(25) << std::left << property << std::endl;
#endif
  if(property=="components")
    {
      value->set_vector(components);
    } 
  else if(property=="componentsWeights")
    {
      value->set_vector(componentsWeights);
    }
  else if(property=="componentsFractions")
    {
      value->set_vector(componentsFractions);
    }
  else
    {
      std::cerr << "Error: parser> unknown material option \"" << property << "\", or doesn't expect vector type" << std::endl;
      if(bExit)
        {exit(1);}
      else
        {std::rethrow_exception(std::current_exception());} // to be caught by python
    }
}
 
#if __cplusplus >= 201703L
std::list<std::variant<bool, int, double, std::string>> Material::get_value_array(const std::string & property) {
  std::list<std::variant<bool, int, double, std::string>> retval(components.size());
 
  if(property=="components")
  {
    std::copy(components.begin(),components.end(), retval.begin());
  }
  else if(property=="componentsWeights")
  {
    std::copy(componentsWeights.begin(),componentsWeights.end(), retval.begin());
  }
  else if(property=="componentsFractions")
  {
    std::copy(componentsFractions.begin(),componentsFractions.end(), retval.begin());
  }
  else
  {
    std::cerr << "Error: parser> unknown material option \"" << property << "\", or doesn't expect vector type" << std::endl;
  }
 
  return retval;
}
#endif