softsusy::MssmSusy Class Reference

Contains all supersymmetric RPC-MSSM parameters. More...

#include <susy.h>

Inheritance diagram for softsusy::MssmSusy:

Public Member Functions
	MssmSusy ()
	MssmSusy (const MssmSusy &)
	Constructor sets object to be equal to another.
	MssmSusy (const DoubleMatrix &u, const DoubleMatrix &d, const DoubleMatrix &e, const DoubleVector &v, double m, double tb, double h, int mix)
const MssmSusy &	operator= (const MssmSusy &s)
	Default destructor. More...
void	setMssmSusy (const MssmSusy &s)
	sets object to be equal to another
void	setMssmLoops (double l)
void	setHvev (double h)
	Sets DRbar running Higgs vev.
void	setSomePars (const MssmSusy &s)
	Copies Yukawa matrices and gauge couplings from s only.
void	setYukawaElement (yukawa, int, int, double)
	Sets one element of a Yukawa matrix.
void	setYukawaMatrix (yukawa, const DoubleMatrix &)
	Sets whole Yukawa matrix.
void	setGaugeCoupling (int, double)
	Set a single gauge coupling.
void	setAllGauge (const DoubleVector &)
	Set all gauge couplings.
void	setSusyMu (double)
	Sets superpotential mu parameter.
void	setTanb (double)
	Sets tan beta.
void	setMssmApprox (int l, int t)
	Set loops/thresholds.
void	setMssmApprox (const Approx &a)
void	set (const DoubleVector &)
void	set (const DoubleVector &, int &k)
void	setMixing (double mix)
	Sets quark mixing parameter.
int	displayMssmLoops () const
double	displayHvev () const
	Returns DRbar running Higgs vev.
const MssmSusy &	displayMssmSusy () const
	Returns whole object as a const.
double	displayYukawaElement (yukawa, int, int) const
	Returns a single Yukawa matrix element.
const DoubleMatrix &	displayYukawaMatrix (yukawa) const
	Returns a whole Yukawa matrix.
double	displayGaugeCoupling (int) const
	Returns a single gauge coupling.
DoubleVector	displayGauge () const
	Returns all gauge couplings.
double	displaySusyMu () const
	Returns superpotential mu parameter.
double	displayTanb () const
	Returns tan beta.
Approx	displayMssmApprox () const
const DoubleVector	display () const
int	displayMixing () const
	Returns quark mixing parameter.
void	getMasses (QedQcd &r, double vev) const
void	quarkMixing (const DoubleMatrix &CKM, int mix)
void	setDiagYukawas (const QedQcd &, double vev)
void	getQuarkMixedYukawas (const QedQcd &r, const DoubleMatrix &CKM, int mix, double vev)
void	diagQuarkBasis (DoubleMatrix &vdl, DoubleMatrix &vdr, DoubleMatrix &vul, DoubleMatrix &vur) const
MssmSusy	beta (sBrevity &) const
void	getOneLpAnom (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gDD, DoubleMatrix &gUU, double &gH1H1, double &gH2H2, sBrevity &a) const
void	getTwoLpAnom (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gDD, DoubleMatrix &gUU, double &gH1H1, double &gH2H2, sBrevity &a) const
void	getThreeLpAnom (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gDD, DoubleMatrix &gUU, double &gH1H1, double &gH2H2, sBrevity &a) const
	By Bednyakov, see arXiv:1009.5455. More...
void	anomalousDimension (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, DoubleVector &dg, double &gH1H1, double &gH2H2, sBrevity &a) const

Detailed Description

Contains all supersymmetric RPC-MSSM parameters.

Constructor & Destructor Documentation

MssmSusy() [1/2]

softsusy::MssmSusy::MssmSusy

(

)

Constructor fills object with zeroes by default

MssmSusy() [2/2]

softsusy::MssmSusy::MssmSusy	(	const DoubleMatrix &	u,
		const DoubleMatrix &	d,
		const DoubleMatrix &	e,
		const DoubleVector &	v,
		double	m,
		double	tb,
		double	h,
		int	mix
	)

