Mathematica Six gluon scattering amplitude
There are some files to download
integrals.math ;
treefunctions.math ;
egzdata.math ;
function_simple.math ;
ghfunction.math ;
kinematic6pointEGZ
The mathematica code to run these looks like
Npoints=6
<< kinematic6pointEGZ
<< function_simple.math
<< integrals.math
<< treefunctions.math
<< egzdata.math
<< ghfunction.math
<< Apppppp.math
<< Amppppp.math
<< Ammpppp.math
<< Ampmppp.math
<< Amppmpp.math
<< Ammmppp.math
<< Ammpmpp.math
<< Ampmpmp.math
After this has been run Eval[AmmppppN1] will numerically
evaluate the AmmpppppN1 amplitude at the kinematic point
of Ellis Geile and Zandereghi. AN1[1] gives the EGZ
point for comparisom.
Ammpppp = Amplitude 1;
Ampmppp = Amplitude 2;
Amppmpp = Amplitude 3;
Ammmppp = Amplitude 4;
Ammpmpp = Amplitude 5;
Ampmpmp = Amplitude 6;
In no sense is any of this code claimed to be optimal !!!!
The code run on Mathematica 5 and 6.The code was originally written in
5. The analytic continuation
needed for the integrals needed some minor adjustments before working
on 6...