Contains all supersymmetric NMSSM parameters, incorporating R_p NMSSM. More...

#include <nmssmsoftsusy.h>

Inheritance diagram for softsusy::NmssmSoftsusy:

Public Member Functions
	NmssmSoftsusy ()
	Default constructor fills object with zeroes.

	NmssmSoftsusy (const NmssmSusy &)
	Constructor sets SUSY parameters only from another object.

	NmssmSoftsusy (const NmssmSoftsusy &)
	Constructor copies another object.

	NmssmSoftsusy (const MssmSoftsusy &)
	Constructor copies an MSSM object, needed for test vs mssm.

	NmssmSoftsusy (const NmssmSusyPars &ss, const SoftParsNmssm &s, const MssmSoftsusy &sp)

const NmssmSoftsusy &	operator= (const NmssmSoftsusy &s)
	Set all data in the object equal to another.

double	displayTadpoleSMs () const
	displays t_s/s tadpole

double	displayTadpoleSMs1loop () const
	displays t_2/v_s tadpole @1 loop

double	displaySoftAkappa () const

double	displaySoftAlambda () const

NmssmSusy	displayNmssmSusy () const

MssmSoftsusy	convertToMssm () const

bool	displayGUTlambda () const

bool	displayGUTkappa () const

bool	displayGUTmuPrime () const

bool	displayGUTxiF () const

bool	displayGUTxiS () const

bool	displayGUTsVev () const

bool	displayMICROMEGAS () const

bool	displayNMSDECAY () const

double	displayDrBarHiggsAccuracy () const
	DH: returns the achieved accuracy in Higgs mass diagonalisation.

double	displayPhysHiggsAccuracy () const

void	P1SfSfCouplings (DoubleMatrix &lTP1Lr, DoubleMatrix &lBP1Lr, DoubleMatrix &lTauP1Lr) const

void	P2SfSfCouplings (DoubleMatrix &lTP2Lr, DoubleMatrix &lBP2Lr, DoubleMatrix &lTauP2Lr) const
	PA: obtains NMSSM P2-sfermion-sfermion couplings.

void	P3SfSfCouplings (DoubleMatrix &lTP3Lr, DoubleMatrix &lBP3Lr, DoubleMatrix &lTauP3Lr) const
	PA: obtains NMSSM P3-sfermion-sfermion couplings.

void	H1SfSfCouplings (DoubleMatrix &lTS1Lr, DoubleMatrix &lBS1Lr, DoubleMatrix &lTauS1Lr, double gmzOcthW, double mu, double cosb, double v1) const
	PA: obtains NMSSM H1-sfermion-sfermion couplings.

void	H2SfSfCouplings (DoubleMatrix &lTS1Lr, DoubleMatrix &lBS1Lr, DoubleMatrix &lTauS1Lr, double gmzOcthW, double mu, double sinb) const
	PA: obtains NMSSM H2-sfermion-sfermion couplings.

void	SSfSfCouplings (DoubleMatrix &lTS3Lr, DoubleMatrix &lBS3Lr, DoubleMatrix &lTauS3Lr, double lam) const
	PA: obtains NMSSM S-sfermion-sfermion couplings.

double	doCalcTadSSfermions (DoubleMatrix lTS3Lr, DoubleMatrix lBS3Lr, DoubleMatrix lTauS3Lr, double q, double s) const
	PA:routine to calculate sfermiom contributions to (16 \pi^2) ts / s.

void	assignHiggs (DoubleVector &higgsm, DoubleVector &higgsa, DoubleVector &higgsc) const
	PA: for loop corrections, helps adding Higgs corrections in a tidy way.

double	doCalcTad1Higgs (double q, double costhDRbar, double g, double tanb) const
	PA: NMSSM routine to obtain Higgs loop parts of (16 \pi^2) t1/v1.

double	doCalcTad2Higgs (double q, double costhDRbar, double g, double tanb) const
	PA: NMSSM routine to obtain Higgs loop parts of (16 \pi^2) t2/v2.

double	doCalcTadSHiggs (double q, double tanb) const
	PA: NMSSM routine to obtain Higgs loop parts of (16 \pi^2) ts/s.

double	doCalcTad1Neutralinos (double q, double costhDRbar, double g, double tanb) const
	PA: NMSSM routine to obtain Neutralino loop parts of (16 \pi^2) t1/v1.

double	doCalcTad2Neutralinos (double q, double costhDRbar, double g, double sinb) const
	PA: NMSSM routine to obtain Neutralino loop parts of (16 \pi^2) t2/v2.

double	doCalcTadSNeutralinos (double q, double lam, double kap) const
	PA: NMSSM routine to obtain Neutralino loop parts of (16 \pi^2) tS/s.

double	doCalcTadSCharginos (double q, double lam) const
	PA: NMSSM routine to obtain Chargino loop parts of (16 \pi^2) tS/s.

virtual double	doCalcTadpole1oneLoop (double mt, double sinthDRbar) const
	Does the calculation of one-loop pieces of \( t_1 / v_1 \).

virtual double	doCalcTadpole2oneLoop (double mt, double sinthDRbar) const
	Does the calculation of one-loop pieces of \( t_2 / v_2 \).

virtual double	doCalcTadpoleSoneLoop (double mt, double sinthDRbar) const
	Does the calculation of one-loop pieces of \( t_s / s \).

virtual void	calcTadpole1Ms1loop (double mt, double sinthDRbar)

virtual void	calcTadpole2Ms1loop (double mt, double sinthDRbar)

virtual void	calcTadpoleSMs1loop (double mt, double sinthDRbar)

void	doTadpoles (double mt, double sinthDRbar)

virtual void	calcDrBarPars ()
	Organises tree-level calculation of all sparticle masses and mixings.

DoubleMatrix	addStopHiggs (double p, double mt, DoubleMatrix &higgs)

virtual void	addStopCorrection (double p, DoubleMatrix &mass, double mt)

DoubleMatrix	addSupHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

virtual void	addSupCorrection (DoubleMatrix &mass, int family)

DoubleMatrix	addSbotHiggs (double p, double mt, DoubleMatrix &higgs)

virtual void	addSbotCorrection (double p, DoubleMatrix &mass, double mt)

DoubleMatrix	addSdownHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

virtual void	addSdownCorrection (DoubleMatrix &mass, int family)

double	addSnuTauHiggs (double p, double &higgs)

virtual void	addSnuTauCorrection (double &mass)

double	addSnuHiggs (double p, int family, double &higgs)

virtual void	addSnuCorrection (double &mass, int family)

DoubleMatrix	addStauHiggs (double p, double mtau, DoubleMatrix &higgs)

virtual void	addStauCorrection (double p, DoubleMatrix &mass, double mtau)

DoubleMatrix	addSlepHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

virtual void	addSlepCorrection (DoubleMatrix &mass, int family)

void	getNeutPassarinoVeltman (double p, double q, DoubleMatrix &b0fn, DoubleMatrix &b1fn)

void	addNeutLoopHiggs (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)

virtual void	addNeutralinoLoop (double p, DoubleMatrix &)
	Adds the loop corrections on to an input tree-level neutralino mass.

virtual void	treeUpSquark (DoubleMatrix &mass, double mtrun, double pizztMS, double sinthDRbarMS, int family)

virtual void	treeDownSquark (DoubleMatrix &mass, double mbrun, double pizztMS, double sinthDRbarMS, int family)

virtual void	treeChargedSlepton (DoubleMatrix &mass, double mTrun, double pizztMS, double sinthDRbarMS, int family)

void	treeCharginos (DoubleMatrix &mass, double beta, double mw)
	LCT: new routine to allocate NMSSM chargino masses,.

void	treeNeutralinos (DoubleMatrix &mass, double beta, double mz, double mw, double sinthDRbar)
	LCT: new routine for NMSSM neutralino masses,.

void	calcDrBarGauginos (double beta, double mw, double mz, double sinth, drBarPars &eg)

void	treeHiggsAlt (DoubleMatrix &mS, DoubleMatrix &mPpr, DoubleMatrix &mP2, double &mHpm, double beta) const
	PA: fills tree level CP even and CP odd Higgs mass matrices.

void	treeHiggs (DoubleMatrix &mS, DoubleMatrix &mPpr, DoubleMatrix &mP2, double &mHpm, double beta) const
	PA: fills tree level CP even and CP odd Higgs mass matrices.

void	calcDrBarHiggs (double beta, double mz2, double mw2, double sinthDRbar, drBarPars &eg)

DoubleVector	beta () const
	Beta functions of NMSSM.

NmssmSoftsusy	beta2 () const
	Beta functions of NMSSM.

const DoubleVector	display () const
	Return contents of object in a vector: for RG evolution.

void	set (const DoubleVector &y)

void	setNmssmSusy (const NmssmSusy &n)

void	setGUTlambda (bool i)

void	setGUTkappa (bool i)

