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.27.tar.gz (29.4 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.27-py3-none-any.whl (11.8 MB view details)

Uploaded Python 3

File details

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

File metadata

  • Download URL: TDCRPy-0.0.27.tar.gz
  • Upload date:
  • Size: 29.4 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.27.tar.gz
Algorithm Hash digest
SHA256 463073c8bc4982c425888840d1044b25303f11e0bbb19191dfac737d51c53c30
MD5 d68fb8237bd861a53bc0638c3701e6e6
BLAKE2b-256 9015b493f849c8cf8d947478f64f48cb49028bb41c847a60d91c21ddfbd800fb

See more details on using hashes here.

File details

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

File metadata

  • Download URL: TDCRPy-0.0.27-py3-none-any.whl
  • Upload date:
  • Size: 11.8 MB
  • 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.27-py3-none-any.whl
Algorithm Hash digest
SHA256 80199507229576f439dd9ce7d56cc519fbc443cc10415a6d915cdb21d665ed94
MD5 578322a60ddd3fd4a0cc75cb7f88aa0f
BLAKE2b-256 30609e8969b5b21f775908c871f0d4c90777757523c663efba034c05b879ad7d

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