/* This file calculates the eta wavefunction. */

/**************************************************************************/
/*                                                                        */
/* IMPORTANT : all dimensional quantities in this program                 */
/*             are in units GeV                                           */  
/*                                                                        */
/**************************************************************************/


#include <eta_wavefunction.h>

/* the function eta_distribution_ampl() returns the 
   asymptotic eta distribution amplitude 
   as function of the quark longitudinal momentum fraction x */

double eta_distribution_ampl(double x)
{
     return 1./(2. * sqrt(6.)) * F_ETA * 6. * x * (1. - x);     
}


/* The function eta_wavef_impactpar() returns the eta wavefunction 
   as function of the quark longitudinal momentum fraction x and 
   as function of the transverse distance b (in GeV^-1)
   between the quark and antiquark. 
   To model the transverse momentum dependence, a gaussian ansatz is used. */

double eta_wavef_impactpar(double x, double b)
{
     double eta_avmom_sqr;

     eta_avmom_sqr = 8. * pow(PI * F_ETA,2.);  /* average squared transverse 
                                                quark momentum in eta 
                                                wavefunction in units GeV^2 */
   
     return 4. * PI * eta_distribution_ampl(x) 
            * exp(- .5 * eta_avmom_sqr * x * (1. - x) * pow(b,2.));        
}


/* The function eta_wavef_momspace() returns the eta wavefunction 
   in momentum space as function of 
   the quark longitudinal momentum fraction x and as function of 
   the transverse momentum k_perp between the quark and antiquark. 
   To model the transverse momentum dependence, a gaussian ansatz is used. */

double eta_wavef_momspace(double x, double k_perp)
{
     double eta_avmom_sqr;

     eta_avmom_sqr = 8. * pow(PI * F_ETA,2.);  /* average squared transverse 
                                                quark momentum in eta 
                                                wavefunction in units GeV^2 */
 
     return 8. * pow(PI,2.) / eta_avmom_sqr
            * eta_distribution_ampl(x) / (x * (1. - x)) 
            * exp(- .5 * pow(k_perp,2.)/(eta_avmom_sqr * x * (1. - x)) );
        
}