TDCRPy 0.0.6__py3-none-any.whl → 0.0.8__py3-none-any.whl

TDCRPy/TDCRoptimize.py

@@ -0,0 +1,99 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Wed Jul  5 10:04:53 2023
+
+@author: romain.coulon
+"""
+
+import numpy as np
+import TDCRPy as td
+import scipy.optimize as opt
+import sys, time
+sys.path.insert(1, 'G:\Python_modules\BIPM_RI_PyModules')
+import TDCRcalculation as tc 
+
+
+
+def eff(TDCR_measure, TAB, TBC, TAC, Rad, kB):
+    
+    N=50000
+    L=1
+    # r=opt.minimize_scalar(td.TDCRPy, args=(TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "res", "sym"), method='golden')#,options={'xatol': 1e-7, 'disp': True, 'maxiter':100})
+    # r=opt.minimize_scalar(td.TDCRPy, args=(TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "res", "sym"), method='brent')
+    r=opt.minimize_scalar(td.TDCRPy, args=(TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "res", "sym"), method='bounded', bounds=[0.5, 2])
+    
+    # r=opt.minimize(td.TDCRPy, L, args=(TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "res", "sym"), method='nelder-mead',options={'xatol': 1e-7, 'disp': True, 'maxiter':100})
+    L=r.x
+    print(r)
+    #L=(L*0.995, L*1.021, L*0.988)           # Free paramete in keV-1
+    # L = (L, L, L)
+    # r=opt.minimize(td.TDCRPy, L, args=(TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "res", "asym"), method='nelder-mead',options={'xatol': 1e-7, 'disp': True, 'maxiter':100})
+    # L=r.x
+    # print(r)
+    out=td.TDCRPy(L,TDCR_measure, TAB, TBC, TAC, Rad, "1", N, kB, 0.98, 1000, "eff", "sym")
+    return np.mean(L), L, out[2], out[2], out[3]
+    
+def TicTocGenerator():
+    """
+    Generator that returns time differences
+    """
+    ti = 0           # initial time
+    tf = time.time() # final time
+    while True:
+        ti = tf
+        tf = time.time()
+        yield tf-ti # returns the time difference
+
+TicToc = TicTocGenerator() # create an instance of the TicTocGen generator
+
+# This will be the main function through which we define both tic() and toc()
+def toc(tempBool=True):
+    """
+    Prints the time difference yielded by generator instance TicToc
+    """
+    
+    tempTimeInterval = next(TicToc)
+    if tempBool:
+        print( "Elapsed time: %f seconds.\n" %tempTimeInterval )
+
+def tic():
+    """
+    Records a time in TicToc, marks the beginning of a time interval
+    """
+    toc(False)
+
+
+"""
+TEST
+"""
+
+Rad="Co-60"    # list of radionuclides (Na-24)
+kB =1.0e-5       # Birks constant in cm/keV  
+TDCR_measure = 0.977667386529166        # Measured TDCR value
+TAB = 0.992232838598821
+TBC = 0.992343419459002
+TAC = 0.99275350064608
+
+
+
+tic()
+F1, FFF, eff1, eff2 = tc.I2calc(TDCR_measure, TAB, TBC, TAC, Rad, kB)
+toc()
+
+print("/nanalytical model")
+print("free parameter = ", F1)
+print("free parameters = ", FFF)
+print("Double count rate efficiency (sym) = ", eff1)
+print("Double count rate efficiency (asym) = ", eff2)
+
+
+tic()
+F1, FFF, eff1, eff2, u = eff(TDCR_measure, TAB, TBC, TAC, Rad, kB)
+toc()
+
+print("/nstochastic model")
+print("free parameter = ", F1)
+print("free parameters = ", FFF)
+print("Double count rate efficiency (sym) = ", eff1)
+print("Double count rate efficiency (asym) = ", eff2)
+print("u(eff) = ", u)

TDCRPy/__init__.py

@@ -0,0 +1,3 @@
+from TDCRPy import TDCRPy
+from TDCRPy import TDCR_model_lib
+from TDCRPy import TDCRoptimize

TDCRPy/decay.py

@@ -0,0 +1,29 @@
+"""
+This program aims to calculate the relative composition of radionuclides
+decaying in chain using the Bateman model implemented in the radioactivedecay module form (MIT)
+
+Romain Coulon, Jialin Hu
+Bureau International des Poids et Mesures
+"""
+
+# pip install radioactivedecay
+import radioactivedecay as rd
+import matplotlib.pyplot as plt
+
+rad = "Ra-223" # Radionuclide
+A0 = 1         # Inital activity
+unit = 'Bq'    # unit for the activity
+decaytime = 2  # decay time
+unitt = 'd'    # unit for the time 
+
+rad_t0 = rd.Inventory({rad: A0}, unit)
+rad_t1 = rad_t0.decay(decaytime, unitt)
+rad_t1.activities(unit)
+print(rad_t1)
+
+decayGraph = rad_t0.plot(decaytime, unitt, yunits=unit)
+plt.savefig('decayGraph.png')
+
+nuc = rd.Nuclide(rad)
+decayChain = nuc.plot()
+plt.savefig('decayChain.png')

TDCRPy/test.py

@@ -0,0 +1,12 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Thu Jul 13 19:15:23 2023
+
+@author: romain.coulon
+"""
+import TDCRPy
+
+TDCRPy.TDCRPy(1, 0.5, 0.5, 0.5, 0.5, "H-3", "1", 10, 0.00001, 1,10, "eff", "sym", Display=True)
+
+
+TDCRPy.TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, RHO, nE, mode, mode2)

