HauserFeshbachDecay.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 <memory>
#include <ostream>
 
#include "marley/ExitChannel.hh"
#include "marley/Parity.hh"
 
namespace marley {
 
  // Forward-declare the StructureDatabase class
  class StructureDatabase;
 
  class HauserFeshbachDecay {
 
    public:
 
      HauserFeshbachDecay( const marley::Particle& compound_nucleus,
        double Exi, int twoJi, marley::Parity Pi,
        marley::StructureDatabase& sdb );
 
      bool do_decay( double& Exf, int& twoJf, marley::Parity& Pf,
        marley::Particle& emitted_particle, marley::Particle& residual_nucleus,
        marley::Generator& gen );
 
      void print( std::ostream& out ) const;
 
      inline std::vector< std::unique_ptr<marley::ExitChannel> >&
        exit_channels();
 
      inline const std::vector< std::unique_ptr<marley::ExitChannel> >&
        exit_channels() const;
 
      const std::unique_ptr<marley::ExitChannel>& sample_exit_channel(
        marley::Generator& gen) const;
 
    private:
 
      void build_exit_channels( marley::StructureDatabase& sdb );
 
      const marley::Particle compound_nucleus_;
      const double Exi_; 
      const int twoJi_; 
      const marley::Parity Pi_; 
 
      double total_width_ = 0.;
 
      std::vector< std::unique_ptr<marley::ExitChannel> > exit_channels_;
  };
 
  // Inline function definitions
  inline std::vector<std::unique_ptr<marley::ExitChannel> >&
    HauserFeshbachDecay::exit_channels() { return exit_channels_; }
 
  inline const std::vector<std::unique_ptr<marley::ExitChannel> >&
    HauserFeshbachDecay::exit_channels() const { return exit_channels_; }
}
 
inline std::ostream& operator<<( std::ostream& out,
  const marley::HauserFeshbachDecay& hfd )
{
  hfd.print(out);
  return out;
}