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

tdcrpy/TDCRPy.py

@@ -19,6 +19,7 @@ import numpy as np
 from tqdm import tqdm
 
 def relaxAtom(daughter_relax,particle_vec,energy_vec,rad,Display=False,uncData=False):
+         
     for i_part in range(len(particle_vec)):
         relaxation = False
         if "Atom_K" in particle_vec[i_part] or "Atom_L" in particle_vec[i_part] or "Atom_M" in particle_vec[i_part]:
@@ -40,22 +41,37 @@ def relaxAtom(daughter_relax,particle_vec,energy_vec,rad,Display=False,uncData=F
                 particle_vec.append(tf)
                 energy_vec.append(ef)
                 relaxation = False
-            elif tf == "Auger K":
+            elif tf == "Auger KLL":
                 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("Atom_L")
-                energy_vec.append(0)
+                energy_vec.append(ef1)
                 relaxation = True
-            elif tf == "Auger L":
+            elif tf == "Auger KLX":
+                particle_vec[i_part] = "Atom_L"
+                particle_vec.append("Atom_M")
+                particle_vec.append(tf)
+                energy_vec.append(0)
+                energy_vec.append(ef)
+                relaxation = True    
+            elif tf == "Auger KXY":
                 particle_vec[i_part] = "Atom_M"
+                particle_vec.append("Atom_M")
                 particle_vec.append(tf)
+                energy_vec.append(0)
                 energy_vec.append(ef)
+                relaxation = False    
+            elif tf == "Auger L":
+                particle_vec[i_part] = "Atom_M"
                 particle_vec.append("Atom_M")
+                particle_vec.append(tf)
                 energy_vec.append(0)
+                energy_vec.append(ef)
                 relaxation = False
             else:
-                if Display: print("\t\t untermined x or Auger")
+                if Display: print(f"\t\t untermined x or Auger = {tf}")
                 relaxation = False
     return particle_vec, energy_vec
 
@@ -471,7 +487,67 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
 
             daughter_relax = DaughterVec[index_rad][iDaughter]
             particle_vec, energy_vec = relaxAtom(daughter_relax,particle_vec,energy_vec,Rad[index_rad],Display=Display,uncData=uncData)
+            
+            # for i_part in range(len(particle_vec)):
+            #     relaxation = False
+            #     if "Atom_K" in particle_vec[i_part] or "Atom_L" in particle_vec[i_part] or "Atom_M" in particle_vec[i_part]:
+            #         relaxation = True
+            #     while relaxation:
+            #         tf,ef = tl.relaxation_atom(daughter_relax,Rad[index_rad],particle_vec[i_part],uncData=uncData)
+                    
+            #         if tf == "XKA":
+            #             particle_vec[i_part] = "Atom_L"
+            #             particle_vec.append(tf)
+            #             energy_vec.append(ef)
+            #             relaxation = True
+            #         elif tf == "XKB":
+            #             particle_vec[i_part] = "Atom_M"
+            #             particle_vec.append(tf)
+            #             energy_vec.append(ef)
+            #             relaxation = False
+            #         elif tf == "XL":
+            #             particle_vec[i_part] = "Atom_M"
+            #             particle_vec.append(tf)
+            #             energy_vec.append(ef)
+            #             relaxation = False
+            #         elif tf == "Auger KLL":
+            #             particle_vec[i_part] = "Atom_L"
+            #             tf1,ef1 = tl.relaxation_atom(daughter_relax,Rad[index_rad],'Atom_L',uncData=uncData)
+            #             particle_vec.append(tf)
+            #             particle_vec.append(tf1)
+            #             energy_vec.append(ef)
+            #             energy_vec.append(ef1)
+            #             relaxation = True
+            #         elif tf == "Auger KLX":
+            #             particle_vec[i_part] = "Atom_L"
+            #             particle_vec.append("Atom_M")
+            #             particle_vec.append(tf)
+            #             energy_vec.append(0)
+            #             energy_vec.append(ef)
+            #             relaxation = True    
+            #         elif tf == "Auger KXY":
+            #             particle_vec[i_part] = "Atom_M"
+            #             particle_vec.append("Atom_M")
+            #             particle_vec.append(tf)
+            #             energy_vec.append(0)
+            #             energy_vec.append(ef)
+            #             relaxation = False    
+            #         elif tf == "Auger L":
+            #             particle_vec[i_part] = "Atom_M"
+            #             particle_vec.append("Atom_M")
+            #             particle_vec.append(tf)
+            #             energy_vec.append(0)
+            #             energy_vec.append(ef)
+            #             relaxation = False
+            #         else:
+            #             if Display: print("\t\t untermined x or Auger")
+            #             relaxation = False
+            #         e_sum += ef
+                    
+
+
             ## evenement normal
+
             if Display:
                 print("\n\t ATOMIC RECOMBINATION--Prompt")
                 for i, p in enumerate(particle_vec):
@@ -485,7 +561,69 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
             
             ## evenement retardee             
             if evenement != 1:
+                
+                if Display:print("\n\t ATOMIC RECOMBINATION--Delay\n\t Summary of the atomic relaxation")
+                
                 particle_vec2, energy_vec2 = relaxAtom(daughter_relax,particle_vec2,energy_vec2,Rad[index_rad],Display=Display,uncData=uncData)
+                
+                
+                # for i_part in range(len(particle_vec2)):
+                #     relaxation = False
+                #     if "Atom_K" in particle_vec2[i_part] or "Atom_L" in particle_vec2[i_part] or "Atom_M" in particle_vec2[i_part]:
+                #         relaxation = True
+                #     while relaxation:
+                #         tf,ef = tl.relaxation_atom(daughter_relax,Rad[index_rad],particle_vec2[i_part],uncData=uncData)
+                        
+                #         if tf == "XKA":
+                #             particle_vec2[i_part] = "Atom_L"
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(ef)
+                #             relaxation = True
+                #         elif tf == "XKB":
+                #             particle_vec2[i_part] = "Atom_M"
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(ef)
+                #             relaxation = False
+                #         elif tf == "XL":
+                #             particle_vec2[i_part] = "Atom_M"
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(ef)
+                #             relaxation = False
+                #         elif tf == "Auger KLL":
+                #             particle_vec2[i_part] = "Atom_L"
+                #             tf1,ef1 = tl.relaxation_atom(daughter_relax,Rad[index_rad],'Atom_L',uncData=uncData)
+                #             particle_vec2.append(tf)
+                #             particle_vec2.append(tf1)
+                #             energy_vec2.append(ef)
+                #             energy_vec2.append(ef1)
+                #             relaxation = True
+                #         elif tf == "Auger KLX":
+                #             particle_vec2[i_part] = "Atom_L"
+                #             particle_vec2.append("Atom_M")
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(0)
+                #             energy_vec2.append(ef)
+                #             relaxation = True 
+                #         elif tf == "Auger KXY":
+                #             particle_vec2[i_part] = "Atom_M"
+                #             particle_vec2.append("Atom_M")
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(0)
+                #             energy_vec2.append(ef)
+                #             relaxation = False         
+                #         elif tf == "Auger L":
+                #             particle_vec2[i_part] = "Atom_M"
+                #             particle_vec2.append(tf)
+                #             energy_vec2.append(ef)
+                #             relaxation = False
+                #         else:
+                #             if Display: print("\t\t x ray or Auger electron from X shell")
+                #             relaxation = False
+                #         e_sum2 += ef
+                
+
+                
+
                 if Display:
                     print("\n\t ATOMIC RECOMBINATION--Delay")
                     for i, p in enumerate(particle_vec2):
@@ -762,3 +900,17 @@ def TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=Fals
                 return mean_efficiency_S, std_efficiency_S, mean_efficiency_D, std_efficiency_D, mean_efficiency_T, std_efficiency_T
         if mode =="dis":
             return efficiency_S, efficiency_D, efficiency_T    
+# L = 1
+# TD = 0.977667386529166
+# TAB = 0.992232838598821
+# TBC = 0.992343419459002
+# TAC = 0.99275350064608
+# Rad="Fe-55"
+# pmf_1="1"
+# N = 10
+# kB =1.0e-5
+# V = 10
+# mode = "eff"
+# mode2 = "sym"
+
+# out = TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2, Display=True, barp=False,uncData=False)

tdcrpy/TDCR_model_lib.py

@@ -81,6 +81,8 @@ with importlib.resources.as_file(files('tdcrpy').joinpath('MCNP-MATRIX')) as dat
     fe4 = data_path / 'matrice/fichier/matrice_16ml-beta-_1_200k.txt' # electron-16ml-1-200keV-niveau 0
     fe5 = data_path / 'matrice/fichier/matrice_16ml-beta-_200_2000k.txt' # electron-16ml-200-2000keV-niveau 1
     fe6 = data_path / 'matrice/fichier/matrice_16ml-beta-_2000_10000k.txt' # electron-16ml-2000-10000keV-niveau 2
+    fe7 = data_path / 'matrice/fichier/matrice_13ml-beta-_1_200k.txt' # electron-13ml-1-200keV-niveau 0
+    fe8 = data_path / 'matrice/fichier/matrice_13ml-beta-_200_2000k.txt' # electron-13ml-200-2000keV-niveau 1
     fe = data_path / 'matrice/fichier/E_depose.txt' # electron-10ml-énergie-niveau 'e'   
 
 # import beta spectra calculated for the analytical model (BetaShape + MCNP6 calculation) 
@@ -996,9 +998,15 @@ Matrice10_p_3 = read_matrice(fp3,2)
 Matrice16_p_1 = read_matrice(fp4,0)
 Matrice16_p_2 = read_matrice(fp5,1)
 Matrice16_p_3 = read_matrice(fp6,2)
+
+#Matrice13_p_1 = read_matrice(fp7,0)
+
 Matrice13_p_1 = read_matrice(fp7,0)
 Matrice13_p_2 = read_matrice(fp8,1)
 Matrice13_p_3 = read_matrice(fp9,2)
+Matrice13_e_1 = read_matrice(fe7,0)
+Matrice13_e_2 = read_matrice(fe8,1)
+
 Matrice_e = read_matrice(fe,'e')
 
 Matrice10_e_1 = read_matrice(fe1,0)
@@ -1009,7 +1017,7 @@ Matrice16_e_2 = read_matrice(fe5,1)
 Matrice16_e_3 = read_matrice(fe6,2)
 #Matrice_e = read_matrice(fe,'e')
 
-def energie_dep_gamma(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p_2,matrice10_3=Matrice10_p_3,matrice16_1=Matrice16_p_1,matrice16_2=Matrice16_p_2,matrice16_3=Matrice16_p_3,matrice13_1=Matrice13_p_1,ed=Matrice_e):
+def energie_dep_gamma(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p_2,matrice10_3=Matrice10_p_3,matrice16_1=Matrice16_p_1,matrice16_2=Matrice16_p_2,matrice16_3=Matrice16_p_3,matrice13_1=Matrice13_p_1,matrice13_2=Matrice13_p_2,matrice13_3=Matrice13_p_3,ed=Matrice_e):
     """ This function samples the energy deposited by a x or gamma rays in the scintillator using response calculated by the Monte-Carlo code MCNP6. 
     
     Parameters
@@ -1053,6 +1061,8 @@ def energie_dep_gamma(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p
         index = int((e_inci-200)/2)
         if v == 10: 
             matrice = matrice10_2
+        elif v ==13:
+            matrice = matrice13_2    
         elif v == 16:
             matrice = matrice16_2
         e = ed[:,1]
@@ -1061,6 +1071,8 @@ def energie_dep_gamma(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p
         index = (int(e_inci)-2000)//10
         if v == 10: 
             matrice = matrice10_3
+        elif v ==13:
+            matrice = matrice13_3    
         elif v == 16:
             matrice = matrice16_3
         e = ed[:,2]
@@ -1072,7 +1084,7 @@ def energie_dep_gamma(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p
     if result  > e_inci: result = e_inci
     return result
 
-def energie_dep_gamma2(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p_2,matrice10_3=Matrice10_p_3,matrice16_1=Matrice16_p_1,matrice16_2=Matrice16_p_2,matrice16_3=Matrice16_p_3,matrice13_1=Matrice13_p_1, ed=Matrice_e):
+def energie_dep_gamma2(e_inci,v,matrice10_1=Matrice10_p_1,matrice10_2=Matrice10_p_2,matrice10_3=Matrice10_p_3,matrice16_1=Matrice16_p_1,matrice16_2=Matrice16_p_2,matrice16_3=Matrice16_p_3,matrice13_1=Matrice13_p_1,matrice13_2=Matrice13_p_2,matrice13_3=Matrice13_p_3,ed=Matrice_e):
     """ This function samples the energy deposited by a x or gamma rays in the scintillator using response calculated by the Monte-Carlo code MCNP6. 
     
     Parameters
@@ -1328,7 +1340,7 @@ def transf_name(rad):
 
 
 
-def readEShape(rad, *, z=z_ensdf):
+def read_ENSDF(rad, *, z=z_ensdf):
     """ This function reads the ENSDF zip files and format the data to be processed by TDCRPy. 
 
     Parameters
@@ -1348,7 +1360,10 @@ def readEShape(rad, *, z=z_ensdf):
         comprise all transtion probabilities of the daughter nucleus.
     Type : list
         comprise all type of transition of the daughter nucleus.
-
+    Prob_augerK :list
+        comprise all probabilities of electron Auger KLL, KLM, KXY. 
+    Type_augerK : list
+        comprise all types of electron Auger KLL, KLM, KXY.
     """
        
     name = rad + '.txt'
@@ -1380,17 +1395,21 @@ def readEShape(rad, *, z=z_ensdf):
         if 'P' in p:
             index_end.append(i)
             posi.append(i)
-
+        
     Energy = []           # enregistrer les résultats (énergie) complètes
     energy = []           # enregistrer les résultats (énergie) d'une fille
     Type = []             # enregistrer les résultats (type de transition) complètes
     type_ = []            # enregistrer les résultats (type de transition) d'une fille
     Prob = []             # enregistrer les résultats (proba de transition) complètes
     prob = []             # enregistrer les résultats (proba de transition) d'une fille
-    incertitude = []
-    Incertitude = []
-    prob_str = []
-    Prob_str = []
+    incertitude = []      # enregistrer les résultats (incertitude) d'une fille 
+    Incertitude = []      # enregistrer les résultats (incertitude) complètes
+    prob_str = []         # enregistrer les résultats (proba en str) d'une fille
+    Prob_str = []         # enregistrer les résultats (proba en str) complètes
+    Prob_augerK= []       # enregistrer les résultats (prob d'électron auger K) complètes
+    prob_augerK = []      # enregistrer les résultats (prob d'électron auger K) d'une fille
+    Type_augerK = []
+    type_augerK = []
     
     for i in range(len(posi)-1):
         start = posi[i]+1
@@ -1399,8 +1418,10 @@ def readEShape(rad, *, z=z_ensdf):
         e = []                # enregistrer les résultats (énergie) d'un bloc
         prob_b = []           # enregistrer les résultats (proba) d'un bloc
         type_b = []           # enregistrer les résultats (type) d'un bloc
-        incertitude_b = []
-        prob_str_b = []
+        incertitude_b = []    # enregistrer les résultats (incertitude) d'un bloc
+        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 = []
         
         if start==end:        # sauter les lignes blaches et continues
             continue
@@ -1418,45 +1439,52 @@ def readEShape(rad, *, z=z_ensdf):
                     e.append(p1[2])
                     type_b.append(p1[-2])
                     incertitude_b.append(int(p1[-3]))
-                    #incertitude_b.append(p1[-3])
                     prob_str_b.append(p1[3])
                 continue 
-            elif '|]' in p1:                # traiter un bloc qui comprend |]
-                if len(p1)>6:               # repérer la ligne qui comprend la proba
-                    prob_str_b.append(p1[4])
-                    prob_b.append(float(p1[4]))
-                    incertitude_b.append(int(p1[5]))
-                    #incertitude_b.append(p1[5])
-                e.append(float(p1[2]))      # enregistrer les valeurs d'énergie
-                if 'AUGER' in p1:           # traiter le cas d'auger et |] 
-                    if 'K' in p1[-2]:       # Auger K
-                        type_b.append('Auger K')
-                    else:
-                        print('erreur')
-                elif 'X' in p1[-1]:         # Rayon X
-                    type_b.append(p1[-1][0:3])
+            elif '|]' in p1:                   # traiter un bloc qui comprend |]
+                if 'AUGER' in p1:              # block of electron Auger 
+                    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])
+                        prob_b.append(float(p1[5]))        # enregistrer la proba totale du bloc
+                        incertitude_b.append(int(p1[6]))   # enregistrer l'incertitude pour la proba totale
+                        prob_augerK_b.append(float(p1[4])) # enregistrer la proba d'un type d'électron auger associé cette ligne 
+                        type_augerK_b.append(p1[-2])
+                    elif len(p1)>6 and 'K' in p1[-2]:      # la ligne sans proba totale mais ayant une proba pour un type d'électron Auger K (KLL,KLM ou KXY)
+                        prob_augerK_b.append(float(p1[4]))
+                        type_augerK_b.append(p1[-2])
+                elif 'X' in p1[-1]:        # le bloc pour rayon X
+                    if len(p1)>6:          # la ligne  avec la proba totale
+                        prob_str_b.append(p1[4])   
+                        prob_b.append(float(p1[4]))      # enregistrer la proba totale
+                        incertitude_b.append(int(p1[5])) # enregistrer l'incetitude pour la proba totale
+                        type_b.append(p1[-1][0:3])       # enregistrer le type de rayon X
+                else:print('erreur')   
+                e.append(float(p1[2]))      # enregistrer les valeurs d'énergie  
             else:                           # traiter le cas sans |] ni (total)
                 if len(p1)==4 and 'X' in p1[-1]:   # le cas de rayon X sans |] ni (total) ni proba
                     continue                       # sauter cette ligne
                 elif len(p1)==5 and 'L' in p1:     # le cas de Auger L sans |] ni (total) ni proba
                     continue                       # sauter cette ligne
