HepLib
Loading...
Searching...
No Matches
Init.cpp
Go to the documentation of this file.
1
6#include "BASIC.h"
7#include "HEP.h"
8#include "QGRAF.h"
9#include "QCD.h"
10#include "IBP.h"
11#include "SD.h"
12#include "DE.h"
13#include "AMF.h"
14#include <cstdlib>
15
16#include "cln/cln.h"
17#include <dlfcn.h>
18#include <string.h>
19#include <string>
20
21namespace HepLib {
22
23
24 namespace {
25 string install_prefix() {
26 static string wdir = "";
27 if(wdir == "") {
28 string path;
29 Dl_info dl_info;
30 dladdr((void*)install_prefix, &dl_info);
31 path = dl_info.dli_fname;
32 wdir = path.substr(0, path.find_last_of('/'));
33 path = wdir;
34 wdir = path.substr(0, path.find_last_of('/'));
35 }
36 return wdir;
37 }
38 }
39
40 // Initialization Unit
41 // to make sure the initialization order is correct,
42 // i.e., ordered in one compilation unit
43
44 // FROM class
45
46 // similar to GINAC_DECLARE_UNARCHIVER/GINAC_IMPLEMENT_REGISTERED_CLASS_OPT
47 GiNaC::registered_class_info Symbol::reg_info =
48 GiNaC::registered_class_info(GiNaC::registered_class_options("Symbol", "symbol", typeid(Symbol)));
49
50 GiNaC::registered_class_info iSymbol::reg_info =
51 GiNaC::registered_class_info(GiNaC::registered_class_options("iSymbol", "symbol", typeid(iSymbol)));
52
53 GiNaC::registered_class_info XIntegral::reg_info =
54 GiNaC::registered_class_info(GiNaC::registered_class_options("XIntegral", "basic", typeid(XIntegral)).print_func<print_dflt>(&XIntegral::print));
55
56 GiNaC::registered_class_info Index::reg_info =
57 GiNaC::registered_class_info(GiNaC::registered_class_options("Index", "basic", typeid(Index)).print_func<print_context>(&Index::print));
58
59 GiNaC::registered_class_info Vector::reg_info =
60 GiNaC::registered_class_info(GiNaC::registered_class_options("Vector", "basic", typeid(Vector)).print_func<print_context>(&Vector::print));
61
62 GiNaC::registered_class_info SUNT::reg_info =
63 GiNaC::registered_class_info(GiNaC::registered_class_options("SUNT", "basic", typeid(SUNT)).print_func<print_dflt>(&SUNT::print).print_func<FormFormat>(&SUNT::form_print).print_func<FCFormat>(&SUNT::fc_print));
64
65 GiNaC::registered_class_info SUNF::reg_info =
66 GiNaC::registered_class_info(GiNaC::registered_class_options("SUNF", "basic", typeid(SUNF)).print_func<print_dflt>(&SUNF::print).print_func<FormFormat>(&SUNF::form_print).print_func<FCFormat>(&SUNF::fc_print));
67
68 GiNaC::registered_class_info SUNF4::reg_info =
69 GiNaC::registered_class_info(GiNaC::registered_class_options("SUNF4", "basic", typeid(SUNF4)).print_func<print_dflt>(&SUNF4::print).print_func<FormFormat>(&SUNF4::form_print).print_func<FCFormat>(&SUNF4::fc_print));
70
71 GiNaC::registered_class_info DGamma::reg_info =
72 GiNaC::registered_class_info(GiNaC::registered_class_options("DGamma", "basic", typeid(DGamma)).print_func<print_dflt>(&DGamma::print).print_func<FormFormat>(&DGamma::form_print).print_func<FCFormat>(&DGamma::fc_print));
73
74 GiNaC::registered_class_info AsGamma::reg_info =
75 GiNaC::registered_class_info(GiNaC::registered_class_options("AsGamma", "basic", typeid(AsGamma)).print_func<print_dflt>(&AsGamma::print).print_func<FormFormat>(&AsGamma::form_print).print_func<FCFormat>(&AsGamma::fc_print));
76
77 GiNaC::registered_class_info Eps::reg_info =
78 GiNaC::registered_class_info(GiNaC::registered_class_options("Eps", "basic", typeid(Eps)).print_func<print_dflt>(&Eps::print).print_func<FormFormat>(&Eps::form_print).print_func<FCFormat>(&Eps::fc_print));
79
80 GiNaC::registered_class_info Pair::reg_info =
81 GiNaC::registered_class_info(GiNaC::registered_class_options("Pair", "basic", typeid(Pair)).print_func<print_dflt>(&Pair::print).print_func<FormFormat>(&Pair::form_print).print_func<FCFormat>(&Pair::fc_print));
82
83 // FROM BASIC
84
85 string Version = "1.5 @2024-04-28";
86 exmap Symbol::vmap;
87 std::map<std::string, ex> Symbol::Table; // alias as symtab in parser
88 std::map<std::string, ex> iSymbol::Table; // alias as symtab in parser
89 GiNaC::exmap Index::Dimension; // default dimension is d, others specified here.
90
91 unsigned nopat = GiNaC::subs_options::no_pattern;
92
93 ex w = wild();
94 ex w0 = wild(0);
95 ex w1 = wild(1);
96 ex w2 = wild(2);
97 ex w3 = wild(3);
98 ex w4 = wild(4);
99 ex w5 = wild(5);
100 ex w6 = wild(6);
101 ex w7 = wild(7);
102 ex w8 = wild(8);
103 ex w9 = wild(9);
104
105 // normal/factor options
106 int _o_ = 0;
107 const int o_none = (_o_++);
108 const int o_normal = (_o_++);
109 const int o_fermat = (_o_++);
110 const int o_fermatfD = (_o_++);
111 const int o_fermatN = (_o_++);
112 const int o_form = (_o_++);
113 const int o_flint = (_o_++);
114 const int o_flintf = (_o_++);
115 const int o_flintfD = (_o_++);
116 const int o_normal_fermat = (_o_++);
117 const int o_normal_form = (_o_++);
118 const int o_fermat_form = (_o_++);
119
120 // Symbols
121 const Symbol NA("NA");
122 const Symbol NF("NF");
123 const Symbol TF("TF");
124 const Symbol CA("CA");
125 const Symbol CF("CF");
126 const Symbol gs("gs");
127 const Symbol as("as");
128 const Symbol mu("mu");
129 const Symbol nL("nL");
130 const Symbol nH("nH");
131 const Symbol eps("eps");
132 const Symbol vs("vs");
133 const Symbol vz("vz");
134 const Symbol epz("epz");
135 const Symbol NaN("NaN");
136 const Symbol ep("ep");
137 const Symbol d("d");
138 const Symbol iet("iet");
139 const iSymbol iEpsilon("iEpsilon");
140
141 const ex iEpsilonN = I*pow(ex(10), -50);
142 int Verbose = 0;
143 string PRE = " ";
144 bool Debug = false;
145 bool In_GiNaC_Parallel = false;
147 map<string, int> GiNaC_Parallel_NP;
148 map<string, int> GiNaC_Parallel_Verb;
150 map<string, int> GiNaC_Parallel_NB;
151 map<string, bool> GiNaC_Parallel_RM;
152 map<string, string> GiNaC_Parallel_PRE;
153 map<string, bool> GiNaC_Parallel_ReWR;
155 map<ex,long long,ex_is_less> fermat_weight;
156 bool using_cache = true;
157 long long cache_limit = -1;
158 int NNDigits = 100;
159
160 bool GMat_using_cache = true;
161
163
164 string InstallPrefix = install_prefix();
165 string INC_FLAGS = "-I'"+InstallPrefix+"/include' " + "@INC_FLAGS@";
166 string LIB_FLAGS = "-L'"+InstallPrefix+"/lib' -Wl,-rpath,'"+InstallPrefix+"'/lib " + "@LIB_FLAGS@";
167
168 int Fermat::buffer_size = 1024*128;
169 int Form::buffer_size = 1024*128;
170
171 bool SD::SecDec::use_dlclose = true;
172 string SD::SecDec::cpp = "g++ -w";
173
174 SD::CppFormat::_init::_init() {
175 set_print_func<numeric, CppFormat>(CppFormat::print_numeric);
176 }
177 SD::CppFormat::_init SD::CppFormat::CppFormat_init;
179 set_print_func<numeric, ExFormat>(ExFormat::print_numeric);
180 }
181 SD::ExFormat::_init SD::ExFormat::ExFormat_init;
183
184 exmap SP_map;
185 map<ex,string,ex_is_less> QGRAF::LineTeX; // key is the filed
186 map<ex,string,ex_is_less> QGRAF::VerTeX; // key is the fileds in vertex
187 map<ex,string,ex_is_less> QGRAF::InOutTeX; // key is the id, id<0
188
189 FCFormat::_init::_init() {
190 set_print_func<ncmul, FCFormat>(FCFormat::ncmul_print);
191 }
192 FCFormat::_init FCFormat::FCFormat_init;
193
194 FormFormat::_init::_init() {
195 set_print_func<power, FormFormat>(FormFormat::power_print);
196 }
197 FormFormat::_init FormFormat::FormFormat_init;
199
200 const int form_trace_auto = 0;
201 const int form_trace_all = 1;
202 const int form_trace_each_all = 2;
203 const int form_trace_each_each = 3;
205
206 const int form_expand_none = 0;
207 const int form_expand_tr = 1;
208 const int form_expand_ci = 2;
209 const int form_expand_li = 3;
210 const int form_expand_all = 4;
213 bool form_using_su3 = false;
214 bool form_using_dim4 = false;
215 bool form_using_gamma5 = false;
216
217 QGRAF::Process::_init::_init() {
218 auto A = Symbol("A");
219 auto q = Symbol("q");
220 auto qbar = Symbol("qbar");
221 auto l = Symbol("l");
222 auto lbar = Symbol("lbar");
223 auto gh = Symbol("gh");
224 auto ghbar = Symbol("ghbar");
225 auto g = Symbol("g");
226 auto Q = Symbol("Q");
227 auto Qbar = Symbol("Qbar");
228 auto C = Symbol("C");
229 auto Cbar = Symbol("Cbar");
230 auto B = Symbol("B");
231 auto Bbar = Symbol("Bbar");
232 auto n = Symbol("n");
233 auto nbar = Symbol("nbar");
234 auto e = Symbol("e");
235 auto G = Symbol("G"); // graviton
236
237 LineTeX[q] = "fermion, edge label=$q$";
238 LineTeX[qbar] = "anti fermion, edge label=$q$";
239 LineTeX[l] = "fermion, edge label=$l$";
240 LineTeX[lbar] = "anti fermion, edge label=$l$";
241 LineTeX[gh] = "ghost, edge label=$\\chi$";
242 LineTeX[ghbar] = "ghost, edge label=$\\chi$";
243 LineTeX[g] = "gluon, edge label=$g$";
244 LineTeX[A] = "photon, edge label=$\\gamma$";
245 LineTeX[G] = "graviton, edge label=$\\gamma$";
246 LineTeX[Q] = "fermion, edge label=$Q$";
247 LineTeX[Qbar] = "anti fermion, edge label=$Q$";
248 LineTeX[C] = "fermion, edge label=$c$";
249 LineTeX[Cbar] = "anti fermion, edge label=$c$";
250 LineTeX[B] = "fermion, edge label=$b$";
251 LineTeX[Bbar] = "anti fermion, edge label=$b$";
252 LineTeX[n] = "double distance=1.5pt";
253 LineTeX[nbar] = "double distance=1.5pt";
254 LineTeX[e] = "color=white";
255
256 VerTeX[lst{Qbar, e, nbar}] = "[crossed dot]";
257 VerTeX[lst{nbar, e, g}] = "[crossed dot]";
258
259 {
260 Symbol U("U"); // U-quark
261 Symbol Ubar("Ubar"); // anti U-quark
262 Symbol D("D"); // D-quark
263 Symbol Dbar("Dbar"); // anti D-quark
264 Symbol C("C"); // C-quark
265 Symbol Cbar("Cbar"); // anti C-quark
266 Symbol S("S"); // S-quark
267 Symbol Sbar("Sbar"); // anti S-quark
268 Symbol T("T"); // T-quark
269 Symbol Tbar("Tbar"); // anti T-quark
270 Symbol B("B"); // B-quark
271 Symbol Bbar("Bbar"); // anti B-quark
272
273 Symbol g("g"); // gluon
274 Symbol gh("gh"); // gluon ghost
275 Symbol ghbar("ghbar"); // anti gluon ghost
276
277 Symbol A("A"); // photon
278 Symbol Wm("Wm"); // W-
279 Symbol Wp("Wp"); // W+
280 Symbol Z("Z"); // Z
281
282 Symbol ghA("SW"); // photon ghost
283 Symbol ghAbar("ghAbar"); // anti photon ghost
284 Symbol ghWm("ghWm"); // W- ghost
285 Symbol ghWmbar("ghWmbar"); // anti W- ghost
286 Symbol ghWp("ghWp"); // W+ ghost
287 Symbol ghWpbar("ghWpbar"); // anti W+ ghost
288 Symbol ghZ("ghZ"); // Z ghost
289 Symbol ghZbar("ghZbar"); // anti Z ghost
290
291 Symbol em("em"); // e-
292 Symbol ep("ep"); // e+
293 Symbol ne("ne"); // e-neutrino
294 Symbol nebar("nebar"); // anti e-neutrino
295 Symbol mum("mum"); // mu-
296 Symbol mup("mup"); // mu+
297 Symbol nmu("nmu"); // mu-neutrino
298 Symbol nmubar("nmubar"); // anti mu-neutrino
299 Symbol taum("taum"); // tau-
300 Symbol taup("taup"); // tau+
301 Symbol ntau("ntau"); // tau-neutrino
302 Symbol ntaubar("ntaubar"); // anti tau-neutrino
303
304 Symbol chi("chi"); // Z goldstone
305 Symbol phim("phim"); // W- goldstone
306 Symbol phip("phip"); // W+ goldstone
307 Symbol H("H"); // higgs
308
309 lst fs = { U, D, C, S, T, B, em, mum, taum, ne, nmu, ntau};
310 lst fs2 = {Ubar, Dbar, Cbar, Sbar, Tbar, Bbar, ep, mup, taup, nebar, nmubar, ntaubar};
311 string fsi[] = { "$u$", "$d$", "$c$", "$s$", "$t$", "$b$", "$e$", "$\\mu$", "$\\tau$", "$\\nu_e$", "$\\nu_\\mu$", "$\\nu_\\tau$" };
312 for(int i=0; i<fs.nops(); i++) {
313 auto si = fsi[i];
314 LineTeX[fs.op(i)] = "fermion ,edge label="+si;
315 LineTeX[fs2.op(i)] = "anti fermion ,edge label="+si;
316 }
317
318 LineTeX[Wp] = "photon, edge label=$W$";
319 LineTeX[Wm] = "photon, edge label=$W$";
320 LineTeX[Z] = "photon, edge label=$Z$";
321
322 LineTeX[H]="scalar, edge label=$H$";
323 LineTeX[chi]="scalar, edge label=$\\chi$";
324 LineTeX[phim]="charged scalar, edge label=$\\phi$";
325 LineTeX[phip]="anti charged scalar, edge label=$\\phi$";
326
327 lst ghs = {gh, ghA, ghWm, ghWp, ghZ};
328 lst ghs2 = {ghbar, ghAbar, ghWmbar, ghWpbar, ghZbar};
329 string ghsi[] = {"$c$", "$c_\\gamma$", "$c_-$", "$c_+$", "$c_Z$"};
330
331 for(int i=0; i<ghs.nops(); i++) {
332 auto si = ghsi[i];
333 LineTeX[ghs.op(i)] = "ghost ,edge label="+si;
334 LineTeX[ghs2.op(i)] = "anti ghost ,edge label="+si;
335 }
336 }
337 }
338 QGRAF::Process::_init QGRAF::Process::Process_init;
339
340 // FROM IBP
341 int FIRE::Version = 6;
342 exmap MapPreSP;
343
344 string UKIRA::KArgs = "";
345 string KIRA::KArgs = "";
346
347 // FROM QCD
348 int QCD::FF::cur_mode = 0; // 0 - gluon, 1 - quark, 2 - anti-quark
349
350 // FROM DE
351 slong error_pass_dp = 100;
352
353 // FROM AMF
354 int DEX::Threads = 8;
355
356 // gamma matrix
357 ex DGamma::gi = GAS(1);
358 ex DGamma::g5 = GAS(5);
359 ex DGamma::C = GAS(int('c'));
360
361 // global init
362 std::stack<cln::float_format_t> cln_prec_stack;
363 std::stack<long> digits_stack;
364 _global_init::_init::_init() {
365 ostringstream oss;
366 string path = InstallPrefix + "/bin";
367 if(dir_exists(path)) oss << path;
368
369 auto opath = getenv("PATH");
370 if(opath != NULL) oss << ":" << opath;
371 setenv("PATH", oss.str().c_str(), true);
372
373 // similar to GINAC_DECLARE_UNARCHIVER/GINAC_BIND_UNARCHIVER
374 GiNaC::unarchive_table_t table;
375 table.insert(std::string("Symbol"), []()->GiNaC::basic*{ return new Symbol(); });
376 table.insert(std::string("iSymbol"), []()->GiNaC::basic*{ return new iSymbol(); });
377 table.insert(std::string("XIntegral"), []()->GiNaC::basic*{ return new XIntegral(); });
378 table.insert(std::string("Index"), []()->GiNaC::basic*{ return new Index(); });
379 table.insert(std::string("Vector"), []()->GiNaC::basic*{ return new Vector(); });
380 table.insert(std::string("SUNT"), []()->GiNaC::basic*{ return new SUNT(); });
381 table.insert(std::string("SUNF"), []()->GiNaC::basic*{ return new SUNF(); });
382 table.insert(std::string("SUNF4"), []()->GiNaC::basic*{ return new SUNF4(); });
383 table.insert(std::string("DGamma"), []()->GiNaC::basic*{ return new DGamma(); });
384 table.insert(std::string("AsGamma"), []()->GiNaC::basic*{ return new AsGamma(); });
385 table.insert(std::string("Eps"), []()->GiNaC::basic*{ return new Eps(); });
386 table.insert(std::string("Pair"), []()->GiNaC::basic*{ return new Pair(); });
387
388 // CLN configurations
389 cln::cl_inhibit_floating_point_underflow = true;
390 Digits = 100;
391 set_precision(100);
392 }
393 _global_init::_init _global_init::init_object;
394
395}
396
397std::string HepLib::QGRAF::Process::Style = R"EOF(
398<prologue>
399
400<diagram>
401((<sign><symmetry_factor>)*
402<in_loop>InField(<field>,<field_index>,<momentum>)*
403<end>
404<back><out_loop>OutField(<field>,<field_index>,<momentum>)*
405<end>
406<back><propagator_loop>Propagator(
407<back>Field(<field>,<field_index>),
408<back>Field(<dual-field>,<dual-field_index>),
409<back><momentum>)*
410<end>
411<back><vertex_loop>Vertex(
412<back><ray_loop>Field(<field>,<field_index>,<momentum>),
413<back><end><back>)*
414<end><back><back>)
415,
416<epilogue>
417
418<exit>
419
420)EOF";
421
422std::string HepLib::QGRAF::Models::GluonFF = R"EOF(
423[ model = 'Gluon FF Model' ]
424[q, qbar, -]
425[Q, Qbar, -]
426[gh, ghbar, -]
427[g, g, +, notadpole]
428[e, e, +, external]
429[n, nbar, +]
430[qbar, q, g; QCD='+1']
431[Qbar, Q, g; QCD='+1']
432[g, g, g, g; QCD='+2']
433[g, g, g; QCD='+1']
434[ghbar, gh, g; QCD='+1']
435[nbar, n, g; QCD='+1']
436[nbar, e, g; QCD='+0']
437%[n, e, g; QCD='+0']
438)EOF";
439
440std::string HepLib::QGRAF::Models::QuarkFF = R"EOF(
441[ model = 'Quark FF Model' ]
442[q, qbar, -]
443[Q, Qbar, -]
444[e, ebar, -, external]
445[gh, ghbar, -]
446[g, g, +, notadpole]
447[n, nbar, +]
448[qbar, q, g; QCD='+1']
449[Qbar, Q, g; QCD='+1']
450[g, g, g, g; QCD='+2']
451[g, g, g; QCD='+1']
452[ghbar, gh, g; QCD='+1']
453[nbar, n, g; QCD='+1']
454[qbar, e, nbar; QCD='0']
455[ebar, q, nbar; QCD='0']
456[Qbar, e, nbar; QCD='0']
457[ebar, Q, nbar; QCD='0']
458)EOF";
459
460std::string HepLib::QGRAF::Models::SM = R"EOF(
461% From Boehm, Denner, Joos @ Feynman Gauge
462[ model = 'Standard Model' ]
463
464%------------------------------------
465% Propagators
466%------------------------------------
467
468% quarks
469[U, Ubar, -]
470[D, Dbar, -]
471[C, Cbar, -]
472[S, Sbar, -]
473[T, Tbar, -]
474[B, Bbar, -]
475
476% gluon and its ghost:
477[g, g, +, notadpole]
478[gh, ghbar, -]
479
480% V: EWSM gauge bosons
481[A, A, +, notadpole]
482[Wm, Wp, +]
483[Z, Z, +]
484
485% leptons
486[em, ep, -]
487[ne, nebar, -]
488[mum, mup, -]
489[nmu, nmubar, -]
490[taum, taup, -]
491[ntau, ntaubar, -]
492
493% G: Faddeev-Popov Ghosts
494[ghA, ghAbar, -]
495[ghZ, ghZbar, -]
496[ghWm, ghWmbar, -]
497[ghWp, ghWpbar, -]
498
499% S: scalars
500[H, H, +]
501[chi, chi, +]
502[phim, phip, +]
503
504%------------------------------------
505% Vertices
506%------------------------------------
507
508% QCD
509[Ubar, U, g; QCD='+1', QED='0']
510[Dbar, D, g; QCD='+1', QED='0']
511[Sbar, S, g; QCD='+1', QED='0']
512[Cbar, C, g; QCD='+1', QED='0']
513[Bbar, B, g; QCD='+1', QED='0']
514[Tbar, T, g; QCD='+1', QED='0']
515
516[g, g, g, g; QCD='+2', QED='0']
517[g, g, g; QCD='+1', QED='0']
518[ghbar, gh, g; QCD='+1', QED='0']
519
520% VVVV
521[Wp, Wm, Z, Z; QCD='0', QED='+2']
522[Wp, Wm, A, Z; QCD='0', QED='+2']
523[Wp, Wm, A, A; QCD='0', QED='+2']
524[Wp, Wp, Wm, Wm; QCD='0', QED='+2']
525
526% VVV
527[Wp, Wm, A; QCD='0', QED='+1']
528[Wp, Wm, Z; QCD='0', QED='+1']
529
530% SSSS
531[H, H, H, H; QCD='0', QED='+2']
532[chi, chi, chi, chi; QCD='0', QED='+2']
533[H, H, chi, chi; QCD='0', QED='+2']
534[H, H, phim, phip; QCD='0', QED='+2']
535[chi, chi, phim, phip; QCD='0', QED='+2']
536[phim, phip, phim, phip; QCD='0', QED='+2']
537
538% SSS
539[H, H, H; QCD='0', QED='+1']
540[H, chi, chi; QCD='0', QED='+1']
541[H, phim, phip; QCD='0', QED='+1']
542
543% VVSS
544[Z, Z, H, H; QCD='0', QED='+2']
545[Z, Z, chi, chi; QCD='0', QED='+2']
546[Wp, Wm, H, H; QCD='0', QED='+2']
547[Wp, Wm, chi, chi; QCD='0', QED='+2']
548[Wp, Wm, phip, phim; QCD='0', QED='+2']
549
550[A, A, phip, phim; QCD='0', QED='+2']
551[Z, A, phip, phim; QCD='0', QED='+2']
552[Z, Z, phip, phim; QCD='0', QED='+2']
553
554[Wp, A, phim, H; QCD='0', QED='+2']
555[Wm, A, phip, H; QCD='0', QED='+2']
556[Wp, A, phim, chi; QCD='0', QED='+2']
557[Wm, A, phip, chi; QCD='0', QED='+2']
558[Wp, Z, phim, H; QCD='0', QED='+2']
559[Wm, Z, phip, H; QCD='0', QED='+2']
560[Wp, Z, phim, chi; QCD='0', QED='+2']
561[Wm, Z, phip, chi; QCD='0', QED='+2']
562
563% VSS
564[Z, chi, H; QCD='0', QED='+1']
565[A, phip, phim; QCD='0', QED='+1']
566[Z, phip, phim; QCD='0', QED='+1']
567[Wp, phim, H; QCD='0', QED='+1']
568[Wm, phip, H; QCD='0', QED='+1']
569[Wp, phim, chi; QCD='0', QED='+1']
570[Wm, phip, chi; QCD='0', QED='+1']
571
572% SVV
573[H, Z, Z; QCD='0', QED='+1']
574[H, Wp, Wm; QCD='0', QED='+1']
575[phip, Wm, A; QCD='0', QED='+1']
576[phim, Wp, A; QCD='0', QED='+1']
577[phip, Wm, Z; QCD='0', QED='+1']
578[phim, Wp, Z; QCD='0', QED='+1']
579
580% ffV
581[Ubar, U, A; QCD='0', QED='+1']
582[Dbar, D, A; QCD='0', QED='+1']
583[Sbar, S, A; QCD='0', QED='+1']
584[Cbar, C, A; QCD='0', QED='+1']
585[Bbar, B, A; QCD='0', QED='+1']
586[Tbar, T, A; QCD='0', QED='+1']
587
588[ep, em, A; QCD='0', QED='+1']
589[mup, mum, A; QCD='0', QED='+1']
590[taup, taum, A; QCD='0', QED='+1']
591
592[Ubar, U, Z; QCD='0', QED='+1']
593[Dbar, D, Z; QCD='0', QED='+1']
594[Sbar, S, Z; QCD='0', QED='+1']
595[Cbar, C, Z; QCD='0', QED='+1']
596[Bbar, B, Z; QCD='0', QED='+1']
597[Tbar, T, Z; QCD='0', QED='+1']
598
599[ep, em, Z; QCD='0', QED='+1']
600[mup, mum, Z; QCD='0', QED='+1']
601[taup, taum, Z; QCD='0', QED='+1']
602
603[nebar, ne, Z; QCD='0', QED='+1']
604[nmubar, nmu, Z; QCD='0', QED='+1']
605[ntaubar, ntau, Z; QCD='0', QED='+1']
606
607[Dbar, U, Wm; QCD='0', QED='+1']
608[Dbar, C, Wm; QCD='0', QED='+1']
609[Dbar, T, Wm; QCD='0', QED='+1']
610
611[Sbar, U, Wm; QCD='0', QED='+1']
612[Sbar, C, Wm; QCD='0', QED='+1']
613[Sbar, T, Wm; QCD='0', QED='+1']
614
615[Bbar, U, Wm; QCD='0', QED='+1']
616[Bbar, C, Wm; QCD='0', QED='+1']
617[Bbar, T, Wm; QCD='0', QED='+1']
618
619[Ubar, D, Wp; QCD='0', QED='+1']
620[Ubar, S, Wp; QCD='0', QED='+1']
621[Ubar, B, Wp; QCD='0', QED='+1']
622
623[Cbar, D, Wp; QCD='0', QED='+1']
624[Cbar, S, Wp; QCD='0', QED='+1']
625[Cbar, B, Wp; QCD='0', QED='+1']
626
627[Tbar, D, Wp; QCD='0', QED='+1']
628[Tbar, S, Wp; QCD='0', QED='+1']
629[Tbar, B, Wp; QCD='0', QED='+1']
630
631[ep, ne, Wm; QCD='0', QED='+1']
632[mup, nmu, Wm; QCD='0', QED='+1']
633[taup, ntau, Wm; QCD='0', QED='+1']
634
635[nebar, em, Wp; QCD='0', QED='+1']
636[nmubar, mum, Wp; QCD='0', QED='+1']
637[ntaubar, taum, Wp; QCD='0', QED='+1']
638
639% Sff
640[Ubar, U, H; QCD='0', QED='+1']
641[Dbar, D, H; QCD='0', QED='+1']
642[Sbar, S, H; QCD='0', QED='+1']
643[Cbar, C, H; QCD='0', QED='+1']
644[Bbar, B, H; QCD='0', QED='+1']
645[Tbar, T, H; QCD='0', QED='+1']
646
647[ep, em, H; QCD='0', QED='+1']
648[mup, mum, H; QCD='0', QED='+1']
649[taup, taum, H; QCD='0', QED='+1']
650
651[Ubar, U, chi; QCD='0', QED='+1']
652[Dbar, D, chi; QCD='0', QED='+1']
653[Sbar, S, chi; QCD='0', QED='+1']
654[Cbar, C, chi; QCD='0', QED='+1']
655[Bbar, B, chi; QCD='0', QED='+1']
656[Tbar, T, chi; QCD='0', QED='+1']
657
658[ep, em, chi; QCD='0', QED='+1']
659[mup, mum, chi; QCD='0', QED='+1']
660[taup, taum, chi; QCD='0', QED='+1']
661
662[Dbar, U, phim; QCD='0', QED='+1']
663[Dbar, C, phim; QCD='0', QED='+1']
664[Dbar, T, phim; QCD='0', QED='+1']
665
666[Sbar, U, phim; QCD='0', QED='+1']
667[Sbar, C, phim; QCD='0', QED='+1']
668[Sbar, T, phim; QCD='0', QED='+1']
669
670[Bbar, U, phim; QCD='0', QED='+1']
671[Bbar, C, phim; QCD='0', QED='+1']
672[Bbar, T, phim; QCD='0', QED='+1']
673
674[Ubar, D, phip; QCD='0', QED='+1']
675[Ubar, S, phip; QCD='0', QED='+1']
676[Ubar, B, phip; QCD='0', QED='+1']
677
678[Cbar, D, phip; QCD='0', QED='+1']
679[Cbar, S, phip; QCD='0', QED='+1']
680[Cbar, B, phip; QCD='0', QED='+1']
681
682[Tbar, D, phip; QCD='0', QED='+1']
683[Tbar, S, phip; QCD='0', QED='+1']
684[Tbar, B, phip; QCD='0', QED='+1']
685
686[ep, ne, phim; QCD='0', QED='+1']
687[mup, nmu, phim; QCD='0', QED='+1']
688[taup, ntau, phim; QCD='0', QED='+1']
689
690[nebar, em, phip; QCD='0', QED='+1']
691[nmubar, mum, phip; QCD='0', QED='+1']
692[ntaubar, taum, phip; QCD='0', QED='+1']
693
694% VGG
695[ghWpbar, ghWp, A; QCD='0', QED='+1']
696[ghWmbar, ghWm, A; QCD='0', QED='+1']
697[ghAbar, ghWm, Wp; QCD='0', QED='+1']
698[ghAbar, ghWp, Wm; QCD='0', QED='+1']
699[ghWmbar, ghA, Wm; QCD='0', QED='+1']
700[ghWpbar, ghA, Wp; QCD='0', QED='+1']
701[ghWpbar, ghWp, Z; QCD='0', QED='+1']
702[ghWmbar, ghWm, Z; QCD='0', QED='+1']
703[ghZbar, ghWm, Wp; QCD='0', QED='+1']
704[ghZbar, ghWp, Wm; QCD='0', QED='+1']
705[ghWmbar, ghZ, Wm; QCD='0', QED='+1']
706[ghWpbar, ghZ, Wp; QCD='0', QED='+1']
707
708% SGG
709[ghZbar, ghZ, H; QCD='0', QED='+1']
710[ghWpbar, ghWp, H; QCD='0', QED='+1']
711[ghWmbar, ghWm, H; QCD='0', QED='+1']
712[ghWpbar, ghWp, chi; QCD='0', QED='+1']
713[ghWmbar, ghWm, chi; QCD='0', QED='+1']
714[ghWpbar, ghA, phip; QCD='0', QED='+1']
715[ghWmbar, ghA, phim; QCD='0', QED='+1']
716[ghWpbar, ghZ, phip; QCD='0', QED='+1']
717[ghWmbar, ghZ, phim; QCD='0', QED='+1']
718[ghZbar, ghWp, phim; QCD='0', QED='+1']
719[ghZbar, ghWm, phip; QCD='0', QED='+1']
720
721)EOF";
722
723std::string HepLib::QGRAF::Models::QCD = R"EOF(
724[ model = 'qcd Model' ]
725%------------------------------------
726% Propagators
727%------------------------------------
728% quark
729[q, qbar, -]
730[Q, Qbar, -]
731% gluon and its ghost:
732[gh, ghbar, -]
733[g, g, +, notadpole]
734% external photon
735[A, A, +, external]
736%------------------------------------
737% Vertices
738%------------------------------------
739% QCD
740[qbar, q, g; QCD='+1']
741[Qbar, Q, g; QCD='+1']
742[g, g, g, g; QCD='+2']
743[g, g, g; QCD='+1']
744[ghbar, gh, g; QCD='+1']
745% external
746[qbar, q, A; QCD='+0']
747[Qbar, Q, A; QCD='+0']
748)EOF";
Basic header file.
HEP header file.
IBP header file.
QGRAF header file.
SecDec header file.
class for AntiSymmetric Gamma object
Definition HEP.h:496
void form_print(const FormFormat &c, unsigned level=0) const
Definition AsGamma.cpp:150
void fc_print(const FCFormat &c, unsigned level=0) const
Definition AsGamma.cpp:154
void print(const print_dflt &c, unsigned level=0) const
Definition AsGamma.cpp:135
class for Dirac Gamma object
Definition HEP.h:437
static ex C
Definition HEP.h:487
void form_print(const FormFormat &c, unsigned level=0) const
Definition DGamma.cpp:139
void print(const print_dflt &c, unsigned level=0) const
Definition DGamma.cpp:117
void fc_print(const FCFormat &c, unsigned level=0) const
Definition DGamma.cpp:150
static ex g5
Definition HEP.h:486
static ex gi
Definition HEP.h:485
class for Levi-Civita object https://onlinelibrary.wiley.com/doi/pdf/10.1002/9783527630097....
Definition HEP.h:386
void fc_print(const FCFormat &c, unsigned level=0) const
Definition Eps.cpp:121
void form_print(const FormFormat &c, unsigned level=0) const
Eps in FORM format https://onlinelibrary.wiley.com/doi/pdf/10.1002/9783527630097.app3 to make Tr(g5,...
Definition Eps.cpp:112
void print(const print_dflt &c, unsigned level=0) const
Definition Eps.cpp:96
class for FCFormat Output
Definition HEP.h:71
static void ncmul_print(const ncmul &p, const FCFormat &c, unsigned level=0)
Definition Basic.cpp:163
static int Version
Definition IBP.h:86
static int buffer_size
Definition BASIC.h:800
static void power_print(const power &p, const FormFormat &c, unsigned level=0)
Definition Basic.cpp:55
static int buffer_size
Definition BASIC.h:822
class for index object
Definition HEP.h:104
void print(const print_context &c, unsigned level=0) const
Definition Basic.cpp:198
static GiNaC::exmap Dimension
Definition HEP.h:109
static string KArgs
Definition IBP.h:97
class for Mathematica Format Output
Definition BASIC.h:843
class for Pair object
Definition HEP.h:322
void form_print(const FormFormat &c, unsigned level=0) const
FormFormat print function.
Definition Pair.cpp:102
void print(const print_dflt &c, unsigned level=0) const
default print function
Definition Pair.cpp:93
void fc_print(const FCFormat &c, unsigned level=0) const
FCFormat print function.
Definition Pair.cpp:113
static string QuarkFF
Definition QGRAF.h:163
static string SM
Definition QGRAF.h:164
static string QCD
Definition QGRAF.h:165
static string GluonFF
Definition QGRAF.h:162
static string Style
Definition QGRAF.h:77
static bool use_dlclose
Definition SD.h:416
static string cpp
Definition SD.h:417
class for SUNF4 object
Definition HEP.h:278
void form_print(const FormFormat &c, unsigned level=0) const
print the Form Format
Definition Basic.cpp:549
void fc_print(const FCFormat &c, unsigned level=0) const
Definition Basic.cpp:540
void print(const print_dflt &c, unsigned level=0) const
normal priint
Definition Basic.cpp:536
class for SUNF object
Definition HEP.h:233
void print(const print_dflt &c, unsigned level=0) const
Definition Basic.cpp:437
void form_print(const FormFormat &c, unsigned level=0) const
Definition Basic.cpp:441
void fc_print(const FCFormat &c, unsigned level=0) const
Definition Basic.cpp:445
class for SUNT object
Definition HEP.h:189
void form_print(const FormFormat &c, unsigned level=0) const
Definition Basic.cpp:334
void fc_print(const FCFormat &c, unsigned level=0) const
Definition Basic.cpp:345
void print(const print_dflt &c, unsigned level=0) const
Definition Basic.cpp:349
class extended to GiNaC symbol class, represent a positive symbol
Definition BASIC.h:113
static std::map< std::string, ex > Table
Definition BASIC.h:160
static exmap vmap
Definition BASIC.h:167
static string KArgs
Definition IBP.h:121
class for vector object
Definition HEP.h:149
void print(const print_context &c, unsigned level=0) const
Definition Basic.cpp:275
XIntegral Class, preface to SecDec.
Definition BASIC.h:740
void print(const print_dflt &c, unsigned level=0) const
Definition BASIC.cpp:1674
class extended to GiNaC symbol class, pure imaginary symbol
Definition BASIC.h:183
static std::map< std::string, ex > Table
Definition BASIC.h:230
map< ex, string, ex_is_less > LineTeX
Definition Init.cpp:185
map< ex, string, ex_is_less > VerTeX
Definition Init.cpp:186
map< ex, string, ex_is_less > InOutTeX
Definition Init.cpp:187
int VEO_Digits
Definition Init.cpp:182
HepLib namespace.
Definition BASIC.cpp:17
string INC_FLAGS
Definition Init.cpp:165
bool form_using_dim4
Definition Init.cpp:214
const Symbol NA
ex w9
Definition BASIC.h:499
ex w8
Definition BASIC.h:499
const int form_expand_tr
Definition Init.cpp:207
const iSymbol iEpsilon
ex w7
Definition BASIC.h:499
MMAFormat mout
ex w5
Definition BASIC.h:499
int GiNaC_Parallel_Process
Definition Init.cpp:146
int _o_
Definition Init.cpp:106
const Symbol gs
const int o_form
Definition Init.cpp:112
int NNDigits
Definition Init.cpp:158
const int o_fermat
Definition Init.cpp:109
const int o_none
Definition Init.cpp:107
const int o_flint
Definition Init.cpp:113
exmap MapPreSP
Definition Init.cpp:342
const int form_expand_none
Definition Init.cpp:206
string Version
Definition Init.cpp:85
ex w0
Definition BASIC.h:499
const Symbol NF
map< string, bool > GiNaC_Parallel_RM
Definition Init.cpp:151
exmap SP_map
Definition Init.cpp:184
map< ex, long long, ex_is_less > fermat_weight
Definition Init.cpp:155
const Symbol ep
const int o_fermatfD
Definition Init.cpp:110
int form_expand_mode
Definition Init.cpp:211
const Symbol eps
bool Apart_using_fermat
Definition Init.cpp:212
map< string, int > GiNaC_Parallel_NB
Definition Init.cpp:150
const int form_trace_all
Definition Init.cpp:201
ex GAS(const ex &expr, unsigned rl)
function similar to GAD/GSD in FeynClac
Definition DGamma.cpp:280
const int o_fermat_form
Definition Init.cpp:118
ex w4
Definition BASIC.h:499
const int form_expand_ci
Definition Init.cpp:208
slong error_pass_dp
Definition Init.cpp:351
int form_trace_mode
Definition Init.cpp:204
bool form_using_su3
Definition Init.cpp:213
string LIB_FLAGS
Definition Init.cpp:166
bool using_cache
Definition Init.cpp:156
bool Debug
Definition Init.cpp:144
ex w
Definition Init.cpp:93
const int form_trace_each_all
Definition Init.cpp:202
const Symbol vs
const Symbol d
int fermat_using_array
Definition Init.cpp:154
map< string, int > GiNaC_Parallel_NP
Definition Init.cpp:147
const int form_expand_all
Definition Init.cpp:210
map< string, int > GiNaC_Parallel_Verb
Definition Init.cpp:148
map< string, bool > GiNaC_Parallel_ReWR
Definition Init.cpp:153
const int o_normal
Definition Init.cpp:108
const Symbol mu
const Symbol CA
const Symbol TF
const int o_flintf
Definition Init.cpp:114
const Symbol CF
const int o_flintfD
Definition Init.cpp:115
const int o_normal_fermat
Definition Init.cpp:116
string InstallPrefix
Definition Init.cpp:164
const ex iEpsilonN
Definition Init.cpp:141
FCFormat fcout
const Symbol epz
bool In_GiNaC_Parallel
Definition Init.cpp:145
const Symbol as
bool form_using_gamma5
Definition Init.cpp:215
int Verbose
Definition Init.cpp:142
const int form_trace_each_each
Definition Init.cpp:203
const int form_trace_auto
Definition Init.cpp:200
bool dir_exists(string dir)
Definition BASIC.h:297
const int o_normal_form
Definition Init.cpp:117
long long cache_limit
Definition Init.cpp:157
const Symbol nL
const Symbol nH
ex w1
Definition BASIC.h:499
unsigned nopat
Definition Init.cpp:91
ex w3
Definition BASIC.h:499
void set_precision(long prec, bool push)
Definition BASIC.cpp:2314
const Symbol iet
ex w6
Definition BASIC.h:499
map< string, string > GiNaC_Parallel_PRE
Definition Init.cpp:152
const int form_expand_li
Definition Init.cpp:209
const int o_fermatN
Definition Init.cpp:111
std::stack< cln::float_format_t > cln_prec_stack
Definition Init.cpp:362
ex w2
Definition BASIC.h:499
std::stack< long > digits_stack
Definition Init.cpp:363
bool GMat_using_cache
Definition Init.cpp:160
const Symbol vz
string PRE
Definition Init.cpp:143
const Symbol NaN
int GiNaC_Parallel_Batch
Definition Init.cpp:149