1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
|
/*******************************************************************************
* McXtrace instrument definition URL=http://www.mcxtrace.org
*
* Instrument: Test Isotropic_Sqw output
*
* %Identification
* Written by: E. Farhi
* Date: Sept 21st 2022
* Origin: DTU
* %INSTRUMENT_SITE: Tests_samples
*
* Test output of Isotropic_Sqw on a spherical monitor.
*
* %Description
* A test instrument for testing Isotropic_Sqw output on a spherical monitor.
*
* %Example: E0=12 Detector: Sph_mon_I=3.63728e-21
*
* %Parameters
* E0: [keV] Energy emitted from source
* dE: [keV] Energy spread at the source
* L1: [m] Source-sample distance
* sqw_coh: [str] Sqw material definition
*
* %End
*******************************************************************************/
DEFINE INSTRUMENT Test_Sqw(E0=12.0, dE=1e-6, L1=10,
string sqw_coh="Rb_liq_coh.sqw")
TRACE
REMOVABLE COMPONENT Origin = Progress_bar()
AT (0,0,0) ABSOLUTE
/* source with constant flux */
REMOVABLE COMPONENT Source = Source_flat(
radius = 1e-7, dist = L1, focus_xw = 1e-7, focus_yh = 1e-7,
E0=E0, dE=dE)
AT (0, 0, 0) RELATIVE Origin
COMPONENT sample_cradle = Arm()
AT (0, 0, L1) RELATIVE PREVIOUS
/* The inelastic sample */
COMPONENT Sqw = Isotropic_Sqw(radius=0.0005, yheight=0.001, Sqw_coh=sqw_coh,
verbose=3, p_interact=0.95)
AT (0, 0, 0) RELATIVE sample_cradle
EXTEND %{
if (!SCATTERED) ABSORB;
%}
COMPONENT Sph_mon = PSD_monitor_4PI(nx=100,ny=100, radius=1, filename="Sphere")
AT (0, 0, 0) RELATIVE PREVIOUS
//COMPONENT Sqw_mon = Monitor_Sqw(index=-2, radius=2, yheight=2, bins=128)
// AT (0, 0, 0) RELATIVE PREVIOUS
END
|
