Version 12 (modified by 4 years ago) ( diff ) | ,
---|
Standard Model Effective Theory at One-Loop in QCD
Céline Degrande, Gauthier Durieux, Fabio Maltoni, Ken Mimasu, Eleni Vryonidou & Cen Zhang
The implementation is based on the Warsaw basis of dimension-six SMEFT operators, after canonical normalization. Electroweak input parameters are taken to be GF, MZ, MW. The CKM matrix is approximated as a unit matrix, and a U(2)q x U(2)u x U(3)d x (U(1)l x U(1)e)3 flavor symmetry is enforced. It forbids all fermion masses and Yukawa couplings except that only of the top quark. The model therefore implements the five-flavor scheme for PDFs.
A new coupling order, NP=2
, is assigned to SMEFT interactions.
The cutoff scale Lambda
takes a default value of 1 TeV-2 and can be modified along with the Wilson coefficients in the param_card
.
Operators definitions, normalisations and coefficient names in the UFO model are specified in definitions.pdf.
The notations and normalizations of top-quark operator coefficients comply with the LHC TOP WG standards of 1802.07237.
Note however that the flavor symmetry enforced here is slightly more restrictive than the baseline assumption there (see the dim6top page for more information).
This model has been validated at tree level against the dim6top
implementation (see 1906.12310 and the comparison details).
Current implementation
UFO model: SMEFTatNLO_v1.0.tar.gz
The current implementation imposes CP conservation.
In the quark sector, it focuses primarily on top-quark interactions.
The light-quark current operator, qqHDH, uuHDH, ddHDH, with coefficients cpq3i
, cpqMi
, cpu
, cpd
are however included.
The triple-gluon operator, with coefficient cG
, is currently not available (see the loop-capable GGG implementation).
Vertices including more than four scalars or four leptons are not included.
Scalar and tensor QQll
operators, with coefficients ctlS3
, ctlT3
, and cblS3
, break our flavor symmetry assumption and are not available for one-loop computations.
Top-quark flavor-changing interactions, not compatible with the imposed flavor symmetry, are not included (see the loop-capable TopFCNC implementation).
Unlike prescribed by the LHC TOP WG, the top quark chromomagnetic-dipole operator coefficient ctG
is normalized with a factor of the strong coupling, gS.
This normalization factor temporarily ensures compatibility with the 2.X.X series of MadGraph5_aMC@NLO but may be dropped in the future.
As with every other appearance of this coupling in MadGraph5_aMC@NLO, its value is renormalisation-group evolved to the QCD renormalization scale (set in the run_card).
Counterterms required for one-loop computations are currently included up to five points.
The unitary gauge (default) is recommended when computing anomalous quark-loop amplitudes like ggZ
, gggZ
, ggZH
and ggff
.
MadGraph5_aMC@NLO does not evolve operator coefficients which are therefore kept at fixed scale mueft
distinguished from the QCD renormalization scale MUR
.
We recommend to use fixed renormalization and factorization scales (in the run_card
), and to set mueft
equal to those (in the param_card
).
The 3.0.3-neworders
development branch (tarball) of MG is required for NLO predictions involving four-fermion operators and (in general) H2G2 with coefficient cpG
not normalized with any power of gS.
It also allows for a better control over coupling orders and, in particular, for the separate computation of linear and quadratic EFT contributions at NLO, in fixed order mode.
A branch allowing for the separate computation of different orders in event-generation mode (with matching to parton shower) is being validated.
The 2.X.X series of MadGraph5_aMC@NLO can handle bosonic and two-fermion operators at one-loop.
Version updates
- 2018/12/20 - v0.1: First version upload, 4F and cG operators at LO pending validation; a few minor convention tweaks required to match
dim6top
exactly.decays.py
missing. - 2019/04/03 - v0.1: Added definitions.pdf document and uploaded a new version with a fix for restrict_default.dat
- 2019/08/12 - v0.1: Uploaded a new version matching
dim6top
operator conventions, also some bugfixes and gs normalisation forOtG
- 2020/08/24 - v1.0: Official release including notably four-quark operators at NLO.
Usage notes
Restriction cards
Because of the mixture of LO/NLO compatible operators included in the model, restriction cards must be used to access the SMEFT interactions.
Default loading of the model
> import model SMEFTatNLO_vX.X
will load the pure SM without any effective operators.
The LO
restriction card should be used when importing the model for LO generation:
> import model SMEFTatNLO_vX.X-LO
For NLO QCD generation, the NLO
restriction card should be used when importing the model:
> import model SMEFTatNLO_vX.X-NLO
This invokes a restricted set of operators for which the required counterterms are implemented.
Coupling orders
We recommend specifying the full QCD
, QED
and NP
orders for process generation.
For example:
> generate p p > t t~ QCD=2 QED=0 NP=2 [QCD]
generates top-quark pair production at NLO QCD, including the QCD-induced SM and the SMEFT contributions.
Excluding operators
We recommend avoiding setting values of Wilson coefficients to 0 when computing at NLO using MadGraph5_aMC@NLO.
Operators should either be removed explicitly with restriction cards or set to a very small non-zero value, e.g., 1e-5
Plugin for b-quark Yukawa coupling and operator (ymb
and cbp
)
A plugin-like modification to the model including the bbh (SM+SMEFT), bbhh and bbhhh interactions has been implemented to account for the Higgs coupling to bottom quarks.
It can only be used at LO.
A configuration.py
file is included in the UFO model with a bottomYukawa
flag set to False
by default.
Setting it to True
restores the SM & SMEFT bottom Yukawa parameters (ymb
and cbp
), the bbh(h)(h) vertices, and corresponding couplings.
The bottom mass parameters, MB
, is not restored which has a percent effect on the h > b b~
partial width.
The corresponding Goldstone boson interactions are not included, such that the extended model can only be used in unitary gauge (default).
Generation recipes for validated processes
Among many others, the following processes are supported at the one-loop level.
Gauge invariance (see help check
in MadGraph5_aMC@NLO) and pole cancellation have been checked explicitly for those.
For complicated processes and in case of doubts, please contact the authors.
Widths should be set to zero to ensure gauge invariance.
QCD
> p p > j j QED=0 QCD=2 NP=2 [QCD]
Drell Yan
> p p > mu+ mu- QCD=0 QED=2 NP=2 [QCD] > p p > mu+ vm QCD=0 QED=2 NP=2 [QCD] > p p > W+ j $$ t QCD=1 QED=1 NP=2 [QCD] > p p > W- j $$ t~ QCD=1 QED=1 NP=2 [QCD] > p p > Z j QCD=1 QED=1 NP=2 [QCD]
Multi-boson production
quark-initiated
> p p > W+ W- QED=2 QCD=0 NP=2 [QCD] > p p > W+ Z QED=2 QCD=0 NP=2 [QCD] > p p > Z Z QED=2 QCD=0 NP=2 [QCD]
loop-induced
> g g > W+ W- QED=2 QCD=2 NP=2 [QCD] > g g > Z Z QED=2 QCD=2 NP=2 [QCD] > g g > W+ W- Z QED=3 QCD=2 NP=2 [QCD] > g g > Z Z Z QED=3 QCD=2 NP=2 [QCD]
Higgs production
loop-induced
> g g > H QED=1 QCD=2 NP=2 [QCD] > g g > H H QED=2 QCD=2 NP=2 [QCD] > g g > H H H QED=3 QCD=2 NP=2 [QCD] > g g > H j QED=1 QCD=3 NP=2 [QCD]
Top quark production
> e+ e- > t t~ QED=2 QCD=0 NP=2 [QCD] > p p > t t~ QED=0 QCD=2 NP=2 [QCD] > p p > t t~ h QED=1 QCD=2 NP=2 [QCD] > p p > t t~ Z QED=1 QCD=2 NP=2 [QCD] > p p > t t~ W+ QED=1 QCD=2 NP=2 [QCD] > p p > t W- $$ t~ QED=1 QCD=1 NP=2 [QCD] > p p > t W- j $$ t~ QED=1 QCD=2 NP=2 [QCD] > p p > t j $$ W- QED=2 QCD=0 NP=2 [QCD] > p p > t h j $$ W- QED=3 QCD=0 NP=2 [QCD] > p p > t Z j $$ W- QED=3 QCD=0 NP=2 [QCD] > p p > t a j $$ W- QED=3 QCD=0 NP=2 [QCD]
When generating one of the last four processes (tj
,thj
,tZj
,taj
) with the cQq83
operator coefficient, all loops including a gluon have to be allowed.
This can be achieved with the following modification of MadGraph5_aMC@NLO:
=== modified file 'madgraph/loop/loop_diagram_generation.py' --- madgraph/loop/loop_diagram_generation.py 2020-03-11 09:28:14 +0000 +++ madgraph/loop/loop_diagram_generation.py 2020-04-03 21:08:18 +0000 @@ -384,7 +384,7 @@ # By default the user filter does nothing if filter is not set, # if you want to turn it on and edit it by hand, then set the # variable edit_filter_manually to True - edit_filter_manually = False + edit_filter_manually = True if not edit_filter_manually and filter in [None,'None']: return if isinstance(filter,str) and filter.lower() == 'true': @@ -415,6 +415,10 @@ raise InvalidCmd("The user-defined filter '%s' did not"%filter+ " returned the following error:\n > %s"%str(e)) + # requires a gluon to run in all loops + if 21 not in diag.get_loop_lines_pdgs(): + valid_diag = False + # if any([abs(pdg) not in range(1,7) for pdg in diag.get_loop_lines_pdgs()]): # valid_diag = False @@ -538,7 +542,7 @@ if valid_diag: newloopselection.append(diag) - self['loop_diagrams']=newloopselection + #self['loop_diagrams']=newloopselection # To monitor what are the diagrams filtered, simply comment the line # directly above and uncomment the two directly below. # self['loop_diagrams'] = base_objects.DiagramList(
Analytic validation
The following loop computations of amplitudes relevant for several processes have been cross-checked analytically:
- ttbar: tt, gg, ggg, gtt, ggtt
- single top/decay: tbW, 4f
- ttV: ttV, ggV, gggV, gttV
- ttH: ggh, gggh, htt, ghtt
Attachments (6)
-
SMEFTatNLO_v1.0.tar.gz
(167.9 KB
) - added by 4 years ago.
v1.0 uploaded 2020/09/09
-
SMEFTatNLO_v1.0.1.tar.gz
(168.2 KB
) - added by 4 years ago.
python3 compatible, BR information in restrictions for easier Madspin use, no4q restriction for MG v2
-
SMEFTatNLO_v1.0.2.tar.gz
(168.4 KB
) - added by 4 years ago.
distinguished slightly MZ, MU_R and mueft in restriction cards
-
SMEFTatNLO_v1.0.3.tar.gz
(168.1 KB
) - added by 4 years ago.
Fixed typo in b-yukawa add-on
-
tree-loop_instructions.txt
(7.9 KB
) - added by 3 years ago.
Instructions for computing SMEFT dependence from mixed tree/loop processes
-
definitions.pdf
(208.6 KB
) - added by 2 months ago.
Notation, conventions, field redefinitions & input scheme
Download all attachments as: .zip