Main Page | Class Hierarchy | Class List | File List | Class Members | File Members

MssmSoftsusy Class Reference

#include <softsusy.h>

Inheritance diagram for MssmSoftsusy:

SoftParsMssm MssmSusy RGE MssmSoftsusy2 List of all members.

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 MssmSoftsusyoperator= (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 $ t_1 / v_1 $.
double doCalcTadpole2oneLoop (double mt, double sinthDRbar)
 Does the calculation of one-loop pieces of $ t_2 / v_2 $.
virtual void calcTadpole1Ms1loop (double mt, double sinthDRbar)
 Calculates and sets the one-loop pieces of $ t_1 / v_1 $.
virtual void calcTadpole2Ms1loop (double mt, double sinthDRbar)
 Calculates then sets the one-loop pieces of $ t_2 / v_2 $.
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)

Detailed Description

Contains all supersymmetric MSSM parameters, and is basically the SOFTSUSY kernel that incorporates everything


Constructor & Destructor Documentation

MssmSoftsusy::MssmSoftsusy const SoftParsMssm s,
const sPhysical sp,
double  mu,
int  l,
int  t,
double  hv
[inline]
 

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.


Member Function Documentation

void MssmSoftsusy::addSbotCorrection DoubleMatrix mass,
double  mt
 

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

void MssmSoftsusy::addSdownCorrection DoubleMatrix mass,
int  family
 

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

void MssmSoftsusy::addSlepCorrection DoubleMatrix mass,
int  family
 

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

void MssmSoftsusy::addSnuCorrection double &  mass,
int  family
 

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

void MssmSoftsusy::addSnuTauCorrection double &  mass  ) 
 

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

void MssmSoftsusy::addSquarkCorrection DoubleMatrix 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

void MssmSoftsusy::addStauCorrection DoubleMatrix mass,
double  mtau
 

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

void MssmSoftsusy::addStopCorrection double  p,
DoubleMatrix mass,
double  mt
 

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

void MssmSoftsusy::addSupCorrection DoubleMatrix mass,
int  family
 

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

double MssmSoftsusy::calcMs  )  const
 

Calculates the best scale at which to do symmetry breaking: $ M_{SUSY}=\sqrt{m_{{\tilde t}_1 {\tilde t}_2}} $. Should only be called after calcDrBarPars.

double MssmSoftsusy::calcRunningMb  )  const
 

Applies approximate 1-loop SUSY corrections to mb(MZ) in order to return the DRbar running value

double MssmSoftsusy::calcRunningMt  ) 
 

Applies 1-loop SUSY and 2-loop QCD corrections to pole mt in order to return the DRbar running value at the current scale

double MssmSoftsusy::calcRunningMtau  )  const
 

Applies approximate 1-loop SUSY corrections to pole mtau in order to return the DRbar running value at the current scale

void MssmSoftsusy::charginos int  accuracy,
double  piwwt
[virtual]
 

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)

double MssmSoftsusy::deltaVb double  outrho,
double  outsin,
double  alphaDRbar,
double  pizztMZ
const
 

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

drBarPars MssmSoftsusy::displayDrBarPars  )  const [inline]
 

Displays tree-level masses and mixings of sparticles and third generation fermions

double MssmSoftsusy::displayMw  )  const [inline]
 

Returns predicted pole MW

double MssmSoftsusy::displayTadpole2Ms  )  const [inline]
 

displays t_2/v_2 tadpole

void MssmSoftsusy::doChargedSleptons double  mT,
double  pizztMS,
double  sinthDRbarMS,
int  accuracy
 

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

void MssmSoftsusy::doDownSquarks double  mb,
double  pizztMS,
double  sinthDRbarMS,
int  accuracy,
double  mt
 

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

void MssmSoftsusy::doSnu double  pizztMS,
int  accuracy = 0
 

Organises calculation of all sneutrino masses, pizztMS=Z self energy at Q=M_SUSY

void MssmSoftsusy::doUfb3 double  mgut  )  [inline]
 

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

void MssmSoftsusy::doUpSquarks double  mt,
double  pizztMS,
double  sinthDRbarMS,
int  accuracy
 

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

double MssmSoftsusy::dR double  outrho,
double  outsin,
double  alphaDRbar,
double  pizztMZ,
double  piwwt0
 

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

double MssmSoftsusy::dRho double  outrho,
double  outsin,
double  alphaDRbar,
double  pizztMZ,
double  piwwtMW
 

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

DoubleVector MssmSoftsusy::fineTune void(*)(MssmSoftsusy &, const DoubleVector &)  boundaryCondition,
const DoubleVector bcPars,
double  MX
 

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.

double MssmSoftsusy::getVev double  pizzt  )  const
 

Input for this one (saves time, possibly) is to give the self-energy of the Z at the current scale

void MssmSoftsusy::higgs int  accuracy,
double  piwwt,
double  pizzt
[virtual]
 

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

void MssmSoftsusy::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
 

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.

void MssmSoftsusy::isawigInterface764 char  fnamein[80],
char  fnameout[80],
char  fnamesoft[80]
const
 

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

double MssmSoftsusy::it1par int  numPar,
const DoubleVector bcPars
 

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

void MssmSoftsusy::iterateMu double &  munew,
int  sgnMu,
double  mt,
int  maxTries,
double  pizztMS,
double  sinthDRbar,
double  tol,
int &  err
 

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

