The LAMP Application

Generated by IDLdoc

single page | use frames     summary     class     fields     routine details     file attributes

lamp_mac/IN13/

elascan_cyl.pro


******************************************************************************* Evaluation of correction coefficients in the case of a cylindrical sample. Both absorption and self shielding are taken into account. Self-shielding factors for each region are computed using a Montecarlo method by adapting the routines of R.M.Dimeo and John Copley. Call w6=elascan_cyl (sample_w=1, background_w=2, Rin=10.0, Rout=11.0, t_sample=0.90, t_container=0.99,$ thick1=0.1, thick2=0.1, twice=0.5, f=1.0, nmc=10000, corr_w=20) The corrected data will then be stored in w6. Additionally an ASCII file called 'correction_factors.dat' is written. This file contains the self-shielding factors Assc, Acsc, and Acc applied in the correction. Another ASCII file named 'backthrough.dat' containing the attenuation factor applied to the neutrons passing a second time through the sample is also written. If the background workspace contains more than one temperature measurement, the subtracted background is that measured at the closest T at which the sample has been measured. A file 'temperatures.dat' is also written. It contains a list of the sample temperatures and the corresponding background temperature used in the subraction. MAG (04/5/2009) *******************************************************************************

Routines

insect, nann, rann, x0, y0, u, v, dism, disp, nint

ninsect, nann, rann, x0, y0, nphi, u, v, dism, disp, nint

nout_calc, nann, dism, disp, sigma, nphi, arg_out

calc_area, rann, d, area

result = elascan_cyl(w_data, background_w=background_w, Rin=Rin, Rout=Rout, t_sample=t_sample, t_container=t_container, thick1=thick1, thick2=thick2, twice=twice, f=f, nmc=nmc, corr_w=corr_w)

Routine details

top insect

insect, nann, rann, x0, y0, u, v, dism, disp, nint

** ****** ** Calculates the intersections of a directed line with a set of ** concentric circles centered at the origin. The number of circles ** is nann and their radii, in increasing order, are (rann(i),i=1,nann). ** The line goes through the point (x0,y0) and its direction is (u,v). ** The calculated distances to intersection with circle k are dism(k), ** disp(k), such that dism(k).le.disp(k). If there is no intersection ** both quantities are set equal to zero. The number of intersections is ** nint.

Parameters

nann
rann
x0
y0
u
v
dism
disp
nint

top ninsect

ninsect, nann, rann, x0, y0, nphi, u, v, dism, disp, nint

** ******* ** Calculates the intersections of nphi directed lines with a set of ** concentric circles centered at the origin. The number of circles ** is nann and their radii, in increasing order, are (rann(i),i=1,nann). ** The lines go through the point (x0,y0) and their directions are (u,v). ** The calculated distances to intersection with circle k are dism(k), ** disp(k), such that dism(k).le.disp(k). If there is no intersection ** both quantities are set equal to the same value. ** The numbers of intersections are nint.

Parameters

nann
rann
x0
y0
nphi
u
v
dism
disp
nint

top nout_calc

nout_calc, nann, dism, disp, sigma, nphi, arg_out

** ********* ** calculates the attenuation from scattering point to all detectors

Parameters

nann
dism
disp
sigma
nphi
arg_out

top calc_area

calc_area, rann, d, area

** ********* ** elascan_cyc stuff

Parameters

rann
d
area

top elascan_cyl

result = elascan_cyl(w_data, background_w=background_w, Rin=Rin, Rout=Rout, t_sample=t_sample, t_container=t_container, thick1=thick1, thick2=thick2, twice=twice, f=f, nmc=nmc, corr_w=corr_w)

******* *********** ** User callable

Parameters

w_data

Keywords

background_w
Rin
Rout
t_sample
t_container
thick1
thick2
twice
f
nmc
corr_w

File attributes

Modification date: Tue Oct 1 14:41:21 2013
Lines: 762