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 = 44856
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);
([7511994.438764137, 7524714.251180007, 7663512.596040253, 7670021.247671681, 7816598.143085433, 7835769.7772962265, 7979982.3110476835, 7980125.960075133, 7982791.762461589, 7986032.438756305, 8020807.401621959, 8022799.136136739, 8023107.8125, 8025905.397175101, 8026915.257017463, 8030090.138173257, 8031103.540254865, 8089049.0718760025, 8249746.756191713, 8252264.997113337, 8255072.922176732, 8255892.1875, 8256577.1465149755, 8256587.5256487, 8257142.236506797, 8257330.554292351, 8257830.640800148, 8258163.460169264, 8258788.018673835, 8259104.374610131, 8260189.874552345, 8272926.5625, 8273277.291004267, 8274578.338574811, 8277908.362797634, 8336421.486539121, 8337612.877495819, 8350759.618045778, 8354312.834384247, 8357745.624854878, 8359827.717854511, 8380065.7690770095, 8381441.966595228, 8381459.7542880485, 8409574.767731868, 8411241.934179636, 8416237.77967936, 8416725.0, 8416920.98479833, 8416935.228550505, 8417626.295036018, 8417672.389763096, 8417789.001593085, 8417907.264953557, 8418300.714582127, 8418307.560962096, 8418489.267836154, 8418841.063777242, 8425861.971529922, 8428445.191830635, 8428943.75, 8429735.9375, 8429813.192020077, 8431456.25, 8434030.264133275, 8464554.6875, 8534005.093962153, 8534419.213517658, 8534453.125, 8536866.226679666, 8542094.566427754, 8544471.079070464, 8546908.624879142, 8546985.9375, 8547646.07694954, 8548016.332303612, 8548225.0, 8548232.8125, 8549240.625, 8549262.5, 8549818.75, 8550395.131995067, 8550404.665383007, 8550463.17534325, 8550512.5, 8551022.808437543, 8551351.5625, 8552488.708679345, 8552879.035034383, 8553831.150416285, 8554135.9375, 8562685.284997128, 8569036.82496311, 8578095.16853073, 8579261.997269824, 8579910.197564537, 8595372.144737542, 8595949.49510039, 8596462.5, 8602854.757026348, 8607884.456879377, 8608018.75, 8609402.464673292, 8610669.171848688, 8611558.92648298, 8611620.75318187, 8612290.406779796, 8612315.625, 8612351.334997073, 8612706.499267103, 8613188.515211971, 8613199.616186924, 8613297.90738576, 8613362.802937884, 8613363.357564364, 8613600.0, 8613717.1875, 8613738.28238248, 8613742.412870677, 8613786.90719506, 8613787.3017975, 8613854.6875, 8613910.084709782, 8613999.113649476, 8614343.75, 8614384.279243175, 8614525.97120216, 8614663.381494949, 8614782.8125, 8614841.321152367, 8614932.807839878, 8614949.143193642, 8614965.620116746, 8615093.75, 8615103.579029128, 8615124.377665512, 8615246.693137148, 8615249.396210266, 8615262.178390488, 8615489.412217893, 8615654.6875, 8615737.218003113, 8615737.392035311, 8615758.705107851, 8615758.876098841, 8615913.076826198, 8615951.376320064, 8616007.3236781, 8616703.125, 8617410.9375, 8617871.875, 8618381.81112805, 8618576.247101543, 8618656.292360937, 8619285.654390903, 8619370.30020895, 8620042.1875, 8623062.494641583, 8623521.545157457, 8623647.938356832, 8623728.125, 8623770.3125, 8623797.566122606, 8623818.350859115, 8623856.25, 8623972.02232252, 8624457.634315496, 8624761.880636804, 8625215.768610466, 8625382.8125, 8625612.05912903, 8625792.978472287, 8625807.083330713, 8626141.964876402, 8627016.977066694, 8627043.060718704, 8628707.42769838, 8628975.0, 8629382.94238769, 8630022.894339817, 8630812.5, 8630829.196460636, 8630980.470345436, 8631461.639831116, 8631468.703714218, 8631617.1875, 8631997.224492498, 8632134.326745875, 8632258.160870727, 8632766.426331496, 8632830.550297821, 8633178.125, 8633271.42778151, 8633909.569495976, 8637116.295937605, 8638081.44803915, 8651179.293122534, 8654557.8125, 8670053.85536857, 8671216.543671478, 8676026.397675917, 8676337.811600909, 8716775.505052783, 8719102.734519605, 8720179.295401938, 8721021.047507659, 8723129.781801835, 8723145.020404292, 8723148.453775939, 8723152.066858718, 8727215.788388353, 8728272.444663884, 8730925.0, 8731470.511431323, 8731914.0625, 8734476.349598901, 8735344.740518281, 8744670.858763248, 8753516.207885277, 8758767.1875, 8761346.705210192, 8763147.577769356, 8770045.33668595, 8771447.560867937, 8771572.444676254, 8773333.355762077, 8773776.142397564, 8776197.53189379, 8776481.984622996, 8776831.157204732, 8785850.0, 8791864.613916242, 8792063.951923419, 8799517.19988214, 8805971.30203664, 8807072.313911775, 8808417.630789492, 8814299.87507118, 8818516.28262518, 8821808.826255612, 8825059.922166223, 8827521.875, 8831166.65054598, 8832341.460469697, 8832958.45050792, 8833175.65218612, 8835056.717621816, 8835717.293800538, 8837629.679743072, 8840938.888242802, 8875111.368117586, 8875329.949392375, 8877592.032420017, 8878891.71880474, 8882435.841954865, 8900099.37734584, 8902416.123285383, 8903712.045099033, 8906514.0625, 8907529.514379837, 8909401.5625, 8912454.456821257, 8922733.18986681, 8946416.692863038, 8956775.65647051, 8957194.757950138, 8957420.3125, 8958704.6875, 8959426.278988296, 8959432.8125, 8959941.695013845, 8960118.637017524, 8960146.855230782, 8960846.751190849, 8960902.74375993, 8961503.125, 8962357.410280151, 8962485.92860404, 8962497.659011362, 8963306.25, 8963481.879409714, 8964108.50242617, 8964158.117490832, 8964578.11765334, 8964843.75, 8965911.22211285, 8966096.101272533, 8966444.407634456, 8969396.875, 8969423.228384286, 8970835.910862241, 8973070.206947157, 8977778.125, 8977787.5, 8986294.165972074, 8991777.833191644, 8997665.739591582, 9010139.442712756, 9012970.146661965, 9013177.745263617, 9015066.561409254, 9016649.872575872, 9016721.36641545, 9019457.286526728, 9034545.3125, 9037081.25, 9037225.0, 9040285.497748554, 9040845.960485343, 9041738.942518681, 9041896.425513085, 9043179.6875, 9043974.748226441, 9043989.0625, 9044444.957643317, 9044461.366492078, 9044524.936310196, 9044597.52624323, 9044620.16108115, 9045406.707675476, 9045793.469659824, 9046053.125, 9046189.674164029, 9047256.249667743, 9048894.74560999, 9051094.983828166, 9051291.298377952, 9053321.63159823, 9077229.736168267, 9090318.87039474, 9090457.06581558, 9090462.122916646, 9098930.139509039, 9137588.224984389, 9149808.842375506, 9150521.712700075, 9150834.074316613, 9150918.475813331, 9151946.706362689, 9152068.75, 9152179.861265179, 9152300.996938737, 9152314.56079095, 9152678.072035365, 9153173.717268305, 9153667.1875, 9153928.63388652, 9153967.1875, 9154938.972644608, 9155072.858582472, 9155521.069265207, 9156023.006006416, 9157276.299588341, 9158581.25, 9160970.216301078, 9161585.9375, 9162461.806606093, 9163656.372816183, 9164701.5625, 9190160.456886144, 9193004.109875487, 9213958.45950062, 9219147.334570622, 9220453.68267678, 9279674.874929326, 9287127.530551728, 9288202.537971891, 9300590.037887314, 9328296.875, 9334680.584207732, 9359992.1875, 9360829.986859873, 9398759.938939068, 9429241.626658764, 9429321.206973096, 9438687.589304954, 9521836.870668586, 9523229.627508594, 9533112.5, 10512589.207012761, 10554468.75, 10563035.116365047, 10565385.100110287, 10569866.956230111, 10571418.75, 10574215.720435908, 10578849.334028428, 10589434.375, 10591551.628898447, 10591975.0, 10592326.36727791, 10592339.0625, 10598674.005347274, 10600404.390860729, 10602646.758549646, 10607150.495469848, 10625528.125, 10628012.041131504, 10628752.809736993, 10629258.577169394, 10630759.031355716, 10631774.210508008, 10631824.384319989, 10631953.587268531, 10632340.169412516, 10632424.707066022, 10632726.179140715, 10632791.076285092, 10642750.489182293, 10644779.6875, 10645085.715408614, 10650275.327836044, 10654718.75, 10654920.149894543, 10656109.375, 10656448.87847935, 10656795.124014432, 10657087.390104204, 10658310.949796386, 10658872.150514072, 10661254.561802037, 10663232.8125, 10664381.674715627, 10666398.403934073, 10670551.5625, 10676626.554226434, 10681323.3186173, 10682543.75, 10682667.15974926, 10687468.75, 10691714.845306503, 10694026.024461113, 10694103.939607162, 10694623.411319029, 10694820.513359793, 10698680.192850238, 10699429.763324285, 10702825.0, 10704954.32573644, 10706633.61852851, 10710846.031393297, 10713856.248380642, 10714017.1875, 10716558.44468772, 10728311.591001164, 10731604.99616564, 10731831.60983176, 10735890.687838878, 10751178.110617671, 10760266.065602865, 10764922.621207267, 10767170.799600445, 10767967.913947891, 10768845.93575407, 10769932.416080069, 10771372.746281223, 10771717.332403574, 10773880.578319928, 10774112.31454117, 10774263.116013166, 10774604.708625203, 10774702.091064379, 10774823.4375, 10774994.33401119, 10775859.375, 10776422.500718705, 10777193.75, 10777601.652660763, 10777831.273516225, 10778052.247836204, 10778730.638224918, 10779286.624262106, 10779577.685177404, 10780545.121239139, 10781330.489433218, 10781629.251054093, 10781898.298463125, 10782473.4375, 10783199.014287952, 10784429.688706294, 10785277.132382698, 10785499.129120575, 10786238.72452815, 10788664.95533935, 10788872.434438216, 10790507.93366308, 10790510.65950864, 10791692.1875, 10791714.0625, 10792355.692794392, 10793723.4375, 10793744.262554567, 10795698.935196437, 10795758.425931498, 10796240.625, 10796541.693910245, 10796835.9375, 10798254.400445042, 10799002.464319466, 10799591.13814658, 10802864.502320467, 10804459.651133824, 10807345.491815303, 581170913.1844199], [50.328068056951324, 60.398532092165496, 15.320657035561709, 20.31048519469123, 16.480424800203696, 46.249911070203694, 13.322070241917467, 20.528785849846457, 51.74505713274392, 22.852357191616377, 22.43072586144577, 34.77428973500821, 82.2199240427296, 6.003686589662683, 86.03030838033906, 13.362269861645197, 25.54955364357109, 11.121716843215353, 9.232924646731888, 15.521889442946064, 13.284274216525699, 48.31199053174171, 17.758855190025, 109.52055496045905, 11.205181430983417, 68.70316096948625, 20.46880853803444, 34.82855409223275, 12.048327221791629, 23.58747606258909, 7.352634710864262, 57.726188696377, 43.78546856388622, 15.82323783037916, 18.629731504504942, 5.3838283217014435, 54.18497100215185, 15.956551831550387, 15.276825090697876, 91.38280337498827, 5.9047161250799896, 15.520652659698245, 7.525247599759364, 9.67411682869759, 14.124421160320601, 24.36590261625969, 44.39451627321419, 52.331618579073805, 5.865457607012142, 28.770665852739757, 66.99133532916181, 17.348939555910775, 37.37777652939182, 100.67098633541994, 10.0187694782275, 53.14916486090236, 5.6255057334628, 5.638788464319082, 13.248743682554691, 6.0756399181230165, 45.898360323292366, 92.70363917462052, 8.689892810117302, 33.73288090489229, 12.166944757818229, 31.385450519340115, 26.390109661808303, 22.900170512326238, 43.58895438068009, 95.00218013245296, 110.7370273613008, 33.56863320256011, 90.7827999957133, 37.09394499615103, 5.872794404416293, 12.84102787030653, 32.34444828700197, 34.73604189027679, 70.6413157090247, 30.313809329650372, 64.27263338303965, 22.758611052959594, 39.868569051555035, 5.751130348026007, 37.50582646822388, 20.01700003336809, 42.85772572811644, 28.117391550252897, 20.096023801288442, 14.575058522069838, 46.70370924273643, 29.608958642396157, 33.19432060396436, 5.199769467575753, 12.294999390534905, 7.930553356205287, 9.57039407993327, 13.312647160108996, 74.34219550558633, 15.213302899810325, 14.398643294505403, 114.03935862773851, 20.030133171485712, 5.175303996585863, 5.589703807663831, 5.554857425776802, 28.647506762412824, 65.3467055531007, 14.139917450767909, 46.91735694781806, 21.03399987065619, 10.576476266085963, 74.30013284232001, 35.09785188213027, 8.798886272008794, 34.32940189557161, 32.154762480038634, 10.093943763473865, 15.876228419275016, 12.849426999191294, 69.70195497842366, 47.37104216881056, 49.12859983236488, 23.578892160780228, 176.72054025027586, 10.981260264808128, 20.628264391111077, 14.687587776155183, 45.91616596251496, 297.460874478453, 76.83657246292978, 12.615072219934028, 7.601472447751454, 46.37262506238982, 24.530627678958997, 25.15972700719692, 25.38695998121268, 12.645853630572404, 14.603298482310304, 6.835392906456927, 94.18138748180539, 7.504332119425871, 19.944481229226234, 13.677448192002867, 21.504068766646547, 12.167345680673757, 90.41265094498094, 19.55808231959865, 44.19853244359339, 71.12774074212426, 42.464985755938535, 32.70937836323238, 49.63178291385325, 41.71863255731489, 29.18748914393074, 107.7118300964306, 32.16639833073681, 67.08477316239757, 23.37688511059067, 13.32839330642796, 99.13105908517312, 37.65376161369413, 50.25448684859603, 22.561901523634127, 63.64972468687135, 47.81675566538379, 10.860164920725984, 35.33367499030955, 19.790880202348475, 46.338982716841954, 42.662263062684744, 23.22235988317157, 53.912292381915265, 24.79213488040973, 181.74083144153943, 19.792119681507224, 129.7855958779023, 47.51602161131981, 15.460823187244353, 10.373052960464223, 35.112309510735216, 6.998089780737286, 19.09356425902716, 59.21528694095759, 89.32528662966278, 116.58551787179333, 16.4649720465157, 11.05311096184837, 48.76998821501355, 6.374043225238401, 7.368362206589455, 114.50693930735278, 13.472964651265261, 29.614133711226316, 39.12011634063197, 10.684183300958876, 25.377754430653667, 59.19578296672296, 18.716173192602596, 18.193182631378633, 73.66832479984811, 91.52638152635228, 40.60241271094999, 10.56769029926467, 102.65920229793942, 45.50199260861257, 25.216113825845298, 23.852235874871756, 7.609135039109869, 31.304345033803255, 5.541063190300903, 76.26519413278717, 96.43384526889822, 8.394653728847054, 68.693664965791, 37.938817374951064, 12.558212097533055, 6.576273611332274, 7.283594897442266, 64.42058301400417, 94.08314851042282, 23.21236841542395, 11.786261320766208, 18.410947241496814, 12.152514433721123, 5.560040765634326, 22.328950745071985, 21.570666475153942, 209.33430700308458, 14.599232609453523, 38.92044408711971, 56.31404916147065, 10.959158040364345, 23.527388730255616, 66.31850964934809, 154.21640151756475, 72.89277280939717, 5.80479406964202, 162.6163435112652, 26.20002550718632, 39.321984757514585, 79.02578870823726, 75.81329465847408, 33.6477429322571, 35.75962269784276, 6.594429528058296, 17.08393186893085, 26.09821773140537, 20.282687623570943, 7.076385950159119, 72.54197635431804, 9.563120709297078, 142.83025904899597, 48.19879465918795, 32.85980949040089, 20.351008462406543, 17.255065499483596, 10.484056363360159, 52.438552969787196, 33.20242629760345, 31.771437785307846, 18.14878349498447, 6.841258349528391, 24.253478359687545, 7.859281009328642, 8.013588289153386, 44.39488987082383, 31.775169722247732, 34.755397342498256, 38.41418428551571, 10.43629037565637, 5.805667559774141, 24.172662595969804, 5.92869551513691, 10.579873512320338, 51.578317935760175, 16.014713003961347, 62.960716657993046, 19.102174182383187, 44.82301817071179, 26.555952519558026, 37.89281768863913, 36.8650805534283, 14.91358113887406, 40.27846275040789, 17.25431310626381, 12.63407636456151, 5.749930164610591, 56.62096334531406, 24.985497217290145, 20.950610529697936, 22.947379337319095, 32.289954046390626, 102.93597622338007, 8.519522212111287, 48.38285160888417, 15.397039756071834, 16.620079264970183, 10.878021435034736, 8.095480899111903, 6.966054443052523, 10.934209476172665, 10.4160576611638, 77.04325559539498, 42.55515199923024, 49.927169445301146, 53.93388257359005, 25.401285251863747, 33.41297316644591, 26.359181610946894, 48.70489248792724, 34.22188283697487, 66.80359550448303, 40.64604246272408, 102.03832231727938, 16.874686886338928, 6.0815009856004325, 33.13437983504423, 29.97330671903738, 17.029390959014986, 34.850445927015386, 42.304398936393184, 13.665085930449184, 11.127187502900522, 40.35805891392496, 98.9493376680813, 18.32759448431328, 7.517513486375431, 12.264580876272118, 244.96397401856026, 12.697661199252455, 49.598862785289455, 37.03484873258748, 12.48785815588475, 18.498591488866364, 94.00248719033476, 54.73957244057469, 24.086168028180932, 42.289056632812894, 55.334816370780054, 9.482784305450277, 20.86310682226369, 12.00746199323635, 14.74340411800854, 5.4092781060972, 51.670873584717185, 38.05039205666405, 49.59715445083617, 5.528735409908835, 63.7521678618302, 93.66959477216432, 95.29584456347949, 75.64584618605387, 36.063090301972714, 26.757147855948144, 93.02064292562606, 28.28889348114337, 8.35261025611603, 34.88797073998119, 21.688692041848917, 59.978690865156956, 28.037124871275235, 27.679223907481543, 8.605279373411676, 6.25152486225165, 20.599974877631087, 17.281170116699997, 13.074701119876817, 57.474582677452425, 193.71165548095456, 33.14114727215967, 55.368377127369065, 50.82079763155406, 20.906913002736335, 5.226825497533155, 40.45699339254274, 58.02600040027646, 11.03521346954943, 83.72766863944238, 24.396282927979392, 30.814502313676527, 56.8904158694406, 54.45756607425079, 66.09872227252241, 55.50685037831816, 6.782710357115575, 10.225099627579414, 60.092186018432386, 11.233102502368599, 41.32264670790712, 5.694445187115626, 81.50412490092096, 7.822396411937982, 10.693673436595938, 7.561655675574393, 5.4834422581810705, 74.0452144205008, 20.025439435585746, 10.668267230329064, 57.17078872555496, 7.033955946907739, 7.908988187168939, 6.026331523957529, 25.502923524652758, 26.273748863053296, 83.72357718112576, 13.657267584049574, 25.186168438106197, 6.865262415394904, 67.49957219605956, 17.139671011299555, 7.726276144955801, 35.441955692766555, 10.882712824476945, 92.42617327826977, 74.56477230101801, 6.243674706273813, 26.526210876819928, 23.622613615593274, 16.107053361617524, 48.71242901121086, 53.78562846230359, 125.24823644458476, 27.877483265500235, 68.1234387650275, 13.581350920093673, 26.37725401203678, 76.48171906550473, 9.308048072304231, 36.1823663160837, 14.149291668370731, 8.914437767200056, 5.148302087828111, 9.309363364407627, 46.29291111483302, 85.37167661380425, 6.778739303343675, 72.35142314648891, 33.27594071019088, 5.87612789531654, 47.7366205284443, 43.21878166973605, 70.52607586495826, 6.8966409729215155, 5.31170948289714, 12.951209112036125, 16.516736003490163, 25.638214993373765, 61.889713023683356, 9.73398434246546, 12.535930829581032, 17.495715326811556, 18.688750823986616, 18.002789912226195, 34.534630801109884, 10.195578023299383, 6.103298661963161, 19.935937447424422, 14.278799953186446, 17.957111043462874, 5.245736928644391, 12.068313866417643, 66.9469605922906, 12.182306430703658, 57.51899420990013, 17.17443558978755, 87.1062881202403, 11.271620224160008, 5.0989193479318375, 21.67981941118466, 31.65936503344659, 18.416952642436524, 13.493935777683754, 25.246765918266043, 120.56172827559253, 19.132531676776217, 9.216543013422678, 41.51888580370814, 6.884393078869794, 75.61687630358101, 19.587905340381027, 69.79478695240294, 32.68321589229032, 41.44812339125839, 42.843580190300486, 55.80002092303579, 8.14484540054404, 48.024214673335756, 28.640861441264715, 6.982297636945267, 89.92790162203482, 71.89747650339811, 12.006276565150916, 7.841271376016547, 55.23834838707515, 44.3629838954569, 50.34129885145923, 78.77344686120732, 12.9065168623008, 89.14749874152588, 58.025643460560666, 66.50499265557453, 61.609691808296844, 18.401232610736947])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([7511994.438764137, 7524714.251180007, 7663512.596040253, 7670021.247671681, 7816598.143085433, 7835769.7772962265, 7979982.3110476835, 7980125.960075133, 7982791.762461589, 7986032.438756305, 8020807.401621959, 8022799.136136739, 8023107.8125, 8025905.397175101, 8026915.257017463, 8030090.138173257, 8031103.540254865, 8089049.0718760025, 8249746.756191713, 8252264.997113337, 8255072.922176732, 8255892.1875, 8256577.1465149755, 8256587.5256487, 8257142.236506797, 8257330.554292351, 8257830.640800148, 8258163.460169264, 8258788.018673835, 8259104.374610131, 8260189.874552345, 8272926.5625, 8273277.291004267, 8274578.338574811, 8277908.362797634, 8336421.486539121, 8337612.877495819, 8350759.618045778, 8354312.834384247, 8357745.624854878, 8359827.717854511, 8380065.7690770095, 8381441.966595228, 8381459.7542880485, 8409574.767731868, 8411241.934179636, 8416237.77967936, 8416725.0, 8416920.98479833, 8416935.228550505, 8417626.295036018, 8417672.389763096, 8417789.001593085, 8417907.264953557, 8418300.714582127, 8418307.560962096, 8418489.267836154, 8418841.063777242, 8425861.971529922, 8428445.191830635, 8428943.75, 8429735.9375, 8429813.192020077, 8431456.25, 8434030.264133275, 8464554.6875, 8534005.093962153, 8534419.213517658, 8534453.125, 8536866.226679666, 8542094.566427754, 8544471.079070464, 8546908.624879142, 8546985.9375, 8547646.07694954, 8548016.332303612, 8548225.0, 8548232.8125, 8549240.625, 8549262.5, 8549818.75, 8550395.131995067, 8550404.665383007, 8550463.17534325, 8550512.5, 8551022.808437543, 8551351.5625, 8552488.708679345, 8552879.035034383, 8553831.150416285, 8554135.9375, 8562685.284997128, 8569036.82496311, 8578095.16853073, 8579261.997269824, 8579910.197564537, 8595372.144737542, 8595949.49510039, 8596462.5, 8602854.757026348, 8607884.456879377, 8608018.75, 8609402.464673292, 8610669.171848688, 8611558.92648298, 8611620.75318187, 8612290.406779796, 8612315.625, 8612351.334997073, 8612706.499267103, 8613188.515211971, 8613199.616186924, 8613297.90738576, 8613362.802937884, 8613363.357564364, 8613600.0, 8613717.1875, 8613738.28238248, 8613742.412870677, 8613786.90719506, 8613787.3017975, 8613854.6875, 8613910.084709782, 8613999.113649476, 8614343.75, 8614384.279243175, 8614525.97120216, 8614663.381494949, 8614782.8125, 8614841.321152367, 8614932.807839878, 8614949.143193642, 8614965.620116746, 8615093.75, 8615103.579029128, 8615124.377665512, 8615246.693137148, 8615249.396210266, 8615262.178390488, 8615489.412217893, 8615654.6875, 8615737.218003113, 8615737.392035311, 8615758.705107851, 8615758.876098841, 8615913.076826198, 8615951.376320064, 8616007.3236781, 8616703.125, 8617410.9375, 8617871.875, 8618381.81112805, 8618576.247101543, 8618656.292360937, 8619285.654390903, 8619370.30020895, 8620042.1875, 8623062.494641583, 8623521.545157457, 8623647.938356832, 8623728.125, 8623770.3125, 8623797.566122606, 8623818.350859115, 8623856.25, 8623972.02232252, 8624457.634315496, 8624761.880636804, 8625215.768610466, 8625382.8125, 8625612.05912903, 8625792.978472287, 8625807.083330713, 8626141.964876402, 8627016.977066694, 8627043.060718704, 8628707.42769838, 8628975.0, 8629382.94238769, 8630022.894339817, 8630812.5, 8630829.196460636, 8630980.470345436, 8631461.639831116, 8631468.703714218, 8631617.1875, 8631997.224492498, 8632134.326745875, 8632258.160870727, 8632766.426331496, 8632830.550297821, 8633178.125, 8633271.42778151, 8633909.569495976, 8637116.295937605, 8638081.44803915, 8651179.293122534, 8654557.8125, 8670053.85536857, 8671216.543671478, 8676026.397675917, 8676337.811600909, 8716775.505052783, 8719102.734519605, 8720179.295401938, 8721021.047507659, 8723129.781801835, 8723145.020404292, 8723148.453775939, 8723152.066858718, 8727215.788388353, 8728272.444663884, 8730925.0, 8731470.511431323, 8731914.0625, 8734476.349598901, 8735344.740518281, 8744670.858763248, 8753516.207885277, 8758767.1875, 8761346.705210192, 8763147.577769356, 8770045.33668595, 8771447.560867937, 8771572.444676254, 8773333.355762077, 8773776.142397564, 8776197.53189379, 8776481.984622996, 8776831.157204732, 8785850.0, 8791864.613916242, 8792063.951923419, 8799517.19988214, 8805971.30203664, 8807072.313911775, 8808417.630789492, 8814299.87507118, 8818516.28262518, 8821808.826255612, 8825059.922166223, 8827521.875, 8831166.65054598, 8832341.460469697, 8832958.45050792, 8833175.65218612, 8835056.717621816, 8835717.293800538, 8837629.679743072, 8840938.888242802, 8875111.368117586, 8875329.949392375, 8877592.032420017, 8878891.71880474, 8882435.841954865, 8900099.37734584, 8902416.123285383, 8903712.045099033, 8906514.0625, 8907529.514379837, 8909401.5625, 8912454.456821257, 8922733.18986681, 8946416.692863038, 8956775.65647051, 8957194.757950138, 8957420.3125, 8958704.6875, 8959426.278988296, 8959432.8125, 8959941.695013845, 8960118.637017524, 8960146.855230782, 8960846.751190849, 8960902.74375993, 8961503.125, 8962357.410280151, 8962485.92860404, 8962497.659011362, 8963306.25, 8963481.879409714, 8964108.50242617, 8964158.117490832, 8964578.11765334, 8964843.75, 8965911.22211285, 8966096.101272533, 8966444.407634456, 8969396.875, 8969423.228384286, 8970835.910862241, 8973070.206947157, 8977778.125, 8977787.5, 8986294.165972074, 8991777.833191644, 8997665.739591582, 9010139.442712756, 9012970.146661965, 9013177.745263617, 9015066.561409254, 9016649.872575872, 9016721.36641545, 9019457.286526728, 9034545.3125, 9037081.25, 9037225.0, 9040285.497748554, 9040845.960485343, 9041738.942518681, 9041896.425513085, 9043179.6875, 9043974.748226441, 9043989.0625, 9044444.957643317, 9044461.366492078, 9044524.936310196, 9044597.52624323, 9044620.16108115, 9045406.707675476, 9045793.469659824, 9046053.125, 9046189.674164029, 9047256.249667743, 9048894.74560999, 9051094.983828166, 9051291.298377952, 9053321.63159823, 9077229.736168267, 9090318.87039474, 9090457.06581558, 9090462.122916646, 9098930.139509039, 9137588.224984389, 9149808.842375506, 9150521.712700075, 9150834.074316613, 9150918.475813331, 9151946.706362689, 9152068.75, 9152179.861265179, 9152300.996938737, 9152314.56079095, 9152678.072035365, 9153173.717268305, 9153667.1875, 9153928.63388652, 9153967.1875, 9154938.972644608, 9155072.858582472, 9155521.069265207, 9156023.006006416, 9157276.299588341, 9158581.25, 9160970.216301078, 9161585.9375, 9162461.806606093, 9163656.372816183, 9164701.5625, 9190160.456886144, 9193004.109875487, 9213958.45950062, 9219147.334570622, 9220453.68267678, 9279674.874929326, 9287127.530551728, 9288202.537971891, 9300590.037887314, 9328296.875, 9334680.584207732, 9359992.1875, 9360829.986859873, 9398759.938939068, 9429241.626658764, 9429321.206973096, 9438687.589304954, 9521836.870668586, 9523229.627508594, 9533112.5, 10512589.207012761, 10554468.75, 10563035.116365047, 10565385.100110287, 10569866.956230111, 10571418.75, 10574215.720435908, 10578849.334028428, 10589434.375, 10591551.628898447, 10591975.0, 10592326.36727791, 10592339.0625, 10598674.005347274, 10600404.390860729, 10602646.758549646, 10607150.495469848, 10625528.125, 10628012.041131504, 10628752.809736993, 10629258.577169394, 10630759.031355716, 10631774.210508008, 10631824.384319989, 10631953.587268531, 10632340.169412516, 10632424.707066022, 10632726.179140715, 10632791.076285092, 10642750.489182293, 10644779.6875, 10645085.715408614, 10650275.327836044, 10654718.75, 10654920.149894543, 10656109.375, 10656448.87847935, 10656795.124014432, 10657087.390104204, 10658310.949796386, 10658872.150514072, 10661254.561802037, 10663232.8125, 10664381.674715627, 10666398.403934073, 10670551.5625, 10676626.554226434, 10681323.3186173, 10682543.75, 10682667.15974926, 10687468.75, 10691714.845306503, 10694026.024461113, 10694103.939607162, 10694623.411319029, 10694820.513359793, 10698680.192850238, 10699429.763324285, 10702825.0, 10704954.32573644, 10706633.61852851, 10710846.031393297, 10713856.248380642, 10714017.1875, 10716558.44468772, 10728311.591001164, 10731604.99616564, 10731831.60983176, 10735890.687838878, 10751178.110617671, 10760266.065602865, 10764922.621207267, 10767170.799600445, 10767967.913947891, 10768845.93575407, 10769932.416080069, 10771372.746281223, 10771717.332403574, 10773880.578319928, 10774112.31454117, 10774263.116013166, 10774604.708625203, 10774702.091064379, 10774823.4375, 10774994.33401119, 10775859.375, 10776422.500718705, 10777193.75, 10777601.652660763, 10777831.273516225, 10778052.247836204, 10778730.638224918, 10779286.624262106, 10779577.685177404, 10780545.121239139, 10781330.489433218, 10781629.251054093, 10781898.298463125, 10782473.4375, 10783199.014287952, 10784429.688706294, 10785277.132382698, 10785499.129120575, 10786238.72452815, 10788664.95533935, 10788872.434438216, 10790507.93366308, 10790510.65950864, 10791692.1875, 10791714.0625, 10792355.692794392, 10793723.4375, 10793744.262554567, 10795698.935196437, 10795758.425931498, 10796240.625, 10796541.693910245, 10796835.9375, 10798254.400445042, 10799002.464319466, 10799591.13814658, 10802864.502320467, 10804459.651133824, 10807345.491815303, 581170913.1844199], [50.328068056951324, 60.398532092165496, 15.320657035561709, 20.31048519469123, 16.480424800203696, 46.249911070203694, 13.322070241917467, 20.528785849846457, 51.74505713274392, 22.852357191616377, 22.43072586144577, 34.77428973500821, 82.2199240427296, 6.003686589662683, 86.03030838033906, 13.362269861645197, 25.54955364357109, 11.121716843215353, 9.232924646731888, 15.521889442946064, 13.284274216525699, 48.31199053174171, 17.758855190025, 109.52055496045905, 11.205181430983417, 68.70316096948625, 20.46880853803444, 34.82855409223275, 12.048327221791629, 23.58747606258909, 7.352634710864262, 57.726188696377, 43.78546856388622, 15.82323783037916, 18.629731504504942, 5.3838283217014435, 54.18497100215185, 15.956551831550387, 15.276825090697876, 91.38280337498827, 5.9047161250799896, 15.520652659698245, 7.525247599759364, 9.67411682869759, 14.124421160320601, 24.36590261625969, 44.39451627321419, 52.331618579073805, 5.865457607012142, 28.770665852739757, 66.99133532916181, 17.348939555910775, 37.37777652939182, 100.67098633541994, 10.0187694782275, 53.14916486090236, 5.6255057334628, 5.638788464319082, 13.248743682554691, 6.0756399181230165, 45.898360323292366, 92.70363917462052, 8.689892810117302, 33.73288090489229, 12.166944757818229, 31.385450519340115, 26.390109661808303, 22.900170512326238, 43.58895438068009, 95.00218013245296, 110.7370273613008, 33.56863320256011, 90.7827999957133, 37.09394499615103, 5.872794404416293, 12.84102787030653, 32.34444828700197, 34.73604189027679, 70.6413157090247, 30.313809329650372, 64.27263338303965, 22.758611052959594, 39.868569051555035, 5.751130348026007, 37.50582646822388, 20.01700003336809, 42.85772572811644, 28.117391550252897, 20.096023801288442, 14.575058522069838, 46.70370924273643, 29.608958642396157, 33.19432060396436, 5.199769467575753, 12.294999390534905, 7.930553356205287, 9.57039407993327, 13.312647160108996, 74.34219550558633, 15.213302899810325, 14.398643294505403, 114.03935862773851, 20.030133171485712, 5.175303996585863, 5.589703807663831, 5.554857425776802, 28.647506762412824, 65.3467055531007, 14.139917450767909, 46.91735694781806, 21.03399987065619, 10.576476266085963, 74.30013284232001, 35.09785188213027, 8.798886272008794, 34.32940189557161, 32.154762480038634, 10.093943763473865, 15.876228419275016, 12.849426999191294, 69.70195497842366, 47.37104216881056, 49.12859983236488, 23.578892160780228, 176.72054025027586, 10.981260264808128, 20.628264391111077, 14.687587776155183, 45.91616596251496, 297.460874478453, 76.83657246292978, 12.615072219934028, 7.601472447751454, 46.37262506238982, 24.530627678958997, 25.15972700719692, 25.38695998121268, 12.645853630572404, 14.603298482310304, 6.835392906456927, 94.18138748180539, 7.504332119425871, 19.944481229226234, 13.677448192002867, 21.504068766646547, 12.167345680673757, 90.41265094498094, 19.55808231959865, 44.19853244359339, 71.12774074212426, 42.464985755938535, 32.70937836323238, 49.63178291385325, 41.71863255731489, 29.18748914393074, 107.7118300964306, 32.16639833073681, 67.08477316239757, 23.37688511059067, 13.32839330642796, 99.13105908517312, 37.65376161369413, 50.25448684859603, 22.561901523634127, 63.64972468687135, 47.81675566538379, 10.860164920725984, 35.33367499030955, 19.790880202348475, 46.338982716841954, 42.662263062684744, 23.22235988317157, 53.912292381915265, 24.79213488040973, 181.74083144153943, 19.792119681507224, 129.7855958779023, 47.51602161131981, 15.460823187244353, 10.373052960464223, 35.112309510735216, 6.998089780737286, 19.09356425902716, 59.21528694095759, 89.32528662966278, 116.58551787179333, 16.4649720465157, 11.05311096184837, 48.76998821501355, 6.374043225238401, 7.368362206589455, 114.50693930735278, 13.472964651265261, 29.614133711226316, 39.12011634063197, 10.684183300958876, 25.377754430653667, 59.19578296672296, 18.716173192602596, 18.193182631378633, 73.66832479984811, 91.52638152635228, 40.60241271094999, 10.56769029926467, 102.65920229793942, 45.50199260861257, 25.216113825845298, 23.852235874871756, 7.609135039109869, 31.304345033803255, 5.541063190300903, 76.26519413278717, 96.43384526889822, 8.394653728847054, 68.693664965791, 37.938817374951064, 12.558212097533055, 6.576273611332274, 7.283594897442266, 64.42058301400417, 94.08314851042282, 23.21236841542395, 11.786261320766208, 18.410947241496814, 12.152514433721123, 5.560040765634326, 22.328950745071985, 21.570666475153942, 209.33430700308458, 14.599232609453523, 38.92044408711971, 56.31404916147065, 10.959158040364345, 23.527388730255616, 66.31850964934809, 154.21640151756475, 72.89277280939717, 5.80479406964202, 162.6163435112652, 26.20002550718632, 39.321984757514585, 79.02578870823726, 75.81329465847408, 33.6477429322571, 35.75962269784276, 6.594429528058296, 17.08393186893085, 26.09821773140537, 20.282687623570943, 7.076385950159119, 72.54197635431804, 9.563120709297078, 142.83025904899597, 48.19879465918795, 32.85980949040089, 20.351008462406543, 17.255065499483596, 10.484056363360159, 52.438552969787196, 33.20242629760345, 31.771437785307846, 18.14878349498447, 6.841258349528391, 24.253478359687545, 7.859281009328642, 8.013588289153386, 44.39488987082383, 31.775169722247732, 34.755397342498256, 38.41418428551571, 10.43629037565637, 5.805667559774141, 24.172662595969804, 5.92869551513691, 10.579873512320338, 51.578317935760175, 16.014713003961347, 62.960716657993046, 19.102174182383187, 44.82301817071179, 26.555952519558026, 37.89281768863913, 36.8650805534283, 14.91358113887406, 40.27846275040789, 17.25431310626381, 12.63407636456151, 5.749930164610591, 56.62096334531406, 24.985497217290145, 20.950610529697936, 22.947379337319095, 32.289954046390626, 102.93597622338007, 8.519522212111287, 48.38285160888417, 15.397039756071834, 16.620079264970183, 10.878021435034736, 8.095480899111903, 6.966054443052523, 10.934209476172665, 10.4160576611638, 77.04325559539498, 42.55515199923024, 49.927169445301146, 53.93388257359005, 25.401285251863747, 33.41297316644591, 26.359181610946894, 48.70489248792724, 34.22188283697487, 66.80359550448303, 40.64604246272408, 102.03832231727938, 16.874686886338928, 6.0815009856004325, 33.13437983504423, 29.97330671903738, 17.029390959014986, 34.850445927015386, 42.304398936393184, 13.665085930449184, 11.127187502900522, 40.35805891392496, 98.9493376680813, 18.32759448431328, 7.517513486375431, 12.264580876272118, 244.96397401856026, 12.697661199252455, 49.598862785289455, 37.03484873258748, 12.48785815588475, 18.498591488866364, 94.00248719033476, 54.73957244057469, 24.086168028180932, 42.289056632812894, 55.334816370780054, 9.482784305450277, 20.86310682226369, 12.00746199323635, 14.74340411800854, 5.4092781060972, 51.670873584717185, 38.05039205666405, 49.59715445083617, 5.528735409908835, 63.7521678618302, 93.66959477216432, 95.29584456347949, 75.64584618605387, 36.063090301972714, 26.757147855948144, 93.02064292562606, 28.28889348114337, 8.35261025611603, 34.88797073998119, 21.688692041848917, 59.978690865156956, 28.037124871275235, 27.679223907481543, 8.605279373411676, 6.25152486225165, 20.599974877631087, 17.281170116699997, 13.074701119876817, 57.474582677452425, 193.71165548095456, 33.14114727215967, 55.368377127369065, 50.82079763155406, 20.906913002736335, 5.226825497533155, 40.45699339254274, 58.02600040027646, 11.03521346954943, 83.72766863944238, 24.396282927979392, 30.814502313676527, 56.8904158694406, 54.45756607425079, 66.09872227252241, 55.50685037831816, 6.782710357115575, 10.225099627579414, 60.092186018432386, 11.233102502368599, 41.32264670790712, 5.694445187115626, 81.50412490092096, 7.822396411937982, 10.693673436595938, 7.561655675574393, 5.4834422581810705, 74.0452144205008, 20.025439435585746, 10.668267230329064, 57.17078872555496, 7.033955946907739, 7.908988187168939, 6.026331523957529, 25.502923524652758, 26.273748863053296, 83.72357718112576, 13.657267584049574, 25.186168438106197, 6.865262415394904, 67.49957219605956, 17.139671011299555, 7.726276144955801, 35.441955692766555, 10.882712824476945, 92.42617327826977, 74.56477230101801, 6.243674706273813, 26.526210876819928, 23.622613615593274, 16.107053361617524, 48.71242901121086, 53.78562846230359, 125.24823644458476, 27.877483265500235, 68.1234387650275, 13.581350920093673, 26.37725401203678, 76.48171906550473, 9.308048072304231, 36.1823663160837, 14.149291668370731, 8.914437767200056, 5.148302087828111, 9.309363364407627, 46.29291111483302, 85.37167661380425, 6.778739303343675, 72.35142314648891, 33.27594071019088, 5.87612789531654, 47.7366205284443, 43.21878166973605, 70.52607586495826, 6.8966409729215155, 5.31170948289714, 12.951209112036125, 16.516736003490163, 25.638214993373765, 61.889713023683356, 9.73398434246546, 12.535930829581032, 17.495715326811556, 18.688750823986616, 18.002789912226195, 34.534630801109884, 10.195578023299383, 6.103298661963161, 19.935937447424422, 14.278799953186446, 17.957111043462874, 5.245736928644391, 12.068313866417643, 66.9469605922906, 12.182306430703658, 57.51899420990013, 17.17443558978755, 87.1062881202403, 11.271620224160008, 5.0989193479318375, 21.67981941118466, 31.65936503344659, 18.416952642436524, 13.493935777683754, 25.246765918266043, 120.56172827559253, 19.132531676776217, 9.216543013422678, 41.51888580370814, 6.884393078869794, 75.61687630358101, 19.587905340381027, 69.79478695240294, 32.68321589229032, 41.44812339125839, 42.843580190300486, 55.80002092303579, 8.14484540054404, 48.024214673335756, 28.640861441264715, 6.982297636945267, 89.92790162203482, 71.89747650339811, 12.006276565150916, 7.841271376016547, 55.23834838707515, 44.3629838954569, 50.34129885145923, 78.77344686120732, 12.9065168623008, 89.14749874152588, 58.025643460560666, 66.50499265557453, 61.609691808296844, 18.401232610736947])
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);
([7511994.438764137, 7524714.251180007, 7663512.596040253, 7670021.247671681, 7816598.143085433, 7835769.7772962265, 7979982.3110476835, 7980125.960075133, 7982791.762461589, 7986032.438756305, 8020807.401621959, 8022799.136136739, 8023107.8125, 8025905.397175101, 8026915.257017463, 8030090.138173257, 8031103.540254865, 8089049.0718760025, 8249746.756191713, 8252264.997113337, 8255072.922176732, 8255892.1875, 8256577.1465149755, 8256587.5256487, 8257142.236506797, 8257330.554292351, 8257830.640800148, 8258163.460169264, 8258788.018673835, 8259104.374610131, 8260189.874552345, 8272926.5625, 8273277.291004267, 8274578.338574811, 8277908.362797634, 8336421.486539121, 8337612.877495819, 8350759.618045778, 8354312.834384247, 8357745.624854878, 8359827.717854511, 8380065.7690770095, 8381441.966595228, 8381459.7542880485, 8409574.767731868, 8411241.934179636, 8416237.77967936, 8416725.0, 8416920.98479833, 8416935.228550505, 8417626.295036018, 8417672.389763096, 8417789.001593085, 8417907.264953557, 8418300.714582127, 8418307.560962096, 8418489.267836154, 8418841.063777242, 8425861.971529922, 8428445.191830635, 8428943.75, 8429735.9375, 8429813.192020077, 8431456.25, 8434030.264133275, 8464554.6875, 8534005.093962153, 8534419.213517658, 8534453.125, 8536866.226679666, 8542094.566427754, 8544471.079070464, 8546908.624879142, 8546985.9375, 8547646.07694954, 8548016.332303612, 8548225.0, 8548232.8125, 8549240.625, 8549262.5, 8549818.75, 8550395.131995067, 8550404.665383007, 8550463.17534325, 8550512.5, 8551022.808437543, 8551351.5625, 8552488.708679345, 8552879.035034383, 8553831.150416285, 8554135.9375, 8562685.284997128, 8569036.82496311, 8578095.16853073, 8579261.997269824, 8579910.197564537, 8595372.144737542, 8595949.49510039, 8596462.5, 8602854.757026348, 8607884.456879377, 8608018.75, 8609402.464673292, 8610669.171848688, 8611558.92648298, 8611620.75318187, 8612290.406779796, 8612315.625, 8612351.334997073, 8612706.499267103, 8613188.515211971, 8613199.616186924, 8613297.90738576, 8613362.802937884, 8613363.357564364, 8613600.0, 8613717.1875, 8613738.28238248, 8613742.412870677, 8613786.90719506, 8613787.3017975, 8613854.6875, 8613910.084709782, 8613999.113649476, 8614343.75, 8614384.279243175, 8614525.97120216, 8614663.381494949, 8614782.8125, 8614841.321152367, 8614932.807839878, 8614949.143193642, 8614965.620116746, 8615093.75, 8615103.579029128, 8615124.377665512, 8615246.693137148, 8615249.396210266, 8615262.178390488, 8615489.412217893, 8615654.6875, 8615737.218003113, 8615737.392035311, 8615758.705107851, 8615758.876098841, 8615913.076826198, 8615951.376320064, 8616007.3236781, 8616703.125, 8617410.9375, 8617871.875, 8618381.81112805, 8618576.247101543, 8618656.292360937, 8619285.654390903, 8619370.30020895, 8620042.1875, 8623062.494641583, 8623521.545157457, 8623647.938356832, 8623728.125, 8623770.3125, 8623797.566122606, 8623818.350859115, 8623856.25, 8623972.02232252, 8624457.634315496, 8624761.880636804, 8625215.768610466, 8625382.8125, 8625612.05912903, 8625792.978472287, 8625807.083330713, 8626141.964876402, 8627016.977066694, 8627043.060718704, 8628707.42769838, 8628975.0, 8629382.94238769, 8630022.894339817, 8630812.5, 8630829.196460636, 8630980.470345436, 8631461.639831116, 8631468.703714218, 8631617.1875, 8631997.224492498, 8632134.326745875, 8632258.160870727, 8632766.426331496, 8632830.550297821, 8633178.125, 8633271.42778151, 8633909.569495976, 8637116.295937605, 8638081.44803915, 8651179.293122534, 8654557.8125, 8670053.85536857, 8671216.543671478, 8676026.397675917, 8676337.811600909, 8716775.505052783, 8719102.734519605, 8720179.295401938, 8721021.047507659, 8723129.781801835, 8723145.020404292, 8723148.453775939, 8723152.066858718, 8727215.788388353, 8728272.444663884, 8730925.0, 8731470.511431323, 8731914.0625, 8734476.349598901, 8735344.740518281, 8744670.858763248, 8753516.207885277, 8758767.1875, 8761346.705210192, 8763147.577769356, 8770045.33668595, 8771447.560867937, 8771572.444676254, 8773333.355762077, 8773776.142397564, 8776197.53189379, 8776481.984622996, 8776831.157204732, 8785850.0, 8791864.613916242, 8792063.951923419, 8799517.19988214, 8805971.30203664, 8807072.313911775, 8808417.630789492, 8814299.87507118, 8818516.28262518, 8821808.826255612, 8825059.922166223, 8827521.875, 8831166.65054598, 8832341.460469697, 8832958.45050792, 8833175.65218612, 8835056.717621816, 8835717.293800538, 8837629.679743072, 8840938.888242802, 8875111.368117586, 8875329.949392375, 8877592.032420017, 8878891.71880474, 8882435.841954865, 8900099.37734584, 8902416.123285383, 8903712.045099033, 8906514.0625, 8907529.514379837, 8909401.5625, 8912454.456821257, 8922733.18986681, 8946416.692863038, 8956775.65647051, 8957194.757950138, 8957420.3125, 8958704.6875, 8959426.278988296, 8959432.8125, 8959941.695013845, 8960118.637017524, 8960146.855230782, 8960846.751190849, 8960902.74375993, 8961503.125, 8962357.410280151, 8962485.92860404, 8962497.659011362, 8963306.25, 8963481.879409714, 8964108.50242617, 8964158.117490832, 8964578.11765334, 8964843.75, 8965911.22211285, 8966096.101272533, 8966444.407634456, 8969396.875, 8969423.228384286, 8970835.910862241, 8973070.206947157, 8977778.125, 8977787.5, 8986294.165972074, 8991777.833191644, 8997665.739591582, 9010139.442712756, 9012970.146661965, 9013177.745263617, 9015066.561409254, 9016649.872575872, 9016721.36641545, 9019457.286526728, 9034545.3125, 9037081.25, 9037225.0, 9040285.497748554, 9040845.960485343, 9041738.942518681, 9041896.425513085, 9043179.6875, 9043974.748226441, 9043989.0625, 9044444.957643317, 9044461.366492078, 9044524.936310196, 9044597.52624323, 9044620.16108115, 9045406.707675476, 9045793.469659824, 9046053.125, 9046189.674164029, 9047256.249667743, 9048894.74560999, 9051094.983828166, 9051291.298377952, 9053321.63159823, 9077229.736168267, 9090318.87039474, 9090457.06581558, 9090462.122916646, 9098930.139509039, 9137588.224984389, 9149808.842375506, 9150521.712700075, 9150834.074316613, 9150918.475813331, 9151946.706362689, 9152068.75, 9152179.861265179, 9152300.996938737, 9152314.56079095, 9152678.072035365, 9153173.717268305, 9153667.1875, 9153928.63388652, 9153967.1875, 9154938.972644608, 9155072.858582472, 9155521.069265207, 9156023.006006416, 9157276.299588341, 9158581.25, 9160970.216301078, 9161585.9375, 9162461.806606093, 9163656.372816183, 9164701.5625, 9190160.456886144, 9193004.109875487, 9213958.45950062, 9219147.334570622, 9220453.68267678, 9279674.874929326, 9287127.530551728, 9288202.537971891, 9300590.037887314, 9328296.875, 9334680.584207732, 9359992.1875, 9360829.986859873, 9398759.938939068, 9429241.626658764, 9429321.206973096, 9438687.589304954, 9521836.870668586, 9523229.627508594, 9533112.5, 10512589.207012761, 10554468.75, 10563035.116365047, 10565385.100110287, 10569866.956230111, 10571418.75, 10574215.720435908, 10578849.334028428, 10589434.375, 10591551.628898447, 10591975.0, 10592326.36727791, 10592339.0625, 10598674.005347274, 10600404.390860729, 10602646.758549646, 10607150.495469848, 10625528.125, 10628012.041131504, 10628752.809736993, 10629258.577169394, 10630759.031355716, 10631774.210508008, 10631824.384319989, 10631953.587268531, 10632340.169412516, 10632424.707066022, 10632726.179140715, 10632791.076285092, 10642750.489182293, 10644779.6875, 10645085.715408614, 10650275.327836044, 10654718.75, 10654920.149894543, 10656109.375, 10656448.87847935, 10656795.124014432, 10657087.390104204, 10658310.949796386, 10658872.150514072, 10661254.561802037, 10663232.8125, 10664381.674715627, 10666398.403934073, 10670551.5625, 10676626.554226434, 10681323.3186173, 10682543.75, 10682667.15974926, 10687468.75, 10691714.845306503, 10694026.024461113, 10694103.939607162, 10694623.411319029, 10694820.513359793, 10698680.192850238, 10699429.763324285, 10702825.0, 10704954.32573644, 10706633.61852851, 10710846.031393297, 10713856.248380642, 10714017.1875, 10716558.44468772, 10728311.591001164, 10731604.99616564, 10731831.60983176, 10735890.687838878, 10751178.110617671, 10760266.065602865, 10764922.621207267, 10767170.799600445, 10767967.913947891, 10768845.93575407, 10769932.416080069, 10771372.746281223, 10771717.332403574, 10773880.578319928, 10774112.31454117, 10774263.116013166, 10774604.708625203, 10774702.091064379, 10774823.4375, 10774994.33401119, 10775859.375, 10776422.500718705, 10777193.75, 10777601.652660763, 10777831.273516225, 10778052.247836204, 10778730.638224918, 10779286.624262106, 10779577.685177404, 10780545.121239139, 10781330.489433218, 10781629.251054093, 10781898.298463125, 10782473.4375, 10783199.014287952, 10784429.688706294, 10785277.132382698, 10785499.129120575, 10786238.72452815, 10788664.95533935, 10788872.434438216, 10790507.93366308, 10790510.65950864, 10791692.1875, 10791714.0625, 10792355.692794392, 10793723.4375, 10793744.262554567, 10795698.935196437, 10795758.425931498, 10796240.625, 10796541.693910245, 10796835.9375, 10798254.400445042, 10799002.464319466, 10799591.13814658, 10802864.502320467, 10804459.651133824, 10807345.491815303, 581170913.1844199], [50.328068056951324, 60.398532092165496, 15.320657035561709, 20.31048519469123, 16.480424800203696, 46.249911070203694, 13.322070241917467, 20.528785849846457, 51.74505713274392, 22.852357191616377, 22.43072586144577, 34.77428973500821, 82.2199240427296, 6.003686589662683, 86.03030838033906, 13.362269861645197, 25.54955364357109, 11.121716843215353, 9.232924646731888, 15.521889442946064, 13.284274216525699, 48.31199053174171, 17.758855190025, 109.52055496045905, 11.205181430983417, 68.70316096948625, 20.46880853803444, 34.82855409223275, 12.048327221791629, 23.58747606258909, 7.352634710864262, 57.726188696377, 43.78546856388622, 15.82323783037916, 18.629731504504942, 5.3838283217014435, 54.18497100215185, 15.956551831550387, 15.276825090697876, 91.38280337498827, 5.9047161250799896, 15.520652659698245, 7.525247599759364, 9.67411682869759, 14.124421160320601, 24.36590261625969, 44.39451627321419, 52.331618579073805, 5.865457607012142, 28.770665852739757, 66.99133532916181, 17.348939555910775, 37.37777652939182, 100.67098633541994, 10.0187694782275, 53.14916486090236, 5.6255057334628, 5.638788464319082, 13.248743682554691, 6.0756399181230165, 45.898360323292366, 92.70363917462052, 8.689892810117302, 33.73288090489229, 12.166944757818229, 31.385450519340115, 26.390109661808303, 22.900170512326238, 43.58895438068009, 95.00218013245296, 110.7370273613008, 33.56863320256011, 90.7827999957133, 37.09394499615103, 5.872794404416293, 12.84102787030653, 32.34444828700197, 34.73604189027679, 70.6413157090247, 30.313809329650372, 64.27263338303965, 22.758611052959594, 39.868569051555035, 5.751130348026007, 37.50582646822388, 20.01700003336809, 42.85772572811644, 28.117391550252897, 20.096023801288442, 14.575058522069838, 46.70370924273643, 29.608958642396157, 33.19432060396436, 5.199769467575753, 12.294999390534905, 7.930553356205287, 9.57039407993327, 13.312647160108996, 74.34219550558633, 15.213302899810325, 14.398643294505403, 114.03935862773851, 20.030133171485712, 5.175303996585863, 5.589703807663831, 5.554857425776802, 28.647506762412824, 65.3467055531007, 14.139917450767909, 46.91735694781806, 21.03399987065619, 10.576476266085963, 74.30013284232001, 35.09785188213027, 8.798886272008794, 34.32940189557161, 32.154762480038634, 10.093943763473865, 15.876228419275016, 12.849426999191294, 69.70195497842366, 47.37104216881056, 49.12859983236488, 23.578892160780228, 176.72054025027586, 10.981260264808128, 20.628264391111077, 14.687587776155183, 45.91616596251496, 297.460874478453, 76.83657246292978, 12.615072219934028, 7.601472447751454, 46.37262506238982, 24.530627678958997, 25.15972700719692, 25.38695998121268, 12.645853630572404, 14.603298482310304, 6.835392906456927, 94.18138748180539, 7.504332119425871, 19.944481229226234, 13.677448192002867, 21.504068766646547, 12.167345680673757, 90.41265094498094, 19.55808231959865, 44.19853244359339, 71.12774074212426, 42.464985755938535, 32.70937836323238, 49.63178291385325, 41.71863255731489, 29.18748914393074, 107.7118300964306, 32.16639833073681, 67.08477316239757, 23.37688511059067, 13.32839330642796, 99.13105908517312, 37.65376161369413, 50.25448684859603, 22.561901523634127, 63.64972468687135, 47.81675566538379, 10.860164920725984, 35.33367499030955, 19.790880202348475, 46.338982716841954, 42.662263062684744, 23.22235988317157, 53.912292381915265, 24.79213488040973, 181.74083144153943, 19.792119681507224, 129.7855958779023, 47.51602161131981, 15.460823187244353, 10.373052960464223, 35.112309510735216, 6.998089780737286, 19.09356425902716, 59.21528694095759, 89.32528662966278, 116.58551787179333, 16.4649720465157, 11.05311096184837, 48.76998821501355, 6.374043225238401, 7.368362206589455, 114.50693930735278, 13.472964651265261, 29.614133711226316, 39.12011634063197, 10.684183300958876, 25.377754430653667, 59.19578296672296, 18.716173192602596, 18.193182631378633, 73.66832479984811, 91.52638152635228, 40.60241271094999, 10.56769029926467, 102.65920229793942, 45.50199260861257, 25.216113825845298, 23.852235874871756, 7.609135039109869, 31.304345033803255, 5.541063190300903, 76.26519413278717, 96.43384526889822, 8.394653728847054, 68.693664965791, 37.938817374951064, 12.558212097533055, 6.576273611332274, 7.283594897442266, 64.42058301400417, 94.08314851042282, 23.21236841542395, 11.786261320766208, 18.410947241496814, 12.152514433721123, 5.560040765634326, 22.328950745071985, 21.570666475153942, 209.33430700308458, 14.599232609453523, 38.92044408711971, 56.31404916147065, 10.959158040364345, 23.527388730255616, 66.31850964934809, 154.21640151756475, 72.89277280939717, 5.80479406964202, 162.6163435112652, 26.20002550718632, 39.321984757514585, 79.02578870823726, 75.81329465847408, 33.6477429322571, 35.75962269784276, 6.594429528058296, 17.08393186893085, 26.09821773140537, 20.282687623570943, 7.076385950159119, 72.54197635431804, 9.563120709297078, 142.83025904899597, 48.19879465918795, 32.85980949040089, 20.351008462406543, 17.255065499483596, 10.484056363360159, 52.438552969787196, 33.20242629760345, 31.771437785307846, 18.14878349498447, 6.841258349528391, 24.253478359687545, 7.859281009328642, 8.013588289153386, 44.39488987082383, 31.775169722247732, 34.755397342498256, 38.41418428551571, 10.43629037565637, 5.805667559774141, 24.172662595969804, 5.92869551513691, 10.579873512320338, 51.578317935760175, 16.014713003961347, 62.960716657993046, 19.102174182383187, 44.82301817071179, 26.555952519558026, 37.89281768863913, 36.8650805534283, 14.91358113887406, 40.27846275040789, 17.25431310626381, 12.63407636456151, 5.749930164610591, 56.62096334531406, 24.985497217290145, 20.950610529697936, 22.947379337319095, 32.289954046390626, 102.93597622338007, 8.519522212111287, 48.38285160888417, 15.397039756071834, 16.620079264970183, 10.878021435034736, 8.095480899111903, 6.966054443052523, 10.934209476172665, 10.4160576611638, 77.04325559539498, 42.55515199923024, 49.927169445301146, 53.93388257359005, 25.401285251863747, 33.41297316644591, 26.359181610946894, 48.70489248792724, 34.22188283697487, 66.80359550448303, 40.64604246272408, 102.03832231727938, 16.874686886338928, 6.0815009856004325, 33.13437983504423, 29.97330671903738, 17.029390959014986, 34.850445927015386, 42.304398936393184, 13.665085930449184, 11.127187502900522, 40.35805891392496, 98.9493376680813, 18.32759448431328, 7.517513486375431, 12.264580876272118, 244.96397401856026, 12.697661199252455, 49.598862785289455, 37.03484873258748, 12.48785815588475, 18.498591488866364, 94.00248719033476, 54.73957244057469, 24.086168028180932, 42.289056632812894, 55.334816370780054, 9.482784305450277, 20.86310682226369, 12.00746199323635, 14.74340411800854, 5.4092781060972, 51.670873584717185, 38.05039205666405, 49.59715445083617, 5.528735409908835, 63.7521678618302, 93.66959477216432, 95.29584456347949, 75.64584618605387, 36.063090301972714, 26.757147855948144, 93.02064292562606, 28.28889348114337, 8.35261025611603, 34.88797073998119, 21.688692041848917, 59.978690865156956, 28.037124871275235, 27.679223907481543, 8.605279373411676, 6.25152486225165, 20.599974877631087, 17.281170116699997, 13.074701119876817, 57.474582677452425, 193.71165548095456, 33.14114727215967, 55.368377127369065, 50.82079763155406, 20.906913002736335, 5.226825497533155, 40.45699339254274, 58.02600040027646, 11.03521346954943, 83.72766863944238, 24.396282927979392, 30.814502313676527, 56.8904158694406, 54.45756607425079, 66.09872227252241, 55.50685037831816, 6.782710357115575, 10.225099627579414, 60.092186018432386, 11.233102502368599, 41.32264670790712, 5.694445187115626, 81.50412490092096, 7.822396411937982, 10.693673436595938, 7.561655675574393, 5.4834422581810705, 74.0452144205008, 20.025439435585746, 10.668267230329064, 57.17078872555496, 7.033955946907739, 7.908988187168939, 6.026331523957529, 25.502923524652758, 26.273748863053296, 83.72357718112576, 13.657267584049574, 25.186168438106197, 6.865262415394904, 67.49957219605956, 17.139671011299555, 7.726276144955801, 35.441955692766555, 10.882712824476945, 92.42617327826977, 74.56477230101801, 6.243674706273813, 26.526210876819928, 23.622613615593274, 16.107053361617524, 48.71242901121086, 53.78562846230359, 125.24823644458476, 27.877483265500235, 68.1234387650275, 13.581350920093673, 26.37725401203678, 76.48171906550473, 9.308048072304231, 36.1823663160837, 14.149291668370731, 8.914437767200056, 5.148302087828111, 9.309363364407627, 46.29291111483302, 85.37167661380425, 6.778739303343675, 72.35142314648891, 33.27594071019088, 5.87612789531654, 47.7366205284443, 43.21878166973605, 70.52607586495826, 6.8966409729215155, 5.31170948289714, 12.951209112036125, 16.516736003490163, 25.638214993373765, 61.889713023683356, 9.73398434246546, 12.535930829581032, 17.495715326811556, 18.688750823986616, 18.002789912226195, 34.534630801109884, 10.195578023299383, 6.103298661963161, 19.935937447424422, 14.278799953186446, 17.957111043462874, 5.245736928644391, 12.068313866417643, 66.9469605922906, 12.182306430703658, 57.51899420990013, 17.17443558978755, 87.1062881202403, 11.271620224160008, 5.0989193479318375, 21.67981941118466, 31.65936503344659, 18.416952642436524, 13.493935777683754, 25.246765918266043, 120.56172827559253, 19.132531676776217, 9.216543013422678, 41.51888580370814, 6.884393078869794, 75.61687630358101, 19.587905340381027, 69.79478695240294, 32.68321589229032, 41.44812339125839, 42.843580190300486, 55.80002092303579, 8.14484540054404, 48.024214673335756, 28.640861441264715, 6.982297636945267, 89.92790162203482, 71.89747650339811, 12.006276565150916, 7.841271376016547, 55.23834838707515, 44.3629838954569, 50.34129885145923, 78.77344686120732, 12.9065168623008, 89.14749874152588, 58.025643460560666, 66.50499265557453, 61.609691808296844, 18.401232610736947])
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)