void MssmSoftsusy::itLowsoft int  maxTries,
double &  mx,
int  sgnMu,
double  tol,
double  tanb,
void(*)(MssmSoftsusy &, const DoubleVector &)  boundaryCondition,
const DoubleVector pars,
bool  gaugeUnification,
bool  ewsbBCscale
 

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).

double MssmSoftsusy::lowOrg void(*)(MssmSoftsusy &, const DoubleVector &)  boundaryCondition,
double  mxGuess,
const DoubleVector pars,
int  sgnMu,
double  tanb,
const QedQcd oneset,
bool  gaugeUnification,
bool  ewsbBCscale = false
 

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 $\sqrt{m_{{\tilde t}_1} m_{{\tilde t}_2}} $, ie like in the "pheno MSSM".

int MssmSoftsusy::lsp double &  mass,
int &  posi,
int &  posj
const
 

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.

void MssmSoftsusy::methodBoundaryCondition const DoubleVector pars  )  [virtual]
 

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

void MssmSoftsusy::neutralinos int  accuracy,
double  piwwt,
double  pizzt
[virtual]
 

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

int MssmSoftsusy::nlsp double &  mass,
int &  posi,
int &  posj
const
 

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.

void MssmSoftsusy::physical int  accuracy  )  [virtual]
 

Calculates physical sparticle masses to accuracy number of loops. Should be called at M_{SUSY}.

double MssmSoftsusy::piAA double  p,
double  Q
const [virtual]
 

Calculates transverse part of A^0 self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::piHpHm double  p,
double  Q
const [virtual]
 

Calculates transverse part of H^+H^- self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::pis1s1 double  p,
double  q
const
 

Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::pis1s2 double  p,
double  q
const
 

Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::pis2s2 double  p,
double  q
const
 

Calculates transverse part of Higgs self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::piWWT double  p,
double  Q,
bool  usePoleMt = false
const [virtual]
 

Calculates transverse part of W self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::piZGT double  p,
double  Q
const
 

Calculates Z gamma self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::piZZT double  p,
double  Q,
bool  usePoleMt = false
const [virtual]
 

Calculates transverse part of Z self-energy: for p=external momentum, Q=renormalisation scale

double MssmSoftsusy::predMzsq double &  tanb  )  const
 

Predicts value of MZ(pole) from values of soft parameters and mu that we have. tanb=tan beta is also predicted

double MssmSoftsusy::predTanb  )  const
 

This is a check: predicts tan beta from the values of soft parameters and mu that we have

double MssmSoftsusy::qcdSusythresh double  alphasMSbar,
double  Q
const
 

Does SUSY (and other) threshold corrections to alphaS Input alphas in MSbar and it returns it in DRbar scheme. From hep-ph/9606211

double MssmSoftsusy::qedSusythresh double  alphaEm,
double  Q
const
 

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....

double MssmSoftsusy::realMinMs  )  const
 

You should evaluate this at a scale MSusy average of stops. Returns depth of electroweak minimum

void MssmSoftsusy::rewsb int  sgnMu,
double  mt
[virtual]
 

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.

int MssmSoftsusy::rewsbM3sq double  ,
double & 
const [virtual]
 

Tree-level REWSB calculation, returning m3sq at correct value consistent with mu

int MssmSoftsusy::rewsbMu int  sgnMu,
double &  mu
const [virtual]
 

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

void MssmSoftsusy::rewsbTreeLevel int  sgnMu  )  [virtual]
 

Organises tree-level rewsb: call it at the low scale M_{SUSY} IO parameters: sgnMu is +/-1 (desired sign of mu)

Reimplemented in MssmSoftsusy2.

void MssmSoftsusy::rhohat double &  outrho,
double &  outsin,
double  alphaMZDRbar,
double  pizztMZ,
double  piwwt0,
double  piwwtMW,
double  tol,
int  maxTries
[virtual]
 

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

void MssmSoftsusy::sCkmRotation DoubleMatrix deltaULL,
DoubleMatrix deltaURR,
DoubleMatrix deltaULR,
DoubleMatrix deltaDLL,
DoubleMatrix deltaDRR,
DoubleMatrix deltaDLR
const
 

all matrices in inputs should be 3x3: they're output as the quark FCNC matrices.

void MssmSoftsusy::sparticleThresholdCorrections double  tb  )  [virtual]
 

For an input tan beta=tb, sets gauge and Yukawa couplings according to the tree-level spectrum and data set

void MssmSoftsusy::ssrunInterface764 char  fname[80],
char  softfname[80]
const
 

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

void MssmSoftsusy::ssrunInterface764Inside char  fname[80],
fstream & 
const
 

prints a file into fname which acts as an input to isassrun: fstream should be opened before calling

void MssmSoftsusy::treeChargedSlepton DoubleMatrix mass,
double  mTrun,
double  pizztMS,
double  sinthDRbarMS,
int  family
 

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

void MssmSoftsusy::treeDownSquark DoubleMatrix mass,
double  mbrun,
double  pizztMS,
double  sinthDRbarMS,
int  family
 

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

void MssmSoftsusy::treeSnu double &  mSq,
double  pizztMS,
int  family
 

Organises calculation of all sneutrino masses, pizztMS=Z self energy at Q=M_SUSY, mSq=mass of sneutrino, family=generation of sneutrino

void MssmSoftsusy::treeUpSquark DoubleMatrix mass,
double  mtrun,
double  pizztMS,
double  sinthDRbarMS,
int  family
 

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

double MssmSoftsusy::ufb3sl double   ) 
 

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.


The documentation for this class was generated from the following files:
Generated on Fri Aug 26 16:37:13 2005 for SOFTSUSY by  doxygen 1.3.9.1