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

([3492620.530665972,
  3582442.0812735655,
  3585565.319624986,
  3593032.2745178808,
  3611177.580869366,
  3614108.073302657,
  3616160.892630776,
  3693601.88804027,
  3732662.441962715,
  3803672.662098161,
  3856413.315365188,
  3924180.8969368823,
  3946479.039064631,
  3947820.0847265646,
  3948395.2593749003,
  3982932.808839554,
  3985380.2526885793,
  4004219.513812185,
  4004237.410135303,
  4004249.4211029704,
  4004315.625,
  4004326.5519919926,
  4032283.745416332,
  4048985.846178447,
  4063593.6444698013,
  4065810.262041953,
  4070489.615134211,
  4071332.619565964,
  4103654.834075178,
  4108865.1873446065,
  4111239.255548117,
  4115083.610033039,
  4140226.5625,
  4142648.6840628483,
  4157467.665319404,
  4158065.6581533193,
  4176712.3629468787,
  4176718.0829887926,
  4183423.4375,
  4183431.5333238444,
  4207255.635936901,
  4230417.674915464,
  4240376.5625,
  4245120.833868845,
  4246034.375,
  4278065.625,
  4278338.187965376,
  4296439.785176602,
  4314606.67847943,
  4337496.875,
  4338624.447477073,
  4344571.082869163,
  4368767.926705537,
  4369570.858679403,
  4370025.0,
  4370050.406730455,
  4370078.468745263,
  4382167.1875,
  4385253.618855921,
  4394957.270319626,
  4396710.300909533,
  4400208.5068542315,
  4401721.317354548,
  4403828.125,
  4404592.014905767,
  4405027.614335727,
  4405306.25,
  4411033.476109745,
  4432089.0625,
  4444100.594728054,
  4444638.828062689,
  4444921.463977997,
  4444945.3125,
  4445785.816980325,
  4453308.859898969,
  4467647.402255139,
  4467762.413638079,
  4507940.901302447,
  4513187.907767356,
  4513339.0625,
  4515787.049710502,
  4528356.492906535,
  4541787.5,
  4544634.458786175,
  4545223.42308966,
  4545577.350766959,
  4547352.408610105,
  4547697.608782832,
  4558034.321257396,
  4558949.32115245,
  4559228.836067757,
  4566832.8125,
  4566995.3125,
  4567151.0378156435,
  4567252.79983913,
  4567267.586596807,
  4567650.0,
  4567653.125,
  4568517.733738784,
  4568700.115373324,
  4568706.25,
  4588297.150567199,
  4588322.324956934,
  4606041.277647058,
  4607407.272407189,
  4609157.94107543,
  4622816.727336632,
  4624307.8125,
  4631048.34565392,
  4638880.898431857,
  4661785.110310686,
  4662986.689381507,
  4692118.131158823,
  4700800.0,
  4708326.0583390845,
  4715356.25,
  4718447.077253479,
  4718652.705737127,
  4723722.554191789,
  4738390.625,
  4746995.3125,
  4753088.111919016,
  4761115.5090038525,
  4761337.427545323,
  4762190.138718305,
  4771840.625,
  4778800.644308123,
  4780726.5625,
  4792370.3125,
  4792425.67511725,
  4793339.0625,
  4793453.402516461,
  4793502.388632653,
  4804800.0,
  4810887.5,
  4817915.314200792,
  4823158.910894705,
  4843160.150324343,
  4850747.3125097,
  4852696.875,
  4870146.727989322,
  4874169.022973611,
  4874746.601720019,
  4885514.8859380605,
  4889250.432388543,
  4900423.4375,
  4913836.337753746,
  4920548.32970231,
  4924215.132579058,
  4936441.705124193,
  4936867.781244053,
  4971395.955267249,
  4997933.071737891,
  5004409.9010975715,
  5006145.683023668,
  5014930.233848478,
  5015036.574655041,
  5016085.9375,
  5041387.378949439,
  5043382.149309411,
  5047707.686714367,
  5059738.900431105,
  5063029.740341798,
  5072741.591584372,
  5074093.582756976,
  5074668.611125128,
  5074720.107495489,
  5083679.904929095,
  5087359.006914243,
  5088588.141981034,
  5090221.2414149875,
  5090597.138631836,
  5090893.75,
  5093431.25,
  5098627.692084908,
  5104219.234343314,
  5105066.562350477,
  5110621.086165099,
  5115426.383465704,
  5144596.4875148125,
  5147997.910924866,
  5149588.509453688,
  5151425.0,
  5168914.855442729,
  5173954.6875,
  5202698.4375,
  5202700.0,
  5202898.840043367,
  5202918.736277276,
  5205207.531347487,
  5215830.12677419,
  5246987.987188586,
  5305675.148945387,
  5340362.5,
  5341082.8125,
  5342370.3125,
  5344550.0,
  5356625.0,
  5360156.113141704,
  5386041.450878613,
  5393534.375,
  5393540.545964344,
  5398086.782300612,
  5419719.49582887,
  5429926.08184323,
  5440662.895200803,
  5455821.875,
  5479009.375,
  5479677.034038522,
  5523718.320207222,
  5561756.25,
  5568504.477091591,
  5575988.719570688,
  5601538.788035155,
  5601559.112705029,
  5614991.1128464,
  5631440.625,
  5633206.119030191,
  5633585.9375,
  5633979.6875,
  5677488.938883729,
  5678219.195279469,
  5682221.5841641305,
  5690881.25,
  5707609.81274805,
  5731804.6875,
  5742124.622823982,
  5774245.3125,
  5777536.395349934,
  5782454.376485919,
  5792843.969078138,
  5797463.898883178,
  5822959.356834526,
  5826571.875,
  5833190.984185385,
  5833355.207311864,
  5835251.343775273,
  5835509.922134748,
  5835565.932978031,
  5875859.590684232,
  5876278.717555765,
  5892808.722669422,
  5903393.025912761,
  5911587.5,
  5956122.610491333,
  5964578.125,
  5981715.560970378,
  5990124.825824275,
  6041239.131424672,
  6066414.365561938,
  6137851.095705542,
  6139587.3817310715,
  6157208.400800219,
  6165635.91745936,
  6190257.4107712675,
  6194263.533878264,
  6195880.800181757,
  6211297.072798556,
  6211332.096785197,
  6211990.438800751,
  6214761.580308128,
  6234113.528688682,
  6238683.020853411,
  6243307.8125,
  6253770.3125,
  6284650.0,
  6289653.125,
  6291201.004408521,
  6291447.750224567,
  6299438.434041375,
  6324689.962954273,
  6347230.277480394,
  6364087.5,
  6364272.317130318,
  6384127.970561172,
  6384132.858884388,
  6384133.642231257,
  6391434.375,
  6391584.357495111,
  6391587.1726484755,
  6391875.356729987,
  6391925.323242055,
  6392496.988696254,
  6392859.992813377,
  6392877.7097100215,
  6398806.074027216,
  6411259.5695956135,
  6430553.318833017,
  6432774.412347517,
  6433266.607067543,
  6433924.773417076,
  6442751.922437774,
  6445746.875,
  6451976.251533856,
  6454548.4375,
  6454566.611227143,
  6455097.499999296,
  6456683.77042303,
  6464203.234935578,
  6464219.293255772,
  6465548.4375,
  6466196.555328,
  6466489.722449631,
  6466492.178129724,
  6467284.375,
  6472930.053967671,
  6475051.789275058,
  6475484.375,
  6475532.1163698565,
  6477843.628759042,
  6478815.625,
  6479753.125,
  6492304.6875,
  6492357.8125,
  6496264.0625,
  6498062.5,
  6508464.0625,
  6511036.838520094,
  6511122.878825806,
  6516896.834101093,
  6517370.3125,
  6517588.563651572,
  6523037.357605769,
  6529828.133206868,
  6535794.850179621,
  6536899.519315411,
  6540122.982247955,
  6540238.93716455,
  6540271.531899382,
  6541496.523669068,
  6552219.496800933,
  6552865.625,
  6552875.193282497,
  6556293.117438827,
  6577967.1875,
  6592498.127299166,
  6606261.908734295,
  6607526.446309192,
  6615879.514238279,
  6633326.17841807,
  6639489.156529045,
  6645651.44032431,
  6648528.125,
  6653554.6875,
  6655173.55230372,
  6655863.3431097735,
  6658726.443185236,
  6666032.000636315,
  6680556.022608152,
  6693474.988524131,
  6694570.3125,
  6699944.184130341,
  6703962.396695631,
  6704160.086962452,
  6704943.384933773,
  6721360.531430565,
  6731435.725928807,
  6732862.69270501,
  6742396.875,
  6750157.069466202,
  6751917.165264096,
  6752992.1875,
  6757478.125,
  6764751.5625,
  6765453.823780311,
  6767496.875,
  6795340.220823396,
  6799912.5,
  6816396.705970912,
  6821595.3125,
  6822638.597248429,
  6822917.611695682,
  6825792.1875,
  6844693.75,
  6852496.289600013,
  6855094.333095724,
  6855462.326195688,
  6855793.707262602,
  6855918.749325935,
  6856098.6254318245,
  6856279.586205223,
  6860275.0,
  6864036.312593633,
  6864861.337947928,
  6865712.5,
  6879417.591902216,
  6882189.66631164,
  6886814.092133719,
  6910634.375,
  6911928.027049641,
  6912199.444276428,
  6914246.355529393,
  6914421.875,
  6915046.663080895,
  6917353.679171382,
  6918871.875,
  6932687.224579017,
  6936685.234625235,
  6938915.753499886,
  6943875.271528397,
  6965887.5,
  6968153.317805092,
  6968426.5625,
  6968652.680989516,
  6969326.062819984,
  6969536.718524376,
  6969542.9068637015,
  6969830.83864261,
  6971314.380071097,
  6972119.961167927,
  6972120.722745279,
  6972320.3125,
  6976373.975332953,
  6982231.25,
  6986494.762410174,
  6987039.0625,
  6987539.979333383,
  7029499.667211107,
  7033105.880360608,
  7042340.625,
  7056637.597838184,
  7057428.125,
  7057438.777256733,
  7058022.902303999,
  7062387.314995076,
  7064029.6875,
  7064623.4375,
  7065604.621902001,
  7066596.855095408,
  7067418.75,
  7067548.651237835,
  7069804.6875,
  7070966.939162453,
  7076671.888951619,
  7081841.653132696,
  7083164.688490692,
  7085546.563091655,
  7090439.0625,
  7098923.4375,
  7100643.75,
  7100961.251880867,
  7101021.46857797,
  7101040.954066272,
  7101042.1875,
  7101070.304486073,
  7101967.1875,
  7119114.494449377,
  7119142.706738644,
  7119993.471642685,
  7120156.25,
  7120387.5,
  7120397.504038648,
  7120754.079197214,
  7122098.4375,
  7122749.094465627,
  7126332.8125,
  7126464.0625,
  7132846.875,
  7139929.6875,
  7143292.800249938,
  7143351.014692584,
  7145806.227342033,
  7145815.967921798,
  7146032.8125,
  7146334.7647872185,
  7146368.75,
  7146378.458172594,
  7146648.4375,
  7146990.626153024,
  7147481.904331879,
  7147978.171787495,
  7150636.668999937,
  7153875.0,
  7154065.625,
  7154959.598534582,
  7155929.6875,
  7159628.979901439,
  7214795.150746081,
  7215348.250305075,
  7223568.466567843,
  7223722.411912227,
  7231227.361970977,
  7236016.491392838,
  7239184.120389741,
  7271118.0936836535,
  7271159.375,
  7272755.62827271,
  7272766.169483935,
  7273498.859188631,
  7274816.197291385,
  7275573.065958184,
  7288536.183781996,
  7290547.95598082,
  7290587.24506547,
  7292195.15796612,
  7296486.459053229,
  7296558.822409572,
  7301244.276066935,
  7302043.75,
  7302115.896212602,
  7326172.3981050905,
  7326918.511018772,
  7327214.0625,
  7330618.75,
  7347440.095602483,
  7352156.121754101,
  7357736.494890946,
  7367009.34389197,
  7367735.221630025,
  7368079.446083401,
  7368216.153364367,
  7373823.4375,
  7374358.449346021,
  7375007.8125,
  7376081.066671932,
  7383918.75,
  7386906.20501683,
  7397754.6875,
  7401430.681537078,
  7427151.233944612,
  7428040.397663485,
  7428450.0,
  7428722.885781166,
  7429996.55943085,
  7431465.195108443,
  7447559.097592855,
  7447602.757488114,
  7470672.5777004585,
  7481824.845586427,
  7483934.347211462,
  7487076.5625,
  7488839.9415765675,
  7500006.25,
  7500836.395134269,
  7514646.24465726,
  7517896.189150857,
  7522315.458061876,
  7530896.607276471,
  7533835.239396883,
  7549100.281566813,
  7559627.97182353,
  7565146.83661196,
  7568092.264301567,
  7570412.432027111,
  7572172.26159283,
  7572197.711950465,
  7572737.893777727,
  7576065.207715668,
  7576670.153965108,
  7586321.875,
  7588195.251707606,
  7595362.505461409,
  7595401.498747828,
  7596086.066848376,
  7612506.443378305,
  7618515.39118069,
  7619398.139725847,
  7653207.162622215,
  7660807.8125,
  7676146.044969409,
  7683127.238787973,
  7695823.09776887,
  7703776.732833039,
  7706629.749099611,
  7709252.139933984,
  7725446.6373939235,
  7726557.8125,
  7734461.113694955,
  7734773.867217522,
  7735946.731765266,
  7738141.560988794,
  7740449.117850486,
  7745000.0,
  7745744.558603443,
  7745809.266897999,
  7746845.3125,
  7746900.66416836,
  7750683.946292246,
  7750848.4375,
  7751427.15533355,
  7753053.5451840265,
  7753597.785976197,
  7753674.897091295,
  7756391.152311827,
  7761139.0625,
  7770390.625,
  7770892.1875,
  7776121.875,
  7776493.75,
  7776756.988846817,
  7778085.9375,
  7778681.25,
  7784767.258994736,
  7784774.440559127,
  7786667.625403491,
  7786679.790152364,
  7789170.3125,
  7793947.635603636,
  7822791.632372629,
  7838968.242121083,
  7839329.6875,
  7843628.125,
  7844465.612131976,
  7845749.106219658,
  7849756.25,
  7852133.362998628,
  7852907.548539679,
  7857587.54657736,
  7880150.0,
  7883243.75,
  7892205.61200743,
  7898414.000005828,
  7898436.123579119,
  7898706.790486198,
  7899232.685738603,
  7904143.75,
  7905120.151713699,
  7916242.036869683,
  7918820.3160350695,
  7918825.0,
  7919669.331555442,
  7919935.799800518,
  7920408.145476777,
  7920415.625,
  7926922.095130529,
  7927137.682291421,
  7935427.341019714,
  7938127.131000082,
  7939946.569443002,
  7940570.25281298,
  7944178.125,
  7949877.715906706,
  7956116.897001985,
  7982109.375,
  7983366.395240506,
  7985389.1502094725,
  7996799.295779672,
  7997866.8535613185,
  8034257.8125,
  8034809.375,
  8034994.515342919,
  8036389.790021534,
  8038466.381966146,
  8039290.557657879,
  8041657.8125,
  8042293.066106248,
  8042412.3953181375,
  8043581.25,
  8044189.808884249,
  8047664.878539154,
  8055109.375,
  8063413.854806508,
  8066609.37238781,
  8075780.392069062,
  8077238.2149775075,
  8078544.450445712,
  8079368.333808584,
  8079539.648461029,
  8079745.3125,
  8080160.9375,
  8086657.087519399,
  8087416.855439277,
  8096735.036823418,
  8105136.900707128,
  8108277.5630520005,
  8113895.07974848,
  8113996.398662962,
  8118676.8131856285,
  8140610.9375,
  8153160.736983476,
  8155768.167216057,
  8156762.243277466,
  8156828.874551024,
  8159345.695481911,
  8159704.982094423,
  8159831.874398971,
  8160810.9375,
  8178712.5,
  8178891.6389876995,
  8179357.479943236,
  8183305.013713312,
  8186688.088251388,
  8205784.375,
  8245112.476645995,
  8247057.4976698505,
  8249771.352389682,
  8275849.775021633,
  8275957.662306401,
  8276030.980215641,
  8287891.443281613,
  8296356.25,
  8305550.15672711,
  8322825.654079661,
  8343010.827270233,
  8355466.45611102,
  8365987.469621887,
  8366748.4375,
  8383557.393473591,
  8387384.375,
  8403367.818729887,
  8404104.118490715,
  8408797.275258327,
  8426721.997055119,
  8426736.750988714,
  8434635.9375,
  8436923.13259699,
  8447150.862155644,
  8447308.546684612,
  8448869.037115455,
  8461332.586004118,
  8467905.518682485,
  8469360.4143886,
  8470163.745655095,
  8470257.2851717,
  8471307.824198764,
  8471408.97647027,
  8475051.274818866,
  8478850.0,
  8489110.685261818,
  8490554.5176276,
  8491959.375,
  8503795.40132972,
  8504709.375,
  8505718.782695398,
  8509527.431921665,
  8509702.743742788,
  8509704.696343774,
  8509709.139474593,
  8511585.9375,
  8513908.901903262,
  8528454.6875,
  8529515.64879859,
  8530264.0625,
  8532033.919870757,
  8542618.07049368,
  8542937.812946761,
  8552277.870632626,
  8558145.150324287,
  8560827.576800201,
  8563912.689026097,
  8595963.441127613,
  8600698.4375,
  8604509.124597415,
  8616246.41874446,
  8617727.177632144,
  8716803.191746434,
  8718575.515214734,
  8743674.898386894,
  8744463.578823779,
  8801959.05638304,
  8802316.852881536,
  8816921.244708845,
  8835612.984189479,
  8836453.125,
  8842844.514679093,
  8844840.625,
  8878528.341871582,
  8881756.25,
  8882220.3125,
  8929487.5,
  8934526.871778646,
  8950714.0625,
  8953850.353456711,
  8957453.609139135,
  8959626.22739169,
  8959795.3125,
  8960596.875,
  8965852.975186393,
  8968903.782620886,
  8970583.53594613,
  8973733.97914366,
  8978084.569004256,
  8979235.9375,
  8981777.353465248,
  8987239.9605755,
  8989512.815359963,
  8997095.3125,
  8997666.92170722,
  8997679.792797567,
  8999127.642454976,
  8999901.43034099,
  8999902.37703619,
  9000360.867239624,
  9000675.623208854,
  9012753.090119759,
  9016464.0625,
  9032564.0625,
  9033107.631302927,
  9067721.875,
  9070687.636026215,
  9072659.375,
  9074527.868409682,
  9085459.375,
  9092251.565174889,
  9092685.9375,
  9093028.682758693,
  9108881.218591005,
  9108893.7380893,
  9109812.5,
  9110204.6875,
  9110928.988959128,
  9111068.371907692,
  9111504.6875,
  9113209.375,
  9122298.252772884,
  9122722.512120308,
  9124526.5625,
  9125286.98357668,
  9125445.3125,
  9126162.5,
  9126802.206854604,
  9126842.1875,
  9129084.375,
  9131002.156021588,
  9131031.25,
  9131746.456449954,
  9134075.0,
  9134684.863656877,
  9138520.3125,
  9139216.227082277,
  9139795.085988956,
  9168178.125,
  9365160.318142092,
  9368031.25,
  9453001.5625],
 [15.614769966327133,
  37.1155661752404,
  13.892675712334887,
  20.241280469855834,
  21.063886871005202,
  11.561390086842907,
  62.59705518981089,
  8.466572728202882,
  14.724818400001537,
  36.736823126277464,
  12.696594241595042,
  13.747330507512231,
  63.16803435687988,
  65.12943017927766,
  55.52803137916615,
  39.473196921062566,
  23.650138656438582,
  28.85526830724879,
  16.71224874559791,
  20.740581028718143,
  97.29568450537086,
  59.028458918001974,
  29.59376884549594,
  9.255361172765275,
  54.49549892691492,
  66.94533111533197,
  9.846912369276689,
  10.838960265892576,
  137.02952433588655,
  100.59065489273542,
  15.7321768071244,
  10.11101478338409,
  51.32999251980287,
  5.045723598982221,
  5.306078524144217,
  24.92533260956835,
  75.10472290907104,
  12.890788075776735,
  49.0619809609241,
  6.5324412271671255,
  8.86106871565045,
  11.985834640184864,
  54.03718192212205,
  22.259530610820686,
  43.73564648322419,
  53.84838796365811,
  18.425042132310182,
  6.606676085773855,
  6.852174197851404,
  55.19922074208506,
  8.904173097776319,
  59.68584064725847,
  25.2829581277543,
  17.26891637906644,
  39.19146190815111,
  11.82911222496202,
  40.192983023911104,
  29.219065600560405,
  32.33581906739688,
  6.985867298033824,
  12.966902049470034,
  18.46815003973684,
  13.170893515533804,
  37.315019188323355,
  22.09894389414234,
  85.12608820679435,
  57.78325019931273,
  12.351955590685527,
  30.321818163672166,
  9.89057408230384,
  8.166375418081824,
  22.429128646269906,
  69.10336682472672,
  60.832179637845826,
  55.4057796795294,
  7.8057588135914315,
  22.150615242967262,
  9.576496073515967,
  71.06586551816848,
  104.2069923174784,
  14.525028134304282,
  15.955914680695997,
  44.30549979451328,
  46.906871472977436,
  14.631394093849668,
  19.133685562949957,
  12.770058374115244,
  18.74013905033442,
  6.260432171533588,
  12.861534608013159,
  18.582180018618086,
  79.32597696451268,
  49.788922400085085,
  13.258946002628806,
  36.34987955035869,
  15.453648470308137,
  56.01513576025013,
  56.57802468723617,
  12.104289185775832,
  191.72722170656115,
  45.148517697050735,
  25.0560839378684,
  16.98738677150621,
  14.115762999429755,
  10.376843506954224,
  15.455322342465514,
  21.899914704531447,
  29.24345675181684,
  26.831549118509326,
  70.43448575113271,
  30.20493553074924,
  18.868525731359064,
  20.384218229991372,
  42.37443682306643,
  21.451693058278334,
  55.49891763335433,
  33.83468312838282,
  8.07664690542092,
  28.983502769845,
  35.12579419774115,
  42.749059836789236,
  6.5874928950143605,
  45.14358337721294,
  14.67235505187391,
  30.103936395629724,
  106.2046557148694,
  51.11090250077219,
  32.42336091918229,
  60.252942815786966,
  24.920637316337665,
  95.31749870474144,
  9.577626043854025,
  19.42010652968349,
  36.637775658002134,
  43.378501388559236,
  92.87806174940461,
  38.23834302070315,
  19.0382642717784,
  59.756636032487116,
  72.49994725636412,
  25.92369526389607,
  11.275175673088256,
  24.11265099090807,
  19.399358338450764,
  6.134248414886759,
  71.62383245784048,
  11.816809855207078,
  9.241887621316186,
  21.569040032597886,
  13.374575119250867,
  6.888129305730314,
  30.191015873992004,
  54.32482614634452,
  17.41430837926818,
  6.1274970187611535,
  7.272878724738667,
  18.640481526646383,
  57.77403819906828,
  178.2087407296472,
  21.384764144749912,
  120.23896322770342,
  14.357779945648911,
  15.146816720241523,
  41.90352234949461,
  25.60882627327911,
  45.09960575412909,
  25.04507587472712,
  114.83155242167213,
  6.692180620390819,
  8.657987590190976,
  20.25656181925672,
  16.051525015351807,
  47.73468606226955,
  30.37558368317511,
  74.51203722740736,
  138.58255717857824,
  20.329349015973744,
  11.423477200525113,
  25.429523495426874,
  7.4688558278231545,
  8.922594759955786,
  8.903803042217083,
  130.81485424892517,
  18.124077186263115,
  74.62908220448168,
  52.60595312984717,
  54.20316714004254,
  78.15027000541264,
  38.106281991515786,
  48.60783821814187,
  11.891545942018722,
  17.19260620285983,
  83.93993021879407,
  81.24562335492317,
  187.6732339990674,
  41.56161234324317,
  43.09521258987494,
  50.701982387489124,
  23.800246406813173,
  13.804847946441853,
  72.18935478335115,
  11.006376061748195,
  19.565433417559536,
  26.73720479114686,
  21.68825414086999,
  5.974577115263925,
  31.790842810438928,
  42.22744027248871,
  17.105538022898187,
  49.0002138704361,
  59.818196732797844,
  83.35148476766825,
  23.217442306887218,
  24.082491556793475,
  30.41818903465886,
  104.91605934163945,
  53.08921356835019,
  26.18003708235912,
  32.404860511710915,
  91.72925653536242,
  26.29959689888458,
  16.962067089454333,
  33.48082550129065,
  118.15282494567504,
  167.96461690623929,
  46.07797365178146,
  22.812563466394543,
  219.5682676981942,
  23.184456195895123,
  23.611064226445656,
  5.102792608902199,
  73.22300232501627,
  64.73154076261136,
  54.55749497962934,
  5.354382169414058,
  17.37797619653832,
  32.92329519364646,
  17.537392266493853,
  11.645821800988376,
  15.832678224039043,
  84.52557778446868,
  12.92735129875925,
  10.157585505504422,
  74.70163628972608,
  72.22133960407947,
  57.6178778241228,
  11.029077861433978,
  25.503037441037886,
  5.819097764764022,
  16.23783808909282,
  21.83086440518341,
  123.17392115032798,
  29.616898679535318,
  9.314514975077797,
  5.328341815390537,
  87.63899685968752,
  31.701758846490645,
  42.184393593515225,
  46.748058475859374,
  52.81879919119535,
  93.7476097062621,
  21.16231721422618,
  16.979431698592435,
  29.267180313387556,
  39.28886108362691,
  33.427321198739264,
  139.11216565050213,
  13.169713912713586,
  5.433021377535052,
  151.63388231503578,
  10.229716872252059,
  7.285684956647565,
  103.19273936553219,
  22.533919494451325,
  34.79268730307008,
  11.200990938961816,
  36.73843471608245,
  44.58660611735263,
  39.732487474268446,
  8.120736012487809,
  121.77130222554614,
  33.01803531676436,
  6.281670738032824,
  7.843096545547782,
  22.38190456009417,
  110.92680870249507,
  15.738991499722053,
  20.820211172519667,
  20.661121200768935,
  57.5342239659596,
  137.42663219062536,
  16.517269056668876,
  74.04171029572589,
  6.702972943278148,
  211.18472754692783,
  8.003660565241498,
  12.268758038421767,
  25.647584810419737,
  9.426122626518538,
  17.514523329617575,
  57.18105031931751,
  16.059723036716527,
  40.46699357846643,
  28.34084889602217,
  32.77556313244104,
  104.65795209437849,
  22.429708683220895,
  47.80507906155397,
  17.686752789101448,
  9.230948792139776,
  93.70197663950312,
  85.4058752709449,
  77.00890290063982,
  52.89874182584843,
  37.87262801940629,
  144.54841721466283,
  51.08606870346067,
  21.209218939866595,
  29.298845925913593,
  102.99168130436092,
  84.52575992909023,
  105.77023454626985,
  25.709602800863735,
  121.3541794723686,
  29.088259994771544,
  8.047751647976714,
  58.503977207232765,
  36.78543761901659,
  33.70050751915862,
  5.909382543038365,
  18.8322551422548,
  69.08253167817767,
  35.46836429076977,
  136.90944721156797,
  53.60901716560561,
  103.03247809672973,
  139.73376895827113,
  25.306971868091786,
  21.06640343366051,
  9.02646639453157,
  98.94484438452403,
  26.32317251459509,
  33.971344917693855,
  64.10678228944015,
  125.10978370975555,
  107.48054325867645,
  26.597075968483608,
  24.32661660305525,
  44.67003024912576,
  11.1098165331589,
  89.02173642931777,
  16.902122266950656,
  75.82484670715624,
  83.06291495837716,
  5.035445093365805,
  48.40588333561513,
  10.812362896600913,
  9.499840748353096,
  31.789836924689762,
  105.8424488909115,
  14.857428517569838,
  87.8587700337127,
  42.756548569085425,
  61.44966277342754,
  12.429274468526797,
  78.04413831728195,
  5.592068706900726,
  101.73320624257796,
  84.15046236160833,
  50.23066155073736,
  18.923879426664016,
  11.862489629562008,
  37.06145111980516,
  31.499396268797277,
  6.975625174082379,
  8.62103577198507,
  14.317997139623957,
  14.781078694066746,
  34.346442066486944,
  22.689161776299247,
  7.0950765072860325,
  98.60640424052767,
  24.509887633028924,
  5.975054477249125,
  47.01351712067179,
  21.999729449569212,
  43.53668709539912,
  11.174344295336557,
  51.30840586204079,
  26.727349972655723,
  94.39102148331817,
  136.1257846525645,
  121.6408165983493,
  68.10479088366853,
  12.120544644856292,
  34.20277965063915,
  24.069865312011032,
  19.68285087389928,
  15.455005834525268,
  26.02364905207331,
  63.53270967768129,
  11.443067776574193,
  44.031286006183386,
  77.79269286293427,
  8.038905602635422,
  10.092418434217018,
  54.50888425602657,
  19.115829229354006,
  16.307896317582237,
  23.11646736866344,
  11.835805509931365,
  58.4208122945276,
  68.00264592152232,
  89.3055734145209,
  27.101168332401908,
  155.97445267401434,
  29.75750298426049,
  23.34075377987811,
  36.744082775375,
  123.6318940719091,
  25.333469751902417,
  63.12972683790582,
  25.649069995635703,
  32.480386006698666,
  40.943026369701954,
  68.73748110048398,
  46.93382425767469,
  27.290876316285612,
  11.096736525128136,
  83.04987298005474,
  15.823966867468592,
  38.523956579812904,
  72.66613647423546,
  11.14958433004702,
  21.156984439916876,
  6.34966576766438,
  45.58479520067375,
  32.057942678391235,
  46.7559274063442,
  98.28327272431125,
  16.29242147955353,
  22.480303727191078,
  34.55987189377694,
  29.552415073491265,
  38.01892592085635,
  98.41053956294138,
  11.877403044262572,
  71.79875517686352,
  7.505570473516513,
  40.51240594244258,
  29.30391403623791,
  30.845117844725422,
  102.85009012718496,
  50.86417940135735,
  21.972454116436467,
  76.13169016056835,
  97.45184125151427,
  99.49262479041113,
  54.422195965783885,
  12.242722896037312,
  73.16502773717922,
  33.41847808438594,
  28.507048199199936,
  67.10307868019618,
  22.078330501946507,
  204.34060600189628,
  11.691726991217637,
  88.35785954543536,
  20.749933532233122,
  31.9053429174757,
  33.467916987392506,
  10.037342858862516,
  72.05016479621312,
  50.654780784956976,
  78.5131746117998,
  91.70320522875704,
  19.64524028333861,
  25.693339753124917,
  39.013703458430754,
  5.33230240839648,
  18.707554922753562,
  73.84301921763249,
  26.4735725319478,
  24.348380639024608,
  139.57120572624285,
  45.85875800464462,
  6.967077432692455,
  6.146612087806292,
  11.857680192540462,
  20.851026189626353,
  157.93452591762576,
  83.07569543024667,
  7.000377524745289,
  10.75131256110481,
  52.49332646719358,
  25.536282365750004,
  20.984048434332085,
  5.307618469395774,
  61.66760754352518,
  5.35130183963476,
  74.05510952219996,
  5.521593362303416,
  85.61279378116981,
  145.19718946965094,
  21.20477264098052,
  6.738170534903634,
  5.645642976416612,
  34.98525051499954,
  253.0693348104893,
  10.77022032463575,
  63.290962328003914,
  78.6864083019725,
  119.33723308582427,
  51.381496507523124,
  33.633728270848664,
  33.540758597402984,
  86.15963548236772,
  29.244205983852503,
  8.892929680989385,
  23.39271274291585,
  24.725519906211517,
  38.54723292202948,
  13.186597545117456,
  23.55388387601943,
  22.22218962295973,
  24.043688715304707,
  75.7002835845927,
  77.46936925556162,
  34.95529504158023,
  27.986491478131114,
  41.38583741093936,
  168.69557015698499,
  51.330263082445136,
  8.544375716606,
  55.96632684865462,
  40.311210427630584,
  152.40173820332296,
  102.69393423195073,
  6.018391537242518,
  8.441306876419134,
  25.826756138396977,
  59.06497186292303,
  23.8222860786381,
  14.688606950860617,
  23.69184112206612,
  75.21439034755198,
  16.68983889948522,
  33.13411226442511,
  25.74321488759006,
  35.77688236233212,
  27.37701095088852,
  14.962366695162334,
  11.029499055834753,
  11.33609026871668,
  29.2449815769993,
  19.89243378046452,
  29.70986797661356,
  37.60482869983124,
  53.9068681593587,
  18.73091024211425,
  17.945932179552194,
  23.25616034807022,
  69.75187109237203,
  88.28838559047345,
  12.168609014198669,
  24.5818564188825,
  135.97955763187986,
  34.857019511838644,
  25.118165812364182,
  25.983493960498468,
  64.57200957966609,
  7.31256762020256,
  50.8477883103327,
  149.12229114760015,
  91.12360221342247,
  114.18278604092498,
  18.266763271835366,
  100.59778892481017,
  53.131207121831174,
  5.8591868981474295,
  9.713898889591693,
  8.835226283605628,
  32.445699287235826,
  33.70785473836754,
  72.6733288660623,
  68.90181621460378,
  91.90797461317592,
  105.97438914514262,
  106.49761238431915,
  18.775832338881457,
  103.16810133161476,
  51.0456249303767,
  35.854274526286865,
  19.78811839232225,
  16.922345301407262,
  15.347760266040202,
  71.89215367167961,
  32.16294301542492,
  8.013323662949785,
  7.429325407003637,
  43.09599128429415,
  39.06964310511996,
  28.272739595426216,
  26.588887944434944,
  65.52310324642417,
  5.485901274879486,
  29.961672290030542,
  43.956022194676756,
  29.09513891203491,
  57.47637549785293,
  87.90738239882754,
  35.88854007995246,
  15.701588175643574,
  12.09254340776239,
  68.22088344017432,
  30.405941217924497,
  10.88535589570661,
  25.866681264108188,
  28.354262025418628,
  73.55908828875974,
  17.743854888325586,
  14.441413162295266,
  34.89158735205282,
  79.18490323745337,
  11.490147812174165,
  21.972413593526667,
  12.61022157469751,
  66.5434867530931,
  35.65521726021844,
  35.573676638594314,
  56.166800836335106,
  21.346736277777737,
  36.80926426349184,
  121.19892488632408,
  12.591058563287843,
  5.712362902679958,
  51.96120504166045,
  29.851696317108235,
  81.24886287541852,
  77.46826200203334,
  76.9738457016018,
  22.482486632821992,
  44.010552091686144,
  27.259774117504442,
  63.60943484805226,
  70.20885095903978,
  111.88123674925092,
  53.68296710106666,
  6.610824803963359,
  79.72302145645614,
  108.52097776614347,
  89.99058027164054,
  8.285196438592337,
  119.72463814425151,
  18.609120844397495,
  5.752835120853492,
  99.54640388560067,
  17.770747664219364,
  63.15756993429879,
  127.42564224494123,
  10.156559748925629,
  134.37339674596353,
  18.251358853754205,
  20.55666365653747,
  59.50585546202748,
  20.336207060461767,
  10.456722595089174,
  15.991617823358645,
  42.14789229497926,
  17.405793380259134,
  98.60968083612752,
  108.9609508028696,
  5.028105310650802,
  30.14464100216142,
  106.03182502988724,
  5.747683174569619,
  76.75675658404526,
  87.13427072108821,
  13.195117341786863,
  7.487228825224888,
  144.28260138806618,
  106.84787936914017,
  30.539046495257235,
  78.0191970641961,
  8.126822988204411,
  46.40494212189582,
  24.75879695184923,
  40.042964297482555,
  5.037026685636164,
  110.48223772846353,
  50.59145327444017,
  57.77892632496538,
  249.2546921113635,
  40.09054385555049,
  8.781104630583034,
  52.3050205672632,
  52.30861883394773,
  5.444064360118548,
  40.88697821082237,
  48.99971081339312,
  8.892545279309934,
  16.723446677254188,
  7.64798764561481,
  36.74418018413055,
  53.06731000634828,
  32.28435134531958,
  19.703073959517653,
  42.21498596844757,
  9.586577701233379,
  24.737370069913993,
  39.17806036478852,
  73.37350670264945,
  68.00743471370616,
  44.8090884229112,
  39.91812550577923,
  6.426381676063272,
  13.947752123306596,
  41.817917884904254,
  24.38062669217548,
  128.06747247482542,
  72.83960768606997,
  56.55022160234503,
  73.72885583509289,
  8.30394320840856,
  19.761441519044602,
  6.689705923744102,
  14.975325230820555,
  10.778495354050087,
  128.2785998587404,
  21.765739383088576,
  107.96631902737053,
  53.04970611270958,
  134.25320360569648,
  38.873360476262484,
  30.470622215831224,
  19.577423782356846,
  75.17296550079249,
  25.686821762365195,
  28.795067882788995,
  5.828100233270641,
  81.1293065986106,
  40.0129442549944,
  24.405651247752306,
  5.900633080357883,
  70.02639580439629,
  15.202168508948457,
  22.815517438430888,
  8.233563763525833,
  80.82560617130285,
  40.01493574352517,
  9.063632043623945,
  19.78912551333547,
  17.175829002480473,
  82.66097900120657,
  18.96069815747877,
  32.290137914057965,
  11.484500112572016,
  31.90179012572043,
  100.50298998192858,
  117.73661364747552,
  123.43148312259714,
  58.21147976787574,
  8.482774726718528,
  5.738278628476706,
  10.641417115782785,
  79.45292133932595,
  63.395531297118666,
  25.87951917415535,
  15.739635204042116,
  12.02412149393823,
  5.905657727872323,
  49.93939095351203,
  73.43586526424619,
  7.346147092678332,
  19.996004017435702,
  16.297166521478413,
  68.10973712466347,
  20.7791888913556,
  74.26908225524363,
  21.669008205573935,
  90.67294036640165,
  8.770035245608964,
  27.207922160090327,
  59.06418362688859,
  5.4602714288084755,
  80.74891449325989,
  53.734605509756385,
  25.396963004938815,
  87.4023107488662,
  13.133148347753895,
  57.989542583577425,
  24.321543029841674,
  34.080013057012046,
  23.032670177587498,
  60.782494176893785,
  13.918222930219114,
  11.078964878958166,
  43.727566256015535,
  59.00982334483704,
  64.19268576754467,
  42.99093540147429,
  71.72355378448395,
  44.66086994206202,
  47.29178287896077,
  10.85074189734108,
  52.34692814906521,
  30.24969039753773,
  21.330818493425465,
  76.1754792877278,
  50.3165533086051,
  12.310247719617765,
  91.16205018518662,
  29.652968860364304,
  17.817151456560303,
  55.536675217166824,
  112.96779596297503,
  81.20440966609611,
  5.738579593055229,
  55.20565650611302,
  6.889970630558052,
  10.791198613732218,
  29.837414806711134,
  20.379709612445208,
  41.44132575692947,
  41.72337927013647])

