Hyperiso/MesonMixingWilson_8cpp_source.html

#include "MesonMixingWilson.h"


C_mix_bd_1::C_mix_bd_1() : WilsonCoefficient("C_BD_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = {

        {

            {"WPARAM_MATCH_SM", LhaID(2, 1)},  // x_t

            {ParameterType::SM, "SMINPUTS", 2},                   // G_F

            {ParameterType::SM, "MASS", 24},                      // M_W

            {ParameterType::SM, "VCKM", LhaID(2, 0)},             // V_td

            {ParameterType::SM, "VCKM", LhaID(2, 2)},             // V_tb

        },

        compute_LO,

        get_lhaid_from_name(QCDOrder::LO)

    };

}


complex_t C_mix_bd_1::compute_LO(const ParamSrc& src) {

    double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});;

    double G_F = src.get_val(ParameterType::SM, "SMINPUTS", 2);

    double M_W = src.get_val(ParameterType::SM, "MASS", 24);

    complex_t V_td = src.get_val(ParameterType::SM, "VCKM", {2, 0});

    complex_t V_tb = src.get_val(ParameterType::SM, "VCKM", {2, 2});

    return pow(G_F * M_W * std::conj(V_td) * V_tb / (2 * PI), 2) * S0(xt);

}


C_mix_bd_1_tilde::C_mix_bd_1_tilde(): WilsonCoefficient("CT_BD_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_2::C_mix_bd_2(): WilsonCoefficient("C_BD_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_2_tilde::C_mix_bd_2_tilde(): WilsonCoefficient("CT_BD_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_3::C_mix_bd_3(): WilsonCoefficient("C_BD_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_3_tilde::C_mix_bd_3_tilde(): WilsonCoefficient("CT_BD_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_4::C_mix_bd_4(): WilsonCoefficient("C_BD_4", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bd_5::C_mix_bd_5(): WilsonCoefficient("C_BD_5", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_1::C_mix_bs_1() : WilsonCoefficient("C_BS_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = {

        {

            {"WPARAM_MATCH_SM", LhaID(2, 1)},  // x_t

            {ParameterType::SM, "SMINPUTS", 2},                   // G_F

            {ParameterType::SM, "MASS", 24},                      // M_W

            {ParameterType::SM, "VCKM", LhaID(2, 1)},             // V_ts

            {ParameterType::SM, "VCKM", LhaID(2, 2)},             // V_tb

        },

        compute_LO,

        get_lhaid_from_name(QCDOrder::LO)

    };

}


complex_t C_mix_bs_1::compute_LO(const ParamSrc& src) {

    double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});;

    double G_F = src.get_val(ParameterType::SM, "SMINPUTS", 2);

    double M_W = src.get_val(ParameterType::SM, "MASS", 24);

    complex_t V_ts = src.get_val(ParameterType::SM, "VCKM", {2, 1});

    complex_t V_tb = src.get_val(ParameterType::SM, "VCKM", {2, 2});

    return pow(G_F * M_W * std::conj(V_ts) * V_tb / (2 * PI), 2) * S0(xt);

}


