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
import matplotlib.cm as cm
#from matplotlib.mlab import griddata
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" t3pa2cls_XV_pc - zkousim vyzobat vysoke energie (jednotlive interakce) - funkce single_interaction. Dale delam prumernou velikost stopy stopy interakce pro danou energii - funkce size_of_interactions_average. Dale delam spektra 2 casti vyboje podle zadaneho casu - primarne pro double breakdown, tj. funkce energy_spectra_doublebreakdown
#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 = 44395 #test discharge for which the notebook will definitely work
shot_no = 45028
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
calib[j][i]=float(word[j]) #j = radek, i = sloupec0 - pouze pro stavajici kalibraci - verze XV
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))
ClmNo.append(int(int(matrix_index[i]))//int(256)) #ver XV - bude pro novou kalibraci - nyn9 to bere 3patnou kalibraci jednotlivych pixelu (cislovani radku a sloupcu stejne jako v pixetu)
RowNo.append(int(int(matrix_index[i]))%int(256)) ##ver XV - bude super pro novou kalibraci - odpovida radkum a sloupcum v pixetu. Nyni bere spatne kalibrace pixelu (viz fce load_calib)
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 Timewalk_parameters_Si():
#SI - korekce na TimeWalk - parametry
A = -1.21988
B = 4.33638
C = 29.5075
D = 1
sigma_A=0.7013
sigma_B=0.1366
sigma_C=4.753
sigma_D=0
return A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D
def Timewalk(E,T):
i=0
A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
while i < len(T):
timewalk=float(0)
timewalk = (C / (E[i] - B)) + A
if(timewalk<0):
timewalk=0
T[i]=T[i]-timewalk
i=i+1
return E,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=[]
N_tot=[]
eventNo=0 #for high_energy_event function
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))
N_tot.append(len(TCl)) #new in ver. XV
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')
#eventNo=single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path)
i=i+1
t3pa_cls.close
return T_first, E_tot, N_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, drawstyle='steps-mid')
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.plot(T_hit, count, drawstyle='steps-mid')
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.fill_between(xmean, count,step='mid') #pre mid post
plt.plot(xmean, count, drawstyle='steps-mid')
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, drawstyle='steps-mid')
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].fill_between(x2, y2,step='mid') #pre mid post
ax[1].plot(x2, y2, drawstyle='steps-mid')
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
def single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path): # ver XV
#tuto funkci volam ve funkci file_t3pa_cls_new(file_t3pa_cls,T):
#This function is calling in function file_t3pa_cls_new(file_t3pa_cls,T)
high_E = 1500 #[keV] - zadani, jaka energie je pro me vysoka
if(Etot>=high_E or len(ECl)>15):
#if (len(ECl)>15): #len(ECl) = pocet pixelu dane interakce
print("Total energy of track: "+str(Etot)+" keV\n")
A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
######################
#SI - rekonstrukce z(t) - parametry
a=0.0193947
b=-0.0273934
sigma_a=0.0004163
sigma_b=0.009832
#######################
time=[]
sigma_time=[]
time_min = []
time_max=[]
z_rec=[]
z_rec_min=[]
z_rec_max=[]
sigma_z_rec=[]
for i in range(0, len(TCl)):
time.append(TCl[i] - Tfirst)
#chyba urceni casu = chyba timewalk ->
sigma_time.append(np.sqrt(((sigma_C / (pow((ECl[i] - B), D)))*(sigma_C / (pow((ECl[i] - B), D)))) + (((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))*((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))) + (((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))* ((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))) + (sigma_A*sigma_A)))
time_min.append(time[i] - sigma_time[i])
time_max.append(time[i] + sigma_time[i])
#//
#//rekonstrukce hloubky interakce
z_rec.append(a * time[i] + (b))
sigma_z_rec.append(np.sqrt(((time[i]*time[i])*(sigma_a*sigma_a)) + (sigma_b*sigma_b) + (a*a)*(sigma_time[i]*sigma_time[i])))
z_rec_min.append(z_rec[i] - sigma_z_rec[i])
z_rec_max.append(z_rec[i] + sigma_z_rec[i])
#//
single_interaction_file=str(path)+'event'+str(eventNo)+'.txt'
single_interaction_2D_fig=str(path)+'2D_event'+str(eventNo)+'.png'
single_interaction_3D_fig=str(path)+'3D_event'+str(eventNo)+'.png'
with open(single_interaction_file, "w", encoding="utf-8") as f:
f.write('#\n')
f.write('# x_mean y_mean z_rec z_rec_min z_rec_max sigma_z_rec\n')
f.write('#E_tot='+str(Etot)+'\n')
f.write('\n')
for m in range(0,len(TCl)):
#print(str(RowNoCl[m])+' '+str(ClmNoCl[m])+' '+str(z_rec[m])+' '+str(z_rec_min[m])+' '+str(z_rec_max[m])+' '+str(sigma_z_rec[m])+' '+str(ECl[m])+'\n')
f.write(str(RowNoCl[m])+' '+str(ClmNoCl[m])+' '+str(z_rec[m])+' '+str(z_rec_min[m])+' '+str(z_rec_max[m])+' '+str(sigma_z_rec[m])+' '+str(ECl[m])+'\n')
f.close
############################### plot
#binwidth=1
if ((max(RowNoCl)-min(RowNoCl))>=(max(ClmNoCl)-min(ClmNoCl))):
rang=max(RowNoCl)-min(RowNoCl)
#bins=(range((min(RowNoCl)-1), (max(RowNoCl)+3), binwidth), range((min(ClmNoCl)-1), (min(ClmNoCl)+rang+3), binwidth))
x_min=min(RowNoCl)-1
x_max=max(RowNoCl)+2
y_min=min(ClmNoCl)-1
y_max=min(ClmNoCl)+rang+2
else:
rang = max(ClmNoCl)-min(ClmNoCl)
#bins=(range((min(RowNoCl)-1), (min(RowNoCl)+rang+3), binwidth), range((min(ClmNoCl)-1), (max(ClmNoCl)+3), binwidth))
x_min=min(RowNoCl)-1
x_max=min(RowNoCl)+rang+2
y_min=min(ClmNoCl)-1
y_max=max(ClmNoCl)+2
###############$$$$$$$$$$$$$$$$$$$$
rows, cols = (x_max-x_min, x_max-x_min)
#arrECl = [[0]*cols]*rows
arrECl= [[0] * cols for _ in range(rows)]
for i in range (x_min,x_max):
for j in range (y_min,y_max):
for k in range(0,len(RowNoCl)):
if ((i==RowNoCl[k]) and (j==ClmNoCl[k])):
#print(str(i)+" "+str(j)+" "+str(RowNoCl[k])+" "+str(ClmNoCl[k])+"\n")
arrECl[i-x_min][j-y_min]=ECl[k]
arrECl=np.transpose(arrECl)
arrECl[arrECl == 0] = np.nan
#2D plot
plt.imshow(arrECl,cmap='jet', interpolation='none', origin=[(x_min-0.2),(y_min-0.2)], aspect='equal', extent=[x_min-0.2,x_max,y_min-0.2,y_max])
###############$$$$$$$$$$$$$$$$$$$$
cb=plt.colorbar(pad=0.005)
plt.tick_params(axis="x",direction="in")
plt.tick_params(axis="y",direction="in")
cb.ax.tick_params(axis='y', direction='in')
cb.set_label('Energy in pixel [keV]')
plt.xlabel('x [pixel]')
plt.ylabel('y [pixel]')
plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
plt.savefig(single_interaction_2D_fig, dpi = 1000)
plt.show()
plt.close()
############################### plot
############## 3D plot
# Creating figure
fig = plt.figure(figsize = (16, 9))
ax = plt.axes(projection ="3d")
# Add x, y gridlines
ax.grid(b = True, color ='grey', #
linestyle ='-.', linewidth = 0.3,
alpha = 0.2)
# Creating plot
sctt = ax.scatter3D(RowNoCl, ClmNoCl, z_rec,
alpha = 0.8,
c = ECl,
cmap = 'jet',
s = 4000, #size of bubbles in plot
marker ='o')
plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
ax.set_xlabel('x [pixel]', fontweight ='bold')
ax.set_ylabel('y [pixel]', fontweight ='bold')
ax.set_zlabel('z [mm]', fontweight ='bold')
cb=fig.colorbar(sctt, ax = ax, shrink = 0.5,aspect = 10,pad=0.00005) #
cb.ax.tick_params(axis='y', direction='in')
cb.set_label('Energy in pixel [keV]', fontweight ='bold')
plt.savefig(single_interaction_3D_fig, dpi = 1000)
# show plot
plt.show()
plt.close()
############## 3D plot
eventNo=eventNo+1
return eventNo
def size_of_interactions_average(Etot, Ntot,file_hist): #ver XV
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=[]
count=[]
nop=[] #number of pixels
nop_mean=[] #mean number of pixels in dE
sigma=[] #root mean square error of the mean - stredni kvadraticka chyba prumeru
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.append(0)
nop.append(0)
nop_mean.append(0)
sigma.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
nop[E_index]=nop[E_index]+Ntot[i] #ver XV
for l in range(0,pocet):
if nop[l]!=0:
nop_mean[l]=float(nop[l])/float(count[l])
#####Computing sigma
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])):
if nop_mean[E_index]!=0:
sigma[E_index]=sigma[E_index]+((float(Ntot[i])-float(nop_mean[E_index]))*(float(Ntot[i])-float(nop_mean[E_index])))
for l in range(0,pocet):
if nop[l]!=0 and count[l]>1:
sigma[l]=np.sqrt(sigma[l]/(float(count[l])*(float(count[l])-1)))
if nop[l]!=0 and count[l]==1:
sigma[l]=10
#####Computing sigma
#####write to file
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('# x_mean [keV] mean_size_of_interaction [pixel] sigma_mean_size_of_interaction\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(nop_mean[m])+' '+str(sigma[m])+'\n')
hist.close
#####write to file
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
#plt.fill_between(xmean, nop,step='mid',color='blue', alpha=0.3) #pre mid post
#plt.fill_between(xmean, nop,step='mid') #pre mid post
plt.plot(xmean, nop_mean, drawstyle='steps-mid')
plt.errorbar(xmean, nop_mean, yerr = sigma,fmt ='.', linewidth=0.7)
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Average size of interaction [pixel]')
plt.xlim(0,500)
plt.ylim(0,17)
#ax.set_yscale('log') #log scale y
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
plt.savefig(figure_size_interaction, dpi = 1000)
return
def energy_spectra_doublebreakdown(Etot, Tfirst, file_hist):
#this function separates discharge into two parts acording to time Tsep and creates energy spectrum for each part.
#This is for discharges with double brakdown
Tsep=14636000 # [ns]
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_first=[]
count_second=[]
for l in range(0,pocet):
count_first.append(0)
count_second.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])):
if(Tfirst[i]<Tsep):
count_first[E_index]=count_first[E_index]+1
else:
count_second[E_index]=count_second[E_index]+1
caption = '# x = energy in keV, dE= '+str(dE)+' keV'
with open(file_hist, "w", encoding="utf-8") as hist:
hist.write('#\n')
hist.write('#'+str(caption)+'\n')
hist.write('# x_mean count_before_2nd_breakdown count_2nd_breakdown\n')
hist.write('\n')
for m in range(0,len(xmean)):
hist.write(str(xmean[m])+' '+str(count_first[m])+' '+str(count_second[m])+'\n')
hist.close
plt.rcParams.update({'font.size': 20})
fig, ax = plt.subplots(figsize =(10, 7))
#plt.fill_between(xmean, count_first,step='mid') #pre mid post
plt.plot(xmean, count_first, label='Before 2nd breakdown', drawstyle='steps-mid')
plt.plot(xmean, count_second, label='2nd breakdown', drawstyle='steps-mid')
plt.legend(loc='best')
plt.title(detector+', #'+str(shot_no))
plt.xlabel('Energy [keV]')
plt.ylabel('Count')
plt.xlim(0,500)
ax.set_yscale('log') #log scale y
plt.savefig(figure_energy_spectra_doublebreakdown, 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'
size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction.txt'
file_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown.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'
figure_size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction'
figure_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown'
#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)
E,T=Timewalk(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, Ntot = 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);
([7341949.345040617, 7349617.1875, 7383835.9375, 7393606.25, 7397523.257230354, 7445434.712954623, 7468078.458700638, 7590140.625, 7652841.934960744, 7652864.652044484, 7661440.2490893705, 7687321.802820627, 7690312.09380085, 7692894.330288865, 7745943.105911829, 7750797.360693115, 7756640.625, 7782889.0625, 7801279.2775708325, 7819941.042632043, 7822439.6355259335, 7836283.059348688, 7840826.5625, 7841531.25, 7844158.80428424, 7846588.273996023, 7848017.1875, 7852425.4068670105, 7852495.657854303, 7853495.3125, 7854843.397894762, 7863713.753118767, 7890918.37737375, 7898289.484058416, 7906446.146062636, 7942191.005106731, 7942281.25, 7942853.863910955, 7968112.160764721, 7969794.821085451, 7979844.091446953, 7985475.01446025, 7987010.9375, 8034182.721701814, 8036230.930784949, 8041966.674677169, 8042772.023527585, 8044907.176177609, 8045836.646008059, 8047148.4375, 8051003.238079173, 8065501.418240794, 8066104.64356547, 8101684.375, 8112638.348822645, 8114563.742646959, 8116540.651892431, 8117282.735175535, 8130018.498716827, 8134465.59919873, 8138842.736031749, 8139909.682297353, 8139931.322275146, 8140210.622996046, 8149181.222058955, 8150476.5625, 8152705.231142673, 8156530.836047882, 8186527.871475717, 8208695.339193121, 8213703.3227339415, 8226792.077911775, 8227906.596938308, 8237363.329474506, 8243156.25, 8249196.369103571, 8266282.997052491, 8272882.8125, 8287762.229380851, 8290214.681895636, 8291995.3125, 8305216.716505616, 8323184.362702257, 8324557.8125, 8327392.430370559, 8333623.874711879, 8345696.293689768, 8356584.044487525, 8356648.5906556435, 8356681.6637081485, 8375821.3815600155, 8382733.265861373, 8385429.248776955, 8394089.175476993, 8397533.986457879, 8399313.0454553, 8407952.949725354, 8425510.240162669, 8427825.849235635, 8442885.27843374, 8461856.341659928, 8464745.3125, 8470284.375, 8471206.71784104, 8471648.4375, 8513756.25, 8525624.588513607, 8540905.44721093, 8541380.407720702, 8542445.3125, 8545526.148771456, 8545727.246535618, 8551882.676993605, 8560948.013776204, 8561000.15712115, 8565946.875, 8578140.625, 8580568.63134677, 8580911.433435109, 8582014.157119911, 8600857.8125, 8603903.125, 8603962.174638461, 8619068.637103815, 8622017.244629288, 8626520.3125, 8628428.347521422, 8630178.907587143, 8632025.661892401, 8638678.79559247, 8640939.014642905, 8642656.25, 8661514.0625, 8661572.331479887, 8730926.784285095, 8731486.49317093, 8743176.011245739, 8751563.121227873, 8751636.42294095, 8788279.591452878, 8791042.377614316, 8794171.895300547, 8819184.309256323, 8821583.071259731, 8822595.604458127, 8838859.47066396, 8849970.800846713, 8852687.791912725, 8852709.481128328, 8853015.565010209, 8855471.07453266, 8857821.875, 8860770.289181355, 8870293.75, 8870311.858334515, 8876267.774277605, 8885243.295197597, 8888593.75, 8890747.79732472, 8917227.979842346, 8920241.794157727, 8930942.86470136, 8940040.555770405, 8959222.915736785, 8969429.6875, 8972732.8125, 8972741.79836481, 8990086.848588351, 8992370.3125, 9023863.620321712, 9033924.805646012, 9092160.653829597, 9095844.000911955, 9163353.494224448, 9188182.137801256, 9258283.678236874, 9286639.502158325, 9295834.375, 9299207.407273667, 9343423.261276508, 9360266.18421911, 9410804.47186099, 9410820.22457622, 9432051.879411984, 9473175.195609378, 9476643.75, 9508075.655324738, 9511889.708365962, 9613874.511384714, 9625040.558625018, 9625164.20291628, 9638977.373930046, 9666140.045569753, 9671676.866775503, 9673257.8125, 9679532.717368454, 9680727.452920264, 9681359.375, 9705052.956186488, 9707407.226047626, 9709420.702200953, 9939953.835931828, 9947037.500886794, 9947582.143217256, 9947616.254266793, 9947713.037218357, 10098756.722034426, 10114392.1875, 10140178.54214567, 10168673.030256838, 10245400.45599091, 10250363.579958308, 10250829.060798533, 10341917.016452713, 10341929.948088102, 10421804.62438267, 10456844.32915452, 10476242.949504135, 10493001.642198984, 10575321.064422602, 10598941.13536112, 10669774.496317685, 10734832.8125, 10794357.8125, 10824920.083542883, 10835825.0, 10856231.25, 10860109.301672088, 10860204.501909243, 10931707.8125, 10951379.102108477, 10987262.5, 10994413.248427745, 10999055.153642982, 11006610.9375, 11037250.353905322, 11096368.75, 11128470.59811617, 11212479.359960804, 11212480.35541019, 11268659.365610491, 11382613.633514645, 11412055.952972887, 11485456.690445792, 11510798.231705148, 11679995.96798333, 11758745.3125, 11759409.285680968, 11786719.579370689, 11907971.875, 11947923.474890647, 12702501.55572402, 12703397.614465257, 12770825.85181278, 12770833.96459695, 12903378.91299109, 13161077.118573125, 13362243.018636232, 13812197.390679702, 14130356.279429603, 15073546.476703227, 15703447.85121397, 15967063.132790694, 15975445.285061106, 15983375.0, 16000438.781963458, 16003400.546629013, 16004225.778088594, 16006234.375, 16012029.6875, 16013205.188377704, 16014734.181089526, 16019921.52185321, 16020284.375, 16020746.650079168, 16021970.430315385, 16024862.201077795, 16027269.944995206, 16027942.69561429, 16028143.628686234, 16028153.763264634, 16029173.688740635, 16029797.669739647, 16030219.688424889, 16033219.394997867, 16033374.035032988, 16033861.266319076, 16034277.264519773, 16044696.792563789, 16047613.596988507, 16049559.375, 16053158.701715937, 16057323.4375, 16061267.1875, 16076061.924175762, 16089310.384396724, 16091113.047003359, 16094590.01395407, 16101693.746823965, 16102161.248012103, 16103727.077608466, 16105836.81197619, 16111206.882711843, 16117304.795089623, 16123011.420360623, 16123178.125, 16124209.062846482, 16126799.98532281, 16131457.48823512, 16159658.292866107, 16159670.459857572, 16170748.37696521, 16178816.696698017, 16180820.151172316, 16181253.198424583, 16182625.0, 16182793.581977692, 16185856.441664541, 16188151.471513303, 16190857.8125, 16190873.534141613, 16192794.857072463, 16193136.99130448, 16193942.820979508, 16194149.577890957, 16194628.672531167, 16196087.070173824, 16197832.8125, 16197846.298173515, 16198335.408474391, 16198353.739516448, 16199219.238548953, 16201475.934689812, 16202189.939807765, 16202498.294068523, 16202552.388944034, 16202765.796936864, 16202922.949627612, 16203028.125, 16203054.361786857, 16203276.478806395, 16203492.1875, 16203548.256782122, 16203582.047625408, 16203622.925476385, 16203639.447967488, 16203903.738182325, 16204134.375, 16204141.513762644, 16204304.563938085, 16204310.049756793, 16204352.060637334, 16204420.3125, 16204462.208588608, 16204814.957411963, 16204821.875, 16204829.436457071, 16204851.5625, 16205209.81185865, 16205223.328967564, 16205364.56513149, 16205386.99308477, 16205494.877499951, 16205549.169438977, 16205585.395114595, 16205632.672856508, 16205645.111212235, 16205649.58967072, 16205675.28440221, 16205694.04788572, 16205695.403249042, 16205731.17083987, 16205760.587224947, 16205779.095125193, 16205788.597543662, 16205902.998795146, 16205915.12934558, 16205945.745684996, 16206019.71295666, 16206026.413997134, 16206113.145155473, 16206173.34115582, 16206324.664453605, 16206366.802089235, 16206410.614937702, 16206432.408029495, 16206436.292292394, 16206457.124816453, 16206571.200812276, 16206602.660124427, 16206608.82308006, 16206619.632093782, 16206720.3125, 16206721.556366265, 16206736.821106702, 16206744.747951135, 16206854.485657817, 16206862.298999459, 16207050.144014046, 16207055.151596809, 16207058.075438708, 16207062.734065278, 16207088.636386957, 16207097.452865582, 16207107.455154173, 16207110.662589548, 16207120.465474168, 16207151.122075, 16207154.6384436, 16207156.25, 16207159.54104702, 16207176.5625, 16207181.25, 16207352.933212722, 16207372.764958574, 16207382.96349704, 16207385.53057189, 16207415.03539173, 16207417.088251065, 16207423.4375, 16207429.6875, 16207440.368919, 16207443.298272546, 16207446.291629799, 16207461.812393589, 16207483.534085589, 16207530.825421393, 16207611.600039585, 16207624.50085341, 16207668.75, 16207681.81290489, 16207744.172887065, 16207803.02782234, 16207812.5, 16207839.201977473, 16207923.373718219, 16207930.980923822, 16207934.99277117, 16207943.75, 16208186.662612237, 16208288.780427814, 16208329.4307924, 16208353.125, 16208365.414823158, 16208368.18666921, 16208399.749228103, 16208425.08292508, 16208432.206033265, 16208434.702143043, 16208459.357214125, 16208482.308883851, 16208561.11001401, 16208565.139014997, 16208575.0, 16208585.86209186, 16208662.142132761, 16208669.510600239, 16208681.997399189, 16208691.23561566, 16208706.255704407, 16208709.102339998, 16208746.991572937, 16208828.849811837, 16208888.42539675, 16208902.028963516, 16208910.032445226, 16208921.434395267, 16209000.10386588, 16209008.928752225, 16209010.9375, 16209140.625, 16209171.856984854, 16209174.444345938, 16209185.917648746, 16209228.891844293, 16209236.442323564, 16209240.625, 16209264.0625, 16209272.274624782, 16209289.7936305, 16209346.840639606, 16209370.224296663, 16209404.861381564, 16209420.3125, 16209424.626553623, 16209490.304140162, 16209519.802903652, 16209544.757913612, 16209607.8125, 16209611.184333507, 16209623.603543133, 16209638.075376598, 16209649.329802962, 16209660.9375, 16209670.62542617, 16209680.683055287, 16209746.403555525, 16209837.552143583, 16209845.643559169, 16209862.37236405, 16209874.260214683, 16209875.0, 16209886.866938068, 16210000.013060093, 16210028.742586, 16210082.079162626, 16210198.688931327, 16210348.134684464, 16210400.726020444, 16210431.25, 16210441.09427573, 16210481.986191746, 16210529.850021752, 16210629.6875, 16210653.125, 16210709.375, 16210782.145125387, 16210928.187477548, 16211014.049246574, 16211250.0, 16211253.35827926, 16211274.098878387, 16211313.474744432, 16211343.812054962, 16211432.8125, 16211556.478196135, 16211939.0625, 16212947.216962518, 16215405.194514401, 16227408.388777623, 16228032.466202414, 16228560.9375, 16228776.138028476, 16237829.56004968, 16241242.1875, 16270035.46062934, 519606755.8849367], [20.075061542102226, 66.77501880346156, 71.30921418376734, 33.81409237767746, 7.074857048482096, 11.70050516397771, 15.034745588959952, 29.53884074984992, 37.839303192928384, 12.193867000567472, 31.521251638669202, 46.03448496680552, 13.590501857318467, 29.004085638048977, 27.964998497182975, 9.77887430555006, 49.72566040428604, 97.72299844866785, 31.545103971469203, 5.414683372178093, 11.769160038570622, 26.153907584027376, 47.03098597585904, 61.75344767832768, 105.92749539690658, 22.15326427365694, 71.73132744910212, 7.75733595205109, 5.839983218473605, 53.975902026678945, 19.50700929876801, 52.06737274794605, 22.865349671609543, 51.69173578569769, 7.93537467089529, 16.619537116783736, 45.10979161006365, 13.706026387527274, 5.035873157674501, 15.732565242874298, 16.424995062909233, 14.996913058459718, 53.026874073126315, 7.57072436982771, 9.075326368024315, 21.366132287659628, 6.2861479236927185, 19.730165168515622, 25.171883864739012, 50.09488959457041, 33.40238010010333, 24.26566169758266, 6.261663468686628, 72.56414211260218, 46.89414481711042, 60.36833824809907, 13.49815650227515, 7.101272679523748, 14.063069587278918, 64.33326313971097, 17.54564790495636, 17.400875238925167, 33.81211326632837, 25.641419217359164, 45.959500570241744, 32.60037403722289, 6.338417172140698, 5.513321849800691, 144.77660776556993, 6.09087947205094, 15.418885785530962, 25.341659901650218, 35.50541418600724, 71.27605716641206, 39.36457352126293, 22.356707018424903, 65.42515335689728, 68.2575927594729, 24.133439268042913, 47.76671146708033, 33.40476111405811, 126.96022669926693, 28.28381792259089, 74.6627860206372, 25.777967968611875, 40.406298006575696, 27.580978772542707, 17.60683442732823, 24.149425848278767, 29.955654198862227, 19.937313845276584, 7.245954106744644, 22.12695387466793, 15.390349713919543, 61.72575615106754, 13.837186854878645, 90.43150041250664, 19.72717557096077, 102.87459481374802, 10.233225902257045, 5.820793588919748, 29.095434121483848, 87.87406886971763, 8.550509651572053, 83.28020402227706, 39.59866491728035, 69.78283933231832, 32.517980802263274, 49.67115507706114, 38.368564654080096, 6.046106477957383, 18.39865703019591, 36.12005183428919, 175.4001918731995, 5.207443454605349, 31.64684812583591, 31.19065878790927, 22.562619910795718, 17.753333098417098, 64.35348765113199, 53.65973496730773, 32.259431565601645, 6.22694821099251, 10.955711877036798, 11.218193327512111, 106.67174755282105, 11.494296869717177, 72.70579783879532, 17.660134385819255, 18.309139839309328, 40.16919320368989, 64.21528181476219, 34.92093284491044, 45.82482037527302, 17.625145695633268, 5.85627194580403, 31.890368612354873, 11.782869552397484, 17.182824492849548, 35.08887946940321, 26.87711449273235, 63.38896079721869, 60.07059578532748, 42.165243682810924, 80.25058160964015, 38.874103814829205, 5.6416442625130125, 36.82559897347932, 15.36206204366398, 32.96566193837316, 28.06836206156396, 124.79263185035319, 61.12221867682287, 55.398488785627364, 27.36035219431543, 19.49389610603575, 100.58937753809113, 57.60619097843332, 55.93203579373673, 32.523086258791146, 18.624159832796064, 18.353005213084543, 37.94032440378395, 21.27871339131735, 78.75600929069337, 98.59560620718261, 84.97863253178672, 25.842246402921493, 39.33559768338594, 22.089969157499866, 61.39854032849083, 73.23072089807604, 15.992659867152057, 48.40359176537403, 132.7818624073178, 19.731786470195498, 5.735322272878144, 64.09168706998275, 94.6211113483372, 153.28901569714634, 9.853697188911143, 24.891663336309133, 7.571743723751527, 19.33341907526585, 15.743460218600939, 104.48986857739114, 18.208845155676997, 6.899683133084568, 21.60742561754838, 56.75489258231201, 6.791867150028416, 154.6530107752742, 21.49315371002927, 17.041341153779733, 82.43660062200432, 32.7731206655763, 49.141084537808126, 81.92493486169869, 68.91706542036114, 8.880306959568047, 61.38512717000971, 17.137335378912173, 117.29381605237185, 75.63676544934401, 18.040954334545884, 5.644472144757256, 15.570333763844427, 40.8827060879917, 9.710922911234558, 27.66903185465518, 28.530312013073367, 62.17721225319873, 66.97693539459904, 52.66803061030485, 40.440967609082264, 6.750970765136743, 17.72924435091622, 18.94133570498829, 80.11074487024763, 200.4687435638146, 9.803750163539398, 27.86735799574362, 36.03563045323571, 111.40771153564074, 14.135183274055848, 28.705418623390603, 49.86140889434661, 27.153670018917143, 25.3311248821529, 44.19954400207433, 13.098108576588519, 88.86851163005997, 28.662640324088002, 72.7412717629813, 33.08871361322505, 8.48303612720662, 134.53924851447078, 34.726314967854066, 13.82456864405277, 5.999761167384706, 14.90584138162558, 10.517429905455879, 84.07107132390601, 5.554845547520138, 46.29439608249341, 32.29971596611013, 46.244324574568296, 175.48145717355803, 51.9993133909239, 74.11136151973685, 84.24985240924052, 35.28629817457856, 66.82152995479836, 73.41750955849491, 15.380760736444621, 6.372167792201014, 23.103004882660066, 50.34380366030459, 22.609226339947014, 15.054880423896535, 16.243776038191832, 64.02658086600162, 5.650037852453839, 37.13493187479041, 58.45772323423543, 39.80891410031519, 6.8777746120311525, 12.609221508719802, 36.937759811044494, 29.134428388594007, 17.269744148317944, 14.25056688889304, 99.88607411522486, 30.778874450872873, 59.417979753494855, 7.6464126573571995, 15.618594619492722, 13.704182926784632, 5.418260638441956, 31.29959852148355, 34.161478163168084, 8.424061696299587, 29.54350733939834, 20.338667587948287, 12.311634057422053, 7.834665144646622, 52.561913777985836, 74.36932901112435, 26.9940876502495, 72.67222274453808, 37.157181032647024, 19.922720651328536, 54.473683572012504, 55.0145931336484, 44.45171381886674, 32.222463101184005, 17.53661635240796, 27.22952767156285, 34.72988744855627, 46.150935049855974, 6.14335567859943, 19.80230646197144, 12.285238657815746, 44.86213354232169, 29.948216647436638, 81.73477797450509, 15.561876071655293, 7.597937806432295, 13.834568928806593, 9.773532103381003, 21.43456418034761, 11.034387205159277, 21.58534871972576, 26.628846719737574, 64.62822227149175, 37.675655258466946, 25.596882135503552, 12.701804409578568, 80.94868902357908, 57.38661800315644, 15.323116509784889, 17.136891162234, 37.646783322989066, 39.13811489760208, 130.23304259792008, 27.929973354612947, 22.2228984705647, 62.898221684290874, 21.40617996414875, 5.6738967712777795, 9.91133951037937, 46.92876369831072, 12.989123212956686, 33.67937179737288, 10.910678631548743, 19.422507816416054, 6.289168527409597, 19.420527062482538, 56.731348456075985, 15.816455044057177, 14.631823894047551, 38.86662298407216, 13.815589868011408, 19.20345732481752, 26.591824484836582, 9.680090161869535, 44.11647870041121, 49.02297210020268, 50.73596560627708, 38.32281255368909, 6.3563869225383405, 17.25422819848931, 47.02723714357127, 130.56758244432976, 8.824234517119574, 64.21033530136208, 6.71440178167111, 52.37628249537762, 11.886876112065666, 113.5208912850465, 21.414626888600296, 6.66635568855736, 13.50766208408557, 12.5035490192224, 21.080452662489915, 33.95416016943703, 14.226053488610662, 14.24286990115944, 30.258168724214315, 5.547714671452123, 5.247692290593898, 64.9643512407398, 21.162318443445116, 20.618582634751206, 30.05812513821806, 117.61713861795866, 8.040776772372965, 134.31850932376992, 76.03539171963423, 14.40416549157068, 12.311973797609056, 25.113596484836293, 13.80020186926111, 13.651231320402653, 13.839777217641448, 15.226517806220984, 25.135747668609483, 8.810559784443628, 74.9187770416099, 7.443330189210997, 37.476861465854704, 6.95338673735558, 101.67118240654064, 77.40504268766605, 19.876676478015327, 6.402090871412261, 25.09113810038808, 55.46435473173046, 17.22753825294696, 5.215409166688883, 12.470010538079146, 26.696496578698653, 10.521151123570466, 7.282117998295218, 57.67719575072533, 36.54931471150268, 16.768363719918703, 8.066640717072982, 9.290264843750602, 89.3349235119259, 15.614569824684533, 97.65430615834937, 45.21188735764461, 9.174335520454225, 10.217388888177933, 9.462428806807115, 17.918982783953684, 18.16725921679945, 5.806673010592574, 51.62910037636868, 36.66838772817987, 5.238034194377177, 45.85564508243577, 78.48594679646877, 6.184064728162417, 24.591855045419287, 22.279975467008267, 6.596219987437982, 5.777948306581928, 33.757052308011886, 12.138520191598102, 9.716549205716694, 38.411936037576105, 32.17676677481561, 11.352949711372947, 7.5803306584315475, 6.4457225145365316, 12.253374468235686, 41.587137707384564, 23.374457897237974, 56.052981970018116, 59.096403009494, 58.58714815112546, 18.182677600021197, 13.155881101199387, 6.448486164988477, 51.03462439844092, 5.66921817148592, 9.624237885748267, 6.086251593232763, 22.591174941851676, 9.481661519838582, 48.09399364410018, 65.79830953190535, 20.891795440486142, 13.732933120846967, 7.903629232348305, 41.265466380423526, 79.20365857734731, 26.426247906280242, 15.568340752967961, 6.7870047665510675, 24.399169476286723, 12.965607821430828, 6.562787057628997, 23.558294761021227, 42.59169759231611, 5.294306419879466, 8.980577149229907, 32.8987019549005, 56.66240458079238, 28.17321310933923, 35.90856066221975, 68.42235815537283, 6.92690785216354, 62.19593855137672, 58.01074279079167, 71.04271192827615, 16.72015310948123, 7.890965366385145, 9.302517535880357, 5.9317608789624225, 5.101893962485382, 73.1577357838217, 5.90753024694414, 7.4911691974029635, 21.3979027250536, 26.380956496829693, 74.82327331059268, 22.771014108389085, 39.826113003446906, 22.205780406768014, 19.94817751189301, 42.15099726849792, 8.459293281230426, 13.146364855603132, 92.7986506850038, 5.099483888018506, 24.831513003166272, 72.20748546458901, 12.714039006344375, 50.75934390203771, 56.254828465852796, 7.607968136013016, 12.252981159405593, 5.115132740620922, 15.94402985077545, 13.90060361165151, 58.704611503952044, 41.84505547880708, 111.12550442527078, 10.042515331068806, 6.287952529110162, 64.8233749777597, 29.650362774678168, 96.35057326621563, 29.243951489237286, 19.548042522869395, 30.881475707153676, 32.622616087897434, 9.536763523497997, 25.323146319989455, 10.318476100292251, 6.496863087005234, 49.018398368786585, 5.544259074573224, 90.1785461420169, 32.037353455483014, 11.217234774533058, 28.110425575302507, 58.906534948122726, 34.30405590000561, 17.533230097891433, 41.59450720529936, 66.04597200938152, 73.98664303404462, 22.484793953950756])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([7341949.345040617, 7349617.1875, 7383835.9375, 7393606.25, 7397523.257230354, 7445434.712954623, 7468078.458700638, 7590140.625, 7652841.934960744, 7652864.652044484, 7661440.2490893705, 7687321.802820627, 7690312.09380085, 7692894.330288865, 7745943.105911829, 7750797.360693115, 7756640.625, 7782889.0625, 7801279.2775708325, 7819941.042632043, 7822439.6355259335, 7836283.059348688, 7840826.5625, 7841531.25, 7844158.80428424, 7846588.273996023, 7848017.1875, 7852425.4068670105, 7852495.657854303, 7853495.3125, 7854843.397894762, 7863713.753118767, 7890918.37737375, 7898289.484058416, 7906446.146062636, 7942191.005106731, 7942281.25, 7942853.863910955, 7968112.160764721, 7969794.821085451, 7979844.091446953, 7985475.01446025, 7987010.9375, 8034182.721701814, 8036230.930784949, 8041966.674677169, 8042772.023527585, 8044907.176177609, 8045836.646008059, 8047148.4375, 8051003.238079173, 8065501.418240794, 8066104.64356547, 8101684.375, 8112638.348822645, 8114563.742646959, 8116540.651892431, 8117282.735175535, 8130018.498716827, 8134465.59919873, 8138842.736031749, 8139909.682297353, 8139931.322275146, 8140210.622996046, 8149181.222058955, 8150476.5625, 8152705.231142673, 8156530.836047882, 8186527.871475717, 8208695.339193121, 8213703.3227339415, 8226792.077911775, 8227906.596938308, 8237363.329474506, 8243156.25, 8249196.369103571, 8266282.997052491, 8272882.8125, 8287762.229380851, 8290214.681895636, 8291995.3125, 8305216.716505616, 8323184.362702257, 8324557.8125, 8327392.430370559, 8333623.874711879, 8345696.293689768, 8356584.044487525, 8356648.5906556435, 8356681.6637081485, 8375821.3815600155, 8382733.265861373, 8385429.248776955, 8394089.175476993, 8397533.986457879, 8399313.0454553, 8407952.949725354, 8425510.240162669, 8427825.849235635, 8442885.27843374, 8461856.341659928, 8464745.3125, 8470284.375, 8471206.71784104, 8471648.4375, 8513756.25, 8525624.588513607, 8540905.44721093, 8541380.407720702, 8542445.3125, 8545526.148771456, 8545727.246535618, 8551882.676993605, 8560948.013776204, 8561000.15712115, 8565946.875, 8578140.625, 8580568.63134677, 8580911.433435109, 8582014.157119911, 8600857.8125, 8603903.125, 8603962.174638461, 8619068.637103815, 8622017.244629288, 8626520.3125, 8628428.347521422, 8630178.907587143, 8632025.661892401, 8638678.79559247, 8640939.014642905, 8642656.25, 8661514.0625, 8661572.331479887, 8730926.784285095, 8731486.49317093, 8743176.011245739, 8751563.121227873, 8751636.42294095, 8788279.591452878, 8791042.377614316, 8794171.895300547, 8819184.309256323, 8821583.071259731, 8822595.604458127, 8838859.47066396, 8849970.800846713, 8852687.791912725, 8852709.481128328, 8853015.565010209, 8855471.07453266, 8857821.875, 8860770.289181355, 8870293.75, 8870311.858334515, 8876267.774277605, 8885243.295197597, 8888593.75, 8890747.79732472, 8917227.979842346, 8920241.794157727, 8930942.86470136, 8940040.555770405, 8959222.915736785, 8969429.6875, 8972732.8125, 8972741.79836481, 8990086.848588351, 8992370.3125, 9023863.620321712, 9033924.805646012, 9092160.653829597, 9095844.000911955, 9163353.494224448, 9188182.137801256, 9258283.678236874, 9286639.502158325, 9295834.375, 9299207.407273667, 9343423.261276508, 9360266.18421911, 9410804.47186099, 9410820.22457622, 9432051.879411984, 9473175.195609378, 9476643.75, 9508075.655324738, 9511889.708365962, 9613874.511384714, 9625040.558625018, 9625164.20291628, 9638977.373930046, 9666140.045569753, 9671676.866775503, 9673257.8125, 9679532.717368454, 9680727.452920264, 9681359.375, 9705052.956186488, 9707407.226047626, 9709420.702200953, 9939953.835931828, 9947037.500886794, 9947582.143217256, 9947616.254266793, 9947713.037218357, 10098756.722034426, 10114392.1875, 10140178.54214567, 10168673.030256838, 10245400.45599091, 10250363.579958308, 10250829.060798533, 10341917.016452713, 10341929.948088102, 10421804.62438267, 10456844.32915452, 10476242.949504135, 10493001.642198984, 10575321.064422602, 10598941.13536112, 10669774.496317685, 10734832.8125, 10794357.8125, 10824920.083542883, 10835825.0, 10856231.25, 10860109.301672088, 10860204.501909243, 10931707.8125, 10951379.102108477, 10987262.5, 10994413.248427745, 10999055.153642982, 11006610.9375, 11037250.353905322, 11096368.75, 11128470.59811617, 11212479.359960804, 11212480.35541019, 11268659.365610491, 11382613.633514645, 11412055.952972887, 11485456.690445792, 11510798.231705148, 11679995.96798333, 11758745.3125, 11759409.285680968, 11786719.579370689, 11907971.875, 11947923.474890647, 12702501.55572402, 12703397.614465257, 12770825.85181278, 12770833.96459695, 12903378.91299109, 13161077.118573125, 13362243.018636232, 13812197.390679702, 14130356.279429603, 15073546.476703227, 15703447.85121397, 15967063.132790694, 15975445.285061106, 15983375.0, 16000438.781963458, 16003400.546629013, 16004225.778088594, 16006234.375, 16012029.6875, 16013205.188377704, 16014734.181089526, 16019921.52185321, 16020284.375, 16020746.650079168, 16021970.430315385, 16024862.201077795, 16027269.944995206, 16027942.69561429, 16028143.628686234, 16028153.763264634, 16029173.688740635, 16029797.669739647, 16030219.688424889, 16033219.394997867, 16033374.035032988, 16033861.266319076, 16034277.264519773, 16044696.792563789, 16047613.596988507, 16049559.375, 16053158.701715937, 16057323.4375, 16061267.1875, 16076061.924175762, 16089310.384396724, 16091113.047003359, 16094590.01395407, 16101693.746823965, 16102161.248012103, 16103727.077608466, 16105836.81197619, 16111206.882711843, 16117304.795089623, 16123011.420360623, 16123178.125, 16124209.062846482, 16126799.98532281, 16131457.48823512, 16159658.292866107, 16159670.459857572, 16170748.37696521, 16178816.696698017, 16180820.151172316, 16181253.198424583, 16182625.0, 16182793.581977692, 16185856.441664541, 16188151.471513303, 16190857.8125, 16190873.534141613, 16192794.857072463, 16193136.99130448, 16193942.820979508, 16194149.577890957, 16194628.672531167, 16196087.070173824, 16197832.8125, 16197846.298173515, 16198335.408474391, 16198353.739516448, 16199219.238548953, 16201475.934689812, 16202189.939807765, 16202498.294068523, 16202552.388944034, 16202765.796936864, 16202922.949627612, 16203028.125, 16203054.361786857, 16203276.478806395, 16203492.1875, 16203548.256782122, 16203582.047625408, 16203622.925476385, 16203639.447967488, 16203903.738182325, 16204134.375, 16204141.513762644, 16204304.563938085, 16204310.049756793, 16204352.060637334, 16204420.3125, 16204462.208588608, 16204814.957411963, 16204821.875, 16204829.436457071, 16204851.5625, 16205209.81185865, 16205223.328967564, 16205364.56513149, 16205386.99308477, 16205494.877499951, 16205549.169438977, 16205585.395114595, 16205632.672856508, 16205645.111212235, 16205649.58967072, 16205675.28440221, 16205694.04788572, 16205695.403249042, 16205731.17083987, 16205760.587224947, 16205779.095125193, 16205788.597543662, 16205902.998795146, 16205915.12934558, 16205945.745684996, 16206019.71295666, 16206026.413997134, 16206113.145155473, 16206173.34115582, 16206324.664453605, 16206366.802089235, 16206410.614937702, 16206432.408029495, 16206436.292292394, 16206457.124816453, 16206571.200812276, 16206602.660124427, 16206608.82308006, 16206619.632093782, 16206720.3125, 16206721.556366265, 16206736.821106702, 16206744.747951135, 16206854.485657817, 16206862.298999459, 16207050.144014046, 16207055.151596809, 16207058.075438708, 16207062.734065278, 16207088.636386957, 16207097.452865582, 16207107.455154173, 16207110.662589548, 16207120.465474168, 16207151.122075, 16207154.6384436, 16207156.25, 16207159.54104702, 16207176.5625, 16207181.25, 16207352.933212722, 16207372.764958574, 16207382.96349704, 16207385.53057189, 16207415.03539173, 16207417.088251065, 16207423.4375, 16207429.6875, 16207440.368919, 16207443.298272546, 16207446.291629799, 16207461.812393589, 16207483.534085589, 16207530.825421393, 16207611.600039585, 16207624.50085341, 16207668.75, 16207681.81290489, 16207744.172887065, 16207803.02782234, 16207812.5, 16207839.201977473, 16207923.373718219, 16207930.980923822, 16207934.99277117, 16207943.75, 16208186.662612237, 16208288.780427814, 16208329.4307924, 16208353.125, 16208365.414823158, 16208368.18666921, 16208399.749228103, 16208425.08292508, 16208432.206033265, 16208434.702143043, 16208459.357214125, 16208482.308883851, 16208561.11001401, 16208565.139014997, 16208575.0, 16208585.86209186, 16208662.142132761, 16208669.510600239, 16208681.997399189, 16208691.23561566, 16208706.255704407, 16208709.102339998, 16208746.991572937, 16208828.849811837, 16208888.42539675, 16208902.028963516, 16208910.032445226, 16208921.434395267, 16209000.10386588, 16209008.928752225, 16209010.9375, 16209140.625, 16209171.856984854, 16209174.444345938, 16209185.917648746, 16209228.891844293, 16209236.442323564, 16209240.625, 16209264.0625, 16209272.274624782, 16209289.7936305, 16209346.840639606, 16209370.224296663, 16209404.861381564, 16209420.3125, 16209424.626553623, 16209490.304140162, 16209519.802903652, 16209544.757913612, 16209607.8125, 16209611.184333507, 16209623.603543133, 16209638.075376598, 16209649.329802962, 16209660.9375, 16209670.62542617, 16209680.683055287, 16209746.403555525, 16209837.552143583, 16209845.643559169, 16209862.37236405, 16209874.260214683, 16209875.0, 16209886.866938068, 16210000.013060093, 16210028.742586, 16210082.079162626, 16210198.688931327, 16210348.134684464, 16210400.726020444, 16210431.25, 16210441.09427573, 16210481.986191746, 16210529.850021752, 16210629.6875, 16210653.125, 16210709.375, 16210782.145125387, 16210928.187477548, 16211014.049246574, 16211250.0, 16211253.35827926, 16211274.098878387, 16211313.474744432, 16211343.812054962, 16211432.8125, 16211556.478196135, 16211939.0625, 16212947.216962518, 16215405.194514401, 16227408.388777623, 16228032.466202414, 16228560.9375, 16228776.138028476, 16237829.56004968, 16241242.1875, 16270035.46062934, 519606755.8849367], [20.075061542102226, 66.77501880346156, 71.30921418376734, 33.81409237767746, 7.074857048482096, 11.70050516397771, 15.034745588959952, 29.53884074984992, 37.839303192928384, 12.193867000567472, 31.521251638669202, 46.03448496680552, 13.590501857318467, 29.004085638048977, 27.964998497182975, 9.77887430555006, 49.72566040428604, 97.72299844866785, 31.545103971469203, 5.414683372178093, 11.769160038570622, 26.153907584027376, 47.03098597585904, 61.75344767832768, 105.92749539690658, 22.15326427365694, 71.73132744910212, 7.75733595205109, 5.839983218473605, 53.975902026678945, 19.50700929876801, 52.06737274794605, 22.865349671609543, 51.69173578569769, 7.93537467089529, 16.619537116783736, 45.10979161006365, 13.706026387527274, 5.035873157674501, 15.732565242874298, 16.424995062909233, 14.996913058459718, 53.026874073126315, 7.57072436982771, 9.075326368024315, 21.366132287659628, 6.2861479236927185, 19.730165168515622, 25.171883864739012, 50.09488959457041, 33.40238010010333, 24.26566169758266, 6.261663468686628, 72.56414211260218, 46.89414481711042, 60.36833824809907, 13.49815650227515, 7.101272679523748, 14.063069587278918, 64.33326313971097, 17.54564790495636, 17.400875238925167, 33.81211326632837, 25.641419217359164, 45.959500570241744, 32.60037403722289, 6.338417172140698, 5.513321849800691, 144.77660776556993, 6.09087947205094, 15.418885785530962, 25.341659901650218, 35.50541418600724, 71.27605716641206, 39.36457352126293, 22.356707018424903, 65.42515335689728, 68.2575927594729, 24.133439268042913, 47.76671146708033, 33.40476111405811, 126.96022669926693, 28.28381792259089, 74.6627860206372, 25.777967968611875, 40.406298006575696, 27.580978772542707, 17.60683442732823, 24.149425848278767, 29.955654198862227, 19.937313845276584, 7.245954106744644, 22.12695387466793, 15.390349713919543, 61.72575615106754, 13.837186854878645, 90.43150041250664, 19.72717557096077, 102.87459481374802, 10.233225902257045, 5.820793588919748, 29.095434121483848, 87.87406886971763, 8.550509651572053, 83.28020402227706, 39.59866491728035, 69.78283933231832, 32.517980802263274, 49.67115507706114, 38.368564654080096, 6.046106477957383, 18.39865703019591, 36.12005183428919, 175.4001918731995, 5.207443454605349, 31.64684812583591, 31.19065878790927, 22.562619910795718, 17.753333098417098, 64.35348765113199, 53.65973496730773, 32.259431565601645, 6.22694821099251, 10.955711877036798, 11.218193327512111, 106.67174755282105, 11.494296869717177, 72.70579783879532, 17.660134385819255, 18.309139839309328, 40.16919320368989, 64.21528181476219, 34.92093284491044, 45.82482037527302, 17.625145695633268, 5.85627194580403, 31.890368612354873, 11.782869552397484, 17.182824492849548, 35.08887946940321, 26.87711449273235, 63.38896079721869, 60.07059578532748, 42.165243682810924, 80.25058160964015, 38.874103814829205, 5.6416442625130125, 36.82559897347932, 15.36206204366398, 32.96566193837316, 28.06836206156396, 124.79263185035319, 61.12221867682287, 55.398488785627364, 27.36035219431543, 19.49389610603575, 100.58937753809113, 57.60619097843332, 55.93203579373673, 32.523086258791146, 18.624159832796064, 18.353005213084543, 37.94032440378395, 21.27871339131735, 78.75600929069337, 98.59560620718261, 84.97863253178672, 25.842246402921493, 39.33559768338594, 22.089969157499866, 61.39854032849083, 73.23072089807604, 15.992659867152057, 48.40359176537403, 132.7818624073178, 19.731786470195498, 5.735322272878144, 64.09168706998275, 94.6211113483372, 153.28901569714634, 9.853697188911143, 24.891663336309133, 7.571743723751527, 19.33341907526585, 15.743460218600939, 104.48986857739114, 18.208845155676997, 6.899683133084568, 21.60742561754838, 56.75489258231201, 6.791867150028416, 154.6530107752742, 21.49315371002927, 17.041341153779733, 82.43660062200432, 32.7731206655763, 49.141084537808126, 81.92493486169869, 68.91706542036114, 8.880306959568047, 61.38512717000971, 17.137335378912173, 117.29381605237185, 75.63676544934401, 18.040954334545884, 5.644472144757256, 15.570333763844427, 40.8827060879917, 9.710922911234558, 27.66903185465518, 28.530312013073367, 62.17721225319873, 66.97693539459904, 52.66803061030485, 40.440967609082264, 6.750970765136743, 17.72924435091622, 18.94133570498829, 80.11074487024763, 200.4687435638146, 9.803750163539398, 27.86735799574362, 36.03563045323571, 111.40771153564074, 14.135183274055848, 28.705418623390603, 49.86140889434661, 27.153670018917143, 25.3311248821529, 44.19954400207433, 13.098108576588519, 88.86851163005997, 28.662640324088002, 72.7412717629813, 33.08871361322505, 8.48303612720662, 134.53924851447078, 34.726314967854066, 13.82456864405277, 5.999761167384706, 14.90584138162558, 10.517429905455879, 84.07107132390601, 5.554845547520138, 46.29439608249341, 32.29971596611013, 46.244324574568296, 175.48145717355803, 51.9993133909239, 74.11136151973685, 84.24985240924052, 35.28629817457856, 66.82152995479836, 73.41750955849491, 15.380760736444621, 6.372167792201014, 23.103004882660066, 50.34380366030459, 22.609226339947014, 15.054880423896535, 16.243776038191832, 64.02658086600162, 5.650037852453839, 37.13493187479041, 58.45772323423543, 39.80891410031519, 6.8777746120311525, 12.609221508719802, 36.937759811044494, 29.134428388594007, 17.269744148317944, 14.25056688889304, 99.88607411522486, 30.778874450872873, 59.417979753494855, 7.6464126573571995, 15.618594619492722, 13.704182926784632, 5.418260638441956, 31.29959852148355, 34.161478163168084, 8.424061696299587, 29.54350733939834, 20.338667587948287, 12.311634057422053, 7.834665144646622, 52.561913777985836, 74.36932901112435, 26.9940876502495, 72.67222274453808, 37.157181032647024, 19.922720651328536, 54.473683572012504, 55.0145931336484, 44.45171381886674, 32.222463101184005, 17.53661635240796, 27.22952767156285, 34.72988744855627, 46.150935049855974, 6.14335567859943, 19.80230646197144, 12.285238657815746, 44.86213354232169, 29.948216647436638, 81.73477797450509, 15.561876071655293, 7.597937806432295, 13.834568928806593, 9.773532103381003, 21.43456418034761, 11.034387205159277, 21.58534871972576, 26.628846719737574, 64.62822227149175, 37.675655258466946, 25.596882135503552, 12.701804409578568, 80.94868902357908, 57.38661800315644, 15.323116509784889, 17.136891162234, 37.646783322989066, 39.13811489760208, 130.23304259792008, 27.929973354612947, 22.2228984705647, 62.898221684290874, 21.40617996414875, 5.6738967712777795, 9.91133951037937, 46.92876369831072, 12.989123212956686, 33.67937179737288, 10.910678631548743, 19.422507816416054, 6.289168527409597, 19.420527062482538, 56.731348456075985, 15.816455044057177, 14.631823894047551, 38.86662298407216, 13.815589868011408, 19.20345732481752, 26.591824484836582, 9.680090161869535, 44.11647870041121, 49.02297210020268, 50.73596560627708, 38.32281255368909, 6.3563869225383405, 17.25422819848931, 47.02723714357127, 130.56758244432976, 8.824234517119574, 64.21033530136208, 6.71440178167111, 52.37628249537762, 11.886876112065666, 113.5208912850465, 21.414626888600296, 6.66635568855736, 13.50766208408557, 12.5035490192224, 21.080452662489915, 33.95416016943703, 14.226053488610662, 14.24286990115944, 30.258168724214315, 5.547714671452123, 5.247692290593898, 64.9643512407398, 21.162318443445116, 20.618582634751206, 30.05812513821806, 117.61713861795866, 8.040776772372965, 134.31850932376992, 76.03539171963423, 14.40416549157068, 12.311973797609056, 25.113596484836293, 13.80020186926111, 13.651231320402653, 13.839777217641448, 15.226517806220984, 25.135747668609483, 8.810559784443628, 74.9187770416099, 7.443330189210997, 37.476861465854704, 6.95338673735558, 101.67118240654064, 77.40504268766605, 19.876676478015327, 6.402090871412261, 25.09113810038808, 55.46435473173046, 17.22753825294696, 5.215409166688883, 12.470010538079146, 26.696496578698653, 10.521151123570466, 7.282117998295218, 57.67719575072533, 36.54931471150268, 16.768363719918703, 8.066640717072982, 9.290264843750602, 89.3349235119259, 15.614569824684533, 97.65430615834937, 45.21188735764461, 9.174335520454225, 10.217388888177933, 9.462428806807115, 17.918982783953684, 18.16725921679945, 5.806673010592574, 51.62910037636868, 36.66838772817987, 5.238034194377177, 45.85564508243577, 78.48594679646877, 6.184064728162417, 24.591855045419287, 22.279975467008267, 6.596219987437982, 5.777948306581928, 33.757052308011886, 12.138520191598102, 9.716549205716694, 38.411936037576105, 32.17676677481561, 11.352949711372947, 7.5803306584315475, 6.4457225145365316, 12.253374468235686, 41.587137707384564, 23.374457897237974, 56.052981970018116, 59.096403009494, 58.58714815112546, 18.182677600021197, 13.155881101199387, 6.448486164988477, 51.03462439844092, 5.66921817148592, 9.624237885748267, 6.086251593232763, 22.591174941851676, 9.481661519838582, 48.09399364410018, 65.79830953190535, 20.891795440486142, 13.732933120846967, 7.903629232348305, 41.265466380423526, 79.20365857734731, 26.426247906280242, 15.568340752967961, 6.7870047665510675, 24.399169476286723, 12.965607821430828, 6.562787057628997, 23.558294761021227, 42.59169759231611, 5.294306419879466, 8.980577149229907, 32.8987019549005, 56.66240458079238, 28.17321310933923, 35.90856066221975, 68.42235815537283, 6.92690785216354, 62.19593855137672, 58.01074279079167, 71.04271192827615, 16.72015310948123, 7.890965366385145, 9.302517535880357, 5.9317608789624225, 5.101893962485382, 73.1577357838217, 5.90753024694414, 7.4911691974029635, 21.3979027250536, 26.380956496829693, 74.82327331059268, 22.771014108389085, 39.826113003446906, 22.205780406768014, 19.94817751189301, 42.15099726849792, 8.459293281230426, 13.146364855603132, 92.7986506850038, 5.099483888018506, 24.831513003166272, 72.20748546458901, 12.714039006344375, 50.75934390203771, 56.254828465852796, 7.607968136013016, 12.252981159405593, 5.115132740620922, 15.94402985077545, 13.90060361165151, 58.704611503952044, 41.84505547880708, 111.12550442527078, 10.042515331068806, 6.287952529110162, 64.8233749777597, 29.650362774678168, 96.35057326621563, 29.243951489237286, 19.548042522869395, 30.881475707153676, 32.622616087897434, 9.536763523497997, 25.323146319989455, 10.318476100292251, 6.496863087005234, 49.018398368786585, 5.544259074573224, 90.1785461420169, 32.037353455483014, 11.217234774533058, 28.110425575302507, 58.906534948122726, 34.30405590000561, 17.533230097891433, 41.59450720529936, 66.04597200938152, 73.98664303404462, 22.484793953950756])
caption, T_hit,count1 = 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, count1, caption);
([7341949.345040617, 7349617.1875, 7383835.9375, 7393606.25, 7397523.257230354, 7445434.712954623, 7468078.458700638, 7590140.625, 7652841.934960744, 7652864.652044484, 7661440.2490893705, 7687321.802820627, 7690312.09380085, 7692894.330288865, 7745943.105911829, 7750797.360693115, 7756640.625, 7782889.0625, 7801279.2775708325, 7819941.042632043, 7822439.6355259335, 7836283.059348688, 7840826.5625, 7841531.25, 7844158.80428424, 7846588.273996023, 7848017.1875, 7852425.4068670105, 7852495.657854303, 7853495.3125, 7854843.397894762, 7863713.753118767, 7890918.37737375, 7898289.484058416, 7906446.146062636, 7942191.005106731, 7942281.25, 7942853.863910955, 7968112.160764721, 7969794.821085451, 7979844.091446953, 7985475.01446025, 7987010.9375, 8034182.721701814, 8036230.930784949, 8041966.674677169, 8042772.023527585, 8044907.176177609, 8045836.646008059, 8047148.4375, 8051003.238079173, 8065501.418240794, 8066104.64356547, 8101684.375, 8112638.348822645, 8114563.742646959, 8116540.651892431, 8117282.735175535, 8130018.498716827, 8134465.59919873, 8138842.736031749, 8139909.682297353, 8139931.322275146, 8140210.622996046, 8149181.222058955, 8150476.5625, 8152705.231142673, 8156530.836047882, 8186527.871475717, 8208695.339193121, 8213703.3227339415, 8226792.077911775, 8227906.596938308, 8237363.329474506, 8243156.25, 8249196.369103571, 8266282.997052491, 8272882.8125, 8287762.229380851, 8290214.681895636, 8291995.3125, 8305216.716505616, 8323184.362702257, 8324557.8125, 8327392.430370559, 8333623.874711879, 8345696.293689768, 8356584.044487525, 8356648.5906556435, 8356681.6637081485, 8375821.3815600155, 8382733.265861373, 8385429.248776955, 8394089.175476993, 8397533.986457879, 8399313.0454553, 8407952.949725354, 8425510.240162669, 8427825.849235635, 8442885.27843374, 8461856.341659928, 8464745.3125, 8470284.375, 8471206.71784104, 8471648.4375, 8513756.25, 8525624.588513607, 8540905.44721093, 8541380.407720702, 8542445.3125, 8545526.148771456, 8545727.246535618, 8551882.676993605, 8560948.013776204, 8561000.15712115, 8565946.875, 8578140.625, 8580568.63134677, 8580911.433435109, 8582014.157119911, 8600857.8125, 8603903.125, 8603962.174638461, 8619068.637103815, 8622017.244629288, 8626520.3125, 8628428.347521422, 8630178.907587143, 8632025.661892401, 8638678.79559247, 8640939.014642905, 8642656.25, 8661514.0625, 8661572.331479887, 8730926.784285095, 8731486.49317093, 8743176.011245739, 8751563.121227873, 8751636.42294095, 8788279.591452878, 8791042.377614316, 8794171.895300547, 8819184.309256323, 8821583.071259731, 8822595.604458127, 8838859.47066396, 8849970.800846713, 8852687.791912725, 8852709.481128328, 8853015.565010209, 8855471.07453266, 8857821.875, 8860770.289181355, 8870293.75, 8870311.858334515, 8876267.774277605, 8885243.295197597, 8888593.75, 8890747.79732472, 8917227.979842346, 8920241.794157727, 8930942.86470136, 8940040.555770405, 8959222.915736785, 8969429.6875, 8972732.8125, 8972741.79836481, 8990086.848588351, 8992370.3125, 9023863.620321712, 9033924.805646012, 9092160.653829597, 9095844.000911955, 9163353.494224448, 9188182.137801256, 9258283.678236874, 9286639.502158325, 9295834.375, 9299207.407273667, 9343423.261276508, 9360266.18421911, 9410804.47186099, 9410820.22457622, 9432051.879411984, 9473175.195609378, 9476643.75, 9508075.655324738, 9511889.708365962, 9613874.511384714, 9625040.558625018, 9625164.20291628, 9638977.373930046, 9666140.045569753, 9671676.866775503, 9673257.8125, 9679532.717368454, 9680727.452920264, 9681359.375, 9705052.956186488, 9707407.226047626, 9709420.702200953, 9939953.835931828, 9947037.500886794, 9947582.143217256, 9947616.254266793, 9947713.037218357, 10098756.722034426, 10114392.1875, 10140178.54214567, 10168673.030256838, 10245400.45599091, 10250363.579958308, 10250829.060798533, 10341917.016452713, 10341929.948088102, 10421804.62438267, 10456844.32915452, 10476242.949504135, 10493001.642198984, 10575321.064422602, 10598941.13536112, 10669774.496317685, 10734832.8125, 10794357.8125, 10824920.083542883, 10835825.0, 10856231.25, 10860109.301672088, 10860204.501909243, 10931707.8125, 10951379.102108477, 10987262.5, 10994413.248427745, 10999055.153642982, 11006610.9375, 11037250.353905322, 11096368.75, 11128470.59811617, 11212479.359960804, 11212480.35541019, 11268659.365610491, 11382613.633514645, 11412055.952972887, 11485456.690445792, 11510798.231705148, 11679995.96798333, 11758745.3125, 11759409.285680968, 11786719.579370689, 11907971.875, 11947923.474890647, 12702501.55572402, 12703397.614465257, 12770825.85181278, 12770833.96459695, 12903378.91299109, 13161077.118573125, 13362243.018636232, 13812197.390679702, 14130356.279429603, 15073546.476703227, 15703447.85121397, 15967063.132790694, 15975445.285061106, 15983375.0, 16000438.781963458, 16003400.546629013, 16004225.778088594, 16006234.375, 16012029.6875, 16013205.188377704, 16014734.181089526, 16019921.52185321, 16020284.375, 16020746.650079168, 16021970.430315385, 16024862.201077795, 16027269.944995206, 16027942.69561429, 16028143.628686234, 16028153.763264634, 16029173.688740635, 16029797.669739647, 16030219.688424889, 16033219.394997867, 16033374.035032988, 16033861.266319076, 16034277.264519773, 16044696.792563789, 16047613.596988507, 16049559.375, 16053158.701715937, 16057323.4375, 16061267.1875, 16076061.924175762, 16089310.384396724, 16091113.047003359, 16094590.01395407, 16101693.746823965, 16102161.248012103, 16103727.077608466, 16105836.81197619, 16111206.882711843, 16117304.795089623, 16123011.420360623, 16123178.125, 16124209.062846482, 16126799.98532281, 16131457.48823512, 16159658.292866107, 16159670.459857572, 16170748.37696521, 16178816.696698017, 16180820.151172316, 16181253.198424583, 16182625.0, 16182793.581977692, 16185856.441664541, 16188151.471513303, 16190857.8125, 16190873.534141613, 16192794.857072463, 16193136.99130448, 16193942.820979508, 16194149.577890957, 16194628.672531167, 16196087.070173824, 16197832.8125, 16197846.298173515, 16198335.408474391, 16198353.739516448, 16199219.238548953, 16201475.934689812, 16202189.939807765, 16202498.294068523, 16202552.388944034, 16202765.796936864, 16202922.949627612, 16203028.125, 16203054.361786857, 16203276.478806395, 16203492.1875, 16203548.256782122, 16203582.047625408, 16203622.925476385, 16203639.447967488, 16203903.738182325, 16204134.375, 16204141.513762644, 16204304.563938085, 16204310.049756793, 16204352.060637334, 16204420.3125, 16204462.208588608, 16204814.957411963, 16204821.875, 16204829.436457071, 16204851.5625, 16205209.81185865, 16205223.328967564, 16205364.56513149, 16205386.99308477, 16205494.877499951, 16205549.169438977, 16205585.395114595, 16205632.672856508, 16205645.111212235, 16205649.58967072, 16205675.28440221, 16205694.04788572, 16205695.403249042, 16205731.17083987, 16205760.587224947, 16205779.095125193, 16205788.597543662, 16205902.998795146, 16205915.12934558, 16205945.745684996, 16206019.71295666, 16206026.413997134, 16206113.145155473, 16206173.34115582, 16206324.664453605, 16206366.802089235, 16206410.614937702, 16206432.408029495, 16206436.292292394, 16206457.124816453, 16206571.200812276, 16206602.660124427, 16206608.82308006, 16206619.632093782, 16206720.3125, 16206721.556366265, 16206736.821106702, 16206744.747951135, 16206854.485657817, 16206862.298999459, 16207050.144014046, 16207055.151596809, 16207058.075438708, 16207062.734065278, 16207088.636386957, 16207097.452865582, 16207107.455154173, 16207110.662589548, 16207120.465474168, 16207151.122075, 16207154.6384436, 16207156.25, 16207159.54104702, 16207176.5625, 16207181.25, 16207352.933212722, 16207372.764958574, 16207382.96349704, 16207385.53057189, 16207415.03539173, 16207417.088251065, 16207423.4375, 16207429.6875, 16207440.368919, 16207443.298272546, 16207446.291629799, 16207461.812393589, 16207483.534085589, 16207530.825421393, 16207611.600039585, 16207624.50085341, 16207668.75, 16207681.81290489, 16207744.172887065, 16207803.02782234, 16207812.5, 16207839.201977473, 16207923.373718219, 16207930.980923822, 16207934.99277117, 16207943.75, 16208186.662612237, 16208288.780427814, 16208329.4307924, 16208353.125, 16208365.414823158, 16208368.18666921, 16208399.749228103, 16208425.08292508, 16208432.206033265, 16208434.702143043, 16208459.357214125, 16208482.308883851, 16208561.11001401, 16208565.139014997, 16208575.0, 16208585.86209186, 16208662.142132761, 16208669.510600239, 16208681.997399189, 16208691.23561566, 16208706.255704407, 16208709.102339998, 16208746.991572937, 16208828.849811837, 16208888.42539675, 16208902.028963516, 16208910.032445226, 16208921.434395267, 16209000.10386588, 16209008.928752225, 16209010.9375, 16209140.625, 16209171.856984854, 16209174.444345938, 16209185.917648746, 16209228.891844293, 16209236.442323564, 16209240.625, 16209264.0625, 16209272.274624782, 16209289.7936305, 16209346.840639606, 16209370.224296663, 16209404.861381564, 16209420.3125, 16209424.626553623, 16209490.304140162, 16209519.802903652, 16209544.757913612, 16209607.8125, 16209611.184333507, 16209623.603543133, 16209638.075376598, 16209649.329802962, 16209660.9375, 16209670.62542617, 16209680.683055287, 16209746.403555525, 16209837.552143583, 16209845.643559169, 16209862.37236405, 16209874.260214683, 16209875.0, 16209886.866938068, 16210000.013060093, 16210028.742586, 16210082.079162626, 16210198.688931327, 16210348.134684464, 16210400.726020444, 16210431.25, 16210441.09427573, 16210481.986191746, 16210529.850021752, 16210629.6875, 16210653.125, 16210709.375, 16210782.145125387, 16210928.187477548, 16211014.049246574, 16211250.0, 16211253.35827926, 16211274.098878387, 16211313.474744432, 16211343.812054962, 16211432.8125, 16211556.478196135, 16211939.0625, 16212947.216962518, 16215405.194514401, 16227408.388777623, 16228032.466202414, 16228560.9375, 16228776.138028476, 16237829.56004968, 16241242.1875, 16270035.46062934, 519606755.8849367], [20.075061542102226, 66.77501880346156, 71.30921418376734, 33.81409237767746, 7.074857048482096, 11.70050516397771, 15.034745588959952, 29.53884074984992, 37.839303192928384, 12.193867000567472, 31.521251638669202, 46.03448496680552, 13.590501857318467, 29.004085638048977, 27.964998497182975, 9.77887430555006, 49.72566040428604, 97.72299844866785, 31.545103971469203, 5.414683372178093, 11.769160038570622, 26.153907584027376, 47.03098597585904, 61.75344767832768, 105.92749539690658, 22.15326427365694, 71.73132744910212, 7.75733595205109, 5.839983218473605, 53.975902026678945, 19.50700929876801, 52.06737274794605, 22.865349671609543, 51.69173578569769, 7.93537467089529, 16.619537116783736, 45.10979161006365, 13.706026387527274, 5.035873157674501, 15.732565242874298, 16.424995062909233, 14.996913058459718, 53.026874073126315, 7.57072436982771, 9.075326368024315, 21.366132287659628, 6.2861479236927185, 19.730165168515622, 25.171883864739012, 50.09488959457041, 33.40238010010333, 24.26566169758266, 6.261663468686628, 72.56414211260218, 46.89414481711042, 60.36833824809907, 13.49815650227515, 7.101272679523748, 14.063069587278918, 64.33326313971097, 17.54564790495636, 17.400875238925167, 33.81211326632837, 25.641419217359164, 45.959500570241744, 32.60037403722289, 6.338417172140698, 5.513321849800691, 144.77660776556993, 6.09087947205094, 15.418885785530962, 25.341659901650218, 35.50541418600724, 71.27605716641206, 39.36457352126293, 22.356707018424903, 65.42515335689728, 68.2575927594729, 24.133439268042913, 47.76671146708033, 33.40476111405811, 126.96022669926693, 28.28381792259089, 74.6627860206372, 25.777967968611875, 40.406298006575696, 27.580978772542707, 17.60683442732823, 24.149425848278767, 29.955654198862227, 19.937313845276584, 7.245954106744644, 22.12695387466793, 15.390349713919543, 61.72575615106754, 13.837186854878645, 90.43150041250664, 19.72717557096077, 102.87459481374802, 10.233225902257045, 5.820793588919748, 29.095434121483848, 87.87406886971763, 8.550509651572053, 83.28020402227706, 39.59866491728035, 69.78283933231832, 32.517980802263274, 49.67115507706114, 38.368564654080096, 6.046106477957383, 18.39865703019591, 36.12005183428919, 175.4001918731995, 5.207443454605349, 31.64684812583591, 31.19065878790927, 22.562619910795718, 17.753333098417098, 64.35348765113199, 53.65973496730773, 32.259431565601645, 6.22694821099251, 10.955711877036798, 11.218193327512111, 106.67174755282105, 11.494296869717177, 72.70579783879532, 17.660134385819255, 18.309139839309328, 40.16919320368989, 64.21528181476219, 34.92093284491044, 45.82482037527302, 17.625145695633268, 5.85627194580403, 31.890368612354873, 11.782869552397484, 17.182824492849548, 35.08887946940321, 26.87711449273235, 63.38896079721869, 60.07059578532748, 42.165243682810924, 80.25058160964015, 38.874103814829205, 5.6416442625130125, 36.82559897347932, 15.36206204366398, 32.96566193837316, 28.06836206156396, 124.79263185035319, 61.12221867682287, 55.398488785627364, 27.36035219431543, 19.49389610603575, 100.58937753809113, 57.60619097843332, 55.93203579373673, 32.523086258791146, 18.624159832796064, 18.353005213084543, 37.94032440378395, 21.27871339131735, 78.75600929069337, 98.59560620718261, 84.97863253178672, 25.842246402921493, 39.33559768338594, 22.089969157499866, 61.39854032849083, 73.23072089807604, 15.992659867152057, 48.40359176537403, 132.7818624073178, 19.731786470195498, 5.735322272878144, 64.09168706998275, 94.6211113483372, 153.28901569714634, 9.853697188911143, 24.891663336309133, 7.571743723751527, 19.33341907526585, 15.743460218600939, 104.48986857739114, 18.208845155676997, 6.899683133084568, 21.60742561754838, 56.75489258231201, 6.791867150028416, 154.6530107752742, 21.49315371002927, 17.041341153779733, 82.43660062200432, 32.7731206655763, 49.141084537808126, 81.92493486169869, 68.91706542036114, 8.880306959568047, 61.38512717000971, 17.137335378912173, 117.29381605237185, 75.63676544934401, 18.040954334545884, 5.644472144757256, 15.570333763844427, 40.8827060879917, 9.710922911234558, 27.66903185465518, 28.530312013073367, 62.17721225319873, 66.97693539459904, 52.66803061030485, 40.440967609082264, 6.750970765136743, 17.72924435091622, 18.94133570498829, 80.11074487024763, 200.4687435638146, 9.803750163539398, 27.86735799574362, 36.03563045323571, 111.40771153564074, 14.135183274055848, 28.705418623390603, 49.86140889434661, 27.153670018917143, 25.3311248821529, 44.19954400207433, 13.098108576588519, 88.86851163005997, 28.662640324088002, 72.7412717629813, 33.08871361322505, 8.48303612720662, 134.53924851447078, 34.726314967854066, 13.82456864405277, 5.999761167384706, 14.90584138162558, 10.517429905455879, 84.07107132390601, 5.554845547520138, 46.29439608249341, 32.29971596611013, 46.244324574568296, 175.48145717355803, 51.9993133909239, 74.11136151973685, 84.24985240924052, 35.28629817457856, 66.82152995479836, 73.41750955849491, 15.380760736444621, 6.372167792201014, 23.103004882660066, 50.34380366030459, 22.609226339947014, 15.054880423896535, 16.243776038191832, 64.02658086600162, 5.650037852453839, 37.13493187479041, 58.45772323423543, 39.80891410031519, 6.8777746120311525, 12.609221508719802, 36.937759811044494, 29.134428388594007, 17.269744148317944, 14.25056688889304, 99.88607411522486, 30.778874450872873, 59.417979753494855, 7.6464126573571995, 15.618594619492722, 13.704182926784632, 5.418260638441956, 31.29959852148355, 34.161478163168084, 8.424061696299587, 29.54350733939834, 20.338667587948287, 12.311634057422053, 7.834665144646622, 52.561913777985836, 74.36932901112435, 26.9940876502495, 72.67222274453808, 37.157181032647024, 19.922720651328536, 54.473683572012504, 55.0145931336484, 44.45171381886674, 32.222463101184005, 17.53661635240796, 27.22952767156285, 34.72988744855627, 46.150935049855974, 6.14335567859943, 19.80230646197144, 12.285238657815746, 44.86213354232169, 29.948216647436638, 81.73477797450509, 15.561876071655293, 7.597937806432295, 13.834568928806593, 9.773532103381003, 21.43456418034761, 11.034387205159277, 21.58534871972576, 26.628846719737574, 64.62822227149175, 37.675655258466946, 25.596882135503552, 12.701804409578568, 80.94868902357908, 57.38661800315644, 15.323116509784889, 17.136891162234, 37.646783322989066, 39.13811489760208, 130.23304259792008, 27.929973354612947, 22.2228984705647, 62.898221684290874, 21.40617996414875, 5.6738967712777795, 9.91133951037937, 46.92876369831072, 12.989123212956686, 33.67937179737288, 10.910678631548743, 19.422507816416054, 6.289168527409597, 19.420527062482538, 56.731348456075985, 15.816455044057177, 14.631823894047551, 38.86662298407216, 13.815589868011408, 19.20345732481752, 26.591824484836582, 9.680090161869535, 44.11647870041121, 49.02297210020268, 50.73596560627708, 38.32281255368909, 6.3563869225383405, 17.25422819848931, 47.02723714357127, 130.56758244432976, 8.824234517119574, 64.21033530136208, 6.71440178167111, 52.37628249537762, 11.886876112065666, 113.5208912850465, 21.414626888600296, 6.66635568855736, 13.50766208408557, 12.5035490192224, 21.080452662489915, 33.95416016943703, 14.226053488610662, 14.24286990115944, 30.258168724214315, 5.547714671452123, 5.247692290593898, 64.9643512407398, 21.162318443445116, 20.618582634751206, 30.05812513821806, 117.61713861795866, 8.040776772372965, 134.31850932376992, 76.03539171963423, 14.40416549157068, 12.311973797609056, 25.113596484836293, 13.80020186926111, 13.651231320402653, 13.839777217641448, 15.226517806220984, 25.135747668609483, 8.810559784443628, 74.9187770416099, 7.443330189210997, 37.476861465854704, 6.95338673735558, 101.67118240654064, 77.40504268766605, 19.876676478015327, 6.402090871412261, 25.09113810038808, 55.46435473173046, 17.22753825294696, 5.215409166688883, 12.470010538079146, 26.696496578698653, 10.521151123570466, 7.282117998295218, 57.67719575072533, 36.54931471150268, 16.768363719918703, 8.066640717072982, 9.290264843750602, 89.3349235119259, 15.614569824684533, 97.65430615834937, 45.21188735764461, 9.174335520454225, 10.217388888177933, 9.462428806807115, 17.918982783953684, 18.16725921679945, 5.806673010592574, 51.62910037636868, 36.66838772817987, 5.238034194377177, 45.85564508243577, 78.48594679646877, 6.184064728162417, 24.591855045419287, 22.279975467008267, 6.596219987437982, 5.777948306581928, 33.757052308011886, 12.138520191598102, 9.716549205716694, 38.411936037576105, 32.17676677481561, 11.352949711372947, 7.5803306584315475, 6.4457225145365316, 12.253374468235686, 41.587137707384564, 23.374457897237974, 56.052981970018116, 59.096403009494, 58.58714815112546, 18.182677600021197, 13.155881101199387, 6.448486164988477, 51.03462439844092, 5.66921817148592, 9.624237885748267, 6.086251593232763, 22.591174941851676, 9.481661519838582, 48.09399364410018, 65.79830953190535, 20.891795440486142, 13.732933120846967, 7.903629232348305, 41.265466380423526, 79.20365857734731, 26.426247906280242, 15.568340752967961, 6.7870047665510675, 24.399169476286723, 12.965607821430828, 6.562787057628997, 23.558294761021227, 42.59169759231611, 5.294306419879466, 8.980577149229907, 32.8987019549005, 56.66240458079238, 28.17321310933923, 35.90856066221975, 68.42235815537283, 6.92690785216354, 62.19593855137672, 58.01074279079167, 71.04271192827615, 16.72015310948123, 7.890965366385145, 9.302517535880357, 5.9317608789624225, 5.101893962485382, 73.1577357838217, 5.90753024694414, 7.4911691974029635, 21.3979027250536, 26.380956496829693, 74.82327331059268, 22.771014108389085, 39.826113003446906, 22.205780406768014, 19.94817751189301, 42.15099726849792, 8.459293281230426, 13.146364855603132, 92.7986506850038, 5.099483888018506, 24.831513003166272, 72.20748546458901, 12.714039006344375, 50.75934390203771, 56.254828465852796, 7.607968136013016, 12.252981159405593, 5.115132740620922, 15.94402985077545, 13.90060361165151, 58.704611503952044, 41.84505547880708, 111.12550442527078, 10.042515331068806, 6.287952529110162, 64.8233749777597, 29.650362774678168, 96.35057326621563, 29.243951489237286, 19.548042522869395, 30.881475707153676, 32.622616087897434, 9.536763523497997, 25.323146319989455, 10.318476100292251, 6.496863087005234, 49.018398368786585, 5.544259074573224, 90.1785461420169, 32.037353455483014, 11.217234774533058, 28.110425575302507, 58.906534948122726, 34.30405590000561, 17.533230097891433, 41.59450720529936, 66.04597200938152, 73.98664303404462, 22.484793953950756])
Detected energies during the discharge + Energy spectrum
multiplot(icon_fig, T_int_first,E,xmean,count)
size_of_interactions_average(Etot, Ntot,size_interaction)
#energy_spectra_doublebreakdown(Etot, T_first, file_energy_spectra_doublebreakdown)