The basic properties of the Higgs boson in the BGL-like 3HDM under consideration, such as its physical couplings, decay rates and branching ratios, have been computed numerically in SPheno-4.0.4 for each input parameter-space point in a dedicated numerical scan. A set of input
points has been prepared in a separate special routine using an inversion procedure in the scalar sector in order to reproduce the exact Higgs alignment limit, with fixed VEV \(v\simeq 246\)  GeV and Higgs mass \(m_h\simeq 125\)  GeV values, as detailed above. Then, a candidate point has been chosen to satisfy the electroweak precision bounds, the boundedness-from-below and unitarity constraints discussed in the previous subsection. On top of that, in order to take into account the latest data from the LHC on the properties of the observed 125 GeV state we used HiggsSignals-2.2.3beta [1]}, whereas limits from the LEP, Tevatron and LHC direct searches for new CP-even and CP-odd scalars were taken into account through the use of HiggsBounds-5.3.2beta [2]}, [3]}, [4]}.
The latter two packages are linked to SPheno-4.0.4 such that they get all the necessary inputs containing masses and decay widths of the SM and BSM Higgs states for each parameter-space point. A point that passes all these constraints is considered to be valid, at least, from the point of view of the scalar sector constraints.
