-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathTHDM_2gm2__1loop_8cpp_source.html
162 lines (160 loc) · 18.9 KB
/
THDM_2gm2__1loop_8cpp_source.html
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
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en-US">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=11"/>
<meta name="generator" content="Doxygen 1.9.8"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>GM2Calc: /home/avoigt/research/GM2Calc/src/THDM/gm2_1loop.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 id="projectrow">
<td id="projectalign">
<div id="projectname">GM2Calc<span id="projectnumber"> 2.3.0</span>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.9.8 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
var searchBox = new SearchBox("searchBox", "search/",'.html');
/* @license-end */
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
$(function() {
initMenu('',true,false,'search.php','Search');
$(document).ready(function() { init_search(); });
});
/* @license-end */
</script>
<div id="main-nav"></div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:d3d9a9a6595521f9666a5e94cc830dab83b65699&dn=expat.txt MIT */
$(document).ready(function() { init_codefold(0); });
/* @license-end */
</script>
<!-- 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">
<div id="MSearchResults">
<div class="SRPage">
<div id="SRIndex">
<div id="SRResults"></div>
<div class="SRStatus" id="Loading">Loading...</div>
<div class="SRStatus" id="Searching">Searching...</div>
<div class="SRStatus" id="NoMatches">No Matches</div>
</div>
</div>
</div>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_68267d1309a1af8e8297ef4c3efbcdba.html">src</a></li><li class="navelem"><a class="el" href="dir_1672695a47b377b4e14cff0be868b5f0.html">THDM</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle"><div class="title">gm2_1loop.cpp</div></div>
</div><!--header-->
<div class="contents">
<a href="THDM_2gm2__1loop_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a id="l00001" name="l00001"></a><span class="lineno"> 1</span><span class="comment">// ====================================================================</span></div>
<div class="line"><a id="l00002" name="l00002"></a><span class="lineno"> 2</span><span class="comment">// This file is part of GM2Calc.</span></div>
<div class="line"><a id="l00003" name="l00003"></a><span class="lineno"> 3</span><span class="comment">//</span></div>
<div class="line"><a id="l00004" name="l00004"></a><span class="lineno"> 4</span><span class="comment">// GM2Calc is free software: you can redistribute it and/or modify</span></div>
<div class="line"><a id="l00005" name="l00005"></a><span class="lineno"> 5</span><span class="comment">// it under the terms of the GNU General Public License as published</span></div>
<div class="line"><a id="l00006" name="l00006"></a><span class="lineno"> 6</span><span class="comment">// by the Free Software Foundation, either version 3 of the License,</span></div>
<div class="line"><a id="l00007" name="l00007"></a><span class="lineno"> 7</span><span class="comment">// or (at your option) any later version.</span></div>
<div class="line"><a id="l00008" name="l00008"></a><span class="lineno"> 8</span><span class="comment">//</span></div>
<div class="line"><a id="l00009" name="l00009"></a><span class="lineno"> 9</span><span class="comment">// GM2Calc is distributed in the hope that it will be useful, but</span></div>
<div class="line"><a id="l00010" name="l00010"></a><span class="lineno"> 10</span><span class="comment">// WITHOUT ANY WARRANTY; without even the implied warranty of</span></div>
<div class="line"><a id="l00011" name="l00011"></a><span class="lineno"> 11</span><span class="comment">// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU</span></div>
<div class="line"><a id="l00012" name="l00012"></a><span class="lineno"> 12</span><span class="comment">// General Public License for more details.</span></div>
<div class="line"><a id="l00013" name="l00013"></a><span class="lineno"> 13</span><span class="comment">//</span></div>
<div class="line"><a id="l00014" name="l00014"></a><span class="lineno"> 14</span><span class="comment">// You should have received a copy of the GNU General Public License</span></div>
<div class="line"><a id="l00015" name="l00015"></a><span class="lineno"> 15</span><span class="comment">// along with GM2Calc. If not, see</span></div>
<div class="line"><a id="l00016" name="l00016"></a><span class="lineno"> 16</span><span class="comment">// <http://www.gnu.org/licenses/>.</span></div>
<div class="line"><a id="l00017" name="l00017"></a><span class="lineno"> 17</span><span class="comment">// ====================================================================</span></div>
<div class="line"><a id="l00018" name="l00018"></a><span class="lineno"> 18</span> </div>
<div class="line"><a id="l00019" name="l00019"></a><span class="lineno"> 19</span><span class="preprocessor">#include "<a class="code" href="gm2__1loop_8hpp.html">gm2calc/gm2_1loop.hpp</a>"</span></div>
<div class="line"><a id="l00020" name="l00020"></a><span class="lineno"> 20</span><span class="preprocessor">#include "<a class="code" href="THDM_8hpp.html">gm2calc/THDM.hpp</a>"</span></div>
<div class="line"><a id="l00021" name="l00021"></a><span class="lineno"> 21</span><span class="preprocessor">#include "<a class="code" href="THDM_2gm2__1loop__helpers_8hpp.html">THDM/gm2_1loop_helpers.hpp</a>"</span></div>
<div class="line"><a id="l00022" name="l00022"></a><span class="lineno"> 22</span><span class="preprocessor">#include <cmath></span></div>
<div class="line"><a id="l00023" name="l00023"></a><span class="lineno"> 23</span><span class="comment"></span> </div>
<div class="line"><a id="l00024" name="l00024"></a><span class="lineno"> 24</span><span class="comment">/**</span></div>
<div class="line"><a id="l00025" name="l00025"></a><span class="lineno"> 25</span><span class="comment"> * \file gm2_1loop.cpp</span></div>
<div class="line"><a id="l00026" name="l00026"></a><span class="lineno"> 26</span><span class="comment"> *</span></div>
<div class="line"><a id="l00027" name="l00027"></a><span class="lineno"> 27</span><span class="comment"> * Contains functions necessary to calculate the THDM</span></div>
<div class="line"><a id="l00028" name="l00028"></a><span class="lineno"> 28</span><span class="comment"> * contributions for g-2 at the 1-loop level.</span></div>
<div class="line"><a id="l00029" name="l00029"></a><span class="lineno"> 29</span><span class="comment"> */</span></div>
<div class="line"><a id="l00030" name="l00030"></a><span class="lineno"> 30</span> </div>
<div class="line"><a id="l00031" name="l00031"></a><span class="lineno"> 31</span><span class="keyword">namespace </span><a class="code hl_namespace" href="namespacegm2calc.html">gm2calc</a> {</div>
<div class="line"><a id="l00032" name="l00032"></a><span class="lineno"> 32</span><span class="comment"></span> </div>
<div class="line"><a id="l00033" name="l00033"></a><span class="lineno"> 33</span><span class="comment">/**</span></div>
<div class="line"><a id="l00034" name="l00034"></a><span class="lineno"> 34</span><span class="comment"> * Calculates full 1-loop contribution to a_mu in the general THDM.</span></div>
<div class="line"><a id="l00035" name="l00035"></a><span class="lineno"> 35</span><span class="comment"> *</span></div>
<div class="line"><a id="l00036" name="l00036"></a><span class="lineno"> 36</span><span class="comment"> * @param model THDM model parameters, masses and mixings</span></div>
<div class="line"><a id="l00037" name="l00037"></a><span class="lineno"> 37</span><span class="comment"> * @return 1-loop contribution to a_mu</span></div>
<div class="line"><a id="l00038" name="l00038"></a><span class="lineno"> 38</span><span class="comment"> */</span></div>
<div class="foldopen" id="foldopen00039" data-start="{" data-end="}">
<div class="line"><a id="l00039" name="l00039"></a><span class="lineno"><a class="line" href="namespacegm2calc.html#a3116e8571a608a48bfb4bb49595ff662"> 39</a></span><span class="keywordtype">double</span> <a class="code hl_function" href="namespacegm2calc.html#a9969e35b9b08b4d8570bdf76a4c40586">calculate_amu_1loop</a>(<span class="keyword">const</span> <a class="code hl_class" href="classgm2calc_1_1THDM.html">THDM</a>& <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>)</div>
<div class="line"><a id="l00040" name="l00040"></a><span class="lineno"> 40</span>{</div>
<div class="line"><a id="l00041" name="l00041"></a><span class="lineno"> 41</span> <a class="code hl_struct" href="structgm2calc_1_1thdm_1_1THDM__1L__parameters.html">thdm::THDM_1L_parameters</a> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>;</div>
<div class="line"><a id="l00042" name="l00042"></a><span class="lineno"> 42</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.<a class="code hl_variable" href="structgm2calc_1_1thdm_1_1THDM__1L__parameters.html#afa73bdda97d4ae3bbc8fc00d761b91d0">alpha_em</a> = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_alpha_em();</div>
<div class="line"><a id="l00043" name="l00043"></a><span class="lineno"> 43</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mm = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MFe(1);</div>
<div class="line"><a id="l00044" name="l00044"></a><span class="lineno"> 44</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mw = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MVWm();</div>
<div class="line"><a id="l00045" name="l00045"></a><span class="lineno"> 45</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mz = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MVZ();</div>
<div class="line"><a id="l00046" name="l00046"></a><span class="lineno"> 46</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mhSM = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_sm().get_mh();</div>
<div class="line"><a id="l00047" name="l00047"></a><span class="lineno"> 47</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mA = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MAh(1);</div>
<div class="line"><a id="l00048" name="l00048"></a><span class="lineno"> 48</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mHp = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MHm(1);</div>
<div class="line"><a id="l00049" name="l00049"></a><span class="lineno"> 49</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.ml = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MFe();</div>
<div class="line"><a id="l00050" name="l00050"></a><span class="lineno"> 50</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mv = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_MFv();</div>
<div class="line"><a id="l00051" name="l00051"></a><span class="lineno"> 51</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.mh = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_Mhh();</div>
<div class="line"><a id="l00052" name="l00052"></a><span class="lineno"> 52</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.ylh = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_ylh();</div>
<div class="line"><a id="l00053" name="l00053"></a><span class="lineno"> 53</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.ylH = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_ylH();</div>
<div class="line"><a id="l00054" name="l00054"></a><span class="lineno"> 54</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.ylA = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_ylA();</div>
<div class="line"><a id="l00055" name="l00055"></a><span class="lineno"> 55</span> <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>.ylHp = <a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">model</a>.get_ylHp();</div>
<div class="line"><a id="l00056" name="l00056"></a><span class="lineno"> 56</span> </div>
<div class="line"><a id="l00057" name="l00057"></a><span class="lineno"> 57</span> <span class="keywordflow">return</span> <a class="code hl_function" href="namespacegm2calc_1_1thdm.html#a85aa785d56e074ea0827bc9014186e7d">thdm::amu1L</a>(<a class="code hl_function" href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">pars</a>);</div>
<div class="line"><a id="l00058" name="l00058"></a><span class="lineno"> 58</span>}</div>
</div>
<div class="line"><a id="l00059" name="l00059"></a><span class="lineno"> 59</span> </div>
<div class="line"><a id="l00060" name="l00060"></a><span class="lineno"> 60</span>} <span class="comment">// namespace gm2calc</span></div>
<div class="ttc" id="aTHDM_2gm2__1loop__helpers_8hpp_html"><div class="ttname"><a href="THDM_2gm2__1loop__helpers_8hpp.html">gm2_1loop_helpers.hpp</a></div></div>
<div class="ttc" id="aTHDM_8hpp_html"><div class="ttname"><a href="THDM_8hpp.html">THDM.hpp</a></div></div>
<div class="ttc" id="aclassgm2calc_1_1THDM_html"><div class="ttname"><a href="classgm2calc_1_1THDM.html">gm2calc::THDM</a></div><div class="ttdoc">Contains routines to determine the THDM parameters.</div><div class="ttdef"><b>Definition</b> <a href="THDM_8hpp_source.html#l00098">THDM.hpp:98</a></div></div>
<div class="ttc" id="agm2__1loop_8hpp_html"><div class="ttname"><a href="gm2__1loop_8hpp.html">gm2_1loop.hpp</a></div></div>
<div class="ttc" id="anamespacegm2calc_1_1thdm_html_a85aa785d56e074ea0827bc9014186e7d"><div class="ttname"><a href="namespacegm2calc_1_1thdm.html#a85aa785d56e074ea0827bc9014186e7d">gm2calc::thdm::amu1L</a></div><div class="ttdeci">double amu1L(const THDM_1L_parameters &pars) noexcept</div><div class="ttdoc">Full (CP-conserving) 1-loop contribution.</div><div class="ttdef"><b>Definition</b> <a href="gm2__1loop__H_8cpp_source.html#l00128">gm2_1loop_H.cpp:128</a></div></div>
<div class="ttc" id="anamespacegm2calc_html"><div class="ttname"><a href="namespacegm2calc.html">gm2calc</a></div><div class="ttdef"><b>Definition</b> <a href="gm2__config__options_8hpp_source.html#l00022">gm2_config_options.hpp:22</a></div></div>
<div class="ttc" id="anamespacegm2calc_html_a5cc5ec26f0e3675025a5edf971e9a84b"><div class="ttname"><a href="namespacegm2calc.html#a5cc5ec26f0e3675025a5edf971e9a84b">gm2calc::svd_eigen</a></div><div class="ttdeci">void svd_eigen(const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real,(((M)<(N)) ?(M) :(N)), 1 > &s, Eigen::Matrix< Scalar, M, M > *u, Eigen::Matrix< Scalar, N, N > *vh)</div><div class="ttdef"><b>Definition</b> <a href="gm2__linalg_8hpp_source.html#l00039">gm2_linalg.hpp:40</a></div></div>
<div class="ttc" id="anamespacegm2calc_html_a9969e35b9b08b4d8570bdf76a4c40586"><div class="ttname"><a href="namespacegm2calc.html#a9969e35b9b08b4d8570bdf76a4c40586">gm2calc::calculate_amu_1loop</a></div><div class="ttdeci">double calculate_amu_1loop(const MSSMNoFV_onshell &model)</div><div class="ttdoc">Calculates full 1-loop SUSY contribution to (g-2), Eq (45) of arXiv:hep-ph/0609168.</div><div class="ttdef"><b>Definition</b> <a href="MSSMNoFV_2gm2__1loop_8cpp_source.html#l00074">gm2_1loop.cpp:74</a></div></div>
<div class="ttc" id="astructgm2calc_1_1thdm_1_1THDM__1L__parameters_html"><div class="ttname"><a href="structgm2calc_1_1thdm_1_1THDM__1L__parameters.html">gm2calc::thdm::THDM_1L_parameters</a></div><div class="ttdoc">parameters to be passed to the 1-loop contribution functions</div><div class="ttdef"><b>Definition</b> <a href="THDM_2gm2__1loop__helpers_8hpp_source.html#l00029">gm2_1loop_helpers.hpp:29</a></div></div>
<div class="ttc" id="astructgm2calc_1_1thdm_1_1THDM__1L__parameters_html_afa73bdda97d4ae3bbc8fc00d761b91d0"><div class="ttname"><a href="structgm2calc_1_1thdm_1_1THDM__1L__parameters.html#afa73bdda97d4ae3bbc8fc00d761b91d0">gm2calc::thdm::THDM_1L_parameters::alpha_em</a></div><div class="ttdeci">double alpha_em</div><div class="ttdoc">electromagnetic coupling</div><div class="ttdef"><b>Definition</b> <a href="THDM_2gm2__1loop__helpers_8hpp_source.html#l00030">gm2_1loop_helpers.hpp:30</a></div></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated by <a href="https://www.doxygen.org/index.html"><img class="footer" src="doxygen.svg" width="104" height="31" alt="doxygen"/></a> 1.9.8
</small></address>
</body>
</html>