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 = 44865
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);
([6743402.047899941, 6798381.602915615, 6830177.945460082, 7077994.625212399, 7158199.763187827, 7282689.0625, 7290009.911801992, 7292031.267858106, 7311199.491457405, 7322373.542119066, 7325259.086153767, 7592617.9120263895, 7619414.0625, 7650001.252156334, 7818962.844233887, 7892179.3971398175, 7993259.776386856, 7993287.783796568, 8061903.125, 8099346.142670247, 8147526.255549074, 8191301.684435399, 8198802.344454739, 8243929.6875, 8371743.8449094165, 8559257.8125, 8565608.889043516, 8565949.27115848, 8568086.754444497, 8666300.048107455, 8673133.781371502, 8673821.875, 8675409.375, 8675463.81485712, 8675915.625, 8677207.594032718, 8688659.316214968, 8696385.935082313, 8795041.94489316, 8806132.534069385, 8819115.625, 8821127.23927495, 8828546.705870736, 8829089.592235776, 8829220.676252192, 8829486.473690713, 8829860.193799892, 8830262.433225669, 8831203.125, 8831284.06113106, 8831308.088159239, 8833877.79661033, 8833913.354696944, 8837296.30220899, 8862834.578185035, 8875394.782318195, 8877307.709458087, 8883792.59997623, 8885794.745985812, 8887077.432332234, 8888476.144987103, 8889162.532126136, 8889599.69769525, 8890515.625, 8915554.421248462, 8916862.5, 8917558.547532491, 8919004.91653882, 8920257.707815751, 8921026.523484185, 8921047.204317106, 8921781.25, 8921909.127004456, 8922135.82795166, 8922645.249725224, 8922682.223801712, 8922697.036171935, 8922722.321188247, 8922774.338740181, 8922932.971570926, 8923010.224360788, 8923025.262222094, 8923128.92335266, 8923187.21360701, 8923301.38280371, 8923522.226092564, 8923804.6875, 8924259.684525425, 8924300.0, 8924482.8125, 8924681.478319036, 8924844.117997319, 8924882.814891769, 8924952.352818947, 8925832.459919259, 8926000.839751719, 8928321.875, 8931793.671313053, 8931845.1362293, 9005839.434806848, 9010785.613939468, 9037091.566264534, 9079485.978516554, 9113899.469612477, 9129703.78641302, 9139654.0289821, 9141640.625, 9150921.875, 9155134.217347628, 9163389.32856721, 9173979.265661016, 9227791.578400794, 9229456.25, 9231723.4375, 9233507.8125, 9236501.5625, 9237008.292103166, 9238279.027579477, 9262815.625, 9275746.875, 9278123.859842023, 9290193.365947057, 9290581.21718795, 9293156.286744399, 9321256.25, 9326245.896215344, 9361627.503131231, 9363169.825398516, 9393771.668867664, 9394259.14901842, 9394879.326302538, 9395200.344477037, 9395526.145135935, 9395536.79565682, 9395728.77775762, 9396083.838617431, 9396297.465110872, 9397680.187797802, 9399192.966858855, 9442300.0, 9448465.869956706, 9449515.625, 9449869.31972607, 9450339.547300255, 9450360.931962002, 9451599.106036844, 9451673.152522145, 9451799.269565966, 9451939.555731645, 9452510.216497613, 9452637.258051476, 9452704.6875, 9453697.93623441, 9453769.84482801, 9455366.549606072, 9455778.125, 9455916.153884234, 9456256.112203613, 9456472.50418881, 9456614.028758012, 9457790.084885208, 9458837.79887042, 9471321.49382651, 9486562.793315124, 9487148.133209087, 9487582.8125, 9488648.308595479, 9489118.789938202, 9489733.056391979, 9489857.46482766, 9489869.261599962, 9490036.282599531, 9491150.0, 9491831.25, 9493293.724922678, 9493450.600779254, 9494657.512572423, 9494882.408332054, 9495812.5, 9497845.150561508, 9498379.216677101, 9511203.896424336, 9512483.749137603, 9512506.216105068, 9512577.043258594, 9516327.24100952, 9517293.75, 9520831.509894272, 9520838.100977363, 9520887.285703335, 9528815.049227403, 9541395.956358386, 9547688.11173419, 9549129.44297439, 9549323.074040601, 9549505.7887901, 9549530.055542015, 9549634.130297204, 9550012.5, 9550361.336093327, 9551120.3125, 9551311.873281494, 9551546.875, 9552032.8125, 9552039.718232624, 9552156.472794136, 9589113.481928477, 9595620.3125, 9607625.464300469, 9632885.67079013, 9667444.472582666, 9670919.70969289, 9673465.504699042, 9675728.38105047, 9676754.758390075, 9677471.566221675, 9677641.649909917, 9677748.295663418, 9677812.179909376, 9678812.615135925, 9679031.162448786, 9679152.133787747, 9679246.205583414, 9680095.3125, 9680585.689703846, 9681243.220494457, 9688088.874592062, 9690758.949480427, 9693903.621325847, 9697162.155219382, 9697215.625, 9697216.35321655, 9738272.905085005, 9755259.960072728, 9761356.064686337, 9785026.695268262, 9813548.4375, 9841014.732487738, 9844292.041361872, 9844482.576824209, 9845581.38642707, 9847800.0, 9848280.67154471, 9848316.891865669, 9850382.019590857, 9851938.06102836, 9857828.681101697, 9861078.021794207, 9864912.417904591, 9865471.800352743, 9865520.913177595, 9867945.564919032, 9868532.8125, 9869421.274775546, 9870286.863985201, 9870940.969802054, 9871790.262409143, 9871904.307981879, 9872725.247330341, 9874207.354109738, 9874225.0, 9874258.334242392, 9874870.624736222, 9875300.0, 9876315.355329145, 9876456.25, 9877753.965418749, 9879672.929685753, 9881135.248643547, 9885026.649101434, 9887716.91539338, 9905766.364989016, 9918983.555378642, 9919048.707737828, 9921847.546116821, 9927783.084713886, 9929404.6875, 9944340.625, 9944347.817351045, 9952073.251897117, 9962678.125, 9963928.125, 9969267.342571348, 9970889.888190618, 9972737.5, 9973657.228670252, 9975851.5625, 9976811.151044188, 9979032.8125, 9980025.725869425, 9980478.221923847, 9982718.59516545, 9984211.652536118, 9997553.125, 9997617.1875, 9999187.321820846, 10001132.8125, 10001358.386192853, 10001573.16223938, 10001589.285662187, 10002332.8125, 10002433.307225011, 10035748.874088014, 10067467.887844063, 10081029.325270029, 10081437.864537299, 10082385.97451424, 10084497.20998952, 10092770.365978505, 10197839.40582428, 10282825.04219514, 10421987.5, 10551493.531040112, 10559381.248354742, 10645693.75, 10645712.483247329, 10646298.084241884, 10654005.623125719, 10663141.278514914, 10663745.219320936, 10680872.697373778, 10683051.563900454, 10683396.863421177, 10685053.073056784, 10716763.601832153, 10720984.375, 10721437.237904074, 10721481.84350176, 10722329.6875, 10722676.188178983, 10723466.005365878, 10730434.652022725, 10747680.888836363, 10747683.467603015, 10752866.650712062, 10754037.523062088, 10760308.497057369, 10779490.625, 11016616.564249625, 11017715.973472947, 11025346.652957622, 11025805.277069805, 11034679.6875, 11046143.006089957, 11048287.491699563, 11051360.727217836, 11051699.564975418, 11053482.595030405, 11054983.59835332, 11055007.132792914, 11062893.897731729, 11064693.102195082, 11067666.872203907, 11068452.695268223, 11069261.552363615, 11069541.43331458, 11070129.6875, 11078495.3125, 11079040.509224212, 11079910.094702065, 11079941.578790737, 11081456.25, 11081509.861383995, 11084193.563882072, 11084314.48643726, 11088002.530185007, 11088805.33647758, 11094739.948669706, 11096678.125, 11097071.863696026, 11097179.65920913, 11099583.971130332, 11100087.5, 11101771.157947425, 11103187.00023111, 11107595.845638173, 11109972.192455528, 11111576.5625, 11112555.827417359, 11114137.419940941, 11114871.470160896, 11117139.0625, 11117921.875, 11118027.616111327, 11118238.095704166, 11119807.259347642, 11120008.678784892, 11120142.211577652, 11123608.217645627, 11221038.253414359, 11224935.422531042, 11227424.675069243, 11233421.535170479, 11243921.40356013, 11244980.222639656, 11245903.556423357, 11246008.37408917, 11246260.751260443, 11246783.880057987, 11247253.105925048, 11247388.761919646, 11248786.09665295, 11248786.599453673, 11248953.721059587, 11249387.27468009, 11249557.202176694, 11249780.585204218, 11249830.926156651, 11249923.499364562, 11250673.4375, 11250926.688847262, 11251054.945157798, 11251100.0, 11251136.504450615, 11251288.413683478, 11251642.1875, 11251689.234674538, 11251873.4375, 11251969.15617041, 11252055.153007187, 11252314.302013319, 11252725.361950723, 11252806.25, 11253150.0, 11253470.892550323, 11253715.215855062, 11253962.222325455, 11254337.588814521, 11254529.6875, 11254547.028175132, 11254753.125, 11254809.313665938, 11254824.89800638, 11255043.755572073, 11255087.807511289, 11255312.026655622, 11255341.25025888, 11255354.053990059, 11255478.125, 11255943.746229323, 11257713.00638248, 11258909.058294, 11259077.645152293, 11259133.89665048, 11259226.714597218, 11259404.53347181, 11259458.467007907, 11259509.375, 11259540.625, 11259549.16405323, 11259929.6875, 11260004.6875, 11260180.027693259, 11260332.237797504, 11260428.125, 11260445.3125, 11260452.04862387, 11260523.190944107, 11260563.988868613, 11260667.1875, 11260690.625, 11260942.169599727, 11260947.130013453, 11260955.401367022, 11260993.75, 11261171.432129772, 11262102.016553279, 11262173.05804178, 11262306.957485372, 11262389.0625, 11262413.235523844, 11262451.5625, 11262582.231918441, 11262621.442216007, 11262706.93482574, 11262717.098555, 11262783.211645022, 11262902.169929666, 11262916.03506481, 11262981.25, 11263214.405164868, 11263391.079617755, 11263479.168062484, 11263535.32757945, 11263547.486920984, 11263643.007283958, 11263673.236488672, 11263673.4375, 11263720.3125, 11263926.5625, 11263981.527958298, 11263986.472287355, 11264203.125, 11264241.716405364, 11264510.03419306, 11264737.933939505, 11265021.231589755, 11265325.0, 11265350.49425418, 11265734.066763027, 11265845.120571652, 11266047.576086733, 11266104.191502687, 11266123.145713773, 11266205.738694292, 11266508.367286494, 11266551.5625, 11266688.976482674, 11266749.200895464, 11266934.998156419, 11267168.75, 11267518.265620317, 11267670.030333158, 11267837.1363112, 11267954.594134295, 11268284.524032107, 11268362.5, 11268856.25, 11269122.941708872, 11269447.578494463, 11269635.644399654, 11270027.84429153, 11270029.972819598, 11270209.738215491, 11270540.625, 11270912.5, 11270936.545741724, 11270953.694681376, 11271130.923449667, 11271573.296061149, 11271637.176577523, 11271660.9375, 11271697.017787416, 11271758.690429682, 11272092.605731577, 11272304.79168194, 11272404.6875, 11272722.531271096, 11272748.4375, 11272758.521600097, 11272879.540220736, 11273017.880245928, 11273539.0625, 11273747.37430062, 11273860.686995288, 11274009.382949473, 11274219.567241652, 11274394.887295455, 11275032.825021883, 11275064.091528695, 11275096.875, 11275519.695452455, 11275656.021861097, 11275784.375, 11275843.138644956, 11275849.286556767, 11276034.042185161, 11276076.528115215, 11276458.79493328, 11276852.366278702, 11277007.179686116, 11277195.112426566, 11277610.1557915, 11277665.013559746, 11277814.764031444, 11278172.644000186, 11278181.25, 11278524.588094536, 11279001.5625, 11279029.475510178, 11279175.348068882, 11279477.413664622, 11280029.6875, 11280276.942770876, 11280307.13760438, 11280514.0625, 11280612.905062014, 11281131.590039296, 11281212.859105937, 11281338.0254681, 11281939.817494188, 11282253.822148083, 11282336.208407873, 11282506.770783184, 11282890.300723907, 11283038.545070168, 11283476.5625, 11283593.75, 11283860.537690671, 11284084.685542293, 11284192.825551659, 11284390.492906354, 11284569.620365875, 11284636.886096908, 11284692.306522377, 11284793.258998405, 11285271.88230377, 11285352.090823002, 11285438.05660037, 11285508.18881671, 11285784.375, 11285937.766805628, 11286274.963756703, 11286323.51819735, 11286358.865631655, 11286564.916849026, 11286599.70387126, 11286833.56318754, 11287138.385939635, 11287366.314209733, 11287371.19734061, 11287404.059194844, 11287863.857593374, 11287970.46836263, 11288402.829492657, 11288639.261809329, 11289400.696880508, 11289941.01555833, 11290562.698667949, 11290684.252035806, 11290841.003296647, 11290843.265150448, 11290893.16664481, 11290932.8125, 11291570.277706416, 11291614.0625, 11292246.237497848, 11292673.137717148, 11292836.15532412, 11292890.625, 11293444.09745987, 11294413.966865301, 11294434.375, 11294505.420791708, 11294600.378517732, 11295079.199046861, 11295170.3125, 11295379.073242508, 11298315.434642745, 11298508.191243798, 11298673.740197107, 11298751.363568773, 11298751.386331443, 11298869.479272414, 11298884.2235384, 11299041.57737906, 11299295.3125, 11299351.225879807, 11299458.993460516, 11299513.968991028, 11299763.922689699, 11300493.1337465, 11300568.43669597, 11300592.110096894, 11300617.023694413, 11300625.0, 11300800.892578958, 11300831.743397139, 11300859.013451736, 11301268.75, 11301432.745423332, 11301527.880178917, 11301624.840088332, 11301831.69121393, 11301915.625, 11301924.584974673, 11302039.610722868, 11302134.375, 11302412.443471853, 11302415.25509626, 11302803.125, 11302905.902223945, 11302983.247437404, 11303037.5, 11303161.182508612, 11303168.75, 11303195.623084964, 11303306.904897362, 11303314.342610575, 11303530.968579883, 11304043.75, 11304339.98514398, 11304681.801956072, 11304805.05935668, 11304808.827799581, 11305057.8125, 11305139.0625, 11305198.77722968, 11305401.5625, 11305468.500346933, 11305502.99305331, 11305945.987687774, 11305960.5916688, 11306022.521296643, 11306342.1875, 11306514.9141842, 11307116.655777067, 11307154.6875, 11307454.352491166, 11307526.97690535, 11307893.195201531, 11308052.705176538, 11308634.69348757, 11308806.767206019, 11308979.6875, 11309031.25, 11309154.865051925, 11309256.538312335, 11309265.625, 11309285.729680665, 11309630.99134203, 11309749.769177794, 11309964.0625, 11310441.109902557, 11310474.31848592, 11310761.92346772, 11311055.70757198, 11311094.08640303, 11311220.146440255, 11311239.949898243, 11311400.687187698, 11311829.626127277, 11312009.156366255, 11312634.53075563, 11312836.223302364, 11313051.087343905, 11313079.654864563, 11313247.215488696, 11313331.025665622, 11313492.1875, 11313854.733233238, 11314059.27190385, 11314092.1875, 11314126.482705917, 11314179.627428213, 11314233.009305712, 11314252.218157006, 11314373.4375, 11314429.490514433, 11314443.75, 11314546.875, 11314596.279876657, 11314676.714352708, 11314685.053204091, 11314808.761223745, 11314813.657872908, 11315124.739334922, 11315214.0625, 11315249.6099362, 11315284.039979445, 11315345.691438667, 11315346.875, 11315490.181759084, 11315657.804461807, 11315731.088141968, 11315746.875, 11315757.503182415, 11315885.428975312, 11316336.740289332, 11316342.58624331, 11316685.557988869, 11316762.282451184, 11316893.805212267, 11316953.606789399, 11317383.996715065, 11317508.025802061, 11317564.327525081, 11317939.0625, 11318050.499148585, 11318467.936054807, 11318557.37125582, 11318573.49468094, 11318650.209601376, 11319079.602021655, 11319407.87428653, 11319507.238675147, 11319531.084129144, 11320023.86007344, 11320290.59555386, 11320511.747091891, 11320885.730573282, 11320900.608852735, 11321083.982138049, 11321138.775027994, 11321326.597959816, 11321679.424123392, 11321871.487594794, 11321880.860088378, 11321933.110629348, 11322084.983885553, 11322146.38832061, 11322150.733748997, 11322204.6875, 11322498.282508334, 11322503.784577508, 11322565.477443617, 11322567.851905419, 11322635.44322009, 11322764.756154217, 11323065.309172725, 11323591.256887442, 11323717.1875, 11323784.315162808, 11323842.232603664, 11324790.25787179, 11324912.5, 11325099.600231985, 11325108.947533477, 11325198.191597437, 11325373.4375, 11325510.223137869, 11325921.875, 11326120.514350247, 11326199.188358214, 11326356.25, 11326502.16888725, 11326563.641463842, 11326720.750321187, 11326844.780205414, 11326869.130916703, 11327060.347661238, 11327192.15255655, 11327250.387432713, 11327740.948140536, 11327852.374336302, 11327885.9375, 11327893.653375106, 11327936.22266685, 11328158.623090526, 11328416.383855976, 11329900.771297958, 11330163.080988491, 11331813.801299619, 11331863.34358504, 11332116.895992663, 11332218.75, 11332723.630529251, 11332747.550000314, 11332892.350011595, 11332959.497017317, 11333099.261118433, 11333450.431331258, 11333453.125, 11334217.281076934, 11335158.161430348, 11335610.710293462, 11336081.25, 11336962.133895578, 11337790.858682366, 11338104.6875, 11338134.959872566, 11338138.312320285, 11339244.549027823, 11339499.902185172, 11339782.712953677, 11339825.09242571, 11339981.25, 11340043.756974308, 11340237.35518677, 11341339.92281948, 11346467.0043743, 11346700.24524572, 11346721.717291264, 11346832.657236839, 11347754.6875, 11347891.779960668, 11349281.25, 11350529.760357166, 11351874.886813449, 11353292.356322328, 11358829.312220389], [7.2551750148893595, 119.84383793028702, 152.02153261831268, 19.808276643067572, 24.59289029283166, 75.44705128695676, 46.338376023028076, 15.010015849639794, 112.40892167580526, 103.13657596892365, 50.390367303494905, 18.675349424312092, 81.20371827475279, 65.7895125410652, 7.732674171024905, 34.590663695327244, 49.46531620331742, 6.532518491591088, 32.615170170776594, 277.3564422309868, 5.913112831998846, 11.32380294461667, 39.34793985805817, 66.2155880584084, 76.81471423241454, 38.065188071998584, 5.676513030506049, 29.76094965996174, 29.753206151926378, 11.203743746378201, 71.29376680013323, 44.32412303061886, 67.46423853372754, 60.60973557263923, 59.006462212079015, 31.362549711884537, 27.413182174151263, 55.77323803075808, 10.768551040237826, 13.97680004256283, 41.73985906466917, 23.710292347494633, 51.300035993910456, 37.507123316706796, 48.84942818801527, 17.473073411300025, 10.135148885218918, 11.024914003806693, 29.64407847211484, 5.791116009965499, 59.08793678344231, 9.068354007916032, 5.5760511720110095, 52.40690793900788, 62.617475702191236, 21.197213093694447, 53.37694701490977, 36.30469338997641, 10.34434817823673, 61.482063613761994, 19.532311981345952, 11.17829607640673, 28.461576141567573, 31.444104205013755, 24.191621095843363, 47.33730737762052, 7.32907334274893, 23.39136475102014, 14.556970002147189, 60.89516338060983, 11.68466485457646, 49.11482084706725, 14.07362199461672, 9.240460343601601, 11.030983253017087, 25.774176741193806, 28.852850342375653, 63.67769920769554, 46.991416162662645, 20.62470980609154, 19.60136020340172, 16.0449717292396, 7.348416712332136, 30.71572821789934, 35.3094471586983, 16.47295404779002, 31.10215081601645, 94.07277693719828, 40.681733736930596, 52.99265209715478, 15.889549802406743, 39.26115242454695, 22.407837426703807, 26.801561126573198, 23.700781614822958, 19.525853086653193, 40.83170295956382, 87.34567969664721, 20.758636079561356, 17.665059866615703, 106.5287325853642, 6.967192788291198, 69.05740554490093, 13.24358516875305, 7.306893724973408, 26.476468120396298, 102.83142283729399, 49.23885817331956, 106.9361180896004, 21.92483444268057, 8.998487611282384, 61.10910517115368, 87.94384146554613, 71.43822733051148, 64.46489362784186, 88.13874858531274, 9.772767843772336, 7.965927848615289, 82.35413302391424, 56.31000655868047, 23.546036465511783, 13.65522569575368, 17.978265394948675, 14.240967628299853, 57.58020766290863, 22.890574934969663, 36.581951273712434, 10.443088003732775, 30.46546975797324, 16.340322314119383, 20.337500878021267, 16.440020075645418, 48.27188570409356, 5.225850714693417, 18.19212282761256, 21.13768324069304, 12.195052264775367, 78.55501320985148, 12.612168809081307, 62.39134498822649, 5.073377977682033, 29.19176356431536, 95.60283427197572, 17.179242174737347, 56.684480917974426, 93.73591249677267, 60.29913467728069, 39.65856475307234, 67.95927301137624, 79.16173016718072, 27.44467852157828, 46.096106101218176, 15.466094554430779, 13.415465531118151, 52.58182594937276, 31.218346706927125, 34.872521257131105, 9.217544436979303, 18.040457987006494, 76.67742221449961, 21.1020546026178, 60.57571209033055, 6.220211789558939, 8.448026798879964, 25.081485800102083, 48.75010732841402, 61.07373678633748, 6.771506200059151, 7.187074053548297, 49.64300134651855, 72.2239102053945, 34.22013662407855, 51.4982727928737, 45.35508767258757, 9.310288516846636, 5.362938467121785, 49.63054499825393, 38.273206736868644, 68.62634572841412, 106.42520308636469, 13.359964691940469, 10.081659345754693, 25.38723218198098, 14.81387268275413, 9.773925351541466, 158.56677975490172, 31.023392791961015, 5.962362818948428, 6.585615383199932, 14.766166506176322, 12.919333212818835, 7.8539843355734, 18.440293109212696, 24.48619354534041, 106.32836849468333, 78.51954787606529, 6.839277112307645, 5.6913594651509465, 95.32132024557711, 91.49995466502558, 46.69425747577368, 84.16756655876631, 51.07290128967238, 43.00975297043062, 53.29026245175248, 7.181275461351208, 54.93063701619354, 49.36148856187985, 8.548379858986479, 24.185499512661647, 43.094068821741814, 60.36292559246701, 42.71850066793812, 19.677826608111953, 15.21877337740057, 23.6392546254641, 32.8542468698288, 32.91196058839547, 163.70939304199646, 9.430692259681553, 26.905444956890424, 55.11795460534596, 25.909424152652996, 38.05071895487867, 84.11182712928706, 28.46226264904394, 25.2965681138856, 19.02536901466238, 24.775383421298994, 10.628413133526731, 37.86470464906046, 6.918152509144184, 21.175725025632328, 9.880498084302417, 25.335264902750144, 15.472917318338405, 76.62611054746354, 18.305139742686595, 10.906716425885788, 48.29591315864786, 32.92226940591586, 45.39191067241269, 13.116191400219549, 57.545044459428716, 6.614900554404024, 13.633592411358546, 8.60430684218417, 106.50570169238556, 6.747226790082193, 114.90516636494061, 6.585563989063977, 46.61321218395169, 80.40342349446172, 78.31252445233294, 118.48910607478383, 17.552962032910003, 16.609109604169326, 24.26649104001345, 27.242095360639134, 69.5306202430092, 59.982152272939956, 13.001035750187206, 9.610149749401783, 54.49936463506427, 24.146042700358212, 51.59479189230557, 19.548502452657754, 5.675183990048977, 19.79555978680541, 30.92082965212629, 119.68803313406201, 45.50634495702023, 104.98674740063788, 6.3019555762509185, 18.312610111152228, 38.436439219004384, 34.40920683110593, 66.97578722502051, 13.008605977643963, 9.17924605938643, 39.82781182655966, 90.40151056564521, 15.567454799194547, 54.88851466835766, 109.778967357964, 45.55228653536312, 103.16902728388949, 15.823100894304012, 92.04889486393427, 6.005213659493739, 15.324271188839777, 59.63360253971818, 120.18920007265005, 82.06780795442725, 116.05354956530776, 6.728448298669941, 68.14266059235064, 19.341764251702333, 31.139112584737425, 15.866270053194649, 50.19979684355037, 26.323544629419544, 13.54367226133342, 42.8479630990655, 205.59059661332435, 25.20377929278292, 7.616680189020896, 23.671797779643743, 11.212338831504209, 8.47687944455657, 15.104814063464582, 34.52565074950944, 24.844221247988646, 93.21532413212995, 127.01382525345319, 40.436345571186294, 13.746745261358983, 20.314412541689293, 43.03359713194706, 100.8655735489136, 7.930469845882459, 101.66725511027254, 11.109656462647274, 5.934894726301829, 25.26918680187261, 31.01138930364673, 68.63025337769808, 7.832994023771672, 32.31559942097935, 39.5237027350131, 7.746857558561615, 16.114141317665776, 74.71932300107517, 12.333475791618511, 88.13191934343703, 19.44072166197472, 27.060761835328442, 98.12964802837895, 84.22530683267914, 7.734332076229037, 26.476535331596235, 7.831365593886226, 39.28223611726116, 12.704239603460275, 29.872406821388168, 24.968655669004526, 57.30471456331142, 17.034767237429357, 35.57120838143452, 12.859418060690793, 92.95441294732346, 26.700818855772003, 13.862069295263328, 5.233275081900455, 18.47783280835758, 70.24979391166111, 38.80160729123631, 55.07211710229137, 7.561893500837759, 18.641812509219353, 13.037861357397487, 30.18266787559698, 9.779567925906827, 12.76285912570457, 16.847813704553797, 64.63228231756085, 61.65047984476783, 6.954333105957183, 70.05161198863958, 23.301149835319805, 33.54194190420254, 5.59208510843894, 50.49665250263954, 5.662593596654664, 71.83215771744301, 17.45524505327815, 11.663022673407148, 69.61689931094872, 22.301781015730366, 7.100564398430099, 7.018948502166853, 40.96460613358928, 50.467720131430426, 21.404886892057064, 5.179471281368416, 13.182793956969649, 61.11997083923525, 37.880067285766714, 13.741604199185788, 7.3346630888231035, 56.54344289210581, 44.64074542939467, 61.15431016779998, 13.404952074256649, 37.088881956617215, 73.2336476943661, 45.60204712095868, 57.01947359508412, 18.823896215807075, 11.09802174933678, 23.743331274510318, 52.46389288836986, 12.348351290027548, 17.832800658052665, 5.692566550993013, 13.126749616478339, 5.766375814361699, 8.122439458269653, 11.225801816480574, 62.24644519952106, 15.445994568772756, 27.286855151902678, 48.24307000688432, 9.86168477927395, 10.24532938646214, 55.76439556790817, 78.22219250663505, 104.25951066294515, 25.258589423447926, 5.305686329531366, 33.7910799404562, 5.0476966913461165, 42.95029267443668, 40.72128290271927, 9.875271299938847, 7.4616696319151, 7.050323616381267, 6.7813681844031235, 117.00112311533869, 6.794477002095908, 93.10035706667966, 27.36726988931169, 18.236984138372307, 63.60044383695316, 6.120672321000391, 8.051181268988454, 25.43896272996389, 7.389125290979001, 43.9352636719009, 9.328216887246377, 118.19174877250668, 74.38679937171128, 35.07545161993774, 60.28757185266415, 15.554755528282385, 9.204473200443676, 76.95709074491481, 40.39158235913895, 58.760368808331336, 27.607424208482566, 62.62766159425012, 33.68121683336272, 43.479220042337985, 6.9781528693209145, 81.1941523999904, 68.21588594414116, 31.177081243559574, 14.078250094386643, 37.08228531952298, 62.878279417492244, 31.87910102294951, 5.397790162435022, 28.878874303665455, 63.546361626100214, 82.05195259778432, 46.49932285379889, 17.009644338127693, 65.30168994812132, 143.56890410675777, 35.521932838890514, 70.17353343848521, 52.87154829530011, 8.884418660666768, 7.453864284746731, 8.685276689138812, 102.35087113921301, 7.7542760816463145, 15.580175986879201, 50.74846257921592, 41.99699763321069, 11.709172070754896, 81.38331835406731, 28.057981243189598, 30.02283244473316, 17.780326468529026, 157.43515106932023, 38.59133467379219, 119.26915340769791, 40.9004231073895, 52.409174656482755, 9.578037414195288, 12.080942338989598, 57.65093914669735, 81.16590505104288, 85.40007875415573, 22.881114349167007, 26.930949681629084, 62.13686803014808, 9.787481818188544, 13.883826434670878, 61.81420022948039, 18.51386304072372, 7.433182218365194, 10.700327403959104, 10.412650918565499, 9.849128159384023, 50.269543890817964, 56.529775446271586, 95.87629464684606, 99.01673716032785, 37.60479313757054, 27.40149000198346, 23.981240040930775, 9.041718266848186, 26.805518785993684, 11.368927592830385, 45.96694086148856, 46.75996384262662, 49.11760854457064, 12.439751559786743, 40.27051856568281, 14.019792185382276, 44.00598180970635, 13.470362104809503, 32.72146062626551, 49.21827258130136, 12.844982742989215, 12.152525041876098, 82.64709255835143, 26.012051334013425, 20.766591601853285, 37.21053307604652, 15.515188493808616, 31.096948631745988, 8.520862582224797, 90.41354314205464, 52.70205513636013, 6.902002171450184, 42.20099822988445, 10.012776914026603, 34.35019136764661, 12.415902662002086, 46.006166599420276, 9.792568403398839, 132.43080375430898, 6.76509465758733, 19.63665898119953, 13.535669864585604, 16.90707177595369, 7.12910338273986, 87.24127014257014, 46.91265534611848, 9.145671677697651, 42.45596967658501, 5.66892394263845, 6.628968396886373, 110.46249707851658, 38.611911341949565, 20.729920802614416, 19.249692766052565, 24.618567885093498, 14.230080156891988, 19.07842113173955, 10.289889557600318, 44.1304769620009, 7.339418970411206, 42.43491677392248, 22.41933310063579, 72.34338140465739, 7.5283228334533, 25.84790064286748, 6.381079578615405, 111.43637184664681, 7.225134688878932, 6.954666379161702, 96.86224645499112, 30.785141655415682, 8.474971762122934, 11.739584380918211, 55.056540076223214, 5.314992444920363, 18.48708417838469, 75.6664027935307, 12.05258108684769, 10.656044497309647, 10.404997800876043, 56.05349345075088, 67.96456817486609, 81.22657688654961, 49.42016733799891, 12.296689268872404, 26.161879412843955, 45.38508459961883, 20.427432493791112, 22.517226964177972, 21.58333624988548, 11.139142289263425, 59.060892753653555, 6.147958673166586, 8.496137062902688, 90.30390066336565, 39.600732521789794, 67.31049055069118, 124.35738100255621, 21.400152285708067, 19.64085404725156, 121.12299237559941, 68.20313316995285, 46.19435799477614, 15.134919452917416, 19.886781820820794, 27.162238637196495, 14.214072938293056, 115.14906277809592, 7.045879585977212, 142.95353742968172, 31.60500151064182, 16.67321112133154, 11.45311682700216, 8.222604877953065, 16.61028887905804, 21.49238866585931, 10.323817461630282, 59.707043868910574, 106.6933307669526, 42.45785193152026, 10.25874787538329, 23.753515835725985, 15.842270834389836, 54.247421554061184, 11.718015915453655, 9.251826900667709, 55.555421571855824, 5.182798693038229, 11.287845030306674, 21.609064889339955, 28.68418149825625, 58.57646450613494, 14.262417005490608, 16.612572203656686, 16.236087052979443, 7.07012245466578, 36.74755115014632, 5.56955315814885, 38.91744749703327, 78.85738379606225, 66.55885619672107, 39.64570920071233, 59.29291459700674, 43.05116230352876, 57.616583933068526, 13.018545934843285, 38.0480565816235, 24.925011809751684, 48.79848066987436, 48.274437799688435, 19.76507371158679, 8.565946758021937, 59.157588927921495, 57.14227355593507, 9.275322842239122, 14.083832907023066, 14.365127431573358, 16.940142017461092, 33.65000424197482, 22.38482569043266, 9.8423764229321, 103.30964749356153, 93.482737580894, 70.28695282389145, 36.95272628687471, 13.763226292274851, 66.32050179751066, 105.52022887233039, 20.58372196333329, 57.306558616504354, 7.2055606193971675, 6.901695760628813, 29.114219935320538, 5.196310986913392, 64.93439841265513, 66.54164242492251, 8.601749250514457, 16.577497833620654, 21.034827730068592, 71.15904863888585, 71.14756778662918, 83.33258945489067, 40.09744323758643, 10.758686591572488, 32.71729187462118, 72.60928034296013, 23.182448164356757, 20.02688139994615, 35.64640238564675, 19.619730168894044, 10.387210959701576, 108.49878733131999, 75.39827125131204, 8.518593311533348, 50.09760170113625, 5.783556515088281, 7.7226393825025, 17.362973213161744, 80.6800037352902, 45.02780319648969, 6.515288961328326, 21.02737688071963, 54.36451013450317, 57.52928577866641, 24.29359490456843, 49.28102430125691, 35.8859259513479, 8.56080309083324, 12.854971881557413, 8.88725896384838, 21.0800477677967, 98.2668911820237, 14.344215517868456, 20.26298700390288, 7.314621461467949, 54.44047864778453, 45.24597043884262, 24.081153968623944, 11.605888653834558, 81.92891659436373, 11.077075817245964, 72.47645741895147, 14.150248660347113, 38.98791894640125, 8.310910804301987, 14.562595175466456, 113.34002985790485, 45.54586004938276, 8.255063789030922, 15.748738231301685, 12.334676676031608, 37.13449911745443, 25.16227956950337, 76.85488255560675, 65.08240181335337, 13.072862840680108, 6.511161887246643, 36.836853509896514, 32.28105111670073, 9.011200476779752, 67.99261973751265, 43.56414154407412, 22.664659544292277, 10.641535346023423, 34.64376026085962, 34.85294630448037, 13.484231145440194, 28.36688965478537, 58.00923280980163, 35.763347543902626, 10.676355788751687, 21.408481350254558, 6.935108643188424, 11.81394661943359, 22.301879646778744, 5.798056980720284, 94.08282554710377, 22.67539113883788, 22.799743419864125, 11.478324417670079, 53.32099708503726, 30.313557087185988, 19.386367771690736, 18.844755698317343, 44.82455210759767, 15.373962501989634, 9.515152632133866, 16.939112837300648, 28.396269204357306, 10.335418544876822, 89.17929630357061, 16.840582885523023, 14.764919650394848, 19.293637630127673, 10.818970722833296, 17.912239148805387, 22.802790503154316, 23.11994332941388, 24.98230979382431, 24.23010562508092, 10.571739601893785, 13.638015383570616, 8.450786848419284, 5.422272727496316, 30.803919297691557, 6.16787773870539, 58.70960199565301, 21.757659586333236, 6.354290042601948, 79.01645583162893, 33.34361337283414, 10.43402929171163, 18.463413928011903, 51.398250308694465, 14.286148944932568, 39.54018858552626, 14.718239656426421, 7.133358728257916, 9.039139006726538, 101.67541212273028, 22.554844581646957, 81.56092029014815, 98.24038281892368, 91.74922363819401, 7.3638018403555074, 49.324132034224064, 52.30603773350151, 12.546541696075952, 55.03255187031403, 5.110580889802128, 13.547640273474753, 5.978462608040431, 140.66736056961003, 6.8420134193939175, 10.315950757587592, 47.18125856105616, 5.635843547192914, 6.538968180383621, 21.88992975546214, 28.648955299399073, 26.749893324709912, 69.50432749701974, 51.68692180174859, 12.564852483727488, 10.081517978439976, 84.29318403235067, 14.031376121318083, 8.789471837751982, 6.70667215533237, 44.29845720927677, 7.67457766137185, 80.51883119795062, 5.404723250860584, 22.0490516942453, 12.716189178806728, 33.87083847411847, 28.595637073850124, 5.167546674236084, 8.480139111408706, 24.72735125141671, 44.190865180027814, 15.528885568827413, 11.513728945059556, 34.379631842273476, 17.764093391162096, 110.36408081759643, 10.112704532840677, 5.6086556117359905, 14.575191366941128, 9.410796500229058, 45.32471755036951, 5.245340531125099, 6.464627593235395, 41.34674346917381, 25.368012276144974, 26.56679033128959, 78.43609293008124, 101.68991413483468, 88.73750329161878, 19.75121208094216, 77.54669677668076, 21.242898363177446, 26.471433912750108, 51.41048584252594, 88.97567150220965])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([6743402.047899941, 6798381.602915615, 6830177.945460082, 7077994.625212399, 7158199.763187827, 7282689.0625, 7290009.911801992, 7292031.267858106, 7311199.491457405, 7322373.542119066, 7325259.086153767, 7592617.9120263895, 7619414.0625, 7650001.252156334, 7818962.844233887, 7892179.3971398175, 7993259.776386856, 7993287.783796568, 8061903.125, 8099346.142670247, 8147526.255549074, 8191301.684435399, 8198802.344454739, 8243929.6875, 8371743.8449094165, 8559257.8125, 8565608.889043516, 8565949.27115848, 8568086.754444497, 8666300.048107455, 8673133.781371502, 8673821.875, 8675409.375, 8675463.81485712, 8675915.625, 8677207.594032718, 8688659.316214968, 8696385.935082313, 8795041.94489316, 8806132.534069385, 8819115.625, 8821127.23927495, 8828546.705870736, 8829089.592235776, 8829220.676252192, 8829486.473690713, 8829860.193799892, 8830262.433225669, 8831203.125, 8831284.06113106, 8831308.088159239, 8833877.79661033, 8833913.354696944, 8837296.30220899, 8862834.578185035, 8875394.782318195, 8877307.709458087, 8883792.59997623, 8885794.745985812, 8887077.432332234, 8888476.144987103, 8889162.532126136, 8889599.69769525, 8890515.625, 8915554.421248462, 8916862.5, 8917558.547532491, 8919004.91653882, 8920257.707815751, 8921026.523484185, 8921047.204317106, 8921781.25, 8921909.127004456, 8922135.82795166, 8922645.249725224, 8922682.223801712, 8922697.036171935, 8922722.321188247, 8922774.338740181, 8922932.971570926, 8923010.224360788, 8923025.262222094, 8923128.92335266, 8923187.21360701, 8923301.38280371, 8923522.226092564, 8923804.6875, 8924259.684525425, 8924300.0, 8924482.8125, 8924681.478319036, 8924844.117997319, 8924882.814891769, 8924952.352818947, 8925832.459919259, 8926000.839751719, 8928321.875, 8931793.671313053, 8931845.1362293, 9005839.434806848, 9010785.613939468, 9037091.566264534, 9079485.978516554, 9113899.469612477, 9129703.78641302, 9139654.0289821, 9141640.625, 9150921.875, 9155134.217347628, 9163389.32856721, 9173979.265661016, 9227791.578400794, 9229456.25, 9231723.4375, 9233507.8125, 9236501.5625, 9237008.292103166, 9238279.027579477, 9262815.625, 9275746.875, 9278123.859842023, 9290193.365947057, 9290581.21718795, 9293156.286744399, 9321256.25, 9326245.896215344, 9361627.503131231, 9363169.825398516, 9393771.668867664, 9394259.14901842, 9394879.326302538, 9395200.344477037, 9395526.145135935, 9395536.79565682, 9395728.77775762, 9396083.838617431, 9396297.465110872, 9397680.187797802, 9399192.966858855, 9442300.0, 9448465.869956706, 9449515.625, 9449869.31972607, 9450339.547300255, 9450360.931962002, 9451599.106036844, 9451673.152522145, 9451799.269565966, 9451939.555731645, 9452510.216497613, 9452637.258051476, 9452704.6875, 9453697.93623441, 9453769.84482801, 9455366.549606072, 9455778.125, 9455916.153884234, 9456256.112203613, 9456472.50418881, 9456614.028758012, 9457790.084885208, 9458837.79887042, 9471321.49382651, 9486562.793315124, 9487148.133209087, 9487582.8125, 9488648.308595479, 9489118.789938202, 9489733.056391979, 9489857.46482766, 9489869.261599962, 9490036.282599531, 9491150.0, 9491831.25, 9493293.724922678, 9493450.600779254, 9494657.512572423, 9494882.408332054, 9495812.5, 9497845.150561508, 9498379.216677101, 9511203.896424336, 9512483.749137603, 9512506.216105068, 9512577.043258594, 9516327.24100952, 9517293.75, 9520831.509894272, 9520838.100977363, 9520887.285703335, 9528815.049227403, 9541395.956358386, 9547688.11173419, 9549129.44297439, 9549323.074040601, 9549505.7887901, 9549530.055542015, 9549634.130297204, 9550012.5, 9550361.336093327, 9551120.3125, 9551311.873281494, 9551546.875, 9552032.8125, 9552039.718232624, 9552156.472794136, 9589113.481928477, 9595620.3125, 9607625.464300469, 9632885.67079013, 9667444.472582666, 9670919.70969289, 9673465.504699042, 9675728.38105047, 9676754.758390075, 9677471.566221675, 9677641.649909917, 9677748.295663418, 9677812.179909376, 9678812.615135925, 9679031.162448786, 9679152.133787747, 9679246.205583414, 9680095.3125, 9680585.689703846, 9681243.220494457, 9688088.874592062, 9690758.949480427, 9693903.621325847, 9697162.155219382, 9697215.625, 9697216.35321655, 9738272.905085005, 9755259.960072728, 9761356.064686337, 9785026.695268262, 9813548.4375, 9841014.732487738, 9844292.041361872, 9844482.576824209, 9845581.38642707, 9847800.0, 9848280.67154471, 9848316.891865669, 9850382.019590857, 9851938.06102836, 9857828.681101697, 9861078.021794207, 9864912.417904591, 9865471.800352743, 9865520.913177595, 9867945.564919032, 9868532.8125, 9869421.274775546, 9870286.863985201, 9870940.969802054, 9871790.262409143, 9871904.307981879, 9872725.247330341, 9874207.354109738, 9874225.0, 9874258.334242392, 9874870.624736222, 9875300.0, 9876315.355329145, 9876456.25, 9877753.965418749, 9879672.929685753, 9881135.248643547, 9885026.649101434, 9887716.91539338, 9905766.364989016, 9918983.555378642, 9919048.707737828, 9921847.546116821, 9927783.084713886, 9929404.6875, 9944340.625, 9944347.817351045, 9952073.251897117, 9962678.125, 9963928.125, 9969267.342571348, 9970889.888190618, 9972737.5, 9973657.228670252, 9975851.5625, 9976811.151044188, 9979032.8125, 9980025.725869425, 9980478.221923847, 9982718.59516545, 9984211.652536118, 9997553.125, 9997617.1875, 9999187.321820846, 10001132.8125, 10001358.386192853, 10001573.16223938, 10001589.285662187, 10002332.8125, 10002433.307225011, 10035748.874088014, 10067467.887844063, 10081029.325270029, 10081437.864537299, 10082385.97451424, 10084497.20998952, 10092770.365978505, 10197839.40582428, 10282825.04219514, 10421987.5, 10551493.531040112, 10559381.248354742, 10645693.75, 10645712.483247329, 10646298.084241884, 10654005.623125719, 10663141.278514914, 10663745.219320936, 10680872.697373778, 10683051.563900454, 10683396.863421177, 10685053.073056784, 10716763.601832153, 10720984.375, 10721437.237904074, 10721481.84350176, 10722329.6875, 10722676.188178983, 10723466.005365878, 10730434.652022725, 10747680.888836363, 10747683.467603015, 10752866.650712062, 10754037.523062088, 10760308.497057369, 10779490.625, 11016616.564249625, 11017715.973472947, 11025346.652957622, 11025805.277069805, 11034679.6875, 11046143.006089957, 11048287.491699563, 11051360.727217836, 11051699.564975418, 11053482.595030405, 11054983.59835332, 11055007.132792914, 11062893.897731729, 11064693.102195082, 11067666.872203907, 11068452.695268223, 11069261.552363615, 11069541.43331458, 11070129.6875, 11078495.3125, 11079040.509224212, 11079910.094702065, 11079941.578790737, 11081456.25, 11081509.861383995, 11084193.563882072, 11084314.48643726, 11088002.530185007, 11088805.33647758, 11094739.948669706, 11096678.125, 11097071.863696026, 11097179.65920913, 11099583.971130332, 11100087.5, 11101771.157947425, 11103187.00023111, 11107595.845638173, 11109972.192455528, 11111576.5625, 11112555.827417359, 11114137.419940941, 11114871.470160896, 11117139.0625, 11117921.875, 11118027.616111327, 11118238.095704166, 11119807.259347642, 11120008.678784892, 11120142.211577652, 11123608.217645627, 11221038.253414359, 11224935.422531042, 11227424.675069243, 11233421.535170479, 11243921.40356013, 11244980.222639656, 11245903.556423357, 11246008.37408917, 11246260.751260443, 11246783.880057987, 11247253.105925048, 11247388.761919646, 11248786.09665295, 11248786.599453673, 11248953.721059587, 11249387.27468009, 11249557.202176694, 11249780.585204218, 11249830.926156651, 11249923.499364562, 11250673.4375, 11250926.688847262, 11251054.945157798, 11251100.0, 11251136.504450615, 11251288.413683478, 11251642.1875, 11251689.234674538, 11251873.4375, 11251969.15617041, 11252055.153007187, 11252314.302013319, 11252725.361950723, 11252806.25, 11253150.0, 11253470.892550323, 11253715.215855062, 11253962.222325455, 11254337.588814521, 11254529.6875, 11254547.028175132, 11254753.125, 11254809.313665938, 11254824.89800638, 11255043.755572073, 11255087.807511289, 11255312.026655622, 11255341.25025888, 11255354.053990059, 11255478.125, 11255943.746229323, 11257713.00638248, 11258909.058294, 11259077.645152293, 11259133.89665048, 11259226.714597218, 11259404.53347181, 11259458.467007907, 11259509.375, 11259540.625, 11259549.16405323, 11259929.6875, 11260004.6875, 11260180.027693259, 11260332.237797504, 11260428.125, 11260445.3125, 11260452.04862387, 11260523.190944107, 11260563.988868613, 11260667.1875, 11260690.625, 11260942.169599727, 11260947.130013453, 11260955.401367022, 11260993.75, 11261171.432129772, 11262102.016553279, 11262173.05804178, 11262306.957485372, 11262389.0625, 11262413.235523844, 11262451.5625, 11262582.231918441, 11262621.442216007, 11262706.93482574, 11262717.098555, 11262783.211645022, 11262902.169929666, 11262916.03506481, 11262981.25, 11263214.405164868, 11263391.079617755, 11263479.168062484, 11263535.32757945, 11263547.486920984, 11263643.007283958, 11263673.236488672, 11263673.4375, 11263720.3125, 11263926.5625, 11263981.527958298, 11263986.472287355, 11264203.125, 11264241.716405364, 11264510.03419306, 11264737.933939505, 11265021.231589755, 11265325.0, 11265350.49425418, 11265734.066763027, 11265845.120571652, 11266047.576086733, 11266104.191502687, 11266123.145713773, 11266205.738694292, 11266508.367286494, 11266551.5625, 11266688.976482674, 11266749.200895464, 11266934.998156419, 11267168.75, 11267518.265620317, 11267670.030333158, 11267837.1363112, 11267954.594134295, 11268284.524032107, 11268362.5, 11268856.25, 11269122.941708872, 11269447.578494463, 11269635.644399654, 11270027.84429153, 11270029.972819598, 11270209.738215491, 11270540.625, 11270912.5, 11270936.545741724, 11270953.694681376, 11271130.923449667, 11271573.296061149, 11271637.176577523, 11271660.9375, 11271697.017787416, 11271758.690429682, 11272092.605731577, 11272304.79168194, 11272404.6875, 11272722.531271096, 11272748.4375, 11272758.521600097, 11272879.540220736, 11273017.880245928, 11273539.0625, 11273747.37430062, 11273860.686995288, 11274009.382949473, 11274219.567241652, 11274394.887295455, 11275032.825021883, 11275064.091528695, 11275096.875, 11275519.695452455, 11275656.021861097, 11275784.375, 11275843.138644956, 11275849.286556767, 11276034.042185161, 11276076.528115215, 11276458.79493328, 11276852.366278702, 11277007.179686116, 11277195.112426566, 11277610.1557915, 11277665.013559746, 11277814.764031444, 11278172.644000186, 11278181.25, 11278524.588094536, 11279001.5625, 11279029.475510178, 11279175.348068882, 11279477.413664622, 11280029.6875, 11280276.942770876, 11280307.13760438, 11280514.0625, 11280612.905062014, 11281131.590039296, 11281212.859105937, 11281338.0254681, 11281939.817494188, 11282253.822148083, 11282336.208407873, 11282506.770783184, 11282890.300723907, 11283038.545070168, 11283476.5625, 11283593.75, 11283860.537690671, 11284084.685542293, 11284192.825551659, 11284390.492906354, 11284569.620365875, 11284636.886096908, 11284692.306522377, 11284793.258998405, 11285271.88230377, 11285352.090823002, 11285438.05660037, 11285508.18881671, 11285784.375, 11285937.766805628, 11286274.963756703, 11286323.51819735, 11286358.865631655, 11286564.916849026, 11286599.70387126, 11286833.56318754, 11287138.385939635, 11287366.314209733, 11287371.19734061, 11287404.059194844, 11287863.857593374, 11287970.46836263, 11288402.829492657, 11288639.261809329, 11289400.696880508, 11289941.01555833, 11290562.698667949, 11290684.252035806, 11290841.003296647, 11290843.265150448, 11290893.16664481, 11290932.8125, 11291570.277706416, 11291614.0625, 11292246.237497848, 11292673.137717148, 11292836.15532412, 11292890.625, 11293444.09745987, 11294413.966865301, 11294434.375, 11294505.420791708, 11294600.378517732, 11295079.199046861, 11295170.3125, 11295379.073242508, 11298315.434642745, 11298508.191243798, 11298673.740197107, 11298751.363568773, 11298751.386331443, 11298869.479272414, 11298884.2235384, 11299041.57737906, 11299295.3125, 11299351.225879807, 11299458.993460516, 11299513.968991028, 11299763.922689699, 11300493.1337465, 11300568.43669597, 11300592.110096894, 11300617.023694413, 11300625.0, 11300800.892578958, 11300831.743397139, 11300859.013451736, 11301268.75, 11301432.745423332, 11301527.880178917, 11301624.840088332, 11301831.69121393, 11301915.625, 11301924.584974673, 11302039.610722868, 11302134.375, 11302412.443471853, 11302415.25509626, 11302803.125, 11302905.902223945, 11302983.247437404, 11303037.5, 11303161.182508612, 11303168.75, 11303195.623084964, 11303306.904897362, 11303314.342610575, 11303530.968579883, 11304043.75, 11304339.98514398, 11304681.801956072, 11304805.05935668, 11304808.827799581, 11305057.8125, 11305139.0625, 11305198.77722968, 11305401.5625, 11305468.500346933, 11305502.99305331, 11305945.987687774, 11305960.5916688, 11306022.521296643, 11306342.1875, 11306514.9141842, 11307116.655777067, 11307154.6875, 11307454.352491166, 11307526.97690535, 11307893.195201531, 11308052.705176538, 11308634.69348757, 11308806.767206019, 11308979.6875, 11309031.25, 11309154.865051925, 11309256.538312335, 11309265.625, 11309285.729680665, 11309630.99134203, 11309749.769177794, 11309964.0625, 11310441.109902557, 11310474.31848592, 11310761.92346772, 11311055.70757198, 11311094.08640303, 11311220.146440255, 11311239.949898243, 11311400.687187698, 11311829.626127277, 11312009.156366255, 11312634.53075563, 11312836.223302364, 11313051.087343905, 11313079.654864563, 11313247.215488696, 11313331.025665622, 11313492.1875, 11313854.733233238, 11314059.27190385, 11314092.1875, 11314126.482705917, 11314179.627428213, 11314233.009305712, 11314252.218157006, 11314373.4375, 11314429.490514433, 11314443.75, 11314546.875, 11314596.279876657, 11314676.714352708, 11314685.053204091, 11314808.761223745, 11314813.657872908, 11315124.739334922, 11315214.0625, 11315249.6099362, 11315284.039979445, 11315345.691438667, 11315346.875, 11315490.181759084, 11315657.804461807, 11315731.088141968, 11315746.875, 11315757.503182415, 11315885.428975312, 11316336.740289332, 11316342.58624331, 11316685.557988869, 11316762.282451184, 11316893.805212267, 11316953.606789399, 11317383.996715065, 11317508.025802061, 11317564.327525081, 11317939.0625, 11318050.499148585, 11318467.936054807, 11318557.37125582, 11318573.49468094, 11318650.209601376, 11319079.602021655, 11319407.87428653, 11319507.238675147, 11319531.084129144, 11320023.86007344, 11320290.59555386, 11320511.747091891, 11320885.730573282, 11320900.608852735, 11321083.982138049, 11321138.775027994, 11321326.597959816, 11321679.424123392, 11321871.487594794, 11321880.860088378, 11321933.110629348, 11322084.983885553, 11322146.38832061, 11322150.733748997, 11322204.6875, 11322498.282508334, 11322503.784577508, 11322565.477443617, 11322567.851905419, 11322635.44322009, 11322764.756154217, 11323065.309172725, 11323591.256887442, 11323717.1875, 11323784.315162808, 11323842.232603664, 11324790.25787179, 11324912.5, 11325099.600231985, 11325108.947533477, 11325198.191597437, 11325373.4375, 11325510.223137869, 11325921.875, 11326120.514350247, 11326199.188358214, 11326356.25, 11326502.16888725, 11326563.641463842, 11326720.750321187, 11326844.780205414, 11326869.130916703, 11327060.347661238, 11327192.15255655, 11327250.387432713, 11327740.948140536, 11327852.374336302, 11327885.9375, 11327893.653375106, 11327936.22266685, 11328158.623090526, 11328416.383855976, 11329900.771297958, 11330163.080988491, 11331813.801299619, 11331863.34358504, 11332116.895992663, 11332218.75, 11332723.630529251, 11332747.550000314, 11332892.350011595, 11332959.497017317, 11333099.261118433, 11333450.431331258, 11333453.125, 11334217.281076934, 11335158.161430348, 11335610.710293462, 11336081.25, 11336962.133895578, 11337790.858682366, 11338104.6875, 11338134.959872566, 11338138.312320285, 11339244.549027823, 11339499.902185172, 11339782.712953677, 11339825.09242571, 11339981.25, 11340043.756974308, 11340237.35518677, 11341339.92281948, 11346467.0043743, 11346700.24524572, 11346721.717291264, 11346832.657236839, 11347754.6875, 11347891.779960668, 11349281.25, 11350529.760357166, 11351874.886813449, 11353292.356322328, 11358829.312220389], [7.2551750148893595, 119.84383793028702, 152.02153261831268, 19.808276643067572, 24.59289029283166, 75.44705128695676, 46.338376023028076, 15.010015849639794, 112.40892167580526, 103.13657596892365, 50.390367303494905, 18.675349424312092, 81.20371827475279, 65.7895125410652, 7.732674171024905, 34.590663695327244, 49.46531620331742, 6.532518491591088, 32.615170170776594, 277.3564422309868, 5.913112831998846, 11.32380294461667, 39.34793985805817, 66.2155880584084, 76.81471423241454, 38.065188071998584, 5.676513030506049, 29.76094965996174, 29.753206151926378, 11.203743746378201, 71.29376680013323, 44.32412303061886, 67.46423853372754, 60.60973557263923, 59.006462212079015, 31.362549711884537, 27.413182174151263, 55.77323803075808, 10.768551040237826, 13.97680004256283, 41.73985906466917, 23.710292347494633, 51.300035993910456, 37.507123316706796, 48.84942818801527, 17.473073411300025, 10.135148885218918, 11.024914003806693, 29.64407847211484, 5.791116009965499, 59.08793678344231, 9.068354007916032, 5.5760511720110095, 52.40690793900788, 62.617475702191236, 21.197213093694447, 53.37694701490977, 36.30469338997641, 10.34434817823673, 61.482063613761994, 19.532311981345952, 11.17829607640673, 28.461576141567573, 31.444104205013755, 24.191621095843363, 47.33730737762052, 7.32907334274893, 23.39136475102014, 14.556970002147189, 60.89516338060983, 11.68466485457646, 49.11482084706725, 14.07362199461672, 9.240460343601601, 11.030983253017087, 25.774176741193806, 28.852850342375653, 63.67769920769554, 46.991416162662645, 20.62470980609154, 19.60136020340172, 16.0449717292396, 7.348416712332136, 30.71572821789934, 35.3094471586983, 16.47295404779002, 31.10215081601645, 94.07277693719828, 40.681733736930596, 52.99265209715478, 15.889549802406743, 39.26115242454695, 22.407837426703807, 26.801561126573198, 23.700781614822958, 19.525853086653193, 40.83170295956382, 87.34567969664721, 20.758636079561356, 17.665059866615703, 106.5287325853642, 6.967192788291198, 69.05740554490093, 13.24358516875305, 7.306893724973408, 26.476468120396298, 102.83142283729399, 49.23885817331956, 106.9361180896004, 21.92483444268057, 8.998487611282384, 61.10910517115368, 87.94384146554613, 71.43822733051148, 64.46489362784186, 88.13874858531274, 9.772767843772336, 7.965927848615289, 82.35413302391424, 56.31000655868047, 23.546036465511783, 13.65522569575368, 17.978265394948675, 14.240967628299853, 57.58020766290863, 22.890574934969663, 36.581951273712434, 10.443088003732775, 30.46546975797324, 16.340322314119383, 20.337500878021267, 16.440020075645418, 48.27188570409356, 5.225850714693417, 18.19212282761256, 21.13768324069304, 12.195052264775367, 78.55501320985148, 12.612168809081307, 62.39134498822649, 5.073377977682033, 29.19176356431536, 95.60283427197572, 17.179242174737347, 56.684480917974426, 93.73591249677267, 60.29913467728069, 39.65856475307234, 67.95927301137624, 79.16173016718072, 27.44467852157828, 46.096106101218176, 15.466094554430779, 13.415465531118151, 52.58182594937276, 31.218346706927125, 34.872521257131105, 9.217544436979303, 18.040457987006494, 76.67742221449961, 21.1020546026178, 60.57571209033055, 6.220211789558939, 8.448026798879964, 25.081485800102083, 48.75010732841402, 61.07373678633748, 6.771506200059151, 7.187074053548297, 49.64300134651855, 72.2239102053945, 34.22013662407855, 51.4982727928737, 45.35508767258757, 9.310288516846636, 5.362938467121785, 49.63054499825393, 38.273206736868644, 68.62634572841412, 106.42520308636469, 13.359964691940469, 10.081659345754693, 25.38723218198098, 14.81387268275413, 9.773925351541466, 158.56677975490172, 31.023392791961015, 5.962362818948428, 6.585615383199932, 14.766166506176322, 12.919333212818835, 7.8539843355734, 18.440293109212696, 24.48619354534041, 106.32836849468333, 78.51954787606529, 6.839277112307645, 5.6913594651509465, 95.32132024557711, 91.49995466502558, 46.69425747577368, 84.16756655876631, 51.07290128967238, 43.00975297043062, 53.29026245175248, 7.181275461351208, 54.93063701619354, 49.36148856187985, 8.548379858986479, 24.185499512661647, 43.094068821741814, 60.36292559246701, 42.71850066793812, 19.677826608111953, 15.21877337740057, 23.6392546254641, 32.8542468698288, 32.91196058839547, 163.70939304199646, 9.430692259681553, 26.905444956890424, 55.11795460534596, 25.909424152652996, 38.05071895487867, 84.11182712928706, 28.46226264904394, 25.2965681138856, 19.02536901466238, 24.775383421298994, 10.628413133526731, 37.86470464906046, 6.918152509144184, 21.175725025632328, 9.880498084302417, 25.335264902750144, 15.472917318338405, 76.62611054746354, 18.305139742686595, 10.906716425885788, 48.29591315864786, 32.92226940591586, 45.39191067241269, 13.116191400219549, 57.545044459428716, 6.614900554404024, 13.633592411358546, 8.60430684218417, 106.50570169238556, 6.747226790082193, 114.90516636494061, 6.585563989063977, 46.61321218395169, 80.40342349446172, 78.31252445233294, 118.48910607478383, 17.552962032910003, 16.609109604169326, 24.26649104001345, 27.242095360639134, 69.5306202430092, 59.982152272939956, 13.001035750187206, 9.610149749401783, 54.49936463506427, 24.146042700358212, 51.59479189230557, 19.548502452657754, 5.675183990048977, 19.79555978680541, 30.92082965212629, 119.68803313406201, 45.50634495702023, 104.98674740063788, 6.3019555762509185, 18.312610111152228, 38.436439219004384, 34.40920683110593, 66.97578722502051, 13.008605977643963, 9.17924605938643, 39.82781182655966, 90.40151056564521, 15.567454799194547, 54.88851466835766, 109.778967357964, 45.55228653536312, 103.16902728388949, 15.823100894304012, 92.04889486393427, 6.005213659493739, 15.324271188839777, 59.63360253971818, 120.18920007265005, 82.06780795442725, 116.05354956530776, 6.728448298669941, 68.14266059235064, 19.341764251702333, 31.139112584737425, 15.866270053194649, 50.19979684355037, 26.323544629419544, 13.54367226133342, 42.8479630990655, 205.59059661332435, 25.20377929278292, 7.616680189020896, 23.671797779643743, 11.212338831504209, 8.47687944455657, 15.104814063464582, 34.52565074950944, 24.844221247988646, 93.21532413212995, 127.01382525345319, 40.436345571186294, 13.746745261358983, 20.314412541689293, 43.03359713194706, 100.8655735489136, 7.930469845882459, 101.66725511027254, 11.109656462647274, 5.934894726301829, 25.26918680187261, 31.01138930364673, 68.63025337769808, 7.832994023771672, 32.31559942097935, 39.5237027350131, 7.746857558561615, 16.114141317665776, 74.71932300107517, 12.333475791618511, 88.13191934343703, 19.44072166197472, 27.060761835328442, 98.12964802837895, 84.22530683267914, 7.734332076229037, 26.476535331596235, 7.831365593886226, 39.28223611726116, 12.704239603460275, 29.872406821388168, 24.968655669004526, 57.30471456331142, 17.034767237429357, 35.57120838143452, 12.859418060690793, 92.95441294732346, 26.700818855772003, 13.862069295263328, 5.233275081900455, 18.47783280835758, 70.24979391166111, 38.80160729123631, 55.07211710229137, 7.561893500837759, 18.641812509219353, 13.037861357397487, 30.18266787559698, 9.779567925906827, 12.76285912570457, 16.847813704553797, 64.63228231756085, 61.65047984476783, 6.954333105957183, 70.05161198863958, 23.301149835319805, 33.54194190420254, 5.59208510843894, 50.49665250263954, 5.662593596654664, 71.83215771744301, 17.45524505327815, 11.663022673407148, 69.61689931094872, 22.301781015730366, 7.100564398430099, 7.018948502166853, 40.96460613358928, 50.467720131430426, 21.404886892057064, 5.179471281368416, 13.182793956969649, 61.11997083923525, 37.880067285766714, 13.741604199185788, 7.3346630888231035, 56.54344289210581, 44.64074542939467, 61.15431016779998, 13.404952074256649, 37.088881956617215, 73.2336476943661, 45.60204712095868, 57.01947359508412, 18.823896215807075, 11.09802174933678, 23.743331274510318, 52.46389288836986, 12.348351290027548, 17.832800658052665, 5.692566550993013, 13.126749616478339, 5.766375814361699, 8.122439458269653, 11.225801816480574, 62.24644519952106, 15.445994568772756, 27.286855151902678, 48.24307000688432, 9.86168477927395, 10.24532938646214, 55.76439556790817, 78.22219250663505, 104.25951066294515, 25.258589423447926, 5.305686329531366, 33.7910799404562, 5.0476966913461165, 42.95029267443668, 40.72128290271927, 9.875271299938847, 7.4616696319151, 7.050323616381267, 6.7813681844031235, 117.00112311533869, 6.794477002095908, 93.10035706667966, 27.36726988931169, 18.236984138372307, 63.60044383695316, 6.120672321000391, 8.051181268988454, 25.43896272996389, 7.389125290979001, 43.9352636719009, 9.328216887246377, 118.19174877250668, 74.38679937171128, 35.07545161993774, 60.28757185266415, 15.554755528282385, 9.204473200443676, 76.95709074491481, 40.39158235913895, 58.760368808331336, 27.607424208482566, 62.62766159425012, 33.68121683336272, 43.479220042337985, 6.9781528693209145, 81.1941523999904, 68.21588594414116, 31.177081243559574, 14.078250094386643, 37.08228531952298, 62.878279417492244, 31.87910102294951, 5.397790162435022, 28.878874303665455, 63.546361626100214, 82.05195259778432, 46.49932285379889, 17.009644338127693, 65.30168994812132, 143.56890410675777, 35.521932838890514, 70.17353343848521, 52.87154829530011, 8.884418660666768, 7.453864284746731, 8.685276689138812, 102.35087113921301, 7.7542760816463145, 15.580175986879201, 50.74846257921592, 41.99699763321069, 11.709172070754896, 81.38331835406731, 28.057981243189598, 30.02283244473316, 17.780326468529026, 157.43515106932023, 38.59133467379219, 119.26915340769791, 40.9004231073895, 52.409174656482755, 9.578037414195288, 12.080942338989598, 57.65093914669735, 81.16590505104288, 85.40007875415573, 22.881114349167007, 26.930949681629084, 62.13686803014808, 9.787481818188544, 13.883826434670878, 61.81420022948039, 18.51386304072372, 7.433182218365194, 10.700327403959104, 10.412650918565499, 9.849128159384023, 50.269543890817964, 56.529775446271586, 95.87629464684606, 99.01673716032785, 37.60479313757054, 27.40149000198346, 23.981240040930775, 9.041718266848186, 26.805518785993684, 11.368927592830385, 45.96694086148856, 46.75996384262662, 49.11760854457064, 12.439751559786743, 40.27051856568281, 14.019792185382276, 44.00598180970635, 13.470362104809503, 32.72146062626551, 49.21827258130136, 12.844982742989215, 12.152525041876098, 82.64709255835143, 26.012051334013425, 20.766591601853285, 37.21053307604652, 15.515188493808616, 31.096948631745988, 8.520862582224797, 90.41354314205464, 52.70205513636013, 6.902002171450184, 42.20099822988445, 10.012776914026603, 34.35019136764661, 12.415902662002086, 46.006166599420276, 9.792568403398839, 132.43080375430898, 6.76509465758733, 19.63665898119953, 13.535669864585604, 16.90707177595369, 7.12910338273986, 87.24127014257014, 46.91265534611848, 9.145671677697651, 42.45596967658501, 5.66892394263845, 6.628968396886373, 110.46249707851658, 38.611911341949565, 20.729920802614416, 19.249692766052565, 24.618567885093498, 14.230080156891988, 19.07842113173955, 10.289889557600318, 44.1304769620009, 7.339418970411206, 42.43491677392248, 22.41933310063579, 72.34338140465739, 7.5283228334533, 25.84790064286748, 6.381079578615405, 111.43637184664681, 7.225134688878932, 6.954666379161702, 96.86224645499112, 30.785141655415682, 8.474971762122934, 11.739584380918211, 55.056540076223214, 5.314992444920363, 18.48708417838469, 75.6664027935307, 12.05258108684769, 10.656044497309647, 10.404997800876043, 56.05349345075088, 67.96456817486609, 81.22657688654961, 49.42016733799891, 12.296689268872404, 26.161879412843955, 45.38508459961883, 20.427432493791112, 22.517226964177972, 21.58333624988548, 11.139142289263425, 59.060892753653555, 6.147958673166586, 8.496137062902688, 90.30390066336565, 39.600732521789794, 67.31049055069118, 124.35738100255621, 21.400152285708067, 19.64085404725156, 121.12299237559941, 68.20313316995285, 46.19435799477614, 15.134919452917416, 19.886781820820794, 27.162238637196495, 14.214072938293056, 115.14906277809592, 7.045879585977212, 142.95353742968172, 31.60500151064182, 16.67321112133154, 11.45311682700216, 8.222604877953065, 16.61028887905804, 21.49238866585931, 10.323817461630282, 59.707043868910574, 106.6933307669526, 42.45785193152026, 10.25874787538329, 23.753515835725985, 15.842270834389836, 54.247421554061184, 11.718015915453655, 9.251826900667709, 55.555421571855824, 5.182798693038229, 11.287845030306674, 21.609064889339955, 28.68418149825625, 58.57646450613494, 14.262417005490608, 16.612572203656686, 16.236087052979443, 7.07012245466578, 36.74755115014632, 5.56955315814885, 38.91744749703327, 78.85738379606225, 66.55885619672107, 39.64570920071233, 59.29291459700674, 43.05116230352876, 57.616583933068526, 13.018545934843285, 38.0480565816235, 24.925011809751684, 48.79848066987436, 48.274437799688435, 19.76507371158679, 8.565946758021937, 59.157588927921495, 57.14227355593507, 9.275322842239122, 14.083832907023066, 14.365127431573358, 16.940142017461092, 33.65000424197482, 22.38482569043266, 9.8423764229321, 103.30964749356153, 93.482737580894, 70.28695282389145, 36.95272628687471, 13.763226292274851, 66.32050179751066, 105.52022887233039, 20.58372196333329, 57.306558616504354, 7.2055606193971675, 6.901695760628813, 29.114219935320538, 5.196310986913392, 64.93439841265513, 66.54164242492251, 8.601749250514457, 16.577497833620654, 21.034827730068592, 71.15904863888585, 71.14756778662918, 83.33258945489067, 40.09744323758643, 10.758686591572488, 32.71729187462118, 72.60928034296013, 23.182448164356757, 20.02688139994615, 35.64640238564675, 19.619730168894044, 10.387210959701576, 108.49878733131999, 75.39827125131204, 8.518593311533348, 50.09760170113625, 5.783556515088281, 7.7226393825025, 17.362973213161744, 80.6800037352902, 45.02780319648969, 6.515288961328326, 21.02737688071963, 54.36451013450317, 57.52928577866641, 24.29359490456843, 49.28102430125691, 35.8859259513479, 8.56080309083324, 12.854971881557413, 8.88725896384838, 21.0800477677967, 98.2668911820237, 14.344215517868456, 20.26298700390288, 7.314621461467949, 54.44047864778453, 45.24597043884262, 24.081153968623944, 11.605888653834558, 81.92891659436373, 11.077075817245964, 72.47645741895147, 14.150248660347113, 38.98791894640125, 8.310910804301987, 14.562595175466456, 113.34002985790485, 45.54586004938276, 8.255063789030922, 15.748738231301685, 12.334676676031608, 37.13449911745443, 25.16227956950337, 76.85488255560675, 65.08240181335337, 13.072862840680108, 6.511161887246643, 36.836853509896514, 32.28105111670073, 9.011200476779752, 67.99261973751265, 43.56414154407412, 22.664659544292277, 10.641535346023423, 34.64376026085962, 34.85294630448037, 13.484231145440194, 28.36688965478537, 58.00923280980163, 35.763347543902626, 10.676355788751687, 21.408481350254558, 6.935108643188424, 11.81394661943359, 22.301879646778744, 5.798056980720284, 94.08282554710377, 22.67539113883788, 22.799743419864125, 11.478324417670079, 53.32099708503726, 30.313557087185988, 19.386367771690736, 18.844755698317343, 44.82455210759767, 15.373962501989634, 9.515152632133866, 16.939112837300648, 28.396269204357306, 10.335418544876822, 89.17929630357061, 16.840582885523023, 14.764919650394848, 19.293637630127673, 10.818970722833296, 17.912239148805387, 22.802790503154316, 23.11994332941388, 24.98230979382431, 24.23010562508092, 10.571739601893785, 13.638015383570616, 8.450786848419284, 5.422272727496316, 30.803919297691557, 6.16787773870539, 58.70960199565301, 21.757659586333236, 6.354290042601948, 79.01645583162893, 33.34361337283414, 10.43402929171163, 18.463413928011903, 51.398250308694465, 14.286148944932568, 39.54018858552626, 14.718239656426421, 7.133358728257916, 9.039139006726538, 101.67541212273028, 22.554844581646957, 81.56092029014815, 98.24038281892368, 91.74922363819401, 7.3638018403555074, 49.324132034224064, 52.30603773350151, 12.546541696075952, 55.03255187031403, 5.110580889802128, 13.547640273474753, 5.978462608040431, 140.66736056961003, 6.8420134193939175, 10.315950757587592, 47.18125856105616, 5.635843547192914, 6.538968180383621, 21.88992975546214, 28.648955299399073, 26.749893324709912, 69.50432749701974, 51.68692180174859, 12.564852483727488, 10.081517978439976, 84.29318403235067, 14.031376121318083, 8.789471837751982, 6.70667215533237, 44.29845720927677, 7.67457766137185, 80.51883119795062, 5.404723250860584, 22.0490516942453, 12.716189178806728, 33.87083847411847, 28.595637073850124, 5.167546674236084, 8.480139111408706, 24.72735125141671, 44.190865180027814, 15.528885568827413, 11.513728945059556, 34.379631842273476, 17.764093391162096, 110.36408081759643, 10.112704532840677, 5.6086556117359905, 14.575191366941128, 9.410796500229058, 45.32471755036951, 5.245340531125099, 6.464627593235395, 41.34674346917381, 25.368012276144974, 26.56679033128959, 78.43609293008124, 101.68991413483468, 88.73750329161878, 19.75121208094216, 77.54669677668076, 21.242898363177446, 26.471433912750108, 51.41048584252594, 88.97567150220965])
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);
([6743402.047899941, 6798381.602915615, 6830177.945460082, 7077994.625212399, 7158199.763187827, 7282689.0625, 7290009.911801992, 7292031.267858106, 7311199.491457405, 7322373.542119066, 7325259.086153767, 7592617.9120263895, 7619414.0625, 7650001.252156334, 7818962.844233887, 7892179.3971398175, 7993259.776386856, 7993287.783796568, 8061903.125, 8099346.142670247, 8147526.255549074, 8191301.684435399, 8198802.344454739, 8243929.6875, 8371743.8449094165, 8559257.8125, 8565608.889043516, 8565949.27115848, 8568086.754444497, 8666300.048107455, 8673133.781371502, 8673821.875, 8675409.375, 8675463.81485712, 8675915.625, 8677207.594032718, 8688659.316214968, 8696385.935082313, 8795041.94489316, 8806132.534069385, 8819115.625, 8821127.23927495, 8828546.705870736, 8829089.592235776, 8829220.676252192, 8829486.473690713, 8829860.193799892, 8830262.433225669, 8831203.125, 8831284.06113106, 8831308.088159239, 8833877.79661033, 8833913.354696944, 8837296.30220899, 8862834.578185035, 8875394.782318195, 8877307.709458087, 8883792.59997623, 8885794.745985812, 8887077.432332234, 8888476.144987103, 8889162.532126136, 8889599.69769525, 8890515.625, 8915554.421248462, 8916862.5, 8917558.547532491, 8919004.91653882, 8920257.707815751, 8921026.523484185, 8921047.204317106, 8921781.25, 8921909.127004456, 8922135.82795166, 8922645.249725224, 8922682.223801712, 8922697.036171935, 8922722.321188247, 8922774.338740181, 8922932.971570926, 8923010.224360788, 8923025.262222094, 8923128.92335266, 8923187.21360701, 8923301.38280371, 8923522.226092564, 8923804.6875, 8924259.684525425, 8924300.0, 8924482.8125, 8924681.478319036, 8924844.117997319, 8924882.814891769, 8924952.352818947, 8925832.459919259, 8926000.839751719, 8928321.875, 8931793.671313053, 8931845.1362293, 9005839.434806848, 9010785.613939468, 9037091.566264534, 9079485.978516554, 9113899.469612477, 9129703.78641302, 9139654.0289821, 9141640.625, 9150921.875, 9155134.217347628, 9163389.32856721, 9173979.265661016, 9227791.578400794, 9229456.25, 9231723.4375, 9233507.8125, 9236501.5625, 9237008.292103166, 9238279.027579477, 9262815.625, 9275746.875, 9278123.859842023, 9290193.365947057, 9290581.21718795, 9293156.286744399, 9321256.25, 9326245.896215344, 9361627.503131231, 9363169.825398516, 9393771.668867664, 9394259.14901842, 9394879.326302538, 9395200.344477037, 9395526.145135935, 9395536.79565682, 9395728.77775762, 9396083.838617431, 9396297.465110872, 9397680.187797802, 9399192.966858855, 9442300.0, 9448465.869956706, 9449515.625, 9449869.31972607, 9450339.547300255, 9450360.931962002, 9451599.106036844, 9451673.152522145, 9451799.269565966, 9451939.555731645, 9452510.216497613, 9452637.258051476, 9452704.6875, 9453697.93623441, 9453769.84482801, 9455366.549606072, 9455778.125, 9455916.153884234, 9456256.112203613, 9456472.50418881, 9456614.028758012, 9457790.084885208, 9458837.79887042, 9471321.49382651, 9486562.793315124, 9487148.133209087, 9487582.8125, 9488648.308595479, 9489118.789938202, 9489733.056391979, 9489857.46482766, 9489869.261599962, 9490036.282599531, 9491150.0, 9491831.25, 9493293.724922678, 9493450.600779254, 9494657.512572423, 9494882.408332054, 9495812.5, 9497845.150561508, 9498379.216677101, 9511203.896424336, 9512483.749137603, 9512506.216105068, 9512577.043258594, 9516327.24100952, 9517293.75, 9520831.509894272, 9520838.100977363, 9520887.285703335, 9528815.049227403, 9541395.956358386, 9547688.11173419, 9549129.44297439, 9549323.074040601, 9549505.7887901, 9549530.055542015, 9549634.130297204, 9550012.5, 9550361.336093327, 9551120.3125, 9551311.873281494, 9551546.875, 9552032.8125, 9552039.718232624, 9552156.472794136, 9589113.481928477, 9595620.3125, 9607625.464300469, 9632885.67079013, 9667444.472582666, 9670919.70969289, 9673465.504699042, 9675728.38105047, 9676754.758390075, 9677471.566221675, 9677641.649909917, 9677748.295663418, 9677812.179909376, 9678812.615135925, 9679031.162448786, 9679152.133787747, 9679246.205583414, 9680095.3125, 9680585.689703846, 9681243.220494457, 9688088.874592062, 9690758.949480427, 9693903.621325847, 9697162.155219382, 9697215.625, 9697216.35321655, 9738272.905085005, 9755259.960072728, 9761356.064686337, 9785026.695268262, 9813548.4375, 9841014.732487738, 9844292.041361872, 9844482.576824209, 9845581.38642707, 9847800.0, 9848280.67154471, 9848316.891865669, 9850382.019590857, 9851938.06102836, 9857828.681101697, 9861078.021794207, 9864912.417904591, 9865471.800352743, 9865520.913177595, 9867945.564919032, 9868532.8125, 9869421.274775546, 9870286.863985201, 9870940.969802054, 9871790.262409143, 9871904.307981879, 9872725.247330341, 9874207.354109738, 9874225.0, 9874258.334242392, 9874870.624736222, 9875300.0, 9876315.355329145, 9876456.25, 9877753.965418749, 9879672.929685753, 9881135.248643547, 9885026.649101434, 9887716.91539338, 9905766.364989016, 9918983.555378642, 9919048.707737828, 9921847.546116821, 9927783.084713886, 9929404.6875, 9944340.625, 9944347.817351045, 9952073.251897117, 9962678.125, 9963928.125, 9969267.342571348, 9970889.888190618, 9972737.5, 9973657.228670252, 9975851.5625, 9976811.151044188, 9979032.8125, 9980025.725869425, 9980478.221923847, 9982718.59516545, 9984211.652536118, 9997553.125, 9997617.1875, 9999187.321820846, 10001132.8125, 10001358.386192853, 10001573.16223938, 10001589.285662187, 10002332.8125, 10002433.307225011, 10035748.874088014, 10067467.887844063, 10081029.325270029, 10081437.864537299, 10082385.97451424, 10084497.20998952, 10092770.365978505, 10197839.40582428, 10282825.04219514, 10421987.5, 10551493.531040112, 10559381.248354742, 10645693.75, 10645712.483247329, 10646298.084241884, 10654005.623125719, 10663141.278514914, 10663745.219320936, 10680872.697373778, 10683051.563900454, 10683396.863421177, 10685053.073056784, 10716763.601832153, 10720984.375, 10721437.237904074, 10721481.84350176, 10722329.6875, 10722676.188178983, 10723466.005365878, 10730434.652022725, 10747680.888836363, 10747683.467603015, 10752866.650712062, 10754037.523062088, 10760308.497057369, 10779490.625, 11016616.564249625, 11017715.973472947, 11025346.652957622, 11025805.277069805, 11034679.6875, 11046143.006089957, 11048287.491699563, 11051360.727217836, 11051699.564975418, 11053482.595030405, 11054983.59835332, 11055007.132792914, 11062893.897731729, 11064693.102195082, 11067666.872203907, 11068452.695268223, 11069261.552363615, 11069541.43331458, 11070129.6875, 11078495.3125, 11079040.509224212, 11079910.094702065, 11079941.578790737, 11081456.25, 11081509.861383995, 11084193.563882072, 11084314.48643726, 11088002.530185007, 11088805.33647758, 11094739.948669706, 11096678.125, 11097071.863696026, 11097179.65920913, 11099583.971130332, 11100087.5, 11101771.157947425, 11103187.00023111, 11107595.845638173, 11109972.192455528, 11111576.5625, 11112555.827417359, 11114137.419940941, 11114871.470160896, 11117139.0625, 11117921.875, 11118027.616111327, 11118238.095704166, 11119807.259347642, 11120008.678784892, 11120142.211577652, 11123608.217645627, 11221038.253414359, 11224935.422531042, 11227424.675069243, 11233421.535170479, 11243921.40356013, 11244980.222639656, 11245903.556423357, 11246008.37408917, 11246260.751260443, 11246783.880057987, 11247253.105925048, 11247388.761919646, 11248786.09665295, 11248786.599453673, 11248953.721059587, 11249387.27468009, 11249557.202176694, 11249780.585204218, 11249830.926156651, 11249923.499364562, 11250673.4375, 11250926.688847262, 11251054.945157798, 11251100.0, 11251136.504450615, 11251288.413683478, 11251642.1875, 11251689.234674538, 11251873.4375, 11251969.15617041, 11252055.153007187, 11252314.302013319, 11252725.361950723, 11252806.25, 11253150.0, 11253470.892550323, 11253715.215855062, 11253962.222325455, 11254337.588814521, 11254529.6875, 11254547.028175132, 11254753.125, 11254809.313665938, 11254824.89800638, 11255043.755572073, 11255087.807511289, 11255312.026655622, 11255341.25025888, 11255354.053990059, 11255478.125, 11255943.746229323, 11257713.00638248, 11258909.058294, 11259077.645152293, 11259133.89665048, 11259226.714597218, 11259404.53347181, 11259458.467007907, 11259509.375, 11259540.625, 11259549.16405323, 11259929.6875, 11260004.6875, 11260180.027693259, 11260332.237797504, 11260428.125, 11260445.3125, 11260452.04862387, 11260523.190944107, 11260563.988868613, 11260667.1875, 11260690.625, 11260942.169599727, 11260947.130013453, 11260955.401367022, 11260993.75, 11261171.432129772, 11262102.016553279, 11262173.05804178, 11262306.957485372, 11262389.0625, 11262413.235523844, 11262451.5625, 11262582.231918441, 11262621.442216007, 11262706.93482574, 11262717.098555, 11262783.211645022, 11262902.169929666, 11262916.03506481, 11262981.25, 11263214.405164868, 11263391.079617755, 11263479.168062484, 11263535.32757945, 11263547.486920984, 11263643.007283958, 11263673.236488672, 11263673.4375, 11263720.3125, 11263926.5625, 11263981.527958298, 11263986.472287355, 11264203.125, 11264241.716405364, 11264510.03419306, 11264737.933939505, 11265021.231589755, 11265325.0, 11265350.49425418, 11265734.066763027, 11265845.120571652, 11266047.576086733, 11266104.191502687, 11266123.145713773, 11266205.738694292, 11266508.367286494, 11266551.5625, 11266688.976482674, 11266749.200895464, 11266934.998156419, 11267168.75, 11267518.265620317, 11267670.030333158, 11267837.1363112, 11267954.594134295, 11268284.524032107, 11268362.5, 11268856.25, 11269122.941708872, 11269447.578494463, 11269635.644399654, 11270027.84429153, 11270029.972819598, 11270209.738215491, 11270540.625, 11270912.5, 11270936.545741724, 11270953.694681376, 11271130.923449667, 11271573.296061149, 11271637.176577523, 11271660.9375, 11271697.017787416, 11271758.690429682, 11272092.605731577, 11272304.79168194, 11272404.6875, 11272722.531271096, 11272748.4375, 11272758.521600097, 11272879.540220736, 11273017.880245928, 11273539.0625, 11273747.37430062, 11273860.686995288, 11274009.382949473, 11274219.567241652, 11274394.887295455, 11275032.825021883, 11275064.091528695, 11275096.875, 11275519.695452455, 11275656.021861097, 11275784.375, 11275843.138644956, 11275849.286556767, 11276034.042185161, 11276076.528115215, 11276458.79493328, 11276852.366278702, 11277007.179686116, 11277195.112426566, 11277610.1557915, 11277665.013559746, 11277814.764031444, 11278172.644000186, 11278181.25, 11278524.588094536, 11279001.5625, 11279029.475510178, 11279175.348068882, 11279477.413664622, 11280029.6875, 11280276.942770876, 11280307.13760438, 11280514.0625, 11280612.905062014, 11281131.590039296, 11281212.859105937, 11281338.0254681, 11281939.817494188, 11282253.822148083, 11282336.208407873, 11282506.770783184, 11282890.300723907, 11283038.545070168, 11283476.5625, 11283593.75, 11283860.537690671, 11284084.685542293, 11284192.825551659, 11284390.492906354, 11284569.620365875, 11284636.886096908, 11284692.306522377, 11284793.258998405, 11285271.88230377, 11285352.090823002, 11285438.05660037, 11285508.18881671, 11285784.375, 11285937.766805628, 11286274.963756703, 11286323.51819735, 11286358.865631655, 11286564.916849026, 11286599.70387126, 11286833.56318754, 11287138.385939635, 11287366.314209733, 11287371.19734061, 11287404.059194844, 11287863.857593374, 11287970.46836263, 11288402.829492657, 11288639.261809329, 11289400.696880508, 11289941.01555833, 11290562.698667949, 11290684.252035806, 11290841.003296647, 11290843.265150448, 11290893.16664481, 11290932.8125, 11291570.277706416, 11291614.0625, 11292246.237497848, 11292673.137717148, 11292836.15532412, 11292890.625, 11293444.09745987, 11294413.966865301, 11294434.375, 11294505.420791708, 11294600.378517732, 11295079.199046861, 11295170.3125, 11295379.073242508, 11298315.434642745, 11298508.191243798, 11298673.740197107, 11298751.363568773, 11298751.386331443, 11298869.479272414, 11298884.2235384, 11299041.57737906, 11299295.3125, 11299351.225879807, 11299458.993460516, 11299513.968991028, 11299763.922689699, 11300493.1337465, 11300568.43669597, 11300592.110096894, 11300617.023694413, 11300625.0, 11300800.892578958, 11300831.743397139, 11300859.013451736, 11301268.75, 11301432.745423332, 11301527.880178917, 11301624.840088332, 11301831.69121393, 11301915.625, 11301924.584974673, 11302039.610722868, 11302134.375, 11302412.443471853, 11302415.25509626, 11302803.125, 11302905.902223945, 11302983.247437404, 11303037.5, 11303161.182508612, 11303168.75, 11303195.623084964, 11303306.904897362, 11303314.342610575, 11303530.968579883, 11304043.75, 11304339.98514398, 11304681.801956072, 11304805.05935668, 11304808.827799581, 11305057.8125, 11305139.0625, 11305198.77722968, 11305401.5625, 11305468.500346933, 11305502.99305331, 11305945.987687774, 11305960.5916688, 11306022.521296643, 11306342.1875, 11306514.9141842, 11307116.655777067, 11307154.6875, 11307454.352491166, 11307526.97690535, 11307893.195201531, 11308052.705176538, 11308634.69348757, 11308806.767206019, 11308979.6875, 11309031.25, 11309154.865051925, 11309256.538312335, 11309265.625, 11309285.729680665, 11309630.99134203, 11309749.769177794, 11309964.0625, 11310441.109902557, 11310474.31848592, 11310761.92346772, 11311055.70757198, 11311094.08640303, 11311220.146440255, 11311239.949898243, 11311400.687187698, 11311829.626127277, 11312009.156366255, 11312634.53075563, 11312836.223302364, 11313051.087343905, 11313079.654864563, 11313247.215488696, 11313331.025665622, 11313492.1875, 11313854.733233238, 11314059.27190385, 11314092.1875, 11314126.482705917, 11314179.627428213, 11314233.009305712, 11314252.218157006, 11314373.4375, 11314429.490514433, 11314443.75, 11314546.875, 11314596.279876657, 11314676.714352708, 11314685.053204091, 11314808.761223745, 11314813.657872908, 11315124.739334922, 11315214.0625, 11315249.6099362, 11315284.039979445, 11315345.691438667, 11315346.875, 11315490.181759084, 11315657.804461807, 11315731.088141968, 11315746.875, 11315757.503182415, 11315885.428975312, 11316336.740289332, 11316342.58624331, 11316685.557988869, 11316762.282451184, 11316893.805212267, 11316953.606789399, 11317383.996715065, 11317508.025802061, 11317564.327525081, 11317939.0625, 11318050.499148585, 11318467.936054807, 11318557.37125582, 11318573.49468094, 11318650.209601376, 11319079.602021655, 11319407.87428653, 11319507.238675147, 11319531.084129144, 11320023.86007344, 11320290.59555386, 11320511.747091891, 11320885.730573282, 11320900.608852735, 11321083.982138049, 11321138.775027994, 11321326.597959816, 11321679.424123392, 11321871.487594794, 11321880.860088378, 11321933.110629348, 11322084.983885553, 11322146.38832061, 11322150.733748997, 11322204.6875, 11322498.282508334, 11322503.784577508, 11322565.477443617, 11322567.851905419, 11322635.44322009, 11322764.756154217, 11323065.309172725, 11323591.256887442, 11323717.1875, 11323784.315162808, 11323842.232603664, 11324790.25787179, 11324912.5, 11325099.600231985, 11325108.947533477, 11325198.191597437, 11325373.4375, 11325510.223137869, 11325921.875, 11326120.514350247, 11326199.188358214, 11326356.25, 11326502.16888725, 11326563.641463842, 11326720.750321187, 11326844.780205414, 11326869.130916703, 11327060.347661238, 11327192.15255655, 11327250.387432713, 11327740.948140536, 11327852.374336302, 11327885.9375, 11327893.653375106, 11327936.22266685, 11328158.623090526, 11328416.383855976, 11329900.771297958, 11330163.080988491, 11331813.801299619, 11331863.34358504, 11332116.895992663, 11332218.75, 11332723.630529251, 11332747.550000314, 11332892.350011595, 11332959.497017317, 11333099.261118433, 11333450.431331258, 11333453.125, 11334217.281076934, 11335158.161430348, 11335610.710293462, 11336081.25, 11336962.133895578, 11337790.858682366, 11338104.6875, 11338134.959872566, 11338138.312320285, 11339244.549027823, 11339499.902185172, 11339782.712953677, 11339825.09242571, 11339981.25, 11340043.756974308, 11340237.35518677, 11341339.92281948, 11346467.0043743, 11346700.24524572, 11346721.717291264, 11346832.657236839, 11347754.6875, 11347891.779960668, 11349281.25, 11350529.760357166, 11351874.886813449, 11353292.356322328, 11358829.312220389], [7.2551750148893595, 119.84383793028702, 152.02153261831268, 19.808276643067572, 24.59289029283166, 75.44705128695676, 46.338376023028076, 15.010015849639794, 112.40892167580526, 103.13657596892365, 50.390367303494905, 18.675349424312092, 81.20371827475279, 65.7895125410652, 7.732674171024905, 34.590663695327244, 49.46531620331742, 6.532518491591088, 32.615170170776594, 277.3564422309868, 5.913112831998846, 11.32380294461667, 39.34793985805817, 66.2155880584084, 76.81471423241454, 38.065188071998584, 5.676513030506049, 29.76094965996174, 29.753206151926378, 11.203743746378201, 71.29376680013323, 44.32412303061886, 67.46423853372754, 60.60973557263923, 59.006462212079015, 31.362549711884537, 27.413182174151263, 55.77323803075808, 10.768551040237826, 13.97680004256283, 41.73985906466917, 23.710292347494633, 51.300035993910456, 37.507123316706796, 48.84942818801527, 17.473073411300025, 10.135148885218918, 11.024914003806693, 29.64407847211484, 5.791116009965499, 59.08793678344231, 9.068354007916032, 5.5760511720110095, 52.40690793900788, 62.617475702191236, 21.197213093694447, 53.37694701490977, 36.30469338997641, 10.34434817823673, 61.482063613761994, 19.532311981345952, 11.17829607640673, 28.461576141567573, 31.444104205013755, 24.191621095843363, 47.33730737762052, 7.32907334274893, 23.39136475102014, 14.556970002147189, 60.89516338060983, 11.68466485457646, 49.11482084706725, 14.07362199461672, 9.240460343601601, 11.030983253017087, 25.774176741193806, 28.852850342375653, 63.67769920769554, 46.991416162662645, 20.62470980609154, 19.60136020340172, 16.0449717292396, 7.348416712332136, 30.71572821789934, 35.3094471586983, 16.47295404779002, 31.10215081601645, 94.07277693719828, 40.681733736930596, 52.99265209715478, 15.889549802406743, 39.26115242454695, 22.407837426703807, 26.801561126573198, 23.700781614822958, 19.525853086653193, 40.83170295956382, 87.34567969664721, 20.758636079561356, 17.665059866615703, 106.5287325853642, 6.967192788291198, 69.05740554490093, 13.24358516875305, 7.306893724973408, 26.476468120396298, 102.83142283729399, 49.23885817331956, 106.9361180896004, 21.92483444268057, 8.998487611282384, 61.10910517115368, 87.94384146554613, 71.43822733051148, 64.46489362784186, 88.13874858531274, 9.772767843772336, 7.965927848615289, 82.35413302391424, 56.31000655868047, 23.546036465511783, 13.65522569575368, 17.978265394948675, 14.240967628299853, 57.58020766290863, 22.890574934969663, 36.581951273712434, 10.443088003732775, 30.46546975797324, 16.340322314119383, 20.337500878021267, 16.440020075645418, 48.27188570409356, 5.225850714693417, 18.19212282761256, 21.13768324069304, 12.195052264775367, 78.55501320985148, 12.612168809081307, 62.39134498822649, 5.073377977682033, 29.19176356431536, 95.60283427197572, 17.179242174737347, 56.684480917974426, 93.73591249677267, 60.29913467728069, 39.65856475307234, 67.95927301137624, 79.16173016718072, 27.44467852157828, 46.096106101218176, 15.466094554430779, 13.415465531118151, 52.58182594937276, 31.218346706927125, 34.872521257131105, 9.217544436979303, 18.040457987006494, 76.67742221449961, 21.1020546026178, 60.57571209033055, 6.220211789558939, 8.448026798879964, 25.081485800102083, 48.75010732841402, 61.07373678633748, 6.771506200059151, 7.187074053548297, 49.64300134651855, 72.2239102053945, 34.22013662407855, 51.4982727928737, 45.35508767258757, 9.310288516846636, 5.362938467121785, 49.63054499825393, 38.273206736868644, 68.62634572841412, 106.42520308636469, 13.359964691940469, 10.081659345754693, 25.38723218198098, 14.81387268275413, 9.773925351541466, 158.56677975490172, 31.023392791961015, 5.962362818948428, 6.585615383199932, 14.766166506176322, 12.919333212818835, 7.8539843355734, 18.440293109212696, 24.48619354534041, 106.32836849468333, 78.51954787606529, 6.839277112307645, 5.6913594651509465, 95.32132024557711, 91.49995466502558, 46.69425747577368, 84.16756655876631, 51.07290128967238, 43.00975297043062, 53.29026245175248, 7.181275461351208, 54.93063701619354, 49.36148856187985, 8.548379858986479, 24.185499512661647, 43.094068821741814, 60.36292559246701, 42.71850066793812, 19.677826608111953, 15.21877337740057, 23.6392546254641, 32.8542468698288, 32.91196058839547, 163.70939304199646, 9.430692259681553, 26.905444956890424, 55.11795460534596, 25.909424152652996, 38.05071895487867, 84.11182712928706, 28.46226264904394, 25.2965681138856, 19.02536901466238, 24.775383421298994, 10.628413133526731, 37.86470464906046, 6.918152509144184, 21.175725025632328, 9.880498084302417, 25.335264902750144, 15.472917318338405, 76.62611054746354, 18.305139742686595, 10.906716425885788, 48.29591315864786, 32.92226940591586, 45.39191067241269, 13.116191400219549, 57.545044459428716, 6.614900554404024, 13.633592411358546, 8.60430684218417, 106.50570169238556, 6.747226790082193, 114.90516636494061, 6.585563989063977, 46.61321218395169, 80.40342349446172, 78.31252445233294, 118.48910607478383, 17.552962032910003, 16.609109604169326, 24.26649104001345, 27.242095360639134, 69.5306202430092, 59.982152272939956, 13.001035750187206, 9.610149749401783, 54.49936463506427, 24.146042700358212, 51.59479189230557, 19.548502452657754, 5.675183990048977, 19.79555978680541, 30.92082965212629, 119.68803313406201, 45.50634495702023, 104.98674740063788, 6.3019555762509185, 18.312610111152228, 38.436439219004384, 34.40920683110593, 66.97578722502051, 13.008605977643963, 9.17924605938643, 39.82781182655966, 90.40151056564521, 15.567454799194547, 54.88851466835766, 109.778967357964, 45.55228653536312, 103.16902728388949, 15.823100894304012, 92.04889486393427, 6.005213659493739, 15.324271188839777, 59.63360253971818, 120.18920007265005, 82.06780795442725, 116.05354956530776, 6.728448298669941, 68.14266059235064, 19.341764251702333, 31.139112584737425, 15.866270053194649, 50.19979684355037, 26.323544629419544, 13.54367226133342, 42.8479630990655, 205.59059661332435, 25.20377929278292, 7.616680189020896, 23.671797779643743, 11.212338831504209, 8.47687944455657, 15.104814063464582, 34.52565074950944, 24.844221247988646, 93.21532413212995, 127.01382525345319, 40.436345571186294, 13.746745261358983, 20.314412541689293, 43.03359713194706, 100.8655735489136, 7.930469845882459, 101.66725511027254, 11.109656462647274, 5.934894726301829, 25.26918680187261, 31.01138930364673, 68.63025337769808, 7.832994023771672, 32.31559942097935, 39.5237027350131, 7.746857558561615, 16.114141317665776, 74.71932300107517, 12.333475791618511, 88.13191934343703, 19.44072166197472, 27.060761835328442, 98.12964802837895, 84.22530683267914, 7.734332076229037, 26.476535331596235, 7.831365593886226, 39.28223611726116, 12.704239603460275, 29.872406821388168, 24.968655669004526, 57.30471456331142, 17.034767237429357, 35.57120838143452, 12.859418060690793, 92.95441294732346, 26.700818855772003, 13.862069295263328, 5.233275081900455, 18.47783280835758, 70.24979391166111, 38.80160729123631, 55.07211710229137, 7.561893500837759, 18.641812509219353, 13.037861357397487, 30.18266787559698, 9.779567925906827, 12.76285912570457, 16.847813704553797, 64.63228231756085, 61.65047984476783, 6.954333105957183, 70.05161198863958, 23.301149835319805, 33.54194190420254, 5.59208510843894, 50.49665250263954, 5.662593596654664, 71.83215771744301, 17.45524505327815, 11.663022673407148, 69.61689931094872, 22.301781015730366, 7.100564398430099, 7.018948502166853, 40.96460613358928, 50.467720131430426, 21.404886892057064, 5.179471281368416, 13.182793956969649, 61.11997083923525, 37.880067285766714, 13.741604199185788, 7.3346630888231035, 56.54344289210581, 44.64074542939467, 61.15431016779998, 13.404952074256649, 37.088881956617215, 73.2336476943661, 45.60204712095868, 57.01947359508412, 18.823896215807075, 11.09802174933678, 23.743331274510318, 52.46389288836986, 12.348351290027548, 17.832800658052665, 5.692566550993013, 13.126749616478339, 5.766375814361699, 8.122439458269653, 11.225801816480574, 62.24644519952106, 15.445994568772756, 27.286855151902678, 48.24307000688432, 9.86168477927395, 10.24532938646214, 55.76439556790817, 78.22219250663505, 104.25951066294515, 25.258589423447926, 5.305686329531366, 33.7910799404562, 5.0476966913461165, 42.95029267443668, 40.72128290271927, 9.875271299938847, 7.4616696319151, 7.050323616381267, 6.7813681844031235, 117.00112311533869, 6.794477002095908, 93.10035706667966, 27.36726988931169, 18.236984138372307, 63.60044383695316, 6.120672321000391, 8.051181268988454, 25.43896272996389, 7.389125290979001, 43.9352636719009, 9.328216887246377, 118.19174877250668, 74.38679937171128, 35.07545161993774, 60.28757185266415, 15.554755528282385, 9.204473200443676, 76.95709074491481, 40.39158235913895, 58.760368808331336, 27.607424208482566, 62.62766159425012, 33.68121683336272, 43.479220042337985, 6.9781528693209145, 81.1941523999904, 68.21588594414116, 31.177081243559574, 14.078250094386643, 37.08228531952298, 62.878279417492244, 31.87910102294951, 5.397790162435022, 28.878874303665455, 63.546361626100214, 82.05195259778432, 46.49932285379889, 17.009644338127693, 65.30168994812132, 143.56890410675777, 35.521932838890514, 70.17353343848521, 52.87154829530011, 8.884418660666768, 7.453864284746731, 8.685276689138812, 102.35087113921301, 7.7542760816463145, 15.580175986879201, 50.74846257921592, 41.99699763321069, 11.709172070754896, 81.38331835406731, 28.057981243189598, 30.02283244473316, 17.780326468529026, 157.43515106932023, 38.59133467379219, 119.26915340769791, 40.9004231073895, 52.409174656482755, 9.578037414195288, 12.080942338989598, 57.65093914669735, 81.16590505104288, 85.40007875415573, 22.881114349167007, 26.930949681629084, 62.13686803014808, 9.787481818188544, 13.883826434670878, 61.81420022948039, 18.51386304072372, 7.433182218365194, 10.700327403959104, 10.412650918565499, 9.849128159384023, 50.269543890817964, 56.529775446271586, 95.87629464684606, 99.01673716032785, 37.60479313757054, 27.40149000198346, 23.981240040930775, 9.041718266848186, 26.805518785993684, 11.368927592830385, 45.96694086148856, 46.75996384262662, 49.11760854457064, 12.439751559786743, 40.27051856568281, 14.019792185382276, 44.00598180970635, 13.470362104809503, 32.72146062626551, 49.21827258130136, 12.844982742989215, 12.152525041876098, 82.64709255835143, 26.012051334013425, 20.766591601853285, 37.21053307604652, 15.515188493808616, 31.096948631745988, 8.520862582224797, 90.41354314205464, 52.70205513636013, 6.902002171450184, 42.20099822988445, 10.012776914026603, 34.35019136764661, 12.415902662002086, 46.006166599420276, 9.792568403398839, 132.43080375430898, 6.76509465758733, 19.63665898119953, 13.535669864585604, 16.90707177595369, 7.12910338273986, 87.24127014257014, 46.91265534611848, 9.145671677697651, 42.45596967658501, 5.66892394263845, 6.628968396886373, 110.46249707851658, 38.611911341949565, 20.729920802614416, 19.249692766052565, 24.618567885093498, 14.230080156891988, 19.07842113173955, 10.289889557600318, 44.1304769620009, 7.339418970411206, 42.43491677392248, 22.41933310063579, 72.34338140465739, 7.5283228334533, 25.84790064286748, 6.381079578615405, 111.43637184664681, 7.225134688878932, 6.954666379161702, 96.86224645499112, 30.785141655415682, 8.474971762122934, 11.739584380918211, 55.056540076223214, 5.314992444920363, 18.48708417838469, 75.6664027935307, 12.05258108684769, 10.656044497309647, 10.404997800876043, 56.05349345075088, 67.96456817486609, 81.22657688654961, 49.42016733799891, 12.296689268872404, 26.161879412843955, 45.38508459961883, 20.427432493791112, 22.517226964177972, 21.58333624988548, 11.139142289263425, 59.060892753653555, 6.147958673166586, 8.496137062902688, 90.30390066336565, 39.600732521789794, 67.31049055069118, 124.35738100255621, 21.400152285708067, 19.64085404725156, 121.12299237559941, 68.20313316995285, 46.19435799477614, 15.134919452917416, 19.886781820820794, 27.162238637196495, 14.214072938293056, 115.14906277809592, 7.045879585977212, 142.95353742968172, 31.60500151064182, 16.67321112133154, 11.45311682700216, 8.222604877953065, 16.61028887905804, 21.49238866585931, 10.323817461630282, 59.707043868910574, 106.6933307669526, 42.45785193152026, 10.25874787538329, 23.753515835725985, 15.842270834389836, 54.247421554061184, 11.718015915453655, 9.251826900667709, 55.555421571855824, 5.182798693038229, 11.287845030306674, 21.609064889339955, 28.68418149825625, 58.57646450613494, 14.262417005490608, 16.612572203656686, 16.236087052979443, 7.07012245466578, 36.74755115014632, 5.56955315814885, 38.91744749703327, 78.85738379606225, 66.55885619672107, 39.64570920071233, 59.29291459700674, 43.05116230352876, 57.616583933068526, 13.018545934843285, 38.0480565816235, 24.925011809751684, 48.79848066987436, 48.274437799688435, 19.76507371158679, 8.565946758021937, 59.157588927921495, 57.14227355593507, 9.275322842239122, 14.083832907023066, 14.365127431573358, 16.940142017461092, 33.65000424197482, 22.38482569043266, 9.8423764229321, 103.30964749356153, 93.482737580894, 70.28695282389145, 36.95272628687471, 13.763226292274851, 66.32050179751066, 105.52022887233039, 20.58372196333329, 57.306558616504354, 7.2055606193971675, 6.901695760628813, 29.114219935320538, 5.196310986913392, 64.93439841265513, 66.54164242492251, 8.601749250514457, 16.577497833620654, 21.034827730068592, 71.15904863888585, 71.14756778662918, 83.33258945489067, 40.09744323758643, 10.758686591572488, 32.71729187462118, 72.60928034296013, 23.182448164356757, 20.02688139994615, 35.64640238564675, 19.619730168894044, 10.387210959701576, 108.49878733131999, 75.39827125131204, 8.518593311533348, 50.09760170113625, 5.783556515088281, 7.7226393825025, 17.362973213161744, 80.6800037352902, 45.02780319648969, 6.515288961328326, 21.02737688071963, 54.36451013450317, 57.52928577866641, 24.29359490456843, 49.28102430125691, 35.8859259513479, 8.56080309083324, 12.854971881557413, 8.88725896384838, 21.0800477677967, 98.2668911820237, 14.344215517868456, 20.26298700390288, 7.314621461467949, 54.44047864778453, 45.24597043884262, 24.081153968623944, 11.605888653834558, 81.92891659436373, 11.077075817245964, 72.47645741895147, 14.150248660347113, 38.98791894640125, 8.310910804301987, 14.562595175466456, 113.34002985790485, 45.54586004938276, 8.255063789030922, 15.748738231301685, 12.334676676031608, 37.13449911745443, 25.16227956950337, 76.85488255560675, 65.08240181335337, 13.072862840680108, 6.511161887246643, 36.836853509896514, 32.28105111670073, 9.011200476779752, 67.99261973751265, 43.56414154407412, 22.664659544292277, 10.641535346023423, 34.64376026085962, 34.85294630448037, 13.484231145440194, 28.36688965478537, 58.00923280980163, 35.763347543902626, 10.676355788751687, 21.408481350254558, 6.935108643188424, 11.81394661943359, 22.301879646778744, 5.798056980720284, 94.08282554710377, 22.67539113883788, 22.799743419864125, 11.478324417670079, 53.32099708503726, 30.313557087185988, 19.386367771690736, 18.844755698317343, 44.82455210759767, 15.373962501989634, 9.515152632133866, 16.939112837300648, 28.396269204357306, 10.335418544876822, 89.17929630357061, 16.840582885523023, 14.764919650394848, 19.293637630127673, 10.818970722833296, 17.912239148805387, 22.802790503154316, 23.11994332941388, 24.98230979382431, 24.23010562508092, 10.571739601893785, 13.638015383570616, 8.450786848419284, 5.422272727496316, 30.803919297691557, 6.16787773870539, 58.70960199565301, 21.757659586333236, 6.354290042601948, 79.01645583162893, 33.34361337283414, 10.43402929171163, 18.463413928011903, 51.398250308694465, 14.286148944932568, 39.54018858552626, 14.718239656426421, 7.133358728257916, 9.039139006726538, 101.67541212273028, 22.554844581646957, 81.56092029014815, 98.24038281892368, 91.74922363819401, 7.3638018403555074, 49.324132034224064, 52.30603773350151, 12.546541696075952, 55.03255187031403, 5.110580889802128, 13.547640273474753, 5.978462608040431, 140.66736056961003, 6.8420134193939175, 10.315950757587592, 47.18125856105616, 5.635843547192914, 6.538968180383621, 21.88992975546214, 28.648955299399073, 26.749893324709912, 69.50432749701974, 51.68692180174859, 12.564852483727488, 10.081517978439976, 84.29318403235067, 14.031376121318083, 8.789471837751982, 6.70667215533237, 44.29845720927677, 7.67457766137185, 80.51883119795062, 5.404723250860584, 22.0490516942453, 12.716189178806728, 33.87083847411847, 28.595637073850124, 5.167546674236084, 8.480139111408706, 24.72735125141671, 44.190865180027814, 15.528885568827413, 11.513728945059556, 34.379631842273476, 17.764093391162096, 110.36408081759643, 10.112704532840677, 5.6086556117359905, 14.575191366941128, 9.410796500229058, 45.32471755036951, 5.245340531125099, 6.464627593235395, 41.34674346917381, 25.368012276144974, 26.56679033128959, 78.43609293008124, 101.68991413483468, 88.73750329161878, 19.75121208094216, 77.54669677668076, 21.242898363177446, 26.471433912750108, 51.41048584252594, 88.97567150220965])
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)