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 = 43866
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);
([3862600.0, 3979267.1875, 4087034.375, 4133768.75, 4318565.625, 4321907.8125, 4343157.8125, 4355064.0625, 4355626.5625, 4364679.6875, 4381337.5, 4384401.5625, 4396214.0625, 4476168.75, 4486793.75, 4493960.9375, 4517585.9375, 4524073.4375, 4569754.6875, 4575548.4375, 4580760.9375, 4584581.25, 4585990.625, 4595654.6875, 4612481.25, 4614454.6875, 4621801.5625, 4622539.0625, 4623907.8125, 4623979.6875, 4653740.625, 4663212.5, 4713223.4375, 4757078.125, 4758531.25, 4797179.6875, 4805151.5625, 4807101.5625, 4808800.0, 4809303.125, 4811125.0, 4817148.4375, 4844342.1875, 4848968.75, 4850312.5, 4888259.375, 4911323.4375, 4914681.25, 4919051.5625, 4934954.6875, 4939806.25, 4940601.5625, 4940656.25, 4955415.625, 4955451.5625, 4957335.9375, 4961868.75, 4966589.0625, 4981678.125, 4992470.3125, 4996915.625, 5001065.625, 5011359.375, 5014604.6875, 5015118.75, 5015368.75, 5016535.9375, 5016578.125, 5016731.25, 5017637.5, 5018109.375, 5019284.375, 5019651.5625, 5034879.6875, 5036162.5, 5036553.125, 5037040.625, 5038918.75, 5039106.25, 5039489.0625, 5039770.3125, 5041675.0, 5047054.6875, 5054171.875, 5054234.375, 5064928.125, 5066206.25, 5068498.4375, 5068953.125, 5069032.8125, 5070529.6875, 5070543.75, 5070584.375, 5071076.5625, 5072645.3125, 5072998.4375, 5073195.3125, 5092003.125, 5095356.25, 5105918.75, 5106700.0, 5107995.3125, 5109415.625, 5109445.3125, 5109992.1875, 5110078.125, 5110135.9375, 5110370.3125, 5110489.0625, 5110615.625, 5110673.4375, 5112773.4375, 5112804.6875, 5115187.5, 5116045.3125, 5116140.625, 5116140.625, 5116909.375, 5119240.625, 5124053.125, 5125056.25, 5125128.125, 5125537.5, 5126073.4375, 5126140.625, 5127953.125, 5134668.75, 5135028.125, 5136856.25, 5137165.625, 5137223.4375, 5137500.0, 5138432.8125, 5138853.125, 5138860.9375, 5138868.75, 5139048.4375, 5139057.8125, 5139212.5, 5139331.25, 5139509.375, 5139676.5625, 5139982.8125, 5140035.9375, 5140112.5, 5140334.375, 5140776.5625, 5141051.5625, 5141521.875, 5142256.25, 5143010.9375, 5143107.8125, 5145390.625, 5147039.0625, 5147098.4375, 5169242.1875, 5178893.75, 5186226.5625, 5194907.8125, 5213712.5, 5215975.0, 5216146.875, 5216181.25, 5216420.3125, 5216434.375, 5217042.1875, 5233975.0, 5234034.375, 5255371.875, 5282454.6875, 5285067.1875, 5289193.75, 5291700.0, 5292940.625, 5293376.5625, 5302750.0, 5342226.5625, 5348767.1875, 5372692.1875, 5412470.3125, 5418387.5, 5479260.9375, 5503931.25, 5517609.375, 5538609.375, 5587468.75, 5589815.625, 5597568.75, 5597589.0625, 5635546.875, 5640087.5, 5741042.1875, 5741078.125, 5751882.8125, 5797834.375, 5864878.125, 5897360.9375, 5953207.8125, 5957812.5, 5960584.375, 5983121.875, 5989935.9375, 5991009.375, 6001268.75, 6010135.9375, 6014906.25, 6040656.25, 6058089.0625, 6077487.5, 6100967.1875, 6113545.3125, 6118284.375, 6130020.3125, 6135818.75, 6171043.75, 6173962.5, 6179062.5, 6183042.1875, 6189400.0, 6196764.0625, 6198107.8125, 6207862.5, 6215082.8125, 6217312.5, 6229481.25, 6250023.4375, 6250032.8125, 6257992.1875, 6258051.5625, 6260210.9375, 6267571.875, 6304251.5625, 6309482.8125, 6309543.75, 6319576.5625, 6322145.3125, 6322181.25, 6324040.625, 6326746.875, 6336509.375, 6342437.5, 6352259.375, 6352775.0, 6362039.0625, 6362248.4375, 6368281.25, 6369189.0625, 6371087.5, 6373585.9375, 6380932.8125, 6394240.625, 6399618.75, 6410535.9375, 6410545.3125, 6413240.625, 6416323.4375, 6418545.3125, 6428387.5, 6430470.3125, 6433681.25, 6435090.625, 6435167.1875, 6436484.375, 6444385.9375, 6447400.0, 6448682.8125, 6448700.0, 6454200.0, 6457529.6875, 6460623.4375, 6467598.4375, 6470059.375, 6471812.5, 6481778.125, 6481903.125, 6485557.8125, 6487559.375, 6489975.0, 6490014.0625, 6491114.0625, 6492940.625, 6493853.125, 6496859.375, 6498262.5, 6502225.0, 6508539.0625, 6513260.9375, 6513512.5, 6514801.5625, 6518090.625, 6524323.4375, 6525020.3125, 6525940.625, 6528212.5, 6529785.9375, 6536190.625, 6538281.25, 6539114.0625, 6540150.0, 6540151.5625, 6540615.625, 6540617.1875, 6540714.0625, 6542154.6875, 6546798.4375, 6546917.1875, 6546948.4375, 6547179.6875, 6547817.1875, 6548545.3125, 6548573.4375, 6551942.1875, 6552953.125, 6557951.5625, 6561510.9375, 6562268.75, 6562459.375, 6563432.8125, 6564185.9375, 6567400.0, 6568051.5625, 6571945.3125, 6576687.5, 6585045.3125, 6585960.9375, 6585993.75, 6586079.6875, 6591973.4375, 6593929.6875, 6594459.375, 6594867.1875, 6605275.0, 6606850.0, 6615773.4375, 6623020.3125, 6625071.875, 6625198.4375, 6628996.875, 6629696.875, 6630925.0, 6634560.9375, 6643695.3125, 6645893.75, 6665757.8125, 6666648.4375, 6671643.75, 6673315.625, 6673625.0, 6680246.875, 6686115.625, 6689448.4375, 6690026.5625, 6710073.4375, 6711153.125, 6724657.8125, 6725851.5625, 6733646.875, 6733985.9375, 6736645.3125, 6738150.0, 6738395.3125, 6738909.375, 6739279.6875, 6739962.5, 6740268.75, 6740739.0625, 6741057.8125, 6742370.3125, 6742417.1875, 6743034.375, 6744151.5625, 6762662.5, 6764714.0625, 6766734.375, 6771920.3125, 6772339.0625, 6772350.0, 6772787.5, 6778729.6875, 6780512.5, 6780789.0625, 6793123.4375, 6809100.0, 6817171.875, 6824787.5, 6826395.3125, 6830545.3125, 6837921.875, 6844735.9375, 6861098.4375, 6861168.75, 6875943.75, 6889281.25, 6891171.875, 6893664.0625, 6893681.25, 6894607.8125, 6904993.75, 6913417.1875, 6917462.5, 6934946.875, 6935445.3125, 6940257.8125, 6960084.375, 6960428.125, 6960650.0, 6962103.125, 6962539.0625, 6964318.75, 6968901.5625, 6978320.3125, 6978865.625, 6982576.5625, 6988735.9375, 6991893.75, 6995650.0, 6996182.8125, 6996939.0625, 6997507.8125, 6997509.375, 7000198.4375, 7000220.3125, 7008189.0625, 7012421.875, 7012714.0625, 7013043.75, 7013500.0, 7014400.0, 7018082.8125, 7029426.5625, 7030426.5625, 7030912.5, 7030971.875, 7031157.8125, 7031432.8125, 7047139.0625, 7047982.8125, 7047984.375, 7048273.4375, 7048384.375, 7048868.75, 7050751.5625, 7052115.625, 7058409.375, 7061468.75, 7068865.625, 7076778.125, 7092384.375, 7110717.1875, 7120039.0625, 7130389.0625, 7130400.0, 7137335.9375, 7144703.125, 7144757.8125, 7147253.125, 7154360.9375, 7155145.3125, 7165195.3125, 7176370.3125, 7193917.1875, 7197621.875, 7198029.6875, 7205201.5625, 7226781.25, 7234670.3125, 7239090.625, 7243629.6875, 7247364.0625, 7248442.1875, 7267871.875, 7273971.875, 7288228.125, 7309335.9375, 7311043.75, 7311104.6875, 7314506.25, 7315379.6875, 7316137.5, 7337059.375, 7337065.625, 7390439.0625, 7435539.0625, 7449290.625, 7450670.3125, 7463337.5, 7463343.75, 7470698.4375, 7470970.3125, 7472193.75, 7476443.75, 7478492.1875, 7478557.8125, 7510765.625, 7518442.1875, 7542168.75, 7570048.4375, 7595409.375, 7596476.5625, 7597328.125, 7625543.75, 7650909.375, 7721596.875, 7722310.9375, 7727210.9375, 7764625.0, 7787993.75, 7791193.75, 7803498.4375, 7807184.375, 7810235.9375, 7812321.875, 7817117.1875, 7829871.875, 7846398.4375, 7847503.125, 7847651.5625, 7849807.8125, 7855193.75, 7855264.0625, 7860251.5625, 7892929.6875, 7898376.5625, 7901004.6875, 7924506.25, 7940000.0, 7958260.9375, 7959453.125, 7982820.3125, 7988326.5625, 7991940.625, 7991996.875, 7993918.75, 8005654.6875, 8014804.6875, 8018310.9375, 8019214.0625, 8037101.5625, 8037159.375, 8038393.75, 8043084.375, 8044379.6875, 8048943.75, 8050307.8125, 8054737.5, 8056631.25, 8057229.6875, 8101037.5, 8104162.5, 8113107.8125, 8129426.5625, 8132382.8125, 8135295.3125, 8169570.3125, 8169612.5, 8171153.125, 8172485.9375, 8190481.25, 8243443.75, 8249378.125, 8249789.0625, 8267040.625, 8311570.3125, 8313985.9375, 8315790.625, 8335160.9375, 8336242.1875, 8336887.5, 8361595.3125, 8364807.8125, 8370550.0, 8371848.4375, 8373148.4375, 8381481.25, 8391876.5625, 8422446.875, 8424526.5625, 8427985.9375, 8428779.6875, 8433364.0625, 8438934.375, 8519918.75, 8519971.875, 8525903.125, 8726356.25, 8763387.5, 8929082.8125, 8971112.5, 9008476.5625, 9011539.0625, 9046957.8125, 9068823.4375, 9069871.875, 9074462.5, 9090970.3125, 9103814.0625, 9112545.3125, 9130581.25, 9138062.5, 9148682.8125, 9191271.875, 9191565.625, 9193410.9375, 9219437.5, 9223810.9375, 9225489.0625, 9243531.25, 9247415.625, 9247484.375, 9249042.1875, 9259990.625, 9260784.375, 9265714.0625, 9266857.8125, 9267235.9375, 9267825.0, 9324629.6875, 9325446.875, 9325659.375, 9330354.6875, 9330698.4375, 9351007.8125, 9356687.5, 9371173.4375, 9377645.3125, 9381504.6875, 9383565.625, 9384579.6875, 9409228.125, 9409454.6875, 9409468.75, 9409582.8125, 9417356.25, 9447523.4375, 9471201.5625, 9471673.4375, 9476043.75, 9477607.8125, 9477642.1875, 9490775.0, 9493270.3125, 9525168.75, 9847228.125, 9851789.0625, 9972010.9375, 9974273.4375, 9996565.625, 10016776.5625, 10032243.75, 10032246.875, 10033557.8125, 10079851.5625, 10083731.25, 10094964.0625, 10099923.4375, 10100248.4375, 10107257.8125, 10107665.625, 10108107.8125, 10108204.6875, 10108870.3125, 10109131.25, 10109329.6875, 10110521.875, 10110553.125, 10110634.375, 10110637.5, 10111701.5625, 10111823.4375, 10112857.8125, 10116873.4375, 10117326.5625, 10127357.8125, 10137668.75, 10490954.6875, 535812517.1875], [61.730579369789076, 48.223354137396974, 15.20526725270402, 13.981023110054315, 29.87760757312946, 8.8068870672802, 53.287516968949696, 20.773058536750135, 36.16696175622107, 7.278502931360082, 69.8195187626727, 32.06916357656553, 91.40144064554639, 17.306456813279414, 71.82575058541576, 17.26115223846553, 36.82005593534155, 58.28714650826787, 29.247214909458286, 16.556835375852387, 94.37830766487349, 63.40329793728067, 53.301260756556076, 28.515792196900513, 22.679851731987373, 54.11880073029788, 24.45068907703016, 15.43737716502371, 7.577281942084358, 23.264514576495277, 10.97482042204417, 103.72941902316018, 69.91365078042362, 12.015182099028792, 63.58525176070736, 32.53758904004177, 103.56127250443365, 17.25103969670652, 21.622900519852916, 65.4284915646054, 13.610255899543073, 74.07951269563867, 78.9746427990302, 95.12001068498806, 10.451888871851619, 12.906966405783828, 7.515684239308004, 13.20426808065657, 53.90467035386394, 6.34446665664832, 41.593783170951454, 59.05137608731858, 5.985211047012849, 27.22913630424391, 8.493763433298101, 13.577839784342036, 24.18640466361277, 41.139901008758045, 91.63931870205877, 43.05841521024962, 78.51352714510219, 12.596699906364691, 11.862274207652536, 126.63405666656021, 9.94438284457542, 5.6382213073732546, 21.3045729749904, 7.03082262538515, 55.699376512973146, 21.134967891258377, 64.71416432393306, 115.73982724103377, 6.180512024523153, 54.19827351758506, 5.519729879265535, 39.58878192973954, 93.13651725991707, 47.32969561538789, 27.189196783452815, 31.387098532659103, 59.020388054702046, 79.2507291931169, 47.066529697195385, 143.9945885186026, 8.725691715282995, 65.70089970795175, 8.368557040809657, 38.2257745238806, 109.04133013783127, 56.61446660663913, 14.260438451569263, 17.885764170992623, 9.36848833732462, 11.759256207683247, 51.54628884884034, 15.839446681348281, 14.648683500482521, 81.51910588398297, 84.37002493283056, 9.245702555999216, 20.89688449327884, 6.171814498927349, 9.599322371419133, 10.40409590451822, 58.574494745199154, 5.030764341000801, 68.47480835064042, 15.000266253040886, 53.68535395660619, 14.800630748177747, 41.01064346390145, 9.12712243076058, 6.839159486730435, 26.842019915616735, 42.20707261622759, 35.93148929941992, 155.18396593260155, 10.60221257145808, 12.169251424823665, 13.410961640024007, 31.189443956598364, 47.304515957635275, 16.95262957251609, 31.73403784197447, 29.886110775353764, 66.07933477518047, 9.458399817358995, 23.655391710780297, 39.712076609853554, 20.806342781595227, 22.44993739040573, 43.285001593626326, 17.874312048511108, 9.468639434257565, 77.21973196699481, 10.42224688025424, 12.167689740481698, 48.56171347889771, 32.131744712693646, 5.128008457959243, 5.172687735865324, 62.01907078642327, 10.897371785278276, 42.28373698904865, 21.797581674781497, 9.445183980373956, 88.44222765667955, 7.239103518220504, 66.28352958259907, 83.54178917247957, 41.44306119551484, 81.89117104918886, 69.77496186655, 18.563519036306776, 9.191454933722552, 45.00750389329174, 88.01979931049114, 12.507170616932315, 18.234274268110326, 5.992280817301156, 10.829434963609172, 40.87166846994279, 32.95465470852147, 61.47347200341765, 9.153427359979961, 106.66700754916832, 53.13847538353467, 5.238805293416969, 8.23966495765203, 71.76374562771295, 5.113119931181222, 85.30577594787123, 46.95979836680726, 10.101147287581615, 100.06790626236443, 103.62546590057887, 5.288061662575674, 70.66371335672221, 88.92104396831037, 27.55024671969938, 17.675900795679045, 24.99696370170548, 11.807711105911824, 63.02792756097283, 7.925805942846159, 129.73757803980965, 115.01321023545769, 35.93698464243087, 19.147855659516466, 32.471437573622886, 12.360062693130613, 71.60958832773719, 46.863618105184614, 79.35380088287634, 28.483111377874053, 19.731372525701865, 13.338696102599023, 15.583814341811001, 44.84439561526891, 72.48995876389458, 218.95124201606626, 88.29097742282828, 110.02368492764231, 68.71600603673434, 59.05321111590324, 18.255761081425028, 8.273092663150978, 6.444423622082854, 11.968570930337682, 7.724159656255928, 27.841057381835878, 104.33712467205673, 11.462106339133992, 12.84168828848509, 65.37598566482066, 71.89401455411569, 5.414974034604237, 21.98536502795923, 82.22866062989874, 52.64603610257542, 21.338303581956367, 26.64318439502845, 77.04417609358892, 11.794104617377744, 109.53248715277338, 23.15473535517404, 7.420073722467204, 14.87005110966682, 5.244533443558159, 53.84804212457197, 10.91531745084978, 64.8133654698495, 31.565172727717496, 5.350243591334018, 19.988785946398895, 24.503686371718356, 9.716804283948901, 10.730386066795244, 27.459666727968482, 13.205529583902004, 18.418677884163525, 72.50527529739027, 30.365570785184836, 75.77208667303097, 24.075399494717622, 102.1657487479346, 13.802320103881083, 77.90057282903221, 16.838557146261092, 14.393087009053406, 11.6074338682422, 6.0718285224190325, 8.56149839962373, 7.329834011664996, 41.09463026792563, 121.48962804856919, 85.25947126929407, 60.08825528347694, 45.61652574287363, 50.84051923943387, 36.91125650911725, 32.11498207648948, 83.22881082285933, 14.517755401762706, 51.11155282468238, 25.330047256903022, 48.85817502138973, 45.68152032082069, 30.181022107928314, 72.53806656229108, 24.290739955686664, 59.2228569346157, 15.157472683568715, 24.412246767777454, 8.708342343333262, 101.05161961808629, 12.6818279486582, 39.08604279829844, 21.246873279726923, 13.028283091784859, 58.094650754922924, 21.337914665854708, 128.93703562507434, 24.23800095610355, 9.118804817790405, 6.047818101148177, 26.084261441212526, 45.333068429167, 17.76154547381708, 69.45959712216444, 63.98430500054071, 36.31862808164117, 28.9640825109195, 42.22738900890218, 41.18681810685485, 56.30445527598968, 24.23015793237512, 8.580936491118582, 10.710011393072676, 31.319274681951722, 36.740677176641576, 79.22652428895694, 83.85801107450524, 123.02659888141588, 25.45291106510391, 6.6003120484145885, 25.37959993260264, 62.296225374618864, 18.65264706888196, 42.86541123427597, 56.23518116085906, 77.91505951911773, 35.61965400930281, 12.986764137230292, 16.61274538307657, 7.682075335730902, 127.9640678428251, 28.936709998234544, 34.342105849768764, 30.15412218623032, 65.09603266959236, 15.688650737004744, 126.58777187957705, 9.0005751614303, 55.93454020344289, 34.73161008094415, 5.04338359456866, 11.715795681911231, 157.82052711413485, 19.504252519933612, 19.972752961652617, 11.723023883172063, 128.7009087298802, 51.44214591967797, 12.854145827417124, 90.00071724468377, 18.20750916385242, 42.302904960299955, 53.78025738594092, 9.140328519765541, 19.34818832214175, 39.9207867897276, 39.62879582296684, 37.60952911487856, 126.83127582860953, 8.33835221239202, 22.431608624716453, 21.086272045742575, 22.746047707759455, 40.491007644839705, 5.430735863960645, 101.00769044521091, 72.10993127852295, 22.75785109615057, 28.220112517481425, 14.327136913342457, 5.896760635786537, 9.271446841842039, 30.8769158144217, 7.153191961043938, 39.861171752483294, 15.11653859729901, 64.74602381584143, 27.01840532883351, 8.726637264477676, 56.952328243439396, 17.56655036647452, 12.236078558099925, 5.952318647711503, 5.534676756156285, 14.752929394495096, 15.372911106930875, 25.14682043715115, 12.31346510394708, 45.474799318379695, 28.102158067394786, 90.00540719219475, 21.658676528384206, 72.072381201717, 30.487563945580888, 24.555115297861207, 23.014581948932403, 37.98698636976693, 54.64118323186963, 107.0656874706244, 86.59438558268022, 30.829890313231704, 63.09626140804717, 7.8333189604034485, 10.346840813722324, 22.179946091522975, 13.615936569082141, 37.047275528887056, 19.446964474825027, 8.858553387425975, 11.439775069551928, 7.880245870814844, 15.578015561114382, 9.361808580993966, 19.767045034667653, 20.454542637760525, 76.58506325151703, 74.37606273264524, 10.225619584479105, 106.64777289345308, 7.999693516985066, 61.221941295517894, 36.07623861258619, 56.763285252860356, 25.196490665040216, 20.8713167473563, 36.10083995059748, 57.839942705430914, 20.60996143235243, 43.55151332708882, 23.0945368068699, 63.70677679557721, 8.290141509220348, 44.399001186510695, 41.881703554806734, 35.93480651895555, 39.79009598533669, 42.573157210624125, 26.632371319994455, 109.88974184229326, 11.024083676058078, 117.82100047708973, 154.12989079672278, 24.426219759491794, 152.7302157002492, 19.316585517140265, 6.8391741773212855, 8.445286318466541, 11.494699619130143, 63.40773352512696, 77.53625675997054, 16.722336522321818, 67.4793147193011, 114.94476157414326, 26.338395498068074, 7.622782401857462, 73.7995709200581, 65.88101933712124, 16.870487486741073, 92.83743446798024, 229.69474019658446, 39.29683955621515, 22.012955646277167, 28.10693353945151, 113.30646744647655, 164.09290827294507, 11.433611625337837, 19.390946013142216, 16.81448483606927, 5.082604255421736, 15.397183127648422, 33.63220620396156, 26.876544950873253, 68.28937473644041, 63.745136354713814, 28.634646864644168, 46.85764740630577, 52.830909040636044, 26.858852487381952, 66.11557651941142, 60.96823927614043, 59.85881440455127, 6.1570238904939565, 5.079603237050903, 6.572933635204267, 141.3910471895039, 17.833839867158765, 5.733418509969664, 33.56693838215604, 132.48797259861647, 19.801668201657108, 30.4397560883312, 12.84805224530508, 82.05945727395628, 28.35071327818005, 34.01331663929978, 29.320477133011078, 24.42253488219739, 155.64623068083208, 132.53883722411933, 30.295453693633927, 19.74057419197362, 8.792385057924562, 35.72537290810086, 12.594218922151232, 6.259744927709634, 16.602841157629637, 5.032665288337748, 34.99186781037157, 119.0094341636794, 6.399587376466778, 21.224829518155637, 64.48197170158335, 11.61595974111695, 14.56866216750861, 122.39595888332506, 17.995454509535463, 45.49510380443273, 67.59272937590914, 23.108238644303803, 6.1382265446787665, 53.79902104229021, 72.61577948932128, 52.45230115471795, 5.5055059120451375, 41.97164469007612, 27.51582335218675, 96.46941663923958, 88.3883180024952, 36.15140748755064, 185.5076619031931, 110.95142861150724, 76.2890693392739, 76.09063633221376, 5.166983389575908, 7.359926138544784, 55.164228081777544, 36.39429755335558, 48.90130337246897, 131.55657843021274, 21.733859869742002, 20.338290749005065, 30.07421291123829, 57.09408547942293, 5.701676502283257, 60.135476252272916, 32.491046287927, 54.05816027912783, 11.323303579211828, 49.32710295752727, 49.51388961188801, 112.8132908393257, 22.910951138450866, 5.070490242199125, 26.900182701594943, 53.73237526311319, 107.85546321330995, 23.70808576366379, 12.595883924333544, 35.67117078210234, 7.0640448999035454, 33.5982876391179, 24.836732500846797, 21.088197932078568, 30.260062375058325, 47.7631293809177, 66.08738590276428, 58.731198018757524, 146.00597035953444, 8.25580950200558, 70.34729137355163, 28.497720448015038, 36.63226136615714, 31.222535988338773, 5.665674900894059, 19.90198655424581, 15.856665956534949, 19.47795910997626, 10.921485093674768, 23.40790027024161, 5.355158333292589, 8.022416314295864, 87.3450577773408, 57.67415011156858, 13.80271535250499, 50.169500874922186, 11.041904057194687, 21.230512664101596, 151.44941274291034, 57.71101445413079, 16.417042152107918, 38.09155941628943, 55.6548870462807, 9.416910045766192, 5.0531352479069715, 39.45245823581494, 92.6197017735598, 17.032492279860367, 48.513704753347426, 13.138012817709667, 16.21513101353392, 43.686168398547125, 64.34531165015599, 36.24362298285726, 7.115637965153101, 30.700536713015218, 158.45750899636462, 90.39756825912099, 87.6456846868272, 6.175872922082695, 25.660783736871988, 13.818208392420315, 83.45359045919835, 42.33272960116086, 71.73893665242136, 43.11459858431879, 148.47960480981152, 5.59308165567319, 133.72717279455182, 61.369871928484116, 42.16265595010627, 19.46270464793236, 7.937189834498064, 130.78433580343727, 46.10557380232831, 8.023480231312336, 60.546793549042384, 14.09444739040739, 88.38143796926721, 11.602622602592488, 27.07879043195964, 11.822160457638445, 35.597386885408355, 123.65749981067931, 31.572016262239657, 60.97174393740824, 16.169729212878526, 25.029493975786217, 78.09176721354613, 12.20208574313815, 75.14029509746577, 30.16579318329695, 16.87403508488182, 17.35639238534709, 19.363942356541926, 6.797366901439606, 10.536594218599937, 91.40471737348084, 14.578262980450925, 45.22101724713076, 136.21667749142208, 6.401599496395278, 19.564644025293237, 57.00573166859808, 162.70103716014182, 84.0528027076331, 49.85870402683705, 25.14827671244069, 33.399246446039825, 36.813844454554356, 53.887075476995136, 6.016720100034734, 8.399397618488644, 44.01538259014254, 52.12184027879909, 64.24457395444753, 17.742198116445476, 21.505004734043997, 22.992561906928316, 56.76952034428838, 5.86958070003278, 13.661911845359016, 27.285096729424538, 15.53904600872319, 53.21185964708374, 13.373586004609324, 49.84541988889123, 8.233810708136899, 27.847489858416072, 11.845559703041914, 12.99774655888604, 58.151546800382306, 64.71085777608471, 24.892115376477363, 83.97508281992219, 15.76647395019558, 11.813324859478813, 13.638367506235745, 5.23169709021714, 57.73523655604864, 107.30746304912473])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3862600.0, 3979267.1875, 4087034.375, 4133768.75, 4318565.625, 4321907.8125, 4343157.8125, 4355064.0625, 4355626.5625, 4364679.6875, 4381337.5, 4384401.5625, 4396214.0625, 4476168.75, 4486793.75, 4493960.9375, 4517585.9375, 4524073.4375, 4569754.6875, 4575548.4375, 4580760.9375, 4584581.25, 4585990.625, 4595654.6875, 4612481.25, 4614454.6875, 4621801.5625, 4622539.0625, 4623907.8125, 4623979.6875, 4653740.625, 4663212.5, 4713223.4375, 4757078.125, 4758531.25, 4797179.6875, 4805151.5625, 4807101.5625, 4808800.0, 4809303.125, 4811125.0, 4817148.4375, 4844342.1875, 4848968.75, 4850312.5, 4888259.375, 4911323.4375, 4914681.25, 4919051.5625, 4934954.6875, 4939806.25, 4940601.5625, 4940656.25, 4955415.625, 4955451.5625, 4957335.9375, 4961868.75, 4966589.0625, 4981678.125, 4992470.3125, 4996915.625, 5001065.625, 5011359.375, 5014604.6875, 5015118.75, 5015368.75, 5016535.9375, 5016578.125, 5016731.25, 5017637.5, 5018109.375, 5019284.375, 5019651.5625, 5034879.6875, 5036162.5, 5036553.125, 5037040.625, 5038918.75, 5039106.25, 5039489.0625, 5039770.3125, 5041675.0, 5047054.6875, 5054171.875, 5054234.375, 5064928.125, 5066206.25, 5068498.4375, 5068953.125, 5069032.8125, 5070529.6875, 5070543.75, 5070584.375, 5071076.5625, 5072645.3125, 5072998.4375, 5073195.3125, 5092003.125, 5095356.25, 5105918.75, 5106700.0, 5107995.3125, 5109415.625, 5109445.3125, 5109992.1875, 5110078.125, 5110135.9375, 5110370.3125, 5110489.0625, 5110615.625, 5110673.4375, 5112773.4375, 5112804.6875, 5115187.5, 5116045.3125, 5116140.625, 5116140.625, 5116909.375, 5119240.625, 5124053.125, 5125056.25, 5125128.125, 5125537.5, 5126073.4375, 5126140.625, 5127953.125, 5134668.75, 5135028.125, 5136856.25, 5137165.625, 5137223.4375, 5137500.0, 5138432.8125, 5138853.125, 5138860.9375, 5138868.75, 5139048.4375, 5139057.8125, 5139212.5, 5139331.25, 5139509.375, 5139676.5625, 5139982.8125, 5140035.9375, 5140112.5, 5140334.375, 5140776.5625, 5141051.5625, 5141521.875, 5142256.25, 5143010.9375, 5143107.8125, 5145390.625, 5147039.0625, 5147098.4375, 5169242.1875, 5178893.75, 5186226.5625, 5194907.8125, 5213712.5, 5215975.0, 5216146.875, 5216181.25, 5216420.3125, 5216434.375, 5217042.1875, 5233975.0, 5234034.375, 5255371.875, 5282454.6875, 5285067.1875, 5289193.75, 5291700.0, 5292940.625, 5293376.5625, 5302750.0, 5342226.5625, 5348767.1875, 5372692.1875, 5412470.3125, 5418387.5, 5479260.9375, 5503931.25, 5517609.375, 5538609.375, 5587468.75, 5589815.625, 5597568.75, 5597589.0625, 5635546.875, 5640087.5, 5741042.1875, 5741078.125, 5751882.8125, 5797834.375, 5864878.125, 5897360.9375, 5953207.8125, 5957812.5, 5960584.375, 5983121.875, 5989935.9375, 5991009.375, 6001268.75, 6010135.9375, 6014906.25, 6040656.25, 6058089.0625, 6077487.5, 6100967.1875, 6113545.3125, 6118284.375, 6130020.3125, 6135818.75, 6171043.75, 6173962.5, 6179062.5, 6183042.1875, 6189400.0, 6196764.0625, 6198107.8125, 6207862.5, 6215082.8125, 6217312.5, 6229481.25, 6250023.4375, 6250032.8125, 6257992.1875, 6258051.5625, 6260210.9375, 6267571.875, 6304251.5625, 6309482.8125, 6309543.75, 6319576.5625, 6322145.3125, 6322181.25, 6324040.625, 6326746.875, 6336509.375, 6342437.5, 6352259.375, 6352775.0, 6362039.0625, 6362248.4375, 6368281.25, 6369189.0625, 6371087.5, 6373585.9375, 6380932.8125, 6394240.625, 6399618.75, 6410535.9375, 6410545.3125, 6413240.625, 6416323.4375, 6418545.3125, 6428387.5, 6430470.3125, 6433681.25, 6435090.625, 6435167.1875, 6436484.375, 6444385.9375, 6447400.0, 6448682.8125, 6448700.0, 6454200.0, 6457529.6875, 6460623.4375, 6467598.4375, 6470059.375, 6471812.5, 6481778.125, 6481903.125, 6485557.8125, 6487559.375, 6489975.0, 6490014.0625, 6491114.0625, 6492940.625, 6493853.125, 6496859.375, 6498262.5, 6502225.0, 6508539.0625, 6513260.9375, 6513512.5, 6514801.5625, 6518090.625, 6524323.4375, 6525020.3125, 6525940.625, 6528212.5, 6529785.9375, 6536190.625, 6538281.25, 6539114.0625, 6540150.0, 6540151.5625, 6540615.625, 6540617.1875, 6540714.0625, 6542154.6875, 6546798.4375, 6546917.1875, 6546948.4375, 6547179.6875, 6547817.1875, 6548545.3125, 6548573.4375, 6551942.1875, 6552953.125, 6557951.5625, 6561510.9375, 6562268.75, 6562459.375, 6563432.8125, 6564185.9375, 6567400.0, 6568051.5625, 6571945.3125, 6576687.5, 6585045.3125, 6585960.9375, 6585993.75, 6586079.6875, 6591973.4375, 6593929.6875, 6594459.375, 6594867.1875, 6605275.0, 6606850.0, 6615773.4375, 6623020.3125, 6625071.875, 6625198.4375, 6628996.875, 6629696.875, 6630925.0, 6634560.9375, 6643695.3125, 6645893.75, 6665757.8125, 6666648.4375, 6671643.75, 6673315.625, 6673625.0, 6680246.875, 6686115.625, 6689448.4375, 6690026.5625, 6710073.4375, 6711153.125, 6724657.8125, 6725851.5625, 6733646.875, 6733985.9375, 6736645.3125, 6738150.0, 6738395.3125, 6738909.375, 6739279.6875, 6739962.5, 6740268.75, 6740739.0625, 6741057.8125, 6742370.3125, 6742417.1875, 6743034.375, 6744151.5625, 6762662.5, 6764714.0625, 6766734.375, 6771920.3125, 6772339.0625, 6772350.0, 6772787.5, 6778729.6875, 6780512.5, 6780789.0625, 6793123.4375, 6809100.0, 6817171.875, 6824787.5, 6826395.3125, 6830545.3125, 6837921.875, 6844735.9375, 6861098.4375, 6861168.75, 6875943.75, 6889281.25, 6891171.875, 6893664.0625, 6893681.25, 6894607.8125, 6904993.75, 6913417.1875, 6917462.5, 6934946.875, 6935445.3125, 6940257.8125, 6960084.375, 6960428.125, 6960650.0, 6962103.125, 6962539.0625, 6964318.75, 6968901.5625, 6978320.3125, 6978865.625, 6982576.5625, 6988735.9375, 6991893.75, 6995650.0, 6996182.8125, 6996939.0625, 6997507.8125, 6997509.375, 7000198.4375, 7000220.3125, 7008189.0625, 7012421.875, 7012714.0625, 7013043.75, 7013500.0, 7014400.0, 7018082.8125, 7029426.5625, 7030426.5625, 7030912.5, 7030971.875, 7031157.8125, 7031432.8125, 7047139.0625, 7047982.8125, 7047984.375, 7048273.4375, 7048384.375, 7048868.75, 7050751.5625, 7052115.625, 7058409.375, 7061468.75, 7068865.625, 7076778.125, 7092384.375, 7110717.1875, 7120039.0625, 7130389.0625, 7130400.0, 7137335.9375, 7144703.125, 7144757.8125, 7147253.125, 7154360.9375, 7155145.3125, 7165195.3125, 7176370.3125, 7193917.1875, 7197621.875, 7198029.6875, 7205201.5625, 7226781.25, 7234670.3125, 7239090.625, 7243629.6875, 7247364.0625, 7248442.1875, 7267871.875, 7273971.875, 7288228.125, 7309335.9375, 7311043.75, 7311104.6875, 7314506.25, 7315379.6875, 7316137.5, 7337059.375, 7337065.625, 7390439.0625, 7435539.0625, 7449290.625, 7450670.3125, 7463337.5, 7463343.75, 7470698.4375, 7470970.3125, 7472193.75, 7476443.75, 7478492.1875, 7478557.8125, 7510765.625, 7518442.1875, 7542168.75, 7570048.4375, 7595409.375, 7596476.5625, 7597328.125, 7625543.75, 7650909.375, 7721596.875, 7722310.9375, 7727210.9375, 7764625.0, 7787993.75, 7791193.75, 7803498.4375, 7807184.375, 7810235.9375, 7812321.875, 7817117.1875, 7829871.875, 7846398.4375, 7847503.125, 7847651.5625, 7849807.8125, 7855193.75, 7855264.0625, 7860251.5625, 7892929.6875, 7898376.5625, 7901004.6875, 7924506.25, 7940000.0, 7958260.9375, 7959453.125, 7982820.3125, 7988326.5625, 7991940.625, 7991996.875, 7993918.75, 8005654.6875, 8014804.6875, 8018310.9375, 8019214.0625, 8037101.5625, 8037159.375, 8038393.75, 8043084.375, 8044379.6875, 8048943.75, 8050307.8125, 8054737.5, 8056631.25, 8057229.6875, 8101037.5, 8104162.5, 8113107.8125, 8129426.5625, 8132382.8125, 8135295.3125, 8169570.3125, 8169612.5, 8171153.125, 8172485.9375, 8190481.25, 8243443.75, 8249378.125, 8249789.0625, 8267040.625, 8311570.3125, 8313985.9375, 8315790.625, 8335160.9375, 8336242.1875, 8336887.5, 8361595.3125, 8364807.8125, 8370550.0, 8371848.4375, 8373148.4375, 8381481.25, 8391876.5625, 8422446.875, 8424526.5625, 8427985.9375, 8428779.6875, 8433364.0625, 8438934.375, 8519918.75, 8519971.875, 8525903.125, 8726356.25, 8763387.5, 8929082.8125, 8971112.5, 9008476.5625, 9011539.0625, 9046957.8125, 9068823.4375, 9069871.875, 9074462.5, 9090970.3125, 9103814.0625, 9112545.3125, 9130581.25, 9138062.5, 9148682.8125, 9191271.875, 9191565.625, 9193410.9375, 9219437.5, 9223810.9375, 9225489.0625, 9243531.25, 9247415.625, 9247484.375, 9249042.1875, 9259990.625, 9260784.375, 9265714.0625, 9266857.8125, 9267235.9375, 9267825.0, 9324629.6875, 9325446.875, 9325659.375, 9330354.6875, 9330698.4375, 9351007.8125, 9356687.5, 9371173.4375, 9377645.3125, 9381504.6875, 9383565.625, 9384579.6875, 9409228.125, 9409454.6875, 9409468.75, 9409582.8125, 9417356.25, 9447523.4375, 9471201.5625, 9471673.4375, 9476043.75, 9477607.8125, 9477642.1875, 9490775.0, 9493270.3125, 9525168.75, 9847228.125, 9851789.0625, 9972010.9375, 9974273.4375, 9996565.625, 10016776.5625, 10032243.75, 10032246.875, 10033557.8125, 10079851.5625, 10083731.25, 10094964.0625, 10099923.4375, 10100248.4375, 10107257.8125, 10107665.625, 10108107.8125, 10108204.6875, 10108870.3125, 10109131.25, 10109329.6875, 10110521.875, 10110553.125, 10110634.375, 10110637.5, 10111701.5625, 10111823.4375, 10112857.8125, 10116873.4375, 10117326.5625, 10127357.8125, 10137668.75, 10490954.6875, 535812517.1875], [61.730579369789076, 48.223354137396974, 15.20526725270402, 13.981023110054315, 29.87760757312946, 8.8068870672802, 53.287516968949696, 20.773058536750135, 36.16696175622107, 7.278502931360082, 69.8195187626727, 32.06916357656553, 91.40144064554639, 17.306456813279414, 71.82575058541576, 17.26115223846553, 36.82005593534155, 58.28714650826787, 29.247214909458286, 16.556835375852387, 94.37830766487349, 63.40329793728067, 53.301260756556076, 28.515792196900513, 22.679851731987373, 54.11880073029788, 24.45068907703016, 15.43737716502371, 7.577281942084358, 23.264514576495277, 10.97482042204417, 103.72941902316018, 69.91365078042362, 12.015182099028792, 63.58525176070736, 32.53758904004177, 103.56127250443365, 17.25103969670652, 21.622900519852916, 65.4284915646054, 13.610255899543073, 74.07951269563867, 78.9746427990302, 95.12001068498806, 10.451888871851619, 12.906966405783828, 7.515684239308004, 13.20426808065657, 53.90467035386394, 6.34446665664832, 41.593783170951454, 59.05137608731858, 5.985211047012849, 27.22913630424391, 8.493763433298101, 13.577839784342036, 24.18640466361277, 41.139901008758045, 91.63931870205877, 43.05841521024962, 78.51352714510219, 12.596699906364691, 11.862274207652536, 126.63405666656021, 9.94438284457542, 5.6382213073732546, 21.3045729749904, 7.03082262538515, 55.699376512973146, 21.134967891258377, 64.71416432393306, 115.73982724103377, 6.180512024523153, 54.19827351758506, 5.519729879265535, 39.58878192973954, 93.13651725991707, 47.32969561538789, 27.189196783452815, 31.387098532659103, 59.020388054702046, 79.2507291931169, 47.066529697195385, 143.9945885186026, 8.725691715282995, 65.70089970795175, 8.368557040809657, 38.2257745238806, 109.04133013783127, 56.61446660663913, 14.260438451569263, 17.885764170992623, 9.36848833732462, 11.759256207683247, 51.54628884884034, 15.839446681348281, 14.648683500482521, 81.51910588398297, 84.37002493283056, 9.245702555999216, 20.89688449327884, 6.171814498927349, 9.599322371419133, 10.40409590451822, 58.574494745199154, 5.030764341000801, 68.47480835064042, 15.000266253040886, 53.68535395660619, 14.800630748177747, 41.01064346390145, 9.12712243076058, 6.839159486730435, 26.842019915616735, 42.20707261622759, 35.93148929941992, 155.18396593260155, 10.60221257145808, 12.169251424823665, 13.410961640024007, 31.189443956598364, 47.304515957635275, 16.95262957251609, 31.73403784197447, 29.886110775353764, 66.07933477518047, 9.458399817358995, 23.655391710780297, 39.712076609853554, 20.806342781595227, 22.44993739040573, 43.285001593626326, 17.874312048511108, 9.468639434257565, 77.21973196699481, 10.42224688025424, 12.167689740481698, 48.56171347889771, 32.131744712693646, 5.128008457959243, 5.172687735865324, 62.01907078642327, 10.897371785278276, 42.28373698904865, 21.797581674781497, 9.445183980373956, 88.44222765667955, 7.239103518220504, 66.28352958259907, 83.54178917247957, 41.44306119551484, 81.89117104918886, 69.77496186655, 18.563519036306776, 9.191454933722552, 45.00750389329174, 88.01979931049114, 12.507170616932315, 18.234274268110326, 5.992280817301156, 10.829434963609172, 40.87166846994279, 32.95465470852147, 61.47347200341765, 9.153427359979961, 106.66700754916832, 53.13847538353467, 5.238805293416969, 8.23966495765203, 71.76374562771295, 5.113119931181222, 85.30577594787123, 46.95979836680726, 10.101147287581615, 100.06790626236443, 103.62546590057887, 5.288061662575674, 70.66371335672221, 88.92104396831037, 27.55024671969938, 17.675900795679045, 24.99696370170548, 11.807711105911824, 63.02792756097283, 7.925805942846159, 129.73757803980965, 115.01321023545769, 35.93698464243087, 19.147855659516466, 32.471437573622886, 12.360062693130613, 71.60958832773719, 46.863618105184614, 79.35380088287634, 28.483111377874053, 19.731372525701865, 13.338696102599023, 15.583814341811001, 44.84439561526891, 72.48995876389458, 218.95124201606626, 88.29097742282828, 110.02368492764231, 68.71600603673434, 59.05321111590324, 18.255761081425028, 8.273092663150978, 6.444423622082854, 11.968570930337682, 7.724159656255928, 27.841057381835878, 104.33712467205673, 11.462106339133992, 12.84168828848509, 65.37598566482066, 71.89401455411569, 5.414974034604237, 21.98536502795923, 82.22866062989874, 52.64603610257542, 21.338303581956367, 26.64318439502845, 77.04417609358892, 11.794104617377744, 109.53248715277338, 23.15473535517404, 7.420073722467204, 14.87005110966682, 5.244533443558159, 53.84804212457197, 10.91531745084978, 64.8133654698495, 31.565172727717496, 5.350243591334018, 19.988785946398895, 24.503686371718356, 9.716804283948901, 10.730386066795244, 27.459666727968482, 13.205529583902004, 18.418677884163525, 72.50527529739027, 30.365570785184836, 75.77208667303097, 24.075399494717622, 102.1657487479346, 13.802320103881083, 77.90057282903221, 16.838557146261092, 14.393087009053406, 11.6074338682422, 6.0718285224190325, 8.56149839962373, 7.329834011664996, 41.09463026792563, 121.48962804856919, 85.25947126929407, 60.08825528347694, 45.61652574287363, 50.84051923943387, 36.91125650911725, 32.11498207648948, 83.22881082285933, 14.517755401762706, 51.11155282468238, 25.330047256903022, 48.85817502138973, 45.68152032082069, 30.181022107928314, 72.53806656229108, 24.290739955686664, 59.2228569346157, 15.157472683568715, 24.412246767777454, 8.708342343333262, 101.05161961808629, 12.6818279486582, 39.08604279829844, 21.246873279726923, 13.028283091784859, 58.094650754922924, 21.337914665854708, 128.93703562507434, 24.23800095610355, 9.118804817790405, 6.047818101148177, 26.084261441212526, 45.333068429167, 17.76154547381708, 69.45959712216444, 63.98430500054071, 36.31862808164117, 28.9640825109195, 42.22738900890218, 41.18681810685485, 56.30445527598968, 24.23015793237512, 8.580936491118582, 10.710011393072676, 31.319274681951722, 36.740677176641576, 79.22652428895694, 83.85801107450524, 123.02659888141588, 25.45291106510391, 6.6003120484145885, 25.37959993260264, 62.296225374618864, 18.65264706888196, 42.86541123427597, 56.23518116085906, 77.91505951911773, 35.61965400930281, 12.986764137230292, 16.61274538307657, 7.682075335730902, 127.9640678428251, 28.936709998234544, 34.342105849768764, 30.15412218623032, 65.09603266959236, 15.688650737004744, 126.58777187957705, 9.0005751614303, 55.93454020344289, 34.73161008094415, 5.04338359456866, 11.715795681911231, 157.82052711413485, 19.504252519933612, 19.972752961652617, 11.723023883172063, 128.7009087298802, 51.44214591967797, 12.854145827417124, 90.00071724468377, 18.20750916385242, 42.302904960299955, 53.78025738594092, 9.140328519765541, 19.34818832214175, 39.9207867897276, 39.62879582296684, 37.60952911487856, 126.83127582860953, 8.33835221239202, 22.431608624716453, 21.086272045742575, 22.746047707759455, 40.491007644839705, 5.430735863960645, 101.00769044521091, 72.10993127852295, 22.75785109615057, 28.220112517481425, 14.327136913342457, 5.896760635786537, 9.271446841842039, 30.8769158144217, 7.153191961043938, 39.861171752483294, 15.11653859729901, 64.74602381584143, 27.01840532883351, 8.726637264477676, 56.952328243439396, 17.56655036647452, 12.236078558099925, 5.952318647711503, 5.534676756156285, 14.752929394495096, 15.372911106930875, 25.14682043715115, 12.31346510394708, 45.474799318379695, 28.102158067394786, 90.00540719219475, 21.658676528384206, 72.072381201717, 30.487563945580888, 24.555115297861207, 23.014581948932403, 37.98698636976693, 54.64118323186963, 107.0656874706244, 86.59438558268022, 30.829890313231704, 63.09626140804717, 7.8333189604034485, 10.346840813722324, 22.179946091522975, 13.615936569082141, 37.047275528887056, 19.446964474825027, 8.858553387425975, 11.439775069551928, 7.880245870814844, 15.578015561114382, 9.361808580993966, 19.767045034667653, 20.454542637760525, 76.58506325151703, 74.37606273264524, 10.225619584479105, 106.64777289345308, 7.999693516985066, 61.221941295517894, 36.07623861258619, 56.763285252860356, 25.196490665040216, 20.8713167473563, 36.10083995059748, 57.839942705430914, 20.60996143235243, 43.55151332708882, 23.0945368068699, 63.70677679557721, 8.290141509220348, 44.399001186510695, 41.881703554806734, 35.93480651895555, 39.79009598533669, 42.573157210624125, 26.632371319994455, 109.88974184229326, 11.024083676058078, 117.82100047708973, 154.12989079672278, 24.426219759491794, 152.7302157002492, 19.316585517140265, 6.8391741773212855, 8.445286318466541, 11.494699619130143, 63.40773352512696, 77.53625675997054, 16.722336522321818, 67.4793147193011, 114.94476157414326, 26.338395498068074, 7.622782401857462, 73.7995709200581, 65.88101933712124, 16.870487486741073, 92.83743446798024, 229.69474019658446, 39.29683955621515, 22.012955646277167, 28.10693353945151, 113.30646744647655, 164.09290827294507, 11.433611625337837, 19.390946013142216, 16.81448483606927, 5.082604255421736, 15.397183127648422, 33.63220620396156, 26.876544950873253, 68.28937473644041, 63.745136354713814, 28.634646864644168, 46.85764740630577, 52.830909040636044, 26.858852487381952, 66.11557651941142, 60.96823927614043, 59.85881440455127, 6.1570238904939565, 5.079603237050903, 6.572933635204267, 141.3910471895039, 17.833839867158765, 5.733418509969664, 33.56693838215604, 132.48797259861647, 19.801668201657108, 30.4397560883312, 12.84805224530508, 82.05945727395628, 28.35071327818005, 34.01331663929978, 29.320477133011078, 24.42253488219739, 155.64623068083208, 132.53883722411933, 30.295453693633927, 19.74057419197362, 8.792385057924562, 35.72537290810086, 12.594218922151232, 6.259744927709634, 16.602841157629637, 5.032665288337748, 34.99186781037157, 119.0094341636794, 6.399587376466778, 21.224829518155637, 64.48197170158335, 11.61595974111695, 14.56866216750861, 122.39595888332506, 17.995454509535463, 45.49510380443273, 67.59272937590914, 23.108238644303803, 6.1382265446787665, 53.79902104229021, 72.61577948932128, 52.45230115471795, 5.5055059120451375, 41.97164469007612, 27.51582335218675, 96.46941663923958, 88.3883180024952, 36.15140748755064, 185.5076619031931, 110.95142861150724, 76.2890693392739, 76.09063633221376, 5.166983389575908, 7.359926138544784, 55.164228081777544, 36.39429755335558, 48.90130337246897, 131.55657843021274, 21.733859869742002, 20.338290749005065, 30.07421291123829, 57.09408547942293, 5.701676502283257, 60.135476252272916, 32.491046287927, 54.05816027912783, 11.323303579211828, 49.32710295752727, 49.51388961188801, 112.8132908393257, 22.910951138450866, 5.070490242199125, 26.900182701594943, 53.73237526311319, 107.85546321330995, 23.70808576366379, 12.595883924333544, 35.67117078210234, 7.0640448999035454, 33.5982876391179, 24.836732500846797, 21.088197932078568, 30.260062375058325, 47.7631293809177, 66.08738590276428, 58.731198018757524, 146.00597035953444, 8.25580950200558, 70.34729137355163, 28.497720448015038, 36.63226136615714, 31.222535988338773, 5.665674900894059, 19.90198655424581, 15.856665956534949, 19.47795910997626, 10.921485093674768, 23.40790027024161, 5.355158333292589, 8.022416314295864, 87.3450577773408, 57.67415011156858, 13.80271535250499, 50.169500874922186, 11.041904057194687, 21.230512664101596, 151.44941274291034, 57.71101445413079, 16.417042152107918, 38.09155941628943, 55.6548870462807, 9.416910045766192, 5.0531352479069715, 39.45245823581494, 92.6197017735598, 17.032492279860367, 48.513704753347426, 13.138012817709667, 16.21513101353392, 43.686168398547125, 64.34531165015599, 36.24362298285726, 7.115637965153101, 30.700536713015218, 158.45750899636462, 90.39756825912099, 87.6456846868272, 6.175872922082695, 25.660783736871988, 13.818208392420315, 83.45359045919835, 42.33272960116086, 71.73893665242136, 43.11459858431879, 148.47960480981152, 5.59308165567319, 133.72717279455182, 61.369871928484116, 42.16265595010627, 19.46270464793236, 7.937189834498064, 130.78433580343727, 46.10557380232831, 8.023480231312336, 60.546793549042384, 14.09444739040739, 88.38143796926721, 11.602622602592488, 27.07879043195964, 11.822160457638445, 35.597386885408355, 123.65749981067931, 31.572016262239657, 60.97174393740824, 16.169729212878526, 25.029493975786217, 78.09176721354613, 12.20208574313815, 75.14029509746577, 30.16579318329695, 16.87403508488182, 17.35639238534709, 19.363942356541926, 6.797366901439606, 10.536594218599937, 91.40471737348084, 14.578262980450925, 45.22101724713076, 136.21667749142208, 6.401599496395278, 19.564644025293237, 57.00573166859808, 162.70103716014182, 84.0528027076331, 49.85870402683705, 25.14827671244069, 33.399246446039825, 36.813844454554356, 53.887075476995136, 6.016720100034734, 8.399397618488644, 44.01538259014254, 52.12184027879909, 64.24457395444753, 17.742198116445476, 21.505004734043997, 22.992561906928316, 56.76952034428838, 5.86958070003278, 13.661911845359016, 27.285096729424538, 15.53904600872319, 53.21185964708374, 13.373586004609324, 49.84541988889123, 8.233810708136899, 27.847489858416072, 11.845559703041914, 12.99774655888604, 58.151546800382306, 64.71085777608471, 24.892115376477363, 83.97508281992219, 15.76647395019558, 11.813324859478813, 13.638367506235745, 5.23169709021714, 57.73523655604864, 107.30746304912473])
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);
([3862600.0, 3979267.1875, 4087034.375, 4133768.75, 4318565.625, 4321907.8125, 4343157.8125, 4355064.0625, 4355626.5625, 4364679.6875, 4381337.5, 4384401.5625, 4396214.0625, 4476168.75, 4486793.75, 4493960.9375, 4517585.9375, 4524073.4375, 4569754.6875, 4575548.4375, 4580760.9375, 4584581.25, 4585990.625, 4595654.6875, 4612481.25, 4614454.6875, 4621801.5625, 4622539.0625, 4623907.8125, 4623979.6875, 4653740.625, 4663212.5, 4713223.4375, 4757078.125, 4758531.25, 4797179.6875, 4805151.5625, 4807101.5625, 4808800.0, 4809303.125, 4811125.0, 4817148.4375, 4844342.1875, 4848968.75, 4850312.5, 4888259.375, 4911323.4375, 4914681.25, 4919051.5625, 4934954.6875, 4939806.25, 4940601.5625, 4940656.25, 4955415.625, 4955451.5625, 4957335.9375, 4961868.75, 4966589.0625, 4981678.125, 4992470.3125, 4996915.625, 5001065.625, 5011359.375, 5014604.6875, 5015118.75, 5015368.75, 5016535.9375, 5016578.125, 5016731.25, 5017637.5, 5018109.375, 5019284.375, 5019651.5625, 5034879.6875, 5036162.5, 5036553.125, 5037040.625, 5038918.75, 5039106.25, 5039489.0625, 5039770.3125, 5041675.0, 5047054.6875, 5054171.875, 5054234.375, 5064928.125, 5066206.25, 5068498.4375, 5068953.125, 5069032.8125, 5070529.6875, 5070543.75, 5070584.375, 5071076.5625, 5072645.3125, 5072998.4375, 5073195.3125, 5092003.125, 5095356.25, 5105918.75, 5106700.0, 5107995.3125, 5109415.625, 5109445.3125, 5109992.1875, 5110078.125, 5110135.9375, 5110370.3125, 5110489.0625, 5110615.625, 5110673.4375, 5112773.4375, 5112804.6875, 5115187.5, 5116045.3125, 5116140.625, 5116140.625, 5116909.375, 5119240.625, 5124053.125, 5125056.25, 5125128.125, 5125537.5, 5126073.4375, 5126140.625, 5127953.125, 5134668.75, 5135028.125, 5136856.25, 5137165.625, 5137223.4375, 5137500.0, 5138432.8125, 5138853.125, 5138860.9375, 5138868.75, 5139048.4375, 5139057.8125, 5139212.5, 5139331.25, 5139509.375, 5139676.5625, 5139982.8125, 5140035.9375, 5140112.5, 5140334.375, 5140776.5625, 5141051.5625, 5141521.875, 5142256.25, 5143010.9375, 5143107.8125, 5145390.625, 5147039.0625, 5147098.4375, 5169242.1875, 5178893.75, 5186226.5625, 5194907.8125, 5213712.5, 5215975.0, 5216146.875, 5216181.25, 5216420.3125, 5216434.375, 5217042.1875, 5233975.0, 5234034.375, 5255371.875, 5282454.6875, 5285067.1875, 5289193.75, 5291700.0, 5292940.625, 5293376.5625, 5302750.0, 5342226.5625, 5348767.1875, 5372692.1875, 5412470.3125, 5418387.5, 5479260.9375, 5503931.25, 5517609.375, 5538609.375, 5587468.75, 5589815.625, 5597568.75, 5597589.0625, 5635546.875, 5640087.5, 5741042.1875, 5741078.125, 5751882.8125, 5797834.375, 5864878.125, 5897360.9375, 5953207.8125, 5957812.5, 5960584.375, 5983121.875, 5989935.9375, 5991009.375, 6001268.75, 6010135.9375, 6014906.25, 6040656.25, 6058089.0625, 6077487.5, 6100967.1875, 6113545.3125, 6118284.375, 6130020.3125, 6135818.75, 6171043.75, 6173962.5, 6179062.5, 6183042.1875, 6189400.0, 6196764.0625, 6198107.8125, 6207862.5, 6215082.8125, 6217312.5, 6229481.25, 6250023.4375, 6250032.8125, 6257992.1875, 6258051.5625, 6260210.9375, 6267571.875, 6304251.5625, 6309482.8125, 6309543.75, 6319576.5625, 6322145.3125, 6322181.25, 6324040.625, 6326746.875, 6336509.375, 6342437.5, 6352259.375, 6352775.0, 6362039.0625, 6362248.4375, 6368281.25, 6369189.0625, 6371087.5, 6373585.9375, 6380932.8125, 6394240.625, 6399618.75, 6410535.9375, 6410545.3125, 6413240.625, 6416323.4375, 6418545.3125, 6428387.5, 6430470.3125, 6433681.25, 6435090.625, 6435167.1875, 6436484.375, 6444385.9375, 6447400.0, 6448682.8125, 6448700.0, 6454200.0, 6457529.6875, 6460623.4375, 6467598.4375, 6470059.375, 6471812.5, 6481778.125, 6481903.125, 6485557.8125, 6487559.375, 6489975.0, 6490014.0625, 6491114.0625, 6492940.625, 6493853.125, 6496859.375, 6498262.5, 6502225.0, 6508539.0625, 6513260.9375, 6513512.5, 6514801.5625, 6518090.625, 6524323.4375, 6525020.3125, 6525940.625, 6528212.5, 6529785.9375, 6536190.625, 6538281.25, 6539114.0625, 6540150.0, 6540151.5625, 6540615.625, 6540617.1875, 6540714.0625, 6542154.6875, 6546798.4375, 6546917.1875, 6546948.4375, 6547179.6875, 6547817.1875, 6548545.3125, 6548573.4375, 6551942.1875, 6552953.125, 6557951.5625, 6561510.9375, 6562268.75, 6562459.375, 6563432.8125, 6564185.9375, 6567400.0, 6568051.5625, 6571945.3125, 6576687.5, 6585045.3125, 6585960.9375, 6585993.75, 6586079.6875, 6591973.4375, 6593929.6875, 6594459.375, 6594867.1875, 6605275.0, 6606850.0, 6615773.4375, 6623020.3125, 6625071.875, 6625198.4375, 6628996.875, 6629696.875, 6630925.0, 6634560.9375, 6643695.3125, 6645893.75, 6665757.8125, 6666648.4375, 6671643.75, 6673315.625, 6673625.0, 6680246.875, 6686115.625, 6689448.4375, 6690026.5625, 6710073.4375, 6711153.125, 6724657.8125, 6725851.5625, 6733646.875, 6733985.9375, 6736645.3125, 6738150.0, 6738395.3125, 6738909.375, 6739279.6875, 6739962.5, 6740268.75, 6740739.0625, 6741057.8125, 6742370.3125, 6742417.1875, 6743034.375, 6744151.5625, 6762662.5, 6764714.0625, 6766734.375, 6771920.3125, 6772339.0625, 6772350.0, 6772787.5, 6778729.6875, 6780512.5, 6780789.0625, 6793123.4375, 6809100.0, 6817171.875, 6824787.5, 6826395.3125, 6830545.3125, 6837921.875, 6844735.9375, 6861098.4375, 6861168.75, 6875943.75, 6889281.25, 6891171.875, 6893664.0625, 6893681.25, 6894607.8125, 6904993.75, 6913417.1875, 6917462.5, 6934946.875, 6935445.3125, 6940257.8125, 6960084.375, 6960428.125, 6960650.0, 6962103.125, 6962539.0625, 6964318.75, 6968901.5625, 6978320.3125, 6978865.625, 6982576.5625, 6988735.9375, 6991893.75, 6995650.0, 6996182.8125, 6996939.0625, 6997507.8125, 6997509.375, 7000198.4375, 7000220.3125, 7008189.0625, 7012421.875, 7012714.0625, 7013043.75, 7013500.0, 7014400.0, 7018082.8125, 7029426.5625, 7030426.5625, 7030912.5, 7030971.875, 7031157.8125, 7031432.8125, 7047139.0625, 7047982.8125, 7047984.375, 7048273.4375, 7048384.375, 7048868.75, 7050751.5625, 7052115.625, 7058409.375, 7061468.75, 7068865.625, 7076778.125, 7092384.375, 7110717.1875, 7120039.0625, 7130389.0625, 7130400.0, 7137335.9375, 7144703.125, 7144757.8125, 7147253.125, 7154360.9375, 7155145.3125, 7165195.3125, 7176370.3125, 7193917.1875, 7197621.875, 7198029.6875, 7205201.5625, 7226781.25, 7234670.3125, 7239090.625, 7243629.6875, 7247364.0625, 7248442.1875, 7267871.875, 7273971.875, 7288228.125, 7309335.9375, 7311043.75, 7311104.6875, 7314506.25, 7315379.6875, 7316137.5, 7337059.375, 7337065.625, 7390439.0625, 7435539.0625, 7449290.625, 7450670.3125, 7463337.5, 7463343.75, 7470698.4375, 7470970.3125, 7472193.75, 7476443.75, 7478492.1875, 7478557.8125, 7510765.625, 7518442.1875, 7542168.75, 7570048.4375, 7595409.375, 7596476.5625, 7597328.125, 7625543.75, 7650909.375, 7721596.875, 7722310.9375, 7727210.9375, 7764625.0, 7787993.75, 7791193.75, 7803498.4375, 7807184.375, 7810235.9375, 7812321.875, 7817117.1875, 7829871.875, 7846398.4375, 7847503.125, 7847651.5625, 7849807.8125, 7855193.75, 7855264.0625, 7860251.5625, 7892929.6875, 7898376.5625, 7901004.6875, 7924506.25, 7940000.0, 7958260.9375, 7959453.125, 7982820.3125, 7988326.5625, 7991940.625, 7991996.875, 7993918.75, 8005654.6875, 8014804.6875, 8018310.9375, 8019214.0625, 8037101.5625, 8037159.375, 8038393.75, 8043084.375, 8044379.6875, 8048943.75, 8050307.8125, 8054737.5, 8056631.25, 8057229.6875, 8101037.5, 8104162.5, 8113107.8125, 8129426.5625, 8132382.8125, 8135295.3125, 8169570.3125, 8169612.5, 8171153.125, 8172485.9375, 8190481.25, 8243443.75, 8249378.125, 8249789.0625, 8267040.625, 8311570.3125, 8313985.9375, 8315790.625, 8335160.9375, 8336242.1875, 8336887.5, 8361595.3125, 8364807.8125, 8370550.0, 8371848.4375, 8373148.4375, 8381481.25, 8391876.5625, 8422446.875, 8424526.5625, 8427985.9375, 8428779.6875, 8433364.0625, 8438934.375, 8519918.75, 8519971.875, 8525903.125, 8726356.25, 8763387.5, 8929082.8125, 8971112.5, 9008476.5625, 9011539.0625, 9046957.8125, 9068823.4375, 9069871.875, 9074462.5, 9090970.3125, 9103814.0625, 9112545.3125, 9130581.25, 9138062.5, 9148682.8125, 9191271.875, 9191565.625, 9193410.9375, 9219437.5, 9223810.9375, 9225489.0625, 9243531.25, 9247415.625, 9247484.375, 9249042.1875, 9259990.625, 9260784.375, 9265714.0625, 9266857.8125, 9267235.9375, 9267825.0, 9324629.6875, 9325446.875, 9325659.375, 9330354.6875, 9330698.4375, 9351007.8125, 9356687.5, 9371173.4375, 9377645.3125, 9381504.6875, 9383565.625, 9384579.6875, 9409228.125, 9409454.6875, 9409468.75, 9409582.8125, 9417356.25, 9447523.4375, 9471201.5625, 9471673.4375, 9476043.75, 9477607.8125, 9477642.1875, 9490775.0, 9493270.3125, 9525168.75, 9847228.125, 9851789.0625, 9972010.9375, 9974273.4375, 9996565.625, 10016776.5625, 10032243.75, 10032246.875, 10033557.8125, 10079851.5625, 10083731.25, 10094964.0625, 10099923.4375, 10100248.4375, 10107257.8125, 10107665.625, 10108107.8125, 10108204.6875, 10108870.3125, 10109131.25, 10109329.6875, 10110521.875, 10110553.125, 10110634.375, 10110637.5, 10111701.5625, 10111823.4375, 10112857.8125, 10116873.4375, 10117326.5625, 10127357.8125, 10137668.75, 10490954.6875, 535812517.1875], [61.730579369789076, 48.223354137396974, 15.20526725270402, 13.981023110054315, 29.87760757312946, 8.8068870672802, 53.287516968949696, 20.773058536750135, 36.16696175622107, 7.278502931360082, 69.8195187626727, 32.06916357656553, 91.40144064554639, 17.306456813279414, 71.82575058541576, 17.26115223846553, 36.82005593534155, 58.28714650826787, 29.247214909458286, 16.556835375852387, 94.37830766487349, 63.40329793728067, 53.301260756556076, 28.515792196900513, 22.679851731987373, 54.11880073029788, 24.45068907703016, 15.43737716502371, 7.577281942084358, 23.264514576495277, 10.97482042204417, 103.72941902316018, 69.91365078042362, 12.015182099028792, 63.58525176070736, 32.53758904004177, 103.56127250443365, 17.25103969670652, 21.622900519852916, 65.4284915646054, 13.610255899543073, 74.07951269563867, 78.9746427990302, 95.12001068498806, 10.451888871851619, 12.906966405783828, 7.515684239308004, 13.20426808065657, 53.90467035386394, 6.34446665664832, 41.593783170951454, 59.05137608731858, 5.985211047012849, 27.22913630424391, 8.493763433298101, 13.577839784342036, 24.18640466361277, 41.139901008758045, 91.63931870205877, 43.05841521024962, 78.51352714510219, 12.596699906364691, 11.862274207652536, 126.63405666656021, 9.94438284457542, 5.6382213073732546, 21.3045729749904, 7.03082262538515, 55.699376512973146, 21.134967891258377, 64.71416432393306, 115.73982724103377, 6.180512024523153, 54.19827351758506, 5.519729879265535, 39.58878192973954, 93.13651725991707, 47.32969561538789, 27.189196783452815, 31.387098532659103, 59.020388054702046, 79.2507291931169, 47.066529697195385, 143.9945885186026, 8.725691715282995, 65.70089970795175, 8.368557040809657, 38.2257745238806, 109.04133013783127, 56.61446660663913, 14.260438451569263, 17.885764170992623, 9.36848833732462, 11.759256207683247, 51.54628884884034, 15.839446681348281, 14.648683500482521, 81.51910588398297, 84.37002493283056, 9.245702555999216, 20.89688449327884, 6.171814498927349, 9.599322371419133, 10.40409590451822, 58.574494745199154, 5.030764341000801, 68.47480835064042, 15.000266253040886, 53.68535395660619, 14.800630748177747, 41.01064346390145, 9.12712243076058, 6.839159486730435, 26.842019915616735, 42.20707261622759, 35.93148929941992, 155.18396593260155, 10.60221257145808, 12.169251424823665, 13.410961640024007, 31.189443956598364, 47.304515957635275, 16.95262957251609, 31.73403784197447, 29.886110775353764, 66.07933477518047, 9.458399817358995, 23.655391710780297, 39.712076609853554, 20.806342781595227, 22.44993739040573, 43.285001593626326, 17.874312048511108, 9.468639434257565, 77.21973196699481, 10.42224688025424, 12.167689740481698, 48.56171347889771, 32.131744712693646, 5.128008457959243, 5.172687735865324, 62.01907078642327, 10.897371785278276, 42.28373698904865, 21.797581674781497, 9.445183980373956, 88.44222765667955, 7.239103518220504, 66.28352958259907, 83.54178917247957, 41.44306119551484, 81.89117104918886, 69.77496186655, 18.563519036306776, 9.191454933722552, 45.00750389329174, 88.01979931049114, 12.507170616932315, 18.234274268110326, 5.992280817301156, 10.829434963609172, 40.87166846994279, 32.95465470852147, 61.47347200341765, 9.153427359979961, 106.66700754916832, 53.13847538353467, 5.238805293416969, 8.23966495765203, 71.76374562771295, 5.113119931181222, 85.30577594787123, 46.95979836680726, 10.101147287581615, 100.06790626236443, 103.62546590057887, 5.288061662575674, 70.66371335672221, 88.92104396831037, 27.55024671969938, 17.675900795679045, 24.99696370170548, 11.807711105911824, 63.02792756097283, 7.925805942846159, 129.73757803980965, 115.01321023545769, 35.93698464243087, 19.147855659516466, 32.471437573622886, 12.360062693130613, 71.60958832773719, 46.863618105184614, 79.35380088287634, 28.483111377874053, 19.731372525701865, 13.338696102599023, 15.583814341811001, 44.84439561526891, 72.48995876389458, 218.95124201606626, 88.29097742282828, 110.02368492764231, 68.71600603673434, 59.05321111590324, 18.255761081425028, 8.273092663150978, 6.444423622082854, 11.968570930337682, 7.724159656255928, 27.841057381835878, 104.33712467205673, 11.462106339133992, 12.84168828848509, 65.37598566482066, 71.89401455411569, 5.414974034604237, 21.98536502795923, 82.22866062989874, 52.64603610257542, 21.338303581956367, 26.64318439502845, 77.04417609358892, 11.794104617377744, 109.53248715277338, 23.15473535517404, 7.420073722467204, 14.87005110966682, 5.244533443558159, 53.84804212457197, 10.91531745084978, 64.8133654698495, 31.565172727717496, 5.350243591334018, 19.988785946398895, 24.503686371718356, 9.716804283948901, 10.730386066795244, 27.459666727968482, 13.205529583902004, 18.418677884163525, 72.50527529739027, 30.365570785184836, 75.77208667303097, 24.075399494717622, 102.1657487479346, 13.802320103881083, 77.90057282903221, 16.838557146261092, 14.393087009053406, 11.6074338682422, 6.0718285224190325, 8.56149839962373, 7.329834011664996, 41.09463026792563, 121.48962804856919, 85.25947126929407, 60.08825528347694, 45.61652574287363, 50.84051923943387, 36.91125650911725, 32.11498207648948, 83.22881082285933, 14.517755401762706, 51.11155282468238, 25.330047256903022, 48.85817502138973, 45.68152032082069, 30.181022107928314, 72.53806656229108, 24.290739955686664, 59.2228569346157, 15.157472683568715, 24.412246767777454, 8.708342343333262, 101.05161961808629, 12.6818279486582, 39.08604279829844, 21.246873279726923, 13.028283091784859, 58.094650754922924, 21.337914665854708, 128.93703562507434, 24.23800095610355, 9.118804817790405, 6.047818101148177, 26.084261441212526, 45.333068429167, 17.76154547381708, 69.45959712216444, 63.98430500054071, 36.31862808164117, 28.9640825109195, 42.22738900890218, 41.18681810685485, 56.30445527598968, 24.23015793237512, 8.580936491118582, 10.710011393072676, 31.319274681951722, 36.740677176641576, 79.22652428895694, 83.85801107450524, 123.02659888141588, 25.45291106510391, 6.6003120484145885, 25.37959993260264, 62.296225374618864, 18.65264706888196, 42.86541123427597, 56.23518116085906, 77.91505951911773, 35.61965400930281, 12.986764137230292, 16.61274538307657, 7.682075335730902, 127.9640678428251, 28.936709998234544, 34.342105849768764, 30.15412218623032, 65.09603266959236, 15.688650737004744, 126.58777187957705, 9.0005751614303, 55.93454020344289, 34.73161008094415, 5.04338359456866, 11.715795681911231, 157.82052711413485, 19.504252519933612, 19.972752961652617, 11.723023883172063, 128.7009087298802, 51.44214591967797, 12.854145827417124, 90.00071724468377, 18.20750916385242, 42.302904960299955, 53.78025738594092, 9.140328519765541, 19.34818832214175, 39.9207867897276, 39.62879582296684, 37.60952911487856, 126.83127582860953, 8.33835221239202, 22.431608624716453, 21.086272045742575, 22.746047707759455, 40.491007644839705, 5.430735863960645, 101.00769044521091, 72.10993127852295, 22.75785109615057, 28.220112517481425, 14.327136913342457, 5.896760635786537, 9.271446841842039, 30.8769158144217, 7.153191961043938, 39.861171752483294, 15.11653859729901, 64.74602381584143, 27.01840532883351, 8.726637264477676, 56.952328243439396, 17.56655036647452, 12.236078558099925, 5.952318647711503, 5.534676756156285, 14.752929394495096, 15.372911106930875, 25.14682043715115, 12.31346510394708, 45.474799318379695, 28.102158067394786, 90.00540719219475, 21.658676528384206, 72.072381201717, 30.487563945580888, 24.555115297861207, 23.014581948932403, 37.98698636976693, 54.64118323186963, 107.0656874706244, 86.59438558268022, 30.829890313231704, 63.09626140804717, 7.8333189604034485, 10.346840813722324, 22.179946091522975, 13.615936569082141, 37.047275528887056, 19.446964474825027, 8.858553387425975, 11.439775069551928, 7.880245870814844, 15.578015561114382, 9.361808580993966, 19.767045034667653, 20.454542637760525, 76.58506325151703, 74.37606273264524, 10.225619584479105, 106.64777289345308, 7.999693516985066, 61.221941295517894, 36.07623861258619, 56.763285252860356, 25.196490665040216, 20.8713167473563, 36.10083995059748, 57.839942705430914, 20.60996143235243, 43.55151332708882, 23.0945368068699, 63.70677679557721, 8.290141509220348, 44.399001186510695, 41.881703554806734, 35.93480651895555, 39.79009598533669, 42.573157210624125, 26.632371319994455, 109.88974184229326, 11.024083676058078, 117.82100047708973, 154.12989079672278, 24.426219759491794, 152.7302157002492, 19.316585517140265, 6.8391741773212855, 8.445286318466541, 11.494699619130143, 63.40773352512696, 77.53625675997054, 16.722336522321818, 67.4793147193011, 114.94476157414326, 26.338395498068074, 7.622782401857462, 73.7995709200581, 65.88101933712124, 16.870487486741073, 92.83743446798024, 229.69474019658446, 39.29683955621515, 22.012955646277167, 28.10693353945151, 113.30646744647655, 164.09290827294507, 11.433611625337837, 19.390946013142216, 16.81448483606927, 5.082604255421736, 15.397183127648422, 33.63220620396156, 26.876544950873253, 68.28937473644041, 63.745136354713814, 28.634646864644168, 46.85764740630577, 52.830909040636044, 26.858852487381952, 66.11557651941142, 60.96823927614043, 59.85881440455127, 6.1570238904939565, 5.079603237050903, 6.572933635204267, 141.3910471895039, 17.833839867158765, 5.733418509969664, 33.56693838215604, 132.48797259861647, 19.801668201657108, 30.4397560883312, 12.84805224530508, 82.05945727395628, 28.35071327818005, 34.01331663929978, 29.320477133011078, 24.42253488219739, 155.64623068083208, 132.53883722411933, 30.295453693633927, 19.74057419197362, 8.792385057924562, 35.72537290810086, 12.594218922151232, 6.259744927709634, 16.602841157629637, 5.032665288337748, 34.99186781037157, 119.0094341636794, 6.399587376466778, 21.224829518155637, 64.48197170158335, 11.61595974111695, 14.56866216750861, 122.39595888332506, 17.995454509535463, 45.49510380443273, 67.59272937590914, 23.108238644303803, 6.1382265446787665, 53.79902104229021, 72.61577948932128, 52.45230115471795, 5.5055059120451375, 41.97164469007612, 27.51582335218675, 96.46941663923958, 88.3883180024952, 36.15140748755064, 185.5076619031931, 110.95142861150724, 76.2890693392739, 76.09063633221376, 5.166983389575908, 7.359926138544784, 55.164228081777544, 36.39429755335558, 48.90130337246897, 131.55657843021274, 21.733859869742002, 20.338290749005065, 30.07421291123829, 57.09408547942293, 5.701676502283257, 60.135476252272916, 32.491046287927, 54.05816027912783, 11.323303579211828, 49.32710295752727, 49.51388961188801, 112.8132908393257, 22.910951138450866, 5.070490242199125, 26.900182701594943, 53.73237526311319, 107.85546321330995, 23.70808576366379, 12.595883924333544, 35.67117078210234, 7.0640448999035454, 33.5982876391179, 24.836732500846797, 21.088197932078568, 30.260062375058325, 47.7631293809177, 66.08738590276428, 58.731198018757524, 146.00597035953444, 8.25580950200558, 70.34729137355163, 28.497720448015038, 36.63226136615714, 31.222535988338773, 5.665674900894059, 19.90198655424581, 15.856665956534949, 19.47795910997626, 10.921485093674768, 23.40790027024161, 5.355158333292589, 8.022416314295864, 87.3450577773408, 57.67415011156858, 13.80271535250499, 50.169500874922186, 11.041904057194687, 21.230512664101596, 151.44941274291034, 57.71101445413079, 16.417042152107918, 38.09155941628943, 55.6548870462807, 9.416910045766192, 5.0531352479069715, 39.45245823581494, 92.6197017735598, 17.032492279860367, 48.513704753347426, 13.138012817709667, 16.21513101353392, 43.686168398547125, 64.34531165015599, 36.24362298285726, 7.115637965153101, 30.700536713015218, 158.45750899636462, 90.39756825912099, 87.6456846868272, 6.175872922082695, 25.660783736871988, 13.818208392420315, 83.45359045919835, 42.33272960116086, 71.73893665242136, 43.11459858431879, 148.47960480981152, 5.59308165567319, 133.72717279455182, 61.369871928484116, 42.16265595010627, 19.46270464793236, 7.937189834498064, 130.78433580343727, 46.10557380232831, 8.023480231312336, 60.546793549042384, 14.09444739040739, 88.38143796926721, 11.602622602592488, 27.07879043195964, 11.822160457638445, 35.597386885408355, 123.65749981067931, 31.572016262239657, 60.97174393740824, 16.169729212878526, 25.029493975786217, 78.09176721354613, 12.20208574313815, 75.14029509746577, 30.16579318329695, 16.87403508488182, 17.35639238534709, 19.363942356541926, 6.797366901439606, 10.536594218599937, 91.40471737348084, 14.578262980450925, 45.22101724713076, 136.21667749142208, 6.401599496395278, 19.564644025293237, 57.00573166859808, 162.70103716014182, 84.0528027076331, 49.85870402683705, 25.14827671244069, 33.399246446039825, 36.813844454554356, 53.887075476995136, 6.016720100034734, 8.399397618488644, 44.01538259014254, 52.12184027879909, 64.24457395444753, 17.742198116445476, 21.505004734043997, 22.992561906928316, 56.76952034428838, 5.86958070003278, 13.661911845359016, 27.285096729424538, 15.53904600872319, 53.21185964708374, 13.373586004609324, 49.84541988889123, 8.233810708136899, 27.847489858416072, 11.845559703041914, 12.99774655888604, 58.151546800382306, 64.71085777608471, 24.892115376477363, 83.97508281992219, 15.76647395019558, 11.813324859478813, 13.638367506235745, 5.23169709021714, 57.73523655604864, 107.30746304912473])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)