Energy spectrum

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

([3492620.530665972,
  3582442.0812735655,
  3585565.319624986,
  3593032.2745178808,
  3611177.580869366,
  3614108.073302657,
  3616160.892630776,
  3693601.88804027,
  3732662.441962715,
  3803672.662098161,
  3856413.315365188,
  3924180.8969368823,
  3946479.039064631,
  3947820.0847265646,
  3948395.2593749003,
  3982932.808839554,
  3985380.2526885793,
  4004219.513812185,
  4004237.410135303,
  4004249.4211029704,
  4004315.625,
  4004326.5519919926,
  4032283.745416332,
  4048985.846178447,
  4063593.6444698013,
  4065810.262041953,
  4070489.615134211,
  4071332.619565964,
  4103654.834075178,
  4108865.1873446065,
  4111239.255548117,
  4115083.610033039,
  4140226.5625,
  4142648.6840628483,
  4157467.665319404,
  4158065.6581533193,
  4176712.3629468787,
  4176718.0829887926,
  4183423.4375,
  4183431.5333238444,
  4207255.635936901,
  4230417.674915464,
  4240376.5625,
  4245120.833868845,
  4246034.375,
  4278065.625,
  4278338.187965376,
  4296439.785176602,
  4314606.67847943,
  4337496.875,
  4338624.447477073,
  4344571.082869163,
  4368767.926705537,
  4369570.858679403,
  4370025.0,
  4370050.406730455,
  4370078.468745263,
  4382167.1875,
  4385253.618855921,
  4394957.270319626,
  4396710.300909533,
  4400208.5068542315,
  4401721.317354548,
  4403828.125,
  4404592.014905767,
  4405027.614335727,
  4405306.25,
  4411033.476109745,
  4432089.0625,
  4444100.594728054,
  4444638.828062689,
  4444921.463977997,
  4444945.3125,
  4445785.816980325,
  4453308.859898969,
  4467647.402255139,
  4467762.413638079,
  4507940.901302447,
  4513187.907767356,
  4513339.0625,
  4515787.049710502,
  4528356.492906535,
  4541787.5,
  4544634.458786175,
  4545223.42308966,
  4545577.350766959,
  4547352.408610105,
  4547697.608782832,
  4558034.321257396,
  4558949.32115245,
  4559228.836067757,
  4566832.8125,
  4566995.3125,
  4567151.0378156435,
  4567252.79983913,
  4567267.586596807,
  4567650.0,
  4567653.125,
  4568517.733738784,
  4568700.115373324,
  4568706.25,
  4588297.150567199,
  4588322.324956934,
  4606041.277647058,
  4607407.272407189,
  4609157.94107543,
  4622816.727336632,
  4624307.8125,
  4631048.34565392,
  4638880.898431857,
  4661785.110310686,
  4662986.689381507,
  4692118.131158823,
  4700800.0,
  4708326.0583390845,
  4715356.25,
  4718447.077253479,
  4718652.705737127,
  4723722.554191789,
  4738390.625,
  4746995.3125,
  4753088.111919016,
  4761115.5090038525,
  4761337.427545323,
  4762190.138718305,
  4771840.625,
  4778800.644308123,
  4780726.5625,
  4792370.3125,
  4792425.67511725,
  4793339.0625,
  4793453.402516461,
  4793502.388632653,
  4804800.0,
  4810887.5,
  4817915.314200792,
  4823158.910894705,
  4843160.150324343,
  4850747.3125097,
  4852696.875,
  4870146.727989322,
  4874169.022973611,
  4874746.601720019,
  4885514.8859380605,
  4889250.432388543,
  4900423.4375,
  4913836.337753746,
  4920548.32970231,
  4924215.132579058,
  4936441.705124193,
  4936867.781244053,
  4971395.955267249,
  4997933.071737891,
  5004409.9010975715,
  5006145.683023668,
  5014930.233848478,
  5015036.574655041,
  5016085.9375,
  5041387.378949439,
  5043382.149309411,
  5047707.686714367,
  5059738.900431105,
  5063029.740341798,
  5072741.591584372,
  5074093.582756976,
  5074668.611125128,
  5074720.107495489,
  5083679.904929095,
  5087359.006914243,
  5088588.141981034,
  5090221.2414149875,
  5090597.138631836,
  5090893.75,
  5093431.25,
  5098627.692084908,
  5104219.234343314,
  5105066.562350477,
  5110621.086165099,
  5115426.383465704,
  5144596.4875148125,
  5147997.910924866,
  5149588.509453688,
  5151425.0,
  5168914.855442729,
  5173954.6875,
  5202698.4375,
  5202700.0,
  5202898.840043367,
  5202918.736277276,
  5205207.531347487,
  5215830.12677419,
  5246987.987188586,
  5305675.148945387,
  5340362.5,
  5341082.8125,
  5342370.3125,
  5344550.0,
  5356625.0,
  5360156.113141704,
  5386041.450878613,
  5393534.375,
  5393540.545964344,
  5398086.782300612,
  5419719.49582887,
  5429926.08184323,
  5440662.895200803,
  5455821.875,
  5479009.375,
  5479677.034038522,
  5523718.320207222,
  5561756.25,
  5568504.477091591,
  5575988.719570688,
  5601538.788035155,
  5601559.112705029,
  5614991.1128464,
  5631440.625,
  5633206.119030191,
  5633585.9375,
  5633979.6875,
  5677488.938883729,
  5678219.195279469,
  5682221.5841641305,
  5690881.25,
  5707609.81274805,
  5731804.6875,
  5742124.622823982,
  5774245.3125,
  5777536.395349934,
  5782454.376485919,
  5792843.969078138,
  5797463.898883178,
  5822959.356834526,
  5826571.875,
  5833190.984185385,
  5833355.207311864,
  5835251.343775273,
  5835509.922134748,
  5835565.932978031,
  5875859.590684232,
  5876278.717555765,
  5892808.722669422,
  5903393.025912761,
  5911587.5,
  5956122.610491333,
  5964578.125,
  5981715.560970378,
  5990124.825824275,
  6041239.131424672,
  6066414.365561938,
  6137851.095705542,
  6139587.3817310715,
  6157208.400800219,
  6165635.91745936,
  6190257.4107712675,
  6194263.533878264,
  6195880.800181757,
  6211297.072798556,
  6211332.096785197,
  6211990.438800751,
  6214761.580308128,
  6234113.528688682,
  6238683.020853411,
  6243307.8125,
  6253770.3125,
  6284650.0,
  6289653.125,
  6291201.004408521,
  6291447.750224567,
  6299438.434041375,
  6324689.962954273,
  6347230.277480394,
  6364087.5,
  6364272.317130318,
  6384127.970561172,
  6384132.858884388,
  6384133.642231257,
  6391434.375,
  6391584.357495111,
  6391587.1726484755,
  6391875.356729987,
  6391925.323242055,
  6392496.988696254,
  6392859.992813377,
  6392877.7097100215,
  6398806.074027216,
  6411259.5695956135,
  6430553.318833017,
  6432774.412347517,
  6433266.607067543,
  6433924.773417076,
  6442751.922437774,
  6445746.875,
  6451976.251533856,
  6454548.4375,
  6454566.611227143,
  6455097.499999296,
  6456683.77042303,
  6464203.234935578,
  6464219.293255772,
  6465548.4375,
  6466196.555328,
  6466489.722449631,
  6466492.178129724,
  6467284.375,
  6472930.053967671,
  6475051.789275058,
  6475484.375,
  6475532.1163698565,
  6477843.628759042,
  6478815.625,
  6479753.125,
  6492304.6875,
  6492357.8125,
  6496264.0625,
  6498062.5,
  6508464.0625,
  6511036.838520094,
  6511122.878825806,
  6516896.834101093,
  6517370.3125,
  6517588.563651572,
  6523037.357605769,
  6529828.133206868,
  6535794.850179621,
  6536899.519315411,
  6540122.982247955,
  6540238.93716455,
  6540271.531899382,
  6541496.523669068,
  6552219.496800933,
  6552865.625,
  6552875.193282497,
  6556293.117438827,
  6577967.1875,
  6592498.127299166,
  6606261.908734295,
  6607526.446309192,
  6615879.514238279,
  6633326.17841807,
  6639489.156529045,
  6645651.44032431,
  6648528.125,
  6653554.6875,
  6655173.55230372,
  6655863.3431097735,
  6658726.443185236,
  6666032.000636315,
  6680556.022608152,
  6693474.988524131,
  6694570.3125,
  6699944.184130341,
  6703962.396695631,
  6704160.086962452,
  6704943.384933773,
  6721360.531430565,
  6731435.725928807,
  6732862.69270501,
  6742396.875,
  6750157.069466202,
  6751917.165264096,
  6752992.1875,
  6757478.125,
  6764751.5625,
  6765453.823780311,
  6767496.875,
  6795340.220823396,
  6799912.5,
  6816396.705970912,
  6821595.3125,
  6822638.597248429,
  6822917.611695682,
  6825792.1875,
  6844693.75,
  6852496.289600013,
  6855094.333095724,
  6855462.326195688,
  6855793.707262602,
  6855918.749325935,
  6856098.6254318245,
  6856279.586205223,
  6860275.0,
  6864036.312593633,
  6864861.337947928,
  6865712.5,
  6879417.591902216,
  6882189.66631164,
  6886814.092133719,
  6910634.375,
  6911928.027049641,
  6912199.444276428,
  6914246.355529393,
  6914421.875,
  6915046.663080895,
  6917353.679171382,
  6918871.875,
  6932687.224579017,
  6936685.234625235,
  6938915.753499886,
  6943875.271528397,
  6965887.5,
  6968153.317805092,
  6968426.5625,
  6968652.680989516,
  6969326.062819984,
  6969536.718524376,
  6969542.9068637015,
  6969830.83864261,
  6971314.380071097,
  6972119.961167927,
  6972120.722745279,
  6972320.3125,
  6976373.975332953,
  6982231.25,
  6986494.762410174,
  6987039.0625,
  6987539.979333383,
  7029499.667211107,
  7033105.880360608,
  7042340.625,
  7056637.597838184,
  7057428.125,
  7057438.777256733,
  7058022.902303999,
  7062387.314995076,
  7064029.6875,
  7064623.4375,
  7065604.621902001,
  7066596.855095408,
  7067418.75,
  7067548.651237835,
  7069804.6875,
  7070966.939162453,
  7076671.888951619,
  7081841.653132696,
  7083164.688490692,
  7085546.563091655,
  7090439.0625,
  7098923.4375,
  7100643.75,
  7100961.251880867,
  7101021.46857797,
  7101040.954066272,
  7101042.1875,
  7101070.304486073,
  7101967.1875,
  7119114.494449377,
  7119142.706738644,
  7119993.471642685,
  7120156.25,
  7120387.5,
  7120397.504038648,
  7120754.079197214,
  7122098.4375,
  7122749.094465627,
  7126332.8125,
  7126464.0625,
  7132846.875,
  7139929.6875,
  7143292.800249938,
  7143351.014692584,
  7145806.227342033,
  7145815.967921798,
  7146032.8125,
  7146334.7647872185,
  7146368.75,
  7146378.458172594,
  7146648.4375,
  7146990.626153024,
  7147481.904331879,
  7147978.171787495,
  7150636.668999937,
  7153875.0,
  7154065.625,
  7154959.598534582,
  7155929.6875,
  7159628.979901439,
  7214795.150746081,
  7215348.250305075,
  7223568.466567843,
  7223722.411912227,
  7231227.361970977,
  7236016.491392838,
  7239184.120389741,
  7271118.0936836535,
  7271159.375,
  7272755.62827271,
  7272766.169483935,
  7273498.859188631,
  7274816.197291385,
  7275573.065958184,
  7288536.183781996,
  7290547.95598082,
  7290587.24506547,
  7292195.15796612,
  7296486.459053229,
  7296558.822409572,
  7301244.276066935,
  7302043.75,
  7302115.896212602,
  7326172.3981050905,
  7326918.511018772,
  7327214.0625,
  7330618.75,
  7347440.095602483,
  7352156.121754101,
  7357736.494890946,
  7367009.34389197,
  7367735.221630025,
  7368079.446083401,
  7368216.153364367,
  7373823.4375,
  7374358.449346021,
  7375007.8125,
  7376081.066671932,
  7383918.75,
  7386906.20501683,
  7397754.6875,
  7401430.681537078,
  7427151.233944612,
  7428040.397663485,
  7428450.0,
  7428722.885781166,
  7429996.55943085,
  7431465.195108443,
  7447559.097592855,
  7447602.757488114,
  7470672.5777004585,
  7481824.845586427,
  7483934.347211462,
  7487076.5625,
  7488839.9415765675,
  7500006.25,
  7500836.395134269,
  7514646.24465726,
  7517896.189150857,
  7522315.458061876,
  7530896.607276471,
  7533835.239396883,
  7549100.281566813,
  7559627.97182353,
  7565146.83661196,
  7568092.264301567,
  7570412.432027111,
  7572172.26159283,
  7572197.711950465,
  7572737.893777727,
  7576065.207715668,
  7576670.153965108,
  7586321.875,
  7588195.251707606,
  7595362.505461409,
  7595401.498747828,
  7596086.066848376,
  7612506.443378305,
  7618515.39118069,
  7619398.139725847,
  7653207.162622215,
  7660807.8125,
  7676146.044969409,
  7683127.238787973,
  7695823.09776887,
  7703776.732833039,
  7706629.749099611,
  7709252.139933984,
  7725446.6373939235,
  7726557.8125,
  7734461.113694955,
  7734773.867217522,
  7735946.731765266,
  7738141.560988794,
  7740449.117850486,
  7745000.0,
  7745744.558603443,
  7745809.266897999,
  7746845.3125,
  7746900.66416836,
  7750683.946292246,
  7750848.4375,
  7751427.15533355,
  7753053.5451840265,
  7753597.785976197,
  7753674.897091295,
  7756391.152311827,
  7761139.0625,
  7770390.625,
  7770892.1875,
  7776121.875,
  7776493.75,
  7776756.988846817,
  7778085.9375,
  7778681.25,
  7784767.258994736,
  7784774.440559127,
  7786667.625403491,
  7786679.790152364,
  7789170.3125,
  7793947.635603636,
  7822791.632372629,
  7838968.242121083,
  7839329.6875,
  7843628.125,
  7844465.612131976,
  7845749.106219658,
  7849756.25,
  7852133.362998628,
  7852907.548539679,
  7857587.54657736,
  7880150.0,
  7883243.75,
  7892205.61200743,
  7898414.000005828,
  7898436.123579119,
  7898706.790486198,
  7899232.685738603,
  7904143.75,
  7905120.151713699,
  7916242.036869683,
  7918820.3160350695,
  7918825.0,
  7919669.331555442,
  7919935.799800518,
  7920408.145476777,
  7920415.625,
  7926922.095130529,
  7927137.682291421,
  7935427.341019714,
  7938127.131000082,
  7939946.569443002,
  7940570.25281298,
  7944178.125,
  7949877.715906706,
  7956116.897001985,
  7982109.375,
  7983366.395240506,
  7985389.1502094725,
  7996799.295779672,
  7997866.8535613185,
  8034257.8125,
  8034809.375,
  8034994.515342919,
  8036389.790021534,
  8038466.381966146,
  8039290.557657879,
  8041657.8125,
  8042293.066106248,
  8042412.3953181375,
  8043581.25,
  8044189.808884249,
  8047664.878539154,
  8055109.375,
  8063413.854806508,
  8066609.37238781,
  8075780.392069062,
  8077238.2149775075,
  8078544.450445712,
  8079368.333808584,
  8079539.648461029,
  8079745.3125,
  8080160.9375,
  8086657.087519399,
  8087416.855439277,
  8096735.036823418,
  8105136.900707128,
  8108277.5630520005,
  8113895.07974848,
  8113996.398662962,
  8118676.8131856285,
  8140610.9375,
  8153160.736983476,
  8155768.167216057,
  8156762.243277466,
  8156828.874551024,
  8159345.695481911,
  8159704.982094423,
  8159831.874398971,
  8160810.9375,
  8178712.5,
  8178891.6389876995,
  8179357.479943236,
  8183305.013713312,
  8186688.088251388,
  8205784.375,
  8245112.476645995,
  8247057.4976698505,
  8249771.352389682,
  8275849.775021633,
  8275957.662306401,
  8276030.980215641,
  8287891.443281613,
  8296356.25,
  8305550.15672711,
  8322825.654079661,
  8343010.827270233,
  8355466.45611102,
  8365987.469621887,
  8366748.4375,
  8383557.393473591,
  8387384.375,
  8403367.818729887,
  8404104.118490715,
  8408797.275258327,
  8426721.997055119,
  8426736.750988714,
  8434635.9375,
  8436923.13259699,
  8447150.862155644,
  8447308.546684612,
  8448869.037115455,
  8461332.586004118,
  8467905.518682485,
  8469360.4143886,
  8470163.745655095,
  8470257.2851717,
  8471307.824198764,
  8471408.97647027,
  8475051.274818866,
  8478850.0,
  8489110.685261818,
  8490554.5176276,
  8491959.375,
  8503795.40132972,
  8504709.375,
  8505718.782695398,
  8509527.431921665,
  8509702.743742788,
  8509704.696343774,
  8509709.139474593,
  8511585.9375,
  8513908.901903262,
  8528454.6875,
  8529515.64879859,
  8530264.0625,
  8532033.919870757,
  8542618.07049368,
  8542937.812946761,
  8552277.870632626,
  8558145.150324287,
  8560827.576800201,
  8563912.689026097,
  8595963.441127613,
  8600698.4375,
  8604509.124597415,
  8616246.41874446,
  8617727.177632144,
  8716803.191746434,
  8718575.515214734,
  8743674.898386894,
  8744463.578823779,
  8801959.05638304,
  8802316.852881536,
  8816921.244708845,
  8835612.984189479,
  8836453.125,
  8842844.514679093,
  8844840.625,
  8878528.341871582,
  8881756.25,
  8882220.3125,
  8929487.5,
  8934526.871778646,
  8950714.0625,
  8953850.353456711,
  8957453.609139135,
  8959626.22739169,
  8959795.3125,
  8960596.875,
  8965852.975186393,
  8968903.782620886,
  8970583.53594613,
  8973733.97914366,
  8978084.569004256,
  8979235.9375,
  8981777.353465248,
  8987239.9605755,
  8989512.815359963,
  8997095.3125,
  8997666.92170722,
  8997679.792797567,
  8999127.642454976,
  8999901.43034099,
  8999902.37703619,
  9000360.867239624,
  9000675.623208854,
  9012753.090119759,
  9016464.0625,
  9032564.0625,
  9033107.631302927,
  9067721.875,
  9070687.636026215,
  9072659.375,
  9074527.868409682,
  9085459.375,
  9092251.565174889,
  9092685.9375,
  9093028.682758693,
  9108881.218591005,
  9108893.7380893,
  9109812.5,
  9110204.6875,
  9110928.988959128,
  9111068.371907692,
  9111504.6875,
  9113209.375,
  9122298.252772884,
  9122722.512120308,
  9124526.5625,
  9125286.98357668,
  9125445.3125,
  9126162.5,
  9126802.206854604,
  9126842.1875,
  9129084.375,
  9131002.156021588,
  9131031.25,
  9131746.456449954,
  9134075.0,
  9134684.863656877,
  9138520.3125,
  9139216.227082277,
  9139795.085988956,
  9168178.125,
  9365160.318142092,
  9368031.25,
  9453001.5625],
 [15.614769966327133,
  37.1155661752404,
  13.892675712334887,
  20.241280469855834,
  21.063886871005202,
  11.561390086842907,
  62.59705518981089,
  8.466572728202882,
  14.724818400001537,
  36.736823126277464,
  12.696594241595042,
  13.747330507512231,
  63.16803435687988,
  65.12943017927766,
  55.52803137916615,
  39.473196921062566,
  23.650138656438582,
  28.85526830724879,
  16.71224874559791,
  20.740581028718143,
  97.29568450537086,
  59.028458918001974,
  29.59376884549594,
  9.255361172765275,
  54.49549892691492,
  66.94533111533197,
  9.846912369276689,
  10.838960265892576,
  137.02952433588655,
  100.59065489273542,
  15.7321768071244,
  10.11101478338409,
  51.32999251980287,
  5.045723598982221,
  5.306078524144217,
  24.92533260956835,
  75.10472290907104,
  12.890788075776735,
  49.0619809609241,
  6.5324412271671255,
  8.86106871565045,
  11.985834640184864,
  54.03718192212205,
  22.259530610820686,
  43.73564648322419,
  53.84838796365811,
  18.425042132310182,
  6.606676085773855,
  6.852174197851404,
  55.19922074208506,
  8.904173097776319,
  59.68584064725847,
  25.2829581277543,
  17.26891637906644,
  39.19146190815111,
  11.82911222496202,
  40.192983023911104,
  29.219065600560405,
  32.33581906739688,
  6.985867298033824,
  12.966902049470034,
  18.46815003973684,
  13.170893515533804,
  37.315019188323355,
  22.09894389414234,
  85.12608820679435,
  57.78325019931273,
  12.351955590685527,
  30.321818163672166,
  9.89057408230384,
  8.166375418081824,
  22.429128646269906,
  69.10336682472672,
  60.832179637845826,
  55.4057796795294,
  7.8057588135914315,
  22.150615242967262,
  9.576496073515967,
  71.06586551816848,
  104.2069923174784,
  14.525028134304282,
  15.955914680695997,
  44.30549979451328,
  46.906871472977436,
  14.631394093849668,
  19.133685562949957,
  12.770058374115244,
  18.74013905033442,
  6.260432171533588,
  12.861534608013159,
  18.582180018618086,
  79.32597696451268,
  49.788922400085085,
  13.258946002628806,
  36.34987955035869,
  15.453648470308137,
  56.01513576025013,
  56.57802468723617,
  12.104289185775832,
  191.72722170656115,
  45.148517697050735,
  25.0560839378684,
  16.98738677150621,
  14.115762999429755,
  10.376843506954224,
  15.455322342465514,
  21.899914704531447,
  29.24345675181684,
  26.831549118509326,
  70.43448575113271,
  30.20493553074924,
  18.868525731359064,
  20.384218229991372,
  42.37443682306643,
  21.451693058278334,
  55.49891763335433,
  33.83468312838282,
  8.07664690542092,
  28.983502769845,
  35.12579419774115,
  42.749059836789236,
  6.5874928950143605,
  45.14358337721294,
  14.67235505187391,
  30.103936395629724,
  106.2046557148694,
  51.11090250077219,
  32.42336091918229,
  60.252942815786966,
  24.920637316337665,
  95.31749870474144,
  9.577626043854025,
  19.42010652968349,
  36.637775658002134,
  43.378501388559236,
  92.87806174940461,
  38.23834302070315,
  19.0382642717784,
  59.756636032487116,
  72.49994725636412,
  25.92369526389607,
  11.275175673088256,
  24.11265099090807,
  19.399358338450764,
  6.134248414886759,
  71.62383245784048,
  11.816809855207078,
  9.241887621316186,
  21.569040032597886,
  13.374575119250867,
  6.888129305730314,
  30.191015873992004,
  54.32482614634452,
  17.41430837926818,
  6.1274970187611535,
  7.272878724738667,
  18.640481526646383,
  57.77403819906828,
  178.2087407296472,
  21.384764144749912,
  120.23896322770342,
  14.357779945648911,
  15.146816720241523,
  41.90352234949461,
  25.60882627327911,
  45.09960575412909,
  25.04507587472712,
  114.83155242167213,
  6.692180620390819,
  8.657987590190976,
  20.25656181925672,
  16.051525015351807,
  47.73468606226955,
  30.37558368317511,
  74.51203722740736,
  138.58255717857824,
  20.329349015973744,
  11.423477200525113,
  25.429523495426874,
  7.4688558278231545,
  8.922594759955786,
  8.903803042217083,
  130.81485424892517,
  18.124077186263115,
  74.62908220448168,
  52.60595312984717,
  54.20316714004254,
  78.15027000541264,
  38.106281991515786,
  48.60783821814187,
  11.891545942018722,
  17.19260620285983,
  83.93993021879407,
  81.24562335492317,
  187.6732339990674,
  41.56161234324317,
  43.09521258987494,
  50.701982387489124,
  23.800246406813173,
  13.804847946441853,
  72.18935478335115,
  11.006376061748195,
  19.565433417559536,
  26.73720479114686,
  21.68825414086999,
  5.974577115263925,
  31.790842810438928,
  42.22744027248871,
  17.105538022898187,
  49.0002138704361,
  59.818196732797844,
  83.35148476766825,
  23.217442306887218,
  24.082491556793475,
  30.41818903465886,
  104.91605934163945,
  53.08921356835019,
  26.18003708235912,
  32.404860511710915,
  91.72925653536242,
  26.29959689888458,
  16.962067089454333,
  33.48082550129065,
  118.15282494567504,
  167.96461690623929,
  46.07797365178146,
  22.812563466394543,
  219.5682676981942,
  23.184456195895123,
  23.611064226445656,
  5.102792608902199,
  73.22300232501627,
  64.73154076261136,
  54.55749497962934,
  5.354382169414058,
  17.37797619653832,
  32.92329519364646,
  17.537392266493853,
  11.645821800988376,
  15.832678224039043,
  84.52557778446868,
  12.92735129875925,
  10.157585505504422,
  74.70163628972608,
  72.22133960407947,
  57.6178778241228,
  11.029077861433978,
  25.503037441037886,
  5.819097764764022,
  16.23783808909282,
  21.83086440518341,
  123.17392115032798,
  29.616898679535318,
  9.314514975077797,
  5.328341815390537,
  87.63899685968752,
  31.701758846490645,
  42.184393593515225,
  46.748058475859374,
  52.81879919119535,
  93.7476097062621,
  21.16231721422618,
  16.979431698592435,
  29.267180313387556,
  39.28886108362691,
  33.427321198739264,
  139.11216565050213,
  13.169713912713586,
  5.433021377535052,
  151.63388231503578,
  10.229716872252059,
  7.285684956647565,
  103.19273936553219,
  22.533919494451325,
  34.79268730307008,
  11.200990938961816,
  36.73843471608245,
  44.58660611735263,
  39.732487474268446,
  8.120736012487809,
  121.77130222554614,
  33.01803531676436,
  6.281670738032824,
  7.843096545547782,
  22.38190456009417,
  110.92680870249507,
  15.738991499722053,
  20.820211172519667,
  20.661121200768935,
  57.5342239659596,
  137.42663219062536,
  16.517269056668876,
  74.04171029572589,
  6.702972943278148,
  211.18472754692783,
  8.003660565241498,
  12.268758038421767,
  25.647584810419737,
  9.426122626518538,
  17.514523329617575,
  57.18105031931751,
  16.059723036716527,
  40.46699357846643,
  28.34084889602217,
  32.77556313244104,
  104.65795209437849,
  22.429708683220895,
  47.80507906155397,
  17.686752789101448,
  9.230948792139776,
  93.70197663950312,
  85.4058752709449,
  77.00890290063982,
  52.89874182584843,
  37.87262801940629,
  144.54841721466283,
  51.08606870346067,
  21.209218939866595,
  29.298845925913593,
  102.99168130436092,
  84.52575992909023,
  105.77023454626985,
  25.709602800863735,
  121.3541794723686,
  29.088259994771544,
  8.047751647976714,
  58.503977207232765,
  36.78543761901659,
  33.70050751915862,
  5.909382543038365,
  18.8322551422548,
  69.08253167817767,
  35.46836429076977,
  136.90944721156797,
  53.60901716560561,
  103.03247809672973,
  139.73376895827113,
  25.306971868091786,
  21.06640343366051,
  9.02646639453157,
  98.94484438452403,
  26.32317251459509,
  33.971344917693855,
  64.10678228944015,
  125.10978370975555,
  107.48054325867645,
  26.597075968483608,
  24.32661660305525,
  44.67003024912576,
  11.1098165331589,
  89.02173642931777,
  16.902122266950656,
  75.82484670715624,
  83.06291495837716,
  5.035445093365805,
  48.40588333561513,
  10.812362896600913,
  9.499840748353096,
  31.789836924689762,
  105.8424488909115,
  14.857428517569838,
  87.8587700337127,
  42.756548569085425,
  61.44966277342754,
  12.429274468526797,
  78.04413831728195,
  5.592068706900726,
  101.73320624257796,
  84.15046236160833,
  50.23066155073736,
  18.923879426664016,
  11.862489629562008,
  37.06145111980516,
  31.499396268797277,
  6.975625174082379,
  8.62103577198507,
  14.317997139623957,
  14.781078694066746,
  34.346442066486944,
  22.689161776299247,
  7.0950765072860325,
  98.60640424052767,
  24.509887633028924,
  5.975054477249125,
  47.01351712067179,
  21.999729449569212,
  43.53668709539912,
  11.174344295336557,
  51.30840586204079,
  26.727349972655723,
  94.39102148331817,
  136.1257846525645,
  121.6408165983493,
  68.10479088366853,
  12.120544644856292,
  34.20277965063915,
  24.069865312011032,
  19.68285087389928,
  15.455005834525268,
  26.02364905207331,
  63.53270967768129,
  11.443067776574193,
  44.031286006183386,
  77.79269286293427,
  8.038905602635422,
  10.092418434217018,
  54.50888425602657,
  19.115829229354006,
  16.307896317582237,
  23.11646736866344,
  11.835805509931365,
  58.4208122945276,
  68.00264592152232,
  89.3055734145209,
  27.101168332401908,
  155.97445267401434,
  29.75750298426049,
  23.34075377987811,
  36.744082775375,
  123.6318940719091,
  25.333469751902417,
  63.12972683790582,
  25.649069995635703,
  32.480386006698666,
  40.943026369701954,
  68.73748110048398,
  46.93382425767469,
  27.290876316285612,
  11.096736525128136,
  83.04987298005474,
  15.823966867468592,
  38.523956579812904,
  72.66613647423546,
  11.14958433004702,
  21.156984439916876,
  6.34966576766438,
  45.58479520067375,
  32.057942678391235,
  46.7559274063442,
  98.28327272431125,
  16.29242147955353,
  22.480303727191078,
  34.55987189377694,
  29.552415073491265,
  38.01892592085635,
  98.41053956294138,
  11.877403044262572,
  71.79875517686352,
  7.505570473516513,
  40.51240594244258,
  29.30391403623791,
  30.845117844725422,
  102.85009012718496,
  50.86417940135735,
  21.972454116436467,
  76.13169016056835,
  97.45184125151427,
  99.49262479041113,
  54.422195965783885,
  12.242722896037312,
  73.16502773717922,
  33.41847808438594,
  28.507048199199936,
  67.10307868019618,
  22.078330501946507,
  204.34060600189628,
  11.691726991217637,
  88.35785954543536,
  20.749933532233122,
  31.9053429174757,
  33.467916987392506,
  10.037342858862516,
  72.05016479621312,
  50.654780784956976,
  78.5131746117998,
  91.70320522875704,
  19.64524028333861,
  25.693339753124917,
  39.013703458430754,
  5.33230240839648,
  18.707554922753562,
  73.84301921763249,
  26.4735725319478,
  24.348380639024608,
  139.57120572624285,
  45.85875800464462,
  6.967077432692455,
  6.146612087806292,
  11.857680192540462,
  20.851026189626353,
  157.93452591762576,
  83.07569543024667,
  7.000377524745289,
  10.75131256110481,
  52.49332646719358,
  25.536282365750004,
  20.984048434332085,
  5.307618469395774,
  61.66760754352518,
  5.35130183963476,
  74.05510952219996,
  5.521593362303416,
  85.61279378116981,
  145.19718946965094,
  21.20477264098052,
  6.738170534903634,
  5.645642976416612,
  34.98525051499954,
  253.0693348104893,
  10.77022032463575,
  63.290962328003914,
  78.6864083019725,
  119.33723308582427,
  51.381496507523124,
  33.633728270848664,
  33.540758597402984,
  86.15963548236772,
  29.244205983852503,
  8.892929680989385,
  23.39271274291585,
  24.725519906211517,
  38.54723292202948,
  13.186597545117456,
  23.55388387601943,
  22.22218962295973,
  24.043688715304707,
  75.7002835845927,
  77.46936925556162,
  34.95529504158023,
  27.986491478131114,
  41.38583741093936,
  168.69557015698499,
  51.330263082445136,
  8.544375716606,
  55.96632684865462,
  40.311210427630584,
  152.40173820332296,
  102.69393423195073,
  6.018391537242518,
  8.441306876419134,
  25.826756138396977,
  59.06497186292303,
  23.8222860786381,
  14.688606950860617,
  23.69184112206612,
  75.21439034755198,
  16.68983889948522,
  33.13411226442511,
  25.74321488759006,
  35.77688236233212,
  27.37701095088852,
  14.962366695162334,
  11.029499055834753,
  11.33609026871668,
  29.2449815769993,
  19.89243378046452,
  29.70986797661356,
  37.60482869983124,
  53.9068681593587,
  18.73091024211425,
  17.945932179552194,
  23.25616034807022,
  69.75187109237203,
  88.28838559047345,
  12.168609014198669,
  24.5818564188825,
  135.97955763187986,
  34.857019511838644,
  25.118165812364182,
  25.983493960498468,
  64.57200957966609,
  7.31256762020256,
  50.8477883103327,
  149.12229114760015,
  91.12360221342247,
  114.18278604092498,
  18.266763271835366,
  100.59778892481017,
  53.131207121831174,
  5.8591868981474295,
  9.713898889591693,
  8.835226283605628,
  32.445699287235826,
  33.70785473836754,
  72.6733288660623,
  68.90181621460378,
  91.90797461317592,
  105.97438914514262,
  106.49761238431915,
  18.775832338881457,
  103.16810133161476,
  51.0456249303767,
  35.854274526286865,
  19.78811839232225,
  16.922345301407262,
  15.347760266040202,
  71.89215367167961,
  32.16294301542492,
  8.013323662949785,
  7.429325407003637,
  43.09599128429415,
  39.06964310511996,
  28.272739595426216,
  26.588887944434944,
  65.52310324642417,
  5.485901274879486,
  29.961672290030542,
  43.956022194676756,
  29.09513891203491,
  57.47637549785293,
  87.90738239882754,
  35.88854007995246,
  15.701588175643574,
  12.09254340776239,
  68.22088344017432,
  30.405941217924497,
  10.88535589570661,
  25.866681264108188,
  28.354262025418628,
  73.55908828875974,
  17.743854888325586,
  14.441413162295266,
  34.89158735205282,
  79.18490323745337,
  11.490147812174165,
  21.972413593526667,
  12.61022157469751,
  66.5434867530931,
  35.65521726021844,
  35.573676638594314,
  56.166800836335106,
  21.346736277777737,
  36.80926426349184,
  121.19892488632408,
  12.591058563287843,
  5.712362902679958,
  51.96120504166045,
  29.851696317108235,
  81.24886287541852,
  77.46826200203334,
  76.9738457016018,
  22.482486632821992,
  44.010552091686144,
  27.259774117504442,
  63.60943484805226,
  70.20885095903978,
  111.88123674925092,
  53.68296710106666,
  6.610824803963359,
  79.72302145645614,
  108.52097776614347,
  89.99058027164054,
  8.285196438592337,
  119.72463814425151,
  18.609120844397495,
  5.752835120853492,
  99.54640388560067,
  17.770747664219364,
  63.15756993429879,
  127.42564224494123,
  10.156559748925629,
  134.37339674596353,
  18.251358853754205,
  20.55666365653747,
  59.50585546202748,
  20.336207060461767,
  10.456722595089174,
  15.991617823358645,
  42.14789229497926,
  17.405793380259134,
  98.60968083612752,
  108.9609508028696,
  5.028105310650802,
  30.14464100216142,
  106.03182502988724,
  5.747683174569619,
  76.75675658404526,
  87.13427072108821,
  13.195117341786863,
  7.487228825224888,
  144.28260138806618,
  106.84787936914017,
  30.539046495257235,
  78.0191970641961,
  8.126822988204411,
  46.40494212189582,
  24.75879695184923,
  40.042964297482555,
  5.037026685636164,
  110.48223772846353,
  50.59145327444017,
  57.77892632496538,
  249.2546921113635,
  40.09054385555049,
  8.781104630583034,
  52.3050205672632,
  52.30861883394773,
  5.444064360118548,
  40.88697821082237,
  48.99971081339312,
  8.892545279309934,
  16.723446677254188,
  7.64798764561481,
  36.74418018413055,
  53.06731000634828,
  32.28435134531958,
  19.703073959517653,
  42.21498596844757,
  9.586577701233379,
  24.737370069913993,
  39.17806036478852,
  73.37350670264945,
  68.00743471370616,
  44.8090884229112,
  39.91812550577923,
  6.426381676063272,
  13.947752123306596,
  41.817917884904254,
  24.38062669217548,
  128.06747247482542,
  72.83960768606997,
  56.55022160234503,
  73.72885583509289,
  8.30394320840856,
  19.761441519044602,
  6.689705923744102,
  14.975325230820555,
  10.778495354050087,
  128.2785998587404,
  21.765739383088576,
  107.96631902737053,
  53.04970611270958,
  134.25320360569648,
  38.873360476262484,
  30.470622215831224,
  19.577423782356846,
  75.17296550079249,
  25.686821762365195,
  28.795067882788995,
  5.828100233270641,
  81.1293065986106,
  40.0129442549944,
  24.405651247752306,
  5.900633080357883,
  70.02639580439629,
  15.202168508948457,
  22.815517438430888,
  8.233563763525833,
  80.82560617130285,
  40.01493574352517,
  9.063632043623945,
  19.78912551333547,
  17.175829002480473,
  82.66097900120657,
  18.96069815747877,
  32.290137914057965,
  11.484500112572016,
  31.90179012572043,
  100.50298998192858,
  117.73661364747552,
  123.43148312259714,
  58.21147976787574,
  8.482774726718528,
  5.738278628476706,
  10.641417115782785,
  79.45292133932595,
  63.395531297118666,
  25.87951917415535,
  15.739635204042116,
  12.02412149393823,
  5.905657727872323,
  49.93939095351203,
  73.43586526424619,
  7.346147092678332,
  19.996004017435702,
  16.297166521478413,
  68.10973712466347,
  20.7791888913556,
  74.26908225524363,
  21.669008205573935,
  90.67294036640165,
  8.770035245608964,
  27.207922160090327,
  59.06418362688859,
  5.4602714288084755,
  80.74891449325989,
  53.734605509756385,
  25.396963004938815,
  87.4023107488662,
  13.133148347753895,
  57.989542583577425,
  24.321543029841674,
  34.080013057012046,
  23.032670177587498,
  60.782494176893785,
  13.918222930219114,
  11.078964878958166,
  43.727566256015535,
  59.00982334483704,
  64.19268576754467,
  42.99093540147429,
  71.72355378448395,
  44.66086994206202,
  47.29178287896077,
  10.85074189734108,
  52.34692814906521,
  30.24969039753773,
  21.330818493425465,
  76.1754792877278,
  50.3165533086051,
  12.310247719617765,
  91.16205018518662,
  29.652968860364304,
  17.817151456560303,
  55.536675217166824,
  112.96779596297503,
  81.20440966609611,
  5.738579593055229,
  55.20565650611302,
  6.889970630558052,
  10.791198613732218,
  29.837414806711134,
  20.379709612445208,
  41.44132575692947,
  41.72337927013647])

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