-                else:                       # traiter le cas sans |] ni (total) mais complet
+                else:                         # traiter le cas sans |] ni (total) mais complet
                     e.append(float(p1[2]))         # enregistrer énergie
                     prob_b.append(float(p1[3]))    # enregistrer proba
                     prob_str_b.append(p1[3])
-                    incertitude_b.append(int(p1[4]))
-                    #incertitude_b.append(p1[4])
+                    incertitude_b.append(int(p1[4])) # enregistrer l'incertitude
                     if 'L' in p1:
                         type_b.append('Auger L')   # enregistrer type Auger L
                     else:
                         type_b.append(p1[-1][0:3]) # enregistrer type Rayon X
-
+        
         if len(prob_b)==1 and len(e)>1:            # calculer la valeur moyenne et l'enregistrer au cas où |] compris et valeurs complètes
             energy.append(np.mean(e))
             prob.append(prob_b[0])
             type_.append(type_b[0])
             incertitude.append(incertitude_b[0])
             prob_str.append(prob_str_b[0])
+            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
         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])
@@ -1464,18 +1492,27 @@ def readEShape(rad, *, z=z_ensdf):
                 prob_str.append(prob_str_b[i])
                 type_.append(type_b[i])
                 incertitude.append(incertitude_b[i])
-        if end in index_end or end+1 in index_end: # enregistrer les résultats à la fin d'une fille
+            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
+        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)
             Type.append(type_)
