OpticalModel.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 <iostream>
#include <complex>
 
namespace marley {

  class OpticalModel {
 
    public:

      OpticalModel( int Z, int A ) : Z_( Z ), A_( A ) {}
 
      virtual ~OpticalModel() = default;

      virtual std::complex< double > optical_model_potential( double r,
        double fragment_KE_lab, int fragment_pdg, int two_j, int l, int two_s,
        int target_charge = 0 ) = 0;

      virtual double transmission_coefficient( double total_KE_CM,
        int fragment_pdg, int two_j, int l, int two_s,
        int target_charge = 0 ) = 0;

      virtual double total_cross_section( double fragment_KE_lab,
        int fragment_pdg, int two_s, size_t l_max, int target_charge = 0 ) = 0;

      inline int Z() const;

      inline int A() const;

      virtual void print( std::ostream& out ) const = 0;
 
    protected:
 
      // Nuclear atomic and mass numbers
      int Z_, A_;
  };
 
  // Inline function definitions
  inline int OpticalModel::Z() const { return Z_; }
 
  inline int OpticalModel::A() const { return A_; }
}
 
inline std::ostream& operator<<( std::ostream& out,
  const marley::OpticalModel& om )
{
  om.print( out );
  return out;
}