Constructor given Yukawa matrices u,d,e, gauge couplings v, mu parameter=m, tan beta=tb, renormalisation scale MU, number of loops in RG evolution l and thresholds parameter t

Member Function Documentation

anomalousDimension()

void softsusy::MssmSusy::anomalousDimension	(	DoubleMatrix &	gEE,
		DoubleMatrix &	gLL,
		DoubleMatrix &	gQQ,
		DoubleMatrix &	gUU,
		DoubleMatrix &	gDD,
		DoubleVector &	dg,
		double &	gH1H1,
		double &	gH2H2,
		sBrevity &	a
	)		const

Constants for gauge running

Constants not set yet

sBrevity a contains all of the shortcutted matrices etc;

For calculational brevity

1 loop contributions:

2 loop contributions:

3 loop contributions:

from http://www.liv.ac.uk/~dij/betas/allgennb.log only strong coupling constant + third family yukawa corrections are added Comments: rruiz

beta()

MssmSusy softsusy::MssmSusy::beta

(

sBrevity &

a

)

const

Outputs derivatives (DRbar scheme) in the form of ds. a contains the matrices calculated that are handy for computation. W= LL Y^E H1 ER + QL Y^D H1 DR + QL Y^U H2 UR + smu H2 H1 is the superpotential. Consistent with Allanach, Dedes, Dreiner hep-ph/9902251 and Barger, Berger and Ohmann hep-ph/9209232, 9311269 EXCEPT for the sign of smu, which is opposite. These equations are also valid for W= - LL Y^E H1 ER - QL Y^D H1 DR + QL Y^U H2 UR + smu H2 H1, the new SOFTSUSY convention

Wave function renormalisations: convention for g**(i, j) is that i is the LOWER index and j the upper in our paper hep-ph/9902251

keep this option in order to interface with RPVSUSY

To keep this a const function

contain derivatives of up, down quarks and leptons

mu parameter derivatives

RGEs of SUSY parameters

Following is from hep-ph/9308335: scalar H anomalous dimensions (as opposed to the chiral superfield one - see hep-ph/0111209). Additional contribution from Feynman gauge running at two-loops of tan beta: we need this to link up with BPMZ: hep-ph/0112251

1/(16 pi^2)^2

One-loop RGEs in Feynman gauge

Two-loop pieces

Contains all susy derivatives:

diagQuarkBasis()

void softsusy::MssmSusy::diagQuarkBasis	(	DoubleMatrix &	vdl,
		DoubleMatrix &	vdr,
		DoubleMatrix &	vul,
		DoubleMatrix &	vur
	)		const

Rotates to quark mass basis, returning the mixing matrices defined as \((Y_U)_{diag}\) = vul \(Y_U\) vur^+,
\((Y_D)_{diag}\) = vdl \(Y_D\) vdr^+
All matrices should be 3 by 3

display()

const DoubleVector softsusy::MssmSusy::display

(

)

const

Calculate beta functions of SUSY preserving parameters of RPC MSSM Returns all parameters as elements of a vector

getMasses()

void softsusy::MssmSusy::getMasses	(	QedQcd &	r,
		double	vev
	)		const

outputs object QedQcd & r valid at 1 GeV from SUSY data at mt, from diagonal elements of Yukawa couplings and Higgs VEV vev.

getOneLpAnom()

void softsusy::MssmSusy::getOneLpAnom	(	DoubleMatrix &	gEE,
		DoubleMatrix &	gLL,
		DoubleMatrix &	gQQ,
		DoubleMatrix &	gDD,
		DoubleMatrix &	gUU,
		double &	gH1H1,
		double &	gH2H2,
		sBrevity &	a
	)		const

