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

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"

#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 = 36529  #test discharge for which the notebook will definitely work
shot_no = 44524
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
    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))
    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 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=[]
        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))
            
            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') 
            i=i+1
    t3pa_cls.close
    return T_first, E_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)
    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.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.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)
    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].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

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'

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

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)

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

([5678503.125,
  5860962.5,
  5872554.6875,
  5884178.125,
  5888800.0,
  6060045.3125,
  6062376.5625,
  6062757.8125,
  6072537.5,
  6098928.125,
  6114381.25,
  6149687.5,
  6150795.3125,
  6164462.5,
  6183070.3125,
  6183078.125,
  6183126.5625,
  6183957.8125,
  6187428.125,
  6187665.625,
  6191129.6875,
  6191195.3125,
  6225615.625,
  6288882.8125,
  6289965.625,
  6308125.0,
  6314806.25,
  6319070.3125,
  6320145.3125,
  6322445.3125,
  6323484.375,
  6325187.5,
  6325307.8125,
  6326343.75,
  6327729.6875,
  6328107.8125,
  6328295.3125,
  6351643.75,
  6354384.375,
  6363343.75,
  6368209.375,
  6368776.5625,
  6379029.6875,
  6383468.75,
  6384890.625,
  6387459.375,
  6387607.8125,
  6388387.5,
  6392968.75,
  6401684.375,
  6403248.4375,
  6403882.8125,
  6404075.0,
  6404154.6875,
  6404998.4375,
  6406692.1875,
  6406868.75,
  6407050.0,
  6407400.0,
  6408132.8125,
  6408265.625,
  6408631.25,
  6408812.5,
  6413734.375,
  6414057.8125,
  6415409.375,
  6415915.625,
  6415917.1875,
  6418026.5625,
  6431112.5,
  6445409.375,
  6449481.25,
  6451356.25,
  6458643.75,
  6495289.0625,
  6497846.875,
  6498734.375,
  6499989.0625,
  6500110.9375,
  6501004.6875,
  6503451.5625,
  6503509.375,
  6503714.0625,
  6505385.9375,
  6507204.6875,
  6509682.8125,
  6513817.1875,
  6514254.6875,
  6514712.5,
  6518403.125,
  6518662.5,
  6526100.0,
  6528420.3125,
  6528634.375,
  6528989.0625,
  6529520.3125,
  6531948.4375,
  6536903.125,
  6536912.5,
  6539871.875,
  6541196.875,
  6544579.6875,
  6547056.25,
  6554398.4375,
  6560371.875,
  6560704.6875,
  6562389.0625,
  6562392.1875,
  6562395.3125,
  6564975.0,
  6566864.0625,
  6567648.4375,
  6569500.0,
  6572759.375,
  6575773.4375,
  6581865.625,
  6585104.6875,
  6585753.125,
  6587948.4375,
  6588715.625,
  6593353.125,
  6602042.1875,
  6604985.9375,
  6606676.5625,
  6609346.875,
  6609371.875,
  6610631.25,
  6611006.25,
  6611206.25,
  6611212.5,
  6611412.5,
  6611510.9375,
  6611664.0625,
  6611909.375,
  6612621.875,
  6612631.25,
  6613570.3125,
  6614082.8125,
  6614246.875,
  6614270.3125,
  6615162.5,
  6615198.4375,
  6615228.125,
  6615229.6875,
  6615293.75,
  6615837.5,
  6616484.375,
  6616518.75,
  6617845.3125,
  6618192.1875,
  6618314.0625,
  6618432.8125,
  6618553.125,
  6618834.375,
  6619000.0,
  6619382.8125,
  6619385.9375,
  6619445.3125,
  6619698.4375,
  6619957.8125,
  6619964.0625,
  6620306.25,
  6620470.3125,
  6620510.9375,
  6620579.6875,
  6620581.25,
  6620750.0,
  6621093.75,
  6621170.3125,
  6621343.75,
  6621367.1875,
  6621379.6875,
  6621400.0,
  6621406.25,
  6621492.1875,
  6621885.9375,
  6621889.0625,
  6622054.6875,
  6622654.6875,
  6622775.0,
  6623178.125,
  6623179.6875,
  6623465.625,
  6623529.6875,
  6623551.5625,
  6623648.4375,
  6623704.6875,
  6624196.875,
  6624471.875,
  6624584.375,
  6624640.625,
  6624750.0,
  6624975.0,
  6625257.8125,
  6625257.8125,
  6625315.625,
  6625353.125,
  6625510.9375,
  6625515.625,
  6625965.625,
  6626201.5625,
  6626478.125,
  6626506.25,
  6626592.1875,
  6626923.4375,
  6627325.0,
  6627342.1875,
  6627926.5625,
  6628662.5,
  6628750.0,
  6629064.0625,
  6629412.5,
  6629768.75,
  6630779.6875,
  6631170.3125,
  6631545.3125,
  6634039.0625,
  6634050.0,
  6634210.9375,
  6635254.6875,
  6635893.75,
  6636296.875,
  6636431.25,
  6637026.5625,
  6637601.5625,
  6638598.4375,
  6639932.8125,
  6646700.0,
  6646770.3125,
  6654910.9375,
  6655240.625,
  6655378.125,
  6656523.4375,
  6658023.4375,
  6662957.8125,
  6663925.0,
  6664485.9375,
  6668675.0,
  6670096.875,
  6672360.9375,
  6674323.4375,
  6675815.625,
  6676489.0625,
  6678357.8125,
  6678648.4375,
  6680739.0625,
  6683292.1875,
  6683803.125,
  6687245.3125,
  6687429.6875,
  6688012.5,
  6690857.8125,
  6697565.625,
  6698504.6875,
  6703725.0,
  6705039.0625,
  6710067.1875,
  6710565.625,
  6721034.375,
  6722598.4375,
  6723471.875,
  6723831.25,
  6725176.5625,
  6727451.5625,
  6727943.75,
  6728465.625,
  6730565.625,
  6731046.875,
  6733153.125,
  6733410.9375,
  6733471.875,
  6733742.1875,
  6734029.6875,
  6734193.75,
  6735429.6875,
  6735775.0,
  6735975.0,
  6736114.0625,
  6737390.625,
  6737484.375,
  6737603.125,
  6738626.5625,
  6739487.5,
  6739673.4375,
  6740064.0625,
  6740576.5625,
  6741485.9375,
  6741565.625,
  6742414.0625,
  6742414.0625,
  6742623.4375,
  6742801.5625,
  6743048.4375,
  6743175.0,
  6744157.8125,
  6744782.8125,
  6745332.8125,
  6745620.3125,
  6746023.4375,
  6746356.25,
  6747031.25,
  6747898.4375,
  6748695.3125,
  6749314.0625,
  6749348.4375,
  6749381.25,
  6750228.125,
  6752625.0,
  6753278.125,
  6754379.6875,
  6757500.0,
  6758064.0625,
  6758256.25,
  6760420.3125,
  6760815.625,
  6760851.5625,
  6761587.5,
  6763392.1875,
  6763410.9375,
  6764268.75,
  6765159.375,
  6766542.1875,
  6767042.1875,
  6769521.875,
  6777809.375,
  6779679.6875,
  6780801.5625,
  6780850.0,
  6781104.6875,
  6782050.0,
  6782209.375,
  6782543.75,
  6782792.1875,
  6783254.6875,
  6783643.75,
  6784407.8125,
  6786764.0625,
  6786784.375,
  6786982.8125,
  6787457.8125,
  6787948.4375,
  6788417.1875,
  6788537.5,
  6788556.25,
  6788739.0625,
  6788881.25,
  6788934.375,
  6789032.8125,
  6789140.625,
  6789268.75,
  6789973.4375,
  6790004.6875,
  6790459.375,
  6790929.6875,
  6791135.9375,
  6791154.6875,
  6791239.0625,
  6792198.4375,
  6792254.6875,
  6792759.375,
  6792898.4375,
  6793490.625,
  6793732.8125,
  6793926.5625,
  6794034.375,
  6794114.0625,
  6794187.5,
  6794189.0625,
  6794207.8125,
  6794559.375,
  6794575.0,
  6794739.0625,
  6794751.5625,
  6795009.375,
  6795168.75,
  6795190.625,
  6795210.9375,
  6795228.125,
  6795243.75,
  6795290.625,
  6795678.125,
  6795764.0625,
  6796340.625,
  6796482.8125,
  6796726.5625,
  6797367.1875,
  6797665.625,
  6797901.5625,
  6798106.25,
  6798120.3125,
  6798428.125,
  6798550.0,
  6798664.0625,
  6798951.5625,
  6799101.5625,
  6799198.4375,
  6799215.625,
  6799923.4375,
  6800028.125,
  6800048.4375,
  6800150.0,
  6800357.8125,
  6800451.5625,
  6800528.125,
  6800873.4375,
  6800907.8125,
  6801228.125,
  6801301.5625,
  6801434.375,
  6801617.1875,
  6801932.8125,
  6802395.3125,
  6802421.875,
  6802481.25,
  6802651.5625,
  6802707.8125,
  6802732.8125,
  6802787.5,
  6802820.3125,
  6802964.0625,
  6803150.0,
  6803214.0625,
  6803259.375,
  6803351.5625,
  6803457.8125,
  6803504.6875,
  6803598.4375,
  6803653.125,
  6803689.0625,
  6803703.125,
  6803775.0,
  6803798.4375,
  6804164.0625,
  6804168.75,
  6804193.75,
  6804204.6875,
  6804207.8125,
  6804210.9375,
  6804214.0625,
  6804220.3125,
  6804248.4375,
  6804256.25,
  6804259.375,
  6804512.5,
  6804514.0625,
  6804543.75,
  6804635.9375,
  6804653.125,
  6804793.75,
  6804853.125,
  6804917.1875,
  6804945.3125,
  6805145.3125,
  6805218.75,
  6805304.6875,
  6805348.4375,
  6805385.9375,
  6805387.5,
  6805501.5625,
  6805539.0625,
  6805578.125,
  6805607.8125,
  6805720.3125,
  6805750.0,
  6805773.4375,
  6805773.4375,
  6805901.5625,
  6805937.5,
  6806021.875,
  6806051.5625,
  6806189.0625,
  6806209.375,
  6806298.4375,
  6806320.3125,
  6806321.875,
  6806618.75,
  6807000.0,
  6807101.5625,
  6807142.1875,
  6807285.9375,
  6807290.625,
  6807475.0,
  6807482.8125,
  6807551.5625,
  6807798.4375,
  6807928.125,
  6807979.6875,
  6808868.75,
  6809035.9375,
  6809103.125,
  6809131.25,
  6809871.875,
  6810701.5625,
  6810985.9375,
  6812331.25,
  6814312.5,
  6841392.1875,
  6842440.625,
  6842443.75,
  6843315.625,
  6843464.0625,
  6843492.1875,
  6846635.9375,
  6847753.125,
  6847776.5625,
  6876957.8125,
  6878031.25,
  6878120.3125,
  6879685.9375,
  6881803.125,
  6882735.9375,
  6888564.0625,
  6888709.375,
  6891264.0625,
  6892821.875,
  6893875.0,
  6895110.9375,
  6898823.4375,
  6898892.1875,
  6900846.875,
  6902585.9375,
  6903382.8125,
  6904778.125,
  6905640.625,
  6906429.6875,
  6907373.4375,
  6907781.25,
  6907792.1875,
  6907843.75,
  6908192.1875,
  6908265.625,
  6910601.5625,
  6912381.25,
  6912834.375,
  6914184.375,
  6914898.4375,
  6915556.25,
  6916575.0,
  6916762.5,
  6916864.0625,
  6917006.25,
  6917273.4375,
  6917315.625,
  6917476.5625,
  6918867.1875,
  6920043.75,
  6920695.3125,
  6920700.0,
  6921462.5,
  6923875.0,
  6925793.75,
  6925960.9375,
  6926371.875,
  6926898.4375,
  6929328.125,
  6931801.5625,
  6932778.125,
  6935623.4375,
  6936379.6875,
  6952595.3125,
  6954873.4375,
  6958770.3125,
  6960568.75,
  6960575.0,
  6961759.375,
  6963076.5625,
  6963746.875,
  6965035.9375,
  6968595.3125,
  6968618.75,
  6970173.4375,
  6970673.4375,
  6973121.875,
  6975589.0625,
  6975589.0625,
  6979523.4375,
  6986015.625,
  6990473.4375,
  6992821.875,
  6992875.0,
  6995557.8125,
  7000578.125,
  7003475.0,
  7009964.0625,
  7010067.1875,
  7012059.375,
  7013076.5625,
  7017481.25,
  7019320.3125,
  7019371.875,
  7021131.25,
  7023870.3125,
  7025142.1875,
  7025893.75,
  7029600.0,
  7029676.5625,
  7030925.0,
  7036026.5625,
  7040059.375,
  7045996.875,
  7050882.8125,
  7051054.6875,
  7059025.0,
  7060748.4375,
  7065354.6875,
  7065585.9375,
  7065696.875,
  7066231.25,
  7066287.5,
  7066312.5,
  7066521.875,
  7066557.8125,
  7066704.6875,
  7067154.6875,
  7068790.625,
  7069056.25,
  7069173.4375,
  7069310.9375,
  7069462.5,
  7069679.6875,
  7069912.5,
  7070106.25,
  7070953.125,
  7071235.9375,
  7071382.8125,
  7071450.0,
  7072067.1875,
  7072995.3125,
  7073070.3125,
  7075135.9375,
  7075754.6875,
  7077654.6875,
  7080306.25,
  7083682.8125,
  7087618.75,
  7089595.3125,
  7089709.375,
  7089978.125,
  7091476.5625,
  7091921.875,
  7092782.8125,
  7093464.0625,
  7093693.75,
  7097728.125,
  7097762.5,
  7099778.125,
  7100476.5625,
  7101195.3125,
  7101676.5625,
  7110992.1875,
  7114667.1875,
  7114687.5,
  7135595.3125,
  7135773.4375,
  7136393.75,
  7136567.1875,
  7137306.25,
  7138226.5625,
  7139157.8125,
  7140085.9375,
  7140645.3125,
  7140681.25,
  7141696.875,
  7143968.75,
  7144054.6875,
  7144100.0,
  7144278.125,
  7144456.25,
  7144517.1875,
  7144782.8125,
  7145368.75,
  7146046.875,
  7146190.625,
  7146235.9375,
  7146409.375,
  7147031.25,
  7147806.25,
  7148392.1875,
  7153271.875,
  7153467.1875,
  7154043.75,
  7154815.625,
  7154976.5625,
  7155234.375,
  7155742.1875,
  7155750.0,
  7155759.375,
  7155807.8125,
  7157684.375,
  7160137.5,
  7163331.25,
  7166090.625,
  7170282.8125,
  7178064.0625,
  7181048.4375,
  7181742.1875,
  7181750.0,
  7185629.6875,
  7188131.25,
  7190898.4375,
  7192676.5625,
  7194932.8125,
  7199528.125,
  7200376.5625,
  7202873.4375,
  7206435.9375,
  7233745.3125,
  7238456.25,
  7240337.5,
  7245131.25,
  7245793.75,
  7246970.3125,
  7267000.0,
  7268742.1875,
  7270745.3125,
  7272240.625,
  7274007.8125,
  7274582.8125,
  7277534.375,
  7278976.5625,
  7282473.4375,
  7291060.9375,
  7298778.125,
  7351706.25,
  7360393.75,
  7363103.125,
  7363831.25,
  7367239.0625,
  7374210.9375,
  7379806.25,
  7381143.75,
  7386112.5,
  7390381.25,
  7394418.75,
  7395090.625,
  7395090.625,
  7406987.5,
  7411218.75,
  7420918.75,
  7421657.8125,
  7431370.3125,
  7431389.0625,
  7438256.25,
  7446681.25,
  7449803.125,
  7451342.1875,
  7465565.625,
  7502345.3125,
  7506750.0,
  7509904.6875,
  7517221.875,
  7525262.5,
  7538918.75,
  7548454.6875,
  7548785.9375,
  7550428.125,
  7550450.0,
  7551910.9375,
  7564014.0625,
  7577807.8125,
  7601307.8125,
  7604479.6875,
  7606565.625,
  7608357.8125,
  7609064.0625,
  7629371.875,
  7640804.6875,
  7644239.0625,
  7645853.125,
  7659410.9375,
  7663459.375,
  7667162.5,
  7685168.75,
  7686807.8125,
  7688215.625,
  7688368.75,
  7690368.75,
  7788470.3125,
  7791253.125,
  7796962.5,
  7800895.3125,
  7817806.25,
  7866092.1875,
  7881675.0,
  7901190.625,
  7902153.125,
  7903989.0625,
  7904062.5,
  7905389.0625,
  7907606.25,
  7909485.9375,
  7912476.5625,
  7912518.75,
  7915131.25,
  7949098.4375,
  7967892.1875,
  7967943.75,
  7969754.6875,
  7973856.25,
  7977109.375,
  7979151.5625,
  7984501.5625,
  8003348.4375,
  8033926.5625,
  8041746.875,
  8042765.625,
  8218093.75,
  8222915.625,
  8225373.4375,
  8232354.6875,
  8254564.0625,
  8257664.0625,
  8264907.8125,
  8267310.9375,
  8268610.9375,
  8268626.5625,
  8270878.125,
  8270910.9375,
  8273314.0625,
  8274415.625,
  8277448.4375,
  8277653.125,
  8278329.6875,
  8374012.5,
  8402148.4375,
  8406387.5,
  8409414.0625,
  8409714.0625,
  8465943.75,
  8577912.5,
  8577976.5625,
  8582915.625,
  8590328.125,
  8710014.0625,
  8711340.625,
  8711365.625,
  8715759.375,
  8806160.9375,
  8810582.8125,
  8858293.75,
  8950881.25,
  8975806.25,
  8981576.5625,
  8994664.0625,
  9012298.4375,
  9018550.0,
  9018968.75,
  9025235.9375,
  9057146.875,
  9057575.0,
  9167290.625,
  9167465.625,
  9167543.75,
  9167726.5625,
  9280321.875,
  9284970.3125,
  9307468.75,
  9319876.5625,
  9323296.875,
  9332540.625,
  9358239.0625,
  9368064.0625,
  9369412.5,
  9388545.3125,
  9390225.0,
  9391110.9375,
  9393790.625,
  9394431.25,
  9546664.0625,
  9549668.75,
  9550220.3125,
  9551631.25,
  9551732.8125,
  9552384.375,
  9552435.9375,
  9554171.875,
  9556809.375,
  9604979.6875,
  9704865.625,
  9750707.8125,
  9763701.5625,
  9842953.125,
  9858634.375,
  9883668.75,
  9892403.125,
  9898618.75,
  9905582.8125,
  9931401.5625,
  9935009.375,
  9971907.8125,
  10094093.75,
  10104489.0625,
  10105764.0625,
  10119478.125,
  10141587.5,
  10142315.625,
  10148984.375,
  10260142.1875,
  10262071.875,
  10976793.75,
  10985671.875,
  11081867.1875,
  11085589.0625,
  11088032.8125,
  11089739.0625,
  11090273.4375,
  11090878.125,
  11090971.875,
  11091028.125,
  11091190.625,
  11091320.3125,
  11091643.75,
  11091657.8125,
  11091701.5625,
  11093714.0625,
  11094746.875,
  11094756.25,
  11094843.75,
  11094887.5,
  11095435.9375,
  11095676.5625,
  11095984.375,
  11097200.0,
  11097270.3125,
  11098148.4375,
  11098903.125,
  11099368.75,
  11099498.4375,
  11099810.9375,
  11099860.9375,
  11099898.4375,
  11099960.9375,
  11099976.5625,
  11100250.0,
  11100485.9375,
  11100831.25,
  11101481.25,
  11101570.3125,
  11101578.125,
  11102368.75,
  11102415.625,
  11102873.4375,
  11104481.25,
  11106231.25,
  11108228.125,
  11108443.75,
  11109021.875,
  11110164.0625,
  11111443.75,
  11111971.875,
  11112696.875,
  11114500.0,
  11115729.6875,
  11116270.3125,
  11116929.6875,
  11117157.8125,
  11117815.625,
  11117992.1875,
  11120812.5,
  11120823.4375,
  11122151.5625,
  11126929.6875,
  11127943.75,
  11143981.25,
  11145668.75,
  11147581.25],
 [23.09646681816278,
  76.41414913220474,
  25.075078733203807,
  17.898604521814896,
  78.76138453820957,
  147.54782317156085,
  43.84313436386716,
  7.833122647288747,
  60.29336420415457,
  5.44352088317474,
  6.935434419250196,
  54.49208325453253,
  83.42235272965655,
  9.755770621741968,
  13.464362350621332,
  37.383263790445454,
  9.751771452360021,
  7.773770450052657,
  46.7573095722054,
  10.430067941002834,
  39.55250744823101,
  17.758857026542117,
  40.68133884254476,
  33.74940800963735,
  30.580128199506913,
  16.964085514229957,
  15.004966896612386,
  245.18515574827353,
  5.890444415388376,
  19.136731114907494,
  27.68792197258647,
  5.596824457932807,
  39.000931588622926,
  74.6455818728687,
  6.765714189037943,
  47.51222106614824,
  12.06725468525347,
  16.677661528556477,
  16.16999994114641,
  8.062978836147936,
  46.162819240992185,
  156.74331799869452,
  37.71091157723907,
  56.48503749092269,
  10.869049450635831,
  24.166134979983333,
  52.34594265703679,
  11.252806776483052,
  39.815794913069055,
  41.64119877348424,
  35.1740284401459,
  27.870468537171902,
  7.306225789525635,
  25.447705441255916,
  37.29166752073966,
  34.10814381034519,
  31.114001235457323,
  12.941470497674922,
  77.27955928476669,
  108.98125412119299,
  5.3398275472173165,
  18.527871245900037,
  38.805576006593455,
  79.99212316262845,
  35.516098161421354,
  85.50813448089428,
  84.08170406323603,
  24.55550101636636,
  22.58398745702308,
  24.10580114876118,
  8.048760514120213,
  59.26493170503068,
  39.50766036043862,
  45.987494532412576,
  55.63748291742418,
  72.91849444647515,
  95.68263776290124,
  86.77436429018091,
  24.685054403488337,
  94.61352608197888,
  53.57752154467659,
  83.65786433103156,
  81.46964795179817,
  11.8095624608373,
  54.5482839923181,
  52.745290908921895,
  5.52930637892141,
  14.481070457927395,
  8.65333440859475,
  130.18674058854342,
  107.31760502008677,
  55.66153402761941,
  21.4402597057734,
  37.45138916764656,
  11.450338501938521,
  5.600274801692124,
  45.03289155919002,
  13.019672836360414,
  14.508268209836615,
  70.09213361995104,
  7.699263183584653,
  88.89142974947379,
  34.300862620141665,
  38.14559947008174,
  27.92221445976094,
  29.48743251472165,
  88.71025095302657,
  19.077015125877832,
  10.310533446307305,
  12.778054730716493,
  23.816304618104432,
  25.20174802558589,
  105.84118914771959,
  11.18597591530499,
  6.360335623868345,
  13.680083523795059,
  26.010808352740497,
  21.435728426955674,
  20.24348663800454,
  7.287026655368721,
  18.97657204196281,
  25.090317154221403,
  122.50238002421648,
  28.209934201111654,
  13.250750964895458,
  158.27526495782521,
  106.32741086267562,
  53.27582112658458,
  18.821001625418006,
  21.80317910608514,
  13.709757200410335,
  34.35407720956256,
  63.463368243833386,
  60.746622542136336,
  20.58591512859365,
  49.542736088842595,
  8.227872078351501,
  88.09822791070869,
  9.881491038395996,
  53.702460440828695,
  5.3447145148315665,
  106.1402034063197,
  94.95346518361814,
  54.03440814621889,
  83.66897597288896,
  96.81380570974585,
  6.673972891229778,
  41.65819208481658,
  14.549226478082966,
  11.548226351795945,
  47.0255136529825,
  16.328130398289087,
  23.87731238925581,
  115.73847989490514,
  90.61983541701595,
  40.20101668185603,
  117.93539228475228,
  55.79221911734366,
  61.181829751242844,
  31.960025575937426,
  7.3531650387301255,
  34.3891722865601,
  147.2563781752445,
  51.429841456300935,
  16.63375120613278,
  6.0477098742011215,
  80.78014513560368,
  99.4661951827968,
  25.38416878533699,
  58.47022260621292,
  330.2268603532193,
  49.67101210244098,
  81.81598277577862,
  22.249064939666027,
  26.51450836481815,
  176.52379229725503,
  80.31796248175466,
  37.34719957276074,
  14.979334571526724,
  100.93020368272109,
  13.163255659860852,
  16.21049502265484,
  12.722743943066144,
  71.87459248461144,
  10.581285617606289,
  14.807699984013801,
  8.13267345698736,
  131.14444740991922,
  20.491821853036573,
  18.225103235688298,
  22.372282947273654,
  7.468647539550206,
  8.473255416139036,
  7.830244444720974,
  58.25810435812028,
  71.42919666632234,
  8.02490941236178,
  53.131349293729045,
  158.97531573670602,
  109.23348209066081,
  12.68398497328271,
  11.855992948508753,
  17.23002872297242,
  35.56269833299826,
  27.91433198748642,
  66.14264690405385,
  76.36977542702137,
  16.744557449002546,
  75.00794447623903,
  7.471494869889003,
  7.516320133033916,
  7.641250526344109,
  41.01519066018683,
  28.59302217934698,
  27.135047561329532,
  15.490584425037806,
  38.41889344663898,
  10.474118457382808,
  63.1097333017173,
  13.771638454431148,
  196.35440403837347,
  45.33554632990689,
  43.170288930191234,
  30.50494422418303,
  65.11519956445282,
  89.87990784332499,
  68.59200654282664,
  69.17145479141683,
  30.60723710460382,
  19.216705459006064,
  55.33718636806174,
  21.904487595471856,
  85.29922293685348,
  11.306723591833897,
  80.83224547548106,
  27.27600346659839,
  43.283596900698974,
  19.52004102043005,
  66.86170386999859,
  18.55791504706697,
  44.71218593961408,
  76.15520357708716,
  180.40320086682667,
  57.061742459971526,
  27.6705878581027,
  68.33743480511626,
  18.620400431268088,
  24.530800057901743,
  45.41470309582157,
  58.47919510094424,
  25.125724069295355,
  60.330054943416016,
  24.995603050199712,
  17.815204255937715,
  27.004487503945562,
  72.27093023267976,
  27.9191763044152,
  49.5886643752733,
  90.17286374720061,
  91.15992024378313,
  103.96258110876433,
  5.282378293816456,
  43.65068488987147,
  80.06988595427131,
  12.08540678604747,
  63.60527539957998,
  13.278586783970654,
  14.36800305439691,
  18.69044990012204,
  9.038911259198743,
  5.043591136343069,
  60.89134810469868,
  126.12918335024405,
  47.29182710158029,
  11.256055453485402,
  10.08873642581787,
  17.177918251110125,
  11.551871204318578,
  79.49388357412488,
  46.83315449200037,
  91.81541689983051,
  57.182059266192965,
  38.18307005211219,
  8.573750701242474,
  68.7897489082834,
  49.55515559462687,
  54.6617441482015,
  13.684773019959055,
  8.597130172921673,
  56.43386582842407,
  121.76610873983014,
  93.3820218085413,
  114.77288350816902,
  33.57555467030754,
  68.70577305279268,
  126.18799902077353,
  11.806736893432724,
  67.5717274962382,
  31.09733312297488,
  24.9993178571767,
  50.91404483677396,
  74.99383888751169,
  17.164033408826647,
  53.625767602512845,
  5.077356712305322,
  17.048460788630074,
  24.490395523338314,
  32.51869845102124,
  7.40936319877423,
  14.008185744776586,
  58.83469407136514,
  15.96515050247702,
  140.4294180233351,
  15.00029204586558,
  58.959776751222755,
  29.00786598454628,
  22.385333494226465,
  68.08628641581663,
  88.3513608283919,
  57.52898471488527,
  9.170096934300823,
  59.202925898044924,
  10.383975078825948,
  15.908029817261188,
  68.44877639548497,
  13.964203147775955,
  38.12228485606426,
  14.657511901910114,
  55.24969493735518,
  28.418871179080792,
  13.634985284050115,
  20.593114874826995,
  7.268482365145665,
  13.79962092497581,
  8.765552540738394,
  7.052060997673832,
  83.56842636929363,
  44.714391618498134,
  38.14198033587367,
  14.593740943570765,
  14.15180293486147,
  126.14168179685748,
  47.96128619383243,
  26.788844991694212,
  34.126671958694054,
  138.51132064101247,
  33.89650422085187,
  139.78208201277732,
  17.92959290678197,
  14.058029814699449,
  11.700523940850971,
  63.92208444707154,
  68.67228810727593,
  67.17039380611489,
  82.80194814505042,
  65.55138104825598,
  159.2600351333081,
  15.186617889429188,
  19.564137007425693,
  16.56356673971569,
  7.4023145860666855,
  66.03785765602143,
  6.717644962731825,
  56.43136262214742,
  176.94409330818468,
  59.981433873743,
  89.21498829623785,
  21.48360170085487,
  6.803592364634948,
  15.452444273169021,
  67.33265802241313,
  75.33599995109903,
  29.369952890623452,
  53.35715545937019,
  33.1669456365453,
  88.0331667248147,
  82.28566960822343,
  72.11796308355532,
  22.64400155187669,
  19.100835933485165,
  16.23567369086813,
  11.608350448985593,
  7.232920380501984,
  48.77432665652317,
  18.870389674506736,
  26.612544181833155,
  76.42600867198205,
  14.759024019388667,
  31.133951694881247,
  23.202970867937008,
  8.68879198888946,
  112.10675524404846,
  5.249294819257113,
  73.16996781306247,
  10.088316178286297,
  12.626010192571085,
  71.06954982983893,
  96.4113934677713,
  25.716591665227995,
  74.70676345939859,
  15.061637140305567,
  31.349563703591844,
  14.765575628876116,
  41.37783741683598,
  6.607019620649585,
  20.016470084221307,
  21.39283167792491,
  54.415168874197455,
  89.48726452015974,
  8.974007576842856,
  10.37863948876639,
  10.010938030597671,
  160.44211687263984,
  5.425689625397781,
  5.988969731469766,
  14.868070754044982,
  16.703555716568903,
  5.032782901148494,
  25.747973838410505,
  12.732413670650539,
  175.34566492353284,
  23.424524974862887,
  32.0514130049495,
  91.65272460705333,
  39.36689932641691,
  14.38267800892593,
  26.6598596835178,
  33.63610288230363,
  98.20603715132502,
  65.26401332348384,
  313.33811760509525,
  61.039257054260936,
  89.77645657215288,
  19.995126261665384,
  38.79491955904897,
  26.78565655374219,
  8.227283651834314,
  15.077980884520434,
  50.93695406713049,
  56.412173704381,
  5.379946272545311,
  42.65603231855647,
  9.03805733257106,
  16.795191646795338,
  25.13183788510992,
  181.1554443409089,
  56.27861656223156,
  25.60811747633166,
  5.848399199904395,
  11.124725613811325,
  23.88846446546136,
  29.68447464326836,
  16.284857881132073,
  53.90884739056179,
  5.12946536475925,
  9.087724984729412,
  31.370536497249784,
  67.21793660523171,
  11.087284574123233,
  69.57536470147195,
  10.71050659714387,
  15.026406033500642,
  114.17159311802381,
  29.037351510597386,
  19.881878828036456,
  78.2986779092786,
  24.661532508463115,
  156.115212201316,
  60.461577464050656,
  46.51681115564065,
  64.58622582448633,
  9.989703034361797,
  53.15762941177178,
  85.27307199795871,
  176.305296323644,
  31.837548026754792,
  22.34985351742204,
  43.37568616079928,
  6.460963445932313,
  31.50834246245548,
  20.58552746852196,
  44.94705353616304,
  37.17461838773451,
  5.099194077332632,
  32.90993039607598,
  104.9916000738912,
  90.29404954860809,
  21.284921249870045,
  63.362470157772734,
  61.21107144910174,
  39.99791073761996,
  29.648512739081998,
  115.00463624065695,
  36.22922519753022,
  42.537829192352824,
  64.45626154113465,
  54.51225385741079,
  80.82654142283819,
  56.97474121039772,
  6.2683575364450475,
  6.3022993995514724,
  34.43252960671323,
  17.084033084814408,
  11.510704211004917,
  13.846648250539005,
  25.24238478363363,
  51.31691896485763,
  10.627951732475035,
  51.23120075689179,
  13.937854078563294,
  20.77498925745532,
  90.58011909951925,
  86.28657754369108,
  28.166031570200573,
  43.34354334706103,
  63.755564071280325,
  8.440272586614352,
  47.423728401763235,
  15.60628441250115,
  6.462352017126078,
  46.34345583843195,
  42.644513503265884,
  67.05489588593181,
  57.700090619293356,
  50.23355206617134,
  29.93019455834159,
  42.323728366451064,
  324.6762747630948,
  6.62776076943185,
  90.55236982052737,
  5.192636666306374,
  74.55620479146661,
  29.349448639369495,
  26.334210243547137,
  90.41760447507804,
  14.649971147491229,
  6.477199593589102,
  59.8287600317738,
  47.55215898458397,
  208.2547132477722,
  11.25065545428929,
  7.618536229633273,
  29.21729821751388,
  5.396782012222748,
  12.447092500610562,
  24.526174011643924,
  102.91069207739035,
  76.19809814420171,
  67.7368598307686,
  66.39835136692548,
  28.66129661975702,
  15.87962279485578,
  26.23462793271085,
  40.88067496199566,
  9.299442518713496,
  15.304711307160304,
  15.15215809710493,
  160.08728455182893,
  32.315523665824294,
  14.641212781246002,
  22.689494937996766,
  51.39737824257823,
  24.507053897445232,
  6.615755445004151,
  41.53936272775266,
  32.64512891996418,
  6.499348054041406,
  39.956518034349926,
  55.220384188716416,
  7.855200405854162,
  26.23293667230579,
  94.02840278444549,
  137.5732317944795,
  11.74106718209256,
  53.216812116956646,
  9.121092394572802,
  26.978770443785013,
  23.976403344690638,
  131.13900161001277,
  5.512709807354648,
  83.50587695776316,
  120.75754111154025,
  21.432627649464347,
  15.938218510478407,
  5.50792898869451,
  8.349560441435239,
  46.04330878224397,
  42.5215676405118,
  204.2566448446522,
  15.540972179187387,
  21.829801157422736,
  75.73574552508313,
  30.100486366045565,
  18.725531276530628,
  63.65354681442656,
  179.41647082102958,
  71.29779825656281,
  18.026845898146256,
  50.580100325133216,
  25.374202899652573,
  117.99260613097047,
  29.569000140727198,
  83.50072556539898,
  20.17129171681023,
  5.152129416571529,
  12.418697704452997,
  12.918247060528927,
  117.13813439275386,
  18.818056661350806,
  44.77227416705334,
  10.863722794633293,
  17.257761057982442,
  20.913078397867118,
  46.23452095016029,
  14.500311047434908,
  51.84250759072827,
  26.786096373728334,
  8.854665134358058,
  75.59907570488885,
  5.7794721015084685,
  16.311281090981943,
  44.1448549007064,
  11.90691444290629,
  66.32186209786246,
  8.88679758115977,
  5.3769170970019236,
  15.825348266967726,
  248.7647104338708,
  95.88248605413318,
  42.2519594346983,
  59.32768348766648,
  130.97140295118984,
  58.570285304433995,
  22.916988222888275,
  69.45381813591031,
  39.744556718555536,
  23.630047931521005,
  16.20378871730779,
  49.85046416041436,
  72.04361428474628,
  100.64556626023996,
  104.25108810520656,
  7.702191769034201,
  54.38835104227716,
  14.39549290924738,
  39.19447563840014,
  5.249218109261389,
  8.54203346361912,
  75.08793963255036,
  6.603855332714815,
  48.617068927368884,
  46.24467225265641,
  7.23011467546573,
  5.292418075307008,
  30.24723549342636,
  20.507620622784927,
  45.39737759814344,
  15.11899561294872,
  17.971645091223593,
  9.73349556867445,
  14.198727066882569,
  5.80524972819613,
  35.03792022298571,
  148.5710923952991,
  37.08055407342612,
  5.8737625917498555,
  29.335296710939083,
  45.148205970740364,
  48.20780377183974,
  13.897105916045245,
  105.22865824518311,
  62.41919324789416,
  71.19244883367816,
  18.461070889653755,
  47.67862493246133,
  20.242123879078814,
  10.880988443715635,
  31.141345074159133,
  5.637329293642736,
  23.265373497688167,
  9.93942213205608,
  70.24366398277022,
  70.64072638104335,
  11.627387288026377,
  29.40552434894623,
  11.79661355602266,
  25.17408269921721,
  8.028697841877118,
  95.86410563706423,
  54.15152188844783,
  7.6522072900293345,
  89.27910220284586,
  25.878933171395296,
  113.09481307963549,
  70.26634073319221,
  106.3503912943476,
  17.972447492576638,
  5.385817299348203,
  13.07458698314503,
  14.009872863302448,
  9.615292714087882,
  10.646363787476574,
  17.54145887959419,
  47.361354483070414,
  50.97301100206579,
  34.55933999788362,
  7.997463485011918,
  37.48030616404397,
  19.552018062377808,
  63.41710539045099,
  80.13770821137499,
  9.539568409049592,
  20.72629194435229,
  106.13315063164029,
  8.774496474086412,
  83.06297317549364,
  10.997196061981022,
  9.877489272470132,
  31.168933551461365,
  29.943901691072227,
  12.856009314396326,
  5.38110965223008,
  40.57188921262078,
  27.064820768877112,
  163.84810559956765,
  48.919726586645695,
  6.1713507826946055,
  16.65902539830648,
  148.78754953733747,
  12.420142059606889,
  16.92961780800558,
  38.43958774119711,
  81.94086042832038,
  19.83292741911609,
  53.37252801495728,
  48.740876131409756,
  75.47926639866053,
  139.69361338466956,
  17.578225571721887,
  11.340479901121961,
  6.194325530092807,
  21.090540718314756,
  46.88065943311076,
  58.030598924109796,
  185.99291227923194,
  6.604815087221636,
  11.34559824691836,
  10.77523594242167,
  134.63759007171868,
  22.303078670593184,
  11.462877499945172,
  111.20880094189025,
  49.33848601370531,
  107.42201470004261,
  37.451382911049514,
  29.716720293362197,
  89.94955215673986,
  6.072013284829785,
  28.963532023674382,
  126.91197967481754,
  60.25867905104224,
  20.691830293424896,
  78.29741749805024,
  54.37888426891168,
  16.937010273912406,
  7.352485570199575,
  97.59141527917566,
  54.79435777258159,
  108.27722485514072,
  92.26586728032673,
  11.579327749792661,
  12.431137668497986,
  13.577408393609055,
  40.9028747991367,
  41.16533159228178,
  112.71383667687874,
  11.31823441804757,
  72.9319093451652,
  16.638607330865103,
  10.440241915268214,
  155.98499119284438,
  26.8434550999985,
  32.54190314571381,
  10.519969282958632,
  118.51697667439582,
  29.154088198990237,
  40.72790912104638,
  53.04413414909574,
  60.907404841461506,
  146.38469734345168,
  8.642480698497081,
  71.43106764356533,
  28.402211411925478,
  22.065729394112545,
  68.10017307000676,
  77.00978165031549,
  5.309758889323243,
  27.99778135881106,
  18.661261356024767,
  54.69931075962305,
  33.81723448596374,
  82.41399349664184,
  35.47818686275496,
  19.08572738311426,
  32.41250777936446,
  34.35461860032439,
  24.415892744744156,
  5.13062813839644,
  13.612852055873017,
  52.009854991578194,
  5.580271875296687,
  26.442735264009773,
  5.746522042028261,
  20.328272984819076,
  21.965844840627856,
  11.929193471490276,
  20.79543070707279,
  22.89528470495999,
  12.183623574882127,
  24.3092891425887,
  74.3320937831429,
  28.127506546643282,
  19.483343445062122,
  182.37556629039452,
  32.25122617212112,
  15.600633867317907,
  57.24508359765869,
  8.048476474559743,
  20.242070402730736,
  56.63173754070698,
  5.036582668575584,
  5.152819356371734,
  18.94931258753779,
  10.112816446706383,
  65.31845825774111,
  8.885312023415915,
  8.299693222890182,
  30.474495787567868,
  5.634890417683859,
  15.578056641766537,
  26.38535609145219,
  17.812532591067523,
  15.266133509541646,
  6.314812793814502,
  20.4580238254229,
  33.10733152014791,
  34.98019841220422,
  6.223366604155347,
  26.163352467987494,
  18.177070756327915,
  5.047114330222604,
  17.10749424453423,
  13.416464862116142,
  15.544373022705157,
  63.20163433849206,
  16.11500586051435,
  38.872233052806855,
  7.277746090901515,
  8.757433423977686,
  40.600452373231555,
  7.0815913205529295,
  5.925992373342192,
  24.62537968757753,
  57.854075101505906,
  20.175740873767914,
  28.87070208056962,
  64.47387747461397,
  88.24292921784107,
  14.097615983920917,
  33.13578221402966,
  20.563481770435782,
  27.18087927802395,
  5.606927873887151,
  16.412155775802216,
  6.451790461334643,
  21.935676144929356,
  12.14082664236699,
  22.36388594216339,
  8.340142717409163,
  6.678983011082577,
  117.98343041161252,
  6.693610119808251,
  8.712969274987193,
  41.39149943080047,
  5.562127579715948,
  56.637391453645485,
  88.54937524550088,
  14.837974107430604,
  15.53116409073412,
  16.928024137461136,
  27.77850471320413,
  19.149004144307735,
  22.151958147466146,
  152.45963455967888,
  14.129868939481765,
  7.965612117307012,
  47.785447491745174,
  42.94388536101974,
  5.861437667978293,
  86.69090644946715,
  46.266655787238484,
  14.680345605269244,
  50.56865604680539,
  66.77499383058839,
  26.720263285608375,
  8.972851083138643,
  38.63317139744809,
  21.170853121561734,
  18.170342667268365,
  16.85524444517276,
  76.31174879923557,
  13.18518358807878,
  6.3255234517095325,
  28.336549432247992,
  39.86597397288527,
  19.315255818557585,
  14.917933649184565,
  14.661277350816956,
  5.198887787396429,
  12.353963742261262,
  52.24058896466551,
  13.706679082816027,
  17.293808976705286,
  58.1097379060424,
  12.107441245541743,
  10.174894450405562,
  41.79909314140049,
  33.18711371552374,
  6.532810932434207,
  5.783015978335344,
  49.658516952763996,
  5.1686063468196695,
  71.031184499903,
  18.244887258021098,
  70.92163024034326,
  34.84218192673624,
  120.18083984850456,
  15.594978685148265,
  9.485102212893365,
  15.134936863091614,
  10.859802037020275,
  8.687640683123133,
  84.05623191754816,
  26.995388524964685,
  23.41861402244887,
  55.981966176233556,
  27.94251376559162,
  16.093502188542093,
  17.552646316599713,
  19.87104944265511,
  50.110761261535,
  99.36883421197973,
  6.525270257340417,
  6.794619366257612,
  16.753376032615233,
  16.055384113524283,
  56.36045723468583,
  63.768470484820625,
  62.894571415145165,
  16.63593416169257,
  26.267111242820263,
  14.683615690575087,
  53.32980836243981,
  15.647127897861955,
  11.060981384999367,
  6.573082482109846])

