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 = 44632
shot = shot_no
identifier='H03-W0051_shot_'+str(shot)+'_450V'
detector = 'H03-W0051'
ds = np.DataSource('/tmp') # temporary storage for downloaded files
scalars_URL = 'http://golem.fjfi.cvut.cz/shots/{shot_no}/Diagnostics/PlasmaDetection/Results/{name}'
def get_scalar(shot_no, name):
return float(ds.open(scalars_URL.format(shot_no=shot_no, name=name)).read())
t_plasma_start = get_scalar(shot_no, 't_plasma_start')
t_plasma_end = get_scalar(shot_no, 't_plasma_end')
is_plasma = get_scalar(shot_no, 'b_plasma')
def get_file(shot, identifier):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/H03/{identifier}.t3pa'
url = URL.format(shot=shot, identifier=identifier)
try:
file_name_t3pa=url
with urllib.request.urlopen(file_name_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
ft3pa.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_name_t3pa
def get_file_calib(name_calib):
#Pick the discharge to analyse
URL = 'http://golem.fjfi.cvut.cz/shots/{shot}/Diagnostics/TimepixDetector/calib_matrix_H03/{name_calib}.txt'
url = URL.format(shot=shot, name_calib=name_calib)
#print(url)
try:
file_calib=url
with urllib.request.urlopen(file_calib) as calib:
line = calib.readline()
line = line.decode('utf‐8')
calib.close
except HTTPError:
print('File not found at %s. Aborting notebook execution.' % url)
raise StopExecution
return file_calib
def load_calib(file_calib):
with urllib.request.urlopen(file_calib) as fc:
calib=[] #vytvoreni 1D pole
for i in range(0,256): #tj. rozsah 0-255
temp = [] # docasne pole
for j in range(0,256):
temp.append(0) #naplneni docasneho pole 0
calib.append(temp) #naplneni pole a[] docasnym polem temp
for i in range(0,256): #nacteni calib matice do pole calib
line = fc.readline()
line = line.decode('utf‐8')
word=line.strip().split(' ')
for j in range(0,256):
calib[i][j]=float(word[j]) #i = radek, j = sloupec0
fc.close
return calib
def load_t3pa_file(file_t3pa):
index=[]
matrix_index=[]
ToA=[]
ToT=[]
FToA=[]
overflow=[]
pocet_udalosti = 0
with urllib.request.urlopen(file_t3pa) as ft3pa:
line = ft3pa.readline()
line = line.decode('utf‐8')
while True:
line = ft3pa.readline()
line = line.decode('utf‐8')
word=line.strip().split('\t') #v t3pa souboru je oddelovac \t
if line == '':
break
index.append(word[0])
matrix_index.append(word[1])
ToA.append(float(word[2]))
ToT.append(float(word[3]))
FToA.append(float(word[4]))
overflow.append(float(word[5]))
pocet_udalosti = pocet_udalosti + 1
ft3pa.close
return index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti
def noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti): #tuto fci nemus9m explicitn2 volat - volam ji v fci load_t3pa
pocet=int(0) #pocet sumicich pixelu
konst=int(len(index)/1000)+1
noise_matrix_index=[]
for i in range(0,konst):
pom = [] # pomocne pole
k=0 #pomocna promenna - udava, kolik je v czklu ve skutecnosti udalosti - aby nebyla chyba 'list index out of range'
for j in range(0,1001):
if i*1000+j>=len(index):
break
pom.append(matrix_index[i*1000+j])
k=k+1
for m in range(0,k):
count=int(0) #pocet vvyskytu stejneho matrix index behem 1000 udalosti
index_=int(-1) #budu testovat, jestli pixel na ktery koukam je sumici (abych ho nezapocital 2x)
for p in range(0,pocet):
#index=int(p)
if pom[m]==noise_matrix_index[p]:
index_=p #pixel na ktery jsem uz koukal a byl sumici
break
if index_ >=0 and pom[m]==noise_matrix_index[index_]:
continue
for l in range(0,k):
if pom[m]==pom[l]:
count=count+1
####podminka na sumici pixely
if count>=50: #kdyz se pixel vyskytne behem tisice udalosti vicekrat nez toto cislo, je sumici
noise_matrix_index.append(pom[m])
#noise_matrix_index[pocet]=pom[i]
pocet=pocet+1
pom.clear()
pocet_udalosti=len(index)
for n in range (0,pocet_udalosti):
for o in range(0,len(noise_matrix_index)):
if n >=pocet_udalosti:
break
if(matrix_index[n]==noise_matrix_index[o]):
del matrix_index[n]
del index[n]
del ToA[n]
del ToT[n]
del FToA[n]
del overflow[n]
pocet_udalosti=pocet_udalosti-1
continue
return pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow
def t3pa_data(pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow):
#rovnou vyhodim sumici pixely
pocet_udalosti,index, matrix_index, ToA, ToT, FToA, overflow=noise(index, matrix_index, ToA, ToT, FToA, overflow, pocet_udalosti)
RowNo=[]
ClmNo=[]
for i in range(0,len(matrix_index)):
#RowNo.append(int(int(matrix_index[i]))//int(256))
#ClmNo.append(int(int(matrix_index[i]))%int(256))
ClmNo.append(int(int(matrix_index[i]))//int(256)) #ver XV
RowNo.append(int(int(matrix_index[i]))%int(256)) #ver XV
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);
([5176359.375, 5278197.853278245, 5280165.688778647, 5400470.240403823, 5401085.482786234, 5415904.6875, 5421084.375, 5421088.869105486, 5426535.709674103, 5457611.751633061, 5528860.9375, 5553021.875, 5554602.984113993, 5576199.208337134, 5576243.73766927, 5577168.753350776, 5584078.125, 5589565.36641185, 5605721.994601503, 5608462.117545506, 5624286.546441099, 5642694.575049502, 5643606.008462296, 5698043.505273514, 5717114.638477524, 5733592.813033938, 5755443.75, 5796592.416548447, 6674254.6875, 6698932.8125, 6778443.830400588, 6788427.690742854, 6797823.4375, 6798121.875, 6801829.319750751, 6806044.2668590825, 6806149.070278282, 6808251.5625, 6809051.5625, 6831962.848255071, 6835660.572420684, 6865896.13556999, 6869433.299348795, 6871527.8141121445, 6881832.8125, 6897929.761116174, 6909714.138681498, 6909720.109257517, 6910810.111330104, 6914277.03586287, 6920757.8125, 6924763.951216762, 6929736.086144281, 6931215.625, 6933735.9375, 6934251.162219503, 6935482.16439149, 6944615.625, 6956656.5159969255, 6975016.22892419, 6979668.042256201, 6981323.2284298055, 7000471.694890829, 7002713.87768687, 7006943.702376271, 7011939.150588079, 7015379.649322781, 7019842.473972929, 7019972.30364368, 7022898.528069976, 7025258.801919314, 7025926.342328861, 7028357.512768078, 7029591.957604919, 7032109.43066144, 7034251.5625, 7042021.489379221, 7043996.875, 7052110.9375, 7053718.64135675, 7054299.6798326075, 7057560.276968299, 7059480.162380217, 7060811.77259277, 7062554.8266791515, 7063186.948298047, 7064594.960696833, 7065919.940128694, 7069554.6875, 7069676.118294685, 7070306.413759812, 7071677.179196502, 7072155.999276852, 7075406.902509678, 7075411.784833743, 7076497.912771927, 7077178.748001068, 7077952.966489455, 7079237.673905073, 7079435.519231649, 7082707.399083979, 7086461.382597969, 7087673.872517844, 7088955.940033788, 7089304.323921151, 7090448.405209036, 7090507.1838376885, 7092968.75, 7093257.8125, 7094675.0, 7099199.818232614, 7106135.854509128, 7111083.314558238, 7115721.303995722, 7118423.677889852, 7123196.875, 7123208.406824104, 7127319.811735158, 7131226.096245381, 7131531.772857061, 7136480.382267867, 7137679.6875, 7141068.807337026, 7143090.625, 7153533.3499095, 7153923.904151928, 7166920.3125, 7168064.0625, 7174831.045751801, 7177429.08311779, 7178807.046372351, 7179002.052617889, 7179455.518584996, 7180581.900047146, 7181371.875, 7181386.889180274, 7184402.406058726, 7184737.5, 7186859.375, 7194644.2726303805, 7198633.841741391, 7204351.553327817, 7209897.522612095, 7214095.3125, 7215396.318173606, 7217396.377757388, 7220332.055667191, 7231883.6093235295, 7244774.711913046, 7248542.446329602, 7260450.163488039, 7272235.935841728, 7273263.680975504, 7277593.75, 7283084.375, 7296131.25, 7297169.541810198, 7311802.158685105, 7316646.78175492, 7318259.375, 7320298.637116808, 7322367.630368029, 7322416.177918096, 7322933.77513504, 7324717.896715794, 7330125.8831184255, 7334932.654152049, 7341517.1875, 7341560.9375, 7343160.9375, 7360978.702128716, 7366147.549165035, 7370115.919074487, 7374157.101263578, 7383192.030639198, 7383194.673123452, 7385265.625, 7396125.443678394, 7403463.434886904, 7410735.4358977135, 7417703.125, 7417775.572758404, 7419913.934088181, 7425251.870221939, 7429350.0, 7433251.842459021, 7433646.076638063, 7433744.028405825, 7445418.598797604, 7449007.8125, 7451329.410092838, 7455553.595693597, 7457114.222968031, 7461385.9375, 7475319.68637211, 7481364.0625, 7481763.982619709, 7481813.865356766, 7482945.922651834, 7484432.059944403, 7485832.222239301, 7487999.924963288, 7491304.563842744, 7492060.786520454, 7499286.03838995, 7518776.15149659, 7518847.394523694, 7531252.610251022, 7534707.595795229, 7539795.765466375, 7551064.858113601, 7568617.1875, 7575163.949322625, 7582510.202324729, 7585976.190737988, 7588896.875, 7591698.788442281, 7643392.13328214, 7646962.113922395, 7649611.20102011, 7671263.345297058, 7691327.42148156, 7710009.7152959, 7718834.78317065, 7746721.604688283, 7757493.601930725, 7798804.389315552, 7826642.1875, 8013757.956788702, 8077939.0625, 8130068.75, 8151998.4375, 8152100.821255108, 8305703.125, 8415407.628299989, 8896535.23400131, 9789375.117735244, 9852244.69609916, 10405777.701116268, 10467177.466736976, 10520561.20661951, 10564724.606593326, 10579152.649226941, 10605637.947142633, 10654496.875, 10660319.040646847, 10699251.216312686, 10769570.542770673, 10815071.875, 10859384.487486435, 10950733.825077353, 10966074.20728845, 11023655.363445703, 11023715.16195357, 11040703.955105523, 11125486.327487502, 11125498.92967076, 11144811.185865572, 11147740.310647408, 11162409.787489574, 11210626.687112886, 11215238.07866137, 11263245.3125, 11284743.452293929, 11313321.858444536, 11325637.5, 11359357.45875484, 11375595.3125, 11458699.151188212, 11466123.110315347, 11509106.319945596, 11521186.707100084, 11532137.817812612, 11548196.593452753, 11551781.25, 11588440.49719893, 11611078.125, 11613115.625, 11664132.315269183, 11677733.127305765, 11678416.712696565, 11704965.625, 11734585.9375, 11742937.5, 11760975.095458027, 11791107.182924503, 11819824.554404428, 11821983.946317345, 11834204.6875, 11876764.65200226, 11878132.8125, 11896499.105642298, 11908143.10512788, 11911215.625, 11934916.84448819, 11938329.6875, 11980809.83617919, 12002333.511217179, 12006390.098545287, 12075260.9375, 12085235.9375, 12098443.331469974, 12123302.852603862, 12127043.173896488, 12178870.005816469, 12181725.60153862, 12226050.670633523, 12236124.731765738, 12248092.163772088, 12292051.78613527, 12303845.887398228, 12324060.170074416, 12344226.518092453, 12353876.5625, 12371062.958381189, 12425622.78022515, 12497515.135917256, 12501923.430186063, 12547076.634274248, 12549834.273760656, 12557903.607911628, 12587467.117977092, 12611262.672294455, 12669396.564752629, 12671462.5, 12671571.542608751, 12671679.6875, 12680839.0625, 12696104.6875, 12751219.675828837, 12766687.202837577, 12776142.237474922, 12778873.4375, 12826359.960741932, 12864698.4375, 12904932.8125, 12906733.581488125, 12910429.6875, 12912159.787484478, 12923238.818687357, 12995738.796300104, 13053302.653925, 13066264.807913376, 13105453.880543645, 13134266.018862609, 13143799.737385549, 13183479.222262494, 13187346.401760308, 13192376.5625, 13208086.42013418, 13212172.933846084, 13220994.886442848, 13222825.784947636, 13237331.31228769, 13237516.411811348, 13256157.520928895, 13280460.252598384, 13297696.531420432, 13303355.648557529, 13326878.44834986, 13332312.5, 13336596.636951024, 13342112.14194723, 13374319.610202905, 13374331.545349503, 13385670.683836093, 13393732.777929438, 13417545.586610572, 13417989.48557362, 13419925.36773336, 13431107.156304933, 13457082.835724477, 13457948.455262294, 13460867.06422994, 13468678.920015655, 13469367.1875, 13479339.64240139, 13496035.9375, 13498180.081452075, 13498304.393954448, 13503491.208253877, 13517572.431725089, 13525687.417931704, 13551563.193934044, 13554934.612031162, 13555589.0625, 13556738.831478262, 13566140.834861567, 13567221.875, 13567562.5, 13589725.133138914, 13610787.5, 13620124.955963738, 13622606.991826974, 13623079.667350292, 13623342.761467515, 13625411.163905488, 13626034.515583176, 13626268.198861016, 13629561.960817927, 13673682.871533323, 13684306.250068264, 13687718.648866095, 13697279.505606005, 13705484.375, 13709692.576316193, 13714731.25, 13714959.905672263, 13760138.161014616, 13763331.122411804, 13790038.820078395, 13844170.786839887, 13844212.099482466, 13851028.059902573, 13861474.230245875, 13866315.41452956, 13874507.096958136, 13878731.097244963, 13885627.857059715, 13891667.919879634, 13897539.0625, 13911759.65099801, 13914119.696912253, 13914785.603167066, 13931431.76309551, 13931579.257231202, 13938797.45668302, 13947502.799169827, 13949985.9375, 13958570.3125, 13986615.876107745, 14014105.935385259, 14028500.0, 14039227.882173218, 14047759.375, 14067078.284353012, 14067671.924098322, 14070810.9375, 14099833.348685153, 14105589.0625, 14131776.885052864, 14137525.0, 14151551.020488773, 14167248.4375, 14168974.754252618, 14342466.294668982, 14373752.398337016, 14379141.647287978, 14427713.974990902, 14457939.87276123, 14470221.875, 14477551.588301234, 14581937.5, 14741059.375, 14792668.502805809, 14798075.0, 14798323.438254615, 14824510.566071574, 14833496.875, 14841736.615851704, 14845901.421011275, 14875315.63068344, 14887735.735675896, 14945123.165863588, 15004876.166372702, 15004889.0625, 15709889.0625, 17069522.929147452, 17299114.43715722, 17314383.625766963, 17361295.910388537, 17367556.019878175, 17387342.1875, 17395277.70673881, 17398797.85850792, 17407473.228699677, 17414694.656451683, 17416327.638028655, 17420373.801409166, 17422911.452424914, 17423957.26303541, 17424162.5, 17427920.263256446, 17442637.40244545, 17446295.018775973, 17448063.02536805, 17468964.468700606, 17484452.14979487, 17486437.43446301, 17490357.8125, 17531687.429321755, 17538085.58961876, 17538571.68998395, 17539127.473311298, 17551430.671352863, 17555147.580319773, 17555221.125538804, 17560066.481300715, 17638278.13989996, 336124202.52591705], [52.184091793816975, 8.000041187308591, 9.385920209364965, 77.21459060921242, 13.451809496003616, 43.98308450743192, 60.23730819009052, 25.21519704467307, 64.7237789546834, 7.926865929694985, 78.60484411621515, 94.81025112142927, 26.020857607538858, 75.07670718512074, 6.761047328106072, 6.764175709236198, 71.8087308575958, 26.053869816392087, 23.891888445962778, 22.751698467997162, 17.912791088750918, 72.83497014239761, 9.135191989141818, 10.765580431748635, 14.853145766194656, 17.469459759503238, 49.13899516612646, 15.892850198626467, 38.7274022119345, 30.560074181090208, 36.88767128470816, 68.76682116345336, 38.69266575111246, 34.81017535350598, 37.87378118530908, 25.512392810314353, 8.652170698578146, 58.655004228252636, 44.47235895041948, 23.918681892455595, 8.10251314853253, 55.3568096164301, 5.077857726426501, 23.61265360079565, 41.871087723345184, 35.906487138577276, 42.573057911584826, 69.7804047375858, 12.513487384248945, 35.35152305286264, 68.3518035472774, 26.50308289151804, 7.657899339515086, 41.94752465297758, 43.16289641100993, 5.383389028934666, 7.971208992949851, 74.35215335918157, 6.817798672226604, 19.537419597565545, 47.26980633224066, 61.5424491779041, 61.8605849044191, 8.191204680954392, 8.261552645615602, 43.93884734552492, 8.266491230969137, 5.253234450333963, 9.722202254776686, 31.9580250930399, 7.40833052542899, 9.151925285738606, 41.27927267039169, 20.651467407546853, 98.69326749268916, 53.940736265008965, 58.323210022053516, 48.694891030978816, 61.66753433040779, 20.568140004217593, 39.001284184421394, 20.02842232641656, 17.124030170985115, 33.03286184746271, 8.361574945300097, 54.03135260434296, 6.223750749466156, 32.3130619137347, 61.31263869750808, 22.068343457244712, 15.604719001519141, 8.642043201307008, 49.213885642234985, 7.857615874556166, 74.18758433981594, 44.619238477536314, 8.646007830935265, 10.888665373703887, 61.74781857076299, 97.99409051793238, 49.01627555413962, 11.902828525069237, 5.73501456777269, 7.494847833514609, 7.029048778422687, 11.077606293895254, 52.743806836157226, 81.13299197501658, 32.80565299543323, 31.761954000640998, 10.855001284849589, 5.609469325706591, 12.014982295037415, 18.236332756694406, 57.4427381877077, 43.104259973214944, 5.992889706771412, 20.12548812585788, 30.19946596342023, 12.056767990520193, 25.256964455685118, 34.44416496309503, 60.25892926598113, 50.356589879113265, 7.809905893993661, 6.860388754497189, 38.54637552048089, 91.88919376986767, 5.52324916898747, 13.048978691046507, 8.757967799445296, 17.209035291686625, 7.358514322706778, 90.08050797698263, 76.41467856161766, 41.83278485501796, 71.43122507136941, 53.13854357297615, 42.92728499135152, 94.97302282927551, 64.8076135989443, 58.50915208481255, 24.774067079441593, 54.48480224663971, 10.35622240368168, 41.25627812221201, 36.22668018989738, 12.65290490414397, 49.34417320040431, 16.02923147841358, 15.60354964030981, 18.06206820529519, 13.66267288753798, 71.70025475295485, 36.747746656235, 44.438369402779124, 7.917139669181635, 12.207786136474148, 14.597627270569857, 37.43312247211789, 15.761159513684962, 9.736221770882924, 8.602342486289114, 13.060615042411076, 45.02799979352094, 48.30361665740414, 9.201006345892084, 58.14908977436397, 58.59782787130499, 46.526617203023136, 8.617326574934927, 15.659634155054272, 11.620733668318712, 24.14285253719602, 7.547472258387531, 41.891624092312604, 34.85510991184974, 38.15685175688678, 18.344434210221618, 58.818678487691294, 53.92573661049925, 12.158900893640647, 14.473655982495778, 7.204762095808517, 47.94566211746255, 6.820715618970731, 10.073520186051885, 5.0462670083552394, 56.566095771800846, 49.03081850929036, 24.04368849379824, 60.764844750655875, 15.590571730825879, 38.25927975787228, 10.272298962603921, 38.24395100063418, 14.64554024552145, 6.456625077138597, 47.508720384112486, 19.29631078394827, 13.085459023507754, 21.1919615338223, 10.939771290777413, 23.127726029199525, 5.537997166782251, 7.017446017268324, 5.307409338260002, 8.931040135257902, 32.16699490617025, 9.746219083037845, 55.3375939693899, 46.2642658228088, 35.207911056264216, 29.543076148661687, 10.59241976783132, 59.79628938965598, 11.72445214287882, 47.35757002528918, 23.021690707204293, 11.56620064749119, 12.768099004114301, 6.948033338091945, 61.19772817297196, 22.172583575672864, 14.002439573958537, 8.209797758167843, 13.914981745130515, 49.72575065726217, 19.445830219808148, 60.428355480353936, 56.86217398186792, 67.60956231710293, 5.0292699982072415, 59.006292352770764, 7.537947171963403, 20.325167496226612, 15.410088770745258, 5.887464941993215, 12.338858046378842, 20.047377305069343, 81.86918185981216, 132.54048527442228, 28.040579601643188, 16.972140463292074, 73.56137613476818, 6.118854592998611, 23.178162690019924, 11.507777730828455, 103.05126997654735, 61.186759651763346, 134.48507149223772, 18.997824839564935, 18.344649820215526, 5.033826140470363, 5.761706479540241, 16.670267532491412, 28.013217893018993, 9.55062717581921, 9.079029987607063, 89.32639157372306, 34.79003542728762, 6.343185267566669, 57.1566900297415, 54.9977730858003, 61.317479699017554, 106.84170565411802, 16.505296361276468, 64.0709718739057, 18.600225466707236, 29.13347960370785, 23.928186152365367, 48.00395318223027, 8.462110141758568, 10.887202007300303, 83.48217031311705, 29.042827752447042, 59.93287465100778, 36.845837363068775, 66.86863134167604, 5.9672977243176355, 29.268332280668442, 59.69843478290925, 47.58290269258862, 43.645542405509644, 19.980506041524862, 12.98464637556844, 22.053541699749136, 22.235305407442336, 59.0851599454004, 8.625025349713995, 56.05223598739221, 12.361836912108249, 55.10758371821497, 89.26059860528872, 13.777596267339112, 87.03942224782665, 17.048550292954346, 33.651029682766506, 33.189481897577565, 101.33393949041167, 94.48584221926666, 52.93001693308808, 9.11122933780715, 40.84896474373514, 8.130793802350777, 6.889850494537435, 31.227424900259454, 41.94995628880258, 9.311420150305832, 19.85847316123365, 28.181087768076807, 19.184373300162132, 140.97667171733235, 43.2898869379084, 25.449554899675178, 12.930072575517297, 55.8284285755228, 6.485945317900089, 57.50123456354396, 26.69638656594714, 15.130446122025557, 89.25411227497271, 17.004226417821275, 7.494752782109766, 75.03805109188852, 7.487284513900317, 56.332493027363924, 146.24187644762105, 44.423127282205684, 38.48555775136522, 48.318903598690234, 27.38251307416874, 126.97836279159674, 5.8586297845842354, 29.193796181336754, 64.77422062014851, 33.27478271396353, 65.30178166897642, 11.840075859156377, 34.0615166126331, 5.520299008501564, 157.07416641740588, 18.822380728519843, 12.557259039282332, 7.228983636397862, 24.24033307126636, 11.79015976871047, 76.70270523874484, 64.66575186249148, 12.491081370478739, 94.81221160865641, 22.26385672636749, 10.096826970464756, 107.53764764441608, 81.24905468655976, 10.700622679252872, 19.827657186761947, 39.89251460109547, 46.88683872330879, 16.336029584743645, 38.477600565063284, 120.05844374135889, 68.75158764823911, 19.68746293387615, 26.677591950366313, 7.743301881493542, 8.268380052414974, 40.96479151929429, 24.983465745515268, 7.741885293083955, 27.66861467558177, 25.272940123516097, 15.0096459394469, 21.671177740293412, 99.40908313834669, 74.79435617150136, 5.48721949363613, 101.33053217283268, 16.690740662212026, 49.36982611587534, 12.180726844509678, 30.00738740003372, 20.422881777833773, 93.14777401499248, 6.524901161546767, 171.74241681634726, 14.128469658645328, 5.543353218247235, 47.167758694575284, 46.30041241425348, 14.61841245994605, 45.749860262099276, 39.48385026730473, 61.61657488015735, 38.24385519010955, 7.280972126193216, 31.549198138622124, 8.362346024784586, 88.40424285992091, 22.426877267652724, 16.77422040981621, 14.941769621837086, 82.57119249341278, 7.538748413944334, 57.09653881993552, 7.750191786431132, 62.02152536110484, 17.440878204648627, 73.08682067420676, 14.476525187934236, 32.892847334636016, 18.31049395518196, 27.87418301483283, 20.375755861403228, 10.105609876396027, 40.884518279347944, 26.604967154373206, 21.62200747307059, 7.509120195767349, 67.32107298197641, 70.47724724423693, 24.81604035344633, 7.394796111771666, 53.328115998419875, 47.546614720263264, 47.7948738566419, 6.582312226585585, 6.226891447147788, 57.80126089261969, 42.792565534391905, 20.86151152904177, 5.3865812856304025, 88.75432634994053, 8.162209353421876, 71.14950041393638, 15.58578770292082, 8.312713125528532, 44.80347889735715, 41.46920786919761, 64.12426501695356, 23.43185598629856, 80.90459048113037, 6.659916575271539, 73.05965010832362, 63.32880390176372, 69.77479226795391, 150.14806259193728, 16.846913084055895, 42.63555828677881, 19.29871414564431, 29.871432313734346, 11.168276705846576, 87.25722845872136, 53.834056621319235, 46.165435852153536, 37.24586493886418, 14.944386441314396, 10.592862260645717, 33.65395675453615, 67.89050482372113, 14.428317350873861, 5.533356440522385, 27.024031747349987, 24.119937030166, 13.429312804892243, 60.26240745232952, 51.364845759858255, 54.6766621212557, 5.842231723744949, 18.326326624570772, 17.844100445670385, 50.79873514861414, 46.90751059176918, 22.349174142263625, 13.114789321874598, 104.13415307285484, 12.91806246828783, 75.03929809004569, 7.740382916554944, 5.740332090304366, 10.365277117312147, 38.00611466447616, 111.0554599568406, 46.668514814600194, 23.831240950249413, 59.43215443542239, 11.828104269398732, 31.262038521614144, 34.16269539010579, 49.1887419035139, 11.019000488036555, 74.24313639111887, 53.54759397848636, 20.102565845710767, 29.218945824328877, 57.5651938389466, 5.370575437886704, 12.79469428453429, 7.125373919670654, 20.55853576324594])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([5176359.375, 5278197.853278245, 5280165.688778647, 5400470.240403823, 5401085.482786234, 5415904.6875, 5421084.375, 5421088.869105486, 5426535.709674103, 5457611.751633061, 5528860.9375, 5553021.875, 5554602.984113993, 5576199.208337134, 5576243.73766927, 5577168.753350776, 5584078.125, 5589565.36641185, 5605721.994601503, 5608462.117545506, 5624286.546441099, 5642694.575049502, 5643606.008462296, 5698043.505273514, 5717114.638477524, 5733592.813033938, 5755443.75, 5796592.416548447, 6674254.6875, 6698932.8125, 6778443.830400588, 6788427.690742854, 6797823.4375, 6798121.875, 6801829.319750751, 6806044.2668590825, 6806149.070278282, 6808251.5625, 6809051.5625, 6831962.848255071, 6835660.572420684, 6865896.13556999, 6869433.299348795, 6871527.8141121445, 6881832.8125, 6897929.761116174, 6909714.138681498, 6909720.109257517, 6910810.111330104, 6914277.03586287, 6920757.8125, 6924763.951216762, 6929736.086144281, 6931215.625, 6933735.9375, 6934251.162219503, 6935482.16439149, 6944615.625, 6956656.5159969255, 6975016.22892419, 6979668.042256201, 6981323.2284298055, 7000471.694890829, 7002713.87768687, 7006943.702376271, 7011939.150588079, 7015379.649322781, 7019842.473972929, 7019972.30364368, 7022898.528069976, 7025258.801919314, 7025926.342328861, 7028357.512768078, 7029591.957604919, 7032109.43066144, 7034251.5625, 7042021.489379221, 7043996.875, 7052110.9375, 7053718.64135675, 7054299.6798326075, 7057560.276968299, 7059480.162380217, 7060811.77259277, 7062554.8266791515, 7063186.948298047, 7064594.960696833, 7065919.940128694, 7069554.6875, 7069676.118294685, 7070306.413759812, 7071677.179196502, 7072155.999276852, 7075406.902509678, 7075411.784833743, 7076497.912771927, 7077178.748001068, 7077952.966489455, 7079237.673905073, 7079435.519231649, 7082707.399083979, 7086461.382597969, 7087673.872517844, 7088955.940033788, 7089304.323921151, 7090448.405209036, 7090507.1838376885, 7092968.75, 7093257.8125, 7094675.0, 7099199.818232614, 7106135.854509128, 7111083.314558238, 7115721.303995722, 7118423.677889852, 7123196.875, 7123208.406824104, 7127319.811735158, 7131226.096245381, 7131531.772857061, 7136480.382267867, 7137679.6875, 7141068.807337026, 7143090.625, 7153533.3499095, 7153923.904151928, 7166920.3125, 7168064.0625, 7174831.045751801, 7177429.08311779, 7178807.046372351, 7179002.052617889, 7179455.518584996, 7180581.900047146, 7181371.875, 7181386.889180274, 7184402.406058726, 7184737.5, 7186859.375, 7194644.2726303805, 7198633.841741391, 7204351.553327817, 7209897.522612095, 7214095.3125, 7215396.318173606, 7217396.377757388, 7220332.055667191, 7231883.6093235295, 7244774.711913046, 7248542.446329602, 7260450.163488039, 7272235.935841728, 7273263.680975504, 7277593.75, 7283084.375, 7296131.25, 7297169.541810198, 7311802.158685105, 7316646.78175492, 7318259.375, 7320298.637116808, 7322367.630368029, 7322416.177918096, 7322933.77513504, 7324717.896715794, 7330125.8831184255, 7334932.654152049, 7341517.1875, 7341560.9375, 7343160.9375, 7360978.702128716, 7366147.549165035, 7370115.919074487, 7374157.101263578, 7383192.030639198, 7383194.673123452, 7385265.625, 7396125.443678394, 7403463.434886904, 7410735.4358977135, 7417703.125, 7417775.572758404, 7419913.934088181, 7425251.870221939, 7429350.0, 7433251.842459021, 7433646.076638063, 7433744.028405825, 7445418.598797604, 7449007.8125, 7451329.410092838, 7455553.595693597, 7457114.222968031, 7461385.9375, 7475319.68637211, 7481364.0625, 7481763.982619709, 7481813.865356766, 7482945.922651834, 7484432.059944403, 7485832.222239301, 7487999.924963288, 7491304.563842744, 7492060.786520454, 7499286.03838995, 7518776.15149659, 7518847.394523694, 7531252.610251022, 7534707.595795229, 7539795.765466375, 7551064.858113601, 7568617.1875, 7575163.949322625, 7582510.202324729, 7585976.190737988, 7588896.875, 7591698.788442281, 7643392.13328214, 7646962.113922395, 7649611.20102011, 7671263.345297058, 7691327.42148156, 7710009.7152959, 7718834.78317065, 7746721.604688283, 7757493.601930725, 7798804.389315552, 7826642.1875, 8013757.956788702, 8077939.0625, 8130068.75, 8151998.4375, 8152100.821255108, 8305703.125, 8415407.628299989, 8896535.23400131, 9789375.117735244, 9852244.69609916, 10405777.701116268, 10467177.466736976, 10520561.20661951, 10564724.606593326, 10579152.649226941, 10605637.947142633, 10654496.875, 10660319.040646847, 10699251.216312686, 10769570.542770673, 10815071.875, 10859384.487486435, 10950733.825077353, 10966074.20728845, 11023655.363445703, 11023715.16195357, 11040703.955105523, 11125486.327487502, 11125498.92967076, 11144811.185865572, 11147740.310647408, 11162409.787489574, 11210626.687112886, 11215238.07866137, 11263245.3125, 11284743.452293929, 11313321.858444536, 11325637.5, 11359357.45875484, 11375595.3125, 11458699.151188212, 11466123.110315347, 11509106.319945596, 11521186.707100084, 11532137.817812612, 11548196.593452753, 11551781.25, 11588440.49719893, 11611078.125, 11613115.625, 11664132.315269183, 11677733.127305765, 11678416.712696565, 11704965.625, 11734585.9375, 11742937.5, 11760975.095458027, 11791107.182924503, 11819824.554404428, 11821983.946317345, 11834204.6875, 11876764.65200226, 11878132.8125, 11896499.105642298, 11908143.10512788, 11911215.625, 11934916.84448819, 11938329.6875, 11980809.83617919, 12002333.511217179, 12006390.098545287, 12075260.9375, 12085235.9375, 12098443.331469974, 12123302.852603862, 12127043.173896488, 12178870.005816469, 12181725.60153862, 12226050.670633523, 12236124.731765738, 12248092.163772088, 12292051.78613527, 12303845.887398228, 12324060.170074416, 12344226.518092453, 12353876.5625, 12371062.958381189, 12425622.78022515, 12497515.135917256, 12501923.430186063, 12547076.634274248, 12549834.273760656, 12557903.607911628, 12587467.117977092, 12611262.672294455, 12669396.564752629, 12671462.5, 12671571.542608751, 12671679.6875, 12680839.0625, 12696104.6875, 12751219.675828837, 12766687.202837577, 12776142.237474922, 12778873.4375, 12826359.960741932, 12864698.4375, 12904932.8125, 12906733.581488125, 12910429.6875, 12912159.787484478, 12923238.818687357, 12995738.796300104, 13053302.653925, 13066264.807913376, 13105453.880543645, 13134266.018862609, 13143799.737385549, 13183479.222262494, 13187346.401760308, 13192376.5625, 13208086.42013418, 13212172.933846084, 13220994.886442848, 13222825.784947636, 13237331.31228769, 13237516.411811348, 13256157.520928895, 13280460.252598384, 13297696.531420432, 13303355.648557529, 13326878.44834986, 13332312.5, 13336596.636951024, 13342112.14194723, 13374319.610202905, 13374331.545349503, 13385670.683836093, 13393732.777929438, 13417545.586610572, 13417989.48557362, 13419925.36773336, 13431107.156304933, 13457082.835724477, 13457948.455262294, 13460867.06422994, 13468678.920015655, 13469367.1875, 13479339.64240139, 13496035.9375, 13498180.081452075, 13498304.393954448, 13503491.208253877, 13517572.431725089, 13525687.417931704, 13551563.193934044, 13554934.612031162, 13555589.0625, 13556738.831478262, 13566140.834861567, 13567221.875, 13567562.5, 13589725.133138914, 13610787.5, 13620124.955963738, 13622606.991826974, 13623079.667350292, 13623342.761467515, 13625411.163905488, 13626034.515583176, 13626268.198861016, 13629561.960817927, 13673682.871533323, 13684306.250068264, 13687718.648866095, 13697279.505606005, 13705484.375, 13709692.576316193, 13714731.25, 13714959.905672263, 13760138.161014616, 13763331.122411804, 13790038.820078395, 13844170.786839887, 13844212.099482466, 13851028.059902573, 13861474.230245875, 13866315.41452956, 13874507.096958136, 13878731.097244963, 13885627.857059715, 13891667.919879634, 13897539.0625, 13911759.65099801, 13914119.696912253, 13914785.603167066, 13931431.76309551, 13931579.257231202, 13938797.45668302, 13947502.799169827, 13949985.9375, 13958570.3125, 13986615.876107745, 14014105.935385259, 14028500.0, 14039227.882173218, 14047759.375, 14067078.284353012, 14067671.924098322, 14070810.9375, 14099833.348685153, 14105589.0625, 14131776.885052864, 14137525.0, 14151551.020488773, 14167248.4375, 14168974.754252618, 14342466.294668982, 14373752.398337016, 14379141.647287978, 14427713.974990902, 14457939.87276123, 14470221.875, 14477551.588301234, 14581937.5, 14741059.375, 14792668.502805809, 14798075.0, 14798323.438254615, 14824510.566071574, 14833496.875, 14841736.615851704, 14845901.421011275, 14875315.63068344, 14887735.735675896, 14945123.165863588, 15004876.166372702, 15004889.0625, 15709889.0625, 17069522.929147452, 17299114.43715722, 17314383.625766963, 17361295.910388537, 17367556.019878175, 17387342.1875, 17395277.70673881, 17398797.85850792, 17407473.228699677, 17414694.656451683, 17416327.638028655, 17420373.801409166, 17422911.452424914, 17423957.26303541, 17424162.5, 17427920.263256446, 17442637.40244545, 17446295.018775973, 17448063.02536805, 17468964.468700606, 17484452.14979487, 17486437.43446301, 17490357.8125, 17531687.429321755, 17538085.58961876, 17538571.68998395, 17539127.473311298, 17551430.671352863, 17555147.580319773, 17555221.125538804, 17560066.481300715, 17638278.13989996, 336124202.52591705], [52.184091793816975, 8.000041187308591, 9.385920209364965, 77.21459060921242, 13.451809496003616, 43.98308450743192, 60.23730819009052, 25.21519704467307, 64.7237789546834, 7.926865929694985, 78.60484411621515, 94.81025112142927, 26.020857607538858, 75.07670718512074, 6.761047328106072, 6.764175709236198, 71.8087308575958, 26.053869816392087, 23.891888445962778, 22.751698467997162, 17.912791088750918, 72.83497014239761, 9.135191989141818, 10.765580431748635, 14.853145766194656, 17.469459759503238, 49.13899516612646, 15.892850198626467, 38.7274022119345, 30.560074181090208, 36.88767128470816, 68.76682116345336, 38.69266575111246, 34.81017535350598, 37.87378118530908, 25.512392810314353, 8.652170698578146, 58.655004228252636, 44.47235895041948, 23.918681892455595, 8.10251314853253, 55.3568096164301, 5.077857726426501, 23.61265360079565, 41.871087723345184, 35.906487138577276, 42.573057911584826, 69.7804047375858, 12.513487384248945, 35.35152305286264, 68.3518035472774, 26.50308289151804, 7.657899339515086, 41.94752465297758, 43.16289641100993, 5.383389028934666, 7.971208992949851, 74.35215335918157, 6.817798672226604, 19.537419597565545, 47.26980633224066, 61.5424491779041, 61.8605849044191, 8.191204680954392, 8.261552645615602, 43.93884734552492, 8.266491230969137, 5.253234450333963, 9.722202254776686, 31.9580250930399, 7.40833052542899, 9.151925285738606, 41.27927267039169, 20.651467407546853, 98.69326749268916, 53.940736265008965, 58.323210022053516, 48.694891030978816, 61.66753433040779, 20.568140004217593, 39.001284184421394, 20.02842232641656, 17.124030170985115, 33.03286184746271, 8.361574945300097, 54.03135260434296, 6.223750749466156, 32.3130619137347, 61.31263869750808, 22.068343457244712, 15.604719001519141, 8.642043201307008, 49.213885642234985, 7.857615874556166, 74.18758433981594, 44.619238477536314, 8.646007830935265, 10.888665373703887, 61.74781857076299, 97.99409051793238, 49.01627555413962, 11.902828525069237, 5.73501456777269, 7.494847833514609, 7.029048778422687, 11.077606293895254, 52.743806836157226, 81.13299197501658, 32.80565299543323, 31.761954000640998, 10.855001284849589, 5.609469325706591, 12.014982295037415, 18.236332756694406, 57.4427381877077, 43.104259973214944, 5.992889706771412, 20.12548812585788, 30.19946596342023, 12.056767990520193, 25.256964455685118, 34.44416496309503, 60.25892926598113, 50.356589879113265, 7.809905893993661, 6.860388754497189, 38.54637552048089, 91.88919376986767, 5.52324916898747, 13.048978691046507, 8.757967799445296, 17.209035291686625, 7.358514322706778, 90.08050797698263, 76.41467856161766, 41.83278485501796, 71.43122507136941, 53.13854357297615, 42.92728499135152, 94.97302282927551, 64.8076135989443, 58.50915208481255, 24.774067079441593, 54.48480224663971, 10.35622240368168, 41.25627812221201, 36.22668018989738, 12.65290490414397, 49.34417320040431, 16.02923147841358, 15.60354964030981, 18.06206820529519, 13.66267288753798, 71.70025475295485, 36.747746656235, 44.438369402779124, 7.917139669181635, 12.207786136474148, 14.597627270569857, 37.43312247211789, 15.761159513684962, 9.736221770882924, 8.602342486289114, 13.060615042411076, 45.02799979352094, 48.30361665740414, 9.201006345892084, 58.14908977436397, 58.59782787130499, 46.526617203023136, 8.617326574934927, 15.659634155054272, 11.620733668318712, 24.14285253719602, 7.547472258387531, 41.891624092312604, 34.85510991184974, 38.15685175688678, 18.344434210221618, 58.818678487691294, 53.92573661049925, 12.158900893640647, 14.473655982495778, 7.204762095808517, 47.94566211746255, 6.820715618970731, 10.073520186051885, 5.0462670083552394, 56.566095771800846, 49.03081850929036, 24.04368849379824, 60.764844750655875, 15.590571730825879, 38.25927975787228, 10.272298962603921, 38.24395100063418, 14.64554024552145, 6.456625077138597, 47.508720384112486, 19.29631078394827, 13.085459023507754, 21.1919615338223, 10.939771290777413, 23.127726029199525, 5.537997166782251, 7.017446017268324, 5.307409338260002, 8.931040135257902, 32.16699490617025, 9.746219083037845, 55.3375939693899, 46.2642658228088, 35.207911056264216, 29.543076148661687, 10.59241976783132, 59.79628938965598, 11.72445214287882, 47.35757002528918, 23.021690707204293, 11.56620064749119, 12.768099004114301, 6.948033338091945, 61.19772817297196, 22.172583575672864, 14.002439573958537, 8.209797758167843, 13.914981745130515, 49.72575065726217, 19.445830219808148, 60.428355480353936, 56.86217398186792, 67.60956231710293, 5.0292699982072415, 59.006292352770764, 7.537947171963403, 20.325167496226612, 15.410088770745258, 5.887464941993215, 12.338858046378842, 20.047377305069343, 81.86918185981216, 132.54048527442228, 28.040579601643188, 16.972140463292074, 73.56137613476818, 6.118854592998611, 23.178162690019924, 11.507777730828455, 103.05126997654735, 61.186759651763346, 134.48507149223772, 18.997824839564935, 18.344649820215526, 5.033826140470363, 5.761706479540241, 16.670267532491412, 28.013217893018993, 9.55062717581921, 9.079029987607063, 89.32639157372306, 34.79003542728762, 6.343185267566669, 57.1566900297415, 54.9977730858003, 61.317479699017554, 106.84170565411802, 16.505296361276468, 64.0709718739057, 18.600225466707236, 29.13347960370785, 23.928186152365367, 48.00395318223027, 8.462110141758568, 10.887202007300303, 83.48217031311705, 29.042827752447042, 59.93287465100778, 36.845837363068775, 66.86863134167604, 5.9672977243176355, 29.268332280668442, 59.69843478290925, 47.58290269258862, 43.645542405509644, 19.980506041524862, 12.98464637556844, 22.053541699749136, 22.235305407442336, 59.0851599454004, 8.625025349713995, 56.05223598739221, 12.361836912108249, 55.10758371821497, 89.26059860528872, 13.777596267339112, 87.03942224782665, 17.048550292954346, 33.651029682766506, 33.189481897577565, 101.33393949041167, 94.48584221926666, 52.93001693308808, 9.11122933780715, 40.84896474373514, 8.130793802350777, 6.889850494537435, 31.227424900259454, 41.94995628880258, 9.311420150305832, 19.85847316123365, 28.181087768076807, 19.184373300162132, 140.97667171733235, 43.2898869379084, 25.449554899675178, 12.930072575517297, 55.8284285755228, 6.485945317900089, 57.50123456354396, 26.69638656594714, 15.130446122025557, 89.25411227497271, 17.004226417821275, 7.494752782109766, 75.03805109188852, 7.487284513900317, 56.332493027363924, 146.24187644762105, 44.423127282205684, 38.48555775136522, 48.318903598690234, 27.38251307416874, 126.97836279159674, 5.8586297845842354, 29.193796181336754, 64.77422062014851, 33.27478271396353, 65.30178166897642, 11.840075859156377, 34.0615166126331, 5.520299008501564, 157.07416641740588, 18.822380728519843, 12.557259039282332, 7.228983636397862, 24.24033307126636, 11.79015976871047, 76.70270523874484, 64.66575186249148, 12.491081370478739, 94.81221160865641, 22.26385672636749, 10.096826970464756, 107.53764764441608, 81.24905468655976, 10.700622679252872, 19.827657186761947, 39.89251460109547, 46.88683872330879, 16.336029584743645, 38.477600565063284, 120.05844374135889, 68.75158764823911, 19.68746293387615, 26.677591950366313, 7.743301881493542, 8.268380052414974, 40.96479151929429, 24.983465745515268, 7.741885293083955, 27.66861467558177, 25.272940123516097, 15.0096459394469, 21.671177740293412, 99.40908313834669, 74.79435617150136, 5.48721949363613, 101.33053217283268, 16.690740662212026, 49.36982611587534, 12.180726844509678, 30.00738740003372, 20.422881777833773, 93.14777401499248, 6.524901161546767, 171.74241681634726, 14.128469658645328, 5.543353218247235, 47.167758694575284, 46.30041241425348, 14.61841245994605, 45.749860262099276, 39.48385026730473, 61.61657488015735, 38.24385519010955, 7.280972126193216, 31.549198138622124, 8.362346024784586, 88.40424285992091, 22.426877267652724, 16.77422040981621, 14.941769621837086, 82.57119249341278, 7.538748413944334, 57.09653881993552, 7.750191786431132, 62.02152536110484, 17.440878204648627, 73.08682067420676, 14.476525187934236, 32.892847334636016, 18.31049395518196, 27.87418301483283, 20.375755861403228, 10.105609876396027, 40.884518279347944, 26.604967154373206, 21.62200747307059, 7.509120195767349, 67.32107298197641, 70.47724724423693, 24.81604035344633, 7.394796111771666, 53.328115998419875, 47.546614720263264, 47.7948738566419, 6.582312226585585, 6.226891447147788, 57.80126089261969, 42.792565534391905, 20.86151152904177, 5.3865812856304025, 88.75432634994053, 8.162209353421876, 71.14950041393638, 15.58578770292082, 8.312713125528532, 44.80347889735715, 41.46920786919761, 64.12426501695356, 23.43185598629856, 80.90459048113037, 6.659916575271539, 73.05965010832362, 63.32880390176372, 69.77479226795391, 150.14806259193728, 16.846913084055895, 42.63555828677881, 19.29871414564431, 29.871432313734346, 11.168276705846576, 87.25722845872136, 53.834056621319235, 46.165435852153536, 37.24586493886418, 14.944386441314396, 10.592862260645717, 33.65395675453615, 67.89050482372113, 14.428317350873861, 5.533356440522385, 27.024031747349987, 24.119937030166, 13.429312804892243, 60.26240745232952, 51.364845759858255, 54.6766621212557, 5.842231723744949, 18.326326624570772, 17.844100445670385, 50.79873514861414, 46.90751059176918, 22.349174142263625, 13.114789321874598, 104.13415307285484, 12.91806246828783, 75.03929809004569, 7.740382916554944, 5.740332090304366, 10.365277117312147, 38.00611466447616, 111.0554599568406, 46.668514814600194, 23.831240950249413, 59.43215443542239, 11.828104269398732, 31.262038521614144, 34.16269539010579, 49.1887419035139, 11.019000488036555, 74.24313639111887, 53.54759397848636, 20.102565845710767, 29.218945824328877, 57.5651938389466, 5.370575437886704, 12.79469428453429, 7.125373919670654, 20.55853576324594])
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);
([5176359.375, 5278197.853278245, 5280165.688778647, 5400470.240403823, 5401085.482786234, 5415904.6875, 5421084.375, 5421088.869105486, 5426535.709674103, 5457611.751633061, 5528860.9375, 5553021.875, 5554602.984113993, 5576199.208337134, 5576243.73766927, 5577168.753350776, 5584078.125, 5589565.36641185, 5605721.994601503, 5608462.117545506, 5624286.546441099, 5642694.575049502, 5643606.008462296, 5698043.505273514, 5717114.638477524, 5733592.813033938, 5755443.75, 5796592.416548447, 6674254.6875, 6698932.8125, 6778443.830400588, 6788427.690742854, 6797823.4375, 6798121.875, 6801829.319750751, 6806044.2668590825, 6806149.070278282, 6808251.5625, 6809051.5625, 6831962.848255071, 6835660.572420684, 6865896.13556999, 6869433.299348795, 6871527.8141121445, 6881832.8125, 6897929.761116174, 6909714.138681498, 6909720.109257517, 6910810.111330104, 6914277.03586287, 6920757.8125, 6924763.951216762, 6929736.086144281, 6931215.625, 6933735.9375, 6934251.162219503, 6935482.16439149, 6944615.625, 6956656.5159969255, 6975016.22892419, 6979668.042256201, 6981323.2284298055, 7000471.694890829, 7002713.87768687, 7006943.702376271, 7011939.150588079, 7015379.649322781, 7019842.473972929, 7019972.30364368, 7022898.528069976, 7025258.801919314, 7025926.342328861, 7028357.512768078, 7029591.957604919, 7032109.43066144, 7034251.5625, 7042021.489379221, 7043996.875, 7052110.9375, 7053718.64135675, 7054299.6798326075, 7057560.276968299, 7059480.162380217, 7060811.77259277, 7062554.8266791515, 7063186.948298047, 7064594.960696833, 7065919.940128694, 7069554.6875, 7069676.118294685, 7070306.413759812, 7071677.179196502, 7072155.999276852, 7075406.902509678, 7075411.784833743, 7076497.912771927, 7077178.748001068, 7077952.966489455, 7079237.673905073, 7079435.519231649, 7082707.399083979, 7086461.382597969, 7087673.872517844, 7088955.940033788, 7089304.323921151, 7090448.405209036, 7090507.1838376885, 7092968.75, 7093257.8125, 7094675.0, 7099199.818232614, 7106135.854509128, 7111083.314558238, 7115721.303995722, 7118423.677889852, 7123196.875, 7123208.406824104, 7127319.811735158, 7131226.096245381, 7131531.772857061, 7136480.382267867, 7137679.6875, 7141068.807337026, 7143090.625, 7153533.3499095, 7153923.904151928, 7166920.3125, 7168064.0625, 7174831.045751801, 7177429.08311779, 7178807.046372351, 7179002.052617889, 7179455.518584996, 7180581.900047146, 7181371.875, 7181386.889180274, 7184402.406058726, 7184737.5, 7186859.375, 7194644.2726303805, 7198633.841741391, 7204351.553327817, 7209897.522612095, 7214095.3125, 7215396.318173606, 7217396.377757388, 7220332.055667191, 7231883.6093235295, 7244774.711913046, 7248542.446329602, 7260450.163488039, 7272235.935841728, 7273263.680975504, 7277593.75, 7283084.375, 7296131.25, 7297169.541810198, 7311802.158685105, 7316646.78175492, 7318259.375, 7320298.637116808, 7322367.630368029, 7322416.177918096, 7322933.77513504, 7324717.896715794, 7330125.8831184255, 7334932.654152049, 7341517.1875, 7341560.9375, 7343160.9375, 7360978.702128716, 7366147.549165035, 7370115.919074487, 7374157.101263578, 7383192.030639198, 7383194.673123452, 7385265.625, 7396125.443678394, 7403463.434886904, 7410735.4358977135, 7417703.125, 7417775.572758404, 7419913.934088181, 7425251.870221939, 7429350.0, 7433251.842459021, 7433646.076638063, 7433744.028405825, 7445418.598797604, 7449007.8125, 7451329.410092838, 7455553.595693597, 7457114.222968031, 7461385.9375, 7475319.68637211, 7481364.0625, 7481763.982619709, 7481813.865356766, 7482945.922651834, 7484432.059944403, 7485832.222239301, 7487999.924963288, 7491304.563842744, 7492060.786520454, 7499286.03838995, 7518776.15149659, 7518847.394523694, 7531252.610251022, 7534707.595795229, 7539795.765466375, 7551064.858113601, 7568617.1875, 7575163.949322625, 7582510.202324729, 7585976.190737988, 7588896.875, 7591698.788442281, 7643392.13328214, 7646962.113922395, 7649611.20102011, 7671263.345297058, 7691327.42148156, 7710009.7152959, 7718834.78317065, 7746721.604688283, 7757493.601930725, 7798804.389315552, 7826642.1875, 8013757.956788702, 8077939.0625, 8130068.75, 8151998.4375, 8152100.821255108, 8305703.125, 8415407.628299989, 8896535.23400131, 9789375.117735244, 9852244.69609916, 10405777.701116268, 10467177.466736976, 10520561.20661951, 10564724.606593326, 10579152.649226941, 10605637.947142633, 10654496.875, 10660319.040646847, 10699251.216312686, 10769570.542770673, 10815071.875, 10859384.487486435, 10950733.825077353, 10966074.20728845, 11023655.363445703, 11023715.16195357, 11040703.955105523, 11125486.327487502, 11125498.92967076, 11144811.185865572, 11147740.310647408, 11162409.787489574, 11210626.687112886, 11215238.07866137, 11263245.3125, 11284743.452293929, 11313321.858444536, 11325637.5, 11359357.45875484, 11375595.3125, 11458699.151188212, 11466123.110315347, 11509106.319945596, 11521186.707100084, 11532137.817812612, 11548196.593452753, 11551781.25, 11588440.49719893, 11611078.125, 11613115.625, 11664132.315269183, 11677733.127305765, 11678416.712696565, 11704965.625, 11734585.9375, 11742937.5, 11760975.095458027, 11791107.182924503, 11819824.554404428, 11821983.946317345, 11834204.6875, 11876764.65200226, 11878132.8125, 11896499.105642298, 11908143.10512788, 11911215.625, 11934916.84448819, 11938329.6875, 11980809.83617919, 12002333.511217179, 12006390.098545287, 12075260.9375, 12085235.9375, 12098443.331469974, 12123302.852603862, 12127043.173896488, 12178870.005816469, 12181725.60153862, 12226050.670633523, 12236124.731765738, 12248092.163772088, 12292051.78613527, 12303845.887398228, 12324060.170074416, 12344226.518092453, 12353876.5625, 12371062.958381189, 12425622.78022515, 12497515.135917256, 12501923.430186063, 12547076.634274248, 12549834.273760656, 12557903.607911628, 12587467.117977092, 12611262.672294455, 12669396.564752629, 12671462.5, 12671571.542608751, 12671679.6875, 12680839.0625, 12696104.6875, 12751219.675828837, 12766687.202837577, 12776142.237474922, 12778873.4375, 12826359.960741932, 12864698.4375, 12904932.8125, 12906733.581488125, 12910429.6875, 12912159.787484478, 12923238.818687357, 12995738.796300104, 13053302.653925, 13066264.807913376, 13105453.880543645, 13134266.018862609, 13143799.737385549, 13183479.222262494, 13187346.401760308, 13192376.5625, 13208086.42013418, 13212172.933846084, 13220994.886442848, 13222825.784947636, 13237331.31228769, 13237516.411811348, 13256157.520928895, 13280460.252598384, 13297696.531420432, 13303355.648557529, 13326878.44834986, 13332312.5, 13336596.636951024, 13342112.14194723, 13374319.610202905, 13374331.545349503, 13385670.683836093, 13393732.777929438, 13417545.586610572, 13417989.48557362, 13419925.36773336, 13431107.156304933, 13457082.835724477, 13457948.455262294, 13460867.06422994, 13468678.920015655, 13469367.1875, 13479339.64240139, 13496035.9375, 13498180.081452075, 13498304.393954448, 13503491.208253877, 13517572.431725089, 13525687.417931704, 13551563.193934044, 13554934.612031162, 13555589.0625, 13556738.831478262, 13566140.834861567, 13567221.875, 13567562.5, 13589725.133138914, 13610787.5, 13620124.955963738, 13622606.991826974, 13623079.667350292, 13623342.761467515, 13625411.163905488, 13626034.515583176, 13626268.198861016, 13629561.960817927, 13673682.871533323, 13684306.250068264, 13687718.648866095, 13697279.505606005, 13705484.375, 13709692.576316193, 13714731.25, 13714959.905672263, 13760138.161014616, 13763331.122411804, 13790038.820078395, 13844170.786839887, 13844212.099482466, 13851028.059902573, 13861474.230245875, 13866315.41452956, 13874507.096958136, 13878731.097244963, 13885627.857059715, 13891667.919879634, 13897539.0625, 13911759.65099801, 13914119.696912253, 13914785.603167066, 13931431.76309551, 13931579.257231202, 13938797.45668302, 13947502.799169827, 13949985.9375, 13958570.3125, 13986615.876107745, 14014105.935385259, 14028500.0, 14039227.882173218, 14047759.375, 14067078.284353012, 14067671.924098322, 14070810.9375, 14099833.348685153, 14105589.0625, 14131776.885052864, 14137525.0, 14151551.020488773, 14167248.4375, 14168974.754252618, 14342466.294668982, 14373752.398337016, 14379141.647287978, 14427713.974990902, 14457939.87276123, 14470221.875, 14477551.588301234, 14581937.5, 14741059.375, 14792668.502805809, 14798075.0, 14798323.438254615, 14824510.566071574, 14833496.875, 14841736.615851704, 14845901.421011275, 14875315.63068344, 14887735.735675896, 14945123.165863588, 15004876.166372702, 15004889.0625, 15709889.0625, 17069522.929147452, 17299114.43715722, 17314383.625766963, 17361295.910388537, 17367556.019878175, 17387342.1875, 17395277.70673881, 17398797.85850792, 17407473.228699677, 17414694.656451683, 17416327.638028655, 17420373.801409166, 17422911.452424914, 17423957.26303541, 17424162.5, 17427920.263256446, 17442637.40244545, 17446295.018775973, 17448063.02536805, 17468964.468700606, 17484452.14979487, 17486437.43446301, 17490357.8125, 17531687.429321755, 17538085.58961876, 17538571.68998395, 17539127.473311298, 17551430.671352863, 17555147.580319773, 17555221.125538804, 17560066.481300715, 17638278.13989996, 336124202.52591705], [52.184091793816975, 8.000041187308591, 9.385920209364965, 77.21459060921242, 13.451809496003616, 43.98308450743192, 60.23730819009052, 25.21519704467307, 64.7237789546834, 7.926865929694985, 78.60484411621515, 94.81025112142927, 26.020857607538858, 75.07670718512074, 6.761047328106072, 6.764175709236198, 71.8087308575958, 26.053869816392087, 23.891888445962778, 22.751698467997162, 17.912791088750918, 72.83497014239761, 9.135191989141818, 10.765580431748635, 14.853145766194656, 17.469459759503238, 49.13899516612646, 15.892850198626467, 38.7274022119345, 30.560074181090208, 36.88767128470816, 68.76682116345336, 38.69266575111246, 34.81017535350598, 37.87378118530908, 25.512392810314353, 8.652170698578146, 58.655004228252636, 44.47235895041948, 23.918681892455595, 8.10251314853253, 55.3568096164301, 5.077857726426501, 23.61265360079565, 41.871087723345184, 35.906487138577276, 42.573057911584826, 69.7804047375858, 12.513487384248945, 35.35152305286264, 68.3518035472774, 26.50308289151804, 7.657899339515086, 41.94752465297758, 43.16289641100993, 5.383389028934666, 7.971208992949851, 74.35215335918157, 6.817798672226604, 19.537419597565545, 47.26980633224066, 61.5424491779041, 61.8605849044191, 8.191204680954392, 8.261552645615602, 43.93884734552492, 8.266491230969137, 5.253234450333963, 9.722202254776686, 31.9580250930399, 7.40833052542899, 9.151925285738606, 41.27927267039169, 20.651467407546853, 98.69326749268916, 53.940736265008965, 58.323210022053516, 48.694891030978816, 61.66753433040779, 20.568140004217593, 39.001284184421394, 20.02842232641656, 17.124030170985115, 33.03286184746271, 8.361574945300097, 54.03135260434296, 6.223750749466156, 32.3130619137347, 61.31263869750808, 22.068343457244712, 15.604719001519141, 8.642043201307008, 49.213885642234985, 7.857615874556166, 74.18758433981594, 44.619238477536314, 8.646007830935265, 10.888665373703887, 61.74781857076299, 97.99409051793238, 49.01627555413962, 11.902828525069237, 5.73501456777269, 7.494847833514609, 7.029048778422687, 11.077606293895254, 52.743806836157226, 81.13299197501658, 32.80565299543323, 31.761954000640998, 10.855001284849589, 5.609469325706591, 12.014982295037415, 18.236332756694406, 57.4427381877077, 43.104259973214944, 5.992889706771412, 20.12548812585788, 30.19946596342023, 12.056767990520193, 25.256964455685118, 34.44416496309503, 60.25892926598113, 50.356589879113265, 7.809905893993661, 6.860388754497189, 38.54637552048089, 91.88919376986767, 5.52324916898747, 13.048978691046507, 8.757967799445296, 17.209035291686625, 7.358514322706778, 90.08050797698263, 76.41467856161766, 41.83278485501796, 71.43122507136941, 53.13854357297615, 42.92728499135152, 94.97302282927551, 64.8076135989443, 58.50915208481255, 24.774067079441593, 54.48480224663971, 10.35622240368168, 41.25627812221201, 36.22668018989738, 12.65290490414397, 49.34417320040431, 16.02923147841358, 15.60354964030981, 18.06206820529519, 13.66267288753798, 71.70025475295485, 36.747746656235, 44.438369402779124, 7.917139669181635, 12.207786136474148, 14.597627270569857, 37.43312247211789, 15.761159513684962, 9.736221770882924, 8.602342486289114, 13.060615042411076, 45.02799979352094, 48.30361665740414, 9.201006345892084, 58.14908977436397, 58.59782787130499, 46.526617203023136, 8.617326574934927, 15.659634155054272, 11.620733668318712, 24.14285253719602, 7.547472258387531, 41.891624092312604, 34.85510991184974, 38.15685175688678, 18.344434210221618, 58.818678487691294, 53.92573661049925, 12.158900893640647, 14.473655982495778, 7.204762095808517, 47.94566211746255, 6.820715618970731, 10.073520186051885, 5.0462670083552394, 56.566095771800846, 49.03081850929036, 24.04368849379824, 60.764844750655875, 15.590571730825879, 38.25927975787228, 10.272298962603921, 38.24395100063418, 14.64554024552145, 6.456625077138597, 47.508720384112486, 19.29631078394827, 13.085459023507754, 21.1919615338223, 10.939771290777413, 23.127726029199525, 5.537997166782251, 7.017446017268324, 5.307409338260002, 8.931040135257902, 32.16699490617025, 9.746219083037845, 55.3375939693899, 46.2642658228088, 35.207911056264216, 29.543076148661687, 10.59241976783132, 59.79628938965598, 11.72445214287882, 47.35757002528918, 23.021690707204293, 11.56620064749119, 12.768099004114301, 6.948033338091945, 61.19772817297196, 22.172583575672864, 14.002439573958537, 8.209797758167843, 13.914981745130515, 49.72575065726217, 19.445830219808148, 60.428355480353936, 56.86217398186792, 67.60956231710293, 5.0292699982072415, 59.006292352770764, 7.537947171963403, 20.325167496226612, 15.410088770745258, 5.887464941993215, 12.338858046378842, 20.047377305069343, 81.86918185981216, 132.54048527442228, 28.040579601643188, 16.972140463292074, 73.56137613476818, 6.118854592998611, 23.178162690019924, 11.507777730828455, 103.05126997654735, 61.186759651763346, 134.48507149223772, 18.997824839564935, 18.344649820215526, 5.033826140470363, 5.761706479540241, 16.670267532491412, 28.013217893018993, 9.55062717581921, 9.079029987607063, 89.32639157372306, 34.79003542728762, 6.343185267566669, 57.1566900297415, 54.9977730858003, 61.317479699017554, 106.84170565411802, 16.505296361276468, 64.0709718739057, 18.600225466707236, 29.13347960370785, 23.928186152365367, 48.00395318223027, 8.462110141758568, 10.887202007300303, 83.48217031311705, 29.042827752447042, 59.93287465100778, 36.845837363068775, 66.86863134167604, 5.9672977243176355, 29.268332280668442, 59.69843478290925, 47.58290269258862, 43.645542405509644, 19.980506041524862, 12.98464637556844, 22.053541699749136, 22.235305407442336, 59.0851599454004, 8.625025349713995, 56.05223598739221, 12.361836912108249, 55.10758371821497, 89.26059860528872, 13.777596267339112, 87.03942224782665, 17.048550292954346, 33.651029682766506, 33.189481897577565, 101.33393949041167, 94.48584221926666, 52.93001693308808, 9.11122933780715, 40.84896474373514, 8.130793802350777, 6.889850494537435, 31.227424900259454, 41.94995628880258, 9.311420150305832, 19.85847316123365, 28.181087768076807, 19.184373300162132, 140.97667171733235, 43.2898869379084, 25.449554899675178, 12.930072575517297, 55.8284285755228, 6.485945317900089, 57.50123456354396, 26.69638656594714, 15.130446122025557, 89.25411227497271, 17.004226417821275, 7.494752782109766, 75.03805109188852, 7.487284513900317, 56.332493027363924, 146.24187644762105, 44.423127282205684, 38.48555775136522, 48.318903598690234, 27.38251307416874, 126.97836279159674, 5.8586297845842354, 29.193796181336754, 64.77422062014851, 33.27478271396353, 65.30178166897642, 11.840075859156377, 34.0615166126331, 5.520299008501564, 157.07416641740588, 18.822380728519843, 12.557259039282332, 7.228983636397862, 24.24033307126636, 11.79015976871047, 76.70270523874484, 64.66575186249148, 12.491081370478739, 94.81221160865641, 22.26385672636749, 10.096826970464756, 107.53764764441608, 81.24905468655976, 10.700622679252872, 19.827657186761947, 39.89251460109547, 46.88683872330879, 16.336029584743645, 38.477600565063284, 120.05844374135889, 68.75158764823911, 19.68746293387615, 26.677591950366313, 7.743301881493542, 8.268380052414974, 40.96479151929429, 24.983465745515268, 7.741885293083955, 27.66861467558177, 25.272940123516097, 15.0096459394469, 21.671177740293412, 99.40908313834669, 74.79435617150136, 5.48721949363613, 101.33053217283268, 16.690740662212026, 49.36982611587534, 12.180726844509678, 30.00738740003372, 20.422881777833773, 93.14777401499248, 6.524901161546767, 171.74241681634726, 14.128469658645328, 5.543353218247235, 47.167758694575284, 46.30041241425348, 14.61841245994605, 45.749860262099276, 39.48385026730473, 61.61657488015735, 38.24385519010955, 7.280972126193216, 31.549198138622124, 8.362346024784586, 88.40424285992091, 22.426877267652724, 16.77422040981621, 14.941769621837086, 82.57119249341278, 7.538748413944334, 57.09653881993552, 7.750191786431132, 62.02152536110484, 17.440878204648627, 73.08682067420676, 14.476525187934236, 32.892847334636016, 18.31049395518196, 27.87418301483283, 20.375755861403228, 10.105609876396027, 40.884518279347944, 26.604967154373206, 21.62200747307059, 7.509120195767349, 67.32107298197641, 70.47724724423693, 24.81604035344633, 7.394796111771666, 53.328115998419875, 47.546614720263264, 47.7948738566419, 6.582312226585585, 6.226891447147788, 57.80126089261969, 42.792565534391905, 20.86151152904177, 5.3865812856304025, 88.75432634994053, 8.162209353421876, 71.14950041393638, 15.58578770292082, 8.312713125528532, 44.80347889735715, 41.46920786919761, 64.12426501695356, 23.43185598629856, 80.90459048113037, 6.659916575271539, 73.05965010832362, 63.32880390176372, 69.77479226795391, 150.14806259193728, 16.846913084055895, 42.63555828677881, 19.29871414564431, 29.871432313734346, 11.168276705846576, 87.25722845872136, 53.834056621319235, 46.165435852153536, 37.24586493886418, 14.944386441314396, 10.592862260645717, 33.65395675453615, 67.89050482372113, 14.428317350873861, 5.533356440522385, 27.024031747349987, 24.119937030166, 13.429312804892243, 60.26240745232952, 51.364845759858255, 54.6766621212557, 5.842231723744949, 18.326326624570772, 17.844100445670385, 50.79873514861414, 46.90751059176918, 22.349174142263625, 13.114789321874598, 104.13415307285484, 12.91806246828783, 75.03929809004569, 7.740382916554944, 5.740332090304366, 10.365277117312147, 38.00611466447616, 111.0554599568406, 46.668514814600194, 23.831240950249413, 59.43215443542239, 11.828104269398732, 31.262038521614144, 34.16269539010579, 49.1887419035139, 11.019000488036555, 74.24313639111887, 53.54759397848636, 20.102565845710767, 29.218945824328877, 57.5651938389466, 5.370575437886704, 12.79469428453429, 7.125373919670654, 20.55853576324594])
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)