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"
 
// HepMC3 includes
#include "HepMC3/GenParticle.h"
 
namespace marley {
 
  // Forward-declare the StructureDatabase class
  class StructureDatabase;

  class HauserFeshbachDecay {
 
    public:

      HauserFeshbachDecay( const std::shared_ptr< HepMC3::GenParticle >&
        compound_nucleus, double Exi, int twoJi, marley::Parity Pi,
        marley::StructureDatabase& sdb );

      const marley::ExitChannel& do_decay( double& Exf, int& twoJf,
        marley::Parity& Pf, std::shared_ptr< HepMC3::GenParticle >&
        emitted_particle, std::shared_ptr< HepMC3::GenParticle >&
        residual_nucleus, int& qIon, 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;
 
      inline double total_width() const;
 
    private:

      void build_exit_channels( marley::StructureDatabase& sdb );

      const std::shared_ptr< HepMC3::GenParticle > 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 double HauserFeshbachDecay::total_width() const
    { return total_width_; }
}

inline std::ostream& operator<<( std::ostream& out,
  const marley::HauserFeshbachDecay& hfd )
{
  hfd.print(out);
  return out;
}