On photon statistics parametrized by a non-central Wishart random matrix

In order to tackle parameter estimation of photocounting distributions, polykays of acting intensities are proposed as a new tool for computing photon statistics. As unbiased estimators of cumulants, polykays are computationally feasible thanks to a symbolic method recently developed in dealing with sequences of moments. This method includes the so-called method of moments for random matrices and results to be particularly suited to deal with convolutions or random summations of random vectors. The overall photocounting effect on a deterministic number of pixels is introduced. A random number of pixels is also considered. The role played by spectral statistics of random matrices is highlighted in approximating the overall photocounting distribution when acting intensities are modeled by a non-central Wishart random matrix. Generalized complete Bell polynomials are used in order to compute joint moments and joint cumulants of multivariate photocounters. Multivariate polykays can be successfully employed in order to approximate the multivariate Poisson–Mandel transform. Open problems are addressed at the end of the paper.