ContinuumNuclearReaction.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
 
// Standard library includes
#include <memory>
 
// MARLEY includes
#include "marley/NuclearReaction.hh"
#include "marley/TabulatedXSec.hh"
 
namespace marley {

  class ContinuumNuclearReaction : public NuclearReaction {
 
    public:

      enum class SubContinuumMode { IGNORE, MIRROR, ACCUMULATE };
 
      ContinuumNuclearReaction( Reaction::ProcessType pt, int pdg_a,
        int pdg_b, int pdg_c, int pdg_d, int q_d,
        const std::shared_ptr<TabulatedXSec>& txsec,
        const std::string& source_file );
 
      virtual std::shared_ptr< HepMC3::GenEvent > create_event(
        int particle_id_a, double KEa, marley::Generator& gen ) const override;
 
      virtual double total_xs( int pdg_a, double KEa ) const override;
 
      inline const TabulatedXSec& get_tabulated_xsec() const
        { return *xsec_; }

      static inline SubContinuumMode sub_continuum_mode() { return sc_mode_; }

      static inline void set_sub_continuum_mode( SubContinuumMode scm )
        { sc_mode_ = scm; }
 
      static SubContinuumMode sub_continuum_mode_from_string(
        const std::string& sc_mode_str );
 
      static std::string string_from_sub_continuum_mode(
        const SubContinuumMode sc_mode );
 
    protected:
 
      virtual void set_description() override;

      bool reassign_sub_continuum( double& w, const double ctl,
        const double KEa ) const;

      double get_Ec_from_Ex( const double Ex, const double cos_theta,
        const double KEa, double* jacobian = nullptr) const;

      std::shared_ptr< TabulatedXSec > xsec_;

      static SubContinuumMode sc_mode_;

      static const std::map< SubContinuumMode, std::string > sc_mode_string_map_;
  };
 
}