complex parameters not exported to UFO in feynrules

Asked by Ilaria Brivio

Hi,

I seem to have a very weird problem with Feynrules. I realized that none of the UFO models exported for SMEFTsim-set A
http://feynrules.irmp.ucl.ac.be/wiki/SMEFT
contain the complex conjugates of complex Wilson coefficients.

I am running some tests with the latest version of Feynrules 2.3.29 and the operator OuW in the U3^5 symmetric model.

The Lagrangian is just

cuW OuW + Conjugate[cuW] HC[OuW]

with these definitions:

[from SMEFTsim_A_operators_U35.fr]

OuW:= Block[{sp1,sp2,mu,nu,ii,jj1,jj2,kk,cc,ff1,ff2},
 feynmangaugerules = If[Not[FeynmanGauge], {G0|GP|GPbar ->0}, {}];
 ExpandIndices[QLmbar[sp1,ii,ff1,cc].uR[sp2,ff2,cc] Conjugate[yu[ff2, ff1]] sigmaT[mu,nu,sp1,sp2] PauliSigma[kk,ii,jj1] Eps[jj1,jj2] Phibar[jj2] FS[Wi,mu,nu,kk] , FlavorExpand->{SU2W,SU2D}]/.feynmangaugerules/.redefCtoZero
 ]//OptimizeIndex;

[from d6_parameters_U35.fr]

  cuWPh == {ParameterType -> External, Value -> 0, ComplexParameter -> False},

  cuWAbs == {
    ParameterType -> External,
    Value -> 1,
    InteractionOrder -> {NP, 1}
  },

  cuW == {
    ParameterType -> Internal,
    Value -> cuWAbs Exp[I cuWPh],
    ComplexParameter -> True,
    InteractionOrder -> {NP, 1},
    TeX -> Subscript[c,uW]
  },

when I look at the model in Mathematica I see that cuW is complex and the vertices given by FeynmanRules are correct.

However, when I do ExportUFO something goes wrong (see files attached).
cuW is actually treated as a real quantity in the couplings, although the parameter.py file says that it's a complex quantity.

The problem appears to be solved if I write explicitly cuWAbs Exp[I cuWPh] in the Lagrangian, but expanding all the complex coefficients would make the Lagrangian manipulation quite heavy, so I'd prefer to keep the notation compact if possible. (besides, SMEFTsim-setB is coded in a similar way but it doesn't appear to have the same issue, which is quite weird)

Is there something fishy I am doing wrong?

Thanks a lot!
Ilaria

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Ilaria Brivio
Solved:
Last query:
Last reply:
Revision history for this message
Ilaria Brivio (ilariab) said :
#1

This is the couplings.py I get for the exported UFO

# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.0.1 for Linux x86 (64-bit) (September 21, 2016)
# Date: Thu 18 Jan 2018 11:50:52

from object_library import all_couplings, Coupling

from function_library import complexconjugate, re, im, csc, sec, acsc, asec, cot

GC_1 = Coupling(name = 'GC_1',
                value = 'cuW*complex(0,1)*yc',
                order = {'NP':1,'QED':1})

GC_2 = Coupling(name = 'GC_2',
                value = '(cth*cuW*complex(0,1)*yc)/cmath.sqrt(2)',
                order = {'NP':1,'QED':1})

GC_3 = Coupling(name = 'GC_3',
                value = '-(cuW*ee*complex(0,1)*yc)',
                order = {'NP':1,'QED':2})

GC_4 = Coupling(name = 'GC_4',
                value = 'cuW*ee*complex(0,1)*yc',
                order = {'NP':1,'QED':2})

GC_5 = Coupling(name = 'GC_5',
                value = '-((cuW*ee*complex(0,1)*yc)/(sth*cmath.sqrt(2)))',
                order = {'NP':1,'QED':2})

GC_6 = Coupling(name = 'GC_6',
                value = '-((cth*cuW*ee*complex(0,1)*yc)/sth)',
                order = {'NP':1,'QED':2})

GC_7 = Coupling(name = 'GC_7',
                value = '(cth*cuW*ee*complex(0,1)*yc)/sth',
                order = {'NP':1,'QED':2})

GC_8 = Coupling(name = 'GC_8',
                value = '(cuW*complex(0,1)*sth*yc)/cmath.sqrt(2)',
                order = {'NP':1,'QED':1})

GC_9 = Coupling(name = 'GC_9',
                value = 'cuW*complex(0,1)*vevhat*yc',
                order = {'NP':1})

GC_10 = Coupling(name = 'GC_10',
                 value = '(cth*cuW*complex(0,1)*vevhat*yc)/cmath.sqrt(2)',
                 order = {'NP':1})

GC_11 = Coupling(name = 'GC_11',
                 value = '-(cuW*ee*complex(0,1)*vevhat*yc)',
                 order = {'NP':1,'QED':1})

GC_12 = Coupling(name = 'GC_12',
                 value = 'cuW*ee*complex(0,1)*vevhat*yc',
                 order = {'NP':1,'QED':1})

GC_13 = Coupling(name = 'GC_13',
                 value = '-((cuW*ee*complex(0,1)*vevhat*yc)/(sth*cmath.sqrt(2)))',
                 order = {'NP':1,'QED':1})

GC_14 = Coupling(name = 'GC_14',
                 value = '-((cth*cuW*ee*complex(0,1)*vevhat*yc)/sth)',
                 order = {'NP':1,'QED':1})

GC_15 = Coupling(name = 'GC_15',
                 value = '(cth*cuW*ee*complex(0,1)*vevhat*yc)/sth',
                 order = {'NP':1,'QED':1})

GC_16 = Coupling(name = 'GC_16',
                 value = '(cuW*complex(0,1)*sth*vevhat*yc)/cmath.sqrt(2)',
                 order = {'NP':1})

GC_17 = Coupling(name = 'GC_17',
                 value = 'cuW*complex(0,1)*yt',
                 order = {'NP':1,'QED':1})

GC_18 = Coupling(name = 'GC_18',
                 value = '(cth*cuW*complex(0,1)*yt)/cmath.sqrt(2)',
                 order = {'NP':1,'QED':1})

GC_19 = Coupling(name = 'GC_19',
                 value = '-(cuW*ee*complex(0,1)*yt)',
                 order = {'NP':1,'QED':2})

