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 = 44530
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)
<ipython-input-1-fe971fb7f7b0>:25: RuntimeWarning: invalid value encountered in scalar divide 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)))) <ipython-input-1-fe971fb7f7b0>:26: RuntimeWarning: invalid value encountered in scalar divide 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))))
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);
([2924346.875, 2944565.625, 2986623.4375, 3002956.25, 3147323.4375, 3397809.375, 3436142.1875, 3999750.0, 4792262.5, 5266640.625, 5416700.0, 5455485.9375, 5455493.75, 5456928.125, 5734228.125, 5736831.25, 5759906.25, 5788390.625, 5790696.875, 5795334.375, 5797781.25, 5839612.5, 5879600.0, 5883587.5, 5910326.5625, 5915831.25, 5917845.3125, 5924237.5, 5925865.625, 5926953.125, 5927317.1875, 5927843.75, 5929675.0, 5931595.3125, 5934132.8125, 5934760.9375, 5936929.6875, 5937835.9375, 5939276.5625, 5943960.9375, 5944914.0625, 5945410.9375, 5946307.8125, 5956275.0, 5956276.5625, 5958031.25, 5959417.1875, 5959895.3125, 5960125.0, 5960753.125, 5966514.0625, 5970970.3125, 5972035.9375, 5977026.5625, 5977026.5625, 5978792.1875, 5978793.75, 5979534.375, 5980356.25, 5980540.625, 5980615.625, 5980959.375, 5981001.5625, 5981025.0, 5981032.8125, 5981273.4375, 5981423.4375, 5981506.25, 5981570.3125, 5981609.375, 5981657.8125, 5981925.0, 5982148.4375, 5982321.875, 5982434.375, 5982484.375, 5982528.125, 5982728.125, 5983053.125, 5983329.6875, 5983426.5625, 5983896.875, 5983953.125, 5984093.75, 5984282.8125, 5984553.125, 5985385.9375, 5986021.875, 5986590.625, 5990492.1875, 5990950.0, 5990971.875, 6009821.875, 6015390.625, 6023539.0625, 6024298.4375, 6054309.375, 6055198.4375, 6058721.875, 6059100.0, 6063656.25, 6066542.1875, 6078581.25, 6084646.875, 6084657.8125, 6088356.25, 6091434.375, 6093745.3125, 6093850.0, 6096214.0625, 6098381.25, 6099982.8125, 6100262.5, 6124200.0, 6126929.6875, 6128434.375, 6128501.5625, 6128693.75, 6132237.5, 6132362.5, 6132404.6875, 6132954.6875, 6134525.0, 6138635.9375, 6142756.25, 6170331.25, 6175581.25, 6175920.3125, 6175925.0, 6176109.375, 6176148.4375, 6176251.5625, 6176275.0, 6177120.3125, 6177282.8125, 6177342.1875, 6177582.8125, 6177603.125, 6178493.75, 6179835.9375, 6180748.4375, 6180907.8125, 6181957.8125, 6182935.9375, 6183243.75, 6183395.3125, 6184270.3125, 6185128.125, 6186509.375, 6194025.0, 6209742.1875, 6210400.0, 6217031.25, 6218221.875, 6262400.0, 6267309.375, 6271317.1875, 6285198.4375, 6286095.3125, 6287734.375, 6291410.9375, 6292990.625, 6298095.3125, 6301600.0, 6303882.8125, 6306959.375, 6336956.25, 6338228.125, 6341228.125, 6342273.4375, 6343179.6875, 6343450.0, 6343943.75, 6344215.625, 6344232.8125, 6344571.875, 6344598.4375, 6344617.1875, 6344620.3125, 6344890.625, 6345023.4375, 6345162.5, 6345195.3125, 6345231.25, 6345573.4375, 6345728.125, 6345757.8125, 6345910.9375, 6346014.0625, 6346209.375, 6346900.0, 6347385.9375, 6347401.5625, 6347420.3125, 6347785.9375, 6348076.5625, 6348728.125, 6348964.0625, 6349156.25, 6349192.1875, 6349276.5625, 6349921.875, 6350020.3125, 6351696.875, 6352292.1875, 6352331.25, 6352590.625, 6352612.5, 6353323.4375, 6355182.8125, 6355465.625, 6355529.6875, 6356089.0625, 6356146.875, 6356576.5625, 6356970.3125, 6358621.875, 6359306.25, 6359323.4375, 6360632.8125, 6362714.0625, 6372925.0, 6373746.875, 6377125.0, 6387075.0, 6387101.5625, 6395089.0625, 6447707.8125, 6481442.1875, 6483025.0, 6484126.5625, 6485400.0, 6486600.0, 6487826.5625, 6495279.6875, 6496150.0, 6496159.375, 6496362.5, 6496417.1875, 6496687.5, 6497096.875, 6497334.375, 6497967.1875, 6498026.5625, 6498237.5, 6498279.6875, 6498581.25, 6498909.375, 6499729.6875, 6499948.4375, 6500839.0625, 6508262.5, 6512271.875, 6514573.4375, 6516168.75, 6516629.6875, 6518032.8125, 6523879.6875, 6527639.0625, 6529360.9375, 6541267.1875, 6542935.9375, 6544412.5, 6545184.375, 6547835.9375, 6548075.0, 6550118.75, 6550818.75, 6551314.0625, 6551706.25, 6551853.125, 6552904.6875, 6553110.9375, 6553231.25, 6553465.625, 6553667.1875, 6553776.5625, 6554157.8125, 6554596.875, 6554623.4375, 6554843.75, 6555184.375, 6555467.1875, 6556028.125, 6556085.9375, 6556203.125, 6556267.1875, 6556943.75, 6556965.625, 6558034.375, 6562751.5625, 6566706.25, 6567043.75, 6581642.1875, 6584039.0625, 6587303.125, 6588934.375, 6590175.0, 6590446.875, 6593926.5625, 6596121.875, 6598346.875, 6598392.1875, 6599023.4375, 6599153.125, 6599817.1875, 6600389.0625, 6600815.625, 6601582.8125, 6607121.875, 6616378.125, 6616556.25, 6616979.6875, 6617123.4375, 6618862.5, 6623275.0, 6642670.3125, 6647689.0625, 6652178.125, 6674725.0, 6683126.5625, 6695668.75, 6703696.875, 6704498.4375, 6719276.5625, 6719456.25, 6719482.8125, 6719523.4375, 6719896.875, 6720107.8125, 6722787.5, 6725245.3125, 6725745.3125, 6726210.9375, 6729340.625, 6730289.0625, 6733443.75, 6733454.6875, 6733928.125, 6734906.25, 6734996.875, 6735070.3125, 6736481.25, 6737221.875, 6737251.5625, 6737345.3125, 6737987.5, 6742279.6875, 6745264.0625, 6745907.8125, 6746176.5625, 6747793.75, 6747860.9375, 6749214.0625, 6750520.3125, 6751846.875, 6751926.5625, 6752326.5625, 6753300.0, 6753306.25, 6759739.0625, 6763539.0625, 6765523.4375, 6767343.75, 6768370.3125, 6769764.0625, 6769864.0625, 6770679.6875, 6776676.5625, 6778278.125, 6779393.75, 6779853.125, 6781718.75, 6781939.0625, 6782900.0, 6783450.0, 6784067.1875, 6785901.5625, 6785929.6875, 6787056.25, 6787090.625, 6787775.0, 6788421.875, 6788440.625, 6788642.1875, 6789384.375, 6789495.3125, 6789673.4375, 6790017.1875, 6791360.9375, 6791470.3125, 6791640.625, 6791792.1875, 6791828.125, 6791956.25, 6792179.6875, 6792260.9375, 6792781.25, 6793235.9375, 6793946.875, 6794045.3125, 6795615.625, 6809879.6875, 6811578.125, 6812414.0625, 6813021.875, 6815412.5, 6819867.1875, 6820987.5, 6821856.25, 6821895.3125, 6822035.9375, 6823390.625, 6824665.625, 6824876.5625, 6825195.3125, 6826151.5625, 6826357.8125, 6826454.6875, 6826839.0625, 6826850.0, 6827001.5625, 6827818.75, 6827881.25, 6828200.0, 6828453.125, 6828470.3125, 6829500.0, 6831670.3125, 6832375.0, 6833915.625, 6834368.75, 6836112.5, 6840454.6875, 6840464.0625, 6840509.375, 6840510.9375, 6840934.375, 6840937.5, 6841137.5, 6841140.625, 6841253.125, 6842260.9375, 6844643.75, 6847510.9375, 6848787.5, 6851209.375, 6852473.4375, 6853573.4375, 6853928.125, 6858198.4375, 6858671.875, 6859639.0625, 6859907.8125, 6862125.0, 6863965.625, 6864787.5, 6864960.9375, 6865084.375, 6865190.625, 6865206.25, 6865842.1875, 6868729.6875, 6869650.0, 6873051.5625, 6873943.75, 6874078.125, 6876496.875, 6876517.1875, 6878250.0, 6879478.125, 6879828.125, 6882109.375, 6883210.9375, 6893704.6875, 6903237.5, 6907682.8125, 6909642.1875, 6910242.1875, 6910378.125, 6910415.625, 6910529.6875, 6914932.8125, 6915995.3125, 6916112.5, 6916350.0, 6916507.8125, 6917054.6875, 6917259.375, 6917795.3125, 6918117.1875, 6918154.6875, 6918909.375, 6918934.375, 6918957.8125, 6919028.125, 6919071.875, 6919606.25, 6919785.9375, 6919873.4375, 6919948.4375, 6919978.125, 6920162.5, 6920246.875, 6920254.6875, 6920275.0, 6920364.0625, 6920423.4375, 6920693.75, 6920917.1875, 6920948.4375, 6921317.1875, 6921323.4375, 6921364.0625, 6921839.0625, 6921879.6875, 6921887.5, 6921920.3125, 6921942.1875, 6921964.0625, 6921968.75, 6922018.75, 6922120.3125, 6922253.125, 6922284.375, 6922292.1875, 6922354.6875, 6922448.4375, 6922482.8125, 6922507.8125, 6922600.0, 6922673.4375, 6922675.0, 6922767.1875, 6922793.75, 6922796.875, 6922907.8125, 6922909.375, 6922978.125, 6923031.25, 6923357.8125, 6923382.8125, 6923385.9375, 6923448.4375, 6923642.1875, 6923706.25, 6923723.4375, 6923778.125, 6923943.75, 6924182.8125, 6924273.4375, 6924296.875, 6924346.875, 6924639.0625, 6924685.9375, 6924710.9375, 6924765.625, 6925482.8125, 6925582.8125, 6925870.3125, 6926137.5, 6926567.1875, 6926606.25, 6926610.9375, 6926665.625, 6926670.3125, 6926801.5625, 6926859.375, 6926965.625, 6927110.9375, 6927850.0, 6928290.625, 6928543.75, 6928693.75, 6929810.9375, 6929990.625, 6931276.5625, 6931575.0, 6931692.1875, 6931715.625, 6931985.9375, 6932751.5625, 6933606.25, 6934704.6875, 6934714.0625, 6936109.375, 6936382.8125, 6937206.25, 6937353.125, 6937642.1875, 6937823.4375, 6938221.875, 6938401.5625, 6938945.3125, 6940976.5625, 6941468.75, 6942489.0625, 6943125.0, 6943398.4375, 6943618.75, 6943878.125, 6944743.75, 6945048.4375, 6945275.0, 6945312.5, 6945375.0, 6945657.8125, 6946314.0625, 6946690.625, 6946729.6875, 6948473.4375, 6950375.0, 6950990.625, 6951054.6875, 6953428.125, 6953579.6875, 6954956.25, 6955893.75, 6956270.3125, 6957575.0, 6958329.6875, 6958557.8125, 6958882.8125, 6961229.6875, 6961531.25, 6962248.4375, 6963715.625, 6966657.8125, 6972007.8125, 6972071.875, 6972503.125, 6972506.25, 6973884.375, 6976331.25, 6977882.8125, 6980239.0625, 6981967.1875, 6982017.1875, 6985701.5625, 6986907.8125, 6996946.875, 7003584.375, 7021931.25, 7022676.5625, 7024148.4375, 7024573.4375, 7024801.5625, 7026842.1875, 7027198.4375, 7027773.4375, 7028279.6875, 7029885.9375, 7030929.6875, 7031110.9375, 7040485.9375, 7041504.6875, 7043698.4375, 7043746.875, 7044601.5625, 7050225.0, 7050925.0, 7053329.6875, 7054293.75, 7054600.0, 7056348.4375, 7058498.4375, 7058601.5625, 7060473.4375, 7062290.625, 7063640.625, 7065537.5, 7066434.375, 7066515.625, 7067143.75, 7068031.25, 7068701.5625, 7068787.5, 7069237.5, 7069501.5625, 7069520.3125, 7069556.25, 7069771.875, 7072853.125, 7074782.8125, 7075415.625, 7076504.6875, 7077718.75, 7079539.0625, 7081946.875, 7082650.0, 7083671.875, 7083726.5625, 7084162.5, 7084609.375, 7087012.5, 7087178.125, 7087289.0625, 7108465.625, 7109739.0625, 7110051.5625, 7112256.25, 7118868.75, 7120237.5, 7120610.9375, 7121957.8125, 7122028.125, 7122396.875, 7126531.25, 7126537.5, 7127698.4375, 7127884.375, 7127981.25, 7128231.25, 7128260.9375, 7128593.75, 7131671.875, 7132656.25, 7134743.75, 7134743.75, 7138676.5625, 7142254.6875, 7144043.75, 7144675.0, 7144685.9375, 7145096.875, 7145753.125, 7146346.875, 7146564.0625, 7146718.75, 7146743.75, 7147043.75, 7147646.875, 7148246.875, 7148701.5625, 7149018.75, 7149393.75, 7150721.875, 7150729.6875, 7151068.75, 7151268.75, 7152187.5, 7152265.625, 7152637.5, 7152928.125, 7153235.9375, 7154014.0625, 7154770.3125, 7156054.6875, 7156275.0, 7156303.125, 7157048.4375, 7157450.0, 7157885.9375, 7158431.25, 7158657.8125, 7159223.4375, 7159434.375, 7159506.25, 7159804.6875, 7160067.1875, 7160132.8125, 7160151.5625, 7160165.625, 7160301.5625, 7160365.625, 7160467.1875, 7160656.25, 7160710.9375, 7160714.0625, 7160921.875, 7161073.4375, 7161112.5, 7161284.375, 7161420.3125, 7161517.1875, 7161712.5, 7162232.8125, 7162370.3125, 7162492.1875, 7162581.25, 7162696.875, 7162795.3125, 7162846.875, 7163475.0, 7165593.75, 7167217.1875, 7167607.8125, 7167757.8125, 7169803.125, 7170253.125, 7170739.0625, 7176346.875, 7176409.375, 7176529.6875, 7176837.5, 7178376.5625, 7179979.6875, 7188281.25, 7192954.6875, 7196901.5625, 7205226.5625, 7206895.3125, 7208975.0, 7209975.0, 7210984.375, 7211687.5, 7212275.0, 7213364.0625, 7213571.875, 7214796.875, 7215300.0, 7216334.375, 7216740.625, 7216760.9375, 7216767.1875, 7217109.375, 7218109.375, 7218309.375, 7218712.5, 7219104.6875, 7219184.375, 7221693.75, 7221932.8125, 7222025.0, 7222228.125, 7222354.6875, 7222914.0625, 7223493.75, 7223623.4375, 7223759.375, 7224009.375, 7224017.1875, 7224107.8125, 7224135.9375, 7224290.625, 7225051.5625, 7225400.0, 7226271.875, 7226517.1875, 7226937.5, 7227267.1875, 7227476.5625, 7228623.4375, 7228885.9375, 7229714.0625, 7230523.4375, 7231045.3125, 7231375.0, 7231400.0, 7231829.6875, 7231848.4375, 7232993.75, 7233007.8125, 7233868.75, 7233898.4375, 7234257.8125, 7234362.5, 7235712.5, 7237325.0, 7238143.75, 7238920.3125, 7245321.875, 7247782.8125, 7249503.125, 7249884.375, 7250350.0, 7250675.0, 7250742.1875, 7250753.125, 7253918.75, 7254842.1875, 7258203.125, 7259095.3125, 7260581.25, 7260595.3125, 7263329.6875, 7264870.3125, 7267870.3125, 7270634.375, 7270667.1875, 7271054.6875, 7275165.625, 7287562.5, 7295409.375, 7295478.125, 7307853.125, 7312629.6875, 7314253.125, 7315639.0625, 7316850.0, 7316942.1875, 7317525.0, 7319484.375, 7320073.4375, 7322621.875, 7322870.3125, 7322926.5625, 7324837.5, 7325293.75, 7326329.6875, 7326343.75, 7326862.5, 7327868.75, 7329995.3125, 7330507.8125, 7330781.25, 7330864.0625, 7330965.625, 7330970.3125, 7331039.0625, 7331103.125, 7331132.8125, 7331268.75, 7331450.0, 7331490.625, 7331496.875, 7331556.25, 7331595.3125, 7331656.25, 7331831.25, 7331931.25, 7332026.5625, 7332139.0625, 7332150.0, 7332170.3125, 7332189.0625, 7332245.3125, 7332296.875, 7332376.5625, 7332570.3125, 7332795.3125, 7333228.125, 7333245.3125, 7333573.4375, 7349715.625, 7404610.9375, 7422159.375, 7423446.875, 7424810.9375, 7432853.125, 7433401.5625, 7444418.75, 7458937.5, 7459159.375, 7499614.0625, 7507142.1875, 7514735.9375, 7516075.0, 7526598.4375, 7530245.3125, 7533860.9375, 7537785.9375, 7539229.6875, 7546154.6875, 7547800.0, 7547820.3125, 7548732.8125, 7552317.1875, 7558478.125, 7558610.9375, 7563135.9375, 7563310.9375, 7566651.5625, 7566848.4375, 7582910.9375, 7586281.25, 7587695.3125, 7588598.4375, 7589415.625, 7589984.375, 7591601.5625, 7605498.4375, 7618710.9375, 7620423.4375, 7621089.0625, 7621237.5, 7621331.25, 7621521.875, 7622914.0625, 7624485.9375, 7624565.625, 7624603.125, 7624728.125, 7625617.1875, 7627723.4375, 7630253.125, 7630934.375, 7631220.3125, 7632337.5, 7688767.1875, 7689118.75, 7691126.5625, 7693764.0625, 7694540.625, 7697353.125, 7697673.4375, 7698129.6875, 7698875.0, 7702300.0, 7702573.4375, ...], [10.713051006161061, 46.205314157275645, 8.201989650568361, 39.191171277572906, 73.81780751787149, 16.290522894242418, 11.639361768222534, 9.738681771399184, 5.9691967392463745, 6.732654842965873, 29.352996094753685, 40.12125904893643, 12.017598496462186, 68.34107171205457, 37.147371835394026, 94.99008762392955, 80.77086210229642, 49.65768933980893, 51.86516697791328, 8.451240349931041, 22.88988921125111, 25.353349388135815, 16.118568391287695, 189.04782623262463, 79.7423608681323, 14.734670586719915, 15.748960378227371, 69.25263500929465, 33.96008278542762, 18.718211988398764, 11.07439989541346, 59.49617827599503, 31.518008305291136, 21.048199636278362, 83.00129900642807, 28.759538032206706, 47.89007404220147, 25.234124368296882, 5.203005334323263, 69.95991192049884, 12.353724370063453, 52.0919660819248, 65.9823739755005, 17.073080500997563, 9.330219102883145, 162.02498344023223, 21.573660930648437, 9.578592551035683, 45.58965015493776, 5.919388677843027, 20.3561706082224, 57.65207190589465, 9.592643509130072, 19.757382372328532, 65.12601775090488, 17.362042864366654, 23.134872640463982, 13.25350392272108, 14.741929745440705, 34.22820653433416, 5.896698899301772, 107.76521471397601, 70.70603547451226, 59.029288992108995, 24.00099199192772, 132.08803994529163, 24.140164981141535, 18.29872189918449, 142.06813615802804, 12.073728390312809, 31.390276475147832, 11.261277487943248, 8.418631101847945, 56.18341850578175, 38.867383485099715, 32.15888864659335, 7.751113626922416, 9.663634067209365, 77.12737849232339, 28.988637712743635, 15.517399424267868, 40.397779223349886, 98.97913608797337, 16.700759144972917, 7.09302357784072, 81.95291949975585, 13.134130942625509, 25.140228873830974, 15.603974882966293, 136.36749136434074, 51.12327276518053, 31.85876848629799, 18.884702434330613, 40.79936324334394, 22.10295502501878, 15.157816701191905, 37.634204285703255, 17.00835975278213, 6.601492450544085, 21.783283143651737, 121.20824694432605, 40.91605215524861, 14.163149939540903, 59.87050462976426, 25.903360670210766, 135.3827979296443, 15.440225188024213, 55.26969908813271, 38.38513297567471, 30.966463387672768, 44.336218752106795, 7.400767357764462, 10.425683430908393, 120.02755428289045, 9.707186410030868, 55.64219926303227, 93.75080651201856, 26.256721822094605, 48.52290975704314, 264.47327012391105, 32.149619814385154, 16.595032312372005, 55.10612886869767, 30.482279644392683, 6.738721400627129, 17.14912541764162, 31.036769729117083, 57.67941455871323, 32.54465771788227, 18.73624352390624, 48.473051302894035, 5.124025258869424, 38.24449133353264, 99.12449408006285, 23.77179586760235, 17.196367781490746, 19.911650030215483, 10.064222237947954, 73.4298316205631, 11.705310106359512, 9.407745729499934, 32.96836339263979, 50.949771078810436, 12.942495050154493, 48.43369718198619, 15.323182683336192, 17.70932693795569, 13.08826764799719, 55.167752461682426, 54.798207918968025, 7.101952636019721, 14.430295221770141, 5.851312221704075, 18.53686804482017, 13.294598980874488, 54.064717740056295, 45.18571726350412, 68.88784114637386, 69.26914843107113, 33.093727749161204, 31.724061398414996, 118.1231016667246, 47.017646445387115, 8.837330458194552, 32.98303372226831, 7.641316819198258, 22.192540683773082, 56.05726272069935, 54.02239980189213, 82.56865532363777, 19.772048089181393, 55.792751228464994, 83.05263997823434, 11.280276607121538, 67.12335933863798, 40.31528544501474, 82.24858064103832, 50.43358457325405, 22.59420590158868, 131.01569162943122, 20.047218641106586, 71.7253569837988, 15.739605471539642, 5.134116227523508, 17.837474786412624, 47.78823245688457, 7.525876678995922, 15.96981254426212, 77.54406542175451, 33.132638452788974, 13.121468741569801, 32.901354497422886, 11.39277828244078, 12.546123099899562, 9.023400581507552, 30.31271131870807, 150.78537421664808, 25.534395563008463, 13.674889303847308, 20.91758858731326, 5.052926636921982, 29.507661278011287, 11.394459606525716, 155.34331456051294, 42.64896115966768, 15.229340705011019, 9.203719082034143, 10.539028624710923, 17.103138739597494, 22.135482322808517, 19.800250168423556, 51.4368507962477, 19.868799525911314, 10.954782941387526, 22.475431895625636, 28.548947615854463, 12.145965228984089, 74.74050633683366, 29.566695151328346, 6.946035610338883, 142.1554719517788, 24.78160114181562, 45.85085389639405, 34.358122359079275, 79.92529124611775, 20.46638653961609, 9.906649331132638, 19.81841686857041, 5.726247656544598, 24.54430034455945, 6.379397206152645, 40.510494070799304, 51.281405098744315, 20.680554179442353, 296.1861297905942, 19.05118378806145, 68.51594384504914, 28.355964016934664, 5.350764833392888, 92.17253780908504, 45.49380559753511, 18.73108758539995, 205.4873664774685, 120.89439055706512, 34.69866804561807, 44.12746711471958, 12.917345050915397, 24.67971460869807, 8.626418306704283, 191.42206980347584, 14.039874578054588, 52.07264398030967, 7.6894169588307815, 37.85302745629882, 24.838355098441056, 22.75420992262668, 13.437049805430775, 9.268988559588795, 14.124137656345068, 18.109515646692206, 16.175644694466502, 12.894940220278274, 66.68462418310563, 5.669314021639514, 11.447520965683479, 13.213117412481651, 51.085618237203676, 20.18919466302434, 99.28238726810977, 14.586143524237917, 283.1449188438106, 18.001978262989546, 77.64648024003074, 26.226342544448393, 24.99042037967624, 313.19531193083225, 141.0131251480738, 27.857801846489703, 18.234903529750422, 16.316697374267136, 231.4565092097347, 179.81987208164455, 28.28995158924046, 54.514717849607706, 40.52120878621977, 28.797343190549597, 27.48995766890571, 74.13655450852944, 45.25076554609947, 55.99258027750915, 47.77164415650097, 36.63632405448321, 32.242423537409266, 23.522335622840085, 52.83308057403912, 328.32384395425595, 11.939908875386346, 14.136301356810709, 58.37861561479734, 79.26308643772694, 70.3351997959965, 49.626451791977395, 46.958960186247666, 110.96451829340843, 40.381943307575874, 22.40470757081033, 10.47793037456153, 29.682704915094867, 93.68212454764652, 50.85295863018345, 56.456409672097834, 95.84940216200175, 14.95015360259729, 379.70048444526566, 6.529676474082922, 32.27906559007781, 82.92487483240097, 87.96660119122541, 12.762061189935608, 59.41596727412542, 36.15318989393144, 30.41785250202821, 68.63049143630359, 33.10722025645674, 22.57892658336801, 17.70560157250341, 7.625062255641228, 8.26166392855743, 25.807839237669526, 86.50703294937978, 27.60413260739759, 111.51652688823657, 8.061709977128482, 14.899546882613032, 67.63536240221649, 6.8062329426472346, 42.982243879826186, 20.070293216136466, 6.527426905694235, 37.95197918099651, 77.30050769500922, 22.375874523374698, 58.24987553488152, 35.84912053376168, 21.16291869976248, 5.03502414482645, 26.155160213158435, 109.96334972003814, 7.62634753375584, 116.07139562465346, 38.828670279264045, 6.029506614058144, 42.55839983747294, 56.565471128398514, 20.106555526303847, 33.77668339732171, 61.79541534894802, 8.585750934981327, 15.0622448201121, 6.000270182131832, 8.420863653511608, 7.155596986445784, 25.645215785354807, 26.121538511683315, 6.793379708982833, 50.1641123465557, 5.127070269467935, 12.330319894325836, 31.687069188073227, 122.11994585244659, 60.545315244154764, 76.76749326871932, 50.685468111187674, 143.35267532111848, 39.38121490414949, 37.02799323642322, 135.979358651345, 33.97263260345547, 32.794849967143264, 27.154328943240433, 20.925360369176627, 34.50314110537914, 37.475828074527946, 28.056586067061883, 28.066003356594525, 24.017830000825622, 27.32049126712908, 21.71515180283756, 41.97998668558744, 34.79487773341015, 35.420491214506185, 6.305705988921438, 78.397551373021, 70.50643435520176, 11.034800962326024, 10.668849590083159, 79.272262314767, 16.24944480851766, 196.76980416042454, 21.46195123958132, 111.24738326437202, 46.64736898553465, 183.57719852993728, 113.66004491997514, 97.97878328659418, 31.68111211200847, 12.044697999352595, 11.106111628398567, 66.32695527335767, 25.59211031956925, 70.09018695575686, 28.710383871100476, 14.868243461216393, 6.000420592400463, 24.198040164545528, 57.1730844450471, 184.90810273916534, 20.318946394143953, 36.68504803572285, 44.294204931928164, 53.1834723229263, 104.49827218279741, 115.73453240930353, 12.612518523384017, 64.51194759560335, 81.96107970734764, 255.96977045809686, 49.38730799671227, 54.19476850669163, 231.82897518128703, 36.54336095702415, 13.026670412491947, 10.417742924261002, 181.37464157412927, 58.259611082480276, 126.26774134640492, 84.37277195539801, 24.69831180324078, 37.008529070798524, 65.66973085176188, 36.98518035476906, 52.51536119677562, 26.93605394127892, 16.202097655007115, 6.300877373928769, 389.7069815251737, 31.765554356555302, 7.140733162062209, 47.67390041423606, 47.811796044031475, 23.630395525902887, 40.20396416511805, 12.833807220158466, 22.4965966451724, 6.439573416346811, 185.7347891447746, 72.60552141514613, 19.53454465026214, 15.978238808239569, 13.063469704766954, 6.265258188180271, 34.1107527538008, 63.19854403276283, 6.616648342812059, 113.60911820641822, 6.566208184461597, 75.11241150499096, 60.77075102224196, 42.061599035983775, 46.103264391408565, 87.4182232995216, 141.07774795394573, 151.1320118820747, 481.2416278638292, 42.42607186288484, 34.345017876046555, 54.75459215206463, 72.45773191237308, 100.84666572906005, 23.549525770417777, 8.489577528931774, 24.551929978171806, 12.888766658763322, 49.68381203692404, 15.422285564937944, 125.65384138465679, 30.397966951469712, 27.270698420745653, 27.099048518054474, 156.8463684831965, 7.915983211087909, 187.60012400997834, 96.46589167967049, 71.89903790577651, 9.139245069329956, 127.91713457575737, 52.49748303051388, 27.90779716937122, 8.576253750437305, 65.16999107734897, 69.96772077381678, 43.934527506734725, 50.06803511995554, 42.41229435511866, 11.337792535260055, 92.1153742574229, 17.179003256449953, 23.41436647679125, 13.112990431266166, 64.63599102137371, 20.255723945897284, 75.35798583234478, 102.06125894702285, 18.580745060513834, 26.5161773152523, 23.313901383660426, 192.1034565216484, 46.511801533349434, 72.1324346430014, 75.45665074422266, 34.93598624567049, 79.79150228418723, 10.345334337714803, 38.01825854766095, 124.34573559411038, 43.4923312979913, 5.054635033502185, 49.767462485557246, 8.229730595017402, 68.04143004524683, 28.810070982190513, 21.83451257561652, 6.903511458030656, 16.736564749167666, 49.70731867312349, 15.308672873281234, 21.189865449223387, 70.34183429676627, 43.39422817619453, 70.40887313326093, 37.3121000955872, 47.596328662981335, 52.855897425360155, 44.192992087668685, 6.377305485498518, 79.68320926206587, 32.99096599656137, 209.18475897163356, 32.55436768167911, 41.30594383100017, 34.61974854568331, 18.951698090338887, 41.55385547276799, 57.223889772184755, 47.51120842333313, 123.41441812335937, 6.251562586175794, 39.33257676771784, 130.35141222254927, 135.69912137809462, 131.99242029176696, 70.12218105111806, 45.49992224806314, 11.77690941529406, 5.83595202666188, 9.237351289526366, 37.542151875532205, 11.765170166005976, 28.461890034519204, 5.828350662928304, 23.065277487636983, 8.198526763505894, 18.860360721560703, 105.36819942755818, 39.181707898487595, 65.26223606272829, 29.398127665671414, 16.499785571056965, 110.5558867030779, 19.64329093700004, 44.21781087804614, 71.3331350597563, 38.66510254643331, 5.839161097082865, 6.669106137391587, 326.37849793012384, 26.606931290969335, 19.86139772610958, 22.537071039979068, 88.37353029003185, 11.476423063622484, 22.261595573223605, 11.580915956600872, 9.85048496065002, 15.639935324056893, 12.280741534153266, 23.521632785283376, 20.465935167624316, 13.632020121429225, 26.39937274650152, 31.255116337012236, 145.08090353834717, 140.2831203992128, 66.71878119903745, 32.15599793364798, 25.804201896748015, 111.47960093246029, 40.642368267183244, 145.9095509207118, 28.88035101898782, 49.722527268790884, 6.872663031826666, 21.17337126008074, 15.90655224561645, 42.806297764758504, 14.328425763040805, 14.496560346618661, 17.673396803032606, 19.23121050859685, 24.897494331191595, 32.52686848933328, 80.14849275562675, 63.374577430291154, 76.91140856192244, 27.893827261715444, 63.465544434511216, 69.95905019568605, 26.2949500155077, 9.119469239515121, 5.557785809486354, 16.44533598579733, 54.88177842253924, 14.341521226574372, 251.0304916112949, 17.63651405487081, 139.08284683084594, 13.838233465237819, 17.46843995361195, 202.38405201465744, 53.53642578573328, 68.98051038271804, 15.747750093205322, 27.764991297744718, 30.29204731258981, 53.5370907479409, 25.85808266064117, 57.78752997500231, 173.52764229644686, 35.334995101282075, 27.37846126567299, 10.502481413938463, 62.186478613708296, 86.59540085177463, 328.95088187580143, 104.7686347599838, 42.075123306795035, 34.89136589201932, 62.56467527308959, 21.992481583954355, 52.384306875634444, 6.898234352915194, 53.054096533400184, 37.004748286637835, 57.48910312922262, 34.49272234952973, 40.655072862253874, 153.23832597706482, 168.47460004661588, 153.27092334282594, 14.374096624582021, 17.47522131062833, 143.58984890130242, 56.75249100435342, 111.69987917328646, 16.956956576584297, 26.53271236146316, 33.09078923924028, 71.77839587525273, 5.55127303476317, 5.445040312390084, 7.458309920179268, 60.53099909465615, 14.606821962128965, 289.4249005870328, 121.55731411637484, 25.88391820732727, 8.125076584329253, 18.85500784395236, 6.235379320506899, 38.67121915893703, 41.83502730627808, 5.173751914269115, 6.871221080242857, 86.80841677891489, 8.836705979776102, 23.809406110862543, 5.164961056748544, 150.78142407126924, 103.29076974180371, 38.90818641242245, 55.487481329300515, 33.982875039335845, 30.944194456726795, 13.296989899127201, 11.22075936182022, 5.745755269611203, 8.433031025626377, 225.20154284735193, 110.78413961479244, 186.7946777581286, 20.214631961667536, 21.639784958095436, 8.053923690355866, 31.57656468171556, 14.931975096864418, 16.4062582202943, 16.3396321839484, 31.467426023937065, 89.13562022617785, 139.5754481616569, 108.80033218047721, 11.475977284580075, 37.02669229079955, 25.01959486495189, 61.346937726608175, 31.112180757678907, 29.26337756784772, 12.15341106885423, 58.788523868141525, 41.82524109223579, 46.73677208802073, 29.49497932166495, 13.230050691870435, 70.49728900846466, 122.6009624982872, 10.175442916362238, 7.541080455753765, 125.19608917609159, 7.096837181987375, 34.57975638827786, 48.99719827741131, 20.136884995886334, 28.108610525583657, 47.54586000802074, 5.739948487806487, 288.8359434216485, 20.80451714897945, 22.580838688193786, 96.91243391844982, 26.530617682188282, 46.690585571946265, 45.27015687932401, 211.274087895192, 58.07436142327462, 25.48902365489932, 38.07328197113766, 28.779674924885253, 80.58950622076149, 6.83147885470799, 21.278179439616775, 78.39583959034182, 70.28697519892509, 218.7195807252162, 51.07795831732477, 23.46650137533203, 33.12777459382858, 8.804583687427268, 80.24445735634535, 13.644637212905879, 38.909463491319926, 61.848413750218775, 169.96402675900413, 56.42003923261203, 66.2534085845368, 46.5832559148559, 12.011346834521214, 43.61758057242759, 199.22996960407045, 6.580800752144334, 5.030624536721584, 122.46700338561193, 36.57596409469311, 8.741658173310547, 18.15357477040057, 96.05667750502765, 61.48159568342763, 24.39197526465823, 44.33280770260135, 84.72307413420737, 87.84117481893159, 42.96536766418469, 36.654240988005704, 5.327699324012194, 57.62838360659262, 111.35263958569983, 26.299305105467557, 9.936778144116541, 61.069909701123386, 50.48315224668051, 26.432696377515615, 21.956056230169764, 46.12688829050798, 192.232185543653, 53.98895843549289, 117.58406511752338, 47.157211422830784, 26.835375552712343, 34.05657477980806, 25.52157426455279, 19.858693098137227, 60.295374911785345, 17.6809760633518, 11.044201874916757, 62.01950067362743, 85.96438612221726, 19.940851847160328, 46.1364693115263, 126.10836485018214, 80.99269685157304, 190.76372824917453, 6.23722355704399, 76.77517684781239, 406.3867758872526, 101.45198344937525, 29.274656666931108, 26.768638631398904, 68.09728441020489, 115.18131491261298, 34.9214361382527, 19.373091706535668, 14.796044290538244, 8.712554247472482, 31.44831113221531, 51.967921962611, 88.87906553398393, 85.03526935983585, 51.51420000559395, 51.05064703489565, 63.55823802365108, 13.754074412960426, 22.325844656050656, 61.24893113655137, 35.268784703835294, 62.90902100836937, 39.629975388239416, 25.502124489530342, 59.02010039311995, 94.00603166354365, 84.41621604192699, 6.13023169476253, 40.27523507904435, 25.129034785069255, 9.128792590726567, 56.35050920148815, 6.575522898248985, 15.282615937155223, 20.82855333306505, 70.94604373267254, 127.19018506462595, 65.52241710444846, 25.550397126110965, 68.16639292822502, 56.603025582534755, 36.105150965970275, 97.26428449003525, 10.603097642580483, 78.9361266641682, 45.86227197308742, 5.324028114150204, 84.57367181752093, 64.04487260358093, 22.902032410225743, 16.69681658522037, 17.59842258252129, 75.52904057316563, 6.708320852032042, 64.05437088213958, 8.42227866180616, 24.755531231071302, 112.26872133137002, 296.1840148640123, 5.79017954419084, 18.52637662503884, 56.711067341812424, 5.085310704130503, 74.49077463939767, 25.697846124798815, 60.89298273185257, 64.60763011398305, 49.448039589481034, 27.342507289003276, 81.20909486182096, 60.13197453543258, 23.876764430505887, 18.844443423359284, 94.9558563009197, 56.74727530891187, 54.38466495621606, 156.03950602306492, 18.08013873755965, 9.403597119548579, 20.86511772100988, 48.612463731101826, 40.68873780723625, 13.177970389348966, 19.710136076767288, 33.02868895676784, 36.37187648588037, 5.197196798069252, 134.23758526834223, 5.194657368466044, 50.359583428338496, 6.130650320560632, 20.621866376959694, 7.827290815186451, 56.58968508150255, 27.838495016141508, 32.64924536459766, 57.10647258336654, 95.47082935762646, 65.0469683102413, 87.17718458077739, 143.45069910904334, 55.25145213608673, 8.075166702332739, 13.984158600434942, 28.163170454450643, 39.742279846785465, 16.447943599481203, 26.344085638217773, 7.9014542191118045, 33.6196909709104, 8.378331185072366, 12.025680777069136, 72.14012793695497, 26.296930324574895, 27.435283377491807, 24.655539542425203, 38.21687445841762, 49.04456986496268, 12.107499539665593, 10.415604643372966, 83.50004305336114, 22.172590894844113, 88.5271409191658, 9.259784857972262, 183.70537593482115, 104.36363148185106, 7.696088935860922, 39.52729238475217, 12.140113259435068, 16.91736667687468, 97.6734639670456, 9.875356756250646, 17.175779146125794, 6.035715432654942, 131.42689031845754, 13.29891426661904, 11.490635711230624, 26.231679848766746, 13.480727851579072, 72.92436678297061, 29.83397458290599, 93.3119589020483, 214.57346518298095, 38.792456305738035, 64.08595584596128, 6.231261770904037, 26.81023192452331, 52.37849212813234, 338.7829110772877, 80.64909782121588, 13.691368987871545, 58.784159352082575, 41.33406431604816, 60.22604450369563, 123.10683518221651, 68.60351925005239, 104.52289419445616, 20.204207069586374, 8.31501541416413, 37.85999763025312, 41.22648052905725, 38.77118722406425, 7.010181070636918, 43.968786026681556, 24.041916461648107, 30.730799642842165, 18.54029988339476, 12.008210614430977, 44.94171361638127, 30.52742582369338, 17.616801392060374, 59.86572233242134, 34.440894620651065, 22.871166407260183, 15.655471628935373, 17.372029099047424, 46.52590805969684, 73.01809164742689, ...])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([2924346.875, 2944565.625, 2986623.4375, 3002956.25, 3147323.4375, 3397809.375, 3436142.1875, 3999750.0, 4792262.5, 5266640.625, 5416700.0, 5455485.9375, 5455493.75, 5456928.125, 5734228.125, 5736831.25, 5759906.25, 5788390.625, 5790696.875, 5795334.375, 5797781.25, 5839612.5, 5879600.0, 5883587.5, 5910326.5625, 5915831.25, 5917845.3125, 5924237.5, 5925865.625, 5926953.125, 5927317.1875, 5927843.75, 5929675.0, 5931595.3125, 5934132.8125, 5934760.9375, 5936929.6875, 5937835.9375, 5939276.5625, 5943960.9375, 5944914.0625, 5945410.9375, 5946307.8125, 5956275.0, 5956276.5625, 5958031.25, 5959417.1875, 5959895.3125, 5960125.0, 5960753.125, 5966514.0625, 5970970.3125, 5972035.9375, 5977026.5625, 5977026.5625, 5978792.1875, 5978793.75, 5979534.375, 5980356.25, 5980540.625, 5980615.625, 5980959.375, 5981001.5625, 5981025.0, 5981032.8125, 5981273.4375, 5981423.4375, 5981506.25, 5981570.3125, 5981609.375, 5981657.8125, 5981925.0, 5982148.4375, 5982321.875, 5982434.375, 5982484.375, 5982528.125, 5982728.125, 5983053.125, 5983329.6875, 5983426.5625, 5983896.875, 5983953.125, 5984093.75, 5984282.8125, 5984553.125, 5985385.9375, 5986021.875, 5986590.625, 5990492.1875, 5990950.0, 5990971.875, 6009821.875, 6015390.625, 6023539.0625, 6024298.4375, 6054309.375, 6055198.4375, 6058721.875, 6059100.0, 6063656.25, 6066542.1875, 6078581.25, 6084646.875, 6084657.8125, 6088356.25, 6091434.375, 6093745.3125, 6093850.0, 6096214.0625, 6098381.25, 6099982.8125, 6100262.5, 6124200.0, 6126929.6875, 6128434.375, 6128501.5625, 6128693.75, 6132237.5, 6132362.5, 6132404.6875, 6132954.6875, 6134525.0, 6138635.9375, 6142756.25, 6170331.25, 6175581.25, 6175920.3125, 6175925.0, 6176109.375, 6176148.4375, 6176251.5625, 6176275.0, 6177120.3125, 6177282.8125, 6177342.1875, 6177582.8125, 6177603.125, 6178493.75, 6179835.9375, 6180748.4375, 6180907.8125, 6181957.8125, 6182935.9375, 6183243.75, 6183395.3125, 6184270.3125, 6185128.125, 6186509.375, 6194025.0, 6209742.1875, 6210400.0, 6217031.25, 6218221.875, 6262400.0, 6267309.375, 6271317.1875, 6285198.4375, 6286095.3125, 6287734.375, 6291410.9375, 6292990.625, 6298095.3125, 6301600.0, 6303882.8125, 6306959.375, 6336956.25, 6338228.125, 6341228.125, 6342273.4375, 6343179.6875, 6343450.0, 6343943.75, 6344215.625, 6344232.8125, 6344571.875, 6344598.4375, 6344617.1875, 6344620.3125, 6344890.625, 6345023.4375, 6345162.5, 6345195.3125, 6345231.25, 6345573.4375, 6345728.125, 6345757.8125, 6345910.9375, 6346014.0625, 6346209.375, 6346900.0, 6347385.9375, 6347401.5625, 6347420.3125, 6347785.9375, 6348076.5625, 6348728.125, 6348964.0625, 6349156.25, 6349192.1875, 6349276.5625, 6349921.875, 6350020.3125, 6351696.875, 6352292.1875, 6352331.25, 6352590.625, 6352612.5, 6353323.4375, 6355182.8125, 6355465.625, 6355529.6875, 6356089.0625, 6356146.875, 6356576.5625, 6356970.3125, 6358621.875, 6359306.25, 6359323.4375, 6360632.8125, 6362714.0625, 6372925.0, 6373746.875, 6377125.0, 6387075.0, 6387101.5625, 6395089.0625, 6447707.8125, 6481442.1875, 6483025.0, 6484126.5625, 6485400.0, 6486600.0, 6487826.5625, 6495279.6875, 6496150.0, 6496159.375, 6496362.5, 6496417.1875, 6496687.5, 6497096.875, 6497334.375, 6497967.1875, 6498026.5625, 6498237.5, 6498279.6875, 6498581.25, 6498909.375, 6499729.6875, 6499948.4375, 6500839.0625, 6508262.5, 6512271.875, 6514573.4375, 6516168.75, 6516629.6875, 6518032.8125, 6523879.6875, 6527639.0625, 6529360.9375, 6541267.1875, 6542935.9375, 6544412.5, 6545184.375, 6547835.9375, 6548075.0, 6550118.75, 6550818.75, 6551314.0625, 6551706.25, 6551853.125, 6552904.6875, 6553110.9375, 6553231.25, 6553465.625, 6553667.1875, 6553776.5625, 6554157.8125, 6554596.875, 6554623.4375, 6554843.75, 6555184.375, 6555467.1875, 6556028.125, 6556085.9375, 6556203.125, 6556267.1875, 6556943.75, 6556965.625, 6558034.375, 6562751.5625, 6566706.25, 6567043.75, 6581642.1875, 6584039.0625, 6587303.125, 6588934.375, 6590175.0, 6590446.875, 6593926.5625, 6596121.875, 6598346.875, 6598392.1875, 6599023.4375, 6599153.125, 6599817.1875, 6600389.0625, 6600815.625, 6601582.8125, 6607121.875, 6616378.125, 6616556.25, 6616979.6875, 6617123.4375, 6618862.5, 6623275.0, 6642670.3125, 6647689.0625, 6652178.125, 6674725.0, 6683126.5625, 6695668.75, 6703696.875, 6704498.4375, 6719276.5625, 6719456.25, 6719482.8125, 6719523.4375, 6719896.875, 6720107.8125, 6722787.5, 6725245.3125, 6725745.3125, 6726210.9375, 6729340.625, 6730289.0625, 6733443.75, 6733454.6875, 6733928.125, 6734906.25, 6734996.875, 6735070.3125, 6736481.25, 6737221.875, 6737251.5625, 6737345.3125, 6737987.5, 6742279.6875, 6745264.0625, 6745907.8125, 6746176.5625, 6747793.75, 6747860.9375, 6749214.0625, 6750520.3125, 6751846.875, 6751926.5625, 6752326.5625, 6753300.0, 6753306.25, 6759739.0625, 6763539.0625, 6765523.4375, 6767343.75, 6768370.3125, 6769764.0625, 6769864.0625, 6770679.6875, 6776676.5625, 6778278.125, 6779393.75, 6779853.125, 6781718.75, 6781939.0625, 6782900.0, 6783450.0, 6784067.1875, 6785901.5625, 6785929.6875, 6787056.25, 6787090.625, 6787775.0, 6788421.875, 6788440.625, 6788642.1875, 6789384.375, 6789495.3125, 6789673.4375, 6790017.1875, 6791360.9375, 6791470.3125, 6791640.625, 6791792.1875, 6791828.125, 6791956.25, 6792179.6875, 6792260.9375, 6792781.25, 6793235.9375, 6793946.875, 6794045.3125, 6795615.625, 6809879.6875, 6811578.125, 6812414.0625, 6813021.875, 6815412.5, 6819867.1875, 6820987.5, 6821856.25, 6821895.3125, 6822035.9375, 6823390.625, 6824665.625, 6824876.5625, 6825195.3125, 6826151.5625, 6826357.8125, 6826454.6875, 6826839.0625, 6826850.0, 6827001.5625, 6827818.75, 6827881.25, 6828200.0, 6828453.125, 6828470.3125, 6829500.0, 6831670.3125, 6832375.0, 6833915.625, 6834368.75, 6836112.5, 6840454.6875, 6840464.0625, 6840509.375, 6840510.9375, 6840934.375, 6840937.5, 6841137.5, 6841140.625, 6841253.125, 6842260.9375, 6844643.75, 6847510.9375, 6848787.5, 6851209.375, 6852473.4375, 6853573.4375, 6853928.125, 6858198.4375, 6858671.875, 6859639.0625, 6859907.8125, 6862125.0, 6863965.625, 6864787.5, 6864960.9375, 6865084.375, 6865190.625, 6865206.25, 6865842.1875, 6868729.6875, 6869650.0, 6873051.5625, 6873943.75, 6874078.125, 6876496.875, 6876517.1875, 6878250.0, 6879478.125, 6879828.125, 6882109.375, 6883210.9375, 6893704.6875, 6903237.5, 6907682.8125, 6909642.1875, 6910242.1875, 6910378.125, 6910415.625, 6910529.6875, 6914932.8125, 6915995.3125, 6916112.5, 6916350.0, 6916507.8125, 6917054.6875, 6917259.375, 6917795.3125, 6918117.1875, 6918154.6875, 6918909.375, 6918934.375, 6918957.8125, 6919028.125, 6919071.875, 6919606.25, 6919785.9375, 6919873.4375, 6919948.4375, 6919978.125, 6920162.5, 6920246.875, 6920254.6875, 6920275.0, 6920364.0625, 6920423.4375, 6920693.75, 6920917.1875, 6920948.4375, 6921317.1875, 6921323.4375, 6921364.0625, 6921839.0625, 6921879.6875, 6921887.5, 6921920.3125, 6921942.1875, 6921964.0625, 6921968.75, 6922018.75, 6922120.3125, 6922253.125, 6922284.375, 6922292.1875, 6922354.6875, 6922448.4375, 6922482.8125, 6922507.8125, 6922600.0, 6922673.4375, 6922675.0, 6922767.1875, 6922793.75, 6922796.875, 6922907.8125, 6922909.375, 6922978.125, 6923031.25, 6923357.8125, 6923382.8125, 6923385.9375, 6923448.4375, 6923642.1875, 6923706.25, 6923723.4375, 6923778.125, 6923943.75, 6924182.8125, 6924273.4375, 6924296.875, 6924346.875, 6924639.0625, 6924685.9375, 6924710.9375, 6924765.625, 6925482.8125, 6925582.8125, 6925870.3125, 6926137.5, 6926567.1875, 6926606.25, 6926610.9375, 6926665.625, 6926670.3125, 6926801.5625, 6926859.375, 6926965.625, 6927110.9375, 6927850.0, 6928290.625, 6928543.75, 6928693.75, 6929810.9375, 6929990.625, 6931276.5625, 6931575.0, 6931692.1875, 6931715.625, 6931985.9375, 6932751.5625, 6933606.25, 6934704.6875, 6934714.0625, 6936109.375, 6936382.8125, 6937206.25, 6937353.125, 6937642.1875, 6937823.4375, 6938221.875, 6938401.5625, 6938945.3125, 6940976.5625, 6941468.75, 6942489.0625, 6943125.0, 6943398.4375, 6943618.75, 6943878.125, 6944743.75, 6945048.4375, 6945275.0, 6945312.5, 6945375.0, 6945657.8125, 6946314.0625, 6946690.625, 6946729.6875, 6948473.4375, 6950375.0, 6950990.625, 6951054.6875, 6953428.125, 6953579.6875, 6954956.25, 6955893.75, 6956270.3125, 6957575.0, 6958329.6875, 6958557.8125, 6958882.8125, 6961229.6875, 6961531.25, 6962248.4375, 6963715.625, 6966657.8125, 6972007.8125, 6972071.875, 6972503.125, 6972506.25, 6973884.375, 6976331.25, 6977882.8125, 6980239.0625, 6981967.1875, 6982017.1875, 6985701.5625, 6986907.8125, 6996946.875, 7003584.375, 7021931.25, 7022676.5625, 7024148.4375, 7024573.4375, 7024801.5625, 7026842.1875, 7027198.4375, 7027773.4375, 7028279.6875, 7029885.9375, 7030929.6875, 7031110.9375, 7040485.9375, 7041504.6875, 7043698.4375, 7043746.875, 7044601.5625, 7050225.0, 7050925.0, 7053329.6875, 7054293.75, 7054600.0, 7056348.4375, 7058498.4375, 7058601.5625, 7060473.4375, 7062290.625, 7063640.625, 7065537.5, 7066434.375, 7066515.625, 7067143.75, 7068031.25, 7068701.5625, 7068787.5, 7069237.5, 7069501.5625, 7069520.3125, 7069556.25, 7069771.875, 7072853.125, 7074782.8125, 7075415.625, 7076504.6875, 7077718.75, 7079539.0625, 7081946.875, 7082650.0, 7083671.875, 7083726.5625, 7084162.5, 7084609.375, 7087012.5, 7087178.125, 7087289.0625, 7108465.625, 7109739.0625, 7110051.5625, 7112256.25, 7118868.75, 7120237.5, 7120610.9375, 7121957.8125, 7122028.125, 7122396.875, 7126531.25, 7126537.5, 7127698.4375, 7127884.375, 7127981.25, 7128231.25, 7128260.9375, 7128593.75, 7131671.875, 7132656.25, 7134743.75, 7134743.75, 7138676.5625, 7142254.6875, 7144043.75, 7144675.0, 7144685.9375, 7145096.875, 7145753.125, 7146346.875, 7146564.0625, 7146718.75, 7146743.75, 7147043.75, 7147646.875, 7148246.875, 7148701.5625, 7149018.75, 7149393.75, 7150721.875, 7150729.6875, 7151068.75, 7151268.75, 7152187.5, 7152265.625, 7152637.5, 7152928.125, 7153235.9375, 7154014.0625, 7154770.3125, 7156054.6875, 7156275.0, 7156303.125, 7157048.4375, 7157450.0, 7157885.9375, 7158431.25, 7158657.8125, 7159223.4375, 7159434.375, 7159506.25, 7159804.6875, 7160067.1875, 7160132.8125, 7160151.5625, 7160165.625, 7160301.5625, 7160365.625, 7160467.1875, 7160656.25, 7160710.9375, 7160714.0625, 7160921.875, 7161073.4375, 7161112.5, 7161284.375, 7161420.3125, 7161517.1875, 7161712.5, 7162232.8125, 7162370.3125, 7162492.1875, 7162581.25, 7162696.875, 7162795.3125, 7162846.875, 7163475.0, 7165593.75, 7167217.1875, 7167607.8125, 7167757.8125, 7169803.125, 7170253.125, 7170739.0625, 7176346.875, 7176409.375, 7176529.6875, 7176837.5, 7178376.5625, 7179979.6875, 7188281.25, 7192954.6875, 7196901.5625, 7205226.5625, 7206895.3125, 7208975.0, 7209975.0, 7210984.375, 7211687.5, 7212275.0, 7213364.0625, 7213571.875, 7214796.875, 7215300.0, 7216334.375, 7216740.625, 7216760.9375, 7216767.1875, 7217109.375, 7218109.375, 7218309.375, 7218712.5, 7219104.6875, 7219184.375, 7221693.75, 7221932.8125, 7222025.0, 7222228.125, 7222354.6875, 7222914.0625, 7223493.75, 7223623.4375, 7223759.375, 7224009.375, 7224017.1875, 7224107.8125, 7224135.9375, 7224290.625, 7225051.5625, 7225400.0, 7226271.875, 7226517.1875, 7226937.5, 7227267.1875, 7227476.5625, 7228623.4375, 7228885.9375, 7229714.0625, 7230523.4375, 7231045.3125, 7231375.0, 7231400.0, 7231829.6875, 7231848.4375, 7232993.75, 7233007.8125, 7233868.75, 7233898.4375, 7234257.8125, 7234362.5, 7235712.5, 7237325.0, 7238143.75, 7238920.3125, 7245321.875, 7247782.8125, 7249503.125, 7249884.375, 7250350.0, 7250675.0, 7250742.1875, 7250753.125, 7253918.75, 7254842.1875, 7258203.125, 7259095.3125, 7260581.25, 7260595.3125, 7263329.6875, 7264870.3125, 7267870.3125, 7270634.375, 7270667.1875, 7271054.6875, 7275165.625, 7287562.5, 7295409.375, 7295478.125, 7307853.125, 7312629.6875, 7314253.125, 7315639.0625, 7316850.0, 7316942.1875, 7317525.0, 7319484.375, 7320073.4375, 7322621.875, 7322870.3125, 7322926.5625, 7324837.5, 7325293.75, 7326329.6875, 7326343.75, 7326862.5, 7327868.75, 7329995.3125, 7330507.8125, 7330781.25, 7330864.0625, 7330965.625, 7330970.3125, 7331039.0625, 7331103.125, 7331132.8125, 7331268.75, 7331450.0, 7331490.625, 7331496.875, 7331556.25, 7331595.3125, 7331656.25, 7331831.25, 7331931.25, 7332026.5625, 7332139.0625, 7332150.0, 7332170.3125, 7332189.0625, 7332245.3125, 7332296.875, 7332376.5625, 7332570.3125, 7332795.3125, 7333228.125, 7333245.3125, 7333573.4375, 7349715.625, 7404610.9375, 7422159.375, 7423446.875, 7424810.9375, 7432853.125, 7433401.5625, 7444418.75, 7458937.5, 7459159.375, 7499614.0625, 7507142.1875, 7514735.9375, 7516075.0, 7526598.4375, 7530245.3125, 7533860.9375, 7537785.9375, 7539229.6875, 7546154.6875, 7547800.0, 7547820.3125, 7548732.8125, 7552317.1875, 7558478.125, 7558610.9375, 7563135.9375, 7563310.9375, 7566651.5625, 7566848.4375, 7582910.9375, 7586281.25, 7587695.3125, 7588598.4375, 7589415.625, 7589984.375, 7591601.5625, 7605498.4375, 7618710.9375, 7620423.4375, 7621089.0625, 7621237.5, 7621331.25, 7621521.875, 7622914.0625, 7624485.9375, 7624565.625, 7624603.125, 7624728.125, 7625617.1875, 7627723.4375, 7630253.125, 7630934.375, 7631220.3125, 7632337.5, 7688767.1875, 7689118.75, 7691126.5625, 7693764.0625, 7694540.625, 7697353.125, 7697673.4375, 7698129.6875, 7698875.0, 7702300.0, 7702573.4375, ...], [10.713051006161061, 46.205314157275645, 8.201989650568361, 39.191171277572906, 73.81780751787149, 16.290522894242418, 11.639361768222534, 9.738681771399184, 5.9691967392463745, 6.732654842965873, 29.352996094753685, 40.12125904893643, 12.017598496462186, 68.34107171205457, 37.147371835394026, 94.99008762392955, 80.77086210229642, 49.65768933980893, 51.86516697791328, 8.451240349931041, 22.88988921125111, 25.353349388135815, 16.118568391287695, 189.04782623262463, 79.7423608681323, 14.734670586719915, 15.748960378227371, 69.25263500929465, 33.96008278542762, 18.718211988398764, 11.07439989541346, 59.49617827599503, 31.518008305291136, 21.048199636278362, 83.00129900642807, 28.759538032206706, 47.89007404220147, 25.234124368296882, 5.203005334323263, 69.95991192049884, 12.353724370063453, 52.0919660819248, 65.9823739755005, 17.073080500997563, 9.330219102883145, 162.02498344023223, 21.573660930648437, 9.578592551035683, 45.58965015493776, 5.919388677843027, 20.3561706082224, 57.65207190589465, 9.592643509130072, 19.757382372328532, 65.12601775090488, 17.362042864366654, 23.134872640463982, 13.25350392272108, 14.741929745440705, 34.22820653433416, 5.896698899301772, 107.76521471397601, 70.70603547451226, 59.029288992108995, 24.00099199192772, 132.08803994529163, 24.140164981141535, 18.29872189918449, 142.06813615802804, 12.073728390312809, 31.390276475147832, 11.261277487943248, 8.418631101847945, 56.18341850578175, 38.867383485099715, 32.15888864659335, 7.751113626922416, 9.663634067209365, 77.12737849232339, 28.988637712743635, 15.517399424267868, 40.397779223349886, 98.97913608797337, 16.700759144972917, 7.09302357784072, 81.95291949975585, 13.134130942625509, 25.140228873830974, 15.603974882966293, 136.36749136434074, 51.12327276518053, 31.85876848629799, 18.884702434330613, 40.79936324334394, 22.10295502501878, 15.157816701191905, 37.634204285703255, 17.00835975278213, 6.601492450544085, 21.783283143651737, 121.20824694432605, 40.91605215524861, 14.163149939540903, 59.87050462976426, 25.903360670210766, 135.3827979296443, 15.440225188024213, 55.26969908813271, 38.38513297567471, 30.966463387672768, 44.336218752106795, 7.400767357764462, 10.425683430908393, 120.02755428289045, 9.707186410030868, 55.64219926303227, 93.75080651201856, 26.256721822094605, 48.52290975704314, 264.47327012391105, 32.149619814385154, 16.595032312372005, 55.10612886869767, 30.482279644392683, 6.738721400627129, 17.14912541764162, 31.036769729117083, 57.67941455871323, 32.54465771788227, 18.73624352390624, 48.473051302894035, 5.124025258869424, 38.24449133353264, 99.12449408006285, 23.77179586760235, 17.196367781490746, 19.911650030215483, 10.064222237947954, 73.4298316205631, 11.705310106359512, 9.407745729499934, 32.96836339263979, 50.949771078810436, 12.942495050154493, 48.43369718198619, 15.323182683336192, 17.70932693795569, 13.08826764799719, 55.167752461682426, 54.798207918968025, 7.101952636019721, 14.430295221770141, 5.851312221704075, 18.53686804482017, 13.294598980874488, 54.064717740056295, 45.18571726350412, 68.88784114637386, 69.26914843107113, 33.093727749161204, 31.724061398414996, 118.1231016667246, 47.017646445387115, 8.837330458194552, 32.98303372226831, 7.641316819198258, 22.192540683773082, 56.05726272069935, 54.02239980189213, 82.56865532363777, 19.772048089181393, 55.792751228464994, 83.05263997823434, 11.280276607121538, 67.12335933863798, 40.31528544501474, 82.24858064103832, 50.43358457325405, 22.59420590158868, 131.01569162943122, 20.047218641106586, 71.7253569837988, 15.739605471539642, 5.134116227523508, 17.837474786412624, 47.78823245688457, 7.525876678995922, 15.96981254426212, 77.54406542175451, 33.132638452788974, 13.121468741569801, 32.901354497422886, 11.39277828244078, 12.546123099899562, 9.023400581507552, 30.31271131870807, 150.78537421664808, 25.534395563008463, 13.674889303847308, 20.91758858731326, 5.052926636921982, 29.507661278011287, 11.394459606525716, 155.34331456051294, 42.64896115966768, 15.229340705011019, 9.203719082034143, 10.539028624710923, 17.103138739597494, 22.135482322808517, 19.800250168423556, 51.4368507962477, 19.868799525911314, 10.954782941387526, 22.475431895625636, 28.548947615854463, 12.145965228984089, 74.74050633683366, 29.566695151328346, 6.946035610338883, 142.1554719517788, 24.78160114181562, 45.85085389639405, 34.358122359079275, 79.92529124611775, 20.46638653961609, 9.906649331132638, 19.81841686857041, 5.726247656544598, 24.54430034455945, 6.379397206152645, 40.510494070799304, 51.281405098744315, 20.680554179442353, 296.1861297905942, 19.05118378806145, 68.51594384504914, 28.355964016934664, 5.350764833392888, 92.17253780908504, 45.49380559753511, 18.73108758539995, 205.4873664774685, 120.89439055706512, 34.69866804561807, 44.12746711471958, 12.917345050915397, 24.67971460869807, 8.626418306704283, 191.42206980347584, 14.039874578054588, 52.07264398030967, 7.6894169588307815, 37.85302745629882, 24.838355098441056, 22.75420992262668, 13.437049805430775, 9.268988559588795, 14.124137656345068, 18.109515646692206, 16.175644694466502, 12.894940220278274, 66.68462418310563, 5.669314021639514, 11.447520965683479, 13.213117412481651, 51.085618237203676, 20.18919466302434, 99.28238726810977, 14.586143524237917, 283.1449188438106, 18.001978262989546, 77.64648024003074, 26.226342544448393, 24.99042037967624, 313.19531193083225, 141.0131251480738, 27.857801846489703, 18.234903529750422, 16.316697374267136, 231.4565092097347, 179.81987208164455, 28.28995158924046, 54.514717849607706, 40.52120878621977, 28.797343190549597, 27.48995766890571, 74.13655450852944, 45.25076554609947, 55.99258027750915, 47.77164415650097, 36.63632405448321, 32.242423537409266, 23.522335622840085, 52.83308057403912, 328.32384395425595, 11.939908875386346, 14.136301356810709, 58.37861561479734, 79.26308643772694, 70.3351997959965, 49.626451791977395, 46.958960186247666, 110.96451829340843, 40.381943307575874, 22.40470757081033, 10.47793037456153, 29.682704915094867, 93.68212454764652, 50.85295863018345, 56.456409672097834, 95.84940216200175, 14.95015360259729, 379.70048444526566, 6.529676474082922, 32.27906559007781, 82.92487483240097, 87.96660119122541, 12.762061189935608, 59.41596727412542, 36.15318989393144, 30.41785250202821, 68.63049143630359, 33.10722025645674, 22.57892658336801, 17.70560157250341, 7.625062255641228, 8.26166392855743, 25.807839237669526, 86.50703294937978, 27.60413260739759, 111.51652688823657, 8.061709977128482, 14.899546882613032, 67.63536240221649, 6.8062329426472346, 42.982243879826186, 20.070293216136466, 6.527426905694235, 37.95197918099651, 77.30050769500922, 22.375874523374698, 58.24987553488152, 35.84912053376168, 21.16291869976248, 5.03502414482645, 26.155160213158435, 109.96334972003814, 7.62634753375584, 116.07139562465346, 38.828670279264045, 6.029506614058144, 42.55839983747294, 56.565471128398514, 20.106555526303847, 33.77668339732171, 61.79541534894802, 8.585750934981327, 15.0622448201121, 6.000270182131832, 8.420863653511608, 7.155596986445784, 25.645215785354807, 26.121538511683315, 6.793379708982833, 50.1641123465557, 5.127070269467935, 12.330319894325836, 31.687069188073227, 122.11994585244659, 60.545315244154764, 76.76749326871932, 50.685468111187674, 143.35267532111848, 39.38121490414949, 37.02799323642322, 135.979358651345, 33.97263260345547, 32.794849967143264, 27.154328943240433, 20.925360369176627, 34.50314110537914, 37.475828074527946, 28.056586067061883, 28.066003356594525, 24.017830000825622, 27.32049126712908, 21.71515180283756, 41.97998668558744, 34.79487773341015, 35.420491214506185, 6.305705988921438, 78.397551373021, 70.50643435520176, 11.034800962326024, 10.668849590083159, 79.272262314767, 16.24944480851766, 196.76980416042454, 21.46195123958132, 111.24738326437202, 46.64736898553465, 183.57719852993728, 113.66004491997514, 97.97878328659418, 31.68111211200847, 12.044697999352595, 11.106111628398567, 66.32695527335767, 25.59211031956925, 70.09018695575686, 28.710383871100476, 14.868243461216393, 6.000420592400463, 24.198040164545528, 57.1730844450471, 184.90810273916534, 20.318946394143953, 36.68504803572285, 44.294204931928164, 53.1834723229263, 104.49827218279741, 115.73453240930353, 12.612518523384017, 64.51194759560335, 81.96107970734764, 255.96977045809686, 49.38730799671227, 54.19476850669163, 231.82897518128703, 36.54336095702415, 13.026670412491947, 10.417742924261002, 181.37464157412927, 58.259611082480276, 126.26774134640492, 84.37277195539801, 24.69831180324078, 37.008529070798524, 65.66973085176188, 36.98518035476906, 52.51536119677562, 26.93605394127892, 16.202097655007115, 6.300877373928769, 389.7069815251737, 31.765554356555302, 7.140733162062209, 47.67390041423606, 47.811796044031475, 23.630395525902887, 40.20396416511805, 12.833807220158466, 22.4965966451724, 6.439573416346811, 185.7347891447746, 72.60552141514613, 19.53454465026214, 15.978238808239569, 13.063469704766954, 6.265258188180271, 34.1107527538008, 63.19854403276283, 6.616648342812059, 113.60911820641822, 6.566208184461597, 75.11241150499096, 60.77075102224196, 42.061599035983775, 46.103264391408565, 87.4182232995216, 141.07774795394573, 151.1320118820747, 481.2416278638292, 42.42607186288484, 34.345017876046555, 54.75459215206463, 72.45773191237308, 100.84666572906005, 23.549525770417777, 8.489577528931774, 24.551929978171806, 12.888766658763322, 49.68381203692404, 15.422285564937944, 125.65384138465679, 30.397966951469712, 27.270698420745653, 27.099048518054474, 156.8463684831965, 7.915983211087909, 187.60012400997834, 96.46589167967049, 71.89903790577651, 9.139245069329956, 127.91713457575737, 52.49748303051388, 27.90779716937122, 8.576253750437305, 65.16999107734897, 69.96772077381678, 43.934527506734725, 50.06803511995554, 42.41229435511866, 11.337792535260055, 92.1153742574229, 17.179003256449953, 23.41436647679125, 13.112990431266166, 64.63599102137371, 20.255723945897284, 75.35798583234478, 102.06125894702285, 18.580745060513834, 26.5161773152523, 23.313901383660426, 192.1034565216484, 46.511801533349434, 72.1324346430014, 75.45665074422266, 34.93598624567049, 79.79150228418723, 10.345334337714803, 38.01825854766095, 124.34573559411038, 43.4923312979913, 5.054635033502185, 49.767462485557246, 8.229730595017402, 68.04143004524683, 28.810070982190513, 21.83451257561652, 6.903511458030656, 16.736564749167666, 49.70731867312349, 15.308672873281234, 21.189865449223387, 70.34183429676627, 43.39422817619453, 70.40887313326093, 37.3121000955872, 47.596328662981335, 52.855897425360155, 44.192992087668685, 6.377305485498518, 79.68320926206587, 32.99096599656137, 209.18475897163356, 32.55436768167911, 41.30594383100017, 34.61974854568331, 18.951698090338887, 41.55385547276799, 57.223889772184755, 47.51120842333313, 123.41441812335937, 6.251562586175794, 39.33257676771784, 130.35141222254927, 135.69912137809462, 131.99242029176696, 70.12218105111806, 45.49992224806314, 11.77690941529406, 5.83595202666188, 9.237351289526366, 37.542151875532205, 11.765170166005976, 28.461890034519204, 5.828350662928304, 23.065277487636983, 8.198526763505894, 18.860360721560703, 105.36819942755818, 39.181707898487595, 65.26223606272829, 29.398127665671414, 16.499785571056965, 110.5558867030779, 19.64329093700004, 44.21781087804614, 71.3331350597563, 38.66510254643331, 5.839161097082865, 6.669106137391587, 326.37849793012384, 26.606931290969335, 19.86139772610958, 22.537071039979068, 88.37353029003185, 11.476423063622484, 22.261595573223605, 11.580915956600872, 9.85048496065002, 15.639935324056893, 12.280741534153266, 23.521632785283376, 20.465935167624316, 13.632020121429225, 26.39937274650152, 31.255116337012236, 145.08090353834717, 140.2831203992128, 66.71878119903745, 32.15599793364798, 25.804201896748015, 111.47960093246029, 40.642368267183244, 145.9095509207118, 28.88035101898782, 49.722527268790884, 6.872663031826666, 21.17337126008074, 15.90655224561645, 42.806297764758504, 14.328425763040805, 14.496560346618661, 17.673396803032606, 19.23121050859685, 24.897494331191595, 32.52686848933328, 80.14849275562675, 63.374577430291154, 76.91140856192244, 27.893827261715444, 63.465544434511216, 69.95905019568605, 26.2949500155077, 9.119469239515121, 5.557785809486354, 16.44533598579733, 54.88177842253924, 14.341521226574372, 251.0304916112949, 17.63651405487081, 139.08284683084594, 13.838233465237819, 17.46843995361195, 202.38405201465744, 53.53642578573328, 68.98051038271804, 15.747750093205322, 27.764991297744718, 30.29204731258981, 53.5370907479409, 25.85808266064117, 57.78752997500231, 173.52764229644686, 35.334995101282075, 27.37846126567299, 10.502481413938463, 62.186478613708296, 86.59540085177463, 328.95088187580143, 104.7686347599838, 42.075123306795035, 34.89136589201932, 62.56467527308959, 21.992481583954355, 52.384306875634444, 6.898234352915194, 53.054096533400184, 37.004748286637835, 57.48910312922262, 34.49272234952973, 40.655072862253874, 153.23832597706482, 168.47460004661588, 153.27092334282594, 14.374096624582021, 17.47522131062833, 143.58984890130242, 56.75249100435342, 111.69987917328646, 16.956956576584297, 26.53271236146316, 33.09078923924028, 71.77839587525273, 5.55127303476317, 5.445040312390084, 7.458309920179268, 60.53099909465615, 14.606821962128965, 289.4249005870328, 121.55731411637484, 25.88391820732727, 8.125076584329253, 18.85500784395236, 6.235379320506899, 38.67121915893703, 41.83502730627808, 5.173751914269115, 6.871221080242857, 86.80841677891489, 8.836705979776102, 23.809406110862543, 5.164961056748544, 150.78142407126924, 103.29076974180371, 38.90818641242245, 55.487481329300515, 33.982875039335845, 30.944194456726795, 13.296989899127201, 11.22075936182022, 5.745755269611203, 8.433031025626377, 225.20154284735193, 110.78413961479244, 186.7946777581286, 20.214631961667536, 21.639784958095436, 8.053923690355866, 31.57656468171556, 14.931975096864418, 16.4062582202943, 16.3396321839484, 31.467426023937065, 89.13562022617785, 139.5754481616569, 108.80033218047721, 11.475977284580075, 37.02669229079955, 25.01959486495189, 61.346937726608175, 31.112180757678907, 29.26337756784772, 12.15341106885423, 58.788523868141525, 41.82524109223579, 46.73677208802073, 29.49497932166495, 13.230050691870435, 70.49728900846466, 122.6009624982872, 10.175442916362238, 7.541080455753765, 125.19608917609159, 7.096837181987375, 34.57975638827786, 48.99719827741131, 20.136884995886334, 28.108610525583657, 47.54586000802074, 5.739948487806487, 288.8359434216485, 20.80451714897945, 22.580838688193786, 96.91243391844982, 26.530617682188282, 46.690585571946265, 45.27015687932401, 211.274087895192, 58.07436142327462, 25.48902365489932, 38.07328197113766, 28.779674924885253, 80.58950622076149, 6.83147885470799, 21.278179439616775, 78.39583959034182, 70.28697519892509, 218.7195807252162, 51.07795831732477, 23.46650137533203, 33.12777459382858, 8.804583687427268, 80.24445735634535, 13.644637212905879, 38.909463491319926, 61.848413750218775, 169.96402675900413, 56.42003923261203, 66.2534085845368, 46.5832559148559, 12.011346834521214, 43.61758057242759, 199.22996960407045, 6.580800752144334, 5.030624536721584, 122.46700338561193, 36.57596409469311, 8.741658173310547, 18.15357477040057, 96.05667750502765, 61.48159568342763, 24.39197526465823, 44.33280770260135, 84.72307413420737, 87.84117481893159, 42.96536766418469, 36.654240988005704, 5.327699324012194, 57.62838360659262, 111.35263958569983, 26.299305105467557, 9.936778144116541, 61.069909701123386, 50.48315224668051, 26.432696377515615, 21.956056230169764, 46.12688829050798, 192.232185543653, 53.98895843549289, 117.58406511752338, 47.157211422830784, 26.835375552712343, 34.05657477980806, 25.52157426455279, 19.858693098137227, 60.295374911785345, 17.6809760633518, 11.044201874916757, 62.01950067362743, 85.96438612221726, 19.940851847160328, 46.1364693115263, 126.10836485018214, 80.99269685157304, 190.76372824917453, 6.23722355704399, 76.77517684781239, 406.3867758872526, 101.45198344937525, 29.274656666931108, 26.768638631398904, 68.09728441020489, 115.18131491261298, 34.9214361382527, 19.373091706535668, 14.796044290538244, 8.712554247472482, 31.44831113221531, 51.967921962611, 88.87906553398393, 85.03526935983585, 51.51420000559395, 51.05064703489565, 63.55823802365108, 13.754074412960426, 22.325844656050656, 61.24893113655137, 35.268784703835294, 62.90902100836937, 39.629975388239416, 25.502124489530342, 59.02010039311995, 94.00603166354365, 84.41621604192699, 6.13023169476253, 40.27523507904435, 25.129034785069255, 9.128792590726567, 56.35050920148815, 6.575522898248985, 15.282615937155223, 20.82855333306505, 70.94604373267254, 127.19018506462595, 65.52241710444846, 25.550397126110965, 68.16639292822502, 56.603025582534755, 36.105150965970275, 97.26428449003525, 10.603097642580483, 78.9361266641682, 45.86227197308742, 5.324028114150204, 84.57367181752093, 64.04487260358093, 22.902032410225743, 16.69681658522037, 17.59842258252129, 75.52904057316563, 6.708320852032042, 64.05437088213958, 8.42227866180616, 24.755531231071302, 112.26872133137002, 296.1840148640123, 5.79017954419084, 18.52637662503884, 56.711067341812424, 5.085310704130503, 74.49077463939767, 25.697846124798815, 60.89298273185257, 64.60763011398305, 49.448039589481034, 27.342507289003276, 81.20909486182096, 60.13197453543258, 23.876764430505887, 18.844443423359284, 94.9558563009197, 56.74727530891187, 54.38466495621606, 156.03950602306492, 18.08013873755965, 9.403597119548579, 20.86511772100988, 48.612463731101826, 40.68873780723625, 13.177970389348966, 19.710136076767288, 33.02868895676784, 36.37187648588037, 5.197196798069252, 134.23758526834223, 5.194657368466044, 50.359583428338496, 6.130650320560632, 20.621866376959694, 7.827290815186451, 56.58968508150255, 27.838495016141508, 32.64924536459766, 57.10647258336654, 95.47082935762646, 65.0469683102413, 87.17718458077739, 143.45069910904334, 55.25145213608673, 8.075166702332739, 13.984158600434942, 28.163170454450643, 39.742279846785465, 16.447943599481203, 26.344085638217773, 7.9014542191118045, 33.6196909709104, 8.378331185072366, 12.025680777069136, 72.14012793695497, 26.296930324574895, 27.435283377491807, 24.655539542425203, 38.21687445841762, 49.04456986496268, 12.107499539665593, 10.415604643372966, 83.50004305336114, 22.172590894844113, 88.5271409191658, 9.259784857972262, 183.70537593482115, 104.36363148185106, 7.696088935860922, 39.52729238475217, 12.140113259435068, 16.91736667687468, 97.6734639670456, 9.875356756250646, 17.175779146125794, 6.035715432654942, 131.42689031845754, 13.29891426661904, 11.490635711230624, 26.231679848766746, 13.480727851579072, 72.92436678297061, 29.83397458290599, 93.3119589020483, 214.57346518298095, 38.792456305738035, 64.08595584596128, 6.231261770904037, 26.81023192452331, 52.37849212813234, 338.7829110772877, 80.64909782121588, 13.691368987871545, 58.784159352082575, 41.33406431604816, 60.22604450369563, 123.10683518221651, 68.60351925005239, 104.52289419445616, 20.204207069586374, 8.31501541416413, 37.85999763025312, 41.22648052905725, 38.77118722406425, 7.010181070636918, 43.968786026681556, 24.041916461648107, 30.730799642842165, 18.54029988339476, 12.008210614430977, 44.94171361638127, 30.52742582369338, 17.616801392060374, 59.86572233242134, 34.440894620651065, 22.871166407260183, 15.655471628935373, 17.372029099047424, 46.52590805969684, 73.01809164742689, ...])
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);
([2924346.875, 2944565.625, 2986623.4375, 3002956.25, 3147323.4375, 3397809.375, 3436142.1875, 3999750.0, 4792262.5, 5266640.625, 5416700.0, 5455485.9375, 5455493.75, 5456928.125, 5734228.125, 5736831.25, 5759906.25, 5788390.625, 5790696.875, 5795334.375, 5797781.25, 5839612.5, 5879600.0, 5883587.5, 5910326.5625, 5915831.25, 5917845.3125, 5924237.5, 5925865.625, 5926953.125, 5927317.1875, 5927843.75, 5929675.0, 5931595.3125, 5934132.8125, 5934760.9375, 5936929.6875, 5937835.9375, 5939276.5625, 5943960.9375, 5944914.0625, 5945410.9375, 5946307.8125, 5956275.0, 5956276.5625, 5958031.25, 5959417.1875, 5959895.3125, 5960125.0, 5960753.125, 5966514.0625, 5970970.3125, 5972035.9375, 5977026.5625, 5977026.5625, 5978792.1875, 5978793.75, 5979534.375, 5980356.25, 5980540.625, 5980615.625, 5980959.375, 5981001.5625, 5981025.0, 5981032.8125, 5981273.4375, 5981423.4375, 5981506.25, 5981570.3125, 5981609.375, 5981657.8125, 5981925.0, 5982148.4375, 5982321.875, 5982434.375, 5982484.375, 5982528.125, 5982728.125, 5983053.125, 5983329.6875, 5983426.5625, 5983896.875, 5983953.125, 5984093.75, 5984282.8125, 5984553.125, 5985385.9375, 5986021.875, 5986590.625, 5990492.1875, 5990950.0, 5990971.875, 6009821.875, 6015390.625, 6023539.0625, 6024298.4375, 6054309.375, 6055198.4375, 6058721.875, 6059100.0, 6063656.25, 6066542.1875, 6078581.25, 6084646.875, 6084657.8125, 6088356.25, 6091434.375, 6093745.3125, 6093850.0, 6096214.0625, 6098381.25, 6099982.8125, 6100262.5, 6124200.0, 6126929.6875, 6128434.375, 6128501.5625, 6128693.75, 6132237.5, 6132362.5, 6132404.6875, 6132954.6875, 6134525.0, 6138635.9375, 6142756.25, 6170331.25, 6175581.25, 6175920.3125, 6175925.0, 6176109.375, 6176148.4375, 6176251.5625, 6176275.0, 6177120.3125, 6177282.8125, 6177342.1875, 6177582.8125, 6177603.125, 6178493.75, 6179835.9375, 6180748.4375, 6180907.8125, 6181957.8125, 6182935.9375, 6183243.75, 6183395.3125, 6184270.3125, 6185128.125, 6186509.375, 6194025.0, 6209742.1875, 6210400.0, 6217031.25, 6218221.875, 6262400.0, 6267309.375, 6271317.1875, 6285198.4375, 6286095.3125, 6287734.375, 6291410.9375, 6292990.625, 6298095.3125, 6301600.0, 6303882.8125, 6306959.375, 6336956.25, 6338228.125, 6341228.125, 6342273.4375, 6343179.6875, 6343450.0, 6343943.75, 6344215.625, 6344232.8125, 6344571.875, 6344598.4375, 6344617.1875, 6344620.3125, 6344890.625, 6345023.4375, 6345162.5, 6345195.3125, 6345231.25, 6345573.4375, 6345728.125, 6345757.8125, 6345910.9375, 6346014.0625, 6346209.375, 6346900.0, 6347385.9375, 6347401.5625, 6347420.3125, 6347785.9375, 6348076.5625, 6348728.125, 6348964.0625, 6349156.25, 6349192.1875, 6349276.5625, 6349921.875, 6350020.3125, 6351696.875, 6352292.1875, 6352331.25, 6352590.625, 6352612.5, 6353323.4375, 6355182.8125, 6355465.625, 6355529.6875, 6356089.0625, 6356146.875, 6356576.5625, 6356970.3125, 6358621.875, 6359306.25, 6359323.4375, 6360632.8125, 6362714.0625, 6372925.0, 6373746.875, 6377125.0, 6387075.0, 6387101.5625, 6395089.0625, 6447707.8125, 6481442.1875, 6483025.0, 6484126.5625, 6485400.0, 6486600.0, 6487826.5625, 6495279.6875, 6496150.0, 6496159.375, 6496362.5, 6496417.1875, 6496687.5, 6497096.875, 6497334.375, 6497967.1875, 6498026.5625, 6498237.5, 6498279.6875, 6498581.25, 6498909.375, 6499729.6875, 6499948.4375, 6500839.0625, 6508262.5, 6512271.875, 6514573.4375, 6516168.75, 6516629.6875, 6518032.8125, 6523879.6875, 6527639.0625, 6529360.9375, 6541267.1875, 6542935.9375, 6544412.5, 6545184.375, 6547835.9375, 6548075.0, 6550118.75, 6550818.75, 6551314.0625, 6551706.25, 6551853.125, 6552904.6875, 6553110.9375, 6553231.25, 6553465.625, 6553667.1875, 6553776.5625, 6554157.8125, 6554596.875, 6554623.4375, 6554843.75, 6555184.375, 6555467.1875, 6556028.125, 6556085.9375, 6556203.125, 6556267.1875, 6556943.75, 6556965.625, 6558034.375, 6562751.5625, 6566706.25, 6567043.75, 6581642.1875, 6584039.0625, 6587303.125, 6588934.375, 6590175.0, 6590446.875, 6593926.5625, 6596121.875, 6598346.875, 6598392.1875, 6599023.4375, 6599153.125, 6599817.1875, 6600389.0625, 6600815.625, 6601582.8125, 6607121.875, 6616378.125, 6616556.25, 6616979.6875, 6617123.4375, 6618862.5, 6623275.0, 6642670.3125, 6647689.0625, 6652178.125, 6674725.0, 6683126.5625, 6695668.75, 6703696.875, 6704498.4375, 6719276.5625, 6719456.25, 6719482.8125, 6719523.4375, 6719896.875, 6720107.8125, 6722787.5, 6725245.3125, 6725745.3125, 6726210.9375, 6729340.625, 6730289.0625, 6733443.75, 6733454.6875, 6733928.125, 6734906.25, 6734996.875, 6735070.3125, 6736481.25, 6737221.875, 6737251.5625, 6737345.3125, 6737987.5, 6742279.6875, 6745264.0625, 6745907.8125, 6746176.5625, 6747793.75, 6747860.9375, 6749214.0625, 6750520.3125, 6751846.875, 6751926.5625, 6752326.5625, 6753300.0, 6753306.25, 6759739.0625, 6763539.0625, 6765523.4375, 6767343.75, 6768370.3125, 6769764.0625, 6769864.0625, 6770679.6875, 6776676.5625, 6778278.125, 6779393.75, 6779853.125, 6781718.75, 6781939.0625, 6782900.0, 6783450.0, 6784067.1875, 6785901.5625, 6785929.6875, 6787056.25, 6787090.625, 6787775.0, 6788421.875, 6788440.625, 6788642.1875, 6789384.375, 6789495.3125, 6789673.4375, 6790017.1875, 6791360.9375, 6791470.3125, 6791640.625, 6791792.1875, 6791828.125, 6791956.25, 6792179.6875, 6792260.9375, 6792781.25, 6793235.9375, 6793946.875, 6794045.3125, 6795615.625, 6809879.6875, 6811578.125, 6812414.0625, 6813021.875, 6815412.5, 6819867.1875, 6820987.5, 6821856.25, 6821895.3125, 6822035.9375, 6823390.625, 6824665.625, 6824876.5625, 6825195.3125, 6826151.5625, 6826357.8125, 6826454.6875, 6826839.0625, 6826850.0, 6827001.5625, 6827818.75, 6827881.25, 6828200.0, 6828453.125, 6828470.3125, 6829500.0, 6831670.3125, 6832375.0, 6833915.625, 6834368.75, 6836112.5, 6840454.6875, 6840464.0625, 6840509.375, 6840510.9375, 6840934.375, 6840937.5, 6841137.5, 6841140.625, 6841253.125, 6842260.9375, 6844643.75, 6847510.9375, 6848787.5, 6851209.375, 6852473.4375, 6853573.4375, 6853928.125, 6858198.4375, 6858671.875, 6859639.0625, 6859907.8125, 6862125.0, 6863965.625, 6864787.5, 6864960.9375, 6865084.375, 6865190.625, 6865206.25, 6865842.1875, 6868729.6875, 6869650.0, 6873051.5625, 6873943.75, 6874078.125, 6876496.875, 6876517.1875, 6878250.0, 6879478.125, 6879828.125, 6882109.375, 6883210.9375, 6893704.6875, 6903237.5, 6907682.8125, 6909642.1875, 6910242.1875, 6910378.125, 6910415.625, 6910529.6875, 6914932.8125, 6915995.3125, 6916112.5, 6916350.0, 6916507.8125, 6917054.6875, 6917259.375, 6917795.3125, 6918117.1875, 6918154.6875, 6918909.375, 6918934.375, 6918957.8125, 6919028.125, 6919071.875, 6919606.25, 6919785.9375, 6919873.4375, 6919948.4375, 6919978.125, 6920162.5, 6920246.875, 6920254.6875, 6920275.0, 6920364.0625, 6920423.4375, 6920693.75, 6920917.1875, 6920948.4375, 6921317.1875, 6921323.4375, 6921364.0625, 6921839.0625, 6921879.6875, 6921887.5, 6921920.3125, 6921942.1875, 6921964.0625, 6921968.75, 6922018.75, 6922120.3125, 6922253.125, 6922284.375, 6922292.1875, 6922354.6875, 6922448.4375, 6922482.8125, 6922507.8125, 6922600.0, 6922673.4375, 6922675.0, 6922767.1875, 6922793.75, 6922796.875, 6922907.8125, 6922909.375, 6922978.125, 6923031.25, 6923357.8125, 6923382.8125, 6923385.9375, 6923448.4375, 6923642.1875, 6923706.25, 6923723.4375, 6923778.125, 6923943.75, 6924182.8125, 6924273.4375, 6924296.875, 6924346.875, 6924639.0625, 6924685.9375, 6924710.9375, 6924765.625, 6925482.8125, 6925582.8125, 6925870.3125, 6926137.5, 6926567.1875, 6926606.25, 6926610.9375, 6926665.625, 6926670.3125, 6926801.5625, 6926859.375, 6926965.625, 6927110.9375, 6927850.0, 6928290.625, 6928543.75, 6928693.75, 6929810.9375, 6929990.625, 6931276.5625, 6931575.0, 6931692.1875, 6931715.625, 6931985.9375, 6932751.5625, 6933606.25, 6934704.6875, 6934714.0625, 6936109.375, 6936382.8125, 6937206.25, 6937353.125, 6937642.1875, 6937823.4375, 6938221.875, 6938401.5625, 6938945.3125, 6940976.5625, 6941468.75, 6942489.0625, 6943125.0, 6943398.4375, 6943618.75, 6943878.125, 6944743.75, 6945048.4375, 6945275.0, 6945312.5, 6945375.0, 6945657.8125, 6946314.0625, 6946690.625, 6946729.6875, 6948473.4375, 6950375.0, 6950990.625, 6951054.6875, 6953428.125, 6953579.6875, 6954956.25, 6955893.75, 6956270.3125, 6957575.0, 6958329.6875, 6958557.8125, 6958882.8125, 6961229.6875, 6961531.25, 6962248.4375, 6963715.625, 6966657.8125, 6972007.8125, 6972071.875, 6972503.125, 6972506.25, 6973884.375, 6976331.25, 6977882.8125, 6980239.0625, 6981967.1875, 6982017.1875, 6985701.5625, 6986907.8125, 6996946.875, 7003584.375, 7021931.25, 7022676.5625, 7024148.4375, 7024573.4375, 7024801.5625, 7026842.1875, 7027198.4375, 7027773.4375, 7028279.6875, 7029885.9375, 7030929.6875, 7031110.9375, 7040485.9375, 7041504.6875, 7043698.4375, 7043746.875, 7044601.5625, 7050225.0, 7050925.0, 7053329.6875, 7054293.75, 7054600.0, 7056348.4375, 7058498.4375, 7058601.5625, 7060473.4375, 7062290.625, 7063640.625, 7065537.5, 7066434.375, 7066515.625, 7067143.75, 7068031.25, 7068701.5625, 7068787.5, 7069237.5, 7069501.5625, 7069520.3125, 7069556.25, 7069771.875, 7072853.125, 7074782.8125, 7075415.625, 7076504.6875, 7077718.75, 7079539.0625, 7081946.875, 7082650.0, 7083671.875, 7083726.5625, 7084162.5, 7084609.375, 7087012.5, 7087178.125, 7087289.0625, 7108465.625, 7109739.0625, 7110051.5625, 7112256.25, 7118868.75, 7120237.5, 7120610.9375, 7121957.8125, 7122028.125, 7122396.875, 7126531.25, 7126537.5, 7127698.4375, 7127884.375, 7127981.25, 7128231.25, 7128260.9375, 7128593.75, 7131671.875, 7132656.25, 7134743.75, 7134743.75, 7138676.5625, 7142254.6875, 7144043.75, 7144675.0, 7144685.9375, 7145096.875, 7145753.125, 7146346.875, 7146564.0625, 7146718.75, 7146743.75, 7147043.75, 7147646.875, 7148246.875, 7148701.5625, 7149018.75, 7149393.75, 7150721.875, 7150729.6875, 7151068.75, 7151268.75, 7152187.5, 7152265.625, 7152637.5, 7152928.125, 7153235.9375, 7154014.0625, 7154770.3125, 7156054.6875, 7156275.0, 7156303.125, 7157048.4375, 7157450.0, 7157885.9375, 7158431.25, 7158657.8125, 7159223.4375, 7159434.375, 7159506.25, 7159804.6875, 7160067.1875, 7160132.8125, 7160151.5625, 7160165.625, 7160301.5625, 7160365.625, 7160467.1875, 7160656.25, 7160710.9375, 7160714.0625, 7160921.875, 7161073.4375, 7161112.5, 7161284.375, 7161420.3125, 7161517.1875, 7161712.5, 7162232.8125, 7162370.3125, 7162492.1875, 7162581.25, 7162696.875, 7162795.3125, 7162846.875, 7163475.0, 7165593.75, 7167217.1875, 7167607.8125, 7167757.8125, 7169803.125, 7170253.125, 7170739.0625, 7176346.875, 7176409.375, 7176529.6875, 7176837.5, 7178376.5625, 7179979.6875, 7188281.25, 7192954.6875, 7196901.5625, 7205226.5625, 7206895.3125, 7208975.0, 7209975.0, 7210984.375, 7211687.5, 7212275.0, 7213364.0625, 7213571.875, 7214796.875, 7215300.0, 7216334.375, 7216740.625, 7216760.9375, 7216767.1875, 7217109.375, 7218109.375, 7218309.375, 7218712.5, 7219104.6875, 7219184.375, 7221693.75, 7221932.8125, 7222025.0, 7222228.125, 7222354.6875, 7222914.0625, 7223493.75, 7223623.4375, 7223759.375, 7224009.375, 7224017.1875, 7224107.8125, 7224135.9375, 7224290.625, 7225051.5625, 7225400.0, 7226271.875, 7226517.1875, 7226937.5, 7227267.1875, 7227476.5625, 7228623.4375, 7228885.9375, 7229714.0625, 7230523.4375, 7231045.3125, 7231375.0, 7231400.0, 7231829.6875, 7231848.4375, 7232993.75, 7233007.8125, 7233868.75, 7233898.4375, 7234257.8125, 7234362.5, 7235712.5, 7237325.0, 7238143.75, 7238920.3125, 7245321.875, 7247782.8125, 7249503.125, 7249884.375, 7250350.0, 7250675.0, 7250742.1875, 7250753.125, 7253918.75, 7254842.1875, 7258203.125, 7259095.3125, 7260581.25, 7260595.3125, 7263329.6875, 7264870.3125, 7267870.3125, 7270634.375, 7270667.1875, 7271054.6875, 7275165.625, 7287562.5, 7295409.375, 7295478.125, 7307853.125, 7312629.6875, 7314253.125, 7315639.0625, 7316850.0, 7316942.1875, 7317525.0, 7319484.375, 7320073.4375, 7322621.875, 7322870.3125, 7322926.5625, 7324837.5, 7325293.75, 7326329.6875, 7326343.75, 7326862.5, 7327868.75, 7329995.3125, 7330507.8125, 7330781.25, 7330864.0625, 7330965.625, 7330970.3125, 7331039.0625, 7331103.125, 7331132.8125, 7331268.75, 7331450.0, 7331490.625, 7331496.875, 7331556.25, 7331595.3125, 7331656.25, 7331831.25, 7331931.25, 7332026.5625, 7332139.0625, 7332150.0, 7332170.3125, 7332189.0625, 7332245.3125, 7332296.875, 7332376.5625, 7332570.3125, 7332795.3125, 7333228.125, 7333245.3125, 7333573.4375, 7349715.625, 7404610.9375, 7422159.375, 7423446.875, 7424810.9375, 7432853.125, 7433401.5625, 7444418.75, 7458937.5, 7459159.375, 7499614.0625, 7507142.1875, 7514735.9375, 7516075.0, 7526598.4375, 7530245.3125, 7533860.9375, 7537785.9375, 7539229.6875, 7546154.6875, 7547800.0, 7547820.3125, 7548732.8125, 7552317.1875, 7558478.125, 7558610.9375, 7563135.9375, 7563310.9375, 7566651.5625, 7566848.4375, 7582910.9375, 7586281.25, 7587695.3125, 7588598.4375, 7589415.625, 7589984.375, 7591601.5625, 7605498.4375, 7618710.9375, 7620423.4375, 7621089.0625, 7621237.5, 7621331.25, 7621521.875, 7622914.0625, 7624485.9375, 7624565.625, 7624603.125, 7624728.125, 7625617.1875, 7627723.4375, 7630253.125, 7630934.375, 7631220.3125, 7632337.5, 7688767.1875, 7689118.75, 7691126.5625, 7693764.0625, 7694540.625, 7697353.125, 7697673.4375, 7698129.6875, 7698875.0, 7702300.0, 7702573.4375, ...], [10.713051006161061, 46.205314157275645, 8.201989650568361, 39.191171277572906, 73.81780751787149, 16.290522894242418, 11.639361768222534, 9.738681771399184, 5.9691967392463745, 6.732654842965873, 29.352996094753685, 40.12125904893643, 12.017598496462186, 68.34107171205457, 37.147371835394026, 94.99008762392955, 80.77086210229642, 49.65768933980893, 51.86516697791328, 8.451240349931041, 22.88988921125111, 25.353349388135815, 16.118568391287695, 189.04782623262463, 79.7423608681323, 14.734670586719915, 15.748960378227371, 69.25263500929465, 33.96008278542762, 18.718211988398764, 11.07439989541346, 59.49617827599503, 31.518008305291136, 21.048199636278362, 83.00129900642807, 28.759538032206706, 47.89007404220147, 25.234124368296882, 5.203005334323263, 69.95991192049884, 12.353724370063453, 52.0919660819248, 65.9823739755005, 17.073080500997563, 9.330219102883145, 162.02498344023223, 21.573660930648437, 9.578592551035683, 45.58965015493776, 5.919388677843027, 20.3561706082224, 57.65207190589465, 9.592643509130072, 19.757382372328532, 65.12601775090488, 17.362042864366654, 23.134872640463982, 13.25350392272108, 14.741929745440705, 34.22820653433416, 5.896698899301772, 107.76521471397601, 70.70603547451226, 59.029288992108995, 24.00099199192772, 132.08803994529163, 24.140164981141535, 18.29872189918449, 142.06813615802804, 12.073728390312809, 31.390276475147832, 11.261277487943248, 8.418631101847945, 56.18341850578175, 38.867383485099715, 32.15888864659335, 7.751113626922416, 9.663634067209365, 77.12737849232339, 28.988637712743635, 15.517399424267868, 40.397779223349886, 98.97913608797337, 16.700759144972917, 7.09302357784072, 81.95291949975585, 13.134130942625509, 25.140228873830974, 15.603974882966293, 136.36749136434074, 51.12327276518053, 31.85876848629799, 18.884702434330613, 40.79936324334394, 22.10295502501878, 15.157816701191905, 37.634204285703255, 17.00835975278213, 6.601492450544085, 21.783283143651737, 121.20824694432605, 40.91605215524861, 14.163149939540903, 59.87050462976426, 25.903360670210766, 135.3827979296443, 15.440225188024213, 55.26969908813271, 38.38513297567471, 30.966463387672768, 44.336218752106795, 7.400767357764462, 10.425683430908393, 120.02755428289045, 9.707186410030868, 55.64219926303227, 93.75080651201856, 26.256721822094605, 48.52290975704314, 264.47327012391105, 32.149619814385154, 16.595032312372005, 55.10612886869767, 30.482279644392683, 6.738721400627129, 17.14912541764162, 31.036769729117083, 57.67941455871323, 32.54465771788227, 18.73624352390624, 48.473051302894035, 5.124025258869424, 38.24449133353264, 99.12449408006285, 23.77179586760235, 17.196367781490746, 19.911650030215483, 10.064222237947954, 73.4298316205631, 11.705310106359512, 9.407745729499934, 32.96836339263979, 50.949771078810436, 12.942495050154493, 48.43369718198619, 15.323182683336192, 17.70932693795569, 13.08826764799719, 55.167752461682426, 54.798207918968025, 7.101952636019721, 14.430295221770141, 5.851312221704075, 18.53686804482017, 13.294598980874488, 54.064717740056295, 45.18571726350412, 68.88784114637386, 69.26914843107113, 33.093727749161204, 31.724061398414996, 118.1231016667246, 47.017646445387115, 8.837330458194552, 32.98303372226831, 7.641316819198258, 22.192540683773082, 56.05726272069935, 54.02239980189213, 82.56865532363777, 19.772048089181393, 55.792751228464994, 83.05263997823434, 11.280276607121538, 67.12335933863798, 40.31528544501474, 82.24858064103832, 50.43358457325405, 22.59420590158868, 131.01569162943122, 20.047218641106586, 71.7253569837988, 15.739605471539642, 5.134116227523508, 17.837474786412624, 47.78823245688457, 7.525876678995922, 15.96981254426212, 77.54406542175451, 33.132638452788974, 13.121468741569801, 32.901354497422886, 11.39277828244078, 12.546123099899562, 9.023400581507552, 30.31271131870807, 150.78537421664808, 25.534395563008463, 13.674889303847308, 20.91758858731326, 5.052926636921982, 29.507661278011287, 11.394459606525716, 155.34331456051294, 42.64896115966768, 15.229340705011019, 9.203719082034143, 10.539028624710923, 17.103138739597494, 22.135482322808517, 19.800250168423556, 51.4368507962477, 19.868799525911314, 10.954782941387526, 22.475431895625636, 28.548947615854463, 12.145965228984089, 74.74050633683366, 29.566695151328346, 6.946035610338883, 142.1554719517788, 24.78160114181562, 45.85085389639405, 34.358122359079275, 79.92529124611775, 20.46638653961609, 9.906649331132638, 19.81841686857041, 5.726247656544598, 24.54430034455945, 6.379397206152645, 40.510494070799304, 51.281405098744315, 20.680554179442353, 296.1861297905942, 19.05118378806145, 68.51594384504914, 28.355964016934664, 5.350764833392888, 92.17253780908504, 45.49380559753511, 18.73108758539995, 205.4873664774685, 120.89439055706512, 34.69866804561807, 44.12746711471958, 12.917345050915397, 24.67971460869807, 8.626418306704283, 191.42206980347584, 14.039874578054588, 52.07264398030967, 7.6894169588307815, 37.85302745629882, 24.838355098441056, 22.75420992262668, 13.437049805430775, 9.268988559588795, 14.124137656345068, 18.109515646692206, 16.175644694466502, 12.894940220278274, 66.68462418310563, 5.669314021639514, 11.447520965683479, 13.213117412481651, 51.085618237203676, 20.18919466302434, 99.28238726810977, 14.586143524237917, 283.1449188438106, 18.001978262989546, 77.64648024003074, 26.226342544448393, 24.99042037967624, 313.19531193083225, 141.0131251480738, 27.857801846489703, 18.234903529750422, 16.316697374267136, 231.4565092097347, 179.81987208164455, 28.28995158924046, 54.514717849607706, 40.52120878621977, 28.797343190549597, 27.48995766890571, 74.13655450852944, 45.25076554609947, 55.99258027750915, 47.77164415650097, 36.63632405448321, 32.242423537409266, 23.522335622840085, 52.83308057403912, 328.32384395425595, 11.939908875386346, 14.136301356810709, 58.37861561479734, 79.26308643772694, 70.3351997959965, 49.626451791977395, 46.958960186247666, 110.96451829340843, 40.381943307575874, 22.40470757081033, 10.47793037456153, 29.682704915094867, 93.68212454764652, 50.85295863018345, 56.456409672097834, 95.84940216200175, 14.95015360259729, 379.70048444526566, 6.529676474082922, 32.27906559007781, 82.92487483240097, 87.96660119122541, 12.762061189935608, 59.41596727412542, 36.15318989393144, 30.41785250202821, 68.63049143630359, 33.10722025645674, 22.57892658336801, 17.70560157250341, 7.625062255641228, 8.26166392855743, 25.807839237669526, 86.50703294937978, 27.60413260739759, 111.51652688823657, 8.061709977128482, 14.899546882613032, 67.63536240221649, 6.8062329426472346, 42.982243879826186, 20.070293216136466, 6.527426905694235, 37.95197918099651, 77.30050769500922, 22.375874523374698, 58.24987553488152, 35.84912053376168, 21.16291869976248, 5.03502414482645, 26.155160213158435, 109.96334972003814, 7.62634753375584, 116.07139562465346, 38.828670279264045, 6.029506614058144, 42.55839983747294, 56.565471128398514, 20.106555526303847, 33.77668339732171, 61.79541534894802, 8.585750934981327, 15.0622448201121, 6.000270182131832, 8.420863653511608, 7.155596986445784, 25.645215785354807, 26.121538511683315, 6.793379708982833, 50.1641123465557, 5.127070269467935, 12.330319894325836, 31.687069188073227, 122.11994585244659, 60.545315244154764, 76.76749326871932, 50.685468111187674, 143.35267532111848, 39.38121490414949, 37.02799323642322, 135.979358651345, 33.97263260345547, 32.794849967143264, 27.154328943240433, 20.925360369176627, 34.50314110537914, 37.475828074527946, 28.056586067061883, 28.066003356594525, 24.017830000825622, 27.32049126712908, 21.71515180283756, 41.97998668558744, 34.79487773341015, 35.420491214506185, 6.305705988921438, 78.397551373021, 70.50643435520176, 11.034800962326024, 10.668849590083159, 79.272262314767, 16.24944480851766, 196.76980416042454, 21.46195123958132, 111.24738326437202, 46.64736898553465, 183.57719852993728, 113.66004491997514, 97.97878328659418, 31.68111211200847, 12.044697999352595, 11.106111628398567, 66.32695527335767, 25.59211031956925, 70.09018695575686, 28.710383871100476, 14.868243461216393, 6.000420592400463, 24.198040164545528, 57.1730844450471, 184.90810273916534, 20.318946394143953, 36.68504803572285, 44.294204931928164, 53.1834723229263, 104.49827218279741, 115.73453240930353, 12.612518523384017, 64.51194759560335, 81.96107970734764, 255.96977045809686, 49.38730799671227, 54.19476850669163, 231.82897518128703, 36.54336095702415, 13.026670412491947, 10.417742924261002, 181.37464157412927, 58.259611082480276, 126.26774134640492, 84.37277195539801, 24.69831180324078, 37.008529070798524, 65.66973085176188, 36.98518035476906, 52.51536119677562, 26.93605394127892, 16.202097655007115, 6.300877373928769, 389.7069815251737, 31.765554356555302, 7.140733162062209, 47.67390041423606, 47.811796044031475, 23.630395525902887, 40.20396416511805, 12.833807220158466, 22.4965966451724, 6.439573416346811, 185.7347891447746, 72.60552141514613, 19.53454465026214, 15.978238808239569, 13.063469704766954, 6.265258188180271, 34.1107527538008, 63.19854403276283, 6.616648342812059, 113.60911820641822, 6.566208184461597, 75.11241150499096, 60.77075102224196, 42.061599035983775, 46.103264391408565, 87.4182232995216, 141.07774795394573, 151.1320118820747, 481.2416278638292, 42.42607186288484, 34.345017876046555, 54.75459215206463, 72.45773191237308, 100.84666572906005, 23.549525770417777, 8.489577528931774, 24.551929978171806, 12.888766658763322, 49.68381203692404, 15.422285564937944, 125.65384138465679, 30.397966951469712, 27.270698420745653, 27.099048518054474, 156.8463684831965, 7.915983211087909, 187.60012400997834, 96.46589167967049, 71.89903790577651, 9.139245069329956, 127.91713457575737, 52.49748303051388, 27.90779716937122, 8.576253750437305, 65.16999107734897, 69.96772077381678, 43.934527506734725, 50.06803511995554, 42.41229435511866, 11.337792535260055, 92.1153742574229, 17.179003256449953, 23.41436647679125, 13.112990431266166, 64.63599102137371, 20.255723945897284, 75.35798583234478, 102.06125894702285, 18.580745060513834, 26.5161773152523, 23.313901383660426, 192.1034565216484, 46.511801533349434, 72.1324346430014, 75.45665074422266, 34.93598624567049, 79.79150228418723, 10.345334337714803, 38.01825854766095, 124.34573559411038, 43.4923312979913, 5.054635033502185, 49.767462485557246, 8.229730595017402, 68.04143004524683, 28.810070982190513, 21.83451257561652, 6.903511458030656, 16.736564749167666, 49.70731867312349, 15.308672873281234, 21.189865449223387, 70.34183429676627, 43.39422817619453, 70.40887313326093, 37.3121000955872, 47.596328662981335, 52.855897425360155, 44.192992087668685, 6.377305485498518, 79.68320926206587, 32.99096599656137, 209.18475897163356, 32.55436768167911, 41.30594383100017, 34.61974854568331, 18.951698090338887, 41.55385547276799, 57.223889772184755, 47.51120842333313, 123.41441812335937, 6.251562586175794, 39.33257676771784, 130.35141222254927, 135.69912137809462, 131.99242029176696, 70.12218105111806, 45.49992224806314, 11.77690941529406, 5.83595202666188, 9.237351289526366, 37.542151875532205, 11.765170166005976, 28.461890034519204, 5.828350662928304, 23.065277487636983, 8.198526763505894, 18.860360721560703, 105.36819942755818, 39.181707898487595, 65.26223606272829, 29.398127665671414, 16.499785571056965, 110.5558867030779, 19.64329093700004, 44.21781087804614, 71.3331350597563, 38.66510254643331, 5.839161097082865, 6.669106137391587, 326.37849793012384, 26.606931290969335, 19.86139772610958, 22.537071039979068, 88.37353029003185, 11.476423063622484, 22.261595573223605, 11.580915956600872, 9.85048496065002, 15.639935324056893, 12.280741534153266, 23.521632785283376, 20.465935167624316, 13.632020121429225, 26.39937274650152, 31.255116337012236, 145.08090353834717, 140.2831203992128, 66.71878119903745, 32.15599793364798, 25.804201896748015, 111.47960093246029, 40.642368267183244, 145.9095509207118, 28.88035101898782, 49.722527268790884, 6.872663031826666, 21.17337126008074, 15.90655224561645, 42.806297764758504, 14.328425763040805, 14.496560346618661, 17.673396803032606, 19.23121050859685, 24.897494331191595, 32.52686848933328, 80.14849275562675, 63.374577430291154, 76.91140856192244, 27.893827261715444, 63.465544434511216, 69.95905019568605, 26.2949500155077, 9.119469239515121, 5.557785809486354, 16.44533598579733, 54.88177842253924, 14.341521226574372, 251.0304916112949, 17.63651405487081, 139.08284683084594, 13.838233465237819, 17.46843995361195, 202.38405201465744, 53.53642578573328, 68.98051038271804, 15.747750093205322, 27.764991297744718, 30.29204731258981, 53.5370907479409, 25.85808266064117, 57.78752997500231, 173.52764229644686, 35.334995101282075, 27.37846126567299, 10.502481413938463, 62.186478613708296, 86.59540085177463, 328.95088187580143, 104.7686347599838, 42.075123306795035, 34.89136589201932, 62.56467527308959, 21.992481583954355, 52.384306875634444, 6.898234352915194, 53.054096533400184, 37.004748286637835, 57.48910312922262, 34.49272234952973, 40.655072862253874, 153.23832597706482, 168.47460004661588, 153.27092334282594, 14.374096624582021, 17.47522131062833, 143.58984890130242, 56.75249100435342, 111.69987917328646, 16.956956576584297, 26.53271236146316, 33.09078923924028, 71.77839587525273, 5.55127303476317, 5.445040312390084, 7.458309920179268, 60.53099909465615, 14.606821962128965, 289.4249005870328, 121.55731411637484, 25.88391820732727, 8.125076584329253, 18.85500784395236, 6.235379320506899, 38.67121915893703, 41.83502730627808, 5.173751914269115, 6.871221080242857, 86.80841677891489, 8.836705979776102, 23.809406110862543, 5.164961056748544, 150.78142407126924, 103.29076974180371, 38.90818641242245, 55.487481329300515, 33.982875039335845, 30.944194456726795, 13.296989899127201, 11.22075936182022, 5.745755269611203, 8.433031025626377, 225.20154284735193, 110.78413961479244, 186.7946777581286, 20.214631961667536, 21.639784958095436, 8.053923690355866, 31.57656468171556, 14.931975096864418, 16.4062582202943, 16.3396321839484, 31.467426023937065, 89.13562022617785, 139.5754481616569, 108.80033218047721, 11.475977284580075, 37.02669229079955, 25.01959486495189, 61.346937726608175, 31.112180757678907, 29.26337756784772, 12.15341106885423, 58.788523868141525, 41.82524109223579, 46.73677208802073, 29.49497932166495, 13.230050691870435, 70.49728900846466, 122.6009624982872, 10.175442916362238, 7.541080455753765, 125.19608917609159, 7.096837181987375, 34.57975638827786, 48.99719827741131, 20.136884995886334, 28.108610525583657, 47.54586000802074, 5.739948487806487, 288.8359434216485, 20.80451714897945, 22.580838688193786, 96.91243391844982, 26.530617682188282, 46.690585571946265, 45.27015687932401, 211.274087895192, 58.07436142327462, 25.48902365489932, 38.07328197113766, 28.779674924885253, 80.58950622076149, 6.83147885470799, 21.278179439616775, 78.39583959034182, 70.28697519892509, 218.7195807252162, 51.07795831732477, 23.46650137533203, 33.12777459382858, 8.804583687427268, 80.24445735634535, 13.644637212905879, 38.909463491319926, 61.848413750218775, 169.96402675900413, 56.42003923261203, 66.2534085845368, 46.5832559148559, 12.011346834521214, 43.61758057242759, 199.22996960407045, 6.580800752144334, 5.030624536721584, 122.46700338561193, 36.57596409469311, 8.741658173310547, 18.15357477040057, 96.05667750502765, 61.48159568342763, 24.39197526465823, 44.33280770260135, 84.72307413420737, 87.84117481893159, 42.96536766418469, 36.654240988005704, 5.327699324012194, 57.62838360659262, 111.35263958569983, 26.299305105467557, 9.936778144116541, 61.069909701123386, 50.48315224668051, 26.432696377515615, 21.956056230169764, 46.12688829050798, 192.232185543653, 53.98895843549289, 117.58406511752338, 47.157211422830784, 26.835375552712343, 34.05657477980806, 25.52157426455279, 19.858693098137227, 60.295374911785345, 17.6809760633518, 11.044201874916757, 62.01950067362743, 85.96438612221726, 19.940851847160328, 46.1364693115263, 126.10836485018214, 80.99269685157304, 190.76372824917453, 6.23722355704399, 76.77517684781239, 406.3867758872526, 101.45198344937525, 29.274656666931108, 26.768638631398904, 68.09728441020489, 115.18131491261298, 34.9214361382527, 19.373091706535668, 14.796044290538244, 8.712554247472482, 31.44831113221531, 51.967921962611, 88.87906553398393, 85.03526935983585, 51.51420000559395, 51.05064703489565, 63.55823802365108, 13.754074412960426, 22.325844656050656, 61.24893113655137, 35.268784703835294, 62.90902100836937, 39.629975388239416, 25.502124489530342, 59.02010039311995, 94.00603166354365, 84.41621604192699, 6.13023169476253, 40.27523507904435, 25.129034785069255, 9.128792590726567, 56.35050920148815, 6.575522898248985, 15.282615937155223, 20.82855333306505, 70.94604373267254, 127.19018506462595, 65.52241710444846, 25.550397126110965, 68.16639292822502, 56.603025582534755, 36.105150965970275, 97.26428449003525, 10.603097642580483, 78.9361266641682, 45.86227197308742, 5.324028114150204, 84.57367181752093, 64.04487260358093, 22.902032410225743, 16.69681658522037, 17.59842258252129, 75.52904057316563, 6.708320852032042, 64.05437088213958, 8.42227866180616, 24.755531231071302, 112.26872133137002, 296.1840148640123, 5.79017954419084, 18.52637662503884, 56.711067341812424, 5.085310704130503, 74.49077463939767, 25.697846124798815, 60.89298273185257, 64.60763011398305, 49.448039589481034, 27.342507289003276, 81.20909486182096, 60.13197453543258, 23.876764430505887, 18.844443423359284, 94.9558563009197, 56.74727530891187, 54.38466495621606, 156.03950602306492, 18.08013873755965, 9.403597119548579, 20.86511772100988, 48.612463731101826, 40.68873780723625, 13.177970389348966, 19.710136076767288, 33.02868895676784, 36.37187648588037, 5.197196798069252, 134.23758526834223, 5.194657368466044, 50.359583428338496, 6.130650320560632, 20.621866376959694, 7.827290815186451, 56.58968508150255, 27.838495016141508, 32.64924536459766, 57.10647258336654, 95.47082935762646, 65.0469683102413, 87.17718458077739, 143.45069910904334, 55.25145213608673, 8.075166702332739, 13.984158600434942, 28.163170454450643, 39.742279846785465, 16.447943599481203, 26.344085638217773, 7.9014542191118045, 33.6196909709104, 8.378331185072366, 12.025680777069136, 72.14012793695497, 26.296930324574895, 27.435283377491807, 24.655539542425203, 38.21687445841762, 49.04456986496268, 12.107499539665593, 10.415604643372966, 83.50004305336114, 22.172590894844113, 88.5271409191658, 9.259784857972262, 183.70537593482115, 104.36363148185106, 7.696088935860922, 39.52729238475217, 12.140113259435068, 16.91736667687468, 97.6734639670456, 9.875356756250646, 17.175779146125794, 6.035715432654942, 131.42689031845754, 13.29891426661904, 11.490635711230624, 26.231679848766746, 13.480727851579072, 72.92436678297061, 29.83397458290599, 93.3119589020483, 214.57346518298095, 38.792456305738035, 64.08595584596128, 6.231261770904037, 26.81023192452331, 52.37849212813234, 338.7829110772877, 80.64909782121588, 13.691368987871545, 58.784159352082575, 41.33406431604816, 60.22604450369563, 123.10683518221651, 68.60351925005239, 104.52289419445616, 20.204207069586374, 8.31501541416413, 37.85999763025312, 41.22648052905725, 38.77118722406425, 7.010181070636918, 43.968786026681556, 24.041916461648107, 30.730799642842165, 18.54029988339476, 12.008210614430977, 44.94171361638127, 30.52742582369338, 17.616801392060374, 59.86572233242134, 34.440894620651065, 22.871166407260183, 15.655471628935373, 17.372029099047424, 46.52590805969684, 73.01809164742689, ...])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)