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
only if centerfocus==1
[deg] angular center position in the horizontal plane
horizontal coverage as center and full width
-180 and 180 if debyescherrer==1
[deg] vertical analyzer coverage between -90 and 90
only if centerfocus==1
[deg] angular vertical position
vertical coverage as center and full width
only if debyescherrer==0
[m] absolute height
(0/1) bend analyzer following a Debye-Scherrer ring along scattering angle tt (twotheta) (phi =
width of square crystal facettes arranged on the spherical surface (set 0 to disable). Default: 0 Warning: Facettation will fail at the poles along +-y axis.
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
sometimes also called Q...
[A^-1] Scattering vector of the reflex
Alternatively to tau: backscattered wavelength
Plateau width of the Ewald/Darwin curve (see
randomized for each trajectory, full width
 Relative variation of tau
[0/1] Use Ewald curve if 1, Darwin curve if 0. Default: 1
The width will be transformed and the Gaussian is actually calculated in k-space.
[meV] Width of the Gaussian reflectivity curve in meV (corresponding to the energy resolution).
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. Default: 5*sigma OR 10*2*dtauovertau*E0
If set to 1, absorb all neutrons that missed the monochromator/analyzer surface. Default: 0
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