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 = 43936
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);
([5593895.3125, 5848604.6875, 5991420.3125, 5991457.8125, 6033123.4375, 6252290.625, 6325731.25, 6723187.5, 7058406.25, 7075506.25, 7096498.4375, 7144004.6875, 7144901.5625, 7351868.75, 7464667.1875, 7485085.9375, 7838417.1875, 7894279.6875, 7895776.5625, 7937912.5, 7988645.3125, 7988650.0, 7988751.5625, 7989456.25, 8031082.8125, 8036070.3125, 8036851.5625, 8036903.125, 8043421.875, 8083312.5, 8085968.75, 8111518.75, 8331220.3125, 8332556.25, 8333475.0, 8335776.5625, 8338275.0, 8339084.375, 8339089.0625, 8339548.4375, 8343417.1875, 8346607.8125, 8349054.6875, 8352450.0, 8352742.1875, 8352748.4375, 8356268.75, 8356965.625, 8357035.9375, 8357364.0625, 8357523.4375, 8358268.75, 8358459.375, 8358818.75, 8359064.0625, 8362582.8125, 8365468.75, 8366692.1875, 8367253.125, 8367451.5625, 8369973.4375, 8374234.375, 8374492.1875, 8382121.875, 8383051.5625, 8394087.5, 8397353.125, 8398828.125, 8399879.6875, 8560350.0, 8562370.3125, 8568525.0, 8588934.375, 8589625.0, 8597015.625, 8597246.875, 8597325.0, 8597395.3125, 8597509.375, 8597515.625, 8597676.5625, 8597693.75, 8597773.4375, 8597793.75, 8597793.75, 8597929.6875, 8597971.875, 8598359.375, 8598571.875, 8598793.75, 8599925.0, 8599967.1875, 8600562.5, 8601332.8125, 8601501.5625, 8601715.625, 8602268.75, 8602387.5, 8602512.5, 8602610.9375, 8602789.0625, 8603528.125, 8603623.4375, 8603831.25, 8603845.3125, 8603954.6875, 8604090.625, 8604492.1875, 8604570.3125, 8604593.75, 8604860.9375, 8605042.1875, 8605314.0625, 8605314.0625, 8605412.5, 8605598.4375, 8605678.125, 8605745.3125, 8605973.4375, 8606071.875, 8606131.25, 8606209.375, 8606476.5625, 8607101.5625, 8607170.3125, 8607743.75, 8607984.375, 8608126.5625, 8608376.5625, 8608420.3125, 8608643.75, 8608732.8125, 8608751.5625, 8608839.0625, 8609157.8125, 8609217.1875, 8609248.4375, 8609382.8125, 8609473.4375, 8609559.375, 8609728.125, 8609860.9375, 8609925.0, 8610054.6875, 8610226.5625, 8610373.4375, 8610478.125, 8610514.0625, 8610521.875, 8610570.3125, 8610607.8125, 8610760.9375, 8610820.3125, 8610862.5, 8611165.625, 8611376.5625, 8611543.75, 8611812.5, 8611954.6875, 8611981.25, 8612028.125, 8612139.0625, 8612292.1875, 8612434.375, 8612534.375, 8612551.5625, 8612556.25, 8612695.3125, 8612850.0, 8612900.0, 8613046.875, 8613390.625, 8613934.375, 8614112.5, 8614123.4375, 8614309.375, 8614465.625, 8614657.8125, 8614934.375, 8615059.375, 8615140.625, 8615143.75, 8615310.9375, 8615440.625, 8615479.6875, 8615593.75, 8615596.875, 8615651.5625, 8615771.875, 8615776.5625, 8615784.375, 8615976.5625, 8615979.6875, 8616081.25, 8616118.75, 8616346.875, 8616454.6875, 8616462.5, 8617076.5625, 8617235.9375, 8617242.1875, 8617259.375, 8617282.8125, 8617535.9375, 8617715.625, 8617870.3125, 8617928.125, 8617951.5625, 8617954.6875, 8618084.375, 8618239.0625, 8618285.9375, 8618287.5, 8619643.75, 8620910.9375, 8621507.8125, 8621828.125, 8622807.8125, 8623387.5, 8623506.25, 8623564.0625, 8623589.0625, 8623643.75, 8623851.5625, 8623942.1875, 8624025.0, 8624106.25, 8624159.375, 8624160.9375, 8624243.75, 8624310.9375, 8624312.5, 8624457.8125, 8624510.9375, 8624684.375, 8624834.375, 8624901.5625, 8624934.375, 8624956.25, 8624957.8125, 8625057.8125, 8625181.25, 8625253.125, 8625276.5625, 8625315.625, 8625437.5, 8625504.6875, 8625843.75, 8626071.875, 8626560.9375, 8627154.6875, 8627845.3125, 8627909.375, 8628589.0625, 8631904.6875, 8634373.4375, 8635718.75, 8638992.1875, 8639064.0625, 8639131.25, 8639229.6875, 8639239.0625, 8639295.3125, 8639509.375, 8639554.6875, 8639557.8125, 8639650.0, 8639660.9375, 8639681.25, 8639704.6875, 8639909.375, 8639957.8125, 8639968.75, 8639968.75, 8640026.5625, 8640054.6875, 8640056.25, 8640070.3125, 8640078.125, 8640096.875, 8640101.5625, 8640114.0625, 8640121.875, 8640148.4375, 8640346.875, 8640551.5625, 8640821.875, 8640853.125, 8640864.0625, 8640918.75, 8641004.6875, 8641307.8125, 8641389.0625, 8641429.6875, 8641551.5625, 8641689.0625, 8641690.625, 8641732.8125, 8641734.375, 8641784.375, 8641801.5625, 8641803.125, 8641809.375, 8641856.25, 8641871.875, 8641934.375, 8641950.0, 8641962.5, 8642003.125, 8642121.875, 8642135.9375, 8642307.8125, 8642381.25, 8642425.0, 8642548.4375, 8642706.25, 8642795.3125, 8642837.5, 8642857.8125, 8642873.4375, 8642931.25, 8643020.3125, 8643071.875, 8643275.0, 8643392.1875, 8643485.9375, 8643503.125, 8643525.0, 8643614.0625, 8643629.6875, 8643707.8125, 8643712.5, 8643714.0625, 8643725.0, 8643778.125, 8643850.0, 8643873.4375, 8643876.5625, 8643940.625, 8643942.1875, 8643993.75, 8644101.5625, 8644153.125, 8644196.875, 8644204.6875, 8644389.0625, 8644764.0625, 8645134.375, 8645153.125, 8645296.875, 8645540.625, 8645634.375, 8645687.5, 8645875.0, 8645884.375, 8646034.375, 8646095.3125, 8646123.4375, 8646132.8125, 8646295.3125, 8646315.625, 8646350.0, 8646592.1875, 8646828.125, 8646903.125, 8647328.125, 8647479.6875, 8647510.9375, 8647725.0, 8648160.9375, 8648223.4375, 8648718.75, 8648848.4375, 8648975.0, 8649401.5625, 8649403.125, 8649503.125, 8649507.8125, 8650290.625, 8650360.9375, 8650614.0625, 8650679.6875, 8650746.875, 8651656.25, 8652153.125, 8652303.125, 8652901.5625, 8653271.875, 8653300.0, 8653795.3125, 8654400.0, 8654492.1875, 8654612.5, 8654826.5625, 8654914.0625, 8655090.625, 8655118.75, 8655179.6875, 8655368.75, 8655568.75, 8655643.75, 8655645.3125, 8655673.4375, 8655712.5, 8655748.4375, 8655782.8125, 8655796.875, 8655801.5625, 8655803.125, 8655812.5, 8655820.3125, 8655832.8125, 8655851.5625, 8655860.9375, 8655929.6875, 8655951.5625, 8655965.625, 8656018.75, 8656031.25, 8656070.3125, 8656070.3125, 8656073.4375, 8656109.375, 8656131.25, 8656135.9375, 8656203.125, 8656214.0625, 8656226.5625, 8656228.125, 8656276.5625, 8656290.625, 8656321.875, 8656350.0, 8656351.5625, 8656359.375, 8656375.0, 8656378.125, 8656384.375, 8656392.1875, 8656435.9375, 8656439.0625, 8656443.75, 8656453.125, 8656454.6875, 8656476.5625, 8656479.6875, 8656501.5625, 8656506.25, 8656507.8125, 8656512.5, 8656534.375, 8656537.5, 8656543.75, 8656546.875, 8656548.4375, 8656548.4375, 8656556.25, 8656567.1875, 8656581.25, 8656589.0625, 8656592.1875, 8656600.0, 8656612.5, 8656614.0625, 8656626.5625, 8656639.0625, 8656650.0, 8656660.9375, 8656687.5, 8656696.875, 8656712.5, 8656717.1875, 8656728.125, 8656745.3125, 8656748.4375, 8656767.1875, 8656771.875, 8656798.4375, 8656820.3125, 8656846.875, 8656848.4375, 8656867.1875, 8656900.0, 8656904.6875, 8656956.25, 8656957.8125, 8656962.5, 8656981.25, 8656984.375, 8656984.375, 8656989.0625, 8657048.4375, 8657064.0625, 8657075.0, 8657075.0, 8657106.25, 8657143.75, 8657160.9375, 8657192.1875, 8657214.0625, 8657231.25, 8657248.4375, 8657257.8125, 8657284.375, 8657287.5, 8657318.75, 8657345.3125, 8657354.6875, 8657378.125, 8657395.3125, 8657442.1875, 8657456.25, 8657510.9375, 8657551.5625, 8657551.5625, 8657585.9375, 8657595.3125, 8657703.125, 8657706.25, 8657728.125, 8657831.25, 8657860.9375, 8657882.8125, 8657926.5625, 8658018.75, 8658021.875, 8658071.875, 8658106.25, 8658106.25, 8658134.375, 8658159.375, 8658181.25, 8658184.375, 8658189.0625, 8658196.875, 8658198.4375, 8658198.4375, 8658212.5, 8658221.875, 8658248.4375, 8658314.0625, 8658320.3125, 8658354.6875, 8658373.4375, 8658385.9375, 8658392.1875, 8658451.5625, 8658535.9375, 8658537.5, 8658607.8125, 8658612.5, 8658676.5625, 8658676.5625, 8658690.625, 8658726.5625, 8658771.875, 8658789.0625, 8658790.625, 8658798.4375, 8658871.875, 8658885.9375, 8658887.5, 8658921.875, 8658970.3125, 8659007.8125, 8659164.0625, 8659228.125, 8659240.625, 8659242.1875, 8659350.0, 8659401.5625, 8659603.125, 8659629.6875, 8659685.9375, 8659740.625, 8660125.0, 8660170.3125, 8660189.0625, 8660321.875, 8660390.625, 8660539.0625, 8660704.6875, 8661618.75, 8661754.6875, 8661856.25, 8662339.0625, 8662414.0625, 8662467.1875, 8662523.4375, 8662707.8125, 8662884.375, 8662892.1875, 8662892.1875, 8662904.6875, 8662917.1875, 8662957.8125, 8662959.375, 8662981.25, 8663121.875, 8663214.0625, 8663284.375, 8663306.25, 8663317.1875, 8663321.875, 8663554.6875, 8663564.0625, 8663570.3125, 8663717.1875, 8663740.625, 8663806.25, 8663814.0625, 8663842.1875, 8663871.875, 8663907.8125, 8663937.5, 8664001.5625, 8664003.125, 8664014.0625, 8664082.8125, 8664095.3125, 8664101.5625, 8664121.875, 8664139.0625, 8664146.875, 8664157.8125, 8664171.875, 8664200.0, 8664207.8125, 8664218.75, 8664273.4375, 8664375.0, 8664375.0, 8664392.1875, 8664393.75, 8664423.4375, 8664457.8125, 8664459.375, 8664464.0625, 8664467.1875, 8664473.4375, 8664481.25, 8664515.625, 8664532.8125, 8664534.375, 8664628.125, 8664645.3125, 8664682.8125, 8664703.125, 8664703.125, 8664740.625, 8664778.125, 8664825.0, 8664828.125, 8664842.1875, 8664856.25, 8664909.375, 8664992.1875, 8664995.3125, 8665000.0, 8665021.875, 8665059.375, 8665114.0625, 8665132.8125, 8665140.625, 8665179.6875, 8665262.5, 8665293.75, 8665337.5, 8665423.4375, 8665442.1875, 8665478.125, 8665503.125, 8665520.3125, 8665525.0, 8665551.5625, 8665556.25, 8665560.9375, 8665606.25, 8665609.375, 8665628.125, 8665657.8125, 8665732.8125, 8665751.5625, 8665787.5, 8665792.1875, 8665801.5625, 8665810.9375, 8665825.0, 8665828.125, 8665831.25, 8665837.5, 8665839.0625, 8665840.625, 8665845.3125, 8665848.4375, 8665851.5625, 8665860.9375, 8665878.125, 8665900.0, 8665906.25, 8665909.375, 8665935.9375, 8665948.4375, 8665948.4375, 8665962.5, 8665964.0625, 8665965.625, 8665973.4375, 8665976.5625, 8666023.4375, 8666035.9375, 8666062.5, 8666168.75, 8666179.6875, 8666184.375, 8666214.0625, 8666221.875, 8666254.6875, 8666262.5, 8666315.625, 8666328.125, 8666368.75, 8666412.5, 8666417.1875, 8666426.5625, 8666429.6875, 8666432.8125, 8666437.5, 8666440.625, 8666442.1875, 8666442.1875, 8666445.3125, 8666453.125, 8666465.625, 8666468.75, 8666496.875, 8666525.0, 8666534.375, 8666545.3125, 8666551.5625, 8666562.5, 8666582.8125, 8666598.4375, 8666603.125, 8666628.125, 8666632.8125, 8666651.5625, 8666701.5625, 8666718.75, 8666721.875, 8666721.875, 8666725.0, 8666743.75, 8666748.4375, 8666754.6875, 8666784.375, 8666807.8125, 8666834.375, 8666839.0625, 8666842.1875, 8666879.6875, 8666920.3125, 8666981.25, 8666985.9375, 8667065.625, 8667065.625, 8667198.4375, 8667201.5625, 8667248.4375, 8667307.8125, 8667567.1875, 8667756.25, 8671078.125, 8673585.9375], [20.022511874343287, 78.46195819126174, 16.4882159166492, 8.3806210227965, 22.655731992662222, 32.7894390233988, 19.992122006824662, 266.3683665547642, 37.971975006894496, 10.479829773658906, 69.07416418976844, 21.18869271607582, 211.06424103209028, 55.74804504420147, 47.64310231793246, 99.0119021928493, 57.440755801619154, 24.558960119136263, 14.00753740105926, 44.73793927065504, 11.942165738710592, 109.35471652356631, 9.115883724207817, 263.72326672636524, 83.74119530890762, 29.438975596620367, 27.96239366700433, 5.575536345699506, 16.55432859967531, 38.05377785665454, 9.517965289059601, 16.898529382988475, 54.086802474986435, 35.035608059892716, 65.39603760722872, 82.51020713418367, 27.12399671080831, 52.15975136926483, 57.12103036787454, 80.35785084335491, 109.76880160587254, 32.347070633102255, 43.79900953922901, 43.001750458424894, 121.31265616031341, 20.550560465803418, 29.589548108661166, 55.18851215289648, 12.639908647972959, 23.31526852045658, 110.79713508812934, 49.80955969255578, 74.46663672128123, 65.97957744967745, 30.964319590318677, 9.090094185837403, 11.307750037033616, 13.188127574503481, 103.05033685647533, 86.36100595153673, 69.76472482982308, 26.248142158722764, 6.722830443276934, 32.04335015358831, 11.671207498930526, 90.95072179078244, 29.672800258648127, 17.713107124165145, 34.586679194393454, 54.28498840904236, 24.533642372448206, 134.8606402561906, 7.0586037825888175, 12.59370809350691, 7.599111804700165, 51.108265454936515, 10.574791549462693, 11.627843477375071, 15.633348816853314, 15.456471071424817, 7.336746640943726, 27.238365400571542, 35.16720754684803, 27.32120359113039, 75.01455418802523, 20.246664143806306, 66.01829339722121, 15.354966589129722, 87.80361657410519, 6.567474863368416, 11.207176214421898, 8.827632784939132, 31.266688100606867, 9.047613354473457, 13.77135809467871, 18.392743426337734, 23.68670557234559, 33.47930157405227, 126.18021487531374, 50.54205196762329, 17.167018492401144, 15.107186850625812, 97.11346589073548, 41.092014407007355, 68.52099390329809, 16.809344891127335, 62.09238777153939, 9.081398378053724, 63.7713158004524, 85.05809963823498, 32.89963711215362, 6.28032830461624, 10.239256394195609, 84.25618049742822, 74.19170638390517, 20.625928059181273, 20.658027950333345, 18.65693408190566, 5.119125621013834, 57.32187041853211, 66.61046769498147, 47.714857279115904, 16.85534833666673, 30.463125340182778, 30.78821050302901, 36.34171075723825, 84.02079413917467, 32.60556176778896, 51.326718497721664, 59.72526238174085, 12.888860951436305, 27.429538697979616, 17.475360841031318, 13.235281348158273, 12.719850189976258, 59.231161340322366, 7.279705078240602, 15.017377153977073, 10.215416068621655, 28.180066957549535, 45.31755918376416, 11.73443880818547, 50.987476867285295, 106.18768313861088, 20.471439204566476, 9.42814307637266, 64.50711046847144, 61.38188200713779, 10.576642261587493, 92.61795480621976, 21.370780776051433, 49.483394211513556, 8.242707347926583, 71.82373137823643, 6.630914744746697, 22.852199563994994, 45.957270068714855, 13.845829546365433, 42.10796933532797, 11.40769711508683, 51.86207872478932, 31.56648948070083, 91.91959275006454, 12.664885596122, 66.64504524126973, 13.06646029002879, 17.741370755701762, 19.18840193266069, 14.656030646639255, 51.869781721360766, 76.54557599956274, 106.9839653635072, 59.070445608763066, 52.27023388168509, 43.571569797211836, 136.81591654951663, 29.209323085934436, 82.29972974349566, 8.513164735397329, 5.241121877643787, 27.070803994854682, 35.18780783075982, 26.667265187721117, 19.425925438821704, 108.90220960543692, 30.088448532601333, 77.88906824306977, 23.921220020601616, 19.211913475387956, 7.822580381785889, 79.3586105724116, 41.62256024242084, 16.731965433250142, 28.216761811564197, 48.745783913234575, 144.5740276545098, 108.36258700257906, 30.314148750774542, 94.16480367452326, 12.979280073092603, 6.579631609900123, 10.102405849151372, 11.755195295016643, 25.719677945668785, 14.000364461736062, 6.217876755306642, 8.347533360605425, 21.245435036259742, 60.14353331864014, 62.749066036269625, 113.09301290638473, 49.613538963498385, 29.933507753127913, 7.995748444568275, 5.679634866878216, 20.485660837258248, 33.78087741531586, 80.38556927958545, 5.699437769570864, 19.7616461800589, 58.41416928432731, 47.755126419368544, 42.451297655611675, 214.70992030312587, 100.9259985312806, 9.679266967878336, 64.36074394078048, 6.773238148044947, 92.54697354725297, 7.924083919981538, 8.174139286676452, 36.934720968456915, 92.93565708099888, 5.459213106317348, 6.090518003376605, 19.742142815854663, 19.803496751812347, 36.70184950543227, 27.913369717333985, 21.82117476984962, 5.125129303538489, 64.71258211212823, 51.78035997042599, 7.5458986114752795, 41.53948156003127, 20.10877436216406, 81.23550556147207, 44.467730693883595, 78.93483251147221, 16.393016965493437, 13.522621849028953, 16.318113515006356, 24.25150428408444, 69.69893291524218, 5.196451320385723, 76.60252628496639, 5.9382938197097594, 19.671632592932, 43.111224875237454, 18.504120477561496, 104.2334352893701, 28.70239436985423, 72.895920708373, 6.704828785895089, 5.219837702203113, 6.416149740213459, 14.852233857396223, 25.92187242281656, 65.9153753285124, 81.11533441930361, 12.678802646645423, 9.715760583687038, 21.40041399152576, 116.12946758162937, 30.737168083352625, 72.95060327693193, 13.737962967167903, 8.76030894904448, 33.01551040340161, 60.044895172334215, 61.58875211956814, 16.226434645763486, 23.83426538308519, 54.18206561852311, 89.00336895620045, 32.24268484031113, 45.52678353352946, 11.496322074438652, 24.851635574066623, 10.666430260161512, 28.83009946244532, 37.92612504920996, 64.4979782865023, 21.466793226517524, 24.554756428897164, 25.622392735159483, 24.396844781605573, 98.02827419184628, 26.786308615280976, 12.321624106663517, 113.4858965020849, 13.107961964124069, 36.84146593628902, 33.71028143484317, 66.78912462706931, 5.043443635812554, 8.739713791442218, 102.19965101063667, 28.457927867923743, 17.70694410850836, 91.59411401199382, 11.817212413967663, 48.606839961114005, 6.3391261383011575, 54.82502246307992, 17.651432804703134, 11.640600052826082, 7.519789398829752, 99.17138732617657, 13.809206423930991, 5.209696634599435, 15.868694825980015, 49.56455393489706, 30.926531444804127, 50.140494833105784, 5.662951429002725, 44.05934191867396, 73.3787904252479, 61.72763036999872, 17.115134546835463, 10.39227027119012, 80.00101962420351, 65.31042364056391, 27.795164747480506, 47.04273403509437, 66.44677002429101, 39.82537274026105, 20.35476358797571, 12.471259104111294, 7.408576455650318, 59.04874617140051, 22.3652553331502, 44.863944752851054, 21.8947129662407, 30.163699905595173, 54.01620857639806, 9.159586802774069, 31.885884725163745, 54.76488390261464, 61.65274490997159, 76.12458457793046, 71.02273743655472, 47.21128125138487, 17.928830725902557, 57.173070909319726, 29.944412868816933, 47.69260718562149, 8.582214411950929, 38.27455340358486, 82.84764914210425, 14.44497035508794, 11.558789004578951, 21.929138931610876, 89.60991337647583, 35.84241031802008, 204.76226066873977, 89.11441562687088, 130.3395686104191, 41.90226735708519, 10.906542596737676, 27.191609552778452, 21.4923950034853, 6.083369277025365, 6.3193460510604345, 14.628246260451183, 73.9696172451684, 17.762397614125252, 35.41798218095028, 25.170586494985553, 78.79448036693064, 51.267141346087925, 6.052412467583229, 32.01826712268636, 45.069618259105916, 15.2514955309576, 6.992247778645817, 19.575726083963634, 26.178468038852913, 77.0864992101433, 70.56467898256747, 11.69493680767102, 59.845506257241155, 100.93188083223943, 25.084907670659074, 126.3759158689946, 17.614106472524064, 5.283415402562491, 11.726348864065061, 15.937684882968709, 33.01024471737816, 89.19911701546434, 7.329428646013629, 6.442150934690001, 128.9611954871741, 28.776619523078274, 17.92625518166424, 21.607997670398746, 14.387768291917036, 11.675769363633774, 7.351455462336332, 96.88181622775627, 13.74798904204543, 42.598834825656375, 104.04222694035619, 15.224770233498267, 15.52351975734394, 13.24305784421485, 51.76451885824973, 105.9699117227311, 89.13416976708402, 50.937843533638976, 18.81046268119141, 41.21616643218582, 92.40899212100695, 10.103177935836095, 15.242602103859008, 6.525388242595328, 82.35647351597595, 25.57638234280975, 10.535537271580983, 9.902477610004253, 11.747016092488686, 32.19044362918392, 33.43058284987549, 65.65440270258851, 35.653220762959165, 94.25097240053438, 21.305260105431543, 117.94558302278995, 18.48559958881206, 33.68338644520263, 11.038231427715784, 22.026480489066678, 17.939826558337337, 20.600248528157383, 14.158737370029208, 25.575586794874734, 9.876561232705079, 109.54988620799597, 11.855428854986165, 28.68600718683858, 45.31653979728256, 77.65318337014287, 5.404245803455997, 58.560427337419554, 82.65493326315645, 103.55044672549838, 192.2821972775673, 10.89848038006748, 18.964574429713057, 29.806038730470995, 10.639721078601042, 5.853382900776596, 87.88459209376487, 12.94285883472737, 90.91315672034716, 20.448903597920555, 31.54120673670293, 18.622901973714768, 31.62538244245834, 62.76904861141467, 59.62812744885474, 17.709529335054118, 19.549264266853896, 83.0881184526006, 82.67009498700453, 92.2683732744003, 96.00571268258807, 140.65669593745503, 58.60451998251467, 84.46395574417433, 6.106678285131179, 149.13056839647655, 5.473634742081172, 24.147986864349097, 5.8269289365751, 8.72255948939593, 46.51641383412958, 19.388191873079457, 35.143627831512525, 57.81423672304632, 25.989744880207525, 28.428504908674135, 5.351621713915876, 5.962861094229473, 21.16035834130944, 12.824110878297748, 63.5973958510727, 6.162788945821083, 139.16828769483578, 84.80622533477869, 23.792076193277655, 17.556555539486602, 34.43751476099532, 34.635072595277705, 53.88655854257308, 5.070238021880104, 37.39054469587327, 77.12971427373495, 49.47619860833129, 35.99489257206772, 9.227171433471742, 16.175188402730786, 101.9419040258488, 66.31183869636145, 160.30263755650012, 82.75280998625504, 56.197924765544116, 84.60058145531889, 51.05360236113959, 8.461708011225136, 64.02861912700928, 24.469496950698353, 35.698299676169235, 120.1103699768657, 44.078659527620744, 8.810057758912102, 49.80594052248908, 60.59631290449353, 6.800482558910093, 268.1098685238079, 14.123992341776995, 18.01220299323216, 63.75094663903423, 65.77015331141875, 16.003008459238306, 45.261756339645, 59.46130310450647, 101.90110577656549, 30.95439635989883, 36.22264920973203, 5.834729884022688, 21.739053503173196, 53.55693304202912, 27.310036785786142, 15.864366766213793, 6.244712889682742, 17.81070332950543, 20.150770196073786, 19.897369770775857, 5.960205226816818, 12.235611204295266, 9.760620692507182, 10.301595669906776, 51.79188856497396, 6.661198377973035, 242.04154665597727, 70.62927489336107, 27.39152614941209, 11.669839141179615, 7.247172383162028, 8.027153950050806, 12.274792664247807, 140.94173602388585, 24.769180940669287, 78.95423260980462, 13.746892510407667, 50.58976834799397, 62.2663696708432, 6.589222718023403, 34.655756993546916, 17.661469204087584, 25.93921681623126, 41.59191780335905, 25.853977914090713, 5.674000142756098, 7.064583546260622, 16.96776611735011, 38.07506326797337, 27.894119883742512, 50.22611678101204, 22.331732811492387, 80.33384899338502, 7.453864284746731, 21.74563383116587, 43.087305145997405, 11.705325889875924, 5.13034697034523, 94.72400170687536, 31.022614469216222, 24.70087671786948, 6.535184344715774, 54.64346073784952, 36.984548369203935, 50.50368400665933, 14.312099607754234, 39.01142474040529, 164.2699329319394, 37.637661537495575, 62.92612670945739, 47.420675574632476, 59.56291967940864, 8.475325069754271, 7.362716790387138, 54.53587353043423, 9.981078612570554, 5.103098438423599, 111.31558223360828, 30.371546320290552, 135.83829324653965, 63.44798249026432, 61.4235018682927, 11.644028396324392, 103.32205672599883, 67.50617579043089, 47.201487949399336, 138.43995005943023, 33.25335871897556, 18.93688409068891, 29.39651862728524, 24.604227722655885, 66.2290043237164, 11.747061364281853, 72.00783120900645, 15.262258631532301, 194.59683263383826, 9.193671949082251, 14.333904367738203, 110.52435057405344, 7.311619127819326, 9.062697903193513, 35.14221476652692, 34.064974665538685, 6.156190680673775, 15.595471085729466, 79.67883339761102, 58.74197575479402, 52.70885025237563, 59.03949256911913, 49.539647018308855, 13.483474484748847, 5.678545366574067, 33.05941412032325, 7.894199239509674, 9.139836909229233, 28.44085182709599, 72.42039556863406, 25.130035427524316, 77.42351348495234, 65.07863038560795, 7.6947478396579445, 85.26641872112141, 84.02046636732139, 17.851658077672322, 20.39860524071606, 20.705606421693677, 66.6453174313181, 7.044914289454934, 43.39808779809452, 34.22995280665725, 12.424332790675388, 79.52004863645571, 94.01327923456759, 19.29874028501852, 74.16212097065974, 31.35814098147169, 5.633906664560219, 30.196130092479482, 22.200253961690052, 15.838434389782737, 54.607942114334925, 49.162520490964496, 13.847197336615569, 17.742091935065215, 10.676323571027847, 36.88807937831369, 9.400586886918543, 81.54197641079415, 74.4075568262742, 34.74525316187774, 8.194758611285586, 16.03138352423043, 9.891802850090373, 37.770360668300526, 25.91490705308511, 10.743371107468722, 112.78705592702772, 52.37018176569424, 99.20865864447491, 49.03835005170571, 65.55276210701308, 28.29871712038704, 75.47751554798941, 13.327607723974037, 8.95543719925393, 57.345375546623316, 57.43902798655975, 12.503607212202237, 18.351478482514604, 12.631854985408651, 54.16636366913231, 31.244070562358342, 49.815666605032014, 16.8438836565649, 7.05275400135345, 68.46363078734322, 84.41995261979798, 41.8908489682004, 76.14719386890287, 70.32925946034982, 7.258801841416832, 24.341853894424776, 64.7618300446885, 15.389104949024048, 29.24554684715484, 49.9566024019079, 45.16250262091817, 19.52449808983004, 11.352137036961382, 10.475539076057933, 108.72392502529694, 9.65589915384454, 10.153989750920672, 10.354607964319086, 25.54150175583259, 43.64366558125047, 8.104351799400517, 14.77639117722654, 20.891162055165875, 6.362054930396035, 51.19176794989235, 85.1224388289672, 45.59532921829374, 53.711103501495344, 9.035469961495378, 11.495398441954519, 23.638127990669535, 9.77640617647872, 10.090924376665647, 67.77431591196958, 44.64970585376821, 20.252866110965815, 6.917982172198404, 23.898553461192414, 71.79213771030712, 5.399232020239222, 99.50113926316546, 111.32782521083351, 13.366886557108431, 41.30704166040016, 17.584387282847256, 34.68783448152784, 25.042383812973835, 13.914366077272486, 99.62416896322168, 99.96587868132904, 27.583286254274224, 16.73111097122793, 6.1622996883229675, 12.387770139742969, 12.025053960250817, 29.23278296781624, 119.47019801700915, 31.856924793195535, 25.872075362000384, 40.61798042680902, 16.66448934095897, 43.989406796388835, 31.652349458941533, 62.517473557477764, 60.86169483950462, 59.114561262462814, 21.04670805195296, 72.20016327831769])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([5593895.3125, 5848604.6875, 5991420.3125, 5991457.8125, 6033123.4375, 6252290.625, 6325731.25, 6723187.5, 7058406.25, 7075506.25, 7096498.4375, 7144004.6875, 7144901.5625, 7351868.75, 7464667.1875, 7485085.9375, 7838417.1875, 7894279.6875, 7895776.5625, 7937912.5, 7988645.3125, 7988650.0, 7988751.5625, 7989456.25, 8031082.8125, 8036070.3125, 8036851.5625, 8036903.125, 8043421.875, 8083312.5, 8085968.75, 8111518.75, 8331220.3125, 8332556.25, 8333475.0, 8335776.5625, 8338275.0, 8339084.375, 8339089.0625, 8339548.4375, 8343417.1875, 8346607.8125, 8349054.6875, 8352450.0, 8352742.1875, 8352748.4375, 8356268.75, 8356965.625, 8357035.9375, 8357364.0625, 8357523.4375, 8358268.75, 8358459.375, 8358818.75, 8359064.0625, 8362582.8125, 8365468.75, 8366692.1875, 8367253.125, 8367451.5625, 8369973.4375, 8374234.375, 8374492.1875, 8382121.875, 8383051.5625, 8394087.5, 8397353.125, 8398828.125, 8399879.6875, 8560350.0, 8562370.3125, 8568525.0, 8588934.375, 8589625.0, 8597015.625, 8597246.875, 8597325.0, 8597395.3125, 8597509.375, 8597515.625, 8597676.5625, 8597693.75, 8597773.4375, 8597793.75, 8597793.75, 8597929.6875, 8597971.875, 8598359.375, 8598571.875, 8598793.75, 8599925.0, 8599967.1875, 8600562.5, 8601332.8125, 8601501.5625, 8601715.625, 8602268.75, 8602387.5, 8602512.5, 8602610.9375, 8602789.0625, 8603528.125, 8603623.4375, 8603831.25, 8603845.3125, 8603954.6875, 8604090.625, 8604492.1875, 8604570.3125, 8604593.75, 8604860.9375, 8605042.1875, 8605314.0625, 8605314.0625, 8605412.5, 8605598.4375, 8605678.125, 8605745.3125, 8605973.4375, 8606071.875, 8606131.25, 8606209.375, 8606476.5625, 8607101.5625, 8607170.3125, 8607743.75, 8607984.375, 8608126.5625, 8608376.5625, 8608420.3125, 8608643.75, 8608732.8125, 8608751.5625, 8608839.0625, 8609157.8125, 8609217.1875, 8609248.4375, 8609382.8125, 8609473.4375, 8609559.375, 8609728.125, 8609860.9375, 8609925.0, 8610054.6875, 8610226.5625, 8610373.4375, 8610478.125, 8610514.0625, 8610521.875, 8610570.3125, 8610607.8125, 8610760.9375, 8610820.3125, 8610862.5, 8611165.625, 8611376.5625, 8611543.75, 8611812.5, 8611954.6875, 8611981.25, 8612028.125, 8612139.0625, 8612292.1875, 8612434.375, 8612534.375, 8612551.5625, 8612556.25, 8612695.3125, 8612850.0, 8612900.0, 8613046.875, 8613390.625, 8613934.375, 8614112.5, 8614123.4375, 8614309.375, 8614465.625, 8614657.8125, 8614934.375, 8615059.375, 8615140.625, 8615143.75, 8615310.9375, 8615440.625, 8615479.6875, 8615593.75, 8615596.875, 8615651.5625, 8615771.875, 8615776.5625, 8615784.375, 8615976.5625, 8615979.6875, 8616081.25, 8616118.75, 8616346.875, 8616454.6875, 8616462.5, 8617076.5625, 8617235.9375, 8617242.1875, 8617259.375, 8617282.8125, 8617535.9375, 8617715.625, 8617870.3125, 8617928.125, 8617951.5625, 8617954.6875, 8618084.375, 8618239.0625, 8618285.9375, 8618287.5, 8619643.75, 8620910.9375, 8621507.8125, 8621828.125, 8622807.8125, 8623387.5, 8623506.25, 8623564.0625, 8623589.0625, 8623643.75, 8623851.5625, 8623942.1875, 8624025.0, 8624106.25, 8624159.375, 8624160.9375, 8624243.75, 8624310.9375, 8624312.5, 8624457.8125, 8624510.9375, 8624684.375, 8624834.375, 8624901.5625, 8624934.375, 8624956.25, 8624957.8125, 8625057.8125, 8625181.25, 8625253.125, 8625276.5625, 8625315.625, 8625437.5, 8625504.6875, 8625843.75, 8626071.875, 8626560.9375, 8627154.6875, 8627845.3125, 8627909.375, 8628589.0625, 8631904.6875, 8634373.4375, 8635718.75, 8638992.1875, 8639064.0625, 8639131.25, 8639229.6875, 8639239.0625, 8639295.3125, 8639509.375, 8639554.6875, 8639557.8125, 8639650.0, 8639660.9375, 8639681.25, 8639704.6875, 8639909.375, 8639957.8125, 8639968.75, 8639968.75, 8640026.5625, 8640054.6875, 8640056.25, 8640070.3125, 8640078.125, 8640096.875, 8640101.5625, 8640114.0625, 8640121.875, 8640148.4375, 8640346.875, 8640551.5625, 8640821.875, 8640853.125, 8640864.0625, 8640918.75, 8641004.6875, 8641307.8125, 8641389.0625, 8641429.6875, 8641551.5625, 8641689.0625, 8641690.625, 8641732.8125, 8641734.375, 8641784.375, 8641801.5625, 8641803.125, 8641809.375, 8641856.25, 8641871.875, 8641934.375, 8641950.0, 8641962.5, 8642003.125, 8642121.875, 8642135.9375, 8642307.8125, 8642381.25, 8642425.0, 8642548.4375, 8642706.25, 8642795.3125, 8642837.5, 8642857.8125, 8642873.4375, 8642931.25, 8643020.3125, 8643071.875, 8643275.0, 8643392.1875, 8643485.9375, 8643503.125, 8643525.0, 8643614.0625, 8643629.6875, 8643707.8125, 8643712.5, 8643714.0625, 8643725.0, 8643778.125, 8643850.0, 8643873.4375, 8643876.5625, 8643940.625, 8643942.1875, 8643993.75, 8644101.5625, 8644153.125, 8644196.875, 8644204.6875, 8644389.0625, 8644764.0625, 8645134.375, 8645153.125, 8645296.875, 8645540.625, 8645634.375, 8645687.5, 8645875.0, 8645884.375, 8646034.375, 8646095.3125, 8646123.4375, 8646132.8125, 8646295.3125, 8646315.625, 8646350.0, 8646592.1875, 8646828.125, 8646903.125, 8647328.125, 8647479.6875, 8647510.9375, 8647725.0, 8648160.9375, 8648223.4375, 8648718.75, 8648848.4375, 8648975.0, 8649401.5625, 8649403.125, 8649503.125, 8649507.8125, 8650290.625, 8650360.9375, 8650614.0625, 8650679.6875, 8650746.875, 8651656.25, 8652153.125, 8652303.125, 8652901.5625, 8653271.875, 8653300.0, 8653795.3125, 8654400.0, 8654492.1875, 8654612.5, 8654826.5625, 8654914.0625, 8655090.625, 8655118.75, 8655179.6875, 8655368.75, 8655568.75, 8655643.75, 8655645.3125, 8655673.4375, 8655712.5, 8655748.4375, 8655782.8125, 8655796.875, 8655801.5625, 8655803.125, 8655812.5, 8655820.3125, 8655832.8125, 8655851.5625, 8655860.9375, 8655929.6875, 8655951.5625, 8655965.625, 8656018.75, 8656031.25, 8656070.3125, 8656070.3125, 8656073.4375, 8656109.375, 8656131.25, 8656135.9375, 8656203.125, 8656214.0625, 8656226.5625, 8656228.125, 8656276.5625, 8656290.625, 8656321.875, 8656350.0, 8656351.5625, 8656359.375, 8656375.0, 8656378.125, 8656384.375, 8656392.1875, 8656435.9375, 8656439.0625, 8656443.75, 8656453.125, 8656454.6875, 8656476.5625, 8656479.6875, 8656501.5625, 8656506.25, 8656507.8125, 8656512.5, 8656534.375, 8656537.5, 8656543.75, 8656546.875, 8656548.4375, 8656548.4375, 8656556.25, 8656567.1875, 8656581.25, 8656589.0625, 8656592.1875, 8656600.0, 8656612.5, 8656614.0625, 8656626.5625, 8656639.0625, 8656650.0, 8656660.9375, 8656687.5, 8656696.875, 8656712.5, 8656717.1875, 8656728.125, 8656745.3125, 8656748.4375, 8656767.1875, 8656771.875, 8656798.4375, 8656820.3125, 8656846.875, 8656848.4375, 8656867.1875, 8656900.0, 8656904.6875, 8656956.25, 8656957.8125, 8656962.5, 8656981.25, 8656984.375, 8656984.375, 8656989.0625, 8657048.4375, 8657064.0625, 8657075.0, 8657075.0, 8657106.25, 8657143.75, 8657160.9375, 8657192.1875, 8657214.0625, 8657231.25, 8657248.4375, 8657257.8125, 8657284.375, 8657287.5, 8657318.75, 8657345.3125, 8657354.6875, 8657378.125, 8657395.3125, 8657442.1875, 8657456.25, 8657510.9375, 8657551.5625, 8657551.5625, 8657585.9375, 8657595.3125, 8657703.125, 8657706.25, 8657728.125, 8657831.25, 8657860.9375, 8657882.8125, 8657926.5625, 8658018.75, 8658021.875, 8658071.875, 8658106.25, 8658106.25, 8658134.375, 8658159.375, 8658181.25, 8658184.375, 8658189.0625, 8658196.875, 8658198.4375, 8658198.4375, 8658212.5, 8658221.875, 8658248.4375, 8658314.0625, 8658320.3125, 8658354.6875, 8658373.4375, 8658385.9375, 8658392.1875, 8658451.5625, 8658535.9375, 8658537.5, 8658607.8125, 8658612.5, 8658676.5625, 8658676.5625, 8658690.625, 8658726.5625, 8658771.875, 8658789.0625, 8658790.625, 8658798.4375, 8658871.875, 8658885.9375, 8658887.5, 8658921.875, 8658970.3125, 8659007.8125, 8659164.0625, 8659228.125, 8659240.625, 8659242.1875, 8659350.0, 8659401.5625, 8659603.125, 8659629.6875, 8659685.9375, 8659740.625, 8660125.0, 8660170.3125, 8660189.0625, 8660321.875, 8660390.625, 8660539.0625, 8660704.6875, 8661618.75, 8661754.6875, 8661856.25, 8662339.0625, 8662414.0625, 8662467.1875, 8662523.4375, 8662707.8125, 8662884.375, 8662892.1875, 8662892.1875, 8662904.6875, 8662917.1875, 8662957.8125, 8662959.375, 8662981.25, 8663121.875, 8663214.0625, 8663284.375, 8663306.25, 8663317.1875, 8663321.875, 8663554.6875, 8663564.0625, 8663570.3125, 8663717.1875, 8663740.625, 8663806.25, 8663814.0625, 8663842.1875, 8663871.875, 8663907.8125, 8663937.5, 8664001.5625, 8664003.125, 8664014.0625, 8664082.8125, 8664095.3125, 8664101.5625, 8664121.875, 8664139.0625, 8664146.875, 8664157.8125, 8664171.875, 8664200.0, 8664207.8125, 8664218.75, 8664273.4375, 8664375.0, 8664375.0, 8664392.1875, 8664393.75, 8664423.4375, 8664457.8125, 8664459.375, 8664464.0625, 8664467.1875, 8664473.4375, 8664481.25, 8664515.625, 8664532.8125, 8664534.375, 8664628.125, 8664645.3125, 8664682.8125, 8664703.125, 8664703.125, 8664740.625, 8664778.125, 8664825.0, 8664828.125, 8664842.1875, 8664856.25, 8664909.375, 8664992.1875, 8664995.3125, 8665000.0, 8665021.875, 8665059.375, 8665114.0625, 8665132.8125, 8665140.625, 8665179.6875, 8665262.5, 8665293.75, 8665337.5, 8665423.4375, 8665442.1875, 8665478.125, 8665503.125, 8665520.3125, 8665525.0, 8665551.5625, 8665556.25, 8665560.9375, 8665606.25, 8665609.375, 8665628.125, 8665657.8125, 8665732.8125, 8665751.5625, 8665787.5, 8665792.1875, 8665801.5625, 8665810.9375, 8665825.0, 8665828.125, 8665831.25, 8665837.5, 8665839.0625, 8665840.625, 8665845.3125, 8665848.4375, 8665851.5625, 8665860.9375, 8665878.125, 8665900.0, 8665906.25, 8665909.375, 8665935.9375, 8665948.4375, 8665948.4375, 8665962.5, 8665964.0625, 8665965.625, 8665973.4375, 8665976.5625, 8666023.4375, 8666035.9375, 8666062.5, 8666168.75, 8666179.6875, 8666184.375, 8666214.0625, 8666221.875, 8666254.6875, 8666262.5, 8666315.625, 8666328.125, 8666368.75, 8666412.5, 8666417.1875, 8666426.5625, 8666429.6875, 8666432.8125, 8666437.5, 8666440.625, 8666442.1875, 8666442.1875, 8666445.3125, 8666453.125, 8666465.625, 8666468.75, 8666496.875, 8666525.0, 8666534.375, 8666545.3125, 8666551.5625, 8666562.5, 8666582.8125, 8666598.4375, 8666603.125, 8666628.125, 8666632.8125, 8666651.5625, 8666701.5625, 8666718.75, 8666721.875, 8666721.875, 8666725.0, 8666743.75, 8666748.4375, 8666754.6875, 8666784.375, 8666807.8125, 8666834.375, 8666839.0625, 8666842.1875, 8666879.6875, 8666920.3125, 8666981.25, 8666985.9375, 8667065.625, 8667065.625, 8667198.4375, 8667201.5625, 8667248.4375, 8667307.8125, 8667567.1875, 8667756.25, 8671078.125, 8673585.9375], [20.022511874343287, 78.46195819126174, 16.4882159166492, 8.3806210227965, 22.655731992662222, 32.7894390233988, 19.992122006824662, 266.3683665547642, 37.971975006894496, 10.479829773658906, 69.07416418976844, 21.18869271607582, 211.06424103209028, 55.74804504420147, 47.64310231793246, 99.0119021928493, 57.440755801619154, 24.558960119136263, 14.00753740105926, 44.73793927065504, 11.942165738710592, 109.35471652356631, 9.115883724207817, 263.72326672636524, 83.74119530890762, 29.438975596620367, 27.96239366700433, 5.575536345699506, 16.55432859967531, 38.05377785665454, 9.517965289059601, 16.898529382988475, 54.086802474986435, 35.035608059892716, 65.39603760722872, 82.51020713418367, 27.12399671080831, 52.15975136926483, 57.12103036787454, 80.35785084335491, 109.76880160587254, 32.347070633102255, 43.79900953922901, 43.001750458424894, 121.31265616031341, 20.550560465803418, 29.589548108661166, 55.18851215289648, 12.639908647972959, 23.31526852045658, 110.79713508812934, 49.80955969255578, 74.46663672128123, 65.97957744967745, 30.964319590318677, 9.090094185837403, 11.307750037033616, 13.188127574503481, 103.05033685647533, 86.36100595153673, 69.76472482982308, 26.248142158722764, 6.722830443276934, 32.04335015358831, 11.671207498930526, 90.95072179078244, 29.672800258648127, 17.713107124165145, 34.586679194393454, 54.28498840904236, 24.533642372448206, 134.8606402561906, 7.0586037825888175, 12.59370809350691, 7.599111804700165, 51.108265454936515, 10.574791549462693, 11.627843477375071, 15.633348816853314, 15.456471071424817, 7.336746640943726, 27.238365400571542, 35.16720754684803, 27.32120359113039, 75.01455418802523, 20.246664143806306, 66.01829339722121, 15.354966589129722, 87.80361657410519, 6.567474863368416, 11.207176214421898, 8.827632784939132, 31.266688100606867, 9.047613354473457, 13.77135809467871, 18.392743426337734, 23.68670557234559, 33.47930157405227, 126.18021487531374, 50.54205196762329, 17.167018492401144, 15.107186850625812, 97.11346589073548, 41.092014407007355, 68.52099390329809, 16.809344891127335, 62.09238777153939, 9.081398378053724, 63.7713158004524, 85.05809963823498, 32.89963711215362, 6.28032830461624, 10.239256394195609, 84.25618049742822, 74.19170638390517, 20.625928059181273, 20.658027950333345, 18.65693408190566, 5.119125621013834, 57.32187041853211, 66.61046769498147, 47.714857279115904, 16.85534833666673, 30.463125340182778, 30.78821050302901, 36.34171075723825, 84.02079413917467, 32.60556176778896, 51.326718497721664, 59.72526238174085, 12.888860951436305, 27.429538697979616, 17.475360841031318, 13.235281348158273, 12.719850189976258, 59.231161340322366, 7.279705078240602, 15.017377153977073, 10.215416068621655, 28.180066957549535, 45.31755918376416, 11.73443880818547, 50.987476867285295, 106.18768313861088, 20.471439204566476, 9.42814307637266, 64.50711046847144, 61.38188200713779, 10.576642261587493, 92.61795480621976, 21.370780776051433, 49.483394211513556, 8.242707347926583, 71.82373137823643, 6.630914744746697, 22.852199563994994, 45.957270068714855, 13.845829546365433, 42.10796933532797, 11.40769711508683, 51.86207872478932, 31.56648948070083, 91.91959275006454, 12.664885596122, 66.64504524126973, 13.06646029002879, 17.741370755701762, 19.18840193266069, 14.656030646639255, 51.869781721360766, 76.54557599956274, 106.9839653635072, 59.070445608763066, 52.27023388168509, 43.571569797211836, 136.81591654951663, 29.209323085934436, 82.29972974349566, 8.513164735397329, 5.241121877643787, 27.070803994854682, 35.18780783075982, 26.667265187721117, 19.425925438821704, 108.90220960543692, 30.088448532601333, 77.88906824306977, 23.921220020601616, 19.211913475387956, 7.822580381785889, 79.3586105724116, 41.62256024242084, 16.731965433250142, 28.216761811564197, 48.745783913234575, 144.5740276545098, 108.36258700257906, 30.314148750774542, 94.16480367452326, 12.979280073092603, 6.579631609900123, 10.102405849151372, 11.755195295016643, 25.719677945668785, 14.000364461736062, 6.217876755306642, 8.347533360605425, 21.245435036259742, 60.14353331864014, 62.749066036269625, 113.09301290638473, 49.613538963498385, 29.933507753127913, 7.995748444568275, 5.679634866878216, 20.485660837258248, 33.78087741531586, 80.38556927958545, 5.699437769570864, 19.7616461800589, 58.41416928432731, 47.755126419368544, 42.451297655611675, 214.70992030312587, 100.9259985312806, 9.679266967878336, 64.36074394078048, 6.773238148044947, 92.54697354725297, 7.924083919981538, 8.174139286676452, 36.934720968456915, 92.93565708099888, 5.459213106317348, 6.090518003376605, 19.742142815854663, 19.803496751812347, 36.70184950543227, 27.913369717333985, 21.82117476984962, 5.125129303538489, 64.71258211212823, 51.78035997042599, 7.5458986114752795, 41.53948156003127, 20.10877436216406, 81.23550556147207, 44.467730693883595, 78.93483251147221, 16.393016965493437, 13.522621849028953, 16.318113515006356, 24.25150428408444, 69.69893291524218, 5.196451320385723, 76.60252628496639, 5.9382938197097594, 19.671632592932, 43.111224875237454, 18.504120477561496, 104.2334352893701, 28.70239436985423, 72.895920708373, 6.704828785895089, 5.219837702203113, 6.416149740213459, 14.852233857396223, 25.92187242281656, 65.9153753285124, 81.11533441930361, 12.678802646645423, 9.715760583687038, 21.40041399152576, 116.12946758162937, 30.737168083352625, 72.95060327693193, 13.737962967167903, 8.76030894904448, 33.01551040340161, 60.044895172334215, 61.58875211956814, 16.226434645763486, 23.83426538308519, 54.18206561852311, 89.00336895620045, 32.24268484031113, 45.52678353352946, 11.496322074438652, 24.851635574066623, 10.666430260161512, 28.83009946244532, 37.92612504920996, 64.4979782865023, 21.466793226517524, 24.554756428897164, 25.622392735159483, 24.396844781605573, 98.02827419184628, 26.786308615280976, 12.321624106663517, 113.4858965020849, 13.107961964124069, 36.84146593628902, 33.71028143484317, 66.78912462706931, 5.043443635812554, 8.739713791442218, 102.19965101063667, 28.457927867923743, 17.70694410850836, 91.59411401199382, 11.817212413967663, 48.606839961114005, 6.3391261383011575, 54.82502246307992, 17.651432804703134, 11.640600052826082, 7.519789398829752, 99.17138732617657, 13.809206423930991, 5.209696634599435, 15.868694825980015, 49.56455393489706, 30.926531444804127, 50.140494833105784, 5.662951429002725, 44.05934191867396, 73.3787904252479, 61.72763036999872, 17.115134546835463, 10.39227027119012, 80.00101962420351, 65.31042364056391, 27.795164747480506, 47.04273403509437, 66.44677002429101, 39.82537274026105, 20.35476358797571, 12.471259104111294, 7.408576455650318, 59.04874617140051, 22.3652553331502, 44.863944752851054, 21.8947129662407, 30.163699905595173, 54.01620857639806, 9.159586802774069, 31.885884725163745, 54.76488390261464, 61.65274490997159, 76.12458457793046, 71.02273743655472, 47.21128125138487, 17.928830725902557, 57.173070909319726, 29.944412868816933, 47.69260718562149, 8.582214411950929, 38.27455340358486, 82.84764914210425, 14.44497035508794, 11.558789004578951, 21.929138931610876, 89.60991337647583, 35.84241031802008, 204.76226066873977, 89.11441562687088, 130.3395686104191, 41.90226735708519, 10.906542596737676, 27.191609552778452, 21.4923950034853, 6.083369277025365, 6.3193460510604345, 14.628246260451183, 73.9696172451684, 17.762397614125252, 35.41798218095028, 25.170586494985553, 78.79448036693064, 51.267141346087925, 6.052412467583229, 32.01826712268636, 45.069618259105916, 15.2514955309576, 6.992247778645817, 19.575726083963634, 26.178468038852913, 77.0864992101433, 70.56467898256747, 11.69493680767102, 59.845506257241155, 100.93188083223943, 25.084907670659074, 126.3759158689946, 17.614106472524064, 5.283415402562491, 11.726348864065061, 15.937684882968709, 33.01024471737816, 89.19911701546434, 7.329428646013629, 6.442150934690001, 128.9611954871741, 28.776619523078274, 17.92625518166424, 21.607997670398746, 14.387768291917036, 11.675769363633774, 7.351455462336332, 96.88181622775627, 13.74798904204543, 42.598834825656375, 104.04222694035619, 15.224770233498267, 15.52351975734394, 13.24305784421485, 51.76451885824973, 105.9699117227311, 89.13416976708402, 50.937843533638976, 18.81046268119141, 41.21616643218582, 92.40899212100695, 10.103177935836095, 15.242602103859008, 6.525388242595328, 82.35647351597595, 25.57638234280975, 10.535537271580983, 9.902477610004253, 11.747016092488686, 32.19044362918392, 33.43058284987549, 65.65440270258851, 35.653220762959165, 94.25097240053438, 21.305260105431543, 117.94558302278995, 18.48559958881206, 33.68338644520263, 11.038231427715784, 22.026480489066678, 17.939826558337337, 20.600248528157383, 14.158737370029208, 25.575586794874734, 9.876561232705079, 109.54988620799597, 11.855428854986165, 28.68600718683858, 45.31653979728256, 77.65318337014287, 5.404245803455997, 58.560427337419554, 82.65493326315645, 103.55044672549838, 192.2821972775673, 10.89848038006748, 18.964574429713057, 29.806038730470995, 10.639721078601042, 5.853382900776596, 87.88459209376487, 12.94285883472737, 90.91315672034716, 20.448903597920555, 31.54120673670293, 18.622901973714768, 31.62538244245834, 62.76904861141467, 59.62812744885474, 17.709529335054118, 19.549264266853896, 83.0881184526006, 82.67009498700453, 92.2683732744003, 96.00571268258807, 140.65669593745503, 58.60451998251467, 84.46395574417433, 6.106678285131179, 149.13056839647655, 5.473634742081172, 24.147986864349097, 5.8269289365751, 8.72255948939593, 46.51641383412958, 19.388191873079457, 35.143627831512525, 57.81423672304632, 25.989744880207525, 28.428504908674135, 5.351621713915876, 5.962861094229473, 21.16035834130944, 12.824110878297748, 63.5973958510727, 6.162788945821083, 139.16828769483578, 84.80622533477869, 23.792076193277655, 17.556555539486602, 34.43751476099532, 34.635072595277705, 53.88655854257308, 5.070238021880104, 37.39054469587327, 77.12971427373495, 49.47619860833129, 35.99489257206772, 9.227171433471742, 16.175188402730786, 101.9419040258488, 66.31183869636145, 160.30263755650012, 82.75280998625504, 56.197924765544116, 84.60058145531889, 51.05360236113959, 8.461708011225136, 64.02861912700928, 24.469496950698353, 35.698299676169235, 120.1103699768657, 44.078659527620744, 8.810057758912102, 49.80594052248908, 60.59631290449353, 6.800482558910093, 268.1098685238079, 14.123992341776995, 18.01220299323216, 63.75094663903423, 65.77015331141875, 16.003008459238306, 45.261756339645, 59.46130310450647, 101.90110577656549, 30.95439635989883, 36.22264920973203, 5.834729884022688, 21.739053503173196, 53.55693304202912, 27.310036785786142, 15.864366766213793, 6.244712889682742, 17.81070332950543, 20.150770196073786, 19.897369770775857, 5.960205226816818, 12.235611204295266, 9.760620692507182, 10.301595669906776, 51.79188856497396, 6.661198377973035, 242.04154665597727, 70.62927489336107, 27.39152614941209, 11.669839141179615, 7.247172383162028, 8.027153950050806, 12.274792664247807, 140.94173602388585, 24.769180940669287, 78.95423260980462, 13.746892510407667, 50.58976834799397, 62.2663696708432, 6.589222718023403, 34.655756993546916, 17.661469204087584, 25.93921681623126, 41.59191780335905, 25.853977914090713, 5.674000142756098, 7.064583546260622, 16.96776611735011, 38.07506326797337, 27.894119883742512, 50.22611678101204, 22.331732811492387, 80.33384899338502, 7.453864284746731, 21.74563383116587, 43.087305145997405, 11.705325889875924, 5.13034697034523, 94.72400170687536, 31.022614469216222, 24.70087671786948, 6.535184344715774, 54.64346073784952, 36.984548369203935, 50.50368400665933, 14.312099607754234, 39.01142474040529, 164.2699329319394, 37.637661537495575, 62.92612670945739, 47.420675574632476, 59.56291967940864, 8.475325069754271, 7.362716790387138, 54.53587353043423, 9.981078612570554, 5.103098438423599, 111.31558223360828, 30.371546320290552, 135.83829324653965, 63.44798249026432, 61.4235018682927, 11.644028396324392, 103.32205672599883, 67.50617579043089, 47.201487949399336, 138.43995005943023, 33.25335871897556, 18.93688409068891, 29.39651862728524, 24.604227722655885, 66.2290043237164, 11.747061364281853, 72.00783120900645, 15.262258631532301, 194.59683263383826, 9.193671949082251, 14.333904367738203, 110.52435057405344, 7.311619127819326, 9.062697903193513, 35.14221476652692, 34.064974665538685, 6.156190680673775, 15.595471085729466, 79.67883339761102, 58.74197575479402, 52.70885025237563, 59.03949256911913, 49.539647018308855, 13.483474484748847, 5.678545366574067, 33.05941412032325, 7.894199239509674, 9.139836909229233, 28.44085182709599, 72.42039556863406, 25.130035427524316, 77.42351348495234, 65.07863038560795, 7.6947478396579445, 85.26641872112141, 84.02046636732139, 17.851658077672322, 20.39860524071606, 20.705606421693677, 66.6453174313181, 7.044914289454934, 43.39808779809452, 34.22995280665725, 12.424332790675388, 79.52004863645571, 94.01327923456759, 19.29874028501852, 74.16212097065974, 31.35814098147169, 5.633906664560219, 30.196130092479482, 22.200253961690052, 15.838434389782737, 54.607942114334925, 49.162520490964496, 13.847197336615569, 17.742091935065215, 10.676323571027847, 36.88807937831369, 9.400586886918543, 81.54197641079415, 74.4075568262742, 34.74525316187774, 8.194758611285586, 16.03138352423043, 9.891802850090373, 37.770360668300526, 25.91490705308511, 10.743371107468722, 112.78705592702772, 52.37018176569424, 99.20865864447491, 49.03835005170571, 65.55276210701308, 28.29871712038704, 75.47751554798941, 13.327607723974037, 8.95543719925393, 57.345375546623316, 57.43902798655975, 12.503607212202237, 18.351478482514604, 12.631854985408651, 54.16636366913231, 31.244070562358342, 49.815666605032014, 16.8438836565649, 7.05275400135345, 68.46363078734322, 84.41995261979798, 41.8908489682004, 76.14719386890287, 70.32925946034982, 7.258801841416832, 24.341853894424776, 64.7618300446885, 15.389104949024048, 29.24554684715484, 49.9566024019079, 45.16250262091817, 19.52449808983004, 11.352137036961382, 10.475539076057933, 108.72392502529694, 9.65589915384454, 10.153989750920672, 10.354607964319086, 25.54150175583259, 43.64366558125047, 8.104351799400517, 14.77639117722654, 20.891162055165875, 6.362054930396035, 51.19176794989235, 85.1224388289672, 45.59532921829374, 53.711103501495344, 9.035469961495378, 11.495398441954519, 23.638127990669535, 9.77640617647872, 10.090924376665647, 67.77431591196958, 44.64970585376821, 20.252866110965815, 6.917982172198404, 23.898553461192414, 71.79213771030712, 5.399232020239222, 99.50113926316546, 111.32782521083351, 13.366886557108431, 41.30704166040016, 17.584387282847256, 34.68783448152784, 25.042383812973835, 13.914366077272486, 99.62416896322168, 99.96587868132904, 27.583286254274224, 16.73111097122793, 6.1622996883229675, 12.387770139742969, 12.025053960250817, 29.23278296781624, 119.47019801700915, 31.856924793195535, 25.872075362000384, 40.61798042680902, 16.66448934095897, 43.989406796388835, 31.652349458941533, 62.517473557477764, 60.86169483950462, 59.114561262462814, 21.04670805195296, 72.20016327831769])
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);
([5593895.3125, 5848604.6875, 5991420.3125, 5991457.8125, 6033123.4375, 6252290.625, 6325731.25, 6723187.5, 7058406.25, 7075506.25, 7096498.4375, 7144004.6875, 7144901.5625, 7351868.75, 7464667.1875, 7485085.9375, 7838417.1875, 7894279.6875, 7895776.5625, 7937912.5, 7988645.3125, 7988650.0, 7988751.5625, 7989456.25, 8031082.8125, 8036070.3125, 8036851.5625, 8036903.125, 8043421.875, 8083312.5, 8085968.75, 8111518.75, 8331220.3125, 8332556.25, 8333475.0, 8335776.5625, 8338275.0, 8339084.375, 8339089.0625, 8339548.4375, 8343417.1875, 8346607.8125, 8349054.6875, 8352450.0, 8352742.1875, 8352748.4375, 8356268.75, 8356965.625, 8357035.9375, 8357364.0625, 8357523.4375, 8358268.75, 8358459.375, 8358818.75, 8359064.0625, 8362582.8125, 8365468.75, 8366692.1875, 8367253.125, 8367451.5625, 8369973.4375, 8374234.375, 8374492.1875, 8382121.875, 8383051.5625, 8394087.5, 8397353.125, 8398828.125, 8399879.6875, 8560350.0, 8562370.3125, 8568525.0, 8588934.375, 8589625.0, 8597015.625, 8597246.875, 8597325.0, 8597395.3125, 8597509.375, 8597515.625, 8597676.5625, 8597693.75, 8597773.4375, 8597793.75, 8597793.75, 8597929.6875, 8597971.875, 8598359.375, 8598571.875, 8598793.75, 8599925.0, 8599967.1875, 8600562.5, 8601332.8125, 8601501.5625, 8601715.625, 8602268.75, 8602387.5, 8602512.5, 8602610.9375, 8602789.0625, 8603528.125, 8603623.4375, 8603831.25, 8603845.3125, 8603954.6875, 8604090.625, 8604492.1875, 8604570.3125, 8604593.75, 8604860.9375, 8605042.1875, 8605314.0625, 8605314.0625, 8605412.5, 8605598.4375, 8605678.125, 8605745.3125, 8605973.4375, 8606071.875, 8606131.25, 8606209.375, 8606476.5625, 8607101.5625, 8607170.3125, 8607743.75, 8607984.375, 8608126.5625, 8608376.5625, 8608420.3125, 8608643.75, 8608732.8125, 8608751.5625, 8608839.0625, 8609157.8125, 8609217.1875, 8609248.4375, 8609382.8125, 8609473.4375, 8609559.375, 8609728.125, 8609860.9375, 8609925.0, 8610054.6875, 8610226.5625, 8610373.4375, 8610478.125, 8610514.0625, 8610521.875, 8610570.3125, 8610607.8125, 8610760.9375, 8610820.3125, 8610862.5, 8611165.625, 8611376.5625, 8611543.75, 8611812.5, 8611954.6875, 8611981.25, 8612028.125, 8612139.0625, 8612292.1875, 8612434.375, 8612534.375, 8612551.5625, 8612556.25, 8612695.3125, 8612850.0, 8612900.0, 8613046.875, 8613390.625, 8613934.375, 8614112.5, 8614123.4375, 8614309.375, 8614465.625, 8614657.8125, 8614934.375, 8615059.375, 8615140.625, 8615143.75, 8615310.9375, 8615440.625, 8615479.6875, 8615593.75, 8615596.875, 8615651.5625, 8615771.875, 8615776.5625, 8615784.375, 8615976.5625, 8615979.6875, 8616081.25, 8616118.75, 8616346.875, 8616454.6875, 8616462.5, 8617076.5625, 8617235.9375, 8617242.1875, 8617259.375, 8617282.8125, 8617535.9375, 8617715.625, 8617870.3125, 8617928.125, 8617951.5625, 8617954.6875, 8618084.375, 8618239.0625, 8618285.9375, 8618287.5, 8619643.75, 8620910.9375, 8621507.8125, 8621828.125, 8622807.8125, 8623387.5, 8623506.25, 8623564.0625, 8623589.0625, 8623643.75, 8623851.5625, 8623942.1875, 8624025.0, 8624106.25, 8624159.375, 8624160.9375, 8624243.75, 8624310.9375, 8624312.5, 8624457.8125, 8624510.9375, 8624684.375, 8624834.375, 8624901.5625, 8624934.375, 8624956.25, 8624957.8125, 8625057.8125, 8625181.25, 8625253.125, 8625276.5625, 8625315.625, 8625437.5, 8625504.6875, 8625843.75, 8626071.875, 8626560.9375, 8627154.6875, 8627845.3125, 8627909.375, 8628589.0625, 8631904.6875, 8634373.4375, 8635718.75, 8638992.1875, 8639064.0625, 8639131.25, 8639229.6875, 8639239.0625, 8639295.3125, 8639509.375, 8639554.6875, 8639557.8125, 8639650.0, 8639660.9375, 8639681.25, 8639704.6875, 8639909.375, 8639957.8125, 8639968.75, 8639968.75, 8640026.5625, 8640054.6875, 8640056.25, 8640070.3125, 8640078.125, 8640096.875, 8640101.5625, 8640114.0625, 8640121.875, 8640148.4375, 8640346.875, 8640551.5625, 8640821.875, 8640853.125, 8640864.0625, 8640918.75, 8641004.6875, 8641307.8125, 8641389.0625, 8641429.6875, 8641551.5625, 8641689.0625, 8641690.625, 8641732.8125, 8641734.375, 8641784.375, 8641801.5625, 8641803.125, 8641809.375, 8641856.25, 8641871.875, 8641934.375, 8641950.0, 8641962.5, 8642003.125, 8642121.875, 8642135.9375, 8642307.8125, 8642381.25, 8642425.0, 8642548.4375, 8642706.25, 8642795.3125, 8642837.5, 8642857.8125, 8642873.4375, 8642931.25, 8643020.3125, 8643071.875, 8643275.0, 8643392.1875, 8643485.9375, 8643503.125, 8643525.0, 8643614.0625, 8643629.6875, 8643707.8125, 8643712.5, 8643714.0625, 8643725.0, 8643778.125, 8643850.0, 8643873.4375, 8643876.5625, 8643940.625, 8643942.1875, 8643993.75, 8644101.5625, 8644153.125, 8644196.875, 8644204.6875, 8644389.0625, 8644764.0625, 8645134.375, 8645153.125, 8645296.875, 8645540.625, 8645634.375, 8645687.5, 8645875.0, 8645884.375, 8646034.375, 8646095.3125, 8646123.4375, 8646132.8125, 8646295.3125, 8646315.625, 8646350.0, 8646592.1875, 8646828.125, 8646903.125, 8647328.125, 8647479.6875, 8647510.9375, 8647725.0, 8648160.9375, 8648223.4375, 8648718.75, 8648848.4375, 8648975.0, 8649401.5625, 8649403.125, 8649503.125, 8649507.8125, 8650290.625, 8650360.9375, 8650614.0625, 8650679.6875, 8650746.875, 8651656.25, 8652153.125, 8652303.125, 8652901.5625, 8653271.875, 8653300.0, 8653795.3125, 8654400.0, 8654492.1875, 8654612.5, 8654826.5625, 8654914.0625, 8655090.625, 8655118.75, 8655179.6875, 8655368.75, 8655568.75, 8655643.75, 8655645.3125, 8655673.4375, 8655712.5, 8655748.4375, 8655782.8125, 8655796.875, 8655801.5625, 8655803.125, 8655812.5, 8655820.3125, 8655832.8125, 8655851.5625, 8655860.9375, 8655929.6875, 8655951.5625, 8655965.625, 8656018.75, 8656031.25, 8656070.3125, 8656070.3125, 8656073.4375, 8656109.375, 8656131.25, 8656135.9375, 8656203.125, 8656214.0625, 8656226.5625, 8656228.125, 8656276.5625, 8656290.625, 8656321.875, 8656350.0, 8656351.5625, 8656359.375, 8656375.0, 8656378.125, 8656384.375, 8656392.1875, 8656435.9375, 8656439.0625, 8656443.75, 8656453.125, 8656454.6875, 8656476.5625, 8656479.6875, 8656501.5625, 8656506.25, 8656507.8125, 8656512.5, 8656534.375, 8656537.5, 8656543.75, 8656546.875, 8656548.4375, 8656548.4375, 8656556.25, 8656567.1875, 8656581.25, 8656589.0625, 8656592.1875, 8656600.0, 8656612.5, 8656614.0625, 8656626.5625, 8656639.0625, 8656650.0, 8656660.9375, 8656687.5, 8656696.875, 8656712.5, 8656717.1875, 8656728.125, 8656745.3125, 8656748.4375, 8656767.1875, 8656771.875, 8656798.4375, 8656820.3125, 8656846.875, 8656848.4375, 8656867.1875, 8656900.0, 8656904.6875, 8656956.25, 8656957.8125, 8656962.5, 8656981.25, 8656984.375, 8656984.375, 8656989.0625, 8657048.4375, 8657064.0625, 8657075.0, 8657075.0, 8657106.25, 8657143.75, 8657160.9375, 8657192.1875, 8657214.0625, 8657231.25, 8657248.4375, 8657257.8125, 8657284.375, 8657287.5, 8657318.75, 8657345.3125, 8657354.6875, 8657378.125, 8657395.3125, 8657442.1875, 8657456.25, 8657510.9375, 8657551.5625, 8657551.5625, 8657585.9375, 8657595.3125, 8657703.125, 8657706.25, 8657728.125, 8657831.25, 8657860.9375, 8657882.8125, 8657926.5625, 8658018.75, 8658021.875, 8658071.875, 8658106.25, 8658106.25, 8658134.375, 8658159.375, 8658181.25, 8658184.375, 8658189.0625, 8658196.875, 8658198.4375, 8658198.4375, 8658212.5, 8658221.875, 8658248.4375, 8658314.0625, 8658320.3125, 8658354.6875, 8658373.4375, 8658385.9375, 8658392.1875, 8658451.5625, 8658535.9375, 8658537.5, 8658607.8125, 8658612.5, 8658676.5625, 8658676.5625, 8658690.625, 8658726.5625, 8658771.875, 8658789.0625, 8658790.625, 8658798.4375, 8658871.875, 8658885.9375, 8658887.5, 8658921.875, 8658970.3125, 8659007.8125, 8659164.0625, 8659228.125, 8659240.625, 8659242.1875, 8659350.0, 8659401.5625, 8659603.125, 8659629.6875, 8659685.9375, 8659740.625, 8660125.0, 8660170.3125, 8660189.0625, 8660321.875, 8660390.625, 8660539.0625, 8660704.6875, 8661618.75, 8661754.6875, 8661856.25, 8662339.0625, 8662414.0625, 8662467.1875, 8662523.4375, 8662707.8125, 8662884.375, 8662892.1875, 8662892.1875, 8662904.6875, 8662917.1875, 8662957.8125, 8662959.375, 8662981.25, 8663121.875, 8663214.0625, 8663284.375, 8663306.25, 8663317.1875, 8663321.875, 8663554.6875, 8663564.0625, 8663570.3125, 8663717.1875, 8663740.625, 8663806.25, 8663814.0625, 8663842.1875, 8663871.875, 8663907.8125, 8663937.5, 8664001.5625, 8664003.125, 8664014.0625, 8664082.8125, 8664095.3125, 8664101.5625, 8664121.875, 8664139.0625, 8664146.875, 8664157.8125, 8664171.875, 8664200.0, 8664207.8125, 8664218.75, 8664273.4375, 8664375.0, 8664375.0, 8664392.1875, 8664393.75, 8664423.4375, 8664457.8125, 8664459.375, 8664464.0625, 8664467.1875, 8664473.4375, 8664481.25, 8664515.625, 8664532.8125, 8664534.375, 8664628.125, 8664645.3125, 8664682.8125, 8664703.125, 8664703.125, 8664740.625, 8664778.125, 8664825.0, 8664828.125, 8664842.1875, 8664856.25, 8664909.375, 8664992.1875, 8664995.3125, 8665000.0, 8665021.875, 8665059.375, 8665114.0625, 8665132.8125, 8665140.625, 8665179.6875, 8665262.5, 8665293.75, 8665337.5, 8665423.4375, 8665442.1875, 8665478.125, 8665503.125, 8665520.3125, 8665525.0, 8665551.5625, 8665556.25, 8665560.9375, 8665606.25, 8665609.375, 8665628.125, 8665657.8125, 8665732.8125, 8665751.5625, 8665787.5, 8665792.1875, 8665801.5625, 8665810.9375, 8665825.0, 8665828.125, 8665831.25, 8665837.5, 8665839.0625, 8665840.625, 8665845.3125, 8665848.4375, 8665851.5625, 8665860.9375, 8665878.125, 8665900.0, 8665906.25, 8665909.375, 8665935.9375, 8665948.4375, 8665948.4375, 8665962.5, 8665964.0625, 8665965.625, 8665973.4375, 8665976.5625, 8666023.4375, 8666035.9375, 8666062.5, 8666168.75, 8666179.6875, 8666184.375, 8666214.0625, 8666221.875, 8666254.6875, 8666262.5, 8666315.625, 8666328.125, 8666368.75, 8666412.5, 8666417.1875, 8666426.5625, 8666429.6875, 8666432.8125, 8666437.5, 8666440.625, 8666442.1875, 8666442.1875, 8666445.3125, 8666453.125, 8666465.625, 8666468.75, 8666496.875, 8666525.0, 8666534.375, 8666545.3125, 8666551.5625, 8666562.5, 8666582.8125, 8666598.4375, 8666603.125, 8666628.125, 8666632.8125, 8666651.5625, 8666701.5625, 8666718.75, 8666721.875, 8666721.875, 8666725.0, 8666743.75, 8666748.4375, 8666754.6875, 8666784.375, 8666807.8125, 8666834.375, 8666839.0625, 8666842.1875, 8666879.6875, 8666920.3125, 8666981.25, 8666985.9375, 8667065.625, 8667065.625, 8667198.4375, 8667201.5625, 8667248.4375, 8667307.8125, 8667567.1875, 8667756.25, 8671078.125, 8673585.9375], [20.022511874343287, 78.46195819126174, 16.4882159166492, 8.3806210227965, 22.655731992662222, 32.7894390233988, 19.992122006824662, 266.3683665547642, 37.971975006894496, 10.479829773658906, 69.07416418976844, 21.18869271607582, 211.06424103209028, 55.74804504420147, 47.64310231793246, 99.0119021928493, 57.440755801619154, 24.558960119136263, 14.00753740105926, 44.73793927065504, 11.942165738710592, 109.35471652356631, 9.115883724207817, 263.72326672636524, 83.74119530890762, 29.438975596620367, 27.96239366700433, 5.575536345699506, 16.55432859967531, 38.05377785665454, 9.517965289059601, 16.898529382988475, 54.086802474986435, 35.035608059892716, 65.39603760722872, 82.51020713418367, 27.12399671080831, 52.15975136926483, 57.12103036787454, 80.35785084335491, 109.76880160587254, 32.347070633102255, 43.79900953922901, 43.001750458424894, 121.31265616031341, 20.550560465803418, 29.589548108661166, 55.18851215289648, 12.639908647972959, 23.31526852045658, 110.79713508812934, 49.80955969255578, 74.46663672128123, 65.97957744967745, 30.964319590318677, 9.090094185837403, 11.307750037033616, 13.188127574503481, 103.05033685647533, 86.36100595153673, 69.76472482982308, 26.248142158722764, 6.722830443276934, 32.04335015358831, 11.671207498930526, 90.95072179078244, 29.672800258648127, 17.713107124165145, 34.586679194393454, 54.28498840904236, 24.533642372448206, 134.8606402561906, 7.0586037825888175, 12.59370809350691, 7.599111804700165, 51.108265454936515, 10.574791549462693, 11.627843477375071, 15.633348816853314, 15.456471071424817, 7.336746640943726, 27.238365400571542, 35.16720754684803, 27.32120359113039, 75.01455418802523, 20.246664143806306, 66.01829339722121, 15.354966589129722, 87.80361657410519, 6.567474863368416, 11.207176214421898, 8.827632784939132, 31.266688100606867, 9.047613354473457, 13.77135809467871, 18.392743426337734, 23.68670557234559, 33.47930157405227, 126.18021487531374, 50.54205196762329, 17.167018492401144, 15.107186850625812, 97.11346589073548, 41.092014407007355, 68.52099390329809, 16.809344891127335, 62.09238777153939, 9.081398378053724, 63.7713158004524, 85.05809963823498, 32.89963711215362, 6.28032830461624, 10.239256394195609, 84.25618049742822, 74.19170638390517, 20.625928059181273, 20.658027950333345, 18.65693408190566, 5.119125621013834, 57.32187041853211, 66.61046769498147, 47.714857279115904, 16.85534833666673, 30.463125340182778, 30.78821050302901, 36.34171075723825, 84.02079413917467, 32.60556176778896, 51.326718497721664, 59.72526238174085, 12.888860951436305, 27.429538697979616, 17.475360841031318, 13.235281348158273, 12.719850189976258, 59.231161340322366, 7.279705078240602, 15.017377153977073, 10.215416068621655, 28.180066957549535, 45.31755918376416, 11.73443880818547, 50.987476867285295, 106.18768313861088, 20.471439204566476, 9.42814307637266, 64.50711046847144, 61.38188200713779, 10.576642261587493, 92.61795480621976, 21.370780776051433, 49.483394211513556, 8.242707347926583, 71.82373137823643, 6.630914744746697, 22.852199563994994, 45.957270068714855, 13.845829546365433, 42.10796933532797, 11.40769711508683, 51.86207872478932, 31.56648948070083, 91.91959275006454, 12.664885596122, 66.64504524126973, 13.06646029002879, 17.741370755701762, 19.18840193266069, 14.656030646639255, 51.869781721360766, 76.54557599956274, 106.9839653635072, 59.070445608763066, 52.27023388168509, 43.571569797211836, 136.81591654951663, 29.209323085934436, 82.29972974349566, 8.513164735397329, 5.241121877643787, 27.070803994854682, 35.18780783075982, 26.667265187721117, 19.425925438821704, 108.90220960543692, 30.088448532601333, 77.88906824306977, 23.921220020601616, 19.211913475387956, 7.822580381785889, 79.3586105724116, 41.62256024242084, 16.731965433250142, 28.216761811564197, 48.745783913234575, 144.5740276545098, 108.36258700257906, 30.314148750774542, 94.16480367452326, 12.979280073092603, 6.579631609900123, 10.102405849151372, 11.755195295016643, 25.719677945668785, 14.000364461736062, 6.217876755306642, 8.347533360605425, 21.245435036259742, 60.14353331864014, 62.749066036269625, 113.09301290638473, 49.613538963498385, 29.933507753127913, 7.995748444568275, 5.679634866878216, 20.485660837258248, 33.78087741531586, 80.38556927958545, 5.699437769570864, 19.7616461800589, 58.41416928432731, 47.755126419368544, 42.451297655611675, 214.70992030312587, 100.9259985312806, 9.679266967878336, 64.36074394078048, 6.773238148044947, 92.54697354725297, 7.924083919981538, 8.174139286676452, 36.934720968456915, 92.93565708099888, 5.459213106317348, 6.090518003376605, 19.742142815854663, 19.803496751812347, 36.70184950543227, 27.913369717333985, 21.82117476984962, 5.125129303538489, 64.71258211212823, 51.78035997042599, 7.5458986114752795, 41.53948156003127, 20.10877436216406, 81.23550556147207, 44.467730693883595, 78.93483251147221, 16.393016965493437, 13.522621849028953, 16.318113515006356, 24.25150428408444, 69.69893291524218, 5.196451320385723, 76.60252628496639, 5.9382938197097594, 19.671632592932, 43.111224875237454, 18.504120477561496, 104.2334352893701, 28.70239436985423, 72.895920708373, 6.704828785895089, 5.219837702203113, 6.416149740213459, 14.852233857396223, 25.92187242281656, 65.9153753285124, 81.11533441930361, 12.678802646645423, 9.715760583687038, 21.40041399152576, 116.12946758162937, 30.737168083352625, 72.95060327693193, 13.737962967167903, 8.76030894904448, 33.01551040340161, 60.044895172334215, 61.58875211956814, 16.226434645763486, 23.83426538308519, 54.18206561852311, 89.00336895620045, 32.24268484031113, 45.52678353352946, 11.496322074438652, 24.851635574066623, 10.666430260161512, 28.83009946244532, 37.92612504920996, 64.4979782865023, 21.466793226517524, 24.554756428897164, 25.622392735159483, 24.396844781605573, 98.02827419184628, 26.786308615280976, 12.321624106663517, 113.4858965020849, 13.107961964124069, 36.84146593628902, 33.71028143484317, 66.78912462706931, 5.043443635812554, 8.739713791442218, 102.19965101063667, 28.457927867923743, 17.70694410850836, 91.59411401199382, 11.817212413967663, 48.606839961114005, 6.3391261383011575, 54.82502246307992, 17.651432804703134, 11.640600052826082, 7.519789398829752, 99.17138732617657, 13.809206423930991, 5.209696634599435, 15.868694825980015, 49.56455393489706, 30.926531444804127, 50.140494833105784, 5.662951429002725, 44.05934191867396, 73.3787904252479, 61.72763036999872, 17.115134546835463, 10.39227027119012, 80.00101962420351, 65.31042364056391, 27.795164747480506, 47.04273403509437, 66.44677002429101, 39.82537274026105, 20.35476358797571, 12.471259104111294, 7.408576455650318, 59.04874617140051, 22.3652553331502, 44.863944752851054, 21.8947129662407, 30.163699905595173, 54.01620857639806, 9.159586802774069, 31.885884725163745, 54.76488390261464, 61.65274490997159, 76.12458457793046, 71.02273743655472, 47.21128125138487, 17.928830725902557, 57.173070909319726, 29.944412868816933, 47.69260718562149, 8.582214411950929, 38.27455340358486, 82.84764914210425, 14.44497035508794, 11.558789004578951, 21.929138931610876, 89.60991337647583, 35.84241031802008, 204.76226066873977, 89.11441562687088, 130.3395686104191, 41.90226735708519, 10.906542596737676, 27.191609552778452, 21.4923950034853, 6.083369277025365, 6.3193460510604345, 14.628246260451183, 73.9696172451684, 17.762397614125252, 35.41798218095028, 25.170586494985553, 78.79448036693064, 51.267141346087925, 6.052412467583229, 32.01826712268636, 45.069618259105916, 15.2514955309576, 6.992247778645817, 19.575726083963634, 26.178468038852913, 77.0864992101433, 70.56467898256747, 11.69493680767102, 59.845506257241155, 100.93188083223943, 25.084907670659074, 126.3759158689946, 17.614106472524064, 5.283415402562491, 11.726348864065061, 15.937684882968709, 33.01024471737816, 89.19911701546434, 7.329428646013629, 6.442150934690001, 128.9611954871741, 28.776619523078274, 17.92625518166424, 21.607997670398746, 14.387768291917036, 11.675769363633774, 7.351455462336332, 96.88181622775627, 13.74798904204543, 42.598834825656375, 104.04222694035619, 15.224770233498267, 15.52351975734394, 13.24305784421485, 51.76451885824973, 105.9699117227311, 89.13416976708402, 50.937843533638976, 18.81046268119141, 41.21616643218582, 92.40899212100695, 10.103177935836095, 15.242602103859008, 6.525388242595328, 82.35647351597595, 25.57638234280975, 10.535537271580983, 9.902477610004253, 11.747016092488686, 32.19044362918392, 33.43058284987549, 65.65440270258851, 35.653220762959165, 94.25097240053438, 21.305260105431543, 117.94558302278995, 18.48559958881206, 33.68338644520263, 11.038231427715784, 22.026480489066678, 17.939826558337337, 20.600248528157383, 14.158737370029208, 25.575586794874734, 9.876561232705079, 109.54988620799597, 11.855428854986165, 28.68600718683858, 45.31653979728256, 77.65318337014287, 5.404245803455997, 58.560427337419554, 82.65493326315645, 103.55044672549838, 192.2821972775673, 10.89848038006748, 18.964574429713057, 29.806038730470995, 10.639721078601042, 5.853382900776596, 87.88459209376487, 12.94285883472737, 90.91315672034716, 20.448903597920555, 31.54120673670293, 18.622901973714768, 31.62538244245834, 62.76904861141467, 59.62812744885474, 17.709529335054118, 19.549264266853896, 83.0881184526006, 82.67009498700453, 92.2683732744003, 96.00571268258807, 140.65669593745503, 58.60451998251467, 84.46395574417433, 6.106678285131179, 149.13056839647655, 5.473634742081172, 24.147986864349097, 5.8269289365751, 8.72255948939593, 46.51641383412958, 19.388191873079457, 35.143627831512525, 57.81423672304632, 25.989744880207525, 28.428504908674135, 5.351621713915876, 5.962861094229473, 21.16035834130944, 12.824110878297748, 63.5973958510727, 6.162788945821083, 139.16828769483578, 84.80622533477869, 23.792076193277655, 17.556555539486602, 34.43751476099532, 34.635072595277705, 53.88655854257308, 5.070238021880104, 37.39054469587327, 77.12971427373495, 49.47619860833129, 35.99489257206772, 9.227171433471742, 16.175188402730786, 101.9419040258488, 66.31183869636145, 160.30263755650012, 82.75280998625504, 56.197924765544116, 84.60058145531889, 51.05360236113959, 8.461708011225136, 64.02861912700928, 24.469496950698353, 35.698299676169235, 120.1103699768657, 44.078659527620744, 8.810057758912102, 49.80594052248908, 60.59631290449353, 6.800482558910093, 268.1098685238079, 14.123992341776995, 18.01220299323216, 63.75094663903423, 65.77015331141875, 16.003008459238306, 45.261756339645, 59.46130310450647, 101.90110577656549, 30.95439635989883, 36.22264920973203, 5.834729884022688, 21.739053503173196, 53.55693304202912, 27.310036785786142, 15.864366766213793, 6.244712889682742, 17.81070332950543, 20.150770196073786, 19.897369770775857, 5.960205226816818, 12.235611204295266, 9.760620692507182, 10.301595669906776, 51.79188856497396, 6.661198377973035, 242.04154665597727, 70.62927489336107, 27.39152614941209, 11.669839141179615, 7.247172383162028, 8.027153950050806, 12.274792664247807, 140.94173602388585, 24.769180940669287, 78.95423260980462, 13.746892510407667, 50.58976834799397, 62.2663696708432, 6.589222718023403, 34.655756993546916, 17.661469204087584, 25.93921681623126, 41.59191780335905, 25.853977914090713, 5.674000142756098, 7.064583546260622, 16.96776611735011, 38.07506326797337, 27.894119883742512, 50.22611678101204, 22.331732811492387, 80.33384899338502, 7.453864284746731, 21.74563383116587, 43.087305145997405, 11.705325889875924, 5.13034697034523, 94.72400170687536, 31.022614469216222, 24.70087671786948, 6.535184344715774, 54.64346073784952, 36.984548369203935, 50.50368400665933, 14.312099607754234, 39.01142474040529, 164.2699329319394, 37.637661537495575, 62.92612670945739, 47.420675574632476, 59.56291967940864, 8.475325069754271, 7.362716790387138, 54.53587353043423, 9.981078612570554, 5.103098438423599, 111.31558223360828, 30.371546320290552, 135.83829324653965, 63.44798249026432, 61.4235018682927, 11.644028396324392, 103.32205672599883, 67.50617579043089, 47.201487949399336, 138.43995005943023, 33.25335871897556, 18.93688409068891, 29.39651862728524, 24.604227722655885, 66.2290043237164, 11.747061364281853, 72.00783120900645, 15.262258631532301, 194.59683263383826, 9.193671949082251, 14.333904367738203, 110.52435057405344, 7.311619127819326, 9.062697903193513, 35.14221476652692, 34.064974665538685, 6.156190680673775, 15.595471085729466, 79.67883339761102, 58.74197575479402, 52.70885025237563, 59.03949256911913, 49.539647018308855, 13.483474484748847, 5.678545366574067, 33.05941412032325, 7.894199239509674, 9.139836909229233, 28.44085182709599, 72.42039556863406, 25.130035427524316, 77.42351348495234, 65.07863038560795, 7.6947478396579445, 85.26641872112141, 84.02046636732139, 17.851658077672322, 20.39860524071606, 20.705606421693677, 66.6453174313181, 7.044914289454934, 43.39808779809452, 34.22995280665725, 12.424332790675388, 79.52004863645571, 94.01327923456759, 19.29874028501852, 74.16212097065974, 31.35814098147169, 5.633906664560219, 30.196130092479482, 22.200253961690052, 15.838434389782737, 54.607942114334925, 49.162520490964496, 13.847197336615569, 17.742091935065215, 10.676323571027847, 36.88807937831369, 9.400586886918543, 81.54197641079415, 74.4075568262742, 34.74525316187774, 8.194758611285586, 16.03138352423043, 9.891802850090373, 37.770360668300526, 25.91490705308511, 10.743371107468722, 112.78705592702772, 52.37018176569424, 99.20865864447491, 49.03835005170571, 65.55276210701308, 28.29871712038704, 75.47751554798941, 13.327607723974037, 8.95543719925393, 57.345375546623316, 57.43902798655975, 12.503607212202237, 18.351478482514604, 12.631854985408651, 54.16636366913231, 31.244070562358342, 49.815666605032014, 16.8438836565649, 7.05275400135345, 68.46363078734322, 84.41995261979798, 41.8908489682004, 76.14719386890287, 70.32925946034982, 7.258801841416832, 24.341853894424776, 64.7618300446885, 15.389104949024048, 29.24554684715484, 49.9566024019079, 45.16250262091817, 19.52449808983004, 11.352137036961382, 10.475539076057933, 108.72392502529694, 9.65589915384454, 10.153989750920672, 10.354607964319086, 25.54150175583259, 43.64366558125047, 8.104351799400517, 14.77639117722654, 20.891162055165875, 6.362054930396035, 51.19176794989235, 85.1224388289672, 45.59532921829374, 53.711103501495344, 9.035469961495378, 11.495398441954519, 23.638127990669535, 9.77640617647872, 10.090924376665647, 67.77431591196958, 44.64970585376821, 20.252866110965815, 6.917982172198404, 23.898553461192414, 71.79213771030712, 5.399232020239222, 99.50113926316546, 111.32782521083351, 13.366886557108431, 41.30704166040016, 17.584387282847256, 34.68783448152784, 25.042383812973835, 13.914366077272486, 99.62416896322168, 99.96587868132904, 27.583286254274224, 16.73111097122793, 6.1622996883229675, 12.387770139742969, 12.025053960250817, 29.23278296781624, 119.47019801700915, 31.856924793195535, 25.872075362000384, 40.61798042680902, 16.66448934095897, 43.989406796388835, 31.652349458941533, 62.517473557477764, 60.86169483950462, 59.114561262462814, 21.04670805195296, 72.20016327831769])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)