This paper describes a fast, robust and efficient method to calculate the scattered field due to multiple PEC cylinders located at random positions excited by a plane wave. The cylinder collocation follows stochastic geometry principles while including multiple electromagnetic interactions. By using this method, we study the statistics of scattered field due to Poisson point process distributed multiple cylinder samples.