Hyperiso 1.0.3
Modular flavour-physics calculations, Wilson coefficients and statistical inference
Loading...
Searching...
No Matches
DChargedCurrentWilsonTHDM.cpp
Go to the documentation of this file.
1#include "DChargedCurrentWilsonTHDM.h"
2
3C_S1_cs_THDM::C_S1_cs_THDM()
4 : WilsonCoefficient("C_S1_cs_THDM", GroupMapper::str(WGroup::CC_cs, ScaleType::MATCHING))
5{
6 matching_info[QCDOrder::LO] = {
7 {
8 {"WPARAM_SI_BSM", 8}, // ld
9 {ParameterType::SM, "MASS", 15}, // m_tau
10 {ParameterType::BSM, "MASS", 37}, // m_H
11 {ParameterType::BSM, "YE", {3, 3}},// l_tau
12 {ParameterType::SM, "QCD", {5, 1}} // m_b
13 },
14 compute_LO,
15 get_lhaid_from_name(QCDOrder::LO)
16 };
17}
18
19double C_S1_cs_THDM::compute_LO(const ParamSrc& src) {
20 double ld = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", 8);
21 double mH = src.get_val(ParameterType::BSM, "MASS", 37);
22 double m_tau = src.get_val(ParameterType::SM, "MASS", 15);
23 double l_tau = src.get_val(ParameterType::BSM, "YE", {3, 3});
24 double m_b = src.get_val(ParameterType::SM, "QCD", {5, 1});
25
26 return -m_b * m_tau * ld * l_tau / std::pow(mH, 2);
27}
28
29C_S2_cs_THDM::C_S2_cs_THDM()
30 : WilsonCoefficient("C_S2_cs_THDM", GroupMapper::str(WGroup::CC_cs, ScaleType::MATCHING))
31{
32 matching_info[QCDOrder::LO] = {
33 {
34 {"WPARAM_SI_BSM", 7}, // lu
35 {ParameterType::SM, "MASS", 15}, // m_tau
36 {ParameterType::SM, "MASS", 4}, // m_c (charm)
37 {ParameterType::BSM, "MASS", 37}, // m_H
38 {ParameterType::BSM, "YE", {3, 3}} // l_tau
39 },
40 compute_LO,
41 get_lhaid_from_name(QCDOrder::LO)
42 };
43}
44
45double C_S2_cs_THDM::compute_LO(const ParamSrc& src) {
46 double lu = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", 7);
47 double mH = src.get_val(ParameterType::BSM, "MASS", 37);
48 double m_c = src.get_val(ParameterType::SM, "MASS", 4); // m_c
49 double m_tau = src.get_val(ParameterType::SM, "MASS", 15); // m_tau
50 double l_tau = src.get_val(ParameterType::BSM, "YE", {3, 3}); // l_tau
51
52 return -m_c * m_tau * lu * l_tau / std::pow(mH, 2);
53}
54
55
56
57C_S1_cd_THDM::C_S1_cd_THDM()
58 : WilsonCoefficient("C_S1_cd_THDM", GroupMapper::str(WGroup::CC_cd, ScaleType::MATCHING))
59{
60 matching_info[QCDOrder::LO] = {
61 {
62 {"WPARAM_SI_BSM", 8}, // ld
63 {ParameterType::SM, "MASS", 15}, // m_tau
64 {ParameterType::BSM, "MASS", 37}, // m_H
65 {ParameterType::BSM, "YE", {3, 3}},// l_tau
66 {ParameterType::SM, "QCD", {5, 1}} // m_b
67 },
68 compute_LO,
69 get_lhaid_from_name(QCDOrder::LO)
70 };
71}
72
73double C_S1_cd_THDM::compute_LO(const ParamSrc& src) {
74 double ld = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", 8);
75 double mH = src.get_val(ParameterType::BSM, "MASS", 37);
76 double m_tau = src.get_val(ParameterType::SM, "MASS", 15);
77 double l_tau = src.get_val(ParameterType::BSM, "YE", {3, 3});
78 double m_b = src.get_val(ParameterType::SM, "QCD", {5, 1});
79
80 return -m_b * m_tau * ld * l_tau / std::pow(mH, 2);
81}
82
83C_S2_cd_THDM::C_S2_cd_THDM()
84 : WilsonCoefficient("C_S2_cd_THDM", GroupMapper::str(WGroup::CC_cd, ScaleType::MATCHING))
85{
86 matching_info[QCDOrder::LO] = {
87 {
88 {"WPARAM_SI_BSM", 7}, // lu
89 {ParameterType::SM, "MASS", 15}, // m_tau
90 {ParameterType::SM, "MASS", 4}, // m_c (charm)
91 {ParameterType::BSM, "MASS", 37}, // m_H
92 {ParameterType::BSM, "YE", {3, 3}} // l_tau
93 },
94 compute_LO,
95 get_lhaid_from_name(QCDOrder::LO)
96 };
97}
98
99double C_S2_cd_THDM::compute_LO(const ParamSrc& src) {
100 double lu = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", 7);
101 double mH = src.get_val(ParameterType::BSM, "MASS", 37);
102 double m_c = src.get_val(ParameterType::SM, "MASS", 4); // m_c
103 double m_tau = src.get_val(ParameterType::SM, "MASS", 15); // m_tau
104 double l_tau = src.get_val(ParameterType::BSM, "YE", {3, 3}); // l_tau
105
106 return -m_c * m_tau * lu * l_tau / std::pow(mH, 2);
107}
ScaleType
ScaleType
Definition GeneralEnum.h:666
WGroup
WGroup
Definition GeneralEnum.h:618
C_S1_cd_THDM::compute_LO
static double compute_LO(const ParamSrc &src)
Definition DChargedCurrentWilsonTHDM.cpp:73
C_S1_cd_THDM::C_S1_cd_THDM
C_S1_cd_THDM()
Definition DChargedCurrentWilsonTHDM.cpp:57
C_S1_cs_THDM::compute_LO
static double compute_LO(const ParamSrc &src)
Definition DChargedCurrentWilsonTHDM.cpp:19
C_S1_cs_THDM::C_S1_cs_THDM
C_S1_cs_THDM()
Definition DChargedCurrentWilsonTHDM.cpp:3
C_S2_cd_THDM::C_S2_cd_THDM
C_S2_cd_THDM()
Definition DChargedCurrentWilsonTHDM.cpp:83
C_S2_cd_THDM::compute_LO
static double compute_LO(const ParamSrc &src)
Definition DChargedCurrentWilsonTHDM.cpp:99
C_S2_cs_THDM::C_S2_cs_THDM
C_S2_cs_THDM()
Definition DChargedCurrentWilsonTHDM.cpp:29
C_S2_cs_THDM::compute_LO
static double compute_LO(const ParamSrc &src)
Definition DChargedCurrentWilsonTHDM.cpp:45
GroupMapper
Mapper for WGroup <-> WGroupId <-> optional external string.
Definition wgroup_ids.hpp:65
ParamSrc
Lightweight view over a set of source parameters keyed by ParamId.
Definition SourcesView.h:249
ParamSrc::get_val
scalar_t get_val(const ParamId &id) const
Retrieves the current value of a parameter.
Definition SourceView.cpp:111
WilsonCoefficient
Abstract base class representing a Wilson coefficient and its matching information.
Definition Wilson.h:153
WilsonCoefficient::get_lhaid_from_name
LhaID get_lhaid_from_name(QCDOrder order)
Computes the LhaID directly from the coefficient base name and mapping conventions.
Definition Wilson.cpp:61
WilsonCoefficient::matching_info
std::map< QCDOrder, MatchingInfo > matching_info
Matching metadata indexed by QCD order.
Definition Wilson.h:330