Real valued R-Parity violating MSSM object: use for UDD couplings. More...

#include <rpvsoft.h>

Inheritance diagram for RpvSoftsusy:

Public Member Functions
	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
const 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)
virtual DoubleVector	calculateSneutrinoVevs (const DoubleVector &sneutrinoVevs, double tol, 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.
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.
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)
void	drbarRpv (ostream &out)
	Does the RPV coupling SLHA1 outputs.
bool	leptonNumberViolation () const

Detailed Description

Real valued R-Parity violating MSSM object: use for UDD couplings.

Member Function Documentation

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

Returns value of m3sq consistent with sneutrino vevs given (which should be a length-3 vector) IO parameters: sgnMu=+/-1, sign of Higgs potential mu parameter, v1 and v2 are the two Higgs doublet VEVs and must be input and snuSq is the sum of squares of sneutrino VEVs

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

Returns value of mu consistent with sneutrino vevs given (which should be a length-3 vector) IO parameters: sgnMu=+/-1, sign of Higgs potential mu parameter, v1 and v2 are the two Higgs doublet VEVs and must be input

DoubleMatrix RpvSoftsusy::chargedLeptons ( double vev )

Calculates charged lepton-chargino mass matrix including sneutrino VEVs.

Calculates the charged lepton mass matrix once the leptonic Yukawas have been set. IO parameters: vev is um of slepton+Higgs Vevs in quadrature

void RpvSoftsusy::higgs	(	int	accuracy,
		double	piwwtMS,
		double	pizztMS
	)		`[virtual]`

Calculates pole Higgs masses and mixings: old feynHiggsfast calculation 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

Reimplemented from MssmSoftsusy.

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

Writes a file for input to ISAWIG, which will write HERWIG information on a file called herwigInputFile, ISAJET output in isajetOutputFile. The input file for ISAWIG is called softOutputFile. The RPV couplings will not be output if they are smaller than eps.

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

Iterates the charged lepton mass matrix in order to get pole lepton masses and sneutrino VEVs etc correct. IO: yeOld is input and a better approximation to it is output. vev is sum of slepton+Higgs Vevs in quadrature tol=desired fractional accuracy, maxTries=maximum number of iterations after which it bails out with zeroes in yeOld and an error flag in err!=0. mtau=the running DRbar tau mass

Save the charged lepton mixing matrices

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

Iterative solution to electroweak symmetry breaking boundary conditions. IO parameters: mu, m3sq are input and output Higgs potential parameters, iterated along with the sneutrino Vevs. Input sgnMu is the sign of mu, maxTries is the maximum number of iterations, tol is the desired fractional accuracy and mt is the DRbar running top mass

void RpvSoftsusy::methodBoundaryCondition ( const DoubleVector & v ) [virtual]

This is used to set both GUT-scale RPV parameters v(i=4,...) and the SUSY breaking ones: v(1)=m0, v(2)=m12, v(3)=a0 where IO parameters: m0=scalar mass, m12=gaugino mass, a0=trilinear coupling

Reimplemented from MssmSoftsusy.

void RpvSoftsusy::neutralinos	(	int	accuracy,
		double	piwwtMS,
		double	pizztMS
	)		`[virtual]`

Calculates neutralino masses. Uses RPC approximation for the 4 neutralinos then diagonalises the full 7 by 7 case for the neutrinos.

Reimplemented from MssmSoftsusy.

const RpvSoftsusy & RpvSoftsusy::operator= ( const RpvSoftsusy & s )

All data in object set equal to another one.

RP violating MSSM object including all real RPV MSSM couplings.

rpvsoft.cpp

Project: SOFTSUSY
Authors: Ben Allanach, Markus Bernhardt
Manual: hep-ph/0104145, Comp. Phys. Comm. 143 (2002) 305
Webpage: http://hepforge.cedar.ac.uk/softsusy/

void RpvSoftsusy::outputNonZeroInputs ( ostream & out )

For SLHA output: prints out the non-zero values of the inputted SUSY RPV couplings

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

Performs radiative electroweak symmetry breaking. IO parameters: sgnMu is +/-1, the sign of mu mt is the RUNNING DRbar top mass

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

Reimplemented from MssmSoftsusy.

void RpvSoftsusy::rotateAwayVevs ( DoubleVector & snVevs )

Input a set of values for sneutrino VEVs and it returns a more accurate set - the next step in the iteration.

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

Derivative of anomalous dimensions of fields with respect to renormalisation scale of RPV parts of RPV MSSM.

void RpvSoftsusy::rpvDisplay ( DoubleVector & parameters ) const

returns double vector of 9 (untouched) free entries followed by 112 RPV parameters the 9 entries are left free for general sugraBcs, amsbBcs, If you need *more parameters left, reserve the first pos (optional argument)

void RpvSoftsusy::rpvSet ( const DoubleVector & parameters ) [virtual]

sets RPV parameters from entry 10 onwards using vector produced by RpvSoftsusy::rpvDisplay() the 9 entries are left free for general sugraBcs, amsbBcs, ... parameters

Reimplemented from MssmSoftsusy.

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

Ytilde quantities are for calculational brevity in beta functions. They are all outputs here. This function returns the RPV parts (only) of the ytildes

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

Outputs with close to SLHA1 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).

implementing positive-only chargino masses in SLHA

void RpvSoftsusy::sparticleThresholdCorrections ( double tb ) [virtual]

Calculates DRbar values of gauge and Yukawa couplings from data, depending upon the spectrum and input value of tan beta=tb pars contains RPV SUSY parameters for the point in question

Reimplemented from MssmSoftsusy.

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

Theoretical boundary condition upon SUSY breaking terms. SUSY RPV parameters must be set before this is applied IO parameters: m0=scalar mass, m12=gaugino mass, a0=trilinear coupling

Reimplemented from SoftParsMssm.

int RpvSoftsusy::usefulVevs	(	double	vSM,
		const DoubleVector &	sneutrinoVevs,
		double &	snuSq,
		double &	v1,
		double &	v2
	)		const

IO parameters: vSM = total VEVs of Higgs+sneutrinos (added in.

v2 are the two Higgs doublet VEVs and snuSq is the sum of squares of sneutrino VEVs: all arguments are outputs except for the The VEVs that come out are in SOFTSUSY notation

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

rpvsoft.h
rpvsoft.cpp

Public Member Functions

Detailed Description

Member Function Documentation