GeneralTHDMMatching.cpp File Reference

Go to the source code of this file.

Classes
struct	__f_params

Functions
double	__fPS2_integrand (double x, void *p)
double	__fS2_integrand (double x, void *p)
double	gpseudoscalar2 (double r)
double	gscalar2 (double r)

Function Documentation

__fPS2_integrand()

double __fPS2_integrand	(	double	x,
		void *	p
	)

Definition at line 123 of file GeneralTHDMMatching.cpp.

                                          {
  __f_params &params= *reinterpret_cast<__f_params *>(p);
  double integ = log(x*(1.0-x)/params.a) / (x*(1.0-x)-params.a);
  return integ;
}

__fS2_integrand()

double __fS2_integrand	(	double	x,
		void *	p
	)

Definition at line 117 of file GeneralTHDMMatching.cpp.

                                         {
  __f_params &params= *reinterpret_cast<__f_params *>(p);
  double integ = (2.0*x*(1.0-x)-1.0) * log(x*(1.0-x)/params.a) / (x*(1.0-x)-params.a);
  return integ;
}

gpseudoscalar2()

double gpseudoscalar2

(

double

r

)

Definition at line 150 of file GeneralTHDMMatching.cpp.

{
    __f_params params;
    params.a=r;
 
    double result;
    gsl_integration_glfixed_table * w
        = gsl_integration_glfixed_table_alloc(100);
    gsl_function F;
 
    F.function = &__fPS2_integrand;
    F.params = reinterpret_cast<void *>(&params);
 
    result = gsl_integration_glfixed (&F, 0, 1, w);
 
    gsl_integration_glfixed_table_free (w);
 
  return result;
}

gscalar2()

double gscalar2

(

double

r

)

Definition at line 130 of file GeneralTHDMMatching.cpp.

{
    __f_params params;
    params.a=r;
 
    double result;
    gsl_integration_glfixed_table * w
        = gsl_integration_glfixed_table_alloc(100);
    gsl_function F;
 
    F.function = &__fS2_integrand;
    F.params = reinterpret_cast<void *>(&params);
 
    result = gsl_integration_glfixed (&F, 0, 1, w);
 
    gsl_integration_glfixed_table_free (w);
    
  return result;
}