#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <cstdlib>
#include <cmath>
#include "def.h"
#include "utils.h"
#include "numerics.h"
#include "physpars.h"
#include "lowe.h"
#include "softpars.h"
#include "twoloophiggs.h"
Classes | |
class | MssmSoftsusy |
Functions | |
void | generalBcs (MssmSoftsusy &m, const DoubleVector &inputParameters) |
void | printShortInitialise () |
Prints out header line for print-short output. | |
ostream & | operator<< (ostream &, const MssmSoftsusy &) |
Formatted output. | |
double | sumTol (const MssmSoftsusy &in, const MssmSoftsusy &out) |
Calculates fractional difference in Drbar masses between in and out. | |
double | ccbSqrt (double f) |
returns the square root of the absolute value of the argument | |
void | recogLsp (int temp, int j) |
double | rho2 (double r) |
Two-loop Standard Model corrections to rho parameter. | |
MssmSusy | guessAtSusyMt (double tanb, const QedQcd &oneset) |
double | minimufb3 (double) |
double | ufb3fn (double mu, double htau, double h2, int family, const MssmSoftsusy &temp) |
double | getQhat (double inminTol, double eR, double h2, double Lisq, double mx, MssmSoftsusy &temp) |
void | extendedSugraBcs (MssmSoftsusy &m, const DoubleVector &inputParameters) |
non-universal mSUGRA boundary conditions | |
void | sugraBcs (MssmSoftsusy &m, const DoubleVector &inputParameters) |
universal mSUGRA boundary conditions | |
void | amsbBcs (MssmSoftsusy &m, const DoubleVector &inputParameters) |
Adds 2-loop AMSB boundary conditions onto m. | |
void | gmsbBcs (MssmSoftsusy &m, const DoubleVector &inputParameters) |
One-loop GMSB boundary conditions. | |
void | generalBcs2 (MssmSoftsusy &m, const DoubleVector &inputParameters) |
void | initialise () |
initialises factors (gfL,R) for loop calculations | |
double | gEff (double x) |
function used for calculating sin theta_eff | |
double | fEff (double x) |
function used for calculating sin theta_eff |
Revision 1.3 2005/07/26 11:00:42 allanach Added method nlsp to object MssmSoftsusy by Bernhardt
Revision 1.2 2005/07/15 15:10:06 allanach Added a routine to calculate sin^2 theta_eff
Revision 1.1.1.1 2004/11/19 16:18:31 allanach
Revision 1.65 2004/09/28 10:43:12 allanach Name-change to les Houches function - overloaded standard and file output functions.
Revision 1.64 2004/09/28 10:39:07 allanach Added Les Houches output function for output to a file
Revision 1.63 2004/04/20 13:55:08 allanach Calculation of tadpoles split up to allow calculation at different scales. The situation with mA^2(MZ) < 0 is handled by putting mApole and pretending that's the DRbar mass.
Revision 1.62 2004/03/21 20:43:05 allanach Added alternative electroweak symmetry breaking conditions. Added possibility of EWSB=SUSY breaking boundary condition scale in SUSY Les Houches Accord.
Revision 1.61 2004/03/17 13:23:06 allanach Solved 1-loop problem in alternative EWSB conditions. 2-loop problem remains.
Revision 1.60 2004/02/12 14:17:24 allanach Added tau^2 corrections and ht hb to higgs masses and Higgs vev calculation
Revision 1.59 2004/01/28 14:10:38 allanach Bug-fixes: for calculation sin theta_w, in the W and Z self-energies, the pole top mass is used in the top loop since the BPMZ SM two-loop corrections have assume that. Some errors were untrapped before, which is changed and gluino mass is allowed to be negative (as it is eg in AMSB).
Revision 1.58 2004/01/15 13:54:55 allanach New heaer style implemented
Revision 1.57 2003/12/02 19:15:53 allanach Added 3-family trilinear information, and added non-universal input
Revision 1.56 2003/11/06 14:41:47 allanach Les Houches interface changed to incorporate running mb
Revision 1.55 2003/10/27 15:50:39 allanach Taken old Higgs routine out of MssmSoftsusy, but added it to RpvSoft (there are currently problems using the new routine there -- needs to be fixed)
Revision 1.52 2003/10/24 10:07:42 allanach Any calculation of sin theta_w (DRbar) is now made from the EW gauge couplings
Revision 1.51 2003/08/19 14:38:59 allanach Altered so that mx in arguments to lowOrg is unchanged. It's used as the initial guess and lowOrg returns a double number as the calculated mgut.
Revision 1.50 2003/08/19 14:26:22 allanach Changing lowOrg to be more sensible about gauge unification. Should now be called with POSITIVE mgut and a flag for gauge unification.
Revision 1.49 2003/08/12 15:19:54 allanach Bug-fixed problem flags, iterateMu. Made iterateMu less precise (was sometimes causing problems)
Revision 1.48 2003/08/04 16:47:39 allanach Added function maxMass to return mass of heaviest superparticle, efficiency savings and 2-loop squark/gluino corrections to mt included
Revision 1.47 2003/07/30 13:39:02 allanach Corrected bug in piAA. Started to add piH+H-.
Revision 1.46 2003/07/28 12:11:37 allanach More error trapping, and rearranging rpvsoftsusy to use correct Higgs VEV (which is sometimes called at MZ)
Revision 1.45 2003/07/25 13:39:15 allanach Trapped errors properly rather than exiting
Revision 1.44 2003/07/24 14:55:28 allanach Implemented les Houches input and output properly in the usual command-line interface
Revision 1.43 2003/07/22 10:26:01 allanach gdL and gdR used in loop corrections are also 1-loop values now. This completes the DR-bar-ness of all 1-loop corrections in SOFTSUSY
Revision 1.42 2003/07/22 09:20:04 allanach Added display functions for running MW/MZ
Revision 1.41 2003/07/21 16:08:55 allanach displayHiggsVevMs now only called at MSUSY, and is checked @ that scale.
Revision 1.40 2003/07/21 14:00:18 allanach MZ fully implemented as an input now. Kept MZ as the central PDG 2002 value, for defaults etc
Revision 1.39 2003/07/18 15:39:14 allanach Added prediction of MW to definition of a softsusy object
Revision 1.38 2003/07/18 14:39:20 allanach Implemented MW as a global variable (in preparation for predicting it), and also speed corrections in getVev and rhohat: allowing input of self-energies to remove their calculation several times
Revision 1.37 2003/07/16 11:07:06 allanach Changed isajet number to 764
Revision 1.35 2003/06/05 09:17:19 allanach Started coding Les Houches Discord
Revision 1.34 2003/05/20 15:19:40 allanach doxygen comment style implemented
Revision 1.22 2002/11/27 18:19:16 allanach New higgs mass calculation and bugfix in piZZ (included neutrino self energy diagrams)
Revision 1.21 2002/11/26 13:25:59 allanach New mb calculation
Revision 1.20 2002/11/19 17:00:00 allanach Added FCNC routine
Revision 1.18 2002/10/22 13:12:10 allanach Introduced new problem flag for infra-red quasi fixed points
Revision 1.17 2002/10/14 14:17:30 allanach Added "runto" command in softpoint.x to get micromegas inputs to a different scale
Revision 1.16 2002/10/10 13:37:12 allanach Added physical parameters to Micromegas output
Revision 1.15 2002/10/09 17:18:51 allanach Added new softpoint option to specify mu and mA instead of mH1, mH2. Some fine-tuning still necessary
Revision 1.14 2002/09/20 15:37:10 allanach Adding quark-mixing routines
Revision 1.12 2002/09/04 13:59:45 allanach Added gauge unification possibility
Revision 1.11 2002/08/12 16:29:49 allanach display and set whole-object routines added
Revision 1.10 2002/07/30 12:57:32 allanach SOFTSUSY1.5
Revision 1.9 2002/06/14 16:26:30 allanach Switches included for 2-loop running of scalar masses, and calulating mt at mt.
Revision 1.7 2002/04/26 15:14:44 allanach Deleted all translation routines and defined boundary conditions within softsusy.h and softsusy.cpp
Revision 1.6 2002/04/14 13:50:41 allanach Now use V=m3^2 H1 H2 instead of V=mu B H1 H2. It's more natural!
Revision 1.5 2002/04/12 16:51:27 allanach Added display/set functions to work automatically
Revision 1.4 2002/04/12 06:24:50 allanach Code maintenance - returning a subobject made simpler
Revision 1.11 2002/02/04 15:07:22 allanach Latest version
Revision 1.10 2001/10/31 09:12:08 allanach Altered so that an early no rho convergence can disappear
Revision 1.9 2001/10/23 12:53:58 allanach Easier user interface utilised
Revision 1.8 2001/10/03 13:34:16 allanach Changed name of complex.h to avoid conflict with STD libraries
Revision 1.7 2001/09/28 13:59:58 allanach Split rhohat determination up into sensible pieces - small bugs fixed there.
Revision 1.6 2001/07/30 14:08:47 allanach Added ISAWIG and SSRUN interface
Revision 1.5 2001/07/26 14:20:34 allanach Added isajet interface
Revision 1.4 2001/07/18 15:54:50 allanach Put MIXING switch into def.h
Revision 1.3 2001/07/18 14:42:51 allanach Added proper header info
|
Allows user to specify a boundary condition where ALL SUSY breaking parameters are specified in inputParameters |
|
Sets all soft parameters in m except for mh1sq or mh2sq: it is intended for the case where mu and M_A^0(pole) is specified |
|
For UFB-3direction, returns scale at which one-loop corrections are smallest. IO parameters: inminTol is fractional accuracy with which minimum is found, eR is value of RH selectron field, h2 is value of H2 field, Lisq=|L_i|^2 slepton VEV value, mx=high BC-scale |
|
Provides the first guess at a SUSY object at mt, inputting tanb and oneset (should be at MZ) - it's very crude, doesn't take radiative corrections into account etc. oneset provides low energy data and tanb=tan beta |
|
For a given trial value of the log of field H2, gives the value of the potential at the minimum. The following global variables must be set before it is called: unificationScale=high BC scale, minTol=fractional accuracy with which minimum is found |
|
Prints out the identity of the LSP to standard output. temp, j are defined from lsp function |
|
Given mu paramer, tau Yukawa htau, family number examined, finds height of potential for temp at |H_2|=h2. |