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

([2128322.8820732716,
  2673853.9038674626,
  2786856.25,
  2789778.125,
  2805233.259980844,
  2807010.716995365,
  2808506.6130725285,
  2812117.506482652,
  2874049.8859884213,
  2889719.5982088023,
  2889974.2945079044,
  2890028.125,
  2891280.9293882987,
  2892351.5625,
  2893143.3147037867,
  2905490.5032702447,
  2905515.9204744357,
  2906018.734324037,
  2948343.75,
  3034213.041952173,
  3038439.0625,
  3141742.1875,
  3424040.625,
  6785328.776244561,
  6785493.381396405,
  6791056.25,
  6795218.742225575,
  6803876.5625,
  6861424.907240535,
  6866769.9989933055,
  6866830.420718658,
  6879946.775391034,
  6940407.8125,
  6953758.444825136,
  6977915.625,
  6978375.856584645,
  6982520.137476653,
  6982803.634441398,
  6987272.607001323,
  6988256.8744002385,
  6988266.546827361,
  6988965.625,
  6990647.289649813,
  6992557.203741236,
  6993033.774311789,
  7020286.637365565,
  7110217.1875,
  7114679.595704689,
  7115817.552414775,
  7116227.676238994,
  7117364.95151519,
  7118310.367652139,
  7118339.0625,
  7122178.9844065495,
  7123494.955652764,
  7125337.107539448,
  7125688.541932925,
  7129987.5,
  7130548.4375,
  7146026.817144993,
  7157582.864702709,
  7169130.601346257,
  7172898.599309346,
  7173441.9109135885,
  7174609.7948041605,
  7176262.816393783,
  7188563.677035553,
  7192779.81081335,
  7193623.024825048,
  7193984.352505192,
  7199272.54899749,
  7199509.375,
  7210189.0625,
  7219604.912121246,
  7220298.4375,
  7220809.009714929,
  7255282.293938559,
  7269896.415108262,
  7292217.873084736,
  7301620.525656571,
  7303248.4375,
  7306081.167012519,
  7308306.903424871,
  7308351.252823675,
  7310126.5625,
  7310192.1875,
  7312481.25,
  7318070.146381164,
  7318290.625,
  7320417.1875,
  7331122.39212852,
  7351126.5625,
  7352783.587341864,
  7352798.284965201,
  7357039.0625,
  7366504.6875,
  7368967.886011667,
  7369167.9398481585,
  7369285.676346161,
  7369331.25,
  7369440.2057869695,
  7369593.75,
  7370225.326871223,
  7370403.450339849,
  7370856.366219008,
  7371115.00902001,
  7371560.9375,
  7372068.883017307,
  7375343.5475491695,
  7375645.3125,
  7376091.551598742,
  7376977.478443644,
  7377316.100927785,
  7377719.397872175,
  7379546.887680742,
  7380329.11404829,
  7380764.708721365,
  7380780.7633840125,
  7380978.125,
  7381117.873272525,
  7381319.525419512,
  7381475.772936901,
  7381608.512732085,
  7381697.844236114,
  7382160.9375,
  7382195.188103089,
  7382382.979949364,
  7383914.557640953,
  7384037.262309488,
  7384590.995170809,
  7385517.137191394,
  7386523.834975215,
  7386632.8125,
  7407451.204093161,
  7411083.48540908,
  7416076.917212375,
  7420263.841231703,
  7421003.125,
  7423408.618539594,
  7424423.0467876475,
  7424773.4375,
  7443801.285092924,
  7473906.030115981,
  7498054.6875,
  7498323.4375,
  7514310.9375,
  7514313.62134014,
  7515864.0625,
  7517616.54779265,
  7517636.413010648,
  7517884.767175848,
  7518367.767904604,
  7518534.561223161,
  7518955.476446981,
  7519572.428332865,
  7520449.367957458,
  7520524.0528527545,
  7523556.444266042,
  7524612.658927334,
  7529778.125,
  7531141.88947555,
  7531305.553590942,
  7532726.177962356,
  7532863.096114185,
  7532898.4375,
  7533681.611589373,
  7538795.601631367,
  7540046.136425193,
  7553094.853595726,
  7553971.859383369,
  7566477.775048121,
  7668799.085138062,
  7671643.07313268,
  7676275.0,
  7676429.536808817,
  7676441.053957746,
  7676862.5,
  7677190.317446775,
  7682425.0,
  7683807.8125,
  7683960.109159015,
  7684184.7994470615,
  7684206.194942015,
  7686345.07821117,
  7687282.027107207,
  7690796.875,
  7690810.614015361,
  7691735.302900683,
  7692922.4904334545,
  7692935.517091191,
  7704117.347873407,
  7705287.5,
  7705876.371888888,
  7706700.0,
  7708684.375,
  7708724.897220039,
  7709552.185082504,
  7710551.5625,
  7711298.249613774,
  7712076.298431986,
  7712229.6875,
  7714553.485871672,
  7716235.752709214,
  7719565.511932308,
  7720832.781053785,
  7720887.5,
  7723117.1875,
  7725433.831847788,
  7735924.301922041,
  7747629.945353442,
  7748082.513033068,
  7750192.142803754,
  7751469.642833331,
  7751663.953054529,
  7754494.266689285,
  7754900.442947881,
  7755684.375,
  7755784.56515661,
  7755860.9375,
  7757521.875,
  7761043.75,
  7761525.0,
  7762597.610975244,
  7763178.125,
  7763250.856852085,
  7763937.710940844,
  7764038.509513045,
  7764114.408345873,
  7764550.0,
  7764590.688630665,
  7764998.4375,
  7766364.043550493,
  7767175.0,
  7767673.4375,
  7768175.148553596,
  7769952.683060873,
  7770311.449855896,
  7770827.02526446,
  7773127.937071002,
  7775690.625,
  7783923.269901901,
  7786251.211772154,
  7787733.52179854,
  7788352.344947424,
  7809994.5073372945,
  7821059.375,
  7822560.559423309,
  7822659.139663249,
  7825477.844137313,
  7825526.378763232,
  7825692.1875,
  7825704.079392292,
  7829315.625,
  7833160.406130971,
  7833259.375,
  7833314.871079315,
  7834111.786234456,
  7834416.500473576,
  7834965.364753085,
  7835912.519376425,
  7836031.917813412,
  7838401.5625,
  7863032.110494282,
  7921702.291830664,
  7945260.601712729,
  7947456.784712147,
  7955254.923686325,
  7957453.125,
  7968525.714854295,
  7975208.703757318,
  7978595.168187187,
  8015402.231421224,
  8015760.449504381,
  8016258.667728237,
  8016607.5604044385,
  8020241.9054473275,
  8022636.682392288,
  8023439.0625,
  8024208.105269995,
  8024979.1824363535,
  8025038.3967448,
  8025178.929093858,
  8034829.380865554,
  8040323.564832087,
  8040382.665801466,
  8041511.289273224,
  8045439.8724204935,
  8045841.311765584,
  8046823.34052544,
  8048256.93844812,
  8048334.327577891,
  8048555.779237429,
  8048890.5710173445,
  8049136.474058713,
  8049313.746090505,
  8049513.7644822635,
  8049647.805680556,
  8049743.568639341,
  8049929.533299217,
  8049941.820588816,
  8049960.9375,
  8049969.91470732,
  8050437.5,
  8050465.413352307,
  8050519.471912093,
  8050848.400838206,
  8051066.9303468885,
  8051665.456136617,
  8052010.173352713,
  8052112.4118886795,
  8052207.185729658,
  8052296.875,
  8052500.458025544,
  8052754.698346187,
  8053609.662566652,
  8056897.520079125,
  8057525.195532655,
  8057626.824388124,
  8058221.875,
  8058278.744517987,
  8062337.5,
  8064962.433896966,
  8071439.0625,
  8117371.352970065,
  8121045.683174125,
  8131168.75,
  8132333.0812669145,
  8133817.124324756,
  8137193.75,
  8140355.6252022395,
  8143071.975665295,
  8177413.946792574,
  8182925.0,
  8183262.266429158,
  8183282.801933412,
  8184017.1875,
  8262919.676514626,
  8319524.453658248,
  8320258.112596558,
  8587742.383229863,
  8587789.0625,
  8596391.256393256,
  8608374.578842437,
  8608943.83626232,
  8610839.122841472,
  8611495.241446324,
  8611543.606846513,
  8619893.678702816,
  8634428.929455204,
  8665534.375,
  8666536.712971518,
  8668197.883039132,
  8859658.38476883,
  8881511.71554166,
  8883381.561294029,
  8900606.18624703,
  8900914.153240364,
  8902285.112518042,
  8902397.981962793,
  8902476.125288535,
  8903878.070447821,
  8904249.635932965,
  8904428.693151668,
  8906205.826134952,
  8908557.003073901,
  8935698.369310148,
  8935764.0625,
  8936295.3125,
  8946430.985996012,
  8947640.816003395,
  8948671.678208986,
  8948756.252352752,
  8948870.647372551,
  8949218.535801215,
  8949322.570562083,
  8949475.429906635,
  8953466.860424375,
  8956999.85178064,
  8957852.23299286,
  8962763.333784113,
  8962821.720186707,
  8963335.9375,
  8992293.75,
  8997832.8125,
  8997859.365292648,
  9010415.625,
  9017079.6875,
  9066584.21633308,
  9067332.8125,
  9067412.688677158,
  9072793.062565196,
  9075806.25,
  9080371.875,
  9083409.375,
  9153968.440669978,
  9165490.625,
  9168406.25,
  9179781.288412115,
  9189612.5,
  9192460.363105677,
  9194167.1875,
  9194639.767033448,
  9202335.3049427,
  9249010.281356322,
  9249443.899958383,
  9250288.036185153,
  9251842.404319722,
  9252202.637876833,
  9252948.4375,
  9255107.906843672,
  9255286.727617305,
  9255894.164045552,
  9255990.625,
  9257420.48283648,
  9257495.3125,
  9260731.593960673,
  9266870.254935114,
  9269437.747614633,
  9270153.577840896,
  9271154.156907385,
  9271160.328719957,
  9271168.563589526,
  9274792.769897306,
  9280541.542912329,
  9281732.089715682,
  9329502.853522174,
  9338001.538835473,
  9338771.977189163,
  9340776.099574013,
  9341239.663460804,
  9342315.625,
  9409865.91240456,
  9429447.17349302,
  9444318.75,
  9481693.814526606,
  9508645.3125,
  9510992.89959578,
  9514419.547320418,
  9516304.6875,
  9519653.89717113,
  9521137.5,
  9527107.8125,
  9595102.958110666,
  9711871.360816404,
  9815831.592904603,
  9823874.118673626,
  9823880.209583577,
  9824207.524427542,
  9825361.428761894,
  9830106.798444437,
  9831970.894541385,
  9840715.625,
  9843538.27984846,
  9892337.337896347,
  9938704.17520631,
  9939038.498756967,
  9957040.625,
  9980320.431322752,
  9980928.993972246,
  9982367.740320219,
  9992629.70397702,
  9996173.253948873,
  10058179.216325019,
  10082157.8125,
  10099350.0,
  10149936.94229971,
  10150868.385338144,
  10155579.555075698,
  10159236.54191782,
  10164091.137727635,
  10164912.817033634,
  10166479.98664615,
  10166712.5,
  10167534.375,
  10168466.908327503,
  10169152.688641166,
  10173309.542516291,
  10174112.5,
  10175453.869309856,
  10176014.34157418,
  10176147.88502453,
  10176196.875,
  10177354.075021673,
  10177476.5625,
  10177632.8125,
  10178531.025056299,
  10178925.333659176,
  10179859.858418666,
  10180989.232086599,
  10181584.375,
  10181921.982898805,
  10181993.75,
  10182011.960467963,
  10183334.375,
  10184765.292858748,
  10185068.061137648,
  10185526.5625,
  10185625.438219737,
  10185627.08548888,
  10187477.90595186,
  10187887.5,
  10188219.054379314,
  10191548.4375,
  10197010.9375,
  10197971.567508955,
  10198071.158335626,
  10201736.259498581,
  10206708.190323729,
  10210881.25,
  10214397.337914497,
  10215226.742799269,
  10215227.496320343,
  10246248.4375,
  10261672.217736252,
  10264516.783159353,
  10269250.577283144,
  10285512.084199846,
  10286591.187746808,
  10287396.875,
  10287406.25,
  10290715.53321069,
  10292374.219148992,
  10345453.125,
  10366834.088538708,
  10408426.451024199,
  10415284.04161117,
  10415443.75,
  10415888.323784096,
  10417464.0625,
  10418431.353151234,
  10420387.354281852,
  10420738.669505117,
  10423329.789523406,
  10424124.473360265,
  10424710.9375,
  10424845.3125,
  10424987.710396757,
  10425013.556489943,
  10425446.875,
  10429316.855651218,
  10429819.336790143,
  10443600.345360486,
  10444951.5625,
  10446738.61328377,
  10449353.958790153,
  10450143.916567834,
  10455825.27547293,
  10498625.917974658,
  10516085.9375,
  10522646.593404826,
  10533317.558245003,
  10533732.62736665,
  10537687.5,
  10543342.1875,
  10543373.686645012,
  10554498.837182047,
  10558617.956379369,
  10563234.375,
  10564486.725717241,
  10565403.125,
  10565564.93853149,
  10565673.016238809,
  10565896.875,
  10565998.4375,
  10566718.41956189,
  10566728.578391833,
  10566806.25,
  10566965.625,
  10567043.392290866,
  10567181.25,
  10567260.216435133,
  10567451.5625,
  10567487.644836836,
  10567582.8125,
  10567643.394951439,
  10567842.039170349,
  10567868.75,
  10568107.741681658,
  10568240.494677877,
  10568347.67987977,
  10568351.660458645,
  10568357.618779495,
  10568368.677976139,
  10568385.9375,
  10568680.686633037,
  10568684.375,
  10568885.821263582,
  10569118.65168374,
  10569152.401719801,
  10569177.285681259,
  10569251.5625,
  10569361.887776734,
  10569459.375,
  10569520.3125,
  10569578.178846367,
  10569633.942680601,
  10569811.97602166,
  10569895.3125,
  10570029.6875,
  10570227.644355712,
  10570479.6875,
  10570507.453489214,
  10570573.056735102,
  10570868.75,
  10570881.315848306,
  10571049.428281816,
  10571300.0,
  10571733.477127165,
  10571848.850952918,
  10572145.3125,
  10572149.64550355,
  10572169.68631037,
  10572219.241805673,
  10572751.5625,
  10573114.849994272,
  10573329.136818837,
  10573368.040424915,
  10573884.375,
  10573885.991845291,
  10574179.028357456,
  10575708.493846755,
  10575820.53458742,
  10575939.11864759,
  10575956.25,
  10576259.27274382,
  10576654.6875,
  10577139.412696933,
  10577323.4375,
  10577675.0,
  10577730.830906747,
  10580065.869354483,
  10582627.71610995,
  10582940.391953839,
  10583024.784083473,
  10583235.940249423,
  10583451.5625,
  10583580.134203324,
  10583585.9375,
  10583620.189527338,
  10583891.566072026,
  10583893.774212722,
  10584210.865249392,
  10584312.750949088,
  10584386.848524414,
  10584391.203458857,
  10584410.854561605,
  10584426.5625,
  10584464.0625,
  10584664.341340942,
  10584858.169555966,
  10584931.922686916,
  10584995.284722965,
  10585044.654981168,
  10585318.405018348,
  10585360.254112825,
  10585482.8125,
  10585548.867458332,
  10585781.25,
  10585783.46990954,
  10585918.75,
  10586100.584047938,
  10586141.974293912,
  10586145.488775168,
  10586187.5,
  10586230.686924592,
  10586390.625,
  10586408.631720155,
  10586418.75,
  10586517.1875,
  10586555.16350034,
  10586812.109829774,
  10586925.553527204,
  10587115.625,
  10587312.87950314,
  10587501.5625,
  10587567.177092817,
  10587593.492776426,
  10587665.082314996,
  10587703.02778036,
  10587718.75,
  10587928.510295745,
  10588153.125,
  10588181.28257486,
  10588193.708831562,
  10588214.884376239,
  10588704.6875,
  10588765.625,
  10589155.912788745,
  10589476.261585321,
  10589704.589205928,
  10589765.625,
  10589773.4375,
  10590460.9375,
  10590932.290547805,
  10591484.375,
  10592675.599312136,
  10593053.992193552,
  10593148.037111022,
  10594687.197271833,
  10594689.0625,
  10595043.75,
  10596071.113808136,
  10596081.456567012,
  10599826.5625,
  10613496.742654076,
  10613535.9375,
  10660493.61125328,
  11266312.259867236,
  11270329.6875,
  11283058.054572864,
  11312946.205440858,
  11464919.360253304,
  11504846.195103101,
  11842039.1701484,
  11853501.70948552,
  12222704.818192063,
  12421846.875,
  13021132.129239254,
  13021236.963152034,
  14088962.023301644,
  15352404.085066361,
  16915983.181643557,
  18774067.1875,
  21136857.264631905,
  21150732.392295297,
  21150897.419354163,
  21151325.0,
  21153409.375,
  21154900.0,
  21155055.155825462,
  21156345.110535234,
  21156879.909279555,
  21158929.6821773,
  21159081.25,
  21161334.33923666,
  21161639.829834018,
  21161872.86379478,
  21171122.908343907,
  21171507.8125,
  21185470.3125],
 [44.1151146399842,
  19.06426387887476,
  72.9903032619747,
  38.13472018705268,
  6.325439695638338,
  24.822227514762076,
  11.746959293914506,
  8.462785199678626,
  10.954056047980188,
  12.775120256909076,
  10.179018807807488,
  80.61816726765272,
  30.76843308656368,
  35.224487691026255,
  59.15985230392127,
  14.4969808328599,
  6.11937455439975,
  28.21834154736074,
  58.469715920556624,
  7.271594547236137,
  39.89527397712712,
  77.67457694140487,
  83.48073164201625,
  12.325121847512246,
  37.14355983439159,
  60.76458705799514,
  57.60419889689573,
  58.90327392228315,
  55.142436694849835,
  5.448469472028692,
  5.266140685345057,
  26.699204358902662,
  60.43746248368338,
  12.284410691492242,
  54.28663991460195,
  19.65862071861211,
  14.313882741610369,
  9.802818611938054,
  12.5036899595841,
  56.865680112598156,
  12.956609693402438,
  72.57305452599724,
  18.08948475784833,
  14.45857897813693,
  13.058492048759996,
  56.59026424879305,
  55.5293496600556,
  26.83241979143611,
  9.660274476863322,
  5.277423733475989,
  19.921068089425823,
  13.13873509542886,
  47.81287930275692,
  27.574392241768237,
  9.046857268249944,
  5.6822278719798485,
  8.926881392463072,
  44.51478965713822,
  34.64802089011277,
  5.961892628678742,
  16.354084309577956,
  52.67755841132305,
  15.466787815070466,
  9.107993885711476,
  11.85406913600058,
  6.828359929736517,
  10.574297837662817,
  15.43331884419696,
  22.410809872894255,
  28.08726405099809,
  18.331705734777255,
  61.69834655165389,
  95.62651882812845,
  15.872847209055095,
  67.69293541037939,
  44.90747192526904,
  21.97033749124155,
  6.417337334516956,
  5.5673057473178735,
  53.182527603473076,
  52.789231939786504,
  5.339086024129061,
  18.196465446112608,
  8.129334092487895,
  59.35673075711627,
  42.681243451635034,
  88.56324384782026,
  104.62034746196505,
  63.05397070058917,
  48.0120416650076,
  12.04521380241369,
  28.705100479201604,
  12.601697867579386,
  6.072379288653127,
  65.81251228832417,
  43.855912782894464,
  11.542550740059664,
  47.16632222383809,
  19.154453525451583,
  55.40464740080633,
  24.66976290069403,
  106.93358363554164,
  11.680191733130869,
  103.20938263705304,
  14.409855284907763,
  6.394679956988659,
  49.697252317956526,
  36.87084176635848,
  25.082256671901476,
  74.84550407404369,
  12.968642171618674,
  20.14591987843851,
  5.640927607687405,
  38.44459780334823,
  9.342148523193076,
  13.1292821919657,
  5.997729949451351,
  21.627660035575317,
  39.39605867880935,
  12.400504983076265,
  31.825092144975738,
  82.28018200319113,
  12.432498111721994,
  13.077675289573866,
  45.401253421468454,
  43.34713211722807,
  65.99974675542694,
  25.24529936553993,
  21.57063169696749,
  41.82884419973743,
  37.86790300012656,
  8.508581803636059,
  39.97695041433258,
  16.86400992609002,
  45.05263228286865,
  14.284770459610836,
  34.9029385521692,
  29.498099299066983,
  44.6242773481123,
  9.021528827598978,
  33.47237665238689,
  24.043689620182455,
  56.52747540026464,
  32.37333376757789,
  64.53368214734593,
  69.62323855926952,
  19.3573075450932,
  69.18362975044886,
  6.189623192982941,
  17.127524836659706,
  94.71629865961143,
  17.73684762004018,
  128.74327013315232,
  5.758878386471944,
  12.118823143131959,
  8.84862856773068,
  5.649018814768355,
  11.445569218297525,
  6.10478367403138,
  37.13639519310876,
  47.819700209322924,
  19.734631054615683,
  37.30336617237097,
  23.47551075640914,
  30.07025744852214,
  16.84755776168545,
  15.08784426023395,
  12.716919264076385,
  8.971344247998388,
  108.70073457902065,
  70.31945360752536,
  17.823422729232263,
  27.922835638237142,
  38.30954713210878,
  45.710242536650156,
  13.28270205813922,
  110.16489039160074,
  38.80967442941125,
  62.45460963814571,
  32.25537420516064,
  5.655993084918314,
  7.764322597989998,
  14.213397081976634,
  38.76950877128941,
  19.05133554295636,
  78.7293859284863,
  5.911721051721169,
  20.339693383724885,
  16.25350287319941,
  16.545055094874876,
  16.598355608792932,
  46.509058977045704,
  9.175268653739458,
  87.88679300656369,
  75.78015652481228,
  26.645590883972027,
  5.564079130979073,
  75.26612782309063,
  70.80270410025568,
  7.053724254883401,
  39.548391287150906,
  6.313894561559217,
  25.343081561235234,
  8.227676788332154,
  27.917361251380587,
  44.031519762669234,
  47.62848283814515,
  8.018818709501337,
  77.50765631026522,
  16.02471019408454,
  87.5737807581747,
  14.773841686147545,
  19.952559522320428,
  65.96123218190318,
  12.044329315181297,
  11.898659173965536,
  29.5375383703038,
  7.1723515961332005,
  85.26475400441294,
  102.01775663013956,
  42.78794565295557,
  34.76736928452512,
  12.512683336939471,
  110.02825299799426,
  19.660748849851885,
  29.870095277623705,
  5.152679894997583,
  16.446819854709723,
  40.744659035337634,
  11.228643843107745,
  49.965925034787745,
  73.34811212073528,
  65.04944831015943,
  50.3945203782646,
  26.554241111965332,
  22.09252417201001,
  12.035614029447801,
  8.547456559337043,
  14.270531638196326,
  29.294444274851116,
  25.60338234327652,
  48.27599220935607,
  55.84730174540409,
  70.1200800220691,
  8.732255198758045,
  34.2173478441145,
  13.672847124999334,
  59.998747703778605,
  73.7752923924916,
  25.358856134644583,
  74.90341900558244,
  10.293895288654406,
  93.72194249286376,
  75.81884456317172,
  73.40699059816055,
  32.177150977121016,
  12.776389041989697,
  19.810395727676376,
  10.74391241859486,
  27.713833555461225,
  20.05444870105609,
  32.1465109091896,
  7.947235684971142,
  35.77949498082458,
  37.37167673741854,
  27.620099492935076,
  19.945952482713597,
  137.97105921619627,
  10.019067323920998,
  12.880306976007637,
  62.547279902002046,
  7.309140683844345,
  77.80728321381923,
  12.792094782112349,
  69.68048274756882,
  5.237319199379071,
  10.052132103849646,
  56.6528496199913,
  5.352054937020157,
  32.43299770226985,
  5.226886544443705,
  28.497173999539996,
  10.680012884976959,
  15.965543716271927,
  38.411749146696785,
  16.4763527377916,
  16.84595594224027,
  53.745381150535614,
  78.04525656715913,
  53.535962147417116,
  24.11018810374957,
  21.78980303610099,
  11.044363929260559,
  36.861906921835,
  17.267103874493166,
  9.091514360962233,
  7.3896592883333945,
  25.3945043178869,
  25.81073607046493,
  22.932084522689678,
  111.24845415770794,
  70.5817401917607,
  31.95795808170676,
  24.948974607958476,
  6.004840238454423,
  7.111847620016791,
  8.155116062592676,
  25.584005990752164,
  46.05357217033031,
  9.257994997089728,
  5.667996941309796,
  74.81569555068364,
  22.74395239701665,
  11.144702517986893,
  16.163918129957743,
  58.445715835220625,
  70.70819859710596,
  37.86236301120175,
  28.637458810901727,
  12.25708673720708,
  34.925675696221234,
  78.35091797027128,
  75.8754487392097,
  6.408258008946374,
  12.953418872988227,
  52.60236289268757,
  9.569493099430568,
  5.33976154816605,
  31.146579840221357,
  20.33239872212456,
  26.468911006709785,
  71.89694592944358,
  44.65323010397874,
  14.552860842324495,
  11.111230616288038,
  79.6473702368225,
  20.23597078053509,
  7.416905853498759,
  6.3058728422776325,
  21.744129771305836,
  69.46254431765593,
  31.070871907119106,
  14.540261200490216,
  93.1089906854249,
  5.944677707199027,
  47.70924175827307,
  7.084322631769302,
  5.3394845099884725,
  63.58578075136934,
  88.63714682439617,
  19.039342417512334,
  7.414296949653421,
  23.44119228206218,
  10.921669192712823,
  16.277486191274637,
  11.029497843892338,
  5.94734689058128,
  12.516180182975104,
  20.015444147434213,
  7.011076823634419,
  31.817788417424886,
  10.60264287948333,
  17.28167714313644,
  6.422637911135516,
  120.12891584833949,
  11.022768610164345,
  46.08537567168395,
  51.88826173624381,
  24.22169463624742,
  26.73061451819018,
  92.18299088748196,
  14.95048459754247,
  38.327774785568856,
  46.61623283478523,
  40.42659496383164,
  5.356872034598011,
  100.77168584833755,
  9.209142916726053,
  18.308480791907638,
  19.479335078374206,
  5.696975610512134,
  36.41198809006641,
  52.039281505551386,
  65.6319387126819,
  9.323208525588003,
  72.47052907164043,
  50.936769141401044,
  41.113564597624965,
  50.68197049373533,
  5.102187912127682,
  18.608548102569607,
  68.73205538791473,
  48.80228820514137,
  70.69101785032288,
  54.61711063185947,
  37.663163279959825,
  75.57348890239108,
  8.306995287967535,
  43.63043132243486,
  30.936921061156582,
  118.11881886514102,
  6.360512981174484,
  12.97709501116696,
  20.065126038915896,
  6.7938184861573925,
  17.473040884504,
  52.702457743103146,
  31.980294675116177,
  67.87911720947356,
  15.313724198302006,
  5.429563640523038,
  16.79555780845051,
  85.10854728562053,
  15.633090702010843,
  60.15107090139588,
  11.711559948534903,
  45.44521950597776,
  51.97451049250542,
  27.17248980605158,
  55.94857937069234,
  45.875253225879774,
  8.190400442147737,
  12.178940989205541,
  7.972786132607971,
  7.3611880950442625,
  32.61733806122776,
  86.18945987039602,
  6.5032304853187375,
  80.76789131992854,
  21.590644063066907,
  49.73562227932904,
  16.65894442596359,
  53.996355074428614,
  38.270144423115774,
  193.64953557646763,
  48.35176764591674,
  72.09628287322178,
  12.654069756665395,
  38.8339941460539,
  12.595519237052462,
  122.17652320991567,
  32.78708495522601,
  6.246340876534137,
  21.35276860440632,
  103.65006610807802,
  18.37950320939907,
  65.30875245199718,
  31.777682584271428,
  11.9934698260187,
  77.92956892869215,
  36.49721293314736,
  71.239265142184,
  122.0610606632733,
  35.907863728312236,
  6.204573074331645,
  5.1524370709084275,
  55.804340067416746,
  20.779583193009863,
  19.75781849101845,
  17.571055407258257,
  52.55299288463493,
  50.14998879565123,
  35.65022850058403,
  85.41052275778333,
  56.2555650788846,
  57.23849621100226,
  9.040988273203602,
  5.393435592483144,
  17.884595057155213,
  57.504636863031024,
  69.78241276972112,
  11.629865209189063,
  60.89575977819283,
  30.90328360084,
  29.8121935448994,
  8.987816759050844,
  15.62090695596031,
  49.25848111915589,
  46.702426614008445,
  16.12379318123065,
  117.69004611416153,
  37.48535648654719,
  10.288642889780611,
  79.41222840551767,
  33.23594943738304,
  30.341678915904232,
  16.3865490921623,
  37.26010689837779,
  13.09119309603918,
  52.93914847664452,
  15.369360924522296,
  47.46616114642142,
  13.219065535564372,
  94.96129393064626,
  9.065508855893341,
  107.57836860958712,
  86.15077401788794,
  9.73163239302457,
  17.396314484373416,
  14.167533245697808,
  32.058409728159546,
  6.825834002638735,
  51.239449661593525,
  58.93740336626033,
  23.655523237322512,
  35.64049152404362,
  16.32943919920143,
  149.95813784939372,
  75.92476748778904,
  14.806688492056985,
  11.42172709769792,
  12.986917651710812,
  58.892394085191164,
  40.04194519137185,
  22.50355454347949,
  17.71787837931663,
  7.01627802521345,
  7.820297110945039,
  125.06472974834055,
  47.76856722361478,
  7.57037303484295,
  19.08473940808929,
  41.835368283748984,
  7.502814439472981,
  9.238889874409926,
  23.33341351072511,
  29.93371572184674,
  12.716583438589957,
  52.81800808660188,
  97.58754923380447,
  104.53354352277216,
  13.628983579329898,
  15.345176395034061,
  13.253673447037574,
  49.44141229447902,
  86.91245433298742,
  81.4023266407402,
  24.244521008207688,
  54.89592184591685,
  13.623843696980469,
  17.775820838582288,
  26.735643861058588,
  31.893504129394383,
  62.15543177348213,
  12.740465708216751,
  90.58736324751223,
  13.36420624141576,
  63.514639157202254,
  29.343776749341284,
  9.104132199686298,
  24.864152586811628,
  110.25589949515893,
  80.38591094106349,
  42.76975445077802,
  26.92601703749842,
  5.139137018060453,
  18.991205300705037,
  30.743551390884154,
  10.100506165017265,
  188.28953587846502,
  26.5083642002,
  7.457664106588232,
  85.2037456826129,
  60.294774072359175,
  37.186989379124206,
  7.2459628458718335,
  39.30215180743267,
  64.30323022022141,
  23.040552951353664,
  69.11603749455873,
  54.62361005109108,
  76.65400538607057,
  82.76549022260227,
  110.53504336507139,
  130.87130999124744,
  6.072808879926113,
  58.111851209136226,
  105.35752406971719,
  34.38036317554156,
  28.174263923786683,
  32.93411117660727,
  32.99055313746904,
  82.34374965142545,
  68.38648816663668,
  13.155785748994665,
  91.90680382008566,
  31.983977102081052,
  12.576389497423712,
  91.49274019929285,
  8.710000511982177,
  124.78674027659044,
  21.191379185396155,
  58.310930454959674,
  48.11106634708665,
  78.71454323431442,
  8.990780491054972,
  21.257192502627564,
  44.14320460383532,
  68.64433277544835,
  14.795440048944869,
  52.42760322465816,
  5.684284310337724,
  13.664912508120777,
  80.15294075129698,
  59.25054773553046,
  20.806311849834053,
  84.1034830302795,
  29.654837220431865,
  97.93362790216196,
  77.60043279710577,
  23.078720834992595,
  20.320338310963102,
  8.564982134652046,
  74.02373323676741,
  7.915265259407697,
  8.905475779391224,
  89.78927083052886,
  33.58490252823504,
  15.15277463915615,
  7.966275210505147,
  9.982631873396429,
  48.51825091615752,
  24.3351324715449,
  37.26790851629909,
  14.565572918921177,
  81.06153852630564,
  45.169041940064815,
  66.2566467434483,
  42.746235222276084,
  121.97616513811549,
  8.420166007439002,
  6.4248491233608345,
  36.864613302282265,
  140.99979581251912,
  97.79007539472119,
  120.14393530791497,
  16.9697638552666,
  71.91385799543811,
  5.241716525991229,
  7.983194593430584,
  9.351960656227808,
  61.69359033606149,
  27.516682857566728,
  51.06542241085249,
  66.79630259813428,
  14.634537460080129,
  90.59091303014901,
  61.085293292397424,
  11.593442606607278,
  7.218584255612209,
  7.310268912187566,
  11.084565318484033,
  20.0536050738007,
  30.29319862748777,
  54.913324375700626,
  60.3459176569474,
  70.39077208628738,
  36.229832084454166,
  29.14664949548409,
  73.83968586297199,
  25.61621942376874,
  16.977445889745475,
  15.658833202617416,
  35.54911488093509,
  23.770979608342763,
  35.20370765171865,
  53.877788441314955,
  70.65449756901506,
  46.24701145989481,
  7.252463696345999,
  7.467963201015802,
  24.255745544995037,
  54.7784164460617,
  119.0141470367416,
  70.16071427657803,
  28.32061750196084,
  12.370195205868335,
  75.99137627421038,
  7.096127943969323,
  70.97090832941112,
  38.81718835026143,
  63.86188918029771,
  87.56821433702686,
  27.735560469253016,
  12.712045199423539,
  79.77510687315201,
  35.228424982198234,
  34.70263446960864,
  71.83231182535239,
  55.88773199103073,
  36.908583024142715,
  30.894621351864586,
  81.93839162962121,
  36.04662453303813,
  50.3789414816418,
  39.620746577144274,
  5.881058026258764,
  31.25339257157878,
  74.79105159945803,
  34.99223608786868,
  67.50951682876378,
  19.231094599960265,
  15.791986493508077,
  57.56086448587226,
  26.157807521580434,
  100.69714851916748,
  79.72215670357868,
  7.5186349885285635,
  61.68177699135296,
  11.528407073010875,
  56.22668864926588,
  5.051980027559639,
  5.882305076554136,
  6.280895447041716,
  11.3654990365045,
  58.890535397130265,
  64.39763481596512,
  12.850678556387784,
  5.032611640767161,
  21.728739094415914,
  40.9365438830743,
  9.664339593854667,
  47.954509770718566,
  21.03970765732501,
  22.327828817994206,
  25.229019999770614,
  54.23928613320853,
  64.52049644226265,
  186.17141128534433,
  45.186552877895515,
  29.806307565536624,
  8.40744664974339,
  58.31608778728173,
  72.80306627412551,
  30.654815742712234,
  16.34176916151307,
  13.801535634857345,
  26.63197255242998,
  35.26472329384794,
  56.45371909211843])

