I am simulating a simple detector, imaging a very faint astronomical source (so that some dozens of (uncorrelated) photons are detected in every frame)
Currently I am first sampling from Poisson-distribution to determine how many photons are to arrive in a given frame.
Then for every photon I have to determine the pixel it will hit: I have for every pixel a number, that gives the probability, that if a photon comes, it will hit that pixel. (so that the sum over all pixels is 1.) I have tools to sample N hits from this PSF (the pixelwise hit probabilities are proportional to the Point Spread Function of the telescope)
It would however be more convenient, if I could simply multiply the PSF with the expected number of photons per frame, and sample for each pixel from this Poisson distribution.
Are the two methods equivalent? If not, which of them is more correct?