script for convert .t3pa files to .t3pa_cls
.t3pa file example:
Index $\quad$ Matrix $\quad$ Index $\quad$ ToA $\quad$ ToT $\quad$ FToA $\quad$ Overflow
0 $\quad$ 4574 $\quad$ 832 $\quad$ 29 $\quad$ 6 $\quad$ 0
1 $\quad$ 4831 $\quad$ 832 $\quad$ 35 $\quad$ 7 $\quad$ 0
2 $\quad$ 4575 $\quad$ 832 $\quad$ 100 $\quad$ 8 $\quad$ 0
3 $\quad$ 31031 $\quad$ 1745 $\quad$ 22 $\quad$ 11 $\quad$ 0
.
.
.
.t3pa_cls file example:
% Index $\quad$ Matrix Index $\quad$ [ RowNo, ClmNo ] $\quad$ ToA $\quad$ FToA $\quad$ ( ToA_in_ns ) $\quad$ ToT ( ToT_in_keV ) $\quad$ Overflow
# 1, $\quad$ Nunmasked = 3, $\quad$ Nmasked = 0, $\quad$ Ntot = 3 # Tfirst = 2.0787500000000000e+04 ns, $\quad$ Tlast = 2.0790625000000000e+04 ns, $\quad$ dT = 3.125000 ns, $\quad$ Etot = 64.428148 keV
2 $\quad$ 4575 $\quad$ [ 17, 223 ] $\quad$ 832 $\quad$ 8 $\quad$ ( 2.0787500000000000e+04 ns ) $\quad$ 100 $\quad$ ( 37.867914 keV ) $\quad$ 0
1 $\quad$ 4831 $\quad$ [ 18, 223 ] $\quad$ 832 $\quad$ 7 $\quad$ ( 2.0789062500000000e+04 ns ) $\quad$ 35 $\quad$ ( 14.733453 keV ) $\quad$ 0
0 $\quad$ 4574 $\quad$ [ 17, 222 ] $\quad$ 832 $\quad$ 6 $\quad$ ( 2.0790625000000000e+04 ns ) $\quad$ 29 $\quad$ ( 11.826781 keV ) $\quad$ 0
# 2, $\quad$ Nunmasked = 3, $\quad$ Nmasked = 0, $\quad$ Ntot = 3 # Tfirst = 4.3601562500000000e+04 ns, $\quad$ Tlast = 4.3607812500000000e+04 ns, $\quad$ dT = 6.250000 ns, $\quad$ Etot = 63.577435 keV
5 $\quad$ 30775 $\quad$ [ 120, 55 ] $\quad$ 1745 $\quad$ 15 $\quad$ ( 4.3601562500000000e+04 ns ) $\quad$ 99 $\quad$ ( 37.617059 keV ) $\quad$ 0
4 $\quad$ 30776 $\quad$ [ 120, 56 ] $\quad$ 1745 $\quad$ 13 $\quad$ ( 4.3604687500000000e+04 ns ) $\quad$ 44 $\quad$ ( 14.715446 keV ) $\quad$ 0
3 $\quad$ 31031 $\quad$ [ 121, 55 ] $\quad$ 1745 $\quad$ 11 $\quad$ ( 4.3607812500000000e+04 ns ) $\quad$2 2 $\quad$ ( 11.244929 keV ) $\quad$ 0
.
.
.
import numpy as np
import math
#import pandas as pd
import matplotlib.pyplot as plt
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"
#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 = 36529 #test discharge for which the notebook will definitely work
shot_no = 44419
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
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))
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 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=[]
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))
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')
i=i+1
t3pa_cls.close
return T_first, E_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)
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.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.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)
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].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
#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'
#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'
#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)
#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 = 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);
([3470703.125, 3480925.0, 3598635.9375, 3653518.75, 3739956.25, 3801859.375, 3801917.1875, 3816970.3125, 3831207.8125, 3905751.5625, 3939035.9375, 3939779.6875, 3952146.875, 3995578.125, 4061498.4375, 4099175.0, 4335209.375, 4352259.375, 4366034.375, 4380714.0625, 4403509.375, 4422071.875, 4440512.5, 4440539.0625, 4444476.5625, 4453784.375, 4465164.0625, 4491889.0625, 4491940.625, 4498448.4375, 4502106.25, 4511367.1875, 4515553.125, 4519835.9375, 4529248.4375, 4538576.5625, 4551434.375, 4560790.625, 4561012.5, 4561021.875, 4566389.0625, 4567720.3125, 4598221.875, 4623303.125, 4639548.4375, 4657598.4375, 4722226.5625, 4728223.4375, 4735293.75, 4750329.6875, 4753018.75, 4816007.8125, 4831603.125, 4831620.3125, 4858807.8125, 4904787.5, 4938901.5625, 5004753.125, 5042281.25, 5838914.0625, 5838939.0625, 5868325.0, 5878576.5625, 5884629.6875, 5909895.3125, 5910268.75, 5913950.0, 5920875.0, 5922079.6875, 5923128.125, 5924895.3125, 5926420.3125, 5931232.8125, 5931289.0625, 5931648.4375, 5937810.9375, 5950868.75, 5952243.75, 5956689.0625, 5975262.5, 6010259.375, 6034350.0, 6034367.1875, 6035717.1875, 6045203.125, 6045373.4375, 6107350.0, 6107462.5, 6116495.3125, 6117214.0625, 6118854.6875, 6119615.625, 6122517.1875, 6122528.125, 6123198.4375, 6137473.4375, 6140943.75, 6153803.125, 6154179.6875, 6164606.25, 6172214.0625, 6178732.8125, 6179531.25, 6181407.8125, 6205753.125, 6211875.0, 6213450.0, 6225665.625, 6231765.625, 6236398.4375, 6236400.0, 6236404.6875, 6260998.4375, 6270893.75, 6276895.3125, 6284298.4375, 6285343.75, 6290525.0, 6292803.125, 6294256.25, 6294520.3125, 6295495.3125, 6295501.5625, 6295828.125, 6296404.6875, 6296918.75, 6297087.5, 6298154.6875, 6298265.625, 6301901.5625, 6310104.6875, 6311875.0, 6328570.3125, 6328709.375, 6329470.3125, 6337890.625, 6338879.6875, 6348040.625, 6348187.5, 6354000.0, 6359185.9375, 6362120.3125, 6365779.6875, 6365871.875, 6369507.8125, 6373795.3125, 6374860.9375, 6374959.375, 6375948.4375, 6381031.25, 6384264.0625, 6384392.1875, 6385601.5625, 6386307.8125, 6386879.6875, 6390520.3125, 6395662.5, 6396964.0625, 6399181.25, 6401507.8125, 6407946.875, 6408001.5625, 6410643.75, 6410648.4375, 6411650.0, 6412393.75, 6412451.5625, 6416378.125, 6416635.9375, 6418126.5625, 6419359.375, 6419471.875, 6422156.25, 6426595.3125, 6431257.8125, 6431275.0, 6443267.1875, 6444823.4375, 6445287.5, 6450365.625, 6453176.5625, 6453195.3125, 6453445.3125, 6453525.0, 6453953.125, 6454562.5, 6455629.6875, 6456620.3125, 6462695.3125, 6462704.6875, 6464318.75, 6465598.4375, 6469204.6875, 6473532.8125, 6475532.8125, 6477868.75, 6477882.8125, 6479181.25, 6481104.6875, 6481543.75, 6482462.5, 6483679.6875, 6485120.3125, 6490718.75, 6490878.125, 6493290.625, 6495804.6875, 6503553.125, 6510087.5, 6512015.625, 6512403.125, 6513850.0, 6514632.8125, 6516128.125, 6518473.4375, 6521920.3125, 6522589.0625, 6524535.9375, 6528795.3125, 6529454.6875, 6531757.8125, 6532001.5625, 6532228.125, 6533046.875, 6533890.625, 6533993.75, 6534092.1875, 6534104.6875, 6534398.4375, 6534485.9375, 6534540.625, 6534870.3125, 6535714.0625, 6536657.8125, 6537370.3125, 6538042.1875, 6538471.875, 6538646.875, 6538751.5625, 6538987.5, 6541532.8125, 6541592.1875, 6541934.375, 6542043.75, 6543900.0, 6543946.875, 6548914.0625, 6550000.0, 6555868.75, 6556376.5625, 6556914.0625, 6559601.5625, 6559757.8125, 6562078.125, 6562137.5, 6562659.375, 6563259.375, 6563328.125, 6563348.4375, 6563589.0625, 6563610.9375, 6564532.8125, 6565235.9375, 6570531.25, 6571701.5625, 6572684.375, 6573951.5625, 6573985.9375, 6575054.6875, 6575329.6875, 6575429.6875, 6576720.3125, 6577021.875, 6578181.25, 6578235.9375, 6579962.5, 6580465.625, 6580521.875, 6580626.5625, 6582337.5, 6583095.3125, 6586312.5, 6587685.9375, 6587742.1875, 6590193.75, 6590210.9375, 6590318.75, 6591201.5625, 6593070.3125, 6594659.375, 6594965.625, 6597537.5, 6599832.8125, 6600621.875, 6600692.1875, 6601557.8125, 6601810.9375, 6602834.375, 6605082.8125, 6605095.3125, 6605914.0625, 6607081.25, 6607164.0625, 6607885.9375, 6609381.25, 6609635.9375, 6617843.75, 6622465.625, 6626889.0625, 6628828.125, 6631379.6875, 6633635.9375, 6634189.0625, 6637800.0, 6638217.1875, 6639678.125, 6648617.1875, 6648714.0625, 6649259.375, 6649587.5, 6654384.375, 6656885.9375, 6675267.1875, 6676045.3125, 6679459.375, 6688209.375, 6688637.5, 6688857.8125, 6692514.0625, 6694685.9375, 6700328.125, 6700721.875, 6701975.0, 6703598.4375, 6707439.0625, 6716284.375, 6717662.5, 6722821.875, 6726289.0625, 6727328.125, 6729115.625, 6729201.5625, 6733485.9375, 6733845.3125, 6734176.5625, 6734270.3125, 6735014.0625, 6737450.0, 6737739.0625, 6737745.3125, 6738331.25, 6739204.6875, 6740117.1875, 6743189.0625, 6744381.25, 6746832.8125, 6748928.125, 6758453.125, 6765581.25, 6766345.3125, 6767529.6875, 6769968.75, 6773750.0, 6775525.0, 6778418.75, 6785853.125, 6788176.5625, 6789664.0625, 6791737.5, 6797053.125, 6798537.5, 6802812.5, 6806085.9375, 6807973.4375, 6808314.0625, 6809484.375, 6810926.5625, 6810946.875, 6817709.375, 6820435.9375, 6821660.9375, 6824609.375, 6824776.5625, 6827606.25, 6827667.1875, 6830176.5625, 6830321.875, 6832028.125, 6832303.125, 6832750.0, 6832782.8125, 6833887.5, 6834196.875, 6834754.6875, 6835962.5, 6837501.5625, 6839189.0625, 6839371.875, 6840176.5625, 6841854.6875, 6841890.625, 6841989.0625, 6842385.9375, 6843281.25, 6843492.1875, 6844131.25, 6844210.9375, 6844256.25, 6845025.0, 6845179.6875, 6846879.6875, 6853315.625, 6877454.6875, 6877467.1875, 6879728.125, 6880060.9375, 6883640.625, 6885059.375, 6885073.4375, 6891931.25, 6892937.5, 6893018.75, 6895690.625, 6897728.125, 6918784.375, 6918787.5, 6920029.6875, 6921639.0625, 6922851.5625, 6925350.0, 6925621.875, 6927460.9375, 6932770.3125, 6934234.375, 6940609.375, 6941910.9375, 6941950.0, 6942523.4375, 6943695.3125, 6945593.75, 6945820.3125, 6947032.8125, 6947296.875, 6952042.1875, 6960453.125, 6977898.4375, 6980271.875, 6981100.0, 6981117.1875, 6981534.375, 6984995.3125, 6986426.5625, 6986484.375, 6986551.5625, 6987107.8125, 6987556.25, 6987846.875, 6987859.375, 6988218.75, 6988842.1875, 6988959.375, 6988985.9375, 6989051.5625, 6989057.8125, 6989379.6875, 6989409.375, 6990923.4375, 6991404.6875, 6991453.125, 6991573.4375, 6991942.1875, 6992051.5625, 6993012.5, 6993681.25, 6993881.25, 6994371.875, 6994445.3125, 6994925.0, 6995023.4375, 6995107.8125, 6995309.375, 6995473.4375, 6996628.125, 6996981.25, 6997337.5, 6997748.4375, 6997798.4375, 6997889.0625, 6998454.6875, 6998489.0625, 6998540.625, 6998542.1875, 6999332.8125, 6999451.5625, 6999521.875, 7001189.0625, 7001204.6875, 7001207.8125, 7001390.625, 7003003.125, 7003445.3125, 7005314.0625, 7005335.9375, 7009192.1875, 7016895.3125, 7020592.1875, 7022682.8125, 7032073.4375, 7032923.4375, 7035607.8125, 7040003.125, 7042409.375, 7043978.125, 7045521.875, 7045550.0, 7045654.6875, 7047379.6875, 7051001.5625, 7052907.8125, 7063760.9375, 7080603.125, 7088414.0625, 7097507.8125, 7111150.0, 7118500.0, 7122682.8125, 7123173.4375, 7123182.8125, 7129781.25, 7130006.25, 7130442.1875, 7150948.4375, 7151957.8125, 7152945.3125, 7153857.8125, 7157454.6875, 7158378.125, 7158442.1875, 7160437.5, 7167154.6875, 7167593.75, 7168675.0, 7169735.9375, 7169871.875, 7170834.375, 7176632.8125, 7180800.0, 7180803.125, 7186521.875, 7186787.5, 7190051.5625, 7191351.5625, 7192045.3125, 7192112.5, 7196346.875, 7199087.5, 7199139.0625, 7199753.125, 7200848.4375, 7201506.25, 7201720.3125, 7202951.5625, 7202957.8125, 7202971.875, 7202975.0, 7203603.125, 7203700.0, 7203759.375, 7204253.125, 7206335.9375, 7206928.125, 7209689.0625, 7222415.625, 7224040.625, 7228606.25, 7229360.9375, 7230101.5625, 7232987.5, 7235026.5625, 7236795.3125, 7255696.875, 7257245.3125, 7261845.3125, 7262132.8125, 7264621.875, 7267503.125, 7268381.25, 7269864.0625, 7271757.8125, 7278056.25, 7282176.5625, 7286479.6875, 7301062.5, 7306542.1875, 7308003.125, 7312517.1875, 7313971.875, 7318020.3125, 7327485.9375, 7332106.25, 7344387.5, 7353012.5, 7355571.875, 7370806.25, 7374682.8125, 7378696.875, 7378742.1875, 7380267.1875, 7384067.1875, 7384193.75, 7387600.0, 7391070.3125, 7394920.3125, 7396345.3125, 7397000.0, 7398423.4375, 7399346.875, 7399812.5, 7399842.1875, 7399860.9375, 7400406.25, 7400467.1875, 7401298.4375, 7401298.4375, 7401443.75, 7401820.3125, 7401942.1875, 7402331.25, 7402517.1875, 7402585.9375, 7403037.5, 7403042.1875, 7403589.0625, 7403900.0, 7404025.0, 7404831.25, 7405089.0625, 7405331.25, 7405571.875, 7406167.1875, 7406365.625, 7406421.875, 7406443.75, 7406571.875, 7406625.0, 7406628.125, 7406768.75, 7406898.4375, 7406985.9375, 7407554.6875, 7407560.9375, 7409626.5625, 7409710.9375, 7410323.4375, 7410696.875, 7410720.3125, 7410920.3125, 7411446.875, 7412246.875, 7412610.9375, 7413185.9375, 7413284.375, 7413585.9375, 7415032.8125, 7415071.875, 7418085.9375, 7418142.1875, 7418146.875, 7418939.0625, 7420471.875, 7420496.875, 7421318.75, 7421784.375, 7422079.6875, 7422989.0625, 7423114.0625, 7423535.9375, 7423620.3125, 7423739.0625, 7424296.875, 7424368.75, 7424629.6875, 7424981.25, 7425031.25, 7425132.8125, 7425151.5625, 7425226.5625, 7425331.25, 7425989.0625, 7427351.5625, 7428912.5, 7429190.625, 7430017.1875, 7430881.25, 7430885.9375, 7430940.625, 7431781.25, 7431784.375, 7431870.3125, 7432498.4375, 7433459.375, 7433464.0625, 7437148.4375, 7438045.3125, 7438800.0, 7439235.9375, 7439407.8125, 7442600.0, 7444151.5625, 7451131.25, 7452737.5, 7453160.9375, 7453285.9375, 7453464.0625, 7454062.5, 7454131.25, 7455009.375, 7455282.8125, 7455856.25, 7456696.875, 7457203.125, 7459103.125, 7459659.375, 7459793.75, 7463287.5, 7464278.125, 7464387.5, 7466350.0, 7469775.0, 7469829.6875, 7470676.5625, 7471781.25, 7479200.0, 7479704.6875, 7481104.6875, 7482906.25, 7488018.75, 7490910.9375, 7492776.5625, 7504732.8125, 7506354.6875, 7506360.9375, 7507032.8125, 7509010.9375, 7517221.875, 7520279.6875, 7539343.75, 7542398.4375, 7542471.875, 7542475.0, 7542571.875, 7542679.6875, 7542687.5, 7542807.8125, 7543168.75, 7543812.5, 7543864.0625, 7545325.0, 7545520.3125, 7545715.625, 7546445.3125, 7546448.4375, 7546454.6875, 7547245.3125, 7547459.375, 7547767.1875, 7548000.0, 7548103.125, 7548165.625, 7548456.25, 7548637.5, 7548821.875, 7548914.0625, 7549021.875, 7549201.5625, 7549326.5625, 7549345.3125, 7549782.8125, 7550331.25, 7550503.125, 7550539.0625, 7551343.75, 7551512.5, 7552079.6875, 7552206.25, 7552321.875, 7552971.875, 7553295.3125, 7554565.625, 7554662.5, 7554762.5, 7554907.8125, 7554945.3125, 7555431.25, 7556512.5, 7558326.5625, 7559665.625, 7570739.0625, 7571868.75, 7575484.375, 7575648.4375, 7575720.3125, 7575804.6875, 7577882.8125, 7577887.5, 7579265.625, 7579795.3125, 7579984.375, 7580779.6875, 7581193.75, 7581590.625, 7581696.875, 7582870.3125, 7583057.8125, 7584373.4375, 7584904.6875, 7587115.625, 7594201.5625, 7598643.75, 7602770.3125, 7603598.4375, 7605089.0625, 7605134.375, 7605914.0625, 7606571.875, 7607732.8125, 7607981.25, 7616067.1875, 7623285.9375, 7625106.25, 7626532.8125, 7628393.75, 7630475.0, 7631223.4375, 7631621.875, 7649718.75, 7658464.0625, 7664843.75, 7677440.625, 7679860.9375, 7682253.125, 7685754.6875, 7696390.625, 7697290.625, 7700509.375, 7700785.9375, 7700865.625, 7706804.6875, 7712309.375, 7713662.5, 7719706.25, 7720771.875, 7757917.1875, 7759796.875, 7760037.5, 7760629.6875, 7762906.25, 7763150.0, 7763242.1875, 7763879.6875, 7764178.125, 7764412.5, 7764415.625, 7764895.3125, 7766042.1875, 7766298.4375, 7767123.4375, 7767359.375, 7767546.875, 7767732.8125, 7767743.75, 7767879.6875, 7768082.8125, 7768085.9375, 7768504.6875, 7768779.6875, 7769778.125, 7769887.5, 7769901.5625, 7769945.3125, 7769968.75, 7770259.375, 7770329.6875, 7770342.1875, 7770381.25, 7770573.4375, 7770579.6875, 7770607.8125, 7770623.4375, 7770917.1875, 7770925.0, 7771493.75, 7771523.4375, 7771650.0, 7771703.125, 7772079.6875, 7776275.0, 7781475.0, 7781792.1875, 7783450.0, 7787321.875, 7787773.4375, 7788987.5, 7789923.4375, 7790475.0, 7792484.375, 7792971.875, 7799373.4375, 7801475.0, 7802045.3125, 7802546.875, 7802665.625, 7803865.625, 7803998.4375, 7805984.375, 7808557.8125, 7813931.25, 7828271.875, 7862162.5, 7863454.6875, 7864765.625, 7865290.625, 7866617.1875, 7866756.25, 7866757.8125, 7868154.6875, 7868167.1875, 7868465.625, 7868681.25, 7869381.25, 7869670.3125, 7869812.5, 7869987.5, 7870084.375, 7870825.0, 7870915.625, 7870923.4375, 7871015.625, 7872685.9375, 7881601.5625, 7883218.75, 7884246.875, 7892050.0, 7893835.9375, 7893890.625, 7894515.625, 7894953.125, 7895953.125, 7899523.4375, 7899595.3125, 7900356.25, 7902392.1875, 7903342.1875, 7918568.75, 7918996.875, 7919084.375, 7919198.4375, 7929390.625, 7930764.0625, 7939103.125, 7939123.4375, 7941054.6875, 7942542.1875, 7942778.125, 7944251.5625, 7944532.8125, 7945579.6875, 7946056.25, 7947525.0, 7947751.5625, 7952901.5625, 7952968.75, 7953270.3125, 7954014.0625, 7959406.25, 7963553.125, 7975454.6875, 7975643.75, 7990943.75, 7991581.25, 7992701.5625, 7992789.0625, 7993678.125, 7993996.875, 7995153.125, 7995551.5625, 7997607.8125, 7997617.1875, 7997660.9375, 7998934.375, 8000775.0, 8002462.5, 8006862.5, 8033621.875, 8034043.75, 8039560.9375, 8039589.0625, 8042185.9375, 8048268.75, 8048740.625, ...], [45.39907160677232, 58.43472597827107, 54.13986925637357, 89.12066462801411, 42.31772412446346, 44.334811280756355, 5.039176232203347, 5.998134463936753, 45.86931391332472, 84.75723619386078, 17.776354316266083, 13.483571108427403, 26.221404666351518, 122.39943040280984, 6.617917723442894, 6.467919871230439, 5.631474168321943, 39.62706428287824, 15.071489453073063, 84.36847826439555, 22.102515067820722, 6.216946875875842, 9.868371161867207, 10.322389149980804, 62.79481252107678, 93.3164113725934, 17.54321205886551, 12.152756737713236, 5.093415425987701, 6.376425591099246, 8.405526394965008, 74.97857949686866, 44.42105796148373, 29.4333807339959, 7.44665754383101, 29.082651627926793, 36.500974862091454, 12.555573518695828, 10.942724627448971, 7.588849913883986, 19.265618039863664, 9.809752109315776, 10.8106488535208, 80.32804657115567, 39.30421997334413, 16.6046771494948, 10.846970809067779, 16.143030242617993, 9.621170738191072, 6.45128784867267, 29.07133178912585, 10.544364280910047, 10.960144356161543, 59.14702833810563, 9.847484428910874, 41.4362725105399, 19.777389623817008, 26.1335112794383, 6.59164146767698, 13.939314326181142, 95.82585115343834, 23.38383925029865, 48.384512544701074, 78.15949706271037, 48.56458390501364, 5.475348409468644, 20.003753260607663, 26.05833028013529, 22.2631641647619, 59.94786198705246, 19.383723204735603, 51.56502116682048, 37.29317231216768, 6.1591816033042175, 181.55652163621295, 66.91443961230476, 31.526858619254003, 75.95028961562377, 65.21755936792128, 16.49099155461073, 55.00423227520606, 58.93448126855243, 52.57303396798196, 5.217553170091159, 6.297038156735923, 55.762434798435486, 122.59244627234175, 5.567883196850918, 36.64121541488077, 9.123455006903669, 6.155750364387428, 115.11362335820226, 14.155945289610052, 13.896946396998004, 109.52134482320002, 6.521761124976892, 28.894246811617883, 17.701096686631846, 6.841739455429855, 115.98313440439134, 87.99904748472544, 104.12694747318481, 38.645243162261444, 7.5688279260069775, 18.211069104169813, 122.5099663949049, 56.258755716947455, 42.31201243482917, 85.38430680996511, 28.627017323532996, 78.52249036185356, 22.257552213644875, 8.337888736273742, 10.884909902192174, 73.42403188280775, 18.94452891252995, 64.69291889600817, 53.75845680785268, 200.4307381448822, 15.403083707772161, 42.345565772351115, 10.792066111075812, 48.952598545524296, 82.04082184310545, 17.459621910125318, 22.578644341523503, 32.518121959953405, 18.596584589046813, 23.562605590381438, 66.2884170327356, 7.150933068576006, 74.30785344515007, 10.707748350875635, 59.1317663370279, 8.145896436792253, 39.19595319609624, 22.69215664318009, 6.756174477758828, 86.82839859316316, 23.55519733864917, 24.66577642177038, 8.245412884681722, 41.888715122311275, 5.177433677294963, 19.40547904648074, 5.395565708163598, 7.059264146471218, 5.156112256451142, 24.265766943295723, 25.07708309080054, 45.60033010662322, 17.213412925676487, 88.4263384452878, 62.19277057868756, 18.004235835022033, 9.285092982644406, 37.843745776663624, 21.201363741903585, 6.189170000128272, 17.632239923206257, 12.158797066862332, 5.387424293221695, 19.956927791612976, 31.544750338157037, 19.833414538623572, 65.70790924322222, 5.9092602237596985, 9.053405579304599, 88.24950473840443, 7.046788741057286, 52.259797094185, 32.91496672789286, 5.5568869593321235, 31.668638883872763, 14.004768803728018, 12.447933234196881, 91.88842608447418, 50.395568154099, 53.13890747925048, 12.856285612819297, 89.23664603720526, 58.09053308802784, 21.94958266811906, 17.63833934598756, 65.29770703729992, 69.71905195484634, 30.90831948911661, 15.933674677695823, 22.296600553194374, 60.34363931774433, 10.461005516721935, 5.0454657469302715, 50.0406825842092, 24.414766266230497, 9.893918438928216, 15.744856432955881, 10.643970707191222, 21.10846471513952, 136.64393338770643, 6.412187549294856, 22.571910458958776, 81.54459522258487, 16.673840323265726, 6.0985001814627715, 30.4960606464605, 163.44339035452413, 18.918455278445634, 44.91867856895598, 27.909485983714127, 29.521661742604593, 23.275619798221083, 57.28065002286976, 29.385702597457005, 7.4102310088436685, 102.6445590302597, 12.474542145873293, 21.56091987854063, 118.42552625068336, 153.20306438310783, 42.855443346989965, 96.45129324467587, 48.67620800629379, 53.447686966167716, 75.11213513887965, 45.31277737228109, 71.88347432463344, 23.646482817622456, 26.69004483286118, 160.01336347064196, 60.671913085868574, 6.388936207090076, 86.339830885314, 68.67242066656685, 20.332071743757655, 192.90851180289417, 13.50086638863289, 18.73622133159905, 61.92237559855329, 25.657698742525973, 57.16394346391178, 7.231209516791493, 94.66997077007252, 65.28848137767659, 115.5350612809463, 105.97551059442645, 85.70309810791288, 76.99501375011621, 21.532455954106926, 186.96168205187627, 84.84958542742243, 19.59833539313485, 130.41424353371184, 20.606194671347275, 44.53396692833216, 5.147449868012342, 27.402908760491336, 41.23654180808661, 6.320732704602066, 5.066039091525922, 27.528909248391273, 11.940316350267715, 5.64966048804249, 15.34483684382827, 68.79196299466794, 49.4414407200757, 25.297074013286498, 25.80380395072883, 49.925189759728056, 203.41652578947088, 41.37034528561626, 19.413489063201066, 22.740068262720936, 69.91825689622371, 6.491557075161738, 5.245993923447846, 40.1234117878206, 61.13955880638159, 105.04512716729712, 24.02470375118437, 95.43278300318732, 38.86067762456453, 52.8405476598883, 5.597005252311804, 5.1062186712719075, 120.85197042249324, 37.04775804224042, 148.4827380106409, 6.054074193955454, 65.56498508992377, 60.06850119040066, 66.09593712981703, 6.289839946194938, 97.59979343563244, 67.75257948420965, 12.55574537977661, 19.579557538956394, 42.243455950252496, 93.94876955573942, 10.214295549138312, 86.26770136285887, 11.239175216115399, 93.89453134919016, 48.632434070424566, 56.31239966064271, 16.974253018025315, 13.284059282926744, 55.94175496813692, 45.328864203068946, 33.46019774074443, 33.67656916228863, 166.85197878594138, 35.078030617766686, 20.302544005692127, 13.688339736074413, 15.04452041373881, 173.10406082439712, 119.04374525554634, 5.536420164809132, 28.910807613527968, 11.798533416214562, 71.23856367948319, 36.56810007910176, 5.183970751278581, 6.271037943667835, 16.58279971874704, 99.16234110249987, 5.29110979192944, 63.82699222193669, 165.34392507816844, 56.558922170957445, 8.417177091023483, 60.70816415202323, 6.358755406015505, 10.224438558589181, 126.75165885212118, 33.45551341054375, 87.51603664614346, 41.34328960950465, 38.45863954065044, 30.71793378594299, 52.573576853859514, 12.405968439363845, 9.936772881970768, 11.516754586909748, 53.208626157814926, 5.1254533998877525, 33.38268387540789, 56.57641282180743, 25.377644201192204, 31.566780382558477, 54.60241753769744, 164.83431835179906, 6.041046825041573, 18.417918709146356, 87.33145765011676, 28.277067543023254, 50.450172829026606, 6.138004097805336, 10.930642203226215, 20.977485891403944, 21.271989860957593, 25.08747440137797, 46.577723697459646, 15.287785020210238, 15.68805122087268, 92.15269797873349, 71.33130413015576, 18.96887726153514, 91.48311912760275, 36.46875108484006, 38.0887057852149, 21.824647857830204, 30.153747873777245, 44.71174270140521, 100.81145934145079, 102.36991611454101, 11.274615133291446, 5.556098077682061, 11.931118094582532, 6.526406230879605, 81.51127708683586, 63.55272371555312, 36.29085351670314, 140.76784046351992, 5.469384465884944, 29.803411402502242, 43.557315029470246, 45.30427071432581, 8.089787249377078, 73.45350093798532, 31.186414730630744, 76.10207636595962, 51.928865988721824, 36.146434187043766, 72.97701412916159, 11.08585042027072, 40.21400977253449, 14.2150281952447, 14.058412170920274, 11.115480864398485, 16.201867928201565, 69.17968872420569, 11.281210534829352, 67.65173082092512, 49.242126947805446, 43.83393473092084, 7.7995722527754, 152.52156616805956, 19.984556308480105, 77.18302296339232, 8.048554268843452, 23.02039655196422, 128.0673188399014, 50.196296468777554, 9.34252793461419, 47.299820719892935, 18.782294796546264, 133.53990843514185, 67.38985519011622, 58.03749151043067, 21.14784692288888, 50.73945281382011, 13.460610761028006, 49.68708529325806, 35.935099630540044, 29.221446354052762, 34.65648861643882, 16.469711879910875, 97.20530485075014, 84.1989207251092, 25.445292421066412, 74.19879793961121, 31.65260505276899, 118.60926307237679, 6.201124605319172, 55.74620815433049, 21.715132998828715, 6.208759179210996, 11.716829216942017, 9.12402337275018, 22.578604468202663, 36.38174943106439, 59.61004846525074, 13.47185660808665, 87.03466393497298, 11.6279752627118, 5.350525528549183, 7.239744481714719, 51.47016199737295, 10.47602973334048, 17.09946966572913, 123.54978959282262, 35.51582230486843, 31.024608105432446, 8.78835583754656, 33.85602790624039, 51.527600122273405, 6.439514952237363, 15.531036282208104, 32.40325296215359, 20.829721709435912, 68.92696691691638, 5.171471823328056, 22.822942443731623, 61.68168187061277, 12.159356070165803, 60.373039087197085, 108.03092352171295, 146.72266425352174, 58.78546962661646, 42.316889095073556, 50.37956971798566, 38.51733141671543, 14.262183383064622, 5.050412223802015, 10.391315907239585, 92.98133332472621, 109.01681276723616, 107.57677757156662, 65.17331094454556, 38.45714045876865, 26.489926954124385, 75.78606146481745, 62.62642188514866, 107.75808847282715, 47.201472066537214, 15.224220158645599, 88.61261378065052, 5.246507166208237, 9.196183800333909, 25.5433698764445, 74.38245817399869, 72.59973650038027, 77.72701736203068, 87.41163007485763, 35.9562353589701, 52.08143432446701, 109.53785943024516, 13.030440405708367, 17.487146313041478, 41.042352685427005, 97.87246164702474, 62.34486879837736, 38.31870887799064, 20.280999746874148, 85.14945489347565, 7.635625081324512, 34.161714483319514, 22.400516908770786, 52.82574428391334, 65.77786237932874, 156.86033516204452, 72.79349761874583, 142.3151378353625, 233.5330024470549, 5.2203407063935074, 123.7779831063507, 33.88265626281833, 23.46587028439319, 36.80403613385671, 8.606679572351961, 97.66796020796627, 87.40365545968612, 18.464249177237146, 88.64058987993567, 24.724413704303355, 32.01993004030046, 12.135625574283186, 111.03745894821475, 28.560903146677177, 60.71702151980888, 58.63033881120785, 22.9922407513341, 21.695210465384772, 16.84382828732733, 34.46430373014289, 72.2521372167779, 91.41723436965684, 211.40120643536588, 43.2932570111605, 48.99300365345403, 95.91261646263202, 9.177647166185327, 128.49438608132056, 60.1128605065641, 19.970274137891288, 304.59736704998477, 12.167539500027882, 21.75076481220733, 19.406917777809213, 79.89871699143055, 57.34471705950752, 45.917912871055954, 80.7559504153946, 185.39896762273395, 13.748685547157848, 20.276649402702382, 122.71099181489973, 40.816159299512805, 57.61074094123394, 16.470648868376493, 19.71684519728231, 10.767044093380509, 47.6203582457846, 14.555605887664182, 77.06464152766388, 12.696578649140196, 77.73410738845169, 8.325066835398632, 72.35507476471018, 22.715771897096168, 162.58787937064534, 19.180345691539177, 58.073722015515756, 11.951487697543904, 63.381894735702424, 7.548881960412486, 13.14133130234406, 68.71958501549082, 53.83584026335273, 8.473532988363099, 22.326460173759187, 5.705549334331559, 78.17408898892896, 25.468255842330855, 10.71092731355036, 13.632211425926862, 8.568961975309552, 180.76506336449532, 83.53662257697995, 109.51898888462628, 7.270365270042803, 18.002409288274954, 121.28398442527717, 48.42043038639474, 25.80537217468133, 158.72491843572791, 54.77437089312585, 75.68020319506877, 9.183297854807824, 51.80007553063447, 8.204838522401333, 37.36375274347866, 29.08008488932845, 11.119943796007474, 211.37403816321228, 37.49607251184844, 27.29079211988199, 16.019604834639615, 16.81871168730808, 162.24795618616378, 348.0208824018029, 11.512538419311905, 84.95693028645478, 62.65182410166494, 61.814011075056726, 7.05996500998451, 21.846587849008756, 86.28542315174298, 78.0221493605246, 9.00742927464562, 21.360821169463694, 18.402261284517184, 39.446510400864625, 16.429311389551888, 43.05576062745309, 39.23774738371048, 17.43963231424453, 13.926713281780472, 27.72959521880666, 113.64483357506909, 159.2781336004997, 7.976879425568566, 92.74080512230546, 246.9537892285316, 19.58411955664895, 116.95838285876593, 16.608827621611027, 23.82849815299224, 25.912249373524652, 17.482252030058426, 31.653362635410506, 25.01550224285874, 20.95775182040264, 101.06927365944459, 30.617737634293213, 9.04249233137478, 106.335675908925, 73.53989227980753, 6.010301376745381, 12.649533435369984, 146.9303231401622, 110.04507905263563, 43.677188294266, 29.286999113927482, 16.213836664877828, 33.240068874474346, 6.18920465537694, 19.130284831495686, 136.2562310080708, 126.7274729281073, 22.95495896703929, 7.004520810990987, 140.78092181156836, 28.085399735984392, 45.43943228508579, 5.648259787593622, 69.62042488994938, 47.739556490555444, 114.19971017845471, 71.53522721889205, 44.46902136650733, 9.425323360214371, 57.9810514966439, 114.66597044116139, 41.12324250059271, 23.175964011569718, 79.5861969635528, 53.0532927746423, 62.727365927343456, 109.94989438830285, 15.658905834381851, 21.778542632880804, 36.36632946188011, 10.222091554119853, 6.7873618683553, 28.880471070308236, 131.67131667382372, 52.619429656445, 28.06603059178011, 157.21840007160353, 50.803075046710006, 19.95689149122119, 21.988561441577644, 12.551705730635698, 5.436004994134631, 48.716992764783406, 13.180159338230013, 156.03150992438563, 123.65860071903353, 58.321931811172945, 5.922546506948526, 13.924990471995743, 38.9451241884922, 6.259461496906064, 30.37976501613859, 34.78347690602414, 33.518193600280554, 5.772827310837642, 88.51492599116492, 9.565603052899677, 11.40817185140431, 12.52468523380964, 32.391526589132496, 59.77838631147118, 5.10018779349841, 14.397051408508585, 22.85011628409739, 106.19166152193783, 19.915839772133424, 44.988250154192315, 34.354753470506274, 69.4366667891856, 80.49454288386778, 65.04169189006475, 120.95514196717421, 611.4128352883533, 32.93152462355539, 23.068341812207123, 5.997137705745988, 34.75113389208721, 20.97452648198206, 20.03801649051698, 25.641992317704887, 95.16720337214147, 22.45960127705453, 437.88582059585286, 11.675223798831418, 58.15562503068316, 11.444024528109768, 6.937374818795225, 26.374443818529542, 73.81957967964568, 103.24728354463329, 20.168593988614965, 58.23224495807202, 34.39429696608186, 94.81145069717866, 21.742586468418516, 30.913765264350925, 10.226784405376707, 23.107715632397078, 20.5675082270712, 79.52161818918337, 90.48449923903593, 8.227089622147698, 6.952172153816562, 63.16990399274361, 16.044133783858697, 54.84976493540472, 8.812064684851697, 41.791993607787646, 6.144525099978994, 92.93738921275857, 15.908261090293172, 13.29835014140322, 9.59569761540489, 43.256685568109944, 7.024742463942281, 35.92521821170954, 14.576421070416197, 81.1560630541061, 33.79438771842482, 7.3684933456447, 24.431637367536588, 57.874373230253326, 67.57188236297732, 47.66323954987702, 90.8041140108025, 51.439682039556146, 51.731592215243325, 45.06556246986837, 10.254086518435326, 22.054003786842458, 8.081482865107146, 8.98044614607415, 66.80082676830175, 10.113797372636842, 6.9422301854809545, 83.7398824194271, 10.335155359737767, 22.02517719826753, 9.96656889405723, 56.579637635076196, 57.52168270548949, 22.0480743143437, 5.258909977878666, 144.47759392387496, 79.09108007978206, 69.24270188279641, 35.89004833355028, 60.26363541158082, 6.183840691371277, 31.97344266911736, 25.597662899510457, 67.9329868498281, 5.28261375676724, 20.734021386319046, 6.096347921549088, 67.14936939666518, 10.882703339478807, 5.169400099528197, 19.023574036903337, 9.303207850954688, 26.6600230183077, 21.860610439885672, 23.442237148723738, 127.29885380051086, 17.4648295190347, 42.68881929337756, 11.700428001468545, 29.921093116833646, 5.696499699919835, 19.2443317893341, 182.5170317848979, 29.351825875242476, 7.384217203800981, 158.02435935390545, 93.6796375627975, 49.734829366015376, 26.218616400446056, 70.48462904966239, 21.90022916982915, 5.884769877879415, 16.967991716379746, 22.625678562515667, 19.949479680015866, 23.960426107520373, 107.73796671430168, 45.33241198334421, 16.551612949877065, 44.9481301203236, 26.995265260177415, 40.769819419450286, 29.661874115851987, 38.128562972575686, 5.164584629287279, 25.23612596315981, 9.034669612732687, 105.35188320443694, 11.572188841072608, 62.41001496027766, 17.05613234304937, 15.67957286764548, 42.22906164301985, 41.40982766779813, 7.678404078917541, 105.6860288700245, 34.8983117569608, 13.12299864359369, 14.113481194007631, 87.43656831514963, 7.007152878075235, 76.97407925764693, 9.732250889280134, 26.929013880208228, 76.31602589461515, 77.79864999731775, 97.61357800441296, 21.620174771384054, 99.81059895227935, 71.19413335457321, 23.88995214378473, 21.27003325796731, 92.23502632421537, 11.054757102750854, 117.32421406600638, 24.871945528717053, 15.74115208184082, 107.90607263392289, 32.92274559824446, 5.141751749210712, 14.788473040988752, 55.17857078267191, 6.308394520941128, 6.6929168546123465, 40.176013510789836, 31.472872141328132, 161.9471794212171, 67.06244420388275, 61.39961217649998, 101.56502637727972, 46.73385749013866, 29.978349784088167, 55.30800203219855, 134.9125742941927, 24.704969679476264, 31.644982610625135, 25.42918615315442, 6.700093784072196, 18.612093252725106, 26.88599035105503, 49.864008982291196, 83.35506395713085, 67.23198892872176, 38.30576017927969, 23.82450356142916, 33.90428826276521, 18.172804751843756, 173.90643258183582, 101.22909867546974, 12.303051119172679, 18.976830191959877, 68.59231350290744, 153.606440785948, 85.63466852403086, 25.407117093849344, 64.1794818355354, 31.67049761975297, 23.752084294829487, 109.40375691974054, 246.93007837495756, 13.933803505673657, 36.39785808745186, 41.745722473090936, 34.62831860015093, 18.93191654541656, 17.59793771899404, 25.68149534796574, 43.347524211242, 74.88563508298893, 45.22936225302462, 66.62927528803118, 118.8088923752862, 8.722983434509487, 116.95347218458855, 41.768967289424516, 54.55988050220941, 20.278028923471513, 48.20969421488407, 32.31265642337569, 15.350711904886683, 85.67224521712531, 84.24951964685421, 5.161379431657956, 76.28379423604815, 33.1024970761187, 133.15944206051464, 91.09716623163561, 43.26042252568403, 11.712283702598027, 16.75597725330885, 33.45915741767293, 139.43996894880084, 19.158661069567664, 5.032455309033085, 29.58190433628818, 60.302959311488415, 21.451901482309, 7.383592000411834, 53.79639980374819, 125.50334648992018, 110.05042240838097, 13.80468840917342, 11.233641694111, 136.54139488580356, 30.58340829934098, 57.71066906076764, 11.252004605354923, 47.82266105891013, 26.46997064619461, 5.082945307235388, 122.16335575401973, 8.339207306507818, 11.987078167692127, 57.96118894037936, 50.94346112312344, 7.032659541177036, 60.41548987421088, 160.19421742491545, 42.77496277006378, 60.93280250334797, 10.970897776940634, 96.23301908766499, 63.046752910794055, 58.62633110801654, 9.827740970297738, 18.809408384890244, 20.119658655171765, 29.974775599673983, 64.8488788100482, 17.026645795901803, 44.96103701925489, 40.194098688191, 104.24190175524134, 30.897882771881072, 35.152632493071216, 85.48149784118453, 67.45284002160504, 9.141389828315232, ...])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3470703.125, 3480925.0, 3598635.9375, 3653518.75, 3739956.25, 3801859.375, 3801917.1875, 3816970.3125, 3831207.8125, 3905751.5625, 3939035.9375, 3939779.6875, 3952146.875, 3995578.125, 4061498.4375, 4099175.0, 4335209.375, 4352259.375, 4366034.375, 4380714.0625, 4403509.375, 4422071.875, 4440512.5, 4440539.0625, 4444476.5625, 4453784.375, 4465164.0625, 4491889.0625, 4491940.625, 4498448.4375, 4502106.25, 4511367.1875, 4515553.125, 4519835.9375, 4529248.4375, 4538576.5625, 4551434.375, 4560790.625, 4561012.5, 4561021.875, 4566389.0625, 4567720.3125, 4598221.875, 4623303.125, 4639548.4375, 4657598.4375, 4722226.5625, 4728223.4375, 4735293.75, 4750329.6875, 4753018.75, 4816007.8125, 4831603.125, 4831620.3125, 4858807.8125, 4904787.5, 4938901.5625, 5004753.125, 5042281.25, 5838914.0625, 5838939.0625, 5868325.0, 5878576.5625, 5884629.6875, 5909895.3125, 5910268.75, 5913950.0, 5920875.0, 5922079.6875, 5923128.125, 5924895.3125, 5926420.3125, 5931232.8125, 5931289.0625, 5931648.4375, 5937810.9375, 5950868.75, 5952243.75, 5956689.0625, 5975262.5, 6010259.375, 6034350.0, 6034367.1875, 6035717.1875, 6045203.125, 6045373.4375, 6107350.0, 6107462.5, 6116495.3125, 6117214.0625, 6118854.6875, 6119615.625, 6122517.1875, 6122528.125, 6123198.4375, 6137473.4375, 6140943.75, 6153803.125, 6154179.6875, 6164606.25, 6172214.0625, 6178732.8125, 6179531.25, 6181407.8125, 6205753.125, 6211875.0, 6213450.0, 6225665.625, 6231765.625, 6236398.4375, 6236400.0, 6236404.6875, 6260998.4375, 6270893.75, 6276895.3125, 6284298.4375, 6285343.75, 6290525.0, 6292803.125, 6294256.25, 6294520.3125, 6295495.3125, 6295501.5625, 6295828.125, 6296404.6875, 6296918.75, 6297087.5, 6298154.6875, 6298265.625, 6301901.5625, 6310104.6875, 6311875.0, 6328570.3125, 6328709.375, 6329470.3125, 6337890.625, 6338879.6875, 6348040.625, 6348187.5, 6354000.0, 6359185.9375, 6362120.3125, 6365779.6875, 6365871.875, 6369507.8125, 6373795.3125, 6374860.9375, 6374959.375, 6375948.4375, 6381031.25, 6384264.0625, 6384392.1875, 6385601.5625, 6386307.8125, 6386879.6875, 6390520.3125, 6395662.5, 6396964.0625, 6399181.25, 6401507.8125, 6407946.875, 6408001.5625, 6410643.75, 6410648.4375, 6411650.0, 6412393.75, 6412451.5625, 6416378.125, 6416635.9375, 6418126.5625, 6419359.375, 6419471.875, 6422156.25, 6426595.3125, 6431257.8125, 6431275.0, 6443267.1875, 6444823.4375, 6445287.5, 6450365.625, 6453176.5625, 6453195.3125, 6453445.3125, 6453525.0, 6453953.125, 6454562.5, 6455629.6875, 6456620.3125, 6462695.3125, 6462704.6875, 6464318.75, 6465598.4375, 6469204.6875, 6473532.8125, 6475532.8125, 6477868.75, 6477882.8125, 6479181.25, 6481104.6875, 6481543.75, 6482462.5, 6483679.6875, 6485120.3125, 6490718.75, 6490878.125, 6493290.625, 6495804.6875, 6503553.125, 6510087.5, 6512015.625, 6512403.125, 6513850.0, 6514632.8125, 6516128.125, 6518473.4375, 6521920.3125, 6522589.0625, 6524535.9375, 6528795.3125, 6529454.6875, 6531757.8125, 6532001.5625, 6532228.125, 6533046.875, 6533890.625, 6533993.75, 6534092.1875, 6534104.6875, 6534398.4375, 6534485.9375, 6534540.625, 6534870.3125, 6535714.0625, 6536657.8125, 6537370.3125, 6538042.1875, 6538471.875, 6538646.875, 6538751.5625, 6538987.5, 6541532.8125, 6541592.1875, 6541934.375, 6542043.75, 6543900.0, 6543946.875, 6548914.0625, 6550000.0, 6555868.75, 6556376.5625, 6556914.0625, 6559601.5625, 6559757.8125, 6562078.125, 6562137.5, 6562659.375, 6563259.375, 6563328.125, 6563348.4375, 6563589.0625, 6563610.9375, 6564532.8125, 6565235.9375, 6570531.25, 6571701.5625, 6572684.375, 6573951.5625, 6573985.9375, 6575054.6875, 6575329.6875, 6575429.6875, 6576720.3125, 6577021.875, 6578181.25, 6578235.9375, 6579962.5, 6580465.625, 6580521.875, 6580626.5625, 6582337.5, 6583095.3125, 6586312.5, 6587685.9375, 6587742.1875, 6590193.75, 6590210.9375, 6590318.75, 6591201.5625, 6593070.3125, 6594659.375, 6594965.625, 6597537.5, 6599832.8125, 6600621.875, 6600692.1875, 6601557.8125, 6601810.9375, 6602834.375, 6605082.8125, 6605095.3125, 6605914.0625, 6607081.25, 6607164.0625, 6607885.9375, 6609381.25, 6609635.9375, 6617843.75, 6622465.625, 6626889.0625, 6628828.125, 6631379.6875, 6633635.9375, 6634189.0625, 6637800.0, 6638217.1875, 6639678.125, 6648617.1875, 6648714.0625, 6649259.375, 6649587.5, 6654384.375, 6656885.9375, 6675267.1875, 6676045.3125, 6679459.375, 6688209.375, 6688637.5, 6688857.8125, 6692514.0625, 6694685.9375, 6700328.125, 6700721.875, 6701975.0, 6703598.4375, 6707439.0625, 6716284.375, 6717662.5, 6722821.875, 6726289.0625, 6727328.125, 6729115.625, 6729201.5625, 6733485.9375, 6733845.3125, 6734176.5625, 6734270.3125, 6735014.0625, 6737450.0, 6737739.0625, 6737745.3125, 6738331.25, 6739204.6875, 6740117.1875, 6743189.0625, 6744381.25, 6746832.8125, 6748928.125, 6758453.125, 6765581.25, 6766345.3125, 6767529.6875, 6769968.75, 6773750.0, 6775525.0, 6778418.75, 6785853.125, 6788176.5625, 6789664.0625, 6791737.5, 6797053.125, 6798537.5, 6802812.5, 6806085.9375, 6807973.4375, 6808314.0625, 6809484.375, 6810926.5625, 6810946.875, 6817709.375, 6820435.9375, 6821660.9375, 6824609.375, 6824776.5625, 6827606.25, 6827667.1875, 6830176.5625, 6830321.875, 6832028.125, 6832303.125, 6832750.0, 6832782.8125, 6833887.5, 6834196.875, 6834754.6875, 6835962.5, 6837501.5625, 6839189.0625, 6839371.875, 6840176.5625, 6841854.6875, 6841890.625, 6841989.0625, 6842385.9375, 6843281.25, 6843492.1875, 6844131.25, 6844210.9375, 6844256.25, 6845025.0, 6845179.6875, 6846879.6875, 6853315.625, 6877454.6875, 6877467.1875, 6879728.125, 6880060.9375, 6883640.625, 6885059.375, 6885073.4375, 6891931.25, 6892937.5, 6893018.75, 6895690.625, 6897728.125, 6918784.375, 6918787.5, 6920029.6875, 6921639.0625, 6922851.5625, 6925350.0, 6925621.875, 6927460.9375, 6932770.3125, 6934234.375, 6940609.375, 6941910.9375, 6941950.0, 6942523.4375, 6943695.3125, 6945593.75, 6945820.3125, 6947032.8125, 6947296.875, 6952042.1875, 6960453.125, 6977898.4375, 6980271.875, 6981100.0, 6981117.1875, 6981534.375, 6984995.3125, 6986426.5625, 6986484.375, 6986551.5625, 6987107.8125, 6987556.25, 6987846.875, 6987859.375, 6988218.75, 6988842.1875, 6988959.375, 6988985.9375, 6989051.5625, 6989057.8125, 6989379.6875, 6989409.375, 6990923.4375, 6991404.6875, 6991453.125, 6991573.4375, 6991942.1875, 6992051.5625, 6993012.5, 6993681.25, 6993881.25, 6994371.875, 6994445.3125, 6994925.0, 6995023.4375, 6995107.8125, 6995309.375, 6995473.4375, 6996628.125, 6996981.25, 6997337.5, 6997748.4375, 6997798.4375, 6997889.0625, 6998454.6875, 6998489.0625, 6998540.625, 6998542.1875, 6999332.8125, 6999451.5625, 6999521.875, 7001189.0625, 7001204.6875, 7001207.8125, 7001390.625, 7003003.125, 7003445.3125, 7005314.0625, 7005335.9375, 7009192.1875, 7016895.3125, 7020592.1875, 7022682.8125, 7032073.4375, 7032923.4375, 7035607.8125, 7040003.125, 7042409.375, 7043978.125, 7045521.875, 7045550.0, 7045654.6875, 7047379.6875, 7051001.5625, 7052907.8125, 7063760.9375, 7080603.125, 7088414.0625, 7097507.8125, 7111150.0, 7118500.0, 7122682.8125, 7123173.4375, 7123182.8125, 7129781.25, 7130006.25, 7130442.1875, 7150948.4375, 7151957.8125, 7152945.3125, 7153857.8125, 7157454.6875, 7158378.125, 7158442.1875, 7160437.5, 7167154.6875, 7167593.75, 7168675.0, 7169735.9375, 7169871.875, 7170834.375, 7176632.8125, 7180800.0, 7180803.125, 7186521.875, 7186787.5, 7190051.5625, 7191351.5625, 7192045.3125, 7192112.5, 7196346.875, 7199087.5, 7199139.0625, 7199753.125, 7200848.4375, 7201506.25, 7201720.3125, 7202951.5625, 7202957.8125, 7202971.875, 7202975.0, 7203603.125, 7203700.0, 7203759.375, 7204253.125, 7206335.9375, 7206928.125, 7209689.0625, 7222415.625, 7224040.625, 7228606.25, 7229360.9375, 7230101.5625, 7232987.5, 7235026.5625, 7236795.3125, 7255696.875, 7257245.3125, 7261845.3125, 7262132.8125, 7264621.875, 7267503.125, 7268381.25, 7269864.0625, 7271757.8125, 7278056.25, 7282176.5625, 7286479.6875, 7301062.5, 7306542.1875, 7308003.125, 7312517.1875, 7313971.875, 7318020.3125, 7327485.9375, 7332106.25, 7344387.5, 7353012.5, 7355571.875, 7370806.25, 7374682.8125, 7378696.875, 7378742.1875, 7380267.1875, 7384067.1875, 7384193.75, 7387600.0, 7391070.3125, 7394920.3125, 7396345.3125, 7397000.0, 7398423.4375, 7399346.875, 7399812.5, 7399842.1875, 7399860.9375, 7400406.25, 7400467.1875, 7401298.4375, 7401298.4375, 7401443.75, 7401820.3125, 7401942.1875, 7402331.25, 7402517.1875, 7402585.9375, 7403037.5, 7403042.1875, 7403589.0625, 7403900.0, 7404025.0, 7404831.25, 7405089.0625, 7405331.25, 7405571.875, 7406167.1875, 7406365.625, 7406421.875, 7406443.75, 7406571.875, 7406625.0, 7406628.125, 7406768.75, 7406898.4375, 7406985.9375, 7407554.6875, 7407560.9375, 7409626.5625, 7409710.9375, 7410323.4375, 7410696.875, 7410720.3125, 7410920.3125, 7411446.875, 7412246.875, 7412610.9375, 7413185.9375, 7413284.375, 7413585.9375, 7415032.8125, 7415071.875, 7418085.9375, 7418142.1875, 7418146.875, 7418939.0625, 7420471.875, 7420496.875, 7421318.75, 7421784.375, 7422079.6875, 7422989.0625, 7423114.0625, 7423535.9375, 7423620.3125, 7423739.0625, 7424296.875, 7424368.75, 7424629.6875, 7424981.25, 7425031.25, 7425132.8125, 7425151.5625, 7425226.5625, 7425331.25, 7425989.0625, 7427351.5625, 7428912.5, 7429190.625, 7430017.1875, 7430881.25, 7430885.9375, 7430940.625, 7431781.25, 7431784.375, 7431870.3125, 7432498.4375, 7433459.375, 7433464.0625, 7437148.4375, 7438045.3125, 7438800.0, 7439235.9375, 7439407.8125, 7442600.0, 7444151.5625, 7451131.25, 7452737.5, 7453160.9375, 7453285.9375, 7453464.0625, 7454062.5, 7454131.25, 7455009.375, 7455282.8125, 7455856.25, 7456696.875, 7457203.125, 7459103.125, 7459659.375, 7459793.75, 7463287.5, 7464278.125, 7464387.5, 7466350.0, 7469775.0, 7469829.6875, 7470676.5625, 7471781.25, 7479200.0, 7479704.6875, 7481104.6875, 7482906.25, 7488018.75, 7490910.9375, 7492776.5625, 7504732.8125, 7506354.6875, 7506360.9375, 7507032.8125, 7509010.9375, 7517221.875, 7520279.6875, 7539343.75, 7542398.4375, 7542471.875, 7542475.0, 7542571.875, 7542679.6875, 7542687.5, 7542807.8125, 7543168.75, 7543812.5, 7543864.0625, 7545325.0, 7545520.3125, 7545715.625, 7546445.3125, 7546448.4375, 7546454.6875, 7547245.3125, 7547459.375, 7547767.1875, 7548000.0, 7548103.125, 7548165.625, 7548456.25, 7548637.5, 7548821.875, 7548914.0625, 7549021.875, 7549201.5625, 7549326.5625, 7549345.3125, 7549782.8125, 7550331.25, 7550503.125, 7550539.0625, 7551343.75, 7551512.5, 7552079.6875, 7552206.25, 7552321.875, 7552971.875, 7553295.3125, 7554565.625, 7554662.5, 7554762.5, 7554907.8125, 7554945.3125, 7555431.25, 7556512.5, 7558326.5625, 7559665.625, 7570739.0625, 7571868.75, 7575484.375, 7575648.4375, 7575720.3125, 7575804.6875, 7577882.8125, 7577887.5, 7579265.625, 7579795.3125, 7579984.375, 7580779.6875, 7581193.75, 7581590.625, 7581696.875, 7582870.3125, 7583057.8125, 7584373.4375, 7584904.6875, 7587115.625, 7594201.5625, 7598643.75, 7602770.3125, 7603598.4375, 7605089.0625, 7605134.375, 7605914.0625, 7606571.875, 7607732.8125, 7607981.25, 7616067.1875, 7623285.9375, 7625106.25, 7626532.8125, 7628393.75, 7630475.0, 7631223.4375, 7631621.875, 7649718.75, 7658464.0625, 7664843.75, 7677440.625, 7679860.9375, 7682253.125, 7685754.6875, 7696390.625, 7697290.625, 7700509.375, 7700785.9375, 7700865.625, 7706804.6875, 7712309.375, 7713662.5, 7719706.25, 7720771.875, 7757917.1875, 7759796.875, 7760037.5, 7760629.6875, 7762906.25, 7763150.0, 7763242.1875, 7763879.6875, 7764178.125, 7764412.5, 7764415.625, 7764895.3125, 7766042.1875, 7766298.4375, 7767123.4375, 7767359.375, 7767546.875, 7767732.8125, 7767743.75, 7767879.6875, 7768082.8125, 7768085.9375, 7768504.6875, 7768779.6875, 7769778.125, 7769887.5, 7769901.5625, 7769945.3125, 7769968.75, 7770259.375, 7770329.6875, 7770342.1875, 7770381.25, 7770573.4375, 7770579.6875, 7770607.8125, 7770623.4375, 7770917.1875, 7770925.0, 7771493.75, 7771523.4375, 7771650.0, 7771703.125, 7772079.6875, 7776275.0, 7781475.0, 7781792.1875, 7783450.0, 7787321.875, 7787773.4375, 7788987.5, 7789923.4375, 7790475.0, 7792484.375, 7792971.875, 7799373.4375, 7801475.0, 7802045.3125, 7802546.875, 7802665.625, 7803865.625, 7803998.4375, 7805984.375, 7808557.8125, 7813931.25, 7828271.875, 7862162.5, 7863454.6875, 7864765.625, 7865290.625, 7866617.1875, 7866756.25, 7866757.8125, 7868154.6875, 7868167.1875, 7868465.625, 7868681.25, 7869381.25, 7869670.3125, 7869812.5, 7869987.5, 7870084.375, 7870825.0, 7870915.625, 7870923.4375, 7871015.625, 7872685.9375, 7881601.5625, 7883218.75, 7884246.875, 7892050.0, 7893835.9375, 7893890.625, 7894515.625, 7894953.125, 7895953.125, 7899523.4375, 7899595.3125, 7900356.25, 7902392.1875, 7903342.1875, 7918568.75, 7918996.875, 7919084.375, 7919198.4375, 7929390.625, 7930764.0625, 7939103.125, 7939123.4375, 7941054.6875, 7942542.1875, 7942778.125, 7944251.5625, 7944532.8125, 7945579.6875, 7946056.25, 7947525.0, 7947751.5625, 7952901.5625, 7952968.75, 7953270.3125, 7954014.0625, 7959406.25, 7963553.125, 7975454.6875, 7975643.75, 7990943.75, 7991581.25, 7992701.5625, 7992789.0625, 7993678.125, 7993996.875, 7995153.125, 7995551.5625, 7997607.8125, 7997617.1875, 7997660.9375, 7998934.375, 8000775.0, 8002462.5, 8006862.5, 8033621.875, 8034043.75, 8039560.9375, 8039589.0625, 8042185.9375, 8048268.75, 8048740.625, ...], [45.39907160677232, 58.43472597827107, 54.13986925637357, 89.12066462801411, 42.31772412446346, 44.334811280756355, 5.039176232203347, 5.998134463936753, 45.86931391332472, 84.75723619386078, 17.776354316266083, 13.483571108427403, 26.221404666351518, 122.39943040280984, 6.617917723442894, 6.467919871230439, 5.631474168321943, 39.62706428287824, 15.071489453073063, 84.36847826439555, 22.102515067820722, 6.216946875875842, 9.868371161867207, 10.322389149980804, 62.79481252107678, 93.3164113725934, 17.54321205886551, 12.152756737713236, 5.093415425987701, 6.376425591099246, 8.405526394965008, 74.97857949686866, 44.42105796148373, 29.4333807339959, 7.44665754383101, 29.082651627926793, 36.500974862091454, 12.555573518695828, 10.942724627448971, 7.588849913883986, 19.265618039863664, 9.809752109315776, 10.8106488535208, 80.32804657115567, 39.30421997334413, 16.6046771494948, 10.846970809067779, 16.143030242617993, 9.621170738191072, 6.45128784867267, 29.07133178912585, 10.544364280910047, 10.960144356161543, 59.14702833810563, 9.847484428910874, 41.4362725105399, 19.777389623817008, 26.1335112794383, 6.59164146767698, 13.939314326181142, 95.82585115343834, 23.38383925029865, 48.384512544701074, 78.15949706271037, 48.56458390501364, 5.475348409468644, 20.003753260607663, 26.05833028013529, 22.2631641647619, 59.94786198705246, 19.383723204735603, 51.56502116682048, 37.29317231216768, 6.1591816033042175, 181.55652163621295, 66.91443961230476, 31.526858619254003, 75.95028961562377, 65.21755936792128, 16.49099155461073, 55.00423227520606, 58.93448126855243, 52.57303396798196, 5.217553170091159, 6.297038156735923, 55.762434798435486, 122.59244627234175, 5.567883196850918, 36.64121541488077, 9.123455006903669, 6.155750364387428, 115.11362335820226, 14.155945289610052, 13.896946396998004, 109.52134482320002, 6.521761124976892, 28.894246811617883, 17.701096686631846, 6.841739455429855, 115.98313440439134, 87.99904748472544, 104.12694747318481, 38.645243162261444, 7.5688279260069775, 18.211069104169813, 122.5099663949049, 56.258755716947455, 42.31201243482917, 85.38430680996511, 28.627017323532996, 78.52249036185356, 22.257552213644875, 8.337888736273742, 10.884909902192174, 73.42403188280775, 18.94452891252995, 64.69291889600817, 53.75845680785268, 200.4307381448822, 15.403083707772161, 42.345565772351115, 10.792066111075812, 48.952598545524296, 82.04082184310545, 17.459621910125318, 22.578644341523503, 32.518121959953405, 18.596584589046813, 23.562605590381438, 66.2884170327356, 7.150933068576006, 74.30785344515007, 10.707748350875635, 59.1317663370279, 8.145896436792253, 39.19595319609624, 22.69215664318009, 6.756174477758828, 86.82839859316316, 23.55519733864917, 24.66577642177038, 8.245412884681722, 41.888715122311275, 5.177433677294963, 19.40547904648074, 5.395565708163598, 7.059264146471218, 5.156112256451142, 24.265766943295723, 25.07708309080054, 45.60033010662322, 17.213412925676487, 88.4263384452878, 62.19277057868756, 18.004235835022033, 9.285092982644406, 37.843745776663624, 21.201363741903585, 6.189170000128272, 17.632239923206257, 12.158797066862332, 5.387424293221695, 19.956927791612976, 31.544750338157037, 19.833414538623572, 65.70790924322222, 5.9092602237596985, 9.053405579304599, 88.24950473840443, 7.046788741057286, 52.259797094185, 32.91496672789286, 5.5568869593321235, 31.668638883872763, 14.004768803728018, 12.447933234196881, 91.88842608447418, 50.395568154099, 53.13890747925048, 12.856285612819297, 89.23664603720526, 58.09053308802784, 21.94958266811906, 17.63833934598756, 65.29770703729992, 69.71905195484634, 30.90831948911661, 15.933674677695823, 22.296600553194374, 60.34363931774433, 10.461005516721935, 5.0454657469302715, 50.0406825842092, 24.414766266230497, 9.893918438928216, 15.744856432955881, 10.643970707191222, 21.10846471513952, 136.64393338770643, 6.412187549294856, 22.571910458958776, 81.54459522258487, 16.673840323265726, 6.0985001814627715, 30.4960606464605, 163.44339035452413, 18.918455278445634, 44.91867856895598, 27.909485983714127, 29.521661742604593, 23.275619798221083, 57.28065002286976, 29.385702597457005, 7.4102310088436685, 102.6445590302597, 12.474542145873293, 21.56091987854063, 118.42552625068336, 153.20306438310783, 42.855443346989965, 96.45129324467587, 48.67620800629379, 53.447686966167716, 75.11213513887965, 45.31277737228109, 71.88347432463344, 23.646482817622456, 26.69004483286118, 160.01336347064196, 60.671913085868574, 6.388936207090076, 86.339830885314, 68.67242066656685, 20.332071743757655, 192.90851180289417, 13.50086638863289, 18.73622133159905, 61.92237559855329, 25.657698742525973, 57.16394346391178, 7.231209516791493, 94.66997077007252, 65.28848137767659, 115.5350612809463, 105.97551059442645, 85.70309810791288, 76.99501375011621, 21.532455954106926, 186.96168205187627, 84.84958542742243, 19.59833539313485, 130.41424353371184, 20.606194671347275, 44.53396692833216, 5.147449868012342, 27.402908760491336, 41.23654180808661, 6.320732704602066, 5.066039091525922, 27.528909248391273, 11.940316350267715, 5.64966048804249, 15.34483684382827, 68.79196299466794, 49.4414407200757, 25.297074013286498, 25.80380395072883, 49.925189759728056, 203.41652578947088, 41.37034528561626, 19.413489063201066, 22.740068262720936, 69.91825689622371, 6.491557075161738, 5.245993923447846, 40.1234117878206, 61.13955880638159, 105.04512716729712, 24.02470375118437, 95.43278300318732, 38.86067762456453, 52.8405476598883, 5.597005252311804, 5.1062186712719075, 120.85197042249324, 37.04775804224042, 148.4827380106409, 6.054074193955454, 65.56498508992377, 60.06850119040066, 66.09593712981703, 6.289839946194938, 97.59979343563244, 67.75257948420965, 12.55574537977661, 19.579557538956394, 42.243455950252496, 93.94876955573942, 10.214295549138312, 86.26770136285887, 11.239175216115399, 93.89453134919016, 48.632434070424566, 56.31239966064271, 16.974253018025315, 13.284059282926744, 55.94175496813692, 45.328864203068946, 33.46019774074443, 33.67656916228863, 166.85197878594138, 35.078030617766686, 20.302544005692127, 13.688339736074413, 15.04452041373881, 173.10406082439712, 119.04374525554634, 5.536420164809132, 28.910807613527968, 11.798533416214562, 71.23856367948319, 36.56810007910176, 5.183970751278581, 6.271037943667835, 16.58279971874704, 99.16234110249987, 5.29110979192944, 63.82699222193669, 165.34392507816844, 56.558922170957445, 8.417177091023483, 60.70816415202323, 6.358755406015505, 10.224438558589181, 126.75165885212118, 33.45551341054375, 87.51603664614346, 41.34328960950465, 38.45863954065044, 30.71793378594299, 52.573576853859514, 12.405968439363845, 9.936772881970768, 11.516754586909748, 53.208626157814926, 5.1254533998877525, 33.38268387540789, 56.57641282180743, 25.377644201192204, 31.566780382558477, 54.60241753769744, 164.83431835179906, 6.041046825041573, 18.417918709146356, 87.33145765011676, 28.277067543023254, 50.450172829026606, 6.138004097805336, 10.930642203226215, 20.977485891403944, 21.271989860957593, 25.08747440137797, 46.577723697459646, 15.287785020210238, 15.68805122087268, 92.15269797873349, 71.33130413015576, 18.96887726153514, 91.48311912760275, 36.46875108484006, 38.0887057852149, 21.824647857830204, 30.153747873777245, 44.71174270140521, 100.81145934145079, 102.36991611454101, 11.274615133291446, 5.556098077682061, 11.931118094582532, 6.526406230879605, 81.51127708683586, 63.55272371555312, 36.29085351670314, 140.76784046351992, 5.469384465884944, 29.803411402502242, 43.557315029470246, 45.30427071432581, 8.089787249377078, 73.45350093798532, 31.186414730630744, 76.10207636595962, 51.928865988721824, 36.146434187043766, 72.97701412916159, 11.08585042027072, 40.21400977253449, 14.2150281952447, 14.058412170920274, 11.115480864398485, 16.201867928201565, 69.17968872420569, 11.281210534829352, 67.65173082092512, 49.242126947805446, 43.83393473092084, 7.7995722527754, 152.52156616805956, 19.984556308480105, 77.18302296339232, 8.048554268843452, 23.02039655196422, 128.0673188399014, 50.196296468777554, 9.34252793461419, 47.299820719892935, 18.782294796546264, 133.53990843514185, 67.38985519011622, 58.03749151043067, 21.14784692288888, 50.73945281382011, 13.460610761028006, 49.68708529325806, 35.935099630540044, 29.221446354052762, 34.65648861643882, 16.469711879910875, 97.20530485075014, 84.1989207251092, 25.445292421066412, 74.19879793961121, 31.65260505276899, 118.60926307237679, 6.201124605319172, 55.74620815433049, 21.715132998828715, 6.208759179210996, 11.716829216942017, 9.12402337275018, 22.578604468202663, 36.38174943106439, 59.61004846525074, 13.47185660808665, 87.03466393497298, 11.6279752627118, 5.350525528549183, 7.239744481714719, 51.47016199737295, 10.47602973334048, 17.09946966572913, 123.54978959282262, 35.51582230486843, 31.024608105432446, 8.78835583754656, 33.85602790624039, 51.527600122273405, 6.439514952237363, 15.531036282208104, 32.40325296215359, 20.829721709435912, 68.92696691691638, 5.171471823328056, 22.822942443731623, 61.68168187061277, 12.159356070165803, 60.373039087197085, 108.03092352171295, 146.72266425352174, 58.78546962661646, 42.316889095073556, 50.37956971798566, 38.51733141671543, 14.262183383064622, 5.050412223802015, 10.391315907239585, 92.98133332472621, 109.01681276723616, 107.57677757156662, 65.17331094454556, 38.45714045876865, 26.489926954124385, 75.78606146481745, 62.62642188514866, 107.75808847282715, 47.201472066537214, 15.224220158645599, 88.61261378065052, 5.246507166208237, 9.196183800333909, 25.5433698764445, 74.38245817399869, 72.59973650038027, 77.72701736203068, 87.41163007485763, 35.9562353589701, 52.08143432446701, 109.53785943024516, 13.030440405708367, 17.487146313041478, 41.042352685427005, 97.87246164702474, 62.34486879837736, 38.31870887799064, 20.280999746874148, 85.14945489347565, 7.635625081324512, 34.161714483319514, 22.400516908770786, 52.82574428391334, 65.77786237932874, 156.86033516204452, 72.79349761874583, 142.3151378353625, 233.5330024470549, 5.2203407063935074, 123.7779831063507, 33.88265626281833, 23.46587028439319, 36.80403613385671, 8.606679572351961, 97.66796020796627, 87.40365545968612, 18.464249177237146, 88.64058987993567, 24.724413704303355, 32.01993004030046, 12.135625574283186, 111.03745894821475, 28.560903146677177, 60.71702151980888, 58.63033881120785, 22.9922407513341, 21.695210465384772, 16.84382828732733, 34.46430373014289, 72.2521372167779, 91.41723436965684, 211.40120643536588, 43.2932570111605, 48.99300365345403, 95.91261646263202, 9.177647166185327, 128.49438608132056, 60.1128605065641, 19.970274137891288, 304.59736704998477, 12.167539500027882, 21.75076481220733, 19.406917777809213, 79.89871699143055, 57.34471705950752, 45.917912871055954, 80.7559504153946, 185.39896762273395, 13.748685547157848, 20.276649402702382, 122.71099181489973, 40.816159299512805, 57.61074094123394, 16.470648868376493, 19.71684519728231, 10.767044093380509, 47.6203582457846, 14.555605887664182, 77.06464152766388, 12.696578649140196, 77.73410738845169, 8.325066835398632, 72.35507476471018, 22.715771897096168, 162.58787937064534, 19.180345691539177, 58.073722015515756, 11.951487697543904, 63.381894735702424, 7.548881960412486, 13.14133130234406, 68.71958501549082, 53.83584026335273, 8.473532988363099, 22.326460173759187, 5.705549334331559, 78.17408898892896, 25.468255842330855, 10.71092731355036, 13.632211425926862, 8.568961975309552, 180.76506336449532, 83.53662257697995, 109.51898888462628, 7.270365270042803, 18.002409288274954, 121.28398442527717, 48.42043038639474, 25.80537217468133, 158.72491843572791, 54.77437089312585, 75.68020319506877, 9.183297854807824, 51.80007553063447, 8.204838522401333, 37.36375274347866, 29.08008488932845, 11.119943796007474, 211.37403816321228, 37.49607251184844, 27.29079211988199, 16.019604834639615, 16.81871168730808, 162.24795618616378, 348.0208824018029, 11.512538419311905, 84.95693028645478, 62.65182410166494, 61.814011075056726, 7.05996500998451, 21.846587849008756, 86.28542315174298, 78.0221493605246, 9.00742927464562, 21.360821169463694, 18.402261284517184, 39.446510400864625, 16.429311389551888, 43.05576062745309, 39.23774738371048, 17.43963231424453, 13.926713281780472, 27.72959521880666, 113.64483357506909, 159.2781336004997, 7.976879425568566, 92.74080512230546, 246.9537892285316, 19.58411955664895, 116.95838285876593, 16.608827621611027, 23.82849815299224, 25.912249373524652, 17.482252030058426, 31.653362635410506, 25.01550224285874, 20.95775182040264, 101.06927365944459, 30.617737634293213, 9.04249233137478, 106.335675908925, 73.53989227980753, 6.010301376745381, 12.649533435369984, 146.9303231401622, 110.04507905263563, 43.677188294266, 29.286999113927482, 16.213836664877828, 33.240068874474346, 6.18920465537694, 19.130284831495686, 136.2562310080708, 126.7274729281073, 22.95495896703929, 7.004520810990987, 140.78092181156836, 28.085399735984392, 45.43943228508579, 5.648259787593622, 69.62042488994938, 47.739556490555444, 114.19971017845471, 71.53522721889205, 44.46902136650733, 9.425323360214371, 57.9810514966439, 114.66597044116139, 41.12324250059271, 23.175964011569718, 79.5861969635528, 53.0532927746423, 62.727365927343456, 109.94989438830285, 15.658905834381851, 21.778542632880804, 36.36632946188011, 10.222091554119853, 6.7873618683553, 28.880471070308236, 131.67131667382372, 52.619429656445, 28.06603059178011, 157.21840007160353, 50.803075046710006, 19.95689149122119, 21.988561441577644, 12.551705730635698, 5.436004994134631, 48.716992764783406, 13.180159338230013, 156.03150992438563, 123.65860071903353, 58.321931811172945, 5.922546506948526, 13.924990471995743, 38.9451241884922, 6.259461496906064, 30.37976501613859, 34.78347690602414, 33.518193600280554, 5.772827310837642, 88.51492599116492, 9.565603052899677, 11.40817185140431, 12.52468523380964, 32.391526589132496, 59.77838631147118, 5.10018779349841, 14.397051408508585, 22.85011628409739, 106.19166152193783, 19.915839772133424, 44.988250154192315, 34.354753470506274, 69.4366667891856, 80.49454288386778, 65.04169189006475, 120.95514196717421, 611.4128352883533, 32.93152462355539, 23.068341812207123, 5.997137705745988, 34.75113389208721, 20.97452648198206, 20.03801649051698, 25.641992317704887, 95.16720337214147, 22.45960127705453, 437.88582059585286, 11.675223798831418, 58.15562503068316, 11.444024528109768, 6.937374818795225, 26.374443818529542, 73.81957967964568, 103.24728354463329, 20.168593988614965, 58.23224495807202, 34.39429696608186, 94.81145069717866, 21.742586468418516, 30.913765264350925, 10.226784405376707, 23.107715632397078, 20.5675082270712, 79.52161818918337, 90.48449923903593, 8.227089622147698, 6.952172153816562, 63.16990399274361, 16.044133783858697, 54.84976493540472, 8.812064684851697, 41.791993607787646, 6.144525099978994, 92.93738921275857, 15.908261090293172, 13.29835014140322, 9.59569761540489, 43.256685568109944, 7.024742463942281, 35.92521821170954, 14.576421070416197, 81.1560630541061, 33.79438771842482, 7.3684933456447, 24.431637367536588, 57.874373230253326, 67.57188236297732, 47.66323954987702, 90.8041140108025, 51.439682039556146, 51.731592215243325, 45.06556246986837, 10.254086518435326, 22.054003786842458, 8.081482865107146, 8.98044614607415, 66.80082676830175, 10.113797372636842, 6.9422301854809545, 83.7398824194271, 10.335155359737767, 22.02517719826753, 9.96656889405723, 56.579637635076196, 57.52168270548949, 22.0480743143437, 5.258909977878666, 144.47759392387496, 79.09108007978206, 69.24270188279641, 35.89004833355028, 60.26363541158082, 6.183840691371277, 31.97344266911736, 25.597662899510457, 67.9329868498281, 5.28261375676724, 20.734021386319046, 6.096347921549088, 67.14936939666518, 10.882703339478807, 5.169400099528197, 19.023574036903337, 9.303207850954688, 26.6600230183077, 21.860610439885672, 23.442237148723738, 127.29885380051086, 17.4648295190347, 42.68881929337756, 11.700428001468545, 29.921093116833646, 5.696499699919835, 19.2443317893341, 182.5170317848979, 29.351825875242476, 7.384217203800981, 158.02435935390545, 93.6796375627975, 49.734829366015376, 26.218616400446056, 70.48462904966239, 21.90022916982915, 5.884769877879415, 16.967991716379746, 22.625678562515667, 19.949479680015866, 23.960426107520373, 107.73796671430168, 45.33241198334421, 16.551612949877065, 44.9481301203236, 26.995265260177415, 40.769819419450286, 29.661874115851987, 38.128562972575686, 5.164584629287279, 25.23612596315981, 9.034669612732687, 105.35188320443694, 11.572188841072608, 62.41001496027766, 17.05613234304937, 15.67957286764548, 42.22906164301985, 41.40982766779813, 7.678404078917541, 105.6860288700245, 34.8983117569608, 13.12299864359369, 14.113481194007631, 87.43656831514963, 7.007152878075235, 76.97407925764693, 9.732250889280134, 26.929013880208228, 76.31602589461515, 77.79864999731775, 97.61357800441296, 21.620174771384054, 99.81059895227935, 71.19413335457321, 23.88995214378473, 21.27003325796731, 92.23502632421537, 11.054757102750854, 117.32421406600638, 24.871945528717053, 15.74115208184082, 107.90607263392289, 32.92274559824446, 5.141751749210712, 14.788473040988752, 55.17857078267191, 6.308394520941128, 6.6929168546123465, 40.176013510789836, 31.472872141328132, 161.9471794212171, 67.06244420388275, 61.39961217649998, 101.56502637727972, 46.73385749013866, 29.978349784088167, 55.30800203219855, 134.9125742941927, 24.704969679476264, 31.644982610625135, 25.42918615315442, 6.700093784072196, 18.612093252725106, 26.88599035105503, 49.864008982291196, 83.35506395713085, 67.23198892872176, 38.30576017927969, 23.82450356142916, 33.90428826276521, 18.172804751843756, 173.90643258183582, 101.22909867546974, 12.303051119172679, 18.976830191959877, 68.59231350290744, 153.606440785948, 85.63466852403086, 25.407117093849344, 64.1794818355354, 31.67049761975297, 23.752084294829487, 109.40375691974054, 246.93007837495756, 13.933803505673657, 36.39785808745186, 41.745722473090936, 34.62831860015093, 18.93191654541656, 17.59793771899404, 25.68149534796574, 43.347524211242, 74.88563508298893, 45.22936225302462, 66.62927528803118, 118.8088923752862, 8.722983434509487, 116.95347218458855, 41.768967289424516, 54.55988050220941, 20.278028923471513, 48.20969421488407, 32.31265642337569, 15.350711904886683, 85.67224521712531, 84.24951964685421, 5.161379431657956, 76.28379423604815, 33.1024970761187, 133.15944206051464, 91.09716623163561, 43.26042252568403, 11.712283702598027, 16.75597725330885, 33.45915741767293, 139.43996894880084, 19.158661069567664, 5.032455309033085, 29.58190433628818, 60.302959311488415, 21.451901482309, 7.383592000411834, 53.79639980374819, 125.50334648992018, 110.05042240838097, 13.80468840917342, 11.233641694111, 136.54139488580356, 30.58340829934098, 57.71066906076764, 11.252004605354923, 47.82266105891013, 26.46997064619461, 5.082945307235388, 122.16335575401973, 8.339207306507818, 11.987078167692127, 57.96118894037936, 50.94346112312344, 7.032659541177036, 60.41548987421088, 160.19421742491545, 42.77496277006378, 60.93280250334797, 10.970897776940634, 96.23301908766499, 63.046752910794055, 58.62633110801654, 9.827740970297738, 18.809408384890244, 20.119658655171765, 29.974775599673983, 64.8488788100482, 17.026645795901803, 44.96103701925489, 40.194098688191, 104.24190175524134, 30.897882771881072, 35.152632493071216, 85.48149784118453, 67.45284002160504, 9.141389828315232, ...])
caption, T_hit,count = 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, count, caption);
([3470703.125, 3480925.0, 3598635.9375, 3653518.75, 3739956.25, 3801859.375, 3801917.1875, 3816970.3125, 3831207.8125, 3905751.5625, 3939035.9375, 3939779.6875, 3952146.875, 3995578.125, 4061498.4375, 4099175.0, 4335209.375, 4352259.375, 4366034.375, 4380714.0625, 4403509.375, 4422071.875, 4440512.5, 4440539.0625, 4444476.5625, 4453784.375, 4465164.0625, 4491889.0625, 4491940.625, 4498448.4375, 4502106.25, 4511367.1875, 4515553.125, 4519835.9375, 4529248.4375, 4538576.5625, 4551434.375, 4560790.625, 4561012.5, 4561021.875, 4566389.0625, 4567720.3125, 4598221.875, 4623303.125, 4639548.4375, 4657598.4375, 4722226.5625, 4728223.4375, 4735293.75, 4750329.6875, 4753018.75, 4816007.8125, 4831603.125, 4831620.3125, 4858807.8125, 4904787.5, 4938901.5625, 5004753.125, 5042281.25, 5838914.0625, 5838939.0625, 5868325.0, 5878576.5625, 5884629.6875, 5909895.3125, 5910268.75, 5913950.0, 5920875.0, 5922079.6875, 5923128.125, 5924895.3125, 5926420.3125, 5931232.8125, 5931289.0625, 5931648.4375, 5937810.9375, 5950868.75, 5952243.75, 5956689.0625, 5975262.5, 6010259.375, 6034350.0, 6034367.1875, 6035717.1875, 6045203.125, 6045373.4375, 6107350.0, 6107462.5, 6116495.3125, 6117214.0625, 6118854.6875, 6119615.625, 6122517.1875, 6122528.125, 6123198.4375, 6137473.4375, 6140943.75, 6153803.125, 6154179.6875, 6164606.25, 6172214.0625, 6178732.8125, 6179531.25, 6181407.8125, 6205753.125, 6211875.0, 6213450.0, 6225665.625, 6231765.625, 6236398.4375, 6236400.0, 6236404.6875, 6260998.4375, 6270893.75, 6276895.3125, 6284298.4375, 6285343.75, 6290525.0, 6292803.125, 6294256.25, 6294520.3125, 6295495.3125, 6295501.5625, 6295828.125, 6296404.6875, 6296918.75, 6297087.5, 6298154.6875, 6298265.625, 6301901.5625, 6310104.6875, 6311875.0, 6328570.3125, 6328709.375, 6329470.3125, 6337890.625, 6338879.6875, 6348040.625, 6348187.5, 6354000.0, 6359185.9375, 6362120.3125, 6365779.6875, 6365871.875, 6369507.8125, 6373795.3125, 6374860.9375, 6374959.375, 6375948.4375, 6381031.25, 6384264.0625, 6384392.1875, 6385601.5625, 6386307.8125, 6386879.6875, 6390520.3125, 6395662.5, 6396964.0625, 6399181.25, 6401507.8125, 6407946.875, 6408001.5625, 6410643.75, 6410648.4375, 6411650.0, 6412393.75, 6412451.5625, 6416378.125, 6416635.9375, 6418126.5625, 6419359.375, 6419471.875, 6422156.25, 6426595.3125, 6431257.8125, 6431275.0, 6443267.1875, 6444823.4375, 6445287.5, 6450365.625, 6453176.5625, 6453195.3125, 6453445.3125, 6453525.0, 6453953.125, 6454562.5, 6455629.6875, 6456620.3125, 6462695.3125, 6462704.6875, 6464318.75, 6465598.4375, 6469204.6875, 6473532.8125, 6475532.8125, 6477868.75, 6477882.8125, 6479181.25, 6481104.6875, 6481543.75, 6482462.5, 6483679.6875, 6485120.3125, 6490718.75, 6490878.125, 6493290.625, 6495804.6875, 6503553.125, 6510087.5, 6512015.625, 6512403.125, 6513850.0, 6514632.8125, 6516128.125, 6518473.4375, 6521920.3125, 6522589.0625, 6524535.9375, 6528795.3125, 6529454.6875, 6531757.8125, 6532001.5625, 6532228.125, 6533046.875, 6533890.625, 6533993.75, 6534092.1875, 6534104.6875, 6534398.4375, 6534485.9375, 6534540.625, 6534870.3125, 6535714.0625, 6536657.8125, 6537370.3125, 6538042.1875, 6538471.875, 6538646.875, 6538751.5625, 6538987.5, 6541532.8125, 6541592.1875, 6541934.375, 6542043.75, 6543900.0, 6543946.875, 6548914.0625, 6550000.0, 6555868.75, 6556376.5625, 6556914.0625, 6559601.5625, 6559757.8125, 6562078.125, 6562137.5, 6562659.375, 6563259.375, 6563328.125, 6563348.4375, 6563589.0625, 6563610.9375, 6564532.8125, 6565235.9375, 6570531.25, 6571701.5625, 6572684.375, 6573951.5625, 6573985.9375, 6575054.6875, 6575329.6875, 6575429.6875, 6576720.3125, 6577021.875, 6578181.25, 6578235.9375, 6579962.5, 6580465.625, 6580521.875, 6580626.5625, 6582337.5, 6583095.3125, 6586312.5, 6587685.9375, 6587742.1875, 6590193.75, 6590210.9375, 6590318.75, 6591201.5625, 6593070.3125, 6594659.375, 6594965.625, 6597537.5, 6599832.8125, 6600621.875, 6600692.1875, 6601557.8125, 6601810.9375, 6602834.375, 6605082.8125, 6605095.3125, 6605914.0625, 6607081.25, 6607164.0625, 6607885.9375, 6609381.25, 6609635.9375, 6617843.75, 6622465.625, 6626889.0625, 6628828.125, 6631379.6875, 6633635.9375, 6634189.0625, 6637800.0, 6638217.1875, 6639678.125, 6648617.1875, 6648714.0625, 6649259.375, 6649587.5, 6654384.375, 6656885.9375, 6675267.1875, 6676045.3125, 6679459.375, 6688209.375, 6688637.5, 6688857.8125, 6692514.0625, 6694685.9375, 6700328.125, 6700721.875, 6701975.0, 6703598.4375, 6707439.0625, 6716284.375, 6717662.5, 6722821.875, 6726289.0625, 6727328.125, 6729115.625, 6729201.5625, 6733485.9375, 6733845.3125, 6734176.5625, 6734270.3125, 6735014.0625, 6737450.0, 6737739.0625, 6737745.3125, 6738331.25, 6739204.6875, 6740117.1875, 6743189.0625, 6744381.25, 6746832.8125, 6748928.125, 6758453.125, 6765581.25, 6766345.3125, 6767529.6875, 6769968.75, 6773750.0, 6775525.0, 6778418.75, 6785853.125, 6788176.5625, 6789664.0625, 6791737.5, 6797053.125, 6798537.5, 6802812.5, 6806085.9375, 6807973.4375, 6808314.0625, 6809484.375, 6810926.5625, 6810946.875, 6817709.375, 6820435.9375, 6821660.9375, 6824609.375, 6824776.5625, 6827606.25, 6827667.1875, 6830176.5625, 6830321.875, 6832028.125, 6832303.125, 6832750.0, 6832782.8125, 6833887.5, 6834196.875, 6834754.6875, 6835962.5, 6837501.5625, 6839189.0625, 6839371.875, 6840176.5625, 6841854.6875, 6841890.625, 6841989.0625, 6842385.9375, 6843281.25, 6843492.1875, 6844131.25, 6844210.9375, 6844256.25, 6845025.0, 6845179.6875, 6846879.6875, 6853315.625, 6877454.6875, 6877467.1875, 6879728.125, 6880060.9375, 6883640.625, 6885059.375, 6885073.4375, 6891931.25, 6892937.5, 6893018.75, 6895690.625, 6897728.125, 6918784.375, 6918787.5, 6920029.6875, 6921639.0625, 6922851.5625, 6925350.0, 6925621.875, 6927460.9375, 6932770.3125, 6934234.375, 6940609.375, 6941910.9375, 6941950.0, 6942523.4375, 6943695.3125, 6945593.75, 6945820.3125, 6947032.8125, 6947296.875, 6952042.1875, 6960453.125, 6977898.4375, 6980271.875, 6981100.0, 6981117.1875, 6981534.375, 6984995.3125, 6986426.5625, 6986484.375, 6986551.5625, 6987107.8125, 6987556.25, 6987846.875, 6987859.375, 6988218.75, 6988842.1875, 6988959.375, 6988985.9375, 6989051.5625, 6989057.8125, 6989379.6875, 6989409.375, 6990923.4375, 6991404.6875, 6991453.125, 6991573.4375, 6991942.1875, 6992051.5625, 6993012.5, 6993681.25, 6993881.25, 6994371.875, 6994445.3125, 6994925.0, 6995023.4375, 6995107.8125, 6995309.375, 6995473.4375, 6996628.125, 6996981.25, 6997337.5, 6997748.4375, 6997798.4375, 6997889.0625, 6998454.6875, 6998489.0625, 6998540.625, 6998542.1875, 6999332.8125, 6999451.5625, 6999521.875, 7001189.0625, 7001204.6875, 7001207.8125, 7001390.625, 7003003.125, 7003445.3125, 7005314.0625, 7005335.9375, 7009192.1875, 7016895.3125, 7020592.1875, 7022682.8125, 7032073.4375, 7032923.4375, 7035607.8125, 7040003.125, 7042409.375, 7043978.125, 7045521.875, 7045550.0, 7045654.6875, 7047379.6875, 7051001.5625, 7052907.8125, 7063760.9375, 7080603.125, 7088414.0625, 7097507.8125, 7111150.0, 7118500.0, 7122682.8125, 7123173.4375, 7123182.8125, 7129781.25, 7130006.25, 7130442.1875, 7150948.4375, 7151957.8125, 7152945.3125, 7153857.8125, 7157454.6875, 7158378.125, 7158442.1875, 7160437.5, 7167154.6875, 7167593.75, 7168675.0, 7169735.9375, 7169871.875, 7170834.375, 7176632.8125, 7180800.0, 7180803.125, 7186521.875, 7186787.5, 7190051.5625, 7191351.5625, 7192045.3125, 7192112.5, 7196346.875, 7199087.5, 7199139.0625, 7199753.125, 7200848.4375, 7201506.25, 7201720.3125, 7202951.5625, 7202957.8125, 7202971.875, 7202975.0, 7203603.125, 7203700.0, 7203759.375, 7204253.125, 7206335.9375, 7206928.125, 7209689.0625, 7222415.625, 7224040.625, 7228606.25, 7229360.9375, 7230101.5625, 7232987.5, 7235026.5625, 7236795.3125, 7255696.875, 7257245.3125, 7261845.3125, 7262132.8125, 7264621.875, 7267503.125, 7268381.25, 7269864.0625, 7271757.8125, 7278056.25, 7282176.5625, 7286479.6875, 7301062.5, 7306542.1875, 7308003.125, 7312517.1875, 7313971.875, 7318020.3125, 7327485.9375, 7332106.25, 7344387.5, 7353012.5, 7355571.875, 7370806.25, 7374682.8125, 7378696.875, 7378742.1875, 7380267.1875, 7384067.1875, 7384193.75, 7387600.0, 7391070.3125, 7394920.3125, 7396345.3125, 7397000.0, 7398423.4375, 7399346.875, 7399812.5, 7399842.1875, 7399860.9375, 7400406.25, 7400467.1875, 7401298.4375, 7401298.4375, 7401443.75, 7401820.3125, 7401942.1875, 7402331.25, 7402517.1875, 7402585.9375, 7403037.5, 7403042.1875, 7403589.0625, 7403900.0, 7404025.0, 7404831.25, 7405089.0625, 7405331.25, 7405571.875, 7406167.1875, 7406365.625, 7406421.875, 7406443.75, 7406571.875, 7406625.0, 7406628.125, 7406768.75, 7406898.4375, 7406985.9375, 7407554.6875, 7407560.9375, 7409626.5625, 7409710.9375, 7410323.4375, 7410696.875, 7410720.3125, 7410920.3125, 7411446.875, 7412246.875, 7412610.9375, 7413185.9375, 7413284.375, 7413585.9375, 7415032.8125, 7415071.875, 7418085.9375, 7418142.1875, 7418146.875, 7418939.0625, 7420471.875, 7420496.875, 7421318.75, 7421784.375, 7422079.6875, 7422989.0625, 7423114.0625, 7423535.9375, 7423620.3125, 7423739.0625, 7424296.875, 7424368.75, 7424629.6875, 7424981.25, 7425031.25, 7425132.8125, 7425151.5625, 7425226.5625, 7425331.25, 7425989.0625, 7427351.5625, 7428912.5, 7429190.625, 7430017.1875, 7430881.25, 7430885.9375, 7430940.625, 7431781.25, 7431784.375, 7431870.3125, 7432498.4375, 7433459.375, 7433464.0625, 7437148.4375, 7438045.3125, 7438800.0, 7439235.9375, 7439407.8125, 7442600.0, 7444151.5625, 7451131.25, 7452737.5, 7453160.9375, 7453285.9375, 7453464.0625, 7454062.5, 7454131.25, 7455009.375, 7455282.8125, 7455856.25, 7456696.875, 7457203.125, 7459103.125, 7459659.375, 7459793.75, 7463287.5, 7464278.125, 7464387.5, 7466350.0, 7469775.0, 7469829.6875, 7470676.5625, 7471781.25, 7479200.0, 7479704.6875, 7481104.6875, 7482906.25, 7488018.75, 7490910.9375, 7492776.5625, 7504732.8125, 7506354.6875, 7506360.9375, 7507032.8125, 7509010.9375, 7517221.875, 7520279.6875, 7539343.75, 7542398.4375, 7542471.875, 7542475.0, 7542571.875, 7542679.6875, 7542687.5, 7542807.8125, 7543168.75, 7543812.5, 7543864.0625, 7545325.0, 7545520.3125, 7545715.625, 7546445.3125, 7546448.4375, 7546454.6875, 7547245.3125, 7547459.375, 7547767.1875, 7548000.0, 7548103.125, 7548165.625, 7548456.25, 7548637.5, 7548821.875, 7548914.0625, 7549021.875, 7549201.5625, 7549326.5625, 7549345.3125, 7549782.8125, 7550331.25, 7550503.125, 7550539.0625, 7551343.75, 7551512.5, 7552079.6875, 7552206.25, 7552321.875, 7552971.875, 7553295.3125, 7554565.625, 7554662.5, 7554762.5, 7554907.8125, 7554945.3125, 7555431.25, 7556512.5, 7558326.5625, 7559665.625, 7570739.0625, 7571868.75, 7575484.375, 7575648.4375, 7575720.3125, 7575804.6875, 7577882.8125, 7577887.5, 7579265.625, 7579795.3125, 7579984.375, 7580779.6875, 7581193.75, 7581590.625, 7581696.875, 7582870.3125, 7583057.8125, 7584373.4375, 7584904.6875, 7587115.625, 7594201.5625, 7598643.75, 7602770.3125, 7603598.4375, 7605089.0625, 7605134.375, 7605914.0625, 7606571.875, 7607732.8125, 7607981.25, 7616067.1875, 7623285.9375, 7625106.25, 7626532.8125, 7628393.75, 7630475.0, 7631223.4375, 7631621.875, 7649718.75, 7658464.0625, 7664843.75, 7677440.625, 7679860.9375, 7682253.125, 7685754.6875, 7696390.625, 7697290.625, 7700509.375, 7700785.9375, 7700865.625, 7706804.6875, 7712309.375, 7713662.5, 7719706.25, 7720771.875, 7757917.1875, 7759796.875, 7760037.5, 7760629.6875, 7762906.25, 7763150.0, 7763242.1875, 7763879.6875, 7764178.125, 7764412.5, 7764415.625, 7764895.3125, 7766042.1875, 7766298.4375, 7767123.4375, 7767359.375, 7767546.875, 7767732.8125, 7767743.75, 7767879.6875, 7768082.8125, 7768085.9375, 7768504.6875, 7768779.6875, 7769778.125, 7769887.5, 7769901.5625, 7769945.3125, 7769968.75, 7770259.375, 7770329.6875, 7770342.1875, 7770381.25, 7770573.4375, 7770579.6875, 7770607.8125, 7770623.4375, 7770917.1875, 7770925.0, 7771493.75, 7771523.4375, 7771650.0, 7771703.125, 7772079.6875, 7776275.0, 7781475.0, 7781792.1875, 7783450.0, 7787321.875, 7787773.4375, 7788987.5, 7789923.4375, 7790475.0, 7792484.375, 7792971.875, 7799373.4375, 7801475.0, 7802045.3125, 7802546.875, 7802665.625, 7803865.625, 7803998.4375, 7805984.375, 7808557.8125, 7813931.25, 7828271.875, 7862162.5, 7863454.6875, 7864765.625, 7865290.625, 7866617.1875, 7866756.25, 7866757.8125, 7868154.6875, 7868167.1875, 7868465.625, 7868681.25, 7869381.25, 7869670.3125, 7869812.5, 7869987.5, 7870084.375, 7870825.0, 7870915.625, 7870923.4375, 7871015.625, 7872685.9375, 7881601.5625, 7883218.75, 7884246.875, 7892050.0, 7893835.9375, 7893890.625, 7894515.625, 7894953.125, 7895953.125, 7899523.4375, 7899595.3125, 7900356.25, 7902392.1875, 7903342.1875, 7918568.75, 7918996.875, 7919084.375, 7919198.4375, 7929390.625, 7930764.0625, 7939103.125, 7939123.4375, 7941054.6875, 7942542.1875, 7942778.125, 7944251.5625, 7944532.8125, 7945579.6875, 7946056.25, 7947525.0, 7947751.5625, 7952901.5625, 7952968.75, 7953270.3125, 7954014.0625, 7959406.25, 7963553.125, 7975454.6875, 7975643.75, 7990943.75, 7991581.25, 7992701.5625, 7992789.0625, 7993678.125, 7993996.875, 7995153.125, 7995551.5625, 7997607.8125, 7997617.1875, 7997660.9375, 7998934.375, 8000775.0, 8002462.5, 8006862.5, 8033621.875, 8034043.75, 8039560.9375, 8039589.0625, 8042185.9375, 8048268.75, 8048740.625, ...], [45.39907160677232, 58.43472597827107, 54.13986925637357, 89.12066462801411, 42.31772412446346, 44.334811280756355, 5.039176232203347, 5.998134463936753, 45.86931391332472, 84.75723619386078, 17.776354316266083, 13.483571108427403, 26.221404666351518, 122.39943040280984, 6.617917723442894, 6.467919871230439, 5.631474168321943, 39.62706428287824, 15.071489453073063, 84.36847826439555, 22.102515067820722, 6.216946875875842, 9.868371161867207, 10.322389149980804, 62.79481252107678, 93.3164113725934, 17.54321205886551, 12.152756737713236, 5.093415425987701, 6.376425591099246, 8.405526394965008, 74.97857949686866, 44.42105796148373, 29.4333807339959, 7.44665754383101, 29.082651627926793, 36.500974862091454, 12.555573518695828, 10.942724627448971, 7.588849913883986, 19.265618039863664, 9.809752109315776, 10.8106488535208, 80.32804657115567, 39.30421997334413, 16.6046771494948, 10.846970809067779, 16.143030242617993, 9.621170738191072, 6.45128784867267, 29.07133178912585, 10.544364280910047, 10.960144356161543, 59.14702833810563, 9.847484428910874, 41.4362725105399, 19.777389623817008, 26.1335112794383, 6.59164146767698, 13.939314326181142, 95.82585115343834, 23.38383925029865, 48.384512544701074, 78.15949706271037, 48.56458390501364, 5.475348409468644, 20.003753260607663, 26.05833028013529, 22.2631641647619, 59.94786198705246, 19.383723204735603, 51.56502116682048, 37.29317231216768, 6.1591816033042175, 181.55652163621295, 66.91443961230476, 31.526858619254003, 75.95028961562377, 65.21755936792128, 16.49099155461073, 55.00423227520606, 58.93448126855243, 52.57303396798196, 5.217553170091159, 6.297038156735923, 55.762434798435486, 122.59244627234175, 5.567883196850918, 36.64121541488077, 9.123455006903669, 6.155750364387428, 115.11362335820226, 14.155945289610052, 13.896946396998004, 109.52134482320002, 6.521761124976892, 28.894246811617883, 17.701096686631846, 6.841739455429855, 115.98313440439134, 87.99904748472544, 104.12694747318481, 38.645243162261444, 7.5688279260069775, 18.211069104169813, 122.5099663949049, 56.258755716947455, 42.31201243482917, 85.38430680996511, 28.627017323532996, 78.52249036185356, 22.257552213644875, 8.337888736273742, 10.884909902192174, 73.42403188280775, 18.94452891252995, 64.69291889600817, 53.75845680785268, 200.4307381448822, 15.403083707772161, 42.345565772351115, 10.792066111075812, 48.952598545524296, 82.04082184310545, 17.459621910125318, 22.578644341523503, 32.518121959953405, 18.596584589046813, 23.562605590381438, 66.2884170327356, 7.150933068576006, 74.30785344515007, 10.707748350875635, 59.1317663370279, 8.145896436792253, 39.19595319609624, 22.69215664318009, 6.756174477758828, 86.82839859316316, 23.55519733864917, 24.66577642177038, 8.245412884681722, 41.888715122311275, 5.177433677294963, 19.40547904648074, 5.395565708163598, 7.059264146471218, 5.156112256451142, 24.265766943295723, 25.07708309080054, 45.60033010662322, 17.213412925676487, 88.4263384452878, 62.19277057868756, 18.004235835022033, 9.285092982644406, 37.843745776663624, 21.201363741903585, 6.189170000128272, 17.632239923206257, 12.158797066862332, 5.387424293221695, 19.956927791612976, 31.544750338157037, 19.833414538623572, 65.70790924322222, 5.9092602237596985, 9.053405579304599, 88.24950473840443, 7.046788741057286, 52.259797094185, 32.91496672789286, 5.5568869593321235, 31.668638883872763, 14.004768803728018, 12.447933234196881, 91.88842608447418, 50.395568154099, 53.13890747925048, 12.856285612819297, 89.23664603720526, 58.09053308802784, 21.94958266811906, 17.63833934598756, 65.29770703729992, 69.71905195484634, 30.90831948911661, 15.933674677695823, 22.296600553194374, 60.34363931774433, 10.461005516721935, 5.0454657469302715, 50.0406825842092, 24.414766266230497, 9.893918438928216, 15.744856432955881, 10.643970707191222, 21.10846471513952, 136.64393338770643, 6.412187549294856, 22.571910458958776, 81.54459522258487, 16.673840323265726, 6.0985001814627715, 30.4960606464605, 163.44339035452413, 18.918455278445634, 44.91867856895598, 27.909485983714127, 29.521661742604593, 23.275619798221083, 57.28065002286976, 29.385702597457005, 7.4102310088436685, 102.6445590302597, 12.474542145873293, 21.56091987854063, 118.42552625068336, 153.20306438310783, 42.855443346989965, 96.45129324467587, 48.67620800629379, 53.447686966167716, 75.11213513887965, 45.31277737228109, 71.88347432463344, 23.646482817622456, 26.69004483286118, 160.01336347064196, 60.671913085868574, 6.388936207090076, 86.339830885314, 68.67242066656685, 20.332071743757655, 192.90851180289417, 13.50086638863289, 18.73622133159905, 61.92237559855329, 25.657698742525973, 57.16394346391178, 7.231209516791493, 94.66997077007252, 65.28848137767659, 115.5350612809463, 105.97551059442645, 85.70309810791288, 76.99501375011621, 21.532455954106926, 186.96168205187627, 84.84958542742243, 19.59833539313485, 130.41424353371184, 20.606194671347275, 44.53396692833216, 5.147449868012342, 27.402908760491336, 41.23654180808661, 6.320732704602066, 5.066039091525922, 27.528909248391273, 11.940316350267715, 5.64966048804249, 15.34483684382827, 68.79196299466794, 49.4414407200757, 25.297074013286498, 25.80380395072883, 49.925189759728056, 203.41652578947088, 41.37034528561626, 19.413489063201066, 22.740068262720936, 69.91825689622371, 6.491557075161738, 5.245993923447846, 40.1234117878206, 61.13955880638159, 105.04512716729712, 24.02470375118437, 95.43278300318732, 38.86067762456453, 52.8405476598883, 5.597005252311804, 5.1062186712719075, 120.85197042249324, 37.04775804224042, 148.4827380106409, 6.054074193955454, 65.56498508992377, 60.06850119040066, 66.09593712981703, 6.289839946194938, 97.59979343563244, 67.75257948420965, 12.55574537977661, 19.579557538956394, 42.243455950252496, 93.94876955573942, 10.214295549138312, 86.26770136285887, 11.239175216115399, 93.89453134919016, 48.632434070424566, 56.31239966064271, 16.974253018025315, 13.284059282926744, 55.94175496813692, 45.328864203068946, 33.46019774074443, 33.67656916228863, 166.85197878594138, 35.078030617766686, 20.302544005692127, 13.688339736074413, 15.04452041373881, 173.10406082439712, 119.04374525554634, 5.536420164809132, 28.910807613527968, 11.798533416214562, 71.23856367948319, 36.56810007910176, 5.183970751278581, 6.271037943667835, 16.58279971874704, 99.16234110249987, 5.29110979192944, 63.82699222193669, 165.34392507816844, 56.558922170957445, 8.417177091023483, 60.70816415202323, 6.358755406015505, 10.224438558589181, 126.75165885212118, 33.45551341054375, 87.51603664614346, 41.34328960950465, 38.45863954065044, 30.71793378594299, 52.573576853859514, 12.405968439363845, 9.936772881970768, 11.516754586909748, 53.208626157814926, 5.1254533998877525, 33.38268387540789, 56.57641282180743, 25.377644201192204, 31.566780382558477, 54.60241753769744, 164.83431835179906, 6.041046825041573, 18.417918709146356, 87.33145765011676, 28.277067543023254, 50.450172829026606, 6.138004097805336, 10.930642203226215, 20.977485891403944, 21.271989860957593, 25.08747440137797, 46.577723697459646, 15.287785020210238, 15.68805122087268, 92.15269797873349, 71.33130413015576, 18.96887726153514, 91.48311912760275, 36.46875108484006, 38.0887057852149, 21.824647857830204, 30.153747873777245, 44.71174270140521, 100.81145934145079, 102.36991611454101, 11.274615133291446, 5.556098077682061, 11.931118094582532, 6.526406230879605, 81.51127708683586, 63.55272371555312, 36.29085351670314, 140.76784046351992, 5.469384465884944, 29.803411402502242, 43.557315029470246, 45.30427071432581, 8.089787249377078, 73.45350093798532, 31.186414730630744, 76.10207636595962, 51.928865988721824, 36.146434187043766, 72.97701412916159, 11.08585042027072, 40.21400977253449, 14.2150281952447, 14.058412170920274, 11.115480864398485, 16.201867928201565, 69.17968872420569, 11.281210534829352, 67.65173082092512, 49.242126947805446, 43.83393473092084, 7.7995722527754, 152.52156616805956, 19.984556308480105, 77.18302296339232, 8.048554268843452, 23.02039655196422, 128.0673188399014, 50.196296468777554, 9.34252793461419, 47.299820719892935, 18.782294796546264, 133.53990843514185, 67.38985519011622, 58.03749151043067, 21.14784692288888, 50.73945281382011, 13.460610761028006, 49.68708529325806, 35.935099630540044, 29.221446354052762, 34.65648861643882, 16.469711879910875, 97.20530485075014, 84.1989207251092, 25.445292421066412, 74.19879793961121, 31.65260505276899, 118.60926307237679, 6.201124605319172, 55.74620815433049, 21.715132998828715, 6.208759179210996, 11.716829216942017, 9.12402337275018, 22.578604468202663, 36.38174943106439, 59.61004846525074, 13.47185660808665, 87.03466393497298, 11.6279752627118, 5.350525528549183, 7.239744481714719, 51.47016199737295, 10.47602973334048, 17.09946966572913, 123.54978959282262, 35.51582230486843, 31.024608105432446, 8.78835583754656, 33.85602790624039, 51.527600122273405, 6.439514952237363, 15.531036282208104, 32.40325296215359, 20.829721709435912, 68.92696691691638, 5.171471823328056, 22.822942443731623, 61.68168187061277, 12.159356070165803, 60.373039087197085, 108.03092352171295, 146.72266425352174, 58.78546962661646, 42.316889095073556, 50.37956971798566, 38.51733141671543, 14.262183383064622, 5.050412223802015, 10.391315907239585, 92.98133332472621, 109.01681276723616, 107.57677757156662, 65.17331094454556, 38.45714045876865, 26.489926954124385, 75.78606146481745, 62.62642188514866, 107.75808847282715, 47.201472066537214, 15.224220158645599, 88.61261378065052, 5.246507166208237, 9.196183800333909, 25.5433698764445, 74.38245817399869, 72.59973650038027, 77.72701736203068, 87.41163007485763, 35.9562353589701, 52.08143432446701, 109.53785943024516, 13.030440405708367, 17.487146313041478, 41.042352685427005, 97.87246164702474, 62.34486879837736, 38.31870887799064, 20.280999746874148, 85.14945489347565, 7.635625081324512, 34.161714483319514, 22.400516908770786, 52.82574428391334, 65.77786237932874, 156.86033516204452, 72.79349761874583, 142.3151378353625, 233.5330024470549, 5.2203407063935074, 123.7779831063507, 33.88265626281833, 23.46587028439319, 36.80403613385671, 8.606679572351961, 97.66796020796627, 87.40365545968612, 18.464249177237146, 88.64058987993567, 24.724413704303355, 32.01993004030046, 12.135625574283186, 111.03745894821475, 28.560903146677177, 60.71702151980888, 58.63033881120785, 22.9922407513341, 21.695210465384772, 16.84382828732733, 34.46430373014289, 72.2521372167779, 91.41723436965684, 211.40120643536588, 43.2932570111605, 48.99300365345403, 95.91261646263202, 9.177647166185327, 128.49438608132056, 60.1128605065641, 19.970274137891288, 304.59736704998477, 12.167539500027882, 21.75076481220733, 19.406917777809213, 79.89871699143055, 57.34471705950752, 45.917912871055954, 80.7559504153946, 185.39896762273395, 13.748685547157848, 20.276649402702382, 122.71099181489973, 40.816159299512805, 57.61074094123394, 16.470648868376493, 19.71684519728231, 10.767044093380509, 47.6203582457846, 14.555605887664182, 77.06464152766388, 12.696578649140196, 77.73410738845169, 8.325066835398632, 72.35507476471018, 22.715771897096168, 162.58787937064534, 19.180345691539177, 58.073722015515756, 11.951487697543904, 63.381894735702424, 7.548881960412486, 13.14133130234406, 68.71958501549082, 53.83584026335273, 8.473532988363099, 22.326460173759187, 5.705549334331559, 78.17408898892896, 25.468255842330855, 10.71092731355036, 13.632211425926862, 8.568961975309552, 180.76506336449532, 83.53662257697995, 109.51898888462628, 7.270365270042803, 18.002409288274954, 121.28398442527717, 48.42043038639474, 25.80537217468133, 158.72491843572791, 54.77437089312585, 75.68020319506877, 9.183297854807824, 51.80007553063447, 8.204838522401333, 37.36375274347866, 29.08008488932845, 11.119943796007474, 211.37403816321228, 37.49607251184844, 27.29079211988199, 16.019604834639615, 16.81871168730808, 162.24795618616378, 348.0208824018029, 11.512538419311905, 84.95693028645478, 62.65182410166494, 61.814011075056726, 7.05996500998451, 21.846587849008756, 86.28542315174298, 78.0221493605246, 9.00742927464562, 21.360821169463694, 18.402261284517184, 39.446510400864625, 16.429311389551888, 43.05576062745309, 39.23774738371048, 17.43963231424453, 13.926713281780472, 27.72959521880666, 113.64483357506909, 159.2781336004997, 7.976879425568566, 92.74080512230546, 246.9537892285316, 19.58411955664895, 116.95838285876593, 16.608827621611027, 23.82849815299224, 25.912249373524652, 17.482252030058426, 31.653362635410506, 25.01550224285874, 20.95775182040264, 101.06927365944459, 30.617737634293213, 9.04249233137478, 106.335675908925, 73.53989227980753, 6.010301376745381, 12.649533435369984, 146.9303231401622, 110.04507905263563, 43.677188294266, 29.286999113927482, 16.213836664877828, 33.240068874474346, 6.18920465537694, 19.130284831495686, 136.2562310080708, 126.7274729281073, 22.95495896703929, 7.004520810990987, 140.78092181156836, 28.085399735984392, 45.43943228508579, 5.648259787593622, 69.62042488994938, 47.739556490555444, 114.19971017845471, 71.53522721889205, 44.46902136650733, 9.425323360214371, 57.9810514966439, 114.66597044116139, 41.12324250059271, 23.175964011569718, 79.5861969635528, 53.0532927746423, 62.727365927343456, 109.94989438830285, 15.658905834381851, 21.778542632880804, 36.36632946188011, 10.222091554119853, 6.7873618683553, 28.880471070308236, 131.67131667382372, 52.619429656445, 28.06603059178011, 157.21840007160353, 50.803075046710006, 19.95689149122119, 21.988561441577644, 12.551705730635698, 5.436004994134631, 48.716992764783406, 13.180159338230013, 156.03150992438563, 123.65860071903353, 58.321931811172945, 5.922546506948526, 13.924990471995743, 38.9451241884922, 6.259461496906064, 30.37976501613859, 34.78347690602414, 33.518193600280554, 5.772827310837642, 88.51492599116492, 9.565603052899677, 11.40817185140431, 12.52468523380964, 32.391526589132496, 59.77838631147118, 5.10018779349841, 14.397051408508585, 22.85011628409739, 106.19166152193783, 19.915839772133424, 44.988250154192315, 34.354753470506274, 69.4366667891856, 80.49454288386778, 65.04169189006475, 120.95514196717421, 611.4128352883533, 32.93152462355539, 23.068341812207123, 5.997137705745988, 34.75113389208721, 20.97452648198206, 20.03801649051698, 25.641992317704887, 95.16720337214147, 22.45960127705453, 437.88582059585286, 11.675223798831418, 58.15562503068316, 11.444024528109768, 6.937374818795225, 26.374443818529542, 73.81957967964568, 103.24728354463329, 20.168593988614965, 58.23224495807202, 34.39429696608186, 94.81145069717866, 21.742586468418516, 30.913765264350925, 10.226784405376707, 23.107715632397078, 20.5675082270712, 79.52161818918337, 90.48449923903593, 8.227089622147698, 6.952172153816562, 63.16990399274361, 16.044133783858697, 54.84976493540472, 8.812064684851697, 41.791993607787646, 6.144525099978994, 92.93738921275857, 15.908261090293172, 13.29835014140322, 9.59569761540489, 43.256685568109944, 7.024742463942281, 35.92521821170954, 14.576421070416197, 81.1560630541061, 33.79438771842482, 7.3684933456447, 24.431637367536588, 57.874373230253326, 67.57188236297732, 47.66323954987702, 90.8041140108025, 51.439682039556146, 51.731592215243325, 45.06556246986837, 10.254086518435326, 22.054003786842458, 8.081482865107146, 8.98044614607415, 66.80082676830175, 10.113797372636842, 6.9422301854809545, 83.7398824194271, 10.335155359737767, 22.02517719826753, 9.96656889405723, 56.579637635076196, 57.52168270548949, 22.0480743143437, 5.258909977878666, 144.47759392387496, 79.09108007978206, 69.24270188279641, 35.89004833355028, 60.26363541158082, 6.183840691371277, 31.97344266911736, 25.597662899510457, 67.9329868498281, 5.28261375676724, 20.734021386319046, 6.096347921549088, 67.14936939666518, 10.882703339478807, 5.169400099528197, 19.023574036903337, 9.303207850954688, 26.6600230183077, 21.860610439885672, 23.442237148723738, 127.29885380051086, 17.4648295190347, 42.68881929337756, 11.700428001468545, 29.921093116833646, 5.696499699919835, 19.2443317893341, 182.5170317848979, 29.351825875242476, 7.384217203800981, 158.02435935390545, 93.6796375627975, 49.734829366015376, 26.218616400446056, 70.48462904966239, 21.90022916982915, 5.884769877879415, 16.967991716379746, 22.625678562515667, 19.949479680015866, 23.960426107520373, 107.73796671430168, 45.33241198334421, 16.551612949877065, 44.9481301203236, 26.995265260177415, 40.769819419450286, 29.661874115851987, 38.128562972575686, 5.164584629287279, 25.23612596315981, 9.034669612732687, 105.35188320443694, 11.572188841072608, 62.41001496027766, 17.05613234304937, 15.67957286764548, 42.22906164301985, 41.40982766779813, 7.678404078917541, 105.6860288700245, 34.8983117569608, 13.12299864359369, 14.113481194007631, 87.43656831514963, 7.007152878075235, 76.97407925764693, 9.732250889280134, 26.929013880208228, 76.31602589461515, 77.79864999731775, 97.61357800441296, 21.620174771384054, 99.81059895227935, 71.19413335457321, 23.88995214378473, 21.27003325796731, 92.23502632421537, 11.054757102750854, 117.32421406600638, 24.871945528717053, 15.74115208184082, 107.90607263392289, 32.92274559824446, 5.141751749210712, 14.788473040988752, 55.17857078267191, 6.308394520941128, 6.6929168546123465, 40.176013510789836, 31.472872141328132, 161.9471794212171, 67.06244420388275, 61.39961217649998, 101.56502637727972, 46.73385749013866, 29.978349784088167, 55.30800203219855, 134.9125742941927, 24.704969679476264, 31.644982610625135, 25.42918615315442, 6.700093784072196, 18.612093252725106, 26.88599035105503, 49.864008982291196, 83.35506395713085, 67.23198892872176, 38.30576017927969, 23.82450356142916, 33.90428826276521, 18.172804751843756, 173.90643258183582, 101.22909867546974, 12.303051119172679, 18.976830191959877, 68.59231350290744, 153.606440785948, 85.63466852403086, 25.407117093849344, 64.1794818355354, 31.67049761975297, 23.752084294829487, 109.40375691974054, 246.93007837495756, 13.933803505673657, 36.39785808745186, 41.745722473090936, 34.62831860015093, 18.93191654541656, 17.59793771899404, 25.68149534796574, 43.347524211242, 74.88563508298893, 45.22936225302462, 66.62927528803118, 118.8088923752862, 8.722983434509487, 116.95347218458855, 41.768967289424516, 54.55988050220941, 20.278028923471513, 48.20969421488407, 32.31265642337569, 15.350711904886683, 85.67224521712531, 84.24951964685421, 5.161379431657956, 76.28379423604815, 33.1024970761187, 133.15944206051464, 91.09716623163561, 43.26042252568403, 11.712283702598027, 16.75597725330885, 33.45915741767293, 139.43996894880084, 19.158661069567664, 5.032455309033085, 29.58190433628818, 60.302959311488415, 21.451901482309, 7.383592000411834, 53.79639980374819, 125.50334648992018, 110.05042240838097, 13.80468840917342, 11.233641694111, 136.54139488580356, 30.58340829934098, 57.71066906076764, 11.252004605354923, 47.82266105891013, 26.46997064619461, 5.082945307235388, 122.16335575401973, 8.339207306507818, 11.987078167692127, 57.96118894037936, 50.94346112312344, 7.032659541177036, 60.41548987421088, 160.19421742491545, 42.77496277006378, 60.93280250334797, 10.970897776940634, 96.23301908766499, 63.046752910794055, 58.62633110801654, 9.827740970297738, 18.809408384890244, 20.119658655171765, 29.974775599673983, 64.8488788100482, 17.026645795901803, 44.96103701925489, 40.194098688191, 104.24190175524134, 30.897882771881072, 35.152632493071216, 85.48149784118453, 67.45284002160504, 9.141389828315232, ...])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)