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
|
/*******************************************************************************
*
* McStas, neutron ray-tracing package
* Copyright (C) 1997-2008, All rights reserved
* Risoe National Laboratory, Roskilde, Denmark
* Institut Laue Langevin, Grenoble, France
*
* Instrument: Test_absorption
*
* %Identification
* Written by: Mads Bertelsen
* Date: September 2015
* Origin: University of Copenhagen
* %INSTRUMENT_SITE: Tests_union
*
* Test of Union material with only absorption
*
* %Description
* Testing an absorber made using the Union components. An absorber is made
* without processes and just the make_material component, and requires
* a few special cases.
*
* %Example: Detector: screen_I=5.91451e-06
*
* %End
*******************************************************************************/
DEFINE INSTRUMENT Test_absorption()
DECLARE
%{
%}
TRACE
COMPONENT init = Union_init()
AT (0,0,0) ABSOLUTE
COMPONENT absorber_material = Union_make_material(my_absorption=3,absorber=1)
AT (0,0,0) ABSOLUTE
COMPONENT a1 = Progress_bar()
AT (0,0,0) ABSOLUTE
COMPONENT source = Source_div(
xwidth=0.15, yheight=0.15,focus_aw=0.01, focus_ah=0.01,
E0 = 10,
dE = 0)
AT (0,0,0) RELATIVE a1
ROTATED (0,0,0) RELATIVE a1
COMPONENT beam_center = Arm()
AT (0,0,1.4) RELATIVE a1
ROTATED (0,0,0) RELATIVE a1
COMPONENT cylinder = Union_cylinder(
radius = 0.05, yheight = 0.1,
priority=6, material_string="absorber_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 screen = PSD_monitor(xwidth=0.2,yheight=0.2,nx=200,ny=200,filename="absoprtion_picture.dat",restore_neutron=1)
AT (0,0,2) RELATIVE a1
COMPONENT stop = Union_stop()
AT (0,0,0) ABSOLUTE
END
|
