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 = 44868
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);
([5650010.149478754, 5834519.560168161, 5971456.154268411, 6083502.057448806, 6083783.929331466, 6084460.160275835, 6084930.685319969, 6085101.5625, 6135571.067988711, 6223906.842998926, 6344690.519016334, 6345945.8464447465, 6364166.678359792, 6444564.91446551, 6453035.8840658935, 6457277.852698897, 6457948.400619329, 6493423.4375, 6496719.381703425, 6528702.164729113, 6573195.3125, 6577434.375, 6628519.521133806, 6629720.3125, 6643955.915323365, 6721243.654233312, 6814796.875, 6814850.379035697, 6815310.9375, 6815325.268840411, 6817760.462303553, 6818134.064159564, 6818890.16687808, 6820594.623221774, 6986460.98362148, 6994957.8125, 6995196.769803997, 7002671.992094324, 7002809.375, 7007407.121657274, 7050857.511957153, 7063616.3680581, 7072725.149314771, 7081831.25, 7119694.159932093, 7123394.776982735, 7384442.1875, 7384475.0, 7436651.048234851, 7439820.134561317, 7517861.0886544995, 7595343.453200425, 7596008.984146238, 7596285.767276969, 7596716.815584948, 7597082.8125, 7598532.8125, 7602332.216093347, 7668487.728431678, 7672903.024584496, 7692127.105163648, 7713419.705374903, 7729113.021075299, 7739262.11025961, 7819827.290910277, 7859977.539442188, 7860738.325808369, 7863337.5, 7864018.097725859, 7905474.479832735, 7906905.637523285, 7907353.3615061045, 7921457.8125, 7944801.490511613, 7954916.387447286, 7964179.6875, 7969528.658138428, 7983981.797035061, 8122887.033496739, 8126253.312186745, 8200253.125, 8202865.250079742, 8228515.625, 8233602.87081225, 8233821.905540149, 8235215.118984976, 8235481.173556576, 8236922.449291308, 8237110.417560471, 8260482.8125, 8260904.990987003, 8260927.177773141, 8268232.3104852885, 8311040.558997082, 8311956.478988514, 8312572.319912597, 8312592.577646264, 8312655.680472765, 8315403.5521477135, 8315476.125473669, 8316313.829021708, 8326244.527067874, 8328307.8125, 8328413.549167401, 8329590.625, 8333193.660836596, 8337834.515885758, 8337942.1875, 8354527.440929546, 8365472.3354978785, 8374350.66994341, 8405806.25, 8420987.266310066, 8421052.401879825, 8422068.527561206, 8423051.549649833, 8427248.708010983, 8429481.027023643, 8432340.193654068, 8586691.456939548, 8737489.0625, 8762463.524828007, 8762514.579860989, 8767907.178109517, 8797334.375, 8823213.179879434, 8826295.040119182, 8849899.75649391, 8853120.246765753, 8871692.263507934, 8887887.508107219, 8891101.36245869, 8892801.587294618, 8892942.561309246, 8893767.864124345, 9099995.982392494, 9100006.062344579, 9143124.098537764, 9144404.269933697, 9144695.543227602, 9156366.207049087, 9159982.007957745, 9160906.965342935, 9161570.63368976, 9166040.296261396, 9166138.032359548, 9166659.375, 9166762.733817402, 9166826.062188804, 9168073.830701645, 9181113.356845982, 9191728.125, 9199630.453027934, 9204932.8125, 9206657.8125, 9208170.669990703, 9218820.3125, 9218859.06625597, 9219511.413419334, 9219611.101817088, 9219939.610085592, 9220351.66341563, 9223617.711374458, 9224421.29983405, 9226549.351477593, 9231147.49240129, 9232210.635178663, 9233712.5, 9235599.277793182, 9235734.95456352, 9252806.25, 9284781.136957949, 9288838.818928365, 9291866.651848529, 9293026.495422162, 9293049.798592588, 9294226.514245996, 9296115.877822813, 9296902.477708274, 9327290.23917403, 9331209.375, 9332060.69599856, 9332346.954932164, 9343106.404287294, 9344685.9375, 9356163.888151158, 9356571.130663658, 9358762.195290742, 9358764.645003187, 9389498.681202574, 9389681.25, 9390076.007268371, 9390237.574764578, 9391071.40672936, 9394732.12015899, 9395495.365252065, 9395539.163294032, 9395598.571587518, 9396287.871441193, 9400418.75, 9400426.5625, 9412322.061584583, 9446074.02351651, 9448266.28909525, 9449813.948676571, 9450755.092652123, 9457984.478180861, 9459356.25, 9478613.738147855, 9486556.25, 9487014.436455939, 9488039.37916052, 9488436.117195396, 9488721.166956985, 9489556.25, 9489805.208458984, 9490484.99550014, 9490625.873310875, 9490656.48221062, 9491073.21967388, 9491412.019936442, 9493596.894412303, 9493777.95081446, 9494777.422491811, 9495330.93439376, 9618505.169648716, 9639570.25366738, 9655098.288104456, 9674640.625, 9713300.869181545, 9728658.515217142, 9747364.410821691, 9750892.1875, 9756418.966568131, 9760202.441081824, 9760915.625, 9761201.5625, 9761340.654041583, 9762650.450902902, 9763651.5625, 9763932.8125, 9764363.94284901, 9764444.434632944, 9766351.560765345, 9767585.58696794, 9768889.716021651, 9784950.557352988, 9788562.223371498, 9793321.246096184, 9794737.119604545, 9796882.749256954, 9799859.220631717, 9799963.51876261, 9801450.781419005, 9801798.237598484, 9801915.625, 9805306.25, 9843769.394504279, 9858477.743810415, 9876021.875, 9940694.863522, 9942585.9375, 9942735.149072006, 9943527.040403629, 9943626.691930879, 9944053.411087243, 9944384.321161749, 9944505.38775203, 9945035.9375, 9945531.356021112, 9947254.836832307, 9949030.233178798, 9954179.124266963, 9954194.419675427, 9954239.409144122, 9954462.041654106, 9954463.339943308, 9961524.672751375, 9961790.36692223, 9961827.382252315, 9964817.1875, 9966451.638499176, 9966659.375, 9967326.5625, 9967632.103537498, 9968206.25, 9968290.205549547, 9968314.0625, 9968319.125107434, 9968752.59313484, 9969400.0, 9969421.875, 9969430.0985406, 9969941.164164066, 9970297.048098054, 9970311.516139532, 9970312.24131377, 9970342.1875, 9970446.431130908, 9970488.557121152, 9970617.083615229, 9970645.1363668, 9970745.375105875, 9970923.4375, 9971308.85948774, 9971756.278095528, 9972339.410111476, 9972407.818205886, 9972458.652366625, 9972596.705590542, 9974597.88905751, 9981389.0625, 9982980.73652565, 9985831.25, 9987032.8125, 9988915.625, 9988941.805764593, 9989247.582947334, 9989849.324026976, 9989858.384305863, 9990108.244304767, 9990123.373963622, 9990475.81743776, 9990483.316946441, 9990508.749511708, 9990760.9375, 9990763.050561476, 9990886.420268143, 9991059.817627193, 9991263.524638884, 9991459.58101752, 9991723.881299783, 9991843.994976593, 9991920.3125, 9992104.6875, 9992108.119142484, 9992264.70196745, 9992575.889160749, 9992619.353131874, 9992736.895012658, 9992928.168449724, 9992979.447692685, 9993005.026850259, 9993207.513846269, 9993232.8125, 9993307.422317876, 9993338.083432285, 9993484.5273462, 9993527.15688437, 9993656.25, 9993694.711612845, 9993721.59390002, 9993825.0, 9993858.76887652, 9993895.975861868, 9993953.003687872, 9993985.547273131, 9994023.88060655, 9994037.232658917, 9994090.501672875, 9994244.611784525, 9994353.25298394, 9994370.307084437, 9994378.125, 9994396.798699781, 9994398.4375, 9994422.140434366, 9994495.3125, 9994601.615951195, 9994644.597473709, 9994701.129912514, 9994760.369899498, 9994837.444660623, 9994894.164360018, 9994952.493571095, 9994965.090379981, 9995050.750730887, 9995254.777568452, 9995414.875678767, 9995434.73172232, 9995571.875, 9995835.1747051, 9995872.59743038, 9995920.917156462, 9996091.165238474, 9996122.620778972, 9996290.284067478, 9996725.0, 9996854.376378296, 9996921.324627195, 9997101.5625, 9997106.435018685, 9997114.042391866, 9997276.052033849, 9997318.326068342, 9998217.113290736, 9999420.164131673, 9999449.35758812, 9999955.071490226, 10000391.639353033, 10004228.125, 10010262.715491254, 10039639.0625, 10046313.25966646, 10052862.5, 10056390.364279294, 10056451.330323461, 10059515.625, 10061307.67890372, 10061959.250825306, 10063582.876356196, 10063601.409340363, 10064185.98605043, 10064382.560337849, 10065175.495672893, 10065393.42662671, 10065732.491640521, 10065739.160890747, 10065819.113149567, 10066163.652399134, 10066849.265524585, 10066896.934770873, 10067168.020069718, 10067555.357688386, 10067696.709877292, 10067821.795268534, 10068589.4791701, 10069055.9474011, 10069071.43725129, 10069245.910038985, 10069434.223772554, 10069527.45949128, 10070489.0625, 10070874.351849819, 10070885.679216096, 10071139.0625, 10071240.625, 10071601.020251315, 10072188.329063639, 10073552.373758392, 10074698.4375, 10074828.196189726, 10075592.705895163, 10075603.303326415, 10077473.397365302, 10077623.721591633, 10078073.855335947, 10083098.4375, 10083678.125, 10084057.840163851, 10084145.19574886, 10084167.44806772, 10084194.704202894, 10084343.532960378, 10085386.9143225, 10085447.625853777, 10085514.270569243, 10086857.435750296, 10086919.989592819, 10087332.068862079, 10088026.018983832, 10088787.5, 10089093.850596668, 10091479.195466435, 10091695.775936661, 10096618.75, 10097392.890450584, 10101815.625, 10103651.69833498, 10104093.38616116, 10104804.460814303, 10104851.5625, 10105229.572318457, 10105251.5625, 10105457.456109663, 10105639.768308673, 10105868.67873328, 10105902.22574194, 10106099.362332057, 10106349.744540399, 10106656.25, 10107463.14733648, 10107875.0, 10109033.422640847, 10109119.244667992, 10109580.70760156, 10109823.4375, 10109909.778691443, 10110187.318484277, 10110241.898935623, 10110497.655495794, 10110746.212709107, 10110924.230056213, 10111154.44214466, 10111203.280068876, 10111308.632147295, 10111374.280137293, 10111498.4375, 10111507.8125, 10111571.711850654, 10111657.8125, 10111862.5, 10112075.23490237, 10112075.377685847, 10112228.467436923, 10112382.29095917, 10112391.184264895, 10112573.63165682, 10112645.3750339, 10112792.1875, 10112876.794640781, 10113126.959686395, 10113251.5625, 10113283.856731955, 10113429.38398803, 10113574.471596982, 10114024.477661908, 10114046.999779532, 10114109.658655295, 10114147.233683009, 10114219.706106583, 10114345.207396135, 10114550.509708354, 10115045.825920217, 10115126.480962075, 10115234.460313393, 10115260.9375, 10115503.480005974, 10115796.56752703, 10116177.517830309, 10116262.635382336, 10116957.41898073, 10117485.9375, 10117835.761294806, 10118457.688804802, 10118720.96239377, 10120289.0625, 10120310.093685366, 10120597.814746432, 10120870.3125, 10146997.350923626, 10147626.30618556, 10149435.9375, 10149837.053968014, 10157785.9375, 10159394.603572048, 10159699.553988999, 10160205.94996986, 10160216.082186824, 10160489.0625, 10161382.263196547, 10161517.976007499, 10170836.508026848, 10183509.32610211, 10187556.969227718, 10188854.976091778, 10214745.3125, 10269737.5, 10285723.247905832, 10314731.492623065, 10322283.621778898, 10328345.69444538, 10330182.904631207, 10343612.372284446, 10344185.902296538, 10349583.463945003, 10350767.85369421, 10351103.97550238, 10358825.631777976, 10363223.210733505, 10383834.375, 10383903.316167818, 10388675.542464778, 10388681.25, 10453865.876860924, 10455998.555143412, 10457791.129829785, 10458467.72717686, 10458475.659819195, 10458498.914832234, 10458671.031893939, 10459163.503636247, 10459405.434902653, 10462026.035933854, 10462442.1875, 10475159.282629902, 10481757.644591527, 10482030.198632324, 10482558.347195486, 10484132.562542448, 10484289.61474013, 10485018.00854844, 10485802.180381123, 10501087.729559405, 10501792.177487034, 10517109.027022203, 10519802.32645364, 10524867.857694944, 10525643.458613945, 10546022.916404609, 10546355.973651404, 10546775.16587195, 10568174.299104834, 10599535.19204174, 10600677.340253225, 10602503.971500603, 10605445.9652286, 10621562.408395305, 10621798.09031073, 10622500.0, 10622726.95481798, 10630045.556617822, 10633609.13088851, 10635503.05269953, 10641754.761830518, 10641808.026231151, 10641841.30953866, 10660127.776480662, 10661268.393868212, 10665020.783303784, 10669052.833793676, 10670915.182281567, 10673023.401022632, 10673490.625, 10688121.097369825, 10688280.069494756, 10689694.853554057, 10689807.696382383, 10690450.86001738, 10692183.118455727, 10738619.225335443, 10740126.220666409, 10740348.86496705, 10741774.502098123, 10742623.19992798, 10746667.067536397, 10746681.726369005, 10747959.226036163, 10748542.1875, 10760357.216469266, 10760705.478093358, 10760706.12284493, 10762257.68513528, 10763692.741520805, 10764397.138284575, 10765473.280233582, 10766566.623278372, 10767637.663944054, 10768873.617257113, 10770480.912157293, 10772139.404482273, 10773106.12470411, 10801356.116821123, 10810680.31847346, 10818057.82757061, 10855675.119504781, 10866752.657224385, 10869820.272533592, 10869825.105518147, 10879636.888652379, 10889498.4375, 10957316.014677003, 10958726.999761188, 10959498.4375, 10960975.0, 10962308.527164439, 10963472.288977932, 10964507.65110128, 10967456.730959859, 10968744.25554954, 10986413.407669174, 11299365.381522806, 11301955.626901818, 11303869.577273088, 11304218.696072562, 11304447.071286488, 11304654.544997586, 11304854.6875, 11305487.5, 11305925.320724953, 11306151.320066718, 11306480.323739847, 11306571.941971555, 11307924.157067388, 11308049.508752769, 11308244.767764298, 11308260.146917308, 11308965.311162245, 11309408.932880666, 11310209.375, 11310726.030202148, 11311353.88224646, 11311949.663972778, 11312306.25, 11312491.230177203, 11312514.051928347, 11314136.948383793, 11314631.638811775, 11315104.6875, 11315527.635495024, 11315791.607703395, 11316262.113818437, 11316498.4375, 11316556.137962764, 11316639.621954117, 11316877.010217378, 11317138.209989619, 11317562.423130447, 11317742.1875, 11317797.313866088, 11317854.338882864, 11318868.296620138, 11318883.235805903, 11318892.09450813, 11318952.809462044, 11319114.0625, 11319315.625, 11319867.1875, 11320009.189201688, 11320108.498884661, 11320634.788229933, 11320771.607512483, 11321099.215910107, 11321482.471323596, 11321572.811270539, 11321632.8125, 11321682.029790265, 11321968.75, 11322571.46730476, 11322980.679594887, 11322985.9375, 11323009.10437657, 11323111.781291237, 11323589.283322692, 11324180.513131185, 11324210.202272784, 11324705.549954517, 11324733.337303597, 11325375.896355465, 11325583.484296152, 11326021.875, 11326170.981356645, 11326555.425578406, 11326578.24763215, 11326934.375, 11327150.979013652, 11328200.76864632, 11328212.993038382, 11328608.806487942, 11329303.363750942, 11330243.75, 11330381.535327647, 11330713.059969002, 11330759.143100874, 11331496.74701107, 11331866.700925088, 11332270.70851687, 11332309.885536298, 11332541.716856275, 11332583.95454559, 11332624.03006431, 11332627.27789879, 11332996.64112055, 11333486.340197692, 11333699.299344273, 11333730.746825922, 11333815.968681837, 11333909.955000104, 11334237.5, 11334249.703637168, 11334283.274657663, 11334420.057393488, 11334540.625, 11334964.009698339, 11335645.053888284, 11335749.431585792, 11335822.9727387, 11336065.774919787, 11336319.708718732, 11336481.697861573, 11336979.554605737, 11337072.860621884, 11337146.60124551, 11337338.00072358, 11337425.034825966, 11337498.4375, 11337644.776351845, 11337901.116662104, 11337903.306146994, 11338095.440329913, 11338645.32412605, 11338662.156389475, 11339012.5, 11339067.608608138, 11339285.724739525, 11339360.9375, 11339452.626922332, 11339490.924699571, 11339562.5, 11339581.697393106, 11339769.685084667, 11339853.125, 11339943.120961653, 11339944.52839517, 11340029.77130868, 11340152.895697717, 11340161.480879674, 11340293.480582705, 11340665.688704155, 11340739.0625, 11341152.17201815, 11341167.842062626, 11341202.56663367, 11341265.993473988, 11341504.745027874, 11341505.067141082, 11341590.17264421, 11341862.809718678, 11341984.06117406, 11342032.021605086, 11342241.64630901, 11342264.777063593, 11342274.826835481, 11342311.992790975, 11342432.8125, 11342614.974061603, 11343007.492282942, 11343102.234121814, 11343145.59069278, 11343154.341160363, 11343314.204841234, 11343371.875, 11343418.584700938, 11343662.967322119, 11343852.793219458, 11343990.625, 11344025.8375701, 11344203.125, 11344258.427627277, 11344298.4375, 11344323.4375, 11344761.321885636, 11344792.1875, 11344839.856737968, 11344876.5625, 11344994.58869172, 11345162.41240205, 11345370.3125, 11345594.51278295, 11345703.193269689, 11345719.304799426, 11345840.077570412, 11345943.101627449, 11345959.375, 11346131.25, 11346466.715251409, 11346707.8125, 11346767.070324406, 11347060.789394204, 11347365.982702896, 11347460.53183657, 11347510.261861237, 11347673.706304615, 11347725.0, 11348195.669537522, 11348232.378096377, 11348482.564315438, 11348830.730562022, 11348889.577241935, 11348903.125, 11348992.123475896, 11349071.739978774, 11349182.8125, 11349206.683818387, 11349341.734150752, 11349427.965271123, 11349446.729496295, 11350298.4375, 11350414.510555848, 11350623.77389559, 11350801.923099872, 11351658.773174861, 11351670.998425651, 11351731.25, 11351977.889988754, 11352292.35176889, 11352641.288735265, 11353291.287283603, 11353715.202113796, 11353857.825136425, 11354026.5625, 11354624.845397815, 11355235.9375, 11356275.029038105, 11356576.363725629, 11356862.593597626, 11357164.0625, 11358264.0625, 11358537.40964914, 11359179.6875, 11360395.126303079, 11360848.762718486, 11362070.62716068, 11363596.591489922, 11364048.342244009, 11364312.403429532, 11366046.994555661, 11366349.435707312, 11366373.4375, 11370281.832267223, 11373352.085883295], [95.58548775922932, 6.01320814014058, 26.765112195213412, 9.788181407113093, 5.678969627754286, 21.11154063307883, 90.73018489566181, 49.897952180760534, 25.50329536951411, 22.741834911686922, 8.231315455462765, 7.269258430875213, 81.9088647427757, 17.3187151941055, 11.045200535918978, 30.66751517463912, 27.815392897821855, 42.52925686639747, 12.283079923645813, 17.830656342819303, 37.84698385358875, 123.6001816899272, 41.51264996786801, 28.56746001451353, 46.86364705168485, 26.764513833847776, 78.2581655073064, 5.174284090745851, 29.773427653344168, 16.075801224118013, 10.458170872290466, 43.14641012019057, 39.749096594798566, 14.040568931586025, 73.98951768059624, 32.27034338585554, 14.55515864855923, 46.749646624558295, 53.42289563900929, 5.664157719912628, 10.68833136577545, 149.96907136114402, 46.98130250958803, 38.804211890181676, 35.0173583891157, 13.22981389641947, 51.80859772125987, 51.82465103430784, 31.781066675546914, 25.446056359838387, 17.814735506626903, 38.1774822203444, 56.86897106336503, 58.2602287757091, 7.474042082945834, 53.479486689804745, 55.73951817543015, 38.18324240259955, 24.967019355341698, 33.566006764085145, 19.199306741890528, 12.343809839921356, 9.814625594487591, 20.35658623960603, 33.32650955788597, 61.69900234816275, 19.417606529172566, 69.6825829271247, 22.2065373772513, 65.32627937681207, 36.061117113495655, 9.539722649028153, 69.05124291439263, 7.577406587036659, 10.071634792171823, 44.46341184716933, 8.590178370972566, 13.372913300220752, 23.28309007305858, 100.11373812404942, 51.01383026491613, 77.61887752503485, 46.212337902206954, 30.867307781519802, 11.175780931783569, 19.57170228450253, 78.89369102134931, 9.869290046391244, 21.296472417326097, 48.80659837119738, 6.306318633333746, 49.90161562691786, 24.922796391755746, 27.283648041756706, 11.505543740904177, 6.325095405214135, 29.019430609902617, 68.15491816240488, 65.94902961050941, 64.89453855761684, 53.61127262161861, 19.051046917230867, 29.481889498911784, 46.42234492087289, 28.815407357444617, 70.88415823893075, 110.47118024275929, 65.74147429598145, 42.6364094170515, 30.212018785359984, 12.36577024392671, 104.32107943374065, 18.30233207765047, 63.91274109107497, 24.794755863970785, 111.12885106443973, 11.578605699492982, 14.154683222497548, 64.35268703272439, 6.17910839145654, 63.01349029833551, 43.825318769784396, 5.078426671529605, 44.24566149825455, 56.03139671960262, 138.50269417923613, 96.23813617210975, 198.07239738402617, 14.696743887930015, 15.23935242686958, 14.972500645475694, 61.19026072131363, 21.637365401738528, 160.06320391432413, 41.045711239217376, 33.114224119959104, 9.17761519570235, 5.566046576201569, 13.557628486144786, 8.412478744502902, 39.35664263790514, 42.061365482654054, 32.87476692605246, 16.325497913605545, 13.820910730421968, 12.076010612986154, 35.10785848654111, 6.112756425421403, 13.325193823119935, 264.74807263510854, 12.79601653400451, 40.14764222771586, 83.9460908223016, 116.20172193394596, 74.43921570907858, 52.616671677346375, 224.84678198052134, 28.194637745983336, 19.276055945969517, 7.663769565687397, 12.107044275119144, 51.61185008389951, 26.629857361099774, 13.124792858251917, 29.515648633614475, 78.75987726144922, 123.33062458263726, 36.029687422746925, 75.55449663972391, 5.362181333375631, 31.088195404092065, 126.92449351060091, 35.861032510124005, 176.60015005288042, 27.264481560403773, 37.84110985049236, 27.6048200080266, 16.001465160736533, 90.18361792058003, 71.8734032667486, 58.26359680918382, 10.770101308431876, 35.026196825688075, 15.564104172285408, 39.85817288464503, 6.729191290195267, 66.36906532250043, 7.4966261094217685, 12.179161693437239, 20.942137764356083, 37.339336876278814, 90.24750834357965, 6.276690770610132, 75.39521299448896, 5.76446945789797, 6.774084005075873, 9.41810963680264, 16.436306260368188, 21.787419941102513, 66.54629681361676, 59.21481421731108, 106.62541280526696, 5.85865134769586, 16.80204722086571, 7.5625820223898925, 25.38565067096753, 33.23370654126729, 92.50312835932648, 8.122192309905532, 31.40507364905291, 8.494753137229836, 11.661576689890573, 9.488112687438981, 6.181703289711271, 29.173348164512305, 30.17042610418095, 38.425251737339636, 33.4520839141225, 15.907180055332487, 114.69670582323175, 51.21357437909144, 61.6031484299749, 96.74937378926643, 48.86236825031578, 5.236397153462196, 29.576542829449565, 11.036976191641012, 63.41734779131698, 29.35995614199963, 26.564578341477137, 18.524977876812027, 22.545814380852043, 31.706457136129927, 7.1795687932909615, 43.25854617750927, 75.03600460804275, 109.93453165680683, 73.91349444542976, 6.856993140199748, 76.24621826559901, 45.795743089345635, 27.560391720092053, 5.2212223837923375, 47.116742998281744, 72.18998378183838, 100.63466415232244, 7.276098116403006, 10.721200019199284, 101.42385544333827, 9.02921615668905, 26.897493166125052, 30.38868642867213, 13.207834777908795, 6.930226797753219, 64.68278419449726, 57.702700457687484, 46.894195711599025, 18.138632257571757, 33.151425054729835, 67.16275861407048, 10.491652017175737, 38.15903646208226, 48.08883722746971, 8.556574467817198, 68.35135883896513, 8.443889917671592, 23.233219970620603, 5.568059358134129, 35.912014140724324, 7.149603088663292, 17.475357913009596, 14.779538202096411, 20.8846859106485, 7.864843757600646, 6.1249045561263475, 21.9189355475237, 19.527351285713518, 9.069220041293361, 24.301429480282415, 12.706998779915972, 60.547611765134675, 5.618262478328589, 107.61871791756307, 72.23315411574018, 19.11288459556006, 102.54776417174521, 21.736156647862202, 50.67225888077736, 33.85651913419626, 21.181010519665563, 59.275636181897724, 37.086590796509356, 5.9760190035442555, 11.266195025571879, 8.380285959610037, 12.171116055817329, 10.746084677311847, 105.59855072567083, 30.210810115698635, 60.7040415939613, 32.589259070148984, 14.31013979260201, 99.01210311918832, 49.66003759874583, 122.93518456879336, 23.498936969959317, 12.053544208142327, 16.652544789868394, 10.159255604005269, 23.23785843650893, 27.22507472905824, 60.96400267947714, 39.38426921864413, 54.20304615966998, 50.763985140826506, 39.02303532532699, 16.28661535070851, 12.920814411473119, 12.868620574685957, 6.342250044523681, 125.49674680821963, 9.278251253417247, 17.613568130270536, 8.572656141268673, 21.140787086893837, 47.11218711614873, 12.113138810269797, 83.97328998623063, 139.01300278627, 56.04949557296476, 5.719996171185611, 31.238889458630997, 21.12007042858298, 50.72593412951427, 43.18730863509642, 9.604882570581308, 30.716882985835497, 10.216453890052513, 56.4812540649851, 20.506048510788954, 15.109746615864445, 35.9843912378238, 135.9979704949151, 55.89875107009698, 78.75195589152028, 36.308875491691225, 31.3656281842295, 66.19876849655732, 9.890214795466772, 41.18078433587508, 79.62144003337747, 6.708655394277812, 53.89189228736277, 6.396096098922511, 15.060218313243986, 41.02967145674825, 22.662803281844567, 27.86958923202867, 12.889665676830223, 14.491395193118471, 37.99854475407462, 5.950629174574408, 14.92090793129184, 52.076502216920126, 27.101347081706365, 70.63730989691489, 16.059914905282625, 67.17754230553906, 36.068851712370616, 8.792086688270468, 10.512769497115887, 5.8914540618017694, 11.042295654218115, 8.518566440764177, 5.320740686105919, 19.05787470946104, 145.24454428424082, 89.56547786858984, 19.320882084038857, 28.521089308763077, 51.94591047240775, 19.219051881378086, 16.186385872354094, 98.44645748190469, 17.45006624830529, 14.32915017019144, 10.633580295009876, 137.56833595382767, 13.874923312536906, 8.015503428408152, 31.01012853103994, 15.69694805849591, 22.797877385056463, 33.93925493859352, 26.301537205196865, 27.138130102837255, 35.077320820453394, 7.973761450050758, 7.748286814051039, 7.416325404286163, 61.25760425657081, 9.520511545841945, 40.85927720050594, 41.81559362738416, 82.85331128932124, 24.265791017400865, 5.125570959490278, 32.69628643760449, 45.060795463443206, 26.29047169950258, 5.464514009237976, 33.452444669066345, 12.344969207041158, 78.17332370301033, 69.49499402274404, 5.448056088717443, 7.491169318257705, 5.043205903871608, 38.91842410249734, 7.461353237116094, 19.434708462344425, 35.01865395012505, 32.357277564802914, 28.176869925346246, 64.28941478223892, 33.2438187724189, 7.761228369551284, 23.7176009968896, 13.499833596128612, 67.30577703002992, 51.863512898213585, 7.379052570959108, 42.78061405216009, 33.02101214031829, 35.842196481456746, 32.385305422371175, 57.078431816481405, 161.7281037254038, 12.729167612891814, 12.686877954553186, 60.65902570679487, 105.92133582940134, 91.66793008080981, 25.0713097564618, 176.58207556185698, 114.4495993609547, 43.91672995556732, 42.648661816976094, 73.25333899556912, 15.04800994535002, 6.252559538475853, 6.816666256715689, 10.2936674881045, 72.83785717590803, 20.678977775034696, 5.656979058769342, 5.20875549974933, 59.56930418211217, 5.036144013547315, 89.4410596320624, 8.910550662556759, 79.61249534958054, 158.953522487227, 38.89553123657973, 11.891976055232082, 32.51902827313742, 42.50173455133175, 51.768105209041416, 8.359739541032159, 18.792758365908327, 10.7909520127425, 56.516792451076626, 26.438357377844845, 28.797765205316107, 21.638547727792364, 35.770771310280026, 52.69510322082316, 18.260673353459126, 32.004226143666656, 9.12435144047788, 45.110148939382505, 6.593591838807247, 78.58023606595353, 25.536544332342814, 78.92193065172485, 10.373993851514802, 61.302791779656886, 5.303722029610844, 51.56131299212668, 46.790762187440635, 88.99318852706148, 72.30880894227006, 118.47055344671813, 31.26929724355679, 15.56744423898236, 82.38581898791122, 7.939362542543894, 58.940910914592614, 53.38361153894469, 119.63238262541284, 79.01991319109841, 54.6157482179556, 15.919406153040544, 16.60717126278108, 16.88564322394499, 47.02352888244826, 31.234239467077636, 12.011899666194765, 25.326128321192545, 78.42571378639701, 15.906863184942546, 5.413878316242108, 45.97958177036469, 29.543411563263188, 98.42420682483484, 35.379980968371264, 31.17923256835921, 25.216557353781717, 9.58334736594936, 8.48582207539188, 87.96417575924038, 5.5280011701501826, 58.27533174749955, 36.43183876523967, 27.009726848809734, 15.27697001935811, 71.15503015452435, 16.492520622570293, 107.3734503402859, 20.4867343855342, 68.157248743764, 182.2920442956025, 66.0601889106289, 25.47226777996774, 58.505783107017486, 81.73120136805935, 29.366158085604415, 31.619177054837813, 33.46372787313438, 28.95193994373235, 9.769084962223468, 21.368649730129807, 56.44857157268002, 92.78314822686065, 72.82989824442224, 24.489155650354768, 22.049123494390514, 66.89450693555474, 16.169955120572553, 54.41324920566595, 30.768457406630272, 10.100833009975094, 86.9054797155474, 55.36070078676495, 7.324283903256295, 45.31854305179026, 55.578670411829435, 72.54372313386452, 25.27149230156136, 131.6503218582221, 12.682202717497331, 11.782251008041401, 11.274833560107759, 21.322399422146734, 116.99424615061356, 8.663985970419603, 9.96630826055913, 64.53353624486206, 124.7523792763555, 10.790837936212586, 92.500615809082, 6.802281452408287, 19.284540057105303, 61.70470238007596, 26.39444736863409, 17.236326423447043, 44.79871466273643, 7.270930360500436, 68.90404064376018, 11.965891461129416, 7.324334194813069, 67.49881262634817, 5.14680346494961, 163.47564258288008, 76.13589337236677, 14.600750460898963, 9.193351011699438, 6.870028675659468, 12.080755602584505, 24.33069006786257, 17.56761288808604, 20.398174901294222, 12.253608036877603, 8.411821081548393, 28.328305173835357, 20.99847034933253, 29.406571864624915, 12.366319859910263, 81.65464948808815, 21.285212332799194, 17.63073038266638, 52.64771918123475, 19.198062660930095, 69.37996136091381, 25.438912771670708, 12.770999764020226, 9.951880963839702, 26.83568947122827, 23.16611558438287, 72.345140270634, 11.801791005427297, 40.159147240199, 59.11758679972437, 86.12883343215259, 71.35747968356017, 6.521125654130583, 5.221219882149466, 31.047729688916185, 13.738128978017885, 66.77000888528026, 5.9144404740612355, 27.16643538991025, 26.672457857994292, 28.708948422867486, 98.46497337421395, 16.62919843951085, 43.91720938313228, 10.950931840026188, 80.42849592859929, 25.21774000490931, 23.123175771779668, 51.11435363076505, 5.480417979560089, 17.525905081018315, 14.107230599485103, 8.224141269738592, 17.132286034146492, 9.208543931185579, 47.70344446213627, 59.03297498091047, 27.424505263211962, 32.958495060804424, 6.467309336384388, 7.2751225349049, 11.565783437945672, 21.67948954489315, 61.59516238869725, 8.375219133435062, 15.67398105776758, 68.76331881299625, 31.980735984186868, 14.709822313106596, 18.435105410889683, 109.59137669367658, 14.999264351983173, 56.28980978265738, 27.453195065361776, 17.539507408194122, 7.149468195061383, 46.46736579755781, 40.184085918192665, 11.796806184412539, 25.42049532496941, 102.45735908167347, 46.40520703331625, 5.868821945680249, 17.24827685447288, 52.683223064364164, 17.157811201253367, 20.823404588067127, 81.67412686537327, 24.50062980827529, 20.34715011691696, 8.778536578778262, 33.03300697173112, 68.34388407751757, 58.252814893778876, 112.12276566229554, 52.883954058297086, 67.92978988816171, 60.956247748653674, 16.27514095951107, 7.1391684319463655, 36.919492080456116, 21.580860398078933, 22.205908193750876, 24.931550920773425, 45.31653199472129, 8.9835968035357, 32.00745729517433, 21.176850823873586, 6.6127307054415425, 23.301199739960364, 138.29785540594224, 47.55823656643202, 85.09085088756143, 10.362627903095994, 5.138899690693089, 94.37581346522359, 36.18550111165656, 22.16499015027219, 129.27161890574607, 70.12184719594458, 14.531005585558676, 5.178598961409752, 115.45123448990353, 6.604462138580024, 81.59225584292057, 31.16867336193774, 5.214748963641918, 6.69584799579283, 20.526880215232783, 21.937115231024713, 7.097219609729197, 34.71690450382481, 31.529646715218522, 79.12929204666194, 29.50968445039194, 10.849201693992526, 64.29661454383167, 6.848907499410452, 14.011390419950876, 47.99501695574824, 137.23777012297907, 20.319993441084975, 59.65511042711247, 10.091033069158863, 35.61584726199297, 66.78061299889751, 20.81839160101983, 105.43718216639098, 14.001452663280446, 27.15013416921883, 18.710976137781948, 12.720981814676621, 26.176425749219916, 12.809640990324379, 6.873006365426899, 12.892937941431997, 7.865738779360716, 53.56929158231249, 18.297664094188, 14.926604661189653, 18.242438316215804, 62.07744899181844, 20.698834591709506, 12.587378855915636, 11.997001316953039, 72.20114083284584, 87.6907879172194, 49.535303951469395, 8.443455646315252, 30.77925384189016, 57.85321765818091, 7.557593056045657, 113.86429564489136, 22.984427442071848, 30.312589646906297, 69.91433367623092, 64.07142047762008, 17.521213902119353, 8.704961254820354, 31.747084376410644, 12.000699987176372, 19.70058260056357, 117.96735908866785, 5.254867144015603, 6.1505649303763485, 46.67235919656056, 32.4258540925215, 69.63705165692747, 6.061966156429703, 66.24250539716236, 27.521675195333707, 8.154395348976237, 13.142501228150062, 10.47145290066898, 15.545476288124364, 20.516742290278145, 30.00316521699822, 89.1655863490199, 20.75901348127942, 19.156241762018105, 73.2217002040056, 5.7089780152054725, 59.62651006765546, 31.943543617359417, 5.3268730825534405, 51.2153143539848, 11.33356983528869, 9.341253040666418, 7.483514153169173, 46.604079096058314, 9.714804682221267, 62.30186206231812, 42.87629386283293, 21.51229797945035, 47.78283384141179, 44.2744473674061, 25.46388188398381, 6.64439504679244, 71.5468006258162, 7.390538192955694, 19.06079285878362, 5.344801724438608, 24.697862697005466, 22.91386395661366, 43.111995734406435, 8.320554913189627, 33.061064462201195, 24.926030854373025, 18.20387650508765, 26.524144354786714, 16.560343921427915, 8.317235131842342, 9.674457726642478, 39.840771764066275, 87.8738376103924, 91.72415911110124, 8.74161871871524, 8.912199984277576, 15.286799238946447, 25.518402246944248, 41.69141167957219, 132.55597243374925, 12.930834864720376, 22.81406691005065, 52.63481115341691, 6.820346238981218, 32.06798901368958, 9.454732399542298, 37.73975498520255, 46.69377487428331, 77.96678160580464, 23.35310383909254, 67.34168947595097, 47.76360623474593, 53.90657301214392, 5.994826541823762, 55.161183763386354, 76.27979351579964, 6.549083870893113, 34.36354576044011, 7.9181949147842365, 64.08952834056191, 49.07352311084684, 10.993491515691165, 33.08664536182323, 10.072008994255157, 15.208266772590372, 83.08232848376191, 21.842763589696645, 60.444153522845774, 34.822462065763595, 37.657679340168386, 21.774475513034307, 48.34254379271037, 39.91819688998364, 25.90641389445405, 17.36252283630768, 5.0615930481436395, 27.173108426884962, 52.435442466146256, 33.88215387300342, 81.27900880248808, 5.155366586274538, 24.43597287863193, 5.378980868360515, 5.8530744537476895, 62.92176614430332, 129.0556306474168, 26.114720303349543, 43.85820173891202, 6.323609522200777, 80.45940114719376, 60.930013467182775, 23.7790416929426, 43.69631823556113, 11.908571830062808, 45.11868435813659, 14.857690501311453, 65.1255332165766, 7.314242139581658, 145.77144785234924, 24.617964669026865, 5.044323075070157, 48.69530913986638, 52.42387066076883, 13.542323634861507, 9.342110889768392, 31.305135271550906, 16.744092944203725, 48.256385416231005, 14.970547589171723, 30.625229406975112, 8.33670136818492, 50.69497731318576, 67.97346984870039, 39.25612033985715, 59.07382116471256, 25.322073977004504, 15.280623116392896, 6.121444034932263, 13.960828120380441, 50.61813338766851, 64.51275905420002, 5.620692663652788, 6.65584700129912, 38.02226174977707, 12.178669851864264, 58.865179229923214])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([5650010.149478754, 5834519.560168161, 5971456.154268411, 6083502.057448806, 6083783.929331466, 6084460.160275835, 6084930.685319969, 6085101.5625, 6135571.067988711, 6223906.842998926, 6344690.519016334, 6345945.8464447465, 6364166.678359792, 6444564.91446551, 6453035.8840658935, 6457277.852698897, 6457948.400619329, 6493423.4375, 6496719.381703425, 6528702.164729113, 6573195.3125, 6577434.375, 6628519.521133806, 6629720.3125, 6643955.915323365, 6721243.654233312, 6814796.875, 6814850.379035697, 6815310.9375, 6815325.268840411, 6817760.462303553, 6818134.064159564, 6818890.16687808, 6820594.623221774, 6986460.98362148, 6994957.8125, 6995196.769803997, 7002671.992094324, 7002809.375, 7007407.121657274, 7050857.511957153, 7063616.3680581, 7072725.149314771, 7081831.25, 7119694.159932093, 7123394.776982735, 7384442.1875, 7384475.0, 7436651.048234851, 7439820.134561317, 7517861.0886544995, 7595343.453200425, 7596008.984146238, 7596285.767276969, 7596716.815584948, 7597082.8125, 7598532.8125, 7602332.216093347, 7668487.728431678, 7672903.024584496, 7692127.105163648, 7713419.705374903, 7729113.021075299, 7739262.11025961, 7819827.290910277, 7859977.539442188, 7860738.325808369, 7863337.5, 7864018.097725859, 7905474.479832735, 7906905.637523285, 7907353.3615061045, 7921457.8125, 7944801.490511613, 7954916.387447286, 7964179.6875, 7969528.658138428, 7983981.797035061, 8122887.033496739, 8126253.312186745, 8200253.125, 8202865.250079742, 8228515.625, 8233602.87081225, 8233821.905540149, 8235215.118984976, 8235481.173556576, 8236922.449291308, 8237110.417560471, 8260482.8125, 8260904.990987003, 8260927.177773141, 8268232.3104852885, 8311040.558997082, 8311956.478988514, 8312572.319912597, 8312592.577646264, 8312655.680472765, 8315403.5521477135, 8315476.125473669, 8316313.829021708, 8326244.527067874, 8328307.8125, 8328413.549167401, 8329590.625, 8333193.660836596, 8337834.515885758, 8337942.1875, 8354527.440929546, 8365472.3354978785, 8374350.66994341, 8405806.25, 8420987.266310066, 8421052.401879825, 8422068.527561206, 8423051.549649833, 8427248.708010983, 8429481.027023643, 8432340.193654068, 8586691.456939548, 8737489.0625, 8762463.524828007, 8762514.579860989, 8767907.178109517, 8797334.375, 8823213.179879434, 8826295.040119182, 8849899.75649391, 8853120.246765753, 8871692.263507934, 8887887.508107219, 8891101.36245869, 8892801.587294618, 8892942.561309246, 8893767.864124345, 9099995.982392494, 9100006.062344579, 9143124.098537764, 9144404.269933697, 9144695.543227602, 9156366.207049087, 9159982.007957745, 9160906.965342935, 9161570.63368976, 9166040.296261396, 9166138.032359548, 9166659.375, 9166762.733817402, 9166826.062188804, 9168073.830701645, 9181113.356845982, 9191728.125, 9199630.453027934, 9204932.8125, 9206657.8125, 9208170.669990703, 9218820.3125, 9218859.06625597, 9219511.413419334, 9219611.101817088, 9219939.610085592, 9220351.66341563, 9223617.711374458, 9224421.29983405, 9226549.351477593, 9231147.49240129, 9232210.635178663, 9233712.5, 9235599.277793182, 9235734.95456352, 9252806.25, 9284781.136957949, 9288838.818928365, 9291866.651848529, 9293026.495422162, 9293049.798592588, 9294226.514245996, 9296115.877822813, 9296902.477708274, 9327290.23917403, 9331209.375, 9332060.69599856, 9332346.954932164, 9343106.404287294, 9344685.9375, 9356163.888151158, 9356571.130663658, 9358762.195290742, 9358764.645003187, 9389498.681202574, 9389681.25, 9390076.007268371, 9390237.574764578, 9391071.40672936, 9394732.12015899, 9395495.365252065, 9395539.163294032, 9395598.571587518, 9396287.871441193, 9400418.75, 9400426.5625, 9412322.061584583, 9446074.02351651, 9448266.28909525, 9449813.948676571, 9450755.092652123, 9457984.478180861, 9459356.25, 9478613.738147855, 9486556.25, 9487014.436455939, 9488039.37916052, 9488436.117195396, 9488721.166956985, 9489556.25, 9489805.208458984, 9490484.99550014, 9490625.873310875, 9490656.48221062, 9491073.21967388, 9491412.019936442, 9493596.894412303, 9493777.95081446, 9494777.422491811, 9495330.93439376, 9618505.169648716, 9639570.25366738, 9655098.288104456, 9674640.625, 9713300.869181545, 9728658.515217142, 9747364.410821691, 9750892.1875, 9756418.966568131, 9760202.441081824, 9760915.625, 9761201.5625, 9761340.654041583, 9762650.450902902, 9763651.5625, 9763932.8125, 9764363.94284901, 9764444.434632944, 9766351.560765345, 9767585.58696794, 9768889.716021651, 9784950.557352988, 9788562.223371498, 9793321.246096184, 9794737.119604545, 9796882.749256954, 9799859.220631717, 9799963.51876261, 9801450.781419005, 9801798.237598484, 9801915.625, 9805306.25, 9843769.394504279, 9858477.743810415, 9876021.875, 9940694.863522, 9942585.9375, 9942735.149072006, 9943527.040403629, 9943626.691930879, 9944053.411087243, 9944384.321161749, 9944505.38775203, 9945035.9375, 9945531.356021112, 9947254.836832307, 9949030.233178798, 9954179.124266963, 9954194.419675427, 9954239.409144122, 9954462.041654106, 9954463.339943308, 9961524.672751375, 9961790.36692223, 9961827.382252315, 9964817.1875, 9966451.638499176, 9966659.375, 9967326.5625, 9967632.103537498, 9968206.25, 9968290.205549547, 9968314.0625, 9968319.125107434, 9968752.59313484, 9969400.0, 9969421.875, 9969430.0985406, 9969941.164164066, 9970297.048098054, 9970311.516139532, 9970312.24131377, 9970342.1875, 9970446.431130908, 9970488.557121152, 9970617.083615229, 9970645.1363668, 9970745.375105875, 9970923.4375, 9971308.85948774, 9971756.278095528, 9972339.410111476, 9972407.818205886, 9972458.652366625, 9972596.705590542, 9974597.88905751, 9981389.0625, 9982980.73652565, 9985831.25, 9987032.8125, 9988915.625, 9988941.805764593, 9989247.582947334, 9989849.324026976, 9989858.384305863, 9990108.244304767, 9990123.373963622, 9990475.81743776, 9990483.316946441, 9990508.749511708, 9990760.9375, 9990763.050561476, 9990886.420268143, 9991059.817627193, 9991263.524638884, 9991459.58101752, 9991723.881299783, 9991843.994976593, 9991920.3125, 9992104.6875, 9992108.119142484, 9992264.70196745, 9992575.889160749, 9992619.353131874, 9992736.895012658, 9992928.168449724, 9992979.447692685, 9993005.026850259, 9993207.513846269, 9993232.8125, 9993307.422317876, 9993338.083432285, 9993484.5273462, 9993527.15688437, 9993656.25, 9993694.711612845, 9993721.59390002, 9993825.0, 9993858.76887652, 9993895.975861868, 9993953.003687872, 9993985.547273131, 9994023.88060655, 9994037.232658917, 9994090.501672875, 9994244.611784525, 9994353.25298394, 9994370.307084437, 9994378.125, 9994396.798699781, 9994398.4375, 9994422.140434366, 9994495.3125, 9994601.615951195, 9994644.597473709, 9994701.129912514, 9994760.369899498, 9994837.444660623, 9994894.164360018, 9994952.493571095, 9994965.090379981, 9995050.750730887, 9995254.777568452, 9995414.875678767, 9995434.73172232, 9995571.875, 9995835.1747051, 9995872.59743038, 9995920.917156462, 9996091.165238474, 9996122.620778972, 9996290.284067478, 9996725.0, 9996854.376378296, 9996921.324627195, 9997101.5625, 9997106.435018685, 9997114.042391866, 9997276.052033849, 9997318.326068342, 9998217.113290736, 9999420.164131673, 9999449.35758812, 9999955.071490226, 10000391.639353033, 10004228.125, 10010262.715491254, 10039639.0625, 10046313.25966646, 10052862.5, 10056390.364279294, 10056451.330323461, 10059515.625, 10061307.67890372, 10061959.250825306, 10063582.876356196, 10063601.409340363, 10064185.98605043, 10064382.560337849, 10065175.495672893, 10065393.42662671, 10065732.491640521, 10065739.160890747, 10065819.113149567, 10066163.652399134, 10066849.265524585, 10066896.934770873, 10067168.020069718, 10067555.357688386, 10067696.709877292, 10067821.795268534, 10068589.4791701, 10069055.9474011, 10069071.43725129, 10069245.910038985, 10069434.223772554, 10069527.45949128, 10070489.0625, 10070874.351849819, 10070885.679216096, 10071139.0625, 10071240.625, 10071601.020251315, 10072188.329063639, 10073552.373758392, 10074698.4375, 10074828.196189726, 10075592.705895163, 10075603.303326415, 10077473.397365302, 10077623.721591633, 10078073.855335947, 10083098.4375, 10083678.125, 10084057.840163851, 10084145.19574886, 10084167.44806772, 10084194.704202894, 10084343.532960378, 10085386.9143225, 10085447.625853777, 10085514.270569243, 10086857.435750296, 10086919.989592819, 10087332.068862079, 10088026.018983832, 10088787.5, 10089093.850596668, 10091479.195466435, 10091695.775936661, 10096618.75, 10097392.890450584, 10101815.625, 10103651.69833498, 10104093.38616116, 10104804.460814303, 10104851.5625, 10105229.572318457, 10105251.5625, 10105457.456109663, 10105639.768308673, 10105868.67873328, 10105902.22574194, 10106099.362332057, 10106349.744540399, 10106656.25, 10107463.14733648, 10107875.0, 10109033.422640847, 10109119.244667992, 10109580.70760156, 10109823.4375, 10109909.778691443, 10110187.318484277, 10110241.898935623, 10110497.655495794, 10110746.212709107, 10110924.230056213, 10111154.44214466, 10111203.280068876, 10111308.632147295, 10111374.280137293, 10111498.4375, 10111507.8125, 10111571.711850654, 10111657.8125, 10111862.5, 10112075.23490237, 10112075.377685847, 10112228.467436923, 10112382.29095917, 10112391.184264895, 10112573.63165682, 10112645.3750339, 10112792.1875, 10112876.794640781, 10113126.959686395, 10113251.5625, 10113283.856731955, 10113429.38398803, 10113574.471596982, 10114024.477661908, 10114046.999779532, 10114109.658655295, 10114147.233683009, 10114219.706106583, 10114345.207396135, 10114550.509708354, 10115045.825920217, 10115126.480962075, 10115234.460313393, 10115260.9375, 10115503.480005974, 10115796.56752703, 10116177.517830309, 10116262.635382336, 10116957.41898073, 10117485.9375, 10117835.761294806, 10118457.688804802, 10118720.96239377, 10120289.0625, 10120310.093685366, 10120597.814746432, 10120870.3125, 10146997.350923626, 10147626.30618556, 10149435.9375, 10149837.053968014, 10157785.9375, 10159394.603572048, 10159699.553988999, 10160205.94996986, 10160216.082186824, 10160489.0625, 10161382.263196547, 10161517.976007499, 10170836.508026848, 10183509.32610211, 10187556.969227718, 10188854.976091778, 10214745.3125, 10269737.5, 10285723.247905832, 10314731.492623065, 10322283.621778898, 10328345.69444538, 10330182.904631207, 10343612.372284446, 10344185.902296538, 10349583.463945003, 10350767.85369421, 10351103.97550238, 10358825.631777976, 10363223.210733505, 10383834.375, 10383903.316167818, 10388675.542464778, 10388681.25, 10453865.876860924, 10455998.555143412, 10457791.129829785, 10458467.72717686, 10458475.659819195, 10458498.914832234, 10458671.031893939, 10459163.503636247, 10459405.434902653, 10462026.035933854, 10462442.1875, 10475159.282629902, 10481757.644591527, 10482030.198632324, 10482558.347195486, 10484132.562542448, 10484289.61474013, 10485018.00854844, 10485802.180381123, 10501087.729559405, 10501792.177487034, 10517109.027022203, 10519802.32645364, 10524867.857694944, 10525643.458613945, 10546022.916404609, 10546355.973651404, 10546775.16587195, 10568174.299104834, 10599535.19204174, 10600677.340253225, 10602503.971500603, 10605445.9652286, 10621562.408395305, 10621798.09031073, 10622500.0, 10622726.95481798, 10630045.556617822, 10633609.13088851, 10635503.05269953, 10641754.761830518, 10641808.026231151, 10641841.30953866, 10660127.776480662, 10661268.393868212, 10665020.783303784, 10669052.833793676, 10670915.182281567, 10673023.401022632, 10673490.625, 10688121.097369825, 10688280.069494756, 10689694.853554057, 10689807.696382383, 10690450.86001738, 10692183.118455727, 10738619.225335443, 10740126.220666409, 10740348.86496705, 10741774.502098123, 10742623.19992798, 10746667.067536397, 10746681.726369005, 10747959.226036163, 10748542.1875, 10760357.216469266, 10760705.478093358, 10760706.12284493, 10762257.68513528, 10763692.741520805, 10764397.138284575, 10765473.280233582, 10766566.623278372, 10767637.663944054, 10768873.617257113, 10770480.912157293, 10772139.404482273, 10773106.12470411, 10801356.116821123, 10810680.31847346, 10818057.82757061, 10855675.119504781, 10866752.657224385, 10869820.272533592, 10869825.105518147, 10879636.888652379, 10889498.4375, 10957316.014677003, 10958726.999761188, 10959498.4375, 10960975.0, 10962308.527164439, 10963472.288977932, 10964507.65110128, 10967456.730959859, 10968744.25554954, 10986413.407669174, 11299365.381522806, 11301955.626901818, 11303869.577273088, 11304218.696072562, 11304447.071286488, 11304654.544997586, 11304854.6875, 11305487.5, 11305925.320724953, 11306151.320066718, 11306480.323739847, 11306571.941971555, 11307924.157067388, 11308049.508752769, 11308244.767764298, 11308260.146917308, 11308965.311162245, 11309408.932880666, 11310209.375, 11310726.030202148, 11311353.88224646, 11311949.663972778, 11312306.25, 11312491.230177203, 11312514.051928347, 11314136.948383793, 11314631.638811775, 11315104.6875, 11315527.635495024, 11315791.607703395, 11316262.113818437, 11316498.4375, 11316556.137962764, 11316639.621954117, 11316877.010217378, 11317138.209989619, 11317562.423130447, 11317742.1875, 11317797.313866088, 11317854.338882864, 11318868.296620138, 11318883.235805903, 11318892.09450813, 11318952.809462044, 11319114.0625, 11319315.625, 11319867.1875, 11320009.189201688, 11320108.498884661, 11320634.788229933, 11320771.607512483, 11321099.215910107, 11321482.471323596, 11321572.811270539, 11321632.8125, 11321682.029790265, 11321968.75, 11322571.46730476, 11322980.679594887, 11322985.9375, 11323009.10437657, 11323111.781291237, 11323589.283322692, 11324180.513131185, 11324210.202272784, 11324705.549954517, 11324733.337303597, 11325375.896355465, 11325583.484296152, 11326021.875, 11326170.981356645, 11326555.425578406, 11326578.24763215, 11326934.375, 11327150.979013652, 11328200.76864632, 11328212.993038382, 11328608.806487942, 11329303.363750942, 11330243.75, 11330381.535327647, 11330713.059969002, 11330759.143100874, 11331496.74701107, 11331866.700925088, 11332270.70851687, 11332309.885536298, 11332541.716856275, 11332583.95454559, 11332624.03006431, 11332627.27789879, 11332996.64112055, 11333486.340197692, 11333699.299344273, 11333730.746825922, 11333815.968681837, 11333909.955000104, 11334237.5, 11334249.703637168, 11334283.274657663, 11334420.057393488, 11334540.625, 11334964.009698339, 11335645.053888284, 11335749.431585792, 11335822.9727387, 11336065.774919787, 11336319.708718732, 11336481.697861573, 11336979.554605737, 11337072.860621884, 11337146.60124551, 11337338.00072358, 11337425.034825966, 11337498.4375, 11337644.776351845, 11337901.116662104, 11337903.306146994, 11338095.440329913, 11338645.32412605, 11338662.156389475, 11339012.5, 11339067.608608138, 11339285.724739525, 11339360.9375, 11339452.626922332, 11339490.924699571, 11339562.5, 11339581.697393106, 11339769.685084667, 11339853.125, 11339943.120961653, 11339944.52839517, 11340029.77130868, 11340152.895697717, 11340161.480879674, 11340293.480582705, 11340665.688704155, 11340739.0625, 11341152.17201815, 11341167.842062626, 11341202.56663367, 11341265.993473988, 11341504.745027874, 11341505.067141082, 11341590.17264421, 11341862.809718678, 11341984.06117406, 11342032.021605086, 11342241.64630901, 11342264.777063593, 11342274.826835481, 11342311.992790975, 11342432.8125, 11342614.974061603, 11343007.492282942, 11343102.234121814, 11343145.59069278, 11343154.341160363, 11343314.204841234, 11343371.875, 11343418.584700938, 11343662.967322119, 11343852.793219458, 11343990.625, 11344025.8375701, 11344203.125, 11344258.427627277, 11344298.4375, 11344323.4375, 11344761.321885636, 11344792.1875, 11344839.856737968, 11344876.5625, 11344994.58869172, 11345162.41240205, 11345370.3125, 11345594.51278295, 11345703.193269689, 11345719.304799426, 11345840.077570412, 11345943.101627449, 11345959.375, 11346131.25, 11346466.715251409, 11346707.8125, 11346767.070324406, 11347060.789394204, 11347365.982702896, 11347460.53183657, 11347510.261861237, 11347673.706304615, 11347725.0, 11348195.669537522, 11348232.378096377, 11348482.564315438, 11348830.730562022, 11348889.577241935, 11348903.125, 11348992.123475896, 11349071.739978774, 11349182.8125, 11349206.683818387, 11349341.734150752, 11349427.965271123, 11349446.729496295, 11350298.4375, 11350414.510555848, 11350623.77389559, 11350801.923099872, 11351658.773174861, 11351670.998425651, 11351731.25, 11351977.889988754, 11352292.35176889, 11352641.288735265, 11353291.287283603, 11353715.202113796, 11353857.825136425, 11354026.5625, 11354624.845397815, 11355235.9375, 11356275.029038105, 11356576.363725629, 11356862.593597626, 11357164.0625, 11358264.0625, 11358537.40964914, 11359179.6875, 11360395.126303079, 11360848.762718486, 11362070.62716068, 11363596.591489922, 11364048.342244009, 11364312.403429532, 11366046.994555661, 11366349.435707312, 11366373.4375, 11370281.832267223, 11373352.085883295], [95.58548775922932, 6.01320814014058, 26.765112195213412, 9.788181407113093, 5.678969627754286, 21.11154063307883, 90.73018489566181, 49.897952180760534, 25.50329536951411, 22.741834911686922, 8.231315455462765, 7.269258430875213, 81.9088647427757, 17.3187151941055, 11.045200535918978, 30.66751517463912, 27.815392897821855, 42.52925686639747, 12.283079923645813, 17.830656342819303, 37.84698385358875, 123.6001816899272, 41.51264996786801, 28.56746001451353, 46.86364705168485, 26.764513833847776, 78.2581655073064, 5.174284090745851, 29.773427653344168, 16.075801224118013, 10.458170872290466, 43.14641012019057, 39.749096594798566, 14.040568931586025, 73.98951768059624, 32.27034338585554, 14.55515864855923, 46.749646624558295, 53.42289563900929, 5.664157719912628, 10.68833136577545, 149.96907136114402, 46.98130250958803, 38.804211890181676, 35.0173583891157, 13.22981389641947, 51.80859772125987, 51.82465103430784, 31.781066675546914, 25.446056359838387, 17.814735506626903, 38.1774822203444, 56.86897106336503, 58.2602287757091, 7.474042082945834, 53.479486689804745, 55.73951817543015, 38.18324240259955, 24.967019355341698, 33.566006764085145, 19.199306741890528, 12.343809839921356, 9.814625594487591, 20.35658623960603, 33.32650955788597, 61.69900234816275, 19.417606529172566, 69.6825829271247, 22.2065373772513, 65.32627937681207, 36.061117113495655, 9.539722649028153, 69.05124291439263, 7.577406587036659, 10.071634792171823, 44.46341184716933, 8.590178370972566, 13.372913300220752, 23.28309007305858, 100.11373812404942, 51.01383026491613, 77.61887752503485, 46.212337902206954, 30.867307781519802, 11.175780931783569, 19.57170228450253, 78.89369102134931, 9.869290046391244, 21.296472417326097, 48.80659837119738, 6.306318633333746, 49.90161562691786, 24.922796391755746, 27.283648041756706, 11.505543740904177, 6.325095405214135, 29.019430609902617, 68.15491816240488, 65.94902961050941, 64.89453855761684, 53.61127262161861, 19.051046917230867, 29.481889498911784, 46.42234492087289, 28.815407357444617, 70.88415823893075, 110.47118024275929, 65.74147429598145, 42.6364094170515, 30.212018785359984, 12.36577024392671, 104.32107943374065, 18.30233207765047, 63.91274109107497, 24.794755863970785, 111.12885106443973, 11.578605699492982, 14.154683222497548, 64.35268703272439, 6.17910839145654, 63.01349029833551, 43.825318769784396, 5.078426671529605, 44.24566149825455, 56.03139671960262, 138.50269417923613, 96.23813617210975, 198.07239738402617, 14.696743887930015, 15.23935242686958, 14.972500645475694, 61.19026072131363, 21.637365401738528, 160.06320391432413, 41.045711239217376, 33.114224119959104, 9.17761519570235, 5.566046576201569, 13.557628486144786, 8.412478744502902, 39.35664263790514, 42.061365482654054, 32.87476692605246, 16.325497913605545, 13.820910730421968, 12.076010612986154, 35.10785848654111, 6.112756425421403, 13.325193823119935, 264.74807263510854, 12.79601653400451, 40.14764222771586, 83.9460908223016, 116.20172193394596, 74.43921570907858, 52.616671677346375, 224.84678198052134, 28.194637745983336, 19.276055945969517, 7.663769565687397, 12.107044275119144, 51.61185008389951, 26.629857361099774, 13.124792858251917, 29.515648633614475, 78.75987726144922, 123.33062458263726, 36.029687422746925, 75.55449663972391, 5.362181333375631, 31.088195404092065, 126.92449351060091, 35.861032510124005, 176.60015005288042, 27.264481560403773, 37.84110985049236, 27.6048200080266, 16.001465160736533, 90.18361792058003, 71.8734032667486, 58.26359680918382, 10.770101308431876, 35.026196825688075, 15.564104172285408, 39.85817288464503, 6.729191290195267, 66.36906532250043, 7.4966261094217685, 12.179161693437239, 20.942137764356083, 37.339336876278814, 90.24750834357965, 6.276690770610132, 75.39521299448896, 5.76446945789797, 6.774084005075873, 9.41810963680264, 16.436306260368188, 21.787419941102513, 66.54629681361676, 59.21481421731108, 106.62541280526696, 5.85865134769586, 16.80204722086571, 7.5625820223898925, 25.38565067096753, 33.23370654126729, 92.50312835932648, 8.122192309905532, 31.40507364905291, 8.494753137229836, 11.661576689890573, 9.488112687438981, 6.181703289711271, 29.173348164512305, 30.17042610418095, 38.425251737339636, 33.4520839141225, 15.907180055332487, 114.69670582323175, 51.21357437909144, 61.6031484299749, 96.74937378926643, 48.86236825031578, 5.236397153462196, 29.576542829449565, 11.036976191641012, 63.41734779131698, 29.35995614199963, 26.564578341477137, 18.524977876812027, 22.545814380852043, 31.706457136129927, 7.1795687932909615, 43.25854617750927, 75.03600460804275, 109.93453165680683, 73.91349444542976, 6.856993140199748, 76.24621826559901, 45.795743089345635, 27.560391720092053, 5.2212223837923375, 47.116742998281744, 72.18998378183838, 100.63466415232244, 7.276098116403006, 10.721200019199284, 101.42385544333827, 9.02921615668905, 26.897493166125052, 30.38868642867213, 13.207834777908795, 6.930226797753219, 64.68278419449726, 57.702700457687484, 46.894195711599025, 18.138632257571757, 33.151425054729835, 67.16275861407048, 10.491652017175737, 38.15903646208226, 48.08883722746971, 8.556574467817198, 68.35135883896513, 8.443889917671592, 23.233219970620603, 5.568059358134129, 35.912014140724324, 7.149603088663292, 17.475357913009596, 14.779538202096411, 20.8846859106485, 7.864843757600646, 6.1249045561263475, 21.9189355475237, 19.527351285713518, 9.069220041293361, 24.301429480282415, 12.706998779915972, 60.547611765134675, 5.618262478328589, 107.61871791756307, 72.23315411574018, 19.11288459556006, 102.54776417174521, 21.736156647862202, 50.67225888077736, 33.85651913419626, 21.181010519665563, 59.275636181897724, 37.086590796509356, 5.9760190035442555, 11.266195025571879, 8.380285959610037, 12.171116055817329, 10.746084677311847, 105.59855072567083, 30.210810115698635, 60.7040415939613, 32.589259070148984, 14.31013979260201, 99.01210311918832, 49.66003759874583, 122.93518456879336, 23.498936969959317, 12.053544208142327, 16.652544789868394, 10.159255604005269, 23.23785843650893, 27.22507472905824, 60.96400267947714, 39.38426921864413, 54.20304615966998, 50.763985140826506, 39.02303532532699, 16.28661535070851, 12.920814411473119, 12.868620574685957, 6.342250044523681, 125.49674680821963, 9.278251253417247, 17.613568130270536, 8.572656141268673, 21.140787086893837, 47.11218711614873, 12.113138810269797, 83.97328998623063, 139.01300278627, 56.04949557296476, 5.719996171185611, 31.238889458630997, 21.12007042858298, 50.72593412951427, 43.18730863509642, 9.604882570581308, 30.716882985835497, 10.216453890052513, 56.4812540649851, 20.506048510788954, 15.109746615864445, 35.9843912378238, 135.9979704949151, 55.89875107009698, 78.75195589152028, 36.308875491691225, 31.3656281842295, 66.19876849655732, 9.890214795466772, 41.18078433587508, 79.62144003337747, 6.708655394277812, 53.89189228736277, 6.396096098922511, 15.060218313243986, 41.02967145674825, 22.662803281844567, 27.86958923202867, 12.889665676830223, 14.491395193118471, 37.99854475407462, 5.950629174574408, 14.92090793129184, 52.076502216920126, 27.101347081706365, 70.63730989691489, 16.059914905282625, 67.17754230553906, 36.068851712370616, 8.792086688270468, 10.512769497115887, 5.8914540618017694, 11.042295654218115, 8.518566440764177, 5.320740686105919, 19.05787470946104, 145.24454428424082, 89.56547786858984, 19.320882084038857, 28.521089308763077, 51.94591047240775, 19.219051881378086, 16.186385872354094, 98.44645748190469, 17.45006624830529, 14.32915017019144, 10.633580295009876, 137.56833595382767, 13.874923312536906, 8.015503428408152, 31.01012853103994, 15.69694805849591, 22.797877385056463, 33.93925493859352, 26.301537205196865, 27.138130102837255, 35.077320820453394, 7.973761450050758, 7.748286814051039, 7.416325404286163, 61.25760425657081, 9.520511545841945, 40.85927720050594, 41.81559362738416, 82.85331128932124, 24.265791017400865, 5.125570959490278, 32.69628643760449, 45.060795463443206, 26.29047169950258, 5.464514009237976, 33.452444669066345, 12.344969207041158, 78.17332370301033, 69.49499402274404, 5.448056088717443, 7.491169318257705, 5.043205903871608, 38.91842410249734, 7.461353237116094, 19.434708462344425, 35.01865395012505, 32.357277564802914, 28.176869925346246, 64.28941478223892, 33.2438187724189, 7.761228369551284, 23.7176009968896, 13.499833596128612, 67.30577703002992, 51.863512898213585, 7.379052570959108, 42.78061405216009, 33.02101214031829, 35.842196481456746, 32.385305422371175, 57.078431816481405, 161.7281037254038, 12.729167612891814, 12.686877954553186, 60.65902570679487, 105.92133582940134, 91.66793008080981, 25.0713097564618, 176.58207556185698, 114.4495993609547, 43.91672995556732, 42.648661816976094, 73.25333899556912, 15.04800994535002, 6.252559538475853, 6.816666256715689, 10.2936674881045, 72.83785717590803, 20.678977775034696, 5.656979058769342, 5.20875549974933, 59.56930418211217, 5.036144013547315, 89.4410596320624, 8.910550662556759, 79.61249534958054, 158.953522487227, 38.89553123657973, 11.891976055232082, 32.51902827313742, 42.50173455133175, 51.768105209041416, 8.359739541032159, 18.792758365908327, 10.7909520127425, 56.516792451076626, 26.438357377844845, 28.797765205316107, 21.638547727792364, 35.770771310280026, 52.69510322082316, 18.260673353459126, 32.004226143666656, 9.12435144047788, 45.110148939382505, 6.593591838807247, 78.58023606595353, 25.536544332342814, 78.92193065172485, 10.373993851514802, 61.302791779656886, 5.303722029610844, 51.56131299212668, 46.790762187440635, 88.99318852706148, 72.30880894227006, 118.47055344671813, 31.26929724355679, 15.56744423898236, 82.38581898791122, 7.939362542543894, 58.940910914592614, 53.38361153894469, 119.63238262541284, 79.01991319109841, 54.6157482179556, 15.919406153040544, 16.60717126278108, 16.88564322394499, 47.02352888244826, 31.234239467077636, 12.011899666194765, 25.326128321192545, 78.42571378639701, 15.906863184942546, 5.413878316242108, 45.97958177036469, 29.543411563263188, 98.42420682483484, 35.379980968371264, 31.17923256835921, 25.216557353781717, 9.58334736594936, 8.48582207539188, 87.96417575924038, 5.5280011701501826, 58.27533174749955, 36.43183876523967, 27.009726848809734, 15.27697001935811, 71.15503015452435, 16.492520622570293, 107.3734503402859, 20.4867343855342, 68.157248743764, 182.2920442956025, 66.0601889106289, 25.47226777996774, 58.505783107017486, 81.73120136805935, 29.366158085604415, 31.619177054837813, 33.46372787313438, 28.95193994373235, 9.769084962223468, 21.368649730129807, 56.44857157268002, 92.78314822686065, 72.82989824442224, 24.489155650354768, 22.049123494390514, 66.89450693555474, 16.169955120572553, 54.41324920566595, 30.768457406630272, 10.100833009975094, 86.9054797155474, 55.36070078676495, 7.324283903256295, 45.31854305179026, 55.578670411829435, 72.54372313386452, 25.27149230156136, 131.6503218582221, 12.682202717497331, 11.782251008041401, 11.274833560107759, 21.322399422146734, 116.99424615061356, 8.663985970419603, 9.96630826055913, 64.53353624486206, 124.7523792763555, 10.790837936212586, 92.500615809082, 6.802281452408287, 19.284540057105303, 61.70470238007596, 26.39444736863409, 17.236326423447043, 44.79871466273643, 7.270930360500436, 68.90404064376018, 11.965891461129416, 7.324334194813069, 67.49881262634817, 5.14680346494961, 163.47564258288008, 76.13589337236677, 14.600750460898963, 9.193351011699438, 6.870028675659468, 12.080755602584505, 24.33069006786257, 17.56761288808604, 20.398174901294222, 12.253608036877603, 8.411821081548393, 28.328305173835357, 20.99847034933253, 29.406571864624915, 12.366319859910263, 81.65464948808815, 21.285212332799194, 17.63073038266638, 52.64771918123475, 19.198062660930095, 69.37996136091381, 25.438912771670708, 12.770999764020226, 9.951880963839702, 26.83568947122827, 23.16611558438287, 72.345140270634, 11.801791005427297, 40.159147240199, 59.11758679972437, 86.12883343215259, 71.35747968356017, 6.521125654130583, 5.221219882149466, 31.047729688916185, 13.738128978017885, 66.77000888528026, 5.9144404740612355, 27.16643538991025, 26.672457857994292, 28.708948422867486, 98.46497337421395, 16.62919843951085, 43.91720938313228, 10.950931840026188, 80.42849592859929, 25.21774000490931, 23.123175771779668, 51.11435363076505, 5.480417979560089, 17.525905081018315, 14.107230599485103, 8.224141269738592, 17.132286034146492, 9.208543931185579, 47.70344446213627, 59.03297498091047, 27.424505263211962, 32.958495060804424, 6.467309336384388, 7.2751225349049, 11.565783437945672, 21.67948954489315, 61.59516238869725, 8.375219133435062, 15.67398105776758, 68.76331881299625, 31.980735984186868, 14.709822313106596, 18.435105410889683, 109.59137669367658, 14.999264351983173, 56.28980978265738, 27.453195065361776, 17.539507408194122, 7.149468195061383, 46.46736579755781, 40.184085918192665, 11.796806184412539, 25.42049532496941, 102.45735908167347, 46.40520703331625, 5.868821945680249, 17.24827685447288, 52.683223064364164, 17.157811201253367, 20.823404588067127, 81.67412686537327, 24.50062980827529, 20.34715011691696, 8.778536578778262, 33.03300697173112, 68.34388407751757, 58.252814893778876, 112.12276566229554, 52.883954058297086, 67.92978988816171, 60.956247748653674, 16.27514095951107, 7.1391684319463655, 36.919492080456116, 21.580860398078933, 22.205908193750876, 24.931550920773425, 45.31653199472129, 8.9835968035357, 32.00745729517433, 21.176850823873586, 6.6127307054415425, 23.301199739960364, 138.29785540594224, 47.55823656643202, 85.09085088756143, 10.362627903095994, 5.138899690693089, 94.37581346522359, 36.18550111165656, 22.16499015027219, 129.27161890574607, 70.12184719594458, 14.531005585558676, 5.178598961409752, 115.45123448990353, 6.604462138580024, 81.59225584292057, 31.16867336193774, 5.214748963641918, 6.69584799579283, 20.526880215232783, 21.937115231024713, 7.097219609729197, 34.71690450382481, 31.529646715218522, 79.12929204666194, 29.50968445039194, 10.849201693992526, 64.29661454383167, 6.848907499410452, 14.011390419950876, 47.99501695574824, 137.23777012297907, 20.319993441084975, 59.65511042711247, 10.091033069158863, 35.61584726199297, 66.78061299889751, 20.81839160101983, 105.43718216639098, 14.001452663280446, 27.15013416921883, 18.710976137781948, 12.720981814676621, 26.176425749219916, 12.809640990324379, 6.873006365426899, 12.892937941431997, 7.865738779360716, 53.56929158231249, 18.297664094188, 14.926604661189653, 18.242438316215804, 62.07744899181844, 20.698834591709506, 12.587378855915636, 11.997001316953039, 72.20114083284584, 87.6907879172194, 49.535303951469395, 8.443455646315252, 30.77925384189016, 57.85321765818091, 7.557593056045657, 113.86429564489136, 22.984427442071848, 30.312589646906297, 69.91433367623092, 64.07142047762008, 17.521213902119353, 8.704961254820354, 31.747084376410644, 12.000699987176372, 19.70058260056357, 117.96735908866785, 5.254867144015603, 6.1505649303763485, 46.67235919656056, 32.4258540925215, 69.63705165692747, 6.061966156429703, 66.24250539716236, 27.521675195333707, 8.154395348976237, 13.142501228150062, 10.47145290066898, 15.545476288124364, 20.516742290278145, 30.00316521699822, 89.1655863490199, 20.75901348127942, 19.156241762018105, 73.2217002040056, 5.7089780152054725, 59.62651006765546, 31.943543617359417, 5.3268730825534405, 51.2153143539848, 11.33356983528869, 9.341253040666418, 7.483514153169173, 46.604079096058314, 9.714804682221267, 62.30186206231812, 42.87629386283293, 21.51229797945035, 47.78283384141179, 44.2744473674061, 25.46388188398381, 6.64439504679244, 71.5468006258162, 7.390538192955694, 19.06079285878362, 5.344801724438608, 24.697862697005466, 22.91386395661366, 43.111995734406435, 8.320554913189627, 33.061064462201195, 24.926030854373025, 18.20387650508765, 26.524144354786714, 16.560343921427915, 8.317235131842342, 9.674457726642478, 39.840771764066275, 87.8738376103924, 91.72415911110124, 8.74161871871524, 8.912199984277576, 15.286799238946447, 25.518402246944248, 41.69141167957219, 132.55597243374925, 12.930834864720376, 22.81406691005065, 52.63481115341691, 6.820346238981218, 32.06798901368958, 9.454732399542298, 37.73975498520255, 46.69377487428331, 77.96678160580464, 23.35310383909254, 67.34168947595097, 47.76360623474593, 53.90657301214392, 5.994826541823762, 55.161183763386354, 76.27979351579964, 6.549083870893113, 34.36354576044011, 7.9181949147842365, 64.08952834056191, 49.07352311084684, 10.993491515691165, 33.08664536182323, 10.072008994255157, 15.208266772590372, 83.08232848376191, 21.842763589696645, 60.444153522845774, 34.822462065763595, 37.657679340168386, 21.774475513034307, 48.34254379271037, 39.91819688998364, 25.90641389445405, 17.36252283630768, 5.0615930481436395, 27.173108426884962, 52.435442466146256, 33.88215387300342, 81.27900880248808, 5.155366586274538, 24.43597287863193, 5.378980868360515, 5.8530744537476895, 62.92176614430332, 129.0556306474168, 26.114720303349543, 43.85820173891202, 6.323609522200777, 80.45940114719376, 60.930013467182775, 23.7790416929426, 43.69631823556113, 11.908571830062808, 45.11868435813659, 14.857690501311453, 65.1255332165766, 7.314242139581658, 145.77144785234924, 24.617964669026865, 5.044323075070157, 48.69530913986638, 52.42387066076883, 13.542323634861507, 9.342110889768392, 31.305135271550906, 16.744092944203725, 48.256385416231005, 14.970547589171723, 30.625229406975112, 8.33670136818492, 50.69497731318576, 67.97346984870039, 39.25612033985715, 59.07382116471256, 25.322073977004504, 15.280623116392896, 6.121444034932263, 13.960828120380441, 50.61813338766851, 64.51275905420002, 5.620692663652788, 6.65584700129912, 38.02226174977707, 12.178669851864264, 58.865179229923214])
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);
([5650010.149478754, 5834519.560168161, 5971456.154268411, 6083502.057448806, 6083783.929331466, 6084460.160275835, 6084930.685319969, 6085101.5625, 6135571.067988711, 6223906.842998926, 6344690.519016334, 6345945.8464447465, 6364166.678359792, 6444564.91446551, 6453035.8840658935, 6457277.852698897, 6457948.400619329, 6493423.4375, 6496719.381703425, 6528702.164729113, 6573195.3125, 6577434.375, 6628519.521133806, 6629720.3125, 6643955.915323365, 6721243.654233312, 6814796.875, 6814850.379035697, 6815310.9375, 6815325.268840411, 6817760.462303553, 6818134.064159564, 6818890.16687808, 6820594.623221774, 6986460.98362148, 6994957.8125, 6995196.769803997, 7002671.992094324, 7002809.375, 7007407.121657274, 7050857.511957153, 7063616.3680581, 7072725.149314771, 7081831.25, 7119694.159932093, 7123394.776982735, 7384442.1875, 7384475.0, 7436651.048234851, 7439820.134561317, 7517861.0886544995, 7595343.453200425, 7596008.984146238, 7596285.767276969, 7596716.815584948, 7597082.8125, 7598532.8125, 7602332.216093347, 7668487.728431678, 7672903.024584496, 7692127.105163648, 7713419.705374903, 7729113.021075299, 7739262.11025961, 7819827.290910277, 7859977.539442188, 7860738.325808369, 7863337.5, 7864018.097725859, 7905474.479832735, 7906905.637523285, 7907353.3615061045, 7921457.8125, 7944801.490511613, 7954916.387447286, 7964179.6875, 7969528.658138428, 7983981.797035061, 8122887.033496739, 8126253.312186745, 8200253.125, 8202865.250079742, 8228515.625, 8233602.87081225, 8233821.905540149, 8235215.118984976, 8235481.173556576, 8236922.449291308, 8237110.417560471, 8260482.8125, 8260904.990987003, 8260927.177773141, 8268232.3104852885, 8311040.558997082, 8311956.478988514, 8312572.319912597, 8312592.577646264, 8312655.680472765, 8315403.5521477135, 8315476.125473669, 8316313.829021708, 8326244.527067874, 8328307.8125, 8328413.549167401, 8329590.625, 8333193.660836596, 8337834.515885758, 8337942.1875, 8354527.440929546, 8365472.3354978785, 8374350.66994341, 8405806.25, 8420987.266310066, 8421052.401879825, 8422068.527561206, 8423051.549649833, 8427248.708010983, 8429481.027023643, 8432340.193654068, 8586691.456939548, 8737489.0625, 8762463.524828007, 8762514.579860989, 8767907.178109517, 8797334.375, 8823213.179879434, 8826295.040119182, 8849899.75649391, 8853120.246765753, 8871692.263507934, 8887887.508107219, 8891101.36245869, 8892801.587294618, 8892942.561309246, 8893767.864124345, 9099995.982392494, 9100006.062344579, 9143124.098537764, 9144404.269933697, 9144695.543227602, 9156366.207049087, 9159982.007957745, 9160906.965342935, 9161570.63368976, 9166040.296261396, 9166138.032359548, 9166659.375, 9166762.733817402, 9166826.062188804, 9168073.830701645, 9181113.356845982, 9191728.125, 9199630.453027934, 9204932.8125, 9206657.8125, 9208170.669990703, 9218820.3125, 9218859.06625597, 9219511.413419334, 9219611.101817088, 9219939.610085592, 9220351.66341563, 9223617.711374458, 9224421.29983405, 9226549.351477593, 9231147.49240129, 9232210.635178663, 9233712.5, 9235599.277793182, 9235734.95456352, 9252806.25, 9284781.136957949, 9288838.818928365, 9291866.651848529, 9293026.495422162, 9293049.798592588, 9294226.514245996, 9296115.877822813, 9296902.477708274, 9327290.23917403, 9331209.375, 9332060.69599856, 9332346.954932164, 9343106.404287294, 9344685.9375, 9356163.888151158, 9356571.130663658, 9358762.195290742, 9358764.645003187, 9389498.681202574, 9389681.25, 9390076.007268371, 9390237.574764578, 9391071.40672936, 9394732.12015899, 9395495.365252065, 9395539.163294032, 9395598.571587518, 9396287.871441193, 9400418.75, 9400426.5625, 9412322.061584583, 9446074.02351651, 9448266.28909525, 9449813.948676571, 9450755.092652123, 9457984.478180861, 9459356.25, 9478613.738147855, 9486556.25, 9487014.436455939, 9488039.37916052, 9488436.117195396, 9488721.166956985, 9489556.25, 9489805.208458984, 9490484.99550014, 9490625.873310875, 9490656.48221062, 9491073.21967388, 9491412.019936442, 9493596.894412303, 9493777.95081446, 9494777.422491811, 9495330.93439376, 9618505.169648716, 9639570.25366738, 9655098.288104456, 9674640.625, 9713300.869181545, 9728658.515217142, 9747364.410821691, 9750892.1875, 9756418.966568131, 9760202.441081824, 9760915.625, 9761201.5625, 9761340.654041583, 9762650.450902902, 9763651.5625, 9763932.8125, 9764363.94284901, 9764444.434632944, 9766351.560765345, 9767585.58696794, 9768889.716021651, 9784950.557352988, 9788562.223371498, 9793321.246096184, 9794737.119604545, 9796882.749256954, 9799859.220631717, 9799963.51876261, 9801450.781419005, 9801798.237598484, 9801915.625, 9805306.25, 9843769.394504279, 9858477.743810415, 9876021.875, 9940694.863522, 9942585.9375, 9942735.149072006, 9943527.040403629, 9943626.691930879, 9944053.411087243, 9944384.321161749, 9944505.38775203, 9945035.9375, 9945531.356021112, 9947254.836832307, 9949030.233178798, 9954179.124266963, 9954194.419675427, 9954239.409144122, 9954462.041654106, 9954463.339943308, 9961524.672751375, 9961790.36692223, 9961827.382252315, 9964817.1875, 9966451.638499176, 9966659.375, 9967326.5625, 9967632.103537498, 9968206.25, 9968290.205549547, 9968314.0625, 9968319.125107434, 9968752.59313484, 9969400.0, 9969421.875, 9969430.0985406, 9969941.164164066, 9970297.048098054, 9970311.516139532, 9970312.24131377, 9970342.1875, 9970446.431130908, 9970488.557121152, 9970617.083615229, 9970645.1363668, 9970745.375105875, 9970923.4375, 9971308.85948774, 9971756.278095528, 9972339.410111476, 9972407.818205886, 9972458.652366625, 9972596.705590542, 9974597.88905751, 9981389.0625, 9982980.73652565, 9985831.25, 9987032.8125, 9988915.625, 9988941.805764593, 9989247.582947334, 9989849.324026976, 9989858.384305863, 9990108.244304767, 9990123.373963622, 9990475.81743776, 9990483.316946441, 9990508.749511708, 9990760.9375, 9990763.050561476, 9990886.420268143, 9991059.817627193, 9991263.524638884, 9991459.58101752, 9991723.881299783, 9991843.994976593, 9991920.3125, 9992104.6875, 9992108.119142484, 9992264.70196745, 9992575.889160749, 9992619.353131874, 9992736.895012658, 9992928.168449724, 9992979.447692685, 9993005.026850259, 9993207.513846269, 9993232.8125, 9993307.422317876, 9993338.083432285, 9993484.5273462, 9993527.15688437, 9993656.25, 9993694.711612845, 9993721.59390002, 9993825.0, 9993858.76887652, 9993895.975861868, 9993953.003687872, 9993985.547273131, 9994023.88060655, 9994037.232658917, 9994090.501672875, 9994244.611784525, 9994353.25298394, 9994370.307084437, 9994378.125, 9994396.798699781, 9994398.4375, 9994422.140434366, 9994495.3125, 9994601.615951195, 9994644.597473709, 9994701.129912514, 9994760.369899498, 9994837.444660623, 9994894.164360018, 9994952.493571095, 9994965.090379981, 9995050.750730887, 9995254.777568452, 9995414.875678767, 9995434.73172232, 9995571.875, 9995835.1747051, 9995872.59743038, 9995920.917156462, 9996091.165238474, 9996122.620778972, 9996290.284067478, 9996725.0, 9996854.376378296, 9996921.324627195, 9997101.5625, 9997106.435018685, 9997114.042391866, 9997276.052033849, 9997318.326068342, 9998217.113290736, 9999420.164131673, 9999449.35758812, 9999955.071490226, 10000391.639353033, 10004228.125, 10010262.715491254, 10039639.0625, 10046313.25966646, 10052862.5, 10056390.364279294, 10056451.330323461, 10059515.625, 10061307.67890372, 10061959.250825306, 10063582.876356196, 10063601.409340363, 10064185.98605043, 10064382.560337849, 10065175.495672893, 10065393.42662671, 10065732.491640521, 10065739.160890747, 10065819.113149567, 10066163.652399134, 10066849.265524585, 10066896.934770873, 10067168.020069718, 10067555.357688386, 10067696.709877292, 10067821.795268534, 10068589.4791701, 10069055.9474011, 10069071.43725129, 10069245.910038985, 10069434.223772554, 10069527.45949128, 10070489.0625, 10070874.351849819, 10070885.679216096, 10071139.0625, 10071240.625, 10071601.020251315, 10072188.329063639, 10073552.373758392, 10074698.4375, 10074828.196189726, 10075592.705895163, 10075603.303326415, 10077473.397365302, 10077623.721591633, 10078073.855335947, 10083098.4375, 10083678.125, 10084057.840163851, 10084145.19574886, 10084167.44806772, 10084194.704202894, 10084343.532960378, 10085386.9143225, 10085447.625853777, 10085514.270569243, 10086857.435750296, 10086919.989592819, 10087332.068862079, 10088026.018983832, 10088787.5, 10089093.850596668, 10091479.195466435, 10091695.775936661, 10096618.75, 10097392.890450584, 10101815.625, 10103651.69833498, 10104093.38616116, 10104804.460814303, 10104851.5625, 10105229.572318457, 10105251.5625, 10105457.456109663, 10105639.768308673, 10105868.67873328, 10105902.22574194, 10106099.362332057, 10106349.744540399, 10106656.25, 10107463.14733648, 10107875.0, 10109033.422640847, 10109119.244667992, 10109580.70760156, 10109823.4375, 10109909.778691443, 10110187.318484277, 10110241.898935623, 10110497.655495794, 10110746.212709107, 10110924.230056213, 10111154.44214466, 10111203.280068876, 10111308.632147295, 10111374.280137293, 10111498.4375, 10111507.8125, 10111571.711850654, 10111657.8125, 10111862.5, 10112075.23490237, 10112075.377685847, 10112228.467436923, 10112382.29095917, 10112391.184264895, 10112573.63165682, 10112645.3750339, 10112792.1875, 10112876.794640781, 10113126.959686395, 10113251.5625, 10113283.856731955, 10113429.38398803, 10113574.471596982, 10114024.477661908, 10114046.999779532, 10114109.658655295, 10114147.233683009, 10114219.706106583, 10114345.207396135, 10114550.509708354, 10115045.825920217, 10115126.480962075, 10115234.460313393, 10115260.9375, 10115503.480005974, 10115796.56752703, 10116177.517830309, 10116262.635382336, 10116957.41898073, 10117485.9375, 10117835.761294806, 10118457.688804802, 10118720.96239377, 10120289.0625, 10120310.093685366, 10120597.814746432, 10120870.3125, 10146997.350923626, 10147626.30618556, 10149435.9375, 10149837.053968014, 10157785.9375, 10159394.603572048, 10159699.553988999, 10160205.94996986, 10160216.082186824, 10160489.0625, 10161382.263196547, 10161517.976007499, 10170836.508026848, 10183509.32610211, 10187556.969227718, 10188854.976091778, 10214745.3125, 10269737.5, 10285723.247905832, 10314731.492623065, 10322283.621778898, 10328345.69444538, 10330182.904631207, 10343612.372284446, 10344185.902296538, 10349583.463945003, 10350767.85369421, 10351103.97550238, 10358825.631777976, 10363223.210733505, 10383834.375, 10383903.316167818, 10388675.542464778, 10388681.25, 10453865.876860924, 10455998.555143412, 10457791.129829785, 10458467.72717686, 10458475.659819195, 10458498.914832234, 10458671.031893939, 10459163.503636247, 10459405.434902653, 10462026.035933854, 10462442.1875, 10475159.282629902, 10481757.644591527, 10482030.198632324, 10482558.347195486, 10484132.562542448, 10484289.61474013, 10485018.00854844, 10485802.180381123, 10501087.729559405, 10501792.177487034, 10517109.027022203, 10519802.32645364, 10524867.857694944, 10525643.458613945, 10546022.916404609, 10546355.973651404, 10546775.16587195, 10568174.299104834, 10599535.19204174, 10600677.340253225, 10602503.971500603, 10605445.9652286, 10621562.408395305, 10621798.09031073, 10622500.0, 10622726.95481798, 10630045.556617822, 10633609.13088851, 10635503.05269953, 10641754.761830518, 10641808.026231151, 10641841.30953866, 10660127.776480662, 10661268.393868212, 10665020.783303784, 10669052.833793676, 10670915.182281567, 10673023.401022632, 10673490.625, 10688121.097369825, 10688280.069494756, 10689694.853554057, 10689807.696382383, 10690450.86001738, 10692183.118455727, 10738619.225335443, 10740126.220666409, 10740348.86496705, 10741774.502098123, 10742623.19992798, 10746667.067536397, 10746681.726369005, 10747959.226036163, 10748542.1875, 10760357.216469266, 10760705.478093358, 10760706.12284493, 10762257.68513528, 10763692.741520805, 10764397.138284575, 10765473.280233582, 10766566.623278372, 10767637.663944054, 10768873.617257113, 10770480.912157293, 10772139.404482273, 10773106.12470411, 10801356.116821123, 10810680.31847346, 10818057.82757061, 10855675.119504781, 10866752.657224385, 10869820.272533592, 10869825.105518147, 10879636.888652379, 10889498.4375, 10957316.014677003, 10958726.999761188, 10959498.4375, 10960975.0, 10962308.527164439, 10963472.288977932, 10964507.65110128, 10967456.730959859, 10968744.25554954, 10986413.407669174, 11299365.381522806, 11301955.626901818, 11303869.577273088, 11304218.696072562, 11304447.071286488, 11304654.544997586, 11304854.6875, 11305487.5, 11305925.320724953, 11306151.320066718, 11306480.323739847, 11306571.941971555, 11307924.157067388, 11308049.508752769, 11308244.767764298, 11308260.146917308, 11308965.311162245, 11309408.932880666, 11310209.375, 11310726.030202148, 11311353.88224646, 11311949.663972778, 11312306.25, 11312491.230177203, 11312514.051928347, 11314136.948383793, 11314631.638811775, 11315104.6875, 11315527.635495024, 11315791.607703395, 11316262.113818437, 11316498.4375, 11316556.137962764, 11316639.621954117, 11316877.010217378, 11317138.209989619, 11317562.423130447, 11317742.1875, 11317797.313866088, 11317854.338882864, 11318868.296620138, 11318883.235805903, 11318892.09450813, 11318952.809462044, 11319114.0625, 11319315.625, 11319867.1875, 11320009.189201688, 11320108.498884661, 11320634.788229933, 11320771.607512483, 11321099.215910107, 11321482.471323596, 11321572.811270539, 11321632.8125, 11321682.029790265, 11321968.75, 11322571.46730476, 11322980.679594887, 11322985.9375, 11323009.10437657, 11323111.781291237, 11323589.283322692, 11324180.513131185, 11324210.202272784, 11324705.549954517, 11324733.337303597, 11325375.896355465, 11325583.484296152, 11326021.875, 11326170.981356645, 11326555.425578406, 11326578.24763215, 11326934.375, 11327150.979013652, 11328200.76864632, 11328212.993038382, 11328608.806487942, 11329303.363750942, 11330243.75, 11330381.535327647, 11330713.059969002, 11330759.143100874, 11331496.74701107, 11331866.700925088, 11332270.70851687, 11332309.885536298, 11332541.716856275, 11332583.95454559, 11332624.03006431, 11332627.27789879, 11332996.64112055, 11333486.340197692, 11333699.299344273, 11333730.746825922, 11333815.968681837, 11333909.955000104, 11334237.5, 11334249.703637168, 11334283.274657663, 11334420.057393488, 11334540.625, 11334964.009698339, 11335645.053888284, 11335749.431585792, 11335822.9727387, 11336065.774919787, 11336319.708718732, 11336481.697861573, 11336979.554605737, 11337072.860621884, 11337146.60124551, 11337338.00072358, 11337425.034825966, 11337498.4375, 11337644.776351845, 11337901.116662104, 11337903.306146994, 11338095.440329913, 11338645.32412605, 11338662.156389475, 11339012.5, 11339067.608608138, 11339285.724739525, 11339360.9375, 11339452.626922332, 11339490.924699571, 11339562.5, 11339581.697393106, 11339769.685084667, 11339853.125, 11339943.120961653, 11339944.52839517, 11340029.77130868, 11340152.895697717, 11340161.480879674, 11340293.480582705, 11340665.688704155, 11340739.0625, 11341152.17201815, 11341167.842062626, 11341202.56663367, 11341265.993473988, 11341504.745027874, 11341505.067141082, 11341590.17264421, 11341862.809718678, 11341984.06117406, 11342032.021605086, 11342241.64630901, 11342264.777063593, 11342274.826835481, 11342311.992790975, 11342432.8125, 11342614.974061603, 11343007.492282942, 11343102.234121814, 11343145.59069278, 11343154.341160363, 11343314.204841234, 11343371.875, 11343418.584700938, 11343662.967322119, 11343852.793219458, 11343990.625, 11344025.8375701, 11344203.125, 11344258.427627277, 11344298.4375, 11344323.4375, 11344761.321885636, 11344792.1875, 11344839.856737968, 11344876.5625, 11344994.58869172, 11345162.41240205, 11345370.3125, 11345594.51278295, 11345703.193269689, 11345719.304799426, 11345840.077570412, 11345943.101627449, 11345959.375, 11346131.25, 11346466.715251409, 11346707.8125, 11346767.070324406, 11347060.789394204, 11347365.982702896, 11347460.53183657, 11347510.261861237, 11347673.706304615, 11347725.0, 11348195.669537522, 11348232.378096377, 11348482.564315438, 11348830.730562022, 11348889.577241935, 11348903.125, 11348992.123475896, 11349071.739978774, 11349182.8125, 11349206.683818387, 11349341.734150752, 11349427.965271123, 11349446.729496295, 11350298.4375, 11350414.510555848, 11350623.77389559, 11350801.923099872, 11351658.773174861, 11351670.998425651, 11351731.25, 11351977.889988754, 11352292.35176889, 11352641.288735265, 11353291.287283603, 11353715.202113796, 11353857.825136425, 11354026.5625, 11354624.845397815, 11355235.9375, 11356275.029038105, 11356576.363725629, 11356862.593597626, 11357164.0625, 11358264.0625, 11358537.40964914, 11359179.6875, 11360395.126303079, 11360848.762718486, 11362070.62716068, 11363596.591489922, 11364048.342244009, 11364312.403429532, 11366046.994555661, 11366349.435707312, 11366373.4375, 11370281.832267223, 11373352.085883295], [95.58548775922932, 6.01320814014058, 26.765112195213412, 9.788181407113093, 5.678969627754286, 21.11154063307883, 90.73018489566181, 49.897952180760534, 25.50329536951411, 22.741834911686922, 8.231315455462765, 7.269258430875213, 81.9088647427757, 17.3187151941055, 11.045200535918978, 30.66751517463912, 27.815392897821855, 42.52925686639747, 12.283079923645813, 17.830656342819303, 37.84698385358875, 123.6001816899272, 41.51264996786801, 28.56746001451353, 46.86364705168485, 26.764513833847776, 78.2581655073064, 5.174284090745851, 29.773427653344168, 16.075801224118013, 10.458170872290466, 43.14641012019057, 39.749096594798566, 14.040568931586025, 73.98951768059624, 32.27034338585554, 14.55515864855923, 46.749646624558295, 53.42289563900929, 5.664157719912628, 10.68833136577545, 149.96907136114402, 46.98130250958803, 38.804211890181676, 35.0173583891157, 13.22981389641947, 51.80859772125987, 51.82465103430784, 31.781066675546914, 25.446056359838387, 17.814735506626903, 38.1774822203444, 56.86897106336503, 58.2602287757091, 7.474042082945834, 53.479486689804745, 55.73951817543015, 38.18324240259955, 24.967019355341698, 33.566006764085145, 19.199306741890528, 12.343809839921356, 9.814625594487591, 20.35658623960603, 33.32650955788597, 61.69900234816275, 19.417606529172566, 69.6825829271247, 22.2065373772513, 65.32627937681207, 36.061117113495655, 9.539722649028153, 69.05124291439263, 7.577406587036659, 10.071634792171823, 44.46341184716933, 8.590178370972566, 13.372913300220752, 23.28309007305858, 100.11373812404942, 51.01383026491613, 77.61887752503485, 46.212337902206954, 30.867307781519802, 11.175780931783569, 19.57170228450253, 78.89369102134931, 9.869290046391244, 21.296472417326097, 48.80659837119738, 6.306318633333746, 49.90161562691786, 24.922796391755746, 27.283648041756706, 11.505543740904177, 6.325095405214135, 29.019430609902617, 68.15491816240488, 65.94902961050941, 64.89453855761684, 53.61127262161861, 19.051046917230867, 29.481889498911784, 46.42234492087289, 28.815407357444617, 70.88415823893075, 110.47118024275929, 65.74147429598145, 42.6364094170515, 30.212018785359984, 12.36577024392671, 104.32107943374065, 18.30233207765047, 63.91274109107497, 24.794755863970785, 111.12885106443973, 11.578605699492982, 14.154683222497548, 64.35268703272439, 6.17910839145654, 63.01349029833551, 43.825318769784396, 5.078426671529605, 44.24566149825455, 56.03139671960262, 138.50269417923613, 96.23813617210975, 198.07239738402617, 14.696743887930015, 15.23935242686958, 14.972500645475694, 61.19026072131363, 21.637365401738528, 160.06320391432413, 41.045711239217376, 33.114224119959104, 9.17761519570235, 5.566046576201569, 13.557628486144786, 8.412478744502902, 39.35664263790514, 42.061365482654054, 32.87476692605246, 16.325497913605545, 13.820910730421968, 12.076010612986154, 35.10785848654111, 6.112756425421403, 13.325193823119935, 264.74807263510854, 12.79601653400451, 40.14764222771586, 83.9460908223016, 116.20172193394596, 74.43921570907858, 52.616671677346375, 224.84678198052134, 28.194637745983336, 19.276055945969517, 7.663769565687397, 12.107044275119144, 51.61185008389951, 26.629857361099774, 13.124792858251917, 29.515648633614475, 78.75987726144922, 123.33062458263726, 36.029687422746925, 75.55449663972391, 5.362181333375631, 31.088195404092065, 126.92449351060091, 35.861032510124005, 176.60015005288042, 27.264481560403773, 37.84110985049236, 27.6048200080266, 16.001465160736533, 90.18361792058003, 71.8734032667486, 58.26359680918382, 10.770101308431876, 35.026196825688075, 15.564104172285408, 39.85817288464503, 6.729191290195267, 66.36906532250043, 7.4966261094217685, 12.179161693437239, 20.942137764356083, 37.339336876278814, 90.24750834357965, 6.276690770610132, 75.39521299448896, 5.76446945789797, 6.774084005075873, 9.41810963680264, 16.436306260368188, 21.787419941102513, 66.54629681361676, 59.21481421731108, 106.62541280526696, 5.85865134769586, 16.80204722086571, 7.5625820223898925, 25.38565067096753, 33.23370654126729, 92.50312835932648, 8.122192309905532, 31.40507364905291, 8.494753137229836, 11.661576689890573, 9.488112687438981, 6.181703289711271, 29.173348164512305, 30.17042610418095, 38.425251737339636, 33.4520839141225, 15.907180055332487, 114.69670582323175, 51.21357437909144, 61.6031484299749, 96.74937378926643, 48.86236825031578, 5.236397153462196, 29.576542829449565, 11.036976191641012, 63.41734779131698, 29.35995614199963, 26.564578341477137, 18.524977876812027, 22.545814380852043, 31.706457136129927, 7.1795687932909615, 43.25854617750927, 75.03600460804275, 109.93453165680683, 73.91349444542976, 6.856993140199748, 76.24621826559901, 45.795743089345635, 27.560391720092053, 5.2212223837923375, 47.116742998281744, 72.18998378183838, 100.63466415232244, 7.276098116403006, 10.721200019199284, 101.42385544333827, 9.02921615668905, 26.897493166125052, 30.38868642867213, 13.207834777908795, 6.930226797753219, 64.68278419449726, 57.702700457687484, 46.894195711599025, 18.138632257571757, 33.151425054729835, 67.16275861407048, 10.491652017175737, 38.15903646208226, 48.08883722746971, 8.556574467817198, 68.35135883896513, 8.443889917671592, 23.233219970620603, 5.568059358134129, 35.912014140724324, 7.149603088663292, 17.475357913009596, 14.779538202096411, 20.8846859106485, 7.864843757600646, 6.1249045561263475, 21.9189355475237, 19.527351285713518, 9.069220041293361, 24.301429480282415, 12.706998779915972, 60.547611765134675, 5.618262478328589, 107.61871791756307, 72.23315411574018, 19.11288459556006, 102.54776417174521, 21.736156647862202, 50.67225888077736, 33.85651913419626, 21.181010519665563, 59.275636181897724, 37.086590796509356, 5.9760190035442555, 11.266195025571879, 8.380285959610037, 12.171116055817329, 10.746084677311847, 105.59855072567083, 30.210810115698635, 60.7040415939613, 32.589259070148984, 14.31013979260201, 99.01210311918832, 49.66003759874583, 122.93518456879336, 23.498936969959317, 12.053544208142327, 16.652544789868394, 10.159255604005269, 23.23785843650893, 27.22507472905824, 60.96400267947714, 39.38426921864413, 54.20304615966998, 50.763985140826506, 39.02303532532699, 16.28661535070851, 12.920814411473119, 12.868620574685957, 6.342250044523681, 125.49674680821963, 9.278251253417247, 17.613568130270536, 8.572656141268673, 21.140787086893837, 47.11218711614873, 12.113138810269797, 83.97328998623063, 139.01300278627, 56.04949557296476, 5.719996171185611, 31.238889458630997, 21.12007042858298, 50.72593412951427, 43.18730863509642, 9.604882570581308, 30.716882985835497, 10.216453890052513, 56.4812540649851, 20.506048510788954, 15.109746615864445, 35.9843912378238, 135.9979704949151, 55.89875107009698, 78.75195589152028, 36.308875491691225, 31.3656281842295, 66.19876849655732, 9.890214795466772, 41.18078433587508, 79.62144003337747, 6.708655394277812, 53.89189228736277, 6.396096098922511, 15.060218313243986, 41.02967145674825, 22.662803281844567, 27.86958923202867, 12.889665676830223, 14.491395193118471, 37.99854475407462, 5.950629174574408, 14.92090793129184, 52.076502216920126, 27.101347081706365, 70.63730989691489, 16.059914905282625, 67.17754230553906, 36.068851712370616, 8.792086688270468, 10.512769497115887, 5.8914540618017694, 11.042295654218115, 8.518566440764177, 5.320740686105919, 19.05787470946104, 145.24454428424082, 89.56547786858984, 19.320882084038857, 28.521089308763077, 51.94591047240775, 19.219051881378086, 16.186385872354094, 98.44645748190469, 17.45006624830529, 14.32915017019144, 10.633580295009876, 137.56833595382767, 13.874923312536906, 8.015503428408152, 31.01012853103994, 15.69694805849591, 22.797877385056463, 33.93925493859352, 26.301537205196865, 27.138130102837255, 35.077320820453394, 7.973761450050758, 7.748286814051039, 7.416325404286163, 61.25760425657081, 9.520511545841945, 40.85927720050594, 41.81559362738416, 82.85331128932124, 24.265791017400865, 5.125570959490278, 32.69628643760449, 45.060795463443206, 26.29047169950258, 5.464514009237976, 33.452444669066345, 12.344969207041158, 78.17332370301033, 69.49499402274404, 5.448056088717443, 7.491169318257705, 5.043205903871608, 38.91842410249734, 7.461353237116094, 19.434708462344425, 35.01865395012505, 32.357277564802914, 28.176869925346246, 64.28941478223892, 33.2438187724189, 7.761228369551284, 23.7176009968896, 13.499833596128612, 67.30577703002992, 51.863512898213585, 7.379052570959108, 42.78061405216009, 33.02101214031829, 35.842196481456746, 32.385305422371175, 57.078431816481405, 161.7281037254038, 12.729167612891814, 12.686877954553186, 60.65902570679487, 105.92133582940134, 91.66793008080981, 25.0713097564618, 176.58207556185698, 114.4495993609547, 43.91672995556732, 42.648661816976094, 73.25333899556912, 15.04800994535002, 6.252559538475853, 6.816666256715689, 10.2936674881045, 72.83785717590803, 20.678977775034696, 5.656979058769342, 5.20875549974933, 59.56930418211217, 5.036144013547315, 89.4410596320624, 8.910550662556759, 79.61249534958054, 158.953522487227, 38.89553123657973, 11.891976055232082, 32.51902827313742, 42.50173455133175, 51.768105209041416, 8.359739541032159, 18.792758365908327, 10.7909520127425, 56.516792451076626, 26.438357377844845, 28.797765205316107, 21.638547727792364, 35.770771310280026, 52.69510322082316, 18.260673353459126, 32.004226143666656, 9.12435144047788, 45.110148939382505, 6.593591838807247, 78.58023606595353, 25.536544332342814, 78.92193065172485, 10.373993851514802, 61.302791779656886, 5.303722029610844, 51.56131299212668, 46.790762187440635, 88.99318852706148, 72.30880894227006, 118.47055344671813, 31.26929724355679, 15.56744423898236, 82.38581898791122, 7.939362542543894, 58.940910914592614, 53.38361153894469, 119.63238262541284, 79.01991319109841, 54.6157482179556, 15.919406153040544, 16.60717126278108, 16.88564322394499, 47.02352888244826, 31.234239467077636, 12.011899666194765, 25.326128321192545, 78.42571378639701, 15.906863184942546, 5.413878316242108, 45.97958177036469, 29.543411563263188, 98.42420682483484, 35.379980968371264, 31.17923256835921, 25.216557353781717, 9.58334736594936, 8.48582207539188, 87.96417575924038, 5.5280011701501826, 58.27533174749955, 36.43183876523967, 27.009726848809734, 15.27697001935811, 71.15503015452435, 16.492520622570293, 107.3734503402859, 20.4867343855342, 68.157248743764, 182.2920442956025, 66.0601889106289, 25.47226777996774, 58.505783107017486, 81.73120136805935, 29.366158085604415, 31.619177054837813, 33.46372787313438, 28.95193994373235, 9.769084962223468, 21.368649730129807, 56.44857157268002, 92.78314822686065, 72.82989824442224, 24.489155650354768, 22.049123494390514, 66.89450693555474, 16.169955120572553, 54.41324920566595, 30.768457406630272, 10.100833009975094, 86.9054797155474, 55.36070078676495, 7.324283903256295, 45.31854305179026, 55.578670411829435, 72.54372313386452, 25.27149230156136, 131.6503218582221, 12.682202717497331, 11.782251008041401, 11.274833560107759, 21.322399422146734, 116.99424615061356, 8.663985970419603, 9.96630826055913, 64.53353624486206, 124.7523792763555, 10.790837936212586, 92.500615809082, 6.802281452408287, 19.284540057105303, 61.70470238007596, 26.39444736863409, 17.236326423447043, 44.79871466273643, 7.270930360500436, 68.90404064376018, 11.965891461129416, 7.324334194813069, 67.49881262634817, 5.14680346494961, 163.47564258288008, 76.13589337236677, 14.600750460898963, 9.193351011699438, 6.870028675659468, 12.080755602584505, 24.33069006786257, 17.56761288808604, 20.398174901294222, 12.253608036877603, 8.411821081548393, 28.328305173835357, 20.99847034933253, 29.406571864624915, 12.366319859910263, 81.65464948808815, 21.285212332799194, 17.63073038266638, 52.64771918123475, 19.198062660930095, 69.37996136091381, 25.438912771670708, 12.770999764020226, 9.951880963839702, 26.83568947122827, 23.16611558438287, 72.345140270634, 11.801791005427297, 40.159147240199, 59.11758679972437, 86.12883343215259, 71.35747968356017, 6.521125654130583, 5.221219882149466, 31.047729688916185, 13.738128978017885, 66.77000888528026, 5.9144404740612355, 27.16643538991025, 26.672457857994292, 28.708948422867486, 98.46497337421395, 16.62919843951085, 43.91720938313228, 10.950931840026188, 80.42849592859929, 25.21774000490931, 23.123175771779668, 51.11435363076505, 5.480417979560089, 17.525905081018315, 14.107230599485103, 8.224141269738592, 17.132286034146492, 9.208543931185579, 47.70344446213627, 59.03297498091047, 27.424505263211962, 32.958495060804424, 6.467309336384388, 7.2751225349049, 11.565783437945672, 21.67948954489315, 61.59516238869725, 8.375219133435062, 15.67398105776758, 68.76331881299625, 31.980735984186868, 14.709822313106596, 18.435105410889683, 109.59137669367658, 14.999264351983173, 56.28980978265738, 27.453195065361776, 17.539507408194122, 7.149468195061383, 46.46736579755781, 40.184085918192665, 11.796806184412539, 25.42049532496941, 102.45735908167347, 46.40520703331625, 5.868821945680249, 17.24827685447288, 52.683223064364164, 17.157811201253367, 20.823404588067127, 81.67412686537327, 24.50062980827529, 20.34715011691696, 8.778536578778262, 33.03300697173112, 68.34388407751757, 58.252814893778876, 112.12276566229554, 52.883954058297086, 67.92978988816171, 60.956247748653674, 16.27514095951107, 7.1391684319463655, 36.919492080456116, 21.580860398078933, 22.205908193750876, 24.931550920773425, 45.31653199472129, 8.9835968035357, 32.00745729517433, 21.176850823873586, 6.6127307054415425, 23.301199739960364, 138.29785540594224, 47.55823656643202, 85.09085088756143, 10.362627903095994, 5.138899690693089, 94.37581346522359, 36.18550111165656, 22.16499015027219, 129.27161890574607, 70.12184719594458, 14.531005585558676, 5.178598961409752, 115.45123448990353, 6.604462138580024, 81.59225584292057, 31.16867336193774, 5.214748963641918, 6.69584799579283, 20.526880215232783, 21.937115231024713, 7.097219609729197, 34.71690450382481, 31.529646715218522, 79.12929204666194, 29.50968445039194, 10.849201693992526, 64.29661454383167, 6.848907499410452, 14.011390419950876, 47.99501695574824, 137.23777012297907, 20.319993441084975, 59.65511042711247, 10.091033069158863, 35.61584726199297, 66.78061299889751, 20.81839160101983, 105.43718216639098, 14.001452663280446, 27.15013416921883, 18.710976137781948, 12.720981814676621, 26.176425749219916, 12.809640990324379, 6.873006365426899, 12.892937941431997, 7.865738779360716, 53.56929158231249, 18.297664094188, 14.926604661189653, 18.242438316215804, 62.07744899181844, 20.698834591709506, 12.587378855915636, 11.997001316953039, 72.20114083284584, 87.6907879172194, 49.535303951469395, 8.443455646315252, 30.77925384189016, 57.85321765818091, 7.557593056045657, 113.86429564489136, 22.984427442071848, 30.312589646906297, 69.91433367623092, 64.07142047762008, 17.521213902119353, 8.704961254820354, 31.747084376410644, 12.000699987176372, 19.70058260056357, 117.96735908866785, 5.254867144015603, 6.1505649303763485, 46.67235919656056, 32.4258540925215, 69.63705165692747, 6.061966156429703, 66.24250539716236, 27.521675195333707, 8.154395348976237, 13.142501228150062, 10.47145290066898, 15.545476288124364, 20.516742290278145, 30.00316521699822, 89.1655863490199, 20.75901348127942, 19.156241762018105, 73.2217002040056, 5.7089780152054725, 59.62651006765546, 31.943543617359417, 5.3268730825534405, 51.2153143539848, 11.33356983528869, 9.341253040666418, 7.483514153169173, 46.604079096058314, 9.714804682221267, 62.30186206231812, 42.87629386283293, 21.51229797945035, 47.78283384141179, 44.2744473674061, 25.46388188398381, 6.64439504679244, 71.5468006258162, 7.390538192955694, 19.06079285878362, 5.344801724438608, 24.697862697005466, 22.91386395661366, 43.111995734406435, 8.320554913189627, 33.061064462201195, 24.926030854373025, 18.20387650508765, 26.524144354786714, 16.560343921427915, 8.317235131842342, 9.674457726642478, 39.840771764066275, 87.8738376103924, 91.72415911110124, 8.74161871871524, 8.912199984277576, 15.286799238946447, 25.518402246944248, 41.69141167957219, 132.55597243374925, 12.930834864720376, 22.81406691005065, 52.63481115341691, 6.820346238981218, 32.06798901368958, 9.454732399542298, 37.73975498520255, 46.69377487428331, 77.96678160580464, 23.35310383909254, 67.34168947595097, 47.76360623474593, 53.90657301214392, 5.994826541823762, 55.161183763386354, 76.27979351579964, 6.549083870893113, 34.36354576044011, 7.9181949147842365, 64.08952834056191, 49.07352311084684, 10.993491515691165, 33.08664536182323, 10.072008994255157, 15.208266772590372, 83.08232848376191, 21.842763589696645, 60.444153522845774, 34.822462065763595, 37.657679340168386, 21.774475513034307, 48.34254379271037, 39.91819688998364, 25.90641389445405, 17.36252283630768, 5.0615930481436395, 27.173108426884962, 52.435442466146256, 33.88215387300342, 81.27900880248808, 5.155366586274538, 24.43597287863193, 5.378980868360515, 5.8530744537476895, 62.92176614430332, 129.0556306474168, 26.114720303349543, 43.85820173891202, 6.323609522200777, 80.45940114719376, 60.930013467182775, 23.7790416929426, 43.69631823556113, 11.908571830062808, 45.11868435813659, 14.857690501311453, 65.1255332165766, 7.314242139581658, 145.77144785234924, 24.617964669026865, 5.044323075070157, 48.69530913986638, 52.42387066076883, 13.542323634861507, 9.342110889768392, 31.305135271550906, 16.744092944203725, 48.256385416231005, 14.970547589171723, 30.625229406975112, 8.33670136818492, 50.69497731318576, 67.97346984870039, 39.25612033985715, 59.07382116471256, 25.322073977004504, 15.280623116392896, 6.121444034932263, 13.960828120380441, 50.61813338766851, 64.51275905420002, 5.620692663652788, 6.65584700129912, 38.02226174977707, 12.178669851864264, 58.865179229923214])
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)