example_slha.py.in

Go to the documentation of this file.

#!/usr/bin/env python
 
from __future__ import print_function
from gm2_python_interface import *
 
cppyy.include(os.path.join("gm2calc","gm2_1loop.hpp"))
cppyy.include(os.path.join("gm2calc","gm2_2loop.hpp"))
cppyy.include(os.path.join("gm2calc","gm2_uncertainty.hpp"))
cppyy.include(os.path.join("gm2calc","gm2_error.hpp"))
cppyy.include(os.path.join("gm2calc","MSSMNoFV_onshell.hpp"))
 
cppyy.load_library("libgm2calc")
    
# Load data types
from cppyy.gbl import gm2calc
from cppyy.gbl.gm2calc import Error
 
def setup():
    # load data types
    from cppyy.gbl import std
    from cppyy.gbl import Eigen
    # Shorten class call name
    from cppyy.gbl.Eigen import Matrix3d
    
    model = gm2calc.MSSMNoFV_onshell()
    
    Pi = 3.141592653589793
    # Outer Matrix3cd is to convert type from CwiseNullaryOp to Matrix
    UnitMatrix = Matrix3d(Matrix3d().Identity())
    # __mul__ is not defined for Matrix type
    UnitMatrix.__imul__(7000*7000)
    
    # fill SM parameters
    model.set_alpha_MZ(0.0077552)               # 1L
    model.set_alpha_thompson(0.00729735)        # 2L
    model.set_g3(std.sqrt(4. * Pi * 0.1184))    # 2L
    model.get_physical().MFt   = 173.34         # 2L
    model.get_physical().MFb   = 4.18           # 2L, mb(mb) MS-bar
    model.get_physical().MFm   = 0.1056583715   # 1L
    model.get_physical().MFtau = 1.777          # 2L
    model.get_physical().MVWm  = 80.385         # 1L
    model.get_physical().MVZ   = 91.1876        # 1L
    
    # fill pole masses
    model.get_physical().MSvmL   =  5.18860573e+02 # 1L
    model.get_physical().MSm[0]  =  5.05095249e+02 # 1L
    model.get_physical().MSm[1]  =  5.25187016e+02 # 1L
    model.get_physical().MChi[0] =  2.01611468e+02 # 1L
    model.get_physical().MChi[1] =  4.10040273e+02 # 1L
    model.get_physical().MChi[2] = -5.16529941e+02 # 1L
    model.get_physical().MChi[3] =  5.45628749e+02 # 1L
    model.get_physical().MCha[0] =  4.09989890e+02 # 1L
    model.get_physical().MCha[1] =  5.46057190e+02 # 1L
    model.get_physical().MAh[1]  =  1.50000000e+03 # 2L
    
    # fill DR-bar parameters
    model.set_TB(40)                        # 1L
    model.set_Mu(500)                       # initial guess
    model.set_MassB(200)                    # initial guess
    model.set_MassWB(400)                   # initial guess
    model.set_MassG(2000)                   # 2L
    model.set_mq2(UnitMatrix)               # 2L
    model.set_ml2(0, 0, 500 * 500)          # 2L
    model.set_ml2(1, 1, 500 * 500)          # irrelevant
    model.set_ml2(2, 2, 500 * 500)          # 2L
    model.set_md2(UnitMatrix)               # 2L
    model.set_mu2(UnitMatrix)               # 2L
    model.set_me2(0, 0, 500 * 500)          # 2L
    model.set_me2(1, 1, 500 * 500)          # initial guess
    model.set_me2(2, 2, 500 * 500)          # 2L
    model.set_Au(2, 2, 0)                   # 2L
    model.set_Ad(2, 2, 0)                   # 2L
    model.set_Ae(1, 1, 0)                   # 1L
    model.set_Ae(2, 2, 0)                   # 2L
    model.set_scale(1000)                   # 2L
    
    # convert DR-bar parameters to on-shell
    model.convert_to_onshell()
    
    # check for warnings
    if model.get_problems().have_warning():
        print(model.get_problems().get_warnings())
    
    # check for problems
    if model.get_problems().have_problem():
        print(model.get_problems().get_problems())
    
    return model
 
 
try:
    model = setup()
    amu = gm2calc.calculate_amu_1loop(model) + gm2calc.calculate_amu_2loop(model)
    delta_amu = gm2calc.calculate_uncertainty_amu_2loop(model)
    print("amu =",amu,"+-",delta_amu)
except gm2calc.Error as e:
    print(e.what())