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 :-)