~maddevelopers/mg5amcnlo/2.9.4

« back to all changes in this revision

Viewing changes to tests/input_files/sm_with_custom_propa/propagator.py

  • Committer: olivier Mattelaer
  • Date: 2013-04-01 15:53:26 UTC
  • mfrom: (240.1.21 1.5.9)
  • Revision ID: olivier.mattelaer@uclouvain.be-20130401155326-tngbopy89f9jg0tw
Pass to 1.5.9

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
# This file was automatically created by FeynRules 1.7.69
 
2
# Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
 
3
# Date: Mon 1 Oct 2012 14:58:26
 
4
 
 
5
from object_library import all_propagator, Propagator
 
6
 
 
7
# define only once the denominator since this is always the same
 
8
denominator = "P('mu', id) * P('mu', id) - Mass(id) * Mass(id) + complex(0,1) * Mass(id) * Width(id)"
 
9
 
 
10
# propagator for the scalar
 
11
S = Propagator(name = "S",
 
12
               spin = 1,
 
13
               numerator = "complex(0,1)",
 
14
               denominator = denominator
 
15
               )
 
16
 
 
17
# propagator for the incoming fermion # the one for the outcomming is computed on the flight
 
18
F = Propagator(name = "F",
 
19
                spin = 2,
 
20
                numerator = "complex(0,1) * (Gamma('mu', s1, s2) * P('mu', id) + Mass(id) * Identity(s1, s2))",
 
21
                denominator = denominator
 
22
               )
 
23
 
 
24
# massive vector in the unitary gauge, can't be use for massless particles
 
25
V1 = Propagator(name = "V1",
 
26
                spin = 3,
 
27
                numerator = "complex(0,1) * (-1 * Metric(l1, l2) + Metric(l1,'mu')* P('mu', id) * P(l2, id) / Mass(id)**2 ",
 
28
                denominator = denominator
 
29
               )
 
30
 
 
31
# massless vector and massive vector in unitary gauge
 
32
V2 = Propagator(name = "V2",
 
33
                spin = 3,
 
34
                expression = "complex(0,-1) * Metric(l1, l2)",
 
35
                denominator = denominator
 
36
               )
 
37
 
 
38