def simulate_particles_for_density(p):
        """Get number of particles in detector given a particle density

        :param p: particle density in number per detector area.
        :return: random number from Poisson distribution.

        """
        return random.poisson(p)