Updates an RPV object with loop-corrected neutrino data: use for L violation. More...

#include <rpvneut.h>

Inheritance diagram for softsusy::RpvNeutrino:

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

	RpvNeutrino (const RpvNeutrino &)
	Constructor initialises object equal to another one.

const RpvNeutrino &	operator= (const RpvNeutrino &)
	All data in object set equal to another one.

const DoubleVector &	displayPhysNuMasses () const
	Display functions to read data.

const DoubleMatrix &	displayUpmns () const

const DoubleMatrix &	displayPhysNeutMix () const

bool	displayInvertedOutput () const

const DoubleVector &	displayCPEMasses () const

const DoubleVector &	displayCPOMasses () const

const DoubleMatrix &	displayCPEMassSq () const

const DoubleMatrix &	displayCPOMassSq () const

const DoubleMatrix &	displayCPEMix () const

const DoubleMatrix &	displayCPOMix () const

void	setInvertedOutput ()

void	setNormalOutput ()

void	setCPEMasses (const DoubleVector &a)

void	setCPOMasses (const DoubleVector &a)

void	setCPEMassSq (const DoubleMatrix &a)

void	setCPOMassSq (const DoubleMatrix &a)

void	setCPEMix (const DoubleMatrix &a)

void	setCPOMix (const DoubleMatrix &a)

void	setUpmns (const DoubleMatrix &a)

DoubleVector	calculateSneutrinoVevs (const DoubleVector &sneutrinoVevs, double tol, double snuSq, double v1, double v2)
	replaces tree-level routine in rpvsoft.cpp - now includes 1-loop tadpoles More...

DoubleVector	calculateSneutrinoTadpoles (double sinthDRbar)

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

void	calculateLNVNeutralScalarMassMatrix (double sinthDRbar)

DoubleMatrix	calculateLNVSleptonMassMatrix (double sinthDRbar)

DoubleMatrix	calculateLNVDownSquarkMassMatrix (double sinthDRbar)
	down squark mass matrix - order: q2_i, d_j

DoubleMatrix	calculateLNVUpSquarkMassMatrix (double sinthDRbar)
	up squark mass matrix - order: q1_i, u_j

DoubleMatrix	neutrinoSeesaw ()

void	neutrinoData ()
	returns neutrino mixing angles and sets mass eigenstates to one loop. More...

void	AnalyticDiagonalise (DoubleMatrix &Mass, DoubleMatrix &MassMix, DoubleMatrix &MassMix0, DoubleMatrix &MassMix1, DoubleMatrix &MassMix2)

void	Positivise (DoubleMatrix &A)

void	FermionMassSym (DoubleMatrix &MM, ComplexMatrix &Z, DoubleVector &MassVector, DoubleVector &PositiveMassVector, DoubleVector &SquaredVector, DoubleMatrix &mix)

void	FermionMassAsym (DoubleMatrix &MM, ComplexMatrix &M, ComplexMatrix &P, DoubleVector &MassVector, DoubleVector &PositiveMassVector, DoubleVector &SquaredVector, DoubleMatrix &mixM, DoubleMatrix &mixP)

void	ScalarMass (DoubleMatrix &MM, ComplexMatrix &Z, DoubleVector &MassVector, DoubleVector &PositiveMassVector, DoubleVector &SquaredVector, DoubleMatrix &mix)

void	neutrino (DoubleVector &neutrino, DoubleMatrix &mixMnu, DoubleMatrix &mNeutrino)

void	sparticleThresholdCorrections (double tb, const DoubleVector &pars)

virtual void	physical (int accuracy)

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

void	spinfoSLHA (ostream &out)
	SPINFO block of SLHA output.

void	headerSLHA (ostream &out)
	Header part of SLHA output.

void	modselSLHA (ostream &out, const char model[])
	modsel block of SLHA output

void	massSLHA (ostream &out)
	mass block of SLHA2 output

void	inomixingSLHA (ostream &out)
	neutralino/chargino mixing for RPV More...

void	rvhmixSLHA (ostream &out)
	CP-even higgs mixing for RPV SLHA2. More...

void	rvamixSLHA (ostream &out)
	CP-odd higgs mixing for RPV SLHA2. More...

virtual void	selmixSLHA (ostream &out)
	charged slepton-higgs mixing for RPV SLHA2 More...

virtual void	snumixSLHA (ostream &out)
	uncharged slepton-higgs mixing for RPV SLHA2

virtual void	hmixSLHA (ostream &out)
	higgs part of DRbar couplings

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

virtual void	sleptonsSLHA (ostream &out)
	slepton part of mass block output of SLHA2

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

	RpvSoftsusy (const RpvSoftsusy &)
	Constructor initialises object equal to another one.

const RpvSoftsusy &	operator= (const RpvSoftsusy &)
	All data in object set equal to another one.

const DoubleVector	display () const
	Displays all RGE parameters in a double vector.

void	set (const DoubleVector &v)
	Sets all RGE parameters from elements of v.

DoubleVector	beta () const
	Beta functions of RPV MSSM.

