Here the underlying code will be in C++, with SLHA writing and reading performed via python using PySLHA (See https://pypi.org/project/pyslha/ & http://insectnation.org/projects/pyslha.html)
Create two scripts to handle the model creation, which consists of the steps:
- a
- Allow specification of relevant parameters in a straightforward C++ script
- Pass these parameters to FeynHiggs (FH) and SUSY-HIT (SH), then run both
- b
- Extract output from FH and SH, combining relevant blocks
- Higgs masses and decays from FH
- Sparticle properties and decay from SH
- Define Sushi block
- Is there a parameter dependance, or will we use the same block for all inputs?
- Extract output from FH and SH, combining relevant blocks