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
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
/*******************************************************************************
*
* McStas, neutron ray-tracing package
* Copyright 1997-2002, All rights reserved
* Risoe National Laboratory, Roskilde, Denmark
* Institut Laue Langevin, Grenoble, France
*
* %I
* Written by: Kim Lefmann
* Date: May 7, 2001
* Version: $Revision$
Origin: Risoe
* Release: McStas 1.6
*
* Rectangular 1D PSD, measuring intensity vs. vertical position, x
*
* %D
*
* Example: PSDlin_monitor(nx=20, filename="Output.x", xmin=-0.1, xmax=0.1, ymin=-0.1, ymax=0.1)
*
* %P
* INPUT PARAMETERS:
*
* xmin: Lower x bound of detector opening [m]
* xmax: Upper x bound of detector opening [m]
* ymin: Lower y bound of detector opening [m]
* ymax: Upper y bound of detector opening [m]
* xwidth: Width of detector. Overrides xmin, xmax [m]
* yheight: Height of detector. Overrides ymin, ymax [m]
* nx: Number of x bins [1]
* filename: Name of file in which to store the detector image [string]
* restore_neutron: If set, the monitor does not influence the neutron state [1]
* nowritefile: [1] If set, monitor will skip writing to disk
*
* OUTPUT PARAMETERS:
*
* PSDlin_N: Array of neutron counts
* PSDlin_p: Array of neutron weight counts
* PSDlin_p2: Array of second moments
*
* %E
******************************************************************************/
DEFINE COMPONENT PSDlin_monitor
DEFINITION PARAMETERS (nx=20)
SETTING PARAMETERS (string filename=0, xmin=-0.05, xmax=0.05, ymin=-0.05, ymax=0.05,
xwidth=0, yheight=0, restore_neutron=0, int nowritefile=0)
OUTPUT PARAMETERS (PSDlin_N, PSDlin_p, PSDlin_p2)
/* Neutron parameters: (x,y,z,vx,vy,vz,t,sx,sy,sz,p) */
DECLARE
%{
double PSDlin_N[nx];
double PSDlin_p[nx];
double PSDlin_p2[nx];
%}
INITIALIZE
%{
int i;
if (xwidth > 0) { xmax = xwidth/2; xmin = -xmax; }
if (yheight > 0) { ymax = yheight/2; ymin = -ymax; }
if ((xmin >= xmax) || (ymin >= ymax)) {
printf("PSDlin_monitor: %s: Null detection area !\n"
"ERROR (xwidth,yheight,xmin,xmax,ymin,ymax). Exiting",
NAME_CURRENT_COMP);
exit(0);
}
for (i=0; i<nx; i++)
{
PSDlin_N[i] = 0;
PSDlin_p[i] = 0;
PSDlin_p2[i] = 0;
}
%}
TRACE
%{
int i;
PROP_Z0;
if (x>xmin && x<xmax && y>ymin && y<ymax)
{
i = floor(nx*(x-xmin)/(xmax-xmin)); /* Bin number */
if((i >= nx) || (i<0))
{
printf("FATAL ERROR: wrong positioning in linear PSD. i= %i \n",i);
exit(1);
}
PSDlin_N[i]++;
PSDlin_p[i] += p;
PSDlin_p2[i] += p*p;
}
if (restore_neutron) {
RESTORE_NEUTRON(INDEX_CURRENT_COMP, x, y, z, vx, vy, vz, t, sx, sy, sz, p);
}
%}
SAVE
%{
if (!nowritefile) {
DETECTOR_OUT_1D(
"Linear PSD monitor",
"Test",
"TEST I",
"x", xmin, xmax, nx,
&PSDlin_N[0],&PSDlin_p[0],&PSDlin_p2[0],
filename);
}
%}
MCDISPLAY
%{
multiline(5, (double)xmin, (double)ymin, 0.0,
(double)xmax, (double)ymin, 0.0,
(double)xmax, (double)ymax, 0.0,
(double)xmin, (double)ymax, 0.0,
(double)xmin, (double)ymin, 0.0);
%}
END
|
