A new mathematical approach for the calculation of the photon detection statistics in single-molecule detection experiments is presented. The 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 single-molecule detection experiments with B-phycoerythrin.