a Code for the Combination of Indirect and Direct Constraints on High Energy Physics Models Logo
THDMquantities.cpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2015 HEPfit Collaboration
3 *
4 *
5 * For the licensing terms see doc/COPYING.
6 */
7
8#include "THDMquantities.h"
9#include "THDM.h"
10
12: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
13{}
14
16{
17 return myTHDM->gettanb();
18}
19
21: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
22{}
23
25{
26 return sqrt(myTHDM->getmHl2());
27}
28
30: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
31{}
32
34{
35 return myTHDM->getmHh();
36}
37
38
40: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
41{}
42
44{
45 return myTHDM->getmA();
46}
47
48
50: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
51{}
52
54{
55 return myTHDM->getmHp();
56}
57
58
60: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
61{}
62
64{
65 return sqrt(myTHDM->getmHl2()) - myTHDM->getmA();
66}
67
68
70: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
71{}
72
74{
75 return myTHDM->getmA() - sqrt(myTHDM->getmHl2());
76}
77
78
80: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
81{}
82
84{
85 return sqrt(myTHDM->getmHl2()) - myTHDM->getmHp();
86}
87
88
90: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
91{}
92
94{
95 return myTHDM->getmHp() - sqrt(myTHDM->getmHl2());
96}
97
98
100: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
101{}
102
104{
105 return myTHDM->getmHh() - myTHDM->getmA();
106}
107
108
110: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
111{}
112
114{
115 return myTHDM->getmA() - myTHDM->getmHh();
116}
117
118
120: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
121{}
122
124{
125 return myTHDM->getmHh() - myTHDM->getmHp();
126}
127
128
130: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
131{}
132
134{
135 return myTHDM->getmHp() - myTHDM->getmHh();
136}
137
138
140: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
141{}
142
144{
145 return myTHDM->getmA() - myTHDM->getmHp();
146}
147
148
150: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
151{}
152
154{
155 return myTHDM->getmHp() - myTHDM->getmA();
156}
157
158
160: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
161{}
162
164{
165 double mHl2=myTHDM->getmHl2();
166 double tanb=myTHDM->gettanb();
167 double mHh2=myTHDM->getmHh2();
168 double sina=myTHDM->getsina();
169 double cosa=myTHDM->getcosa();
170 double m12_2=myTHDM->getm12_2();
171
172 return (-(cosa*cosa*mHh2) - mHl2*sina*sina + 2.0*m12_2*tanb + cosa*(mHl2 - mHh2)*sina*tanb)/2.;
173}
174
175
177: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
178{}
179
181{
182 double mHl2=myTHDM->getmHl2();
183 double tanb=myTHDM->gettanb();
184 double mHh2=myTHDM->getmHh2();
185 double sina=myTHDM->getsina();
186 double cosa=myTHDM->getcosa();
187 double m12_2=myTHDM->getm12_2();
188
189 return (2.0*m12_2 + sina*cosa*(mHl2 - mHh2) - (cosa*cosa*mHl2 + sina*sina*mHh2)*tanb)/(2.0*tanb);
190}
191
192
194: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
195{}
196
198{
199 double mHl2=myTHDM->getmHl2();
200 double vev=myTHDM->v();
201 double tanb=myTHDM->gettanb();
202 double cosb=myTHDM->getcosb();
203 double mHh2=myTHDM->getmHh2();
204 double sina=myTHDM->getsina();
205 double cosa=myTHDM->getcosa();
206 double m12_2=myTHDM->getm12_2();
207
208 return (mHh2*cosa*cosa+mHl2*sina*sina-m12_2*tanb)/(vev*vev*cosb*cosb);
209}
210
211
213: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
214{}
215
217{
218 double mHl2=myTHDM->getmHl2();
219 double vev=myTHDM->v();
220 double tanb=myTHDM->gettanb();
221 double sinb=myTHDM->getsinb();
222 double mHh2=myTHDM->getmHh2();
223 double sina=myTHDM->getsina();
224 double cosa=myTHDM->getcosa();
225 double m12_2=myTHDM->getm12_2();
226
227 return (mHh2*sina*sina+mHl2*cosa*cosa-m12_2/tanb)/(vev*vev*sinb*sinb);
228}
229
230
232: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
233{}
234
236{
237 double mHl2=myTHDM->getmHl2();
238 double vev=myTHDM->v();
239 double sinb=myTHDM->getsinb();
240 double cosb=myTHDM->getcosb();
241 double mHh2=myTHDM->getmHh2();
242 double mHp2=myTHDM->getmHp2();
243 double sina=myTHDM->getsina();
244 double cosa=myTHDM->getcosa();
245 double m12_2=myTHDM->getm12_2();
246
247 return ((mHh2-mHl2)*cosa*sina+2.0*mHp2*sinb*cosb-m12_2)/(vev*vev*sinb*cosb);
248}
249
250
252: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
253{}
254
256{
257 double mA2=myTHDM->getmA2();
258 double mHp2=myTHDM->getmHp2();
259 double vev=myTHDM->v();
260 double sinb=myTHDM->getsinb();
261 double cosb=myTHDM->getcosb();
262 double m12_2=myTHDM->getm12_2();
263
264 return ((mA2-2.0*mHp2)*sinb*cosb+m12_2)/(vev*vev*sinb*cosb);
265}
266
267
269: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
270{}
271
273{
274 double mA2=myTHDM->getmA2();
275 double vev=myTHDM->v();
276 double sinb=myTHDM->getsinb();
277 double cosb=myTHDM->getcosb();
278 double m12_2=myTHDM->getm12_2();
279
280 return (m12_2-mA2*sinb*cosb)/(vev*vev*sinb*cosb);
281}
282
283
285: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
286{}
287
289{
290 double mHl2=myTHDM->getmHl2();
291 double vev=myTHDM->v();
292 double sinb=myTHDM->getsinb();
293 double cosb=myTHDM->getcosb();
294 double mHh2=myTHDM->getmHh2();
295 double sina=myTHDM->getsina();
296 double cosa=myTHDM->getcosa();
297 double m12_2=myTHDM->getm12_2();
298
299 return (m12_2+(mHh2-mHl2)*cosa*sina)/(vev*vev*sinb*cosb);
300}
301
302
304: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
305{}
306
308{
309 double mHl2=myTHDM->getmHl2();
310 double vev=myTHDM->v();
311 double sinb=myTHDM->getsinb();
312 double cosb=myTHDM->getcosb();
313 double sina=myTHDM->getsina();
314 double cosa=myTHDM->getcosa();
315 double bma=myTHDM->getbma();
316 double m12_2=myTHDM->getm12_2();
317
318 return 3.0*( -mHl2*(cosa*cosa*cosa/sinb - sina*sina*sina/cosb)
319 +cos(bma)*m12_2*(cosa*cosa/(sinb*sinb) - sina*sina/(cosb*cosb)) )/vev;
320}
321
322
324: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
325{}
326
328{
329 double mHl2=myTHDM->getmHl2();
330 double mHh2=myTHDM->getmHh2();
331 double vev=myTHDM->v();
332 double sinb=myTHDM->getsinb();
333 double cosb=myTHDM->getcosb();
334 double sina=myTHDM->getsina();
335 double cosa=myTHDM->getcosa();
336 double bma=myTHDM->getbma();
337 double m12_2=myTHDM->getm12_2();
338
339 return -cos(bma)*(sinb*cosb*m12_2
340 +sina*cosa*(sinb*cosb*(2.0*mHl2 + mHh2) - 3.0*m12_2))/(sinb*sinb*cosb*cosb*vev);
341}
342
343
345: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
346{}
347
349{
350 double mHl2=myTHDM->getmHl2();
351 double mHh2=myTHDM->getmHh2();
352 double vev=myTHDM->v();
353 double sinb=myTHDM->getsinb();
354 double cosb=myTHDM->getcosb();
355 double sina=myTHDM->getsina();
356 double cosa=myTHDM->getcosa();
357 double bma=myTHDM->getbma();
358 double m12_2=myTHDM->getm12_2();
359
360 return -sin(bma)*(sinb*cosb*m12_2
361 -sina*cosa*(sinb*cosb*(mHl2 + 2.0*mHh2) - 3.0*m12_2))/(sinb*sinb*cosb*cosb*vev);
362}
363
364
366: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
367{}
368
370{
371 double mHh2=myTHDM->getmHh2();
372 double vev=myTHDM->v();
373 double sinb=myTHDM->getsinb();
374 double cosb=myTHDM->getcosb();
375 double sina=myTHDM->getsina();
376 double cosa=myTHDM->getcosa();
377 double bma=myTHDM->getbma();
378 double m12_2=myTHDM->getm12_2();
379
380 return -3.0*( (cosa*cosa*cosa*sinb + sina*sina*sina*cosb)*sinb*cosb*mHh2
381 -(cosa*cosa*sinb*sinb - cosb*cosb*sina*sina)*m12_2*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
382}
383
384
386: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
387{}
388
390{
391 double mHl2=myTHDM->getmHl2();
392 double mA2=myTHDM->getmA2();
393 double vev=myTHDM->v();
394 double sinb=myTHDM->getsinb();
395 double cosb=myTHDM->getcosb();
396 double sina=myTHDM->getsina();
397 double cosa=myTHDM->getcosa();
398 double bma=myTHDM->getbma();
399 double m12_2=myTHDM->getm12_2();
400
401 return ( (cosa*cosb-sina*sinb)*(m12_2-sinb*cosb*mHl2)
402 +(mHl2-2.0*mA2)*sinb*sinb*cosb*cosb*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
403}
404
405
407: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
408{}
409
411{
412 double mHh2=myTHDM->getmHh2();
413 double mA2=myTHDM->getmA2();
414 double vev=myTHDM->v();
415 double sinb=myTHDM->getsinb();
416 double cosb=myTHDM->getcosb();
417 double sina=myTHDM->getsina();
418 double cosa=myTHDM->getcosa();
419 double bma=myTHDM->getbma();
420 double m12_2=myTHDM->getm12_2();
421
422 return ( (cosa*sinb+cosb*sina)*(m12_2-sinb*cosb*mHh2)
423 +(mHh2-2.0*mA2)*sinb*sinb*cosb*cosb*cos(bma) )/(sinb*sinb*cosb*cosb*vev);
424}
425
426
428: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
429{}
430
432{
433 double mHl2=myTHDM->getmHl2();
434 double mHp2=myTHDM->getmHp2();
435 double vev=myTHDM->v();
436 double sinb=myTHDM->getsinb();
437 double cosb=myTHDM->getcosb();
438 double sina=myTHDM->getsina();
439 double cosa=myTHDM->getcosa();
440 double bma=myTHDM->getbma();
441 double m12_2=myTHDM->getm12_2();
442
443 return ( (cosa*cosb-sina*sinb)*(m12_2-sinb*cosb*mHl2)
444 +(mHl2-2.0*mHp2)*sinb*sinb*cosb*cosb*sin(bma) )/(sinb*sinb*cosb*cosb*vev);
445}
446
447
449: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
450{}
451
453{
454 double mHh2=myTHDM->getmHh2();
455 double mHp2=myTHDM->getmHp2();
456 double vev=myTHDM->v();
457 double sinb=myTHDM->getsinb();
458 double cosb=myTHDM->getcosb();
459 double sina=myTHDM->getsina();
460 double cosa=myTHDM->getcosa();
461 double bma=myTHDM->getbma();
462 double m12_2=myTHDM->getm12_2();
463
464 return ( (cosa*sinb+cosb*sina)*(m12_2-sinb*cosb*mHh2)
465 +(mHh2-2.0*mHp2)*sinb*sinb*cosb*cosb*cos(bma) )/(sinb*sinb*cosb*cosb*vev);
466}
467
468
470: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
471{}
472
474{
475 double mHl2=myTHDM->getmHl2();
476 double mHh2=myTHDM->getmHh2();
477 double bma=myTHDM->getbma();
478
479 return 0.25*(-mHl2-mHh2+(mHl2-mHh2)*cos(2.0*bma));
480}
481
482
484: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
485{}
486
488{
489 double mHl2=myTHDM->getmHl2();
490 double sinb=myTHDM->getsinb();
491 double cosb=myTHDM->getcosb();
492 double sin2b=2.0*sinb*cosb;
493 double mHh2=myTHDM->getmHh2();
494 double bma=myTHDM->getbma();
495 double sina=myTHDM->getsina();
496 double cosa=myTHDM->getcosa();
497 double sin2a=2.0*sina*cosa;
498 double m12_2=myTHDM->getm12_2();
499
500
501 return 0.125*((mHl2-mHh2)*(3.0*sin2a-sin2b*cos(2.0*bma)
502 -sin(2.0*bma)*(cosb*cosb-sinb*sinb))
503 +2.0*(8.0*m12_2-(mHl2+mHh2)*sin2b))/sin2b;
504}
505
506
508: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
509{}
510
512{
513 double mHl2=myTHDM->getmHl2();
514 double mHh2=myTHDM->getmHh2();
515 double bma=myTHDM->getbma();
516
517 return 0.25*(mHh2-mHl2)*sin(2.0*bma);
518}
519
520
522: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
523{}
524
526{
527 double mHl2=myTHDM->getmHl2();
528 double vev=myTHDM->v();
529 double mHh2=myTHDM->getmHh2();
530 double bma=myTHDM->getbma();
531
532 return (mHl2+mHh2+(mHh2-mHl2)*cos(2.0*bma))/(2.0*vev*vev);
533}
534
535
537: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
538{}
539
541{
542 double mHl2=myTHDM->getmHl2();
543 double vev=myTHDM->v();
544 double sinb=myTHDM->getsinb();
545 double cosb=myTHDM->getcosb();
546 double sin2b=2.0*sinb*cosb;
547 double cos2b=cosb*cosb-sinb*sinb;
548 double mHh2=myTHDM->getmHh2();
549 double bma=myTHDM->getbma();
550 double m12_2=myTHDM->getm12_2();
551
552 return (-128.0*m12_2*cos2b*cos2b
553 +4.0*(mHl2+mHh2)*(5.0+3.0*(cos2b*cos2b-sin2b*sin2b))*sin2b
554 +2.0*(mHl2-mHh2)*((6.0+10.0*(cos2b*cos2b-sin2b*sin2b))*cos(2.0*bma)
555 +16.0*sin2b*cos2b*sin(2.0*bma))*sin2b)
556 /(16.0*sin2b*sin2b*sin2b*vev*vev);
557}
558
559
561: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
562{}
563
565{
566 double mHl2=myTHDM->getmHl2();
567 double vev=myTHDM->v();
568 double sinb=myTHDM->getsinb();
569 double cosb=myTHDM->getcosb();
570 double sin2b=2.0*sinb*cosb;
571 double cos2b=cosb*cosb-sinb*sinb;
572 double mHh2=myTHDM->getmHh2();
573 double mHp2=myTHDM->getmHp2();
574 double bma=myTHDM->getbma();
575 double sina=myTHDM->getsina();
576 double cosa=myTHDM->getcosa();
577 double sin2a=2.0*sina*cosa;
578 double m12_2=myTHDM->getm12_2();
579
580 return ((mHh2-mHl2)*(3.0*sin2a-sin(2.0*bma)*cos2b-cos(2.0*bma)*sin2b)
581 +2.0*(-8.0*m12_2+(mHh2+mHl2+4.0*mHp2)*sin2b))/(4.0*sin2b*vev*vev);
582}
583
584
586: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
587{}
588
590{
591 double mHl2=myTHDM->getmHl2();
592 double vev=myTHDM->v();
593 double mHh2=myTHDM->getmHh2();
594 double mA2=myTHDM->getmA2();
595 double mHp2=myTHDM->getmHp2();
596 double bma=myTHDM->getbma();
597
598 return (2.0*mA2+mHl2+mHh2-4.0*mHp2+(mHl2-mHh2)*cos(2.0*bma))/(2.0*vev*vev);
599}
600
601
603: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
604{}
605
607{
608 double mHl2=myTHDM->getmHl2();
609 double vev=myTHDM->v();
610 double mHh2=myTHDM->getmHh2();
611 double mA2=myTHDM->getmA2();
612 double bma=myTHDM->getbma();
613
614 return (-2.0*mA2+mHl2+mHh2+(mHl2-mHh2)*cos(2.0*bma))/(2.0*vev*vev);
615}
616
617
619: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
620{}
621
623{
624 double mHl2=myTHDM->getmHl2();
625 double vev=myTHDM->v();
626 double mHh2=myTHDM->getmHh2();
627 double bma=myTHDM->getbma();
628
629 return ((mHl2-mHh2)*sin(2.0*bma))/(2.0*vev*vev);
630}
631
632
634: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
635{}
636
638{
639 double mHl2=myTHDM->getmHl2();
640 double vev=myTHDM->v();
641 double sinb=myTHDM->getsinb();
642 double cosb=myTHDM->getcosb();
643 double sin2b=2.0*sinb*cosb;
644 double cos2b=cosb*cosb-sinb*sinb;
645 double mHh2=myTHDM->getmHh2();
646 double bma=myTHDM->getbma();
647 double sina=myTHDM->getsina();
648 double cosa=myTHDM->getcosa();
649 double m12_2=myTHDM->getm12_2();
650
651 return (4.0*cos2b*(-4.0*m12_2+(mHl2+mHh2)*sin2b)
652 +(mHl2-mHh2)*(3.0*(cosa*cosa-sina*sina)
653 +cos(2.0*bma)*cos2b-sin(2.0*bma)*sin2b)*sin2b)
654 /(4.0*sin2b*sin2b*vev*vev);
655}
656
657
659: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
660{}
661
663{
664 double mHl2=myTHDM->getmHl2();
665 double sinb=myTHDM->getsinb();
666 double cosb=myTHDM->getcosb();
667 double sin2b=2.0*sinb*cosb;
668 double mHh2=myTHDM->getmHh2();
669 double sina=myTHDM->getsina();
670 double cosa=myTHDM->getcosa();
671 double sin2a=2.0*sina*cosa;
672 double m12_2=myTHDM->getm12_2();
673
674 return -((mHh2+mHl2)*sin2b+(mHh2-mHl2)*sin2a-4.0*m12_2)/(4.0*sin2b);
675}
676
677
679: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
680{}
681
683{
684 double m12_2=myTHDM->getm12_2();
685
686 return -m12_2;
687}
688
689
691: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
692{}
693
695{
696 double mHl2=myTHDM->getmHl2();
697 double sinb=myTHDM->getsinb();
698 double cosb=myTHDM->getcosb();
699 double sin2b=2.0*sinb*cosb;
700 double cos2b=cosb*cosb-sinb*sinb;
701 double mHh2=myTHDM->getmHh2();
702 double bma=myTHDM->getbma();
703 double m12_2=myTHDM->getm12_2();
704
705 return ((mHl2-mHh2)*sin(2.0*bma)-4.0*m12_2*cos2b)/(4.0*sin2b);
706}
707
708
710: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
711{}
712
714{
715 double mHl2=myTHDM->getmHl2();
716 double vev=myTHDM->v();
717 double sinb=myTHDM->getsinb();
718 double cosb=myTHDM->getcosb();
719 double sin2b=2.0*sinb*cosb;
720 double cos2b=cosb*cosb-sinb*sinb;
721 double mHh2=myTHDM->getmHh2();
722 double mHp2=myTHDM->getmHp2();
723 double sina=myTHDM->getsina();
724 double cosa=myTHDM->getcosa();
725 double sin2a=2.0*sina*cosa;
726 double cos2a=cosa*cosa-sina*sina;
727 double m12_2=myTHDM->getm12_2();
728
729 return ((mHh2+mHl2-(mHh2-mHl2)*(cos2a*cos2b-sin2a*sin2b)
730 -4.0*m12_2/sin2b)/(sin2b*sin2b)
731 +2.0*mHp2)/(4.0*vev*vev);
732}
733
734
736: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
737{}
738
740{
741 double mHl2=myTHDM->getmHl2();
742 double vev=myTHDM->v();
743 double tanb=myTHDM->gettanb();
744 double sinb=myTHDM->getsinb();
745 double cosb=myTHDM->getcosb();
746 double sin2b=2.0*sinb*cosb;
747 double cos2b=cosb*cosb-sinb*sinb;
748 double mHh2=myTHDM->getmHh2();
749 double sina=myTHDM->getsina();
750 double cosa=myTHDM->getcosa();
751 double cos2a=cosa*cosa-sina*sina;
752 double m12_2=myTHDM->getm12_2();
753
754 return ((mHh2+mHl2)*(0.5*tanb-0.5/tanb)
755 +(mHh2-mHl2)*cos2a/sin2b
756 +4.0*m12_2*cos2b/(sin2b*sin2b))/(4.0*vev*vev*sin2b);
757}
758
759
761: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
762{}
763
765{
766 double vev=myTHDM->v();
767 double sinb=myTHDM->getsinb();
768 double cosb=myTHDM->getcosb();
769 double sin2b=2.0*sinb*cosb;
770 double mHp2=myTHDM->getmHp2();
771 double m12_2=myTHDM->getm12_2();
772
773 return (m12_2/sin2b-0.5*mHp2)/(vev*vev);
774}
775
776
778: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
779{}
780
782{
783 double vev=myTHDM->v();
784 double mA2=myTHDM->getmA2();
785 double mHp2=myTHDM->getmHp2();
786
787 return (mA2-mHp2)/(2.0*vev*vev);
788}
789
790
792: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
793{}
794
796{
797 double mHl2=myTHDM->getmHl2();
798 double vev=myTHDM->v();
799 double sinb=myTHDM->getsinb();
800 double cosb=myTHDM->getcosb();
801 double sin2b=2.0*sinb*cosb;
802 double cos2b=cosb*cosb-sinb*sinb;
803 double mHh2=myTHDM->getmHh2();
804 double mHp2=myTHDM->getmHp2();
805 double bma=myTHDM->getbma();
806 double m12_2=myTHDM->getm12_2();
807
808 return ((mHh2+mHl2-(mHh2-mHl2)*cos(2.0*bma)
809 -4.0*m12_2*cos2b*cos2b/sin2b)/(sin2b*sin2b)
810 -2.0*mHp2)/(4.0*vev*vev);
811}
812
813
814
816: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
817{}
818
820{
821 double mHl2=myTHDM->getmHl2();
822 double vev=myTHDM->v();
823 double tanb=myTHDM->gettanb();
824 double sinb=myTHDM->getsinb();
825 double cosb=myTHDM->getcosb();
826 double sin2b=2.0*sinb*cosb;
827 double cos2b=cosb*cosb-sinb*sinb;
828 double mHh2=myTHDM->getmHh2();
829 double sina=myTHDM->getsina();
830 double cosa=myTHDM->getcosa();
831 double sin2a=2.0*sina*cosa;
832 double cos2a=cosa*cosa-sina*sina;
833 double m12_2=myTHDM->getm12_2();
834
835 return (mHh2+mHl2+4.0*m12_2*cos2b/sin2b
836 +(mHh2-mHl2)*(cos2a-sin2a/tanb))/(4.0*vev*vev*cosb*cosb);
837}
838
839
841: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
842{}
843
845{
846 double mHl2=myTHDM->getmHl2();
847 double vev=myTHDM->v();
848 double tanb=myTHDM->gettanb();
849 double sinb=myTHDM->getsinb();
850 double cosb=myTHDM->getcosb();
851 double sin2b=2.0*sinb*cosb;
852 double cos2b=cosb*cosb-sinb*sinb;
853 double mHh2=myTHDM->getmHh2();
854 double sina=myTHDM->getsina();
855 double cosa=myTHDM->getcosa();
856 double sin2a=2.0*sina*cosa;
857 double cos2a=cosa*cosa-sina*sina;
858 double m12_2=myTHDM->getm12_2();
859
860 return (mHh2+mHl2-4.0*m12_2*cos2b/sin2b
861 -(mHh2-mHl2)*(cos2a+sin2a*tanb))/(4.0*vev*vev*sinb*sinb);
862}
863
864
866: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
867{}
868
870{
871 double mHl2=myTHDM->getmHl2();
872 double vev=myTHDM->v();
873 double sinb=myTHDM->getsinb();
874 double cosb=myTHDM->getcosb();
875 double sin2b=2.0*sinb*cosb;
876 double mHh2=myTHDM->getmHh2();
877 double sina=myTHDM->getsina();
878 double cosa=myTHDM->getcosa();
879 double sin2a=2.0*sina*cosa;
880 double m12_2=myTHDM->getm12_2();
881
882 return ((mHh2-mHl2)*sin2a-2.0*m12_2)/(2.0*vev*vev*sin2b);
883}
884
885
887: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
888{}
889
891{
892 double vev=myTHDM->v();
893 double mHp2=myTHDM->getmHp2();
894
895 return (2.0*mHp2)/(vev*vev);
896}
897
898
900: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
901{}
902
904{
905 double vev=myTHDM->v();
906 double sinb=myTHDM->getsinb();
907 double cosb=myTHDM->getcosb();
908 double sin2b=2.0*sinb*cosb;
909 double m12_2=myTHDM->getm12_2();
910
911 return 4.0*m12_2/(vev*vev*sin2b);
912}
913
914
916: ThObservable(SM_i), myTHDM(static_cast<const THDM*> (&SM_i))
917{}
918
920{
921 double vev=myTHDM->v();
922 double mA2=myTHDM->getmA2();
923
924 return 2.0*mA2/(vev*vev);
925}
double computeThValue()
The quartic coupling .
const THDM * myTHDM
E11_THDM(const StandardModel &SM_i)
E11_THDM constructor.
double computeThValue()
The quartic coupling .
E22_THDM(const StandardModel &SM_i)
E22_THDM constructor.
const THDM * myTHDM
E33_THDM(const StandardModel &SM_i)
E33_THDM constructor.
double computeThValue()
The quartic coupling .
const THDM * myTHDM
double computeThValue()
The quartic coupling .
const THDM * myTHDM
HHlambda1(const StandardModel &SM_i)
HHlambda1 constructor.
double computeThValue()
The quartic coupling .
HHlambda2(const StandardModel &SM_i)
HHlambda2 constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
HHlambda3(const StandardModel &SM_i)
HHlambda3 constructor.
const THDM * myTHDM
HHlambda4(const StandardModel &SM_i)
HHlambda4 constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
const THDM * myTHDM
double computeThValue()
The quartic coupling .
HHlambda5(const StandardModel &SM_i)
HHlambda5 constructor.
const THDM * myTHDM
HHlambda6(const StandardModel &SM_i)
HHlambda6 constructor.
double computeThValue()
The quartic coupling .
A model class for the Standard Model.
const double v() const
The Higgs vacuum expectation value.
A base class for symmetric Two-Higgs-Doublet models.
Definition: THDM.h:120
double getmHp() const
A method get the charged Higgs mass.
Definition: THDM.h:471
double getmHh() const
A method get the mass of the "non-125 GeV" neutral scalar Higgs.
Definition: THDM.h:394
double gettanb() const
A method get .
Definition: THDM.h:283
double getbma() const
A method get .
Definition: THDM.h:307
double getmHh2() const
A method get the squared mass of the "non-125 GeV" neutral scalar Higgs.
Definition: THDM.h:365
double getsinb() const
A method get .
Definition: THDM.h:291
double getmHl2() const
A method get the squared mass of the lighter neutral scalar Higgs.
Definition: THDM.h:339
double getmHp2() const
A method get the squared charged Higgs mass.
Definition: THDM.h:457
double getmA() const
A method get the mass of the pseudoscalar Higgs A.
Definition: THDM.h:437
double getm12_2() const
A method get the parameter of the Higgs potential .
Definition: THDM.h:491
double getcosb() const
A method get .
Definition: THDM.h:299
double getcosa() const
A method get .
Definition: THDM.h:323
double getsina() const
A method get .
Definition: THDM.h:331
double getmA2() const
A method get the squared mass of the pseudoscalar Higgs A.
Definition: THDM.h:423
A class for a model prediction of an observable.
Definition: ThObservable.h:25
Y1_THDM(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
double computeThValue()
The quadratic coupling .
double computeThValue()
The quadratic coupling .
const THDM * myTHDM
Y2_THDM(const StandardModel &SM_i)
Constructor.
Y3_THDM(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quadratic coupling .
const THDM * myTHDM
double computeThValue()
The quartic coupling .
const THDM * myTHDM
Z1_THDM(const StandardModel &SM_i)
Z1_THDM constructor.
double computeThValue()
The quartic coupling .
Z2_THDM(const StandardModel &SM_i)
Z2_THDM constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
Z3_THDM(const StandardModel &SM_i)
Z3_THDM constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
Z4_THDM(const StandardModel &SM_i)
Z4_THDM constructor.
const THDM * myTHDM
const THDM * myTHDM
Z5_THDM(const StandardModel &SM_i)
Z5_THDM constructor.
double computeThValue()
The quartic coupling .
const THDM * myTHDM
Z6_THDM(const StandardModel &SM_i)
Z6_THDM constructor.
double computeThValue()
The quartic coupling .
Z7_THDM(const StandardModel &SM_i)
Z7_THDM constructor.
double computeThValue()
The quartic coupling .
const THDM * myTHDM
eta00_THDM(const StandardModel &SM_i)
eta00_THDM constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
const THDM * myTHDM
double computeThValue()
The quartic coupling .
eta3_THDM(const StandardModel &SM_i)
eta3_THDM constructor.
g_HhAA(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
const THDM * myTHDM
const THDM * myTHDM
double computeThValue()
The triple Higgs coupling .
g_HhHhHh(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
g_HhHpHm(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
const THDM * myTHDM
g_hAA(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
const THDM * myTHDM
double computeThValue()
The triple Higgs coupling .
g_hHhHh(const StandardModel &SM_i)
Constructor.
g_hHpHm(const StandardModel &SM_i)
Constructor.
double computeThValue()
The triple Higgs coupling .
const THDM * myTHDM
g_hhHh(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
double computeThValue()
The triple Higgs coupling .
g_hhh(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
double computeThValue()
The triple Higgs coupling .
lambda1(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
const THDM * myTHDM
double computeThValue()
The quartic coupling .
lambda2(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
lambda345(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling combination .
const THDM * myTHDM
lambda3(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
double computeThValue()
The quartic coupling .
const THDM * myTHDM
lambda4(const StandardModel &SM_i)
Constructor.
lambda5(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
const THDM * myTHDM
m11_2(const StandardModel &SM_i)
Constructor.
const THDM * myTHDM
double computeThValue()
The quartic coupling .
m22_2(const StandardModel &SM_i)
Constructor.
double computeThValue()
The quartic coupling .
const THDM * myTHDM
const THDM * myTHDM
mass_mA(const StandardModel &SM_i)
mass_mA constructor.
double computeThValue()
double computeThValue()
mass_mHh(const StandardModel &SM_i)
mass_mHh constructor.
const THDM * myTHDM
mass_mHl(const StandardModel &SM_i)
mass_mHl constructor.
double computeThValue()
const THDM * myTHDM
double computeThValue()
mass_mHp(const StandardModel &SM_i)
mass_mHp constructor.
const THDM * myTHDM
massdifference_mAmmHh(const StandardModel &SM_i)
massdifference_mAmmHh constructor.
massdifference_mAmmHl(const StandardModel &SM_i)
massdifference_mAmmHl constructor.
massdifference_mAmmHp(const StandardModel &SM_i)
massdifference_mAmmHp constructor.
massdifference_mHhmmA(const StandardModel &SM_i)
massdifference_mHhmmA constructor.
massdifference_mHhmmHp(const StandardModel &SM_i)
massdifference_mHhmmHp constructor.
massdifference_mHlmmA(const StandardModel &SM_i)
massdifference_mHlmmA constructor.
massdifference_mHlmmHp(const StandardModel &SM_i)
massdifference_mHlmmHp constructor.
massdifference_mHpmmA(const StandardModel &SM_i)
massdifference_mHpmmA constructor.
massdifference_mHpmmHh(const StandardModel &SM_i)
massdifference_mHpmmHh constructor.
massdifference_mHpmmHl(const StandardModel &SM_i)
massdifference_mHpmmHl constructor.
tanbeta(const StandardModel &SM_i)
tanbeta constructor.
double computeThValue()
const THDM * myTHDM
double computeThValue()
The quadratic coupling .
const THDM * myTHDM
xi0_THDM(const StandardModel &SM_i)
xi0_THDM constructor.
double computeThValue()
The quadratic coupling .
xi1_THDM(const StandardModel &SM_i)
xi1_THDM constructor.
const THDM * myTHDM
double computeThValue()
The quadratic coupling .
xi3_THDM(const StandardModel &SM_i)
xi3_THDM constructor.
const THDM * myTHDM