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
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
/*******************************************************************************
*
* McStas, neutron ray-tracing package
* Copyright (C) 1997-2008, All rights reserved
* Risoe National Laboratory, Roskilde, Denmark
* Institut Laue Langevin, Grenoble, France
*
* Instrument: Union_test_powder
*
* %Identification
* Written by: Mads Bertelsen
* Date: September 2015
* Origin: University of Copenhagen
* %INSTRUMENT_SITE: Tests_union
*
* Simple test instrument for powder process in Union framework.
*
* %Description
* Simple test instrument for powder sample simulated using Union components.
*
* %Example: Detector: m4pi_I=4.23557e-06
*
* %End
*******************************************************************************/
DEFINE INSTRUMENT Test_powder()
DECLARE
%{
%}
TRACE
COMPONENT init = Union_init()
AT (0,0,0) ABSOLUTE
COMPONENT test_incoherent = Incoherent_process(sigma=1.5,packing_factor=1,unit_cell_volume=66.4)
AT (0,0,0) ABSOLUTE
COMPONENT test_powder = Powder_process(reflections="Al.laz")
AT (0,0,0) ABSOLUTE
COMPONENT test_material = Union_make_material(my_absorption=4)
AT (0,0,0) ABSOLUTE
COMPONENT a1 = Progress_bar()
AT (0,0,0) ABSOLUTE
COMPONENT source = Source_div(
xwidth=0.0005, yheight=0.0005,focus_aw=3, focus_ah=2,
E0 = 100,
dE = 0)
AT (0,0,0) RELATIVE a1 ROTATED (0,0,0) RELATIVE a1
// Sample position
COMPONENT beam_center = Arm()
AT (0,0,1.4) RELATIVE a1
ROTATED (0,0,0) RELATIVE a1
// V1
COMPONENT test_powder_sample = Union_cylinder(radius=0.04,yheight=0.03,priority=5,material_string="test_material")
AT (0,0,0) RELATIVE beam_center
ROTATED (0,0,0) RELATIVE beam_center
COMPONENT test_sample = Union_master()
AT(0,0,0) RELATIVE beam_center
ROTATED(0,0,0) RELATIVE beam_center
COMPONENT m4pi = PSD_monitor_4PI(radius=1, nx=180, ny=180, filename="Events",restore_neutron=1)
AT (0, 0, 0) RELATIVE beam_center
COMPONENT stop = Union_stop()
AT (0,0,0) ABSOLUTE
END
|