RpvSoftsusy	beta2 () const
	Beta functions of RPV MSSM.

const RpvSoftsusy &	displayRpvSoftsusy () const
	Returns object as a const.

void	setRpvSoftsusy (const RpvSoftsusy &s)
	Sets whole object equal to another

void	rpvDisplay (DoubleVector &parameters) const

virtual void	rpvSet (const DoubleVector &parameters)

const DoubleVector &	displaySneutrinoVevs () const
	Returns the vacuum expectation values of sneutrinos.

const DoubleVector &	displayNeutrinoMasses () const
	Returns the neutrino masses.

const DoubleMatrix &	displayNeutralMixing () const
	Returns the full neutral mixing matrix.

const DoubleMatrix &	displayUch () const
	5 x 5 charged fermion mixing

const DoubleMatrix &	displayVch () const
	5 x 5 charged fermion mixing

int	displayNumRpcBcs () const
	Displays number of parameters in R-parity conserving BCs.

void	setNumRpcBcs (int pos)
	Sets number of parameters in R-parity conserving BCs.

void	setSneutrinoVevs (DoubleVector &v)
	Set vacuum expectation values of sneutrinos.

void	setNeutrinoMasses (DoubleVector &v)
	Set neutrino masses.

void	setNeutralMixing (DoubleMatrix &v)
	Set neutral fermion mixing matrix.

void	setUch (DoubleMatrix &v)
	Set charged fermion mixing matrix.

void	setVch (DoubleMatrix &v)
	Set charged fermion mixing matrix.

void	check (const DoubleVector &sneutrinoVevs) const
	Checks to what extent REWSB conditions are satisfied.

void	rpvAnomalousDimension (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, double &gH1H1, double &gH2H2, DoubleVector &gH1L) const
	Anomalous dimensions of fields in RPV MSSM.

void	rpvAnomalousDeriv (DoubleMatrix &gEE, DoubleMatrix &gLL, DoubleMatrix &gQQ, DoubleMatrix &gUU, DoubleMatrix &gDD, double &gH1H1, double &gH2H2, DoubleVector &gH1L) const

void	rpvyTildes (DoubleMatrix &ye, DoubleMatrix &yd, Tensor &letilde, Tensor &ldtilde, Tensor &lutilde) const

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

void	iterateRewsb (double &mu, double &m3sq, DoubleVector &sneutrinoVevs, int sgnMu, int &numTries, int maxTries, double tol, double mt, double oldMu, double eps)

double	calculateMu (const DoubleVector &sneutrinoVevs, int sgnMu, double v1, double v2)

virtual void	higgs (int accuracy, double piwwtMS, double pizztMS)

double	calculateM3sq (const DoubleVector &sneutrinoVevs, double snuSq, double v1, double v2)

void	rotateAwayVevs (DoubleVector &snVevs)

