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

([8194426.5625,
  9007979.6875,
  9007980.788781097,
  9008015.567974329,
  9048940.625,
  9207851.5625,
  9222218.75,
  9224647.086230844,
  9225658.813289803,
  9300041.692603217,
  9321231.26734774,
  9381733.714107273,
  9439306.25,
  9752058.698369922,
  9752248.458814077,
  9949970.742946515,
  9951515.763801731,
  9952592.02227511,
  9954007.8125,
  9979901.443027167,
  9980293.0188882,
  9983306.297564307,
  9997508.916296797,
  10012672.370592983,
  10020371.607496593,
  10024610.753136452,
  10084867.1875,
  10155739.712103916,
  10160677.892555565,
  10160683.964612478,
  10163671.364678029,
  10172400.002939677,
  10176477.930598183,
  10185125.0,
  10186471.875,
  10189105.90446893,
  10189122.234375035,
  10190985.777817901,
  10192462.055379925,
  10208813.569724252,
  10211635.06068229,
  10211656.25,
  10212315.625,
  10212321.64251549,
  10213287.035756813,
  10213426.71581303,
  10213446.938607415,
  10217683.343096478,
  10218348.46623591,
  10218354.296484858,
  10222039.0625,
  10222967.181086302,
  10226560.9375,
  10226571.875,
  10226670.3125,
  10227224.577280194,
  10227726.510403039,
  10228258.188632103,
  10228268.620362328,
  10228625.0,
  10229437.487831319,
  10239182.052403504,
  10242891.673760792,
  10246754.6875,
  10246756.25,
  10247077.222481124,
  10248730.93154603,
  10255358.12880965,
  10256565.835860334,
  10257321.37043249,
  10329926.290634952,
  10331160.79814444,
  10331510.469889743,
  10335000.64140006,
  10336010.9375,
  10336084.375,
  10336988.970066717,
  10337535.9375,
  10337805.323144559,
  10337877.27304477,
  10346157.8125,
  10347713.281582655,
  10348215.678977918,
  10350388.883519586,
  10361455.961740501,
  10365614.13425318,
  10368029.910331147,
  10369880.18940451,
  10375499.863919131,
  10375662.830386726,
  10375785.9375,
  10379495.65879957,
  10381309.28501451,
  10383887.5,
  10383897.971979002,
  10389103.125,
  10396035.9375,
  10397480.374793146,
  10404267.99218849,
  10404319.239253357,
  10404496.32994227,
  10405163.724547233,
  10405798.4375,
  10406348.840215039,
  10406399.327740619,
  10406834.463060094,
  10407250.0,
  10407287.24484313,
  10408000.0,
  10409433.12570805,
  10409947.627913915,
  10410156.998947991,
  10410730.159676611,
  10410932.987844327,
  10413623.4375,
  10415888.317510428,
  10415892.856994668,
  10458247.417821053,
  10468743.296435447,
  10478046.554902913,
  10485881.80461816,
  10485990.625,
  10485994.170084119,
  10487242.61766642,
  10487470.904151736,
  10488639.0625,
  10490499.463557629,
  10493322.501128238,
  10502476.562339796,
  10519979.6875,
  10531096.71966659,
  10544617.1875,
  10567410.9375,
  10568839.0625,
  10573619.559459034,
  10573620.3125,
  10585144.468153704,
  10585417.1875,
  10591295.3125,
  10591563.28767183,
  10591740.625,
  10597931.27866446,
  10598174.753579237,
  10598704.747670537,
  10599034.910545133,
  10602705.906525807,
  10602914.066767057,
  10603428.125,
  10603832.8125,
  10603839.0625,
  10603950.0,
  10604618.37074226,
  10604981.25,
  10605607.8125,
  10605849.460008204,
  10606156.72160564,
  10607432.764858946,
  10608663.301518561,
  10609726.623548063,
  10609818.471214741,
  10610360.104812052,
  10610679.358630372,
  10610798.711701337,
  10610901.5625,
  10610974.520186258,
  10611125.0,
  10611155.318204867,
  10611221.875,
  10611276.5625,
  10612032.136345292,
  10612305.670112545,
  10613735.9375,
  10617029.429298695,
  10617707.902967792,
  10617709.013776207,
  10618025.859266253,
  10618030.808463143,
  10618389.540149285,
  10618733.433692362,
  10619320.3125,
  10619760.041487088,
  10620042.1875,
  10620482.8125,
  10620904.245523456,
  10620918.693427663,
  10621302.345684633,
  10621531.95838207,
  10621537.651053509,
  10621581.718308015,
  10621623.4375,
  10621746.875,
  10622053.540397298,
  10622203.00199093,
  10622490.625,
  10623146.221323734,
  10623150.676031128,
  10623317.306764124,
  10624888.492823606,
  10627462.901868073,
  10627750.572135545,
  10628219.80696334,
  10628414.0625,
  10629635.9375,
  10630165.294636628,
  10630193.75,
  10630830.3244337,
  10630904.6875,
  10631321.046748666,
  10631584.250007391,
  10631708.365214951,
  10631858.179277526,
  10632210.9375,
  10632470.710617732,
  10632704.195582036,
  10633071.875,
  10633225.662233654,
  10634573.4375,
  10646257.941417972,
  10647519.350374313,
  10649684.58193129,
  10650589.877306392,
  10650893.666895663,
  10650897.894471098,
  10651593.209459284,
  10652615.625,
  10654652.239153497,
  10655073.101043368,
  10655666.256528901,
  10656784.375,
  10661214.39199662,
  10663162.5,
  10663237.5,
  10663318.75,
  10664620.750345187,
  10664845.902922507,
  10665345.88729036,
  10665358.48814202,
  10665363.23798249,
  10665442.1875,
  10665827.522226505,
  10666101.678479921,
  10666102.162396852,
  10666323.4375,
  10666704.473170944,
  10666884.44053268,
  10666980.837394796,
  10666987.203726197,
  10667035.837948762,
  10667085.581458196,
  10667095.995219715,
  10667097.607153941,
  10667106.95328706,
  10667178.463881353,
  10667314.572177563,
  10667326.5625,
  10667641.809596898,
  10667670.3125,
  10667779.6875,
  10667800.0,
  10668045.147534115,
  10668198.4375,
  10668267.855773026,
  10668307.77531367,
  10668731.25827751,
  10668779.6875,
  10668841.683653843,
  10668934.200368706,
  10669082.761082401,
  10669284.402462967,
  10669330.688826755,
  10669388.870983042,
  10669628.955177063,
  10669840.287639402,
  10670051.508158518,
  10670923.4375,
  10670948.556673035,
  10671042.08900885,
  10671414.0625,
  10671583.594266461,
  10672114.976239199,
  10672806.25,
  10673075.0,
  10673099.02772723,
  10673718.75,
  10674424.723862963,
  10674846.215777814,
  10675068.75,
  10675311.464289881,
  10675768.29475451,
  10676151.398616692,
  10676241.521586902,
  10676295.3125,
  10676431.18793412,
  10677302.071385961,
  10677557.8125,
  10677801.5625,
  10677974.927239602,
  10677983.920668779,
  10677993.75,
  10678178.125,
  10678225.08399677,
  10678247.735952746,
  10678787.5,
  10678986.782303922,
  10679021.332575934,
  10679249.257874658,
  10679287.5,
  10679332.204993406,
  10679384.375,
  10679460.839701138,
  10680052.637895998,
  10684150.214178625,
  10684433.408242043,
  10684593.75,
  10685437.417640617,
  10687610.000355273,
  10689395.410292612,
  10693050.983483374,
  10695584.062223101,
  10698153.801184675,
  10698417.1875,
  10699159.149815282,
  10700074.411455559,
  10700093.75,
  10702228.243623562,
  10703377.601866776,
  10703675.930593522,
  10708579.382195435,
  10716253.125,
  10716897.029633995,
  10717238.620679403,
  10724123.4375,
  10724691.540826252,
  10725653.205679426,
  10726043.75,
  10726333.878310405,
  10726709.375,
  10727562.470892569,
  10728194.930847412,
  10729293.75,
  10729406.581059577,
  10729498.855642261,
  10729548.4375,
  10729829.6875,
  10729932.45185938,
  10729936.91494653,
  10729957.8125,
  10730022.343092557,
  10730109.706127003,
  10730210.594816752,
  10730469.954139352,
  10730619.829356158,
  10730629.63012114,
  10730650.136280507,
  10730949.616914203,
  10730994.304800713,
  10731255.691437561,
  10731546.875,
  10731749.483427815,
  10732935.495574178,
  10732960.720114041,
  10733059.605167136,
  10733125.395247223,
  10733255.906637331,
  10733288.18837584,
  10733372.61213823,
  10733376.352230968,
  10733454.661006652,
  10733465.734903507,
  10733587.961363181,
  10733687.472192083,
  10733815.473487189,
  10733826.035608636,
  10733967.959578259,
  10734204.53229983,
  10734271.890595283,
  10734535.740755893,
  10734557.010542015,
  10734563.989835024,
  10734698.35166462,
  10734950.0,
  10734983.761578597,
  10735127.649069538,
  10735565.459862305,
  10736184.375,
  10736194.876853177,
  10736328.382379571,
  10736607.981930826,
  10737511.672253339,
  10737551.5625,
  10737692.1875,
  10738206.25,
  10738556.606714938,
  10738656.509106094,
  10738981.398952805,
  10739422.992568307,
  10739457.8125,
  10739870.3125,
  10740257.652749237,
  10740350.76157476,
  10740422.662799142,
  10740749.041236116,
  10741681.78597816,
  10742285.9375,
  10742452.70571754,
  10744815.635573402,
  10744995.176785918,
  10746163.911209399,
  10748233.958437288,
  10748466.078470403,
  10748770.140936641,
  10749543.571246717,
  10749619.959137792,
  10749652.616260264,
  10749836.627291061,
  10750446.875,
  10750914.611101342,
  10751867.062037077,
  10752054.282003747,
  10752130.472185422,
  10752621.654243039,
  10753334.634885427,
  10754821.593598707,
  10756781.438443037,
  10757650.215437913,
  10757717.243570117,
  10757734.375,
  10761685.287727175,
  10761787.368747976,
  10762811.297617273,
  10763409.442316184,
  10784325.0,
  10785611.776125701,
  10790179.533992535,
  10800853.736558627,
  10801025.70014421,
  10801224.244481876,
  10801814.791647192,
  10804018.717560792,
  10807697.787808374,
  10807805.271445373,
  10810526.556379255,
  10821350.0,
  10823167.18693609,
  10830950.0,
  10831001.64710043,
  10831690.133419774,
  10833378.615434831,
  10849751.06873595,
  10849760.904473815,
  10869719.665163906,
  10883703.64894028,
  10917773.4375,
  10922868.096013507,
  10944681.426996509,
  10946166.54840526,
  10946276.5625,
  10946641.181216683,
  10946674.830690654,
  10949240.657753805,
  10949733.980566131,
  10950429.6875,
  10950707.265768236,
  10950790.623566462,
  10950809.375,
  10951138.409132151,
  10951328.635135477,
  10951362.5,
  10951438.888354348,
  10951744.967070037,
  10951809.05738127,
  10952185.9375,
  10952224.67525353,
  10952755.97506146,
  10952832.8125,
  10952913.567811498,
  10953097.931360053,
  10953344.109932166,
  10953513.671547595,
  10953585.9375,
  10953857.808962967,
  10953881.25,
  10953887.5,
  10954103.125,
  10954140.027678743,
  10954667.638055988,
  10955087.9499275,
  10955151.417637186,
  10955789.562155727,
  10956045.3125,
  10956146.116911612,
  10956183.064540097,
  10956238.241192807,
  10956475.552347831,
  10956545.3125,
  10956636.727818947,
  10959150.38932686,
  10959563.266369937,
  10960616.138274461,
  10961279.568197519,
  10961479.914005838,
  10961935.9375,
  10962015.312400375,
  10962158.352514409,
  10962309.831483819,
  10962422.80298365,
  10963235.9375,
  10963247.589060167,
  10963428.19313659,
  10964258.745098965,
  10965071.875,
  10970921.153333273,
  10970940.195213242,
  10980378.125,
  10987440.81026531,
  10987465.625,
  10998293.430567272,
  10998457.8125,
  11003134.375,
  11009473.512380058,
  11017965.625,
  11019489.040944865,
  11020406.245042948,
  11023040.420225363,
  11027090.56817726,
  11030069.338191263,
  11031638.820132403,
  11033573.4375,
  11033797.438723268,
  11034015.625,
  11036008.264716007,
  11036121.717488794,
  11047920.3125,
  11047936.234150114,
  11048089.304512007,
  11048518.544181656,
  11049419.806924181,
  11049706.25,
  11050223.4375,
  11050456.124238543,
  11050695.686735896,
  11050795.3125,
  11051185.579969808,
  11052187.5,
  11052459.375,
  11052498.4375,
  11052608.585341772,
  11052628.125,
  11052704.6875,
  11052807.794660637,
  11053052.010298248,
  11053057.8125,
  11053079.6875,
  11053420.3125,
  11053537.5,
  11053551.432155235,
  11053630.802711086,
  11053636.469264077,
  11053725.0,
  11053751.293748796,
  11053754.201753967,
  11054132.8125,
  11064935.9375,
  11066828.469645403,
  11069199.852247389,
  11070449.814942306,
  11070477.475705775,
  11072197.895028273,
  11073664.0625,
  11073696.928858122,
  11074027.797770614,
  11074724.74118652,
  11074772.514279012,
  11075249.253313374,
  11075279.075294733,
  11075294.666258208,
  11075624.921386193,
  11075840.625,
  11076190.74257594,
  11076355.490347927,
  11077287.09093583,
  11079744.392452262,
  11081800.0,
  11083407.741341138,
  11100244.93374215,
  11101520.3125,
  11102357.290362218,
  11107054.673681838,
  11107383.333312264,
  11108311.356817525,
  11109194.69750157,
  11111381.684249843,
  11118648.667844908,
  11120078.208567759,
  11120192.1875,
  11121236.025419911,
  11121564.746058064,
  11122366.65699521,
  11122967.1875,
  11123949.138291277,
  11126325.601857217,
  11126709.375,
  11127000.16201752,
  11127323.4375,
  11127749.04798044,
  11127944.628785584,
  11128120.3125,
  11128544.465129457,
  11128713.18315504,
  11130294.407833057,
  11131531.25,
  11131816.763923656,
  11132039.64816561,
  11134118.75,
  11134706.406815937,
  11135050.0,
  11140265.625,
  11147363.739480121,
  11147376.196770534,
  11148091.66939472,
  11150190.415576423,
  11150499.535508212,
  11167139.795033026,
  11167146.733910596,
  11167181.207791744,
  11168387.958579028,
  11170831.25,
  11173257.8125,
  11173660.238696704,
  11207580.720249677,
  11208721.875,
  11217254.257382715,
  11217560.900176646,
  11219270.3125,
  11256710.771725755,
  11259235.662703464,
  11259769.066427844,
  11260432.158499919,
  11264139.104121426,
  11264947.645925201,
  11265092.216840269,
  11265641.60255622,
  11265812.5,
  11266009.312730523,
  11266062.158604996,
  11266139.0625,
  11266493.75,
  11266506.211218331,
  11266541.131823104,
  11266789.0625,
  11267146.533449108,
  11267248.4375,
  11267360.94005182,
  11267464.273117088,
  11267889.0625,
  11268016.517799657,
  11268050.0,
  11268100.470255062,
  11268243.516179841,
  11268710.301103912,
  11268909.814778535,
  11273924.728893256,
  11285058.507062903,
  11301622.869773796,
  11301678.125,
  11305879.008617057,
  11306688.891365174,
  11317179.6875,
  11317521.875,
  11325264.725981757,
  11327557.504718628,
  11335565.388846623,
  11335606.46004931,
  11343465.810245508,
  11343683.688706458,
  11358159.375,
  11361792.1875,
  11364869.84936274,
  11386646.995896783,
  11391195.3125,
  11613991.801342947,
  11644575.721773047,
  11825167.704729134,
  11826532.618146012,
  12359810.9375,
  12405788.561093709,
  12412255.992851308,
  12766435.9375,
  12775486.38418735,
  12776460.666151168,
  12777679.6875,
  12779045.3125,
  12780132.562330503,
  12781445.3125,
  12789219.27077604,
  12789290.3691226,
  12791514.493317692,
  12791669.926699363,
  12793248.4375,
  12793939.012875333,
  12797828.125,
  12800243.972104227,
  12800286.710933585,
  12804194.611044167,
  12833462.757222446,
  12847538.959554236,
  12852471.875,
  13076147.750639915,
  13137595.408848872,
  13164774.863675997,
  13174558.425561663,
  13203401.5625,
  13203407.67259257,
  13206251.808167202,
  13209997.076022692,
  13210596.875,
  14081595.109291937,
  14099023.430746213,
  14138059.809210049,
  14175562.858426778],
 [39.50316505699135,
  67.98118973901671,
  10.475974835742623,
  9.283645768105531,
  43.19733067911364,
  46.72084718618287,
  112.76887399498642,
  14.338767177944092,
  97.918332693213,
  10.433252092058375,
  11.154931133489224,
  5.578499140480954,
  171.0161099825164,
  8.818071195696492,
  12.843343636464379,
  16.882440591201906,
  11.351822441213697,
  19.4165965061271,
  72.28675643873743,
  14.504883258177252,
  5.317478152842318,
  15.125955356278041,
  21.915192867842556,
  28.43948359473771,
  7.052864853486413,
  52.67705388113174,
  110.18496295559889,
  12.321574972962075,
  54.28176176914436,
  33.524651031298205,
  13.29786520492926,
  7.122224623379808,
  25.200326373046824,
  100.9278902654627,
  53.47690291703544,
  45.68603105929498,
  28.34435113107704,
  63.55581010773996,
  22.06392500416911,
  59.25098737433268,
  137.2246871549226,
  69.04375388938725,
  36.47629596348966,
  37.628643560581224,
  55.25368581672657,
  9.464492020581726,
  15.189623174582835,
  12.07243313710575,
  112.73610992701094,
  22.65383545202508,
  78.30760580888762,
  74.52571466998572,
  109.36107687799321,
  49.25151524499393,
  76.76951799222392,
  97.38296466841138,
  184.67112861261216,
  21.820369773355335,
  14.469388563229382,
  56.27002224380116,
  40.07203909668597,
  5.370190083145916,
  8.032385437302798,
  42.95557843724215,
  46.03756987483408,
  18.239280315574085,
  23.517845529202095,
  11.660938696150717,
  141.82845458917473,
  30.12287385241698,
  100.12519723385515,
  10.916654295819153,
  20.525468473093618,
  18.118620291298082,
  45.249649132104395,
  81.52535769530303,
  164.840737797297,
  47.690690827365245,
  20.243315694287496,
  46.82265063992616,
  30.640038596597066,
  6.930264112494245,
  6.495586163170042,
  52.35241462809125,
  75.29059316127328,
  11.241744909791224,
  5.170583868778721,
  20.8779091691581,
  36.532671435122296,
  59.1264381076774,
  65.51480269991409,
  100.28263274581711,
  5.163271165077303,
  100.08339375069251,
  28.530200005352768,
  59.42226366510195,
  169.16689145982713,
  18.420521638537195,
  124.48521262540953,
  5.049886296354172,
  21.05510015536838,
  71.58755458122793,
  135.62392647896442,
  11.821480264401828,
  156.56124919286145,
  11.268197800474285,
  84.19985025652363,
  32.06203287820941,
  40.72358733787902,
  21.890162543640233,
  85.05769072728724,
  7.8986112063784795,
  43.61816050384583,
  21.148538191844224,
  120.1479857856996,
  12.701678672839142,
  18.30187995691229,
  9.836832189598512,
  21.969172175576638,
  87.53097201744308,
  109.78303601419351,
  30.754174396294367,
  16.827406308327546,
  33.69657480567423,
  46.84345867585989,
  54.49256816428344,
  21.13711115856213,
  63.914465018185616,
  5.098453072358372,
  54.14416381720697,
  54.155580063974526,
  43.50012159871536,
  38.40719992467651,
  59.46819460795218,
  36.12469615871113,
  64.77537698628942,
  16.540021706797532,
  118.93985599089858,
  119.37355715310139,
  12.631508817806465,
  61.52752982278932,
  11.172808747355564,
  167.78914028142808,
  8.318654898119922,
  12.082482934963709,
  14.320003277902503,
  13.189125169058485,
  64.20040837410554,
  60.743866726123784,
  148.0040344779177,
  63.506012846898585,
  33.037747486236604,
  83.61165944756007,
  66.68941854937256,
  5.507434673401909,
  12.241326015228964,
  92.07402479353142,
  19.232676850345797,
  183.04346327815165,
  69.11463909836347,
  18.712274366461852,
  77.40474405681056,
  5.834551274349883,
  63.39302136388497,
  59.17137936773298,
  71.56796665841725,
  12.280943448568983,
  45.72530776167984,
  46.31745114280217,
  34.630257348748444,
  29.896252614242968,
  57.603192790750796,
  60.88308765824151,
  14.543146064676076,
  7.477613189232802,
  31.339899846522307,
  10.501221232944795,
  45.89334633945878,
  12.380655616575535,
  32.80777126786844,
  39.94217339735157,
  58.657165281011466,
  56.15204557435317,
  22.09212438479612,
  62.22700927973831,
  38.225569100157756,
  73.45287949372364,
  15.55030579217439,
  8.550790715402066,
  62.156793227722886,
  48.694236095658574,
  8.519181382170569,
  127.49896016940642,
  74.35725229200243,
  10.239584531634653,
  79.31756973372542,
  27.187897341703405,
  39.676251016094014,
  25.270125138171736,
  6.580681302260346,
  21.438046325948587,
  68.52869301077271,
  56.32787557445083,
  6.699296911303664,
  84.81195249402639,
  6.156937611144566,
  79.34649857951442,
  7.328835452646465,
  7.558647047478901,
  9.84699546056425,
  122.46400322140755,
  90.37138020116767,
  123.23393382180436,
  21.5741035135771,
  31.973075449484124,
  211.05273636199007,
  52.3756682219781,
  54.37979379497091,
  181.80183785015723,
  15.793606803403478,
  12.695270854836156,
  27.480866527807308,
  21.074340742478125,
  54.16684174716094,
  88.261347702655,
  107.61382329947793,
  23.29596714529758,
  18.055367811545647,
  117.366133581356,
  86.57085999908834,
  58.5329977855276,
  118.53085534949352,
  114.28349954453392,
  16.922032320229476,
  17.798091778896588,
  20.458006591773838,
  9.976475147979615,
  18.769727652286633,
  99.9455425558037,
  43.89973990013543,
  49.80085948127204,
  5.74603120059032,
  38.130784164585314,
  105.11395066089355,
  11.231707206535946,
  9.005299157995132,
  13.92092632862266,
  14.575173908021776,
  167.46296337851243,
  65.60731307079084,
  52.59176147794144,
  16.77276589992591,
  16.41230238264284,
  32.77197772102608,
  39.08331626137966,
  22.80415572628247,
  141.91655220764213,
  91.92467701756958,
  121.07781002347892,
  57.666724434261084,
  31.38139508354997,
  13.836561093882093,
  27.8096839030851,
  10.374361731516549,
  43.2046311574842,
  82.70019922836426,
  6.729136485409227,
  21.365479383067,
  7.128689597436717,
  58.103539549638526,
  14.258546238283651,
  5.498563743302872,
  23.284961250720276,
  17.010685648832077,
  30.370300503328743,
  27.873587333053564,
  39.794891099579914,
  317.264460074431,
  42.869925292417975,
  20.12726928224934,
  31.087947989101828,
  92.85253580039058,
  17.796892604748905,
  51.86295698276669,
  47.60891241664507,
  20.039357856911572,
  64.65907003187658,
  17.418322391910742,
  8.974009552725779,
  41.25246632208623,
  19.9836427031011,
  87.60460094520714,
  205.70680637198802,
  8.575357914864341,
  91.70076528392454,
  39.94930390076782,
  38.724564142595966,
  6.211454123512151,
  34.50196643718585,
  66.19265898201763,
  43.368771562997715,
  7.367787491327854,
  70.64014478191929,
  16.748967683683865,
  21.080085331452803,
  37.02097284894502,
  93.29296954000951,
  56.57272601070615,
  63.11233906260891,
  26.7299243347594,
  21.62271658412902,
  20.60244840400898,
  12.206855918780509,
  158.21937990948868,
  43.81581669546195,
  5.9957804594264195,
  14.49804175605524,
  29.045476364368568,
  42.95685976086347,
  12.3794300179355,
  60.002422912892186,
  24.75588069153058,
  68.96586736966843,
  46.49526127280067,
  8.350319691507895,
  19.31364446519322,
  32.58273834981632,
  51.60347229964887,
  30.842982512175833,
  21.72821663263724,
  13.488260944579789,
  68.16291121039966,
  12.94152362513063,
  24.13371368218136,
  52.7024518326302,
  8.94399700941504,
  44.07779010452224,
  7.113821162289917,
  22.760919223852838,
  183.68952898995354,
  11.687854875288638,
  54.21955842375287,
  38.989546208267456,
  85.80772985010569,
  49.79673443711015,
  117.93790648148806,
  29.66464002724236,
  9.761263370426112,
  280.14832419444963,
  23.22041559779727,
  33.95434962496035,
  10.448093811092093,
  13.438438024661767,
  5.133425479576933,
  16.710918313127475,
  34.18430914018503,
  13.168468712223998,
  35.18378824343972,
  6.409052293536936,
  21.975517311805895,
  24.866679081739058,
  54.986971361983436,
  11.807327722191385,
  32.16828886480468,
  12.406246384690967,
  138.62882808050097,
  53.38215312488828,
  61.69675895424119,
  7.643477533354556,
  44.01396304080855,
  27.98612411810494,
  6.248243738001401,
  21.228974309955674,
  30.272422260672542,
  25.795128827537994,
  89.55463739983735,
  17.22852174774867,
  49.03696515142512,
  11.015995072244392,
  34.53788227779839,
  93.07063806943887,
  143.12699157188348,
  188.0788107248651,
  19.09171758762296,
  94.594987882581,
  217.69143758541898,
  48.021933137036754,
  42.433374335887734,
  5.427325995685693,
  69.10741624106771,
  58.44649630496571,
  77.27462336593153,
  79.92324475179814,
  6.115464905356458,
  15.541360360602457,
  119.83512328313321,
  245.2513389074015,
  54.95601040747751,
  39.79674531032294,
  36.33063957014039,
  19.13021494192061,
  17.880355775712836,
  7.269851314316384,
  36.67046067112894,
  196.83566563410648,
  11.144274027383975,
  121.45688832239892,
  115.09905926256451,
  42.95695223082459,
  17.006472493453117,
  55.235415287515856,
  36.350266621736125,
  23.092233267582078,
  75.94126158127837,
  79.54774983702579,
  67.52821268696961,
  107.95182017905461,
  26.269449616442465,
  86.69673658557966,
  8.750238064077863,
  53.143423608527215,
  18.534601510029795,
  29.47390563761767,
  86.54984289177627,
  15.83157506605092,
  5.292822371868798,
  64.86804594195122,
  45.62181094600685,
  26.175474528346665,
  35.136595532772105,
  105.84798670794297,
  42.31068596256322,
  12.752055955854289,
  35.01353804943281,
  29.912663654740633,
  148.96922531686255,
  21.397040671964866,
  5.754789353631851,
  27.898663322647938,
  6.960535101717085,
  17.758430239961733,
  208.73987030044356,
  121.08638988630685,
  66.87419870337315,
  28.533081436975788,
  81.88307268943217,
  10.437432965647504,
  31.655140345207386,
  21.555533522184323,
  26.339029500836943,
  29.830742328708673,
  40.917244378492335,
  51.328837124411564,
  100.72647055894237,
  249.71683219370942,
  45.29884271245434,
  85.72728235851778,
  14.296925893606659,
  33.436782034711484,
  89.22098673132109,
  93.00418892387427,
  106.04860353705973,
  24.35947589657486,
  240.5133756211648,
  55.29936135366497,
  12.924757640611173,
  5.558362076355664,
  32.09336785343787,
  72.61274505494005,
  39.933268901371434,
  10.665068051913394,
  37.44563436079499,
  13.833096750077926,
  24.076234157346146,
  51.36190317645799,
  52.06298262482526,
  34.962808104015224,
  6.126346234412228,
  10.56706923952123,
  69.42156034995841,
  5.8137186022315595,
  55.870555121813126,
  31.267887602902707,
  77.45468142604233,
  55.1674483705808,
  19.488894849414848,
  5.140235838273593,
  14.416684870112073,
  6.568376930671347,
  61.745195099021714,
  19.104353814473367,
  9.554015103434526,
  154.364746163956,
  180.50940237057165,
  43.88335866126284,
  67.15078573173517,
  56.702802705225466,
  6.171593579121499,
  133.64086025580562,
  26.370345058017897,
  15.881353712258855,
  44.291779230488956,
  23.59112212455028,
  185.12404191948178,
  11.924984305919258,
  50.81014347660295,
  29.608125076867246,
  8.704095647400791,
  63.85135895188379,
  10.267796839308943,
  239.83471601720842,
  10.160366551818761,
  13.522546916263728,
  117.65038163483892,
  15.69802687563399,
  54.29769705298167,
  23.505649161650165,
  82.71773438798344,
  112.15022448191117,
  11.246801638266403,
  38.897067836405846,
  28.105241657299036,
  77.08391804322174,
  78.32327295153378,
  140.5529760462183,
  17.784402725463,
  56.88055790759188,
  49.53040158816948,
  5.361616764259064,
  95.5865679863165,
  11.916664361816641,
  31.277784732601013,
  76.3292078122316,
  43.04791942118595,
  69.79503048522363,
  205.1453442959711,
  21.437658215864364,
  69.69174512797319,
  167.59064327992328,
  10.936663265682077,
  90.32961923125609,
  59.6360171846313,
  43.36172722626316,
  80.81973251501296,
  55.502445334401315,
  65.24663440316296,
  12.597069867649909,
  122.7851219197621,
  36.38648019411924,
  92.96538678721858,
  16.97568893827716,
  94.20382537397045,
  46.014881344188026,
  100.87989361553686,
  48.45975791076177,
  10.929903575409403,
  6.419188511802526,
  74.73375876292599,
  33.53142057881886,
  14.00738332081079,
  102.88951099268658,
  104.57337625198284,
  178.82023307639827,
  9.640875386187458,
  14.406743329393132,
  25.33909073688116,
  12.934952611126516,
  29.659763248731572,
  67.68060034522952,
  8.315265241217455,
  6.699890069788223,
  61.14502615351799,
  153.16604304872877,
  10.13174754877868,
  70.67247077433476,
  48.6727632232121,
  57.693098935846436,
  87.8705213561676,
  14.74996893220896,
  8.762262413819911,
  266.8341285953496,
  53.94273713016208,
  78.2918588799144,
  14.677048703049477,
  65.10229005857747,
  48.482791836267424,
  92.51818144712378,
  21.74348126099796,
  8.582633039813405,
  16.823354130900178,
  13.021753325625962,
  120.57941575875051,
  6.112385156712894,
  15.269893477015982,
  49.53842392320138,
  84.99749220461472,
  65.111582861139,
  21.19410190848156,
  83.18432435823185,
  65.96495571799105,
  5.746163229296661,
  36.53985710851037,
  15.597226634384171,
  30.121563948284003,
  83.80534800357908,
  83.05571949720654,
  105.40613805787011,
  17.63034772766711,
  19.503482933788554,
  172.20884847222757,
  172.76590816379965,
  62.66060381272337,
  5.745073480472719,
  46.26407272278802,
  184.27969168430545,
  106.1502881619002,
  77.12238463178318,
  23.461081744352935,
  48.0592746327572,
  8.928640150612145,
  122.39405626496612,
  21.854780752492147,
  12.50489211144153,
  15.511759110799773,
  87.41363178193195,
  9.751791601320102,
  51.72895020235779,
  57.09813386267075,
  5.269978964581719,
  61.99640944060233,
  81.24362552023702,
  38.160456006768655,
  14.801133306894108,
  34.71918848391615,
  61.126036860280095,
  60.02210479015393,
  74.7360472403996,
  20.08311870618091,
  15.102559950606805,
  36.2013445158081,
  228.8605933688175,
  7.999712780848196,
  51.6585053053881,
  85.09709289105838,
  106.94011337740989,
  111.70531486061282,
  46.83885217464058,
  312.4528765779274,
  8.574102087018632,
  297.0373019280596,
  23.23411036707691,
  35.65677939817195,
  14.951244695236474,
  6.110278859062512,
  84.16856601252927,
  8.822793188947056,
  96.46873951833324,
  6.1384060074468305,
  22.520152641449176,
  29.356014436532597,
  15.333481833557244,
  27.378530491661877,
  12.419924308306046,
  13.144309721098136,
  42.666640096519295,
  27.359965317166278,
  25.54930992101975,
  52.456059678950716,
  110.96406293338795,
  18.262249301253334,
  56.926206554031815,
  24.602055552335177,
  33.08469173311845,
  11.430152314651679,
  69.76586902768939,
  90.53347811252846,
  79.10642867431459,
  7.004805931076852,
  88.26067278520506,
  46.553460625703856,
  195.13321801985725,
  12.586149015353122,
  12.04541644021504,
  6.725315838947938,
  39.16008383245138,
  26.882421123351886,
  53.01498674088756,
  61.55531170779169,
  18.506147656527652,
  24.123752240435177,
  36.220924431059395,
  77.00875205535822,
  24.408833387887526,
  74.55698943557425,
  42.030491949670974,
  5.6906645331502475,
  16.88442088971928,
  17.72218341082624,
  149.495764297482,
  27.579697480683027,
  94.41426366208414,
  6.5224809174786085,
  22.789855823658115,
  59.75870854853966,
  16.021789474290898,
  62.13039931733603,
  32.372088360356216,
  93.0881293434344,
  93.50188585235915,
  31.13445806268396,
  7.292504726327672,
  73.75925724298725,
  14.433778254666748,
  5.83238393698658,
  13.82565030106506,
  80.39383503016572,
  6.72360460858169,
  60.33725145198902,
  11.8817623427884,
  36.811620496808466])