void	setGUTmuPrime (bool i)

void	setGUTxiF (bool i)

void	setGUTxiS (bool i)

void	setGUTsVev (bool i)

void	setMICROMEGAS (int i)

void	setNMSDECAY (int i)

double	VhAtMin (double v1, double v2, double s, double mHu2, double mHd2, double mSsq)

virtual int	rewsbMu (int sgnMu, double &mu) const

virtual int	rewsbSvev (int sgnMu, double &svev) const

virtual int	rewsbM3sq (double mu, double &m3sq) const

virtual int	rewsbKap (double &kap) const
	PA:: In case of Z3 invariance EWSB outputs kappa instead.

virtual int	rewsbXiS (double &xiS) const
	PA: third EWSB condition (for the singlet Higgs field)

virtual int	rewsbSvevNoZ3 (double &xiS) const

virtual int	rewsbmSsq (double &mSsq) const

virtual int	rewsbmH1sq (double &mH1sq) const
	PA: for non universal Higgs.

virtual int	rewsbmH2sq (double &mH211sq) const
	PA: for non universal Higgs.

virtual void	rewsbTreeLevel (int sgnMu)
	PA: Imposes EWSB at the tree level.

void	iterateMu (double &munew, int sgnMu, double mt, int maxTries, double pizztMS, double sinthDRbar, double tol, int &err)
	PA: finds mu iteratively in the case where we use EWSB to swap. More...

void	iterateSvev (double &sold, int sgnMu, double mt, int maxTries, double pizzMS, double sinthDRbar, double tol, int &err)

void	rewsb (int sgnMu, double mt, double muOld=-6.66e66, double eps=0.)
	ie (mu, m3sq, XiS) --> (mZ, tb, s) More...

virtual void	physical (int accuracy)

