McStas logo McStas - A neutron ray-trace simulation package ILL; PSI; ESS Niels Bohr Institute DTU Physics NEXMAP


About McStas
 Conditions of use
 Project funding

 Other Downloads (share)

Mailing list

Search web/mailinglist

 Wiki (GitHub)
 McStas manual


Report bugs (GitHub)

Code-repository (GitHub)

McStas: PSD_Detector

[ Identification | Description | Input parameters | Links ]

The PSD_Detector Component

Position-sensitive gas-filled detector with gaseous thermal-neutron converter (box, cylinder or 'banana').


  • Site:
  • Author: Thorwald van Vuure
  • Origin: ILL
  • Date: Feb 21st, 2005


An n times m pixel Position Sensitive Detector (PSD), box, cylinder or banana
filled with a mixture of thermal-neutron converter gas and stopping gas.
This model implements wires detection, including charge drift, parallax, etc.
The default is a simple 1D/2D position detector. However, setting
the parameter 'type' to "events" will save the signal as an event file which
contains all neutron parameters at detection points, including time. This is
especially suited for TOF detectors. Please use Histogrammer afterwards.
The gas creation event (neutron absorption in gas) is flagged as SCATTERED==1,
and the charge detection (signal) is flagged as SCATTERED==2.

A flat rectangular box of given depth is simulated, or alternatively
a cylinder (when giving radius), or a horizontal cylindrical 'banana'
detector (when giving the width along the arc).
Box position is at its center, the 'banana' and cylinder have their
position at focus point, symetric in angular range for the 'banana'.

The converter gas can be specified as He-3, BF3 or N2.
The filling pressure of converter gas must be specified along
with the zdepth: this gives an absorption probability (e.g. ~90% for
0.18 nm neutrons at zdepth 3 cm and He-3 pressure 5 bar).
The stopping gas can be specified as one of the following gases:
CF4, C3H8, CO2, Xe/TMA or He, or any for which a table of stopping
powers is available. NB: each stopping gas table has specific
data on a pair of particles (from thermal-neutron absorptions in
either He-3, BF3 or N2) in a specific gas. Unusual choices like N2
as a converter and He as a stopping gas probably require the
calculation of a new table.
Tip: to do a simulation in pure He-3, specify 0 bar for the stopping gas.
The pressures of the two gases together put a lower limit on the
achievable spatial resolution, no matter what pixel size is
Example: 1 bar CF4 gives ~2.5 mm spatial resolution, 3 bar CF4 gives
~1 mm. C3H8 requires more pressure for the same resolution, then
Xe/TMA, then CO2 and finally He.

Taken into account are the following effects:
- (rectangularly distributed) error in the recorded position due to
the range of the neutron capture products in the gas
- angular dependence of absorption efficiency
- border effect (events outside the sensitive volume tend to be
counted in the border pixels)
- parallax effect (with, if desired, an electrostatic lens to
partially correct for this - use "LensOn=1")
- wall effect
To correctly take the wall effect into account, a value must be
specified for the energy threshold for acceptance of an event as a
valid neutron event. This normally serves to discriminate from
gammas. Because these are not simulated in McStas, a value of 0 keV
can be chosen for this threshold if one is interested in seeing the
maximum number of neutrons; however, a more realistic value would be
100 keV. This implies seeing the maximum number of neutrons as well,
except in geometries with dominating wall effect.
The border effect is a considerably complex phenomenon: it depends
on the precise electric field configuration near the border. It is
simulated here simply by introducing a dead zone, gas-filled,
bordering the sensitive volume. Events in there can still cause
signals in the detector. The dead zone can cause a shadow in the
sensitive volume, just as in reality.
The algorithm simply adds all events on the first 'pixel' in the
dead zone to the border pixels. This is a coarse approximation,
because the physical effect depends on the orientation of the
tracks and the energy threshold setting. The bottom line is that
tracking the position of the Center Of Gravity of the charge
deposition in the detector does not allow for accurate modeling of
the border effect, and therefore the border pixels should be
ignored, just as in a real detector.
Note that specifying 0 width of the dead zone next to the
sensitive volume, apart from being unrealistic, does not allow
accurate simulation of the border effect: in this case, there will
be a relative lack of events on the border pixels due to the wall
This component determines the position of the Center Of Gravity of
the charge cloud in the detector. It does not simulate independent
Manufacturing imperfections such as wire diameter variations and
spread in electronic amplifier component properties are not simulated.
This should allow an experimenter to identify these manufacturing
imperfections in his physical machine and correct for them.

