script for convert .t3pa files to .t3pa_cls
.t3pa file example:
Index Matrix Index ToA ToT FToA Overflow
0 4574 832 29 6 0
1 4831 832 35 7 0
2 4575 832 100 8 0
3 31031 1745 22 11 0
.
.
.
.t3pa_cls file example:
% Index Matrix Index [ RowNo, ClmNo ] ToA FToA ( ToA_in_ns ) ToT ( ToT_in_keV ) Overflow
# 1, Nunmasked = 3, Nmasked = 0, Ntot = 3 # Tfirst = 2.0787500000000000e+04 ns, Tlast = 2.0790625000000000e+04 ns, dT = 3.125000 ns, Etot = 64.428148 keV
2 4575 [ 17, 223 ] 832 8 ( 2.0787500000000000e+04 ns ) 100 ( 37.867914 keV ) 0
1 4831 [ 18, 223 ] 832 7 ( 2.0789062500000000e+04 ns ) 35 ( 14.733453 keV ) 0
0 4574 [ 17, 222 ] 832 6 ( 2.0790625000000000e+04 ns ) 29 ( 11.826781 keV ) 0
# 2, Nunmasked = 3, Nmasked = 0, Ntot = 3 # Tfirst = 4.3601562500000000e+04 ns, Tlast = 4.3607812500000000e+04 ns, dT = 6.250000 ns, Etot = 63.577435 keV
5 30775 [ 120, 55 ] 1745 15 ( 4.3601562500000000e+04 ns ) 99 ( 37.617059 keV ) 0
4 30776 [ 120, 56 ] 1745 13 ( 4.3604687500000000e+04 ns ) 44 ( 14.715446 keV ) 0
3 31031 [ 121, 55 ] 1745 11 ( 4.3607812500000000e+04 ns ) 2 2 ( 11.244929 keV ) 0
.
.
.
import numpy as np
import math
#import pandas as pd
import matplotlib.pyplot as plt
import matplotlib.cm as cm
#from matplotlib.mlab import griddata
from urllib.error import HTTPError # recognise the error stemming from missing data
#import urllib
import urllib.request
t3pa2cls_XII - upravena fce energy(a, b, c, t, ToT, pocet_udalosti, RowNo, ClmNo) - nyni je se pocita i s pripadem "nan" t3pa2cls_XV_pc - zkousim vyzobat vysoke energie (jednotlive interakce) - funkce single_interaction. Dale delam prumernou velikost stopy stopy interakce pro danou energii - funkce size_of_interactions_average. Dale delam spektra 2 casti vyboje podle zadaneho casu - primarne pro double breakdown, tj. funkce energy_spectra_doublebreakdown
#Define an exception which will be raised if the data is missing and stop the notebook execution
class StopExecution(Exception):
def _render_traceback_(self):
pass
#shot_no = 44395 #test discharge for which the notebook will definitely work
shot_no = 45912
shot = shot_no
identifier='H03-W0051_shot_'+str(shot)+'_450V'
detector = 'H03-W0051'
ds = np.DataSource('/tmp') # temporary storage for downloaded files
scalars_URL = 'http://golem.fjfi.cvut.cz/shots/{shot_no}/Diagnostics/PlasmaDetection/Results/{name}'
def get_scalar(shot_no, name):
return float(ds.open(scalars_URL.format(shot_no=shot_no, name=name)).read())
t_plasma_start = get_scalar(shot_no, 't_plasma_start')
t_plasma_end = get_scalar(shot_no, 't_plasma_end')
is_plasma = get_scalar(shot_no, 'b_plasma')
def get_file(shot, identifier):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/H03/{identifier}.t3pa'
url = URL.format(shot=shot, identifier=identifier)
try:
file_name_t3pa=url
with urllib.request.urlopen(file_name_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
ft3pa.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_name_t3pa
def get_file_calib(name_calib):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/calib_matrix_H03/{name_calib}.txt'
url = URL.format(shot=shot, name_calib=name_calib)
#print(url)
try:
file_calib=url
with urllib.request.urlopen(file_calib) as calib:
line = calib.readline()
line = line.decode('utf‐8')
calib.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_calib
def load_calib(file_calib):
with urllib.request.urlopen(file_calib) as fc:
calib=[] #vytvoreni 1D pole
for i in range(0,256): #tj. rozsah 0-255
temp = [] # docasne pole
for j in range(0,256):
temp.append(0) #naplneni docasneho pole 0
calib.append(temp) #naplneni pole a[] docasnym polem temp
for i in range(0,256): #nacteni calib matice do pole calib
line = fc.readline()
line = line.decode('utf‐8')
word=line.strip().split(' ')
for j in range(0,256):
#calib[i][j]=float(word[j]) #i = radek, j = sloupec0
calib[j][i]=float(word[j]) #j = radek, i = sloupec0 - pouze pro stavajici kalibraci - verze XV
fc.close
return calib
def load_t3pa_file(file_t3pa):
index=[]
matrix_index=[]
ToA=[]
ToT=[]
FToA=[]
overflow=[]
pocet_udalosti = 0
with urllib.request.urlopen(file_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
while True:
line = ft3pa.readline()
line = line.decode('utf‐8')
word=line.strip().split('\t') #v t3pa souboru je oddelovac \t
if line == '':
break
index.append(word[0])
matrix_index.append(word[1])
ToA.append(float(word[2]))
ToT.append(float(word[3]))
FToA.append(float(word[4]))
overflow.append(float(word[5]))
pocet_udalosti = pocet_udalosti + 1
ft3pa.close
return index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti
def noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti): #tuto fci nemus9m explicitn2 volat - volam ji v fci load_t3pa
pocet=int(0) #pocet sumicich pixelu
konst=int(len(index)/1000)+1
noise_matrix_index=[]
for i in range(0,konst):
pom = [] # pomocne pole
k=0 #pomocna promenna - udava, kolik je v czklu ve skutecnosti udalosti - aby nebyla chyba 'list index out of range'
for j in range(0,1001):
if i*1000+j>=len(index):
break
pom.append(matrix_index[i*1000+j])
k=k+1
for m in range(0,k):
count=int(0) #pocet vvyskytu stejneho matrix index behem 1000 udalosti
index_=int(-1) #budu testovat, jestli pixel na ktery koukam je sumici (abych ho nezapocital 2x)
for p in range(0,pocet):
#index=int(p)
if pom[m]==noise_matrix_index[p]:
index_=p #pixel na ktery jsem uz koukal a byl sumici
break
if index_ >=0 and pom[m]==noise_matrix_index[index_]:
continue
for l in range(0,k):
if pom[m]==pom[l]:
count=count+1
####podminka na sumici pixely
if count>=50: #kdyz se pixel vyskytne behem tisice udalosti vicekrat nez toto cislo, je sumici
noise_matrix_index.append(pom[m])
#noise_matrix_index[pocet]=pom[i]
pocet=pocet+1
pom.clear()
pocet_udalosti=len(index)
for n in range (0,pocet_udalosti):
for o in range(0,len(noise_matrix_index)):
if n >=pocet_udalosti:
break
if(matrix_index[n]==noise_matrix_index[o]):
del matrix_index[n]
del index[n]
del ToA[n]
del ToT[n]
del FToA[n]
del overflow[n]
pocet_udalosti=pocet_udalosti-1
continue
return pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow
def t3pa_data(pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow):
#rovnou vyhodim sumici pixely
pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow=noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti)
RowNo=[]
ClmNo=[]
for i in range(0,len(matrix_index)):
#RowNo.append(int(int(matrix_index[i]))//int(256))
#ClmNo.append(int(int(matrix_index[i]))%int(256))
ClmNo.append(int(int(matrix_index[i]))//int(256)) #ver XV - bude pro novou kalibraci - nyn9 to bere 3patnou kalibraci jednotlivych pixelu (cislovani radku a sloupcu stejne jako v pixetu)
RowNo.append(int(int(matrix_index[i]))%int(256)) ##ver XV - bude super pro novou kalibraci - odpovida radkum a sloupcum v pixetu. Nyni bere spatne kalibrace pixelu (viz fce load_calib)
return index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti, RowNo, ClmNo
def hit_map(detector,hit_map_fig,RowNo,ClmNo):
plt.hist2d(RowNo,ClmNo,bins=(256,256),cmap='Blues')
cb=plt.colorbar()
cb.set_label('Counts in pixel')
plt.xlabel('x [pixel]')
plt.ylabel('y [pixel]')
plt.title(detector)
plt.savefig(hit_map_fig, dpi = 1000)
return
def energy(a, b, c, t, ToT, pocet_udalosti, RowNo, ClmNo):
E=[] #energy in keV
#for i in range (0,pocet_udalosti):
pom=0
for i in range (0,len(ToT)):
Sqrt=float(0.0)
e1=float(0.0)
e2=float(0.0)
# promenna sqrt je vnitrek odmocniny
Sqrt = (((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i])))*(((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i])))) + (float(4)*float(a[RowNo[i]][ClmNo[i]])*float(c[RowNo[i]][ClmNo[i]]))) #zmena oproti verzi VI
if float(Sqrt)<float(0):
E.append(float(0))
else:
'''
V kalibracni matici a se obcas vyskytne 0 -> ve vypoctu energie
je tim padem deleni nulou -> energie diverguje. Jak to vyresit?
zatim polozim energii = 0 (kdyz a=0), pak se uvidi
nakonec udelam limitu vyrazu energie pro a->0 (L'hopital)
'''
if a[RowNo[i]][ClmNo[i]]==0:
e1=((float(t[RowNo[i]][ClmNo[i]]))/float(2)) + ((((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i]))*(float(t[RowNo[i]][ClmNo[i]]))) - 2*(float(c[RowNo[i]][ClmNo[i]])))/(float(2)*np.sqrt(float(Sqrt))))
e2=((float(t[RowNo[i]][ClmNo[i]]))/float(2)) - ((((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i]))*(float(t[RowNo[i]][ClmNo[i]]))) - 2*(float(c[RowNo[i]][ClmNo[i]])))/(float(2)*np.sqrt(float(Sqrt))))
else:
e1=((-(float(b[RowNo[i]][ClmNo[i]]) - (float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]]))-float(ToT[i])))+np.sqrt(float(Sqrt)))/(float(2)*float(a[RowNo[i]][ClmNo[i]]))
e2=((-(float(b[RowNo[i]][ClmNo[i]]) - (float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]]))-float(ToT[i])))-np.sqrt(float(Sqrt)))/(float(2)*float(a[RowNo[i]][ClmNo[i]]))
if a[RowNo[i]][ClmNo[i]]<0:
e1=-1
e2=-1
if math.isnan(e1):
e1=-1
if math.isnan(e2):
e2=-1
if e1<0 and e2<0:
E.append(float(0))
if e1>=0 and e1>e2:
E.append(float(e1))
if e2>=0 and e2>e1:
E.append(float(e2))
if e1>=0 and e2==e1:
E.append(float(e1))
return E
def Time(ToA, FToA, pocet_udalosti, RowNo, ClmNo):
T=[] #time in ns
for i in range (0,pocet_udalosti):
Time=float(0.0)
Time=(float(ToA[i])-((float(FToA[i])/float(16))))*float(25)
T.append(float(Time))
return T
def Timewalk_parameters_Si():
#SI - korekce na TimeWalk - parametry
A = -1.21988
B = 4.33638
C = 29.5075
D = 1
sigma_A=0.7013
sigma_B=0.1366
sigma_C=4.753
sigma_D=0
return A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D
def Timewalk(E,T):
i=0
A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
while i < len(T):
timewalk=float(0)
timewalk = (C / (E[i] - B)) + A
if(timewalk<0):
timewalk=0
T[i]=T[i]-timewalk
i=i+1
return E,T
def remove_interactions_with_zero_energy(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T):
i=0
treshold=5.015347
while i < len(T):
if E[i]<treshold: #E[i] < energy treshold
index.pop(i)
matrix_index.pop(i)
ToA.pop(i)
ToT.pop(i)
FToA.pop(i)
overflow.pop(i)
RowNo.pop(i)
ClmNo.pop(i)
E.pop(i)
T.pop(i)
continue
i=i+1
return index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T
def clustering_new(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T):
dT=float(50)
indexCl, TCl,ECl, matrix_indexCl, ToACl,ToTCl,FToACl,RowNoCl,ClmNoCl,overflowCl=[],[],[],[],[],[],[],[],[],[]
StartLastElem=len(T)-1
indexCl.append(int(index[StartLastElem]))
TCl.append(float(T[StartLastElem]))
ECl.append(float(E[StartLastElem]))
matrix_indexCl.append(int(matrix_index[StartLastElem]))
RowNoCl.append(int(RowNo[StartLastElem]))
ClmNoCl.append(int(ClmNo[StartLastElem]))
ToACl.append(float(ToA[StartLastElem]))
ToTCl.append(float(ToT[StartLastElem]))
FToACl.append(float(FToA[StartLastElem]))
overflowCl.append(float(overflow[StartLastElem]))
del index[StartLastElem]
del T[StartLastElem]
del E[StartLastElem]
del matrix_index[StartLastElem]
del RowNo[StartLastElem]
del ClmNo[StartLastElem]
del ToA[StartLastElem]
del ToT[StartLastElem]
del FToA[StartLastElem]
del overflow[StartLastElem]
j=1
pom=float(TCl[0]+dT)
while(j >0):
if(len(T) == 0):
break
k=0
j=0
while (k<=(len(TCl)-1)):
i=len(T)-1
if(len(T) == 0):
break
pocet_sousedu=0 #pocet sousednich pixelu - mohou byt maximalne 4
delka=0
# verze X
count=0 #pomocna promanna, kterou urcuji, ze se ma nasledujici cyklus while projit jeste jednou, pokud je i = -1
while(float(T[i])<=(pom)):
delka=delka+1
if(((((int(RowNoCl[k]))==(int(RowNo[i])+1))or((int(RowNoCl[k]))==(int(RowNo[i])-1))) and ((int(ClmNoCl[k]))==(int(ClmNo[i])))) or (((int(RowNoCl[k]))==(int(RowNo[i]))) and (((int(ClmNoCl[k]))==(int(ClmNo[i])+1))or((int(ClmNoCl[k]))==(int(ClmNo[i])-1))))):
#beru jen pixely, které mají společnou jednu stranu.
#pixely, kter0 spolu sousedí přes roh neuvažuji
indexCl.append(int(index[i]))
TCl.append(float(T[i]))
ECl.append(float(E[i]))
matrix_indexCl.append(int(matrix_index[i]))
RowNoCl.append(int(RowNo[i]))
ClmNoCl.append(int(ClmNo[i]))
ToACl.append(float(ToA[i]))
ToTCl.append(float(ToT[i]))
FToACl.append(float(FToA[i]))
overflowCl.append(float(overflow[i]))
# Removes i-th Row
del index[i]
del T[i]
del E[i]
del matrix_index[i]
del RowNo[i]
del ClmNo[i]
del ToA[i]
del ToT[i]
del FToA[i]
del overflow[i]
j=j+1
i=len(T)-1
pocet_sousedu=pocet_sousedu+1
if(len(T) == 0):
break
if(pocet_sousedu==4):
break
continue
i=i-1
if(i==-1): # verze X
count=count+1
if(i<0 and len(T)>0): # verze X
i=0
if(count>1):
break
if(i>=len(T)):
break
k=k+1
if(len(TCl)>2):
indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl = insertionSort(indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl)
return T, indexCl,TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl
def insertionSort(indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl):
# Function to do insertion sort
# Traverse through 1 to len(arr)
for i in range(1, len(TCl)):
key = TCl[i]
# Move elements of arr[0..i-1], that are
# greater than key, to one position ahead
# of their current position
#ostatni
key1 = indexCl[i]
key2 = ECl[i]
key3 = matrix_indexCl[i]
key4 = RowNoCl[i]
key5 = ClmNoCl[i]
key6 = ToACl[i]
key7 = ToTCl[i]
key8 = FToACl[i]
key9 = overflowCl[i]
j = i-1
while j >= 0 and key < TCl[j] :
TCl[j + 1] = TCl[j]
#ostatni
indexCl[j + 1] = indexCl[j]
ECl[j + 1] = ECl[j]
matrix_indexCl[j + 1] = matrix_indexCl[j]
RowNoCl[j + 1] = RowNoCl[j]
ClmNoCl[j + 1] = ClmNoCl[j]
ToACl[j + 1] = ToACl[j]
ToTCl[j + 1] = ToTCl[j]
FToACl[j + 1] = FToACl[j]
overflowCl[j + 1] = overflowCl[j]
j -= 1
TCl[j + 1] = key
#ostatni
indexCl[j + 1] = key1
ECl[j + 1] = key2
matrix_indexCl[j + 1] = key3
RowNoCl[j + 1] =key4
ClmNoCl[j + 1] = key5
ToACl[j + 1] = key6
ToTCl[j + 1] = key7
FToACl[j + 1] = key8
overflowCl [j + 1] = key9
return indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl
def file_t3pa_cls_new(file_t3pa_cls,T):
with open(file_t3pa_cls, "w", encoding="utf-8") as t3pa_cls:
t3pa_cls.write('%\n')
t3pa_cls.write('% Index Matrix Index [ RowNo, ClmNo ] ToA FToA ( ToA_in_ns ) ToT ( ToT_in_keV ) Overflow\n')
t3pa_cls.write('\n')
i=1
T_first=[]
E_tot=[]
N_tot=[]
eventNo=0 #for high_energy_event function
while(len(T) > 0):
T, indexCl,TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl = clustering_new(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T)
Tfirst=float(TCl[0])
Tlast=float(TCl[len(TCl)-1])
dT=Tlast-Tfirst
Etot=float(0)
for k in range(0,len(TCl)):
Etot=Etot+float(ECl[k])
T_first.append(float(Tfirst))
dT=Tlast-Tfirst
E_tot.append(float(Etot))
N_tot.append(len(TCl)) #new in ver. XV
t3pa_cls.write('# '+str(i)+', Nunmasked = '+str(len(TCl))+', Nmasked = 0, Ntot = '+str(len(TCl))+'\n')
t3pa_cls.write('# Tfirst = '+str(Tfirst)+' ns, Tlast = '+str(Tlast)+' ns, dT = '+str(dT)+' ns, Etot = '+str(Etot)+' keV\n')
for j in range(0,len(TCl)):
t3pa_cls.write(str(indexCl[j])+' '+str(matrix_indexCl[j])+' [ '+str(RowNoCl[j])+', '+str(ClmNoCl[j])+' ] '+str(ToACl[j])+' '+str(FToACl[j])+' ( '+str(TCl[j])+' ns ) '+str(ToTCl[j])+' ( '+str(ECl[j])+' keV ) '+str(overflowCl[j])+'\n')
t3pa_cls.write('\n')
#eventNo=single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path)
i=i+1
t3pa_cls.close
return T_first, E_tot, N_tot
def energy_spectrum_in_time(Tfirst, Etot): #dela histogram - energie zaznamenana v case
pom = 0
dt=100 #(ns) time width of 1 bin
T_first=0 #cas, kdy prisel trigger a yacalo mereni
T_last=(max(Tfirst)) #posledni z Tfirst
Delta_T = T_last - T_first
poc = int(int(Delta_T) / float(dt)) + 1 #pocet casovych oken
T_int_first=[] #cas
E=[] #energie
for i in range(0,poc):
T_int_first.append((i*dt) + dt/2)
E.append(0)
#XII
for j in range(0,len(Tfirst)):
time_index=0
time_index=int(((Tfirst[j]-T_first)/dt))
if float(Tfirst[j]-T_first) >= (T_int_first[time_index] - dt / 2) and float(Tfirst[j]-T_first) < (T_int_first[time_index] + dt / 2):
E[time_index]=float(E[time_index])+float(Etot[j])
pom=pom+1
for l in range(0,len(T_int_first)):
T_int_first[l]=T_int_first[l]+T_first
caption, T_int_first = energy_in_time_hist(T_int_first, E, figure_E_in_time_hist, t_plasma_start, t_plasma_end, is_plasma, dt)
return dt, caption, T_int_first, E
def energy_in_time_hist(T_int_first, E,figure_E_in_time_hist, t_plasma_start, t_plasma_end, is_plasma, dt):
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
for k in range(0,len(T_int_first)):
T_int_first[k] = T_int_first[k] / 1000000
plt.plot(T_int_first, E, drawstyle='steps-mid')
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Time [ms]')
plt.ylabel('Energy [keV]')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
plt.axvline(t, color='k', linestyle='--')
plt.xlim([0, (t_plasma_start + t_plasma_end)])
else:
plt.xlim(0,)
plt.ylim(0,) #10 000 keV
plt.savefig(figure_E_in_time_hist, dpi = 1000)
caption = '# x = time in ms, count = energy in keV, dT= '+str(dt)+' ns'
return caption, T_int_first
def hits_in_time_hist_new(T, dt, t_plasma_start, t_plasma_end, is_plasma,figure_count_in_time_hist):
pom = 0
T_first=0 #cas, kdy prisel trigger a yacalo mereni
T_last=(max(T)) #posledni z Tfirst
Delta_T = T_last - T_first
poc = int(int(Delta_T) / float(dt)) + 1 #pocet casovych oken
T_hit=[] #cas
count=[] #energie
for i in range(0,poc):
T_hit.append((i*dt) + dt/2)
count.append(0)
for j in range(0,len(T)):
time_index=0
time_index=int(((T[j]-T_first)/dt))
k=time_index
for j in range(0,len(T)):
time_index=0
time_index=int(((T[j]-T_first)/dt))
if float(T[j]-T_first) >= (T_hit[time_index] - dt / 2) and float(T[j]-T_first) < (T_hit[time_index] + dt / 2):
count[time_index] = count[time_index] + 1
pom=pom+1
for l in range(0,len(T_hit)):
T_hit[l]=T_hit[l]+T_first
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
for k in range(0,len(T_hit)):
T_hit[k] = T_hit[k] / 1000000
#plt.plot(T_hit, count)
plt.plot(T_hit, count, drawstyle='steps-mid')
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Time [ms]')
plt.ylabel('Count')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
plt.axvline(t, color='k', linestyle='--')
plt.xlim([0, (t_plasma_start + t_plasma_end)])
else:
plt.xlim(0,)
plt.ylim(0,) #10 000 keV
plt.savefig(figure_count_in_time_hist, dpi = 1000)
caption = '# x = time in ms, dT= '+str(dt)+' ns'
return caption, T_hit,count
def energy_spectrum(Etot):
E_min=0
dE=5 #keV
E_max=max(Etot)
pocet=(E_max//dE) + 3
pocet=int(pocet)
E_max=float(dE*pocet)
xle=[]
xre=[]
xmean=[]
for p in range (0,pocet):
xle.append(E_min + (p * (E_max - E_min)) / pocet)
xre.append(xle[p]+dE)
xmean.append((xle[p] + xre[p]) / 2)
count=[]
for l in range(0,pocet):
count.append(0)
#XII
for i in range(0,len(Etot)):
E_index=int(((Etot[i]-E_min)/dE))
if ((xle[E_index] <= Etot[i]) and (Etot[i] < xre[E_index])):
count[E_index]=count[E_index]+1
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
#ax.hist(Etot, bins = xle)
plt.fill_between(xmean, count,step='mid') #pre mid post
plt.plot(xmean, count, drawstyle='steps-mid')
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Count')
plt.xlim(0,)
ax.set_yscale('log') #log scale y
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
plt.savefig(figure_E_hist, dpi = 1000)
return caption, xmean,count, xle, Etot
def hist_file(file_hist, xmean, count, caption ):
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('#\n')
hist.write('#'+str(caption)+'\n')
hist.write('# x_mean count\n')
hist.write('\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(count[m])+'\n')
hist.close
return T_first, E_tot
def multiplot(icon_fig, x1,y1,x2,y2):
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(nrows=2,figsize =(10, 7))
ax[0].plot(x1, y1, drawstyle='steps-mid')
ax[0].set_xlabel('Time [ms]')
ax[0].set_ylabel('Energy [keV]')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
ax[0].axvline(t, color='k', linestyle='--')
ax[0].set_xlim([0, (t_plasma_start + t_plasma_end)])
else:
ax[0].set_xlim(0,)
ax[0].set_ylim(0,) #keV
#ax[1].hist(y2, bins = x2)
ax[1].fill_between(x2, y2,step='mid') #pre mid post
ax[1].plot(x2, y2, drawstyle='steps-mid')
ax[1].set_xlabel('Energy [keV]')
ax[1].set_ylabel('Count')
ax[1].set_xlim(0,)
#ax[1].set_ylim(0,)
ax[1].set_yscale('log') #log scale y
fig.subplots_adjust(hspace=0.3)
plt.savefig(icon_fig, dpi = 1000)
return
def single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path): # ver XV
#tuto funkci volam ve funkci file_t3pa_cls_new(file_t3pa_cls,T):
#This function is calling in function file_t3pa_cls_new(file_t3pa_cls,T)
high_E = 1500 #[keV] - zadani, jaka energie je pro me vysoka
if(Etot>=high_E or len(ECl)>15):
#if (len(ECl)>15): #len(ECl) = pocet pixelu dane interakce
print("Total energy of track: "+str(Etot)+" keV\n")
A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
######################
#SI - rekonstrukce z(t) - parametry
a=0.0193947
b=-0.0273934
sigma_a=0.0004163
sigma_b=0.009832
#######################
time=[]
sigma_time=[]
time_min = []
time_max=[]
z_rec=[]
z_rec_min=[]
z_rec_max=[]
sigma_z_rec=[]
for i in range(0, len(TCl)):
time.append(TCl[i] - Tfirst)
#chyba urceni casu = chyba timewalk ->
sigma_time.append(np.sqrt(((sigma_C / (pow((ECl[i] - B), D)))*(sigma_C / (pow((ECl[i] - B), D)))) + (((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))*((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))) + (((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))* ((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))) + (sigma_A*sigma_A)))
time_min.append(time[i] - sigma_time[i])
time_max.append(time[i] + sigma_time[i])
#//
#//rekonstrukce hloubky interakce
z_rec.append(a * time[i] + (b))
sigma_z_rec.append(np.sqrt(((time[i]*time[i])*(sigma_a*sigma_a)) + (sigma_b*sigma_b) + (a*a)*(sigma_time[i]*sigma_time[i])))
z_rec_min.append(z_rec[i] - sigma_z_rec[i])
z_rec_max.append(z_rec[i] + sigma_z_rec[i])
#//
single_interaction_file=str(path)+'event'+str(eventNo)+'.txt'
single_interaction_2D_fig=str(path)+'2D_event'+str(eventNo)+'.png'
single_interaction_3D_fig=str(path)+'3D_event'+str(eventNo)+'.png'
with open(single_interaction_file, "w", encoding="utf-8") as f:
f.write('#\n')
f.write('# x_mean y_mean z_rec z_rec_min z_rec_max sigma_z_rec\n')
f.write('#E_tot='+str(Etot)+'\n')
f.write('\n')
for m in range(0,len(TCl)):
#print(str(RowNoCl[m])+' '+str(ClmNoCl[m])+' '+str(z_rec[m])+' '+str(z_rec_min[m])+' '+str(z_rec_max[m])+' '+str(sigma_z_rec[m])+' '+str(ECl[m])+'\n')
f.write(str(RowNoCl[m])+' '+str(ClmNoCl[m])+' '+str(z_rec[m])+' '+str(z_rec_min[m])+' '+str(z_rec_max[m])+' '+str(sigma_z_rec[m])+' '+str(ECl[m])+'\n')
f.close
############################### plot
#binwidth=1
if ((max(RowNoCl)-min(RowNoCl))>=(max(ClmNoCl)-min(ClmNoCl))):
rang=max(RowNoCl)-min(RowNoCl)
#bins=(range((min(RowNoCl)-1), (max(RowNoCl)+3), binwidth), range((min(ClmNoCl)-1), (min(ClmNoCl)+rang+3), binwidth))
x_min=min(RowNoCl)-1
x_max=max(RowNoCl)+2
y_min=min(ClmNoCl)-1
y_max=min(ClmNoCl)+rang+2
else:
rang = max(ClmNoCl)-min(ClmNoCl)
#bins=(range((min(RowNoCl)-1), (min(RowNoCl)+rang+3), binwidth), range((min(ClmNoCl)-1), (max(ClmNoCl)+3), binwidth))
x_min=min(RowNoCl)-1
x_max=min(RowNoCl)+rang+2
y_min=min(ClmNoCl)-1
y_max=max(ClmNoCl)+2
###############$$$$$$$$$$$$$$$$$$$$
rows, cols = (x_max-x_min, x_max-x_min)
#arrECl = [[0]*cols]*rows
arrECl= [[0] * cols for _ in range(rows)]
for i in range (x_min,x_max):
for j in range (y_min,y_max):
for k in range(0,len(RowNoCl)):
if ((i==RowNoCl[k]) and (j==ClmNoCl[k])):
#print(str(i)+" "+str(j)+" "+str(RowNoCl[k])+" "+str(ClmNoCl[k])+"\n")
arrECl[i-x_min][j-y_min]=ECl[k]
arrECl=np.transpose(arrECl)
arrECl[arrECl == 0] = np.nan
#2D plot
plt.imshow(arrECl,cmap='jet', interpolation='none', origin=[(x_min-0.2),(y_min-0.2)], aspect='equal', extent=[x_min-0.2,x_max,y_min-0.2,y_max])
###############$$$$$$$$$$$$$$$$$$$$
cb=plt.colorbar(pad=0.005)
plt.tick_params(axis="x",direction="in")
plt.tick_params(axis="y",direction="in")
cb.ax.tick_params(axis='y', direction='in')
cb.set_label('Energy in pixel [keV]')
plt.xlabel('x [pixel]')
plt.ylabel('y [pixel]')
plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
plt.savefig(single_interaction_2D_fig, dpi = 1000)
plt.show()
plt.close()
############################### plot
############## 3D plot
# Creating figure
fig = plt.figure(figsize = (16, 9))
ax = plt.axes(projection ="3d")
# Add x, y gridlines
ax.grid(b = True, color ='grey', #
linestyle ='-.', linewidth = 0.3,
alpha = 0.2)
# Creating plot
sctt = ax.scatter3D(RowNoCl, ClmNoCl, z_rec,
alpha = 0.8,
c = ECl,
cmap = 'jet',
s = 4000, #size of bubbles in plot
marker ='o')
plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
ax.set_xlabel('x [pixel]', fontweight ='bold')
ax.set_ylabel('y [pixel]', fontweight ='bold')
ax.set_zlabel('z [mm]', fontweight ='bold')
cb=fig.colorbar(sctt, ax = ax, shrink = 0.5,aspect = 10,pad=0.00005) #
cb.ax.tick_params(axis='y', direction='in')
cb.set_label('Energy in pixel [keV]', fontweight ='bold')
plt.savefig(single_interaction_3D_fig, dpi = 1000)
# show plot
plt.show()
plt.close()
############## 3D plot
eventNo=eventNo+1
return eventNo
def size_of_interactions_average(Etot, Ntot,file_hist): #ver XV
E_min=0
dE=5 #keV
E_max=max(Etot)
pocet=(E_max//dE) + 3
pocet=int(pocet)
E_max=float(dE*pocet)
xle=[]
xre=[]
xmean=[]
count=[]
nop=[] #number of pixels
nop_mean=[] #mean number of pixels in dE
sigma=[] #root mean square error of the mean - stredni kvadraticka chyba prumeru
for p in range (0,pocet):
xle.append(E_min + (p * (E_max - E_min)) / pocet)
xre.append(xle[p]+dE)
xmean.append((xle[p] + xre[p]) / 2)
count.append(0)
nop.append(0)
nop_mean.append(0)
sigma.append(0)
#XII
for i in range(0,len(Etot)):
E_index=int(((Etot[i]-E_min)/dE))
if ((xle[E_index] <= Etot[i]) and (Etot[i] < xre[E_index])):
count[E_index]=count[E_index]+1
nop[E_index]=nop[E_index]+Ntot[i] #ver XV
for l in range(0,pocet):
if nop[l]!=0:
nop_mean[l]=float(nop[l])/float(count[l])
#####Computing sigma
for i in range(0,len(Etot)):
E_index=int(((Etot[i]-E_min)/dE))
if ((xle[E_index] <= Etot[i]) and (Etot[i] < xre[E_index])):
if nop_mean[E_index]!=0:
sigma[E_index]=sigma[E_index]+((float(Ntot[i])-float(nop_mean[E_index]))*(float(Ntot[i])-float(nop_mean[E_index])))
for l in range(0,pocet):
if nop[l]!=0 and count[l]>1:
sigma[l]=np.sqrt(sigma[l]/(float(count[l])*(float(count[l])-1)))
if nop[l]!=0 and count[l]==1:
sigma[l]=10
#####Computing sigma
#####write to file
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('# x_mean [keV] mean_size_of_interaction [pixel] sigma_mean_size_of_interaction\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(nop_mean[m])+' '+str(sigma[m])+'\n')
hist.close
#####write to file
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
#plt.fill_between(xmean, nop,step='mid',color='blue', alpha=0.3) #pre mid post
#plt.fill_between(xmean, nop,step='mid') #pre mid post
plt.plot(xmean, nop_mean, drawstyle='steps-mid')
plt.errorbar(xmean, nop_mean, yerr = sigma,fmt ='.', linewidth=0.7)
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Average size of interaction [pixel]')
plt.xlim(0,500)
plt.ylim(0,17)
#ax.set_yscale('log') #log scale y
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
plt.savefig(figure_size_interaction, dpi = 1000)
return
def energy_spectra_doublebreakdown(Etot, Tfirst, file_hist):
#this function separates discharge into two parts acording to time Tsep and creates energy spectrum for each part.
#This is for discharges with double brakdown
Tsep=14636000 # [ns]
E_min=0
dE=5 #keV
E_max=max(Etot)
pocet=(E_max//dE) + 3
pocet=int(pocet)
E_max=float(dE*pocet)
xle=[]
xre=[]
xmean=[]
for p in range (0,pocet):
xle.append(E_min + (p * (E_max - E_min)) / pocet)
xre.append(xle[p]+dE)
xmean.append((xle[p] + xre[p]) / 2)
count_first=[]
count_second=[]
for l in range(0,pocet):
count_first.append(0)
count_second.append(0)
#XII
for i in range(0,len(Etot)):
E_index=int(((Etot[i]-E_min)/dE))
if ((xle[E_index] <= Etot[i]) and (Etot[i] < xre[E_index])):
if(Tfirst[i]<Tsep):
count_first[E_index]=count_first[E_index]+1
else:
count_second[E_index]=count_second[E_index]+1
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('#\n')
hist.write('#'+str(caption)+'\n')
hist.write('# x_mean count_before_2nd_breakdown count_2nd_breakdown\n')
hist.write('\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(count_first[m])+' '+str(count_second[m])+'\n')
hist.close
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
#plt.fill_between(xmean, count_first,step='mid') #pre mid post
plt.plot(xmean, count_first, label='Before 2nd breakdown', drawstyle='steps-mid')
plt.plot(xmean, count_second, label='2nd breakdown', drawstyle='steps-mid')
plt.legend(loc='best')
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Count')
plt.xlim(0,500)
ax.set_yscale('log') #log scale y
plt.savefig(figure_energy_spectra_doublebreakdown, dpi = 1000)
return
#soubory, ktere ctu
#read files
t3pa=get_file(shot, identifier)
name_calib='caliba'
caliba=get_file_calib(name_calib)
name_calib='calibb'
calibb=get_file_calib(name_calib)
name_calib='calibc'
calibc=get_file_calib(name_calib)
name_calib='calibt'
calibt=get_file_calib(name_calib)
#vytvorene soubory:
#created files
t3pa_cls= 'H03-W0051_shot_'+str(shot)+'_450V.t3pa_cls'
E_hist= 'H03-W0051_shot_'+str(shot)+'_450V_E_hist.txt'
E_in_time_hist= 'H03-W0051_shot_'+str(shot)+'_450V_discharge_energy.txt'
count_in_time_hist= 'H03-W0051_shot_'+str(shot)+'_450V_discharge_hits.txt'
size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction.txt'
file_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown.txt'
#created figures
icon_fig='icon-fig'
figure_E_in_time_hist='discharge_energy'
figure_count_in_time_hist='discharge_hits'
figure_E_hist='Energy_spectrum'
hit_map_fig='hit-map'
figure_size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction'
figure_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown'
#nactu jednotlive kalibracni matice - abych to nemusel delat v kazde funkci
a=load_calib(caliba)
b=load_calib(calibb)
c=load_calib(calibc)
t=load_calib(calibt)
#nactu a urcim jednotlive hodnoty - abych to nemusel delat v kazde funkci
index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti = load_t3pa_file(t3pa)
index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti, RowNo, ClmNo = t3pa_data(pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow)
raw data
#hit map
hit_map(detector,hit_map_fig,RowNo,ClmNo)
Energy and time calculation from raw data.
E=energy(a, b, c, t, ToT, pocet_udalosti, RowNo, ClmNo)
T=Time(ToA, FToA, pocet_udalosti, RowNo, ClmNo)
index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T = remove_interactions_with_zero_energy(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T)
E,T=Timewalk(E,T)
#sort by time
T, index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E = (list(t) for t in zip(*sorted(zip(T, index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E), reverse=True))) #serazeni od nejvetsiho po nejmensi
T_pom=T.copy()
#save to file
T_first, E_tot, Ntot = file_t3pa_cls_new(t3pa_cls,T)
dt, caption, T_int_first, E = energy_spectrum_in_time(T_first, E_tot)
hist_file(E_in_time_hist, T_int_first, E, caption);
([4183873.4375, 4200589.0625, 4245645.486060796, 4245648.321076619, 4273347.460452169, 4314272.98618653, 4463036.308541708, 4516373.4375, 4519480.333225946, 4520039.226822126, 4520225.0, 4520815.422263485, 4525556.211428195, 4525936.918065219, 4526217.1875, 4526593.75, 4529406.065605748, 4529410.086621125, 4580491.832632525, 4586693.342413341, 4595065.625, 4596572.697327951, 4604054.484354057, 4711086.660961539, 4722958.058993217, 4810845.3125, 4814384.375, 4891162.5, 6695396.9985546, 6866860.352375621, 6868261.685191906, 6994003.125, 7099131.25, 7129463.278797314, 7130623.073756609, 7130945.010654807, 8185245.3125, 8525754.387882605, 8628820.776659146, 8696595.68875797, 8701679.491570348, 8703296.875, 8703743.628683295, 8705085.336426139, 8705454.924390363, 8706293.053784892, 8706748.4375, 8730331.25, 8731529.476926165, 8748865.30920586, 8749673.4375, 8750856.25, 8779617.1875, 8790217.1875, 8825184.74532947, 8834315.625, 8853562.445120098, 8863321.875, 8864176.697293386, 8886742.76603319, 8892519.140832322, 8904692.1875, 8904913.658927286, 8905095.020341542, 8905349.964492438, 8905650.0, 8918130.442457212, 8919474.243177438, 8922011.282267557, 8928745.3125, 8931389.0625, 8937252.011676298, 8947414.0625, 8952057.307835227, 8965704.685666695, 8968389.0625, 8969540.953686716, 8970240.882425923, 8970963.154934406, 8971078.125, 8975503.980992258, 8975887.5, 8976480.68812859, 8976954.5977973, 8977042.565043097, 8980778.213154148, 8980976.5625, 8980985.456670586, 8981047.939277556, 8981551.095397804, 8981751.29438005, 8981798.4375, 8982078.030210437, 8982135.129054297, 8982614.469613882, 8982882.12367435, 8983096.929754166, 8983256.25, 8983643.557566142, 8983773.4375, 8987637.05927988, 8988124.254688518, 8988876.295442553, 8996076.24896725, 8997781.09847637, 8999818.443383655, 9000604.424873259, 9001122.357356636, 9001125.894374518, 9001335.9375, 9001542.16464767, 9004281.3693377, 9005364.941423027, 9005503.562604766, 9005923.390236454, 9006085.420970289, 9006089.0625, 9006169.007326934, 9006280.44515653, 9006329.517238103, 9006368.950798277, 9006621.719176397, 9006647.152944328, 9006790.625, 9006796.875, 9006884.923797358, 9007031.321908507, 9007034.372597642, 9007046.104904564, 9007155.040208349, 9007167.259587096, 9007182.30970794, 9007227.141970402, 9007312.964309137, 9007439.402320221, 9007560.507849539, 9007670.558947109, 9007743.023725435, 9007760.9375, 9007764.0625, 9007819.531306375, 9008027.550515702, 9008043.839343544, 9008055.380151317, 9008062.5, 9008104.6875, 9008143.505249066, 9008199.465843676, 9008345.143207641, 9008404.6875, 9008411.211793315, 9008517.924231231, 9008560.837918967, 9008656.505350163, 9008714.0625, 9008840.750927528, 9008893.75, 9009074.880919412, 9009133.625789277, 9009182.45264877, 9009201.682209646, 9009321.996323926, 9009331.541482333, 9009363.358914142, 9009468.169689404, 9009491.236132098, 9009671.800921798, 9009751.5625, 9009754.405385567, 9009768.69761471, 9009918.39377025, 9009929.048365531, 9009946.984452512, 9009953.125, 9010026.5625, 9010035.73518393, 9010084.06247413, 9010100.137702378, 9010101.757499326, 9010116.035248041, 9010130.758633042, 9010135.885259666, 9010141.002152685, 9010151.5625, 9010167.15521591, 9010210.9375, 9010212.5, 9010222.428020593, 9010326.869958924, 9010334.384608725, 9010350.248461494, 9010354.845121, 9010377.016465543, 9010404.987460053, 9010425.872245127, 9010455.531909695, 9010478.18827834, 9010510.9375, 9010589.786125397, 9010597.601751687, 9010634.375, 9010642.1875, 9010657.8125, 9010668.75, 9010671.440412492, 9010747.83882921, 9010783.42131511, 9010807.8125, 9010813.912468296, 9010868.673984608, 9010883.92879088, 9010905.421934988, 9010912.5, 9010929.6875, 9010981.25, 9011003.32062525, 9011030.585197488, 9011042.1875, 9011107.8125, 9011150.0, 9011244.019793529, 9011268.424141327, 9011268.631759254, 9011300.524493318, 9011330.166647157, 9011347.483890036, 9011403.785153715, 9011485.086929744, 9011520.473087369, 9011555.373735217, 9011576.543449743, 9011648.38375643, 9011740.095098844, 9011748.940079143, 9011795.3125, 9011848.4375, 9011875.0, 9011901.270979028, 9011943.324927993, 9011951.5625, 9012105.207295107, 9012110.623040995, 9012123.06182649, 9012215.392783487, 9012225.0, 9012557.8125, 9012846.238316586, 9012985.653734386, 9013038.183354102, 9013056.555102466, 9013261.361907978, 9013425.0, 9015011.177805066, 9017948.338090392, 9019175.928783504, 9023421.875, 9024487.5, 9030298.4375, 9061534.749423213, 9063212.96730573, 9063230.76319783, 9063521.57554675, 9063547.998211809, 9063607.018780848, 9063667.1875, 9063702.865812365, 9063825.550099153, 9064021.875, 9064248.59936124, 9064320.3125, 9064390.070606042, 9064467.140083848, 9064715.280235466, 9064726.532983446, 9065384.767233357, 9066142.158094348, 9066198.905409403, 9066262.5, 9066264.0625, 9066506.25, 9066533.510633184, 9066636.528919457, 9066854.911846513, 9066909.48110164, 9066939.0625, 9067093.723861238, 9067145.3125, 9067179.6875, 9067247.905065252, 9067314.434433706, 9067383.916610664, 9067411.974206004, 9067421.875, 9067463.288606301, 9067533.20841679, 9067562.292744435, 9067566.478803854, 9067648.344025167, 9067799.012766793, 9067979.6875, 9068041.77729005, 9068117.158770138, 9068486.447622718, 9068523.950253056, 9068767.1875, 9068840.62600095, 9070351.527994512, 9070593.115576735, 9071427.604652686, 9071523.114482097, 9071887.5, 9071983.578473588, 9071989.0625, 9072018.812990367, 9072409.375, 9073620.3125, 9073904.609205527, 9079628.08675259, 9079909.080042688, 9080418.531281415, 9081234.571802793, 9081265.625, 9081870.721764483, 9085774.34244485, 9086248.4375, 9096315.625, 9109050.98628012, 9109132.8125, 9109457.781630544, 9109699.717481546, 9111366.180808999, 9113138.88884435, 9115875.0, 9131820.13874599, 9147349.445457477, 9159259.375, 9160656.140578115, 9165311.122044193, 9167652.231229564, 9171349.105878329, 9172940.356509535, 9173229.024531271, 9173242.222816251, 9178196.364252988, 9178356.25, 9178383.955048038, 9178865.625, 9179410.422848966, 9179777.13656335, 9179841.026321342, 9179868.113863625, 9180475.929992128, 9182109.375, 9182704.849608814, 9183585.483051082, 9184177.752318073, 9184590.385287365, 9185301.033749372, 9185478.363726486, 9195449.057354204, 9200950.00733304, 9200998.4375, 9201441.400513928, 9202943.2137763, 9203593.352078538, 9207244.239531096, 9207614.192748923, 9207905.513074651, 9215696.660701513, 9217120.3125, 9217249.351318192, 9221911.965753097, 9223551.845500352, 9224036.63207275, 9229127.409765972, 9230625.0, 9230674.197268892, 9232825.15064998, 9236192.24933635, 9236565.625, 9237231.951013215, 9237239.461285375, 9237421.875, 9238530.063670369, 9238602.326439526, 9238670.843280602, 9240362.405934347, 9240542.58173975, 9240676.542443758, 9241408.518504225, 9241526.5625, 9242412.5, 9242884.375, 9243115.425720569, 9243726.98515286, 9244298.4375, 9248301.5625, 9248686.374043081, 9249471.22349186, 9250530.433762833, 9250697.830740603, 9250924.338311244, 9251329.68326511, 9251943.54174985, 9252479.6875, 9253012.5, 9253053.023174386, 9253705.585383918, 9256217.219854118, 9257523.4375, 9258080.845368464, 9258095.3125, 9258101.199564034, 9259207.19113179, 9262720.089933092, 9263687.5, 9264140.551223606, 9264433.998159565, 9264440.453284826, 9265517.404212292, 9269219.783151336, 9271054.32465812, 9271418.897103176, 9273351.000810036, 9273642.1875, 9275489.53549253, 9283621.840782374, 9286755.262242535, 9290206.17938272, 9290893.366560556, 9291002.993811766, 9292003.091235429, 9292214.491689194, 9293006.25, 9293141.8697256, 9293164.915213862, 9293335.9375, 9293754.885646014, 9294297.924650637, 9294316.579789124, 9294490.625, 9294793.535139622, 9294964.0625, 9295060.699740462, 9295157.783238593, 9295595.3125, 9295821.593704496, 9295830.711258667, 9296057.8125, 9296547.728200698, 9296898.4375, 9296917.903748617, 9297036.101665298, 9297085.906624135, 9297206.25, 9297495.529910361, 9297789.0625, 9298067.155095145, 9298157.45720855, 9299418.985953556, 9300103.014735801, 9307010.15055015, 9307911.454167297, 9308871.875, 9314456.6044502, 9315615.625, 9315801.649536349, 9316337.91821196, 9316378.125, 9316900.563317673, 9317332.84588332, 9317701.900590431, 9317969.43420567, 9318466.029479621, 9318790.625, 9320993.62328599, 9330527.674047578, 9331041.438932614, 9331164.0625, 9331577.921955507, 9331852.5725616, 9332437.5, 9332986.884632403, 9334005.506223986, 9334146.875, 9334386.362836745, 9334438.206567025, 9334624.66590418, 9334994.884104067, 9335202.262524545, 9335471.218261806, 9336109.685655365, 9336599.564251667, 9337186.842468543, 9337300.0, 9337470.89909258, 9337773.628468541, 9338570.3125, 9338996.445688896, 9339157.411625843, 9339460.739591464, 9339620.072334826, 9340141.672891175, 9340228.24405665, 9340231.919399926, 9345273.127970248, 9345720.744557677, 9347929.92298877, 9349310.9375, 9353172.466441734, 9354263.507174898, 9357708.789581047, 9357737.5, 9358298.4375, 9359907.737752508, 9361405.221549703, 9365042.1875, 9366331.25, 9369052.985822245, 9369233.443840804, 9375553.125, 9378942.134353943, 9380621.875, 9381745.188794693, 9382939.980990546, 9383186.52889745, 9383194.724825159, 9383218.75, 9383400.42650697, 9383411.039667604, 9383467.1875, 9383518.053417074, 9383518.842151439, 9383548.4375, 9383571.210749764, 9383574.991915066, 9383624.505602453, 9383626.450180985, 9383649.002135746, 9383649.89990699, 9383696.17235585, 9383705.526467063, 9383818.234928422, 9383901.681584785, 9383919.80316099, 9383956.860374065, 9383959.09677875, 9383990.026998747, 9384007.108439788, 9384019.46886308, 9384068.651329206, 9384074.753014965, 9384084.375, 9384087.5, 9384098.736740664, 9384112.5, 9384129.736935025, 9384133.363028757, 9384133.690564953, 9384135.9375, 9384157.68631216, 9384213.243030472, 9384218.05240922, 9384225.0, 9384265.154551305, 9384312.5, 9384326.505917348, 9384341.50534859, 9384356.247788457, 9384357.8125, 9384360.349823713, 9384368.558008399, 9384411.21654247, 9384426.5625, 9384466.588870592, 9384488.239121938, 9384507.438372737, 9384523.37776976, 9384539.782972625, 9384552.334091563, 9384603.127578856, 9384617.1875, 9384664.501116205, 9384697.992875282, 9384710.9375, 9384728.293236211, 9384767.1875, 9384775.346470393, 9384785.584538156, 9384794.777358059, 9384819.076719295, 9384823.4375, 9384837.571506204, 9384851.937856372, 9384873.915697863, 9384874.869817106, 9384892.1875, 9384895.3125, 9384895.891132612, 9384897.16701084, 9384946.875, 9385003.913663138, 9385007.8125, 9385008.5348641, 9385009.15462192, 9385018.396866033, 9385029.6875, 9385058.33446472, 9385102.959184108, 9385116.42490072, 9385121.875, 9385131.42161185, 9385144.988563828, 9385146.231073383, 9385159.57488966, 9385181.872342804, 9385187.339553459, 9385200.713205244, 9385259.851077123, 9385262.399786592, 9385337.326968024, 9385345.232064677, 9385353.081587017, 9385374.495866934, 9385380.08281271, 9385404.6875, 9385409.735008465, 9385423.876734253, 9385426.5625, 9385430.572978152, 9385431.728160482, 9385441.63424752, 9385496.875, 9385505.130497972, 9385507.426505128, 9385535.9375, 9385567.710546523, 9385568.75, 9385590.625, 9385625.0, 9385634.375, 9385638.90198532, 9385639.308358315, 9385647.869562797, 9385699.483459262, 9385748.141833985, 9385750.880192107, 9385760.646970747, 9385771.875, 9385783.480476495, 9385796.875, 9385801.433802659, 9385802.884297783, 9385818.75, 9385826.5625, 9385839.265898142, 9385868.691069163, 9385875.284486553, 9385877.48968759, 9385906.578791907, 9385910.126563957, 9385914.0625, 9385929.324142477, 9385973.288077647, 9385977.899757428, 9386001.24668379, 9386042.881670369, 9386054.649106922, 9386090.625, 9386090.625, 9386120.120026069, 9386151.379667, 9386178.125, 9386226.055003846, 9386246.813125495, 9386312.1611892, 9386314.015267078, 9386343.156216921, 9386375.0, 9386403.711666008, 9386417.1875, 9386422.647698505, 9386453.218032125, 9386537.5, 9386673.4375, 9386690.625, 9386758.824125618, 9386826.5625, 9386971.706984648, 9387017.277290575, 9387033.722721951, 9387100.711479243, 9387189.20231424, 9387253.04151081, 9387298.458652133, 9389617.94632632, 9392402.283774145, 9394149.975087762, 9401037.95764324, 9410630.017878864, 9411945.670703603, 9412532.328363247, 9413253.571260372, 9414887.256361376, 9420804.6875, 9421105.89619378, 9421107.69383939, 9421124.681637691, 9433381.917005433, 9434098.58426839, 9434345.726710493, 9434370.3125, 9435955.104388442, 9436375.83997479, 9436376.497665443, 9442838.822168848, 9446480.700714383, 9446504.6875, 9450931.25, 9464278.976143913, 9465028.441702409, 9465869.01443272, 9465938.003995422, 9465971.875, 9466097.81185165, 9469341.595293041, 9469390.47860369, 9469809.704222381, 9469833.488171494, 9469838.467163462, 9470042.174673326, 9470124.51300583, 9470309.375, 9472103.674934672, 9481967.1875, 9484439.36304169, 9484616.671628812, 9485519.390108004, 9486049.022658782, 9487785.9375, 9489252.328555718, 9507384.952241799, 9513228.125, 9523429.60364719, 9524304.6875, 9524514.0625, 9526527.689923856, 9526564.507349363, 9526589.376933398, 9527387.5, 9527800.37000523, 9532643.75, 9533010.844629444, 9533819.592208184, 9534866.235476708, 9535036.700544922, 9535350.0, 9545030.321029177, 9545100.99637735, 9545227.134004977, 9546106.384401431, 9546498.165281378, 9547285.666247137, 9547300.0, 9548478.125, 9548645.310451858, 9548740.899207123, 9550031.34111165, 9551707.8125, 9552605.401446996, 9552605.851353938, 9553657.484965513, 9562054.56344031, 9562371.875, 9563629.6875, 9563693.67078937, 9565815.664124388, 9575085.257596007, 9577381.176103098, 9577963.906861138, 9584071.278416738, 9589389.82176476, 9590979.09258505, 9591040.386781741, 9592120.612660326, 9606948.56807796, 9606976.02125493, 9613446.85710067, 9621459.255553292, 9621836.785819778, 9626922.235126052, 9627074.65278923, 9628367.633027893, 9635205.239818064, 9635745.740496589, 9648885.9375, 9652495.3125, 9655807.8125, 9656141.845308196, 9656593.75, 9656609.08282424, 9656742.004677035, 9656922.422898443, 9657005.998099044, 9657039.0625, 9657192.051306926, 9657221.875, 9657260.929231947, 9657265.773876911, 9657266.149383234, 9657268.577441236, 9657370.081132371, 9657382.8125, 9657561.061781563, 9657573.616184356, 9657628.125, 9657641.35431838, 9657642.971053936, 9657643.031103864, 9657648.336379016, 9657654.6875, 9657671.165542131, 9657715.544216901, 9657799.472992625, 9657801.424135646, 9657859.375, 9657885.9375, 9657999.132716931, 9657999.7062979, 9658016.846393414, 9658033.42021039, 9658134.375, 9658158.699094893, 9658163.337940188, 9658173.2949372, 9658176.853542052, 9658183.25571347, 9658187.5, 9658197.272923008, 9658214.800368657, 9658223.992459789, 9658239.0625, 9658246.359825965, 9658255.63293468, 9658275.0, 9658278.144343846, 9658286.9834116, 9658301.5625, 9658321.875, 9658329.473593608, 9658332.6413758, 9658387.258843498, 9658395.265774716, 9658413.95845659, 9658416.071465353, 9658416.668557458, 9658418.657321084, 9658424.619076567, 9658464.0625, 9658468.75, 9658493.402294464, 9658522.828267436, 9658543.726739371, 9658546.93213835, 9658553.291084366, 9658559.375, 9658562.3143369, 9658568.29072699, 9658587.195126064, 9658619.73635291, 9658629.529785713, 9658660.389166312, 9658671.404018918, 9658689.0625, 9658690.625, 9658691.846906701, 9658700.96590776, 9658707.144443003, 9658708.848314336, 9658710.999314388, 9658747.402172104, 9658760.641002866, 9658787.5, 9658792.1875, 9658807.353921818, 9658810.98537383, 9658829.093074469, 9658843.644490767, 9658846.66323796, 9658926.5625, 9658952.560560957, 9659020.3125, 9659024.045865737, 9659039.0625, 9659063.717655929, 9659067.086802756, 9659178.612032536, 9659201.552739644, 9659215.22642859, 9659221.115147658, 9659290.341234531, 9659306.968595747, 9659320.759863622, 9659370.3125, 9659521.875, 9659546.967779955, 9659574.92076158, 9659595.3125, 9659601.781027209, 9659621.09103631, 9659626.146711342, 9659645.89125797, 9659650.742401293, 9659704.284199376, 9659743.267198969, 9659760.9375, 9659764.963497555, 9659801.15869171, 9659811.529305374, 9659834.550973622, 9659930.808898399, 9659936.577534705, 9659966.015420724, 9659970.3125, 9660145.255293535, 9660218.759104874, 9660311.362844318, 9660378.516949484, 9660434.375, 9660536.748454707, 9660584.847487042, 9660593.346921958, 9660708.422831573, 9661235.9375, 9661400.0, 9661465.625, 9662012.452187384, 9662189.992382219, 10426407.582043407, 10467749.378067244, 10540810.9375, 10543097.839392321, 10671911.547558827, 10717692.571378345, 10745472.12049515, 10748170.3125, 10748569.88804872, 10749989.264167875, 10750039.728189787, 10751219.836592996, 10751967.1875, 10753414.0625, 10753906.187087854, 10754763.858625816, 10754979.143617032, 10754991.08121027, 10755305.881249318, 10755352.31359794, 10755411.486346485, 10755852.02375274, 10756001.013996137, 10757195.068804407, 10758076.861832762, 10758547.463982968, 10758949.2931958, 10759239.392921679, 10759328.125, 10761004.213081356, 10761106.25, 10761150.930300463, 10761359.176743908, 10761412.5, 10761633.29782466, 10761872.19917027, 10762755.550410097, 10762799.915923996, 10763397.472332638, 10763780.936076377, 10763886.167680684, 10765110.823693387, 10767200.504363976, 10767678.791481636, 10768155.453256307, 10769157.785072489, 10769243.75, 10770129.6875, 10770204.13605203, 10770440.90994106, 10771073.179616235, 10771683.454448944, 10774063.398531575, 10774303.203784268, 10774348.979502395, 10774814.11509693, 10776589.0625, 10777628.109538162, 10791905.43411666, 10792365.118187916, 10802035.311484888, 10812466.119662225, ...], [69.77826560890375, 91.55267732843626, 9.48260086091381, 58.110938815141836, 12.18531375387578, 10.488654754911623, 13.757151043791694, 29.42451301356898, 18.14652392422138, 17.013698643403096, 50.135531786664764, 25.078090569469914, 66.15028514486116, 43.26247340860597, 45.22565552233857, 61.12638228105293, 45.25321925344876, 7.321984280432101, 42.97378178584038, 98.00823560613298, 86.50762485022932, 12.713119248037552, 25.07212283452826, 18.66793346000179, 21.849202220640354, 70.85218833703803, 59.9940620473035, 49.7914583735647, 5.467094748059871, 5.771918492542404, 117.35535142021838, 54.07610709568627, 78.81090848751654, 53.76496256509341, 37.820336784745145, 7.497595782644844, 50.19741674938565, 17.637126278536385, 30.183405516532815, 56.94435682094335, 87.81304566695417, 108.2196887180205, 80.03361576365157, 103.455102875924, 91.15444897749427, 25.252121254335393, 42.45770497107402, 53.12210812504012, 5.602499948103115, 23.551068338210577, 36.17739737275985, 49.25835630793095, 63.514911442422274, 54.27549545299744, 9.665480726269346, 67.51404859695468, 23.846389566681204, 69.93098190327427, 8.359168216362772, 17.725468508667205, 36.07742534822684, 41.05005251652426, 148.12882075413344, 21.306393983357516, 27.841073594507304, 70.31016320435519, 18.890571742234965, 5.760026617818309, 5.21229544370629, 53.237428798964054, 39.721450213582074, 96.65502152448519, 111.77604791726147, 31.021309875172797, 36.52866313107864, 56.55636521487171, 59.162723263305374, 5.545706046516259, 9.95423923825268, 51.959870332086965, 61.37472980352535, 75.55777427819783, 5.671908389441759, 7.098069513803744, 11.773988836274457, 15.288505632155287, 59.95914584710735, 11.052700217418792, 31.69626150940286, 46.58064869086881, 7.509058905494638, 35.28736434585281, 51.28882097167621, 18.884092905126163, 104.26727288438829, 59.37053935191965, 15.154396030690538, 64.14409251367024, 14.255488021367157, 102.19269292896851, 13.491385710057397, 82.96419980923928, 21.101392055129082, 74.09458481242889, 14.393799648652077, 6.4400771223857545, 24.240168078439353, 8.55926391178432, 19.96530706759422, 80.76747958587946, 28.080431140001963, 84.60793485446379, 30.188562014509337, 11.888317918385392, 74.8438662612712, 13.281001884818322, 88.96939309636943, 9.558897321027768, 22.023623449562486, 25.56263039213629, 5.719657638938461, 54.395002641971715, 16.11847539616461, 35.95064151530201, 56.586228774140984, 6.339117604009785, 23.179439585480512, 21.396556550613766, 12.642559883862226, 8.482338739056356, 22.101307208431166, 10.423322274980073, 96.30515458100007, 102.65855846653412, 16.416944106737866, 69.92368849194668, 21.000230826399413, 27.355706835256676, 29.72659693089, 46.331622054894396, 5.083385012586598, 50.34345415570628, 39.901953554923, 8.690298301484715, 114.58837202716028, 69.25984415193179, 24.483093639604494, 13.233463316965059, 35.07533709740332, 38.601597606352726, 7.195471338449034, 12.638646440861073, 26.699677797946727, 7.190233169668652, 44.342754045220715, 15.296342450602864, 60.96775148683466, 98.14002113055867, 5.660679524503315, 5.908668401920862, 11.320122081198331, 5.466998104572451, 151.76435792944645, 29.650971082229056, 20.56394780745952, 17.926491000605367, 11.013858865676111, 87.82463792104352, 61.88881908172403, 10.852508559356604, 108.35878614638794, 27.713020986304016, 47.74661144678452, 113.33298046853396, 48.62698687851665, 25.084222448780782, 8.127945293613735, 21.572574124053528, 11.446825413902843, 11.835810773849381, 20.518112986889427, 9.28761822442976, 182.44865673331793, 72.09392643970429, 6.757078370448362, 37.53434549947089, 43.862184228860706, 12.11818223152035, 83.99641048909596, 6.488598538170047, 24.50078737351882, 22.414197838065576, 25.69164853956341, 5.965960585589358, 33.79061639437341, 39.394399632736864, 15.188309675610437, 58.27062158239014, 13.719541369234069, 12.491839798950364, 42.30669312942717, 33.1458119567504, 28.527440687348765, 80.5503299019194, 10.510184917078977, 26.112351842571165, 17.912004140621757, 38.85215284886852, 25.87608937406097, 11.010932863726408, 71.51963132412364, 57.70011033443569, 40.5382917285101, 77.36689761857832, 48.06874733154836, 73.7649713873403, 6.021590453166321, 35.28392573072457, 34.406308164104175, 28.84738541322323, 5.724146255085033, 5.596267675305432, 9.233385842603509, 9.818104164251224, 31.90194260980565, 32.38416776090794, 23.785957854336807, 99.63791219796497, 8.37336435543623, 48.714538973921876, 95.50427948263514, 27.504528514206136, 21.808520984198363, 17.279396450051234, 61.16769341308694, 83.41373918909656, 113.53813997056459, 23.859657102399453, 54.29486651058716, 39.92188747251334, 25.54224565423921, 23.567788374344858, 57.76254494149095, 14.124588759434356, 52.83092875570432, 51.82801483736395, 10.259721238681559, 74.780426753657, 25.54930961553815, 16.247641292194967, 9.718041601855553, 62.315855927097004, 27.46356148824193, 14.575677104897352, 20.255432540843824, 69.65170981410218, 69.30552318905077, 110.70010281924421, 69.9128692736494, 92.51425395888629, 126.56316952612504, 146.5513798027724, 22.120894035622314, 6.927348354020946, 56.689154665369486, 5.196868511760627, 23.239373645660116, 92.00801747574113, 9.472121651408102, 29.708178673980022, 20.967126254218723, 35.00287511331956, 78.58197679064716, 27.95378145571374, 6.5503867777059686, 11.082052860565552, 34.48474351626413, 46.40021736496004, 37.05272406409659, 41.40087410083035, 23.102873416494667, 18.05653922382973, 8.408888953239286, 6.784875430665667, 32.67927096230147, 7.114597481842072, 78.85094021770959, 65.19907024630844, 21.175535193090656, 113.62528944644407, 7.445453559130028, 13.255953163892276, 92.4478821070038, 8.752828319539226, 62.47474399185407, 14.206312089854476, 19.636527769139082, 57.42691269989697, 58.033160042405996, 89.96255001832894, 52.39957507407625, 39.612089620570245, 28.292023645302557, 6.870381315887813, 92.8890883839816, 22.844380161163226, 37.685547258308354, 22.696896710220308, 17.280384425796413, 5.51760618041778, 65.01247399981318, 18.970086697214956, 29.614079755906488, 30.88702771837743, 30.19455806839862, 38.9873339125901, 41.43169696885861, 34.56333666122596, 8.136524010449943, 24.84766113885159, 52.50079941140823, 42.19547848925977, 92.33245360661022, 20.053301029792472, 107.04747762077136, 93.6232831082136, 50.64458534228507, 53.95880776982894, 14.825140580137498, 19.205732972225306, 6.879785543755798, 9.188760695686844, 40.16983503560883, 33.22950870243038, 28.475210282620367, 67.63021431687542, 109.71960166266305, 14.786023503875086, 45.81227220384923, 41.93690318060415, 96.82889999071415, 51.82011710029744, 15.077848967137742, 16.42255082775159, 30.579413936515525, 42.88533437386914, 71.50555989274369, 8.931110209888228, 17.698367734420394, 16.72897683872372, 7.9765774331894574, 31.0869308261036, 47.57458007513725, 36.20160207781897, 5.381380488049234, 59.32357859321691, 6.450159446990361, 13.247988546154946, 7.185652170532421, 5.861313367581496, 6.764971436503363, 74.22611078713354, 23.825951508019, 6.987285139046644, 74.02427294602852, 95.7766083669544, 18.2134675453204, 27.469551590912438, 10.808488974648137, 112.7176341480447, 7.384341120742907, 35.651358806444996, 7.197885467963935, 45.074030655385855, 40.399279351448925, 70.86625677576966, 6.008419338358145, 20.698578336306483, 22.088825287063038, 54.46707939245381, 39.440459336501974, 5.7326166801739395, 33.67928179453982, 7.7452045615218115, 86.0593088256797, 6.874310394906235, 26.793549312090107, 25.462368443591785, 30.069386269026083, 12.44534051936803, 33.328928890411404, 29.46860812153487, 45.51092059449774, 25.128617533951314, 41.471124129808146, 28.708237989294137, 41.33286033423291, 16.9150463847392, 20.104087066646258, 6.366323736958948, 28.620024205295294, 8.828264785054541, 8.284339114761906, 14.203029672649995, 42.657924876798106, 55.486539704723974, 6.25441858712848, 7.963832727428369, 30.67098534775709, 42.615415609572125, 25.780350975092084, 32.76428229249645, 22.97878734531394, 29.304694439916226, 9.15004322425323, 64.79143367109447, 14.667571965803345, 13.676500641371305, 25.540463303240664, 16.85400792516345, 24.08524227352808, 22.97989556067756, 8.36593063420316, 53.173289804208565, 65.77447399368772, 160.4507280206815, 85.25189631004315, 14.660170047425588, 11.019092754691924, 36.64234042926972, 52.166359098954416, 131.44005705698146, 71.95922374674139, 109.27659166187777, 8.12538994839134, 24.207479589659716, 54.49864234514353, 8.394215456708961, 13.29099188085044, 87.25574374227118, 57.14132893996474, 29.99940454200873, 49.87084521061469, 16.446007253148593, 14.831139371859233, 76.65810104222564, 13.967785142642757, 21.115148364671953, 48.64735902389387, 6.62970676344421, 71.31736722240434, 12.468233825043775, 6.513137231132816, 61.69115976618602, 54.69081165622625, 32.963092682298196, 30.384666757174628, 27.89930970396219, 23.069261429950387, 5.9602205989615396, 34.83156137260096, 7.910088134013308, 17.359875859538437, 62.45234644551324, 16.48973797589025, 48.70373538989847, 11.266531131759114, 80.75471036608904, 71.77061776585, 9.857928625729334, 66.08535707538118, 49.848278498878656, 6.357694400854463, 87.06081404735941, 70.6627637770713, 32.03186611556872, 10.489117919340877, 153.46168328427834, 80.82039310971457, 14.22023405971005, 13.184687899079929, 73.29301168875114, 31.14870292714661, 19.363197210734917, 63.565408336594935, 11.864680640969805, 18.551292596396358, 7.036312773519994, 22.238418985422385, 18.506631252375602, 12.91634100634451, 70.94886765954244, 94.91807868698646, 7.381557464355967, 110.69531224531683, 106.34017512554217, 15.723031820776056, 29.409838891722575, 82.77878549863067, 22.54241075026417, 79.56421189125302, 122.9980516510281, 59.604921514839035, 5.714377217766062, 5.081274727558386, 23.629770790920563, 16.891040997691864, 11.51688395950215, 33.68587371120256, 17.8043541830509, 45.32292703450922, 20.681318307821012, 51.75598295007495, 63.377263088492185, 14.664060544467723, 17.436911749376428, 48.8104286884791, 31.467838187166503, 26.048113668271522, 26.273024039091503, 89.67299082683479, 27.515402816751955, 88.3901826370805, 26.29814146026113, 45.49511398434288, 44.90524951118547, 26.532478692975445, 55.2428955488178, 49.793403544275755, 50.87569320749922, 62.03314406369785, 78.97821755623272, 15.304778900270604, 62.201247421478456, 19.99745237598238, 38.52839964422239, 21.54917034216446, 103.51565854414545, 17.64156608598151, 145.4593856524358, 26.657779712734758, 19.519720315460518, 127.06404601446914, 75.36540517033468, 8.934913916819317, 9.90697976111298, 6.7092045550091095, 8.87224021382223, 38.18152125325919, 12.474097126212813, 10.976902718818442, 23.49248926413497, 44.047635883158094, 81.85640687175498, 5.726070875819977, 28.746029952540848, 30.876854950818323, 17.55746724957502, 62.144187615252555, 44.528774948730245, 14.481407248395923, 6.168933455879589, 12.815617864883375, 59.153207612370885, 13.691497122321188, 101.9284440798304, 22.347845293845147, 55.238201054080335, 28.481461562771454, 70.35314407269958, 13.092168474508384, 51.36845124333529, 41.721283775918565, 93.46683385249949, 20.800513050628126, 18.777776585600243, 5.90747328807188, 27.396136088327772, 18.647158631583117, 74.14988856292834, 14.95134793238071, 46.412787311015656, 8.532993988326437, 8.064661307903851, 29.40156838792065, 5.623419667498884, 59.53769571011693, 15.85178728808059, 23.097007536979238, 13.230820041497793, 11.67991406003155, 32.793827500641655, 23.680092722443156, 28.487469785378735, 34.22660152391056, 60.51935085050796, 64.2696727230055, 76.57524488209904, 9.873797679732139, 9.591154737291168, 95.49066297627802, 12.633821111497921, 46.6190565844058, 18.660298562787762, 31.123500360639106, 69.3944357263281, 45.91405222458597, 20.39776673615765, 20.38151942815262, 10.113691739717485, 84.38185358806999, 15.638609193024253, 8.122394366764263, 20.168226361937545, 27.822636818986375, 76.25945900789, 25.713568021890932, 18.59689613864327, 17.130666596408332, 26.688964911396912, 77.54017381426794, 9.264304194180871, 73.02081844012453, 59.18853490417307, 82.09590632966399, 48.51829706190693, 13.80238125734703, 43.93863612471406, 33.20971224524206, 19.89200826635866, 7.785847241239852, 48.60158334400493, 40.13477631564674, 6.855292216660307, 6.428238907641042, 83.47370864726389, 21.87070321216279, 110.16524089744877, 100.9008443695833, 33.15516884757169, 54.6271857424183, 14.282822637060761, 27.06758609068604, 5.583378603177801, 21.329669333758147, 23.80626062898029, 19.848780269837988, 29.171047211904916, 117.41566427668782, 80.26145077099804, 56.06546696714246, 8.219672710772471, 6.7163851879354945, 44.76630773488396, 30.87254847041948, 111.62958020851474, 60.05537100024193, 16.149333841698024, 60.05275947766313, 17.78581064773958, 6.367064299652842, 77.2479352983438, 16.642855823419023, 14.401007980535407, 23.182032338039065, 23.550792186263195, 26.186780864506463, 103.17134443433095, 34.59514545535535, 40.7988887733481, 133.09142995178235, 54.6440877361341, 46.96550380001123, 6.410374507676841, 5.24341684598228, 13.790303806704257, 16.756583799575893, 38.84844856884651, 54.653141264019645, 103.03755278251344, 29.421713828057616, 19.01905487800498, 11.276303739527185, 63.65566606747419, 35.09023204787191, 41.8347777387739, 25.522313355936795, 49.3888453456664, 7.543579133817603, 5.465633693719902, 5.271355276151485, 18.585010653915838, 41.987020887939195, 33.39544699610464, 5.098872764668381, 6.613008175899193, 33.50148939358008, 11.088989365801158, 11.927247712382789, 11.686608334973291, 25.161520452523277, 21.65281772680007, 7.773031395135736, 24.498405631440903, 58.52275196722579, 19.98488787771426, 6.468649639701314, 13.852649673103077, 39.530055844834955, 11.365135482619827, 133.76187876718976, 84.7384813813763, 6.554490104850635, 7.938191585915324, 5.809198587663047, 104.97410782690235, 13.19306113167598, 67.41521014959426, 132.89172246295493, 37.32232045216823, 106.96845039403092, 166.7920621873063, 7.796296673563342, 62.25816501533568, 10.272871640949646, 94.88455433207378, 10.906338731637183, 26.51085657130392, 12.378305105214896, 61.133276583946135, 73.06133011517392, 21.623828901995033, 130.73248346426675, 12.111854335398787, 64.58310030372296, 7.202761413681687, 6.66470913698336, 7.300536287094914, 13.66260938967418, 31.36215478048632, 69.08257636658377, 29.91356266487371, 147.20001365691377, 43.287191073634915, 35.98030459908736, 51.818181167478414, 23.306662844659694, 16.959743958596363, 20.320186323756786, 82.49483306967863, 23.87802323874675, 33.21173986036255, 84.65118678111139, 10.161947438066907, 5.9943931676093705, 18.94886324706177, 63.275757943643846, 113.37240814371873, 86.82889063994705, 17.682905557180604, 5.715365398316607, 51.9313193529557, 24.125025749365097, 67.10846952277934, 86.33003899555197, 11.124506932028083, 7.705518793252499, 11.273170528706327, 40.778511012024815, 92.25678238781433, 22.567473511396813, 52.80326704158034, 14.488835629977112, 66.7529617386063, 49.069312322289875, 56.8709807779707, 5.246241624991588, 19.868408961669093, 175.26049278450168, 20.236829332667504, 20.580859377899355, 18.921560029395177, 56.64419763233434, 98.26725037637073, 11.631693985190036, 11.338132534469736, 49.851144234986236, 21.342413660703365, 53.72032388418809, 19.593143817440968, 11.741479634264495, 106.513350608544, 5.305050581674112, 51.21543539192766, 30.365448796487218, 62.6179015062853, 82.33111529571401, 38.38677243626686, 19.280299440029765, 62.04096657765081, 10.699792804440605, 58.03534618647714, 25.20306004959948, 79.94184810960904, 124.90520617156518, 9.218839372630129, 47.05321074380154, 49.68121313762114, 46.17582538454992, 17.404379228090576, 25.52761709067381, 6.718178458050407, 55.74669613189269, 7.154509276863198, 11.418980654067878, 72.39913812446233, 25.71048634471073, 12.621917317343861, 5.225460690407179, 29.16689130934563, 52.82118755875637, 64.37664063591447, 60.43871007587107, 13.252682827118646, 83.25972817088461, 104.5110256013063, 84.77956184298004, 18.47399121285312, 23.831523349094752, 10.633346382443696, 17.90510779735961, 50.32705636572011, 6.386111902768931, 43.140198657118546, 22.913514130900406, 5.474398605351051, 25.20467608790601, 37.21375192816757, 37.60414558583007, 8.719542550059671, 150.90357192146126, 31.65445329142686, 8.031721305098555, 13.016472551202959, 37.49072403505779, 81.54250882004574, 64.140469946995, 39.446260516502456, 79.92310192098643, 14.184403899686975, 125.00160504645694, 60.27709883940308, 14.766355922354824, 7.566035487226929, 11.90548248748581, 6.20378697940843, 26.817275631840847, 22.769311485704314, 51.53068194112272, 32.716439998202574, 13.763438673215541, 6.970011162992154, 28.072632619209237, 12.682655047810279, 15.946165161397058, 66.28638209950796, 25.33004430235787, 58.43077602174903, 23.68868289961794, 10.332587592674747, 5.336546280758336, 31.041908306121343, 86.92095281556595, 57.38300761848095, 30.09479958875559, 133.33951786138735, 35.010125498291586, 65.72176873422268, 52.59428477568801, 15.1824701921035, 26.04619362604328, 157.36656400576214, 87.87778969766059, 74.81181180318475, 28.94930427900464, 5.535408529767064, 13.074668224043917, 6.2539594389111315, 9.1585429445389, 66.5214701444361, 25.4385661763878, 61.47616848217645, 12.2334460682809, 37.95016324449648, 70.91680563317215, 26.680775348307794, 30.700660686446525, 28.333233915575754, 22.568314438068622, 5.497304982047925, 22.044268297984786, 60.21149368833324, 6.325423978666722, 158.06260499281035, 56.245033666554306, 5.465762899336625, 27.049859792403083, 54.16250030105062, 24.866674301939153, 7.3423362838986455, 9.002948551042019, 17.72683426307021, 10.049375478318348, 8.084229999503211, 56.06073486591909, 56.174026562532624, 42.193191136692036, 6.685481002078481, 17.806589890005085, 19.37101629501525, 22.101477452675635, 18.109836629526605, 94.71065431426054, 41.26115423990972, 231.2090986539324, 7.606534083473084, 9.718960799781552, 9.686370344131078, 53.22109012810866, 14.871360304667427, 36.4910188154177, 40.67924926450571, 14.843721351013137, 32.70807733900208, 90.70543992461603, 77.86839180027894, 14.762348830490495, 94.54933088113195, 9.143855656068174, 30.957854703081367, 45.91616662957094, 53.07797685442267, 5.422319615225528, 16.638852230915695, 5.7225621590299305, 81.55801831222949, 69.34543000426277, 116.174233207947, 5.052943386685845, 6.20888682727358, 57.07011504416968, 61.30565866995507, 27.338937424105687, 7.531127503764302, 58.26624235704098, 6.56106597321113, 22.71922041011763, 18.862919954016824, 5.217634064397359, 17.04895311596002, 13.19513978403063, 61.39846085789381, 5.726076903558251, 65.71865766160028, 12.793227880158955, 16.37063839248726, 58.10246394639079, 10.459158889853171, 41.577994702543215, 61.210382709045156, 81.16742583011197, 59.327075090281234, 11.98169076938157, 45.577654238132844, 53.06708763869012, 45.96337653185522, 48.113244558365565, 23.574501193124455, 28.12615448274729, 10.954360261608782, 65.07561663177394, 7.863496720619498, 12.580716471574828, 36.3830052370592, 64.13341800277311, 29.623305339205345, 14.192461473662146, 8.443234665134106, 8.15475500420226, 5.423595598713941, 7.964121563174872, 18.01102258969498, 23.920927179573255, 15.145847000260014, 50.91316388665784, 11.10362465561797, 10.054041913523104, 59.318671319467505, 61.497372171578625, 6.866450655096844, ...])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([4183873.4375, 4200589.0625, 4245645.486060796, 4245648.321076619, 4273347.460452169, 4314272.98618653, 4463036.308541708, 4516373.4375, 4519480.333225946, 4520039.226822126, 4520225.0, 4520815.422263485, 4525556.211428195, 4525936.918065219, 4526217.1875, 4526593.75, 4529406.065605748, 4529410.086621125, 4580491.832632525, 4586693.342413341, 4595065.625, 4596572.697327951, 4604054.484354057, 4711086.660961539, 4722958.058993217, 4810845.3125, 4814384.375, 4891162.5, 6695396.9985546, 6866860.352375621, 6868261.685191906, 6994003.125, 7099131.25, 7129463.278797314, 7130623.073756609, 7130945.010654807, 8185245.3125, 8525754.387882605, 8628820.776659146, 8696595.68875797, 8701679.491570348, 8703296.875, 8703743.628683295, 8705085.336426139, 8705454.924390363, 8706293.053784892, 8706748.4375, 8730331.25, 8731529.476926165, 8748865.30920586, 8749673.4375, 8750856.25, 8779617.1875, 8790217.1875, 8825184.74532947, 8834315.625, 8853562.445120098, 8863321.875, 8864176.697293386, 8886742.76603319, 8892519.140832322, 8904692.1875, 8904913.658927286, 8905095.020341542, 8905349.964492438, 8905650.0, 8918130.442457212, 8919474.243177438, 8922011.282267557, 8928745.3125, 8931389.0625, 8937252.011676298, 8947414.0625, 8952057.307835227, 8965704.685666695, 8968389.0625, 8969540.953686716, 8970240.882425923, 8970963.154934406, 8971078.125, 8975503.980992258, 8975887.5, 8976480.68812859, 8976954.5977973, 8977042.565043097, 8980778.213154148, 8980976.5625, 8980985.456670586, 8981047.939277556, 8981551.095397804, 8981751.29438005, 8981798.4375, 8982078.030210437, 8982135.129054297, 8982614.469613882, 8982882.12367435, 8983096.929754166, 8983256.25, 8983643.557566142, 8983773.4375, 8987637.05927988, 8988124.254688518, 8988876.295442553, 8996076.24896725, 8997781.09847637, 8999818.443383655, 9000604.424873259, 9001122.357356636, 9001125.894374518, 9001335.9375, 9001542.16464767, 9004281.3693377, 9005364.941423027, 9005503.562604766, 9005923.390236454, 9006085.420970289, 9006089.0625, 9006169.007326934, 9006280.44515653, 9006329.517238103, 9006368.950798277, 9006621.719176397, 9006647.152944328, 9006790.625, 9006796.875, 9006884.923797358, 9007031.321908507, 9007034.372597642, 9007046.104904564, 9007155.040208349, 9007167.259587096, 9007182.30970794, 9007227.141970402, 9007312.964309137, 9007439.402320221, 9007560.507849539, 9007670.558947109, 9007743.023725435, 9007760.9375, 9007764.0625, 9007819.531306375, 9008027.550515702, 9008043.839343544, 9008055.380151317, 9008062.5, 9008104.6875, 9008143.505249066, 9008199.465843676, 9008345.143207641, 9008404.6875, 9008411.211793315, 9008517.924231231, 9008560.837918967, 9008656.505350163, 9008714.0625, 9008840.750927528, 9008893.75, 9009074.880919412, 9009133.625789277, 9009182.45264877, 9009201.682209646, 9009321.996323926, 9009331.541482333, 9009363.358914142, 9009468.169689404, 9009491.236132098, 9009671.800921798, 9009751.5625, 9009754.405385567, 9009768.69761471, 9009918.39377025, 9009929.048365531, 9009946.984452512, 9009953.125, 9010026.5625, 9010035.73518393, 9010084.06247413, 9010100.137702378, 9010101.757499326, 9010116.035248041, 9010130.758633042, 9010135.885259666, 9010141.002152685, 9010151.5625, 9010167.15521591, 9010210.9375, 9010212.5, 9010222.428020593, 9010326.869958924, 9010334.384608725, 9010350.248461494, 9010354.845121, 9010377.016465543, 9010404.987460053, 9010425.872245127, 9010455.531909695, 9010478.18827834, 9010510.9375, 9010589.786125397, 9010597.601751687, 9010634.375, 9010642.1875, 9010657.8125, 9010668.75, 9010671.440412492, 9010747.83882921, 9010783.42131511, 9010807.8125, 9010813.912468296, 9010868.673984608, 9010883.92879088, 9010905.421934988, 9010912.5, 9010929.6875, 9010981.25, 9011003.32062525, 9011030.585197488, 9011042.1875, 9011107.8125, 9011150.0, 9011244.019793529, 9011268.424141327, 9011268.631759254, 9011300.524493318, 9011330.166647157, 9011347.483890036, 9011403.785153715, 9011485.086929744, 9011520.473087369, 9011555.373735217, 9011576.543449743, 9011648.38375643, 9011740.095098844, 9011748.940079143, 9011795.3125, 9011848.4375, 9011875.0, 9011901.270979028, 9011943.324927993, 9011951.5625, 9012105.207295107, 9012110.623040995, 9012123.06182649, 9012215.392783487, 9012225.0, 9012557.8125, 9012846.238316586, 9012985.653734386, 9013038.183354102, 9013056.555102466, 9013261.361907978, 9013425.0, 9015011.177805066, 9017948.338090392, 9019175.928783504, 9023421.875, 9024487.5, 9030298.4375, 9061534.749423213, 9063212.96730573, 9063230.76319783, 9063521.57554675, 9063547.998211809, 9063607.018780848, 9063667.1875, 9063702.865812365, 9063825.550099153, 9064021.875, 9064248.59936124, 9064320.3125, 9064390.070606042, 9064467.140083848, 9064715.280235466, 9064726.532983446, 9065384.767233357, 9066142.158094348, 9066198.905409403, 9066262.5, 9066264.0625, 9066506.25, 9066533.510633184, 9066636.528919457, 9066854.911846513, 9066909.48110164, 9066939.0625, 9067093.723861238, 9067145.3125, 9067179.6875, 9067247.905065252, 9067314.434433706, 9067383.916610664, 9067411.974206004, 9067421.875, 9067463.288606301, 9067533.20841679, 9067562.292744435, 9067566.478803854, 9067648.344025167, 9067799.012766793, 9067979.6875, 9068041.77729005, 9068117.158770138, 9068486.447622718, 9068523.950253056, 9068767.1875, 9068840.62600095, 9070351.527994512, 9070593.115576735, 9071427.604652686, 9071523.114482097, 9071887.5, 9071983.578473588, 9071989.0625, 9072018.812990367, 9072409.375, 9073620.3125, 9073904.609205527, 9079628.08675259, 9079909.080042688, 9080418.531281415, 9081234.571802793, 9081265.625, 9081870.721764483, 9085774.34244485, 9086248.4375, 9096315.625, 9109050.98628012, 9109132.8125, 9109457.781630544, 9109699.717481546, 9111366.180808999, 9113138.88884435, 9115875.0, 9131820.13874599, 9147349.445457477, 9159259.375, 9160656.140578115, 9165311.122044193, 9167652.231229564, 9171349.105878329, 9172940.356509535, 9173229.024531271, 9173242.222816251, 9178196.364252988, 9178356.25, 9178383.955048038, 9178865.625, 9179410.422848966, 9179777.13656335, 9179841.026321342, 9179868.113863625, 9180475.929992128, 9182109.375, 9182704.849608814, 9183585.483051082, 9184177.752318073, 9184590.385287365, 9185301.033749372, 9185478.363726486, 9195449.057354204, 9200950.00733304, 9200998.4375, 9201441.400513928, 9202943.2137763, 9203593.352078538, 9207244.239531096, 9207614.192748923, 9207905.513074651, 9215696.660701513, 9217120.3125, 9217249.351318192, 9221911.965753097, 9223551.845500352, 9224036.63207275, 9229127.409765972, 9230625.0, 9230674.197268892, 9232825.15064998, 9236192.24933635, 9236565.625, 9237231.951013215, 9237239.461285375, 9237421.875, 9238530.063670369, 9238602.326439526, 9238670.843280602, 9240362.405934347, 9240542.58173975, 9240676.542443758, 9241408.518504225, 9241526.5625, 9242412.5, 9242884.375, 9243115.425720569, 9243726.98515286, 9244298.4375, 9248301.5625, 9248686.374043081, 9249471.22349186, 9250530.433762833, 9250697.830740603, 9250924.338311244, 9251329.68326511, 9251943.54174985, 9252479.6875, 9253012.5, 9253053.023174386, 9253705.585383918, 9256217.219854118, 9257523.4375, 9258080.845368464, 9258095.3125, 9258101.199564034, 9259207.19113179, 9262720.089933092, 9263687.5, 9264140.551223606, 9264433.998159565, 9264440.453284826, 9265517.404212292, 9269219.783151336, 9271054.32465812, 9271418.897103176, 9273351.000810036, 9273642.1875, 9275489.53549253, 9283621.840782374, 9286755.262242535, 9290206.17938272, 9290893.366560556, 9291002.993811766, 9292003.091235429, 9292214.491689194, 9293006.25, 9293141.8697256, 9293164.915213862, 9293335.9375, 9293754.885646014, 9294297.924650637, 9294316.579789124, 9294490.625, 9294793.535139622, 9294964.0625, 9295060.699740462, 9295157.783238593, 9295595.3125, 9295821.593704496, 9295830.711258667, 9296057.8125, 9296547.728200698, 9296898.4375, 9296917.903748617, 9297036.101665298, 9297085.906624135, 9297206.25, 9297495.529910361, 9297789.0625, 9298067.155095145, 9298157.45720855, 9299418.985953556, 9300103.014735801, 9307010.15055015, 9307911.454167297, 9308871.875, 9314456.6044502, 9315615.625, 9315801.649536349, 9316337.91821196, 9316378.125, 9316900.563317673, 9317332.84588332, 9317701.900590431, 9317969.43420567, 9318466.029479621, 9318790.625, 9320993.62328599, 9330527.674047578, 9331041.438932614, 9331164.0625, 9331577.921955507, 9331852.5725616, 9332437.5, 9332986.884632403, 9334005.506223986, 9334146.875, 9334386.362836745, 9334438.206567025, 9334624.66590418, 9334994.884104067, 9335202.262524545, 9335471.218261806, 9336109.685655365, 9336599.564251667, 9337186.842468543, 9337300.0, 9337470.89909258, 9337773.628468541, 9338570.3125, 9338996.445688896, 9339157.411625843, 9339460.739591464, 9339620.072334826, 9340141.672891175, 9340228.24405665, 9340231.919399926, 9345273.127970248, 9345720.744557677, 9347929.92298877, 9349310.9375, 9353172.466441734, 9354263.507174898, 9357708.789581047, 9357737.5, 9358298.4375, 9359907.737752508, 9361405.221549703, 9365042.1875, 9366331.25, 9369052.985822245, 9369233.443840804, 9375553.125, 9378942.134353943, 9380621.875, 9381745.188794693, 9382939.980990546, 9383186.52889745, 9383194.724825159, 9383218.75, 9383400.42650697, 9383411.039667604, 9383467.1875, 9383518.053417074, 9383518.842151439, 9383548.4375, 9383571.210749764, 9383574.991915066, 9383624.505602453, 9383626.450180985, 9383649.002135746, 9383649.89990699, 9383696.17235585, 9383705.526467063, 9383818.234928422, 9383901.681584785, 9383919.80316099, 9383956.860374065, 9383959.09677875, 9383990.026998747, 9384007.108439788, 9384019.46886308, 9384068.651329206, 9384074.753014965, 9384084.375, 9384087.5, 9384098.736740664, 9384112.5, 9384129.736935025, 9384133.363028757, 9384133.690564953, 9384135.9375, 9384157.68631216, 9384213.243030472, 9384218.05240922, 9384225.0, 9384265.154551305, 9384312.5, 9384326.505917348, 9384341.50534859, 9384356.247788457, 9384357.8125, 9384360.349823713, 9384368.558008399, 9384411.21654247, 9384426.5625, 9384466.588870592, 9384488.239121938, 9384507.438372737, 9384523.37776976, 9384539.782972625, 9384552.334091563, 9384603.127578856, 9384617.1875, 9384664.501116205, 9384697.992875282, 9384710.9375, 9384728.293236211, 9384767.1875, 9384775.346470393, 9384785.584538156, 9384794.777358059, 9384819.076719295, 9384823.4375, 9384837.571506204, 9384851.937856372, 9384873.915697863, 9384874.869817106, 9384892.1875, 9384895.3125, 9384895.891132612, 9384897.16701084, 9384946.875, 9385003.913663138, 9385007.8125, 9385008.5348641, 9385009.15462192, 9385018.396866033, 9385029.6875, 9385058.33446472, 9385102.959184108, 9385116.42490072, 9385121.875, 9385131.42161185, 9385144.988563828, 9385146.231073383, 9385159.57488966, 9385181.872342804, 9385187.339553459, 9385200.713205244, 9385259.851077123, 9385262.399786592, 9385337.326968024, 9385345.232064677, 9385353.081587017, 9385374.495866934, 9385380.08281271, 9385404.6875, 9385409.735008465, 9385423.876734253, 9385426.5625, 9385430.572978152, 9385431.728160482, 9385441.63424752, 9385496.875, 9385505.130497972, 9385507.426505128, 9385535.9375, 9385567.710546523, 9385568.75, 9385590.625, 9385625.0, 9385634.375, 9385638.90198532, 9385639.308358315, 9385647.869562797, 9385699.483459262, 9385748.141833985, 9385750.880192107, 9385760.646970747, 9385771.875, 9385783.480476495, 9385796.875, 9385801.433802659, 9385802.884297783, 9385818.75, 9385826.5625, 9385839.265898142, 9385868.691069163, 9385875.284486553, 9385877.48968759, 9385906.578791907, 9385910.126563957, 9385914.0625, 9385929.324142477, 9385973.288077647, 9385977.899757428, 9386001.24668379, 9386042.881670369, 9386054.649106922, 9386090.625, 9386090.625, 9386120.120026069, 9386151.379667, 9386178.125, 9386226.055003846, 9386246.813125495, 9386312.1611892, 9386314.015267078, 9386343.156216921, 9386375.0, 9386403.711666008, 9386417.1875, 9386422.647698505, 9386453.218032125, 9386537.5, 9386673.4375, 9386690.625, 9386758.824125618, 9386826.5625, 9386971.706984648, 9387017.277290575, 9387033.722721951, 9387100.711479243, 9387189.20231424, 9387253.04151081, 9387298.458652133, 9389617.94632632, 9392402.283774145, 9394149.975087762, 9401037.95764324, 9410630.017878864, 9411945.670703603, 9412532.328363247, 9413253.571260372, 9414887.256361376, 9420804.6875, 9421105.89619378, 9421107.69383939, 9421124.681637691, 9433381.917005433, 9434098.58426839, 9434345.726710493, 9434370.3125, 9435955.104388442, 9436375.83997479, 9436376.497665443, 9442838.822168848, 9446480.700714383, 9446504.6875, 9450931.25, 9464278.976143913, 9465028.441702409, 9465869.01443272, 9465938.003995422, 9465971.875, 9466097.81185165, 9469341.595293041, 9469390.47860369, 9469809.704222381, 9469833.488171494, 9469838.467163462, 9470042.174673326, 9470124.51300583, 9470309.375, 9472103.674934672, 9481967.1875, 9484439.36304169, 9484616.671628812, 9485519.390108004, 9486049.022658782, 9487785.9375, 9489252.328555718, 9507384.952241799, 9513228.125, 9523429.60364719, 9524304.6875, 9524514.0625, 9526527.689923856, 9526564.507349363, 9526589.376933398, 9527387.5, 9527800.37000523, 9532643.75, 9533010.844629444, 9533819.592208184, 9534866.235476708, 9535036.700544922, 9535350.0, 9545030.321029177, 9545100.99637735, 9545227.134004977, 9546106.384401431, 9546498.165281378, 9547285.666247137, 9547300.0, 9548478.125, 9548645.310451858, 9548740.899207123, 9550031.34111165, 9551707.8125, 9552605.401446996, 9552605.851353938, 9553657.484965513, 9562054.56344031, 9562371.875, 9563629.6875, 9563693.67078937, 9565815.664124388, 9575085.257596007, 9577381.176103098, 9577963.906861138, 9584071.278416738, 9589389.82176476, 9590979.09258505, 9591040.386781741, 9592120.612660326, 9606948.56807796, 9606976.02125493, 9613446.85710067, 9621459.255553292, 9621836.785819778, 9626922.235126052, 9627074.65278923, 9628367.633027893, 9635205.239818064, 9635745.740496589, 9648885.9375, 9652495.3125, 9655807.8125, 9656141.845308196, 9656593.75, 9656609.08282424, 9656742.004677035, 9656922.422898443, 9657005.998099044, 9657039.0625, 9657192.051306926, 9657221.875, 9657260.929231947, 9657265.773876911, 9657266.149383234, 9657268.577441236, 9657370.081132371, 9657382.8125, 9657561.061781563, 9657573.616184356, 9657628.125, 9657641.35431838, 9657642.971053936, 9657643.031103864, 9657648.336379016, 9657654.6875, 9657671.165542131, 9657715.544216901, 9657799.472992625, 9657801.424135646, 9657859.375, 9657885.9375, 9657999.132716931, 9657999.7062979, 9658016.846393414, 9658033.42021039, 9658134.375, 9658158.699094893, 9658163.337940188, 9658173.2949372, 9658176.853542052, 9658183.25571347, 9658187.5, 9658197.272923008, 9658214.800368657, 9658223.992459789, 9658239.0625, 9658246.359825965, 9658255.63293468, 9658275.0, 9658278.144343846, 9658286.9834116, 9658301.5625, 9658321.875, 9658329.473593608, 9658332.6413758, 9658387.258843498, 9658395.265774716, 9658413.95845659, 9658416.071465353, 9658416.668557458, 9658418.657321084, 9658424.619076567, 9658464.0625, 9658468.75, 9658493.402294464, 9658522.828267436, 9658543.726739371, 9658546.93213835, 9658553.291084366, 9658559.375, 9658562.3143369, 9658568.29072699, 9658587.195126064, 9658619.73635291, 9658629.529785713, 9658660.389166312, 9658671.404018918, 9658689.0625, 9658690.625, 9658691.846906701, 9658700.96590776, 9658707.144443003, 9658708.848314336, 9658710.999314388, 9658747.402172104, 9658760.641002866, 9658787.5, 9658792.1875, 9658807.353921818, 9658810.98537383, 9658829.093074469, 9658843.644490767, 9658846.66323796, 9658926.5625, 9658952.560560957, 9659020.3125, 9659024.045865737, 9659039.0625, 9659063.717655929, 9659067.086802756, 9659178.612032536, 9659201.552739644, 9659215.22642859, 9659221.115147658, 9659290.341234531, 9659306.968595747, 9659320.759863622, 9659370.3125, 9659521.875, 9659546.967779955, 9659574.92076158, 9659595.3125, 9659601.781027209, 9659621.09103631, 9659626.146711342, 9659645.89125797, 9659650.742401293, 9659704.284199376, 9659743.267198969, 9659760.9375, 9659764.963497555, 9659801.15869171, 9659811.529305374, 9659834.550973622, 9659930.808898399, 9659936.577534705, 9659966.015420724, 9659970.3125, 9660145.255293535, 9660218.759104874, 9660311.362844318, 9660378.516949484, 9660434.375, 9660536.748454707, 9660584.847487042, 9660593.346921958, 9660708.422831573, 9661235.9375, 9661400.0, 9661465.625, 9662012.452187384, 9662189.992382219, 10426407.582043407, 10467749.378067244, 10540810.9375, 10543097.839392321, 10671911.547558827, 10717692.571378345, 10745472.12049515, 10748170.3125, 10748569.88804872, 10749989.264167875, 10750039.728189787, 10751219.836592996, 10751967.1875, 10753414.0625, 10753906.187087854, 10754763.858625816, 10754979.143617032, 10754991.08121027, 10755305.881249318, 10755352.31359794, 10755411.486346485, 10755852.02375274, 10756001.013996137, 10757195.068804407, 10758076.861832762, 10758547.463982968, 10758949.2931958, 10759239.392921679, 10759328.125, 10761004.213081356, 10761106.25, 10761150.930300463, 10761359.176743908, 10761412.5, 10761633.29782466, 10761872.19917027, 10762755.550410097, 10762799.915923996, 10763397.472332638, 10763780.936076377, 10763886.167680684, 10765110.823693387, 10767200.504363976, 10767678.791481636, 10768155.453256307, 10769157.785072489, 10769243.75, 10770129.6875, 10770204.13605203, 10770440.90994106, 10771073.179616235, 10771683.454448944, 10774063.398531575, 10774303.203784268, 10774348.979502395, 10774814.11509693, 10776589.0625, 10777628.109538162, 10791905.43411666, 10792365.118187916, 10802035.311484888, 10812466.119662225, ...], [69.77826560890375, 91.55267732843626, 9.48260086091381, 58.110938815141836, 12.18531375387578, 10.488654754911623, 13.757151043791694, 29.42451301356898, 18.14652392422138, 17.013698643403096, 50.135531786664764, 25.078090569469914, 66.15028514486116, 43.26247340860597, 45.22565552233857, 61.12638228105293, 45.25321925344876, 7.321984280432101, 42.97378178584038, 98.00823560613298, 86.50762485022932, 12.713119248037552, 25.07212283452826, 18.66793346000179, 21.849202220640354, 70.85218833703803, 59.9940620473035, 49.7914583735647, 5.467094748059871, 5.771918492542404, 117.35535142021838, 54.07610709568627, 78.81090848751654, 53.76496256509341, 37.820336784745145, 7.497595782644844, 50.19741674938565, 17.637126278536385, 30.183405516532815, 56.94435682094335, 87.81304566695417, 108.2196887180205, 80.03361576365157, 103.455102875924, 91.15444897749427, 25.252121254335393, 42.45770497107402, 53.12210812504012, 5.602499948103115, 23.551068338210577, 36.17739737275985, 49.25835630793095, 63.514911442422274, 54.27549545299744, 9.665480726269346, 67.51404859695468, 23.846389566681204, 69.93098190327427, 8.359168216362772, 17.725468508667205, 36.07742534822684, 41.05005251652426, 148.12882075413344, 21.306393983357516, 27.841073594507304, 70.31016320435519, 18.890571742234965, 5.760026617818309, 5.21229544370629, 53.237428798964054, 39.721450213582074, 96.65502152448519, 111.77604791726147, 31.021309875172797, 36.52866313107864, 56.55636521487171, 59.162723263305374, 5.545706046516259, 9.95423923825268, 51.959870332086965, 61.37472980352535, 75.55777427819783, 5.671908389441759, 7.098069513803744, 11.773988836274457, 15.288505632155287, 59.95914584710735, 11.052700217418792, 31.69626150940286, 46.58064869086881, 7.509058905494638, 35.28736434585281, 51.28882097167621, 18.884092905126163, 104.26727288438829, 59.37053935191965, 15.154396030690538, 64.14409251367024, 14.255488021367157, 102.19269292896851, 13.491385710057397, 82.96419980923928, 21.101392055129082, 74.09458481242889, 14.393799648652077, 6.4400771223857545, 24.240168078439353, 8.55926391178432, 19.96530706759422, 80.76747958587946, 28.080431140001963, 84.60793485446379, 30.188562014509337, 11.888317918385392, 74.8438662612712, 13.281001884818322, 88.96939309636943, 9.558897321027768, 22.023623449562486, 25.56263039213629, 5.719657638938461, 54.395002641971715, 16.11847539616461, 35.95064151530201, 56.586228774140984, 6.339117604009785, 23.179439585480512, 21.396556550613766, 12.642559883862226, 8.482338739056356, 22.101307208431166, 10.423322274980073, 96.30515458100007, 102.65855846653412, 16.416944106737866, 69.92368849194668, 21.000230826399413, 27.355706835256676, 29.72659693089, 46.331622054894396, 5.083385012586598, 50.34345415570628, 39.901953554923, 8.690298301484715, 114.58837202716028, 69.25984415193179, 24.483093639604494, 13.233463316965059, 35.07533709740332, 38.601597606352726, 7.195471338449034, 12.638646440861073, 26.699677797946727, 7.190233169668652, 44.342754045220715, 15.296342450602864, 60.96775148683466, 98.14002113055867, 5.660679524503315, 5.908668401920862, 11.320122081198331, 5.466998104572451, 151.76435792944645, 29.650971082229056, 20.56394780745952, 17.926491000605367, 11.013858865676111, 87.82463792104352, 61.88881908172403, 10.852508559356604, 108.35878614638794, 27.713020986304016, 47.74661144678452, 113.33298046853396, 48.62698687851665, 25.084222448780782, 8.127945293613735, 21.572574124053528, 11.446825413902843, 11.835810773849381, 20.518112986889427, 9.28761822442976, 182.44865673331793, 72.09392643970429, 6.757078370448362, 37.53434549947089, 43.862184228860706, 12.11818223152035, 83.99641048909596, 6.488598538170047, 24.50078737351882, 22.414197838065576, 25.69164853956341, 5.965960585589358, 33.79061639437341, 39.394399632736864, 15.188309675610437, 58.27062158239014, 13.719541369234069, 12.491839798950364, 42.30669312942717, 33.1458119567504, 28.527440687348765, 80.5503299019194, 10.510184917078977, 26.112351842571165, 17.912004140621757, 38.85215284886852, 25.87608937406097, 11.010932863726408, 71.51963132412364, 57.70011033443569, 40.5382917285101, 77.36689761857832, 48.06874733154836, 73.7649713873403, 6.021590453166321, 35.28392573072457, 34.406308164104175, 28.84738541322323, 5.724146255085033, 5.596267675305432, 9.233385842603509, 9.818104164251224, 31.90194260980565, 32.38416776090794, 23.785957854336807, 99.63791219796497, 8.37336435543623, 48.714538973921876, 95.50427948263514, 27.504528514206136, 21.808520984198363, 17.279396450051234, 61.16769341308694, 83.41373918909656, 113.53813997056459, 23.859657102399453, 54.29486651058716, 39.92188747251334, 25.54224565423921, 23.567788374344858, 57.76254494149095, 14.124588759434356, 52.83092875570432, 51.82801483736395, 10.259721238681559, 74.780426753657, 25.54930961553815, 16.247641292194967, 9.718041601855553, 62.315855927097004, 27.46356148824193, 14.575677104897352, 20.255432540843824, 69.65170981410218, 69.30552318905077, 110.70010281924421, 69.9128692736494, 92.51425395888629, 126.56316952612504, 146.5513798027724, 22.120894035622314, 6.927348354020946, 56.689154665369486, 5.196868511760627, 23.239373645660116, 92.00801747574113, 9.472121651408102, 29.708178673980022, 20.967126254218723, 35.00287511331956, 78.58197679064716, 27.95378145571374, 6.5503867777059686, 11.082052860565552, 34.48474351626413, 46.40021736496004, 37.05272406409659, 41.40087410083035, 23.102873416494667, 18.05653922382973, 8.408888953239286, 6.784875430665667, 32.67927096230147, 7.114597481842072, 78.85094021770959, 65.19907024630844, 21.175535193090656, 113.62528944644407, 7.445453559130028, 13.255953163892276, 92.4478821070038, 8.752828319539226, 62.47474399185407, 14.206312089854476, 19.636527769139082, 57.42691269989697, 58.033160042405996, 89.96255001832894, 52.39957507407625, 39.612089620570245, 28.292023645302557, 6.870381315887813, 92.8890883839816, 22.844380161163226, 37.685547258308354, 22.696896710220308, 17.280384425796413, 5.51760618041778, 65.01247399981318, 18.970086697214956, 29.614079755906488, 30.88702771837743, 30.19455806839862, 38.9873339125901, 41.43169696885861, 34.56333666122596, 8.136524010449943, 24.84766113885159, 52.50079941140823, 42.19547848925977, 92.33245360661022, 20.053301029792472, 107.04747762077136, 93.6232831082136, 50.64458534228507, 53.95880776982894, 14.825140580137498, 19.205732972225306, 6.879785543755798, 9.188760695686844, 40.16983503560883, 33.22950870243038, 28.475210282620367, 67.63021431687542, 109.71960166266305, 14.786023503875086, 45.81227220384923, 41.93690318060415, 96.82889999071415, 51.82011710029744, 15.077848967137742, 16.42255082775159, 30.579413936515525, 42.88533437386914, 71.50555989274369, 8.931110209888228, 17.698367734420394, 16.72897683872372, 7.9765774331894574, 31.0869308261036, 47.57458007513725, 36.20160207781897, 5.381380488049234, 59.32357859321691, 6.450159446990361, 13.247988546154946, 7.185652170532421, 5.861313367581496, 6.764971436503363, 74.22611078713354, 23.825951508019, 6.987285139046644, 74.02427294602852, 95.7766083669544, 18.2134675453204, 27.469551590912438, 10.808488974648137, 112.7176341480447, 7.384341120742907, 35.651358806444996, 7.197885467963935, 45.074030655385855, 40.399279351448925, 70.86625677576966, 6.008419338358145, 20.698578336306483, 22.088825287063038, 54.46707939245381, 39.440459336501974, 5.7326166801739395, 33.67928179453982, 7.7452045615218115, 86.0593088256797, 6.874310394906235, 26.793549312090107, 25.462368443591785, 30.069386269026083, 12.44534051936803, 33.328928890411404, 29.46860812153487, 45.51092059449774, 25.128617533951314, 41.471124129808146, 28.708237989294137, 41.33286033423291, 16.9150463847392, 20.104087066646258, 6.366323736958948, 28.620024205295294, 8.828264785054541, 8.284339114761906, 14.203029672649995, 42.657924876798106, 55.486539704723974, 6.25441858712848, 7.963832727428369, 30.67098534775709, 42.615415609572125, 25.780350975092084, 32.76428229249645, 22.97878734531394, 29.304694439916226, 9.15004322425323, 64.79143367109447, 14.667571965803345, 13.676500641371305, 25.540463303240664, 16.85400792516345, 24.08524227352808, 22.97989556067756, 8.36593063420316, 53.173289804208565, 65.77447399368772, 160.4507280206815, 85.25189631004315, 14.660170047425588, 11.019092754691924, 36.64234042926972, 52.166359098954416, 131.44005705698146, 71.95922374674139, 109.27659166187777, 8.12538994839134, 24.207479589659716, 54.49864234514353, 8.394215456708961, 13.29099188085044, 87.25574374227118, 57.14132893996474, 29.99940454200873, 49.87084521061469, 16.446007253148593, 14.831139371859233, 76.65810104222564, 13.967785142642757, 21.115148364671953, 48.64735902389387, 6.62970676344421, 71.31736722240434, 12.468233825043775, 6.513137231132816, 61.69115976618602, 54.69081165622625, 32.963092682298196, 30.384666757174628, 27.89930970396219, 23.069261429950387, 5.9602205989615396, 34.83156137260096, 7.910088134013308, 17.359875859538437, 62.45234644551324, 16.48973797589025, 48.70373538989847, 11.266531131759114, 80.75471036608904, 71.77061776585, 9.857928625729334, 66.08535707538118, 49.848278498878656, 6.357694400854463, 87.06081404735941, 70.6627637770713, 32.03186611556872, 10.489117919340877, 153.46168328427834, 80.82039310971457, 14.22023405971005, 13.184687899079929, 73.29301168875114, 31.14870292714661, 19.363197210734917, 63.565408336594935, 11.864680640969805, 18.551292596396358, 7.036312773519994, 22.238418985422385, 18.506631252375602, 12.91634100634451, 70.94886765954244, 94.91807868698646, 7.381557464355967, 110.69531224531683, 106.34017512554217, 15.723031820776056, 29.409838891722575, 82.77878549863067, 22.54241075026417, 79.56421189125302, 122.9980516510281, 59.604921514839035, 5.714377217766062, 5.081274727558386, 23.629770790920563, 16.891040997691864, 11.51688395950215, 33.68587371120256, 17.8043541830509, 45.32292703450922, 20.681318307821012, 51.75598295007495, 63.377263088492185, 14.664060544467723, 17.436911749376428, 48.8104286884791, 31.467838187166503, 26.048113668271522, 26.273024039091503, 89.67299082683479, 27.515402816751955, 88.3901826370805, 26.29814146026113, 45.49511398434288, 44.90524951118547, 26.532478692975445, 55.2428955488178, 49.793403544275755, 50.87569320749922, 62.03314406369785, 78.97821755623272, 15.304778900270604, 62.201247421478456, 19.99745237598238, 38.52839964422239, 21.54917034216446, 103.51565854414545, 17.64156608598151, 145.4593856524358, 26.657779712734758, 19.519720315460518, 127.06404601446914, 75.36540517033468, 8.934913916819317, 9.90697976111298, 6.7092045550091095, 8.87224021382223, 38.18152125325919, 12.474097126212813, 10.976902718818442, 23.49248926413497, 44.047635883158094, 81.85640687175498, 5.726070875819977, 28.746029952540848, 30.876854950818323, 17.55746724957502, 62.144187615252555, 44.528774948730245, 14.481407248395923, 6.168933455879589, 12.815617864883375, 59.153207612370885, 13.691497122321188, 101.9284440798304, 22.347845293845147, 55.238201054080335, 28.481461562771454, 70.35314407269958, 13.092168474508384, 51.36845124333529, 41.721283775918565, 93.46683385249949, 20.800513050628126, 18.777776585600243, 5.90747328807188, 27.396136088327772, 18.647158631583117, 74.14988856292834, 14.95134793238071, 46.412787311015656, 8.532993988326437, 8.064661307903851, 29.40156838792065, 5.623419667498884, 59.53769571011693, 15.85178728808059, 23.097007536979238, 13.230820041497793, 11.67991406003155, 32.793827500641655, 23.680092722443156, 28.487469785378735, 34.22660152391056, 60.51935085050796, 64.2696727230055, 76.57524488209904, 9.873797679732139, 9.591154737291168, 95.49066297627802, 12.633821111497921, 46.6190565844058, 18.660298562787762, 31.123500360639106, 69.3944357263281, 45.91405222458597, 20.39776673615765, 20.38151942815262, 10.113691739717485, 84.38185358806999, 15.638609193024253, 8.122394366764263, 20.168226361937545, 27.822636818986375, 76.25945900789, 25.713568021890932, 18.59689613864327, 17.130666596408332, 26.688964911396912, 77.54017381426794, 9.264304194180871, 73.02081844012453, 59.18853490417307, 82.09590632966399, 48.51829706190693, 13.80238125734703, 43.93863612471406, 33.20971224524206, 19.89200826635866, 7.785847241239852, 48.60158334400493, 40.13477631564674, 6.855292216660307, 6.428238907641042, 83.47370864726389, 21.87070321216279, 110.16524089744877, 100.9008443695833, 33.15516884757169, 54.6271857424183, 14.282822637060761, 27.06758609068604, 5.583378603177801, 21.329669333758147, 23.80626062898029, 19.848780269837988, 29.171047211904916, 117.41566427668782, 80.26145077099804, 56.06546696714246, 8.219672710772471, 6.7163851879354945, 44.76630773488396, 30.87254847041948, 111.62958020851474, 60.05537100024193, 16.149333841698024, 60.05275947766313, 17.78581064773958, 6.367064299652842, 77.2479352983438, 16.642855823419023, 14.401007980535407, 23.182032338039065, 23.550792186263195, 26.186780864506463, 103.17134443433095, 34.59514545535535, 40.7988887733481, 133.09142995178235, 54.6440877361341, 46.96550380001123, 6.410374507676841, 5.24341684598228, 13.790303806704257, 16.756583799575893, 38.84844856884651, 54.653141264019645, 103.03755278251344, 29.421713828057616, 19.01905487800498, 11.276303739527185, 63.65566606747419, 35.09023204787191, 41.8347777387739, 25.522313355936795, 49.3888453456664, 7.543579133817603, 5.465633693719902, 5.271355276151485, 18.585010653915838, 41.987020887939195, 33.39544699610464, 5.098872764668381, 6.613008175899193, 33.50148939358008, 11.088989365801158, 11.927247712382789, 11.686608334973291, 25.161520452523277, 21.65281772680007, 7.773031395135736, 24.498405631440903, 58.52275196722579, 19.98488787771426, 6.468649639701314, 13.852649673103077, 39.530055844834955, 11.365135482619827, 133.76187876718976, 84.7384813813763, 6.554490104850635, 7.938191585915324, 5.809198587663047, 104.97410782690235, 13.19306113167598, 67.41521014959426, 132.89172246295493, 37.32232045216823, 106.96845039403092, 166.7920621873063, 7.796296673563342, 62.25816501533568, 10.272871640949646, 94.88455433207378, 10.906338731637183, 26.51085657130392, 12.378305105214896, 61.133276583946135, 73.06133011517392, 21.623828901995033, 130.73248346426675, 12.111854335398787, 64.58310030372296, 7.202761413681687, 6.66470913698336, 7.300536287094914, 13.66260938967418, 31.36215478048632, 69.08257636658377, 29.91356266487371, 147.20001365691377, 43.287191073634915, 35.98030459908736, 51.818181167478414, 23.306662844659694, 16.959743958596363, 20.320186323756786, 82.49483306967863, 23.87802323874675, 33.21173986036255, 84.65118678111139, 10.161947438066907, 5.9943931676093705, 18.94886324706177, 63.275757943643846, 113.37240814371873, 86.82889063994705, 17.682905557180604, 5.715365398316607, 51.9313193529557, 24.125025749365097, 67.10846952277934, 86.33003899555197, 11.124506932028083, 7.705518793252499, 11.273170528706327, 40.778511012024815, 92.25678238781433, 22.567473511396813, 52.80326704158034, 14.488835629977112, 66.7529617386063, 49.069312322289875, 56.8709807779707, 5.246241624991588, 19.868408961669093, 175.26049278450168, 20.236829332667504, 20.580859377899355, 18.921560029395177, 56.64419763233434, 98.26725037637073, 11.631693985190036, 11.338132534469736, 49.851144234986236, 21.342413660703365, 53.72032388418809, 19.593143817440968, 11.741479634264495, 106.513350608544, 5.305050581674112, 51.21543539192766, 30.365448796487218, 62.6179015062853, 82.33111529571401, 38.38677243626686, 19.280299440029765, 62.04096657765081, 10.699792804440605, 58.03534618647714, 25.20306004959948, 79.94184810960904, 124.90520617156518, 9.218839372630129, 47.05321074380154, 49.68121313762114, 46.17582538454992, 17.404379228090576, 25.52761709067381, 6.718178458050407, 55.74669613189269, 7.154509276863198, 11.418980654067878, 72.39913812446233, 25.71048634471073, 12.621917317343861, 5.225460690407179, 29.16689130934563, 52.82118755875637, 64.37664063591447, 60.43871007587107, 13.252682827118646, 83.25972817088461, 104.5110256013063, 84.77956184298004, 18.47399121285312, 23.831523349094752, 10.633346382443696, 17.90510779735961, 50.32705636572011, 6.386111902768931, 43.140198657118546, 22.913514130900406, 5.474398605351051, 25.20467608790601, 37.21375192816757, 37.60414558583007, 8.719542550059671, 150.90357192146126, 31.65445329142686, 8.031721305098555, 13.016472551202959, 37.49072403505779, 81.54250882004574, 64.140469946995, 39.446260516502456, 79.92310192098643, 14.184403899686975, 125.00160504645694, 60.27709883940308, 14.766355922354824, 7.566035487226929, 11.90548248748581, 6.20378697940843, 26.817275631840847, 22.769311485704314, 51.53068194112272, 32.716439998202574, 13.763438673215541, 6.970011162992154, 28.072632619209237, 12.682655047810279, 15.946165161397058, 66.28638209950796, 25.33004430235787, 58.43077602174903, 23.68868289961794, 10.332587592674747, 5.336546280758336, 31.041908306121343, 86.92095281556595, 57.38300761848095, 30.09479958875559, 133.33951786138735, 35.010125498291586, 65.72176873422268, 52.59428477568801, 15.1824701921035, 26.04619362604328, 157.36656400576214, 87.87778969766059, 74.81181180318475, 28.94930427900464, 5.535408529767064, 13.074668224043917, 6.2539594389111315, 9.1585429445389, 66.5214701444361, 25.4385661763878, 61.47616848217645, 12.2334460682809, 37.95016324449648, 70.91680563317215, 26.680775348307794, 30.700660686446525, 28.333233915575754, 22.568314438068622, 5.497304982047925, 22.044268297984786, 60.21149368833324, 6.325423978666722, 158.06260499281035, 56.245033666554306, 5.465762899336625, 27.049859792403083, 54.16250030105062, 24.866674301939153, 7.3423362838986455, 9.002948551042019, 17.72683426307021, 10.049375478318348, 8.084229999503211, 56.06073486591909, 56.174026562532624, 42.193191136692036, 6.685481002078481, 17.806589890005085, 19.37101629501525, 22.101477452675635, 18.109836629526605, 94.71065431426054, 41.26115423990972, 231.2090986539324, 7.606534083473084, 9.718960799781552, 9.686370344131078, 53.22109012810866, 14.871360304667427, 36.4910188154177, 40.67924926450571, 14.843721351013137, 32.70807733900208, 90.70543992461603, 77.86839180027894, 14.762348830490495, 94.54933088113195, 9.143855656068174, 30.957854703081367, 45.91616662957094, 53.07797685442267, 5.422319615225528, 16.638852230915695, 5.7225621590299305, 81.55801831222949, 69.34543000426277, 116.174233207947, 5.052943386685845, 6.20888682727358, 57.07011504416968, 61.30565866995507, 27.338937424105687, 7.531127503764302, 58.26624235704098, 6.56106597321113, 22.71922041011763, 18.862919954016824, 5.217634064397359, 17.04895311596002, 13.19513978403063, 61.39846085789381, 5.726076903558251, 65.71865766160028, 12.793227880158955, 16.37063839248726, 58.10246394639079, 10.459158889853171, 41.577994702543215, 61.210382709045156, 81.16742583011197, 59.327075090281234, 11.98169076938157, 45.577654238132844, 53.06708763869012, 45.96337653185522, 48.113244558365565, 23.574501193124455, 28.12615448274729, 10.954360261608782, 65.07561663177394, 7.863496720619498, 12.580716471574828, 36.3830052370592, 64.13341800277311, 29.623305339205345, 14.192461473662146, 8.443234665134106, 8.15475500420226, 5.423595598713941, 7.964121563174872, 18.01102258969498, 23.920927179573255, 15.145847000260014, 50.91316388665784, 11.10362465561797, 10.054041913523104, 59.318671319467505, 61.497372171578625, 6.866450655096844, ...])
caption, T_hit,count1 = hits_in_time_hist_new(T_pom, dt, t_plasma_start, t_plasma_end, is_plasma, figure_count_in_time_hist)
hist_file(count_in_time_hist, T_hit, count1, caption);
([4183873.4375, 4200589.0625, 4245645.486060796, 4245648.321076619, 4273347.460452169, 4314272.98618653, 4463036.308541708, 4516373.4375, 4519480.333225946, 4520039.226822126, 4520225.0, 4520815.422263485, 4525556.211428195, 4525936.918065219, 4526217.1875, 4526593.75, 4529406.065605748, 4529410.086621125, 4580491.832632525, 4586693.342413341, 4595065.625, 4596572.697327951, 4604054.484354057, 4711086.660961539, 4722958.058993217, 4810845.3125, 4814384.375, 4891162.5, 6695396.9985546, 6866860.352375621, 6868261.685191906, 6994003.125, 7099131.25, 7129463.278797314, 7130623.073756609, 7130945.010654807, 8185245.3125, 8525754.387882605, 8628820.776659146, 8696595.68875797, 8701679.491570348, 8703296.875, 8703743.628683295, 8705085.336426139, 8705454.924390363, 8706293.053784892, 8706748.4375, 8730331.25, 8731529.476926165, 8748865.30920586, 8749673.4375, 8750856.25, 8779617.1875, 8790217.1875, 8825184.74532947, 8834315.625, 8853562.445120098, 8863321.875, 8864176.697293386, 8886742.76603319, 8892519.140832322, 8904692.1875, 8904913.658927286, 8905095.020341542, 8905349.964492438, 8905650.0, 8918130.442457212, 8919474.243177438, 8922011.282267557, 8928745.3125, 8931389.0625, 8937252.011676298, 8947414.0625, 8952057.307835227, 8965704.685666695, 8968389.0625, 8969540.953686716, 8970240.882425923, 8970963.154934406, 8971078.125, 8975503.980992258, 8975887.5, 8976480.68812859, 8976954.5977973, 8977042.565043097, 8980778.213154148, 8980976.5625, 8980985.456670586, 8981047.939277556, 8981551.095397804, 8981751.29438005, 8981798.4375, 8982078.030210437, 8982135.129054297, 8982614.469613882, 8982882.12367435, 8983096.929754166, 8983256.25, 8983643.557566142, 8983773.4375, 8987637.05927988, 8988124.254688518, 8988876.295442553, 8996076.24896725, 8997781.09847637, 8999818.443383655, 9000604.424873259, 9001122.357356636, 9001125.894374518, 9001335.9375, 9001542.16464767, 9004281.3693377, 9005364.941423027, 9005503.562604766, 9005923.390236454, 9006085.420970289, 9006089.0625, 9006169.007326934, 9006280.44515653, 9006329.517238103, 9006368.950798277, 9006621.719176397, 9006647.152944328, 9006790.625, 9006796.875, 9006884.923797358, 9007031.321908507, 9007034.372597642, 9007046.104904564, 9007155.040208349, 9007167.259587096, 9007182.30970794, 9007227.141970402, 9007312.964309137, 9007439.402320221, 9007560.507849539, 9007670.558947109, 9007743.023725435, 9007760.9375, 9007764.0625, 9007819.531306375, 9008027.550515702, 9008043.839343544, 9008055.380151317, 9008062.5, 9008104.6875, 9008143.505249066, 9008199.465843676, 9008345.143207641, 9008404.6875, 9008411.211793315, 9008517.924231231, 9008560.837918967, 9008656.505350163, 9008714.0625, 9008840.750927528, 9008893.75, 9009074.880919412, 9009133.625789277, 9009182.45264877, 9009201.682209646, 9009321.996323926, 9009331.541482333, 9009363.358914142, 9009468.169689404, 9009491.236132098, 9009671.800921798, 9009751.5625, 9009754.405385567, 9009768.69761471, 9009918.39377025, 9009929.048365531, 9009946.984452512, 9009953.125, 9010026.5625, 9010035.73518393, 9010084.06247413, 9010100.137702378, 9010101.757499326, 9010116.035248041, 9010130.758633042, 9010135.885259666, 9010141.002152685, 9010151.5625, 9010167.15521591, 9010210.9375, 9010212.5, 9010222.428020593, 9010326.869958924, 9010334.384608725, 9010350.248461494, 9010354.845121, 9010377.016465543, 9010404.987460053, 9010425.872245127, 9010455.531909695, 9010478.18827834, 9010510.9375, 9010589.786125397, 9010597.601751687, 9010634.375, 9010642.1875, 9010657.8125, 9010668.75, 9010671.440412492, 9010747.83882921, 9010783.42131511, 9010807.8125, 9010813.912468296, 9010868.673984608, 9010883.92879088, 9010905.421934988, 9010912.5, 9010929.6875, 9010981.25, 9011003.32062525, 9011030.585197488, 9011042.1875, 9011107.8125, 9011150.0, 9011244.019793529, 9011268.424141327, 9011268.631759254, 9011300.524493318, 9011330.166647157, 9011347.483890036, 9011403.785153715, 9011485.086929744, 9011520.473087369, 9011555.373735217, 9011576.543449743, 9011648.38375643, 9011740.095098844, 9011748.940079143, 9011795.3125, 9011848.4375, 9011875.0, 9011901.270979028, 9011943.324927993, 9011951.5625, 9012105.207295107, 9012110.623040995, 9012123.06182649, 9012215.392783487, 9012225.0, 9012557.8125, 9012846.238316586, 9012985.653734386, 9013038.183354102, 9013056.555102466, 9013261.361907978, 9013425.0, 9015011.177805066, 9017948.338090392, 9019175.928783504, 9023421.875, 9024487.5, 9030298.4375, 9061534.749423213, 9063212.96730573, 9063230.76319783, 9063521.57554675, 9063547.998211809, 9063607.018780848, 9063667.1875, 9063702.865812365, 9063825.550099153, 9064021.875, 9064248.59936124, 9064320.3125, 9064390.070606042, 9064467.140083848, 9064715.280235466, 9064726.532983446, 9065384.767233357, 9066142.158094348, 9066198.905409403, 9066262.5, 9066264.0625, 9066506.25, 9066533.510633184, 9066636.528919457, 9066854.911846513, 9066909.48110164, 9066939.0625, 9067093.723861238, 9067145.3125, 9067179.6875, 9067247.905065252, 9067314.434433706, 9067383.916610664, 9067411.974206004, 9067421.875, 9067463.288606301, 9067533.20841679, 9067562.292744435, 9067566.478803854, 9067648.344025167, 9067799.012766793, 9067979.6875, 9068041.77729005, 9068117.158770138, 9068486.447622718, 9068523.950253056, 9068767.1875, 9068840.62600095, 9070351.527994512, 9070593.115576735, 9071427.604652686, 9071523.114482097, 9071887.5, 9071983.578473588, 9071989.0625, 9072018.812990367, 9072409.375, 9073620.3125, 9073904.609205527, 9079628.08675259, 9079909.080042688, 9080418.531281415, 9081234.571802793, 9081265.625, 9081870.721764483, 9085774.34244485, 9086248.4375, 9096315.625, 9109050.98628012, 9109132.8125, 9109457.781630544, 9109699.717481546, 9111366.180808999, 9113138.88884435, 9115875.0, 9131820.13874599, 9147349.445457477, 9159259.375, 9160656.140578115, 9165311.122044193, 9167652.231229564, 9171349.105878329, 9172940.356509535, 9173229.024531271, 9173242.222816251, 9178196.364252988, 9178356.25, 9178383.955048038, 9178865.625, 9179410.422848966, 9179777.13656335, 9179841.026321342, 9179868.113863625, 9180475.929992128, 9182109.375, 9182704.849608814, 9183585.483051082, 9184177.752318073, 9184590.385287365, 9185301.033749372, 9185478.363726486, 9195449.057354204, 9200950.00733304, 9200998.4375, 9201441.400513928, 9202943.2137763, 9203593.352078538, 9207244.239531096, 9207614.192748923, 9207905.513074651, 9215696.660701513, 9217120.3125, 9217249.351318192, 9221911.965753097, 9223551.845500352, 9224036.63207275, 9229127.409765972, 9230625.0, 9230674.197268892, 9232825.15064998, 9236192.24933635, 9236565.625, 9237231.951013215, 9237239.461285375, 9237421.875, 9238530.063670369, 9238602.326439526, 9238670.843280602, 9240362.405934347, 9240542.58173975, 9240676.542443758, 9241408.518504225, 9241526.5625, 9242412.5, 9242884.375, 9243115.425720569, 9243726.98515286, 9244298.4375, 9248301.5625, 9248686.374043081, 9249471.22349186, 9250530.433762833, 9250697.830740603, 9250924.338311244, 9251329.68326511, 9251943.54174985, 9252479.6875, 9253012.5, 9253053.023174386, 9253705.585383918, 9256217.219854118, 9257523.4375, 9258080.845368464, 9258095.3125, 9258101.199564034, 9259207.19113179, 9262720.089933092, 9263687.5, 9264140.551223606, 9264433.998159565, 9264440.453284826, 9265517.404212292, 9269219.783151336, 9271054.32465812, 9271418.897103176, 9273351.000810036, 9273642.1875, 9275489.53549253, 9283621.840782374, 9286755.262242535, 9290206.17938272, 9290893.366560556, 9291002.993811766, 9292003.091235429, 9292214.491689194, 9293006.25, 9293141.8697256, 9293164.915213862, 9293335.9375, 9293754.885646014, 9294297.924650637, 9294316.579789124, 9294490.625, 9294793.535139622, 9294964.0625, 9295060.699740462, 9295157.783238593, 9295595.3125, 9295821.593704496, 9295830.711258667, 9296057.8125, 9296547.728200698, 9296898.4375, 9296917.903748617, 9297036.101665298, 9297085.906624135, 9297206.25, 9297495.529910361, 9297789.0625, 9298067.155095145, 9298157.45720855, 9299418.985953556, 9300103.014735801, 9307010.15055015, 9307911.454167297, 9308871.875, 9314456.6044502, 9315615.625, 9315801.649536349, 9316337.91821196, 9316378.125, 9316900.563317673, 9317332.84588332, 9317701.900590431, 9317969.43420567, 9318466.029479621, 9318790.625, 9320993.62328599, 9330527.674047578, 9331041.438932614, 9331164.0625, 9331577.921955507, 9331852.5725616, 9332437.5, 9332986.884632403, 9334005.506223986, 9334146.875, 9334386.362836745, 9334438.206567025, 9334624.66590418, 9334994.884104067, 9335202.262524545, 9335471.218261806, 9336109.685655365, 9336599.564251667, 9337186.842468543, 9337300.0, 9337470.89909258, 9337773.628468541, 9338570.3125, 9338996.445688896, 9339157.411625843, 9339460.739591464, 9339620.072334826, 9340141.672891175, 9340228.24405665, 9340231.919399926, 9345273.127970248, 9345720.744557677, 9347929.92298877, 9349310.9375, 9353172.466441734, 9354263.507174898, 9357708.789581047, 9357737.5, 9358298.4375, 9359907.737752508, 9361405.221549703, 9365042.1875, 9366331.25, 9369052.985822245, 9369233.443840804, 9375553.125, 9378942.134353943, 9380621.875, 9381745.188794693, 9382939.980990546, 9383186.52889745, 9383194.724825159, 9383218.75, 9383400.42650697, 9383411.039667604, 9383467.1875, 9383518.053417074, 9383518.842151439, 9383548.4375, 9383571.210749764, 9383574.991915066, 9383624.505602453, 9383626.450180985, 9383649.002135746, 9383649.89990699, 9383696.17235585, 9383705.526467063, 9383818.234928422, 9383901.681584785, 9383919.80316099, 9383956.860374065, 9383959.09677875, 9383990.026998747, 9384007.108439788, 9384019.46886308, 9384068.651329206, 9384074.753014965, 9384084.375, 9384087.5, 9384098.736740664, 9384112.5, 9384129.736935025, 9384133.363028757, 9384133.690564953, 9384135.9375, 9384157.68631216, 9384213.243030472, 9384218.05240922, 9384225.0, 9384265.154551305, 9384312.5, 9384326.505917348, 9384341.50534859, 9384356.247788457, 9384357.8125, 9384360.349823713, 9384368.558008399, 9384411.21654247, 9384426.5625, 9384466.588870592, 9384488.239121938, 9384507.438372737, 9384523.37776976, 9384539.782972625, 9384552.334091563, 9384603.127578856, 9384617.1875, 9384664.501116205, 9384697.992875282, 9384710.9375, 9384728.293236211, 9384767.1875, 9384775.346470393, 9384785.584538156, 9384794.777358059, 9384819.076719295, 9384823.4375, 9384837.571506204, 9384851.937856372, 9384873.915697863, 9384874.869817106, 9384892.1875, 9384895.3125, 9384895.891132612, 9384897.16701084, 9384946.875, 9385003.913663138, 9385007.8125, 9385008.5348641, 9385009.15462192, 9385018.396866033, 9385029.6875, 9385058.33446472, 9385102.959184108, 9385116.42490072, 9385121.875, 9385131.42161185, 9385144.988563828, 9385146.231073383, 9385159.57488966, 9385181.872342804, 9385187.339553459, 9385200.713205244, 9385259.851077123, 9385262.399786592, 9385337.326968024, 9385345.232064677, 9385353.081587017, 9385374.495866934, 9385380.08281271, 9385404.6875, 9385409.735008465, 9385423.876734253, 9385426.5625, 9385430.572978152, 9385431.728160482, 9385441.63424752, 9385496.875, 9385505.130497972, 9385507.426505128, 9385535.9375, 9385567.710546523, 9385568.75, 9385590.625, 9385625.0, 9385634.375, 9385638.90198532, 9385639.308358315, 9385647.869562797, 9385699.483459262, 9385748.141833985, 9385750.880192107, 9385760.646970747, 9385771.875, 9385783.480476495, 9385796.875, 9385801.433802659, 9385802.884297783, 9385818.75, 9385826.5625, 9385839.265898142, 9385868.691069163, 9385875.284486553, 9385877.48968759, 9385906.578791907, 9385910.126563957, 9385914.0625, 9385929.324142477, 9385973.288077647, 9385977.899757428, 9386001.24668379, 9386042.881670369, 9386054.649106922, 9386090.625, 9386090.625, 9386120.120026069, 9386151.379667, 9386178.125, 9386226.055003846, 9386246.813125495, 9386312.1611892, 9386314.015267078, 9386343.156216921, 9386375.0, 9386403.711666008, 9386417.1875, 9386422.647698505, 9386453.218032125, 9386537.5, 9386673.4375, 9386690.625, 9386758.824125618, 9386826.5625, 9386971.706984648, 9387017.277290575, 9387033.722721951, 9387100.711479243, 9387189.20231424, 9387253.04151081, 9387298.458652133, 9389617.94632632, 9392402.283774145, 9394149.975087762, 9401037.95764324, 9410630.017878864, 9411945.670703603, 9412532.328363247, 9413253.571260372, 9414887.256361376, 9420804.6875, 9421105.89619378, 9421107.69383939, 9421124.681637691, 9433381.917005433, 9434098.58426839, 9434345.726710493, 9434370.3125, 9435955.104388442, 9436375.83997479, 9436376.497665443, 9442838.822168848, 9446480.700714383, 9446504.6875, 9450931.25, 9464278.976143913, 9465028.441702409, 9465869.01443272, 9465938.003995422, 9465971.875, 9466097.81185165, 9469341.595293041, 9469390.47860369, 9469809.704222381, 9469833.488171494, 9469838.467163462, 9470042.174673326, 9470124.51300583, 9470309.375, 9472103.674934672, 9481967.1875, 9484439.36304169, 9484616.671628812, 9485519.390108004, 9486049.022658782, 9487785.9375, 9489252.328555718, 9507384.952241799, 9513228.125, 9523429.60364719, 9524304.6875, 9524514.0625, 9526527.689923856, 9526564.507349363, 9526589.376933398, 9527387.5, 9527800.37000523, 9532643.75, 9533010.844629444, 9533819.592208184, 9534866.235476708, 9535036.700544922, 9535350.0, 9545030.321029177, 9545100.99637735, 9545227.134004977, 9546106.384401431, 9546498.165281378, 9547285.666247137, 9547300.0, 9548478.125, 9548645.310451858, 9548740.899207123, 9550031.34111165, 9551707.8125, 9552605.401446996, 9552605.851353938, 9553657.484965513, 9562054.56344031, 9562371.875, 9563629.6875, 9563693.67078937, 9565815.664124388, 9575085.257596007, 9577381.176103098, 9577963.906861138, 9584071.278416738, 9589389.82176476, 9590979.09258505, 9591040.386781741, 9592120.612660326, 9606948.56807796, 9606976.02125493, 9613446.85710067, 9621459.255553292, 9621836.785819778, 9626922.235126052, 9627074.65278923, 9628367.633027893, 9635205.239818064, 9635745.740496589, 9648885.9375, 9652495.3125, 9655807.8125, 9656141.845308196, 9656593.75, 9656609.08282424, 9656742.004677035, 9656922.422898443, 9657005.998099044, 9657039.0625, 9657192.051306926, 9657221.875, 9657260.929231947, 9657265.773876911, 9657266.149383234, 9657268.577441236, 9657370.081132371, 9657382.8125, 9657561.061781563, 9657573.616184356, 9657628.125, 9657641.35431838, 9657642.971053936, 9657643.031103864, 9657648.336379016, 9657654.6875, 9657671.165542131, 9657715.544216901, 9657799.472992625, 9657801.424135646, 9657859.375, 9657885.9375, 9657999.132716931, 9657999.7062979, 9658016.846393414, 9658033.42021039, 9658134.375, 9658158.699094893, 9658163.337940188, 9658173.2949372, 9658176.853542052, 9658183.25571347, 9658187.5, 9658197.272923008, 9658214.800368657, 9658223.992459789, 9658239.0625, 9658246.359825965, 9658255.63293468, 9658275.0, 9658278.144343846, 9658286.9834116, 9658301.5625, 9658321.875, 9658329.473593608, 9658332.6413758, 9658387.258843498, 9658395.265774716, 9658413.95845659, 9658416.071465353, 9658416.668557458, 9658418.657321084, 9658424.619076567, 9658464.0625, 9658468.75, 9658493.402294464, 9658522.828267436, 9658543.726739371, 9658546.93213835, 9658553.291084366, 9658559.375, 9658562.3143369, 9658568.29072699, 9658587.195126064, 9658619.73635291, 9658629.529785713, 9658660.389166312, 9658671.404018918, 9658689.0625, 9658690.625, 9658691.846906701, 9658700.96590776, 9658707.144443003, 9658708.848314336, 9658710.999314388, 9658747.402172104, 9658760.641002866, 9658787.5, 9658792.1875, 9658807.353921818, 9658810.98537383, 9658829.093074469, 9658843.644490767, 9658846.66323796, 9658926.5625, 9658952.560560957, 9659020.3125, 9659024.045865737, 9659039.0625, 9659063.717655929, 9659067.086802756, 9659178.612032536, 9659201.552739644, 9659215.22642859, 9659221.115147658, 9659290.341234531, 9659306.968595747, 9659320.759863622, 9659370.3125, 9659521.875, 9659546.967779955, 9659574.92076158, 9659595.3125, 9659601.781027209, 9659621.09103631, 9659626.146711342, 9659645.89125797, 9659650.742401293, 9659704.284199376, 9659743.267198969, 9659760.9375, 9659764.963497555, 9659801.15869171, 9659811.529305374, 9659834.550973622, 9659930.808898399, 9659936.577534705, 9659966.015420724, 9659970.3125, 9660145.255293535, 9660218.759104874, 9660311.362844318, 9660378.516949484, 9660434.375, 9660536.748454707, 9660584.847487042, 9660593.346921958, 9660708.422831573, 9661235.9375, 9661400.0, 9661465.625, 9662012.452187384, 9662189.992382219, 10426407.582043407, 10467749.378067244, 10540810.9375, 10543097.839392321, 10671911.547558827, 10717692.571378345, 10745472.12049515, 10748170.3125, 10748569.88804872, 10749989.264167875, 10750039.728189787, 10751219.836592996, 10751967.1875, 10753414.0625, 10753906.187087854, 10754763.858625816, 10754979.143617032, 10754991.08121027, 10755305.881249318, 10755352.31359794, 10755411.486346485, 10755852.02375274, 10756001.013996137, 10757195.068804407, 10758076.861832762, 10758547.463982968, 10758949.2931958, 10759239.392921679, 10759328.125, 10761004.213081356, 10761106.25, 10761150.930300463, 10761359.176743908, 10761412.5, 10761633.29782466, 10761872.19917027, 10762755.550410097, 10762799.915923996, 10763397.472332638, 10763780.936076377, 10763886.167680684, 10765110.823693387, 10767200.504363976, 10767678.791481636, 10768155.453256307, 10769157.785072489, 10769243.75, 10770129.6875, 10770204.13605203, 10770440.90994106, 10771073.179616235, 10771683.454448944, 10774063.398531575, 10774303.203784268, 10774348.979502395, 10774814.11509693, 10776589.0625, 10777628.109538162, 10791905.43411666, 10792365.118187916, 10802035.311484888, 10812466.119662225, ...], [69.77826560890375, 91.55267732843626, 9.48260086091381, 58.110938815141836, 12.18531375387578, 10.488654754911623, 13.757151043791694, 29.42451301356898, 18.14652392422138, 17.013698643403096, 50.135531786664764, 25.078090569469914, 66.15028514486116, 43.26247340860597, 45.22565552233857, 61.12638228105293, 45.25321925344876, 7.321984280432101, 42.97378178584038, 98.00823560613298, 86.50762485022932, 12.713119248037552, 25.07212283452826, 18.66793346000179, 21.849202220640354, 70.85218833703803, 59.9940620473035, 49.7914583735647, 5.467094748059871, 5.771918492542404, 117.35535142021838, 54.07610709568627, 78.81090848751654, 53.76496256509341, 37.820336784745145, 7.497595782644844, 50.19741674938565, 17.637126278536385, 30.183405516532815, 56.94435682094335, 87.81304566695417, 108.2196887180205, 80.03361576365157, 103.455102875924, 91.15444897749427, 25.252121254335393, 42.45770497107402, 53.12210812504012, 5.602499948103115, 23.551068338210577, 36.17739737275985, 49.25835630793095, 63.514911442422274, 54.27549545299744, 9.665480726269346, 67.51404859695468, 23.846389566681204, 69.93098190327427, 8.359168216362772, 17.725468508667205, 36.07742534822684, 41.05005251652426, 148.12882075413344, 21.306393983357516, 27.841073594507304, 70.31016320435519, 18.890571742234965, 5.760026617818309, 5.21229544370629, 53.237428798964054, 39.721450213582074, 96.65502152448519, 111.77604791726147, 31.021309875172797, 36.52866313107864, 56.55636521487171, 59.162723263305374, 5.545706046516259, 9.95423923825268, 51.959870332086965, 61.37472980352535, 75.55777427819783, 5.671908389441759, 7.098069513803744, 11.773988836274457, 15.288505632155287, 59.95914584710735, 11.052700217418792, 31.69626150940286, 46.58064869086881, 7.509058905494638, 35.28736434585281, 51.28882097167621, 18.884092905126163, 104.26727288438829, 59.37053935191965, 15.154396030690538, 64.14409251367024, 14.255488021367157, 102.19269292896851, 13.491385710057397, 82.96419980923928, 21.101392055129082, 74.09458481242889, 14.393799648652077, 6.4400771223857545, 24.240168078439353, 8.55926391178432, 19.96530706759422, 80.76747958587946, 28.080431140001963, 84.60793485446379, 30.188562014509337, 11.888317918385392, 74.8438662612712, 13.281001884818322, 88.96939309636943, 9.558897321027768, 22.023623449562486, 25.56263039213629, 5.719657638938461, 54.395002641971715, 16.11847539616461, 35.95064151530201, 56.586228774140984, 6.339117604009785, 23.179439585480512, 21.396556550613766, 12.642559883862226, 8.482338739056356, 22.101307208431166, 10.423322274980073, 96.30515458100007, 102.65855846653412, 16.416944106737866, 69.92368849194668, 21.000230826399413, 27.355706835256676, 29.72659693089, 46.331622054894396, 5.083385012586598, 50.34345415570628, 39.901953554923, 8.690298301484715, 114.58837202716028, 69.25984415193179, 24.483093639604494, 13.233463316965059, 35.07533709740332, 38.601597606352726, 7.195471338449034, 12.638646440861073, 26.699677797946727, 7.190233169668652, 44.342754045220715, 15.296342450602864, 60.96775148683466, 98.14002113055867, 5.660679524503315, 5.908668401920862, 11.320122081198331, 5.466998104572451, 151.76435792944645, 29.650971082229056, 20.56394780745952, 17.926491000605367, 11.013858865676111, 87.82463792104352, 61.88881908172403, 10.852508559356604, 108.35878614638794, 27.713020986304016, 47.74661144678452, 113.33298046853396, 48.62698687851665, 25.084222448780782, 8.127945293613735, 21.572574124053528, 11.446825413902843, 11.835810773849381, 20.518112986889427, 9.28761822442976, 182.44865673331793, 72.09392643970429, 6.757078370448362, 37.53434549947089, 43.862184228860706, 12.11818223152035, 83.99641048909596, 6.488598538170047, 24.50078737351882, 22.414197838065576, 25.69164853956341, 5.965960585589358, 33.79061639437341, 39.394399632736864, 15.188309675610437, 58.27062158239014, 13.719541369234069, 12.491839798950364, 42.30669312942717, 33.1458119567504, 28.527440687348765, 80.5503299019194, 10.510184917078977, 26.112351842571165, 17.912004140621757, 38.85215284886852, 25.87608937406097, 11.010932863726408, 71.51963132412364, 57.70011033443569, 40.5382917285101, 77.36689761857832, 48.06874733154836, 73.7649713873403, 6.021590453166321, 35.28392573072457, 34.406308164104175, 28.84738541322323, 5.724146255085033, 5.596267675305432, 9.233385842603509, 9.818104164251224, 31.90194260980565, 32.38416776090794, 23.785957854336807, 99.63791219796497, 8.37336435543623, 48.714538973921876, 95.50427948263514, 27.504528514206136, 21.808520984198363, 17.279396450051234, 61.16769341308694, 83.41373918909656, 113.53813997056459, 23.859657102399453, 54.29486651058716, 39.92188747251334, 25.54224565423921, 23.567788374344858, 57.76254494149095, 14.124588759434356, 52.83092875570432, 51.82801483736395, 10.259721238681559, 74.780426753657, 25.54930961553815, 16.247641292194967, 9.718041601855553, 62.315855927097004, 27.46356148824193, 14.575677104897352, 20.255432540843824, 69.65170981410218, 69.30552318905077, 110.70010281924421, 69.9128692736494, 92.51425395888629, 126.56316952612504, 146.5513798027724, 22.120894035622314, 6.927348354020946, 56.689154665369486, 5.196868511760627, 23.239373645660116, 92.00801747574113, 9.472121651408102, 29.708178673980022, 20.967126254218723, 35.00287511331956, 78.58197679064716, 27.95378145571374, 6.5503867777059686, 11.082052860565552, 34.48474351626413, 46.40021736496004, 37.05272406409659, 41.40087410083035, 23.102873416494667, 18.05653922382973, 8.408888953239286, 6.784875430665667, 32.67927096230147, 7.114597481842072, 78.85094021770959, 65.19907024630844, 21.175535193090656, 113.62528944644407, 7.445453559130028, 13.255953163892276, 92.4478821070038, 8.752828319539226, 62.47474399185407, 14.206312089854476, 19.636527769139082, 57.42691269989697, 58.033160042405996, 89.96255001832894, 52.39957507407625, 39.612089620570245, 28.292023645302557, 6.870381315887813, 92.8890883839816, 22.844380161163226, 37.685547258308354, 22.696896710220308, 17.280384425796413, 5.51760618041778, 65.01247399981318, 18.970086697214956, 29.614079755906488, 30.88702771837743, 30.19455806839862, 38.9873339125901, 41.43169696885861, 34.56333666122596, 8.136524010449943, 24.84766113885159, 52.50079941140823, 42.19547848925977, 92.33245360661022, 20.053301029792472, 107.04747762077136, 93.6232831082136, 50.64458534228507, 53.95880776982894, 14.825140580137498, 19.205732972225306, 6.879785543755798, 9.188760695686844, 40.16983503560883, 33.22950870243038, 28.475210282620367, 67.63021431687542, 109.71960166266305, 14.786023503875086, 45.81227220384923, 41.93690318060415, 96.82889999071415, 51.82011710029744, 15.077848967137742, 16.42255082775159, 30.579413936515525, 42.88533437386914, 71.50555989274369, 8.931110209888228, 17.698367734420394, 16.72897683872372, 7.9765774331894574, 31.0869308261036, 47.57458007513725, 36.20160207781897, 5.381380488049234, 59.32357859321691, 6.450159446990361, 13.247988546154946, 7.185652170532421, 5.861313367581496, 6.764971436503363, 74.22611078713354, 23.825951508019, 6.987285139046644, 74.02427294602852, 95.7766083669544, 18.2134675453204, 27.469551590912438, 10.808488974648137, 112.7176341480447, 7.384341120742907, 35.651358806444996, 7.197885467963935, 45.074030655385855, 40.399279351448925, 70.86625677576966, 6.008419338358145, 20.698578336306483, 22.088825287063038, 54.46707939245381, 39.440459336501974, 5.7326166801739395, 33.67928179453982, 7.7452045615218115, 86.0593088256797, 6.874310394906235, 26.793549312090107, 25.462368443591785, 30.069386269026083, 12.44534051936803, 33.328928890411404, 29.46860812153487, 45.51092059449774, 25.128617533951314, 41.471124129808146, 28.708237989294137, 41.33286033423291, 16.9150463847392, 20.104087066646258, 6.366323736958948, 28.620024205295294, 8.828264785054541, 8.284339114761906, 14.203029672649995, 42.657924876798106, 55.486539704723974, 6.25441858712848, 7.963832727428369, 30.67098534775709, 42.615415609572125, 25.780350975092084, 32.76428229249645, 22.97878734531394, 29.304694439916226, 9.15004322425323, 64.79143367109447, 14.667571965803345, 13.676500641371305, 25.540463303240664, 16.85400792516345, 24.08524227352808, 22.97989556067756, 8.36593063420316, 53.173289804208565, 65.77447399368772, 160.4507280206815, 85.25189631004315, 14.660170047425588, 11.019092754691924, 36.64234042926972, 52.166359098954416, 131.44005705698146, 71.95922374674139, 109.27659166187777, 8.12538994839134, 24.207479589659716, 54.49864234514353, 8.394215456708961, 13.29099188085044, 87.25574374227118, 57.14132893996474, 29.99940454200873, 49.87084521061469, 16.446007253148593, 14.831139371859233, 76.65810104222564, 13.967785142642757, 21.115148364671953, 48.64735902389387, 6.62970676344421, 71.31736722240434, 12.468233825043775, 6.513137231132816, 61.69115976618602, 54.69081165622625, 32.963092682298196, 30.384666757174628, 27.89930970396219, 23.069261429950387, 5.9602205989615396, 34.83156137260096, 7.910088134013308, 17.359875859538437, 62.45234644551324, 16.48973797589025, 48.70373538989847, 11.266531131759114, 80.75471036608904, 71.77061776585, 9.857928625729334, 66.08535707538118, 49.848278498878656, 6.357694400854463, 87.06081404735941, 70.6627637770713, 32.03186611556872, 10.489117919340877, 153.46168328427834, 80.82039310971457, 14.22023405971005, 13.184687899079929, 73.29301168875114, 31.14870292714661, 19.363197210734917, 63.565408336594935, 11.864680640969805, 18.551292596396358, 7.036312773519994, 22.238418985422385, 18.506631252375602, 12.91634100634451, 70.94886765954244, 94.91807868698646, 7.381557464355967, 110.69531224531683, 106.34017512554217, 15.723031820776056, 29.409838891722575, 82.77878549863067, 22.54241075026417, 79.56421189125302, 122.9980516510281, 59.604921514839035, 5.714377217766062, 5.081274727558386, 23.629770790920563, 16.891040997691864, 11.51688395950215, 33.68587371120256, 17.8043541830509, 45.32292703450922, 20.681318307821012, 51.75598295007495, 63.377263088492185, 14.664060544467723, 17.436911749376428, 48.8104286884791, 31.467838187166503, 26.048113668271522, 26.273024039091503, 89.67299082683479, 27.515402816751955, 88.3901826370805, 26.29814146026113, 45.49511398434288, 44.90524951118547, 26.532478692975445, 55.2428955488178, 49.793403544275755, 50.87569320749922, 62.03314406369785, 78.97821755623272, 15.304778900270604, 62.201247421478456, 19.99745237598238, 38.52839964422239, 21.54917034216446, 103.51565854414545, 17.64156608598151, 145.4593856524358, 26.657779712734758, 19.519720315460518, 127.06404601446914, 75.36540517033468, 8.934913916819317, 9.90697976111298, 6.7092045550091095, 8.87224021382223, 38.18152125325919, 12.474097126212813, 10.976902718818442, 23.49248926413497, 44.047635883158094, 81.85640687175498, 5.726070875819977, 28.746029952540848, 30.876854950818323, 17.55746724957502, 62.144187615252555, 44.528774948730245, 14.481407248395923, 6.168933455879589, 12.815617864883375, 59.153207612370885, 13.691497122321188, 101.9284440798304, 22.347845293845147, 55.238201054080335, 28.481461562771454, 70.35314407269958, 13.092168474508384, 51.36845124333529, 41.721283775918565, 93.46683385249949, 20.800513050628126, 18.777776585600243, 5.90747328807188, 27.396136088327772, 18.647158631583117, 74.14988856292834, 14.95134793238071, 46.412787311015656, 8.532993988326437, 8.064661307903851, 29.40156838792065, 5.623419667498884, 59.53769571011693, 15.85178728808059, 23.097007536979238, 13.230820041497793, 11.67991406003155, 32.793827500641655, 23.680092722443156, 28.487469785378735, 34.22660152391056, 60.51935085050796, 64.2696727230055, 76.57524488209904, 9.873797679732139, 9.591154737291168, 95.49066297627802, 12.633821111497921, 46.6190565844058, 18.660298562787762, 31.123500360639106, 69.3944357263281, 45.91405222458597, 20.39776673615765, 20.38151942815262, 10.113691739717485, 84.38185358806999, 15.638609193024253, 8.122394366764263, 20.168226361937545, 27.822636818986375, 76.25945900789, 25.713568021890932, 18.59689613864327, 17.130666596408332, 26.688964911396912, 77.54017381426794, 9.264304194180871, 73.02081844012453, 59.18853490417307, 82.09590632966399, 48.51829706190693, 13.80238125734703, 43.93863612471406, 33.20971224524206, 19.89200826635866, 7.785847241239852, 48.60158334400493, 40.13477631564674, 6.855292216660307, 6.428238907641042, 83.47370864726389, 21.87070321216279, 110.16524089744877, 100.9008443695833, 33.15516884757169, 54.6271857424183, 14.282822637060761, 27.06758609068604, 5.583378603177801, 21.329669333758147, 23.80626062898029, 19.848780269837988, 29.171047211904916, 117.41566427668782, 80.26145077099804, 56.06546696714246, 8.219672710772471, 6.7163851879354945, 44.76630773488396, 30.87254847041948, 111.62958020851474, 60.05537100024193, 16.149333841698024, 60.05275947766313, 17.78581064773958, 6.367064299652842, 77.2479352983438, 16.642855823419023, 14.401007980535407, 23.182032338039065, 23.550792186263195, 26.186780864506463, 103.17134443433095, 34.59514545535535, 40.7988887733481, 133.09142995178235, 54.6440877361341, 46.96550380001123, 6.410374507676841, 5.24341684598228, 13.790303806704257, 16.756583799575893, 38.84844856884651, 54.653141264019645, 103.03755278251344, 29.421713828057616, 19.01905487800498, 11.276303739527185, 63.65566606747419, 35.09023204787191, 41.8347777387739, 25.522313355936795, 49.3888453456664, 7.543579133817603, 5.465633693719902, 5.271355276151485, 18.585010653915838, 41.987020887939195, 33.39544699610464, 5.098872764668381, 6.613008175899193, 33.50148939358008, 11.088989365801158, 11.927247712382789, 11.686608334973291, 25.161520452523277, 21.65281772680007, 7.773031395135736, 24.498405631440903, 58.52275196722579, 19.98488787771426, 6.468649639701314, 13.852649673103077, 39.530055844834955, 11.365135482619827, 133.76187876718976, 84.7384813813763, 6.554490104850635, 7.938191585915324, 5.809198587663047, 104.97410782690235, 13.19306113167598, 67.41521014959426, 132.89172246295493, 37.32232045216823, 106.96845039403092, 166.7920621873063, 7.796296673563342, 62.25816501533568, 10.272871640949646, 94.88455433207378, 10.906338731637183, 26.51085657130392, 12.378305105214896, 61.133276583946135, 73.06133011517392, 21.623828901995033, 130.73248346426675, 12.111854335398787, 64.58310030372296, 7.202761413681687, 6.66470913698336, 7.300536287094914, 13.66260938967418, 31.36215478048632, 69.08257636658377, 29.91356266487371, 147.20001365691377, 43.287191073634915, 35.98030459908736, 51.818181167478414, 23.306662844659694, 16.959743958596363, 20.320186323756786, 82.49483306967863, 23.87802323874675, 33.21173986036255, 84.65118678111139, 10.161947438066907, 5.9943931676093705, 18.94886324706177, 63.275757943643846, 113.37240814371873, 86.82889063994705, 17.682905557180604, 5.715365398316607, 51.9313193529557, 24.125025749365097, 67.10846952277934, 86.33003899555197, 11.124506932028083, 7.705518793252499, 11.273170528706327, 40.778511012024815, 92.25678238781433, 22.567473511396813, 52.80326704158034, 14.488835629977112, 66.7529617386063, 49.069312322289875, 56.8709807779707, 5.246241624991588, 19.868408961669093, 175.26049278450168, 20.236829332667504, 20.580859377899355, 18.921560029395177, 56.64419763233434, 98.26725037637073, 11.631693985190036, 11.338132534469736, 49.851144234986236, 21.342413660703365, 53.72032388418809, 19.593143817440968, 11.741479634264495, 106.513350608544, 5.305050581674112, 51.21543539192766, 30.365448796487218, 62.6179015062853, 82.33111529571401, 38.38677243626686, 19.280299440029765, 62.04096657765081, 10.699792804440605, 58.03534618647714, 25.20306004959948, 79.94184810960904, 124.90520617156518, 9.218839372630129, 47.05321074380154, 49.68121313762114, 46.17582538454992, 17.404379228090576, 25.52761709067381, 6.718178458050407, 55.74669613189269, 7.154509276863198, 11.418980654067878, 72.39913812446233, 25.71048634471073, 12.621917317343861, 5.225460690407179, 29.16689130934563, 52.82118755875637, 64.37664063591447, 60.43871007587107, 13.252682827118646, 83.25972817088461, 104.5110256013063, 84.77956184298004, 18.47399121285312, 23.831523349094752, 10.633346382443696, 17.90510779735961, 50.32705636572011, 6.386111902768931, 43.140198657118546, 22.913514130900406, 5.474398605351051, 25.20467608790601, 37.21375192816757, 37.60414558583007, 8.719542550059671, 150.90357192146126, 31.65445329142686, 8.031721305098555, 13.016472551202959, 37.49072403505779, 81.54250882004574, 64.140469946995, 39.446260516502456, 79.92310192098643, 14.184403899686975, 125.00160504645694, 60.27709883940308, 14.766355922354824, 7.566035487226929, 11.90548248748581, 6.20378697940843, 26.817275631840847, 22.769311485704314, 51.53068194112272, 32.716439998202574, 13.763438673215541, 6.970011162992154, 28.072632619209237, 12.682655047810279, 15.946165161397058, 66.28638209950796, 25.33004430235787, 58.43077602174903, 23.68868289961794, 10.332587592674747, 5.336546280758336, 31.041908306121343, 86.92095281556595, 57.38300761848095, 30.09479958875559, 133.33951786138735, 35.010125498291586, 65.72176873422268, 52.59428477568801, 15.1824701921035, 26.04619362604328, 157.36656400576214, 87.87778969766059, 74.81181180318475, 28.94930427900464, 5.535408529767064, 13.074668224043917, 6.2539594389111315, 9.1585429445389, 66.5214701444361, 25.4385661763878, 61.47616848217645, 12.2334460682809, 37.95016324449648, 70.91680563317215, 26.680775348307794, 30.700660686446525, 28.333233915575754, 22.568314438068622, 5.497304982047925, 22.044268297984786, 60.21149368833324, 6.325423978666722, 158.06260499281035, 56.245033666554306, 5.465762899336625, 27.049859792403083, 54.16250030105062, 24.866674301939153, 7.3423362838986455, 9.002948551042019, 17.72683426307021, 10.049375478318348, 8.084229999503211, 56.06073486591909, 56.174026562532624, 42.193191136692036, 6.685481002078481, 17.806589890005085, 19.37101629501525, 22.101477452675635, 18.109836629526605, 94.71065431426054, 41.26115423990972, 231.2090986539324, 7.606534083473084, 9.718960799781552, 9.686370344131078, 53.22109012810866, 14.871360304667427, 36.4910188154177, 40.67924926450571, 14.843721351013137, 32.70807733900208, 90.70543992461603, 77.86839180027894, 14.762348830490495, 94.54933088113195, 9.143855656068174, 30.957854703081367, 45.91616662957094, 53.07797685442267, 5.422319615225528, 16.638852230915695, 5.7225621590299305, 81.55801831222949, 69.34543000426277, 116.174233207947, 5.052943386685845, 6.20888682727358, 57.07011504416968, 61.30565866995507, 27.338937424105687, 7.531127503764302, 58.26624235704098, 6.56106597321113, 22.71922041011763, 18.862919954016824, 5.217634064397359, 17.04895311596002, 13.19513978403063, 61.39846085789381, 5.726076903558251, 65.71865766160028, 12.793227880158955, 16.37063839248726, 58.10246394639079, 10.459158889853171, 41.577994702543215, 61.210382709045156, 81.16742583011197, 59.327075090281234, 11.98169076938157, 45.577654238132844, 53.06708763869012, 45.96337653185522, 48.113244558365565, 23.574501193124455, 28.12615448274729, 10.954360261608782, 65.07561663177394, 7.863496720619498, 12.580716471574828, 36.3830052370592, 64.13341800277311, 29.623305339205345, 14.192461473662146, 8.443234665134106, 8.15475500420226, 5.423595598713941, 7.964121563174872, 18.01102258969498, 23.920927179573255, 15.145847000260014, 50.91316388665784, 11.10362465561797, 10.054041913523104, 59.318671319467505, 61.497372171578625, 6.866450655096844, ...])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xmean,count)
size_of_interactions_average(Etot, Ntot,size_interaction)
#energy_spectra_doublebreakdown(Etot, T_first, file_energy_spectra_doublebreakdown)