complex parameters not exported to UFO in feynrules
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://
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_
OuW:= Block[{
feynmangaugerules = If[Not[
ExpandIndices[
]//OptimizeIndex;
[from d6_parameters_
cuWPh == {ParameterType -> External, Value -> 0, ComplexParameter -> False},
cuWAbs == {
ParameterType -> External,
Value -> 1,
Interaction
},
cuW == {
ParameterType -> Internal,
Value -> cuWAbs Exp[I cuWPh],
ComplexPara
Interaction
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
- Assignee:
- No assignee Edit question
- Solved by:
- Ilaria Brivio
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#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',
GC_2 = Coupling(name = 'GC_2',
GC_3 = Coupling(name = 'GC_3',
GC_4 = Coupling(name = 'GC_4',
GC_5 = Coupling(name = 'GC_5',
GC_6 = Coupling(name = 'GC_6',
GC_7 = Coupling(name = 'GC_7',
GC_8 = Coupling(name = 'GC_8',
GC_9 = Coupling(name = 'GC_9',
GC_10 = Coupling(name = 'GC_10',
GC_11 = Coupling(name = 'GC_11',
GC_12 = Coupling(name = 'GC_12',
GC_13 = Coupling(name = 'GC_13',
GC_14 = Coupling(name = 'GC_14',
GC_15 = Coupling(name = 'GC_15',
GC_16 = Coupling(name = 'GC_16',
GC_17 = Coupling(name = 'GC_17',
GC_18 = Coupling(name = 'GC_18',
GC_19 = Coupling(name = 'GC_19',
GC_20 = Coupling(name = 'GC_20',
GC_21 = Coupling(name = 'GC_21',
GC_22 = Coupling(name = 'GC_22',
GC_23 = Coupling(name = 'GC_23',
GC_24 = Coupling(name = 'GC_24',
GC_25 = Coupling(name = 'GC_25',
GC_26 = Coupling(name = 'GC_26',
GC_27 = Coupling(name = 'GC_27',
GC_28 = Coupling(name = 'GC_28',
GC_29 = Coupling(name = 'GC_29',
GC_30 = Coupling(name = 'GC_30',
GC_31 = Coupling(name = 'GC_31',
GC_32 = Coupling(name = 'GC_32',
GC_33 = Coupling(name = 'GC_33',
GC_34 = Coupling(name = 'GC_34',
GC_35 = Coupling(name = 'GC_35',
GC_36 = Coupling(name = 'GC_36',
GC_37 = Coupling(name = 'GC_37',
GC_38 = Coupling(name = 'GC_38',
GC_39 = Coupling(name = 'GC_39',
GC_40 = Coupling(name = 'GC_40',
GC_41 = Coupling(name = 'GC_41',
GC_42 = Coupling(name = 'GC_42',
GC_43 = Coupling(name = 'GC_43',
GC_44 = Coupling(name = 'GC_44',
GC_45 = Coupling(name = 'GC_45',
GC_46 = Coupling(name = 'GC_46',
GC_47 = Coupling(name = 'GC_47',
Revision history for this message
|
#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',
# User-defined parameters.
cabi = Parameter(name = 'cabi',
CKMlambda = Parameter(name = 'CKMlambda',
CKMA = Parameter(name = 'CKMA',
CKMrho = Parameter(name = 'CKMrho',
CKMeta = Parameter(name = 'CKMeta',
aEW = Parameter(name = 'aEW',
Gf = Parameter(name = 'Gf',
type = 'real',
aS = Parameter(name = 'aS',
type = 'real',
ymdo = Parameter(name = 'ymdo',
ymup = Parameter(name = 'ymup',
yms = Parameter(name = 'yms',
ymc = Parameter(name = 'ymc',
ymb = Parameter(name = 'ymb',
ymt = Parameter(name = 'ymt',
yme = Parameter(name = 'yme',
ymm = Parameter(name = 'ymm',
ymtau = Parameter(name = 'ymtau',
ceWPh = Parameter(name = 'ceWPh',
ceBPh = Parameter(name = 'ceBPh',
cuGPh = Parameter(name = 'cuGPh',
cuWPh = Parameter(name = 'cuWPh',
cuBPh = Parameter(name = 'cuBPh',
cdGPh = Parameter(name = 'cdGPh',
cdWPh = Parameter(name = 'cdWPh',
cdBPh = Parameter(name = 'cdBPh',
cHudPh = Parameter(name = 'cHudPh',
ceHPh = Parameter(name = 'ceHPh',
cuHPh = Parameter(name = 'cuHPh',
cdHPh = Parameter(name = 'cdHPh',
cledqPh = Parameter(name = 'cledqPh',
cquqd1Ph = Parameter(name = 'cquqd1Ph',
cquqd8Ph = Parameter(name = 'cquqd8Ph',
clequ1Ph = Parameter(name = 'clequ1Ph',
clequ3Ph = Parameter(name = 'clequ3Ph',
LambdaSMEFT = Parameter(name = 'LambdaSMEFT',
cG = Parameter(name = 'cG',
type = 'real',
cGtil = Parameter(name = 'cGtil',
cW = Parameter(name = 'cW',
type = 'real',
cWtil = Parameter(name = 'cWtil',
cH = Parameter(name = 'cH',
type = 'real',
cHbox = Parameter(name = 'cHbox',
cHD = Parameter(name = 'cHD',
cHG = Parameter(name = 'cHG',
cHGtil = Parameter(name = 'cHGtil',
cHW = Parameter(name = 'cHW',
cHWtil = Parameter(name = 'cHWtil',
cHB = Parameter(name = 'cHB',
cHBtil = Parameter(name = 'cHBtil',
cHWB = Parameter(name = 'cHWB',
cHWBtil = Parameter(name = 'cHWBtil',
ceHAbs = Parameter(name = 'ceHAbs',
cuHAbs = Parameter(name = 'cuHAbs',
cdHAbs = Parameter(name = 'cdHAbs',
ceWAbs = Parameter(name = 'ceWAbs',
ceBAbs = Parameter(name = 'ceBAbs',
cuGAbs = Parameter(name = 'cuGAbs',
cuWAbs = Parameter(name = 'cuWAbs',
cuBAbs = Parameter(name = 'cuBAbs',
cdGAbs = Parameter(name = 'cdGAbs',
cdWAbs = Parameter(name = 'cdWAbs',
cdBAbs = Parameter(name = 'cdBAbs',
cHl1 = Parameter(name = 'cHl1',
cHl3 = Parameter(name = 'cHl3',
cHe = Parameter(name = 'cHe',
cHq1 = Parameter(name = 'cHq1',
cHq3 = Parameter(name = 'cHq3',
cHu = Parameter(name = 'cHu',
cHd = Parameter(name = 'cHd',
cHudAbs = Parameter(name = 'cHudAbs',
cll = Parameter(name = 'cll',
cll1 = Parameter(name = 'cll1',
cqq1 = Parameter(name = 'cqq1',
cqq11 = Parameter(name = 'cqq11',
cqq3 = Parameter(name = 'cqq3',
cqq31 = Parameter(name = 'cqq31',
clq1 = Parameter(name = 'clq1',
clq3 = Parameter(name = 'clq3',
cee = Parameter(name = 'cee',
cuu = Parameter(name = 'cuu',
cuu1 = Parameter(name = 'cuu1',
cdd = Parameter(name = 'cdd',
cdd1 = Parameter(name = 'cdd1',
ceu = Parameter(name = 'ceu',
ced = Parameter(name = 'ced',
cud1 = Parameter(name = 'cud1',
cud8 = Parameter(name = 'cud8',
cle = Parameter(name = 'cle',
clu = Parameter(name = 'clu',
cld = Parameter(name = 'cld',
cqe = Parameter(name = 'cqe',
cqu1 = Parameter(name = 'cqu1',
cqu8 = Parameter(name = 'cqu8',
cqd1 = Parameter(name = 'cqd1',
cqd8 = Parameter(name = 'cqd8',
cledqAbs = Parameter(name = 'cledqAbs',
cquqd1Abs = Parameter(name = 'cquqd1Abs',
cquqd8Abs = Parameter(name = 'cquqd8Abs',
clequ1Abs = Parameter(name = 'clequ1Abs',
clequ3Abs = Parameter(name = 'clequ3Abs',
Hloopcouplings = Parameter(name = 'Hloopcouplings',
MZ = Parameter(name = 'MZ',
type = 'real',
Me = Parameter(name = 'Me',
type = 'real',
MMU = Parameter(name = 'MMU',
MTA = Parameter(name = 'MTA',
MU = Parameter(name = 'MU',
type = 'real',
MC = Parameter(name = 'MC',
type = 'real',
MT = Parameter(name = 'MT',
type = 'real',
MD = Parameter(name = 'MD',
type = 'real',
MS = Parameter(name = 'MS',
type = 'real',
MB = Parameter(name = 'MB',
type = 'real',
MH = Parameter(name = 'MH',
type = 'real',
WZ = Parameter(name = 'WZ',
type = 'real',
WW = Parameter(name = 'WW',
type = 'real',
WT = Parameter(name = 'WT',
type = 'real',
WH = Parameter(name = 'WH',
type = 'real',
ceH = Parameter(name = 'ceH',
cuH = Parameter(name = 'cuH',
cdH = Parameter(name = 'cdH',
ceW = Parameter(name = 'ceW',
ceB = Parameter(name = 'ceB',
cuG = Parameter(name = 'cuG',
cuW = Parameter(name = 'cuW',
cuB = Parameter(name = 'cuB',
cdG = Parameter(name = 'cdG',
cdW = Parameter(name = 'cdW',
cdB = Parameter(name = 'cdB',
cHud = Parameter(name = 'cHud',
cledq = Parameter(name = 'cledq',
cquqd1 = Parameter(name = 'cquqd1',
cquqd8 = Parameter(name = 'cquqd8',
clequ1 = Parameter(name = 'clequ1',
clequ3 = Parameter(name = 'clequ3',
vevhat = Parameter(name = 'vevhat',
lam = Parameter(name = 'lam',
G = Parameter(name = 'G',
type = 'real',
value = '2*cmath.
ee = Parameter(name = 'ee',
type = 'real',
CKM1x1 = Parameter(name = 'CKM1x1',
CKM1x2 = Parameter(name = 'CKM1x2',
CKM1x3 = Parameter(name = 'CKM1x3',
CKM2x1 = Parameter(name = 'CKM2x1',
CKM2x2 = Parameter(name = 'CKM2x2',
CKM2x3 = Parameter(name = 'CKM2x3',
CKM3x1 = Parameter(name = 'CKM3x1',
CKM3x2 = Parameter(name = 'CKM3x2',
CKM3x3 = Parameter(name = 'CKM3x3',
sth2 = Parameter(name = 'sth2',
cth = Parameter(name = 'cth',
gHgg = Parameter(name = 'gHgg',
sth = Parameter(name = 'sth',
dGf = Parameter(name = 'dGf',
dMH2 = Parameter(name = 'dMH2',
yb = Parameter(name = 'yb',
type = 'real',
yc = Parameter(name = 'yc',
type = 'real',
ydo = Parameter(name = 'ydo',
ye = Parameter(name = 'ye',
type = 'real',
ym = Parameter(name = 'ym',
type = 'real',
ys = Parameter(name = 'ys',
type = 'real',
yt = Parameter(name = 'yt',
type = 'real',
ytau = Parameter(name = 'ytau',
yup = Parameter(name = 'yup',
barlam = Parameter(name = 'barlam',
MW0 = Parameter(name = 'MW0',
vevT = Parameter(name = 'vevT',
dMZ2 = Parameter(name = 'dMZ2',
g1 = Parameter(name = 'g1',
type = 'real',
gw = Parameter(name = 'gw',
type = 'real',
gHaa = Parameter(name = 'gHaa',
gHza = Parameter(name = 'gHza',
vev = Parameter(name = 'vev',
dg1 = Parameter(name = 'dg1',
dgw = Parameter(name = 'dgw',
g1sh = Parameter(name = 'g1sh',
gwsh = Parameter(name = 'gwsh',
dMW = Parameter(name = 'dMW',
dsth2 = Parameter(name = 'dsth2',
MW = Parameter(name = 'MW',
type = 'real',
Revision history for this message
|
#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
|
#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
|
#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:/
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:/
>
> 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
|
#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