SHRiMPS/Eikonals/Eikonal_Creator.H

0001 #ifndef SHRIMPS_Eikonals_Eikonal_Creator_H
0002 #define SHRIMPS_Eikonals_Eikonal_Creator_H
0003
0004 #include "SHRiMPS/Eikonals/Omega_ik.H"
0005 #include "SHRiMPS/Eikonals/Form_Factors.H"
0006
0007 namespace SHRIMPS {
0008   /*!
0009     \class Eikonal_Creator
0010     \brief Initialises and fills the single channel eikonals,
0011     \f$\Omega_{ik}\f$, represented by SHRIMPS::Omega_ik.
0012
0013     To this end, the eikonal creator solves the coupled differential
0014     equations for the two terms \f$Omega_{i(k)}\f$ and \f$\Omega_{(i)k}\f$,
0015     individually implemented in the class SHRIMPS::Eikonal_Contributor.
0016     Having them at hand, this class also convolutes them by integrating out
0017     the individual impact parameters \f$\vec b_\perp^{(1,2)}\f$.
0018   */
0019
0020   class Eikonal_Creator {
0021   private:
0022     Form_Factor      * p_ff1, * p_ff2;
0023     Eikonal_Parameters m_params;
0024     int                m_Bsteps,m_ff1steps, m_ff2steps;
0025
0026     void FillBYGrids(Eikonal_Contributor * omegai,
0027                Eikonal_Contributor * omegak);
0028     void CreateImpactParameterGrid(Omega_ik * Omega_ik);
0029     void FixGridAndBorders(DEQ_Solver * solver,int & steps,
0030                const double & ff1,const double & ff2);
0031
0032   public:
0033     Eikonal_Creator();
0034     ~Eikonal_Creator() {}
0035
0036     void       SetFormFactors(Form_Factor * ff1,Form_Factor * ff2);
0037     Omega_ik * InitialiseEikonal();
0038   };
0039 }
0040
0041 #endif