int	usefulVevs (double vSM, const DoubleVector &sneutrinoVevs, double &snuSq, double &v1, double &v2) const
	IO parameters: vSM = total VEVs of Higgs+sneutrinos (added in. More...

void	standardSugra (double m0, double m12, double a0)

DoubleMatrix	neutralinoMassMatrix () const
	Returns the 7 by 7 general RPV neutralino mass matrix.

void	neutralinos (int accuracy, double piwwtMS, double pizztMS)

void	methodBoundaryCondition (const DoubleVector &v)

virtual void	sparticleThresholdCorrections (double tb)

DoubleMatrix	chargedLeptons (double vev)
	Calculates charged lepton-chargino mass matrix including sneutrino VEVs. More...

void	iterateChargedLeptons (double vev, DoubleMatrix &yeOld, double tol, int maxTries, int &err, double mtau)

void	isawigInterface764 (char herwigInputFile[80], char isajetOutputFile[80], char softOutputFile[80], double eps=0.) const

void	outputNonZeroInputs (ostream &out)

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

void	drbarRpv (ostream &out)
	Does the RPV coupling SLHA1 outputs.

bool	leptonNumberViolation () const

Public Member Functions inherited from softsusy::FlavourMssmSoftsusy
	FlavourMssmSoftsusy ()
	constructor functions

	FlavourMssmSoftsusy (const FlavourMssmSoftsusy &s)

const FlavourMssmSoftsusy &	operator= (const FlavourMssmSoftsusy &s)

double	displayTheta12 () const
	standard set and display functions

double	displayTheta13 () const

double	displayTheta23 () const

double	displayDelta () const

double	displayThetaB12 () const

double	displayThetaB13 () const

double	displayThetaB23 () const

double	displayMnuE () const

double	displayMnuMu () const

double	displayPoleMe () const

double	displayPoleMmu () const

double	displayMd2GeV () const

double	displayMu2GeV () const

double	displayMs2GeV () const

double	displayMcMc () const

double	displayMnuTau () const

const flavourPhysical &	displayFlavourPhysical () const

const FlavourMssmSoftsusy &	displayFlavourSoftsusy () const

void	setFlavourSoftsusy (const FlavourMssmSoftsusy &s)

void	setTheta12 (double d)

void	setTheta13 (double d)

void	setTheta23 (double d)

void	setDelta (double d)

void	setThetaB12 (double d)

void	setThetaB13 (double d)

void	setThetaB23 (double d)

void	setMnuE (double d)

void	setMnuMu (double d)

void	setMnuTau (double d)

void	setPoleMe (double d)

void	setPoleMmu (double d)

void	setMd2GeV (double d)

void	setMu2GeV (double d)

void	setMs2GeV (double d)

void	setMcMc (double d)

void	setFlavourPhys (const flavourPhysical &s)
	Sets all squark flavour violating physical parameters.

void	setInitialData (const QedQcd &r)
	Sets the initial data from a un-run QedQcd object.

const DoubleMatrix	displayCkm () const

const ComplexMatrix	displayComplexCkm () const
	same as displayCkm, except the phase is properly incorporated

const DoubleMatrix	displayMns () const
	constructs a MNS matrix out of the angles stored within

void	modselSLHA (ostream &out, const char model[])
	MODSEL block of SLHA2.

void	sminputsSLHA (ostream &out)
	SM inputs block of SLHA2.

void	vckminSLHA (ostream &out)
	VCKMIN inputs block of SLHA2. More...

void	minparSLHA (ostream &out, const char model[], const DoubleVector &pars, double tanb, int sgnMu, bool ewsbBCscale)
	Input SUSY parameters part of SLHA2. More...

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

virtual void	sfermionsSLHA (ostream &out)
	sfermion part of mass block output of SLHA2 More...

virtual void	sfermionmixSLHA (ostream &out)
	sfermion mixing output part of SLHA2

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

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

virtual void	doUpSquarks (double mt, double pizztMS, double sinthDRbarMS, int accuracy)
	flavour violating squark mass/mixing calculation More...

virtual void	doDownSquarks (double mb, double pizztMS, double sinthDRbarMS, int accuracy, double mt)
	flavour violating squark mass/mixing calculation More...

virtual void	doChargedSleptons (double mT, double pizztMS, double sinthDRbarMS, int accuracy)
	flavour violating slepton mass/mixing calculation More...

virtual void	doSnu (double pizztMS, int accuracy=0)
	flavour violating slepton mass/mixing calculation More...

virtual void	doQuarkMixing (DoubleMatrix &mDon, DoubleMatrix &mUpq)

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

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

void	setAngles (double lambdaW, double aCkm, double rhobar, double etabar)

void	sCkm (DoubleMatrix &deltaULL, DoubleMatrix &deltaURR, DoubleMatrix &deltaULR, DoubleMatrix &deltaDLL, DoubleMatrix &deltaDRR, DoubleMatrix &deltaDLR) const

void	sCkmAnalytic (double m32, double &deltaULL23, double &deltaULR23, double &deltaULL12, double &deltaDLL23, double &deltaDLR23, double &deltaDLL12) const

void	sCkmRunning (DoubleMatrix &deltaULL, DoubleMatrix &deltaURR, DoubleMatrix &deltaULR, DoubleMatrix &deltaDLL, DoubleMatrix &deltaDRR, DoubleMatrix &deltaDLR)

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 void	calcTadpole1Ms1loop (double mt, double sinthDRbar)
	checked More...

virtual void	calcTadpole2Ms1loop (double mt, double sinthDRbar)

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)

virtual void	addStopCorrection (double p, DoubleMatrix &mass, double mt)
	As in BPMZ appendix. More...

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)

virtual void	addSdownCorrection (DoubleMatrix &mass, int family)

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)

virtual void	addSbotCorrection (double p, DoubleMatrix &mass, double mb)
	Checked 16.08.2005. More...

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)

virtual void	addSlepCorrection (DoubleMatrix &mass, int family)
	16.09.05 checked. More...

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)

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

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)

virtual void	addSupCorrection (DoubleMatrix &mass, int family)
	As in BPMZ appendix, INCLUDING weak boson loops. More...

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)

virtual void	addSnuTauCorrection (double &mass)

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)

virtual void	addSnuCorrection (double &mass, int family)
	Found+fixed bug 7/09/06. Thanks to J Kersten. More...

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

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

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	calcRunMtHiggs () const
	Obtains (16 Pi^2) / mt times 1-loop Higgs corrections.

virtual double	calcRunMtNeutralinos () const
	Obtains (16 Pi^2) / mt times 1-loop Neutralino 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	calcRunMtauNeutralinos (double mTauSMMZ) const

virtual double	calcRunMtauHiggs () const
	Obtains (1 / mtau) times 1-loop Higgs corrections.

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	calcRunMbHiggs () const
	Obtains (1 / mb) times 1-loop Higgs corrections.

