Contains all supersymmetric RPC-MSSM parameters and RGEs. More...
#include <susy.h>
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 MU, int l, int t, double h) | |
const MssmSusy & | operator= (const MssmSusy &s) |
Default destructor. | |
void | setSusy (const MssmSusy &s) |
sets object to be equal to another | |
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 | set (const DoubleVector &) |
Sets all RGE parameters to elements of vector. | |
void | setTanb (double) |
Sets tan beta. | |
double | displayHvev () const |
Returns DRbar running Higgs vev. | |
const MssmSusy & | displaySusy () 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. | |
const DoubleVector | display () const |
Returns all parameters as elements of a vector. | |
double | displayTanb () const |
Returns tan beta. | |
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) |
DoubleVector | beta () const |
Calculate beta functions of SUSY preserving parameters of RPC MSSM. | |
MssmSusy | beta (sBrevity &) const |
Calculate beta functions of SUSY preserving parameters of RPC MSSM. | |
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. | |
void | anomalousDimension (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, DoubleVector &dg, double &gH1H1, double &gH2H2, sBrevity &a) const |
void | diagQuarkBasis (DoubleMatrix &vdl, DoubleMatrix &vdr, DoubleMatrix &vul, DoubleMatrix &vur) const |
Contains all supersymmetric RPC-MSSM parameters and RGEs.
Constructor fills object with zeroes by default
MssmSusy::MssmSusy | ( | const DoubleMatrix & | u, |
const DoubleMatrix & | d, | ||
const DoubleMatrix & | e, | ||
const DoubleVector & | v, | ||
double | m, | ||
double | tb, | ||
double | MU, | ||
int | l, | ||
int | t, | ||
double | h | ||
) |
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
void MssmSusy::anomalousDimension | ( | DoubleMatrix & | gEE, |
DoubleMatrix & | gLL, | ||
DoubleMatrix & | gQQ, | ||
DoubleMatrix & | gUU, | ||
DoubleMatrix & | gDD, | ||
DoubleVector & | dg, | ||
double & | gH1H1, | ||
double & | gH2H2, | ||
sBrevity & | a | ||
) | const |
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
MssmSusy MssmSusy::beta | ( | sBrevity & | a | ) | const |
Calculate beta functions of SUSY preserving parameters of RPC MSSM.
One-loop RGEs in Feynman gauge
Two-loop pieces
void MssmSusy::diagQuarkBasis | ( | DoubleMatrix & | vdl, |
DoubleMatrix & | vdr, | ||
DoubleMatrix & | vul, | ||
DoubleMatrix & | vur | ||
) | const |
Rotates to quark mass basis, returning the mixing matrices defined as = vul
vur^+,
= vdl
vdr^+ All matrices should be 3 by 3
void 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.
void 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.
void 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
void 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.
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)
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
void 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.
Default destructor.
sets object to be equal to another
void 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
void MssmSusy::setDiagYukawas | ( | const QedQcd & | r, |
double | vev | ||
) |
Sets diagonal Yukawa couplings according to data in QedQcd input and Higgs VEV parameter vev=