Hyperiso 1.0.3
Modular flavour-physics calculations, Wilson coefficients and statistical inference
Loading...
Searching...
No Matches
speeeed.cpp
Go to the documentation of this file.
1#include "ObservableInterface.h"
2#include "HyperisoMaster.h"
3#include "ParameterSetter.h"
4
5#include <chrono>
6#include <fstream>
7#include <iomanip>
8#include <iostream>
9#include <optional>
10#include <sstream>
11#include <string>
12#include <vector>
13
14struct CsvRow {
15 int step;
16 double mu_b;
17 std::string observable;
18 std::optional<double> q2_min;
19 std::optional<double> q2_max;
20 double value;
21};
22
23int main() {
24 HyperisoMaster hyp;
25 HyperisoConfig config;
26 config.model = Model::SM;
27 hyp.init("lha/si_input.flha", config);
28
29 QCDOrder order = QCDOrder::NNLO;
30 ObservableInterface oi;
31 ParameterSetter ps;
32
33 oi.add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::F_L_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
34 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::A_T_2_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
35 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_2_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
36 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_3_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
37 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_4_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
38 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_5_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
39 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_6_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
40 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_8_B0__KSTAR0_MU_MU), {1.1, 2}}, QCDOrder::NNLO, true)
41
42 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::F_L_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
43 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::A_T_2_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
44 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_2_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
45 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_3_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
46 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_4_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
47 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_5_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
48 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_6_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
49 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_8_B0__KSTAR0_MU_MU), {2, 4.3}}, QCDOrder::NNLO, true)
50
51 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::F_L_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
52 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::A_T_2_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
53 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_2_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
54 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_3_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
55 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_4_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
56 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_5_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
57 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_6_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
58 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_8_B0__KSTAR0_MU_MU), {4.3, 6}}, QCDOrder::NNLO, true)
59
60 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::F_L_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
61 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::A_T_2_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
62 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_2_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
63 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_3_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
64 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_4_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
65 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_5_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
66 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_6_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true)
67 .add_observable(BinnedObservableId{ObservableMapper::to_id(Observables::P_PRIME_8_B0__KSTAR0_MU_MU), {14.18, 16}}, QCDOrder::NNLO, true);
68
69 std::vector<CsvRow> rows;
70 rows.reserve(16 * 32);
71
72 auto start = std::chrono::steady_clock::now();
73 int step = 0;
74
75 for (double mu_b = 1.0; mu_b < 5.0; mu_b += 0.05) {
76 ++step;
77 ps.mutate({ParameterType::WILSON, "B_SCALE", 1}, mu_b);
78 const auto results = oi.compute_all();
79
80 for (const auto& group : results) {
81 const std::string obs_name = group.first.str();
82 for (const auto& val : group.second) {
83 CsvRow row;
84 row.step = step;
85 row.mu_b = mu_b;
86 row.observable = obs_name;
87 if (val.bin.has_value()) {
88 row.q2_min = val.bin->first;
89 row.q2_max = val.bin->second;
90 }
91 row.value = val.value;
92 rows.push_back(std::move(row));
93 }
94 }
95 }
96
97 auto stop = std::chrono::steady_clock::now();
98 const auto us = std::chrono::duration_cast<std::chrono::microseconds>(stop - start).count();
99
100 std::ofstream csv("hyperiso_bkstarll_scan.csv");
101 csv << "step,mu_b,observable,q2_min,q2_max,value\n";
102 csv << std::setprecision(17);
103 for (const auto& row : rows) {
104 csv << row.step << ','
105 << row.mu_b << ','
106 << row.observable << ',';
107 if (row.q2_min.has_value()) {
108 csv << *row.q2_min << ',' << *row.q2_max;
109 } else {
110 csv << ',';
111 }
112 csv << ',' << row.value << '\n';
113 }
114 csv.close();
115
116 std::cout << "Temps calcul seul : " << us * 1e-6 << " s\n";
117 std::cout << "num of steps : " << step << '\n';
118 std::cout << "CSV : hyperiso_bkstarll_scan.csv\n";
119 return 0;
120}
Observables::P_2_B0__KSTAR0_MU_MU
@ P_2_B0__KSTAR0_MU_MU
Observables::P_PRIME_8_B0__KSTAR0_MU_MU
@ P_PRIME_8_B0__KSTAR0_MU_MU
Observables::A_T_2_B0__KSTAR0_MU_MU
@ A_T_2_B0__KSTAR0_MU_MU
Observables::P_PRIME_5_B0__KSTAR0_MU_MU
@ P_PRIME_5_B0__KSTAR0_MU_MU
Observables::F_L_B0__KSTAR0_MU_MU
@ F_L_B0__KSTAR0_MU_MU
Observables::P_3_B0__KSTAR0_MU_MU
@ P_3_B0__KSTAR0_MU_MU
Observables::P_PRIME_4_B0__KSTAR0_MU_MU
@ P_PRIME_4_B0__KSTAR0_MU_MU
Observables::P_PRIME_6_B0__KSTAR0_MU_MU
@ P_PRIME_6_B0__KSTAR0_MU_MU
QCDOrder
QCDOrder
Definition GeneralEnum.h:592
HyperisoMaster.h
High-level helpers for initializing and monitoring the Hyperiso framework.
ObservableInterface.h
High-level, user-facing entry point to compute flavor observables.
GenericMapperWithExt< ObservableTag, Observables, LhaID, std::hash< LhaID >, observable_mapping, observable_flha_mapping >::to_id
static IdOf< ObservableTag > to_id(Observables e)
Converts an enum value to an IdOf<Tag>.
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
ObservableInterface::compute_all
std::map< ObservableId, std::vector< ObservableValue > > compute_all()
Compute all currently registered observables.
Definition ObservableInterface.cpp:284
ObservableInterface::add_observable
ObservableInterface & add_observable(Observables obs, QCDOrder order, bool add_dependencies=false)
Add an observable to the manager (enum API).
Definition ObservableInterface.cpp:91
ParameterSetter
Concrete mutator that directly sets the value of parameters.
Definition ParameterSetter.h:35
ParameterSetter::mutate
void mutate(const ParamId &pid, scalar_t value) override
Sets the value of a parameter.
Definition ParameterSetter.cpp:3
main
int main()
Definition speeeed.cpp:23
BinnedObservableId
Identifies an observable together with a numerical bin.
Definition BinnedObservableId.h:285
CsvRow::mu_b
double mu_b
Definition speeeed.cpp:16
CsvRow::step
int step
Definition speeeed.cpp:15
CsvRow::observable
std::string observable
Definition speeeed.cpp:17
CsvRow::q2_max
std::optional< double > q2_max
Definition speeeed.cpp:19
CsvRow::q2_min
std::optional< double > q2_min
Definition speeeed.cpp:18
CsvRow::value
double value
Definition speeeed.cpp:20
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