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

([4928112.4989503175,
  5036230.994222343,
  5073247.653318141,
  5157092.537210761,
  5191763.327334791,
  5300340.535763368,
  5317382.550383241,
  5321083.964715635,
  5355676.5625,
  5438311.969999445,
  5451291.820822522,
  6778494.16440898,
  6785040.20378968,
  6961798.4375,
  6961824.3864494115,
  6965671.875,
  6966445.3125,
  6996359.375,
  6998190.625,
  7042068.75,
  7067759.136083348,
  7085209.378132706,
  7139805.228617017,
  7160751.5625,
  7173569.216348684,
  7174431.809798662,
  7183017.109251836,
  7267945.3125,
  7270080.18207427,
  7277856.25,
  7297820.795334424,
  7308457.590299777,
  7325564.0625,
  7338440.044851796,
  7339704.23392173,
  7383595.006030557,
  7408210.089322634,
  7439252.530490686,
  7449281.994137895,
  7450372.632883947,
  7450440.53282894,
  7453479.6875,
  7460293.943414853,
  7461414.0625,
  7462551.5625,
  7464193.726921775,
  7464526.5625,
  7465321.813099671,
  7467959.614432444,
  7501741.049957041,
  7510363.752011988,
  7514546.034548267,
  7520335.081267191,
  7520371.272001179,
  7522658.281426924,
  7524769.998808046,
  7528399.2505654795,
  7529044.037948253,
  7541440.625,
  7546482.8125,
  7546489.0625,
  7547263.463943135,
  7553357.8791160025,
  7561035.9375,
  7577926.50891812,
  7589500.577856496,
  7589963.907652801,
  7591550.633801611,
  7592200.0,
  7592793.412688895,
  7594035.9375,
  7595455.167362315,
  7597749.088604033,
  7597774.573422763,
  7598457.582963947,
  7600586.047011132,
  7600948.586543177,
  7601541.609075746,
  7607356.241776778,
  7608885.987471846,
  7617818.75,
  7623774.413744111,
  7623906.25,
  7628007.522701366,
  7631127.952310165,
  7642312.933400372,
  7642623.4375,
  7646555.927471407,
  7659704.424472292,
  7660423.246718515,
  7661290.781656286,
  7667187.519684521,
  7669906.820943252,
  7670069.640706044,
  7676299.3879393805,
  7677092.784260475,
  7693406.840187491,
  7693458.300954373,
  7695421.661761349,
  7696093.095569686,
  7698765.346409517,
  7713493.75,
  7721913.727597748,
  7722556.980149264,
  7724259.525576042,
  7724515.863845837,
  7724843.75,
  7727261.023425729,
  7727271.356537071,
  7727415.625,
  7727877.194644847,
  7728186.571872057,
  7728515.680084369,
  7729054.871836707,
  7729200.301527607,
  7730979.611331726,
  7731413.585203696,
  7733510.9375,
  7734096.009670198,
  7740732.768387361,
  7741889.0625,
  7741920.126632115,
  7744231.594028677,
  7747357.8125,
  7747757.514116983,
  7749698.4375,
  7751104.6875,
  7752864.455034001,
  7756960.9375,
  7764938.664784911,
  7766651.5625,
  7767120.3125,
  7767438.239700456,
  7767926.5625,
  7768345.144835219,
  7768367.540090685,
  7768904.336399147,
  7768928.964429455,
  7768935.9375,
  7769650.0,
  7770121.04240849,
  7773571.440123174,
  7773998.411369783,
  7774192.110536828,
  7774225.0,
  7774447.749575863,
  7774453.125,
  7775539.0625,
  7777068.605323072,
  7777170.3125,
  7782460.9375,
  7782693.75,
  7783134.910419227,
  7783312.794751378,
  7783749.427500115,
  7784467.062865479,
  7787787.064645441,
  7795668.816698004,
  7795856.25,
  7796655.232118821,
  7797149.858760677,
  7802854.6875,
  7807122.39877029,
  7807973.4375,
  7826756.25,
  7826756.25,
  7842487.447523741,
  7844718.535868823,
  7851614.102293992,
  7865723.294864476,
  7870084.375,
  7874134.848562917,
  7877729.353606684,
  7878374.032340517,
  7890621.875,
  7902452.33385658,
  7921685.9375,
  7921687.5,
  7924925.46120395,
  7938648.198609276,
  7959710.47465506,
  7985421.918472736,
  7992271.766758991,
  7995504.696577904,
  7997695.11200486,
  8000735.9375,
  8023782.600891208,
  8026537.5,
  8026938.284448174,
  8027025.55993296,
  8027387.111576463,
  8027482.8125,
  8028878.228104227,
  8031037.455439818,
  8031115.537842666,
  8033304.051546997,
  8039536.7103334535,
  8043325.73759618,
  8044762.24497101,
  8045268.983490752,
  8045851.5625,
  8048845.415036554,
  8049198.030699028,
  8049636.389559176,
  8050270.3125,
  8050626.5625,
  8050758.249043236,
  8050821.624179483,
  8051181.404633995,
  8051940.625,
  8052812.131277288,
  8053056.706105952,
  8054601.5625,
  8057075.833765266,
  8063735.326451124,
  8066506.25,
  8076059.375,
  8076200.156499177,
  8078877.447056814,
  8079407.79223353,
  8079756.25,
  8081642.1875,
  8085040.697201559,
  8089092.171798846,
  8089559.375,
  8089950.870734111,
  8090286.385033847,
  8090292.187302922,
  8092241.566301165,
  8094739.0625,
  8104568.75,
  8122196.365773838,
  8124370.317114247,
  8145745.3125,
  8153418.084821994,
  8155963.30454557,
  8181686.231968025,
  8186329.011787595,
  8196593.75,
  8199386.539502798,
  8200132.462763735,
  8200161.615539064,
  8201437.270873369,
  8201957.7435368,
  8204474.977320864,
  8207768.75,
  8208186.36743079,
  8208777.336027159,
  8213970.3125,
  8218710.171837039,
  8220077.517073625,
  8220664.455636583,
  8220664.916431453,
  8233194.623744854,
  8240549.360318952,
  8242492.1875,
  8251937.5,
  8270687.059177668,
  8279017.572944678,
  8287757.8125,
  8294750.347112862,
  8297646.971576843,
  8300795.3125,
  8302445.316396724,
  8305570.725915888,
  8309044.852922054,
  8311505.633486733,
  8312116.216577061,
  8312432.117639927,
  8313599.753057187,
  8313651.069343612,
  8314070.3125,
  8314099.423960631,
  8315051.476407381,
  8315123.4375,
  8316479.6875,
  8316499.716329852,
  8317019.467708244,
  8317121.875,
  8317151.5625,
  8317175.684282461,
  8318581.16909885,
  8319086.0302196145,
  8319989.0625,
  8319990.625,
  8320053.968010188,
  8320442.108191726,
  8321211.738216368,
  8322905.739160612,
  8324197.45631981,
  8326223.3470792435,
  8327634.02880927,
  8340076.558805409,
  8344870.6528990865,
  8345062.327534265,
  8345539.155398426,
  8345935.26695746,
  8346270.3125,
  8348530.291463457,
  8350254.5892514,
  8351448.960669997,
  8352077.085393564,
  8353981.25,
  8354491.011043467,
  8354655.271274155,
  8354942.961162562,
  8357061.836866527,
  8359675.327558865,
  8361828.597066413,
  8362179.213305795,
  8362197.206458042,
  8362717.1875,
  8364066.463802592,
  8365800.331625924,
  8366057.262007171,
  8366451.360100323,
  8371052.421910109,
  8371646.3900579205,
  8372870.3125,
  8373288.062428359,
  8373307.968666479,
  8376871.875,
  8377524.114762574,
  8377525.984602657,
  8385687.1533857025,
  8392364.0625,
  8393432.418561727,
  8393451.884448973,
  8393855.74306448,
  8395354.010603042,
  8395506.202094821,
  8398071.654072808,
  8399175.584507551,
  8399606.085573912,
  8400477.59843692,
  8402713.808092939,
  8403599.736644963,
  8404043.75,
  8405123.196955431,
  8406294.64174241,
  8406318.75,
  8406332.8125,
  8410882.629978953,
  8413853.098848302,
  8414165.443707367,
  8414224.22646879,
  8416485.9375,
  8417261.647422258,
  8451413.919805456,
  8457365.400593042,
  8475653.91058664,
  8478498.4375,
  8482995.555825783,
  8485418.75,
  8505571.881481376,
  8513292.524653705,
  8513342.124890037,
  8516474.98144081,
  8520664.052955315,
  8556786.386712432,
  8556881.25,
  8566668.574877169,
  8571082.239563277,
  8574679.6875,
  8584252.705169203,
  8594757.8125,
  8598044.548606426,
  8601592.1875,
  8615725.749367699,
  8617676.989804832,
  8617726.395943604,
  8623713.128053546,
  8625082.919685105,
  8625810.646208726,
  8626141.722661693,
  8626164.0625,
  8626190.625,
  8626262.240024853,
  8627054.216234744,
  8628262.370622981,
  8638746.937816445,
  8645048.4375,
  8656671.875,
  8657466.76109016,
  8666741.42635224,
  8668084.375,
  8674427.819685213,
  8688781.629497789,
  8689140.83092631,
  8689870.3125,
  8690026.5625,
  8690426.978738997,
  8690464.0625,
  8690515.305800967,
  8690779.494858652,
  8691157.910657002,
  8695473.283444319,
  8697453.125,
  8699466.692589542,
  8699593.394955164,
  8700452.849068219,
  8706508.862549247,
  8707676.5625,
  8707938.553769981,
  8708731.245737484,
  8709403.079703972,
  8712004.6875,
  8713247.66956024,
  8714228.125,
  8716215.202061305,
  8721128.125,
  8725490.625,
  8728900.0,
  8731607.02811587,
  8738446.875,
  8768687.891110368,
  8781291.571391275,
  8783153.966033926,
  8796893.498754788,
  8797083.088780563,
  8800035.590136496,
  8800468.54824382,
  8800645.3125,
  8800786.639872646,
  8803204.53056855,
  8806076.585097393,
  8806084.272478458,
  8806703.66117119,
  8807775.0,
  8808114.0625,
  8808526.279198207,
  8809034.375,
  8826511.516911607,
  8828482.717357265,
  8828489.595247788,
  8841285.9375,
  8855067.1875,
  8859428.125,
  8864448.841614436,
  8867828.028597258,
  8868305.222270196,
  8868792.1875,
  8904137.5,
  8927442.397498403,
  8939105.686883226,
  8957125.290068386,
  8959234.375,
  8974734.053802043,
  8986701.533118851,
  8987067.1875,
  8996024.995150551,
  8996598.4375,
  9003556.240705365,
  9036503.125,
  9037611.849858824,
  9041590.625,
  9042186.222912796,
  9046012.176946638,
  9077160.120997218,
  9077220.47831315,
  9100035.9375,
  9111443.176914165,
  9158772.432460576,
  9178772.941857023,
  9219050.940352,
  9227913.652766291,
  9235752.275383314,
  9239427.7062371,
  9263566.711817905,
  9270237.33245167,
  9274719.662006212,
  9327288.211935423,
  9328164.0625,
  9399050.868746758,
  9399291.926408896,
  9601913.09858774,
  9602331.25,
  9644308.34937074,
  9649944.110416917,
  9757497.128231078,
  9806924.910490336,
  9891566.249271,
  9899257.730456265,
  9899263.942708107,
  9951846.424792135,
  9955604.389635475,
  10079533.060011435,
  10081965.172663918,
  10095268.090938635,
  10096637.727239292,
  10136351.82051684,
  10136363.62075997,
  10148972.433423938,
  10149559.282204537,
  10152686.589196257,
  10208654.366364075,
  10212271.875,
  10277493.75,
  10366323.96732527,
  10366986.022250822,
  10389824.1905386,
  10449226.54390808,
  10481341.530068649,
  10647001.5625,
  10655455.784229299,
  10712977.082335524,
  10802095.210919635,
  10890344.431180144,
  10906858.300711997,
  10927722.133428497,
  10929320.72471968,
  10967161.364347601,
  10972091.964679774,
  10974415.290450843,
  11032509.308349539,
  11034975.615031825,
  11099318.75,
  11145931.23604667,
  11145971.925563348,
  11152262.48499993,
  11167701.5625,
  11200131.755921412,
  11263233.711107356,
  11300817.373708166,
  11351764.62703371,
  11386057.086578716,
  11403738.595086506,
  11450519.71053722,
  11559345.3125,
  11571109.375,
  11683712.169916466,
  11689170.052798927,
  11718090.625,
  11796484.375,
  11798926.250507787,
  11799317.33510752,
  11855752.820890762,
  11855753.125,
  11861334.375,
  11904612.5,
  11905260.803323846,
  11940785.999630323,
  11957873.630855177,
  11968992.015292574,
  11986109.607981514,
  11988990.556150988,
  12000330.865965478,
  12000721.419723637,
  12003637.5,
  12013730.526579404,
  12025076.923100617,
  12031860.9375,
  12044717.114371343,
  12044832.262686346,
  12044866.258667452,
  12051963.891360208,
  12053504.222162308,
  12065428.700254148,
  12100260.457968684,
  12100428.875834126,
  12100521.613606047,
  12101031.944963839,
  12101641.736556266,
  12102864.0625,
  12103403.125,
  12106042.1875,
  12106053.524926012,
  12107975.885402422,
  12111437.5,
  12112353.532229342,
  12113627.02966085,
  12113981.661584968,
  12127799.623681167,
  12129737.890908703,
  12130976.18025973,
  12131888.912029162,
  12132723.133956652,
  12136144.492694957,
  12137367.321974728,
  12156985.9375,
  12178909.156252343,
  12180044.999346372,
  12182261.808479087,
  12206553.75725092,
  12213708.988271683,
  12231860.9375,
  12233003.931692105,
  12233992.1875,
  12235307.877462111,
  12236860.9375,
  12243556.25,
  12246220.337330703,
  12246220.98635424,
  12250039.0625,
  12251293.35837731,
  12270825.417941313,
  12315375.06569629,
  12316543.962705411,
  12323460.32709706,
  12324313.384600127,
  12334077.192554666,
  12336724.173568392,
  12337765.625,
  12341471.169248201,
  12353645.882447833,
  12360157.411962524,
  12360671.875,
  12361101.481186591,
  12365635.9375,
  12365693.50405218,
  12368671.148806043,
  12369346.875,
  12369726.052549252,
  12373270.876579523,
  12385782.544451762,
  12388402.343412558,
  12388435.263011193,
  12388565.625,
  12392239.753562465,
  12392659.682422675,
  12407177.137471808,
  12408483.046715355,
  12410066.679712268,
  12412177.893814208,
  12413194.53484666,
  12415309.487116607,
  12454878.125,
  12454878.125,
  12454884.113512365,
  12464404.17052088,
  12475781.25,
  12482887.5,
  12491459.017914796,
  12496264.890647713,
  12496564.0625,
  12497086.08039192,
  12510701.5625,
  12518944.5534123,
  12558012.231662245,
  12602168.205628183,
  12617085.9375,
  12617875.282408794,
  12622629.427306714,
  12625505.196172269,
  12652141.976454712,
  12670826.5625,
  12675375.64480423,
  12724159.375,
  12729757.8125,
  12731730.560226863,
  12761296.875,
  12762276.5625,
  12766442.1875,
  12767004.6875,
  12767637.776918728,
  12768873.4375,
  12773937.5,
  12775829.787915375,
  12776082.8125,
  12799215.625,
  12801985.67810539,
  12801996.664054416,
  12809708.845568748,
  12809841.287454512,
  12810412.108728705,
  12812492.666752644,
  12814362.4708218,
  12841928.84920282,
  12842684.698368622,
  12842960.9375,
  12842982.630621236,
  12843374.478768272,
  12843565.961579485,
  12851479.622186773,
  12854739.106796453,
  12860861.165745437,
  12861768.575807039,
  12879271.875,
  12883087.5,
  12896689.752215723,
  12898294.711687803,
  12910569.939069323,
  12926728.088782512,
  12944428.125,
  12945717.1875,
  12975956.235638129,
  12976982.8125,
  12977962.23576663,
  12978343.274225801,
  12978397.880479556,
  12979221.115958042,
  12979722.942118015,
  12988841.654664552,
  12996532.646194644,
  13003348.417378366,
  13005721.85358186,
  13020053.125,
  13026039.548362201,
  13047723.868260026,
  13049032.8125,
  13061197.773720844,
  13092529.575770812,
  13101117.1875,
  13139054.500913948,
  13141177.844707327,
  13142220.703160131,
  13144174.71781907,
  13149697.356564377,
  13160144.310996471,
  13161296.277499402,
  13178764.54825073,
  13180466.886455594,
  13184084.375,
  13184859.357688708,
  13192573.4375,
  13217912.5,
  13218251.5625,
  13228065.24626964,
  13250096.424802622,
  13250944.436122917,
  13251824.859801963,
  13252399.885734472,
  13254709.375,
  13258878.149024663,
  13264669.392847164,
  13291125.531940714,
  13318093.55893475,
  13436100.675369803,
  13479318.05163303,
  13512222.468342818,
  13514127.296096109,
  13549889.0625,
  13601174.125953637,
  13609099.9257834,
  13615434.23017301,
  13616828.125,
  13618156.113376312,
  13649542.772527138,
  13690618.791648932,
  13693544.741738344,
  13756411.159074683,
  13757289.912697056,
  13757604.33187155,
  13758779.005601779,
  13808909.375,
  13860483.071516888,
  13905126.390576802,
  13909102.015869262,
  13930563.39813776,
  13984470.167927423,
  13985185.9375,
  13986047.26349431,
  13994370.3125,
  14034506.25,
  14038898.12181526,
  14040118.75,
  14084625.0,
  14172919.38951191,
  14180237.5,
  14214701.28976263,
  14214884.375,
  14272776.78720475,
  14282157.629674047,
  14285638.420834383,
  14294968.540596422,
  14295910.437899055,
  14296434.53983261,
  14298643.038561702,
  14298790.096308548,
  14299470.3125,
  14301478.144789785,
  14302721.343413487,
  14304503.691940747,
  14304682.8125,
  14306492.1875,
  14308017.529457469,
  14309509.524124196,
  14351395.3125,
  14437879.802880725,
  14438519.805870693,
  14438562.677019713,
  14567965.807816833,
  14568380.335978758,
  14568715.625,
  14569565.060232587,
  14570002.829817912,
  14570965.139149569,
  14571254.6875,
  14571320.3125,
  14572818.508692514,
  14572823.4375,
  14572846.468839142,
  14572974.911198325,
  14573925.0,
  14575331.689392272,
  14576285.69654126,
  14576450.0441637,
  14576607.233619459,
  14577712.5,
  14578392.567612343,
  14578404.021282103,
  14578774.770167405,
  14579862.5,
  14579896.490098722,
  14579900.044450147,
  14580976.5625,
  14581135.9375,
  14581193.481323069,
  14581856.477213265,
  14582081.25,
  14582388.566443786,
  14582567.593741134,
  14582857.726646248,
  14583358.131934201,
  14583628.611107716,
  14583633.412135994,
  14585243.69024646,
  14586834.396304809,
  14951506.504111039,
  17830421.857919473,
  18750417.165826026,
  18751571.54319177,
  18760962.66690401,
  18766547.148982238,
  18779312.5,
  18781648.4375,
  18787685.9375,
  18791538.252323695,
  18791776.335534506,
  18791812.03376177,
  18791909.65048032,
  18793443.77453665,
  18797495.3125,
  18810867.508635696,
  18815743.75,
  102069009.375,
  102069040.625],
 [51.019727640144815,
  24.332549071475135,
  17.3993505790785,
  23.14311441386546,
  19.429381360668685,
  77.3385001613488,
  8.152664579600076,
  58.37078312040285,
  30.14550252492934,
  21.19895951299102,
  13.70664402501979,
  6.3210152915565105,
  94.72844267427124,
  37.2287108097906,
  13.025456301338632,
  108.1242463352786,
  50.00846480260808,
  66.96792301048697,
  54.69352784626684,
  85.43212882556722,
  33.071288895178554,
  23.495822414529655,
  13.551066366682026,
  64.49609518319588,
  7.781090731727417,
  5.486317898114985,
  46.04963558183315,
  38.32516084582863,
  52.748205806823364,
  72.58054141015513,
  17.16826309144363,
  33.00985311964993,
  34.37930580246373,
  109.24430228640362,
  16.253474263773896,
  50.68560737194869,
  5.753779302985908,
  27.654414848409456,
  103.96772557390341,
  50.2938702629432,
  5.8071917286905945,
  39.4497157115031,
  9.500482202494188,
  41.911703494531245,
  49.5956397582072,
  28.07611590737734,
  33.94616425981533,
  16.604391903071047,
  17.187940123130844,
  87.39315844690233,
  37.78370264227462,
  45.275673183013,
  8.699063203928658,
  8.868758347081377,
  60.39102764820626,
  23.577407407889883,
  12.691150566274143,
  90.51423539736149,
  38.11950288442491,
  76.11552355717849,
  116.32668504918468,
  7.993531997877487,
  15.201646652661815,
  92.24377710007931,
  19.27012507283209,
  17.72135863534364,
  10.893999760769823,
  8.65281675252923,
  98.1665540285828,
  23.285563504799256,
  37.07183535066272,
  8.557757016071566,
  8.663769570181527,
  23.617375709220394,
  7.0623931036725125,
  15.37601623390135,
  9.46068948486128,
  23.84301158856338,
  37.5487554006468,
  30.453815985959146,
  60.65921914809008,
  33.74751468190908,
  65.14038141773567,
  7.501680854975588,
  10.868099725358467,
  11.880200447835133,
  91.91207700294838,
  13.83988038088738,
  63.47924516550697,
  64.5213682036561,
  37.652879218916944,
  15.01707217921971,
  159.9492055892424,
  126.51231125879093,
  128.5432534119197,
  43.00794538026564,
  17.796648601868746,
  5.0778876430670365,
  24.92609179782862,
  10.238694148450953,
  86.9861887201475,
  113.97306082189817,
  5.1576346291734545,
  34.532664230909205,
  14.740611947767835,
  12.656350278708338,
  98.36547058097557,
  38.21521333248496,
  35.74932402482466,
  84.18974668802073,
  18.05929768248392,
  26.247122385519162,
  34.26502670903773,
  24.12420041879301,
  8.452737341847591,
  27.10366468398521,
  13.388659198874077,
  30.530054930213357,
  6.6022294621061794,
  24.209458760469037,
  105.36106381599684,
  7.073435009116808,
  7.214921408001459,
  52.90917582199029,
  121.58097718475993,
  81.66669052107775,
  46.29084546423743,
  17.777123146793024,
  91.8272917407803,
  29.63514990032996,
  60.5866002291116,
  146.131701741978,
  122.1088564795227,
  39.62207490400502,
  7.543701344727414,
  82.07562207406878,
  23.119231456509876,
  10.1587533115775,
  64.49473409828599,
  62.6931274471124,
  84.76796646393606,
  22.168305240012284,
  41.639322451673316,
  27.08970949779508,
  150.19195903398287,
  10.199413647942263,
  35.925234655776514,
  32.36118581894197,
  25.960615229470005,
  77.95496600398273,
  44.164314970193146,
  35.706978293407325,
  123.56040923114173,
  16.198078217869536,
  98.96520163382988,
  14.486828091488016,
  60.545482158006024,
  15.201974732971163,
  64.94936796020671,
  17.52255030025914,
  75.36788066337975,
  50.426726015037644,
  12.058613132253763,
  51.58431372510711,
  29.20118300661746,
  53.823589298265176,
  88.49258621636241,
  10.808726491800996,
  40.11516331888909,
  7.084455188617361,
  71.97939260395822,
  24.33171868404032,
  83.38245595902278,
  12.20496382123465,
  112.14656948161337,
  32.290526419026996,
  41.54608261342369,
  39.74803956805752,
  8.546518946424731,
  24.564027504530284,
  13.428970051555387,
  92.29854904916911,
  14.544394094778951,
  77.53345364600796,
  10.828143405788646,
  29.89399343235969,
  14.191961316929504,
  55.866106359313335,
  79.02966305943522,
  17.613411792059154,
  44.527727939597426,
  49.546152580837735,
  15.349617343751005,
  5.611583711422576,
  72.04395976547586,
  41.558567489368286,
  16.79906812359159,
  36.33295070169201,
  24.342699166062744,
  11.51339443961825,
  89.5195259055385,
  6.280240627168779,
  36.81996218688203,
  73.1749559687757,
  126.99589923332962,
  38.517865131338944,
  9.729979054268567,
  36.43206622709621,
  15.565585541765214,
  57.21913455550318,
  20.291031351152636,
  15.573674641968893,
  61.963341078082934,
  35.417597386158754,
  20.452516666547517,
  43.00292569181489,
  38.50699852823647,
  6.795157854379685,
  19.88445646071125,
  9.314325319610372,
  48.31806337084415,
  58.174144359128334,
  6.7780071671085995,
  37.588320043095436,
  105.82759213376322,
  58.91739354815732,
  11.90755762656985,
  9.331236539073581,
  88.75385233394404,
  60.935934874920704,
  51.55023971996729,
  28.595339460497833,
  64.62751397591268,
  46.04628958123016,
  12.895336789285572,
  6.93551065607059,
  8.448687432769841,
  96.57061404180752,
  47.4306482397842,
  16.56090360167681,
  157.3154259000464,
  18.358585003645082,
  24.700330971229143,
  58.38571643423502,
  14.855766069187162,
  48.845474273326026,
  73.2327603277362,
  115.47797534579206,
  73.81302802249368,
  31.502367548705436,
  13.039870698759033,
  16.686523897456144,
  5.576011305599249,
  12.837201418512956,
  81.36422136917308,
  52.58972869870623,
  89.50069064146214,
  28.343329062268463,
  16.64689156734508,
  74.77187712224105,
  28.187122455315027,
  5.821160852366742,
  103.52608515239372,
  14.956382739687166,
  58.46735927881736,
  10.478071806390014,
  34.12020169574375,
  67.05192816380818,
  19.747087915256788,
  14.077005818968322,
  52.48166922318646,
  67.5353250314676,
  31.507465028817442,
  6.473679617656323,
  57.58220980417575,
  89.70819536246653,
  169.9679634212193,
  25.568063759888183,
  31.168903674726693,
  43.60001106766382,
  9.98580541850406,
  27.020826201050195,
  49.31325780256245,
  80.1229596977723,
  289.73522540853793,
  19.55137396986149,
  82.05906374272587,
  8.760847838521935,
  8.933838923800112,
  19.54144890443359,
  6.473009789596844,
  23.17812158002227,
  40.16695747824426,
  62.45388966852825,
  69.32183867517249,
  9.411209049278996,
  196.72001821701687,
  75.2678657439432,
  17.88176923711083,
  26.722283778334162,
  8.58407423040281,
  17.39576355138891,
  39.08570705675115,
  5.97374468412474,
  17.298529727484897,
  19.02610162508597,
  68.96127839366913,
  6.83071188075995,
  36.652314223790974,
  73.75273602875131,
  27.96032714660041,
  63.92405419090463,
  12.75248054857289,
  35.75672133830594,
  77.52945493155681,
  25.083002825376514,
  41.85084021719274,
  21.644637769226446,
  39.86130097415324,
  41.34075203161305,
  6.795089692056224,
  41.61269582831005,
  28.765265057211426,
  13.11764922052998,
  19.74670148021197,
  43.419051781772204,
  74.61037120517116,
  126.69806054880688,
  45.60304118376206,
  15.036963766866599,
  27.611222225785276,
  101.21917185698308,
  170.53958271885588,
  5.69637279703897,
  16.43692663407476,
  128.7327693257842,
  30.550187030731397,
  37.06429330512138,
  7.059623543790945,
  19.94354879647016,
  77.25456995278871,
  59.02931626383904,
  25.377080217568487,
  50.08824868493095,
  62.04801495197008,
  28.653259257672243,
  61.94219606766775,
  18.574306011650485,
  32.19927203024752,
  89.04931689868317,
  17.973802436227004,
  34.557040364572316,
  44.11583721527938,
  51.89619477419498,
  6.26801758963527,
  38.00504178398391,
  6.751071017579627,
  20.425796054471213,
  53.13186437973502,
  16.983349901710593,
  81.14156649513676,
  9.187590237442256,
  15.74676066448239,
  46.161293410853844,
  46.88810232052171,
  43.339594208611864,
  19.210809412392116,
  57.350078578848866,
  49.055033447814246,
  9.720886433989314,
  34.592547228103086,
  104.39216011951194,
  9.423709051788068,
  8.138319032084125,
  79.17389568270941,
  46.43790197491268,
  46.22913153699679,
  93.01034920214863,
  13.40543875425305,
  138.54015206748417,
  8.320077423236954,
  36.22645147310939,
  70.50010515108235,
  19.530338328696665,
  10.115334088168769,
  62.181391195169525,
  45.52540581362601,
  49.76706700360043,
  16.504915815964594,
  95.00847027972847,
  43.336059585196594,
  6.55438182335269,
  81.50710199820907,
  23.50855982177313,
  141.13981153732965,
  53.10011727803297,
  10.89515347225526,
  62.70856628020768,
  21.544021814837333,
  13.470764836514936,
  13.984676876953692,
  49.422938205509695,
  78.7000462090955,
  10.415875445354137,
  102.2857716787643,
  27.65922151046059,
  42.3716698825728,
  8.75676749649084,
  85.77456054525874,
  22.29788730141916,
  32.29827523355036,
  85.9671628412062,
  41.56562198283246,
  72.42087838360187,
  30.61119760110114,
  9.685556521975242,
  14.755553475462332,
  26.679464583617175,
  29.65451511680259,
  17.718140344932806,
  13.764468904179129,
  21.997521845675937,
  60.50736150319451,
  12.437255890573878,
  56.765748864475405,
  38.30464209403545,
  43.520394003074514,
  63.68306246500674,
  64.97735108662387,
  40.04895538926286,
  23.966407554350187,
  48.23722491285364,
  69.93308792635801,
  26.775155563154645,
  6.336938358068589,
  48.59318691944191,
  174.550736392183,
  32.16997885619041,
  9.698196482138107,
  14.586371975303496,
  17.464770863258533,
  113.51018827508985,
  45.95399153613879,
  5.406563215586613,
  8.00966667388881,
  11.613536822909879,
  98.60064988813923,
  35.3923526281148,
  102.09263664234726,
  47.052096592033905,
  34.06027987811051,
  71.91287033111117,
  21.989748455857416,
  68.63860813273082,
  20.115613083047016,
  32.53027553404104,
  59.52971022420215,
  6.700680913799637,
  18.826533941191574,
  5.5411824603603,
  58.29512146299198,
  12.07121727270243,
  26.36180633733093,
  21.536674253680363,
  42.78802666586461,
  22.443431894950514,
  14.166873197157742,
  38.449109994219825,
  10.45755413273848,
  35.7064020383829,
  69.01427769538401,
  31.991424973693558,
  63.19702396379522,
  5.149513407433575,
  23.95635927513694,
  17.848425332107094,
  78.3588695956492,
  85.16665870520886,
  14.92361440479509,
  5.408243866311126,
  10.990621561758923,
  5.2196236460076015,
  6.474306608463468,
  45.62299201348876,
  5.589613802676459,
  61.91126842023762,
  135.98364894085228,
  21.9821258498567,
  20.04070188231627,
  8.410515552246348,
  78.65153901800285,
  118.7516431882548,
  17.604402922112953,
  8.238107655226662,
  9.950777944615773,
  9.073864876870976,
  104.75841511247097,
  43.28001990608811,
  5.741328130287127,
  12.326915406338511,
  28.82744169361292,
  11.098770257185098,
  21.726454783274562,
  79.44354172075879,
  6.040966413788864,
  17.378112574962653,
  13.364271904134796,
  65.03925902352591,
  88.69967013226335,
  5.62850281495974,
  8.517298132653174,
  33.88675751084155,
  10.796414269669068,
  30.46945185098557,
  5.067686295655706,
  12.24755992251158,
  42.77836747743613,
  55.66029067334087,
  11.20767114841728,
  158.96983506909697,
  35.772958690158,
  9.79925628108159,
  53.194216105796144,
  17.066650233275546,
  336.2615048333812,
  78.73451572037436,
  90.67608730380516,
  58.13111133272538,
  72.515178390632,
  90.4349291476602,
  5.235999902591054,
  19.277354771767495,
  98.264141736135,
  55.038779150926324,
  5.092364745993398,
  15.53563922832179,
  49.238163393636356,
  60.88776225140054,
  36.36509806061902,
  79.92001444914862,
  152.65141385148377,
  43.59354606300749,
  41.02063172511431,
  14.323391972178356,
  6.297089594210837,
  79.39886730544517,
  22.448726883922575,
  16.815763898354145,
  189.36239127126672,
  28.3736407235359,
  16.52060297249311,
  30.501089602136993,
  27.157186216090633,
  13.191657653892507,
  16.681653252502226,
  25.54923416537571,
  83.00748969813924,
  6.342720348008545,
  210.96699658732354,
  27.331477538838218,
  296.0380689426962,
  64.11204694029546,
  8.061686986837433,
  32.721722216808104,
  73.0644645775102,
  71.531021667596,
  11.816188405022837,
  100.0884314945292,
  67.81183538021493,
  104.24922352092497,
  108.54316587093628,
  16.782626760796155,
  13.678042995397185,
  136.41347889139655,
  22.75416084621835,
  25.8691869168368,
  5.3869952720990035,
  24.251599208985276,
  15.480094363000273,
  76.00058134829126,
  24.847246620879215,
  9.07994409991257,
  19.591710175419557,
  34.841003179280875,
  9.024494425834693,
  65.44328544121508,
  28.200095601828465,
  42.887748280603496,
  112.74417763732102,
  101.71799031821259,
  60.895941456786225,
  26.60800007532731,
  68.88621083963525,
  32.343477556961105,
  8.089797274311929,
  8.520690344654893,
  36.65291372550179,
  9.51797545996134,
  52.841322147188706,
  107.16929798435424,
  37.079805535883835,
  12.512894379241473,
  68.62619829987139,
  60.28862685863818,
  240.58579977733874,
  89.5225707452022,
  51.1643796027628,
  40.20264865868681,
  361.0020765450235,
  5.159672925853372,
  19.49895894680683,
  43.75040927438075,
  21.394407234518805,
  12.140944186867133,
  24.167645551223593,
  70.27009267328918,
  19.912808346496266,
  46.178156138228466,
  5.5675870871275634,
  18.295071351530858,
  17.703866855820305,
  65.47418421892348,
  18.113773588390636,
  74.89261844354837,
  19.108348713502235,
  16.129959487475794,
  34.057766171151066,
  68.70058283187828,
  24.255473221831394,
  55.260548064446,
  51.18031720735816,
  111.21770153649427,
  218.32283639086896,
  26.419700422123206,
  106.41117905135975,
  5.715912782116807,
  58.20947271200813,
  177.26046924695652,
  14.008693986519727,
  19.855863704175423,
  43.95842429372984,
  5.299891110022268,
  39.94726581426675,
  5.177379979994608,
  39.26319285419514,
  29.847025205830793,
  21.14429392689243,
  33.68065862474706,
  45.38890702766213,
  18.390976220322482,
  112.66360130545196,
  40.622440047755205,
  56.935298826492335,
  124.03833796125484,
  11.59001346102021,
  110.03837508930822,
  192.89958454632801,
  15.338575825302389,
  84.86314336068006,
  62.589779733874565,
  24.283656729927518,
  5.798535246139598,
  10.390055786107476,
  24.60512657208503,
  52.35341631273806,
  5.306124199339385,
  27.9601791132417,
  8.710662849406965,
  5.727651826563023,
  57.482683705852224,
  31.181358453263883,
  39.74587185312034,
  47.311358027807074,
  8.252337973574988,
  30.315260904229632,
  185.9137621271206,
  76.28900853866611,
  38.64456171959087,
  32.84541818743825,
  12.409204573314007,
  65.68369510997394,
  10.59020727466267,
  85.67716011457915,
  57.88003468659261,
  46.17763545910871,
  24.758403228998503,
  30.343499575477843,
  30.028883968605626,
  29.31365858752018,
  240.87515441295776,
  7.349987229521624,
  8.040903474206418,
  70.495483000691,
  80.04154274069904,
  37.5077812334381,
  14.86049712352624,
  127.51809157802892,
  44.11336346831556,
  42.708811701121235,
  51.84370794664624,
  6.957251596519924,
  7.5633209075184435,
  154.92547772479304,
  65.10510588279544,
  13.970938815513385,
  16.67373622336668,
  13.965002395384998,
  9.774732875713882,
  50.40053174167759,
  20.572659858812724,
  18.5722970477008,
  50.46364338028776,
  86.24458798055095,
  28.18677420682364,
  42.38023466856238,
  137.65555416761748,
  73.24549361865107,
  105.24085406270217,
  10.490086743301685,
  18.412659211927984,
  126.67968922692351,
  190.3264668268053,
  30.780790589445022,
  36.614259470710245,
  89.06058403203807,
  7.807671189160624,
  6.24335672835445,
  56.79672720877181,
  102.91848094951155,
  17.875086605125155,
  7.892021451760044,
  61.01470381868949,
  12.406418082477062,
  6.748779702580254,
  25.958237456887023,
  69.1698950704019,
  45.01203071874414,
  42.944023825377656,
  72.36518828557253,
  8.891312813508083,
  30.96455731915713,
  7.94269144532792,
  23.065254478573618,
  20.788734966781774,
  45.73047587764724,
  215.98970605848322,
  44.93623625691771,
  35.31738814408766,
  30.31201231925092,
  98.93881085236103,
  78.75158026739878,
  23.4595628758251,
  76.77988180322536,
  70.39598984811391,
  66.70019618152155,
  113.14505199289167,
  54.70124195187026,
  52.94349720964987,
  42.33583923855222,
  24.10534463564322,
  72.64840700400005,
  6.522902179507845,
  7.5384245499817455,
  8.841998333315564,
  35.910716564492105,
  7.432011450374535,
  5.953889844981608,
  10.172147895998716,
  21.211584023493852,
  132.64282875677196,
  46.95613120107641,
  13.240362551425536,
  84.05597771212109,
  95.4263989017441,
  47.14345365812882,
  130.13108768075568,
  15.542089658349264,
  78.26597808965136,
  17.58474488272077,
  44.439362343366284,
  9.485707184580157,
  35.39383327926229,
  12.319144334878755,
  73.70620791928121,
  22.92211956799052,
  6.057931523397967,
  31.50087497424812,
  46.71624497856206,
  68.94485331671018,
  303.6892269669984,
  67.70411086818183,
  7.462657720691568,
  123.31614311389951,
  42.44801957799693,
  5.554792542432959,
  24.535392526438656,
  7.619289424757403,
  20.740731402955074,
  55.12146533469664,
  26.2799560566004,
  79.80855636979358,
  56.83126275754098,
  59.354000627724,
  22.723620925129964,
  6.494081828629084,
  64.80388256890106,
  74.0678067262305,
  24.159269801317635,
  25.231914728561563,
  109.73274350698438,
  174.26285065951365,
  11.828181356511971,
  270.53476515284547,
  42.70299646726895,
  5.557147349491229,
  9.89570990903025,
  181.50940305007032,
  15.0233399560187,
  6.299846323086897,
  28.19122365866761,
  23.714868175546975,
  33.38075789436762,
  15.618265247795502,
  116.10276102300307,
  123.09693618944496,
  61.43580605315465,
  75.58702888738935,
  18.45953911076787,
  108.60358975513593,
  10.469569536729665,
  14.61200007351098,
  7.612136344332092,
  58.25369299131724,
  11.669723714965864,
  39.93866567719314,
  64.81639125204522,
  37.3190142815898])

