QCD BGF model for FeynArts
Load FeynRules
FR$Parallel = False;
$FeynRulesPath = FileNameJoin[{$UserBaseDirectory, "Applications", "FeynRules"}];
<< 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, "QCDBGF.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 in the background field formalism loaded.
Generate Feynman rules
fRules = FeynmanRules[LQCD]
Starting Feynman rule calculation.
Expanding the Lagrangian...
Collecting the different structures that enter the vertex.
17 possible non-zero vertices have been found -> starting the computation: FeynRulesˋFR$FeynmanRules / 17.
17 vertices obtained.
BBB123BBG123ghG†ghGB123BBBB1234ghG†ghGBB1234BGG123BBBG1234ghG†ghGG123ghG†ghGBG1234BBGG1234GGG123BGGG1234GGGG1234dq−dqB123uq−uqB123dq−dqG123uq−uqG123−gsfa1,a2,a3ημ1,μ2p1μ3+gsfa1,a2,a3ημ1,μ2p2μ3+gsfa1,a2,a3ημ1,μ3p1μ2−gsfa1,a2,a3ημ1,μ3p3μ2−gsfa1,a2,a3ημ2,μ3p2μ1+gsfa1,a2,a3ημ2,μ3p3μ1GaugeXi(G)gsfa1,a2,a3ημ1,μ3p2μ2−GaugeXi(G)gsfa1,a2,a3ημ2,μ3p1μ1−gsfa1,a2,a3ημ1,μ2p1μ3+gsfa1,a2,a3ημ1,μ2p2μ3+gsfa1,a2,a3ημ1,μ3p1μ2−gsfa1,a2,a3ημ1,μ3p3μ2−gsfa1,a2,a3ημ2,μ3p2μ1+gsfa1,a2,a3ημ2,μ3p3μ1gsfa3,a1,a2p1μ3−gsfa3,a1,a2p2μ3igs2ημ1,μ4ημ2,μ3fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a2,Gluon$1fa3,a4,Gluon$1+igs2ημ1,μ3ημ2,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ3ημ2,μ4fa1,a2,Gluon$1fa3,a4,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1igs2ημ3,μ4fa3,a2,Gluon$1fa4,a1,Gluon$1+igs2ημ3,μ4fa3,a1,Gluon$1fa4,a2,Gluon$1GaugeXi(G)gsfa1,a2,a3ημ1,μ2p3μ3−GaugeXi(G)gsfa1,a2,a3ημ1,μ3p2μ2−gsfa1,a2,a3ημ1,μ2p1μ3+gsfa1,a2,a3ημ1,μ2p2μ3+gsfa1,a2,a3ημ1,μ3p1μ2−gsfa1,a2,a3ημ1,μ3p3μ2−gsfa1,a2,a3ημ2,μ3p2μ1+gsfa1,a2,a3ημ2,μ3p3μ1igs2ημ1,μ4ημ2,μ3fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a2,Gluon$1fa3,a4,Gluon$1+igs2ημ1,μ3ημ2,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ3ημ2,μ4fa1,a2,Gluon$1fa3,a4,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1gsfa3,a1,a2p1μ3igs2ημ3,μ4fa3,a1,Gluon$1fa4,a2,Gluon$1−GaugeXi(G)igs2ημ1,μ4ημ2,μ3fa1,a4,Gluon$1fa2,a3,Gluon$1−GaugeXi(G)igs2ημ1,μ3ημ2,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a2,Gluon$1fa3,a4,Gluon$1+igs2ημ1,μ3ημ2,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ3ημ2,μ4fa1,a2,Gluon$1fa3,a4,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1−gsfa1,a2,a3ημ1,μ2p1μ3+gsfa1,a2,a3ημ1,μ2p2μ3+gsfa1,a2,a3ημ1,μ3p1μ2−gsfa1,a2,a3ημ1,μ3p3μ2−gsfa1,a2,a3ημ2,μ3p2μ1+gsfa1,a2,a3ημ2,μ3p3μ1igs2ημ1,μ4ημ2,μ3fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a2,Gluon$1fa3,a4,Gluon$1+igs2ημ1,μ3ημ2,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ3ημ2,μ4fa1,a2,Gluon$1fa3,a4,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1igs2ημ1,μ4ημ2,μ3fa1,a3,Gluon$1fa2,a4,Gluon$1+igs2ημ1,μ4ημ2,μ3fa1,a2,Gluon$1fa3,a4,Gluon$1+igs2ημ1,μ3ημ2,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ3ημ2,μ4fa1,a2,Gluon$1fa3,a4,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a4,Gluon$1fa2,a3,Gluon$1−igs2ημ1,μ2ημ3,μ4fa1,a3,Gluon$1fa2,a4,Gluon$1−igsδf1,f2Tm1,m2a3γs1,s2μ3−igsδf1,f2Tm1,m2a3γs1,s2μ3−igsδf1,f2Tm1,m2a3γs1,s2μ3−igsδf1,f2Tm1,m2a3γs1,s2μ3
Create FeynArts model
SetDirectory[FileNameJoin[{$UserBaseDirectory, "Applications", "FeynCalc", "FeynArts", "Models"}]];
WriteFeynArtsOutput[LQCD, Output -> "QCDBGF", CouplingRename -> False,SelectParticles -> {
{ghG, ghGbar, B}, {ghG, ghGbar, B, B}, {B, G, G},
{ghG, ghGbar, G}, {ghG, ghGbar, B, G}, {B, B, G, G},
{G, G, G}, {B, G, G, G}, {G, G, G, G}, {uqbar, uq, G}, {dqbar, dq, G},
{uqbar, uq, B}, {dqbar, dq, B}}];
- - - 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...
Selecting specified field content. Warning! Only mass eigenstates should be selected!
Neglecting all terms with more than 4 particles.
Neglecting all terms with less than 3 particles.
Collecting the different structures that enter the vertex.
13 possible non-zero vertices have been found -> starting the computation: FeynRulesˋFR$FeynmanRules / 13.
13 vertices obtained.
mytimecheck,after LGC
Writing FeynArts model file into directory QCDBGF
Writing FeynArts generic file on QCDBGF.gen.