GC_20 = Coupling(name = 'GC_20',
                 value = 'cuW*ee*complex(0,1)*yt',
                 order = {'NP':1,'QED':2})

GC_21 = Coupling(name = 'GC_21',
                 value = '-((cuW*ee*complex(0,1)*yt)/(sth*cmath.sqrt(2)))',
                 order = {'NP':1,'QED':2})

GC_22 = Coupling(name = 'GC_22',
                 value = '-((cth*cuW*ee*complex(0,1)*yt)/sth)',
                 order = {'NP':1,'QED':2})

GC_23 = Coupling(name = 'GC_23',
                 value = '(cth*cuW*ee*complex(0,1)*yt)/sth',
                 order = {'NP':1,'QED':2})

GC_24 = Coupling(name = 'GC_24',
                 value = '(cuW*complex(0,1)*sth*yt)/cmath.sqrt(2)',
                 order = {'NP':1,'QED':1})

GC_25 = Coupling(name = 'GC_25',
                 value = 'cuW*complex(0,1)*vevhat*yt',
                 order = {'NP':1})

GC_26 = Coupling(name = 'GC_26',
                 value = '(cth*cuW*complex(0,1)*vevhat*yt)/cmath.sqrt(2)',
                 order = {'NP':1})

GC_27 = Coupling(name = 'GC_27',
                 value = '-(cuW*ee*complex(0,1)*vevhat*yt)',
                 order = {'NP':1,'QED':1})

GC_28 = Coupling(name = 'GC_28',
                 value = 'cuW*ee*complex(0,1)*vevhat*yt',
                 order = {'NP':1,'QED':1})

GC_29 = Coupling(name = 'GC_29',
                 value = '-((cuW*ee*complex(0,1)*vevhat*yt)/(sth*cmath.sqrt(2)))',
                 order = {'NP':1,'QED':1})

GC_30 = Coupling(name = 'GC_30',
                 value = '-((cth*cuW*ee*complex(0,1)*vevhat*yt)/sth)',
                 order = {'NP':1,'QED':1})

GC_31 = Coupling(name = 'GC_31',
                 value = '(cth*cuW*ee*complex(0,1)*vevhat*yt)/sth',
                 order = {'NP':1,'QED':1})

GC_32 = Coupling(name = 'GC_32',
                 value = '(cuW*complex(0,1)*sth*vevhat*yt)/cmath.sqrt(2)',
                 order = {'NP':1})

GC_33 = Coupling(name = 'GC_33',
                 value = 'cuW*complex(0,1)*yup',
                 order = {'NP':1,'QED':1})

GC_34 = Coupling(name = 'GC_34',
                 value = '(cth*cuW*complex(0,1)*yup)/cmath.sqrt(2)',
                 order = {'NP':1,'QED':1})

GC_35 = Coupling(name = 'GC_35',
                 value = '-(cuW*ee*complex(0,1)*yup)',
                 order = {'NP':1,'QED':2})

GC_36 = Coupling(name = 'GC_36',
                 value = 'cuW*ee*complex(0,1)*yup',
                 order = {'NP':1,'QED':2})

GC_37 = Coupling(name = 'GC_37',
                 value = '-((cuW*ee*complex(0,1)*yup)/(sth*cmath.sqrt(2)))',
                 order = {'NP':1,'QED':2})

GC_38 = Coupling(name = 'GC_38',
                 value = '(cth*cuW*ee*complex(0,1)*yup)/sth',
                 order = {'NP':1,'QED':2})

GC_39 = Coupling(name = 'GC_39',
                 value = '(cuW*complex(0,1)*sth*yup)/cmath.sqrt(2)',
                 order = {'NP':1,'QED':1})

GC_40 = Coupling(name = 'GC_40',
                 value = 'cuW*complex(0,1)*vevhat*yup',
                 order = {'NP':1})

GC_41 = Coupling(name = 'GC_41',
                 value = '(cth*cuW*complex(0,1)*vevhat*yup)/cmath.sqrt(2)',
                 order = {'NP':1})

GC_42 = Coupling(name = 'GC_42',
                 value = '-(cuW*ee*complex(0,1)*vevhat*yup)',
                 order = {'NP':1,'QED':1})

GC_43 = Coupling(name = 'GC_43',
                 value = 'cuW*ee*complex(0,1)*vevhat*yup',
                 order = {'NP':1,'QED':1})

GC_44 = Coupling(name = 'GC_44',
                 value = '-((cuW*ee*complex(0,1)*vevhat*yup)/(sth*cmath.sqrt(2)))',
                 order = {'NP':1,'QED':1})

GC_45 = Coupling(name = 'GC_45',
                 value = '-((cth*cuW*ee*complex(0,1)*vevhat*yup)/sth)',
                 order = {'NP':1,'QED':1})

GC_46 = Coupling(name = 'GC_46',
                 value = '(cth*cuW*ee*complex(0,1)*vevhat*yup)/sth',
                 order = {'NP':1,'QED':1})

GC_47 = Coupling(name = 'GC_47',
                 value = '(cuW*complex(0,1)*sth*vevhat*yup)/cmath.sqrt(2)',
                 order = {'NP':1})

Revision history for this message
Ilaria Brivio (ilariab) said :
#2

and the parameters.py

# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.0.1 for Linux x86 (64-bit) (September 21, 2016)
# Date: Thu 18 Jan 2018 11:50:52

from object_library import all_parameters, Parameter

from function_library import complexconjugate, re, im, csc, sec, acsc, asec, cot

# This is a default parameter object representing 0.
ZERO = Parameter(name = 'ZERO',
                 nature = 'internal',
                 type = 'real',
                 value = '0.0',
                 texname = '0')

# User-defined parameters.
cabi = Parameter(name = 'cabi',
                 nature = 'external',
                 type = 'real',
                 value = 0.227736,
                 texname = '\\theta _c',
                 lhablock = 'CKMBLOCK',
                 lhacode = [ 1 ])

CKMlambda = Parameter(name = 'CKMlambda',
                      nature = 'external',
                      type = 'real',
                      value = 0.22506,
                      texname = '\\text{CKMlambda}',
                      lhablock = 'CKMBLOCK',
                      lhacode = [ 2 ])

CKMA = Parameter(name = 'CKMA',
                 nature = 'external',
                 type = 'real',
                 value = 0.811,
                 texname = '\\text{CKMA}',
                 lhablock = 'CKMBLOCK',
                 lhacode = [ 3 ])

