Skip to content

Commit

Permalink
Merge remote-tracking branch 'organization/development' into feature-…
Browse files Browse the repository at this point in the history 
…HSSUSY-3L
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Jun 12, 2018
2 parents f57b231 + 8ee3a1d commit bcc4087
Show file tree
Hide file tree
Showing 6 changed files with 48 additions and 31 deletions.
13 changes: 7 additions & 6 deletions doc/slha_input.dox
View file
Expand Up @@ -35,7 +35,7 @@ Block FlexibleSUSY
22 0 # EFT index of SM-like Higgs in the BSM model (0 = lightest Higgs)
23 1 # calculate BSM pole masses
24 123111321 # individual threshold correction loop orders
25 0 # ren. scheme for Higgs 3L corrections (0 = DR, 1 = MDR)
25 0 # ren. scheme for Higgs 3L corrections (0 = DR', 1 = MDR', 2 = H3m)
26 1 # Higgs 3-loop corrections O(alpha_t alpha_s^2)
27 1 # Higgs 3-loop corrections O(alpha_b alpha_s^2)
28 1 # Higgs 3-loop corrections O(alpha_t^2 alpha_s)
Expand Down Expand Up @@ -77,7 +77,7 @@ index | description | possible values
22 | EFT Higgs index | any integer >= 0 | 0 (= lightest)
23 | calculate pole masses of BSM particles | 0 (no) or 1 (yes) | 1 (= yes)
24 | individual threshold corrections | positive integer | 123111321
25 | ren. scheme for higgs 3L corrections | 0 (DR-bar) or 1 (MDR-bar) | 0 (= DR-bar)
25 | ren. scheme for higgs 3L corrections | 0 (DR'), 1 (MDR'), 2 (H3m) | 0 (= DR')
26 | higgs 3-loop correction O(at as^2) | 0, 1 | 1 (= enabled)
27 | higgs 3-loop correction O(ab as^2) | 0, 1 | 1 (= enabled)
28 | higgs 3-loop correction O(at^2 as) | 0, 1 | 1 (= enabled)
Expand Down Expand Up @@ -469,9 +469,10 @@ flag `FlexibleSUSY[26]` to `1`. To enable the 3-loop correction
`FlexibleSUSY[27]` to `1`.

The 3-loop corrections from [arxiv:hep-ph/1005.5709] can be calculated
in the \f$\overline{DR}\f$ or \f$\overline{MDR}\f$ scheme. To use the
\f$\overline{DR}\f$ scheme, set `FlexibleSUSY[25]` to `0`. To use the
\f$\overline{MDR}\f$ scheme, set `FlexibleSUSY[25]` to `1`.
in the \f$\overline{DR}'\f$, \f$\overline{MDR}'\f$ or H3m scheme. To
use the \f$\overline{DR}'\f$ scheme, set `FlexibleSUSY[25]` to `0`.
To use the \f$\overline{MDR}'\f$ scheme, set `FlexibleSUSY[25]` to
`1`. To use the H3m scheme, set `FlexibleSUSY[25]` to `2`.

We recommend to set the following model file options to enable the
3-loop Higgs pole mass corrections in the MSSM:
Expand Down Expand Up @@ -499,7 +500,7 @@ Block FlexibleSUSY
10 1 # Higgs 2-loop corrections O((alpha_t + alpha_b)^2)
11 1 # Higgs 2-loop corrections O(alpha_tau^2)
24 123111221 # individual threshold correction loop orders
25 0 # ren. scheme for Higgs 3L corrections (0 = DR, 1 = MDR)
25 0 # ren. scheme for Higgs 3L corrections (0 = DR', 1 = MDR', 2 = H3m)
26 1 # Higgs 3-loop corrections O(alpha_t alpha_s^2)
27 1 # Higgs 3-loop corrections O(alpha_b alpha_s^2)
28 1 # Higgs 3-loop corrections O(alpha_t^2 alpha_s)
Expand Down
54 changes: 35 additions & 19 deletions meta/SelfEnergies.m
View file
Expand Up @@ -1433,44 +1433,60 @@ therefore not be accessed in the form Glu(gO2).
#endif
try {
const auto mdrScheme = HIGGS_3LOOP_MDR_SCHEME;
const auto ren_scheme = HIGGS_3LOOP_SCHEME;
const bool verbose = false;
himalaya::HierarchyCalculator hc(pars, verbose);
if (HIGGS_3LOOP_CORRECTION_AT_AS_AS) {
const auto hier = hc.calculateDMh3L(false, mdrScheme);
VERBOSE_MSG(\"Himalaya top (hierarchy, uncertainties) = (\"
<< hier.getSuitableHierarchy() << \", {\"
<< hier.getExpUncertainty(1) << \", \"
<< hier.getExpUncertainty(2) << \", \"
<< hier.getExpUncertainty(3) << \"})\");
#if Himalaya_VERSION_MAJOR < 2
const auto hier = hc.calculateDMh3L(false, ren_scheme);
#else
const auto hier = hc.calculateDMh3L(false);
#endif
// calculate the 3-loop corrections
self_energy_3l += - hier.getDMh(3);
if (mdrScheme) {
// calculate the 1- and 2-loop shift DR -> MDR
#if Himalaya_VERSION_MAJOR < 2
if (ren_scheme) {
// calculate shift DR -> MDR
self_energy_3l += - hier.getDRToMDRShift();
}
#else
if (ren_scheme == 1) {
// calculate shift DR' -> MDR'
self_energy_3l += - hier.getDRbarPrimeToMDRbarPrimeShift();
} else if (ren_scheme == 1) {
// calculate shift DR' -> H3m
self_energy_3l += - hier.getDRbarPrimeToH3mShift();
}
#endif
}
if (HIGGS_3LOOP_CORRECTION_AB_AS_AS) {
const auto hier = hc.calculateDMh3L(true, mdrScheme);
VERBOSE_MSG(\"Himalaya bottom (hierarchy, uncertainties) = (\"
<< hier.getSuitableHierarchy() << \", {\"
<< hier.getExpUncertainty(1) << \", \"
<< hier.getExpUncertainty(2) << \", \"
<< hier.getExpUncertainty(3) << \"})\");
#if Himalaya_VERSION_MAJOR < 2
const auto hier = hc.calculateDMh3L(true, ren_scheme);
#else
const auto hier = hc.calculateDMh3L(false);
#endif
// calculate the 3-loop corrections
self_energy_3l += - hier.getDMh(3);
if (mdrScheme) {
// calculate the 1- and 2-loop shift DR -> MDR
#if Himalaya_VERSION_MAJOR < 2
if (ren_scheme) {
// calculate the shift DR -> MDR
self_energy_3l += - hier.getDRToMDRShift();
}
#else
if (ren_scheme == 1) {
// calculate shift DR' -> MDR'
self_energy_3l += - hier.getDRbarPrimeToMDRbarPrimeShift();
} else if (ren_scheme == 1) {
// calculate shift DR' -> H3m
self_energy_3l += - hier.getDRbarPrimeToH3mShift();
}
#endif
}
} catch (const std::exception& e) {
VERBOSE_MSG(e.what());
Expand Down
2 changes: 1 addition & 1 deletion src/loop_corrections.hpp
View file
Expand Up @@ -39,7 +39,7 @@ struct Loop_corrections {
bool higgs_at_at_as{true};
bool higgs_at_at_at{true};
bool higgs_at_as_as_as{true};
int higgs_3L_mdr_scheme{1};
int higgs_3L_scheme{1};
int top_qcd{1}; ///< top pole mass QCD corrections
};

Expand Down
6 changes: 3 additions & 3 deletions src/spectrum_generator_settings.cpp
View file
Expand Up @@ -125,7 +125,7 @@ void Spectrum_generator_settings::set(const Spectrum_generator_settings::Setting
* | eft_higgs_index | any integer >= 0 | 0 (= lightest) |
* | calculate_bsm_masses | 0 (no) or 1 (yes) | 1 (= yes) |
* | threshold_corrections | positive integer | 123111321 |
* | higgs_3loop_ren_scheme_atb_as2 | 0 (DR-bar) or 1 (MDR-bar) | 0 (= MDR-bar) |
* | higgs_3loop_ren_scheme_atb_as2 | 0 (DR'), 1 (MDR'), 2 (H3m) | 0 (= DR') |
* | higgs_3loop_correction_at_as2 | 0, 1 | 1 (= enabled) |
* | higgs_3loop_correction_ab_as2 | 0, 1 | 1 (= enabled) |
* | higgs_3loop_correction_at2_as | 0, 1 | 1 (= enabled) |
Expand Down Expand Up @@ -180,7 +180,7 @@ Loop_corrections Spectrum_generator_settings::get_loop_corrections() const
loop_corrections.higgs_at_at_as = get(higgs_3loop_correction_at2_as);
loop_corrections.higgs_at_at_at = get(higgs_3loop_correction_at3);
loop_corrections.higgs_at_as_as_as = get(higgs_4loop_correction_at_as3);
loop_corrections.higgs_3L_mdr_scheme = get(higgs_3loop_ren_scheme_atb_as2);
loop_corrections.higgs_3L_scheme = get(higgs_3loop_ren_scheme_atb_as2);
loop_corrections.top_qcd = get(top_pole_qcd_corrections);

return loop_corrections;
Expand All @@ -198,7 +198,7 @@ void Spectrum_generator_settings::set_loop_corrections(
set(higgs_3loop_correction_at2_as, loop_corrections.higgs_at_at_as);
set(higgs_3loop_correction_at3 , loop_corrections.higgs_at_at_at);
set(higgs_4loop_correction_at_as3, loop_corrections.higgs_at_as_as_as);
set(higgs_3loop_ren_scheme_atb_as2, loop_corrections.higgs_3L_mdr_scheme);
set(higgs_3loop_ren_scheme_atb_as2, loop_corrections.higgs_3L_scheme);
set(top_pole_qcd_corrections, loop_corrections.top_qcd);
}

Expand Down
2 changes: 1 addition & 1 deletion templates/LesHouches.in
View file
Expand Up @@ -26,7 +26,7 @@ Block FlexibleSUSY
22 0 # EFT index of SM-like Higgs in the BSM model
23 1 # calculate BSM pole masses
24 123111321 # individual threshold correction loop orders
25 0 # ren. scheme for Higgs 3L corrections (0 = DR, 1 = MDR)
25 0 # ren. scheme for Higgs 3L corrections (0 = DR', 1 = MDR', 2 = H3m)
26 1 # Higgs 3-loop corrections O(alpha_t alpha_s^2)
27 1 # Higgs 3-loop corrections O(alpha_b alpha_s^2)
28 1 # Higgs 3-loop corrections O(alpha_t^2 alpha_s)
Expand Down
2 changes: 1 addition & 1 deletion templates/mass_eigenstates.cpp.in
View file
Expand Up @@ -76,7 +76,7 @@ namespace flexiblesusy {
#define TOP_POLE_QCD_CORRECTION loop_corrections.top_qcd
#define HIGGS_3LOOP_CORRECTION_AT_AS_AS loop_corrections.higgs_at_as_as
#define HIGGS_3LOOP_CORRECTION_AB_AS_AS loop_corrections.higgs_ab_as_as
#define HIGGS_3LOOP_MDR_SCHEME loop_corrections.higgs_3L_mdr_scheme
#define HIGGS_3LOOP_SCHEME loop_corrections.higgs_3L_scheme
#define HIGGS_3LOOP_CORRECTION_AT_AT_AS loop_corrections.higgs_at_at_as
#define HIGGS_3LOOP_CORRECTION_AT_AT_AT loop_corrections.higgs_at_at_at
#define HIGGS_4LOOP_CORRECTION_AT_AS_AS_AS loop_corrections.higgs_at_as_as_as
Expand Down

0 comments on commit bcc4087

Please sign in to comment.