Energy spectrum

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

([4928112.4989503175,
  5036230.994222343,
  5073247.653318141,
  5157092.537210761,
  5191763.327334791,
  5300340.535763368,
  5317382.550383241,
  5321083.964715635,
  5355676.5625,
  5438311.969999445,
  5451291.820822522,
  6778494.16440898,
  6785040.20378968,
  6961798.4375,
  6961824.3864494115,
  6965671.875,
  6966445.3125,
  6996359.375,
  6998190.625,
  7042068.75,
  7067759.136083348,
  7085209.378132706,
  7139805.228617017,
  7160751.5625,
  7173569.216348684,
  7174431.809798662,
  7183017.109251836,
  7267945.3125,
  7270080.18207427,
  7277856.25,
  7297820.795334424,
  7308457.590299777,
  7325564.0625,
  7338440.044851796,
  7339704.23392173,
  7383595.006030557,
  7408210.089322634,
  7439252.530490686,
  7449281.994137895,
  7450372.632883947,
  7450440.53282894,
  7453479.6875,
  7460293.943414853,
  7461414.0625,
  7462551.5625,
  7464193.726921775,
  7464526.5625,
  7465321.813099671,
  7467959.614432444,
  7501741.049957041,
  7510363.752011988,
  7514546.034548267,
  7520335.081267191,
  7520371.272001179,
  7522658.281426924,
  7524769.998808046,
  7528399.2505654795,
  7529044.037948253,
  7541440.625,
  7546482.8125,
  7546489.0625,
  7547263.463943135,
  7553357.8791160025,
  7561035.9375,
  7577926.50891812,
  7589500.577856496,
  7589963.907652801,
  7591550.633801611,
  7592200.0,
  7592793.412688895,
  7594035.9375,
  7595455.167362315,
  7597749.088604033,
  7597774.573422763,
  7598457.582963947,
  7600586.047011132,
  7600948.586543177,
  7601541.609075746,
  7607356.241776778,
  7608885.987471846,
  7617818.75,
  7623774.413744111,
  7623906.25,
  7628007.522701366,
  7631127.952310165,
  7642312.933400372,
  7642623.4375,
  7646555.927471407,
  7659704.424472292,
  7660423.246718515,
  7661290.781656286,
  7667187.519684521,
  7669906.820943252,
  7670069.640706044,
  7676299.3879393805,
  7677092.784260475,
  7693406.840187491,
  7693458.300954373,
  7695421.661761349,
  7696093.095569686,
  7698765.346409517,
  7713493.75,
  7721913.727597748,
  7722556.980149264,
  7724259.525576042,
  7724515.863845837,
  7724843.75,
  7727261.023425729,
  7727271.356537071,
  7727415.625,
  7727877.194644847,
  7728186.571872057,
  7728515.680084369,
  7729054.871836707,
  7729200.301527607,
  7730979.611331726,
  7731413.585203696,
  7733510.9375,
  7734096.009670198,
  7740732.768387361,
  7741889.0625,
  7741920.126632115,
  7744231.594028677,
  7747357.8125,
  7747757.514116983,
  7749698.4375,
  7751104.6875,
  7752864.455034001,
  7756960.9375,
  7764938.664784911,
  7766651.5625,
  7767120.3125,
  7767438.239700456,
  7767926.5625,
  7768345.144835219,
  7768367.540090685,
  7768904.336399147,
  7768928.964429455,
  7768935.9375,
  7769650.0,
  7770121.04240849,
  7773571.440123174,
  7773998.411369783,
  7774192.110536828,
  7774225.0,
  7774447.749575863,
  7774453.125,
  7775539.0625,
  7777068.605323072,
  7777170.3125,
  7782460.9375,
  7782693.75,
  7783134.910419227,
  7783312.794751378,
  7783749.427500115,
  7784467.062865479,
  7787787.064645441,
  7795668.816698004,
  7795856.25,
  7796655.232118821,
  7797149.858760677,
  7802854.6875,
  7807122.39877029,
  7807973.4375,
  7826756.25,
  7826756.25,
  7842487.447523741,
  7844718.535868823,
  7851614.102293992,
  7865723.294864476,
  7870084.375,
  7874134.848562917,
  7877729.353606684,
  7878374.032340517,
  7890621.875,
  7902452.33385658,
  7921685.9375,
  7921687.5,
  7924925.46120395,
  7938648.198609276,
  7959710.47465506,
  7985421.918472736,
  7992271.766758991,
  7995504.696577904,
  7997695.11200486,
  8000735.9375,
  8023782.600891208,
  8026537.5,
  8026938.284448174,
  8027025.55993296,
  8027387.111576463,
  8027482.8125,
  8028878.228104227,
  8031037.455439818,
  8031115.537842666,
  8033304.051546997,
  8039536.7103334535,
  8043325.73759618,
  8044762.24497101,
  8045268.983490752,
  8045851.5625,
  8048845.415036554,
  8049198.030699028,
  8049636.389559176,
  8050270.3125,
  8050626.5625,
  8050758.249043236,
  8050821.624179483,
  8051181.404633995,
  8051940.625,
  8052812.131277288,
  8053056.706105952,
  8054601.5625,
  8057075.833765266,
  8063735.326451124,
  8066506.25,
  8076059.375,
  8076200.156499177,
  8078877.447056814,
  8079407.79223353,
  8079756.25,
  8081642.1875,
  8085040.697201559,
  8089092.171798846,
  8089559.375,
  8089950.870734111,
  8090286.385033847,
  8090292.187302922,
  8092241.566301165,
  8094739.0625,
  8104568.75,
  8122196.365773838,
  8124370.317114247,
  8145745.3125,
  8153418.084821994,
  8155963.30454557,
  8181686.231968025,
  8186329.011787595,
  8196593.75,
  8199386.539502798,
  8200132.462763735,
  8200161.615539064,
  8201437.270873369,
  8201957.7435368,
  8204474.977320864,
  8207768.75,
  8208186.36743079,
  8208777.336027159,
  8213970.3125,
  8218710.171837039,
  8220077.517073625,
  8220664.455636583,
  8220664.916431453,
  8233194.623744854,
  8240549.360318952,
  8242492.1875,
  8251937.5,
  8270687.059177668,
  8279017.572944678,
  8287757.8125,
  8294750.347112862,
  8297646.971576843,
  8300795.3125,
  8302445.316396724,
  8305570.725915888,
  8309044.852922054,
  8311505.633486733,
  8312116.216577061,
  8312432.117639927,
  8313599.753057187,
  8313651.069343612,
  8314070.3125,
  8314099.423960631,
  8315051.476407381,
  8315123.4375,
  8316479.6875,
  8316499.716329852,
  8317019.467708244,
  8317121.875,
  8317151.5625,
  8317175.684282461,
  8318581.16909885,
  8319086.0302196145,
  8319989.0625,
  8319990.625,
  8320053.968010188,
  8320442.108191726,
  8321211.738216368,
  8322905.739160612,
  8324197.45631981,
  8326223.3470792435,
  8327634.02880927,
  8340076.558805409,
  8344870.6528990865,
  8345062.327534265,
  8345539.155398426,
  8345935.26695746,
  8346270.3125,
  8348530.291463457,
  8350254.5892514,
  8351448.960669997,
  8352077.085393564,
  8353981.25,
  8354491.011043467,
  8354655.271274155,
  8354942.961162562,
  8357061.836866527,
  8359675.327558865,
  8361828.597066413,
  8362179.213305795,
  8362197.206458042,
  8362717.1875,
  8364066.463802592,
  8365800.331625924,
  8366057.262007171,
  8366451.360100323,
  8371052.421910109,
  8371646.3900579205,
  8372870.3125,
  8373288.062428359,
  8373307.968666479,
  8376871.875,
  8377524.114762574,
  8377525.984602657,
  8385687.1533857025,
  8392364.0625,
  8393432.418561727,
  8393451.884448973,
  8393855.74306448,
  8395354.010603042,
  8395506.202094821,
  8398071.654072808,
  8399175.584507551,
  8399606.085573912,
  8400477.59843692,
  8402713.808092939,
  8403599.736644963,
  8404043.75,
  8405123.196955431,
  8406294.64174241,
  8406318.75,
  8406332.8125,
  8410882.629978953,
  8413853.098848302,
  8414165.443707367,
  8414224.22646879,
  8416485.9375,
  8417261.647422258,
  8451413.919805456,
  8457365.400593042,
  8475653.91058664,
  8478498.4375,
  8482995.555825783,
  8485418.75,
  8505571.881481376,
  8513292.524653705,
  8513342.124890037,
  8516474.98144081,
  8520664.052955315,
  8556786.386712432,
  8556881.25,
  8566668.574877169,
  8571082.239563277,
  8574679.6875,
  8584252.705169203,
  8594757.8125,
  8598044.548606426,
  8601592.1875,
  8615725.749367699,
  8617676.989804832,
  8617726.395943604,
  8623713.128053546,
  8625082.919685105,
  8625810.646208726,
  8626141.722661693,
  8626164.0625,
  8626190.625,
  8626262.240024853,
  8627054.216234744,
  8628262.370622981,
  8638746.937816445,
  8645048.4375,
  8656671.875,
  8657466.76109016,
  8666741.42635224,
  8668084.375,
  8674427.819685213,
  8688781.629497789,
  8689140.83092631,
  8689870.3125,
  8690026.5625,
  8690426.978738997,
  8690464.0625,
  8690515.305800967,
  8690779.494858652,
  8691157.910657002,
  8695473.283444319,
  8697453.125,
  8699466.692589542,
  8699593.394955164,
  8700452.849068219,
  8706508.862549247,
  8707676.5625,
  8707938.553769981,
  8708731.245737484,
  8709403.079703972,
  8712004.6875,
  8713247.66956024,
  8714228.125,
  8716215.202061305,
  8721128.125,
  8725490.625,
  8728900.0,
  8731607.02811587,
  8738446.875,
  8768687.891110368,
  8781291.571391275,
  8783153.966033926,
  8796893.498754788,
  8797083.088780563,
  8800035.590136496,
  8800468.54824382,
  8800645.3125,
  8800786.639872646,
  8803204.53056855,
  8806076.585097393,
  8806084.272478458,
  8806703.66117119,
  8807775.0,
  8808114.0625,
  8808526.279198207,
  8809034.375,
  8826511.516911607,
  8828482.717357265,
  8828489.595247788,
  8841285.9375,
  8855067.1875,
  8859428.125,
  8864448.841614436,
  8867828.028597258,
  8868305.222270196,
  8868792.1875,
  8904137.5,
  8927442.397498403,
  8939105.686883226,
  8957125.290068386,
  8959234.375,
  8974734.053802043,
  8986701.533118851,
  8987067.1875,
  8996024.995150551,
  8996598.4375,
  9003556.240705365,
  9036503.125,
  9037611.849858824,
  9041590.625,
  9042186.222912796,
  9046012.176946638,
  9077160.120997218,
  9077220.47831315,
  9100035.9375,
  9111443.176914165,
  9158772.432460576,
  9178772.941857023,
  9219050.940352,
  9227913.652766291,
  9235752.275383314,
  9239427.7062371,
  9263566.711817905,
  9270237.33245167,
  9274719.662006212,
  9327288.211935423,
  9328164.0625,
  9399050.868746758,
  9399291.926408896,
  9601913.09858774,
  9602331.25,
  9644308.34937074,
  9649944.110416917,
  9757497.128231078,
  9806924.910490336,
  9891566.249271,
  9899257.730456265,
  9899263.942708107,
  9951846.424792135,
  9955604.389635475,
  10079533.060011435,
  10081965.172663918,
  10095268.090938635,
  10096637.727239292,
  10136351.82051684,
  10136363.62075997,
  10148972.433423938,
  10149559.282204537,
  10152686.589196257,
  10208654.366364075,
  10212271.875,
  10277493.75,
  10366323.96732527,
  10366986.022250822,
  10389824.1905386,
  10449226.54390808,
  10481341.530068649,
  10647001.5625,
  10655455.784229299,
  10712977.082335524,
  10802095.210919635,
  10890344.431180144,
  10906858.300711997,
  10927722.133428497,
  10929320.72471968,
  10967161.364347601,
  10972091.964679774,
  10974415.290450843,
  11032509.308349539,
  11034975.615031825,
  11099318.75,
  11145931.23604667,
  11145971.925563348,
  11152262.48499993,
  11167701.5625,
  11200131.755921412,
  11263233.711107356,
  11300817.373708166,
  11351764.62703371,
  11386057.086578716,
  11403738.595086506,
  11450519.71053722,
  11559345.3125,
  11571109.375,
  11683712.169916466,
  11689170.052798927,
  11718090.625,
  11796484.375,
  11798926.250507787,
  11799317.33510752,
  11855752.820890762,
  11855753.125,
  11861334.375,
  11904612.5,
  11905260.803323846,
  11940785.999630323,
  11957873.630855177,
  11968992.015292574,
  11986109.607981514,
  11988990.556150988,
  12000330.865965478,
  12000721.419723637,
  12003637.5,
  12013730.526579404,
  12025076.923100617,
  12031860.9375,
  12044717.114371343,
  12044832.262686346,
  12044866.258667452,
  12051963.891360208,
  12053504.222162308,
  12065428.700254148,
  12100260.457968684,
  12100428.875834126,
  12100521.613606047,
  12101031.944963839,
  12101641.736556266,
  12102864.0625,
  12103403.125,
  12106042.1875,
  12106053.524926012,
  12107975.885402422,
  12111437.5,
  12112353.532229342,
  12113627.02966085,
  12113981.661584968,
  12127799.623681167,
  12129737.890908703,
  12130976.18025973,
  12131888.912029162,
  12132723.133956652,
  12136144.492694957,
  12137367.321974728,
  12156985.9375,
  12178909.156252343,
  12180044.999346372,
  12182261.808479087,
  12206553.75725092,
  12213708.988271683,
  12231860.9375,
  12233003.931692105,
  12233992.1875,
  12235307.877462111,
  12236860.9375,
  12243556.25,
  12246220.337330703,
  12246220.98635424,
  12250039.0625,
  12251293.35837731,
  12270825.417941313,
  12315375.06569629,
  12316543.962705411,
  12323460.32709706,
  12324313.384600127,
  12334077.192554666,
  12336724.173568392,
  12337765.625,
  12341471.169248201,
  12353645.882447833,
  12360157.411962524,
  12360671.875,
  12361101.481186591,
  12365635.9375,
  12365693.50405218,
  12368671.148806043,
  12369346.875,
  12369726.052549252,
  12373270.876579523,
  12385782.544451762,
  12388402.343412558,
  12388435.263011193,
  12388565.625,
  12392239.753562465,
  12392659.682422675,
  12407177.137471808,
  12408483.046715355,
  12410066.679712268,
  12412177.893814208,
  12413194.53484666,
  12415309.487116607,
  12454878.125,
  12454878.125,
  12454884.113512365,
  12464404.17052088,
  12475781.25,
  12482887.5,
  12491459.017914796,
  12496264.890647713,
  12496564.0625,
  12497086.08039192,
  12510701.5625,
  12518944.5534123,
  12558012.231662245,
  12602168.205628183,
  12617085.9375,
  12617875.282408794,
  12622629.427306714,
  12625505.196172269,
  12652141.976454712,
  12670826.5625,
  12675375.64480423,
  12724159.375,
  12729757.8125,
  12731730.560226863,
  12761296.875,
  12762276.5625,
  12766442.1875,
  12767004.6875,
  12767637.776918728,
  12768873.4375,
  12773937.5,
  12775829.787915375,
  12776082.8125,
  12799215.625,
  12801985.67810539,
  12801996.664054416,
  12809708.845568748,
  12809841.287454512,
  12810412.108728705,
  12812492.666752644,
  12814362.4708218,
  12841928.84920282,
  12842684.698368622,
  12842960.9375,
  12842982.630621236,
  12843374.478768272,
  12843565.961579485,
  12851479.622186773,
  12854739.106796453,
  12860861.165745437,
  12861768.575807039,
  12879271.875,
  12883087.5,
  12896689.752215723,
  12898294.711687803,
  12910569.939069323,
  12926728.088782512,
  12944428.125,
  12945717.1875,
  12975956.235638129,
  12976982.8125,
  12977962.23576663,
  12978343.274225801,
  12978397.880479556,
  12979221.115958042,
  12979722.942118015,
  12988841.654664552,
  12996532.646194644,
  13003348.417378366,
  13005721.85358186,
  13020053.125,
  13026039.548362201,
  13047723.868260026,
  13049032.8125,
  13061197.773720844,
  13092529.575770812,
  13101117.1875,
  13139054.500913948,
  13141177.844707327,
  13142220.703160131,
  13144174.71781907,
  13149697.356564377,
  13160144.310996471,
  13161296.277499402,
  13178764.54825073,
  13180466.886455594,
  13184084.375,
  13184859.357688708,
  13192573.4375,
  13217912.5,
  13218251.5625,
  13228065.24626964,
  13250096.424802622,
  13250944.436122917,
  13251824.859801963,
  13252399.885734472,
  13254709.375,
  13258878.149024663,
  13264669.392847164,
  13291125.531940714,
  13318093.55893475,
  13436100.675369803,
  13479318.05163303,
  13512222.468342818,
  13514127.296096109,
  13549889.0625,
  13601174.125953637,
  13609099.9257834,
  13615434.23017301,
  13616828.125,
  13618156.113376312,
  13649542.772527138,
  13690618.791648932,
  13693544.741738344,
  13756411.159074683,
  13757289.912697056,
  13757604.33187155,
  13758779.005601779,
  13808909.375,
  13860483.071516888,
  13905126.390576802,
  13909102.015869262,
  13930563.39813776,
  13984470.167927423,
  13985185.9375,
  13986047.26349431,
  13994370.3125,
  14034506.25,
  14038898.12181526,
  14040118.75,
  14084625.0,
  14172919.38951191,
  14180237.5,
  14214701.28976263,
  14214884.375,
  14272776.78720475,
  14282157.629674047,
  14285638.420834383,
  14294968.540596422,
  14295910.437899055,
  14296434.53983261,
  14298643.038561702,
  14298790.096308548,
  14299470.3125,
  14301478.144789785,
  14302721.343413487,
  14304503.691940747,
  14304682.8125,
  14306492.1875,
  14308017.529457469,
  14309509.524124196,
  14351395.3125,
  14437879.802880725,
  14438519.805870693,
  14438562.677019713,
  14567965.807816833,
  14568380.335978758,
  14568715.625,
  14569565.060232587,
  14570002.829817912,
  14570965.139149569,
  14571254.6875,
  14571320.3125,
  14572818.508692514,
  14572823.4375,
  14572846.468839142,
  14572974.911198325,
  14573925.0,
  14575331.689392272,
  14576285.69654126,
  14576450.0441637,
  14576607.233619459,
  14577712.5,
  14578392.567612343,
  14578404.021282103,
  14578774.770167405,
  14579862.5,
  14579896.490098722,
  14579900.044450147,
  14580976.5625,
  14581135.9375,
  14581193.481323069,
  14581856.477213265,
  14582081.25,
  14582388.566443786,
  14582567.593741134,
  14582857.726646248,
  14583358.131934201,
  14583628.611107716,
  14583633.412135994,
  14585243.69024646,
  14586834.396304809,
  14951506.504111039,
  17830421.857919473,
  18750417.165826026,
  18751571.54319177,
  18760962.66690401,
  18766547.148982238,
  18779312.5,
  18781648.4375,
  18787685.9375,
  18791538.252323695,
  18791776.335534506,
  18791812.03376177,
  18791909.65048032,
  18793443.77453665,
  18797495.3125,
  18810867.508635696,
  18815743.75,
  102069009.375,
  102069040.625],
 [51.019727640144815,
  24.332549071475135,
  17.3993505790785,
  23.14311441386546,
  19.429381360668685,
  77.3385001613488,
  8.152664579600076,
  58.37078312040285,
  30.14550252492934,
  21.19895951299102,
  13.70664402501979,
  6.3210152915565105,
  94.72844267427124,
  37.2287108097906,
  13.025456301338632,
  108.1242463352786,
  50.00846480260808,
  66.96792301048697,
  54.69352784626684,
  85.43212882556722,
  33.071288895178554,
  23.495822414529655,
  13.551066366682026,
  64.49609518319588,
  7.781090731727417,
  5.486317898114985,
  46.04963558183315,
  38.32516084582863,
  52.748205806823364,
  72.58054141015513,
  17.16826309144363,
  33.00985311964993,
  34.37930580246373,
  109.24430228640362,
  16.253474263773896,
  50.68560737194869,
  5.753779302985908,
  27.654414848409456,
  103.96772557390341,
  50.2938702629432,
  5.8071917286905945,
  39.4497157115031,
  9.500482202494188,
  41.911703494531245,
  49.5956397582072,
  28.07611590737734,
  33.94616425981533,
  16.604391903071047,
  17.187940123130844,
  87.39315844690233,
  37.78370264227462,
  45.275673183013,
  8.699063203928658,
  8.868758347081377,
  60.39102764820626,
  23.577407407889883,
  12.691150566274143,
  90.51423539736149,
  38.11950288442491,
  76.11552355717849,
  116.32668504918468,
  7.993531997877487,
  15.201646652661815,
  92.24377710007931,
  19.27012507283209,
  17.72135863534364,
  10.893999760769823,
  8.65281675252923,
  98.1665540285828,
  23.285563504799256,
  37.07183535066272,
  8.557757016071566,
  8.663769570181527,
  23.617375709220394,
  7.0623931036725125,
  15.37601623390135,
  9.46068948486128,
  23.84301158856338,
  37.5487554006468,
  30.453815985959146,
  60.65921914809008,
  33.74751468190908,
  65.14038141773567,
  7.501680854975588,
  10.868099725358467,
  11.880200447835133,
  91.91207700294838,
  13.83988038088738,
  63.47924516550697,
  64.5213682036561,
  37.652879218916944,
  15.01707217921971,
  159.9492055892424,
  126.51231125879093,
  128.5432534119197,
  43.00794538026564,
  17.796648601868746,
  5.0778876430670365,
  24.92609179782862,
  10.238694148450953,
  86.9861887201475,
  113.97306082189817,
  5.1576346291734545,
  34.532664230909205,
  14.740611947767835,
  12.656350278708338,
  98.36547058097557,
  38.21521333248496,
  35.74932402482466,
  84.18974668802073,
  18.05929768248392,
  26.247122385519162,
  34.26502670903773,
  24.12420041879301,
  8.452737341847591,
  27.10366468398521,
  13.388659198874077,
  30.530054930213357,
  6.6022294621061794,
  24.209458760469037,
  105.36106381599684,
  7.073435009116808,
  7.214921408001459,
  52.90917582199029,
  121.58097718475993,
  81.66669052107775,
  46.29084546423743,
  17.777123146793024,
  91.8272917407803,
  29.63514990032996,
  60.5866002291116,
  146.131701741978,
  122.1088564795227,
  39.62207490400502,
  7.543701344727414,
  82.07562207406878,
  23.119231456509876,
  10.1587533115775,
  64.49473409828599,
  62.6931274471124,
  84.76796646393606,
  22.168305240012284,
  41.639322451673316,
  27.08970949779508,
  150.19195903398287,
  10.199413647942263,
  35.925234655776514,
  32.36118581894197,
  25.960615229470005,
  77.95496600398273,
  44.164314970193146,
  35.706978293407325,
  123.56040923114173,
  16.198078217869536,
  98.96520163382988,
  14.486828091488016,
  60.545482158006024,
  15.201974732971163,
  64.94936796020671,
  17.52255030025914,
  75.36788066337975,
  50.426726015037644,
  12.058613132253763,
  51.58431372510711,
  29.20118300661746,
  53.823589298265176,
  88.49258621636241,
  10.808726491800996,
  40.11516331888909,
  7.084455188617361,
  71.97939260395822,
  24.33171868404032,
  83.38245595902278,
  12.20496382123465,
  112.14656948161337,
  32.290526419026996,
  41.54608261342369,
  39.74803956805752,
  8.546518946424731,
  24.564027504530284,
  13.428970051555387,
  92.29854904916911,
  14.544394094778951,
  77.53345364600796,
  10.828143405788646,
  29.89399343235969,
  14.191961316929504,
  55.866106359313335,
  79.02966305943522,
  17.613411792059154,
  44.527727939597426,
  49.546152580837735,
  15.349617343751005,
  5.611583711422576,
  72.04395976547586,
  41.558567489368286,
  16.79906812359159,
  36.33295070169201,
  24.342699166062744,
  11.51339443961825,
  89.5195259055385,
  6.280240627168779,
  36.81996218688203,
  73.1749559687757,
  126.99589923332962,
  38.517865131338944,
  9.729979054268567,
  36.43206622709621,
  15.565585541765214,
  57.21913455550318,
  20.291031351152636,
  15.573674641968893,
  61.963341078082934,
  35.417597386158754,
  20.452516666547517,
  43.00292569181489,
  38.50699852823647,
  6.795157854379685,
  19.88445646071125,
  9.314325319610372,
  48.31806337084415,
  58.174144359128334,
  6.7780071671085995,
  37.588320043095436,
  105.82759213376322,
  58.91739354815732,
  11.90755762656985,
  9.331236539073581,
  88.75385233394404,
  60.935934874920704,
  51.55023971996729,
  28.595339460497833,
  64.62751397591268,
  46.04628958123016,
  12.895336789285572,
  6.93551065607059,
  8.448687432769841,
  96.57061404180752,
  47.4306482397842,
  16.56090360167681,
  157.3154259000464,
  18.358585003645082,
  24.700330971229143,
  58.38571643423502,
  14.855766069187162,
  48.845474273326026,
  73.2327603277362,
  115.47797534579206,
  73.81302802249368,
  31.502367548705436,
  13.039870698759033,
  16.686523897456144,
  5.576011305599249,
  12.837201418512956,
  81.36422136917308,
  52.58972869870623,
  89.50069064146214,
  28.343329062268463,
  16.64689156734508,
  74.77187712224105,
  28.187122455315027,
  5.821160852366742,
  103.52608515239372,
  14.956382739687166,
  58.46735927881736,
  10.478071806390014,
  34.12020169574375,
  67.05192816380818,
  19.747087915256788,
  14.077005818968322,
  52.48166922318646,
  67.5353250314676,
  31.507465028817442,
  6.473679617656323,
  57.58220980417575,
  89.70819536246653,
  169.9679634212193,
  25.568063759888183,
  31.168903674726693,
  43.60001106766382,
  9.98580541850406,
  27.020826201050195,
  49.31325780256245,
  80.1229596977723,
  289.73522540853793,
  19.55137396986149,
  82.05906374272587,
  8.760847838521935,
  8.933838923800112,
  19.54144890443359,
  6.473009789596844,
  23.17812158002227,
  40.16695747824426,
  62.45388966852825,
  69.32183867517249,
  9.411209049278996,
  196.72001821701687,
  75.2678657439432,
  17.88176923711083,
  26.722283778334162,
  8.58407423040281,
  17.39576355138891,
  39.08570705675115,
  5.97374468412474,
  17.298529727484897,
  19.02610162508597,
  68.96127839366913,
  6.83071188075995,
  36.652314223790974,
  73.75273602875131,
  27.96032714660041,
  63.92405419090463,
  12.75248054857289,
  35.75672133830594,
  77.52945493155681,
  25.083002825376514,
  41.85084021719274,
  21.644637769226446,
  39.86130097415324,
  41.34075203161305,
  6.795089692056224,
  41.61269582831005,
  28.765265057211426,
  13.11764922052998,
  19.74670148021197,
  43.419051781772204,
  74.61037120517116,
  126.69806054880688,
  45.60304118376206,
  15.036963766866599,
  27.611222225785276,
  101.21917185698308,
  170.53958271885588,
  5.69637279703897,
  16.43692663407476,
  128.7327693257842,
  30.550187030731397,
  37.06429330512138,
  7.059623543790945,
  19.94354879647016,
  77.25456995278871,
  59.02931626383904,
  25.377080217568487,
  50.08824868493095,
  62.04801495197008,
  28.653259257672243,
  61.94219606766775,
  18.574306011650485,
  32.19927203024752,
  89.04931689868317,
  17.973802436227004,
  34.557040364572316,
  44.11583721527938,
  51.89619477419498,
  6.26801758963527,
  38.00504178398391,
  6.751071017579627,
  20.425796054471213,
  53.13186437973502,
  16.983349901710593,
  81.14156649513676,
  9.187590237442256,
  15.74676066448239,
  46.161293410853844,
  46.88810232052171,
  43.339594208611864,
  19.210809412392116,
  57.350078578848866,
  49.055033447814246,
  9.720886433989314,
  34.592547228103086,
  104.39216011951194,
  9.423709051788068,
  8.138319032084125,
  79.17389568270941,
  46.43790197491268,
  46.22913153699679,
  93.01034920214863,
  13.40543875425305,
  138.54015206748417,
  8.320077423236954,
  36.22645147310939,
  70.50010515108235,
  19.530338328696665,
  10.115334088168769,
  62.181391195169525,
  45.52540581362601,
  49.76706700360043,
  16.504915815964594,
  95.00847027972847,
  43.336059585196594,
  6.55438182335269,
  81.50710199820907,
  23.50855982177313,
  141.13981153732965,
  53.10011727803297,
  10.89515347225526,
  62.70856628020768,
  21.544021814837333,
  13.470764836514936,
  13.984676876953692,
  49.422938205509695,
  78.7000462090955,
  10.415875445354137,
  102.2857716787643,
  27.65922151046059,
  42.3716698825728,
  8.75676749649084,
  85.77456054525874,
  22.29788730141916,
  32.29827523355036,
  85.9671628412062,
  41.56562198283246,
  72.42087838360187,
  30.61119760110114,
  9.685556521975242,
  14.755553475462332,
  26.679464583617175,
  29.65451511680259,
  17.718140344932806,
  13.764468904179129,
  21.997521845675937,
  60.50736150319451,
  12.437255890573878,
  56.765748864475405,
  38.30464209403545,
  43.520394003074514,
  63.68306246500674,
  64.97735108662387,
  40.04895538926286,
  23.966407554350187,
  48.23722491285364,
  69.93308792635801,
  26.775155563154645,
  6.336938358068589,
  48.59318691944191,
  174.550736392183,
  32.16997885619041,
  9.698196482138107,
  14.586371975303496,
  17.464770863258533,
  113.51018827508985,
  45.95399153613879,
  5.406563215586613,
  8.00966667388881,
  11.613536822909879,
  98.60064988813923,
  35.3923526281148,
  102.09263664234726,
  47.052096592033905,
  34.06027987811051,
  71.91287033111117,
  21.989748455857416,
  68.63860813273082,
  20.115613083047016,
  32.53027553404104,
  59.52971022420215,
  6.700680913799637,
  18.826533941191574,
  5.5411824603603,
  58.29512146299198,
  12.07121727270243,
  26.36180633733093,
  21.536674253680363,
  42.78802666586461,
  22.443431894950514,
  14.166873197157742,
  38.449109994219825,
  10.45755413273848,
  35.7064020383829,
  69.01427769538401,
  31.991424973693558,
  63.19702396379522,
  5.149513407433575,
  23.95635927513694,
  17.848425332107094,
  78.3588695956492,
  85.16665870520886,
  14.92361440479509,
  5.408243866311126,
  10.990621561758923,
  5.2196236460076015,
  6.474306608463468,
  45.62299201348876,
  5.589613802676459,
  61.91126842023762,
  135.98364894085228,
  21.9821258498567,
  20.04070188231627,
  8.410515552246348,
  78.65153901800285,
  118.7516431882548,
  17.604402922112953,
  8.238107655226662,
  9.950777944615773,
  9.073864876870976,
  104.75841511247097,
  43.28001990608811,
  5.741328130287127,
  12.326915406338511,
  28.82744169361292,
  11.098770257185098,
  21.726454783274562,
  79.44354172075879,
  6.040966413788864,
  17.378112574962653,
  13.364271904134796,
  65.03925902352591,
  88.69967013226335,
  5.62850281495974,
  8.517298132653174,
  33.88675751084155,
  10.796414269669068,
  30.46945185098557,
  5.067686295655706,
  12.24755992251158,
  42.77836747743613,
  55.66029067334087,
  11.20767114841728,
  158.96983506909697,
  35.772958690158,
  9.79925628108159,
  53.194216105796144,
  17.066650233275546,
  336.2615048333812,
  78.73451572037436,
  90.67608730380516,
  58.13111133272538,
  72.515178390632,
  90.4349291476602,
  5.235999902591054,
  19.277354771767495,
  98.264141736135,
  55.038779150926324,
  5.092364745993398,
  15.53563922832179,
  49.238163393636356,
  60.88776225140054,
  36.36509806061902,
  79.92001444914862,
  152.65141385148377,
  43.59354606300749,
  41.02063172511431,
  14.323391972178356,
  6.297089594210837,
  79.39886730544517,
  22.448726883922575,
  16.815763898354145,
  189.36239127126672,
  28.3736407235359,
  16.52060297249311,
  30.501089602136993,
  27.157186216090633,
  13.191657653892507,
  16.681653252502226,
  25.54923416537571,
  83.00748969813924,
  6.342720348008545,
  210.96699658732354,
  27.331477538838218,
  296.0380689426962,
  64.11204694029546,
  8.061686986837433,
  32.721722216808104,
  73.0644645775102,
  71.531021667596,
  11.816188405022837,
  100.0884314945292,
  67.81183538021493,
  104.24922352092497,
  108.54316587093628,
  16.782626760796155,
  13.678042995397185,
  136.41347889139655,
  22.75416084621835,
  25.8691869168368,
  5.3869952720990035,
  24.251599208985276,
  15.480094363000273,
  76.00058134829126,
  24.847246620879215,
  9.07994409991257,
  19.591710175419557,
  34.841003179280875,
  9.024494425834693,
  65.44328544121508,
  28.200095601828465,
  42.887748280603496,
  112.74417763732102,
  101.71799031821259,
  60.895941456786225,
  26.60800007532731,
  68.88621083963525,
  32.343477556961105,
  8.089797274311929,
  8.520690344654893,
  36.65291372550179,
  9.51797545996134,
  52.841322147188706,
  107.16929798435424,
  37.079805535883835,
  12.512894379241473,
  68.62619829987139,
  60.28862685863818,
  240.58579977733874,
  89.5225707452022,
  51.1643796027628,
  40.20264865868681,
  361.0020765450235,
  5.159672925853372,
  19.49895894680683,
  43.75040927438075,
  21.394407234518805,
  12.140944186867133,
  24.167645551223593,
  70.27009267328918,
  19.912808346496266,
  46.178156138228466,
  5.5675870871275634,
  18.295071351530858,
  17.703866855820305,
  65.47418421892348,
  18.113773588390636,
  74.89261844354837,
  19.108348713502235,
  16.129959487475794,
  34.057766171151066,
  68.70058283187828,
  24.255473221831394,
  55.260548064446,
  51.18031720735816,
  111.21770153649427,
  218.32283639086896,
  26.419700422123206,
  106.41117905135975,
  5.715912782116807,
  58.20947271200813,
  177.26046924695652,
  14.008693986519727,
  19.855863704175423,
  43.95842429372984,
  5.299891110022268,
  39.94726581426675,
  5.177379979994608,
  39.26319285419514,
  29.847025205830793,
  21.14429392689243,
  33.68065862474706,
  45.38890702766213,
  18.390976220322482,
  112.66360130545196,
  40.622440047755205,
  56.935298826492335,
  124.03833796125484,
  11.59001346102021,
  110.03837508930822,
  192.89958454632801,
  15.338575825302389,
  84.86314336068006,
  62.589779733874565,
  24.283656729927518,
  5.798535246139598,
  10.390055786107476,
  24.60512657208503,
  52.35341631273806,
  5.306124199339385,
  27.9601791132417,
  8.710662849406965,
  5.727651826563023,
  57.482683705852224,
  31.181358453263883,
  39.74587185312034,
  47.311358027807074,
  8.252337973574988,
  30.315260904229632,
  185.9137621271206,
  76.28900853866611,
  38.64456171959087,
  32.84541818743825,
  12.409204573314007,
  65.68369510997394,
  10.59020727466267,
  85.67716011457915,
  57.88003468659261,
  46.17763545910871,
  24.758403228998503,
  30.343499575477843,
  30.028883968605626,
  29.31365858752018,
  240.87515441295776,
  7.349987229521624,
  8.040903474206418,
  70.495483000691,
  80.04154274069904,
  37.5077812334381,
  14.86049712352624,
  127.51809157802892,
  44.11336346831556,
  42.708811701121235,
  51.84370794664624,
  6.957251596519924,
  7.5633209075184435,
  154.92547772479304,
  65.10510588279544,
  13.970938815513385,
  16.67373622336668,
  13.965002395384998,
  9.774732875713882,
  50.40053174167759,
  20.572659858812724,
  18.5722970477008,
  50.46364338028776,
  86.24458798055095,
  28.18677420682364,
  42.38023466856238,
  137.65555416761748,
  73.24549361865107,
  105.24085406270217,
  10.490086743301685,
  18.412659211927984,
  126.67968922692351,
  190.3264668268053,
  30.780790589445022,
  36.614259470710245,
  89.06058403203807,
  7.807671189160624,
  6.24335672835445,
  56.79672720877181,
  102.91848094951155,
  17.875086605125155,
  7.892021451760044,
  61.01470381868949,
  12.406418082477062,
  6.748779702580254,
  25.958237456887023,
  69.1698950704019,
  45.01203071874414,
  42.944023825377656,
  72.36518828557253,
  8.891312813508083,
  30.96455731915713,
  7.94269144532792,
  23.065254478573618,
  20.788734966781774,
  45.73047587764724,
  215.98970605848322,
  44.93623625691771,
  35.31738814408766,
  30.31201231925092,
  98.93881085236103,
  78.75158026739878,
  23.4595628758251,
  76.77988180322536,
  70.39598984811391,
  66.70019618152155,
  113.14505199289167,
  54.70124195187026,
  52.94349720964987,
  42.33583923855222,
  24.10534463564322,
  72.64840700400005,
  6.522902179507845,
  7.5384245499817455,
  8.841998333315564,
  35.910716564492105,
  7.432011450374535,
  5.953889844981608,
  10.172147895998716,
  21.211584023493852,
  132.64282875677196,
  46.95613120107641,
  13.240362551425536,
  84.05597771212109,
  95.4263989017441,
  47.14345365812882,
  130.13108768075568,
  15.542089658349264,
  78.26597808965136,
  17.58474488272077,
  44.439362343366284,
  9.485707184580157,
  35.39383327926229,
  12.319144334878755,
  73.70620791928121,
  22.92211956799052,
  6.057931523397967,
  31.50087497424812,
  46.71624497856206,
  68.94485331671018,
  303.6892269669984,
  67.70411086818183,
  7.462657720691568,
  123.31614311389951,
  42.44801957799693,
  5.554792542432959,
  24.535392526438656,
  7.619289424757403,
  20.740731402955074,
  55.12146533469664,
  26.2799560566004,
  79.80855636979358,
  56.83126275754098,
  59.354000627724,
  22.723620925129964,
  6.494081828629084,
  64.80388256890106,
  74.0678067262305,
  24.159269801317635,
  25.231914728561563,
  109.73274350698438,
  174.26285065951365,
  11.828181356511971,
  270.53476515284547,
  42.70299646726895,
  5.557147349491229,
  9.89570990903025,
  181.50940305007032,
  15.0233399560187,
  6.299846323086897,
  28.19122365866761,
  23.714868175546975,
  33.38075789436762,
  15.618265247795502,
  116.10276102300307,
  123.09693618944496,
  61.43580605315465,
  75.58702888738935,
  18.45953911076787,
  108.60358975513593,
  10.469569536729665,
  14.61200007351098,
  7.612136344332092,
  58.25369299131724,
  11.669723714965864,
  39.93866567719314,
  64.81639125204522,
  37.3190142815898])

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