CKMrho = Parameter(name = 'CKMrho',
                   nature = 'external',
                   type = 'real',
                   value = 0.124,
                   texname = '\\text{CKMrho}',
                   lhablock = 'CKMBLOCK',
                   lhacode = [ 4 ])

CKMeta = Parameter(name = 'CKMeta',
                   nature = 'external',
                   type = 'real',
                   value = 0.356,
                   texname = '\\text{CKMeta}',
                   lhablock = 'CKMBLOCK',
                   lhacode = [ 5 ])

aEW = Parameter(name = 'aEW',
                nature = 'external',
                type = 'real',
                value = 0.007815552950371238,
                texname = '\\alpha _{\\text{EW}}',
                lhablock = 'SMINPUTS',
                lhacode = [ 1 ])

Gf = Parameter(name = 'Gf',
               nature = 'external',
               type = 'real',
               value = 0.000011663787,
               texname = 'G_f',
               lhablock = 'SMINPUTS',
               lhacode = [ 2 ])

aS = Parameter(name = 'aS',
               nature = 'external',
               type = 'real',
               value = 0.1181,
               texname = '\\alpha _s',
               lhablock = 'SMINPUTS',
               lhacode = [ 3 ])

ymdo = Parameter(name = 'ymdo',
                 nature = 'external',
                 type = 'real',
                 value = 0.0047,
                 texname = '\\text{ymdo}',
                 lhablock = 'YUKAWA',
                 lhacode = [ 1 ])

ymup = Parameter(name = 'ymup',
                 nature = 'external',
                 type = 'real',
                 value = 0.0022,
                 texname = '\\text{ymup}',
                 lhablock = 'YUKAWA',
                 lhacode = [ 2 ])

yms = Parameter(name = 'yms',
                nature = 'external',
                type = 'real',
                value = 0.096,
                texname = '\\text{yms}',
                lhablock = 'YUKAWA',
                lhacode = [ 3 ])

ymc = Parameter(name = 'ymc',
                nature = 'external',
                type = 'real',
                value = 1.28,
                texname = '\\text{ymc}',
                lhablock = 'YUKAWA',
                lhacode = [ 4 ])

ymb = Parameter(name = 'ymb',
                nature = 'external',
                type = 'real',
                value = 4.18,
                texname = '\\text{ymb}',
                lhablock = 'YUKAWA',
                lhacode = [ 5 ])

ymt = Parameter(name = 'ymt',
                nature = 'external',
                type = 'real',
                value = 173.2,
                texname = '\\text{ymt}',
                lhablock = 'YUKAWA',
                lhacode = [ 6 ])

yme = Parameter(name = 'yme',
                nature = 'external',
                type = 'real',
                value = 0.000511,
                texname = '\\text{yme}',
                lhablock = 'YUKAWA',
                lhacode = [ 11 ])

ymm = Parameter(name = 'ymm',
                nature = 'external',
                type = 'real',
                value = 0.10566,
                texname = '\\text{ymm}',
                lhablock = 'YUKAWA',
                lhacode = [ 13 ])

ymtau = Parameter(name = 'ymtau',
                  nature = 'external',
                  type = 'real',
                  value = 1.777,
                  texname = '\\text{ymtau}',
                  lhablock = 'YUKAWA',
                  lhacode = [ 15 ])

