Hyperiso 1.0.3
Modular flavour-physics calculations, Wilson coefficients and statistical inference
Loading...
Searching...
No Matches
special_THDM.cpp
Go to the documentation of this file.
1#include <cmath>
2#include "special_THDM.h"
3
4
5double F1SP(double xt, double xH)
6{
7 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F1SP(0.9998,1.0002);
8
9 if(fabs(1.-xt)<1.e-5) return F1SP(0.9999,xH);
10 if(fabs(1.-xH)<1.e-5) return F1SP(xt,0.9999);
11
12 if(fabs(1.-xt/xH)<1.e-5) return F1SP(xH*0.9998,xH);
13
14 return 1./4./(xH-xt)*(xt*log(xt)/(xt-1.)-xH*log(xH)/(xH-1.));
15}
16
17
18/*--------------------------------------------------------------------*/
19
20double F2SP(double xt, double xH)
21{
22 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F2SP(0.9998,1.0002);
23
24 if(fabs(1.-xt)<1.e-5) return F2SP(0.9999,xH);
25 if(fabs(1.-xH)<1.e-5) return F2SP(xt,0.9999);
26
27 if(fabs(1.-xt/xH)<1.e-5) return F2SP(xH*0.9998,xH);
28
29 return 1./8./(xH-xt)*(xH/(xH-1.)+xt*xt*log(xt)/(xt-1.)/(xH-xt)-xH*(xH*xt+xH-2.*xt)/(xH-1.)/(xH-1.)/(xH-xt)*log(xH));
30}
31
32/*--------------------------------------------------------------------*/
33
34double F3SP(double xt, double xH)
35{
36 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F3SP(0.9998,1.0002);
37
38 if(fabs(1.-xt)<1.e-5) return F3SP(0.9999,xH);
39 if(fabs(1.-xH)<1.e-5) return F3SP(xt,0.9999);
40
41 if(fabs(1.-xt/xH)<1.e-5) return F3SP(xH*0.9998,xH);
42
43 return 1./8./(xH-xt)*((xH-xt)/(xH-1.)/(xt-1.)+xt*(xt-2.)*log(xt)/(xt-1.)/(xt-1.)-xH*(xH-2.)/(xH-1.)/(xH-1.)*log(xH));
44}
45
46/*--------------------------------------------------------------------*/
47
48double F4SP(double xt, double xH)
49{
50 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F4SP(0.9998,1.0002);
51
52 if(fabs(1.-xt)<1.e-5) return F4SP(0.9999,xH);
53 if(fabs(1.-xH)<1.e-5) return F4SP(xt,0.9999);
54
55 if(fabs(1.-xt/xH)<1.e-5) return F4SP(xH*0.9998,xH);
56
57 return xt/(xH-xt)*(1.-xH/(xH-xt)*log(xH/xt));
58}
59
60/*--------------------------------------------------------------------*/
61
62double F5SP(double xt, double xH)
63{
64 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F5SP(0.9998,1.0002);
65
66 if(fabs(1.-xt)<1.e-5) return F5SP(0.9999,xH);
67 if(fabs(1.-xH)<1.e-5) return F5SP(xt,0.9999);
68
69 if(fabs(1.-xt/xH)<1.e-5) return F5SP(xH*0.9998,xH);
70
71 return xt/2./(xH-xt)/(xH-xt)*((xH+xt)/2.-xH*xt/(xH-xt)*log(xH/xt));
72}
73
74/*--------------------------------------------------------------------*/
75
76double F6SP(double xt, double xH)
77{
78 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F6SP(0.9998,1.0002);
79
80 if(fabs(1.-xt)<1.e-5) return F6SP(0.9999,xH);
81 if(fabs(1.-xH)<1.e-5) return F6SP(xt,0.9999);
82
83 if(fabs(1.-xt/xH)<1.e-5) return F6SP(xH*0.9998,xH);
84
85 return 1./2./(xH-xt)*(-xH+xt+xH*log(xH)-xt*log(xt));
86}
87
88/*--------------------------------------------------------------------*/
89
90double F7SP(double xt, double xH)
91{
92 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F7SP(0.9998,1.0002);
93
94 if(fabs(1.-xt)<1.e-5) return F7SP(0.9999,xH);
95 if(fabs(1.-xH)<1.e-5) return F7SP(xt,0.9999);
96
97 if(fabs(1.-xt/xH)<1.e-5) return F7SP(xH*0.9998,xH);
98
99 return 1./2./(xH-xt)*(xt-xH*xt/(xH-xt)*log(xH/xt));
100}
101
102/*--------------------------------------------------------------------*/
103
104double F8SP(double xt, double xH)
105{
106 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F8SP(0.9998,1.0002);
107
108 if(fabs(1.-xt)<1.e-5) return F8SP(0.9999,xH);
109 if(fabs(1.-xH)<1.e-5) return F8SP(xt,0.9999);
110
111 if(fabs(1.-xt/xH)<1.e-5) return F8SP(xH*0.9998,xH);
112
113 return 1./2./(xH-xt)*(xH-xH*xH*log(xH)/(xH-xt)+xt*(2.*xH-xt)*log(xt)/(xH-xt));
114}
115
116/*--------------------------------------------------------------------*/
117
118double F9SP(double xt, double xH)
119{
120 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F9SP(0.9998,1.0002);
121
122 if(fabs(1.-xt)<1.e-5) return F9SP(0.9999,xH);
123 if(fabs(1.-xH)<1.e-5) return F9SP(xt,0.9999);
124
125 if(fabs(1.-xt/xH)<1.e-5) return F9SP(xH*0.9998,xH);
126
127 return 1./4./(xH-xt)/(xH-xt)*(xt*(3.*xH-xt)/2.-xH*xH*xt/(xH-xt)*log(xH/xt));
128}
129
130/*--------------------------------------------------------------------*/
131
132double F10SP(double xt, double xH)
133{
134 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F10SP(0.9998,1.0002);
135
136 if(fabs(1.-xt)<1.e-5) return F10SP(0.9999,xH);
137 if(fabs(1.-xH)<1.e-5) return F10SP(xt,0.9999);
138
139 if(fabs(1.-xt/xH)<1.e-5) return F10SP(xH*0.9998,xH);
140
141 return 1./4./(xH-xt)/(xH-xt)*(xt*(xH-3.*xt)/2.-xH*xt*(xH-2.*xt)/(xH-xt)*log(xH/xt));
142}
143
144/*--------------------------------------------------------------------*/
145
146double F11SP(double xt, double xH)
147{
148 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F11SP(0.9998,1.0002);
149
150 if(fabs(1.-xt)<1.e-5) return F11SP(0.9999,xH);
151 if(fabs(1.-xH)<1.e-5) return F11SP(xt,0.9999);
152
153 if(fabs(1.-xt/xH)<1.e-5) return F11SP(xH*0.9998,xH);
154
155 double xt2=xt*xt;
156 double xt3=xt*xt2;
157
158 return 1./2./(xH-xt)*(xt*(xt2-3.*xH*xt+9.*xH-5.*xt-2.)/4./(xt-1.)/(xt-1.)+xH*(xH*xt-3.*xH+2.*xt)/2./(xH-1.)/(xH-xt)*log(xH)
159 +(xH*xH*(-2.*xt3+6.*xt2-9.*xt+2.)+3.*xH*xt2*(xt2-2.*xt+3.)-xt2*(2.*xt3-2.*xt2+3.*xt+1.))*log(xt)/2./pow(xt-1.,3.)/(xH-xt));
160}
161
162/*--------------------------------------------------------------------*/
163
164double F12SP(double xt, double xH)
165{
166 if((fabs(1.-xt)<1.e-5)&&(fabs(1.-xH)<1.e-5)) return F12SP(0.9998,1.0002);
167
168 if(fabs(1.-xt)<1.e-5) return F12SP(0.9999,xH);
169 if(fabs(1.-xH)<1.e-5) return F12SP(xt,0.9999);
170
171 if(fabs(1.-xt/xH)<1.e-5) return F12SP(xH*0.9998,xH);
172
173 double xt2=xt*xt;
174 double xt3=xt*xt2;
175
176 return 1./2./(xH-xt)*((xt2+xt-8.)*(xH-xt)/4./(xt-1.)/(xt-1.)-xH*(xH+2.)*log(xH)/2./(xH-1.)
177 +(xH*(xt3-3.*xt2+3.*xt+2.)+3.*(xt-2.)*xt2)*log(xt)/2./pow(xt-1.,3.));
178}
179
180/*--------------------------------------------------------------------*/
181
182double CSc_2HDM(double xH, double xt, double lu, double ld, double ll)
183{
184 if((fabs(1.-xH)<1.e-5)&&(fabs(1.-xt)<1.e-5)) return CSc_2HDM(0.9998,1.0002,lu,ld,ll);
185
186 if(fabs(1.-xH)<1.e-5) return CSc_2HDM(0.9999,xt,lu,ld,ll);
187 if(fabs(1.-xt)<1.e-5) return CSc_2HDM(xH,0.9999,lu,ld,ll);
188
189 if(fabs(1.-xH/xt)<1.e-5) return CSc_2HDM(xt*0.9998,xt,lu,ld,ll);
190
191 return xt/8./(xH-xt)*(2.*ld*ll*(1./(xH-1.)*log(xH)-1./(xt-1.)*log(xt))
192 +lu*ll*(1./(xH-1.)+xH/(xH-xt)/(xt-1.)*log(xt)-xH*(2.*xH-xt-1.)/(xH-xt)/(xH-1.)/(xH-1.)*log(xH))
193 -ll*lu*((xt-xH)/(xH-1.)/(xt-1.)+xt/(xt-1.)/(xt-1.)*log(xt)-xH/(xH-1.)/(xH-1.)*log(xH)));
194}
195
196/*--------------------------------------------------------------------*/
197
198double CPc_2HDM(double xH, double xt, double lu, double ld, double ll, double sw2)
199{
200 if((fabs(1.-xH)<1.e-5)&&(fabs(1.-xt)<1.e-5)) return CPc_2HDM(0.9998,1.0002,lu,ld,ll,sw2);
201
202 if(fabs(1.-xH)<1.e-5) return CPc_2HDM(0.9999,xt,lu,ld,ll,sw2);
203 if(fabs(1.-xt)<1.e-5) return CPc_2HDM(xH,0.9999,lu,ld,ll,sw2);
204
205 if(fabs(1.-xH/xt)<1.e-5) return CPc_2HDM(xt*0.9998,xt,lu,ld,ll,sw2);
206
207 return -xt/8./(xH-xt)*(2.*ld*ll*(1./(xH-1.)*log(xH)-1./(xt-1.)*log(xt))
208 +lu*ll*(1./(xH-1.)+xH/(xH-xt)/(xt-1.)*log(xt)-xH*(2.*xH-xt-1.)/(xH-xt)/(xH-1.)/(xH-1.)*log(xH))
209 +ll*lu*((xt-xH)/(xH-1.)/(xt-1.)+xt/(xt-1.)/(xt-1.)*log(xt)-xH/(xH-1.)/(xH-1.)*log(xH)))
210 +xt/4./(xH-xt)/(xH-xt)*(-ld*lu*(-(xt+xH)/2.+xt*xH/(xH-xt)*log(xH/xt))
211 +lu*lu/6./(xH-xt)*((xH*xH-8.*xH*xt-17.*xt*xt)/6.+xt*xt*(3.*xH+xt)/(xH-xt)*log(xH/xt)))
212 +sw2*xt/6./(xH-xt)/(xH-xt)*(-ld*lu*((5.*xt-3.*xH)/2.+xH*(2.*xH-3.*xt)/(xH-xt)*log(xH/xt))
213 +lu*lu/6./(xH-xt)*((4.*xH*xH*xH-12.*xH*xH*xt+xH*xt*xt+3.*xt*xt*xt)/(xH-xt)*log(xH/xt)-(17.*xH*xH-64.*xH*xt+71.*xt*xt)/6.))
214 +(1.-sw2)*lu*lu*xt*xt/4./(xH-xt)/(xH-xt)*(xH*log(xH/xt)+xt-xH);
215}
pow
scalar_t pow(const scalar_t &base, const scalar_t &exp)
Definition scalar.cpp:75
F6SP
double F6SP(double xt, double xH)
Definition special_THDM.cpp:76
CSc_2HDM
double CSc_2HDM(double xH, double xt, double lu, double ld, double ll)
Definition special_THDM.cpp:182
F7SP
double F7SP(double xt, double xH)
Definition special_THDM.cpp:90
F1SP
double F1SP(double xt, double xH)
Definition special_THDM.cpp:5
F10SP
double F10SP(double xt, double xH)
Definition special_THDM.cpp:132
F5SP
double F5SP(double xt, double xH)
Definition special_THDM.cpp:62
F4SP
double F4SP(double xt, double xH)
Definition special_THDM.cpp:48
F12SP
double F12SP(double xt, double xH)
Definition special_THDM.cpp:164
CPc_2HDM
double CPc_2HDM(double xH, double xt, double lu, double ld, double ll, double sw2)
Definition special_THDM.cpp:198
F3SP
double F3SP(double xt, double xH)
Definition special_THDM.cpp:34
F9SP
double F9SP(double xt, double xH)
Definition special_THDM.cpp:118
F8SP
double F8SP(double xt, double xH)
Definition special_THDM.cpp:104
F2SP
double F2SP(double xt, double xH)
Definition special_THDM.cpp:20
F11SP
double F11SP(double xt, double xH)
Definition special_THDM.cpp:146
special_THDM.h