Spherical detector

This shape of detector represents a portion of a sphere, defined by the range of $\varphi$ and $\alpha$ angles as shown in the plot below. The sphere has its center located at the origin of the sample coordinate system.

A GISAS or off-specular simulation creates a spherical detector by default. To set its parameters the method setDetectorParameters should be used

setDetectorParameters(n_phi, phi_min, phi_max, n_alpha, alpha_min, alpha_max)
"""
Sets detector parameters using angle ranges
 
n_phi     : number of phi-axis bins
phi_min   : low edge of first phi-bin
phi_max   : upper edge of last phi-bin
n_alpha   : number of alpha-axis bins
alpha_min : low edge of first alpha-bin
alpha_max : upper edge of last alpha-bin
"""

All angles are expressed in radians. For convenience, a special deg multiplier can be used to define the detector range directly in degrees. For example, the following code snippet

simulation = ScatteringSimulation()
simulation.setDetectorParameters(20, -1.0*deg, 1.0*deg, 10, 0.0*deg, 1.0*deg)

will create a detector plane with a total number of bins equal to 200 and with a fixed bin size equal to 0.1*deg along both directions

Here, the vertical and horizontal lines denote bin boundaries while blue dots show the bin centers. During the simulation, the bin intensity will be calculated for values of $\varphi_f$ and $\alpha_f$ corresponding to the bin centers and then normalized to the bin area.

Alternative way to create SphericalDetector

The detector can be initialized separately and then later assigned to the simulation object as shown in the script below. The results will be exactly the same as in the previous example

detector = SphericalDetector(20, -1.0*deg, 1.0*deg, 10, 0.0*deg, 1.0*deg)
simulation = ScatteringSimulation()
simulation.setDetector(detector)

Spherical detector with variable bin size

It is possible to create a spherical detector with a non-equidistant binning. This is done by creating custom axes and using them as SphericalDetector initializers.

There are two basic axis types in BornAgain: EquiDivision and GenericScale. Their signatures are shown in the code snippet below.

EquiDivision(name, nbins, start, end)
"""
Constructs an axis with fixed bin size
 
name  : The name of the axis
nbins : Number of axis bins
start : Low edge of first bin
end   : Upper edge of last bin
"""
 
GenericScale(name, bin_boundaries)
"""
Constructs an axis with variable bin size
 
name           : The name of the axis
bin_boundaries : Array of even size containing lower and upper edges for each bin.
"""

In the following example we create a SphericalDetector with two axes: an x-axis with variable bin size and a y-axis with fixed bin size. The code and resulting detector plane are shown below.

xedges = [-1.0*deg, -0.7*deg, -0.7*deg, -0.5*deg, -0.5*deg, -0.3*deg, -0.3*deg, -0.2*deg, -0.2*deg, -0.1*deg, -0.1*deg, -0.025*deg, -0.025*deg, 0.025*deg, 0.025*deg, 0.1*deg, 0.1*deg, 0.2*deg, 0.2*deg, 0.3*deg, 0.3*deg, 0.5*deg, 0.5*deg, 0.7*deg, 0.7*deg, 1.0*deg]
xaxis = GenericScale("my-x-axis", xedges)
 
yaxis = EquiDivision("my-y-axis", 10, 0.0, 10.)
 
detector = SphericalDetector()
detector.setDetectorAxes(xaxis, yaxis)
 
simulation = ScatteringSimulation()
simulation.setDetector(detector)