Energy spectrum

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

([5678503.125,
  5860962.5,
  5872554.6875,
  5884178.125,
  5888800.0,
  6060045.3125,
  6062376.5625,
  6062757.8125,
  6072537.5,
  6098928.125,
  6114381.25,
  6149687.5,
  6150795.3125,
  6164462.5,
  6183070.3125,
  6183078.125,
  6183126.5625,
  6183957.8125,
  6187428.125,
  6187665.625,
  6191129.6875,
  6191195.3125,
  6225615.625,
  6288882.8125,
  6289965.625,
  6308125.0,
  6314806.25,
  6319070.3125,
  6320145.3125,
  6322445.3125,
  6323484.375,
  6325187.5,
  6325307.8125,
  6326343.75,
  6327729.6875,
  6328107.8125,
  6328295.3125,
  6351643.75,
  6354384.375,
  6363343.75,
  6368209.375,
  6368776.5625,
  6379029.6875,
  6383468.75,
  6384890.625,
  6387459.375,
  6387607.8125,
  6388387.5,
  6392968.75,
  6401684.375,
  6403248.4375,
  6403882.8125,
  6404075.0,
  6404154.6875,
  6404998.4375,
  6406692.1875,
  6406868.75,
  6407050.0,
  6407400.0,
  6408132.8125,
  6408265.625,
  6408631.25,
  6408812.5,
  6413734.375,
  6414057.8125,
  6415409.375,
  6415915.625,
  6415917.1875,
  6418026.5625,
  6431112.5,
  6445409.375,
  6449481.25,
  6451356.25,
  6458643.75,
  6495289.0625,
  6497846.875,
  6498734.375,
  6499989.0625,
  6500110.9375,
  6501004.6875,
  6503451.5625,
  6503509.375,
  6503714.0625,
  6505385.9375,
  6507204.6875,
  6509682.8125,
  6513817.1875,
  6514254.6875,
  6514712.5,
  6518403.125,
  6518662.5,
  6526100.0,
  6528420.3125,
  6528634.375,
  6528989.0625,
  6529520.3125,
  6531948.4375,
  6536903.125,
  6536912.5,
  6539871.875,
  6541196.875,
  6544579.6875,
  6547056.25,
  6554398.4375,
  6560371.875,
  6560704.6875,
  6562389.0625,
  6562392.1875,
  6562395.3125,
  6564975.0,
  6566864.0625,
  6567648.4375,
  6569500.0,
  6572759.375,
  6575773.4375,
  6581865.625,
  6585104.6875,
  6585753.125,
  6587948.4375,
  6588715.625,
  6593353.125,
  6602042.1875,
  6604985.9375,
  6606676.5625,
  6609346.875,
  6609371.875,
  6610631.25,
  6611006.25,
  6611206.25,
  6611212.5,
  6611412.5,
  6611510.9375,
  6611664.0625,
  6611909.375,
  6612621.875,
  6612631.25,
  6613570.3125,
  6614082.8125,
  6614246.875,
  6614270.3125,
  6615162.5,
  6615198.4375,
  6615228.125,
  6615229.6875,
  6615293.75,
  6615837.5,
  6616484.375,
  6616518.75,
  6617845.3125,
  6618192.1875,
  6618314.0625,
  6618432.8125,
  6618553.125,
  6618834.375,
  6619000.0,
  6619382.8125,
  6619385.9375,
  6619445.3125,
  6619698.4375,
  6619957.8125,
  6619964.0625,
  6620306.25,
  6620470.3125,
  6620510.9375,
  6620579.6875,
  6620581.25,
  6620750.0,
  6621093.75,
  6621170.3125,
  6621343.75,
  6621367.1875,
  6621379.6875,
  6621400.0,
  6621406.25,
  6621492.1875,
  6621885.9375,
  6621889.0625,
  6622054.6875,
  6622654.6875,
  6622775.0,
  6623178.125,
  6623179.6875,
  6623465.625,
  6623529.6875,
  6623551.5625,
  6623648.4375,
  6623704.6875,
  6624196.875,
  6624471.875,
  6624584.375,
  6624640.625,
  6624750.0,
  6624975.0,
  6625257.8125,
  6625257.8125,
  6625315.625,
  6625353.125,
  6625510.9375,
  6625515.625,
  6625965.625,
  6626201.5625,
  6626478.125,
  6626506.25,
  6626592.1875,
  6626923.4375,
  6627325.0,
  6627342.1875,
  6627926.5625,
  6628662.5,
  6628750.0,
  6629064.0625,
  6629412.5,
  6629768.75,
  6630779.6875,
  6631170.3125,
  6631545.3125,
  6634039.0625,
  6634050.0,
  6634210.9375,
  6635254.6875,
  6635893.75,
  6636296.875,
  6636431.25,
  6637026.5625,
  6637601.5625,
  6638598.4375,
  6639932.8125,
  6646700.0,
  6646770.3125,
  6654910.9375,
  6655240.625,
  6655378.125,
  6656523.4375,
  6658023.4375,
  6662957.8125,
  6663925.0,
  6664485.9375,
  6668675.0,
  6670096.875,
  6672360.9375,
  6674323.4375,
  6675815.625,
  6676489.0625,
  6678357.8125,
  6678648.4375,
  6680739.0625,
  6683292.1875,
  6683803.125,
  6687245.3125,
  6687429.6875,
  6688012.5,
  6690857.8125,
  6697565.625,
  6698504.6875,
  6703725.0,
  6705039.0625,
  6710067.1875,
  6710565.625,
  6721034.375,
  6722598.4375,
  6723471.875,
  6723831.25,
  6725176.5625,
  6727451.5625,
  6727943.75,
  6728465.625,
  6730565.625,
  6731046.875,
  6733153.125,
  6733410.9375,
  6733471.875,
  6733742.1875,
  6734029.6875,
  6734193.75,
  6735429.6875,
  6735775.0,
  6735975.0,
  6736114.0625,
  6737390.625,
  6737484.375,
  6737603.125,
  6738626.5625,
  6739487.5,
  6739673.4375,
  6740064.0625,
  6740576.5625,
  6741485.9375,
  6741565.625,
  6742414.0625,
  6742414.0625,
  6742623.4375,
  6742801.5625,
  6743048.4375,
  6743175.0,
  6744157.8125,
  6744782.8125,
  6745332.8125,
  6745620.3125,
  6746023.4375,
  6746356.25,
  6747031.25,
  6747898.4375,
  6748695.3125,
  6749314.0625,
  6749348.4375,
  6749381.25,
  6750228.125,
  6752625.0,
  6753278.125,
  6754379.6875,
  6757500.0,
  6758064.0625,
  6758256.25,
  6760420.3125,
  6760815.625,
  6760851.5625,
  6761587.5,
  6763392.1875,
  6763410.9375,
  6764268.75,
  6765159.375,
  6766542.1875,
  6767042.1875,
  6769521.875,
  6777809.375,
  6779679.6875,
  6780801.5625,
  6780850.0,
  6781104.6875,
  6782050.0,
  6782209.375,
  6782543.75,
  6782792.1875,
  6783254.6875,
  6783643.75,
  6784407.8125,
  6786764.0625,
  6786784.375,
  6786982.8125,
  6787457.8125,
  6787948.4375,
  6788417.1875,
  6788537.5,
  6788556.25,
  6788739.0625,
  6788881.25,
  6788934.375,
  6789032.8125,
  6789140.625,
  6789268.75,
  6789973.4375,
  6790004.6875,
  6790459.375,
  6790929.6875,
  6791135.9375,
  6791154.6875,
  6791239.0625,
  6792198.4375,
  6792254.6875,
  6792759.375,
  6792898.4375,
  6793490.625,
  6793732.8125,
  6793926.5625,
  6794034.375,
  6794114.0625,
  6794187.5,
  6794189.0625,
  6794207.8125,
  6794559.375,
  6794575.0,
  6794739.0625,
  6794751.5625,
  6795009.375,
  6795168.75,
  6795190.625,
  6795210.9375,
  6795228.125,
  6795243.75,
  6795290.625,
  6795678.125,
  6795764.0625,
  6796340.625,
  6796482.8125,
  6796726.5625,
  6797367.1875,
  6797665.625,
  6797901.5625,
  6798106.25,
  6798120.3125,
  6798428.125,
  6798550.0,
  6798664.0625,
  6798951.5625,
  6799101.5625,
  6799198.4375,
  6799215.625,
  6799923.4375,
  6800028.125,
  6800048.4375,
  6800150.0,
  6800357.8125,
  6800451.5625,
  6800528.125,
  6800873.4375,
  6800907.8125,
  6801228.125,
  6801301.5625,
  6801434.375,
  6801617.1875,
  6801932.8125,
  6802395.3125,
  6802421.875,
  6802481.25,
  6802651.5625,
  6802707.8125,
  6802732.8125,
  6802787.5,
  6802820.3125,
  6802964.0625,
  6803150.0,
  6803214.0625,
  6803259.375,
  6803351.5625,
  6803457.8125,
  6803504.6875,
  6803598.4375,
  6803653.125,
  6803689.0625,
  6803703.125,
  6803775.0,
  6803798.4375,
  6804164.0625,
  6804168.75,
  6804193.75,
  6804204.6875,
  6804207.8125,
  6804210.9375,
  6804214.0625,
  6804220.3125,
  6804248.4375,
  6804256.25,
  6804259.375,
  6804512.5,
  6804514.0625,
  6804543.75,
  6804635.9375,
  6804653.125,
  6804793.75,
  6804853.125,
  6804917.1875,
  6804945.3125,
  6805145.3125,
  6805218.75,
  6805304.6875,
  6805348.4375,
  6805385.9375,
  6805387.5,
  6805501.5625,
  6805539.0625,
  6805578.125,
  6805607.8125,
  6805720.3125,
  6805750.0,
  6805773.4375,
  6805773.4375,
  6805901.5625,
  6805937.5,
  6806021.875,
  6806051.5625,
  6806189.0625,
  6806209.375,
  6806298.4375,
  6806320.3125,
  6806321.875,
  6806618.75,
  6807000.0,
  6807101.5625,
  6807142.1875,
  6807285.9375,
  6807290.625,
  6807475.0,
  6807482.8125,
  6807551.5625,
  6807798.4375,
  6807928.125,
  6807979.6875,
  6808868.75,
  6809035.9375,
  6809103.125,
  6809131.25,
  6809871.875,
  6810701.5625,
  6810985.9375,
  6812331.25,
  6814312.5,
  6841392.1875,
  6842440.625,
  6842443.75,
  6843315.625,
  6843464.0625,
  6843492.1875,
  6846635.9375,
  6847753.125,
  6847776.5625,
  6876957.8125,
  6878031.25,
  6878120.3125,
  6879685.9375,
  6881803.125,
  6882735.9375,
  6888564.0625,
  6888709.375,
  6891264.0625,
  6892821.875,
  6893875.0,
  6895110.9375,
  6898823.4375,
  6898892.1875,
  6900846.875,
  6902585.9375,
  6903382.8125,
  6904778.125,
  6905640.625,
  6906429.6875,
  6907373.4375,
  6907781.25,
  6907792.1875,
  6907843.75,
  6908192.1875,
  6908265.625,
  6910601.5625,
  6912381.25,
  6912834.375,
  6914184.375,
  6914898.4375,
  6915556.25,
  6916575.0,
  6916762.5,
  6916864.0625,
  6917006.25,
  6917273.4375,
  6917315.625,
  6917476.5625,
  6918867.1875,
  6920043.75,
  6920695.3125,
  6920700.0,
  6921462.5,
  6923875.0,
  6925793.75,
  6925960.9375,
  6926371.875,
  6926898.4375,
  6929328.125,
  6931801.5625,
  6932778.125,
  6935623.4375,
  6936379.6875,
  6952595.3125,
  6954873.4375,
  6958770.3125,
  6960568.75,
  6960575.0,
  6961759.375,
  6963076.5625,
  6963746.875,
  6965035.9375,
  6968595.3125,
  6968618.75,
  6970173.4375,
  6970673.4375,
  6973121.875,
  6975589.0625,
  6975589.0625,
  6979523.4375,
  6986015.625,
  6990473.4375,
  6992821.875,
  6992875.0,
  6995557.8125,
  7000578.125,
  7003475.0,
  7009964.0625,
  7010067.1875,
  7012059.375,
  7013076.5625,
  7017481.25,
  7019320.3125,
  7019371.875,
  7021131.25,
  7023870.3125,
  7025142.1875,
  7025893.75,
  7029600.0,
  7029676.5625,
  7030925.0,
  7036026.5625,
  7040059.375,
  7045996.875,
  7050882.8125,
  7051054.6875,
  7059025.0,
  7060748.4375,
  7065354.6875,
  7065585.9375,
  7065696.875,
  7066231.25,
  7066287.5,
  7066312.5,
  7066521.875,
  7066557.8125,
  7066704.6875,
  7067154.6875,
  7068790.625,
  7069056.25,
  7069173.4375,
  7069310.9375,
  7069462.5,
  7069679.6875,
  7069912.5,
  7070106.25,
  7070953.125,
  7071235.9375,
  7071382.8125,
  7071450.0,
  7072067.1875,
  7072995.3125,
  7073070.3125,
  7075135.9375,
  7075754.6875,
  7077654.6875,
  7080306.25,
  7083682.8125,
  7087618.75,
  7089595.3125,
  7089709.375,
  7089978.125,
  7091476.5625,
  7091921.875,
  7092782.8125,
  7093464.0625,
  7093693.75,
  7097728.125,
  7097762.5,
  7099778.125,
  7100476.5625,
  7101195.3125,
  7101676.5625,
  7110992.1875,
  7114667.1875,
  7114687.5,
  7135595.3125,
  7135773.4375,
  7136393.75,
  7136567.1875,
  7137306.25,
  7138226.5625,
  7139157.8125,
  7140085.9375,
  7140645.3125,
  7140681.25,
  7141696.875,
  7143968.75,
  7144054.6875,
  7144100.0,
  7144278.125,
  7144456.25,
  7144517.1875,
  7144782.8125,
  7145368.75,
  7146046.875,
  7146190.625,
  7146235.9375,
  7146409.375,
  7147031.25,
  7147806.25,
  7148392.1875,
  7153271.875,
  7153467.1875,
  7154043.75,
  7154815.625,
  7154976.5625,
  7155234.375,
  7155742.1875,
  7155750.0,
  7155759.375,
  7155807.8125,
  7157684.375,
  7160137.5,
  7163331.25,
  7166090.625,
  7170282.8125,
  7178064.0625,
  7181048.4375,
  7181742.1875,
  7181750.0,
  7185629.6875,
  7188131.25,
  7190898.4375,
  7192676.5625,
  7194932.8125,
  7199528.125,
  7200376.5625,
  7202873.4375,
  7206435.9375,
  7233745.3125,
  7238456.25,
  7240337.5,
  7245131.25,
  7245793.75,
  7246970.3125,
  7267000.0,
  7268742.1875,
  7270745.3125,
  7272240.625,
  7274007.8125,
  7274582.8125,
  7277534.375,
  7278976.5625,
  7282473.4375,
  7291060.9375,
  7298778.125,
  7351706.25,
  7360393.75,
  7363103.125,
  7363831.25,
  7367239.0625,
  7374210.9375,
  7379806.25,
  7381143.75,
  7386112.5,
  7390381.25,
  7394418.75,
  7395090.625,
  7395090.625,
  7406987.5,
  7411218.75,
  7420918.75,
  7421657.8125,
  7431370.3125,
  7431389.0625,
  7438256.25,
  7446681.25,
  7449803.125,
  7451342.1875,
  7465565.625,
  7502345.3125,
  7506750.0,
  7509904.6875,
  7517221.875,
  7525262.5,
  7538918.75,
  7548454.6875,
  7548785.9375,
  7550428.125,
  7550450.0,
  7551910.9375,
  7564014.0625,
  7577807.8125,
  7601307.8125,
  7604479.6875,
  7606565.625,
  7608357.8125,
  7609064.0625,
  7629371.875,
  7640804.6875,
  7644239.0625,
  7645853.125,
  7659410.9375,
  7663459.375,
  7667162.5,
  7685168.75,
  7686807.8125,
  7688215.625,
  7688368.75,
  7690368.75,
  7788470.3125,
  7791253.125,
  7796962.5,
  7800895.3125,
  7817806.25,
  7866092.1875,
  7881675.0,
  7901190.625,
  7902153.125,
  7903989.0625,
  7904062.5,
  7905389.0625,
  7907606.25,
  7909485.9375,
  7912476.5625,
  7912518.75,
  7915131.25,
  7949098.4375,
  7967892.1875,
  7967943.75,
  7969754.6875,
  7973856.25,
  7977109.375,
  7979151.5625,
  7984501.5625,
  8003348.4375,
  8033926.5625,
  8041746.875,
  8042765.625,
  8218093.75,
  8222915.625,
  8225373.4375,
  8232354.6875,
  8254564.0625,
  8257664.0625,
  8264907.8125,
  8267310.9375,
  8268610.9375,
  8268626.5625,
  8270878.125,
  8270910.9375,
  8273314.0625,
  8274415.625,
  8277448.4375,
  8277653.125,
  8278329.6875,
  8374012.5,
  8402148.4375,
  8406387.5,
  8409414.0625,
  8409714.0625,
  8465943.75,
  8577912.5,
  8577976.5625,
  8582915.625,
  8590328.125,
  8710014.0625,
  8711340.625,
  8711365.625,
  8715759.375,
  8806160.9375,
  8810582.8125,
  8858293.75,
  8950881.25,
  8975806.25,
  8981576.5625,
  8994664.0625,
  9012298.4375,
  9018550.0,
  9018968.75,
  9025235.9375,
  9057146.875,
  9057575.0,
  9167290.625,
  9167465.625,
  9167543.75,
  9167726.5625,
  9280321.875,
  9284970.3125,
  9307468.75,
  9319876.5625,
  9323296.875,
  9332540.625,
  9358239.0625,
  9368064.0625,
  9369412.5,
  9388545.3125,
  9390225.0,
  9391110.9375,
  9393790.625,
  9394431.25,
  9546664.0625,
  9549668.75,
  9550220.3125,
  9551631.25,
  9551732.8125,
  9552384.375,
  9552435.9375,
  9554171.875,
  9556809.375,
  9604979.6875,
  9704865.625,
  9750707.8125,
  9763701.5625,
  9842953.125,
  9858634.375,
  9883668.75,
  9892403.125,
  9898618.75,
  9905582.8125,
  9931401.5625,
  9935009.375,
  9971907.8125,
  10094093.75,
  10104489.0625,
  10105764.0625,
  10119478.125,
  10141587.5,
  10142315.625,
  10148984.375,
  10260142.1875,
  10262071.875,
  10976793.75,
  10985671.875,
  11081867.1875,
  11085589.0625,
  11088032.8125,
  11089739.0625,
  11090273.4375,
  11090878.125,
  11090971.875,
  11091028.125,
  11091190.625,
  11091320.3125,
  11091643.75,
  11091657.8125,
  11091701.5625,
  11093714.0625,
  11094746.875,
  11094756.25,
  11094843.75,
  11094887.5,
  11095435.9375,
  11095676.5625,
  11095984.375,
  11097200.0,
  11097270.3125,
  11098148.4375,
  11098903.125,
  11099368.75,
  11099498.4375,
  11099810.9375,
  11099860.9375,
  11099898.4375,
  11099960.9375,
  11099976.5625,
  11100250.0,
  11100485.9375,
  11100831.25,
  11101481.25,
  11101570.3125,
  11101578.125,
  11102368.75,
  11102415.625,
  11102873.4375,
  11104481.25,
  11106231.25,
  11108228.125,
  11108443.75,
  11109021.875,
  11110164.0625,
  11111443.75,
  11111971.875,
  11112696.875,
  11114500.0,
  11115729.6875,
  11116270.3125,
  11116929.6875,
  11117157.8125,
  11117815.625,
  11117992.1875,
  11120812.5,
  11120823.4375,
  11122151.5625,
  11126929.6875,
  11127943.75,
  11143981.25,
  11145668.75,
  11147581.25],
 [23.09646681816278,
  76.41414913220474,
  25.075078733203807,
  17.898604521814896,
  78.76138453820957,
  147.54782317156085,
  43.84313436386716,
  7.833122647288747,
  60.29336420415457,
  5.44352088317474,
  6.935434419250196,
  54.49208325453253,
  83.42235272965655,
  9.755770621741968,
  13.464362350621332,
  37.383263790445454,
  9.751771452360021,
  7.773770450052657,
  46.7573095722054,
  10.430067941002834,
  39.55250744823101,
  17.758857026542117,
  40.68133884254476,
  33.74940800963735,
  30.580128199506913,
  16.964085514229957,
  15.004966896612386,
  245.18515574827353,
  5.890444415388376,
  19.136731114907494,
  27.68792197258647,
  5.596824457932807,
  39.000931588622926,
  74.6455818728687,
  6.765714189037943,
  47.51222106614824,
  12.06725468525347,
  16.677661528556477,
  16.16999994114641,
  8.062978836147936,
  46.162819240992185,
  156.74331799869452,
  37.71091157723907,
  56.48503749092269,
  10.869049450635831,
  24.166134979983333,
  52.34594265703679,
  11.252806776483052,
  39.815794913069055,
  41.64119877348424,
  35.1740284401459,
  27.870468537171902,
  7.306225789525635,
  25.447705441255916,
  37.29166752073966,
  34.10814381034519,
  31.114001235457323,
  12.941470497674922,
  77.27955928476669,
  108.98125412119299,
  5.3398275472173165,
  18.527871245900037,
  38.805576006593455,
  79.99212316262845,
  35.516098161421354,
  85.50813448089428,
  84.08170406323603,
  24.55550101636636,
  22.58398745702308,
  24.10580114876118,
  8.048760514120213,
  59.26493170503068,
  39.50766036043862,
  45.987494532412576,
  55.63748291742418,
  72.91849444647515,
  95.68263776290124,
  86.77436429018091,
  24.685054403488337,
  94.61352608197888,
  53.57752154467659,
  83.65786433103156,
  81.46964795179817,
  11.8095624608373,
  54.5482839923181,
  52.745290908921895,
  5.52930637892141,
  14.481070457927395,
  8.65333440859475,
  130.18674058854342,
  107.31760502008677,
  55.66153402761941,
  21.4402597057734,
  37.45138916764656,
  11.450338501938521,
  5.600274801692124,
  45.03289155919002,
  13.019672836360414,
  14.508268209836615,
  70.09213361995104,
  7.699263183584653,
  88.89142974947379,
  34.300862620141665,
  38.14559947008174,
  27.92221445976094,
  29.48743251472165,
  88.71025095302657,
  19.077015125877832,
  10.310533446307305,
  12.778054730716493,
  23.816304618104432,
  25.20174802558589,
  105.84118914771959,
  11.18597591530499,
  6.360335623868345,
  13.680083523795059,
  26.010808352740497,
  21.435728426955674,
  20.24348663800454,
  7.287026655368721,
  18.97657204196281,
  25.090317154221403,
  122.50238002421648,
  28.209934201111654,
  13.250750964895458,
  158.27526495782521,
  106.32741086267562,
  53.27582112658458,
  18.821001625418006,
  21.80317910608514,
  13.709757200410335,
  34.35407720956256,
  63.463368243833386,
  60.746622542136336,
  20.58591512859365,
  49.542736088842595,
  8.227872078351501,
  88.09822791070869,
  9.881491038395996,
  53.702460440828695,
  5.3447145148315665,
  106.1402034063197,
  94.95346518361814,
  54.03440814621889,
  83.66897597288896,
  96.81380570974585,
  6.673972891229778,
  41.65819208481658,
  14.549226478082966,
  11.548226351795945,
  47.0255136529825,
  16.328130398289087,
  23.87731238925581,
  115.73847989490514,
  90.61983541701595,
  40.20101668185603,
  117.93539228475228,
  55.79221911734366,
  61.181829751242844,
  31.960025575937426,
  7.3531650387301255,
  34.3891722865601,
  147.2563781752445,
  51.429841456300935,
  16.63375120613278,
  6.0477098742011215,
  80.78014513560368,
  99.4661951827968,
  25.38416878533699,
  58.47022260621292,
  330.2268603532193,
  49.67101210244098,
  81.81598277577862,
  22.249064939666027,
  26.51450836481815,
  176.52379229725503,
  80.31796248175466,
  37.34719957276074,
  14.979334571526724,
  100.93020368272109,
  13.163255659860852,
  16.21049502265484,
  12.722743943066144,
  71.87459248461144,
  10.581285617606289,
  14.807699984013801,
  8.13267345698736,
  131.14444740991922,
  20.491821853036573,
  18.225103235688298,
  22.372282947273654,
  7.468647539550206,
  8.473255416139036,
  7.830244444720974,
  58.25810435812028,
  71.42919666632234,
  8.02490941236178,
  53.131349293729045,
  158.97531573670602,
  109.23348209066081,
  12.68398497328271,
  11.855992948508753,
  17.23002872297242,
  35.56269833299826,
  27.91433198748642,
  66.14264690405385,
  76.36977542702137,
  16.744557449002546,
  75.00794447623903,
  7.471494869889003,
  7.516320133033916,
  7.641250526344109,
  41.01519066018683,
  28.59302217934698,
  27.135047561329532,
  15.490584425037806,
  38.41889344663898,
  10.474118457382808,
  63.1097333017173,
  13.771638454431148,
  196.35440403837347,
  45.33554632990689,
  43.170288930191234,
  30.50494422418303,
  65.11519956445282,
  89.87990784332499,
  68.59200654282664,
  69.17145479141683,
  30.60723710460382,
  19.216705459006064,
  55.33718636806174,
  21.904487595471856,
  85.29922293685348,
  11.306723591833897,
  80.83224547548106,
  27.27600346659839,
  43.283596900698974,
  19.52004102043005,
  66.86170386999859,
  18.55791504706697,
  44.71218593961408,
  76.15520357708716,
  180.40320086682667,
  57.061742459971526,
  27.6705878581027,
  68.33743480511626,
  18.620400431268088,
  24.530800057901743,
  45.41470309582157,
  58.47919510094424,
  25.125724069295355,
  60.330054943416016,
  24.995603050199712,
  17.815204255937715,
  27.004487503945562,
  72.27093023267976,
  27.9191763044152,
  49.5886643752733,
  90.17286374720061,
  91.15992024378313,
  103.96258110876433,
  5.282378293816456,
  43.65068488987147,
  80.06988595427131,
  12.08540678604747,
  63.60527539957998,
  13.278586783970654,
  14.36800305439691,
  18.69044990012204,
  9.038911259198743,
  5.043591136343069,
  60.89134810469868,
  126.12918335024405,
  47.29182710158029,
  11.256055453485402,
  10.08873642581787,
  17.177918251110125,
  11.551871204318578,
  79.49388357412488,
  46.83315449200037,
  91.81541689983051,
  57.182059266192965,
  38.18307005211219,
  8.573750701242474,
  68.7897489082834,
  49.55515559462687,
  54.6617441482015,
  13.684773019959055,
  8.597130172921673,
  56.43386582842407,
  121.76610873983014,
  93.3820218085413,
  114.77288350816902,
  33.57555467030754,
  68.70577305279268,
  126.18799902077353,
  11.806736893432724,
  67.5717274962382,
  31.09733312297488,
  24.9993178571767,
  50.91404483677396,
  74.99383888751169,
  17.164033408826647,
  53.625767602512845,
  5.077356712305322,
  17.048460788630074,
  24.490395523338314,
  32.51869845102124,
  7.40936319877423,
  14.008185744776586,
  58.83469407136514,
  15.96515050247702,
  140.4294180233351,
  15.00029204586558,
  58.959776751222755,
  29.00786598454628,
  22.385333494226465,
  68.08628641581663,
  88.3513608283919,
  57.52898471488527,
  9.170096934300823,
  59.202925898044924,
  10.383975078825948,
  15.908029817261188,
  68.44877639548497,
  13.964203147775955,
  38.12228485606426,
  14.657511901910114,
  55.24969493735518,
  28.418871179080792,
  13.634985284050115,
  20.593114874826995,
  7.268482365145665,
  13.79962092497581,
  8.765552540738394,
  7.052060997673832,
  83.56842636929363,
  44.714391618498134,
  38.14198033587367,
  14.593740943570765,
  14.15180293486147,
  126.14168179685748,
  47.96128619383243,
  26.788844991694212,
  34.126671958694054,
  138.51132064101247,
  33.89650422085187,
  139.78208201277732,
  17.92959290678197,
  14.058029814699449,
  11.700523940850971,
  63.92208444707154,
  68.67228810727593,
  67.17039380611489,
  82.80194814505042,
  65.55138104825598,
  159.2600351333081,
  15.186617889429188,
  19.564137007425693,
  16.56356673971569,
  7.4023145860666855,
  66.03785765602143,
  6.717644962731825,
  56.43136262214742,
  176.94409330818468,
  59.981433873743,
  89.21498829623785,
  21.48360170085487,
  6.803592364634948,
  15.452444273169021,
  67.33265802241313,
  75.33599995109903,
  29.369952890623452,
  53.35715545937019,
  33.1669456365453,
  88.0331667248147,
  82.28566960822343,
  72.11796308355532,
  22.64400155187669,
  19.100835933485165,
  16.23567369086813,
  11.608350448985593,
  7.232920380501984,
  48.77432665652317,
  18.870389674506736,
  26.612544181833155,
  76.42600867198205,
  14.759024019388667,
  31.133951694881247,
  23.202970867937008,
  8.68879198888946,
  112.10675524404846,
  5.249294819257113,
  73.16996781306247,
  10.088316178286297,
  12.626010192571085,
  71.06954982983893,
  96.4113934677713,
  25.716591665227995,
  74.70676345939859,
  15.061637140305567,
  31.349563703591844,
  14.765575628876116,
  41.37783741683598,
  6.607019620649585,
  20.016470084221307,
  21.39283167792491,
  54.415168874197455,
  89.48726452015974,
  8.974007576842856,
  10.37863948876639,
  10.010938030597671,
  160.44211687263984,
  5.425689625397781,
  5.988969731469766,
  14.868070754044982,
  16.703555716568903,
  5.032782901148494,
  25.747973838410505,
  12.732413670650539,
  175.34566492353284,
  23.424524974862887,
  32.0514130049495,
  91.65272460705333,
  39.36689932641691,
  14.38267800892593,
  26.6598596835178,
  33.63610288230363,
  98.20603715132502,
  65.26401332348384,
  313.33811760509525,
  61.039257054260936,
  89.77645657215288,
  19.995126261665384,
  38.79491955904897,
  26.78565655374219,
  8.227283651834314,
  15.077980884520434,
  50.93695406713049,
  56.412173704381,
  5.379946272545311,
  42.65603231855647,
  9.03805733257106,
  16.795191646795338,
  25.13183788510992,
  181.1554443409089,
  56.27861656223156,
  25.60811747633166,
  5.848399199904395,
  11.124725613811325,
  23.88846446546136,
  29.68447464326836,
  16.284857881132073,
  53.90884739056179,
  5.12946536475925,
  9.087724984729412,
  31.370536497249784,
  67.21793660523171,
  11.087284574123233,
  69.57536470147195,
  10.71050659714387,
  15.026406033500642,
  114.17159311802381,
  29.037351510597386,
  19.881878828036456,
  78.2986779092786,
  24.661532508463115,
  156.115212201316,
  60.461577464050656,
  46.51681115564065,
  64.58622582448633,
  9.989703034361797,
  53.15762941177178,
  85.27307199795871,
  176.305296323644,
  31.837548026754792,
  22.34985351742204,
  43.37568616079928,
  6.460963445932313,
  31.50834246245548,
  20.58552746852196,
  44.94705353616304,
  37.17461838773451,
  5.099194077332632,
  32.90993039607598,
  104.9916000738912,
  90.29404954860809,
  21.284921249870045,
  63.362470157772734,
  61.21107144910174,
  39.99791073761996,
  29.648512739081998,
  115.00463624065695,
  36.22922519753022,
  42.537829192352824,
  64.45626154113465,
  54.51225385741079,
  80.82654142283819,
  56.97474121039772,
  6.2683575364450475,
  6.3022993995514724,
  34.43252960671323,
  17.084033084814408,
  11.510704211004917,
  13.846648250539005,
  25.24238478363363,
  51.31691896485763,
  10.627951732475035,
  51.23120075689179,
  13.937854078563294,
  20.77498925745532,
  90.58011909951925,
  86.28657754369108,
  28.166031570200573,
  43.34354334706103,
  63.755564071280325,
  8.440272586614352,
  47.423728401763235,
  15.60628441250115,
  6.462352017126078,
  46.34345583843195,
  42.644513503265884,
  67.05489588593181,
  57.700090619293356,
  50.23355206617134,
  29.93019455834159,
  42.323728366451064,
  324.6762747630948,
  6.62776076943185,
  90.55236982052737,
  5.192636666306374,
  74.55620479146661,
  29.349448639369495,
  26.334210243547137,
  90.41760447507804,
  14.649971147491229,
  6.477199593589102,
  59.8287600317738,
  47.55215898458397,
  208.2547132477722,
  11.25065545428929,
  7.618536229633273,
  29.21729821751388,
  5.396782012222748,
  12.447092500610562,
  24.526174011643924,
  102.91069207739035,
  76.19809814420171,
  67.7368598307686,
  66.39835136692548,
  28.66129661975702,
  15.87962279485578,
  26.23462793271085,
  40.88067496199566,
  9.299442518713496,
  15.304711307160304,
  15.15215809710493,
  160.08728455182893,
  32.315523665824294,
  14.641212781246002,
  22.689494937996766,
  51.39737824257823,
  24.507053897445232,
  6.615755445004151,
  41.53936272775266,
  32.64512891996418,
  6.499348054041406,
  39.956518034349926,
  55.220384188716416,
  7.855200405854162,
  26.23293667230579,
  94.02840278444549,
  137.5732317944795,
  11.74106718209256,
  53.216812116956646,
  9.121092394572802,
  26.978770443785013,
  23.976403344690638,
  131.13900161001277,
  5.512709807354648,
  83.50587695776316,
  120.75754111154025,
  21.432627649464347,
  15.938218510478407,
  5.50792898869451,
  8.349560441435239,
  46.04330878224397,
  42.5215676405118,
  204.2566448446522,
  15.540972179187387,
  21.829801157422736,
  75.73574552508313,
  30.100486366045565,
  18.725531276530628,
  63.65354681442656,
  179.41647082102958,
  71.29779825656281,
  18.026845898146256,
  50.580100325133216,
  25.374202899652573,
  117.99260613097047,
  29.569000140727198,
  83.50072556539898,
  20.17129171681023,
  5.152129416571529,
  12.418697704452997,
  12.918247060528927,
  117.13813439275386,
  18.818056661350806,
  44.77227416705334,
  10.863722794633293,
  17.257761057982442,
  20.913078397867118,
  46.23452095016029,
  14.500311047434908,
  51.84250759072827,
  26.786096373728334,
  8.854665134358058,
  75.59907570488885,
  5.7794721015084685,
  16.311281090981943,
  44.1448549007064,
  11.90691444290629,
  66.32186209786246,
  8.88679758115977,
  5.3769170970019236,
  15.825348266967726,
  248.7647104338708,
  95.88248605413318,
  42.2519594346983,
  59.32768348766648,
  130.97140295118984,
  58.570285304433995,
  22.916988222888275,
  69.45381813591031,
  39.744556718555536,
  23.630047931521005,
  16.20378871730779,
  49.85046416041436,
  72.04361428474628,
  100.64556626023996,
  104.25108810520656,
  7.702191769034201,
  54.38835104227716,
  14.39549290924738,
  39.19447563840014,
  5.249218109261389,
  8.54203346361912,
  75.08793963255036,
  6.603855332714815,
  48.617068927368884,
  46.24467225265641,
  7.23011467546573,
  5.292418075307008,
  30.24723549342636,
  20.507620622784927,
  45.39737759814344,
  15.11899561294872,
  17.971645091223593,
  9.73349556867445,
  14.198727066882569,
  5.80524972819613,
  35.03792022298571,
  148.5710923952991,
  37.08055407342612,
  5.8737625917498555,
  29.335296710939083,
  45.148205970740364,
  48.20780377183974,
  13.897105916045245,
  105.22865824518311,
  62.41919324789416,
  71.19244883367816,
  18.461070889653755,
  47.67862493246133,
  20.242123879078814,
  10.880988443715635,
  31.141345074159133,
  5.637329293642736,
  23.265373497688167,
  9.93942213205608,
  70.24366398277022,
  70.64072638104335,
  11.627387288026377,
  29.40552434894623,
  11.79661355602266,
  25.17408269921721,
  8.028697841877118,
  95.86410563706423,
  54.15152188844783,
  7.6522072900293345,
  89.27910220284586,
  25.878933171395296,
  113.09481307963549,
  70.26634073319221,
  106.3503912943476,
  17.972447492576638,
  5.385817299348203,
  13.07458698314503,
  14.009872863302448,
  9.615292714087882,
  10.646363787476574,
  17.54145887959419,
  47.361354483070414,
  50.97301100206579,
  34.55933999788362,
  7.997463485011918,
  37.48030616404397,
  19.552018062377808,
  63.41710539045099,
  80.13770821137499,
  9.539568409049592,
  20.72629194435229,
  106.13315063164029,
  8.774496474086412,
  83.06297317549364,
  10.997196061981022,
  9.877489272470132,
  31.168933551461365,
  29.943901691072227,
  12.856009314396326,
  5.38110965223008,
  40.57188921262078,
  27.064820768877112,
  163.84810559956765,
  48.919726586645695,
  6.1713507826946055,
  16.65902539830648,
  148.78754953733747,
  12.420142059606889,
  16.92961780800558,
  38.43958774119711,
  81.94086042832038,
  19.83292741911609,
  53.37252801495728,
  48.740876131409756,
  75.47926639866053,
  139.69361338466956,
  17.578225571721887,
  11.340479901121961,
  6.194325530092807,
  21.090540718314756,
  46.88065943311076,
  58.030598924109796,
  185.99291227923194,
  6.604815087221636,
  11.34559824691836,
  10.77523594242167,
  134.63759007171868,
  22.303078670593184,
  11.462877499945172,
  111.20880094189025,
  49.33848601370531,
  107.42201470004261,
  37.451382911049514,
  29.716720293362197,
  89.94955215673986,
  6.072013284829785,
  28.963532023674382,
  126.91197967481754,
  60.25867905104224,
  20.691830293424896,
  78.29741749805024,
  54.37888426891168,
  16.937010273912406,
  7.352485570199575,
  97.59141527917566,
  54.79435777258159,
  108.27722485514072,
  92.26586728032673,
  11.579327749792661,
  12.431137668497986,
  13.577408393609055,
  40.9028747991367,
  41.16533159228178,
  112.71383667687874,
  11.31823441804757,
  72.9319093451652,
  16.638607330865103,
  10.440241915268214,
  155.98499119284438,
  26.8434550999985,
  32.54190314571381,
  10.519969282958632,
  118.51697667439582,
  29.154088198990237,
  40.72790912104638,
  53.04413414909574,
  60.907404841461506,
  146.38469734345168,
  8.642480698497081,
  71.43106764356533,
  28.402211411925478,
  22.065729394112545,
  68.10017307000676,
  77.00978165031549,
  5.309758889323243,
  27.99778135881106,
  18.661261356024767,
  54.69931075962305,
  33.81723448596374,
  82.41399349664184,
  35.47818686275496,
  19.08572738311426,
  32.41250777936446,
  34.35461860032439,
  24.415892744744156,
  5.13062813839644,
  13.612852055873017,
  52.009854991578194,
  5.580271875296687,
  26.442735264009773,
  5.746522042028261,
  20.328272984819076,
  21.965844840627856,
  11.929193471490276,
  20.79543070707279,
  22.89528470495999,
  12.183623574882127,
  24.3092891425887,
  74.3320937831429,
  28.127506546643282,
  19.483343445062122,
  182.37556629039452,
  32.25122617212112,
  15.600633867317907,
  57.24508359765869,
  8.048476474559743,
  20.242070402730736,
  56.63173754070698,
  5.036582668575584,
  5.152819356371734,
  18.94931258753779,
  10.112816446706383,
  65.31845825774111,
  8.885312023415915,
  8.299693222890182,
  30.474495787567868,
  5.634890417683859,
  15.578056641766537,
  26.38535609145219,
  17.812532591067523,
  15.266133509541646,
  6.314812793814502,
  20.4580238254229,
  33.10733152014791,
  34.98019841220422,
  6.223366604155347,
  26.163352467987494,
  18.177070756327915,
  5.047114330222604,
  17.10749424453423,
  13.416464862116142,
  15.544373022705157,
  63.20163433849206,
  16.11500586051435,
  38.872233052806855,
  7.277746090901515,
  8.757433423977686,
  40.600452373231555,
  7.0815913205529295,
  5.925992373342192,
  24.62537968757753,
  57.854075101505906,
  20.175740873767914,
  28.87070208056962,
  64.47387747461397,
  88.24292921784107,
  14.097615983920917,
  33.13578221402966,
  20.563481770435782,
  27.18087927802395,
  5.606927873887151,
  16.412155775802216,
  6.451790461334643,
  21.935676144929356,
  12.14082664236699,
  22.36388594216339,
  8.340142717409163,
  6.678983011082577,
  117.98343041161252,
  6.693610119808251,
  8.712969274987193,
  41.39149943080047,
  5.562127579715948,
  56.637391453645485,
  88.54937524550088,
  14.837974107430604,
  15.53116409073412,
  16.928024137461136,
  27.77850471320413,
  19.149004144307735,
  22.151958147466146,
  152.45963455967888,
  14.129868939481765,
  7.965612117307012,
  47.785447491745174,
  42.94388536101974,
  5.861437667978293,
  86.69090644946715,
  46.266655787238484,
  14.680345605269244,
  50.56865604680539,
  66.77499383058839,
  26.720263285608375,
  8.972851083138643,
  38.63317139744809,
  21.170853121561734,
  18.170342667268365,
  16.85524444517276,
  76.31174879923557,
  13.18518358807878,
  6.3255234517095325,
  28.336549432247992,
  39.86597397288527,
  19.315255818557585,
  14.917933649184565,
  14.661277350816956,
  5.198887787396429,
  12.353963742261262,
  52.24058896466551,
  13.706679082816027,
  17.293808976705286,
  58.1097379060424,
  12.107441245541743,
  10.174894450405562,
  41.79909314140049,
  33.18711371552374,
  6.532810932434207,
  5.783015978335344,
  49.658516952763996,
  5.1686063468196695,
  71.031184499903,
  18.244887258021098,
  70.92163024034326,
  34.84218192673624,
  120.18083984850456,
  15.594978685148265,
  9.485102212893365,
  15.134936863091614,
  10.859802037020275,
  8.687640683123133,
  84.05623191754816,
  26.995388524964685,
  23.41861402244887,
  55.981966176233556,
  27.94251376559162,
  16.093502188542093,
  17.552646316599713,
  19.87104944265511,
  50.110761261535,
  99.36883421197973,
  6.525270257340417,
  6.794619366257612,
  16.753376032615233,
  16.055384113524283,
  56.36045723468583,
  63.768470484820625,
  62.894571415145165,
  16.63593416169257,
  26.267111242820263,
  14.683615690575087,
  53.32980836243981,
  15.647127897861955,
  11.060981384999367,
  6.573082482109846])

