A new mathematical approach for the calculation of the photon detectio
n statistics in single-molecule detection experiments is presented. Th
e approach is based on a path integral representation of the involved
probability distributions. A numerical Monte Carlo sampling method for
computing the involved path integrals is given. The general ideas are
applied to the modeling of measured burst size distributions in singl
e-molecule detection experiments with B-phycoerythrin.