-            Incertitude.append(incertitude)
+            Incertitude.append(incertitude)        
             Prob_str.append(prob_str)
+            Prob_augerK.append(prob_augerK)
+            Type_augerK.append(type_augerK)
             energy = []
             prob = []
             type_ = []   
             incertitude = []
             prob_str = []
-    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  
+
 
 
 def incer(prob,incer):
@@ -1486,12 +1523,12 @@ def incer(prob,incer):
     prob : list of str
         probability (str) of rayon X and Auger electron.
     incer : list of int
-        uncertainty of the probability.
+        standard uncertainty of the probability.
 
     Returns
     -------
     incertitude : list of float
-        DESCRIPTION.
+        absolute uncertainty of proba
 
     '''
     incertitude = []
@@ -1544,15 +1581,16 @@ def relaxation_atom(daugther,rad,lacune='defaut',uncData=False):
     Energy : corresponding energy in keV.
 
     """
-    daug_name,Energy,Prob,Type,Incertitude,f = readEShape(rad)  # tirer les vecteurs de rad d'Ensdf 
-    incertitude = incer(f,Incertitude)
+    daug_name,Energy,Prob,Type,Incertitude,prob_str,Prob_K,Type_K = 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
     
     Energie = np.array(Energy[index_daug])                  # tirer le vecteur d'énergie
     probability = np.array(Prob[index_daug])                # tirer le vecteur de proba
     type_transi = Type[index_daug]                # tirer le vecteur de type
-
+    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])
     
     if len(probability) > 0:                      # le cas où le vecteur de proba/energie/type n'est pas vide
@@ -1610,6 +1648,12 @@ def relaxation_atom(daugther,rad,lacune='defaut',uncData=False):
                 index_fin = sampling(prob_2)        # sample in probability of transition
                 type_fin = type_2[index_fin]        # type of transition     
                 energie_fin = energy_2[index_fin]   # energy of the transition
+                if 'Auger K' == type_fin:
+                    prob_AugerK = normalise(prob_augerK)
+                    index_K = sampling(prob_AugerK)
+                    type_K = type_augerK[index_K]
+                    type_fin = 'Auger K' + type_K[1:]
+                
             else:
                 # print("pas de transition de rayon X ni d'électron Auger pour cette lacune: ",lacune)
                 type_fin = 'NON'