SOFTSUSY
4.1
|
Soft SUSY breaking parameters and beta functions. More...
#include <softpars.h>
Public Member Functions | |
MssmSoftPars () | |
Default constructor fills object with zeroes. | |
MssmSoftPars (const MssmSoftPars &) | |
Constructor sets all parameters equal to those in another object. More... | |
const MssmSoftPars & | operator= (const MssmSoftPars &s) |
Sets all parameters equal to those in another object. | |
MssmSoftPars (const DoubleVector &mG, const DoubleMatrix &aU, const DoubleMatrix &aD, const DoubleMatrix &aE, const DoubleMatrix &mQl, const DoubleMatrix &mUr, const DoubleMatrix &mDr, const DoubleMatrix &mLl, const DoubleMatrix &mEr, double m3sq, double mH1sq, double mH2sq, double mGravitino) | |
const MssmSoftPars & | displayMssmSoftPars () const |
Returns whole object as a const. | |
const DoubleMatrix & | displayTrilinear (trilinears) const |
Return a trilinear coupling matrix. | |
double | displayTrilinear (trilinears, int i, int j) const |
Return a trilinear element. | |
const DoubleMatrix & | displaySoftMassSquared (softMasses) const |
Return a soft mass squared matrix. | |
double | displaySoftMassSquared (softMasses, int i, int j) const |
Return a soft mass squared element. | |
double | displayGravitino () const |
Returns the gravitino mass. | |
double | displayM3Squared () const |
Return \( m_3^2\). | |
double | displayMh1Squared () const |
Return \(m_{H_1}^2\). | |
double | displayMh2Squared () const |
Return \(m_{H_2}^2\)=mH2sq. | |
DoubleVector | displayGaugino () const |
Return \(M_{G_i}\). | |
double | displayGaugino (int i) const |
void | display (DoubleVector &v, int &k) const |
const DoubleVector | display2 () const |
void | setM32 (double) |
Sets gravitino mass. | |
void | setSoftPars (MssmSoftPars const &) |
Sets whole thing equal to another object. | |
void | setSoftMassElement (softMasses, int, int, double) |
Set one element of a soft mass squared matrix. | |
void | setSoftMassMatrix (softMasses, const DoubleMatrix &) |
Set whole of a soft mass squared matrix. | |
void | setTrilinearMatrix (trilinears, const DoubleMatrix &) |
Set whole of a trilinear SUSY breaking parameter matrix. | |
void | setTrilinearElement (trilinears k, int i, int j, double a) |
Set one element of a trilinear SUSY breaking parameter matrix. | |
void | setGauginoMass (int, double) |
Set one gaugino mass. | |
void | setAllGauginos (const DoubleVector &) |
Set all gaugino masses. | |
void | setM3Squared (double) |
Sets \( m_3^2\). | |
void | setMh1Squared (double) |
Sets \( m_{H_1}^2\). | |
void | setMh2Squared (double) |
void | set (const DoubleVector &y) |
Sets all SUSY breaking parameters only - not SUSY ones. | |
DoubleVector | beta (const MssmSusy &) const |
Returns double vector containing numerical beta functions of parameters. | |
MssmSoftPars | beta2 (const MssmSusy &xx) const |
Returns numerical beta functions of parameters | |
MssmSoftPars | beta2 (const MssmSusy &xx, sBrevity &a) const |
Returns numerical beta functions of parameters and Brevity. More... | |
void | anomalousDeriv (const MssmSusy &xx, DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, double &gH1H1, double &gH2H2) const |
void | yTildes (const MssmSusy &xx, DoubleMatrix &yu, DoubleMatrix &yd, DoubleMatrix &ye) const |
Ytilde quantities are for calculational brevity in beta functions. | |
void | u1R_PQflip (MssmSusy &s) |
void | universal (MssmSusy &s, double m0, double m12, double a0, double mu, double m3sq) |
void | addAmsb (const MssmSusy &, double m32) |
void | universalTrilinears (const MssmSusy &a, double a0) |
Sets all SUSY breaking trilinear couplings to a0. | |
void | minimalGmsb (const MssmSusy &xx, int n5, double LAMBDA, double mMess, double cgrav) |
void | standardSugra (const MssmSusy &xx, double m0, double m12, double a0) |
Reads in universal boundary conditions at the current scale: m0, M1/2, A0. | |
void | universalScalars (double m0) |
Sets all flavour-diagonal SUSY breaking scalar masses to m0. | |
void | universalGauginos (double m12) |
Sets all flavour-diagonal SUSY breaking gaugino masses to m12. | |
void | inputSoftParsOnly () |
Reads in soft SUSY breaking parameters from a file. | |
Soft SUSY breaking parameters and beta functions.
softsusy::MssmSoftPars::MssmSoftPars | ( | const MssmSoftPars & | s | ) |
Constructor sets all parameters equal to those in another object.
Constructor fills SUSY conserving parts with another object, all SUSY breaking parameters set to zero
softsusy::MssmSoftPars::MssmSoftPars | ( | const DoubleVector & | mG, |
const DoubleMatrix & | aU, | ||
const DoubleMatrix & | aD, | ||
const DoubleMatrix & | aE, | ||
const DoubleMatrix & | mQl, | ||
const DoubleMatrix & | mUr, | ||
const DoubleMatrix & | mDr, | ||
const DoubleMatrix & | mLl, | ||
const DoubleMatrix & | mEr, | ||
double | m3sq, | ||
double | mH1sq, | ||
double | mH2sq, | ||
double | mGravitino | ||
) |
Constructor sets RPC SUSY parameters to s, gaugino masses to mG, trilinears to aU, aD, aE for au, ad, ae trilnears respectively, \(m_Q^2\)=mQl, \(m_U^2\)=mUr, \(m_D^2\)=mDr, \(m_L^2\)=mLl, \(m_E^2\)=mEr, \( m_3^2\)=m3sq, \(m_{H_1}^2\)=mH1sq, \(m_{H_2}^2\)=mH2sq, mu parameter, number of loops=l, and threshold parameter=t
void softsusy::MssmSoftPars::addAmsb | ( | const MssmSusy & | run, |
double | m32 | ||
) |
Give it a SUSY object and a value of M3/2, and it will return a soft object with AMSB soft breaking terms. Note that the sleptons will be tachyonic, ie nothing has been done to fix that problem. Note that in the following, we are neglecting all Yukawa couplings except that of the third family.
three family version – will be out by some two-loop terms. See hep-ph/9904378
void softsusy::MssmSoftPars::anomalousDeriv | ( | const MssmSusy & | xx, |
DoubleMatrix & | gEE, | ||
DoubleMatrix & | gLL, | ||
DoubleMatrix & | gQQ, | ||
DoubleMatrix & | gUU, | ||
DoubleMatrix & | gDD, | ||
double & | gH1H1, | ||
double & | gH2H2 | ||
) | const |
Returns derivatives of anomalous dimensions of fields with respect to renormalisation scale in MSSM for: RH leptons, LH leptons, LH quarks, RH up quarks, RH down quarks, H1 and H2 respectively
MssmSoftPars softsusy::MssmSoftPars::beta2 | ( | const MssmSusy & | xx, |
sBrevity & | a | ||
) | const |
Returns numerical beta functions of parameters and Brevity.
< dominant 3rd family approximation
< fully mixed case – slower
new dominant 3-family version
Old full 3-family version
The following are valid in the third-family approximation – they are much faster, and should be a good approximation to the no-mixed case
In the mixed case, we need to use the slower full 3-family expressions
void softsusy::MssmSoftPars::display | ( | DoubleVector & | v, |
int & | k | ||
) | const |
Obligatory function that outputs contents in v, k is incremented on output by the number of inputs.
double softsusy::MssmSoftPars::displayGaugino | ( | int | i | ) | const |
Return \(M_{G_i}\)
void softsusy::MssmSoftPars::minimalGmsb | ( | const MssmSusy & | xx, |
int | n5, | ||
double | LAMBDA, | ||
double | mMess, | ||
double | cgrav | ||
) |
Boundary conditions to be applied at messenger scale for Gauge mediated SUSY breaking (see hep-ph/9703211 for example), n5 is the number of 5-plets, mMess is the messenger scale and lambda is the GMSB scale
hep-ph/9801271
There is a relative minus in the mGMSB conditions for gaugino masses, since these equations are for L=-M/2 gaugino gaugino. See hep-ph/9801271: BCA 27/7/12
void softsusy::MssmSoftPars::setMh2Squared | ( | double | f | ) |
Sets \( m_{H_2}^2\)
void softsusy::MssmSoftPars::u1R_PQflip | ( | MssmSusy & | s | ) |
This flips the signs of various parameters (including mu) with no physical effect - apparently, through a U(1)_{R-PQ} transformation. See hep-ph/0312378
void softsusy::MssmSoftPars::universal | ( | MssmSusy & | s, |
double | m0, | ||
double | m12, | ||
double | a0, | ||
double | mu, | ||
double | m3sq | ||
) |
Reads in universal boundary conditions at the current scale: m0, M1/2, A0, B-parameter and mu