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 = 48240
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);
([3910014.8194018146, 4365871.041709561, 4366648.189549966, 4432115.625, 4456751.305011168, 4586187.5, 4675888.037006122, 4724910.922729095, 4725035.211758849, 4725908.737295396, 4764932.8125, 4897414.818716402, 4898139.329151454, 4899593.607904949, 4962067.1875, 4993750.919991304, 4995000.0, 4996474.053462202, 5000587.600701952, 5016280.019841695, 5089785.784367688, 5111960.9375, 5117663.099922496, 5138038.515552952, 5147295.3125, 5173257.713381661, 5189553.125, 5272208.651500754, 7858878.118125062, 7998937.183893107, 7999271.875, 8008642.21425294, 8140826.5625, 8232093.029546252, 8318254.6875, 8329914.0625, 8330007.732865723, 8331063.2071464, 8333396.7204927895, 8334824.156808111, 8335331.25, 8343903.125, 8351388.481924267, 8352501.5625, 8352670.454579472, 8358210.203884542, 8426044.126635674, 8461484.221274178, 8478923.603848565, 8519017.09176515, 8519671.163631419, 8521453.50796835, 8521743.187961867, 8521950.0, 8522193.141182508, 8522765.625, 8523133.438325534, 8523599.651927931, 8524160.620010953, 8524376.351478698, 8524385.9375, 8524450.38444078, 8524485.376575157, 8526485.359116502, 8526504.6875, 8527680.192819232, 8528430.99899828, 8528771.33624739, 8528825.44824642, 8535415.59023353, 8537355.384617068, 8538428.125, 8540335.584436208, 8542365.486523516, 8543185.396108765, 8545460.67703779, 8562026.5625, 8566509.375, 8569592.599586722, 8580070.66442481, 8625198.46348207, 8682684.375, 8688032.8125, 8689419.988341711, 8702893.325506413, 8711207.59329776, 8714309.800957955, 8715674.778663142, 8716381.555926824, 8716432.352189507, 8716797.244935399, 8717168.47455672, 8717377.700179731, 8717397.139809752, 8719137.388298994, 8719156.25, 8719828.350238103, 8719842.100319996, 8720666.233044147, 8721428.835282575, 8723305.145912651, 8731212.5, 8743274.07983356, 8745881.16577824, 8746935.789313221, 8748029.549046723, 8748348.4375, 8748386.886841381, 8748434.27611158, 8752620.3125, 8752626.451647926, 8780590.625, 8785692.321825508, 8811787.5, 8833841.627961911, 8834132.976984683, 8834940.625, 8843214.0625, 8843798.4375, 8843978.639688369, 8845145.792431833, 8850379.226035608, 8903631.25, 8906150.0, 8907231.25, 8907286.394763708, 8922509.375, 8923197.898707053, 8932870.3125, 8934614.0625, 8935296.593524855, 8935996.206659578, 8936816.290710649, 8936818.94084929, 8938843.452153, 8939674.78728012, 8939681.217692744, 8940223.4375, 8940226.465782732, 8940475.0, 8940801.541903665, 8940842.162362223, 8940964.569668746, 8941298.116985053, 8941391.42301674, 8941741.547489034, 8942217.1875, 8942303.125, 8942547.782576567, 8942575.720424745, 8942590.106948173, 8942618.75, 8943098.336435124, 8943220.392540742, 8944040.044485684, 8945493.75, 8953315.625, 8954069.196523713, 8954075.931590796, 8954633.512970926, 8955016.118982531, 8955032.206794204, 8957075.752796663, 8957185.515313458, 8957507.414243396, 8958693.957597852, 8960132.8125, 8976781.660278516, 8976805.114621326, 8978153.915198157, 8978576.529812444, 8979989.0625, 8979995.3125, 8980849.119752444, 8981021.250518287, 8981613.940275062, 8983196.367137792, 8983217.471941829, 8988609.371579766, 8990137.967347275, 8995624.655706385, 8995888.105106452, 8996627.04051476, 9001272.497743383, 9007156.978826828, 9008153.125, 9008340.729593217, 9008618.75, 9008909.706133831, 9009090.625, 9009458.477815984, 9009523.614459397, 9009653.125, 9009974.365987148, 9012239.819104508, 9012655.364676045, 9013360.340155289, 9015226.5625, 9015297.410609175, 9015531.25, 9015546.84275286, 9016079.17015203, 9016456.25, 9016796.875, 9016939.0625, 9018065.47048876, 9019007.598853597, 9019424.781650782, 9020216.204898566, 9020700.219281388, 9021750.0, 9022009.375, 9022039.07890427, 9022117.054477932, 9022137.428465886, 9022166.893257631, 9022683.500176037, 9023096.11375978, 9023111.215038115, 9023174.094177118, 9023442.564468732, 9023861.97372736, 9024143.75, 9025093.75, 9025234.138274176, 9025281.25, 9025393.75, 9025399.585635705, 9025634.375, 9025823.4375, 9025824.3051855, 9025962.5, 9026859.829160407, 9027274.890437229, 9027672.454261405, 9027911.069156326, 9028061.852266202, 9028183.787916828, 9028326.5625, 9028390.098379996, 9028460.9375, 9029149.50941549, 9029163.729770934, 9029282.072274571, 9029302.516967375, 9029634.375, 9029941.853473887, 9030044.33999454, 9030186.37782236, 9030281.25, 9030500.0, 9030503.084601222, 9030520.997055272, 9031378.303729942, 9032110.656052893, 9032149.973909745, 9033028.125, 9033934.487218443, 9034633.83690556, 9034652.957940388, 9034899.967225933, 9035402.147926139, 9036098.50580103, 9036178.125, 9038883.443018883, 9039474.471852457, 9040592.1875, 9041704.594084738, 9042593.75, 9046849.481746728, 9051306.141803468, 9052698.150654785, 9052714.877641186, 9052753.570555959, 9056847.325566074, 9062510.9375, 9063555.71338701, 9064218.75, 9065071.132197255, 9065638.770624513, 9066702.78622392, 9067486.701920865, 9067615.324827421, 9068552.479030035, 9069988.11945227, 9069993.731729507, 9071403.720326673, 9077095.218276614, 9077495.3125, 9079097.853132205, 9079171.875, 9079740.625, 9083656.25, 9086332.8125, 9086332.8125, 9086787.12871061, 9086810.9375, 9089397.075729573, 9089400.940688163, 9090137.319029123, 9090449.658358864, 9091783.114446089, 9092626.502024684, 9093303.125, 9093757.052094826, 9097567.1875, 9100103.428038187, 9101284.339614797, 9102750.957168976, 9105148.592938282, 9105280.947507199, 9107571.274942506, 9110897.773587795, 9116896.875, 9177286.006820817, 9178843.02244361, 9179099.066186946, 9179265.625, 9179705.40038077, 9184450.984484904, 9187837.109803181, 9198187.346512605, 9199600.0, 9208436.538005788, 9215647.58424489, 9215656.349891137, 9220460.592340775, 9233321.095769065, 9243669.121055583, 9255972.01304002, 9256781.524706552, 9256958.579869783, 9261357.196251675, 9265636.392539795, 9265651.5625, 9266696.8470426, 9268285.614791855, 9271367.805880785, 9283673.4375, 9298368.75, 9298392.857950756, 9313939.0625, 9318899.676206494, 9319383.217727654, 9320184.679209884, 9322546.172988059, 9326092.102253003, 9327159.375, 9327232.8125, 9331745.418387065, 9332896.875, 9335329.6875, 9342043.011454243, 9342114.0625, 9342346.758064944, 9342512.5, 9342928.125, 9344628.054153837, 9347010.9375, 9347018.26855879, 9347328.438777333, 9348664.0625, 9348967.261186296, 9349613.654678162, 9351163.772366062, 9353159.709445283, 9353695.130763272, 9354211.458132006, 9355114.0625, 9359686.29699372, 9370001.5625, 9375216.905901369, 9376018.75, 9378689.21165212, 9385837.553814525, 9386630.004669448, 9390265.720974179, 9391996.875, 9392456.093642805, 9396526.980268562, 9396584.330165122, 9397627.740121828, 9399044.44253145, 9400444.377489066, 9402468.75, 9404839.0625, 9405517.1875, 9405838.55837837, 9411281.41896581, 9411321.021800742, 9414594.16525168, 9416834.197788605, 9417010.9375, 9417518.75, 9417782.604005229, 9418923.233139755, 9419249.801805234, 9420269.086537804, 9421682.8125, 9422029.560480203, 9422380.041298838, 9422544.85604374, 9422851.309094522, 9428275.925191343, 9439181.094588943, 9446658.828845475, 9454511.05135791, 9454546.430704959, 9455671.378693778, 9456116.39656712, 9456397.468078075, 9456695.3125, 9456853.125, 9457146.77241501, 9459726.213059133, 9460774.386727327, 9461456.25, 9462055.985219043, 9462303.125, 9463882.8125, 9465464.021324731, 9492736.819528794, 9492779.333749047, 9508375.550402498, 9509081.856690725, 9513296.875, 9517017.754068272, 9518942.87787002, 9526969.779341122, 9534720.3125, 9555951.343560228, 9562494.045772677, 9567825.112379488, 9568079.286456928, 9568466.94948145, 9570962.5, 9570974.599218486, 9571044.70856867, 9572650.188685182, 9573213.878311653, 9573415.625, 9576210.83247661, 9583905.22582162, 9586254.638998972, 9587317.1875, 9587349.22339871, 9587553.392005583, 9589885.9375, 9590143.75, 9593087.5, 9603923.650641385, 9603982.965480449, 9604390.680842701, 9604546.875, 9606937.238704486, 9609308.850961689, 9609516.284902824, 9611308.823663715, 9613093.75, 9624877.759948839, 9628488.58698483, 9630603.673003042, 9630767.874217683, 9632660.9375, 9632876.838467829, 9652425.38970571, 9665323.025330754, 9679040.569819259, 9680973.498929648, 9681725.575091952, 9684459.375, 9684514.672573822, 9684690.628494522, 9685429.077217486, 9685441.636011573, 9686747.68076135, 9689137.5, 9689807.00624082, 9689820.919587554, 9691253.029709218, 9692564.080210293, 9692689.781769462, 9695335.059501812, 9695995.892970309, 9696205.681000935, 9696382.8125, 9696424.074754816, 9698119.834255315, 9698979.514743427, 9699112.510830484, 9699218.226893734, 9701579.83995475, 9702131.217521526, 9705942.1875, 9706640.441950958, 9710110.938596912, 9711336.244502613, 9754081.25, 9765270.150989942, 9768143.763976626, 9801535.397468517, 9804400.914825995, 9814025.153936988, 9817267.84673101, 9817740.324264565, 9817784.53537451, 9818118.75, 9818284.439746203, 9818298.129739815, 9818518.75401193, 9819407.417905157, 9819928.172236305, 9820739.0625, 9822303.445409037, 9822758.170460366, 9823079.964889828, 9823281.339533838, 9823735.809455218, 9823939.853611143, 9824493.070789374, 9824848.306342283, 9824912.095960679, 9825208.922525693, 9826026.147432694, 9826245.085697332, 9826267.001211539, 9826732.050847774, 9827350.0, 9827810.761088971, 9831209.961995056, 9840444.318977376, 9842151.5625, 9842727.491664417, 9854598.074789742, 9855262.771473255, 9856192.09117025, 9856197.579067843, 9872913.511049548, 9886643.50121541, 9890120.844563497, 9891354.6875, 9892342.61468205, 9898347.202098096, 9900263.739561766, 9907278.125, 9939453.994210582, 9949859.689833472, 9951289.923148287, 9984108.625222363, 9984126.5625, 10007128.125, 10007372.137489181, 10009089.614021001, 10013912.33083273, 10014116.727380464, 10015307.68547173, 10018446.105300602, 10019328.68269172, 10024282.734712167, 10026312.695555907, 10026445.702386715, 10027767.861249648, 10028839.0625, 10029917.1875, 10035609.375, 10036142.132504603, 10036410.9375, 10037533.624144655, 10061839.0625, 10063397.979330145, 10072938.61562099, 10073888.988652343, 10076427.512934284, 10076628.125, 10084051.19274822, 10116061.035375463, 10126977.886553975, 10149259.375, 10157478.169020737, 10163592.765153438, 10164237.883640427, 10164788.829404725, 10164839.416025998, 10170884.375, 10170892.13330737, 10173927.544992948, 10196498.98989321, 10203669.42634784, 10203674.379509365, 10204014.135250777, 10205352.778397124, 10207341.37849955, 10209378.812626692, 10211544.916665554, 10213093.75, 10213858.372417396, 10213872.430699797, 10214065.625, 10218065.474249193, 10237555.724615576, 10238848.987491501, 10240943.75, 10296221.41665913, 10296309.907359548, 10298082.8125, 10298438.268401809, 10299260.9375, 10300863.32957696, 10301972.91339047, 10302721.875, 10302995.495858198, 10303339.421607787, 10303477.78248735, 10303543.75, 10303545.470577497, 10303741.73410811, 10303925.0, 10305243.511840178, 10305698.052348431, 10306057.959726714, 10314070.616445133, 10324254.14542219, 10372448.194981279, 10380000.242500465, 10386225.140979253, 10387783.551169975, 10388081.25, 10389148.379518785, 10389757.86504857, 10390478.125, 10390523.4375, 10390762.788446251, 10392996.875, 10396074.99820794, 10396432.63452368, 10396526.63375702, 10396935.14668103, 10397268.75, 10397306.840768661, 10397718.299602026, 10397768.73538385, 10398174.170509558, 10398295.3125, 10399236.613046618, 10399598.82372481, 10400141.430366283, 10400161.725282285, 10400202.685622443, 10404360.728222607, 10432733.368654316, 10435409.375, 10435624.09256947, 10435714.418971682, 10449332.565323243, 10454811.76591008, 10516926.722230682, 10522560.343609666, 10522574.863961294, 10522577.445805307, 10524158.330515428, 10528937.5, 10530455.449537195, 10532096.079064209, 10534935.9375, 10538225.303373078, 10542774.469408182, 10542916.484300751, 10542943.475709692, 10544071.541011686, 10544528.155959854, 10545671.397630995, 10545790.06601156, 10545815.513190813, 10546058.815893244, 10548848.593949974, 10564817.987001283, 10568182.168388236, 10582660.9375, 10588945.3125, 10589417.934475303, 10589520.639396966, 10589757.8125, 10590440.405347642, 10591008.014278261, 10591416.879572209, 10591537.5, 10591761.368099602, 10591817.835589144, 10591835.9375, 10592053.125, 10592286.277655097, 10592458.73411632, 10592472.488433879, 10592759.375, 10593278.559160301, 10594103.125, 10594131.594740573, 10594145.168887274, 10594284.25667788, 10594623.157352129, 10595312.08726856, 10595661.06267492, 10596050.0, 10596196.875, 10596215.993059194, 10596243.22360108, 10596926.5625, 10597052.59136572, 10598710.9227283, 10599052.84061991, 10599985.416991087, 10609299.083491072, 10636650.811206527, 10636653.125, 10637817.382119207, 10638064.0625, 10640930.697161864, 10642357.277146436, 10643373.1084644, 10644029.060808377, 10645199.428642595, 10647466.784150237, 10668005.966715604, 10670204.6875, 10670398.4375, 10672816.55273429, 10673974.64164251, 10674109.26085461, 10674489.0625, 10688870.3125, 10697141.90810875, 10697733.29197915, 10698082.8125, 10698412.5, 10698554.351246353, 10699064.3394213, 10699954.191797096, 10700435.9375, 10721221.875, 10725882.082000937, 10727024.737788638, 10729781.798060134, 10730305.034707038, 10731214.0625, 10731894.777215632, 10733339.901713032, 10738946.07258351, 10739159.283389304, 10740117.1875, 10740879.45022794, 10741406.398486733, 10746292.164069183, 10749454.483458579, 10752243.75, 10770221.166200245, 10772271.862388028, 10772978.125, 10772987.5, 10773065.625, 10773746.63209447, 10774124.357733347, 10774435.66191616, 10774984.375, 10775233.165574025, 10775287.996270292, 10775420.3125, 10775515.625, 10775524.93598372, 10775842.000110881, 10775871.875, 10776688.587699244, 10776810.036570204, 10776953.125, 10777241.288671251, 10777275.84324734, 10777395.898087712, 10777742.1875, 10777834.379076317, 10778701.440969745, 10779051.442084434, 10779072.02657176, 10779553.125, 10779578.125, 10779583.61533177, 10779954.626848776, 10779995.097503902, 10780004.41989751, 10780092.1875, 10780421.799909264, 10780706.25, 10780872.910663037, 10780884.327216959, 10780944.361719858, 10781143.75, 10781166.578915285, 10781292.1875, 10781373.070245069, 10781506.983192934, 10781625.0, 10781878.045459421, 10782285.235903585, 10782326.115524577, 10782425.41484338, 10782654.218193593, 10783016.704770388, 10783053.045161724, 10783215.145005118, 10783225.0, 10783269.317205181, 10783338.072656099, 10783341.93095099, 10783640.625, 10783951.5625, 10783977.47172084, 10783984.78956332, 10785651.5625, 10786607.8125, 10786811.991781265, 10787028.034716027, 10787391.489974169, 10788310.77878921, 10788553.125, 10788565.471452907, 10798530.938055158, 10800907.817632362, 10802299.403113931, 10803359.375, 10804462.5, 10808588.674578888, 10808612.316061666, 10819733.614552058, 10822183.7059848, 10832126.537001822, 10836026.627184968, 10837050.790442552, 10841107.163319295, 10844948.135790141, 10847010.586759742, 10847312.545251215, 10864289.010278149, 10864658.09354318, 10865036.024376439, 10865141.129382402, 10866945.380682684, 10868675.593782062, 10876874.661304872, 10922466.858693598, 10953160.196094112, 10953997.715523146, 10958459.507378101, 10959701.860012108, 10960353.715128489, 10960354.6875, 10961534.50683899, 10962338.203074178, 10962462.656191004, 10962866.249768902, 10963769.149754304, 10964746.421447467, 10966654.146663897, 10967927.723956537, 10968253.125, 10970155.528168725, 10971875.0, 10973522.891605634, 10977303.84571479, 10981945.3125, 10988909.871059148, 10992332.473321991, 10995832.39009761, 10996101.117891068, 10996200.0, 11016298.367706403, 11017429.611613475, 11017741.530830815, 11018125.015003186, 11019336.846298076, 11021028.275282966, 11023297.95419125, 11032236.952595845, 11035385.9375, 11063679.574795235, 11065327.91903799, 11080902.638586815, 11082199.964970492, 11083810.99171903, 11085159.375, 11086743.629448319, 11087113.35734961, 11089652.251600103, 11090215.925553828, 11090673.805488931, 11158989.232536165, 11160354.018419227, 11167437.5, 11168110.359819282, 11169225.0, 11177554.6875, 11181131.141750481, 11187539.338001879, 11189665.625, 11191737.236550096, 11206631.56392146, 11211884.647520201, 11214126.620531777, 11295855.643720938, 11295979.6875, 11296564.923850661, 11296744.506728716, 11297624.610689642, 11313015.414485484, 11314030.550396627, 11314727.885854278, 11340253.396617088, 11395433.91478568, 11403448.114480177, 11492980.7912614, 11493967.76805489, 11647892.585424425, 16228155.364985706, 16966407.405190874, 16966407.735301036, 16966409.26151319, 16966411.043223858, 16966419.23521041, 16966424.215175677, 16966426.399346806, 16966426.474611294, 17019036.096217223, 17026679.6875, 17028646.630960144, 17028732.8125, 17029206.70469579, 17029578.92315375, 17029942.104518156, 17030474.429124456, 17031115.54155804, 17033077.043456983, 17033225.618404295, 17033288.33227493, 17034658.879324015, 17034967.1875, 17035715.3996179, 17035835.29785716, 17035968.75, 17038375.0, 17039714.693723567, 17040379.539999146, 17043613.4145505, 17043743.75, 17045100.0, 17048337.5], [10.065459696394544, 12.326574375071807, 8.159669690351624, 67.90080682193609, 12.673214011517171, 78.83340792819155, 22.15804401485982, 28.235848172362726, 24.376230438396476, 6.024202519815038, 41.934350152355144, 6.6125498325183, 27.336288191385698, 7.552640313244892, 57.18113545269022, 6.189297189141456, 31.063538671754323, 17.95679062663734, 15.339751575277347, 33.718717044999835, 76.87883709580312, 58.841636439642684, 12.215641638910325, 21.037222353227897, 61.34892485615109, 94.39492082299645, 34.968065838312704, 19.519983544344502, 53.713227745639436, 13.859576650916187, 69.1274456075273, 17.01002492809936, 57.045100904849875, 19.543815337145745, 29.255540507501614, 54.02718570924705, 23.285083349427197, 18.555261184776526, 164.49361023456058, 96.27009658228715, 96.51793925946248, 87.73290846818036, 42.94959745352294, 50.0104046091421, 9.454500010236732, 6.178756896108188, 5.302864789419282, 60.99320441044666, 123.04106404617539, 5.394826747899507, 56.02721212423233, 34.221435534876804, 13.159289941202506, 36.18866569165773, 20.47218165547233, 70.4145057115702, 23.511384517176914, 20.119371644724723, 49.29620062267678, 115.97819889059197, 76.28279020347254, 71.18256387459354, 20.906135036125214, 26.05840568469945, 78.00236695436391, 5.1772996169510295, 10.756801979318887, 117.23148808335588, 20.3920821489323, 27.85495734193822, 18.48687371385993, 33.46518069902873, 13.747328481721174, 57.1806240353857, 138.5696484394292, 141.78809681086875, 38.62134949131011, 88.10219637428074, 9.705975125328253, 20.2363178140899, 5.815905221489664, 122.35462834937749, 60.6242781811777, 39.31723922460627, 17.324307946503897, 36.66445185487366, 16.858542597302872, 5.197819362108079, 11.642109655699372, 5.033871245317018, 11.759753899138556, 160.5072730740715, 10.522827436614566, 6.529419460729219, 59.70910964443793, 75.55180408811208, 11.499017128810143, 14.619315347952877, 24.06474157228263, 5.027469132053725, 7.247404026985484, 69.4093628626539, 119.9425404591059, 208.23390962816254, 147.80559523951146, 275.8375934686521, 59.76773777879205, 20.433967735746013, 9.24916410158467, 91.07225708151792, 14.535181632232376, 34.506380310638164, 15.47946640011782, 37.77053765525138, 40.039300725286665, 7.1653715294197005, 64.64485516763472, 157.8528911155728, 99.88735124733068, 26.31051016322639, 232.62944169638956, 10.475661248583219, 194.03580927282368, 63.47737823014859, 64.0177193969828, 6.85836347983801, 48.723145962181015, 27.116895715423034, 113.48950163865842, 30.865936465142692, 23.990290725278545, 122.33336364120123, 106.70438197321933, 25.556181183450228, 83.28954850096201, 108.73852395216578, 8.269566281733962, 29.717047787788154, 14.585252213038318, 38.30902162389954, 28.123613319398967, 52.873992016953096, 46.40061248258097, 10.661139281559956, 19.2063892367583, 64.67726967439827, 89.10508649879478, 101.67991738707937, 92.73435444361691, 69.43333169784071, 84.0124461749939, 97.86917353449624, 26.674556670267737, 83.1773630397125, 8.001728409237012, 31.864799635065236, 99.00394800279389, 16.968792435542497, 47.83386924890513, 24.238586174245164, 19.916205608514428, 118.04414115898064, 17.698062327434872, 15.937448260462798, 188.2653959157293, 7.177113463957441, 41.58661621733429, 6.848276235933935, 5.734491269605495, 35.510178982449574, 56.23593188966842, 205.3470599882475, 71.46467667378529, 152.57855942250472, 27.710727269360188, 33.43903012921928, 8.93597250843404, 42.720581694753115, 57.31326136866095, 22.731773109874915, 17.95266389864998, 16.301969352994643, 67.51546792021543, 12.263950379629946, 12.496519965979616, 81.46236944322362, 11.29522535927037, 29.213805194518436, 9.628022280951724, 223.57620509587545, 18.274314696791553, 30.03376882115046, 81.87519684522694, 12.973413642526028, 32.15820699348023, 130.45946612629496, 33.942662110412485, 145.635519367409, 7.809169911195373, 69.77107079458565, 27.902277642309922, 142.2131935205544, 79.71398832422027, 53.8838063839399, 36.62525566514588, 32.328151625181405, 6.066160129017725, 5.197894504185542, 9.875113520031144, 12.762697001482076, 31.3254844012218, 315.80796899854676, 5.533877222599556, 19.538112312733475, 86.63581119026658, 23.82456708047306, 23.01881321963276, 24.259030202183915, 11.591118065352681, 6.595205822988215, 66.64062821332942, 6.407640974682039, 62.00215721879985, 66.05933862480246, 40.426794174372844, 32.61087589486444, 137.3036201013023, 8.078970779357652, 44.05283632857133, 69.42436241361507, 139.57198843343193, 32.20323066491779, 479.60358455066694, 49.27135503473501, 9.874451174467843, 121.73531173395843, 69.58737526305512, 5.670386154241651, 137.93884700119983, 6.989574262155711, 105.48470569001742, 52.41244420290899, 6.698849362435621, 7.355789738177282, 13.039597944474542, 80.22688052604067, 28.910954195080766, 17.795463977158093, 16.935344048199802, 47.471275247959596, 110.49523320644775, 75.89619257546168, 5.264717220175217, 70.27631937158931, 14.227596071218645, 43.38274343953543, 65.2639489197569, 15.38720946124976, 10.379315583679716, 34.57712925055541, 27.892364691119386, 98.08434940176893, 158.15084275439068, 150.40985011948044, 9.92824594925841, 13.24961199017178, 125.10850388964877, 9.253644913489001, 61.74350832904565, 21.312926313502608, 168.06902980265556, 39.46650572302364, 29.242406965757322, 25.59309003571559, 91.99012681649751, 42.626109500489285, 13.226877837918321, 55.56789156222786, 99.15619445953797, 60.110690088661094, 23.26775324595096, 18.95882675822094, 17.839463591663954, 20.150888427470214, 17.978768316925354, 16.430335044549018, 5.992982945100992, 63.76914484638368, 85.97126911898198, 20.6908408946346, 37.06738855120284, 89.10996031127964, 28.782809981411603, 110.3937149153673, 39.61050126644795, 279.91825738110845, 101.6387698961279, 8.395657668498057, 5.769360846555785, 25.400363662995677, 23.23301823929342, 11.634916322668293, 48.356174667413796, 165.50736595347345, 38.58038692391099, 118.14936549096926, 68.21676231015972, 63.24843990277795, 10.297236993788147, 7.6604414773303, 16.300830201272582, 20.549848923035157, 20.00026263552206, 59.12997577296272, 11.237816453166808, 6.626457273464917, 18.0372646034879, 45.47213889630851, 338.66284978512766, 86.72734122303476, 16.0370270911384, 25.821890802518194, 73.13304793216244, 22.355939496397795, 6.165096365903978, 15.336425672218802, 23.190540028702983, 6.173524488040721, 6.127554851299917, 28.564689700970696, 90.70776920376426, 53.378220714067396, 13.018559156914483, 11.922167191350615, 37.06442693598332, 27.983291031244303, 128.23797468473612, 17.749849216711628, 159.7344748869917, 64.07675349860621, 11.804764732056647, 79.57103909975811, 19.893197669174945, 145.63421374644633, 9.60259526611847, 110.47286501984465, 14.117279831033876, 78.1129848380767, 106.9369372691047, 14.934163914400377, 72.66787981357452, 66.46051691694319, 127.10221219197, 73.44920303039733, 36.37537185311975, 41.780700034748875, 66.14160531688321, 121.67854605858096, 39.640005971909794, 72.44008334460138, 11.656337448957867, 67.27050074006516, 85.51835470166353, 28.25052885260542, 37.73617390642234, 16.390418778497622, 14.291284854072707, 121.2407542868341, 174.3094470280064, 5.477410906521663, 73.72028038674958, 62.22240867066296, 33.095014466396385, 11.36912875067891, 9.377324724586217, 7.2073987790956755, 5.9478073414078505, 87.03264983686122, 179.72046613126642, 27.141477073248826, 69.83435110743261, 26.66025927609796, 18.455824200760375, 18.029648430908413, 31.628242956647103, 42.58030262632378, 31.476992650902265, 8.938655988111936, 15.627165870220361, 6.003651321862843, 19.59925403533, 25.457045490189284, 134.9695500164637, 36.329553137335836, 24.994473797092187, 42.71499633458385, 8.184477586045038, 59.60599125811213, 58.54091170179636, 84.93361830342135, 24.634654765351552, 26.78600116319455, 24.364745470245747, 155.23175199740942, 80.98047339367035, 34.432593171301136, 15.393998790065904, 27.447998714878107, 6.668309169395795, 7.340203695421902, 69.35328828634582, 29.17519662993057, 159.0314774688759, 77.184407289088, 44.92183480620913, 93.92492811851584, 50.38908036339462, 14.01998398779195, 61.21445045635694, 95.34837126174128, 61.93788664363602, 10.208108917672467, 13.745266441722075, 136.30732278817152, 8.635766017497307, 46.41499155280363, 46.55713767002439, 51.63752346378877, 31.051016118110184, 56.081423998232935, 18.739290541337922, 16.202950061627664, 14.217057326011915, 38.58277585856864, 76.46731269646023, 108.82405792402594, 26.15026590113674, 20.51541097116636, 5.956007547241795, 30.916344647362077, 149.4839133046695, 30.502535883437528, 12.088132908065273, 27.600288351746094, 57.08060228599658, 6.931248623199689, 8.43300841421391, 105.8130201421919, 116.19300313285602, 73.3297633313933, 15.821300020237416, 8.369167356200602, 15.158714989085617, 39.48771082669069, 20.1309020098496, 16.39901706671206, 18.238767426921566, 20.99583613804402, 102.11129709086629, 77.65320227576704, 29.06666785773202, 74.12603212756545, 69.41010720163604, 52.90979270018388, 34.81196920966274, 47.12058290840053, 50.67416763950388, 11.042537421217501, 103.12659872798778, 12.163743139299722, 71.60845435651234, 12.237057728643467, 14.954845654235575, 142.2898355441057, 39.21564365064953, 8.764197323283433, 53.57715230122585, 14.367513716039007, 8.636014621737234, 19.00258351596674, 214.7182086708191, 33.9593156473441, 64.03086381255758, 128.00294858978313, 10.341310038781456, 35.60166793607071, 62.889541442550914, 7.4389627404447465, 100.67797950943924, 5.9403492647604885, 27.936928293769004, 93.60228502830421, 56.60210916048675, 52.440015388456644, 25.36916149703492, 6.267336964182166, 142.86993435731887, 77.48843123356143, 25.697109967647346, 89.42970172773758, 201.86661988401644, 52.417022386164064, 5.622617448234552, 30.615128794190746, 51.54192424453679, 32.03288908526014, 41.73177361975442, 15.194171078163473, 48.56539518339554, 61.89643205659379, 121.60747512517852, 5.292548587274148, 47.52382290839143, 16.32169601344293, 19.109159741396557, 9.577508153934511, 6.5012186363983036, 9.225431071071085, 19.154820962656537, 70.57575143040913, 26.176998967952287, 32.18661949449039, 17.151401885736995, 36.39151453257347, 5.197681925558865, 69.89341217299983, 40.84599490468364, 246.90893509464712, 89.09501314793873, 28.72664120378201, 131.95434371899222, 91.21018374741894, 6.64332676985908, 10.604448881328048, 55.26382310678627, 14.58663186725983, 61.20125539035429, 97.20093871987211, 24.464620053201173, 8.589519297009351, 76.7025913549202, 84.63074616163107, 7.725988783089798, 32.044679196995645, 104.69502620111811, 6.183407453262966, 301.9032276027648, 68.23134702754746, 19.31740992465969, 85.22836757300185, 39.63781040147598, 78.05571346886943, 5.361182281697263, 70.42191858798284, 10.477379553265456, 26.24395950377896, 7.917570078999962, 283.29435942754867, 8.245865770958975, 116.73378282189418, 16.66974794586078, 16.28322935215401, 97.02527452853384, 93.2037563919223, 44.12266641989035, 7.5834671629781525, 79.47055430307769, 72.70418532662393, 88.76014787255907, 15.469562506365564, 93.0984403950143, 27.144503354023733, 66.3639181817211, 45.391404108765286, 106.26864580433606, 15.396465654520192, 58.746573210664195, 72.30458406472137, 9.368904645548543, 17.720125862958486, 8.50043622373016, 24.644709058910614, 5.729632911506516, 72.183788932704, 32.56472921881566, 20.730464268369563, 52.89002266113788, 178.0563791885188, 47.14942698220819, 29.621451336517723, 185.26046426300917, 18.880115192690663, 12.404593359390443, 22.599198691572123, 72.6465322524311, 383.32416430876106, 99.09549650019648, 95.86189969155802, 69.70951041318037, 43.17127263607583, 11.624529203698046, 104.55858651647182, 52.89531522158302, 12.076010612986154, 83.28068194113732, 10.078280218566665, 45.40544323654357, 19.44671083833787, 6.663196572492929, 63.24347467103084, 9.491409282598319, 138.49925508672797, 7.483882707337109, 38.384597204178505, 15.58032004392454, 38.53712379414376, 62.94225979247324, 33.62344592642402, 17.71585488700972, 48.043733934347266, 73.04478714127971, 13.212263432792328, 24.51384567024685, 15.95405697925217, 90.35402324414531, 15.7985270107426, 70.28664657480216, 48.26671124703293, 9.376234752126175, 52.724639014308906, 56.00612547973738, 16.168089650258885, 41.089293710678675, 7.6025896541126095, 86.44607878882387, 15.220246285385699, 50.23777935004686, 32.31150764119185, 17.800217993589715, 58.98211207271988, 28.2388441793577, 34.49001661001151, 55.849097805765346, 18.341996323764768, 17.79121592841217, 12.672977159194982, 6.174040860273673, 69.53184609556901, 24.983161194525028, 56.25727018970835, 100.9915571042396, 139.08525438361283, 5.0754364483683325, 141.6603488874823, 96.13179946636197, 15.587407655527041, 5.669975746392192, 14.447163606509626, 8.8163261594362, 150.98989757725772, 53.61217849057084, 115.73764120992095, 18.974374273025973, 29.537627376195985, 16.239331841585514, 13.243035902456594, 93.90975460366869, 13.989858136777631, 23.326084413973316, 11.176446342871463, 362.39255354598663, 10.353568338503006, 35.148138151397305, 87.66175691216887, 8.371080529492417, 10.113239952736825, 6.96183792526149, 135.0409280951333, 146.46341861609227, 125.41377121128748, 84.93206909539391, 59.28490254738445, 5.920502442664166, 114.07319500294477, 13.884781643939055, 37.164875501287774, 9.724125692667407, 23.428802564744636, 31.334570732072123, 53.9847987377896, 11.704551526353407, 83.99686123483983, 17.940913422183446, 41.28985228945841, 6.8533933084314445, 81.20282521392664, 27.286396885378352, 25.97749289303358, 26.3864870672791, 24.007681170223993, 65.35252547563826, 40.42881319245795, 52.45053160245735, 76.07737497511714, 5.972112315864063, 60.594771615975816, 50.15551984015993, 7.420998016902543, 28.2358327752387, 13.958097751026655, 13.270225723141753, 5.749142163068124, 16.965340944839255, 58.98720347255606, 33.15442047556916, 125.52647676316052, 47.98119672721745, 32.49420622901363, 35.62810355048529, 6.644525704252957, 8.890893967153294, 133.73462616344085, 37.27390272299792, 67.94663672529002, 99.82435922776703, 30.909821652988278, 59.95325817888908, 26.455524223267073, 59.75990833382976, 88.82381543047235, 42.35341226897886, 171.47926504204983, 110.1623147324641, 66.08631837690804, 23.298440263118618, 98.03271895407283, 21.536073437454448, 41.13782277900193, 31.35816424224633, 30.78355607556629, 88.62757950465914, 17.542859718756944, 28.183194195654025, 36.158604397184504, 39.19304850130103, 23.696178202837846, 301.12955990407653, 122.82182654189265, 67.70194188011621, 24.58649720745498, 5.622311939638599, 6.264242584517828, 25.059082512659224, 113.62676198197086, 185.30195965643145, 117.47350733752184, 36.40414025937308, 62.10297631053583, 32.97330868378769, 47.65373627319608, 86.61934269209016, 42.72068977812477, 50.109052070575046, 16.482853237803578, 145.58825427444586, 40.33743802688613, 82.02927139611825, 27.31915552584696, 44.00442719645317, 47.50431005778305, 21.748213995091653, 9.961345667156532, 51.047474997846194, 46.29669297096342, 24.73272789478049, 17.7684646626583, 29.043442417272896, 9.33485231090715, 19.62129069750685, 10.944565181586293, 11.37318681516897, 51.95682191914705, 55.581556185500595, 52.22554472494025, 7.581447488853529, 24.90087336217376, 22.60715365755828, 65.6785166311473, 39.7517630450412, 32.55035642590353, 100.4412471799809, 27.613464713728156, 17.93017086130963, 45.08593786324791, 7.397017692495815, 61.417531761955246, 54.891105686862744, 40.45820534823959, 54.73987180374691, 14.646763949633657, 27.627511104933678, 104.75154875608781, 119.04327473969308, 5.276807539656052, 13.373595747806007, 5.06744735110097, 13.381171313498967, 204.94114190167687, 20.272921093612005, 68.55689119311319, 14.046215243124289, 99.85135948458216, 45.555620558410055, 8.8340225818876, 6.554102487354428, 40.62616220040944, 31.4115098150593, 118.0483098395159, 26.85837008469606, 25.22798335550291, 5.161681720697702, 40.47379749839956, 14.386867990922868, 30.97194167445604, 48.854923975740114, 44.35617078161087, 36.06529695493804, 54.651992854470954, 10.571059908188062, 61.12192188125745, 16.7178819683983, 185.83646826845256, 49.31376504909399, 190.1516808460131, 81.67649128028869, 109.88853086326847, 19.667635845824883, 99.23130428068568, 16.148365717740667, 80.82221162082917, 12.789701790520379, 5.100169953670094, 104.94380869176253, 59.44527635606772, 72.54733814399464, 117.3186700377615, 30.88603038998292, 174.0046456231751, 12.75661265184694, 74.77876702923287, 62.92098599606597, 17.79628633064022, 57.86146682378606, 112.41108283832766, 6.60325717437112, 9.467058199607353, 65.30352443305922, 123.5737953727717, 21.936016525223863, 34.37039282683617, 46.15017625309061, 56.629963079487304, 134.1043870097884, 47.73172656438225, 21.047178710776805, 12.478254968234028, 39.162956909214635, 17.242486315825293, 23.26287267313329, 29.174765264560158, 148.24773162620716, 84.310106893874, 14.681998524766712, 38.51802091207599, 6.639125960566017, 33.27205361564582, 57.056555491201664, 6.103868004718433, 10.447885066075768, 48.79289887500329, 113.74595070618177, 6.7412126123920775, 25.03116931522865, 249.31740185213908, 27.8500276297756, 15.152274044739151, 127.32745324234878, 154.50425460083176, 19.664709440802774, 8.688018078395125, 9.599161313075442, 55.073017372851055, 9.479439374646871, 6.186210584107579, 101.40786959246094, 20.751680855190685, 45.99453428874781, 95.46101296402868, 14.544364039424249, 16.106996110291963, 54.28799929660118, 24.229122620054678, 8.459866719167694, 9.572978789467573, 55.37952011887972, 57.742413330089306, 73.50700039661027, 83.16175028438187, 27.1766141658076, 119.10774957498592, 24.96530416807874, 41.11463218909563, 7.059975158322338, 176.53832848186468, 78.05949088325752, 35.11174603679402, 21.914822164344994, 178.25477306200423, 13.466039325466113, 8.680577160422564, 197.06124588064873, 131.70261997563264, 139.8993133124308, 82.90671994441763, 15.50703962848715, 5.223698068289605, 5.240601799538259, 5.201188546917007, 9.469312847282577, 37.02017232787812, 10.820932184672952, 69.97585055411137, 56.501639957575364, 15.337162307062552, 8.243177265278701, 15.317181192462161, 5.807831964779528, 11.777096421987164, 5.649197026536696, 19.467615982095005, 21.53634330521458, 79.3231112237834, 24.753091962443133, 6.189630701765054, 89.14052766737274, 103.74618285607986, 12.282053091461476, 107.53580062706234, 17.723495552468894, 47.96542030644419, 50.29343427466593, 114.90962770801391])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3910014.8194018146, 4365871.041709561, 4366648.189549966, 4432115.625, 4456751.305011168, 4586187.5, 4675888.037006122, 4724910.922729095, 4725035.211758849, 4725908.737295396, 4764932.8125, 4897414.818716402, 4898139.329151454, 4899593.607904949, 4962067.1875, 4993750.919991304, 4995000.0, 4996474.053462202, 5000587.600701952, 5016280.019841695, 5089785.784367688, 5111960.9375, 5117663.099922496, 5138038.515552952, 5147295.3125, 5173257.713381661, 5189553.125, 5272208.651500754, 7858878.118125062, 7998937.183893107, 7999271.875, 8008642.21425294, 8140826.5625, 8232093.029546252, 8318254.6875, 8329914.0625, 8330007.732865723, 8331063.2071464, 8333396.7204927895, 8334824.156808111, 8335331.25, 8343903.125, 8351388.481924267, 8352501.5625, 8352670.454579472, 8358210.203884542, 8426044.126635674, 8461484.221274178, 8478923.603848565, 8519017.09176515, 8519671.163631419, 8521453.50796835, 8521743.187961867, 8521950.0, 8522193.141182508, 8522765.625, 8523133.438325534, 8523599.651927931, 8524160.620010953, 8524376.351478698, 8524385.9375, 8524450.38444078, 8524485.376575157, 8526485.359116502, 8526504.6875, 8527680.192819232, 8528430.99899828, 8528771.33624739, 8528825.44824642, 8535415.59023353, 8537355.384617068, 8538428.125, 8540335.584436208, 8542365.486523516, 8543185.396108765, 8545460.67703779, 8562026.5625, 8566509.375, 8569592.599586722, 8580070.66442481, 8625198.46348207, 8682684.375, 8688032.8125, 8689419.988341711, 8702893.325506413, 8711207.59329776, 8714309.800957955, 8715674.778663142, 8716381.555926824, 8716432.352189507, 8716797.244935399, 8717168.47455672, 8717377.700179731, 8717397.139809752, 8719137.388298994, 8719156.25, 8719828.350238103, 8719842.100319996, 8720666.233044147, 8721428.835282575, 8723305.145912651, 8731212.5, 8743274.07983356, 8745881.16577824, 8746935.789313221, 8748029.549046723, 8748348.4375, 8748386.886841381, 8748434.27611158, 8752620.3125, 8752626.451647926, 8780590.625, 8785692.321825508, 8811787.5, 8833841.627961911, 8834132.976984683, 8834940.625, 8843214.0625, 8843798.4375, 8843978.639688369, 8845145.792431833, 8850379.226035608, 8903631.25, 8906150.0, 8907231.25, 8907286.394763708, 8922509.375, 8923197.898707053, 8932870.3125, 8934614.0625, 8935296.593524855, 8935996.206659578, 8936816.290710649, 8936818.94084929, 8938843.452153, 8939674.78728012, 8939681.217692744, 8940223.4375, 8940226.465782732, 8940475.0, 8940801.541903665, 8940842.162362223, 8940964.569668746, 8941298.116985053, 8941391.42301674, 8941741.547489034, 8942217.1875, 8942303.125, 8942547.782576567, 8942575.720424745, 8942590.106948173, 8942618.75, 8943098.336435124, 8943220.392540742, 8944040.044485684, 8945493.75, 8953315.625, 8954069.196523713, 8954075.931590796, 8954633.512970926, 8955016.118982531, 8955032.206794204, 8957075.752796663, 8957185.515313458, 8957507.414243396, 8958693.957597852, 8960132.8125, 8976781.660278516, 8976805.114621326, 8978153.915198157, 8978576.529812444, 8979989.0625, 8979995.3125, 8980849.119752444, 8981021.250518287, 8981613.940275062, 8983196.367137792, 8983217.471941829, 8988609.371579766, 8990137.967347275, 8995624.655706385, 8995888.105106452, 8996627.04051476, 9001272.497743383, 9007156.978826828, 9008153.125, 9008340.729593217, 9008618.75, 9008909.706133831, 9009090.625, 9009458.477815984, 9009523.614459397, 9009653.125, 9009974.365987148, 9012239.819104508, 9012655.364676045, 9013360.340155289, 9015226.5625, 9015297.410609175, 9015531.25, 9015546.84275286, 9016079.17015203, 9016456.25, 9016796.875, 9016939.0625, 9018065.47048876, 9019007.598853597, 9019424.781650782, 9020216.204898566, 9020700.219281388, 9021750.0, 9022009.375, 9022039.07890427, 9022117.054477932, 9022137.428465886, 9022166.893257631, 9022683.500176037, 9023096.11375978, 9023111.215038115, 9023174.094177118, 9023442.564468732, 9023861.97372736, 9024143.75, 9025093.75, 9025234.138274176, 9025281.25, 9025393.75, 9025399.585635705, 9025634.375, 9025823.4375, 9025824.3051855, 9025962.5, 9026859.829160407, 9027274.890437229, 9027672.454261405, 9027911.069156326, 9028061.852266202, 9028183.787916828, 9028326.5625, 9028390.098379996, 9028460.9375, 9029149.50941549, 9029163.729770934, 9029282.072274571, 9029302.516967375, 9029634.375, 9029941.853473887, 9030044.33999454, 9030186.37782236, 9030281.25, 9030500.0, 9030503.084601222, 9030520.997055272, 9031378.303729942, 9032110.656052893, 9032149.973909745, 9033028.125, 9033934.487218443, 9034633.83690556, 9034652.957940388, 9034899.967225933, 9035402.147926139, 9036098.50580103, 9036178.125, 9038883.443018883, 9039474.471852457, 9040592.1875, 9041704.594084738, 9042593.75, 9046849.481746728, 9051306.141803468, 9052698.150654785, 9052714.877641186, 9052753.570555959, 9056847.325566074, 9062510.9375, 9063555.71338701, 9064218.75, 9065071.132197255, 9065638.770624513, 9066702.78622392, 9067486.701920865, 9067615.324827421, 9068552.479030035, 9069988.11945227, 9069993.731729507, 9071403.720326673, 9077095.218276614, 9077495.3125, 9079097.853132205, 9079171.875, 9079740.625, 9083656.25, 9086332.8125, 9086332.8125, 9086787.12871061, 9086810.9375, 9089397.075729573, 9089400.940688163, 9090137.319029123, 9090449.658358864, 9091783.114446089, 9092626.502024684, 9093303.125, 9093757.052094826, 9097567.1875, 9100103.428038187, 9101284.339614797, 9102750.957168976, 9105148.592938282, 9105280.947507199, 9107571.274942506, 9110897.773587795, 9116896.875, 9177286.006820817, 9178843.02244361, 9179099.066186946, 9179265.625, 9179705.40038077, 9184450.984484904, 9187837.109803181, 9198187.346512605, 9199600.0, 9208436.538005788, 9215647.58424489, 9215656.349891137, 9220460.592340775, 9233321.095769065, 9243669.121055583, 9255972.01304002, 9256781.524706552, 9256958.579869783, 9261357.196251675, 9265636.392539795, 9265651.5625, 9266696.8470426, 9268285.614791855, 9271367.805880785, 9283673.4375, 9298368.75, 9298392.857950756, 9313939.0625, 9318899.676206494, 9319383.217727654, 9320184.679209884, 9322546.172988059, 9326092.102253003, 9327159.375, 9327232.8125, 9331745.418387065, 9332896.875, 9335329.6875, 9342043.011454243, 9342114.0625, 9342346.758064944, 9342512.5, 9342928.125, 9344628.054153837, 9347010.9375, 9347018.26855879, 9347328.438777333, 9348664.0625, 9348967.261186296, 9349613.654678162, 9351163.772366062, 9353159.709445283, 9353695.130763272, 9354211.458132006, 9355114.0625, 9359686.29699372, 9370001.5625, 9375216.905901369, 9376018.75, 9378689.21165212, 9385837.553814525, 9386630.004669448, 9390265.720974179, 9391996.875, 9392456.093642805, 9396526.980268562, 9396584.330165122, 9397627.740121828, 9399044.44253145, 9400444.377489066, 9402468.75, 9404839.0625, 9405517.1875, 9405838.55837837, 9411281.41896581, 9411321.021800742, 9414594.16525168, 9416834.197788605, 9417010.9375, 9417518.75, 9417782.604005229, 9418923.233139755, 9419249.801805234, 9420269.086537804, 9421682.8125, 9422029.560480203, 9422380.041298838, 9422544.85604374, 9422851.309094522, 9428275.925191343, 9439181.094588943, 9446658.828845475, 9454511.05135791, 9454546.430704959, 9455671.378693778, 9456116.39656712, 9456397.468078075, 9456695.3125, 9456853.125, 9457146.77241501, 9459726.213059133, 9460774.386727327, 9461456.25, 9462055.985219043, 9462303.125, 9463882.8125, 9465464.021324731, 9492736.819528794, 9492779.333749047, 9508375.550402498, 9509081.856690725, 9513296.875, 9517017.754068272, 9518942.87787002, 9526969.779341122, 9534720.3125, 9555951.343560228, 9562494.045772677, 9567825.112379488, 9568079.286456928, 9568466.94948145, 9570962.5, 9570974.599218486, 9571044.70856867, 9572650.188685182, 9573213.878311653, 9573415.625, 9576210.83247661, 9583905.22582162, 9586254.638998972, 9587317.1875, 9587349.22339871, 9587553.392005583, 9589885.9375, 9590143.75, 9593087.5, 9603923.650641385, 9603982.965480449, 9604390.680842701, 9604546.875, 9606937.238704486, 9609308.850961689, 9609516.284902824, 9611308.823663715, 9613093.75, 9624877.759948839, 9628488.58698483, 9630603.673003042, 9630767.874217683, 9632660.9375, 9632876.838467829, 9652425.38970571, 9665323.025330754, 9679040.569819259, 9680973.498929648, 9681725.575091952, 9684459.375, 9684514.672573822, 9684690.628494522, 9685429.077217486, 9685441.636011573, 9686747.68076135, 9689137.5, 9689807.00624082, 9689820.919587554, 9691253.029709218, 9692564.080210293, 9692689.781769462, 9695335.059501812, 9695995.892970309, 9696205.681000935, 9696382.8125, 9696424.074754816, 9698119.834255315, 9698979.514743427, 9699112.510830484, 9699218.226893734, 9701579.83995475, 9702131.217521526, 9705942.1875, 9706640.441950958, 9710110.938596912, 9711336.244502613, 9754081.25, 9765270.150989942, 9768143.763976626, 9801535.397468517, 9804400.914825995, 9814025.153936988, 9817267.84673101, 9817740.324264565, 9817784.53537451, 9818118.75, 9818284.439746203, 9818298.129739815, 9818518.75401193, 9819407.417905157, 9819928.172236305, 9820739.0625, 9822303.445409037, 9822758.170460366, 9823079.964889828, 9823281.339533838, 9823735.809455218, 9823939.853611143, 9824493.070789374, 9824848.306342283, 9824912.095960679, 9825208.922525693, 9826026.147432694, 9826245.085697332, 9826267.001211539, 9826732.050847774, 9827350.0, 9827810.761088971, 9831209.961995056, 9840444.318977376, 9842151.5625, 9842727.491664417, 9854598.074789742, 9855262.771473255, 9856192.09117025, 9856197.579067843, 9872913.511049548, 9886643.50121541, 9890120.844563497, 9891354.6875, 9892342.61468205, 9898347.202098096, 9900263.739561766, 9907278.125, 9939453.994210582, 9949859.689833472, 9951289.923148287, 9984108.625222363, 9984126.5625, 10007128.125, 10007372.137489181, 10009089.614021001, 10013912.33083273, 10014116.727380464, 10015307.68547173, 10018446.105300602, 10019328.68269172, 10024282.734712167, 10026312.695555907, 10026445.702386715, 10027767.861249648, 10028839.0625, 10029917.1875, 10035609.375, 10036142.132504603, 10036410.9375, 10037533.624144655, 10061839.0625, 10063397.979330145, 10072938.61562099, 10073888.988652343, 10076427.512934284, 10076628.125, 10084051.19274822, 10116061.035375463, 10126977.886553975, 10149259.375, 10157478.169020737, 10163592.765153438, 10164237.883640427, 10164788.829404725, 10164839.416025998, 10170884.375, 10170892.13330737, 10173927.544992948, 10196498.98989321, 10203669.42634784, 10203674.379509365, 10204014.135250777, 10205352.778397124, 10207341.37849955, 10209378.812626692, 10211544.916665554, 10213093.75, 10213858.372417396, 10213872.430699797, 10214065.625, 10218065.474249193, 10237555.724615576, 10238848.987491501, 10240943.75, 10296221.41665913, 10296309.907359548, 10298082.8125, 10298438.268401809, 10299260.9375, 10300863.32957696, 10301972.91339047, 10302721.875, 10302995.495858198, 10303339.421607787, 10303477.78248735, 10303543.75, 10303545.470577497, 10303741.73410811, 10303925.0, 10305243.511840178, 10305698.052348431, 10306057.959726714, 10314070.616445133, 10324254.14542219, 10372448.194981279, 10380000.242500465, 10386225.140979253, 10387783.551169975, 10388081.25, 10389148.379518785, 10389757.86504857, 10390478.125, 10390523.4375, 10390762.788446251, 10392996.875, 10396074.99820794, 10396432.63452368, 10396526.63375702, 10396935.14668103, 10397268.75, 10397306.840768661, 10397718.299602026, 10397768.73538385, 10398174.170509558, 10398295.3125, 10399236.613046618, 10399598.82372481, 10400141.430366283, 10400161.725282285, 10400202.685622443, 10404360.728222607, 10432733.368654316, 10435409.375, 10435624.09256947, 10435714.418971682, 10449332.565323243, 10454811.76591008, 10516926.722230682, 10522560.343609666, 10522574.863961294, 10522577.445805307, 10524158.330515428, 10528937.5, 10530455.449537195, 10532096.079064209, 10534935.9375, 10538225.303373078, 10542774.469408182, 10542916.484300751, 10542943.475709692, 10544071.541011686, 10544528.155959854, 10545671.397630995, 10545790.06601156, 10545815.513190813, 10546058.815893244, 10548848.593949974, 10564817.987001283, 10568182.168388236, 10582660.9375, 10588945.3125, 10589417.934475303, 10589520.639396966, 10589757.8125, 10590440.405347642, 10591008.014278261, 10591416.879572209, 10591537.5, 10591761.368099602, 10591817.835589144, 10591835.9375, 10592053.125, 10592286.277655097, 10592458.73411632, 10592472.488433879, 10592759.375, 10593278.559160301, 10594103.125, 10594131.594740573, 10594145.168887274, 10594284.25667788, 10594623.157352129, 10595312.08726856, 10595661.06267492, 10596050.0, 10596196.875, 10596215.993059194, 10596243.22360108, 10596926.5625, 10597052.59136572, 10598710.9227283, 10599052.84061991, 10599985.416991087, 10609299.083491072, 10636650.811206527, 10636653.125, 10637817.382119207, 10638064.0625, 10640930.697161864, 10642357.277146436, 10643373.1084644, 10644029.060808377, 10645199.428642595, 10647466.784150237, 10668005.966715604, 10670204.6875, 10670398.4375, 10672816.55273429, 10673974.64164251, 10674109.26085461, 10674489.0625, 10688870.3125, 10697141.90810875, 10697733.29197915, 10698082.8125, 10698412.5, 10698554.351246353, 10699064.3394213, 10699954.191797096, 10700435.9375, 10721221.875, 10725882.082000937, 10727024.737788638, 10729781.798060134, 10730305.034707038, 10731214.0625, 10731894.777215632, 10733339.901713032, 10738946.07258351, 10739159.283389304, 10740117.1875, 10740879.45022794, 10741406.398486733, 10746292.164069183, 10749454.483458579, 10752243.75, 10770221.166200245, 10772271.862388028, 10772978.125, 10772987.5, 10773065.625, 10773746.63209447, 10774124.357733347, 10774435.66191616, 10774984.375, 10775233.165574025, 10775287.996270292, 10775420.3125, 10775515.625, 10775524.93598372, 10775842.000110881, 10775871.875, 10776688.587699244, 10776810.036570204, 10776953.125, 10777241.288671251, 10777275.84324734, 10777395.898087712, 10777742.1875, 10777834.379076317, 10778701.440969745, 10779051.442084434, 10779072.02657176, 10779553.125, 10779578.125, 10779583.61533177, 10779954.626848776, 10779995.097503902, 10780004.41989751, 10780092.1875, 10780421.799909264, 10780706.25, 10780872.910663037, 10780884.327216959, 10780944.361719858, 10781143.75, 10781166.578915285, 10781292.1875, 10781373.070245069, 10781506.983192934, 10781625.0, 10781878.045459421, 10782285.235903585, 10782326.115524577, 10782425.41484338, 10782654.218193593, 10783016.704770388, 10783053.045161724, 10783215.145005118, 10783225.0, 10783269.317205181, 10783338.072656099, 10783341.93095099, 10783640.625, 10783951.5625, 10783977.47172084, 10783984.78956332, 10785651.5625, 10786607.8125, 10786811.991781265, 10787028.034716027, 10787391.489974169, 10788310.77878921, 10788553.125, 10788565.471452907, 10798530.938055158, 10800907.817632362, 10802299.403113931, 10803359.375, 10804462.5, 10808588.674578888, 10808612.316061666, 10819733.614552058, 10822183.7059848, 10832126.537001822, 10836026.627184968, 10837050.790442552, 10841107.163319295, 10844948.135790141, 10847010.586759742, 10847312.545251215, 10864289.010278149, 10864658.09354318, 10865036.024376439, 10865141.129382402, 10866945.380682684, 10868675.593782062, 10876874.661304872, 10922466.858693598, 10953160.196094112, 10953997.715523146, 10958459.507378101, 10959701.860012108, 10960353.715128489, 10960354.6875, 10961534.50683899, 10962338.203074178, 10962462.656191004, 10962866.249768902, 10963769.149754304, 10964746.421447467, 10966654.146663897, 10967927.723956537, 10968253.125, 10970155.528168725, 10971875.0, 10973522.891605634, 10977303.84571479, 10981945.3125, 10988909.871059148, 10992332.473321991, 10995832.39009761, 10996101.117891068, 10996200.0, 11016298.367706403, 11017429.611613475, 11017741.530830815, 11018125.015003186, 11019336.846298076, 11021028.275282966, 11023297.95419125, 11032236.952595845, 11035385.9375, 11063679.574795235, 11065327.91903799, 11080902.638586815, 11082199.964970492, 11083810.99171903, 11085159.375, 11086743.629448319, 11087113.35734961, 11089652.251600103, 11090215.925553828, 11090673.805488931, 11158989.232536165, 11160354.018419227, 11167437.5, 11168110.359819282, 11169225.0, 11177554.6875, 11181131.141750481, 11187539.338001879, 11189665.625, 11191737.236550096, 11206631.56392146, 11211884.647520201, 11214126.620531777, 11295855.643720938, 11295979.6875, 11296564.923850661, 11296744.506728716, 11297624.610689642, 11313015.414485484, 11314030.550396627, 11314727.885854278, 11340253.396617088, 11395433.91478568, 11403448.114480177, 11492980.7912614, 11493967.76805489, 11647892.585424425, 16228155.364985706, 16966407.405190874, 16966407.735301036, 16966409.26151319, 16966411.043223858, 16966419.23521041, 16966424.215175677, 16966426.399346806, 16966426.474611294, 17019036.096217223, 17026679.6875, 17028646.630960144, 17028732.8125, 17029206.70469579, 17029578.92315375, 17029942.104518156, 17030474.429124456, 17031115.54155804, 17033077.043456983, 17033225.618404295, 17033288.33227493, 17034658.879324015, 17034967.1875, 17035715.3996179, 17035835.29785716, 17035968.75, 17038375.0, 17039714.693723567, 17040379.539999146, 17043613.4145505, 17043743.75, 17045100.0, 17048337.5], [10.065459696394544, 12.326574375071807, 8.159669690351624, 67.90080682193609, 12.673214011517171, 78.83340792819155, 22.15804401485982, 28.235848172362726, 24.376230438396476, 6.024202519815038, 41.934350152355144, 6.6125498325183, 27.336288191385698, 7.552640313244892, 57.18113545269022, 6.189297189141456, 31.063538671754323, 17.95679062663734, 15.339751575277347, 33.718717044999835, 76.87883709580312, 58.841636439642684, 12.215641638910325, 21.037222353227897, 61.34892485615109, 94.39492082299645, 34.968065838312704, 19.519983544344502, 53.713227745639436, 13.859576650916187, 69.1274456075273, 17.01002492809936, 57.045100904849875, 19.543815337145745, 29.255540507501614, 54.02718570924705, 23.285083349427197, 18.555261184776526, 164.49361023456058, 96.27009658228715, 96.51793925946248, 87.73290846818036, 42.94959745352294, 50.0104046091421, 9.454500010236732, 6.178756896108188, 5.302864789419282, 60.99320441044666, 123.04106404617539, 5.394826747899507, 56.02721212423233, 34.221435534876804, 13.159289941202506, 36.18866569165773, 20.47218165547233, 70.4145057115702, 23.511384517176914, 20.119371644724723, 49.29620062267678, 115.97819889059197, 76.28279020347254, 71.18256387459354, 20.906135036125214, 26.05840568469945, 78.00236695436391, 5.1772996169510295, 10.756801979318887, 117.23148808335588, 20.3920821489323, 27.85495734193822, 18.48687371385993, 33.46518069902873, 13.747328481721174, 57.1806240353857, 138.5696484394292, 141.78809681086875, 38.62134949131011, 88.10219637428074, 9.705975125328253, 20.2363178140899, 5.815905221489664, 122.35462834937749, 60.6242781811777, 39.31723922460627, 17.324307946503897, 36.66445185487366, 16.858542597302872, 5.197819362108079, 11.642109655699372, 5.033871245317018, 11.759753899138556, 160.5072730740715, 10.522827436614566, 6.529419460729219, 59.70910964443793, 75.55180408811208, 11.499017128810143, 14.619315347952877, 24.06474157228263, 5.027469132053725, 7.247404026985484, 69.4093628626539, 119.9425404591059, 208.23390962816254, 147.80559523951146, 275.8375934686521, 59.76773777879205, 20.433967735746013, 9.24916410158467, 91.07225708151792, 14.535181632232376, 34.506380310638164, 15.47946640011782, 37.77053765525138, 40.039300725286665, 7.1653715294197005, 64.64485516763472, 157.8528911155728, 99.88735124733068, 26.31051016322639, 232.62944169638956, 10.475661248583219, 194.03580927282368, 63.47737823014859, 64.0177193969828, 6.85836347983801, 48.723145962181015, 27.116895715423034, 113.48950163865842, 30.865936465142692, 23.990290725278545, 122.33336364120123, 106.70438197321933, 25.556181183450228, 83.28954850096201, 108.73852395216578, 8.269566281733962, 29.717047787788154, 14.585252213038318, 38.30902162389954, 28.123613319398967, 52.873992016953096, 46.40061248258097, 10.661139281559956, 19.2063892367583, 64.67726967439827, 89.10508649879478, 101.67991738707937, 92.73435444361691, 69.43333169784071, 84.0124461749939, 97.86917353449624, 26.674556670267737, 83.1773630397125, 8.001728409237012, 31.864799635065236, 99.00394800279389, 16.968792435542497, 47.83386924890513, 24.238586174245164, 19.916205608514428, 118.04414115898064, 17.698062327434872, 15.937448260462798, 188.2653959157293, 7.177113463957441, 41.58661621733429, 6.848276235933935, 5.734491269605495, 35.510178982449574, 56.23593188966842, 205.3470599882475, 71.46467667378529, 152.57855942250472, 27.710727269360188, 33.43903012921928, 8.93597250843404, 42.720581694753115, 57.31326136866095, 22.731773109874915, 17.95266389864998, 16.301969352994643, 67.51546792021543, 12.263950379629946, 12.496519965979616, 81.46236944322362, 11.29522535927037, 29.213805194518436, 9.628022280951724, 223.57620509587545, 18.274314696791553, 30.03376882115046, 81.87519684522694, 12.973413642526028, 32.15820699348023, 130.45946612629496, 33.942662110412485, 145.635519367409, 7.809169911195373, 69.77107079458565, 27.902277642309922, 142.2131935205544, 79.71398832422027, 53.8838063839399, 36.62525566514588, 32.328151625181405, 6.066160129017725, 5.197894504185542, 9.875113520031144, 12.762697001482076, 31.3254844012218, 315.80796899854676, 5.533877222599556, 19.538112312733475, 86.63581119026658, 23.82456708047306, 23.01881321963276, 24.259030202183915, 11.591118065352681, 6.595205822988215, 66.64062821332942, 6.407640974682039, 62.00215721879985, 66.05933862480246, 40.426794174372844, 32.61087589486444, 137.3036201013023, 8.078970779357652, 44.05283632857133, 69.42436241361507, 139.57198843343193, 32.20323066491779, 479.60358455066694, 49.27135503473501, 9.874451174467843, 121.73531173395843, 69.58737526305512, 5.670386154241651, 137.93884700119983, 6.989574262155711, 105.48470569001742, 52.41244420290899, 6.698849362435621, 7.355789738177282, 13.039597944474542, 80.22688052604067, 28.910954195080766, 17.795463977158093, 16.935344048199802, 47.471275247959596, 110.49523320644775, 75.89619257546168, 5.264717220175217, 70.27631937158931, 14.227596071218645, 43.38274343953543, 65.2639489197569, 15.38720946124976, 10.379315583679716, 34.57712925055541, 27.892364691119386, 98.08434940176893, 158.15084275439068, 150.40985011948044, 9.92824594925841, 13.24961199017178, 125.10850388964877, 9.253644913489001, 61.74350832904565, 21.312926313502608, 168.06902980265556, 39.46650572302364, 29.242406965757322, 25.59309003571559, 91.99012681649751, 42.626109500489285, 13.226877837918321, 55.56789156222786, 99.15619445953797, 60.110690088661094, 23.26775324595096, 18.95882675822094, 17.839463591663954, 20.150888427470214, 17.978768316925354, 16.430335044549018, 5.992982945100992, 63.76914484638368, 85.97126911898198, 20.6908408946346, 37.06738855120284, 89.10996031127964, 28.782809981411603, 110.3937149153673, 39.61050126644795, 279.91825738110845, 101.6387698961279, 8.395657668498057, 5.769360846555785, 25.400363662995677, 23.23301823929342, 11.634916322668293, 48.356174667413796, 165.50736595347345, 38.58038692391099, 118.14936549096926, 68.21676231015972, 63.24843990277795, 10.297236993788147, 7.6604414773303, 16.300830201272582, 20.549848923035157, 20.00026263552206, 59.12997577296272, 11.237816453166808, 6.626457273464917, 18.0372646034879, 45.47213889630851, 338.66284978512766, 86.72734122303476, 16.0370270911384, 25.821890802518194, 73.13304793216244, 22.355939496397795, 6.165096365903978, 15.336425672218802, 23.190540028702983, 6.173524488040721, 6.127554851299917, 28.564689700970696, 90.70776920376426, 53.378220714067396, 13.018559156914483, 11.922167191350615, 37.06442693598332, 27.983291031244303, 128.23797468473612, 17.749849216711628, 159.7344748869917, 64.07675349860621, 11.804764732056647, 79.57103909975811, 19.893197669174945, 145.63421374644633, 9.60259526611847, 110.47286501984465, 14.117279831033876, 78.1129848380767, 106.9369372691047, 14.934163914400377, 72.66787981357452, 66.46051691694319, 127.10221219197, 73.44920303039733, 36.37537185311975, 41.780700034748875, 66.14160531688321, 121.67854605858096, 39.640005971909794, 72.44008334460138, 11.656337448957867, 67.27050074006516, 85.51835470166353, 28.25052885260542, 37.73617390642234, 16.390418778497622, 14.291284854072707, 121.2407542868341, 174.3094470280064, 5.477410906521663, 73.72028038674958, 62.22240867066296, 33.095014466396385, 11.36912875067891, 9.377324724586217, 7.2073987790956755, 5.9478073414078505, 87.03264983686122, 179.72046613126642, 27.141477073248826, 69.83435110743261, 26.66025927609796, 18.455824200760375, 18.029648430908413, 31.628242956647103, 42.58030262632378, 31.476992650902265, 8.938655988111936, 15.627165870220361, 6.003651321862843, 19.59925403533, 25.457045490189284, 134.9695500164637, 36.329553137335836, 24.994473797092187, 42.71499633458385, 8.184477586045038, 59.60599125811213, 58.54091170179636, 84.93361830342135, 24.634654765351552, 26.78600116319455, 24.364745470245747, 155.23175199740942, 80.98047339367035, 34.432593171301136, 15.393998790065904, 27.447998714878107, 6.668309169395795, 7.340203695421902, 69.35328828634582, 29.17519662993057, 159.0314774688759, 77.184407289088, 44.92183480620913, 93.92492811851584, 50.38908036339462, 14.01998398779195, 61.21445045635694, 95.34837126174128, 61.93788664363602, 10.208108917672467, 13.745266441722075, 136.30732278817152, 8.635766017497307, 46.41499155280363, 46.55713767002439, 51.63752346378877, 31.051016118110184, 56.081423998232935, 18.739290541337922, 16.202950061627664, 14.217057326011915, 38.58277585856864, 76.46731269646023, 108.82405792402594, 26.15026590113674, 20.51541097116636, 5.956007547241795, 30.916344647362077, 149.4839133046695, 30.502535883437528, 12.088132908065273, 27.600288351746094, 57.08060228599658, 6.931248623199689, 8.43300841421391, 105.8130201421919, 116.19300313285602, 73.3297633313933, 15.821300020237416, 8.369167356200602, 15.158714989085617, 39.48771082669069, 20.1309020098496, 16.39901706671206, 18.238767426921566, 20.99583613804402, 102.11129709086629, 77.65320227576704, 29.06666785773202, 74.12603212756545, 69.41010720163604, 52.90979270018388, 34.81196920966274, 47.12058290840053, 50.67416763950388, 11.042537421217501, 103.12659872798778, 12.163743139299722, 71.60845435651234, 12.237057728643467, 14.954845654235575, 142.2898355441057, 39.21564365064953, 8.764197323283433, 53.57715230122585, 14.367513716039007, 8.636014621737234, 19.00258351596674, 214.7182086708191, 33.9593156473441, 64.03086381255758, 128.00294858978313, 10.341310038781456, 35.60166793607071, 62.889541442550914, 7.4389627404447465, 100.67797950943924, 5.9403492647604885, 27.936928293769004, 93.60228502830421, 56.60210916048675, 52.440015388456644, 25.36916149703492, 6.267336964182166, 142.86993435731887, 77.48843123356143, 25.697109967647346, 89.42970172773758, 201.86661988401644, 52.417022386164064, 5.622617448234552, 30.615128794190746, 51.54192424453679, 32.03288908526014, 41.73177361975442, 15.194171078163473, 48.56539518339554, 61.89643205659379, 121.60747512517852, 5.292548587274148, 47.52382290839143, 16.32169601344293, 19.109159741396557, 9.577508153934511, 6.5012186363983036, 9.225431071071085, 19.154820962656537, 70.57575143040913, 26.176998967952287, 32.18661949449039, 17.151401885736995, 36.39151453257347, 5.197681925558865, 69.89341217299983, 40.84599490468364, 246.90893509464712, 89.09501314793873, 28.72664120378201, 131.95434371899222, 91.21018374741894, 6.64332676985908, 10.604448881328048, 55.26382310678627, 14.58663186725983, 61.20125539035429, 97.20093871987211, 24.464620053201173, 8.589519297009351, 76.7025913549202, 84.63074616163107, 7.725988783089798, 32.044679196995645, 104.69502620111811, 6.183407453262966, 301.9032276027648, 68.23134702754746, 19.31740992465969, 85.22836757300185, 39.63781040147598, 78.05571346886943, 5.361182281697263, 70.42191858798284, 10.477379553265456, 26.24395950377896, 7.917570078999962, 283.29435942754867, 8.245865770958975, 116.73378282189418, 16.66974794586078, 16.28322935215401, 97.02527452853384, 93.2037563919223, 44.12266641989035, 7.5834671629781525, 79.47055430307769, 72.70418532662393, 88.76014787255907, 15.469562506365564, 93.0984403950143, 27.144503354023733, 66.3639181817211, 45.391404108765286, 106.26864580433606, 15.396465654520192, 58.746573210664195, 72.30458406472137, 9.368904645548543, 17.720125862958486, 8.50043622373016, 24.644709058910614, 5.729632911506516, 72.183788932704, 32.56472921881566, 20.730464268369563, 52.89002266113788, 178.0563791885188, 47.14942698220819, 29.621451336517723, 185.26046426300917, 18.880115192690663, 12.404593359390443, 22.599198691572123, 72.6465322524311, 383.32416430876106, 99.09549650019648, 95.86189969155802, 69.70951041318037, 43.17127263607583, 11.624529203698046, 104.55858651647182, 52.89531522158302, 12.076010612986154, 83.28068194113732, 10.078280218566665, 45.40544323654357, 19.44671083833787, 6.663196572492929, 63.24347467103084, 9.491409282598319, 138.49925508672797, 7.483882707337109, 38.384597204178505, 15.58032004392454, 38.53712379414376, 62.94225979247324, 33.62344592642402, 17.71585488700972, 48.043733934347266, 73.04478714127971, 13.212263432792328, 24.51384567024685, 15.95405697925217, 90.35402324414531, 15.7985270107426, 70.28664657480216, 48.26671124703293, 9.376234752126175, 52.724639014308906, 56.00612547973738, 16.168089650258885, 41.089293710678675, 7.6025896541126095, 86.44607878882387, 15.220246285385699, 50.23777935004686, 32.31150764119185, 17.800217993589715, 58.98211207271988, 28.2388441793577, 34.49001661001151, 55.849097805765346, 18.341996323764768, 17.79121592841217, 12.672977159194982, 6.174040860273673, 69.53184609556901, 24.983161194525028, 56.25727018970835, 100.9915571042396, 139.08525438361283, 5.0754364483683325, 141.6603488874823, 96.13179946636197, 15.587407655527041, 5.669975746392192, 14.447163606509626, 8.8163261594362, 150.98989757725772, 53.61217849057084, 115.73764120992095, 18.974374273025973, 29.537627376195985, 16.239331841585514, 13.243035902456594, 93.90975460366869, 13.989858136777631, 23.326084413973316, 11.176446342871463, 362.39255354598663, 10.353568338503006, 35.148138151397305, 87.66175691216887, 8.371080529492417, 10.113239952736825, 6.96183792526149, 135.0409280951333, 146.46341861609227, 125.41377121128748, 84.93206909539391, 59.28490254738445, 5.920502442664166, 114.07319500294477, 13.884781643939055, 37.164875501287774, 9.724125692667407, 23.428802564744636, 31.334570732072123, 53.9847987377896, 11.704551526353407, 83.99686123483983, 17.940913422183446, 41.28985228945841, 6.8533933084314445, 81.20282521392664, 27.286396885378352, 25.97749289303358, 26.3864870672791, 24.007681170223993, 65.35252547563826, 40.42881319245795, 52.45053160245735, 76.07737497511714, 5.972112315864063, 60.594771615975816, 50.15551984015993, 7.420998016902543, 28.2358327752387, 13.958097751026655, 13.270225723141753, 5.749142163068124, 16.965340944839255, 58.98720347255606, 33.15442047556916, 125.52647676316052, 47.98119672721745, 32.49420622901363, 35.62810355048529, 6.644525704252957, 8.890893967153294, 133.73462616344085, 37.27390272299792, 67.94663672529002, 99.82435922776703, 30.909821652988278, 59.95325817888908, 26.455524223267073, 59.75990833382976, 88.82381543047235, 42.35341226897886, 171.47926504204983, 110.1623147324641, 66.08631837690804, 23.298440263118618, 98.03271895407283, 21.536073437454448, 41.13782277900193, 31.35816424224633, 30.78355607556629, 88.62757950465914, 17.542859718756944, 28.183194195654025, 36.158604397184504, 39.19304850130103, 23.696178202837846, 301.12955990407653, 122.82182654189265, 67.70194188011621, 24.58649720745498, 5.622311939638599, 6.264242584517828, 25.059082512659224, 113.62676198197086, 185.30195965643145, 117.47350733752184, 36.40414025937308, 62.10297631053583, 32.97330868378769, 47.65373627319608, 86.61934269209016, 42.72068977812477, 50.109052070575046, 16.482853237803578, 145.58825427444586, 40.33743802688613, 82.02927139611825, 27.31915552584696, 44.00442719645317, 47.50431005778305, 21.748213995091653, 9.961345667156532, 51.047474997846194, 46.29669297096342, 24.73272789478049, 17.7684646626583, 29.043442417272896, 9.33485231090715, 19.62129069750685, 10.944565181586293, 11.37318681516897, 51.95682191914705, 55.581556185500595, 52.22554472494025, 7.581447488853529, 24.90087336217376, 22.60715365755828, 65.6785166311473, 39.7517630450412, 32.55035642590353, 100.4412471799809, 27.613464713728156, 17.93017086130963, 45.08593786324791, 7.397017692495815, 61.417531761955246, 54.891105686862744, 40.45820534823959, 54.73987180374691, 14.646763949633657, 27.627511104933678, 104.75154875608781, 119.04327473969308, 5.276807539656052, 13.373595747806007, 5.06744735110097, 13.381171313498967, 204.94114190167687, 20.272921093612005, 68.55689119311319, 14.046215243124289, 99.85135948458216, 45.555620558410055, 8.8340225818876, 6.554102487354428, 40.62616220040944, 31.4115098150593, 118.0483098395159, 26.85837008469606, 25.22798335550291, 5.161681720697702, 40.47379749839956, 14.386867990922868, 30.97194167445604, 48.854923975740114, 44.35617078161087, 36.06529695493804, 54.651992854470954, 10.571059908188062, 61.12192188125745, 16.7178819683983, 185.83646826845256, 49.31376504909399, 190.1516808460131, 81.67649128028869, 109.88853086326847, 19.667635845824883, 99.23130428068568, 16.148365717740667, 80.82221162082917, 12.789701790520379, 5.100169953670094, 104.94380869176253, 59.44527635606772, 72.54733814399464, 117.3186700377615, 30.88603038998292, 174.0046456231751, 12.75661265184694, 74.77876702923287, 62.92098599606597, 17.79628633064022, 57.86146682378606, 112.41108283832766, 6.60325717437112, 9.467058199607353, 65.30352443305922, 123.5737953727717, 21.936016525223863, 34.37039282683617, 46.15017625309061, 56.629963079487304, 134.1043870097884, 47.73172656438225, 21.047178710776805, 12.478254968234028, 39.162956909214635, 17.242486315825293, 23.26287267313329, 29.174765264560158, 148.24773162620716, 84.310106893874, 14.681998524766712, 38.51802091207599, 6.639125960566017, 33.27205361564582, 57.056555491201664, 6.103868004718433, 10.447885066075768, 48.79289887500329, 113.74595070618177, 6.7412126123920775, 25.03116931522865, 249.31740185213908, 27.8500276297756, 15.152274044739151, 127.32745324234878, 154.50425460083176, 19.664709440802774, 8.688018078395125, 9.599161313075442, 55.073017372851055, 9.479439374646871, 6.186210584107579, 101.40786959246094, 20.751680855190685, 45.99453428874781, 95.46101296402868, 14.544364039424249, 16.106996110291963, 54.28799929660118, 24.229122620054678, 8.459866719167694, 9.572978789467573, 55.37952011887972, 57.742413330089306, 73.50700039661027, 83.16175028438187, 27.1766141658076, 119.10774957498592, 24.96530416807874, 41.11463218909563, 7.059975158322338, 176.53832848186468, 78.05949088325752, 35.11174603679402, 21.914822164344994, 178.25477306200423, 13.466039325466113, 8.680577160422564, 197.06124588064873, 131.70261997563264, 139.8993133124308, 82.90671994441763, 15.50703962848715, 5.223698068289605, 5.240601799538259, 5.201188546917007, 9.469312847282577, 37.02017232787812, 10.820932184672952, 69.97585055411137, 56.501639957575364, 15.337162307062552, 8.243177265278701, 15.317181192462161, 5.807831964779528, 11.777096421987164, 5.649197026536696, 19.467615982095005, 21.53634330521458, 79.3231112237834, 24.753091962443133, 6.189630701765054, 89.14052766737274, 103.74618285607986, 12.282053091461476, 107.53580062706234, 17.723495552468894, 47.96542030644419, 50.29343427466593, 114.90962770801391])
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);
([3910014.8194018146, 4365871.041709561, 4366648.189549966, 4432115.625, 4456751.305011168, 4586187.5, 4675888.037006122, 4724910.922729095, 4725035.211758849, 4725908.737295396, 4764932.8125, 4897414.818716402, 4898139.329151454, 4899593.607904949, 4962067.1875, 4993750.919991304, 4995000.0, 4996474.053462202, 5000587.600701952, 5016280.019841695, 5089785.784367688, 5111960.9375, 5117663.099922496, 5138038.515552952, 5147295.3125, 5173257.713381661, 5189553.125, 5272208.651500754, 7858878.118125062, 7998937.183893107, 7999271.875, 8008642.21425294, 8140826.5625, 8232093.029546252, 8318254.6875, 8329914.0625, 8330007.732865723, 8331063.2071464, 8333396.7204927895, 8334824.156808111, 8335331.25, 8343903.125, 8351388.481924267, 8352501.5625, 8352670.454579472, 8358210.203884542, 8426044.126635674, 8461484.221274178, 8478923.603848565, 8519017.09176515, 8519671.163631419, 8521453.50796835, 8521743.187961867, 8521950.0, 8522193.141182508, 8522765.625, 8523133.438325534, 8523599.651927931, 8524160.620010953, 8524376.351478698, 8524385.9375, 8524450.38444078, 8524485.376575157, 8526485.359116502, 8526504.6875, 8527680.192819232, 8528430.99899828, 8528771.33624739, 8528825.44824642, 8535415.59023353, 8537355.384617068, 8538428.125, 8540335.584436208, 8542365.486523516, 8543185.396108765, 8545460.67703779, 8562026.5625, 8566509.375, 8569592.599586722, 8580070.66442481, 8625198.46348207, 8682684.375, 8688032.8125, 8689419.988341711, 8702893.325506413, 8711207.59329776, 8714309.800957955, 8715674.778663142, 8716381.555926824, 8716432.352189507, 8716797.244935399, 8717168.47455672, 8717377.700179731, 8717397.139809752, 8719137.388298994, 8719156.25, 8719828.350238103, 8719842.100319996, 8720666.233044147, 8721428.835282575, 8723305.145912651, 8731212.5, 8743274.07983356, 8745881.16577824, 8746935.789313221, 8748029.549046723, 8748348.4375, 8748386.886841381, 8748434.27611158, 8752620.3125, 8752626.451647926, 8780590.625, 8785692.321825508, 8811787.5, 8833841.627961911, 8834132.976984683, 8834940.625, 8843214.0625, 8843798.4375, 8843978.639688369, 8845145.792431833, 8850379.226035608, 8903631.25, 8906150.0, 8907231.25, 8907286.394763708, 8922509.375, 8923197.898707053, 8932870.3125, 8934614.0625, 8935296.593524855, 8935996.206659578, 8936816.290710649, 8936818.94084929, 8938843.452153, 8939674.78728012, 8939681.217692744, 8940223.4375, 8940226.465782732, 8940475.0, 8940801.541903665, 8940842.162362223, 8940964.569668746, 8941298.116985053, 8941391.42301674, 8941741.547489034, 8942217.1875, 8942303.125, 8942547.782576567, 8942575.720424745, 8942590.106948173, 8942618.75, 8943098.336435124, 8943220.392540742, 8944040.044485684, 8945493.75, 8953315.625, 8954069.196523713, 8954075.931590796, 8954633.512970926, 8955016.118982531, 8955032.206794204, 8957075.752796663, 8957185.515313458, 8957507.414243396, 8958693.957597852, 8960132.8125, 8976781.660278516, 8976805.114621326, 8978153.915198157, 8978576.529812444, 8979989.0625, 8979995.3125, 8980849.119752444, 8981021.250518287, 8981613.940275062, 8983196.367137792, 8983217.471941829, 8988609.371579766, 8990137.967347275, 8995624.655706385, 8995888.105106452, 8996627.04051476, 9001272.497743383, 9007156.978826828, 9008153.125, 9008340.729593217, 9008618.75, 9008909.706133831, 9009090.625, 9009458.477815984, 9009523.614459397, 9009653.125, 9009974.365987148, 9012239.819104508, 9012655.364676045, 9013360.340155289, 9015226.5625, 9015297.410609175, 9015531.25, 9015546.84275286, 9016079.17015203, 9016456.25, 9016796.875, 9016939.0625, 9018065.47048876, 9019007.598853597, 9019424.781650782, 9020216.204898566, 9020700.219281388, 9021750.0, 9022009.375, 9022039.07890427, 9022117.054477932, 9022137.428465886, 9022166.893257631, 9022683.500176037, 9023096.11375978, 9023111.215038115, 9023174.094177118, 9023442.564468732, 9023861.97372736, 9024143.75, 9025093.75, 9025234.138274176, 9025281.25, 9025393.75, 9025399.585635705, 9025634.375, 9025823.4375, 9025824.3051855, 9025962.5, 9026859.829160407, 9027274.890437229, 9027672.454261405, 9027911.069156326, 9028061.852266202, 9028183.787916828, 9028326.5625, 9028390.098379996, 9028460.9375, 9029149.50941549, 9029163.729770934, 9029282.072274571, 9029302.516967375, 9029634.375, 9029941.853473887, 9030044.33999454, 9030186.37782236, 9030281.25, 9030500.0, 9030503.084601222, 9030520.997055272, 9031378.303729942, 9032110.656052893, 9032149.973909745, 9033028.125, 9033934.487218443, 9034633.83690556, 9034652.957940388, 9034899.967225933, 9035402.147926139, 9036098.50580103, 9036178.125, 9038883.443018883, 9039474.471852457, 9040592.1875, 9041704.594084738, 9042593.75, 9046849.481746728, 9051306.141803468, 9052698.150654785, 9052714.877641186, 9052753.570555959, 9056847.325566074, 9062510.9375, 9063555.71338701, 9064218.75, 9065071.132197255, 9065638.770624513, 9066702.78622392, 9067486.701920865, 9067615.324827421, 9068552.479030035, 9069988.11945227, 9069993.731729507, 9071403.720326673, 9077095.218276614, 9077495.3125, 9079097.853132205, 9079171.875, 9079740.625, 9083656.25, 9086332.8125, 9086332.8125, 9086787.12871061, 9086810.9375, 9089397.075729573, 9089400.940688163, 9090137.319029123, 9090449.658358864, 9091783.114446089, 9092626.502024684, 9093303.125, 9093757.052094826, 9097567.1875, 9100103.428038187, 9101284.339614797, 9102750.957168976, 9105148.592938282, 9105280.947507199, 9107571.274942506, 9110897.773587795, 9116896.875, 9177286.006820817, 9178843.02244361, 9179099.066186946, 9179265.625, 9179705.40038077, 9184450.984484904, 9187837.109803181, 9198187.346512605, 9199600.0, 9208436.538005788, 9215647.58424489, 9215656.349891137, 9220460.592340775, 9233321.095769065, 9243669.121055583, 9255972.01304002, 9256781.524706552, 9256958.579869783, 9261357.196251675, 9265636.392539795, 9265651.5625, 9266696.8470426, 9268285.614791855, 9271367.805880785, 9283673.4375, 9298368.75, 9298392.857950756, 9313939.0625, 9318899.676206494, 9319383.217727654, 9320184.679209884, 9322546.172988059, 9326092.102253003, 9327159.375, 9327232.8125, 9331745.418387065, 9332896.875, 9335329.6875, 9342043.011454243, 9342114.0625, 9342346.758064944, 9342512.5, 9342928.125, 9344628.054153837, 9347010.9375, 9347018.26855879, 9347328.438777333, 9348664.0625, 9348967.261186296, 9349613.654678162, 9351163.772366062, 9353159.709445283, 9353695.130763272, 9354211.458132006, 9355114.0625, 9359686.29699372, 9370001.5625, 9375216.905901369, 9376018.75, 9378689.21165212, 9385837.553814525, 9386630.004669448, 9390265.720974179, 9391996.875, 9392456.093642805, 9396526.980268562, 9396584.330165122, 9397627.740121828, 9399044.44253145, 9400444.377489066, 9402468.75, 9404839.0625, 9405517.1875, 9405838.55837837, 9411281.41896581, 9411321.021800742, 9414594.16525168, 9416834.197788605, 9417010.9375, 9417518.75, 9417782.604005229, 9418923.233139755, 9419249.801805234, 9420269.086537804, 9421682.8125, 9422029.560480203, 9422380.041298838, 9422544.85604374, 9422851.309094522, 9428275.925191343, 9439181.094588943, 9446658.828845475, 9454511.05135791, 9454546.430704959, 9455671.378693778, 9456116.39656712, 9456397.468078075, 9456695.3125, 9456853.125, 9457146.77241501, 9459726.213059133, 9460774.386727327, 9461456.25, 9462055.985219043, 9462303.125, 9463882.8125, 9465464.021324731, 9492736.819528794, 9492779.333749047, 9508375.550402498, 9509081.856690725, 9513296.875, 9517017.754068272, 9518942.87787002, 9526969.779341122, 9534720.3125, 9555951.343560228, 9562494.045772677, 9567825.112379488, 9568079.286456928, 9568466.94948145, 9570962.5, 9570974.599218486, 9571044.70856867, 9572650.188685182, 9573213.878311653, 9573415.625, 9576210.83247661, 9583905.22582162, 9586254.638998972, 9587317.1875, 9587349.22339871, 9587553.392005583, 9589885.9375, 9590143.75, 9593087.5, 9603923.650641385, 9603982.965480449, 9604390.680842701, 9604546.875, 9606937.238704486, 9609308.850961689, 9609516.284902824, 9611308.823663715, 9613093.75, 9624877.759948839, 9628488.58698483, 9630603.673003042, 9630767.874217683, 9632660.9375, 9632876.838467829, 9652425.38970571, 9665323.025330754, 9679040.569819259, 9680973.498929648, 9681725.575091952, 9684459.375, 9684514.672573822, 9684690.628494522, 9685429.077217486, 9685441.636011573, 9686747.68076135, 9689137.5, 9689807.00624082, 9689820.919587554, 9691253.029709218, 9692564.080210293, 9692689.781769462, 9695335.059501812, 9695995.892970309, 9696205.681000935, 9696382.8125, 9696424.074754816, 9698119.834255315, 9698979.514743427, 9699112.510830484, 9699218.226893734, 9701579.83995475, 9702131.217521526, 9705942.1875, 9706640.441950958, 9710110.938596912, 9711336.244502613, 9754081.25, 9765270.150989942, 9768143.763976626, 9801535.397468517, 9804400.914825995, 9814025.153936988, 9817267.84673101, 9817740.324264565, 9817784.53537451, 9818118.75, 9818284.439746203, 9818298.129739815, 9818518.75401193, 9819407.417905157, 9819928.172236305, 9820739.0625, 9822303.445409037, 9822758.170460366, 9823079.964889828, 9823281.339533838, 9823735.809455218, 9823939.853611143, 9824493.070789374, 9824848.306342283, 9824912.095960679, 9825208.922525693, 9826026.147432694, 9826245.085697332, 9826267.001211539, 9826732.050847774, 9827350.0, 9827810.761088971, 9831209.961995056, 9840444.318977376, 9842151.5625, 9842727.491664417, 9854598.074789742, 9855262.771473255, 9856192.09117025, 9856197.579067843, 9872913.511049548, 9886643.50121541, 9890120.844563497, 9891354.6875, 9892342.61468205, 9898347.202098096, 9900263.739561766, 9907278.125, 9939453.994210582, 9949859.689833472, 9951289.923148287, 9984108.625222363, 9984126.5625, 10007128.125, 10007372.137489181, 10009089.614021001, 10013912.33083273, 10014116.727380464, 10015307.68547173, 10018446.105300602, 10019328.68269172, 10024282.734712167, 10026312.695555907, 10026445.702386715, 10027767.861249648, 10028839.0625, 10029917.1875, 10035609.375, 10036142.132504603, 10036410.9375, 10037533.624144655, 10061839.0625, 10063397.979330145, 10072938.61562099, 10073888.988652343, 10076427.512934284, 10076628.125, 10084051.19274822, 10116061.035375463, 10126977.886553975, 10149259.375, 10157478.169020737, 10163592.765153438, 10164237.883640427, 10164788.829404725, 10164839.416025998, 10170884.375, 10170892.13330737, 10173927.544992948, 10196498.98989321, 10203669.42634784, 10203674.379509365, 10204014.135250777, 10205352.778397124, 10207341.37849955, 10209378.812626692, 10211544.916665554, 10213093.75, 10213858.372417396, 10213872.430699797, 10214065.625, 10218065.474249193, 10237555.724615576, 10238848.987491501, 10240943.75, 10296221.41665913, 10296309.907359548, 10298082.8125, 10298438.268401809, 10299260.9375, 10300863.32957696, 10301972.91339047, 10302721.875, 10302995.495858198, 10303339.421607787, 10303477.78248735, 10303543.75, 10303545.470577497, 10303741.73410811, 10303925.0, 10305243.511840178, 10305698.052348431, 10306057.959726714, 10314070.616445133, 10324254.14542219, 10372448.194981279, 10380000.242500465, 10386225.140979253, 10387783.551169975, 10388081.25, 10389148.379518785, 10389757.86504857, 10390478.125, 10390523.4375, 10390762.788446251, 10392996.875, 10396074.99820794, 10396432.63452368, 10396526.63375702, 10396935.14668103, 10397268.75, 10397306.840768661, 10397718.299602026, 10397768.73538385, 10398174.170509558, 10398295.3125, 10399236.613046618, 10399598.82372481, 10400141.430366283, 10400161.725282285, 10400202.685622443, 10404360.728222607, 10432733.368654316, 10435409.375, 10435624.09256947, 10435714.418971682, 10449332.565323243, 10454811.76591008, 10516926.722230682, 10522560.343609666, 10522574.863961294, 10522577.445805307, 10524158.330515428, 10528937.5, 10530455.449537195, 10532096.079064209, 10534935.9375, 10538225.303373078, 10542774.469408182, 10542916.484300751, 10542943.475709692, 10544071.541011686, 10544528.155959854, 10545671.397630995, 10545790.06601156, 10545815.513190813, 10546058.815893244, 10548848.593949974, 10564817.987001283, 10568182.168388236, 10582660.9375, 10588945.3125, 10589417.934475303, 10589520.639396966, 10589757.8125, 10590440.405347642, 10591008.014278261, 10591416.879572209, 10591537.5, 10591761.368099602, 10591817.835589144, 10591835.9375, 10592053.125, 10592286.277655097, 10592458.73411632, 10592472.488433879, 10592759.375, 10593278.559160301, 10594103.125, 10594131.594740573, 10594145.168887274, 10594284.25667788, 10594623.157352129, 10595312.08726856, 10595661.06267492, 10596050.0, 10596196.875, 10596215.993059194, 10596243.22360108, 10596926.5625, 10597052.59136572, 10598710.9227283, 10599052.84061991, 10599985.416991087, 10609299.083491072, 10636650.811206527, 10636653.125, 10637817.382119207, 10638064.0625, 10640930.697161864, 10642357.277146436, 10643373.1084644, 10644029.060808377, 10645199.428642595, 10647466.784150237, 10668005.966715604, 10670204.6875, 10670398.4375, 10672816.55273429, 10673974.64164251, 10674109.26085461, 10674489.0625, 10688870.3125, 10697141.90810875, 10697733.29197915, 10698082.8125, 10698412.5, 10698554.351246353, 10699064.3394213, 10699954.191797096, 10700435.9375, 10721221.875, 10725882.082000937, 10727024.737788638, 10729781.798060134, 10730305.034707038, 10731214.0625, 10731894.777215632, 10733339.901713032, 10738946.07258351, 10739159.283389304, 10740117.1875, 10740879.45022794, 10741406.398486733, 10746292.164069183, 10749454.483458579, 10752243.75, 10770221.166200245, 10772271.862388028, 10772978.125, 10772987.5, 10773065.625, 10773746.63209447, 10774124.357733347, 10774435.66191616, 10774984.375, 10775233.165574025, 10775287.996270292, 10775420.3125, 10775515.625, 10775524.93598372, 10775842.000110881, 10775871.875, 10776688.587699244, 10776810.036570204, 10776953.125, 10777241.288671251, 10777275.84324734, 10777395.898087712, 10777742.1875, 10777834.379076317, 10778701.440969745, 10779051.442084434, 10779072.02657176, 10779553.125, 10779578.125, 10779583.61533177, 10779954.626848776, 10779995.097503902, 10780004.41989751, 10780092.1875, 10780421.799909264, 10780706.25, 10780872.910663037, 10780884.327216959, 10780944.361719858, 10781143.75, 10781166.578915285, 10781292.1875, 10781373.070245069, 10781506.983192934, 10781625.0, 10781878.045459421, 10782285.235903585, 10782326.115524577, 10782425.41484338, 10782654.218193593, 10783016.704770388, 10783053.045161724, 10783215.145005118, 10783225.0, 10783269.317205181, 10783338.072656099, 10783341.93095099, 10783640.625, 10783951.5625, 10783977.47172084, 10783984.78956332, 10785651.5625, 10786607.8125, 10786811.991781265, 10787028.034716027, 10787391.489974169, 10788310.77878921, 10788553.125, 10788565.471452907, 10798530.938055158, 10800907.817632362, 10802299.403113931, 10803359.375, 10804462.5, 10808588.674578888, 10808612.316061666, 10819733.614552058, 10822183.7059848, 10832126.537001822, 10836026.627184968, 10837050.790442552, 10841107.163319295, 10844948.135790141, 10847010.586759742, 10847312.545251215, 10864289.010278149, 10864658.09354318, 10865036.024376439, 10865141.129382402, 10866945.380682684, 10868675.593782062, 10876874.661304872, 10922466.858693598, 10953160.196094112, 10953997.715523146, 10958459.507378101, 10959701.860012108, 10960353.715128489, 10960354.6875, 10961534.50683899, 10962338.203074178, 10962462.656191004, 10962866.249768902, 10963769.149754304, 10964746.421447467, 10966654.146663897, 10967927.723956537, 10968253.125, 10970155.528168725, 10971875.0, 10973522.891605634, 10977303.84571479, 10981945.3125, 10988909.871059148, 10992332.473321991, 10995832.39009761, 10996101.117891068, 10996200.0, 11016298.367706403, 11017429.611613475, 11017741.530830815, 11018125.015003186, 11019336.846298076, 11021028.275282966, 11023297.95419125, 11032236.952595845, 11035385.9375, 11063679.574795235, 11065327.91903799, 11080902.638586815, 11082199.964970492, 11083810.99171903, 11085159.375, 11086743.629448319, 11087113.35734961, 11089652.251600103, 11090215.925553828, 11090673.805488931, 11158989.232536165, 11160354.018419227, 11167437.5, 11168110.359819282, 11169225.0, 11177554.6875, 11181131.141750481, 11187539.338001879, 11189665.625, 11191737.236550096, 11206631.56392146, 11211884.647520201, 11214126.620531777, 11295855.643720938, 11295979.6875, 11296564.923850661, 11296744.506728716, 11297624.610689642, 11313015.414485484, 11314030.550396627, 11314727.885854278, 11340253.396617088, 11395433.91478568, 11403448.114480177, 11492980.7912614, 11493967.76805489, 11647892.585424425, 16228155.364985706, 16966407.405190874, 16966407.735301036, 16966409.26151319, 16966411.043223858, 16966419.23521041, 16966424.215175677, 16966426.399346806, 16966426.474611294, 17019036.096217223, 17026679.6875, 17028646.630960144, 17028732.8125, 17029206.70469579, 17029578.92315375, 17029942.104518156, 17030474.429124456, 17031115.54155804, 17033077.043456983, 17033225.618404295, 17033288.33227493, 17034658.879324015, 17034967.1875, 17035715.3996179, 17035835.29785716, 17035968.75, 17038375.0, 17039714.693723567, 17040379.539999146, 17043613.4145505, 17043743.75, 17045100.0, 17048337.5], [10.065459696394544, 12.326574375071807, 8.159669690351624, 67.90080682193609, 12.673214011517171, 78.83340792819155, 22.15804401485982, 28.235848172362726, 24.376230438396476, 6.024202519815038, 41.934350152355144, 6.6125498325183, 27.336288191385698, 7.552640313244892, 57.18113545269022, 6.189297189141456, 31.063538671754323, 17.95679062663734, 15.339751575277347, 33.718717044999835, 76.87883709580312, 58.841636439642684, 12.215641638910325, 21.037222353227897, 61.34892485615109, 94.39492082299645, 34.968065838312704, 19.519983544344502, 53.713227745639436, 13.859576650916187, 69.1274456075273, 17.01002492809936, 57.045100904849875, 19.543815337145745, 29.255540507501614, 54.02718570924705, 23.285083349427197, 18.555261184776526, 164.49361023456058, 96.27009658228715, 96.51793925946248, 87.73290846818036, 42.94959745352294, 50.0104046091421, 9.454500010236732, 6.178756896108188, 5.302864789419282, 60.99320441044666, 123.04106404617539, 5.394826747899507, 56.02721212423233, 34.221435534876804, 13.159289941202506, 36.18866569165773, 20.47218165547233, 70.4145057115702, 23.511384517176914, 20.119371644724723, 49.29620062267678, 115.97819889059197, 76.28279020347254, 71.18256387459354, 20.906135036125214, 26.05840568469945, 78.00236695436391, 5.1772996169510295, 10.756801979318887, 117.23148808335588, 20.3920821489323, 27.85495734193822, 18.48687371385993, 33.46518069902873, 13.747328481721174, 57.1806240353857, 138.5696484394292, 141.78809681086875, 38.62134949131011, 88.10219637428074, 9.705975125328253, 20.2363178140899, 5.815905221489664, 122.35462834937749, 60.6242781811777, 39.31723922460627, 17.324307946503897, 36.66445185487366, 16.858542597302872, 5.197819362108079, 11.642109655699372, 5.033871245317018, 11.759753899138556, 160.5072730740715, 10.522827436614566, 6.529419460729219, 59.70910964443793, 75.55180408811208, 11.499017128810143, 14.619315347952877, 24.06474157228263, 5.027469132053725, 7.247404026985484, 69.4093628626539, 119.9425404591059, 208.23390962816254, 147.80559523951146, 275.8375934686521, 59.76773777879205, 20.433967735746013, 9.24916410158467, 91.07225708151792, 14.535181632232376, 34.506380310638164, 15.47946640011782, 37.77053765525138, 40.039300725286665, 7.1653715294197005, 64.64485516763472, 157.8528911155728, 99.88735124733068, 26.31051016322639, 232.62944169638956, 10.475661248583219, 194.03580927282368, 63.47737823014859, 64.0177193969828, 6.85836347983801, 48.723145962181015, 27.116895715423034, 113.48950163865842, 30.865936465142692, 23.990290725278545, 122.33336364120123, 106.70438197321933, 25.556181183450228, 83.28954850096201, 108.73852395216578, 8.269566281733962, 29.717047787788154, 14.585252213038318, 38.30902162389954, 28.123613319398967, 52.873992016953096, 46.40061248258097, 10.661139281559956, 19.2063892367583, 64.67726967439827, 89.10508649879478, 101.67991738707937, 92.73435444361691, 69.43333169784071, 84.0124461749939, 97.86917353449624, 26.674556670267737, 83.1773630397125, 8.001728409237012, 31.864799635065236, 99.00394800279389, 16.968792435542497, 47.83386924890513, 24.238586174245164, 19.916205608514428, 118.04414115898064, 17.698062327434872, 15.937448260462798, 188.2653959157293, 7.177113463957441, 41.58661621733429, 6.848276235933935, 5.734491269605495, 35.510178982449574, 56.23593188966842, 205.3470599882475, 71.46467667378529, 152.57855942250472, 27.710727269360188, 33.43903012921928, 8.93597250843404, 42.720581694753115, 57.31326136866095, 22.731773109874915, 17.95266389864998, 16.301969352994643, 67.51546792021543, 12.263950379629946, 12.496519965979616, 81.46236944322362, 11.29522535927037, 29.213805194518436, 9.628022280951724, 223.57620509587545, 18.274314696791553, 30.03376882115046, 81.87519684522694, 12.973413642526028, 32.15820699348023, 130.45946612629496, 33.942662110412485, 145.635519367409, 7.809169911195373, 69.77107079458565, 27.902277642309922, 142.2131935205544, 79.71398832422027, 53.8838063839399, 36.62525566514588, 32.328151625181405, 6.066160129017725, 5.197894504185542, 9.875113520031144, 12.762697001482076, 31.3254844012218, 315.80796899854676, 5.533877222599556, 19.538112312733475, 86.63581119026658, 23.82456708047306, 23.01881321963276, 24.259030202183915, 11.591118065352681, 6.595205822988215, 66.64062821332942, 6.407640974682039, 62.00215721879985, 66.05933862480246, 40.426794174372844, 32.61087589486444, 137.3036201013023, 8.078970779357652, 44.05283632857133, 69.42436241361507, 139.57198843343193, 32.20323066491779, 479.60358455066694, 49.27135503473501, 9.874451174467843, 121.73531173395843, 69.58737526305512, 5.670386154241651, 137.93884700119983, 6.989574262155711, 105.48470569001742, 52.41244420290899, 6.698849362435621, 7.355789738177282, 13.039597944474542, 80.22688052604067, 28.910954195080766, 17.795463977158093, 16.935344048199802, 47.471275247959596, 110.49523320644775, 75.89619257546168, 5.264717220175217, 70.27631937158931, 14.227596071218645, 43.38274343953543, 65.2639489197569, 15.38720946124976, 10.379315583679716, 34.57712925055541, 27.892364691119386, 98.08434940176893, 158.15084275439068, 150.40985011948044, 9.92824594925841, 13.24961199017178, 125.10850388964877, 9.253644913489001, 61.74350832904565, 21.312926313502608, 168.06902980265556, 39.46650572302364, 29.242406965757322, 25.59309003571559, 91.99012681649751, 42.626109500489285, 13.226877837918321, 55.56789156222786, 99.15619445953797, 60.110690088661094, 23.26775324595096, 18.95882675822094, 17.839463591663954, 20.150888427470214, 17.978768316925354, 16.430335044549018, 5.992982945100992, 63.76914484638368, 85.97126911898198, 20.6908408946346, 37.06738855120284, 89.10996031127964, 28.782809981411603, 110.3937149153673, 39.61050126644795, 279.91825738110845, 101.6387698961279, 8.395657668498057, 5.769360846555785, 25.400363662995677, 23.23301823929342, 11.634916322668293, 48.356174667413796, 165.50736595347345, 38.58038692391099, 118.14936549096926, 68.21676231015972, 63.24843990277795, 10.297236993788147, 7.6604414773303, 16.300830201272582, 20.549848923035157, 20.00026263552206, 59.12997577296272, 11.237816453166808, 6.626457273464917, 18.0372646034879, 45.47213889630851, 338.66284978512766, 86.72734122303476, 16.0370270911384, 25.821890802518194, 73.13304793216244, 22.355939496397795, 6.165096365903978, 15.336425672218802, 23.190540028702983, 6.173524488040721, 6.127554851299917, 28.564689700970696, 90.70776920376426, 53.378220714067396, 13.018559156914483, 11.922167191350615, 37.06442693598332, 27.983291031244303, 128.23797468473612, 17.749849216711628, 159.7344748869917, 64.07675349860621, 11.804764732056647, 79.57103909975811, 19.893197669174945, 145.63421374644633, 9.60259526611847, 110.47286501984465, 14.117279831033876, 78.1129848380767, 106.9369372691047, 14.934163914400377, 72.66787981357452, 66.46051691694319, 127.10221219197, 73.44920303039733, 36.37537185311975, 41.780700034748875, 66.14160531688321, 121.67854605858096, 39.640005971909794, 72.44008334460138, 11.656337448957867, 67.27050074006516, 85.51835470166353, 28.25052885260542, 37.73617390642234, 16.390418778497622, 14.291284854072707, 121.2407542868341, 174.3094470280064, 5.477410906521663, 73.72028038674958, 62.22240867066296, 33.095014466396385, 11.36912875067891, 9.377324724586217, 7.2073987790956755, 5.9478073414078505, 87.03264983686122, 179.72046613126642, 27.141477073248826, 69.83435110743261, 26.66025927609796, 18.455824200760375, 18.029648430908413, 31.628242956647103, 42.58030262632378, 31.476992650902265, 8.938655988111936, 15.627165870220361, 6.003651321862843, 19.59925403533, 25.457045490189284, 134.9695500164637, 36.329553137335836, 24.994473797092187, 42.71499633458385, 8.184477586045038, 59.60599125811213, 58.54091170179636, 84.93361830342135, 24.634654765351552, 26.78600116319455, 24.364745470245747, 155.23175199740942, 80.98047339367035, 34.432593171301136, 15.393998790065904, 27.447998714878107, 6.668309169395795, 7.340203695421902, 69.35328828634582, 29.17519662993057, 159.0314774688759, 77.184407289088, 44.92183480620913, 93.92492811851584, 50.38908036339462, 14.01998398779195, 61.21445045635694, 95.34837126174128, 61.93788664363602, 10.208108917672467, 13.745266441722075, 136.30732278817152, 8.635766017497307, 46.41499155280363, 46.55713767002439, 51.63752346378877, 31.051016118110184, 56.081423998232935, 18.739290541337922, 16.202950061627664, 14.217057326011915, 38.58277585856864, 76.46731269646023, 108.82405792402594, 26.15026590113674, 20.51541097116636, 5.956007547241795, 30.916344647362077, 149.4839133046695, 30.502535883437528, 12.088132908065273, 27.600288351746094, 57.08060228599658, 6.931248623199689, 8.43300841421391, 105.8130201421919, 116.19300313285602, 73.3297633313933, 15.821300020237416, 8.369167356200602, 15.158714989085617, 39.48771082669069, 20.1309020098496, 16.39901706671206, 18.238767426921566, 20.99583613804402, 102.11129709086629, 77.65320227576704, 29.06666785773202, 74.12603212756545, 69.41010720163604, 52.90979270018388, 34.81196920966274, 47.12058290840053, 50.67416763950388, 11.042537421217501, 103.12659872798778, 12.163743139299722, 71.60845435651234, 12.237057728643467, 14.954845654235575, 142.2898355441057, 39.21564365064953, 8.764197323283433, 53.57715230122585, 14.367513716039007, 8.636014621737234, 19.00258351596674, 214.7182086708191, 33.9593156473441, 64.03086381255758, 128.00294858978313, 10.341310038781456, 35.60166793607071, 62.889541442550914, 7.4389627404447465, 100.67797950943924, 5.9403492647604885, 27.936928293769004, 93.60228502830421, 56.60210916048675, 52.440015388456644, 25.36916149703492, 6.267336964182166, 142.86993435731887, 77.48843123356143, 25.697109967647346, 89.42970172773758, 201.86661988401644, 52.417022386164064, 5.622617448234552, 30.615128794190746, 51.54192424453679, 32.03288908526014, 41.73177361975442, 15.194171078163473, 48.56539518339554, 61.89643205659379, 121.60747512517852, 5.292548587274148, 47.52382290839143, 16.32169601344293, 19.109159741396557, 9.577508153934511, 6.5012186363983036, 9.225431071071085, 19.154820962656537, 70.57575143040913, 26.176998967952287, 32.18661949449039, 17.151401885736995, 36.39151453257347, 5.197681925558865, 69.89341217299983, 40.84599490468364, 246.90893509464712, 89.09501314793873, 28.72664120378201, 131.95434371899222, 91.21018374741894, 6.64332676985908, 10.604448881328048, 55.26382310678627, 14.58663186725983, 61.20125539035429, 97.20093871987211, 24.464620053201173, 8.589519297009351, 76.7025913549202, 84.63074616163107, 7.725988783089798, 32.044679196995645, 104.69502620111811, 6.183407453262966, 301.9032276027648, 68.23134702754746, 19.31740992465969, 85.22836757300185, 39.63781040147598, 78.05571346886943, 5.361182281697263, 70.42191858798284, 10.477379553265456, 26.24395950377896, 7.917570078999962, 283.29435942754867, 8.245865770958975, 116.73378282189418, 16.66974794586078, 16.28322935215401, 97.02527452853384, 93.2037563919223, 44.12266641989035, 7.5834671629781525, 79.47055430307769, 72.70418532662393, 88.76014787255907, 15.469562506365564, 93.0984403950143, 27.144503354023733, 66.3639181817211, 45.391404108765286, 106.26864580433606, 15.396465654520192, 58.746573210664195, 72.30458406472137, 9.368904645548543, 17.720125862958486, 8.50043622373016, 24.644709058910614, 5.729632911506516, 72.183788932704, 32.56472921881566, 20.730464268369563, 52.89002266113788, 178.0563791885188, 47.14942698220819, 29.621451336517723, 185.26046426300917, 18.880115192690663, 12.404593359390443, 22.599198691572123, 72.6465322524311, 383.32416430876106, 99.09549650019648, 95.86189969155802, 69.70951041318037, 43.17127263607583, 11.624529203698046, 104.55858651647182, 52.89531522158302, 12.076010612986154, 83.28068194113732, 10.078280218566665, 45.40544323654357, 19.44671083833787, 6.663196572492929, 63.24347467103084, 9.491409282598319, 138.49925508672797, 7.483882707337109, 38.384597204178505, 15.58032004392454, 38.53712379414376, 62.94225979247324, 33.62344592642402, 17.71585488700972, 48.043733934347266, 73.04478714127971, 13.212263432792328, 24.51384567024685, 15.95405697925217, 90.35402324414531, 15.7985270107426, 70.28664657480216, 48.26671124703293, 9.376234752126175, 52.724639014308906, 56.00612547973738, 16.168089650258885, 41.089293710678675, 7.6025896541126095, 86.44607878882387, 15.220246285385699, 50.23777935004686, 32.31150764119185, 17.800217993589715, 58.98211207271988, 28.2388441793577, 34.49001661001151, 55.849097805765346, 18.341996323764768, 17.79121592841217, 12.672977159194982, 6.174040860273673, 69.53184609556901, 24.983161194525028, 56.25727018970835, 100.9915571042396, 139.08525438361283, 5.0754364483683325, 141.6603488874823, 96.13179946636197, 15.587407655527041, 5.669975746392192, 14.447163606509626, 8.8163261594362, 150.98989757725772, 53.61217849057084, 115.73764120992095, 18.974374273025973, 29.537627376195985, 16.239331841585514, 13.243035902456594, 93.90975460366869, 13.989858136777631, 23.326084413973316, 11.176446342871463, 362.39255354598663, 10.353568338503006, 35.148138151397305, 87.66175691216887, 8.371080529492417, 10.113239952736825, 6.96183792526149, 135.0409280951333, 146.46341861609227, 125.41377121128748, 84.93206909539391, 59.28490254738445, 5.920502442664166, 114.07319500294477, 13.884781643939055, 37.164875501287774, 9.724125692667407, 23.428802564744636, 31.334570732072123, 53.9847987377896, 11.704551526353407, 83.99686123483983, 17.940913422183446, 41.28985228945841, 6.8533933084314445, 81.20282521392664, 27.286396885378352, 25.97749289303358, 26.3864870672791, 24.007681170223993, 65.35252547563826, 40.42881319245795, 52.45053160245735, 76.07737497511714, 5.972112315864063, 60.594771615975816, 50.15551984015993, 7.420998016902543, 28.2358327752387, 13.958097751026655, 13.270225723141753, 5.749142163068124, 16.965340944839255, 58.98720347255606, 33.15442047556916, 125.52647676316052, 47.98119672721745, 32.49420622901363, 35.62810355048529, 6.644525704252957, 8.890893967153294, 133.73462616344085, 37.27390272299792, 67.94663672529002, 99.82435922776703, 30.909821652988278, 59.95325817888908, 26.455524223267073, 59.75990833382976, 88.82381543047235, 42.35341226897886, 171.47926504204983, 110.1623147324641, 66.08631837690804, 23.298440263118618, 98.03271895407283, 21.536073437454448, 41.13782277900193, 31.35816424224633, 30.78355607556629, 88.62757950465914, 17.542859718756944, 28.183194195654025, 36.158604397184504, 39.19304850130103, 23.696178202837846, 301.12955990407653, 122.82182654189265, 67.70194188011621, 24.58649720745498, 5.622311939638599, 6.264242584517828, 25.059082512659224, 113.62676198197086, 185.30195965643145, 117.47350733752184, 36.40414025937308, 62.10297631053583, 32.97330868378769, 47.65373627319608, 86.61934269209016, 42.72068977812477, 50.109052070575046, 16.482853237803578, 145.58825427444586, 40.33743802688613, 82.02927139611825, 27.31915552584696, 44.00442719645317, 47.50431005778305, 21.748213995091653, 9.961345667156532, 51.047474997846194, 46.29669297096342, 24.73272789478049, 17.7684646626583, 29.043442417272896, 9.33485231090715, 19.62129069750685, 10.944565181586293, 11.37318681516897, 51.95682191914705, 55.581556185500595, 52.22554472494025, 7.581447488853529, 24.90087336217376, 22.60715365755828, 65.6785166311473, 39.7517630450412, 32.55035642590353, 100.4412471799809, 27.613464713728156, 17.93017086130963, 45.08593786324791, 7.397017692495815, 61.417531761955246, 54.891105686862744, 40.45820534823959, 54.73987180374691, 14.646763949633657, 27.627511104933678, 104.75154875608781, 119.04327473969308, 5.276807539656052, 13.373595747806007, 5.06744735110097, 13.381171313498967, 204.94114190167687, 20.272921093612005, 68.55689119311319, 14.046215243124289, 99.85135948458216, 45.555620558410055, 8.8340225818876, 6.554102487354428, 40.62616220040944, 31.4115098150593, 118.0483098395159, 26.85837008469606, 25.22798335550291, 5.161681720697702, 40.47379749839956, 14.386867990922868, 30.97194167445604, 48.854923975740114, 44.35617078161087, 36.06529695493804, 54.651992854470954, 10.571059908188062, 61.12192188125745, 16.7178819683983, 185.83646826845256, 49.31376504909399, 190.1516808460131, 81.67649128028869, 109.88853086326847, 19.667635845824883, 99.23130428068568, 16.148365717740667, 80.82221162082917, 12.789701790520379, 5.100169953670094, 104.94380869176253, 59.44527635606772, 72.54733814399464, 117.3186700377615, 30.88603038998292, 174.0046456231751, 12.75661265184694, 74.77876702923287, 62.92098599606597, 17.79628633064022, 57.86146682378606, 112.41108283832766, 6.60325717437112, 9.467058199607353, 65.30352443305922, 123.5737953727717, 21.936016525223863, 34.37039282683617, 46.15017625309061, 56.629963079487304, 134.1043870097884, 47.73172656438225, 21.047178710776805, 12.478254968234028, 39.162956909214635, 17.242486315825293, 23.26287267313329, 29.174765264560158, 148.24773162620716, 84.310106893874, 14.681998524766712, 38.51802091207599, 6.639125960566017, 33.27205361564582, 57.056555491201664, 6.103868004718433, 10.447885066075768, 48.79289887500329, 113.74595070618177, 6.7412126123920775, 25.03116931522865, 249.31740185213908, 27.8500276297756, 15.152274044739151, 127.32745324234878, 154.50425460083176, 19.664709440802774, 8.688018078395125, 9.599161313075442, 55.073017372851055, 9.479439374646871, 6.186210584107579, 101.40786959246094, 20.751680855190685, 45.99453428874781, 95.46101296402868, 14.544364039424249, 16.106996110291963, 54.28799929660118, 24.229122620054678, 8.459866719167694, 9.572978789467573, 55.37952011887972, 57.742413330089306, 73.50700039661027, 83.16175028438187, 27.1766141658076, 119.10774957498592, 24.96530416807874, 41.11463218909563, 7.059975158322338, 176.53832848186468, 78.05949088325752, 35.11174603679402, 21.914822164344994, 178.25477306200423, 13.466039325466113, 8.680577160422564, 197.06124588064873, 131.70261997563264, 139.8993133124308, 82.90671994441763, 15.50703962848715, 5.223698068289605, 5.240601799538259, 5.201188546917007, 9.469312847282577, 37.02017232787812, 10.820932184672952, 69.97585055411137, 56.501639957575364, 15.337162307062552, 8.243177265278701, 15.317181192462161, 5.807831964779528, 11.777096421987164, 5.649197026536696, 19.467615982095005, 21.53634330521458, 79.3231112237834, 24.753091962443133, 6.189630701765054, 89.14052766737274, 103.74618285607986, 12.282053091461476, 107.53580062706234, 17.723495552468894, 47.96542030644419, 50.29343427466593, 114.90962770801391])
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)