a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
THDMWpositivity.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2017 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "THDMWpositivity.h"
9#include "THDMW.h"
10#include "THDMWcache.h"
11
12THDMWpositivity1::THDMWpositivity1(const StandardModel& SM_i)
13: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
14{}
15
16double THDMWpositivity1::computeThValue()
17{
18 double mu1=myTHDMW.getTHDMW_mu1();
19 double mu2=myTHDMW.getTHDMW_mu2();
20 double mu3=myTHDMW.getTHDMW_mu3();
21 double mu4=myTHDMW.getTHDMW_mu4();
22 double mu5=myTHDMW.getTHDMW_mu5();
23 double mu6=myTHDMW.getTHDMW_mu6();
24
25 return mu1+mu2+2.0*(mu3+mu4+mu5)+mu6;
26}
27
28
29
30THDMWpositivity2::THDMWpositivity2(const StandardModel& SM_i)
31: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
32{}
33
34double THDMWpositivity2::computeThValue()
35{
36 double mu1=myTHDMW.getTHDMW_mu1();
37 double mu2=myTHDMW.getTHDMW_mu2();
38 double mu3=myTHDMW.getTHDMW_mu3();
39 double mu4=myTHDMW.getTHDMW_mu4();
40
41 return mu1+mu2+mu3+mu4;
42}
43
44
45
46THDMWpositivity3::THDMWpositivity3(const StandardModel& SM_i)
47: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
48{}
49
50double THDMWpositivity3::computeThValue()
51{
52 double mu1=myTHDMW.getTHDMW_mu1();
53 double mu2=myTHDMW.getTHDMW_mu2();
54 double mu3=myTHDMW.getTHDMW_mu3();
55 double mu4=myTHDMW.getTHDMW_mu4();
56 double mu6=myTHDMW.getTHDMW_mu6();
57
58 return 14.0*(mu1+mu2) + 24.0*(mu3+mu4) + 5.0*mu6 - 3.0*fabs(2.0*(mu1+mu2)-mu6);
59}
60
61
62
63THDMWpositivity4::THDMWpositivity4(const StandardModel& SM_i)
64: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
65{}
66
67double THDMWpositivity4::computeThValue()
68{
69 double mu1=myTHDMW.getTHDMW_mu1();
70 double mu2=myTHDMW.getTHDMW_mu2();
71 double mu3=myTHDMW.getTHDMW_mu3();
72 double mu4=myTHDMW.getTHDMW_mu4();
73 double mu5=myTHDMW.getTHDMW_mu5();
74 double mu6=myTHDMW.getTHDMW_mu6();
75
76 return 5.0*(mu1+mu2+mu6) + 6.0*(2.0*mu3+mu4+mu5) - fabs(mu1+mu2+mu6);
77}
78
79
80
81THDMWpositivity5::THDMWpositivity5(const StandardModel& SM_i)
82: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
83{}
84
85double THDMWpositivity5::computeThValue()
86{
87 double lambda1=myTHDMW.getTHDMW_lambda1();
88 double mu1=myTHDMW.getTHDMW_mu1();
89 double mu2=myTHDMW.getTHDMW_mu2();
90 double mu3=myTHDMW.getTHDMW_mu3();
91 double mu4=myTHDMW.getTHDMW_mu4();
92 double mu5=myTHDMW.getTHDMW_mu5();
93 double mu6=myTHDMW.getTHDMW_mu6();
94 double muAtimes2=mu1+mu2+2.0*(mu3+mu4+mu5)+mu6;
95 double nu1=myTHDMW.getTHDMW_nu1();
96 if (lambda1>0 && muAtimes2>0)
97 {
98 return sqrt(lambda1*muAtimes2)+nu1;
99 }
100 else
101 {
102 return -1.;
103 }
104}
105
106
107
108THDMWpositivity6::THDMWpositivity6(const StandardModel& SM_i)
109: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
110{}
111
112double THDMWpositivity6::computeThValue()
113{
114 double lambda1=myTHDMW.getTHDMW_lambda1();
115 double mu1=myTHDMW.getTHDMW_mu1();
116 double mu2=myTHDMW.getTHDMW_mu2();
117 double mu3=myTHDMW.getTHDMW_mu3();
118 double mu4=myTHDMW.getTHDMW_mu4();
119 double mu5=myTHDMW.getTHDMW_mu5();
120 double mu6=myTHDMW.getTHDMW_mu6();
121 double muAtimes2=mu1+mu2+2.0*(mu3+mu4+mu5)+mu6;
122 double nu1=myTHDMW.getTHDMW_nu1();
123 double nu2=myTHDMW.getTHDMW_nu2();
124 double nu3=myTHDMW.getTHDMW_nu3();
125 if (lambda1>0 && muAtimes2>0)
126 {
127 return sqrt(lambda1*muAtimes2)+nu1+nu2-2.0*fabs(nu3);
128 }
129 else
130 {
131 return -1.;
132 }
133}
134
135
136
137THDMWpositivity7::THDMWpositivity7(const StandardModel& SM_i)
138: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
139{}
140
141double THDMWpositivity7::computeThValue()
142{
143 double lambda1=myTHDMW.getTHDMW_lambda1();
144 double mu1=myTHDMW.getTHDMW_mu1();
145 double mu2=myTHDMW.getTHDMW_mu2();
146 double mu3=myTHDMW.getTHDMW_mu3();
147 double mu4=myTHDMW.getTHDMW_mu4();
148 double mu5=myTHDMW.getTHDMW_mu5();
149 double mu6=myTHDMW.getTHDMW_mu6();
150 double muA=0.5*(mu1+mu2+mu6)+mu3+mu4+mu5;
151 double nu1=myTHDMW.getTHDMW_nu1();
152 double nu2=myTHDMW.getTHDMW_nu2();
153 double nu3=myTHDMW.getTHDMW_nu3();
154 double nu4=myTHDMW.getTHDMW_nu4();
155 double nu5=myTHDMW.getTHDMW_nu5();
156 if (lambda1>0 && muA>0)
157 {
158 return lambda1 + 0.5*muA + nu1 + nu2 + 2.0*nu3 - fabs(nu4+nu5)/sqrt(3.0);
159 }
160 else
161 {
162 return -1.;
163 }
164}
165
166
167
168THDMWpositivity8::THDMWpositivity8(const StandardModel& SM_i)
169: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
170{}
171
172double THDMWpositivity8::computeThValue()
173{
174 double lambda1=myTHDMW.getTHDMW_lambda1();
175 double lambda2=myTHDMW.getTHDMW_lambda2();
176 double lambda3=myTHDMW.getTHDMW_lambda3();
177
178 if (lambda1>0 && lambda2>0)
179 {
180 return lambda3+sqrt(lambda1*lambda2);
181 }
182 else
183 {
184 return -1.;
185 }
186}
187
188
189
190THDMWpositivity9::THDMWpositivity9(const StandardModel& SM_i)
191: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
192{}
193
194double THDMWpositivity9::computeThValue()
195{
196 double lambda1=myTHDMW.getTHDMW_lambda1();
197 double lambda2=myTHDMW.getTHDMW_lambda2();
198 double lambda3=myTHDMW.getTHDMW_lambda3();
199 double lambda4=myTHDMW.getTHDMW_lambda4();
200 double lambda5=myTHDMW.getTHDMW_lambda5();
201
202 if (lambda1>0 && lambda2>0)
203 {
204 return lambda3+lambda4-fabs(lambda5)+sqrt(lambda1*lambda2);
205 }
206 else
207 {
208 return -1.;
209 }
210}
211
212
213
214THDMWpositivity10::THDMWpositivity10(const StandardModel& SM_i)
215: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
216{}
217
218double THDMWpositivity10::computeThValue()
219{
220 double lambda2=myTHDMW.getTHDMW_lambda2();
221 double mu1=myTHDMW.getTHDMW_mu1();
222 double mu2=myTHDMW.getTHDMW_mu2();
223 double mu3=myTHDMW.getTHDMW_mu3();
224 double mu4=myTHDMW.getTHDMW_mu4();
225 double mu5=myTHDMW.getTHDMW_mu5();
226 double mu6=myTHDMW.getTHDMW_mu6();
227 double muAtimes2=mu1+mu2+2.0*(mu3+mu4+mu5)+mu6;
228 double omega1=myTHDMW.getTHDMW_omega1();
229 if (lambda2>0 && muAtimes2>0)
230 {
231 return sqrt(lambda2*muAtimes2)+omega1;
232 }
233 else
234 {
235 return -1.;
236 }
237}
238
239
240
241THDMWpositivity11::THDMWpositivity11(const StandardModel& SM_i)
242: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
243{}
244
245double THDMWpositivity11::computeThValue()
246{
247 double lambda2=myTHDMW.getTHDMW_lambda2();
248 double mu1=myTHDMW.getTHDMW_mu1();
249 double mu2=myTHDMW.getTHDMW_mu2();
250 double mu3=myTHDMW.getTHDMW_mu3();
251 double mu4=myTHDMW.getTHDMW_mu4();
252 double mu5=myTHDMW.getTHDMW_mu5();
253 double mu6=myTHDMW.getTHDMW_mu6();
254 double muAtimes2=mu1+mu2+2.0*(mu3+mu4+mu5)+mu6;
255 double omega1=myTHDMW.getTHDMW_omega1();
256 double omega2=myTHDMW.getTHDMW_omega2();
257 double omega3=myTHDMW.getTHDMW_omega3();
258 if (lambda2>0 && muAtimes2>0)
259 {
260 return sqrt(lambda2*muAtimes2)+omega1+omega2-2.0*fabs(omega3);
261 }
262 else
263 {
264 return -1.;
265 }
266}
267
268
269
270THDMWpositivity12::THDMWpositivity12(const StandardModel& SM_i)
271: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
272{}
273
274double THDMWpositivity12::computeThValue()
275{
276 double lambda2=myTHDMW.getTHDMW_lambda2();
277 double mu1=myTHDMW.getTHDMW_mu1();
278 double mu2=myTHDMW.getTHDMW_mu2();
279 double mu3=myTHDMW.getTHDMW_mu3();
280 double mu4=myTHDMW.getTHDMW_mu4();
281 double mu5=myTHDMW.getTHDMW_mu5();
282 double mu6=myTHDMW.getTHDMW_mu6();
283 double muA=0.5*(mu1+mu2+mu6)+mu3+mu4+mu5;
284 double omega1=myTHDMW.getTHDMW_omega1();
285 double omega2=myTHDMW.getTHDMW_omega2();
286 double omega3=myTHDMW.getTHDMW_omega3();
287 double omega4=myTHDMW.getTHDMW_omega4();
288 double omega5=myTHDMW.getTHDMW_omega5();
289 if (lambda2>0 && muA>0)
290 {
291 return lambda2 + 0.5*muA + omega1 + omega2 + 2.0*omega3 - fabs(omega4+omega5)/sqrt(3.0);
292 }
293 else
294 {
295 return -1.;
296 }
297}
298
299
300
301THDMWpositiveMassSquares::THDMWpositiveMassSquares(const StandardModel& SM_i)
302: ThObservable(SM_i),myTHDMW(static_cast<const THDMW&> (SM_i))
303{}
304
305double THDMWpositiveMassSquares::computeThValue()
306{
307 double a = myTHDMW.getMyTHDMWCache()->setOtherParameters();
308 return a;
309}
StandardModel
A model class for the Standard Model.
Definition: StandardModel.h:509
THDMW
A base class for symmetric Two-Higgs-Doublet-Manohar-Wise models.
Definition: THDMW.h:233
THDMW::getTHDMW_nu3
const double getTHDMW_nu3() const
A getter for .
Definition: THDMW.h:558
THDMW::getTHDMW_omega5
const double getTHDMW_omega5() const
A getter for .
Definition: THDMW.h:646
THDMW::getTHDMW_lambda2
const double getTHDMW_lambda2() const
A getter for .
Definition: THDMW.h:417
THDMW::getTHDMW_nu4
const double getTHDMW_nu4() const
A getter for .
Definition: THDMW.h:570
THDMW::getTHDMW_nu2
const double getTHDMW_nu2() const
A getter for .
Definition: THDMW.h:550
THDMW::getMyTHDMWCache
THDMWcache * getMyTHDMWCache() const
A method get the THDMWCache.
Definition: THDMW.h:318
THDMW::getTHDMW_mu2
const double getTHDMW_mu2() const
A getter for .
Definition: THDMW.h:488
THDMW::getTHDMW_mu3
const double getTHDMW_mu3() const
A getter for .
Definition: THDMW.h:501
THDMW::getTHDMW_lambda1
const double getTHDMW_lambda1() const
A getter for .
Definition: THDMW.h:409
THDMW::getTHDMW_mu5
const double getTHDMW_mu5() const
A getter for .
Definition: THDMW.h:516
THDMW::getTHDMW_mu1
const double getTHDMW_mu1() const
A getter for .
Definition: THDMW.h:480
THDMW::getTHDMW_omega4
const double getTHDMW_omega4() const
A getter for .
Definition: THDMW.h:633
THDMW::getTHDMW_mu4
const double getTHDMW_mu4() const
A getter for .
Definition: THDMW.h:508
THDMW::getTHDMW_nu5
const double getTHDMW_nu5() const
A getter for .
Definition: THDMW.h:578
THDMW::getTHDMW_omega3
const double getTHDMW_omega3() const
A getter for .
Definition: THDMW.h:617
THDMW::getTHDMW_lambda3
const double getTHDMW_lambda3() const
A getter for .
Definition: THDMW.h:430
THDMW::getTHDMW_omega1
const double getTHDMW_omega1() const
A getter for .
Definition: THDMW.h:591
THDMW::getTHDMW_omega2
const double getTHDMW_omega2() const
A getter for .
Definition: THDMW.h:604
THDMW::getTHDMW_mu6
const double getTHDMW_mu6() const
A getter for .
Definition: THDMW.h:529
THDMW::getTHDMW_nu1
const double getTHDMW_nu1() const
A getter for .
Definition: THDMW.h:542
THDMW::getTHDMW_lambda5
const double getTHDMW_lambda5() const
A getter for .
Definition: THDMW.h:456
THDMW::getTHDMW_lambda4
const double getTHDMW_lambda4() const
A getter for .
Definition: THDMW.h:443
THDMWcache::setOtherParameters
double setOtherParameters()
Definition: THDMWcache.cpp:2737
THDMWpositiveMassSquares::THDMWpositiveMassSquares
THDMWpositiveMassSquares(const StandardModel &SM_i)
THDMWpositiveMassSquares constructor.
Definition: THDMWpositivity.cpp:301
THDMWpositivity10::THDMWpositivity10
THDMWpositivity10(const StandardModel &SM_i)
THDMWpositivity10 constructor.
Definition: THDMWpositivity.cpp:214
THDMWpositivity10::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:218
THDMWpositivity11::THDMWpositivity11
THDMWpositivity11(const StandardModel &SM_i)
THDMWpositivity11 constructor.
Definition: THDMWpositivity.cpp:241
THDMWpositivity11::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:237
THDMWpositivity12::THDMWpositivity12
THDMWpositivity12(const StandardModel &SM_i)
THDMWpositivity12 constructor.
Definition: THDMWpositivity.cpp:270
THDMWpositivity12::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:256
THDMWpositivity1::THDMWpositivity1
THDMWpositivity1(const StandardModel &SM_i)
THDMWpositivity1 constructor.
Definition: THDMWpositivity.cpp:12
THDMWpositivity1::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:44
THDMWpositivity2::THDMWpositivity2
THDMWpositivity2(const StandardModel &SM_i)
THDMWpositivity2 constructor.
Definition: THDMWpositivity.cpp:30
THDMWpositivity2::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:66
THDMWpositivity3::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:85
THDMWpositivity3::THDMWpositivity3
THDMWpositivity3(const StandardModel &SM_i)
THDMWpositivity3 constructor.
Definition: THDMWpositivity.cpp:46
THDMWpositivity4::THDMWpositivity4
THDMWpositivity4(const StandardModel &SM_i)
THDMWpositivity4 constructor.
Definition: THDMWpositivity.cpp:63
THDMWpositivity4::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:104
THDMWpositivity5::THDMWpositivity5
THDMWpositivity5(const StandardModel &SM_i)
THDMWpositivity5 constructor.
Definition: THDMWpositivity.cpp:81
THDMWpositivity5::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:123
THDMWpositivity6::THDMWpositivity6
THDMWpositivity6(const StandardModel &SM_i)
THDMWpositivity6 constructor.
Definition: THDMWpositivity.cpp:108
THDMWpositivity6::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:142
THDMWpositivity7::THDMWpositivity7
THDMWpositivity7(const StandardModel &SM_i)
THDMWpositivity7 constructor.
Definition: THDMWpositivity.cpp:137
THDMWpositivity7::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:161
THDMWpositivity8::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:180
THDMWpositivity8::THDMWpositivity8
THDMWpositivity8(const StandardModel &SM_i)
THDMWpositivity8 constructor.
Definition: THDMWpositivity.cpp:168
THDMWpositivity9::THDMWpositivity9
THDMWpositivity9(const StandardModel &SM_i)
THDMWpositivity9 constructor.
Definition: THDMWpositivity.cpp:190
THDMWpositivity9::myTHDMW
const THDMW & myTHDMW
Definition: THDMWpositivity.h:199
ThObservable
A class for a model prediction of an observable.
Definition: ThObservable.h:25
lambda1
parameter of the Higgs potential
Definition: THDMquantities.h:382
lambda2
An observable class for the quartic Higgs potential coupling .
Definition: THDMquantities.h:405
lambda3
An observable class for the quartic Higgs potential coupling .
Definition: THDMquantities.h:428
lambda4
An observable class for the quartic Higgs potential coupling .
Definition: THDMquantities.h:451
lambda5
An observable class for the quartic Higgs potential coupling .
Definition: THDMquantities.h:474