t3pa2cls

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 = 44859
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}'

Functions

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 

File names

#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'

Calling functions

#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)

Hit map

raw data

#hit map
hit_map(detector,hit_map_fig,RowNo,ClmNo)

Energy and time calculation

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)

Detected energies during the discharge

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);

([5578977.992637578,
  5613539.99644256,
  5621984.067356959,
  5625185.491974772,
  5630368.1422281265,
  5632342.76281048,
  5700317.071617094,
  5743468.085524118,
  5744836.701108547,
  5774398.334295918,
  5774767.373556603,
  5774883.6945538055,
  5777203.405797355,
  5793085.617795732,
  5821571.743378873,
  5855827.613748997,
  5925625.695451189,
  5933944.544405776,
  5934801.5625,
  5940192.840524625,
  5946482.516304447,
  5969773.381161929,
  5981136.763201982,
  5984070.243328004,
  6013602.588212368,
  6013635.187115682,
  6031473.283657751,
  6037470.24548906,
  6042843.5027544135,
  6064294.474320638,
  6094475.0,
  6094905.9829924675,
  6095592.947819843,
  6095976.71663147,
  6096367.1875,
  6100604.671760609,
  6127496.660020787,
  6127514.66433957,
  6131121.586458964,
  6135403.121850532,
  6136735.932735263,
  6141650.962756152,
  6142076.310749517,
  6145591.865585281,
  6145618.491309015,
  6149554.421623171,
  6150099.041201457,
  6151184.375,
  6153137.020367047,
  6155257.734810135,
  6159103.125,
  6159982.962820449,
  6164579.583200394,
  6165053.125,
  6165408.029470348,
  6167709.375,
  6175208.476366167,
  6182233.18007594,
  6184032.704570248,
  6185279.554791648,
  6185349.042431714,
  6187329.873797793,
  6187452.158667843,
  6191708.011044823,
  6192119.60121315,
  6192157.165202618,
  6199457.507883905,
  6218398.34190127,
  6218690.570040212,
  6218746.875,
  6220405.099065951,
  6220965.103675448,
  6222523.60000546,
  6223349.563405839,
  6224387.228177192,
  6224443.575621733,
  6226151.5625,
  6228858.775595083,
  6230673.432989561,
  6253256.072526208,
  6265052.452044818,
  6274830.950682911,
  6277836.20900789,
  6314249.280928683,
  6315299.579220695,
  6315736.783741842,
  6316443.123659642,
  6316910.9375,
  6317419.647893409,
  6318677.500032619,
  6319442.44569348,
  6319450.113932871,
  6319712.830529361,
  6320429.6875,
  6322545.811595063,
  6323816.559097993,
  6323824.821548586,
  6331484.375,
  6355723.014527185,
  6357360.588598309,
  6363796.557840736,
  6395125.521991844,
  6395641.637876963,
  6400400.543865813,
  6401275.067899309,
  6423478.595334258,
  6426760.9375,
  6434744.811132539,
  6442036.355002249,
  6444163.76806014,
  6452397.865784013,
  6462296.591062106,
  6640179.5698098745,
  6924175.0,
  6932273.239346904,
  6978174.583763218,
  6985594.689982541,
  6986289.0625,
  7009900.567918236,
  7049085.756867959,
  7050059.154200402,
  7053084.837842204,
  7054303.048294704,
  7063516.622922211,
  7073588.617017695,
  7200593.640329458,
  7204766.782160621,
  7240366.622016894,
  7317253.268836167,
  7364354.6875,
  7365233.114811757,
  7367649.523492468,
  7493069.116851672,
  7503121.868162566,
  7503184.560158351,
  7531570.1224347735,
  7573246.109365999,
  7614677.878747302,
  7642043.848840267,
  7642056.026815211,
  7740341.077747443,
  7837390.625,
  7871174.495460565,
  7872129.6875,
  7908156.925440548,
  7918716.881316441,
  7919718.632875429,
  7922573.709157177,
  7923678.177498328,
  7924698.879963058,
  7925407.8125,
  7926249.4093552325,
  7927146.96573762,
  7931012.345965633,
  7939362.115827618,
  7941796.346337401,
  7942101.001009151,
  7943187.988806127,
  7946514.697059679,
  7948926.203486298,
  7949316.016096984,
  7977941.478919659,
  7997917.289809064,
  7998522.693319561,
  8001081.420265849,
  8004648.631118895,
  8004875.722019136,
  8005660.933549159,
  8006579.534724715,
  8006725.0,
  8006839.0625,
  8006954.982699717,
  8008182.244234026,
  8010172.892383554,
  8025226.5625,
  8027709.559995733,
  8027712.273284319,
  8028437.681708809,
  8028718.75,
  8028848.4375,
  8029351.5625,
  8030138.62628424,
  8030164.0625,
  8030282.884151183,
  8030343.75,
  8030579.296380469,
  8030788.566129072,
  8031312.5,
  8032257.8125,
  8039522.919105169,
  8045755.033987068,
  8045815.030558283,
  8051632.177019472,
  8055023.17636913,
  8055675.966119183,
  8055677.96209952,
  8058149.7678323435,
  8058416.82410652,
  8059667.893819579,
  8059670.3125,
  8063958.710953828,
  8072553.726499772,
  8094066.422983345,
  8129896.6949967425,
  8141425.052780143,
  8142815.625,
  8158021.623618728,
  8167913.766938603,
  8174902.733161797,
  8181806.800345287,
  8187138.001237432,
  8189534.981143247,
  8189853.418130742,
  8190473.4375,
  8190840.761745352,
  8257782.8125,
  8262367.272545,
  8263788.017466279,
  8268400.745466446,
  8268609.375,
  8268650.712523266,
  8271246.875,
  8274221.875,
  8274650.0,
  8274695.3125,
  8274698.026401529,
  8274734.146565861,
  8274752.336769589,
  8274771.139793693,
  8274815.314464465,
  8274909.577474334,
  8276568.75,
  8277005.824703687,
  8278117.1875,
  8292855.987967529,
  8329567.989503332,
  8369545.9785145745,
  8370513.431092635,
  8371887.753969119,
  8374074.405507117,
  8374076.643530111,
  8382171.457617436,
  8422064.237096388,
  8422096.821402773,
  8448323.72457198,
  8448454.343248893,
  8452726.160929244,
  8454011.583374472,
  8454603.548884489,
  8454943.75,
  8455778.728842134,
  8482536.891196294,
  8482783.0513977,
  8484149.487759253,
  8485708.153880486,
  8487173.382844176,
  8488025.0,
  8488662.5,
  8495114.019980155,
  8495364.513272393,
  8495590.435909467,
  8496115.522262499,
  8502990.527542615,
  8505265.625,
  8506783.675371563,
  8506812.495052168,
  8509096.150287287,
  8510243.75,
  8510427.786932219,
  8511181.25,
  8512575.838905351,
  8513537.51068289,
  8518229.666317971,
  8518587.5,
  8518845.066480987,
  8538177.352235738,
  8538651.698696647,
  8541294.02065175,
  8550751.333775036,
  8550788.678866152,
  8552459.256242624,
  8552723.455720842,
  8553490.139605707,
  8563675.928755844,
  8568046.706448179,
  8569368.75,
  8570442.8937564,
  8572202.725852901,
  8576441.562753767,
  8578176.267249139,
  8579606.241359796,
  8580919.879156416,
  8583482.621314054,
  8585365.625,
  8585791.663102852,
  8585818.44421204,
  8586276.206291188,
  8586568.75,
  8589244.446631467,
  8589526.5625,
  8591485.654464064,
  8600761.594623223,
  8606515.061278822,
  8610595.910508148,
  8622675.934620805,
  8625634.908882054,
  8625814.485315993,
  8626881.52629706,
  8633973.971662201,
  8639836.424902648,
  8664127.054110654,
  8664489.716332259,
  8664549.756569194,
  8664722.522404205,
  8665114.948647423,
  8665252.087675111,
  8665749.782820417,
  8666228.125,
  8666942.1875,
  8667709.605118735,
  8668166.404484307,
  8668435.85852967,
  8669265.625,
  8669269.705527764,
  8669285.9375,
  8670121.88575467,
  8670287.5,
  8670294.517241456,
  8670517.89229598,
  8671609.375,
  8672554.104710143,
  8672950.339512775,
  8673110.460658483,
  8673431.25,
  8674107.46158853,
  8674760.951248473,
  8675641.758034034,
  8677970.237481065,
  8678623.8448693,
  8683148.323229674,
  8691526.043281795,
  8707683.95163068,
  8713364.0625,
  8725007.569254583,
  8725190.625,
  8725530.76560904,
  8725793.054255715,
  8727321.962750567,
  8728206.780048713,
  8735364.953771826,
  8737135.380865166,
  8741326.743886946,
  8752984.99766059,
  8754162.15317943,
  8761408.526926685,
  8792685.360675436,
  8792758.041077849,
  8793653.832615513,
  8793996.875,
  8794022.188779937,
  8794154.6875,
  8795154.6875,
  8795201.5625,
  8797156.25,
  8798655.251237357,
  8802349.879995534,
  8802656.63022709,
  8802989.522512034,
  8803828.891342113,
  8804015.738557713,
  8804023.4375,
  8804145.943872036,
  8804186.860842645,
  8804192.780479454,
  8804219.813627753,
  8804569.712109134,
  8806134.96419869,
  8806203.889612991,
  8806254.6875,
  8806975.906851102,
  8807864.134671347,
  8808102.117476793,
  8919282.607320037,
  8927556.573354574,
  8953757.731735593,
  8954361.746108068,
  8954362.31208534,
  8954415.39249643,
  8955966.60381047,
  8958651.5625,
  8962332.555338075,
  8971112.321127882,
  8977868.233932119,
  8980165.436467165,
  8983996.92816116,
  8984599.17041027,
  8985237.5,
  8985248.935334267,
  8985299.896004494,
  8988940.625,
  8989925.0,
  8989927.606602695,
  8991127.35511087,
  8991185.500495778,
  8993061.513141083,
  8996273.120763056,
  9004222.812932868,
  9007014.95784737,
  9028259.238826694,
  9030218.242034517,
  9032344.708575925,
  9035695.65798077,
  9036945.3125,
  9037044.923101418,
  9037098.69224798,
  9038609.375,
  9039574.616101801,
  9039973.000686746,
  9040032.284541702,
  9041434.707256664,
  9041732.359365825,
  9042234.361240216,
  9042501.5625,
  9042831.284776246,
  9042887.5,
  9043064.575609773,
  9043202.546242805,
  9043281.25,
  9043592.066514777,
  9044362.991936032,
  9044486.512079779,
  9044854.017704967,
  9045270.3125,
  9045492.95490584,
  9045494.87490307,
  9045720.127382657,
  9045908.082004739,
  9046732.13088164,
  9046779.482426535,
  9048838.858840046,
  9048903.789355332,
  9049677.784857232,
  9052394.930674184,
  9053596.875,
  9053828.930609608,
  9062087.914366238,
  9067817.172224868,
  9070774.077437205,
  9107027.744620299,
  9109396.695817403,
  9129387.419966767,
  9129628.125,
  9133234.420686716,
  9137449.976875912,
  9157097.79604507,
  9174989.89398588,
  9178357.8125,
  9180308.896938387,
  9205695.091304233,
  9227506.414677413,
  9233628.954627037,
  9279602.294272551,
  9289142.182902094,
  9289290.504566358,
  9299359.361218933,
  9303648.4375,
  9307773.782939412,
  9309326.5625,
  9309650.800261535,
  9349581.598764736,
  9351163.369059006,
  9358130.934509596,
  9360859.375,
  9361247.1715868,
  9365434.482828114,
  9386427.518439142,
  9387284.375,
  9387340.625,
  9387694.203454724,
  9387923.469339475,
  9388744.119543774,
  9389051.614938,
  9389479.370882215,
  9389733.16428069,
  9390684.004603539,
  9391330.75013092,
  9411608.573555443,
  9412030.666474579,
  9412792.944299256,
  9412845.3125,
  9412851.308301268,
  9413583.376862906,
  9415168.94505935,
  9415216.51633216,
  9415425.442093844,
  9415459.87429085,
  9415565.631651614,
  9415673.838264123,
  9415698.997529536,
  9415723.4375,
  9415790.500855843,
  9415873.614106368,
  9416124.885779316,
  9416262.24589332,
  9416313.714787347,
  9416316.916951617,
  9416423.613059886,
  9416474.545621373,
  9416500.829499535,
  9416602.546723524,
  9416828.932083739,
  9417016.571739135,
  9417062.328590088,
  9418001.5625,
  9418397.90640377,
  9418610.835803391,
  9418635.250859113,
  9418748.08130043,
  9432352.294239253,
  9550716.457110025,
  9552684.303133195,
  9553108.909647174,
  9654869.71805141,
  9745972.988154227,
  9748710.47078961,
  9764723.150676485,
  9765188.735763712,
  9791240.625,
  9798694.121990992,
  9799732.710472522,
  9800043.708696079,
  9802784.171198906,
  9802805.990969595,
  9803242.267307982,
  9850576.5625,
  9850640.550180549,
  9872050.018251656,
  9905679.190650804,
  9909968.176003896,
  9911297.66728433,
  9913121.840636093,
  9916705.572066901,
  9948380.833003595,
  9965781.261695607,
  9972839.846731849,
  9973690.479978068,
  9973754.6875,
  9974052.746241923,
  9974161.851776069,
  9974490.047096923,
  9974545.950481862,
  9974584.68219173,
  9974604.730765797,
  9974646.83061326,
  9974682.513091875,
  9975133.547908572,
  9975212.208521167,
  9975220.735421192,
  9975229.476164848,
  9975235.811844306,
  9975765.021756265,
  9975856.501909183,
  9976007.419921272,
  9976057.80337224,
  9976065.266040968,
  9976886.342057634,
  9977187.243001143,
  9977210.144756453,
  9977274.089917224,
  9977470.3125,
  9977504.185187697,
  9978905.05714484,
  9979243.707714986,
  9982509.31431189,
  9985512.5,
  9992113.654304542,
  9993020.555281322,
  9996762.471096417,
  9996787.381402094,
  10002379.336804705,
  10024478.125,
  10027391.971597543,
  10035632.623142052,
  10037956.857204715,
  10038008.361725478,
  10040890.37532692,
  10041982.8125,
  10065853.38125877,
  10091403.838716848,
  10092895.3125,
  10092963.321090968,
  10095299.720654031,
  10097807.304352304,
  10097856.253012737,
  10098667.1875,
  10099229.6875,
  10099909.014946127,
  10100182.936398685,
  10100247.188947475,
  10100470.178203166,
  10101104.057642618,
  10101162.825340625,
  10101854.6875,
  10105054.78391054,
  10131292.45049253,
  10161616.416436505,
  10161626.269829186,
  10161729.834783094,
  10169688.595853087,
  10170138.014750518,
  10179129.221269872,
  10179714.788008817,
  10184406.21191207,
  10188168.196301162,
  10188721.875,
  10190524.772325411,
  10191598.302842084,
  10191712.078098169,
  10199184.375,
  10199235.857941551,
  10208423.172623368,
  10208909.375,
  10210021.875,
  10211332.8125,
  10226490.625,
  10232635.438595347,
  10235212.585390566,
  10237767.544899004,
  10270728.125,
  10283002.062005082,
  10290050.0,
  10292251.5625,
  10311263.88779654,
  10386710.098369364,
  10391426.332413316,
  10466773.854691127,
  10477916.968040938,
  10480777.584047493,
  10490746.993626973,
  10503769.415995158,
  10520482.016492838,
  10541145.250601705,
  10541158.407772304,
  10548296.875,
  10573725.0,
  10773496.875,
  10774728.905083535,
  10774751.5625,
  10776792.206136825,
  10777665.020240204,
  10777668.13461999,
  10778846.71205643,
  10780182.317867313,
  10780895.3125,
  10781879.34564069,
  10782364.0625,
  10782616.685812997,
  10782751.512144936,
  10784220.47316309,
  10785047.461768515,
  10785567.649202945,
  10785963.83995435,
  10786879.554680398,
  10789251.242923414,
  10790342.008619381,
  10791541.585273359,
  10794124.449462075,
  10794867.909633797,
  10794883.177574959,
  10794896.152109306,
  10798670.88748433,
  10798675.251460712,
  10800201.774210311,
  10800576.472235681,
  10800866.467850756,
  10801371.854027035,
  10801583.580403736,
  10801655.203027135,
  10804001.787265554,
  10804373.7500902,
  10804408.038444053,
  10804666.718057577,
  10805223.262449266,
  10805324.392792817,
  10806161.656789578,
  10806224.997911619,
  10807836.321072653,
  10808036.644279381,
  10808415.625,
  10808581.873706874,
  10808873.4375,
  10808939.833839264,
  10808980.868800774,
  10809117.556962334,
  10809218.567149037,
  10809468.243870266,
  10809663.455287416,
  10809883.819216076,
  10810860.650136447,
  10811066.631674938,
  10811229.246688634,
  10811451.55775868,
  10812146.451787293,
  10812288.74001905,
  10814082.73614174,
  10814429.6875,
  10814470.794442,
  10814634.175857754,
  10814822.984567305,
  10815081.21234314,
  10815575.934037207,
  10816437.5,
  10817226.336476753,
  10817231.25,
  10817241.190187208,
  10817316.994471196,
  10817685.920652514,
  10817708.61376789,
  10818650.44962985,
  10818785.0410984,
  10818949.73800539,
  10819542.563609054,
  10819618.57517276,
  10820283.905993085,
  10821180.385645043,
  10821563.474445011,
  10821962.366196258,
  10822690.405543769,
  10823060.742860926,
  10824867.355729358,
  10825224.736878874,
  10825530.41632611,
  10825575.0,
  10825796.616920745,
  10825852.074789431,
  10826027.925790953,
  10826181.140819779,
  10826309.928126782,
  10826581.25,
  10826678.736373166,
  10827726.356233587,
  10828785.467177365,
  10828997.179867625,
  10829662.646082487,
  10829965.531124143,
  10830422.094331188,
  10830576.343721554,
  10831079.266695555,
  10831707.462456234,
  10831860.775297178,
  10832075.0,
  10832606.842710985,
  10833504.725712355,
  10834208.209129078,
  10834461.755489081,
  10835198.16407682,
  10835381.25,
  10835405.036137272,
  10835788.706193324,
  10836144.412098594,
  10836521.80788944,
  10836907.84599894,
  10836979.299859766,
  10837014.0625,
  10837414.803512307,
  10838227.097552575,
  10838352.524900647,
  10838645.084898634,
  10838932.537461119,
  10839782.707676722,
  10839966.237466928,
  10840163.452415587,
  10840168.75,
  10840573.001947701,
  10840835.183819659,
  10841051.5625,
  10841305.140630526,
  10841315.26370543,
  10841873.732887015,
  10841873.955943927,
  10842635.052716026,
  10843007.96124787,
  10845167.493300281,
  10846154.474631691,
  10846417.784582015,
  10846660.9375,
  10850361.280366791,
  10851604.490230924,
  10851920.3125,
  10852854.71921794,
  10852988.403083177,
  10854441.81241589,
  10854476.873962175,
  10854951.012120415,
  10855042.432351183,
  10855117.968278056,
  10856024.359284345,
  10856636.050576696,
  10856885.9375,
  10857321.755933389,
  10858837.696414525,
  10859606.116575606,
  445897168.75],
 [10.926932224398904,
  20.299863187360607,
  128.03679543984842,
  22.054290030078057,
  92.0134442811754,
  8.616197603256872,
  6.076130750777888,
  99.513277837992,
  62.942603518181855,
  14.562212679960686,
  69.58291917703299,
  7.9567872446315615,
  6.820890981208564,
  67.7001006400315,
  6.73751087261893,
  117.25589082591692,
  186.77098044642872,
  15.759888942893621,
  111.42139242152703,
  209.23843133454497,
  8.135918106382611,
  72.9941144126226,
  25.883289989482332,
  36.68592379563776,
  64.15916495810745,
  53.54018954297891,
  10.895464639696927,
  26.62252597598342,
  72.73778335531173,
  8.7107392510694,
  49.865399405062874,
  112.20082055767892,
  82.64480473827645,
  6.286870166924756,
  48.66363464808074,
  17.15804067148037,
  9.156008998167568,
  48.4507440634721,
  59.66783805117294,
  28.462945030734435,
  28.43112284737409,
  13.060927451972757,
  24.3872693557942,
  42.247653860659256,
  6.061604517463069,
  24.196628613696685,
  12.223600475285949,
  87.33087495537146,
  112.04535780831122,
  99.35487052097429,
  34.3612146902758,
  69.82759419328755,
  93.7767087179709,
  80.87655749057232,
  65.84386632806763,
  46.10461411743419,
  95.29772568023517,
  54.90189286240101,
  90.08621857118047,
  26.15196044445469,
  5.422114749261968,
  20.46881669574627,
  9.892079739524835,
  64.04483853197218,
  19.616002262174746,
  22.426241203402252,
  23.69195598151165,
  26.767377412322464,
  27.48242570625202,
  78.17984072236509,
  16.782526933300517,
  5.776388073977832,
  17.16245783712184,
  22.14981813755204,
  41.21862258940446,
  32.30695594607249,
  52.61387621741336,
  16.04861373266941,
  55.8379853597047,
  25.453079411126218,
  78.58967712147752,
  23.759468818820626,
  6.116796188298128,
  59.45093622154312,
  6.682311490039749,
  55.73381667475616,
  65.93507541908905,
  74.76998507261524,
  7.963836959850351,
  12.996342339540242,
  115.28563924400606,
  11.310586744054634,
  43.415073325361774,
  89.11735884359675,
  7.794309767515815,
  8.851141962612628,
  14.302330732534013,
  36.26773275395881,
  130.4621411981585,
  23.145561783174436,
  139.05611076736474,
  17.39055387202637,
  13.192164245260043,
  21.52157689636502,
  44.57887992328304,
  137.21496266899408,
  29.43156867413411,
  21.479469924041624,
  11.849662875070816,
  78.45085521222896,
  78.00925312492446,
  64.21466078080728,
  25.522473652892575,
  44.0564159458759,
  7.070319545796551,
  22.37146180037758,
  21.563094880394516,
  107.55994206406913,
  145.031435332794,
  14.350142723374486,
  61.36056062936864,
  97.6706598756864,
  36.30821344342823,
  5.229613973701861,
  95.05866320511922,
  10.960505086983037,
  10.54819837019146,
  10.34560973780322,
  29.46945492806583,
  102.24459712685095,
  41.12275630692974,
  6.6719635584127674,
  110.69055943650886,
  228.96045833314355,
  6.51651349414073,
  10.843073802257575,
  10.110261060396558,
  7.516536059594127,
  7.147678427959863,
  11.498108948890188,
  77.84171886428932,
  58.28900402519001,
  39.79299893595178,
  96.20372095280756,
  69.1724976588664,
  23.750907548738354,
  6.252513071890919,
  44.64535189066154,
  9.376191502703486,
  34.17020809979182,
  39.8741001375511,
  20.63414102952904,
  8.335150925109934,
  25.813337165548567,
  10.576002151515361,
  5.324126930341678,
  97.83932753890085,
  24.294792464321212,
  93.7812136750128,
  39.116434751511186,
  16.675990078357046,
  88.01702599711518,
  32.66879312901625,
  89.10860079246407,
  27.594867061153096,
  31.01760277897975,
  7.887068169157723,
  14.92647202087665,
  15.615817247997839,
  51.789081444436036,
  33.959268727929306,
  22.91985033134695,
  51.64026536488472,
  8.01791625522131,
  35.40306319928818,
  156.88821685805033,
  9.146787521338487,
  6.290432215728766,
  125.86227009929267,
  32.06036861177417,
  53.00435486744217,
  6.043881677488489,
  42.60987894384768,
  120.96697446291563,
  53.56839170152092,
  5.445225985469572,
  10.431395613628844,
  58.17035926230734,
  42.99555130513441,
  6.6642455857204475,
  11.168018450438,
  25.532552010210814,
  19.689447108991835,
  17.524387391947783,
  19.55212414753965,
  7.079278505288826,
  24.657682728096812,
  29.22295533164956,
  6.003370524791343,
  28.59439664529859,
  153.8087198798159,
  17.86646735726233,
  100.59607538761418,
  136.46283028327977,
  6.273889867406289,
  183.8744551224653,
  43.926633821223035,
  23.807604716859107,
  52.572638966366405,
  107.54586394540995,
  53.073939969998364,
  40.95730382658491,
  14.297051575578589,
  36.705622648118975,
  15.489658484792963,
  45.544627970500144,
  108.25080523661097,
  12.045894086567241,
  18.82275814450937,
  40.447185168777445,
  72.69234609727621,
  28.957450863253055,
  89.09193763462116,
  69.26511707333464,
  52.12826840805946,
  105.97248009747663,
  24.71006773504476,
  6.615724045093962,
  38.95273399486335,
  23.617091228379515,
  5.542988620937467,
  31.7991050186689,
  8.995942330408488,
  176.5495347256891,
  6.234650523586943,
  19.236323271931933,
  5.312113059171831,
  12.980005640377659,
  140.47715832597345,
  109.97914872695404,
  25.48527871948557,
  27.884424147593194,
  6.289512289891066,
  50.08977361861045,
  11.608163153690981,
  26.462111074065266,
  13.927865004684602,
  7.3022543197622864,
  8.524219738475079,
  52.07943080778101,
  21.279786557519145,
  17.66432269755832,
  11.522845463577987,
  8.93283532971391,
  34.50895901041959,
  39.80426715013291,
  62.17654697789925,
  49.88482558293417,
  6.08368866025267,
  16.991811347639135,
  9.176475566574483,
  37.33761558860229,
  74.97984286531754,
  48.88550122841634,
  10.185810200964234,
  13.711900607875883,
  10.156867299279295,
  29.16122290640686,
  13.79255490299566,
  68.50916047843575,
  12.752727222747184,
  11.144446005659084,
  7.115894656315024,
  32.209759154141395,
  50.18522031979376,
  19.144592673404816,
  33.926451187125316,
  84.67245875804198,
  22.2356207471094,
  10.576712785890933,
  14.50739094047482,
  9.346857592635573,
  21.640047878045056,
  15.956878803880418,
  6.246135368968968,
  47.482792030189735,
  328.98981815647835,
  19.261139174294765,
  45.33088099576746,
  55.933661687859725,
  11.114228412551258,
  13.512386339955809,
  25.24787608992044,
  59.14961944476913,
  22.905120363684794,
  13.891397854826806,
  65.3949650866028,
  83.27717189883589,
  158.9455942157528,
  56.2937382677257,
  19.606899925577643,
  82.55026200982032,
  20.880156925777406,
  6.152575774114602,
  52.87620633233132,
  12.079102617513728,
  20.505580177970533,
  16.82753595845946,
  117.21772669193572,
  30.173683254735415,
  17.217423679258687,
  7.858171713084091,
  14.08831159433832,
  5.565348336016015,
  12.552443839528934,
  25.691355953270538,
  53.39862287539856,
  39.94515272603766,
  50.77449085279253,
  5.251631814585299,
  51.36063314871676,
  27.054548006621,
  42.93135107541445,
  13.042320824040967,
  38.94119007439107,
  11.144558754842182,
  92.68322931856662,
  18.97929379309364,
  6.912872776632882,
  37.33817091681863,
  57.29560793633324,
  6.541663270948435,
  8.9583217551817,
  74.32972366732852,
  14.140485442797875,
  33.78701765844463,
  19.717480382119135,
  27.123872689534732,
  61.7588220920513,
  14.523146327413878,
  23.922619236500225,
  10.524709933429127,
  44.506801835712075,
  9.133859605102945,
  41.774104081097335,
  15.637719791498562,
  7.95000450316414,
  20.683192917795644,
  46.471162797043924,
  68.29503440851497,
  19.331465650822423,
  80.75747116634749,
  52.67749964093652,
  5.595141066061806,
  8.70433261905094,
  21.941461113579024,
  55.64705553365039,
  5.753322795152172,
  50.88161940674756,
  93.53654113240425,
  46.27715075107495,
  47.241021412927665,
  36.17471792541111,
  52.43536290417473,
  36.81032822712951,
  5.929022504706373,
  16.620152567450507,
  17.042161528676985,
  5.147225244121349,
  15.392754998128405,
  60.542797456236634,
  18.05436732871967,
  37.33168933760188,
  5.859198364572408,
  18.644646682734045,
  20.546879513972723,
  37.8523181725523,
  13.023912759924865,
  35.43682901855933,
  20.069275406330537,
  15.223918110283714,
  26.907083211472383,
  26.481562071959825,
  79.84973667457616,
  27.02321109646978,
  80.06084596224814,
  25.296468013396886,
  5.283959008515075,
  43.99361177737693,
  50.237691467659026,
  52.71513265995131,
  68.36976420458642,
  46.37473661916552,
  25.28726829855669,
  119.89988444841401,
  18.73400665211349,
  35.943513482928225,
  36.210555347685855,
  5.806325352253628,
  85.56712408895677,
  40.5039646438503,
  21.311519652857456,
  12.643042285563753,
  8.987343587627201,
  15.894972587846649,
  119.70250995747693,
  7.391952283345661,
  14.312142661494214,
  38.33810008921233,
  57.22146334905344,
  26.112882090615813,
  29.08788315134973,
  60.32957979414643,
  10.569114193689918,
  78.82062341725094,
  49.20380341159883,
  22.735121178256865,
  27.29899539570856,
  171.60564229959652,
  47.03787626215502,
  84.64026715964413,
  28.255436751363685,
  31.041271486414434,
  26.03138330688587,
  30.69018505212393,
  17.33945898218055,
  20.741856365916416,
  34.39656128436948,
  32.25041486248019,
  17.219253544890268,
  12.162679834210863,
  76.09708704849265,
  47.073834808572336,
  72.099192350003,
  60.254810665237564,
  33.57402852192815,
  11.576817532025064,
  19.85440546798221,
  67.56277601313558,
  25.064635587018305,
  15.443547837581601,
  65.57509063286506,
  22.758926589522268,
  105.46137886924613,
  37.55674567362179,
  100.4918702667481,
  28.22609170658376,
  23.58628829630893,
  21.014816762304413,
  119.55624673828223,
  120.86386619604761,
  32.18507175160181,
  171.79164374499774,
  60.88162981203118,
  15.342634910298637,
  31.16039991158481,
  35.260241571303915,
  27.344143147044896,
  43.74920005398438,
  11.395248413345847,
  117.22010438751501,
  7.8912403139854925,
  14.92401327293268,
  5.393889840367488,
  8.279303087768238,
  59.922764827281384,
  56.19924026232134,
  43.09742253336949,
  7.920815783515331,
  6.312291357021164,
  42.67192747405723,
  58.365192078672614,
  38.32677803151704,
  5.300337379615727,
  42.7261427461295,
  22.148412869993265,
  36.16459499434433,
  69.33811249279483,
  9.746703484072992,
  16.57628123928538,
  12.012023985878013,
  8.314503455211783,
  50.11866377685668,
  74.12072579375187,
  16.836359810931757,
  82.43930612589126,
  42.20531373361757,
  10.323610651744039,
  8.731902879751532,
  63.06134191757923,
  14.053730547292425,
  9.859008783918828,
  122.35038950135372,
  34.984785318756984,
  5.175787148806842,
  26.90752023822015,
  14.966923011582885,
  9.694934265212874,
  12.132592897577792,
  63.48480586647386,
  14.488540589712892,
  86.33819451925297,
  19.916519037233,
  14.054025129256633,
  30.713518193495936,
  59.61974300670012,
  15.655726451629446,
  27.127994680681066,
  19.44611177009812,
  20.74624220488592,
  17.273493357238173,
  29.407363102571487,
  41.41151004164073,
  58.922833632484,
  16.99623868604042,
  111.67970585870904,
  6.6337524708251125,
  72.15417495616964,
  56.763741009605944,
  19.466336565440052,
  5.4586903667915525,
  13.42194879163562,
  36.429114896919046,
  39.90569003842483,
  58.50488771446582,
  8.140508898336048,
  77.11988440481217,
  66.49638457088298,
  66.1187092494653,
  37.1842015976012,
  27.733046819811918,
  25.062580631470595,
  55.668730445523884,
  9.39980948752519,
  49.902978072151924,
  5.339364413536215,
  28.698457873769797,
  41.09262571917477,
  20.785397823265516,
  24.20355309785876,
  14.812128233324046,
  8.817184413313417,
  9.00205739034952,
  14.986275771160003,
  12.623494804779385,
  57.01322277354186,
  55.52008643593597,
  48.865659551273254,
  20.131807390406497,
  20.749650525036113,
  5.026300861271223,
  8.45599248133933,
  92.45525209630476,
  37.667431354283536,
  143.7688430488426,
  18.751578768346857,
  10.700749338077976,
  30.781118280975498,
  14.192862165299847,
  93.09828835541049,
  20.521512926586272,
  51.44049529138689,
  22.63607317752428,
  11.113469436106252,
  62.12643074009492,
  20.954039679420067,
  34.070810908660675,
  18.99759175034995,
  7.857577025138012,
  43.47849580956395,
  63.89474244831867,
  38.90770493838187,
  77.43728837141978,
  5.459167752419965,
  127.29906349863414,
  5.383095064172092,
  7.693472978290682,
  11.082827198674643,
  68.0748354022724,
  57.21124166674078,
  44.14997484735749,
  5.602210525198436,
  19.275917747106114,
  41.09196018570641,
  22.87034700814646,
  42.44229379230722,
  69.71207864950152,
  21.93270739540813,
  67.79745487091058,
  46.26751252021307,
  68.89592531627014,
  90.61133318937571,
  21.41220584214014,
  27.256701390250292,
  36.30318168645079,
  79.34608610252965,
  9.044447224935139,
  5.950268487147574,
  40.38240431743288,
  60.89803801676996,
  28.10544854767372,
  74.76619215455344,
  50.86775087383714,
  85.08042335151902,
  10.830533182773305,
  23.115580720485625,
  31.723429041988723,
  5.716196052743241,
  8.96571407769283,
  70.34579197806306,
  8.966016800407095,
  10.030751580815174,
  47.5688196455782,
  13.181344835298118,
  44.28130498181924,
  7.522916290595509,
  7.555249636135615,
  102.01260700723242,
  48.192499493770384,
  27.044265862519545,
  25.85161728501255,
  38.873623105381476,
  81.51277375987193,
  87.99466326315533,
  31.240924251921584,
  13.329046932719358,
  9.404664704760718,
  14.435792218680437,
  72.0443409836828,
  7.794475927166077,
  77.56027135258952,
  30.830992614557633,
  25.495027606236977,
  6.917058026196925,
  130.56705327642095,
  26.983162966713323,
  31.667837752358892,
  29.600001058099423,
  5.714349567598257,
  6.904173212543638,
  74.86827333964791,
  27.932664063219345,
  28.85069525791109,
  49.38916189379995,
  69.25670652391995,
  62.779149276582444,
  52.595044306505464,
  89.15607907736654,
  9.34721153441597,
  24.34719660791221,
  20.41448328978249,
  16.04441135592343,
  10.433584803268053,
  41.329080110011425,
  9.058154687872289,
  107.64005027072655,
  21.476287962797628,
  27.566332335837902,
  17.96942457506728,
  58.09634975946812,
  17.05141931783922,
  39.30135963240701,
  26.150166260130817,
  11.777513349536502,
  21.855166173760516,
  5.504549418459996,
  13.189733326364223,
  5.86941663139027,
  16.543155521487638,
  19.524739642591005,
  8.615995155354518,
  5.209876062915035,
  20.88846153496624,
  29.353452219197568,
  38.6339454632909,
  14.86216835735905,
  28.08579087417601,
  15.660489209010343,
  11.498195620167337,
  11.654182408129442,
  7.687054459033075,
  123.77114121959707,
  15.601458594924091,
  6.970487333764528,
  7.885902320216253,
  5.068858326898298,
  88.46459860141735,
  16.206535031841753,
  28.775431676906226,
  12.767356441026058,
  51.7642611240329,
  12.593596633462976,
  20.907895896471196,
  25.17439843477087,
  25.372131542085338,
  6.666500213041673,
  20.486355234984117,
  70.74549648886186,
  71.72758345826381,
  5.584017223572758,
  15.88439821731767,
  9.327397335839784,
  6.422734138365126,
  65.84960692681398,
  27.100327724312436,
  129.7550856984676,
  17.163285132520386,
  62.78321675176749,
  13.5993367452964,
  43.29506660258941,
  120.87570479279883,
  59.97229167923042,
  14.144739312122685,
  61.08029609878962,
  17.64487537239363,
  13.968588444332157,
  38.30842664824547,
  8.76121376048288,
  30.38382556535369,
  19.355231618011253,
  31.61707215931181,
  101.70973825363468,
  25.493149793521606,
  29.778798011146367,
  43.53481729816229,
  10.31811299822585,
  69.4129713158506,
  14.166196728180852,
  25.196826886038085,
  33.36760750559982,
  14.025261656043163,
  16.14768328979134,
  32.477808164586584,
  74.26724377052109,
  15.434845843483465,
  8.183968653324403,
  26.538158010686168,
  6.879197658658789,
  60.398272696129474,
  5.310353242136236,
  62.89231453628379,
  24.24569009623564,
  9.603213525856418,
  8.365369051724706,
  14.595376930320295,
  7.180325947575875,
  13.569150077236507,
  13.548306605005875,
  13.756401744802467,
  25.686402971722902,
  41.400325223071654,
  6.345104632887465,
  7.617117158524453,
  100.20620548845756,
  46.074423169882195,
  6.060119770741594,
  50.729143167703825,
  28.291206832572946,
  10.612986744057762,
  76.65606729944457,
  20.066558174302738,
  15.070747671027846,
  17.471548155181615,
  52.7980943118385,
  18.791149765462272,
  6.145565298520987,
  13.060010486809109,
  14.139596835161989,
  32.43691765651124,
  122.06776761250734,
  8.639387517868306,
  85.18182413679119,
  42.787166476670336,
  36.11863781779974,
  7.924546047289133,
  32.740198027067784,
  19.378199770780938,
  32.69629215008728,
  23.530498211989027,
  140.55604676099674,
  18.35643278368189,
  5.168840362449723,
  11.641880772650373,
  24.931413911878312,
  37.8375095301648,
  36.51837841597658,
  7.214595175924872,
  89.05749480457942,
  77.16388397722504,
  40.043820630679434,
  20.037731511293774,
  53.447818244001375,
  25.590491600348273,
  10.364150230970486,
  55.569553515506676,
  6.010412243570199,
  18.929854136344957,
  5.100688294790798,
  54.429358634056,
  5.162597844839081,
  6.803363126210574,
  43.24178418444137,
  53.158592393886494])

