In [25]:
#hit map
hit_map(detector,hit_map_fig,RowNo,ClmNo)
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 = 44409
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);
([3326887.5, 3336750.0, 3430026.5625, 3457040.625, 3465812.5, 3468035.9375, 3469092.1875, 3470846.875, 3471004.6875, 3471117.1875, 3472431.25, 3473471.875, 3476746.875, 3477715.625, 3489785.9375, 3502885.9375, 3505543.75, 3512221.875, 3514873.4375, 3533598.4375, 3535301.5625, 3536723.4375, 3562912.5, 3577376.5625, 3580067.1875, 3580279.6875, 3586718.75, 3587570.3125, 3595275.0, 3602159.375, 3612182.8125, 3629662.5, 3637182.8125, 3637484.375, 3640037.5, 3642481.25, 3645756.25, 3666034.375, 3668587.5, 3669707.8125, 3673342.1875, 3677042.1875, 3686990.625, 3699889.0625, 3712414.0625, 3713782.8125, 3732732.8125, 3735007.8125, 3737292.1875, 3740648.4375, 3741112.5, 3742118.75, 3743825.0, 3744315.625, 3748340.625, 3751559.375, 3752273.4375, 3755212.5, 3755807.8125, 3758048.4375, 3767800.0, 3775487.5, 3775782.8125, 3777957.8125, 3778107.8125, 3779579.6875, 3781940.625, 3782976.5625, 3783531.25, 3785453.125, 3785595.3125, 3791892.1875, 3792293.75, 3792343.75, 3794779.6875, 3794807.8125, 3794885.9375, 3795123.4375, 3795201.5625, 3795425.0, 3796775.0, 3797351.5625, 3797432.8125, 3797748.4375, 3797817.1875, 3798037.5, 3798153.125, 3798317.1875, 3798760.9375, 3799110.9375, 3799131.25, 3799132.8125, 3799193.75, 3799443.75, 3800071.875, 3800340.625, 3800468.75, 3801090.625, 3801528.125, 3801806.25, 3802296.875, 3802773.4375, 3802925.0, 3803112.5, 3803467.1875, 3804554.6875, 3805026.5625, 3805328.125, 3806145.3125, 3806717.1875, 3826245.3125, 3833881.25, 3837154.6875, 3855659.375, 3862954.6875, 3890884.375, 3895787.5, 3906654.6875, 3913764.0625, 3949573.4375, 4017015.625, 4022571.875, 4030139.0625, 4057093.75, 4068350.0, 4080389.0625, 4081840.625, 4128017.1875, 4148826.5625, 4148854.6875, 4153829.6875, 4161646.875, 4166392.1875, 4180204.6875, 4191534.375, 4193710.9375, 4194979.6875, 4196157.8125, 4196228.125, 4199753.125, 4220971.875, 4222853.125, 4223217.1875, 4239856.25, 4258600.0, 4262967.1875, 4267545.3125, 4288862.5, 4294617.1875, 4294910.9375, 4303500.0, 4305481.25, 4305985.9375, 4308029.6875, 4310221.875, 4310242.1875, 4318239.0625, 4324614.0625, 4327668.75, 4331854.6875, 4337943.75, 4340582.8125, 4347409.375, 4367134.375, 4381517.1875, 4383151.5625, 4387140.625, 4398503.125, 4400676.5625, 4404487.5, 4404500.0, 4427906.25, 4433817.1875, 4435356.25, 4435442.1875, 4440901.5625, 4458620.3125, 4459748.4375, 4463070.3125, 4463907.8125, 4465964.0625, 4465967.1875, 4469293.75, 4473443.75, 4490329.6875, 4497471.875, 4500485.9375, 4501984.375, 4508179.6875, 4511987.5, 4514610.9375, 4514971.875, 4522603.125, 4529471.875, 4530121.875, 4530643.75, 4531968.75, 4533301.5625, 4536485.9375, 4538485.9375, 4544220.3125, 4544756.25, 4545223.4375, 4545935.9375, 4547050.0, 4547217.1875, 4548378.125, 4549014.0625, 4549131.25, 4550893.75, 4554931.25, 4563829.6875, 4569568.75, 4572073.4375, 4572143.75, 4572332.8125, 4574929.6875, 4574951.5625, 4575776.5625, 4579729.6875, 4582498.4375, 4585431.25, 4595532.8125, 4596210.9375, 4600998.4375, 4608214.0625, 4612792.1875, 4640540.625, 4646134.375, 4654810.9375, 4667539.0625, 4668748.4375, 4670723.4375, 4676304.6875, 4701053.125, 4708215.625, 4711590.625, 4714379.6875, 4715809.375, 4722696.875, 4723876.5625, 4724245.3125, 4724629.6875, 4724970.3125, 4725031.25, 4731400.0, 4732053.125, 4733025.0, 4733084.375, 4735140.625, 4740709.375, 4741070.3125, 4743775.0, 4747467.1875, 4760954.6875, 4760976.5625, 4770456.25, 4783200.0, 4790409.375, 4793589.0625, 4799028.125, 4814804.6875, 4820009.375, 4821987.5, 4822035.9375, 4823054.6875, 4826410.9375, 4827120.3125, 4827203.125, 4830959.375, 4830973.4375, 4832904.6875, 4833489.0625, 4844934.375, 4845545.3125, 4845893.75, 4849990.625, 4853407.8125, 4855856.25, 4858214.0625, 4863420.3125, 4865295.3125, 4866053.125, 4870170.3125, 4874395.3125, 4877328.125, 4880217.1875, 4880818.75, 4896176.5625, 4896715.625, 4902978.125, 4911889.0625, 4915531.25, 4917796.875, 4923159.375, 4928773.4375, 4935479.6875, 4943054.6875, 4945256.25, 4945687.5, 4947821.875, 4947912.5, 4952792.1875, 4953240.625, 4953317.1875, 4953471.875, 4955229.6875, 4958467.1875, 4958510.9375, 4964034.375, 4964064.0625, 4964178.125, 4969896.875, 4983842.1875, 4990490.625, 4992773.4375, 5001396.875, 5003035.9375, 5006504.6875, 5008365.625, 5008385.9375, 5015668.75, 5023470.3125, 5025200.0, 5028593.75, 5030446.875, 5032739.0625, 5045187.5, 5050703.125, 5054545.3125, 5054554.6875, 5054651.5625, 5057193.75, 5059929.6875, 5071673.4375, 5077125.0, 5080265.625, 5089387.5, 5094989.0625, 5107109.375, 5128735.9375, 5145960.9375, 5146853.125, 5153967.1875, 5158428.125, 5170628.125, 5197831.25, 5208023.4375, 5225720.3125, 5227214.0625, 5240971.875, 5267218.75, 5272814.0625, 5273248.4375, 5273268.75, 5398412.5, 5415448.4375, 5435407.8125, 5446914.0625, 5450184.375, 5453114.0625, 5454279.6875, 5458439.0625, 5459515.625, 5466901.5625, 5470607.8125, 5478843.75, 5478865.625, 5480057.8125, 5482807.8125, 5499718.75, 5515796.875, 5530859.375, 5531109.375, 5535317.1875, 5545223.4375, 5550046.875, 5554468.75, 5556876.5625, 5556904.6875, 5559548.4375, 5580939.0625, 5581068.75, 5582159.375, 5583506.25, 5585168.75, 5586623.4375, 5610334.375, 5617987.5, 5623110.9375, 5625445.3125, 5627564.0625, 5647457.8125, 5651701.5625, 5651725.0, 5654246.875, 5668335.9375, 5672264.0625, 5680914.0625, 5682459.375, 5689290.625, 5690996.875, 5706857.8125, 5707810.9375, 5721684.375, 5742173.4375, 5781223.4375, 5781284.375, 5799792.1875, 5805968.75, 5816207.8125, 5889990.625, 5895464.0625, 5941703.125, 5979157.8125, 5990748.4375, 6019798.4375, 6033767.1875, 6060998.4375, 6064048.4375, 6101187.5, 6110467.1875, 6130032.8125, 6132146.875, 6158523.4375, 6176171.875, 6211081.25, 6251334.375, 6269348.4375, 6277125.0, 6293648.4375, 6297175.0, 6300803.125, 6311553.125, 6321376.5625, 6337212.5, 6359751.5625, 6359839.0625, 6366684.375, 6368276.5625, 6392606.25, 6404895.3125, 6408560.9375, 6412671.875, 6412967.1875, 6432392.1875, 6434925.0, 6449143.75, 6453175.0, 6453189.0625, 6457040.625, 6464821.875, 6465721.875, 6467864.0625, 6469910.9375, 6481212.5, 6481223.4375, 6484256.25, 6490796.875, 6491153.125, 6492528.125, 6494725.0, 6498693.75, 6501935.9375, 6501946.875, 6504720.3125, 6512354.6875, 6520545.3125, 6520851.5625, 6522270.3125, 6523392.1875, 6527471.875, 6534904.6875, 6552546.875, 6557964.0625, 6570459.375, 6572651.5625, 6577250.0, 6584509.375, 6601348.4375, 6646229.6875, 6669004.6875, 6671432.8125, 6682762.5, 6700582.8125, 6702760.9375, 6746168.75, 6772123.4375, 6773493.75, 6778526.5625, 6781509.375, 6786373.4375, 6790454.6875, 6798782.8125, 6800720.3125, 6803750.0, 6809440.625, 6809489.0625, 6832151.5625, 6832212.5, 6878804.6875, 6880190.625, 6988062.5, 6989648.4375, 6991495.3125, 6997343.75, 7062939.0625, 7088926.5625, 7443489.0625, 7454237.5, 8496096.875, 8612093.75], [50.944314495704695, 15.796553761775797, 39.599775479332536, 8.137706685417054, 63.074565691841194, 80.21789700366425, 15.960937936626587, 12.512999259053608, 29.837789083558576, 53.120041011142696, 27.89621550330027, 47.63394783322977, 18.300284868169413, 50.57820505054154, 35.690030261814336, 67.36186904267599, 7.196413843423989, 65.38448733491555, 24.998896924739103, 6.509668300256832, 19.315350304061933, 5.439443571260751, 25.632901477304635, 66.16595145584803, 19.278302018446666, 14.15874738738331, 6.522623815425454, 107.53968102319064, 33.41513750724508, 5.316264066312792, 46.10586863436758, 91.16589071864006, 48.09769939310871, 7.158275626253481, 54.74312561294476, 6.794411236218008, 13.644069024890666, 9.542268560378988, 24.600628551372115, 8.980340001393829, 10.803338798417208, 40.45007002599705, 68.50365582357375, 16.959660412291317, 21.905485907458093, 8.268913915283843, 62.64292542871192, 91.47650348548203, 58.34117945242116, 25.060641910279976, 37.35763724901154, 40.520881028062504, 26.911704247992077, 5.121250559007611, 23.60418085712736, 5.700602386763834, 17.9621753510238, 59.68005251310863, 69.0006913012932, 10.955843315049043, 25.84869330680094, 94.43485138712042, 33.13091245400119, 23.792743755176126, 73.61615820046723, 45.69454257566443, 8.750624063703121, 71.5191867657087, 64.94882434073011, 5.711309894727182, 20.481644961295867, 38.950150823930734, 18.088430464499684, 5.647520900683906, 12.135772125361198, 13.752119500296295, 13.091904364232594, 27.168080232946714, 113.1331687422565, 12.86995546178879, 12.14858966370281, 44.947831610102234, 43.854046163781256, 78.75315976097336, 93.75752775515157, 52.20668232804204, 55.90487138715627, 16.09285818625129, 19.798294202270533, 60.54741029516812, 17.389647197718674, 15.9963206441607, 54.395748676954206, 23.103191481878326, 51.80352019490101, 59.57832523208183, 11.51942720465644, 28.086517632765965, 6.065950813993589, 10.13062257170731, 10.574647561438637, 21.320669548016888, 19.745649394005632, 61.228638772155804, 72.73013554124239, 5.2719012231358695, 32.57270000951705, 7.253053753444399, 7.280814963669682, 42.912864000026964, 67.61316065124427, 5.62951948229284, 7.085303980947569, 59.43974368411783, 17.31304622136468, 64.10192608274598, 11.492817322970527, 33.20667022877087, 15.024348292525623, 9.834152479081345, 9.56967415166326, 64.06353291501641, 79.69358726818837, 7.549817268686985, 73.97081699540749, 17.74799126293028, 49.266121924875584, 17.13425397130002, 27.214492930305425, 12.195626065470634, 94.37692473064298, 16.373007131649477, 42.14259295651776, 7.98458684063542, 75.55815453060413, 35.12623972643939, 49.602479270763496, 17.224787612594323, 5.169139747263604, 26.830700834371942, 25.67510155413373, 27.48320485520895, 44.32946137990008, 75.48996115018434, 58.447575040326555, 36.518753572604936, 5.521493061442035, 8.031471558821353, 27.821059716650893, 13.163987644669419, 25.786962165195888, 41.54410680529542, 21.374984613534597, 34.703267100647466, 88.53354610613425, 12.17820705532583, 45.08392716131503, 13.34068538552218, 10.626892004130326, 40.55367665331279, 5.029691580161599, 89.72498654448499, 26.63211685324429, 12.71011993697828, 17.96824894569834, 37.943651597656164, 70.22964903144896, 73.1276765720511, 106.92512423577564, 20.20393569463335, 15.249400727842824, 10.410655364025292, 50.12537437218167, 161.61939059722428, 5.21016275423511, 10.861592297274491, 13.98605953305139, 25.778956749256366, 14.351774903446644, 19.987153139850143, 22.08495909861419, 28.94689530639569, 73.53783070486082, 83.43072176960183, 78.51153103982357, 31.492972990266875, 5.470208713948725, 16.366000167455127, 33.47860789500325, 71.33295277617917, 15.776243049260728, 22.907238928496348, 83.02688974523193, 26.25698520129269, 5.786110069010675, 36.033767298774734, 8.339859353035278, 13.805229401299217, 6.588777894823672, 31.577525446498296, 72.44380337767733, 9.46734512134162, 50.183104440415796, 6.366286970021599, 23.376396962840523, 62.504425608304146, 19.384430652404177, 37.407907875019774, 16.118298272153197, 20.34051371970572, 6.8365284522596435, 11.223147344283218, 9.508576428819588, 78.63040007514823, 10.766912906583594, 61.598593402889236, 26.181738301365808, 78.04871585780177, 89.40377328229651, 40.29039994151957, 129.79098658785202, 63.755752766138386, 86.46792473340524, 125.78626764771685, 76.2125002110638, 25.498601563515233, 28.312596114138803, 30.941215692162736, 9.054369042689874, 14.164660251426868, 18.693462145242975, 19.899878014321146, 43.26646387948695, 27.972438283688014, 29.18301157608849, 54.57432964556627, 38.62558218029787, 65.49394062729058, 37.691739056808075, 18.710073706045748, 37.102507653289834, 99.38100862261118, 103.78187255075929, 105.32317128029207, 6.067994961857043, 121.59654713008884, 25.010919430937037, 13.252717642964797, 5.401295107188826, 6.249232438686556, 54.48134949199029, 9.790394190622914, 103.45283563733504, 37.51866728377649, 14.745446678476293, 59.74435588506283, 14.983347717471712, 9.514996551443158, 41.802267844644426, 19.66102646435677, 16.46179074001237, 66.92915642180607, 128.02455324041892, 72.49916501972811, 68.96663773029925, 72.90284714656786, 26.720954706378436, 11.33293797430555, 36.34627591652319, 60.40303690511706, 29.9650978750554, 117.99923337949089, 17.094445723828184, 91.70197062341235, 44.722346253657754, 36.51871247738389, 61.322620210664596, 5.438962070653053, 38.671147273977304, 97.95481064266868, 33.55581502171796, 65.953701445872, 6.723190745924475, 73.77858508677758, 15.71509270840959, 15.355024075580742, 62.652882573216765, 36.625712796389536, 6.672648004620804, 58.04449545218593, 15.370729463586787, 8.467784766132016, 100.66027502822826, 32.041961162279236, 17.438030165431037, 130.96867778309647, 41.12916289466618, 52.23345039876758, 123.86519551561767, 154.35490746532048, 11.658505625421697, 25.648847646213042, 71.56051894896395, 20.18543790349799, 77.60325287608424, 49.88799803276582, 10.627805648995272, 102.72709699499185, 18.495906855980817, 7.939345881218487, 10.741123252350581, 42.85624952299829, 34.96980933190555, 44.087608556740314, 42.82915315862963, 37.48159910805486, 7.110450032533722, 19.765006074314748, 20.785542498566095, 60.679918771112995, 106.13644172889829, 55.999691995288416, 31.898404678909362, 89.3954371006658, 55.01484533186828, 41.01534224610968, 23.07885957068808, 11.384386466817384, 8.50711056549528, 23.983318871884126, 15.096644262654838, 5.1155136764734035, 18.736755183604288, 7.855463055278699, 30.66548143806213, 38.60271901350903, 14.957982285149768, 27.11355877172525, 38.69326356259688, 83.38412336903227, 32.73809522327452, 14.112279019704951, 10.725791158170066, 47.124592597325226, 32.76567063421831, 114.23063205732157, 34.81246903405556, 53.16611174253132, 32.47627953177644, 33.81104503044368, 7.137695223752519, 27.84019008140716, 30.715824413239577, 61.91498337181847, 49.61342053035381, 13.640448684487664, 27.102532530390054, 20.126295418323753, 36.92291451479147, 11.564345437419183, 27.340257804324402, 34.38113355977191, 8.275789781157172, 52.96589467592817, 30.413434174258548, 63.35069299412984, 89.93052660655312, 11.097756967341901, 14.182264921613376, 25.582335321899002, 27.156526344122497, 41.35811313786884, 43.90833221232265, 16.230177861795134, 47.81085930995156, 6.876853992204028, 31.672761754490352, 5.955097910899297, 166.85813640451374, 11.411397247443487, 44.70645638228816, 18.07786856143491, 5.185435539395201, 65.95546909937444, 9.971028261650556, 118.39103987322837, 49.890096393814595, 69.41531786768181, 26.21951758573789, 21.23867454555524, 10.143994802986827, 32.51560401289473, 21.79621245669433, 48.00168823356475, 24.296253328550815, 74.27095379095478, 16.575567528892083, 73.627128048453, 85.16554642018879, 20.307248070382904, 210.8618319081322, 8.035702849446732, 12.086441242558092, 53.98103473962383, 39.99103512153062, 103.5195688754545, 17.35749940991072, 5.31059899485489, 7.2174439282374365, 23.5509437563909, 20.837841185910676, 125.89647693372366, 20.757245732220515, 16.995482224641496, 75.2894909688611, 76.38747288304455, 32.4561234393083, 165.50694487951301, 37.91771556902011, 75.95210886056685, 94.47482444903409, 37.51500454799752, 24.919236587801556, 23.86275841033838, 13.079024537964218, 66.35252802786141, 8.658965112661603, 74.07376097052445, 10.066941470899478, 14.029687147044973, 59.70204416174195, 12.460485992299345, 7.613495581663946, 10.829249025980673, 5.041197028084741, 8.570950591301235, 35.99194379275371, 5.26242857463036, 36.665630139726744, 117.88856411398243, 15.177450688591072, 18.95123630595777, 52.34620010855623, 80.47857381717395, 24.345119376864314, 53.46429305952985, 32.82346072044773, 20.55074504902428, 14.062462460633181, 7.690274984551698, 15.708740461407524, 122.16682973358604, 7.976266346969618, 19.726846351661077, 8.043416970914338, 85.00391104918818, 52.21749771569018, 131.72901701947387, 85.97567327941024, 47.05270627025925, 63.407819479049095, 32.019609506160634, 38.73526774545978, 47.40688913661187, 9.179524584071384, 13.146208003173333, 123.41685647243351, 59.04266508708121, 15.513365620683317, 5.974535652039651, 43.15826456774729, 5.803372783228023, 65.90218942643276, 19.813920478810942, 14.049400316570729, 103.34250304794251, 54.05178421181283, 51.1319095357435, 55.80989774829781, 100.2944279675149, 49.86491054820108, 12.781548653137662, 56.052358061165634, 8.221534051912824, 78.5618487473546, 39.09805077074164, 52.5747485104642, 12.250288110175825, 70.33840919877106, 10.984700483099775, 16.793637521933956, 28.57460566821389, 41.07797487023091, 55.059771295601095, 19.462488708631753, 5.593877153424963, 14.954774831355817, 6.033095841148297, 30.259375928780567, 5.48931152161992, 54.34734269218338, 33.779584732135454, 66.9756497119167, 18.67147411604182, 22.764441232765726, 5.198432572722792, 5.075672511810246, 69.42741829288413, 58.49251865541673, 10.939096943049515, 83.81919927495807, 33.220221302536594])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3326887.5, 3336750.0, 3430026.5625, 3457040.625, 3465812.5, 3468035.9375, 3469092.1875, 3470846.875, 3471004.6875, 3471117.1875, 3472431.25, 3473471.875, 3476746.875, 3477715.625, 3489785.9375, 3502885.9375, 3505543.75, 3512221.875, 3514873.4375, 3533598.4375, 3535301.5625, 3536723.4375, 3562912.5, 3577376.5625, 3580067.1875, 3580279.6875, 3586718.75, 3587570.3125, 3595275.0, 3602159.375, 3612182.8125, 3629662.5, 3637182.8125, 3637484.375, 3640037.5, 3642481.25, 3645756.25, 3666034.375, 3668587.5, 3669707.8125, 3673342.1875, 3677042.1875, 3686990.625, 3699889.0625, 3712414.0625, 3713782.8125, 3732732.8125, 3735007.8125, 3737292.1875, 3740648.4375, 3741112.5, 3742118.75, 3743825.0, 3744315.625, 3748340.625, 3751559.375, 3752273.4375, 3755212.5, 3755807.8125, 3758048.4375, 3767800.0, 3775487.5, 3775782.8125, 3777957.8125, 3778107.8125, 3779579.6875, 3781940.625, 3782976.5625, 3783531.25, 3785453.125, 3785595.3125, 3791892.1875, 3792293.75, 3792343.75, 3794779.6875, 3794807.8125, 3794885.9375, 3795123.4375, 3795201.5625, 3795425.0, 3796775.0, 3797351.5625, 3797432.8125, 3797748.4375, 3797817.1875, 3798037.5, 3798153.125, 3798317.1875, 3798760.9375, 3799110.9375, 3799131.25, 3799132.8125, 3799193.75, 3799443.75, 3800071.875, 3800340.625, 3800468.75, 3801090.625, 3801528.125, 3801806.25, 3802296.875, 3802773.4375, 3802925.0, 3803112.5, 3803467.1875, 3804554.6875, 3805026.5625, 3805328.125, 3806145.3125, 3806717.1875, 3826245.3125, 3833881.25, 3837154.6875, 3855659.375, 3862954.6875, 3890884.375, 3895787.5, 3906654.6875, 3913764.0625, 3949573.4375, 4017015.625, 4022571.875, 4030139.0625, 4057093.75, 4068350.0, 4080389.0625, 4081840.625, 4128017.1875, 4148826.5625, 4148854.6875, 4153829.6875, 4161646.875, 4166392.1875, 4180204.6875, 4191534.375, 4193710.9375, 4194979.6875, 4196157.8125, 4196228.125, 4199753.125, 4220971.875, 4222853.125, 4223217.1875, 4239856.25, 4258600.0, 4262967.1875, 4267545.3125, 4288862.5, 4294617.1875, 4294910.9375, 4303500.0, 4305481.25, 4305985.9375, 4308029.6875, 4310221.875, 4310242.1875, 4318239.0625, 4324614.0625, 4327668.75, 4331854.6875, 4337943.75, 4340582.8125, 4347409.375, 4367134.375, 4381517.1875, 4383151.5625, 4387140.625, 4398503.125, 4400676.5625, 4404487.5, 4404500.0, 4427906.25, 4433817.1875, 4435356.25, 4435442.1875, 4440901.5625, 4458620.3125, 4459748.4375, 4463070.3125, 4463907.8125, 4465964.0625, 4465967.1875, 4469293.75, 4473443.75, 4490329.6875, 4497471.875, 4500485.9375, 4501984.375, 4508179.6875, 4511987.5, 4514610.9375, 4514971.875, 4522603.125, 4529471.875, 4530121.875, 4530643.75, 4531968.75, 4533301.5625, 4536485.9375, 4538485.9375, 4544220.3125, 4544756.25, 4545223.4375, 4545935.9375, 4547050.0, 4547217.1875, 4548378.125, 4549014.0625, 4549131.25, 4550893.75, 4554931.25, 4563829.6875, 4569568.75, 4572073.4375, 4572143.75, 4572332.8125, 4574929.6875, 4574951.5625, 4575776.5625, 4579729.6875, 4582498.4375, 4585431.25, 4595532.8125, 4596210.9375, 4600998.4375, 4608214.0625, 4612792.1875, 4640540.625, 4646134.375, 4654810.9375, 4667539.0625, 4668748.4375, 4670723.4375, 4676304.6875, 4701053.125, 4708215.625, 4711590.625, 4714379.6875, 4715809.375, 4722696.875, 4723876.5625, 4724245.3125, 4724629.6875, 4724970.3125, 4725031.25, 4731400.0, 4732053.125, 4733025.0, 4733084.375, 4735140.625, 4740709.375, 4741070.3125, 4743775.0, 4747467.1875, 4760954.6875, 4760976.5625, 4770456.25, 4783200.0, 4790409.375, 4793589.0625, 4799028.125, 4814804.6875, 4820009.375, 4821987.5, 4822035.9375, 4823054.6875, 4826410.9375, 4827120.3125, 4827203.125, 4830959.375, 4830973.4375, 4832904.6875, 4833489.0625, 4844934.375, 4845545.3125, 4845893.75, 4849990.625, 4853407.8125, 4855856.25, 4858214.0625, 4863420.3125, 4865295.3125, 4866053.125, 4870170.3125, 4874395.3125, 4877328.125, 4880217.1875, 4880818.75, 4896176.5625, 4896715.625, 4902978.125, 4911889.0625, 4915531.25, 4917796.875, 4923159.375, 4928773.4375, 4935479.6875, 4943054.6875, 4945256.25, 4945687.5, 4947821.875, 4947912.5, 4952792.1875, 4953240.625, 4953317.1875, 4953471.875, 4955229.6875, 4958467.1875, 4958510.9375, 4964034.375, 4964064.0625, 4964178.125, 4969896.875, 4983842.1875, 4990490.625, 4992773.4375, 5001396.875, 5003035.9375, 5006504.6875, 5008365.625, 5008385.9375, 5015668.75, 5023470.3125, 5025200.0, 5028593.75, 5030446.875, 5032739.0625, 5045187.5, 5050703.125, 5054545.3125, 5054554.6875, 5054651.5625, 5057193.75, 5059929.6875, 5071673.4375, 5077125.0, 5080265.625, 5089387.5, 5094989.0625, 5107109.375, 5128735.9375, 5145960.9375, 5146853.125, 5153967.1875, 5158428.125, 5170628.125, 5197831.25, 5208023.4375, 5225720.3125, 5227214.0625, 5240971.875, 5267218.75, 5272814.0625, 5273248.4375, 5273268.75, 5398412.5, 5415448.4375, 5435407.8125, 5446914.0625, 5450184.375, 5453114.0625, 5454279.6875, 5458439.0625, 5459515.625, 5466901.5625, 5470607.8125, 5478843.75, 5478865.625, 5480057.8125, 5482807.8125, 5499718.75, 5515796.875, 5530859.375, 5531109.375, 5535317.1875, 5545223.4375, 5550046.875, 5554468.75, 5556876.5625, 5556904.6875, 5559548.4375, 5580939.0625, 5581068.75, 5582159.375, 5583506.25, 5585168.75, 5586623.4375, 5610334.375, 5617987.5, 5623110.9375, 5625445.3125, 5627564.0625, 5647457.8125, 5651701.5625, 5651725.0, 5654246.875, 5668335.9375, 5672264.0625, 5680914.0625, 5682459.375, 5689290.625, 5690996.875, 5706857.8125, 5707810.9375, 5721684.375, 5742173.4375, 5781223.4375, 5781284.375, 5799792.1875, 5805968.75, 5816207.8125, 5889990.625, 5895464.0625, 5941703.125, 5979157.8125, 5990748.4375, 6019798.4375, 6033767.1875, 6060998.4375, 6064048.4375, 6101187.5, 6110467.1875, 6130032.8125, 6132146.875, 6158523.4375, 6176171.875, 6211081.25, 6251334.375, 6269348.4375, 6277125.0, 6293648.4375, 6297175.0, 6300803.125, 6311553.125, 6321376.5625, 6337212.5, 6359751.5625, 6359839.0625, 6366684.375, 6368276.5625, 6392606.25, 6404895.3125, 6408560.9375, 6412671.875, 6412967.1875, 6432392.1875, 6434925.0, 6449143.75, 6453175.0, 6453189.0625, 6457040.625, 6464821.875, 6465721.875, 6467864.0625, 6469910.9375, 6481212.5, 6481223.4375, 6484256.25, 6490796.875, 6491153.125, 6492528.125, 6494725.0, 6498693.75, 6501935.9375, 6501946.875, 6504720.3125, 6512354.6875, 6520545.3125, 6520851.5625, 6522270.3125, 6523392.1875, 6527471.875, 6534904.6875, 6552546.875, 6557964.0625, 6570459.375, 6572651.5625, 6577250.0, 6584509.375, 6601348.4375, 6646229.6875, 6669004.6875, 6671432.8125, 6682762.5, 6700582.8125, 6702760.9375, 6746168.75, 6772123.4375, 6773493.75, 6778526.5625, 6781509.375, 6786373.4375, 6790454.6875, 6798782.8125, 6800720.3125, 6803750.0, 6809440.625, 6809489.0625, 6832151.5625, 6832212.5, 6878804.6875, 6880190.625, 6988062.5, 6989648.4375, 6991495.3125, 6997343.75, 7062939.0625, 7088926.5625, 7443489.0625, 7454237.5, 8496096.875, 8612093.75], [50.944314495704695, 15.796553761775797, 39.599775479332536, 8.137706685417054, 63.074565691841194, 80.21789700366425, 15.960937936626587, 12.512999259053608, 29.837789083558576, 53.120041011142696, 27.89621550330027, 47.63394783322977, 18.300284868169413, 50.57820505054154, 35.690030261814336, 67.36186904267599, 7.196413843423989, 65.38448733491555, 24.998896924739103, 6.509668300256832, 19.315350304061933, 5.439443571260751, 25.632901477304635, 66.16595145584803, 19.278302018446666, 14.15874738738331, 6.522623815425454, 107.53968102319064, 33.41513750724508, 5.316264066312792, 46.10586863436758, 91.16589071864006, 48.09769939310871, 7.158275626253481, 54.74312561294476, 6.794411236218008, 13.644069024890666, 9.542268560378988, 24.600628551372115, 8.980340001393829, 10.803338798417208, 40.45007002599705, 68.50365582357375, 16.959660412291317, 21.905485907458093, 8.268913915283843, 62.64292542871192, 91.47650348548203, 58.34117945242116, 25.060641910279976, 37.35763724901154, 40.520881028062504, 26.911704247992077, 5.121250559007611, 23.60418085712736, 5.700602386763834, 17.9621753510238, 59.68005251310863, 69.0006913012932, 10.955843315049043, 25.84869330680094, 94.43485138712042, 33.13091245400119, 23.792743755176126, 73.61615820046723, 45.69454257566443, 8.750624063703121, 71.5191867657087, 64.94882434073011, 5.711309894727182, 20.481644961295867, 38.950150823930734, 18.088430464499684, 5.647520900683906, 12.135772125361198, 13.752119500296295, 13.091904364232594, 27.168080232946714, 113.1331687422565, 12.86995546178879, 12.14858966370281, 44.947831610102234, 43.854046163781256, 78.75315976097336, 93.75752775515157, 52.20668232804204, 55.90487138715627, 16.09285818625129, 19.798294202270533, 60.54741029516812, 17.389647197718674, 15.9963206441607, 54.395748676954206, 23.103191481878326, 51.80352019490101, 59.57832523208183, 11.51942720465644, 28.086517632765965, 6.065950813993589, 10.13062257170731, 10.574647561438637, 21.320669548016888, 19.745649394005632, 61.228638772155804, 72.73013554124239, 5.2719012231358695, 32.57270000951705, 7.253053753444399, 7.280814963669682, 42.912864000026964, 67.61316065124427, 5.62951948229284, 7.085303980947569, 59.43974368411783, 17.31304622136468, 64.10192608274598, 11.492817322970527, 33.20667022877087, 15.024348292525623, 9.834152479081345, 9.56967415166326, 64.06353291501641, 79.69358726818837, 7.549817268686985, 73.97081699540749, 17.74799126293028, 49.266121924875584, 17.13425397130002, 27.214492930305425, 12.195626065470634, 94.37692473064298, 16.373007131649477, 42.14259295651776, 7.98458684063542, 75.55815453060413, 35.12623972643939, 49.602479270763496, 17.224787612594323, 5.169139747263604, 26.830700834371942, 25.67510155413373, 27.48320485520895, 44.32946137990008, 75.48996115018434, 58.447575040326555, 36.518753572604936, 5.521493061442035, 8.031471558821353, 27.821059716650893, 13.163987644669419, 25.786962165195888, 41.54410680529542, 21.374984613534597, 34.703267100647466, 88.53354610613425, 12.17820705532583, 45.08392716131503, 13.34068538552218, 10.626892004130326, 40.55367665331279, 5.029691580161599, 89.72498654448499, 26.63211685324429, 12.71011993697828, 17.96824894569834, 37.943651597656164, 70.22964903144896, 73.1276765720511, 106.92512423577564, 20.20393569463335, 15.249400727842824, 10.410655364025292, 50.12537437218167, 161.61939059722428, 5.21016275423511, 10.861592297274491, 13.98605953305139, 25.778956749256366, 14.351774903446644, 19.987153139850143, 22.08495909861419, 28.94689530639569, 73.53783070486082, 83.43072176960183, 78.51153103982357, 31.492972990266875, 5.470208713948725, 16.366000167455127, 33.47860789500325, 71.33295277617917, 15.776243049260728, 22.907238928496348, 83.02688974523193, 26.25698520129269, 5.786110069010675, 36.033767298774734, 8.339859353035278, 13.805229401299217, 6.588777894823672, 31.577525446498296, 72.44380337767733, 9.46734512134162, 50.183104440415796, 6.366286970021599, 23.376396962840523, 62.504425608304146, 19.384430652404177, 37.407907875019774, 16.118298272153197, 20.34051371970572, 6.8365284522596435, 11.223147344283218, 9.508576428819588, 78.63040007514823, 10.766912906583594, 61.598593402889236, 26.181738301365808, 78.04871585780177, 89.40377328229651, 40.29039994151957, 129.79098658785202, 63.755752766138386, 86.46792473340524, 125.78626764771685, 76.2125002110638, 25.498601563515233, 28.312596114138803, 30.941215692162736, 9.054369042689874, 14.164660251426868, 18.693462145242975, 19.899878014321146, 43.26646387948695, 27.972438283688014, 29.18301157608849, 54.57432964556627, 38.62558218029787, 65.49394062729058, 37.691739056808075, 18.710073706045748, 37.102507653289834, 99.38100862261118, 103.78187255075929, 105.32317128029207, 6.067994961857043, 121.59654713008884, 25.010919430937037, 13.252717642964797, 5.401295107188826, 6.249232438686556, 54.48134949199029, 9.790394190622914, 103.45283563733504, 37.51866728377649, 14.745446678476293, 59.74435588506283, 14.983347717471712, 9.514996551443158, 41.802267844644426, 19.66102646435677, 16.46179074001237, 66.92915642180607, 128.02455324041892, 72.49916501972811, 68.96663773029925, 72.90284714656786, 26.720954706378436, 11.33293797430555, 36.34627591652319, 60.40303690511706, 29.9650978750554, 117.99923337949089, 17.094445723828184, 91.70197062341235, 44.722346253657754, 36.51871247738389, 61.322620210664596, 5.438962070653053, 38.671147273977304, 97.95481064266868, 33.55581502171796, 65.953701445872, 6.723190745924475, 73.77858508677758, 15.71509270840959, 15.355024075580742, 62.652882573216765, 36.625712796389536, 6.672648004620804, 58.04449545218593, 15.370729463586787, 8.467784766132016, 100.66027502822826, 32.041961162279236, 17.438030165431037, 130.96867778309647, 41.12916289466618, 52.23345039876758, 123.86519551561767, 154.35490746532048, 11.658505625421697, 25.648847646213042, 71.56051894896395, 20.18543790349799, 77.60325287608424, 49.88799803276582, 10.627805648995272, 102.72709699499185, 18.495906855980817, 7.939345881218487, 10.741123252350581, 42.85624952299829, 34.96980933190555, 44.087608556740314, 42.82915315862963, 37.48159910805486, 7.110450032533722, 19.765006074314748, 20.785542498566095, 60.679918771112995, 106.13644172889829, 55.999691995288416, 31.898404678909362, 89.3954371006658, 55.01484533186828, 41.01534224610968, 23.07885957068808, 11.384386466817384, 8.50711056549528, 23.983318871884126, 15.096644262654838, 5.1155136764734035, 18.736755183604288, 7.855463055278699, 30.66548143806213, 38.60271901350903, 14.957982285149768, 27.11355877172525, 38.69326356259688, 83.38412336903227, 32.73809522327452, 14.112279019704951, 10.725791158170066, 47.124592597325226, 32.76567063421831, 114.23063205732157, 34.81246903405556, 53.16611174253132, 32.47627953177644, 33.81104503044368, 7.137695223752519, 27.84019008140716, 30.715824413239577, 61.91498337181847, 49.61342053035381, 13.640448684487664, 27.102532530390054, 20.126295418323753, 36.92291451479147, 11.564345437419183, 27.340257804324402, 34.38113355977191, 8.275789781157172, 52.96589467592817, 30.413434174258548, 63.35069299412984, 89.93052660655312, 11.097756967341901, 14.182264921613376, 25.582335321899002, 27.156526344122497, 41.35811313786884, 43.90833221232265, 16.230177861795134, 47.81085930995156, 6.876853992204028, 31.672761754490352, 5.955097910899297, 166.85813640451374, 11.411397247443487, 44.70645638228816, 18.07786856143491, 5.185435539395201, 65.95546909937444, 9.971028261650556, 118.39103987322837, 49.890096393814595, 69.41531786768181, 26.21951758573789, 21.23867454555524, 10.143994802986827, 32.51560401289473, 21.79621245669433, 48.00168823356475, 24.296253328550815, 74.27095379095478, 16.575567528892083, 73.627128048453, 85.16554642018879, 20.307248070382904, 210.8618319081322, 8.035702849446732, 12.086441242558092, 53.98103473962383, 39.99103512153062, 103.5195688754545, 17.35749940991072, 5.31059899485489, 7.2174439282374365, 23.5509437563909, 20.837841185910676, 125.89647693372366, 20.757245732220515, 16.995482224641496, 75.2894909688611, 76.38747288304455, 32.4561234393083, 165.50694487951301, 37.91771556902011, 75.95210886056685, 94.47482444903409, 37.51500454799752, 24.919236587801556, 23.86275841033838, 13.079024537964218, 66.35252802786141, 8.658965112661603, 74.07376097052445, 10.066941470899478, 14.029687147044973, 59.70204416174195, 12.460485992299345, 7.613495581663946, 10.829249025980673, 5.041197028084741, 8.570950591301235, 35.99194379275371, 5.26242857463036, 36.665630139726744, 117.88856411398243, 15.177450688591072, 18.95123630595777, 52.34620010855623, 80.47857381717395, 24.345119376864314, 53.46429305952985, 32.82346072044773, 20.55074504902428, 14.062462460633181, 7.690274984551698, 15.708740461407524, 122.16682973358604, 7.976266346969618, 19.726846351661077, 8.043416970914338, 85.00391104918818, 52.21749771569018, 131.72901701947387, 85.97567327941024, 47.05270627025925, 63.407819479049095, 32.019609506160634, 38.73526774545978, 47.40688913661187, 9.179524584071384, 13.146208003173333, 123.41685647243351, 59.04266508708121, 15.513365620683317, 5.974535652039651, 43.15826456774729, 5.803372783228023, 65.90218942643276, 19.813920478810942, 14.049400316570729, 103.34250304794251, 54.05178421181283, 51.1319095357435, 55.80989774829781, 100.2944279675149, 49.86491054820108, 12.781548653137662, 56.052358061165634, 8.221534051912824, 78.5618487473546, 39.09805077074164, 52.5747485104642, 12.250288110175825, 70.33840919877106, 10.984700483099775, 16.793637521933956, 28.57460566821389, 41.07797487023091, 55.059771295601095, 19.462488708631753, 5.593877153424963, 14.954774831355817, 6.033095841148297, 30.259375928780567, 5.48931152161992, 54.34734269218338, 33.779584732135454, 66.9756497119167, 18.67147411604182, 22.764441232765726, 5.198432572722792, 5.075672511810246, 69.42741829288413, 58.49251865541673, 10.939096943049515, 83.81919927495807, 33.220221302536594])
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);
([3326887.5, 3336750.0, 3430026.5625, 3457040.625, 3465812.5, 3468035.9375, 3469092.1875, 3470846.875, 3471004.6875, 3471117.1875, 3472431.25, 3473471.875, 3476746.875, 3477715.625, 3489785.9375, 3502885.9375, 3505543.75, 3512221.875, 3514873.4375, 3533598.4375, 3535301.5625, 3536723.4375, 3562912.5, 3577376.5625, 3580067.1875, 3580279.6875, 3586718.75, 3587570.3125, 3595275.0, 3602159.375, 3612182.8125, 3629662.5, 3637182.8125, 3637484.375, 3640037.5, 3642481.25, 3645756.25, 3666034.375, 3668587.5, 3669707.8125, 3673342.1875, 3677042.1875, 3686990.625, 3699889.0625, 3712414.0625, 3713782.8125, 3732732.8125, 3735007.8125, 3737292.1875, 3740648.4375, 3741112.5, 3742118.75, 3743825.0, 3744315.625, 3748340.625, 3751559.375, 3752273.4375, 3755212.5, 3755807.8125, 3758048.4375, 3767800.0, 3775487.5, 3775782.8125, 3777957.8125, 3778107.8125, 3779579.6875, 3781940.625, 3782976.5625, 3783531.25, 3785453.125, 3785595.3125, 3791892.1875, 3792293.75, 3792343.75, 3794779.6875, 3794807.8125, 3794885.9375, 3795123.4375, 3795201.5625, 3795425.0, 3796775.0, 3797351.5625, 3797432.8125, 3797748.4375, 3797817.1875, 3798037.5, 3798153.125, 3798317.1875, 3798760.9375, 3799110.9375, 3799131.25, 3799132.8125, 3799193.75, 3799443.75, 3800071.875, 3800340.625, 3800468.75, 3801090.625, 3801528.125, 3801806.25, 3802296.875, 3802773.4375, 3802925.0, 3803112.5, 3803467.1875, 3804554.6875, 3805026.5625, 3805328.125, 3806145.3125, 3806717.1875, 3826245.3125, 3833881.25, 3837154.6875, 3855659.375, 3862954.6875, 3890884.375, 3895787.5, 3906654.6875, 3913764.0625, 3949573.4375, 4017015.625, 4022571.875, 4030139.0625, 4057093.75, 4068350.0, 4080389.0625, 4081840.625, 4128017.1875, 4148826.5625, 4148854.6875, 4153829.6875, 4161646.875, 4166392.1875, 4180204.6875, 4191534.375, 4193710.9375, 4194979.6875, 4196157.8125, 4196228.125, 4199753.125, 4220971.875, 4222853.125, 4223217.1875, 4239856.25, 4258600.0, 4262967.1875, 4267545.3125, 4288862.5, 4294617.1875, 4294910.9375, 4303500.0, 4305481.25, 4305985.9375, 4308029.6875, 4310221.875, 4310242.1875, 4318239.0625, 4324614.0625, 4327668.75, 4331854.6875, 4337943.75, 4340582.8125, 4347409.375, 4367134.375, 4381517.1875, 4383151.5625, 4387140.625, 4398503.125, 4400676.5625, 4404487.5, 4404500.0, 4427906.25, 4433817.1875, 4435356.25, 4435442.1875, 4440901.5625, 4458620.3125, 4459748.4375, 4463070.3125, 4463907.8125, 4465964.0625, 4465967.1875, 4469293.75, 4473443.75, 4490329.6875, 4497471.875, 4500485.9375, 4501984.375, 4508179.6875, 4511987.5, 4514610.9375, 4514971.875, 4522603.125, 4529471.875, 4530121.875, 4530643.75, 4531968.75, 4533301.5625, 4536485.9375, 4538485.9375, 4544220.3125, 4544756.25, 4545223.4375, 4545935.9375, 4547050.0, 4547217.1875, 4548378.125, 4549014.0625, 4549131.25, 4550893.75, 4554931.25, 4563829.6875, 4569568.75, 4572073.4375, 4572143.75, 4572332.8125, 4574929.6875, 4574951.5625, 4575776.5625, 4579729.6875, 4582498.4375, 4585431.25, 4595532.8125, 4596210.9375, 4600998.4375, 4608214.0625, 4612792.1875, 4640540.625, 4646134.375, 4654810.9375, 4667539.0625, 4668748.4375, 4670723.4375, 4676304.6875, 4701053.125, 4708215.625, 4711590.625, 4714379.6875, 4715809.375, 4722696.875, 4723876.5625, 4724245.3125, 4724629.6875, 4724970.3125, 4725031.25, 4731400.0, 4732053.125, 4733025.0, 4733084.375, 4735140.625, 4740709.375, 4741070.3125, 4743775.0, 4747467.1875, 4760954.6875, 4760976.5625, 4770456.25, 4783200.0, 4790409.375, 4793589.0625, 4799028.125, 4814804.6875, 4820009.375, 4821987.5, 4822035.9375, 4823054.6875, 4826410.9375, 4827120.3125, 4827203.125, 4830959.375, 4830973.4375, 4832904.6875, 4833489.0625, 4844934.375, 4845545.3125, 4845893.75, 4849990.625, 4853407.8125, 4855856.25, 4858214.0625, 4863420.3125, 4865295.3125, 4866053.125, 4870170.3125, 4874395.3125, 4877328.125, 4880217.1875, 4880818.75, 4896176.5625, 4896715.625, 4902978.125, 4911889.0625, 4915531.25, 4917796.875, 4923159.375, 4928773.4375, 4935479.6875, 4943054.6875, 4945256.25, 4945687.5, 4947821.875, 4947912.5, 4952792.1875, 4953240.625, 4953317.1875, 4953471.875, 4955229.6875, 4958467.1875, 4958510.9375, 4964034.375, 4964064.0625, 4964178.125, 4969896.875, 4983842.1875, 4990490.625, 4992773.4375, 5001396.875, 5003035.9375, 5006504.6875, 5008365.625, 5008385.9375, 5015668.75, 5023470.3125, 5025200.0, 5028593.75, 5030446.875, 5032739.0625, 5045187.5, 5050703.125, 5054545.3125, 5054554.6875, 5054651.5625, 5057193.75, 5059929.6875, 5071673.4375, 5077125.0, 5080265.625, 5089387.5, 5094989.0625, 5107109.375, 5128735.9375, 5145960.9375, 5146853.125, 5153967.1875, 5158428.125, 5170628.125, 5197831.25, 5208023.4375, 5225720.3125, 5227214.0625, 5240971.875, 5267218.75, 5272814.0625, 5273248.4375, 5273268.75, 5398412.5, 5415448.4375, 5435407.8125, 5446914.0625, 5450184.375, 5453114.0625, 5454279.6875, 5458439.0625, 5459515.625, 5466901.5625, 5470607.8125, 5478843.75, 5478865.625, 5480057.8125, 5482807.8125, 5499718.75, 5515796.875, 5530859.375, 5531109.375, 5535317.1875, 5545223.4375, 5550046.875, 5554468.75, 5556876.5625, 5556904.6875, 5559548.4375, 5580939.0625, 5581068.75, 5582159.375, 5583506.25, 5585168.75, 5586623.4375, 5610334.375, 5617987.5, 5623110.9375, 5625445.3125, 5627564.0625, 5647457.8125, 5651701.5625, 5651725.0, 5654246.875, 5668335.9375, 5672264.0625, 5680914.0625, 5682459.375, 5689290.625, 5690996.875, 5706857.8125, 5707810.9375, 5721684.375, 5742173.4375, 5781223.4375, 5781284.375, 5799792.1875, 5805968.75, 5816207.8125, 5889990.625, 5895464.0625, 5941703.125, 5979157.8125, 5990748.4375, 6019798.4375, 6033767.1875, 6060998.4375, 6064048.4375, 6101187.5, 6110467.1875, 6130032.8125, 6132146.875, 6158523.4375, 6176171.875, 6211081.25, 6251334.375, 6269348.4375, 6277125.0, 6293648.4375, 6297175.0, 6300803.125, 6311553.125, 6321376.5625, 6337212.5, 6359751.5625, 6359839.0625, 6366684.375, 6368276.5625, 6392606.25, 6404895.3125, 6408560.9375, 6412671.875, 6412967.1875, 6432392.1875, 6434925.0, 6449143.75, 6453175.0, 6453189.0625, 6457040.625, 6464821.875, 6465721.875, 6467864.0625, 6469910.9375, 6481212.5, 6481223.4375, 6484256.25, 6490796.875, 6491153.125, 6492528.125, 6494725.0, 6498693.75, 6501935.9375, 6501946.875, 6504720.3125, 6512354.6875, 6520545.3125, 6520851.5625, 6522270.3125, 6523392.1875, 6527471.875, 6534904.6875, 6552546.875, 6557964.0625, 6570459.375, 6572651.5625, 6577250.0, 6584509.375, 6601348.4375, 6646229.6875, 6669004.6875, 6671432.8125, 6682762.5, 6700582.8125, 6702760.9375, 6746168.75, 6772123.4375, 6773493.75, 6778526.5625, 6781509.375, 6786373.4375, 6790454.6875, 6798782.8125, 6800720.3125, 6803750.0, 6809440.625, 6809489.0625, 6832151.5625, 6832212.5, 6878804.6875, 6880190.625, 6988062.5, 6989648.4375, 6991495.3125, 6997343.75, 7062939.0625, 7088926.5625, 7443489.0625, 7454237.5, 8496096.875, 8612093.75], [50.944314495704695, 15.796553761775797, 39.599775479332536, 8.137706685417054, 63.074565691841194, 80.21789700366425, 15.960937936626587, 12.512999259053608, 29.837789083558576, 53.120041011142696, 27.89621550330027, 47.63394783322977, 18.300284868169413, 50.57820505054154, 35.690030261814336, 67.36186904267599, 7.196413843423989, 65.38448733491555, 24.998896924739103, 6.509668300256832, 19.315350304061933, 5.439443571260751, 25.632901477304635, 66.16595145584803, 19.278302018446666, 14.15874738738331, 6.522623815425454, 107.53968102319064, 33.41513750724508, 5.316264066312792, 46.10586863436758, 91.16589071864006, 48.09769939310871, 7.158275626253481, 54.74312561294476, 6.794411236218008, 13.644069024890666, 9.542268560378988, 24.600628551372115, 8.980340001393829, 10.803338798417208, 40.45007002599705, 68.50365582357375, 16.959660412291317, 21.905485907458093, 8.268913915283843, 62.64292542871192, 91.47650348548203, 58.34117945242116, 25.060641910279976, 37.35763724901154, 40.520881028062504, 26.911704247992077, 5.121250559007611, 23.60418085712736, 5.700602386763834, 17.9621753510238, 59.68005251310863, 69.0006913012932, 10.955843315049043, 25.84869330680094, 94.43485138712042, 33.13091245400119, 23.792743755176126, 73.61615820046723, 45.69454257566443, 8.750624063703121, 71.5191867657087, 64.94882434073011, 5.711309894727182, 20.481644961295867, 38.950150823930734, 18.088430464499684, 5.647520900683906, 12.135772125361198, 13.752119500296295, 13.091904364232594, 27.168080232946714, 113.1331687422565, 12.86995546178879, 12.14858966370281, 44.947831610102234, 43.854046163781256, 78.75315976097336, 93.75752775515157, 52.20668232804204, 55.90487138715627, 16.09285818625129, 19.798294202270533, 60.54741029516812, 17.389647197718674, 15.9963206441607, 54.395748676954206, 23.103191481878326, 51.80352019490101, 59.57832523208183, 11.51942720465644, 28.086517632765965, 6.065950813993589, 10.13062257170731, 10.574647561438637, 21.320669548016888, 19.745649394005632, 61.228638772155804, 72.73013554124239, 5.2719012231358695, 32.57270000951705, 7.253053753444399, 7.280814963669682, 42.912864000026964, 67.61316065124427, 5.62951948229284, 7.085303980947569, 59.43974368411783, 17.31304622136468, 64.10192608274598, 11.492817322970527, 33.20667022877087, 15.024348292525623, 9.834152479081345, 9.56967415166326, 64.06353291501641, 79.69358726818837, 7.549817268686985, 73.97081699540749, 17.74799126293028, 49.266121924875584, 17.13425397130002, 27.214492930305425, 12.195626065470634, 94.37692473064298, 16.373007131649477, 42.14259295651776, 7.98458684063542, 75.55815453060413, 35.12623972643939, 49.602479270763496, 17.224787612594323, 5.169139747263604, 26.830700834371942, 25.67510155413373, 27.48320485520895, 44.32946137990008, 75.48996115018434, 58.447575040326555, 36.518753572604936, 5.521493061442035, 8.031471558821353, 27.821059716650893, 13.163987644669419, 25.786962165195888, 41.54410680529542, 21.374984613534597, 34.703267100647466, 88.53354610613425, 12.17820705532583, 45.08392716131503, 13.34068538552218, 10.626892004130326, 40.55367665331279, 5.029691580161599, 89.72498654448499, 26.63211685324429, 12.71011993697828, 17.96824894569834, 37.943651597656164, 70.22964903144896, 73.1276765720511, 106.92512423577564, 20.20393569463335, 15.249400727842824, 10.410655364025292, 50.12537437218167, 161.61939059722428, 5.21016275423511, 10.861592297274491, 13.98605953305139, 25.778956749256366, 14.351774903446644, 19.987153139850143, 22.08495909861419, 28.94689530639569, 73.53783070486082, 83.43072176960183, 78.51153103982357, 31.492972990266875, 5.470208713948725, 16.366000167455127, 33.47860789500325, 71.33295277617917, 15.776243049260728, 22.907238928496348, 83.02688974523193, 26.25698520129269, 5.786110069010675, 36.033767298774734, 8.339859353035278, 13.805229401299217, 6.588777894823672, 31.577525446498296, 72.44380337767733, 9.46734512134162, 50.183104440415796, 6.366286970021599, 23.376396962840523, 62.504425608304146, 19.384430652404177, 37.407907875019774, 16.118298272153197, 20.34051371970572, 6.8365284522596435, 11.223147344283218, 9.508576428819588, 78.63040007514823, 10.766912906583594, 61.598593402889236, 26.181738301365808, 78.04871585780177, 89.40377328229651, 40.29039994151957, 129.79098658785202, 63.755752766138386, 86.46792473340524, 125.78626764771685, 76.2125002110638, 25.498601563515233, 28.312596114138803, 30.941215692162736, 9.054369042689874, 14.164660251426868, 18.693462145242975, 19.899878014321146, 43.26646387948695, 27.972438283688014, 29.18301157608849, 54.57432964556627, 38.62558218029787, 65.49394062729058, 37.691739056808075, 18.710073706045748, 37.102507653289834, 99.38100862261118, 103.78187255075929, 105.32317128029207, 6.067994961857043, 121.59654713008884, 25.010919430937037, 13.252717642964797, 5.401295107188826, 6.249232438686556, 54.48134949199029, 9.790394190622914, 103.45283563733504, 37.51866728377649, 14.745446678476293, 59.74435588506283, 14.983347717471712, 9.514996551443158, 41.802267844644426, 19.66102646435677, 16.46179074001237, 66.92915642180607, 128.02455324041892, 72.49916501972811, 68.96663773029925, 72.90284714656786, 26.720954706378436, 11.33293797430555, 36.34627591652319, 60.40303690511706, 29.9650978750554, 117.99923337949089, 17.094445723828184, 91.70197062341235, 44.722346253657754, 36.51871247738389, 61.322620210664596, 5.438962070653053, 38.671147273977304, 97.95481064266868, 33.55581502171796, 65.953701445872, 6.723190745924475, 73.77858508677758, 15.71509270840959, 15.355024075580742, 62.652882573216765, 36.625712796389536, 6.672648004620804, 58.04449545218593, 15.370729463586787, 8.467784766132016, 100.66027502822826, 32.041961162279236, 17.438030165431037, 130.96867778309647, 41.12916289466618, 52.23345039876758, 123.86519551561767, 154.35490746532048, 11.658505625421697, 25.648847646213042, 71.56051894896395, 20.18543790349799, 77.60325287608424, 49.88799803276582, 10.627805648995272, 102.72709699499185, 18.495906855980817, 7.939345881218487, 10.741123252350581, 42.85624952299829, 34.96980933190555, 44.087608556740314, 42.82915315862963, 37.48159910805486, 7.110450032533722, 19.765006074314748, 20.785542498566095, 60.679918771112995, 106.13644172889829, 55.999691995288416, 31.898404678909362, 89.3954371006658, 55.01484533186828, 41.01534224610968, 23.07885957068808, 11.384386466817384, 8.50711056549528, 23.983318871884126, 15.096644262654838, 5.1155136764734035, 18.736755183604288, 7.855463055278699, 30.66548143806213, 38.60271901350903, 14.957982285149768, 27.11355877172525, 38.69326356259688, 83.38412336903227, 32.73809522327452, 14.112279019704951, 10.725791158170066, 47.124592597325226, 32.76567063421831, 114.23063205732157, 34.81246903405556, 53.16611174253132, 32.47627953177644, 33.81104503044368, 7.137695223752519, 27.84019008140716, 30.715824413239577, 61.91498337181847, 49.61342053035381, 13.640448684487664, 27.102532530390054, 20.126295418323753, 36.92291451479147, 11.564345437419183, 27.340257804324402, 34.38113355977191, 8.275789781157172, 52.96589467592817, 30.413434174258548, 63.35069299412984, 89.93052660655312, 11.097756967341901, 14.182264921613376, 25.582335321899002, 27.156526344122497, 41.35811313786884, 43.90833221232265, 16.230177861795134, 47.81085930995156, 6.876853992204028, 31.672761754490352, 5.955097910899297, 166.85813640451374, 11.411397247443487, 44.70645638228816, 18.07786856143491, 5.185435539395201, 65.95546909937444, 9.971028261650556, 118.39103987322837, 49.890096393814595, 69.41531786768181, 26.21951758573789, 21.23867454555524, 10.143994802986827, 32.51560401289473, 21.79621245669433, 48.00168823356475, 24.296253328550815, 74.27095379095478, 16.575567528892083, 73.627128048453, 85.16554642018879, 20.307248070382904, 210.8618319081322, 8.035702849446732, 12.086441242558092, 53.98103473962383, 39.99103512153062, 103.5195688754545, 17.35749940991072, 5.31059899485489, 7.2174439282374365, 23.5509437563909, 20.837841185910676, 125.89647693372366, 20.757245732220515, 16.995482224641496, 75.2894909688611, 76.38747288304455, 32.4561234393083, 165.50694487951301, 37.91771556902011, 75.95210886056685, 94.47482444903409, 37.51500454799752, 24.919236587801556, 23.86275841033838, 13.079024537964218, 66.35252802786141, 8.658965112661603, 74.07376097052445, 10.066941470899478, 14.029687147044973, 59.70204416174195, 12.460485992299345, 7.613495581663946, 10.829249025980673, 5.041197028084741, 8.570950591301235, 35.99194379275371, 5.26242857463036, 36.665630139726744, 117.88856411398243, 15.177450688591072, 18.95123630595777, 52.34620010855623, 80.47857381717395, 24.345119376864314, 53.46429305952985, 32.82346072044773, 20.55074504902428, 14.062462460633181, 7.690274984551698, 15.708740461407524, 122.16682973358604, 7.976266346969618, 19.726846351661077, 8.043416970914338, 85.00391104918818, 52.21749771569018, 131.72901701947387, 85.97567327941024, 47.05270627025925, 63.407819479049095, 32.019609506160634, 38.73526774545978, 47.40688913661187, 9.179524584071384, 13.146208003173333, 123.41685647243351, 59.04266508708121, 15.513365620683317, 5.974535652039651, 43.15826456774729, 5.803372783228023, 65.90218942643276, 19.813920478810942, 14.049400316570729, 103.34250304794251, 54.05178421181283, 51.1319095357435, 55.80989774829781, 100.2944279675149, 49.86491054820108, 12.781548653137662, 56.052358061165634, 8.221534051912824, 78.5618487473546, 39.09805077074164, 52.5747485104642, 12.250288110175825, 70.33840919877106, 10.984700483099775, 16.793637521933956, 28.57460566821389, 41.07797487023091, 55.059771295601095, 19.462488708631753, 5.593877153424963, 14.954774831355817, 6.033095841148297, 30.259375928780567, 5.48931152161992, 54.34734269218338, 33.779584732135454, 66.9756497119167, 18.67147411604182, 22.764441232765726, 5.198432572722792, 5.075672511810246, 69.42741829288413, 58.49251865541673, 10.939096943049515, 83.81919927495807, 33.220221302536594])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xle,Etot)