virtual double	calcRunMbNeutralinos () const
	Obtains (1 / mb) times 1-loop neutralino-suqrk corrections. More...

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	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. 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	addNeutralinoLoop (double p, DoubleMatrix &)
	Adds the loop corrections on to an input tree-level neutralino mass. More...

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)

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

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

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	rewsbTreeLevel (int sgnMu)

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	piZZTHiggs (double p, double Q, double thetaWDRbar) const
	Calculates Higgs contribution to the transverse part of Z self-energy: More...

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	piZZTNeutralinos (double p, double Q, double thetaWDRbar) const

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

virtual double	piZZT (double p, double Q, bool usePoleMt=false) const
	Transverse part of Z self-energy: has been checked. More...

virtual double	piWWTHiggs (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

virtual double	piWWTgauginos (double p, double Q, double thetaWDRbar) const

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

virtual void	getNeutralinoCharginoHpmCoup (ComplexMatrix &apph1, ComplexMatrix &apph2, ComplexMatrix &bpph1, ComplexMatrix &bpph2) const
	LCT: Returns trilinear neutralino-chargino-hpm coupling in unrotated basis. More...

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

virtual double	piHpHm (double p, double Q) const
	LCT: Charged Higgs self-energy. More...

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 double	h1s2Mix ()
	gets h1 mixing element with Hu.

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

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

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)

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	neutralinoCharginoMSLHA (ostream &out)
	neutralino and charigno part of mass block of SLHA More...

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

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

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

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

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

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

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

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

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)

Public Member Functions inherited from softsusy::RpvSoftPars
	RpvSoftPars ()
	bilinear RPV parameters: 1..3 More...

	RpvSoftPars (const RpvSoftPars &)
	Constructor sets object with another as copy.

	RpvSoftPars (const DoubleVector &m, const Tensor &hu, const Tensor &hd, const Tensor &he, const DoubleVector &d)

const RpvSoftPars &	operator= (const RpvSoftPars &s)
	Sets contents of operator equal to those of another.

void	setDr (const DoubleVector &V)
	Sets bilinear RPV D parameters equal to a vector.

void	setD (int, double)
	Sets an individual bilinear RPV D parameter.

void	setMh1lSquared (int i, double f)
	Sets a single $ m_{H_1 L_i}^2 $.

void	setMh1lSquared (const DoubleVector &v)
	Sets all $ m_{H_1 L_i}^2 $.

void	setHr (RpvCouplings, int, int, int, double)
	Sets an element of $ h_U $=hu, $ h_D $=hd or $ h_E $=he.

void	setHr (RpvCouplings, const Tensor &)
	Sets all elements of $ h_U $=hu, $ h_D $=hd or $ h_E $=he.

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

void	setRpvSoftPars (const RpvSoftPars &s)
	Sets all parameters in object equal to those in another.

const RpvSoftPars &	displayRpvSoft () const
	Returns object as a const.

const Tensor &	displayHr (RpvCouplings) const
	Returns one of RPV $ h_{U,D,E} $ parameter tensors.

const DoubleVector &	displayMh1lSquared () const
	Returns vector of $ m_{H_1 L_i}^2 $.

const DoubleVector &	displayDr () const
	Returns vector of $D$ parameters.

void	checkOut (double tol) const

void	display (DoubleVector &v, int &k) const

Public Member Functions inherited from softsusy::RpvSusyPars
	RpvSusyPars ()
	bilinear RPV parameters More...

	RpvSusyPars (const RpvSusyPars &)
	Constructor copies another object.

	RpvSusyPars (const Tensor &hu, const Tensor &hd, const Tensor &he, const DoubleVector &d)

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

void	setLam (int i, int j, int k, double f)
	Conventional couplings: lambda_{ijk}.

void	setLamPrime (int i, int j, int k, double f)
	Conventional couplings: lambda'_{ijk}.

void	setLamPrimePrime (int i, int j, int k, double f)
	Conventional couplings: lambda'_{ijk}.

void	setKappa (int, double)
	Set one kappa value.

void	setLambda (RpvCouplings, int, int, int, double)
	Set a single element of a lambda tensor.

void	setKappa (const DoubleVector &k)
	Set all kappa values.

void	setLambda (const RpvCouplings &r, const Tensor &t)
	Set all of one tensor.

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

void	setRpvSusyPars (const RpvSusyPars &s)
	sets all parameters to be those in s

const Tensor &	displayLambda (RpvCouplings) const
	Returns one of the tensors LU, LD or LE.

const DoubleVector &	displayKappa () const
	Returns all kappa values in a vector.

double	displayLam (int i, int j, int k) const
	Conventional couplings: lambda_{ijk}.

double	displayLamPrime (int i, int j, int k) const
	Conventional couplings: lambda'_{ijk}.

double	displayLamPrimePrime (int i, int j, int k) const
	Conventional couplings: lambda''_{ijk}.

const RpvSusyPars &	displayRpvSusy () const
	Returns all object as a const.

void	display (DoubleVector &v, int &k) const

void	checkOut (double tol) const

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

Updates an RPV object with loop-corrected neutrino data: use for L violation.

Member Function Documentation

◆ AnalyticDiagonalise()

