subroutine esensmap Loop over image pixels: Loop over input images: IF detection mask == 1 Accumulate background in 3 x 3 pixel detection cell Determine exposure in detection cell END IF END Loop Evaluate cumulative Poissonian distribution Subtract background and divide by exposure Correct for fraction of PSF outside of the detection cell IF exposure in pixel /= 0: Write result to output sensitivity map END Loop end subroutine esensmap