Monte Carlo calculation of angular transmission functions for small angle scattering

Abstract

We define an angular transmission function η in the center of mass system. The convolution of the differential cross section σ with η yields the signal in the laboratory system. For the case of elastic small angle scattering by spherically symmetric potentials we calculate η by a Monte Carlo method. Random positions are taken in the beam defining collimators, resulting in a trajectory with a deflection angle at the scattering centre. These deflection angles are transformed to the c.m. system with the small angle tranformation formulae. From the distribution we calculate η as a histogram and the central moments of η. The function η depends on the velocity ratio and on the mass ratio of the scattering partners. We store the results in such a way that the central moments can be calculated afterwards for all mass and velocity ratios. By using the central moments the convolution integral can be reduced to a simple weighted sum of σ-values at different scattering angles. The r.m.s. deviations of the central moments scale with N 1 2 , with N the number of Monte Carlo trajectories. A typical deviation is 1% in the second order moment for N = 2 × 10 4, increasing slightly with increasing order of the moments. This method of calculation gives a large degree of freedom for optimisation of the collimation geometry. The use of an angular transmission function defined in the center of mass system gives a good insight in the experimental reflection of the physical events. As an example we apply the method to the case of small angle scattering of Ar as a primary beam by Kr as a secondary beam and the inverse configuration of Kr by Ar.