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)