a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
THDMMatching.cpp File Reference

Go to the source code of this file.

Classes

struct  __f_params
 

Functions

double __fPS_integrand (double x, void *p)
 
double __fS_integrand (double x, void *p)
 
double gpseudoscalar (double r)
 
double gscalar (double r)
 

Function Documentation

◆ __fPS_integrand()

double __fPS_integrand ( double  x,
void *  p 
)

Definition at line 499 of file THDMMatching.cpp.

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

◆ __fS_integrand()

double __fS_integrand ( double  x,
void *  p 
)

Definition at line 493 of file THDMMatching.cpp.

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

◆ gpseudoscalar()

double gpseudoscalar ( double  r)

Definition at line 526 of file THDMMatching.cpp.

527{
528 __f_params params;
529 params.a=r;
530
531 double result;
532 gsl_integration_glfixed_table * w
533 = gsl_integration_glfixed_table_alloc(100);
534 gsl_function F;
535
536 F.function = &__fPS_integrand;
537 F.params = reinterpret_cast<void *>(&params);
538
539 result = gsl_integration_glfixed (&F, 0, 1, w);
540
541 gsl_integration_glfixed_table_free (w);
542
543 return result;
544}
double __fPS_integrand(double x, void *p)

◆ gscalar()

double gscalar ( double  r)

Definition at line 505 of file THDMMatching.cpp.

506{
507 __f_params params;
508 params.a=r;
509
510 double result;
511 gsl_integration_glfixed_table * w
512 = gsl_integration_glfixed_table_alloc(100);
513 gsl_function F;
514
515 F.function = &__fS_integrand;
516 F.params = reinterpret_cast<void *>(&params);
517
518 result = gsl_integration_glfixed (&F, 0, 1, w);
519
520 gsl_integration_glfixed_table_free (w);
521
522 return result;
523}
double __fS_integrand(double x, void *p)