c --------------------------------------------------------------
c
c FeynHiggs
c =========
c
c Calculation of the masses of the neutral CP-even
c Higgs bosons in the MSSM
c
c Authors: Sven Heinemeyer (two-loop part)
c Andreas Dabelstein (one-loop part)
c
c Based on hep-ph/9803277, hep-ph/9807423, hep-ph/9812472,
c hep-ph/9903404, hep-ph/9910283
c by S. Heinemeyer, W. Hollik, G. Weiglein
c and on hep-ph/0001002
c by M. Carena, H. Haber, S. Heinemeyer, W. Hollik,
c C. Wagner and G. Weiglein
c
c In case of problems or questions,
c contact Sven Heinemeyer
c email: Sven.Heinemeyer@desy.de
c
c FeynHiggs homepage:
c http://www-itp.physik.uni-karlsruhe.de/feynhiggs
c
c --------------------------------------------------------------
This program provides the subroutine FeynHiggsSub in which the masses
of the neutral CP-even Higgs bosons in the MSSM are calculated.
Besides this also the mixing angle in the CP-even sector,
the charged Higgs mass as well as the SUSY contribution to the rho
parameter are eveluated.
All *.f files are needed for the computation.
FeynHiggs.f is the front-end which can be manipulated by the user at will.
What happenes in FeynHiggs.f? (FH == FeynHiggs)
1.) Start the program with 'FeynHiggs.exe'
Enter '0'
2.) Now the options are set (SETUP).
SETUP: ( * indicates one input)
======
* depth of calculation ?
1: full 1-loop + 2-loop QCD
2: same as 1, but in addition with
running mt at 2-loop
3: same as 2, but in addition with
two-loop Yukawa term included
(3 gives most accurate result)
[3 is recommended.]
* Select input:
1: On-Shell: Msusy_top_L,R, Xt, Msusy_bot_L,R, Xb
2: On-Shell: MSt1, MSt2, stt, MSb1, MSb2, stb
3: MSbar : MSt1, MSt2, stt, MSb1, MSb2, stb
[New feature: FH contains the evaluation for the Higgs masses worked out
in the On-Shell scheme. If your parameters are OS parameters, you can
choose 1 (unphysical Sfermion parameters) or 2 (physical Sfermion
parameters). If your parameters are MSbar parameters (e.g. the output from
some renormalization group equations) you can opt for 3. Then you can
enter the MSbar parameters supplemented with the scale at which they are
given.]
* mt in the OS stop mass matrix at 2-loop ?
--> irrelevant if input = 2, 3, 4 has been chosen
1: mt = pole mass
2: mt = running mass
[If 'Select input' has been chosen to be 1, you can specify what top quark
mass should be used to evaluate the Stop masses.]
* Are you interested in the 'fast and approximate'
result from the short approximation formula?
1: no, I prefer the full result only
2: yes, I'd like to get this additional information
[Do you want (e.g. for comparison) also the approximate result
(FeynHiggsFast) as explained in hep-ph/9903404? Then choose 2.]
* Limit for Delta rho = 1.3 * 10^-3 ? (0 = ok)
[The program also performs a check on \Delta\rho (see hep-ph/9612363
and hep-ph/9710438) where also the leading two-loop corrections
in the case of heavy gluino are implemented.
The program gives a warning if \Delta\rho exceeds the limit which
can be specified here (this indicates that this choice of sfermion
masses is experimentally disfavored.)
What additional information do you want?
* Test output of your input parameters? (0: yes // 1: no)
* Sfermion masses etc? (0: yes // 1: no)
* Neutralino/Chargino masses etc?(0: yes // 1: no)
[As a check you can get more information about the SUSY parameters.]
This ends the setup. After one calculation the user can choose whether to
change the setup or to perform the next calculation.
A possible input would look like:
3,1,1,2,0,1,1,1
3) Depending on the choice of 'Select input' you get three different types
of input masks. With the above choice you will get:
INPUT:
======
tan(beta)
Msusy_top_L
Msusy_top_R (0: Msusy_top_R = Msusy_top_L)
Msusy_bot_R (0: Msusy_bot_R = Msusy_top_L)
Xt == MtLR
Xb == MbLR (1: Ab = At)
Mtop(pole) (0: Mtop = 175 // 1: Mtop = 174.3)
Mbottom (0: Mbot = 4.5 // 1: Mbot = 2.97)
Mgluino (0: Mgl = 500 // 1: Mgl = Msusy_top_L)
Mu (1: Mu = -200 // 2: Mu = -Msusy_top_L)
M2 (1: M2 = 400 // 2: M2 = Msusy_top_L)
MA
Accuracy for the one-loop calculation:
1 = top only, 2 = top/bottom only, 3 = all
Now the user has to specify the MSSM parameters.
a) tan(beta) == v2/v1 (ratio of the two vevs).
b) input parameters of the stop mass matrix:
Msusy_top_L = Soft Susy breaking term in the upper left part of the
stop mass matrix
Msusy_top_R = SSB term in the lower right part of the stop mass matrix
c) input parameter of the sbottom mass matrix:
Msusy_bot_L is fixed via the SU(2) invariance to
Msusy_bot_L == Msusy_top_L. Only Msusy_bot_R is free to choose.
d) Xt defines the off-diagonal entry in the stop mass matrix.
It reads: Mtop * Xt == Mtop * (At - mu/tan(beta))
e) Same for the sbottom mass matrix: Xb == (Ab - mu*tan(beta))
In order to get reasonable values for Ab, one can choose '1'.
This fixes Ab = At, where At is calculated from Xt.
f) Mtop(pole): the top quark pole mass
g) Mbottom: the bottom quark mass.
The choice '1' (Mbot = 2.97) corresponds to the running bottom mass
at the scale Mtop. This choice incorporates automatically the leading
two-loop corrections of the bottom/sbottom sector into the calculation.
h) Mgluino is the gluino mass (surprise!)
i) Mu is the Higgs mixing parameter.
Convention in the stop mass matrix: Xt = At - mu/tan(beta)
Convention in the chargino mass matrix: +mu
Convention in the neutralino mass matrix: -mu
j) M2 is the SSB term in the chargino/neutralino sector.
k) MA is the mass of the CP-odd Higgs boson A.
l) The user can choose the accuracy of the one-loop part of the calculation:
A possible input would be:
3,1000,0,0,2000,1,1,1,800,200,1,500,3
Then the output should look like:
-------------------------------------------------
-------------------------------------------------
Performing the calculation...
using running mt for two-loop contribution: 166.665430375750
including two-loop Yukawa term
-------------------------------------------------
-------------------------------------------------
----------------------------------------------------
----------------------------------------------------
The results: light Higgs heavy Higgs alpha
----------------------------------------------------
mh-tree : 72.50753 503.0484 -0.3378600
----------------------------------------------------
mh-1loop:
--> BEST : 133.1238 504.7614 -0.3560725
short form.: 131.4415 504.6244 -0.3560664
Yuk-approx : 138.6594 504.9646
----------------------------------------------------
mh-2loop:
--> BEST : 116.3357 504.0734 -0.3489271
short form.: 115.6372 504.0988 -0.3500472
Yuk-approx : 122.9116 504.2478
----------------------------------------------------
----------------------------------------------------
charged Higgs, tree: 506.4213
1-loop: 506.1459
----------------------------------------------------
----------------------------------------------------
Delta rho 1-loop : 6.068775119477837E-005
Delta rho 2-loop (gluon) : 7.993995015988343E-006
Delta rho 2-loop (gluino): 0.000000000000000E+000
Delta rho total : 6.868174621076671E-005
----------------------------------------------------
collected WARNINGS:
----------------------------------------------------
The output for the Higgses should be self-explanatory. :-)
Delta rho gives only the SUSY contribution. The gluino-exchange part of
Delta rho can (only when it is negligible small!) be numerically instable.
In this case it is automatically set to zero.
The 'collected WARNINGS' contain possible warnings that can appear
in the Higgs output. Be aware! If it says something like
'Higgs sector not ok at 2-loop' then the result is NOT OK. Normally this
happens only if the parameter point under consideration is unphysical
(try e.g. Xt = 3 Msusy_top_L ...)
4.) The next round:
What next? (0: new calculation // 1: setup)
'0' takes you to
INPUT:
======
'1' takes you to
SETUP: ( * indicates one input)
======
As pointed out before, the program FeynHiggs.f is a front-end only.
So you can change everything at will in FeynHiggs, but some things always
have to be done: the SM parameters have to be set, the SUSY parameters have
to be set, including the four SSB terms in the stop and sbottom
mass matrices (msusytl, msusytr, msusybl, msusybr)
(To specify this a bit more: if the phsical stop masses are entered
directly, for the sake of [internal] simplicity, from these masses the
unphysical parameters are calculated in the subroutine def3:
call def3(ttb,mst2,delmst,stt,msusytl,msusytr,msusybl,msusybr,mtlr)
This automatically specifies the four SSB terms.)
You also have to specify the 'setup' which are explained here in (2)
The 'switches' set there are transferred to the subroutines via common blocks.
The subroutine FeynHiggsSub gives You the one- and two-loop values for
the Higgs-boson masses.
mh1: 1-loop value for the light Higgs
mh2: 1-loop value for the heavy Higgs
mh12: 2-loop value for the light Higgs
mh22: 2-loop value for the heavy Higgs
The tree-level masses are given by mlh (light Higgs) and mhh (heavy Higgs)
which are transferred via a common block.
Final note:
Please don't be too picky about how the program is written. It has been
developed over several years, the one-loop code was taken over from
A. Dabelstein, everything afterwards has been implemented into this
code. It is planned to rewrite the whole code within 2000.
Nevertheless the existing code works :-)