Hyperiso 1.0.3
Modular flavour-physics calculations, Wilson coefficients and statistical inference
Loading...
Searching...
No Matches
main_mixing.cpp
Go to the documentation of this file.
1#include "HyperisoMaster.h"
2#include "WilsonInterface.h"
3#include "ParameterSetter.h"
4#include "BlockProxy.h"
5#include <iostream>
6
7int main(){
8 HyperisoMaster hyp = HyperisoMaster();
9 HyperisoConfig config;
10 config.model = Model::SM;
11
12 hyp.init("lha/testInput.flha", config);
13 LOG_INFO("HyperisoMaster initialized");
14
15 auto wi = WilsonInterface(); // Initialize interface and build the required groups
16 LOG_INFO("WilsonInterface created");
17
18 WilsonBuildConfig config_k{{WGroup::MESON_MIXING}, 81, 4.18, QCDOrder::LO};
19
20 std::cout << "here" << std::endl;
21 wi.build(config_k);
22
23
24
25 LOG_INFO("Interface built");
26
27 BlockProxy().log_all_blocks(ParameterType::WILSON);
28
29 LOG_INFO("C_BD_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BD_1, QCDOrder::LO, ContributionType::SM));
30 LOG_INFO("CT_BD_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BD_1, QCDOrder::LO, ContributionType::SM));
31 LOG_INFO("C_BD_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BD_2, QCDOrder::LO, ContributionType::SM));
32 LOG_INFO("CT_BD_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BD_2, QCDOrder::LO, ContributionType::SM));
33 LOG_INFO("C_BD_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BD_3, QCDOrder::LO, ContributionType::SM));
34 LOG_INFO("CT_BD_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BD_3, QCDOrder::LO, ContributionType::SM));
35 LOG_INFO("C_BD_4(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BD_4, QCDOrder::LO, ContributionType::SM));
36 LOG_INFO("C_BD_5(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BD_5, QCDOrder::LO, ContributionType::SM));
37
38 LOG_INFO("C_BS_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BS_1, QCDOrder::LO, ContributionType::SM));
39 LOG_INFO("CT_BS_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BS_1, QCDOrder::LO, ContributionType::SM));
40 LOG_INFO("C_BS_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BS_2, QCDOrder::LO, ContributionType::SM));
41 LOG_INFO("CT_BS_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BS_2, QCDOrder::LO, ContributionType::SM));
42 LOG_INFO("C_BS_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BS_3, QCDOrder::LO, ContributionType::SM));
43 LOG_INFO("CT_BS_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_BS_3, QCDOrder::LO, ContributionType::SM));
44 LOG_INFO("C_BS_4(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BS_4, QCDOrder::LO, ContributionType::SM));
45 LOG_INFO("C_BS_5(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_BS_5, QCDOrder::LO, ContributionType::SM));
46
47 LOG_INFO("C_SD_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_SD_1, QCDOrder::LO, ContributionType::SM));
48 LOG_INFO("CT_SD_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_SD_1, QCDOrder::LO, ContributionType::SM));
49 LOG_INFO("C_SD_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_SD_2, QCDOrder::LO, ContributionType::SM));
50 LOG_INFO("CT_SD_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_SD_2, QCDOrder::LO, ContributionType::SM));
51 LOG_INFO("C_SD_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_SD_3, QCDOrder::LO, ContributionType::SM));
52 LOG_INFO("CT_SD_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_SD_3, QCDOrder::LO, ContributionType::SM));
53 LOG_INFO("C_SD_4(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_SD_4, QCDOrder::LO, ContributionType::SM));
54 LOG_INFO("C_SD_5(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_SD_5, QCDOrder::LO, ContributionType::SM));
55
56 LOG_INFO("C_CU_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_CU_1, QCDOrder::LO, ContributionType::SM));
57 LOG_INFO("CT_CU_1(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_CU_1, QCDOrder::LO, ContributionType::SM));
58 LOG_INFO("C_CU_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_CU_2, QCDOrder::LO, ContributionType::SM));
59 LOG_INFO("CT_CU_2(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_CU_2, QCDOrder::LO, ContributionType::SM));
60 LOG_INFO("C_CU_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_CU_3, QCDOrder::LO, ContributionType::SM));
61 LOG_INFO("CT_CU_3(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::CT_CU_3, QCDOrder::LO, ContributionType::SM));
62 LOG_INFO("C_CU_4(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_CU_4, QCDOrder::LO, ContributionType::SM));
63 LOG_INFO("C_CU_5(mu_W) at LO =", wi.getM(WGroup::MESON_MIXING, WCoef::C_CU_5, QCDOrder::LO, ContributionType::SM));
64
65
66
67 LOG_INFO("C_BD_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BD_1, QCDOrder::LO, ContributionType::SM));
68 LOG_INFO("CT_BD_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BD_1, QCDOrder::LO, ContributionType::SM));
69 LOG_INFO("C_BD_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BD_2, QCDOrder::LO, ContributionType::SM));
70 LOG_INFO("CT_BD_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BD_2, QCDOrder::LO, ContributionType::SM));
71 LOG_INFO("C_BD_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BD_3, QCDOrder::LO, ContributionType::SM));
72 LOG_INFO("CT_BD_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BD_3, QCDOrder::LO, ContributionType::SM));
73 LOG_INFO("C_BD_4(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BD_4, QCDOrder::LO, ContributionType::SM));
74 LOG_INFO("C_BD_5(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BD_5, QCDOrder::LO, ContributionType::SM));
75
76 LOG_INFO("C_BS_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BS_1, QCDOrder::LO, ContributionType::SM));
77 LOG_INFO("CT_BS_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BS_1, QCDOrder::LO, ContributionType::SM));
78 LOG_INFO("C_BS_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BS_2, QCDOrder::LO, ContributionType::SM));
79 LOG_INFO("CT_BS_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BS_2, QCDOrder::LO, ContributionType::SM));
80 LOG_INFO("C_BS_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BS_3, QCDOrder::LO, ContributionType::SM));
81 LOG_INFO("CT_BS_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_BS_3, QCDOrder::LO, ContributionType::SM));
82 LOG_INFO("C_BS_4(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BS_4, QCDOrder::LO, ContributionType::SM));
83 LOG_INFO("C_BS_5(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_BS_5, QCDOrder::LO, ContributionType::SM));
84
85 LOG_INFO("C_SD_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_SD_1, QCDOrder::LO, ContributionType::SM));
86 LOG_INFO("CT_SD_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_SD_1, QCDOrder::LO, ContributionType::SM));
87 LOG_INFO("C_SD_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_SD_2, QCDOrder::LO, ContributionType::SM));
88 LOG_INFO("CT_SD_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_SD_2, QCDOrder::LO, ContributionType::SM));
89 LOG_INFO("C_SD_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_SD_3, QCDOrder::LO, ContributionType::SM));
90 LOG_INFO("CT_SD_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_SD_3, QCDOrder::LO, ContributionType::SM));
91 LOG_INFO("C_SD_4(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_SD_4, QCDOrder::LO, ContributionType::SM));
92 LOG_INFO("C_SD_5(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_SD_5, QCDOrder::LO, ContributionType::SM));
93
94 LOG_INFO("C_CU_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_CU_1, QCDOrder::LO, ContributionType::SM));
95 LOG_INFO("CT_CU_1(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_CU_1, QCDOrder::LO, ContributionType::SM));
96 LOG_INFO("C_CU_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_CU_2, QCDOrder::LO, ContributionType::SM));
97 LOG_INFO("CT_CU_2(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_CU_2, QCDOrder::LO, ContributionType::SM));
98 LOG_INFO("C_CU_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_CU_3, QCDOrder::LO, ContributionType::SM));
99 LOG_INFO("CT_CU_3(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::CT_CU_3, QCDOrder::LO, ContributionType::SM));
100 LOG_INFO("C_CU_4(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_CU_4, QCDOrder::LO, ContributionType::SM));
101 LOG_INFO("C_CU_5(mu_h) at LO =", wi.getR(WGroup::MESON_MIXING, WCoef::C_CU_5, QCDOrder::LO, ContributionType::SM));
102
103 // LOG_INFO("CQ1(mu_h) at LO =", wi.getR(WGroup::BScalar, WCoef::CQ1_MU, QCDOrder::LO));
104 // LOG_INFO("CQ2(mu_h) at LO =", wi.getR(WGroup::BScalar, WCoef::CQ2_MU, QCDOrder::LO));
105 // LOG_INFO("C10(mu_h) at LO =", wi.getR(WGroup::B, WCoef::C10, QCDOrder::LO));
106
107 // LOG_INFO("C7(mu_h) at LO =", wi.getR(WGroup::B, WCoef::C7, QCDOrder::LO));
108 // LOG_INFO("C7(mu_h) at NLO =", wi.getR(WGroup::B, WCoef::C7, QCDOrder::NLO));
109 // LOG_INFO("C7(mu_h) at NNLO =", wi.getR(WGroup::B, WCoef::C7, QCDOrder::NNLO));
110
111 // LOG_INFO("C7(mu_h) full =", wi.getFR(WGroup::B, WCoef::C7, QCDOrder::NNLO));
112
113 // LOG_INFO("C7(mu_W) at LO =", wi.getM(WGroup::B, WCoef::C7, QCDOrder::LO));
114 // LOG_INFO("C7(mu_W) at NLO =", wi.getM(WGroup::B, WCoef::C7, QCDOrder::NLO));
115 // LOG_INFO("C7(mu_W) at NNLO =", wi.getM(WGroup::B, WCoef::C7, QCDOrder::NNLO));
116
117 // LOG_INFO("C7(mu_W) full =", wi.getFM(WGroup::B, WCoef::C7, QCDOrder::NNLO));
118
119 return 0;
120}
BlockProxy.h
Facade combining block existence/logging (BlockProvider) and block listing (APIAdapter).
WGroup::MESON_MIXING
@ MESON_MIXING
HyperisoMaster.h
High-level helpers for initializing and monitoring the Hyperiso framework.
LOG_INFO
#define LOG_INFO(...)
Macro for logging informational messages.
Definition Logger.h:39
BlockProxy
Unified proxy for block queries and diagnostics.
Definition BlockProxy.h:51
BlockProxy::log_all_blocks
void log_all_blocks(ParameterType pt) override
Logs all blocks for a given ParameterType.
Definition BlockProxy.cpp:7
HyperisoMaster
High-level interface to initialize and monitor the main framework configuration.
Definition HyperisoMaster.h:61
HyperisoMaster::init
void init(const std::string &lhaFile, HyperisoConfig config)
Initializes Hyperiso using a LHA file and a full Config object.
Definition HyperisoMaster.cpp:185
WilsonInterface
User-facing API to build and query Wilson coefficients at matching and hadronic scales.
Definition WilsonInterface.h:76
main
int main()
Definition main_mixing.cpp:7
HyperisoConfig
Configuration object controlling model, input flags and optional MARTY resources.
Definition Config.h:24
HyperisoConfig::model
Model model
Current model.
Definition Config.h:33
WilsonBuildConfig
Configuration for building sets of Wilson coefficients.
Definition Configs.h:32