TDCRPy 1.9.6__py3-none-any.whl → 1.9.8__py3-none-any.whl

{TDCRPy-1.9.6.dist-info → TDCRPy-1.9.8.dist-info}/METADATA

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

{TDCRPy-1.9.6.dist-info → TDCRPy-1.9.8.dist-info}/RECORD

@@ -1,6 +1,6 @@
-tdcrpy/TDCRPy.py,sha256=qYl8BEI2lLm2oPGrnrxiFya7Y_aQU6dtyutUYHay4Yc,53858
+tdcrpy/TDCRPy.py,sha256=kl0UM0HEY7D95M1c806IkYkXl7du1BMtWUA7jTiOGVs,53859
 tdcrpy/TDCRPy1.py,sha256=QTBZh5B5JWnGB0BQfD-cFmwA9W080OD4sG-aj50-ejo,38106
-tdcrpy/TDCR_model_lib.py,sha256=YfOTuOqywl9gMqP7FlJo6I28obn4Qyv3zUHTRW7HZYs,104753
+tdcrpy/TDCR_model_lib.py,sha256=Lfl2xIXbleN-W1t3ad6ki6wKgttXNpww7KITcb3SgPg,104751
 tdcrpy/TDCRoptimize.py,sha256=c2XIGveeLdVYYek4Rg6dygMvVA2xIrIkMb3L-_jUucM,6496
 tdcrpy/__init__.py,sha256=vQslGLsoZPIceaitnSHOqN6lUdjEyJ3YhfJ6tYdXt-s,127
 tdcrpy/config.toml,sha256=rHuuM520HZquGrEd8LAVUdvnzZGtk90pigudy-MXnm8,1389
@@ -1068,7 +1068,7 @@ tdcrpy/Quenching/TandataUG.txt,sha256=AkUGz2dAZ6g3bQIIacR5f92aa8gmiklgRbN4Dyz-12
 tdcrpy/Quenching/alpha_toulene.txt,sha256=rJB9xCLM3GLg91EKfhu8sGwro281ySq6R0G5mXMt3f8,1553
 tdcrpy/Quenching/inputVecteurAlpha.txt,sha256=ee4emplX51vhomC7cGNP_eE5eqS9By5PaSSaa4hh1eg,178652
 tdcrpy/Quenching/inputVecteurElectron.txt,sha256=5NBRBaiB7op1SJ6v952w9gFKqh0x83pn8zYcoUpy33Y,181498
-tdcrpy/decayData/All-nuclides_BetaShape.zip,sha256=JfbU66QYnYIki7nmwbXBnf7MDOa7JDBqFv4hEsEQQyM,7130930
+tdcrpy/decayData/All-nuclides_BetaShape.zip,sha256=-qTbYuzgTJv9AyJ2x5zX-9ndWBY3SuNyS2LBDP5iSwg,7633072
 tdcrpy/decayData/All-nuclides_Ensdf.zip,sha256=GL-X_IG6x4ykGrWYaMiZrjfMXOvYPM07G-Cp8WCPOu8,437902
 tdcrpy/decayData/All-nuclides_PenNuc.zip,sha256=DmhS2CYegWhL8qggc_MO6FgyiXE5qX73dOoH-SSQXoo,388954
 tdcrpy/decayData/atom-ENDF-VII0.zip,sha256=oP4Z7qWweWF4j12QAIfvLFLlqeY1RomXwXD2iFgQCqo,1624444
@@ -1080,8 +1080,8 @@ tdcrpy/docs/_build/html/source/modules.html,sha256=Jf-qxVBId0UgpwyvYuyjtMNG-ezPO
 tdcrpy/docs/_build/html/source/tdcrpy.html,sha256=-38lHMNFB22p1tWJEeN3yDqfDiCYE304vxDamO1-iRc,3779
 tdcrpy/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 tdcrpy/test/test_tdcrpy.py,sha256=MTMiE0N7RiyE8Rom2Yu8kCDR03pJga9D6_o7wr7leVA,4675
-TDCRPy-1.9.6.dist-info/LICENCE.md,sha256=ZTpWyGU3qv_iwEpgvCijoCuCYpOPpyzJCgOk46WpUKU,1066
-TDCRPy-1.9.6.dist-info/METADATA,sha256=tA4rnXeUBedoSDn3z2B1bPcY7GQmE3CSq1kX27q4Ego,15995
-TDCRPy-1.9.6.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
-TDCRPy-1.9.6.dist-info/top_level.txt,sha256=f4vzFFcKSEnonAACs0ZXuRczmroLLqtPTqXFymU_VU0,14
-TDCRPy-1.9.6.dist-info/RECORD,,
+TDCRPy-1.9.8.dist-info/LICENCE.md,sha256=ZTpWyGU3qv_iwEpgvCijoCuCYpOPpyzJCgOk46WpUKU,1066
+TDCRPy-1.9.8.dist-info/METADATA,sha256=qqblBiwN8n0Dj8GWVq4HO5_ZUm-nzcI9G7BtX8o_OIo,15995
+TDCRPy-1.9.8.dist-info/WHEEL,sha256=GJ7t_kWBFywbagK5eo9IoUwLW6oyOeTKmQ-9iHFVNxQ,92
+TDCRPy-1.9.8.dist-info/top_level.txt,sha256=f4vzFFcKSEnonAACs0ZXuRczmroLLqtPTqXFymU_VU0,14
+TDCRPy-1.9.8.dist-info/RECORD,,

tdcrpy/TDCRPy.py

@@ -179,9 +179,9 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
         if mode == "eff":
             return out[0], 0, out[1], 0, out[2], 0
     elif (not Smodel) and (not Rad in radListPureBeta):
-        print("cannot be processed by the analytical model.")
+        # print("cannot be processed by the analytical model.")
         # print(f"Analytical model used for {Rad}")
-        out=tl.modelAnalytical(L,TD,TAB,TBC,TAC,Rad,kB,V,mode,mode2,10000)
+        out=tl.modelAnalytical(L,TD,TAB,TBC,TAC,Rad,kB,V,mode,mode2,1000)
         if mode == "res":
             return out
         if mode == "eff":

tdcrpy/TDCR_model_lib.py

@@ -824,7 +824,7 @@ def readBetaShape(rad,mode,level,z=z_betashape):
         else:
             p.append(p0 * (e[k]-e[k-1]))
 
-    e=(np.asarray(e[:-1])+np.asarray(e[1:]))/2
+    e=(np.asarray(e[:-1])+np.asarray(e[1:]))/2 # deal with the continuity hypothesis
     p.pop(-1)
     p /= sum(np.asarray(p)) # normalization
     p = list(p); e = list(e)
@@ -2524,9 +2524,8 @@ def modelAnalytical(L,TD,TAB,TBC,TAC,rad,kB,V,mode,mode2,ne):
     e, p = readBetaSpectra(rad)
     em=np.empty(len(e))
     for i, ei in enumerate(e):
-        # ed = energie_dep_beta2(ei,V)
-        em[i] = E_quench_e(ei*1e3,ei*1e3,kB*1e3,ne)*1e-3
-        # em[i] = Em_e(ei*1e3,ei*1e3,kB*1e3,ne)*1e-3
+        #em[i] = E_quench_e(ei*1e3,ei*1e3,kB*1e3,ne)*1e-3
+        em[i] = Em_e(ei*1e3,ei*1e3,kB*1e3,ne)*1e-3
         
         
     if mode2=="sym":