Example
An example script that demonstrates the usage of thdmTools can be found in the example folder.
This script initializes a 2HDM parameter point and applies a set of theoretical and experimental constraints. The results are printed to the terminal. See below for the source code.
import numpy as np
from thdmtools.parampoint import Parampoint
TB = 2.0
dc = {
'type': 2,
'tb': TB,
'alpha': np.arctan(TB) + 1.01 * np.pi / 2.,
'mHh': 500.0,
'mHl': 125.09,
'mA': 600.0,
'mHp': 700.0,
'm12sq': 500.0**2} # 156000.0}
pt = Parampoint(dc)
print(pt.check_vacuum_stability())
print(pt.check_perturbative_unitarity())
print(pt.check_flavour_constraints())
print(pt.check_ewpo_constraints())
print(pt.check_ewpo_fit())
print(pt.check_collider_constraints())
print(pt.reshb)
print(pt.chisqhs)
print(pt.delta_chisqhs)