C_mix_bs_1_tilde::C_mix_bs_1_tilde() : WilsonCoefficient("CT_BS_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_2::C_mix_bs_2() : WilsonCoefficient("C_BS_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_2_tilde::C_mix_bs_2_tilde() : WilsonCoefficient("CT_BS_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_3::C_mix_bs_3() : WilsonCoefficient("C_BS_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_3_tilde::C_mix_bs_3_tilde() : WilsonCoefficient("CT_BS_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_4::C_mix_bs_4() : WilsonCoefficient("C_BS_4", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_bs_5::C_mix_bs_5() : WilsonCoefficient("C_BS_5", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_1::C_mix_sd_1() : WilsonCoefficient("C_SD_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = {

        {

            {"WPARAM_MATCH_SM", LhaID(2, 1)},  // x_t

            {ParameterType::SM, "SMINPUTS", 2},                   // G_F

            {ParameterType::SM, "MASS", 24},                      // M_W

            {ParameterType::SM, "VCKM", LhaID(2, 0)},             // V_td

            {ParameterType::SM, "VCKM", LhaID(2, 1)},             // V_ts

        },

        compute_LO,

        get_lhaid_from_name(QCDOrder::LO)

    };

}


complex_t C_mix_sd_1::compute_LO(const ParamSrc& src) {

    double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", {2, 1});;

    double G_F = src.get_val(ParameterType::SM, "SMINPUTS", 2);

    double M_W = src.get_val(ParameterType::SM, "MASS", 24);

    complex_t V_td = src.get_val(ParameterType::SM, "VCKM", {2, 0});

    complex_t V_ts = src.get_val(ParameterType::SM, "VCKM", {2, 1});

    return pow(G_F * M_W * abs(std::conj(V_td) * V_ts) / (2 * PI), 2) * S0(xt);

}


C_mix_sd_1_tilde::C_mix_sd_1_tilde() : WilsonCoefficient("CT_SD_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_2::C_mix_sd_2() : WilsonCoefficient("C_SD_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_2_tilde::C_mix_sd_2_tilde() : WilsonCoefficient("CT_SD_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_3::C_mix_sd_3() : WilsonCoefficient("C_SD_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_3_tilde::C_mix_sd_3_tilde() : WilsonCoefficient("CT_SD_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_4::C_mix_sd_4() : WilsonCoefficient("C_SD_4", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_sd_5::C_mix_sd_5() : WilsonCoefficient("C_SD_5", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_1::C_mix_cu_1() : WilsonCoefficient("C_CU_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_1_tilde::C_mix_cu_1_tilde() : WilsonCoefficient("CT_CU_1", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_2::C_mix_cu_2() : WilsonCoefficient("C_CU_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_2_tilde::C_mix_cu_2_tilde() : WilsonCoefficient("CT_CU_2", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_3::C_mix_cu_3() : WilsonCoefficient("C_CU_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_3_tilde::C_mix_cu_3_tilde() : WilsonCoefficient("CT_CU_3", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_4::C_mix_cu_4() : WilsonCoefficient("C_CU_4", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


C_mix_cu_5::C_mix_cu_5() : WilsonCoefficient("C_CU_5", GroupMapper::str(WGroup::MESON_MIXING, ScaleType::MATCHING)) {

    matching_info[QCDOrder::LO] = MatchingInfo(get_lhaid_from_name(QCDOrder::LO));

}


QCDOrder::LO
@ LO

ParameterType::SM
@ SM

ParameterType::WILSON
@ WILSON

ScaleType
ScaleType
Definition GeneralEnum.h:666

ScaleType::MATCHING
@ MATCHING

WGroup
WGroup
Definition GeneralEnum.h:618

WGroup::MESON_MIXING
@ MESON_MIXING

MesonMixingWilson.h

complex_t
std::complex< double > complex_t
Convenience alias for std::complex<double>.
Definition Utils.h:35

C_mix_bd_1_tilde::C_mix_bd_1_tilde
C_mix_bd_1_tilde()
Definition MesonMixingWilson.cpp:26

C_mix_bd_1::C_mix_bd_1
C_mix_bd_1()
Definition MesonMixingWilson.cpp:3

C_mix_bd_1::compute_LO
static complex_t compute_LO(const ParamSrc &src)
Definition MesonMixingWilson.cpp:17

C_mix_bd_2_tilde::C_mix_bd_2_tilde
C_mix_bd_2_tilde()
Definition MesonMixingWilson.cpp:34

C_mix_bd_2::C_mix_bd_2
C_mix_bd_2()
Definition MesonMixingWilson.cpp:30

C_mix_bd_3_tilde::C_mix_bd_3_tilde
C_mix_bd_3_tilde()
Definition MesonMixingWilson.cpp:42

C_mix_bd_3::C_mix_bd_3
C_mix_bd_3()
Definition MesonMixingWilson.cpp:38

C_mix_bd_4::C_mix_bd_4
C_mix_bd_4()
Definition MesonMixingWilson.cpp:46

C_mix_bd_5::C_mix_bd_5
C_mix_bd_5()
Definition MesonMixingWilson.cpp:50

C_mix_bs_1_tilde::C_mix_bs_1_tilde
C_mix_bs_1_tilde()
Definition MesonMixingWilson.cpp:77

C_mix_bs_1::C_mix_bs_1
C_mix_bs_1()
Definition MesonMixingWilson.cpp:54

C_mix_bs_1::compute_LO
static complex_t compute_LO(const ParamSrc &src)
Definition MesonMixingWilson.cpp:68

C_mix_bs_2_tilde::C_mix_bs_2_tilde
C_mix_bs_2_tilde()
Definition MesonMixingWilson.cpp:85

C_mix_bs_2::C_mix_bs_2
C_mix_bs_2()
Definition MesonMixingWilson.cpp:81

C_mix_bs_3_tilde::C_mix_bs_3_tilde
C_mix_bs_3_tilde()
Definition MesonMixingWilson.cpp:93

C_mix_bs_3::C_mix_bs_3
C_mix_bs_3()
Definition MesonMixingWilson.cpp:89

C_mix_bs_4::C_mix_bs_4
C_mix_bs_4()
Definition MesonMixingWilson.cpp:97

C_mix_bs_5::C_mix_bs_5
C_mix_bs_5()
Definition MesonMixingWilson.cpp:101

C_mix_cu_1_tilde::C_mix_cu_1_tilde
C_mix_cu_1_tilde()
Definition MesonMixingWilson.cpp:160

C_mix_cu_1::C_mix_cu_1
C_mix_cu_1()
Definition MesonMixingWilson.cpp:156

C_mix_cu_2_tilde::C_mix_cu_2_tilde
C_mix_cu_2_tilde()
Definition MesonMixingWilson.cpp:168

C_mix_cu_2::C_mix_cu_2
C_mix_cu_2()
Definition MesonMixingWilson.cpp:164

C_mix_cu_3_tilde::C_mix_cu_3_tilde
C_mix_cu_3_tilde()
Definition MesonMixingWilson.cpp:176

C_mix_cu_3::C_mix_cu_3
C_mix_cu_3()
Definition MesonMixingWilson.cpp:172

C_mix_cu_4::C_mix_cu_4
C_mix_cu_4()
Definition MesonMixingWilson.cpp:180

C_mix_cu_5::C_mix_cu_5
C_mix_cu_5()
Definition MesonMixingWilson.cpp:184

C_mix_sd_1_tilde::C_mix_sd_1_tilde
C_mix_sd_1_tilde()
Definition MesonMixingWilson.cpp:128

C_mix_sd_1::compute_LO
static complex_t compute_LO(const ParamSrc &src)
Definition MesonMixingWilson.cpp:119

C_mix_sd_1::C_mix_sd_1
C_mix_sd_1()
Definition MesonMixingWilson.cpp:105

C_mix_sd_2_tilde::C_mix_sd_2_tilde
C_mix_sd_2_tilde()
Definition MesonMixingWilson.cpp:136

C_mix_sd_2::C_mix_sd_2
C_mix_sd_2()
Definition MesonMixingWilson.cpp:132

C_mix_sd_3_tilde::C_mix_sd_3_tilde
C_mix_sd_3_tilde()
Definition MesonMixingWilson.cpp:144

C_mix_sd_3::C_mix_sd_3
C_mix_sd_3()
Definition MesonMixingWilson.cpp:140

C_mix_sd_4::C_mix_sd_4
C_mix_sd_4()
Definition MesonMixingWilson.cpp:148

C_mix_sd_5::C_mix_sd_5
C_mix_sd_5()
Definition MesonMixingWilson.cpp:152

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

PI
constexpr double PI
Definition constants.h:7

pow
scalar_t pow(const scalar_t &base, const scalar_t &exp)
Definition scalar.cpp:75

S0
double S0(double x)
Wilson special function S0 depending on xt.
Definition special_SM.cpp:1159

LhaID
Represents an identifier of a LHA element, possibly containing several sub-ids.
Definition LhaID.h:56

MatchingInfo
Matching metadata for a Wilson coefficient at a specific QCD order.
Definition Wilson.h:62