add script to calculate background contamination

This commit adds a script to calculate the background contamination using a
method inspired by the bifurcated analysis method used in SNO. The method works
by looking at the distribution of several high level variables (radius, udotr,
psi, and reconstructed z position) for events tagged by the different data
cleaning cuts and assuming that any background events which sneak past the data
cleaning cuts will have a similar distribution (for certain backgrounds this is
assumed and for others I will actually test this assumption. For more details
see the unidoc). Then, by looking at the distribution of these high level
variables for all the untagged events we can use a maximum likelihood fit to
determine the residual contamination.

There are also a few other updates to the plot-energy script:

- add a --dc command line argument to plot corner plots for the high level
  variables used in the contamination analysis

- add a fudge factor to the Ockham factor of 100 per extra particle

- fix a bug by correctly setting the final kinetic energy to the sum of the
  individual kinetic energies instead of just the first particle

- fix calculation of prompt events by applying at the run level

0 files changed, 0 insertions, 0 deletions