Gambit::FlavBit::LoopFunctions Namespace Reference

Functions
double	G1 (const double x)
double	G1 (const double a, const double b, const double c)
double	MFVV (const double a, const double b)
double	B1 (const double a, const double b, const double Q)
double	B0 (const double a, const double b, const double Q)
double	C0 (const double a, const double b, const double c)
double	C00 (const double a, const double b, const double c, const double Q)
double	B02C00C0 (const double a, const double b, const double c, const double Q)
double	D0 (const double a, const double b, const double c, const double d)
double	D27 (const double a, const double b, const double c, const double d)
double	IC0D0 (const double a, const double b, const double c, const double d)
double	E (const double x, const double y)

Function Documentation

B0()

double Gambit::FlavBit::LoopFunctions::B0	(	const double	a,
		const double	b,
		const double	Q
	)

Definition at line 69 of file flav_loop_functions.hpp.

References Gambit::Scanner::pow().

Referenced by B02C00C0(), mt2_bisect::mt2::nsols(), mt2_bisect::mt2::nsols_massless(), mt2w_bisect::mt2w::teco(), Gambit::FlavBit::Penguins::VZw10LL(), and Gambit::FlavBit::Penguins::VZw8LL().

      {
        // TODO: behaviour when a = 0 and b = 0 undefined
        if(a == 0 and b == 0)
          return 0;
        else if(a == b)
          return -log(b / pow(Q,2));
        else if(a == 0)
          return 1. - log(b / pow(Q,2));
        else if(b == 0)
          return 1. - log(a) - log(1./pow(Q,2));
        else
          return 1. - log(b / pow(Q,2)) + 1./(b-a) * a * log(a/b);
      }

Here is the call graph for this function:

Here is the caller graph for this function:

B02C00C0()

double Gambit::FlavBit::LoopFunctions::B02C00C0	(	const double	a,
		const double	b,
		const double	c,
		const double	Q
	)

Definition at line 149 of file flav_loop_functions.hpp.

References B0(), C0(), C00(), and Gambit::Scanner::pow().

Referenced by Gambit::FlavBit::Penguins::VZw10LL().

      {
        if(a == 0 and b == 0)
          return 0.25*(1.0 - 2.0*log(c) - 2.0*log(1 / pow(Q,2)));
        else
          return B0(a,b,Q) - 2*C00(a,b,c,Q) + C0(a,b,c)*c;
      }

Here is the call graph for this function:

Here is the caller graph for this function:

B1()

double Gambit::FlavBit::LoopFunctions::B1	(	const double	a,
		const double	b,
		const double	Q
	)

Definition at line 59 of file flav_loop_functions.hpp.

References Gambit::Scanner::pow().

