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

added parity operator to baryon fields

Browse Source
This commit is contained in:
ferben 2019-02-08 15:49:52 +00:00
parent d26a5dce12
commit c4d27ee30f
2 changed files with 34 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Hadrons/Modules/MDistil/BContraction.hpp
View File

@ -26,6 +26,7 @@ public:
std::string, two, std::string, two,
std::string, three, std::string, three,
std::string, output, std::string, output,
int, parity,
std::vector<std::string>, mom); std::vector<std::string>, mom);
}; };
@ -104,6 +105,8 @@ void TBContraction<FImpl>::execute(void)
int N_2 = two.size(); int N_2 = two.size();
int N_3 = three.size(); int N_3 = three.size();
int parity = par().parity;
LOG(Message) << "Computing distillation baryon fields" << std::endl; LOG(Message) << "Computing distillation baryon fields" << std::endl;
LOG(Message) << "One: '" << par().one << "' Two: '" << par().two << "' Three: '" << par().three << "'" << std::endl; LOG(Message) << "One: '" << par().one << "' Two: '" << par().two << "' Three: '" << par().three << "'" << std::endl;
LOG(Message) << "Momenta:" << std::endl; LOG(Message) << "Momenta:" << std::endl;
@ -138,9 +141,13 @@ void TBContraction<FImpl>::execute(void)
SpinVector tmp222; SpinVector tmp222;
SpinVector tmp111; SpinVector tmp111;
assert(parity == 1 || parity == -1);
std::vector<std::vector<int>> epsilon = {{0,1,2},{1,2,0},{2,0,1},{0,2,1},{2,1,0},{1,0,2}}; std::vector<std::vector<int>> epsilon = {{0,1,2},{1,2,0},{2,0,1},{0,2,1},{2,1,0},{1,0,2}};
std::vector<int> epsilon_sgn = {1,1,1,-1,-1,-1}; std::vector<int> epsilon_sgn = {1,1,1,-1,-1,-1};
Gamma g4(Gamma::Algebra::GammaT);
gamma12_ = { gamma12_ = {
Gamma::Algebra::Identity, // I Gamma::Algebra::Identity, // I
Gamma::Algebra::Gamma5, // gamma_5 Gamma::Algebra::Gamma5, // gamma_5
@ -227,8 +234,11 @@ void TBContraction<FImpl>::execute(void)
tmp33s()(is)() = tmp33[sU]()(is)(epsilon[ie][2]); tmp33s()(is)() = tmp33[sU]()(is)(epsilon[ie][2]);
} }
tmp333 = Gamma(gamma23_[0])*tmp33s; tmp333 = Gamma(gamma23_[0])*tmp33s;
tmp111 = Gamma(gamma12_[0])*tmp11s;
tmp222 = g4*tmp111;
tmp111 = 0.5*(double)parity*(tmp111 + tmp222); // P_\pm * ...
diquark2 = factor23[0]*innerProduct(tmp22s,tmp333); diquark2 = factor23[0]*innerProduct(tmp22s,tmp333);
BField2[Bindex]+=(double)epsilon_sgn[ie]*tmp11s*diquark2; BField2[Bindex]+=(double)epsilon_sgn[ie]*tmp111*diquark2;
} }
} }
} }

28
tests/hadrons/Test_hadrons_distil.cc
View File

@ -225,19 +225,36 @@ void test_MesonField(Application &application)
application.createModule<MContraction::A2AMesonField>("DistilMesonField",A2AMesonFieldPar); application.createModule<MContraction::A2AMesonField>("DistilMesonField",A2AMesonFieldPar);
} }
///////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////
// BaryonFields // BaryonFields - phiphiphi
///////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////
void test_BaryonField(Application &application) void test_BaryonFieldPhi(Application &application)
{ {
// DistilVectors parameters // DistilVectors parameters
MDistil::BContraction::Par BContractionPar; MDistil::BContraction::Par BContractionPar;
BContractionPar.one="DistilVecs_phi"; BContractionPar.one="DistilVecs_phi";
BContractionPar.two="DistilVecs_phi"; BContractionPar.two="DistilVecs_phi";
BContractionPar.three="DistilVecs_phi"; BContractionPar.three="DistilVecs_phi";
BContractionPar.output="BaryonField"; BContractionPar.output="BaryonFieldPhi";
BContractionPar.parity=1;
BContractionPar.mom={"0 0 0"}; BContractionPar.mom={"0 0 0"};
application.createModule<MDistil::BContraction>("BaryonField",BContractionPar); application.createModule<MDistil::BContraction>("BaryonFieldPhi",BContractionPar);
}
/////////////////////////////////////////////////////////////
// BaryonFields - rhorhorho
/////////////////////////////////////////////////////////////
void test_BaryonFieldRho(Application &application)
{
// DistilVectors parameters
MDistil::BContraction::Par BContractionPar;
BContractionPar.one="DistilVecs_rho";
BContractionPar.two="DistilVecs_rho";
BContractionPar.three="DistilVecs_rho";
BContractionPar.output="BaryonFieldRho";
BContractionPar.parity=1;
BContractionPar.mom={"0 0 0"};
application.createModule<MDistil::BContraction>("BaryonFieldRho",BContractionPar);
} }
bool bNumber( int &ri, const char * & pstr, bool bGobbleWhiteSpace = true ) bool bNumber( int &ri, const char * & pstr, bool bGobbleWhiteSpace = true )
@ -565,7 +582,8 @@ int main(int argc, char *argv[])
test_LapEvec( application ); test_LapEvec( application );
test_Perambulators( application ); test_Perambulators( application );
test_DistilVectors( application ); test_DistilVectors( application );
test_BaryonField( application ); test_BaryonFieldPhi( application );
test_BaryonFieldRho( application );
break; break;
} }
LOG(Message) << "====== XML creation for test " << iTestNum << " complete ======" << std::endl; LOG(Message) << "====== XML creation for test " << iTestNum << " complete ======" << std::endl;