-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathexample.cpp
78 lines (69 loc) · 1.88 KB
/
example.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
// ====================================================================
// This file is part of Himalaya.
//
// Himalaya is licenced under the GNU General Public License (GNU GPL)
// version 3.
// ====================================================================
#include "himalaya/HierarchyCalculator.hpp"
#include <cmath>
#include <iostream>
himalaya::Parameters setup_point(double MS, double tb, double xt)
{
himalaya::Parameters pars;
const double MS2 = MS*MS;
const double Xt = xt*MS;
const double beta = std::atan(tb);
pars.scale = MS;
pars.mu = MS;
pars.g1 = 0.46;
pars.g2 = 0.65;
pars.g3 = 1.166;
pars.vd = 246*std::cos(beta);
pars.vu = 246*std::sin(beta);
pars.mq2 << MS2, 0, 0,
0, MS2, 0,
0, 0, MS2;
pars.md2 << MS2, 0, 0,
0, MS2, 0,
0, 0, MS2;
pars.mu2 << MS2, 0, 0,
0, MS2, 0,
0, 0, MS2;
pars.ml2 << MS2, 0, 0,
0, MS2, 0,
0, 0, MS2;
pars.me2 << MS2, 0, 0,
0, MS2, 0,
0, 0, MS2;
pars.Au << 0, 0, 0,
0, 0, 0,
0, 0, Xt + pars.mu/tb;
pars.Ad << 0, 0, 0,
0, 0, 0,
0, 0, 0;
pars.Ae << 0, 0, 0,
0, 0, 0,
0, 0, 0;
pars.Yu << 0, 0, 0, 0, 0, 0, 0, 0, 0.862;
pars.Yd << 0, 0, 0, 0 ,0 ,0 ,0 ,0, 0.133;
pars.Ye << 0, 0, 0, 0, 0, 0, 0, 0, 0.101;
pars.MA = MS;
pars.M1 = MS;
pars.M2 = MS;
pars.MG = MS;
pars.validate(true);
return pars;
}
int main()
{
try {
const auto point = setup_point(2000., 20., std::sqrt(6.));
himalaya::HierarchyCalculator hc(point);
// calculate the 3-loop corrections O(α_t*α_s^2)
const auto ho = hc.calculateDMh3L(false);
std::cout << ho;
} catch (const std::exception& e) {
std::cerr << e.what() << '\n';
}
return 0;
}