A Coulombscattering and energyabsorbing element simulating material in the beam path.
Parallel capable? : yes
GPU capable? : yes
Backtracking capable? : no
Parameter Name  Units  Type  Default  Description 
L  M  double  0.0  length 
LEFFECTIVE  M  double  0.0  effective length (used if L=0) 
XO  M  double  0.0  radiation length 
ENERGY_DECAY  long  0  If nonzero, then particles will lose energy due to material using a simple exponential model. 

ENERGY_STRAGGLE  long  0  Use simpleminded energy straggling model coupled with ENERGY_DECAY=1? 

NUCLEAR_BREMSSTRAHLUNG  long  0  Model energy loss to nuclear bremsstrahlung? If enabled, set ENERGY_DECAY=0 to disable simpler model. 

ELECTRON_RECOIL  long  0  If nonzero, electron recoil during Coulomb scattering is included (results in energy change). 

Z  long  0  Atomic number 

A  AMU  double  0.0  Atomic mass 
RHO  KG∕M^{3}  double  0.0  Density 
PRESSURE  PASCAL  double  0.0  Pressure. Used with temperature and atomic mass to compute density for ideal gas. 
TEMPERATURE  K  double  0.0  Temperature. Used with pressure and atomic mass to compute density for ideal gas. 
MULTIPLICITY  long  1  Atoms per gas molecule. 

PLIMIT  double  0.05  Probability cutoff for each slice 

WIDTH  M  double  0.0  Full width of slots. If 0, no slots are present. 
SPACING  M  double  0.0  Centertocenter spacing of slots. If 0, no slots are present. 
MATTER continued
A Coulombscattering and energyabsorbing element simulating material in the beam path.
Parameter Name  Units  Type  Default  Description 
TILT  RAD  double  0.0  Tilt of slot array about the longitudinal axis. 
CENTER  M  double  0.0  Position of center of slot array in rotated frame. 
N_SLOTS  long  0  Number of empty slots in material. If <=0, an infinite array is assumed. 

START_PASS  long  1  If nonnegative, pass on which to start interaction with beam. 

END_PASS  long  1  If nonnegative, pass on which to end interaction with beam. 

GROUP  string  NULL  Optionally used to assign an element to a group, with a userdefined name. Group names will appear in the parameter output file in the column ElementGroup 

This element is based on section 3.3.1 of the Handbook of Accelerator Physics and Engineering, specifically, the subsections Single Coulomb scattering of spin particles, Multiple Coulomb scattering through small angles, and Radiation length. There are two aspects to this element: scattering and energy loss.
Scattering. The multiple Coulomb scattering formula is used whenever the thickness of the material is greater than 0.001X_{o}, where X_{o} is the radiation length. (Note that this is inaccurate for materials thicker than 100X_{o}.) For this regime, the user need only specify the material thickness (L) and the radiation length (XO).
For materials thinner than 0.001X_{o}, the user must specify additional parameters, namely, the atomic number (Z), atomic mass (A), and mass density (RHO) of the material. Note that the density is given in units of kg∕m^{3}. (Multiply by 10^{3} to convert g∕cm^{3} to kg∕m^{3}.) In addition, the simulation parameter PLIMIT may be modified.
To understand this parameter, one must understand how elegant simulates the thin materials. First, it computes the expected number of scattering events per particle, E = σ_{T }nL = , where n is the number density of the material, L is the thickness of the material, K_{1} = ()^{2}, and K_{2} = , with r_{e} the classical electron radius and α the fine structure constant. The material is then broken into N slices, where N = E∕P_{limit}. For each slice, each simulation particle has a probability E∕N of scattering. If scattering occurs, the location within the slice is computed using a uniform distribution over the slice thickness.
For each scatter that occurs, the scattering angle, θ is computed using the cumulative probability distribution F(θ > θ_{o}) = . This can be solved for θ_{o}, giving θ_{o} = . For each scatter, F is chosen from a uniform random distribution on [0,1].
Energy loss. There are two ways to compute energy loss in materials, using a simple minded approach and using the bremsstrahlung cross section. The latter is recommended, but the former is kept for backward compatibility.
Slotted absorber. If the WIDTH and SPACING parameters are set to nonzero values, then a slotted absorber is simulated. The number of slots is by default infinite, but can be limited by setting N_SLOTS to a positive value; in this case, the slot array is centered about the transverse coordinate given by the CENTER parameter.
Note that the simulation contains a simplification in that particles cannot leave or enter the material through the side of the slot. I.e., if a particle is inside (outside) the material when it hits the front face of the object, it is assumed to remain inside (outside) until it has passed the object. For long objects, breaking the simulation up into multiple MATTER elements is suggested if a slotted arrangement is being simulated.
Onesided scrapers. One sided scrapers may be modeled using the SCRAPER element. It uses the same materialmodeling algorithm as described here.
MAXAMP