portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

Hadrons: some updates in WeakMesonDecayKl2

Browse Source
This commit is contained in:
Vera Guelpers 2019-03-06 15:27:59 +00:00
parent 0b426bf9f6
commit d3935ae7fc
1 changed files with 13 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
Hadrons/Modules/MContraction/WeakMesonDecayKl2.hpp
View File

@ -34,7 +34,6 @@ See the full license in the file "LICENSE" in the top level distribution directo
#include <Hadrons/Global.hpp> #include <Hadrons/Global.hpp>
#include <Hadrons/Module.hpp> #include <Hadrons/Module.hpp>
#include <Hadrons/ModuleFactory.hpp> #include <Hadrons/ModuleFactory.hpp>
#include <Hadrons/Modules/MContraction/WeakHamiltonianNonEye.hpp>
BEGIN_HADRONS_NAMESPACE BEGIN_HADRONS_NAMESPACE
@ -84,14 +83,14 @@ class TWeakMesonDecayKl2: public Module<WeakMesonDecayKl2Par>
{ {
public: public:
FERM_TYPE_ALIASES(FImpl,); FERM_TYPE_ALIASES(FImpl,);
class Result: Serializable class Metadata: Serializable
{ {
public: public:
GRID_SERIALIZABLE_CLASS_MEMBERS(Result, GRID_SERIALIZABLE_CLASS_MEMBERS(Metadata,
int, spinidx1, int, spinidx1,
int, spinidx2, int, spinidx2);
std::vector<Complex>, corr);
}; };
typedef Correlator<Metadata> Result;
public: public:
// constructor // constructor
TWeakMesonDecayKl2(const std::string name); TWeakMesonDecayKl2(const std::string name);
@ -107,7 +106,7 @@ protected:
virtual void execute(void); virtual void execute(void);
}; };
MODULE_REGISTER_TMP(WeakMesonDecayKl2, ARG(TWeakMesonDecayKl2<FIMPL>), MContraction); MODULE_REGISTER_TMP(WeakMesonDecayKl2, TWeakMesonDecayKl2<FIMPL>, MContraction);
/****************************************************************************** /******************************************************************************
* TWeakMesonDecayKl2 implementation * * TWeakMesonDecayKl2 implementation *
@ -167,18 +166,13 @@ void TWeakMesonDecayKl2<FImpl>::execute(void)
envGetTmp(PropagatorField, prop_buf); envGetTmp(PropagatorField, prop_buf);
std::vector<Result> result; std::vector<Result> result;
result.resize(Ns*Ns); Result r;
for (unsigned int i = 0; i < Ns*Ns ; ++i)
{
result[i].corr.resize(nt);
}
for (unsigned int mu = 0; mu < 4; ++mu) for (unsigned int mu = 0; mu < 4; ++mu)
{ {
c = zero; c = zero;
//hadronic part: trace(q1*adj(q2)*g5*gL[mu]) //hadronic part: trace(q1*adj(q2)*g5*gL[mu])
c = trace(MAKE_CW_SUBDIAG(q1, q2, GammaL(Gamma::gmu[mu]))); c = trace(q1*adj(q2)*g5*GammaL(Gamma::gmu[mu]));
prop_buf = 1.; prop_buf = 1.;
//multiply lepton part //multiply lepton part
res += c * prop_buf * GammaL(Gamma::gmu[mu]) * lepton; res += c * prop_buf * GammaL(Gamma::gmu[mu]) * lepton;
@ -193,12 +187,15 @@ void TWeakMesonDecayKl2<FImpl>::execute(void)
sliceSum(buf, res_summed, Tp); sliceSum(buf, res_summed, Tp);
r.corr.clear();
for (unsigned int t = 0; t < nt; ++t) for (unsigned int t = 0; t < nt; ++t)
{ {
result[i].corr[t] = TensorRemove(res_summed[t]); r.corr.push_back(TensorRemove(res_summed[t]));
} }
result[i].spinidx1 = s1;
result[i].spinidx2 = s2; r.info.spinidx1 = s1;
r.info.spinidx2 = s2;
result.push_back(r);
i+=1; i+=1;
} }