void softsusy::RpvNeutrino::AnalyticDiagonalise	(	DoubleMatrix &	Mass,
		DoubleMatrix &	MassMix,
		DoubleMatrix &	MassMix0,
		DoubleMatrix &	MassMix1,
		DoubleMatrix &	MassMix2
	)

routines defined specifically for "neutrino" below. Given a DoubleMatrix (Mass), perform analytic diagonalisation, assuming all off diagonal elements are much smaller than the difference between any diagonal entries. The resulting mixing matrices are returned used in RpvNeutrino::neutrino

Unit mixing matrix

Corrections of order (OffD/D) to the off diagonal entries

Corrections of order (OffD^2/D^2) to the diagonal entries

Off diagonal terms required to diagonalize completely to order OffD^2/D

◆ calculateLNVNeutralScalarMassMatrix()

void softsusy::RpvNeutrino::calculateLNVNeutralScalarMassMatrix ( double sinthDRbar )

LNV neutral scalar mass matrices (CPE/CPO) in Feynman gauge (xi=1), c.f. paper Allanach Kom

< Identity matrix

lower 3x3 part (sneutrinos)

◆ calculateLNVSleptonMassMatrix()

DoubleMatrix softsusy::RpvNeutrino::calculateLNVSleptonMassMatrix ( double sinthDRbar )

charged higgs - slepton mass matrix in xi=1 (feynman) gauge order: Hu1, Hd2, tildeLi, tildeEbarj gauge can be changed in the code

xi=1 Feynman gauge for loop calculations (m_G+=m+W); xi=0 gauge - m_G+=0

◆ calculateSneutrinoTadpoles()

DoubleVector softsusy::RpvNeutrino::calculateSneutrinoTadpoles ( double sinthDRbar )

next 3 routines calculate snu_i,h_d^0, h_u^0 tadpoles for minimisation of one-loop scalar potential snu_i tadpoles are also used for corrections to sneutrino vevs procedure: get full LNV mass matrices from routines below, for 1. CPE/CPO scalars 2. charged higgs-sleptons

up- and down-type squarks 4. quarks 5. gauge bosons
charginos-charged leptons 7. neutralinos-neutrinos then diagonalise to obtain mass eigenvalues and mixings then set-up couplings snu_i XY (or corresponding for h_d^0, h_u^0) in format matrix(X,Y) for 1. to 7. Use mixing matrices to rotate couplings to mass eigenbasis calculate tadpoles: for scalars: Sum (rotatedCoupling* a0(massEigenvalue))

for fermions: Sum (rotatedCoupling* massEigenvalue * a0(massEigenvalue)) / (tadpoles calculated with DRbar parameters and in the Feynman gauge (xi=1))

CPE/CPO scalar couplings

snu-higgs1-higg1 (down-type)

snu-higgs2-higg2 (up-type)

charged higgs-slepton mass matrix - order: Hu1, Hd2, tildeLi, tildeEbarj

slepton couplings

squark mass matrices

squark couplings

initialise

Quarks

root2 in last term?? (not in softsusy manual - I think we discussed this issue last year and agreed that it was missing in the manual)

Weak bosons

chargino-charged lepton mass matrix in basis (W-, Hd-, l-) M (W+, hu+, ebar+)

chargino couplings

neutralino-neutrino mass matrix

Swap rows and columns to the format (Bino, Wino, Hu, Hd, nui)

neutralino couplings

◆ calculateSneutrinoVevs()

DoubleVector softsusy::RpvNeutrino::calculateSneutrinoVevs	(	const DoubleVector &	sneutrinoVevs,
		double	tol,
		double	snuSq,
		double	v1,
		double	v2
	)

virtual

replaces tree-level routine in rpvsoft.cpp - now includes 1-loop tadpoles

An adaption of the original rpvsoftsusy version (Input a set of values for sneutrino VEVs and it returns a more accurate set)

Identity matrix

tree level sneutrino vevs (as before)

adding the 1-loop correction

Reimplemented from softsusy::RpvSoftsusy.

◆ doCalcTadpole1oneLoop()

double softsusy::RpvNeutrino::doCalcTadpole1oneLoop	(	double	mt,
		double	sinthDRbar
	)

virtual

Does the calculation of one-loop pieces of $ t_1 / v_1 $.

From hep-ph/9606211's appendix. It should be done at MSusy to minimize the 1-loop contributions. Only call if you've calculated drbarpars. inputs are running top/bottom masses: call at MSusy only

CPE/CPO scalar masses in order (Hu,Hd,snu_i)

neutral scalar couplings

higgs1-higgs2-higg2 (up-type)

charged higgs-slepton mass matrix - order: Hu1, Hd2, tildeLi, tildeEbarj

slepton couplings

squark mass matrices

squark couplings

if I include kk=1,2 3 and not only kk=3 as in softsusy.cpp,
result changes by 0.6%. Keep this??

SM quarks

Weak bosons

chargino-charged lepton mass matrix in basis (W-, Hd-, l-) M (W+, hu+, ebar+)

chargino couplings

neutralino-neutrino mass matrix

neutralino couplings

PA: stop, sbottom, stau, couplings in the left right basis

