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
76
77
78
79
80
81
82
83
|
/*******************************************************************************
*
* McStas, neutron ray-tracing package
* Copyright (C) 1997-2008, All rights reserved
* Risoe National Laboratory, Roskilde, Denmark
* Institut Laue Langevin, Grenoble, France
*
* Component: Beamstop
*
* %ID
*
* Written by: Peter Willendrup
* Date: November 2008
* Origin: Risoe
*
* Box-shaped absorbing slab.
*
* %D
* Infinitely absorbing slab of material.
*
* Example: Absorber(xmin=-0.01, xmax=0.01, ymin=-0.05, ymax=0.05, zmin=-0.2, zmax=0)
*
* %P
* INPUT PARAMETERS
*
* xmin: [m] Lower x bound
* xmax: [m] Upper x bound
* ymin: [m] Lower y bound
* ymax: [m] Upper y bound
* zmin: [m] Lower z bound
* zmax: [m] Upper z bound
*
* %END
*******************************************************************************/
DEFINE COMPONENT Absorber
SETTING PARAMETERS (xmin=0, xmax=0, ymin=0, ymax=0, zmin=0, zmax=0)
DECLARE
%{
double xw;
double yh;
double zt;
double xm;
double ym;
double zm;
%}
INITIALIZE
%{
xw = xmax-xmin;
yh = ymax-ymin;
zt = zmax-zmin;
xm = (xmax+xmin)/2;
ym = (ymax+ymin)/2;
zm = (zmax+zmin)/2;
if (xw == 0 || yh == 0 || zt == 0)
{ fprintf(stderr,"Absorber: %s: Error: Slab volume is zero!\n", NAME_CURRENT_COMP); exit(-1); }
%}
TRACE
%{
double t0,t1,xp,yp,zp;
xp = x - xm;
yp = y - ym;
zp = z - zm;
if (box_intersect(&t0, &t1, xp, yp, zp, vx, vy, vz, xw, yh, zt)) {
if (t1>=0) {
PROP_DT((t1+t0)/2);
SCATTER;
ABSORB;
}
}
%}
MCDISPLAY
%{
box(xm, ym, zm, xw, yh, zt,0, 0, 1, 0);
%}
END
|
