Based on a perfect crystal component by: Miguel A. Gonzalez, A. Dianoux June 2013 (ILL)
- BUGFIX: correct neutron energy shift in Doppler mode
- added option 'debyescherrer' to select analyzer geometry
- added option 'facette' to approximate analyzer sphere by small, flat crystals
- inital release
Author: Markus Appel
Origin: ILL / FAU Erlangen-Nuernberg
Perfect crystal component, primarily for use as monochromator and analyzer in
backscattering spectrometers. Reflection from a single Bragg reflex of a flat or
spherical surface is simulated using a Darwin, Ewald or Gaussian profile. Doppler
movement of the monochromator is supported as well.
Orientation of the crystal surface is *different* from other monochromator components!
Gravitational energy shift for tall analyzers should work in principle, but it not tested yet.
See the parameter description on how to define the geometry and properties.
 Website for Backscattering Spectroscopy: http://www.ill.eu/sites/BS-review/index.htm
IN16B (ILL) Si111 large angle analyzers (approximate dimensions):
PerfectCrystal(radius=2, lambda=6.2708, sigma=0.24e-3,
ttmin=40, ttmax=165, phimin=-45, phimax=+45, centerfocus=1)
IN16B (ILL) Si111 Doppler monochromator:
PerfectCrystal(radius=2.165, lambda=6.2708, sigma=0.24e-3,
width = 0.500, height = 0.250, centerfocus=0,
speed = 4.7, amplitude = 0.075, exclusive=1)
WARNING: This is a contributed Component.
Parameters in boldface are required;
the others are optional.
analyzer coverage angle in horizontal xz-plane between -180 and 180
angular center position in the horizontal plane (only if centerfocus==1)
horizontal coverage as center and full width
vertical analyzer coverage between -90 and 90 (-180 and 180 if debyescherrer==1)
angular vertical position (only if centerfocus==1)
vertical coverage as center and full width
absolute height (only if debyescherrer==0)
bend analyzer following a Debye-Scherrer ring along scattering angle tt (twotheta)
(phi = [-180...180] in that case). Default: 0
width of square crystal facettes arranged on the spherical surface (set 0 to disable). Default: 0
random misalignemt of each facette. Default: 0
Component origin is the center of the analyzer sphere if 1,
if set to 0 the origin is on the analyzer surface. Default: 0
Radius of curvature, set to 0 for a flat surface. Default: 0
Scattering vector of the reflex (sometimes also called Q...)
Alternatively to tau: backscattered wavelength
Peak reflectivity. Default: 1
Plateau width of the Ewald/Darwin curve
randomized for each trajectory, full width
Relative variation of tau
Use Ewald curve if 1, Darwin curve if 0. Default: 1 (Ewald)
Width of the Gaussian reflectivity curve in meV (corresponding to
the energy resolution).
(The width will be transformed and the Gaussian is actually calculated in k-space.)
Simply reflect all neutrons if 1. Good for debugging/visualization. Default: 0
Maximum Doppler speed. The actual monochromator velocity is randomized between
-speed and +speed. Default: 0
Amplitude of the Doppler movement. Default: 0
Optimize Doppler phase for better MC efficiency if set to 1.
*WARNING:* Experimental option! Always compare to a simulation
without smartphase. Better do not use smartwidth with
Ewald/Darwin curves due to their endless tails. Default: 0
Half width of the possible energy reflection window for smartphase.
If set to 1, absorb all neutrons that missed the monochromator/analyzer surface.
Monte-Carlo probability of transmitting an event through the
monochromator/analyzer surface. (Events with R=1.0 for DarwinE/Ewald curves
are always reflected!). Default: 0