DSZHiggs.h

Go to the documentation of this file.

// ====================================================================
// This file is part of Himalaya.
//
// Himalaya is licenced under the GNU General Public License (GNU GPL)
// version 3.
// ====================================================================

#pragma once

#ifdef __cplusplus
extern "C" {
#endif

/**
 * @file DSZHiggs.h
 * @brief Contains two-loop routines from Slavich et al.
 */

/**
 * Two-loop O(a_t a_s) corrections to the CP-even Higgs mass matrix.
 * Routine written by P. Slavich (e-mail: slavich@pd.infn.it).
 * Based on G. Degrassi, P. Slavich and F. Zwirner,
 * Nucl. Phys. B611 (2001) 403 [hep-ph/0105096].
 *
 * Last update:  13/12/2001.
 *
 * I/O PARAMETERS:
 * t = m_top^2, g = m_gluino^2, T1 = m_stop1^2, T2 = m_stop2^2,
 * st = sin(theta_stop), ct = cos(theta_stop), q = Q^2 (ren. scale),
 * mu = Higgs mixing parameter, tanb = tan(beta), vv = v^2,
 * OS = renormalization scheme for 1-loop (0 = DRbar, 1 = On-Shell),
 * Sij = 2-loop corrections to the CP-even Higgs mass matrix elements.
 */
int dszhiggs_(double * t, double * g, double * T1, double * T2,
              double * st, double * ct, double * q, double * mu,
              double * tanb, double *vv, double * gs, int * OS,
              double * S11, double * S22, double * S12);

/**
 * Two-loop O(a_t^2 + at ab + ab^2) corrections to the CP-even Higgs masses
 * Written by P. Slavich (e-mail: slavich@mppmu.mpg.de).
 * Based on A. Dedes, G. Degrassi and P. Slavich, hep-ph/0305127.
 *
 * Last update:  13/05/2003.
 *
 * I/O PARAMETERS:
 * t = m_top^2, b = m_bot^2, A0 = m_A^2, T1 = m_stop1^2, T2 = m_stop2^2,
 * B1 = m_sbot1^2, B2 = m_sbot2^2, st = sin(theta_stop),
 * ct = cos(theta_stop), sb = sin(theta_sbot), cb = cos(theta_sbot),
 * q = Q^2 (ren. scale), mu = Higgs mixing parameter, tanb = tan(beta),
 * vv = v^2,
 * Sij = 2-loop corrections to the CP-even Higgs mass matrix elements,
 * atat = if 1, returns only O(at^2) corrections
 *
 * Notice: we assume that the 1-loop part is computed in terms of
 *         running (DRbar) parameters, evaluated at the scale Q. The
 *         parameters in the bottom/sbottom sector should be computed
 *         in terms of the "resummed" bottom Yukawa coupling.
 */
int ddshiggs_(double * t, double * b, double * A0, double * T1,
              double * T2, double * B1, double * B2, double * st,
              double * ct, double * sb, double * cb, double * q,
              double * mu, double * tanb, double * vv,
              double * S11, double * S12, double * S22, int * atat);

/**
 * Two-loop O(a_tau^2) corrections to the CP-even Higgs mass matrix.
 * Routine written by P. Slavich (e-mail: slavich@mppmu.mpg.de).
 * Based on A. Brignole, G. Degrassi, P. Slavich and F. Zwirner,
 * hep-ph/0112177 with appropriate substitutions for top -> tau.
 *
 * Last update:  19/06/2003.
 *
 * I/O PARAMETERS:
 * t = m_tau^2, A0 = m_A^2, BL = m_snu^2, T1 = m_stau1^2, T2 = m_stau2^2,
 * st = sin(theta_stau), ct = cos(theta_stau), q = Q^2 (ren. scale),
 * mu = Higgs mixing parameter, tb = tan(beta), vv = v^2,
 * OS = renormalization scheme for 1-loop (0 = DRbar, 1 = On-Shell),
 * Sij = 2-loop corrections to the CP-even Higgs mass matrix elements.
 */
int tausqhiggs_(double * t, double * A0, double * BL, double * T1,
                double * T2, double * st, double * ct, double * q,
                double * mu, double * tanb, double * vv, int * OS,
                double * S11, double * S22, double * S12);

/**
 * Two-loop O(a_tau * a_bottom) corrections to the CP-even Higgs mass
 * matrix.  Routine written by P. Slavich (e-mail:
 * slavich@mppmu.mpg.de).  Based on Allanach et.al., JHEP 0409 (2004)
 * 044 [hep-ph/0406166].
 *
 * I/O PARAMETERS:
 * t = m_tau^2, b = m_b^2,
 * T1 = m_stau1^2, T2 = m_stau2^2, B1 = m_sbottom1^2, B2 = m_sbottom2^2,
 * st = sin(theta_stau), ct = cos(theta_stau),
 * sb = sin(theta_sbottom), cb = cos(theta_sbottom),
 * q = Q^2 (ren. scale), mu = Higgs mixing parameter, tb = tan(beta),
 * vv = v^2, Sij = 2-loop corrections to the CP-even Higgs mass matrix
 * elements.
 */
int taubot_(double * t, double * b,
            double * T1, double * T2, double * B1, double * B2,
            double * st, double * ct, double * sb, double * cb,
            double * q, double * mu, double * tanb, double * vv,
            double * S11, double * S22, double * S12);

/**
 * Two-loop O(a_t a_s) corrections to the CP-odd Higgs mass in the
 * DRbar scheme.  Written by P. Slavich (e-mail: slavich@pd.infn.it).
 * Based on G. Degrassi, P. Slavich and F. Zwirner,
 * Nucl. Phys. B611 (2001) 403 [hep-ph/0105096].
 *
 * Last update:  13/12/2001.
 *
 * I/O PARAMETERS:
 * t = m_top^2, g = m_gluino^2, T1 = m_stop1^2, T2 = m_stop2^2,
 * st = sin(theta_stop), ct = cos(theta_stop), q = Q^2 (ren. scale),
 * mu = Higgs mixing parameter, tanb = tan(beta), vv = v^2,
 * gs = strong coupling,
 * dma = 2-loop corrections to the CP-odd Higgs mass.
 */
int dszodd_(double * t, double * g, double * T1,
            double * T2, double * st, double * ct,
            double * q, double * mu, double * tanb,
            double * vv, double * gs, double * dma);

/**
 * Two-loop O(a_t^2 + at ab + ab^2) corrections to the CP-odd Higgs mass
 * Written by P. Slavich (e-mail: slavich@mppmu.mpg.de).
 * Based on A. Dedes, G. Degrassi and P. Slavich, hep-ph/0305127.
 *
 * Last update:  13/05/2003.
 *
 * I/O PARAMETERS:
 * t = m_top^2, b = m_bot^2, A0 = m_A^2, T1 = m_stop1^2, T2 = m_stop2^2,
 * B1 = m_sbot1^2, B2 = m_sbot2^2, st = sin(theta_stop),
 * ct = cos(theta_stop), sb = sin(theta_sbot), cb = cos(theta_sbot),
 * q = Q^2 (ren. scale), mu = Higgs mixing parameter, tanb = tan(beta),
 * vv = v^2,
 * dma = 2-loop corrections to the CP-odd Higgs mass.
 * atat = if 1, returns only O(at^2) corrections
 *
 * Notice: we assume that the 1-loop part is computed in terms of
 *         running (DRbar) parameters, evaluated at the scale Q. The
 *         parameters in the bottom/sbottom sector should be computed
 *         in term of the "resummed" bottom Yukawa coupling.
*/
int ddsodd_(double * t, double * b, double * A0, double * T1,
            double * T2, double * B1, double * B2, double * st,
            double * ct, double * sb, double * cb, double * q,
            double * mu, double * tanb, double * vv, double * dma, int * atat);

/**
 * Two-loop O(a_tau^2) corrections to the CP-odd Higgs mass in the
 * DRbar scheme. Written by P. Slavich (e-mail: slavich@mppmu.mpg.de).
 * Based on A. Brignole, G. Degrassi, P. Slavich and F. Zwirner,
 * hep-ph/0112177 with appropriate substitutions for top -> tau.
 *
 * Last update:  19/06/2003.
 *
 * I/O PARAMETERS:
 * t = m_tau^2, A0 = m_A^2, BL = m_snu^2, T1 = m_stau1^2, T2 = m_stau2^2,
 * st = sin(theta_stau), ct = cos(theta_stau), q = Q^2 (ren. scale),
 * mu = Higgs mixing parameter, tb = tan(beta), vv = v^2,
 * dma = 2-loop corrections to the CP-odd Higgs mass.
 */
int tausqodd_(double * t, double * A0, double * BL, double * T1,
              double * T2, double * st, double * ct, double * q,
              double * mu, double * tanb, double * vv,
              double * dma);

/**
 * Two-loop O(a_tau * a_bottom) corrections to the CP-odd Higgs mass
 * matrix.  Routine written by P. Slavich (e-mail:
 * slavich@mppmu.mpg.de).  Based on Allanach et.al., JHEP 0409 (2004)
 * 044 [hep-ph/0406166].
 *
 * I/O PARAMETERS:
 * t = m_tau^2, b = m_b^2,
 * T1 = m_stau1^2, T2 = m_stau2^2, B1 = m_sbottom1^2, B2 = m_sbottom2^2,
 * st = sin(theta_stau), ct = cos(theta_stau),
 * sb = sin(theta_sbottom), cb = cos(theta_sbottom),
 * q = Q^2 (ren. scale), mu = Higgs mixing parameter, tb = tan(beta),
 * vv = v^2, dma = 2-loop corrections to the CP-odd Higgs mass.
 */
int taubotodd_(double * t, double * b,
            double * T1, double * T2, double * B1, double * B2,
            double * st, double * ct, double * sb, double * cb,
            double * q, double * mu, double * tanb, double * vv,
            double * dma);

#ifdef __cplusplus
}
#endif