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 232 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 226 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 257 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 238 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;
}