CoulombCorrector.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 <map>
#include <string>
 
namespace marley {

  class CoulombCorrector {
 
    public:

      enum class CoulombMode { NO_CORRECTION, FERMI_FUNCTION, EMA, MEMA,
        FERMI_AND_EMA, FERMI_AND_MEMA };

      CoulombCorrector( int pdg_c, int pdg_d,
        CoulombMode mode = CoulombMode::FERMI_AND_MEMA );

      double fermi_function( double beta_c ) const;

      double coulomb_correction_factor( double beta_rel_cd ) const;

      double ema_factor( double beta_rel_cd, bool& ok,
        bool modified_ema ) const;

      inline CoulombMode coulomb_mode() const
        { return coulomb_mode_; }

      void set_coulomb_mode( CoulombMode mode );

      static CoulombMode coulomb_mode_from_string( const std::string& str );

      static std::string string_from_coulomb_mode( CoulombMode mode );
 
    protected:

      static std::map< CoulombMode, std::string > coulomb_mode_string_map_;
 
      int pdg_c_; 
      double mc_; 
 
      int Zf_; 
      int Af_; 

      CoulombMode coulomb_mode_;
  };
 
}