def get_diagram_one(ax, fermion_style, boson_style, vertex_style):
D = Diagram(ax)
w = 0.75
xy0 = [0.5 - w/2, 0.25]
v1 = D.vertex(xy0, **vertex_style)
v2 = D.vertex(v1.xy, dx=w, **vertex_style)
G = D.line(v1, v2, **fermion_style)
B = D.line(v1, v2, **boson_style)
# In case the axes get smaller (you have more diagrams), you might want to change the scale
D.scale(1.0)
D.plot()
return D
""" VBH === Vector Boson Fusion. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10., 10.)) ax = fig.add_axes([0, 0, 1, 1], frameon=False) diagram = Diagram(ax) diagram.text(.5, 0.9, "Vector Boson Fusion (VBF)", fontsize=40) in1 = diagram.vertex(xy=(.1, .8), marker='') in2 = diagram.vertex(xy=(.1, .2), marker='') v1 = diagram.vertex(xy=(.4, .7)) v2 = diagram.vertex(xy=(.4, .3)) v3 = diagram.vertex(xy=(.6, .5)) out1 = diagram.vertex(xy=(.9, .8), marker='') out2 = diagram.vertex(xy=(.9, .2), marker='') higgsout = diagram.vertex(xy=(.9, .5)) q1 = diagram.line(in1, v1, arrow=False) q2 = diagram.line(in2, v2, arrow=False) wz1 = diagram.line(v1, v3, style='wiggly') wz2 = diagram.line(v2, v3, style='wiggly') higgs = diagram.line(v3, higgsout, style='dashed', arrow=False) q3 = diagram.line(v1, out1, arrow=False) q4 = diagram.line(v2, out2, arrow=False)
"""
VH
==
Higgs Strahlung.
"""
import matplotlib.pyplot as plt
from feynman import Diagram
fig = plt.figure(figsize=(10.,10.))
ax = fig.add_axes([0,0,1,1], frameon=False)
diagram = Diagram(ax)
diagram.text(.4,0.9,"Associated Vector Boson", fontsize=40)
diagram.text(.6,0.83,"(VH or 'Higgs Strahlung')", fontsize=40)
in1 = diagram.vertex(xy=(.1,.75), marker='')
in2= diagram.vertex(xy=(.1,.25), marker='')
v1 = diagram.vertex(xy=(.35,.5))
v2 = diagram.vertex(xy=(.65,.5))
higgsout = diagram.vertex(xy=(.9,.75))
out1 = diagram.vertex(xy=(.9,.25),marker='')
q1 = diagram.line(in1, v1)
q2 = diagram.line(v1, in2)
wz1 = diagram.line(v1, v2, style='wiggly')
wz2 = diagram.line(v2, out1, style='wiggly')
higgs = diagram.line(v2, higgsout, arrow=False, style='dashed')
q1.text("q",fontsize=30)
q2.text(r"$\bar{\mathrm{q}}$",fontsize=30)
diagram.text(0.5,0.55,"$Z/W^\pm$",fontsize=30)
"""
Common neutrino nucleon interactions
"""
import matplotlib.pyplot as plt
from feynman import Diagram
fig = plt.figure(figsize=(10, 5))
ax1 = fig.add_subplot(121, frameon=False)
diagram = Diagram(ax1)
diagram.text(.5, 0.9, "Charged-Current", fontsize=25)
in1 = diagram.vertex(xy=(0.05, 0.75), marker='')
in2 = diagram.vertex(xy=(0.05, 0.25), marker='')
v1 = diagram.vertex(xy=(.45, .6))
v2 = diagram.vertex(xy=(.45, .4))
out1 = diagram.vertex(xy=(0.85, 0.75), marker='')
out2 = diagram.vertex(xy=(0.85, 0.25), marker='')
nu1 = diagram.line(in1, v1) #incoming neutrino
N = diagram.line(in2, v2) #incoming neucleon
W = diagram.line(v1, v2, style='wiggly') #W mediator
nu2 = diagram.line(v1, out1) #outgoing neutrino
X = diagram.line(v2, out2) #outgoing shower
diagram.text(0.10, 0.68, "$\\nu_{\ell}$", fontsize=30)
diagram.text(0.10, 0.32, "$N$", fontsize=30)
diagram.text(0.57, 0.5, "$W^{\pm}$", fontsize=30)
diagram.text(0.81, 0.68, "$\ell$", fontsize=30)
diagram.text(0.81, 0.32, "$X$", fontsize=30)
ax2 = fig.add_subplot(122, frameon=False)
diagram2 = Diagram(ax2)
diagram2.text(.5, 0.9, "Neutral-Current", fontsize=25)
=====
Doubly Charged Higgs Production
"""
import matplotlib
matplotlib.rcParams['mathtext.fontset'] = 'stix'
matplotlib.rcParams['font.family'] = 'STIXGeneral'
from feynman import Diagram
fig = matplotlib.pyplot.figure(figsize=(10., 10.))
ax = fig.add_axes([0, 0, 1, 1], frameon=False)
diagram = Diagram(ax)
diagram.text(.4, 0.9, "Doubly Charged Higgs Production", fontsize=40)
in1 = diagram.vertex(xy=(.1, .75), marker='')
in2 = diagram.vertex(xy=(.1, .25), marker='')
v1 = diagram.vertex(xy=(.35, .5))
v2 = diagram.vertex(xy=(.65, .5))
higgsplusout = diagram.vertex(xy=(.8, .7))
higgsminusout = diagram.vertex(xy=(.8, .3))
l1plus = diagram.vertex(xy=(.95, .8), marker='')
l2plus = diagram.vertex(xy=(.95, .6), marker='')
l1minus = diagram.vertex(xy=(.95, .4), marker='')
l2minus = diagram.vertex(xy=(.95, .2), marker='')
lw = 5
q1 = diagram.line(v1,
in1,
color='blue',
lw=lw,
""" ggFZZ EFT ========= The gluon-gluon fusion. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10., 10.)) ax = fig.add_axes([0, 0, 1, 1], frameon=False) diagram = Diagram(ax) #diagram.text(.5,0.9,"Gluon-Gluon Fusion (ggF)",fontsize=40) in1 = diagram.vertex(xy=(.05, .7), marker='') in2 = diagram.vertex(xy=(.05, .3), marker='') v1 = diagram.vertex(xy=(.25, .7)) v2 = diagram.vertex(xy=(.25, .3)) v3 = diagram.vertex(xy=(.45, .5)) higgsout = diagram.vertex(xy=(.60, .5)) zout1 = diagram.vertex(xy=(.85, .7), marker='') zout2 = diagram.vertex(xy=(.85, .3), marker='') gluon_style = dict(style='linear loopy', xamp=.025, yamp=.035, nloops=4) g1 = diagram.line(in1, v1, **gluon_style) g2 = diagram.line(in2, v2, **gluon_style) t1 = diagram.line(v1, v2) t2 = diagram.line(v2, v3) t3 = diagram.line(v3, v1) higgs = diagram.line(v3, higgsout, arrow=False, style='dashed')
"""
ggF EFT
=========
The gluon-gluon fusion.
"""
import matplotlib.pyplot as plt
from feynman import Diagram
fig = plt.figure(figsize=(10.,10.))
ax = fig.add_axes([0,0,1,1], frameon=False)
diagram = Diagram(ax)
in1 = diagram.vertex(xy=(.1,.6), marker='')
in2= diagram.vertex(xy=(.1,.4), marker='')
v1 = diagram.vertex(xy=(.4,.6))
v2 = diagram.vertex(xy=(.4,.4))
v3 = diagram.vertex(xy=(.6,.5))
v4 = diagram.vertex(xy=(.34,.5), marker='')
higgsout = diagram.vertex(xy=(.9,.5))
epsilon = diagram.operator([v4,v3], c=1.1)
epsilon.text("Effective \n coupling", fontsize=30)
gluon_up_style = dict(style='linear loopy', xamp=.025, yamp=.035, nloops=7)
gluon_down_style = dict(style='linear loopy', xamp=.025, yamp=-.035, nloops=7)
g1 = diagram.line(in1, v1, **gluon_up_style)
g2 = diagram.line(in2, v2, **gluon_down_style)
higgs = diagram.line(v3, higgsout, arrow=False, style='dashed')
side = 0.3
gammalen = side * np.sqrt(3) / 2
linlen = 0.3
txtpad = 0.2
tail_marker = 'o'
W_style = dict(style = 'double wiggly', nwiggles=2)
v_style = dict(style = 'simple wiggly', nwiggles=2)
G_style = dict(style = 'double', arrow=True, arrow_param={'width':0.05})
D = Diagram(ax)
# Left hand size
xy = [0.4, y0]
v11 = D.vertex(xy, dy= side/2)
v12 = D.vertex(xy, dy=-side/2)
v13 = D.vertex(xy, dx=gammalen)
gamma0 = D.operator([v11,v12,v13])
gamma0.text("$\Gamma$")
# Symbol
D.text(v13.x + txtpad, y0, "=")
# Create a three-vertex dot.
chunkdist = .03
chunklen = .03
chunkstyle=dict(arrow=False, linewidth=6.)
v20 = D.vertex([v13.x + 2 * txtpad, y0])
v210 = D.vertex(v20.xy, angle=0., radius=chunkdist, marker='')
v220 = D.vertex(v20.xy, angle=1./3, radius=chunkdist, marker='')
D.x0 = 0.2
D.y0 = sum(D.ax.get_ylim()) * .35
# Various size
opwidth = 1.
linlen = 2.
txtpad = .8
wiggle_amplitude=.1
# Line styles
Ph_style = dict(style='elliptic loopy', ellipse_spread=.6, xamp=.10, yamp=-.15, nloops=15)
DW_style = dict(style='circular loopy', circle_radius=.7, xamp=.10, yamp=.15, nloops=18)
G_style = dict(style='simple', arrow=True, arrow_param={'width':0.15, 'length': .3})
# Item 1
v11 = D.vertex([D.x0, D.y0])
v12 = D.vertex(v11.xy, dx=opwidth)
Sigma = D.operator([v11, v12])
Sigma.text("$\Sigma^{ep}$")
# Symbol
D.text(v12.x + txtpad, D.y0, "=")
# Item 3
v21 = D.vertex([v12.x + 2 * txtpad, D.y0 - 0.3])
v22 = D.vertex(v21.xy, dx=linlen)
G = D.line(v21, v22, **G_style)
Ph = D.line(v21, v22, **Ph_style)
# Symbol
D.text(v22.x + txtpad, D.y0, "+")
""" VBH tau-tau =========== Vector Boson Fusion. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10., 10.)) ax = fig.add_axes([0, 0, 1, 1], frameon=False) diagram = Diagram(ax) #diagram.text(.5,0.9,r"Vector Boson Fusion (VBF) Higgs $\rightarrow\tau\tau$",fontsize=40) in1 = diagram.vertex(xy=(.1, .8), marker='') in2 = diagram.vertex(xy=(.1, .2), marker='') v1 = diagram.vertex(xy=(.3, .7)) v2 = diagram.vertex(xy=(.3, .3)) v3 = diagram.vertex(xy=(.5, .5)) out1 = diagram.vertex(xy=(.9, .8), marker='') out2 = diagram.vertex(xy=(.9, .2), marker='') higgsf = diagram.vertex(xy=(.7, .5)) tau1 = diagram.vertex(xy=(.9, .7), marker='') tau2 = diagram.vertex(xy=(.9, .3), marker='') q1 = diagram.line(in1, v1, arrow=False) q2 = diagram.line(in2, v2, arrow=False) wz1 = diagram.line(v1, v3, style='wiggly') wz2 = diagram.line(v2, v3, style='wiggly') higgs = diagram.line(v3, higgsf, style='dashed', arrow=False) q3 = diagram.line(v1, out1, arrow=False)
""" LFV === The LFV diagram. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10.,10.)) ax = fig.add_axes([0,0,1,1], frameon=False) diagram = Diagram(ax) in1 = diagram.vertex(xy=(.1,.5)) in2= diagram.vertex(xy=(.4,.5)) v1 = diagram.vertex(xy=(.65,.65)) v2 = diagram.vertex(xy=(.65,.35)) out1 = diagram.vertex(xy=(.9,.65),marker='') out2 = diagram.vertex(xy=(.9,.35),marker='') higgs = diagram.line(in1, in2, arrow=False, style='dashed') nu1 = diagram.line(v1, in2) nu2 = diagram.line(in2, v2) w = diagram.line(v1, v2, style='wiggly') lep = diagram.line(out1, v1) tau = diagram.line(v2, out2) nu1.text(r"$\nu_\ell$",fontsize=40) nu2.text(r"$\nu_\tau$",fontsize=40) lep.text(r"$\ell^+$",fontsize=40) tau.text(r"$\tau^-$",fontsize=40)
# Set up the figure and ax
fig = plt.figure(figsize=(8,2))
ax = fig.add_axes([0,0,1,1], frameon=False)
# It is best to keep the xlim/ylim ratio the same as the figure aspect ratio.
ax.set_xlim(0, 1)
ax.set_ylim(0, 0.25)
y0 = sum(ax.get_ylim()) / 2
# Initialize diagram with the ax
D = Diagram(ax)
# Polarizability operator
v11 = D.vertex([0.1, y0])
v12 = D.vertex(v11.xy, dx=0.15)
P = D.operator([v11, v12], c=1.3) # c is the excentricity of the patch
P.text("$P$")
# Symbols
D.text(.35, y0, "=", fontsize=30)
# Propagator lines
G_style = dict(style='double elliptic',
ellipse_excentricity=-1.2, ellipse_spread=.3,
arrow=True, arrow_param={'width':0.03})
v21 = D.vertex([0.45, y0])
v22 = D.vertex(v21.xy, dx=0.4)
ax.set_xlim(0, 1.0)
ax.set_ylim(0, .15)
l = 0.15 # Length of the propagator
txt_l = 0.05 # Padding around the symbol
op_l = 0.08 # Size of the operator
G_style = dict(arrow=True, arrow_param={'width':0.02, 'length': 0.05}, style = 'double')
G0_style = dict(arrow=True, arrow_param={'width':0.02, 'length': 0.05}, style = 'simple')
text_prop = dict(y=0.02, fontsize=20)
D = Diagram(ax)
# Left hand side
v11 = D.vertex(xy=[0.05, 0.06])
v12 = D.vertex(v11.xy, dx=l)
G = D.line(v11, v12, **G_style)
G.text("$G$", **text_prop)
# Symbol
D.text(v12.x + txt_l, v12.y, "=")
# First term
v21 = D.vertex(v12.xy, dx=2*txt_l)
v22 = D.vertex(v21.xy, dx=l)
G0 = D.line(v21, v22, **G0_style)
G0.text("$G_0$", **text_prop)
# Symbol
D.text(v22.x + txt_l, v22.y, "+")
===== Doubly Charged Higgs Production """ import matplotlib matplotlib.rcParams['mathtext.fontset'] = 'stix' matplotlib.rcParams['font.family'] = 'STIXGeneral' from feynman import Diagram fig = matplotlib.pyplot.figure(figsize=(10.,10.)) ax = fig.add_axes([0,0,1,1], frameon=False) diagram = Diagram(ax) diagram.text(.4,0.9,"Doubly Charged Higgs Production", fontsize=40) q1in = diagram.vertex(xy=(.1,.75), marker='') q2in= diagram.vertex(xy=(.1,.25), marker='') v1 = diagram.vertex(xy=(.3,.75)) v2 = diagram.vertex(xy=(.3,.25)) vmerge = diagram.vertex(xy=(.6,.5)) higgsout = diagram.vertex(xy=(.8,.5)) q1pout = diagram.vertex(xy=(.95,.75), marker='') q2pout = diagram.vertex(xy=(.95,.25), marker='') l1out = diagram.vertex(xy=(.95,.62), marker='') l2out = diagram.vertex(xy=(.95,.38), marker='') lw = 5 # Quarks q1 = diagram.line(q1in, v1, color='blue', lw=lw, arrow_param=dict(color='blue', length=0.08, width=0.02)) q2 = diagram.line(q2in, v2, color='blue', lw=lw, arrow_param=dict(color='blue', length=0.08, width=0.02)) q1p = diagram.line(v1, q1pout, color='blue', lw=lw, arrow_param=dict(color='blue', length=0.08, width=0.02))
""" VBH tau-tau =========== Vector Boson Fusion. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10.,10.)) ax = fig.add_axes([0,0,1,1], frameon=False) diagram = Diagram(ax) #diagram.text(.5,0.9,r"Vector Boson Fusion (VBF) Higgs $\rightarrow\tau\tau$",fontsize=40) in1 = diagram.vertex(xy=(.1,.8), marker='') in2= diagram.vertex(xy=(.1,.2), marker='') v1 = diagram.vertex(xy=(.3,.7)) v2 = diagram.vertex(xy=(.3,.3)) v3 = diagram.vertex(xy=(.5,.5)) out1 = diagram.vertex(xy=(.9,.8), marker='') out2 = diagram.vertex(xy=(.9,.2), marker='') higgsf = diagram.vertex(xy=(.7,.5)) tau1 = diagram.vertex(xy=(.9,.7), marker='') tau2 = diagram.vertex(xy=(.9,.3), marker='') q1 = diagram.line(in1, v1, arrow=False) q2 = diagram.line(in2, v2, arrow=False) wz1 = diagram.line(v1, v3, style='wiggly') wz2 = diagram.line(v2, v3, style='wiggly') higgs = diagram.line(v3, higgsf, style='dashed', arrow=False) q3 = diagram.line(v1, out1, arrow=False)
A diagram containing double lines.
"""
import matplotlib.pyplot as plt
from feynman import Diagram
fig = plt.figure(figsize=(8,2))
ax = fig.add_axes([0,0,1,1], frameon=False)
ax.set_xlim(0, 1)
ax.set_ylim(0,.25)
# Sigma operator
D = Diagram(ax)
v11 = D.vertex([.1, .08])
v12 = D.vertex(v11.xy, dx=.15)
Sigma = D.operator([v11, v12])
Sigma.text("$\Sigma$")
# Symbols
D.text(v12.x+.1, v12.y, "=")
# GW convolution
v21 = D.vertex(v12.xy, dxy=[0.2, -0.04])
v22 = D.vertex(v21.xy, dx=0.4)
l21 = D.line(v21, v22, style='double', arrow=True)
# Specifying the number of wiggles and the amplitude of the wiggles
"""
Fock exchange
=============
A simple diagram composed of vertices and lines.
"""
import matplotlib.pyplot as plt
from feynman import Diagram
# If no Axes is given, a new one is initialized.
diagram = Diagram()
# Create four vertices.
v1 = diagram.vertex(xy=(.1, .5), marker='')
v2 = diagram.vertex(v1.xy, dx=.2)
v3 = diagram.vertex(v2.xy, dx=.4)
v4 = diagram.vertex(v3.xy, dx=.2, marker='')
# Create four lines.
# By default, 'simple' lines have arrows
# and others flavours such as 'wiggly' or 'loopy' don't.
l12 = diagram.line(v1, v2)
l23 = diagram.line(v2, v3)
l34 = diagram.line(v3, v4, arrow=True)
w23 = diagram.line(v2, v3, style='wiggly elliptic')
# Add labels.
l12.text("p")
w23.text("q")
l23.text("p - q")
Gamma_side = .4
Gamma_height = Gamma_side * np.sqrt(3) / 2
# Positon of the first vertex
x0 = 0.2
y0 = sum(d.ax.get_ylim()) / 2
# Define line styles
G_style = dict(
style='double elliptic',
ellipse_excentricity=-1.2, ellipse_spread=.3,
arrow=True, arrow_param={'width':0.05},
)
# Draw the diagram
v01 = d.vertex([x0, y0])
v02 = d.vertex(v01.xy, dx=opwidth)
P = d.operator([v01,v02], c=1.3)
P.text("$P$")
d.text(v02.x+textpad, y0, "=", fontsize=30)
v21 = d.vertex(v02.xy, dx=.4)
v22 = d.vertex(v21.xy, dx=linlen, dy=Gamma_side/2)
v23 = d.vertex(v21.xy, dx=linlen, dy=-Gamma_side/2)
v24 = d.vertex(v21.xy, dx=linlen+Gamma_height)
l21 = d.line(v22, v21, **G_style)
l21 = d.line(v21, v23, **G_style)
Gamma = d.operator([v22,v23, v24])
""" ggFZZ EFT ========= The gluon-gluon fusion. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10.,10.)) ax = fig.add_axes([0,0,1,1], frameon=False) diagram = Diagram(ax) #diagram.text(.5,0.9,"Gluon-Gluon Fusion (ggF)",fontsize=40) in1 = diagram.vertex(xy=(.05,.7), marker='') in2= diagram.vertex(xy=(.05,.3), marker='') v1 = diagram.vertex(xy=(.25,.7)) v2 = diagram.vertex(xy=(.25,.3)) v3 = diagram.vertex(xy=(.45,.5)) higgsout = diagram.vertex(xy=(.60,.5)) zout1 = diagram.vertex(xy=(.85,.7), marker='') zout2 = diagram.vertex(xy=(.85,.3), marker='') gluon_style = dict(style='linear loopy', xamp=.025, yamp=.035, nloops=4) g1 = diagram.line(in1, v1, **gluon_style) g2 = diagram.line(in2, v2, **gluon_style) t1 = diagram.line(v1, v2) t2 = diagram.line(v2, v3) t3 = diagram.line(v3, v1) higgs = diagram.line(v3, higgsout, arrow=False, style='dashed')
""" ttH === The ttH diagram. """ import matplotlib.pyplot as plt from feynman import Diagram fig = plt.figure(figsize=(10.,10.)) ax = fig.add_axes([0,0,1,1], frameon=False) diagram = Diagram(ax) diagram.text(.5,0.9,"Associated Top Pair (ttH)", fontsize=40) in1 = diagram.vertex(xy=(.1,.8), marker='') in2= diagram.vertex(xy=(.1,.2), marker='') v1 = diagram.vertex(xy=(.4,.7)) v2 = diagram.vertex(xy=(.4,.3)) v3 = diagram.vertex(xy=(.6,.5)) out1 = diagram.vertex(xy=(.9,.8), marker='') out2 = diagram.vertex(xy=(.9,.2), marker='') higgsout = diagram.vertex(xy=(.9,.5)) g1 = diagram.line(in1, v1, style='loopy',nloops=7,yamp=0.04) g2 = diagram.line(in2, v2, style='loopy',nloops=7,yamp=0.04) t1 = diagram.line(v3, v1, arrow = True) t2 = diagram.line(v2, v3, arrow = True) higgs = diagram.line(v3, higgsout, arrow=False, style='dashed') t3 = diagram.line(v1, out1, arrow=True) t4 = diagram.line(out2, v2, arrow=True)
"""
Fock exchange
=============
A simple diagram composed of vertices and lines.
"""
import matplotlib.pyplot as plt
from feynman import Diagram
# If no Axes is given, a new one is initialized.
diagram = Diagram()
# Create four vertices.
v1 = diagram.vertex(xy=(.1,.5), marker='')
v2 = diagram.vertex(v1.xy, dx=.2)
v3 = diagram.vertex(v2.xy, dx=.4)
v4 = diagram.vertex(v3.xy, dx=.2, marker='')
# Create four lines.
# By default, 'simple' lines have arrows
# and others flavours such as 'wiggly' or 'loopy' don't.
l12 = diagram.line(v1, v2)
l23 = diagram.line(v2, v3)
l34 = diagram.line(v3, v4, arrow=True)
w23 = diagram.line(v2, v3, style='wiggly elliptic')
# Add labels.
l12.text("p")
w23.text("q")
l23.text("p - q")
