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 = 44439
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);
([3705604.6875, 6423890.625, 6500756.25, 6513976.5625, 6794501.5625, 7090943.75, 7102685.9375, 7202932.8125, 7302295.3125, 7313960.9375, 7326523.4375, 7361387.5, 7535378.125, 7646403.125, 7696426.5625, 7709309.375, 7712450.0, 7715162.5, 7730945.3125, 7791773.4375, 7837214.0625, 7841965.625, 7842064.0625, 7848854.6875, 7870062.5, 7882487.5, 7934584.375, 7934942.1875, 7935275.0, 7968071.875, 7974945.3125, 7990709.375, 7998029.6875, 7999793.75, 8000640.625, 8001064.0625, 8002039.0625, 8003228.125, 8003254.6875, 8004153.125, 8004223.4375, 8004232.8125, 8032606.25, 8038804.6875, 8044009.375, 8067870.3125, 8069576.5625, 8070392.1875, 8072015.625, 8072076.5625, 8072134.375, 8079217.1875, 8099459.375, 8099460.9375, 8100382.8125, 8100540.625, 8117115.625, 8135281.25, 8137257.8125, 8138187.5, 8142634.375, 8162804.6875, 8167034.375, 8171234.375, 8176718.75, 8176734.375, 8177679.6875, 8179432.8125, 8184225.0, 8196209.375, 8197962.5, 8198471.875, 8199068.75, 8199381.25, 8199450.0, 8200035.9375, 8201050.0, 8201942.1875, 8202917.1875, 8205037.5, 8222787.5, 8223082.8125, 8228965.625, 8230635.9375, 8231132.8125, 8231281.25, 8234160.9375, 8244215.625, 8244221.875, 8245184.375, 8245387.5, 8245595.3125, 8247450.0, 8248010.9375, 8250584.375, 8250634.375, 8253532.8125, 8253695.3125, 8253762.5, 8253915.625, 8254303.125, 8254896.875, 8259273.4375, 8260720.3125, 8263546.875, 8272223.4375, 8279895.3125, 8281201.5625, 8281942.1875, 8285032.8125, 8291051.5625, 8296437.5, 8303520.3125, 8312259.375, 8312553.125, 8313365.625, 8314039.0625, 8315084.375, 8320967.1875, 8321356.25, 8330812.5, 8344367.1875, 8352362.5, 8353245.3125, 8353407.8125, 8353850.0, 8354371.875, 8354785.9375, 8355100.0, 8355523.4375, 8355531.25, 8355810.9375, 8355817.1875, 8356234.375, 8356581.25, 8357290.625, 8357590.625, 8358098.4375, 8359235.9375, 8359728.125, 8359735.9375, 8359810.9375, 8360296.875, 8368237.5, 8373996.875, 8379579.6875, 8403715.625, 8411807.8125, 8411870.3125, 8411882.8125, 8414796.875, 8417959.375, 8428298.4375, 8437976.5625, 8439992.1875, 8440532.8125, 8441332.8125, 8441337.5, 8445678.125, 8446182.8125, 8446273.4375, 8446281.25, 8447825.0, 8448601.5625, 8449657.8125, 8450968.75, 8452254.6875, 8452484.375, 8452984.375, 8453623.4375, 8464078.125, 8465335.9375, 8466023.4375, 8466051.5625, 8466254.6875, 8466278.125, 8467796.875, 8467920.3125, 8468101.5625, 8468282.8125, 8468625.0, 8469062.5, 8469207.8125, 8469251.5625, 8469889.0625, 8470368.75, 8470989.0625, 8471123.4375, 8471471.875, 8473667.1875, 8478482.8125, 8499490.625, 8504246.875, 8507578.125, 8508534.375, 8509046.875, 8509196.875, 8509568.75, 8509812.5, 8510259.375, 8510823.4375, 8512421.875, 8513431.25, 8519854.6875, 8531000.0, 8532071.875, 8534100.0, 8534706.25, 8534857.8125, 8535060.9375, 8535234.375, 8535242.1875, 8535328.125, 8536035.9375, 8536109.375, 8537093.75, 8537965.625, 8540003.125, 8541484.375, 8542129.6875, 8543498.4375, 8553581.25, 8562779.6875, 8565059.375, 8567484.375, 8570851.5625, 8573003.125, 8573556.25, 8574856.25, 8575293.75, 8575860.9375, 8575989.0625, 8576050.0, 8576101.5625, 8577500.0, 8578662.5, 8579892.1875, 8580604.6875, 8582126.5625, 8582167.1875, 8587206.25, 8588489.0625, 8590443.75, 8610356.25, 8638378.125, 8639889.0625, 8644953.125, 8645210.9375, 8645900.0, 8645900.0, 8646054.6875, 8646954.6875, 8647520.3125, 8648146.875, 8648215.625, 8648367.1875, 8648393.75, 8649575.0, 8650079.6875, 8651812.5, 8653095.3125, 8653823.4375, 8654046.875, 8655251.5625, 8655796.875, 8658664.0625, 8668089.0625, 8678429.6875, 8678451.5625, 8678456.25, 8679614.0625, 8680557.8125, 8682000.0, 8682310.9375, 8682328.125, 8682381.25, 8683881.25, 8687270.3125, 8687559.375, 8688534.375, 8689162.5, 8689168.75, 8693081.25, 8693425.0, 8704051.5625, 8707806.25, 8709673.4375, 8714892.1875, 8715429.6875, 8715456.25, 8716992.1875, 8719417.1875, 8719742.1875, 8723953.125, 8724154.6875, 8725073.4375, 8726704.6875, 8727390.625, 8727929.6875, 8736792.1875, 8738900.0, 8739376.5625, 8740020.3125, 8740739.0625, 8740756.25, 8759081.25, 8759553.125, 8763004.6875, 8764178.125, 8764562.5, 8764637.5, 8764709.375, 8764740.625, 8764851.5625, 8764945.3125, 8765257.8125, 8765273.4375, 8765832.8125, 8766326.5625, 8766956.25, 8769465.625, 8784559.375, 8785243.75, 8786615.625, 8788962.5, 8789932.8125, 8791257.8125, 8792246.875, 8792800.0, 8792884.375, 8793276.5625, 8798778.125, 8800532.8125, 8800568.75, 8801943.75, 8803034.375, 8803264.0625, 8803975.0, 8805159.375, 8818435.9375, 8822832.8125, 8825456.25, 8835157.8125, 8835281.25, 8840689.0625, 8851539.0625, 8852985.9375, 8854314.0625, 8858096.875, 8858700.0, 8858707.8125, 8858987.5, 8859803.125, 8870996.875, 8910817.1875, 8914404.6875, 8919392.1875, 8921954.6875, 8922271.875, 8922279.6875, 8925782.8125, 8926048.4375, 8926601.5625, 8926934.375, 8927957.8125, 8928184.375, 8928443.75, 8929007.8125, 8929214.0625, 8929362.5, 8929915.625, 8930304.6875, 8930553.125, 8930953.125, 8931014.0625, 8931046.875, 8931276.5625, 8931398.4375, 8931682.8125, 8932106.25, 8932226.5625, 8932337.5, 8932810.9375, 8932820.3125, 8934115.625, 8934584.375, 8935342.1875, 8935382.8125, 8935726.5625, 8935753.125, 8936329.6875, 8936339.0625, 8936540.625, 8936837.5, 8937060.9375, 8937307.8125, 8937628.125, 8938035.9375, 8938592.1875, 8939223.4375, 8939276.5625, 8947126.5625, 8948243.75, 8948312.5, 8948328.125, 8951134.375, 8952315.625, 8952662.5, 8954704.6875, 8955128.125, 8955151.5625, 8957035.9375, 8959267.1875, 8961725.0, 8962639.0625, 8963487.5, 8971059.375, 8972078.125, 8973175.0, 8975200.0, 8975264.0625, 8975703.125, 8977512.5, 8977548.4375, 8977918.75, 8985459.375, 8985473.4375, 8986131.25, 8986459.375, 8986596.875, 8986629.6875, 8987832.8125, 8988731.25, 8988835.9375, 8991071.875, 8992198.4375, 8992881.25, 8995573.4375, 9004615.625, 9030792.1875, 9030960.9375, 9031095.3125, 9031282.8125, 9032473.4375, 9036792.1875, 9036856.25, 9039098.4375, 9039409.375, 9039762.5, 9039989.0625, 9040539.0625, 9040915.625, 9040976.5625, 9041345.3125, 9041409.375, 9043701.5625, 9044026.5625, 9044132.8125, 9046395.3125, 9046598.4375, 9046812.5, 9046817.1875, 9046912.5, 9047639.0625, 9049565.625, 9049846.875, 9050868.75, 9051934.375, 9053914.0625, 9054840.625, 9056235.9375, 9058725.0, 9061704.6875, 9069498.4375, 9069506.25, 9071943.75, 9074067.1875, 9075142.1875, 9076225.0, 9079643.75, 9081614.0625, 9081675.0, 9081676.5625, 9082278.125, 9082650.0, 9083378.125, 9084345.3125, 9084979.6875, 9085076.5625, 9085618.75, 9086126.5625, 9086673.4375, 9086779.6875, 9087821.875, 9088042.1875, 9088301.5625, 9088359.375, 9088392.1875, 9088423.4375, 9088679.6875, 9088706.25, 9089032.8125, 9089089.0625, 9089090.625, 9089173.4375, 9089217.1875, 9089218.75, 9089262.5, 9089264.0625, 9089270.3125, 9089279.6875, 9089284.375, 9089296.875, 9089362.5, 9089504.6875, 9089829.6875, 9090100.0, 9090242.1875, 9090256.25, 9090464.0625, 9090473.4375, 9090717.1875, 9090735.9375, 9091303.125, 9091325.0, 9091785.9375, 9092210.9375, 9092214.0625, 9092325.0, 9092954.6875, 9093354.6875, 9093379.6875, 9093604.6875, 9093629.6875, 9093687.5, 9093778.125, 9093931.25, 9094551.5625, 9094818.75, 9094921.875, 9095182.8125, 9095276.5625, 9095471.875, 9095534.375, 9095685.9375, 9095809.375, 9096037.5, 9096096.875, 9096195.3125, 9096242.1875, 9096242.1875, 9096648.4375, 9096678.125, 9096814.0625, 9096831.25, 9096970.3125, 9097176.5625, 9097214.0625, 9097537.5, 9097664.0625, 9097728.125, 9097748.4375, 9098431.25, 9098443.75, 9098864.0625, 9098931.25, 9099078.125, 9099678.125, 9099709.375, 9099714.0625, 9105495.3125, 9108160.9375, 9119932.8125, 9122787.5, 9124354.6875, 9124364.0625, 9124396.875, 9124523.4375, 9124770.3125, 9124803.125, 9125487.5, 9125600.0, 9125693.75, 9126878.125, 9128679.6875, 9128707.8125, 9128709.375, 9129215.625, 9129317.1875, 9130765.625, 9131331.25, 9131907.8125, 9132014.0625, 9132043.75, 9133329.6875, 9134289.0625, 9134609.375, 9134832.8125, 9144650.0, 9153195.3125, 9153267.1875, 9154398.4375, 9156373.4375, 9157345.3125, 9160698.4375, 9160732.8125, 9161307.8125, 9163648.4375, 9164257.8125, 9164760.9375, 9165667.1875, 9165689.0625, 9166060.9375, 9166140.625, 9166407.8125, 9166967.1875, 9169387.5, 9169910.9375, 9170482.8125, 9170926.5625, 9171160.9375, 9172639.0625, 9172945.3125, 9173367.1875, 9173400.0, 9173781.25, 9175565.625, 9175992.1875, 9176331.25, 9179114.0625, 9184820.3125, 9184832.8125, 9185929.6875, 9188104.6875, 9188690.625, 9189268.75, 9189432.8125, 9189615.625, 9190246.875, 9190251.5625, 9190262.5, 9190501.5625, 9191315.625, 9191387.5, 9191881.25, 9192189.0625, 9192607.8125, 9192660.9375, 9194310.9375, 9198184.375, 9198498.4375, 9202181.25, 9205237.5, 9206037.5, 9210162.5, 9210389.0625, 9216592.1875, 9217500.0, 9221870.3125, 9223168.75, 9240759.375, 9241048.4375, 9241375.0, 9241795.3125, 9242664.0625, 9243109.375, 9243143.75, 9244673.4375, 9244793.75, 9244818.75, 9245759.375, 9248735.9375, 9249584.375, 9249587.5, 9250393.75, 9251046.875, 9251048.4375, 9251504.6875, 9252232.8125, 9252876.5625, 9252981.25, 9254890.625, 9255895.3125, 9265015.625, 9266023.4375, 9279628.125, 9287020.3125, 9287171.875, 9290442.1875, 9290993.75, 9291664.0625, 9293692.1875, 9293726.5625, 9294456.25, 9294507.8125, 9294514.0625, 9298482.8125, 9298946.875, 9299182.8125, 9299525.0, 9300209.375, 9301351.5625, 9302185.9375, 9303482.8125, 9306239.0625, 9306371.875, 9309070.3125, 9317134.375, 9348784.375, 9352528.125, 9383178.125, 9387189.0625, 9388006.25, 9388368.75, 9388903.125, 9389192.1875, 9389398.4375, 9390550.0, 9398114.0625, 9399178.125, 9400410.9375, 9407659.375, 9408317.1875, 9413340.625, 9414045.3125, 9414671.875, 9415359.375, 9420593.75, 9451454.6875, 9454281.25, 9472517.1875, 9475528.125, 9490748.4375, 9491090.625, 9491407.8125, 9491739.0625, 9492835.9375, 9493862.5, 9494696.875, 9519956.25, 9520500.0, 9520517.1875, 9520606.25, 9520757.8125, 9522117.1875, 9522321.875, 9525976.5625, 9526635.9375, 9527129.6875, 9527682.8125, 9529068.75, 9537310.9375, 9547309.375, 9549134.375, 9551140.625, 9555470.3125, 9556670.3125, 9577035.9375, 9627831.25, 9633131.25, 9633567.1875, 9634851.5625, 9634942.1875, 9634996.875, 9636200.0, 9636203.125, 9674834.375, 9686006.25, 9765542.1875, 9765554.6875, 9765781.25, 9765817.1875, 9766620.3125, 9768617.1875, 9770923.4375, 9771557.8125, 9775157.8125, 9859306.25, 10024181.25, 10025139.0625, 10030987.5, 10253185.9375, 10254517.1875, 10254604.6875, 10256017.1875, 10256078.125, 10256331.25, 10256784.375, 10256810.9375, 10256853.125, 10259625.0, 10273998.4375, 10294403.125, 10294453.125, 10296456.25, 10317260.9375, 10325451.5625, 10333953.125, 10420651.5625, 10480359.375, 10553034.375, 10554129.6875, 10555826.5625, 10580101.5625, 10601931.25, 10602167.1875, 10607145.3125, 10673568.75, 10689646.875, 10734523.4375, 10742792.1875, 11032814.0625, 11034987.5, 11067242.1875, 11070414.0625, 11073853.125, 11076217.1875, 11093540.625, 11128373.4375], [32.924657097325294, 116.6798834281576, 81.00691152143972, 13.81606582078742, 41.72555211612796, 27.067201900344223, 9.661675720852736, 82.78728982993094, 14.823666307008384, 50.47189156931825, 73.32101371470245, 16.13505519609505, 14.418534181308267, 20.83983263389171, 5.6965500646658445, 8.894401020402597, 58.60041078339818, 11.856012650062224, 57.13348803974097, 61.19822632708832, 6.018440701337877, 15.15612161580969, 5.039071960146389, 28.72195793469428, 23.10672267943091, 97.58608264704937, 75.35854863886624, 9.659575726354321, 5.685506942610424, 20.335945054826034, 58.75481869670872, 93.3671116648696, 17.897319648338065, 58.727334773427245, 53.03302349377744, 80.39872580958976, 100.51427166015104, 80.13425660560868, 21.60709772063372, 69.38740188222866, 74.5877110724207, 12.78891806195076, 9.145712155769688, 115.848675321599, 61.38400223940031, 17.525835980327578, 6.607019620649585, 37.88480676717241, 17.77627576050788, 6.198439810677378, 71.16212715782888, 13.925571693954117, 14.53270219022515, 56.47847311257261, 14.231431015597083, 5.191186985993551, 61.42556379724776, 16.260243983812476, 6.241256862918884, 7.906379731635614, 8.549370819691367, 64.65193001122084, 30.189325033217266, 5.758292771471893, 15.115693467049823, 74.04248132686806, 55.20489449149838, 5.212654762827162, 6.571187188756432, 89.76614159046558, 10.465239039799483, 5.029637402935534, 53.66821510682772, 24.985309575968472, 20.83597270999763, 60.34876815587375, 16.30591490556375, 42.063002244972274, 75.90222048543963, 10.728107875943651, 8.947732010090657, 5.908348804966057, 87.09635711889646, 8.024669039121248, 62.78872949930695, 18.035970787139966, 36.6274253304104, 89.57669413445683, 64.25627901150489, 125.8307875237152, 12.03887369343171, 60.13972059924853, 12.465507224304247, 125.15352480115966, 88.92581524380074, 18.23269855016324, 97.77919819247317, 164.12279828797415, 5.500703622250263, 51.410707353117, 12.308662345715513, 5.551288558591052, 7.729951864461873, 16.779813778676992, 114.01548200973278, 6.5225869783925114, 30.922506918927596, 62.74326095199561, 23.335546081325806, 74.26539906834553, 40.56839744776385, 47.17016024195928, 5.87290474759791, 20.994805317749, 27.642249403970858, 13.73910176767776, 21.1696653546308, 9.488463220362247, 50.98497923456754, 16.613098978222425, 68.44332311093979, 50.416608986945874, 46.39505885148927, 80.97850407183961, 5.856558114020477, 5.226520871294633, 27.797452185396622, 59.124549090495805, 25.49913342130373, 17.582139815441657, 12.727219598469084, 211.41610012978518, 17.15872374385065, 40.74508636717269, 23.595943065991705, 34.332486437982745, 12.870203049543697, 28.44004240532065, 46.15426647738029, 57.41312617612931, 74.0944622169626, 108.3079182877103, 128.15072494588716, 22.762698457694675, 59.56549031715537, 77.67081983330007, 71.28092089875014, 109.1983634950865, 5.087629448471466, 5.02612513931732, 172.80589021724447, 78.51851359915393, 74.45847793068981, 59.61830626984404, 42.12325779692252, 80.06486479088169, 59.05093839089593, 43.18937772131305, 68.07164278428874, 35.42035311134909, 17.079249324906534, 7.61054281595935, 62.004903548381776, 84.10756466593455, 115.36061851669324, 52.9216767045551, 96.82843774883669, 20.884818239873937, 7.668040580038916, 7.11332217615482, 14.344377386407078, 38.745995673451446, 5.48425647205571, 18.60014119719927, 24.366283839026533, 143.428561693925, 48.59463022082683, 11.138488383048657, 11.01047232136609, 7.041889125072349, 5.1902418741117105, 31.636465234104293, 66.70689770171029, 98.21874010969626, 10.942707948514915, 87.98031492065158, 78.20688103665884, 75.80383253753945, 56.17196284850187, 45.62371635788635, 74.53985454164697, 100.36954080425524, 14.67279338980616, 21.47146830699488, 11.602682034196775, 32.61574729027931, 87.73579343352041, 14.670771706702002, 20.43916203688921, 6.006154518156957, 21.416024844655446, 20.244868802468304, 6.70066476499348, 75.10459147983843, 7.1079328292833415, 33.78828705240565, 68.33043595502669, 7.088111775935066, 5.2592237910222615, 41.957385930903285, 16.951986633866152, 32.30809348890219, 72.6679734717716, 61.608331974283146, 66.31096207653296, 43.56543120750834, 35.0035442328987, 17.938597863334167, 10.3130314100056, 97.39706789239654, 22.851046801484483, 12.129908465551022, 37.49006667836005, 16.3947930434081, 15.934215058800195, 151.14226456567974, 12.914239952692064, 16.337265247714317, 11.844890387450622, 6.944662971350165, 6.08410190696702, 25.40423052074012, 43.533659711533936, 8.805757463629188, 80.8808191305475, 27.18252032850513, 125.39988663636795, 24.19919465587248, 5.887121302681618, 165.90263302406132, 34.508508498647906, 20.561151933922357, 36.86235548799108, 90.92139993504784, 50.05667601548069, 25.487998106876354, 141.48940988255507, 7.243799439184941, 6.587818403883618, 121.46000486333148, 55.15767294459172, 21.486340336232416, 25.48974125085641, 23.477542479703597, 14.320472620269658, 21.497624785132977, 15.566253335702699, 44.19846895731496, 5.167972977641077, 8.081175676856867, 126.0309991132097, 14.978635612315523, 16.28236732110962, 7.583503821865062, 82.17689781359329, 41.66946614246198, 14.4587787266855, 100.85049292882078, 30.768355277168403, 22.96530113646394, 16.35121506637152, 39.62119117255041, 8.525017105057323, 60.40707931082138, 10.87866388263598, 27.565435484921743, 24.571254282781453, 21.588018015068272, 128.2668210825618, 27.624151243146777, 114.22089447212235, 11.800058439264006, 16.53051006309701, 30.692140727846795, 41.678786439275264, 76.49965600994615, 30.84402499776383, 45.286576270859285, 20.842401000776654, 49.414366725282974, 11.70703496129263, 13.022546636047572, 11.838428320057156, 27.90432821271198, 47.3357336198286, 9.649760613378662, 33.6895291039139, 42.50795531437848, 11.675433502415846, 173.7987759447381, 104.9719710209032, 40.1259550382697, 10.10100506682189, 26.551873548215895, 8.996573066915651, 32.41341909816228, 33.071746696685985, 38.46286975728948, 28.650557694252885, 79.23470504918554, 121.83991317282812, 45.60153871703153, 24.44779696699885, 39.434282884618405, 74.24008279882129, 19.20779174963493, 51.38379060714112, 64.4015378613884, 43.19510912778543, 59.084060065371446, 38.81403666317835, 18.83579110389391, 49.96896346212341, 55.92919405936321, 57.788837123061114, 49.59490519928747, 96.29733125882511, 55.169546460358234, 12.437623525684874, 15.726839210117292, 98.98227815791182, 9.82112699442779, 40.228218061072845, 46.81303675477902, 110.09165585815043, 100.2597246318238, 31.03538933586202, 15.317744582333278, 37.43139119613992, 59.52038818983135, 48.22109401350293, 32.717217153901444, 25.865032447816056, 13.208018311005814, 5.9003606704665055, 27.379948933004183, 14.479266667070963, 75.57768193779717, 18.589234915344335, 15.722342295291924, 32.85661532047301, 99.89236781071969, 37.424265835090324, 92.30584422484525, 92.05991887316418, 15.61320573818369, 8.63095863393564, 8.760347537697715, 33.476510617181674, 18.827789038991863, 26.888165929433242, 30.42175663256922, 16.069880451532686, 5.262989843630826, 40.03971070819298, 14.134046835277053, 17.020048795183868, 50.50387825877301, 151.22889611269392, 63.46524483034014, 16.14903298780898, 5.667255884978681, 9.629352873382768, 44.689850219089266, 17.669481666425575, 55.332010720627196, 13.73043771923356, 19.473972904421682, 41.86750954546676, 66.08167162159023, 10.620321132711709, 76.37485386571336, 65.41713063108222, 11.298951269996039, 41.64482405005603, 84.65469006238222, 81.02903411316058, 27.52793745658734, 42.474287791121014, 12.715564225568837, 15.652578630488666, 37.14309259239211, 27.71625870172194, 15.152537691637114, 23.096684496603803, 39.686267491266236, 34.016309444265524, 14.873539009205832, 18.445964371100416, 189.78456389753975, 39.92180615912983, 8.261552645615602, 15.588276082733929, 60.55771075958826, 29.076066288314937, 108.79209425868126, 20.710337521655884, 47.36582564176112, 54.92282267726806, 25.410236199615987, 30.261546145513798, 24.772628264234005, 88.90794795481821, 11.64244701970497, 192.50614707089505, 28.675964727111598, 13.378745903290769, 12.875199809248988, 70.8817727454168, 14.232152521302606, 67.66386061235895, 131.5717429906109, 10.723698552120933, 85.4214145050205, 27.11142764702241, 29.30161052447525, 6.792385579804992, 14.465028927842702, 86.2864948315138, 26.886914076650385, 29.165727144989788, 36.382368188891824, 123.3290671831657, 35.47817831868368, 47.65733211960614, 64.52371765098714, 9.119815875592273, 77.32845423359748, 19.839212661948206, 130.86225844963167, 16.41030622824398, 44.63710634732309, 33.40569040456431, 58.37489605488815, 27.812340684483356, 5.024411214800042, 41.80996621495609, 19.28292084063712, 14.097641579230096, 40.59118815383231, 42.9961004748717, 7.6652682804176635, 68.15295699729796, 104.18122702871335, 70.46224013378512, 68.64033638212265, 14.388731608272009, 172.00391572644023, 50.801326739344915, 22.300684320065606, 30.302611871915524, 20.740584225982126, 15.984587492122163, 65.79103634635933, 61.28101592649081, 8.969244101743984, 56.69863098235196, 36.39553605323646, 19.938333112624754, 40.852170448050344, 31.234427155665674, 5.8190971116115, 182.8420777150173, 11.111837386432626, 63.616456844031106, 56.61250372517559, 11.663793021300107, 10.089796506861942, 123.35194978335716, 79.62773183223466, 17.83898689457889, 7.998607243290227, 43.07773779308858, 67.78000143328207, 43.154529947723304, 26.75180031867474, 52.47549176457697, 40.472853778637734, 96.37876472137663, 45.39333692737364, 22.713574577667814, 58.752260654513776, 21.406971700910937, 50.74241546664142, 61.12786118575792, 27.19324336266716, 23.699518848411273, 116.91549173960323, 66.37872713013465, 32.71818471615478, 34.809660207489486, 10.637679965990744, 74.71396575972813, 13.976734496700999, 132.71265320995658, 30.94502928732678, 11.661705287161318, 6.033893429643892, 24.661209020386643, 26.283522644497964, 89.35037805921596, 35.47240006511114, 50.42968566431969, 36.208889858415404, 10.019641602200524, 65.26607823632742, 24.94490476700947, 24.42131428353687, 10.208898612059109, 15.088051621424992, 36.67242457886821, 47.77566024903621, 36.07453992963525, 50.71645487604973, 7.5307774819693245, 39.04070543560128, 64.52297220800261, 26.81748508265912, 40.16312478517369, 15.454442383661968, 115.37027191545279, 134.6126938270227, 64.0158972814034, 8.261571179151076, 13.56233045772898, 118.5113676453322, 99.12068873961482, 116.65802870113792, 60.79419175324312, 8.85095042406844, 33.0001338894392, 29.848726291123196, 80.978850633698, 37.42312878248814, 41.49936258411806, 6.208724910462915, 9.038345578817697, 5.378216745630152, 12.599176265682507, 9.576960779908998, 37.9768887904705, 20.57193689307942, 7.47054296424512, 28.972645932052327, 51.946150292117785, 16.455350704646413, 6.622425127852205, 121.03723275264161, 74.04103502061801, 115.62040635983799, 55.15835805526446, 18.153447245211975, 26.568445284829195, 17.804512556224864, 54.87089187103154, 79.43299725586613, 73.00508114607749, 8.05672486198217, 68.66837907266587, 131.2828158163505, 39.748842547621585, 26.00035424728346, 36.53023358958911, 45.923752593590265, 65.54674609903553, 16.936684315190387, 151.90447416488368, 19.232006465704693, 88.39092742095102, 28.570599775977055, 96.83386437787198, 68.94188119201912, 39.2016901752217, 8.950548092181926, 158.03791325754122, 71.23842625613656, 37.035100116729225, 78.22005215695403, 5.103273678674295, 5.352198469536739, 14.744225916779666, 8.905076504710934, 76.36010568699902, 43.02168149222666, 27.268192235020923, 30.474329458927773, 18.826413672841, 25.436450950583133, 28.37893010368078, 79.82336219502548, 5.190963802034682, 6.830287051783128, 13.836966409664848, 51.437551372011285, 57.42260720854609, 7.360613115738645, 22.89453082897128, 23.85191904698461, 97.29391762377679, 21.49236196980726, 70.16292413629075, 40.08199733977131, 63.94249628045913, 29.131857014537644, 25.198277011720318, 39.598212776084644, 62.81178636082504, 25.823389893061904, 66.04055978992595, 100.76757322412215, 33.63814292770292, 12.537869524972136, 83.65746169930932, 111.68595237956785, 19.039552547928842, 262.2673539374696, 90.5290369595815, 25.205566693605142, 22.51922230158157, 25.228548720088604, 90.22342466673786, 31.618407060578242, 101.39848925022315, 58.4626720075059, 57.05113638378787, 73.68164603567162, 175.76699501964458, 30.60912336745204, 112.53885813476577, 52.01525922659552, 20.742798833424285, 19.4055430654669, 99.59278737399248, 7.1539607964270076, 57.52491950535498, 62.91495774737473, 68.55294170451401, 5.065912136159971, 11.70429541438887, 44.03004171602744, 54.01662279964817, 7.6301049148870925, 24.156903634909927, 92.74577462411881, 93.96752856229226, 125.42717094001031, 89.06086857117145, 14.092339971336317, 8.366229186720703, 29.86870122150375, 81.16301792941056, 26.270465451492335, 22.588621813012146, 15.0348306778534, 12.135946702055199, 7.362079397250393, 65.48394931434993, 8.351663089339961, 24.475895806110724, 112.68030759982324, 24.218606158559993, 15.468194403853039, 6.183517318434525, 22.38633697814639, 24.00400358928337, 61.040960658143604, 17.95748969298244, 21.95179224030898, 12.931370348004553, 12.74110154457865, 48.93929917469359, 227.67706996602868, 112.51998434093974, 41.582941685496486, 96.74943666517736, 67.81621497977122, 50.690291204949254, 57.622419757448895, 47.072283883215206, 65.72980877496182, 26.610746086471913, 21.397502792242683, 37.52089472206346, 25.20553236446617, 16.676651704488503, 27.99152455666368, 60.36288306624489, 13.01485261568247, 5.61294970333641, 14.577076388617769, 9.580380276606364, 113.32975505283501, 52.87055021224667, 126.0882688787489, 89.92300040829763, 46.63568302929067, 38.79100329453242, 19.98976518116541, 6.888931911550446, 9.979114224176064, 101.72348327452988, 5.5675429672881185, 72.26115133243339, 43.41732545177426, 108.14565597231343, 19.329126093552816, 21.28846306528019, 12.220944754374417, 13.078644239520068, 137.24892675503747, 8.109864531215313, 70.02273076282208, 60.26695669772552, 118.82262364836316, 24.724064203774613, 5.555967513944521, 28.328050168583378, 93.28746462161298, 52.38185693816583, 149.3918214312369, 46.401212688278754, 141.4068641054947, 245.12840350386492, 32.47302767718489, 18.900967615763058, 24.08344403997618, 6.110924514752332, 63.57555515550892, 61.110101710307035, 32.34003870050973, 119.41474780228708, 19.919834974370858, 101.66238449718038, 8.440986096295617, 25.686160201138723, 25.757665396587274, 31.33793909630443, 30.22003949943895, 57.872915494078455, 59.295995251688666, 23.80224850446178, 9.473432820587302, 53.53458669037609, 210.0824999246817, 19.351889690919258, 72.09739342401426, 19.192115014275775, 6.184879174974784, 86.94955722549035, 16.64295476996734, 84.18898657715172, 108.72345167435368, 13.616465639399365, 5.252983170080755, 21.79300048815033, 81.54846902119787, 26.33760758170059, 207.0198146956072, 12.186649544505443, 8.141526326198806, 9.70631881222407, 23.72574446910234, 57.8872899125006, 25.86148399056148, 29.51829299553497, 67.48981852624637, 19.41387128134355, 25.769864105974953, 33.74541501654854, 15.632634922910334, 44.73308627702366, 14.159321979927359, 37.398354260490876, 104.86826328920962, 16.53196829735078, 69.97091518926179, 19.801065962390993, 75.03072688972566, 16.9822742979065, 6.5496748682950825, 5.4626251441583005, 14.577347472508464, 17.221316960507394, 37.1947214574466, 89.31249921380503, 26.03248583382198, 10.199416811522697, 49.295336627075514, 22.663647664048113, 90.0617620885011, 13.385534696093494, 15.71123891577353, 5.6769395033954675, 69.79158888822076, 99.11028403624456, 126.74564138490084, 35.62904257849863, 8.716734065067033, 103.79829180561171, 23.318639176450983, 9.968603535685736, 55.64024346618522, 76.21506179656281, 17.400712255662746, 72.86711976442741, 7.787611774742526, 17.554050112901276, 53.93162568684195, 48.86700727739273])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3705604.6875, 6423890.625, 6500756.25, 6513976.5625, 6794501.5625, 7090943.75, 7102685.9375, 7202932.8125, 7302295.3125, 7313960.9375, 7326523.4375, 7361387.5, 7535378.125, 7646403.125, 7696426.5625, 7709309.375, 7712450.0, 7715162.5, 7730945.3125, 7791773.4375, 7837214.0625, 7841965.625, 7842064.0625, 7848854.6875, 7870062.5, 7882487.5, 7934584.375, 7934942.1875, 7935275.0, 7968071.875, 7974945.3125, 7990709.375, 7998029.6875, 7999793.75, 8000640.625, 8001064.0625, 8002039.0625, 8003228.125, 8003254.6875, 8004153.125, 8004223.4375, 8004232.8125, 8032606.25, 8038804.6875, 8044009.375, 8067870.3125, 8069576.5625, 8070392.1875, 8072015.625, 8072076.5625, 8072134.375, 8079217.1875, 8099459.375, 8099460.9375, 8100382.8125, 8100540.625, 8117115.625, 8135281.25, 8137257.8125, 8138187.5, 8142634.375, 8162804.6875, 8167034.375, 8171234.375, 8176718.75, 8176734.375, 8177679.6875, 8179432.8125, 8184225.0, 8196209.375, 8197962.5, 8198471.875, 8199068.75, 8199381.25, 8199450.0, 8200035.9375, 8201050.0, 8201942.1875, 8202917.1875, 8205037.5, 8222787.5, 8223082.8125, 8228965.625, 8230635.9375, 8231132.8125, 8231281.25, 8234160.9375, 8244215.625, 8244221.875, 8245184.375, 8245387.5, 8245595.3125, 8247450.0, 8248010.9375, 8250584.375, 8250634.375, 8253532.8125, 8253695.3125, 8253762.5, 8253915.625, 8254303.125, 8254896.875, 8259273.4375, 8260720.3125, 8263546.875, 8272223.4375, 8279895.3125, 8281201.5625, 8281942.1875, 8285032.8125, 8291051.5625, 8296437.5, 8303520.3125, 8312259.375, 8312553.125, 8313365.625, 8314039.0625, 8315084.375, 8320967.1875, 8321356.25, 8330812.5, 8344367.1875, 8352362.5, 8353245.3125, 8353407.8125, 8353850.0, 8354371.875, 8354785.9375, 8355100.0, 8355523.4375, 8355531.25, 8355810.9375, 8355817.1875, 8356234.375, 8356581.25, 8357290.625, 8357590.625, 8358098.4375, 8359235.9375, 8359728.125, 8359735.9375, 8359810.9375, 8360296.875, 8368237.5, 8373996.875, 8379579.6875, 8403715.625, 8411807.8125, 8411870.3125, 8411882.8125, 8414796.875, 8417959.375, 8428298.4375, 8437976.5625, 8439992.1875, 8440532.8125, 8441332.8125, 8441337.5, 8445678.125, 8446182.8125, 8446273.4375, 8446281.25, 8447825.0, 8448601.5625, 8449657.8125, 8450968.75, 8452254.6875, 8452484.375, 8452984.375, 8453623.4375, 8464078.125, 8465335.9375, 8466023.4375, 8466051.5625, 8466254.6875, 8466278.125, 8467796.875, 8467920.3125, 8468101.5625, 8468282.8125, 8468625.0, 8469062.5, 8469207.8125, 8469251.5625, 8469889.0625, 8470368.75, 8470989.0625, 8471123.4375, 8471471.875, 8473667.1875, 8478482.8125, 8499490.625, 8504246.875, 8507578.125, 8508534.375, 8509046.875, 8509196.875, 8509568.75, 8509812.5, 8510259.375, 8510823.4375, 8512421.875, 8513431.25, 8519854.6875, 8531000.0, 8532071.875, 8534100.0, 8534706.25, 8534857.8125, 8535060.9375, 8535234.375, 8535242.1875, 8535328.125, 8536035.9375, 8536109.375, 8537093.75, 8537965.625, 8540003.125, 8541484.375, 8542129.6875, 8543498.4375, 8553581.25, 8562779.6875, 8565059.375, 8567484.375, 8570851.5625, 8573003.125, 8573556.25, 8574856.25, 8575293.75, 8575860.9375, 8575989.0625, 8576050.0, 8576101.5625, 8577500.0, 8578662.5, 8579892.1875, 8580604.6875, 8582126.5625, 8582167.1875, 8587206.25, 8588489.0625, 8590443.75, 8610356.25, 8638378.125, 8639889.0625, 8644953.125, 8645210.9375, 8645900.0, 8645900.0, 8646054.6875, 8646954.6875, 8647520.3125, 8648146.875, 8648215.625, 8648367.1875, 8648393.75, 8649575.0, 8650079.6875, 8651812.5, 8653095.3125, 8653823.4375, 8654046.875, 8655251.5625, 8655796.875, 8658664.0625, 8668089.0625, 8678429.6875, 8678451.5625, 8678456.25, 8679614.0625, 8680557.8125, 8682000.0, 8682310.9375, 8682328.125, 8682381.25, 8683881.25, 8687270.3125, 8687559.375, 8688534.375, 8689162.5, 8689168.75, 8693081.25, 8693425.0, 8704051.5625, 8707806.25, 8709673.4375, 8714892.1875, 8715429.6875, 8715456.25, 8716992.1875, 8719417.1875, 8719742.1875, 8723953.125, 8724154.6875, 8725073.4375, 8726704.6875, 8727390.625, 8727929.6875, 8736792.1875, 8738900.0, 8739376.5625, 8740020.3125, 8740739.0625, 8740756.25, 8759081.25, 8759553.125, 8763004.6875, 8764178.125, 8764562.5, 8764637.5, 8764709.375, 8764740.625, 8764851.5625, 8764945.3125, 8765257.8125, 8765273.4375, 8765832.8125, 8766326.5625, 8766956.25, 8769465.625, 8784559.375, 8785243.75, 8786615.625, 8788962.5, 8789932.8125, 8791257.8125, 8792246.875, 8792800.0, 8792884.375, 8793276.5625, 8798778.125, 8800532.8125, 8800568.75, 8801943.75, 8803034.375, 8803264.0625, 8803975.0, 8805159.375, 8818435.9375, 8822832.8125, 8825456.25, 8835157.8125, 8835281.25, 8840689.0625, 8851539.0625, 8852985.9375, 8854314.0625, 8858096.875, 8858700.0, 8858707.8125, 8858987.5, 8859803.125, 8870996.875, 8910817.1875, 8914404.6875, 8919392.1875, 8921954.6875, 8922271.875, 8922279.6875, 8925782.8125, 8926048.4375, 8926601.5625, 8926934.375, 8927957.8125, 8928184.375, 8928443.75, 8929007.8125, 8929214.0625, 8929362.5, 8929915.625, 8930304.6875, 8930553.125, 8930953.125, 8931014.0625, 8931046.875, 8931276.5625, 8931398.4375, 8931682.8125, 8932106.25, 8932226.5625, 8932337.5, 8932810.9375, 8932820.3125, 8934115.625, 8934584.375, 8935342.1875, 8935382.8125, 8935726.5625, 8935753.125, 8936329.6875, 8936339.0625, 8936540.625, 8936837.5, 8937060.9375, 8937307.8125, 8937628.125, 8938035.9375, 8938592.1875, 8939223.4375, 8939276.5625, 8947126.5625, 8948243.75, 8948312.5, 8948328.125, 8951134.375, 8952315.625, 8952662.5, 8954704.6875, 8955128.125, 8955151.5625, 8957035.9375, 8959267.1875, 8961725.0, 8962639.0625, 8963487.5, 8971059.375, 8972078.125, 8973175.0, 8975200.0, 8975264.0625, 8975703.125, 8977512.5, 8977548.4375, 8977918.75, 8985459.375, 8985473.4375, 8986131.25, 8986459.375, 8986596.875, 8986629.6875, 8987832.8125, 8988731.25, 8988835.9375, 8991071.875, 8992198.4375, 8992881.25, 8995573.4375, 9004615.625, 9030792.1875, 9030960.9375, 9031095.3125, 9031282.8125, 9032473.4375, 9036792.1875, 9036856.25, 9039098.4375, 9039409.375, 9039762.5, 9039989.0625, 9040539.0625, 9040915.625, 9040976.5625, 9041345.3125, 9041409.375, 9043701.5625, 9044026.5625, 9044132.8125, 9046395.3125, 9046598.4375, 9046812.5, 9046817.1875, 9046912.5, 9047639.0625, 9049565.625, 9049846.875, 9050868.75, 9051934.375, 9053914.0625, 9054840.625, 9056235.9375, 9058725.0, 9061704.6875, 9069498.4375, 9069506.25, 9071943.75, 9074067.1875, 9075142.1875, 9076225.0, 9079643.75, 9081614.0625, 9081675.0, 9081676.5625, 9082278.125, 9082650.0, 9083378.125, 9084345.3125, 9084979.6875, 9085076.5625, 9085618.75, 9086126.5625, 9086673.4375, 9086779.6875, 9087821.875, 9088042.1875, 9088301.5625, 9088359.375, 9088392.1875, 9088423.4375, 9088679.6875, 9088706.25, 9089032.8125, 9089089.0625, 9089090.625, 9089173.4375, 9089217.1875, 9089218.75, 9089262.5, 9089264.0625, 9089270.3125, 9089279.6875, 9089284.375, 9089296.875, 9089362.5, 9089504.6875, 9089829.6875, 9090100.0, 9090242.1875, 9090256.25, 9090464.0625, 9090473.4375, 9090717.1875, 9090735.9375, 9091303.125, 9091325.0, 9091785.9375, 9092210.9375, 9092214.0625, 9092325.0, 9092954.6875, 9093354.6875, 9093379.6875, 9093604.6875, 9093629.6875, 9093687.5, 9093778.125, 9093931.25, 9094551.5625, 9094818.75, 9094921.875, 9095182.8125, 9095276.5625, 9095471.875, 9095534.375, 9095685.9375, 9095809.375, 9096037.5, 9096096.875, 9096195.3125, 9096242.1875, 9096242.1875, 9096648.4375, 9096678.125, 9096814.0625, 9096831.25, 9096970.3125, 9097176.5625, 9097214.0625, 9097537.5, 9097664.0625, 9097728.125, 9097748.4375, 9098431.25, 9098443.75, 9098864.0625, 9098931.25, 9099078.125, 9099678.125, 9099709.375, 9099714.0625, 9105495.3125, 9108160.9375, 9119932.8125, 9122787.5, 9124354.6875, 9124364.0625, 9124396.875, 9124523.4375, 9124770.3125, 9124803.125, 9125487.5, 9125600.0, 9125693.75, 9126878.125, 9128679.6875, 9128707.8125, 9128709.375, 9129215.625, 9129317.1875, 9130765.625, 9131331.25, 9131907.8125, 9132014.0625, 9132043.75, 9133329.6875, 9134289.0625, 9134609.375, 9134832.8125, 9144650.0, 9153195.3125, 9153267.1875, 9154398.4375, 9156373.4375, 9157345.3125, 9160698.4375, 9160732.8125, 9161307.8125, 9163648.4375, 9164257.8125, 9164760.9375, 9165667.1875, 9165689.0625, 9166060.9375, 9166140.625, 9166407.8125, 9166967.1875, 9169387.5, 9169910.9375, 9170482.8125, 9170926.5625, 9171160.9375, 9172639.0625, 9172945.3125, 9173367.1875, 9173400.0, 9173781.25, 9175565.625, 9175992.1875, 9176331.25, 9179114.0625, 9184820.3125, 9184832.8125, 9185929.6875, 9188104.6875, 9188690.625, 9189268.75, 9189432.8125, 9189615.625, 9190246.875, 9190251.5625, 9190262.5, 9190501.5625, 9191315.625, 9191387.5, 9191881.25, 9192189.0625, 9192607.8125, 9192660.9375, 9194310.9375, 9198184.375, 9198498.4375, 9202181.25, 9205237.5, 9206037.5, 9210162.5, 9210389.0625, 9216592.1875, 9217500.0, 9221870.3125, 9223168.75, 9240759.375, 9241048.4375, 9241375.0, 9241795.3125, 9242664.0625, 9243109.375, 9243143.75, 9244673.4375, 9244793.75, 9244818.75, 9245759.375, 9248735.9375, 9249584.375, 9249587.5, 9250393.75, 9251046.875, 9251048.4375, 9251504.6875, 9252232.8125, 9252876.5625, 9252981.25, 9254890.625, 9255895.3125, 9265015.625, 9266023.4375, 9279628.125, 9287020.3125, 9287171.875, 9290442.1875, 9290993.75, 9291664.0625, 9293692.1875, 9293726.5625, 9294456.25, 9294507.8125, 9294514.0625, 9298482.8125, 9298946.875, 9299182.8125, 9299525.0, 9300209.375, 9301351.5625, 9302185.9375, 9303482.8125, 9306239.0625, 9306371.875, 9309070.3125, 9317134.375, 9348784.375, 9352528.125, 9383178.125, 9387189.0625, 9388006.25, 9388368.75, 9388903.125, 9389192.1875, 9389398.4375, 9390550.0, 9398114.0625, 9399178.125, 9400410.9375, 9407659.375, 9408317.1875, 9413340.625, 9414045.3125, 9414671.875, 9415359.375, 9420593.75, 9451454.6875, 9454281.25, 9472517.1875, 9475528.125, 9490748.4375, 9491090.625, 9491407.8125, 9491739.0625, 9492835.9375, 9493862.5, 9494696.875, 9519956.25, 9520500.0, 9520517.1875, 9520606.25, 9520757.8125, 9522117.1875, 9522321.875, 9525976.5625, 9526635.9375, 9527129.6875, 9527682.8125, 9529068.75, 9537310.9375, 9547309.375, 9549134.375, 9551140.625, 9555470.3125, 9556670.3125, 9577035.9375, 9627831.25, 9633131.25, 9633567.1875, 9634851.5625, 9634942.1875, 9634996.875, 9636200.0, 9636203.125, 9674834.375, 9686006.25, 9765542.1875, 9765554.6875, 9765781.25, 9765817.1875, 9766620.3125, 9768617.1875, 9770923.4375, 9771557.8125, 9775157.8125, 9859306.25, 10024181.25, 10025139.0625, 10030987.5, 10253185.9375, 10254517.1875, 10254604.6875, 10256017.1875, 10256078.125, 10256331.25, 10256784.375, 10256810.9375, 10256853.125, 10259625.0, 10273998.4375, 10294403.125, 10294453.125, 10296456.25, 10317260.9375, 10325451.5625, 10333953.125, 10420651.5625, 10480359.375, 10553034.375, 10554129.6875, 10555826.5625, 10580101.5625, 10601931.25, 10602167.1875, 10607145.3125, 10673568.75, 10689646.875, 10734523.4375, 10742792.1875, 11032814.0625, 11034987.5, 11067242.1875, 11070414.0625, 11073853.125, 11076217.1875, 11093540.625, 11128373.4375], [32.924657097325294, 116.6798834281576, 81.00691152143972, 13.81606582078742, 41.72555211612796, 27.067201900344223, 9.661675720852736, 82.78728982993094, 14.823666307008384, 50.47189156931825, 73.32101371470245, 16.13505519609505, 14.418534181308267, 20.83983263389171, 5.6965500646658445, 8.894401020402597, 58.60041078339818, 11.856012650062224, 57.13348803974097, 61.19822632708832, 6.018440701337877, 15.15612161580969, 5.039071960146389, 28.72195793469428, 23.10672267943091, 97.58608264704937, 75.35854863886624, 9.659575726354321, 5.685506942610424, 20.335945054826034, 58.75481869670872, 93.3671116648696, 17.897319648338065, 58.727334773427245, 53.03302349377744, 80.39872580958976, 100.51427166015104, 80.13425660560868, 21.60709772063372, 69.38740188222866, 74.5877110724207, 12.78891806195076, 9.145712155769688, 115.848675321599, 61.38400223940031, 17.525835980327578, 6.607019620649585, 37.88480676717241, 17.77627576050788, 6.198439810677378, 71.16212715782888, 13.925571693954117, 14.53270219022515, 56.47847311257261, 14.231431015597083, 5.191186985993551, 61.42556379724776, 16.260243983812476, 6.241256862918884, 7.906379731635614, 8.549370819691367, 64.65193001122084, 30.189325033217266, 5.758292771471893, 15.115693467049823, 74.04248132686806, 55.20489449149838, 5.212654762827162, 6.571187188756432, 89.76614159046558, 10.465239039799483, 5.029637402935534, 53.66821510682772, 24.985309575968472, 20.83597270999763, 60.34876815587375, 16.30591490556375, 42.063002244972274, 75.90222048543963, 10.728107875943651, 8.947732010090657, 5.908348804966057, 87.09635711889646, 8.024669039121248, 62.78872949930695, 18.035970787139966, 36.6274253304104, 89.57669413445683, 64.25627901150489, 125.8307875237152, 12.03887369343171, 60.13972059924853, 12.465507224304247, 125.15352480115966, 88.92581524380074, 18.23269855016324, 97.77919819247317, 164.12279828797415, 5.500703622250263, 51.410707353117, 12.308662345715513, 5.551288558591052, 7.729951864461873, 16.779813778676992, 114.01548200973278, 6.5225869783925114, 30.922506918927596, 62.74326095199561, 23.335546081325806, 74.26539906834553, 40.56839744776385, 47.17016024195928, 5.87290474759791, 20.994805317749, 27.642249403970858, 13.73910176767776, 21.1696653546308, 9.488463220362247, 50.98497923456754, 16.613098978222425, 68.44332311093979, 50.416608986945874, 46.39505885148927, 80.97850407183961, 5.856558114020477, 5.226520871294633, 27.797452185396622, 59.124549090495805, 25.49913342130373, 17.582139815441657, 12.727219598469084, 211.41610012978518, 17.15872374385065, 40.74508636717269, 23.595943065991705, 34.332486437982745, 12.870203049543697, 28.44004240532065, 46.15426647738029, 57.41312617612931, 74.0944622169626, 108.3079182877103, 128.15072494588716, 22.762698457694675, 59.56549031715537, 77.67081983330007, 71.28092089875014, 109.1983634950865, 5.087629448471466, 5.02612513931732, 172.80589021724447, 78.51851359915393, 74.45847793068981, 59.61830626984404, 42.12325779692252, 80.06486479088169, 59.05093839089593, 43.18937772131305, 68.07164278428874, 35.42035311134909, 17.079249324906534, 7.61054281595935, 62.004903548381776, 84.10756466593455, 115.36061851669324, 52.9216767045551, 96.82843774883669, 20.884818239873937, 7.668040580038916, 7.11332217615482, 14.344377386407078, 38.745995673451446, 5.48425647205571, 18.60014119719927, 24.366283839026533, 143.428561693925, 48.59463022082683, 11.138488383048657, 11.01047232136609, 7.041889125072349, 5.1902418741117105, 31.636465234104293, 66.70689770171029, 98.21874010969626, 10.942707948514915, 87.98031492065158, 78.20688103665884, 75.80383253753945, 56.17196284850187, 45.62371635788635, 74.53985454164697, 100.36954080425524, 14.67279338980616, 21.47146830699488, 11.602682034196775, 32.61574729027931, 87.73579343352041, 14.670771706702002, 20.43916203688921, 6.006154518156957, 21.416024844655446, 20.244868802468304, 6.70066476499348, 75.10459147983843, 7.1079328292833415, 33.78828705240565, 68.33043595502669, 7.088111775935066, 5.2592237910222615, 41.957385930903285, 16.951986633866152, 32.30809348890219, 72.6679734717716, 61.608331974283146, 66.31096207653296, 43.56543120750834, 35.0035442328987, 17.938597863334167, 10.3130314100056, 97.39706789239654, 22.851046801484483, 12.129908465551022, 37.49006667836005, 16.3947930434081, 15.934215058800195, 151.14226456567974, 12.914239952692064, 16.337265247714317, 11.844890387450622, 6.944662971350165, 6.08410190696702, 25.40423052074012, 43.533659711533936, 8.805757463629188, 80.8808191305475, 27.18252032850513, 125.39988663636795, 24.19919465587248, 5.887121302681618, 165.90263302406132, 34.508508498647906, 20.561151933922357, 36.86235548799108, 90.92139993504784, 50.05667601548069, 25.487998106876354, 141.48940988255507, 7.243799439184941, 6.587818403883618, 121.46000486333148, 55.15767294459172, 21.486340336232416, 25.48974125085641, 23.477542479703597, 14.320472620269658, 21.497624785132977, 15.566253335702699, 44.19846895731496, 5.167972977641077, 8.081175676856867, 126.0309991132097, 14.978635612315523, 16.28236732110962, 7.583503821865062, 82.17689781359329, 41.66946614246198, 14.4587787266855, 100.85049292882078, 30.768355277168403, 22.96530113646394, 16.35121506637152, 39.62119117255041, 8.525017105057323, 60.40707931082138, 10.87866388263598, 27.565435484921743, 24.571254282781453, 21.588018015068272, 128.2668210825618, 27.624151243146777, 114.22089447212235, 11.800058439264006, 16.53051006309701, 30.692140727846795, 41.678786439275264, 76.49965600994615, 30.84402499776383, 45.286576270859285, 20.842401000776654, 49.414366725282974, 11.70703496129263, 13.022546636047572, 11.838428320057156, 27.90432821271198, 47.3357336198286, 9.649760613378662, 33.6895291039139, 42.50795531437848, 11.675433502415846, 173.7987759447381, 104.9719710209032, 40.1259550382697, 10.10100506682189, 26.551873548215895, 8.996573066915651, 32.41341909816228, 33.071746696685985, 38.46286975728948, 28.650557694252885, 79.23470504918554, 121.83991317282812, 45.60153871703153, 24.44779696699885, 39.434282884618405, 74.24008279882129, 19.20779174963493, 51.38379060714112, 64.4015378613884, 43.19510912778543, 59.084060065371446, 38.81403666317835, 18.83579110389391, 49.96896346212341, 55.92919405936321, 57.788837123061114, 49.59490519928747, 96.29733125882511, 55.169546460358234, 12.437623525684874, 15.726839210117292, 98.98227815791182, 9.82112699442779, 40.228218061072845, 46.81303675477902, 110.09165585815043, 100.2597246318238, 31.03538933586202, 15.317744582333278, 37.43139119613992, 59.52038818983135, 48.22109401350293, 32.717217153901444, 25.865032447816056, 13.208018311005814, 5.9003606704665055, 27.379948933004183, 14.479266667070963, 75.57768193779717, 18.589234915344335, 15.722342295291924, 32.85661532047301, 99.89236781071969, 37.424265835090324, 92.30584422484525, 92.05991887316418, 15.61320573818369, 8.63095863393564, 8.760347537697715, 33.476510617181674, 18.827789038991863, 26.888165929433242, 30.42175663256922, 16.069880451532686, 5.262989843630826, 40.03971070819298, 14.134046835277053, 17.020048795183868, 50.50387825877301, 151.22889611269392, 63.46524483034014, 16.14903298780898, 5.667255884978681, 9.629352873382768, 44.689850219089266, 17.669481666425575, 55.332010720627196, 13.73043771923356, 19.473972904421682, 41.86750954546676, 66.08167162159023, 10.620321132711709, 76.37485386571336, 65.41713063108222, 11.298951269996039, 41.64482405005603, 84.65469006238222, 81.02903411316058, 27.52793745658734, 42.474287791121014, 12.715564225568837, 15.652578630488666, 37.14309259239211, 27.71625870172194, 15.152537691637114, 23.096684496603803, 39.686267491266236, 34.016309444265524, 14.873539009205832, 18.445964371100416, 189.78456389753975, 39.92180615912983, 8.261552645615602, 15.588276082733929, 60.55771075958826, 29.076066288314937, 108.79209425868126, 20.710337521655884, 47.36582564176112, 54.92282267726806, 25.410236199615987, 30.261546145513798, 24.772628264234005, 88.90794795481821, 11.64244701970497, 192.50614707089505, 28.675964727111598, 13.378745903290769, 12.875199809248988, 70.8817727454168, 14.232152521302606, 67.66386061235895, 131.5717429906109, 10.723698552120933, 85.4214145050205, 27.11142764702241, 29.30161052447525, 6.792385579804992, 14.465028927842702, 86.2864948315138, 26.886914076650385, 29.165727144989788, 36.382368188891824, 123.3290671831657, 35.47817831868368, 47.65733211960614, 64.52371765098714, 9.119815875592273, 77.32845423359748, 19.839212661948206, 130.86225844963167, 16.41030622824398, 44.63710634732309, 33.40569040456431, 58.37489605488815, 27.812340684483356, 5.024411214800042, 41.80996621495609, 19.28292084063712, 14.097641579230096, 40.59118815383231, 42.9961004748717, 7.6652682804176635, 68.15295699729796, 104.18122702871335, 70.46224013378512, 68.64033638212265, 14.388731608272009, 172.00391572644023, 50.801326739344915, 22.300684320065606, 30.302611871915524, 20.740584225982126, 15.984587492122163, 65.79103634635933, 61.28101592649081, 8.969244101743984, 56.69863098235196, 36.39553605323646, 19.938333112624754, 40.852170448050344, 31.234427155665674, 5.8190971116115, 182.8420777150173, 11.111837386432626, 63.616456844031106, 56.61250372517559, 11.663793021300107, 10.089796506861942, 123.35194978335716, 79.62773183223466, 17.83898689457889, 7.998607243290227, 43.07773779308858, 67.78000143328207, 43.154529947723304, 26.75180031867474, 52.47549176457697, 40.472853778637734, 96.37876472137663, 45.39333692737364, 22.713574577667814, 58.752260654513776, 21.406971700910937, 50.74241546664142, 61.12786118575792, 27.19324336266716, 23.699518848411273, 116.91549173960323, 66.37872713013465, 32.71818471615478, 34.809660207489486, 10.637679965990744, 74.71396575972813, 13.976734496700999, 132.71265320995658, 30.94502928732678, 11.661705287161318, 6.033893429643892, 24.661209020386643, 26.283522644497964, 89.35037805921596, 35.47240006511114, 50.42968566431969, 36.208889858415404, 10.019641602200524, 65.26607823632742, 24.94490476700947, 24.42131428353687, 10.208898612059109, 15.088051621424992, 36.67242457886821, 47.77566024903621, 36.07453992963525, 50.71645487604973, 7.5307774819693245, 39.04070543560128, 64.52297220800261, 26.81748508265912, 40.16312478517369, 15.454442383661968, 115.37027191545279, 134.6126938270227, 64.0158972814034, 8.261571179151076, 13.56233045772898, 118.5113676453322, 99.12068873961482, 116.65802870113792, 60.79419175324312, 8.85095042406844, 33.0001338894392, 29.848726291123196, 80.978850633698, 37.42312878248814, 41.49936258411806, 6.208724910462915, 9.038345578817697, 5.378216745630152, 12.599176265682507, 9.576960779908998, 37.9768887904705, 20.57193689307942, 7.47054296424512, 28.972645932052327, 51.946150292117785, 16.455350704646413, 6.622425127852205, 121.03723275264161, 74.04103502061801, 115.62040635983799, 55.15835805526446, 18.153447245211975, 26.568445284829195, 17.804512556224864, 54.87089187103154, 79.43299725586613, 73.00508114607749, 8.05672486198217, 68.66837907266587, 131.2828158163505, 39.748842547621585, 26.00035424728346, 36.53023358958911, 45.923752593590265, 65.54674609903553, 16.936684315190387, 151.90447416488368, 19.232006465704693, 88.39092742095102, 28.570599775977055, 96.83386437787198, 68.94188119201912, 39.2016901752217, 8.950548092181926, 158.03791325754122, 71.23842625613656, 37.035100116729225, 78.22005215695403, 5.103273678674295, 5.352198469536739, 14.744225916779666, 8.905076504710934, 76.36010568699902, 43.02168149222666, 27.268192235020923, 30.474329458927773, 18.826413672841, 25.436450950583133, 28.37893010368078, 79.82336219502548, 5.190963802034682, 6.830287051783128, 13.836966409664848, 51.437551372011285, 57.42260720854609, 7.360613115738645, 22.89453082897128, 23.85191904698461, 97.29391762377679, 21.49236196980726, 70.16292413629075, 40.08199733977131, 63.94249628045913, 29.131857014537644, 25.198277011720318, 39.598212776084644, 62.81178636082504, 25.823389893061904, 66.04055978992595, 100.76757322412215, 33.63814292770292, 12.537869524972136, 83.65746169930932, 111.68595237956785, 19.039552547928842, 262.2673539374696, 90.5290369595815, 25.205566693605142, 22.51922230158157, 25.228548720088604, 90.22342466673786, 31.618407060578242, 101.39848925022315, 58.4626720075059, 57.05113638378787, 73.68164603567162, 175.76699501964458, 30.60912336745204, 112.53885813476577, 52.01525922659552, 20.742798833424285, 19.4055430654669, 99.59278737399248, 7.1539607964270076, 57.52491950535498, 62.91495774737473, 68.55294170451401, 5.065912136159971, 11.70429541438887, 44.03004171602744, 54.01662279964817, 7.6301049148870925, 24.156903634909927, 92.74577462411881, 93.96752856229226, 125.42717094001031, 89.06086857117145, 14.092339971336317, 8.366229186720703, 29.86870122150375, 81.16301792941056, 26.270465451492335, 22.588621813012146, 15.0348306778534, 12.135946702055199, 7.362079397250393, 65.48394931434993, 8.351663089339961, 24.475895806110724, 112.68030759982324, 24.218606158559993, 15.468194403853039, 6.183517318434525, 22.38633697814639, 24.00400358928337, 61.040960658143604, 17.95748969298244, 21.95179224030898, 12.931370348004553, 12.74110154457865, 48.93929917469359, 227.67706996602868, 112.51998434093974, 41.582941685496486, 96.74943666517736, 67.81621497977122, 50.690291204949254, 57.622419757448895, 47.072283883215206, 65.72980877496182, 26.610746086471913, 21.397502792242683, 37.52089472206346, 25.20553236446617, 16.676651704488503, 27.99152455666368, 60.36288306624489, 13.01485261568247, 5.61294970333641, 14.577076388617769, 9.580380276606364, 113.32975505283501, 52.87055021224667, 126.0882688787489, 89.92300040829763, 46.63568302929067, 38.79100329453242, 19.98976518116541, 6.888931911550446, 9.979114224176064, 101.72348327452988, 5.5675429672881185, 72.26115133243339, 43.41732545177426, 108.14565597231343, 19.329126093552816, 21.28846306528019, 12.220944754374417, 13.078644239520068, 137.24892675503747, 8.109864531215313, 70.02273076282208, 60.26695669772552, 118.82262364836316, 24.724064203774613, 5.555967513944521, 28.328050168583378, 93.28746462161298, 52.38185693816583, 149.3918214312369, 46.401212688278754, 141.4068641054947, 245.12840350386492, 32.47302767718489, 18.900967615763058, 24.08344403997618, 6.110924514752332, 63.57555515550892, 61.110101710307035, 32.34003870050973, 119.41474780228708, 19.919834974370858, 101.66238449718038, 8.440986096295617, 25.686160201138723, 25.757665396587274, 31.33793909630443, 30.22003949943895, 57.872915494078455, 59.295995251688666, 23.80224850446178, 9.473432820587302, 53.53458669037609, 210.0824999246817, 19.351889690919258, 72.09739342401426, 19.192115014275775, 6.184879174974784, 86.94955722549035, 16.64295476996734, 84.18898657715172, 108.72345167435368, 13.616465639399365, 5.252983170080755, 21.79300048815033, 81.54846902119787, 26.33760758170059, 207.0198146956072, 12.186649544505443, 8.141526326198806, 9.70631881222407, 23.72574446910234, 57.8872899125006, 25.86148399056148, 29.51829299553497, 67.48981852624637, 19.41387128134355, 25.769864105974953, 33.74541501654854, 15.632634922910334, 44.73308627702366, 14.159321979927359, 37.398354260490876, 104.86826328920962, 16.53196829735078, 69.97091518926179, 19.801065962390993, 75.03072688972566, 16.9822742979065, 6.5496748682950825, 5.4626251441583005, 14.577347472508464, 17.221316960507394, 37.1947214574466, 89.31249921380503, 26.03248583382198, 10.199416811522697, 49.295336627075514, 22.663647664048113, 90.0617620885011, 13.385534696093494, 15.71123891577353, 5.6769395033954675, 69.79158888822076, 99.11028403624456, 126.74564138490084, 35.62904257849863, 8.716734065067033, 103.79829180561171, 23.318639176450983, 9.968603535685736, 55.64024346618522, 76.21506179656281, 17.400712255662746, 72.86711976442741, 7.787611774742526, 17.554050112901276, 53.93162568684195, 48.86700727739273])
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);
([3705604.6875, 6423890.625, 6500756.25, 6513976.5625, 6794501.5625, 7090943.75, 7102685.9375, 7202932.8125, 7302295.3125, 7313960.9375, 7326523.4375, 7361387.5, 7535378.125, 7646403.125, 7696426.5625, 7709309.375, 7712450.0, 7715162.5, 7730945.3125, 7791773.4375, 7837214.0625, 7841965.625, 7842064.0625, 7848854.6875, 7870062.5, 7882487.5, 7934584.375, 7934942.1875, 7935275.0, 7968071.875, 7974945.3125, 7990709.375, 7998029.6875, 7999793.75, 8000640.625, 8001064.0625, 8002039.0625, 8003228.125, 8003254.6875, 8004153.125, 8004223.4375, 8004232.8125, 8032606.25, 8038804.6875, 8044009.375, 8067870.3125, 8069576.5625, 8070392.1875, 8072015.625, 8072076.5625, 8072134.375, 8079217.1875, 8099459.375, 8099460.9375, 8100382.8125, 8100540.625, 8117115.625, 8135281.25, 8137257.8125, 8138187.5, 8142634.375, 8162804.6875, 8167034.375, 8171234.375, 8176718.75, 8176734.375, 8177679.6875, 8179432.8125, 8184225.0, 8196209.375, 8197962.5, 8198471.875, 8199068.75, 8199381.25, 8199450.0, 8200035.9375, 8201050.0, 8201942.1875, 8202917.1875, 8205037.5, 8222787.5, 8223082.8125, 8228965.625, 8230635.9375, 8231132.8125, 8231281.25, 8234160.9375, 8244215.625, 8244221.875, 8245184.375, 8245387.5, 8245595.3125, 8247450.0, 8248010.9375, 8250584.375, 8250634.375, 8253532.8125, 8253695.3125, 8253762.5, 8253915.625, 8254303.125, 8254896.875, 8259273.4375, 8260720.3125, 8263546.875, 8272223.4375, 8279895.3125, 8281201.5625, 8281942.1875, 8285032.8125, 8291051.5625, 8296437.5, 8303520.3125, 8312259.375, 8312553.125, 8313365.625, 8314039.0625, 8315084.375, 8320967.1875, 8321356.25, 8330812.5, 8344367.1875, 8352362.5, 8353245.3125, 8353407.8125, 8353850.0, 8354371.875, 8354785.9375, 8355100.0, 8355523.4375, 8355531.25, 8355810.9375, 8355817.1875, 8356234.375, 8356581.25, 8357290.625, 8357590.625, 8358098.4375, 8359235.9375, 8359728.125, 8359735.9375, 8359810.9375, 8360296.875, 8368237.5, 8373996.875, 8379579.6875, 8403715.625, 8411807.8125, 8411870.3125, 8411882.8125, 8414796.875, 8417959.375, 8428298.4375, 8437976.5625, 8439992.1875, 8440532.8125, 8441332.8125, 8441337.5, 8445678.125, 8446182.8125, 8446273.4375, 8446281.25, 8447825.0, 8448601.5625, 8449657.8125, 8450968.75, 8452254.6875, 8452484.375, 8452984.375, 8453623.4375, 8464078.125, 8465335.9375, 8466023.4375, 8466051.5625, 8466254.6875, 8466278.125, 8467796.875, 8467920.3125, 8468101.5625, 8468282.8125, 8468625.0, 8469062.5, 8469207.8125, 8469251.5625, 8469889.0625, 8470368.75, 8470989.0625, 8471123.4375, 8471471.875, 8473667.1875, 8478482.8125, 8499490.625, 8504246.875, 8507578.125, 8508534.375, 8509046.875, 8509196.875, 8509568.75, 8509812.5, 8510259.375, 8510823.4375, 8512421.875, 8513431.25, 8519854.6875, 8531000.0, 8532071.875, 8534100.0, 8534706.25, 8534857.8125, 8535060.9375, 8535234.375, 8535242.1875, 8535328.125, 8536035.9375, 8536109.375, 8537093.75, 8537965.625, 8540003.125, 8541484.375, 8542129.6875, 8543498.4375, 8553581.25, 8562779.6875, 8565059.375, 8567484.375, 8570851.5625, 8573003.125, 8573556.25, 8574856.25, 8575293.75, 8575860.9375, 8575989.0625, 8576050.0, 8576101.5625, 8577500.0, 8578662.5, 8579892.1875, 8580604.6875, 8582126.5625, 8582167.1875, 8587206.25, 8588489.0625, 8590443.75, 8610356.25, 8638378.125, 8639889.0625, 8644953.125, 8645210.9375, 8645900.0, 8645900.0, 8646054.6875, 8646954.6875, 8647520.3125, 8648146.875, 8648215.625, 8648367.1875, 8648393.75, 8649575.0, 8650079.6875, 8651812.5, 8653095.3125, 8653823.4375, 8654046.875, 8655251.5625, 8655796.875, 8658664.0625, 8668089.0625, 8678429.6875, 8678451.5625, 8678456.25, 8679614.0625, 8680557.8125, 8682000.0, 8682310.9375, 8682328.125, 8682381.25, 8683881.25, 8687270.3125, 8687559.375, 8688534.375, 8689162.5, 8689168.75, 8693081.25, 8693425.0, 8704051.5625, 8707806.25, 8709673.4375, 8714892.1875, 8715429.6875, 8715456.25, 8716992.1875, 8719417.1875, 8719742.1875, 8723953.125, 8724154.6875, 8725073.4375, 8726704.6875, 8727390.625, 8727929.6875, 8736792.1875, 8738900.0, 8739376.5625, 8740020.3125, 8740739.0625, 8740756.25, 8759081.25, 8759553.125, 8763004.6875, 8764178.125, 8764562.5, 8764637.5, 8764709.375, 8764740.625, 8764851.5625, 8764945.3125, 8765257.8125, 8765273.4375, 8765832.8125, 8766326.5625, 8766956.25, 8769465.625, 8784559.375, 8785243.75, 8786615.625, 8788962.5, 8789932.8125, 8791257.8125, 8792246.875, 8792800.0, 8792884.375, 8793276.5625, 8798778.125, 8800532.8125, 8800568.75, 8801943.75, 8803034.375, 8803264.0625, 8803975.0, 8805159.375, 8818435.9375, 8822832.8125, 8825456.25, 8835157.8125, 8835281.25, 8840689.0625, 8851539.0625, 8852985.9375, 8854314.0625, 8858096.875, 8858700.0, 8858707.8125, 8858987.5, 8859803.125, 8870996.875, 8910817.1875, 8914404.6875, 8919392.1875, 8921954.6875, 8922271.875, 8922279.6875, 8925782.8125, 8926048.4375, 8926601.5625, 8926934.375, 8927957.8125, 8928184.375, 8928443.75, 8929007.8125, 8929214.0625, 8929362.5, 8929915.625, 8930304.6875, 8930553.125, 8930953.125, 8931014.0625, 8931046.875, 8931276.5625, 8931398.4375, 8931682.8125, 8932106.25, 8932226.5625, 8932337.5, 8932810.9375, 8932820.3125, 8934115.625, 8934584.375, 8935342.1875, 8935382.8125, 8935726.5625, 8935753.125, 8936329.6875, 8936339.0625, 8936540.625, 8936837.5, 8937060.9375, 8937307.8125, 8937628.125, 8938035.9375, 8938592.1875, 8939223.4375, 8939276.5625, 8947126.5625, 8948243.75, 8948312.5, 8948328.125, 8951134.375, 8952315.625, 8952662.5, 8954704.6875, 8955128.125, 8955151.5625, 8957035.9375, 8959267.1875, 8961725.0, 8962639.0625, 8963487.5, 8971059.375, 8972078.125, 8973175.0, 8975200.0, 8975264.0625, 8975703.125, 8977512.5, 8977548.4375, 8977918.75, 8985459.375, 8985473.4375, 8986131.25, 8986459.375, 8986596.875, 8986629.6875, 8987832.8125, 8988731.25, 8988835.9375, 8991071.875, 8992198.4375, 8992881.25, 8995573.4375, 9004615.625, 9030792.1875, 9030960.9375, 9031095.3125, 9031282.8125, 9032473.4375, 9036792.1875, 9036856.25, 9039098.4375, 9039409.375, 9039762.5, 9039989.0625, 9040539.0625, 9040915.625, 9040976.5625, 9041345.3125, 9041409.375, 9043701.5625, 9044026.5625, 9044132.8125, 9046395.3125, 9046598.4375, 9046812.5, 9046817.1875, 9046912.5, 9047639.0625, 9049565.625, 9049846.875, 9050868.75, 9051934.375, 9053914.0625, 9054840.625, 9056235.9375, 9058725.0, 9061704.6875, 9069498.4375, 9069506.25, 9071943.75, 9074067.1875, 9075142.1875, 9076225.0, 9079643.75, 9081614.0625, 9081675.0, 9081676.5625, 9082278.125, 9082650.0, 9083378.125, 9084345.3125, 9084979.6875, 9085076.5625, 9085618.75, 9086126.5625, 9086673.4375, 9086779.6875, 9087821.875, 9088042.1875, 9088301.5625, 9088359.375, 9088392.1875, 9088423.4375, 9088679.6875, 9088706.25, 9089032.8125, 9089089.0625, 9089090.625, 9089173.4375, 9089217.1875, 9089218.75, 9089262.5, 9089264.0625, 9089270.3125, 9089279.6875, 9089284.375, 9089296.875, 9089362.5, 9089504.6875, 9089829.6875, 9090100.0, 9090242.1875, 9090256.25, 9090464.0625, 9090473.4375, 9090717.1875, 9090735.9375, 9091303.125, 9091325.0, 9091785.9375, 9092210.9375, 9092214.0625, 9092325.0, 9092954.6875, 9093354.6875, 9093379.6875, 9093604.6875, 9093629.6875, 9093687.5, 9093778.125, 9093931.25, 9094551.5625, 9094818.75, 9094921.875, 9095182.8125, 9095276.5625, 9095471.875, 9095534.375, 9095685.9375, 9095809.375, 9096037.5, 9096096.875, 9096195.3125, 9096242.1875, 9096242.1875, 9096648.4375, 9096678.125, 9096814.0625, 9096831.25, 9096970.3125, 9097176.5625, 9097214.0625, 9097537.5, 9097664.0625, 9097728.125, 9097748.4375, 9098431.25, 9098443.75, 9098864.0625, 9098931.25, 9099078.125, 9099678.125, 9099709.375, 9099714.0625, 9105495.3125, 9108160.9375, 9119932.8125, 9122787.5, 9124354.6875, 9124364.0625, 9124396.875, 9124523.4375, 9124770.3125, 9124803.125, 9125487.5, 9125600.0, 9125693.75, 9126878.125, 9128679.6875, 9128707.8125, 9128709.375, 9129215.625, 9129317.1875, 9130765.625, 9131331.25, 9131907.8125, 9132014.0625, 9132043.75, 9133329.6875, 9134289.0625, 9134609.375, 9134832.8125, 9144650.0, 9153195.3125, 9153267.1875, 9154398.4375, 9156373.4375, 9157345.3125, 9160698.4375, 9160732.8125, 9161307.8125, 9163648.4375, 9164257.8125, 9164760.9375, 9165667.1875, 9165689.0625, 9166060.9375, 9166140.625, 9166407.8125, 9166967.1875, 9169387.5, 9169910.9375, 9170482.8125, 9170926.5625, 9171160.9375, 9172639.0625, 9172945.3125, 9173367.1875, 9173400.0, 9173781.25, 9175565.625, 9175992.1875, 9176331.25, 9179114.0625, 9184820.3125, 9184832.8125, 9185929.6875, 9188104.6875, 9188690.625, 9189268.75, 9189432.8125, 9189615.625, 9190246.875, 9190251.5625, 9190262.5, 9190501.5625, 9191315.625, 9191387.5, 9191881.25, 9192189.0625, 9192607.8125, 9192660.9375, 9194310.9375, 9198184.375, 9198498.4375, 9202181.25, 9205237.5, 9206037.5, 9210162.5, 9210389.0625, 9216592.1875, 9217500.0, 9221870.3125, 9223168.75, 9240759.375, 9241048.4375, 9241375.0, 9241795.3125, 9242664.0625, 9243109.375, 9243143.75, 9244673.4375, 9244793.75, 9244818.75, 9245759.375, 9248735.9375, 9249584.375, 9249587.5, 9250393.75, 9251046.875, 9251048.4375, 9251504.6875, 9252232.8125, 9252876.5625, 9252981.25, 9254890.625, 9255895.3125, 9265015.625, 9266023.4375, 9279628.125, 9287020.3125, 9287171.875, 9290442.1875, 9290993.75, 9291664.0625, 9293692.1875, 9293726.5625, 9294456.25, 9294507.8125, 9294514.0625, 9298482.8125, 9298946.875, 9299182.8125, 9299525.0, 9300209.375, 9301351.5625, 9302185.9375, 9303482.8125, 9306239.0625, 9306371.875, 9309070.3125, 9317134.375, 9348784.375, 9352528.125, 9383178.125, 9387189.0625, 9388006.25, 9388368.75, 9388903.125, 9389192.1875, 9389398.4375, 9390550.0, 9398114.0625, 9399178.125, 9400410.9375, 9407659.375, 9408317.1875, 9413340.625, 9414045.3125, 9414671.875, 9415359.375, 9420593.75, 9451454.6875, 9454281.25, 9472517.1875, 9475528.125, 9490748.4375, 9491090.625, 9491407.8125, 9491739.0625, 9492835.9375, 9493862.5, 9494696.875, 9519956.25, 9520500.0, 9520517.1875, 9520606.25, 9520757.8125, 9522117.1875, 9522321.875, 9525976.5625, 9526635.9375, 9527129.6875, 9527682.8125, 9529068.75, 9537310.9375, 9547309.375, 9549134.375, 9551140.625, 9555470.3125, 9556670.3125, 9577035.9375, 9627831.25, 9633131.25, 9633567.1875, 9634851.5625, 9634942.1875, 9634996.875, 9636200.0, 9636203.125, 9674834.375, 9686006.25, 9765542.1875, 9765554.6875, 9765781.25, 9765817.1875, 9766620.3125, 9768617.1875, 9770923.4375, 9771557.8125, 9775157.8125, 9859306.25, 10024181.25, 10025139.0625, 10030987.5, 10253185.9375, 10254517.1875, 10254604.6875, 10256017.1875, 10256078.125, 10256331.25, 10256784.375, 10256810.9375, 10256853.125, 10259625.0, 10273998.4375, 10294403.125, 10294453.125, 10296456.25, 10317260.9375, 10325451.5625, 10333953.125, 10420651.5625, 10480359.375, 10553034.375, 10554129.6875, 10555826.5625, 10580101.5625, 10601931.25, 10602167.1875, 10607145.3125, 10673568.75, 10689646.875, 10734523.4375, 10742792.1875, 11032814.0625, 11034987.5, 11067242.1875, 11070414.0625, 11073853.125, 11076217.1875, 11093540.625, 11128373.4375], [32.924657097325294, 116.6798834281576, 81.00691152143972, 13.81606582078742, 41.72555211612796, 27.067201900344223, 9.661675720852736, 82.78728982993094, 14.823666307008384, 50.47189156931825, 73.32101371470245, 16.13505519609505, 14.418534181308267, 20.83983263389171, 5.6965500646658445, 8.894401020402597, 58.60041078339818, 11.856012650062224, 57.13348803974097, 61.19822632708832, 6.018440701337877, 15.15612161580969, 5.039071960146389, 28.72195793469428, 23.10672267943091, 97.58608264704937, 75.35854863886624, 9.659575726354321, 5.685506942610424, 20.335945054826034, 58.75481869670872, 93.3671116648696, 17.897319648338065, 58.727334773427245, 53.03302349377744, 80.39872580958976, 100.51427166015104, 80.13425660560868, 21.60709772063372, 69.38740188222866, 74.5877110724207, 12.78891806195076, 9.145712155769688, 115.848675321599, 61.38400223940031, 17.525835980327578, 6.607019620649585, 37.88480676717241, 17.77627576050788, 6.198439810677378, 71.16212715782888, 13.925571693954117, 14.53270219022515, 56.47847311257261, 14.231431015597083, 5.191186985993551, 61.42556379724776, 16.260243983812476, 6.241256862918884, 7.906379731635614, 8.549370819691367, 64.65193001122084, 30.189325033217266, 5.758292771471893, 15.115693467049823, 74.04248132686806, 55.20489449149838, 5.212654762827162, 6.571187188756432, 89.76614159046558, 10.465239039799483, 5.029637402935534, 53.66821510682772, 24.985309575968472, 20.83597270999763, 60.34876815587375, 16.30591490556375, 42.063002244972274, 75.90222048543963, 10.728107875943651, 8.947732010090657, 5.908348804966057, 87.09635711889646, 8.024669039121248, 62.78872949930695, 18.035970787139966, 36.6274253304104, 89.57669413445683, 64.25627901150489, 125.8307875237152, 12.03887369343171, 60.13972059924853, 12.465507224304247, 125.15352480115966, 88.92581524380074, 18.23269855016324, 97.77919819247317, 164.12279828797415, 5.500703622250263, 51.410707353117, 12.308662345715513, 5.551288558591052, 7.729951864461873, 16.779813778676992, 114.01548200973278, 6.5225869783925114, 30.922506918927596, 62.74326095199561, 23.335546081325806, 74.26539906834553, 40.56839744776385, 47.17016024195928, 5.87290474759791, 20.994805317749, 27.642249403970858, 13.73910176767776, 21.1696653546308, 9.488463220362247, 50.98497923456754, 16.613098978222425, 68.44332311093979, 50.416608986945874, 46.39505885148927, 80.97850407183961, 5.856558114020477, 5.226520871294633, 27.797452185396622, 59.124549090495805, 25.49913342130373, 17.582139815441657, 12.727219598469084, 211.41610012978518, 17.15872374385065, 40.74508636717269, 23.595943065991705, 34.332486437982745, 12.870203049543697, 28.44004240532065, 46.15426647738029, 57.41312617612931, 74.0944622169626, 108.3079182877103, 128.15072494588716, 22.762698457694675, 59.56549031715537, 77.67081983330007, 71.28092089875014, 109.1983634950865, 5.087629448471466, 5.02612513931732, 172.80589021724447, 78.51851359915393, 74.45847793068981, 59.61830626984404, 42.12325779692252, 80.06486479088169, 59.05093839089593, 43.18937772131305, 68.07164278428874, 35.42035311134909, 17.079249324906534, 7.61054281595935, 62.004903548381776, 84.10756466593455, 115.36061851669324, 52.9216767045551, 96.82843774883669, 20.884818239873937, 7.668040580038916, 7.11332217615482, 14.344377386407078, 38.745995673451446, 5.48425647205571, 18.60014119719927, 24.366283839026533, 143.428561693925, 48.59463022082683, 11.138488383048657, 11.01047232136609, 7.041889125072349, 5.1902418741117105, 31.636465234104293, 66.70689770171029, 98.21874010969626, 10.942707948514915, 87.98031492065158, 78.20688103665884, 75.80383253753945, 56.17196284850187, 45.62371635788635, 74.53985454164697, 100.36954080425524, 14.67279338980616, 21.47146830699488, 11.602682034196775, 32.61574729027931, 87.73579343352041, 14.670771706702002, 20.43916203688921, 6.006154518156957, 21.416024844655446, 20.244868802468304, 6.70066476499348, 75.10459147983843, 7.1079328292833415, 33.78828705240565, 68.33043595502669, 7.088111775935066, 5.2592237910222615, 41.957385930903285, 16.951986633866152, 32.30809348890219, 72.6679734717716, 61.608331974283146, 66.31096207653296, 43.56543120750834, 35.0035442328987, 17.938597863334167, 10.3130314100056, 97.39706789239654, 22.851046801484483, 12.129908465551022, 37.49006667836005, 16.3947930434081, 15.934215058800195, 151.14226456567974, 12.914239952692064, 16.337265247714317, 11.844890387450622, 6.944662971350165, 6.08410190696702, 25.40423052074012, 43.533659711533936, 8.805757463629188, 80.8808191305475, 27.18252032850513, 125.39988663636795, 24.19919465587248, 5.887121302681618, 165.90263302406132, 34.508508498647906, 20.561151933922357, 36.86235548799108, 90.92139993504784, 50.05667601548069, 25.487998106876354, 141.48940988255507, 7.243799439184941, 6.587818403883618, 121.46000486333148, 55.15767294459172, 21.486340336232416, 25.48974125085641, 23.477542479703597, 14.320472620269658, 21.497624785132977, 15.566253335702699, 44.19846895731496, 5.167972977641077, 8.081175676856867, 126.0309991132097, 14.978635612315523, 16.28236732110962, 7.583503821865062, 82.17689781359329, 41.66946614246198, 14.4587787266855, 100.85049292882078, 30.768355277168403, 22.96530113646394, 16.35121506637152, 39.62119117255041, 8.525017105057323, 60.40707931082138, 10.87866388263598, 27.565435484921743, 24.571254282781453, 21.588018015068272, 128.2668210825618, 27.624151243146777, 114.22089447212235, 11.800058439264006, 16.53051006309701, 30.692140727846795, 41.678786439275264, 76.49965600994615, 30.84402499776383, 45.286576270859285, 20.842401000776654, 49.414366725282974, 11.70703496129263, 13.022546636047572, 11.838428320057156, 27.90432821271198, 47.3357336198286, 9.649760613378662, 33.6895291039139, 42.50795531437848, 11.675433502415846, 173.7987759447381, 104.9719710209032, 40.1259550382697, 10.10100506682189, 26.551873548215895, 8.996573066915651, 32.41341909816228, 33.071746696685985, 38.46286975728948, 28.650557694252885, 79.23470504918554, 121.83991317282812, 45.60153871703153, 24.44779696699885, 39.434282884618405, 74.24008279882129, 19.20779174963493, 51.38379060714112, 64.4015378613884, 43.19510912778543, 59.084060065371446, 38.81403666317835, 18.83579110389391, 49.96896346212341, 55.92919405936321, 57.788837123061114, 49.59490519928747, 96.29733125882511, 55.169546460358234, 12.437623525684874, 15.726839210117292, 98.98227815791182, 9.82112699442779, 40.228218061072845, 46.81303675477902, 110.09165585815043, 100.2597246318238, 31.03538933586202, 15.317744582333278, 37.43139119613992, 59.52038818983135, 48.22109401350293, 32.717217153901444, 25.865032447816056, 13.208018311005814, 5.9003606704665055, 27.379948933004183, 14.479266667070963, 75.57768193779717, 18.589234915344335, 15.722342295291924, 32.85661532047301, 99.89236781071969, 37.424265835090324, 92.30584422484525, 92.05991887316418, 15.61320573818369, 8.63095863393564, 8.760347537697715, 33.476510617181674, 18.827789038991863, 26.888165929433242, 30.42175663256922, 16.069880451532686, 5.262989843630826, 40.03971070819298, 14.134046835277053, 17.020048795183868, 50.50387825877301, 151.22889611269392, 63.46524483034014, 16.14903298780898, 5.667255884978681, 9.629352873382768, 44.689850219089266, 17.669481666425575, 55.332010720627196, 13.73043771923356, 19.473972904421682, 41.86750954546676, 66.08167162159023, 10.620321132711709, 76.37485386571336, 65.41713063108222, 11.298951269996039, 41.64482405005603, 84.65469006238222, 81.02903411316058, 27.52793745658734, 42.474287791121014, 12.715564225568837, 15.652578630488666, 37.14309259239211, 27.71625870172194, 15.152537691637114, 23.096684496603803, 39.686267491266236, 34.016309444265524, 14.873539009205832, 18.445964371100416, 189.78456389753975, 39.92180615912983, 8.261552645615602, 15.588276082733929, 60.55771075958826, 29.076066288314937, 108.79209425868126, 20.710337521655884, 47.36582564176112, 54.92282267726806, 25.410236199615987, 30.261546145513798, 24.772628264234005, 88.90794795481821, 11.64244701970497, 192.50614707089505, 28.675964727111598, 13.378745903290769, 12.875199809248988, 70.8817727454168, 14.232152521302606, 67.66386061235895, 131.5717429906109, 10.723698552120933, 85.4214145050205, 27.11142764702241, 29.30161052447525, 6.792385579804992, 14.465028927842702, 86.2864948315138, 26.886914076650385, 29.165727144989788, 36.382368188891824, 123.3290671831657, 35.47817831868368, 47.65733211960614, 64.52371765098714, 9.119815875592273, 77.32845423359748, 19.839212661948206, 130.86225844963167, 16.41030622824398, 44.63710634732309, 33.40569040456431, 58.37489605488815, 27.812340684483356, 5.024411214800042, 41.80996621495609, 19.28292084063712, 14.097641579230096, 40.59118815383231, 42.9961004748717, 7.6652682804176635, 68.15295699729796, 104.18122702871335, 70.46224013378512, 68.64033638212265, 14.388731608272009, 172.00391572644023, 50.801326739344915, 22.300684320065606, 30.302611871915524, 20.740584225982126, 15.984587492122163, 65.79103634635933, 61.28101592649081, 8.969244101743984, 56.69863098235196, 36.39553605323646, 19.938333112624754, 40.852170448050344, 31.234427155665674, 5.8190971116115, 182.8420777150173, 11.111837386432626, 63.616456844031106, 56.61250372517559, 11.663793021300107, 10.089796506861942, 123.35194978335716, 79.62773183223466, 17.83898689457889, 7.998607243290227, 43.07773779308858, 67.78000143328207, 43.154529947723304, 26.75180031867474, 52.47549176457697, 40.472853778637734, 96.37876472137663, 45.39333692737364, 22.713574577667814, 58.752260654513776, 21.406971700910937, 50.74241546664142, 61.12786118575792, 27.19324336266716, 23.699518848411273, 116.91549173960323, 66.37872713013465, 32.71818471615478, 34.809660207489486, 10.637679965990744, 74.71396575972813, 13.976734496700999, 132.71265320995658, 30.94502928732678, 11.661705287161318, 6.033893429643892, 24.661209020386643, 26.283522644497964, 89.35037805921596, 35.47240006511114, 50.42968566431969, 36.208889858415404, 10.019641602200524, 65.26607823632742, 24.94490476700947, 24.42131428353687, 10.208898612059109, 15.088051621424992, 36.67242457886821, 47.77566024903621, 36.07453992963525, 50.71645487604973, 7.5307774819693245, 39.04070543560128, 64.52297220800261, 26.81748508265912, 40.16312478517369, 15.454442383661968, 115.37027191545279, 134.6126938270227, 64.0158972814034, 8.261571179151076, 13.56233045772898, 118.5113676453322, 99.12068873961482, 116.65802870113792, 60.79419175324312, 8.85095042406844, 33.0001338894392, 29.848726291123196, 80.978850633698, 37.42312878248814, 41.49936258411806, 6.208724910462915, 9.038345578817697, 5.378216745630152, 12.599176265682507, 9.576960779908998, 37.9768887904705, 20.57193689307942, 7.47054296424512, 28.972645932052327, 51.946150292117785, 16.455350704646413, 6.622425127852205, 121.03723275264161, 74.04103502061801, 115.62040635983799, 55.15835805526446, 18.153447245211975, 26.568445284829195, 17.804512556224864, 54.87089187103154, 79.43299725586613, 73.00508114607749, 8.05672486198217, 68.66837907266587, 131.2828158163505, 39.748842547621585, 26.00035424728346, 36.53023358958911, 45.923752593590265, 65.54674609903553, 16.936684315190387, 151.90447416488368, 19.232006465704693, 88.39092742095102, 28.570599775977055, 96.83386437787198, 68.94188119201912, 39.2016901752217, 8.950548092181926, 158.03791325754122, 71.23842625613656, 37.035100116729225, 78.22005215695403, 5.103273678674295, 5.352198469536739, 14.744225916779666, 8.905076504710934, 76.36010568699902, 43.02168149222666, 27.268192235020923, 30.474329458927773, 18.826413672841, 25.436450950583133, 28.37893010368078, 79.82336219502548, 5.190963802034682, 6.830287051783128, 13.836966409664848, 51.437551372011285, 57.42260720854609, 7.360613115738645, 22.89453082897128, 23.85191904698461, 97.29391762377679, 21.49236196980726, 70.16292413629075, 40.08199733977131, 63.94249628045913, 29.131857014537644, 25.198277011720318, 39.598212776084644, 62.81178636082504, 25.823389893061904, 66.04055978992595, 100.76757322412215, 33.63814292770292, 12.537869524972136, 83.65746169930932, 111.68595237956785, 19.039552547928842, 262.2673539374696, 90.5290369595815, 25.205566693605142, 22.51922230158157, 25.228548720088604, 90.22342466673786, 31.618407060578242, 101.39848925022315, 58.4626720075059, 57.05113638378787, 73.68164603567162, 175.76699501964458, 30.60912336745204, 112.53885813476577, 52.01525922659552, 20.742798833424285, 19.4055430654669, 99.59278737399248, 7.1539607964270076, 57.52491950535498, 62.91495774737473, 68.55294170451401, 5.065912136159971, 11.70429541438887, 44.03004171602744, 54.01662279964817, 7.6301049148870925, 24.156903634909927, 92.74577462411881, 93.96752856229226, 125.42717094001031, 89.06086857117145, 14.092339971336317, 8.366229186720703, 29.86870122150375, 81.16301792941056, 26.270465451492335, 22.588621813012146, 15.0348306778534, 12.135946702055199, 7.362079397250393, 65.48394931434993, 8.351663089339961, 24.475895806110724, 112.68030759982324, 24.218606158559993, 15.468194403853039, 6.183517318434525, 22.38633697814639, 24.00400358928337, 61.040960658143604, 17.95748969298244, 21.95179224030898, 12.931370348004553, 12.74110154457865, 48.93929917469359, 227.67706996602868, 112.51998434093974, 41.582941685496486, 96.74943666517736, 67.81621497977122, 50.690291204949254, 57.622419757448895, 47.072283883215206, 65.72980877496182, 26.610746086471913, 21.397502792242683, 37.52089472206346, 25.20553236446617, 16.676651704488503, 27.99152455666368, 60.36288306624489, 13.01485261568247, 5.61294970333641, 14.577076388617769, 9.580380276606364, 113.32975505283501, 52.87055021224667, 126.0882688787489, 89.92300040829763, 46.63568302929067, 38.79100329453242, 19.98976518116541, 6.888931911550446, 9.979114224176064, 101.72348327452988, 5.5675429672881185, 72.26115133243339, 43.41732545177426, 108.14565597231343, 19.329126093552816, 21.28846306528019, 12.220944754374417, 13.078644239520068, 137.24892675503747, 8.109864531215313, 70.02273076282208, 60.26695669772552, 118.82262364836316, 24.724064203774613, 5.555967513944521, 28.328050168583378, 93.28746462161298, 52.38185693816583, 149.3918214312369, 46.401212688278754, 141.4068641054947, 245.12840350386492, 32.47302767718489, 18.900967615763058, 24.08344403997618, 6.110924514752332, 63.57555515550892, 61.110101710307035, 32.34003870050973, 119.41474780228708, 19.919834974370858, 101.66238449718038, 8.440986096295617, 25.686160201138723, 25.757665396587274, 31.33793909630443, 30.22003949943895, 57.872915494078455, 59.295995251688666, 23.80224850446178, 9.473432820587302, 53.53458669037609, 210.0824999246817, 19.351889690919258, 72.09739342401426, 19.192115014275775, 6.184879174974784, 86.94955722549035, 16.64295476996734, 84.18898657715172, 108.72345167435368, 13.616465639399365, 5.252983170080755, 21.79300048815033, 81.54846902119787, 26.33760758170059, 207.0198146956072, 12.186649544505443, 8.141526326198806, 9.70631881222407, 23.72574446910234, 57.8872899125006, 25.86148399056148, 29.51829299553497, 67.48981852624637, 19.41387128134355, 25.769864105974953, 33.74541501654854, 15.632634922910334, 44.73308627702366, 14.159321979927359, 37.398354260490876, 104.86826328920962, 16.53196829735078, 69.97091518926179, 19.801065962390993, 75.03072688972566, 16.9822742979065, 6.5496748682950825, 5.4626251441583005, 14.577347472508464, 17.221316960507394, 37.1947214574466, 89.31249921380503, 26.03248583382198, 10.199416811522697, 49.295336627075514, 22.663647664048113, 90.0617620885011, 13.385534696093494, 15.71123891577353, 5.6769395033954675, 69.79158888822076, 99.11028403624456, 126.74564138490084, 35.62904257849863, 8.716734065067033, 103.79829180561171, 23.318639176450983, 9.968603535685736, 55.64024346618522, 76.21506179656281, 17.400712255662746, 72.86711976442741, 7.787611774742526, 17.554050112901276, 53.93162568684195, 48.86700727739273])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)