([3492620.530665972,
  3582442.0812735655,
  3585565.319624986,
  3593032.2745178808,
  3611177.580869366,
  3614108.073302657,
  3616160.892630776,
  3693601.88804027,
  3732662.441962715,
  3803672.662098161,
  3856413.315365188,
  3924180.8969368823,
  3946479.039064631,
  3947820.0847265646,
  3948395.2593749003,
  3982932.808839554,
  3985380.2526885793,
  4004219.513812185,
  4004237.410135303,
  4004249.4211029704,
  4004315.625,
  4004326.5519919926,
  4032283.745416332,
  4048985.846178447,
  4063593.6444698013,
  4065810.262041953,
  4070489.615134211,
  4071332.619565964,
  4103654.834075178,
  4108865.1873446065,
  4111239.255548117,
  4115083.610033039,
  4140226.5625,
  4142648.6840628483,
  4157467.665319404,
  4158065.6581533193,
  4176712.3629468787,
  4176718.0829887926,
  4183423.4375,
  4183431.5333238444,
  4207255.635936901,
  4230417.674915464,
  4240376.5625,
  4245120.833868845,
  4246034.375,
  4278065.625,
  4278338.187965376,
  4296439.785176602,
  4314606.67847943,
  4337496.875,
  4338624.447477073,
  4344571.082869163,
  4368767.926705537,
  4369570.858679403,
  4370025.0,
  4370050.406730455,
  4370078.468745263,
  4382167.1875,
  4385253.618855921,
  4394957.270319626,
  4396710.300909533,
  4400208.5068542315,
  4401721.317354548,
  4403828.125,
  4404592.014905767,
  4405027.614335727,
  4405306.25,
  4411033.476109745,
  4432089.0625,
  4444100.594728054,
  4444638.828062689,
  4444921.463977997,
  4444945.3125,
  4445785.816980325,
  4453308.859898969,
  4467647.402255139,
  4467762.413638079,
  4507940.901302447,
  4513187.907767356,
  4513339.0625,
  4515787.049710502,
  4528356.492906535,
  4541787.5,
  4544634.458786175,
  4545223.42308966,
  4545577.350766959,
  4547352.408610105,
  4547697.608782832,
  4558034.321257396,
  4558949.32115245,
  4559228.836067757,
  4566832.8125,
  4566995.3125,
  4567151.0378156435,
  4567252.79983913,
  4567267.586596807,
  4567650.0,
  4567653.125,
  4568517.733738784,
  4568700.115373324,
  4568706.25,
  4588297.150567199,
  4588322.324956934,
  4606041.277647058,
  4607407.272407189,
  4609157.94107543,
  4622816.727336632,
  4624307.8125,
  4631048.34565392,
  4638880.898431857,
  4661785.110310686,
  4662986.689381507,
  4692118.131158823,
  4700800.0,
  4708326.0583390845,
  4715356.25,
  4718447.077253479,
  4718652.705737127,
  4723722.554191789,
  4738390.625,
  4746995.3125,
  4753088.111919016,
  4761115.5090038525,
  4761337.427545323,
  4762190.138718305,
  4771840.625,
  4778800.644308123,
  4780726.5625,
  4792370.3125,
  4792425.67511725,
  4793339.0625,
  4793453.402516461,
  4793502.388632653,
  4804800.0,
  4810887.5,
  4817915.314200792,
  4823158.910894705,
  4843160.150324343,
  4850747.3125097,
  4852696.875,
  4870146.727989322,
  4874169.022973611,
  4874746.601720019,
  4885514.8859380605,
  4889250.432388543,
  4900423.4375,
  4913836.337753746,
  4920548.32970231,
  4924215.132579058,
  4936441.705124193,
  4936867.781244053,
  4971395.955267249,
  4997933.071737891,
  5004409.9010975715,
  5006145.683023668,
  5014930.233848478,
  5015036.574655041,
  5016085.9375,
  5041387.378949439,
  5043382.149309411,
  5047707.686714367,
  5059738.900431105,
  5063029.740341798,
  5072741.591584372,
  5074093.582756976,
  5074668.611125128,
  5074720.107495489,
  5083679.904929095,
  5087359.006914243,
  5088588.141981034,
  5090221.2414149875,
  5090597.138631836,
  5090893.75,
  5093431.25,
  5098627.692084908,
  5104219.234343314,
  5105066.562350477,
  5110621.086165099,
  5115426.383465704,
  5144596.4875148125,
  5147997.910924866,
  5149588.509453688,
  5151425.0,
  5168914.855442729,
  5173954.6875,
  5202698.4375,
  5202700.0,
  5202898.840043367,
  5202918.736277276,
  5205207.531347487,
  5215830.12677419,
  5246987.987188586,
  5305675.148945387,
  5340362.5,
  5341082.8125,
  5342370.3125,
  5344550.0,
  5356625.0,
  5360156.113141704,
  5386041.450878613,
  5393534.375,
  5393540.545964344,
  5398086.782300612,
  5419719.49582887,
  5429926.08184323,
  5440662.895200803,
  5455821.875,
  5479009.375,
  5479677.034038522,
  5523718.320207222,
  5561756.25,
  5568504.477091591,
  5575988.719570688,
  5601538.788035155,
  5601559.112705029,
  5614991.1128464,
  5631440.625,
  5633206.119030191,
  5633585.9375,
  5633979.6875,
  5677488.938883729,
  5678219.195279469,
  5682221.5841641305,
  5690881.25,
  5707609.81274805,
  5731804.6875,
  5742124.622823982,
  5774245.3125,
  5777536.395349934,
  5782454.376485919,
  5792843.969078138,
  5797463.898883178,
  5822959.356834526,
  5826571.875,
  5833190.984185385,
  5833355.207311864,
  5835251.343775273,
  5835509.922134748,
  5835565.932978031,
  5875859.590684232,
  5876278.717555765,
  5892808.722669422,
  5903393.025912761,
  5911587.5,
  5956122.610491333,
  5964578.125,
  5981715.560970378,
  5990124.825824275,
  6041239.131424672,
  6066414.365561938,
  6137851.095705542,
  6139587.3817310715,
  6157208.400800219,
  6165635.91745936,
  6190257.4107712675,
  6194263.533878264,
  6195880.800181757,
  6211297.072798556,
  6211332.096785197,
  6211990.438800751,
  6214761.580308128,
  6234113.528688682,
  6238683.020853411,
  6243307.8125,
  6253770.3125,
  6284650.0,
  6289653.125,
  6291201.004408521,
  6291447.750224567,
  6299438.434041375,
  6324689.962954273,
  6347230.277480394,
  6364087.5,
  6364272.317130318,
  6384127.970561172,
  6384132.858884388,
  6384133.642231257,
  6391434.375,
  6391584.357495111,
  6391587.1726484755,
  6391875.356729987,
  6391925.323242055,
  6392496.988696254,
  6392859.992813377,
  6392877.7097100215,
  6398806.074027216,
  6411259.5695956135,
  6430553.318833017,
  6432774.412347517,
  6433266.607067543,
  6433924.773417076,
  6442751.922437774,
  6445746.875,
  6451976.251533856,
  6454548.4375,
  6454566.611227143,
  6455097.499999296,
  6456683.77042303,
  6464203.234935578,
  6464219.293255772,
  6465548.4375,
  6466196.555328,
  6466489.722449631,
  6466492.178129724,
  6467284.375,
  6472930.053967671,
  6475051.789275058,
  6475484.375,
  6475532.1163698565,
  6477843.628759042,
  6478815.625,
  6479753.125,
  6492304.6875,
  6492357.8125,
  6496264.0625,
  6498062.5,
  6508464.0625,
  6511036.838520094,
  6511122.878825806,
  6516896.834101093,
  6517370.3125,
  6517588.563651572,
  6523037.357605769,
  6529828.133206868,
  6535794.850179621,
  6536899.519315411,
  6540122.982247955,
  6540238.93716455,
  6540271.531899382,
  6541496.523669068,
  6552219.496800933,
  6552865.625,
  6552875.193282497,
  6556293.117438827,
  6577967.1875,
  6592498.127299166,
  6606261.908734295,
  6607526.446309192,
  6615879.514238279,
  6633326.17841807,
  6639489.156529045,
  6645651.44032431,
  6648528.125,
  6653554.6875,
  6655173.55230372,
  6655863.3431097735,
  6658726.443185236,
  6666032.000636315,
  6680556.022608152,
  6693474.988524131,
  6694570.3125,
  6699944.184130341,
  6703962.396695631,
  6704160.086962452,
  6704943.384933773,
  6721360.531430565,
  6731435.725928807,
  6732862.69270501,
  6742396.875,
  6750157.069466202,
  6751917.165264096,
  6752992.1875,
  6757478.125,
  6764751.5625,
  6765453.823780311,
  6767496.875,
  6795340.220823396,
  6799912.5,
  6816396.705970912,
  6821595.3125,
  6822638.597248429,
  6822917.611695682,
  6825792.1875,
  6844693.75,
  6852496.289600013,
  6855094.333095724,
  6855462.326195688,
  6855793.707262602,
  6855918.749325935,
  6856098.6254318245,
  6856279.586205223,
  6860275.0,
  6864036.312593633,
  6864861.337947928,
  6865712.5,
  6879417.591902216,
  6882189.66631164,
  6886814.092133719,
  6910634.375,
  6911928.027049641,
  6912199.444276428,
  6914246.355529393,
  6914421.875,
  6915046.663080895,
  6917353.679171382,
  6918871.875,
  6932687.224579017,
  6936685.234625235,
  6938915.753499886,
  6943875.271528397,
  6965887.5,
  6968153.317805092,
  6968426.5625,
  6968652.680989516,
  6969326.062819984,
  6969536.718524376,
  6969542.9068637015,
  6969830.83864261,
  6971314.380071097,
  6972119.961167927,
  6972120.722745279,
  6972320.3125,
  6976373.975332953,
  6982231.25,
  6986494.762410174,
  6987039.0625,
  6987539.979333383,
  7029499.667211107,
  7033105.880360608,
  7042340.625,
  7056637.597838184,
  7057428.125,
  7057438.777256733,
  7058022.902303999,
  7062387.314995076,
  7064029.6875,
  7064623.4375,
  7065604.621902001,
  7066596.855095408,
  7067418.75,
  7067548.651237835,
  7069804.6875,
  7070966.939162453,
  7076671.888951619,
  7081841.653132696,
  7083164.688490692,
  7085546.563091655,
  7090439.0625,
  7098923.4375,
  7100643.75,
  7100961.251880867,
  7101021.46857797,
  7101040.954066272,
  7101042.1875,
  7101070.304486073,
  7101967.1875,
  7119114.494449377,
  7119142.706738644,
  7119993.471642685,
  7120156.25,
  7120387.5,
  7120397.504038648,
  7120754.079197214,
  7122098.4375,
  7122749.094465627,
  7126332.8125,
  7126464.0625,
  7132846.875,
  7139929.6875,
  7143292.800249938,
  7143351.014692584,
  7145806.227342033,
  7145815.967921798,
  7146032.8125,
  7146334.7647872185,
  7146368.75,
  7146378.458172594,
  7146648.4375,
  7146990.626153024,
  7147481.904331879,
  7147978.171787495,
  7150636.668999937,
  7153875.0,
  7154065.625,
  7154959.598534582,
  7155929.6875,
  7159628.979901439,
  7214795.150746081,
  7215348.250305075,
  7223568.466567843,
  7223722.411912227,
  7231227.361970977,
  7236016.491392838,
  7239184.120389741,
  7271118.0936836535,
  7271159.375,
  7272755.62827271,
  7272766.169483935,
  7273498.859188631,
  7274816.197291385,
  7275573.065958184,
  7288536.183781996,
  7290547.95598082,
  7290587.24506547,
  7292195.15796612,
  7296486.459053229,
  7296558.822409572,
  7301244.276066935,
  7302043.75,
  7302115.896212602,
  7326172.3981050905,
  7326918.511018772,
  7327214.0625,
  7330618.75,
  7347440.095602483,
  7352156.121754101,
  7357736.494890946,
  7367009.34389197,
  7367735.221630025,
  7368079.446083401,
  7368216.153364367,
  7373823.4375,
  7374358.449346021,
  7375007.8125,
  7376081.066671932,
  7383918.75,
  7386906.20501683,
  7397754.6875,
  7401430.681537078,
  7427151.233944612,
  7428040.397663485,
  7428450.0,
  7428722.885781166,
  7429996.55943085,
  7431465.195108443,
  7447559.097592855,
  7447602.757488114,
  7470672.5777004585,
  7481824.845586427,
  7483934.347211462,
  7487076.5625,
  7488839.9415765675,
  7500006.25,
  7500836.395134269,
  7514646.24465726,
  7517896.189150857,
  7522315.458061876,
  7530896.607276471,
  7533835.239396883,
  7549100.281566813,
  7559627.97182353,
  7565146.83661196,
  7568092.264301567,
  7570412.432027111,
  7572172.26159283,
  7572197.711950465,
  7572737.893777727,
  7576065.207715668,
  7576670.153965108,
  7586321.875,
  7588195.251707606,
  7595362.505461409,
  7595401.498747828,
  7596086.066848376,
  7612506.443378305,
  7618515.39118069,
  7619398.139725847,
  7653207.162622215,
  7660807.8125,
  7676146.044969409,
  7683127.238787973,
  7695823.09776887,
  7703776.732833039,
  7706629.749099611,
  7709252.139933984,
  7725446.6373939235,
  7726557.8125,
  7734461.113694955,
  7734773.867217522,
  7735946.731765266,
  7738141.560988794,
  7740449.117850486,
  7745000.0,
  7745744.558603443,
  7745809.266897999,
  7746845.3125,
  7746900.66416836,
  7750683.946292246,
  7750848.4375,
  7751427.15533355,
  7753053.5451840265,
  7753597.785976197,
  7753674.897091295,
  7756391.152311827,
  7761139.0625,
  7770390.625,
  7770892.1875,
  7776121.875,
  7776493.75,
  7776756.988846817,
  7778085.9375,
  7778681.25,
  7784767.258994736,
  7784774.440559127,
  7786667.625403491,
  7786679.790152364,
  7789170.3125,
  7793947.635603636,
  7822791.632372629,
  7838968.242121083,
  7839329.6875,
  7843628.125,
  7844465.612131976,
  7845749.106219658,
  7849756.25,
  7852133.362998628,
  7852907.548539679,
  7857587.54657736,
  7880150.0,
  7883243.75,
  7892205.61200743,
  7898414.000005828,
  7898436.123579119,
  7898706.790486198,
  7899232.685738603,
  7904143.75,
  7905120.151713699,
  7916242.036869683,
  7918820.3160350695,
  7918825.0,
  7919669.331555442,
  7919935.799800518,
  7920408.145476777,
  7920415.625,
  7926922.095130529,
  7927137.682291421,
  7935427.341019714,
  7938127.131000082,
  7939946.569443002,
  7940570.25281298,
  7944178.125,
  7949877.715906706,
  7956116.897001985,
  7982109.375,
  7983366.395240506,
  7985389.1502094725,
  7996799.295779672,
  7997866.8535613185,
  8034257.8125,
  8034809.375,
  8034994.515342919,
  8036389.790021534,
  8038466.381966146,
  8039290.557657879,
  8041657.8125,
  8042293.066106248,
  8042412.3953181375,
  8043581.25,
  8044189.808884249,
  8047664.878539154,
  8055109.375,
  8063413.854806508,
  8066609.37238781,
  8075780.392069062,
  8077238.2149775075,
  8078544.450445712,
  8079368.333808584,
  8079539.648461029,
  8079745.3125,
  8080160.9375,
  8086657.087519399,
  8087416.855439277,
  8096735.036823418,
  8105136.900707128,
  8108277.5630520005,
  8113895.07974848,
  8113996.398662962,
  8118676.8131856285,
  8140610.9375,
  8153160.736983476,
  8155768.167216057,
  8156762.243277466,
  8156828.874551024,
  8159345.695481911,
  8159704.982094423,
  8159831.874398971,
  8160810.9375,
  8178712.5,
  8178891.6389876995,
  8179357.479943236,
  8183305.013713312,
  8186688.088251388,
  8205784.375,
  8245112.476645995,
  8247057.4976698505,
  8249771.352389682,
  8275849.775021633,
  8275957.662306401,
  8276030.980215641,
  8287891.443281613,
  8296356.25,
  8305550.15672711,
  8322825.654079661,
  8343010.827270233,
  8355466.45611102,
  8365987.469621887,
  8366748.4375,
  8383557.393473591,
  8387384.375,
  8403367.818729887,
  8404104.118490715,
  8408797.275258327,
  8426721.997055119,
  8426736.750988714,
  8434635.9375,
  8436923.13259699,
  8447150.862155644,
  8447308.546684612,
  8448869.037115455,
  8461332.586004118,
  8467905.518682485,
  8469360.4143886,
  8470163.745655095,
  8470257.2851717,
  8471307.824198764,
  8471408.97647027,
  8475051.274818866,
  8478850.0,
  8489110.685261818,
  8490554.5176276,
  8491959.375,
  8503795.40132972,
  8504709.375,
  8505718.782695398,
  8509527.431921665,
  8509702.743742788,
  8509704.696343774,
  8509709.139474593,
  8511585.9375,
  8513908.901903262,
  8528454.6875,
  8529515.64879859,
  8530264.0625,
  8532033.919870757,
  8542618.07049368,
  8542937.812946761,
  8552277.870632626,
  8558145.150324287,
  8560827.576800201,
  8563912.689026097,
  8595963.441127613,
  8600698.4375,
  8604509.124597415,
  8616246.41874446,
  8617727.177632144,
  8716803.191746434,
  8718575.515214734,
  8743674.898386894,
  8744463.578823779,
  8801959.05638304,
  8802316.852881536,
  8816921.244708845,
  8835612.984189479,
  8836453.125,
  8842844.514679093,
  8844840.625,
  8878528.341871582,
  8881756.25,
  8882220.3125,
  8929487.5,
  8934526.871778646,
  8950714.0625,
  8953850.353456711,
  8957453.609139135,
  8959626.22739169,
  8959795.3125,
  8960596.875,
  8965852.975186393,
  8968903.782620886,
  8970583.53594613,
  8973733.97914366,
  8978084.569004256,
  8979235.9375,
  8981777.353465248,
  8987239.9605755,
  8989512.815359963,
  8997095.3125,
  8997666.92170722,
  8997679.792797567,
  8999127.642454976,
  8999901.43034099,
  8999902.37703619,
  9000360.867239624,
  9000675.623208854,
  9012753.090119759,
  9016464.0625,
  9032564.0625,
  9033107.631302927,
  9067721.875,
  9070687.636026215,
  9072659.375,
  9074527.868409682,
  9085459.375,
  9092251.565174889,
  9092685.9375,
  9093028.682758693,
  9108881.218591005,
  9108893.7380893,
  9109812.5,
  9110204.6875,
  9110928.988959128,
  9111068.371907692,
  9111504.6875,
  9113209.375,
  9122298.252772884,
  9122722.512120308,
  9124526.5625,
  9125286.98357668,
  9125445.3125,
  9126162.5,
  9126802.206854604,
  9126842.1875,
  9129084.375,
  9131002.156021588,
  9131031.25,
  9131746.456449954,
  9134075.0,
  9134684.863656877,
  9138520.3125,
  9139216.227082277,
  9139795.085988956,
  9168178.125,
  9365160.318142092,
  9368031.25,
  9453001.5625],
 [15.614769966327133,
  37.1155661752404,
  13.892675712334887,
  20.241280469855834,
  21.063886871005202,
  11.561390086842907,
  62.59705518981089,
  8.466572728202882,
  14.724818400001537,
  36.736823126277464,
  12.696594241595042,
  13.747330507512231,
  63.16803435687988,
  65.12943017927766,
  55.52803137916615,
  39.473196921062566,
  23.650138656438582,
  28.85526830724879,
  16.71224874559791,
  20.740581028718143,
  97.29568450537086,
  59.028458918001974,
  29.59376884549594,
  9.255361172765275,
  54.49549892691492,
  66.94533111533197,
  9.846912369276689,
  10.838960265892576,
  137.02952433588655,
  100.59065489273542,
  15.7321768071244,
  10.11101478338409,
  51.32999251980287,
  5.045723598982221,
  5.306078524144217,
  24.92533260956835,
  75.10472290907104,
  12.890788075776735,
  49.0619809609241,
  6.5324412271671255,
  8.86106871565045,
  11.985834640184864,
  54.03718192212205,
  22.259530610820686,
  43.73564648322419,
  53.84838796365811,
  18.425042132310182,
  6.606676085773855,
  6.852174197851404,
  55.19922074208506,
  8.904173097776319,
  59.68584064725847,
  25.2829581277543,
  17.26891637906644,
  39.19146190815111,
  11.82911222496202,
  40.192983023911104,
  29.219065600560405,
  32.33581906739688,
  6.985867298033824,
  12.966902049470034,
  18.46815003973684,
  13.170893515533804,
  37.315019188323355,
  22.09894389414234,
  85.12608820679435,
  57.78325019931273,
  12.351955590685527,
  30.321818163672166,
  9.89057408230384,
  8.166375418081824,
  22.429128646269906,
  69.10336682472672,
  60.832179637845826,
  55.4057796795294,
  7.8057588135914315,
  22.150615242967262,
  9.576496073515967,
  71.06586551816848,
  104.2069923174784,
  14.525028134304282,
  15.955914680695997,
  44.30549979451328,
  46.906871472977436,
  14.631394093849668,
  19.133685562949957,
  12.770058374115244,
  18.74013905033442,
  6.260432171533588,
  12.861534608013159,
  18.582180018618086,
  79.32597696451268,
  49.788922400085085,
  13.258946002628806,
  36.34987955035869,
  15.453648470308137,
  56.01513576025013,
  56.57802468723617,
  12.104289185775832,
  191.72722170656115,
  45.148517697050735,
  25.0560839378684,
  16.98738677150621,
  14.115762999429755,
  10.376843506954224,
  15.455322342465514,
  21.899914704531447,
  29.24345675181684,
  26.831549118509326,
  70.43448575113271,
  30.20493553074924,
  18.868525731359064,
  20.384218229991372,
  42.37443682306643,
  21.451693058278334,
  55.49891763335433,
  33.83468312838282,
  8.07664690542092,
  28.983502769845,
  35.12579419774115,
  42.749059836789236,
  6.5874928950143605,
  45.14358337721294,
  14.67235505187391,
  30.103936395629724,
  106.2046557148694,
  51.11090250077219,
  32.42336091918229,
  60.252942815786966,
  24.920637316337665,
  95.31749870474144,
  9.577626043854025,
  19.42010652968349,
  36.637775658002134,
  43.378501388559236,
  92.87806174940461,
  38.23834302070315,
  19.0382642717784,
  59.756636032487116,
  72.49994725636412,
  25.92369526389607,
  11.275175673088256,
  24.11265099090807,
  19.399358338450764,
  6.134248414886759,
  71.62383245784048,
  11.816809855207078,
  9.241887621316186,
  21.569040032597886,
  13.374575119250867,
  6.888129305730314,
  30.191015873992004,
  54.32482614634452,
  17.41430837926818,
  6.1274970187611535,
  7.272878724738667,
  18.640481526646383,
  57.77403819906828,
  178.2087407296472,
  21.384764144749912,
  120.23896322770342,
  14.357779945648911,
  15.146816720241523,
  41.90352234949461,
  25.60882627327911,
  45.09960575412909,
  25.04507587472712,
  114.83155242167213,
  6.692180620390819,
  8.657987590190976,
  20.25656181925672,
  16.051525015351807,
  47.73468606226955,
  30.37558368317511,
  74.51203722740736,
  138.58255717857824,
  20.329349015973744,
  11.423477200525113,
  25.429523495426874,
  7.4688558278231545,
  8.922594759955786,
  8.903803042217083,
  130.81485424892517,
  18.124077186263115,
  74.62908220448168,
  52.60595312984717,
  54.20316714004254,
  78.15027000541264,
  38.106281991515786,
  48.60783821814187,
  11.891545942018722,
  17.19260620285983,
  83.93993021879407,
  81.24562335492317,
  187.6732339990674,
  41.56161234324317,
  43.09521258987494,
  50.701982387489124,
  23.800246406813173,
  13.804847946441853,
  72.18935478335115,
  11.006376061748195,
  19.565433417559536,
  26.73720479114686,
  21.68825414086999,
  5.974577115263925,
  31.790842810438928,
  42.22744027248871,
  17.105538022898187,
  49.0002138704361,
  59.818196732797844,
  83.35148476766825,
  23.217442306887218,
  24.082491556793475,
  30.41818903465886,
  104.91605934163945,
  53.08921356835019,
  26.18003708235912,
  32.404860511710915,
  91.72925653536242,
  26.29959689888458,
  16.962067089454333,
  33.48082550129065,
  118.15282494567504,
  167.96461690623929,
  46.07797365178146,
  22.812563466394543,
  219.5682676981942,
  23.184456195895123,
  23.611064226445656,
  5.102792608902199,
  73.22300232501627,
  64.73154076261136,
  54.55749497962934,
  5.354382169414058,
  17.37797619653832,
  32.92329519364646,
  17.537392266493853,
  11.645821800988376,
  15.832678224039043,
  84.52557778446868,
  12.92735129875925,
  10.157585505504422,
  74.70163628972608,
  72.22133960407947,
  57.6178778241228,
  11.029077861433978,
  25.503037441037886,
  5.819097764764022,
  16.23783808909282,
  21.83086440518341,
  123.17392115032798,
  29.616898679535318,
  9.314514975077797,
  5.328341815390537,
  87.63899685968752,
  31.701758846490645,
  42.184393593515225,
  46.748058475859374,
  52.81879919119535,
  93.7476097062621,
  21.16231721422618,
  16.979431698592435,
  29.267180313387556,
  39.28886108362691,
  33.427321198739264,
  139.11216565050213,
  13.169713912713586,
  5.433021377535052,
  151.63388231503578,
  10.229716872252059,
  7.285684956647565,
  103.19273936553219,
  22.533919494451325,
  34.79268730307008,
  11.200990938961816,
  36.73843471608245,
  44.58660611735263,
  39.732487474268446,
  8.120736012487809,
  121.77130222554614,
  33.01803531676436,
  6.281670738032824,
  7.843096545547782,
  22.38190456009417,
  110.92680870249507,
  15.738991499722053,
  20.820211172519667,
  20.661121200768935,
  57.5342239659596,
  137.42663219062536,
  16.517269056668876,
  74.04171029572589,
  6.702972943278148,
  211.18472754692783,
  8.003660565241498,
  12.268758038421767,
  25.647584810419737,
  9.426122626518538,
  17.514523329617575,
  57.18105031931751,
  16.059723036716527,
  40.46699357846643,
  28.34084889602217,
  32.77556313244104,
  104.65795209437849,
  22.429708683220895,
  47.80507906155397,
  17.686752789101448,
  9.230948792139776,
  93.70197663950312,
  85.4058752709449,
  77.00890290063982,
  52.89874182584843,
  37.87262801940629,
  144.54841721466283,
  51.08606870346067,
  21.209218939866595,
  29.298845925913593,
  102.99168130436092,
  84.52575992909023,
  105.77023454626985,
  25.709602800863735,
  121.3541794723686,
  29.088259994771544,
  8.047751647976714,
  58.503977207232765,
  36.78543761901659,
  33.70050751915862,
  5.909382543038365,
  18.8322551422548,
  69.08253167817767,
  35.46836429076977,
  136.90944721156797,
  53.60901716560561,
  103.03247809672973,
  139.73376895827113,
  25.306971868091786,
  21.06640343366051,
  9.02646639453157,
  98.94484438452403,
  26.32317251459509,
  33.971344917693855,
  64.10678228944015,
  125.10978370975555,
  107.48054325867645,
  26.597075968483608,
  24.32661660305525,
  44.67003024912576,
  11.1098165331589,
  89.02173642931777,
  16.902122266950656,
  75.82484670715624,
  83.06291495837716,
  5.035445093365805,
  48.40588333561513,
  10.812362896600913,
  9.499840748353096,
  31.789836924689762,
  105.8424488909115,
  14.857428517569838,
  87.8587700337127,
  42.756548569085425,
  61.44966277342754,
  12.429274468526797,
  78.04413831728195,
  5.592068706900726,
  101.73320624257796,
  84.15046236160833,
  50.23066155073736,
  18.923879426664016,
  11.862489629562008,
  37.06145111980516,
  31.499396268797277,
  6.975625174082379,
  8.62103577198507,
  14.317997139623957,
  14.781078694066746,
  34.346442066486944,
  22.689161776299247,
  7.0950765072860325,
  98.60640424052767,
  24.509887633028924,
  5.975054477249125,
  47.01351712067179,
  21.999729449569212,
  43.53668709539912,
  11.174344295336557,
  51.30840586204079,
  26.727349972655723,
  94.39102148331817,
  136.1257846525645,
  121.6408165983493,
  68.10479088366853,
  12.120544644856292,
  34.20277965063915,
  24.069865312011032,
  19.68285087389928,
  15.455005834525268,
  26.02364905207331,
  63.53270967768129,
  11.443067776574193,
  44.031286006183386,
  77.79269286293427,
  8.038905602635422,
  10.092418434217018,
  54.50888425602657,
  19.115829229354006,
  16.307896317582237,
  23.11646736866344,
  11.835805509931365,
  58.4208122945276,
  68.00264592152232,
  89.3055734145209,
  27.101168332401908,
  155.97445267401434,
  29.75750298426049,
  23.34075377987811,
  36.744082775375,
  123.6318940719091,
  25.333469751902417,
  63.12972683790582,
  25.649069995635703,
  32.480386006698666,
  40.943026369701954,
  68.73748110048398,
  46.93382425767469,
  27.290876316285612,
  11.096736525128136,
  83.04987298005474,
  15.823966867468592,
  38.523956579812904,
  72.66613647423546,
  11.14958433004702,
  21.156984439916876,
  6.34966576766438,
  45.58479520067375,
  32.057942678391235,
  46.7559274063442,
  98.28327272431125,
  16.29242147955353,
  22.480303727191078,
  34.55987189377694,
  29.552415073491265,
  38.01892592085635,
  98.41053956294138,
  11.877403044262572,
  71.79875517686352,
  7.505570473516513,
  40.51240594244258,
  29.30391403623791,
  30.845117844725422,
  102.85009012718496,
  50.86417940135735,
  21.972454116436467,
  76.13169016056835,
  97.45184125151427,
  99.49262479041113,
  54.422195965783885,
  12.242722896037312,
  73.16502773717922,
  33.41847808438594,
  28.507048199199936,
  67.10307868019618,
  22.078330501946507,
  204.34060600189628,
  11.691726991217637,
  88.35785954543536,
  20.749933532233122,
  31.9053429174757,
  33.467916987392506,
  10.037342858862516,
  72.05016479621312,
  50.654780784956976,
  78.5131746117998,
  91.70320522875704,
  19.64524028333861,
  25.693339753124917,
  39.013703458430754,
  5.33230240839648,
  18.707554922753562,
  73.84301921763249,
  26.4735725319478,
  24.348380639024608,
  139.57120572624285,
  45.85875800464462,
  6.967077432692455,
  6.146612087806292,
  11.857680192540462,
  20.851026189626353,
  157.93452591762576,
  83.07569543024667,
  7.000377524745289,
  10.75131256110481,
  52.49332646719358,
  25.536282365750004,
  20.984048434332085,
  5.307618469395774,
  61.66760754352518,
  5.35130183963476,
  74.05510952219996,
  5.521593362303416,
  85.61279378116981,
  145.19718946965094,
  21.20477264098052,
  6.738170534903634,
  5.645642976416612,
  34.98525051499954,
  253.0693348104893,
  10.77022032463575,
  63.290962328003914,
  78.6864083019725,
  119.33723308582427,
  51.381496507523124,
  33.633728270848664,
  33.540758597402984,
  86.15963548236772,
  29.244205983852503,
  8.892929680989385,
  23.39271274291585,
  24.725519906211517,
  38.54723292202948,
  13.186597545117456,
  23.55388387601943,
  22.22218962295973,
  24.043688715304707,
  75.7002835845927,
  77.46936925556162,
  34.95529504158023,
  27.986491478131114,
  41.38583741093936,
  168.69557015698499,
  51.330263082445136,
  8.544375716606,
  55.96632684865462,
  40.311210427630584,
  152.40173820332296,
  102.69393423195073,
  6.018391537242518,
  8.441306876419134,
  25.826756138396977,
  59.06497186292303,
  23.8222860786381,
  14.688606950860617,
  23.69184112206612,
  75.21439034755198,
  16.68983889948522,
  33.13411226442511,
  25.74321488759006,
  35.77688236233212,
  27.37701095088852,
  14.962366695162334,
  11.029499055834753,
  11.33609026871668,
  29.2449815769993,
  19.89243378046452,
  29.70986797661356,
  37.60482869983124,
  53.9068681593587,
  18.73091024211425,
  17.945932179552194,
  23.25616034807022,
  69.75187109237203,
  88.28838559047345,
  12.168609014198669,
  24.5818564188825,
  135.97955763187986,
  34.857019511838644,
  25.118165812364182,
  25.983493960498468,
  64.57200957966609,
  7.31256762020256,
  50.8477883103327,
  149.12229114760015,
  91.12360221342247,
  114.18278604092498,
  18.266763271835366,
  100.59778892481017,
  53.131207121831174,
  5.8591868981474295,
  9.713898889591693,
  8.835226283605628,
  32.445699287235826,
  33.70785473836754,
  72.6733288660623,
  68.90181621460378,
  91.90797461317592,
  105.97438914514262,
  106.49761238431915,
  18.775832338881457,
  103.16810133161476,
  51.0456249303767,
  35.854274526286865,
  19.78811839232225,
  16.922345301407262,
  15.347760266040202,
  71.89215367167961,
  32.16294301542492,
  8.013323662949785,
  7.429325407003637,
  43.09599128429415,
  39.06964310511996,
  28.272739595426216,
  26.588887944434944,
  65.52310324642417,
  5.485901274879486,
  29.961672290030542,
  43.956022194676756,
  29.09513891203491,
  57.47637549785293,
  87.90738239882754,
  35.88854007995246,
  15.701588175643574,
  12.09254340776239,
  68.22088344017432,
  30.405941217924497,
  10.88535589570661,
  25.866681264108188,
  28.354262025418628,
  73.55908828875974,
  17.743854888325586,
  14.441413162295266,
  34.89158735205282,
  79.18490323745337,
  11.490147812174165,
  21.972413593526667,
  12.61022157469751,
  66.5434867530931,
  35.65521726021844,
  35.573676638594314,
  56.166800836335106,
  21.346736277777737,
  36.80926426349184,
  121.19892488632408,
  12.591058563287843,
  5.712362902679958,
  51.96120504166045,
  29.851696317108235,
  81.24886287541852,
  77.46826200203334,
  76.9738457016018,
  22.482486632821992,
  44.010552091686144,
  27.259774117504442,
  63.60943484805226,
  70.20885095903978,
  111.88123674925092,
  53.68296710106666,
  6.610824803963359,
  79.72302145645614,
  108.52097776614347,
  89.99058027164054,
  8.285196438592337,
  119.72463814425151,
  18.609120844397495,
  5.752835120853492,
  99.54640388560067,
  17.770747664219364,
  63.15756993429879,
  127.42564224494123,
  10.156559748925629,
  134.37339674596353,
  18.251358853754205,
  20.55666365653747,
  59.50585546202748,
  20.336207060461767,
  10.456722595089174,
  15.991617823358645,
  42.14789229497926,
  17.405793380259134,
  98.60968083612752,
  108.9609508028696,
  5.028105310650802,
  30.14464100216142,
  106.03182502988724,
  5.747683174569619,
  76.75675658404526,
  87.13427072108821,
  13.195117341786863,
  7.487228825224888,
  144.28260138806618,
  106.84787936914017,
  30.539046495257235,
  78.0191970641961,
  8.126822988204411,
  46.40494212189582,
  24.75879695184923,
  40.042964297482555,
  5.037026685636164,
  110.48223772846353,
  50.59145327444017,
  57.77892632496538,
  249.2546921113635,
  40.09054385555049,
  8.781104630583034,
  52.3050205672632,
  52.30861883394773,
  5.444064360118548,
  40.88697821082237,
  48.99971081339312,
  8.892545279309934,
  16.723446677254188,
  7.64798764561481,
  36.74418018413055,
  53.06731000634828,
  32.28435134531958,
  19.703073959517653,
  42.21498596844757,
  9.586577701233379,
  24.737370069913993,
  39.17806036478852,
  73.37350670264945,
  68.00743471370616,
  44.8090884229112,
  39.91812550577923,
  6.426381676063272,
  13.947752123306596,
  41.817917884904254,
  24.38062669217548,
  128.06747247482542,
  72.83960768606997,
  56.55022160234503,
  73.72885583509289,
  8.30394320840856,
  19.761441519044602,
  6.689705923744102,
  14.975325230820555,
  10.778495354050087,
  128.2785998587404,
  21.765739383088576,
  107.96631902737053,
  53.04970611270958,
  134.25320360569648,
  38.873360476262484,
  30.470622215831224,
  19.577423782356846,
  75.17296550079249,
  25.686821762365195,
  28.795067882788995,
  5.828100233270641,
  81.1293065986106,
  40.0129442549944,
  24.405651247752306,
  5.900633080357883,
  70.02639580439629,
  15.202168508948457,
  22.815517438430888,
  8.233563763525833,
  80.82560617130285,
  40.01493574352517,
  9.063632043623945,
  19.78912551333547,
  17.175829002480473,
  82.66097900120657,
  18.96069815747877,
  32.290137914057965,
  11.484500112572016,
  31.90179012572043,
  100.50298998192858,
  117.73661364747552,
  123.43148312259714,
  58.21147976787574,
  8.482774726718528,
  5.738278628476706,
  10.641417115782785,
  79.45292133932595,
  63.395531297118666,
  25.87951917415535,
  15.739635204042116,
  12.02412149393823,
  5.905657727872323,
  49.93939095351203,
  73.43586526424619,
  7.346147092678332,
  19.996004017435702,
  16.297166521478413,
  68.10973712466347,
  20.7791888913556,
  74.26908225524363,
  21.669008205573935,
  90.67294036640165,
  8.770035245608964,
  27.207922160090327,
  59.06418362688859,
  5.4602714288084755,
  80.74891449325989,
  53.734605509756385,
  25.396963004938815,
  87.4023107488662,
  13.133148347753895,
  57.989542583577425,
  24.321543029841674,
  34.080013057012046,
  23.032670177587498,
  60.782494176893785,
  13.918222930219114,
  11.078964878958166,
  43.727566256015535,
  59.00982334483704,
  64.19268576754467,
  42.99093540147429,
  71.72355378448395,
  44.66086994206202,
  47.29178287896077,
  10.85074189734108,
  52.34692814906521,
  30.24969039753773,
  21.330818493425465,
  76.1754792877278,
  50.3165533086051,
  12.310247719617765,
  91.16205018518662,
  29.652968860364304,
  17.817151456560303,
  55.536675217166824,
  112.96779596297503,
  81.20440966609611,
  5.738579593055229,
  55.20565650611302,
  6.889970630558052,
  10.791198613732218,
  29.837414806711134,
  20.379709612445208,
  41.44132575692947,
  41.72337927013647])

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)