Energy spectrum

caption, xmean,count, xle, Etot = energy_spectrum(E_tot)
hist_file(E_hist, xmean, count, caption);

([5578977.992637578,
  5613539.99644256,
  5621984.067356959,
  5625185.491974772,
  5630368.1422281265,
  5632342.76281048,
  5700317.071617094,
  5743468.085524118,
  5744836.701108547,
  5774398.334295918,
  5774767.373556603,
  5774883.6945538055,
  5777203.405797355,
  5793085.617795732,
  5821571.743378873,
  5855827.613748997,
  5925625.695451189,
  5933944.544405776,
  5934801.5625,
  5940192.840524625,
  5946482.516304447,
  5969773.381161929,
  5981136.763201982,
  5984070.243328004,
  6013602.588212368,
  6013635.187115682,
  6031473.283657751,
  6037470.24548906,
  6042843.5027544135,
  6064294.474320638,
  6094475.0,
  6094905.9829924675,
  6095592.947819843,
  6095976.71663147,
  6096367.1875,
  6100604.671760609,
  6127496.660020787,
  6127514.66433957,
  6131121.586458964,
  6135403.121850532,
  6136735.932735263,
  6141650.962756152,
  6142076.310749517,
  6145591.865585281,
  6145618.491309015,
  6149554.421623171,
  6150099.041201457,
  6151184.375,
  6153137.020367047,
  6155257.734810135,
  6159103.125,
  6159982.962820449,
  6164579.583200394,
  6165053.125,
  6165408.029470348,
  6167709.375,
  6175208.476366167,
  6182233.18007594,
  6184032.704570248,
  6185279.554791648,
  6185349.042431714,
  6187329.873797793,
  6187452.158667843,
  6191708.011044823,
  6192119.60121315,
  6192157.165202618,
  6199457.507883905,
  6218398.34190127,
  6218690.570040212,
  6218746.875,
  6220405.099065951,
  6220965.103675448,
  6222523.60000546,
  6223349.563405839,
  6224387.228177192,
  6224443.575621733,
  6226151.5625,
  6228858.775595083,
  6230673.432989561,
  6253256.072526208,
  6265052.452044818,
  6274830.950682911,
  6277836.20900789,
  6314249.280928683,
  6315299.579220695,
  6315736.783741842,
  6316443.123659642,
  6316910.9375,
  6317419.647893409,
  6318677.500032619,
  6319442.44569348,
  6319450.113932871,
  6319712.830529361,
  6320429.6875,
  6322545.811595063,
  6323816.559097993,
  6323824.821548586,
  6331484.375,
  6355723.014527185,
  6357360.588598309,
  6363796.557840736,
  6395125.521991844,
  6395641.637876963,
  6400400.543865813,
  6401275.067899309,
  6423478.595334258,
  6426760.9375,
  6434744.811132539,
  6442036.355002249,
  6444163.76806014,
  6452397.865784013,
  6462296.591062106,
  6640179.5698098745,
  6924175.0,
  6932273.239346904,
  6978174.583763218,
  6985594.689982541,
  6986289.0625,
  7009900.567918236,
  7049085.756867959,
  7050059.154200402,
  7053084.837842204,
  7054303.048294704,
  7063516.622922211,
  7073588.617017695,
  7200593.640329458,
  7204766.782160621,
  7240366.622016894,
  7317253.268836167,
  7364354.6875,
  7365233.114811757,
  7367649.523492468,
  7493069.116851672,
  7503121.868162566,
  7503184.560158351,
  7531570.1224347735,
  7573246.109365999,
  7614677.878747302,
  7642043.848840267,
  7642056.026815211,
  7740341.077747443,
  7837390.625,
  7871174.495460565,
  7872129.6875,
  7908156.925440548,
  7918716.881316441,
  7919718.632875429,
  7922573.709157177,
  7923678.177498328,
  7924698.879963058,
  7925407.8125,
  7926249.4093552325,
  7927146.96573762,
  7931012.345965633,
  7939362.115827618,
  7941796.346337401,
  7942101.001009151,
  7943187.988806127,
  7946514.697059679,
  7948926.203486298,
  7949316.016096984,
  7977941.478919659,
  7997917.289809064,
  7998522.693319561,
  8001081.420265849,
  8004648.631118895,
  8004875.722019136,
  8005660.933549159,
  8006579.534724715,
  8006725.0,
  8006839.0625,
  8006954.982699717,
  8008182.244234026,
  8010172.892383554,
  8025226.5625,
  8027709.559995733,
  8027712.273284319,
  8028437.681708809,
  8028718.75,
  8028848.4375,
  8029351.5625,
  8030138.62628424,
  8030164.0625,
  8030282.884151183,
  8030343.75,
  8030579.296380469,
  8030788.566129072,
  8031312.5,
  8032257.8125,
  8039522.919105169,
  8045755.033987068,
  8045815.030558283,
  8051632.177019472,
  8055023.17636913,
  8055675.966119183,
  8055677.96209952,
  8058149.7678323435,
  8058416.82410652,
  8059667.893819579,
  8059670.3125,
  8063958.710953828,
  8072553.726499772,
  8094066.422983345,
  8129896.6949967425,
  8141425.052780143,
  8142815.625,
  8158021.623618728,
  8167913.766938603,
  8174902.733161797,
  8181806.800345287,
  8187138.001237432,
  8189534.981143247,
  8189853.418130742,
  8190473.4375,
  8190840.761745352,
  8257782.8125,
  8262367.272545,
  8263788.017466279,
  8268400.745466446,
  8268609.375,
  8268650.712523266,
  8271246.875,
  8274221.875,
  8274650.0,
  8274695.3125,
  8274698.026401529,
  8274734.146565861,
  8274752.336769589,
  8274771.139793693,
  8274815.314464465,
  8274909.577474334,
  8276568.75,
  8277005.824703687,
  8278117.1875,
  8292855.987967529,
  8329567.989503332,
  8369545.9785145745,
  8370513.431092635,
  8371887.753969119,
  8374074.405507117,
  8374076.643530111,
  8382171.457617436,
  8422064.237096388,
  8422096.821402773,
  8448323.72457198,
  8448454.343248893,
  8452726.160929244,
  8454011.583374472,
  8454603.548884489,
  8454943.75,
  8455778.728842134,
  8482536.891196294,
  8482783.0513977,
  8484149.487759253,
  8485708.153880486,
  8487173.382844176,
  8488025.0,
  8488662.5,
  8495114.019980155,
  8495364.513272393,
  8495590.435909467,
  8496115.522262499,
  8502990.527542615,
  8505265.625,
  8506783.675371563,
  8506812.495052168,
  8509096.150287287,
  8510243.75,
  8510427.786932219,
  8511181.25,
  8512575.838905351,
  8513537.51068289,
  8518229.666317971,
  8518587.5,
  8518845.066480987,
  8538177.352235738,
  8538651.698696647,
  8541294.02065175,
  8550751.333775036,
  8550788.678866152,
  8552459.256242624,
  8552723.455720842,
  8553490.139605707,
  8563675.928755844,
  8568046.706448179,
  8569368.75,
  8570442.8937564,
  8572202.725852901,
  8576441.562753767,
  8578176.267249139,
  8579606.241359796,
  8580919.879156416,
  8583482.621314054,
  8585365.625,
  8585791.663102852,
  8585818.44421204,
  8586276.206291188,
  8586568.75,
  8589244.446631467,
  8589526.5625,
  8591485.654464064,
  8600761.594623223,
  8606515.061278822,
  8610595.910508148,
  8622675.934620805,
  8625634.908882054,
  8625814.485315993,
  8626881.52629706,
  8633973.971662201,
  8639836.424902648,
  8664127.054110654,
  8664489.716332259,
  8664549.756569194,
  8664722.522404205,
  8665114.948647423,
  8665252.087675111,
  8665749.782820417,
  8666228.125,
  8666942.1875,
  8667709.605118735,
  8668166.404484307,
  8668435.85852967,
  8669265.625,
  8669269.705527764,
  8669285.9375,
  8670121.88575467,
  8670287.5,
  8670294.517241456,
  8670517.89229598,
  8671609.375,
  8672554.104710143,
  8672950.339512775,
  8673110.460658483,
  8673431.25,
  8674107.46158853,
  8674760.951248473,
  8675641.758034034,
  8677970.237481065,
  8678623.8448693,
  8683148.323229674,
  8691526.043281795,
  8707683.95163068,
  8713364.0625,
  8725007.569254583,
  8725190.625,
  8725530.76560904,
  8725793.054255715,
  8727321.962750567,
  8728206.780048713,
  8735364.953771826,
  8737135.380865166,
  8741326.743886946,
  8752984.99766059,
  8754162.15317943,
  8761408.526926685,
  8792685.360675436,
  8792758.041077849,
  8793653.832615513,
  8793996.875,
  8794022.188779937,
  8794154.6875,
  8795154.6875,
  8795201.5625,
  8797156.25,
  8798655.251237357,
  8802349.879995534,
  8802656.63022709,
  8802989.522512034,
  8803828.891342113,
  8804015.738557713,
  8804023.4375,
  8804145.943872036,
  8804186.860842645,
  8804192.780479454,
  8804219.813627753,
  8804569.712109134,
  8806134.96419869,
  8806203.889612991,
  8806254.6875,
  8806975.906851102,
  8807864.134671347,
  8808102.117476793,
  8919282.607320037,
  8927556.573354574,
  8953757.731735593,
  8954361.746108068,
  8954362.31208534,
  8954415.39249643,
  8955966.60381047,
  8958651.5625,
  8962332.555338075,
  8971112.321127882,
  8977868.233932119,
  8980165.436467165,
  8983996.92816116,
  8984599.17041027,
  8985237.5,
  8985248.935334267,
  8985299.896004494,
  8988940.625,
  8989925.0,
  8989927.606602695,
  8991127.35511087,
  8991185.500495778,
  8993061.513141083,
  8996273.120763056,
  9004222.812932868,
  9007014.95784737,
  9028259.238826694,
  9030218.242034517,
  9032344.708575925,
  9035695.65798077,
  9036945.3125,
  9037044.923101418,
  9037098.69224798,
  9038609.375,
  9039574.616101801,
  9039973.000686746,
  9040032.284541702,
  9041434.707256664,
  9041732.359365825,
  9042234.361240216,
  9042501.5625,
  9042831.284776246,
  9042887.5,
  9043064.575609773,
  9043202.546242805,
  9043281.25,
  9043592.066514777,
  9044362.991936032,
  9044486.512079779,
  9044854.017704967,
  9045270.3125,
  9045492.95490584,
  9045494.87490307,
  9045720.127382657,
  9045908.082004739,
  9046732.13088164,
  9046779.482426535,
  9048838.858840046,
  9048903.789355332,
  9049677.784857232,
  9052394.930674184,
  9053596.875,
  9053828.930609608,
  9062087.914366238,
  9067817.172224868,
  9070774.077437205,
  9107027.744620299,
  9109396.695817403,
  9129387.419966767,
  9129628.125,
  9133234.420686716,
  9137449.976875912,
  9157097.79604507,
  9174989.89398588,
  9178357.8125,
  9180308.896938387,
  9205695.091304233,
  9227506.414677413,
  9233628.954627037,
  9279602.294272551,
  9289142.182902094,
  9289290.504566358,
  9299359.361218933,
  9303648.4375,
  9307773.782939412,
  9309326.5625,
  9309650.800261535,
  9349581.598764736,
  9351163.369059006,
  9358130.934509596,
  9360859.375,
  9361247.1715868,
  9365434.482828114,
  9386427.518439142,
  9387284.375,
  9387340.625,
  9387694.203454724,
  9387923.469339475,
  9388744.119543774,
  9389051.614938,
  9389479.370882215,
  9389733.16428069,
  9390684.004603539,
  9391330.75013092,
  9411608.573555443,
  9412030.666474579,
  9412792.944299256,
  9412845.3125,
  9412851.308301268,
  9413583.376862906,
  9415168.94505935,
  9415216.51633216,
  9415425.442093844,
  9415459.87429085,
  9415565.631651614,
  9415673.838264123,
  9415698.997529536,
  9415723.4375,
  9415790.500855843,
  9415873.614106368,
  9416124.885779316,
  9416262.24589332,
  9416313.714787347,
  9416316.916951617,
  9416423.613059886,
  9416474.545621373,
  9416500.829499535,
  9416602.546723524,
  9416828.932083739,
  9417016.571739135,
  9417062.328590088,
  9418001.5625,
  9418397.90640377,
  9418610.835803391,
  9418635.250859113,
  9418748.08130043,
  9432352.294239253,
  9550716.457110025,
  9552684.303133195,
  9553108.909647174,
  9654869.71805141,
  9745972.988154227,
  9748710.47078961,
  9764723.150676485,
  9765188.735763712,
  9791240.625,
  9798694.121990992,
  9799732.710472522,
  9800043.708696079,
  9802784.171198906,
  9802805.990969595,
  9803242.267307982,
  9850576.5625,
  9850640.550180549,
  9872050.018251656,
  9905679.190650804,
  9909968.176003896,
  9911297.66728433,
  9913121.840636093,
  9916705.572066901,
  9948380.833003595,
  9965781.261695607,
  9972839.846731849,
  9973690.479978068,
  9973754.6875,
  9974052.746241923,
  9974161.851776069,
  9974490.047096923,
  9974545.950481862,
  9974584.68219173,
  9974604.730765797,
  9974646.83061326,
  9974682.513091875,
  9975133.547908572,
  9975212.208521167,
  9975220.735421192,
  9975229.476164848,
  9975235.811844306,
  9975765.021756265,
  9975856.501909183,
  9976007.419921272,
  9976057.80337224,
  9976065.266040968,
  9976886.342057634,
  9977187.243001143,
  9977210.144756453,
  9977274.089917224,
  9977470.3125,
  9977504.185187697,
  9978905.05714484,
  9979243.707714986,
  9982509.31431189,
  9985512.5,
  9992113.654304542,
  9993020.555281322,
  9996762.471096417,
  9996787.381402094,
  10002379.336804705,
  10024478.125,
  10027391.971597543,
  10035632.623142052,
  10037956.857204715,
  10038008.361725478,
  10040890.37532692,
  10041982.8125,
  10065853.38125877,
  10091403.838716848,
  10092895.3125,
  10092963.321090968,
  10095299.720654031,
  10097807.304352304,
  10097856.253012737,
  10098667.1875,
  10099229.6875,
  10099909.014946127,
  10100182.936398685,
  10100247.188947475,
  10100470.178203166,
  10101104.057642618,
  10101162.825340625,
  10101854.6875,
  10105054.78391054,
  10131292.45049253,
  10161616.416436505,
  10161626.269829186,
  10161729.834783094,
  10169688.595853087,
  10170138.014750518,
  10179129.221269872,
  10179714.788008817,
  10184406.21191207,
  10188168.196301162,
  10188721.875,
  10190524.772325411,
  10191598.302842084,
  10191712.078098169,
  10199184.375,
  10199235.857941551,
  10208423.172623368,
  10208909.375,
  10210021.875,
  10211332.8125,
  10226490.625,
  10232635.438595347,
  10235212.585390566,
  10237767.544899004,
  10270728.125,
  10283002.062005082,
  10290050.0,
  10292251.5625,
  10311263.88779654,
  10386710.098369364,
  10391426.332413316,
  10466773.854691127,
  10477916.968040938,
  10480777.584047493,
  10490746.993626973,
  10503769.415995158,
  10520482.016492838,
  10541145.250601705,
  10541158.407772304,
  10548296.875,
  10573725.0,
  10773496.875,
  10774728.905083535,
  10774751.5625,
  10776792.206136825,
  10777665.020240204,
  10777668.13461999,
  10778846.71205643,
  10780182.317867313,
  10780895.3125,
  10781879.34564069,
  10782364.0625,
  10782616.685812997,
  10782751.512144936,
  10784220.47316309,
  10785047.461768515,
  10785567.649202945,
  10785963.83995435,
  10786879.554680398,
  10789251.242923414,
  10790342.008619381,
  10791541.585273359,
  10794124.449462075,
  10794867.909633797,
  10794883.177574959,
  10794896.152109306,
  10798670.88748433,
  10798675.251460712,
  10800201.774210311,
  10800576.472235681,
  10800866.467850756,
  10801371.854027035,
  10801583.580403736,
  10801655.203027135,
  10804001.787265554,
  10804373.7500902,
  10804408.038444053,
  10804666.718057577,
  10805223.262449266,
  10805324.392792817,
  10806161.656789578,
  10806224.997911619,
  10807836.321072653,
  10808036.644279381,
  10808415.625,
  10808581.873706874,
  10808873.4375,
  10808939.833839264,
  10808980.868800774,
  10809117.556962334,
  10809218.567149037,
  10809468.243870266,
  10809663.455287416,
  10809883.819216076,
  10810860.650136447,
  10811066.631674938,
  10811229.246688634,
  10811451.55775868,
  10812146.451787293,
  10812288.74001905,
  10814082.73614174,
  10814429.6875,
  10814470.794442,
  10814634.175857754,
  10814822.984567305,
  10815081.21234314,
  10815575.934037207,
  10816437.5,
  10817226.336476753,
  10817231.25,
  10817241.190187208,
  10817316.994471196,
  10817685.920652514,
  10817708.61376789,
  10818650.44962985,
  10818785.0410984,
  10818949.73800539,
  10819542.563609054,
  10819618.57517276,
  10820283.905993085,
  10821180.385645043,
  10821563.474445011,
  10821962.366196258,
  10822690.405543769,
  10823060.742860926,
  10824867.355729358,
  10825224.736878874,
  10825530.41632611,
  10825575.0,
  10825796.616920745,
  10825852.074789431,
  10826027.925790953,
  10826181.140819779,
  10826309.928126782,
  10826581.25,
  10826678.736373166,
  10827726.356233587,
  10828785.467177365,
  10828997.179867625,
  10829662.646082487,
  10829965.531124143,
  10830422.094331188,
  10830576.343721554,
  10831079.266695555,
  10831707.462456234,
  10831860.775297178,
  10832075.0,
  10832606.842710985,
  10833504.725712355,
  10834208.209129078,
  10834461.755489081,
  10835198.16407682,
  10835381.25,
  10835405.036137272,
  10835788.706193324,
  10836144.412098594,
  10836521.80788944,
  10836907.84599894,
  10836979.299859766,
  10837014.0625,
  10837414.803512307,
  10838227.097552575,
  10838352.524900647,
  10838645.084898634,
  10838932.537461119,
  10839782.707676722,
  10839966.237466928,
  10840163.452415587,
  10840168.75,
  10840573.001947701,
  10840835.183819659,
  10841051.5625,
  10841305.140630526,
  10841315.26370543,
  10841873.732887015,
  10841873.955943927,
  10842635.052716026,
  10843007.96124787,
  10845167.493300281,
  10846154.474631691,
  10846417.784582015,
  10846660.9375,
  10850361.280366791,
  10851604.490230924,
  10851920.3125,
  10852854.71921794,
  10852988.403083177,
  10854441.81241589,
  10854476.873962175,
  10854951.012120415,
  10855042.432351183,
  10855117.968278056,
  10856024.359284345,
  10856636.050576696,
  10856885.9375,
  10857321.755933389,
  10858837.696414525,
  10859606.116575606,
  445897168.75],
 [10.926932224398904,
  20.299863187360607,
  128.03679543984842,
  22.054290030078057,
  92.0134442811754,
  8.616197603256872,
  6.076130750777888,
  99.513277837992,
  62.942603518181855,
  14.562212679960686,
  69.58291917703299,
  7.9567872446315615,
  6.820890981208564,
  67.7001006400315,
  6.73751087261893,
  117.25589082591692,
  186.77098044642872,
  15.759888942893621,
  111.42139242152703,
  209.23843133454497,
  8.135918106382611,
  72.9941144126226,
  25.883289989482332,
  36.68592379563776,
  64.15916495810745,
  53.54018954297891,
  10.895464639696927,
  26.62252597598342,
  72.73778335531173,
  8.7107392510694,
  49.865399405062874,
  112.20082055767892,
  82.64480473827645,
  6.286870166924756,
  48.66363464808074,
  17.15804067148037,
  9.156008998167568,
  48.4507440634721,
  59.66783805117294,
  28.462945030734435,
  28.43112284737409,
  13.060927451972757,
  24.3872693557942,
  42.247653860659256,
  6.061604517463069,
  24.196628613696685,
  12.223600475285949,
  87.33087495537146,
  112.04535780831122,
  99.35487052097429,
  34.3612146902758,
  69.82759419328755,
  93.7767087179709,
  80.87655749057232,
  65.84386632806763,
  46.10461411743419,
  95.29772568023517,
  54.90189286240101,
  90.08621857118047,
  26.15196044445469,
  5.422114749261968,
  20.46881669574627,
  9.892079739524835,
  64.04483853197218,
  19.616002262174746,
  22.426241203402252,
  23.69195598151165,
  26.767377412322464,
  27.48242570625202,
  78.17984072236509,
  16.782526933300517,
  5.776388073977832,
  17.16245783712184,
  22.14981813755204,
  41.21862258940446,
  32.30695594607249,
  52.61387621741336,
  16.04861373266941,
  55.8379853597047,
  25.453079411126218,
  78.58967712147752,
  23.759468818820626,
  6.116796188298128,
  59.45093622154312,
  6.682311490039749,
  55.73381667475616,
  65.93507541908905,
  74.76998507261524,
  7.963836959850351,
  12.996342339540242,
  115.28563924400606,
  11.310586744054634,
  43.415073325361774,
  89.11735884359675,
  7.794309767515815,
  8.851141962612628,
  14.302330732534013,
  36.26773275395881,
  130.4621411981585,
  23.145561783174436,
  139.05611076736474,
  17.39055387202637,
  13.192164245260043,
  21.52157689636502,
  44.57887992328304,
  137.21496266899408,
  29.43156867413411,
  21.479469924041624,
  11.849662875070816,
  78.45085521222896,
  78.00925312492446,
  64.21466078080728,
  25.522473652892575,
  44.0564159458759,
  7.070319545796551,
  22.37146180037758,
  21.563094880394516,
  107.55994206406913,
  145.031435332794,
  14.350142723374486,
  61.36056062936864,
  97.6706598756864,
  36.30821344342823,
  5.229613973701861,
  95.05866320511922,
  10.960505086983037,
  10.54819837019146,
  10.34560973780322,
  29.46945492806583,
  102.24459712685095,
  41.12275630692974,
  6.6719635584127674,
  110.69055943650886,
  228.96045833314355,
  6.51651349414073,
  10.843073802257575,
  10.110261060396558,
  7.516536059594127,
  7.147678427959863,
  11.498108948890188,
  77.84171886428932,
  58.28900402519001,
  39.79299893595178,
  96.20372095280756,
  69.1724976588664,
  23.750907548738354,
  6.252513071890919,
  44.64535189066154,
  9.376191502703486,
  34.17020809979182,
  39.8741001375511,
  20.63414102952904,
  8.335150925109934,
  25.813337165548567,
  10.576002151515361,
  5.324126930341678,
  97.83932753890085,
  24.294792464321212,
  93.7812136750128,
  39.116434751511186,
  16.675990078357046,
  88.01702599711518,
  32.66879312901625,
  89.10860079246407,
  27.594867061153096,
  31.01760277897975,
  7.887068169157723,
  14.92647202087665,
  15.615817247997839,
  51.789081444436036,
  33.959268727929306,
  22.91985033134695,
  51.64026536488472,
  8.01791625522131,
  35.40306319928818,
  156.88821685805033,
  9.146787521338487,
  6.290432215728766,
  125.86227009929267,
  32.06036861177417,
  53.00435486744217,
  6.043881677488489,
  42.60987894384768,
  120.96697446291563,
  53.56839170152092,
  5.445225985469572,
  10.431395613628844,
  58.17035926230734,
  42.99555130513441,
  6.6642455857204475,
  11.168018450438,
  25.532552010210814,
  19.689447108991835,
  17.524387391947783,
  19.55212414753965,
  7.079278505288826,
  24.657682728096812,
  29.22295533164956,
  6.003370524791343,
  28.59439664529859,
  153.8087198798159,
  17.86646735726233,
  100.59607538761418,
  136.46283028327977,
  6.273889867406289,
  183.8744551224653,
  43.926633821223035,
  23.807604716859107,
  52.572638966366405,
  107.54586394540995,
  53.073939969998364,
  40.95730382658491,
  14.297051575578589,
  36.705622648118975,
  15.489658484792963,
  45.544627970500144,
  108.25080523661097,
  12.045894086567241,
  18.82275814450937,
  40.447185168777445,
  72.69234609727621,
  28.957450863253055,
  89.09193763462116,
  69.26511707333464,
  52.12826840805946,
  105.97248009747663,
  24.71006773504476,
  6.615724045093962,
  38.95273399486335,
  23.617091228379515,
  5.542988620937467,
  31.7991050186689,
  8.995942330408488,
  176.5495347256891,
  6.234650523586943,
  19.236323271931933,
  5.312113059171831,
  12.980005640377659,
  140.47715832597345,
  109.97914872695404,
  25.48527871948557,
  27.884424147593194,
  6.289512289891066,
  50.08977361861045,
  11.608163153690981,
  26.462111074065266,
  13.927865004684602,
  7.3022543197622864,
  8.524219738475079,
  52.07943080778101,
  21.279786557519145,
  17.66432269755832,
  11.522845463577987,
  8.93283532971391,
  34.50895901041959,
  39.80426715013291,
  62.17654697789925,
  49.88482558293417,
  6.08368866025267,
  16.991811347639135,
  9.176475566574483,
  37.33761558860229,
  74.97984286531754,
  48.88550122841634,
  10.185810200964234,
  13.711900607875883,
  10.156867299279295,
  29.16122290640686,
  13.79255490299566,
  68.50916047843575,
  12.752727222747184,
  11.144446005659084,
  7.115894656315024,
  32.209759154141395,
  50.18522031979376,
  19.144592673404816,
  33.926451187125316,
  84.67245875804198,
  22.2356207471094,
  10.576712785890933,
  14.50739094047482,
  9.346857592635573,
  21.640047878045056,
  15.956878803880418,
  6.246135368968968,
  47.482792030189735,
  328.98981815647835,
  19.261139174294765,
  45.33088099576746,
  55.933661687859725,
  11.114228412551258,
  13.512386339955809,
  25.24787608992044,
  59.14961944476913,
  22.905120363684794,
  13.891397854826806,
  65.3949650866028,
  83.27717189883589,
  158.9455942157528,
  56.2937382677257,
  19.606899925577643,
  82.55026200982032,
  20.880156925777406,
  6.152575774114602,
  52.87620633233132,
  12.079102617513728,
  20.505580177970533,
  16.82753595845946,
  117.21772669193572,
  30.173683254735415,
  17.217423679258687,
  7.858171713084091,
  14.08831159433832,
  5.565348336016015,
  12.552443839528934,
  25.691355953270538,
  53.39862287539856,
  39.94515272603766,
  50.77449085279253,
  5.251631814585299,
  51.36063314871676,
  27.054548006621,
  42.93135107541445,
  13.042320824040967,
  38.94119007439107,
  11.144558754842182,
  92.68322931856662,
  18.97929379309364,
  6.912872776632882,
  37.33817091681863,
  57.29560793633324,
  6.541663270948435,
  8.9583217551817,
  74.32972366732852,
  14.140485442797875,
  33.78701765844463,
  19.717480382119135,
  27.123872689534732,
  61.7588220920513,
  14.523146327413878,
  23.922619236500225,
  10.524709933429127,
  44.506801835712075,
  9.133859605102945,
  41.774104081097335,
  15.637719791498562,
  7.95000450316414,
  20.683192917795644,
  46.471162797043924,
  68.29503440851497,
  19.331465650822423,
  80.75747116634749,
  52.67749964093652,
  5.595141066061806,
  8.70433261905094,
  21.941461113579024,
  55.64705553365039,
  5.753322795152172,
  50.88161940674756,
  93.53654113240425,
  46.27715075107495,
  47.241021412927665,
  36.17471792541111,
  52.43536290417473,
  36.81032822712951,
  5.929022504706373,
  16.620152567450507,
  17.042161528676985,
  5.147225244121349,
  15.392754998128405,
  60.542797456236634,
  18.05436732871967,
  37.33168933760188,
  5.859198364572408,
  18.644646682734045,
  20.546879513972723,
  37.8523181725523,
  13.023912759924865,
  35.43682901855933,
  20.069275406330537,
  15.223918110283714,
  26.907083211472383,
  26.481562071959825,
  79.84973667457616,
  27.02321109646978,
  80.06084596224814,
  25.296468013396886,
  5.283959008515075,
  43.99361177737693,
  50.237691467659026,
  52.71513265995131,
  68.36976420458642,
  46.37473661916552,
  25.28726829855669,
  119.89988444841401,
  18.73400665211349,
  35.943513482928225,
  36.210555347685855,
  5.806325352253628,
  85.56712408895677,
  40.5039646438503,
  21.311519652857456,
  12.643042285563753,
  8.987343587627201,
  15.894972587846649,
  119.70250995747693,
  7.391952283345661,
  14.312142661494214,
  38.33810008921233,
  57.22146334905344,
  26.112882090615813,
  29.08788315134973,
  60.32957979414643,
  10.569114193689918,
  78.82062341725094,
  49.20380341159883,
  22.735121178256865,
  27.29899539570856,
  171.60564229959652,
  47.03787626215502,
  84.64026715964413,
  28.255436751363685,
  31.041271486414434,
  26.03138330688587,
  30.69018505212393,
  17.33945898218055,
  20.741856365916416,
  34.39656128436948,
  32.25041486248019,
  17.219253544890268,
  12.162679834210863,
  76.09708704849265,
  47.073834808572336,
  72.099192350003,
  60.254810665237564,
  33.57402852192815,
  11.576817532025064,
  19.85440546798221,
  67.56277601313558,
  25.064635587018305,
  15.443547837581601,
  65.57509063286506,
  22.758926589522268,
  105.46137886924613,
  37.55674567362179,
  100.4918702667481,
  28.22609170658376,
  23.58628829630893,
  21.014816762304413,
  119.55624673828223,
  120.86386619604761,
  32.18507175160181,
  171.79164374499774,
  60.88162981203118,
  15.342634910298637,
  31.16039991158481,
  35.260241571303915,
  27.344143147044896,
  43.74920005398438,
  11.395248413345847,
  117.22010438751501,
  7.8912403139854925,
  14.92401327293268,
  5.393889840367488,
  8.279303087768238,
  59.922764827281384,
  56.19924026232134,
  43.09742253336949,
  7.920815783515331,
  6.312291357021164,
  42.67192747405723,
  58.365192078672614,
  38.32677803151704,
  5.300337379615727,
  42.7261427461295,
  22.148412869993265,
  36.16459499434433,
  69.33811249279483,
  9.746703484072992,
  16.57628123928538,
  12.012023985878013,
  8.314503455211783,
  50.11866377685668,
  74.12072579375187,
  16.836359810931757,
  82.43930612589126,
  42.20531373361757,
  10.323610651744039,
  8.731902879751532,
  63.06134191757923,
  14.053730547292425,
  9.859008783918828,
  122.35038950135372,
  34.984785318756984,
  5.175787148806842,
  26.90752023822015,
  14.966923011582885,
  9.694934265212874,
  12.132592897577792,
  63.48480586647386,
  14.488540589712892,
  86.33819451925297,
  19.916519037233,
  14.054025129256633,
  30.713518193495936,
  59.61974300670012,
  15.655726451629446,
  27.127994680681066,
  19.44611177009812,
  20.74624220488592,
  17.273493357238173,
  29.407363102571487,
  41.41151004164073,
  58.922833632484,
  16.99623868604042,
  111.67970585870904,
  6.6337524708251125,
  72.15417495616964,
  56.763741009605944,
  19.466336565440052,
  5.4586903667915525,
  13.42194879163562,
  36.429114896919046,
  39.90569003842483,
  58.50488771446582,
  8.140508898336048,
  77.11988440481217,
  66.49638457088298,
  66.1187092494653,
  37.1842015976012,
  27.733046819811918,
  25.062580631470595,
  55.668730445523884,
  9.39980948752519,
  49.902978072151924,
  5.339364413536215,
  28.698457873769797,
  41.09262571917477,
  20.785397823265516,
  24.20355309785876,
  14.812128233324046,
  8.817184413313417,
  9.00205739034952,
  14.986275771160003,
  12.623494804779385,
  57.01322277354186,
  55.52008643593597,
  48.865659551273254,
  20.131807390406497,
  20.749650525036113,
  5.026300861271223,
  8.45599248133933,
  92.45525209630476,
  37.667431354283536,
  143.7688430488426,
  18.751578768346857,
  10.700749338077976,
  30.781118280975498,
  14.192862165299847,
  93.09828835541049,
  20.521512926586272,
  51.44049529138689,
  22.63607317752428,
  11.113469436106252,
  62.12643074009492,
  20.954039679420067,
  34.070810908660675,
  18.99759175034995,
  7.857577025138012,
  43.47849580956395,
  63.89474244831867,
  38.90770493838187,
  77.43728837141978,
  5.459167752419965,
  127.29906349863414,
  5.383095064172092,
  7.693472978290682,
  11.082827198674643,
  68.0748354022724,
  57.21124166674078,
  44.14997484735749,
  5.602210525198436,
  19.275917747106114,
  41.09196018570641,
  22.87034700814646,
  42.44229379230722,
  69.71207864950152,
  21.93270739540813,
  67.79745487091058,
  46.26751252021307,
  68.89592531627014,
  90.61133318937571,
  21.41220584214014,
  27.256701390250292,
  36.30318168645079,
  79.34608610252965,
  9.044447224935139,
  5.950268487147574,
  40.38240431743288,
  60.89803801676996,
  28.10544854767372,
  74.76619215455344,
  50.86775087383714,
  85.08042335151902,
  10.830533182773305,
  23.115580720485625,
  31.723429041988723,
  5.716196052743241,
  8.96571407769283,
  70.34579197806306,
  8.966016800407095,
  10.030751580815174,
  47.5688196455782,
  13.181344835298118,
  44.28130498181924,
  7.522916290595509,
  7.555249636135615,
  102.01260700723242,
  48.192499493770384,
  27.044265862519545,
  25.85161728501255,
  38.873623105381476,
  81.51277375987193,
  87.99466326315533,
  31.240924251921584,
  13.329046932719358,
  9.404664704760718,
  14.435792218680437,
  72.0443409836828,
  7.794475927166077,
  77.56027135258952,
  30.830992614557633,
  25.495027606236977,
  6.917058026196925,
  130.56705327642095,
  26.983162966713323,
  31.667837752358892,
  29.600001058099423,
  5.714349567598257,
  6.904173212543638,
  74.86827333964791,
  27.932664063219345,
  28.85069525791109,
  49.38916189379995,
  69.25670652391995,
  62.779149276582444,
  52.595044306505464,
  89.15607907736654,
  9.34721153441597,
  24.34719660791221,
  20.41448328978249,
  16.04441135592343,
  10.433584803268053,
  41.329080110011425,
  9.058154687872289,
  107.64005027072655,
  21.476287962797628,
  27.566332335837902,
  17.96942457506728,
  58.09634975946812,
  17.05141931783922,
  39.30135963240701,
  26.150166260130817,
  11.777513349536502,
  21.855166173760516,
  5.504549418459996,
  13.189733326364223,
  5.86941663139027,
  16.543155521487638,
  19.524739642591005,
  8.615995155354518,
  5.209876062915035,
  20.88846153496624,
  29.353452219197568,
  38.6339454632909,
  14.86216835735905,
  28.08579087417601,
  15.660489209010343,
  11.498195620167337,
  11.654182408129442,
  7.687054459033075,
  123.77114121959707,
  15.601458594924091,
  6.970487333764528,
  7.885902320216253,
  5.068858326898298,
  88.46459860141735,
  16.206535031841753,
  28.775431676906226,
  12.767356441026058,
  51.7642611240329,
  12.593596633462976,
  20.907895896471196,
  25.17439843477087,
  25.372131542085338,
  6.666500213041673,
  20.486355234984117,
  70.74549648886186,
  71.72758345826381,
  5.584017223572758,
  15.88439821731767,
  9.327397335839784,
  6.422734138365126,
  65.84960692681398,
  27.100327724312436,
  129.7550856984676,
  17.163285132520386,
  62.78321675176749,
  13.5993367452964,
  43.29506660258941,
  120.87570479279883,
  59.97229167923042,
  14.144739312122685,
  61.08029609878962,
  17.64487537239363,
  13.968588444332157,
  38.30842664824547,
  8.76121376048288,
  30.38382556535369,
  19.355231618011253,
  31.61707215931181,
  101.70973825363468,
  25.493149793521606,
  29.778798011146367,
  43.53481729816229,
  10.31811299822585,
  69.4129713158506,
  14.166196728180852,
  25.196826886038085,
  33.36760750559982,
  14.025261656043163,
  16.14768328979134,
  32.477808164586584,
  74.26724377052109,
  15.434845843483465,
  8.183968653324403,
  26.538158010686168,
  6.879197658658789,
  60.398272696129474,
  5.310353242136236,
  62.89231453628379,
  24.24569009623564,
  9.603213525856418,
  8.365369051724706,
  14.595376930320295,
  7.180325947575875,
  13.569150077236507,
  13.548306605005875,
  13.756401744802467,
  25.686402971722902,
  41.400325223071654,
  6.345104632887465,
  7.617117158524453,
  100.20620548845756,
  46.074423169882195,
  6.060119770741594,
  50.729143167703825,
  28.291206832572946,
  10.612986744057762,
  76.65606729944457,
  20.066558174302738,
  15.070747671027846,
  17.471548155181615,
  52.7980943118385,
  18.791149765462272,
  6.145565298520987,
  13.060010486809109,
  14.139596835161989,
  32.43691765651124,
  122.06776761250734,
  8.639387517868306,
  85.18182413679119,
  42.787166476670336,
  36.11863781779974,
  7.924546047289133,
  32.740198027067784,
  19.378199770780938,
  32.69629215008728,
  23.530498211989027,
  140.55604676099674,
  18.35643278368189,
  5.168840362449723,
  11.641880772650373,
  24.931413911878312,
  37.8375095301648,
  36.51837841597658,
  7.214595175924872,
  89.05749480457942,
  77.16388397722504,
  40.043820630679434,
  20.037731511293774,
  53.447818244001375,
  25.590491600348273,
  10.364150230970486,
  55.569553515506676,
  6.010412243570199,
  18.929854136344957,
  5.100688294790798,
  54.429358634056,
  5.162597844839081,
  6.803363126210574,
  43.24178418444137,
  53.158592393886494])