([4928112.4989503175,
  5036230.994222343,
  5073247.653318141,
  5157092.537210761,
  5191763.327334791,
  5300340.535763368,
  5317382.550383241,
  5321083.964715635,
  5355676.5625,
  5438311.969999445,
  5451291.820822522,
  6778494.16440898,
  6785040.20378968,
  6961798.4375,
  6961824.3864494115,
  6965671.875,
  6966445.3125,
  6996359.375,
  6998190.625,
  7042068.75,
  7067759.136083348,
  7085209.378132706,
  7139805.228617017,
  7160751.5625,
  7173569.216348684,
  7174431.809798662,
  7183017.109251836,
  7267945.3125,
  7270080.18207427,
  7277856.25,
  7297820.795334424,
  7308457.590299777,
  7325564.0625,
  7338440.044851796,
  7339704.23392173,
  7383595.006030557,
  7408210.089322634,
  7439252.530490686,
  7449281.994137895,
  7450372.632883947,
  7450440.53282894,
  7453479.6875,
  7460293.943414853,
  7461414.0625,
  7462551.5625,
  7464193.726921775,
  7464526.5625,
  7465321.813099671,
  7467959.614432444,
  7501741.049957041,
  7510363.752011988,
  7514546.034548267,
  7520335.081267191,
  7520371.272001179,
  7522658.281426924,
  7524769.998808046,
  7528399.2505654795,
  7529044.037948253,
  7541440.625,
  7546482.8125,
  7546489.0625,
  7547263.463943135,
  7553357.8791160025,
  7561035.9375,
  7577926.50891812,
  7589500.577856496,
  7589963.907652801,
  7591550.633801611,
  7592200.0,
  7592793.412688895,
  7594035.9375,
  7595455.167362315,
  7597749.088604033,
  7597774.573422763,
  7598457.582963947,
  7600586.047011132,
  7600948.586543177,
  7601541.609075746,
  7607356.241776778,
  7608885.987471846,
  7617818.75,
  7623774.413744111,
  7623906.25,
  7628007.522701366,
  7631127.952310165,
  7642312.933400372,
  7642623.4375,
  7646555.927471407,
  7659704.424472292,
  7660423.246718515,
  7661290.781656286,
  7667187.519684521,
  7669906.820943252,
  7670069.640706044,
  7676299.3879393805,
  7677092.784260475,
  7693406.840187491,
  7693458.300954373,
  7695421.661761349,
  7696093.095569686,
  7698765.346409517,
  7713493.75,
  7721913.727597748,
  7722556.980149264,
  7724259.525576042,
  7724515.863845837,
  7724843.75,
  7727261.023425729,
  7727271.356537071,
  7727415.625,
  7727877.194644847,
  7728186.571872057,
  7728515.680084369,
  7729054.871836707,
  7729200.301527607,
  7730979.611331726,
  7731413.585203696,
  7733510.9375,
  7734096.009670198,
  7740732.768387361,
  7741889.0625,
  7741920.126632115,
  7744231.594028677,
  7747357.8125,
  7747757.514116983,
  7749698.4375,
  7751104.6875,
  7752864.455034001,
  7756960.9375,
  7764938.664784911,
  7766651.5625,
  7767120.3125,
  7767438.239700456,
  7767926.5625,
  7768345.144835219,
  7768367.540090685,
  7768904.336399147,
  7768928.964429455,
  7768935.9375,
  7769650.0,
  7770121.04240849,
  7773571.440123174,
  7773998.411369783,
  7774192.110536828,
  7774225.0,
  7774447.749575863,
  7774453.125,
  7775539.0625,
  7777068.605323072,
  7777170.3125,
  7782460.9375,
  7782693.75,
  7783134.910419227,
  7783312.794751378,
  7783749.427500115,
  7784467.062865479,
  7787787.064645441,
  7795668.816698004,
  7795856.25,
  7796655.232118821,
  7797149.858760677,
  7802854.6875,
  7807122.39877029,
  7807973.4375,
  7826756.25,
  7826756.25,
  7842487.447523741,
  7844718.535868823,
  7851614.102293992,
  7865723.294864476,
  7870084.375,
  7874134.848562917,
  7877729.353606684,
  7878374.032340517,
  7890621.875,
  7902452.33385658,
  7921685.9375,
  7921687.5,
  7924925.46120395,
  7938648.198609276,
  7959710.47465506,
  7985421.918472736,
  7992271.766758991,
  7995504.696577904,
  7997695.11200486,
  8000735.9375,
  8023782.600891208,
  8026537.5,
  8026938.284448174,
  8027025.55993296,
  8027387.111576463,
  8027482.8125,
  8028878.228104227,
  8031037.455439818,
  8031115.537842666,
  8033304.051546997,
  8039536.7103334535,
  8043325.73759618,
  8044762.24497101,
  8045268.983490752,
  8045851.5625,
  8048845.415036554,
  8049198.030699028,
  8049636.389559176,
  8050270.3125,
  8050626.5625,
  8050758.249043236,
  8050821.624179483,
  8051181.404633995,
  8051940.625,
  8052812.131277288,
  8053056.706105952,
  8054601.5625,
  8057075.833765266,
  8063735.326451124,
  8066506.25,
  8076059.375,
  8076200.156499177,
  8078877.447056814,
  8079407.79223353,
  8079756.25,
  8081642.1875,
  8085040.697201559,
  8089092.171798846,
  8089559.375,
  8089950.870734111,
  8090286.385033847,
  8090292.187302922,
  8092241.566301165,
  8094739.0625,
  8104568.75,
  8122196.365773838,
  8124370.317114247,
  8145745.3125,
  8153418.084821994,
  8155963.30454557,
  8181686.231968025,
  8186329.011787595,
  8196593.75,
  8199386.539502798,
  8200132.462763735,
  8200161.615539064,
  8201437.270873369,
  8201957.7435368,
  8204474.977320864,
  8207768.75,
  8208186.36743079,
  8208777.336027159,
  8213970.3125,
  8218710.171837039,
  8220077.517073625,
  8220664.455636583,
  8220664.916431453,
  8233194.623744854,
  8240549.360318952,
  8242492.1875,
  8251937.5,
  8270687.059177668,
  8279017.572944678,
  8287757.8125,
  8294750.347112862,
  8297646.971576843,
  8300795.3125,
  8302445.316396724,
  8305570.725915888,
  8309044.852922054,
  8311505.633486733,
  8312116.216577061,
  8312432.117639927,
  8313599.753057187,
  8313651.069343612,
  8314070.3125,
  8314099.423960631,
  8315051.476407381,
  8315123.4375,
  8316479.6875,
  8316499.716329852,
  8317019.467708244,
  8317121.875,
  8317151.5625,
  8317175.684282461,
  8318581.16909885,
  8319086.0302196145,
  8319989.0625,
  8319990.625,
  8320053.968010188,
  8320442.108191726,
  8321211.738216368,
  8322905.739160612,
  8324197.45631981,
  8326223.3470792435,
  8327634.02880927,
  8340076.558805409,
  8344870.6528990865,
  8345062.327534265,
  8345539.155398426,
  8345935.26695746,
  8346270.3125,
  8348530.291463457,
  8350254.5892514,
  8351448.960669997,
  8352077.085393564,
  8353981.25,
  8354491.011043467,
  8354655.271274155,
  8354942.961162562,
  8357061.836866527,
  8359675.327558865,
  8361828.597066413,
  8362179.213305795,
  8362197.206458042,
  8362717.1875,
  8364066.463802592,
  8365800.331625924,
  8366057.262007171,
  8366451.360100323,
  8371052.421910109,
  8371646.3900579205,
  8372870.3125,
  8373288.062428359,
  8373307.968666479,
  8376871.875,
  8377524.114762574,
  8377525.984602657,
  8385687.1533857025,
  8392364.0625,
  8393432.418561727,
  8393451.884448973,
  8393855.74306448,
  8395354.010603042,
  8395506.202094821,
  8398071.654072808,
  8399175.584507551,
  8399606.085573912,
  8400477.59843692,
  8402713.808092939,
  8403599.736644963,
  8404043.75,
  8405123.196955431,
  8406294.64174241,
  8406318.75,
  8406332.8125,
  8410882.629978953,
  8413853.098848302,
  8414165.443707367,
  8414224.22646879,
  8416485.9375,
  8417261.647422258,
  8451413.919805456,
  8457365.400593042,
  8475653.91058664,
  8478498.4375,
  8482995.555825783,
  8485418.75,
  8505571.881481376,
  8513292.524653705,
  8513342.124890037,
  8516474.98144081,
  8520664.052955315,
  8556786.386712432,
  8556881.25,
  8566668.574877169,
  8571082.239563277,
  8574679.6875,
  8584252.705169203,
  8594757.8125,
  8598044.548606426,
  8601592.1875,
  8615725.749367699,
  8617676.989804832,
  8617726.395943604,
  8623713.128053546,
  8625082.919685105,
  8625810.646208726,
  8626141.722661693,
  8626164.0625,
  8626190.625,
  8626262.240024853,
  8627054.216234744,
  8628262.370622981,
  8638746.937816445,
  8645048.4375,
  8656671.875,
  8657466.76109016,
  8666741.42635224,
  8668084.375,
  8674427.819685213,
  8688781.629497789,
  8689140.83092631,
  8689870.3125,
  8690026.5625,
  8690426.978738997,
  8690464.0625,
  8690515.305800967,
  8690779.494858652,
  8691157.910657002,
  8695473.283444319,
  8697453.125,
  8699466.692589542,
  8699593.394955164,
  8700452.849068219,
  8706508.862549247,
  8707676.5625,
  8707938.553769981,
  8708731.245737484,
  8709403.079703972,
  8712004.6875,
  8713247.66956024,
  8714228.125,
  8716215.202061305,
  8721128.125,
  8725490.625,
  8728900.0,
  8731607.02811587,
  8738446.875,
  8768687.891110368,
  8781291.571391275,
  8783153.966033926,
  8796893.498754788,
  8797083.088780563,
  8800035.590136496,
  8800468.54824382,
  8800645.3125,
  8800786.639872646,
  8803204.53056855,
  8806076.585097393,
  8806084.272478458,
  8806703.66117119,
  8807775.0,
  8808114.0625,
  8808526.279198207,
  8809034.375,
  8826511.516911607,
  8828482.717357265,
  8828489.595247788,
  8841285.9375,
  8855067.1875,
  8859428.125,
  8864448.841614436,
  8867828.028597258,
  8868305.222270196,
  8868792.1875,
  8904137.5,
  8927442.397498403,
  8939105.686883226,
  8957125.290068386,
  8959234.375,
  8974734.053802043,
  8986701.533118851,
  8987067.1875,
  8996024.995150551,
  8996598.4375,
  9003556.240705365,
  9036503.125,
  9037611.849858824,
  9041590.625,
  9042186.222912796,
  9046012.176946638,
  9077160.120997218,
  9077220.47831315,
  9100035.9375,
  9111443.176914165,
  9158772.432460576,
  9178772.941857023,
  9219050.940352,
  9227913.652766291,
  9235752.275383314,
  9239427.7062371,
  9263566.711817905,
  9270237.33245167,
  9274719.662006212,
  9327288.211935423,
  9328164.0625,
  9399050.868746758,
  9399291.926408896,
  9601913.09858774,
  9602331.25,
  9644308.34937074,
  9649944.110416917,
  9757497.128231078,
  9806924.910490336,
  9891566.249271,
  9899257.730456265,
  9899263.942708107,
  9951846.424792135,
  9955604.389635475,
  10079533.060011435,
  10081965.172663918,
  10095268.090938635,
  10096637.727239292,
  10136351.82051684,
  10136363.62075997,
  10148972.433423938,
  10149559.282204537,
  10152686.589196257,
  10208654.366364075,
  10212271.875,
  10277493.75,
  10366323.96732527,
  10366986.022250822,
  10389824.1905386,
  10449226.54390808,
  10481341.530068649,
  10647001.5625,
  10655455.784229299,
  10712977.082335524,
  10802095.210919635,
  10890344.431180144,
  10906858.300711997,
  10927722.133428497,
  10929320.72471968,
  10967161.364347601,
  10972091.964679774,
  10974415.290450843,
  11032509.308349539,
  11034975.615031825,
  11099318.75,
  11145931.23604667,
  11145971.925563348,
  11152262.48499993,
  11167701.5625,
  11200131.755921412,
  11263233.711107356,
  11300817.373708166,
  11351764.62703371,
  11386057.086578716,
  11403738.595086506,
  11450519.71053722,
  11559345.3125,
  11571109.375,
  11683712.169916466,
  11689170.052798927,
  11718090.625,
  11796484.375,
  11798926.250507787,
  11799317.33510752,
  11855752.820890762,
  11855753.125,
  11861334.375,
  11904612.5,
  11905260.803323846,
  11940785.999630323,
  11957873.630855177,
  11968992.015292574,
  11986109.607981514,
  11988990.556150988,
  12000330.865965478,
  12000721.419723637,
  12003637.5,
  12013730.526579404,
  12025076.923100617,
  12031860.9375,
  12044717.114371343,
  12044832.262686346,
  12044866.258667452,
  12051963.891360208,
  12053504.222162308,
  12065428.700254148,
  12100260.457968684,
  12100428.875834126,
  12100521.613606047,
  12101031.944963839,
  12101641.736556266,
  12102864.0625,
  12103403.125,
  12106042.1875,
  12106053.524926012,
  12107975.885402422,
  12111437.5,
  12112353.532229342,
  12113627.02966085,
  12113981.661584968,
  12127799.623681167,
  12129737.890908703,
  12130976.18025973,
  12131888.912029162,
  12132723.133956652,
  12136144.492694957,
  12137367.321974728,
  12156985.9375,
  12178909.156252343,
  12180044.999346372,
  12182261.808479087,
  12206553.75725092,
  12213708.988271683,
  12231860.9375,
  12233003.931692105,
  12233992.1875,
  12235307.877462111,
  12236860.9375,
  12243556.25,
  12246220.337330703,
  12246220.98635424,
  12250039.0625,
  12251293.35837731,
  12270825.417941313,
  12315375.06569629,
  12316543.962705411,
  12323460.32709706,
  12324313.384600127,
  12334077.192554666,
  12336724.173568392,
  12337765.625,
  12341471.169248201,
  12353645.882447833,
  12360157.411962524,
  12360671.875,
  12361101.481186591,
  12365635.9375,
  12365693.50405218,
  12368671.148806043,
  12369346.875,
  12369726.052549252,
  12373270.876579523,
  12385782.544451762,
  12388402.343412558,
  12388435.263011193,
  12388565.625,
  12392239.753562465,
  12392659.682422675,
  12407177.137471808,
  12408483.046715355,
  12410066.679712268,
  12412177.893814208,
  12413194.53484666,
  12415309.487116607,
  12454878.125,
  12454878.125,
  12454884.113512365,
  12464404.17052088,
  12475781.25,
  12482887.5,
  12491459.017914796,
  12496264.890647713,
  12496564.0625,
  12497086.08039192,
  12510701.5625,
  12518944.5534123,
  12558012.231662245,
  12602168.205628183,
  12617085.9375,
  12617875.282408794,
  12622629.427306714,
  12625505.196172269,
  12652141.976454712,
  12670826.5625,
  12675375.64480423,
  12724159.375,
  12729757.8125,
  12731730.560226863,
  12761296.875,
  12762276.5625,
  12766442.1875,
  12767004.6875,
  12767637.776918728,
  12768873.4375,
  12773937.5,
  12775829.787915375,
  12776082.8125,
  12799215.625,
  12801985.67810539,
  12801996.664054416,
  12809708.845568748,
  12809841.287454512,
  12810412.108728705,
  12812492.666752644,
  12814362.4708218,
  12841928.84920282,
  12842684.698368622,
  12842960.9375,
  12842982.630621236,
  12843374.478768272,
  12843565.961579485,
  12851479.622186773,
  12854739.106796453,
  12860861.165745437,
  12861768.575807039,
  12879271.875,
  12883087.5,
  12896689.752215723,
  12898294.711687803,
  12910569.939069323,
  12926728.088782512,
  12944428.125,
  12945717.1875,
  12975956.235638129,
  12976982.8125,
  12977962.23576663,
  12978343.274225801,
  12978397.880479556,
  12979221.115958042,
  12979722.942118015,
  12988841.654664552,
  12996532.646194644,
  13003348.417378366,
  13005721.85358186,
  13020053.125,
  13026039.548362201,
  13047723.868260026,
  13049032.8125,
  13061197.773720844,
  13092529.575770812,
  13101117.1875,
  13139054.500913948,
  13141177.844707327,
  13142220.703160131,
  13144174.71781907,
  13149697.356564377,
  13160144.310996471,
  13161296.277499402,
  13178764.54825073,
  13180466.886455594,
  13184084.375,
  13184859.357688708,
  13192573.4375,
  13217912.5,
  13218251.5625,
  13228065.24626964,
  13250096.424802622,
  13250944.436122917,
  13251824.859801963,
  13252399.885734472,
  13254709.375,
  13258878.149024663,
  13264669.392847164,
  13291125.531940714,
  13318093.55893475,
  13436100.675369803,
  13479318.05163303,
  13512222.468342818,
  13514127.296096109,
  13549889.0625,
  13601174.125953637,
  13609099.9257834,
  13615434.23017301,
  13616828.125,
  13618156.113376312,
  13649542.772527138,
  13690618.791648932,
  13693544.741738344,
  13756411.159074683,
  13757289.912697056,
  13757604.33187155,
  13758779.005601779,
  13808909.375,
  13860483.071516888,
  13905126.390576802,
  13909102.015869262,
  13930563.39813776,
  13984470.167927423,
  13985185.9375,
  13986047.26349431,
  13994370.3125,
  14034506.25,
  14038898.12181526,
  14040118.75,
  14084625.0,
  14172919.38951191,
  14180237.5,
  14214701.28976263,
  14214884.375,
  14272776.78720475,
  14282157.629674047,
  14285638.420834383,
  14294968.540596422,
  14295910.437899055,
  14296434.53983261,
  14298643.038561702,
  14298790.096308548,
  14299470.3125,
  14301478.144789785,
  14302721.343413487,
  14304503.691940747,
  14304682.8125,
  14306492.1875,
  14308017.529457469,
  14309509.524124196,
  14351395.3125,
  14437879.802880725,
  14438519.805870693,
  14438562.677019713,
  14567965.807816833,
  14568380.335978758,
  14568715.625,
  14569565.060232587,
  14570002.829817912,
  14570965.139149569,
  14571254.6875,
  14571320.3125,
  14572818.508692514,
  14572823.4375,
  14572846.468839142,
  14572974.911198325,
  14573925.0,
  14575331.689392272,
  14576285.69654126,
  14576450.0441637,
  14576607.233619459,
  14577712.5,
  14578392.567612343,
  14578404.021282103,
  14578774.770167405,
  14579862.5,
  14579896.490098722,
  14579900.044450147,
  14580976.5625,
  14581135.9375,
  14581193.481323069,
  14581856.477213265,
  14582081.25,
  14582388.566443786,
  14582567.593741134,
  14582857.726646248,
  14583358.131934201,
  14583628.611107716,
  14583633.412135994,
  14585243.69024646,
  14586834.396304809,
  14951506.504111039,
  17830421.857919473,
  18750417.165826026,
  18751571.54319177,
  18760962.66690401,
  18766547.148982238,
  18779312.5,
  18781648.4375,
  18787685.9375,
  18791538.252323695,
  18791776.335534506,
  18791812.03376177,
  18791909.65048032,
  18793443.77453665,
  18797495.3125,
  18810867.508635696,
  18815743.75,
  102069009.375,
  102069040.625],
 [51.019727640144815,
  24.332549071475135,
  17.3993505790785,
  23.14311441386546,
  19.429381360668685,
  77.3385001613488,
  8.152664579600076,
  58.37078312040285,
  30.14550252492934,
  21.19895951299102,
  13.70664402501979,
  6.3210152915565105,
  94.72844267427124,
  37.2287108097906,
  13.025456301338632,
  108.1242463352786,
  50.00846480260808,
  66.96792301048697,
  54.69352784626684,
  85.43212882556722,
  33.071288895178554,
  23.495822414529655,
  13.551066366682026,
  64.49609518319588,
  7.781090731727417,
  5.486317898114985,
  46.04963558183315,
  38.32516084582863,
  52.748205806823364,
  72.58054141015513,
  17.16826309144363,
  33.00985311964993,
  34.37930580246373,
  109.24430228640362,
  16.253474263773896,
  50.68560737194869,
  5.753779302985908,
  27.654414848409456,
  103.96772557390341,
  50.2938702629432,
  5.8071917286905945,
  39.4497157115031,
  9.500482202494188,
  41.911703494531245,
  49.5956397582072,
  28.07611590737734,
  33.94616425981533,
  16.604391903071047,
  17.187940123130844,
  87.39315844690233,
  37.78370264227462,
  45.275673183013,
  8.699063203928658,
  8.868758347081377,
  60.39102764820626,
  23.577407407889883,
  12.691150566274143,
  90.51423539736149,
  38.11950288442491,
  76.11552355717849,
  116.32668504918468,
  7.993531997877487,
  15.201646652661815,
  92.24377710007931,
  19.27012507283209,
  17.72135863534364,
  10.893999760769823,
  8.65281675252923,
  98.1665540285828,
  23.285563504799256,
  37.07183535066272,
  8.557757016071566,
  8.663769570181527,
  23.617375709220394,
  7.0623931036725125,
  15.37601623390135,
  9.46068948486128,
  23.84301158856338,
  37.5487554006468,
  30.453815985959146,
  60.65921914809008,
  33.74751468190908,
  65.14038141773567,
  7.501680854975588,
  10.868099725358467,
  11.880200447835133,
  91.91207700294838,
  13.83988038088738,
  63.47924516550697,
  64.5213682036561,
  37.652879218916944,
  15.01707217921971,
  159.9492055892424,
  126.51231125879093,
  128.5432534119197,
  43.00794538026564,
  17.796648601868746,
  5.0778876430670365,
  24.92609179782862,
  10.238694148450953,
  86.9861887201475,
  113.97306082189817,
  5.1576346291734545,
  34.532664230909205,
  14.740611947767835,
  12.656350278708338,
  98.36547058097557,
  38.21521333248496,
  35.74932402482466,
  84.18974668802073,
  18.05929768248392,
  26.247122385519162,
  34.26502670903773,
  24.12420041879301,
  8.452737341847591,
  27.10366468398521,
  13.388659198874077,
  30.530054930213357,
  6.6022294621061794,
  24.209458760469037,
  105.36106381599684,
  7.073435009116808,
  7.214921408001459,
  52.90917582199029,
  121.58097718475993,
  81.66669052107775,
  46.29084546423743,
  17.777123146793024,
  91.8272917407803,
  29.63514990032996,
  60.5866002291116,
  146.131701741978,
  122.1088564795227,
  39.62207490400502,
  7.543701344727414,
  82.07562207406878,
  23.119231456509876,
  10.1587533115775,
  64.49473409828599,
  62.6931274471124,
  84.76796646393606,
  22.168305240012284,
  41.639322451673316,
  27.08970949779508,
  150.19195903398287,
  10.199413647942263,
  35.925234655776514,
  32.36118581894197,
  25.960615229470005,
  77.95496600398273,
  44.164314970193146,
  35.706978293407325,
  123.56040923114173,
  16.198078217869536,
  98.96520163382988,
  14.486828091488016,
  60.545482158006024,
  15.201974732971163,
  64.94936796020671,
  17.52255030025914,
  75.36788066337975,
  50.426726015037644,
  12.058613132253763,
  51.58431372510711,
  29.20118300661746,
  53.823589298265176,
  88.49258621636241,
  10.808726491800996,
  40.11516331888909,
  7.084455188617361,
  71.97939260395822,
  24.33171868404032,
  83.38245595902278,
  12.20496382123465,
  112.14656948161337,
  32.290526419026996,
  41.54608261342369,
  39.74803956805752,
  8.546518946424731,
  24.564027504530284,
  13.428970051555387,
  92.29854904916911,
  14.544394094778951,
  77.53345364600796,
  10.828143405788646,
  29.89399343235969,
  14.191961316929504,
  55.866106359313335,
  79.02966305943522,
  17.613411792059154,
  44.527727939597426,
  49.546152580837735,
  15.349617343751005,
  5.611583711422576,
  72.04395976547586,
  41.558567489368286,
  16.79906812359159,
  36.33295070169201,
  24.342699166062744,
  11.51339443961825,
  89.5195259055385,
  6.280240627168779,
  36.81996218688203,
  73.1749559687757,
  126.99589923332962,
  38.517865131338944,
  9.729979054268567,
  36.43206622709621,
  15.565585541765214,
  57.21913455550318,
  20.291031351152636,
  15.573674641968893,
  61.963341078082934,
  35.417597386158754,
  20.452516666547517,
  43.00292569181489,
  38.50699852823647,
  6.795157854379685,
  19.88445646071125,
  9.314325319610372,
  48.31806337084415,
  58.174144359128334,
  6.7780071671085995,
  37.588320043095436,
  105.82759213376322,
  58.91739354815732,
  11.90755762656985,
  9.331236539073581,
  88.75385233394404,
  60.935934874920704,
  51.55023971996729,
  28.595339460497833,
  64.62751397591268,
  46.04628958123016,
  12.895336789285572,
  6.93551065607059,
  8.448687432769841,
  96.57061404180752,
  47.4306482397842,
  16.56090360167681,
  157.3154259000464,
  18.358585003645082,
  24.700330971229143,
  58.38571643423502,
  14.855766069187162,
  48.845474273326026,
  73.2327603277362,
  115.47797534579206,
  73.81302802249368,
  31.502367548705436,
  13.039870698759033,
  16.686523897456144,
  5.576011305599249,
  12.837201418512956,
  81.36422136917308,
  52.58972869870623,
  89.50069064146214,
  28.343329062268463,
  16.64689156734508,
  74.77187712224105,
  28.187122455315027,
  5.821160852366742,
  103.52608515239372,
  14.956382739687166,
  58.46735927881736,
  10.478071806390014,
  34.12020169574375,
  67.05192816380818,
  19.747087915256788,
  14.077005818968322,
  52.48166922318646,
  67.5353250314676,
  31.507465028817442,
  6.473679617656323,
  57.58220980417575,
  89.70819536246653,
  169.9679634212193,
  25.568063759888183,
  31.168903674726693,
  43.60001106766382,
  9.98580541850406,
  27.020826201050195,
  49.31325780256245,
  80.1229596977723,
  289.73522540853793,
  19.55137396986149,
  82.05906374272587,
  8.760847838521935,
  8.933838923800112,
  19.54144890443359,
  6.473009789596844,
  23.17812158002227,
  40.16695747824426,
  62.45388966852825,
  69.32183867517249,
  9.411209049278996,
  196.72001821701687,
  75.2678657439432,
  17.88176923711083,
  26.722283778334162,
  8.58407423040281,
  17.39576355138891,
  39.08570705675115,
  5.97374468412474,
  17.298529727484897,
  19.02610162508597,
  68.96127839366913,
  6.83071188075995,
  36.652314223790974,
  73.75273602875131,
  27.96032714660041,
  63.92405419090463,
  12.75248054857289,
  35.75672133830594,
  77.52945493155681,
  25.083002825376514,
  41.85084021719274,
  21.644637769226446,
  39.86130097415324,
  41.34075203161305,
  6.795089692056224,
  41.61269582831005,
  28.765265057211426,
  13.11764922052998,
  19.74670148021197,
  43.419051781772204,
  74.61037120517116,
  126.69806054880688,
  45.60304118376206,
  15.036963766866599,
  27.611222225785276,
  101.21917185698308,
  170.53958271885588,
  5.69637279703897,
  16.43692663407476,
  128.7327693257842,
  30.550187030731397,
  37.06429330512138,
  7.059623543790945,
  19.94354879647016,
  77.25456995278871,
  59.02931626383904,
  25.377080217568487,
  50.08824868493095,
  62.04801495197008,
  28.653259257672243,
  61.94219606766775,
  18.574306011650485,
  32.19927203024752,
  89.04931689868317,
  17.973802436227004,
  34.557040364572316,
  44.11583721527938,
  51.89619477419498,
  6.26801758963527,
  38.00504178398391,
  6.751071017579627,
  20.425796054471213,
  53.13186437973502,
  16.983349901710593,
  81.14156649513676,
  9.187590237442256,
  15.74676066448239,
  46.161293410853844,
  46.88810232052171,
  43.339594208611864,
  19.210809412392116,
  57.350078578848866,
  49.055033447814246,
  9.720886433989314,
  34.592547228103086,
  104.39216011951194,
  9.423709051788068,
  8.138319032084125,
  79.17389568270941,
  46.43790197491268,
  46.22913153699679,
  93.01034920214863,
  13.40543875425305,
  138.54015206748417,
  8.320077423236954,
  36.22645147310939,
  70.50010515108235,
  19.530338328696665,
  10.115334088168769,
  62.181391195169525,
  45.52540581362601,
  49.76706700360043,
  16.504915815964594,
  95.00847027972847,
  43.336059585196594,
  6.55438182335269,
  81.50710199820907,
  23.50855982177313,
  141.13981153732965,
  53.10011727803297,
  10.89515347225526,
  62.70856628020768,
  21.544021814837333,
  13.470764836514936,
  13.984676876953692,
  49.422938205509695,
  78.7000462090955,
  10.415875445354137,
  102.2857716787643,
  27.65922151046059,
  42.3716698825728,
  8.75676749649084,
  85.77456054525874,
  22.29788730141916,
  32.29827523355036,
  85.9671628412062,
  41.56562198283246,
  72.42087838360187,
  30.61119760110114,
  9.685556521975242,
  14.755553475462332,
  26.679464583617175,
  29.65451511680259,
  17.718140344932806,
  13.764468904179129,
  21.997521845675937,
  60.50736150319451,
  12.437255890573878,
  56.765748864475405,
  38.30464209403545,
  43.520394003074514,
  63.68306246500674,
  64.97735108662387,
  40.04895538926286,
  23.966407554350187,
  48.23722491285364,
  69.93308792635801,
  26.775155563154645,
  6.336938358068589,
  48.59318691944191,
  174.550736392183,
  32.16997885619041,
  9.698196482138107,
  14.586371975303496,
  17.464770863258533,
  113.51018827508985,
  45.95399153613879,
  5.406563215586613,
  8.00966667388881,
  11.613536822909879,
  98.60064988813923,
  35.3923526281148,
  102.09263664234726,
  47.052096592033905,
  34.06027987811051,
  71.91287033111117,
  21.989748455857416,
  68.63860813273082,
  20.115613083047016,
  32.53027553404104,
  59.52971022420215,
  6.700680913799637,
  18.826533941191574,
  5.5411824603603,
  58.29512146299198,
  12.07121727270243,
  26.36180633733093,
  21.536674253680363,
  42.78802666586461,
  22.443431894950514,
  14.166873197157742,
  38.449109994219825,
  10.45755413273848,
  35.7064020383829,
  69.01427769538401,
  31.991424973693558,
  63.19702396379522,
  5.149513407433575,
  23.95635927513694,
  17.848425332107094,
  78.3588695956492,
  85.16665870520886,
  14.92361440479509,
  5.408243866311126,
  10.990621561758923,
  5.2196236460076015,
  6.474306608463468,
  45.62299201348876,
  5.589613802676459,
  61.91126842023762,
  135.98364894085228,
  21.9821258498567,
  20.04070188231627,
  8.410515552246348,
  78.65153901800285,
  118.7516431882548,
  17.604402922112953,
  8.238107655226662,
  9.950777944615773,
  9.073864876870976,
  104.75841511247097,
  43.28001990608811,
  5.741328130287127,
  12.326915406338511,
  28.82744169361292,
  11.098770257185098,
  21.726454783274562,
  79.44354172075879,
  6.040966413788864,
  17.378112574962653,
  13.364271904134796,
  65.03925902352591,
  88.69967013226335,
  5.62850281495974,
  8.517298132653174,
  33.88675751084155,
  10.796414269669068,
  30.46945185098557,
  5.067686295655706,
  12.24755992251158,
  42.77836747743613,
  55.66029067334087,
  11.20767114841728,
  158.96983506909697,
  35.772958690158,
  9.79925628108159,
  53.194216105796144,
  17.066650233275546,
  336.2615048333812,
  78.73451572037436,
  90.67608730380516,
  58.13111133272538,
  72.515178390632,
  90.4349291476602,
  5.235999902591054,
  19.277354771767495,
  98.264141736135,
  55.038779150926324,
  5.092364745993398,
  15.53563922832179,
  49.238163393636356,
  60.88776225140054,
  36.36509806061902,
  79.92001444914862,
  152.65141385148377,
  43.59354606300749,
  41.02063172511431,
  14.323391972178356,
  6.297089594210837,
  79.39886730544517,
  22.448726883922575,
  16.815763898354145,
  189.36239127126672,
  28.3736407235359,
  16.52060297249311,
  30.501089602136993,
  27.157186216090633,
  13.191657653892507,
  16.681653252502226,
  25.54923416537571,
  83.00748969813924,
  6.342720348008545,
  210.96699658732354,
  27.331477538838218,
  296.0380689426962,
  64.11204694029546,
  8.061686986837433,
  32.721722216808104,
  73.0644645775102,
  71.531021667596,
  11.816188405022837,
  100.0884314945292,
  67.81183538021493,
  104.24922352092497,
  108.54316587093628,
  16.782626760796155,
  13.678042995397185,
  136.41347889139655,
  22.75416084621835,
  25.8691869168368,
  5.3869952720990035,
  24.251599208985276,
  15.480094363000273,
  76.00058134829126,
  24.847246620879215,
  9.07994409991257,
  19.591710175419557,
  34.841003179280875,
  9.024494425834693,
  65.44328544121508,
  28.200095601828465,
  42.887748280603496,
  112.74417763732102,
  101.71799031821259,
  60.895941456786225,
  26.60800007532731,
  68.88621083963525,
  32.343477556961105,
  8.089797274311929,
  8.520690344654893,
  36.65291372550179,
  9.51797545996134,
  52.841322147188706,
  107.16929798435424,
  37.079805535883835,
  12.512894379241473,
  68.62619829987139,
  60.28862685863818,
  240.58579977733874,
  89.5225707452022,
  51.1643796027628,
  40.20264865868681,
  361.0020765450235,
  5.159672925853372,
  19.49895894680683,
  43.75040927438075,
  21.394407234518805,
  12.140944186867133,
  24.167645551223593,
  70.27009267328918,
  19.912808346496266,
  46.178156138228466,
  5.5675870871275634,
  18.295071351530858,
  17.703866855820305,
  65.47418421892348,
  18.113773588390636,
  74.89261844354837,
  19.108348713502235,
  16.129959487475794,
  34.057766171151066,
  68.70058283187828,
  24.255473221831394,
  55.260548064446,
  51.18031720735816,
  111.21770153649427,
  218.32283639086896,
  26.419700422123206,
  106.41117905135975,
  5.715912782116807,
  58.20947271200813,
  177.26046924695652,
  14.008693986519727,
  19.855863704175423,
  43.95842429372984,
  5.299891110022268,
  39.94726581426675,
  5.177379979994608,
  39.26319285419514,
  29.847025205830793,
  21.14429392689243,
  33.68065862474706,
  45.38890702766213,
  18.390976220322482,
  112.66360130545196,
  40.622440047755205,
  56.935298826492335,
  124.03833796125484,
  11.59001346102021,
  110.03837508930822,
  192.89958454632801,
  15.338575825302389,
  84.86314336068006,
  62.589779733874565,
  24.283656729927518,
  5.798535246139598,
  10.390055786107476,
  24.60512657208503,
  52.35341631273806,
  5.306124199339385,
  27.9601791132417,
  8.710662849406965,
  5.727651826563023,
  57.482683705852224,
  31.181358453263883,
  39.74587185312034,
  47.311358027807074,
  8.252337973574988,
  30.315260904229632,
  185.9137621271206,
  76.28900853866611,
  38.64456171959087,
  32.84541818743825,
  12.409204573314007,
  65.68369510997394,
  10.59020727466267,
  85.67716011457915,
  57.88003468659261,
  46.17763545910871,
  24.758403228998503,
  30.343499575477843,
  30.028883968605626,
  29.31365858752018,
  240.87515441295776,
  7.349987229521624,
  8.040903474206418,
  70.495483000691,
  80.04154274069904,
  37.5077812334381,
  14.86049712352624,
  127.51809157802892,
  44.11336346831556,
  42.708811701121235,
  51.84370794664624,
  6.957251596519924,
  7.5633209075184435,
  154.92547772479304,
  65.10510588279544,
  13.970938815513385,
  16.67373622336668,
  13.965002395384998,
  9.774732875713882,
  50.40053174167759,
  20.572659858812724,
  18.5722970477008,
  50.46364338028776,
  86.24458798055095,
  28.18677420682364,
  42.38023466856238,
  137.65555416761748,
  73.24549361865107,
  105.24085406270217,
  10.490086743301685,
  18.412659211927984,
  126.67968922692351,
  190.3264668268053,
  30.780790589445022,
  36.614259470710245,
  89.06058403203807,
  7.807671189160624,
  6.24335672835445,
  56.79672720877181,
  102.91848094951155,
  17.875086605125155,
  7.892021451760044,
  61.01470381868949,
  12.406418082477062,
  6.748779702580254,
  25.958237456887023,
  69.1698950704019,
  45.01203071874414,
  42.944023825377656,
  72.36518828557253,
  8.891312813508083,
  30.96455731915713,
  7.94269144532792,
  23.065254478573618,
  20.788734966781774,
  45.73047587764724,
  215.98970605848322,
  44.93623625691771,
  35.31738814408766,
  30.31201231925092,
  98.93881085236103,
  78.75158026739878,
  23.4595628758251,
  76.77988180322536,
  70.39598984811391,
  66.70019618152155,
  113.14505199289167,
  54.70124195187026,
  52.94349720964987,
  42.33583923855222,
  24.10534463564322,
  72.64840700400005,
  6.522902179507845,
  7.5384245499817455,
  8.841998333315564,
  35.910716564492105,
  7.432011450374535,
  5.953889844981608,
  10.172147895998716,
  21.211584023493852,
  132.64282875677196,
  46.95613120107641,
  13.240362551425536,
  84.05597771212109,
  95.4263989017441,
  47.14345365812882,
  130.13108768075568,
  15.542089658349264,
  78.26597808965136,
  17.58474488272077,
  44.439362343366284,
  9.485707184580157,
  35.39383327926229,
  12.319144334878755,
  73.70620791928121,
  22.92211956799052,
  6.057931523397967,
  31.50087497424812,
  46.71624497856206,
  68.94485331671018,
  303.6892269669984,
  67.70411086818183,
  7.462657720691568,
  123.31614311389951,
  42.44801957799693,
  5.554792542432959,
  24.535392526438656,
  7.619289424757403,
  20.740731402955074,
  55.12146533469664,
  26.2799560566004,
  79.80855636979358,
  56.83126275754098,
  59.354000627724,
  22.723620925129964,
  6.494081828629084,
  64.80388256890106,
  74.0678067262305,
  24.159269801317635,
  25.231914728561563,
  109.73274350698438,
  174.26285065951365,
  11.828181356511971,
  270.53476515284547,
  42.70299646726895,
  5.557147349491229,
  9.89570990903025,
  181.50940305007032,
  15.0233399560187,
  6.299846323086897,
  28.19122365866761,
  23.714868175546975,
  33.38075789436762,
  15.618265247795502,
  116.10276102300307,
  123.09693618944496,
  61.43580605315465,
  75.58702888738935,
  18.45953911076787,
  108.60358975513593,
  10.469569536729665,
  14.61200007351098,
  7.612136344332092,
  58.25369299131724,
  11.669723714965864,
  39.93866567719314,
  64.81639125204522,
  37.3190142815898])

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)