TDCRPy/test1.py

@@ -0,0 +1,87 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Mon May 15 15:41:41 2023
+
+@author: romain.coulon
+"""
+
+
+
+import zipfile
+import numpy as np
+
+file = 'decayData/All-nuclides_Ensdf.zip'
+z = zipfile.ZipFile(file)
+
+with z.open('Ag-108m.txt') as f:
+    data = f.readlines()
+    nL = np.size(data)
+    
+    ## fomatage en list de vecteur pour chaque line du fichier ENSDF
+    for i in range(nL):              # boucle dans data
+        data[i] = str(data[i])                  # conversion en string
+        data[i] = data[i].replace("b'",'')       # on enlève les caratères de formatage
+        data[i] = data[i].replace("\\r\\n'",'')  # on enlève les caratères de formatage
+    for i in range(nL):
+        data[i] = data[i].split()
+    
+    for i in range(nL):
+        if i>0 and ("L" in data[i]) and ("AUGER" in data[i]) and ("|]" in data[i-1]):
+            data.insert(i,[data[i][0], 'T'])
+
+    ## Définir des repères
+    index_auger = []          # indice de
+    daugther = []             # liste des filles
+    index_end = []            # indice dernier valeur pour une fille
+    posi = []                 # indice des blocs transitions 
+    for i, p in enumerate(data): # boucle dans le fichier
+        if 'DECAY' in p:
+            daugther.append(p[0])   # liste des filles
+        if 'Auger' in p:
+            index_auger.append(i)   # repère de ligne pour ensuite définir l'ensemble des blocs
+        if len(p)==2:
+            posi.append(i)          # repère de ligne pour ensuite définir les blocs (2)
+        # if 'L' in p
+        #     posi:
+        if 'P' in p:
+            index_end.append(i)     # repère de fin pour ensuite définir l'ensemble des blocs
+            posi.append(i)
+            
+    ## Filtrage des données utiles et formatage
+    energy = []
+    prob = []
+    Type = []
+    Fille = []
+    for i in range(len(posi)-1): # bloucle dans les blocs
+        start = posi[i]+1        # indice debut du bloc
+        end = posi[i+1]          # indice fin du bloc
+        d = data[start:end]      # Le bloc
+        e_b = []                   # Energie des transitions du bloc
+        prob_b = []                # Proba des transitions du bloc
+        Type_b = []                # Types des transitions du bloc
+        Fille_b = []               # Fille des transitions du bloc
+        #e2 = []
+        if start == end:
+            continue
+        if start-1 in index_end:
+            #e = []
+            continue
+        for n, p1 in enumerate(d):
+            #e2 = []
+            if '-' in p1[2]: # traitement des intervals de nombres réels
+                x = p1[2].split('-')
+                p1[2] = round((float(x[1])+float(x[0]))/2,3) # Moyenne
+            if '|]' in p1:        # traitment des accolades 
+                if len(p1)>6:     # repère de la ligne portant la proba pour le groupe 
+                    prob_b.append(float(p1[4]))
+                e_b.append(float(p1[2]))
+            else:
+                e_b.append(float(p1[2]))
+                prob_b.append(float(p1[3]))
+        
+        if len(prob_b)==1 and len(e_b)>1:
+            e_b = [np.mean(e_b)]
+        
+        print(e_b)
+        print(prob_b)
+        print(" ")

TDCRPy-0.0.8.dist-info/METADATA

@@ -0,0 +1,88 @@
+Metadata-Version: 2.1
+Name: TDCRPy
+Version: 0.0.8
+Summary: TDCR model
+Home-page: https://github.com/RomainCoulon/TDCRPy
+Author: RomainCoulon (Romain Coulon)
+Author-email: <romain.coulon@bipm.org>
+License: UNKNOWN
+Keywords: python,TDCR,Monte-Carlo,radionuclide,scintillation,counting
+Platform: UNKNOWN
+Classifier: Development Status :: 2 - Pre-Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Natural Language :: English
+Classifier: Natural Language :: French
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: Unix
+Classifier: Operating System :: MacOS :: MacOS X
+Classifier: Operating System :: Microsoft :: Windows
+Classifier: Topic :: Scientific/Engineering :: Physics
+Description-Content-Type: text/markdown
+License-File: LICENCE.md
+Requires-Dist: numpy
+Requires-Dist: scipy
+Requires-Dist: sys
+Requires-Dist: time
+Requires-Dist: urllib
+Requires-Dist: zipfile
+Requires-Dist: re
+
+# 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]](#1) and a Monte-Carlo sampling allowing to adress complexe decay schemes and radionuclide mixtures. The process is summarized in the figure below.
+
+<img src="./FlowChart.jpg" alt="drawing" width="500"/>
+
+## Nuclear decay 
+
+The code directly reads decay data from the Decay Data Evaluation Project (DDEP) web interface [[2]](#2) that is recommanded to be used by the radionuclide metrology community. The PenNuc format [[3]](#3) is used to simulate decays and the $\beta$ spectra from the BetaShape code [[4]](#4) are used. The BetaShape code estimates accurate $\beta$ spectra by taking the atomic exchange effect and also simulate accurately electron capture decay [[5]](#5). It has been demonstrated to offer significant improvement in the context of liquid scintillation counting [[6]](#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]](#13). 
+
+## Scintillation
+
+The stopping power of electrons between 20 keV and 1000 keV is a mixture of a radiative loss model [[7]](#7) and a collision model [[8]](#8) that has been validated agaisnt the NIST model ESTAR [[9]](#9) recommanded by the ICRU Report 37 [[10]](#10). At low energy - between 10 eV and 20 keV - the model from Tan and Xia [[11]](#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]](#9) recommanded in the ICRU Report 49 [[12]](#12). For energy below 1 keV, an extrapolation is made.
+
+## Statistical model
+
+...
+
+## References
+
+<a id="1">[1]</a> 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
+
+<b id="2">[2]</b> http://www.lnhb.fr/ddep_wg/
+
+<c id="3">[3]</c> 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
+
+<c id="4">[4]</c> 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
+
+<c id="5">[5]</c> 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
+
+<c id="6">[6]</c> K. Kossert, X. Mougeot, Improved activity standardization of <sup>90</sup>Sr/<sup>90</sup>Y by means of liquid scintillation counting, *Applied Radiation and Isotopes* **168**, 2021, 109478, https://doi.org/10.1016/j.apradiso.2020.109478
+
+<c id="7">[7]</c> 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
+
+<c id="8">[8]</c> 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
+
+<c id="9">[9]</c> 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
+
+<c id="10">[10]</c> ICRU Report 37, *Stopping Powers for Electrons and Positrons*
+
+<c id="11">[11]</c> 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
+
+<c id="12">[12]</c> 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/
+
+<c id="13">[13]</c> https://mcnp.lanl.gov/
+
+

TDCRPy-0.0.8.dist-info/RECORD

@@ -0,0 +1,23 @@
+Code/Activity_TDCR.py,sha256=wcllEOJbeEOexAYzcFGQx2luIBhxSU3LtpEalRmNziw,11530
+Code/EfficiencyProfils.py,sha256=-ZUPva1dU7lMRcKbQAewX9QyiTDwCiNQDoaQiw7dOI4,4829
+Code/TDCRPy.py,sha256=igN-xE3ixC27R0aaZHQ832wp8w5UV0iDjx_1cEc8Kvk,23214
+Code/TDCR_model_lib.py,sha256=bF2S2L3L5c4y2Bg6G1tG1PVR7mJJ1L6UU7G0IEHNsTY,79829
+Code/TDCRoptimize.py,sha256=2CssmP_SZSe6C_Bkgei3hon21RPytH4DetMbXi-e6wk,3171
+Code/__init__.py,sha256=IQBC9crikuuGlHfxIdP65CrgSHgl9U8z2jSYduNpWB0,93
+Code/decay.py,sha256=Lnq-kj-6Xs2Rw9w3Y3LlK_ffFip1QzemHdesiug0Wko,817
+Code/test.py,sha256=3r8LxXmQwlHNIefMswwUCy3OoHHKDVK8A3izerGT15U,292
+Code/test1.py,sha256=rGA9d1cgArhAfxRAAdtiz-vPGQ0Ns3HEdN8sCwydHYM,3250
+TDCRPy/Activity_TDCR.py,sha256=wcllEOJbeEOexAYzcFGQx2luIBhxSU3LtpEalRmNziw,11530
+TDCRPy/EfficiencyProfils.py,sha256=-ZUPva1dU7lMRcKbQAewX9QyiTDwCiNQDoaQiw7dOI4,4829
+TDCRPy/TDCRPy.py,sha256=lrsngHHgBxuOdOmUqxEEbMHjpHnSF8ZP-ZUknfPs37c,23207
+TDCRPy/TDCR_model_lib.py,sha256=bF2S2L3L5c4y2Bg6G1tG1PVR7mJJ1L6UU7G0IEHNsTY,79829
+TDCRPy/TDCRoptimize.py,sha256=2CssmP_SZSe6C_Bkgei3hon21RPytH4DetMbXi-e6wk,3171
+TDCRPy/__init__.py,sha256=IQBC9crikuuGlHfxIdP65CrgSHgl9U8z2jSYduNpWB0,93
+TDCRPy/decay.py,sha256=Lnq-kj-6Xs2Rw9w3Y3LlK_ffFip1QzemHdesiug0Wko,817
+TDCRPy/test.py,sha256=3r8LxXmQwlHNIefMswwUCy3OoHHKDVK8A3izerGT15U,292
+TDCRPy/test1.py,sha256=rGA9d1cgArhAfxRAAdtiz-vPGQ0Ns3HEdN8sCwydHYM,3250
+TDCRPy-0.0.8.dist-info/LICENCE.md,sha256=fuYzrZRiOAjJBzA1tsGQwojCgGROArb2Ec48GDTjlWM,1086
+TDCRPy-0.0.8.dist-info/METADATA,sha256=CU_LjClcOCPImFHxekh2IA_RxqPwuLem6ek03GO7j-I,5761
+TDCRPy-0.0.8.dist-info/WHEEL,sha256=2wepM1nk4DS4eFpYrW1TTqPcoGNfHhhO_i5m4cOimbo,92
+TDCRPy-0.0.8.dist-info/top_level.txt,sha256=5dSMF0KOskVmT0QGz9KnqcD9Ccc2Dkq0Sa8FOZk-7Ck,5
+TDCRPy-0.0.8.dist-info/RECORD,,

TDCRPy-0.0.6.dist-info/METADATA

@@ -1,31 +0,0 @@
-Metadata-Version: 2.1
-Name: TDCRPy
-Version: 0.0.6
-Summary: TDCR model
-Home-page: UNKNOWN
-Author: RomainCoulon (Romain Coulon)
-Author-email: <romain.coulon@bipm.org>
-License: UNKNOWN
-Keywords: python,TDCR,Monte-Carlo,radionuclide,scintillation,counting
-Platform: UNKNOWN
-Classifier: Development Status :: 2 - Pre-Alpha
-Classifier: Intended Audience :: Science/Research
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Natural Language :: English
-Classifier: Natural Language :: French
-Classifier: Programming Language :: Python :: 3
-Classifier: Operating System :: Unix
-Classifier: Operating System :: MacOS :: MacOS X
-Classifier: Operating System :: Microsoft :: Windows
-Classifier: Topic :: Scientific/Engineering :: Physics
-License-File: LICENCE.md
-Requires-Dist: numpy
-Requires-Dist: scipy
-Requires-Dist: sys
-Requires-Dist: time
-Requires-Dist: urllib
-Requires-Dist: zipfile
-Requires-Dist: re
-
-UNKNOWN
-

TDCRPy-0.0.6.dist-info/RECORD

@@ -1,14 +0,0 @@
-Code/Activity_TDCR.py,sha256=wcllEOJbeEOexAYzcFGQx2luIBhxSU3LtpEalRmNziw,11530
-Code/EfficiencyProfils.py,sha256=-ZUPva1dU7lMRcKbQAewX9QyiTDwCiNQDoaQiw7dOI4,4829
-Code/TDCRPy.py,sha256=ffOfmaqWdk5PWoVDep8zs_VoQLUhUCQePC0oAbPOYVs,23202
-Code/TDCR_model_lib.py,sha256=vl0pMsGG3FKAmfpnHfjUX8mokay1dvlH99GHYd9BnsI,79115
-Code/TDCRoptimize.py,sha256=2CssmP_SZSe6C_Bkgei3hon21RPytH4DetMbXi-e6wk,3171
-Code/__init__.py,sha256=IQBC9crikuuGlHfxIdP65CrgSHgl9U8z2jSYduNpWB0,93
-Code/decay.py,sha256=Lnq-kj-6Xs2Rw9w3Y3LlK_ffFip1QzemHdesiug0Wko,817
-Code/test.py,sha256=rbvsRbDGzSifUU1dSUA7CHkU58Qt8z8U42A9hGdIMOQ,212
-Code/test1.py,sha256=rGA9d1cgArhAfxRAAdtiz-vPGQ0Ns3HEdN8sCwydHYM,3250
-TDCRPy-0.0.6.dist-info/LICENCE.md,sha256=fuYzrZRiOAjJBzA1tsGQwojCgGROArb2Ec48GDTjlWM,1086
-TDCRPy-0.0.6.dist-info/METADATA,sha256=2i9g8aqPHJVoECr8P8yaquCvQXmXzjyV8cfNhyTnIZY,950
-TDCRPy-0.0.6.dist-info/WHEEL,sha256=2wepM1nk4DS4eFpYrW1TTqPcoGNfHhhO_i5m4cOimbo,92
-TDCRPy-0.0.6.dist-info/top_level.txt,sha256=5dSMF0KOskVmT0QGz9KnqcD9Ccc2Dkq0Sa8FOZk-7Ck,5
-TDCRPy-0.0.6.dist-info/RECORD,,