cctbx/eltbx/neutron.h

#ifndef CCTBX_ELTBX_NEUTRON_H
#define CCTBX_ELTBX_NEUTRON_H

#include <string>
#include <complex>

namespace cctbx { namespace eltbx { namespace neutron {

  namespace detail
  {
    struct raw_record_neutron_news_1992
    {
      const char* label;
      float       bound_coh_scatt_length_real;
      float       bound_coh_scatt_length_imag;
      float       abs_cross_sect; // For 2200 m/s neutrons
    };
  }


  class neutron_news_1992_table
  {
    public:
      neutron_news_1992_table() : record_(0) {}


      explicit
      neutron_news_1992_table(std::string const& label, bool exact=false);


      bool
      is_valid() const { return record_->label != 0; }

      const char*
      label() const
      {
        return record_->label;
      }


      std::complex<float>
      bound_coh_scatt_length() const
      {
        return std::complex<float>(record_->bound_coh_scatt_length_real,
                                   record_->bound_coh_scatt_length_imag);
      }


      float
      bound_coh_scatt_length_real() const
      {
        return record_->bound_coh_scatt_length_real;
      }


      float
      bound_coh_scatt_length_imag() const
      {
        return record_->bound_coh_scatt_length_imag;
      }


      float
      abs_cross_sect() const
      {
        return record_->abs_cross_sect;
      }

    private:
      const detail::raw_record_neutron_news_1992* record_;
      friend class neutron_news_1992_table_iterator;
  };

  class neutron_news_1992_table_iterator
  {
    public:
      neutron_news_1992_table_iterator();


      neutron_news_1992_table
      next();

    private:
      neutron_news_1992_table current_;
  };

}}} // cctbx::eltbx::neutron

#endif // CCTBX_ELTBX_NEUTRON_H

---