Energy spectrum

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

([2128322.8820732716,
  2673853.9038674626,
  2786856.25,
  2789778.125,
  2805233.259980844,
  2807010.716995365,
  2808506.6130725285,
  2812117.506482652,
  2874049.8859884213,
  2889719.5982088023,
  2889974.2945079044,
  2890028.125,
  2891280.9293882987,
  2892351.5625,
  2893143.3147037867,
  2905490.5032702447,
  2905515.9204744357,
  2906018.734324037,
  2948343.75,
  3034213.041952173,
  3038439.0625,
  3141742.1875,
  3424040.625,
  6785328.776244561,
  6785493.381396405,
  6791056.25,
  6795218.742225575,
  6803876.5625,
  6861424.907240535,
  6866769.9989933055,
  6866830.420718658,
  6879946.775391034,
  6940407.8125,
  6953758.444825136,
  6977915.625,
  6978375.856584645,
  6982520.137476653,
  6982803.634441398,
  6987272.607001323,
  6988256.8744002385,
  6988266.546827361,
  6988965.625,
  6990647.289649813,
  6992557.203741236,
  6993033.774311789,
  7020286.637365565,
  7110217.1875,
  7114679.595704689,
  7115817.552414775,
  7116227.676238994,
  7117364.95151519,
  7118310.367652139,
  7118339.0625,
  7122178.9844065495,
  7123494.955652764,
  7125337.107539448,
  7125688.541932925,
  7129987.5,
  7130548.4375,
  7146026.817144993,
  7157582.864702709,
  7169130.601346257,
  7172898.599309346,
  7173441.9109135885,
  7174609.7948041605,
  7176262.816393783,
  7188563.677035553,
  7192779.81081335,
  7193623.024825048,
  7193984.352505192,
  7199272.54899749,
  7199509.375,
  7210189.0625,
  7219604.912121246,
  7220298.4375,
  7220809.009714929,
  7255282.293938559,
  7269896.415108262,
  7292217.873084736,
  7301620.525656571,
  7303248.4375,
  7306081.167012519,
  7308306.903424871,
  7308351.252823675,
  7310126.5625,
  7310192.1875,
  7312481.25,
  7318070.146381164,
  7318290.625,
  7320417.1875,
  7331122.39212852,
  7351126.5625,
  7352783.587341864,
  7352798.284965201,
  7357039.0625,
  7366504.6875,
  7368967.886011667,
  7369167.9398481585,
  7369285.676346161,
  7369331.25,
  7369440.2057869695,
  7369593.75,
  7370225.326871223,
  7370403.450339849,
  7370856.366219008,
  7371115.00902001,
  7371560.9375,
  7372068.883017307,
  7375343.5475491695,
  7375645.3125,
  7376091.551598742,
  7376977.478443644,
  7377316.100927785,
  7377719.397872175,
  7379546.887680742,
  7380329.11404829,
  7380764.708721365,
  7380780.7633840125,
  7380978.125,
  7381117.873272525,
  7381319.525419512,
  7381475.772936901,
  7381608.512732085,
  7381697.844236114,
  7382160.9375,
  7382195.188103089,
  7382382.979949364,
  7383914.557640953,
  7384037.262309488,
  7384590.995170809,
  7385517.137191394,
  7386523.834975215,
  7386632.8125,
  7407451.204093161,
  7411083.48540908,
  7416076.917212375,
  7420263.841231703,
  7421003.125,
  7423408.618539594,
  7424423.0467876475,
  7424773.4375,
  7443801.285092924,
  7473906.030115981,
  7498054.6875,
  7498323.4375,
  7514310.9375,
  7514313.62134014,
  7515864.0625,
  7517616.54779265,
  7517636.413010648,
  7517884.767175848,
  7518367.767904604,
  7518534.561223161,
  7518955.476446981,
  7519572.428332865,
  7520449.367957458,
  7520524.0528527545,
  7523556.444266042,
  7524612.658927334,
  7529778.125,
  7531141.88947555,
  7531305.553590942,
  7532726.177962356,
  7532863.096114185,
  7532898.4375,
  7533681.611589373,
  7538795.601631367,
  7540046.136425193,
  7553094.853595726,
  7553971.859383369,
  7566477.775048121,
  7668799.085138062,
  7671643.07313268,
  7676275.0,
  7676429.536808817,
  7676441.053957746,
  7676862.5,
  7677190.317446775,
  7682425.0,
  7683807.8125,
  7683960.109159015,
  7684184.7994470615,
  7684206.194942015,
  7686345.07821117,
  7687282.027107207,
  7690796.875,
  7690810.614015361,
  7691735.302900683,
  7692922.4904334545,
  7692935.517091191,
  7704117.347873407,
  7705287.5,
  7705876.371888888,
  7706700.0,
  7708684.375,
  7708724.897220039,
  7709552.185082504,
  7710551.5625,
  7711298.249613774,
  7712076.298431986,
  7712229.6875,
  7714553.485871672,
  7716235.752709214,
  7719565.511932308,
  7720832.781053785,
  7720887.5,
  7723117.1875,
  7725433.831847788,
  7735924.301922041,
  7747629.945353442,
  7748082.513033068,
  7750192.142803754,
  7751469.642833331,
  7751663.953054529,
  7754494.266689285,
  7754900.442947881,
  7755684.375,
  7755784.56515661,
  7755860.9375,
  7757521.875,
  7761043.75,
  7761525.0,
  7762597.610975244,
  7763178.125,
  7763250.856852085,
  7763937.710940844,
  7764038.509513045,
  7764114.408345873,
  7764550.0,
  7764590.688630665,
  7764998.4375,
  7766364.043550493,
  7767175.0,
  7767673.4375,
  7768175.148553596,
  7769952.683060873,
  7770311.449855896,
  7770827.02526446,
  7773127.937071002,
  7775690.625,
  7783923.269901901,
  7786251.211772154,
  7787733.52179854,
  7788352.344947424,
  7809994.5073372945,
  7821059.375,
  7822560.559423309,
  7822659.139663249,
  7825477.844137313,
  7825526.378763232,
  7825692.1875,
  7825704.079392292,
  7829315.625,
  7833160.406130971,
  7833259.375,
  7833314.871079315,
  7834111.786234456,
  7834416.500473576,
  7834965.364753085,
  7835912.519376425,
  7836031.917813412,
  7838401.5625,
  7863032.110494282,
  7921702.291830664,
  7945260.601712729,
  7947456.784712147,
  7955254.923686325,
  7957453.125,
  7968525.714854295,
  7975208.703757318,
  7978595.168187187,
  8015402.231421224,
  8015760.449504381,
  8016258.667728237,
  8016607.5604044385,
  8020241.9054473275,
  8022636.682392288,
  8023439.0625,
  8024208.105269995,
  8024979.1824363535,
  8025038.3967448,
  8025178.929093858,
  8034829.380865554,
  8040323.564832087,
  8040382.665801466,
  8041511.289273224,
  8045439.8724204935,
  8045841.311765584,
  8046823.34052544,
  8048256.93844812,
  8048334.327577891,
  8048555.779237429,
  8048890.5710173445,
  8049136.474058713,
  8049313.746090505,
  8049513.7644822635,
  8049647.805680556,
  8049743.568639341,
  8049929.533299217,
  8049941.820588816,
  8049960.9375,
  8049969.91470732,
  8050437.5,
  8050465.413352307,
  8050519.471912093,
  8050848.400838206,
  8051066.9303468885,
  8051665.456136617,
  8052010.173352713,
  8052112.4118886795,
  8052207.185729658,
  8052296.875,
  8052500.458025544,
  8052754.698346187,
  8053609.662566652,
  8056897.520079125,
  8057525.195532655,
  8057626.824388124,
  8058221.875,
  8058278.744517987,
  8062337.5,
  8064962.433896966,
  8071439.0625,
  8117371.352970065,
  8121045.683174125,
  8131168.75,
  8132333.0812669145,
  8133817.124324756,
  8137193.75,
  8140355.6252022395,
  8143071.975665295,
  8177413.946792574,
  8182925.0,
  8183262.266429158,
  8183282.801933412,
  8184017.1875,
  8262919.676514626,
  8319524.453658248,
  8320258.112596558,
  8587742.383229863,
  8587789.0625,
  8596391.256393256,
  8608374.578842437,
  8608943.83626232,
  8610839.122841472,
  8611495.241446324,
  8611543.606846513,
  8619893.678702816,
  8634428.929455204,
  8665534.375,
  8666536.712971518,
  8668197.883039132,
  8859658.38476883,
  8881511.71554166,
  8883381.561294029,
  8900606.18624703,
  8900914.153240364,
  8902285.112518042,
  8902397.981962793,
  8902476.125288535,
  8903878.070447821,
  8904249.635932965,
  8904428.693151668,
  8906205.826134952,
  8908557.003073901,
  8935698.369310148,
  8935764.0625,
  8936295.3125,
  8946430.985996012,
  8947640.816003395,
  8948671.678208986,
  8948756.252352752,
  8948870.647372551,
  8949218.535801215,
  8949322.570562083,
  8949475.429906635,
  8953466.860424375,
  8956999.85178064,
  8957852.23299286,
  8962763.333784113,
  8962821.720186707,
  8963335.9375,
  8992293.75,
  8997832.8125,
  8997859.365292648,
  9010415.625,
  9017079.6875,
  9066584.21633308,
  9067332.8125,
  9067412.688677158,
  9072793.062565196,
  9075806.25,
  9080371.875,
  9083409.375,
  9153968.440669978,
  9165490.625,
  9168406.25,
  9179781.288412115,
  9189612.5,
  9192460.363105677,
  9194167.1875,
  9194639.767033448,
  9202335.3049427,
  9249010.281356322,
  9249443.899958383,
  9250288.036185153,
  9251842.404319722,
  9252202.637876833,
  9252948.4375,
  9255107.906843672,
  9255286.727617305,
  9255894.164045552,
  9255990.625,
  9257420.48283648,
  9257495.3125,
  9260731.593960673,
  9266870.254935114,
  9269437.747614633,
  9270153.577840896,
  9271154.156907385,
  9271160.328719957,
  9271168.563589526,
  9274792.769897306,
  9280541.542912329,
  9281732.089715682,
  9329502.853522174,
  9338001.538835473,
  9338771.977189163,
  9340776.099574013,
  9341239.663460804,
  9342315.625,
  9409865.91240456,
  9429447.17349302,
  9444318.75,
  9481693.814526606,
  9508645.3125,
  9510992.89959578,
  9514419.547320418,
  9516304.6875,
  9519653.89717113,
  9521137.5,
  9527107.8125,
  9595102.958110666,
  9711871.360816404,
  9815831.592904603,
  9823874.118673626,
  9823880.209583577,
  9824207.524427542,
  9825361.428761894,
  9830106.798444437,
  9831970.894541385,
  9840715.625,
  9843538.27984846,
  9892337.337896347,
  9938704.17520631,
  9939038.498756967,
  9957040.625,
  9980320.431322752,
  9980928.993972246,
  9982367.740320219,
  9992629.70397702,
  9996173.253948873,
  10058179.216325019,
  10082157.8125,
  10099350.0,
  10149936.94229971,
  10150868.385338144,
  10155579.555075698,
  10159236.54191782,
  10164091.137727635,
  10164912.817033634,
  10166479.98664615,
  10166712.5,
  10167534.375,
  10168466.908327503,
  10169152.688641166,
  10173309.542516291,
  10174112.5,
  10175453.869309856,
  10176014.34157418,
  10176147.88502453,
  10176196.875,
  10177354.075021673,
  10177476.5625,
  10177632.8125,
  10178531.025056299,
  10178925.333659176,
  10179859.858418666,
  10180989.232086599,
  10181584.375,
  10181921.982898805,
  10181993.75,
  10182011.960467963,
  10183334.375,
  10184765.292858748,
  10185068.061137648,
  10185526.5625,
  10185625.438219737,
  10185627.08548888,
  10187477.90595186,
  10187887.5,
  10188219.054379314,
  10191548.4375,
  10197010.9375,
  10197971.567508955,
  10198071.158335626,
  10201736.259498581,
  10206708.190323729,
  10210881.25,
  10214397.337914497,
  10215226.742799269,
  10215227.496320343,
  10246248.4375,
  10261672.217736252,
  10264516.783159353,
  10269250.577283144,
  10285512.084199846,
  10286591.187746808,
  10287396.875,
  10287406.25,
  10290715.53321069,
  10292374.219148992,
  10345453.125,
  10366834.088538708,
  10408426.451024199,
  10415284.04161117,
  10415443.75,
  10415888.323784096,
  10417464.0625,
  10418431.353151234,
  10420387.354281852,
  10420738.669505117,
  10423329.789523406,
  10424124.473360265,
  10424710.9375,
  10424845.3125,
  10424987.710396757,
  10425013.556489943,
  10425446.875,
  10429316.855651218,
  10429819.336790143,
  10443600.345360486,
  10444951.5625,
  10446738.61328377,
  10449353.958790153,
  10450143.916567834,
  10455825.27547293,
  10498625.917974658,
  10516085.9375,
  10522646.593404826,
  10533317.558245003,
  10533732.62736665,
  10537687.5,
  10543342.1875,
  10543373.686645012,
  10554498.837182047,
  10558617.956379369,
  10563234.375,
  10564486.725717241,
  10565403.125,
  10565564.93853149,
  10565673.016238809,
  10565896.875,
  10565998.4375,
  10566718.41956189,
  10566728.578391833,
  10566806.25,
  10566965.625,
  10567043.392290866,
  10567181.25,
  10567260.216435133,
  10567451.5625,
  10567487.644836836,
  10567582.8125,
  10567643.394951439,
  10567842.039170349,
  10567868.75,
  10568107.741681658,
  10568240.494677877,
  10568347.67987977,
  10568351.660458645,
  10568357.618779495,
  10568368.677976139,
  10568385.9375,
  10568680.686633037,
  10568684.375,
  10568885.821263582,
  10569118.65168374,
  10569152.401719801,
  10569177.285681259,
  10569251.5625,
  10569361.887776734,
  10569459.375,
  10569520.3125,
  10569578.178846367,
  10569633.942680601,
  10569811.97602166,
  10569895.3125,
  10570029.6875,
  10570227.644355712,
  10570479.6875,
  10570507.453489214,
  10570573.056735102,
  10570868.75,
  10570881.315848306,
  10571049.428281816,
  10571300.0,
  10571733.477127165,
  10571848.850952918,
  10572145.3125,
  10572149.64550355,
  10572169.68631037,
  10572219.241805673,
  10572751.5625,
  10573114.849994272,
  10573329.136818837,
  10573368.040424915,
  10573884.375,
  10573885.991845291,
  10574179.028357456,
  10575708.493846755,
  10575820.53458742,
  10575939.11864759,
  10575956.25,
  10576259.27274382,
  10576654.6875,
  10577139.412696933,
  10577323.4375,
  10577675.0,
  10577730.830906747,
  10580065.869354483,
  10582627.71610995,
  10582940.391953839,
  10583024.784083473,
  10583235.940249423,
  10583451.5625,
  10583580.134203324,
  10583585.9375,
  10583620.189527338,
  10583891.566072026,
  10583893.774212722,
  10584210.865249392,
  10584312.750949088,
  10584386.848524414,
  10584391.203458857,
  10584410.854561605,
  10584426.5625,
  10584464.0625,
  10584664.341340942,
  10584858.169555966,
  10584931.922686916,
  10584995.284722965,
  10585044.654981168,
  10585318.405018348,
  10585360.254112825,
  10585482.8125,
  10585548.867458332,
  10585781.25,
  10585783.46990954,
  10585918.75,
  10586100.584047938,
  10586141.974293912,
  10586145.488775168,
  10586187.5,
  10586230.686924592,
  10586390.625,
  10586408.631720155,
  10586418.75,
  10586517.1875,
  10586555.16350034,
  10586812.109829774,
  10586925.553527204,
  10587115.625,
  10587312.87950314,
  10587501.5625,
  10587567.177092817,
  10587593.492776426,
  10587665.082314996,
  10587703.02778036,
  10587718.75,
  10587928.510295745,
  10588153.125,
  10588181.28257486,
  10588193.708831562,
  10588214.884376239,
  10588704.6875,
  10588765.625,
  10589155.912788745,
  10589476.261585321,
  10589704.589205928,
  10589765.625,
  10589773.4375,
  10590460.9375,
  10590932.290547805,
  10591484.375,
  10592675.599312136,
  10593053.992193552,
  10593148.037111022,
  10594687.197271833,
  10594689.0625,
  10595043.75,
  10596071.113808136,
  10596081.456567012,
  10599826.5625,
  10613496.742654076,
  10613535.9375,
  10660493.61125328,
  11266312.259867236,
  11270329.6875,
  11283058.054572864,
  11312946.205440858,
  11464919.360253304,
  11504846.195103101,
  11842039.1701484,
  11853501.70948552,
  12222704.818192063,
  12421846.875,
  13021132.129239254,
  13021236.963152034,
  14088962.023301644,
  15352404.085066361,
  16915983.181643557,
  18774067.1875,
  21136857.264631905,
  21150732.392295297,
  21150897.419354163,
  21151325.0,
  21153409.375,
  21154900.0,
  21155055.155825462,
  21156345.110535234,
  21156879.909279555,
  21158929.6821773,
  21159081.25,
  21161334.33923666,
  21161639.829834018,
  21161872.86379478,
  21171122.908343907,
  21171507.8125,
  21185470.3125],
 [44.1151146399842,
  19.06426387887476,
  72.9903032619747,
  38.13472018705268,
  6.325439695638338,
  24.822227514762076,
  11.746959293914506,
  8.462785199678626,
  10.954056047980188,
  12.775120256909076,
  10.179018807807488,
  80.61816726765272,
  30.76843308656368,
  35.224487691026255,
  59.15985230392127,
  14.4969808328599,
  6.11937455439975,
  28.21834154736074,
  58.469715920556624,
  7.271594547236137,
  39.89527397712712,
  77.67457694140487,
  83.48073164201625,
  12.325121847512246,
  37.14355983439159,
  60.76458705799514,
  57.60419889689573,
  58.90327392228315,
  55.142436694849835,
  5.448469472028692,
  5.266140685345057,
  26.699204358902662,
  60.43746248368338,
  12.284410691492242,
  54.28663991460195,
  19.65862071861211,
  14.313882741610369,
  9.802818611938054,
  12.5036899595841,
  56.865680112598156,
  12.956609693402438,
  72.57305452599724,
  18.08948475784833,
  14.45857897813693,
  13.058492048759996,
  56.59026424879305,
  55.5293496600556,
  26.83241979143611,
  9.660274476863322,
  5.277423733475989,
  19.921068089425823,
  13.13873509542886,
  47.81287930275692,
  27.574392241768237,
  9.046857268249944,
  5.6822278719798485,
  8.926881392463072,
  44.51478965713822,
  34.64802089011277,
  5.961892628678742,
  16.354084309577956,
  52.67755841132305,
  15.466787815070466,
  9.107993885711476,
  11.85406913600058,
  6.828359929736517,
  10.574297837662817,
  15.43331884419696,
  22.410809872894255,
  28.08726405099809,
  18.331705734777255,
  61.69834655165389,
  95.62651882812845,
  15.872847209055095,
  67.69293541037939,
  44.90747192526904,
  21.97033749124155,
  6.417337334516956,
  5.5673057473178735,
  53.182527603473076,
  52.789231939786504,
  5.339086024129061,
  18.196465446112608,
  8.129334092487895,
  59.35673075711627,
  42.681243451635034,
  88.56324384782026,
  104.62034746196505,
  63.05397070058917,
  48.0120416650076,
  12.04521380241369,
  28.705100479201604,
  12.601697867579386,
  6.072379288653127,
  65.81251228832417,
  43.855912782894464,
  11.542550740059664,
  47.16632222383809,
  19.154453525451583,
  55.40464740080633,
  24.66976290069403,
  106.93358363554164,
  11.680191733130869,
  103.20938263705304,
  14.409855284907763,
  6.394679956988659,
  49.697252317956526,
  36.87084176635848,
  25.082256671901476,
  74.84550407404369,
  12.968642171618674,
  20.14591987843851,
  5.640927607687405,
  38.44459780334823,
  9.342148523193076,
  13.1292821919657,
  5.997729949451351,
  21.627660035575317,
  39.39605867880935,
  12.400504983076265,
  31.825092144975738,
  82.28018200319113,
  12.432498111721994,
  13.077675289573866,
  45.401253421468454,
  43.34713211722807,
  65.99974675542694,
  25.24529936553993,
  21.57063169696749,
  41.82884419973743,
  37.86790300012656,
  8.508581803636059,
  39.97695041433258,
  16.86400992609002,
  45.05263228286865,
  14.284770459610836,
  34.9029385521692,
  29.498099299066983,
  44.6242773481123,
  9.021528827598978,
  33.47237665238689,
  24.043689620182455,
  56.52747540026464,
  32.37333376757789,
  64.53368214734593,
  69.62323855926952,
  19.3573075450932,
  69.18362975044886,
  6.189623192982941,
  17.127524836659706,
  94.71629865961143,
  17.73684762004018,
  128.74327013315232,
  5.758878386471944,
  12.118823143131959,
  8.84862856773068,
  5.649018814768355,
  11.445569218297525,
  6.10478367403138,
  37.13639519310876,
  47.819700209322924,
  19.734631054615683,
  37.30336617237097,
  23.47551075640914,
  30.07025744852214,
  16.84755776168545,
  15.08784426023395,
  12.716919264076385,
  8.971344247998388,
  108.70073457902065,
  70.31945360752536,
  17.823422729232263,
  27.922835638237142,
  38.30954713210878,
  45.710242536650156,
  13.28270205813922,
  110.16489039160074,
  38.80967442941125,
  62.45460963814571,
  32.25537420516064,
  5.655993084918314,
  7.764322597989998,
  14.213397081976634,
  38.76950877128941,
  19.05133554295636,
  78.7293859284863,
  5.911721051721169,
  20.339693383724885,
  16.25350287319941,
  16.545055094874876,
  16.598355608792932,
  46.509058977045704,
  9.175268653739458,
  87.88679300656369,
  75.78015652481228,
  26.645590883972027,
  5.564079130979073,
  75.26612782309063,
  70.80270410025568,
  7.053724254883401,
  39.548391287150906,
  6.313894561559217,
  25.343081561235234,
  8.227676788332154,
  27.917361251380587,
  44.031519762669234,
  47.62848283814515,
  8.018818709501337,
  77.50765631026522,
  16.02471019408454,
  87.5737807581747,
  14.773841686147545,
  19.952559522320428,
  65.96123218190318,
  12.044329315181297,
  11.898659173965536,
  29.5375383703038,
  7.1723515961332005,
  85.26475400441294,
  102.01775663013956,
  42.78794565295557,
  34.76736928452512,
  12.512683336939471,
  110.02825299799426,
  19.660748849851885,
  29.870095277623705,
  5.152679894997583,
  16.446819854709723,
  40.744659035337634,
  11.228643843107745,
  49.965925034787745,
  73.34811212073528,
  65.04944831015943,
  50.3945203782646,
  26.554241111965332,
  22.09252417201001,
  12.035614029447801,
  8.547456559337043,
  14.270531638196326,
  29.294444274851116,
  25.60338234327652,
  48.27599220935607,
  55.84730174540409,
  70.1200800220691,
  8.732255198758045,
  34.2173478441145,
  13.672847124999334,
  59.998747703778605,
  73.7752923924916,
  25.358856134644583,
  74.90341900558244,
  10.293895288654406,
  93.72194249286376,
  75.81884456317172,
  73.40699059816055,
  32.177150977121016,
  12.776389041989697,
  19.810395727676376,
  10.74391241859486,
  27.713833555461225,
  20.05444870105609,
  32.1465109091896,
  7.947235684971142,
  35.77949498082458,
  37.37167673741854,
  27.620099492935076,
  19.945952482713597,
  137.97105921619627,
  10.019067323920998,
  12.880306976007637,
  62.547279902002046,
  7.309140683844345,
  77.80728321381923,
  12.792094782112349,
  69.68048274756882,
  5.237319199379071,
  10.052132103849646,
  56.6528496199913,
  5.352054937020157,
  32.43299770226985,
  5.226886544443705,
  28.497173999539996,
  10.680012884976959,
  15.965543716271927,
  38.411749146696785,
  16.4763527377916,
  16.84595594224027,
  53.745381150535614,
  78.04525656715913,
  53.535962147417116,
  24.11018810374957,
  21.78980303610099,
  11.044363929260559,
  36.861906921835,
  17.267103874493166,
  9.091514360962233,
  7.3896592883333945,
  25.3945043178869,
  25.81073607046493,
  22.932084522689678,
  111.24845415770794,
  70.5817401917607,
  31.95795808170676,
  24.948974607958476,
  6.004840238454423,
  7.111847620016791,
  8.155116062592676,
  25.584005990752164,
  46.05357217033031,
  9.257994997089728,
  5.667996941309796,
  74.81569555068364,
  22.74395239701665,
  11.144702517986893,
  16.163918129957743,
  58.445715835220625,
  70.70819859710596,
  37.86236301120175,
  28.637458810901727,
  12.25708673720708,
  34.925675696221234,
  78.35091797027128,
  75.8754487392097,
  6.408258008946374,
  12.953418872988227,
  52.60236289268757,
  9.569493099430568,
  5.33976154816605,
  31.146579840221357,
  20.33239872212456,
  26.468911006709785,
  71.89694592944358,
  44.65323010397874,
  14.552860842324495,
  11.111230616288038,
  79.6473702368225,
  20.23597078053509,
  7.416905853498759,
  6.3058728422776325,
  21.744129771305836,
  69.46254431765593,
  31.070871907119106,
  14.540261200490216,
  93.1089906854249,
  5.944677707199027,
  47.70924175827307,
  7.084322631769302,
  5.3394845099884725,
  63.58578075136934,
  88.63714682439617,
  19.039342417512334,
  7.414296949653421,
  23.44119228206218,
  10.921669192712823,
  16.277486191274637,
  11.029497843892338,
  5.94734689058128,
  12.516180182975104,
  20.015444147434213,
  7.011076823634419,
  31.817788417424886,
  10.60264287948333,
  17.28167714313644,
  6.422637911135516,
  120.12891584833949,
  11.022768610164345,
  46.08537567168395,
  51.88826173624381,
  24.22169463624742,
  26.73061451819018,
  92.18299088748196,
  14.95048459754247,
  38.327774785568856,
  46.61623283478523,
  40.42659496383164,
  5.356872034598011,
  100.77168584833755,
  9.209142916726053,
  18.308480791907638,
  19.479335078374206,
  5.696975610512134,
  36.41198809006641,
  52.039281505551386,
  65.6319387126819,
  9.323208525588003,
  72.47052907164043,
  50.936769141401044,
  41.113564597624965,
  50.68197049373533,
  5.102187912127682,
  18.608548102569607,
  68.73205538791473,
  48.80228820514137,
  70.69101785032288,
  54.61711063185947,
  37.663163279959825,
  75.57348890239108,
  8.306995287967535,
  43.63043132243486,
  30.936921061156582,
  118.11881886514102,
  6.360512981174484,
  12.97709501116696,
  20.065126038915896,
  6.7938184861573925,
  17.473040884504,
  52.702457743103146,
  31.980294675116177,
  67.87911720947356,
  15.313724198302006,
  5.429563640523038,
  16.79555780845051,
  85.10854728562053,
  15.633090702010843,
  60.15107090139588,
  11.711559948534903,
  45.44521950597776,
  51.97451049250542,
  27.17248980605158,
  55.94857937069234,
  45.875253225879774,
  8.190400442147737,
  12.178940989205541,
  7.972786132607971,
  7.3611880950442625,
  32.61733806122776,
  86.18945987039602,
  6.5032304853187375,
  80.76789131992854,
  21.590644063066907,
  49.73562227932904,
  16.65894442596359,
  53.996355074428614,
  38.270144423115774,
  193.64953557646763,
  48.35176764591674,
  72.09628287322178,
  12.654069756665395,
  38.8339941460539,
  12.595519237052462,
  122.17652320991567,
  32.78708495522601,
  6.246340876534137,
  21.35276860440632,
  103.65006610807802,
  18.37950320939907,
  65.30875245199718,
  31.777682584271428,
  11.9934698260187,
  77.92956892869215,
  36.49721293314736,
  71.239265142184,
  122.0610606632733,
  35.907863728312236,
  6.204573074331645,
  5.1524370709084275,
  55.804340067416746,
  20.779583193009863,
  19.75781849101845,
  17.571055407258257,
  52.55299288463493,
  50.14998879565123,
  35.65022850058403,
  85.41052275778333,
  56.2555650788846,
  57.23849621100226,
  9.040988273203602,
  5.393435592483144,
  17.884595057155213,
  57.504636863031024,
  69.78241276972112,
  11.629865209189063,
  60.89575977819283,
  30.90328360084,
  29.8121935448994,
  8.987816759050844,
  15.62090695596031,
  49.25848111915589,
  46.702426614008445,
  16.12379318123065,
  117.69004611416153,
  37.48535648654719,
  10.288642889780611,
  79.41222840551767,
  33.23594943738304,
  30.341678915904232,
  16.3865490921623,
  37.26010689837779,
  13.09119309603918,
  52.93914847664452,
  15.369360924522296,
  47.46616114642142,
  13.219065535564372,
  94.96129393064626,
  9.065508855893341,
  107.57836860958712,
  86.15077401788794,
  9.73163239302457,
  17.396314484373416,
  14.167533245697808,
  32.058409728159546,
  6.825834002638735,
  51.239449661593525,
  58.93740336626033,
  23.655523237322512,
  35.64049152404362,
  16.32943919920143,
  149.95813784939372,
  75.92476748778904,
  14.806688492056985,
  11.42172709769792,
  12.986917651710812,
  58.892394085191164,
  40.04194519137185,
  22.50355454347949,
  17.71787837931663,
  7.01627802521345,
  7.820297110945039,
  125.06472974834055,
  47.76856722361478,
  7.57037303484295,
  19.08473940808929,
  41.835368283748984,
  7.502814439472981,
  9.238889874409926,
  23.33341351072511,
  29.93371572184674,
  12.716583438589957,
  52.81800808660188,
  97.58754923380447,
  104.53354352277216,
  13.628983579329898,
  15.345176395034061,
  13.253673447037574,
  49.44141229447902,
  86.91245433298742,
  81.4023266407402,
  24.244521008207688,
  54.89592184591685,
  13.623843696980469,
  17.775820838582288,
  26.735643861058588,
  31.893504129394383,
  62.15543177348213,
  12.740465708216751,
  90.58736324751223,
  13.36420624141576,
  63.514639157202254,
  29.343776749341284,
  9.104132199686298,
  24.864152586811628,
  110.25589949515893,
  80.38591094106349,
  42.76975445077802,
  26.92601703749842,
  5.139137018060453,
  18.991205300705037,
  30.743551390884154,
  10.100506165017265,
  188.28953587846502,
  26.5083642002,
  7.457664106588232,
  85.2037456826129,
  60.294774072359175,
  37.186989379124206,
  7.2459628458718335,
  39.30215180743267,
  64.30323022022141,
  23.040552951353664,
  69.11603749455873,
  54.62361005109108,
  76.65400538607057,
  82.76549022260227,
  110.53504336507139,
  130.87130999124744,
  6.072808879926113,
  58.111851209136226,
  105.35752406971719,
  34.38036317554156,
  28.174263923786683,
  32.93411117660727,
  32.99055313746904,
  82.34374965142545,
  68.38648816663668,
  13.155785748994665,
  91.90680382008566,
  31.983977102081052,
  12.576389497423712,
  91.49274019929285,
  8.710000511982177,
  124.78674027659044,
  21.191379185396155,
  58.310930454959674,
  48.11106634708665,
  78.71454323431442,
  8.990780491054972,
  21.257192502627564,
  44.14320460383532,
  68.64433277544835,
  14.795440048944869,
  52.42760322465816,
  5.684284310337724,
  13.664912508120777,
  80.15294075129698,
  59.25054773553046,
  20.806311849834053,
  84.1034830302795,
  29.654837220431865,
  97.93362790216196,
  77.60043279710577,
  23.078720834992595,
  20.320338310963102,
  8.564982134652046,
  74.02373323676741,
  7.915265259407697,
  8.905475779391224,
  89.78927083052886,
  33.58490252823504,
  15.15277463915615,
  7.966275210505147,
  9.982631873396429,
  48.51825091615752,
  24.3351324715449,
  37.26790851629909,
  14.565572918921177,
  81.06153852630564,
  45.169041940064815,
  66.2566467434483,
  42.746235222276084,
  121.97616513811549,
  8.420166007439002,
  6.4248491233608345,
  36.864613302282265,
  140.99979581251912,
  97.79007539472119,
  120.14393530791497,
  16.9697638552666,
  71.91385799543811,
  5.241716525991229,
  7.983194593430584,
  9.351960656227808,
  61.69359033606149,
  27.516682857566728,
  51.06542241085249,
  66.79630259813428,
  14.634537460080129,
  90.59091303014901,
  61.085293292397424,
  11.593442606607278,
  7.218584255612209,
  7.310268912187566,
  11.084565318484033,
  20.0536050738007,
  30.29319862748777,
  54.913324375700626,
  60.3459176569474,
  70.39077208628738,
  36.229832084454166,
  29.14664949548409,
  73.83968586297199,
  25.61621942376874,
  16.977445889745475,
  15.658833202617416,
  35.54911488093509,
  23.770979608342763,
  35.20370765171865,
  53.877788441314955,
  70.65449756901506,
  46.24701145989481,
  7.252463696345999,
  7.467963201015802,
  24.255745544995037,
  54.7784164460617,
  119.0141470367416,
  70.16071427657803,
  28.32061750196084,
  12.370195205868335,
  75.99137627421038,
  7.096127943969323,
  70.97090832941112,
  38.81718835026143,
  63.86188918029771,
  87.56821433702686,
  27.735560469253016,
  12.712045199423539,
  79.77510687315201,
  35.228424982198234,
  34.70263446960864,
  71.83231182535239,
  55.88773199103073,
  36.908583024142715,
  30.894621351864586,
  81.93839162962121,
  36.04662453303813,
  50.3789414816418,
  39.620746577144274,
  5.881058026258764,
  31.25339257157878,
  74.79105159945803,
  34.99223608786868,
  67.50951682876378,
  19.231094599960265,
  15.791986493508077,
  57.56086448587226,
  26.157807521580434,
  100.69714851916748,
  79.72215670357868,
  7.5186349885285635,
  61.68177699135296,
  11.528407073010875,
  56.22668864926588,
  5.051980027559639,
  5.882305076554136,
  6.280895447041716,
  11.3654990365045,
  58.890535397130265,
  64.39763481596512,
  12.850678556387784,
  5.032611640767161,
  21.728739094415914,
  40.9365438830743,
  9.664339593854667,
  47.954509770718566,
  21.03970765732501,
  22.327828817994206,
  25.229019999770614,
  54.23928613320853,
  64.52049644226265,
  186.17141128534433,
  45.186552877895515,
  29.806307565536624,
  8.40744664974339,
  58.31608778728173,
  72.80306627412551,
  30.654815742712234,
  16.34176916151307,
  13.801535634857345,
  26.63197255242998,
  35.26472329384794,
  56.45371909211843])

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

