#include <softsusy.h>
Inheritance diagram for MssmSoftsusy:
Public Member Functions | |
MssmSoftsusy () | |
Default constructor fills object with zeroes. | |
MssmSoftsusy (const MssmSusy &) | |
Constructor sets SUSY parameters only from another object. | |
MssmSoftsusy (const MssmSoftsusy &) | |
Constructor copies another object. | |
MssmSoftsusy (const SoftParsMssm &s, const sPhysical &sp, double mu, int l, int t, double hv) | |
const MssmSoftsusy & | operator= (const MssmSoftsusy &s) |
Set all data in the object equal to another. | |
MssmSoftsusy | displayMssmSoft () const |
Displays whole object as a const. | |
sPhysical | displayPhys () const |
Displays physical parameters only. | |
drBarPars | displayDrBarPars () const |
sProblem | displayProblem () const |
Returns any problem flags associated with the object. | |
QedQcd | displayDataSet () const |
Gives the low energy Standard Model data set used for the object. | |
double | displayMinpot () const |
Returns minimum of Higgs potential. | |
double | displayMsusy () const |
Returns Higgs minimisation scale. | |
double | displayMw () const |
double | displayMwRun () const |
Returns DRbar MW, must be calculated by calcDrBarPars first. | |
double | displayMzRun () const |
Returns DRbar MZ, must be calculated by calcDrBarPars first. | |
double | displayTadpole1Ms () const |
displays t_1/v_1 tadpole | |
double | displayTadpole2Ms () const |
MssmSoftsusy | displaySoftsusy () const |
Returns object as a const. | |
double | displayMz () const |
Returns value of pole MZ being used. | |
void | flagIrqfp (bool a) |
Flags Infra-red quasi fixed point breach problem. | |
void | flagNonperturbative (bool a) |
Flags non-perturbative RG evolution. | |
void | flagTachyon (bool a) |
Flags a negative-mass squared scalar (really a CCB problem). | |
void | flagB (bool a) |
Flags problem with Higgs potential minimum. | |
void | flagBadConvergenve (bool a) |
Flags a really bad convergence: nowhere near a solution. | |
void | flagNoConvergence (bool a) |
Flags fact that calculation hasn't acheived required accuracy. | |
void | flagNoMuConvergence (bool a) |
Flags fact that mu sub iteration didn't converge. | |
void | flagNoRhoConvergence (bool a) |
Flags fact that rho parameter sub iteration didn't converge. | |
void | flagMusqwrongsign (bool a) |
Flags point inconsistent with electroweak symmetry breaking. | |
void | flagHiggsufb (bool a) |
Flags an inconsistent Higgs minimum. | |
void | flagAllProblems (bool a) |
Sets all problems equal to either true or false (contained in a). | |
void | setSoftsusy (const MssmSoftsusy &s) |
Sets whole object equal to another. | |
void | setData (const QedQcd &r) |
Sets low energy Standard Model fermion mass and gauge coupling data. | |
void | setMinpot (double) |
Sets potential value at minimum of Higgs potential. | |
void | setMsusy (double) |
Sets scale of Higgs potential minimisation and sparticle mass calculation. | |
void | setMw (double) |
sets pole MW prediction | |
void | setPhys (const sPhysical &s) |
Sets all physical parameters. | |
void | setDrBarPars (const drBarPars &s) |
Sets tree-level DRbar parameters. | |
void | doTadpoles (double mt, double sinthDRbar) |
Does the full 2-loop calculation of both tadpoles and sets them. | |
double | doCalcTadpole1oneLoop (double mt, double sinthDRbar) |
Does the calculation of one-loop pieces of ![]() | |
double | doCalcTadpole2oneLoop (double mt, double sinthDRbar) |
Does the calculation of one-loop pieces of ![]() | |
virtual void | calcTadpole1Ms1loop (double mt, double sinthDRbar) |
Calculates and sets the one-loop pieces of ![]() | |
virtual void | calcTadpole2Ms1loop (double mt, double sinthDRbar) |
Calculates then sets the one-loop pieces of ![]() | |
void | addStopCorrection (double p, DoubleMatrix &mass, double mt) |
void | addSdownCorrection (DoubleMatrix &mass, int family) |
void | addSbotCorrection (DoubleMatrix &mass, double mt) |
void | addSlepCorrection (DoubleMatrix &mass, int family) |
void | addStauCorrection (DoubleMatrix &mass, double mtau) |
void | addSupCorrection (DoubleMatrix &mass, int family) |
void | addSnuTauCorrection (double &mass) |
void | addSnuCorrection (double &mass, int family) |
void | addSquarkCorrection (DoubleMatrix &mass) |
void | doUpSquarks (double mt, double pizztMS, double sinthDRbarMS, int accuracy) |
void | doDownSquarks (double mb, double pizztMS, double sinthDRbarMS, int accuracy, double mt) |
void | doChargedSleptons (double mT, double pizztMS, double sinthDRbarMS, int accuracy) |
void | doSnu (double pizztMS, int accuracy=0) |
void | treeUpSquark (DoubleMatrix &mass, double mtrun, double pizztMS, double sinthDRbarMS, int family) |
void | treeDownSquark (DoubleMatrix &mass, double mbrun, double pizztMS, double sinthDRbarMS, int family) |
void | treeChargedSlepton (DoubleMatrix &mass, double mTrun, double pizztMS, double sinthDRbarMS, int family) |
void | treeSnu (double &mSq, double pizztMS, int family) |
void | calcDrBarPars () |
Organises tree-level calculation of all sparticle masses and mixings. | |
virtual void | sparticleThresholdCorrections (double tb) |
double | qedSusythresh (double alphaEm, double Q) const |
double | qcdSusythresh (double alphasMSbar, double Q) const |
double | calcMs () const |
virtual void | physical (int accuracy) |
double | calcRunningMt () |
double | calcRunningMtau () const |
double | calcRunningMb () const |
double | calcSinthdrbar () const |
Calculates DRbar sin theta_w at the current scale from gauge couplings. | |
double | getVev () const |
Calculates Higgs VEV parameter from gauge couplings and MZ. | |
double | getVev (double pizzt) const |
virtual void | charginos (int accuracy, double piwwt) |
virtual void | addCharginoLoop (DoubleMatrix &) |
Adds the loop corrections on to an input tree-level chargino mass. | |
virtual void | neutralinos (int accuracy, double piwwt, double pizzt) |
virtual void | addNeutralinoLoop (DoubleMatrix &) |
Adds the loop corrections on to an input tree-level neutralino mass. | |
virtual void | gluino (int accuracy) |
Calculates pole gluino mass to 1-loop SUSY corrections. | |
virtual void | higgs (int accuracy, double piwwt, double pizzt) |
virtual int | rewsbMu (int sgnMu, double &mu) const |
virtual int | rewsbM3sq (double, double &) const |
virtual void | rewsb (int sgnMu, double mt) |
virtual void | rewsbTreeLevel (int sgnMu) |
void | iterateMu (double &munew, int sgnMu, double mt, int maxTries, double pizztMS, double sinthDRbar, double tol, int &err) |
double | predTanb () const |
double | predMzsq (double &tanb) const |
DoubleVector | fineTune (void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &bcPars, double MX) |
double | it1par (int numPar, const DoubleVector &bcPars) |
double | ufb3sl (double) |
double | realMinMs () const |
virtual double | piZZT (double p, double Q, bool usePoleMt=false) const |
virtual double | piWWT (double p, double Q, bool usePoleMt=false) const |
virtual double | piHpHm (double p, double Q) const |
double | piZGT (double p, double Q) const |
virtual double | piAA (double p, double Q) const |
double | pis1s1 (double p, double q) const |
double | pis1s2 (double p, double q) const |
double | pis2s2 (double p, double q) const |
double | sinSqThetaEff () |
Calculates sin^2 theta^l_eff. | |
virtual void | rhohat (double &outrho, double &outsin, double alphaMZDRbar, double pizztMZ, double piwwt0, double piwwtMW, double tol, int maxTries) |
double | deltaVb (double outrho, double outsin, double alphaDRbar, double pizztMZ) const |
double | dRho (double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwtMW) |
double | dR (double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwt0) |
double | maxMass () const |
Returns the mass of the heaviest SUSY particle, excluding gravitino. | |
int | lsp (double &mass, int &posi, int &posj) const |
int | nlsp (double &mass, int &posi, int &posj) const |
void | printShort () const |
Prints a list of important sparticle/Higgs masses to standard output. | |
void | printLong () |
Prints a list of all sparticle/Higgs masses to standard output. | |
virtual void | printObj () |
Prints whols object to standard output. | |
double | thet (double a, double b, double c) |
log(max(a^2, b^2, c^2) / Q^2) | |
double | lowOrg (void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false) |
void | itLowsoft (int maxTries, double &mx, int sgnMu, double tol, double tanb, void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &pars, bool gaugeUnification, bool ewsbBCscale) |
virtual void | methodBoundaryCondition (const DoubleVector &pars) |
void | isajetNumbers764 (double &mtopPole, double &mGPole, double &smu, double &mA, double &tanb, double &mq1l, double &mdr, double &mur, double &meL, double &meR, double &mql3, double &mdr3, double &mur3, double &mtauL, double &mtauR, double &at, double &ab, double &atau, double &mq2l, double &msr, double &mcr, double &mmuL, double &mmuR, double &m1, double &m2) const |
void | isajetInterface764 (char fname[80]) const |
prints a file into fname which acts as an input to isajet | |
void | ssrunInterface764Inside (char fname[80], fstream &) const |
void | ssrunInterface764 (char fname[80], char softfname[80]) const |
void | isawigInterface764 (char fnamein[80], char fnameout[80], char fnamesoft[80]) const |
void | microMegasInterface (char fname[80]) const |
Outputs to softoutput a micromegas file (the name of which is input). | |
void | lesHouchesAccordOutput (char model[], const DoubleVector &pars, int sgnMu, double tanb, double qMax, int numPoints, double mb, double mtau, double mgut, bool altEwsb=false) |
Outputs with Les Houches accord conventions to standard output. | |
void | lesHouchesAccordOutput (char model[], const DoubleVector &pars, int sgnMu, double tanb, double qMax, char fileName[], int numPoints, double mb, double mtau, double mgut, bool altEwsb=false) |
Outputs with Les Houches accord conventions to a file "fileName". | |
void | microMegasIn (char fname[80]) |
Inputs a micromegas filename into a MssmSoftsusy object. | |
void | sCkmRotation (DoubleMatrix &deltaULL, DoubleMatrix &deltaURR, DoubleMatrix &deltaULR, DoubleMatrix &deltaDLL, DoubleMatrix &deltaDRR, DoubleMatrix &deltaDLR) const |
void | outputFcncs () const |
Prints out quark FCNC delta parameters (defined by Masiero and co). | |
void | doUfb3 (double mgut) |
|
Sets all parameters from s, sp, mu is the mu superpotential parameter, l is the number of loops used for RG evolution, t is the thresholds accuracy parameter, mg is the gravitino mass, hv is the Higgs VEV parameter. |
|
Adds one-loop corrections to sbottom mass matrix at p=root(mb1 mb2) IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds one-loop corrections to sbottom mass matrix at p=root(mb1 mb2) IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds one-loop corrections to sel_fam mass matrix at p=root(msel1 msel2) IO parameters: mass=tree level mass on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds one-loop corrections to sneutrino mass of family "family" IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added |
|
Adds one-loop corrections to tau sneutrino mass IO parameters: p=external momentum, mass=tree level mass matrix on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds approximate one-loop corrections to squark mass matrix for first two families. IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added |
|
Adds one-loop corrections to stau mass matrix at p=root(mtau1 mtau2) IO parameters: mass=tree level mass on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds one-loop corrections to stop mass matrix IO parameters: p=external momentum, mass=tree level mass matrix on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Adds one-loop corrections to stau mass matrix at p=root(mtau1 mtau2) IO parameters: mass=tree level mass on input, is returned with radiative corrections added, mt=DR bar top mass |
|
Calculates the best scale at which to do symmetry breaking: |
|
Applies approximate 1-loop SUSY corrections to mb(MZ) in order to return the DRbar running value |
|
Applies 1-loop SUSY and 2-loop QCD corrections to pole mt in order to return the DRbar running value at the current scale |
|
Applies approximate 1-loop SUSY corrections to pole mtau in order to return the DRbar running value at the current scale |
|
Calculates pole chargino masses and mixing using approximate 1-loop SUSY corrections. IO parameters: piwwt is the W self-energy at the current, accuracy is the number of loops required (0 or 1 currently) |
|
Calculates delta_v corrections for outrho=DRbar rho parameter, outsin=DRbar sin theta_w, alphaDRbar=alpha(Q) in the DR bar scheme, pizztMZ=self-energy of the Z at MZ |
|
Displays tree-level masses and mixings of sparticles and third generation fermions |
|
Returns predicted pole MW |
|
displays t_2/v_2 tadpole |
|
Organises calculation of all slepton masses. IO parameters: mT=DRbar tau mass, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w, accuracy=0,1: number of loops to add to tree-level squark mass matrix |
|
Organises calculation of all down squark masses. IO parameters: mb=DRbar bottom mass, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w, accuracy=0,1: number of loops to add to tree-level squark mass matrix |
|
Organises calculation of all sneutrino masses, pizztMS=Z self energy at Q=M_SUSY |
|
Sets the minimum of potential to be the difference between the UFB-3 direction minimum and the standard EW breaking minimum. mgut is obviously the high-scale at which boundary conditions are employed |
|
Organises calculation of all up squark masses. IO parameters: mt=DRbar top mass, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w, accuracy=0,1: number of loops to add to tree-level squark mass matrix |
|
Calculates delta r corrections for outrho=DRbar rho parameter, outsin=DRbar sin theta_w, alphaDRbar=alpha(Q) in the DR bar scheme, pizztMZ=self-energy of the Z at p=MZ, pizzt0=self-energy of the W at p=0 |
|
Calculates delta rho corrections for outrho=DRbar rho parameter, outsin=DRbar sin theta_w, alphaDRbar=alpha(Q) in the DR bar scheme, pizztMZ=self-energy of the Z at MZ, piwwtMW=self-energy of the W at p=MW |
|
Calculates fine-tuning for soft parameters and mu, m_3^2, top Yukawa. IO parameters: bcPars should be a vector giving the high-scale SUSY breaking boundary condition parameters, MX is the high-scale, boundaryCondition is the user-supplied function that sets the SUSY breaking BCs. |
|
Input for this one (saves time, possibly) is to give the self-energy of the Z at the current scale |
|
Calculates pole Higgs masses and mixings: full 1-loop SUSY corrections and 2-loop alpha_t (alpha_s + alpha_t) + alpha_s alpha_b effective potential corrections. IO parameters: piwwt is the W self-energy at M_SUSY, accuracy is number of loops (0 or 1) to use and pizzt is the Z self-energy at M_SUSY |
|
Works out how best to fit the isajet numbers to the spectrum. There are problems with the Higgs and sbottoms because ISAJET assumes certain tree-level relations between masses that are broken by SOFTSUSY's higher accuracy. The differences get large for high tan beta around 50, at around 10 they're typically only a percent. Output parameters: mtopPole is the pole mass of the top quark, mGPole is the gluino pole mass, smu is the superpotential Higgs parameter, mA is the pseudoscalar Higgs mass, tanb is tan beta, mq1l is the mass of the LH first family squark, mdr is the mass of the RH first family squark, meL is the LH selectron mass, meR is the RH selectron mass, mql3 is the LH stop mass, mdr3 is the RH sbottom mass, mur3 is the RH stop mass, mtauL is the LH stau mass, mtauR is the RH stau mass, at is the stop trilinear term, ab is the sbottom trilinear term, atau is the stau trilinear term, mq2l, msr, mcr are 2nd family squark masses, mmuL, mmuR are the smu masses, m1 and m2 are the gaugino mass parameters. |
|
Prints a file into fnamesoft that can be input into isawig. fnamein gives the ISASSRUN output file name, fnameout is the filename for the ISAWIG input file |
|
Give it a GUT scale object consistent with rewsb and it'll return the fine tuning of one parameter specified by numPar ht, mu and m3sq at the high |
|
Obtains solution of one-loop effective potential minimisation via iteration technique. Currently includes: all 1-loop SUSY tadpoles, plus 2-loop alpha_t (alpha_t + alpha_s) + alpha_b alpha_s corrections IO parameters: munew=current value of mu for this iteration, sgnMu=desired sign of mu, mt=DRbar mtop, maxTries=maximum number of iterations before it bails out, pizzMS=self-energy of MZ at current scale, sinthDRbar=DRbar value of sin theta_w, tol=desired fractional accuracy on mu, err=error flag: err=1 if no iteration reached, 2 if incorrect rewsb |
|
Main iteration routine: Boundary condition is the theoretical condition on parameters at the high energy scale mx: the parameters themselves are contained within the vector. IO parameters: maxTries is the maximum number of iterations allowed, mx is the GUT scale (negative if you require gauge unification), sgnMu is the desired sign of mu: + or - 1. If mu is 0, mu is set initially as a boundary condition. tanb = desired value of DR bar tan beta(MZ). |
|
Driver calculation to determine all sparticle masses and parameters. Returns low energy softsusy object consistent with BC's m0 etc at MGUT. oneset should be at MZ and contains the SM data to fit the model to. If the running comes into difficulty, eg if a Landau pole is reached, it returns a ZERO object: no result is possible! Boundary condition is the theoretical condition on parameters at the high energy scale mx: the parameters themselves are contained within the vector. IO parameters: sgnMu is the desired sign of mu: + or - 1. If mu is 0, mu is set initially as a boundary condition. tanb = desired value of DR bar tan beta(MZ). boundaryCondition is the function which sets to SUSY BCs at the high scale, mxGuess is the GUT scale gaugeUnification=true if the user requests true gauge unification, in which case mxGuess is the first guess. Returns actual GUT scale. ewsbBCscale = false (or omitted) means that the boundary condition on SUSY breaking is set in the usual way. If it is true, the boundary condition is set to |
|
Returns lsp mass in mass and function return labels which particle is lsp=0 means LSP is neutralino, 1=up squark, 2=down squark, 3=sleptons, 4=charginos, 5=sneutrinos, 6=gluino. posi, posj give the "handedness" (ie 1 or 2) for scalars and family respectively. |
|
Dummy function to allow users to re-define it in user supplied objects. Pars will contain necessary parameters to describe high-scale boundayr conditions on SUSY breaking terms |
|
Calculates pole neutralino masses and mixingusing approximate 1-loop SUSY corrections. IO parameters: piwwt is the W self-energy at M_SUSY, accuracy is the number of loops required (0 or 1 currently), pizzt is the Z self-energy at M_SUSY |
|
Returns nlsp mass in mass and function return labels which particle is nlsp=0 means NLSP is neutralino, 1=up squark, 2=down squark, 3=sleptons, 4=charginos, 5=sneutrinos, 6=gluino. posi, posj give the "handedness" (ie 1 or 2) for scalars and family respectively. |
|
Calculates physical sparticle masses to accuracy number of loops. Should be called at M_{SUSY}. |
|
Calculates transverse part of A^0 self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of H^+H^- self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of W self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates Z gamma self-energy: for p=external momentum, Q=renormalisation scale |
|
Calculates transverse part of Z self-energy: for p=external momentum, Q=renormalisation scale |
|
Predicts value of MZ(pole) from values of soft parameters and mu that we have. tanb=tan beta is also predicted |
|
This is a check: predicts tan beta from the values of soft parameters and mu that we have |
|
Does SUSY (and other) threshold corrections to alphaS Input alphas in MSbar and it returns it in DRbar scheme. From hep-ph/9606211 |
|
Does SUSY (and other) threshold corrections to alphaEm - returns alpha in DRbar scheme at scale Q. From hep-ph/9606211. Input empirical value of alpha at 0 momentum.... |
|
You should evaluate this at a scale MSusy average of stops. Returns depth of electroweak minimum |
|
Organises high accuracy rewsb: call it at the low scale M_{SUSY} IO parameters: sgnMu is +/-1 (desired sign of mu), mt is DRbar top mass Reimplemented in MssmSoftsusy2. |
|
Tree-level REWSB calculation, returning m3sq at correct value consistent with mu |
|
Tree-level REWSB calculation, returning mu at correct value. sgnMu is the required sign (+/- 1). Returns 1 if mu^2<0, indicating an inconsistent minimum |
|
Organises tree-level rewsb: call it at the low scale M_{SUSY} IO parameters: sgnMu is +/-1 (desired sign of mu) Reimplemented in MssmSoftsusy2. |
|
Iterative determination of rho parameter consistent with muon decay constant, MZ and alpha_0. IO parameters: outrho and outsin are the current DRbar values of sin theta_w and the rho parameter. alphaMZDRbar=already calculated value of fine-structure constant, pizztMZ=Z self-energy at Q=MZ, tol=desired piwwt0=self-energy of the W at p=0, piwwtMW=self-energy of the W at p=MW accuracy and maxTries is the maximum number of iterations the routines will allow |
|
all matrices in inputs should be 3x3: they're output as the quark FCNC matrices. |
|
For an input tan beta=tb, sets gauge and Yukawa couplings according to the tree-level spectrum and data set |
|
prints a file into fname which acts as an input to isajet First name input is the name of an OUTPUT file from ssrun, the second name is the name of the interface file for INPUT to ssrun |
|
prints a file into fname which acts as an input to isassrun: fstream should be opened before calling |
|
Returns tree-level down squark mass matrix in "mass". IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added, mTrun=DR bar tau mass, family=generation of slepton, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w |
|
Returns tree-level down squark mass matrix in "mass". IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added, mbrun=DR bar bottom mass, family=generation of squark, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w |
|
Organises calculation of all sneutrino masses, pizztMS=Z self energy at Q=M_SUSY, mSq=mass of sneutrino, family=generation of sneutrino |
|
Returns tree-level up squark mass matrix in "mass". IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added, mtrun=DR bar top mass, family=generation of squark, pizztMS=Z self energy at Q=M_SUSY, sinthDRbarMS=DRbar value of sin theta_w |
|
Input mx the scale up to which you search for minima. Returns minimum value of potential along UFB3 direction. Does ufbs truly properly but takes a long time. |