Example: PSD_Detector(xwidth=0.192, yheight=0.192, nx=64, ny=64,
zdepth=0.03, threshold=100, borderx=-1, bordery=-1,
PressureConv=5, PressureStop=1,
FN_Conv="Gas_tables/He3inHe.table", FN_Stop="Gas_tables/He3inCF4.table",
xChDivRelSigma=0, yChDivRelSigma=0,

Example: PSD_Detector(xwidth=0.26, yheight=0.26, nx=128, ny=128,
zdepth=0.03, threshold=100, borderx=-1, bordery=-1,
PressureConv=5, PressureStop=1,
FN_Conv="Gas_tables/He3inHe.table", FN_Stop="Gas_tables/He3inCF4.table",
xChDivRelSigma=0, yChDivRelSigma=0,

Example: PSD_Detector(radius=0.0127, yheight=0.20, nx=1, ny=128,
threshold=100, borderx=-1,
PressureConv=9.9, PressureStop=0.1,
FN_Conv="Gas_tables/He3inHe.table", FN_Stop="Gas_tables/He3inCO2.table",

Example: PSD_Detector(awidth=1.533, yheight=0.4, nx=640, ny=256,
radius=0.73, zdepth=0.032, dc=0.026, threshold=100,
borderx=-1, bordery=-1,
PressureConv=5, PressureStop=1,
FN_Conv="Gas_tables/He3inHe.table", FN_Stop="Gas_tables/He3inCF4.table",
xChDivRelSigma=0, yChDivRelSigma=0.0037,
LensOn=1, filename="D19.psd")

Input parameters

Parameters in boldface are required; the others are optional.
nx1Number of pixel columns128
ny1Number of pixel rows128
xwidthmWidth of detector opening, in the case of a flat box0
radiusmRadius of detector for cylindrical/banana shape0
awidthmWidth of detector opening along the horizontal arc of the detector, measured along the inside arc of the sensitive volume in the case of a 'banana' detector0
yheightmHeight of detector opening0
zdepthmDepth of the sensitive volume of the detector0.03
thresholdkeVEnergy threshold for acceptance of an event100
PressureConvbarPressure of gaseous thermal-neutron converter (e.g. He-3)5
PressureStopbarPressure of stopping gas1
interpolate1Performs energy interpolation if true1
p_interact1Fraction of statistical events to be treated by component. Set it to 0 to use real absorption probability0
verbose1Gives more information and spectrum0
LensOn1set to 1 to simulate an electrostatic lens according to P. Van Esch, NIM A 540 (2005) pp. 361-367.1
dcmDistance from the entrance window to the cathode plane, where the correction zone of the lens ends for banana shape and LensOn=10
borderxmWidth of (gas-filled) border on the side of the detector-2
borderymHeight of (gas-filled) border on the top/bottom of the detector giving rise to various border effects. Set to -2 to obtain the height of one pixel, -1 for the depth of the detector (default), or specify a non-negativedistance-2
xChDivRelSigma1Sigma of the error in position determination along the x axis (relative to xwidth), uniquely due to charge division readout. Set to 0 in case of independent channel readout along the x dimension0
yChDivRelSigma1Relative sigma due to charge division in y direction, relative to yheight.0
bufsize1Size of neutron output buffer default is 0, i.e. save all - recommended.0
restore_neutron1If set, the monitor does not influence the neutron state0
angledegAngular opening of detector in the case of a 'banana' detector, then overrides awidth parameter0
typestringIf set to 'events' detector signal saves signal into an event file to be read e.g. by Histogrammer. This is to be used for TOF detectors.0
filenametextName of file in which to store the detector image. The name of the component is used if file name is not specified0
FN_ConvtextFilename of the stopping power table of the converter gas"He3inHe.table"
FN_StoptextFilename of the stopping power table of the stopping gas"He3inCF4.table"


[ Identification | Description | Input parameters | Links ]

Generated on 2024-01-03 11:47:43

Last Modified: Tuesday, 16-Apr-2024 20:05:21 CEST
Search website mailinglist archive GitHub repos