Correlation Function Manager

Name pending.

Chris has set up a workflow for computing the finite volume spectrum of an n to m meson system(at least he thinks so). The outline is as follows.

1. Determine the basis of operators you want to include by looking at the non-interacting energy levels. Note all of the operators that appear below the inelastic threshold, then pick at least the next 1 or as many as you wish(this is for use in the variational basis to remove excited state effects.

2. Project the operators to the relevant finite volume symmetry group, $O_h$ for cubic and $D_{4h}$ for elongated boxes.

3. Take the operators and compute all of the relevant wick contractions, ex. $\bar{u}u\bar{d}d -> u\bar{u}d\bar{d}$.

4. The wick contractions give a list of diagrams for each element of the correlation matrix that needs to be computed.

5. Once all relevant diagrams are computed, combine them into correlators, apply the GEVP to the correlation matrix and fit the resultant spectrum to obtain the energy levels.