Referenced by mt2_bisect::mt2::nsols(), mt2_bisect::mt2::nsols_massless(), Gambit::FlavBit::Penguins::Shw2w4LL(), Gambit::FlavBit::Penguins::Shw2w4RR(), mt2w_bisect::mt2w::teco(), Gambit::FlavBit::Penguins::VZw2w4LL(), and Gambit::FlavBit::Penguins::VZw2w4RR().

      {
        if(a == b)
          return 0.5 * log(b / pow(Q,2));
        else if(a == 0)
          return -0.25 + 0.5*log(b / pow(Q,2));
        else
          return -0.5 + 0.5*log(b / pow(Q,2)) - (a*a - b*b + 2.*a*a*log(b/a)) / (4.*pow(a-b,2));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

C0()

double Gambit::FlavBit::LoopFunctions::C0	(	const double	a,
		const double	b,
		const double	c
	)

Definition at line 84 of file flav_loop_functions.hpp.

References b, and Gambit::Scanner::pow().

Referenced by B02C00C0(), IC0D0(), mt2_bisect::mt2::nsols(), mt2_bisect::mt2::nsols_massless(), mt2w_bisect::mt2w::teco(), Gambit::FlavBit::Penguins::VZw10LL(), and Gambit::FlavBit::Penguins::VZw8LL().

      {
        // TODO: behaviour when two paramers are 0 undefined, set it to zero
        if(a == 0 and b == 0 and c == 0)
          return 0;
        if(a == 0 and b == 0) 
          return 0;
        else if(a == 0 and c == 0) 
          return 0;
        else if(b == 0 and c == 0)
          return 0;
        else if(c == 0)
          return C0(a,c,b);
        else if(a == b and b == c) 
          return - 1./(2*c);
        else if(a == b)
          return (-b + c- c*log(c/b)) / pow(b-c,2);
        else if(a == c and b != 0)
          return C0(a,c,b);
        else if(a == c and b == 0)
          return -1./c;
        else if(b == c and a != 0)
          return (a - c + a*log(c/a)) / pow(a-c,2);
        else if(b == c and a == 0)
          return -1./c;
        else if(a == 0)
          return (-log(b) + log(c)) / (b-c);
        else if(b == 0)
          return log(c/a)/(a-c);
        else
          return -1. / (a-b)*(a-c)*(b-c)*( b*(c-a)*log(b/a) + c*(a-b)*log(c/a));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

C00()

double Gambit::FlavBit::LoopFunctions::C00	(	const double	a,
		const double	b,
		const double	c,
		const double	Q
	)

Definition at line 117 of file flav_loop_functions.hpp.

References b, and Gambit::Scanner::pow().

Referenced by B02C00C0(), Gambit::FlavBit::Penguins::VZw10LL(), and Gambit::FlavBit::Penguins::VZw8LL().

      {
        // TODO: behaviour when all three parameters are zero is undefined, set it to zero
        if(a == 0 and b == 0 and c == 0)
          return 0;
        else if(b == 0 and c == 0)
          return 0.125*(3. - 2.*log(a/pow(Q,2)));
        else if(a == 0 and b == 0)
          return 0.125*(3. - 2.*log(c) - 2.*log(1./pow(Q,2)));
        else if(c == 0)
          return C00(a,c,b,Q);
        else if(a == b and b == c)
          return -0.25*log(c/pow(Q,2));
        else if(a == b)
          return - (2.*c*c*log(c/b) + (b-c)*(-b + 3.*c + 2.*(b-c)*log(b/pow(Q,2))))/(8.*pow(b-c,2));
        else if(a == c and b != 0)
          return C00(a,c,b,Q);
        else if(a == c and b == 0)
          return 0.125*(1. - 2.*log(c/pow(Q,2)));
        else if(b == c and a != 0)
          return (2.*(2.*a-c)*c*log(c/a)-(a-c)*(-3.*a+c+2.*(a-c)*log(a/pow(Q,2))))/(8.*pow(a-c,2));
        else if(b == c and a == 0)
          return 0.125*(1. - 2.*log(c) - 2.*log(1./pow(Q,2)));
        else if(a == 0)
          return -(2.*b*log(b) - 2.*c*log(c) + (b-c)*(-3.+2.*log(1./pow(Q,2))))/(8.*(b-c));
        else if(b == 0)
          return (2.*c*log(c/a) - (a-c)*(-3. + 2.*log(a/pow(Q,2))))/(8.*(a-c));
        else
          return 1. / (8.*(a-b)*(a-c)*(b-c)) * ( (c-a)*((a-b)*(2.*log(a/pow(Q,2))-3.)*(b-c) - 2.*b*b*log(b/a)) + 2.*c*c*(b-a)*log(c/a));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

D0()

double Gambit::FlavBit::LoopFunctions::D0	(	const double	a,
		const double	b,
		const double	c,
		const double	d
	)

Definition at line 157 of file flav_loop_functions.hpp.

References b, and Gambit::Scanner::pow().

Referenced by IC0D0(), mt2_bisect::mt2::nsols(), mt2_bisect::mt2::nsols_massless(), mt2w_bisect::mt2w::teco(), and Gambit::FlavBit::Boxes::Vw8lLL().

      {
        //TODO: behaviour when two or more parameters are zero is undefined, set it to zero
        if((!a and !b) or (!b and !c) or (!b and !d) or (!c and !d))
          return 0;
        else if(c == 0)
          return D0(a,c,b,d);
        else if(d == 0)
          return D0(a,d,c,b);
        else if(a == b and b == c and c == d)
          return 1. / (6.*d*d);
        else if(a == b and b == c)
          return D0(a,d,c,d);
        else if(a == b and b == d)
          return D0(a,c,b,d);
        else if(a == c and c == d and b == 0)
          return 1. / (2.*c*c);
        else if(a == c and c == d and b != 0)
          return (-b*b + c*c + 2.*b*c*log(b/c)) / (2.*c*pow(c-b,3));
        else if(b == c and c == d and a == 0)
          return 1. / (2.*d*d);
        else if(b == c and c == d and a != 0)
          return (a*a - d*d + 2.*a*d*log(d/a)) / (2.*d*pow(a-d,3));
        else if(a == d and b == c)
          return (-2*c + 2*d + (c+d)*log(c/d)) / pow(c-d,3);
        else if(a == d and b == 0)
          return (c - d -d*log(c/d)) / (pow(c-d,2)*d);
        else if(a == d)
          return 1./ ((b-d)*(d-c))-(b*log(b/d))/((b-c)*pow(b-d,2))+(c*log(c/d))/((b-c)*pow(c-d,2));
        else if(a == c)
          return D0(a,b,d,c);
        else if(a == b)
          return D0(a,d,c,b);
        else if(b == c)
          return D0(a,d,c,b);
        else if(b == d)
          return D0(a,c,b,d);
        else if(c == d and b == 0)
          return (a - d + d*log(d/a)) / (d*pow(a-d,2));
        else if(c == d and a == 0)
          return (b - d + d*log(d/b)) / (d*pow(b-d,2));
        else if(c == d)
          return (b*pow(a-d,2)*log(b/a) - (a-b)*( (a-d)*(b-d) + (a*b-d*d)*log(d/a) )) / ((a-b)*pow(a-d,2)*pow(b-d,2));
        else if(b == 0)
          return log(c/a)/((a-c)*(c-d)) + log(d/a)/((a-d)*(d-c));
        else if(a == 0)
          return ((d-c)*log(b) + (b-d)*log(c) + (c-b)*log(d))/((b-c)*(b-d)*(c-d));
        else
          return -(b*log(b/a)/((b-a)*(b-c)*(b-d)) + c*log(c/a)/((c-a)*(c-b)*(c-d)) + d*log(d/a)/((d-a)*(d-b)*(d-c)));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

D27()

double Gambit::FlavBit::LoopFunctions::D27	(	const double	a,
		const double	b,
		const double	c,
		const double	d
	)

Definition at line 208 of file flav_loop_functions.hpp.

References b, and Gambit::Scanner::pow().

Referenced by Gambit::FlavBit::Boxes::Vw4dLL(), Gambit::FlavBit::Boxes::Vw4lLL(), and Gambit::FlavBit::Boxes::Vw4uLL().

      {
        //TODO: behaviour when three or more parameters are zero is undefined, set it to zero
        if((!a and !b and !c) or (!a and !b and !d) or (!a and !c and !d) or (!b and !c and !d))
          return 0; 
        if(a == b and b == c and c == d)
           return -1./(12.*d);
        if(a == d and c == d and b == 0)
           return -1. / (8.*d);
        if(a == d and c == d)
           return (3.*b*b - 4.*b*d + d*d - 2.*b*b*log(b/d))/(8.*pow(b-d,3));
        if(b == c and c == d)
           return D27(b,a,c,d);
        if(a == b and b == c)
           return D27(a,d,c,b);
        if(a == b and b == d)
           return D27(a,c,b,d);
        if(a == b and c == d and b == 0)
          return -1./(4.*d);
        if(a == b and c == d and d == 0)
          return -1./(4.*b);
        if(a == b and c == d) 
          return (-b*b + d*d -2.*b*d*log(d/b)) / (4.*pow(b-d,3));
        if(a == c and b == d)
          return D27(a,c,b,d);
        if(a == d and b == c)
          return D27(a,b,d,c);
        if(a == b and b == 0)
          return log(d/c)/(4.*(c-d));
        if(a == b and c == 0)
          return - (b -d +d*log(d/b))/(4.*pow(b-d,2));
        if(a == b and d == 0)
          return D27(a,b,d,c);
        if(a == b)
          return 0.25*(-c*c*log(c/b)/(pow(b-c,2)*(c-d)) + (b*(d-b)/(b-c) + d*d*log(d/b)/(c-d))/pow(b-d,2));
        if(a == c)
          return D27(a,c,b,d);
        if(a == d)
          return D27(a,d,c,b);
        if(b == c and a == 0)
          return (-c+d+d*log(c/d))/(4.*pow(c-d,2));
        if(b == c and c == 0)
          return log(d/a)/(4.*(a-d));
        if(b == c and d == 0)
          return (a-c+a*log(c/a))/(4.*pow(a-c,2));
        if(b == c)
          return (c*(a-d)*(a*(c-2.*d)+c*d)*log(c/a)+(a-c)*(c*(a-d)*(c-d)+(a-c)*d*d*log(d/a)))/(4.*pow(a-c,2)*(a-d)*pow(c-d,2));
        if(b == d)
          return D27(a,b,d,c);
        if(c == d)
          return D27(a,c,d,b);
        if(a == 0)
          return (b*(-c+d)*log(b)+c*(b-d)*log(c)+(-b+c)*d*log(d))/(4.*(b-c)*(b-d)*(c-d));
        if(b == 0)
          return ((c*log(c/a))/((a - c)*(c - d)) + (d*log(d/a))/((a - d)*(-c + d)))/4.;
        if(c == 0)
          return D27(a,c,b,d);
        if(d == 0)
          return D27(a,d,c,b);
        else
          return -0.25*(b*b*log(b/a)/((b-a)*(b-c)*(b-d)) + c*c*log(c/a)/((c-a)*(c-b)*(c-d)) + d*d*log(d/a)/((d-a)*(d-b)*(d-c)));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

E()

double Gambit::FlavBit::LoopFunctions::E	(	const double	x,
		const double	y
	)

Definition at line 280 of file flav_loop_functions.hpp.

References Gambit::Scanner::pow().

Referenced by ben_findsols(), Gambit::DarkBit::boost_dNdE(), Gambit::DarkBit::calc_ALPS1_signal_gas(), Gambit::DarkBit::calc_AxionOscillationTemperature(), Gambit::DarkBit::calc_Haloscope_signal(), Gambit::ColliderBit::calcMT_1l(), Gambit::DarkBit::cascadeMC_fetchSpectra(), Gambit::DarkBit::cascadeMC_gammaSpectra(), Gambit::DarkBit::cascadeMC_Histograms(), Gambit::DarkBit::DecayChain::ChainParticle::E_Lab(), Gambit::DarkBit::DecayChain::Ep4vec(), Gambit::DarkBit::getYield(), Gambit::DarkBit::intersect_parabola_line(), Gambit::PrecisionBit::lnL_GF(), Gambit::DarkBit::RD_oh2_Axions(), Gambit::FlavBit::RHN_RK(), Gambit::FlavBit::RHN_RKstar_0045_11(), Gambit::FlavBit::RHN_RKstar_11_60(), Gambit::DarkBit::SimpleHist::SimpleHist(), Gambit::SMInputs::SMInputs(), Gambit::ColliderBit::sortByPT(), Gambit::DecayBit::stau_1_decays_smallsplit(), and Gambit::DarkBit::TH_ProcessCatalog_WIMP().

      {
        if(x == 0 or y == 0)
          return 0.0;
        if(x == y)
          return (x*(-4.0 + 15.0*x - 12.0*pow(x,2) + pow(x,3) + 6.0*pow(x,2)*log(x)))/ (4.*pow(-1.0 + x,3));
        return x*y*(-3.0/(4.0*(1.0-x)*(1.0 - y)) + ((0.25 - 3.0/(4.0*pow(-1.0 + x,2)) - 3.0/(2.0*(-1.0 + x)))*log(x))/(x - y) + ((0.25 - 3.0/(4.0*pow(-1.0 + y,2)) - 3.0/(2.0*(-1 + y)))*log(y))/(-x + y));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

G1() [1/2]

double Gambit::FlavBit::LoopFunctions::G1

(

const double

x

)

Definition at line 31 of file flav_loop_functions.hpp.

References Gambit::Scanner::pow().

Referenced by Gambit::FlavBit::Penguins::A2L(), and Gambit::FlavBit::Penguins::A2R().

      {
        if(x == 0)
          return -7./12.;
        if(x == 1)
          return -5./12.;
        else
          return (-7. + 33.*x - 57.*x*x + 31.*x*x*x + 6.*x*x*(1. - 3*x)*log(x))/(12.*pow(1.-x,4));
      }

Here is the call graph for this function:

Here is the caller graph for this function:

G1() [2/2]

double Gambit::FlavBit::LoopFunctions::G1	(	const double	a,
		const double	b,
		const double	c
	)

Definition at line 41 of file flav_loop_functions.hpp.

References b.

      {
        if(b == c and b != 0)
          return G1(a/b)/b;
        else
          return 0; // TODO: 2C12 + 2C22 - C1 or 2C12 + C11 - C2
      }

IC0D0()

double Gambit::FlavBit::LoopFunctions::IC0D0	(	const double	a,
		const double	b,
		const double	c,
		const double	d
	)

Definition at line 271 of file flav_loop_functions.hpp.

References C0(), and D0().

Referenced by Gambit::FlavBit::Boxes::Vw4dLL(), and Gambit::FlavBit::Boxes::Vw4lLL().

      {
        return C0(a,b,c) + d*D0(a,b,c,d);
      }

Here is the call graph for this function:

Here is the caller graph for this function:

MFVV()

double Gambit::FlavBit::LoopFunctions::MFVV	(	const double	a,
		const double	b
	)

Definition at line 49 of file flav_loop_functions.hpp.

References b, and Gambit::Scanner::pow().

Referenced by Gambit::FlavBit::Penguins::A1R().

      {
        if(a == b)
          return 1. / (3. * b);
        else if(a == 0)
          return 5. / (9. * b);
        else
          return (6.*a*a*(a-3.*b)*log(a/b) - (a-b)*(5.*a*a - 22.*a*b + 5.*b*b))/(9.*pow(a-b,4));
      }

Here is the call graph for this function:

Here is the caller graph for this function: