TDCRPy 0.0.1__py3-none-any.whl

Code/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)

Code/__init__.py

@@ -0,0 +1 @@
+from TDCRPy import TDCRPy

Code/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')

Code/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.1.dist-info/METADATA

@@ -0,0 +1,21 @@
+Metadata-Version: 2.1
+Name: TDCRPy
+Version: 0.0.1
+Summary: TDCR model
+Home-page: UNKNOWN
+Author: RomainCoulon (Romain Coulon)
+Author-email: <romain.coulon@bipm.org>
+License: UNKNOWN
+Keywords: python,TDCR
+Platform: UNKNOWN
+Classifier: Development Status :: 1 - Planning
+Classifier: Intended Audience :: Developers
+Classifier: Programming Language :: Python :: 3
+Classifier: Operating System :: Unix
+Classifier: Operating System :: MacOS :: MacOS X
+Classifier: Operating System :: Microsoft :: Windows
+License-File: LICENCE
+Requires-Dist: numpy
+
+UNKNOWN
+

TDCRPy-0.0.1.dist-info/RECORD

@@ -0,0 +1,13 @@
+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=SUEsH4EdNiXOzPbsBgLHk3OmEti8M7guimMrV0HpP4o,25
+Code/decay.py,sha256=Lnq-kj-6Xs2Rw9w3Y3LlK_ffFip1QzemHdesiug0Wko,817
+Code/test1.py,sha256=rGA9d1cgArhAfxRAAdtiz-vPGQ0Ns3HEdN8sCwydHYM,3250
+TDCRPy-0.0.1.dist-info/LICENCE,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+TDCRPy-0.0.1.dist-info/METADATA,sha256=WPlPzU3GXr_yyOrvbFyihNPp0ToVOO8wnDfgtwhfswI,574
+TDCRPy-0.0.1.dist-info/WHEEL,sha256=2wepM1nk4DS4eFpYrW1TTqPcoGNfHhhO_i5m4cOimbo,92
+TDCRPy-0.0.1.dist-info/top_level.txt,sha256=5dSMF0KOskVmT0QGz9KnqcD9Ccc2Dkq0Sa8FOZk-7Ck,5
+TDCRPy-0.0.1.dist-info/RECORD,,

TDCRPy-0.0.1.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: bdist_wheel (0.38.4)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

TDCRPy-0.0.1.dist-info/top_level.txt

@@ -0,0 +1 @@
+Code