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 = 44090
shot = shot_no
identifier='H03-W0051_shot_'+str(shot)+'_450V'
detector = 'H03-W0051'
ds = np.DataSource('/tmp') # temporary storage for downloaded files
scalars_URL = 'http://golem.fjfi.cvut.cz/shots/{shot_no}/Diagnostics/PlasmaDetection/Results/{name}'
def get_scalar(shot_no, name):
return float(ds.open(scalars_URL.format(shot_no=shot_no, name=name)).read())
t_plasma_start = get_scalar(shot_no, 't_plasma_start')
t_plasma_end = get_scalar(shot_no, 't_plasma_end')
is_plasma = get_scalar(shot_no, 'b_plasma')
def get_file(shot, identifier):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/H03/{identifier}.t3pa'
url = URL.format(shot=shot, identifier=identifier)
try:
file_name_t3pa=url
with urllib.request.urlopen(file_name_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
ft3pa.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_name_t3pa
def get_file_calib(name_calib):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/calib_matrix_H03/{name_calib}.txt'
url = URL.format(shot=shot, name_calib=name_calib)
#print(url)
try:
file_calib=url
with urllib.request.urlopen(file_calib) as calib:
line = calib.readline()
line = line.decode('utf‐8')
calib.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_calib
def load_calib(file_calib):
with urllib.request.urlopen(file_calib) as fc:
calib=[] #vytvoreni 1D pole
for i in range(0,256): #tj. rozsah 0-255
temp = [] # docasne pole
for j in range(0,256):
temp.append(0) #naplneni docasneho pole 0
calib.append(temp) #naplneni pole a[] docasnym polem temp
for i in range(0,256): #nacteni calib matice do pole calib
line = fc.readline()
line = line.decode('utf‐8')
word=line.strip().split(' ')
for j in range(0,256):
calib[i][j]=float(word[j]) #i = radek, j = sloupec0
fc.close
return calib
def load_t3pa_file(file_t3pa):
index=[]
matrix_index=[]
ToA=[]
ToT=[]
FToA=[]
overflow=[]
pocet_udalosti = 0
with urllib.request.urlopen(file_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
while True:
line = ft3pa.readline()
line = line.decode('utf‐8')
word=line.strip().split('\t') #v t3pa souboru je oddelovac \t
if line == '':
break
index.append(word[0])
matrix_index.append(word[1])
ToA.append(float(word[2]))
ToT.append(float(word[3]))
FToA.append(float(word[4]))
overflow.append(float(word[5]))
pocet_udalosti = pocet_udalosti + 1
ft3pa.close
return index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti
def noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti): #tuto fci nemus9m explicitn2 volat - volam ji v fci load_t3pa
pocet=int(0) #pocet sumicich pixelu
konst=int(len(index)/1000)+1
noise_matrix_index=[]
for i in range(0,konst):
pom = [] # pomocne pole
k=0 #pomocna promenna - udava, kolik je v czklu ve skutecnosti udalosti - aby nebyla chyba 'list index out of range'
for j in range(0,1001):
if i*1000+j>=len(index):
break
pom.append(matrix_index[i*1000+j])
k=k+1
for m in range(0,k):
count=int(0) #pocet vvyskytu stejneho matrix index behem 1000 udalosti
index_=int(-1) #budu testovat, jestli pixel na ktery koukam je sumici (abych ho nezapocital 2x)
for p in range(0,pocet):
#index=int(p)
if pom[m]==noise_matrix_index[p]:
index_=p #pixel na ktery jsem uz koukal a byl sumici
break
if index_ >=0 and pom[m]==noise_matrix_index[index_]:
continue
for l in range(0,k):
if pom[m]==pom[l]:
count=count+1
####podminka na sumici pixely
if count>=50: #kdyz se pixel vyskytne behem tisice udalosti vicekrat nez toto cislo, je sumici
noise_matrix_index.append(pom[m])
#noise_matrix_index[pocet]=pom[i]
pocet=pocet+1
pom.clear()
pocet_udalosti=len(index)
for n in range (0,pocet_udalosti):
for o in range(0,len(noise_matrix_index)):
if n >=pocet_udalosti:
break
if(matrix_index[n]==noise_matrix_index[o]):
del matrix_index[n]
del index[n]
del ToA[n]
del ToT[n]
del FToA[n]
del overflow[n]
pocet_udalosti=pocet_udalosti-1
continue
return pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow
def t3pa_data(pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow):
#rovnou vyhodim sumici pixely
pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow=noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti)
RowNo=[]
ClmNo=[]
for i in range(0,len(matrix_index)):
RowNo.append(int(int(matrix_index[i]))//int(256))
ClmNo.append(int(int(matrix_index[i]))%int(256))
return index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti, RowNo, ClmNo
def hit_map(detector,hit_map_fig,RowNo,ClmNo):
plt.hist2d(RowNo,ClmNo,bins=(256,256),cmap='Blues')
cb=plt.colorbar()
cb.set_label('Counts in pixel')
plt.xlabel('x [pixel]')
plt.ylabel('y [pixel]')
plt.title(detector)
plt.savefig(hit_map_fig, dpi = 1000)
return
def energy(a, b, c, t, ToT, pocet_udalosti, RowNo, ClmNo):
E=[] #energy in keV
#for i in range (0,pocet_udalosti):
pom=0
for i in range (0,len(ToT)):
sqrt=float(0.0)
e1=float(0.0)
e2=float(0.0)
# promenna sqrt je vnitrek odmocniny
sqrt = (((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i])))*(((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i])))) + (float(4)*float(a[RowNo[i]][ClmNo[i]])*float(c[RowNo[i]][ClmNo[i]]))) #zmena oproti verzi VI
if float(sqrt)<float(0):
E.append(float(0))
else:
'''
V kalibracni matici a se obcas vyskytne 0 -> ve vypoctu energie
je tim padem deleni nulou -> energie diverguje. Jak to vyresit?
zatim polozim energii = 0 (kdyz a=0), pak se uvidi
nakonec udelam limitu vyrazu energie pro a->0 (L'hopital)
'''
if a[RowNo[i]][ClmNo[i]]==0:
e1=((float(t[RowNo[i]][ClmNo[i]]))/float(2)) + ((((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i]))*(float(t[RowNo[i]][ClmNo[i]]))) - 2*(float(c[RowNo[i]][ClmNo[i]])))/(float(2)*np.sqrt(float(sqrt))))
e2=((float(t[RowNo[i]][ClmNo[i]]))/float(2)) - ((((float(b[RowNo[i]][ClmNo[i]])+float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]])-float(ToT[i]))*(float(t[RowNo[i]][ClmNo[i]]))) - 2*(float(c[RowNo[i]][ClmNo[i]])))/(float(2)*np.sqrt(float(sqrt))))
else:
e1=((-(float(b[RowNo[i]][ClmNo[i]]) - (float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]]))-float(ToT[i])))+np.sqrt(float(sqrt)))/(float(2)*float(a[RowNo[i]][ClmNo[i]]))
e2=((-(float(b[RowNo[i]][ClmNo[i]]) - (float(a[RowNo[i]][ClmNo[i]])*float(t[RowNo[i]][ClmNo[i]]))-float(ToT[i])))-np.sqrt(float(sqrt)))/(float(2)*float(a[RowNo[i]][ClmNo[i]]))
if a[RowNo[i]][ClmNo[i]]<0:
e1=-1
e2=-1
if math.isnan(e1):
e1=-1
if math.isnan(e2):
e2=-1
if e1<0 and e2<0:
E.append(float(0))
if e1>=0 and e1>e2:
E.append(float(e1))
if e2>=0 and e2>e1:
E.append(float(e2))
if e1>=0 and e2==e1:
E.append(float(e1))
return E
def Time(ToA, FToA, pocet_udalosti, RowNo, ClmNo):
T=[] #time in ns
for i in range (0,pocet_udalosti):
Time=float(0.0)
Time=(float(ToA[i])-((float(FToA[i])/float(16))))*float(25)
T.append(float(Time))
return T
def remove_interactions_with_zero_energy(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T):
i=0
treshold=5.015347
while i < len(T):
if E[i]<treshold: #E[i] < energy treshold
index.pop(i)
matrix_index.pop(i)
ToA.pop(i)
ToT.pop(i)
FToA.pop(i)
overflow.pop(i)
RowNo.pop(i)
ClmNo.pop(i)
E.pop(i)
T.pop(i)
continue
i=i+1
return index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T
def clustering_new(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T):
dT=float(50)
indexCl, TCl,ECl, matrix_indexCl, ToACl,ToTCl,FToACl,RowNoCl,ClmNoCl,overflowCl=[],[],[],[],[],[],[],[],[],[]
StartLastElem=len(T)-1
indexCl.append(int(index[StartLastElem]))
TCl.append(float(T[StartLastElem]))
ECl.append(float(E[StartLastElem]))
matrix_indexCl.append(int(matrix_index[StartLastElem]))
RowNoCl.append(int(RowNo[StartLastElem]))
ClmNoCl.append(int(ClmNo[StartLastElem]))
ToACl.append(float(ToA[StartLastElem]))
ToTCl.append(float(ToT[StartLastElem]))
FToACl.append(float(FToA[StartLastElem]))
overflowCl.append(float(overflow[StartLastElem]))
del index[StartLastElem]
del T[StartLastElem]
del E[StartLastElem]
del matrix_index[StartLastElem]
del RowNo[StartLastElem]
del ClmNo[StartLastElem]
del ToA[StartLastElem]
del ToT[StartLastElem]
del FToA[StartLastElem]
del overflow[StartLastElem]
j=1
pom=float(TCl[0]+dT)
while(j >0):
if(len(T) == 0):
break
k=0
j=0
while (k<=(len(TCl)-1)):
i=len(T)-1
if(len(T) == 0):
break
pocet_sousedu=0 #pocet sousednich pixelu - mohou byt maximalne 4
delka=0
# verze X
count=0 #pomocna promanna, kterou urcuji, ze se ma nasledujici cyklus while projit jeste jednou, pokud je i = -1
while(float(T[i])<=(pom)):
delka=delka+1
if(((((int(RowNoCl[k]))==(int(RowNo[i])+1))or((int(RowNoCl[k]))==(int(RowNo[i])-1))) and ((int(ClmNoCl[k]))==(int(ClmNo[i])))) or (((int(RowNoCl[k]))==(int(RowNo[i]))) and (((int(ClmNoCl[k]))==(int(ClmNo[i])+1))or((int(ClmNoCl[k]))==(int(ClmNo[i])-1))))):
#beru jen pixely, které mají společnou jednu stranu.
#pixely, kter0 spolu sousedí přes roh neuvažuji
indexCl.append(int(index[i]))
TCl.append(float(T[i]))
ECl.append(float(E[i]))
matrix_indexCl.append(int(matrix_index[i]))
RowNoCl.append(int(RowNo[i]))
ClmNoCl.append(int(ClmNo[i]))
ToACl.append(float(ToA[i]))
ToTCl.append(float(ToT[i]))
FToACl.append(float(FToA[i]))
overflowCl.append(float(overflow[i]))
# Removes i-th Row
del index[i]
del T[i]
del E[i]
del matrix_index[i]
del RowNo[i]
del ClmNo[i]
del ToA[i]
del ToT[i]
del FToA[i]
del overflow[i]
j=j+1
i=len(T)-1
pocet_sousedu=pocet_sousedu+1
if(len(T) == 0):
break
if(pocet_sousedu==4):
break
continue
i=i-1
if(i==-1): # verze X
count=count+1
if(i<0 and len(T)>0): # verze X
i=0
if(count>1):
break
if(i>=len(T)):
break
k=k+1
if(len(TCl)>2):
indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl = insertionSort(indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl)
return T, indexCl,TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl
def insertionSort(indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl):
# Function to do insertion sort
# Traverse through 1 to len(arr)
for i in range(1, len(TCl)):
key = TCl[i]
# Move elements of arr[0..i-1], that are
# greater than key, to one position ahead
# of their current position
#ostatni
key1 = indexCl[i]
key2 = ECl[i]
key3 = matrix_indexCl[i]
key4 = RowNoCl[i]
key5 = ClmNoCl[i]
key6 = ToACl[i]
key7 = ToTCl[i]
key8 = FToACl[i]
key9 = overflowCl[i]
j = i-1
while j >= 0 and key < TCl[j] :
TCl[j + 1] = TCl[j]
#ostatni
indexCl[j + 1] = indexCl[j]
ECl[j + 1] = ECl[j]
matrix_indexCl[j + 1] = matrix_indexCl[j]
RowNoCl[j + 1] = RowNoCl[j]
ClmNoCl[j + 1] = ClmNoCl[j]
ToACl[j + 1] = ToACl[j]
ToTCl[j + 1] = ToTCl[j]
FToACl[j + 1] = FToACl[j]
overflowCl[j + 1] = overflowCl[j]
j -= 1
TCl[j + 1] = key
#ostatni
indexCl[j + 1] = key1
ECl[j + 1] = key2
matrix_indexCl[j + 1] = key3
RowNoCl[j + 1] =key4
ClmNoCl[j + 1] = key5
ToACl[j + 1] = key6
ToTCl[j + 1] = key7
FToACl[j + 1] = key8
overflowCl [j + 1] = key9
return indexCl, TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl
def file_t3pa_cls_new(file_t3pa_cls,T):
with open(file_t3pa_cls, "w", encoding="utf-8") as t3pa_cls:
t3pa_cls.write('%\n')
t3pa_cls.write('% Index Matrix Index [ RowNo, ClmNo ] ToA FToA ( ToA_in_ns ) ToT ( ToT_in_keV ) Overflow\n')
t3pa_cls.write('\n')
i=1
T_first=[]
E_tot=[]
while(len(T) > 0):
T, indexCl,TCl, ECl, matrix_indexCl, RowNoCl, ClmNoCl, ToACl, ToTCl, FToACl, overflowCl = clustering_new(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T)
Tfirst=float(TCl[0])
Tlast=float(TCl[len(TCl)-1])
dT=Tlast-Tfirst
Etot=float(0)
for k in range(0,len(TCl)):
Etot=Etot+float(ECl[k])
T_first.append(float(Tfirst))
dT=Tlast-Tfirst
E_tot.append(float(Etot))
t3pa_cls.write('# '+str(i)+', Nunmasked = '+str(len(TCl))+', Nmasked = 0, Ntot = '+str(len(TCl))+'\n')
t3pa_cls.write('# Tfirst = '+str(Tfirst)+' ns, Tlast = '+str(Tlast)+' ns, dT = '+str(dT)+' ns, Etot = '+str(Etot)+' keV\n')
for j in range(0,len(TCl)):
t3pa_cls.write(str(indexCl[j])+' '+str(matrix_indexCl[j])+' [ '+str(RowNoCl[j])+', '+str(ClmNoCl[j])+' ] '+str(ToACl[j])+' '+str(FToACl[j])+' ( '+str(TCl[j])+' ns ) '+str(ToTCl[j])+' ( '+str(ECl[j])+' keV ) '+str(overflowCl[j])+'\n')
t3pa_cls.write('\n')
i=i+1
t3pa_cls.close
return T_first, E_tot
def energy_spectrum_in_time(Tfirst, Etot): #dela histogram - energie zaznamenana v case
pom = 0
dt=100 #(ns) time width of 1 bin
T_first=0 #cas, kdy prisel trigger a yacalo mereni
T_last=(max(Tfirst)) #posledni z Tfirst
Delta_T = T_last - T_first
poc = int(int(Delta_T) / float(dt)) + 1 #pocet casovych oken
T_int_first=[] #cas
E=[] #energie
for i in range(0,poc):
T_int_first.append((i*dt) + dt/2)
E.append(0)
#XII
for j in range(0,len(Tfirst)):
time_index=0
time_index=int(((Tfirst[j]-T_first)/dt))
if float(Tfirst[j]-T_first) >= (T_int_first[time_index] - dt / 2) and float(Tfirst[j]-T_first) < (T_int_first[time_index] + dt / 2):
E[time_index]=float(E[time_index])+float(Etot[j])
pom=pom+1
for l in range(0,len(T_int_first)):
T_int_first[l]=T_int_first[l]+T_first
caption, T_int_first = energy_in_time_hist(T_int_first, E, figure_E_in_time_hist, t_plasma_start, t_plasma_end, is_plasma, dt)
return dt, caption, T_int_first, E
def energy_in_time_hist(T_int_first, E,figure_E_in_time_hist, t_plasma_start, t_plasma_end, is_plasma, dt):
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
for k in range(0,len(T_int_first)):
T_int_first[k] = T_int_first[k] / 1000000
plt.plot(T_int_first, E)
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Time [ms]')
plt.ylabel('Energy [keV]')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
plt.axvline(t, color='k', linestyle='--')
plt.xlim([0, (t_plasma_start + t_plasma_end)])
else:
plt.xlim(0,)
plt.ylim(0,) #10 000 keV
plt.savefig(figure_E_in_time_hist, dpi = 1000)
caption = '# x = time in ms, count = energy in keV, dT= '+str(dt)+' ns'
return caption, T_int_first
def hits_in_time_hist_new(T, dt, t_plasma_start, t_plasma_end, is_plasma,figure_count_in_time_hist):
pom = 0
T_first=0 #cas, kdy prisel trigger a yacalo mereni
T_last=(max(T)) #posledni z Tfirst
Delta_T = T_last - T_first
poc = int(int(Delta_T) / float(dt)) + 1 #pocet casovych oken
T_hit=[] #cas
count=[] #energie
for i in range(0,poc):
T_hit.append((i*dt) + dt/2)
count.append(0)
for j in range(0,len(T)):
time_index=0
time_index=int(((T[j]-T_first)/dt))
k=time_index
for j in range(0,len(T)):
time_index=0
time_index=int(((T[j]-T_first)/dt))
if float(T[j]-T_first) >= (T_hit[time_index] - dt / 2) and float(T[j]-T_first) < (T_hit[time_index] + dt / 2):
count[time_index] = count[time_index] + 1
pom=pom+1
for l in range(0,len(T_hit)):
T_hit[l]=T_hit[l]+T_first
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
for k in range(0,len(T_hit)):
T_hit[k] = T_hit[k] / 1000000
plt.plot(T_hit, count)
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Time [ms]')
plt.ylabel('Count')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
plt.axvline(t, color='k', linestyle='--')
plt.xlim([0, (t_plasma_start + t_plasma_end)])
else:
plt.xlim(0,)
plt.ylim(0,) #10 000 keV
plt.savefig(figure_count_in_time_hist, dpi = 1000)
caption = '# x = time in ms, dT= '+str(dt)+' ns'
return caption, T_hit,count
def energy_spectrum(Etot):
E_min=0
dE=5 #keV
E_max=max(Etot)
pocet=(E_max//dE) + 3
pocet=int(pocet)
E_max=float(dE*pocet)
xle=[]
xre=[]
xmean=[]
for p in range (0,pocet):
xle.append(E_min + (p * (E_max - E_min)) / pocet)
xre.append(xle[p]+dE)
xmean.append((xle[p] + xre[p]) / 2)
count=[]
for l in range(0,pocet):
count.append(0)
#XII
for i in range(0,len(Etot)):
E_index=int(((Etot[i]-E_min)/dE))
if ((xle[E_index] <= Etot[i]) and (Etot[i] < xre[E_index])):
count[E_index]=count[E_index]+1
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
ax.hist(Etot, bins = xle)
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Count')
plt.xlim(0,)
ax.set_yscale('log') #log scale y
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
plt.savefig(figure_E_hist, dpi = 1000)
return caption, xmean,count, xle, Etot
def hist_file(file_hist, xmean, count, caption ):
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('#\n')
hist.write('#'+str(caption)+'\n')
hist.write('# x_mean count\n')
hist.write('\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(count[m])+'\n')
hist.close
return T_first, E_tot
def multiplot(icon_fig, x1,y1,x2,y2):
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(nrows=2,figsize =(10, 7))
ax[0].plot(x1, y1)
ax[0].set_xlabel('Time [ms]')
ax[0].set_ylabel('Energy [keV]')
if is_plasma == 1:
for t in (t_plasma_start, t_plasma_end):
ax[0].axvline(t, color='k', linestyle='--')
ax[0].set_xlim([0, (t_plasma_start + t_plasma_end)])
else:
ax[0].set_xlim(0,)
ax[0].set_ylim(0,) #keV
ax[1].hist(y2, bins = x2)
ax[1].set_xlabel('Energy [keV]')
ax[1].set_ylabel('Count')
ax[1].set_xlim(0,)
#ax[1].set_ylim(0,)
ax[1].set_yscale('log') #log scale y
fig.subplots_adjust(hspace=0.3)
plt.savefig(icon_fig, dpi = 1000)
return
#soubory, ktere ctu
#read files
t3pa=get_file(shot, identifier)
name_calib='caliba'
caliba=get_file_calib(name_calib)
name_calib='calibb'
calibb=get_file_calib(name_calib)
name_calib='calibc'
calibc=get_file_calib(name_calib)
name_calib='calibt'
calibt=get_file_calib(name_calib)
#vytvorene soubory:
#created files
t3pa_cls='H03-W0051_shot_'+str(shot)+'_450V.t3pa_cls'
E_hist='H03-W0051_shot_'+str(shot)+'_450V_E_hist.txt'
E_in_time_hist='H03-W0051_shot_'+str(shot)+'_450V_discharge_energy.txt'
count_in_time_hist= 'H03-W0051_shot_'+str(shot)+'_450V_discharge_hits.txt'
#created figures
icon_fig='icon-fig'
figure_E_in_time_hist='discharge_energy'
figure_count_in_time_hist='discharge_hits'
figure_E_hist='Energy_spectrum'
hit_map_fig='hit-map'
#nactu jednotlive kalibracni matice - abych to nemusel delat v kazde funkci
a=load_calib(caliba)
b=load_calib(calibb)
c=load_calib(calibc)
t=load_calib(calibt)
#nactu a urcim jednotlive hodnoty - abych to nemusel delat v kazde funkci
index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti = load_t3pa_file(t3pa)
index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti, RowNo, ClmNo = t3pa_data(pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow)
raw data
#hit map
hit_map(detector,hit_map_fig,RowNo,ClmNo)
Energy and time calculation from raw data.
E=energy(a, b, c, t, ToT, pocet_udalosti, RowNo, ClmNo)
T=Time(ToA, FToA, pocet_udalosti, RowNo, ClmNo)
index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T = remove_interactions_with_zero_energy(index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E, T)
#sort by time
T, index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E = (list(t) for t in zip(*sorted(zip(T, index, matrix_index, ToA, ToT, FToA, overflow, RowNo, ClmNo, E), reverse=True))) #serazeni od nejvetsiho po nejmensi
T_pom=T.copy()
#save to file
T_first, E_tot = file_t3pa_cls_new(t3pa_cls,T)
dt, caption, T_int_first, E = energy_spectrum_in_time(T_first, E_tot)
hist_file(E_in_time_hist, T_int_first, E, caption);
([7061593.75, 7094107.8125, 7136457.8125, 7165618.75, 7314129.6875, 7369234.375, 7379143.75, 7475423.4375, 7507050.0, 7516884.375, 7516917.1875, 7532898.4375, 7545882.8125, 7665126.5625, 7731723.4375, 7740578.125, 7791290.625, 7796645.3125, 7894650.0, 8063903.125, 8064143.75, 8081240.625, 8149629.6875, 8173607.8125, 8196509.375, 8296281.25, 8298535.9375, 8332868.75, 8335373.4375, 8345445.3125, 8357275.0, 8381796.875, 8384681.25, 8388559.375, 8397385.9375, 8437132.8125, 8444285.9375, 8449117.1875, 8457600.0, 8507312.5, 8540082.8125, 8546862.5, 8572471.875, 8595956.25, 8611620.3125, 8621629.6875, 8648803.125, 8692535.9375, 8801998.4375, 8869753.125, 8869776.5625, 8870157.8125, 8892581.25, 8909712.5, 8912267.1875, 8924081.25, 8937540.625, 8944387.5, 9077096.875, 9079775.0, 9080823.4375, 9105562.5, 9114835.9375, 9126712.5, 9145242.1875, 9164987.5, 9168923.4375, 9181139.0625, 9184987.5, 9190543.75, 9212739.0625, 9226307.8125, 9227689.0625, 9228715.625, 9235398.4375, 9242964.0625, 9283775.0, 9307321.875, 9308928.125, 9308971.875, 9322628.125, 9330110.9375, 9372259.375, 9387946.875, 9390525.0, 9396026.5625, 9420315.625, 9428125.0, 9431637.5, 9432137.5, 9432303.125, 9432321.875, 9435570.3125, 9468221.875, 9471503.125, 9471971.875, 9528312.5, 9538973.4375, 9602695.3125, 9636742.1875, 9663726.5625, 9717195.3125, 9762660.9375, 9762712.5, 9770528.125, 9815509.375, 9819387.5, 9837284.375, 9862785.9375, 9867406.25, 9940407.8125, 9987234.375, 13897032.8125, 13903092.1875, 13928628.125, 13936584.375, 13949195.3125, 13995673.4375, 14015260.9375, 14019679.6875, 14022262.5, 14023420.3125, 14023457.8125, 14024776.5625, 14026320.3125, 14027896.875, 14027903.125, 14027970.3125, 14027996.875, 14030532.8125, 14031292.1875, 14033853.125, 14037579.6875, 14038510.9375, 14040364.0625, 14043387.5, 14048459.375, 14050093.75, 14050621.875, 14050696.875, 14050860.9375, 14052976.5625, 14053023.4375, 14054107.8125, 14054814.0625, 14055162.5, 14057601.5625, 14062215.625, 14062654.6875, 14063395.3125, 14073559.375, 14077432.8125, 14079910.9375, 14082051.5625, 14083595.3125, 14087943.75, 14090359.375, 14091442.1875, 14091701.5625, 14092000.0, 14092485.9375, 14092987.5, 14093596.875, 14096893.75, 14098607.8125, 14099900.0, 14104860.9375, 14104918.75, 14106817.1875, 14107292.1875, 14107464.0625, 14107500.0, 14107642.1875, 14108034.375, 14108104.6875, 14108226.5625, 14108246.875, 14108579.6875, 14108695.3125, 14108748.4375, 14108998.4375, 14109448.4375, 14109643.75, 14109675.0, 14109957.8125, 14110004.6875, 14110053.125, 14110567.1875, 14110571.875, 14110729.6875, 14111112.5, 14111128.125, 14111203.125, 14111695.3125, 14111703.125, 14112042.1875, 14112104.6875, 14112287.5, 14112400.0, 14112401.5625, 14112532.8125, 14113296.875, 14113303.125, 14113503.125, 14114776.5625, 14114965.625, 14115142.1875, 14116695.3125, 14117165.625, 14118179.6875, 14118445.3125, 14118448.4375, 14119423.4375, 14119635.9375, 14119814.0625, 14119906.25, 14120198.4375, 14120276.5625, 14120325.0, 14120854.6875, 14120954.6875, 14121303.125, 14121356.25, 14122217.1875, 14122653.125, 14122839.0625, 14122981.25, 14123003.125, 14123051.5625, 14123315.625, 14123360.9375, 14123365.625, 14123476.5625, 14124060.9375, 14124554.6875, 14124560.9375, 14124871.875, 14124940.625, 14125478.125, 14125637.5, 14125768.75, 14126159.375, 14126185.9375, 14126329.6875, 14126412.5, 14126518.75, 14126742.1875, 14126865.625, 14127070.3125, 14127175.0, 14127382.8125, 14127995.3125, 14128059.375, 14128065.625, 14128264.0625, 14128284.375, 14128451.5625, 14129212.5, 14129837.5, 14129945.3125, 14130034.375, 14130248.4375, 14130300.0, 14130306.25, 14130351.5625, 14130531.25, 14130751.5625, 14131023.4375, 14131317.1875, 14131418.75, 14132281.25, 14132332.8125, 14132343.75, 14132353.125, 14132528.125, 14132748.4375, 14133131.25, 14133604.6875, 14133862.5, 14134568.75, 14134798.4375, 14135043.75, 14135828.125, 14135989.0625, 14136137.5, 14136423.4375, 14136735.9375, 14136982.8125, 14137323.4375, 14137462.5, 14137632.8125, 14137753.125, 14137996.875, 14138070.3125, 14138200.0, 14138242.1875, 14138253.125, 14138320.3125, 14138354.6875, 14138454.6875, 14138465.625, 14138535.9375, 14138562.5, 14138801.5625, 14138868.75, 14139028.125, 14139065.625, 14139078.125, 14139262.5, 14139335.9375, 14139535.9375, 14139651.5625, 14139687.5, 14140417.1875, 14140503.125, 14140587.5, 14140775.0, 14140829.6875, 14140850.0, 14140882.8125, 14140992.1875, 14141443.75, 14142292.1875, 14142407.8125, 14142443.75, 14142807.8125, 14142812.5, 14143096.875, 14143157.8125, 14143176.5625, 14143240.625, 14143267.1875, 14143450.0, 14143456.25, 14143457.8125, 14143484.375, 14143887.5, 14144118.75, 14144193.75, 14144375.0, 14144512.5, 14144529.6875, 14144742.1875, 14144743.75, 14145101.5625, 14145293.75, 14145401.5625, 14145543.75, 14145615.625, 14146081.25, 14146143.75, 14146489.0625, 14146518.75, 14146684.375, 14146764.0625, 14146810.9375, 14146868.75, 14147015.625, 14147353.125, 14147364.0625, 14147365.625, 14147634.375, 14147801.5625, 14147831.25, 14148221.875, 14148665.625, 14148689.0625, 14148896.875, 14148960.9375, 14148968.75, 14148995.3125, 14149040.625, 14149048.4375, 14149171.875, 14149173.4375, 14149204.6875, 14149210.9375, 14149239.0625, 14149292.1875, 14149407.8125, 14149548.4375, 14149626.5625, 14149679.6875, 14149842.1875, 14149842.1875, 14149873.4375, 14149875.0, 14149917.1875, 14149923.4375, 14149931.25, 14150006.25, 14150023.4375, 14150057.8125, 14150059.375, 14150171.875, 14150210.9375, 14150231.25, 14150240.625, 14150279.6875, 14150282.8125, 14150287.5, 14150440.625, 14150471.875, 14150484.375, 14150520.3125, 14150589.0625, 14150615.625, 14150646.875, 14150664.0625, 14150676.5625, 14150723.4375, 14150837.5, 14150920.3125, 14150945.3125, 14151031.25, 14151081.25, 14151201.5625, 14151245.3125, 14151351.5625, 14151476.5625, 14151514.0625, 14151557.8125, 14151609.375, 14151610.9375, 14151623.4375, 14151626.5625, 14151907.8125, 14152131.25, 14152246.875, 14152264.0625, 14152531.25, 14152603.125, 14152712.5, 14152737.5, 14152754.6875, 14152889.0625, 14152982.8125, 14153065.625, 14153101.5625, 14153160.9375, 14153220.3125, 14153245.3125, 14153257.8125, 14153278.125, 14153296.875, 14153523.4375, 14153556.25, 14153709.375, 14153781.25, 14153887.5, 14153992.1875, 14154017.1875, 14154028.125, 14154029.6875, 14154053.125, 14154098.4375, 14154131.25, 14154167.1875, 14154346.875, 14154371.875, 14154412.5, 14154442.1875, 14154468.75, 14154523.4375, 14154728.125, 14154745.3125, 14154840.625, 14155045.3125, 14155142.1875, 14155357.8125, 14155379.6875, 14155414.0625, 14155490.625, 14155507.8125, 14155603.125, 14155615.625, 14155620.3125, 14155671.875, 14155682.8125, 14155704.6875, 14155706.25, 14155707.8125, 14155725.0, 14155745.3125, 14155765.625, 14155767.1875, 14155778.125, 14155781.25, 14155853.125, 14155989.0625, 14156073.4375, 14156081.25, 14156114.0625, 14156339.0625, 14156535.9375, 14156575.0, 14156714.0625, 14156725.0, 14156839.0625, 14156851.5625, 14156985.9375, 14157046.875, 14157051.5625, 14157053.125, 14157354.6875, 14157404.6875, 14157528.125, 14157545.3125, 14157750.0, 14157906.25, 14157981.25, 14158114.0625, 14158123.4375, 14158150.0, 14158296.875, 14158418.75, 14158445.3125, 14158450.0, 14158459.375, 14158490.625, 14158514.0625, 14158632.8125, 14158714.0625, 14158742.1875, 14158765.625, 14158856.25, 14159012.5, 14159321.875, 14159334.375, 14159495.3125, 14159534.375, 14159546.875, 14159564.0625, 14159585.9375, 14159593.75, 14159715.625, 14159753.125, 14159845.3125, 14159854.6875, 14159971.875, 14159992.1875, 14160057.8125, 14160104.6875, 14160114.0625, 14160181.25, 14160351.5625, 14160354.6875, 14160406.25, 14160451.5625, 14160476.5625, 14160726.5625, 14160832.8125, 14160848.4375, 14160876.5625, 14160895.3125, 14160976.5625, 14161010.9375, 14161070.3125, 14161079.6875, 14161150.0, 14161150.0, 14161203.125, 14161239.0625, 14161260.9375, 14161278.125, 14161300.0, 14161309.375, 14161310.9375, 14161403.125, 14161448.4375, 14161462.5, 14161482.8125, 14161489.0625, 14161501.5625, 14161546.875, 14161551.5625, 14161557.8125, 14161568.75, 14161596.875, 14161606.25, 14161642.1875, 14161651.5625, 14161660.9375, 14161743.75, 14161832.8125, 14161857.8125, 14162048.4375, 14162209.375, 14162246.875, 14162329.6875, 14162453.125, 14162492.1875, 14162517.1875, 14162554.6875, 14162623.4375, 14162646.875, 14162690.625, 14162765.625, 14162942.1875, 14162993.75, 14162996.875, 14163006.25, 14163035.9375, 14163151.5625, 14163176.5625, 14163353.125, 14163376.5625, 14163400.0, 14163417.1875, 14163512.5, 14163551.5625, 14163609.375, 14163787.5, 14163898.4375, 14163903.125, 14163906.25, 14163979.6875, 14163981.25, 14163992.1875, 14163996.875, 14164181.25, 14164562.5, 14164750.0, 14164784.375, 14164821.875, 14164829.6875, 14164962.5, 14165006.25, 14165042.1875, 14165078.125, 14165243.75, 14165342.1875, 14165428.125, 14165503.125, 14165590.625, 14165609.375, 14165635.9375, 14165660.9375, 14165681.25, 14165689.0625, 14165709.375, 14165740.625, 14165760.9375, 14165765.625, 14165776.5625, 14165785.9375, 14165804.6875, 14165818.75, 14165859.375, 14165898.4375, 14165907.8125, 14165910.9375, 14165929.6875, 14165976.5625, 14165978.125, 14166018.75, 14166132.8125, 14166157.8125, 14166246.875, 14166440.625, 14166603.125, 14166671.875, 14166792.1875, 14166901.5625, 14166943.75, 14167029.6875, 14167082.8125, 14167084.375, 14167142.1875, 14167143.75, 14167160.9375, 14167364.0625, 14167459.375, 14167664.0625, 14167792.1875, 14167800.0, 14167812.5, 14168106.25, 14168115.625, 14168290.625, 14168300.0, 14168579.6875, 14168857.8125, 14169046.875, 14169220.3125, 14169687.5, 14169921.875, 14169996.875, 14170292.1875, 14170348.4375, 14170384.375, 14171021.875, 14171310.9375, 14171335.9375, 14171637.5, 14171764.0625, 14171826.5625, 14171915.625, 14172003.125, 14172129.6875, 14172282.8125, 14172326.5625, 14172656.25, 14172676.5625, 14172878.125, 14173100.0, 14173310.9375, 14173646.875, 14173726.5625, 14173798.4375, 14174464.0625, 14174809.375, 14174962.5, 14174967.1875, 14174976.5625, 14175090.625, 14175114.0625, 14175153.125, 14175573.4375, 14175637.5, 14175790.625, 14176089.0625, 14176718.75, 14176718.75, 14176867.1875, 14177201.5625, 14177900.0, 14178412.5, 14178857.8125, 14178859.375, 14178993.75, 14179395.3125, 14179567.1875, 14180114.0625, 14180900.0, 14181170.3125, 14182025.0, 14182746.875, 14183140.625, 14184557.8125, 14184934.375, 14184971.875, 14184990.625, 14185576.5625, 14185637.5, 14185712.5, 14185879.6875, 14185942.1875, 14186414.0625, 14186792.1875, 14187131.25, 14187528.125, 14188085.9375, 14188307.8125, 14188537.5, 14188575.0, 14188901.5625, 14189168.75, 14189495.3125, 14189845.3125, 14190823.4375, 14191256.25, 14191310.9375, 14191317.1875, 14192559.375, 14192592.1875, 14193084.375, 14193559.375, 14194429.6875, 14196209.375, 14197817.1875, 14199010.9375, 14199012.5, 14199107.8125, 14199748.4375, 14199756.25, 14200703.125, 14201914.0625, 14202307.8125, 14202560.9375, 14202989.0625, 14203189.0625, 14203348.4375, 14203595.3125, 14203756.25, 14204378.125, 14205239.0625, 14205481.25, 14205678.125, 14205978.125, 14206098.4375, 14206164.0625, 14206626.5625, 14207867.1875, 14210195.3125, 14210546.875, 14211262.5, 14212501.5625, 14214840.625, 14215814.0625, 14215979.6875, 14215981.25, 14216007.8125, 14216457.8125, 14216664.0625, 14221175.0, 14221579.6875, 14223371.875, 14224518.75, 14224739.0625, 14224815.625, 14225807.8125, 14225965.625, 14226720.3125, 14230601.5625, 14232651.5625, 14233571.875, 14234320.3125, 14234446.875, 14240440.625, 14240776.5625, 14243726.5625, 14245967.1875, 14248412.5, 14249862.5, 14250431.25, 14256106.25, 14262431.25, 14270289.0625, 14271584.375, 14282850.0, 14283479.6875, 14288748.4375, 14288756.25, 14316635.9375, 442749226.5625], [14.353494663391347, 473.72552766437803, 64.83849779080148, 36.77691193062272, 45.5215314765082, 9.35358762473517, 138.96490170140314, 10.736181369094693, 63.82109119366702, 55.938173334457915, 88.7681590638431, 19.734325713764772, 35.651572755161176, 46.550626326252385, 20.871126056680342, 21.84359957165835, 47.23218286247899, 6.13684904713674, 9.24638261305465, 30.44021871282233, 47.03775059633008, 6.55245419575673, 58.702711176285426, 23.636303341479437, 18.179892111916878, 89.94716101189168, 5.1297197217189465, 12.358679632605428, 10.71630147318415, 93.0766903499036, 17.52061275270044, 24.991509180526812, 20.46416628164172, 20.60424365329352, 9.121317134698764, 33.91656963226593, 19.029721329047504, 51.49528148779245, 88.70332418188231, 82.62434546059671, 5.2298468191344565, 58.33207553908103, 41.93331390532528, 86.784007584458, 51.30684960822346, 8.61703840863088, 32.63254148345757, 43.17868874281369, 12.753193533769785, 34.34814232164087, 10.167794305449949, 87.83604021181014, 39.235160788821744, 18.255997217943104, 6.01336053078085, 11.0482600586928, 32.090075499029574, 68.50731007923501, 39.51983998604173, 29.065530742111154, 25.36276690784278, 5.921981987503392, 84.37475442262945, 11.488771362701042, 50.02737375445998, 22.722740087727537, 17.169690837376297, 91.45494429636368, 21.335876551298693, 6.5764150127191465, 111.03219718051129, 85.299902567554, 23.945056982174265, 38.52536228253169, 28.972469800237363, 66.00821662089689, 13.676826791699021, 72.9907585653588, 95.70378321655397, 62.12037854049055, 66.7891484470813, 54.47259835167837, 27.28921849345727, 39.35013358710248, 28.2602832444484, 65.7435275050917, 12.801872069724844, 25.61020355370054, 16.63461569573696, 46.28079114633559, 96.29792409196753, 89.07687421935228, 44.22610227911099, 14.941244148153501, 13.057208369360833, 12.69087761508895, 38.26527543577155, 107.25872188508636, 35.72856526079983, 58.60820293278006, 7.372982563358372, 26.732792541294373, 68.17684566975235, 5.354712784245878, 64.41885900111, 66.85545750295871, 28.996827169854154, 112.19906756178189, 88.49173997056705, 28.492001388069177, 14.434911830894555, 14.534011182788714, 80.34683665108685, 29.485995613491788, 10.977925767294726, 105.69469775692211, 17.168300381557057, 20.036805185614934, 6.328084087484737, 39.78822185923844, 22.85319254143005, 44.848943485916244, 51.88871040089554, 20.871193667039336, 26.467093006754947, 15.410974375179872, 89.7989222877052, 5.1492059434781, 5.387654290123565, 16.536214208219764, 26.62462523961121, 33.20454150504682, 10.211351191310994, 19.601035348003556, 27.953888009990216, 140.85742704841465, 87.5724959221242, 7.7884728836066905, 17.888564589548565, 9.29390012359141, 27.901629567809184, 61.82385221825618, 69.82649026251953, 8.898215595327333, 95.92870616255469, 35.27538308689147, 18.52547826745096, 19.67910539562954, 77.94228142861903, 5.71437746038673, 44.47637852710898, 56.36949227037983, 29.022241895616872, 19.432695983014664, 26.510717661558907, 25.756896138374913, 11.326472613546093, 65.36063868588136, 50.949582862337884, 19.770498431721553, 13.538392312171128, 173.75082529786062, 42.85044921383103, 14.150076505973443, 113.0498348797251, 28.173464299292945, 7.611891421922517, 5.09488032341785, 35.222522822861414, 19.14699249790841, 58.64371926551527, 60.664350018351925, 23.22640493206027, 6.795742014943427, 68.45664886767102, 6.519365833077661, 81.72007403771516, 26.181263879217074, 17.219910229944606, 32.86166634729164, 44.27113433493194, 8.566632772941917, 7.517428094158713, 60.515464745670485, 72.78434080480221, 28.199652866668668, 8.960780437403438, 21.635228861273784, 14.910618693956629, 29.148818182271953, 5.216270756522859, 81.2912865049599, 46.614024831305414, 67.18256762792173, 15.867553571281848, 10.270541321095285, 73.82180026059304, 20.058651924123808, 37.071630126810064, 58.4809364584972, 22.550319185008476, 6.933272729737776, 25.587246155338363, 73.06943091049588, 97.59715330110649, 5.712836102823359, 83.06397033333556, 47.11838201349126, 25.303326137019152, 40.30559279508621, 11.340039130667057, 29.153939809226653, 64.22302227576327, 36.7537455405884, 17.5776176592042, 8.118588397117676, 29.570597285957536, 25.396005680513376, 18.822971165666523, 36.453855472548426, 18.710383379088043, 57.55298437809341, 29.590219543437573, 48.267791310782805, 46.37818819078264, 25.73942578247345, 29.593005254134496, 58.61179312681341, 5.744956007708263, 78.00161546103493, 11.199893235750558, 22.72205252054634, 97.31172513930277, 32.675650253293114, 80.89593158843272, 50.725112666314686, 21.639503016235178, 5.234156784139545, 21.2604218195829, 31.40272912710506, 38.89152686193783, 23.576606320267587, 24.12015277532069, 80.56142884208839, 61.783822976485546, 130.96431079441888, 7.225273460899372, 25.27921487972804, 9.823766914230644, 12.560017545618122, 58.076908098202786, 41.12700090295604, 76.07454717282067, 150.85943632080748, 27.65022827389516, 5.57519222573279, 11.76680217485271, 5.640250335941951, 38.65796054701666, 47.792779484209134, 5.357885297869613, 26.029438017087415, 29.604191167731127, 42.795952554645254, 68.58387361761261, 5.901564481523938, 11.45765452471494, 10.330741219667953, 115.23435989695996, 25.440892274317907, 43.92058086375303, 14.68165371458611, 77.07678788115845, 72.06944552585097, 10.3784304246614, 20.516357260276195, 29.66975523850248, 37.390459818246526, 45.10135817742194, 7.822877740003752, 6.9047549546724545, 69.02149245526022, 115.91437765085564, 18.29596293877664, 91.04978680844967, 19.987977949081962, 13.58785447880324, 32.631274458799595, 122.34863721889519, 8.687979854997394, 67.43830196743917, 19.796190578387257, 9.70874289866778, 5.3208183280715105, 26.166139992579883, 52.352664071284345, 65.78788479758539, 32.96582134133165, 22.31088865752372, 53.04019828907692, 7.6305856592876875, 68.29311692352044, 5.198107194523418, 14.925240601590707, 70.30190856292792, 15.257973668040526, 24.523792994991883, 12.128948368450294, 61.880366493780315, 64.35540451634562, 56.12343846276438, 6.036973747837919, 27.986491478131114, 9.635269287039392, 15.384076907805271, 6.379008749669283, 33.151920247776644, 16.08108154619088, 45.126872121236985, 14.26252348097054, 21.745799117344784, 66.19139995038633, 10.092529100972401, 32.880980227759096, 59.18544360724054, 62.37415011395015, 23.700612120865156, 39.22478817111433, 12.128612534526674, 46.7952996734885, 17.240567391596084, 19.563062491699153, 32.77943073174624, 75.9244532855812, 6.2623924026738536, 27.785666873615423, 6.0800116185134, 24.544525550550887, 190.73086764807385, 68.7782447723781, 6.76145788477542, 88.4272289875377, 32.31481695488803, 28.392484217973628, 9.289049871627414, 35.26804737796784, 162.9244149185386, 20.013158130557656, 12.28036558966326, 181.68393035370133, 104.97754878116041, 10.516428442332703, 21.81040952229024, 16.018448455714275, 18.646876630704586, 11.558945118318242, 34.64340098544084, 43.12197024156084, 26.711302433437393, 57.47291746726988, 9.971635292758434, 13.454504760957123, 42.061852668240164, 65.88109421558121, 5.839771119107743, 58.50996494348714, 14.806594037032847, 9.5681734380959, 79.75198020536143, 34.78258862659328, 12.586082027440696, 27.273693806676206, 23.35680763294256, 11.284664052981473, 57.65542306537539, 100.14971778596654, 31.78167166919605, 45.64889588083866, 38.65420264274514, 79.66758193074546, 17.1201419980156, 102.9097270235377, 17.219748229071147, 5.368158579478964, 31.601807003765117, 71.76241143826363, 32.939004211533664, 42.155056073147755, 11.987017610004475, 9.748965266877068, 19.767290591809402, 7.110682854338585, 7.392816442907737, 24.204429172715763, 91.49244085784125, 8.293983059766695, 23.937622027379852, 14.20480637365656, 9.68367346001704, 85.73824785275004, 60.35682557482924, 133.2052298197641, 11.943339722888041, 72.69601261208908, 120.63466956028064, 91.80284339144525, 9.815860612181233, 6.559279010428072, 33.564597357816474, 9.555392088807762, 78.3831574207153, 27.179852850244465, 6.594196462714507, 40.75154635054125, 36.724422283613656, 19.247476621369373, 37.929071005299356, 17.703290882445362, 107.20111832155399, 120.53668625294208, 16.841635354393926, 6.1544185316873605, 135.75325572677988, 74.31783458236261, 47.199725363220104, 49.58537836937867, 57.278543193439255, 53.959106161268366, 77.60096407933997, 58.56103256615294, 144.74222144461513, 85.70036375992598, 176.41859715798293, 11.115712889705602, 54.354030798643635, 30.350359721749584, 6.293897324225599, 55.71879324120669, 29.430641341484563, 107.56090829397996, 16.845050039553357, 40.235935353644614, 22.69808557703132, 76.87226004600889, 15.579405560730098, 10.010674334203065, 19.290120177564052, 38.52634537582876, 6.997355117862078, 13.638581184715129, 48.785429134526986, 42.00743014501606, 101.61683499211131, 73.20579185488225, 27.091800585534948, 30.619262391062065, 42.29559070599094, 6.137766473903482, 10.385163457992718, 43.24332129557821, 8.060733604365055, 50.43517550739337, 18.16350169807365, 9.86721091156452, 11.406037867396963, 17.72729173742572, 103.40576469931088, 75.34640461928947, 5.7576106600359225, 30.240818891386503, 37.46407154749019, 11.442331310233499, 79.65797703991664, 51.40732769005136, 65.57899699966401, 20.5501338259211, 29.623393177108042, 25.338588503117435, 128.6027701241257, 58.325084278925424, 47.76298990475915, 17.595070019309876, 9.867098645014597, 132.67222661305547, 8.921745451046558, 50.526320092458434, 51.58389783260908, 32.022575343000476, 31.963147853396553, 71.49294496650889, 105.55434234269504, 112.65856204740324, 29.25502346049976, 171.86869265024822, 58.31970632752492, 49.932414663227775, 80.57132976094407, 65.61759374109833, 23.22642695072207, 166.34354258819423, 27.542601311789923, 25.772105636994027, 63.93691962545487, 66.90935078348784, 106.14916392544066, 10.939203343101395, 42.124368198549114, 10.844508684040633, 7.241840816275167, 6.73736242483435, 11.593779428282327, 18.335779657925514, 21.41947902885403, 18.753579045500953, 19.35315697538517, 107.13513831793823, 31.550913024453926, 102.6628717816799, 35.88619177362325, 17.235990196897486, 24.141560673491583, 22.806446988252, 62.00804244895558, 21.8257731318311, 8.383265822904729, 11.875105911471099, 6.636418077090405, 5.961896993654193, 12.40819180158402, 72.92229378544869, 16.950758129527294, 12.532837565835399, 21.40140000993174, 5.978316730841685, 163.52628580424195, 18.29264744716334, 9.53987319215665, 23.17433371053665, 27.759873069797564, 25.77991529317198, 35.19983605801254, 38.35084678260363, 10.69775915491015, 15.950459283085388, 16.93716050207172, 13.2905645636859, 37.018675301060775, 27.618674979992512, 57.673680193818264, 60.50167468273238, 23.631718542799774, 57.29710130082808, 26.36064094057911, 15.111439999683835, 32.86672071402299, 7.4691393743382815, 60.48189238985312, 79.35631343185365, 24.277194971053582, 34.47297681181391, 42.79701611969088, 53.08722263706735, 8.30078276625573, 42.66492192149025, 72.4143974575048, 23.492589951121623, 24.503793887336872, 24.5654317567233, 73.56338891035281, 78.6167044099696, 22.832434501843665, 45.346258978838904, 27.835536129408673, 5.908811296577237, 9.81245982135465, 26.736256819297026, 19.553482018962708, 32.033837241875744, 22.219346359091624, 17.374203878533606, 9.454173248858073, 16.905984015553795, 14.604699809172235, 5.892264206598048, 24.750274386990252, 10.129997786648358, 115.03410106358012, 120.59809824278094, 28.445128655074825, 24.804282183490045, 45.60001329839841, 54.379135472663855, 100.36656776423267, 7.1822472978541665, 69.40832478912773, 23.052443759982285, 27.523937880857577, 8.04043585859457, 33.09133066645042, 36.03935778629557, 24.983223811319892, 29.60588377718645, 11.675196360863536, 5.603972321785733, 85.86093185281321, 6.330764863829363, 12.86462784245834, 10.095862764191946, 54.41931073050113, 26.569248425808297, 16.057670851516743, 5.553189605060328, 35.91835265641494, 22.943831240828146, 21.979195849849486, 57.64208553402391, 39.86349284853823, 33.64598782450551, 153.69751159374601, 8.179707912007789, 65.19907640435628, 6.3731746240121545, 8.203519165503396, 18.461192792801985, 9.378688151242594, 31.601214988486028, 17.527842294205264, 15.502400387266192, 15.759471304195365, 192.94054873959666, 81.66392771040061, 31.28089465363296, 7.897973252990759, 76.80174422391055, 12.998234122218102, 27.228943722390547, 10.255578470981368, 8.344127005302429, 31.227121693565774, 112.64753732933204, 84.57942169022655, 6.755584123921043, 129.6651810321321, 35.06480684977999, 10.725893355668381, 119.5733993740958, 15.437543296927409, 5.422740319555255, 5.111246949472533, 32.18588515908512, 40.77242655841336, 90.25513639741051, 7.280682639352037, 77.62422027098206, 22.99121109848209, 67.72342756292294, 16.276133475377677, 32.102843183515155, 10.832739398097864, 66.0874783140498, 57.8059826646064, 19.403919895286432, 90.47193379697455, 9.934088548358066, 5.297742092654492, 74.81000662410375, 14.028357983947314, 6.503026641515817, 38.110746131507874, 31.19567044194315, 28.85526189513476, 46.68178040088819, 25.370980019490304, 29.74618510053752, 5.812437935303385, 114.34868919062774, 28.228757623061927, 41.81577428548081, 57.84271723124805, 22.92408061043053, 27.165065377009284, 40.352581437828285, 56.09387447922988, 50.64836482406433, 7.555244432039499, 6.579163328083184, 17.63436454018527, 49.06731276152528, 45.18078872034652, 22.259135349107297, 16.029458217998748, 51.05830677514964, 11.939508329800633, 6.61063322611364, 11.814231171268109, 10.742024049354702, 13.755892036566769, 14.541561349340286, 42.66995775086826, 11.323124212455566, 68.19026545077588, 24.057057603169547, 107.2010066940721, 29.38064581138545, 5.238439750192551, 15.38065479939424, 26.428753066981713, 100.92082851566539, 18.95662968989847, 106.32004525234811, 32.54977997824817, 49.40870064583637, 13.029481349876233, 10.678344216231519, 45.35456948327778, 65.48779644159757, 24.123998075217237, 69.4666661381233, 5.727448883501773, 19.325350682528594, 72.44816343599858, 42.5988272568806, 73.89845024180056, 51.98895013992934, 9.491066001418263, 46.387717662298286, 7.222415955323188, 34.48454880254624, 65.68751152076723, 5.211427768320859, 27.440242615085037, 39.39707298326284, 9.6765657827956, 14.331049459827428, 13.246368053407936, 15.686035953566362, 20.941839106548784, 10.720418571356598, 13.845165038436733, 7.113718941942257, 63.05466688401277, 69.59607448694868, 29.281196176750544, 8.891762549465556, 5.364152630269463, 131.94940676017814, 103.58620754774593, 51.84897823387635, 13.318800262476152, 8.052412111664024, 13.150209584174325, 62.09884185938603, 9.312112045567218, 53.60266778314205, 15.492859181962173, 6.720785975304419, 15.336933400850091, 100.80972994565191, 40.125024789834015, 97.35213209262277, 79.17677888834599, 15.434523461552672, 72.41510840814257, 76.1290006632271, 25.95245345017412, 25.946580354324322, 67.39171313697904, 84.67664208775805, 7.883125964175795, 25.701161924387222, 15.0305663848643, 65.61402769362405, 18.714518515639014, 7.884252759880994, 66.12187057903228, 7.7422409944461945, 7.930287508880684, 33.63132475150989, 83.28621539989015, 28.71326309822679, 13.555599003069942, 11.198195616775747, 28.83893870711918, 71.99012084775593, 34.27059600174316, 15.491059102396095, 94.85796046303663, 60.80621641633502, 9.298607510593195, 24.418031367950462, 43.07865403523679, 58.993156938069575, 23.980590210919814, 22.608508374697323, 18.85112564381806, 76.78964369657797, 55.06861024881866, 13.293523150838999, 63.095268398131324, 56.17066112459775, 35.573440325706656, 35.01348108909672, 9.488463612785868, 79.62241535886976, 13.769130029172004, 25.56633873415291, 46.014920023969054, 37.90023427118021, 47.74529842289515, 14.821632874934174, 81.4225354428202, 65.16212021741647, 9.3581247125086, 130.8053047235203, 33.23626145194534, 67.33242632089512, 14.458149277836041, 25.679156622015196, 86.3978634807841, 75.85960454048146, 14.125211046946713, 77.32937189092672, 21.577605834222336, 17.301057635941863, 50.93057878173799, 7.845349388155738, 12.252503029805252, 5.8840402398967235, 15.690595495616433, 60.627401180799644, 7.976440232600355, 7.357624827285209, 52.26075744299292, 7.634449117748097, 8.493106532176006, 12.378104995582737])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([7061593.75, 7094107.8125, 7136457.8125, 7165618.75, 7314129.6875, 7369234.375, 7379143.75, 7475423.4375, 7507050.0, 7516884.375, 7516917.1875, 7532898.4375, 7545882.8125, 7665126.5625, 7731723.4375, 7740578.125, 7791290.625, 7796645.3125, 7894650.0, 8063903.125, 8064143.75, 8081240.625, 8149629.6875, 8173607.8125, 8196509.375, 8296281.25, 8298535.9375, 8332868.75, 8335373.4375, 8345445.3125, 8357275.0, 8381796.875, 8384681.25, 8388559.375, 8397385.9375, 8437132.8125, 8444285.9375, 8449117.1875, 8457600.0, 8507312.5, 8540082.8125, 8546862.5, 8572471.875, 8595956.25, 8611620.3125, 8621629.6875, 8648803.125, 8692535.9375, 8801998.4375, 8869753.125, 8869776.5625, 8870157.8125, 8892581.25, 8909712.5, 8912267.1875, 8924081.25, 8937540.625, 8944387.5, 9077096.875, 9079775.0, 9080823.4375, 9105562.5, 9114835.9375, 9126712.5, 9145242.1875, 9164987.5, 9168923.4375, 9181139.0625, 9184987.5, 9190543.75, 9212739.0625, 9226307.8125, 9227689.0625, 9228715.625, 9235398.4375, 9242964.0625, 9283775.0, 9307321.875, 9308928.125, 9308971.875, 9322628.125, 9330110.9375, 9372259.375, 9387946.875, 9390525.0, 9396026.5625, 9420315.625, 9428125.0, 9431637.5, 9432137.5, 9432303.125, 9432321.875, 9435570.3125, 9468221.875, 9471503.125, 9471971.875, 9528312.5, 9538973.4375, 9602695.3125, 9636742.1875, 9663726.5625, 9717195.3125, 9762660.9375, 9762712.5, 9770528.125, 9815509.375, 9819387.5, 9837284.375, 9862785.9375, 9867406.25, 9940407.8125, 9987234.375, 13897032.8125, 13903092.1875, 13928628.125, 13936584.375, 13949195.3125, 13995673.4375, 14015260.9375, 14019679.6875, 14022262.5, 14023420.3125, 14023457.8125, 14024776.5625, 14026320.3125, 14027896.875, 14027903.125, 14027970.3125, 14027996.875, 14030532.8125, 14031292.1875, 14033853.125, 14037579.6875, 14038510.9375, 14040364.0625, 14043387.5, 14048459.375, 14050093.75, 14050621.875, 14050696.875, 14050860.9375, 14052976.5625, 14053023.4375, 14054107.8125, 14054814.0625, 14055162.5, 14057601.5625, 14062215.625, 14062654.6875, 14063395.3125, 14073559.375, 14077432.8125, 14079910.9375, 14082051.5625, 14083595.3125, 14087943.75, 14090359.375, 14091442.1875, 14091701.5625, 14092000.0, 14092485.9375, 14092987.5, 14093596.875, 14096893.75, 14098607.8125, 14099900.0, 14104860.9375, 14104918.75, 14106817.1875, 14107292.1875, 14107464.0625, 14107500.0, 14107642.1875, 14108034.375, 14108104.6875, 14108226.5625, 14108246.875, 14108579.6875, 14108695.3125, 14108748.4375, 14108998.4375, 14109448.4375, 14109643.75, 14109675.0, 14109957.8125, 14110004.6875, 14110053.125, 14110567.1875, 14110571.875, 14110729.6875, 14111112.5, 14111128.125, 14111203.125, 14111695.3125, 14111703.125, 14112042.1875, 14112104.6875, 14112287.5, 14112400.0, 14112401.5625, 14112532.8125, 14113296.875, 14113303.125, 14113503.125, 14114776.5625, 14114965.625, 14115142.1875, 14116695.3125, 14117165.625, 14118179.6875, 14118445.3125, 14118448.4375, 14119423.4375, 14119635.9375, 14119814.0625, 14119906.25, 14120198.4375, 14120276.5625, 14120325.0, 14120854.6875, 14120954.6875, 14121303.125, 14121356.25, 14122217.1875, 14122653.125, 14122839.0625, 14122981.25, 14123003.125, 14123051.5625, 14123315.625, 14123360.9375, 14123365.625, 14123476.5625, 14124060.9375, 14124554.6875, 14124560.9375, 14124871.875, 14124940.625, 14125478.125, 14125637.5, 14125768.75, 14126159.375, 14126185.9375, 14126329.6875, 14126412.5, 14126518.75, 14126742.1875, 14126865.625, 14127070.3125, 14127175.0, 14127382.8125, 14127995.3125, 14128059.375, 14128065.625, 14128264.0625, 14128284.375, 14128451.5625, 14129212.5, 14129837.5, 14129945.3125, 14130034.375, 14130248.4375, 14130300.0, 14130306.25, 14130351.5625, 14130531.25, 14130751.5625, 14131023.4375, 14131317.1875, 14131418.75, 14132281.25, 14132332.8125, 14132343.75, 14132353.125, 14132528.125, 14132748.4375, 14133131.25, 14133604.6875, 14133862.5, 14134568.75, 14134798.4375, 14135043.75, 14135828.125, 14135989.0625, 14136137.5, 14136423.4375, 14136735.9375, 14136982.8125, 14137323.4375, 14137462.5, 14137632.8125, 14137753.125, 14137996.875, 14138070.3125, 14138200.0, 14138242.1875, 14138253.125, 14138320.3125, 14138354.6875, 14138454.6875, 14138465.625, 14138535.9375, 14138562.5, 14138801.5625, 14138868.75, 14139028.125, 14139065.625, 14139078.125, 14139262.5, 14139335.9375, 14139535.9375, 14139651.5625, 14139687.5, 14140417.1875, 14140503.125, 14140587.5, 14140775.0, 14140829.6875, 14140850.0, 14140882.8125, 14140992.1875, 14141443.75, 14142292.1875, 14142407.8125, 14142443.75, 14142807.8125, 14142812.5, 14143096.875, 14143157.8125, 14143176.5625, 14143240.625, 14143267.1875, 14143450.0, 14143456.25, 14143457.8125, 14143484.375, 14143887.5, 14144118.75, 14144193.75, 14144375.0, 14144512.5, 14144529.6875, 14144742.1875, 14144743.75, 14145101.5625, 14145293.75, 14145401.5625, 14145543.75, 14145615.625, 14146081.25, 14146143.75, 14146489.0625, 14146518.75, 14146684.375, 14146764.0625, 14146810.9375, 14146868.75, 14147015.625, 14147353.125, 14147364.0625, 14147365.625, 14147634.375, 14147801.5625, 14147831.25, 14148221.875, 14148665.625, 14148689.0625, 14148896.875, 14148960.9375, 14148968.75, 14148995.3125, 14149040.625, 14149048.4375, 14149171.875, 14149173.4375, 14149204.6875, 14149210.9375, 14149239.0625, 14149292.1875, 14149407.8125, 14149548.4375, 14149626.5625, 14149679.6875, 14149842.1875, 14149842.1875, 14149873.4375, 14149875.0, 14149917.1875, 14149923.4375, 14149931.25, 14150006.25, 14150023.4375, 14150057.8125, 14150059.375, 14150171.875, 14150210.9375, 14150231.25, 14150240.625, 14150279.6875, 14150282.8125, 14150287.5, 14150440.625, 14150471.875, 14150484.375, 14150520.3125, 14150589.0625, 14150615.625, 14150646.875, 14150664.0625, 14150676.5625, 14150723.4375, 14150837.5, 14150920.3125, 14150945.3125, 14151031.25, 14151081.25, 14151201.5625, 14151245.3125, 14151351.5625, 14151476.5625, 14151514.0625, 14151557.8125, 14151609.375, 14151610.9375, 14151623.4375, 14151626.5625, 14151907.8125, 14152131.25, 14152246.875, 14152264.0625, 14152531.25, 14152603.125, 14152712.5, 14152737.5, 14152754.6875, 14152889.0625, 14152982.8125, 14153065.625, 14153101.5625, 14153160.9375, 14153220.3125, 14153245.3125, 14153257.8125, 14153278.125, 14153296.875, 14153523.4375, 14153556.25, 14153709.375, 14153781.25, 14153887.5, 14153992.1875, 14154017.1875, 14154028.125, 14154029.6875, 14154053.125, 14154098.4375, 14154131.25, 14154167.1875, 14154346.875, 14154371.875, 14154412.5, 14154442.1875, 14154468.75, 14154523.4375, 14154728.125, 14154745.3125, 14154840.625, 14155045.3125, 14155142.1875, 14155357.8125, 14155379.6875, 14155414.0625, 14155490.625, 14155507.8125, 14155603.125, 14155615.625, 14155620.3125, 14155671.875, 14155682.8125, 14155704.6875, 14155706.25, 14155707.8125, 14155725.0, 14155745.3125, 14155765.625, 14155767.1875, 14155778.125, 14155781.25, 14155853.125, 14155989.0625, 14156073.4375, 14156081.25, 14156114.0625, 14156339.0625, 14156535.9375, 14156575.0, 14156714.0625, 14156725.0, 14156839.0625, 14156851.5625, 14156985.9375, 14157046.875, 14157051.5625, 14157053.125, 14157354.6875, 14157404.6875, 14157528.125, 14157545.3125, 14157750.0, 14157906.25, 14157981.25, 14158114.0625, 14158123.4375, 14158150.0, 14158296.875, 14158418.75, 14158445.3125, 14158450.0, 14158459.375, 14158490.625, 14158514.0625, 14158632.8125, 14158714.0625, 14158742.1875, 14158765.625, 14158856.25, 14159012.5, 14159321.875, 14159334.375, 14159495.3125, 14159534.375, 14159546.875, 14159564.0625, 14159585.9375, 14159593.75, 14159715.625, 14159753.125, 14159845.3125, 14159854.6875, 14159971.875, 14159992.1875, 14160057.8125, 14160104.6875, 14160114.0625, 14160181.25, 14160351.5625, 14160354.6875, 14160406.25, 14160451.5625, 14160476.5625, 14160726.5625, 14160832.8125, 14160848.4375, 14160876.5625, 14160895.3125, 14160976.5625, 14161010.9375, 14161070.3125, 14161079.6875, 14161150.0, 14161150.0, 14161203.125, 14161239.0625, 14161260.9375, 14161278.125, 14161300.0, 14161309.375, 14161310.9375, 14161403.125, 14161448.4375, 14161462.5, 14161482.8125, 14161489.0625, 14161501.5625, 14161546.875, 14161551.5625, 14161557.8125, 14161568.75, 14161596.875, 14161606.25, 14161642.1875, 14161651.5625, 14161660.9375, 14161743.75, 14161832.8125, 14161857.8125, 14162048.4375, 14162209.375, 14162246.875, 14162329.6875, 14162453.125, 14162492.1875, 14162517.1875, 14162554.6875, 14162623.4375, 14162646.875, 14162690.625, 14162765.625, 14162942.1875, 14162993.75, 14162996.875, 14163006.25, 14163035.9375, 14163151.5625, 14163176.5625, 14163353.125, 14163376.5625, 14163400.0, 14163417.1875, 14163512.5, 14163551.5625, 14163609.375, 14163787.5, 14163898.4375, 14163903.125, 14163906.25, 14163979.6875, 14163981.25, 14163992.1875, 14163996.875, 14164181.25, 14164562.5, 14164750.0, 14164784.375, 14164821.875, 14164829.6875, 14164962.5, 14165006.25, 14165042.1875, 14165078.125, 14165243.75, 14165342.1875, 14165428.125, 14165503.125, 14165590.625, 14165609.375, 14165635.9375, 14165660.9375, 14165681.25, 14165689.0625, 14165709.375, 14165740.625, 14165760.9375, 14165765.625, 14165776.5625, 14165785.9375, 14165804.6875, 14165818.75, 14165859.375, 14165898.4375, 14165907.8125, 14165910.9375, 14165929.6875, 14165976.5625, 14165978.125, 14166018.75, 14166132.8125, 14166157.8125, 14166246.875, 14166440.625, 14166603.125, 14166671.875, 14166792.1875, 14166901.5625, 14166943.75, 14167029.6875, 14167082.8125, 14167084.375, 14167142.1875, 14167143.75, 14167160.9375, 14167364.0625, 14167459.375, 14167664.0625, 14167792.1875, 14167800.0, 14167812.5, 14168106.25, 14168115.625, 14168290.625, 14168300.0, 14168579.6875, 14168857.8125, 14169046.875, 14169220.3125, 14169687.5, 14169921.875, 14169996.875, 14170292.1875, 14170348.4375, 14170384.375, 14171021.875, 14171310.9375, 14171335.9375, 14171637.5, 14171764.0625, 14171826.5625, 14171915.625, 14172003.125, 14172129.6875, 14172282.8125, 14172326.5625, 14172656.25, 14172676.5625, 14172878.125, 14173100.0, 14173310.9375, 14173646.875, 14173726.5625, 14173798.4375, 14174464.0625, 14174809.375, 14174962.5, 14174967.1875, 14174976.5625, 14175090.625, 14175114.0625, 14175153.125, 14175573.4375, 14175637.5, 14175790.625, 14176089.0625, 14176718.75, 14176718.75, 14176867.1875, 14177201.5625, 14177900.0, 14178412.5, 14178857.8125, 14178859.375, 14178993.75, 14179395.3125, 14179567.1875, 14180114.0625, 14180900.0, 14181170.3125, 14182025.0, 14182746.875, 14183140.625, 14184557.8125, 14184934.375, 14184971.875, 14184990.625, 14185576.5625, 14185637.5, 14185712.5, 14185879.6875, 14185942.1875, 14186414.0625, 14186792.1875, 14187131.25, 14187528.125, 14188085.9375, 14188307.8125, 14188537.5, 14188575.0, 14188901.5625, 14189168.75, 14189495.3125, 14189845.3125, 14190823.4375, 14191256.25, 14191310.9375, 14191317.1875, 14192559.375, 14192592.1875, 14193084.375, 14193559.375, 14194429.6875, 14196209.375, 14197817.1875, 14199010.9375, 14199012.5, 14199107.8125, 14199748.4375, 14199756.25, 14200703.125, 14201914.0625, 14202307.8125, 14202560.9375, 14202989.0625, 14203189.0625, 14203348.4375, 14203595.3125, 14203756.25, 14204378.125, 14205239.0625, 14205481.25, 14205678.125, 14205978.125, 14206098.4375, 14206164.0625, 14206626.5625, 14207867.1875, 14210195.3125, 14210546.875, 14211262.5, 14212501.5625, 14214840.625, 14215814.0625, 14215979.6875, 14215981.25, 14216007.8125, 14216457.8125, 14216664.0625, 14221175.0, 14221579.6875, 14223371.875, 14224518.75, 14224739.0625, 14224815.625, 14225807.8125, 14225965.625, 14226720.3125, 14230601.5625, 14232651.5625, 14233571.875, 14234320.3125, 14234446.875, 14240440.625, 14240776.5625, 14243726.5625, 14245967.1875, 14248412.5, 14249862.5, 14250431.25, 14256106.25, 14262431.25, 14270289.0625, 14271584.375, 14282850.0, 14283479.6875, 14288748.4375, 14288756.25, 14316635.9375, 442749226.5625], [14.353494663391347, 473.72552766437803, 64.83849779080148, 36.77691193062272, 45.5215314765082, 9.35358762473517, 138.96490170140314, 10.736181369094693, 63.82109119366702, 55.938173334457915, 88.7681590638431, 19.734325713764772, 35.651572755161176, 46.550626326252385, 20.871126056680342, 21.84359957165835, 47.23218286247899, 6.13684904713674, 9.24638261305465, 30.44021871282233, 47.03775059633008, 6.55245419575673, 58.702711176285426, 23.636303341479437, 18.179892111916878, 89.94716101189168, 5.1297197217189465, 12.358679632605428, 10.71630147318415, 93.0766903499036, 17.52061275270044, 24.991509180526812, 20.46416628164172, 20.60424365329352, 9.121317134698764, 33.91656963226593, 19.029721329047504, 51.49528148779245, 88.70332418188231, 82.62434546059671, 5.2298468191344565, 58.33207553908103, 41.93331390532528, 86.784007584458, 51.30684960822346, 8.61703840863088, 32.63254148345757, 43.17868874281369, 12.753193533769785, 34.34814232164087, 10.167794305449949, 87.83604021181014, 39.235160788821744, 18.255997217943104, 6.01336053078085, 11.0482600586928, 32.090075499029574, 68.50731007923501, 39.51983998604173, 29.065530742111154, 25.36276690784278, 5.921981987503392, 84.37475442262945, 11.488771362701042, 50.02737375445998, 22.722740087727537, 17.169690837376297, 91.45494429636368, 21.335876551298693, 6.5764150127191465, 111.03219718051129, 85.299902567554, 23.945056982174265, 38.52536228253169, 28.972469800237363, 66.00821662089689, 13.676826791699021, 72.9907585653588, 95.70378321655397, 62.12037854049055, 66.7891484470813, 54.47259835167837, 27.28921849345727, 39.35013358710248, 28.2602832444484, 65.7435275050917, 12.801872069724844, 25.61020355370054, 16.63461569573696, 46.28079114633559, 96.29792409196753, 89.07687421935228, 44.22610227911099, 14.941244148153501, 13.057208369360833, 12.69087761508895, 38.26527543577155, 107.25872188508636, 35.72856526079983, 58.60820293278006, 7.372982563358372, 26.732792541294373, 68.17684566975235, 5.354712784245878, 64.41885900111, 66.85545750295871, 28.996827169854154, 112.19906756178189, 88.49173997056705, 28.492001388069177, 14.434911830894555, 14.534011182788714, 80.34683665108685, 29.485995613491788, 10.977925767294726, 105.69469775692211, 17.168300381557057, 20.036805185614934, 6.328084087484737, 39.78822185923844, 22.85319254143005, 44.848943485916244, 51.88871040089554, 20.871193667039336, 26.467093006754947, 15.410974375179872, 89.7989222877052, 5.1492059434781, 5.387654290123565, 16.536214208219764, 26.62462523961121, 33.20454150504682, 10.211351191310994, 19.601035348003556, 27.953888009990216, 140.85742704841465, 87.5724959221242, 7.7884728836066905, 17.888564589548565, 9.29390012359141, 27.901629567809184, 61.82385221825618, 69.82649026251953, 8.898215595327333, 95.92870616255469, 35.27538308689147, 18.52547826745096, 19.67910539562954, 77.94228142861903, 5.71437746038673, 44.47637852710898, 56.36949227037983, 29.022241895616872, 19.432695983014664, 26.510717661558907, 25.756896138374913, 11.326472613546093, 65.36063868588136, 50.949582862337884, 19.770498431721553, 13.538392312171128, 173.75082529786062, 42.85044921383103, 14.150076505973443, 113.0498348797251, 28.173464299292945, 7.611891421922517, 5.09488032341785, 35.222522822861414, 19.14699249790841, 58.64371926551527, 60.664350018351925, 23.22640493206027, 6.795742014943427, 68.45664886767102, 6.519365833077661, 81.72007403771516, 26.181263879217074, 17.219910229944606, 32.86166634729164, 44.27113433493194, 8.566632772941917, 7.517428094158713, 60.515464745670485, 72.78434080480221, 28.199652866668668, 8.960780437403438, 21.635228861273784, 14.910618693956629, 29.148818182271953, 5.216270756522859, 81.2912865049599, 46.614024831305414, 67.18256762792173, 15.867553571281848, 10.270541321095285, 73.82180026059304, 20.058651924123808, 37.071630126810064, 58.4809364584972, 22.550319185008476, 6.933272729737776, 25.587246155338363, 73.06943091049588, 97.59715330110649, 5.712836102823359, 83.06397033333556, 47.11838201349126, 25.303326137019152, 40.30559279508621, 11.340039130667057, 29.153939809226653, 64.22302227576327, 36.7537455405884, 17.5776176592042, 8.118588397117676, 29.570597285957536, 25.396005680513376, 18.822971165666523, 36.453855472548426, 18.710383379088043, 57.55298437809341, 29.590219543437573, 48.267791310782805, 46.37818819078264, 25.73942578247345, 29.593005254134496, 58.61179312681341, 5.744956007708263, 78.00161546103493, 11.199893235750558, 22.72205252054634, 97.31172513930277, 32.675650253293114, 80.89593158843272, 50.725112666314686, 21.639503016235178, 5.234156784139545, 21.2604218195829, 31.40272912710506, 38.89152686193783, 23.576606320267587, 24.12015277532069, 80.56142884208839, 61.783822976485546, 130.96431079441888, 7.225273460899372, 25.27921487972804, 9.823766914230644, 12.560017545618122, 58.076908098202786, 41.12700090295604, 76.07454717282067, 150.85943632080748, 27.65022827389516, 5.57519222573279, 11.76680217485271, 5.640250335941951, 38.65796054701666, 47.792779484209134, 5.357885297869613, 26.029438017087415, 29.604191167731127, 42.795952554645254, 68.58387361761261, 5.901564481523938, 11.45765452471494, 10.330741219667953, 115.23435989695996, 25.440892274317907, 43.92058086375303, 14.68165371458611, 77.07678788115845, 72.06944552585097, 10.3784304246614, 20.516357260276195, 29.66975523850248, 37.390459818246526, 45.10135817742194, 7.822877740003752, 6.9047549546724545, 69.02149245526022, 115.91437765085564, 18.29596293877664, 91.04978680844967, 19.987977949081962, 13.58785447880324, 32.631274458799595, 122.34863721889519, 8.687979854997394, 67.43830196743917, 19.796190578387257, 9.70874289866778, 5.3208183280715105, 26.166139992579883, 52.352664071284345, 65.78788479758539, 32.96582134133165, 22.31088865752372, 53.04019828907692, 7.6305856592876875, 68.29311692352044, 5.198107194523418, 14.925240601590707, 70.30190856292792, 15.257973668040526, 24.523792994991883, 12.128948368450294, 61.880366493780315, 64.35540451634562, 56.12343846276438, 6.036973747837919, 27.986491478131114, 9.635269287039392, 15.384076907805271, 6.379008749669283, 33.151920247776644, 16.08108154619088, 45.126872121236985, 14.26252348097054, 21.745799117344784, 66.19139995038633, 10.092529100972401, 32.880980227759096, 59.18544360724054, 62.37415011395015, 23.700612120865156, 39.22478817111433, 12.128612534526674, 46.7952996734885, 17.240567391596084, 19.563062491699153, 32.77943073174624, 75.9244532855812, 6.2623924026738536, 27.785666873615423, 6.0800116185134, 24.544525550550887, 190.73086764807385, 68.7782447723781, 6.76145788477542, 88.4272289875377, 32.31481695488803, 28.392484217973628, 9.289049871627414, 35.26804737796784, 162.9244149185386, 20.013158130557656, 12.28036558966326, 181.68393035370133, 104.97754878116041, 10.516428442332703, 21.81040952229024, 16.018448455714275, 18.646876630704586, 11.558945118318242, 34.64340098544084, 43.12197024156084, 26.711302433437393, 57.47291746726988, 9.971635292758434, 13.454504760957123, 42.061852668240164, 65.88109421558121, 5.839771119107743, 58.50996494348714, 14.806594037032847, 9.5681734380959, 79.75198020536143, 34.78258862659328, 12.586082027440696, 27.273693806676206, 23.35680763294256, 11.284664052981473, 57.65542306537539, 100.14971778596654, 31.78167166919605, 45.64889588083866, 38.65420264274514, 79.66758193074546, 17.1201419980156, 102.9097270235377, 17.219748229071147, 5.368158579478964, 31.601807003765117, 71.76241143826363, 32.939004211533664, 42.155056073147755, 11.987017610004475, 9.748965266877068, 19.767290591809402, 7.110682854338585, 7.392816442907737, 24.204429172715763, 91.49244085784125, 8.293983059766695, 23.937622027379852, 14.20480637365656, 9.68367346001704, 85.73824785275004, 60.35682557482924, 133.2052298197641, 11.943339722888041, 72.69601261208908, 120.63466956028064, 91.80284339144525, 9.815860612181233, 6.559279010428072, 33.564597357816474, 9.555392088807762, 78.3831574207153, 27.179852850244465, 6.594196462714507, 40.75154635054125, 36.724422283613656, 19.247476621369373, 37.929071005299356, 17.703290882445362, 107.20111832155399, 120.53668625294208, 16.841635354393926, 6.1544185316873605, 135.75325572677988, 74.31783458236261, 47.199725363220104, 49.58537836937867, 57.278543193439255, 53.959106161268366, 77.60096407933997, 58.56103256615294, 144.74222144461513, 85.70036375992598, 176.41859715798293, 11.115712889705602, 54.354030798643635, 30.350359721749584, 6.293897324225599, 55.71879324120669, 29.430641341484563, 107.56090829397996, 16.845050039553357, 40.235935353644614, 22.69808557703132, 76.87226004600889, 15.579405560730098, 10.010674334203065, 19.290120177564052, 38.52634537582876, 6.997355117862078, 13.638581184715129, 48.785429134526986, 42.00743014501606, 101.61683499211131, 73.20579185488225, 27.091800585534948, 30.619262391062065, 42.29559070599094, 6.137766473903482, 10.385163457992718, 43.24332129557821, 8.060733604365055, 50.43517550739337, 18.16350169807365, 9.86721091156452, 11.406037867396963, 17.72729173742572, 103.40576469931088, 75.34640461928947, 5.7576106600359225, 30.240818891386503, 37.46407154749019, 11.442331310233499, 79.65797703991664, 51.40732769005136, 65.57899699966401, 20.5501338259211, 29.623393177108042, 25.338588503117435, 128.6027701241257, 58.325084278925424, 47.76298990475915, 17.595070019309876, 9.867098645014597, 132.67222661305547, 8.921745451046558, 50.526320092458434, 51.58389783260908, 32.022575343000476, 31.963147853396553, 71.49294496650889, 105.55434234269504, 112.65856204740324, 29.25502346049976, 171.86869265024822, 58.31970632752492, 49.932414663227775, 80.57132976094407, 65.61759374109833, 23.22642695072207, 166.34354258819423, 27.542601311789923, 25.772105636994027, 63.93691962545487, 66.90935078348784, 106.14916392544066, 10.939203343101395, 42.124368198549114, 10.844508684040633, 7.241840816275167, 6.73736242483435, 11.593779428282327, 18.335779657925514, 21.41947902885403, 18.753579045500953, 19.35315697538517, 107.13513831793823, 31.550913024453926, 102.6628717816799, 35.88619177362325, 17.235990196897486, 24.141560673491583, 22.806446988252, 62.00804244895558, 21.8257731318311, 8.383265822904729, 11.875105911471099, 6.636418077090405, 5.961896993654193, 12.40819180158402, 72.92229378544869, 16.950758129527294, 12.532837565835399, 21.40140000993174, 5.978316730841685, 163.52628580424195, 18.29264744716334, 9.53987319215665, 23.17433371053665, 27.759873069797564, 25.77991529317198, 35.19983605801254, 38.35084678260363, 10.69775915491015, 15.950459283085388, 16.93716050207172, 13.2905645636859, 37.018675301060775, 27.618674979992512, 57.673680193818264, 60.50167468273238, 23.631718542799774, 57.29710130082808, 26.36064094057911, 15.111439999683835, 32.86672071402299, 7.4691393743382815, 60.48189238985312, 79.35631343185365, 24.277194971053582, 34.47297681181391, 42.79701611969088, 53.08722263706735, 8.30078276625573, 42.66492192149025, 72.4143974575048, 23.492589951121623, 24.503793887336872, 24.5654317567233, 73.56338891035281, 78.6167044099696, 22.832434501843665, 45.346258978838904, 27.835536129408673, 5.908811296577237, 9.81245982135465, 26.736256819297026, 19.553482018962708, 32.033837241875744, 22.219346359091624, 17.374203878533606, 9.454173248858073, 16.905984015553795, 14.604699809172235, 5.892264206598048, 24.750274386990252, 10.129997786648358, 115.03410106358012, 120.59809824278094, 28.445128655074825, 24.804282183490045, 45.60001329839841, 54.379135472663855, 100.36656776423267, 7.1822472978541665, 69.40832478912773, 23.052443759982285, 27.523937880857577, 8.04043585859457, 33.09133066645042, 36.03935778629557, 24.983223811319892, 29.60588377718645, 11.675196360863536, 5.603972321785733, 85.86093185281321, 6.330764863829363, 12.86462784245834, 10.095862764191946, 54.41931073050113, 26.569248425808297, 16.057670851516743, 5.553189605060328, 35.91835265641494, 22.943831240828146, 21.979195849849486, 57.64208553402391, 39.86349284853823, 33.64598782450551, 153.69751159374601, 8.179707912007789, 65.19907640435628, 6.3731746240121545, 8.203519165503396, 18.461192792801985, 9.378688151242594, 31.601214988486028, 17.527842294205264, 15.502400387266192, 15.759471304195365, 192.94054873959666, 81.66392771040061, 31.28089465363296, 7.897973252990759, 76.80174422391055, 12.998234122218102, 27.228943722390547, 10.255578470981368, 8.344127005302429, 31.227121693565774, 112.64753732933204, 84.57942169022655, 6.755584123921043, 129.6651810321321, 35.06480684977999, 10.725893355668381, 119.5733993740958, 15.437543296927409, 5.422740319555255, 5.111246949472533, 32.18588515908512, 40.77242655841336, 90.25513639741051, 7.280682639352037, 77.62422027098206, 22.99121109848209, 67.72342756292294, 16.276133475377677, 32.102843183515155, 10.832739398097864, 66.0874783140498, 57.8059826646064, 19.403919895286432, 90.47193379697455, 9.934088548358066, 5.297742092654492, 74.81000662410375, 14.028357983947314, 6.503026641515817, 38.110746131507874, 31.19567044194315, 28.85526189513476, 46.68178040088819, 25.370980019490304, 29.74618510053752, 5.812437935303385, 114.34868919062774, 28.228757623061927, 41.81577428548081, 57.84271723124805, 22.92408061043053, 27.165065377009284, 40.352581437828285, 56.09387447922988, 50.64836482406433, 7.555244432039499, 6.579163328083184, 17.63436454018527, 49.06731276152528, 45.18078872034652, 22.259135349107297, 16.029458217998748, 51.05830677514964, 11.939508329800633, 6.61063322611364, 11.814231171268109, 10.742024049354702, 13.755892036566769, 14.541561349340286, 42.66995775086826, 11.323124212455566, 68.19026545077588, 24.057057603169547, 107.2010066940721, 29.38064581138545, 5.238439750192551, 15.38065479939424, 26.428753066981713, 100.92082851566539, 18.95662968989847, 106.32004525234811, 32.54977997824817, 49.40870064583637, 13.029481349876233, 10.678344216231519, 45.35456948327778, 65.48779644159757, 24.123998075217237, 69.4666661381233, 5.727448883501773, 19.325350682528594, 72.44816343599858, 42.5988272568806, 73.89845024180056, 51.98895013992934, 9.491066001418263, 46.387717662298286, 7.222415955323188, 34.48454880254624, 65.68751152076723, 5.211427768320859, 27.440242615085037, 39.39707298326284, 9.6765657827956, 14.331049459827428, 13.246368053407936, 15.686035953566362, 20.941839106548784, 10.720418571356598, 13.845165038436733, 7.113718941942257, 63.05466688401277, 69.59607448694868, 29.281196176750544, 8.891762549465556, 5.364152630269463, 131.94940676017814, 103.58620754774593, 51.84897823387635, 13.318800262476152, 8.052412111664024, 13.150209584174325, 62.09884185938603, 9.312112045567218, 53.60266778314205, 15.492859181962173, 6.720785975304419, 15.336933400850091, 100.80972994565191, 40.125024789834015, 97.35213209262277, 79.17677888834599, 15.434523461552672, 72.41510840814257, 76.1290006632271, 25.95245345017412, 25.946580354324322, 67.39171313697904, 84.67664208775805, 7.883125964175795, 25.701161924387222, 15.0305663848643, 65.61402769362405, 18.714518515639014, 7.884252759880994, 66.12187057903228, 7.7422409944461945, 7.930287508880684, 33.63132475150989, 83.28621539989015, 28.71326309822679, 13.555599003069942, 11.198195616775747, 28.83893870711918, 71.99012084775593, 34.27059600174316, 15.491059102396095, 94.85796046303663, 60.80621641633502, 9.298607510593195, 24.418031367950462, 43.07865403523679, 58.993156938069575, 23.980590210919814, 22.608508374697323, 18.85112564381806, 76.78964369657797, 55.06861024881866, 13.293523150838999, 63.095268398131324, 56.17066112459775, 35.573440325706656, 35.01348108909672, 9.488463612785868, 79.62241535886976, 13.769130029172004, 25.56633873415291, 46.014920023969054, 37.90023427118021, 47.74529842289515, 14.821632874934174, 81.4225354428202, 65.16212021741647, 9.3581247125086, 130.8053047235203, 33.23626145194534, 67.33242632089512, 14.458149277836041, 25.679156622015196, 86.3978634807841, 75.85960454048146, 14.125211046946713, 77.32937189092672, 21.577605834222336, 17.301057635941863, 50.93057878173799, 7.845349388155738, 12.252503029805252, 5.8840402398967235, 15.690595495616433, 60.627401180799644, 7.976440232600355, 7.357624827285209, 52.26075744299292, 7.634449117748097, 8.493106532176006, 12.378104995582737])
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);
([7061593.75, 7094107.8125, 7136457.8125, 7165618.75, 7314129.6875, 7369234.375, 7379143.75, 7475423.4375, 7507050.0, 7516884.375, 7516917.1875, 7532898.4375, 7545882.8125, 7665126.5625, 7731723.4375, 7740578.125, 7791290.625, 7796645.3125, 7894650.0, 8063903.125, 8064143.75, 8081240.625, 8149629.6875, 8173607.8125, 8196509.375, 8296281.25, 8298535.9375, 8332868.75, 8335373.4375, 8345445.3125, 8357275.0, 8381796.875, 8384681.25, 8388559.375, 8397385.9375, 8437132.8125, 8444285.9375, 8449117.1875, 8457600.0, 8507312.5, 8540082.8125, 8546862.5, 8572471.875, 8595956.25, 8611620.3125, 8621629.6875, 8648803.125, 8692535.9375, 8801998.4375, 8869753.125, 8869776.5625, 8870157.8125, 8892581.25, 8909712.5, 8912267.1875, 8924081.25, 8937540.625, 8944387.5, 9077096.875, 9079775.0, 9080823.4375, 9105562.5, 9114835.9375, 9126712.5, 9145242.1875, 9164987.5, 9168923.4375, 9181139.0625, 9184987.5, 9190543.75, 9212739.0625, 9226307.8125, 9227689.0625, 9228715.625, 9235398.4375, 9242964.0625, 9283775.0, 9307321.875, 9308928.125, 9308971.875, 9322628.125, 9330110.9375, 9372259.375, 9387946.875, 9390525.0, 9396026.5625, 9420315.625, 9428125.0, 9431637.5, 9432137.5, 9432303.125, 9432321.875, 9435570.3125, 9468221.875, 9471503.125, 9471971.875, 9528312.5, 9538973.4375, 9602695.3125, 9636742.1875, 9663726.5625, 9717195.3125, 9762660.9375, 9762712.5, 9770528.125, 9815509.375, 9819387.5, 9837284.375, 9862785.9375, 9867406.25, 9940407.8125, 9987234.375, 13897032.8125, 13903092.1875, 13928628.125, 13936584.375, 13949195.3125, 13995673.4375, 14015260.9375, 14019679.6875, 14022262.5, 14023420.3125, 14023457.8125, 14024776.5625, 14026320.3125, 14027896.875, 14027903.125, 14027970.3125, 14027996.875, 14030532.8125, 14031292.1875, 14033853.125, 14037579.6875, 14038510.9375, 14040364.0625, 14043387.5, 14048459.375, 14050093.75, 14050621.875, 14050696.875, 14050860.9375, 14052976.5625, 14053023.4375, 14054107.8125, 14054814.0625, 14055162.5, 14057601.5625, 14062215.625, 14062654.6875, 14063395.3125, 14073559.375, 14077432.8125, 14079910.9375, 14082051.5625, 14083595.3125, 14087943.75, 14090359.375, 14091442.1875, 14091701.5625, 14092000.0, 14092485.9375, 14092987.5, 14093596.875, 14096893.75, 14098607.8125, 14099900.0, 14104860.9375, 14104918.75, 14106817.1875, 14107292.1875, 14107464.0625, 14107500.0, 14107642.1875, 14108034.375, 14108104.6875, 14108226.5625, 14108246.875, 14108579.6875, 14108695.3125, 14108748.4375, 14108998.4375, 14109448.4375, 14109643.75, 14109675.0, 14109957.8125, 14110004.6875, 14110053.125, 14110567.1875, 14110571.875, 14110729.6875, 14111112.5, 14111128.125, 14111203.125, 14111695.3125, 14111703.125, 14112042.1875, 14112104.6875, 14112287.5, 14112400.0, 14112401.5625, 14112532.8125, 14113296.875, 14113303.125, 14113503.125, 14114776.5625, 14114965.625, 14115142.1875, 14116695.3125, 14117165.625, 14118179.6875, 14118445.3125, 14118448.4375, 14119423.4375, 14119635.9375, 14119814.0625, 14119906.25, 14120198.4375, 14120276.5625, 14120325.0, 14120854.6875, 14120954.6875, 14121303.125, 14121356.25, 14122217.1875, 14122653.125, 14122839.0625, 14122981.25, 14123003.125, 14123051.5625, 14123315.625, 14123360.9375, 14123365.625, 14123476.5625, 14124060.9375, 14124554.6875, 14124560.9375, 14124871.875, 14124940.625, 14125478.125, 14125637.5, 14125768.75, 14126159.375, 14126185.9375, 14126329.6875, 14126412.5, 14126518.75, 14126742.1875, 14126865.625, 14127070.3125, 14127175.0, 14127382.8125, 14127995.3125, 14128059.375, 14128065.625, 14128264.0625, 14128284.375, 14128451.5625, 14129212.5, 14129837.5, 14129945.3125, 14130034.375, 14130248.4375, 14130300.0, 14130306.25, 14130351.5625, 14130531.25, 14130751.5625, 14131023.4375, 14131317.1875, 14131418.75, 14132281.25, 14132332.8125, 14132343.75, 14132353.125, 14132528.125, 14132748.4375, 14133131.25, 14133604.6875, 14133862.5, 14134568.75, 14134798.4375, 14135043.75, 14135828.125, 14135989.0625, 14136137.5, 14136423.4375, 14136735.9375, 14136982.8125, 14137323.4375, 14137462.5, 14137632.8125, 14137753.125, 14137996.875, 14138070.3125, 14138200.0, 14138242.1875, 14138253.125, 14138320.3125, 14138354.6875, 14138454.6875, 14138465.625, 14138535.9375, 14138562.5, 14138801.5625, 14138868.75, 14139028.125, 14139065.625, 14139078.125, 14139262.5, 14139335.9375, 14139535.9375, 14139651.5625, 14139687.5, 14140417.1875, 14140503.125, 14140587.5, 14140775.0, 14140829.6875, 14140850.0, 14140882.8125, 14140992.1875, 14141443.75, 14142292.1875, 14142407.8125, 14142443.75, 14142807.8125, 14142812.5, 14143096.875, 14143157.8125, 14143176.5625, 14143240.625, 14143267.1875, 14143450.0, 14143456.25, 14143457.8125, 14143484.375, 14143887.5, 14144118.75, 14144193.75, 14144375.0, 14144512.5, 14144529.6875, 14144742.1875, 14144743.75, 14145101.5625, 14145293.75, 14145401.5625, 14145543.75, 14145615.625, 14146081.25, 14146143.75, 14146489.0625, 14146518.75, 14146684.375, 14146764.0625, 14146810.9375, 14146868.75, 14147015.625, 14147353.125, 14147364.0625, 14147365.625, 14147634.375, 14147801.5625, 14147831.25, 14148221.875, 14148665.625, 14148689.0625, 14148896.875, 14148960.9375, 14148968.75, 14148995.3125, 14149040.625, 14149048.4375, 14149171.875, 14149173.4375, 14149204.6875, 14149210.9375, 14149239.0625, 14149292.1875, 14149407.8125, 14149548.4375, 14149626.5625, 14149679.6875, 14149842.1875, 14149842.1875, 14149873.4375, 14149875.0, 14149917.1875, 14149923.4375, 14149931.25, 14150006.25, 14150023.4375, 14150057.8125, 14150059.375, 14150171.875, 14150210.9375, 14150231.25, 14150240.625, 14150279.6875, 14150282.8125, 14150287.5, 14150440.625, 14150471.875, 14150484.375, 14150520.3125, 14150589.0625, 14150615.625, 14150646.875, 14150664.0625, 14150676.5625, 14150723.4375, 14150837.5, 14150920.3125, 14150945.3125, 14151031.25, 14151081.25, 14151201.5625, 14151245.3125, 14151351.5625, 14151476.5625, 14151514.0625, 14151557.8125, 14151609.375, 14151610.9375, 14151623.4375, 14151626.5625, 14151907.8125, 14152131.25, 14152246.875, 14152264.0625, 14152531.25, 14152603.125, 14152712.5, 14152737.5, 14152754.6875, 14152889.0625, 14152982.8125, 14153065.625, 14153101.5625, 14153160.9375, 14153220.3125, 14153245.3125, 14153257.8125, 14153278.125, 14153296.875, 14153523.4375, 14153556.25, 14153709.375, 14153781.25, 14153887.5, 14153992.1875, 14154017.1875, 14154028.125, 14154029.6875, 14154053.125, 14154098.4375, 14154131.25, 14154167.1875, 14154346.875, 14154371.875, 14154412.5, 14154442.1875, 14154468.75, 14154523.4375, 14154728.125, 14154745.3125, 14154840.625, 14155045.3125, 14155142.1875, 14155357.8125, 14155379.6875, 14155414.0625, 14155490.625, 14155507.8125, 14155603.125, 14155615.625, 14155620.3125, 14155671.875, 14155682.8125, 14155704.6875, 14155706.25, 14155707.8125, 14155725.0, 14155745.3125, 14155765.625, 14155767.1875, 14155778.125, 14155781.25, 14155853.125, 14155989.0625, 14156073.4375, 14156081.25, 14156114.0625, 14156339.0625, 14156535.9375, 14156575.0, 14156714.0625, 14156725.0, 14156839.0625, 14156851.5625, 14156985.9375, 14157046.875, 14157051.5625, 14157053.125, 14157354.6875, 14157404.6875, 14157528.125, 14157545.3125, 14157750.0, 14157906.25, 14157981.25, 14158114.0625, 14158123.4375, 14158150.0, 14158296.875, 14158418.75, 14158445.3125, 14158450.0, 14158459.375, 14158490.625, 14158514.0625, 14158632.8125, 14158714.0625, 14158742.1875, 14158765.625, 14158856.25, 14159012.5, 14159321.875, 14159334.375, 14159495.3125, 14159534.375, 14159546.875, 14159564.0625, 14159585.9375, 14159593.75, 14159715.625, 14159753.125, 14159845.3125, 14159854.6875, 14159971.875, 14159992.1875, 14160057.8125, 14160104.6875, 14160114.0625, 14160181.25, 14160351.5625, 14160354.6875, 14160406.25, 14160451.5625, 14160476.5625, 14160726.5625, 14160832.8125, 14160848.4375, 14160876.5625, 14160895.3125, 14160976.5625, 14161010.9375, 14161070.3125, 14161079.6875, 14161150.0, 14161150.0, 14161203.125, 14161239.0625, 14161260.9375, 14161278.125, 14161300.0, 14161309.375, 14161310.9375, 14161403.125, 14161448.4375, 14161462.5, 14161482.8125, 14161489.0625, 14161501.5625, 14161546.875, 14161551.5625, 14161557.8125, 14161568.75, 14161596.875, 14161606.25, 14161642.1875, 14161651.5625, 14161660.9375, 14161743.75, 14161832.8125, 14161857.8125, 14162048.4375, 14162209.375, 14162246.875, 14162329.6875, 14162453.125, 14162492.1875, 14162517.1875, 14162554.6875, 14162623.4375, 14162646.875, 14162690.625, 14162765.625, 14162942.1875, 14162993.75, 14162996.875, 14163006.25, 14163035.9375, 14163151.5625, 14163176.5625, 14163353.125, 14163376.5625, 14163400.0, 14163417.1875, 14163512.5, 14163551.5625, 14163609.375, 14163787.5, 14163898.4375, 14163903.125, 14163906.25, 14163979.6875, 14163981.25, 14163992.1875, 14163996.875, 14164181.25, 14164562.5, 14164750.0, 14164784.375, 14164821.875, 14164829.6875, 14164962.5, 14165006.25, 14165042.1875, 14165078.125, 14165243.75, 14165342.1875, 14165428.125, 14165503.125, 14165590.625, 14165609.375, 14165635.9375, 14165660.9375, 14165681.25, 14165689.0625, 14165709.375, 14165740.625, 14165760.9375, 14165765.625, 14165776.5625, 14165785.9375, 14165804.6875, 14165818.75, 14165859.375, 14165898.4375, 14165907.8125, 14165910.9375, 14165929.6875, 14165976.5625, 14165978.125, 14166018.75, 14166132.8125, 14166157.8125, 14166246.875, 14166440.625, 14166603.125, 14166671.875, 14166792.1875, 14166901.5625, 14166943.75, 14167029.6875, 14167082.8125, 14167084.375, 14167142.1875, 14167143.75, 14167160.9375, 14167364.0625, 14167459.375, 14167664.0625, 14167792.1875, 14167800.0, 14167812.5, 14168106.25, 14168115.625, 14168290.625, 14168300.0, 14168579.6875, 14168857.8125, 14169046.875, 14169220.3125, 14169687.5, 14169921.875, 14169996.875, 14170292.1875, 14170348.4375, 14170384.375, 14171021.875, 14171310.9375, 14171335.9375, 14171637.5, 14171764.0625, 14171826.5625, 14171915.625, 14172003.125, 14172129.6875, 14172282.8125, 14172326.5625, 14172656.25, 14172676.5625, 14172878.125, 14173100.0, 14173310.9375, 14173646.875, 14173726.5625, 14173798.4375, 14174464.0625, 14174809.375, 14174962.5, 14174967.1875, 14174976.5625, 14175090.625, 14175114.0625, 14175153.125, 14175573.4375, 14175637.5, 14175790.625, 14176089.0625, 14176718.75, 14176718.75, 14176867.1875, 14177201.5625, 14177900.0, 14178412.5, 14178857.8125, 14178859.375, 14178993.75, 14179395.3125, 14179567.1875, 14180114.0625, 14180900.0, 14181170.3125, 14182025.0, 14182746.875, 14183140.625, 14184557.8125, 14184934.375, 14184971.875, 14184990.625, 14185576.5625, 14185637.5, 14185712.5, 14185879.6875, 14185942.1875, 14186414.0625, 14186792.1875, 14187131.25, 14187528.125, 14188085.9375, 14188307.8125, 14188537.5, 14188575.0, 14188901.5625, 14189168.75, 14189495.3125, 14189845.3125, 14190823.4375, 14191256.25, 14191310.9375, 14191317.1875, 14192559.375, 14192592.1875, 14193084.375, 14193559.375, 14194429.6875, 14196209.375, 14197817.1875, 14199010.9375, 14199012.5, 14199107.8125, 14199748.4375, 14199756.25, 14200703.125, 14201914.0625, 14202307.8125, 14202560.9375, 14202989.0625, 14203189.0625, 14203348.4375, 14203595.3125, 14203756.25, 14204378.125, 14205239.0625, 14205481.25, 14205678.125, 14205978.125, 14206098.4375, 14206164.0625, 14206626.5625, 14207867.1875, 14210195.3125, 14210546.875, 14211262.5, 14212501.5625, 14214840.625, 14215814.0625, 14215979.6875, 14215981.25, 14216007.8125, 14216457.8125, 14216664.0625, 14221175.0, 14221579.6875, 14223371.875, 14224518.75, 14224739.0625, 14224815.625, 14225807.8125, 14225965.625, 14226720.3125, 14230601.5625, 14232651.5625, 14233571.875, 14234320.3125, 14234446.875, 14240440.625, 14240776.5625, 14243726.5625, 14245967.1875, 14248412.5, 14249862.5, 14250431.25, 14256106.25, 14262431.25, 14270289.0625, 14271584.375, 14282850.0, 14283479.6875, 14288748.4375, 14288756.25, 14316635.9375, 442749226.5625], [14.353494663391347, 473.72552766437803, 64.83849779080148, 36.77691193062272, 45.5215314765082, 9.35358762473517, 138.96490170140314, 10.736181369094693, 63.82109119366702, 55.938173334457915, 88.7681590638431, 19.734325713764772, 35.651572755161176, 46.550626326252385, 20.871126056680342, 21.84359957165835, 47.23218286247899, 6.13684904713674, 9.24638261305465, 30.44021871282233, 47.03775059633008, 6.55245419575673, 58.702711176285426, 23.636303341479437, 18.179892111916878, 89.94716101189168, 5.1297197217189465, 12.358679632605428, 10.71630147318415, 93.0766903499036, 17.52061275270044, 24.991509180526812, 20.46416628164172, 20.60424365329352, 9.121317134698764, 33.91656963226593, 19.029721329047504, 51.49528148779245, 88.70332418188231, 82.62434546059671, 5.2298468191344565, 58.33207553908103, 41.93331390532528, 86.784007584458, 51.30684960822346, 8.61703840863088, 32.63254148345757, 43.17868874281369, 12.753193533769785, 34.34814232164087, 10.167794305449949, 87.83604021181014, 39.235160788821744, 18.255997217943104, 6.01336053078085, 11.0482600586928, 32.090075499029574, 68.50731007923501, 39.51983998604173, 29.065530742111154, 25.36276690784278, 5.921981987503392, 84.37475442262945, 11.488771362701042, 50.02737375445998, 22.722740087727537, 17.169690837376297, 91.45494429636368, 21.335876551298693, 6.5764150127191465, 111.03219718051129, 85.299902567554, 23.945056982174265, 38.52536228253169, 28.972469800237363, 66.00821662089689, 13.676826791699021, 72.9907585653588, 95.70378321655397, 62.12037854049055, 66.7891484470813, 54.47259835167837, 27.28921849345727, 39.35013358710248, 28.2602832444484, 65.7435275050917, 12.801872069724844, 25.61020355370054, 16.63461569573696, 46.28079114633559, 96.29792409196753, 89.07687421935228, 44.22610227911099, 14.941244148153501, 13.057208369360833, 12.69087761508895, 38.26527543577155, 107.25872188508636, 35.72856526079983, 58.60820293278006, 7.372982563358372, 26.732792541294373, 68.17684566975235, 5.354712784245878, 64.41885900111, 66.85545750295871, 28.996827169854154, 112.19906756178189, 88.49173997056705, 28.492001388069177, 14.434911830894555, 14.534011182788714, 80.34683665108685, 29.485995613491788, 10.977925767294726, 105.69469775692211, 17.168300381557057, 20.036805185614934, 6.328084087484737, 39.78822185923844, 22.85319254143005, 44.848943485916244, 51.88871040089554, 20.871193667039336, 26.467093006754947, 15.410974375179872, 89.7989222877052, 5.1492059434781, 5.387654290123565, 16.536214208219764, 26.62462523961121, 33.20454150504682, 10.211351191310994, 19.601035348003556, 27.953888009990216, 140.85742704841465, 87.5724959221242, 7.7884728836066905, 17.888564589548565, 9.29390012359141, 27.901629567809184, 61.82385221825618, 69.82649026251953, 8.898215595327333, 95.92870616255469, 35.27538308689147, 18.52547826745096, 19.67910539562954, 77.94228142861903, 5.71437746038673, 44.47637852710898, 56.36949227037983, 29.022241895616872, 19.432695983014664, 26.510717661558907, 25.756896138374913, 11.326472613546093, 65.36063868588136, 50.949582862337884, 19.770498431721553, 13.538392312171128, 173.75082529786062, 42.85044921383103, 14.150076505973443, 113.0498348797251, 28.173464299292945, 7.611891421922517, 5.09488032341785, 35.222522822861414, 19.14699249790841, 58.64371926551527, 60.664350018351925, 23.22640493206027, 6.795742014943427, 68.45664886767102, 6.519365833077661, 81.72007403771516, 26.181263879217074, 17.219910229944606, 32.86166634729164, 44.27113433493194, 8.566632772941917, 7.517428094158713, 60.515464745670485, 72.78434080480221, 28.199652866668668, 8.960780437403438, 21.635228861273784, 14.910618693956629, 29.148818182271953, 5.216270756522859, 81.2912865049599, 46.614024831305414, 67.18256762792173, 15.867553571281848, 10.270541321095285, 73.82180026059304, 20.058651924123808, 37.071630126810064, 58.4809364584972, 22.550319185008476, 6.933272729737776, 25.587246155338363, 73.06943091049588, 97.59715330110649, 5.712836102823359, 83.06397033333556, 47.11838201349126, 25.303326137019152, 40.30559279508621, 11.340039130667057, 29.153939809226653, 64.22302227576327, 36.7537455405884, 17.5776176592042, 8.118588397117676, 29.570597285957536, 25.396005680513376, 18.822971165666523, 36.453855472548426, 18.710383379088043, 57.55298437809341, 29.590219543437573, 48.267791310782805, 46.37818819078264, 25.73942578247345, 29.593005254134496, 58.61179312681341, 5.744956007708263, 78.00161546103493, 11.199893235750558, 22.72205252054634, 97.31172513930277, 32.675650253293114, 80.89593158843272, 50.725112666314686, 21.639503016235178, 5.234156784139545, 21.2604218195829, 31.40272912710506, 38.89152686193783, 23.576606320267587, 24.12015277532069, 80.56142884208839, 61.783822976485546, 130.96431079441888, 7.225273460899372, 25.27921487972804, 9.823766914230644, 12.560017545618122, 58.076908098202786, 41.12700090295604, 76.07454717282067, 150.85943632080748, 27.65022827389516, 5.57519222573279, 11.76680217485271, 5.640250335941951, 38.65796054701666, 47.792779484209134, 5.357885297869613, 26.029438017087415, 29.604191167731127, 42.795952554645254, 68.58387361761261, 5.901564481523938, 11.45765452471494, 10.330741219667953, 115.23435989695996, 25.440892274317907, 43.92058086375303, 14.68165371458611, 77.07678788115845, 72.06944552585097, 10.3784304246614, 20.516357260276195, 29.66975523850248, 37.390459818246526, 45.10135817742194, 7.822877740003752, 6.9047549546724545, 69.02149245526022, 115.91437765085564, 18.29596293877664, 91.04978680844967, 19.987977949081962, 13.58785447880324, 32.631274458799595, 122.34863721889519, 8.687979854997394, 67.43830196743917, 19.796190578387257, 9.70874289866778, 5.3208183280715105, 26.166139992579883, 52.352664071284345, 65.78788479758539, 32.96582134133165, 22.31088865752372, 53.04019828907692, 7.6305856592876875, 68.29311692352044, 5.198107194523418, 14.925240601590707, 70.30190856292792, 15.257973668040526, 24.523792994991883, 12.128948368450294, 61.880366493780315, 64.35540451634562, 56.12343846276438, 6.036973747837919, 27.986491478131114, 9.635269287039392, 15.384076907805271, 6.379008749669283, 33.151920247776644, 16.08108154619088, 45.126872121236985, 14.26252348097054, 21.745799117344784, 66.19139995038633, 10.092529100972401, 32.880980227759096, 59.18544360724054, 62.37415011395015, 23.700612120865156, 39.22478817111433, 12.128612534526674, 46.7952996734885, 17.240567391596084, 19.563062491699153, 32.77943073174624, 75.9244532855812, 6.2623924026738536, 27.785666873615423, 6.0800116185134, 24.544525550550887, 190.73086764807385, 68.7782447723781, 6.76145788477542, 88.4272289875377, 32.31481695488803, 28.392484217973628, 9.289049871627414, 35.26804737796784, 162.9244149185386, 20.013158130557656, 12.28036558966326, 181.68393035370133, 104.97754878116041, 10.516428442332703, 21.81040952229024, 16.018448455714275, 18.646876630704586, 11.558945118318242, 34.64340098544084, 43.12197024156084, 26.711302433437393, 57.47291746726988, 9.971635292758434, 13.454504760957123, 42.061852668240164, 65.88109421558121, 5.839771119107743, 58.50996494348714, 14.806594037032847, 9.5681734380959, 79.75198020536143, 34.78258862659328, 12.586082027440696, 27.273693806676206, 23.35680763294256, 11.284664052981473, 57.65542306537539, 100.14971778596654, 31.78167166919605, 45.64889588083866, 38.65420264274514, 79.66758193074546, 17.1201419980156, 102.9097270235377, 17.219748229071147, 5.368158579478964, 31.601807003765117, 71.76241143826363, 32.939004211533664, 42.155056073147755, 11.987017610004475, 9.748965266877068, 19.767290591809402, 7.110682854338585, 7.392816442907737, 24.204429172715763, 91.49244085784125, 8.293983059766695, 23.937622027379852, 14.20480637365656, 9.68367346001704, 85.73824785275004, 60.35682557482924, 133.2052298197641, 11.943339722888041, 72.69601261208908, 120.63466956028064, 91.80284339144525, 9.815860612181233, 6.559279010428072, 33.564597357816474, 9.555392088807762, 78.3831574207153, 27.179852850244465, 6.594196462714507, 40.75154635054125, 36.724422283613656, 19.247476621369373, 37.929071005299356, 17.703290882445362, 107.20111832155399, 120.53668625294208, 16.841635354393926, 6.1544185316873605, 135.75325572677988, 74.31783458236261, 47.199725363220104, 49.58537836937867, 57.278543193439255, 53.959106161268366, 77.60096407933997, 58.56103256615294, 144.74222144461513, 85.70036375992598, 176.41859715798293, 11.115712889705602, 54.354030798643635, 30.350359721749584, 6.293897324225599, 55.71879324120669, 29.430641341484563, 107.56090829397996, 16.845050039553357, 40.235935353644614, 22.69808557703132, 76.87226004600889, 15.579405560730098, 10.010674334203065, 19.290120177564052, 38.52634537582876, 6.997355117862078, 13.638581184715129, 48.785429134526986, 42.00743014501606, 101.61683499211131, 73.20579185488225, 27.091800585534948, 30.619262391062065, 42.29559070599094, 6.137766473903482, 10.385163457992718, 43.24332129557821, 8.060733604365055, 50.43517550739337, 18.16350169807365, 9.86721091156452, 11.406037867396963, 17.72729173742572, 103.40576469931088, 75.34640461928947, 5.7576106600359225, 30.240818891386503, 37.46407154749019, 11.442331310233499, 79.65797703991664, 51.40732769005136, 65.57899699966401, 20.5501338259211, 29.623393177108042, 25.338588503117435, 128.6027701241257, 58.325084278925424, 47.76298990475915, 17.595070019309876, 9.867098645014597, 132.67222661305547, 8.921745451046558, 50.526320092458434, 51.58389783260908, 32.022575343000476, 31.963147853396553, 71.49294496650889, 105.55434234269504, 112.65856204740324, 29.25502346049976, 171.86869265024822, 58.31970632752492, 49.932414663227775, 80.57132976094407, 65.61759374109833, 23.22642695072207, 166.34354258819423, 27.542601311789923, 25.772105636994027, 63.93691962545487, 66.90935078348784, 106.14916392544066, 10.939203343101395, 42.124368198549114, 10.844508684040633, 7.241840816275167, 6.73736242483435, 11.593779428282327, 18.335779657925514, 21.41947902885403, 18.753579045500953, 19.35315697538517, 107.13513831793823, 31.550913024453926, 102.6628717816799, 35.88619177362325, 17.235990196897486, 24.141560673491583, 22.806446988252, 62.00804244895558, 21.8257731318311, 8.383265822904729, 11.875105911471099, 6.636418077090405, 5.961896993654193, 12.40819180158402, 72.92229378544869, 16.950758129527294, 12.532837565835399, 21.40140000993174, 5.978316730841685, 163.52628580424195, 18.29264744716334, 9.53987319215665, 23.17433371053665, 27.759873069797564, 25.77991529317198, 35.19983605801254, 38.35084678260363, 10.69775915491015, 15.950459283085388, 16.93716050207172, 13.2905645636859, 37.018675301060775, 27.618674979992512, 57.673680193818264, 60.50167468273238, 23.631718542799774, 57.29710130082808, 26.36064094057911, 15.111439999683835, 32.86672071402299, 7.4691393743382815, 60.48189238985312, 79.35631343185365, 24.277194971053582, 34.47297681181391, 42.79701611969088, 53.08722263706735, 8.30078276625573, 42.66492192149025, 72.4143974575048, 23.492589951121623, 24.503793887336872, 24.5654317567233, 73.56338891035281, 78.6167044099696, 22.832434501843665, 45.346258978838904, 27.835536129408673, 5.908811296577237, 9.81245982135465, 26.736256819297026, 19.553482018962708, 32.033837241875744, 22.219346359091624, 17.374203878533606, 9.454173248858073, 16.905984015553795, 14.604699809172235, 5.892264206598048, 24.750274386990252, 10.129997786648358, 115.03410106358012, 120.59809824278094, 28.445128655074825, 24.804282183490045, 45.60001329839841, 54.379135472663855, 100.36656776423267, 7.1822472978541665, 69.40832478912773, 23.052443759982285, 27.523937880857577, 8.04043585859457, 33.09133066645042, 36.03935778629557, 24.983223811319892, 29.60588377718645, 11.675196360863536, 5.603972321785733, 85.86093185281321, 6.330764863829363, 12.86462784245834, 10.095862764191946, 54.41931073050113, 26.569248425808297, 16.057670851516743, 5.553189605060328, 35.91835265641494, 22.943831240828146, 21.979195849849486, 57.64208553402391, 39.86349284853823, 33.64598782450551, 153.69751159374601, 8.179707912007789, 65.19907640435628, 6.3731746240121545, 8.203519165503396, 18.461192792801985, 9.378688151242594, 31.601214988486028, 17.527842294205264, 15.502400387266192, 15.759471304195365, 192.94054873959666, 81.66392771040061, 31.28089465363296, 7.897973252990759, 76.80174422391055, 12.998234122218102, 27.228943722390547, 10.255578470981368, 8.344127005302429, 31.227121693565774, 112.64753732933204, 84.57942169022655, 6.755584123921043, 129.6651810321321, 35.06480684977999, 10.725893355668381, 119.5733993740958, 15.437543296927409, 5.422740319555255, 5.111246949472533, 32.18588515908512, 40.77242655841336, 90.25513639741051, 7.280682639352037, 77.62422027098206, 22.99121109848209, 67.72342756292294, 16.276133475377677, 32.102843183515155, 10.832739398097864, 66.0874783140498, 57.8059826646064, 19.403919895286432, 90.47193379697455, 9.934088548358066, 5.297742092654492, 74.81000662410375, 14.028357983947314, 6.503026641515817, 38.110746131507874, 31.19567044194315, 28.85526189513476, 46.68178040088819, 25.370980019490304, 29.74618510053752, 5.812437935303385, 114.34868919062774, 28.228757623061927, 41.81577428548081, 57.84271723124805, 22.92408061043053, 27.165065377009284, 40.352581437828285, 56.09387447922988, 50.64836482406433, 7.555244432039499, 6.579163328083184, 17.63436454018527, 49.06731276152528, 45.18078872034652, 22.259135349107297, 16.029458217998748, 51.05830677514964, 11.939508329800633, 6.61063322611364, 11.814231171268109, 10.742024049354702, 13.755892036566769, 14.541561349340286, 42.66995775086826, 11.323124212455566, 68.19026545077588, 24.057057603169547, 107.2010066940721, 29.38064581138545, 5.238439750192551, 15.38065479939424, 26.428753066981713, 100.92082851566539, 18.95662968989847, 106.32004525234811, 32.54977997824817, 49.40870064583637, 13.029481349876233, 10.678344216231519, 45.35456948327778, 65.48779644159757, 24.123998075217237, 69.4666661381233, 5.727448883501773, 19.325350682528594, 72.44816343599858, 42.5988272568806, 73.89845024180056, 51.98895013992934, 9.491066001418263, 46.387717662298286, 7.222415955323188, 34.48454880254624, 65.68751152076723, 5.211427768320859, 27.440242615085037, 39.39707298326284, 9.6765657827956, 14.331049459827428, 13.246368053407936, 15.686035953566362, 20.941839106548784, 10.720418571356598, 13.845165038436733, 7.113718941942257, 63.05466688401277, 69.59607448694868, 29.281196176750544, 8.891762549465556, 5.364152630269463, 131.94940676017814, 103.58620754774593, 51.84897823387635, 13.318800262476152, 8.052412111664024, 13.150209584174325, 62.09884185938603, 9.312112045567218, 53.60266778314205, 15.492859181962173, 6.720785975304419, 15.336933400850091, 100.80972994565191, 40.125024789834015, 97.35213209262277, 79.17677888834599, 15.434523461552672, 72.41510840814257, 76.1290006632271, 25.95245345017412, 25.946580354324322, 67.39171313697904, 84.67664208775805, 7.883125964175795, 25.701161924387222, 15.0305663848643, 65.61402769362405, 18.714518515639014, 7.884252759880994, 66.12187057903228, 7.7422409944461945, 7.930287508880684, 33.63132475150989, 83.28621539989015, 28.71326309822679, 13.555599003069942, 11.198195616775747, 28.83893870711918, 71.99012084775593, 34.27059600174316, 15.491059102396095, 94.85796046303663, 60.80621641633502, 9.298607510593195, 24.418031367950462, 43.07865403523679, 58.993156938069575, 23.980590210919814, 22.608508374697323, 18.85112564381806, 76.78964369657797, 55.06861024881866, 13.293523150838999, 63.095268398131324, 56.17066112459775, 35.573440325706656, 35.01348108909672, 9.488463612785868, 79.62241535886976, 13.769130029172004, 25.56633873415291, 46.014920023969054, 37.90023427118021, 47.74529842289515, 14.821632874934174, 81.4225354428202, 65.16212021741647, 9.3581247125086, 130.8053047235203, 33.23626145194534, 67.33242632089512, 14.458149277836041, 25.679156622015196, 86.3978634807841, 75.85960454048146, 14.125211046946713, 77.32937189092672, 21.577605834222336, 17.301057635941863, 50.93057878173799, 7.845349388155738, 12.252503029805252, 5.8840402398967235, 15.690595495616433, 60.627401180799644, 7.976440232600355, 7.357624827285209, 52.26075744299292, 7.634449117748097, 8.493106532176006, 12.378104995582737])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)