#ifndef SHRIMPS_Cross_Sections_Cross_Sections_H
#define SHRIMPS_Cross_Sections_Cross_Sections_H

#include "SHRiMPS/Cross_Sections/Sigma_Inelastic.H"
#include "SHRiMPS/Cross_Sections/Sigma_Elastic.H"
#include "SHRiMPS/Cross_Sections/Sigma_SD.H"
//#include "SHRiMPS/Cross_Sections/Sigma_DD.H"
#include "SHRiMPS/Eikonals/Omega_ik.H"
#include <list>
#include <map>

namespace SHRIMPS {
  class Cross_Sections {
  private:
    Omega_ik      * p_selected;       
    double          m_xstot, m_slope, m_xsinel, m_xsel, m_xsSD[2], m_xsDD;
    Sigma_Inelastic m_sigma_inelastic;
    Sigma_Elastic   m_sigma_elastic;
    Sigma_SD        m_sigma_SD;
    std::map<run_mode::code,double> m_modemap;
  public:
    Cross_Sections();
    ~Cross_Sections();
    
    void CalculateCrossSections();
    void Test(const std::string & dirname);
    //run_mode::code SelectCollisionMode();

    //std::list<Omega_ik *> * GetEikonals() { return p_eikonals; }    

    Sigma_Inelastic * GetSigmaInelastic() { return (&m_sigma_inelastic); }
    Sigma_Elastic   * GetSigmaElastic()   { return (&m_sigma_elastic); }
    Sigma_SD        * GetSigmaSD()        { return (&m_sigma_SD); }

    const double & SigmaTot()        const { return m_xstot; }
    const double & SigmaSD(size_t i) const { return m_xsSD[i]; }
    const double & SigmaDD()         const { return m_xsDD; }
    const double & SigmaEl()         const { return m_xsel; }
    const double & SigmaInel()       const { return m_xsinel; }
    const double & ElasticSlope()    const { return m_slope; }
  };
}
#endif