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 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 |