QCD model for FeynArts
Load FeynRules
FR$Parallel = False;
$FeynRulesPath = FileNameJoin[{$UserBaseDirectory, "Applications", "FeynRules"}]
<< FeynRules`;
/home/vs/.Mathematica/Applications/FeynRules
- FeynRules -
Version: 2.3.49 (29 September 2021).
Authors: A. Alloul, N. Christensen, C. Degrande, C. Duhr, B. Fuks
$$$$
Please cite:
- Comput.Phys.Commun.185:2250-2300,2014 (arXiv:1310.1921);
- Comput.Phys.Commun.180:1614-1641,2009 (arXiv:0806.4194).
$$$$
http://feynrules.phys.ucl.ac.be
$$$$
The FeynRules palette can be opened using the command FRPalette[].
Load FeynRules model
If[$FrontEnd === Null,
nbDir = DirectoryName[$InputFileName],
nbDir = NotebookDirectory[]
];
frModelPath = FileNameJoin[{nbDir, "QCD.fr"}];
LoadModel[frModelPath];
This model implementation was created by
Vladyslav Shtabovenko
Model Version: 0
For more information, type ModelInformation[].
- Loading particle classes.
- Loading gauge group classes.
- Loading parameter classes.
\nModel QCD loaded.
FR$Loop = True;
SetDirectory[FileNameJoin[{$UserBaseDirectory, "Applications", "FeynCalc", "FeynArts", "Models"}]];
WriteFeynArtsOutput[LQCD, Output -> "QCD", CouplingRename -> False];
- - - FeynRules interface to FeynArts - - -
C. Degrande C. Duhr, 2013
Counterterms: B. Fuks, 2012
Calculating Feynman rules for L1
Starting Feynman rules calculation for L1.
Expanding the Lagrangian...
Collecting the different structures that enter the vertex.
9 possible non-zero vertices have been found -> starting the computation: FeynRulesˋFR$FeynmanRules / 9.
9 vertices obtained.
mytimecheck,after LGC
Writing FeynArts model file into directory QCD
Writing FeynArts generic file on QCD.gen.
input = FileNameJoin[{$UserBaseDirectory, "Applications", "FeynCalc", "FeynArts", "Models", "QCD", "QCD.gen"}];
tmp = Import[input, "Text"] <> "\n";
res = StringJoin[{tmp, "GaugeXi[V[5,___]]=GaugeXi[\"G\"]; \n GaugeXi[U[5,___]]=GaugeXi[\"G\"]; \n"}];
Export[input, res, "Text"];