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 = 44880
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);
([3289616.3830499407, 3303694.5670141024, 3303706.9044706216, 3303708.62705532, 3361407.9583998397, 3372267.817054787, 3401464.450281572, 3494554.9672297654, 3544403.9626828395, 3562595.7390554203, 3588408.0159743354, 3640633.9512020834, 3653202.4190230253, 3665428.678084648, 3725898.413753792, 3745536.409961781, 3754319.835869735, 3768993.058757255, 3784231.6581447837, 3813170.791731474, 3818703.125, 3819553.685008866, 3846665.6453555063, 3868875.0, 3915306.232238509, 3924389.3778110305, 3933269.0407879357, 3948245.09856601, 4017153.2948206244, 4065865.625, 4131079.1985916193, 4134002.792875643, 4134054.4819827727, 4146324.835310962, 4146927.7973075043, 4157772.8509146525, 4157840.2121202634, 4157855.086712655, 4166128.7508831043, 4166895.778472287, 4168572.673913555, 4175763.4946510848, 4179744.7832615515, 4181995.0015037656, 4182065.6993583217, 4182268.2224579006, 4186079.4393615285, 4188457.455351227, 4188625.4633132075, 4189796.4602017277, 4198432.226202912, 4202405.9591939, 4206865.955835955, 4209252.923257171, 4210675.0, 4221150.2505495185, 4225454.900116744, 4227334.518139776, 4238702.260954452, 4239033.334343837, 4243979.72942757, 4246993.925483671, 4257229.6875, 4265934.469385223, 4273098.4375, 4281392.1875, 4291977.451712106, 4301885.9375, 4302605.221500269, 4303754.6875, 4305334.140882799, 4305680.650675773, 4306107.782168999, 4315227.670876264, 4315253.105989039, 4320772.532055834, 4321144.163587471, 4322505.8720118245, 4329816.825886039, 4332003.125, 4332479.6875, 4336782.350282249, 4339399.574673927, 4339697.638694021, 4346062.441295115, 4353772.262266037, 4367696.875, 4371000.0, 4383442.689193379, 4401287.375764157, 4401634.8334376495, 4402617.1875, 4407311.163384037, 4414490.065961275, 4417180.672459226, 4421170.3125, 4433593.75, 4433596.875, 4437880.498377094, 4442429.9356823955, 4458206.25, 4464609.514155789, 4466451.352186747, 4466484.74530662, 4477463.670471125, 4488474.663080906, 4543459.375, 4575419.491041033, 4591650.0, 4603477.212171129, 4606817.421538717, 4633991.718609705, 4687024.084459776, 4785099.650177267, 4793053.330363577, 4799740.625, 4811941.749209431, 4841086.3524588235, 4862875.89835738, 4886606.25, 4898007.515825558, 4904237.180384333, 4919047.558897727, 4941061.7465640595, 4947778.125, 4949101.992831532, 5116764.0625, 5421576.563424182, 5505224.920949464, 5547839.688672303, 5547849.769453231, 5550546.61246778, 5571444.081099645, 5573311.371803439, 5604875.6817831155, 5661937.123181856, 5673700.0, 5681827.608924588, 5706551.548781086, 5735443.876131765, 5742395.281312907, 5754239.030822766, 5754272.028870332, 5754593.632696086, 5755321.5437598955, 5796128.113457928, 5800395.3125, 5806799.823324698, 5813269.846477928, 5815070.248028967, 5815090.812571533, 5817348.587655403, 5817540.939416457, 5817585.395797357, 5819959.656112894, 5821922.787940805, 5833321.266087948, 5835439.384940442, 5838804.109323221, 5839901.042427323, 5841887.850037032, 5846614.249342886, 5849080.505958405, 5851749.98352387, 5855143.27607705, 5855253.50196926, 5856360.9375, 5856904.6875, 5857206.460665059, 5858657.8125, 5859129.6875, 5859493.6399430055, 5860915.625, 5866837.5, 5870728.023045324, 5874533.955911645, 5877444.229755689, 5877829.710389734, 5886926.587098153, 5887062.499905393, 5887134.58187558, 5889012.860171742, 5890648.334333355, 5892566.816433423, 5892569.078513558, 5893661.804914798, 5895718.75, 5900498.184819494, 5902919.586696376, 5903920.215226034, 5904478.392413667, 5905178.75964319, 5907528.125, 5909263.838586813, 5914055.060733484, 5921007.687689695, 5927692.883480396, 5928389.11303896, 5930898.439695728, 5931655.298718919, 5932479.6812824905, 5933501.483703938, 5933837.485952549, 5934342.533183626, 5934454.419233422, 5934970.43986909, 5935554.352698775, 5936195.697668868, 5937976.5625, 5937992.8899850845, 5939317.295001825, 5940470.947406035, 5941362.043840735, 5949695.3125, 5949924.898704369, 5950548.4375, 5950573.622817895, 5951059.375, 5951839.925213284, 5952471.551872035, 5954026.5625, 5958367.417918985, 5959164.9684796315, 5962646.113974172, 5962667.1875, 5965065.144417244, 5966992.709973228, 5966997.6062842505, 5967020.779085337, 5967049.810437342, 5967459.474043333, 5969865.114193529, 5973260.702807805, 5973984.375, 5974008.134318849, 5974076.774725979, 5974688.015285756, 5977938.830444098, 5980260.476401447, 5981235.0268454375, 5984757.047407948, 5985282.720587147, 5987693.937852148, 5987889.079851043, 5988233.575709202, 5993968.75, 5997159.246196752, 6009191.668260643, 6012345.053799067, 6014705.034889674, 6017794.687987047, 6017860.954466893, 6017869.868230353, 6018362.7187450025, 6025099.468489863, 6026748.25129162, 6029698.4375, 6033835.9375, 6034915.446732479, 6034948.193695159, 6035386.325138181, 6035445.770275978, 6036271.875, 6037090.625395155, 6037853.30064235, 6047157.951857143, 6049715.686131479, 6049960.023042742, 6050001.276796241, 6050111.42252504, 6050790.806351585, 6057817.1875, 6059178.125, 6059725.125933288, 6060762.5, 6062453.125, 6065739.0625, 6066374.688657359, 6066985.587188809, 6067050.219859117, 6067504.687435813, 6072065.042190622, 6076798.10204158, 6077040.625, 6081055.656353732, 6086703.125, 6087114.767780678, 6101069.2181322565, 6108098.521106405, 6112348.4375, 6114194.397461384, 6114850.974775548, 6123734.375, 6124299.688264167, 6124912.197712646, 6125829.070609204, 6128445.684862198, 6130656.201244765, 6135457.646843425, 6135724.49364819, 6135728.524719748, 6145148.108845356, 6147923.205752619, 6154363.96950454, 6158850.349198318, 6161921.565706026, 6166757.8125, 6167287.589179943, 6169733.13435848, 6170772.219290323, 6176295.706944693, 6176800.548644639, 6182279.096597498, 6185676.385094928, 6185891.586901235, 6186636.291000028, 6198443.070768497, 6198462.208253725, 6199052.238999337, 6203021.860696968, 6208991.855663654, 6212173.4375, 6233754.715131111, 6240473.4375, 6244115.659093715, 6244608.750181072, 6245601.5625, 6247979.267473207, 6256820.878299428, 6277460.9375, 6317093.744322637, 6318363.602595327, 6319417.105744561, 6323637.5, 6325873.807591194, 6327563.469539941, 6327843.309945661, 6328268.75, 6328307.16036434, 6328345.585117513, 6329215.625, 6331150.167781317, 6331560.9375, 6332701.5625, 6332804.6875, 6333047.823056042, 6334517.783022789, 6336358.787123833, 6337334.375, 6358337.087313737, 6361166.965613411, 6367449.149153101, 6380437.5, 6386279.1977802785, 6386280.308141492, 6386446.043282702, 6389031.686949426, 6390645.009116986, 6396435.536584096, 6397412.420196039, 6399119.464237197, 6402526.505205462, 6405727.470812603, 6406661.198706667, 6406668.778281935, 6408417.1875, 6413128.723050547, 6417146.64863134, 6446437.659766821, 6458698.71588326, 6459535.929230323, 6466610.302083575, 6474251.303918442, 6476085.9375, 6476734.375, 6482265.67246656, 6482708.015784061, 6483517.394575945, 6484635.9375, 6487250.731318553, 6487326.727050012, 6495400.544489233, 6524210.6505563315, 6524670.027462907, 6525123.4375, 6526838.233645814, 6528319.564374161, 6529075.0, 6529465.949582848, 6529952.810811531, 6541050.938262046, 6542268.706984871, 6543000.0, 6545800.313993159, 6554682.692061474, 6559012.264259453, 6562657.983704123, 6574273.325909638, 6581576.163124048, 6583087.284877243, 6583459.316417675, 6595304.726133283, 6596133.328257039, 6597067.4528720025, 6599028.125, 6599416.319157749, 6605305.419027014, 6605988.4041949, 6611066.676840408, 6614205.0548423305, 6626606.25, 6630750.230010365, 6639432.834772261, 6640821.096792355, 6657221.724006354, 6662053.360775384, 6667753.125, 6672033.828855704, 6675252.853577672, 6677019.486909356, 6679177.36643968, 6679564.43942662, 6684409.884115663, 6692256.25, 6699323.4375, 6702493.75, 6716439.743684723, 6718960.9375, 6722282.8125, 6749125.269646801, 6750060.253922819, 6781578.16904868, 6809678.125, 6816781.724816517, 6832516.820785397, 6833057.484971444, 6834786.267557633, 6846960.063786482, 6854871.662572202, 6858669.5077494085, 6861126.078826982, 6863315.249637634, 6863844.309015883, 6864710.586290001, 6873079.6875, 6888504.681932202, 6894310.513748245, 6896878.125, 6923192.2570928, 6926181.754821068, 6941635.214056919, 6949385.9375, 6951900.0, 6955658.747382614, 6960945.3125, 6964371.636372421, 6990299.673243299, 6993711.680432656, 6993746.842152963, 7001763.023415416, 7007029.6875, 7022373.4375, 7025400.740033284, 7025469.608496787, 7035525.83280347, 7036709.708146781, 7041413.570514713, 7046545.271019533, 7048725.0, 7049251.491930821, 7051662.188468126, 7054182.026183254, 7056485.930265008, 7072083.888557982, 7088368.195021688, 7088818.75, 7090670.128102084, 7091463.600759135, 7092691.019895673, 7095576.984524515, 7099834.480588252, 7100436.707604002, 7101899.367351357, 7107732.70770487, 7112933.772735193, 7124803.57551533, 7129161.583280774, 7130523.10667298, 7130543.294759118, 7132214.454927503, 7132735.74705855, 7132758.789284774, 7133041.890095139, 7134105.801049691, 7134111.97947878, 7134697.081540903, 7134818.10311611, 7134880.697842841, 7135036.856000439, 7135576.953527398, 7137508.725500216, 7137577.6224639015, 7142183.393114338, 7142784.116201784, 7143368.090196706, 7146214.686929061, 7174771.839883506, 7179492.065261579, 7180857.8125, 7180989.437804634, 7181790.460382599, 7185202.455272951, 7189741.449947029, 7204196.7224793425, 7215190.802299676, 7224121.431291077, 7230879.370529711, 7249851.3870703215, 7254886.294253225, 7256880.608094226, 7260231.825628557, 7260452.023194743, 7261393.035001027, 7264263.1671164045, 7264314.250142559, 7268682.49685133, 7271042.397571985, 7272323.304372533, 7278461.782847884, 7285739.0625, 7287165.625, 7287970.307663501, 7289515.625, 7289587.5, 7290125.0, 7292579.887027521, 7295330.145272437, 7296136.186589512, 7297409.13124628, 7317139.229395409, 7319046.592480765, 7321564.128949349, 7328741.619418339, 7368325.0, 7376058.434119968, 7394138.896245486, 7395987.3637844445, 7396326.5625, 7399131.963512219, 7399774.573643863, 7399799.29204587, 7399994.744605444, 7401266.887865392, 7403372.4459833875, 7442324.830215419, 7447318.6997430995, 7448412.872578683, 7448599.3626386, 7448830.125831067, 7459831.370481638, 7463375.0, 7469901.15574963, 7475951.749028853, 7476313.497937507, 7478014.0625, 7478490.625, 7481839.0625, 7490710.513538135, 7493799.336985714, 7512636.983155156, 7513703.125, 7560288.201323972, 7563457.8125, 7569170.410417212, 7572234.177227849, 7576099.685209399, 7607892.561496948, 7620949.852031229, 7683170.3125, 7684048.210021862, 7695679.405734168, 7699631.02236792, 7699697.68176256, 7702626.5625, 7722303.680115381, 7725431.825363969, 7727442.25782351, 7786434.375, 7790034.303200654, 7797412.5, 7798145.454415936, 7803120.3125, 7821280.128831979, 7841223.010341579, 7858688.966487167, 7899405.975533073, 7908079.61981421, 7923534.358954104, 7951457.524088077, 7966084.283818429, 7968991.738725209, 7969764.5095685925, 7970616.18734958, 8000632.373788593, 8000898.4375, 8011707.296970498, 8012768.75, 8043410.9375, 8053118.070887132, 8053808.697154209, 8056702.714123152, 8059010.9375, 8059679.676433637, 8078506.196328188, 8081007.1612793645, 8095090.247517933, 8101286.151340823, 8127726.691403712, 8160329.049428431, 8173951.317712264, 8196317.044492019, 8209026.172193555, 8210971.79687322, 8246708.46457428, 8262961.878054241, 8271823.657164918, 8287509.375, 8287601.5625, 8297528.125, 8341496.875, 8365681.120108673, 8365722.956547563, 8367198.715055553, 8368008.0408552, 8368169.2268717075, 8373332.5523564825, 8380617.961785056, 8382192.8950491985, 8382331.009871078, 8382349.853952357, 8383450.125555632, 8384451.122629687, 8390727.440719936, 8391429.595035817, 8391433.376010943, 8393553.470759666, 8393992.172360806, 8394164.0464347, 8403725.603455922, 8405034.193148145, 8412140.00640307, 8414426.08726184, 8419634.573110107, 8433933.333003258, 8450601.5625, 8452876.890527936, 8456736.748251423, 8470040.48705608, 8470296.875, 8472338.33143693, 8472833.006901639, 8474060.568921125, 8475544.818850242, 8476308.642575273, 8477449.966659853, 8481650.914613513, 8482790.406239, 8489335.377942413, 8494129.654625783, 8494137.8962669, 8497600.837970074, 8499541.473605614, 8504205.553811705, 8512641.274911946, 8514124.576563193, 8515378.125, 8530138.165465537, 8539426.5625, 8541384.529661704, 8552110.070445914, 8579324.681289528, 8581282.139916603, 8581331.78789907, 8582965.625, 8590296.494222675, 8590302.947508479, 8596026.5625, 8600251.5625, 8606281.25, 8632326.363060344, 8633648.4375, 8633669.845301637, 8635640.15352166, 8681954.41009322, 8697483.686836494, 8733350.0, 8757856.400109652, 8793468.813630393, 8805935.932094613, 8848099.307259345, 8848370.273104887, 8852517.1875, 8853406.25, 8870984.717912393, 8876682.943980174, 8876967.387249144, 8879285.768579409, 8879412.5, 8881633.676934514, 8882778.125, 8903796.875, 8904585.484266013, 8905578.029146688, 8909803.233810954, 8909851.089518474, 8910978.385799784, 8930865.984135864, 8943120.233602319, 8945735.47362893, 8945751.5625, 8948221.875, 8950684.375, 8950901.011943795, 8956461.562494345, 8961579.6875, 8964436.795387678, 8964688.98091937, 8965792.03532538, 8970058.72365756, 8986254.975925455, 8993529.897136461, 8994241.789232586, 9019545.249114357, 9019743.746494891, 9022372.74302696, 9026669.866807934, 9028458.494580518, 9030086.541759923, 9045727.17471792, 9048846.567883665, 9050650.388532156, 9053207.8125, 9065083.48379432, 9069584.886152979, 9083421.875, 9086801.452734517, 9089983.167441096, 9090667.83408997, 9135931.25, 9154788.540910957, 9246389.0625, 9262965.522973154, 9267862.838519286, 9291503.542262172, 9360618.653265452, 9400183.63610315, 9443471.770797372, 9451215.625, 9457009.924365016, 9460111.02363269, 9461836.018810347, 9469067.688450215, 9473146.711364785, 9502144.707199419, 9502153.051739577, 9528425.673556097, 9540988.88600209, 9658557.8125, 9734686.821164336, 9748865.625, 9766832.8125, 9782686.19552208, 9794991.135267347, 9796087.516187126, 9801835.9375, 9801857.267142529, 9805505.811231393, 9808943.176634561, 9813069.911690751, 9814023.379680492, 9814215.909012, 9815021.489402318, 9816020.3125, 9817443.461948879, 9824382.067596195, 9824570.3125, 9832667.833517255, 9834509.375, 9844351.5625, 9845117.124794083, 9848192.70278282, 9849404.509532757, 9849561.151365537, 9849703.706774313, 9851301.99159263, 9852740.601028964, 9853515.625, 9853698.976033585, 9853871.774844388, 9856639.546954235, 9857753.632811321, 9858404.852083176, 9858967.80744238, 9866631.01706779, 10200285.716260478, 10203940.876564685, 10203954.466497345, 10203971.875, 10203979.091645405, 10206241.962831002, 10206507.431677168, 10206936.18081582, 10206985.903690565, 10207082.8125, 10207179.6875, 10208039.107279541, 10208341.442292217, 10208350.506202934, 10208542.1875, 10208907.474298516, 10208964.530767387, 10209186.61808688, 10209360.87677859, 10209938.112223286, 10210013.031422365, 10210080.73375325, 10211784.794553448, 10213155.580725893, 10213496.558544423, 10213651.526987754, 10213967.547036096, 10214197.424953973, 10214393.595190335, 10214423.393919578, 10214466.9662926, 10215703.973774176, 10215760.061409608, 10215855.445624068, 10216113.747815834, 10216721.477298686, 10216828.798439866, 10217177.305935368, 10217441.950855317, 10217630.852304908, 10218371.830002315, 10218773.4375, 10220193.155714465, 10220224.508995816, 10220234.490792703, 10220881.25, 10221139.226992983, 10222200.663476087, 10224507.062049726, 10226428.584123164, 10228152.841636727, 10229960.789658021, 10231150.76025522, 10249387.257046646, 10259074.357178938, 10260149.032567868, 10263467.951217942, 10270410.9375, 10270416.373344054, 10273151.5625, 10276444.657358345, 10278823.4375, 10281074.277478717, 10284992.724302351, 71543668.75], [44.210644164147894, 27.846974509521036, 5.220422786062373, 5.3706304100837965, 30.802528362090428, 5.363967177980706, 41.18255161895853, 52.2081174609432, 51.69012048864629, 35.93237377173189, 42.86041630930448, 33.10098985656373, 38.29796372041915, 69.44594539646117, 16.149665451704976, 36.33183206792127, 13.390509199029628, 29.719384060336814, 45.40016640077884, 42.86974915662669, 44.55293979737853, 35.4894532611812, 27.720454175393403, 32.091051798830435, 5.068573459091698, 18.37151353835292, 57.47568794275436, 46.145265952271025, 15.630860145996897, 54.09951552716386, 21.60446324045682, 42.904903628683, 6.241577301293607, 11.764817155360968, 9.068883113577602, 29.883174748753426, 22.40854293693177, 11.81483609406427, 64.71026984217056, 42.0910118935903, 6.012136380726729, 30.102758320112994, 23.71449998967597, 60.97198260149402, 11.245957382949038, 65.75087197048313, 91.43340453464448, 8.106329089947515, 66.34150422829208, 22.387332594203922, 47.54997631575006, 86.16977856033377, 22.170654304356674, 25.092588603585877, 39.22409360503969, 92.14138904428381, 11.477139959364774, 20.753873450528804, 10.001175882867992, 57.11266261355941, 27.002208385685336, 67.97180569807406, 78.79409371505454, 18.881683816052277, 83.21544980167798, 44.75461636129157, 26.789189932791096, 40.74717867933141, 76.86474617647269, 61.79933660489296, 18.263199027117274, 10.304478938822841, 56.9472417332081, 36.631112412768694, 5.937751644759673, 43.44201486118911, 6.320905639273739, 22.803172474756543, 13.721735177748569, 40.508227433497254, 59.19515400731169, 36.382231658402155, 13.535321340106497, 39.94691406404384, 11.037170554341765, 16.65625306723172, 49.71196845998314, 108.88820821309281, 5.642415332488307, 10.938916393556935, 9.751651093268586, 52.53647276923455, 5.721286277003084, 5.189409482110752, 115.23603369249324, 65.78574774603757, 55.378352672984946, 34.906274976880844, 139.14706316954187, 15.980104997978856, 31.745598646882947, 13.146948483268272, 7.528901832343341, 32.226680848000726, 60.429772997049085, 23.29033501024202, 42.21305542510233, 24.136387987004895, 48.30179755727117, 54.52097846697564, 85.32389101625012, 159.90462779523733, 108.4442827088075, 68.16732582059338, 6.519772958822442, 47.150891023728036, 13.498291896837197, 11.844800315138638, 51.537446603088725, 88.29919048116699, 22.854335919009284, 7.491588904076591, 9.985389312382974, 20.004837797451174, 34.24851809768077, 27.431144083191185, 44.235827865223875, 9.332184745181168, 27.053145202409453, 18.021286163225973, 52.622352001260396, 10.740734242582041, 11.687928263704798, 7.24049669300364, 18.383577830650825, 81.15847897403148, 53.44922752254976, 33.617333176326525, 26.535482071531852, 6.507046948448335, 35.27709678446164, 39.735721401662914, 6.794619366257612, 32.24780746311441, 119.50434500685101, 72.83275584642422, 77.76701999225023, 28.654011556407667, 21.838879389782196, 34.34538451103127, 18.536446767622472, 6.793858808706219, 27.84371236596501, 8.911836949653253, 58.073703659177305, 20.120525360640436, 82.72407200754405, 16.331593418764292, 19.598327740278833, 87.28025107851832, 11.72277797237795, 8.387917777271184, 65.9514512620011, 11.102050842121374, 8.050910708446224, 48.17069176889395, 30.542284476137766, 100.79835136234841, 6.294186401995925, 36.976565960282706, 56.097859492228764, 10.959930464867906, 44.058364768863726, 72.35152166082298, 69.50166984403525, 18.602127058429893, 82.45905835565492, 106.28772679077855, 25.4062034786718, 6.088089378520371, 15.793358979485912, 16.51844564224598, 23.118472994375125, 22.883514414008562, 27.485456079277885, 29.579341634771176, 44.85806232838513, 14.058591450783448, 56.327585561219294, 10.978990819206693, 62.09955537904838, 12.28937640115151, 73.92961705683211, 52.71271071769997, 78.15277725674827, 121.28884918884194, 18.529041841460437, 22.629655183119343, 14.949885119146868, 23.13444102441966, 14.917864024865011, 18.38803349170303, 35.19150829409274, 11.71473735925323, 6.233889524815095, 42.01975445067001, 111.84436293460266, 5.635714008647698, 85.10599376329716, 28.635609867508315, 8.344256164693682, 18.076942392533507, 13.447740891730083, 67.66971193667446, 61.62834625753729, 29.6969768700478, 15.698256924915693, 31.362448813422656, 8.802363872827373, 27.68897303090911, 144.6837524292166, 22.808963573232948, 52.23874380814196, 67.53688550685831, 57.20125526090664, 51.39702285771252, 28.255098888530785, 27.786382385781593, 11.944749802571216, 25.27196033731165, 11.28612920427231, 22.45537709229452, 69.74724345510245, 42.32799280191649, 15.781817824723383, 57.03754021246817, 79.3346011648395, 23.08585676218783, 5.900231580201749, 5.651153183515341, 33.743693377089755, 10.987017291838338, 37.948924395450256, 9.346117679879173, 67.81888781804079, 74.21743506202517, 100.1782433626228, 5.439912977513531, 6.44728812011987, 11.71788629345039, 12.997497449435436, 94.52051016955855, 20.037090650223476, 33.613206608964944, 69.99106215583319, 25.321902967527684, 85.77415049461906, 39.453360941980655, 34.22149854226853, 53.211352909633746, 16.658167579684655, 5.140407748703948, 49.38174386936166, 8.286786301321563, 70.33231988932548, 6.791650750544605, 20.737834144708682, 18.16151466199128, 34.77763273274769, 5.176813548732988, 79.25918085317332, 38.03142172308502, 72.61339213368892, 24.90251215780132, 31.741901917326146, 83.34453558318724, 105.98409569219098, 23.606928389552966, 6.695528427196868, 26.71921087117419, 61.69825660042816, 63.590043111606604, 192.00083349886418, 31.96317286142182, 17.80341728818045, 31.751169230485942, 10.008609262335414, 128.4263087253355, 76.91065553398778, 32.18550864307337, 18.15774994447192, 53.17382644553258, 33.144570174884784, 37.84467778285767, 22.908055345219353, 142.2538225294796, 5.354845160751018, 17.961942437755273, 44.136180704916775, 82.88743317827695, 5.732678414261903, 72.58837936733187, 24.66356615350088, 34.530806558748836, 16.46350565448162, 13.235833375492856, 46.46984662069872, 5.10021549759303, 11.797869996364003, 6.834244725129301, 31.01565453142038, 9.842809784619199, 30.45220111249258, 30.55100661883002, 7.399555028271431, 16.997835241676178, 86.97599993353967, 61.97000533799286, 130.56083565126616, 40.820569205697204, 73.66173916848525, 138.34560812841892, 15.047882642566679, 55.36860262362, 102.68532834751502, 5.435572808224367, 51.83863037194866, 23.07806620306089, 25.897958376269997, 72.02287279261743, 19.6236304846526, 7.005586219211726, 11.002277928513996, 94.74842568371281, 7.74281225602333, 6.028533380313232, 22.112685542372233, 59.81304513313698, 20.098801560783638, 67.02394097068692, 42.698935906908744, 7.665069888285425, 32.74221539095845, 76.50258825225518, 64.3168028470581, 5.668738086946826, 20.96124466676502, 8.879310908759818, 74.17356499946771, 49.221346376723304, 5.601885656054741, 12.457948188839826, 51.63121273086832, 21.59626816404885, 29.045987854637175, 7.890816446394075, 39.88583941156479, 7.497075055720798, 8.085365535783508, 40.03898491772664, 18.604011788002854, 9.283422761519024, 12.92270949509748, 8.429712938823435, 15.05041388490439, 35.92021909367521, 6.585113692636059, 106.9576828236705, 11.38696590700949, 16.12054070604987, 42.034761131901135, 27.06318718251292, 36.624380191024045, 48.7160439801092, 51.451413288304565, 66.68661090555378, 21.94389277376897, 5.406449796565294, 107.35172135432403, 18.722833433116946, 6.796808441840091, 12.100713145284194, 13.950061446192853, 6.231845370032907, 34.074375621506455, 13.827663430751961, 37.48854967280618, 78.12907124318967, 16.34204937281738, 7.041239730499225, 5.7691920323519845, 32.78749079031385, 66.68091541050704, 7.720893733590304, 60.38101754481355, 13.571053476584124, 7.666355660969796, 26.497969802461107, 27.909465879552656, 10.807319070576611, 39.358219386930244, 70.60350674584946, 9.809221890831216, 75.86860452021843, 69.67399279712482, 6.910481710270753, 6.608223091382152, 12.91243364761048, 13.296946423964002, 72.12041954347082, 95.81114802476702, 15.897205516519364, 5.069301535907275, 68.19741733238823, 10.899620500534606, 20.29707636601043, 50.36399485573122, 36.544727755020666, 31.089975034393127, 51.88320816730719, 6.620960082306551, 16.60329826056348, 103.0130466978617, 51.677612230436395, 68.36184834658779, 108.92189378252965, 18.379577377455302, 83.95296202269797, 70.39482220895694, 8.434511155537102, 76.36968052372765, 30.721151849584732, 85.548763284069, 8.554711980866147, 6.4311020268700405, 71.17698204781308, 16.36885191314452, 12.40288550265065, 42.981381924969966, 6.9248410702123415, 8.953381285004948, 63.970669484443, 75.35321111940877, 36.965498334474475, 69.75690697715339, 113.83458708569849, 8.074519922884239, 30.835810952332757, 9.390949306939632, 11.417259048654538, 12.753091089224462, 70.01529198453608, 77.58476481748784, 14.258689394148014, 32.01306157779783, 69.86504262862046, 152.98945645962425, 5.570257562606189, 38.08495668852991, 12.057896018483525, 39.388620958892005, 29.639652371390405, 5.870843609409775, 41.33112425264589, 19.471718131828435, 113.14562860043229, 30.28685057018614, 33.759515037158664, 43.19155850871684, 126.98580088007857, 39.69923006516617, 17.08600227194007, 11.116416844554246, 8.951381780937968, 19.471593796149726, 39.79317666899776, 77.7614432325422, 42.14668212030068, 6.550829133739846, 11.858324185642934, 119.9066040758278, 19.000123957484696, 10.264522681544916, 130.40520747273678, 151.44998103190397, 165.15772611591757, 18.146878311818064, 7.487810908004603, 13.312356881497095, 5.138978616232775, 79.05825990317436, 220.1090118165187, 31.68594973928165, 13.468067934621027, 8.926914140296864, 15.7918703803895, 6.023316788380455, 10.36196221215797, 41.85958783986794, 5.478803957216898, 20.121027001195277, 63.26543130696367, 17.738124111935367, 19.706277112502534, 19.21503513222622, 38.77375264977928, 27.848397888371412, 60.06617583657638, 98.33604745253103, 16.595032312372005, 54.29437893298605, 43.98012563087998, 6.940189797563369, 25.837025479577797, 6.289861833089044, 6.419714962754574, 73.3843728778506, 5.801113005884095, 7.737575183650098, 72.98670491724789, 5.218339506760482, 20.69468482511607, 82.2307100105246, 6.597012377553708, 22.81397546037877, 10.667743261449564, 98.92222271502045, 8.217409990939213, 6.180652681595996, 68.84986478324159, 45.91825466895254, 11.120155449835211, 77.83786889974296, 57.485588909774506, 35.41370274852639, 15.760764565709696, 7.24722017950947, 15.984274366204746, 8.161749607964719, 15.618228153944655, 23.976632278741146, 11.233184554693864, 8.007397758580833, 32.10472599588564, 24.240127885489287, 82.05744596785105, 62.91635882610656, 37.320387680661106, 21.393476285731097, 22.26060080019255, 7.365793530003927, 64.55801220294916, 16.569304333473944, 14.576787632537119, 14.331588324632364, 48.70466216501629, 45.94493394182712, 5.667360793516271, 193.64722538200405, 9.435398188344804, 55.990889221244814, 9.009628455531075, 99.59417441427951, 248.89490881517486, 58.49349797509953, 60.00593587731011, 113.51263423702898, 68.13797309042805, 20.0077324586111, 106.09741697825359, 56.48700769485229, 18.51547926797228, 46.70328939388363, 117.65117821337749, 20.89872641512945, 7.493217642102088, 27.23958916291044, 25.908574700420953, 72.90226488612129, 54.5042673204575, 10.714110449652576, 60.87411767062735, 5.083866725066565, 53.229607834101415, 12.122483614857835, 12.164307978962817, 15.216499980202325, 46.172993298559625, 22.583757047089883, 62.613575454694306, 34.53071250730285, 145.6851051530556, 16.940777574703723, 7.455718247487381, 69.80790302190725, 36.94350200033803, 7.103772068591833, 38.32458585983661, 53.00645683369863, 26.842950755959173, 13.468563974407907, 11.906653908354128, 22.035188892050712, 81.35952842334021, 69.09125482061091, 6.409204107700083, 46.41560201435009, 38.86366936413045, 86.80871998096978, 231.3304920335484, 64.739500486279, 47.16959200263488, 39.92189922680159, 52.77367418676694, 83.37182640731871, 13.674604067421473, 15.824427472942867, 11.335352002008426, 44.23746601974849, 24.482587419760247, 44.863308102505336, 115.96024503114913, 5.040229801498758, 18.187678304663496, 10.2772970821902, 17.826940692978056, 52.82693640406076, 49.41140425754997, 41.97296080480667, 53.541201917598926, 7.087720714161999, 105.2664781207838, 16.116646690547043, 11.504440787395053, 79.73703577705115, 152.88396943269592, 7.909540662093617, 85.2706667388582, 24.546872932191018, 166.6844498902187, 6.788834349098631, 57.987954497553915, 81.0508993721119, 47.619455260297364, 17.634821348493546, 20.07047778707327, 14.88411817519722, 34.12683557540174, 101.93618134602555, 30.73027005453305, 6.394304293123174, 29.26588043764485, 5.470334382121479, 87.46537837562161, 36.09264715969766, 11.682306495250904, 12.685680079364975, 117.81749563342655, 44.250863102786, 19.461116570062824, 15.7381368489984, 13.700989685471963, 24.284219108385933, 44.58633695147721, 14.815937083369521, 34.43906477389972, 24.847056396466343, 26.383360617262873, 27.89048151628637, 11.809256994446155, 25.63565834795466, 19.595398949195967, 19.73640518150133, 15.193605505002047, 5.591040367122229, 86.1936680047944, 8.672903026843763, 40.02044422039164, 15.56570395327665, 18.475542416106414, 7.040118982673326, 12.876991463972603, 17.722162933238252, 69.28308449362856, 10.580484888075189, 26.722218142078496, 66.7138510916618, 46.5903978171471, 44.79068415460786, 8.183852960310698, 35.33576995824221, 86.73880655156367, 21.145190721057183, 110.44883250182596, 7.3719612906533225, 61.793933543379644, 37.92146227214558, 31.164118324241326, 11.119268333246968, 41.952819094474265, 46.9135566306072, 62.41430657897161, 68.19881229845525, 18.37940514121238, 9.445107798993, 26.896431537907926, 13.685798801704255, 36.34766209720343, 12.814461360682499, 59.742171268612466, 44.397879166051474, 28.26505195740198, 34.67525458571081, 7.643072630024772, 18.83253472872454, 16.038367433587386, 34.48297328842959, 45.51842166208217, 8.967730997798599, 57.65891998864141, 52.23216690684241, 41.53100716119738, 65.4001631657973, 7.846091723473722, 45.8406040195537, 6.017767669667718, 62.286296346126136, 25.842448020291695, 12.929823563688961, 14.48339070410606, 43.3603060183674, 20.35384217613722, 105.60504291941655, 76.90489830765591, 30.929180533842707, 157.87295333665537, 12.392376553609214, 75.805977267255, 5.860711806082114, 53.60520120849155, 5.843298738935159, 61.448032454674824, 12.368722939737864, 17.32825607331559, 62.628001570840034, 7.564014044659709, 109.77055939144033, 77.34727958593783, 87.61115620454642, 83.78687905271478, 65.01005664397863, 5.338437708586707, 248.95259818850028, 303.87017565236636, 26.748026644483023, 111.16887142124719, 84.87159686486002, 58.02941151263842, 75.64934224268958, 15.280294487647078, 8.330048786082735, 17.270155214471746, 25.664298684140483, 10.29727365162832, 7.102325968998719, 7.8664818829635905, 14.308910432711299, 83.76568580083402, 50.18201810232434, 59.15244427711543, 57.62508872574951, 63.13072273567729, 7.798843030457015, 42.238002553252336, 30.082026923529387, 21.052261276204003, 22.12646517535551, 13.765003702484783, 71.28992716923128, 7.582458370072916, 16.147090024041823, 10.574122143882033, 33.11196772016257, 22.369780309119513, 6.938501919480215, 65.90338956437972, 12.313831998677287, 48.517692697396974, 46.30659029611051, 8.253693432901054, 17.351922640560396, 10.860478297805452, 46.42816370119364, 8.620213805734736, 18.686151740465235, 6.758730354333417, 38.20753770521827, 12.088564582721174, 63.87826164965371, 9.777630956207998, 47.30639260284667, 9.47455585814036, 45.04249546431757, 38.62701258247531, 24.81178093509485, 38.972851122717074, 27.12532690816741, 52.64818331909333, 5.8891419706641415, 45.499189224675874, 5.514879119892801, 42.59221923638069, 11.075268526991017, 68.71888508661057, 34.94309316907344, 7.1332727986129925, 5.660917480754918, 97.53049157676335, 54.338705722020975, 13.79214974954378, 52.58119517160878, 52.84459384944259, 6.081807405036744, 59.74401388839681, 8.589126295423446, 5.586108584093594, 38.70532161406103, 7.377871654444003, 6.7019314631375435, 13.158141454921738, 7.219282643617524, 12.111893881070241, 59.18549410714026, 9.294823328470022, 6.452965252628435, 19.596729154146722, 5.26477138194249, 15.366741294527094, 57.21238397558369, 17.040929952848504, 12.373416945880273, 5.755764221903059, 42.28719645501298, 7.4654643087068155, 59.90143665698936, 89.83044399810169, 13.075030509901751, 8.042940643115244, 5.168067121846161, 41.21378411500444, 34.336648571056756, 7.307510794870824, 50.51950491522598, 5.848750855081385, 35.28328692281315, 54.9030741283749, 49.244725560905295, 38.220896723498385, 66.76695033044527, 12.205440893687532, 83.53019688278354, 38.797156124043106, 7.33311917442329, 36.98268087723461, 30.96196745763318, 41.02760466045344, 19.527628198511152, 6.337488447825835, 38.04608050975831])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([3289616.3830499407, 3303694.5670141024, 3303706.9044706216, 3303708.62705532, 3361407.9583998397, 3372267.817054787, 3401464.450281572, 3494554.9672297654, 3544403.9626828395, 3562595.7390554203, 3588408.0159743354, 3640633.9512020834, 3653202.4190230253, 3665428.678084648, 3725898.413753792, 3745536.409961781, 3754319.835869735, 3768993.058757255, 3784231.6581447837, 3813170.791731474, 3818703.125, 3819553.685008866, 3846665.6453555063, 3868875.0, 3915306.232238509, 3924389.3778110305, 3933269.0407879357, 3948245.09856601, 4017153.2948206244, 4065865.625, 4131079.1985916193, 4134002.792875643, 4134054.4819827727, 4146324.835310962, 4146927.7973075043, 4157772.8509146525, 4157840.2121202634, 4157855.086712655, 4166128.7508831043, 4166895.778472287, 4168572.673913555, 4175763.4946510848, 4179744.7832615515, 4181995.0015037656, 4182065.6993583217, 4182268.2224579006, 4186079.4393615285, 4188457.455351227, 4188625.4633132075, 4189796.4602017277, 4198432.226202912, 4202405.9591939, 4206865.955835955, 4209252.923257171, 4210675.0, 4221150.2505495185, 4225454.900116744, 4227334.518139776, 4238702.260954452, 4239033.334343837, 4243979.72942757, 4246993.925483671, 4257229.6875, 4265934.469385223, 4273098.4375, 4281392.1875, 4291977.451712106, 4301885.9375, 4302605.221500269, 4303754.6875, 4305334.140882799, 4305680.650675773, 4306107.782168999, 4315227.670876264, 4315253.105989039, 4320772.532055834, 4321144.163587471, 4322505.8720118245, 4329816.825886039, 4332003.125, 4332479.6875, 4336782.350282249, 4339399.574673927, 4339697.638694021, 4346062.441295115, 4353772.262266037, 4367696.875, 4371000.0, 4383442.689193379, 4401287.375764157, 4401634.8334376495, 4402617.1875, 4407311.163384037, 4414490.065961275, 4417180.672459226, 4421170.3125, 4433593.75, 4433596.875, 4437880.498377094, 4442429.9356823955, 4458206.25, 4464609.514155789, 4466451.352186747, 4466484.74530662, 4477463.670471125, 4488474.663080906, 4543459.375, 4575419.491041033, 4591650.0, 4603477.212171129, 4606817.421538717, 4633991.718609705, 4687024.084459776, 4785099.650177267, 4793053.330363577, 4799740.625, 4811941.749209431, 4841086.3524588235, 4862875.89835738, 4886606.25, 4898007.515825558, 4904237.180384333, 4919047.558897727, 4941061.7465640595, 4947778.125, 4949101.992831532, 5116764.0625, 5421576.563424182, 5505224.920949464, 5547839.688672303, 5547849.769453231, 5550546.61246778, 5571444.081099645, 5573311.371803439, 5604875.6817831155, 5661937.123181856, 5673700.0, 5681827.608924588, 5706551.548781086, 5735443.876131765, 5742395.281312907, 5754239.030822766, 5754272.028870332, 5754593.632696086, 5755321.5437598955, 5796128.113457928, 5800395.3125, 5806799.823324698, 5813269.846477928, 5815070.248028967, 5815090.812571533, 5817348.587655403, 5817540.939416457, 5817585.395797357, 5819959.656112894, 5821922.787940805, 5833321.266087948, 5835439.384940442, 5838804.109323221, 5839901.042427323, 5841887.850037032, 5846614.249342886, 5849080.505958405, 5851749.98352387, 5855143.27607705, 5855253.50196926, 5856360.9375, 5856904.6875, 5857206.460665059, 5858657.8125, 5859129.6875, 5859493.6399430055, 5860915.625, 5866837.5, 5870728.023045324, 5874533.955911645, 5877444.229755689, 5877829.710389734, 5886926.587098153, 5887062.499905393, 5887134.58187558, 5889012.860171742, 5890648.334333355, 5892566.816433423, 5892569.078513558, 5893661.804914798, 5895718.75, 5900498.184819494, 5902919.586696376, 5903920.215226034, 5904478.392413667, 5905178.75964319, 5907528.125, 5909263.838586813, 5914055.060733484, 5921007.687689695, 5927692.883480396, 5928389.11303896, 5930898.439695728, 5931655.298718919, 5932479.6812824905, 5933501.483703938, 5933837.485952549, 5934342.533183626, 5934454.419233422, 5934970.43986909, 5935554.352698775, 5936195.697668868, 5937976.5625, 5937992.8899850845, 5939317.295001825, 5940470.947406035, 5941362.043840735, 5949695.3125, 5949924.898704369, 5950548.4375, 5950573.622817895, 5951059.375, 5951839.925213284, 5952471.551872035, 5954026.5625, 5958367.417918985, 5959164.9684796315, 5962646.113974172, 5962667.1875, 5965065.144417244, 5966992.709973228, 5966997.6062842505, 5967020.779085337, 5967049.810437342, 5967459.474043333, 5969865.114193529, 5973260.702807805, 5973984.375, 5974008.134318849, 5974076.774725979, 5974688.015285756, 5977938.830444098, 5980260.476401447, 5981235.0268454375, 5984757.047407948, 5985282.720587147, 5987693.937852148, 5987889.079851043, 5988233.575709202, 5993968.75, 5997159.246196752, 6009191.668260643, 6012345.053799067, 6014705.034889674, 6017794.687987047, 6017860.954466893, 6017869.868230353, 6018362.7187450025, 6025099.468489863, 6026748.25129162, 6029698.4375, 6033835.9375, 6034915.446732479, 6034948.193695159, 6035386.325138181, 6035445.770275978, 6036271.875, 6037090.625395155, 6037853.30064235, 6047157.951857143, 6049715.686131479, 6049960.023042742, 6050001.276796241, 6050111.42252504, 6050790.806351585, 6057817.1875, 6059178.125, 6059725.125933288, 6060762.5, 6062453.125, 6065739.0625, 6066374.688657359, 6066985.587188809, 6067050.219859117, 6067504.687435813, 6072065.042190622, 6076798.10204158, 6077040.625, 6081055.656353732, 6086703.125, 6087114.767780678, 6101069.2181322565, 6108098.521106405, 6112348.4375, 6114194.397461384, 6114850.974775548, 6123734.375, 6124299.688264167, 6124912.197712646, 6125829.070609204, 6128445.684862198, 6130656.201244765, 6135457.646843425, 6135724.49364819, 6135728.524719748, 6145148.108845356, 6147923.205752619, 6154363.96950454, 6158850.349198318, 6161921.565706026, 6166757.8125, 6167287.589179943, 6169733.13435848, 6170772.219290323, 6176295.706944693, 6176800.548644639, 6182279.096597498, 6185676.385094928, 6185891.586901235, 6186636.291000028, 6198443.070768497, 6198462.208253725, 6199052.238999337, 6203021.860696968, 6208991.855663654, 6212173.4375, 6233754.715131111, 6240473.4375, 6244115.659093715, 6244608.750181072, 6245601.5625, 6247979.267473207, 6256820.878299428, 6277460.9375, 6317093.744322637, 6318363.602595327, 6319417.105744561, 6323637.5, 6325873.807591194, 6327563.469539941, 6327843.309945661, 6328268.75, 6328307.16036434, 6328345.585117513, 6329215.625, 6331150.167781317, 6331560.9375, 6332701.5625, 6332804.6875, 6333047.823056042, 6334517.783022789, 6336358.787123833, 6337334.375, 6358337.087313737, 6361166.965613411, 6367449.149153101, 6380437.5, 6386279.1977802785, 6386280.308141492, 6386446.043282702, 6389031.686949426, 6390645.009116986, 6396435.536584096, 6397412.420196039, 6399119.464237197, 6402526.505205462, 6405727.470812603, 6406661.198706667, 6406668.778281935, 6408417.1875, 6413128.723050547, 6417146.64863134, 6446437.659766821, 6458698.71588326, 6459535.929230323, 6466610.302083575, 6474251.303918442, 6476085.9375, 6476734.375, 6482265.67246656, 6482708.015784061, 6483517.394575945, 6484635.9375, 6487250.731318553, 6487326.727050012, 6495400.544489233, 6524210.6505563315, 6524670.027462907, 6525123.4375, 6526838.233645814, 6528319.564374161, 6529075.0, 6529465.949582848, 6529952.810811531, 6541050.938262046, 6542268.706984871, 6543000.0, 6545800.313993159, 6554682.692061474, 6559012.264259453, 6562657.983704123, 6574273.325909638, 6581576.163124048, 6583087.284877243, 6583459.316417675, 6595304.726133283, 6596133.328257039, 6597067.4528720025, 6599028.125, 6599416.319157749, 6605305.419027014, 6605988.4041949, 6611066.676840408, 6614205.0548423305, 6626606.25, 6630750.230010365, 6639432.834772261, 6640821.096792355, 6657221.724006354, 6662053.360775384, 6667753.125, 6672033.828855704, 6675252.853577672, 6677019.486909356, 6679177.36643968, 6679564.43942662, 6684409.884115663, 6692256.25, 6699323.4375, 6702493.75, 6716439.743684723, 6718960.9375, 6722282.8125, 6749125.269646801, 6750060.253922819, 6781578.16904868, 6809678.125, 6816781.724816517, 6832516.820785397, 6833057.484971444, 6834786.267557633, 6846960.063786482, 6854871.662572202, 6858669.5077494085, 6861126.078826982, 6863315.249637634, 6863844.309015883, 6864710.586290001, 6873079.6875, 6888504.681932202, 6894310.513748245, 6896878.125, 6923192.2570928, 6926181.754821068, 6941635.214056919, 6949385.9375, 6951900.0, 6955658.747382614, 6960945.3125, 6964371.636372421, 6990299.673243299, 6993711.680432656, 6993746.842152963, 7001763.023415416, 7007029.6875, 7022373.4375, 7025400.740033284, 7025469.608496787, 7035525.83280347, 7036709.708146781, 7041413.570514713, 7046545.271019533, 7048725.0, 7049251.491930821, 7051662.188468126, 7054182.026183254, 7056485.930265008, 7072083.888557982, 7088368.195021688, 7088818.75, 7090670.128102084, 7091463.600759135, 7092691.019895673, 7095576.984524515, 7099834.480588252, 7100436.707604002, 7101899.367351357, 7107732.70770487, 7112933.772735193, 7124803.57551533, 7129161.583280774, 7130523.10667298, 7130543.294759118, 7132214.454927503, 7132735.74705855, 7132758.789284774, 7133041.890095139, 7134105.801049691, 7134111.97947878, 7134697.081540903, 7134818.10311611, 7134880.697842841, 7135036.856000439, 7135576.953527398, 7137508.725500216, 7137577.6224639015, 7142183.393114338, 7142784.116201784, 7143368.090196706, 7146214.686929061, 7174771.839883506, 7179492.065261579, 7180857.8125, 7180989.437804634, 7181790.460382599, 7185202.455272951, 7189741.449947029, 7204196.7224793425, 7215190.802299676, 7224121.431291077, 7230879.370529711, 7249851.3870703215, 7254886.294253225, 7256880.608094226, 7260231.825628557, 7260452.023194743, 7261393.035001027, 7264263.1671164045, 7264314.250142559, 7268682.49685133, 7271042.397571985, 7272323.304372533, 7278461.782847884, 7285739.0625, 7287165.625, 7287970.307663501, 7289515.625, 7289587.5, 7290125.0, 7292579.887027521, 7295330.145272437, 7296136.186589512, 7297409.13124628, 7317139.229395409, 7319046.592480765, 7321564.128949349, 7328741.619418339, 7368325.0, 7376058.434119968, 7394138.896245486, 7395987.3637844445, 7396326.5625, 7399131.963512219, 7399774.573643863, 7399799.29204587, 7399994.744605444, 7401266.887865392, 7403372.4459833875, 7442324.830215419, 7447318.6997430995, 7448412.872578683, 7448599.3626386, 7448830.125831067, 7459831.370481638, 7463375.0, 7469901.15574963, 7475951.749028853, 7476313.497937507, 7478014.0625, 7478490.625, 7481839.0625, 7490710.513538135, 7493799.336985714, 7512636.983155156, 7513703.125, 7560288.201323972, 7563457.8125, 7569170.410417212, 7572234.177227849, 7576099.685209399, 7607892.561496948, 7620949.852031229, 7683170.3125, 7684048.210021862, 7695679.405734168, 7699631.02236792, 7699697.68176256, 7702626.5625, 7722303.680115381, 7725431.825363969, 7727442.25782351, 7786434.375, 7790034.303200654, 7797412.5, 7798145.454415936, 7803120.3125, 7821280.128831979, 7841223.010341579, 7858688.966487167, 7899405.975533073, 7908079.61981421, 7923534.358954104, 7951457.524088077, 7966084.283818429, 7968991.738725209, 7969764.5095685925, 7970616.18734958, 8000632.373788593, 8000898.4375, 8011707.296970498, 8012768.75, 8043410.9375, 8053118.070887132, 8053808.697154209, 8056702.714123152, 8059010.9375, 8059679.676433637, 8078506.196328188, 8081007.1612793645, 8095090.247517933, 8101286.151340823, 8127726.691403712, 8160329.049428431, 8173951.317712264, 8196317.044492019, 8209026.172193555, 8210971.79687322, 8246708.46457428, 8262961.878054241, 8271823.657164918, 8287509.375, 8287601.5625, 8297528.125, 8341496.875, 8365681.120108673, 8365722.956547563, 8367198.715055553, 8368008.0408552, 8368169.2268717075, 8373332.5523564825, 8380617.961785056, 8382192.8950491985, 8382331.009871078, 8382349.853952357, 8383450.125555632, 8384451.122629687, 8390727.440719936, 8391429.595035817, 8391433.376010943, 8393553.470759666, 8393992.172360806, 8394164.0464347, 8403725.603455922, 8405034.193148145, 8412140.00640307, 8414426.08726184, 8419634.573110107, 8433933.333003258, 8450601.5625, 8452876.890527936, 8456736.748251423, 8470040.48705608, 8470296.875, 8472338.33143693, 8472833.006901639, 8474060.568921125, 8475544.818850242, 8476308.642575273, 8477449.966659853, 8481650.914613513, 8482790.406239, 8489335.377942413, 8494129.654625783, 8494137.8962669, 8497600.837970074, 8499541.473605614, 8504205.553811705, 8512641.274911946, 8514124.576563193, 8515378.125, 8530138.165465537, 8539426.5625, 8541384.529661704, 8552110.070445914, 8579324.681289528, 8581282.139916603, 8581331.78789907, 8582965.625, 8590296.494222675, 8590302.947508479, 8596026.5625, 8600251.5625, 8606281.25, 8632326.363060344, 8633648.4375, 8633669.845301637, 8635640.15352166, 8681954.41009322, 8697483.686836494, 8733350.0, 8757856.400109652, 8793468.813630393, 8805935.932094613, 8848099.307259345, 8848370.273104887, 8852517.1875, 8853406.25, 8870984.717912393, 8876682.943980174, 8876967.387249144, 8879285.768579409, 8879412.5, 8881633.676934514, 8882778.125, 8903796.875, 8904585.484266013, 8905578.029146688, 8909803.233810954, 8909851.089518474, 8910978.385799784, 8930865.984135864, 8943120.233602319, 8945735.47362893, 8945751.5625, 8948221.875, 8950684.375, 8950901.011943795, 8956461.562494345, 8961579.6875, 8964436.795387678, 8964688.98091937, 8965792.03532538, 8970058.72365756, 8986254.975925455, 8993529.897136461, 8994241.789232586, 9019545.249114357, 9019743.746494891, 9022372.74302696, 9026669.866807934, 9028458.494580518, 9030086.541759923, 9045727.17471792, 9048846.567883665, 9050650.388532156, 9053207.8125, 9065083.48379432, 9069584.886152979, 9083421.875, 9086801.452734517, 9089983.167441096, 9090667.83408997, 9135931.25, 9154788.540910957, 9246389.0625, 9262965.522973154, 9267862.838519286, 9291503.542262172, 9360618.653265452, 9400183.63610315, 9443471.770797372, 9451215.625, 9457009.924365016, 9460111.02363269, 9461836.018810347, 9469067.688450215, 9473146.711364785, 9502144.707199419, 9502153.051739577, 9528425.673556097, 9540988.88600209, 9658557.8125, 9734686.821164336, 9748865.625, 9766832.8125, 9782686.19552208, 9794991.135267347, 9796087.516187126, 9801835.9375, 9801857.267142529, 9805505.811231393, 9808943.176634561, 9813069.911690751, 9814023.379680492, 9814215.909012, 9815021.489402318, 9816020.3125, 9817443.461948879, 9824382.067596195, 9824570.3125, 9832667.833517255, 9834509.375, 9844351.5625, 9845117.124794083, 9848192.70278282, 9849404.509532757, 9849561.151365537, 9849703.706774313, 9851301.99159263, 9852740.601028964, 9853515.625, 9853698.976033585, 9853871.774844388, 9856639.546954235, 9857753.632811321, 9858404.852083176, 9858967.80744238, 9866631.01706779, 10200285.716260478, 10203940.876564685, 10203954.466497345, 10203971.875, 10203979.091645405, 10206241.962831002, 10206507.431677168, 10206936.18081582, 10206985.903690565, 10207082.8125, 10207179.6875, 10208039.107279541, 10208341.442292217, 10208350.506202934, 10208542.1875, 10208907.474298516, 10208964.530767387, 10209186.61808688, 10209360.87677859, 10209938.112223286, 10210013.031422365, 10210080.73375325, 10211784.794553448, 10213155.580725893, 10213496.558544423, 10213651.526987754, 10213967.547036096, 10214197.424953973, 10214393.595190335, 10214423.393919578, 10214466.9662926, 10215703.973774176, 10215760.061409608, 10215855.445624068, 10216113.747815834, 10216721.477298686, 10216828.798439866, 10217177.305935368, 10217441.950855317, 10217630.852304908, 10218371.830002315, 10218773.4375, 10220193.155714465, 10220224.508995816, 10220234.490792703, 10220881.25, 10221139.226992983, 10222200.663476087, 10224507.062049726, 10226428.584123164, 10228152.841636727, 10229960.789658021, 10231150.76025522, 10249387.257046646, 10259074.357178938, 10260149.032567868, 10263467.951217942, 10270410.9375, 10270416.373344054, 10273151.5625, 10276444.657358345, 10278823.4375, 10281074.277478717, 10284992.724302351, 71543668.75], [44.210644164147894, 27.846974509521036, 5.220422786062373, 5.3706304100837965, 30.802528362090428, 5.363967177980706, 41.18255161895853, 52.2081174609432, 51.69012048864629, 35.93237377173189, 42.86041630930448, 33.10098985656373, 38.29796372041915, 69.44594539646117, 16.149665451704976, 36.33183206792127, 13.390509199029628, 29.719384060336814, 45.40016640077884, 42.86974915662669, 44.55293979737853, 35.4894532611812, 27.720454175393403, 32.091051798830435, 5.068573459091698, 18.37151353835292, 57.47568794275436, 46.145265952271025, 15.630860145996897, 54.09951552716386, 21.60446324045682, 42.904903628683, 6.241577301293607, 11.764817155360968, 9.068883113577602, 29.883174748753426, 22.40854293693177, 11.81483609406427, 64.71026984217056, 42.0910118935903, 6.012136380726729, 30.102758320112994, 23.71449998967597, 60.97198260149402, 11.245957382949038, 65.75087197048313, 91.43340453464448, 8.106329089947515, 66.34150422829208, 22.387332594203922, 47.54997631575006, 86.16977856033377, 22.170654304356674, 25.092588603585877, 39.22409360503969, 92.14138904428381, 11.477139959364774, 20.753873450528804, 10.001175882867992, 57.11266261355941, 27.002208385685336, 67.97180569807406, 78.79409371505454, 18.881683816052277, 83.21544980167798, 44.75461636129157, 26.789189932791096, 40.74717867933141, 76.86474617647269, 61.79933660489296, 18.263199027117274, 10.304478938822841, 56.9472417332081, 36.631112412768694, 5.937751644759673, 43.44201486118911, 6.320905639273739, 22.803172474756543, 13.721735177748569, 40.508227433497254, 59.19515400731169, 36.382231658402155, 13.535321340106497, 39.94691406404384, 11.037170554341765, 16.65625306723172, 49.71196845998314, 108.88820821309281, 5.642415332488307, 10.938916393556935, 9.751651093268586, 52.53647276923455, 5.721286277003084, 5.189409482110752, 115.23603369249324, 65.78574774603757, 55.378352672984946, 34.906274976880844, 139.14706316954187, 15.980104997978856, 31.745598646882947, 13.146948483268272, 7.528901832343341, 32.226680848000726, 60.429772997049085, 23.29033501024202, 42.21305542510233, 24.136387987004895, 48.30179755727117, 54.52097846697564, 85.32389101625012, 159.90462779523733, 108.4442827088075, 68.16732582059338, 6.519772958822442, 47.150891023728036, 13.498291896837197, 11.844800315138638, 51.537446603088725, 88.29919048116699, 22.854335919009284, 7.491588904076591, 9.985389312382974, 20.004837797451174, 34.24851809768077, 27.431144083191185, 44.235827865223875, 9.332184745181168, 27.053145202409453, 18.021286163225973, 52.622352001260396, 10.740734242582041, 11.687928263704798, 7.24049669300364, 18.383577830650825, 81.15847897403148, 53.44922752254976, 33.617333176326525, 26.535482071531852, 6.507046948448335, 35.27709678446164, 39.735721401662914, 6.794619366257612, 32.24780746311441, 119.50434500685101, 72.83275584642422, 77.76701999225023, 28.654011556407667, 21.838879389782196, 34.34538451103127, 18.536446767622472, 6.793858808706219, 27.84371236596501, 8.911836949653253, 58.073703659177305, 20.120525360640436, 82.72407200754405, 16.331593418764292, 19.598327740278833, 87.28025107851832, 11.72277797237795, 8.387917777271184, 65.9514512620011, 11.102050842121374, 8.050910708446224, 48.17069176889395, 30.542284476137766, 100.79835136234841, 6.294186401995925, 36.976565960282706, 56.097859492228764, 10.959930464867906, 44.058364768863726, 72.35152166082298, 69.50166984403525, 18.602127058429893, 82.45905835565492, 106.28772679077855, 25.4062034786718, 6.088089378520371, 15.793358979485912, 16.51844564224598, 23.118472994375125, 22.883514414008562, 27.485456079277885, 29.579341634771176, 44.85806232838513, 14.058591450783448, 56.327585561219294, 10.978990819206693, 62.09955537904838, 12.28937640115151, 73.92961705683211, 52.71271071769997, 78.15277725674827, 121.28884918884194, 18.529041841460437, 22.629655183119343, 14.949885119146868, 23.13444102441966, 14.917864024865011, 18.38803349170303, 35.19150829409274, 11.71473735925323, 6.233889524815095, 42.01975445067001, 111.84436293460266, 5.635714008647698, 85.10599376329716, 28.635609867508315, 8.344256164693682, 18.076942392533507, 13.447740891730083, 67.66971193667446, 61.62834625753729, 29.6969768700478, 15.698256924915693, 31.362448813422656, 8.802363872827373, 27.68897303090911, 144.6837524292166, 22.808963573232948, 52.23874380814196, 67.53688550685831, 57.20125526090664, 51.39702285771252, 28.255098888530785, 27.786382385781593, 11.944749802571216, 25.27196033731165, 11.28612920427231, 22.45537709229452, 69.74724345510245, 42.32799280191649, 15.781817824723383, 57.03754021246817, 79.3346011648395, 23.08585676218783, 5.900231580201749, 5.651153183515341, 33.743693377089755, 10.987017291838338, 37.948924395450256, 9.346117679879173, 67.81888781804079, 74.21743506202517, 100.1782433626228, 5.439912977513531, 6.44728812011987, 11.71788629345039, 12.997497449435436, 94.52051016955855, 20.037090650223476, 33.613206608964944, 69.99106215583319, 25.321902967527684, 85.77415049461906, 39.453360941980655, 34.22149854226853, 53.211352909633746, 16.658167579684655, 5.140407748703948, 49.38174386936166, 8.286786301321563, 70.33231988932548, 6.791650750544605, 20.737834144708682, 18.16151466199128, 34.77763273274769, 5.176813548732988, 79.25918085317332, 38.03142172308502, 72.61339213368892, 24.90251215780132, 31.741901917326146, 83.34453558318724, 105.98409569219098, 23.606928389552966, 6.695528427196868, 26.71921087117419, 61.69825660042816, 63.590043111606604, 192.00083349886418, 31.96317286142182, 17.80341728818045, 31.751169230485942, 10.008609262335414, 128.4263087253355, 76.91065553398778, 32.18550864307337, 18.15774994447192, 53.17382644553258, 33.144570174884784, 37.84467778285767, 22.908055345219353, 142.2538225294796, 5.354845160751018, 17.961942437755273, 44.136180704916775, 82.88743317827695, 5.732678414261903, 72.58837936733187, 24.66356615350088, 34.530806558748836, 16.46350565448162, 13.235833375492856, 46.46984662069872, 5.10021549759303, 11.797869996364003, 6.834244725129301, 31.01565453142038, 9.842809784619199, 30.45220111249258, 30.55100661883002, 7.399555028271431, 16.997835241676178, 86.97599993353967, 61.97000533799286, 130.56083565126616, 40.820569205697204, 73.66173916848525, 138.34560812841892, 15.047882642566679, 55.36860262362, 102.68532834751502, 5.435572808224367, 51.83863037194866, 23.07806620306089, 25.897958376269997, 72.02287279261743, 19.6236304846526, 7.005586219211726, 11.002277928513996, 94.74842568371281, 7.74281225602333, 6.028533380313232, 22.112685542372233, 59.81304513313698, 20.098801560783638, 67.02394097068692, 42.698935906908744, 7.665069888285425, 32.74221539095845, 76.50258825225518, 64.3168028470581, 5.668738086946826, 20.96124466676502, 8.879310908759818, 74.17356499946771, 49.221346376723304, 5.601885656054741, 12.457948188839826, 51.63121273086832, 21.59626816404885, 29.045987854637175, 7.890816446394075, 39.88583941156479, 7.497075055720798, 8.085365535783508, 40.03898491772664, 18.604011788002854, 9.283422761519024, 12.92270949509748, 8.429712938823435, 15.05041388490439, 35.92021909367521, 6.585113692636059, 106.9576828236705, 11.38696590700949, 16.12054070604987, 42.034761131901135, 27.06318718251292, 36.624380191024045, 48.7160439801092, 51.451413288304565, 66.68661090555378, 21.94389277376897, 5.406449796565294, 107.35172135432403, 18.722833433116946, 6.796808441840091, 12.100713145284194, 13.950061446192853, 6.231845370032907, 34.074375621506455, 13.827663430751961, 37.48854967280618, 78.12907124318967, 16.34204937281738, 7.041239730499225, 5.7691920323519845, 32.78749079031385, 66.68091541050704, 7.720893733590304, 60.38101754481355, 13.571053476584124, 7.666355660969796, 26.497969802461107, 27.909465879552656, 10.807319070576611, 39.358219386930244, 70.60350674584946, 9.809221890831216, 75.86860452021843, 69.67399279712482, 6.910481710270753, 6.608223091382152, 12.91243364761048, 13.296946423964002, 72.12041954347082, 95.81114802476702, 15.897205516519364, 5.069301535907275, 68.19741733238823, 10.899620500534606, 20.29707636601043, 50.36399485573122, 36.544727755020666, 31.089975034393127, 51.88320816730719, 6.620960082306551, 16.60329826056348, 103.0130466978617, 51.677612230436395, 68.36184834658779, 108.92189378252965, 18.379577377455302, 83.95296202269797, 70.39482220895694, 8.434511155537102, 76.36968052372765, 30.721151849584732, 85.548763284069, 8.554711980866147, 6.4311020268700405, 71.17698204781308, 16.36885191314452, 12.40288550265065, 42.981381924969966, 6.9248410702123415, 8.953381285004948, 63.970669484443, 75.35321111940877, 36.965498334474475, 69.75690697715339, 113.83458708569849, 8.074519922884239, 30.835810952332757, 9.390949306939632, 11.417259048654538, 12.753091089224462, 70.01529198453608, 77.58476481748784, 14.258689394148014, 32.01306157779783, 69.86504262862046, 152.98945645962425, 5.570257562606189, 38.08495668852991, 12.057896018483525, 39.388620958892005, 29.639652371390405, 5.870843609409775, 41.33112425264589, 19.471718131828435, 113.14562860043229, 30.28685057018614, 33.759515037158664, 43.19155850871684, 126.98580088007857, 39.69923006516617, 17.08600227194007, 11.116416844554246, 8.951381780937968, 19.471593796149726, 39.79317666899776, 77.7614432325422, 42.14668212030068, 6.550829133739846, 11.858324185642934, 119.9066040758278, 19.000123957484696, 10.264522681544916, 130.40520747273678, 151.44998103190397, 165.15772611591757, 18.146878311818064, 7.487810908004603, 13.312356881497095, 5.138978616232775, 79.05825990317436, 220.1090118165187, 31.68594973928165, 13.468067934621027, 8.926914140296864, 15.7918703803895, 6.023316788380455, 10.36196221215797, 41.85958783986794, 5.478803957216898, 20.121027001195277, 63.26543130696367, 17.738124111935367, 19.706277112502534, 19.21503513222622, 38.77375264977928, 27.848397888371412, 60.06617583657638, 98.33604745253103, 16.595032312372005, 54.29437893298605, 43.98012563087998, 6.940189797563369, 25.837025479577797, 6.289861833089044, 6.419714962754574, 73.3843728778506, 5.801113005884095, 7.737575183650098, 72.98670491724789, 5.218339506760482, 20.69468482511607, 82.2307100105246, 6.597012377553708, 22.81397546037877, 10.667743261449564, 98.92222271502045, 8.217409990939213, 6.180652681595996, 68.84986478324159, 45.91825466895254, 11.120155449835211, 77.83786889974296, 57.485588909774506, 35.41370274852639, 15.760764565709696, 7.24722017950947, 15.984274366204746, 8.161749607964719, 15.618228153944655, 23.976632278741146, 11.233184554693864, 8.007397758580833, 32.10472599588564, 24.240127885489287, 82.05744596785105, 62.91635882610656, 37.320387680661106, 21.393476285731097, 22.26060080019255, 7.365793530003927, 64.55801220294916, 16.569304333473944, 14.576787632537119, 14.331588324632364, 48.70466216501629, 45.94493394182712, 5.667360793516271, 193.64722538200405, 9.435398188344804, 55.990889221244814, 9.009628455531075, 99.59417441427951, 248.89490881517486, 58.49349797509953, 60.00593587731011, 113.51263423702898, 68.13797309042805, 20.0077324586111, 106.09741697825359, 56.48700769485229, 18.51547926797228, 46.70328939388363, 117.65117821337749, 20.89872641512945, 7.493217642102088, 27.23958916291044, 25.908574700420953, 72.90226488612129, 54.5042673204575, 10.714110449652576, 60.87411767062735, 5.083866725066565, 53.229607834101415, 12.122483614857835, 12.164307978962817, 15.216499980202325, 46.172993298559625, 22.583757047089883, 62.613575454694306, 34.53071250730285, 145.6851051530556, 16.940777574703723, 7.455718247487381, 69.80790302190725, 36.94350200033803, 7.103772068591833, 38.32458585983661, 53.00645683369863, 26.842950755959173, 13.468563974407907, 11.906653908354128, 22.035188892050712, 81.35952842334021, 69.09125482061091, 6.409204107700083, 46.41560201435009, 38.86366936413045, 86.80871998096978, 231.3304920335484, 64.739500486279, 47.16959200263488, 39.92189922680159, 52.77367418676694, 83.37182640731871, 13.674604067421473, 15.824427472942867, 11.335352002008426, 44.23746601974849, 24.482587419760247, 44.863308102505336, 115.96024503114913, 5.040229801498758, 18.187678304663496, 10.2772970821902, 17.826940692978056, 52.82693640406076, 49.41140425754997, 41.97296080480667, 53.541201917598926, 7.087720714161999, 105.2664781207838, 16.116646690547043, 11.504440787395053, 79.73703577705115, 152.88396943269592, 7.909540662093617, 85.2706667388582, 24.546872932191018, 166.6844498902187, 6.788834349098631, 57.987954497553915, 81.0508993721119, 47.619455260297364, 17.634821348493546, 20.07047778707327, 14.88411817519722, 34.12683557540174, 101.93618134602555, 30.73027005453305, 6.394304293123174, 29.26588043764485, 5.470334382121479, 87.46537837562161, 36.09264715969766, 11.682306495250904, 12.685680079364975, 117.81749563342655, 44.250863102786, 19.461116570062824, 15.7381368489984, 13.700989685471963, 24.284219108385933, 44.58633695147721, 14.815937083369521, 34.43906477389972, 24.847056396466343, 26.383360617262873, 27.89048151628637, 11.809256994446155, 25.63565834795466, 19.595398949195967, 19.73640518150133, 15.193605505002047, 5.591040367122229, 86.1936680047944, 8.672903026843763, 40.02044422039164, 15.56570395327665, 18.475542416106414, 7.040118982673326, 12.876991463972603, 17.722162933238252, 69.28308449362856, 10.580484888075189, 26.722218142078496, 66.7138510916618, 46.5903978171471, 44.79068415460786, 8.183852960310698, 35.33576995824221, 86.73880655156367, 21.145190721057183, 110.44883250182596, 7.3719612906533225, 61.793933543379644, 37.92146227214558, 31.164118324241326, 11.119268333246968, 41.952819094474265, 46.9135566306072, 62.41430657897161, 68.19881229845525, 18.37940514121238, 9.445107798993, 26.896431537907926, 13.685798801704255, 36.34766209720343, 12.814461360682499, 59.742171268612466, 44.397879166051474, 28.26505195740198, 34.67525458571081, 7.643072630024772, 18.83253472872454, 16.038367433587386, 34.48297328842959, 45.51842166208217, 8.967730997798599, 57.65891998864141, 52.23216690684241, 41.53100716119738, 65.4001631657973, 7.846091723473722, 45.8406040195537, 6.017767669667718, 62.286296346126136, 25.842448020291695, 12.929823563688961, 14.48339070410606, 43.3603060183674, 20.35384217613722, 105.60504291941655, 76.90489830765591, 30.929180533842707, 157.87295333665537, 12.392376553609214, 75.805977267255, 5.860711806082114, 53.60520120849155, 5.843298738935159, 61.448032454674824, 12.368722939737864, 17.32825607331559, 62.628001570840034, 7.564014044659709, 109.77055939144033, 77.34727958593783, 87.61115620454642, 83.78687905271478, 65.01005664397863, 5.338437708586707, 248.95259818850028, 303.87017565236636, 26.748026644483023, 111.16887142124719, 84.87159686486002, 58.02941151263842, 75.64934224268958, 15.280294487647078, 8.330048786082735, 17.270155214471746, 25.664298684140483, 10.29727365162832, 7.102325968998719, 7.8664818829635905, 14.308910432711299, 83.76568580083402, 50.18201810232434, 59.15244427711543, 57.62508872574951, 63.13072273567729, 7.798843030457015, 42.238002553252336, 30.082026923529387, 21.052261276204003, 22.12646517535551, 13.765003702484783, 71.28992716923128, 7.582458370072916, 16.147090024041823, 10.574122143882033, 33.11196772016257, 22.369780309119513, 6.938501919480215, 65.90338956437972, 12.313831998677287, 48.517692697396974, 46.30659029611051, 8.253693432901054, 17.351922640560396, 10.860478297805452, 46.42816370119364, 8.620213805734736, 18.686151740465235, 6.758730354333417, 38.20753770521827, 12.088564582721174, 63.87826164965371, 9.777630956207998, 47.30639260284667, 9.47455585814036, 45.04249546431757, 38.62701258247531, 24.81178093509485, 38.972851122717074, 27.12532690816741, 52.64818331909333, 5.8891419706641415, 45.499189224675874, 5.514879119892801, 42.59221923638069, 11.075268526991017, 68.71888508661057, 34.94309316907344, 7.1332727986129925, 5.660917480754918, 97.53049157676335, 54.338705722020975, 13.79214974954378, 52.58119517160878, 52.84459384944259, 6.081807405036744, 59.74401388839681, 8.589126295423446, 5.586108584093594, 38.70532161406103, 7.377871654444003, 6.7019314631375435, 13.158141454921738, 7.219282643617524, 12.111893881070241, 59.18549410714026, 9.294823328470022, 6.452965252628435, 19.596729154146722, 5.26477138194249, 15.366741294527094, 57.21238397558369, 17.040929952848504, 12.373416945880273, 5.755764221903059, 42.28719645501298, 7.4654643087068155, 59.90143665698936, 89.83044399810169, 13.075030509901751, 8.042940643115244, 5.168067121846161, 41.21378411500444, 34.336648571056756, 7.307510794870824, 50.51950491522598, 5.848750855081385, 35.28328692281315, 54.9030741283749, 49.244725560905295, 38.220896723498385, 66.76695033044527, 12.205440893687532, 83.53019688278354, 38.797156124043106, 7.33311917442329, 36.98268087723461, 30.96196745763318, 41.02760466045344, 19.527628198511152, 6.337488447825835, 38.04608050975831])
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);
([3289616.3830499407, 3303694.5670141024, 3303706.9044706216, 3303708.62705532, 3361407.9583998397, 3372267.817054787, 3401464.450281572, 3494554.9672297654, 3544403.9626828395, 3562595.7390554203, 3588408.0159743354, 3640633.9512020834, 3653202.4190230253, 3665428.678084648, 3725898.413753792, 3745536.409961781, 3754319.835869735, 3768993.058757255, 3784231.6581447837, 3813170.791731474, 3818703.125, 3819553.685008866, 3846665.6453555063, 3868875.0, 3915306.232238509, 3924389.3778110305, 3933269.0407879357, 3948245.09856601, 4017153.2948206244, 4065865.625, 4131079.1985916193, 4134002.792875643, 4134054.4819827727, 4146324.835310962, 4146927.7973075043, 4157772.8509146525, 4157840.2121202634, 4157855.086712655, 4166128.7508831043, 4166895.778472287, 4168572.673913555, 4175763.4946510848, 4179744.7832615515, 4181995.0015037656, 4182065.6993583217, 4182268.2224579006, 4186079.4393615285, 4188457.455351227, 4188625.4633132075, 4189796.4602017277, 4198432.226202912, 4202405.9591939, 4206865.955835955, 4209252.923257171, 4210675.0, 4221150.2505495185, 4225454.900116744, 4227334.518139776, 4238702.260954452, 4239033.334343837, 4243979.72942757, 4246993.925483671, 4257229.6875, 4265934.469385223, 4273098.4375, 4281392.1875, 4291977.451712106, 4301885.9375, 4302605.221500269, 4303754.6875, 4305334.140882799, 4305680.650675773, 4306107.782168999, 4315227.670876264, 4315253.105989039, 4320772.532055834, 4321144.163587471, 4322505.8720118245, 4329816.825886039, 4332003.125, 4332479.6875, 4336782.350282249, 4339399.574673927, 4339697.638694021, 4346062.441295115, 4353772.262266037, 4367696.875, 4371000.0, 4383442.689193379, 4401287.375764157, 4401634.8334376495, 4402617.1875, 4407311.163384037, 4414490.065961275, 4417180.672459226, 4421170.3125, 4433593.75, 4433596.875, 4437880.498377094, 4442429.9356823955, 4458206.25, 4464609.514155789, 4466451.352186747, 4466484.74530662, 4477463.670471125, 4488474.663080906, 4543459.375, 4575419.491041033, 4591650.0, 4603477.212171129, 4606817.421538717, 4633991.718609705, 4687024.084459776, 4785099.650177267, 4793053.330363577, 4799740.625, 4811941.749209431, 4841086.3524588235, 4862875.89835738, 4886606.25, 4898007.515825558, 4904237.180384333, 4919047.558897727, 4941061.7465640595, 4947778.125, 4949101.992831532, 5116764.0625, 5421576.563424182, 5505224.920949464, 5547839.688672303, 5547849.769453231, 5550546.61246778, 5571444.081099645, 5573311.371803439, 5604875.6817831155, 5661937.123181856, 5673700.0, 5681827.608924588, 5706551.548781086, 5735443.876131765, 5742395.281312907, 5754239.030822766, 5754272.028870332, 5754593.632696086, 5755321.5437598955, 5796128.113457928, 5800395.3125, 5806799.823324698, 5813269.846477928, 5815070.248028967, 5815090.812571533, 5817348.587655403, 5817540.939416457, 5817585.395797357, 5819959.656112894, 5821922.787940805, 5833321.266087948, 5835439.384940442, 5838804.109323221, 5839901.042427323, 5841887.850037032, 5846614.249342886, 5849080.505958405, 5851749.98352387, 5855143.27607705, 5855253.50196926, 5856360.9375, 5856904.6875, 5857206.460665059, 5858657.8125, 5859129.6875, 5859493.6399430055, 5860915.625, 5866837.5, 5870728.023045324, 5874533.955911645, 5877444.229755689, 5877829.710389734, 5886926.587098153, 5887062.499905393, 5887134.58187558, 5889012.860171742, 5890648.334333355, 5892566.816433423, 5892569.078513558, 5893661.804914798, 5895718.75, 5900498.184819494, 5902919.586696376, 5903920.215226034, 5904478.392413667, 5905178.75964319, 5907528.125, 5909263.838586813, 5914055.060733484, 5921007.687689695, 5927692.883480396, 5928389.11303896, 5930898.439695728, 5931655.298718919, 5932479.6812824905, 5933501.483703938, 5933837.485952549, 5934342.533183626, 5934454.419233422, 5934970.43986909, 5935554.352698775, 5936195.697668868, 5937976.5625, 5937992.8899850845, 5939317.295001825, 5940470.947406035, 5941362.043840735, 5949695.3125, 5949924.898704369, 5950548.4375, 5950573.622817895, 5951059.375, 5951839.925213284, 5952471.551872035, 5954026.5625, 5958367.417918985, 5959164.9684796315, 5962646.113974172, 5962667.1875, 5965065.144417244, 5966992.709973228, 5966997.6062842505, 5967020.779085337, 5967049.810437342, 5967459.474043333, 5969865.114193529, 5973260.702807805, 5973984.375, 5974008.134318849, 5974076.774725979, 5974688.015285756, 5977938.830444098, 5980260.476401447, 5981235.0268454375, 5984757.047407948, 5985282.720587147, 5987693.937852148, 5987889.079851043, 5988233.575709202, 5993968.75, 5997159.246196752, 6009191.668260643, 6012345.053799067, 6014705.034889674, 6017794.687987047, 6017860.954466893, 6017869.868230353, 6018362.7187450025, 6025099.468489863, 6026748.25129162, 6029698.4375, 6033835.9375, 6034915.446732479, 6034948.193695159, 6035386.325138181, 6035445.770275978, 6036271.875, 6037090.625395155, 6037853.30064235, 6047157.951857143, 6049715.686131479, 6049960.023042742, 6050001.276796241, 6050111.42252504, 6050790.806351585, 6057817.1875, 6059178.125, 6059725.125933288, 6060762.5, 6062453.125, 6065739.0625, 6066374.688657359, 6066985.587188809, 6067050.219859117, 6067504.687435813, 6072065.042190622, 6076798.10204158, 6077040.625, 6081055.656353732, 6086703.125, 6087114.767780678, 6101069.2181322565, 6108098.521106405, 6112348.4375, 6114194.397461384, 6114850.974775548, 6123734.375, 6124299.688264167, 6124912.197712646, 6125829.070609204, 6128445.684862198, 6130656.201244765, 6135457.646843425, 6135724.49364819, 6135728.524719748, 6145148.108845356, 6147923.205752619, 6154363.96950454, 6158850.349198318, 6161921.565706026, 6166757.8125, 6167287.589179943, 6169733.13435848, 6170772.219290323, 6176295.706944693, 6176800.548644639, 6182279.096597498, 6185676.385094928, 6185891.586901235, 6186636.291000028, 6198443.070768497, 6198462.208253725, 6199052.238999337, 6203021.860696968, 6208991.855663654, 6212173.4375, 6233754.715131111, 6240473.4375, 6244115.659093715, 6244608.750181072, 6245601.5625, 6247979.267473207, 6256820.878299428, 6277460.9375, 6317093.744322637, 6318363.602595327, 6319417.105744561, 6323637.5, 6325873.807591194, 6327563.469539941, 6327843.309945661, 6328268.75, 6328307.16036434, 6328345.585117513, 6329215.625, 6331150.167781317, 6331560.9375, 6332701.5625, 6332804.6875, 6333047.823056042, 6334517.783022789, 6336358.787123833, 6337334.375, 6358337.087313737, 6361166.965613411, 6367449.149153101, 6380437.5, 6386279.1977802785, 6386280.308141492, 6386446.043282702, 6389031.686949426, 6390645.009116986, 6396435.536584096, 6397412.420196039, 6399119.464237197, 6402526.505205462, 6405727.470812603, 6406661.198706667, 6406668.778281935, 6408417.1875, 6413128.723050547, 6417146.64863134, 6446437.659766821, 6458698.71588326, 6459535.929230323, 6466610.302083575, 6474251.303918442, 6476085.9375, 6476734.375, 6482265.67246656, 6482708.015784061, 6483517.394575945, 6484635.9375, 6487250.731318553, 6487326.727050012, 6495400.544489233, 6524210.6505563315, 6524670.027462907, 6525123.4375, 6526838.233645814, 6528319.564374161, 6529075.0, 6529465.949582848, 6529952.810811531, 6541050.938262046, 6542268.706984871, 6543000.0, 6545800.313993159, 6554682.692061474, 6559012.264259453, 6562657.983704123, 6574273.325909638, 6581576.163124048, 6583087.284877243, 6583459.316417675, 6595304.726133283, 6596133.328257039, 6597067.4528720025, 6599028.125, 6599416.319157749, 6605305.419027014, 6605988.4041949, 6611066.676840408, 6614205.0548423305, 6626606.25, 6630750.230010365, 6639432.834772261, 6640821.096792355, 6657221.724006354, 6662053.360775384, 6667753.125, 6672033.828855704, 6675252.853577672, 6677019.486909356, 6679177.36643968, 6679564.43942662, 6684409.884115663, 6692256.25, 6699323.4375, 6702493.75, 6716439.743684723, 6718960.9375, 6722282.8125, 6749125.269646801, 6750060.253922819, 6781578.16904868, 6809678.125, 6816781.724816517, 6832516.820785397, 6833057.484971444, 6834786.267557633, 6846960.063786482, 6854871.662572202, 6858669.5077494085, 6861126.078826982, 6863315.249637634, 6863844.309015883, 6864710.586290001, 6873079.6875, 6888504.681932202, 6894310.513748245, 6896878.125, 6923192.2570928, 6926181.754821068, 6941635.214056919, 6949385.9375, 6951900.0, 6955658.747382614, 6960945.3125, 6964371.636372421, 6990299.673243299, 6993711.680432656, 6993746.842152963, 7001763.023415416, 7007029.6875, 7022373.4375, 7025400.740033284, 7025469.608496787, 7035525.83280347, 7036709.708146781, 7041413.570514713, 7046545.271019533, 7048725.0, 7049251.491930821, 7051662.188468126, 7054182.026183254, 7056485.930265008, 7072083.888557982, 7088368.195021688, 7088818.75, 7090670.128102084, 7091463.600759135, 7092691.019895673, 7095576.984524515, 7099834.480588252, 7100436.707604002, 7101899.367351357, 7107732.70770487, 7112933.772735193, 7124803.57551533, 7129161.583280774, 7130523.10667298, 7130543.294759118, 7132214.454927503, 7132735.74705855, 7132758.789284774, 7133041.890095139, 7134105.801049691, 7134111.97947878, 7134697.081540903, 7134818.10311611, 7134880.697842841, 7135036.856000439, 7135576.953527398, 7137508.725500216, 7137577.6224639015, 7142183.393114338, 7142784.116201784, 7143368.090196706, 7146214.686929061, 7174771.839883506, 7179492.065261579, 7180857.8125, 7180989.437804634, 7181790.460382599, 7185202.455272951, 7189741.449947029, 7204196.7224793425, 7215190.802299676, 7224121.431291077, 7230879.370529711, 7249851.3870703215, 7254886.294253225, 7256880.608094226, 7260231.825628557, 7260452.023194743, 7261393.035001027, 7264263.1671164045, 7264314.250142559, 7268682.49685133, 7271042.397571985, 7272323.304372533, 7278461.782847884, 7285739.0625, 7287165.625, 7287970.307663501, 7289515.625, 7289587.5, 7290125.0, 7292579.887027521, 7295330.145272437, 7296136.186589512, 7297409.13124628, 7317139.229395409, 7319046.592480765, 7321564.128949349, 7328741.619418339, 7368325.0, 7376058.434119968, 7394138.896245486, 7395987.3637844445, 7396326.5625, 7399131.963512219, 7399774.573643863, 7399799.29204587, 7399994.744605444, 7401266.887865392, 7403372.4459833875, 7442324.830215419, 7447318.6997430995, 7448412.872578683, 7448599.3626386, 7448830.125831067, 7459831.370481638, 7463375.0, 7469901.15574963, 7475951.749028853, 7476313.497937507, 7478014.0625, 7478490.625, 7481839.0625, 7490710.513538135, 7493799.336985714, 7512636.983155156, 7513703.125, 7560288.201323972, 7563457.8125, 7569170.410417212, 7572234.177227849, 7576099.685209399, 7607892.561496948, 7620949.852031229, 7683170.3125, 7684048.210021862, 7695679.405734168, 7699631.02236792, 7699697.68176256, 7702626.5625, 7722303.680115381, 7725431.825363969, 7727442.25782351, 7786434.375, 7790034.303200654, 7797412.5, 7798145.454415936, 7803120.3125, 7821280.128831979, 7841223.010341579, 7858688.966487167, 7899405.975533073, 7908079.61981421, 7923534.358954104, 7951457.524088077, 7966084.283818429, 7968991.738725209, 7969764.5095685925, 7970616.18734958, 8000632.373788593, 8000898.4375, 8011707.296970498, 8012768.75, 8043410.9375, 8053118.070887132, 8053808.697154209, 8056702.714123152, 8059010.9375, 8059679.676433637, 8078506.196328188, 8081007.1612793645, 8095090.247517933, 8101286.151340823, 8127726.691403712, 8160329.049428431, 8173951.317712264, 8196317.044492019, 8209026.172193555, 8210971.79687322, 8246708.46457428, 8262961.878054241, 8271823.657164918, 8287509.375, 8287601.5625, 8297528.125, 8341496.875, 8365681.120108673, 8365722.956547563, 8367198.715055553, 8368008.0408552, 8368169.2268717075, 8373332.5523564825, 8380617.961785056, 8382192.8950491985, 8382331.009871078, 8382349.853952357, 8383450.125555632, 8384451.122629687, 8390727.440719936, 8391429.595035817, 8391433.376010943, 8393553.470759666, 8393992.172360806, 8394164.0464347, 8403725.603455922, 8405034.193148145, 8412140.00640307, 8414426.08726184, 8419634.573110107, 8433933.333003258, 8450601.5625, 8452876.890527936, 8456736.748251423, 8470040.48705608, 8470296.875, 8472338.33143693, 8472833.006901639, 8474060.568921125, 8475544.818850242, 8476308.642575273, 8477449.966659853, 8481650.914613513, 8482790.406239, 8489335.377942413, 8494129.654625783, 8494137.8962669, 8497600.837970074, 8499541.473605614, 8504205.553811705, 8512641.274911946, 8514124.576563193, 8515378.125, 8530138.165465537, 8539426.5625, 8541384.529661704, 8552110.070445914, 8579324.681289528, 8581282.139916603, 8581331.78789907, 8582965.625, 8590296.494222675, 8590302.947508479, 8596026.5625, 8600251.5625, 8606281.25, 8632326.363060344, 8633648.4375, 8633669.845301637, 8635640.15352166, 8681954.41009322, 8697483.686836494, 8733350.0, 8757856.400109652, 8793468.813630393, 8805935.932094613, 8848099.307259345, 8848370.273104887, 8852517.1875, 8853406.25, 8870984.717912393, 8876682.943980174, 8876967.387249144, 8879285.768579409, 8879412.5, 8881633.676934514, 8882778.125, 8903796.875, 8904585.484266013, 8905578.029146688, 8909803.233810954, 8909851.089518474, 8910978.385799784, 8930865.984135864, 8943120.233602319, 8945735.47362893, 8945751.5625, 8948221.875, 8950684.375, 8950901.011943795, 8956461.562494345, 8961579.6875, 8964436.795387678, 8964688.98091937, 8965792.03532538, 8970058.72365756, 8986254.975925455, 8993529.897136461, 8994241.789232586, 9019545.249114357, 9019743.746494891, 9022372.74302696, 9026669.866807934, 9028458.494580518, 9030086.541759923, 9045727.17471792, 9048846.567883665, 9050650.388532156, 9053207.8125, 9065083.48379432, 9069584.886152979, 9083421.875, 9086801.452734517, 9089983.167441096, 9090667.83408997, 9135931.25, 9154788.540910957, 9246389.0625, 9262965.522973154, 9267862.838519286, 9291503.542262172, 9360618.653265452, 9400183.63610315, 9443471.770797372, 9451215.625, 9457009.924365016, 9460111.02363269, 9461836.018810347, 9469067.688450215, 9473146.711364785, 9502144.707199419, 9502153.051739577, 9528425.673556097, 9540988.88600209, 9658557.8125, 9734686.821164336, 9748865.625, 9766832.8125, 9782686.19552208, 9794991.135267347, 9796087.516187126, 9801835.9375, 9801857.267142529, 9805505.811231393, 9808943.176634561, 9813069.911690751, 9814023.379680492, 9814215.909012, 9815021.489402318, 9816020.3125, 9817443.461948879, 9824382.067596195, 9824570.3125, 9832667.833517255, 9834509.375, 9844351.5625, 9845117.124794083, 9848192.70278282, 9849404.509532757, 9849561.151365537, 9849703.706774313, 9851301.99159263, 9852740.601028964, 9853515.625, 9853698.976033585, 9853871.774844388, 9856639.546954235, 9857753.632811321, 9858404.852083176, 9858967.80744238, 9866631.01706779, 10200285.716260478, 10203940.876564685, 10203954.466497345, 10203971.875, 10203979.091645405, 10206241.962831002, 10206507.431677168, 10206936.18081582, 10206985.903690565, 10207082.8125, 10207179.6875, 10208039.107279541, 10208341.442292217, 10208350.506202934, 10208542.1875, 10208907.474298516, 10208964.530767387, 10209186.61808688, 10209360.87677859, 10209938.112223286, 10210013.031422365, 10210080.73375325, 10211784.794553448, 10213155.580725893, 10213496.558544423, 10213651.526987754, 10213967.547036096, 10214197.424953973, 10214393.595190335, 10214423.393919578, 10214466.9662926, 10215703.973774176, 10215760.061409608, 10215855.445624068, 10216113.747815834, 10216721.477298686, 10216828.798439866, 10217177.305935368, 10217441.950855317, 10217630.852304908, 10218371.830002315, 10218773.4375, 10220193.155714465, 10220224.508995816, 10220234.490792703, 10220881.25, 10221139.226992983, 10222200.663476087, 10224507.062049726, 10226428.584123164, 10228152.841636727, 10229960.789658021, 10231150.76025522, 10249387.257046646, 10259074.357178938, 10260149.032567868, 10263467.951217942, 10270410.9375, 10270416.373344054, 10273151.5625, 10276444.657358345, 10278823.4375, 10281074.277478717, 10284992.724302351, 71543668.75], [44.210644164147894, 27.846974509521036, 5.220422786062373, 5.3706304100837965, 30.802528362090428, 5.363967177980706, 41.18255161895853, 52.2081174609432, 51.69012048864629, 35.93237377173189, 42.86041630930448, 33.10098985656373, 38.29796372041915, 69.44594539646117, 16.149665451704976, 36.33183206792127, 13.390509199029628, 29.719384060336814, 45.40016640077884, 42.86974915662669, 44.55293979737853, 35.4894532611812, 27.720454175393403, 32.091051798830435, 5.068573459091698, 18.37151353835292, 57.47568794275436, 46.145265952271025, 15.630860145996897, 54.09951552716386, 21.60446324045682, 42.904903628683, 6.241577301293607, 11.764817155360968, 9.068883113577602, 29.883174748753426, 22.40854293693177, 11.81483609406427, 64.71026984217056, 42.0910118935903, 6.012136380726729, 30.102758320112994, 23.71449998967597, 60.97198260149402, 11.245957382949038, 65.75087197048313, 91.43340453464448, 8.106329089947515, 66.34150422829208, 22.387332594203922, 47.54997631575006, 86.16977856033377, 22.170654304356674, 25.092588603585877, 39.22409360503969, 92.14138904428381, 11.477139959364774, 20.753873450528804, 10.001175882867992, 57.11266261355941, 27.002208385685336, 67.97180569807406, 78.79409371505454, 18.881683816052277, 83.21544980167798, 44.75461636129157, 26.789189932791096, 40.74717867933141, 76.86474617647269, 61.79933660489296, 18.263199027117274, 10.304478938822841, 56.9472417332081, 36.631112412768694, 5.937751644759673, 43.44201486118911, 6.320905639273739, 22.803172474756543, 13.721735177748569, 40.508227433497254, 59.19515400731169, 36.382231658402155, 13.535321340106497, 39.94691406404384, 11.037170554341765, 16.65625306723172, 49.71196845998314, 108.88820821309281, 5.642415332488307, 10.938916393556935, 9.751651093268586, 52.53647276923455, 5.721286277003084, 5.189409482110752, 115.23603369249324, 65.78574774603757, 55.378352672984946, 34.906274976880844, 139.14706316954187, 15.980104997978856, 31.745598646882947, 13.146948483268272, 7.528901832343341, 32.226680848000726, 60.429772997049085, 23.29033501024202, 42.21305542510233, 24.136387987004895, 48.30179755727117, 54.52097846697564, 85.32389101625012, 159.90462779523733, 108.4442827088075, 68.16732582059338, 6.519772958822442, 47.150891023728036, 13.498291896837197, 11.844800315138638, 51.537446603088725, 88.29919048116699, 22.854335919009284, 7.491588904076591, 9.985389312382974, 20.004837797451174, 34.24851809768077, 27.431144083191185, 44.235827865223875, 9.332184745181168, 27.053145202409453, 18.021286163225973, 52.622352001260396, 10.740734242582041, 11.687928263704798, 7.24049669300364, 18.383577830650825, 81.15847897403148, 53.44922752254976, 33.617333176326525, 26.535482071531852, 6.507046948448335, 35.27709678446164, 39.735721401662914, 6.794619366257612, 32.24780746311441, 119.50434500685101, 72.83275584642422, 77.76701999225023, 28.654011556407667, 21.838879389782196, 34.34538451103127, 18.536446767622472, 6.793858808706219, 27.84371236596501, 8.911836949653253, 58.073703659177305, 20.120525360640436, 82.72407200754405, 16.331593418764292, 19.598327740278833, 87.28025107851832, 11.72277797237795, 8.387917777271184, 65.9514512620011, 11.102050842121374, 8.050910708446224, 48.17069176889395, 30.542284476137766, 100.79835136234841, 6.294186401995925, 36.976565960282706, 56.097859492228764, 10.959930464867906, 44.058364768863726, 72.35152166082298, 69.50166984403525, 18.602127058429893, 82.45905835565492, 106.28772679077855, 25.4062034786718, 6.088089378520371, 15.793358979485912, 16.51844564224598, 23.118472994375125, 22.883514414008562, 27.485456079277885, 29.579341634771176, 44.85806232838513, 14.058591450783448, 56.327585561219294, 10.978990819206693, 62.09955537904838, 12.28937640115151, 73.92961705683211, 52.71271071769997, 78.15277725674827, 121.28884918884194, 18.529041841460437, 22.629655183119343, 14.949885119146868, 23.13444102441966, 14.917864024865011, 18.38803349170303, 35.19150829409274, 11.71473735925323, 6.233889524815095, 42.01975445067001, 111.84436293460266, 5.635714008647698, 85.10599376329716, 28.635609867508315, 8.344256164693682, 18.076942392533507, 13.447740891730083, 67.66971193667446, 61.62834625753729, 29.6969768700478, 15.698256924915693, 31.362448813422656, 8.802363872827373, 27.68897303090911, 144.6837524292166, 22.808963573232948, 52.23874380814196, 67.53688550685831, 57.20125526090664, 51.39702285771252, 28.255098888530785, 27.786382385781593, 11.944749802571216, 25.27196033731165, 11.28612920427231, 22.45537709229452, 69.74724345510245, 42.32799280191649, 15.781817824723383, 57.03754021246817, 79.3346011648395, 23.08585676218783, 5.900231580201749, 5.651153183515341, 33.743693377089755, 10.987017291838338, 37.948924395450256, 9.346117679879173, 67.81888781804079, 74.21743506202517, 100.1782433626228, 5.439912977513531, 6.44728812011987, 11.71788629345039, 12.997497449435436, 94.52051016955855, 20.037090650223476, 33.613206608964944, 69.99106215583319, 25.321902967527684, 85.77415049461906, 39.453360941980655, 34.22149854226853, 53.211352909633746, 16.658167579684655, 5.140407748703948, 49.38174386936166, 8.286786301321563, 70.33231988932548, 6.791650750544605, 20.737834144708682, 18.16151466199128, 34.77763273274769, 5.176813548732988, 79.25918085317332, 38.03142172308502, 72.61339213368892, 24.90251215780132, 31.741901917326146, 83.34453558318724, 105.98409569219098, 23.606928389552966, 6.695528427196868, 26.71921087117419, 61.69825660042816, 63.590043111606604, 192.00083349886418, 31.96317286142182, 17.80341728818045, 31.751169230485942, 10.008609262335414, 128.4263087253355, 76.91065553398778, 32.18550864307337, 18.15774994447192, 53.17382644553258, 33.144570174884784, 37.84467778285767, 22.908055345219353, 142.2538225294796, 5.354845160751018, 17.961942437755273, 44.136180704916775, 82.88743317827695, 5.732678414261903, 72.58837936733187, 24.66356615350088, 34.530806558748836, 16.46350565448162, 13.235833375492856, 46.46984662069872, 5.10021549759303, 11.797869996364003, 6.834244725129301, 31.01565453142038, 9.842809784619199, 30.45220111249258, 30.55100661883002, 7.399555028271431, 16.997835241676178, 86.97599993353967, 61.97000533799286, 130.56083565126616, 40.820569205697204, 73.66173916848525, 138.34560812841892, 15.047882642566679, 55.36860262362, 102.68532834751502, 5.435572808224367, 51.83863037194866, 23.07806620306089, 25.897958376269997, 72.02287279261743, 19.6236304846526, 7.005586219211726, 11.002277928513996, 94.74842568371281, 7.74281225602333, 6.028533380313232, 22.112685542372233, 59.81304513313698, 20.098801560783638, 67.02394097068692, 42.698935906908744, 7.665069888285425, 32.74221539095845, 76.50258825225518, 64.3168028470581, 5.668738086946826, 20.96124466676502, 8.879310908759818, 74.17356499946771, 49.221346376723304, 5.601885656054741, 12.457948188839826, 51.63121273086832, 21.59626816404885, 29.045987854637175, 7.890816446394075, 39.88583941156479, 7.497075055720798, 8.085365535783508, 40.03898491772664, 18.604011788002854, 9.283422761519024, 12.92270949509748, 8.429712938823435, 15.05041388490439, 35.92021909367521, 6.585113692636059, 106.9576828236705, 11.38696590700949, 16.12054070604987, 42.034761131901135, 27.06318718251292, 36.624380191024045, 48.7160439801092, 51.451413288304565, 66.68661090555378, 21.94389277376897, 5.406449796565294, 107.35172135432403, 18.722833433116946, 6.796808441840091, 12.100713145284194, 13.950061446192853, 6.231845370032907, 34.074375621506455, 13.827663430751961, 37.48854967280618, 78.12907124318967, 16.34204937281738, 7.041239730499225, 5.7691920323519845, 32.78749079031385, 66.68091541050704, 7.720893733590304, 60.38101754481355, 13.571053476584124, 7.666355660969796, 26.497969802461107, 27.909465879552656, 10.807319070576611, 39.358219386930244, 70.60350674584946, 9.809221890831216, 75.86860452021843, 69.67399279712482, 6.910481710270753, 6.608223091382152, 12.91243364761048, 13.296946423964002, 72.12041954347082, 95.81114802476702, 15.897205516519364, 5.069301535907275, 68.19741733238823, 10.899620500534606, 20.29707636601043, 50.36399485573122, 36.544727755020666, 31.089975034393127, 51.88320816730719, 6.620960082306551, 16.60329826056348, 103.0130466978617, 51.677612230436395, 68.36184834658779, 108.92189378252965, 18.379577377455302, 83.95296202269797, 70.39482220895694, 8.434511155537102, 76.36968052372765, 30.721151849584732, 85.548763284069, 8.554711980866147, 6.4311020268700405, 71.17698204781308, 16.36885191314452, 12.40288550265065, 42.981381924969966, 6.9248410702123415, 8.953381285004948, 63.970669484443, 75.35321111940877, 36.965498334474475, 69.75690697715339, 113.83458708569849, 8.074519922884239, 30.835810952332757, 9.390949306939632, 11.417259048654538, 12.753091089224462, 70.01529198453608, 77.58476481748784, 14.258689394148014, 32.01306157779783, 69.86504262862046, 152.98945645962425, 5.570257562606189, 38.08495668852991, 12.057896018483525, 39.388620958892005, 29.639652371390405, 5.870843609409775, 41.33112425264589, 19.471718131828435, 113.14562860043229, 30.28685057018614, 33.759515037158664, 43.19155850871684, 126.98580088007857, 39.69923006516617, 17.08600227194007, 11.116416844554246, 8.951381780937968, 19.471593796149726, 39.79317666899776, 77.7614432325422, 42.14668212030068, 6.550829133739846, 11.858324185642934, 119.9066040758278, 19.000123957484696, 10.264522681544916, 130.40520747273678, 151.44998103190397, 165.15772611591757, 18.146878311818064, 7.487810908004603, 13.312356881497095, 5.138978616232775, 79.05825990317436, 220.1090118165187, 31.68594973928165, 13.468067934621027, 8.926914140296864, 15.7918703803895, 6.023316788380455, 10.36196221215797, 41.85958783986794, 5.478803957216898, 20.121027001195277, 63.26543130696367, 17.738124111935367, 19.706277112502534, 19.21503513222622, 38.77375264977928, 27.848397888371412, 60.06617583657638, 98.33604745253103, 16.595032312372005, 54.29437893298605, 43.98012563087998, 6.940189797563369, 25.837025479577797, 6.289861833089044, 6.419714962754574, 73.3843728778506, 5.801113005884095, 7.737575183650098, 72.98670491724789, 5.218339506760482, 20.69468482511607, 82.2307100105246, 6.597012377553708, 22.81397546037877, 10.667743261449564, 98.92222271502045, 8.217409990939213, 6.180652681595996, 68.84986478324159, 45.91825466895254, 11.120155449835211, 77.83786889974296, 57.485588909774506, 35.41370274852639, 15.760764565709696, 7.24722017950947, 15.984274366204746, 8.161749607964719, 15.618228153944655, 23.976632278741146, 11.233184554693864, 8.007397758580833, 32.10472599588564, 24.240127885489287, 82.05744596785105, 62.91635882610656, 37.320387680661106, 21.393476285731097, 22.26060080019255, 7.365793530003927, 64.55801220294916, 16.569304333473944, 14.576787632537119, 14.331588324632364, 48.70466216501629, 45.94493394182712, 5.667360793516271, 193.64722538200405, 9.435398188344804, 55.990889221244814, 9.009628455531075, 99.59417441427951, 248.89490881517486, 58.49349797509953, 60.00593587731011, 113.51263423702898, 68.13797309042805, 20.0077324586111, 106.09741697825359, 56.48700769485229, 18.51547926797228, 46.70328939388363, 117.65117821337749, 20.89872641512945, 7.493217642102088, 27.23958916291044, 25.908574700420953, 72.90226488612129, 54.5042673204575, 10.714110449652576, 60.87411767062735, 5.083866725066565, 53.229607834101415, 12.122483614857835, 12.164307978962817, 15.216499980202325, 46.172993298559625, 22.583757047089883, 62.613575454694306, 34.53071250730285, 145.6851051530556, 16.940777574703723, 7.455718247487381, 69.80790302190725, 36.94350200033803, 7.103772068591833, 38.32458585983661, 53.00645683369863, 26.842950755959173, 13.468563974407907, 11.906653908354128, 22.035188892050712, 81.35952842334021, 69.09125482061091, 6.409204107700083, 46.41560201435009, 38.86366936413045, 86.80871998096978, 231.3304920335484, 64.739500486279, 47.16959200263488, 39.92189922680159, 52.77367418676694, 83.37182640731871, 13.674604067421473, 15.824427472942867, 11.335352002008426, 44.23746601974849, 24.482587419760247, 44.863308102505336, 115.96024503114913, 5.040229801498758, 18.187678304663496, 10.2772970821902, 17.826940692978056, 52.82693640406076, 49.41140425754997, 41.97296080480667, 53.541201917598926, 7.087720714161999, 105.2664781207838, 16.116646690547043, 11.504440787395053, 79.73703577705115, 152.88396943269592, 7.909540662093617, 85.2706667388582, 24.546872932191018, 166.6844498902187, 6.788834349098631, 57.987954497553915, 81.0508993721119, 47.619455260297364, 17.634821348493546, 20.07047778707327, 14.88411817519722, 34.12683557540174, 101.93618134602555, 30.73027005453305, 6.394304293123174, 29.26588043764485, 5.470334382121479, 87.46537837562161, 36.09264715969766, 11.682306495250904, 12.685680079364975, 117.81749563342655, 44.250863102786, 19.461116570062824, 15.7381368489984, 13.700989685471963, 24.284219108385933, 44.58633695147721, 14.815937083369521, 34.43906477389972, 24.847056396466343, 26.383360617262873, 27.89048151628637, 11.809256994446155, 25.63565834795466, 19.595398949195967, 19.73640518150133, 15.193605505002047, 5.591040367122229, 86.1936680047944, 8.672903026843763, 40.02044422039164, 15.56570395327665, 18.475542416106414, 7.040118982673326, 12.876991463972603, 17.722162933238252, 69.28308449362856, 10.580484888075189, 26.722218142078496, 66.7138510916618, 46.5903978171471, 44.79068415460786, 8.183852960310698, 35.33576995824221, 86.73880655156367, 21.145190721057183, 110.44883250182596, 7.3719612906533225, 61.793933543379644, 37.92146227214558, 31.164118324241326, 11.119268333246968, 41.952819094474265, 46.9135566306072, 62.41430657897161, 68.19881229845525, 18.37940514121238, 9.445107798993, 26.896431537907926, 13.685798801704255, 36.34766209720343, 12.814461360682499, 59.742171268612466, 44.397879166051474, 28.26505195740198, 34.67525458571081, 7.643072630024772, 18.83253472872454, 16.038367433587386, 34.48297328842959, 45.51842166208217, 8.967730997798599, 57.65891998864141, 52.23216690684241, 41.53100716119738, 65.4001631657973, 7.846091723473722, 45.8406040195537, 6.017767669667718, 62.286296346126136, 25.842448020291695, 12.929823563688961, 14.48339070410606, 43.3603060183674, 20.35384217613722, 105.60504291941655, 76.90489830765591, 30.929180533842707, 157.87295333665537, 12.392376553609214, 75.805977267255, 5.860711806082114, 53.60520120849155, 5.843298738935159, 61.448032454674824, 12.368722939737864, 17.32825607331559, 62.628001570840034, 7.564014044659709, 109.77055939144033, 77.34727958593783, 87.61115620454642, 83.78687905271478, 65.01005664397863, 5.338437708586707, 248.95259818850028, 303.87017565236636, 26.748026644483023, 111.16887142124719, 84.87159686486002, 58.02941151263842, 75.64934224268958, 15.280294487647078, 8.330048786082735, 17.270155214471746, 25.664298684140483, 10.29727365162832, 7.102325968998719, 7.8664818829635905, 14.308910432711299, 83.76568580083402, 50.18201810232434, 59.15244427711543, 57.62508872574951, 63.13072273567729, 7.798843030457015, 42.238002553252336, 30.082026923529387, 21.052261276204003, 22.12646517535551, 13.765003702484783, 71.28992716923128, 7.582458370072916, 16.147090024041823, 10.574122143882033, 33.11196772016257, 22.369780309119513, 6.938501919480215, 65.90338956437972, 12.313831998677287, 48.517692697396974, 46.30659029611051, 8.253693432901054, 17.351922640560396, 10.860478297805452, 46.42816370119364, 8.620213805734736, 18.686151740465235, 6.758730354333417, 38.20753770521827, 12.088564582721174, 63.87826164965371, 9.777630956207998, 47.30639260284667, 9.47455585814036, 45.04249546431757, 38.62701258247531, 24.81178093509485, 38.972851122717074, 27.12532690816741, 52.64818331909333, 5.8891419706641415, 45.499189224675874, 5.514879119892801, 42.59221923638069, 11.075268526991017, 68.71888508661057, 34.94309316907344, 7.1332727986129925, 5.660917480754918, 97.53049157676335, 54.338705722020975, 13.79214974954378, 52.58119517160878, 52.84459384944259, 6.081807405036744, 59.74401388839681, 8.589126295423446, 5.586108584093594, 38.70532161406103, 7.377871654444003, 6.7019314631375435, 13.158141454921738, 7.219282643617524, 12.111893881070241, 59.18549410714026, 9.294823328470022, 6.452965252628435, 19.596729154146722, 5.26477138194249, 15.366741294527094, 57.21238397558369, 17.040929952848504, 12.373416945880273, 5.755764221903059, 42.28719645501298, 7.4654643087068155, 59.90143665698936, 89.83044399810169, 13.075030509901751, 8.042940643115244, 5.168067121846161, 41.21378411500444, 34.336648571056756, 7.307510794870824, 50.51950491522598, 5.848750855081385, 35.28328692281315, 54.9030741283749, 49.244725560905295, 38.220896723498385, 66.76695033044527, 12.205440893687532, 83.53019688278354, 38.797156124043106, 7.33311917442329, 36.98268087723461, 30.96196745763318, 41.02760466045344, 19.527628198511152, 6.337488447825835, 38.04608050975831])
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)