Herbert Lee, Bruno Sanso, Weining Zhou and David Higdon

12/31/2003 09:00 AM

Applied Mathematics & Statistics

A beam of protons is produced by a linear charged particle accelerator, then focused through the use of successive quadrupoles.

The initial state of the beam is unknown, in terms of particle position and momentum. Wire scans are used to collect data on the current state of the beam as it passes through and beyond the focusing region, and the goal is to infer the initial state from the wire trace data. This setup is that of an inverse problem, in which a computer simulator is used to link an initial state configuration to observable values (wire traces), and then inference is performed for the distribution of the initial state. We present a Bayesian approach for solving this inverse problem.