FeynHiggs
FeynHiggs - the FeynHiggs command-line frontend
FeynHiggs parameterfile [flags]
FeynHiggs calculates the masses and couplings of the Higgs sector in the
MSSM. FeynHiggs provides a convenient way of invoking the subroutines
of the FeynHiggs library from the command line. The input parameters
are read from a file and the output is displayed in a human-readable
form on screen, but can be piped through a utility program to yield ma-
chine-readable data files, e.g. for plotting.
The input parameters for FeynHiggs are read from a file. FeynHiggs at-
tempts to read an SLHA file first and if that fails, falls back onto its
native format (description follows below). Output-wise, the default is
to write out an SLHA file (original filename + .fh-NNN) if the input was
one, and write the output to screen otherwise. Output style can be
changed by appending #SLHA to the filename to force SLHA output, and #FH
to force screen output.
Description of the native FeynHiggs format: comments begin with a # and
continue to the end of the line. Apart from a possible comment, each
line has the format
var lower [upper] [*][/][step]
and declares a loop over var from lower to upper in steps of step. A
"*" in front of step results in a logarithmic stepping, i.e. var is
multiplied by step in each iteration. A prefix of "/" indicates that
step is to be interpreted as the number of steps, rather than the step
size.
Simple assignments of the form name lower (i.e. no loop) are just a spe-
cial case of this since the items in brackets may be omitted.
var is the name of the parameter, which is one of
-------------------+--------------------------------------
invAlfa0, | the inverse of the electromagnetic
invAlfaMZ | coupling constant at 0, MZ
-------------------+--------------------------------------
AlfasMZ | the strong coupling constant at MZ
-------------------+--------------------------------------
GF | the Fermi constant
-------------------+--------------------------------------
MS, MC, MT, MB | the strange, charm, top, bottom masses
-------------------+--------------------------------------
MW, MZ | the W and Z masses
-------------------+--------------------------------------
TB | tan(beta)
-------------------+--------------------------------------
MA0, MHp | the CP-odd and charged Higgs mass
| (only one should be specified)
-------------------+--------------------------------------
MSusy | a generic soft-SUSY breaking mass
| used as a default value for all
| MgSL, MgSE, MgSQ, MgSU, MgSD.
-------------------+--------------------------------------
MgSL, MgSQ | individual soft-SUSY breaking masses
MgSE, MgSU, MgSD | for the slepton and squark doublets
g = 1..3 | and singlets.
| Default: all MgSx = MSusy.
-------------------+--------------------------------------
Abs(MUE), Re(MUE) | the Higgs mixing parameter
Arg(MUE), Im(MUE) |
-------------------+--------------------------------------
Abs(Af), Re(Af) | the sfermion mixing parameters
Arg(Af), Im(Af) |
f = e,mu,tau, |
u,c,t, |
d,s,b |
-------------------+--------------------------------------
Abs(Xf), Re(Xf) | alternative sfermion mixing parameter;
Arg(Xf), Im(Xf) | since Xf and Af are not independent,
f = tau,t,b | only one should be specified
-------------------+--------------------------------------
Abs(M_1), Re(M_1) | the U1 gaugino mass parameter
Arg(M_1), Im(M_1) |
-------------------+--------------------------------------
Abs(M_2), Re(M_2) | the SU2 gaugino mass parameter
Arg(M_2), Im(M_2) |
-------------------+--------------------------------------
Abs(M_3), Re(M_3) | the gluino mass parameter
Arg(M_3), Im(M_3) |
-------------------+--------------------------------------
Qtau, Qt, Qb | the scales at which the stau, stop,
| and sbottom DRbar masses are defined
| Default: all 0 (= on-shell masses)
-------------------+--------------------------------------
CKMlambda | the CKM parameters
CKMA |
CKMrhobar |
CKMetabar | Default: all -1 (= use PDG values)
-------------------+--------------------------------------
Abs(deltaFXYij) | the flavour-violating
Arg(deltaFXYij) | parameters
Re(deltaFXYij) |
Im(deltaFXYij) |
F = L,E,Q,U,D |
XY = LL,LR,RL,RR |
ij = 12,23,13, | Default: all 0
-------------------+--------------------------------------
scalefactor | the ren. scale is MT*scalefactor
-------------------+--------------------------------------
prodSqrts | sqrt(s) for the production x-sections
-------------------+--------------------------------------
Needless to add that complex quantities should be specified EITHER
through Abs and Arg OR through Re and Im. In conflicting cases a warn-
ing is issued and the Abs/Arg value takes precedence.
The flags are given as an optional 8-digit string on the command line,
where the default is 42020110. The positions in the string are as fol-
lows:
#1: mssmpart = 0..4
specifies the scope of the 1-loop part:
0: MT^4-approximation
1: top/stop-sector
2: top/stop- + bottom/sbottom-sector
3: full (s-)quark/lepton-sector
4: full MSSM
recommended: 4
#2: higgsmix = 1..3
determines the mixing in the Higgs sector:
1: All non-diagonal (ND) Self Energies = 0
2: All CP-violating ND Self Energies = 0
3: Full 3x3 mixing in the neutral sector
This means:
2: evaluation in the rMSSM
3: evaluation in the cMSSM
#3: p2approx = 0..4
determines the approximation for the 1-loop result:
0: None, i.e. full determination of the propagator matrices's
poles,
UHiggs is evaluated at p^2 = m^2
1: p^2 = 0 approximation
2: All self energies are calculated at p^2 = m_tree^2
3: imaginary parts of self-energies are discarded
4: as 0, but with UHiggs evaluated at p^2 = 0
recommended: 4
#4: looplevel = 0..2
determines the inclusion of higher-order corrections:
0: tree level
1: one-loop contributions only
2: include various two-loop contributions
recommended: 2
#5: loglevel = 0..3
determines the inclusion of log resummations:
0: no log resummation
1: NLL resummation (for large MSUSY)
2: NLL resummation (for large MCha,MNeu,MGlu,MSUSY)
3: NNLL resummation (for large MCha,MNeu,MGlu,MSUSY)
recommended: 3
#6: runningMT = 0..3
determines which top mass shall be used in the 1-/2-loop correc-
tions:
0: pole MT
1: running MT (SM MSbar 2L)
2: running MT (SM MSbar 1L)
3: running MT (MSSM DRbar)
recommended: 1
#7: botResum = 0..1
determines whether the O(tan beta^n) corrections shall be re-
summed:
0: no resummation
1: 1L resummation
2: 1L resummation for Higgs self-energies, 2L resummation for
Higgs decays
recommended: 2
#8: tlCplxApprox = 0..3
determines how the two-loop corrections are treated in the pres-
ence of complex parameters (cMSSM):
0: all corrections (asat, asab, atat, atab) are computed in
the rMSSM
1: only the cMSSM asat corrections are used
2: the cMSSM asat corrections are combined with the remaining
corrections in the rMSSM
3: the cMSSM a_s a_t corrections are combined with the remaining
corrections, interpolated in the complex phases (cf. FHSelec-
tIPol)
Using the example parameter file that comes with FeynHiggs, the command
FeynHiggs example/var.in
should result in the following Higgs masses and mixings:
| Mh0 = 113.89008624
| MHH = 195.87658697
| MA0 = 200.00000000
| MHp = 213.63955069
| SAeff = -0.35592463 0.00000000
| UHiggsRe = 0.99672973 0.08080741 0.00000000
| -0.08080741 0.99672973 0.00000000
| 0.00000000 0.00000000 1.00000000
and the following widths and branching ratios:
| h0-b-b = 6.018044E-03 0.825117 0.725474
| HH-b-b = 5.165486E-02 0.559950 2.526224E-03
| A0-b-b = 5.896187E-02 0.848461 2.280889E-03
| Hp-t-b = 6.224386E-02 0.781530
Before the results, FeynHiggs lists the parameters of the calculation.
With all the detailed information on widths and branching ratios, the
screen can become rather crowded. To mask off this detailed output,
type
FeynHiggs ... | grep -v %
The output of FeynHiggs can be converted into a format more amenable to
further processing, e.g. plotting, by the "table" utility. This works
as follows:
FeynHiggs ... | table TB Mh0 > var.out
var.out then contains two columns listing TB and Mh0 for each pass of
the loop.
The SPS benchmark scenarios are predefined. The input files can be
found in example/SPS/. To run a scenario, type e.g.
FeynHiggs example/SPS/SPS1a.in
For more information on the SPS, see
http://www.ippp.dur.ac.uk/~georg/sps/
The Les Houches benchmark scenarios for Higgs boson searches at hadron
colliders are predefined. The input files can be found in exam-
ple/LHBMS/. To run a scenario, type e.g.
FeynHiggs example/SPS/mhmax.in
For more information on the Les Houches Benchmark scenarios, see Eur.
Phys. J. C26 (2003) 601 [hep-ph/0202167].
FeynHiggs can handle data files in the SUSY Les Houches Accord (SLHA)
format. If the input file is in SLHA format, the results (masses,
widths, etc.) are equally written out in SLHA format, to a file named as
the source file suffixed by ".fh". For example:
FeynHiggs example/SLHA/TestEXTPAR.spc
produces example/SLHA/TestEXTPAR.spc.fh. FeynHiggs can process SLHA
files in which the input parameters are given (Block EXTPAR) as well as
files in which the masses and mixings is given (Blocks MASS, STOPMIX,
etc.), where in the latter case the input parameters are reconstructed
from the mass spectrum. For instance, try:
FeynHiggs example/SLHA/SphenoRP.spc
For more information on the SLHA, see
http://home.fnal.gov/~skands/slha/
The FeynHiggs debug level can be set using the FHDEBUG environment vari-
able, e.g.
setenv FHDEBUG 1 (in tcsh)
export FHDEBUG=1 (in bash)
The pager spawned to view FeynHiggs output is taken from the environment
variable PAGER and defaults to "less" if PAGER is not set. Should a
program abort occur with FeynHiggs, it is advisable to set PAGER to an
empty string, e.g.
setenv PAGER
This prevents a pager from being started and ensures that the output is
shown right up to the point where the abort occurs.
libFH(1), MFeynHiggs(1)
30-Jan-2018 FEYNHIGGS(1)