def prop_trace(prop_list, aq, sin=True):
"Trace of qslash*inv(prop)*(1/12q^2). Supports two forms of q."
N = len(prop_list)
prop_inv = inv(sum(prop_list)/N)
qslash = dw.slash(aq)
qslash_sin = dw.slash(np.sin(aq))
qsquared = dw.inner(aq)
qsquared_sin = dw.inner(np.sin(aq))
if sin:
return (1./12)*trace(dot(qslash_sin, prop_inv)).imag/qsquared_sin
else:
return (1./12)*trace(dot(qslash, prop_inv)).imag/qsquared
def populate_kinematic_variables(self):
'''Requires that L, T, and a are defined.'''
self.ap = dw.ap(self.p, self.tw, self.L, self.T)
self.ap2 = dw.ap(self.p2, self.tw, self.L, self.T)
self.aq = dw.aq(self.ap, self.ap2) # ap - ap2
self.apSq = dw.inner(self.ap) # (ap)^2
self.mu = dw.mu(self.ap, self.a)