Abstract
Simulation of multiplicity distributions with the Monte Carlo method is difficult because each history is treated individually. In order to accurately model the multiplicity distribution, the intensity and time width of the interrogation pulse must be incorporated into the calculation. This behavior dictates how many photons arrive at the target essentially simultaneously. In order to model the pulse width correctly, a Monte Carlo code system consisting of modified versions of the codes MCNPX and MCNP-PoliMi has been developed in conjunction with a post-processing algorithm to operate on the MCNP-PoliMi output file. The purpose of this subroutine is to assemble the interactions into groups corresponding to the number of interactions which would occur during a given pulse. The resulting multiplicity distributions appear more realistic and capture the higher-order multiplets which are a product of multiple reactions occurring during a single accelerator pulse. Plans are underway to gather relevant experimental data to verify and validate the methodology developed and presented here. This capability will enable the simulation of a large number of materials and detector geometries. Analysis of this information will determine the feasibility of using multiplicity distributions as an identification tool for special nuclear material.