Cylindrical holes in layered system

Simulate the next sample, representing unordered cylindrical holes in two magnetic layers:

Parameters:

• cylindrical holes (air): radius: 5 nm, height: 20 nm
• particle density: 0.002
• thickness layer 1 and 1: 10 nm
• layer 1 and 2: $\delta = 6\cdot 10^{-6}$, $\beta = 2\cdot 10^{-8}$
• layer 1 magnetization: (0, 1e5, 0)
• layer 2 magnetizaiton: (0, -1e5, 0)
• substrate material: $\delta = 3\cdot 10^{-6}$, $\beta = 2\cdot 10^{-8}$
• beam: wavelength = $1$ Å, $\alpha = 0.2$ degree, intensity $1\cdot 10^8$
• detector (200x200): $\alpha$ between 0 and 2 degree, $\phi$ between -1 and 1 degree
• polarization analysis: spin flip from +z to -z