/* This is the header file for the functions that specify the
   parton distribution functions. */

#ifndef parton_distr_h
#define parton_distr_h

#include <default.h>
#include <int.h>
#include <matrix.h>


double up_valence(int param, double x, double scale_sqr);
double norm_up_valence(int param, int parton_npoints, double scale_sqr);
double int_up_valence(double x);
double mom_up_valence(int param, int parton_npoints, double scale_sqr);
double int_mom_up_valence(double x);
double moment_up_valence(int param, int parton_npoints, double scale_sqr,
                         double n);
double int_moment_up_valence(double x);

double down_valence(int param, double x, double scale_sqr);
double norm_down_valence(int param, int parton_npoints, double scale_sqr);
double int_down_valence(double x);
double mom_down_valence(int param, int parton_npoints, double scale_sqr);
double int_mom_down_valence(double x);
double moment_down_valence(int param, int parton_npoints, double scale_sqr,
			   double n);
double int_moment_down_valence(double x);

double up_bar(int param, double x, double scale_sqr);
double mom_up_bar(int param, int parton_npoints, double scale_sqr);
double int_mom_up_bar(double x);

double down_bar(int param, double x, double scale_sqr);
double mom_down_bar(int param, int parton_npoints, double scale_sqr);
double int_mom_down_bar(double x);

double strange_sea(int param, double x, double scale_sqr);
double mom_strange_sea(int param, int parton_npoints, double scale_sqr);
double int_mom_strange_sea(double x);

double sea_distr(int param, double x, double scale_sqr);
double diff_distr(int param, double x, double scale_sqr);
double norm_dbar_min_ubar(int param, int parton_npoints, double scale_sqr);
double int_norm_dbar_min_ubar(double x);
 
double pion_valence(int param, double x, double scale_sqr);
double norm_pion_valence(int param, int parton_npoints, double scale_sqr);
double int_pion_valence(double x);
double mom_pion_valence(int param, int parton_npoints, double scale_sqr);
double int_mom_pion_valence(double x);

double F2_distr(int param, double x, double scale_sqr);

double gluon_distr(int param, double x, double scale_sqr);

/**************************************************************************/
/**************************************************************************/

double pol_up_valence(int param, double x, double scale_sqr);
double norm_pol_up_valence(int param, int parton_npoints, double scale_sqr);
double int_pol_up_valence(double x);

double pol_down_valence(int param, double x, double scale_sqr);
double norm_pol_down_valence(int param, int parton_npoints, double scale_sqr);
double int_pol_down_valence(double x);

double spin_dilution(int param, double x);

double bjorken_sumrule(int param, int parton_npoints, double scale_sqr);
double int_bjorken_sumrule(double x);

double pol_strange_bar(int param, double x, double scale_sqr);
double pol_up_bar(int param, double x, double scale_sqr);
double pol_down_bar(int param, double x, double scale_sqr);

double norm_pol_strange_bar(int param, int parton_npoints, double scale_sqr);
double int_pol_strange_bar(double x);

/**************************************************************************/

double pol_strange_sea(int param, double x, double scale_sqr);
double norm_pol_strange_sea(int param, int parton_npoints, double scale_sqr);
double int_pol_strange_sea(double x);


#endif