Gamma.hh

//
// This file is part of MARLEY (Model of Argon Reaction Low Energy Yields)
//
// MARLEY is free software: you can redistribute it and/or modify it under the
// terms of version 3 of the GNU General Public License as published by the
// Free Software Foundation.
//
// For the full text of the license please see COPYING or
// visit http://opensource.org/licenses/GPL-3.0
//
// Please respect the MCnet academic usage guidelines. See GUIDELINES
// or visit https://www.montecarlonet.org/GUIDELINES for details.
 
#pragma once
#include "marley/IteratorToMember.hh"
 
namespace marley {
 
  class Level;
 
  class Gamma {
 
    public:
 
      Gamma(double energy = 0., double rel_intensity = 0.,
        marley::Level* start_lev = nullptr, marley::Level* end_lev = nullptr);
 
      inline marley::Level* start_level() const;
 
      inline void set_start_level(marley::Level* start_lev);
 
      inline marley::Level* end_level() const;
 
      inline void set_end_level(marley::Level* end_lev);
 
      inline double energy() const;
 
      inline double relative_intensity() const;
 
      template<typename It> inline static marley::IteratorToMember<It, double>
        make_intensity_iterator(It it);
 
    protected:
 
      double energy_; 
      double relative_intensity_; 
 
      marley::Level* start_level_;
 
      marley::Level* end_level_;
  };
 
  // Inline function definitions
  inline marley::Level* Gamma::start_level() const { return start_level_; }
  inline void Gamma::set_start_level(marley::Level* start_lev)
    { start_level_ = start_lev; }
 
  inline marley::Level* Gamma::end_level() const { return end_level_; }
  inline void Gamma::set_end_level(marley::Level* end_lev)
    { end_level_ = end_lev; }
 
  inline double Gamma::energy() const { return energy_; }
  inline double Gamma::relative_intensity() const
    { return relative_intensity_; }
 
  template<typename It> inline marley::IteratorToMember<It, double>
    Gamma::make_intensity_iterator(It it)
  {
    return marley::IteratorToMember<It, double>(it,
      &marley::Gamma::relative_intensity_);
  }
}