% !TEX program = pdflatexmk
% !TEX parameter = -shell-escape
% Author: Izaak Neutelings (October 24)
% Instructions: To compile via command line, run the following twice
%   pdflatex -shell-escape extract/higgs_pp$TAG.tex
\documentclass[11pt,border=4pt]{standalone}
\usepackage{feynmp-auto}
\usepackage{xcolor}
\definecolor{colggH}{rgb}{0,0,1} % gluon fusion (blue)
\definecolor{colkappaf}{rgb}{1,0,0} % kappa_q (red)
\definecolor{colkappaV}{rgb}{.15,.75,.15} % kappa_V (green)
\begin{document}
\begin{fmffile}{feyngraph}
  \fmfframe(4,14)(13,10){ % padding (L,T)(R,B)
  %\color{colggH} % blue text color
  %\fmfcmd{foreground:=(0,0,1);} % blue line color
  \begin{fmfgraph*}(170,100) % dimensions (WH)
    % line style
    \fmfset{wiggly_len}{12} % boson wavelength
    \fmfset{wiggly_slope}{65} % boson slope of waves
    % external vertices
    \fmfstraight
    \fmfleft{i2,i1}
    \fmfright{o4,o3,o2,o1}
    \fmfforce{(0,.71h)}{i1} % exact placement
    \fmfforce{(0,.29h)}{i2} % exact placement
    % internal vertices (exact placement)
    \fmfforce{(.24w,.71h)}{d1} % exact placement
    \fmfforce{(.24w,.29h)}{d2} % exact placement
    \fmfforce{(.84w,.82h)}{v1} % exact placement
    \fmfforce{(.84w,.18h)}{v2} % exact placement
    % gluons
    \fmf{gluon}{i2,d2}
    \fmf{gluon}{d1,i1}
    % top triangle loop
    \fmf{fermion}{d1,t3,d2}
    \fmf{fermion,l.s=left,label=t}{d2,d1}
    % Higgs boson
    \fmf{dashes,t=1.4,l.s=left,label=h}{t3,h}
    \fmf{boson,l.s=left,label=Z}{v2,h,v1}
    \fmf{fermion}{o2,v1,o1}
    \fmf{fermion}{o4,v2,o3}
    % Higss coupling
    \fmfv{decor.shape=circle,decor.filled=full,decor.size=4,f=(1,,0,,0),
          l.d=5,l.a=-60,l=\color{colkappaf}$\kappa_\mathrm{t}$}{t3}
    \fmfv{decor.shape=circle,decor.filled=full,decor.size=4,f=(.15,,.75,,.15),
          l.d=5,l.a=-110,l=\color{colkappaV}$\kappa_\mathrm{Z}$\hspace{4pt}}{h}
    % labels
    \fmflabel{$g$}{i1}
    \fmflabel{$g$}{i2}
    \fmfv{l.a=24,l.d=5,l=$\mu^-$}{o1}
    \fmfv{l.a=-22,l.d=4,l=$\mu^+$}{o2}
    \fmfv{l.a=24,l.d=5,l=$\mu^-$}{o3}
    \fmfv{l.a=-22,l.d=4,l=$\mu^+$}{o4}
  \end{fmfgraph*}
  } % close \fmfframe
\end{fmffile}
\end{document}