  We present a method to integrate predictions from a theoretical model of a
reaction with three bodies in the final state over the region of phase space
covered by a given experiment. The method takes into account the true
experimental acceptance, as well as variations of detector efficiency, and
eliminates the need for a Monte-Carlo simulation of the detector setup. The
method is applicable to kinematically complete experiments. Examples for the
use of this method include several polarization observables in dp break-up at
270 MeV. The calculations are carried out in the Faddeev framework with the CD
Bonn nucleon-nucleon interaction, with or without the inclusion of an
additional three-nucleon force.
