TDCRPy 1.12.6__tar.gz → 1.12.9__tar.gz

{tdcrpy-1.12.6/TDCRPy.egg-info → tdcrpy-1.12.9}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: TDCRPy
-Version: 1.12.6
+Version: 1.12.9
 Summary: TDCR model
 Home-page: https://pypi.org/project/TDCRPy/
 Author: RomainCoulon (Romain Coulon)

{tdcrpy-1.12.6 → tdcrpy-1.12.9/TDCRPy.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: TDCRPy
-Version: 1.12.6
+Version: 1.12.9
 Summary: TDCR model
 Home-page: https://pypi.org/project/TDCRPy/
 Author: RomainCoulon (Romain Coulon)

{tdcrpy-1.12.6 → tdcrpy-1.12.9}/setup.py

@@ -2,7 +2,7 @@ from setuptools import setup, find_packages
 import codecs
 import os
 
-VERSION = "1.12.6"
+VERSION = "1.12.9"
 
 DESCRIPTION = "TDCR model"
 

{tdcrpy-1.12.6 → tdcrpy-1.12.9}/tdcrpy/TDCRPy.py

@@ -166,6 +166,8 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
         file_conf = data_path       
     config.read(file_conf)
     tau=config["Inputs"].getfloat("tau")
+    extDT=config["Inputs"].getfloat("extDT")
+    measTime=config["Inputs"].getfloat("measTime")
     micCorr=config["Inputs"].getboolean("micCorr")
     radListPureBeta = ["H-3", "C-14", "S-35", "Ca-45", "Ni-63", "Sr-89", "Sr-90", "Tc-99", "Pm-147", "Pu-241"]
     if (Rad in radListPureBeta) and (not Smodel):
@@ -694,7 +696,7 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
             ====================
             '''            
             if mode2=="sym":
-                if evenement !=1:
+                if evenement !=1 and t1 > extDT*1e-6 and t1 < measTime*60:
                     p_nosingle = np.exp(-L*np.sum(np.asarray(e_quenching))/3) # probability to have 0 electrons in a PMT
                     p_single = 1-p_nosingle                                    # probability to have at least 1 electrons in a PMT
                     p_nosingle2 = np.exp(-L*np.sum(np.asarray(e_quenching2))/3) # probability to have 0 electrons in a PMT
@@ -727,7 +729,7 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
                     if Display: print("\t\t Efficiency of triple events = ", round(efficiency_T[-1],5))                    
                                     
             elif mode2=="asym":
-                if evenement !=1:
+                if evenement !=1 and t1 > extDT*1e-6 and t1 < measTime*60:
                     pA_nosingle = np.exp(-L[0]*np.sum(np.asarray(e_quenching))/3) # probability to have 0 electrons in a PMT
                     pA_single = 1-pA_nosingle                                    # probability to have at least 1 electrons in a PMT
                     pB_nosingle = np.exp(-L[1]*np.sum(np.asarray(e_quenching))/3) # probability to have 0 electrons in a PMT

{tdcrpy-1.12.6 → tdcrpy-1.12.9}/tdcrpy/TDCR_model_lib.py

@@ -41,6 +41,8 @@ def readParameters(disp=False):
     nE_electron = config["Inputs"].getint("nE_electron")
     nE_alpha = config["Inputs"].getint("nE_alpha")
     tau = config["Inputs"].getint("tau")
+    extDT = config["Inputs"].getint("extDT")
+    measTime = config["Inputs"].getint("measTime")
     RHO = config["Inputs"].getfloat("density")
     Z = config["Inputs"].getfloat("Z")
     A = config["Inputs"].getfloat("A")
@@ -49,6 +51,7 @@ def readParameters(disp=False):
     Einterp_e = config["Inputs"].getfloat("Einterp_e")
     diam_micelle = config["Inputs"].getfloat("diam_micelle")
     fAq = config["Inputs"].getfloat("fAq")
+    micCorr = config["Inputs"].getboolean("micCorr")
     
     if disp:
         print(f"number of integration bins for electrons = {nE_electron}")
@@ -59,13 +62,16 @@ def readParameters(disp=False):
         print(f"depth of spline interp. = {depthSpline}")
         print(f"energy above which interp. in implemented (for alpha) = {Einterp_a} keV")
         print(f"energy above which interp. in implemented (for electron) = {Einterp_e} keV")
+        print(f"activation of the micelle correction = {micCorr}")
         print(f"diameter of micelle = {diam_micelle} nm")
         print(f"acqueous fraction = {fAq}")
         print(f"coincidence resolving time = {tau} ns")
+        print(f"extended dead time = {extDT} µs")
+        print(f"measurement time = {measTime} min")
     
-    return nE_electron, nE_alpha, RHO, Z, A, depthSpline, Einterp_a, Einterp_e, diam_micelle, fAq, tau
+    return nE_electron, nE_alpha, RHO, Z, A, depthSpline, Einterp_a, Einterp_e, diam_micelle, fAq, tau, extDT, measTime, micCorr
 
-nE_electron, nE_alpha, RHO, Z, A, depthSpline, Einterp_a, Einterp_e, diam_micelle, fAq, tau = readParameters()
+nE_electron, nE_alpha, RHO, Z, A, depthSpline, Einterp_a, Einterp_e, diam_micelle, fAq, tau, extDT, measTime, micCorr = readParameters()
 
 def readConfigAsstr():
     path2config = str(config.read(file_conf)[0])
@@ -144,6 +150,24 @@ def modifyTau(x):
     data1 = data0.replace(f"tau = {x0}",f"tau = {x}")
     writeConfifAsstr(data1)
 
+def modifyDeadTime(x):
+    data0 = readConfigAsstr()
+    x0 = readParameters()[11]
+    data1 = data0.replace(f"extDT = {x0}",f"extDT=  {x}")
+    writeConfifAsstr(data1)
+
+def modifyMeasTime(x):
+    data0 = readConfigAsstr()
+    x0 = readParameters()[12]
+    data1 = data0.replace(f"measTime = {x0}",f"measTime = {x}")
+    writeConfifAsstr(data1)
+
+def modifyMicCorr(x):
+    data0 = readConfigAsstr()
+    x0 = readParameters()[13]
+    data1 = data0.replace(f"micCorr = {x0}",f"micCorr= {x}")
+    writeConfifAsstr(data1)
+
 # import PenNuc data
 with importlib.resources.as_file(files('tdcrpy').joinpath('decayData')) as data_path:
 #with importlib.resources.path('tdcrpy', 'decayData') as data_path:

{tdcrpy-1.12.6 → tdcrpy-1.12.9}/tdcrpy/config.toml

@@ -34,6 +34,10 @@ fAq = 0.1
 ## PROPERTIES OF THE COUNTER
 # Coincidence resolving time (ns)
 tau = 50
+# Extended dead time (µs)
+extDT = 10
+# Measurement time (minutes)
+measTime = 20
 
 ## INTERPOLATION TO ESTIMATE THE QUENCHED ENERGY
 # depth (number of indexes on each side of the energy point) on which the spline interpolation is done.