Higgs

Neutralinos

Charginos

Weak bosons

Reimplemented from softsusy::MssmSoftsusy.

◆ doCalcTadpole2oneLoop()

double softsusy::RpvNeutrino::doCalcTadpole2oneLoop	(	double	mt,
		double	sinthDRbar
	)

virtual

Does the calculation of one-loop pieces of $ t_2 / v_2 $.

From hep-ph/9311269's appendix. It should be done at MSusy to minimize the 1-loop contributions. Only call if you've calculated physpars inputs are running top/bottom masses. Call at MSusy

neutral scalar couplings

charged higgs-slepton mass matrix order: Hu1, Hd2, tildeLi, tildeEbarj

slepton couplings

squark mass matrices

squark couplings

SM quarks

Weak bosons

chargino-charged lepton mass matrix in basis (W-, Hd-, l-) M (W+,hu+,ebar+)

chargino couplings

neutralino-neutrino mass matrix

Swap rows and columns to the format (Bino,Wino, Hu, Hd, nui)

neutralino couplings

CHECKED

Sfermion couplings

Neutralinos

Reimplemented from softsusy::MssmSoftsusy.

◆ FermionMassAsym()

void softsusy::RpvNeutrino::FermionMassAsym	(	DoubleMatrix &	MM,
		ComplexMatrix &	M,
		ComplexMatrix &	P,
		DoubleVector &	MassVector,
		DoubleVector &	PositiveMassVector,
		DoubleVector &	SquaredVector,
		DoubleMatrix &	mixM,
		DoubleMatrix &	mixP
	)

given an asymmetric fermion mass matrix, output mixing matrix, mass eigenvector, positive mass eigenvector, squared mass eigenvector used in RpvNeutrino::neutrino

force the largest component of every mass eigenvector in the mixing matrices P and M to be positive. Since an asymmetric matrix is diagonalised by a bi-unitary diagonalisation, the diagonal values can always be made positive by changing phases of the mixing matrices. What we are really doing here is to 'un-do' these phase changes, and compensate this by forcing the mixing matrices to be imaginary. In this case the matrix MM is diagonalised by M.hermitianConjugate * MM * P = diag MM, where all eigenvalues are positive.

◆ FermionMassSym()

void softsusy::RpvNeutrino::FermionMassSym	(	DoubleMatrix &	MM,
		ComplexMatrix &	Z,
		DoubleVector &	MassVector,
		DoubleVector &	PositiveMassVector,
		DoubleVector &	SquaredVector,
		DoubleMatrix &	mix
	)

given a symmetric fermion mass matrix, output mixing matrix, mass eigenvector, positive mass eigenvector, squared mass eigenvector used in RpvNeutrino::neutrino

◆ inomixingSLHA()

void softsusy::RpvNeutrino::inomixingSLHA ( ostream & out )

neutralino/chargino mixing for RPV

need charginos too

We must re-order the rows and columns of the chargino mass matrix. SOFTSUSY is in the basis U = [ wino higgsino e_Ri ], whereas SLHA2 is in the basis U~ = [ e_Ri wino higgsino ].

◆ lesHouchesAccordOutput()

void softsusy::RpvNeutrino::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).

There are some re-orderings of masses required in the SLHA 2 version. CP+ scalars

goldstone boson in last place

charged scalars excluding goldstone

Reimplemented from softsusy::FlavourMssmSoftsusy.

◆ neutrino()

void softsusy::RpvNeutrino::neutrino	(	DoubleVector &	neutrino,
		DoubleMatrix &	mixMnu,
		DoubleMatrix &	mNeutrino
	)

Steve's "neutrino" calculates the full 1 loop 7x7 neutralino-neutrino mass matrix (now in rpvsoft.cpp basis) MH: corrected Higgs vevs; now in rpvsoftsusy basis Return tree + 1 loop neutrino–neutralino masses and mixings. calculation performed in ADD basis where Hd and L_i fields are on equal footing. At the end of the code the results are rotated to the rpvsoftsusy basis

square root of -1.

'zero'

loop integral external momentum scale

no. of colours

gauge fixing parameter (for reference only so do not change !)

OBTAIN TREE LEVEL PARAMETERS FROM MAIN SOFTSUSY CODE

gauge couplings and EW mixing angles

< = sqrt(oneset.displayAlpha(ALPHA)

4. * PI)/g2;

the Yukawa matrices

third family approximation: Useful when cross checked with RPC softsusy yu(1, 1) = 0.; yu(2, 2) = 0.; yd(1, 1) = 0.; yd(2, 2) = 0.; ye(1, 1) = 0.; ye(2, 2) = 0.; the VEVs

various VEVs

the tree level gauge boson msasses

the supersymmetric bilinear parameters

< (susy.h)

< (rpvsusypars.h)

the soft bilinear parameters

< (rpvsusypars.h)

< (softpars.h)

the msass mixing between down-higgs and sleptons

< (rpvsusypars.h)

the up- and down- type higgs squared, tadpoles included