Energy spectrum

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

([8194426.5625,
  9007979.6875,
  9007980.788781097,
  9008015.567974329,
  9048940.625,
  9207851.5625,
  9222218.75,
  9224647.086230844,
  9225658.813289803,
  9300041.692603217,
  9321231.26734774,
  9381733.714107273,
  9439306.25,
  9752058.698369922,
  9752248.458814077,
  9949970.742946515,
  9951515.763801731,
  9952592.02227511,
  9954007.8125,
  9979901.443027167,
  9980293.0188882,
  9983306.297564307,
  9997508.916296797,
  10012672.370592983,
  10020371.607496593,
  10024610.753136452,
  10084867.1875,
  10155739.712103916,
  10160677.892555565,
  10160683.964612478,
  10163671.364678029,
  10172400.002939677,
  10176477.930598183,
  10185125.0,
  10186471.875,
  10189105.90446893,
  10189122.234375035,
  10190985.777817901,
  10192462.055379925,
  10208813.569724252,
  10211635.06068229,
  10211656.25,
  10212315.625,
  10212321.64251549,
  10213287.035756813,
  10213426.71581303,
  10213446.938607415,
  10217683.343096478,
  10218348.46623591,
  10218354.296484858,
  10222039.0625,
  10222967.181086302,
  10226560.9375,
  10226571.875,
  10226670.3125,
  10227224.577280194,
  10227726.510403039,
  10228258.188632103,
  10228268.620362328,
  10228625.0,
  10229437.487831319,
  10239182.052403504,
  10242891.673760792,
  10246754.6875,
  10246756.25,
  10247077.222481124,
  10248730.93154603,
  10255358.12880965,
  10256565.835860334,
  10257321.37043249,
  10329926.290634952,
  10331160.79814444,
  10331510.469889743,
  10335000.64140006,
  10336010.9375,
  10336084.375,
  10336988.970066717,
  10337535.9375,
  10337805.323144559,
  10337877.27304477,
  10346157.8125,
  10347713.281582655,
  10348215.678977918,
  10350388.883519586,
  10361455.961740501,
  10365614.13425318,
  10368029.910331147,
  10369880.18940451,
  10375499.863919131,
  10375662.830386726,
  10375785.9375,
  10379495.65879957,
  10381309.28501451,
  10383887.5,
  10383897.971979002,
  10389103.125,
  10396035.9375,
  10397480.374793146,
  10404267.99218849,
  10404319.239253357,
  10404496.32994227,
  10405163.724547233,
  10405798.4375,
  10406348.840215039,
  10406399.327740619,
  10406834.463060094,
  10407250.0,
  10407287.24484313,
  10408000.0,
  10409433.12570805,
  10409947.627913915,
  10410156.998947991,
  10410730.159676611,
  10410932.987844327,
  10413623.4375,
  10415888.317510428,
  10415892.856994668,
  10458247.417821053,
  10468743.296435447,
  10478046.554902913,
  10485881.80461816,
  10485990.625,
  10485994.170084119,
  10487242.61766642,
  10487470.904151736,
  10488639.0625,
  10490499.463557629,
  10493322.501128238,
  10502476.562339796,
  10519979.6875,
  10531096.71966659,
  10544617.1875,
  10567410.9375,
  10568839.0625,
  10573619.559459034,
  10573620.3125,
  10585144.468153704,
  10585417.1875,
  10591295.3125,
  10591563.28767183,
  10591740.625,
  10597931.27866446,
  10598174.753579237,
  10598704.747670537,
  10599034.910545133,
  10602705.906525807,
  10602914.066767057,
  10603428.125,
  10603832.8125,
  10603839.0625,
  10603950.0,
  10604618.37074226,
  10604981.25,
  10605607.8125,
  10605849.460008204,
  10606156.72160564,
  10607432.764858946,
  10608663.301518561,
  10609726.623548063,
  10609818.471214741,
  10610360.104812052,
  10610679.358630372,
  10610798.711701337,
  10610901.5625,
  10610974.520186258,
  10611125.0,
  10611155.318204867,
  10611221.875,
  10611276.5625,
  10612032.136345292,
  10612305.670112545,
  10613735.9375,
  10617029.429298695,
  10617707.902967792,
  10617709.013776207,
  10618025.859266253,
  10618030.808463143,
  10618389.540149285,
  10618733.433692362,
  10619320.3125,
  10619760.041487088,
  10620042.1875,
  10620482.8125,
  10620904.245523456,
  10620918.693427663,
  10621302.345684633,
  10621531.95838207,
  10621537.651053509,
  10621581.718308015,
  10621623.4375,
  10621746.875,
  10622053.540397298,
  10622203.00199093,
  10622490.625,
  10623146.221323734,
  10623150.676031128,
  10623317.306764124,
  10624888.492823606,
  10627462.901868073,
  10627750.572135545,
  10628219.80696334,
  10628414.0625,
  10629635.9375,
  10630165.294636628,
  10630193.75,
  10630830.3244337,
  10630904.6875,
  10631321.046748666,
  10631584.250007391,
  10631708.365214951,
  10631858.179277526,
  10632210.9375,
  10632470.710617732,
  10632704.195582036,
  10633071.875,
  10633225.662233654,
  10634573.4375,
  10646257.941417972,
  10647519.350374313,
  10649684.58193129,
  10650589.877306392,
  10650893.666895663,
  10650897.894471098,
  10651593.209459284,
  10652615.625,
  10654652.239153497,
  10655073.101043368,
  10655666.256528901,
  10656784.375,
  10661214.39199662,
  10663162.5,
  10663237.5,
  10663318.75,
  10664620.750345187,
  10664845.902922507,
  10665345.88729036,
  10665358.48814202,
  10665363.23798249,
  10665442.1875,
  10665827.522226505,
  10666101.678479921,
  10666102.162396852,
  10666323.4375,
  10666704.473170944,
  10666884.44053268,
  10666980.837394796,
  10666987.203726197,
  10667035.837948762,
  10667085.581458196,
  10667095.995219715,
  10667097.607153941,
  10667106.95328706,
  10667178.463881353,
  10667314.572177563,
  10667326.5625,
  10667641.809596898,
  10667670.3125,
  10667779.6875,
  10667800.0,
  10668045.147534115,
  10668198.4375,
  10668267.855773026,
  10668307.77531367,
  10668731.25827751,
  10668779.6875,
  10668841.683653843,
  10668934.200368706,
  10669082.761082401,
  10669284.402462967,
  10669330.688826755,
  10669388.870983042,
  10669628.955177063,
  10669840.287639402,
  10670051.508158518,
  10670923.4375,
  10670948.556673035,
  10671042.08900885,
  10671414.0625,
  10671583.594266461,
  10672114.976239199,
  10672806.25,
  10673075.0,
  10673099.02772723,
  10673718.75,
  10674424.723862963,
  10674846.215777814,
  10675068.75,
  10675311.464289881,
  10675768.29475451,
  10676151.398616692,
  10676241.521586902,
  10676295.3125,
  10676431.18793412,
  10677302.071385961,
  10677557.8125,
  10677801.5625,
  10677974.927239602,
  10677983.920668779,
  10677993.75,
  10678178.125,
  10678225.08399677,
  10678247.735952746,
  10678787.5,
  10678986.782303922,
  10679021.332575934,
  10679249.257874658,
  10679287.5,
  10679332.204993406,
  10679384.375,
  10679460.839701138,
  10680052.637895998,
  10684150.214178625,
  10684433.408242043,
  10684593.75,
  10685437.417640617,
  10687610.000355273,
  10689395.410292612,
  10693050.983483374,
  10695584.062223101,
  10698153.801184675,
  10698417.1875,
  10699159.149815282,
  10700074.411455559,
  10700093.75,
  10702228.243623562,
  10703377.601866776,
  10703675.930593522,
  10708579.382195435,
  10716253.125,
  10716897.029633995,
  10717238.620679403,
  10724123.4375,
  10724691.540826252,
  10725653.205679426,
  10726043.75,
  10726333.878310405,
  10726709.375,
  10727562.470892569,
  10728194.930847412,
  10729293.75,
  10729406.581059577,
  10729498.855642261,
  10729548.4375,
  10729829.6875,
  10729932.45185938,
  10729936.91494653,
  10729957.8125,
  10730022.343092557,
  10730109.706127003,
  10730210.594816752,
  10730469.954139352,
  10730619.829356158,
  10730629.63012114,
  10730650.136280507,
  10730949.616914203,
  10730994.304800713,
  10731255.691437561,
  10731546.875,
  10731749.483427815,
  10732935.495574178,
  10732960.720114041,
  10733059.605167136,
  10733125.395247223,
  10733255.906637331,
  10733288.18837584,
  10733372.61213823,
  10733376.352230968,
  10733454.661006652,
  10733465.734903507,
  10733587.961363181,
  10733687.472192083,
  10733815.473487189,
  10733826.035608636,
  10733967.959578259,
  10734204.53229983,
  10734271.890595283,
  10734535.740755893,
  10734557.010542015,
  10734563.989835024,
  10734698.35166462,
  10734950.0,
  10734983.761578597,
  10735127.649069538,
  10735565.459862305,
  10736184.375,
  10736194.876853177,
  10736328.382379571,
  10736607.981930826,
  10737511.672253339,
  10737551.5625,
  10737692.1875,
  10738206.25,
  10738556.606714938,
  10738656.509106094,
  10738981.398952805,
  10739422.992568307,
  10739457.8125,
  10739870.3125,
  10740257.652749237,
  10740350.76157476,
  10740422.662799142,
  10740749.041236116,
  10741681.78597816,
  10742285.9375,
  10742452.70571754,
  10744815.635573402,
  10744995.176785918,
  10746163.911209399,
  10748233.958437288,
  10748466.078470403,
  10748770.140936641,
  10749543.571246717,
  10749619.959137792,
  10749652.616260264,
  10749836.627291061,
  10750446.875,
  10750914.611101342,
  10751867.062037077,
  10752054.282003747,
  10752130.472185422,
  10752621.654243039,
  10753334.634885427,
  10754821.593598707,
  10756781.438443037,
  10757650.215437913,
  10757717.243570117,
  10757734.375,
  10761685.287727175,
  10761787.368747976,
  10762811.297617273,
  10763409.442316184,
  10784325.0,
  10785611.776125701,
  10790179.533992535,
  10800853.736558627,
  10801025.70014421,
  10801224.244481876,
  10801814.791647192,
  10804018.717560792,
  10807697.787808374,
  10807805.271445373,
  10810526.556379255,
  10821350.0,
  10823167.18693609,
  10830950.0,
  10831001.64710043,
  10831690.133419774,
  10833378.615434831,
  10849751.06873595,
  10849760.904473815,
  10869719.665163906,
  10883703.64894028,
  10917773.4375,
  10922868.096013507,
  10944681.426996509,
  10946166.54840526,
  10946276.5625,
  10946641.181216683,
  10946674.830690654,
  10949240.657753805,
  10949733.980566131,
  10950429.6875,
  10950707.265768236,
  10950790.623566462,
  10950809.375,
  10951138.409132151,
  10951328.635135477,
  10951362.5,
  10951438.888354348,
  10951744.967070037,
  10951809.05738127,
  10952185.9375,
  10952224.67525353,
  10952755.97506146,
  10952832.8125,
  10952913.567811498,
  10953097.931360053,
  10953344.109932166,
  10953513.671547595,
  10953585.9375,
  10953857.808962967,
  10953881.25,
  10953887.5,
  10954103.125,
  10954140.027678743,
  10954667.638055988,
  10955087.9499275,
  10955151.417637186,
  10955789.562155727,
  10956045.3125,
  10956146.116911612,
  10956183.064540097,
  10956238.241192807,
  10956475.552347831,
  10956545.3125,
  10956636.727818947,
  10959150.38932686,
  10959563.266369937,
  10960616.138274461,
  10961279.568197519,
  10961479.914005838,
  10961935.9375,
  10962015.312400375,
  10962158.352514409,
  10962309.831483819,
  10962422.80298365,
  10963235.9375,
  10963247.589060167,
  10963428.19313659,
  10964258.745098965,
  10965071.875,
  10970921.153333273,
  10970940.195213242,
  10980378.125,
  10987440.81026531,
  10987465.625,
  10998293.430567272,
  10998457.8125,
  11003134.375,
  11009473.512380058,
  11017965.625,
  11019489.040944865,
  11020406.245042948,
  11023040.420225363,
  11027090.56817726,
  11030069.338191263,
  11031638.820132403,
  11033573.4375,
  11033797.438723268,
  11034015.625,
  11036008.264716007,
  11036121.717488794,
  11047920.3125,
  11047936.234150114,
  11048089.304512007,
  11048518.544181656,
  11049419.806924181,
  11049706.25,
  11050223.4375,
  11050456.124238543,
  11050695.686735896,
  11050795.3125,
  11051185.579969808,
  11052187.5,
  11052459.375,
  11052498.4375,
  11052608.585341772,
  11052628.125,
  11052704.6875,
  11052807.794660637,
  11053052.010298248,
  11053057.8125,
  11053079.6875,
  11053420.3125,
  11053537.5,
  11053551.432155235,
  11053630.802711086,
  11053636.469264077,
  11053725.0,
  11053751.293748796,
  11053754.201753967,
  11054132.8125,
  11064935.9375,
  11066828.469645403,
  11069199.852247389,
  11070449.814942306,
  11070477.475705775,
  11072197.895028273,
  11073664.0625,
  11073696.928858122,
  11074027.797770614,
  11074724.74118652,
  11074772.514279012,
  11075249.253313374,
  11075279.075294733,
  11075294.666258208,
  11075624.921386193,
  11075840.625,
  11076190.74257594,
  11076355.490347927,
  11077287.09093583,
  11079744.392452262,
  11081800.0,
  11083407.741341138,
  11100244.93374215,
  11101520.3125,
  11102357.290362218,
  11107054.673681838,
  11107383.333312264,
  11108311.356817525,
  11109194.69750157,
  11111381.684249843,
  11118648.667844908,
  11120078.208567759,
  11120192.1875,
  11121236.025419911,
  11121564.746058064,
  11122366.65699521,
  11122967.1875,
  11123949.138291277,
  11126325.601857217,
  11126709.375,
  11127000.16201752,
  11127323.4375,
  11127749.04798044,
  11127944.628785584,
  11128120.3125,
  11128544.465129457,
  11128713.18315504,
  11130294.407833057,
  11131531.25,
  11131816.763923656,
  11132039.64816561,
  11134118.75,
  11134706.406815937,
  11135050.0,
  11140265.625,
  11147363.739480121,
  11147376.196770534,
  11148091.66939472,
  11150190.415576423,
  11150499.535508212,
  11167139.795033026,
  11167146.733910596,
  11167181.207791744,
  11168387.958579028,
  11170831.25,
  11173257.8125,
  11173660.238696704,
  11207580.720249677,
  11208721.875,
  11217254.257382715,
  11217560.900176646,
  11219270.3125,
  11256710.771725755,
  11259235.662703464,
  11259769.066427844,
  11260432.158499919,
  11264139.104121426,
  11264947.645925201,
  11265092.216840269,
  11265641.60255622,
  11265812.5,
  11266009.312730523,
  11266062.158604996,
  11266139.0625,
  11266493.75,
  11266506.211218331,
  11266541.131823104,
  11266789.0625,
  11267146.533449108,
  11267248.4375,
  11267360.94005182,
  11267464.273117088,
  11267889.0625,
  11268016.517799657,
  11268050.0,
  11268100.470255062,
  11268243.516179841,
  11268710.301103912,
  11268909.814778535,
  11273924.728893256,
  11285058.507062903,
  11301622.869773796,
  11301678.125,
  11305879.008617057,
  11306688.891365174,
  11317179.6875,
  11317521.875,
  11325264.725981757,
  11327557.504718628,
  11335565.388846623,
  11335606.46004931,
  11343465.810245508,
  11343683.688706458,
  11358159.375,
  11361792.1875,
  11364869.84936274,
  11386646.995896783,
  11391195.3125,
  11613991.801342947,
  11644575.721773047,
  11825167.704729134,
  11826532.618146012,
  12359810.9375,
  12405788.561093709,
  12412255.992851308,
  12766435.9375,
  12775486.38418735,
  12776460.666151168,
  12777679.6875,
  12779045.3125,
  12780132.562330503,
  12781445.3125,
  12789219.27077604,
  12789290.3691226,
  12791514.493317692,
  12791669.926699363,
  12793248.4375,
  12793939.012875333,
  12797828.125,
  12800243.972104227,
  12800286.710933585,
  12804194.611044167,
  12833462.757222446,
  12847538.959554236,
  12852471.875,
  13076147.750639915,
  13137595.408848872,
  13164774.863675997,
  13174558.425561663,
  13203401.5625,
  13203407.67259257,
  13206251.808167202,
  13209997.076022692,
  13210596.875,
  14081595.109291937,
  14099023.430746213,
  14138059.809210049,
  14175562.858426778],
 [39.50316505699135,
  67.98118973901671,
  10.475974835742623,
  9.283645768105531,
  43.19733067911364,
  46.72084718618287,
  112.76887399498642,
  14.338767177944092,
  97.918332693213,
  10.433252092058375,
  11.154931133489224,
  5.578499140480954,
  171.0161099825164,
  8.818071195696492,
  12.843343636464379,
  16.882440591201906,
  11.351822441213697,
  19.4165965061271,
  72.28675643873743,
  14.504883258177252,
  5.317478152842318,
  15.125955356278041,
  21.915192867842556,
  28.43948359473771,
  7.052864853486413,
  52.67705388113174,
  110.18496295559889,
  12.321574972962075,
  54.28176176914436,
  33.524651031298205,
  13.29786520492926,
  7.122224623379808,
  25.200326373046824,
  100.9278902654627,
  53.47690291703544,
  45.68603105929498,
  28.34435113107704,
  63.55581010773996,
  22.06392500416911,
  59.25098737433268,
  137.2246871549226,
  69.04375388938725,
  36.47629596348966,
  37.628643560581224,
  55.25368581672657,
  9.464492020581726,
  15.189623174582835,
  12.07243313710575,
  112.73610992701094,
  22.65383545202508,
  78.30760580888762,
  74.52571466998572,
  109.36107687799321,
  49.25151524499393,
  76.76951799222392,
  97.38296466841138,
  184.67112861261216,
  21.820369773355335,
  14.469388563229382,
  56.27002224380116,
  40.07203909668597,
  5.370190083145916,
  8.032385437302798,
  42.95557843724215,
  46.03756987483408,
  18.239280315574085,
  23.517845529202095,
  11.660938696150717,
  141.82845458917473,
  30.12287385241698,
  100.12519723385515,
  10.916654295819153,
  20.525468473093618,
  18.118620291298082,
  45.249649132104395,
  81.52535769530303,
  164.840737797297,
  47.690690827365245,
  20.243315694287496,
  46.82265063992616,
  30.640038596597066,
  6.930264112494245,
  6.495586163170042,
  52.35241462809125,
  75.29059316127328,
  11.241744909791224,
  5.170583868778721,
  20.8779091691581,
  36.532671435122296,
  59.1264381076774,
  65.51480269991409,
  100.28263274581711,
  5.163271165077303,
  100.08339375069251,
  28.530200005352768,
  59.42226366510195,
  169.16689145982713,
  18.420521638537195,
  124.48521262540953,
  5.049886296354172,
  21.05510015536838,
  71.58755458122793,
  135.62392647896442,
  11.821480264401828,
  156.56124919286145,
  11.268197800474285,
  84.19985025652363,
  32.06203287820941,
  40.72358733787902,
  21.890162543640233,
  85.05769072728724,
  7.8986112063784795,
  43.61816050384583,
  21.148538191844224,
  120.1479857856996,
  12.701678672839142,
  18.30187995691229,
  9.836832189598512,
  21.969172175576638,
  87.53097201744308,
  109.78303601419351,
  30.754174396294367,
  16.827406308327546,
  33.69657480567423,
  46.84345867585989,
  54.49256816428344,
  21.13711115856213,
  63.914465018185616,
  5.098453072358372,
  54.14416381720697,
  54.155580063974526,
  43.50012159871536,
  38.40719992467651,
  59.46819460795218,
  36.12469615871113,
  64.77537698628942,
  16.540021706797532,
  118.93985599089858,
  119.37355715310139,
  12.631508817806465,
  61.52752982278932,
  11.172808747355564,
  167.78914028142808,
  8.318654898119922,
  12.082482934963709,
  14.320003277902503,
  13.189125169058485,
  64.20040837410554,
  60.743866726123784,
  148.0040344779177,
  63.506012846898585,
  33.037747486236604,
  83.61165944756007,
  66.68941854937256,
  5.507434673401909,
  12.241326015228964,
  92.07402479353142,
  19.232676850345797,
  183.04346327815165,
  69.11463909836347,
  18.712274366461852,
  77.40474405681056,
  5.834551274349883,
  63.39302136388497,
  59.17137936773298,
  71.56796665841725,
  12.280943448568983,
  45.72530776167984,
  46.31745114280217,
  34.630257348748444,
  29.896252614242968,
  57.603192790750796,
  60.88308765824151,
  14.543146064676076,
  7.477613189232802,
  31.339899846522307,
  10.501221232944795,
  45.89334633945878,
  12.380655616575535,
  32.80777126786844,
  39.94217339735157,
  58.657165281011466,
  56.15204557435317,
  22.09212438479612,
  62.22700927973831,
  38.225569100157756,
  73.45287949372364,
  15.55030579217439,
  8.550790715402066,
  62.156793227722886,
  48.694236095658574,
  8.519181382170569,
  127.49896016940642,
  74.35725229200243,
  10.239584531634653,
  79.31756973372542,
  27.187897341703405,
  39.676251016094014,
  25.270125138171736,
  6.580681302260346,
  21.438046325948587,
  68.52869301077271,
  56.32787557445083,
  6.699296911303664,
  84.81195249402639,
  6.156937611144566,
  79.34649857951442,
  7.328835452646465,
  7.558647047478901,
  9.84699546056425,
  122.46400322140755,
  90.37138020116767,
  123.23393382180436,
  21.5741035135771,
  31.973075449484124,
  211.05273636199007,
  52.3756682219781,
  54.37979379497091,
  181.80183785015723,
  15.793606803403478,
  12.695270854836156,
  27.480866527807308,
  21.074340742478125,
  54.16684174716094,
  88.261347702655,
  107.61382329947793,
  23.29596714529758,
  18.055367811545647,
  117.366133581356,
  86.57085999908834,
  58.5329977855276,
  118.53085534949352,
  114.28349954453392,
  16.922032320229476,
  17.798091778896588,
  20.458006591773838,
  9.976475147979615,
  18.769727652286633,
  99.9455425558037,
  43.89973990013543,
  49.80085948127204,
  5.74603120059032,
  38.130784164585314,
  105.11395066089355,
  11.231707206535946,
  9.005299157995132,
  13.92092632862266,
  14.575173908021776,
  167.46296337851243,
  65.60731307079084,
  52.59176147794144,
  16.77276589992591,
  16.41230238264284,
  32.77197772102608,
  39.08331626137966,
  22.80415572628247,
  141.91655220764213,
  91.92467701756958,
  121.07781002347892,
  57.666724434261084,
  31.38139508354997,
  13.836561093882093,
  27.8096839030851,
  10.374361731516549,
  43.2046311574842,
  82.70019922836426,
  6.729136485409227,
  21.365479383067,
  7.128689597436717,
  58.103539549638526,
  14.258546238283651,
  5.498563743302872,
  23.284961250720276,
  17.010685648832077,
  30.370300503328743,
  27.873587333053564,
  39.794891099579914,
  317.264460074431,
  42.869925292417975,
  20.12726928224934,
  31.087947989101828,
  92.85253580039058,
  17.796892604748905,
  51.86295698276669,
  47.60891241664507,
  20.039357856911572,
  64.65907003187658,
  17.418322391910742,
  8.974009552725779,
  41.25246632208623,
  19.9836427031011,
  87.60460094520714,
  205.70680637198802,
  8.575357914864341,
  91.70076528392454,
  39.94930390076782,
  38.724564142595966,
  6.211454123512151,
  34.50196643718585,
  66.19265898201763,
  43.368771562997715,
  7.367787491327854,
  70.64014478191929,
  16.748967683683865,
  21.080085331452803,
  37.02097284894502,
  93.29296954000951,
  56.57272601070615,
  63.11233906260891,
  26.7299243347594,
  21.62271658412902,
  20.60244840400898,
  12.206855918780509,
  158.21937990948868,
  43.81581669546195,
  5.9957804594264195,
  14.49804175605524,
  29.045476364368568,
  42.95685976086347,
  12.3794300179355,
  60.002422912892186,
  24.75588069153058,
  68.96586736966843,
  46.49526127280067,
  8.350319691507895,
  19.31364446519322,
  32.58273834981632,
  51.60347229964887,
  30.842982512175833,
  21.72821663263724,
  13.488260944579789,
  68.16291121039966,
  12.94152362513063,
  24.13371368218136,
  52.7024518326302,
  8.94399700941504,
  44.07779010452224,
  7.113821162289917,
  22.760919223852838,
  183.68952898995354,
  11.687854875288638,
  54.21955842375287,
  38.989546208267456,
  85.80772985010569,
  49.79673443711015,
  117.93790648148806,
  29.66464002724236,
  9.761263370426112,
  280.14832419444963,
  23.22041559779727,
  33.95434962496035,
  10.448093811092093,
  13.438438024661767,
  5.133425479576933,
  16.710918313127475,
  34.18430914018503,
  13.168468712223998,
  35.18378824343972,
  6.409052293536936,
  21.975517311805895,
  24.866679081739058,
  54.986971361983436,
  11.807327722191385,
  32.16828886480468,
  12.406246384690967,
  138.62882808050097,
  53.38215312488828,
  61.69675895424119,
  7.643477533354556,
  44.01396304080855,
  27.98612411810494,
  6.248243738001401,
  21.228974309955674,
  30.272422260672542,
  25.795128827537994,
  89.55463739983735,
  17.22852174774867,
  49.03696515142512,
  11.015995072244392,
  34.53788227779839,
  93.07063806943887,
  143.12699157188348,
  188.0788107248651,
  19.09171758762296,
  94.594987882581,
  217.69143758541898,
  48.021933137036754,
  42.433374335887734,
  5.427325995685693,
  69.10741624106771,
  58.44649630496571,
  77.27462336593153,
  79.92324475179814,
  6.115464905356458,
  15.541360360602457,
  119.83512328313321,
  245.2513389074015,
  54.95601040747751,
  39.79674531032294,
  36.33063957014039,
  19.13021494192061,
  17.880355775712836,
  7.269851314316384,
  36.67046067112894,
  196.83566563410648,
  11.144274027383975,
  121.45688832239892,
  115.09905926256451,
  42.95695223082459,
  17.006472493453117,
  55.235415287515856,
  36.350266621736125,
  23.092233267582078,
  75.94126158127837,
  79.54774983702579,
  67.52821268696961,
  107.95182017905461,
  26.269449616442465,
  86.69673658557966,
  8.750238064077863,
  53.143423608527215,
  18.534601510029795,
  29.47390563761767,
  86.54984289177627,
  15.83157506605092,
  5.292822371868798,
  64.86804594195122,
  45.62181094600685,
  26.175474528346665,
  35.136595532772105,
  105.84798670794297,
  42.31068596256322,
  12.752055955854289,
  35.01353804943281,
  29.912663654740633,
  148.96922531686255,
  21.397040671964866,
  5.754789353631851,
  27.898663322647938,
  6.960535101717085,
  17.758430239961733,
  208.73987030044356,
  121.08638988630685,
  66.87419870337315,
  28.533081436975788,
  81.88307268943217,
  10.437432965647504,
  31.655140345207386,
  21.555533522184323,
  26.339029500836943,
  29.830742328708673,
  40.917244378492335,
  51.328837124411564,
  100.72647055894237,
  249.71683219370942,
  45.29884271245434,
  85.72728235851778,
  14.296925893606659,
  33.436782034711484,
  89.22098673132109,
  93.00418892387427,
  106.04860353705973,
  24.35947589657486,
  240.5133756211648,
  55.29936135366497,
  12.924757640611173,
  5.558362076355664,
  32.09336785343787,
  72.61274505494005,
  39.933268901371434,
  10.665068051913394,
  37.44563436079499,
  13.833096750077926,
  24.076234157346146,
  51.36190317645799,
  52.06298262482526,
  34.962808104015224,
  6.126346234412228,
  10.56706923952123,
  69.42156034995841,
  5.8137186022315595,
  55.870555121813126,
  31.267887602902707,
  77.45468142604233,
  55.1674483705808,
  19.488894849414848,
  5.140235838273593,
  14.416684870112073,
  6.568376930671347,
  61.745195099021714,
  19.104353814473367,
  9.554015103434526,
  154.364746163956,
  180.50940237057165,
  43.88335866126284,
  67.15078573173517,
  56.702802705225466,
  6.171593579121499,
  133.64086025580562,
  26.370345058017897,
  15.881353712258855,
  44.291779230488956,
  23.59112212455028,
  185.12404191948178,
  11.924984305919258,
  50.81014347660295,
  29.608125076867246,
  8.704095647400791,
  63.85135895188379,
  10.267796839308943,
  239.83471601720842,
  10.160366551818761,
  13.522546916263728,
  117.65038163483892,
  15.69802687563399,
  54.29769705298167,
  23.505649161650165,
  82.71773438798344,
  112.15022448191117,
  11.246801638266403,
  38.897067836405846,
  28.105241657299036,
  77.08391804322174,
  78.32327295153378,
  140.5529760462183,
  17.784402725463,
  56.88055790759188,
  49.53040158816948,
  5.361616764259064,
  95.5865679863165,
  11.916664361816641,
  31.277784732601013,
  76.3292078122316,
  43.04791942118595,
  69.79503048522363,
  205.1453442959711,
  21.437658215864364,
  69.69174512797319,
  167.59064327992328,
  10.936663265682077,
  90.32961923125609,
  59.6360171846313,
  43.36172722626316,
  80.81973251501296,
  55.502445334401315,
  65.24663440316296,
  12.597069867649909,
  122.7851219197621,
  36.38648019411924,
  92.96538678721858,
  16.97568893827716,
  94.20382537397045,
  46.014881344188026,
  100.87989361553686,
  48.45975791076177,
  10.929903575409403,
  6.419188511802526,
  74.73375876292599,
  33.53142057881886,
  14.00738332081079,
  102.88951099268658,
  104.57337625198284,
  178.82023307639827,
  9.640875386187458,
  14.406743329393132,
  25.33909073688116,
  12.934952611126516,
  29.659763248731572,
  67.68060034522952,
  8.315265241217455,
  6.699890069788223,
  61.14502615351799,
  153.16604304872877,
  10.13174754877868,
  70.67247077433476,
  48.6727632232121,
  57.693098935846436,
  87.8705213561676,
  14.74996893220896,
  8.762262413819911,
  266.8341285953496,
  53.94273713016208,
  78.2918588799144,
  14.677048703049477,
  65.10229005857747,
  48.482791836267424,
  92.51818144712378,
  21.74348126099796,
  8.582633039813405,
  16.823354130900178,
  13.021753325625962,
  120.57941575875051,
  6.112385156712894,
  15.269893477015982,
  49.53842392320138,
  84.99749220461472,
  65.111582861139,
  21.19410190848156,
  83.18432435823185,
  65.96495571799105,
  5.746163229296661,
  36.53985710851037,
  15.597226634384171,
  30.121563948284003,
  83.80534800357908,
  83.05571949720654,
  105.40613805787011,
  17.63034772766711,
  19.503482933788554,
  172.20884847222757,
  172.76590816379965,
  62.66060381272337,
  5.745073480472719,
  46.26407272278802,
  184.27969168430545,
  106.1502881619002,
  77.12238463178318,
  23.461081744352935,
  48.0592746327572,
  8.928640150612145,
  122.39405626496612,
  21.854780752492147,
  12.50489211144153,
  15.511759110799773,
  87.41363178193195,
  9.751791601320102,
  51.72895020235779,
  57.09813386267075,
  5.269978964581719,
  61.99640944060233,
  81.24362552023702,
  38.160456006768655,
  14.801133306894108,
  34.71918848391615,
  61.126036860280095,
  60.02210479015393,
  74.7360472403996,
  20.08311870618091,
  15.102559950606805,
  36.2013445158081,
  228.8605933688175,
  7.999712780848196,
  51.6585053053881,
  85.09709289105838,
  106.94011337740989,
  111.70531486061282,
  46.83885217464058,
  312.4528765779274,
  8.574102087018632,
  297.0373019280596,
  23.23411036707691,
  35.65677939817195,
  14.951244695236474,
  6.110278859062512,
  84.16856601252927,
  8.822793188947056,
  96.46873951833324,
  6.1384060074468305,
  22.520152641449176,
  29.356014436532597,
  15.333481833557244,
  27.378530491661877,
  12.419924308306046,
  13.144309721098136,
  42.666640096519295,
  27.359965317166278,
  25.54930992101975,
  52.456059678950716,
  110.96406293338795,
  18.262249301253334,
  56.926206554031815,
  24.602055552335177,
  33.08469173311845,
  11.430152314651679,
  69.76586902768939,
  90.53347811252846,
  79.10642867431459,
  7.004805931076852,
  88.26067278520506,
  46.553460625703856,
  195.13321801985725,
  12.586149015353122,
  12.04541644021504,
  6.725315838947938,
  39.16008383245138,
  26.882421123351886,
  53.01498674088756,
  61.55531170779169,
  18.506147656527652,
  24.123752240435177,
  36.220924431059395,
  77.00875205535822,
  24.408833387887526,
  74.55698943557425,
  42.030491949670974,
  5.6906645331502475,
  16.88442088971928,
  17.72218341082624,
  149.495764297482,
  27.579697480683027,
  94.41426366208414,
  6.5224809174786085,
  22.789855823658115,
  59.75870854853966,
  16.021789474290898,
  62.13039931733603,
  32.372088360356216,
  93.0881293434344,
  93.50188585235915,
  31.13445806268396,
  7.292504726327672,
  73.75925724298725,
  14.433778254666748,
  5.83238393698658,
  13.82565030106506,
  80.39383503016572,
  6.72360460858169,
  60.33725145198902,
  11.8817623427884,
  36.811620496808466])

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

