throughput¶
| author: | Stanislav Stoupin |
|---|---|
| email: | <sstoupin@gmail.com> |
calculate throughput of a multicrystal configuration (dynamical theory of x-ray diffraction for perfect crystals)
SYNOPSIS¶
throughput [options] func dpsi Ec dEx ne input_file
DESCRIPTION¶
A program to calculate throughput of a multicrystal configuration given by input_file using dynamical theory of x-ray diffraction for perfect crystals For a brief summary of options and parameters run:
throughput -h
PARAMETERS¶
| func: | angular divergence distribution function for incoming x-rays:
|
|||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| dpsi: | angular divergence in units of [urad]
|
|||||||||
| Ec: | central energy for energy distribution in units of [keV] (corrected automatically depending on the choice of source and presence of crystals in backscattering configuration in the input file) |
|||||||||
| dEx: | energy half-range in units of [meV] |
|||||||||
| ne: | number of steps in the energy grid to perform the calculation (max 1000) |
OPTIONS¶
| -v, –version: | show version of program. |
||||||||
|---|---|---|---|---|---|---|---|---|---|
| -h, –help: | show summary of options. |
||||||||
| -o F, –output=F: | |||||||||
write results to file F (default to stdout) |
|||||||||
| -w D, –write=D: | write data to file D (default - no action) |
||||||||
| -p, –pi: | \(\pi\) polarization for incident wave (default - \(\sigma\) polarization) |
||||||||
| -a NTH, –angular_scan=NTH: | |||||||||
perform crystal rotation (angular scan) with NTH points |
|||||||||
| -s SRC, –source=SRC: | |||||||||
type of energy distribution for the source:
|
|||||||||
EXAMPLES¶
This is an example input file for calculation of throughput of a multicrystal configuration
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | # Crystal input file for throughput program
# Author: Stanislav Stoupin <sstoupin@aps.anl.gov> 2014
#################################################################################################################
# Example: calculation of throughput of a double-crystal C (diamond) high-heat-load monochromator
# to perform calculation for incident beam with Gaussian divergence of 5 urad rms at 23.7 keV
# in the photon energy range from -5000 meV to 5000 meV with 250 steps run:
# throughput g 5 23.7 5000 250 thru_hhlmC.in
#################################################################################################################
# columns:
# element h k l eta[deg] phi[deg] T[K] dc[mm] sign(+counterclockwise,-clockwise) R/T entity exactbb(1true,0false)
# offset[urad] angular_scan(0,1,2) angular_increment[urad]
## CDDW #########################################################################################################
C 1 1 1 0 0 350 0.3 1 R 1 0 0 0 0
C 1 1 1 0 0 300 0.5 -1 R 2 0 0 0 0
|
This is an example input file for calculation of a rocking curve in a multicrystal configuration
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | # Crystal input file for throughput program
# Author: Stanislav Stoupin <sstoupin@aps.anl.gov> 2014
#################################################################################################################
# Example: calculation of rocking curve of the second crystal of a double-crystal C (diamond) high-heat-load monochromator
# to perform calculation for an incident beam with a Gaussian divergence of 5 urad rms at 23.7 keV
# in the photon energy range from -5000 meV to 5000 meV with 250 steps,
# with angular scan of 100 points (with angular increment being the last parameter in the second line)
# and save data into a file "c2rc_hhlmC_23keV.dat" run the following:
# throughput -a 100 -w c2rc_hhlmC_23keV.dat g 5 23.7 5000 250 c2rc_hhlmC.in
#################################################################################################################
# columns:
# element h k l eta[deg] phi[deg] T[K] dc[mm] sign(+counterclockwise,-clockwise) R/T entity exactbb(1true,0false)
# offset[urad] angular_scan(0,1,2) angular_increment[urad]
## CDDW #########################################################################################################
C 1 1 1 0 0 350 0.3 1 R 1 0 0 0 0
C 1 1 1 0 0 300 0.5 -1 R 2 0 0 1 0.5
|