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 = 44521
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);
([5470700.0, 6122045.3125, 6243959.375, 6290701.5625, 6297417.1875, 6312395.3125, 6315159.375, 6317226.5625, 6317284.375, 6318421.875, 6365495.3125, 6365546.875, 6368198.4375, 6372684.375, 6375050.0, 6376064.0625, 6405745.3125, 6413112.5, 6415925.0, 6417312.5, 6434789.0625, 6461846.875, 6544517.1875, 6557629.6875, 6565242.1875, 6565315.625, 6568340.625, 6576395.3125, 6610321.875, 6610728.125, 6611550.0, 6612453.125, 6613289.0625, 6614284.375, 6614582.8125, 6615031.25, 6615117.1875, 6616489.0625, 6638365.625, 6666957.8125, 6671709.375, 6692260.9375, 6727423.4375, 6747229.6875, 6754425.0, 6770390.625, 6782639.0625, 6786725.0, 6790895.3125, 6790925.0, 6791337.5, 6792334.375, 6792496.875, 6792654.6875, 6793148.4375, 6793245.3125, 6793326.5625, 6794039.0625, 6794529.6875, 6795168.75, 6797026.5625, 6809607.8125, 6841125.0, 6841576.5625, 6842840.625, 6843528.125, 6844050.0, 6844876.5625, 6845046.875, 6845260.9375, 6845435.9375, 6931915.625, 6932621.875, 6939903.125, 6948675.0, 6948954.6875, 6976082.8125, 6977854.6875, 6995010.9375, 6995221.875, 6999707.8125, 7010454.6875, 7024406.25, 7051040.625, 7056859.375, 7087010.9375, 7152462.5, 7191401.5625, 7203760.9375, 7206790.625, 7227307.8125, 7268150.0, 7274557.8125, 7303523.4375, 7483157.8125, 7515953.125, 7532070.3125, 7585101.5625, 7592435.9375, 7617507.8125, 7640865.625, 8158739.0625, 8341992.1875, 8742915.625, 8745429.6875, 8896539.0625, 8897218.75, 8900353.125, 8922768.75, 8952323.4375, 9016885.9375, 9053489.0625, 9056553.125, 9071548.4375, 9081501.5625, 9083954.6875, 9084728.125, 9093759.375, 9184885.9375, 9185790.625, 9186853.125, 9187056.25, 9187242.1875, 9187826.5625, 9187940.625, 9188356.25, 9188840.625, 9189709.375, 9190681.25, 9191962.5, 9211017.1875, 9211059.375, 9212698.4375, 9213515.625, 9214260.9375, 9215653.125, 9216051.5625, 9216325.0, 9216329.6875, 9218117.1875, 9219581.25, 9224764.0625, 9225714.0625, 9231679.6875, 9239634.375, 9241787.5, 9249125.0, 9250290.625, 9251634.375, 9252456.25, 9255898.4375, 9283743.75, 9344678.125, 9345293.75, 9345715.625, 9351218.75, 9358037.5, 9366471.875, 9368329.6875, 9368656.25, 9368768.75, 9368971.875, 9370218.75, 9371384.375, 9372092.1875, 9372323.4375, 9378300.0, 9378612.5, 9380526.5625, 9387962.5, 9388795.3125, 9391978.125, 9392406.25, 9393000.0, 9393389.0625, 9401007.8125, 9402292.1875, 9404040.625, 9411431.25, 9411435.9375, 9423295.3125, 9424525.0, 9447748.4375, 9448029.6875, 9449560.9375, 9449623.4375, 9449703.125, 9450314.0625, 9452034.375, 9452543.75, 9452779.6875, 9453740.625, 9455087.5, 9456306.25, 9456981.25, 9456992.1875, 9495570.3125, 9496193.75, 9497059.375, 9497332.8125, 9497600.0, 9497871.875, 9499110.9375, 9499148.4375, 9500146.875, 9507839.0625, 9508890.625, 9510567.1875, 9511082.8125, 9511800.0, 9522014.0625, 9525356.25, 9526981.25, 9529025.0, 9529437.5, 9538645.3125, 9541935.9375, 9542082.8125, 9543035.9375, 9543720.3125, 9610778.125, 9619928.125, 9638068.75, 9715906.25, 9717585.9375, 9719315.625, 9741928.125, 9742382.8125, 9743175.0, 9743557.8125, 9744164.0625, 9744270.3125, 9747600.0, 9747651.5625, 9749429.6875, 9751496.875, 9767779.6875, 9773337.5, 9776160.9375, 9786034.375, 9786995.3125, 9791762.5, 9796870.3125, 9796893.75, 9806301.5625, 9811043.75, 9890292.1875, 9906764.0625, 9925570.3125, 9931945.3125, 9965332.8125, 9969200.0, 9969675.0, 9969748.4375, 9973417.1875, 9976837.5, 9988881.25, 10009853.125, 10011873.4375, 10013937.5, 10023450.0, 10023475.0, 10023487.5, 10038245.3125, 10100821.875, 10103193.75, 10104520.3125, 10105098.4375, 10115637.5, 10132003.125, 10133939.0625, 10141028.125, 10171753.125, 10204665.625, 10235901.5625, 10257421.875, 10263312.5, 10297973.4375, 10316348.4375, 10421642.1875, 10515720.3125, 10520773.4375, 10521157.8125, 10557134.375, 10616289.0625, 10664932.8125, 11031510.9375, 11081223.4375, 11087146.875, 11148254.6875, 11150075.0, 11154546.875, 11157140.625, 11157390.625, 11158046.875, 11159251.5625, 11163929.6875, 11172446.875, 11178334.375, 11181218.75, 11181789.0625, 11181923.4375, 11182198.4375, 11185790.625, 11186618.75, 11186893.75, 11187134.375, 11187229.6875, 11187509.375, 11187545.3125, 11187996.875, 11188053.125, 11188093.75, 11188289.0625, 11188298.4375, 11188331.25, 11188475.0, 11188684.375, 11189248.4375, 11189285.9375, 11189503.125, 11189548.4375, 11189592.1875, 11189967.1875, 11189973.4375, 11190028.125, 11190104.6875, 11190146.875, 11190362.5, 11190379.6875, 11190412.5, 11190440.625, 11190889.0625, 11190901.5625, 11190912.5, 11190926.5625, 11190970.3125, 11190984.375, 11191029.6875, 11191092.1875, 11191110.9375, 11191132.8125, 11191142.1875, 11191226.5625, 11191315.625, 11191331.25, 11191434.375, 11191521.875, 11191575.0, 11191585.9375, 11191670.3125, 11191746.875, 11191912.5, 11192279.6875, 11192371.875, 11192506.25, 11192553.125, 11192556.25, 11192593.75, 11192596.875, 11192756.25, 11192851.5625, 11192982.8125, 11193021.875, 11193568.75, 11193670.3125, 11193929.6875, 11195054.6875, 11195354.6875, 11195404.6875, 11195704.6875, 11195765.625, 11195878.125, 11196057.8125, 11196545.3125, 11197056.25, 11197087.5, 11197525.0, 11197757.8125, 11197810.9375, 11198140.625, 11198479.6875, 11198828.125, 11198859.375, 11199081.25, 11199135.9375, 11199156.25, 11199279.6875, 11199400.0, 11199479.6875, 11199698.4375, 11199834.375, 11199957.8125, 11200068.75, 11200262.5, 11200337.5, 11200340.625, 11200615.625, 11200709.375, 11200728.125, 11200754.6875, 11200759.375, 11200873.4375, 11201026.5625, 11201135.9375, 11201300.0, 11201323.4375, 11201425.0, 11201487.5, 11201720.3125, 11201721.875, 11201723.4375, 11201868.75, 11202296.875, 11202446.875, 11202575.0, 11202614.0625, 11202729.6875, 11203282.8125, 11203473.4375, 11203918.75, 11203931.25, 11203984.375, 11204048.4375, 11204375.0, 11204503.125, 11204721.875, 11204970.3125, 11204973.4375, 11205012.5, 11205246.875, 11205748.4375, 11205829.6875, 11205895.3125, 11205915.625, 11206426.5625, 11206657.8125, 11207381.25, 11207428.125, 11207442.1875, 11207800.0, 11208667.1875, 11208676.5625, 11209234.375, 11209435.9375, 11209575.0, 11209787.5, 11210067.1875, 11210945.3125, 11211142.1875, 11211165.625, 11211206.25, 11211392.1875, 11211953.125, 11212381.25, 11212504.6875, 11212795.3125, 11213148.4375, 11213168.75, 11213350.0, 11213417.1875, 11213596.875, 11213603.125, 11214114.0625, 11214151.5625, 11214151.5625, 11214232.8125, 11214279.6875, 11214700.0, 11214840.625, 11214898.4375, 11214903.125, 11214906.25, 11215139.0625, 11215178.125, 11215231.25, 11215243.75, 11215276.5625, 11215678.125, 11215817.1875, 11215976.5625, 11215985.9375, 11216157.8125, 11216309.375, 11216551.5625, 11217090.625, 11217496.875, 11217503.125, 11217675.0, 11217793.75, 11217915.625, 11217970.3125, 11218021.875, 11218479.6875, 11218504.6875, 11218801.5625, 11218946.875, 11219442.1875, 11219942.1875, 11219964.0625, 11220021.875, 11220178.125, 11220557.8125, 11220560.9375, 11220590.625, 11220715.625, 11220875.0, 11220925.0, 11220970.3125, 11221028.125, 11221162.5, 11221514.0625, 11221576.5625, 11221579.6875, 11221593.75, 11221625.0, 11221628.125, 11221654.6875, 11221712.5, 11221957.8125, 11222310.9375, 11222510.9375, 11222565.625, 11222570.3125, 11222685.9375, 11222725.0, 11223121.875, 11223135.9375, 11223337.5, 11223454.6875, 11223609.375, 11223725.0, 11223840.625, 11224075.0, 11224110.9375, 11224146.875, 11224329.6875, 11224376.5625, 11224506.25, 11224570.3125, 11224675.0, 11224732.8125, 11224953.125, 11224957.8125, 11225084.375, 11225321.875, 11225578.125, 11225617.1875, 11225681.25, 11225715.625, 11225753.125, 11225798.4375, 11226046.875, 11226140.625, 11226384.375, 11226423.4375, 11226467.1875, 11226467.1875, 11226482.8125, 11226514.0625, 11226540.625, 11226557.8125, 11226567.1875, 11226576.5625, 11226651.5625, 11226818.75, 11226828.125, 11226860.9375, 11226870.3125, 11226992.1875, 11227156.25, 11227181.25, 11227203.125, 11227292.1875, 11227292.1875, 11227478.125, 11227657.8125, 11227671.875, 11227723.4375, 11227740.625, 11227809.375, 11227812.5, 11227842.1875, 11227895.3125, 11227901.5625, 11227907.8125, 11228212.5, 11228301.5625, 11228368.75, 11228375.0, 11228400.0, 11228468.75, 11228610.9375, 11228725.0, 11228742.1875, 11228815.625, 11228875.0, 11228882.8125, 11229112.5, 11229334.375, 11229376.5625, 11229526.5625, 11229543.75, 11229554.6875, 11229739.0625, 11229798.4375, 11229803.125, 11229870.3125, 11229929.6875, 11229985.9375, 11230042.1875, 11230084.375, 11230106.25, 11230178.125, 11230264.0625, 11230359.375, 11230414.0625, 11230500.0, 11230635.9375, 11230651.5625, 11230670.3125, 11230912.5, 11231665.625, 11231750.0, 11231854.6875, 11231867.1875, 11232121.875, 11232148.4375, 11232245.3125, 11232439.0625, 11232779.6875, 11232881.25, 11233065.625, 11233134.375, 11233137.5, 11233325.0, 11233395.3125, 11233398.4375, 11233868.75, 11234004.6875, 11234018.75, 11234459.375, 11234584.375, 11234676.5625, 11235034.375, 11235103.125, 11235625.0, 11235676.5625, 11235748.4375, 11235876.5625, 11235884.375, 11236228.125, 11236296.875, 11236773.4375, 11237193.75, 11237264.0625, 11237368.75, 11237385.9375, 11237579.6875, 11237826.5625, 11238039.0625, 11238207.8125, 11238731.25, 11238868.75, 11239037.5, 11239109.375, 11239160.9375, 11239621.875, 11240298.4375, 11240473.4375, 11240667.1875, 11241101.5625, 11241221.875, 11241307.8125, 11241346.875, 11241737.5, 11241751.5625, 11241756.25, 11241821.875, 11241904.6875, 11241917.1875, 11241943.75, 11242012.5, 11242014.0625, 11242056.25, 11242626.5625, 11242646.875, 11242710.9375, 11243243.75, 11243496.875, 11243656.25, 11243726.5625, 11244084.375, 11244337.5, 11245167.1875, 11245357.8125, 11245587.5, 11245631.25, 11245703.125, 11245731.25, 11245851.5625, 11245975.0, 11246042.1875, 11246170.3125, 11246182.8125, 11246904.6875, 11247129.6875, 11247309.375, 11247362.5, 11247471.875, 11247734.375, 11247900.0, 11248231.25, 11248239.0625, 11248315.625, 11248387.5, 11248528.125, 11248682.8125, 11248753.125, 11248785.9375, 11248892.1875, 11249142.1875, 11249175.0, 11249695.3125, 11249734.375, 11250970.3125, 11251001.5625, 11251384.375, 11251653.125, 11251717.1875, 11251846.875, 11252343.75, 11252778.125, 11253340.625, 11254192.1875, 11254642.1875, 11254710.9375, 11254740.625, 11254776.5625, 11254784.375, 11255078.125, 11255620.3125, 11255815.625, 11255928.125, 11256446.875, 11256446.875, 11256614.0625, 11257279.6875, 11257450.0, 11258204.6875, 11258609.375, 11258634.375, 11259451.5625, 11259792.1875, 11260662.5, 11260715.625, 11260859.375, 11261665.625, 11262460.9375, 11262564.0625, 11262565.625, 11262625.0, 11264050.0, 11264137.5, 11265189.0625, 11265687.5, 11265751.5625, 11265840.625, 11265943.75, 11266532.8125, 11267440.625, 11267560.9375, 11267832.8125, 11268154.6875, 11269114.0625, 11269442.1875, 11269500.0, 11269903.125, 11269984.375, 11270670.3125, 11271157.8125, 11271178.125, 11271323.4375, 11271457.8125, 11271507.8125, 11271621.875, 11271768.75, 11271820.3125, 11271985.9375, 11272401.5625, 11272409.375, 11272700.0, 11272768.75, 11272804.6875, 11272954.6875, 11273239.0625, 11273246.875, 11273371.875, 11273478.125, 11273490.625, 11273871.875, 11274337.5, 11274381.25, 11274671.875, 11274860.9375, 11274910.9375, 11275226.5625, 11275287.5, 11275298.4375, 11275357.8125, 11275526.5625, 11275825.0, 11276135.9375, 11276182.8125, 11276209.375, 11276214.0625, 11276262.5, 11276278.125, 11276431.25, 11276592.1875, 11276671.875, 11276798.4375, 11276882.8125, 11277109.375, 11277114.0625, 11277500.0, 11277615.625, 11277729.6875, 11277864.0625, 11278112.5, 11278159.375, 11278257.8125, 11278439.0625, 11278887.5, 11278912.5, 11278935.9375, 11278951.5625, 11279217.1875, 11279303.125, 11279500.0, 11279790.625, 11280242.1875, 11280390.625, 11280829.6875, 11280932.8125, 11280975.0, 11281148.4375, 11282146.875, 11282423.4375, 11282645.3125, 11283068.75, 11283117.1875, 11283446.875, 11283675.0, 11285146.875, 11285829.6875, 11286868.75, 11286884.375, 11287160.9375, 11287489.0625, 11288742.1875, 11289150.0, 11289871.875, 11290551.5625, 11290565.625, 11290871.875, 11292012.5, 11293326.5625, 11293557.8125, 11294004.6875, 11294284.375, 11297067.1875, 11298284.375, 11301239.0625, 11307773.4375], [19.2154295523539, 89.55258854213393, 14.370355298661053, 57.20399228797166, 44.811780504814074, 26.75388049298675, 6.530170367706963, 46.51528292778501, 5.291261754564691, 11.241879054981549, 24.51955003314039, 5.4735318475161385, 16.98705148080214, 6.854651163154591, 134.54444194779637, 108.56779278433693, 29.617489241934027, 51.53468854959295, 45.763343383599434, 29.438834060151855, 8.984525156225294, 11.413458231150608, 79.11277376253217, 29.332268207312964, 7.754292493196194, 5.1596478858584565, 65.75759256528346, 29.886755584893912, 26.5584010269535, 23.863277590325925, 14.547172927251028, 60.30310941298213, 19.09707671046739, 156.3456596865374, 17.16082850848504, 46.84886998523748, 22.061472761700315, 21.5868778827754, 9.24097056919769, 23.262784910507882, 118.25907569736367, 61.11357374615261, 38.57253171293951, 47.46288157118514, 34.99858159090979, 45.693231652771615, 9.251113344025454, 8.059526900380703, 90.49611075964678, 20.520631561015886, 11.720049670994765, 29.740928249506485, 107.58167149662388, 6.40105327770997, 28.50712426756071, 13.07560101312344, 68.26889573018359, 17.74258765455825, 56.945125949754264, 72.29505848984574, 6.020204843413952, 44.54701547359016, 40.36894174918001, 17.067673224056758, 9.605313045803141, 18.59978060841521, 57.41665724064427, 217.70902154594663, 60.42647272904755, 64.01513099818922, 5.657883631660489, 69.18948929366888, 12.09263699146219, 13.107291736562065, 39.879337497829894, 22.24404132702483, 23.154982776707463, 123.23317597273478, 53.05712160521556, 25.354345409682807, 6.980673241152645, 30.94223972389825, 10.214020951023041, 142.33953166532407, 297.29521644971095, 10.760884240449247, 23.924188692267542, 145.2495245410926, 7.913224832692245, 70.96656852421646, 30.52159985969437, 99.8772153032514, 49.93174090042638, 30.353901470667097, 9.26361765387983, 21.240483298665268, 9.544389337823151, 33.41682334052246, 124.01046154983763, 31.81760139138788, 97.40791875541008, 31.90252501810467, 26.154246510660762, 5.140815620989394, 17.177225782037713, 64.50741966314932, 35.957545738789264, 104.5144784399327, 26.64889928125563, 59.48767625255684, 16.942852235320416, 134.40837950624547, 22.053034040793786, 6.23291908819242, 5.911287509029278, 110.56842971665007, 7.677522869714147, 6.284719157619654, 42.683435128956695, 8.32042025593095, 9.928282408649084, 19.722045606635636, 28.987548108221837, 22.541174969365414, 100.96387202082396, 5.803734727703886, 19.65109117609528, 24.842907083337348, 30.680043183427458, 8.658986238020574, 26.621158222116563, 35.042213976536594, 56.984986104709336, 22.89199804998787, 12.055277195586118, 67.62044480693226, 16.47506936148762, 31.985816272631364, 20.16372603464425, 25.89806132220823, 28.880419122242554, 15.46104680811437, 21.045317789123715, 45.491566144160146, 85.2147434472077, 17.309245418570566, 8.078382959682134, 8.195185905875588, 20.622677176985167, 6.1076938727813745, 137.97297094019163, 40.72741328583843, 105.8852169166256, 59.00638138940108, 37.84499218295711, 18.052755338778198, 12.777847454904615, 18.981400770137018, 166.2523480868819, 119.7609513969797, 114.13651494861878, 16.652970221240064, 46.62280759256032, 195.60053890605388, 46.59376171522019, 84.21447092167637, 35.67032040601055, 54.621348313428555, 22.53155975251148, 43.436918197163806, 15.338022649426478, 38.304689734371735, 15.100770425629005, 44.81912929448974, 6.320434263713759, 28.76550382602135, 107.88435606427258, 8.450202504917593, 84.97320765353425, 37.56270282959764, 58.030205410898795, 5.7709920107508434, 32.06595252128765, 17.993235980718183, 17.6253432354124, 41.00399468709798, 59.93992228451643, 84.16686319144326, 10.351134411316657, 13.5293570573456, 49.50857954099376, 34.1751830579345, 63.41813705918761, 10.368743511179968, 18.738693240509853, 6.657544770175715, 52.444248361771116, 49.51492145254019, 24.429689187005337, 12.706705068052425, 53.697597530578605, 8.236528576630326, 157.20851863249132, 34.22233017568783, 24.64701976296206, 47.981146511860544, 56.64239394103669, 5.7529193433668295, 21.83882126418456, 19.57399568647594, 32.83299838602224, 298.4792889812215, 98.62888764356468, 26.47639052807266, 55.785696132754126, 19.240462877606962, 18.30042911850452, 55.666689248636914, 6.032081601307642, 103.43258385604304, 7.426222274868422, 175.87626436116688, 5.208550669803976, 87.66195913679998, 6.970929251076096, 101.57004136135618, 51.365073943841615, 60.9338213153302, 52.62441405820311, 12.093841119121558, 79.93321541879693, 131.19417342765024, 8.139494005406712, 5.821175516598972, 8.455144136445806, 70.4630275673035, 34.323304711144594, 84.58534687972843, 108.64944103920774, 13.681570178809102, 62.22614937523675, 12.365094712858127, 28.492406438982332, 11.695273511376318, 52.46457819806943, 80.60044291018835, 66.2486926663718, 74.93155719876815, 5.887282941632003, 18.54482833822992, 10.659805712491513, 106.2018419397337, 78.7222965759835, 5.176864716447867, 47.39263194185077, 58.657984790067154, 12.137771994305721, 5.6537688156527866, 57.06383264487278, 15.620701400616198, 48.72667263087499, 6.868133893677933, 7.624943522082688, 6.511652217974223, 28.12847202639719, 6.3308568389160085, 85.81788067101157, 39.5355884853977, 25.799927371278617, 5.635185781668891, 19.413694709855502, 103.66360819193886, 39.27383368194472, 13.696639081188916, 44.29238541410025, 110.26364335749889, 20.176164818304052, 11.964093289182197, 48.53086118336387, 49.204634448661025, 6.970078167026922, 8.935250347429678, 47.064436161810555, 98.01106467781845, 13.380820703223215, 35.649780551146065, 8.416738377800952, 6.425517229849981, 85.7151444578088, 12.068191739864806, 10.415593125635109, 8.972398263117928, 42.37344961171816, 58.2661475358052, 23.319547183968105, 73.68292085277879, 80.81713775262017, 38.13017052214444, 7.693987769674597, 16.377388457263514, 86.45912336994074, 93.21326374404448, 52.01944740977451, 18.291413416405508, 7.514705503132166, 40.84015074842837, 17.22129104662929, 6.398440824001257, 20.405598255080214, 25.33041735696217, 13.445138437343513, 92.95164637472942, 41.36511519869156, 65.44956925553036, 55.155084595589955, 30.00283269197329, 19.612254510784922, 32.83023133034581, 44.758132906000625, 15.308796375793193, 12.413608308115478, 10.68845631053543, 9.245130215340437, 16.671255064750632, 37.374326254738456, 5.777629518547999, 35.04937588090091, 17.980773221899423, 33.13911708660834, 15.058198397742165, 76.7816122129089, 6.377286346592362, 44.489582365049124, 45.18491444488918, 10.724624285581632, 14.070434476963813, 45.13512991270552, 12.609869601272743, 10.595083751588458, 20.89065575264003, 48.26806553866681, 36.10995743514956, 7.7460302051703, 7.505432818216993, 25.27096680329298, 5.128056347923415, 8.081620682470584, 16.787434254518708, 12.548672309875796, 7.309578165500295, 29.020397725107284, 15.819777974452645, 21.842558275844425, 42.749059836789236, 56.04955688135264, 18.431463067382307, 41.86821160239721, 24.012749600134843, 81.63171093430608, 27.620072308716992, 9.997263199271416, 7.754638245339949, 60.24650286862469, 23.835065228706924, 7.401011322532126, 15.703710682886056, 7.13657384749423, 13.657442972863167, 13.416787128507428, 31.59162517362728, 39.78198164977961, 15.291569046090522, 21.48855460426641, 146.21183152002635, 59.78011245708154, 70.81183821481665, 141.7277492793319, 26.403224294035706, 39.51417239946688, 30.897482814382986, 17.552609536530802, 9.58766721154342, 77.56996682126315, 17.147744434704386, 31.058917482281707, 11.290727540660633, 20.224815672530422, 65.82104401212247, 72.92231591405691, 6.585060122686211, 11.125336338174773, 5.866850039796829, 16.118845809924984, 63.95025140890468, 67.25577173679876, 13.384798335304454, 82.94665171031781, 35.619306463418056, 23.205619874504407, 74.20522964493901, 7.340666706822463, 14.749773076930264, 9.271464809005398, 46.48912410617532, 7.948426370465435, 23.756361712600018, 6.088395514391004, 6.302727603136325, 71.53022670204054, 35.021021869095875, 17.623555629125867, 14.411791763744088, 13.628827912187504, 5.404040694168145, 17.53156105658011, 19.99823029524813, 7.754790468958489, 89.70109119222542, 7.556148365213274, 25.965060177981975, 23.426498016563794, 10.14428547043145, 57.42988902025119, 31.815525977940702, 7.743952994812329, 32.109797045428806, 53.94240851334756, 46.25298270850679, 26.772815264614813, 39.50925386853139, 23.77816653558072, 69.79518841739711, 10.550828769485605, 49.40105418719045, 43.02799056819079, 5.908332187910523, 8.389694091899239, 5.411725324815218, 6.659641162254998, 5.148741916051411, 13.018939389113124, 14.369263879145777, 21.146769197012503, 21.872154900812458, 13.688695176942094, 8.448224543756433, 37.24216527716307, 16.28188339153168, 17.928664213940998, 11.430343282894036, 44.70953331731848, 52.62467879165811, 45.493779866664454, 40.00077846891103, 54.48894618210974, 11.324782041286394, 18.9468971704655, 15.624337007099601, 10.07164556284105, 36.136225125879946, 6.383005547928466, 14.901742383627454, 14.925043733038475, 22.620977183698667, 6.983810820614679, 24.85458182474022, 12.64207354586782, 8.959333766573748, 14.072681049693339, 14.727675607944024, 5.897178780864305, 14.793275357159445, 41.300704270154355, 40.57628120280969, 12.853708620908415, 13.453643611525175, 30.722172602768957, 9.926232490786155, 47.14677218625363, 11.911645024717282, 35.48180106788522, 69.73056913729127, 26.493920191066046, 91.13206307859123, 26.142266212372864, 30.708904276228445, 20.39439571720566, 23.20537404952701, 42.49547458507458, 25.64729266426884, 34.1197785723608, 20.262190240739343, 5.189546737127401, 12.572815226237152, 9.430343393089538, 7.0183345759590745, 48.09402447066543, 6.628072789602321, 9.908964510692961, 13.518875043548238, 37.97285072811504, 39.549782492630676, 9.002334403506342, 22.00836846897292, 49.1403915087156, 6.1593746587870095, 20.334718043361036, 66.37898921317844, 87.91118950629291, 26.38949628689685, 75.24588157355609, 35.69825321773843, 14.673983140228355, 14.82047446937802, 72.71234538104842, 29.86007877256362, 51.52256662333684, 29.625169596203836, 95.2098589382308, 10.266281834717775, 49.506976388462604, 22.206844028616224, 10.483578363836186, 15.082327217535068, 42.24666349694151, 51.87424614935948, 51.48758727009644, 19.646884891683953, 5.293078022010689, 11.093344208484977, 5.8194539965376055, 6.939367469405111, 8.352271377889013, 6.284496889469913, 16.104931226655747, 72.69328131479669, 11.427941787757222, 36.74397531453554, 10.347121196824602, 21.014102576829373, 11.505990441010656, 11.939822363503923, 40.83490349167668, 5.728664673603575, 29.26027164895753, 12.51418926142831, 22.6775898377942, 6.079412548335617, 29.302165490976535, 102.70976793125672, 56.84214704002435, 66.2668457954252, 96.72747174152747, 67.9405082300604, 20.27287876012243, 36.542178314746145, 6.151914691269665, 41.244366789223335, 14.441408443638895, 27.744506314230723, 26.239260501988262, 32.229118398702624, 17.697433462978232, 14.39602176586904, 5.231643332867902, 22.082874885854054, 5.402808866845115, 25.878025141512442, 6.040398486905212, 24.636275886575717, 71.04539707230435, 58.04083447124795, 26.003062357681554, 14.4022288742417, 60.91498228126416, 19.45826261863541, 10.123778425960948, 23.445382235071065, 26.434888263639614, 5.612731546844431, 31.64826045373723, 33.16981878834403, 66.95048665925208, 42.33817539051895, 70.92840005186784, 30.607143968110407, 44.112125804838236, 60.86494062425038, 88.76341085264215, 79.31061842136104, 8.908223246334472, 5.223142510423527, 76.9641490504155, 47.7765082859956, 25.943899801221846, 36.53781248430798, 10.5910286757399, 23.748638609542652, 18.089933854231074, 61.70399034704518, 17.456940543601768, 13.18213828624253, 32.6484229293632, 32.43281571002123, 101.61661638600606, 57.345436713391315, 28.68309170424345, 5.2186775478605, 5.681630259765935, 79.74819338676163, 7.451958603634344, 57.63097491347129, 12.277893560714832, 19.52871576319681, 136.8781468931789, 6.037513420892447, 38.47534471979837, 12.983316575943968, 32.907006897292824, 9.578760087295343, 13.116074587375607, 24.775521522643736, 10.755654488319447, 99.96875360790315, 12.85317541431894, 36.782147607281594, 49.43601960509329, 5.220232653968736, 5.457796766727303, 9.43810473661178, 10.067959735892794, 20.204842789321127, 25.958579603808438, 25.722918264879215, 59.08686006389523, 21.52328120958489, 23.031003790454857, 20.876197297286414, 75.99798096131326, 27.09070877294054, 53.797213417670186, 43.047582219957135, 67.87573448830743, 5.745826318142593, 31.418046575271, 47.14878914859085, 88.18478941213253, 5.280473165712254, 15.48063081445228, 5.046605275281482, 7.080052212698353, 60.49974931357806, 6.233840845351417, 186.37660783624187, 44.38737537705964, 39.542584946948445, 54.17550854876919, 6.1934164678063075, 49.28731723008159, 32.38719736798186, 13.160990412669836, 20.857098742305503, 10.29534003529524, 21.745973892790097, 65.47553683411499, 24.142971358875677, 24.731588505718477, 42.32035782285217, 15.249931475779778, 61.86547253028861, 75.76114057251098, 24.65683997342812, 16.17679519295178, 5.424220788653245, 18.696611551758302, 31.427751664361818, 44.59367498670005, 6.338670706173984, 44.25557778986483, 24.021748303272865, 49.89078313695396, 37.01252132827999, 8.923371271224662, 56.207632262218226, 16.490000663356724, 7.684860379061482, 31.537362101845687, 5.593332306509236, 20.07892381470297, 5.256933266731531, 17.066383251529615, 21.4851241315282, 62.75012960561685, 34.878283989950035, 112.29818462854801, 14.517489635698217, 55.35497554948135, 46.61620988033934, 7.995595640237112, 16.367726420642015, 14.633649664447823, 63.244302879746726, 6.323392369985584, 33.083540173689094, 21.04261395869716, 41.78009336545705, 21.754887736277286, 10.652181316482707, 24.86035827570537, 6.611075669615951, 5.486563150483383, 27.136199183780025, 10.339596730452229, 17.89746605688413, 19.63441586536146, 46.02126251324685, 69.83957315389293, 10.38244692954969, 47.937329206654624, 6.357993353579579, 30.606312484337685, 121.39901888296282, 31.131664872198336, 63.713611657327334, 25.21776310494605, 7.627902233338459, 10.657564454179456, 28.642218443231357, 11.613065391992013, 42.34442597735319, 35.32740052263784, 37.40350165919267, 50.29054156536031, 38.10146008129219, 6.1630078659859695, 47.81007741958481, 14.015993735012309, 13.535613551136759, 12.222022180438449, 21.00914396156214, 5.255335921072931, 23.552704328368158, 61.63082233086503, 34.883074170987534, 111.0614040977351, 54.28632388426466, 21.16683783679268, 10.82600413512975, 13.574385218876655, 36.11431673752004, 34.434675240142354, 21.166802252015376, 42.253559781072866, 48.93302433198453, 49.0087195292572, 5.094514081473814, 53.70087328743548, 14.055154591890712, 9.400016264302172, 85.26016705524454, 26.133648807929344, 27.85875123119011, 18.721963963624365, 27.04448873474811, 18.551434270301595, 88.10343861551631, 5.172812180767204, 7.536165016560863, 84.44553839054439, 94.33528250898142, 14.542033743820722, 41.550086840018096, 31.181586621794544, 6.346373239572135, 14.007272419614747, 22.88840235289654, 84.99552667987346, 8.923836385514063, 6.525574964902848, 57.592678054651316, 27.032001811327987, 16.936495837460107, 21.58473914165049, 49.0053320072829, 6.476124317277085, 7.6381218063281855, 50.22313504984391, 23.68894001012083, 10.058284265745536, 21.349926314564343, 86.82137874304998, 93.30187717605122, 22.457234697770264, 57.364868405822314, 7.2176461450538225, 8.356776482822138, 35.53633266949107, 9.123054006519741, 53.326938448449226, 68.20658738474245, 5.794823176706958, 7.392507385370459, 10.039751759120962, 8.464705776114249, 110.1473650536748, 25.12321833681538, 17.041871140599483, 10.207849835930586, 17.52927864815489, 14.882854672939615, 103.74907638309892, 42.436508502657674, 42.20014002723499, 86.9821261483078, 45.190905517084616, 15.226006590605854, 5.162098224343081, 18.05866892224462, 11.65938425894094, 5.23054428556187, 81.92084227205126, 131.48054999218957, 100.2094161952809, 20.443788902006748, 11.148196921283928, 9.42054530875808, 24.27798270549415, 22.479825784191405, 19.724411713573843, 31.035190944182688, 6.574803553942991, 46.662117481836496, 5.20269142198402, 31.689227212255634, 25.90241841883148, 11.126566503321005, 21.87362598558279, 16.451087858115578, 35.35891668003922, 80.75855591424877, 44.207469370612145, 5.414109078510478, 18.427740172277627, 6.946598965966919, 45.90815609864583, 12.734407656729031, 26.06043557087196, 11.96663407508298, 14.279078972841074, 24.26501670706422, 7.03381356372197, 28.901460383439407, 96.98107554255881, 101.47997996304743, 17.901811991650305, 14.37263516794507, 51.314730012017876, 20.51014678250479, 14.465450029298346, 14.981327034442483, 29.03062224571492, 98.26969569482813, 5.266282271203422, 43.33310546741839, 11.109202027010381, 7.422884720861018, 15.797596666543068, 36.472231883671995, 9.192369551553885, 17.831687671194416, 108.38294338321113, 5.806484111987648, 6.13595417273739, 20.202230120484447, 14.425503779032276])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([5470700.0, 6122045.3125, 6243959.375, 6290701.5625, 6297417.1875, 6312395.3125, 6315159.375, 6317226.5625, 6317284.375, 6318421.875, 6365495.3125, 6365546.875, 6368198.4375, 6372684.375, 6375050.0, 6376064.0625, 6405745.3125, 6413112.5, 6415925.0, 6417312.5, 6434789.0625, 6461846.875, 6544517.1875, 6557629.6875, 6565242.1875, 6565315.625, 6568340.625, 6576395.3125, 6610321.875, 6610728.125, 6611550.0, 6612453.125, 6613289.0625, 6614284.375, 6614582.8125, 6615031.25, 6615117.1875, 6616489.0625, 6638365.625, 6666957.8125, 6671709.375, 6692260.9375, 6727423.4375, 6747229.6875, 6754425.0, 6770390.625, 6782639.0625, 6786725.0, 6790895.3125, 6790925.0, 6791337.5, 6792334.375, 6792496.875, 6792654.6875, 6793148.4375, 6793245.3125, 6793326.5625, 6794039.0625, 6794529.6875, 6795168.75, 6797026.5625, 6809607.8125, 6841125.0, 6841576.5625, 6842840.625, 6843528.125, 6844050.0, 6844876.5625, 6845046.875, 6845260.9375, 6845435.9375, 6931915.625, 6932621.875, 6939903.125, 6948675.0, 6948954.6875, 6976082.8125, 6977854.6875, 6995010.9375, 6995221.875, 6999707.8125, 7010454.6875, 7024406.25, 7051040.625, 7056859.375, 7087010.9375, 7152462.5, 7191401.5625, 7203760.9375, 7206790.625, 7227307.8125, 7268150.0, 7274557.8125, 7303523.4375, 7483157.8125, 7515953.125, 7532070.3125, 7585101.5625, 7592435.9375, 7617507.8125, 7640865.625, 8158739.0625, 8341992.1875, 8742915.625, 8745429.6875, 8896539.0625, 8897218.75, 8900353.125, 8922768.75, 8952323.4375, 9016885.9375, 9053489.0625, 9056553.125, 9071548.4375, 9081501.5625, 9083954.6875, 9084728.125, 9093759.375, 9184885.9375, 9185790.625, 9186853.125, 9187056.25, 9187242.1875, 9187826.5625, 9187940.625, 9188356.25, 9188840.625, 9189709.375, 9190681.25, 9191962.5, 9211017.1875, 9211059.375, 9212698.4375, 9213515.625, 9214260.9375, 9215653.125, 9216051.5625, 9216325.0, 9216329.6875, 9218117.1875, 9219581.25, 9224764.0625, 9225714.0625, 9231679.6875, 9239634.375, 9241787.5, 9249125.0, 9250290.625, 9251634.375, 9252456.25, 9255898.4375, 9283743.75, 9344678.125, 9345293.75, 9345715.625, 9351218.75, 9358037.5, 9366471.875, 9368329.6875, 9368656.25, 9368768.75, 9368971.875, 9370218.75, 9371384.375, 9372092.1875, 9372323.4375, 9378300.0, 9378612.5, 9380526.5625, 9387962.5, 9388795.3125, 9391978.125, 9392406.25, 9393000.0, 9393389.0625, 9401007.8125, 9402292.1875, 9404040.625, 9411431.25, 9411435.9375, 9423295.3125, 9424525.0, 9447748.4375, 9448029.6875, 9449560.9375, 9449623.4375, 9449703.125, 9450314.0625, 9452034.375, 9452543.75, 9452779.6875, 9453740.625, 9455087.5, 9456306.25, 9456981.25, 9456992.1875, 9495570.3125, 9496193.75, 9497059.375, 9497332.8125, 9497600.0, 9497871.875, 9499110.9375, 9499148.4375, 9500146.875, 9507839.0625, 9508890.625, 9510567.1875, 9511082.8125, 9511800.0, 9522014.0625, 9525356.25, 9526981.25, 9529025.0, 9529437.5, 9538645.3125, 9541935.9375, 9542082.8125, 9543035.9375, 9543720.3125, 9610778.125, 9619928.125, 9638068.75, 9715906.25, 9717585.9375, 9719315.625, 9741928.125, 9742382.8125, 9743175.0, 9743557.8125, 9744164.0625, 9744270.3125, 9747600.0, 9747651.5625, 9749429.6875, 9751496.875, 9767779.6875, 9773337.5, 9776160.9375, 9786034.375, 9786995.3125, 9791762.5, 9796870.3125, 9796893.75, 9806301.5625, 9811043.75, 9890292.1875, 9906764.0625, 9925570.3125, 9931945.3125, 9965332.8125, 9969200.0, 9969675.0, 9969748.4375, 9973417.1875, 9976837.5, 9988881.25, 10009853.125, 10011873.4375, 10013937.5, 10023450.0, 10023475.0, 10023487.5, 10038245.3125, 10100821.875, 10103193.75, 10104520.3125, 10105098.4375, 10115637.5, 10132003.125, 10133939.0625, 10141028.125, 10171753.125, 10204665.625, 10235901.5625, 10257421.875, 10263312.5, 10297973.4375, 10316348.4375, 10421642.1875, 10515720.3125, 10520773.4375, 10521157.8125, 10557134.375, 10616289.0625, 10664932.8125, 11031510.9375, 11081223.4375, 11087146.875, 11148254.6875, 11150075.0, 11154546.875, 11157140.625, 11157390.625, 11158046.875, 11159251.5625, 11163929.6875, 11172446.875, 11178334.375, 11181218.75, 11181789.0625, 11181923.4375, 11182198.4375, 11185790.625, 11186618.75, 11186893.75, 11187134.375, 11187229.6875, 11187509.375, 11187545.3125, 11187996.875, 11188053.125, 11188093.75, 11188289.0625, 11188298.4375, 11188331.25, 11188475.0, 11188684.375, 11189248.4375, 11189285.9375, 11189503.125, 11189548.4375, 11189592.1875, 11189967.1875, 11189973.4375, 11190028.125, 11190104.6875, 11190146.875, 11190362.5, 11190379.6875, 11190412.5, 11190440.625, 11190889.0625, 11190901.5625, 11190912.5, 11190926.5625, 11190970.3125, 11190984.375, 11191029.6875, 11191092.1875, 11191110.9375, 11191132.8125, 11191142.1875, 11191226.5625, 11191315.625, 11191331.25, 11191434.375, 11191521.875, 11191575.0, 11191585.9375, 11191670.3125, 11191746.875, 11191912.5, 11192279.6875, 11192371.875, 11192506.25, 11192553.125, 11192556.25, 11192593.75, 11192596.875, 11192756.25, 11192851.5625, 11192982.8125, 11193021.875, 11193568.75, 11193670.3125, 11193929.6875, 11195054.6875, 11195354.6875, 11195404.6875, 11195704.6875, 11195765.625, 11195878.125, 11196057.8125, 11196545.3125, 11197056.25, 11197087.5, 11197525.0, 11197757.8125, 11197810.9375, 11198140.625, 11198479.6875, 11198828.125, 11198859.375, 11199081.25, 11199135.9375, 11199156.25, 11199279.6875, 11199400.0, 11199479.6875, 11199698.4375, 11199834.375, 11199957.8125, 11200068.75, 11200262.5, 11200337.5, 11200340.625, 11200615.625, 11200709.375, 11200728.125, 11200754.6875, 11200759.375, 11200873.4375, 11201026.5625, 11201135.9375, 11201300.0, 11201323.4375, 11201425.0, 11201487.5, 11201720.3125, 11201721.875, 11201723.4375, 11201868.75, 11202296.875, 11202446.875, 11202575.0, 11202614.0625, 11202729.6875, 11203282.8125, 11203473.4375, 11203918.75, 11203931.25, 11203984.375, 11204048.4375, 11204375.0, 11204503.125, 11204721.875, 11204970.3125, 11204973.4375, 11205012.5, 11205246.875, 11205748.4375, 11205829.6875, 11205895.3125, 11205915.625, 11206426.5625, 11206657.8125, 11207381.25, 11207428.125, 11207442.1875, 11207800.0, 11208667.1875, 11208676.5625, 11209234.375, 11209435.9375, 11209575.0, 11209787.5, 11210067.1875, 11210945.3125, 11211142.1875, 11211165.625, 11211206.25, 11211392.1875, 11211953.125, 11212381.25, 11212504.6875, 11212795.3125, 11213148.4375, 11213168.75, 11213350.0, 11213417.1875, 11213596.875, 11213603.125, 11214114.0625, 11214151.5625, 11214151.5625, 11214232.8125, 11214279.6875, 11214700.0, 11214840.625, 11214898.4375, 11214903.125, 11214906.25, 11215139.0625, 11215178.125, 11215231.25, 11215243.75, 11215276.5625, 11215678.125, 11215817.1875, 11215976.5625, 11215985.9375, 11216157.8125, 11216309.375, 11216551.5625, 11217090.625, 11217496.875, 11217503.125, 11217675.0, 11217793.75, 11217915.625, 11217970.3125, 11218021.875, 11218479.6875, 11218504.6875, 11218801.5625, 11218946.875, 11219442.1875, 11219942.1875, 11219964.0625, 11220021.875, 11220178.125, 11220557.8125, 11220560.9375, 11220590.625, 11220715.625, 11220875.0, 11220925.0, 11220970.3125, 11221028.125, 11221162.5, 11221514.0625, 11221576.5625, 11221579.6875, 11221593.75, 11221625.0, 11221628.125, 11221654.6875, 11221712.5, 11221957.8125, 11222310.9375, 11222510.9375, 11222565.625, 11222570.3125, 11222685.9375, 11222725.0, 11223121.875, 11223135.9375, 11223337.5, 11223454.6875, 11223609.375, 11223725.0, 11223840.625, 11224075.0, 11224110.9375, 11224146.875, 11224329.6875, 11224376.5625, 11224506.25, 11224570.3125, 11224675.0, 11224732.8125, 11224953.125, 11224957.8125, 11225084.375, 11225321.875, 11225578.125, 11225617.1875, 11225681.25, 11225715.625, 11225753.125, 11225798.4375, 11226046.875, 11226140.625, 11226384.375, 11226423.4375, 11226467.1875, 11226467.1875, 11226482.8125, 11226514.0625, 11226540.625, 11226557.8125, 11226567.1875, 11226576.5625, 11226651.5625, 11226818.75, 11226828.125, 11226860.9375, 11226870.3125, 11226992.1875, 11227156.25, 11227181.25, 11227203.125, 11227292.1875, 11227292.1875, 11227478.125, 11227657.8125, 11227671.875, 11227723.4375, 11227740.625, 11227809.375, 11227812.5, 11227842.1875, 11227895.3125, 11227901.5625, 11227907.8125, 11228212.5, 11228301.5625, 11228368.75, 11228375.0, 11228400.0, 11228468.75, 11228610.9375, 11228725.0, 11228742.1875, 11228815.625, 11228875.0, 11228882.8125, 11229112.5, 11229334.375, 11229376.5625, 11229526.5625, 11229543.75, 11229554.6875, 11229739.0625, 11229798.4375, 11229803.125, 11229870.3125, 11229929.6875, 11229985.9375, 11230042.1875, 11230084.375, 11230106.25, 11230178.125, 11230264.0625, 11230359.375, 11230414.0625, 11230500.0, 11230635.9375, 11230651.5625, 11230670.3125, 11230912.5, 11231665.625, 11231750.0, 11231854.6875, 11231867.1875, 11232121.875, 11232148.4375, 11232245.3125, 11232439.0625, 11232779.6875, 11232881.25, 11233065.625, 11233134.375, 11233137.5, 11233325.0, 11233395.3125, 11233398.4375, 11233868.75, 11234004.6875, 11234018.75, 11234459.375, 11234584.375, 11234676.5625, 11235034.375, 11235103.125, 11235625.0, 11235676.5625, 11235748.4375, 11235876.5625, 11235884.375, 11236228.125, 11236296.875, 11236773.4375, 11237193.75, 11237264.0625, 11237368.75, 11237385.9375, 11237579.6875, 11237826.5625, 11238039.0625, 11238207.8125, 11238731.25, 11238868.75, 11239037.5, 11239109.375, 11239160.9375, 11239621.875, 11240298.4375, 11240473.4375, 11240667.1875, 11241101.5625, 11241221.875, 11241307.8125, 11241346.875, 11241737.5, 11241751.5625, 11241756.25, 11241821.875, 11241904.6875, 11241917.1875, 11241943.75, 11242012.5, 11242014.0625, 11242056.25, 11242626.5625, 11242646.875, 11242710.9375, 11243243.75, 11243496.875, 11243656.25, 11243726.5625, 11244084.375, 11244337.5, 11245167.1875, 11245357.8125, 11245587.5, 11245631.25, 11245703.125, 11245731.25, 11245851.5625, 11245975.0, 11246042.1875, 11246170.3125, 11246182.8125, 11246904.6875, 11247129.6875, 11247309.375, 11247362.5, 11247471.875, 11247734.375, 11247900.0, 11248231.25, 11248239.0625, 11248315.625, 11248387.5, 11248528.125, 11248682.8125, 11248753.125, 11248785.9375, 11248892.1875, 11249142.1875, 11249175.0, 11249695.3125, 11249734.375, 11250970.3125, 11251001.5625, 11251384.375, 11251653.125, 11251717.1875, 11251846.875, 11252343.75, 11252778.125, 11253340.625, 11254192.1875, 11254642.1875, 11254710.9375, 11254740.625, 11254776.5625, 11254784.375, 11255078.125, 11255620.3125, 11255815.625, 11255928.125, 11256446.875, 11256446.875, 11256614.0625, 11257279.6875, 11257450.0, 11258204.6875, 11258609.375, 11258634.375, 11259451.5625, 11259792.1875, 11260662.5, 11260715.625, 11260859.375, 11261665.625, 11262460.9375, 11262564.0625, 11262565.625, 11262625.0, 11264050.0, 11264137.5, 11265189.0625, 11265687.5, 11265751.5625, 11265840.625, 11265943.75, 11266532.8125, 11267440.625, 11267560.9375, 11267832.8125, 11268154.6875, 11269114.0625, 11269442.1875, 11269500.0, 11269903.125, 11269984.375, 11270670.3125, 11271157.8125, 11271178.125, 11271323.4375, 11271457.8125, 11271507.8125, 11271621.875, 11271768.75, 11271820.3125, 11271985.9375, 11272401.5625, 11272409.375, 11272700.0, 11272768.75, 11272804.6875, 11272954.6875, 11273239.0625, 11273246.875, 11273371.875, 11273478.125, 11273490.625, 11273871.875, 11274337.5, 11274381.25, 11274671.875, 11274860.9375, 11274910.9375, 11275226.5625, 11275287.5, 11275298.4375, 11275357.8125, 11275526.5625, 11275825.0, 11276135.9375, 11276182.8125, 11276209.375, 11276214.0625, 11276262.5, 11276278.125, 11276431.25, 11276592.1875, 11276671.875, 11276798.4375, 11276882.8125, 11277109.375, 11277114.0625, 11277500.0, 11277615.625, 11277729.6875, 11277864.0625, 11278112.5, 11278159.375, 11278257.8125, 11278439.0625, 11278887.5, 11278912.5, 11278935.9375, 11278951.5625, 11279217.1875, 11279303.125, 11279500.0, 11279790.625, 11280242.1875, 11280390.625, 11280829.6875, 11280932.8125, 11280975.0, 11281148.4375, 11282146.875, 11282423.4375, 11282645.3125, 11283068.75, 11283117.1875, 11283446.875, 11283675.0, 11285146.875, 11285829.6875, 11286868.75, 11286884.375, 11287160.9375, 11287489.0625, 11288742.1875, 11289150.0, 11289871.875, 11290551.5625, 11290565.625, 11290871.875, 11292012.5, 11293326.5625, 11293557.8125, 11294004.6875, 11294284.375, 11297067.1875, 11298284.375, 11301239.0625, 11307773.4375], [19.2154295523539, 89.55258854213393, 14.370355298661053, 57.20399228797166, 44.811780504814074, 26.75388049298675, 6.530170367706963, 46.51528292778501, 5.291261754564691, 11.241879054981549, 24.51955003314039, 5.4735318475161385, 16.98705148080214, 6.854651163154591, 134.54444194779637, 108.56779278433693, 29.617489241934027, 51.53468854959295, 45.763343383599434, 29.438834060151855, 8.984525156225294, 11.413458231150608, 79.11277376253217, 29.332268207312964, 7.754292493196194, 5.1596478858584565, 65.75759256528346, 29.886755584893912, 26.5584010269535, 23.863277590325925, 14.547172927251028, 60.30310941298213, 19.09707671046739, 156.3456596865374, 17.16082850848504, 46.84886998523748, 22.061472761700315, 21.5868778827754, 9.24097056919769, 23.262784910507882, 118.25907569736367, 61.11357374615261, 38.57253171293951, 47.46288157118514, 34.99858159090979, 45.693231652771615, 9.251113344025454, 8.059526900380703, 90.49611075964678, 20.520631561015886, 11.720049670994765, 29.740928249506485, 107.58167149662388, 6.40105327770997, 28.50712426756071, 13.07560101312344, 68.26889573018359, 17.74258765455825, 56.945125949754264, 72.29505848984574, 6.020204843413952, 44.54701547359016, 40.36894174918001, 17.067673224056758, 9.605313045803141, 18.59978060841521, 57.41665724064427, 217.70902154594663, 60.42647272904755, 64.01513099818922, 5.657883631660489, 69.18948929366888, 12.09263699146219, 13.107291736562065, 39.879337497829894, 22.24404132702483, 23.154982776707463, 123.23317597273478, 53.05712160521556, 25.354345409682807, 6.980673241152645, 30.94223972389825, 10.214020951023041, 142.33953166532407, 297.29521644971095, 10.760884240449247, 23.924188692267542, 145.2495245410926, 7.913224832692245, 70.96656852421646, 30.52159985969437, 99.8772153032514, 49.93174090042638, 30.353901470667097, 9.26361765387983, 21.240483298665268, 9.544389337823151, 33.41682334052246, 124.01046154983763, 31.81760139138788, 97.40791875541008, 31.90252501810467, 26.154246510660762, 5.140815620989394, 17.177225782037713, 64.50741966314932, 35.957545738789264, 104.5144784399327, 26.64889928125563, 59.48767625255684, 16.942852235320416, 134.40837950624547, 22.053034040793786, 6.23291908819242, 5.911287509029278, 110.56842971665007, 7.677522869714147, 6.284719157619654, 42.683435128956695, 8.32042025593095, 9.928282408649084, 19.722045606635636, 28.987548108221837, 22.541174969365414, 100.96387202082396, 5.803734727703886, 19.65109117609528, 24.842907083337348, 30.680043183427458, 8.658986238020574, 26.621158222116563, 35.042213976536594, 56.984986104709336, 22.89199804998787, 12.055277195586118, 67.62044480693226, 16.47506936148762, 31.985816272631364, 20.16372603464425, 25.89806132220823, 28.880419122242554, 15.46104680811437, 21.045317789123715, 45.491566144160146, 85.2147434472077, 17.309245418570566, 8.078382959682134, 8.195185905875588, 20.622677176985167, 6.1076938727813745, 137.97297094019163, 40.72741328583843, 105.8852169166256, 59.00638138940108, 37.84499218295711, 18.052755338778198, 12.777847454904615, 18.981400770137018, 166.2523480868819, 119.7609513969797, 114.13651494861878, 16.652970221240064, 46.62280759256032, 195.60053890605388, 46.59376171522019, 84.21447092167637, 35.67032040601055, 54.621348313428555, 22.53155975251148, 43.436918197163806, 15.338022649426478, 38.304689734371735, 15.100770425629005, 44.81912929448974, 6.320434263713759, 28.76550382602135, 107.88435606427258, 8.450202504917593, 84.97320765353425, 37.56270282959764, 58.030205410898795, 5.7709920107508434, 32.06595252128765, 17.993235980718183, 17.6253432354124, 41.00399468709798, 59.93992228451643, 84.16686319144326, 10.351134411316657, 13.5293570573456, 49.50857954099376, 34.1751830579345, 63.41813705918761, 10.368743511179968, 18.738693240509853, 6.657544770175715, 52.444248361771116, 49.51492145254019, 24.429689187005337, 12.706705068052425, 53.697597530578605, 8.236528576630326, 157.20851863249132, 34.22233017568783, 24.64701976296206, 47.981146511860544, 56.64239394103669, 5.7529193433668295, 21.83882126418456, 19.57399568647594, 32.83299838602224, 298.4792889812215, 98.62888764356468, 26.47639052807266, 55.785696132754126, 19.240462877606962, 18.30042911850452, 55.666689248636914, 6.032081601307642, 103.43258385604304, 7.426222274868422, 175.87626436116688, 5.208550669803976, 87.66195913679998, 6.970929251076096, 101.57004136135618, 51.365073943841615, 60.9338213153302, 52.62441405820311, 12.093841119121558, 79.93321541879693, 131.19417342765024, 8.139494005406712, 5.821175516598972, 8.455144136445806, 70.4630275673035, 34.323304711144594, 84.58534687972843, 108.64944103920774, 13.681570178809102, 62.22614937523675, 12.365094712858127, 28.492406438982332, 11.695273511376318, 52.46457819806943, 80.60044291018835, 66.2486926663718, 74.93155719876815, 5.887282941632003, 18.54482833822992, 10.659805712491513, 106.2018419397337, 78.7222965759835, 5.176864716447867, 47.39263194185077, 58.657984790067154, 12.137771994305721, 5.6537688156527866, 57.06383264487278, 15.620701400616198, 48.72667263087499, 6.868133893677933, 7.624943522082688, 6.511652217974223, 28.12847202639719, 6.3308568389160085, 85.81788067101157, 39.5355884853977, 25.799927371278617, 5.635185781668891, 19.413694709855502, 103.66360819193886, 39.27383368194472, 13.696639081188916, 44.29238541410025, 110.26364335749889, 20.176164818304052, 11.964093289182197, 48.53086118336387, 49.204634448661025, 6.970078167026922, 8.935250347429678, 47.064436161810555, 98.01106467781845, 13.380820703223215, 35.649780551146065, 8.416738377800952, 6.425517229849981, 85.7151444578088, 12.068191739864806, 10.415593125635109, 8.972398263117928, 42.37344961171816, 58.2661475358052, 23.319547183968105, 73.68292085277879, 80.81713775262017, 38.13017052214444, 7.693987769674597, 16.377388457263514, 86.45912336994074, 93.21326374404448, 52.01944740977451, 18.291413416405508, 7.514705503132166, 40.84015074842837, 17.22129104662929, 6.398440824001257, 20.405598255080214, 25.33041735696217, 13.445138437343513, 92.95164637472942, 41.36511519869156, 65.44956925553036, 55.155084595589955, 30.00283269197329, 19.612254510784922, 32.83023133034581, 44.758132906000625, 15.308796375793193, 12.413608308115478, 10.68845631053543, 9.245130215340437, 16.671255064750632, 37.374326254738456, 5.777629518547999, 35.04937588090091, 17.980773221899423, 33.13911708660834, 15.058198397742165, 76.7816122129089, 6.377286346592362, 44.489582365049124, 45.18491444488918, 10.724624285581632, 14.070434476963813, 45.13512991270552, 12.609869601272743, 10.595083751588458, 20.89065575264003, 48.26806553866681, 36.10995743514956, 7.7460302051703, 7.505432818216993, 25.27096680329298, 5.128056347923415, 8.081620682470584, 16.787434254518708, 12.548672309875796, 7.309578165500295, 29.020397725107284, 15.819777974452645, 21.842558275844425, 42.749059836789236, 56.04955688135264, 18.431463067382307, 41.86821160239721, 24.012749600134843, 81.63171093430608, 27.620072308716992, 9.997263199271416, 7.754638245339949, 60.24650286862469, 23.835065228706924, 7.401011322532126, 15.703710682886056, 7.13657384749423, 13.657442972863167, 13.416787128507428, 31.59162517362728, 39.78198164977961, 15.291569046090522, 21.48855460426641, 146.21183152002635, 59.78011245708154, 70.81183821481665, 141.7277492793319, 26.403224294035706, 39.51417239946688, 30.897482814382986, 17.552609536530802, 9.58766721154342, 77.56996682126315, 17.147744434704386, 31.058917482281707, 11.290727540660633, 20.224815672530422, 65.82104401212247, 72.92231591405691, 6.585060122686211, 11.125336338174773, 5.866850039796829, 16.118845809924984, 63.95025140890468, 67.25577173679876, 13.384798335304454, 82.94665171031781, 35.619306463418056, 23.205619874504407, 74.20522964493901, 7.340666706822463, 14.749773076930264, 9.271464809005398, 46.48912410617532, 7.948426370465435, 23.756361712600018, 6.088395514391004, 6.302727603136325, 71.53022670204054, 35.021021869095875, 17.623555629125867, 14.411791763744088, 13.628827912187504, 5.404040694168145, 17.53156105658011, 19.99823029524813, 7.754790468958489, 89.70109119222542, 7.556148365213274, 25.965060177981975, 23.426498016563794, 10.14428547043145, 57.42988902025119, 31.815525977940702, 7.743952994812329, 32.109797045428806, 53.94240851334756, 46.25298270850679, 26.772815264614813, 39.50925386853139, 23.77816653558072, 69.79518841739711, 10.550828769485605, 49.40105418719045, 43.02799056819079, 5.908332187910523, 8.389694091899239, 5.411725324815218, 6.659641162254998, 5.148741916051411, 13.018939389113124, 14.369263879145777, 21.146769197012503, 21.872154900812458, 13.688695176942094, 8.448224543756433, 37.24216527716307, 16.28188339153168, 17.928664213940998, 11.430343282894036, 44.70953331731848, 52.62467879165811, 45.493779866664454, 40.00077846891103, 54.48894618210974, 11.324782041286394, 18.9468971704655, 15.624337007099601, 10.07164556284105, 36.136225125879946, 6.383005547928466, 14.901742383627454, 14.925043733038475, 22.620977183698667, 6.983810820614679, 24.85458182474022, 12.64207354586782, 8.959333766573748, 14.072681049693339, 14.727675607944024, 5.897178780864305, 14.793275357159445, 41.300704270154355, 40.57628120280969, 12.853708620908415, 13.453643611525175, 30.722172602768957, 9.926232490786155, 47.14677218625363, 11.911645024717282, 35.48180106788522, 69.73056913729127, 26.493920191066046, 91.13206307859123, 26.142266212372864, 30.708904276228445, 20.39439571720566, 23.20537404952701, 42.49547458507458, 25.64729266426884, 34.1197785723608, 20.262190240739343, 5.189546737127401, 12.572815226237152, 9.430343393089538, 7.0183345759590745, 48.09402447066543, 6.628072789602321, 9.908964510692961, 13.518875043548238, 37.97285072811504, 39.549782492630676, 9.002334403506342, 22.00836846897292, 49.1403915087156, 6.1593746587870095, 20.334718043361036, 66.37898921317844, 87.91118950629291, 26.38949628689685, 75.24588157355609, 35.69825321773843, 14.673983140228355, 14.82047446937802, 72.71234538104842, 29.86007877256362, 51.52256662333684, 29.625169596203836, 95.2098589382308, 10.266281834717775, 49.506976388462604, 22.206844028616224, 10.483578363836186, 15.082327217535068, 42.24666349694151, 51.87424614935948, 51.48758727009644, 19.646884891683953, 5.293078022010689, 11.093344208484977, 5.8194539965376055, 6.939367469405111, 8.352271377889013, 6.284496889469913, 16.104931226655747, 72.69328131479669, 11.427941787757222, 36.74397531453554, 10.347121196824602, 21.014102576829373, 11.505990441010656, 11.939822363503923, 40.83490349167668, 5.728664673603575, 29.26027164895753, 12.51418926142831, 22.6775898377942, 6.079412548335617, 29.302165490976535, 102.70976793125672, 56.84214704002435, 66.2668457954252, 96.72747174152747, 67.9405082300604, 20.27287876012243, 36.542178314746145, 6.151914691269665, 41.244366789223335, 14.441408443638895, 27.744506314230723, 26.239260501988262, 32.229118398702624, 17.697433462978232, 14.39602176586904, 5.231643332867902, 22.082874885854054, 5.402808866845115, 25.878025141512442, 6.040398486905212, 24.636275886575717, 71.04539707230435, 58.04083447124795, 26.003062357681554, 14.4022288742417, 60.91498228126416, 19.45826261863541, 10.123778425960948, 23.445382235071065, 26.434888263639614, 5.612731546844431, 31.64826045373723, 33.16981878834403, 66.95048665925208, 42.33817539051895, 70.92840005186784, 30.607143968110407, 44.112125804838236, 60.86494062425038, 88.76341085264215, 79.31061842136104, 8.908223246334472, 5.223142510423527, 76.9641490504155, 47.7765082859956, 25.943899801221846, 36.53781248430798, 10.5910286757399, 23.748638609542652, 18.089933854231074, 61.70399034704518, 17.456940543601768, 13.18213828624253, 32.6484229293632, 32.43281571002123, 101.61661638600606, 57.345436713391315, 28.68309170424345, 5.2186775478605, 5.681630259765935, 79.74819338676163, 7.451958603634344, 57.63097491347129, 12.277893560714832, 19.52871576319681, 136.8781468931789, 6.037513420892447, 38.47534471979837, 12.983316575943968, 32.907006897292824, 9.578760087295343, 13.116074587375607, 24.775521522643736, 10.755654488319447, 99.96875360790315, 12.85317541431894, 36.782147607281594, 49.43601960509329, 5.220232653968736, 5.457796766727303, 9.43810473661178, 10.067959735892794, 20.204842789321127, 25.958579603808438, 25.722918264879215, 59.08686006389523, 21.52328120958489, 23.031003790454857, 20.876197297286414, 75.99798096131326, 27.09070877294054, 53.797213417670186, 43.047582219957135, 67.87573448830743, 5.745826318142593, 31.418046575271, 47.14878914859085, 88.18478941213253, 5.280473165712254, 15.48063081445228, 5.046605275281482, 7.080052212698353, 60.49974931357806, 6.233840845351417, 186.37660783624187, 44.38737537705964, 39.542584946948445, 54.17550854876919, 6.1934164678063075, 49.28731723008159, 32.38719736798186, 13.160990412669836, 20.857098742305503, 10.29534003529524, 21.745973892790097, 65.47553683411499, 24.142971358875677, 24.731588505718477, 42.32035782285217, 15.249931475779778, 61.86547253028861, 75.76114057251098, 24.65683997342812, 16.17679519295178, 5.424220788653245, 18.696611551758302, 31.427751664361818, 44.59367498670005, 6.338670706173984, 44.25557778986483, 24.021748303272865, 49.89078313695396, 37.01252132827999, 8.923371271224662, 56.207632262218226, 16.490000663356724, 7.684860379061482, 31.537362101845687, 5.593332306509236, 20.07892381470297, 5.256933266731531, 17.066383251529615, 21.4851241315282, 62.75012960561685, 34.878283989950035, 112.29818462854801, 14.517489635698217, 55.35497554948135, 46.61620988033934, 7.995595640237112, 16.367726420642015, 14.633649664447823, 63.244302879746726, 6.323392369985584, 33.083540173689094, 21.04261395869716, 41.78009336545705, 21.754887736277286, 10.652181316482707, 24.86035827570537, 6.611075669615951, 5.486563150483383, 27.136199183780025, 10.339596730452229, 17.89746605688413, 19.63441586536146, 46.02126251324685, 69.83957315389293, 10.38244692954969, 47.937329206654624, 6.357993353579579, 30.606312484337685, 121.39901888296282, 31.131664872198336, 63.713611657327334, 25.21776310494605, 7.627902233338459, 10.657564454179456, 28.642218443231357, 11.613065391992013, 42.34442597735319, 35.32740052263784, 37.40350165919267, 50.29054156536031, 38.10146008129219, 6.1630078659859695, 47.81007741958481, 14.015993735012309, 13.535613551136759, 12.222022180438449, 21.00914396156214, 5.255335921072931, 23.552704328368158, 61.63082233086503, 34.883074170987534, 111.0614040977351, 54.28632388426466, 21.16683783679268, 10.82600413512975, 13.574385218876655, 36.11431673752004, 34.434675240142354, 21.166802252015376, 42.253559781072866, 48.93302433198453, 49.0087195292572, 5.094514081473814, 53.70087328743548, 14.055154591890712, 9.400016264302172, 85.26016705524454, 26.133648807929344, 27.85875123119011, 18.721963963624365, 27.04448873474811, 18.551434270301595, 88.10343861551631, 5.172812180767204, 7.536165016560863, 84.44553839054439, 94.33528250898142, 14.542033743820722, 41.550086840018096, 31.181586621794544, 6.346373239572135, 14.007272419614747, 22.88840235289654, 84.99552667987346, 8.923836385514063, 6.525574964902848, 57.592678054651316, 27.032001811327987, 16.936495837460107, 21.58473914165049, 49.0053320072829, 6.476124317277085, 7.6381218063281855, 50.22313504984391, 23.68894001012083, 10.058284265745536, 21.349926314564343, 86.82137874304998, 93.30187717605122, 22.457234697770264, 57.364868405822314, 7.2176461450538225, 8.356776482822138, 35.53633266949107, 9.123054006519741, 53.326938448449226, 68.20658738474245, 5.794823176706958, 7.392507385370459, 10.039751759120962, 8.464705776114249, 110.1473650536748, 25.12321833681538, 17.041871140599483, 10.207849835930586, 17.52927864815489, 14.882854672939615, 103.74907638309892, 42.436508502657674, 42.20014002723499, 86.9821261483078, 45.190905517084616, 15.226006590605854, 5.162098224343081, 18.05866892224462, 11.65938425894094, 5.23054428556187, 81.92084227205126, 131.48054999218957, 100.2094161952809, 20.443788902006748, 11.148196921283928, 9.42054530875808, 24.27798270549415, 22.479825784191405, 19.724411713573843, 31.035190944182688, 6.574803553942991, 46.662117481836496, 5.20269142198402, 31.689227212255634, 25.90241841883148, 11.126566503321005, 21.87362598558279, 16.451087858115578, 35.35891668003922, 80.75855591424877, 44.207469370612145, 5.414109078510478, 18.427740172277627, 6.946598965966919, 45.90815609864583, 12.734407656729031, 26.06043557087196, 11.96663407508298, 14.279078972841074, 24.26501670706422, 7.03381356372197, 28.901460383439407, 96.98107554255881, 101.47997996304743, 17.901811991650305, 14.37263516794507, 51.314730012017876, 20.51014678250479, 14.465450029298346, 14.981327034442483, 29.03062224571492, 98.26969569482813, 5.266282271203422, 43.33310546741839, 11.109202027010381, 7.422884720861018, 15.797596666543068, 36.472231883671995, 9.192369551553885, 17.831687671194416, 108.38294338321113, 5.806484111987648, 6.13595417273739, 20.202230120484447, 14.425503779032276])
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);
([5470700.0, 6122045.3125, 6243959.375, 6290701.5625, 6297417.1875, 6312395.3125, 6315159.375, 6317226.5625, 6317284.375, 6318421.875, 6365495.3125, 6365546.875, 6368198.4375, 6372684.375, 6375050.0, 6376064.0625, 6405745.3125, 6413112.5, 6415925.0, 6417312.5, 6434789.0625, 6461846.875, 6544517.1875, 6557629.6875, 6565242.1875, 6565315.625, 6568340.625, 6576395.3125, 6610321.875, 6610728.125, 6611550.0, 6612453.125, 6613289.0625, 6614284.375, 6614582.8125, 6615031.25, 6615117.1875, 6616489.0625, 6638365.625, 6666957.8125, 6671709.375, 6692260.9375, 6727423.4375, 6747229.6875, 6754425.0, 6770390.625, 6782639.0625, 6786725.0, 6790895.3125, 6790925.0, 6791337.5, 6792334.375, 6792496.875, 6792654.6875, 6793148.4375, 6793245.3125, 6793326.5625, 6794039.0625, 6794529.6875, 6795168.75, 6797026.5625, 6809607.8125, 6841125.0, 6841576.5625, 6842840.625, 6843528.125, 6844050.0, 6844876.5625, 6845046.875, 6845260.9375, 6845435.9375, 6931915.625, 6932621.875, 6939903.125, 6948675.0, 6948954.6875, 6976082.8125, 6977854.6875, 6995010.9375, 6995221.875, 6999707.8125, 7010454.6875, 7024406.25, 7051040.625, 7056859.375, 7087010.9375, 7152462.5, 7191401.5625, 7203760.9375, 7206790.625, 7227307.8125, 7268150.0, 7274557.8125, 7303523.4375, 7483157.8125, 7515953.125, 7532070.3125, 7585101.5625, 7592435.9375, 7617507.8125, 7640865.625, 8158739.0625, 8341992.1875, 8742915.625, 8745429.6875, 8896539.0625, 8897218.75, 8900353.125, 8922768.75, 8952323.4375, 9016885.9375, 9053489.0625, 9056553.125, 9071548.4375, 9081501.5625, 9083954.6875, 9084728.125, 9093759.375, 9184885.9375, 9185790.625, 9186853.125, 9187056.25, 9187242.1875, 9187826.5625, 9187940.625, 9188356.25, 9188840.625, 9189709.375, 9190681.25, 9191962.5, 9211017.1875, 9211059.375, 9212698.4375, 9213515.625, 9214260.9375, 9215653.125, 9216051.5625, 9216325.0, 9216329.6875, 9218117.1875, 9219581.25, 9224764.0625, 9225714.0625, 9231679.6875, 9239634.375, 9241787.5, 9249125.0, 9250290.625, 9251634.375, 9252456.25, 9255898.4375, 9283743.75, 9344678.125, 9345293.75, 9345715.625, 9351218.75, 9358037.5, 9366471.875, 9368329.6875, 9368656.25, 9368768.75, 9368971.875, 9370218.75, 9371384.375, 9372092.1875, 9372323.4375, 9378300.0, 9378612.5, 9380526.5625, 9387962.5, 9388795.3125, 9391978.125, 9392406.25, 9393000.0, 9393389.0625, 9401007.8125, 9402292.1875, 9404040.625, 9411431.25, 9411435.9375, 9423295.3125, 9424525.0, 9447748.4375, 9448029.6875, 9449560.9375, 9449623.4375, 9449703.125, 9450314.0625, 9452034.375, 9452543.75, 9452779.6875, 9453740.625, 9455087.5, 9456306.25, 9456981.25, 9456992.1875, 9495570.3125, 9496193.75, 9497059.375, 9497332.8125, 9497600.0, 9497871.875, 9499110.9375, 9499148.4375, 9500146.875, 9507839.0625, 9508890.625, 9510567.1875, 9511082.8125, 9511800.0, 9522014.0625, 9525356.25, 9526981.25, 9529025.0, 9529437.5, 9538645.3125, 9541935.9375, 9542082.8125, 9543035.9375, 9543720.3125, 9610778.125, 9619928.125, 9638068.75, 9715906.25, 9717585.9375, 9719315.625, 9741928.125, 9742382.8125, 9743175.0, 9743557.8125, 9744164.0625, 9744270.3125, 9747600.0, 9747651.5625, 9749429.6875, 9751496.875, 9767779.6875, 9773337.5, 9776160.9375, 9786034.375, 9786995.3125, 9791762.5, 9796870.3125, 9796893.75, 9806301.5625, 9811043.75, 9890292.1875, 9906764.0625, 9925570.3125, 9931945.3125, 9965332.8125, 9969200.0, 9969675.0, 9969748.4375, 9973417.1875, 9976837.5, 9988881.25, 10009853.125, 10011873.4375, 10013937.5, 10023450.0, 10023475.0, 10023487.5, 10038245.3125, 10100821.875, 10103193.75, 10104520.3125, 10105098.4375, 10115637.5, 10132003.125, 10133939.0625, 10141028.125, 10171753.125, 10204665.625, 10235901.5625, 10257421.875, 10263312.5, 10297973.4375, 10316348.4375, 10421642.1875, 10515720.3125, 10520773.4375, 10521157.8125, 10557134.375, 10616289.0625, 10664932.8125, 11031510.9375, 11081223.4375, 11087146.875, 11148254.6875, 11150075.0, 11154546.875, 11157140.625, 11157390.625, 11158046.875, 11159251.5625, 11163929.6875, 11172446.875, 11178334.375, 11181218.75, 11181789.0625, 11181923.4375, 11182198.4375, 11185790.625, 11186618.75, 11186893.75, 11187134.375, 11187229.6875, 11187509.375, 11187545.3125, 11187996.875, 11188053.125, 11188093.75, 11188289.0625, 11188298.4375, 11188331.25, 11188475.0, 11188684.375, 11189248.4375, 11189285.9375, 11189503.125, 11189548.4375, 11189592.1875, 11189967.1875, 11189973.4375, 11190028.125, 11190104.6875, 11190146.875, 11190362.5, 11190379.6875, 11190412.5, 11190440.625, 11190889.0625, 11190901.5625, 11190912.5, 11190926.5625, 11190970.3125, 11190984.375, 11191029.6875, 11191092.1875, 11191110.9375, 11191132.8125, 11191142.1875, 11191226.5625, 11191315.625, 11191331.25, 11191434.375, 11191521.875, 11191575.0, 11191585.9375, 11191670.3125, 11191746.875, 11191912.5, 11192279.6875, 11192371.875, 11192506.25, 11192553.125, 11192556.25, 11192593.75, 11192596.875, 11192756.25, 11192851.5625, 11192982.8125, 11193021.875, 11193568.75, 11193670.3125, 11193929.6875, 11195054.6875, 11195354.6875, 11195404.6875, 11195704.6875, 11195765.625, 11195878.125, 11196057.8125, 11196545.3125, 11197056.25, 11197087.5, 11197525.0, 11197757.8125, 11197810.9375, 11198140.625, 11198479.6875, 11198828.125, 11198859.375, 11199081.25, 11199135.9375, 11199156.25, 11199279.6875, 11199400.0, 11199479.6875, 11199698.4375, 11199834.375, 11199957.8125, 11200068.75, 11200262.5, 11200337.5, 11200340.625, 11200615.625, 11200709.375, 11200728.125, 11200754.6875, 11200759.375, 11200873.4375, 11201026.5625, 11201135.9375, 11201300.0, 11201323.4375, 11201425.0, 11201487.5, 11201720.3125, 11201721.875, 11201723.4375, 11201868.75, 11202296.875, 11202446.875, 11202575.0, 11202614.0625, 11202729.6875, 11203282.8125, 11203473.4375, 11203918.75, 11203931.25, 11203984.375, 11204048.4375, 11204375.0, 11204503.125, 11204721.875, 11204970.3125, 11204973.4375, 11205012.5, 11205246.875, 11205748.4375, 11205829.6875, 11205895.3125, 11205915.625, 11206426.5625, 11206657.8125, 11207381.25, 11207428.125, 11207442.1875, 11207800.0, 11208667.1875, 11208676.5625, 11209234.375, 11209435.9375, 11209575.0, 11209787.5, 11210067.1875, 11210945.3125, 11211142.1875, 11211165.625, 11211206.25, 11211392.1875, 11211953.125, 11212381.25, 11212504.6875, 11212795.3125, 11213148.4375, 11213168.75, 11213350.0, 11213417.1875, 11213596.875, 11213603.125, 11214114.0625, 11214151.5625, 11214151.5625, 11214232.8125, 11214279.6875, 11214700.0, 11214840.625, 11214898.4375, 11214903.125, 11214906.25, 11215139.0625, 11215178.125, 11215231.25, 11215243.75, 11215276.5625, 11215678.125, 11215817.1875, 11215976.5625, 11215985.9375, 11216157.8125, 11216309.375, 11216551.5625, 11217090.625, 11217496.875, 11217503.125, 11217675.0, 11217793.75, 11217915.625, 11217970.3125, 11218021.875, 11218479.6875, 11218504.6875, 11218801.5625, 11218946.875, 11219442.1875, 11219942.1875, 11219964.0625, 11220021.875, 11220178.125, 11220557.8125, 11220560.9375, 11220590.625, 11220715.625, 11220875.0, 11220925.0, 11220970.3125, 11221028.125, 11221162.5, 11221514.0625, 11221576.5625, 11221579.6875, 11221593.75, 11221625.0, 11221628.125, 11221654.6875, 11221712.5, 11221957.8125, 11222310.9375, 11222510.9375, 11222565.625, 11222570.3125, 11222685.9375, 11222725.0, 11223121.875, 11223135.9375, 11223337.5, 11223454.6875, 11223609.375, 11223725.0, 11223840.625, 11224075.0, 11224110.9375, 11224146.875, 11224329.6875, 11224376.5625, 11224506.25, 11224570.3125, 11224675.0, 11224732.8125, 11224953.125, 11224957.8125, 11225084.375, 11225321.875, 11225578.125, 11225617.1875, 11225681.25, 11225715.625, 11225753.125, 11225798.4375, 11226046.875, 11226140.625, 11226384.375, 11226423.4375, 11226467.1875, 11226467.1875, 11226482.8125, 11226514.0625, 11226540.625, 11226557.8125, 11226567.1875, 11226576.5625, 11226651.5625, 11226818.75, 11226828.125, 11226860.9375, 11226870.3125, 11226992.1875, 11227156.25, 11227181.25, 11227203.125, 11227292.1875, 11227292.1875, 11227478.125, 11227657.8125, 11227671.875, 11227723.4375, 11227740.625, 11227809.375, 11227812.5, 11227842.1875, 11227895.3125, 11227901.5625, 11227907.8125, 11228212.5, 11228301.5625, 11228368.75, 11228375.0, 11228400.0, 11228468.75, 11228610.9375, 11228725.0, 11228742.1875, 11228815.625, 11228875.0, 11228882.8125, 11229112.5, 11229334.375, 11229376.5625, 11229526.5625, 11229543.75, 11229554.6875, 11229739.0625, 11229798.4375, 11229803.125, 11229870.3125, 11229929.6875, 11229985.9375, 11230042.1875, 11230084.375, 11230106.25, 11230178.125, 11230264.0625, 11230359.375, 11230414.0625, 11230500.0, 11230635.9375, 11230651.5625, 11230670.3125, 11230912.5, 11231665.625, 11231750.0, 11231854.6875, 11231867.1875, 11232121.875, 11232148.4375, 11232245.3125, 11232439.0625, 11232779.6875, 11232881.25, 11233065.625, 11233134.375, 11233137.5, 11233325.0, 11233395.3125, 11233398.4375, 11233868.75, 11234004.6875, 11234018.75, 11234459.375, 11234584.375, 11234676.5625, 11235034.375, 11235103.125, 11235625.0, 11235676.5625, 11235748.4375, 11235876.5625, 11235884.375, 11236228.125, 11236296.875, 11236773.4375, 11237193.75, 11237264.0625, 11237368.75, 11237385.9375, 11237579.6875, 11237826.5625, 11238039.0625, 11238207.8125, 11238731.25, 11238868.75, 11239037.5, 11239109.375, 11239160.9375, 11239621.875, 11240298.4375, 11240473.4375, 11240667.1875, 11241101.5625, 11241221.875, 11241307.8125, 11241346.875, 11241737.5, 11241751.5625, 11241756.25, 11241821.875, 11241904.6875, 11241917.1875, 11241943.75, 11242012.5, 11242014.0625, 11242056.25, 11242626.5625, 11242646.875, 11242710.9375, 11243243.75, 11243496.875, 11243656.25, 11243726.5625, 11244084.375, 11244337.5, 11245167.1875, 11245357.8125, 11245587.5, 11245631.25, 11245703.125, 11245731.25, 11245851.5625, 11245975.0, 11246042.1875, 11246170.3125, 11246182.8125, 11246904.6875, 11247129.6875, 11247309.375, 11247362.5, 11247471.875, 11247734.375, 11247900.0, 11248231.25, 11248239.0625, 11248315.625, 11248387.5, 11248528.125, 11248682.8125, 11248753.125, 11248785.9375, 11248892.1875, 11249142.1875, 11249175.0, 11249695.3125, 11249734.375, 11250970.3125, 11251001.5625, 11251384.375, 11251653.125, 11251717.1875, 11251846.875, 11252343.75, 11252778.125, 11253340.625, 11254192.1875, 11254642.1875, 11254710.9375, 11254740.625, 11254776.5625, 11254784.375, 11255078.125, 11255620.3125, 11255815.625, 11255928.125, 11256446.875, 11256446.875, 11256614.0625, 11257279.6875, 11257450.0, 11258204.6875, 11258609.375, 11258634.375, 11259451.5625, 11259792.1875, 11260662.5, 11260715.625, 11260859.375, 11261665.625, 11262460.9375, 11262564.0625, 11262565.625, 11262625.0, 11264050.0, 11264137.5, 11265189.0625, 11265687.5, 11265751.5625, 11265840.625, 11265943.75, 11266532.8125, 11267440.625, 11267560.9375, 11267832.8125, 11268154.6875, 11269114.0625, 11269442.1875, 11269500.0, 11269903.125, 11269984.375, 11270670.3125, 11271157.8125, 11271178.125, 11271323.4375, 11271457.8125, 11271507.8125, 11271621.875, 11271768.75, 11271820.3125, 11271985.9375, 11272401.5625, 11272409.375, 11272700.0, 11272768.75, 11272804.6875, 11272954.6875, 11273239.0625, 11273246.875, 11273371.875, 11273478.125, 11273490.625, 11273871.875, 11274337.5, 11274381.25, 11274671.875, 11274860.9375, 11274910.9375, 11275226.5625, 11275287.5, 11275298.4375, 11275357.8125, 11275526.5625, 11275825.0, 11276135.9375, 11276182.8125, 11276209.375, 11276214.0625, 11276262.5, 11276278.125, 11276431.25, 11276592.1875, 11276671.875, 11276798.4375, 11276882.8125, 11277109.375, 11277114.0625, 11277500.0, 11277615.625, 11277729.6875, 11277864.0625, 11278112.5, 11278159.375, 11278257.8125, 11278439.0625, 11278887.5, 11278912.5, 11278935.9375, 11278951.5625, 11279217.1875, 11279303.125, 11279500.0, 11279790.625, 11280242.1875, 11280390.625, 11280829.6875, 11280932.8125, 11280975.0, 11281148.4375, 11282146.875, 11282423.4375, 11282645.3125, 11283068.75, 11283117.1875, 11283446.875, 11283675.0, 11285146.875, 11285829.6875, 11286868.75, 11286884.375, 11287160.9375, 11287489.0625, 11288742.1875, 11289150.0, 11289871.875, 11290551.5625, 11290565.625, 11290871.875, 11292012.5, 11293326.5625, 11293557.8125, 11294004.6875, 11294284.375, 11297067.1875, 11298284.375, 11301239.0625, 11307773.4375], [19.2154295523539, 89.55258854213393, 14.370355298661053, 57.20399228797166, 44.811780504814074, 26.75388049298675, 6.530170367706963, 46.51528292778501, 5.291261754564691, 11.241879054981549, 24.51955003314039, 5.4735318475161385, 16.98705148080214, 6.854651163154591, 134.54444194779637, 108.56779278433693, 29.617489241934027, 51.53468854959295, 45.763343383599434, 29.438834060151855, 8.984525156225294, 11.413458231150608, 79.11277376253217, 29.332268207312964, 7.754292493196194, 5.1596478858584565, 65.75759256528346, 29.886755584893912, 26.5584010269535, 23.863277590325925, 14.547172927251028, 60.30310941298213, 19.09707671046739, 156.3456596865374, 17.16082850848504, 46.84886998523748, 22.061472761700315, 21.5868778827754, 9.24097056919769, 23.262784910507882, 118.25907569736367, 61.11357374615261, 38.57253171293951, 47.46288157118514, 34.99858159090979, 45.693231652771615, 9.251113344025454, 8.059526900380703, 90.49611075964678, 20.520631561015886, 11.720049670994765, 29.740928249506485, 107.58167149662388, 6.40105327770997, 28.50712426756071, 13.07560101312344, 68.26889573018359, 17.74258765455825, 56.945125949754264, 72.29505848984574, 6.020204843413952, 44.54701547359016, 40.36894174918001, 17.067673224056758, 9.605313045803141, 18.59978060841521, 57.41665724064427, 217.70902154594663, 60.42647272904755, 64.01513099818922, 5.657883631660489, 69.18948929366888, 12.09263699146219, 13.107291736562065, 39.879337497829894, 22.24404132702483, 23.154982776707463, 123.23317597273478, 53.05712160521556, 25.354345409682807, 6.980673241152645, 30.94223972389825, 10.214020951023041, 142.33953166532407, 297.29521644971095, 10.760884240449247, 23.924188692267542, 145.2495245410926, 7.913224832692245, 70.96656852421646, 30.52159985969437, 99.8772153032514, 49.93174090042638, 30.353901470667097, 9.26361765387983, 21.240483298665268, 9.544389337823151, 33.41682334052246, 124.01046154983763, 31.81760139138788, 97.40791875541008, 31.90252501810467, 26.154246510660762, 5.140815620989394, 17.177225782037713, 64.50741966314932, 35.957545738789264, 104.5144784399327, 26.64889928125563, 59.48767625255684, 16.942852235320416, 134.40837950624547, 22.053034040793786, 6.23291908819242, 5.911287509029278, 110.56842971665007, 7.677522869714147, 6.284719157619654, 42.683435128956695, 8.32042025593095, 9.928282408649084, 19.722045606635636, 28.987548108221837, 22.541174969365414, 100.96387202082396, 5.803734727703886, 19.65109117609528, 24.842907083337348, 30.680043183427458, 8.658986238020574, 26.621158222116563, 35.042213976536594, 56.984986104709336, 22.89199804998787, 12.055277195586118, 67.62044480693226, 16.47506936148762, 31.985816272631364, 20.16372603464425, 25.89806132220823, 28.880419122242554, 15.46104680811437, 21.045317789123715, 45.491566144160146, 85.2147434472077, 17.309245418570566, 8.078382959682134, 8.195185905875588, 20.622677176985167, 6.1076938727813745, 137.97297094019163, 40.72741328583843, 105.8852169166256, 59.00638138940108, 37.84499218295711, 18.052755338778198, 12.777847454904615, 18.981400770137018, 166.2523480868819, 119.7609513969797, 114.13651494861878, 16.652970221240064, 46.62280759256032, 195.60053890605388, 46.59376171522019, 84.21447092167637, 35.67032040601055, 54.621348313428555, 22.53155975251148, 43.436918197163806, 15.338022649426478, 38.304689734371735, 15.100770425629005, 44.81912929448974, 6.320434263713759, 28.76550382602135, 107.88435606427258, 8.450202504917593, 84.97320765353425, 37.56270282959764, 58.030205410898795, 5.7709920107508434, 32.06595252128765, 17.993235980718183, 17.6253432354124, 41.00399468709798, 59.93992228451643, 84.16686319144326, 10.351134411316657, 13.5293570573456, 49.50857954099376, 34.1751830579345, 63.41813705918761, 10.368743511179968, 18.738693240509853, 6.657544770175715, 52.444248361771116, 49.51492145254019, 24.429689187005337, 12.706705068052425, 53.697597530578605, 8.236528576630326, 157.20851863249132, 34.22233017568783, 24.64701976296206, 47.981146511860544, 56.64239394103669, 5.7529193433668295, 21.83882126418456, 19.57399568647594, 32.83299838602224, 298.4792889812215, 98.62888764356468, 26.47639052807266, 55.785696132754126, 19.240462877606962, 18.30042911850452, 55.666689248636914, 6.032081601307642, 103.43258385604304, 7.426222274868422, 175.87626436116688, 5.208550669803976, 87.66195913679998, 6.970929251076096, 101.57004136135618, 51.365073943841615, 60.9338213153302, 52.62441405820311, 12.093841119121558, 79.93321541879693, 131.19417342765024, 8.139494005406712, 5.821175516598972, 8.455144136445806, 70.4630275673035, 34.323304711144594, 84.58534687972843, 108.64944103920774, 13.681570178809102, 62.22614937523675, 12.365094712858127, 28.492406438982332, 11.695273511376318, 52.46457819806943, 80.60044291018835, 66.2486926663718, 74.93155719876815, 5.887282941632003, 18.54482833822992, 10.659805712491513, 106.2018419397337, 78.7222965759835, 5.176864716447867, 47.39263194185077, 58.657984790067154, 12.137771994305721, 5.6537688156527866, 57.06383264487278, 15.620701400616198, 48.72667263087499, 6.868133893677933, 7.624943522082688, 6.511652217974223, 28.12847202639719, 6.3308568389160085, 85.81788067101157, 39.5355884853977, 25.799927371278617, 5.635185781668891, 19.413694709855502, 103.66360819193886, 39.27383368194472, 13.696639081188916, 44.29238541410025, 110.26364335749889, 20.176164818304052, 11.964093289182197, 48.53086118336387, 49.204634448661025, 6.970078167026922, 8.935250347429678, 47.064436161810555, 98.01106467781845, 13.380820703223215, 35.649780551146065, 8.416738377800952, 6.425517229849981, 85.7151444578088, 12.068191739864806, 10.415593125635109, 8.972398263117928, 42.37344961171816, 58.2661475358052, 23.319547183968105, 73.68292085277879, 80.81713775262017, 38.13017052214444, 7.693987769674597, 16.377388457263514, 86.45912336994074, 93.21326374404448, 52.01944740977451, 18.291413416405508, 7.514705503132166, 40.84015074842837, 17.22129104662929, 6.398440824001257, 20.405598255080214, 25.33041735696217, 13.445138437343513, 92.95164637472942, 41.36511519869156, 65.44956925553036, 55.155084595589955, 30.00283269197329, 19.612254510784922, 32.83023133034581, 44.758132906000625, 15.308796375793193, 12.413608308115478, 10.68845631053543, 9.245130215340437, 16.671255064750632, 37.374326254738456, 5.777629518547999, 35.04937588090091, 17.980773221899423, 33.13911708660834, 15.058198397742165, 76.7816122129089, 6.377286346592362, 44.489582365049124, 45.18491444488918, 10.724624285581632, 14.070434476963813, 45.13512991270552, 12.609869601272743, 10.595083751588458, 20.89065575264003, 48.26806553866681, 36.10995743514956, 7.7460302051703, 7.505432818216993, 25.27096680329298, 5.128056347923415, 8.081620682470584, 16.787434254518708, 12.548672309875796, 7.309578165500295, 29.020397725107284, 15.819777974452645, 21.842558275844425, 42.749059836789236, 56.04955688135264, 18.431463067382307, 41.86821160239721, 24.012749600134843, 81.63171093430608, 27.620072308716992, 9.997263199271416, 7.754638245339949, 60.24650286862469, 23.835065228706924, 7.401011322532126, 15.703710682886056, 7.13657384749423, 13.657442972863167, 13.416787128507428, 31.59162517362728, 39.78198164977961, 15.291569046090522, 21.48855460426641, 146.21183152002635, 59.78011245708154, 70.81183821481665, 141.7277492793319, 26.403224294035706, 39.51417239946688, 30.897482814382986, 17.552609536530802, 9.58766721154342, 77.56996682126315, 17.147744434704386, 31.058917482281707, 11.290727540660633, 20.224815672530422, 65.82104401212247, 72.92231591405691, 6.585060122686211, 11.125336338174773, 5.866850039796829, 16.118845809924984, 63.95025140890468, 67.25577173679876, 13.384798335304454, 82.94665171031781, 35.619306463418056, 23.205619874504407, 74.20522964493901, 7.340666706822463, 14.749773076930264, 9.271464809005398, 46.48912410617532, 7.948426370465435, 23.756361712600018, 6.088395514391004, 6.302727603136325, 71.53022670204054, 35.021021869095875, 17.623555629125867, 14.411791763744088, 13.628827912187504, 5.404040694168145, 17.53156105658011, 19.99823029524813, 7.754790468958489, 89.70109119222542, 7.556148365213274, 25.965060177981975, 23.426498016563794, 10.14428547043145, 57.42988902025119, 31.815525977940702, 7.743952994812329, 32.109797045428806, 53.94240851334756, 46.25298270850679, 26.772815264614813, 39.50925386853139, 23.77816653558072, 69.79518841739711, 10.550828769485605, 49.40105418719045, 43.02799056819079, 5.908332187910523, 8.389694091899239, 5.411725324815218, 6.659641162254998, 5.148741916051411, 13.018939389113124, 14.369263879145777, 21.146769197012503, 21.872154900812458, 13.688695176942094, 8.448224543756433, 37.24216527716307, 16.28188339153168, 17.928664213940998, 11.430343282894036, 44.70953331731848, 52.62467879165811, 45.493779866664454, 40.00077846891103, 54.48894618210974, 11.324782041286394, 18.9468971704655, 15.624337007099601, 10.07164556284105, 36.136225125879946, 6.383005547928466, 14.901742383627454, 14.925043733038475, 22.620977183698667, 6.983810820614679, 24.85458182474022, 12.64207354586782, 8.959333766573748, 14.072681049693339, 14.727675607944024, 5.897178780864305, 14.793275357159445, 41.300704270154355, 40.57628120280969, 12.853708620908415, 13.453643611525175, 30.722172602768957, 9.926232490786155, 47.14677218625363, 11.911645024717282, 35.48180106788522, 69.73056913729127, 26.493920191066046, 91.13206307859123, 26.142266212372864, 30.708904276228445, 20.39439571720566, 23.20537404952701, 42.49547458507458, 25.64729266426884, 34.1197785723608, 20.262190240739343, 5.189546737127401, 12.572815226237152, 9.430343393089538, 7.0183345759590745, 48.09402447066543, 6.628072789602321, 9.908964510692961, 13.518875043548238, 37.97285072811504, 39.549782492630676, 9.002334403506342, 22.00836846897292, 49.1403915087156, 6.1593746587870095, 20.334718043361036, 66.37898921317844, 87.91118950629291, 26.38949628689685, 75.24588157355609, 35.69825321773843, 14.673983140228355, 14.82047446937802, 72.71234538104842, 29.86007877256362, 51.52256662333684, 29.625169596203836, 95.2098589382308, 10.266281834717775, 49.506976388462604, 22.206844028616224, 10.483578363836186, 15.082327217535068, 42.24666349694151, 51.87424614935948, 51.48758727009644, 19.646884891683953, 5.293078022010689, 11.093344208484977, 5.8194539965376055, 6.939367469405111, 8.352271377889013, 6.284496889469913, 16.104931226655747, 72.69328131479669, 11.427941787757222, 36.74397531453554, 10.347121196824602, 21.014102576829373, 11.505990441010656, 11.939822363503923, 40.83490349167668, 5.728664673603575, 29.26027164895753, 12.51418926142831, 22.6775898377942, 6.079412548335617, 29.302165490976535, 102.70976793125672, 56.84214704002435, 66.2668457954252, 96.72747174152747, 67.9405082300604, 20.27287876012243, 36.542178314746145, 6.151914691269665, 41.244366789223335, 14.441408443638895, 27.744506314230723, 26.239260501988262, 32.229118398702624, 17.697433462978232, 14.39602176586904, 5.231643332867902, 22.082874885854054, 5.402808866845115, 25.878025141512442, 6.040398486905212, 24.636275886575717, 71.04539707230435, 58.04083447124795, 26.003062357681554, 14.4022288742417, 60.91498228126416, 19.45826261863541, 10.123778425960948, 23.445382235071065, 26.434888263639614, 5.612731546844431, 31.64826045373723, 33.16981878834403, 66.95048665925208, 42.33817539051895, 70.92840005186784, 30.607143968110407, 44.112125804838236, 60.86494062425038, 88.76341085264215, 79.31061842136104, 8.908223246334472, 5.223142510423527, 76.9641490504155, 47.7765082859956, 25.943899801221846, 36.53781248430798, 10.5910286757399, 23.748638609542652, 18.089933854231074, 61.70399034704518, 17.456940543601768, 13.18213828624253, 32.6484229293632, 32.43281571002123, 101.61661638600606, 57.345436713391315, 28.68309170424345, 5.2186775478605, 5.681630259765935, 79.74819338676163, 7.451958603634344, 57.63097491347129, 12.277893560714832, 19.52871576319681, 136.8781468931789, 6.037513420892447, 38.47534471979837, 12.983316575943968, 32.907006897292824, 9.578760087295343, 13.116074587375607, 24.775521522643736, 10.755654488319447, 99.96875360790315, 12.85317541431894, 36.782147607281594, 49.43601960509329, 5.220232653968736, 5.457796766727303, 9.43810473661178, 10.067959735892794, 20.204842789321127, 25.958579603808438, 25.722918264879215, 59.08686006389523, 21.52328120958489, 23.031003790454857, 20.876197297286414, 75.99798096131326, 27.09070877294054, 53.797213417670186, 43.047582219957135, 67.87573448830743, 5.745826318142593, 31.418046575271, 47.14878914859085, 88.18478941213253, 5.280473165712254, 15.48063081445228, 5.046605275281482, 7.080052212698353, 60.49974931357806, 6.233840845351417, 186.37660783624187, 44.38737537705964, 39.542584946948445, 54.17550854876919, 6.1934164678063075, 49.28731723008159, 32.38719736798186, 13.160990412669836, 20.857098742305503, 10.29534003529524, 21.745973892790097, 65.47553683411499, 24.142971358875677, 24.731588505718477, 42.32035782285217, 15.249931475779778, 61.86547253028861, 75.76114057251098, 24.65683997342812, 16.17679519295178, 5.424220788653245, 18.696611551758302, 31.427751664361818, 44.59367498670005, 6.338670706173984, 44.25557778986483, 24.021748303272865, 49.89078313695396, 37.01252132827999, 8.923371271224662, 56.207632262218226, 16.490000663356724, 7.684860379061482, 31.537362101845687, 5.593332306509236, 20.07892381470297, 5.256933266731531, 17.066383251529615, 21.4851241315282, 62.75012960561685, 34.878283989950035, 112.29818462854801, 14.517489635698217, 55.35497554948135, 46.61620988033934, 7.995595640237112, 16.367726420642015, 14.633649664447823, 63.244302879746726, 6.323392369985584, 33.083540173689094, 21.04261395869716, 41.78009336545705, 21.754887736277286, 10.652181316482707, 24.86035827570537, 6.611075669615951, 5.486563150483383, 27.136199183780025, 10.339596730452229, 17.89746605688413, 19.63441586536146, 46.02126251324685, 69.83957315389293, 10.38244692954969, 47.937329206654624, 6.357993353579579, 30.606312484337685, 121.39901888296282, 31.131664872198336, 63.713611657327334, 25.21776310494605, 7.627902233338459, 10.657564454179456, 28.642218443231357, 11.613065391992013, 42.34442597735319, 35.32740052263784, 37.40350165919267, 50.29054156536031, 38.10146008129219, 6.1630078659859695, 47.81007741958481, 14.015993735012309, 13.535613551136759, 12.222022180438449, 21.00914396156214, 5.255335921072931, 23.552704328368158, 61.63082233086503, 34.883074170987534, 111.0614040977351, 54.28632388426466, 21.16683783679268, 10.82600413512975, 13.574385218876655, 36.11431673752004, 34.434675240142354, 21.166802252015376, 42.253559781072866, 48.93302433198453, 49.0087195292572, 5.094514081473814, 53.70087328743548, 14.055154591890712, 9.400016264302172, 85.26016705524454, 26.133648807929344, 27.85875123119011, 18.721963963624365, 27.04448873474811, 18.551434270301595, 88.10343861551631, 5.172812180767204, 7.536165016560863, 84.44553839054439, 94.33528250898142, 14.542033743820722, 41.550086840018096, 31.181586621794544, 6.346373239572135, 14.007272419614747, 22.88840235289654, 84.99552667987346, 8.923836385514063, 6.525574964902848, 57.592678054651316, 27.032001811327987, 16.936495837460107, 21.58473914165049, 49.0053320072829, 6.476124317277085, 7.6381218063281855, 50.22313504984391, 23.68894001012083, 10.058284265745536, 21.349926314564343, 86.82137874304998, 93.30187717605122, 22.457234697770264, 57.364868405822314, 7.2176461450538225, 8.356776482822138, 35.53633266949107, 9.123054006519741, 53.326938448449226, 68.20658738474245, 5.794823176706958, 7.392507385370459, 10.039751759120962, 8.464705776114249, 110.1473650536748, 25.12321833681538, 17.041871140599483, 10.207849835930586, 17.52927864815489, 14.882854672939615, 103.74907638309892, 42.436508502657674, 42.20014002723499, 86.9821261483078, 45.190905517084616, 15.226006590605854, 5.162098224343081, 18.05866892224462, 11.65938425894094, 5.23054428556187, 81.92084227205126, 131.48054999218957, 100.2094161952809, 20.443788902006748, 11.148196921283928, 9.42054530875808, 24.27798270549415, 22.479825784191405, 19.724411713573843, 31.035190944182688, 6.574803553942991, 46.662117481836496, 5.20269142198402, 31.689227212255634, 25.90241841883148, 11.126566503321005, 21.87362598558279, 16.451087858115578, 35.35891668003922, 80.75855591424877, 44.207469370612145, 5.414109078510478, 18.427740172277627, 6.946598965966919, 45.90815609864583, 12.734407656729031, 26.06043557087196, 11.96663407508298, 14.279078972841074, 24.26501670706422, 7.03381356372197, 28.901460383439407, 96.98107554255881, 101.47997996304743, 17.901811991650305, 14.37263516794507, 51.314730012017876, 20.51014678250479, 14.465450029298346, 14.981327034442483, 29.03062224571492, 98.26969569482813, 5.266282271203422, 43.33310546741839, 11.109202027010381, 7.422884720861018, 15.797596666543068, 36.472231883671995, 9.192369551553885, 17.831687671194416, 108.38294338321113, 5.806484111987648, 6.13595417273739, 20.202230120484447, 14.425503779032276])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)