TDCRPy 1.6.0__py3-none-any.whl → 1.6.2__py3-none-any.whl

{TDCRPy-1.6.0.dist-info → TDCRPy-1.6.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: TDCRPy
-Version: 1.6.0
+Version: 1.6.2
 Summary: TDCR model
 Home-page: https://github.com/RomainCoulon/TDCRPy
 Author: RomainCoulon (Romain Coulon)

{TDCRPy-1.6.0.dist-info → TDCRPy-1.6.2.dist-info}/RECORD

@@ -1,6 +1,6 @@
-tdcrpy/TDCRPy.py,sha256=5imzxtIfyzkRWbfOTOvrbbyMWnjnfPIOFoWoBYN6j00,57307
+tdcrpy/TDCRPy.py,sha256=02OJKQnFA3E1x4txLYNiomiDCvH_pmogosoy2vOyMC0,57647
 tdcrpy/TDCRPy1.py,sha256=ZS4riwGPyITxkg0RHxTO59yzue7uLSFa9S7UrM8NGOI,38777
-tdcrpy/TDCR_model_lib.py,sha256=dJ6fj-QGZazLsyKkuAV9rsNbEkeerjECtD3EOlPPJH0,75156
+tdcrpy/TDCR_model_lib.py,sha256=wkVoZ8Afb-DEVDfYzpxDFZCT2xfzn-PBMw81TG6mYBs,75627
 tdcrpy/TDCRoptimize.py,sha256=Mo5vEULS1_FCEzIVyv0C6Yte6qo072HLRWETNyIs7SU,4552
 tdcrpy/__init__.py,sha256=pxZ5gMTFSEzeGQLzeRjy-Imu7bZIfhg9eqloOq5FTM4,130
 tdcrpy/config.toml,sha256=mWDY-Fm3SxO9zsYzkjt9iQGWwzOH7z4CuCN1I8ASOYA,1542
@@ -68,8 +68,8 @@ tdcrpy/docs/_build/html/source/modules.html,sha256=gJrhuEaMQHkYGewGHPWGmgVsPgVV6
 tdcrpy/docs/_build/html/source/tdcrpy.html,sha256=VZV_LFdrCCIvYBXzsZE9mZx8iGa7ZjhEsXtE8BqqU9g,90373
 tdcrpy/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 tdcrpy/test/test_tdcrpy.py,sha256=dWZpfvarkr8lr_K_zUakt2AOcF0FfCOjBOQPl4vnTWc,4116
-TDCRPy-1.6.0.dist-info/LICENCE.md,sha256=fuYzrZRiOAjJBzA1tsGQwojCgGROArb2Ec48GDTjlWM,1086
-TDCRPy-1.6.0.dist-info/METADATA,sha256=rQ_NR-_ell_qDpnqRQo-mau6Ukl4l8oVRjBqAD6tLMM,15929
-TDCRPy-1.6.0.dist-info/WHEEL,sha256=AtBG6SXL3KF_v0NxLf0ehyVOh0cold-JbJYXNGorC6Q,92
-TDCRPy-1.6.0.dist-info/top_level.txt,sha256=VgPJa9YJ7fi8wrr9qDJPLhu3tK1BcFRIP6Ib4r3BEfs,7
-TDCRPy-1.6.0.dist-info/RECORD,,
+TDCRPy-1.6.2.dist-info/LICENCE.md,sha256=fuYzrZRiOAjJBzA1tsGQwojCgGROArb2Ec48GDTjlWM,1086
+TDCRPy-1.6.2.dist-info/METADATA,sha256=9RbYa5ZGqZvWXlF4JiXYq2BwLNXbWRG6B5eT8zyvsMw,15929
+TDCRPy-1.6.2.dist-info/WHEEL,sha256=AtBG6SXL3KF_v0NxLf0ehyVOh0cold-JbJYXNGorC6Q,92
+TDCRPy-1.6.2.dist-info/top_level.txt,sha256=VgPJa9YJ7fi8wrr9qDJPLhu3tK1BcFRIP6Ib4r3BEfs,7
+TDCRPy-1.6.2.dist-info/RECORD,,

tdcrpy/TDCRPy.py

@@ -45,9 +45,17 @@ def relaxAtom(daughter_relax,particle_vec,energy_vec,rad,Display=False,uncData=F
                 particle_vec[i_part] = "Atom_L"
                 tf1,ef1 = tl.relaxation_atom(daughter_relax,rad,'Atom_L',uncData=uncData)
                 particle_vec.append(tf)
-                particle_vec.append(tf1)
                 energy_vec.append(ef)
+                particle_vec.append(tf1)
                 energy_vec.append(ef1)
+                if tf1 == 'Auger L':
+                    particle_vec.append("Atom_M")
+                    particle_vec.append("Atom_M")
+                    energy_vec.append(0)
+                    energy_vec.append(0)
+                else:
+                    particle_vec.append("Atom_M")
+                    energy_vec.append(0)
                 relaxation = True
             elif tf == "Auger KLX":
                 particle_vec[i_part] = "Atom_L"

tdcrpy/TDCR_model_lib.py

@@ -1410,6 +1410,8 @@ def read_ENSDF(rad, *, z=z_ensdf):
     prob_augerK = []      # enregistrer les résultats (prob d'électron auger K) d'une fille
     Type_augerK = []
     type_augerK = []
+    energie_augerK = []
+    Energie_augerK = []
     
     for i in range(len(posi)-1):
         start = posi[i]+1
@@ -1422,6 +1424,7 @@ def read_ENSDF(rad, *, z=z_ensdf):
         prob_str_b = []       # enregistrer les résultats (proba en str) d'un bloc
         prob_augerK_b = []    # enregistrer les résultats (prob d'électron auger K) d'un bloc
         type_augerK_b = []
+        energie_augerK_b = []
         
         if start==end:        # sauter les lignes blaches et continues
             continue
@@ -1442,7 +1445,8 @@ def read_ENSDF(rad, *, z=z_ensdf):
                     prob_str_b.append(p1[3])
                 continue 
             elif '|]' in p1:                   # traiter un bloc qui comprend |]
-                if 'AUGER' in p1:              # block of electron Auger 
+                if 'AUGER' in p1:              # block of electron Auger
+                    energie_augerK_b.append(float(p1[2]))
                     type_b.append('Auger K')   # for electron Auger, only the block of electron Auger K has |] 
                     if len(p1)>7:               # repérer la ligne qui comprend la proba totale et l'incertitude
                         prob_str_b.append(p1[5])
@@ -1485,6 +1489,7 @@ def read_ENSDF(rad, *, z=z_ensdf):
             if len(prob_augerK_b)!=0:              # si l'électron Auger K, enregistrer les proba dans la liste pour un noyau fils
                 prob_augerK = prob_augerK_b
                 type_augerK = type_augerK_b
+                energie_augerK = energie_augerK_b
         elif len(e)==len(prob_b) and len(e)>=1:    # enregistrer les valeurs au cas où sans |] et valeurs complètes
             for i in range(len(e)):
                 energy.append(e[i])
@@ -1495,6 +1500,7 @@ def read_ENSDF(rad, *, z=z_ensdf):
             if len(prob_augerK_b)!=0:               # si l'électron Auger, enregistrer les proba dans la liste pour un noyau fils 
                 prob_augerK = prob_augerK_b         # pour certains cas spécifiques qui n'ont qu'un seul type d'électron Auger
                 type_augerK = type_augerK_b
+                energie_augerK = energie_augerK_b
         if end in index_end or end+1 in index_end:  # enregistrer les résultats à la fin d'une fille   
             Energy.append(energy)
             Prob.append(prob)
@@ -1503,6 +1509,7 @@ def read_ENSDF(rad, *, z=z_ensdf):
             Prob_str.append(prob_str)
             Prob_augerK.append(prob_augerK)
             Type_augerK.append(type_augerK)
+            Energie_augerK.append(energie_augerK)
             energy = []
             prob = []
             type_ = []   
@@ -1510,9 +1517,9 @@ def read_ENSDF(rad, *, z=z_ensdf):
             prob_str = []
             prob_augerK = []
             type_augerK = []
+            energie_augerK = []
             
-    return  daug_name,Energy,Prob,Type,Incertitude,Prob_str,Prob_augerK,Type_augerK  
-
+    return  daug_name,Energy,Prob,Type,Incertitude,Prob_str,Prob_augerK,Type_augerK,Energie_augerK  
 
 
 def incer(prob,incer):
@@ -1581,7 +1588,7 @@ def relaxation_atom(daugther,rad,lacune='defaut',uncData=False):
     Energy : corresponding energy in keV.
 
     """
-    daug_name,Energy,Prob,Type,Incertitude,prob_str,Prob_K,Type_K = read_ENSDF(rad)  # tirer les vecteurs de rad d'Ensdf 
+    daug_name,Energy,Prob,Type,Incertitude,prob_str,Prob_K,Type_K,Energie_augerK = read_ENSDF(rad)  # tirer les vecteurs de rad d'Ensdf 
     incertitude = incer(prob_str,Incertitude)
 
     index_daug = daug_name.index(daugther)        # repérer l'indice de fille correspondante
@@ -1592,6 +1599,7 @@ def relaxation_atom(daugther,rad,lacune='defaut',uncData=False):
     prob_augerK = Prob_K[index_daug]    # tirer le vecteur de prob 'électron auger K
     type_augerK = Type_K[index_daug]
     u_probability = np.array(incertitude[index_daug])
+    energie_augerK = np.array(Energie_augerK[index_daug])
     
     if len(probability) > 0:                      # le cas où le vecteur de proba/energie/type n'est pas vide
         if 'L' in lacune:                         # traiter le transition de couche L
@@ -1653,6 +1661,7 @@ def relaxation_atom(daugther,rad,lacune='defaut',uncData=False):
                     index_K = sampling(prob_AugerK)
                     type_K = type_augerK[index_K]
                     type_fin = 'Auger K' + type_K[1:]
+                    energie_fin = energie_augerK[index_K]
                 
             else:
                 # print("pas de transition de rayon X ni d'électron Auger pour cette lacune: ",lacune)