the 4x4x3 lamda and 4x3x3 lamda' tensors

the 4x4x3 lambda and 4x3x3 lambda' rpv soft trilinear tensors

DEFINE TREE LEVEL MASS MATRICES up type quarks

down type quarks

up type squarks

down type squarks

neutralino-neutrino

a modified chargino-charged lepton matrix, Hd and Leptons are on equal footing

4x4 lepton doublet - down higgs doublet for CPeven, CPodd, and Charged scalar matrices

4x4 'sneutrino-neutral down higgs matrix for CPeven and CPodd msasses

CP even scalars

CP odd scalars (gauge-independent part)

Charged scalars (gauge-independent part)

gauge DEPENDENT part of the mass matrices (CP odd and charged scalars)

combine the gauge dependent and the gauge independent parts

OBTAIN THE mixing MATRICES AND MASS EIGENVALUES convention as in Dedes et al, hep-ph/0603225

get to the same convention as Dedes

SET ELECTRON AND MUON MASSES TO ZERO for easy comparison with BPMZ

< get to the same notation as Dedes (hep-ph/0603225)

< get to the same notation as Dedes

define the renormalization scales note actual renormalisation needs to be done at msusy otherwise get negative goldstone masses

define neutral fermion feynman rules all external neutrino/neutralino legs are in flavour basis all internal legs are in mass basis
in these feynman rules, 'k' denotes neutral fermion the first index correspond to a scalar, while the next two indices are fermions for vector boson feynman rules, only two indices, and both are fermions CP even neutral scalar – neutral fermion(external) – neutral fermion(internal)

CP odd neutral scalar – neutral fermion(external) – neutral fermion(internal)

in +ve scalar – in -ve LH fermion – neutral fermion

in -ve scalar – in +ve LH fermion – neutral fermion

in -ve scalar down – neutral fermion – in +ve LH down type (anti) quark

out +ve scalar up – neutral fermion – in +ve LH up type quark

in +ve scalar up – neutral fermion – in -ve LH up type (anti) quark

Z – neutral fermion(external) – neutral fermion(internal)

in W+ – in -ve LH fermion – out neutral fermion

in W- – in +ve LH fermion – out neutral fermion

calculate the loop corrections note in rpc softsusy a different qq is used for neutrino mass calculation pp=0 is a suitable choice in fact, pp=0 must be used in the present calculation, for consistency with mia approx. (we only have c0 at zero external momentum) sigma d, scalar contribution

CP even/odd scalar -neutral fermion -neutral fermion because the large cancellation between these two contributions lead to numerical instability in the 3x3 neutrino–neutrino block. that part will later be replaced by a cpe-cpo analytic mia calculation

charged scalar -charged fermion -neutral fermion

down squark -neutral fermion - down quark

up squark -neutral fermion -up quark

CPE+CPO scalar -neutral fermion -neutral fermion

charged scalar -charged fermion -neutral fermion

down squark -neutral fermion - down quark

up squark -neutral fermion -up quark

Sigma D, gauge boson contribution initializsation

Z boson -neutral fermion -neutral fermion

W boson -charged fermion -neutral fermion

Sigma L, gauge boson contribution

Z boson -neutral fermion -neutral fermion

W boson -charged fermion -neutral fermion

mia approximation for the cpe+cpo contribution to the 3x3 neutrino-neutrino block This is a systematic procedure to expand the combined CPE+CPO contributions in powers of RPV parameters (call it B), that is the (2x3) off block-diagonal part of the full 5x5 mass matrix. Currently the approximation is to order B^2, but is generalisable to include higher corrections if desired. The external mass scale is set to zero. CPE MIA approximation

put the h-H, sneutrinos and RPV contributions to different blocks

RPC rotation

the 5x5 mass matrix to be further diagonalised, after the RPC rotation

the 5x5 mixing (RPC part) (h-H mixing, sneutrino mixing)

CPO MIA approximation

put the G-A, sneutrinos and RPV contributions to different blocks

RPC rotation

the 5x5 mass matrix to be further diagonalised, after the RPC rotation

the 5x5 mixing (RPC part) (G-A mixing, sneutrino mixing)

define the full (RPC + RPV using MIA) mixing matrices in different orders of perturbations

define the sneutrino mass^2 differences (CPE minus CPO)

compatibility with wsr2-wsa2 : CHECKED

define the hH, GA, CPE-sneutrino, CPO-sneutrino masses after MIA

Check for Tachyonic states

sum over the loop contributions (the 3x3 block)

sneutrino mixing difference contributions wCPOsneuttot is used here, while using wCPEsneuttot will give practically the same results

sneutrino mass difference contributions

H, h, A cancellations.

combine all the loop corrections and obtain the corrected mass matrix the output loop corrected mass matrix is in flavour basis
combine the sigma D CPE+CPO scalar contribution ONLY

combine all the sigma D scalar contribution except the cpe+cpo contributions (i=1, 2)

combine all the sigma L scalar contribution

combine all the (sigma D+L) gauge boson contribution

Use formal matrix manipulations to calc gauge basis correction