Count of hits during the discharge

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);

([5578977.992637578,
  5613539.99644256,
  5621984.067356959,
  5625185.491974772,
  5630368.1422281265,
  5632342.76281048,
  5700317.071617094,
  5743468.085524118,
  5744836.701108547,
  5774398.334295918,
  5774767.373556603,
  5774883.6945538055,
  5777203.405797355,
  5793085.617795732,
  5821571.743378873,
  5855827.613748997,
  5925625.695451189,
  5933944.544405776,
  5934801.5625,
  5940192.840524625,
  5946482.516304447,
  5969773.381161929,
  5981136.763201982,
  5984070.243328004,
  6013602.588212368,
  6013635.187115682,
  6031473.283657751,
  6037470.24548906,
  6042843.5027544135,
  6064294.474320638,
  6094475.0,
  6094905.9829924675,
  6095592.947819843,
  6095976.71663147,
  6096367.1875,
  6100604.671760609,
  6127496.660020787,
  6127514.66433957,
  6131121.586458964,
  6135403.121850532,
  6136735.932735263,
  6141650.962756152,
  6142076.310749517,
  6145591.865585281,
  6145618.491309015,
  6149554.421623171,
  6150099.041201457,
  6151184.375,
  6153137.020367047,
  6155257.734810135,
  6159103.125,
  6159982.962820449,
  6164579.583200394,
  6165053.125,
  6165408.029470348,
  6167709.375,
  6175208.476366167,
  6182233.18007594,
  6184032.704570248,
  6185279.554791648,
  6185349.042431714,
  6187329.873797793,
  6187452.158667843,
  6191708.011044823,
  6192119.60121315,
  6192157.165202618,
  6199457.507883905,
  6218398.34190127,
  6218690.570040212,
  6218746.875,
  6220405.099065951,
  6220965.103675448,
  6222523.60000546,
  6223349.563405839,
  6224387.228177192,
  6224443.575621733,
  6226151.5625,
  6228858.775595083,
  6230673.432989561,
  6253256.072526208,
  6265052.452044818,
  6274830.950682911,
  6277836.20900789,
  6314249.280928683,
  6315299.579220695,
  6315736.783741842,
  6316443.123659642,
  6316910.9375,
  6317419.647893409,
  6318677.500032619,
  6319442.44569348,
  6319450.113932871,
  6319712.830529361,
  6320429.6875,
  6322545.811595063,
  6323816.559097993,
  6323824.821548586,
  6331484.375,
  6355723.014527185,
  6357360.588598309,
  6363796.557840736,
  6395125.521991844,
  6395641.637876963,
  6400400.543865813,
  6401275.067899309,
  6423478.595334258,
  6426760.9375,
  6434744.811132539,
  6442036.355002249,
  6444163.76806014,
  6452397.865784013,
  6462296.591062106,
  6640179.5698098745,
  6924175.0,
  6932273.239346904,
  6978174.583763218,
  6985594.689982541,
  6986289.0625,
  7009900.567918236,
  7049085.756867959,
  7050059.154200402,
  7053084.837842204,
  7054303.048294704,
  7063516.622922211,
  7073588.617017695,
  7200593.640329458,
  7204766.782160621,
  7240366.622016894,
  7317253.268836167,
  7364354.6875,
  7365233.114811757,
  7367649.523492468,
  7493069.116851672,
  7503121.868162566,
  7503184.560158351,
  7531570.1224347735,
  7573246.109365999,
  7614677.878747302,
  7642043.848840267,
  7642056.026815211,
  7740341.077747443,
  7837390.625,
  7871174.495460565,
  7872129.6875,
  7908156.925440548,
  7918716.881316441,
  7919718.632875429,
  7922573.709157177,
  7923678.177498328,
  7924698.879963058,
  7925407.8125,
  7926249.4093552325,
  7927146.96573762,
  7931012.345965633,
  7939362.115827618,
  7941796.346337401,
  7942101.001009151,
  7943187.988806127,
  7946514.697059679,
  7948926.203486298,
  7949316.016096984,
  7977941.478919659,
  7997917.289809064,
  7998522.693319561,
  8001081.420265849,
  8004648.631118895,
  8004875.722019136,
  8005660.933549159,
  8006579.534724715,
  8006725.0,
  8006839.0625,
  8006954.982699717,
  8008182.244234026,
  8010172.892383554,
  8025226.5625,
  8027709.559995733,
  8027712.273284319,
  8028437.681708809,
  8028718.75,
  8028848.4375,
  8029351.5625,
  8030138.62628424,
  8030164.0625,
  8030282.884151183,
  8030343.75,
  8030579.296380469,
  8030788.566129072,
  8031312.5,
  8032257.8125,
  8039522.919105169,
  8045755.033987068,
  8045815.030558283,
  8051632.177019472,
  8055023.17636913,
  8055675.966119183,
  8055677.96209952,
  8058149.7678323435,
  8058416.82410652,
  8059667.893819579,
  8059670.3125,
  8063958.710953828,
  8072553.726499772,
  8094066.422983345,
  8129896.6949967425,
  8141425.052780143,
  8142815.625,
  8158021.623618728,
  8167913.766938603,
  8174902.733161797,
  8181806.800345287,
  8187138.001237432,
  8189534.981143247,
  8189853.418130742,
  8190473.4375,
  8190840.761745352,
  8257782.8125,
  8262367.272545,
  8263788.017466279,
  8268400.745466446,
  8268609.375,
  8268650.712523266,
  8271246.875,
  8274221.875,
  8274650.0,
  8274695.3125,
  8274698.026401529,
  8274734.146565861,
  8274752.336769589,
  8274771.139793693,
  8274815.314464465,
  8274909.577474334,
  8276568.75,
  8277005.824703687,
  8278117.1875,
  8292855.987967529,
  8329567.989503332,
  8369545.9785145745,
  8370513.431092635,
  8371887.753969119,
  8374074.405507117,
  8374076.643530111,
  8382171.457617436,
  8422064.237096388,
  8422096.821402773,
  8448323.72457198,
  8448454.343248893,
  8452726.160929244,
  8454011.583374472,
  8454603.548884489,
  8454943.75,
  8455778.728842134,
  8482536.891196294,
  8482783.0513977,
  8484149.487759253,
  8485708.153880486,
  8487173.382844176,
  8488025.0,
  8488662.5,
  8495114.019980155,
  8495364.513272393,
  8495590.435909467,
  8496115.522262499,
  8502990.527542615,
  8505265.625,
  8506783.675371563,
  8506812.495052168,
  8509096.150287287,
  8510243.75,
  8510427.786932219,
  8511181.25,
  8512575.838905351,
  8513537.51068289,
  8518229.666317971,
  8518587.5,
  8518845.066480987,
  8538177.352235738,
  8538651.698696647,
  8541294.02065175,
  8550751.333775036,
  8550788.678866152,
  8552459.256242624,
  8552723.455720842,
  8553490.139605707,
  8563675.928755844,
  8568046.706448179,
  8569368.75,
  8570442.8937564,
  8572202.725852901,
  8576441.562753767,
  8578176.267249139,
  8579606.241359796,
  8580919.879156416,
  8583482.621314054,
  8585365.625,
  8585791.663102852,
  8585818.44421204,
  8586276.206291188,
  8586568.75,
  8589244.446631467,
  8589526.5625,
  8591485.654464064,
  8600761.594623223,
  8606515.061278822,
  8610595.910508148,
  8622675.934620805,
  8625634.908882054,
  8625814.485315993,
  8626881.52629706,
  8633973.971662201,
  8639836.424902648,
  8664127.054110654,
  8664489.716332259,
  8664549.756569194,
  8664722.522404205,
  8665114.948647423,
  8665252.087675111,
  8665749.782820417,
  8666228.125,
  8666942.1875,
  8667709.605118735,
  8668166.404484307,
  8668435.85852967,
  8669265.625,
  8669269.705527764,
  8669285.9375,
  8670121.88575467,
  8670287.5,
  8670294.517241456,
  8670517.89229598,
  8671609.375,
  8672554.104710143,
  8672950.339512775,
  8673110.460658483,
  8673431.25,
  8674107.46158853,
  8674760.951248473,
  8675641.758034034,
  8677970.237481065,
  8678623.8448693,
  8683148.323229674,
  8691526.043281795,
  8707683.95163068,
  8713364.0625,
  8725007.569254583,
  8725190.625,
  8725530.76560904,
  8725793.054255715,
  8727321.962750567,
  8728206.780048713,
  8735364.953771826,
  8737135.380865166,
  8741326.743886946,
  8752984.99766059,
  8754162.15317943,
  8761408.526926685,
  8792685.360675436,
  8792758.041077849,
  8793653.832615513,
  8793996.875,
  8794022.188779937,
  8794154.6875,
  8795154.6875,
  8795201.5625,
  8797156.25,
  8798655.251237357,
  8802349.879995534,
  8802656.63022709,
  8802989.522512034,
  8803828.891342113,
  8804015.738557713,
  8804023.4375,
  8804145.943872036,
  8804186.860842645,
  8804192.780479454,
  8804219.813627753,
  8804569.712109134,
  8806134.96419869,
  8806203.889612991,
  8806254.6875,
  8806975.906851102,
  8807864.134671347,
  8808102.117476793,
  8919282.607320037,
  8927556.573354574,
  8953757.731735593,
  8954361.746108068,
  8954362.31208534,
  8954415.39249643,
  8955966.60381047,
  8958651.5625,
  8962332.555338075,
  8971112.321127882,
  8977868.233932119,
  8980165.436467165,
  8983996.92816116,
  8984599.17041027,
  8985237.5,
  8985248.935334267,
  8985299.896004494,
  8988940.625,
  8989925.0,
  8989927.606602695,
  8991127.35511087,
  8991185.500495778,
  8993061.513141083,
  8996273.120763056,
  9004222.812932868,
  9007014.95784737,
  9028259.238826694,
  9030218.242034517,
  9032344.708575925,
  9035695.65798077,
  9036945.3125,
  9037044.923101418,
  9037098.69224798,
  9038609.375,
  9039574.616101801,
  9039973.000686746,
  9040032.284541702,
  9041434.707256664,
  9041732.359365825,
  9042234.361240216,
  9042501.5625,
  9042831.284776246,
  9042887.5,
  9043064.575609773,
  9043202.546242805,
  9043281.25,
  9043592.066514777,
  9044362.991936032,
  9044486.512079779,
  9044854.017704967,
  9045270.3125,
  9045492.95490584,
  9045494.87490307,
  9045720.127382657,
  9045908.082004739,
  9046732.13088164,
  9046779.482426535,
  9048838.858840046,
  9048903.789355332,
  9049677.784857232,
  9052394.930674184,
  9053596.875,
  9053828.930609608,
  9062087.914366238,
  9067817.172224868,
  9070774.077437205,
  9107027.744620299,
  9109396.695817403,
  9129387.419966767,
  9129628.125,
  9133234.420686716,
  9137449.976875912,
  9157097.79604507,
  9174989.89398588,
  9178357.8125,
  9180308.896938387,
  9205695.091304233,
  9227506.414677413,
  9233628.954627037,
  9279602.294272551,
  9289142.182902094,
  9289290.504566358,
  9299359.361218933,
  9303648.4375,
  9307773.782939412,
  9309326.5625,
  9309650.800261535,
  9349581.598764736,
  9351163.369059006,
  9358130.934509596,
  9360859.375,
  9361247.1715868,
  9365434.482828114,
  9386427.518439142,
  9387284.375,
  9387340.625,
  9387694.203454724,
  9387923.469339475,
  9388744.119543774,
  9389051.614938,
  9389479.370882215,
  9389733.16428069,
  9390684.004603539,
  9391330.75013092,
  9411608.573555443,
  9412030.666474579,
  9412792.944299256,
  9412845.3125,
  9412851.308301268,
  9413583.376862906,
  9415168.94505935,
  9415216.51633216,
  9415425.442093844,
  9415459.87429085,
  9415565.631651614,
  9415673.838264123,
  9415698.997529536,
  9415723.4375,
  9415790.500855843,
  9415873.614106368,
  9416124.885779316,
  9416262.24589332,
  9416313.714787347,
  9416316.916951617,
  9416423.613059886,
  9416474.545621373,
  9416500.829499535,
  9416602.546723524,
  9416828.932083739,
  9417016.571739135,
  9417062.328590088,
  9418001.5625,
  9418397.90640377,
  9418610.835803391,
  9418635.250859113,
  9418748.08130043,
  9432352.294239253,
  9550716.457110025,
  9552684.303133195,
  9553108.909647174,
  9654869.71805141,
  9745972.988154227,
  9748710.47078961,
  9764723.150676485,
  9765188.735763712,
  9791240.625,
  9798694.121990992,
  9799732.710472522,
  9800043.708696079,
  9802784.171198906,
  9802805.990969595,
  9803242.267307982,
  9850576.5625,
  9850640.550180549,
  9872050.018251656,
  9905679.190650804,
  9909968.176003896,
  9911297.66728433,
  9913121.840636093,
  9916705.572066901,
  9948380.833003595,
  9965781.261695607,
  9972839.846731849,
  9973690.479978068,
  9973754.6875,
  9974052.746241923,
  9974161.851776069,
  9974490.047096923,
  9974545.950481862,
  9974584.68219173,
  9974604.730765797,
  9974646.83061326,
  9974682.513091875,
  9975133.547908572,
  9975212.208521167,
  9975220.735421192,
  9975229.476164848,
  9975235.811844306,
  9975765.021756265,
  9975856.501909183,
  9976007.419921272,
  9976057.80337224,
  9976065.266040968,
  9976886.342057634,
  9977187.243001143,
  9977210.144756453,
  9977274.089917224,
  9977470.3125,
  9977504.185187697,
  9978905.05714484,
  9979243.707714986,
  9982509.31431189,
  9985512.5,
  9992113.654304542,
  9993020.555281322,
  9996762.471096417,
  9996787.381402094,
  10002379.336804705,
  10024478.125,
  10027391.971597543,
  10035632.623142052,
  10037956.857204715,
  10038008.361725478,
  10040890.37532692,
  10041982.8125,
  10065853.38125877,
  10091403.838716848,
  10092895.3125,
  10092963.321090968,
  10095299.720654031,
  10097807.304352304,
  10097856.253012737,
  10098667.1875,
  10099229.6875,
  10099909.014946127,
  10100182.936398685,
  10100247.188947475,
  10100470.178203166,
  10101104.057642618,
  10101162.825340625,
  10101854.6875,
  10105054.78391054,
  10131292.45049253,
  10161616.416436505,
  10161626.269829186,
  10161729.834783094,
  10169688.595853087,
  10170138.014750518,
  10179129.221269872,
  10179714.788008817,
  10184406.21191207,
  10188168.196301162,
  10188721.875,
  10190524.772325411,
  10191598.302842084,
  10191712.078098169,
  10199184.375,
  10199235.857941551,
  10208423.172623368,
  10208909.375,
  10210021.875,
  10211332.8125,
  10226490.625,
  10232635.438595347,
  10235212.585390566,
  10237767.544899004,
  10270728.125,
  10283002.062005082,
  10290050.0,
  10292251.5625,
  10311263.88779654,
  10386710.098369364,
  10391426.332413316,
  10466773.854691127,
  10477916.968040938,
  10480777.584047493,
  10490746.993626973,
  10503769.415995158,
  10520482.016492838,
  10541145.250601705,
  10541158.407772304,
  10548296.875,
  10573725.0,
  10773496.875,
  10774728.905083535,
  10774751.5625,
  10776792.206136825,
  10777665.020240204,
  10777668.13461999,
  10778846.71205643,
  10780182.317867313,
  10780895.3125,
  10781879.34564069,
  10782364.0625,
  10782616.685812997,
  10782751.512144936,
  10784220.47316309,
  10785047.461768515,
  10785567.649202945,
  10785963.83995435,
  10786879.554680398,
  10789251.242923414,
  10790342.008619381,
  10791541.585273359,
  10794124.449462075,
  10794867.909633797,
  10794883.177574959,
  10794896.152109306,
  10798670.88748433,
  10798675.251460712,
  10800201.774210311,
  10800576.472235681,
  10800866.467850756,
  10801371.854027035,
  10801583.580403736,
  10801655.203027135,
  10804001.787265554,
  10804373.7500902,
  10804408.038444053,
  10804666.718057577,
  10805223.262449266,
  10805324.392792817,
  10806161.656789578,
  10806224.997911619,
  10807836.321072653,
  10808036.644279381,
  10808415.625,
  10808581.873706874,
  10808873.4375,
  10808939.833839264,
  10808980.868800774,
  10809117.556962334,
  10809218.567149037,
  10809468.243870266,
  10809663.455287416,
  10809883.819216076,
  10810860.650136447,
  10811066.631674938,
  10811229.246688634,
  10811451.55775868,
  10812146.451787293,
  10812288.74001905,
  10814082.73614174,
  10814429.6875,
  10814470.794442,
  10814634.175857754,
  10814822.984567305,
  10815081.21234314,
  10815575.934037207,
  10816437.5,
  10817226.336476753,
  10817231.25,
  10817241.190187208,
  10817316.994471196,
  10817685.920652514,
  10817708.61376789,
  10818650.44962985,
  10818785.0410984,
  10818949.73800539,
  10819542.563609054,
  10819618.57517276,
  10820283.905993085,
  10821180.385645043,
  10821563.474445011,
  10821962.366196258,
  10822690.405543769,
  10823060.742860926,
  10824867.355729358,
  10825224.736878874,
  10825530.41632611,
  10825575.0,
  10825796.616920745,
  10825852.074789431,
  10826027.925790953,
  10826181.140819779,
  10826309.928126782,
  10826581.25,
  10826678.736373166,
  10827726.356233587,
  10828785.467177365,
  10828997.179867625,
  10829662.646082487,
  10829965.531124143,
  10830422.094331188,
  10830576.343721554,
  10831079.266695555,
  10831707.462456234,
  10831860.775297178,
  10832075.0,
  10832606.842710985,
  10833504.725712355,
  10834208.209129078,
  10834461.755489081,
  10835198.16407682,
  10835381.25,
  10835405.036137272,
  10835788.706193324,
  10836144.412098594,
  10836521.80788944,
  10836907.84599894,
  10836979.299859766,
  10837014.0625,
  10837414.803512307,
  10838227.097552575,
  10838352.524900647,
  10838645.084898634,
  10838932.537461119,
  10839782.707676722,
  10839966.237466928,
  10840163.452415587,
  10840168.75,
  10840573.001947701,
  10840835.183819659,
  10841051.5625,
  10841305.140630526,
  10841315.26370543,
  10841873.732887015,
  10841873.955943927,
  10842635.052716026,
  10843007.96124787,
  10845167.493300281,
  10846154.474631691,
  10846417.784582015,
  10846660.9375,
  10850361.280366791,
  10851604.490230924,
  10851920.3125,
  10852854.71921794,
  10852988.403083177,
  10854441.81241589,
  10854476.873962175,
  10854951.012120415,
  10855042.432351183,
  10855117.968278056,
  10856024.359284345,
  10856636.050576696,
  10856885.9375,
  10857321.755933389,
  10858837.696414525,
  10859606.116575606,
  445897168.75],
 [10.926932224398904,
  20.299863187360607,
  128.03679543984842,
  22.054290030078057,
  92.0134442811754,
  8.616197603256872,
  6.076130750777888,
  99.513277837992,
  62.942603518181855,
  14.562212679960686,
  69.58291917703299,
  7.9567872446315615,
  6.820890981208564,
  67.7001006400315,
  6.73751087261893,
  117.25589082591692,
  186.77098044642872,
  15.759888942893621,
  111.42139242152703,
  209.23843133454497,
  8.135918106382611,
  72.9941144126226,
  25.883289989482332,
  36.68592379563776,
  64.15916495810745,
  53.54018954297891,
  10.895464639696927,
  26.62252597598342,
  72.73778335531173,
  8.7107392510694,
  49.865399405062874,
  112.20082055767892,
  82.64480473827645,
  6.286870166924756,
  48.66363464808074,
  17.15804067148037,
  9.156008998167568,
  48.4507440634721,
  59.66783805117294,
  28.462945030734435,
  28.43112284737409,
  13.060927451972757,
  24.3872693557942,
  42.247653860659256,
  6.061604517463069,
  24.196628613696685,
  12.223600475285949,
  87.33087495537146,
  112.04535780831122,
  99.35487052097429,
  34.3612146902758,
  69.82759419328755,
  93.7767087179709,
  80.87655749057232,
  65.84386632806763,
  46.10461411743419,
  95.29772568023517,
  54.90189286240101,
  90.08621857118047,
  26.15196044445469,
  5.422114749261968,
  20.46881669574627,
  9.892079739524835,
  64.04483853197218,
  19.616002262174746,
  22.426241203402252,
  23.69195598151165,
  26.767377412322464,
  27.48242570625202,
  78.17984072236509,
  16.782526933300517,
  5.776388073977832,
  17.16245783712184,
  22.14981813755204,
  41.21862258940446,
  32.30695594607249,
  52.61387621741336,
  16.04861373266941,
  55.8379853597047,
  25.453079411126218,
  78.58967712147752,
  23.759468818820626,
  6.116796188298128,
  59.45093622154312,
  6.682311490039749,
  55.73381667475616,
  65.93507541908905,
  74.76998507261524,
  7.963836959850351,
  12.996342339540242,
  115.28563924400606,
  11.310586744054634,
  43.415073325361774,
  89.11735884359675,
  7.794309767515815,
  8.851141962612628,
  14.302330732534013,
  36.26773275395881,
  130.4621411981585,
  23.145561783174436,
  139.05611076736474,
  17.39055387202637,
  13.192164245260043,
  21.52157689636502,
  44.57887992328304,
  137.21496266899408,
  29.43156867413411,
  21.479469924041624,
  11.849662875070816,
  78.45085521222896,
  78.00925312492446,
  64.21466078080728,
  25.522473652892575,
  44.0564159458759,
  7.070319545796551,
  22.37146180037758,
  21.563094880394516,
  107.55994206406913,
  145.031435332794,
  14.350142723374486,
  61.36056062936864,
  97.6706598756864,
  36.30821344342823,
  5.229613973701861,
  95.05866320511922,
  10.960505086983037,
  10.54819837019146,
  10.34560973780322,
  29.46945492806583,
  102.24459712685095,
  41.12275630692974,
  6.6719635584127674,
  110.69055943650886,
  228.96045833314355,
  6.51651349414073,
  10.843073802257575,
  10.110261060396558,
  7.516536059594127,
  7.147678427959863,
  11.498108948890188,
  77.84171886428932,
  58.28900402519001,
  39.79299893595178,
  96.20372095280756,
  69.1724976588664,
  23.750907548738354,
  6.252513071890919,
  44.64535189066154,
  9.376191502703486,
  34.17020809979182,
  39.8741001375511,
  20.63414102952904,
  8.335150925109934,
  25.813337165548567,
  10.576002151515361,
  5.324126930341678,
  97.83932753890085,
  24.294792464321212,
  93.7812136750128,
  39.116434751511186,
  16.675990078357046,
  88.01702599711518,
  32.66879312901625,
  89.10860079246407,
  27.594867061153096,
  31.01760277897975,
  7.887068169157723,
  14.92647202087665,
  15.615817247997839,
  51.789081444436036,
  33.959268727929306,
  22.91985033134695,
  51.64026536488472,
  8.01791625522131,
  35.40306319928818,
  156.88821685805033,
  9.146787521338487,
  6.290432215728766,
  125.86227009929267,
  32.06036861177417,
  53.00435486744217,
  6.043881677488489,
  42.60987894384768,
  120.96697446291563,
  53.56839170152092,
  5.445225985469572,
  10.431395613628844,
  58.17035926230734,
  42.99555130513441,
  6.6642455857204475,
  11.168018450438,
  25.532552010210814,
  19.689447108991835,
  17.524387391947783,
  19.55212414753965,
  7.079278505288826,
  24.657682728096812,
  29.22295533164956,
  6.003370524791343,
  28.59439664529859,
  153.8087198798159,
  17.86646735726233,
  100.59607538761418,
  136.46283028327977,
  6.273889867406289,
  183.8744551224653,
  43.926633821223035,
  23.807604716859107,
  52.572638966366405,
  107.54586394540995,
  53.073939969998364,
  40.95730382658491,
  14.297051575578589,
  36.705622648118975,
  15.489658484792963,
  45.544627970500144,
  108.25080523661097,
  12.045894086567241,
  18.82275814450937,
  40.447185168777445,
  72.69234609727621,
  28.957450863253055,
  89.09193763462116,
  69.26511707333464,
  52.12826840805946,
  105.97248009747663,
  24.71006773504476,
  6.615724045093962,
  38.95273399486335,
  23.617091228379515,
  5.542988620937467,
  31.7991050186689,
  8.995942330408488,
  176.5495347256891,
  6.234650523586943,
  19.236323271931933,
  5.312113059171831,
  12.980005640377659,
  140.47715832597345,
  109.97914872695404,
  25.48527871948557,
  27.884424147593194,
  6.289512289891066,
  50.08977361861045,
  11.608163153690981,
  26.462111074065266,
  13.927865004684602,
  7.3022543197622864,
  8.524219738475079,
  52.07943080778101,
  21.279786557519145,
  17.66432269755832,
  11.522845463577987,
  8.93283532971391,
  34.50895901041959,
  39.80426715013291,
  62.17654697789925,
  49.88482558293417,
  6.08368866025267,
  16.991811347639135,
  9.176475566574483,
  37.33761558860229,
  74.97984286531754,
  48.88550122841634,
  10.185810200964234,
  13.711900607875883,
  10.156867299279295,
  29.16122290640686,
  13.79255490299566,
  68.50916047843575,
  12.752727222747184,
  11.144446005659084,
  7.115894656315024,
  32.209759154141395,
  50.18522031979376,
  19.144592673404816,
  33.926451187125316,
  84.67245875804198,
  22.2356207471094,
  10.576712785890933,
  14.50739094047482,
  9.346857592635573,
  21.640047878045056,
  15.956878803880418,
  6.246135368968968,
  47.482792030189735,
  328.98981815647835,
  19.261139174294765,
  45.33088099576746,
  55.933661687859725,
  11.114228412551258,
  13.512386339955809,
  25.24787608992044,
  59.14961944476913,
  22.905120363684794,
  13.891397854826806,
  65.3949650866028,
  83.27717189883589,
  158.9455942157528,
  56.2937382677257,
  19.606899925577643,
  82.55026200982032,
  20.880156925777406,
  6.152575774114602,
  52.87620633233132,
  12.079102617513728,
  20.505580177970533,
  16.82753595845946,
  117.21772669193572,
  30.173683254735415,
  17.217423679258687,
  7.858171713084091,
  14.08831159433832,
  5.565348336016015,
  12.552443839528934,
  25.691355953270538,
  53.39862287539856,
  39.94515272603766,
  50.77449085279253,
  5.251631814585299,
  51.36063314871676,
  27.054548006621,
  42.93135107541445,
  13.042320824040967,
  38.94119007439107,
  11.144558754842182,
  92.68322931856662,
  18.97929379309364,
  6.912872776632882,
  37.33817091681863,
  57.29560793633324,
  6.541663270948435,
  8.9583217551817,
  74.32972366732852,
  14.140485442797875,
  33.78701765844463,
  19.717480382119135,
  27.123872689534732,
  61.7588220920513,
  14.523146327413878,
  23.922619236500225,
  10.524709933429127,
  44.506801835712075,
  9.133859605102945,
  41.774104081097335,
  15.637719791498562,
  7.95000450316414,
  20.683192917795644,
  46.471162797043924,
  68.29503440851497,
  19.331465650822423,
  80.75747116634749,
  52.67749964093652,
  5.595141066061806,
  8.70433261905094,
  21.941461113579024,
  55.64705553365039,
  5.753322795152172,
  50.88161940674756,
  93.53654113240425,
  46.27715075107495,
  47.241021412927665,
  36.17471792541111,
  52.43536290417473,
  36.81032822712951,
  5.929022504706373,
  16.620152567450507,
  17.042161528676985,
  5.147225244121349,
  15.392754998128405,
  60.542797456236634,
  18.05436732871967,
  37.33168933760188,
  5.859198364572408,
  18.644646682734045,
  20.546879513972723,
  37.8523181725523,
  13.023912759924865,
  35.43682901855933,
  20.069275406330537,
  15.223918110283714,
  26.907083211472383,
  26.481562071959825,
  79.84973667457616,
  27.02321109646978,
  80.06084596224814,
  25.296468013396886,
  5.283959008515075,
  43.99361177737693,
  50.237691467659026,
  52.71513265995131,
  68.36976420458642,
  46.37473661916552,
  25.28726829855669,
  119.89988444841401,
  18.73400665211349,
  35.943513482928225,
  36.210555347685855,
  5.806325352253628,
  85.56712408895677,
  40.5039646438503,
  21.311519652857456,
  12.643042285563753,
  8.987343587627201,
  15.894972587846649,
  119.70250995747693,
  7.391952283345661,
  14.312142661494214,
  38.33810008921233,
  57.22146334905344,
  26.112882090615813,
  29.08788315134973,
  60.32957979414643,
  10.569114193689918,
  78.82062341725094,
  49.20380341159883,
  22.735121178256865,
  27.29899539570856,
  171.60564229959652,
  47.03787626215502,
  84.64026715964413,
  28.255436751363685,
  31.041271486414434,
  26.03138330688587,
  30.69018505212393,
  17.33945898218055,
  20.741856365916416,
  34.39656128436948,
  32.25041486248019,
  17.219253544890268,
  12.162679834210863,
  76.09708704849265,
  47.073834808572336,
  72.099192350003,
  60.254810665237564,
  33.57402852192815,
  11.576817532025064,
  19.85440546798221,
  67.56277601313558,
  25.064635587018305,
  15.443547837581601,
  65.57509063286506,
  22.758926589522268,
  105.46137886924613,
  37.55674567362179,
  100.4918702667481,
  28.22609170658376,
  23.58628829630893,
  21.014816762304413,
  119.55624673828223,
  120.86386619604761,
  32.18507175160181,
  171.79164374499774,
  60.88162981203118,
  15.342634910298637,
  31.16039991158481,
  35.260241571303915,
  27.344143147044896,
  43.74920005398438,
  11.395248413345847,
  117.22010438751501,
  7.8912403139854925,
  14.92401327293268,
  5.393889840367488,
  8.279303087768238,
  59.922764827281384,
  56.19924026232134,
  43.09742253336949,
  7.920815783515331,
  6.312291357021164,
  42.67192747405723,
  58.365192078672614,
  38.32677803151704,
  5.300337379615727,
  42.7261427461295,
  22.148412869993265,
  36.16459499434433,
  69.33811249279483,
  9.746703484072992,
  16.57628123928538,
  12.012023985878013,
  8.314503455211783,
  50.11866377685668,
  74.12072579375187,
  16.836359810931757,
  82.43930612589126,
  42.20531373361757,
  10.323610651744039,
  8.731902879751532,
  63.06134191757923,
  14.053730547292425,
  9.859008783918828,
  122.35038950135372,
  34.984785318756984,
  5.175787148806842,
  26.90752023822015,
  14.966923011582885,
  9.694934265212874,
  12.132592897577792,
  63.48480586647386,
  14.488540589712892,
  86.33819451925297,
  19.916519037233,
  14.054025129256633,
  30.713518193495936,
  59.61974300670012,
  15.655726451629446,
  27.127994680681066,
  19.44611177009812,
  20.74624220488592,
  17.273493357238173,
  29.407363102571487,
  41.41151004164073,
  58.922833632484,
  16.99623868604042,
  111.67970585870904,
  6.6337524708251125,
  72.15417495616964,
  56.763741009605944,
  19.466336565440052,
  5.4586903667915525,
  13.42194879163562,
  36.429114896919046,
  39.90569003842483,
  58.50488771446582,
  8.140508898336048,
  77.11988440481217,
  66.49638457088298,
  66.1187092494653,
  37.1842015976012,
  27.733046819811918,
  25.062580631470595,
  55.668730445523884,
  9.39980948752519,
  49.902978072151924,
  5.339364413536215,
  28.698457873769797,
  41.09262571917477,
  20.785397823265516,
  24.20355309785876,
  14.812128233324046,
  8.817184413313417,
  9.00205739034952,
  14.986275771160003,
  12.623494804779385,
  57.01322277354186,
  55.52008643593597,
  48.865659551273254,
  20.131807390406497,
  20.749650525036113,
  5.026300861271223,
  8.45599248133933,
  92.45525209630476,
  37.667431354283536,
  143.7688430488426,
  18.751578768346857,
  10.700749338077976,
  30.781118280975498,
  14.192862165299847,
  93.09828835541049,
  20.521512926586272,
  51.44049529138689,
  22.63607317752428,
  11.113469436106252,
  62.12643074009492,
  20.954039679420067,
  34.070810908660675,
  18.99759175034995,
  7.857577025138012,
  43.47849580956395,
  63.89474244831867,
  38.90770493838187,
  77.43728837141978,
  5.459167752419965,
  127.29906349863414,
  5.383095064172092,
  7.693472978290682,
  11.082827198674643,
  68.0748354022724,
  57.21124166674078,
  44.14997484735749,
  5.602210525198436,
  19.275917747106114,
  41.09196018570641,
  22.87034700814646,
  42.44229379230722,
  69.71207864950152,
  21.93270739540813,
  67.79745487091058,
  46.26751252021307,
  68.89592531627014,
  90.61133318937571,
  21.41220584214014,
  27.256701390250292,
  36.30318168645079,
  79.34608610252965,
  9.044447224935139,
  5.950268487147574,
  40.38240431743288,
  60.89803801676996,
  28.10544854767372,
  74.76619215455344,
  50.86775087383714,
  85.08042335151902,
  10.830533182773305,
  23.115580720485625,
  31.723429041988723,
  5.716196052743241,
  8.96571407769283,
  70.34579197806306,
  8.966016800407095,
  10.030751580815174,
  47.5688196455782,
  13.181344835298118,
  44.28130498181924,
  7.522916290595509,
  7.555249636135615,
  102.01260700723242,
  48.192499493770384,
  27.044265862519545,
  25.85161728501255,
  38.873623105381476,
  81.51277375987193,
  87.99466326315533,
  31.240924251921584,
  13.329046932719358,
  9.404664704760718,
  14.435792218680437,
  72.0443409836828,
  7.794475927166077,
  77.56027135258952,
  30.830992614557633,
  25.495027606236977,
  6.917058026196925,
  130.56705327642095,
  26.983162966713323,
  31.667837752358892,
  29.600001058099423,
  5.714349567598257,
  6.904173212543638,
  74.86827333964791,
  27.932664063219345,
  28.85069525791109,
  49.38916189379995,
  69.25670652391995,
  62.779149276582444,
  52.595044306505464,
  89.15607907736654,
  9.34721153441597,
  24.34719660791221,
  20.41448328978249,
  16.04441135592343,
  10.433584803268053,
  41.329080110011425,
  9.058154687872289,
  107.64005027072655,
  21.476287962797628,
  27.566332335837902,
  17.96942457506728,
  58.09634975946812,
  17.05141931783922,
  39.30135963240701,
  26.150166260130817,
  11.777513349536502,
  21.855166173760516,
  5.504549418459996,
  13.189733326364223,
  5.86941663139027,
  16.543155521487638,
  19.524739642591005,
  8.615995155354518,
  5.209876062915035,
  20.88846153496624,
  29.353452219197568,
  38.6339454632909,
  14.86216835735905,
  28.08579087417601,
  15.660489209010343,
  11.498195620167337,
  11.654182408129442,
  7.687054459033075,
  123.77114121959707,
  15.601458594924091,
  6.970487333764528,
  7.885902320216253,
  5.068858326898298,
  88.46459860141735,
  16.206535031841753,
  28.775431676906226,
  12.767356441026058,
  51.7642611240329,
  12.593596633462976,
  20.907895896471196,
  25.17439843477087,
  25.372131542085338,
  6.666500213041673,
  20.486355234984117,
  70.74549648886186,
  71.72758345826381,
  5.584017223572758,
  15.88439821731767,
  9.327397335839784,
  6.422734138365126,
  65.84960692681398,
  27.100327724312436,
  129.7550856984676,
  17.163285132520386,
  62.78321675176749,
  13.5993367452964,
  43.29506660258941,
  120.87570479279883,
  59.97229167923042,
  14.144739312122685,
  61.08029609878962,
  17.64487537239363,
  13.968588444332157,
  38.30842664824547,
  8.76121376048288,
  30.38382556535369,
  19.355231618011253,
  31.61707215931181,
  101.70973825363468,
  25.493149793521606,
  29.778798011146367,
  43.53481729816229,
  10.31811299822585,
  69.4129713158506,
  14.166196728180852,
  25.196826886038085,
  33.36760750559982,
  14.025261656043163,
  16.14768328979134,
  32.477808164586584,
  74.26724377052109,
  15.434845843483465,
  8.183968653324403,
  26.538158010686168,
  6.879197658658789,
  60.398272696129474,
  5.310353242136236,
  62.89231453628379,
  24.24569009623564,
  9.603213525856418,
  8.365369051724706,
  14.595376930320295,
  7.180325947575875,
  13.569150077236507,
  13.548306605005875,
  13.756401744802467,
  25.686402971722902,
  41.400325223071654,
  6.345104632887465,
  7.617117158524453,
  100.20620548845756,
  46.074423169882195,
  6.060119770741594,
  50.729143167703825,
  28.291206832572946,
  10.612986744057762,
  76.65606729944457,
  20.066558174302738,
  15.070747671027846,
  17.471548155181615,
  52.7980943118385,
  18.791149765462272,
  6.145565298520987,
  13.060010486809109,
  14.139596835161989,
  32.43691765651124,
  122.06776761250734,
  8.639387517868306,
  85.18182413679119,
  42.787166476670336,
  36.11863781779974,
  7.924546047289133,
  32.740198027067784,
  19.378199770780938,
  32.69629215008728,
  23.530498211989027,
  140.55604676099674,
  18.35643278368189,
  5.168840362449723,
  11.641880772650373,
  24.931413911878312,
  37.8375095301648,
  36.51837841597658,
  7.214595175924872,
  89.05749480457942,
  77.16388397722504,
  40.043820630679434,
  20.037731511293774,
  53.447818244001375,
  25.590491600348273,
  10.364150230970486,
  55.569553515506676,
  6.010412243570199,
  18.929854136344957,
  5.100688294790798,
  54.429358634056,
  5.162597844839081,
  6.803363126210574,
  43.24178418444137,
  53.158592393886494])

Graph for home page

Detected energies during the discharge + Energy spectrum

multiplot(icon_fig, T_int_first,E,xmean,count)

Average size of interactions

size_of_interactions_average(Etot, Ntot,size_interaction)

Energy spectra for parts of doublebreadown discharges

#energy_spectra_doublebreakdown(Etot, T_first, file_energy_spectra_doublebreakdown)