ceWPh = Parameter(name = 'ceWPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{ceWPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 1 ])

ceBPh = Parameter(name = 'ceBPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{ceBPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 2 ])

cuGPh = Parameter(name = 'cuGPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cuGPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 3 ])

cuWPh = Parameter(name = 'cuWPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cuWPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 4 ])

cuBPh = Parameter(name = 'cuBPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cuBPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 5 ])

cdGPh = Parameter(name = 'cdGPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cdGPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 6 ])

cdWPh = Parameter(name = 'cdWPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cdWPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 7 ])

cdBPh = Parameter(name = 'cdBPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cdBPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 8 ])

cHudPh = Parameter(name = 'cHudPh',
                   nature = 'external',
                   type = 'real',
                   value = 0,
                   texname = '\\text{cHudPh}',
                   lhablock = 'FRBlock',
                   lhacode = [ 9 ])

ceHPh = Parameter(name = 'ceHPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{ceHPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 10 ])

cuHPh = Parameter(name = 'cuHPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cuHPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 11 ])

cdHPh = Parameter(name = 'cdHPh',
                  nature = 'external',
                  type = 'real',
                  value = 0,
                  texname = '\\text{cdHPh}',
                  lhablock = 'FRBlock',
                  lhacode = [ 12 ])

cledqPh = Parameter(name = 'cledqPh',
                    nature = 'external',
                    type = 'real',
                    value = 0,
                    texname = '\\text{cledqPh}',
                    lhablock = 'FRBlock',
                    lhacode = [ 13 ])

cquqd1Ph = Parameter(name = 'cquqd1Ph',
                     nature = 'external',
                     type = 'real',
                     value = 0,
                     texname = '\\text{cquqd1Ph}',
                     lhablock = 'FRBlock',
                     lhacode = [ 14 ])

cquqd8Ph = Parameter(name = 'cquqd8Ph',
                     nature = 'external',
                     type = 'real',
                     value = 0,
                     texname = '\\text{cquqd8Ph}',
                     lhablock = 'FRBlock',
                     lhacode = [ 15 ])

clequ1Ph = Parameter(name = 'clequ1Ph',
                     nature = 'external',
                     type = 'real',
                     value = 0,
                     texname = '\\text{clequ1Ph}',
                     lhablock = 'FRBlock',
                     lhacode = [ 16 ])

clequ3Ph = Parameter(name = 'clequ3Ph',
                     nature = 'external',
                     type = 'real',
                     value = 0,
                     texname = '\\text{clequ3Ph}',
                     lhablock = 'FRBlock',
                     lhacode = [ 17 ])

LambdaSMEFT = Parameter(name = 'LambdaSMEFT',
                        nature = 'external',
                        type = 'real',
                        value = 1000,
                        texname = '\\Lambda',
                        lhablock = 'FRBlock',
                        lhacode = [ 18 ])

cG = Parameter(name = 'cG',
               nature = 'external',
               type = 'real',
               value = 1,
               texname = 'c_G',
               lhablock = 'FRBlock',
               lhacode = [ 19 ])

cGtil = Parameter(name = 'cGtil',
                  nature = 'external',
                  type = 'real',
                  value = 1,
                  texname = 'c_{\\tilde{G}}',
                  lhablock = 'FRBlock',
                  lhacode = [ 20 ])

cW = Parameter(name = 'cW',
               nature = 'external',
               type = 'real',
               value = 1,
               texname = 'c_W',
               lhablock = 'FRBlock',
               lhacode = [ 21 ])

cWtil = Parameter(name = 'cWtil',
                  nature = 'external',
                  type = 'real',
                  value = 1,
                  texname = 'c_{\\tilde{W}}',
                  lhablock = 'FRBlock',
                  lhacode = [ 22 ])

cH = Parameter(name = 'cH',
               nature = 'external',
               type = 'real',
               value = 1,
               texname = 'c_H',
               lhablock = 'FRBlock',
               lhacode = [ 23 ])

cHbox = Parameter(name = 'cHbox',
                  nature = 'external',
                  type = 'real',
                  value = 1,
                  texname = 'c_{H \\square }',
                  lhablock = 'FRBlock',
                  lhacode = [ 24 ])

cHD = Parameter(name = 'cHD',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{HD}}',
                lhablock = 'FRBlock',
                lhacode = [ 25 ])

cHG = Parameter(name = 'cHG',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{HG}}',
                lhablock = 'FRBlock',
                lhacode = [ 26 ])

cHGtil = Parameter(name = 'cHGtil',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = 'c_{H \\tilde{G}}',
                   lhablock = 'FRBlock',
                   lhacode = [ 27 ])

cHW = Parameter(name = 'cHW',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{HW}}',
                lhablock = 'FRBlock',
                lhacode = [ 28 ])

cHWtil = Parameter(name = 'cHWtil',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = 'c_{H \\tilde{W}}',
                   lhablock = 'FRBlock',
                   lhacode = [ 29 ])

cHB = Parameter(name = 'cHB',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{HB}}',
                lhablock = 'FRBlock',
                lhacode = [ 30 ])

cHBtil = Parameter(name = 'cHBtil',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = 'c_{H \\tilde{B}}',
                   lhablock = 'FRBlock',
                   lhacode = [ 31 ])

cHWB = Parameter(name = 'cHWB',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = 'c_{\\text{HWB}}',
                 lhablock = 'FRBlock',
                 lhacode = [ 32 ])

cHWBtil = Parameter(name = 'cHWBtil',
                    nature = 'external',
                    type = 'real',
                    value = 1,
                    texname = 'c_{B H \\tilde{W}}',
                    lhablock = 'FRBlock',
                    lhacode = [ 33 ])

ceHAbs = Parameter(name = 'ceHAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{ceHAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 34 ])

cuHAbs = Parameter(name = 'cuHAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cuHAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 35 ])

cdHAbs = Parameter(name = 'cdHAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cdHAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 36 ])

ceWAbs = Parameter(name = 'ceWAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{ceWAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 37 ])

ceBAbs = Parameter(name = 'ceBAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{ceBAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 38 ])

cuGAbs = Parameter(name = 'cuGAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cuGAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 39 ])

cuWAbs = Parameter(name = 'cuWAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cuWAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 40 ])

cuBAbs = Parameter(name = 'cuBAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cuBAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 41 ])

cdGAbs = Parameter(name = 'cdGAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cdGAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 42 ])

cdWAbs = Parameter(name = 'cdWAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cdWAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 43 ])

cdBAbs = Parameter(name = 'cdBAbs',
                   nature = 'external',
                   type = 'real',
                   value = 1,
                   texname = '\\text{cdBAbs}',
                   lhablock = 'FRBlock',
                   lhacode = [ 44 ])

cHl1 = Parameter(name = 'cHl1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{Hl},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 45 ])

cHl3 = Parameter(name = 'cHl3',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{Hl},3]',
                 lhablock = 'FRBlock',
                 lhacode = [ 46 ])

cHe = Parameter(name = 'cHe',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{He}}',
                lhablock = 'FRBlock',
                lhacode = [ 47 ])

cHq1 = Parameter(name = 'cHq1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{Hq},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 48 ])

cHq3 = Parameter(name = 'cHq3',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{Hq},3]',
                 lhablock = 'FRBlock',
                 lhacode = [ 49 ])

cHu = Parameter(name = 'cHu',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{Hu}}',
                lhablock = 'FRBlock',
                lhacode = [ 50 ])

cHd = Parameter(name = 'cHd',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{Hd}}',
                lhablock = 'FRBlock',
                lhacode = [ 51 ])

cHudAbs = Parameter(name = 'cHudAbs',
                    nature = 'external',
                    type = 'real',
                    value = 1,
                    texname = '\\text{cHudAbs}',
                    lhablock = 'FRBlock',
                    lhacode = [ 52 ])

cll = Parameter(name = 'cll',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{ll}}',
                lhablock = 'FRBlock',
                lhacode = [ 53 ])

cll1 = Parameter(name = 'cll1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{ll},\\text{Prime}]',
                 lhablock = 'FRBlock',
                 lhacode = [ 54 ])

cqq1 = Parameter(name = 'cqq1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qq},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 55 ])

cqq11 = Parameter(name = 'cqq11',
                  nature = 'external',
                  type = 'real',
                  value = 1,
                  texname = '\\text{Subsuperscript}[c,\\text{qq},\\text{Prime}]',
                  lhablock = 'FRBlock',
                  lhacode = [ 56 ])

cqq3 = Parameter(name = 'cqq3',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qq},3]',
                 lhablock = 'FRBlock',
                 lhacode = [ 57 ])

cqq31 = Parameter(name = 'cqq31',
                  nature = 'external',
                  type = 'real',
                  value = 1,
                  texname = '\\text{Subsuperscript}[c,\\text{qq},3 \\text{Prime}]',
                  lhablock = 'FRBlock',
                  lhacode = [ 58 ])

clq1 = Parameter(name = 'clq1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{lq},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 59 ])

clq3 = Parameter(name = 'clq3',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{lq},3]',
                 lhablock = 'FRBlock',
                 lhacode = [ 60 ])

cee = Parameter(name = 'cee',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{e^2}',
                lhablock = 'FRBlock',
                lhacode = [ 61 ])

cuu = Parameter(name = 'cuu',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{uu}}',
                lhablock = 'FRBlock',
                lhacode = [ 62 ])

cuu1 = Parameter(name = 'cuu1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{uu},\\text{Prime}]',
                 lhablock = 'FRBlock',
                 lhacode = [ 63 ])

cdd = Parameter(name = 'cdd',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{dd}}',
                lhablock = 'FRBlock',
                lhacode = [ 64 ])

cdd1 = Parameter(name = 'cdd1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{dd},\\text{Prime}]',
                 lhablock = 'FRBlock',
                 lhacode = [ 65 ])

ceu = Parameter(name = 'ceu',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{eu}}',
                lhablock = 'FRBlock',
                lhacode = [ 66 ])

ced = Parameter(name = 'ced',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{ed}}',
                lhablock = 'FRBlock',
                lhacode = [ 67 ])

cud1 = Parameter(name = 'cud1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{ud},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 68 ])

cud8 = Parameter(name = 'cud8',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{ud},8]',
                 lhablock = 'FRBlock',
                 lhacode = [ 69 ])

cle = Parameter(name = 'cle',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{le}}',
                lhablock = 'FRBlock',
                lhacode = [ 70 ])

clu = Parameter(name = 'clu',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{lu}}',
                lhablock = 'FRBlock',
                lhacode = [ 71 ])

cld = Parameter(name = 'cld',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{ld}}',
                lhablock = 'FRBlock',
                lhacode = [ 72 ])

cqe = Parameter(name = 'cqe',
                nature = 'external',
                type = 'real',
                value = 1,
                texname = 'c_{\\text{qe}}',
                lhablock = 'FRBlock',
                lhacode = [ 73 ])

cqu1 = Parameter(name = 'cqu1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qu},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 74 ])

cqu8 = Parameter(name = 'cqu8',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qu},8]',
                 lhablock = 'FRBlock',
                 lhacode = [ 75 ])

cqd1 = Parameter(name = 'cqd1',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qd},1]',
                 lhablock = 'FRBlock',
                 lhacode = [ 76 ])

cqd8 = Parameter(name = 'cqd8',
                 nature = 'external',
                 type = 'real',
                 value = 1,
                 texname = '\\text{Subsuperscript}[c,\\text{qd},8]',
                 lhablock = 'FRBlock',
                 lhacode = [ 77 ])

cledqAbs = Parameter(name = 'cledqAbs',
                     nature = 'external',
                     type = 'real',
                     value = 1,
                     texname = '\\text{cledqAbs}',
                     lhablock = 'FRBlock',
                     lhacode = [ 78 ])

cquqd1Abs = Parameter(name = 'cquqd1Abs',
                      nature = 'external',
                      type = 'real',
                      value = 1,
                      texname = '\\text{cquqd1Abs}',
                      lhablock = 'FRBlock',
                      lhacode = [ 79 ])

cquqd8Abs = Parameter(name = 'cquqd8Abs',
                      nature = 'external',
                      type = 'real',
                      value = 1,
                      texname = '\\text{cquqd8Abs}',
                      lhablock = 'FRBlock',
                      lhacode = [ 80 ])

clequ1Abs = Parameter(name = 'clequ1Abs',
                      nature = 'external',
                      type = 'real',
                      value = 1,
                      texname = '\\text{clequ1Abs}',
                      lhablock = 'FRBlock',
                      lhacode = [ 81 ])

clequ3Abs = Parameter(name = 'clequ3Abs',
                      nature = 'external',
                      type = 'real',
                      value = 1,
                      texname = '\\text{clequ3Abs}',
                      lhablock = 'FRBlock',
                      lhacode = [ 82 ])

Hloopcouplings = Parameter(name = 'Hloopcouplings',
                           nature = 'external',
                           type = 'real',
                           value = 1,
                           texname = '\\eta _{H \\text{loops}}',
                           lhablock = 'FRBlock',
                           lhacode = [ 83 ])

MZ = Parameter(name = 'MZ',
               nature = 'external',
               type = 'real',
               value = 91.1876,
               texname = '\\text{MZ}',
               lhablock = 'MASS',
               lhacode = [ 23 ])

Me = Parameter(name = 'Me',
               nature = 'external',
               type = 'real',
               value = 0.000511,
               texname = '\\text{Me}',
               lhablock = 'MASS',
               lhacode = [ 11 ])

MMU = Parameter(name = 'MMU',
                nature = 'external',
                type = 'real',
                value = 0.10566,
                texname = '\\text{MMU}',
                lhablock = 'MASS',
                lhacode = [ 13 ])

MTA = Parameter(name = 'MTA',
                nature = 'external',
                type = 'real',
                value = 1.777,
                texname = '\\text{MTA}',
                lhablock = 'MASS',
                lhacode = [ 15 ])

MU = Parameter(name = 'MU',
               nature = 'external',
               type = 'real',
               value = 0.0022,
               texname = 'M',
               lhablock = 'MASS',
               lhacode = [ 2 ])

MC = Parameter(name = 'MC',
               nature = 'external',
               type = 'real',
               value = 1.28,
               texname = '\\text{MC}',
               lhablock = 'MASS',
               lhacode = [ 4 ])

MT = Parameter(name = 'MT',
               nature = 'external',
               type = 'real',
               value = 173.2,
               texname = '\\text{MT}',
               lhablock = 'MASS',
               lhacode = [ 6 ])

MD = Parameter(name = 'MD',
               nature = 'external',
               type = 'real',
               value = 0.0047,
               texname = '\\text{MD}',
               lhablock = 'MASS',
               lhacode = [ 1 ])

MS = Parameter(name = 'MS',
               nature = 'external',
               type = 'real',
               value = 0.096,
               texname = '\\text{MS}',
               lhablock = 'MASS',
               lhacode = [ 3 ])

MB = Parameter(name = 'MB',
               nature = 'external',
               type = 'real',
               value = 4.18,
               texname = '\\text{MB}',
               lhablock = 'MASS',
               lhacode = [ 5 ])

MH = Parameter(name = 'MH',
               nature = 'external',
               type = 'real',
               value = 125.09,
               texname = '\\text{MH}',
               lhablock = 'MASS',
               lhacode = [ 25 ])

WZ = Parameter(name = 'WZ',
               nature = 'external',
               type = 'real',
               value = 2.4952,
               texname = '\\text{WZ}',
               lhablock = 'DECAY',
               lhacode = [ 23 ])

WW = Parameter(name = 'WW',
               nature = 'external',
               type = 'real',
               value = 2.085,
               texname = '\\text{WW}',
               lhablock = 'DECAY',
               lhacode = [ 24 ])

WT = Parameter(name = 'WT',
               nature = 'external',
               type = 'real',
               value = 1.50833649,
               texname = '\\text{WT}',
               lhablock = 'DECAY',
               lhacode = [ 6 ])

WH = Parameter(name = 'WH',
               nature = 'external',
               type = 'real',
               value = 0.00407,
               texname = '\\text{WH}',
               lhablock = 'DECAY',
               lhacode = [ 25 ])

ceH = Parameter(name = 'ceH',
                nature = 'internal',
                type = 'complex',
                value = 'ceHAbs*cmath.exp(ceHPh*complex(0,1))',
                texname = 'c_{\\text{eH}}')

cuH = Parameter(name = 'cuH',
                nature = 'internal',
                type = 'complex',
                value = 'cuHAbs*cmath.exp(cuHPh*complex(0,1))',
                texname = 'c_{\\text{uH}}')

cdH = Parameter(name = 'cdH',
                nature = 'internal',
                type = 'complex',
                value = 'cdHAbs*cmath.exp(cdHPh*complex(0,1))',
                texname = 'c_{\\text{dH}}')

ceW = Parameter(name = 'ceW',
                nature = 'internal',
                type = 'complex',
                value = 'ceWAbs*cmath.exp(ceWPh*complex(0,1))',
                texname = 'c_{\\text{eW}}')

ceB = Parameter(name = 'ceB',
                nature = 'internal',
                type = 'complex',
                value = 'ceBAbs*cmath.exp(ceBPh*complex(0,1))',
                texname = 'c_{\\text{eB}}')

cuG = Parameter(name = 'cuG',
                nature = 'internal',
                type = 'complex',
                value = 'cuGAbs*cmath.exp(cuGPh*complex(0,1))',
                texname = 'c_{\\text{uG}}')

cuW = Parameter(name = 'cuW',
                nature = 'internal',
                type = 'complex',
                value = 'cuWAbs*cmath.exp(cuWPh*complex(0,1))',
                texname = 'c_{\\text{uW}}')

cuB = Parameter(name = 'cuB',
                nature = 'internal',
                type = 'complex',
                value = 'cuBAbs*cmath.exp(cuBPh*complex(0,1))',
                texname = 'c_{\\text{uB}}')

cdG = Parameter(name = 'cdG',
                nature = 'internal',
                type = 'complex',
                value = 'cdGAbs*cmath.exp(cdGPh*complex(0,1))',
                texname = 'c_{\\text{dG}}')

cdW = Parameter(name = 'cdW',
                nature = 'internal',
                type = 'complex',
                value = 'cdWAbs*cmath.exp(cdWPh*complex(0,1))',
                texname = 'c_{\\text{dW}}')

cdB = Parameter(name = 'cdB',
                nature = 'internal',
                type = 'complex',
                value = 'cdBAbs*cmath.exp(cdBPh*complex(0,1))',
                texname = 'c_{\\text{dB}}')

cHud = Parameter(name = 'cHud',
                 nature = 'internal',
                 type = 'complex',
                 value = 'cHudAbs*cmath.exp(cHudPh*complex(0,1))',
                 texname = 'c_{\\text{Hud}}')

cledq = Parameter(name = 'cledq',
                  nature = 'internal',
                  type = 'complex',
                  value = 'cledqAbs*cmath.exp(cledqPh*complex(0,1))',
                  texname = 'c_{\\text{ledq}}')

cquqd1 = Parameter(name = 'cquqd1',
                   nature = 'internal',
                   type = 'complex',
                   value = 'cquqd1Abs*cmath.exp(cquqd1Ph*complex(0,1))',
                   texname = '\\text{Subsuperscript}[c,\\text{quqd},1]')

cquqd8 = Parameter(name = 'cquqd8',
                   nature = 'internal',
                   type = 'complex',
                   value = 'cquqd8Abs*cmath.exp(cquqd8Ph*complex(0,1))',
                   texname = '\\text{Subsuperscript}[c,\\text{quqd},8]')

clequ1 = Parameter(name = 'clequ1',
                   nature = 'internal',
                   type = 'complex',
                   value = 'clequ1Abs*cmath.exp(clequ1Ph*complex(0,1))',
                   texname = '\\text{Subsuperscript}[c,\\text{lequ},1]')

clequ3 = Parameter(name = 'clequ3',
                   nature = 'internal',
                   type = 'complex',
                   value = 'clequ3Abs*cmath.exp(clequ3Ph*complex(0,1))',
                   texname = '\\text{Subsuperscript}[c,\\text{lequ},3]')

vevhat = Parameter(name = 'vevhat',
                   nature = 'internal',
                   type = 'real',
                   value = '1/(2**0.25*cmath.sqrt(Gf))',
                   texname = '\\hat{v}')

lam = Parameter(name = 'lam',
                nature = 'internal',
                type = 'real',
                value = '(Gf*MH**2)/cmath.sqrt(2)',
                texname = '\\text{lam}')

G = Parameter(name = 'G',
              nature = 'internal',
              type = 'real',
              value = '2*cmath.sqrt(aS)*cmath.sqrt(cmath.pi)',
              texname = 'G')

ee = Parameter(name = 'ee',
               nature = 'internal',
               type = 'real',
               value = '2*cmath.sqrt(aEW)*cmath.sqrt(cmath.pi)',
               texname = 'e')

CKM1x1 = Parameter(name = 'CKM1x1',
                   nature = 'internal',
                   type = 'complex',
                   value = '1 - CKMlambda**2/2.',
                   texname = '\\text{CKM1x1}')

CKM1x2 = Parameter(name = 'CKM1x2',
                   nature = 'internal',
                   type = 'complex',
                   value = 'CKMlambda',
                   texname = '\\text{CKM1x2}')

CKM1x3 = Parameter(name = 'CKM1x3',
                   nature = 'internal',
                   type = 'complex',
                   value = 'CKMA*CKMlambda**3*(CKMrho - CKMeta*complex(0,1))',
                   texname = '\\text{CKM1x3}')

CKM2x1 = Parameter(name = 'CKM2x1',
                   nature = 'internal',
                   type = 'complex',
                   value = '-CKMlambda',
                   texname = '\\text{CKM2x1}')

CKM2x2 = Parameter(name = 'CKM2x2',
                   nature = 'internal',
                   type = 'complex',
                   value = '1 - CKMlambda**2/2.',
                   texname = '\\text{CKM2x2}')

CKM2x3 = Parameter(name = 'CKM2x3',
                   nature = 'internal',
                   type = 'complex',
                   value = 'CKMA*CKMlambda**2',
                   texname = '\\text{CKM2x3}')

CKM3x1 = Parameter(name = 'CKM3x1',
                   nature = 'internal',
                   type = 'complex',
                   value = 'CKMA*CKMlambda**3*(1 - CKMrho - CKMeta*complex(0,1))',
                   texname = '\\text{CKM3x1}')

CKM3x2 = Parameter(name = 'CKM3x2',
                   nature = 'internal',
                   type = 'complex',
                   value = '-(CKMA*CKMlambda**2)',
                   texname = '\\text{CKM3x2}')

CKM3x3 = Parameter(name = 'CKM3x3',
                   nature = 'internal',
                   type = 'complex',
                   value = '1',
                   texname = '\\text{CKM3x3}')

sth2 = Parameter(name = 'sth2',
                 nature = 'internal',
                 type = 'real',
                 value = '0.5 - cmath.sqrt(1 - (2*aEW*cmath.pi*cmath.sqrt(2))/(Gf*MZ**2))/2.',
                 texname = '\\text{sth2}')

cth = Parameter(name = 'cth',
                nature = 'internal',
                type = 'real',
                value = 'cmath.sqrt(1 - sth2)',
                texname = 'c_{\\theta }')

gHgg = Parameter(name = 'gHgg',
                 nature = 'internal',
                 type = 'real',
                 value = '(G**2*Hloopcouplings*(0.3333333333333333 + (13*MH**6)/(50400.*MT**6) + MH**4/(504.*MT**4) + (7*MH**2)/(360.*MT**2)))/(16.*cmath.pi**2)',
                 texname = 'g_{\\text{HGG}}')

sth = Parameter(name = 'sth',
                nature = 'internal',
                type = 'real',
                value = 'cmath.sqrt(sth2)',
                texname = 's_{\\theta }')

dGf = Parameter(name = 'dGf',
                nature = 'internal',
                type = 'real',
                value = '(vevhat**2*(-(cll + cll1)/(2.*cmath.sqrt(2)) + cHl3*cmath.sqrt(2)))/LambdaSMEFT**2',
                texname = '\\text{dGf}')

dMH2 = Parameter(name = 'dMH2',
                 nature = 'internal',
                 type = 'real',
                 value = '((2*cHbox - cHD/2. - (3*cH)/(2.*lam))*MH**2*vevhat**2)/LambdaSMEFT**2',
                 texname = '\\text{dMH2}')

yb = Parameter(name = 'yb',
               nature = 'internal',
               type = 'real',
               value = '(ymb*cmath.sqrt(2))/vevhat',
               texname = '\\text{yb}')

yc = Parameter(name = 'yc',
               nature = 'internal',
               type = 'real',
               value = '(ymc*cmath.sqrt(2))/vevhat',
               texname = '\\text{yc}')

ydo = Parameter(name = 'ydo',
                nature = 'internal',
                type = 'real',
                value = '(ymdo*cmath.sqrt(2))/vevhat',
                texname = '\\text{ydo}')

ye = Parameter(name = 'ye',
               nature = 'internal',
               type = 'real',
               value = '(yme*cmath.sqrt(2))/vevhat',
               texname = '\\text{ye}')

ym = Parameter(name = 'ym',
               nature = 'internal',
               type = 'real',
               value = '(ymm*cmath.sqrt(2))/vevhat',
               texname = '\\text{ym}')

ys = Parameter(name = 'ys',
               nature = 'internal',
               type = 'real',
               value = '(yms*cmath.sqrt(2))/vevhat',
               texname = '\\text{ys}')

yt = Parameter(name = 'yt',
               nature = 'internal',
               type = 'real',
               value = '(ymt*cmath.sqrt(2))/vevhat',
               texname = '\\text{yt}')

ytau = Parameter(name = 'ytau',
                 nature = 'internal',
                 type = 'real',
                 value = '(ymtau*cmath.sqrt(2))/vevhat',
                 texname = '\\text{ytau}')

yup = Parameter(name = 'yup',
                nature = 'internal',
                type = 'real',
                value = '(ymup*cmath.sqrt(2))/vevhat',
                texname = '\\text{yup}')

barlam = Parameter(name = 'barlam',
                   nature = 'internal',
                   type = 'real',
                   value = 'lam*(1 - dMH2/MH**2 - dGf*cmath.sqrt(2))',
                   texname = '\\text{barlam}')

MW0 = Parameter(name = 'MW0',
                nature = 'internal',
                type = 'real',
                value = 'cth*MZ',
                texname = '\\text{MW0}')

vevT = Parameter(name = 'vevT',
                 nature = 'internal',
                 type = 'real',
                 value = 'vevhat*(1 + dGf/cmath.sqrt(2))',
                 texname = '\\text{vevT}')

dMZ2 = Parameter(name = 'dMZ2',
                 nature = 'internal',
                 type = 'real',
                 value = '(MZ**2*(cHD/2. + 2*cHWB*cth*sth)*vevhat**2)/LambdaSMEFT**2',
                 texname = '\\text{dMZ2}')

g1 = Parameter(name = 'g1',
               nature = 'internal',
               type = 'real',
               value = 'ee/cth',
               texname = 'g_1')

gw = Parameter(name = 'gw',
               nature = 'internal',
               type = 'real',
               value = 'ee/sth',
               texname = 'g_w')

gHaa = Parameter(name = 'gHaa',
                 nature = 'internal',
                 type = 'real',
                 value = '(ee**2*Hloopcouplings*(-1.75 + (4*(0.3333333333333333 + (13*MH**6)/(50400.*MT**6) + MH**4/(504.*MT**4) + (7*MH**2)/(360.*MT**2)))/3. - (29*MH**6)/(16800.*MW0**6) - (19*MH**4)/(1680.*MW0**4) - (11*MH**2)/(120.*MW0**2)))/(8.*cmath.pi**2)',
                 texname = 'g_{\\text{H$\\gamma \\gamma $}}')

gHza = Parameter(name = 'gHza',
                 nature = 'internal',
                 type = 'real',
                 value = '(ee**2*Hloopcouplings*(((0.4583333333333333 + (29*MH**6)/(100800.*MW0**6) + (19*MH**4)/(10080.*MW0**4) + (11*MH**2)/(720.*MW0**2) + (MH**4*MZ**2)/(2100.*MW0**6) + (MH**2*MZ**2)/(280.*MW0**4) + (7*MZ**2)/(180.*MW0**2) + (67*MH**2*MZ**4)/(100800.*MW0**6) + (53*MZ**4)/(10080.*MW0**4) + (43*MZ**6)/(50400.*MW0**6) - (31*cth**2)/(24.*sth**2) - (29*cth**2*MH**6)/(20160.*MW0**6*sth**2) - (19*cth**2*MH**4)/(2016.*MW0**4*sth**2) - (11*cth**2*MH**2)/(144.*MW0**2*sth**2) - (cth**2*MH**4*MZ**2)/(560.*MW0**6*sth**2) - (31*cth**2*MH**2*MZ**2)/(2520.*MW0**4*sth**2) - (cth**2*MZ**2)/(9.*MW0**2*sth**2) - (43*cth**2*MH**2*MZ**4)/(20160.*MW0**6*sth**2) - (17*cth**2*MZ**4)/(1120.*MW0**4*sth**2) - (5*cth**2*MZ**6)/(2016.*MW0**6*sth**2))*sth)/cth + ((0.3333333333333333 + (13*MH**6)/(50400.*MT**6) + MH**4/(504.*MT**4) + (7*MH**2)/(360.*MT**2) + (MH**4*MZ**2)/(2400.*MT**6) + (MH**2*MZ**2)/(315.*MT**4) + (11*MZ**2)/(360.*MT**2) + (29*MH**2*MZ**4)/(50400.*MT**6) + (11*MZ**4)/(2520.*MT**4) + (37*MZ**6)/(50400.*MT**6))*(0.5 - (4*sth**2)/3.))/(cth*sth)))/(4.*cmath.pi**2)',
                 texname = 'g_{\\text{HZ$\\gamma $}}')

vev = Parameter(name = 'vev',
                nature = 'internal',
                type = 'real',
                value = '(1 - (3*cH*vevhat**2)/(8.*lam*LambdaSMEFT**2))*vevT',
                texname = '\\text{vev}')

dg1 = Parameter(name = 'dg1',
                nature = 'internal',
                type = 'real',
                value = '(sth**2*(dMZ2/MZ**2 + (2*cHWB*cth**3*vevhat**2)/(LambdaSMEFT**2*sth) + dGf*cmath.sqrt(2)))/(2.*(1 - 2*sth**2))',
                texname = '\\text{dg1}')

dgw = Parameter(name = 'dgw',
                nature = 'internal',
                type = 'real',
                value = '-(cth**2*(dMZ2/MZ**2 + (2*cHWB*sth**3*vevhat**2)/(cth*LambdaSMEFT**2) + dGf*cmath.sqrt(2)))/(2.*(1 - 2*sth**2))',
                texname = '\\text{dgw}')

g1sh = Parameter(name = 'g1sh',
                 nature = 'internal',
                 type = 'real',
                 value = '(ee*(1 + dg1 - (cHB*vevhat**2)/LambdaSMEFT**2))/cth',
                 texname = 'g_1')

gwsh = Parameter(name = 'gwsh',
                 nature = 'internal',
                 type = 'real',
                 value = '(ee*(1 + dgw - (cHW*vevhat**2)/LambdaSMEFT**2))/sth',
                 texname = 'g_w')

dMW = Parameter(name = 'dMW',
                nature = 'internal',
                type = 'real',
                value = '(MW0*(2*dgw + dGf*cmath.sqrt(2)))/2.',
                texname = '\\text{dMW}')

dsth2 = Parameter(name = 'dsth2',
                  nature = 'internal',
                  type = 'real',
                  value = '2*cth**2*(dg1 - dgw)*sth**2 + (cHWB*cth*sth*(1 - 2*sth**2)*vevhat**2)/LambdaSMEFT**2',
                  texname = '\\text{dsth2}')

MW = Parameter(name = 'MW',
               nature = 'internal',
               type = 'real',
               value = 'dMW + MW0',
               texname = 'M_W')

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

This is actually a question (which is more a request actually) for the FeynRules author.
I can only tell you one potential reason of this.

The UFO model wants to have as input a card following the LesHouchesAccord, (like the SLHA ones for examples).
But in that community-wide agreement, all input parameter are supposed to be a real number. Therefore if you want complex parameter it is usual to have two entry in that file (the real part and the imaginary part).

The next LesHouches is actually in two years, so this is probably not going to change soon.

Cheers,

Olivier

Revision history for this message
Ilaria Brivio (ilariab) said :
#4

Hi Olivier,

thanks. It can be as you say, but the issue is still weird, because this complex parameter is an internal one. in the param card one fixes absolute value and phase, which are real numbers.

The trick works for the CKM and in the other version of SMEFTsim, so I'm just really confused about why it's not interpreted correctly in this case.

Is is possible to pass this on to a Feynrules author?

Thanks again!
Ilaria

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

Hi,

Feynrules decides to not use the launchpad platform for their support to the user.
In the past they were using a ticket system on their website https://feynrules.irmp.ucl.ac.be but it does not seem active.
So my advice would be that you just send an email to Claude Duhr.

Cheers,

Olivier

> On 18 Jan 2018, at 13:27, Ilaria Brivio <email address hidden> wrote:
>
> Question #663218 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/663218
>
> Status: Answered => Open
>
> Ilaria Brivio is still having a problem:
> Hi Olivier,
>
> thanks. It can be as you say, but the issue is still weird, because this
> complex parameter is an internal one. in the param card one fixes
> absolute value and phase, which are real numbers.
>
> The trick works for the CKM and in the other version of SMEFTsim, so I'm
> just really confused about why it's not interpreted correctly in this
> case.
>
> Is is possible to pass this on to a Feynrules author?
>
> Thanks again!
> Ilaria
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Ilaria Brivio (ilariab) said :
#6

Oh, sorry, I thought you could pass it on (I had the impression I had read some FR answers here but maybe that's inaccurate, or just old messages). My bad.
I'll send an email then.

Thanks anyways!
Cheers,
Ilaria