([2128322.8820732716,
  2673853.9038674626,
  2786856.25,
  2789778.125,
  2805233.259980844,
  2807010.716995365,
  2808506.6130725285,
  2812117.506482652,
  2874049.8859884213,
  2889719.5982088023,
  2889974.2945079044,
  2890028.125,
  2891280.9293882987,
  2892351.5625,
  2893143.3147037867,
  2905490.5032702447,
  2905515.9204744357,
  2906018.734324037,
  2948343.75,
  3034213.041952173,
  3038439.0625,
  3141742.1875,
  3424040.625,
  6785328.776244561,
  6785493.381396405,
  6791056.25,
  6795218.742225575,
  6803876.5625,
  6861424.907240535,
  6866769.9989933055,
  6866830.420718658,
  6879946.775391034,
  6940407.8125,
  6953758.444825136,
  6977915.625,
  6978375.856584645,
  6982520.137476653,
  6982803.634441398,
  6987272.607001323,
  6988256.8744002385,
  6988266.546827361,
  6988965.625,
  6990647.289649813,
  6992557.203741236,
  6993033.774311789,
  7020286.637365565,
  7110217.1875,
  7114679.595704689,
  7115817.552414775,
  7116227.676238994,
  7117364.95151519,
  7118310.367652139,
  7118339.0625,
  7122178.9844065495,
  7123494.955652764,
  7125337.107539448,
  7125688.541932925,
  7129987.5,
  7130548.4375,
  7146026.817144993,
  7157582.864702709,
  7169130.601346257,
  7172898.599309346,
  7173441.9109135885,
  7174609.7948041605,
  7176262.816393783,
  7188563.677035553,
  7192779.81081335,
  7193623.024825048,
  7193984.352505192,
  7199272.54899749,
  7199509.375,
  7210189.0625,
  7219604.912121246,
  7220298.4375,
  7220809.009714929,
  7255282.293938559,
  7269896.415108262,
  7292217.873084736,
  7301620.525656571,
  7303248.4375,
  7306081.167012519,
  7308306.903424871,
  7308351.252823675,
  7310126.5625,
  7310192.1875,
  7312481.25,
  7318070.146381164,
  7318290.625,
  7320417.1875,
  7331122.39212852,
  7351126.5625,
  7352783.587341864,
  7352798.284965201,
  7357039.0625,
  7366504.6875,
  7368967.886011667,
  7369167.9398481585,
  7369285.676346161,
  7369331.25,
  7369440.2057869695,
  7369593.75,
  7370225.326871223,
  7370403.450339849,
  7370856.366219008,
  7371115.00902001,
  7371560.9375,
  7372068.883017307,
  7375343.5475491695,
  7375645.3125,
  7376091.551598742,
  7376977.478443644,
  7377316.100927785,
  7377719.397872175,
  7379546.887680742,
  7380329.11404829,
  7380764.708721365,
  7380780.7633840125,
  7380978.125,
  7381117.873272525,
  7381319.525419512,
  7381475.772936901,
  7381608.512732085,
  7381697.844236114,
  7382160.9375,
  7382195.188103089,
  7382382.979949364,
  7383914.557640953,
  7384037.262309488,
  7384590.995170809,
  7385517.137191394,
  7386523.834975215,
  7386632.8125,
  7407451.204093161,
  7411083.48540908,
  7416076.917212375,
  7420263.841231703,
  7421003.125,
  7423408.618539594,
  7424423.0467876475,
  7424773.4375,
  7443801.285092924,
  7473906.030115981,
  7498054.6875,
  7498323.4375,
  7514310.9375,
  7514313.62134014,
  7515864.0625,
  7517616.54779265,
  7517636.413010648,
  7517884.767175848,
  7518367.767904604,
  7518534.561223161,
  7518955.476446981,
  7519572.428332865,
  7520449.367957458,
  7520524.0528527545,
  7523556.444266042,
  7524612.658927334,
  7529778.125,
  7531141.88947555,
  7531305.553590942,
  7532726.177962356,
  7532863.096114185,
  7532898.4375,
  7533681.611589373,
  7538795.601631367,
  7540046.136425193,
  7553094.853595726,
  7553971.859383369,
  7566477.775048121,
  7668799.085138062,
  7671643.07313268,
  7676275.0,
  7676429.536808817,
  7676441.053957746,
  7676862.5,
  7677190.317446775,
  7682425.0,
  7683807.8125,
  7683960.109159015,
  7684184.7994470615,
  7684206.194942015,
  7686345.07821117,
  7687282.027107207,
  7690796.875,
  7690810.614015361,
  7691735.302900683,
  7692922.4904334545,
  7692935.517091191,
  7704117.347873407,
  7705287.5,
  7705876.371888888,
  7706700.0,
  7708684.375,
  7708724.897220039,
  7709552.185082504,
  7710551.5625,
  7711298.249613774,
  7712076.298431986,
  7712229.6875,
  7714553.485871672,
  7716235.752709214,
  7719565.511932308,
  7720832.781053785,
  7720887.5,
  7723117.1875,
  7725433.831847788,
  7735924.301922041,
  7747629.945353442,
  7748082.513033068,
  7750192.142803754,
  7751469.642833331,
  7751663.953054529,
  7754494.266689285,
  7754900.442947881,
  7755684.375,
  7755784.56515661,
  7755860.9375,
  7757521.875,
  7761043.75,
  7761525.0,
  7762597.610975244,
  7763178.125,
  7763250.856852085,
  7763937.710940844,
  7764038.509513045,
  7764114.408345873,
  7764550.0,
  7764590.688630665,
  7764998.4375,
  7766364.043550493,
  7767175.0,
  7767673.4375,
  7768175.148553596,
  7769952.683060873,
  7770311.449855896,
  7770827.02526446,
  7773127.937071002,
  7775690.625,
  7783923.269901901,
  7786251.211772154,
  7787733.52179854,
  7788352.344947424,
  7809994.5073372945,
  7821059.375,
  7822560.559423309,
  7822659.139663249,
  7825477.844137313,
  7825526.378763232,
  7825692.1875,
  7825704.079392292,
  7829315.625,
  7833160.406130971,
  7833259.375,
  7833314.871079315,
  7834111.786234456,
  7834416.500473576,
  7834965.364753085,
  7835912.519376425,
  7836031.917813412,
  7838401.5625,
  7863032.110494282,
  7921702.291830664,
  7945260.601712729,
  7947456.784712147,
  7955254.923686325,
  7957453.125,
  7968525.714854295,
  7975208.703757318,
  7978595.168187187,
  8015402.231421224,
  8015760.449504381,
  8016258.667728237,
  8016607.5604044385,
  8020241.9054473275,
  8022636.682392288,
  8023439.0625,
  8024208.105269995,
  8024979.1824363535,
  8025038.3967448,
  8025178.929093858,
  8034829.380865554,
  8040323.564832087,
  8040382.665801466,
  8041511.289273224,
  8045439.8724204935,
  8045841.311765584,
  8046823.34052544,
  8048256.93844812,
  8048334.327577891,
  8048555.779237429,
  8048890.5710173445,
  8049136.474058713,
  8049313.746090505,
  8049513.7644822635,
  8049647.805680556,
  8049743.568639341,
  8049929.533299217,
  8049941.820588816,
  8049960.9375,
  8049969.91470732,
  8050437.5,
  8050465.413352307,
  8050519.471912093,
  8050848.400838206,
  8051066.9303468885,
  8051665.456136617,
  8052010.173352713,
  8052112.4118886795,
  8052207.185729658,
  8052296.875,
  8052500.458025544,
  8052754.698346187,
  8053609.662566652,
  8056897.520079125,
  8057525.195532655,
  8057626.824388124,
  8058221.875,
  8058278.744517987,
  8062337.5,
  8064962.433896966,
  8071439.0625,
  8117371.352970065,
  8121045.683174125,
  8131168.75,
  8132333.0812669145,
  8133817.124324756,
  8137193.75,
  8140355.6252022395,
  8143071.975665295,
  8177413.946792574,
  8182925.0,
  8183262.266429158,
  8183282.801933412,
  8184017.1875,
  8262919.676514626,
  8319524.453658248,
  8320258.112596558,
  8587742.383229863,
  8587789.0625,
  8596391.256393256,
  8608374.578842437,
  8608943.83626232,
  8610839.122841472,
  8611495.241446324,
  8611543.606846513,
  8619893.678702816,
  8634428.929455204,
  8665534.375,
  8666536.712971518,
  8668197.883039132,
  8859658.38476883,
  8881511.71554166,
  8883381.561294029,
  8900606.18624703,
  8900914.153240364,
  8902285.112518042,
  8902397.981962793,
  8902476.125288535,
  8903878.070447821,
  8904249.635932965,
  8904428.693151668,
  8906205.826134952,
  8908557.003073901,
  8935698.369310148,
  8935764.0625,
  8936295.3125,
  8946430.985996012,
  8947640.816003395,
  8948671.678208986,
  8948756.252352752,
  8948870.647372551,
  8949218.535801215,
  8949322.570562083,
  8949475.429906635,
  8953466.860424375,
  8956999.85178064,
  8957852.23299286,
  8962763.333784113,
  8962821.720186707,
  8963335.9375,
  8992293.75,
  8997832.8125,
  8997859.365292648,
  9010415.625,
  9017079.6875,
  9066584.21633308,
  9067332.8125,
  9067412.688677158,
  9072793.062565196,
  9075806.25,
  9080371.875,
  9083409.375,
  9153968.440669978,
  9165490.625,
  9168406.25,
  9179781.288412115,
  9189612.5,
  9192460.363105677,
  9194167.1875,
  9194639.767033448,
  9202335.3049427,
  9249010.281356322,
  9249443.899958383,
  9250288.036185153,
  9251842.404319722,
  9252202.637876833,
  9252948.4375,
  9255107.906843672,
  9255286.727617305,
  9255894.164045552,
  9255990.625,
  9257420.48283648,
  9257495.3125,
  9260731.593960673,
  9266870.254935114,
  9269437.747614633,
  9270153.577840896,
  9271154.156907385,
  9271160.328719957,
  9271168.563589526,
  9274792.769897306,
  9280541.542912329,
  9281732.089715682,
  9329502.853522174,
  9338001.538835473,
  9338771.977189163,
  9340776.099574013,
  9341239.663460804,
  9342315.625,
  9409865.91240456,
  9429447.17349302,
  9444318.75,
  9481693.814526606,
  9508645.3125,
  9510992.89959578,
  9514419.547320418,
  9516304.6875,
  9519653.89717113,
  9521137.5,
  9527107.8125,
  9595102.958110666,
  9711871.360816404,
  9815831.592904603,
  9823874.118673626,
  9823880.209583577,
  9824207.524427542,
  9825361.428761894,
  9830106.798444437,
  9831970.894541385,
  9840715.625,
  9843538.27984846,
  9892337.337896347,
  9938704.17520631,
  9939038.498756967,
  9957040.625,
  9980320.431322752,
  9980928.993972246,
  9982367.740320219,
  9992629.70397702,
  9996173.253948873,
  10058179.216325019,
  10082157.8125,
  10099350.0,
  10149936.94229971,
  10150868.385338144,
  10155579.555075698,
  10159236.54191782,
  10164091.137727635,
  10164912.817033634,
  10166479.98664615,
  10166712.5,
  10167534.375,
  10168466.908327503,
  10169152.688641166,
  10173309.542516291,
  10174112.5,
  10175453.869309856,
  10176014.34157418,
  10176147.88502453,
  10176196.875,
  10177354.075021673,
  10177476.5625,
  10177632.8125,
  10178531.025056299,
  10178925.333659176,
  10179859.858418666,
  10180989.232086599,
  10181584.375,
  10181921.982898805,
  10181993.75,
  10182011.960467963,
  10183334.375,
  10184765.292858748,
  10185068.061137648,
  10185526.5625,
  10185625.438219737,
  10185627.08548888,
  10187477.90595186,
  10187887.5,
  10188219.054379314,
  10191548.4375,
  10197010.9375,
  10197971.567508955,
  10198071.158335626,
  10201736.259498581,
  10206708.190323729,
  10210881.25,
  10214397.337914497,
  10215226.742799269,
  10215227.496320343,
  10246248.4375,
  10261672.217736252,
  10264516.783159353,
  10269250.577283144,
  10285512.084199846,
  10286591.187746808,
  10287396.875,
  10287406.25,
  10290715.53321069,
  10292374.219148992,
  10345453.125,
  10366834.088538708,
  10408426.451024199,
  10415284.04161117,
  10415443.75,
  10415888.323784096,
  10417464.0625,
  10418431.353151234,
  10420387.354281852,
  10420738.669505117,
  10423329.789523406,
  10424124.473360265,
  10424710.9375,
  10424845.3125,
  10424987.710396757,
  10425013.556489943,
  10425446.875,
  10429316.855651218,
  10429819.336790143,
  10443600.345360486,
  10444951.5625,
  10446738.61328377,
  10449353.958790153,
  10450143.916567834,
  10455825.27547293,
  10498625.917974658,
  10516085.9375,
  10522646.593404826,
  10533317.558245003,
  10533732.62736665,
  10537687.5,
  10543342.1875,
  10543373.686645012,
  10554498.837182047,
  10558617.956379369,
  10563234.375,
  10564486.725717241,
  10565403.125,
  10565564.93853149,
  10565673.016238809,
  10565896.875,
  10565998.4375,
  10566718.41956189,
  10566728.578391833,
  10566806.25,
  10566965.625,
  10567043.392290866,
  10567181.25,
  10567260.216435133,
  10567451.5625,
  10567487.644836836,
  10567582.8125,
  10567643.394951439,
  10567842.039170349,
  10567868.75,
  10568107.741681658,
  10568240.494677877,
  10568347.67987977,
  10568351.660458645,
  10568357.618779495,
  10568368.677976139,
  10568385.9375,
  10568680.686633037,
  10568684.375,
  10568885.821263582,
  10569118.65168374,
  10569152.401719801,
  10569177.285681259,
  10569251.5625,
  10569361.887776734,
  10569459.375,
  10569520.3125,
  10569578.178846367,
  10569633.942680601,
  10569811.97602166,
  10569895.3125,
  10570029.6875,
  10570227.644355712,
  10570479.6875,
  10570507.453489214,
  10570573.056735102,
  10570868.75,
  10570881.315848306,
  10571049.428281816,
  10571300.0,
  10571733.477127165,
  10571848.850952918,
  10572145.3125,
  10572149.64550355,
  10572169.68631037,
  10572219.241805673,
  10572751.5625,
  10573114.849994272,
  10573329.136818837,
  10573368.040424915,
  10573884.375,
  10573885.991845291,
  10574179.028357456,
  10575708.493846755,
  10575820.53458742,
  10575939.11864759,
  10575956.25,
  10576259.27274382,
  10576654.6875,
  10577139.412696933,
  10577323.4375,
  10577675.0,
  10577730.830906747,
  10580065.869354483,
  10582627.71610995,
  10582940.391953839,
  10583024.784083473,
  10583235.940249423,
  10583451.5625,
  10583580.134203324,
  10583585.9375,
  10583620.189527338,
  10583891.566072026,
  10583893.774212722,
  10584210.865249392,
  10584312.750949088,
  10584386.848524414,
  10584391.203458857,
  10584410.854561605,
  10584426.5625,
  10584464.0625,
  10584664.341340942,
  10584858.169555966,
  10584931.922686916,
  10584995.284722965,
  10585044.654981168,
  10585318.405018348,
  10585360.254112825,
  10585482.8125,
  10585548.867458332,
  10585781.25,
  10585783.46990954,
  10585918.75,
  10586100.584047938,
  10586141.974293912,
  10586145.488775168,
  10586187.5,
  10586230.686924592,
  10586390.625,
  10586408.631720155,
  10586418.75,
  10586517.1875,
  10586555.16350034,
  10586812.109829774,
  10586925.553527204,
  10587115.625,
  10587312.87950314,
  10587501.5625,
  10587567.177092817,
  10587593.492776426,
  10587665.082314996,
  10587703.02778036,
  10587718.75,
  10587928.510295745,
  10588153.125,
  10588181.28257486,
  10588193.708831562,
  10588214.884376239,
  10588704.6875,
  10588765.625,
  10589155.912788745,
  10589476.261585321,
  10589704.589205928,
  10589765.625,
  10589773.4375,
  10590460.9375,
  10590932.290547805,
  10591484.375,
  10592675.599312136,
  10593053.992193552,
  10593148.037111022,
  10594687.197271833,
  10594689.0625,
  10595043.75,
  10596071.113808136,
  10596081.456567012,
  10599826.5625,
  10613496.742654076,
  10613535.9375,
  10660493.61125328,
  11266312.259867236,
  11270329.6875,
  11283058.054572864,
  11312946.205440858,
  11464919.360253304,
  11504846.195103101,
  11842039.1701484,
  11853501.70948552,
  12222704.818192063,
  12421846.875,
  13021132.129239254,
  13021236.963152034,
  14088962.023301644,
  15352404.085066361,
  16915983.181643557,
  18774067.1875,
  21136857.264631905,
  21150732.392295297,
  21150897.419354163,
  21151325.0,
  21153409.375,
  21154900.0,
  21155055.155825462,
  21156345.110535234,
  21156879.909279555,
  21158929.6821773,
  21159081.25,
  21161334.33923666,
  21161639.829834018,
  21161872.86379478,
  21171122.908343907,
  21171507.8125,
  21185470.3125],
 [44.1151146399842,
  19.06426387887476,
  72.9903032619747,
  38.13472018705268,
  6.325439695638338,
  24.822227514762076,
  11.746959293914506,
  8.462785199678626,
  10.954056047980188,
  12.775120256909076,
  10.179018807807488,
  80.61816726765272,
  30.76843308656368,
  35.224487691026255,
  59.15985230392127,
  14.4969808328599,
  6.11937455439975,
  28.21834154736074,
  58.469715920556624,
  7.271594547236137,
  39.89527397712712,
  77.67457694140487,
  83.48073164201625,
  12.325121847512246,
  37.14355983439159,
  60.76458705799514,
  57.60419889689573,
  58.90327392228315,
  55.142436694849835,
  5.448469472028692,
  5.266140685345057,
  26.699204358902662,
  60.43746248368338,
  12.284410691492242,
  54.28663991460195,
  19.65862071861211,
  14.313882741610369,
  9.802818611938054,
  12.5036899595841,
  56.865680112598156,
  12.956609693402438,
  72.57305452599724,
  18.08948475784833,
  14.45857897813693,
  13.058492048759996,
  56.59026424879305,
  55.5293496600556,
  26.83241979143611,
  9.660274476863322,
  5.277423733475989,
  19.921068089425823,
  13.13873509542886,
  47.81287930275692,
  27.574392241768237,
  9.046857268249944,
  5.6822278719798485,
  8.926881392463072,
  44.51478965713822,
  34.64802089011277,
  5.961892628678742,
  16.354084309577956,
  52.67755841132305,
  15.466787815070466,
  9.107993885711476,
  11.85406913600058,
  6.828359929736517,
  10.574297837662817,
  15.43331884419696,
  22.410809872894255,
  28.08726405099809,
  18.331705734777255,
  61.69834655165389,
  95.62651882812845,
  15.872847209055095,
  67.69293541037939,
  44.90747192526904,
  21.97033749124155,
  6.417337334516956,
  5.5673057473178735,
  53.182527603473076,
  52.789231939786504,
  5.339086024129061,
  18.196465446112608,
  8.129334092487895,
  59.35673075711627,
  42.681243451635034,
  88.56324384782026,
  104.62034746196505,
  63.05397070058917,
  48.0120416650076,
  12.04521380241369,
  28.705100479201604,
  12.601697867579386,
  6.072379288653127,
  65.81251228832417,
  43.855912782894464,
  11.542550740059664,
  47.16632222383809,
  19.154453525451583,
  55.40464740080633,
  24.66976290069403,
  106.93358363554164,
  11.680191733130869,
  103.20938263705304,
  14.409855284907763,
  6.394679956988659,
  49.697252317956526,
  36.87084176635848,
  25.082256671901476,
  74.84550407404369,
  12.968642171618674,
  20.14591987843851,
  5.640927607687405,
  38.44459780334823,
  9.342148523193076,
  13.1292821919657,
  5.997729949451351,
  21.627660035575317,
  39.39605867880935,
  12.400504983076265,
  31.825092144975738,
  82.28018200319113,
  12.432498111721994,
  13.077675289573866,
  45.401253421468454,
  43.34713211722807,
  65.99974675542694,
  25.24529936553993,
  21.57063169696749,
  41.82884419973743,
  37.86790300012656,
  8.508581803636059,
  39.97695041433258,
  16.86400992609002,
  45.05263228286865,
  14.284770459610836,
  34.9029385521692,
  29.498099299066983,
  44.6242773481123,
  9.021528827598978,
  33.47237665238689,
  24.043689620182455,
  56.52747540026464,
  32.37333376757789,
  64.53368214734593,
  69.62323855926952,
  19.3573075450932,
  69.18362975044886,
  6.189623192982941,
  17.127524836659706,
  94.71629865961143,
  17.73684762004018,
  128.74327013315232,
  5.758878386471944,
  12.118823143131959,
  8.84862856773068,
  5.649018814768355,
  11.445569218297525,
  6.10478367403138,
  37.13639519310876,
  47.819700209322924,
  19.734631054615683,
  37.30336617237097,
  23.47551075640914,
  30.07025744852214,
  16.84755776168545,
  15.08784426023395,
  12.716919264076385,
  8.971344247998388,
  108.70073457902065,
  70.31945360752536,
  17.823422729232263,
  27.922835638237142,
  38.30954713210878,
  45.710242536650156,
  13.28270205813922,
  110.16489039160074,
  38.80967442941125,
  62.45460963814571,
  32.25537420516064,
  5.655993084918314,
  7.764322597989998,
  14.213397081976634,
  38.76950877128941,
  19.05133554295636,
  78.7293859284863,
  5.911721051721169,
  20.339693383724885,
  16.25350287319941,
  16.545055094874876,
  16.598355608792932,
  46.509058977045704,
  9.175268653739458,
  87.88679300656369,
  75.78015652481228,
  26.645590883972027,
  5.564079130979073,
  75.26612782309063,
  70.80270410025568,
  7.053724254883401,
  39.548391287150906,
  6.313894561559217,
  25.343081561235234,
  8.227676788332154,
  27.917361251380587,
  44.031519762669234,
  47.62848283814515,
  8.018818709501337,
  77.50765631026522,
  16.02471019408454,
  87.5737807581747,
  14.773841686147545,
  19.952559522320428,
  65.96123218190318,
  12.044329315181297,
  11.898659173965536,
  29.5375383703038,
  7.1723515961332005,
  85.26475400441294,
  102.01775663013956,
  42.78794565295557,
  34.76736928452512,
  12.512683336939471,
  110.02825299799426,
  19.660748849851885,
  29.870095277623705,
  5.152679894997583,
  16.446819854709723,
  40.744659035337634,
  11.228643843107745,
  49.965925034787745,
  73.34811212073528,
  65.04944831015943,
  50.3945203782646,
  26.554241111965332,
  22.09252417201001,
  12.035614029447801,
  8.547456559337043,
  14.270531638196326,
  29.294444274851116,
  25.60338234327652,
  48.27599220935607,
  55.84730174540409,
  70.1200800220691,
  8.732255198758045,
  34.2173478441145,
  13.672847124999334,
  59.998747703778605,
  73.7752923924916,
  25.358856134644583,
  74.90341900558244,
  10.293895288654406,
  93.72194249286376,
  75.81884456317172,
  73.40699059816055,
  32.177150977121016,
  12.776389041989697,
  19.810395727676376,
  10.74391241859486,
  27.713833555461225,
  20.05444870105609,
  32.1465109091896,
  7.947235684971142,
  35.77949498082458,
  37.37167673741854,
  27.620099492935076,
  19.945952482713597,
  137.97105921619627,
  10.019067323920998,
  12.880306976007637,
  62.547279902002046,
  7.309140683844345,
  77.80728321381923,
  12.792094782112349,
  69.68048274756882,
  5.237319199379071,
  10.052132103849646,
  56.6528496199913,
  5.352054937020157,
  32.43299770226985,
  5.226886544443705,
  28.497173999539996,
  10.680012884976959,
  15.965543716271927,
  38.411749146696785,
  16.4763527377916,
  16.84595594224027,
  53.745381150535614,
  78.04525656715913,
  53.535962147417116,
  24.11018810374957,
  21.78980303610099,
  11.044363929260559,
  36.861906921835,
  17.267103874493166,
  9.091514360962233,
  7.3896592883333945,
  25.3945043178869,
  25.81073607046493,
  22.932084522689678,
  111.24845415770794,
  70.5817401917607,
  31.95795808170676,
  24.948974607958476,
  6.004840238454423,
  7.111847620016791,
  8.155116062592676,
  25.584005990752164,
  46.05357217033031,
  9.257994997089728,
  5.667996941309796,
  74.81569555068364,
  22.74395239701665,
  11.144702517986893,
  16.163918129957743,
  58.445715835220625,
  70.70819859710596,
  37.86236301120175,
  28.637458810901727,
  12.25708673720708,
  34.925675696221234,
  78.35091797027128,
  75.8754487392097,
  6.408258008946374,
  12.953418872988227,
  52.60236289268757,
  9.569493099430568,
  5.33976154816605,
  31.146579840221357,
  20.33239872212456,
  26.468911006709785,
  71.89694592944358,
  44.65323010397874,
  14.552860842324495,
  11.111230616288038,
  79.6473702368225,
  20.23597078053509,
  7.416905853498759,
  6.3058728422776325,
  21.744129771305836,
  69.46254431765593,
  31.070871907119106,
  14.540261200490216,
  93.1089906854249,
  5.944677707199027,
  47.70924175827307,
  7.084322631769302,
  5.3394845099884725,
  63.58578075136934,
  88.63714682439617,
  19.039342417512334,
  7.414296949653421,
  23.44119228206218,
  10.921669192712823,
  16.277486191274637,
  11.029497843892338,
  5.94734689058128,
  12.516180182975104,
  20.015444147434213,
  7.011076823634419,
  31.817788417424886,
  10.60264287948333,
  17.28167714313644,
  6.422637911135516,
  120.12891584833949,
  11.022768610164345,
  46.08537567168395,
  51.88826173624381,
  24.22169463624742,
  26.73061451819018,
  92.18299088748196,
  14.95048459754247,
  38.327774785568856,
  46.61623283478523,
  40.42659496383164,
  5.356872034598011,
  100.77168584833755,
  9.209142916726053,
  18.308480791907638,
  19.479335078374206,
  5.696975610512134,
  36.41198809006641,
  52.039281505551386,
  65.6319387126819,
  9.323208525588003,
  72.47052907164043,
  50.936769141401044,
  41.113564597624965,
  50.68197049373533,
  5.102187912127682,
  18.608548102569607,
  68.73205538791473,
  48.80228820514137,
  70.69101785032288,
  54.61711063185947,
  37.663163279959825,
  75.57348890239108,
  8.306995287967535,
  43.63043132243486,
  30.936921061156582,
  118.11881886514102,
  6.360512981174484,
  12.97709501116696,
  20.065126038915896,
  6.7938184861573925,
  17.473040884504,
  52.702457743103146,
  31.980294675116177,
  67.87911720947356,
  15.313724198302006,
  5.429563640523038,
  16.79555780845051,
  85.10854728562053,
  15.633090702010843,
  60.15107090139588,
  11.711559948534903,
  45.44521950597776,
  51.97451049250542,
  27.17248980605158,
  55.94857937069234,
  45.875253225879774,
  8.190400442147737,
  12.178940989205541,
  7.972786132607971,
  7.3611880950442625,
  32.61733806122776,
  86.18945987039602,
  6.5032304853187375,
  80.76789131992854,
  21.590644063066907,
  49.73562227932904,
  16.65894442596359,
  53.996355074428614,
  38.270144423115774,
  193.64953557646763,
  48.35176764591674,
  72.09628287322178,
  12.654069756665395,
  38.8339941460539,
  12.595519237052462,
  122.17652320991567,
  32.78708495522601,
  6.246340876534137,
  21.35276860440632,
  103.65006610807802,
  18.37950320939907,
  65.30875245199718,
  31.777682584271428,
  11.9934698260187,
  77.92956892869215,
  36.49721293314736,
  71.239265142184,
  122.0610606632733,
  35.907863728312236,
  6.204573074331645,
  5.1524370709084275,
  55.804340067416746,
  20.779583193009863,
  19.75781849101845,
  17.571055407258257,
  52.55299288463493,
  50.14998879565123,
  35.65022850058403,
  85.41052275778333,
  56.2555650788846,
  57.23849621100226,
  9.040988273203602,
  5.393435592483144,
  17.884595057155213,
  57.504636863031024,
  69.78241276972112,
  11.629865209189063,
  60.89575977819283,
  30.90328360084,
  29.8121935448994,
  8.987816759050844,
  15.62090695596031,
  49.25848111915589,
  46.702426614008445,
  16.12379318123065,
  117.69004611416153,
  37.48535648654719,
  10.288642889780611,
  79.41222840551767,
  33.23594943738304,
  30.341678915904232,
  16.3865490921623,
  37.26010689837779,
  13.09119309603918,
  52.93914847664452,
  15.369360924522296,
  47.46616114642142,
  13.219065535564372,
  94.96129393064626,
  9.065508855893341,
  107.57836860958712,
  86.15077401788794,
  9.73163239302457,
  17.396314484373416,
  14.167533245697808,
  32.058409728159546,
  6.825834002638735,
  51.239449661593525,
  58.93740336626033,
  23.655523237322512,
  35.64049152404362,
  16.32943919920143,
  149.95813784939372,
  75.92476748778904,
  14.806688492056985,
  11.42172709769792,
  12.986917651710812,
  58.892394085191164,
  40.04194519137185,
  22.50355454347949,
  17.71787837931663,
  7.01627802521345,
  7.820297110945039,
  125.06472974834055,
  47.76856722361478,
  7.57037303484295,
  19.08473940808929,
  41.835368283748984,
  7.502814439472981,
  9.238889874409926,
  23.33341351072511,
  29.93371572184674,
  12.716583438589957,
  52.81800808660188,
  97.58754923380447,
  104.53354352277216,
  13.628983579329898,
  15.345176395034061,
  13.253673447037574,
  49.44141229447902,
  86.91245433298742,
  81.4023266407402,
  24.244521008207688,
  54.89592184591685,
  13.623843696980469,
  17.775820838582288,
  26.735643861058588,
  31.893504129394383,
  62.15543177348213,
  12.740465708216751,
  90.58736324751223,
  13.36420624141576,
  63.514639157202254,
  29.343776749341284,
  9.104132199686298,
  24.864152586811628,
  110.25589949515893,
  80.38591094106349,
  42.76975445077802,
  26.92601703749842,
  5.139137018060453,
  18.991205300705037,
  30.743551390884154,
  10.100506165017265,
  188.28953587846502,
  26.5083642002,
  7.457664106588232,
  85.2037456826129,
  60.294774072359175,
  37.186989379124206,
  7.2459628458718335,
  39.30215180743267,
  64.30323022022141,
  23.040552951353664,
  69.11603749455873,
  54.62361005109108,
  76.65400538607057,
  82.76549022260227,
  110.53504336507139,
  130.87130999124744,
  6.072808879926113,
  58.111851209136226,
  105.35752406971719,
  34.38036317554156,
  28.174263923786683,
  32.93411117660727,
  32.99055313746904,
  82.34374965142545,
  68.38648816663668,
  13.155785748994665,
  91.90680382008566,
  31.983977102081052,
  12.576389497423712,
  91.49274019929285,
  8.710000511982177,
  124.78674027659044,
  21.191379185396155,
  58.310930454959674,
  48.11106634708665,
  78.71454323431442,
  8.990780491054972,
  21.257192502627564,
  44.14320460383532,
  68.64433277544835,
  14.795440048944869,
  52.42760322465816,
  5.684284310337724,
  13.664912508120777,
  80.15294075129698,
  59.25054773553046,
  20.806311849834053,
  84.1034830302795,
  29.654837220431865,
  97.93362790216196,
  77.60043279710577,
  23.078720834992595,
  20.320338310963102,
  8.564982134652046,
  74.02373323676741,
  7.915265259407697,
  8.905475779391224,
  89.78927083052886,
  33.58490252823504,
  15.15277463915615,
  7.966275210505147,
  9.982631873396429,
  48.51825091615752,
  24.3351324715449,
  37.26790851629909,
  14.565572918921177,
  81.06153852630564,
  45.169041940064815,
  66.2566467434483,
  42.746235222276084,
  121.97616513811549,
  8.420166007439002,
  6.4248491233608345,
  36.864613302282265,
  140.99979581251912,
  97.79007539472119,
  120.14393530791497,
  16.9697638552666,
  71.91385799543811,
  5.241716525991229,
  7.983194593430584,
  9.351960656227808,
  61.69359033606149,
  27.516682857566728,
  51.06542241085249,
  66.79630259813428,
  14.634537460080129,
  90.59091303014901,
  61.085293292397424,
  11.593442606607278,
  7.218584255612209,
  7.310268912187566,
  11.084565318484033,
  20.0536050738007,
  30.29319862748777,
  54.913324375700626,
  60.3459176569474,
  70.39077208628738,
  36.229832084454166,
  29.14664949548409,
  73.83968586297199,
  25.61621942376874,
  16.977445889745475,
  15.658833202617416,
  35.54911488093509,
  23.770979608342763,
  35.20370765171865,
  53.877788441314955,
  70.65449756901506,
  46.24701145989481,
  7.252463696345999,
  7.467963201015802,
  24.255745544995037,
  54.7784164460617,
  119.0141470367416,
  70.16071427657803,
  28.32061750196084,
  12.370195205868335,
  75.99137627421038,
  7.096127943969323,
  70.97090832941112,
  38.81718835026143,
  63.86188918029771,
  87.56821433702686,
  27.735560469253016,
  12.712045199423539,
  79.77510687315201,
  35.228424982198234,
  34.70263446960864,
  71.83231182535239,
  55.88773199103073,
  36.908583024142715,
  30.894621351864586,
  81.93839162962121,
  36.04662453303813,
  50.3789414816418,
  39.620746577144274,
  5.881058026258764,
  31.25339257157878,
  74.79105159945803,
  34.99223608786868,
  67.50951682876378,
  19.231094599960265,
  15.791986493508077,
  57.56086448587226,
  26.157807521580434,
  100.69714851916748,
  79.72215670357868,
  7.5186349885285635,
  61.68177699135296,
  11.528407073010875,
  56.22668864926588,
  5.051980027559639,
  5.882305076554136,
  6.280895447041716,
  11.3654990365045,
  58.890535397130265,
  64.39763481596512,
  12.850678556387784,
  5.032611640767161,
  21.728739094415914,
  40.9365438830743,
  9.664339593854667,
  47.954509770718566,
  21.03970765732501,
  22.327828817994206,
  25.229019999770614,
  54.23928613320853,
  64.52049644226265,
  186.17141128534433,
  45.186552877895515,
  29.806307565536624,
  8.40744664974339,
  58.31608778728173,
  72.80306627412551,
  30.654815742712234,
  16.34176916151307,
  13.801535634857345,
  26.63197255242998,
  35.26472329384794,
  56.45371909211843])

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)