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 = 44928
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);
([6545529.213785079, 6601443.740339845, 6662226.602305766, 6732844.522830824, 6817765.543901257, 6907381.143830056, 7325765.2641742155, 7759829.056142265, 7802295.032793387, 7834572.64073461, 7868227.138216119, 7933558.791688724, 8032593.463836384, 8073361.591505273, 8092861.36892458, 8105679.057481564, 8106128.95183724, 8106993.75, 8108664.7105028825, 8112257.954749868, 8113650.262079768, 8113853.125, 8117096.369119483, 8118103.11295699, 8192141.560404632, 8194721.875, 8206549.852044636, 8208863.197536281, 8209642.047681872, 8211166.8731837515, 8211170.724937779, 8214995.461090431, 8237084.0518939905, 8301608.976428049, 8325546.834221977, 8329060.448176033, 8329093.731634366, 8329432.8025122015, 8330336.141706665, 8331338.3865176905, 8332301.182745698, 8332395.3125, 8332757.747427058, 8332923.852816924, 8335701.858132108, 8337334.874439936, 8337524.285666104, 8337839.834648788, 8338199.886158746, 8340086.02878969, 8345100.904741182, 8358002.910713149, 8358550.211738515, 8360208.245350861, 8360759.375, 8362625.0, 8377450.0, 8435392.255089872, 8435463.715977885, 8441497.624795098, 8441640.45383395, 8444071.2346073, 8444588.599189525, 8445291.668337902, 8445574.794458464, 8447351.888398295, 8463628.047201036, 8470313.079960566, 8486269.05751779, 8490886.429052072, 8557391.842824403, 8563310.238943413, 8588580.784117008, 8590805.300327597, 8591990.304094246, 8594954.340012152, 8618866.668761615, 8618886.231965266, 8624191.689141514, 8625593.439316887, 8626881.764455257, 8649429.223971142, 8650964.538331347, 8653497.23889746, 8702334.351065807, 8732760.790100006, 8755150.106053265, 8757763.728888875, 8759719.464865796, 8760619.811208945, 8763292.18516511, 8763292.1875, 8765500.0, 8767821.64236096, 8767959.026907006, 8795708.370933587, 8797537.5, 8801548.659406807, 8857623.644137666, 8885002.168992119, 8931637.5, 8948024.091380334, 8962002.827475922, 8969404.098924745, 8970245.078404859, 8970995.3125, 8971434.551879466, 8990065.625, 8991261.143546656, 8993468.908154646, 8994562.160221824, 8997750.231978243, 8998873.4375, 9010514.0625, 9019895.802334512, 9027206.25, 9029730.804510111, 9035000.0, 9045142.351123119, 9045557.463825863, 9045923.630364526, 9046205.293393629, 9047670.3125, 9047795.631100819, 9048142.034433374, 9048476.80035415, 9048785.947103929, 9049496.811431251, 9064260.280110447, 9066629.069822852, 9073779.6875, 9079368.75, 9081285.388403768, 9093103.652019558, 9113043.909975333, 9122067.780470518, 9131460.475659497, 9171338.606156098, 9172570.3125, 9179415.625, 9198651.040966596, 9203243.416915242, 9205006.25, 9219766.411712954, 9221012.5, 9225635.534399133, 9225745.3125, 9226133.51447952, 9229023.616649669, 9232575.0, 9243094.852259452, 9247428.312908007, 9289077.112029158, 9314314.726417897, 9356561.940966213, 9360939.0625, 9364352.989204347, 9371840.135847941, 9371862.612425864, 9374944.064178856, 9377998.4375, 9379875.936855169, 9382115.14982822, 9396960.795031626, 9397510.762362551, 9400763.790998338, 9400784.611448744, 9402002.549305363, 9403455.024575658, 9403457.884943636, 9404376.5625, 9405761.46749293, 9405914.0625, 9408407.770410601, 9418831.141779719, 9423271.550803911, 9425867.813507302, 9427284.478111198, 9427374.674525179, 9429082.991766201, 9431965.625, 9432021.875, 9432260.738258172, 9432262.340787336, 9432274.536176618, 9432602.846919704, 9432880.020945568, 9433003.125, 9433318.221633919, 9433404.217507947, 9433414.941164924, 9433717.485008687, 9433965.20944518, 9434406.25, 9434491.66853627, 9434503.813545281, 9434704.6875, 9434822.20995872, 9434828.125, 9434942.1875, 9435025.604198152, 9435509.040992258, 9435520.271206118, 9435630.230430307, 9435657.8125, 9435839.62672939, 9435859.18247504, 9436091.694685213, 9436098.299875695, 9436182.8125, 9436563.029858764, 9436611.011853868, 9436636.519795697, 9436651.971010687, 9436740.558339758, 9436877.38446657, 9437382.813397799, 9437571.875, 9437745.3125, 9437757.8125, 9437759.527289007, 9437777.851352466, 9437967.39339713, 9438185.258683223, 9438515.612502528, 9438753.330132032, 9438861.721077124, 9438911.581195638, 9439590.625, 9439960.62207716, 9440740.1627974, 9441598.622983204, 9442841.666150799, 9444425.292663345, 9444469.119553005, 9446005.14197341, 9446172.844474513, 9447012.5, 9450673.4375, 9454821.875, 9454986.440922506, 9455357.530089725, 9455485.429943906, 9456978.125, 9457161.705962108, 9457779.322191294, 9457830.432026725, 9458060.916366607, 9459089.87700894, 9460407.821022231, 9464142.1875, 9464836.849501545, 9465302.195138903, 9465729.466735832, 9477117.1875, 9478935.739605473, 9479725.312078428, 9495448.4375, 9499450.0, 9499516.037360229, 9500297.589766206, 9501476.185267553, 9501576.5625, 9503620.274694307, 9526613.904179893, 9529504.591730673, 9536462.5, 9537227.686359692, 9539370.678346911, 9541635.74898897, 9543403.125, 9543963.479973812, 9544760.265846293, 9545089.878665935, 9566414.601978991, 9576412.283920541, 9578401.31579933, 9584643.281961033, 9588043.188141929, 9596770.084800536, 9603180.778193986, 9605610.737446869, 9607126.478778133, 9611895.374009397, 9623030.122808771, 9623271.499024883, 9623290.329741282, 9625993.402314162, 9627225.023797566, 9629996.275939662, 9644870.226136304, 9645815.340879427, 9648001.46177157, 9651005.269882688, 9652051.5625, 9655262.5, 9663489.603204867, 9667303.125, 9672744.890929077, 9674814.0625, 9676816.011572262, 9677023.031553837, 9677025.35669642, 9677832.375236703, 9690299.351021303, 9690727.210483301, 9691023.096060744, 9703917.961847419, 9707123.226414874, 9707942.1875, 9719890.625, 9720425.42306955, 9721730.703852907, 9726236.738438789, 9729999.451539755, 9730955.54459186, 9732982.478701603, 9736371.875, 9737835.203441843, 9738841.3534429, 9740083.773235837, 9741128.478713376, 9741784.457737312, 9742604.6875, 9746357.8125, 9747543.75, 9751997.449634438, 9754362.410861436, 9754582.8125, 9760983.70510806, 9764226.251113283, 9765085.236234723, 9766337.86459973, 9766383.748807963, 9769110.686924092, 9769869.84502496, 9770868.75, 9770919.671887003, 9771246.875, 9771643.825479457, 9772663.497492108, 9772757.557257941, 9773170.315507691, 9773367.252160383, 9774263.16861902, 9780753.797575934, 9821015.176074697, 9842382.8125, 9852063.579718066, 9853594.8962652, 9891549.773014551, 9896141.735449243, 9900335.406103969, 9900725.94183369, 9903371.383948818, 9903517.47241881, 9904158.110073226, 9908627.058942651, 9925148.4375, 9926445.0137219, 9939159.375, 9939535.782134289, 9940564.50152192, 9941115.625, 9941650.0, 9942491.67337611, 9942577.70440662, 9944264.227969505, 9947063.212841185, 9953679.904006047, 9958560.9375, 9959659.375, 9959668.436231282, 9965295.3125, 9979183.892351035, 9979589.143452322, 9983087.5, 9991424.049376253, 9993962.972848086, 9994070.3125, 10015067.1875, 10037926.173215073, 10038112.5, 10038951.35580175, 10043943.138549088, 10047739.0625, 10048773.98275896, 10105865.004850542, 10108504.549188336, 10168835.845919825, 10234386.021237548, 10244173.736683426, 10312677.48105694, 10373066.132452358, 10377382.624846734, 10505564.0625, 10521202.80615693, 10529198.051162573, 10578404.975586843, 10626162.821523914, 10644178.550686719, 10656762.5, 10668391.390841668, 10668418.22461328, 10677596.875, 10749830.730311716, 10762576.336158361, 10767298.429648796, 10875609.333053276, 10883840.192975476, 10948213.609016996, 11151757.8125, 11160524.610786699, 11173220.202527693, 11299404.455005264, 11299428.778118704, 11327404.138420654, 11332867.404561844, 11333410.9375, 11333882.085358296, 11352921.041245257, 11353872.58494068, 11354700.312471772, 11355431.770494482, 11356962.5, 11356997.011465302, 11357232.148145527, 11363069.024270117, 11366154.6875, 11368171.233499516, 11371871.192846676, 11400900.0, 11401667.385170722, 11403754.588236444, 11405050.952554448, 11406632.782814851, 11412135.9375, 11416020.3125, 11424591.399408825, 11427416.139031742, 11438531.25, 11438978.138026679, 11440947.77955094, 11441187.331962826, 11448940.57226916, 11448963.207845934, 11455521.601381162, 11457763.37095222, 11458092.1875, 11459920.117868388, 11461113.792870315, 11465185.9375, 11477218.813439531, 11477919.771468744, 11483428.006083459, 11485238.854629679, 11485408.605975455, 11488195.3125, 11488372.636590999, 11489262.5, 11489702.987376181, 11489955.52887039, 11491079.81967093, 11491269.316853989, 11492766.729524102, 11493381.925117178, 11493442.1875, 11493888.08592875, 11494912.09131952, 11495871.43439221, 11496816.759470627, 11497315.91217744, 11500639.483836515, 11503382.8125, 11505615.625, 11506890.210830119, 11508698.954743141, 11510505.86428975, 11511813.079705795, 11513998.138291761, 11514951.203077454, 11517391.415535057, 11518128.125, 11518866.575209213, 11518898.826402118, 11519306.393639904, 11520620.38704898, 11521648.4375, 11521871.875, 11522048.218101023, 11522065.415555563, 11522067.620109076, 11522129.288961964, 11528187.17296505, 11528297.055454506, 11534693.75, 11534714.448949916, 11542151.552611802, 11543082.218131347, 11543309.005644133, 11543379.945199056, 11544802.442480206, 11545403.125, 11545555.560950374, 11545860.9375, 11547273.917260563, 11549045.525223903, 11553188.187430112, 11566474.325492933, 11570828.341470277, 11575864.0625, 11590032.300134156, 11592060.870179474, 11594838.869490411, 11594995.3125, 11595126.662241884, 11596142.14975528, 11601490.807095895, 11603089.0625, 11603600.415995216, 11608178.397806078, 11612622.270998115, 11613879.601216013, 11622389.0625, 11624160.714110233, 11629766.648679698, 11634683.1117809, 11637053.1077956, 11637378.6952132, 11639658.86344426, 11639910.831204617, 11641647.445234919, 11642433.430389158, 11644003.125, 11645330.865115263, 11648181.528639555, 11651168.693455769, 11655637.856510479, 11661164.0625, 11661182.733942738, 11661560.299101988, 11662399.840582158, 11664781.25, 11664840.394941106, 11665268.514318293, 11665538.88459844, 11666531.25, 11666830.473082704, 11666905.311987659, 11667506.191202579, 11670328.157690931, 11672012.5, 11674246.977770058, 11677543.75, 11685315.382099167, 11693904.252602592, 11693913.72523525, 11701722.687697079, 11702433.70238604, 11707979.6875, 11709510.431218628, 11711466.339072473, 11712299.960626382, 11716733.784559725, 11717059.96813094, 11718998.47461315, 11721661.460789476, 11721713.670180703, 11722691.251878558, 11723494.297862254, 11726993.790713605, 11730350.0, 11730352.892464817, 11730870.873112239, 11735086.28858121, 11743232.8125, 11747436.876557792, 11747899.838587401, 11748796.011783522, 11748797.8415802, 11752307.733269854, 11762040.213197922, 11762456.320982656, 11766809.847482182, 11767024.95197246, 11767847.311315097, 11768232.431149941, 11768729.381362313, 11770964.0625, 11771362.5, 11799322.413064573, 11800479.6875, 11802043.551689828, 11803662.83498937, 11803666.97686269, 11804820.41457572, 11805637.5, 11807121.673764924, 11807131.64604544, 11810110.9375, 11810252.516474968, 11813601.011391968, 11817235.285055207, 11818204.711702053, 11818670.581941573, 11819502.456537625, 11819518.532420712, 11819977.192026425, 11823423.137045275, 11826253.584050966, 11826262.64200081, 11836469.770282792, 11837861.227453059, 11838500.0, 11840973.517801823, 11848420.980863173, 11849158.064101465, 11849518.578272313, 11851240.625, 11854638.213618154, 11855264.775113294, 11861203.867439386, 11862718.005050613, 11879123.258855205, 11882607.023091545, 11885013.534939742, 11886403.125, 11891513.913511489, 11894375.785114475, 11895823.4375, 11896949.939463265, 11899438.925213845, 11904428.405943016, 11906252.886308905, 11906279.123036087, 11906279.592363678, 11932480.48582081, 11933103.057154663, 11933129.382322943, 11933607.8125, 11938084.375, 11938124.54868778, 11938337.5, 11939007.8125, 11942148.4375, 11942570.3125, 11943890.625, 11949318.75, 11950425.55961304, 11954123.955127371, 11967712.51602519, 11984147.400599526, 11998102.145341592, 11998440.688169323, 12003535.09044335, 12003643.562584821, 12003880.92697216, 12006603.295409188, 12006666.224627413, 12006818.280884814, 12007656.063976634, 12021760.904552123, 12025167.1875, 12026949.351648882, 12034604.825674282, 12036595.022035923, 12036644.718651596, 12037748.629993232, 12046750.480471171, 12052907.079733664, 12065759.22314647, 12066550.0, 12067487.860439241, 12068535.095491441, 12072198.31479832, 12072784.063159177, 12075505.769797347, 12080189.0625, 12083563.44953464, 12085642.1875, 12093829.193722008, 12095408.424948497, 12096765.275814628, 12097687.283872547, 12099518.427992776, 12104641.462846002, 12106722.425311297, 12107981.419838626, 12109543.675004374, 12109555.5617175, 12111009.40173332, 12118956.25, 12130685.9375, 12156460.963113505, 12158903.861347273, 12160021.875, 12161343.516034434, 12161435.442353878, 12162293.576249557, 12171394.532892589, 12182315.920075202, 12189519.792895334, 12190764.1994611, 12228982.8125, 12250967.39681442, 12257127.38689808, 12260809.874362111, 12277119.85655168, 12283352.550303686, 12284861.559482442, 12294715.167966884, 12304309.03100434, 12304940.136731103, 12307857.8125, 12307876.21857812, 12312274.954218056, 12312325.0, 12319710.193869041, 12333271.139402378, 12347278.705457788, 12350305.702439817, 12350903.774050955, 12352714.57503582, 12363439.0625, 12363475.027313102, 12382578.125, 12385620.09127635, 12395413.514889034, 12396338.908420142, 12404720.189026555, 12414685.262414258, 12429753.50300578, 12430477.592767531, 12491428.997429512, 12491487.455662034, 12505950.0, 12521524.194384303, 12522096.733955272, 12535976.748088233, 12545199.56781789, 12546245.514236271, 12585963.532921545, 12587055.081408978, 12587098.213551808, 12601330.047961889, 12623849.354517985, 12629663.63298652, 12630808.627519049, 12631189.776683358, 12640764.732623996, 12640784.12392131, 12642096.475688936, 12644814.0625, 12648048.011604825, 12650172.536429618, 12650685.788957292, 12650790.625, 12659515.579054026, 12687373.23501867, 12691241.793954123, 12691435.9375, 12694668.75, 12704594.864356427, 12725418.96091594, 12729585.923266191, 12755124.684385259, 12765718.094615197, 12766621.71747, 12767161.154606812, 12767207.4480891, 12778809.546092963, 12784032.8125, 12785545.3125, 12788303.125, 12788700.0, 12788746.79849205, 12789178.669848654, 12789427.616552504, 12800002.374935947, 12806134.212748306, 12826503.49843812, 12834753.649074761, 12849625.910309618, 12853181.06132268, 12856142.974524688, 12859390.97912963, 12861295.3125, 12861796.738231504, 12862332.296256898, 12862547.330079835, 12862806.167456886, 12863685.9381652, 12864459.375, 12865390.625, 12871498.191598183, 12872913.646558296, 12874745.59387956, 12882815.625, 12894935.9375, 12896462.768712457, 12897845.3125, 12920970.234010998, 12922365.25815487, 12923056.320378527, 12924702.32949444, 12927113.808847269, 12938816.046170648, 12955906.25, 12960581.25, 12961335.473904584, 12971468.75, 12982040.969907956, 12984873.4375, 12989683.307689022, 12992551.296673022, 12995479.6875, 13004751.5625, 13006326.5625, 13006348.4375, 13016577.63930378, 13035602.952919273, 13040328.60738103, 13042052.982904142, 13046184.146535974, 13046518.75, 13053731.042401625, 13071697.86593385, 13087968.440281667, 13101690.57579928, 13102119.088868024, 13107442.565876858, 13117316.501595251, 13120861.104912566, 13126672.090813775, 13127017.998469863, 13127883.434762586, 13128695.3125, 13130914.41729843, 13132032.776201544, 13138273.227820558, 13142433.123649413, 13147455.761792997, 13154959.424166335, 13171745.332598332, 13180323.354040999, 13214886.312817538, 13215305.013614299, 13231551.5625, 13249223.076949947, 13252125.051517086, 13263656.25, 13271519.978331707, 13285523.1385595, 13290871.863095077, 13295278.125, 13295300.193188995, 13312300.607757116, 13313914.01212709, 13357484.126843406, 13366547.027387178, 13368765.691720141, 13385956.25, 13386588.673917877, 13398531.333380474, 13406261.603192013, 13415476.465965934, 13482379.010671396, 13490359.342487063, 13522149.751474824, 13524280.220960297, 13524679.633438936, 13555312.792130562, 13573266.601688737, 13614512.5, 13630402.222515922, 13633070.3125, 13663222.470951432, 13678840.625, 13678887.108327253, 13721999.292205464, 13725011.222197175, 13725019.16909765, 13725500.0, 13725791.521667438, 13725820.872642437, 13746900.0, 13747684.375, 13792454.6875, 13794546.984439723, 13808651.80619715, 13822718.75, 13970953.879051354, 13984749.610296266, 13985732.8125, 14049518.555861633, 14074421.473048678, 14197527.130071469, 14239943.910490653, 14382580.759306889, 15372881.202261753, 15377288.481698131, 15392935.885784997, 15418070.241272867, 15481532.8125, 15528969.927964639, 15562996.770062903, 15594319.589288566, 15695528.125, 15862687.369659513, 15908140.806831326, 15934082.8125, 15935221.875, 16029304.6875, 16029330.713078398, 16034444.10512927, 16049114.775341185, 16069709.727653597, 16074374.445150739, 16089023.4375, 16098877.835589826, 16111560.213344913, 16131173.669604735, 16141026.911057696, 16156287.894219257, 16171649.791964475], [10.46005308030074, 14.904816831114386, 21.509579953826076, 5.315571685515921, 14.641153547339282, 14.551713151450613, 7.477746289301568, 20.004992277642685, 178.2872063856559, 59.93966451642283, 51.362576770807834, 13.1035239783946, 34.195371794219966, 6.353521381156771, 11.876391786882824, 49.29720514225406, 132.36260141018587, 53.10963427763747, 18.161256474081533, 7.655510267202092, 11.465416502899068, 46.9817648158792, 21.43463882913759, 40.35769116738798, 18.76656052722761, 40.98496199447338, 19.002676457902197, 28.500643229868082, 37.34412429749422, 36.900183409281084, 6.132065607959456, 35.5600762244857, 9.07236688080641, 65.35378508178178, 13.5901635727614, 47.518673981573706, 162.00679272838605, 11.112131036531967, 24.882819981192128, 8.818512170568809, 22.782783631012762, 33.151835487243325, 8.252462851111957, 22.97597363647257, 45.00266508243338, 15.39863943033076, 7.941796522170276, 19.47197065545811, 39.32959418108189, 20.28250527145475, 10.23476704601325, 61.29166591925012, 6.520803373945005, 6.557154839231342, 156.23966800287673, 35.28726230336729, 67.14434771748074, 9.389215692926406, 13.767004199774586, 10.05756550785566, 29.141779714716762, 43.686308098672754, 6.041208707818841, 10.40283654052718, 39.58016733270071, 7.209839204682123, 7.575339423011128, 9.875177979916275, 21.375157601033006, 26.410988076161814, 10.6285540427579, 96.88147833538343, 8.9662689655524, 15.183393585622014, 23.48732289220043, 88.3852572723313, 115.47603765312692, 16.196714113826587, 107.40044113737153, 80.10397482998351, 32.35797879642434, 8.967979769757086, 48.2424149815306, 5.634498105925283, 5.531920066023435, 25.917548565327237, 5.466336949612975, 100.60185238949083, 121.646224449383, 44.386850200776806, 28.479024807073255, 164.59426231862466, 45.88221007666563, 42.30776035408241, 23.155262788130653, 64.31043972340991, 32.77808818872232, 9.526361378626401, 21.826275585055114, 51.28072073653459, 61.93652129939126, 25.12061613846303, 48.4234809169319, 96.70928611494756, 61.733659847068594, 45.302781941492476, 11.415913640130704, 51.03460121989626, 61.56974874699588, 13.431149996742056, 23.255589646389794, 11.510755110375927, 60.410657109469916, 58.88748526514357, 5.141110712318085, 74.57940042074148, 22.05466600158739, 51.29357868028993, 9.473697022992809, 7.481815424423114, 14.03379984206384, 88.5330791497979, 41.56249942595854, 11.665106763904873, 14.388513035873158, 21.800541174700545, 100.96898633824384, 11.02977753813278, 33.29011745919081, 20.394383983071645, 63.55380257386972, 43.919357458252804, 23.01024424148823, 52.93824705598195, 75.43661023465496, 9.888737247731447, 10.475180869403623, 12.89789137394039, 86.39553385485206, 96.41423892357255, 46.41316530770811, 37.556570428969856, 49.45516408339277, 6.617917066333705, 54.147849956046876, 8.084325092096178, 35.802227558745834, 62.76021087535632, 56.88338260245663, 130.6711572084285, 40.56750600212464, 25.49335405251625, 6.576829047498817, 52.736603576252236, 61.237742416801304, 48.20686758601432, 6.074090595487571, 15.940573885565586, 7.151321960700783, 16.2914429290514, 54.63575849855518, 44.69204191866788, 32.98275736333001, 25.995671096857322, 14.31349780202911, 25.033433646329264, 12.407688834811019, 6.197103995513198, 11.698890112459546, 11.242860867452213, 61.46559166622376, 66.70496444239848, 29.360754381563762, 14.196409422614167, 173.43949869419188, 205.21444123748222, 17.427440157357893, 124.60847524303543, 42.2119156838969, 9.487726679291303, 32.12754485945547, 51.243098127498214, 52.13001059678446, 14.367510538450862, 21.861735304659298, 13.366603473479222, 54.168829738192784, 72.72168227726169, 70.9424489442474, 89.08706179555452, 7.373313596795591, 135.87841700917642, 13.563428602116813, 52.294007630488146, 6.03574452193149, 32.46112106928729, 80.24642603544508, 41.77935585724417, 62.44800718780023, 85.9149071529233, 95.6962602770132, 25.3958146126336, 9.296729245050015, 7.271880166555046, 89.8361150902715, 17.639128628348196, 127.51661966346144, 61.41548376518578, 91.77511047191413, 30.09324106604355, 17.436145515896865, 30.82023086744954, 5.562024721678407, 70.33998509823505, 14.693376565335392, 23.424056867223538, 27.749510271548267, 52.812905721240064, 30.3098036361471, 99.51020235291591, 5.540517540550961, 31.772223427695096, 45.8965926027949, 7.374882944576174, 6.086800553258134, 15.785608154027267, 8.749506438162307, 7.301604934364692, 77.32987370266302, 23.555715244751187, 21.8786219800095, 8.387161530518998, 38.02464967434318, 80.2738366049172, 11.759039832530886, 23.086243794466412, 5.771366899993589, 39.85785672538417, 39.12722844289006, 92.91403045781837, 156.85731100439705, 94.66225747566423, 7.429430069475531, 67.44153409041668, 74.59239359829178, 114.4623992148571, 10.051727391426644, 28.113319219764335, 19.331010790221494, 16.11048730103169, 51.794892973848604, 41.95453607240488, 54.95996210986396, 7.5252960554845805, 35.53186933286451, 99.96477038609784, 9.610001022015487, 51.56846671376123, 28.982923462512446, 86.05114115112606, 18.51965362359951, 25.195128065048458, 53.13692907005788, 62.34127565818207, 59.715944529763384, 9.25171666538587, 51.6489905775533, 46.17152996801633, 18.453163735074376, 25.287592676196066, 43.787470420982025, 83.66012110656628, 46.6402921641026, 19.34364728420873, 12.090490520501456, 6.240278905413386, 19.012242862990725, 6.416142371763144, 30.446169598231187, 14.139277350873845, 27.978337069651904, 71.59450683893728, 66.36487197355532, 116.90839376226157, 59.17809813705439, 29.775401455902266, 70.90357127638957, 40.04808196076903, 7.611867593801878, 98.40674967720241, 26.92596270771244, 14.497536487310022, 16.02166294164095, 144.06385216261046, 48.62399427477424, 73.47949365020442, 8.594823383759417, 74.05095983615143, 13.546102842927178, 46.63661713175448, 69.56341089096935, 21.91789367967564, 11.735111628696554, 22.142625269171518, 32.002270048779195, 42.204814344109835, 22.899857877829845, 19.031111692547118, 59.49666932087761, 44.627803443600456, 67.05204489559316, 17.508628846794508, 30.39458265216645, 43.345792690808075, 21.022931173373195, 29.070446501135432, 29.99797570216571, 52.102735893667, 6.377865142572682, 5.266000795281109, 59.62495941644705, 27.277400204586463, 206.84018002133402, 45.02901706739867, 68.19244353274112, 42.37770458631099, 78.30816841088722, 185.95939386039947, 70.62829096274824, 115.16122749467726, 40.218162398003315, 32.60279394256348, 16.540756214495634, 55.50737314910614, 39.390508513644036, 27.728811880335755, 95.15657115800794, 12.956759891372322, 71.7314885580177, 38.210511165429956, 5.506273196787376, 50.67208622977007, 88.07391667117014, 15.193828212766471, 112.62899174153308, 37.5015304182949, 22.018191016381163, 68.47414725799158, 28.32038456990839, 7.016172240331582, 28.445266863418382, 47.20740459649353, 21.185522844303225, 6.967326293996054, 42.11691045642713, 26.45409229157432, 155.27544494740954, 5.177086943933869, 87.11463274885384, 98.15602438969994, 41.45646777609386, 40.82943453484009, 6.134975346618311, 110.64372977868027, 62.2339636612124, 21.353354503218227, 17.690927850082105, 41.1880027080968, 18.59438752306608, 5.251245528259005, 80.01210067352088, 36.00276588777652, 18.81077143953084, 55.118472665717626, 30.039363902200748, 74.22049446597399, 46.90936199130296, 9.908559986058348, 17.11277845915571, 109.16717911069404, 37.37516937285763, 35.51583663351791, 62.192104250686526, 6.722930732293959, 41.09043166829201, 107.44944307224371, 76.4716792331794, 13.31265237845868, 26.0619587964073, 14.60357168296921, 15.270581374537088, 5.547779206414565, 74.15917302598947, 17.34287689319821, 5.524041923611538, 91.61268523532743, 82.23288021946027, 28.66258016231189, 25.083309343602085, 18.305894510912008, 11.865352898185584, 29.852207572898763, 22.33450854734035, 13.257051405218004, 49.726473677150686, 54.601991080018976, 71.2597066597719, 84.01063941639904, 35.49505355652587, 93.79843519859801, 15.68269031646279, 44.26204534874868, 22.67434708861132, 112.24608451983018, 13.436237123912555, 6.353399137402942, 21.017091219679685, 114.37584594732776, 41.89707069831421, 6.942584154398579, 16.50850163933523, 18.57443256995602, 68.77056582430806, 98.42990465030982, 52.323462314432355, 6.284595118619792, 135.61081354767194, 53.503255212290355, 32.55315815662022, 103.9825048778239, 14.441371419713642, 149.4500938009838, 60.55300531501252, 159.73576478759458, 20.462234155829286, 85.40999915027011, 44.830266160964996, 69.8709751443631, 12.600695424099838, 49.70777220682642, 31.134932741801837, 44.67853744047713, 12.913321184739036, 44.410528568998515, 80.89785451680575, 124.53323675450727, 41.134293649110944, 58.775219974976, 101.97560395560843, 25.196936930288974, 52.84964313625933, 53.20986339075504, 23.607944156480748, 40.20613385199743, 7.09053907774558, 9.161611748811067, 100.40686317953406, 76.8505708126323, 10.070680414025194, 38.680465525076265, 26.07296711513674, 95.11151012711966, 8.357730495497414, 36.605707419413264, 63.8598890055455, 12.377090374751722, 66.75748662706995, 14.968322651239426, 27.804981391383635, 22.10676060245977, 8.993932220826862, 36.9085924589413, 5.6377866380751565, 113.43088218598271, 84.45790557888111, 22.394274274559887, 58.54920632364563, 22.71435657008975, 44.8913671637996, 13.911799463428904, 13.728281485434376, 19.150535145545902, 51.377452720530435, 36.68975611028691, 30.645004890194983, 9.455885651841506, 15.233478004502677, 55.631317699723354, 49.300862500962175, 161.59919467714377, 24.980748227581522, 9.726103225243985, 18.370513047381255, 37.172800474716496, 11.421991218353877, 38.323132887689454, 41.760676636721925, 9.323056115934662, 10.310466763641115, 78.72991463916759, 24.914633917131162, 8.814065689379412, 78.53581305810236, 14.06509920815876, 41.93971447769981, 69.69164099130191, 11.477325140736363, 26.75752412790964, 10.215085828001309, 69.89188717479753, 113.18339962873547, 130.03007456240235, 27.260158710473284, 14.253568695148076, 94.48413618988229, 105.55620901609637, 29.854065329518868, 102.95286200566437, 31.10033772785273, 6.849499239611752, 16.09435204052012, 34.263079763037496, 152.05090160172193, 32.91372536810881, 10.795608878117786, 21.11439496595769, 26.57605755229211, 234.4321639864748, 18.77366329613951, 77.8789585018097, 40.14073192674408, 12.15367125819885, 45.365951030545986, 39.42761410472773, 56.64166862771388, 7.707224747318835, 45.10793019880295, 24.015195898340927, 54.35573414782981, 7.57506984591888, 8.84424750563413, 7.080166757187136, 76.26817418251031, 6.238563157711726, 9.139766552876395, 24.02479033910739, 54.97474384929673, 8.750830569022959, 5.02608325252503, 91.51394733406534, 6.7700506558429225, 56.87981644273842, 29.39628482227317, 28.831415817558295, 71.99769897983538, 24.221276936329392, 6.435490489730579, 10.128202661889697, 12.87691591399845, 42.77654693933018, 119.9671510729133, 30.279461042153734, 41.45253137717987, 20.94816474761589, 12.201378308452913, 11.18621682671237, 7.266155582998123, 5.445175995065317, 14.214353619881635, 8.599226627305159, 39.8440815802911, 42.757769442579985, 24.652540503521674, 48.54345280295502, 22.575182363354173, 90.62945528068869, 67.0776418832133, 39.01717021678313, 41.456834447734494, 83.48234554651404, 14.647882435466503, 33.68568462216086, 147.80540979695195, 6.328797527910687, 47.703484041433114, 144.19888872822813, 8.094708284892777, 33.29828771044275, 57.35404278575308, 33.92769941592709, 13.471009960733731, 74.09015515088883, 5.694252196734011, 11.69505950570191, 24.963533392118325, 69.2989177756227, 107.10465511599463, 8.18295246096663, 26.920217970702648, 80.408428396774, 13.038334115335761, 61.124891169077564, 64.83192814811643, 9.351951560310564, 18.65260142067794, 12.887190726114909, 7.526394040069647, 117.8542058627798, 6.7049697792421235, 5.305480799182547, 15.511858167455491, 158.003421471843, 15.14468356362237, 100.3086861773746, 15.256678399609854, 6.160474786050645, 9.553610452182832, 48.341340959699366, 54.47442227717069, 66.6248936672985, 10.016615835763426, 18.80126186023898, 51.17998097565463, 14.301679867052636, 77.52751029771365, 8.921892657072155, 54.89846921988301, 17.33076715646697, 102.10316010293691, 38.450140532371854, 128.13530916122204, 59.607154563478474, 8.44095068429451, 24.566795988847897, 10.346857269145369, 10.982188901561608, 84.31591920415146, 27.25081613684647, 5.092496767658493, 45.78943359716188, 45.102119009829224, 18.9284950636643, 51.104402998203305, 111.79485160976334, 49.53916811042245, 81.03624590685057, 82.52736583402798, 69.3828527196381, 38.15212096690889, 78.56843966412396, 57.75017315547217, 66.56362324616298, 19.915158544044868, 6.7761837034514665, 12.46642960681854, 26.503109069968744, 23.4609830696894, 200.6357994925697, 146.34723781013432, 26.060271222565984, 82.80220866999161, 27.889096638261737, 69.91641671835687, 187.24765842954773, 14.764638400887188, 23.873318930481776, 10.311096077659993, 62.030587279356716, 11.972088670592324, 12.730767367373948, 14.744926948739588, 75.13660392289218, 16.51979112769956, 45.92522289750906, 26.314558628522253, 103.6013029108133, 110.19349901986118, 128.84037061445736, 15.108716975852179, 107.45754482563311, 13.343119011495226, 23.62001218040819, 35.52648172839175, 65.78981482359202, 54.23635156563837, 19.510966824017967, 54.86017367432424, 15.63093797440105, 285.41646802802813, 45.69000009228574, 15.331499618514034, 52.43298035522011, 34.27608752401803, 15.040043169067104, 6.919990493655217, 71.6797750061305, 24.632552222249586, 31.53428674693425, 69.19490509983791, 94.91574108996396, 8.449035442929048, 13.272672253935308, 180.03725244602123, 38.36112770745725, 15.80421830495635, 12.718044972855013, 57.66727139151401, 124.32413496646772, 5.890872749530303, 141.78366850821695, 49.21254345002892, 24.862206126987317, 127.11754308219128, 108.01070954329697, 191.72434928058928, 72.02571590310596, 29.860304521020527, 5.574188017950052, 33.00205110082356, 17.737171289568295, 34.43730538522372, 5.111448490476592, 87.3112673130835, 81.14294925463486, 25.06552323976976, 30.416920697069553, 24.81200644198517, 110.79433240950424, 5.525648997796797, 82.54634276357329, 34.207482291239515, 7.2244942474555, 50.393069818758974, 289.9863771281537, 9.294192209314993, 59.506221368054995, 148.01703662683747, 14.429850073122667, 7.1711069596050905, 129.74907338994228, 20.36961735410386, 8.920454949366103, 5.552508556854841, 6.06511613050196, 29.70379490382412, 6.641138105415583, 59.772044043957976, 32.0481623617099, 60.902347472461955, 95.8247708695723, 66.64354265054752, 22.42525927588848, 83.02261752987962, 76.47169997903089, 18.248775239295398, 26.681639978982467, 84.20116422305013, 149.20322021001482, 45.67598245106359, 16.11656586360676, 109.41111175694266, 150.36828264807724, 6.419062870226426, 32.98809616617845, 40.07854015769477, 133.3801368369283, 13.186270935428443, 25.758832274902964, 262.3294823261272, 56.257633196179754, 6.2890594725513616, 58.52036261183278, 95.66214036585107, 118.29535592895849, 56.59520466918407, 9.26753851237654, 17.523904749402295, 21.409243826750533, 12.689661624799005, 111.58000988517797, 67.07812222020456, 5.986382592171176, 120.56743583247504, 30.441937892191945, 149.67243871328762, 27.94921410683981, 39.75178348026374, 6.736505556424783, 197.27161225840402, 14.063270527564626, 6.07955731649008, 24.509736987045585, 34.843079116447335, 105.9847749437468, 36.87142771067538, 9.002835603148894, 77.89318863267114, 91.03190717933538, 55.24562061211194, 70.16081829807507, 39.834999718723076, 35.282617266385046, 22.932858627925537, 13.083635004557983, 36.02916338132999, 24.3613847957462, 20.7625135836459, 111.11582473938283, 74.89439280824058, 157.00914179930862, 79.92202177784577, 16.442159872981588, 73.87337579656358, 62.14438589694732, 59.148395646215235, 68.094066942419, 38.53874835114903, 74.70283607632217, 47.317640502741654, 5.778896196151821, 172.2530912279364, 69.49620086919462, 26.001595138787152, 24.709647298132136, 52.415221573384656, 69.41312245627446, 76.84416934341235, 23.627392125530026, 125.59241210530469, 130.8182542375304, 5.973048415159237, 35.23535088861591, 161.4597549767803, 40.03074657455956, 142.01891556901927, 26.863609222982504, 44.89134253814364, 58.96898944948762, 113.93732092138178, 9.160184686044985, 75.63429296358292, 149.2098021182579, 85.2168357443131, 256.4584291439799, 10.999713587102859, 22.915482147833956, 7.725605761112345, 79.79078745708637, 41.31199146266421, 5.04241315796255, 52.74131554358093, 32.29531434318677, 10.690523101746015, 57.08718483670773, 48.972701567198854, 11.443724926437923, 5.9941398408934905, 218.0114386269908, 10.760779063747075, 6.794315663464759, 15.202063307290153, 108.53077834817229, 120.22234106455603, 6.780267788600447, 26.185225598766877, 26.75143978064268, 7.9583949057945675, 17.24768056533221, 33.954093994392, 47.04282213920154, 27.498754487763403, 41.33288818306591, 20.6777671441621, 59.0258496188467, 43.87522345295311, 61.545955207928706, 42.87074019796485, 43.07556484745181, 10.566121665295805, 26.321645473810776, 6.303850598237803, 9.712834439690447, 80.67118243646796, 52.95424603686472, 6.148352706654838, 53.67906388039738, 60.92083116930623, 75.68472294765404, 113.31443379749983, 155.35045909199408, 36.821991516969675, 160.8859167163571, 56.264180861206974, 78.00386262070087, 30.90678716430718, 12.780285019615246, 16.233017842299414, 15.590668832209131, 29.917501354037608, 37.323004688903474, 25.80216880840069, 30.995986736378207, 39.21098722866333, 30.919047224710905, 22.727814482073246, 10.967551429868204, 109.9834275188002, 83.32517340778499, 55.833565750028384, 15.68302399677503, 43.73613658005283, 52.28963701482757, 48.21845633612445, 5.27478532295761, 16.49313810919062, 37.4853504266769, 86.83622237589096, 26.88375651187652, 43.64692058661242, 10.703590272493596, 19.51486118811333, 11.510975767421295, 15.340383725070527, 35.44902153180603, 57.20358509925974])
caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);
([6545529.213785079, 6601443.740339845, 6662226.602305766, 6732844.522830824, 6817765.543901257, 6907381.143830056, 7325765.2641742155, 7759829.056142265, 7802295.032793387, 7834572.64073461, 7868227.138216119, 7933558.791688724, 8032593.463836384, 8073361.591505273, 8092861.36892458, 8105679.057481564, 8106128.95183724, 8106993.75, 8108664.7105028825, 8112257.954749868, 8113650.262079768, 8113853.125, 8117096.369119483, 8118103.11295699, 8192141.560404632, 8194721.875, 8206549.852044636, 8208863.197536281, 8209642.047681872, 8211166.8731837515, 8211170.724937779, 8214995.461090431, 8237084.0518939905, 8301608.976428049, 8325546.834221977, 8329060.448176033, 8329093.731634366, 8329432.8025122015, 8330336.141706665, 8331338.3865176905, 8332301.182745698, 8332395.3125, 8332757.747427058, 8332923.852816924, 8335701.858132108, 8337334.874439936, 8337524.285666104, 8337839.834648788, 8338199.886158746, 8340086.02878969, 8345100.904741182, 8358002.910713149, 8358550.211738515, 8360208.245350861, 8360759.375, 8362625.0, 8377450.0, 8435392.255089872, 8435463.715977885, 8441497.624795098, 8441640.45383395, 8444071.2346073, 8444588.599189525, 8445291.668337902, 8445574.794458464, 8447351.888398295, 8463628.047201036, 8470313.079960566, 8486269.05751779, 8490886.429052072, 8557391.842824403, 8563310.238943413, 8588580.784117008, 8590805.300327597, 8591990.304094246, 8594954.340012152, 8618866.668761615, 8618886.231965266, 8624191.689141514, 8625593.439316887, 8626881.764455257, 8649429.223971142, 8650964.538331347, 8653497.23889746, 8702334.351065807, 8732760.790100006, 8755150.106053265, 8757763.728888875, 8759719.464865796, 8760619.811208945, 8763292.18516511, 8763292.1875, 8765500.0, 8767821.64236096, 8767959.026907006, 8795708.370933587, 8797537.5, 8801548.659406807, 8857623.644137666, 8885002.168992119, 8931637.5, 8948024.091380334, 8962002.827475922, 8969404.098924745, 8970245.078404859, 8970995.3125, 8971434.551879466, 8990065.625, 8991261.143546656, 8993468.908154646, 8994562.160221824, 8997750.231978243, 8998873.4375, 9010514.0625, 9019895.802334512, 9027206.25, 9029730.804510111, 9035000.0, 9045142.351123119, 9045557.463825863, 9045923.630364526, 9046205.293393629, 9047670.3125, 9047795.631100819, 9048142.034433374, 9048476.80035415, 9048785.947103929, 9049496.811431251, 9064260.280110447, 9066629.069822852, 9073779.6875, 9079368.75, 9081285.388403768, 9093103.652019558, 9113043.909975333, 9122067.780470518, 9131460.475659497, 9171338.606156098, 9172570.3125, 9179415.625, 9198651.040966596, 9203243.416915242, 9205006.25, 9219766.411712954, 9221012.5, 9225635.534399133, 9225745.3125, 9226133.51447952, 9229023.616649669, 9232575.0, 9243094.852259452, 9247428.312908007, 9289077.112029158, 9314314.726417897, 9356561.940966213, 9360939.0625, 9364352.989204347, 9371840.135847941, 9371862.612425864, 9374944.064178856, 9377998.4375, 9379875.936855169, 9382115.14982822, 9396960.795031626, 9397510.762362551, 9400763.790998338, 9400784.611448744, 9402002.549305363, 9403455.024575658, 9403457.884943636, 9404376.5625, 9405761.46749293, 9405914.0625, 9408407.770410601, 9418831.141779719, 9423271.550803911, 9425867.813507302, 9427284.478111198, 9427374.674525179, 9429082.991766201, 9431965.625, 9432021.875, 9432260.738258172, 9432262.340787336, 9432274.536176618, 9432602.846919704, 9432880.020945568, 9433003.125, 9433318.221633919, 9433404.217507947, 9433414.941164924, 9433717.485008687, 9433965.20944518, 9434406.25, 9434491.66853627, 9434503.813545281, 9434704.6875, 9434822.20995872, 9434828.125, 9434942.1875, 9435025.604198152, 9435509.040992258, 9435520.271206118, 9435630.230430307, 9435657.8125, 9435839.62672939, 9435859.18247504, 9436091.694685213, 9436098.299875695, 9436182.8125, 9436563.029858764, 9436611.011853868, 9436636.519795697, 9436651.971010687, 9436740.558339758, 9436877.38446657, 9437382.813397799, 9437571.875, 9437745.3125, 9437757.8125, 9437759.527289007, 9437777.851352466, 9437967.39339713, 9438185.258683223, 9438515.612502528, 9438753.330132032, 9438861.721077124, 9438911.581195638, 9439590.625, 9439960.62207716, 9440740.1627974, 9441598.622983204, 9442841.666150799, 9444425.292663345, 9444469.119553005, 9446005.14197341, 9446172.844474513, 9447012.5, 9450673.4375, 9454821.875, 9454986.440922506, 9455357.530089725, 9455485.429943906, 9456978.125, 9457161.705962108, 9457779.322191294, 9457830.432026725, 9458060.916366607, 9459089.87700894, 9460407.821022231, 9464142.1875, 9464836.849501545, 9465302.195138903, 9465729.466735832, 9477117.1875, 9478935.739605473, 9479725.312078428, 9495448.4375, 9499450.0, 9499516.037360229, 9500297.589766206, 9501476.185267553, 9501576.5625, 9503620.274694307, 9526613.904179893, 9529504.591730673, 9536462.5, 9537227.686359692, 9539370.678346911, 9541635.74898897, 9543403.125, 9543963.479973812, 9544760.265846293, 9545089.878665935, 9566414.601978991, 9576412.283920541, 9578401.31579933, 9584643.281961033, 9588043.188141929, 9596770.084800536, 9603180.778193986, 9605610.737446869, 9607126.478778133, 9611895.374009397, 9623030.122808771, 9623271.499024883, 9623290.329741282, 9625993.402314162, 9627225.023797566, 9629996.275939662, 9644870.226136304, 9645815.340879427, 9648001.46177157, 9651005.269882688, 9652051.5625, 9655262.5, 9663489.603204867, 9667303.125, 9672744.890929077, 9674814.0625, 9676816.011572262, 9677023.031553837, 9677025.35669642, 9677832.375236703, 9690299.351021303, 9690727.210483301, 9691023.096060744, 9703917.961847419, 9707123.226414874, 9707942.1875, 9719890.625, 9720425.42306955, 9721730.703852907, 9726236.738438789, 9729999.451539755, 9730955.54459186, 9732982.478701603, 9736371.875, 9737835.203441843, 9738841.3534429, 9740083.773235837, 9741128.478713376, 9741784.457737312, 9742604.6875, 9746357.8125, 9747543.75, 9751997.449634438, 9754362.410861436, 9754582.8125, 9760983.70510806, 9764226.251113283, 9765085.236234723, 9766337.86459973, 9766383.748807963, 9769110.686924092, 9769869.84502496, 9770868.75, 9770919.671887003, 9771246.875, 9771643.825479457, 9772663.497492108, 9772757.557257941, 9773170.315507691, 9773367.252160383, 9774263.16861902, 9780753.797575934, 9821015.176074697, 9842382.8125, 9852063.579718066, 9853594.8962652, 9891549.773014551, 9896141.735449243, 9900335.406103969, 9900725.94183369, 9903371.383948818, 9903517.47241881, 9904158.110073226, 9908627.058942651, 9925148.4375, 9926445.0137219, 9939159.375, 9939535.782134289, 9940564.50152192, 9941115.625, 9941650.0, 9942491.67337611, 9942577.70440662, 9944264.227969505, 9947063.212841185, 9953679.904006047, 9958560.9375, 9959659.375, 9959668.436231282, 9965295.3125, 9979183.892351035, 9979589.143452322, 9983087.5, 9991424.049376253, 9993962.972848086, 9994070.3125, 10015067.1875, 10037926.173215073, 10038112.5, 10038951.35580175, 10043943.138549088, 10047739.0625, 10048773.98275896, 10105865.004850542, 10108504.549188336, 10168835.845919825, 10234386.021237548, 10244173.736683426, 10312677.48105694, 10373066.132452358, 10377382.624846734, 10505564.0625, 10521202.80615693, 10529198.051162573, 10578404.975586843, 10626162.821523914, 10644178.550686719, 10656762.5, 10668391.390841668, 10668418.22461328, 10677596.875, 10749830.730311716, 10762576.336158361, 10767298.429648796, 10875609.333053276, 10883840.192975476, 10948213.609016996, 11151757.8125, 11160524.610786699, 11173220.202527693, 11299404.455005264, 11299428.778118704, 11327404.138420654, 11332867.404561844, 11333410.9375, 11333882.085358296, 11352921.041245257, 11353872.58494068, 11354700.312471772, 11355431.770494482, 11356962.5, 11356997.011465302, 11357232.148145527, 11363069.024270117, 11366154.6875, 11368171.233499516, 11371871.192846676, 11400900.0, 11401667.385170722, 11403754.588236444, 11405050.952554448, 11406632.782814851, 11412135.9375, 11416020.3125, 11424591.399408825, 11427416.139031742, 11438531.25, 11438978.138026679, 11440947.77955094, 11441187.331962826, 11448940.57226916, 11448963.207845934, 11455521.601381162, 11457763.37095222, 11458092.1875, 11459920.117868388, 11461113.792870315, 11465185.9375, 11477218.813439531, 11477919.771468744, 11483428.006083459, 11485238.854629679, 11485408.605975455, 11488195.3125, 11488372.636590999, 11489262.5, 11489702.987376181, 11489955.52887039, 11491079.81967093, 11491269.316853989, 11492766.729524102, 11493381.925117178, 11493442.1875, 11493888.08592875, 11494912.09131952, 11495871.43439221, 11496816.759470627, 11497315.91217744, 11500639.483836515, 11503382.8125, 11505615.625, 11506890.210830119, 11508698.954743141, 11510505.86428975, 11511813.079705795, 11513998.138291761, 11514951.203077454, 11517391.415535057, 11518128.125, 11518866.575209213, 11518898.826402118, 11519306.393639904, 11520620.38704898, 11521648.4375, 11521871.875, 11522048.218101023, 11522065.415555563, 11522067.620109076, 11522129.288961964, 11528187.17296505, 11528297.055454506, 11534693.75, 11534714.448949916, 11542151.552611802, 11543082.218131347, 11543309.005644133, 11543379.945199056, 11544802.442480206, 11545403.125, 11545555.560950374, 11545860.9375, 11547273.917260563, 11549045.525223903, 11553188.187430112, 11566474.325492933, 11570828.341470277, 11575864.0625, 11590032.300134156, 11592060.870179474, 11594838.869490411, 11594995.3125, 11595126.662241884, 11596142.14975528, 11601490.807095895, 11603089.0625, 11603600.415995216, 11608178.397806078, 11612622.270998115, 11613879.601216013, 11622389.0625, 11624160.714110233, 11629766.648679698, 11634683.1117809, 11637053.1077956, 11637378.6952132, 11639658.86344426, 11639910.831204617, 11641647.445234919, 11642433.430389158, 11644003.125, 11645330.865115263, 11648181.528639555, 11651168.693455769, 11655637.856510479, 11661164.0625, 11661182.733942738, 11661560.299101988, 11662399.840582158, 11664781.25, 11664840.394941106, 11665268.514318293, 11665538.88459844, 11666531.25, 11666830.473082704, 11666905.311987659, 11667506.191202579, 11670328.157690931, 11672012.5, 11674246.977770058, 11677543.75, 11685315.382099167, 11693904.252602592, 11693913.72523525, 11701722.687697079, 11702433.70238604, 11707979.6875, 11709510.431218628, 11711466.339072473, 11712299.960626382, 11716733.784559725, 11717059.96813094, 11718998.47461315, 11721661.460789476, 11721713.670180703, 11722691.251878558, 11723494.297862254, 11726993.790713605, 11730350.0, 11730352.892464817, 11730870.873112239, 11735086.28858121, 11743232.8125, 11747436.876557792, 11747899.838587401, 11748796.011783522, 11748797.8415802, 11752307.733269854, 11762040.213197922, 11762456.320982656, 11766809.847482182, 11767024.95197246, 11767847.311315097, 11768232.431149941, 11768729.381362313, 11770964.0625, 11771362.5, 11799322.413064573, 11800479.6875, 11802043.551689828, 11803662.83498937, 11803666.97686269, 11804820.41457572, 11805637.5, 11807121.673764924, 11807131.64604544, 11810110.9375, 11810252.516474968, 11813601.011391968, 11817235.285055207, 11818204.711702053, 11818670.581941573, 11819502.456537625, 11819518.532420712, 11819977.192026425, 11823423.137045275, 11826253.584050966, 11826262.64200081, 11836469.770282792, 11837861.227453059, 11838500.0, 11840973.517801823, 11848420.980863173, 11849158.064101465, 11849518.578272313, 11851240.625, 11854638.213618154, 11855264.775113294, 11861203.867439386, 11862718.005050613, 11879123.258855205, 11882607.023091545, 11885013.534939742, 11886403.125, 11891513.913511489, 11894375.785114475, 11895823.4375, 11896949.939463265, 11899438.925213845, 11904428.405943016, 11906252.886308905, 11906279.123036087, 11906279.592363678, 11932480.48582081, 11933103.057154663, 11933129.382322943, 11933607.8125, 11938084.375, 11938124.54868778, 11938337.5, 11939007.8125, 11942148.4375, 11942570.3125, 11943890.625, 11949318.75, 11950425.55961304, 11954123.955127371, 11967712.51602519, 11984147.400599526, 11998102.145341592, 11998440.688169323, 12003535.09044335, 12003643.562584821, 12003880.92697216, 12006603.295409188, 12006666.224627413, 12006818.280884814, 12007656.063976634, 12021760.904552123, 12025167.1875, 12026949.351648882, 12034604.825674282, 12036595.022035923, 12036644.718651596, 12037748.629993232, 12046750.480471171, 12052907.079733664, 12065759.22314647, 12066550.0, 12067487.860439241, 12068535.095491441, 12072198.31479832, 12072784.063159177, 12075505.769797347, 12080189.0625, 12083563.44953464, 12085642.1875, 12093829.193722008, 12095408.424948497, 12096765.275814628, 12097687.283872547, 12099518.427992776, 12104641.462846002, 12106722.425311297, 12107981.419838626, 12109543.675004374, 12109555.5617175, 12111009.40173332, 12118956.25, 12130685.9375, 12156460.963113505, 12158903.861347273, 12160021.875, 12161343.516034434, 12161435.442353878, 12162293.576249557, 12171394.532892589, 12182315.920075202, 12189519.792895334, 12190764.1994611, 12228982.8125, 12250967.39681442, 12257127.38689808, 12260809.874362111, 12277119.85655168, 12283352.550303686, 12284861.559482442, 12294715.167966884, 12304309.03100434, 12304940.136731103, 12307857.8125, 12307876.21857812, 12312274.954218056, 12312325.0, 12319710.193869041, 12333271.139402378, 12347278.705457788, 12350305.702439817, 12350903.774050955, 12352714.57503582, 12363439.0625, 12363475.027313102, 12382578.125, 12385620.09127635, 12395413.514889034, 12396338.908420142, 12404720.189026555, 12414685.262414258, 12429753.50300578, 12430477.592767531, 12491428.997429512, 12491487.455662034, 12505950.0, 12521524.194384303, 12522096.733955272, 12535976.748088233, 12545199.56781789, 12546245.514236271, 12585963.532921545, 12587055.081408978, 12587098.213551808, 12601330.047961889, 12623849.354517985, 12629663.63298652, 12630808.627519049, 12631189.776683358, 12640764.732623996, 12640784.12392131, 12642096.475688936, 12644814.0625, 12648048.011604825, 12650172.536429618, 12650685.788957292, 12650790.625, 12659515.579054026, 12687373.23501867, 12691241.793954123, 12691435.9375, 12694668.75, 12704594.864356427, 12725418.96091594, 12729585.923266191, 12755124.684385259, 12765718.094615197, 12766621.71747, 12767161.154606812, 12767207.4480891, 12778809.546092963, 12784032.8125, 12785545.3125, 12788303.125, 12788700.0, 12788746.79849205, 12789178.669848654, 12789427.616552504, 12800002.374935947, 12806134.212748306, 12826503.49843812, 12834753.649074761, 12849625.910309618, 12853181.06132268, 12856142.974524688, 12859390.97912963, 12861295.3125, 12861796.738231504, 12862332.296256898, 12862547.330079835, 12862806.167456886, 12863685.9381652, 12864459.375, 12865390.625, 12871498.191598183, 12872913.646558296, 12874745.59387956, 12882815.625, 12894935.9375, 12896462.768712457, 12897845.3125, 12920970.234010998, 12922365.25815487, 12923056.320378527, 12924702.32949444, 12927113.808847269, 12938816.046170648, 12955906.25, 12960581.25, 12961335.473904584, 12971468.75, 12982040.969907956, 12984873.4375, 12989683.307689022, 12992551.296673022, 12995479.6875, 13004751.5625, 13006326.5625, 13006348.4375, 13016577.63930378, 13035602.952919273, 13040328.60738103, 13042052.982904142, 13046184.146535974, 13046518.75, 13053731.042401625, 13071697.86593385, 13087968.440281667, 13101690.57579928, 13102119.088868024, 13107442.565876858, 13117316.501595251, 13120861.104912566, 13126672.090813775, 13127017.998469863, 13127883.434762586, 13128695.3125, 13130914.41729843, 13132032.776201544, 13138273.227820558, 13142433.123649413, 13147455.761792997, 13154959.424166335, 13171745.332598332, 13180323.354040999, 13214886.312817538, 13215305.013614299, 13231551.5625, 13249223.076949947, 13252125.051517086, 13263656.25, 13271519.978331707, 13285523.1385595, 13290871.863095077, 13295278.125, 13295300.193188995, 13312300.607757116, 13313914.01212709, 13357484.126843406, 13366547.027387178, 13368765.691720141, 13385956.25, 13386588.673917877, 13398531.333380474, 13406261.603192013, 13415476.465965934, 13482379.010671396, 13490359.342487063, 13522149.751474824, 13524280.220960297, 13524679.633438936, 13555312.792130562, 13573266.601688737, 13614512.5, 13630402.222515922, 13633070.3125, 13663222.470951432, 13678840.625, 13678887.108327253, 13721999.292205464, 13725011.222197175, 13725019.16909765, 13725500.0, 13725791.521667438, 13725820.872642437, 13746900.0, 13747684.375, 13792454.6875, 13794546.984439723, 13808651.80619715, 13822718.75, 13970953.879051354, 13984749.610296266, 13985732.8125, 14049518.555861633, 14074421.473048678, 14197527.130071469, 14239943.910490653, 14382580.759306889, 15372881.202261753, 15377288.481698131, 15392935.885784997, 15418070.241272867, 15481532.8125, 15528969.927964639, 15562996.770062903, 15594319.589288566, 15695528.125, 15862687.369659513, 15908140.806831326, 15934082.8125, 15935221.875, 16029304.6875, 16029330.713078398, 16034444.10512927, 16049114.775341185, 16069709.727653597, 16074374.445150739, 16089023.4375, 16098877.835589826, 16111560.213344913, 16131173.669604735, 16141026.911057696, 16156287.894219257, 16171649.791964475], [10.46005308030074, 14.904816831114386, 21.509579953826076, 5.315571685515921, 14.641153547339282, 14.551713151450613, 7.477746289301568, 20.004992277642685, 178.2872063856559, 59.93966451642283, 51.362576770807834, 13.1035239783946, 34.195371794219966, 6.353521381156771, 11.876391786882824, 49.29720514225406, 132.36260141018587, 53.10963427763747, 18.161256474081533, 7.655510267202092, 11.465416502899068, 46.9817648158792, 21.43463882913759, 40.35769116738798, 18.76656052722761, 40.98496199447338, 19.002676457902197, 28.500643229868082, 37.34412429749422, 36.900183409281084, 6.132065607959456, 35.5600762244857, 9.07236688080641, 65.35378508178178, 13.5901635727614, 47.518673981573706, 162.00679272838605, 11.112131036531967, 24.882819981192128, 8.818512170568809, 22.782783631012762, 33.151835487243325, 8.252462851111957, 22.97597363647257, 45.00266508243338, 15.39863943033076, 7.941796522170276, 19.47197065545811, 39.32959418108189, 20.28250527145475, 10.23476704601325, 61.29166591925012, 6.520803373945005, 6.557154839231342, 156.23966800287673, 35.28726230336729, 67.14434771748074, 9.389215692926406, 13.767004199774586, 10.05756550785566, 29.141779714716762, 43.686308098672754, 6.041208707818841, 10.40283654052718, 39.58016733270071, 7.209839204682123, 7.575339423011128, 9.875177979916275, 21.375157601033006, 26.410988076161814, 10.6285540427579, 96.88147833538343, 8.9662689655524, 15.183393585622014, 23.48732289220043, 88.3852572723313, 115.47603765312692, 16.196714113826587, 107.40044113737153, 80.10397482998351, 32.35797879642434, 8.967979769757086, 48.2424149815306, 5.634498105925283, 5.531920066023435, 25.917548565327237, 5.466336949612975, 100.60185238949083, 121.646224449383, 44.386850200776806, 28.479024807073255, 164.59426231862466, 45.88221007666563, 42.30776035408241, 23.155262788130653, 64.31043972340991, 32.77808818872232, 9.526361378626401, 21.826275585055114, 51.28072073653459, 61.93652129939126, 25.12061613846303, 48.4234809169319, 96.70928611494756, 61.733659847068594, 45.302781941492476, 11.415913640130704, 51.03460121989626, 61.56974874699588, 13.431149996742056, 23.255589646389794, 11.510755110375927, 60.410657109469916, 58.88748526514357, 5.141110712318085, 74.57940042074148, 22.05466600158739, 51.29357868028993, 9.473697022992809, 7.481815424423114, 14.03379984206384, 88.5330791497979, 41.56249942595854, 11.665106763904873, 14.388513035873158, 21.800541174700545, 100.96898633824384, 11.02977753813278, 33.29011745919081, 20.394383983071645, 63.55380257386972, 43.919357458252804, 23.01024424148823, 52.93824705598195, 75.43661023465496, 9.888737247731447, 10.475180869403623, 12.89789137394039, 86.39553385485206, 96.41423892357255, 46.41316530770811, 37.556570428969856, 49.45516408339277, 6.617917066333705, 54.147849956046876, 8.084325092096178, 35.802227558745834, 62.76021087535632, 56.88338260245663, 130.6711572084285, 40.56750600212464, 25.49335405251625, 6.576829047498817, 52.736603576252236, 61.237742416801304, 48.20686758601432, 6.074090595487571, 15.940573885565586, 7.151321960700783, 16.2914429290514, 54.63575849855518, 44.69204191866788, 32.98275736333001, 25.995671096857322, 14.31349780202911, 25.033433646329264, 12.407688834811019, 6.197103995513198, 11.698890112459546, 11.242860867452213, 61.46559166622376, 66.70496444239848, 29.360754381563762, 14.196409422614167, 173.43949869419188, 205.21444123748222, 17.427440157357893, 124.60847524303543, 42.2119156838969, 9.487726679291303, 32.12754485945547, 51.243098127498214, 52.13001059678446, 14.367510538450862, 21.861735304659298, 13.366603473479222, 54.168829738192784, 72.72168227726169, 70.9424489442474, 89.08706179555452, 7.373313596795591, 135.87841700917642, 13.563428602116813, 52.294007630488146, 6.03574452193149, 32.46112106928729, 80.24642603544508, 41.77935585724417, 62.44800718780023, 85.9149071529233, 95.6962602770132, 25.3958146126336, 9.296729245050015, 7.271880166555046, 89.8361150902715, 17.639128628348196, 127.51661966346144, 61.41548376518578, 91.77511047191413, 30.09324106604355, 17.436145515896865, 30.82023086744954, 5.562024721678407, 70.33998509823505, 14.693376565335392, 23.424056867223538, 27.749510271548267, 52.812905721240064, 30.3098036361471, 99.51020235291591, 5.540517540550961, 31.772223427695096, 45.8965926027949, 7.374882944576174, 6.086800553258134, 15.785608154027267, 8.749506438162307, 7.301604934364692, 77.32987370266302, 23.555715244751187, 21.8786219800095, 8.387161530518998, 38.02464967434318, 80.2738366049172, 11.759039832530886, 23.086243794466412, 5.771366899993589, 39.85785672538417, 39.12722844289006, 92.91403045781837, 156.85731100439705, 94.66225747566423, 7.429430069475531, 67.44153409041668, 74.59239359829178, 114.4623992148571, 10.051727391426644, 28.113319219764335, 19.331010790221494, 16.11048730103169, 51.794892973848604, 41.95453607240488, 54.95996210986396, 7.5252960554845805, 35.53186933286451, 99.96477038609784, 9.610001022015487, 51.56846671376123, 28.982923462512446, 86.05114115112606, 18.51965362359951, 25.195128065048458, 53.13692907005788, 62.34127565818207, 59.715944529763384, 9.25171666538587, 51.6489905775533, 46.17152996801633, 18.453163735074376, 25.287592676196066, 43.787470420982025, 83.66012110656628, 46.6402921641026, 19.34364728420873, 12.090490520501456, 6.240278905413386, 19.012242862990725, 6.416142371763144, 30.446169598231187, 14.139277350873845, 27.978337069651904, 71.59450683893728, 66.36487197355532, 116.90839376226157, 59.17809813705439, 29.775401455902266, 70.90357127638957, 40.04808196076903, 7.611867593801878, 98.40674967720241, 26.92596270771244, 14.497536487310022, 16.02166294164095, 144.06385216261046, 48.62399427477424, 73.47949365020442, 8.594823383759417, 74.05095983615143, 13.546102842927178, 46.63661713175448, 69.56341089096935, 21.91789367967564, 11.735111628696554, 22.142625269171518, 32.002270048779195, 42.204814344109835, 22.899857877829845, 19.031111692547118, 59.49666932087761, 44.627803443600456, 67.05204489559316, 17.508628846794508, 30.39458265216645, 43.345792690808075, 21.022931173373195, 29.070446501135432, 29.99797570216571, 52.102735893667, 6.377865142572682, 5.266000795281109, 59.62495941644705, 27.277400204586463, 206.84018002133402, 45.02901706739867, 68.19244353274112, 42.37770458631099, 78.30816841088722, 185.95939386039947, 70.62829096274824, 115.16122749467726, 40.218162398003315, 32.60279394256348, 16.540756214495634, 55.50737314910614, 39.390508513644036, 27.728811880335755, 95.15657115800794, 12.956759891372322, 71.7314885580177, 38.210511165429956, 5.506273196787376, 50.67208622977007, 88.07391667117014, 15.193828212766471, 112.62899174153308, 37.5015304182949, 22.018191016381163, 68.47414725799158, 28.32038456990839, 7.016172240331582, 28.445266863418382, 47.20740459649353, 21.185522844303225, 6.967326293996054, 42.11691045642713, 26.45409229157432, 155.27544494740954, 5.177086943933869, 87.11463274885384, 98.15602438969994, 41.45646777609386, 40.82943453484009, 6.134975346618311, 110.64372977868027, 62.2339636612124, 21.353354503218227, 17.690927850082105, 41.1880027080968, 18.59438752306608, 5.251245528259005, 80.01210067352088, 36.00276588777652, 18.81077143953084, 55.118472665717626, 30.039363902200748, 74.22049446597399, 46.90936199130296, 9.908559986058348, 17.11277845915571, 109.16717911069404, 37.37516937285763, 35.51583663351791, 62.192104250686526, 6.722930732293959, 41.09043166829201, 107.44944307224371, 76.4716792331794, 13.31265237845868, 26.0619587964073, 14.60357168296921, 15.270581374537088, 5.547779206414565, 74.15917302598947, 17.34287689319821, 5.524041923611538, 91.61268523532743, 82.23288021946027, 28.66258016231189, 25.083309343602085, 18.305894510912008, 11.865352898185584, 29.852207572898763, 22.33450854734035, 13.257051405218004, 49.726473677150686, 54.601991080018976, 71.2597066597719, 84.01063941639904, 35.49505355652587, 93.79843519859801, 15.68269031646279, 44.26204534874868, 22.67434708861132, 112.24608451983018, 13.436237123912555, 6.353399137402942, 21.017091219679685, 114.37584594732776, 41.89707069831421, 6.942584154398579, 16.50850163933523, 18.57443256995602, 68.77056582430806, 98.42990465030982, 52.323462314432355, 6.284595118619792, 135.61081354767194, 53.503255212290355, 32.55315815662022, 103.9825048778239, 14.441371419713642, 149.4500938009838, 60.55300531501252, 159.73576478759458, 20.462234155829286, 85.40999915027011, 44.830266160964996, 69.8709751443631, 12.600695424099838, 49.70777220682642, 31.134932741801837, 44.67853744047713, 12.913321184739036, 44.410528568998515, 80.89785451680575, 124.53323675450727, 41.134293649110944, 58.775219974976, 101.97560395560843, 25.196936930288974, 52.84964313625933, 53.20986339075504, 23.607944156480748, 40.20613385199743, 7.09053907774558, 9.161611748811067, 100.40686317953406, 76.8505708126323, 10.070680414025194, 38.680465525076265, 26.07296711513674, 95.11151012711966, 8.357730495497414, 36.605707419413264, 63.8598890055455, 12.377090374751722, 66.75748662706995, 14.968322651239426, 27.804981391383635, 22.10676060245977, 8.993932220826862, 36.9085924589413, 5.6377866380751565, 113.43088218598271, 84.45790557888111, 22.394274274559887, 58.54920632364563, 22.71435657008975, 44.8913671637996, 13.911799463428904, 13.728281485434376, 19.150535145545902, 51.377452720530435, 36.68975611028691, 30.645004890194983, 9.455885651841506, 15.233478004502677, 55.631317699723354, 49.300862500962175, 161.59919467714377, 24.980748227581522, 9.726103225243985, 18.370513047381255, 37.172800474716496, 11.421991218353877, 38.323132887689454, 41.760676636721925, 9.323056115934662, 10.310466763641115, 78.72991463916759, 24.914633917131162, 8.814065689379412, 78.53581305810236, 14.06509920815876, 41.93971447769981, 69.69164099130191, 11.477325140736363, 26.75752412790964, 10.215085828001309, 69.89188717479753, 113.18339962873547, 130.03007456240235, 27.260158710473284, 14.253568695148076, 94.48413618988229, 105.55620901609637, 29.854065329518868, 102.95286200566437, 31.10033772785273, 6.849499239611752, 16.09435204052012, 34.263079763037496, 152.05090160172193, 32.91372536810881, 10.795608878117786, 21.11439496595769, 26.57605755229211, 234.4321639864748, 18.77366329613951, 77.8789585018097, 40.14073192674408, 12.15367125819885, 45.365951030545986, 39.42761410472773, 56.64166862771388, 7.707224747318835, 45.10793019880295, 24.015195898340927, 54.35573414782981, 7.57506984591888, 8.84424750563413, 7.080166757187136, 76.26817418251031, 6.238563157711726, 9.139766552876395, 24.02479033910739, 54.97474384929673, 8.750830569022959, 5.02608325252503, 91.51394733406534, 6.7700506558429225, 56.87981644273842, 29.39628482227317, 28.831415817558295, 71.99769897983538, 24.221276936329392, 6.435490489730579, 10.128202661889697, 12.87691591399845, 42.77654693933018, 119.9671510729133, 30.279461042153734, 41.45253137717987, 20.94816474761589, 12.201378308452913, 11.18621682671237, 7.266155582998123, 5.445175995065317, 14.214353619881635, 8.599226627305159, 39.8440815802911, 42.757769442579985, 24.652540503521674, 48.54345280295502, 22.575182363354173, 90.62945528068869, 67.0776418832133, 39.01717021678313, 41.456834447734494, 83.48234554651404, 14.647882435466503, 33.68568462216086, 147.80540979695195, 6.328797527910687, 47.703484041433114, 144.19888872822813, 8.094708284892777, 33.29828771044275, 57.35404278575308, 33.92769941592709, 13.471009960733731, 74.09015515088883, 5.694252196734011, 11.69505950570191, 24.963533392118325, 69.2989177756227, 107.10465511599463, 8.18295246096663, 26.920217970702648, 80.408428396774, 13.038334115335761, 61.124891169077564, 64.83192814811643, 9.351951560310564, 18.65260142067794, 12.887190726114909, 7.526394040069647, 117.8542058627798, 6.7049697792421235, 5.305480799182547, 15.511858167455491, 158.003421471843, 15.14468356362237, 100.3086861773746, 15.256678399609854, 6.160474786050645, 9.553610452182832, 48.341340959699366, 54.47442227717069, 66.6248936672985, 10.016615835763426, 18.80126186023898, 51.17998097565463, 14.301679867052636, 77.52751029771365, 8.921892657072155, 54.89846921988301, 17.33076715646697, 102.10316010293691, 38.450140532371854, 128.13530916122204, 59.607154563478474, 8.44095068429451, 24.566795988847897, 10.346857269145369, 10.982188901561608, 84.31591920415146, 27.25081613684647, 5.092496767658493, 45.78943359716188, 45.102119009829224, 18.9284950636643, 51.104402998203305, 111.79485160976334, 49.53916811042245, 81.03624590685057, 82.52736583402798, 69.3828527196381, 38.15212096690889, 78.56843966412396, 57.75017315547217, 66.56362324616298, 19.915158544044868, 6.7761837034514665, 12.46642960681854, 26.503109069968744, 23.4609830696894, 200.6357994925697, 146.34723781013432, 26.060271222565984, 82.80220866999161, 27.889096638261737, 69.91641671835687, 187.24765842954773, 14.764638400887188, 23.873318930481776, 10.311096077659993, 62.030587279356716, 11.972088670592324, 12.730767367373948, 14.744926948739588, 75.13660392289218, 16.51979112769956, 45.92522289750906, 26.314558628522253, 103.6013029108133, 110.19349901986118, 128.84037061445736, 15.108716975852179, 107.45754482563311, 13.343119011495226, 23.62001218040819, 35.52648172839175, 65.78981482359202, 54.23635156563837, 19.510966824017967, 54.86017367432424, 15.63093797440105, 285.41646802802813, 45.69000009228574, 15.331499618514034, 52.43298035522011, 34.27608752401803, 15.040043169067104, 6.919990493655217, 71.6797750061305, 24.632552222249586, 31.53428674693425, 69.19490509983791, 94.91574108996396, 8.449035442929048, 13.272672253935308, 180.03725244602123, 38.36112770745725, 15.80421830495635, 12.718044972855013, 57.66727139151401, 124.32413496646772, 5.890872749530303, 141.78366850821695, 49.21254345002892, 24.862206126987317, 127.11754308219128, 108.01070954329697, 191.72434928058928, 72.02571590310596, 29.860304521020527, 5.574188017950052, 33.00205110082356, 17.737171289568295, 34.43730538522372, 5.111448490476592, 87.3112673130835, 81.14294925463486, 25.06552323976976, 30.416920697069553, 24.81200644198517, 110.79433240950424, 5.525648997796797, 82.54634276357329, 34.207482291239515, 7.2244942474555, 50.393069818758974, 289.9863771281537, 9.294192209314993, 59.506221368054995, 148.01703662683747, 14.429850073122667, 7.1711069596050905, 129.74907338994228, 20.36961735410386, 8.920454949366103, 5.552508556854841, 6.06511613050196, 29.70379490382412, 6.641138105415583, 59.772044043957976, 32.0481623617099, 60.902347472461955, 95.8247708695723, 66.64354265054752, 22.42525927588848, 83.02261752987962, 76.47169997903089, 18.248775239295398, 26.681639978982467, 84.20116422305013, 149.20322021001482, 45.67598245106359, 16.11656586360676, 109.41111175694266, 150.36828264807724, 6.419062870226426, 32.98809616617845, 40.07854015769477, 133.3801368369283, 13.186270935428443, 25.758832274902964, 262.3294823261272, 56.257633196179754, 6.2890594725513616, 58.52036261183278, 95.66214036585107, 118.29535592895849, 56.59520466918407, 9.26753851237654, 17.523904749402295, 21.409243826750533, 12.689661624799005, 111.58000988517797, 67.07812222020456, 5.986382592171176, 120.56743583247504, 30.441937892191945, 149.67243871328762, 27.94921410683981, 39.75178348026374, 6.736505556424783, 197.27161225840402, 14.063270527564626, 6.07955731649008, 24.509736987045585, 34.843079116447335, 105.9847749437468, 36.87142771067538, 9.002835603148894, 77.89318863267114, 91.03190717933538, 55.24562061211194, 70.16081829807507, 39.834999718723076, 35.282617266385046, 22.932858627925537, 13.083635004557983, 36.02916338132999, 24.3613847957462, 20.7625135836459, 111.11582473938283, 74.89439280824058, 157.00914179930862, 79.92202177784577, 16.442159872981588, 73.87337579656358, 62.14438589694732, 59.148395646215235, 68.094066942419, 38.53874835114903, 74.70283607632217, 47.317640502741654, 5.778896196151821, 172.2530912279364, 69.49620086919462, 26.001595138787152, 24.709647298132136, 52.415221573384656, 69.41312245627446, 76.84416934341235, 23.627392125530026, 125.59241210530469, 130.8182542375304, 5.973048415159237, 35.23535088861591, 161.4597549767803, 40.03074657455956, 142.01891556901927, 26.863609222982504, 44.89134253814364, 58.96898944948762, 113.93732092138178, 9.160184686044985, 75.63429296358292, 149.2098021182579, 85.2168357443131, 256.4584291439799, 10.999713587102859, 22.915482147833956, 7.725605761112345, 79.79078745708637, 41.31199146266421, 5.04241315796255, 52.74131554358093, 32.29531434318677, 10.690523101746015, 57.08718483670773, 48.972701567198854, 11.443724926437923, 5.9941398408934905, 218.0114386269908, 10.760779063747075, 6.794315663464759, 15.202063307290153, 108.53077834817229, 120.22234106455603, 6.780267788600447, 26.185225598766877, 26.75143978064268, 7.9583949057945675, 17.24768056533221, 33.954093994392, 47.04282213920154, 27.498754487763403, 41.33288818306591, 20.6777671441621, 59.0258496188467, 43.87522345295311, 61.545955207928706, 42.87074019796485, 43.07556484745181, 10.566121665295805, 26.321645473810776, 6.303850598237803, 9.712834439690447, 80.67118243646796, 52.95424603686472, 6.148352706654838, 53.67906388039738, 60.92083116930623, 75.68472294765404, 113.31443379749983, 155.35045909199408, 36.821991516969675, 160.8859167163571, 56.264180861206974, 78.00386262070087, 30.90678716430718, 12.780285019615246, 16.233017842299414, 15.590668832209131, 29.917501354037608, 37.323004688903474, 25.80216880840069, 30.995986736378207, 39.21098722866333, 30.919047224710905, 22.727814482073246, 10.967551429868204, 109.9834275188002, 83.32517340778499, 55.833565750028384, 15.68302399677503, 43.73613658005283, 52.28963701482757, 48.21845633612445, 5.27478532295761, 16.49313810919062, 37.4853504266769, 86.83622237589096, 26.88375651187652, 43.64692058661242, 10.703590272493596, 19.51486118811333, 11.510975767421295, 15.340383725070527, 35.44902153180603, 57.20358509925974])
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);
([6545529.213785079, 6601443.740339845, 6662226.602305766, 6732844.522830824, 6817765.543901257, 6907381.143830056, 7325765.2641742155, 7759829.056142265, 7802295.032793387, 7834572.64073461, 7868227.138216119, 7933558.791688724, 8032593.463836384, 8073361.591505273, 8092861.36892458, 8105679.057481564, 8106128.95183724, 8106993.75, 8108664.7105028825, 8112257.954749868, 8113650.262079768, 8113853.125, 8117096.369119483, 8118103.11295699, 8192141.560404632, 8194721.875, 8206549.852044636, 8208863.197536281, 8209642.047681872, 8211166.8731837515, 8211170.724937779, 8214995.461090431, 8237084.0518939905, 8301608.976428049, 8325546.834221977, 8329060.448176033, 8329093.731634366, 8329432.8025122015, 8330336.141706665, 8331338.3865176905, 8332301.182745698, 8332395.3125, 8332757.747427058, 8332923.852816924, 8335701.858132108, 8337334.874439936, 8337524.285666104, 8337839.834648788, 8338199.886158746, 8340086.02878969, 8345100.904741182, 8358002.910713149, 8358550.211738515, 8360208.245350861, 8360759.375, 8362625.0, 8377450.0, 8435392.255089872, 8435463.715977885, 8441497.624795098, 8441640.45383395, 8444071.2346073, 8444588.599189525, 8445291.668337902, 8445574.794458464, 8447351.888398295, 8463628.047201036, 8470313.079960566, 8486269.05751779, 8490886.429052072, 8557391.842824403, 8563310.238943413, 8588580.784117008, 8590805.300327597, 8591990.304094246, 8594954.340012152, 8618866.668761615, 8618886.231965266, 8624191.689141514, 8625593.439316887, 8626881.764455257, 8649429.223971142, 8650964.538331347, 8653497.23889746, 8702334.351065807, 8732760.790100006, 8755150.106053265, 8757763.728888875, 8759719.464865796, 8760619.811208945, 8763292.18516511, 8763292.1875, 8765500.0, 8767821.64236096, 8767959.026907006, 8795708.370933587, 8797537.5, 8801548.659406807, 8857623.644137666, 8885002.168992119, 8931637.5, 8948024.091380334, 8962002.827475922, 8969404.098924745, 8970245.078404859, 8970995.3125, 8971434.551879466, 8990065.625, 8991261.143546656, 8993468.908154646, 8994562.160221824, 8997750.231978243, 8998873.4375, 9010514.0625, 9019895.802334512, 9027206.25, 9029730.804510111, 9035000.0, 9045142.351123119, 9045557.463825863, 9045923.630364526, 9046205.293393629, 9047670.3125, 9047795.631100819, 9048142.034433374, 9048476.80035415, 9048785.947103929, 9049496.811431251, 9064260.280110447, 9066629.069822852, 9073779.6875, 9079368.75, 9081285.388403768, 9093103.652019558, 9113043.909975333, 9122067.780470518, 9131460.475659497, 9171338.606156098, 9172570.3125, 9179415.625, 9198651.040966596, 9203243.416915242, 9205006.25, 9219766.411712954, 9221012.5, 9225635.534399133, 9225745.3125, 9226133.51447952, 9229023.616649669, 9232575.0, 9243094.852259452, 9247428.312908007, 9289077.112029158, 9314314.726417897, 9356561.940966213, 9360939.0625, 9364352.989204347, 9371840.135847941, 9371862.612425864, 9374944.064178856, 9377998.4375, 9379875.936855169, 9382115.14982822, 9396960.795031626, 9397510.762362551, 9400763.790998338, 9400784.611448744, 9402002.549305363, 9403455.024575658, 9403457.884943636, 9404376.5625, 9405761.46749293, 9405914.0625, 9408407.770410601, 9418831.141779719, 9423271.550803911, 9425867.813507302, 9427284.478111198, 9427374.674525179, 9429082.991766201, 9431965.625, 9432021.875, 9432260.738258172, 9432262.340787336, 9432274.536176618, 9432602.846919704, 9432880.020945568, 9433003.125, 9433318.221633919, 9433404.217507947, 9433414.941164924, 9433717.485008687, 9433965.20944518, 9434406.25, 9434491.66853627, 9434503.813545281, 9434704.6875, 9434822.20995872, 9434828.125, 9434942.1875, 9435025.604198152, 9435509.040992258, 9435520.271206118, 9435630.230430307, 9435657.8125, 9435839.62672939, 9435859.18247504, 9436091.694685213, 9436098.299875695, 9436182.8125, 9436563.029858764, 9436611.011853868, 9436636.519795697, 9436651.971010687, 9436740.558339758, 9436877.38446657, 9437382.813397799, 9437571.875, 9437745.3125, 9437757.8125, 9437759.527289007, 9437777.851352466, 9437967.39339713, 9438185.258683223, 9438515.612502528, 9438753.330132032, 9438861.721077124, 9438911.581195638, 9439590.625, 9439960.62207716, 9440740.1627974, 9441598.622983204, 9442841.666150799, 9444425.292663345, 9444469.119553005, 9446005.14197341, 9446172.844474513, 9447012.5, 9450673.4375, 9454821.875, 9454986.440922506, 9455357.530089725, 9455485.429943906, 9456978.125, 9457161.705962108, 9457779.322191294, 9457830.432026725, 9458060.916366607, 9459089.87700894, 9460407.821022231, 9464142.1875, 9464836.849501545, 9465302.195138903, 9465729.466735832, 9477117.1875, 9478935.739605473, 9479725.312078428, 9495448.4375, 9499450.0, 9499516.037360229, 9500297.589766206, 9501476.185267553, 9501576.5625, 9503620.274694307, 9526613.904179893, 9529504.591730673, 9536462.5, 9537227.686359692, 9539370.678346911, 9541635.74898897, 9543403.125, 9543963.479973812, 9544760.265846293, 9545089.878665935, 9566414.601978991, 9576412.283920541, 9578401.31579933, 9584643.281961033, 9588043.188141929, 9596770.084800536, 9603180.778193986, 9605610.737446869, 9607126.478778133, 9611895.374009397, 9623030.122808771, 9623271.499024883, 9623290.329741282, 9625993.402314162, 9627225.023797566, 9629996.275939662, 9644870.226136304, 9645815.340879427, 9648001.46177157, 9651005.269882688, 9652051.5625, 9655262.5, 9663489.603204867, 9667303.125, 9672744.890929077, 9674814.0625, 9676816.011572262, 9677023.031553837, 9677025.35669642, 9677832.375236703, 9690299.351021303, 9690727.210483301, 9691023.096060744, 9703917.961847419, 9707123.226414874, 9707942.1875, 9719890.625, 9720425.42306955, 9721730.703852907, 9726236.738438789, 9729999.451539755, 9730955.54459186, 9732982.478701603, 9736371.875, 9737835.203441843, 9738841.3534429, 9740083.773235837, 9741128.478713376, 9741784.457737312, 9742604.6875, 9746357.8125, 9747543.75, 9751997.449634438, 9754362.410861436, 9754582.8125, 9760983.70510806, 9764226.251113283, 9765085.236234723, 9766337.86459973, 9766383.748807963, 9769110.686924092, 9769869.84502496, 9770868.75, 9770919.671887003, 9771246.875, 9771643.825479457, 9772663.497492108, 9772757.557257941, 9773170.315507691, 9773367.252160383, 9774263.16861902, 9780753.797575934, 9821015.176074697, 9842382.8125, 9852063.579718066, 9853594.8962652, 9891549.773014551, 9896141.735449243, 9900335.406103969, 9900725.94183369, 9903371.383948818, 9903517.47241881, 9904158.110073226, 9908627.058942651, 9925148.4375, 9926445.0137219, 9939159.375, 9939535.782134289, 9940564.50152192, 9941115.625, 9941650.0, 9942491.67337611, 9942577.70440662, 9944264.227969505, 9947063.212841185, 9953679.904006047, 9958560.9375, 9959659.375, 9959668.436231282, 9965295.3125, 9979183.892351035, 9979589.143452322, 9983087.5, 9991424.049376253, 9993962.972848086, 9994070.3125, 10015067.1875, 10037926.173215073, 10038112.5, 10038951.35580175, 10043943.138549088, 10047739.0625, 10048773.98275896, 10105865.004850542, 10108504.549188336, 10168835.845919825, 10234386.021237548, 10244173.736683426, 10312677.48105694, 10373066.132452358, 10377382.624846734, 10505564.0625, 10521202.80615693, 10529198.051162573, 10578404.975586843, 10626162.821523914, 10644178.550686719, 10656762.5, 10668391.390841668, 10668418.22461328, 10677596.875, 10749830.730311716, 10762576.336158361, 10767298.429648796, 10875609.333053276, 10883840.192975476, 10948213.609016996, 11151757.8125, 11160524.610786699, 11173220.202527693, 11299404.455005264, 11299428.778118704, 11327404.138420654, 11332867.404561844, 11333410.9375, 11333882.085358296, 11352921.041245257, 11353872.58494068, 11354700.312471772, 11355431.770494482, 11356962.5, 11356997.011465302, 11357232.148145527, 11363069.024270117, 11366154.6875, 11368171.233499516, 11371871.192846676, 11400900.0, 11401667.385170722, 11403754.588236444, 11405050.952554448, 11406632.782814851, 11412135.9375, 11416020.3125, 11424591.399408825, 11427416.139031742, 11438531.25, 11438978.138026679, 11440947.77955094, 11441187.331962826, 11448940.57226916, 11448963.207845934, 11455521.601381162, 11457763.37095222, 11458092.1875, 11459920.117868388, 11461113.792870315, 11465185.9375, 11477218.813439531, 11477919.771468744, 11483428.006083459, 11485238.854629679, 11485408.605975455, 11488195.3125, 11488372.636590999, 11489262.5, 11489702.987376181, 11489955.52887039, 11491079.81967093, 11491269.316853989, 11492766.729524102, 11493381.925117178, 11493442.1875, 11493888.08592875, 11494912.09131952, 11495871.43439221, 11496816.759470627, 11497315.91217744, 11500639.483836515, 11503382.8125, 11505615.625, 11506890.210830119, 11508698.954743141, 11510505.86428975, 11511813.079705795, 11513998.138291761, 11514951.203077454, 11517391.415535057, 11518128.125, 11518866.575209213, 11518898.826402118, 11519306.393639904, 11520620.38704898, 11521648.4375, 11521871.875, 11522048.218101023, 11522065.415555563, 11522067.620109076, 11522129.288961964, 11528187.17296505, 11528297.055454506, 11534693.75, 11534714.448949916, 11542151.552611802, 11543082.218131347, 11543309.005644133, 11543379.945199056, 11544802.442480206, 11545403.125, 11545555.560950374, 11545860.9375, 11547273.917260563, 11549045.525223903, 11553188.187430112, 11566474.325492933, 11570828.341470277, 11575864.0625, 11590032.300134156, 11592060.870179474, 11594838.869490411, 11594995.3125, 11595126.662241884, 11596142.14975528, 11601490.807095895, 11603089.0625, 11603600.415995216, 11608178.397806078, 11612622.270998115, 11613879.601216013, 11622389.0625, 11624160.714110233, 11629766.648679698, 11634683.1117809, 11637053.1077956, 11637378.6952132, 11639658.86344426, 11639910.831204617, 11641647.445234919, 11642433.430389158, 11644003.125, 11645330.865115263, 11648181.528639555, 11651168.693455769, 11655637.856510479, 11661164.0625, 11661182.733942738, 11661560.299101988, 11662399.840582158, 11664781.25, 11664840.394941106, 11665268.514318293, 11665538.88459844, 11666531.25, 11666830.473082704, 11666905.311987659, 11667506.191202579, 11670328.157690931, 11672012.5, 11674246.977770058, 11677543.75, 11685315.382099167, 11693904.252602592, 11693913.72523525, 11701722.687697079, 11702433.70238604, 11707979.6875, 11709510.431218628, 11711466.339072473, 11712299.960626382, 11716733.784559725, 11717059.96813094, 11718998.47461315, 11721661.460789476, 11721713.670180703, 11722691.251878558, 11723494.297862254, 11726993.790713605, 11730350.0, 11730352.892464817, 11730870.873112239, 11735086.28858121, 11743232.8125, 11747436.876557792, 11747899.838587401, 11748796.011783522, 11748797.8415802, 11752307.733269854, 11762040.213197922, 11762456.320982656, 11766809.847482182, 11767024.95197246, 11767847.311315097, 11768232.431149941, 11768729.381362313, 11770964.0625, 11771362.5, 11799322.413064573, 11800479.6875, 11802043.551689828, 11803662.83498937, 11803666.97686269, 11804820.41457572, 11805637.5, 11807121.673764924, 11807131.64604544, 11810110.9375, 11810252.516474968, 11813601.011391968, 11817235.285055207, 11818204.711702053, 11818670.581941573, 11819502.456537625, 11819518.532420712, 11819977.192026425, 11823423.137045275, 11826253.584050966, 11826262.64200081, 11836469.770282792, 11837861.227453059, 11838500.0, 11840973.517801823, 11848420.980863173, 11849158.064101465, 11849518.578272313, 11851240.625, 11854638.213618154, 11855264.775113294, 11861203.867439386, 11862718.005050613, 11879123.258855205, 11882607.023091545, 11885013.534939742, 11886403.125, 11891513.913511489, 11894375.785114475, 11895823.4375, 11896949.939463265, 11899438.925213845, 11904428.405943016, 11906252.886308905, 11906279.123036087, 11906279.592363678, 11932480.48582081, 11933103.057154663, 11933129.382322943, 11933607.8125, 11938084.375, 11938124.54868778, 11938337.5, 11939007.8125, 11942148.4375, 11942570.3125, 11943890.625, 11949318.75, 11950425.55961304, 11954123.955127371, 11967712.51602519, 11984147.400599526, 11998102.145341592, 11998440.688169323, 12003535.09044335, 12003643.562584821, 12003880.92697216, 12006603.295409188, 12006666.224627413, 12006818.280884814, 12007656.063976634, 12021760.904552123, 12025167.1875, 12026949.351648882, 12034604.825674282, 12036595.022035923, 12036644.718651596, 12037748.629993232, 12046750.480471171, 12052907.079733664, 12065759.22314647, 12066550.0, 12067487.860439241, 12068535.095491441, 12072198.31479832, 12072784.063159177, 12075505.769797347, 12080189.0625, 12083563.44953464, 12085642.1875, 12093829.193722008, 12095408.424948497, 12096765.275814628, 12097687.283872547, 12099518.427992776, 12104641.462846002, 12106722.425311297, 12107981.419838626, 12109543.675004374, 12109555.5617175, 12111009.40173332, 12118956.25, 12130685.9375, 12156460.963113505, 12158903.861347273, 12160021.875, 12161343.516034434, 12161435.442353878, 12162293.576249557, 12171394.532892589, 12182315.920075202, 12189519.792895334, 12190764.1994611, 12228982.8125, 12250967.39681442, 12257127.38689808, 12260809.874362111, 12277119.85655168, 12283352.550303686, 12284861.559482442, 12294715.167966884, 12304309.03100434, 12304940.136731103, 12307857.8125, 12307876.21857812, 12312274.954218056, 12312325.0, 12319710.193869041, 12333271.139402378, 12347278.705457788, 12350305.702439817, 12350903.774050955, 12352714.57503582, 12363439.0625, 12363475.027313102, 12382578.125, 12385620.09127635, 12395413.514889034, 12396338.908420142, 12404720.189026555, 12414685.262414258, 12429753.50300578, 12430477.592767531, 12491428.997429512, 12491487.455662034, 12505950.0, 12521524.194384303, 12522096.733955272, 12535976.748088233, 12545199.56781789, 12546245.514236271, 12585963.532921545, 12587055.081408978, 12587098.213551808, 12601330.047961889, 12623849.354517985, 12629663.63298652, 12630808.627519049, 12631189.776683358, 12640764.732623996, 12640784.12392131, 12642096.475688936, 12644814.0625, 12648048.011604825, 12650172.536429618, 12650685.788957292, 12650790.625, 12659515.579054026, 12687373.23501867, 12691241.793954123, 12691435.9375, 12694668.75, 12704594.864356427, 12725418.96091594, 12729585.923266191, 12755124.684385259, 12765718.094615197, 12766621.71747, 12767161.154606812, 12767207.4480891, 12778809.546092963, 12784032.8125, 12785545.3125, 12788303.125, 12788700.0, 12788746.79849205, 12789178.669848654, 12789427.616552504, 12800002.374935947, 12806134.212748306, 12826503.49843812, 12834753.649074761, 12849625.910309618, 12853181.06132268, 12856142.974524688, 12859390.97912963, 12861295.3125, 12861796.738231504, 12862332.296256898, 12862547.330079835, 12862806.167456886, 12863685.9381652, 12864459.375, 12865390.625, 12871498.191598183, 12872913.646558296, 12874745.59387956, 12882815.625, 12894935.9375, 12896462.768712457, 12897845.3125, 12920970.234010998, 12922365.25815487, 12923056.320378527, 12924702.32949444, 12927113.808847269, 12938816.046170648, 12955906.25, 12960581.25, 12961335.473904584, 12971468.75, 12982040.969907956, 12984873.4375, 12989683.307689022, 12992551.296673022, 12995479.6875, 13004751.5625, 13006326.5625, 13006348.4375, 13016577.63930378, 13035602.952919273, 13040328.60738103, 13042052.982904142, 13046184.146535974, 13046518.75, 13053731.042401625, 13071697.86593385, 13087968.440281667, 13101690.57579928, 13102119.088868024, 13107442.565876858, 13117316.501595251, 13120861.104912566, 13126672.090813775, 13127017.998469863, 13127883.434762586, 13128695.3125, 13130914.41729843, 13132032.776201544, 13138273.227820558, 13142433.123649413, 13147455.761792997, 13154959.424166335, 13171745.332598332, 13180323.354040999, 13214886.312817538, 13215305.013614299, 13231551.5625, 13249223.076949947, 13252125.051517086, 13263656.25, 13271519.978331707, 13285523.1385595, 13290871.863095077, 13295278.125, 13295300.193188995, 13312300.607757116, 13313914.01212709, 13357484.126843406, 13366547.027387178, 13368765.691720141, 13385956.25, 13386588.673917877, 13398531.333380474, 13406261.603192013, 13415476.465965934, 13482379.010671396, 13490359.342487063, 13522149.751474824, 13524280.220960297, 13524679.633438936, 13555312.792130562, 13573266.601688737, 13614512.5, 13630402.222515922, 13633070.3125, 13663222.470951432, 13678840.625, 13678887.108327253, 13721999.292205464, 13725011.222197175, 13725019.16909765, 13725500.0, 13725791.521667438, 13725820.872642437, 13746900.0, 13747684.375, 13792454.6875, 13794546.984439723, 13808651.80619715, 13822718.75, 13970953.879051354, 13984749.610296266, 13985732.8125, 14049518.555861633, 14074421.473048678, 14197527.130071469, 14239943.910490653, 14382580.759306889, 15372881.202261753, 15377288.481698131, 15392935.885784997, 15418070.241272867, 15481532.8125, 15528969.927964639, 15562996.770062903, 15594319.589288566, 15695528.125, 15862687.369659513, 15908140.806831326, 15934082.8125, 15935221.875, 16029304.6875, 16029330.713078398, 16034444.10512927, 16049114.775341185, 16069709.727653597, 16074374.445150739, 16089023.4375, 16098877.835589826, 16111560.213344913, 16131173.669604735, 16141026.911057696, 16156287.894219257, 16171649.791964475], [10.46005308030074, 14.904816831114386, 21.509579953826076, 5.315571685515921, 14.641153547339282, 14.551713151450613, 7.477746289301568, 20.004992277642685, 178.2872063856559, 59.93966451642283, 51.362576770807834, 13.1035239783946, 34.195371794219966, 6.353521381156771, 11.876391786882824, 49.29720514225406, 132.36260141018587, 53.10963427763747, 18.161256474081533, 7.655510267202092, 11.465416502899068, 46.9817648158792, 21.43463882913759, 40.35769116738798, 18.76656052722761, 40.98496199447338, 19.002676457902197, 28.500643229868082, 37.34412429749422, 36.900183409281084, 6.132065607959456, 35.5600762244857, 9.07236688080641, 65.35378508178178, 13.5901635727614, 47.518673981573706, 162.00679272838605, 11.112131036531967, 24.882819981192128, 8.818512170568809, 22.782783631012762, 33.151835487243325, 8.252462851111957, 22.97597363647257, 45.00266508243338, 15.39863943033076, 7.941796522170276, 19.47197065545811, 39.32959418108189, 20.28250527145475, 10.23476704601325, 61.29166591925012, 6.520803373945005, 6.557154839231342, 156.23966800287673, 35.28726230336729, 67.14434771748074, 9.389215692926406, 13.767004199774586, 10.05756550785566, 29.141779714716762, 43.686308098672754, 6.041208707818841, 10.40283654052718, 39.58016733270071, 7.209839204682123, 7.575339423011128, 9.875177979916275, 21.375157601033006, 26.410988076161814, 10.6285540427579, 96.88147833538343, 8.9662689655524, 15.183393585622014, 23.48732289220043, 88.3852572723313, 115.47603765312692, 16.196714113826587, 107.40044113737153, 80.10397482998351, 32.35797879642434, 8.967979769757086, 48.2424149815306, 5.634498105925283, 5.531920066023435, 25.917548565327237, 5.466336949612975, 100.60185238949083, 121.646224449383, 44.386850200776806, 28.479024807073255, 164.59426231862466, 45.88221007666563, 42.30776035408241, 23.155262788130653, 64.31043972340991, 32.77808818872232, 9.526361378626401, 21.826275585055114, 51.28072073653459, 61.93652129939126, 25.12061613846303, 48.4234809169319, 96.70928611494756, 61.733659847068594, 45.302781941492476, 11.415913640130704, 51.03460121989626, 61.56974874699588, 13.431149996742056, 23.255589646389794, 11.510755110375927, 60.410657109469916, 58.88748526514357, 5.141110712318085, 74.57940042074148, 22.05466600158739, 51.29357868028993, 9.473697022992809, 7.481815424423114, 14.03379984206384, 88.5330791497979, 41.56249942595854, 11.665106763904873, 14.388513035873158, 21.800541174700545, 100.96898633824384, 11.02977753813278, 33.29011745919081, 20.394383983071645, 63.55380257386972, 43.919357458252804, 23.01024424148823, 52.93824705598195, 75.43661023465496, 9.888737247731447, 10.475180869403623, 12.89789137394039, 86.39553385485206, 96.41423892357255, 46.41316530770811, 37.556570428969856, 49.45516408339277, 6.617917066333705, 54.147849956046876, 8.084325092096178, 35.802227558745834, 62.76021087535632, 56.88338260245663, 130.6711572084285, 40.56750600212464, 25.49335405251625, 6.576829047498817, 52.736603576252236, 61.237742416801304, 48.20686758601432, 6.074090595487571, 15.940573885565586, 7.151321960700783, 16.2914429290514, 54.63575849855518, 44.69204191866788, 32.98275736333001, 25.995671096857322, 14.31349780202911, 25.033433646329264, 12.407688834811019, 6.197103995513198, 11.698890112459546, 11.242860867452213, 61.46559166622376, 66.70496444239848, 29.360754381563762, 14.196409422614167, 173.43949869419188, 205.21444123748222, 17.427440157357893, 124.60847524303543, 42.2119156838969, 9.487726679291303, 32.12754485945547, 51.243098127498214, 52.13001059678446, 14.367510538450862, 21.861735304659298, 13.366603473479222, 54.168829738192784, 72.72168227726169, 70.9424489442474, 89.08706179555452, 7.373313596795591, 135.87841700917642, 13.563428602116813, 52.294007630488146, 6.03574452193149, 32.46112106928729, 80.24642603544508, 41.77935585724417, 62.44800718780023, 85.9149071529233, 95.6962602770132, 25.3958146126336, 9.296729245050015, 7.271880166555046, 89.8361150902715, 17.639128628348196, 127.51661966346144, 61.41548376518578, 91.77511047191413, 30.09324106604355, 17.436145515896865, 30.82023086744954, 5.562024721678407, 70.33998509823505, 14.693376565335392, 23.424056867223538, 27.749510271548267, 52.812905721240064, 30.3098036361471, 99.51020235291591, 5.540517540550961, 31.772223427695096, 45.8965926027949, 7.374882944576174, 6.086800553258134, 15.785608154027267, 8.749506438162307, 7.301604934364692, 77.32987370266302, 23.555715244751187, 21.8786219800095, 8.387161530518998, 38.02464967434318, 80.2738366049172, 11.759039832530886, 23.086243794466412, 5.771366899993589, 39.85785672538417, 39.12722844289006, 92.91403045781837, 156.85731100439705, 94.66225747566423, 7.429430069475531, 67.44153409041668, 74.59239359829178, 114.4623992148571, 10.051727391426644, 28.113319219764335, 19.331010790221494, 16.11048730103169, 51.794892973848604, 41.95453607240488, 54.95996210986396, 7.5252960554845805, 35.53186933286451, 99.96477038609784, 9.610001022015487, 51.56846671376123, 28.982923462512446, 86.05114115112606, 18.51965362359951, 25.195128065048458, 53.13692907005788, 62.34127565818207, 59.715944529763384, 9.25171666538587, 51.6489905775533, 46.17152996801633, 18.453163735074376, 25.287592676196066, 43.787470420982025, 83.66012110656628, 46.6402921641026, 19.34364728420873, 12.090490520501456, 6.240278905413386, 19.012242862990725, 6.416142371763144, 30.446169598231187, 14.139277350873845, 27.978337069651904, 71.59450683893728, 66.36487197355532, 116.90839376226157, 59.17809813705439, 29.775401455902266, 70.90357127638957, 40.04808196076903, 7.611867593801878, 98.40674967720241, 26.92596270771244, 14.497536487310022, 16.02166294164095, 144.06385216261046, 48.62399427477424, 73.47949365020442, 8.594823383759417, 74.05095983615143, 13.546102842927178, 46.63661713175448, 69.56341089096935, 21.91789367967564, 11.735111628696554, 22.142625269171518, 32.002270048779195, 42.204814344109835, 22.899857877829845, 19.031111692547118, 59.49666932087761, 44.627803443600456, 67.05204489559316, 17.508628846794508, 30.39458265216645, 43.345792690808075, 21.022931173373195, 29.070446501135432, 29.99797570216571, 52.102735893667, 6.377865142572682, 5.266000795281109, 59.62495941644705, 27.277400204586463, 206.84018002133402, 45.02901706739867, 68.19244353274112, 42.37770458631099, 78.30816841088722, 185.95939386039947, 70.62829096274824, 115.16122749467726, 40.218162398003315, 32.60279394256348, 16.540756214495634, 55.50737314910614, 39.390508513644036, 27.728811880335755, 95.15657115800794, 12.956759891372322, 71.7314885580177, 38.210511165429956, 5.506273196787376, 50.67208622977007, 88.07391667117014, 15.193828212766471, 112.62899174153308, 37.5015304182949, 22.018191016381163, 68.47414725799158, 28.32038456990839, 7.016172240331582, 28.445266863418382, 47.20740459649353, 21.185522844303225, 6.967326293996054, 42.11691045642713, 26.45409229157432, 155.27544494740954, 5.177086943933869, 87.11463274885384, 98.15602438969994, 41.45646777609386, 40.82943453484009, 6.134975346618311, 110.64372977868027, 62.2339636612124, 21.353354503218227, 17.690927850082105, 41.1880027080968, 18.59438752306608, 5.251245528259005, 80.01210067352088, 36.00276588777652, 18.81077143953084, 55.118472665717626, 30.039363902200748, 74.22049446597399, 46.90936199130296, 9.908559986058348, 17.11277845915571, 109.16717911069404, 37.37516937285763, 35.51583663351791, 62.192104250686526, 6.722930732293959, 41.09043166829201, 107.44944307224371, 76.4716792331794, 13.31265237845868, 26.0619587964073, 14.60357168296921, 15.270581374537088, 5.547779206414565, 74.15917302598947, 17.34287689319821, 5.524041923611538, 91.61268523532743, 82.23288021946027, 28.66258016231189, 25.083309343602085, 18.305894510912008, 11.865352898185584, 29.852207572898763, 22.33450854734035, 13.257051405218004, 49.726473677150686, 54.601991080018976, 71.2597066597719, 84.01063941639904, 35.49505355652587, 93.79843519859801, 15.68269031646279, 44.26204534874868, 22.67434708861132, 112.24608451983018, 13.436237123912555, 6.353399137402942, 21.017091219679685, 114.37584594732776, 41.89707069831421, 6.942584154398579, 16.50850163933523, 18.57443256995602, 68.77056582430806, 98.42990465030982, 52.323462314432355, 6.284595118619792, 135.61081354767194, 53.503255212290355, 32.55315815662022, 103.9825048778239, 14.441371419713642, 149.4500938009838, 60.55300531501252, 159.73576478759458, 20.462234155829286, 85.40999915027011, 44.830266160964996, 69.8709751443631, 12.600695424099838, 49.70777220682642, 31.134932741801837, 44.67853744047713, 12.913321184739036, 44.410528568998515, 80.89785451680575, 124.53323675450727, 41.134293649110944, 58.775219974976, 101.97560395560843, 25.196936930288974, 52.84964313625933, 53.20986339075504, 23.607944156480748, 40.20613385199743, 7.09053907774558, 9.161611748811067, 100.40686317953406, 76.8505708126323, 10.070680414025194, 38.680465525076265, 26.07296711513674, 95.11151012711966, 8.357730495497414, 36.605707419413264, 63.8598890055455, 12.377090374751722, 66.75748662706995, 14.968322651239426, 27.804981391383635, 22.10676060245977, 8.993932220826862, 36.9085924589413, 5.6377866380751565, 113.43088218598271, 84.45790557888111, 22.394274274559887, 58.54920632364563, 22.71435657008975, 44.8913671637996, 13.911799463428904, 13.728281485434376, 19.150535145545902, 51.377452720530435, 36.68975611028691, 30.645004890194983, 9.455885651841506, 15.233478004502677, 55.631317699723354, 49.300862500962175, 161.59919467714377, 24.980748227581522, 9.726103225243985, 18.370513047381255, 37.172800474716496, 11.421991218353877, 38.323132887689454, 41.760676636721925, 9.323056115934662, 10.310466763641115, 78.72991463916759, 24.914633917131162, 8.814065689379412, 78.53581305810236, 14.06509920815876, 41.93971447769981, 69.69164099130191, 11.477325140736363, 26.75752412790964, 10.215085828001309, 69.89188717479753, 113.18339962873547, 130.03007456240235, 27.260158710473284, 14.253568695148076, 94.48413618988229, 105.55620901609637, 29.854065329518868, 102.95286200566437, 31.10033772785273, 6.849499239611752, 16.09435204052012, 34.263079763037496, 152.05090160172193, 32.91372536810881, 10.795608878117786, 21.11439496595769, 26.57605755229211, 234.4321639864748, 18.77366329613951, 77.8789585018097, 40.14073192674408, 12.15367125819885, 45.365951030545986, 39.42761410472773, 56.64166862771388, 7.707224747318835, 45.10793019880295, 24.015195898340927, 54.35573414782981, 7.57506984591888, 8.84424750563413, 7.080166757187136, 76.26817418251031, 6.238563157711726, 9.139766552876395, 24.02479033910739, 54.97474384929673, 8.750830569022959, 5.02608325252503, 91.51394733406534, 6.7700506558429225, 56.87981644273842, 29.39628482227317, 28.831415817558295, 71.99769897983538, 24.221276936329392, 6.435490489730579, 10.128202661889697, 12.87691591399845, 42.77654693933018, 119.9671510729133, 30.279461042153734, 41.45253137717987, 20.94816474761589, 12.201378308452913, 11.18621682671237, 7.266155582998123, 5.445175995065317, 14.214353619881635, 8.599226627305159, 39.8440815802911, 42.757769442579985, 24.652540503521674, 48.54345280295502, 22.575182363354173, 90.62945528068869, 67.0776418832133, 39.01717021678313, 41.456834447734494, 83.48234554651404, 14.647882435466503, 33.68568462216086, 147.80540979695195, 6.328797527910687, 47.703484041433114, 144.19888872822813, 8.094708284892777, 33.29828771044275, 57.35404278575308, 33.92769941592709, 13.471009960733731, 74.09015515088883, 5.694252196734011, 11.69505950570191, 24.963533392118325, 69.2989177756227, 107.10465511599463, 8.18295246096663, 26.920217970702648, 80.408428396774, 13.038334115335761, 61.124891169077564, 64.83192814811643, 9.351951560310564, 18.65260142067794, 12.887190726114909, 7.526394040069647, 117.8542058627798, 6.7049697792421235, 5.305480799182547, 15.511858167455491, 158.003421471843, 15.14468356362237, 100.3086861773746, 15.256678399609854, 6.160474786050645, 9.553610452182832, 48.341340959699366, 54.47442227717069, 66.6248936672985, 10.016615835763426, 18.80126186023898, 51.17998097565463, 14.301679867052636, 77.52751029771365, 8.921892657072155, 54.89846921988301, 17.33076715646697, 102.10316010293691, 38.450140532371854, 128.13530916122204, 59.607154563478474, 8.44095068429451, 24.566795988847897, 10.346857269145369, 10.982188901561608, 84.31591920415146, 27.25081613684647, 5.092496767658493, 45.78943359716188, 45.102119009829224, 18.9284950636643, 51.104402998203305, 111.79485160976334, 49.53916811042245, 81.03624590685057, 82.52736583402798, 69.3828527196381, 38.15212096690889, 78.56843966412396, 57.75017315547217, 66.56362324616298, 19.915158544044868, 6.7761837034514665, 12.46642960681854, 26.503109069968744, 23.4609830696894, 200.6357994925697, 146.34723781013432, 26.060271222565984, 82.80220866999161, 27.889096638261737, 69.91641671835687, 187.24765842954773, 14.764638400887188, 23.873318930481776, 10.311096077659993, 62.030587279356716, 11.972088670592324, 12.730767367373948, 14.744926948739588, 75.13660392289218, 16.51979112769956, 45.92522289750906, 26.314558628522253, 103.6013029108133, 110.19349901986118, 128.84037061445736, 15.108716975852179, 107.45754482563311, 13.343119011495226, 23.62001218040819, 35.52648172839175, 65.78981482359202, 54.23635156563837, 19.510966824017967, 54.86017367432424, 15.63093797440105, 285.41646802802813, 45.69000009228574, 15.331499618514034, 52.43298035522011, 34.27608752401803, 15.040043169067104, 6.919990493655217, 71.6797750061305, 24.632552222249586, 31.53428674693425, 69.19490509983791, 94.91574108996396, 8.449035442929048, 13.272672253935308, 180.03725244602123, 38.36112770745725, 15.80421830495635, 12.718044972855013, 57.66727139151401, 124.32413496646772, 5.890872749530303, 141.78366850821695, 49.21254345002892, 24.862206126987317, 127.11754308219128, 108.01070954329697, 191.72434928058928, 72.02571590310596, 29.860304521020527, 5.574188017950052, 33.00205110082356, 17.737171289568295, 34.43730538522372, 5.111448490476592, 87.3112673130835, 81.14294925463486, 25.06552323976976, 30.416920697069553, 24.81200644198517, 110.79433240950424, 5.525648997796797, 82.54634276357329, 34.207482291239515, 7.2244942474555, 50.393069818758974, 289.9863771281537, 9.294192209314993, 59.506221368054995, 148.01703662683747, 14.429850073122667, 7.1711069596050905, 129.74907338994228, 20.36961735410386, 8.920454949366103, 5.552508556854841, 6.06511613050196, 29.70379490382412, 6.641138105415583, 59.772044043957976, 32.0481623617099, 60.902347472461955, 95.8247708695723, 66.64354265054752, 22.42525927588848, 83.02261752987962, 76.47169997903089, 18.248775239295398, 26.681639978982467, 84.20116422305013, 149.20322021001482, 45.67598245106359, 16.11656586360676, 109.41111175694266, 150.36828264807724, 6.419062870226426, 32.98809616617845, 40.07854015769477, 133.3801368369283, 13.186270935428443, 25.758832274902964, 262.3294823261272, 56.257633196179754, 6.2890594725513616, 58.52036261183278, 95.66214036585107, 118.29535592895849, 56.59520466918407, 9.26753851237654, 17.523904749402295, 21.409243826750533, 12.689661624799005, 111.58000988517797, 67.07812222020456, 5.986382592171176, 120.56743583247504, 30.441937892191945, 149.67243871328762, 27.94921410683981, 39.75178348026374, 6.736505556424783, 197.27161225840402, 14.063270527564626, 6.07955731649008, 24.509736987045585, 34.843079116447335, 105.9847749437468, 36.87142771067538, 9.002835603148894, 77.89318863267114, 91.03190717933538, 55.24562061211194, 70.16081829807507, 39.834999718723076, 35.282617266385046, 22.932858627925537, 13.083635004557983, 36.02916338132999, 24.3613847957462, 20.7625135836459, 111.11582473938283, 74.89439280824058, 157.00914179930862, 79.92202177784577, 16.442159872981588, 73.87337579656358, 62.14438589694732, 59.148395646215235, 68.094066942419, 38.53874835114903, 74.70283607632217, 47.317640502741654, 5.778896196151821, 172.2530912279364, 69.49620086919462, 26.001595138787152, 24.709647298132136, 52.415221573384656, 69.41312245627446, 76.84416934341235, 23.627392125530026, 125.59241210530469, 130.8182542375304, 5.973048415159237, 35.23535088861591, 161.4597549767803, 40.03074657455956, 142.01891556901927, 26.863609222982504, 44.89134253814364, 58.96898944948762, 113.93732092138178, 9.160184686044985, 75.63429296358292, 149.2098021182579, 85.2168357443131, 256.4584291439799, 10.999713587102859, 22.915482147833956, 7.725605761112345, 79.79078745708637, 41.31199146266421, 5.04241315796255, 52.74131554358093, 32.29531434318677, 10.690523101746015, 57.08718483670773, 48.972701567198854, 11.443724926437923, 5.9941398408934905, 218.0114386269908, 10.760779063747075, 6.794315663464759, 15.202063307290153, 108.53077834817229, 120.22234106455603, 6.780267788600447, 26.185225598766877, 26.75143978064268, 7.9583949057945675, 17.24768056533221, 33.954093994392, 47.04282213920154, 27.498754487763403, 41.33288818306591, 20.6777671441621, 59.0258496188467, 43.87522345295311, 61.545955207928706, 42.87074019796485, 43.07556484745181, 10.566121665295805, 26.321645473810776, 6.303850598237803, 9.712834439690447, 80.67118243646796, 52.95424603686472, 6.148352706654838, 53.67906388039738, 60.92083116930623, 75.68472294765404, 113.31443379749983, 155.35045909199408, 36.821991516969675, 160.8859167163571, 56.264180861206974, 78.00386262070087, 30.90678716430718, 12.780285019615246, 16.233017842299414, 15.590668832209131, 29.917501354037608, 37.323004688903474, 25.80216880840069, 30.995986736378207, 39.21098722866333, 30.919047224710905, 22.727814482073246, 10.967551429868204, 109.9834275188002, 83.32517340778499, 55.833565750028384, 15.68302399677503, 43.73613658005283, 52.28963701482757, 48.21845633612445, 5.27478532295761, 16.49313810919062, 37.4853504266769, 86.83622237589096, 26.88375651187652, 43.64692058661242, 10.703590272493596, 19.51486118811333, 11.510975767421295, 15.340383725070527, 35.44902153180603, 57.20358509925974])
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)