Skip to main content

TDCR model

Project description

TDCRPy

TDCRPy is a Python code to calculate detection efficiency of a liquide scintillation counter using 3-photomultiplier tubes. The calculation is based on the photo-physical model called of the Triple-to-Double-Coincidence-Ratio method (TDCR) [1] and a Monte-Carlo sampling allowing to adress complexe decay schemes and radionuclide mixtures. The process is summarized in the figure below.

drawing

Nuclear decay

The code directly reads decay data from the Decay Data Evaluation Project (DDEP) web interface [2] that is recommanded to be used by the radionuclide metrology community. The PenNuc format [3] is used to simulate decays and the $\beta$ spectra from the BetaShape code [4] are used. The BetaShape code estimates accurate $\beta$ spectra by taking the atomic exchange effect and also simulate accurately electron capture decay [5]. It has been demonstrated to offer significant improvement in the context of liquid scintillation counting [6].

Atomic relaxation

The atomic relaxation from missing electrons in the inner-shell following electron capture and internal conversion is simulated by ENSDF data on the DDEP web interface.

Interaction

The interaction of $\gamma$ rays, electrons and positrons are simulated using response kernels calculated by the Monte-Carlo code MCNP6 developped by Los Alamos [13].

Scintillation

The stopping power of electrons between 20 keV and 1000 keV is a mixture of a radiative loss model [7] and a collision model [8] that has been validated agaisnt the NIST model ESTAR [9] recommanded by the ICRU Report 37 [10]. At low energy - between 10 eV and 20 keV - the model from Tan and Xia [11] is implemented.

The stopping power of $\alpha$ particles of energy comprises between 1 keV and 8 MeV comes from the NIST code ASTAR [9] recommanded in the ICRU Report 49 [12]. For energy below 1 keV, an extrapolation is made.

Statistical model

...

References

[1] Ryszard Broda, Krzysztof Pochwalski, Tomasz Radoszewski, Calculation of liquid-scintillation detector efficiency, Applied Radiation and Isotopes 39:2, 1988, 159-164, https://doi.org/10.1016/0883-2889(88)90161-X

[2] http://www.lnhb.fr/ddep_wg/

[3] E. García-Toraño, V. Peyres, F. Salvat, PenNuc: Monte Carlo simulation of the decay of radionuclides, Computer Physics Communications 245, 2019, 106849 https://doi.org/10.1016/j.cpc.2019.08.002

[4] X. Mougeot, Erratum: Reliability of usual assumptions in the calculation of $\beta$ and $\bar{\mu}$ spectra, Physical Review C 91, 2015, 055504, https://doi.org/10.1103/PhysRevC.92.059902

[5] X. Mougeot, Towards high-precision calculation of electron capture decays, Applied Radiation and Isotopes 154, 2019, 108884, https://doi.org/10.1016/j.apradiso.2019.108884

[6] K. Kossert, X. Mougeot, Improved activity standardization of 90Sr/90Y by means of liquid scintillation counting, Applied Radiation and Isotopes 168, 2021, 109478, https://doi.org/10.1016/j.apradiso.2020.109478

[7] S.M. Seltzer, M.R. Berger, M. R., Evaluation of the collision stopping power of elements and compounds for electrons and positrons, Applied Radiation and Isotopes 33:11, 1982, 1189-1218, https://doi.org/10.1016/0020-708x(82)90244-7

[8] M.O. El-Ghossain, Calculations Of Stopping Power, And Range Of Electrons Interaction With Different Material And Human Body Parts, International Journal of Scientific & Technology Research 6:1 2017. https://www.ijstr.org/final-print/jan2017/Calculations-Of-Stopping-Power-And-Range-Of-Electrons-Interaction-With-Different-Material-And-Human-Body-Parts.pdf

[9] M.J. Berger, J.S. Coursey, M.A. Zucker and J. Chang,Stopping-Power & Range Tables for Electrons, Protons, and Helium Ions, NIST Standard Reference Database 124, 2017, https://dx.doi.org/10.18434/T4NC7P

[10] ICRU Report 37, Stopping Powers for Electrons and Positrons

[11] Z. Tan, Y. Xia, Stopping power and mean free path for low-energy electrons in ten scintillators over energy range of 20–20,000 eV, Applied Radiation and Isotopes 70, 2012, 296-300, https://doi.org/10.1016/j.apradiso.2011.08.012

[12] ICRU Report 49, Stopping Power and Ranges for Protons and Alpha Particles, https://www.icru.org/report/stopping-power-and-ranges-for-protons-and-alpha-particles-report-49/

[13] https://mcnp.lanl.gov/

Project details


Release history Release notifications | RSS feed

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

TDCRPy-0.0.13.tar.gz (28.9 kB view details)

Uploaded Source

Built Distribution

If you're not sure about the file name format, learn more about wheel file names.

TDCRPy-0.0.13-py3-none-any.whl (32.4 kB view details)

Uploaded Python 3

File details

Details for the file TDCRPy-0.0.13.tar.gz.

File metadata

  • Download URL: TDCRPy-0.0.13.tar.gz
  • Upload date:
  • Size: 28.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.7.16

File hashes

Hashes for TDCRPy-0.0.13.tar.gz
Algorithm Hash digest
SHA256 e8dfaae8d26bc4a20f15e76e3537ec35c9c3dbb5f6f8b230c4949e9eea5162e4
MD5 12d144e58c73504646b462137bfb9059
BLAKE2b-256 8e803e6e28442d86dcfccef08878766e3ae30af724a15f335762393fead78e1d

See more details on using hashes here.

File details

Details for the file TDCRPy-0.0.13-py3-none-any.whl.

File metadata

  • Download URL: TDCRPy-0.0.13-py3-none-any.whl
  • Upload date:
  • Size: 32.4 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.7.16

File hashes

Hashes for TDCRPy-0.0.13-py3-none-any.whl
Algorithm Hash digest
SHA256 c6a78768280e5ab43128a8ad89c695407f6e31e5c31e91a9ae4ccef199282ed3
MD5 c8f5d3ab28c9320e02b0263120f1f645
BLAKE2b-256 35d2c3617eb315c47a4a94c7ec5cc565cd89b778429bdc17675ef88c6c31b15e

See more details on using hashes here.

Supported by

AWS Cloud computing and Security Sponsor Datadog Monitoring Depot Continuous Integration Fastly CDN Google Download Analytics Pingdom Monitoring Sentry Error logging StatusPage Status page