25 string install_prefix() {
26 static string wdir =
"";
30 dladdr((
void*)install_prefix, &dl_info);
31 path = dl_info.dli_fname;
32 wdir = path.substr(0, path.find_last_of(
'/'));
34 wdir = path.substr(0, path.find_last_of(
'/'));
47 GiNaC::registered_class_info Symbol::reg_info =
48 GiNaC::registered_class_info(GiNaC::registered_class_options(
"Symbol",
"symbol",
typeid(Symbol)));
50 GiNaC::registered_class_info iSymbol::reg_info =
51 GiNaC::registered_class_info(GiNaC::registered_class_options(
"iSymbol",
"symbol",
typeid(iSymbol)));
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));
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));
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));
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));
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));
68 GiNaC::registered_class_info SUNF4::reg_info =
71 GiNaC::registered_class_info DGamma::reg_info =
74 GiNaC::registered_class_info AsGamma::reg_info =
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));
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));
91 unsigned nopat = GiNaC::subs_options::no_pattern;
174 SD::CppFormat::_init::_init() {
175 set_print_func<numeric, CppFormat>(CppFormat::print_numeric);
177 SD::CppFormat::_init SD::CppFormat::CppFormat_init;
179 set_print_func<numeric, ExFormat>(ExFormat::print_numeric);
189 FCFormat::_init::_init() {
192 FCFormat::_init FCFormat::FCFormat_init;
194 FormFormat::_init::_init() {
197 FormFormat::_init FormFormat::FormFormat_init;
217 QGRAF::Process::_init::_init() {
220 auto qbar =
Symbol(
"qbar");
222 auto lbar =
Symbol(
"lbar");
224 auto ghbar =
Symbol(
"ghbar");
227 auto Qbar =
Symbol(
"Qbar");
229 auto Cbar =
Symbol(
"Cbar");
231 auto Bbar =
Symbol(
"Bbar");
233 auto nbar =
Symbol(
"nbar");
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";
256 VerTeX[lst{Qbar, e, nbar}] =
"[crossed dot]";
257 VerTeX[lst{nbar, e, g}] =
"[crossed dot]";
285 Symbol ghWmbar(
"ghWmbar");
287 Symbol ghWpbar(
"ghWpbar");
302 Symbol ntaubar(
"ntaubar");
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++) {
314 LineTeX[fs.op(i)] =
"fermion ,edge label="+si;
315 LineTeX[fs2.op(i)] =
"anti fermion ,edge label="+si;
318 LineTeX[Wp] =
"photon, edge label=$W$";
319 LineTeX[Wm] =
"photon, edge label=$W$";
320 LineTeX[Z] =
"photon, edge label=$Z$";
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$";
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$"};
331 for(
int i=0; i<ghs.nops(); i++) {
333 LineTeX[ghs.op(i)] =
"ghost ,edge label="+si;
334 LineTeX[ghs2.op(i)] =
"anti ghost ,edge label="+si;
338 QGRAF::Process::_init QGRAF::Process::Process_init;
348 int QCD::FF::cur_mode = 0;
354 int DEX::Threads = 8;
364 _global_init::_init::_init() {
369 auto opath = getenv(
"PATH");
370 if(opath != NULL) oss <<
":" << opath;
371 setenv(
"PATH", oss.str().c_str(),
true);
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(); });
389 cln::cl_inhibit_floating_point_underflow =
true;
393 _global_init::_init _global_init::init_object;
401((<sign><symmetry_factor>)*
402<in_loop>InField(<field>,<field_index>,<momentum>)*
404<back><out_loop>OutField(<field>,<field_index>,<momentum>)*
406<back><propagator_loop>Propagator(
407<back>Field(<field>,<field_index>),
408<back>Field(<dual-field>,<dual-field_index>),
411<back><vertex_loop>Vertex(
412<back><ray_loop>Field(<field>,<field_index>,<momentum>),
423[ model = 'Gluon FF Model' ]
430[qbar, q, g; QCD='+1']
431[Qbar, Q, g; QCD='+1']
432[g, g, g, g; QCD='+2']
434[ghbar, gh, g; QCD='+1']
435[nbar, n, g; QCD='+1']
436[nbar, e, g; QCD='+0']
441[ model = 'Quark FF Model' ]
444[e, ebar, -, external]
448[qbar, q, g; QCD='+1']
449[Qbar, Q, g; QCD='+1']
450[g, g, g, g; QCD='+2']
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']
461% From Boehm, Denner, Joos @ Feynman Gauge
462[ model = 'Standard Model' ]
464%------------------------------------
466%------------------------------------
476% gluon and its ghost:
480% V: EWSM gauge bosons
493% G: Faddeev-Popov Ghosts
504%------------------------------------
506%------------------------------------
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']
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']
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']
527[Wp, Wm, A; QCD='0', QED='+1']
528[Wp, Wm, Z; QCD='0', QED='+1']
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']
539[H, H, H; QCD='0', QED='+1']
540[H, chi, chi; QCD='0', QED='+1']
541[H, phim, phip; QCD='0', QED='+1']
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']
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']
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']
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']
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']
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']
588[ep, em, A; QCD='0', QED='+1']
589[mup, mum, A; QCD='0', QED='+1']
590[taup, taum, A; QCD='0', QED='+1']
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']
599[ep, em, Z; QCD='0', QED='+1']
600[mup, mum, Z; QCD='0', QED='+1']
601[taup, taum, Z; QCD='0', QED='+1']
603[nebar, ne, Z; QCD='0', QED='+1']
604[nmubar, nmu, Z; QCD='0', QED='+1']
605[ntaubar, ntau, Z; QCD='0', QED='+1']
607[Dbar, U, Wm; QCD='0', QED='+1']
608[Dbar, C, Wm; QCD='0', QED='+1']
609[Dbar, T, Wm; QCD='0', QED='+1']
611[Sbar, U, Wm; QCD='0', QED='+1']
612[Sbar, C, Wm; QCD='0', QED='+1']
613[Sbar, T, Wm; QCD='0', QED='+1']
615[Bbar, U, Wm; QCD='0', QED='+1']
616[Bbar, C, Wm; QCD='0', QED='+1']
617[Bbar, T, Wm; QCD='0', QED='+1']
619[Ubar, D, Wp; QCD='0', QED='+1']
620[Ubar, S, Wp; QCD='0', QED='+1']
621[Ubar, B, Wp; QCD='0', QED='+1']
623[Cbar, D, Wp; QCD='0', QED='+1']
624[Cbar, S, Wp; QCD='0', QED='+1']
625[Cbar, B, Wp; QCD='0', QED='+1']
627[Tbar, D, Wp; QCD='0', QED='+1']
628[Tbar, S, Wp; QCD='0', QED='+1']
629[Tbar, B, Wp; QCD='0', QED='+1']
631[ep, ne, Wm; QCD='0', QED='+1']
632[mup, nmu, Wm; QCD='0', QED='+1']
633[taup, ntau, Wm; QCD='0', QED='+1']
635[nebar, em, Wp; QCD='0', QED='+1']
636[nmubar, mum, Wp; QCD='0', QED='+1']
637[ntaubar, taum, Wp; QCD='0', QED='+1']
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']
647[ep, em, H; QCD='0', QED='+1']
648[mup, mum, H; QCD='0', QED='+1']
649[taup, taum, H; QCD='0', QED='+1']
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']
658[ep, em, chi; QCD='0', QED='+1']
659[mup, mum, chi; QCD='0', QED='+1']
660[taup, taum, chi; QCD='0', QED='+1']
662[Dbar, U, phim; QCD='0', QED='+1']
663[Dbar, C, phim; QCD='0', QED='+1']
664[Dbar, T, phim; QCD='0', QED='+1']
666[Sbar, U, phim; QCD='0', QED='+1']
667[Sbar, C, phim; QCD='0', QED='+1']
668[Sbar, T, phim; QCD='0', QED='+1']
670[Bbar, U, phim; QCD='0', QED='+1']
671[Bbar, C, phim; QCD='0', QED='+1']
672[Bbar, T, phim; QCD='0', QED='+1']
674[Ubar, D, phip; QCD='0', QED='+1']
675[Ubar, S, phip; QCD='0', QED='+1']
676[Ubar, B, phip; QCD='0', QED='+1']
678[Cbar, D, phip; QCD='0', QED='+1']
679[Cbar, S, phip; QCD='0', QED='+1']
680[Cbar, B, phip; QCD='0', QED='+1']
682[Tbar, D, phip; QCD='0', QED='+1']
683[Tbar, S, phip; QCD='0', QED='+1']
684[Tbar, B, phip; QCD='0', QED='+1']
686[ep, ne, phim; QCD='0', QED='+1']
687[mup, nmu, phim; QCD='0', QED='+1']
688[taup, ntau, phim; QCD='0', QED='+1']
690[nebar, em, phip; QCD='0', QED='+1']
691[nmubar, mum, phip; QCD='0', QED='+1']
692[ntaubar, taum, phip; QCD='0', QED='+1']
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']
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']
724[ model = 'qcd Model' ]
725%------------------------------------
727%------------------------------------
731% gluon and its ghost:
736%------------------------------------
738%------------------------------------
740[qbar, q, g; QCD='+1']
741[Qbar, Q, g; QCD='+1']
742[g, g, g, g; QCD='+2']
744[ghbar, gh, g; QCD='+1']
746[qbar, q, A; QCD='+0']
747[Qbar, Q, A; QCD='+0']
class for AntiSymmetric Gamma object
void form_print(const FormFormat &c, unsigned level=0) const
void fc_print(const FCFormat &c, unsigned level=0) const
void print(const print_dflt &c, unsigned level=0) const
class for Dirac Gamma object
void form_print(const FormFormat &c, unsigned level=0) const
void print(const print_dflt &c, unsigned level=0) const
void fc_print(const FCFormat &c, unsigned level=0) const
class for Levi-Civita object https://onlinelibrary.wiley.com/doi/pdf/10.1002/9783527630097....
void fc_print(const FCFormat &c, unsigned level=0) const
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,...
void print(const print_dflt &c, unsigned level=0) const
void print(const print_context &c, unsigned level=0) const
static GiNaC::exmap Dimension
void form_print(const FormFormat &c, unsigned level=0) const
FormFormat print function.
void print(const print_dflt &c, unsigned level=0) const
default print function
void fc_print(const FCFormat &c, unsigned level=0) const
FCFormat print function.
void form_print(const FormFormat &c, unsigned level=0) const
print the Form Format
void fc_print(const FCFormat &c, unsigned level=0) const
void print(const print_dflt &c, unsigned level=0) const
normal priint
void print(const print_dflt &c, unsigned level=0) const
void form_print(const FormFormat &c, unsigned level=0) const
void fc_print(const FCFormat &c, unsigned level=0) const
void form_print(const FormFormat &c, unsigned level=0) const
void fc_print(const FCFormat &c, unsigned level=0) const
void print(const print_dflt &c, unsigned level=0) const
class extended to GiNaC symbol class, represent a positive symbol
static std::map< std::string, ex > Table
void print(const print_context &c, unsigned level=0) const
XIntegral Class, preface to SecDec.
void print(const print_dflt &c, unsigned level=0) const
class extended to GiNaC symbol class, pure imaginary symbol
static std::map< std::string, ex > Table
map< ex, string, ex_is_less > LineTeX
map< ex, string, ex_is_less > VerTeX
map< ex, string, ex_is_less > InOutTeX
int GiNaC_Parallel_Process
const int form_expand_none
map< string, bool > GiNaC_Parallel_RM
map< ex, long long, ex_is_less > fermat_weight
map< string, int > GiNaC_Parallel_NB
ex GAS(const ex &expr, unsigned rl)
function similar to GAD/GSD in FeynClac
const int form_trace_each_all
map< string, int > GiNaC_Parallel_NP
const int form_expand_all
map< string, int > GiNaC_Parallel_Verb
map< string, bool > GiNaC_Parallel_ReWR
const int o_normal_fermat
const int form_trace_each_each
const int form_trace_auto
bool dir_exists(string dir)
void set_precision(long prec, bool push)
map< string, string > GiNaC_Parallel_PRE
std::stack< cln::float_format_t > cln_prec_stack
std::stack< long > digits_stack