Count of hits during the discharge

caption, T_hit,count = 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, count, caption);

([5678503.125,
  5860962.5,
  5872554.6875,
  5884178.125,
  5888800.0,
  6060045.3125,
  6062376.5625,
  6062757.8125,
  6072537.5,
  6098928.125,
  6114381.25,
  6149687.5,
  6150795.3125,
  6164462.5,
  6183070.3125,
  6183078.125,
  6183126.5625,
  6183957.8125,
  6187428.125,
  6187665.625,
  6191129.6875,
  6191195.3125,
  6225615.625,
  6288882.8125,
  6289965.625,
  6308125.0,
  6314806.25,
  6319070.3125,
  6320145.3125,
  6322445.3125,
  6323484.375,
  6325187.5,
  6325307.8125,
  6326343.75,
  6327729.6875,
  6328107.8125,
  6328295.3125,
  6351643.75,
  6354384.375,
  6363343.75,
  6368209.375,
  6368776.5625,
  6379029.6875,
  6383468.75,
  6384890.625,
  6387459.375,
  6387607.8125,
  6388387.5,
  6392968.75,
  6401684.375,
  6403248.4375,
  6403882.8125,
  6404075.0,
  6404154.6875,
  6404998.4375,
  6406692.1875,
  6406868.75,
  6407050.0,
  6407400.0,
  6408132.8125,
  6408265.625,
  6408631.25,
  6408812.5,
  6413734.375,
  6414057.8125,
  6415409.375,
  6415915.625,
  6415917.1875,
  6418026.5625,
  6431112.5,
  6445409.375,
  6449481.25,
  6451356.25,
  6458643.75,
  6495289.0625,
  6497846.875,
  6498734.375,
  6499989.0625,
  6500110.9375,
  6501004.6875,
  6503451.5625,
  6503509.375,
  6503714.0625,
  6505385.9375,
  6507204.6875,
  6509682.8125,
  6513817.1875,
  6514254.6875,
  6514712.5,
  6518403.125,
  6518662.5,
  6526100.0,
  6528420.3125,
  6528634.375,
  6528989.0625,
  6529520.3125,
  6531948.4375,
  6536903.125,
  6536912.5,
  6539871.875,
  6541196.875,
  6544579.6875,
  6547056.25,
  6554398.4375,
  6560371.875,
  6560704.6875,
  6562389.0625,
  6562392.1875,
  6562395.3125,
  6564975.0,
  6566864.0625,
  6567648.4375,
  6569500.0,
  6572759.375,
  6575773.4375,
  6581865.625,
  6585104.6875,
  6585753.125,
  6587948.4375,
  6588715.625,
  6593353.125,
  6602042.1875,
  6604985.9375,
  6606676.5625,
  6609346.875,
  6609371.875,
  6610631.25,
  6611006.25,
  6611206.25,
  6611212.5,
  6611412.5,
  6611510.9375,
  6611664.0625,
  6611909.375,
  6612621.875,
  6612631.25,
  6613570.3125,
  6614082.8125,
  6614246.875,
  6614270.3125,
  6615162.5,
  6615198.4375,
  6615228.125,
  6615229.6875,
  6615293.75,
  6615837.5,
  6616484.375,
  6616518.75,
  6617845.3125,
  6618192.1875,
  6618314.0625,
  6618432.8125,
  6618553.125,
  6618834.375,
  6619000.0,
  6619382.8125,
  6619385.9375,
  6619445.3125,
  6619698.4375,
  6619957.8125,
  6619964.0625,
  6620306.25,
  6620470.3125,
  6620510.9375,
  6620579.6875,
  6620581.25,
  6620750.0,
  6621093.75,
  6621170.3125,
  6621343.75,
  6621367.1875,
  6621379.6875,
  6621400.0,
  6621406.25,
  6621492.1875,
  6621885.9375,
  6621889.0625,
  6622054.6875,
  6622654.6875,
  6622775.0,
  6623178.125,
  6623179.6875,
  6623465.625,
  6623529.6875,
  6623551.5625,
  6623648.4375,
  6623704.6875,
  6624196.875,
  6624471.875,
  6624584.375,
  6624640.625,
  6624750.0,
  6624975.0,
  6625257.8125,
  6625257.8125,
  6625315.625,
  6625353.125,
  6625510.9375,
  6625515.625,
  6625965.625,
  6626201.5625,
  6626478.125,
  6626506.25,
  6626592.1875,
  6626923.4375,
  6627325.0,
  6627342.1875,
  6627926.5625,
  6628662.5,
  6628750.0,
  6629064.0625,
  6629412.5,
  6629768.75,
  6630779.6875,
  6631170.3125,
  6631545.3125,
  6634039.0625,
  6634050.0,
  6634210.9375,
  6635254.6875,
  6635893.75,
  6636296.875,
  6636431.25,
  6637026.5625,
  6637601.5625,
  6638598.4375,
  6639932.8125,
  6646700.0,
  6646770.3125,
  6654910.9375,
  6655240.625,
  6655378.125,
  6656523.4375,
  6658023.4375,
  6662957.8125,
  6663925.0,
  6664485.9375,
  6668675.0,
  6670096.875,
  6672360.9375,
  6674323.4375,
  6675815.625,
  6676489.0625,
  6678357.8125,
  6678648.4375,
  6680739.0625,
  6683292.1875,
  6683803.125,
  6687245.3125,
  6687429.6875,
  6688012.5,
  6690857.8125,
  6697565.625,
  6698504.6875,
  6703725.0,
  6705039.0625,
  6710067.1875,
  6710565.625,
  6721034.375,
  6722598.4375,
  6723471.875,
  6723831.25,
  6725176.5625,
  6727451.5625,
  6727943.75,
  6728465.625,
  6730565.625,
  6731046.875,
  6733153.125,
  6733410.9375,
  6733471.875,
  6733742.1875,
  6734029.6875,
  6734193.75,
  6735429.6875,
  6735775.0,
  6735975.0,
  6736114.0625,
  6737390.625,
  6737484.375,
  6737603.125,
  6738626.5625,
  6739487.5,
  6739673.4375,
  6740064.0625,
  6740576.5625,
  6741485.9375,
  6741565.625,
  6742414.0625,
  6742414.0625,
  6742623.4375,
  6742801.5625,
  6743048.4375,
  6743175.0,
  6744157.8125,
  6744782.8125,
  6745332.8125,
  6745620.3125,
  6746023.4375,
  6746356.25,
  6747031.25,
  6747898.4375,
  6748695.3125,
  6749314.0625,
  6749348.4375,
  6749381.25,
  6750228.125,
  6752625.0,
  6753278.125,
  6754379.6875,
  6757500.0,
  6758064.0625,
  6758256.25,
  6760420.3125,
  6760815.625,
  6760851.5625,
  6761587.5,
  6763392.1875,
  6763410.9375,
  6764268.75,
  6765159.375,
  6766542.1875,
  6767042.1875,
  6769521.875,
  6777809.375,
  6779679.6875,
  6780801.5625,
  6780850.0,
  6781104.6875,
  6782050.0,
  6782209.375,
  6782543.75,
  6782792.1875,
  6783254.6875,
  6783643.75,
  6784407.8125,
  6786764.0625,
  6786784.375,
  6786982.8125,
  6787457.8125,
  6787948.4375,
  6788417.1875,
  6788537.5,
  6788556.25,
  6788739.0625,
  6788881.25,
  6788934.375,
  6789032.8125,
  6789140.625,
  6789268.75,
  6789973.4375,
  6790004.6875,
  6790459.375,
  6790929.6875,
  6791135.9375,
  6791154.6875,
  6791239.0625,
  6792198.4375,
  6792254.6875,
  6792759.375,
  6792898.4375,
  6793490.625,
  6793732.8125,
  6793926.5625,
  6794034.375,
  6794114.0625,
  6794187.5,
  6794189.0625,
  6794207.8125,
  6794559.375,
  6794575.0,
  6794739.0625,
  6794751.5625,
  6795009.375,
  6795168.75,
  6795190.625,
  6795210.9375,
  6795228.125,
  6795243.75,
  6795290.625,
  6795678.125,
  6795764.0625,
  6796340.625,
  6796482.8125,
  6796726.5625,
  6797367.1875,
  6797665.625,
  6797901.5625,
  6798106.25,
  6798120.3125,
  6798428.125,
  6798550.0,
  6798664.0625,
  6798951.5625,
  6799101.5625,
  6799198.4375,
  6799215.625,
  6799923.4375,
  6800028.125,
  6800048.4375,
  6800150.0,
  6800357.8125,
  6800451.5625,
  6800528.125,
  6800873.4375,
  6800907.8125,
  6801228.125,
  6801301.5625,
  6801434.375,
  6801617.1875,
  6801932.8125,
  6802395.3125,
  6802421.875,
  6802481.25,
  6802651.5625,
  6802707.8125,
  6802732.8125,
  6802787.5,
  6802820.3125,
  6802964.0625,
  6803150.0,
  6803214.0625,
  6803259.375,
  6803351.5625,
  6803457.8125,
  6803504.6875,
  6803598.4375,
  6803653.125,
  6803689.0625,
  6803703.125,
  6803775.0,
  6803798.4375,
  6804164.0625,
  6804168.75,
  6804193.75,
  6804204.6875,
  6804207.8125,
  6804210.9375,
  6804214.0625,
  6804220.3125,
  6804248.4375,
  6804256.25,
  6804259.375,
  6804512.5,
  6804514.0625,
  6804543.75,
  6804635.9375,
  6804653.125,
  6804793.75,
  6804853.125,
  6804917.1875,
  6804945.3125,
  6805145.3125,
  6805218.75,
  6805304.6875,
  6805348.4375,
  6805385.9375,
  6805387.5,
  6805501.5625,
  6805539.0625,
  6805578.125,
  6805607.8125,
  6805720.3125,
  6805750.0,
  6805773.4375,
  6805773.4375,
  6805901.5625,
  6805937.5,
  6806021.875,
  6806051.5625,
  6806189.0625,
  6806209.375,
  6806298.4375,
  6806320.3125,
  6806321.875,
  6806618.75,
  6807000.0,
  6807101.5625,
  6807142.1875,
  6807285.9375,
  6807290.625,
  6807475.0,
  6807482.8125,
  6807551.5625,
  6807798.4375,
  6807928.125,
  6807979.6875,
  6808868.75,
  6809035.9375,
  6809103.125,
  6809131.25,
  6809871.875,
  6810701.5625,
  6810985.9375,
  6812331.25,
  6814312.5,
  6841392.1875,
  6842440.625,
  6842443.75,
  6843315.625,
  6843464.0625,
  6843492.1875,
  6846635.9375,
  6847753.125,
  6847776.5625,
  6876957.8125,
  6878031.25,
  6878120.3125,
  6879685.9375,
  6881803.125,
  6882735.9375,
  6888564.0625,
  6888709.375,
  6891264.0625,
  6892821.875,
  6893875.0,
  6895110.9375,
  6898823.4375,
  6898892.1875,
  6900846.875,
  6902585.9375,
  6903382.8125,
  6904778.125,
  6905640.625,
  6906429.6875,
  6907373.4375,
  6907781.25,
  6907792.1875,
  6907843.75,
  6908192.1875,
  6908265.625,
  6910601.5625,
  6912381.25,
  6912834.375,
  6914184.375,
  6914898.4375,
  6915556.25,
  6916575.0,
  6916762.5,
  6916864.0625,
  6917006.25,
  6917273.4375,
  6917315.625,
  6917476.5625,
  6918867.1875,
  6920043.75,
  6920695.3125,
  6920700.0,
  6921462.5,
  6923875.0,
  6925793.75,
  6925960.9375,
  6926371.875,
  6926898.4375,
  6929328.125,
  6931801.5625,
  6932778.125,
  6935623.4375,
  6936379.6875,
  6952595.3125,
  6954873.4375,
  6958770.3125,
  6960568.75,
  6960575.0,
  6961759.375,
  6963076.5625,
  6963746.875,
  6965035.9375,
  6968595.3125,
  6968618.75,
  6970173.4375,
  6970673.4375,
  6973121.875,
  6975589.0625,
  6975589.0625,
  6979523.4375,
  6986015.625,
  6990473.4375,
  6992821.875,
  6992875.0,
  6995557.8125,
  7000578.125,
  7003475.0,
  7009964.0625,
  7010067.1875,
  7012059.375,
  7013076.5625,
  7017481.25,
  7019320.3125,
  7019371.875,
  7021131.25,
  7023870.3125,
  7025142.1875,
  7025893.75,
  7029600.0,
  7029676.5625,
  7030925.0,
  7036026.5625,
  7040059.375,
  7045996.875,
  7050882.8125,
  7051054.6875,
  7059025.0,
  7060748.4375,
  7065354.6875,
  7065585.9375,
  7065696.875,
  7066231.25,
  7066287.5,
  7066312.5,
  7066521.875,
  7066557.8125,
  7066704.6875,
  7067154.6875,
  7068790.625,
  7069056.25,
  7069173.4375,
  7069310.9375,
  7069462.5,
  7069679.6875,
  7069912.5,
  7070106.25,
  7070953.125,
  7071235.9375,
  7071382.8125,
  7071450.0,
  7072067.1875,
  7072995.3125,
  7073070.3125,
  7075135.9375,
  7075754.6875,
  7077654.6875,
  7080306.25,
  7083682.8125,
  7087618.75,
  7089595.3125,
  7089709.375,
  7089978.125,
  7091476.5625,
  7091921.875,
  7092782.8125,
  7093464.0625,
  7093693.75,
  7097728.125,
  7097762.5,
  7099778.125,
  7100476.5625,
  7101195.3125,
  7101676.5625,
  7110992.1875,
  7114667.1875,
  7114687.5,
  7135595.3125,
  7135773.4375,
  7136393.75,
  7136567.1875,
  7137306.25,
  7138226.5625,
  7139157.8125,
  7140085.9375,
  7140645.3125,
  7140681.25,
  7141696.875,
  7143968.75,
  7144054.6875,
  7144100.0,
  7144278.125,
  7144456.25,
  7144517.1875,
  7144782.8125,
  7145368.75,
  7146046.875,
  7146190.625,
  7146235.9375,
  7146409.375,
  7147031.25,
  7147806.25,
  7148392.1875,
  7153271.875,
  7153467.1875,
  7154043.75,
  7154815.625,
  7154976.5625,
  7155234.375,
  7155742.1875,
  7155750.0,
  7155759.375,
  7155807.8125,
  7157684.375,
  7160137.5,
  7163331.25,
  7166090.625,
  7170282.8125,
  7178064.0625,
  7181048.4375,
  7181742.1875,
  7181750.0,
  7185629.6875,
  7188131.25,
  7190898.4375,
  7192676.5625,
  7194932.8125,
  7199528.125,
  7200376.5625,
  7202873.4375,
  7206435.9375,
  7233745.3125,
  7238456.25,
  7240337.5,
  7245131.25,
  7245793.75,
  7246970.3125,
  7267000.0,
  7268742.1875,
  7270745.3125,
  7272240.625,
  7274007.8125,
  7274582.8125,
  7277534.375,
  7278976.5625,
  7282473.4375,
  7291060.9375,
  7298778.125,
  7351706.25,
  7360393.75,
  7363103.125,
  7363831.25,
  7367239.0625,
  7374210.9375,
  7379806.25,
  7381143.75,
  7386112.5,
  7390381.25,
  7394418.75,
  7395090.625,
  7395090.625,
  7406987.5,
  7411218.75,
  7420918.75,
  7421657.8125,
  7431370.3125,
  7431389.0625,
  7438256.25,
  7446681.25,
  7449803.125,
  7451342.1875,
  7465565.625,
  7502345.3125,
  7506750.0,
  7509904.6875,
  7517221.875,
  7525262.5,
  7538918.75,
  7548454.6875,
  7548785.9375,
  7550428.125,
  7550450.0,
  7551910.9375,
  7564014.0625,
  7577807.8125,
  7601307.8125,
  7604479.6875,
  7606565.625,
  7608357.8125,
  7609064.0625,
  7629371.875,
  7640804.6875,
  7644239.0625,
  7645853.125,
  7659410.9375,
  7663459.375,
  7667162.5,
  7685168.75,
  7686807.8125,
  7688215.625,
  7688368.75,
  7690368.75,
  7788470.3125,
  7791253.125,
  7796962.5,
  7800895.3125,
  7817806.25,
  7866092.1875,
  7881675.0,
  7901190.625,
  7902153.125,
  7903989.0625,
  7904062.5,
  7905389.0625,
  7907606.25,
  7909485.9375,
  7912476.5625,
  7912518.75,
  7915131.25,
  7949098.4375,
  7967892.1875,
  7967943.75,
  7969754.6875,
  7973856.25,
  7977109.375,
  7979151.5625,
  7984501.5625,
  8003348.4375,
  8033926.5625,
  8041746.875,
  8042765.625,
  8218093.75,
  8222915.625,
  8225373.4375,
  8232354.6875,
  8254564.0625,
  8257664.0625,
  8264907.8125,
  8267310.9375,
  8268610.9375,
  8268626.5625,
  8270878.125,
  8270910.9375,
  8273314.0625,
  8274415.625,
  8277448.4375,
  8277653.125,
  8278329.6875,
  8374012.5,
  8402148.4375,
  8406387.5,
  8409414.0625,
  8409714.0625,
  8465943.75,
  8577912.5,
  8577976.5625,
  8582915.625,
  8590328.125,
  8710014.0625,
  8711340.625,
  8711365.625,
  8715759.375,
  8806160.9375,
  8810582.8125,
  8858293.75,
  8950881.25,
  8975806.25,
  8981576.5625,
  8994664.0625,
  9012298.4375,
  9018550.0,
  9018968.75,
  9025235.9375,
  9057146.875,
  9057575.0,
  9167290.625,
  9167465.625,
  9167543.75,
  9167726.5625,
  9280321.875,
  9284970.3125,
  9307468.75,
  9319876.5625,
  9323296.875,
  9332540.625,
  9358239.0625,
  9368064.0625,
  9369412.5,
  9388545.3125,
  9390225.0,
  9391110.9375,
  9393790.625,
  9394431.25,
  9546664.0625,
  9549668.75,
  9550220.3125,
  9551631.25,
  9551732.8125,
  9552384.375,
  9552435.9375,
  9554171.875,
  9556809.375,
  9604979.6875,
  9704865.625,
  9750707.8125,
  9763701.5625,
  9842953.125,
  9858634.375,
  9883668.75,
  9892403.125,
  9898618.75,
  9905582.8125,
  9931401.5625,
  9935009.375,
  9971907.8125,
  10094093.75,
  10104489.0625,
  10105764.0625,
  10119478.125,
  10141587.5,
  10142315.625,
  10148984.375,
  10260142.1875,
  10262071.875,
  10976793.75,
  10985671.875,
  11081867.1875,
  11085589.0625,
  11088032.8125,
  11089739.0625,
  11090273.4375,
  11090878.125,
  11090971.875,
  11091028.125,
  11091190.625,
  11091320.3125,
  11091643.75,
  11091657.8125,
  11091701.5625,
  11093714.0625,
  11094746.875,
  11094756.25,
  11094843.75,
  11094887.5,
  11095435.9375,
  11095676.5625,
  11095984.375,
  11097200.0,
  11097270.3125,
  11098148.4375,
  11098903.125,
  11099368.75,
  11099498.4375,
  11099810.9375,
  11099860.9375,
  11099898.4375,
  11099960.9375,
  11099976.5625,
  11100250.0,
  11100485.9375,
  11100831.25,
  11101481.25,
  11101570.3125,
  11101578.125,
  11102368.75,
  11102415.625,
  11102873.4375,
  11104481.25,
  11106231.25,
  11108228.125,
  11108443.75,
  11109021.875,
  11110164.0625,
  11111443.75,
  11111971.875,
  11112696.875,
  11114500.0,
  11115729.6875,
  11116270.3125,
  11116929.6875,
  11117157.8125,
  11117815.625,
  11117992.1875,
  11120812.5,
  11120823.4375,
  11122151.5625,
  11126929.6875,
  11127943.75,
  11143981.25,
  11145668.75,
  11147581.25],
 [23.09646681816278,
  76.41414913220474,
  25.075078733203807,
  17.898604521814896,
  78.76138453820957,
  147.54782317156085,
  43.84313436386716,
  7.833122647288747,
  60.29336420415457,
  5.44352088317474,
  6.935434419250196,
  54.49208325453253,
  83.42235272965655,
  9.755770621741968,
  13.464362350621332,
  37.383263790445454,
  9.751771452360021,
  7.773770450052657,
  46.7573095722054,
  10.430067941002834,
  39.55250744823101,
  17.758857026542117,
  40.68133884254476,
  33.74940800963735,
  30.580128199506913,
  16.964085514229957,
  15.004966896612386,
  245.18515574827353,
  5.890444415388376,
  19.136731114907494,
  27.68792197258647,
  5.596824457932807,
  39.000931588622926,
  74.6455818728687,
  6.765714189037943,
  47.51222106614824,
  12.06725468525347,
  16.677661528556477,
  16.16999994114641,
  8.062978836147936,
  46.162819240992185,
  156.74331799869452,
  37.71091157723907,
  56.48503749092269,
  10.869049450635831,
  24.166134979983333,
  52.34594265703679,
  11.252806776483052,
  39.815794913069055,
  41.64119877348424,
  35.1740284401459,
  27.870468537171902,
  7.306225789525635,
  25.447705441255916,
  37.29166752073966,
  34.10814381034519,
  31.114001235457323,
  12.941470497674922,
  77.27955928476669,
  108.98125412119299,
  5.3398275472173165,
  18.527871245900037,
  38.805576006593455,
  79.99212316262845,
  35.516098161421354,
  85.50813448089428,
  84.08170406323603,
  24.55550101636636,
  22.58398745702308,
  24.10580114876118,
  8.048760514120213,
  59.26493170503068,
  39.50766036043862,
  45.987494532412576,
  55.63748291742418,
  72.91849444647515,
  95.68263776290124,
  86.77436429018091,
  24.685054403488337,
  94.61352608197888,
  53.57752154467659,
  83.65786433103156,
  81.46964795179817,
  11.8095624608373,
  54.5482839923181,
  52.745290908921895,
  5.52930637892141,
  14.481070457927395,
  8.65333440859475,
  130.18674058854342,
  107.31760502008677,
  55.66153402761941,
  21.4402597057734,
  37.45138916764656,
  11.450338501938521,
  5.600274801692124,
  45.03289155919002,
  13.019672836360414,
  14.508268209836615,
  70.09213361995104,
  7.699263183584653,
  88.89142974947379,
  34.300862620141665,
  38.14559947008174,
  27.92221445976094,
  29.48743251472165,
  88.71025095302657,
  19.077015125877832,
  10.310533446307305,
  12.778054730716493,
  23.816304618104432,
  25.20174802558589,
  105.84118914771959,
  11.18597591530499,
  6.360335623868345,
  13.680083523795059,
  26.010808352740497,
  21.435728426955674,
  20.24348663800454,
  7.287026655368721,
  18.97657204196281,
  25.090317154221403,
  122.50238002421648,
  28.209934201111654,
  13.250750964895458,
  158.27526495782521,
  106.32741086267562,
  53.27582112658458,
  18.821001625418006,
  21.80317910608514,
  13.709757200410335,
  34.35407720956256,
  63.463368243833386,
  60.746622542136336,
  20.58591512859365,
  49.542736088842595,
  8.227872078351501,
  88.09822791070869,
  9.881491038395996,
  53.702460440828695,
  5.3447145148315665,
  106.1402034063197,
  94.95346518361814,
  54.03440814621889,
  83.66897597288896,
  96.81380570974585,
  6.673972891229778,
  41.65819208481658,
  14.549226478082966,
  11.548226351795945,
  47.0255136529825,
  16.328130398289087,
  23.87731238925581,
  115.73847989490514,
  90.61983541701595,
  40.20101668185603,
  117.93539228475228,
  55.79221911734366,
  61.181829751242844,
  31.960025575937426,
  7.3531650387301255,
  34.3891722865601,
  147.2563781752445,
  51.429841456300935,
  16.63375120613278,
  6.0477098742011215,
  80.78014513560368,
  99.4661951827968,
  25.38416878533699,
  58.47022260621292,
  330.2268603532193,
  49.67101210244098,
  81.81598277577862,
  22.249064939666027,
  26.51450836481815,
  176.52379229725503,
  80.31796248175466,
  37.34719957276074,
  14.979334571526724,
  100.93020368272109,
  13.163255659860852,
  16.21049502265484,
  12.722743943066144,
  71.87459248461144,
  10.581285617606289,
  14.807699984013801,
  8.13267345698736,
  131.14444740991922,
  20.491821853036573,
  18.225103235688298,
  22.372282947273654,
  7.468647539550206,
  8.473255416139036,
  7.830244444720974,
  58.25810435812028,
  71.42919666632234,
  8.02490941236178,
  53.131349293729045,
  158.97531573670602,
  109.23348209066081,
  12.68398497328271,
  11.855992948508753,
  17.23002872297242,
  35.56269833299826,
  27.91433198748642,
  66.14264690405385,
  76.36977542702137,
  16.744557449002546,
  75.00794447623903,
  7.471494869889003,
  7.516320133033916,
  7.641250526344109,
  41.01519066018683,
  28.59302217934698,
  27.135047561329532,
  15.490584425037806,
  38.41889344663898,
  10.474118457382808,
  63.1097333017173,
  13.771638454431148,
  196.35440403837347,
  45.33554632990689,
  43.170288930191234,
  30.50494422418303,
  65.11519956445282,
  89.87990784332499,
  68.59200654282664,
  69.17145479141683,
  30.60723710460382,
  19.216705459006064,
  55.33718636806174,
  21.904487595471856,
  85.29922293685348,
  11.306723591833897,
  80.83224547548106,
  27.27600346659839,
  43.283596900698974,
  19.52004102043005,
  66.86170386999859,
  18.55791504706697,
  44.71218593961408,
  76.15520357708716,
  180.40320086682667,
  57.061742459971526,
  27.6705878581027,
  68.33743480511626,
  18.620400431268088,
  24.530800057901743,
  45.41470309582157,
  58.47919510094424,
  25.125724069295355,
  60.330054943416016,
  24.995603050199712,
  17.815204255937715,
  27.004487503945562,
  72.27093023267976,
  27.9191763044152,
  49.5886643752733,
  90.17286374720061,
  91.15992024378313,
  103.96258110876433,
  5.282378293816456,
  43.65068488987147,
  80.06988595427131,
  12.08540678604747,
  63.60527539957998,
  13.278586783970654,
  14.36800305439691,
  18.69044990012204,
  9.038911259198743,
  5.043591136343069,
  60.89134810469868,
  126.12918335024405,
  47.29182710158029,
  11.256055453485402,
  10.08873642581787,
  17.177918251110125,
  11.551871204318578,
  79.49388357412488,
  46.83315449200037,
  91.81541689983051,
  57.182059266192965,
  38.18307005211219,
  8.573750701242474,
  68.7897489082834,
  49.55515559462687,
  54.6617441482015,
  13.684773019959055,
  8.597130172921673,
  56.43386582842407,
  121.76610873983014,
  93.3820218085413,
  114.77288350816902,
  33.57555467030754,
  68.70577305279268,
  126.18799902077353,
  11.806736893432724,
  67.5717274962382,
  31.09733312297488,
  24.9993178571767,
  50.91404483677396,
  74.99383888751169,
  17.164033408826647,
  53.625767602512845,
  5.077356712305322,
  17.048460788630074,
  24.490395523338314,
  32.51869845102124,
  7.40936319877423,
  14.008185744776586,
  58.83469407136514,
  15.96515050247702,
  140.4294180233351,
  15.00029204586558,
  58.959776751222755,
  29.00786598454628,
  22.385333494226465,
  68.08628641581663,
  88.3513608283919,
  57.52898471488527,
  9.170096934300823,
  59.202925898044924,
  10.383975078825948,
  15.908029817261188,
  68.44877639548497,
  13.964203147775955,
  38.12228485606426,
  14.657511901910114,
  55.24969493735518,
  28.418871179080792,
  13.634985284050115,
  20.593114874826995,
  7.268482365145665,
  13.79962092497581,
  8.765552540738394,
  7.052060997673832,
  83.56842636929363,
  44.714391618498134,
  38.14198033587367,
  14.593740943570765,
  14.15180293486147,
  126.14168179685748,
  47.96128619383243,
  26.788844991694212,
  34.126671958694054,
  138.51132064101247,
  33.89650422085187,
  139.78208201277732,
  17.92959290678197,
  14.058029814699449,
  11.700523940850971,
  63.92208444707154,
  68.67228810727593,
  67.17039380611489,
  82.80194814505042,
  65.55138104825598,
  159.2600351333081,
  15.186617889429188,
  19.564137007425693,
  16.56356673971569,
  7.4023145860666855,
  66.03785765602143,
  6.717644962731825,
  56.43136262214742,
  176.94409330818468,
  59.981433873743,
  89.21498829623785,
  21.48360170085487,
  6.803592364634948,
  15.452444273169021,
  67.33265802241313,
  75.33599995109903,
  29.369952890623452,
  53.35715545937019,
  33.1669456365453,
  88.0331667248147,
  82.28566960822343,
  72.11796308355532,
  22.64400155187669,
  19.100835933485165,
  16.23567369086813,
  11.608350448985593,
  7.232920380501984,
  48.77432665652317,
  18.870389674506736,
  26.612544181833155,
  76.42600867198205,
  14.759024019388667,
  31.133951694881247,
  23.202970867937008,
  8.68879198888946,
  112.10675524404846,
  5.249294819257113,
  73.16996781306247,
  10.088316178286297,
  12.626010192571085,
  71.06954982983893,
  96.4113934677713,
  25.716591665227995,
  74.70676345939859,
  15.061637140305567,
  31.349563703591844,
  14.765575628876116,
  41.37783741683598,
  6.607019620649585,
  20.016470084221307,
  21.39283167792491,
  54.415168874197455,
  89.48726452015974,
  8.974007576842856,
  10.37863948876639,
  10.010938030597671,
  160.44211687263984,
  5.425689625397781,
  5.988969731469766,
  14.868070754044982,
  16.703555716568903,
  5.032782901148494,
  25.747973838410505,
  12.732413670650539,
  175.34566492353284,
  23.424524974862887,
  32.0514130049495,
  91.65272460705333,
  39.36689932641691,
  14.38267800892593,
  26.6598596835178,
  33.63610288230363,
  98.20603715132502,
  65.26401332348384,
  313.33811760509525,
  61.039257054260936,
  89.77645657215288,
  19.995126261665384,
  38.79491955904897,
  26.78565655374219,
  8.227283651834314,
  15.077980884520434,
  50.93695406713049,
  56.412173704381,
  5.379946272545311,
  42.65603231855647,
  9.03805733257106,
  16.795191646795338,
  25.13183788510992,
  181.1554443409089,
  56.27861656223156,
  25.60811747633166,
  5.848399199904395,
  11.124725613811325,
  23.88846446546136,
  29.68447464326836,
  16.284857881132073,
  53.90884739056179,
  5.12946536475925,
  9.087724984729412,
  31.370536497249784,
  67.21793660523171,
  11.087284574123233,
  69.57536470147195,
  10.71050659714387,
  15.026406033500642,
  114.17159311802381,
  29.037351510597386,
  19.881878828036456,
  78.2986779092786,
  24.661532508463115,
  156.115212201316,
  60.461577464050656,
  46.51681115564065,
  64.58622582448633,
  9.989703034361797,
  53.15762941177178,
  85.27307199795871,
  176.305296323644,
  31.837548026754792,
  22.34985351742204,
  43.37568616079928,
  6.460963445932313,
  31.50834246245548,
  20.58552746852196,
  44.94705353616304,
  37.17461838773451,
  5.099194077332632,
  32.90993039607598,
  104.9916000738912,
  90.29404954860809,
  21.284921249870045,
  63.362470157772734,
  61.21107144910174,
  39.99791073761996,
  29.648512739081998,
  115.00463624065695,
  36.22922519753022,
  42.537829192352824,
  64.45626154113465,
  54.51225385741079,
  80.82654142283819,
  56.97474121039772,
  6.2683575364450475,
  6.3022993995514724,
  34.43252960671323,
  17.084033084814408,
  11.510704211004917,
  13.846648250539005,
  25.24238478363363,
  51.31691896485763,
  10.627951732475035,
  51.23120075689179,
  13.937854078563294,
  20.77498925745532,
  90.58011909951925,
  86.28657754369108,
  28.166031570200573,
  43.34354334706103,
  63.755564071280325,
  8.440272586614352,
  47.423728401763235,
  15.60628441250115,
  6.462352017126078,
  46.34345583843195,
  42.644513503265884,
  67.05489588593181,
  57.700090619293356,
  50.23355206617134,
  29.93019455834159,
  42.323728366451064,
  324.6762747630948,
  6.62776076943185,
  90.55236982052737,
  5.192636666306374,
  74.55620479146661,
  29.349448639369495,
  26.334210243547137,
  90.41760447507804,
  14.649971147491229,
  6.477199593589102,
  59.8287600317738,
  47.55215898458397,
  208.2547132477722,
  11.25065545428929,
  7.618536229633273,
  29.21729821751388,
  5.396782012222748,
  12.447092500610562,
  24.526174011643924,
  102.91069207739035,
  76.19809814420171,
  67.7368598307686,
  66.39835136692548,
  28.66129661975702,
  15.87962279485578,
  26.23462793271085,
  40.88067496199566,
  9.299442518713496,
  15.304711307160304,
  15.15215809710493,
  160.08728455182893,
  32.315523665824294,
  14.641212781246002,
  22.689494937996766,
  51.39737824257823,
  24.507053897445232,
  6.615755445004151,
  41.53936272775266,
  32.64512891996418,
  6.499348054041406,
  39.956518034349926,
  55.220384188716416,
  7.855200405854162,
  26.23293667230579,
  94.02840278444549,
  137.5732317944795,
  11.74106718209256,
  53.216812116956646,
  9.121092394572802,
  26.978770443785013,
  23.976403344690638,
  131.13900161001277,
  5.512709807354648,
  83.50587695776316,
  120.75754111154025,
  21.432627649464347,
  15.938218510478407,
  5.50792898869451,
  8.349560441435239,
  46.04330878224397,
  42.5215676405118,
  204.2566448446522,
  15.540972179187387,
  21.829801157422736,
  75.73574552508313,
  30.100486366045565,
  18.725531276530628,
  63.65354681442656,
  179.41647082102958,
  71.29779825656281,
  18.026845898146256,
  50.580100325133216,
  25.374202899652573,
  117.99260613097047,
  29.569000140727198,
  83.50072556539898,
  20.17129171681023,
  5.152129416571529,
  12.418697704452997,
  12.918247060528927,
  117.13813439275386,
  18.818056661350806,
  44.77227416705334,
  10.863722794633293,
  17.257761057982442,
  20.913078397867118,
  46.23452095016029,
  14.500311047434908,
  51.84250759072827,
  26.786096373728334,
  8.854665134358058,
  75.59907570488885,
  5.7794721015084685,
  16.311281090981943,
  44.1448549007064,
  11.90691444290629,
  66.32186209786246,
  8.88679758115977,
  5.3769170970019236,
  15.825348266967726,
  248.7647104338708,
  95.88248605413318,
  42.2519594346983,
  59.32768348766648,
  130.97140295118984,
  58.570285304433995,
  22.916988222888275,
  69.45381813591031,
  39.744556718555536,
  23.630047931521005,
  16.20378871730779,
  49.85046416041436,
  72.04361428474628,
  100.64556626023996,
  104.25108810520656,
  7.702191769034201,
  54.38835104227716,
  14.39549290924738,
  39.19447563840014,
  5.249218109261389,
  8.54203346361912,
  75.08793963255036,
  6.603855332714815,
  48.617068927368884,
  46.24467225265641,
  7.23011467546573,
  5.292418075307008,
  30.24723549342636,
  20.507620622784927,
  45.39737759814344,
  15.11899561294872,
  17.971645091223593,
  9.73349556867445,
  14.198727066882569,
  5.80524972819613,
  35.03792022298571,
  148.5710923952991,
  37.08055407342612,
  5.8737625917498555,
  29.335296710939083,
  45.148205970740364,
  48.20780377183974,
  13.897105916045245,
  105.22865824518311,
  62.41919324789416,
  71.19244883367816,
  18.461070889653755,
  47.67862493246133,
  20.242123879078814,
  10.880988443715635,
  31.141345074159133,
  5.637329293642736,
  23.265373497688167,
  9.93942213205608,
  70.24366398277022,
  70.64072638104335,
  11.627387288026377,
  29.40552434894623,
  11.79661355602266,
  25.17408269921721,
  8.028697841877118,
  95.86410563706423,
  54.15152188844783,
  7.6522072900293345,
  89.27910220284586,
  25.878933171395296,
  113.09481307963549,
  70.26634073319221,
  106.3503912943476,
  17.972447492576638,
  5.385817299348203,
  13.07458698314503,
  14.009872863302448,
  9.615292714087882,
  10.646363787476574,
  17.54145887959419,
  47.361354483070414,
  50.97301100206579,
  34.55933999788362,
  7.997463485011918,
  37.48030616404397,
  19.552018062377808,
  63.41710539045099,
  80.13770821137499,
  9.539568409049592,
  20.72629194435229,
  106.13315063164029,
  8.774496474086412,
  83.06297317549364,
  10.997196061981022,
  9.877489272470132,
  31.168933551461365,
  29.943901691072227,
  12.856009314396326,
  5.38110965223008,
  40.57188921262078,
  27.064820768877112,
  163.84810559956765,
  48.919726586645695,
  6.1713507826946055,
  16.65902539830648,
  148.78754953733747,
  12.420142059606889,
  16.92961780800558,
  38.43958774119711,
  81.94086042832038,
  19.83292741911609,
  53.37252801495728,
  48.740876131409756,
  75.47926639866053,
  139.69361338466956,
  17.578225571721887,
  11.340479901121961,
  6.194325530092807,
  21.090540718314756,
  46.88065943311076,
  58.030598924109796,
  185.99291227923194,
  6.604815087221636,
  11.34559824691836,
  10.77523594242167,
  134.63759007171868,
  22.303078670593184,
  11.462877499945172,
  111.20880094189025,
  49.33848601370531,
  107.42201470004261,
  37.451382911049514,
  29.716720293362197,
  89.94955215673986,
  6.072013284829785,
  28.963532023674382,
  126.91197967481754,
  60.25867905104224,
  20.691830293424896,
  78.29741749805024,
  54.37888426891168,
  16.937010273912406,
  7.352485570199575,
  97.59141527917566,
  54.79435777258159,
  108.27722485514072,
  92.26586728032673,
  11.579327749792661,
  12.431137668497986,
  13.577408393609055,
  40.9028747991367,
  41.16533159228178,
  112.71383667687874,
  11.31823441804757,
  72.9319093451652,
  16.638607330865103,
  10.440241915268214,
  155.98499119284438,
  26.8434550999985,
  32.54190314571381,
  10.519969282958632,
  118.51697667439582,
  29.154088198990237,
  40.72790912104638,
  53.04413414909574,
  60.907404841461506,
  146.38469734345168,
  8.642480698497081,
  71.43106764356533,
  28.402211411925478,
  22.065729394112545,
  68.10017307000676,
  77.00978165031549,
  5.309758889323243,
  27.99778135881106,
  18.661261356024767,
  54.69931075962305,
  33.81723448596374,
  82.41399349664184,
  35.47818686275496,
  19.08572738311426,
  32.41250777936446,
  34.35461860032439,
  24.415892744744156,
  5.13062813839644,
  13.612852055873017,
  52.009854991578194,
  5.580271875296687,
  26.442735264009773,
  5.746522042028261,
  20.328272984819076,
  21.965844840627856,
  11.929193471490276,
  20.79543070707279,
  22.89528470495999,
  12.183623574882127,
  24.3092891425887,
  74.3320937831429,
  28.127506546643282,
  19.483343445062122,
  182.37556629039452,
  32.25122617212112,
  15.600633867317907,
  57.24508359765869,
  8.048476474559743,
  20.242070402730736,
  56.63173754070698,
  5.036582668575584,
  5.152819356371734,
  18.94931258753779,
  10.112816446706383,
  65.31845825774111,
  8.885312023415915,
  8.299693222890182,
  30.474495787567868,
  5.634890417683859,
  15.578056641766537,
  26.38535609145219,
  17.812532591067523,
  15.266133509541646,
  6.314812793814502,
  20.4580238254229,
  33.10733152014791,
  34.98019841220422,
  6.223366604155347,
  26.163352467987494,
  18.177070756327915,
  5.047114330222604,
  17.10749424453423,
  13.416464862116142,
  15.544373022705157,
  63.20163433849206,
  16.11500586051435,
  38.872233052806855,
  7.277746090901515,
  8.757433423977686,
  40.600452373231555,
  7.0815913205529295,
  5.925992373342192,
  24.62537968757753,
  57.854075101505906,
  20.175740873767914,
  28.87070208056962,
  64.47387747461397,
  88.24292921784107,
  14.097615983920917,
  33.13578221402966,
  20.563481770435782,
  27.18087927802395,
  5.606927873887151,
  16.412155775802216,
  6.451790461334643,
  21.935676144929356,
  12.14082664236699,
  22.36388594216339,
  8.340142717409163,
  6.678983011082577,
  117.98343041161252,
  6.693610119808251,
  8.712969274987193,
  41.39149943080047,
  5.562127579715948,
  56.637391453645485,
  88.54937524550088,
  14.837974107430604,
  15.53116409073412,
  16.928024137461136,
  27.77850471320413,
  19.149004144307735,
  22.151958147466146,
  152.45963455967888,
  14.129868939481765,
  7.965612117307012,
  47.785447491745174,
  42.94388536101974,
  5.861437667978293,
  86.69090644946715,
  46.266655787238484,
  14.680345605269244,
  50.56865604680539,
  66.77499383058839,
  26.720263285608375,
  8.972851083138643,
  38.63317139744809,
  21.170853121561734,
  18.170342667268365,
  16.85524444517276,
  76.31174879923557,
  13.18518358807878,
  6.3255234517095325,
  28.336549432247992,
  39.86597397288527,
  19.315255818557585,
  14.917933649184565,
  14.661277350816956,
  5.198887787396429,
  12.353963742261262,
  52.24058896466551,
  13.706679082816027,
  17.293808976705286,
  58.1097379060424,
  12.107441245541743,
  10.174894450405562,
  41.79909314140049,
  33.18711371552374,
  6.532810932434207,
  5.783015978335344,
  49.658516952763996,
  5.1686063468196695,
  71.031184499903,
  18.244887258021098,
  70.92163024034326,
  34.84218192673624,
  120.18083984850456,
  15.594978685148265,
  9.485102212893365,
  15.134936863091614,
  10.859802037020275,
  8.687640683123133,
  84.05623191754816,
  26.995388524964685,
  23.41861402244887,
  55.981966176233556,
  27.94251376559162,
  16.093502188542093,
  17.552646316599713,
  19.87104944265511,
  50.110761261535,
  99.36883421197973,
  6.525270257340417,
  6.794619366257612,
  16.753376032615233,
  16.055384113524283,
  56.36045723468583,
  63.768470484820625,
  62.894571415145165,
  16.63593416169257,
  26.267111242820263,
  14.683615690575087,
  53.32980836243981,
  15.647127897861955,
  11.060981384999367,
  6.573082482109846])

Graph for home page

Detected energies during the discharge + Energy spectrum

multiplot(icon_fig, T_int_first,E,xle,Etot)