Commit 74fe337b authored by Despina Sampsoniodou's avatar Despina Sampsoniodou

model, cards, README

parent f1e66553
Pipeline #27 failed with stages
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################
###################################
## INFORMATION FOR FRBLOCK
###################################
Block frblock
18 1.000000e+03 # LambdaSMEFT
19 1.000000e+00 # cG
20 1.000000e+00 # cGtil
21 1.000000e+00 # cW
22 1.000000e+00 # cWtil
24 1.000000e+00 # cHbox
25 1.000000e+00 # cHDD
26 1.000000e+00 # cHG
27 1.000000e+00 # cHGtil
28 1.000000e+00 # cHW
30 1.000000e+00 # cHB
31 1.000000e+00 # cHBtil
32 1.000000e+00 # cHWB
33 1.000000e+00 # cHWBtil
45 1.000000e+00 # cHl1
46 1.000000e+00 # cHl3
47 1.000000e+00 # cHe
48 1.000000e+00 # cHq1
49 1.000000e+00 # cHq3
50 1.000000e+00 # cHu
51 1.000000e+00 # cHd
54 1.000000e+00 # cll1
###################################
## INFORMATION FOR MASS
###################################
Block mass
5 4.180000e+00 # MB
6 1.732000e+02 # MT
23 9.118760e+01 # MZ
25 1.250900e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
4 0.000000 # c : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
15 0.000000 # ta- : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
24 73.187697 # w+ : dMW + MW0
###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
1 7.815553e-03 # aEW
2 1.166379e-05 # Gf
3 1.181000e-01 # aS
###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
5 4.180000e+00 # ymb
6 1.732000e+02 # ymt
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.085000e+00 # WW
DECAY 25 4.070000e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model SMEFTsim_A_U35_alphaScheme_UFO-massless
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > e+ e- mu- vm~ j j QCD=0 QED=6 NP=1
output SMEFTSim_VBS_total2
This diff is collapsed.
<<<<<<< HEAD
# SMEFTSim_Instructions
In order to test the effect of the Dimension 6 operators, we can utilize the SMEFTSim package(https://feynrules.irmp.ucl.ac.be/wiki/SMEFT) using your favourite MC generator. Below you will find instructions for Madgraph standalone:
......@@ -21,12 +20,3 @@ The output folder with the relevant Cards is created.
6.Inside the output folder,
> ./bin/generate_events
=======
# Dimension-6 instructions
>>>>>>> dad4f289907eebb90710ec3e38f960dee80350f5
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.2.0 for Linux x86 (64-bit) (September 11, 2017)
# Date: Fri 9 Mar 2018 12:55:25
from object_library import all_couplings, Coupling
from function_library import complexconjugate, re, im, csc, sec, acsc, asec, cot
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.2.0 for Linux x86 (64-bit) (September 11, 2017)
# Date: Fri 9 Mar 2018 12:29:16
#
# This is the logfile for the model SMEFTsim_A
# Authors: I. Brivio, Y. Jiang, M. Trott
# Model version: 2.1
# Checking the Quantum numbers
* Electric charge defined.
# Checking the Lagrangians
* All Lagrangians are ok.
#
# Particle definitions
#
* No particles removed. All particles correspond to GenInt setup.
# Automatically assigned PDG numbers
* Assigned PDG number 9000001 to particle ghA
* Assigned PDG number 9000002 to particle ghZ
* Assigned PDG number 9000003 to particle ghWp
* Assigned PDG number 9000004 to particle ghWm
# Compulsory PDG codes:
* Class SM leptons complete.
* Class SM neutrinos complete.
* Class SM quarks complete.
* Class SM gauge bosons complete.
#
# Parameter definitions
#
* All parameters are ok.
# Vertices
* Calling FeynmanRules for 1 Lagrangians.
* Number of classes vertices: 82
* Number of flavored vertices: 396
* Saved vertices in InterfaceRun[ 1 ].
* Checked QNumber conservation.
- Quantum number GhostNumber conserved in all vertices.
- Quantum number LeptonNumber conserved in all vertices.
- Quantum number Q conserved in all vertices.
- Quantum number Y conserved in all vertices.
* particles.py written.
* parameters.py written.
#
# Vertex definitions
#
* 467 vertices written.
* vertices.py written.
#
# Lorentz structure definitions
#
* 85 lorentz structures written.
* lorentz.py written.
#
# Coupling definitions
#
>>>> Non positive interaction order QED found for coupling GC_311.
* 1160 couplings written.
* couplings.py written.
#
# Coupling order definitions
#
* 0 couplings orders written.
* coupling_orders.py written.
#
# Decay definitions
#
* 12 decays written.
* decay.py not written
#
# CTCoupling definitions
#
* 0 CTcouplings written.
* CT_couplings.py written.
import particles
import couplings
import lorentz
import parameters
import vertices
import coupling_orders
import write_param_card
import propagators
import function_library
all_particles = particles.all_particles
all_vertices = vertices.all_vertices
all_couplings = couplings.all_couplings
all_lorentz = lorentz.all_lorentz
all_parameters = parameters.all_parameters
all_orders = coupling_orders.all_orders
all_functions = function_library.all_functions
all_propagators = propagators.all_propagators
try:
import decays
except ImportError:
pass
else:
all_decays = decays.all_decays
try:
import form_factors
except ImportError:
pass
else:
all_form_factors = form_factors.all_form_factors
try:
import CT_vertices
except ImportError:
pass
else:
all_CTvertices = CT_vertices.all_CTvertices
gauge = [0]
__author__ = "I. Brivio, Y. Jiang, M. Trott"
__date__ = "09. 03. 2018"
__version__= "2.1"
# This file was automatically created by FeynRules 2.3.29
# Mathematica version: 11.2.0 for Linux x86 (64-bit) (September 11, 2017)
# Date: Fri 9 Mar 2018 12:55:23
from object_library import all_orders, CouplingOrder
QCD = CouplingOrder(name = 'QCD',
expansion_order = 99,
hierarchy = 1)
QED = CouplingOrder(name = 'QED',
expansion_order = 99,
hierarchy = 2)
NP = CouplingOrder(name = 'NP',
expansion_order = 99,
hierarchy = 1)
SMHLOOP = CouplingOrder(name = 'SMHLOOP',
expansion_order = 99,
hierarchy = 0)
This diff is collapsed.
This diff is collapsed.
# This file is part of the UFO.
#
# This file contains definitions for functions that
# are extensions of the cmath library, and correspond
# either to functions that are in cmath, but inconvenient
# to access from there (e.g. z.conjugate()),
# or functions that are simply not defined.
#
#
__date__ = "22 July 2010"
__author__ = "claude.duhr@durham.ac.uk"
import cmath
from object_library import all_functions, Function
#
# shortcuts for functions from cmath
#
complexconjugate = Function(name = 'complexconjugate',
arguments = ('z',),
expression = 'z.conjugate()')
re = Function(name = 're',
arguments = ('z',),
expression = 'z.real')
im = Function(name = 'im',
arguments = ('z',),
expression = 'z.imag')
# New functions (trigonometric)
sec = Function(name = 'sec',
arguments = ('z',),
expression = '1./cmath.cos(z.real)')
asec = Function(name = 'asec',
arguments = ('z',),
expression = 'cmath.acos(1./(z.real))')
csc = Function(name = 'csc',
arguments = ('z',),
expression = '1./cmath.sin(z.real)')
acsc = Function(name = 'acsc',
arguments = ('z',),
expression = 'cmath.asin(1./(z.real))')
cot = Function(name = 'cot',
arguments = ('z',),
expression = '1./cmath.tan(z.real)')
# Heaviside theta function
theta_function = Function(name = 'theta_function',
arguments = ('x','y','z'),
expression = 'y if x else z')
# Auxiliary functions for NLO
cond = Function(name = 'cond',
arguments = ('condition','ExprTrue','ExprFalse'),
expression = '(ExprTrue if condition==0.0 else ExprFalse)')
reglog = Function(name = 'reglog',
arguments = ('z'),
expression = '(0.0 if z==0.0 else cmath.log(z.real))')
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
# This file was automatically created by FeynRules 2.0 (static file)
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
# Date: Mon 1 Oct 2012 14:58:26
from object_library import all_propagators, Propagator
# define only once the denominator since this is always the same
denominator = "P('mu', id) * P('mu', id) - Mass(id) * Mass(id) + complex(0,1) * Mass(id) * Width(id)"
# propagator for the scalar
S = Propagator(name = "S",
numerator = "complex(0,1)",
denominator = denominator
)
# propagator for the incoming fermion # the one for the outcomming is computed on the flight
F = Propagator(name = "F",
numerator = "complex(0,1) * (Gamma('mu', 1, 2) * P('mu', id) + Mass(id) * Identity(1, 2))",
denominator = denominator
)
# massive vector in the unitary gauge, can't be use for massless particles
V1 = Propagator(name = "V1",
numerator = "complex(0,1) * (-1 * Metric(1, 2) + Metric(1,'mu')* P('mu', id) * P(2, id) / Mass(id)**2 ",
denominator = denominator
)
# massless vector and massive vector in unitary gauge
V2 = Propagator(name = "V2",
numerator = "complex(0,-1) * Metric(1, 2)",
denominator = "P('mu', id) * P('mu', id)"
)
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO #####################
######################################################################
###################################
## INFORMATION FOR SMINPUTS
###################################
Block SMINPUTS
1 7.815553e-03 # aEW
2 1.166379e-05 # Gf
3 1.181000e-01 # aS
###################################
## INFORMATION FOR MASS
###################################
Block MASS
1 4.700000e-03 # MD
2 2.200000e-03 # MU
3 9.600000e-02 # MS
4 1.280000e+00 # MC
5 4.180000e+00 # MB
6 1.732000e+02 # MT
11 5.110000e-04 # Me
13 1.056600e-01 # MMU
15 1.777000e+00 # MTA
23 9.118760e+01 # MZ
25 1.250900e+02 # MH
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
22 0.000000 # a : 0.0
24 73.187688 # W+ : dMW + MW0
21 0.000000 # g : 0.0
9000001 0.000000 # ghA : 0.0
9000003 73.187688 # ghWp : dMW + MW0
9000004 73.187688 # ghWm : dMW + MW0
82 0.000000 # ghG : 0.0
12 0.000000 # ve : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
251 73.187688 # G+ : dMW + MW0
9000002 91.187600 # ghZ : MZ
250 91.187600 # G0 : MZ
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00
DECAY 23 2.495200e+00
DECAY 24 2.085000e+00
DECAY 25 4.070000e-03
## Not dependent paramater.
## Those values should be edited following analytical the
## analytical expression. Some generator could simply ignore
## those values and use the analytical expression
DECAY 22 0.000000 # a : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 9000001 0.000000 # ghA : 0.0
DECAY 82 0.000000 # ghG : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 1 0.000000 # d : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 9000002 2.495200 # ghZ : WZ
DECAY 9000003 2.085000 # ghWp : WW
DECAY 9000004 2.085000 # ghWm : WW
DECAY 250 2.495200 # G0 : WZ
DECAY 251 2.085000 # G+ : WW
###################################
## INFORMATION FOR CKMBLOCK
###################################
Block CKMBLOCK
1 2.277360e-01 # cabi
2 2.250600e-01 # CKMlambda
3 8.110000e-01 # CKMA
4 1.240000e-01 # CKMrho
5 0.000000e+00 # CKMeta
###################################
## INFORMATION FOR FRBLOCK
###################################
Block FRBlock
1 0.000000e+00 # ceWPh
2 0.000000e+00 # ceBPh
3 0.000000e+00 # cuGPh
4 0.000000e+00 # cuWPh
5 0.000000e+00 # cuBPh
6 0.000000e+00 # cdGPh
7 0.000000e+00 # cdWPh
8 0.000000e+00 # cdBPh
9 0.000000e+00 # cHudPh
10 0.000000e+00 # ceHPh
11 0.000000e+00 # cuHPh
12 0.000000e+00 # cdHPh
13 0.000000e+00 # cledqPh
14 0.000000e+00 # cquqd1Ph
15 0.000000e+00 # cquqd8Ph
16 0.000000e+00 # clequ1Ph
17 0.000000e+00 # clequ3Ph
18 1.000000e+03 # LambdaSMEFT
19 1.000000e+00 # cG
20 0.000000e+00 # cGtil
21 1.000000e+00 # cW
22 0.000000e+00 # cWtil
23 1.000000e+00 # cH
24 1.000000e+00 # cHbox
25 1.000000e+00 # cHDD
26 1.000000e+00 # cHG
27 0.000000e+00 # cHGtil
28 1.000000e+00 # cHW
29 0.000000e+00 # cHWtil
30 1.000000e+00 # cHB
31 0.000000e+00 # cHBtil
32 1.000000e+00 # cHWB
33 0.000000e+00 # cHWBtil
34 1.000000e+00 # ceHAbs
35 1.000000e+00 # cuHAbs
36 1.000000e+00 # cdHAbs
37 1.000000e+00 # ceWAbs
38 1.000000e+00 # ceBAbs
39 1.000000e+00 # cuGAbs
40 1.000000e+00 # cuWAbs
41 1.000000e+00 # cuBAbs
42 1.000000e+00 # cdGAbs
43 1.000000e+00 # cdWAbs
44 1.000000e+00 # cdBAbs
45 1.000000e+00 # cHl1
46 1.000000e+00 # cHl3
47 1.000000e+00 # cHe
48 1.000000e+00 # cHq1
49 1.000000e+00 # cHq3
50 1.000000e+00 # cHu
51 1.000000e+00 # cHd
52 1.000000e+00 # cHudAbs
53 1.000000e+00 # cll
54 1.000000e+00 # cll1
55 1.000000e+00 # cqq1
56 1.000000e+00 # cqq11
57 1.000000e+00 # cqq3
58 1.000000e+00 # cqq31
59 1.000000e+00 # clq1
60 1.000000e+00 # clq3
61 1.000000e+00 # cee
62 1.000000e+00 # cuu
63 1.000000e+00 # cuu1
64 1.000000e+00 # cdd
65 1.000000e+00 # cdd1
66 1.000000e+00 # ceu
67 1.000000e+00 # ced
68 1.000000e+00 # cud1
69 1.000000e+00 # cud8
70 1.000000e+00 # cle
71 1.000000e+00 # clu
72 1.000000e+00 # cld
73 1.000000e+00 # cqe
74 1.000000e+00 # cqu1
75 1.000000e+00 # cqu8
76 1.000000e+00 # cqd1
77 1.000000e+00 # cqd8
78 1.000000e+00 # cledqAbs
79 1.000000e+00 # cquqd1Abs
80 1.000000e+00 # cquqd8Abs
81 1.000000e+00 # clequ1Abs
82 1.000000e+00 # clequ3Abs
###################################
## INFORMATION FOR YUKAWA
###################################
Block YUKAWA
1 4.700000e-03 # ymdo
2 2.200000e-03 # ymup
3 9.600000e-02 # yms
4 1.280000e+00 # ymc
5 4.180000e+00 # ymb
6 1.732000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
15 1.777000e+00 # ymtau
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================
Block QNUMBERS 9000001 # ghA
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000002 # ghZ
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000003 # ghWp
1 3 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000004 # ghWm
1 -3 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 82 # ghG
1 0 # 3 times electric charge
2 -1 # number of spin states (2S+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 250 # G0
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 251 # G+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY THE UFO #####################
######################################################################
###################################
## INFORMATION FOR SMINPUTS
###################################
Block SMINPUTS
1 7.815553e-03 # aEW
2 1.166379e-05 # Gf
3 1.181000e-01 # aS
###################################
## INFORMATION FOR MASS
###################################
Block MASS
1 4.700000e-03 # MD
2 2.200000e-03 # MU
3 9.600000e-02 # MS
4 1.280000e+00 # MC