([8194426.5625,
  9007979.6875,
  9007980.788781097,
  9008015.567974329,
  9048940.625,
  9207851.5625,
  9222218.75,
  9224647.086230844,
  9225658.813289803,
  9300041.692603217,
  9321231.26734774,
  9381733.714107273,
  9439306.25,
  9752058.698369922,
  9752248.458814077,
  9949970.742946515,
  9951515.763801731,
  9952592.02227511,
  9954007.8125,
  9979901.443027167,
  9980293.0188882,
  9983306.297564307,
  9997508.916296797,
  10012672.370592983,
  10020371.607496593,
  10024610.753136452,
  10084867.1875,
  10155739.712103916,
  10160677.892555565,
  10160683.964612478,
  10163671.364678029,
  10172400.002939677,
  10176477.930598183,
  10185125.0,
  10186471.875,
  10189105.90446893,
  10189122.234375035,
  10190985.777817901,
  10192462.055379925,
  10208813.569724252,
  10211635.06068229,
  10211656.25,
  10212315.625,
  10212321.64251549,
  10213287.035756813,
  10213426.71581303,
  10213446.938607415,
  10217683.343096478,
  10218348.46623591,
  10218354.296484858,
  10222039.0625,
  10222967.181086302,
  10226560.9375,
  10226571.875,
  10226670.3125,
  10227224.577280194,
  10227726.510403039,
  10228258.188632103,
  10228268.620362328,
  10228625.0,
  10229437.487831319,
  10239182.052403504,
  10242891.673760792,
  10246754.6875,
  10246756.25,
  10247077.222481124,
  10248730.93154603,
  10255358.12880965,
  10256565.835860334,
  10257321.37043249,
  10329926.290634952,
  10331160.79814444,
  10331510.469889743,
  10335000.64140006,
  10336010.9375,
  10336084.375,
  10336988.970066717,
  10337535.9375,
  10337805.323144559,
  10337877.27304477,
  10346157.8125,
  10347713.281582655,
  10348215.678977918,
  10350388.883519586,
  10361455.961740501,
  10365614.13425318,
  10368029.910331147,
  10369880.18940451,
  10375499.863919131,
  10375662.830386726,
  10375785.9375,
  10379495.65879957,
  10381309.28501451,
  10383887.5,
  10383897.971979002,
  10389103.125,
  10396035.9375,
  10397480.374793146,
  10404267.99218849,
  10404319.239253357,
  10404496.32994227,
  10405163.724547233,
  10405798.4375,
  10406348.840215039,
  10406399.327740619,
  10406834.463060094,
  10407250.0,
  10407287.24484313,
  10408000.0,
  10409433.12570805,
  10409947.627913915,
  10410156.998947991,
  10410730.159676611,
  10410932.987844327,
  10413623.4375,
  10415888.317510428,
  10415892.856994668,
  10458247.417821053,
  10468743.296435447,
  10478046.554902913,
  10485881.80461816,
  10485990.625,
  10485994.170084119,
  10487242.61766642,
  10487470.904151736,
  10488639.0625,
  10490499.463557629,
  10493322.501128238,
  10502476.562339796,
  10519979.6875,
  10531096.71966659,
  10544617.1875,
  10567410.9375,
  10568839.0625,
  10573619.559459034,
  10573620.3125,
  10585144.468153704,
  10585417.1875,
  10591295.3125,
  10591563.28767183,
  10591740.625,
  10597931.27866446,
  10598174.753579237,
  10598704.747670537,
  10599034.910545133,
  10602705.906525807,
  10602914.066767057,
  10603428.125,
  10603832.8125,
  10603839.0625,
  10603950.0,
  10604618.37074226,
  10604981.25,
  10605607.8125,
  10605849.460008204,
  10606156.72160564,
  10607432.764858946,
  10608663.301518561,
  10609726.623548063,
  10609818.471214741,
  10610360.104812052,
  10610679.358630372,
  10610798.711701337,
  10610901.5625,
  10610974.520186258,
  10611125.0,
  10611155.318204867,
  10611221.875,
  10611276.5625,
  10612032.136345292,
  10612305.670112545,
  10613735.9375,
  10617029.429298695,
  10617707.902967792,
  10617709.013776207,
  10618025.859266253,
  10618030.808463143,
  10618389.540149285,
  10618733.433692362,
  10619320.3125,
  10619760.041487088,
  10620042.1875,
  10620482.8125,
  10620904.245523456,
  10620918.693427663,
  10621302.345684633,
  10621531.95838207,
  10621537.651053509,
  10621581.718308015,
  10621623.4375,
  10621746.875,
  10622053.540397298,
  10622203.00199093,
  10622490.625,
  10623146.221323734,
  10623150.676031128,
  10623317.306764124,
  10624888.492823606,
  10627462.901868073,
  10627750.572135545,
  10628219.80696334,
  10628414.0625,
  10629635.9375,
  10630165.294636628,
  10630193.75,
  10630830.3244337,
  10630904.6875,
  10631321.046748666,
  10631584.250007391,
  10631708.365214951,
  10631858.179277526,
  10632210.9375,
  10632470.710617732,
  10632704.195582036,
  10633071.875,
  10633225.662233654,
  10634573.4375,
  10646257.941417972,
  10647519.350374313,
  10649684.58193129,
  10650589.877306392,
  10650893.666895663,
  10650897.894471098,
  10651593.209459284,
  10652615.625,
  10654652.239153497,
  10655073.101043368,
  10655666.256528901,
  10656784.375,
  10661214.39199662,
  10663162.5,
  10663237.5,
  10663318.75,
  10664620.750345187,
  10664845.902922507,
  10665345.88729036,
  10665358.48814202,
  10665363.23798249,
  10665442.1875,
  10665827.522226505,
  10666101.678479921,
  10666102.162396852,
  10666323.4375,
  10666704.473170944,
  10666884.44053268,
  10666980.837394796,
  10666987.203726197,
  10667035.837948762,
  10667085.581458196,
  10667095.995219715,
  10667097.607153941,
  10667106.95328706,
  10667178.463881353,
  10667314.572177563,
  10667326.5625,
  10667641.809596898,
  10667670.3125,
  10667779.6875,
  10667800.0,
  10668045.147534115,
  10668198.4375,
  10668267.855773026,
  10668307.77531367,
  10668731.25827751,
  10668779.6875,
  10668841.683653843,
  10668934.200368706,
  10669082.761082401,
  10669284.402462967,
  10669330.688826755,
  10669388.870983042,
  10669628.955177063,
  10669840.287639402,
  10670051.508158518,
  10670923.4375,
  10670948.556673035,
  10671042.08900885,
  10671414.0625,
  10671583.594266461,
  10672114.976239199,
  10672806.25,
  10673075.0,
  10673099.02772723,
  10673718.75,
  10674424.723862963,
  10674846.215777814,
  10675068.75,
  10675311.464289881,
  10675768.29475451,
  10676151.398616692,
  10676241.521586902,
  10676295.3125,
  10676431.18793412,
  10677302.071385961,
  10677557.8125,
  10677801.5625,
  10677974.927239602,
  10677983.920668779,
  10677993.75,
  10678178.125,
  10678225.08399677,
  10678247.735952746,
  10678787.5,
  10678986.782303922,
  10679021.332575934,
  10679249.257874658,
  10679287.5,
  10679332.204993406,
  10679384.375,
  10679460.839701138,
  10680052.637895998,
  10684150.214178625,
  10684433.408242043,
  10684593.75,
  10685437.417640617,
  10687610.000355273,
  10689395.410292612,
  10693050.983483374,
  10695584.062223101,
  10698153.801184675,
  10698417.1875,
  10699159.149815282,
  10700074.411455559,
  10700093.75,
  10702228.243623562,
  10703377.601866776,
  10703675.930593522,
  10708579.382195435,
  10716253.125,
  10716897.029633995,
  10717238.620679403,
  10724123.4375,
  10724691.540826252,
  10725653.205679426,
  10726043.75,
  10726333.878310405,
  10726709.375,
  10727562.470892569,
  10728194.930847412,
  10729293.75,
  10729406.581059577,
  10729498.855642261,
  10729548.4375,
  10729829.6875,
  10729932.45185938,
  10729936.91494653,
  10729957.8125,
  10730022.343092557,
  10730109.706127003,
  10730210.594816752,
  10730469.954139352,
  10730619.829356158,
  10730629.63012114,
  10730650.136280507,
  10730949.616914203,
  10730994.304800713,
  10731255.691437561,
  10731546.875,
  10731749.483427815,
  10732935.495574178,
  10732960.720114041,
  10733059.605167136,
  10733125.395247223,
  10733255.906637331,
  10733288.18837584,
  10733372.61213823,
  10733376.352230968,
  10733454.661006652,
  10733465.734903507,
  10733587.961363181,
  10733687.472192083,
  10733815.473487189,
  10733826.035608636,
  10733967.959578259,
  10734204.53229983,
  10734271.890595283,
  10734535.740755893,
  10734557.010542015,
  10734563.989835024,
  10734698.35166462,
  10734950.0,
  10734983.761578597,
  10735127.649069538,
  10735565.459862305,
  10736184.375,
  10736194.876853177,
  10736328.382379571,
  10736607.981930826,
  10737511.672253339,
  10737551.5625,
  10737692.1875,
  10738206.25,
  10738556.606714938,
  10738656.509106094,
  10738981.398952805,
  10739422.992568307,
  10739457.8125,
  10739870.3125,
  10740257.652749237,
  10740350.76157476,
  10740422.662799142,
  10740749.041236116,
  10741681.78597816,
  10742285.9375,
  10742452.70571754,
  10744815.635573402,
  10744995.176785918,
  10746163.911209399,
  10748233.958437288,
  10748466.078470403,
  10748770.140936641,
  10749543.571246717,
  10749619.959137792,
  10749652.616260264,
  10749836.627291061,
  10750446.875,
  10750914.611101342,
  10751867.062037077,
  10752054.282003747,
  10752130.472185422,
  10752621.654243039,
  10753334.634885427,
  10754821.593598707,
  10756781.438443037,
  10757650.215437913,
  10757717.243570117,
  10757734.375,
  10761685.287727175,
  10761787.368747976,
  10762811.297617273,
  10763409.442316184,
  10784325.0,
  10785611.776125701,
  10790179.533992535,
  10800853.736558627,
  10801025.70014421,
  10801224.244481876,
  10801814.791647192,
  10804018.717560792,
  10807697.787808374,
  10807805.271445373,
  10810526.556379255,
  10821350.0,
  10823167.18693609,
  10830950.0,
  10831001.64710043,
  10831690.133419774,
  10833378.615434831,
  10849751.06873595,
  10849760.904473815,
  10869719.665163906,
  10883703.64894028,
  10917773.4375,
  10922868.096013507,
  10944681.426996509,
  10946166.54840526,
  10946276.5625,
  10946641.181216683,
  10946674.830690654,
  10949240.657753805,
  10949733.980566131,
  10950429.6875,
  10950707.265768236,
  10950790.623566462,
  10950809.375,
  10951138.409132151,
  10951328.635135477,
  10951362.5,
  10951438.888354348,
  10951744.967070037,
  10951809.05738127,
  10952185.9375,
  10952224.67525353,
  10952755.97506146,
  10952832.8125,
  10952913.567811498,
  10953097.931360053,
  10953344.109932166,
  10953513.671547595,
  10953585.9375,
  10953857.808962967,
  10953881.25,
  10953887.5,
  10954103.125,
  10954140.027678743,
  10954667.638055988,
  10955087.9499275,
  10955151.417637186,
  10955789.562155727,
  10956045.3125,
  10956146.116911612,
  10956183.064540097,
  10956238.241192807,
  10956475.552347831,
  10956545.3125,
  10956636.727818947,
  10959150.38932686,
  10959563.266369937,
  10960616.138274461,
  10961279.568197519,
  10961479.914005838,
  10961935.9375,
  10962015.312400375,
  10962158.352514409,
  10962309.831483819,
  10962422.80298365,
  10963235.9375,
  10963247.589060167,
  10963428.19313659,
  10964258.745098965,
  10965071.875,
  10970921.153333273,
  10970940.195213242,
  10980378.125,
  10987440.81026531,
  10987465.625,
  10998293.430567272,
  10998457.8125,
  11003134.375,
  11009473.512380058,
  11017965.625,
  11019489.040944865,
  11020406.245042948,
  11023040.420225363,
  11027090.56817726,
  11030069.338191263,
  11031638.820132403,
  11033573.4375,
  11033797.438723268,
  11034015.625,
  11036008.264716007,
  11036121.717488794,
  11047920.3125,
  11047936.234150114,
  11048089.304512007,
  11048518.544181656,
  11049419.806924181,
  11049706.25,
  11050223.4375,
  11050456.124238543,
  11050695.686735896,
  11050795.3125,
  11051185.579969808,
  11052187.5,
  11052459.375,
  11052498.4375,
  11052608.585341772,
  11052628.125,
  11052704.6875,
  11052807.794660637,
  11053052.010298248,
  11053057.8125,
  11053079.6875,
  11053420.3125,
  11053537.5,
  11053551.432155235,
  11053630.802711086,
  11053636.469264077,
  11053725.0,
  11053751.293748796,
  11053754.201753967,
  11054132.8125,
  11064935.9375,
  11066828.469645403,
  11069199.852247389,
  11070449.814942306,
  11070477.475705775,
  11072197.895028273,
  11073664.0625,
  11073696.928858122,
  11074027.797770614,
  11074724.74118652,
  11074772.514279012,
  11075249.253313374,
  11075279.075294733,
  11075294.666258208,
  11075624.921386193,
  11075840.625,
  11076190.74257594,
  11076355.490347927,
  11077287.09093583,
  11079744.392452262,
  11081800.0,
  11083407.741341138,
  11100244.93374215,
  11101520.3125,
  11102357.290362218,
  11107054.673681838,
  11107383.333312264,
  11108311.356817525,
  11109194.69750157,
  11111381.684249843,
  11118648.667844908,
  11120078.208567759,
  11120192.1875,
  11121236.025419911,
  11121564.746058064,
  11122366.65699521,
  11122967.1875,
  11123949.138291277,
  11126325.601857217,
  11126709.375,
  11127000.16201752,
  11127323.4375,
  11127749.04798044,
  11127944.628785584,
  11128120.3125,
  11128544.465129457,
  11128713.18315504,
  11130294.407833057,
  11131531.25,
  11131816.763923656,
  11132039.64816561,
  11134118.75,
  11134706.406815937,
  11135050.0,
  11140265.625,
  11147363.739480121,
  11147376.196770534,
  11148091.66939472,
  11150190.415576423,
  11150499.535508212,
  11167139.795033026,
  11167146.733910596,
  11167181.207791744,
  11168387.958579028,
  11170831.25,
  11173257.8125,
  11173660.238696704,
  11207580.720249677,
  11208721.875,
  11217254.257382715,
  11217560.900176646,
  11219270.3125,
  11256710.771725755,
  11259235.662703464,
  11259769.066427844,
  11260432.158499919,
  11264139.104121426,
  11264947.645925201,
  11265092.216840269,
  11265641.60255622,
  11265812.5,
  11266009.312730523,
  11266062.158604996,
  11266139.0625,
  11266493.75,
  11266506.211218331,
  11266541.131823104,
  11266789.0625,
  11267146.533449108,
  11267248.4375,
  11267360.94005182,
  11267464.273117088,
  11267889.0625,
  11268016.517799657,
  11268050.0,
  11268100.470255062,
  11268243.516179841,
  11268710.301103912,
  11268909.814778535,
  11273924.728893256,
  11285058.507062903,
  11301622.869773796,
  11301678.125,
  11305879.008617057,
  11306688.891365174,
  11317179.6875,
  11317521.875,
  11325264.725981757,
  11327557.504718628,
  11335565.388846623,
  11335606.46004931,
  11343465.810245508,
  11343683.688706458,
  11358159.375,
  11361792.1875,
  11364869.84936274,
  11386646.995896783,
  11391195.3125,
  11613991.801342947,
  11644575.721773047,
  11825167.704729134,
  11826532.618146012,
  12359810.9375,
  12405788.561093709,
  12412255.992851308,
  12766435.9375,
  12775486.38418735,
  12776460.666151168,
  12777679.6875,
  12779045.3125,
  12780132.562330503,
  12781445.3125,
  12789219.27077604,
  12789290.3691226,
  12791514.493317692,
  12791669.926699363,
  12793248.4375,
  12793939.012875333,
  12797828.125,
  12800243.972104227,
  12800286.710933585,
  12804194.611044167,
  12833462.757222446,
  12847538.959554236,
  12852471.875,
  13076147.750639915,
  13137595.408848872,
  13164774.863675997,
  13174558.425561663,
  13203401.5625,
  13203407.67259257,
  13206251.808167202,
  13209997.076022692,
  13210596.875,
  14081595.109291937,
  14099023.430746213,
  14138059.809210049,
  14175562.858426778],
 [39.50316505699135,
  67.98118973901671,
  10.475974835742623,
  9.283645768105531,
  43.19733067911364,
  46.72084718618287,
  112.76887399498642,
  14.338767177944092,
  97.918332693213,
  10.433252092058375,
  11.154931133489224,
  5.578499140480954,
  171.0161099825164,
  8.818071195696492,
  12.843343636464379,
  16.882440591201906,
  11.351822441213697,
  19.4165965061271,
  72.28675643873743,
  14.504883258177252,
  5.317478152842318,
  15.125955356278041,
  21.915192867842556,
  28.43948359473771,
  7.052864853486413,
  52.67705388113174,
  110.18496295559889,
  12.321574972962075,
  54.28176176914436,
  33.524651031298205,
  13.29786520492926,
  7.122224623379808,
  25.200326373046824,
  100.9278902654627,
  53.47690291703544,
  45.68603105929498,
  28.34435113107704,
  63.55581010773996,
  22.06392500416911,
  59.25098737433268,
  137.2246871549226,
  69.04375388938725,
  36.47629596348966,
  37.628643560581224,
  55.25368581672657,
  9.464492020581726,
  15.189623174582835,
  12.07243313710575,
  112.73610992701094,
  22.65383545202508,
  78.30760580888762,
  74.52571466998572,
  109.36107687799321,
  49.25151524499393,
  76.76951799222392,
  97.38296466841138,
  184.67112861261216,
  21.820369773355335,
  14.469388563229382,
  56.27002224380116,
  40.07203909668597,
  5.370190083145916,
  8.032385437302798,
  42.95557843724215,
  46.03756987483408,
  18.239280315574085,
  23.517845529202095,
  11.660938696150717,
  141.82845458917473,
  30.12287385241698,
  100.12519723385515,
  10.916654295819153,
  20.525468473093618,
  18.118620291298082,
  45.249649132104395,
  81.52535769530303,
  164.840737797297,
  47.690690827365245,
  20.243315694287496,
  46.82265063992616,
  30.640038596597066,
  6.930264112494245,
  6.495586163170042,
  52.35241462809125,
  75.29059316127328,
  11.241744909791224,
  5.170583868778721,
  20.8779091691581,
  36.532671435122296,
  59.1264381076774,
  65.51480269991409,
  100.28263274581711,
  5.163271165077303,
  100.08339375069251,
  28.530200005352768,
  59.42226366510195,
  169.16689145982713,
  18.420521638537195,
  124.48521262540953,
  5.049886296354172,
  21.05510015536838,
  71.58755458122793,
  135.62392647896442,
  11.821480264401828,
  156.56124919286145,
  11.268197800474285,
  84.19985025652363,
  32.06203287820941,
  40.72358733787902,
  21.890162543640233,
  85.05769072728724,
  7.8986112063784795,
  43.61816050384583,
  21.148538191844224,
  120.1479857856996,
  12.701678672839142,
  18.30187995691229,
  9.836832189598512,
  21.969172175576638,
  87.53097201744308,
  109.78303601419351,
  30.754174396294367,
  16.827406308327546,
  33.69657480567423,
  46.84345867585989,
  54.49256816428344,
  21.13711115856213,
  63.914465018185616,
  5.098453072358372,
  54.14416381720697,
  54.155580063974526,
  43.50012159871536,
  38.40719992467651,
  59.46819460795218,
  36.12469615871113,
  64.77537698628942,
  16.540021706797532,
  118.93985599089858,
  119.37355715310139,
  12.631508817806465,
  61.52752982278932,
  11.172808747355564,
  167.78914028142808,
  8.318654898119922,
  12.082482934963709,
  14.320003277902503,
  13.189125169058485,
  64.20040837410554,
  60.743866726123784,
  148.0040344779177,
  63.506012846898585,
  33.037747486236604,
  83.61165944756007,
  66.68941854937256,
  5.507434673401909,
  12.241326015228964,
  92.07402479353142,
  19.232676850345797,
  183.04346327815165,
  69.11463909836347,
  18.712274366461852,
  77.40474405681056,
  5.834551274349883,
  63.39302136388497,
  59.17137936773298,
  71.56796665841725,
  12.280943448568983,
  45.72530776167984,
  46.31745114280217,
  34.630257348748444,
  29.896252614242968,
  57.603192790750796,
  60.88308765824151,
  14.543146064676076,
  7.477613189232802,
  31.339899846522307,
  10.501221232944795,
  45.89334633945878,
  12.380655616575535,
  32.80777126786844,
  39.94217339735157,
  58.657165281011466,
  56.15204557435317,
  22.09212438479612,
  62.22700927973831,
  38.225569100157756,
  73.45287949372364,
  15.55030579217439,
  8.550790715402066,
  62.156793227722886,
  48.694236095658574,
  8.519181382170569,
  127.49896016940642,
  74.35725229200243,
  10.239584531634653,
  79.31756973372542,
  27.187897341703405,
  39.676251016094014,
  25.270125138171736,
  6.580681302260346,
  21.438046325948587,
  68.52869301077271,
  56.32787557445083,
  6.699296911303664,
  84.81195249402639,
  6.156937611144566,
  79.34649857951442,
  7.328835452646465,
  7.558647047478901,
  9.84699546056425,
  122.46400322140755,
  90.37138020116767,
  123.23393382180436,
  21.5741035135771,
  31.973075449484124,
  211.05273636199007,
  52.3756682219781,
  54.37979379497091,
  181.80183785015723,
  15.793606803403478,
  12.695270854836156,
  27.480866527807308,
  21.074340742478125,
  54.16684174716094,
  88.261347702655,
  107.61382329947793,
  23.29596714529758,
  18.055367811545647,
  117.366133581356,
  86.57085999908834,
  58.5329977855276,
  118.53085534949352,
  114.28349954453392,
  16.922032320229476,
  17.798091778896588,
  20.458006591773838,
  9.976475147979615,
  18.769727652286633,
  99.9455425558037,
  43.89973990013543,
  49.80085948127204,
  5.74603120059032,
  38.130784164585314,
  105.11395066089355,
  11.231707206535946,
  9.005299157995132,
  13.92092632862266,
  14.575173908021776,
  167.46296337851243,
  65.60731307079084,
  52.59176147794144,
  16.77276589992591,
  16.41230238264284,
  32.77197772102608,
  39.08331626137966,
  22.80415572628247,
  141.91655220764213,
  91.92467701756958,
  121.07781002347892,
  57.666724434261084,
  31.38139508354997,
  13.836561093882093,
  27.8096839030851,
  10.374361731516549,
  43.2046311574842,
  82.70019922836426,
  6.729136485409227,
  21.365479383067,
  7.128689597436717,
  58.103539549638526,
  14.258546238283651,
  5.498563743302872,
  23.284961250720276,
  17.010685648832077,
  30.370300503328743,
  27.873587333053564,
  39.794891099579914,
  317.264460074431,
  42.869925292417975,
  20.12726928224934,
  31.087947989101828,
  92.85253580039058,
  17.796892604748905,
  51.86295698276669,
  47.60891241664507,
  20.039357856911572,
  64.65907003187658,
  17.418322391910742,
  8.974009552725779,
  41.25246632208623,
  19.9836427031011,
  87.60460094520714,
  205.70680637198802,
  8.575357914864341,
  91.70076528392454,
  39.94930390076782,
  38.724564142595966,
  6.211454123512151,
  34.50196643718585,
  66.19265898201763,
  43.368771562997715,
  7.367787491327854,
  70.64014478191929,
  16.748967683683865,
  21.080085331452803,
  37.02097284894502,
  93.29296954000951,
  56.57272601070615,
  63.11233906260891,
  26.7299243347594,
  21.62271658412902,
  20.60244840400898,
  12.206855918780509,
  158.21937990948868,
  43.81581669546195,
  5.9957804594264195,
  14.49804175605524,
  29.045476364368568,
  42.95685976086347,
  12.3794300179355,
  60.002422912892186,
  24.75588069153058,
  68.96586736966843,
  46.49526127280067,
  8.350319691507895,
  19.31364446519322,
  32.58273834981632,
  51.60347229964887,
  30.842982512175833,
  21.72821663263724,
  13.488260944579789,
  68.16291121039966,
  12.94152362513063,
  24.13371368218136,
  52.7024518326302,
  8.94399700941504,
  44.07779010452224,
  7.113821162289917,
  22.760919223852838,
  183.68952898995354,
  11.687854875288638,
  54.21955842375287,
  38.989546208267456,
  85.80772985010569,
  49.79673443711015,
  117.93790648148806,
  29.66464002724236,
  9.761263370426112,
  280.14832419444963,
  23.22041559779727,
  33.95434962496035,
  10.448093811092093,
  13.438438024661767,
  5.133425479576933,
  16.710918313127475,
  34.18430914018503,
  13.168468712223998,
  35.18378824343972,
  6.409052293536936,
  21.975517311805895,
  24.866679081739058,
  54.986971361983436,
  11.807327722191385,
  32.16828886480468,
  12.406246384690967,
  138.62882808050097,
  53.38215312488828,
  61.69675895424119,
  7.643477533354556,
  44.01396304080855,
  27.98612411810494,
  6.248243738001401,
  21.228974309955674,
  30.272422260672542,
  25.795128827537994,
  89.55463739983735,
  17.22852174774867,
  49.03696515142512,
  11.015995072244392,
  34.53788227779839,
  93.07063806943887,
  143.12699157188348,
  188.0788107248651,
  19.09171758762296,
  94.594987882581,
  217.69143758541898,
  48.021933137036754,
  42.433374335887734,
  5.427325995685693,
  69.10741624106771,
  58.44649630496571,
  77.27462336593153,
  79.92324475179814,
  6.115464905356458,
  15.541360360602457,
  119.83512328313321,
  245.2513389074015,
  54.95601040747751,
  39.79674531032294,
  36.33063957014039,
  19.13021494192061,
  17.880355775712836,
  7.269851314316384,
  36.67046067112894,
  196.83566563410648,
  11.144274027383975,
  121.45688832239892,
  115.09905926256451,
  42.95695223082459,
  17.006472493453117,
  55.235415287515856,
  36.350266621736125,
  23.092233267582078,
  75.94126158127837,
  79.54774983702579,
  67.52821268696961,
  107.95182017905461,
  26.269449616442465,
  86.69673658557966,
  8.750238064077863,
  53.143423608527215,
  18.534601510029795,
  29.47390563761767,
  86.54984289177627,
  15.83157506605092,
  5.292822371868798,
  64.86804594195122,
  45.62181094600685,
  26.175474528346665,
  35.136595532772105,
  105.84798670794297,
  42.31068596256322,
  12.752055955854289,
  35.01353804943281,
  29.912663654740633,
  148.96922531686255,
  21.397040671964866,
  5.754789353631851,
  27.898663322647938,
  6.960535101717085,
  17.758430239961733,
  208.73987030044356,
  121.08638988630685,
  66.87419870337315,
  28.533081436975788,
  81.88307268943217,
  10.437432965647504,
  31.655140345207386,
  21.555533522184323,
  26.339029500836943,
  29.830742328708673,
  40.917244378492335,
  51.328837124411564,
  100.72647055894237,
  249.71683219370942,
  45.29884271245434,
  85.72728235851778,
  14.296925893606659,
  33.436782034711484,
  89.22098673132109,
  93.00418892387427,
  106.04860353705973,
  24.35947589657486,
  240.5133756211648,
  55.29936135366497,
  12.924757640611173,
  5.558362076355664,
  32.09336785343787,
  72.61274505494005,
  39.933268901371434,
  10.665068051913394,
  37.44563436079499,
  13.833096750077926,
  24.076234157346146,
  51.36190317645799,
  52.06298262482526,
  34.962808104015224,
  6.126346234412228,
  10.56706923952123,
  69.42156034995841,
  5.8137186022315595,
  55.870555121813126,
  31.267887602902707,
  77.45468142604233,
  55.1674483705808,
  19.488894849414848,
  5.140235838273593,
  14.416684870112073,
  6.568376930671347,
  61.745195099021714,
  19.104353814473367,
  9.554015103434526,
  154.364746163956,
  180.50940237057165,
  43.88335866126284,
  67.15078573173517,
  56.702802705225466,
  6.171593579121499,
  133.64086025580562,
  26.370345058017897,
  15.881353712258855,
  44.291779230488956,
  23.59112212455028,
  185.12404191948178,
  11.924984305919258,
  50.81014347660295,
  29.608125076867246,
  8.704095647400791,
  63.85135895188379,
  10.267796839308943,
  239.83471601720842,
  10.160366551818761,
  13.522546916263728,
  117.65038163483892,
  15.69802687563399,
  54.29769705298167,
  23.505649161650165,
  82.71773438798344,
  112.15022448191117,
  11.246801638266403,
  38.897067836405846,
  28.105241657299036,
  77.08391804322174,
  78.32327295153378,
  140.5529760462183,
  17.784402725463,
  56.88055790759188,
  49.53040158816948,
  5.361616764259064,
  95.5865679863165,
  11.916664361816641,
  31.277784732601013,
  76.3292078122316,
  43.04791942118595,
  69.79503048522363,
  205.1453442959711,
  21.437658215864364,
  69.69174512797319,
  167.59064327992328,
  10.936663265682077,
  90.32961923125609,
  59.6360171846313,
  43.36172722626316,
  80.81973251501296,
  55.502445334401315,
  65.24663440316296,
  12.597069867649909,
  122.7851219197621,
  36.38648019411924,
  92.96538678721858,
  16.97568893827716,
  94.20382537397045,
  46.014881344188026,
  100.87989361553686,
  48.45975791076177,
  10.929903575409403,
  6.419188511802526,
  74.73375876292599,
  33.53142057881886,
  14.00738332081079,
  102.88951099268658,
  104.57337625198284,
  178.82023307639827,
  9.640875386187458,
  14.406743329393132,
  25.33909073688116,
  12.934952611126516,
  29.659763248731572,
  67.68060034522952,
  8.315265241217455,
  6.699890069788223,
  61.14502615351799,
  153.16604304872877,
  10.13174754877868,
  70.67247077433476,
  48.6727632232121,
  57.693098935846436,
  87.8705213561676,
  14.74996893220896,
  8.762262413819911,
  266.8341285953496,
  53.94273713016208,
  78.2918588799144,
  14.677048703049477,
  65.10229005857747,
  48.482791836267424,
  92.51818144712378,
  21.74348126099796,
  8.582633039813405,
  16.823354130900178,
  13.021753325625962,
  120.57941575875051,
  6.112385156712894,
  15.269893477015982,
  49.53842392320138,
  84.99749220461472,
  65.111582861139,
  21.19410190848156,
  83.18432435823185,
  65.96495571799105,
  5.746163229296661,
  36.53985710851037,
  15.597226634384171,
  30.121563948284003,
  83.80534800357908,
  83.05571949720654,
  105.40613805787011,
  17.63034772766711,
  19.503482933788554,
  172.20884847222757,
  172.76590816379965,
  62.66060381272337,
  5.745073480472719,
  46.26407272278802,
  184.27969168430545,
  106.1502881619002,
  77.12238463178318,
  23.461081744352935,
  48.0592746327572,
  8.928640150612145,
  122.39405626496612,
  21.854780752492147,
  12.50489211144153,
  15.511759110799773,
  87.41363178193195,
  9.751791601320102,
  51.72895020235779,
  57.09813386267075,
  5.269978964581719,
  61.99640944060233,
  81.24362552023702,
  38.160456006768655,
  14.801133306894108,
  34.71918848391615,
  61.126036860280095,
  60.02210479015393,
  74.7360472403996,
  20.08311870618091,
  15.102559950606805,
  36.2013445158081,
  228.8605933688175,
  7.999712780848196,
  51.6585053053881,
  85.09709289105838,
  106.94011337740989,
  111.70531486061282,
  46.83885217464058,
  312.4528765779274,
  8.574102087018632,
  297.0373019280596,
  23.23411036707691,
  35.65677939817195,
  14.951244695236474,
  6.110278859062512,
  84.16856601252927,
  8.822793188947056,
  96.46873951833324,
  6.1384060074468305,
  22.520152641449176,
  29.356014436532597,
  15.333481833557244,
  27.378530491661877,
  12.419924308306046,
  13.144309721098136,
  42.666640096519295,
  27.359965317166278,
  25.54930992101975,
  52.456059678950716,
  110.96406293338795,
  18.262249301253334,
  56.926206554031815,
  24.602055552335177,
  33.08469173311845,
  11.430152314651679,
  69.76586902768939,
  90.53347811252846,
  79.10642867431459,
  7.004805931076852,
  88.26067278520506,
  46.553460625703856,
  195.13321801985725,
  12.586149015353122,
  12.04541644021504,
  6.725315838947938,
  39.16008383245138,
  26.882421123351886,
  53.01498674088756,
  61.55531170779169,
  18.506147656527652,
  24.123752240435177,
  36.220924431059395,
  77.00875205535822,
  24.408833387887526,
  74.55698943557425,
  42.030491949670974,
  5.6906645331502475,
  16.88442088971928,
  17.72218341082624,
  149.495764297482,
  27.579697480683027,
  94.41426366208414,
  6.5224809174786085,
  22.789855823658115,
  59.75870854853966,
  16.021789474290898,
  62.13039931733603,
  32.372088360356216,
  93.0881293434344,
  93.50188585235915,
  31.13445806268396,
  7.292504726327672,
  73.75925724298725,
  14.433778254666748,
  5.83238393698658,
  13.82565030106506,
  80.39383503016572,
  6.72360460858169,
  60.33725145198902,
  11.8817623427884,
  36.811620496808466])

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)