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 = 44422
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);
([3087773.4375, 3279492.1875, 3291940.625, 3302131.25, 3332648.4375, 3340603.125, 3341659.375, 3341721.875, 3351781.25, 3373592.1875, 3376653.125, 3379860.9375, 3384843.75, 3396104.6875, 3443282.8125, 3445862.5, 3449479.6875, 3450982.8125, 3453625.0, 3462920.3125, 3480990.625, 3481003.125, 3506332.8125, 3512343.75, 3528484.375, 3547621.875, 3550703.125, 3594795.3125, 3598578.125, 3617646.875, 3662160.9375, 3691170.3125, 3697804.6875, 3703453.125, 3707943.75, 3710864.0625, 3711143.75, 3726654.6875, 3728825.0, 3739898.4375, 3759232.8125, 3776081.25, 3785918.75, 3829915.625, 3838668.75, 3854609.375, 3869617.1875, 3903290.625, 3932415.625, 3938745.3125, 3945310.9375, 3949546.875, 3951993.75, 3957442.1875, 3958454.6875, 3958465.625, 3959054.6875, 3960396.875, 3960943.75, 3961643.75, 3962089.0625, 3963312.5, 3963360.9375, 3963529.6875, 3964828.125, 3967365.625, 3970403.125, 3970807.8125, 3972093.75, 3972198.4375, 3972228.125, 3977357.8125, 3987478.125, 3988395.3125, 3988934.375, 3990846.875, 3994178.125, 3996070.3125, 4001526.5625, 4001946.875, 4002076.5625, 4002365.625, 4003300.0, 4004778.125, 4014709.375, 4014735.9375, 4016078.125, 4017732.8125, 4027700.0, 4041681.25, 4046200.0, 4055504.6875, 4058004.6875, 4066754.6875, 4071087.5, 4083629.6875, 4099378.125, 4101448.4375, 4170231.25, 4172356.25, 4178989.0625, 4187404.6875, 4233329.6875, 4237554.6875, 4257335.9375, 4258609.375, 4260948.4375, 4286045.3125, 4286953.125, 4287385.9375, 4287959.375, 4293464.0625, 4293470.3125, 4298512.5, 4299073.4375, 4299826.5625, 4302045.3125, 4302190.625, 4303270.3125, 4306889.0625, 4313354.6875, 4313610.9375, 4316346.875, 4365762.5, 4367503.125, 4372154.6875, 4372220.3125, 4372660.9375, 4372743.75, 4376043.75, 4382392.1875, 4382400.0, 4402967.1875, 4408334.375, 4411962.5, 4437851.5625, 4439017.1875, 4440671.875, 4448135.9375, 4460140.625, 4462029.6875, 4468298.4375, 4475170.3125, 4476517.1875, 4478134.375, 4487668.75, 4493178.125, 4504848.4375, 4510612.5, 4514323.4375, 4519190.625, 4523985.9375, 4524639.0625, 4535112.5, 4536723.4375, 4551214.0625, 4587389.0625, 4592243.75, 4603082.8125, 4604612.5, 4611210.9375, 4611635.9375, 4611650.0, 4612910.9375, 4614756.25, 4621246.875, 4626018.75, 4648462.5, 4661704.6875, 4670635.9375, 4677840.625, 4685984.375, 4689985.9375, 4703406.25, 4706595.3125, 4726128.125, 4727179.6875, 4729268.75, 4731562.5, 4741643.75, 4741659.375, 4746734.375, 4751593.75, 4754264.0625, 4761575.0, 4767590.625, 4770246.875, 4772746.875, 4777084.375, 4777148.4375, 4781417.1875, 4784057.8125, 4784570.3125, 4786085.9375, 4788357.8125, 4791876.5625, 4792198.4375, 4799470.3125, 4806020.3125, 4818323.4375, 4822328.125, 4823754.6875, 4825445.3125, 4826709.375, 4836514.0625, 4852610.9375, 4861657.8125, 4864784.375, 4870968.75, 4881329.6875, 4881331.25, 4882671.875, 4887182.8125, 4887264.0625, 4887571.875, 4942845.3125, 5202618.75, 5267948.4375, 5304584.375, 5311142.1875, 5431157.8125, 5439004.6875, 5566810.9375, 5613315.625, 5617290.625, 5661426.5625, 5729892.1875, 5777900.0, 5788240.625, 5855723.4375, 5865017.1875, 5885287.5, 5913842.1875, 5945917.1875, 5957729.6875, 5959490.625, 5963306.25, 5968056.25, 5997043.75, 6033198.4375, 6166170.3125, 6226726.5625, 6230910.9375, 6272492.1875, 6276678.125, 6350546.875, 6384792.1875, 6414335.9375, 6415671.875, 6415682.8125, 6444023.4375, 6451354.6875, 6453412.5, 6455237.5, 6455496.875, 6455756.25, 6455810.9375, 6455815.625, 6456259.375, 6456390.625, 6457600.0, 6457610.9375, 6457835.9375, 6458114.0625, 6458634.375, 6458912.5, 6459537.5, 6460218.75, 6460381.25, 6464935.9375, 6482834.375, 6501218.75, 6567126.5625, 6575485.9375, 6575707.8125, 6577029.6875, 6577520.3125, 6580496.875, 6580612.5, 6582296.875, 6595951.5625, 6622120.3125, 6626510.9375, 6626848.4375, 6628718.75, 6628795.3125, 6630976.5625, 6632500.0, 6634203.125, 6635343.75, 6635479.6875, 6636526.5625, 6637414.0625, 6694087.5, 6694162.5, 6727965.625, 6727979.6875, 6728510.9375, 6768018.75, 6769990.625, 6819342.1875, 6824698.4375, 6827067.1875, 6828110.9375, 6920218.75, 6920220.3125, 6924500.0, 6955023.4375, 6955450.0, 6957070.3125, 6984871.875, 6999907.8125, 7009839.0625, 7115892.1875, 7120931.25, 7171362.5, 7177581.25, 7212764.0625, 7233709.375, 7291923.4375, 7337326.5625, 7337984.375, 7358348.4375, 7370460.9375, 7405179.6875, 7479779.6875, 7542878.125, 7546215.625, 7546226.5625, 7579470.3125, 7590853.125, 7688746.875, 7738551.5625, 7801984.375, 7835351.5625, 7835420.3125, 7908001.5625, 7959790.625, 7999307.8125, 8046478.125, 8080765.625, 8094612.5, 8099781.25, 8123092.1875, 8137417.1875, 8251940.625, 8302825.0, 8340545.3125, 8367860.9375, 8370220.3125, 8449735.9375, 8458984.375, 8459814.0625, 8467187.5, 8469817.1875, 8488120.3125, 8493984.375, 8497235.9375, 8498331.25, 8501132.8125, 8504198.4375, 8579548.4375, 8644642.1875, 8657851.5625, 8722443.75, 8791884.375, 8810348.4375, 8823653.125, 8966678.125, 9049198.4375, 9186187.5, 9296284.375, 9325573.4375, 9330907.8125, 9357398.4375, 9360792.1875, 9361407.8125, 9447243.75, 9448639.0625, 9449212.5, 9449701.5625, 9462470.3125, 9489589.0625, 9491907.8125, 9540006.25, 9636868.75, 9727014.0625, 9732673.4375, 9732914.0625, 9735464.0625, 9742296.875, 9744396.875, 9745239.0625, 9803107.8125, 9900350.0, 9908162.5, 9910843.75, 9948620.3125, 9949189.0625, 9951018.75, 9951565.625, 9954076.5625, 9954132.8125, 9975871.875, 9977475.0, 9977601.5625, 9978351.5625, 9978410.9375, 9978426.5625, 9978746.875, 9978934.375, 9979885.9375, 9980848.4375, 9982334.375, 9983507.8125, 9984860.9375, 9984939.0625, 9984970.3125, 9991339.0625, 10018678.125, 10020339.0625, 10020495.3125, 10020675.0, 10021337.5, 10021643.75, 10021660.9375, 10022270.3125, 10022396.875, 10022862.5, 10023589.0625, 10024279.6875, 10024595.3125, 10026085.9375, 10028837.5, 10095937.5, 10107328.125, 10108204.6875, 10127207.8125, 10133090.625, 10133110.9375, 10149882.8125, 10155068.75, 10155465.625, 10159973.4375, 10160479.6875, 10160507.8125, 10175629.6875, 10187621.875, 10187737.5, 10205395.3125, 10205668.75, 10268126.5625, 10268142.1875, 11172354.6875, 11177729.6875, 11177839.0625, 11178792.1875, 11179850.0, 11180498.4375, 11180504.6875, 11183793.75, 11186384.375, 11187345.3125, 11190289.0625, 11190915.625, 11193473.4375, 11195215.625, 11202617.1875, 11206462.5, 11218329.6875, 11295078.125, 11317231.25, 11317300.0, 11318243.75, 11318614.0625, 11319217.1875, 11319470.3125, 11319945.3125, 11320190.625, 11321784.375, 11321896.875, 11321900.0, 11321907.8125, 11322157.8125, 11324617.1875, 11325123.4375, 11325620.3125, 11325760.9375, 11325940.625, 11326442.1875, 11326460.9375, 11326703.125, 11326810.9375, 11327181.25, 11327264.0625, 11327660.9375, 11327950.0, 11328082.8125, 11328645.3125, 11328690.625, 11329137.5, 11329428.125, 11329448.4375, 11329817.1875, 11329920.3125, 11330070.3125, 11331375.0, 11331465.625, 11331715.625, 11331851.5625, 11331929.6875, 11331940.625, 11332281.25, 11332867.1875, 11333100.0, 11333157.8125, 11334282.8125, 11334704.6875, 11334878.125, 11335160.9375, 11335171.875, 11335220.3125, 11335862.5, 11336420.3125, 11336692.1875, 11337010.9375, 11337740.625, 11337904.6875, 11337912.5, 11338114.0625, 11338264.0625, 11338621.875, 11338668.75, 11339717.1875, 11339978.125, 11340881.25, 11341218.75, 11341492.1875, 11341568.75, 11342414.0625, 11342489.0625, 11342684.375, 11342929.6875, 11342934.375, 11343203.125, 11344056.25, 11344065.625, 11344345.3125, 11346003.125, 11346782.8125, 11347150.0, 11347812.5, 11347867.1875, 11347870.3125, 11347881.25, 11348556.25, 11349464.0625, 11349545.3125, 11349548.4375, 11349603.125, 11350007.8125, 11350568.75, 11350673.4375, 11351117.1875, 11351239.0625, 11351496.875, 11352053.125, 11352648.4375, 11352682.8125, 11353907.8125, 11354253.125, 11355803.125, 11355898.4375, 11356171.875, 11356864.0625, 11357542.1875, 11358006.25, 11358673.4375, 11359612.5, 11359800.0, 11360937.5, 11361762.5, 11362464.0625, 11362490.625, 11362646.875, 11363053.125, 11363376.5625, 11363453.125, 11365003.125, 11365451.5625, 11365660.9375, 11365776.5625, 11366259.375, 11367289.0625, 11367437.5, 11367587.5, 11367898.4375, 11368028.125, 11368704.6875, 11369042.1875, 11369257.8125, 11369528.125, 11369960.9375, 11370104.6875, 11370121.875, 11370282.8125, 11370321.875, 11370715.625, 11370957.8125, 11371170.3125, 11371384.375, 11371710.9375, 11372156.25, 11372167.1875, 11372243.75, 11372289.0625, 11372296.875, 11372796.875, 11374464.0625, 11375037.5, 11375103.125, 11375614.0625, 11375701.5625, 11375831.25, 11376528.125, 11376701.5625, 11376862.5, 11377450.0, 11377623.4375, 11377670.3125, 11379512.5, 11380315.625, 11381628.125, 11382507.8125, 11390015.625, 11393490.625, 11400914.0625, 11412100.0], [9.021341576318875, 10.477288474122311, 57.54232546967111, 43.36223247797474, 11.30619164980767, 41.089963124880875, 8.166058461384358, 5.240813951353239, 42.11975825407877, 5.648981120762856, 38.44809445452768, 34.398871321035394, 33.515944875368504, 9.225651296322075, 5.637175833110015, 12.453412421269853, 9.95662535622695, 42.85192698978631, 58.44927274427586, 71.2157030505174, 67.42362637935705, 7.135961360725214, 79.46549018434845, 9.849819364019993, 25.044360851243567, 71.60105463631187, 33.04248052611274, 65.5889921758571, 60.85461379432205, 13.304759103373232, 9.179640706374053, 60.64724193305344, 19.44004131969818, 5.636328907076851, 24.848293754721837, 5.290449908872773, 5.171263084338775, 7.995774964460886, 75.76501538774251, 24.19017595839289, 60.14228394731865, 29.7064543888022, 7.500352011449034, 10.750330995411119, 13.826946834335155, 33.447354321558436, 11.94859398986753, 75.03945225597064, 8.635126547536887, 52.27764754239169, 67.9266378567987, 46.04808813836873, 101.13580485579475, 123.18776831456105, 22.60207153664546, 60.556266716845286, 15.298635025912702, 85.21712537952381, 28.611818329301595, 14.760795640073713, 19.371264829730862, 51.332686155110764, 12.868475116809758, 17.07158809025163, 6.608998110940824, 81.28311107751189, 7.830052982107149, 64.65233504921224, 16.569579362659873, 16.702030866968922, 9.74260796417946, 30.01294975813429, 92.53561265167383, 19.55944798532267, 24.604704762767536, 165.6791276678834, 12.085598395197227, 11.811065604939285, 66.74262162672616, 33.35766781848508, 10.343889980390738, 51.53812635793888, 5.185852330044586, 33.55835004815951, 57.065205701517044, 21.973078643081088, 31.520777244848492, 8.564235694264939, 7.92196593358585, 9.706312248906553, 49.55153439496237, 32.43307708385775, 27.289044132673123, 32.09290456982849, 28.004534571688644, 13.649655348425629, 139.06591816922727, 7.937612684660185, 18.715945570827422, 135.6264237347823, 13.451375773570465, 65.99526794183818, 93.94145740681168, 132.62053297411637, 73.88730511353756, 9.555031159786623, 10.159484030408336, 28.753230155683998, 56.06243456487948, 6.88037298906764, 9.746176409734632, 76.48351885828143, 11.677681085553473, 7.8734026627145, 83.13458928723722, 12.47601459755017, 56.369336043489085, 38.79422547130149, 9.747888887236929, 35.18213256046644, 17.311572216100345, 6.654798398315644, 21.205006035962363, 146.73375709112386, 29.41553357409855, 14.77133786019102, 102.74532226215821, 10.428520389410174, 5.174315570063226, 23.882563835075153, 20.85725712502494, 18.230827555127817, 37.011115066793565, 24.077082261451345, 14.445864546598049, 121.86511622824248, 5.383008799121236, 10.342308830756494, 38.017352474591966, 5.153762521561403, 36.28312707571996, 104.38080187740073, 19.540276420624945, 8.280307167487964, 26.959179897435153, 57.23915855070114, 81.53097542794902, 16.041392747083027, 31.332282256928348, 17.116063325385987, 82.89291172723831, 50.942663206557036, 61.16509563160895, 13.241508996154332, 92.20384760567592, 84.32115200732882, 7.166571345499603, 32.25415245276355, 85.2918498392157, 23.58204239266529, 32.047140800123714, 15.663771838668453, 62.73409213837115, 80.13968246445067, 62.778846064423895, 57.58033072787298, 8.196694625907504, 76.40782581931734, 46.586682489976994, 114.35381760832637, 194.18203056426375, 34.59402256868242, 21.651933506693812, 23.826349217082235, 15.874130483702587, 14.036216805382255, 10.039572677623797, 54.37395246773005, 51.69284486157993, 25.494806326474627, 22.970498996319368, 29.049559360632642, 62.726248735489285, 5.985104622593393, 74.84996325167715, 70.97074195252529, 30.787962678658396, 14.898851547316111, 23.739679174309803, 131.0865084547394, 24.168222574394733, 14.37073089022526, 96.79395563644708, 27.262246127543477, 25.737726054337667, 20.340371735703464, 31.66474509112033, 46.57698551728102, 17.967807193836734, 59.292920229704265, 136.37043121808114, 17.54353663177534, 62.11276542206669, 21.48315950171234, 31.21179598230219, 5.391109844601892, 31.988454329878167, 71.11033831855057, 40.919044566978684, 25.297123589604777, 39.551484206563536, 80.19014551794595, 25.650910056096734, 5.336689840202448, 18.262915567326992, 118.13576015850285, 70.50607742782249, 29.788165855112435, 9.139218642279456, 5.149952051571178, 9.958332934713537, 95.75046259601393, 25.029793268646767, 214.91696354232363, 48.08420860780603, 40.20608745184056, 27.839534613914925, 58.88639165702904, 11.899245397706084, 30.47116918306545, 5.373817131441746, 24.87775343796074, 124.30687177435296, 128.6698593004463, 66.08445045038765, 64.38259088450616, 61.703001646313155, 51.814026213043114, 30.69536301464755, 36.91263224545558, 42.50890956154804, 35.226959027873065, 11.843565757265504, 43.42630273166789, 66.04710627837264, 30.746867034251878, 41.46239555718222, 52.91180541551975, 25.580055007486077, 133.75708739228938, 34.40303160956148, 20.85543539161491, 60.37378443553194, 30.03128253661255, 30.299903994369124, 16.528548427139636, 5.662333461269984, 17.27228407870122, 10.042434678355434, 13.73836268097965, 6.657999652977235, 32.808162388271285, 15.072944133872037, 8.86649357244533, 57.705624552604874, 18.76057017501909, 17.03910002189733, 43.21193784236018, 15.453375549136315, 71.76141430947406, 31.704235575903716, 8.02462762406903, 19.98030092372793, 180.17910151600915, 24.408890791191936, 5.4888252273976885, 89.3463587978369, 71.64351476969134, 55.37447538437138, 72.1857978352336, 42.679712633416365, 8.643517306005004, 5.561395765676327, 53.40116008827213, 17.19664765390621, 88.2193755643768, 10.918543297228254, 102.41681051040266, 27.202836355716975, 14.669160676969513, 102.80900398772879, 24.623023323470782, 112.11353867917133, 62.56346455250921, 86.35286182958929, 7.19712308406061, 7.154950940050308, 43.42000490054383, 66.37194296632163, 36.12332225708188, 43.400400107590535, 11.562255378660073, 214.58225222048264, 58.99461637531479, 74.54915486019145, 104.29956755666134, 103.65203008084146, 134.79304904453852, 12.899219204141861, 11.128346658731921, 5.118775767946365, 144.87369429286144, 6.0330410946785635, 7.084702728148356, 167.84000102526167, 7.691588857052687, 14.53307104785124, 183.54675656824554, 5.585668255883306, 20.06764064855352, 7.267865703719169, 80.23737809738822, 68.62163808254931, 25.999232394033875, 32.890423499539295, 60.279515643102314, 153.59472938346792, 20.483298019145177, 30.65913104945777, 77.22919222226835, 5.18866916473502, 21.917044705032993, 6.861636473787336, 77.52371763856812, 11.836137897767514, 86.5496144850733, 159.8892914441253, 28.392926426256608, 63.80776368308087, 8.897032944558788, 30.77797499672431, 22.684304195349007, 144.33985614999858, 11.257701021853373, 16.728772810005815, 183.94809856801078, 51.53671035664104, 60.40471128759087, 24.422617324156796, 36.242290016240105, 114.84915416001145, 18.281348353839597, 68.93676993387682, 18.21830886137796, 56.28448348867614, 68.90577910965548, 30.56781386306111, 29.583533394671633, 145.5582293654339, 20.39242065177016, 8.280572763948582, 23.860437961548836, 39.609649648667954, 65.42491419499525, 9.55149631981572, 94.89097058992613, 74.4860911168397, 100.31451595517484, 130.13482868631988, 10.512559695399421, 15.566394667900955, 81.39279492012433, 19.498866769125748, 35.731736579455664, 24.511853364985363, 34.08499860664441, 35.24380606237309, 9.240268095219559, 72.1596171400519, 7.785868806569296, 28.49040438915312, 34.6153542417438, 15.348391179493005, 56.34824251396226, 62.26588691460305, 11.290896566965255, 33.65689035365521, 118.13048166901218, 37.19803721435578, 6.604620932998921, 33.293768826754764, 238.91105150854133, 66.3434514520456, 27.173998390458138, 11.690694191586443, 23.567855531781866, 29.055017785531877, 54.299683535906794, 46.55121776490687, 29.22905125735331, 128.5067144673952, 42.7828020651135, 5.463892327023526, 133.58604478477812, 49.7303689759091, 31.66264001385943, 151.69191151898954, 93.28658861785371, 12.286680209557279, 43.96168087819643, 52.44665134370262, 18.30737912454215, 27.262277964126792, 62.80337129942751, 8.973898478598995, 14.160610039373534, 5.152519243478851, 29.39516946480492, 22.404832869442814, 61.74187752392386, 19.850046754710874, 55.992523970153066, 56.59263062667548, 21.2070585932029, 90.29082554855187, 24.817377674570402, 62.26553339617556, 10.179388207095963, 29.840027943612263, 25.148432215436568, 22.34846039962666, 194.80395689194546, 166.4605500712578, 8.695219865185168, 59.59849268858369, 14.212030126126832, 55.33781689202029, 28.84495210871959, 24.043688649642945, 26.631977377903716, 16.55183059250378, 20.222455248867366, 12.116994118988346, 59.85229158134261, 41.377537346054794, 26.64192663953893, 11.224798581029903, 6.977549470257569, 5.096085765150016, 12.268487781308231, 41.4307998743007, 20.090583717254727, 26.747580149438633, 76.85779785887442, 55.03010645074186, 24.90104027020149, 7.086458102945911, 7.569187318280818, 8.874315565176806, 7.695250335589253, 13.515898608904545, 77.18429990996823, 30.062439751363794, 88.56119199983307, 70.46745129449701, 29.8744497364807, 62.85347353867023, 71.29197731684675, 71.34306514086097, 50.62500220655267, 57.85440943410419, 22.638712254399703, 7.4309962488349965, 9.059045607160723, 45.0859702736866, 15.881033120140875, 84.58786933627641, 24.78097920718625, 58.25011522097172, 43.18242141078524, 34.870120130212115, 50.47119521091089, 68.80437733327997, 5.113183277464542, 14.91190639060965, 17.11870141624891, 77.97525458822486, 70.92655817137924, 18.410011486754357, 14.720262230913116, 7.287047973448669, 10.871291718559236, 57.54149552832244, 6.876003622381181, 7.968412210487223, 23.7507424976705, 16.10940164310589, 13.566601774286635, 38.87194449593684, 49.73875969952127, 35.45645998945266, 79.58453677538628, 9.236829289921983, 65.31867028738485, 52.190274251771115, 13.935945985869546, 42.446639098184924, 16.881420423257893, 82.46633068706194, 12.581882276225468, 10.598166632122277, 77.36410583768759, 13.733641874239257, 62.90548551977255, 11.777443829906712, 44.30551951064687, 34.462657662263986, 61.49912204134962, 28.94204199863049, 17.498114058070623, 28.700534699464015, 62.64861627083927, 10.691477407534023, 23.335499218411513, 87.3184803981755, 15.244022403213467, 9.093179826007699, 97.25350982175077, 67.77102157671955, 11.57517828026321, 19.127400663018648, 14.307975113364787, 5.091175813161763, 12.18189411108122, 8.495822580633002, 6.454646232699643, 39.57714115950007, 20.416618871769916, 30.983031334451482, 8.228442871664601, 5.159450684370454, 63.26234204182821, 27.5285880089019, 12.01176807770946, 53.851871123215744, 62.85649559307462, 12.478560575248734, 15.680807593072624, 54.943221485452305, 8.24350925923982, 189.51925173544402, 20.129206532026206, 7.113113983196451, 68.86700903224099, 5.210895440856173, 12.304246313266205, 16.001925737659978, 89.4118103059515, 6.952981295083042, 46.38527652941901, 24.418136025793643, 20.132192508550933, 39.479414730586335, 19.548667864534494, 20.418216381930918, 15.967772450331399, 45.44186675076523, 9.655117980023904, 15.107998916066029, 97.40949718756524, 7.135729478392734, 5.4658136897423555, 7.676627283868108, 16.921767717040005, 8.423557423532772, 6.137994272533657, 16.29234603103157, 29.20874494194414, 13.868315599188143, 34.26206872342423, 24.546954690578893, 54.280971801214115, 86.15760559422289, 13.455040187688954, 77.83820310843062, 24.137809355059577, 11.047177479278954, 78.6599350271893, 26.44893747934936, 32.90956474793774, 8.014399977262284, 39.97975007787519, 18.201640621549206, 65.13049776372034, 52.273259728624254, 76.00032508552368, 40.44337300506608, 72.11470513534665, 63.45351509360373, 13.155007101433918, 19.07850029166463, 11.629352502912301, 20.236715029391107, 73.47750879727143, 7.318210371663687, 6.470235280062387, 26.417195292758343, 6.680119250188373, 87.8891339336148, 9.891237912467195, 73.16290092732898, 23.500219874358503, 12.67375563074878, 73.1175212212965, 43.07414059986069, 5.743926808181919, 31.324349633784266, 8.661170368286205, 20.094764301948093, 13.713711528343914, 65.22718358105725, 13.502316154216112, 14.18886481410034, 6.709677157330401, 32.76343140317322, 133.99406267323621, 47.43691061492423, 42.817696147762305, 56.923591267712936, 6.962283887841462, 18.35591916934167, 18.69666289808284, 48.48546408647252, 55.345573634861836, 6.568993348418585, 25.533279001565415, 15.747297391074174, 5.390049023981349])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3087773.4375, 3279492.1875, 3291940.625, 3302131.25, 3332648.4375, 3340603.125, 3341659.375, 3341721.875, 3351781.25, 3373592.1875, 3376653.125, 3379860.9375, 3384843.75, 3396104.6875, 3443282.8125, 3445862.5, 3449479.6875, 3450982.8125, 3453625.0, 3462920.3125, 3480990.625, 3481003.125, 3506332.8125, 3512343.75, 3528484.375, 3547621.875, 3550703.125, 3594795.3125, 3598578.125, 3617646.875, 3662160.9375, 3691170.3125, 3697804.6875, 3703453.125, 3707943.75, 3710864.0625, 3711143.75, 3726654.6875, 3728825.0, 3739898.4375, 3759232.8125, 3776081.25, 3785918.75, 3829915.625, 3838668.75, 3854609.375, 3869617.1875, 3903290.625, 3932415.625, 3938745.3125, 3945310.9375, 3949546.875, 3951993.75, 3957442.1875, 3958454.6875, 3958465.625, 3959054.6875, 3960396.875, 3960943.75, 3961643.75, 3962089.0625, 3963312.5, 3963360.9375, 3963529.6875, 3964828.125, 3967365.625, 3970403.125, 3970807.8125, 3972093.75, 3972198.4375, 3972228.125, 3977357.8125, 3987478.125, 3988395.3125, 3988934.375, 3990846.875, 3994178.125, 3996070.3125, 4001526.5625, 4001946.875, 4002076.5625, 4002365.625, 4003300.0, 4004778.125, 4014709.375, 4014735.9375, 4016078.125, 4017732.8125, 4027700.0, 4041681.25, 4046200.0, 4055504.6875, 4058004.6875, 4066754.6875, 4071087.5, 4083629.6875, 4099378.125, 4101448.4375, 4170231.25, 4172356.25, 4178989.0625, 4187404.6875, 4233329.6875, 4237554.6875, 4257335.9375, 4258609.375, 4260948.4375, 4286045.3125, 4286953.125, 4287385.9375, 4287959.375, 4293464.0625, 4293470.3125, 4298512.5, 4299073.4375, 4299826.5625, 4302045.3125, 4302190.625, 4303270.3125, 4306889.0625, 4313354.6875, 4313610.9375, 4316346.875, 4365762.5, 4367503.125, 4372154.6875, 4372220.3125, 4372660.9375, 4372743.75, 4376043.75, 4382392.1875, 4382400.0, 4402967.1875, 4408334.375, 4411962.5, 4437851.5625, 4439017.1875, 4440671.875, 4448135.9375, 4460140.625, 4462029.6875, 4468298.4375, 4475170.3125, 4476517.1875, 4478134.375, 4487668.75, 4493178.125, 4504848.4375, 4510612.5, 4514323.4375, 4519190.625, 4523985.9375, 4524639.0625, 4535112.5, 4536723.4375, 4551214.0625, 4587389.0625, 4592243.75, 4603082.8125, 4604612.5, 4611210.9375, 4611635.9375, 4611650.0, 4612910.9375, 4614756.25, 4621246.875, 4626018.75, 4648462.5, 4661704.6875, 4670635.9375, 4677840.625, 4685984.375, 4689985.9375, 4703406.25, 4706595.3125, 4726128.125, 4727179.6875, 4729268.75, 4731562.5, 4741643.75, 4741659.375, 4746734.375, 4751593.75, 4754264.0625, 4761575.0, 4767590.625, 4770246.875, 4772746.875, 4777084.375, 4777148.4375, 4781417.1875, 4784057.8125, 4784570.3125, 4786085.9375, 4788357.8125, 4791876.5625, 4792198.4375, 4799470.3125, 4806020.3125, 4818323.4375, 4822328.125, 4823754.6875, 4825445.3125, 4826709.375, 4836514.0625, 4852610.9375, 4861657.8125, 4864784.375, 4870968.75, 4881329.6875, 4881331.25, 4882671.875, 4887182.8125, 4887264.0625, 4887571.875, 4942845.3125, 5202618.75, 5267948.4375, 5304584.375, 5311142.1875, 5431157.8125, 5439004.6875, 5566810.9375, 5613315.625, 5617290.625, 5661426.5625, 5729892.1875, 5777900.0, 5788240.625, 5855723.4375, 5865017.1875, 5885287.5, 5913842.1875, 5945917.1875, 5957729.6875, 5959490.625, 5963306.25, 5968056.25, 5997043.75, 6033198.4375, 6166170.3125, 6226726.5625, 6230910.9375, 6272492.1875, 6276678.125, 6350546.875, 6384792.1875, 6414335.9375, 6415671.875, 6415682.8125, 6444023.4375, 6451354.6875, 6453412.5, 6455237.5, 6455496.875, 6455756.25, 6455810.9375, 6455815.625, 6456259.375, 6456390.625, 6457600.0, 6457610.9375, 6457835.9375, 6458114.0625, 6458634.375, 6458912.5, 6459537.5, 6460218.75, 6460381.25, 6464935.9375, 6482834.375, 6501218.75, 6567126.5625, 6575485.9375, 6575707.8125, 6577029.6875, 6577520.3125, 6580496.875, 6580612.5, 6582296.875, 6595951.5625, 6622120.3125, 6626510.9375, 6626848.4375, 6628718.75, 6628795.3125, 6630976.5625, 6632500.0, 6634203.125, 6635343.75, 6635479.6875, 6636526.5625, 6637414.0625, 6694087.5, 6694162.5, 6727965.625, 6727979.6875, 6728510.9375, 6768018.75, 6769990.625, 6819342.1875, 6824698.4375, 6827067.1875, 6828110.9375, 6920218.75, 6920220.3125, 6924500.0, 6955023.4375, 6955450.0, 6957070.3125, 6984871.875, 6999907.8125, 7009839.0625, 7115892.1875, 7120931.25, 7171362.5, 7177581.25, 7212764.0625, 7233709.375, 7291923.4375, 7337326.5625, 7337984.375, 7358348.4375, 7370460.9375, 7405179.6875, 7479779.6875, 7542878.125, 7546215.625, 7546226.5625, 7579470.3125, 7590853.125, 7688746.875, 7738551.5625, 7801984.375, 7835351.5625, 7835420.3125, 7908001.5625, 7959790.625, 7999307.8125, 8046478.125, 8080765.625, 8094612.5, 8099781.25, 8123092.1875, 8137417.1875, 8251940.625, 8302825.0, 8340545.3125, 8367860.9375, 8370220.3125, 8449735.9375, 8458984.375, 8459814.0625, 8467187.5, 8469817.1875, 8488120.3125, 8493984.375, 8497235.9375, 8498331.25, 8501132.8125, 8504198.4375, 8579548.4375, 8644642.1875, 8657851.5625, 8722443.75, 8791884.375, 8810348.4375, 8823653.125, 8966678.125, 9049198.4375, 9186187.5, 9296284.375, 9325573.4375, 9330907.8125, 9357398.4375, 9360792.1875, 9361407.8125, 9447243.75, 9448639.0625, 9449212.5, 9449701.5625, 9462470.3125, 9489589.0625, 9491907.8125, 9540006.25, 9636868.75, 9727014.0625, 9732673.4375, 9732914.0625, 9735464.0625, 9742296.875, 9744396.875, 9745239.0625, 9803107.8125, 9900350.0, 9908162.5, 9910843.75, 9948620.3125, 9949189.0625, 9951018.75, 9951565.625, 9954076.5625, 9954132.8125, 9975871.875, 9977475.0, 9977601.5625, 9978351.5625, 9978410.9375, 9978426.5625, 9978746.875, 9978934.375, 9979885.9375, 9980848.4375, 9982334.375, 9983507.8125, 9984860.9375, 9984939.0625, 9984970.3125, 9991339.0625, 10018678.125, 10020339.0625, 10020495.3125, 10020675.0, 10021337.5, 10021643.75, 10021660.9375, 10022270.3125, 10022396.875, 10022862.5, 10023589.0625, 10024279.6875, 10024595.3125, 10026085.9375, 10028837.5, 10095937.5, 10107328.125, 10108204.6875, 10127207.8125, 10133090.625, 10133110.9375, 10149882.8125, 10155068.75, 10155465.625, 10159973.4375, 10160479.6875, 10160507.8125, 10175629.6875, 10187621.875, 10187737.5, 10205395.3125, 10205668.75, 10268126.5625, 10268142.1875, 11172354.6875, 11177729.6875, 11177839.0625, 11178792.1875, 11179850.0, 11180498.4375, 11180504.6875, 11183793.75, 11186384.375, 11187345.3125, 11190289.0625, 11190915.625, 11193473.4375, 11195215.625, 11202617.1875, 11206462.5, 11218329.6875, 11295078.125, 11317231.25, 11317300.0, 11318243.75, 11318614.0625, 11319217.1875, 11319470.3125, 11319945.3125, 11320190.625, 11321784.375, 11321896.875, 11321900.0, 11321907.8125, 11322157.8125, 11324617.1875, 11325123.4375, 11325620.3125, 11325760.9375, 11325940.625, 11326442.1875, 11326460.9375, 11326703.125, 11326810.9375, 11327181.25, 11327264.0625, 11327660.9375, 11327950.0, 11328082.8125, 11328645.3125, 11328690.625, 11329137.5, 11329428.125, 11329448.4375, 11329817.1875, 11329920.3125, 11330070.3125, 11331375.0, 11331465.625, 11331715.625, 11331851.5625, 11331929.6875, 11331940.625, 11332281.25, 11332867.1875, 11333100.0, 11333157.8125, 11334282.8125, 11334704.6875, 11334878.125, 11335160.9375, 11335171.875, 11335220.3125, 11335862.5, 11336420.3125, 11336692.1875, 11337010.9375, 11337740.625, 11337904.6875, 11337912.5, 11338114.0625, 11338264.0625, 11338621.875, 11338668.75, 11339717.1875, 11339978.125, 11340881.25, 11341218.75, 11341492.1875, 11341568.75, 11342414.0625, 11342489.0625, 11342684.375, 11342929.6875, 11342934.375, 11343203.125, 11344056.25, 11344065.625, 11344345.3125, 11346003.125, 11346782.8125, 11347150.0, 11347812.5, 11347867.1875, 11347870.3125, 11347881.25, 11348556.25, 11349464.0625, 11349545.3125, 11349548.4375, 11349603.125, 11350007.8125, 11350568.75, 11350673.4375, 11351117.1875, 11351239.0625, 11351496.875, 11352053.125, 11352648.4375, 11352682.8125, 11353907.8125, 11354253.125, 11355803.125, 11355898.4375, 11356171.875, 11356864.0625, 11357542.1875, 11358006.25, 11358673.4375, 11359612.5, 11359800.0, 11360937.5, 11361762.5, 11362464.0625, 11362490.625, 11362646.875, 11363053.125, 11363376.5625, 11363453.125, 11365003.125, 11365451.5625, 11365660.9375, 11365776.5625, 11366259.375, 11367289.0625, 11367437.5, 11367587.5, 11367898.4375, 11368028.125, 11368704.6875, 11369042.1875, 11369257.8125, 11369528.125, 11369960.9375, 11370104.6875, 11370121.875, 11370282.8125, 11370321.875, 11370715.625, 11370957.8125, 11371170.3125, 11371384.375, 11371710.9375, 11372156.25, 11372167.1875, 11372243.75, 11372289.0625, 11372296.875, 11372796.875, 11374464.0625, 11375037.5, 11375103.125, 11375614.0625, 11375701.5625, 11375831.25, 11376528.125, 11376701.5625, 11376862.5, 11377450.0, 11377623.4375, 11377670.3125, 11379512.5, 11380315.625, 11381628.125, 11382507.8125, 11390015.625, 11393490.625, 11400914.0625, 11412100.0], [9.021341576318875, 10.477288474122311, 57.54232546967111, 43.36223247797474, 11.30619164980767, 41.089963124880875, 8.166058461384358, 5.240813951353239, 42.11975825407877, 5.648981120762856, 38.44809445452768, 34.398871321035394, 33.515944875368504, 9.225651296322075, 5.637175833110015, 12.453412421269853, 9.95662535622695, 42.85192698978631, 58.44927274427586, 71.2157030505174, 67.42362637935705, 7.135961360725214, 79.46549018434845, 9.849819364019993, 25.044360851243567, 71.60105463631187, 33.04248052611274, 65.5889921758571, 60.85461379432205, 13.304759103373232, 9.179640706374053, 60.64724193305344, 19.44004131969818, 5.636328907076851, 24.848293754721837, 5.290449908872773, 5.171263084338775, 7.995774964460886, 75.76501538774251, 24.19017595839289, 60.14228394731865, 29.7064543888022, 7.500352011449034, 10.750330995411119, 13.826946834335155, 33.447354321558436, 11.94859398986753, 75.03945225597064, 8.635126547536887, 52.27764754239169, 67.9266378567987, 46.04808813836873, 101.13580485579475, 123.18776831456105, 22.60207153664546, 60.556266716845286, 15.298635025912702, 85.21712537952381, 28.611818329301595, 14.760795640073713, 19.371264829730862, 51.332686155110764, 12.868475116809758, 17.07158809025163, 6.608998110940824, 81.28311107751189, 7.830052982107149, 64.65233504921224, 16.569579362659873, 16.702030866968922, 9.74260796417946, 30.01294975813429, 92.53561265167383, 19.55944798532267, 24.604704762767536, 165.6791276678834, 12.085598395197227, 11.811065604939285, 66.74262162672616, 33.35766781848508, 10.343889980390738, 51.53812635793888, 5.185852330044586, 33.55835004815951, 57.065205701517044, 21.973078643081088, 31.520777244848492, 8.564235694264939, 7.92196593358585, 9.706312248906553, 49.55153439496237, 32.43307708385775, 27.289044132673123, 32.09290456982849, 28.004534571688644, 13.649655348425629, 139.06591816922727, 7.937612684660185, 18.715945570827422, 135.6264237347823, 13.451375773570465, 65.99526794183818, 93.94145740681168, 132.62053297411637, 73.88730511353756, 9.555031159786623, 10.159484030408336, 28.753230155683998, 56.06243456487948, 6.88037298906764, 9.746176409734632, 76.48351885828143, 11.677681085553473, 7.8734026627145, 83.13458928723722, 12.47601459755017, 56.369336043489085, 38.79422547130149, 9.747888887236929, 35.18213256046644, 17.311572216100345, 6.654798398315644, 21.205006035962363, 146.73375709112386, 29.41553357409855, 14.77133786019102, 102.74532226215821, 10.428520389410174, 5.174315570063226, 23.882563835075153, 20.85725712502494, 18.230827555127817, 37.011115066793565, 24.077082261451345, 14.445864546598049, 121.86511622824248, 5.383008799121236, 10.342308830756494, 38.017352474591966, 5.153762521561403, 36.28312707571996, 104.38080187740073, 19.540276420624945, 8.280307167487964, 26.959179897435153, 57.23915855070114, 81.53097542794902, 16.041392747083027, 31.332282256928348, 17.116063325385987, 82.89291172723831, 50.942663206557036, 61.16509563160895, 13.241508996154332, 92.20384760567592, 84.32115200732882, 7.166571345499603, 32.25415245276355, 85.2918498392157, 23.58204239266529, 32.047140800123714, 15.663771838668453, 62.73409213837115, 80.13968246445067, 62.778846064423895, 57.58033072787298, 8.196694625907504, 76.40782581931734, 46.586682489976994, 114.35381760832637, 194.18203056426375, 34.59402256868242, 21.651933506693812, 23.826349217082235, 15.874130483702587, 14.036216805382255, 10.039572677623797, 54.37395246773005, 51.69284486157993, 25.494806326474627, 22.970498996319368, 29.049559360632642, 62.726248735489285, 5.985104622593393, 74.84996325167715, 70.97074195252529, 30.787962678658396, 14.898851547316111, 23.739679174309803, 131.0865084547394, 24.168222574394733, 14.37073089022526, 96.79395563644708, 27.262246127543477, 25.737726054337667, 20.340371735703464, 31.66474509112033, 46.57698551728102, 17.967807193836734, 59.292920229704265, 136.37043121808114, 17.54353663177534, 62.11276542206669, 21.48315950171234, 31.21179598230219, 5.391109844601892, 31.988454329878167, 71.11033831855057, 40.919044566978684, 25.297123589604777, 39.551484206563536, 80.19014551794595, 25.650910056096734, 5.336689840202448, 18.262915567326992, 118.13576015850285, 70.50607742782249, 29.788165855112435, 9.139218642279456, 5.149952051571178, 9.958332934713537, 95.75046259601393, 25.029793268646767, 214.91696354232363, 48.08420860780603, 40.20608745184056, 27.839534613914925, 58.88639165702904, 11.899245397706084, 30.47116918306545, 5.373817131441746, 24.87775343796074, 124.30687177435296, 128.6698593004463, 66.08445045038765, 64.38259088450616, 61.703001646313155, 51.814026213043114, 30.69536301464755, 36.91263224545558, 42.50890956154804, 35.226959027873065, 11.843565757265504, 43.42630273166789, 66.04710627837264, 30.746867034251878, 41.46239555718222, 52.91180541551975, 25.580055007486077, 133.75708739228938, 34.40303160956148, 20.85543539161491, 60.37378443553194, 30.03128253661255, 30.299903994369124, 16.528548427139636, 5.662333461269984, 17.27228407870122, 10.042434678355434, 13.73836268097965, 6.657999652977235, 32.808162388271285, 15.072944133872037, 8.86649357244533, 57.705624552604874, 18.76057017501909, 17.03910002189733, 43.21193784236018, 15.453375549136315, 71.76141430947406, 31.704235575903716, 8.02462762406903, 19.98030092372793, 180.17910151600915, 24.408890791191936, 5.4888252273976885, 89.3463587978369, 71.64351476969134, 55.37447538437138, 72.1857978352336, 42.679712633416365, 8.643517306005004, 5.561395765676327, 53.40116008827213, 17.19664765390621, 88.2193755643768, 10.918543297228254, 102.41681051040266, 27.202836355716975, 14.669160676969513, 102.80900398772879, 24.623023323470782, 112.11353867917133, 62.56346455250921, 86.35286182958929, 7.19712308406061, 7.154950940050308, 43.42000490054383, 66.37194296632163, 36.12332225708188, 43.400400107590535, 11.562255378660073, 214.58225222048264, 58.99461637531479, 74.54915486019145, 104.29956755666134, 103.65203008084146, 134.79304904453852, 12.899219204141861, 11.128346658731921, 5.118775767946365, 144.87369429286144, 6.0330410946785635, 7.084702728148356, 167.84000102526167, 7.691588857052687, 14.53307104785124, 183.54675656824554, 5.585668255883306, 20.06764064855352, 7.267865703719169, 80.23737809738822, 68.62163808254931, 25.999232394033875, 32.890423499539295, 60.279515643102314, 153.59472938346792, 20.483298019145177, 30.65913104945777, 77.22919222226835, 5.18866916473502, 21.917044705032993, 6.861636473787336, 77.52371763856812, 11.836137897767514, 86.5496144850733, 159.8892914441253, 28.392926426256608, 63.80776368308087, 8.897032944558788, 30.77797499672431, 22.684304195349007, 144.33985614999858, 11.257701021853373, 16.728772810005815, 183.94809856801078, 51.53671035664104, 60.40471128759087, 24.422617324156796, 36.242290016240105, 114.84915416001145, 18.281348353839597, 68.93676993387682, 18.21830886137796, 56.28448348867614, 68.90577910965548, 30.56781386306111, 29.583533394671633, 145.5582293654339, 20.39242065177016, 8.280572763948582, 23.860437961548836, 39.609649648667954, 65.42491419499525, 9.55149631981572, 94.89097058992613, 74.4860911168397, 100.31451595517484, 130.13482868631988, 10.512559695399421, 15.566394667900955, 81.39279492012433, 19.498866769125748, 35.731736579455664, 24.511853364985363, 34.08499860664441, 35.24380606237309, 9.240268095219559, 72.1596171400519, 7.785868806569296, 28.49040438915312, 34.6153542417438, 15.348391179493005, 56.34824251396226, 62.26588691460305, 11.290896566965255, 33.65689035365521, 118.13048166901218, 37.19803721435578, 6.604620932998921, 33.293768826754764, 238.91105150854133, 66.3434514520456, 27.173998390458138, 11.690694191586443, 23.567855531781866, 29.055017785531877, 54.299683535906794, 46.55121776490687, 29.22905125735331, 128.5067144673952, 42.7828020651135, 5.463892327023526, 133.58604478477812, 49.7303689759091, 31.66264001385943, 151.69191151898954, 93.28658861785371, 12.286680209557279, 43.96168087819643, 52.44665134370262, 18.30737912454215, 27.262277964126792, 62.80337129942751, 8.973898478598995, 14.160610039373534, 5.152519243478851, 29.39516946480492, 22.404832869442814, 61.74187752392386, 19.850046754710874, 55.992523970153066, 56.59263062667548, 21.2070585932029, 90.29082554855187, 24.817377674570402, 62.26553339617556, 10.179388207095963, 29.840027943612263, 25.148432215436568, 22.34846039962666, 194.80395689194546, 166.4605500712578, 8.695219865185168, 59.59849268858369, 14.212030126126832, 55.33781689202029, 28.84495210871959, 24.043688649642945, 26.631977377903716, 16.55183059250378, 20.222455248867366, 12.116994118988346, 59.85229158134261, 41.377537346054794, 26.64192663953893, 11.224798581029903, 6.977549470257569, 5.096085765150016, 12.268487781308231, 41.4307998743007, 20.090583717254727, 26.747580149438633, 76.85779785887442, 55.03010645074186, 24.90104027020149, 7.086458102945911, 7.569187318280818, 8.874315565176806, 7.695250335589253, 13.515898608904545, 77.18429990996823, 30.062439751363794, 88.56119199983307, 70.46745129449701, 29.8744497364807, 62.85347353867023, 71.29197731684675, 71.34306514086097, 50.62500220655267, 57.85440943410419, 22.638712254399703, 7.4309962488349965, 9.059045607160723, 45.0859702736866, 15.881033120140875, 84.58786933627641, 24.78097920718625, 58.25011522097172, 43.18242141078524, 34.870120130212115, 50.47119521091089, 68.80437733327997, 5.113183277464542, 14.91190639060965, 17.11870141624891, 77.97525458822486, 70.92655817137924, 18.410011486754357, 14.720262230913116, 7.287047973448669, 10.871291718559236, 57.54149552832244, 6.876003622381181, 7.968412210487223, 23.7507424976705, 16.10940164310589, 13.566601774286635, 38.87194449593684, 49.73875969952127, 35.45645998945266, 79.58453677538628, 9.236829289921983, 65.31867028738485, 52.190274251771115, 13.935945985869546, 42.446639098184924, 16.881420423257893, 82.46633068706194, 12.581882276225468, 10.598166632122277, 77.36410583768759, 13.733641874239257, 62.90548551977255, 11.777443829906712, 44.30551951064687, 34.462657662263986, 61.49912204134962, 28.94204199863049, 17.498114058070623, 28.700534699464015, 62.64861627083927, 10.691477407534023, 23.335499218411513, 87.3184803981755, 15.244022403213467, 9.093179826007699, 97.25350982175077, 67.77102157671955, 11.57517828026321, 19.127400663018648, 14.307975113364787, 5.091175813161763, 12.18189411108122, 8.495822580633002, 6.454646232699643, 39.57714115950007, 20.416618871769916, 30.983031334451482, 8.228442871664601, 5.159450684370454, 63.26234204182821, 27.5285880089019, 12.01176807770946, 53.851871123215744, 62.85649559307462, 12.478560575248734, 15.680807593072624, 54.943221485452305, 8.24350925923982, 189.51925173544402, 20.129206532026206, 7.113113983196451, 68.86700903224099, 5.210895440856173, 12.304246313266205, 16.001925737659978, 89.4118103059515, 6.952981295083042, 46.38527652941901, 24.418136025793643, 20.132192508550933, 39.479414730586335, 19.548667864534494, 20.418216381930918, 15.967772450331399, 45.44186675076523, 9.655117980023904, 15.107998916066029, 97.40949718756524, 7.135729478392734, 5.4658136897423555, 7.676627283868108, 16.921767717040005, 8.423557423532772, 6.137994272533657, 16.29234603103157, 29.20874494194414, 13.868315599188143, 34.26206872342423, 24.546954690578893, 54.280971801214115, 86.15760559422289, 13.455040187688954, 77.83820310843062, 24.137809355059577, 11.047177479278954, 78.6599350271893, 26.44893747934936, 32.90956474793774, 8.014399977262284, 39.97975007787519, 18.201640621549206, 65.13049776372034, 52.273259728624254, 76.00032508552368, 40.44337300506608, 72.11470513534665, 63.45351509360373, 13.155007101433918, 19.07850029166463, 11.629352502912301, 20.236715029391107, 73.47750879727143, 7.318210371663687, 6.470235280062387, 26.417195292758343, 6.680119250188373, 87.8891339336148, 9.891237912467195, 73.16290092732898, 23.500219874358503, 12.67375563074878, 73.1175212212965, 43.07414059986069, 5.743926808181919, 31.324349633784266, 8.661170368286205, 20.094764301948093, 13.713711528343914, 65.22718358105725, 13.502316154216112, 14.18886481410034, 6.709677157330401, 32.76343140317322, 133.99406267323621, 47.43691061492423, 42.817696147762305, 56.923591267712936, 6.962283887841462, 18.35591916934167, 18.69666289808284, 48.48546408647252, 55.345573634861836, 6.568993348418585, 25.533279001565415, 15.747297391074174, 5.390049023981349])
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);
([3087773.4375, 3279492.1875, 3291940.625, 3302131.25, 3332648.4375, 3340603.125, 3341659.375, 3341721.875, 3351781.25, 3373592.1875, 3376653.125, 3379860.9375, 3384843.75, 3396104.6875, 3443282.8125, 3445862.5, 3449479.6875, 3450982.8125, 3453625.0, 3462920.3125, 3480990.625, 3481003.125, 3506332.8125, 3512343.75, 3528484.375, 3547621.875, 3550703.125, 3594795.3125, 3598578.125, 3617646.875, 3662160.9375, 3691170.3125, 3697804.6875, 3703453.125, 3707943.75, 3710864.0625, 3711143.75, 3726654.6875, 3728825.0, 3739898.4375, 3759232.8125, 3776081.25, 3785918.75, 3829915.625, 3838668.75, 3854609.375, 3869617.1875, 3903290.625, 3932415.625, 3938745.3125, 3945310.9375, 3949546.875, 3951993.75, 3957442.1875, 3958454.6875, 3958465.625, 3959054.6875, 3960396.875, 3960943.75, 3961643.75, 3962089.0625, 3963312.5, 3963360.9375, 3963529.6875, 3964828.125, 3967365.625, 3970403.125, 3970807.8125, 3972093.75, 3972198.4375, 3972228.125, 3977357.8125, 3987478.125, 3988395.3125, 3988934.375, 3990846.875, 3994178.125, 3996070.3125, 4001526.5625, 4001946.875, 4002076.5625, 4002365.625, 4003300.0, 4004778.125, 4014709.375, 4014735.9375, 4016078.125, 4017732.8125, 4027700.0, 4041681.25, 4046200.0, 4055504.6875, 4058004.6875, 4066754.6875, 4071087.5, 4083629.6875, 4099378.125, 4101448.4375, 4170231.25, 4172356.25, 4178989.0625, 4187404.6875, 4233329.6875, 4237554.6875, 4257335.9375, 4258609.375, 4260948.4375, 4286045.3125, 4286953.125, 4287385.9375, 4287959.375, 4293464.0625, 4293470.3125, 4298512.5, 4299073.4375, 4299826.5625, 4302045.3125, 4302190.625, 4303270.3125, 4306889.0625, 4313354.6875, 4313610.9375, 4316346.875, 4365762.5, 4367503.125, 4372154.6875, 4372220.3125, 4372660.9375, 4372743.75, 4376043.75, 4382392.1875, 4382400.0, 4402967.1875, 4408334.375, 4411962.5, 4437851.5625, 4439017.1875, 4440671.875, 4448135.9375, 4460140.625, 4462029.6875, 4468298.4375, 4475170.3125, 4476517.1875, 4478134.375, 4487668.75, 4493178.125, 4504848.4375, 4510612.5, 4514323.4375, 4519190.625, 4523985.9375, 4524639.0625, 4535112.5, 4536723.4375, 4551214.0625, 4587389.0625, 4592243.75, 4603082.8125, 4604612.5, 4611210.9375, 4611635.9375, 4611650.0, 4612910.9375, 4614756.25, 4621246.875, 4626018.75, 4648462.5, 4661704.6875, 4670635.9375, 4677840.625, 4685984.375, 4689985.9375, 4703406.25, 4706595.3125, 4726128.125, 4727179.6875, 4729268.75, 4731562.5, 4741643.75, 4741659.375, 4746734.375, 4751593.75, 4754264.0625, 4761575.0, 4767590.625, 4770246.875, 4772746.875, 4777084.375, 4777148.4375, 4781417.1875, 4784057.8125, 4784570.3125, 4786085.9375, 4788357.8125, 4791876.5625, 4792198.4375, 4799470.3125, 4806020.3125, 4818323.4375, 4822328.125, 4823754.6875, 4825445.3125, 4826709.375, 4836514.0625, 4852610.9375, 4861657.8125, 4864784.375, 4870968.75, 4881329.6875, 4881331.25, 4882671.875, 4887182.8125, 4887264.0625, 4887571.875, 4942845.3125, 5202618.75, 5267948.4375, 5304584.375, 5311142.1875, 5431157.8125, 5439004.6875, 5566810.9375, 5613315.625, 5617290.625, 5661426.5625, 5729892.1875, 5777900.0, 5788240.625, 5855723.4375, 5865017.1875, 5885287.5, 5913842.1875, 5945917.1875, 5957729.6875, 5959490.625, 5963306.25, 5968056.25, 5997043.75, 6033198.4375, 6166170.3125, 6226726.5625, 6230910.9375, 6272492.1875, 6276678.125, 6350546.875, 6384792.1875, 6414335.9375, 6415671.875, 6415682.8125, 6444023.4375, 6451354.6875, 6453412.5, 6455237.5, 6455496.875, 6455756.25, 6455810.9375, 6455815.625, 6456259.375, 6456390.625, 6457600.0, 6457610.9375, 6457835.9375, 6458114.0625, 6458634.375, 6458912.5, 6459537.5, 6460218.75, 6460381.25, 6464935.9375, 6482834.375, 6501218.75, 6567126.5625, 6575485.9375, 6575707.8125, 6577029.6875, 6577520.3125, 6580496.875, 6580612.5, 6582296.875, 6595951.5625, 6622120.3125, 6626510.9375, 6626848.4375, 6628718.75, 6628795.3125, 6630976.5625, 6632500.0, 6634203.125, 6635343.75, 6635479.6875, 6636526.5625, 6637414.0625, 6694087.5, 6694162.5, 6727965.625, 6727979.6875, 6728510.9375, 6768018.75, 6769990.625, 6819342.1875, 6824698.4375, 6827067.1875, 6828110.9375, 6920218.75, 6920220.3125, 6924500.0, 6955023.4375, 6955450.0, 6957070.3125, 6984871.875, 6999907.8125, 7009839.0625, 7115892.1875, 7120931.25, 7171362.5, 7177581.25, 7212764.0625, 7233709.375, 7291923.4375, 7337326.5625, 7337984.375, 7358348.4375, 7370460.9375, 7405179.6875, 7479779.6875, 7542878.125, 7546215.625, 7546226.5625, 7579470.3125, 7590853.125, 7688746.875, 7738551.5625, 7801984.375, 7835351.5625, 7835420.3125, 7908001.5625, 7959790.625, 7999307.8125, 8046478.125, 8080765.625, 8094612.5, 8099781.25, 8123092.1875, 8137417.1875, 8251940.625, 8302825.0, 8340545.3125, 8367860.9375, 8370220.3125, 8449735.9375, 8458984.375, 8459814.0625, 8467187.5, 8469817.1875, 8488120.3125, 8493984.375, 8497235.9375, 8498331.25, 8501132.8125, 8504198.4375, 8579548.4375, 8644642.1875, 8657851.5625, 8722443.75, 8791884.375, 8810348.4375, 8823653.125, 8966678.125, 9049198.4375, 9186187.5, 9296284.375, 9325573.4375, 9330907.8125, 9357398.4375, 9360792.1875, 9361407.8125, 9447243.75, 9448639.0625, 9449212.5, 9449701.5625, 9462470.3125, 9489589.0625, 9491907.8125, 9540006.25, 9636868.75, 9727014.0625, 9732673.4375, 9732914.0625, 9735464.0625, 9742296.875, 9744396.875, 9745239.0625, 9803107.8125, 9900350.0, 9908162.5, 9910843.75, 9948620.3125, 9949189.0625, 9951018.75, 9951565.625, 9954076.5625, 9954132.8125, 9975871.875, 9977475.0, 9977601.5625, 9978351.5625, 9978410.9375, 9978426.5625, 9978746.875, 9978934.375, 9979885.9375, 9980848.4375, 9982334.375, 9983507.8125, 9984860.9375, 9984939.0625, 9984970.3125, 9991339.0625, 10018678.125, 10020339.0625, 10020495.3125, 10020675.0, 10021337.5, 10021643.75, 10021660.9375, 10022270.3125, 10022396.875, 10022862.5, 10023589.0625, 10024279.6875, 10024595.3125, 10026085.9375, 10028837.5, 10095937.5, 10107328.125, 10108204.6875, 10127207.8125, 10133090.625, 10133110.9375, 10149882.8125, 10155068.75, 10155465.625, 10159973.4375, 10160479.6875, 10160507.8125, 10175629.6875, 10187621.875, 10187737.5, 10205395.3125, 10205668.75, 10268126.5625, 10268142.1875, 11172354.6875, 11177729.6875, 11177839.0625, 11178792.1875, 11179850.0, 11180498.4375, 11180504.6875, 11183793.75, 11186384.375, 11187345.3125, 11190289.0625, 11190915.625, 11193473.4375, 11195215.625, 11202617.1875, 11206462.5, 11218329.6875, 11295078.125, 11317231.25, 11317300.0, 11318243.75, 11318614.0625, 11319217.1875, 11319470.3125, 11319945.3125, 11320190.625, 11321784.375, 11321896.875, 11321900.0, 11321907.8125, 11322157.8125, 11324617.1875, 11325123.4375, 11325620.3125, 11325760.9375, 11325940.625, 11326442.1875, 11326460.9375, 11326703.125, 11326810.9375, 11327181.25, 11327264.0625, 11327660.9375, 11327950.0, 11328082.8125, 11328645.3125, 11328690.625, 11329137.5, 11329428.125, 11329448.4375, 11329817.1875, 11329920.3125, 11330070.3125, 11331375.0, 11331465.625, 11331715.625, 11331851.5625, 11331929.6875, 11331940.625, 11332281.25, 11332867.1875, 11333100.0, 11333157.8125, 11334282.8125, 11334704.6875, 11334878.125, 11335160.9375, 11335171.875, 11335220.3125, 11335862.5, 11336420.3125, 11336692.1875, 11337010.9375, 11337740.625, 11337904.6875, 11337912.5, 11338114.0625, 11338264.0625, 11338621.875, 11338668.75, 11339717.1875, 11339978.125, 11340881.25, 11341218.75, 11341492.1875, 11341568.75, 11342414.0625, 11342489.0625, 11342684.375, 11342929.6875, 11342934.375, 11343203.125, 11344056.25, 11344065.625, 11344345.3125, 11346003.125, 11346782.8125, 11347150.0, 11347812.5, 11347867.1875, 11347870.3125, 11347881.25, 11348556.25, 11349464.0625, 11349545.3125, 11349548.4375, 11349603.125, 11350007.8125, 11350568.75, 11350673.4375, 11351117.1875, 11351239.0625, 11351496.875, 11352053.125, 11352648.4375, 11352682.8125, 11353907.8125, 11354253.125, 11355803.125, 11355898.4375, 11356171.875, 11356864.0625, 11357542.1875, 11358006.25, 11358673.4375, 11359612.5, 11359800.0, 11360937.5, 11361762.5, 11362464.0625, 11362490.625, 11362646.875, 11363053.125, 11363376.5625, 11363453.125, 11365003.125, 11365451.5625, 11365660.9375, 11365776.5625, 11366259.375, 11367289.0625, 11367437.5, 11367587.5, 11367898.4375, 11368028.125, 11368704.6875, 11369042.1875, 11369257.8125, 11369528.125, 11369960.9375, 11370104.6875, 11370121.875, 11370282.8125, 11370321.875, 11370715.625, 11370957.8125, 11371170.3125, 11371384.375, 11371710.9375, 11372156.25, 11372167.1875, 11372243.75, 11372289.0625, 11372296.875, 11372796.875, 11374464.0625, 11375037.5, 11375103.125, 11375614.0625, 11375701.5625, 11375831.25, 11376528.125, 11376701.5625, 11376862.5, 11377450.0, 11377623.4375, 11377670.3125, 11379512.5, 11380315.625, 11381628.125, 11382507.8125, 11390015.625, 11393490.625, 11400914.0625, 11412100.0], [9.021341576318875, 10.477288474122311, 57.54232546967111, 43.36223247797474, 11.30619164980767, 41.089963124880875, 8.166058461384358, 5.240813951353239, 42.11975825407877, 5.648981120762856, 38.44809445452768, 34.398871321035394, 33.515944875368504, 9.225651296322075, 5.637175833110015, 12.453412421269853, 9.95662535622695, 42.85192698978631, 58.44927274427586, 71.2157030505174, 67.42362637935705, 7.135961360725214, 79.46549018434845, 9.849819364019993, 25.044360851243567, 71.60105463631187, 33.04248052611274, 65.5889921758571, 60.85461379432205, 13.304759103373232, 9.179640706374053, 60.64724193305344, 19.44004131969818, 5.636328907076851, 24.848293754721837, 5.290449908872773, 5.171263084338775, 7.995774964460886, 75.76501538774251, 24.19017595839289, 60.14228394731865, 29.7064543888022, 7.500352011449034, 10.750330995411119, 13.826946834335155, 33.447354321558436, 11.94859398986753, 75.03945225597064, 8.635126547536887, 52.27764754239169, 67.9266378567987, 46.04808813836873, 101.13580485579475, 123.18776831456105, 22.60207153664546, 60.556266716845286, 15.298635025912702, 85.21712537952381, 28.611818329301595, 14.760795640073713, 19.371264829730862, 51.332686155110764, 12.868475116809758, 17.07158809025163, 6.608998110940824, 81.28311107751189, 7.830052982107149, 64.65233504921224, 16.569579362659873, 16.702030866968922, 9.74260796417946, 30.01294975813429, 92.53561265167383, 19.55944798532267, 24.604704762767536, 165.6791276678834, 12.085598395197227, 11.811065604939285, 66.74262162672616, 33.35766781848508, 10.343889980390738, 51.53812635793888, 5.185852330044586, 33.55835004815951, 57.065205701517044, 21.973078643081088, 31.520777244848492, 8.564235694264939, 7.92196593358585, 9.706312248906553, 49.55153439496237, 32.43307708385775, 27.289044132673123, 32.09290456982849, 28.004534571688644, 13.649655348425629, 139.06591816922727, 7.937612684660185, 18.715945570827422, 135.6264237347823, 13.451375773570465, 65.99526794183818, 93.94145740681168, 132.62053297411637, 73.88730511353756, 9.555031159786623, 10.159484030408336, 28.753230155683998, 56.06243456487948, 6.88037298906764, 9.746176409734632, 76.48351885828143, 11.677681085553473, 7.8734026627145, 83.13458928723722, 12.47601459755017, 56.369336043489085, 38.79422547130149, 9.747888887236929, 35.18213256046644, 17.311572216100345, 6.654798398315644, 21.205006035962363, 146.73375709112386, 29.41553357409855, 14.77133786019102, 102.74532226215821, 10.428520389410174, 5.174315570063226, 23.882563835075153, 20.85725712502494, 18.230827555127817, 37.011115066793565, 24.077082261451345, 14.445864546598049, 121.86511622824248, 5.383008799121236, 10.342308830756494, 38.017352474591966, 5.153762521561403, 36.28312707571996, 104.38080187740073, 19.540276420624945, 8.280307167487964, 26.959179897435153, 57.23915855070114, 81.53097542794902, 16.041392747083027, 31.332282256928348, 17.116063325385987, 82.89291172723831, 50.942663206557036, 61.16509563160895, 13.241508996154332, 92.20384760567592, 84.32115200732882, 7.166571345499603, 32.25415245276355, 85.2918498392157, 23.58204239266529, 32.047140800123714, 15.663771838668453, 62.73409213837115, 80.13968246445067, 62.778846064423895, 57.58033072787298, 8.196694625907504, 76.40782581931734, 46.586682489976994, 114.35381760832637, 194.18203056426375, 34.59402256868242, 21.651933506693812, 23.826349217082235, 15.874130483702587, 14.036216805382255, 10.039572677623797, 54.37395246773005, 51.69284486157993, 25.494806326474627, 22.970498996319368, 29.049559360632642, 62.726248735489285, 5.985104622593393, 74.84996325167715, 70.97074195252529, 30.787962678658396, 14.898851547316111, 23.739679174309803, 131.0865084547394, 24.168222574394733, 14.37073089022526, 96.79395563644708, 27.262246127543477, 25.737726054337667, 20.340371735703464, 31.66474509112033, 46.57698551728102, 17.967807193836734, 59.292920229704265, 136.37043121808114, 17.54353663177534, 62.11276542206669, 21.48315950171234, 31.21179598230219, 5.391109844601892, 31.988454329878167, 71.11033831855057, 40.919044566978684, 25.297123589604777, 39.551484206563536, 80.19014551794595, 25.650910056096734, 5.336689840202448, 18.262915567326992, 118.13576015850285, 70.50607742782249, 29.788165855112435, 9.139218642279456, 5.149952051571178, 9.958332934713537, 95.75046259601393, 25.029793268646767, 214.91696354232363, 48.08420860780603, 40.20608745184056, 27.839534613914925, 58.88639165702904, 11.899245397706084, 30.47116918306545, 5.373817131441746, 24.87775343796074, 124.30687177435296, 128.6698593004463, 66.08445045038765, 64.38259088450616, 61.703001646313155, 51.814026213043114, 30.69536301464755, 36.91263224545558, 42.50890956154804, 35.226959027873065, 11.843565757265504, 43.42630273166789, 66.04710627837264, 30.746867034251878, 41.46239555718222, 52.91180541551975, 25.580055007486077, 133.75708739228938, 34.40303160956148, 20.85543539161491, 60.37378443553194, 30.03128253661255, 30.299903994369124, 16.528548427139636, 5.662333461269984, 17.27228407870122, 10.042434678355434, 13.73836268097965, 6.657999652977235, 32.808162388271285, 15.072944133872037, 8.86649357244533, 57.705624552604874, 18.76057017501909, 17.03910002189733, 43.21193784236018, 15.453375549136315, 71.76141430947406, 31.704235575903716, 8.02462762406903, 19.98030092372793, 180.17910151600915, 24.408890791191936, 5.4888252273976885, 89.3463587978369, 71.64351476969134, 55.37447538437138, 72.1857978352336, 42.679712633416365, 8.643517306005004, 5.561395765676327, 53.40116008827213, 17.19664765390621, 88.2193755643768, 10.918543297228254, 102.41681051040266, 27.202836355716975, 14.669160676969513, 102.80900398772879, 24.623023323470782, 112.11353867917133, 62.56346455250921, 86.35286182958929, 7.19712308406061, 7.154950940050308, 43.42000490054383, 66.37194296632163, 36.12332225708188, 43.400400107590535, 11.562255378660073, 214.58225222048264, 58.99461637531479, 74.54915486019145, 104.29956755666134, 103.65203008084146, 134.79304904453852, 12.899219204141861, 11.128346658731921, 5.118775767946365, 144.87369429286144, 6.0330410946785635, 7.084702728148356, 167.84000102526167, 7.691588857052687, 14.53307104785124, 183.54675656824554, 5.585668255883306, 20.06764064855352, 7.267865703719169, 80.23737809738822, 68.62163808254931, 25.999232394033875, 32.890423499539295, 60.279515643102314, 153.59472938346792, 20.483298019145177, 30.65913104945777, 77.22919222226835, 5.18866916473502, 21.917044705032993, 6.861636473787336, 77.52371763856812, 11.836137897767514, 86.5496144850733, 159.8892914441253, 28.392926426256608, 63.80776368308087, 8.897032944558788, 30.77797499672431, 22.684304195349007, 144.33985614999858, 11.257701021853373, 16.728772810005815, 183.94809856801078, 51.53671035664104, 60.40471128759087, 24.422617324156796, 36.242290016240105, 114.84915416001145, 18.281348353839597, 68.93676993387682, 18.21830886137796, 56.28448348867614, 68.90577910965548, 30.56781386306111, 29.583533394671633, 145.5582293654339, 20.39242065177016, 8.280572763948582, 23.860437961548836, 39.609649648667954, 65.42491419499525, 9.55149631981572, 94.89097058992613, 74.4860911168397, 100.31451595517484, 130.13482868631988, 10.512559695399421, 15.566394667900955, 81.39279492012433, 19.498866769125748, 35.731736579455664, 24.511853364985363, 34.08499860664441, 35.24380606237309, 9.240268095219559, 72.1596171400519, 7.785868806569296, 28.49040438915312, 34.6153542417438, 15.348391179493005, 56.34824251396226, 62.26588691460305, 11.290896566965255, 33.65689035365521, 118.13048166901218, 37.19803721435578, 6.604620932998921, 33.293768826754764, 238.91105150854133, 66.3434514520456, 27.173998390458138, 11.690694191586443, 23.567855531781866, 29.055017785531877, 54.299683535906794, 46.55121776490687, 29.22905125735331, 128.5067144673952, 42.7828020651135, 5.463892327023526, 133.58604478477812, 49.7303689759091, 31.66264001385943, 151.69191151898954, 93.28658861785371, 12.286680209557279, 43.96168087819643, 52.44665134370262, 18.30737912454215, 27.262277964126792, 62.80337129942751, 8.973898478598995, 14.160610039373534, 5.152519243478851, 29.39516946480492, 22.404832869442814, 61.74187752392386, 19.850046754710874, 55.992523970153066, 56.59263062667548, 21.2070585932029, 90.29082554855187, 24.817377674570402, 62.26553339617556, 10.179388207095963, 29.840027943612263, 25.148432215436568, 22.34846039962666, 194.80395689194546, 166.4605500712578, 8.695219865185168, 59.59849268858369, 14.212030126126832, 55.33781689202029, 28.84495210871959, 24.043688649642945, 26.631977377903716, 16.55183059250378, 20.222455248867366, 12.116994118988346, 59.85229158134261, 41.377537346054794, 26.64192663953893, 11.224798581029903, 6.977549470257569, 5.096085765150016, 12.268487781308231, 41.4307998743007, 20.090583717254727, 26.747580149438633, 76.85779785887442, 55.03010645074186, 24.90104027020149, 7.086458102945911, 7.569187318280818, 8.874315565176806, 7.695250335589253, 13.515898608904545, 77.18429990996823, 30.062439751363794, 88.56119199983307, 70.46745129449701, 29.8744497364807, 62.85347353867023, 71.29197731684675, 71.34306514086097, 50.62500220655267, 57.85440943410419, 22.638712254399703, 7.4309962488349965, 9.059045607160723, 45.0859702736866, 15.881033120140875, 84.58786933627641, 24.78097920718625, 58.25011522097172, 43.18242141078524, 34.870120130212115, 50.47119521091089, 68.80437733327997, 5.113183277464542, 14.91190639060965, 17.11870141624891, 77.97525458822486, 70.92655817137924, 18.410011486754357, 14.720262230913116, 7.287047973448669, 10.871291718559236, 57.54149552832244, 6.876003622381181, 7.968412210487223, 23.7507424976705, 16.10940164310589, 13.566601774286635, 38.87194449593684, 49.73875969952127, 35.45645998945266, 79.58453677538628, 9.236829289921983, 65.31867028738485, 52.190274251771115, 13.935945985869546, 42.446639098184924, 16.881420423257893, 82.46633068706194, 12.581882276225468, 10.598166632122277, 77.36410583768759, 13.733641874239257, 62.90548551977255, 11.777443829906712, 44.30551951064687, 34.462657662263986, 61.49912204134962, 28.94204199863049, 17.498114058070623, 28.700534699464015, 62.64861627083927, 10.691477407534023, 23.335499218411513, 87.3184803981755, 15.244022403213467, 9.093179826007699, 97.25350982175077, 67.77102157671955, 11.57517828026321, 19.127400663018648, 14.307975113364787, 5.091175813161763, 12.18189411108122, 8.495822580633002, 6.454646232699643, 39.57714115950007, 20.416618871769916, 30.983031334451482, 8.228442871664601, 5.159450684370454, 63.26234204182821, 27.5285880089019, 12.01176807770946, 53.851871123215744, 62.85649559307462, 12.478560575248734, 15.680807593072624, 54.943221485452305, 8.24350925923982, 189.51925173544402, 20.129206532026206, 7.113113983196451, 68.86700903224099, 5.210895440856173, 12.304246313266205, 16.001925737659978, 89.4118103059515, 6.952981295083042, 46.38527652941901, 24.418136025793643, 20.132192508550933, 39.479414730586335, 19.548667864534494, 20.418216381930918, 15.967772450331399, 45.44186675076523, 9.655117980023904, 15.107998916066029, 97.40949718756524, 7.135729478392734, 5.4658136897423555, 7.676627283868108, 16.921767717040005, 8.423557423532772, 6.137994272533657, 16.29234603103157, 29.20874494194414, 13.868315599188143, 34.26206872342423, 24.546954690578893, 54.280971801214115, 86.15760559422289, 13.455040187688954, 77.83820310843062, 24.137809355059577, 11.047177479278954, 78.6599350271893, 26.44893747934936, 32.90956474793774, 8.014399977262284, 39.97975007787519, 18.201640621549206, 65.13049776372034, 52.273259728624254, 76.00032508552368, 40.44337300506608, 72.11470513534665, 63.45351509360373, 13.155007101433918, 19.07850029166463, 11.629352502912301, 20.236715029391107, 73.47750879727143, 7.318210371663687, 6.470235280062387, 26.417195292758343, 6.680119250188373, 87.8891339336148, 9.891237912467195, 73.16290092732898, 23.500219874358503, 12.67375563074878, 73.1175212212965, 43.07414059986069, 5.743926808181919, 31.324349633784266, 8.661170368286205, 20.094764301948093, 13.713711528343914, 65.22718358105725, 13.502316154216112, 14.18886481410034, 6.709677157330401, 32.76343140317322, 133.99406267323621, 47.43691061492423, 42.817696147762305, 56.923591267712936, 6.962283887841462, 18.35591916934167, 18.69666289808284, 48.48546408647252, 55.345573634861836, 6.568993348418585, 25.533279001565415, 15.747297391074174, 5.390049023981349])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)