Outputs one-loop anomlous dimensions gii given matrix inputs. for RH leptons, LH leptons, LH quarks, RH downs, RH ups, H1 and H2 respectively. Note that we use the convention (for matrices in terms of gamma's): gamma^Li_Lj = M_ij for LH fields and gamma^Rj_Ri = M_ij for RH fields (since they are really the complex conjugates of the RH fields). a should already be defined.

For calculational brevity

getQuarkMixedYukawas()

void softsusy::MssmSusy::getQuarkMixedYukawas	(	const QedQcd &	r,
		const DoubleMatrix &	CKM,
		int	mix,
		double	vev
	)

Defines mixed Yukawa matrices from data input in form of CKM matrix and r, vev. If mix=2, all mixing is in down sector mix=1, all mixing is in up sector

getThreeLpAnom()

void softsusy::MssmSusy::getThreeLpAnom	(	DoubleMatrix &	gEE,
		DoubleMatrix &	gLL,
		DoubleMatrix &	gQQ,
		DoubleMatrix &	gDD,
		DoubleMatrix &	gUU,
		double &	gH1H1,
		double &	gH2H2,
		sBrevity &	a
	)		const

By Bednyakov, see arXiv:1009.5455.

Outputs wave function renormalisation for SUSY parameters and gauge beta functions up to 2 loops. Also calculates and outputs a. IO parameters: RH leptons, LH leptons, LH quarks, RH downs, RH ups, H1 and H2 respectively. g^Li_Lj = m_{ij} for LH fields g^Ei_Ej = m_{ji} for RH fields Calculates three-loop anomalous dimensions in the dominant third family approximation and adds them

powers of gauge couplings

< 1/(16 pi^2)^3

< 6 Zeta(3)

full three family

For calculational brevity NB!!! Change notations to that of J&J hep-ph/0408128 (Y->Y^T , etc)

Tr(Y Y^T) = Tr (Y^T Y)

Everything gets the (1/16pi^2)^3 factor at the bottom

aij = ai^j, ai = gi^2

third family approximation

For calculational brevity NB Change notation to that of J&J hep-ph/0408128 (Y->Y^T , etc)

Everything gets the (1/16pi^2)^3 factor at the bottom

third family approximation

getTwoLpAnom()

void softsusy::MssmSusy::getTwoLpAnom	(	DoubleMatrix &	gEE,
		DoubleMatrix &	gLL,
		DoubleMatrix &	gQQ,
		DoubleMatrix &	gDD,
		DoubleMatrix &	gUU,
		double &	gH1H1,
		double &	gH2H2,
		sBrevity &	a
	)		const

Outputs two-loop anomlous dimensions gii given matrix inputs. for RH leptons, LH leptons, LH quarks, RH downs, RH ups, H1 and H2 respectively. Note that we use the convention (for matrices in terms of gamma's): gamma^Li_Lj = M_ij for LH fields and gamma^Rj_Ri = M_ij for RH fields (since they are really the complex conjugates of the RH fields). a should already be defined.

For calculational brevity

Everything gets the (1/16pi^2)^2 factor at the bottom

Two-loop pure gauge anom dimensions

Two-loop pure Yukawa contributions

1/(16 pi^2)^2

operator=()

const MssmSusy & softsusy::MssmSusy::operator=

(

const MssmSusy &

s

)

Default destructor.

sets object to be equal to another

quarkMixing()

void softsusy::MssmSusy::quarkMixing	(	const DoubleMatrix &	CKM,
		int	mix
	)

This turns diagonal Yukawa couplings at MZ into CKM mixed ones Takes diagonal quark Yukawa matrices and mixes them up according to the CKM matrix assuming: mix=2, all mixing is in down sector mix=1, all mixing is in up sector

setDiagYukawas()

void softsusy::MssmSusy::setDiagYukawas	(	const QedQcd &	r,
		double	vev
	)

Sets diagonal Yukawa couplings according to data in QedQcd input and Higgs VEV parameter vev= \(v_1^2+v_2^2\)

Higgs VEVs

---

The documentation for this class was generated from the following files:

• src/susy.h
• src/susy.cpp