Implements the operation $$ L = a*( \Gamma_{ik} M_{ki} ) $$ to give the trace of a quark loop \(M\) multiplied by a gamma structure \(\Gamma\) and scaled by a complex factor \(a\).
loop_spin_project(loop, gamma, reim = "both", stochav = FALSE,
scale_factor = as.complex(1), herm_conj = FALSE)
'raw_cf' container with loop data
4x4 complex matrix
String, one of 'real', 'imag' or 'both'. After the spin projection and trace, the result can be restricted to just the real or imaginary part, if desired. Useful for the cases in which it is clear that only one or the other contains any signal.
Boolean, specifies whether the average over stochastic samples should be performed. This makes the projection much faster but of course prevents the projected loop data to be used for the construction of diagrams with multiple quark loops.
Complex scaling factor to be applied.
Boolean, optionally the loop matrix \(M\) can be hermitian conjugated before the spin projection is performed.
Returns an object of class raw_cf.