virtual void	addChaLoopHiggs (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const

virtual void	addCharginoLoop (double p, DoubleMatrix &)
	Adds the loop corrections on to an input tree-level chargino mass.

virtual void	charginos (int accuracy, double piwwt, sPhysical &phys)

virtual void	neutralinos (int accuracy, double piwwt, double pizzt, sPhysical &phys)
	Calculates pole neutralino masses and mixing. More...

bool	higgs (int accuracy, double piwwt, double pizzt, sPhysical &phys)

void	DegrassiSlavicMix (DoubleMatrix &P) const

virtual double	piZZTHiggs (double p, double q, double thetaWDRbar) const
	Calculates Higgs contribution to the transverse part of Z self-energy:

virtual double	piZZTNeutralinos (double p, double Q, double thetaWDRbar) const
	Calculates neutralino contrib. to the transverse part of Z self-energy:

virtual double	piZZT (double p, double Q, bool usePoleMt=false) const

virtual double	piWWTHiggs (double p, double q, double thetaWDRbar) const
	Calculates Higgs contribution to the transverse part of W self-energy:

virtual double	piWWTgauginos (double p, double Q, double thetaWDRbar) const
	Calculates gaugino contr. to transverse part of W self-energy:

virtual double	piWWT (double p, double Q, bool usePoleMt=false) const

double	pip1p1 (double p, double q) const
	PA: self energy routines for pseudo scalar self energies.

double	pip1p2 (double p, double q) const

double	pip2p2 (double p, double q) const

double	pip1p3 (double p, double q) const

double	pip2p3 (double p, double q) const

double	pip3p3 (double p, double q) const

void	getP1HiggsTriCoup (DoubleMatrix &spp1, DoubleMatrix &hphpp1, double cw2DRbar) const
	PA: Obtains trilnear couplings of P1-Higgs, for use in loop functions.

void	getP2HiggsTriCoup (DoubleMatrix &spp2, DoubleMatrix &hphpp2, double cw2DRbar) const
	PA: Obtains trilnear couplings of P2-Higgs for use in loop functions.

void	getP3HiggsTriCoup (DoubleMatrix &spp2, DoubleMatrix &hphpp2) const
	PA: Obtains trilnear couplings of P3-Higgs for use in loop functions.

void	getS1HiggsTriCoup (DoubleMatrix &sss1, DoubleMatrix &pps1, DoubleMatrix &hphps1, double thetaWDRbar) const
	PA: Obtains trilnear couplings of s1-higgs-higgs for use in loop functions.

void	getS2HiggsTriCoup (DoubleMatrix &sss2, DoubleMatrix &pps2, DoubleMatrix &hphps2, double thetaWDRbar) const
	PA: Obtains trilnear couplings of s2-higgs-higgs for use in loop functions.

void	getS3HiggsTriCoup (DoubleMatrix &sss3, DoubleMatrix &pps3, DoubleMatrix &hphps3) const
	PA: Obtains trilnear couplings of s3-higgs-higgs for use in loop functions.

double	pis1s1Higgs (double p, double q) const
	PA:Calculates (16 Pi^2) times the Higgs contribution to Higgs self-energy:

double	pis1s2Higgs (double p, double q) const

double	pis2s2Higgs (double p, double q) const

double	pis1s3Higgs (double p, double q) const

double	pis2s3Higgs (double p, double q) const

double	pis3s3Higgs (double p, double q) const

void	getS1NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const
	PA: obtains CP even Higgs-Neutralino couplings.

void	getS2NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const

void	getS3NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const

void	getP1NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const
	PA: obtains CP odd Higgs-Neutralino couplings.

void	getP2NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const

void	getP3NeutralinoCoup (ComplexMatrix &aChi, ComplexMatrix &bChi) const

double	pis1s1Neutralinos (double p, double q) const
	PA:Calculates (16 Pi^2) times Neutralino contrib. to Higgs self-energy:

double	pis1s2Neutralinos (double p, double q) const

double	pis2s2Neutralinos (double p, double q) const

double	pis1s3Neutralinos (double p, double q) const

double	pis2s3Neutralinos (double p, double q) const

double	pis3s3Neutralinos (double p, double q) const

double	pis1s3Charginos (double p, double q) const
	PA:Calculates (16 Pi^2) times Chargino contrib. to Higgs self-energy:

double	pis2s3Charginos (double p, double q) const

double	pis3s3Charginos (double p, double q) const

double	pis1s3Sfermions (double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau, DoubleMatrix ls3tt, DoubleMatrix ls3bb, DoubleMatrix ls3tautau) const
	PA:Calculates (16 Pi^2) times sfermion contrib. to Higgs self-energy:

double	pis2s3Sfermions (double p, double q, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau, DoubleMatrix ls3tt, DoubleMatrix ls3bb, DoubleMatrix ls3tautau) const

double	pis3s3Sfermions (double p, double q, DoubleMatrix ls3tt, DoubleMatrix ls3bb, DoubleMatrix ls3tautau) const

double	pis1s1 (double p, double q) const

double	pis1s2 (double p, double q) const

double	pis2s2 (double p, double q) const

double	pis1s3 (double p, double q) const

double	pis2s3 (double p, double q) const

double	pis3s3 (double p, double q) const

virtual double	h1s2Mix ()
	PA: gets h1 mixing element with Hu.

void	getHp1HiggsTriCoup (DoubleMatrix &ahphp1, DoubleMatrix &hhphp1) const
	LCT: Charged Higgs trilinear couplings with CP-odd and CP-even Higgs'.

void	getHp2HiggsTriCoup (DoubleMatrix &ahphp2, DoubleMatrix &hhphp2) const

void	getNeutralinoCharginoHpmCoup (ComplexMatrix &apph1, ComplexMatrix &apph2, ComplexMatrix &bpph1, ComplexMatrix &bpph2) const
	LCT: Charged Higgs trilinear couplings with neutralinos/charginos.

void	getGaugeHiggsHpmCoup (DoubleVector &wmhhp1, DoubleVector &wmahp1, DoubleVector &gamhphp1, DoubleVector &zhphp1, DoubleVector &wmhhp2, DoubleVector &wmahp2, DoubleVector &gamhphp2, DoubleVector &zhphp2) const
	LCT: Returns gauge-Higgs-hpm trilinear couplings with hpm in weak basis.

double	piHpm11Higgs (double p, double q) const

double	piHpm12Higgs (double p, double q) const

double	piHpm22Higgs (double p, double q) const

double	piHpm11Gaugino (double p, double q) const

double	piHpm12Gaugino (double p, double q) const

double	piHpm22Gaugino (double p, double q) const

double	piHpHm (double p, double q) const

virtual double	calcRunMtHiggs () const
	PA: routine to calculate the Higgs loop contributions to running Mt.

virtual double	calcRunMtNeutralinos () const
	PA: routine to calculate the Neutralino loop contributions to running Mt.

virtual double	calcRunMbHiggs () const
	PA: routine to calculate the Higgs loop contributions to running Mb.

virtual double	calcRunMbNeutralinos () const
	PA: routine to calculate the Neutralino loop contributions to running Mb.

virtual double	calcRunMtauHiggs () const
	PA: routine to calculate the Higgs loop contributions to running Mtau.

virtual double	calcRunMtauNeutralinos (double mTauSMMZ) const

NmssmSusyRGE	guessAtSusyMt (double tanb, const DoubleVector &nmpars, const QedQcd &oneset)

double	predTanb (double muSusy=-6.66e66) const

double	predMzsq (double &tanb, double muOld=-6.66e66, double eps=0.)

double	ewsbCondition1TreeLevel () const
	DH: returns the values of the tree-level EWSB conditions.

double	ewsbCondition2TreeLevel () const

double	ewsbConditionSTreeLevel () const

void	ewsbConditions (DoubleVector &values) const

void	predVevs (DoubleVector &vevs, int &err)

DoubleMatrix	fineTune (void(*boundaryCondition)(NmssmSoftsusy &, const DoubleVector &), const DoubleVector &bcPars, double MX, bool doTop=false)

double	calcBayesianNaturalness () const
	Calculates Bayesian naturalness a la arXiv:1709.07895.

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

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

void	printall () const

virtual void	printObj ()
	Prints whole object to standard output.

const NmssmSoftsusy &	displayNmssmSoftsusy () const
	Displays whole object as a const.

virtual void	lesHouchesAccordOutput (ostream &out, const char model[], const DoubleVector &pars, int sgnMu, double tanb, double qMax, int numPoints, bool ewsbBCscale)

void	headerSLHA (ostream &out)
	Proper citation info in SLHA file.

void	modselSLHA (ostream &out, const char model[], double qMax)
	This does the MODSEL block of SLHA.

virtual void	extparSLHA (ostream &out, const DoubleVector &pars, bool ewsbBCscale)
	EXTPAR block of SLHA.

virtual void	higgsMSLHA (ostream &out)
	higgs part of mass block of SLHA

virtual void	neutralinoCharginoMSLHA (ostream &out)
	neutralino and charigno part of mass block of SLHA

virtual void	neutralinoMixingSLHA (ostream &out)
	nmnmix block output of SLHA

void	nmhmixSLHA (ostream &)
	CP even mixing matrix NMHMIX block of SLHA.

void	nmamixSLHA (ostream &)
	CP odd mixing matrix NMHMIX block of SLHA.

void	nmssmrunSLHA (ostream &, const char *blockName="NMSSMRUN")
	NMSSMRUN block of SLHA.

void	extranmssmtoolsSLHA (ostream &)
	extra DRbar parameters block for NMSSMTools decays

void	yukawaMatricesSLHA (ostream &, const char *blockName="")
	extra DRbar Yukawa couplings for NMSSMTools decays

void	extramsoftSLHA (ostream &, const char *blockName="")
	extra DRbar soft parameters for NMSSMTools decays

void	extrahmixSLHA (ostream &, const char *blockName="")
	extra DRbar Higgs mixing parameters for NMSSMTools decays

void	extragaugeSLHA (ostream &, const char *blockName="GAUGE")
	extra DRbar gauge couplings for NMSSMTools decays

void	extrasfermionmixSLHA (ostream &out)
	extra DRbar smuon mixing matrix elements for NMSSMTools decays

virtual void	drbarSLHA (ostream &, int numPoints, double qMax, int n)
	outputs DRbar parameters at scale Q (default is at MSUSY)

virtual void	softsusySLHA (ostream &out)
	SOFTSUSY comments in SLHA.

virtual void	spinfoSLHA (ostream &out)
	This does the SPINFO block of SLHA.

Public Member Functions inherited from softsusy::NmssmSusyPars
	NmssmSusyPars ()
	Constructor fills object with zeroes by default. More...

	NmssmSusyPars (const NmssmSusyPars &)

	NmssmSusyPars (double sv, double lambda, double kappa, double xiF, double mupr)

const NmssmSusyPars &	operator= (const NmssmSusyPars &s)

void	setSvev (double s)
	sets DRbar running singlet vev.

void	setLambda (double l)
	sets the \lambda S H_u H_d coupling

void	setKappa (double k)
	sets the \frac{1}{3} \kappa S^3 coupling

void	setMupr (double m)
	sets the \frac{1}{2} mupr S^2 coupling

void	setXiF (double x)
	sets the xiF S coupling

void	setNmssmSusyPars (const NmssmSusyPars &s)
	sets the whole object

void	setNmssmSusyPars (const DoubleVector &y)

void	set (const DoubleVector &y, int &k)

double	displaySvev () const
	returns DRbar running Singlet Higgs vev

double	displayLambda () const
	returns superpotential parameter lambda

double	displayKappa () const
	returns superpotential parameter lambda

double	displayMupr () const
	returns mupr superpotential parameter

double	displayXiF () const
	returns xiF superpotential parameter

const NmssmSusyPars &	displayNmssmSusyPars () const
	returns whole object

void	display (DoubleVector &y, int &k) const
	whole object in a vector y, starting at index k

void	getOneLpAnom (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gDD, DoubleMatrix &gUU, double &gH1H1, double &gH2H2, double &gSS, nmsBrevity &a, int loops) const

void	getTwoLpAnom (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gDD, DoubleMatrix &gUU, double &gH1H1, double &gH2H2, double &gSS, nmsBrevity &a, int loops) const

void	anomalousDimension (const MssmSusy &s, DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, DoubleVector &dg, double &gH1H1, double &gH2H2, double &gSS, nmsBrevity &a, int loops) const

Public Member Functions inherited from softsusy::SoftParsNmssm
	SoftParsNmssm ()
	Default constructor fills object with zeroes.

	SoftParsNmssm (const SoftParsNmssm &s)

	SoftParsNmssm (double al, double ak, double ms, double msp, double x)

const SoftParsNmssm &	operator= (const SoftParsNmssm &s)

const SoftParsNmssm &	displaySoftParsNmssm () const

double	displayTrialambda () const

double	displayTriakappa () const

double	displayMsSquared () const
	Return \(m_{S}^2\)=mSsq.

double	displayMspSquared () const
	Return \(m_{S'}^2\)=mSpsq.

double	displayXiS () const
	Return xiS i.e \(xiS S\).

double	displaySoftAlambda (double lam) const

double	displaySoftAkappa (double ka) const

void	display (DoubleVector &y, int &k) const
	Whole object output in a doublevector. More...

bool	displayZ3 () const

void	setZ3 (bool z)

void	setTrialambda (double al)
	PA: Set trilinear SUSY breaking parameter alambda.

void	setTriakappa (double ak)
	PA: Set trilinear SUSY breaking parameter akappa.

void	setMsSquared (double f)
	PA: Sets soft scalar mass \(m_S^2 \|S\|^2\).

void	setMspSquared (double f)
	PA: Sets soft scalar mass \(m_S^2 S^2\).

void	setXiS (double f)
	PA: Sets soft breaking \( xiS S\).

void	setSoftParsNmssm (SoftParsNmssm const &s)

void	set (const DoubleVector &, int &k)

void	anomalousDeriv (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, double &gH1H1, double &gH2H2, double &gSS) const

void	addAmsb (double maux, const NmssmSusy &Nms, MssmSoftPars &m)
	adds NMSSM amsb contribution to soft masses. More...

void	standardSugra (double m0, double m12, double a0, const NmssmSusyPars &n, const MssmSusy &m, const MssmSoftPars &msoft)

void	standardsemiSugra (double m0, double m12, double a0, double Al, double Ak, const NmssmSusy &n, const MssmSoftPars &m, double mS=0.0)
	works for Z3 violating NMSSM. More...

void	universalScalars (double m0)
	Sets all flavour-diagonal SUSY breaking scalar masses to m0.

void	semiuniversalScalars (double m0)
	PA: for semi constrained models. More...

void	universalTrilinears (double a0, const NmssmSusy &n)
	Sets all NMSSM SUSY breaking trilinear couplings to a0. More...

void	semiuniversalTrilinears (double a0, double al, double ak, const NmssmSusy &n)

void	inputSoftParsOnly ()
	Reads in soft SUSY breaking parameters from a file.

Public Member Functions inherited from softsusy::MssmSoftsusy
	MssmSoftsusy ()
	Default constructor fills object with zeroes.

	MssmSoftsusy (const MssmSusyRGE &)
	Constructor sets SUSY parameters only from another object.

	MssmSoftsusy (const MssmSoftsusy &)
	Constructor copies another object.

	MssmSoftsusy (const MssmSusy &, const MssmSoftPars &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.

const MssmSoftsusy &	displayMssmSoft () const
	Displays whole object as a const.

const sPhysical &	displayPhys () const
	Displays physical parameters only.

sPhysical	displayPhysUncertainty (TMSSMBoundaryCondition boundaryCondition, double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false, int contributions=DeltaAll) const
	Displays uncertainty estimate of pole masses. More...

sPhysical	displayPhysUncertaintyScaleVariation () const
	Displays uncertainty estimate of pole masses by scale variation only. More...

const drBarPars &	displayDrBarPars () const

const sProblem &	displayProblem () const
	Returns any problem flags associated with the object.

const QedQcd &	displayDataSet () const
	Gives the low energy Standard Model data set used for the object.

double	displaySoftA (trilinears, int, int) const
	Return a trilinear element in "SUGRA style".

double	displayFracDiff () const
	Displays iteration accuracy attained.

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
	displays t_2/v_2 tadpole

double	displayTadpole1Ms1loop () const
	displays t_1/v_1 tadpole@1 loop

double	displayTadpole2Ms1loop () const
	displays t_2/v_2 tadpole@1 loop

double	displayMxBC () const

const MssmSoftsusy &	displaySoftsusy () const
	Returns object as a const. More...

double	displayMz () const
	Returns value of pole MZ being used.

bool	displaySetTbAtMX () const
	Is tan beta set at the user defined SUSY breaking scale?

bool	displayAltEwsb () const

bool	displayAltMt () const

bool	displayAltAlphaS () const

bool	displayAltAlphaEm () const

double	displayPredMzSq () const

double	displayQewsb () const

double	displayMatchingScale () const

void	flagMgutOutOfBounds (bool a)
	Flags weird mgut-type problems.

void	flagIrqfp (bool a)
	Flags Infra-red quasi fixed point breach problem.

void	flagNonperturbative (bool a)
	Flags non-perturbative RG evolution.

void	flagTachyon (tachyonType a)
	Flags a negative-mass squared scalar (really a CCB problem) More...

void	flagTachyonWarning (tachyonType a)

void	flagM3sq (bool a)
	Flags problem with Higgs potential minimum.

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	flagInaccurateHiggsMass (bool a)

void	flagHiggsufb (bool a)
	Flags an inconsistent Higgs minimum.

void	flagNotGlobalMin (bool a)
	LCT: Flags problem if not in global minimum of Higgs potential.

void	flagAllProblems (bool a, tachyonType b)
	Sets all problems equal to either true or false (contained in a)

void	flagProblemThrown (bool a)
	Flags a numerical exception eg number too big/small.

void	setProblem (const sProblem a)
	Sets all problem flags.

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)

void	setMw (double)
	sets pole MW prediction

void	setPhys (const sPhysical &s)
	Sets all physical parameters.

void	setMxBC (double mx)
	Sets the scale at which high-scale boundary conditions are applied.

void	setDrBarPars (const drBarPars &s)
	Sets tree-level DRbar parameters.

void	setSetTbAtMX (bool a)
	Sets the setTbAtMX flag.

void	useAlternativeEwsb ()
	Use alernative EWSB conditions: set mu and MA(pole)

void	useAlternativeMt ()
	Use alernative mt(MSSM,DR-bar) calculation.

void	useAlternativeAlphaS ()
	Use alernative alpha_s(MSSM,DR-bar) calculation.

void	useAlternativeAlphaEm ()
	Use alernative alpha_em(MSSM,DR-bar) calculation.

void	setPredMzSq (double a)
	Set MZ^2 predicted after iteration.

void	set (const DoubleVector &)
	Sets total set of RGE parameters equal to elements of a vector.

void	setQewsb (double q)
	Sets user-set scale qewsb.

void	setTwoLoopAlphasThresholds (bool sw)

void	setTwoLoopMtThresholds (bool sw)
	Switch 2-loop threshold \( O(\alpha_s^2) \) corrections to top mass.

void	setTwoLoopMbMtauThresholds (bool sw)

void	setAllTwoLoopThresholds (bool sw)
	Switch all two-loop thresholds to gauge-Yukawa couplings.

DoubleVector	beta () const
	Returns double vector containing numerical beta functions of parameters.

void	setFracDiff (double fD)
	sets fracDiff, needed for access by NmssmSoftsusy methods

void	H1SfSfCouplings (DoubleMatrix &lTS1Lr, DoubleMatrix &lBS1Lr, DoubleMatrix &lTauS1Lr, double gmzOcthW, double mu, double cosb, double v1) const
	fixes trilnear H1-sfermion-sfermion couplings

void	H2SfSfCouplings (DoubleMatrix &lTS2Lr, DoubleMatrix &lBS2Lr, DoubleMatrix &lTauS2Lr, double gmzOcthW, double mu, double sinb) const
	fixes trilnear H2-sfermion-sfermion couplings

double	doCalcTad1Sfermions (DoubleMatrix lTS1Lr, DoubleMatrix lBS1Lr, DoubleMatrix lTauS1Lr, double costhDRbar) const
	routine to calculate sfermiom contributions to H1 tadpole / v1 More...

double	doCalcTad2Sfermions (DoubleMatrix lTS2Lr, DoubleMatrix lBS2Lr, DoubleMatrix lTauS2Lr, double costhDRbar) const
	routine to calculate sfermiom contributions to (16 \pi^2) t1 / v1 More...

double	doCalcTad1fermions (double q, double v1) const
	fixes trilnear H1-fermion-fermion couplings More...

double	doCalcTad2fermions (double q) const
	fixes trilnear H2-fermion-fermion couplings More...

double	doCalcTad1Higgs (double q, double costhDRbar2, double g, double tanb) const
	one loop H1 tadpole contributions from Higgs bosons in the loops

double	doCalcTad2Higgs (double q, double costhDRbar2, double g, double tanb) const

double	doCalcTad1Neutralinos (double q, double costhDRbar, double g, double cosb) const
	one loop H1 tadpole contributions from Neutralinos in the loops

double	doCalcTad2Neutralinos (double q, double costhDRbar, double g, double sinb) const
	one loop H2 tadpole contributions from Neutralinos in the loops

double	doCalcTad1Charginos (double q, double g, double cosb) const
	one loop H1 tadpole contributions from Charginos in the loops

double	doCalcTad2Charginos (double q, double g, double sinb) const

double	doCalcTad1GaugeBosons (double q, double costhDRbar2, double g, double tanb) const
	one loop H1 tadpole contributions from Charginos in the loops

double	doCalcTad2GaugeBosons (double q, double costhDRbar2, double g, double tanb) const

void	doTadpoles (double mt, double sinthDRbar)
	Does the full 2-loop calculation of both tadpoles and sets them. More...

virtual double	doCalcTadpole1oneLoop (double mt, double sinthDRbar)
	Does the calculation of one-loop pieces of \( t_1 / v_1 \). More...

virtual double	doCalcTadpole2oneLoop (double mt, double sinthDRbar)
	Does the calculation of one-loop pieces of \( t_2 / v_2 \). More...

DoubleMatrix	addStopQCD (double p, double mt, DoubleMatrix &strong)

DoubleMatrix	addStopStop (double p, double mt, DoubleMatrix &stop)

DoubleMatrix	addStopSbottom (double p, double mt, DoubleMatrix &sbottom)

DoubleMatrix	addStopHiggs (double p, double mt, DoubleMatrix &higgs)

DoubleMatrix	addStopEweak (double p, DoubleMatrix &electroweak)

DoubleMatrix	addStopNeutralino (double p, double mt, DoubleMatrix &neutralino)

DoubleMatrix	addStopChargino (double p, DoubleMatrix &chargino)

DoubleMatrix	addSdownQCD (double p1, double p2, int family, DoubleMatrix &strong)

DoubleMatrix	addSdownHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

DoubleMatrix	addSdownEweak (double p1, double p2, int family, DoubleMatrix &electroweak)

DoubleMatrix	addSdownNeutralino (double p1, double p2, int family, DoubleMatrix &neutralino)

DoubleMatrix	addSdownChargino (double p1, double p2, int family, DoubleMatrix &chargino)

DoubleMatrix	addSbotQCD (double p, double mt, DoubleMatrix &strong)

void	addSbotSfermion (double p, double mt, DoubleMatrix &stop, DoubleMatrix &sbottom)

DoubleMatrix	addSbotHiggs (double p, double mt, DoubleMatrix &higgs)

DoubleMatrix	addSbotEweak (double p, DoubleMatrix &electroweak)

DoubleMatrix	addSbotNeutralino (double p, double mt, DoubleMatrix &neutralino)

DoubleMatrix	addSbotChargino (double p, double mt, DoubleMatrix &chargino)

DoubleMatrix	addSlepHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

DoubleMatrix	addSlepEweak (double p1, double p2, int family, DoubleMatrix &electroweak)

void	addSlepGaugino (double p1, double p2, int family, DoubleMatrix &chargino, DoubleMatrix &neutralino)

void	addSlepCorrection (double p, DoubleMatrix &mass, int family)

void	addStauSfermion (double p, double mtau, DoubleMatrix &stop, DoubleMatrix &sbottom)

void	addStauGaugino (double p, double mtau, DoubleMatrix &chargino, DoubleMatrix &neutralino)

DoubleMatrix	addStauEweak (double p, double mtau, DoubleMatrix &electroweak)

DoubleMatrix	addStauHiggs (double p, double mtau, DoubleMatrix &higgs)

DoubleMatrix	addSupQCD (double p1, double p2, int family, DoubleMatrix &strong)

DoubleMatrix	addSupHiggs (double p1, double p2, int family, DoubleMatrix &higgs)

DoubleMatrix	addSupEweak (double p1, double p2, int family, DoubleMatrix &electroweak)

DoubleMatrix	addSupNeutralino (double p1, double p2, int family, DoubleMatrix &neutralino)

DoubleMatrix	addSupChargino (double p1, double p2, int family, DoubleMatrix &chargino)

void	addSnuTauSfermion (double p, double &stop, double &sbottom)

double	addSnuTauHiggs (double p, double &higgs)

double	addSnuTauEweak (double p, double &electroweak)

void	addSnuTauGaugino (double p, double &chargino, double &neutralino)

double	addSnuHiggs (double p, int family, double &higgs)

double	addSnuEweak (double p, int family, double &electroweak)

void	addSnuGaugino (double p, int family, double &chargino, double &neutralino)

void	addSquarkCorrection (DoubleMatrix &mass)
	Now these are calculated at the squark scale. More...

virtual void	doUpSquarks (double mt, double pizztMS, double sinthDRbarMS, int accuracy)

virtual void	doDownSquarks (double mb, double pizztMS, double sinthDRbarMS, int accuracy, double mt)

virtual void	doChargedSleptons (double mT, double pizztMS, double sinthDRbarMS, int accuracy)

virtual void	doSnu (double pizztMS, int accuracy=0)

void	treeSnu (double &mSq, double pizztMS, int family)

void	setNeutCurrCouplings (double sinthDRbar, double &sw2, double &guL, double &gdL, double &geL, double &guR, double &gdR, double &geR)

void	calcDRTrilinears (drBarPars &eg, double vev, double beta)
	sets the Yukawas and Trilinears

void	calcDrBarHiggs (double beta, double mz2, double mw2, double sinthDRbar, drBarPars &eg)

void	treeCharginos (DoubleMatrix &mCh, double beta, double mw)

void	treeNeutralinos (DoubleMatrix &mN, double beta, double mz, double mw, double sinth)

void	calcDrBarGauginos (double beta, double mw, double mz, double sinth, drBarPars &eg)
	masses in eg with the appropriate values on exit.

virtual void	sparticleThresholdCorrections (double tb)

double	qedSusythresh (double alphaEm, double Q) const

double	qcdSusythresh (double alphasMSbar, double Q)

double	calcMs () const
	DRbar pars should be defined already for this. More...

virtual double	calcRunMtQCD () const
	Obtains (16 Pi^2) / mt times 1-loop and 2-loop QCD corrections. More...

virtual double	calcRunMtStopGluino () const
	Obtains (16 Pi^2) / mt times 1-loop stop-Gluino corrections. More...

virtual double	calcRunMtCharginos () const
	Obtains (16 Pi^2) / mt times 1-loop Charginos corrections.

virtual double	calcRunningMt ()

virtual double	calcRunMtauDrBarConv () const
	calculates factor to convert to DrBar for mtau

virtual double	calcRunMtauCharginos (double mTauSMMZ) const

virtual double	calcRunningMtau ()
	Full BPMZ expression. More...

virtual double	calcRunMbDrBarConv () const
	calculates factor to convert to DrBar for mb

virtual double	calcRunMbSquarkGluino () const
	Obtains (1 / mb) times 1-loop sbottom-Gluino corrections.

virtual double	calcRunMbChargino () const
	Obtains (1 / mb) times 1-loop squark-chargino corrections.

virtual double	calcRunningMb ()

double	calcSinthdrbar () const
	Calculates DRbar sin theta_w at the current scale from gauge couplings. More...

double	getVev ()
	Calculates Higgs VEV parameter from gauge couplings and MZ.

double	getVev (double pizzt)

virtual void	charginos (int accuracy, double piwwt)
	checked More...

virtual void	addChaLoopSfermion (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS) const

virtual void	addChaLoopGauge (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS, DoubleMatrix b1pCha, DoubleMatrix b0pCha, DoubleMatrix b1pNeut, DoubleMatrix b0pNeut) const

virtual void	neutralinos (int accuracy, double piwwt, double pizzt)
	mixNeut set to diagonal = mixNeut^T mNeutralino mixNeut: checked More...

void	addNeutLoopSfermion (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)

void	getNeutPassarinoVeltman (double p, double q, DoubleMatrix &b0fn, DoubleMatrix &b1fn)
	LCT: Returns matrix of Passarino-Veltman B0 and B1 functions. More...

void	addNeutLoopGauge (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)

void	addNeutLoopHiggs (double p, DoubleMatrix &sigmaL, DoubleMatrix &sigmaR, DoubleMatrix &sigmaS)

virtual void	gluino (int accuracy)
	Calculates pole gluino mass to 1-loop SUSY corrections. More...

void	calcHiggsAtScale (int accuracy, double &mt, double &sinthDRbar, double &piwwtMS, double &pizztMS, double q=0.)

bool	higgs (int accuracy, double piwwt, double pizzt)

int	rewsbM3sqTree (double, double &) const
	returns 1 if mu < 1.0e-9 More...

void	alternativeEwsb (double mt)
	Fixes mH1^2 and mH2^2 in order to get muCond and mAcond correct. More...

virtual void	rewsb (int sgnMu, double mt, const DoubleVector &pars, double muOld=-6.66e66, double eps=0.)

double	treeLevelMuSq ()

void	iterateMu (double &munew, int sgnMu, double mt, int maxTries, double pizztMS, double sinthDRbar, double tol, int &err)

double	predTanb (double muSusy=-6.66e66) const

double	predMzsq (double &tanb, double muOld=-6.66e66, double eps=0.)

double	deltaEW () const
	Returns fine-tuning as defined in Baer et al arXiv:1212.2655 eq 1.2. More...

double	ewsbCondition1TreeLevel () const
	DH: returns the values of the tree-level EWSB conditions. More...

double	ewsbCondition2TreeLevel () const

void	ewsbConditions (DoubleVector &values) const

void	predVevs (DoubleVector &vevs, int &err)
	DH: solves the EWSB conditions for the Higgs VEVs. More...

DoubleVector	fineTune (void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &bcPars, double MX, bool doTop=false)
	Pass it an object and it'll return the fine tuning parameters. More...

double	it1par (int numPar, const DoubleVector &bcPars)

double	ufb3sl (double)

double	realMinMs () const

double	calcBayesianNaturalness () const
	Calculates Bayesian naturalness a la arXiv:1709.07895.

virtual double	piZZTsfermions (double p, double Q) const
	Calculates sfermion contribution to the transverse part of Z self-energy: More...

virtual double	piZZTfermions (double p, double Q, bool usePoleMt) const
	fermion contribution to the Transverse part of Z self-energy: More...

virtual double	piZZTCharginos (double p, double q, double thetaWDRbar) const

virtual double	piWWTfermions (double p, double Q, bool usePoleMt) const

virtual double	piWWTsfermions (double p, double Q) const

double	piHpHmFermions (double p, double q) const
	LCT: Returns fermion contributions to charged Higgs self-energy. More...

double	piHpHmSfermions (double p, double q, double mu) const
	LCT: Returns sfermion contributions to charged Higgs self-energy. More...

double	piHpHmGauge (double p, double q) const
	LCT: Returns pure gauge contributions to charged Higgs self-energy. More...

double	piHpHmHiggs (double p, double q) const
	LCT: Returns Higgs contributions to charged Higgs self-energy. More...

double	piHpHmGauginos (double p, double q) const

double	piZGT (double p, double Q) const

virtual double	piAA (double p, double Q) const

Complex	pis1s1Sfermions (double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau) const
	Calculates (16 Pi^2) times the sfermion contribution to Higgs. More...

Complex	pis1s2Sfermions (double p, double q, DoubleMatrix ls1tt, DoubleMatrix ls1bb, DoubleMatrix ls1tautau, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const

Complex	pis2s2Sfermions (double p, double q, DoubleMatrix ls2tt, DoubleMatrix ls2bb, DoubleMatrix ls2tautau) const

Complex	pis1s1Fermions (double p, double q) const

Complex	pis2s2Fermions (double p, double q) const

Complex	pis1s1Higgs (double p, double q) const

Complex	pis1s2Higgs (double p, double q) const

Complex	pis2s2Higgs (double p, double q) const

Complex	pis1s1Neutralinos (double p, double q) const

Complex	pis1s2Neutralinos (double p, double q) const

Complex	pis2s2Neutralinos (double p, double q) const

Complex	pis1s1Charginos (double p, double q) const

Complex	pis1s2Charginos (double p, double q) const

Complex	pis2s2Charginos (double p, double q) const

Complex	pis1s1 (double p, double q) const

Complex	pis1s2 (double p, double q) const

Complex	pis2s2 (double p, double q) const
	checked 28.10.02 More...

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)

virtual double	deltaVb (double outrho, double outsin, double alphaDRbar, double pizztMZ) const
	outrho, outsin represent the DRbar values More...

virtual double	dRho (double outrho, double outsin, double alphaDRbar, double pizztMZ, double piwwtMW)

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

int	lsp (double &mass, int &posi, int &posj) const

int	nlsp (double &mass, int &posi, int &posj) const

string	printShort () const
	Prints a list of important sparticle/Higgs masses to standard output. More...

string	printLong ()
	Prints a list of all sparticle/Higgs masses to standard output. More...

string	printLongDrbar ()
	Prints a list of all running sparticle/Higgs masses to standard output. More...

double	thet (double a, double b, double c=0.0)
	log(max(a^2, b^2, c^2) / Q^2)

void	fixedPointIteration (TMSSMBoundaryCondition boundaryCondition, double mxGuess, const DoubleVector &pars, int sgnMu, double tanb, const QedQcd &oneset, bool gaugeUnification, bool ewsbBCscale=false)

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, int sgnMu, double tol, double tanb, void(*boundaryCondition)(MssmSoftsusy &, const DoubleVector &), const DoubleVector &pars, bool gaugeUnification, bool ewsbBCscale)

virtual void	rpvSet (const DoubleVector &parameters)

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 (const char fname[80]) const
	prints a file into fname which acts as an input to isajet More...

void	ssrunInterface764Inside (const char fname[80], fstream &) const

void	ssrunInterface764 (const char fname[80], const char softfname[80]) const

void	isawigInterface764 (const char fnamein[80], const char fnameout[80], const char fnamesoft[80]) const

virtual void	setEwsbConditions (const DoubleVector &inputs)
	Again, another dummy - useful in alternative EWSB conditions sometimes. More...

void	headerSLHA (ostream &out)

void	modselSLHA (ostream &out, const char model[])
	This does the MODSEL block of SLHA. More...

void	sminputsSLHA (ostream &out)
	SMINPUTS block of SLHA.

void	minparSLHA (ostream &out, const char model[], const DoubleVector &pars, double tanb, int sgnMu, bool ewsbBCscale)
	MINPAR block of SLHA. More...

void	massSLHA (ostream &out)
	This does the job of the above method, but outputs the Mass block. More...

virtual void	sfermionsSLHA (ostream &out)
	sfermions part of mass block of SLHA

void	sfermionmixSLHA (ostream &out)
	hmix block output of SLHA

void	inomixingSLHA (ostream &out)
	This does the job of the above method, but outputs the UMIX/VMIX blocks. More...

void	alphaSLHA (ostream &out)
	sfermionic part of mixing blocks

virtual void	hmixSLHA (ostream &out)
	higgs part of mixing blocks

void	gaugeSLHA (ostream &out)
	sfermionic part of mixing blocks

virtual void	yukawasSLHA (ostream &out)
	sfermionic part of mixing blocks

virtual void	msoftSLHA (ostream &out)
	MSOFT block of SLHA.

void	doUfb3 (double mgut)

void	assignHiggs (DoubleVector &higgsm, DoubleVector &higgsc, DoubleVector &dnu, DoubleVector &dnd, DoubleVector &cn, double beta) const

void	assignHiggs (DoubleVector &higgsm, DoubleVector &higgsc) const

void	assignHiggsSfermions (DoubleVector &higgsm, DoubleVector &higgsc, DoubleVector &dnu, DoubleVector &dnd, DoubleVector &cn, double beta) const
	some switches due to BPMZ's different conventions More...

double	smPredictionMW () const
	checked 22/04/06 More...

double	twoLoopGm2 (double amu1Loop) const

double	twoLpMt () const

double	twoLpMb () const

virtual void	doQuarkMixing (DoubleMatrix &mDon, DoubleMatrix &mUpq)
	input diagonal matrices and it'll give you back mixed ones More...

virtual MssmSusyRGE	guessAtSusyMt (double tanb, const QedQcd &oneset)

Public Member Functions inherited from softsusy::MssmSusy
	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

Public Member Functions inherited from softsusy::MssmSoftPars
	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.

Public Member Functions inherited from softsusy::AltEwsbMssm
void	setAltEwsb (double ma, double mu)

void	setAltEwsbMssm (const AltEwsbMssm &s)

double	displayMaCond () const

double	displayMuCond () const

const AltEwsbMssm &	displayAltEwsbMssm () const

void	setMaCond (double maInput)

void	setMuCond (double muInput)

Public Member Functions inherited from softsusy::Approx
const Approx &	displayApprox () const

const Approx &	operator= (const Approx &)

	Approx (const Approx &a)

	Approx (int l, int t)

void	setLoops (int l)
	Set number of loops used.

void	setThresholds (int t)

int	displayLoops () const
	Return number of loops.

int	displayThresholds () const
	Return level of threshold approximation.

Public Member Functions inherited from softsusy::RGE
	RGE ()
	Default constructor fills data with zeroes.

void	setMu (double e)
	Sets renormalisation scale to e.

virtual void	setmu (double f)

void	setPars (int i)
	Set number of parameters in RGE object.

double	displayMu () const
	Return renomalisation scale.

int	howMany () const
	Return number of parameters.

const RGE &	displayRGE () const

virtual int	run (double from, double to, double eps=-1.0)
	Runs parameters from scale "from" to "to, eps is fractional accuracy.

virtual int	runto (double to, double eps=-1.0)
	Runs parameters to scale "to", eps is accuracy.

int	callRK (double x1, double x2, DoubleVector &v, DoubleVector(*derivs)(double, const DoubleVector &), double eps=-1.0)

Additional Inherited Members
Public Types inherited from softsusy::MssmSoftsusy
typedef void(*	TMSSMBoundaryCondition) (MssmSoftsusy &, const DoubleVector &)

Public Attributes inherited from softsusy::MssmSoftsusy
int	included_thresholds
	Flag allowing to choose which two-loop thresholds have to be included.

Protected Member Functions inherited from softsusy::MssmSoftsusy
void	setT1OV1Ms (double t1)

void	setT2OV2Ms (double t2)

void	setT1OV1Ms1loop (double t1)

void	setT2OV2Ms1loop (double t2)

void	check_flags ()
	check consistency of threshold flags More...

DoubleMatrix	calcHiggs3L (bool is_bottom)
	Higgs 3L corrections. More...

Protected Attributes inherited from softsusy::MssmSoftsusy
double	mxBC
	Scale at which SUSY breaking boundary conditions set.

Detailed Description

Contains all supersymmetric NMSSM parameters, incorporating R_p NMSSM.

Constructor & Destructor Documentation

◆ NmssmSoftsusy()

softsusy::NmssmSoftsusy::NmssmSoftsusy	(	const NmssmSusyPars &	ss,
		const SoftParsNmssm &	s,
		const MssmSoftsusy &	sp
	)

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

◆ addChaLoopHiggs()

virtual void softsusy::NmssmSoftsusy::addChaLoopHiggs	(	double	p,
		DoubleMatrix &	sigmaL,
		DoubleMatrix &	sigmaR,
		DoubleMatrix &	sigmaS,
		DoubleMatrix	b1pCha,
		DoubleMatrix	b0pCha,
		DoubleMatrix	b1pNeut,
		DoubleMatrix	b0pNeut
	)		const

virtual

PA: Fills sigmaL, sigmaR and sigmaS with NMSSM Higgs parts of the charginos loop corrections.

Reimplemented from softsusy::MssmSoftsusy.

◆ addNeutLoopHiggs()

void softsusy::NmssmSoftsusy::addNeutLoopHiggs	(	double	p,
		DoubleMatrix &	sigmaL,
		DoubleMatrix &	sigmaR,
		DoubleMatrix &	sigmaS
	)

LCT: Adds Higgs loop corrections to neutralino mass matrix. IO parameters: p=external momentum, {sigmaL,sigmaR,sigmaS}= 4x4 matrices of {left,right,scalar} corrections

◆ addSbotCorrection()

virtual void softsusy::NmssmSoftsusy::addSbotCorrection	(	double	p,
		DoubleMatrix &	mass,
		double	mt
	)

virtual

Adds one-loop corrections to sbottom 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

Reimplemented from softsusy::MssmSoftsusy.

◆ addSbotHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addSbotHiggs	(	double	p,
		double	mt,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop Higgs corrections to sbottom mass matrix IO parameters: p=external momentum, mt=DR bar top mass, outputs higgs=2x2 matrix of self-energies in LR basis

◆ addSdownCorrection()

virtual void softsusy::NmssmSoftsusy::addSdownCorrection	(	DoubleMatrix &	mass,
		int	family
	)

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ addSdownHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addSdownHiggs	(	double	p1,
		double	p2,
		int	family,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop higgs corrections to sdown mass matrix IO parameters: {p1,p2}=external momentum, outputs higgs=2x2 matrix of self-energies in LR basis

◆ addSlepCorrection()

virtual void softsusy::NmssmSoftsusy::addSlepCorrection	(	DoubleMatrix &	mass,
		int	family
	)

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ addSlepHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addSlepHiggs	(	double	p1,
		double	p2,
		int	family,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop Higgs corrections to slepton self energy IO parameters: p=external momentum, family=generation, outputs higgs=2x2 matrix of self-energies

◆ addSnuCorrection()

virtual void softsusy::NmssmSoftsusy::addSnuCorrection	(	double &	mass,
		int	family
	)

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ addSnuHiggs()

double softsusy::NmssmSoftsusy::addSnuHiggs	(	double	p,
		int	family,
		double &	higgs
	)

LCT: Adds one-loop Higgs corrections to sneutrino self energy IO parameters: p=external momentum, family=generation, outputs higgs=self-energy

◆ addSnuTauCorrection()

virtual void softsusy::NmssmSoftsusy::addSnuTauCorrection ( double & mass )

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ addSnuTauHiggs()

double softsusy::NmssmSoftsusy::addSnuTauHiggs	(	double	p,
		double &	higgs
	)

LCT: Adds one-loop Higgs corrections to tau sneutrino self energy IO parameters: p=external momentum, outputs higgs=self-energy

◆ addStauCorrection()

virtual void softsusy::NmssmSoftsusy::addStauCorrection	(	double	p,
		DoubleMatrix &	mass,
		double	mtau
	)

virtual

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, mtau=DR bar tau mass

Reimplemented from softsusy::MssmSoftsusy.

◆ addStauHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addStauHiggs	(	double	p,
		double	mtau,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop electroweak corrections to stau self energy IO parameters: p=external momentum, mtau=DR bar tau mass, outputs electroweak=2x2 matrix of self-energies

◆ addStopCorrection()

virtual void softsusy::NmssmSoftsusy::addStopCorrection	(	double	p,
		DoubleMatrix &	mass,
		double	mt
	)

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ addStopHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addStopHiggs	(	double	p,
		double	mt,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop higgs corrections to stop mass matrix IO parameters: p=external momentum, mt=DR bar top mass, outputs higgs=2x2 matrix of self-energies in LR basis

◆ addSupCorrection()

virtual void softsusy::NmssmSoftsusy::addSupCorrection	(	DoubleMatrix &	mass,
		int	family
	)

virtual

Adds one-loop corrections to sup mass matrix IO parameters: mass=tree level mass matrix on input, is returned with radiative corrections added for a given family

Reimplemented from softsusy::MssmSoftsusy.

◆ addSupHiggs()

DoubleMatrix softsusy::NmssmSoftsusy::addSupHiggs	(	double	p1,
		double	p2,
		int	family,
		DoubleMatrix &	higgs
	)

LCT: Adds one-loop higgs corrections to sup mass matrix IO parameters: {p1,p2}=external momentum, outputs higgs=2x2 matrix of self-energies in LR basis

◆ calcDrBarGauginos()

void softsusy::NmssmSoftsusy::calcDrBarGauginos	(	double	beta,
		double	mw,
		double	mz,
		double	sinth,
		drBarPars &	eg
	)

PA: calls treeCharginos and treeNeutralinos, performs diagonalisation and fills eg with appropriate values.

◆ calcRunMtauNeutralinos()

virtual double softsusy::NmssmSoftsusy::calcRunMtauNeutralinos ( double mTauSMMZ ) const

virtual

PA: routine to calculate the Neutralino loop contributions to running Mtau

Reimplemented from softsusy::MssmSoftsusy.

◆ calcTadpole1Ms1loop()

virtual void softsusy::NmssmSoftsusy::calcTadpole1Ms1loop	(	double	mt,
		double	sinthDRbar
	)

virtual

Calculates and sets the one-loop pieces of \( t_1 / v_1 \): sets both 1-loop and total pieces equal to the one-loop piece

Reimplemented from softsusy::MssmSoftsusy.

◆ calcTadpole2Ms1loop()

virtual void softsusy::NmssmSoftsusy::calcTadpole2Ms1loop	(	double	mt,
		double	sinthDRbar
	)

virtual

Calculates then sets the one-loop pieces of \( t_2 / v_2 \): sets both 1-loop and total pieces equal to the one-loop piece

Reimplemented from softsusy::MssmSoftsusy.

◆ calcTadpoleSMs1loop()

virtual void softsusy::NmssmSoftsusy::calcTadpoleSMs1loop	(	double	mt,
		double	sinthDRbar
	)

virtual

Calculates then sets the one-loop pieces of \( t_s / s \): sets both 1-loop and total pieces equal to the one-loop piece Calculates then sets the one-loop pieces of \( t_s / s \): sets both 1-loop and total pieces equal to the one-loop piece

◆ charginos()

virtual void softsusy::NmssmSoftsusy::charginos	(	int	accuracy,
		double	piwwt,
		sPhysical &	phys
	)

virtual

Calculates pole chargino masses and mixing. IO parameters: piwwt is the W self-energy at the current, accuracy is the number of loops required (0 or 1 currently)

◆ DegrassiSlavicMix()

void softsusy::NmssmSoftsusy::DegrassiSlavicMix ( DoubleMatrix & P ) const

PA:: fixes The CP odd mixing matrix with the conventions Degrassi and Slavich arXiv:0907.4682

◆ doTadpoles()

void softsusy::NmssmSoftsusy::doTadpoles	(	double	mt,
		double	sinthDRbar
	)

PA: calls routines to calculate all three tadpoles and sets them. Currently only works at one loop. Two loop should be added later.

◆ ewsbConditions()

void softsusy::NmssmSoftsusy::ewsbConditions ( DoubleVector & values ) const

DH: returns the EWSB conditions at the current EWSB loop order. Note that this uses the current values of the tadpoles, i.e. it does not calculate or recalculate them.

◆ fineTune()

DoubleMatrix softsusy::NmssmSoftsusy::fineTune	(	void(*)(NmssmSoftsusy &, const DoubleVector &)	boundaryCondition,
		const DoubleVector &	bcPars,
		double	MX,
		bool	doTop = `false`
	)

DH: calculates fine-tuning for soft parameters, using the Barbieri-Giudice measure. 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. If doTop is true, it also calculates the fine tuning associated with the top Yukawa coupling. Returns a matrix in which the first column contains the calculated fine tunings, and the second contains the estimated errors on each.

◆ higgs()

bool softsusy::NmssmSoftsusy::higgs	(	int	accuracy,
		double	piwwt,
		double	pizzt,
		sPhysical &	phys
	)

Calculates pole Higgs masses and mixings: full 1-loop SUSY corrections leading log two loop for general nmssm and full effective potential 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 Returns "true" if there's a tachyon problem

◆ iterateMu()

void softsusy::NmssmSoftsusy::iterateMu	(	double &	munew,
		int	sgnMu,
		double	mt,
		int	maxTries,
		double	pizztMS,
		double	sinthDRbar,
		double	tol,
		int &	err
	)

PA: finds mu iteratively in the case where we use EWSB to swap.

Uses the full one loop tadpole from Degrassi and Slavich. No two loop added yet.

◆ iterateSvev()

void softsusy::NmssmSoftsusy::iterateSvev	(	double &	sold,
		int	sgnMu,
		double	mt,
		int	maxTries,
		double	pizzMS,
		double	sinthDRbar,
		double	tol,
		int &	err
	)

where the EWSB is used to swap (kappa, mS) --> (mZ, tb) and determine s. Uses the full one loop tadpole from Degrassi and Slavich. No two loop added yet.

◆ itLowsoft()

void softsusy::NmssmSoftsusy::itLowsoft	(	int	maxTries,
		int	sgnMu,
		double	tol,
		double	tanb,
		void(*)(NmssmSoftsusy &, const DoubleVector &)	boundaryCondition,
		const DoubleVector &	pars,
		const DoubleVector	nmpars,
		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 when lambda * svev = 0, the desired sign of svev when mu = 0 and the sign of in the quadratic equation for mu when both mu and lambda * svev are non-zero. tanb = desired value of DR bar tan beta(MZ).

◆ lesHouchesAccordOutput()

virtual void softsusy::NmssmSoftsusy::lesHouchesAccordOutput	(	ostream &	out,
		const char	model[],
		const DoubleVector &	pars,
		int	sgnMu,
		double	tanb,
		double	qMax,
		int	numPoints,
		bool	ewsbBCscale
	)

virtual

Outputs with Les Houches accord conventions to standard output. Inputs: out should be something like cout or fout depending on whether you want output in a file or not. model contains what form of model is used for the SUSY breaking terms (eg sugra, gmsb, amsb, nonUniversal). qMax is only relevant if you want a gridded output of running parameters up to some scale qMax. Put numPoints = 1 if you don't want to use this option - then qMaz is immaterial. mb is mb(mb) in the MSbar scheme used to produce the output, whereas mtau is the pole mass used (eg 1.777). mgut is the GUT scale that has been determined, and altEwsb is true if you specified mu and mA as input parameters (not tan beta and mH1, mH2).

Reimplemented from softsusy::MssmSoftsusy.

◆ neutralinos()

virtual void softsusy::NmssmSoftsusy::neutralinos	(	int	accuracy,
		double	piwwt,
		double	pizzt,
		sPhysical &	phys
	)

virtual

Calculates pole neutralino masses and mixing.

accuracy is the number of loops required (0 or 1 currently), pizzt is the Z self-energy at M_SUSY

◆ P1SfSfCouplings()

void softsusy::NmssmSoftsusy::P1SfSfCouplings	(	DoubleMatrix &	lTP1Lr,
		DoubleMatrix &	lBP1Lr,
		DoubleMatrix &	lTauP1Lr
	)		const

PA: obtains NMSSM P1-sfermion-sfermion couplings for 3rd generation sfermions

◆ physical()

virtual void softsusy::NmssmSoftsusy::physical ( int accuracy )

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ piHpHm()

double softsusy::NmssmSoftsusy::piHpHm	(	double	p,
		double	q
	)		const

virtual

LCT: Calculates the charged Higgs self-energy in mass basis: for p=external momentum, q=renormalisation scale

Reimplemented from softsusy::MssmSoftsusy.

◆ piHpm11Gaugino()

double softsusy::NmssmSoftsusy::piHpm11Gaugino	(	double	p,
		double	q
	)		const

LCT: Calculates (16 Pi^2) times gaugino contribution to (i, j) component of charged Higgs self-energy in the weak basis for p=external momentum, q=renormalisation scale

◆ piHpm11Higgs()

double softsusy::NmssmSoftsusy::piHpm11Higgs	(	double	p,
		double	q
	)		const

LCT: Calculates (16 Pi^2) times Higgs contribution to (i, j) component of the charged Higgs self-energy in the weak basis for p=external momentum, q=renormalisation scale

◆ pis1s1()

double softsusy::NmssmSoftsusy::pis1s1	(	double	p,
		double	q
	)		const

Calculates the Higgs self-energy: for p=external momentum, Q=renormalisation scale

◆ piWWT()

virtual double softsusy::NmssmSoftsusy::piWWT	(	double	p,
		double	Q,
		bool	usePoleMt = `false`
	)		const

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ piZZT()

virtual double softsusy::NmssmSoftsusy::piZZT	(	double	p,
		double	Q,
		bool	usePoleMt = `false`
	)		const

virtual

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

Reimplemented from softsusy::MssmSoftsusy.

◆ predVevs()

void softsusy::NmssmSoftsusy::predVevs	(	DoubleVector &	vevs,
		int &	err
	)

DH: solves the EWSB conditions for the Higgs and singlet VEVs, given the current values of all of the other parameters

◆ printall()

void softsusy::NmssmSoftsusy::printall ( ) const

PA: A print method used in development. I find it useful and easier to read than couting the normal display function or calling printlong etc.

◆ rewsb()

void softsusy::NmssmSoftsusy::rewsb	(	int	sgnMu,
		double	mt,
		double	muOld = `-6.66e66`,
		double	eps = `0.`
	)

ie (mu, m3sq, XiS) --> (mZ, tb, s)

PA: organises imposition of EWSB conditions. Currently at full one loop, two loop to be added. Currently works for two cases depending on Z3 switch

◆ rewsbM3sq()

virtual int softsusy::NmssmSoftsusy::rewsbM3sq	(	double	mu,
		double &	m3sq
	)		const

virtual

returns 1 if mu < 1.0e-9 PA: nmssm version for use in Z3 violating case.

Reimplemented from softsusy::MssmSoftsusy.

◆ rewsbmSsq()

virtual int softsusy::NmssmSoftsusy::rewsbmSsq ( double & mSsq ) const

virtual

PA:For Z3 invariant NMSSM when we solve for s, kappa and mS or for non universal Higgs

◆ rewsbMu()

virtual int softsusy::NmssmSoftsusy::rewsbMu	(	int	sgnMu,
		double &	mu
	)		const

virtual

Returns mu from rewsb requirement. returns 1 if there's a problem. Call at MSusy PA: To be used in general Z3 violating nmssm

Reimplemented from softsusy::MssmSoftsusy.

◆ rewsbSvev()

virtual int softsusy::NmssmSoftsusy::rewsbSvev	(	int	sgnMu,
		double &	svev
	)		const

virtual

PA: NMssm rewsb routine which fixes imn much the same way as mu is fixed in the Mssm using mueff = lambda * s / root For use in Z3 constrained version or when other scenarios where mu = 0

◆ rewsbSvevNoZ3()

virtual int softsusy::NmssmSoftsusy::rewsbSvevNoZ3 ( double & xiS ) const

virtual

BCA: fourth EWSB condition (for the singlet Higgs field) new with respect to the MSSM.

◆ set()

void softsusy::NmssmSoftsusy::set ( const DoubleVector & s )

virtual

Sets all parameters in an RGE object to the contents of s. It's an obligatory function for a derived class.

Implements softsusy::RGE.

◆ treeChargedSlepton()

virtual void softsusy::NmssmSoftsusy::treeChargedSlepton	(	DoubleMatrix &	mass,
		double	mTrun,
		double	pizztMS,
		double	sinthDRbarMS,
		int	family
	)

virtual

Returns tree-level down squark mass matrix in "mass" for NMSSM. 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

Reimplemented from softsusy::MssmSoftsusy.

◆ treeDownSquark()

virtual void softsusy::NmssmSoftsusy::treeDownSquark	(	DoubleMatrix &	mass,
		double	mbrun,
		double	pizztMS,
		double	sinthDRbarMS,
		int	family
	)

virtual

Returns tree-level down squark mass matrix in "mass" for NMSSM. 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

Reimplemented from softsusy::MssmSoftsusy.

◆ treeUpSquark()

virtual void softsusy::NmssmSoftsusy::treeUpSquark	(	DoubleMatrix &	mass,
		double	mtrun,
		double	pizztMS,
		double	sinthDRbarMS,
		int	family
	)

virtual

Returns tree-level up squark mass matrix in "mass" for NMSSM. 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

Reimplemented from softsusy::MssmSoftsusy.

◆ VhAtMin()

double softsusy::NmssmSoftsusy::VhAtMin	(	double	v1,
		double	v2,
		double	s,
		double	mHu2,
		double	mHd2,
		double	mSsq
	)

LCT: Returns Higgs potential at minimum. Inputs: {v1,v2,s}=Higgs vevs from EWSB conditions, {mHu2,mHd2,mSsq}=soft masses calculated at scale mu.

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

src/nmssmsoftsusy.h

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ NmssmSoftsusy()

Member Function Documentation

◆ addChaLoopHiggs()

◆ addNeutLoopHiggs()

◆ addSbotCorrection()

◆ addSbotHiggs()

◆ addSdownCorrection()

◆ addSdownHiggs()

◆ addSlepCorrection()

◆ addSlepHiggs()

◆ addSnuCorrection()

◆ addSnuHiggs()

◆ addSnuTauCorrection()

◆ addSnuTauHiggs()

◆ addStauCorrection()

◆ addStauHiggs()

◆ addStopCorrection()

◆ addStopHiggs()

◆ addSupCorrection()

◆ addSupHiggs()

◆ calcDrBarGauginos()

◆ calcRunMtauNeutralinos()

◆ calcTadpole1Ms1loop()

◆ calcTadpole2Ms1loop()

◆ calcTadpoleSMs1loop()

◆ charginos()

◆ DegrassiSlavicMix()

◆ doTadpoles()

◆ ewsbConditions()

◆ fineTune()

◆ higgs()

◆ iterateMu()

◆ iterateSvev()

◆ itLowsoft()

◆ lesHouchesAccordOutput()

◆ neutralinos()

◆ P1SfSfCouplings()

◆ physical()

◆ piHpHm()

◆ piHpm11Gaugino()

◆ piHpm11Higgs()

◆ pis1s1()

◆ piWWT()

◆ piZZT()

◆ predVevs()

◆ printall()

◆ rewsb()

◆ rewsbM3sq()

◆ rewsbmSsq()

◆ rewsbMu()

◆ rewsbSvev()

◆ rewsbSvevNoZ3()

◆ set()

◆ treeChargedSlepton()

◆ treeDownSquark()

◆ treeUpSquark()

◆ VhAtMin()