This is for FULL FLAVOR BASIS LOOP CALCULATION

overwrite the neutrino-neutrino cpecpo part by the mia results

combine all loop contributions

include CPE-CPO cancellations

Calculate the loop corrected neutrino-neutralino masses

rotate to the RPV softsusy basis

obtain loop corrected eigenvalues etc

Main calculation ends here.

◆ neutrinoData()

void softsusy::RpvNeutrino::neutrinoData ( )

returns neutrino mixing angles and sets mass eigenstates to one loop.

in case of IH, swapping of rows in U & NeutrinoMasses needed

Calculate PMNS matrix by multiplying U by the charged lepton mixing matrix (transposed).

We may re-organise the signs in the mixing matrices to keep the diagonal elements positive

◆ neutrinoSeesaw()

DoubleMatrix softsusy::RpvNeutrino::neutrinoSeesaw ( )

using Steve's neutrino routine to obtain 3x3 effective neutrino mass matrix via seesaw

calculate 3x3 effective neutrino mass matrix via seesaw from 7x7 neutralino matrix

Set full 7 by 7 physical mixing matrix: defined at zero external momentum

◆ physical()

void softsusy::RpvNeutrino::physical ( int accuracy )

virtual

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

Organises calculation of physical quantities such as sparticle masses etc Call AT MSusy

Running masses at MSUSY

This part of the code adds higher loop corrections to gluino masses etc from Steve Martin et al

For third family sparticles, we must keep the same mass ordering as SOFTSUSY

Reimplemented from softsusy::MssmSoftsusy.

◆ rvamixSLHA()

void softsusy::RpvNeutrino::rvamixSLHA ( ostream & out )

CP-odd higgs mixing for RPV SLHA2.

Need to find out what basis ordering we should really be in for SLHA2. CP- scalars: find lightest boson that's not the goldstone etc

This is a check to see if the eigenvalues are in the right order cout << slhaMix * CPOscalars * slhaMix.transpose();

H1 <-> H2 difference in positioning of basis vectors

◆ rvhmixSLHA()

void softsusy::RpvNeutrino::rvhmixSLHA ( ostream & out )

CP-even higgs mixing for RPV SLHA2.

This is a check to see if the eigenvalues are in the right order cout << slhaMix * CPEscalars * slhaMix.transpose();

H1 <-> H2 difference in positioning of basis vectors

◆ ScalarMass()

void softsusy::RpvNeutrino::ScalarMass	(	DoubleMatrix &	MM,
		ComplexMatrix &	Z,
		DoubleVector &	MassVector,
		DoubleVector &	PositiveMassVector,
		DoubleVector &	SquaredVector,
		DoubleMatrix &	mix
	)

given squared mass matrix (MM), output mixing matrix, mass eigenvector, positive mass eigenvector, squared mass eigenvector used in RpvNeutrino::neutrino

FORCE the squared mass vector to be positive. This does not change anything if renormalization is done in Msusy. However this is potentially inconsistent if renormalization done at MZ ! This can be traced as the SquaredVector is still negative. At Mu(~MZ), some elements of wsa2, wsr2, wsch2 may become negative.

◆ selmixSLHA()

void softsusy::RpvNeutrino::selmixSLHA ( ostream & out )

virtual

charged slepton-higgs mixing for RPV SLHA2

charged higgs-slepton mass matrix - order: Hu, Hd, tildeLi, tildeEbarj

swap to SLHA2 order: Hu, Hd, tildeLi, tildeEbarj

Don't print out the line corresponding to the goldstone mode: line 8

Reimplemented from softsusy::FlavourMssmSoftsusy.

◆ sparticleThresholdCorrections()

void softsusy::RpvNeutrino::sparticleThresholdCorrections	(	double	tb,
		const DoubleVector &	pars
	)

next two routines are copy of old code which was taken out of the softsusy version when flavorsoft was implemented for objects RpvSoftsusy2, mixing=1 sets CKM mixing in the up-sector, mixing=2 sets CKM mixing in the down-sector pars contains RPV SUSY parameters for the point in question

Set the RPV couplings here if we need to: currently only works with CMSSM

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

src/rpvneut.h
src/rpvneut.cpp

Public Member Functions

Additional Inherited Members

Detailed Description

Member Function Documentation

◆ AnalyticDiagonalise()

◆ calculateLNVNeutralScalarMassMatrix()

◆ calculateLNVSleptonMassMatrix()

◆ calculateSneutrinoTadpoles()

◆ calculateSneutrinoVevs()

◆ doCalcTadpole1oneLoop()

◆ doCalcTadpole2oneLoop()

◆ FermionMassAsym()

◆ FermionMassSym()

◆ inomixingSLHA()

◆ lesHouchesAccordOutput()

◆ neutrino()

◆ neutrinoData()

◆ neutrinoSeesaw()

◆ physical()

◆ rvamixSLHA()

◆ rvhmixSLHA()

◆ ScalarMass()

◆ selmixSLHA()

◆ sparticleThresholdCorrections()