From 98cf20cf0653abed5c11508c60b6279caeeacd74 Mon Sep 17 00:00:00 2001 From: Felix Erben Date: Tue, 9 Jul 2019 17:42:36 +0100 Subject: [PATCH] continued work on baryons --- Hadrons/Modules/MContraction/Nucleon.hpp | 15 +++++++++++---- 1 file changed, 11 insertions(+), 4 deletions(-) diff --git a/Hadrons/Modules/MContraction/Nucleon.hpp b/Hadrons/Modules/MContraction/Nucleon.hpp index 41e1348d..2caa8208 100644 --- a/Hadrons/Modules/MContraction/Nucleon.hpp +++ b/Hadrons/Modules/MContraction/Nucleon.hpp @@ -127,11 +127,12 @@ void TNucleon::execute(void) auto &q2 = envGet(PropagatorField2, par().q2); auto &q3 = envGet(PropagatorField3, par().q2); envGetTmp(LatticeComplex, c); - //envGetTmp(LatticeComplex, quark2); - //envGetTmp(LatticeComplex, quark3); envGetTmp(LatticeComplex, diquark); Result result; - + int nt = env().getDim(Tp); + result.corr.resize(nt); + + std::vector buf; // C = i gamma_2 gamma_4 => C gamma_5 = - i gamma_1 gamma_3 Gamma Cg5(Gamma::Algebra::SigmaXZ); Gamma g4(Gamma::Algebra::GammaT); //needed for parity P_\pm = 0.5*(1 \pm \gamma_4) @@ -179,7 +180,13 @@ void TNucleon::execute(void) } } - // saveResult(par().output, "meson", result); + /* sliceSum(c,buf,Tp); + for (unsigned int t = 0; t < buf.size(); ++t) + { + result.corr[t] = TensorRemove(buf[t]); + }*/ + + // saveResult(par().output, "baryon", result); } END_MODULE_NAMESPACE