can produce same result as original work, might be wrong
put every thing related to diquark momentum in vec_diquark.h