-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathDSZHiggs_8cpp_source.html
More file actions
107 lines (105 loc) · 81.8 KB
/
Copy pathDSZHiggs_8cpp_source.html
File metadata and controls
107 lines (105 loc) · 81.8 KB
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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.13"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Himalaya: /home/avoigt/packages/Himalaya/source/himalaya/mh2l/DSZHiggs.cpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td id="projectalign" style="padding-left: 0.5em;">
<div id="projectname">Himalaya
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.13 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
$(function() {
initMenu('',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
</script>
<div id="main-nav"></div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_b2f33c71d4aa5e7af42a1ca61ff5af1b.html">source</a></li><li class="navelem"><a class="el" href="dir_69bdcaeb8e7a1f530411149fc98d532a.html">himalaya</a></li><li class="navelem"><a class="el" href="dir_6082e19c029e3c02659ecf7e77c14555.html">mh2l</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">DSZHiggs.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="DSZHiggs_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">// ====================================================================</span></div><div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment">// This file is part of Himalaya.</span></div><div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">//</span></div><div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment">// Himalaya is licenced under the GNU General Public License (GNU GPL)</span></div><div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment">// version 3.</span></div><div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment">// ====================================================================</span></div><div class="line"><a name="l00007"></a><span class="lineno"> 7</span> </div><div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="preprocessor">#include "<a class="code" href="DSZHiggs_8hpp.html">himalaya/mh2l/DSZHiggs.hpp</a>"</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="preprocessor">#include "<a class="code" href="DSZHiggs_8h.html">himalaya/mh2l/DSZHiggs.h</a>"</span></div><div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="preprocessor">#include "<a class="code" href="Li2f_8hpp.html">himalaya/mh2l/Li2f.hpp</a>"</span></div><div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="preprocessor">#include "<a class="code" href="CouplingOrders_8hpp.html">himalaya/misc/CouplingOrders.hpp</a>"</span></div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="preprocessor">#include <cmath></span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="preprocessor">#include <iostream></span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="preprocessor">#include <limits></span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="preprocessor">#include <mutex></span></div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="preprocessor">#include <utility></span></div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"></span></div><div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment">/**</span></div><div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment"> * @file DSZHiggs.cpp</span></div><div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment"> * @brief Implementation of C++ wrappers for 2-loop FORTRAN routines.</span></div><div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> */</span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> </div><div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="keyword">namespace </span><a class="code" href="namespacehimalaya.html">himalaya</a> {</div><div class="line"><a name="l00024"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html"> 24</a></span> <span class="keyword">namespace </span>mh2l {</div><div class="line"><a name="l00025"></a><span class="lineno"> 25</span> </div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="keyword">namespace </span>{</div><div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment"></span></div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment">/// locks MSSM fortran functions that are not threadsafe</span></div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment"></span>std::mutex mtx;</div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> </div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="keyword">template</span> <<span class="keyword">typename</span> T> T constexpr sqr(T <a class="code" href="Himalaya__LibraryLink_8m.html#a80c72e70807fc6ba175eba2e03353da8">a</a>) { <span class="keywordflow">return</span> a * <a class="code" href="Himalaya__LibraryLink_8m.html#a80c72e70807fc6ba175eba2e03353da8">a</a>; }</div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="keyword">template</span> <<span class="keyword">typename</span> T> T sqrtabs(T a) { <span class="keywordflow">return</span> std::sqrt(std::abs(a)); }</div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="keyword">template</span> <<span class="keyword">typename</span> T> T logabs(T x) { <span class="keywordflow">return</span> <a class="code" href="namespacehimalaya.html#a99d07fbbb2c5cc71320a23984dc65b58">std::log</a>(std::abs(x)); }</div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment"></span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment">/// limit st -> 0 and mst1 -> mst2</span></div><div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment"></span>Eigen::Matrix<double, 2, 2> delta_mh2_2loop_at_as_st_0_mst1_eq_mst2(</div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> <span class="comment">/* mst22 */</span>,</div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span>  <span class="keywordtype">double</span> <span class="comment">/* sxt */</span>, <span class="keywordtype">double</span> <span class="comment">/* cxt */</span>, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3, <span class="keywordtype">int</span> <span class="comment">/* scheme */</span>)</div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span> {</div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span>  constexpr <span class="keywordtype">double</span> Pi2 = M_PI * M_PI;</div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span>  constexpr <span class="keywordtype">double</span> Pi4 = M_PI * M_PI * M_PI * M_PI;</div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2__eft_1_1threshold__loop__functions.html#a20013baed5c82acb480c5697a992756e">g</a> = sqr(mg);</div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> g2 = sqr(g);</div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> q = scale2;</div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> t = mt2;</div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> tsq = sqr(t);</div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> T = mst12;</div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> Tsq = sqr(mst12);</div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> del = g2 + tsq + Tsq - 2*(g*t + g*T + t*T);</div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> rdel = sqrtabs(del);</div><div class="line"><a name="l00052"></a><span class="lineno"> 52</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> sb = std::sin(std::atan(tanb));</div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> ht2 = 2./vev2*mt2/sqr(sb);</div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span> </div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span> </div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span>  result(0,0) = 0.;</div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span> </div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span>  result(0,1) = (sqr(g3)*ht2*mg*mt2*mu* (-1 + logabs(t/q) -</div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span>  (T*((2*del* (-(logabs(t/g)/T) - (2*(g + t - T +</div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span>  g*sqrtabs(del/g2))* logabs((g + t - T - g*sqrtabs(del/g2))/</div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span>  (2.*g)))/ (g*sqrtabs(del/g2)* (t - T + g*(-1 +</div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span>  sqrtabs(del/g2)))) + (2*logabs((g - t + T - g*sqrtabs(del/g2))/</div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span>  (2.*g)))/(g*sqrtabs(del/g2)) - (2*((g + t - T +</div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span>  g*sqrtabs(del/g2))* logabs((g + t - T + g*sqrtabs(del/g2))/</div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span>  (2.*g)) + (g - t + T - g*sqrtabs(del/g2))* logabs((g - t + T +</div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span>  g*sqrtabs(del/g2))/ (2.*g))))/ (g*sqrtabs(del/g2)* (t - T +</div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span>  g*(-1 + sqrtabs(del/g2))))))/ g2 + (2*(g + t - T)*(Pi2 -</div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span>  3*logabs(t/g)*logabs(T/g) + 6*logabs((g + t - T -</div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span>  g*sqrtabs(del/g2))/(2.*g))*logabs((g - t + T -</div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span>  g*sqrtabs(del/g2))/(2.*g)) - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g + t - T -</div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span>  g*sqrtabs(del/g2))/(2.*g)) - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - t + T -</div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span>  g*sqrtabs(del/g2))/(2.*g))))/ (3.*g2)))/(2.*std::pow(std::abs(del)/g2,1.5))))/</div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span>  (8.*Pi4*T);</div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span> </div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span> </div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span>  result(1,1) = (sqr(g3)*ht2*mt2*(-2 - (8*(g + t)*(-1 +</div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span>  logabs(g/q)))/T + 8*logabs(t/g) + 6*(-1 + logabs(t/q)) +</div><div class="line"><a name="l00080"></a><span class="lineno"> 80</span>  8*sqr(logabs(t/q)) - 4*logabs(T/g) + (4*((g2*T - sqr(t - T)*(2*t</div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span>  + T) + 2*g*t*(t + 5*T))*logabs(t/g) +</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  4*g2*T*logabs(T/g)))/(del*T) + 2*logabs(T/q) -</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  8*sqr(logabs(T/q)) + 5*logabs(Tsq/tsq) + sqr(logabs(Tsq/tsq)) +</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  (4*g2*(g + t - T)*(Pi2 - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g + t - T -</div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span>  g*sqrtabs(del/g2))/ (2.*g)) - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - t + T -</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span>  g*sqrtabs(del/g2))/(2.*g)) - 3*logabs(t/g)*logabs(T/g) +</div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span>  6*logabs((-rdel + g + t - T)/(2.*g))* logabs((-rdel + g - t +</div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span>  T)/(2.*g))))/(3.*std::pow(std::abs(del),1.5)) + (4*g*(g + t - T)*(Pi2 -</div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span>  6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g + t - T - g*sqrtabs(del/g2))/(2.*g)) - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - t</div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span>  + T - g*sqrtabs(del/g2))/ (2.*g)) - 3*logabs(t/g)*logabs(T/g) +</div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span>  6*logabs((g + t - T - g*sqrtabs(del/g2))/(2.*g))* logabs((g - t</div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span>  + T - g*sqrtabs(del/g2))/(2.*g))))/ (3.*del*sqrtabs(del/g2)) +</div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span>  (8*mg*mu*(1/T - logabs(t/q)/T + (g*(g + t - T)* (Pi2 -</div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g + t - T - g*sqrtabs(del/g2))/(2.*g)) - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - t</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span>  + T - g*sqrtabs(del/g2))/ (2.*g)) - 3*logabs(t/g)*logabs(T/g) +</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  6*logabs((-rdel + g + t - T)/(2.*g))*logabs((-rdel + g - t +</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  T)/(2.*g))))/ (3.*std::pow(std::abs(del),1.5)) + (g*(-(logabs(t/g)/T) -</div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  (2*(-(g*logabs(4.)) + (rdel + g + t - T)*logabs((-rdel + g + t -</div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  T)/g) + (-rdel + g - t + T)*logabs((-rdel + g - t + T)/g)))/</div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  (rdel*(rdel - g + t - T)) - 2*(((rdel + g + t - T)* logabs((g +</div><div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  t - T + g*sqrtabs(del/g2))/ (2.*g)))/ (rdel*(t - T + g*(-1 +</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  sqrtabs(del/g2)))) - ((rdel - g - t + T)* logabs((g - t + T +</div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  g*sqrtabs(del/g2))/ (2.*g)))/ (rdel*(-t + T + g*(-1 +</div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  sqrtabs(del/g2)))))))/ rdel))/tanb))/(32.*Pi4);</div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span> </div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span> }</div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span> <span class="comment"></span></div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span> <span class="comment">/// Pietro Slavich implementation</span></div><div class="line"><a name="l00110"></a><span class="lineno"> 110</span> <span class="comment"></span>Eigen::Matrix<double, 2, 2> delta_mh2_2loop_at_as_general(</div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3, <span class="keywordtype">int</span> scheme)</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span> {</div><div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00116"></a><span class="lineno"> 116</span> </div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <a class="code" href="DSZHiggs_8h.html#a88e8c8f0df22394855df1dfc61564c2d">dszhiggs_</a>(&mt2, &mg, &mst12, &mst22, &sxt, &cxt, &scale2, &mu,</div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  &tanb, &vev2, &g3, &scheme,</div><div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  &result(0,0), &result(1,1), &result(0,1));</div><div class="line"><a name="l00120"></a><span class="lineno"> 120</span> </div><div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span> </div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span> }</div><div class="line"><a name="l00125"></a><span class="lineno"> 125</span> </div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span> Eigen::Matrix<double, 2, 2> rotate_by(<span class="keywordtype">double</span> dm2, <span class="keywordtype">double</span> tanb)</div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span> {</div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> tanb2 = sqr(tanb);</div><div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> sinb = tanb / sqrtabs(1. + tanb2);</div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> cosb = 1. / sqrtabs(1. + tanb2);</div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span> </div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00133"></a><span class="lineno"> 133</span> </div><div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  result(0,0) = dm2 * sqr(sinb);</div><div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  result(0,1) = - dm2 * sinb * cosb;</div><div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  result(1,1) = dm2 * sqr(cosb);</div><div class="line"><a name="l00138"></a><span class="lineno"> 138</span> </div><div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00140"></a><span class="lineno"> 140</span> }</div><div class="line"><a name="l00141"></a><span class="lineno"> 141</span> </div><div class="line"><a name="l00142"></a><span class="lineno"> 142</span> <span class="keywordtype">double</span> calc_At(<span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb)</div><div class="line"><a name="l00144"></a><span class="lineno"> 144</span> {</div><div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> s2t = 2*cxt*sxt;</div><div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> Xt = (mst12 - mst22)*s2t/2./sqrtabs(mt2);</div><div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth.html#ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df">At</a> = Xt - mu/tanb;</div><div class="line"><a name="l00148"></a><span class="lineno"> 148</span> </div><div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  <span class="keywordflow">return</span> <a class="code" href="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth.html#ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df">At</a>;</div><div class="line"><a name="l00150"></a><span class="lineno"> 150</span> }</div><div class="line"><a name="l00151"></a><span class="lineno"> 151</span> </div><div class="line"><a name="l00152"></a><span class="lineno"> 152</span> <span class="keywordtype">double</span> <a class="code" href="DSZHiggs_8f.html#a6120448a5a0919fd26c80580dce5c0f3">phi</a>(<span class="keywordtype">double</span> x, <span class="keywordtype">double</span> y, <span class="keywordtype">double</span> z)</div><div class="line"><a name="l00153"></a><span class="lineno"> 153</span> {</div><div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> u = x/z, v = y/z;</div><div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> lambda = sqrtabs(sqr(1 - u - v) - 4*u*v);</div><div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> xp = 0.5 * (1 + (u - v) - lambda);</div><div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> xm = 0.5 * (1 - (u - v) - lambda);</div><div class="line"><a name="l00158"></a><span class="lineno"> 158</span> </div><div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  <span class="keywordflow">return</span> 1./lambda * (2*logabs(xp)*logabs(xm) - logabs(u)*logabs(v) -</div><div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  2*(<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>(xp) + <a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>(xm)) + M_PI*M_PI/3.);</div><div class="line"><a name="l00161"></a><span class="lineno"> 161</span> }</div><div class="line"><a name="l00162"></a><span class="lineno"> 162</span> <span class="comment"></span></div><div class="line"><a name="l00163"></a><span class="lineno"> 163</span> <span class="comment">/// First derivative of phi[t,T,g] w.r.t. T</span></div><div class="line"><a name="l00164"></a><span class="lineno"> 164</span> <span class="comment"></span><span class="keywordtype">double</span> dphi_010(<span class="keywordtype">double</span> t, <span class="keywordtype">double</span> T, <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2__eft_1_1threshold__loop__functions.html#a20013baed5c82acb480c5697a992756e">g</a>)</div><div class="line"><a name="l00165"></a><span class="lineno"> 165</span> {</div><div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  constexpr <span class="keywordtype">double</span> Pi2 = M_PI * M_PI;</div><div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> g2 = sqr(g);</div><div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> abbr = (-4*t*T)/g2 + sqr(1 - t/g - T/g);</div><div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> rabbr = sqrtabs(abbr);</div><div class="line"><a name="l00170"></a><span class="lineno"> 170</span> </div><div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  <span class="keywordflow">return</span> ((g + t - T)*(Pi2 - 6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - rabbr*g + t - T)/(2.*g)) -</div><div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  6*<a class="code" href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">li2</a>((g - rabbr*g - t + T)/(2.*g)) -</div><div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  3*logabs(t/g)*logabs(T/g) + 6*logabs((g - rabbr*g + t -</div><div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  T)/(2.*g))*logabs((g - rabbr*g - t + T)/(2.*g))) + (3*rabbr*g* (</div><div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  rabbr*g*((-1 + rabbr)*g + t - T)*logabs(t/g) +</div><div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  2*T*(-2*g*logabs(4.) + (g + rabbr*g + t - T)*logabs((g - rabbr*g</div><div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  + t - T)/g) + (g + rabbr*g + t - T)*logabs((g + rabbr*g + t -</div><div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  T)/g) + g*logabs((g - rabbr*g - t + T)/g) - rabbr*g*logabs((g -</div><div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  rabbr*g - t + T)/g) - t*logabs((g - rabbr*g - t + T)/g) +</div><div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  T*logabs((g - rabbr*g - t + T)/g) + g*logabs((g + rabbr*g - t +</div><div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  T)/g) - rabbr*g*logabs((g + rabbr*g - t + T)/g) - t*logabs((g +</div><div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  rabbr*g - t + T)/g) + T*logabs((g + rabbr*g - t + T)/g)) ) ) /</div><div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  (T*(g - rabbr*g - t + T)))/(3.*std::pow(std::abs(abbr),1.5)*g2);</div><div class="line"><a name="l00184"></a><span class="lineno"> 184</span> }</div><div class="line"><a name="l00185"></a><span class="lineno"> 185</span> </div><div class="line"><a name="l00186"></a><span class="lineno"> 186</span> } <span class="comment">// anonymous namespace</span></div><div class="line"><a name="l00187"></a><span class="lineno"> 187</span> </div><div class="line"><a name="l00188"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#ad0ce42f717a9918763604fd55fcb4e3d"> 188</a></span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#ad0ce42f717a9918763604fd55fcb4e3d">delta_ma2_2loop_at_as_mst1_eq_mst2</a>(</div><div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3)</div><div class="line"><a name="l00192"></a><span class="lineno"> 192</span> {</div><div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  constexpr <span class="keywordtype">double</span> Pi2 = M_PI * M_PI;</div><div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2__eft_1_1threshold__loop__functions.html#a20013baed5c82acb480c5697a992756e">g</a> = sqr(mg);</div><div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> g2 = sqr(g);</div><div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> q = scale2;</div><div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> q2 = sqr(scale2);</div><div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> t = mt2;</div><div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> T = mst12;</div><div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> sb = std::sin(std::atan(tanb));</div><div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> ht2 = 2./vev2*mt2/sqr(sb);</div><div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth.html#ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df">At</a> = calc_At(mt2, mst12, mst22, sxt, cxt, mu, tanb);</div><div class="line"><a name="l00203"></a><span class="lineno"> 203</span> </div><div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> result = (-2*(g*(2*At*g + 2*At*t - At*T + mg*T + mg*(g</div><div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  - t)*logabs(g/t) - At*T*logabs(g/q)*logabs(t/q) -</div><div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  mg*T*logabs(g/q)*logabs(t/q) - 4*mg*T*logabs(T/q) -</div><div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  2*At*T*sqr(logabs(T/q)) + logabs((g*t)/q2)*(-(At*(g + t - T)) +</div><div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  mg*T + (At + mg)*T*logabs(T/q))) - 2*(At + mg)*(g + t -</div><div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  T)*T*<a class="code" href="DSZHiggs_8f.html#a6120448a5a0919fd26c80580dce5c0f3">phi</a>(t,T,g) + T*(At*(g2 + sqr(t - T) - 2*g*T) + mg*(g2 +</div><div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  sqr(t - T) - 2*g*(t + T)))*dphi_010(t,T,g)))/ (g*T);</div><div class="line"><a name="l00211"></a><span class="lineno"> 211</span> </div><div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> pref = 4*sqr(g3)/sqr(16*Pi2) * ht2*mu*(1./tanb + tanb);</div><div class="line"><a name="l00213"></a><span class="lineno"> 213</span> </div><div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  <span class="keywordflow">return</span> pref * result;</div><div class="line"><a name="l00215"></a><span class="lineno"> 215</span> }</div><div class="line"><a name="l00216"></a><span class="lineno"> 216</span> </div><div class="line"><a name="l00217"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#accc846c9b031c5d40de467d16e2f5054"> 217</a></span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#accc846c9b031c5d40de467d16e2f5054">delta_ma2_2loop_at_as_general</a>(</div><div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3)</div><div class="line"><a name="l00221"></a><span class="lineno"> 221</span> {</div><div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  <span class="keywordtype">double</span> result = 0.;</div><div class="line"><a name="l00223"></a><span class="lineno"> 223</span> </div><div class="line"><a name="l00224"></a><span class="lineno"> 224</span>  <a class="code" href="DSZHiggs_8h.html#ae6e9056d395834b500d50e3219a78fdb">dszodd_</a>(&mt2, &mg, &mst12, &mst22, &sxt, &cxt, &scale2, &mu,</div><div class="line"><a name="l00225"></a><span class="lineno"> 225</span>  &tanb, &vev2, &g3, &result);</div><div class="line"><a name="l00226"></a><span class="lineno"> 226</span> </div><div class="line"><a name="l00227"></a><span class="lineno"> 227</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00228"></a><span class="lineno"> 228</span> }</div><div class="line"><a name="l00229"></a><span class="lineno"> 229</span> <span class="comment"></span></div><div class="line"><a name="l00230"></a><span class="lineno"> 230</span> <span class="comment">/// 2-loop contribution to CP-odd Higgs O(at*as)</span></div><div class="line"><a name="l00231"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a9561889cc35e34a709643cfed97856cc"> 231</a></span> <span class="comment"></span><span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#a9561889cc35e34a709643cfed97856cc">delta_ma2_2loop_at_as</a>(</div><div class="line"><a name="l00232"></a><span class="lineno"> 232</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00233"></a><span class="lineno"> 233</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00234"></a><span class="lineno"> 234</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3)</div><div class="line"><a name="l00235"></a><span class="lineno"> 235</span> {</div><div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  <span class="keywordflow">if</span> (std::abs((mst12 - mst22)/mst12) < 1e-8) {</div><div class="line"><a name="l00237"></a><span class="lineno"> 237</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth.html#ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df">At</a> = calc_At(mt2, mst12, mst22, sxt, cxt, mu, tanb);</div><div class="line"><a name="l00238"></a><span class="lineno"> 238</span> </div><div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <span class="comment">// if At = 0 => mu = 0 => dMA(2L) = 0</span></div><div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  <span class="keywordflow">if</span> (std::abs(At) < std::numeric_limits<double>::epsilon())</div><div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  <span class="keywordflow">return</span> 0.;</div><div class="line"><a name="l00242"></a><span class="lineno"> 242</span> </div><div class="line"><a name="l00243"></a><span class="lineno"> 243</span>  <span class="keywordflow">return</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#ad0ce42f717a9918763604fd55fcb4e3d">delta_ma2_2loop_at_as_mst1_eq_mst2</a>(</div><div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  mt2, mg, mst12, mst22, sxt, cxt, scale2, mu, tanb, vev2, g3);</div><div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  }</div><div class="line"><a name="l00246"></a><span class="lineno"> 246</span> </div><div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  <span class="keywordflow">return</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#accc846c9b031c5d40de467d16e2f5054">delta_ma2_2loop_at_as_general</a>(</div><div class="line"><a name="l00248"></a><span class="lineno"> 248</span>  mt2, mg, mst12, mst22, sxt, cxt, scale2, mu, tanb, vev2, g3);</div><div class="line"><a name="l00249"></a><span class="lineno"> 249</span> }</div><div class="line"><a name="l00250"></a><span class="lineno"> 250</span> </div><div class="line"><a name="l00251"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#acbdc52f7ab258a55f03b9f8bd069784e"> 251</a></span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#acbdc52f7ab258a55f03b9f8bd069784e">delta_ma2_2loop_at_at</a>(</div><div class="line"><a name="l00252"></a><span class="lineno"> 252</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mb2, <span class="keywordtype">double</span> mA2, <span class="keywordtype">double</span> mst12,</div><div class="line"><a name="l00253"></a><span class="lineno"> 253</span>  <span class="keywordtype">double</span> mst22, <span class="keywordtype">double</span> msb12, <span class="keywordtype">double</span> msb22,</div><div class="line"><a name="l00254"></a><span class="lineno"> 254</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> sxb, <span class="keywordtype">double</span> cxb,</div><div class="line"><a name="l00255"></a><span class="lineno"> 255</span>  <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">int</span> atat)</div><div class="line"><a name="l00256"></a><span class="lineno"> 256</span> {</div><div class="line"><a name="l00257"></a><span class="lineno"> 257</span>  <span class="keywordtype">double</span> result = 0.;</div><div class="line"><a name="l00258"></a><span class="lineno"> 258</span> </div><div class="line"><a name="l00259"></a><span class="lineno"> 259</span>  {</div><div class="line"><a name="l00260"></a><span class="lineno"> 260</span>  std::lock_guard<std::mutex> lg(mtx);</div><div class="line"><a name="l00261"></a><span class="lineno"> 261</span> </div><div class="line"><a name="l00262"></a><span class="lineno"> 262</span>  <a class="code" href="DSZHiggs_8h.html#accd0e91b067d839181665628d66d932f">ddsodd_</a>(&mt2, &mb2, &mA2, &mst12, &mst22, &msb12, &msb22,</div><div class="line"><a name="l00263"></a><span class="lineno"> 263</span>  &sxt, &cxt, &sxb, &cxb, &scale2, &mu, &tanb, &vev2, &result, &atat);</div><div class="line"><a name="l00264"></a><span class="lineno"> 264</span>  }</div><div class="line"><a name="l00265"></a><span class="lineno"> 265</span> </div><div class="line"><a name="l00266"></a><span class="lineno"> 266</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00267"></a><span class="lineno"> 267</span> }</div><div class="line"><a name="l00268"></a><span class="lineno"> 268</span> </div><div class="line"><a name="l00269"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#acc1ac4c68ee565a17377005fc22fa312"> 269</a></span> Eigen::Matrix<double, 2, 2> <a class="code" href="namespacehimalaya_1_1mh2l.html#acc1ac4c68ee565a17377005fc22fa312">delta_mh2_2loop_at_as</a>(</div><div class="line"><a name="l00270"></a><span class="lineno"> 270</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> mst12, <span class="keywordtype">double</span> mst22,</div><div class="line"><a name="l00271"></a><span class="lineno"> 271</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00272"></a><span class="lineno"> 272</span>  <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3,</div><div class="line"><a name="l00273"></a><span class="lineno"> 273</span>  <span class="keywordtype">int</span> include_heavy_higgs)</div><div class="line"><a name="l00274"></a><span class="lineno"> 274</span> {</div><div class="line"><a name="l00275"></a><span class="lineno"> 275</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00276"></a><span class="lineno"> 276</span> </div><div class="line"><a name="l00277"></a><span class="lineno"> 277</span>  <span class="keywordflow">if</span> (std::abs((mst12 - mst22)/mst12) < 1e-8) {</div><div class="line"><a name="l00278"></a><span class="lineno"> 278</span>  result = delta_mh2_2loop_at_as_st_0_mst1_eq_mst2(</div><div class="line"><a name="l00279"></a><span class="lineno"> 279</span>  mt2, mg, mst12, mst22, sxt, cxt, scale2, mu, tanb, vev2, g3, 0);</div><div class="line"><a name="l00280"></a><span class="lineno"> 280</span>  } <span class="keywordflow">else</span> {</div><div class="line"><a name="l00281"></a><span class="lineno"> 281</span>  result = delta_mh2_2loop_at_as_general(</div><div class="line"><a name="l00282"></a><span class="lineno"> 282</span>  mt2, mg, mst12, mst22, sxt, cxt, scale2, mu, tanb, vev2, g3, 0);</div><div class="line"><a name="l00283"></a><span class="lineno"> 283</span>  }</div><div class="line"><a name="l00284"></a><span class="lineno"> 284</span> </div><div class="line"><a name="l00285"></a><span class="lineno"> 285</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> dMA = include_heavy_higgs</div><div class="line"><a name="l00286"></a><span class="lineno"> 286</span>  ? <a class="code" href="namespacehimalaya_1_1mh2l.html#a9561889cc35e34a709643cfed97856cc">delta_ma2_2loop_at_as</a>(</div><div class="line"><a name="l00287"></a><span class="lineno"> 287</span>  mt2, mg, mst12, mst22, sxt, cxt, scale2, mu, tanb, vev2, g3)</div><div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  : 0.;</div><div class="line"><a name="l00289"></a><span class="lineno"> 289</span> </div><div class="line"><a name="l00290"></a><span class="lineno"> 290</span>  <span class="keywordflow">return</span> result + rotate_by(dMA, tanb);</div><div class="line"><a name="l00291"></a><span class="lineno"> 291</span> }</div><div class="line"><a name="l00292"></a><span class="lineno"> 292</span> </div><div class="line"><a name="l00293"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#accc07f22b2d75d17656ea35acd713c2b"> 293</a></span> Eigen::Matrix<double, 2, 2> <a class="code" href="namespacehimalaya_1_1mh2l.html#accc07f22b2d75d17656ea35acd713c2b">delta_mh2_2loop_at_at</a>(</div><div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  <span class="keywordtype">double</span> mt2, <span class="keywordtype">double</span> mb2, <span class="keywordtype">double</span> mA2, <span class="keywordtype">double</span> mst12,</div><div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  <span class="keywordtype">double</span> mst22, <span class="keywordtype">double</span> msb12, <span class="keywordtype">double</span> msb22,</div><div class="line"><a name="l00296"></a><span class="lineno"> 296</span>  <span class="keywordtype">double</span> sxt, <span class="keywordtype">double</span> cxt, <span class="keywordtype">double</span> sxb, <span class="keywordtype">double</span> cxb,</div><div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2,</div><div class="line"><a name="l00298"></a><span class="lineno"> 298</span>  <span class="keywordtype">int</span> include_heavy_higgs, <span class="keywordtype">int</span> atat)</div><div class="line"><a name="l00299"></a><span class="lineno"> 299</span> {</div><div class="line"><a name="l00300"></a><span class="lineno"> 300</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00301"></a><span class="lineno"> 301</span> </div><div class="line"><a name="l00302"></a><span class="lineno"> 302</span>  {</div><div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  std::lock_guard<std::mutex> lg(mtx);</div><div class="line"><a name="l00304"></a><span class="lineno"> 304</span> </div><div class="line"><a name="l00305"></a><span class="lineno"> 305</span>  <a class="code" href="DSZHiggs_8h.html#a62881cea97186f48509abe184f224dcb">ddshiggs_</a>(&mt2, &mb2, &mA2, &mst12, &mst22, &msb12, &msb22,</div><div class="line"><a name="l00306"></a><span class="lineno"> 306</span>  &sxt, &cxt, &sxb, &cxb, &scale2, &mu, &tanb, &vev2,</div><div class="line"><a name="l00307"></a><span class="lineno"> 307</span>  &result(0,0), &result(0,1), &result(1,1), &atat);</div><div class="line"><a name="l00308"></a><span class="lineno"> 308</span>  }</div><div class="line"><a name="l00309"></a><span class="lineno"> 309</span> </div><div class="line"><a name="l00310"></a><span class="lineno"> 310</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00311"></a><span class="lineno"> 311</span> </div><div class="line"><a name="l00312"></a><span class="lineno"> 312</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> dMA = include_heavy_higgs</div><div class="line"><a name="l00313"></a><span class="lineno"> 313</span>  ? <a class="code" href="namespacehimalaya_1_1mh2l.html#acbdc52f7ab258a55f03b9f8bd069784e">delta_ma2_2loop_at_at</a>(</div><div class="line"><a name="l00314"></a><span class="lineno"> 314</span>  mt2, mb2, mA2, mst12, mst22, msb12, msb22,</div><div class="line"><a name="l00315"></a><span class="lineno"> 315</span>  sxt, cxt, sxb, cxb, scale2, mu, tanb, vev2, atat)</div><div class="line"><a name="l00316"></a><span class="lineno"> 316</span>  : 0.;</div><div class="line"><a name="l00317"></a><span class="lineno"> 317</span> </div><div class="line"><a name="l00318"></a><span class="lineno"> 318</span>  <span class="keywordflow">return</span> result + rotate_by(dMA, tanb);</div><div class="line"><a name="l00319"></a><span class="lineno"> 319</span> }</div><div class="line"><a name="l00320"></a><span class="lineno"> 320</span> </div><div class="line"><a name="l00321"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a83dba77a0cf5c8e1e1f0dc89c94360b3"> 321</a></span> Eigen::Matrix<double, 2, 2> <a class="code" href="namespacehimalaya_1_1mh2l.html#a83dba77a0cf5c8e1e1f0dc89c94360b3">delta_mh2_2loop_ab_as</a>(</div><div class="line"><a name="l00322"></a><span class="lineno"> 322</span>  <span class="keywordtype">double</span> mb2, <span class="keywordtype">double</span> mg, <span class="keywordtype">double</span> msb12, <span class="keywordtype">double</span> msb22,</div><div class="line"><a name="l00323"></a><span class="lineno"> 323</span>  <span class="keywordtype">double</span> sxb, <span class="keywordtype">double</span> cxb, <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu,</div><div class="line"><a name="l00324"></a><span class="lineno"> 324</span>  <span class="keywordtype">double</span> cotb, <span class="keywordtype">double</span> vev2, <span class="keywordtype">double</span> g3,</div><div class="line"><a name="l00325"></a><span class="lineno"> 325</span>  <span class="keywordtype">int</span> include_heavy_higgs)</div><div class="line"><a name="l00326"></a><span class="lineno"> 326</span> {</div><div class="line"><a name="l00327"></a><span class="lineno"> 327</span>  Eigen::Matrix<double, 2, 2> result(<a class="code" href="namespacehimalaya_1_1mh2l.html#acc1ac4c68ee565a17377005fc22fa312">delta_mh2_2loop_at_as</a>(</div><div class="line"><a name="l00328"></a><span class="lineno"> 328</span>  mb2, mg, msb12, msb22, sxb, cxb, scale2, mu,</div><div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  cotb, vev2, g3, include_heavy_higgs));</div><div class="line"><a name="l00330"></a><span class="lineno"> 330</span> </div><div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  std::swap(result(0,0), result(1,1));</div><div class="line"><a name="l00332"></a><span class="lineno"> 332</span> </div><div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00334"></a><span class="lineno"> 334</span> }</div><div class="line"><a name="l00335"></a><span class="lineno"> 335</span> </div><div class="line"><a name="l00336"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a5ddaca4b030e7354fd06bc149bbaf13a"> 336</a></span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#a5ddaca4b030e7354fd06bc149bbaf13a">delta_ma2_2loop_atau_atau</a>(</div><div class="line"><a name="l00337"></a><span class="lineno"> 337</span>  <span class="keywordtype">double</span> mtau2, <span class="keywordtype">double</span> mA2, <span class="keywordtype">double</span> msv2, <span class="keywordtype">double</span> mstau12,</div><div class="line"><a name="l00338"></a><span class="lineno"> 338</span>  <span class="keywordtype">double</span> mstau22, <span class="keywordtype">double</span> sintau, <span class="keywordtype">double</span> costau, <span class="keywordtype">double</span> scale2,</div><div class="line"><a name="l00339"></a><span class="lineno"> 339</span>  <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2)</div><div class="line"><a name="l00340"></a><span class="lineno"> 340</span> {</div><div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  <span class="keywordtype">double</span> result = 0.;</div><div class="line"><a name="l00342"></a><span class="lineno"> 342</span> </div><div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <a class="code" href="DSZHiggs_8h.html#a20219899e8065151504afdc38c83e481">tausqodd_</a>(&mtau2, &mA2, &msv2, &mstau12, &mstau22, &sintau,</div><div class="line"><a name="l00344"></a><span class="lineno"> 344</span>  &costau, &scale2, &mu, &tanb, &vev2, &result);</div><div class="line"><a name="l00345"></a><span class="lineno"> 345</span> </div><div class="line"><a name="l00346"></a><span class="lineno"> 346</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00347"></a><span class="lineno"> 347</span> }</div><div class="line"><a name="l00348"></a><span class="lineno"> 348</span> </div><div class="line"><a name="l00349"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a21677404a27a0054d6a0117db96a181d"> 349</a></span> Eigen::Matrix<double, 2, 2> <a class="code" href="namespacehimalaya_1_1mh2l.html#a21677404a27a0054d6a0117db96a181d">delta_mh2_2loop_atau_atau</a>(</div><div class="line"><a name="l00350"></a><span class="lineno"> 350</span>  <span class="keywordtype">double</span> mtau2, <span class="keywordtype">double</span> mA2, <span class="keywordtype">double</span> msv2, <span class="keywordtype">double</span> mstau12,</div><div class="line"><a name="l00351"></a><span class="lineno"> 351</span>  <span class="keywordtype">double</span> mstau22, <span class="keywordtype">double</span> sintau, <span class="keywordtype">double</span> costau, <span class="keywordtype">double</span> scale2,</div><div class="line"><a name="l00352"></a><span class="lineno"> 352</span>  <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2,</div><div class="line"><a name="l00353"></a><span class="lineno"> 353</span>  <span class="keywordtype">int</span> include_heavy_higgs)</div><div class="line"><a name="l00354"></a><span class="lineno"> 354</span> {</div><div class="line"><a name="l00355"></a><span class="lineno"> 355</span>  <span class="keywordtype">int</span> scheme = 0;</div><div class="line"><a name="l00356"></a><span class="lineno"> 356</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00357"></a><span class="lineno"> 357</span> </div><div class="line"><a name="l00358"></a><span class="lineno"> 358</span>  <a class="code" href="DSZHiggs_8h.html#afebc9638df87be653b7998cd246cc4a3">tausqhiggs_</a>(&mtau2, &mA2, &msv2, &mstau12, &mstau22, &sintau,</div><div class="line"><a name="l00359"></a><span class="lineno"> 359</span>  &costau, &scale2, &mu, &tanb, &vev2, &scheme,</div><div class="line"><a name="l00360"></a><span class="lineno"> 360</span>  &result(0,0), &result(1,1), &result(0,1));</div><div class="line"><a name="l00361"></a><span class="lineno"> 361</span> </div><div class="line"><a name="l00362"></a><span class="lineno"> 362</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00363"></a><span class="lineno"> 363</span> </div><div class="line"><a name="l00364"></a><span class="lineno"> 364</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> dMA = include_heavy_higgs</div><div class="line"><a name="l00365"></a><span class="lineno"> 365</span>  ? <a class="code" href="namespacehimalaya_1_1mh2l.html#a5ddaca4b030e7354fd06bc149bbaf13a">delta_ma2_2loop_atau_atau</a>(</div><div class="line"><a name="l00366"></a><span class="lineno"> 366</span>  mtau2, mA2, msv2, mstau12, mstau22, sintau, costau, scale2,</div><div class="line"><a name="l00367"></a><span class="lineno"> 367</span>  mu, tanb, vev2)</div><div class="line"><a name="l00368"></a><span class="lineno"> 368</span>  : 0.;</div><div class="line"><a name="l00369"></a><span class="lineno"> 369</span> </div><div class="line"><a name="l00370"></a><span class="lineno"> 370</span>  <span class="keywordflow">return</span> result + rotate_by(dMA, tanb);</div><div class="line"><a name="l00371"></a><span class="lineno"> 371</span> }</div><div class="line"><a name="l00372"></a><span class="lineno"> 372</span> </div><div class="line"><a name="l00373"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a8b600e41354ff1538cbe86668852e468"> 373</a></span> <span class="keywordtype">double</span> <a class="code" href="namespacehimalaya_1_1mh2l.html#a8b600e41354ff1538cbe86668852e468">delta_ma2_2loop_ab_atau</a>(</div><div class="line"><a name="l00374"></a><span class="lineno"> 374</span>  <span class="keywordtype">double</span> mtau2, <span class="keywordtype">double</span> mb2,</div><div class="line"><a name="l00375"></a><span class="lineno"> 375</span>  <span class="keywordtype">double</span> mstau12, <span class="keywordtype">double</span> mstau22, <span class="keywordtype">double</span> msb12, <span class="keywordtype">double</span> msb22,</div><div class="line"><a name="l00376"></a><span class="lineno"> 376</span>  <span class="keywordtype">double</span> sintau, <span class="keywordtype">double</span> costau, <span class="keywordtype">double</span> sxb, <span class="keywordtype">double</span> cxb,</div><div class="line"><a name="l00377"></a><span class="lineno"> 377</span>  <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2)</div><div class="line"><a name="l00378"></a><span class="lineno"> 378</span> {</div><div class="line"><a name="l00379"></a><span class="lineno"> 379</span>  <span class="keywordtype">double</span> result = 0.;</div><div class="line"><a name="l00380"></a><span class="lineno"> 380</span> </div><div class="line"><a name="l00381"></a><span class="lineno"> 381</span>  <a class="code" href="DSZHiggs_8h.html#a9f210e554b67f63390a67ff1b60540e2">taubotodd_</a>(&mtau2, &mb2,</div><div class="line"><a name="l00382"></a><span class="lineno"> 382</span>  &mstau12, &mstau22, &msb12, &msb22,</div><div class="line"><a name="l00383"></a><span class="lineno"> 383</span>  &sintau, &costau, &sxb, &cxb,</div><div class="line"><a name="l00384"></a><span class="lineno"> 384</span>  &scale2, &mu, &tanb, &vev2, &result);</div><div class="line"><a name="l00385"></a><span class="lineno"> 385</span> </div><div class="line"><a name="l00386"></a><span class="lineno"> 386</span>  <span class="keywordflow">return</span> result;</div><div class="line"><a name="l00387"></a><span class="lineno"> 387</span> }</div><div class="line"><a name="l00388"></a><span class="lineno"> 388</span> </div><div class="line"><a name="l00389"></a><span class="lineno"><a class="line" href="namespacehimalaya_1_1mh2l.html#a7827d302c444d22f1b46e1c855566039"> 389</a></span> Eigen::Matrix<double, 2, 2> <a class="code" href="namespacehimalaya_1_1mh2l.html#a7827d302c444d22f1b46e1c855566039">delta_mh2_2loop_ab_atau</a>(</div><div class="line"><a name="l00390"></a><span class="lineno"> 390</span>  <span class="keywordtype">double</span> mtau2, <span class="keywordtype">double</span> mb2,</div><div class="line"><a name="l00391"></a><span class="lineno"> 391</span>  <span class="keywordtype">double</span> mstau12, <span class="keywordtype">double</span> mstau22, <span class="keywordtype">double</span> msb12, <span class="keywordtype">double</span> msb22,</div><div class="line"><a name="l00392"></a><span class="lineno"> 392</span>  <span class="keywordtype">double</span> sintau, <span class="keywordtype">double</span> costau, <span class="keywordtype">double</span> sxb, <span class="keywordtype">double</span> cxb,</div><div class="line"><a name="l00393"></a><span class="lineno"> 393</span>  <span class="keywordtype">double</span> scale2, <span class="keywordtype">double</span> mu, <span class="keywordtype">double</span> tanb, <span class="keywordtype">double</span> vev2,</div><div class="line"><a name="l00394"></a><span class="lineno"> 394</span>  <span class="keywordtype">int</span> include_heavy_higgs)</div><div class="line"><a name="l00395"></a><span class="lineno"> 395</span> {</div><div class="line"><a name="l00396"></a><span class="lineno"> 396</span>  Eigen::Matrix<double, 2, 2> result;</div><div class="line"><a name="l00397"></a><span class="lineno"> 397</span> </div><div class="line"><a name="l00398"></a><span class="lineno"> 398</span>  <a class="code" href="DSZHiggs_8h.html#acbca846367a332a10c3d4c91a8439ab8">taubot_</a>(&mtau2, &mb2,</div><div class="line"><a name="l00399"></a><span class="lineno"> 399</span>  &mstau12, &mstau22, &msb12, &msb22,</div><div class="line"><a name="l00400"></a><span class="lineno"> 400</span>  &sintau, &costau, &sxb, &cxb,</div><div class="line"><a name="l00401"></a><span class="lineno"> 401</span>  &scale2, &mu, &tanb, &vev2,</div><div class="line"><a name="l00402"></a><span class="lineno"> 402</span>  &result(0,0), &result(1,1), &result(0,1));</div><div class="line"><a name="l00403"></a><span class="lineno"> 403</span> </div><div class="line"><a name="l00404"></a><span class="lineno"> 404</span>  result(1,0) = result(0,1);</div><div class="line"><a name="l00405"></a><span class="lineno"> 405</span> </div><div class="line"><a name="l00406"></a><span class="lineno"> 406</span>  <span class="keyword">const</span> <span class="keywordtype">double</span> dMA = include_heavy_higgs</div><div class="line"><a name="l00407"></a><span class="lineno"> 407</span>  ? <a class="code" href="namespacehimalaya_1_1mh2l.html#a8b600e41354ff1538cbe86668852e468">delta_ma2_2loop_ab_atau</a>(</div><div class="line"><a name="l00408"></a><span class="lineno"> 408</span>  mtau2, mb2, mstau12, mstau22, msb12, msb22, sintau, costau, sxb, cxb,</div><div class="line"><a name="l00409"></a><span class="lineno"> 409</span>  scale2, mu, tanb, vev2)</div><div class="line"><a name="l00410"></a><span class="lineno"> 410</span>  : 0.;</div><div class="line"><a name="l00411"></a><span class="lineno"> 411</span> </div><div class="line"><a name="l00412"></a><span class="lineno"> 412</span>  <span class="keywordflow">return</span> result + rotate_by(dMA, tanb);</div><div class="line"><a name="l00413"></a><span class="lineno"> 413</span> }</div><div class="line"><a name="l00414"></a><span class="lineno"> 414</span> </div><div class="line"><a name="l00415"></a><span class="lineno"> 415</span> } <span class="comment">// namespace mh2l</span></div><div class="line"><a name="l00416"></a><span class="lineno"> 416</span> } <span class="comment">// namespace himalaya</span></div><div class="ttc" id="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth_html_ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df"><div class="ttname"><a href="namespacehimalaya_1_1hierarchies_1_1ExpansionDepth.html#ae39843d06a27e228f9fbea913e185c26aa43dcbe26742c23b9c4a1097777101df">himalaya::hierarchies::ExpansionDepth::At</a></div><div class="ttdoc">truncate the expansion depth in At/Ab by one order </div><div class="ttdef"><b>Definition:</b> <a href="HierarchyFlags_8hpp_source.html#l00046">HierarchyFlags.hpp:46</a></div></div>
<div class="ttc" id="Himalaya__LibraryLink_8m_html_a80c72e70807fc6ba175eba2e03353da8"><div class="ttname"><a href="Himalaya__LibraryLink_8m.html#a80c72e70807fc6ba175eba2e03353da8">a</a></div><div class="ttdeci">HimalayaNumericQ [a_] a</div><div class="ttdef"><b>Definition:</b> <a href="Himalaya__LibraryLink_8m_source.html#l00220">Himalaya_LibraryLink.m:220</a></div></div>
<div class="ttc" id="namespacehimalaya_html"><div class="ttname"><a href="namespacehimalaya.html">himalaya</a></div><div class="ttdef"><b>Definition:</b> <a href="H3_8cpp_source.html#l00014">H3.cpp:14</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html_ae6e9056d395834b500d50e3219a78fdb"><div class="ttname"><a href="DSZHiggs_8h.html#ae6e9056d395834b500d50e3219a78fdb">dszodd_</a></div><div class="ttdeci">int dszodd_(double *t, double *g, double *T1, double *T2, double *st, double *ct, double *q, double *mu, double *tanb, double *vv, double *gs, double *dma)</div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_accc07f22b2d75d17656ea35acd713c2b"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#accc07f22b2d75d17656ea35acd713c2b">himalaya::mh2l::delta_mh2_2loop_at_at</a></div><div class="ttdeci">Eigen::Matrix< double, 2, 2 > delta_mh2_2loop_at_at(double mt2, double mb2, double mA2, double mst12, double mst22, double msb12, double msb22, double sxt, double cxt, double sxb, double cxb, double scale2, double mu, double tanb, double vev2, int include_heavy_higgs, int atat)</div><div class="ttdoc">2-loop CP-even Higgs contribution O(at^2) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00293">DSZHiggs.cpp:293</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a5ddaca4b030e7354fd06bc149bbaf13a"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a5ddaca4b030e7354fd06bc149bbaf13a">himalaya::mh2l::delta_ma2_2loop_atau_atau</a></div><div class="ttdeci">double delta_ma2_2loop_atau_atau(double mtau2, double mA2, double msv2, double mstau12, double mstau22, double sintau, double costau, double scale2, double mu, double tanb, double vev2)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00336">DSZHiggs.cpp:336</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a83dba77a0cf5c8e1e1f0dc89c94360b3"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a83dba77a0cf5c8e1e1f0dc89c94360b3">himalaya::mh2l::delta_mh2_2loop_ab_as</a></div><div class="ttdeci">Eigen::Matrix< double, 2, 2 > delta_mh2_2loop_ab_as(double mb2, double mg, double msb12, double msb22, double sxb, double cxb, double scale2, double mu, double cotb, double vev2, double g3, int include_heavy_higgs)</div><div class="ttdoc">2-loop CP-even Higgs contribution O(ab*as) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00321">DSZHiggs.cpp:321</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2__eft_1_1threshold__loop__functions_html_a20013baed5c82acb480c5697a992756e"><div class="ttname"><a href="namespacehimalaya_1_1mh2__eft_1_1threshold__loop__functions.html#a20013baed5c82acb480c5697a992756e">himalaya::mh2_eft::threshold_loop_functions::g</a></div><div class="ttdeci">double g(double r) noexcept</div><div class="ttdef"><b>Definition:</b> <a href="ThresholdLoopFunctions_8cpp_source.html#l00446">ThresholdLoopFunctions.cpp:446</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_accc846c9b031c5d40de467d16e2f5054"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#accc846c9b031c5d40de467d16e2f5054">himalaya::mh2l::delta_ma2_2loop_at_as_general</a></div><div class="ttdeci">double delta_ma2_2loop_at_as_general(double mt2, double mg, double mst12, double mst22, double sxt, double cxt, double scale2, double mu, double tanb, double vev2, double g3)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00217">DSZHiggs.cpp:217</a></div></div>
<div class="ttc" id="Li2f_8hpp_html"><div class="ttname"><a href="Li2f_8hpp.html">Li2f.hpp</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a7827d302c444d22f1b46e1c855566039"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a7827d302c444d22f1b46e1c855566039">himalaya::mh2l::delta_mh2_2loop_ab_atau</a></div><div class="ttdeci">Eigen::Matrix< double, 2, 2 > delta_mh2_2loop_ab_atau(double mtau2, double mb2, double mstau12, double mstau22, double msb12, double msb22, double sintau, double costau, double sxb, double cxb, double scale2, double mu, double tanb, double vev2, int include_heavy_higgs)</div><div class="ttdoc">2-loop CP-even Higgs contribution O(ab*atau) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00389">DSZHiggs.cpp:389</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html_acbca846367a332a10c3d4c91a8439ab8"><div class="ttname"><a href="DSZHiggs_8h.html#acbca846367a332a10c3d4c91a8439ab8">taubot_</a></div><div class="ttdeci">int taubot_(double *t, double *b, double *T1, double *T2, double *B1, double *B2, double *st, double *ct, double *sb, double *cb, double *q, double *mu, double *tanb, double *vv, double *S11, double *S22, double *S12)</div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a21677404a27a0054d6a0117db96a181d"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a21677404a27a0054d6a0117db96a181d">himalaya::mh2l::delta_mh2_2loop_atau_atau</a></div><div class="ttdeci">Eigen::Matrix< double, 2, 2 > delta_mh2_2loop_atau_atau(double mtau2, double mA2, double msv2, double mstau12, double mstau22, double sintau, double costau, double scale2, double mu, double tanb, double vev2, int include_heavy_higgs)</div><div class="ttdoc">2-loop CP-even Higgs contribution O(atau^2) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00349">DSZHiggs.cpp:349</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_acc1ac4c68ee565a17377005fc22fa312"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#acc1ac4c68ee565a17377005fc22fa312">himalaya::mh2l::delta_mh2_2loop_at_as</a></div><div class="ttdeci">Eigen::Matrix< double, 2, 2 > delta_mh2_2loop_at_as(double mt2, double mg, double mst12, double mst22, double sxt, double cxt, double scale2, double mu, double tanb, double vev2, double g3, int include_heavy_higgs)</div><div class="ttdoc">2-loop CP-even Higgs contribution O(at*as) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00269">DSZHiggs.cpp:269</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html_a62881cea97186f48509abe184f224dcb"><div class="ttname"><a href="DSZHiggs_8h.html#a62881cea97186f48509abe184f224dcb">ddshiggs_</a></div><div class="ttdeci">int ddshiggs_(double *t, double *b, double *A0, double *T1, double *T2, double *B1, double *B2, double *st, double *ct, double *sb, double *cb, double *q, double *mu, double *tanb, double *vv, double *S11, double *S12, double *S22, int *atat)</div></div>
<div class="ttc" id="CouplingOrders_8hpp_html"><div class="ttname"><a href="CouplingOrders_8hpp.html">CouplingOrders.hpp</a></div><div class="ttdoc">enum definitions </div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a8b600e41354ff1538cbe86668852e468"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a8b600e41354ff1538cbe86668852e468">himalaya::mh2l::delta_ma2_2loop_ab_atau</a></div><div class="ttdeci">double delta_ma2_2loop_ab_atau(double mtau2, double mb2, double mstau12, double mstau22, double msb12, double msb22, double sintau, double costau, double sxb, double cxb, double scale2, double mu, double tanb, double vev2)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00373">DSZHiggs.cpp:373</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html"><div class="ttname"><a href="DSZHiggs_8h.html">DSZHiggs.h</a></div><div class="ttdoc">Contains two-loop routines from Slavich et al. </div></div>
<div class="ttc" id="DSZHiggs_8h_html_a9f210e554b67f63390a67ff1b60540e2"><div class="ttname"><a href="DSZHiggs_8h.html#a9f210e554b67f63390a67ff1b60540e2">taubotodd_</a></div><div class="ttdeci">int taubotodd_(double *t, double *b, double *T1, double *T2, double *B1, double *B2, double *st, double *ct, double *sb, double *cb, double *q, double *mu, double *tanb, double *vv, double *dma)</div></div>
<div class="ttc" id="DSZHiggs_8h_html_a20219899e8065151504afdc38c83e481"><div class="ttname"><a href="DSZHiggs_8h.html#a20219899e8065151504afdc38c83e481">tausqodd_</a></div><div class="ttdeci">int tausqodd_(double *t, double *A0, double *BL, double *T1, double *T2, double *st, double *ct, double *q, double *mu, double *tanb, double *vv, double *dma)</div></div>
<div class="ttc" id="namespacehimalaya_html_a5747395de1ca13183c7e7abbff826a5e"><div class="ttname"><a href="namespacehimalaya.html#a5747395de1ca13183c7e7abbff826a5e">himalaya::li2</a></div><div class="ttdeci">double li2(double x)</div><div class="ttdoc">real dilogarithm from FORTRAN module </div><div class="ttdef"><b>Definition:</b> <a href="Li2f_8cpp_source.html#l00017">Li2f.cpp:17</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_acbdc52f7ab258a55f03b9f8bd069784e"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#acbdc52f7ab258a55f03b9f8bd069784e">himalaya::mh2l::delta_ma2_2loop_at_at</a></div><div class="ttdeci">double delta_ma2_2loop_at_at(double mt2, double mb2, double mA2, double mst12, double mst22, double msb12, double msb22, double sxt, double cxt, double sxb, double cxb, double scale2, double mu, double tanb, double vev2, int atat)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00251">DSZHiggs.cpp:251</a></div></div>
<div class="ttc" id="DSZHiggs_8hpp_html"><div class="ttname"><a href="DSZHiggs_8hpp.html">DSZHiggs.hpp</a></div><div class="ttdoc">function declarations for 2-loop MSSM Higgs contributions </div></div>
<div class="ttc" id="namespacehimalaya_html_a99d07fbbb2c5cc71320a23984dc65b58"><div class="ttname"><a href="namespacehimalaya.html#a99d07fbbb2c5cc71320a23984dc65b58">himalaya::log</a></div><div class="ttdeci">Complex< T > log(const Complex< T > &z_) noexcept</div><div class="ttdef"><b>Definition:</b> <a href="complex_8hpp_source.html#l00045">complex.hpp:45</a></div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_ad0ce42f717a9918763604fd55fcb4e3d"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#ad0ce42f717a9918763604fd55fcb4e3d">himalaya::mh2l::delta_ma2_2loop_at_as_mst1_eq_mst2</a></div><div class="ttdeci">double delta_ma2_2loop_at_as_mst1_eq_mst2(double mt2, double mg, double mst12, double mst22, double sxt, double cxt, double scale2, double mu, double tanb, double vev2, double g3)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00188">DSZHiggs.cpp:188</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html_afebc9638df87be653b7998cd246cc4a3"><div class="ttname"><a href="DSZHiggs_8h.html#afebc9638df87be653b7998cd246cc4a3">tausqhiggs_</a></div><div class="ttdeci">int tausqhiggs_(double *t, double *A0, double *BL, double *T1, double *T2, double *st, double *ct, double *q, double *mu, double *tanb, double *vv, int *OS, double *S11, double *S22, double *S12)</div></div>
<div class="ttc" id="namespacehimalaya_1_1mh2l_html_a9561889cc35e34a709643cfed97856cc"><div class="ttname"><a href="namespacehimalaya_1_1mh2l.html#a9561889cc35e34a709643cfed97856cc">himalaya::mh2l::delta_ma2_2loop_at_as</a></div><div class="ttdeci">double delta_ma2_2loop_at_as(double mt2, double mg, double mst12, double mst22, double sxt, double cxt, double scale2, double mu, double tanb, double vev2, double g3)</div><div class="ttdoc">2-loop contribution to CP-odd Higgs O(at*as) </div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8cpp_source.html#l00231">DSZHiggs.cpp:231</a></div></div>
<div class="ttc" id="DSZHiggs_8h_html_accd0e91b067d839181665628d66d932f"><div class="ttname"><a href="DSZHiggs_8h.html#accd0e91b067d839181665628d66d932f">ddsodd_</a></div><div class="ttdeci">int ddsodd_(double *t, double *b, double *A0, double *T1, double *T2, double *B1, double *B2, double *st, double *ct, double *sb, double *cb, double *q, double *mu, double *tanb, double *vv, double *dma, int *atat)</div></div>
<div class="ttc" id="DSZHiggs_8h_html_a88e8c8f0df22394855df1dfc61564c2d"><div class="ttname"><a href="DSZHiggs_8h.html#a88e8c8f0df22394855df1dfc61564c2d">dszhiggs_</a></div><div class="ttdeci">int dszhiggs_(double *t, double *g, double *T1, double *T2, double *st, double *ct, double *q, double *mu, double *tanb, double *vv, double *gs, int *OS, double *S11, double *S22, double *S12)</div></div>
<div class="ttc" id="DSZHiggs_8f_html_a6120448a5a0919fd26c80580dce5c0f3"><div class="ttname"><a href="DSZHiggs_8f.html#a6120448a5a0919fd26c80580dce5c0f3">phi</a></div><div class="ttdeci">double precision function phi(x, y, z)</div><div class="ttdef"><b>Definition:</b> <a href="DSZHiggs_8f_source.html#l05161">DSZHiggs.f:5161</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated by  <a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.8.13
</small></address>
</body>
</html>