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

([3095096.875,
  3130134.375,
  3173053.125,
  3184339.0625,
  3236659.375,
  3339829.6875,
  3361368.75,
  3394290.625,
  3478837.5,
  3623385.9375,
  3630409.375,
  3638814.0625,
  3644310.9375,
  3749187.5,
  3800285.9375,
  3825940.625,
  3833795.3125,
  3835598.4375,
  3839553.125,
  3881589.0625,
  3888871.875,
  3936321.875,
  3976903.125,
  4060014.0625,
  4067768.75,
  4092753.125,
  4106103.125,
  4359567.1875,
  4518362.5,
  5543140.625,
  5544182.8125,
  5604010.9375,
  5605882.8125,
  5606962.5,
  5812875.0,
  5815982.8125,
  5816301.5625,
  5817520.3125,
  5822973.4375,
  5824389.0625,
  5827267.1875,
  5828703.125,
  5879342.1875,
  5885393.75,
  5895771.875,
  5913675.0,
  5924465.625,
  5924928.125,
  5952151.5625,
  5955779.6875,
  5969759.375,
  5974773.4375,
  5975303.125,
  5987084.375,
  6001139.0625,
  6002212.5,
  6003457.8125,
  6004635.9375,
  6005606.25,
  6010192.1875,
  6010570.3125,
  6038145.3125,
  6044790.625,
  6054506.25,
  6057182.8125,
  6060404.6875,
  6061629.6875,
  6061817.1875,
  6061825.0,
  6061843.75,
  6062412.5,
  6063631.25,
  6067809.375,
  6069576.5625,
  6139709.375,
  6140701.5625,
  6146695.3125,
  6148635.9375,
  6154395.3125,
  6226862.5,
  6230390.625,
  6234929.6875,
  6271028.125,
  6281798.4375,
  6281946.875,
  6312445.3125,
  6349659.375,
  6354142.1875,
  6354982.8125,
  6360348.4375,
  6366464.0625,
  6429153.125,
  6430348.4375,
  6430471.875,
  6430485.9375,
  6430554.6875,
  6476718.75,
  6481134.375,
  6491523.4375,
  6543426.5625,
  6549515.625,
  6562106.25,
  6574489.0625,
  6574529.6875,
  6585251.5625,
  6594785.9375,
  6597854.6875,
  6652389.0625,
  6660606.25,
  6667795.3125,
  6686635.9375,
  6709715.625,
  6909517.1875,
  6920348.4375,
  6920362.5,
  7085657.8125,
  7090815.625,
  7262906.25,
  7471714.0625,
  7512268.75,
  7623790.625,
  7670892.1875,
  7702435.9375,
  7746726.5625,
  7747443.75,
  7747725.0,
  7821670.3125,
  8293815.625,
  8407873.4375,
  8549767.1875,
  8581109.375,
  8673853.125,
  8674500.0,
  8844775.0,
  8883093.75,
  8958075.0,
  9048804.6875,
  9134260.9375,
  9166535.9375,
  9167731.25,
  9185896.875,
  9252093.75,
  9264071.875,
  9267556.25,
  9276921.875,
  9277060.9375,
  9287301.5625,
  9293276.5625,
  9295248.4375,
  9298473.4375,
  9298865.625,
  9299448.4375,
  9301309.375,
  9301535.9375,
  9302603.125,
  9302654.6875,
  9307309.375,
  9308431.25,
  9308510.9375,
  9322225.0,
  9325959.375,
  9327304.6875,
  9327531.25,
  9328395.3125,
  9330585.9375,
  9331032.8125,
  9331407.8125,
  9332526.5625,
  9332551.5625,
  9338215.625,
  9338270.3125,
  9339307.8125,
  9342260.9375,
  9344514.0625,
  9346723.4375,
  9347164.0625,
  9348017.1875,
  9351326.5625,
  9352309.375,
  9352423.4375,
  9358867.1875,
  9380614.0625,
  9384768.75,
  9385475.0,
  9574679.6875,
  9588153.125,
  9589546.875,
  9596085.9375,
  9597867.1875,
  9598028.125,
  9610746.875,
  9645923.4375,
  9743623.4375,
  9774156.25,
  9807400.0,
  9915476.5625,
  10069742.1875,
  10158950.0,
  10280985.9375,
  10289282.8125,
  10295045.3125,
  10295084.375,
  10301118.75,
  10307501.5625,
  10308054.6875,
  10308654.6875,
  10309025.0,
  10309135.9375,
  10309784.375,
  10309954.6875,
  10310046.875,
  10310456.25,
  10310623.4375,
  10311826.5625,
  10312025.0,
  10312292.1875,
  10312490.625,
  10312698.4375,
  10314018.75,
  10314235.9375,
  10314509.375,
  10314610.9375,
  10314942.1875,
  10317314.0625,
  10317371.875,
  10317496.875,
  10317634.375,
  10317743.75,
  10317978.125,
  10318615.625,
  10319348.4375,
  10319507.8125,
  10320118.75,
  10320273.4375,
  10320300.0,
  10320495.3125,
  10320731.25,
  10320787.5,
  10321564.0625,
  10322139.0625,
  10322340.625,
  10323145.3125,
  10323759.375,
  10324007.8125,
  10324859.375,
  10324929.6875,
  10325692.1875,
  10325932.8125,
  10326228.125,
  10326500.0,
  10327004.6875,
  10327128.125,
  10327271.875,
  10327967.1875,
  10328482.8125,
  10328885.9375,
  10328943.75,
  10329042.1875,
  10329334.375,
  10329971.875,
  10329978.125,
  10331275.0,
  10331809.375,
  10332126.5625,
  10332468.75,
  10332837.5,
  10333168.75,
  10333214.0625,
  10333496.875,
  10334393.75,
  10334576.5625,
  10335740.625,
  10335789.0625,
  10335893.75,
  10336884.375,
  10337393.75,
  10337604.6875,
  10338650.0,
  10339279.6875,
  10339659.375,
  10339803.125,
  10339995.3125,
  10340067.1875,
  10340114.0625,
  10340175.0,
  10340826.5625,
  10341267.1875,
  10342056.25,
  10343078.125,
  10343800.0,
  10343871.875,
  10344368.75,
  10345551.5625,
  10345629.6875,
  10345679.6875,
  10345717.1875,
  10346439.0625,
  10346998.4375,
  10350035.9375,
  10350800.0,
  10351921.875,
  10352040.625,
  10352481.25,
  10352781.25,
  10353268.75,
  10355371.875,
  10355970.3125,
  10356046.875,
  10358014.0625,
  10359557.8125,
  10360114.0625,
  10360432.8125,
  10360439.0625,
  10361782.8125,
  10361954.6875,
  10362314.0625,
  10362467.1875,
  10369871.875,
  10373509.375,
  10373914.0625,
  10388879.6875,
  10477029.6875,
  10498712.5,
  10501132.8125,
  10518271.875,
  10518284.375,
  10524175.0,
  10525876.5625,
  10528479.6875,
  10537621.875,
  10542632.8125,
  10546545.3125,
  10550665.625,
  10551973.4375,
  10568545.3125,
  10570568.75,
  10571379.6875,
  10575551.5625,
  10576414.0625,
  10583564.0625,
  10583575.0,
  10584500.0,
  10584512.5,
  10586576.5625,
  10586842.1875,
  10589575.0,
  10591045.3125,
  10593678.125,
  10593968.75,
  10599493.75,
  10604378.125,
  10606801.5625,
  10610676.5625,
  10620978.125,
  10621756.25,
  10625560.9375,
  10632734.375,
  10634217.1875,
  10637003.125,
  10641867.1875,
  10647748.4375,
  10650081.25,
  10650137.5,
  10651003.125,
  10651812.5,
  10652564.0625,
  10652648.4375,
  10652709.375,
  10653234.375,
  10653278.125,
  10660632.8125,
  10661256.25,
  10664578.125,
  10665196.875,
  10665653.125,
  10665859.375,
  10666350.0,
  10669685.9375,
  10672007.8125,
  10672059.375,
  10672292.1875,
  10674923.4375,
  10675757.8125,
  10677909.375,
  10677921.875,
  10677967.1875,
  10678023.4375,
  10680100.0,
  10683654.6875,
  10683773.4375,
  10687082.8125,
  10688095.3125,
  10689064.0625,
  10689368.75,
  10689915.625,
  10689975.0,
  10691718.75,
  10692067.1875,
  10694081.25,
  10694150.0,
  10694625.0,
  10695029.6875,
  10695581.25,
  10696026.5625,
  10697865.625,
  10698284.375,
  10698354.6875,
  10698550.0,
  10702175.0,
  10702253.125,
  10702845.3125,
  10703668.75,
  10705820.3125,
  10707839.0625,
  10709551.5625,
  10709714.0625,
  10709856.25,
  10711900.0,
  10713340.625,
  10714584.375,
  10715073.4375,
  10718389.0625,
  10722868.75,
  10801353.125,
  10801359.375,
  10805301.5625,
  10806518.75,
  10806573.4375,
  10806756.25,
  10807281.25,
  10810979.6875,
  10812654.6875,
  10821601.5625,
  10832964.0625,
  10836409.375,
  10838545.3125,
  10839285.9375,
  10839290.625,
  10840954.6875,
  10841642.1875,
  10841798.4375,
  10842746.875,
  10843996.875,
  10845215.625,
  10845837.5,
  10847076.5625,
  10847596.875,
  10847832.8125,
  10848335.9375,
  10848575.0,
  10849181.25,
  10849206.25,
  10849225.0,
  10850168.75,
  10850337.5,
  10851575.0,
  10851654.6875,
  10852021.875,
  10852251.5625,
  10852435.9375,
  10852439.0625,
  10852857.8125,
  10852867.1875,
  10853134.375,
  10853429.6875,
  10853762.5,
  10854515.625,
  10854521.875,
  10854643.75,
  10855050.0,
  10855078.125,
  10855346.875,
  10855401.5625,
  10859253.125,
  10859296.875,
  10859881.25,
  10860401.5625,
  10861117.1875,
  10861782.8125,
  10864178.125,
  10864331.25,
  10864904.6875,
  10865728.125,
  10865762.5,
  10866118.75,
  10866176.5625,
  10867034.375,
  10867429.6875,
  10867537.5,
  10867551.5625,
  10869642.1875,
  10869682.8125,
  10870432.8125,
  10871084.375,
  10875185.9375,
  10878128.125,
  10879346.875,
  10879662.5,
  10881098.4375,
  10881107.8125,
  10881153.125,
  10881179.6875,
  10881181.25,
  10881235.9375,
  10881264.0625,
  10881396.875,
  10881445.3125,
  10881451.5625,
  10881667.1875,
  10881771.875,
  10881839.0625,
  10881996.875,
  10882318.75,
  10882718.75,
  10882768.75,
  10882801.5625,
  10882878.125,
  10882926.5625,
  10882973.4375,
  10883081.25,
  10883118.75,
  10883125.0,
  10883190.625,
  10883379.6875,
  10883543.75,
  10883570.3125,
  10883667.1875,
  10883842.1875,
  10883878.125,
  10883901.5625,
  10883932.8125,
  10884385.9375,
  10884482.8125,
  10884726.5625,
  10885081.25,
  10885246.875,
  10885251.5625,
  10885273.4375,
  10885353.125,
  10885370.3125,
  10885437.5,
  10885493.75,
  10885643.75,
  10885656.25,
  10885673.4375,
  10885810.9375,
  10885904.6875,
  10886064.0625,
  10886087.5,
  10886093.75,
  10886135.9375,
  10886206.25,
  10886273.4375,
  10886292.1875,
  10886309.375,
  10886403.125,
  10886456.25,
  10886573.4375,
  10886714.0625,
  10886735.9375,
  10886764.0625,
  10886889.0625,
  10886921.875,
  10886973.4375,
  10886995.3125,
  10887090.625,
  10887093.75,
  10887131.25,
  10887260.9375,
  10887359.375,
  10887360.9375,
  10887437.5,
  10887570.3125,
  10887584.375,
  10887681.25,
  10887714.0625,
  10887718.75,
  10887746.875,
  10887750.0,
  10887770.3125,
  10887784.375,
  10887796.875,
  10887826.5625,
  10887960.9375,
  10888059.375,
  10888104.6875,
  10888156.25,
  10888250.0,
  10888350.0,
  10888368.75,
  10888415.625,
  10888465.625,
  10888537.5,
  10888548.4375,
  10888598.4375,
  10888628.125,
  10888739.0625,
  10888742.1875,
  10888746.875,
  10888860.9375,
  10888881.25,
  10888940.625,
  10888989.0625,
  10889076.5625,
  10889190.625,
  10889231.25,
  10889318.75,
  10889417.1875,
  10889478.125,
  10889490.625,
  10889529.6875,
  10889539.0625,
  10889576.5625,
  10889642.1875,
  10889710.9375,
  10889729.6875,
  10889800.0,
  10889804.6875,
  10889917.1875,
  10890192.1875,
  10890256.25,
  10890264.0625,
  10890296.875,
  10890343.75,
  10890446.875,
  10890506.25,
  10890820.3125,
  10891071.875,
  10891073.4375,
  10891145.3125,
  10891242.1875,
  10891323.4375,
  10891504.6875,
  10891531.25,
  10891567.1875,
  10891675.0,
  10891884.375,
  10892540.625,
  10893507.8125,
  10894684.375,
  10894834.375,
  10895203.125,
  10895268.75,
  10897412.5,
  10899225.0,
  10899317.1875,
  10902109.375,
  10906518.75,
  10911089.0625,
  10912175.0,
  10913253.125,
  10914320.3125,
  10914378.125,
  10914632.8125,
  10914775.0,
  10915095.3125,
  10915164.0625,
  10915228.125,
  10916170.3125,
  10916323.4375,
  10916429.6875,
  10916771.875,
  10917215.625,
  10917495.3125,
  10917556.25,
  10917664.0625,
  10918315.625,
  10918470.3125,
  10918898.4375,
  10919525.0,
  10919698.4375,
  10919989.0625,
  10920093.75,
  10920318.75,
  10920575.0,
  10920596.875,
  10920664.0625,
  10921070.3125,
  10921278.125,
  10921353.125,
  10921356.25,
  10921434.375,
  10921645.3125,
  10922104.6875,
  10922109.375,
  10922129.6875,
  10922198.4375,
  10922257.8125,
  10922264.0625,
  10922293.75,
  10922365.625,
  10922506.25,
  10922531.25,
  10922537.5,
  10922801.5625,
  10922806.25,
  10922835.9375,
  10922892.1875,
  10923278.125,
  10923501.5625,
  10923575.0,
  10923868.75,
  10923948.4375,
  10924093.75,
  10924262.5,
  10924364.0625,
  10924453.125,
  10924465.625,
  10925006.25,
  10925051.5625,
  10925153.125,
  10925348.4375,
  10925429.6875,
  10925596.875,
  10925925.0,
  10926098.4375,
  10926128.125,
  10926151.5625,
  10926454.6875,
  10926654.6875,
  10926770.3125,
  10926889.0625,
  10926915.625,
  10927321.875,
  10927596.875,
  10927726.5625,
  10927981.25,
  10928145.3125,
  10928543.75,
  10928789.0625,
  10929031.25,
  10929031.25,
  10929532.8125,
  10929587.5,
  10930218.75,
  10930218.75,
  10930232.8125,
  10930420.3125,
  10930454.6875,
  10930664.0625,
  10930860.9375,
  10930870.3125,
  10930900.0,
  10930928.125,
  10931104.6875,
  10931218.75,
  10931448.4375,
  10931451.5625,
  10931681.25,
  10931843.75,
  10932384.375,
  10932676.5625,
  10932721.875,
  10932757.8125,
  10932917.1875,
  10934662.5,
  10934687.5,
  10935039.0625,
  10935148.4375,
  10935503.125,
  10936704.6875,
  10937028.125,
  10937239.0625,
  10937546.875,
  10937560.9375,
  10937981.25,
  10938842.1875,
  10939871.875,
  10940214.0625,
  10940606.25,
  10940781.25,
  10940810.9375,
  10940960.9375,
  10941887.5,
  10942423.4375,
  10943606.25,
  10943848.4375,
  10944360.9375,
  10945128.125,
  10945957.8125,
  10946528.125,
  10946817.1875,
  10947159.375,
  10947937.5,
  10949028.125,
  10949226.5625,
  10949809.375,
  10950296.875,
  10950470.3125,
  10950673.4375,
  10951807.8125,
  10952471.875,
  10952725.0,
  10952778.125,
  10952784.375,
  10952873.4375,
  10952900.0,
  10953121.875,
  10954700.0,
  10955804.6875,
  10955865.625,
  10956117.1875,
  10956343.75,
  10956437.5,
  10958182.8125,
  10960417.1875,
  10961435.9375,
  10961564.0625,
  10961642.1875,
  10961832.8125,
  10961914.0625,
  10961915.625,
  10963120.3125,
  10963295.3125,
  10963384.375,
  10963650.0,
  10963768.75,
  10963793.75,
  10964090.625,
  10965314.0625,
  10965384.375,
  10965462.5,
  10966151.5625,
  10967268.75,
  10967918.75,
  10969079.6875,
  10969656.25,
  10969787.5,
  10971337.5,
  10972368.75,
  10972854.6875,
  10973035.9375,
  10973556.25,
  10973951.5625,
  10974559.375,
  10974835.9375,
  10975565.625,
  10975734.375,
  10977785.9375,
  10978518.75,
  10978632.8125,
  10978659.375,
  10978737.5,
  10979903.125,
  10980192.1875,
  10980512.5,
  10981354.6875,
  10982539.0625,
  10982834.375,
  10982890.625,
  10983115.625,
  10983232.8125,
  10983560.9375,
  10983590.625,
  10984081.25,
  10984831.25,
  10984979.6875,
  10985089.0625,
  10985362.5,
  10985415.625,
  10985467.1875,
  10985660.9375,
  10985714.0625,
  10985831.25,
  10986034.375,
  10986057.8125,
  10986114.0625,
  10986234.375,
  10986971.875,
  10987185.9375,
  10987251.5625,
  10987528.125,
  10987803.125,
  10987950.0,
  10988021.875,
  10988181.25,
  10988284.375,
  10988796.875,
  10988881.25,
  10989420.3125,
  10990409.375,
  10990484.375,
  10990981.25,
  10990982.8125,
  10991706.25,
  10992010.9375,
  10992103.125,
  10992479.6875,
  10992546.875,
  10992654.6875,
  10992743.75,
  10992884.375,
  10993004.6875,
  10993120.3125,
  10993168.75,
  10993217.1875,
  10993271.875,
  10993303.125,
  10993303.125,
  10993315.625,
  10993532.8125,
  10994193.75,
  10994250.0,
  10994668.75,
  10994884.375,
  10995017.1875,
  10995057.8125,
  10995295.3125,
  10995485.9375,
  10995782.8125,
  10995939.0625,
  10997456.25,
  10998806.25,
  11000315.625,
  11001193.75,
  11001395.3125,
  11001623.4375,
  11002075.0,
  11002526.5625,
  11002592.1875,
  11002617.1875,
  11002956.25,
  11003267.1875,
  11006328.125,
  11006500.0,
  11006693.75,
  11007612.5,
  11008003.125,
  11008017.1875,
  11008315.625,
  11008356.25,
  11009971.875],
 [22.021009242355184,
  38.4890783953468,
  47.02068742803747,
  28.392128982770146,
  66.74958698675533,
  55.22233102952246,
  5.767541976575615,
  85.98545176225845,
  76.9788880898114,
  25.313246423077054,
  26.001943124255078,
  92.83992457127137,
  22.140125767721898,
  6.775345668204583,
  23.317082582744632,
  29.846308908107645,
  60.91451679256374,
  7.328955250217951,
  10.965439446176706,
  54.96693523938706,
  5.0788987163568,
  49.9400207927214,
  8.494233821597778,
  17.79456726236055,
  26.614098956838923,
  24.077513630222768,
  7.881532955875548,
  51.47549429520111,
  38.09964661220158,
  108.72243320993135,
  48.32018342004417,
  12.352624450128737,
  15.216745508036997,
  47.56984033199902,
  22.11036563746642,
  14.560200383213502,
  13.21185426504061,
  17.00803663695596,
  77.58990132771481,
  18.928524810622605,
  18.26483942357279,
  5.528073577416909,
  27.787439662382923,
  15.147994241430474,
  77.08979221999557,
  44.22948272239949,
  5.82832345829075,
  41.997380502800254,
  121.56549943227482,
  89.5261471073743,
  16.279487112996037,
  18.588269635683083,
  31.178240603922152,
  133.54967171658197,
  46.77596037010173,
  26.653184434819728,
  79.97654101549485,
  65.43011049508247,
  17.955891771497473,
  61.90826001779653,
  20.868631608281472,
  30.946566757831633,
  61.02832664774631,
  165.8110614041866,
  53.76319944199383,
  20.34710226503025,
  26.1027783369538,
  18.50971128519458,
  12.056668303282246,
  17.21388380696328,
  89.09662348260042,
  11.821077989476422,
  200.72604604634267,
  5.4336971123951825,
  16.10064007640934,
  15.351378248519651,
  66.85188253566955,
  94.70977099622877,
  38.94508089319679,
  7.364347189617007,
  56.823631086798514,
  6.9515240150768,
  23.655785158591573,
  109.94492113787433,
  22.57318643421639,
  28.80568689671915,
  36.86961944261854,
  7.283532726832849,
  25.09127887004715,
  46.893205239032824,
  29.074893782244324,
  137.60778568136072,
  45.539803689582556,
  34.73445517965499,
  98.36953253093536,
  5.276151751162502,
  16.73576318457908,
  20.490710855333738,
  11.666800405834781,
  40.666781069542665,
  95.499035707235,
  60.32844883773057,
  53.09497305128083,
  11.138284956021902,
  34.41457539434679,
  18.538386078888358,
  47.81573168259413,
  13.823627868534967,
  88.07231232320218,
  16.841576235918975,
  6.313443538055007,
  46.01188192946974,
  35.4555126037289,
  62.31103758371546,
  18.260874149808338,
  70.50340101649549,
  16.998355474105303,
  32.349993325646366,
  36.35025351582656,
  9.905739734247941,
  6.272346554264856,
  16.548075951133608,
  173.4780343651496,
  14.66564857523588,
  18.176815182647395,
  25.105044093816872,
  16.020130987277902,
  19.486453323053702,
  20.430833529285074,
  27.70935736374711,
  36.63704332839521,
  27.408499881318978,
  19.671181140231766,
  11.382255464408123,
  16.694242426628662,
  72.4963068669911,
  16.224052722840963,
  32.937620370474534,
  13.054374663405506,
  53.593902462809424,
  37.660485575055304,
  124.6842273463259,
  11.973016827946285,
  20.17553459537155,
  24.393674603768687,
  5.041678731527717,
  23.24902946044456,
  76.2990777405718,
  18.977857621516563,
  62.036722606662316,
  7.0193588417710515,
  10.981624497878927,
  8.63548794277709,
  19.284447395602363,
  91.59509239318692,
  5.239379273968614,
  18.311129497436117,
  83.30861089082661,
  68.88100814463934,
  11.500091722865413,
  81.27556251771465,
  38.88334790456944,
  17.95606701292393,
  34.5942353457635,
  11.094001898801896,
  30.218668317589874,
  28.793985019134244,
  35.00334716868556,
  89.27003577678465,
  62.178988330270435,
  5.786125626670979,
  15.4954564536649,
  14.245958827256203,
  12.20793513230758,
  17.964454316839333,
  98.4255103702047,
  17.351265301932973,
  7.86100509754009,
  22.050736113141603,
  26.438877054565015,
  88.08701175903134,
  13.770208485666899,
  88.72686862228527,
  163.30563692716265,
  52.090468440495556,
  21.586448490358404,
  9.474736021000595,
  36.213764981720736,
  85.08699727616452,
  214.53261368846267,
  10.740221575344053,
  35.77077205935591,
  91.18864709509413,
  16.229036704589458,
  10.799802374587344,
  94.71323087772357,
  39.0077658431792,
  5.621626020592023,
  21.458026174122292,
  39.61612410696293,
  23.427202472568087,
  30.099284242926736,
  157.64861930971978,
  30.809233968790622,
  5.8023261510349915,
  13.519787243402202,
  13.73244761866103,
  14.153340818472302,
  17.512550123563887,
  5.799873310174771,
  65.81773487747353,
  31.09619083558779,
  24.440747308417343,
  144.62507467580127,
  65.51660079114662,
  33.451659394009475,
  8.00529218687174,
  59.24554381775309,
  8.748882696933995,
  55.41716628969277,
  5.831643974250857,
  223.5296764355819,
  30.710301430118875,
  8.697083915282567,
  30.55764708558197,
  29.11562981964126,
  7.71821965787399,
  5.332107034424148,
  8.953723696357054,
  12.740229987697337,
  14.068827868213816,
  36.27327213427412,
  32.020303071632846,
  44.60634230266353,
  16.57414431456052,
  6.434576193701606,
  14.907339095605545,
  22.890187732247284,
  5.24721410341351,
  25.26679854568272,
  15.271712538456171,
  36.6513410796253,
  5.383907988892966,
  24.901555215662682,
  25.036365624854696,
  139.8300860633072,
  9.680220528056239,
  6.11022503291402,
  15.076892778596527,
  31.1540587976988,
  22.547681036207095,
  11.572712877285262,
  65.87077213573423,
  37.741516255717784,
  98.86840405537988,
  8.814990455582633,
  41.37817426786875,
  5.204984731342885,
  41.52914443715207,
  106.0944401293016,
  15.021078390882824,
  22.82330875321483,
  16.66503134240758,
  64.0284232252087,
  43.765999967727694,
  56.41239548588396,
  31.73265188069879,
  31.92621215244271,
  39.5202020939635,
  12.743189561458843,
  19.349836969610667,
  17.93846462401773,
  49.29586602947651,
  15.279825888494805,
  34.50729847707833,
  77.14075514074692,
  37.47004128373334,
  10.928065331906891,
  22.9258506804753,
  151.8548261866377,
  50.36576040721054,
  57.160279098744134,
  5.1463006665050255,
  16.246580314279104,
  37.54345238010842,
  36.79520231005443,
  26.596550451061262,
  30.02696603035968,
  30.29493832764156,
  66.75335578822245,
  22.47519437477357,
  20.362495308891447,
  13.782695860583036,
  5.105228635646729,
  11.465829376306463,
  24.689489332707915,
  22.468568125350938,
  18.251894050617672,
  10.41531876261297,
  13.710940199974536,
  52.34867932798088,
  42.540561392271215,
  9.865211756656256,
  19.86159492684419,
  13.716269438775196,
  8.12563030131686,
  6.079465494844198,
  24.220883659054937,
  7.680922032547267,
  30.718974200800645,
  7.282300362708002,
  67.03591453283133,
  44.46697470882711,
  9.25112182635102,
  35.64006940966409,
  7.054792151097475,
  36.9233388624636,
  38.408505668803954,
  31.26842943641544,
  44.97547688189459,
  72.49704185320105,
  18.03348940304471,
  128.7066757451429,
  134.79254782541915,
  9.099341105235155,
  28.009804625469606,
  6.209783846193417,
  65.0208204057012,
  13.855440865426292,
  115.08518137467516,
  37.232730994110085,
  10.031061446947488,
  59.63945759145744,
  76.92252500768866,
  41.165820552689155,
  95.07188602774748,
  24.657813383502127,
  11.137720828219305,
  10.496877095393701,
  25.47161697023405,
  55.64807023954094,
  31.104682557426386,
  13.26363570340616,
  43.488304977597295,
  13.233473167898177,
  15.360079603334214,
  22.848670904501805,
  42.559909705301784,
  34.67602754924622,
  37.97093729504877,
  89.55418966788866,
  23.655749886430065,
  25.256919586640304,
  57.95666255489811,
  26.289583863198995,
  69.03233943083873,
  87.2797337057543,
  6.682150561842103,
  24.833240771015063,
  47.868893005645724,
  31.129884345645216,
  68.4942083232207,
  5.181325306025436,
  37.8727256560143,
  69.26176883543641,
  43.32495886566604,
  43.07275486455848,
  15.867927464700793,
  56.842124237070955,
  11.026938422475894,
  29.742461821466048,
  35.477831521480155,
  109.3976900512591,
  22.06295365220376,
  34.197541967758205,
  16.91099457022679,
  23.679096116763116,
  56.24427229374366,
  13.056530399021842,
  5.245650491351172,
  16.585874732501,
  9.519687281671542,
  5.3110270733573826,
  49.98210644744865,
  10.533285023869338,
  14.03520230293683,
  80.53185550481206,
  142.20907185489466,
  13.621867184182264,
  10.496795945218656,
  52.794487016302504,
  5.81852381302713,
  5.11988686230411,
  13.010602722255353,
  72.03845711836274,
  19.96349716340491,
  7.4510984617700595,
  6.519807299136634,
  76.72752828676161,
  14.955102447042108,
  17.488279106598295,
  13.2240439570794,
  77.17335669046227,
  33.601244916449915,
  43.65337426395233,
  10.535746741680837,
  21.54555258334053,
  20.873519761088723,
  28.238810736462195,
  5.778985931681963,
  7.711172389454429,
  52.85776382538101,
  81.54478801636546,
  9.975727577410039,
  5.156266483097807,
  88.13017486023823,
  17.535101307001124,
  37.6776960318424,
  62.6841803427774,
  14.224764176988328,
  18.411617199313568,
  41.170161100909084,
  15.974723954580107,
  21.634210021222763,
  35.126478197135846,
  29.38502145942212,
  12.113958991867475,
  91.67982697504785,
  38.62457881105571,
  31.782638781353505,
  64.6841965692721,
  11.417150361034581,
  73.87486338456806,
  8.806014972258357,
  18.95679500699656,
  39.14002171153467,
  22.41591445096192,
  65.99529546426692,
  28.67072992745748,
  8.805729919763305,
  65.99059132796356,
  28.11468680341484,
  76.44349714650832,
  66.16629340087987,
  10.152127676087854,
  6.357698827331542,
  84.66973836414083,
  10.592715304784974,
  65.64720746561785,
  47.66557987703953,
  23.33183279362798,
  44.579767139026444,
  16.380984123555233,
  7.905809331005398,
  65.33143082728952,
  38.16439296293021,
  61.836374435680966,
  37.69603815554154,
  22.733831852294873,
  10.911905571870372,
  74.65633083934625,
  12.508210428003846,
  9.725340240774221,
  48.09394253109295,
  97.45586146597684,
  84.18692986328307,
  7.898691548649248,
  10.68649086219551,
  37.46424444986359,
  24.73213961121516,
  41.872491913159266,
  40.62338597392003,
  17.05929120717752,
  12.748262745258055,
  9.773554885743303,
  28.227988516823736,
  5.433260964862621,
  17.68868122704911,
  42.4351797973745,
  48.472047921834985,
  9.856016171300812,
  9.06385286707796,
  18.72618452329767,
  17.19172638298216,
  77.62683033630735,
  5.350408738443339,
  81.43577629676666,
  57.4512005163295,
  65.90236095603541,
  28.558262024427805,
  36.25644384351509,
  114.36925168269248,
  16.345061629878007,
  78.91918042756956,
  5.3961460087826705,
  184.18845408512198,
  5.44601500182603,
  18.178213819118735,
  38.768675227088394,
  57.638076333475645,
  34.11240184579252,
  42.55537570015323,
  14.561482675649726,
  40.89795889076275,
  10.857965566203532,
  87.63891674092557,
  73.51387920871596,
  17.6120719174898,
  5.249086831333948,
  34.32455983600745,
  5.171124406431231,
  9.094198901350126,
  27.623466950636892,
  45.21400875983571,
  28.284419286661436,
  111.75056451168965,
  33.67624455001214,
  6.813122942536255,
  6.112532485850941,
  83.34096537565901,
  97.63111396509836,
  38.36411684364267,
  79.12613931639416,
  19.856801445326067,
  63.479904976588095,
  99.23966045508843,
  5.716553670554265,
  203.4042285228039,
  36.85084191552412,
  13.847658985820932,
  42.40245615647174,
  7.7619145515877115,
  10.125094502875928,
  22.59347082126714,
  35.5995571271542,
  49.83339886896031,
  13.023547661770873,
  86.32743717240116,
  5.365018795285513,
  30.15722298452897,
  53.3130026198923,
  40.11615393095157,
  104.90759889866919,
  22.038872957895677,
  111.06000894916396,
  86.21779869376928,
  16.648411062053647,
  15.122188221791173,
  119.67093414439576,
  56.10826696282055,
  30.4619104345695,
  85.57323752579141,
  26.64526035465237,
  5.512210331251039,
  16.168188756214402,
  50.16706543608371,
  7.862688686332413,
  21.174991404603016,
  29.770891245999692,
  53.83780875022511,
  45.27308263126477,
  18.849458107855693,
  22.682093444403883,
  9.267214533984296,
  84.65023467444472,
  19.710697184538045,
  22.2028086579542,
  29.509347136265312,
  6.933488576430554,
  22.683479401268528,
  5.271866949848501,
  48.45177151309455,
  94.34272152631485,
  50.74603832448135,
  65.8893793756094,
  59.18668021183142,
  6.411250187473843,
  34.330131028621736,
  40.6576972695079,
  11.104617550258228,
  7.887705892111198,
  5.18005252761478,
  19.746549438652668,
  92.05133031035159,
  35.05692637374621,
  45.4773152326746,
  23.73527356299942,
  22.16014330001128,
  7.680815981306755,
  28.249317472032615,
  94.17449764976114,
  12.395297198884663,
  65.55968764439035,
  8.188548098923544,
  19.76390527371317,
  13.782718952804494,
  27.463261022589784,
  7.253873108500993,
  7.098247328973293,
  87.34684372462829,
  66.49412408220236,
  7.283875763722794,
  17.330554785575572,
  36.52366381780125,
  8.01223805647984,
  56.429363548751994,
  54.06117003741237,
  14.770592630022241,
  26.42326911590387,
  59.868183019168306,
  43.23805501909495,
  10.214640248906216,
  59.72851590158683,
  20.659875269530673,
  23.493264175909378,
  31.204148809644746,
  25.198365358576496,
  10.903033004281209,
  51.03699798151449,
  26.935533597535574,
  52.816519902510024,
  28.4479667681646,
  13.573508083331795,
  13.31283556640876,
  59.11371970603487,
  27.41232303341799,
  63.92876644192435,
  16.3301677067414,
  30.17306679133136,
  73.05017821255817,
  43.76804684322735,
  9.377318534420954,
  101.34424601873934,
  41.035859801503136,
  24.265292619418787,
  8.277266186050756,
  34.45577091187768,
  63.05802953976904,
  11.163152759147,
  15.490200850043493,
  7.000659673439083,
  18.65522636836891,
  8.611027181783722,
  15.7071408820859,
  15.671686514707517,
  15.551000979005742,
  55.86499579745119,
  21.053437689644273,
  73.85920695239307,
  21.616976840544154,
  113.0408594774431,
  15.286231190939306,
  24.26820016859054,
  35.35622195433799,
  65.41433752029411,
  59.08459494152116,
  23.314415555225686,
  5.617862029888876,
  12.961658684330763,
  14.75647783470239,
  8.73183091086844,
  29.614177463402136,
  20.342343756181247,
  30.213886107224802,
  5.058864945783083,
  17.524693815422285,
  25.07481566810459,
  9.215964805219468,
  15.434852025089148,
  10.798638114703666,
  24.564514480703515,
  25.390587529801348,
  18.9742667465705,
  18.908219815034375,
  30.34449822091364,
  14.624239766767865,
  46.16688792442261,
  20.307999884402832,
  12.946293858889382,
  48.208410674804554,
  5.062074776884923,
  16.95465162004606,
  16.6640376812401,
  47.45443784378364,
  19.370594172227264,
  29.58239319234457,
  65.03935836284646,
  5.339394518650768,
  5.874130199960952,
  22.421037076783215,
  46.58539084372967,
  7.990716144309684,
  22.979159192473006,
  7.193342301813022,
  56.42877660935584,
  38.63725540578667,
  21.01436858523974,
  14.171886003128202,
  63.64548751945541,
  29.664315790507796,
  74.3046080298979,
  18.35924808820434,
  14.472804136288927,
  60.94876527352899,
  14.602847740925831,
  116.40782153342616,
  8.992212337913951,
  5.9875307768331725,
  78.45442709647688,
  74.62526212553288,
  46.40240576676584,
  12.054653876271365,
  5.285194301886191,
  7.839299582827493,
  63.392201058440435,
  54.77349197330699,
  12.847227483279545,
  71.80080956859267,
  14.586885534611513,
  11.755643620206742,
  16.436451660646064,
  10.46448021601351,
  16.851498424285793,
  8.626307341908115,
  35.91879378318802,
  46.550582331435805,
  8.046991122943412,
  55.17416982013076,
  46.85609323020705,
  123.32736537632918,
  10.134955948997298,
  29.108489257247527,
  21.27870430558702,
  25.34670583857635,
  49.5673741311398,
  71.76938397849558,
  23.111600725871835,
  50.19117509704823,
  8.680596800741602,
  55.06642405893906,
  10.329415999954492,
  90.42979675993865,
  29.71168000827199,
  23.82351004444102,
  35.71374341061582,
  5.510616214255254,
  55.83269388527026,
  5.203303906647725,
  7.079532911234448,
  10.34548884351672,
  23.688498973806844,
  44.93005624181225,
  47.62743080132242,
  19.80910117496984,
  15.33786883264623,
  61.1266382591012,
  43.013515943568244,
  124.00386625107237,
  52.01179770924548,
  81.10372725906613,
  15.249657599687554,
  12.384611721622107,
  17.71189451480888,
  34.728231331536676,
  20.475720516204397,
  65.64446515560198,
  9.566740644794663,
  11.445634567423177,
  68.9017362685597,
  29.98705052870775,
  11.139339005510818,
  13.30437144220927,
  27.846160543068166,
  5.674755432438352,
  12.343400811292362,
  16.94600758369818,
  20.42515646291394,
  29.364399613952926,
  61.024032258373,
  7.105237067849262,
  11.936150474254427,
  20.701798610012098,
  27.923475131570207,
  21.169829402338863,
  10.678265085728777,
  103.51254346744452,
  72.44795854086237,
  10.915504204385975,
  6.776016059219304,
  17.290112236251346,
  14.548129569537364,
  21.262300580920257,
  14.650993960765222,
  19.795831747331114,
  5.928722362066054,
  5.715891480488105,
  20.707832994601986,
  17.599363625216174,
  13.4017068908327,
  14.668583908856826,
  25.27435717097639,
  20.64032799362026,
  81.86736506133546,
  15.225643744904625,
  87.47437020994107,
  29.30794757172552,
  18.102969558140508,
  8.438859460490141,
  56.3849649656277,
  34.13813401389431,
  36.0818578019151,
  86.66110046385121,
  20.708465544551917,
  77.10642400050229,
  44.34881107527772,
  14.713167720706386,
  71.51845277837111,
  22.375394475295273,
  16.28818993360847,
  48.57040585597619,
  41.96772091864355,
  22.283426862356265,
  71.72591491862836,
  20.626442323309714,
  84.91838809504783,
  5.292777944683193,
  13.336394033125114,
  44.48343760206406,
  57.336228013175926,
  5.410824069280008,
  32.36725898468573,
  85.74413940363755,
  13.96525748953117,
  36.63488752594296,
  22.31510641626128,
  32.39794798722611,
  17.56963513315327,
  73.29718922332276,
  28.390550148637203,
  78.737325721106,
  40.13258517613334,
  24.2569644572232,
  5.73677101680262,
  17.007721219121464,
  27.546531131391724,
  17.88566874786028,
  29.04182172084941,
  5.098931286238759,
  5.362676915658964,
  60.05818336903028,
  21.74389677485656,
  37.48173029369633,
  7.21879373844304,
  44.515530418320196,
  142.56983849035865,
  18.138350794524957,
  7.616096618973378,
  14.289807409547906,
  35.20519866100101,
  27.285240800385132,
  14.368589102463558,
  51.62647468520487,
  25.991859649822743,
  9.074751637838316,
  9.046592571905157,
  14.584119126960067,
  5.937328771401298,
  10.045440512891627,
  23.313660511017435,
  17.107995963080118,
  5.818522107576128,
  71.08696593997973,
  25.442802520418088,
  41.45241381870106,
  66.22400949133461,
  146.57638690000948,
  105.88807152064071,
  5.960893475895645,
  6.148202775042389,
  21.192313707992646,
  44.09465620369771,
  25.457030630336206,
  12.425603059532982,
  24.6471178701932,
  35.382632947966556,
  17.186985693617892,
  5.161512798362923,
  11.884989819364487,
  5.0549004826156985,
  6.386668212430329,
  15.697027261721002,
  16.242766529518974,
  10.170805794759305,
  17.497511785217142,
  78.80800858733525,
  88.79806092413611,
  23.144164809664098,
  13.492445974154005,
  50.9021917385543,
  52.65022817104952,
  121.10695939461881,
  44.3338198699621,
  15.139243419728508,
  46.04415661762758,
  7.385940347398921,
  7.797147971087979,
  79.39219995992596,
  21.09274255073534,
  25.550954959178696,
  8.575364572269729,
  69.00936254497469,
  82.56555452336234,
  11.766986380727856,
  67.68651563337193,
  69.35916430244502,
  79.07667980291255,
  25.504299287895186,
  53.09170737629367,
  12.84771979132092,
  64.6783645165732,
  30.29273264034638,
  20.994389289618372,
  101.47276786301427,
  6.200436275231995,
  27.108089089319712,
  22.83534139160889,
  6.223195707648995,
  14.126174705444814,
  9.093728426326399,
  7.587944510358162,
  66.22406423386317,
  71.55801092103763])

Energy spectrum

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

([3095096.875,
  3130134.375,
  3173053.125,
  3184339.0625,
  3236659.375,
  3339829.6875,
  3361368.75,
  3394290.625,
  3478837.5,
  3623385.9375,
  3630409.375,
  3638814.0625,
  3644310.9375,
  3749187.5,
  3800285.9375,
  3825940.625,
  3833795.3125,
  3835598.4375,
  3839553.125,
  3881589.0625,
  3888871.875,
  3936321.875,
  3976903.125,
  4060014.0625,
  4067768.75,
  4092753.125,
  4106103.125,
  4359567.1875,
  4518362.5,
  5543140.625,
  5544182.8125,
  5604010.9375,
  5605882.8125,
  5606962.5,
  5812875.0,
  5815982.8125,
  5816301.5625,
  5817520.3125,
  5822973.4375,
  5824389.0625,
  5827267.1875,
  5828703.125,
  5879342.1875,
  5885393.75,
  5895771.875,
  5913675.0,
  5924465.625,
  5924928.125,
  5952151.5625,
  5955779.6875,
  5969759.375,
  5974773.4375,
  5975303.125,
  5987084.375,
  6001139.0625,
  6002212.5,
  6003457.8125,
  6004635.9375,
  6005606.25,
  6010192.1875,
  6010570.3125,
  6038145.3125,
  6044790.625,
  6054506.25,
  6057182.8125,
  6060404.6875,
  6061629.6875,
  6061817.1875,
  6061825.0,
  6061843.75,
  6062412.5,
  6063631.25,
  6067809.375,
  6069576.5625,
  6139709.375,
  6140701.5625,
  6146695.3125,
  6148635.9375,
  6154395.3125,
  6226862.5,
  6230390.625,
  6234929.6875,
  6271028.125,
  6281798.4375,
  6281946.875,
  6312445.3125,
  6349659.375,
  6354142.1875,
  6354982.8125,
  6360348.4375,
  6366464.0625,
  6429153.125,
  6430348.4375,
  6430471.875,
  6430485.9375,
  6430554.6875,
  6476718.75,
  6481134.375,
  6491523.4375,
  6543426.5625,
  6549515.625,
  6562106.25,
  6574489.0625,
  6574529.6875,
  6585251.5625,
  6594785.9375,
  6597854.6875,
  6652389.0625,
  6660606.25,
  6667795.3125,
  6686635.9375,
  6709715.625,
  6909517.1875,
  6920348.4375,
  6920362.5,
  7085657.8125,
  7090815.625,
  7262906.25,
  7471714.0625,
  7512268.75,
  7623790.625,
  7670892.1875,
  7702435.9375,
  7746726.5625,
  7747443.75,
  7747725.0,
  7821670.3125,
  8293815.625,
  8407873.4375,
  8549767.1875,
  8581109.375,
  8673853.125,
  8674500.0,
  8844775.0,
  8883093.75,
  8958075.0,
  9048804.6875,
  9134260.9375,
  9166535.9375,
  9167731.25,
  9185896.875,
  9252093.75,
  9264071.875,
  9267556.25,
  9276921.875,
  9277060.9375,
  9287301.5625,
  9293276.5625,
  9295248.4375,
  9298473.4375,
  9298865.625,
  9299448.4375,
  9301309.375,
  9301535.9375,
  9302603.125,
  9302654.6875,
  9307309.375,
  9308431.25,
  9308510.9375,
  9322225.0,
  9325959.375,
  9327304.6875,
  9327531.25,
  9328395.3125,
  9330585.9375,
  9331032.8125,
  9331407.8125,
  9332526.5625,
  9332551.5625,
  9338215.625,
  9338270.3125,
  9339307.8125,
  9342260.9375,
  9344514.0625,
  9346723.4375,
  9347164.0625,
  9348017.1875,
  9351326.5625,
  9352309.375,
  9352423.4375,
  9358867.1875,
  9380614.0625,
  9384768.75,
  9385475.0,
  9574679.6875,
  9588153.125,
  9589546.875,
  9596085.9375,
  9597867.1875,
  9598028.125,
  9610746.875,
  9645923.4375,
  9743623.4375,
  9774156.25,
  9807400.0,
  9915476.5625,
  10069742.1875,
  10158950.0,
  10280985.9375,
  10289282.8125,
  10295045.3125,
  10295084.375,
  10301118.75,
  10307501.5625,
  10308054.6875,
  10308654.6875,
  10309025.0,
  10309135.9375,
  10309784.375,
  10309954.6875,
  10310046.875,
  10310456.25,
  10310623.4375,
  10311826.5625,
  10312025.0,
  10312292.1875,
  10312490.625,
  10312698.4375,
  10314018.75,
  10314235.9375,
  10314509.375,
  10314610.9375,
  10314942.1875,
  10317314.0625,
  10317371.875,
  10317496.875,
  10317634.375,
  10317743.75,
  10317978.125,
  10318615.625,
  10319348.4375,
  10319507.8125,
  10320118.75,
  10320273.4375,
  10320300.0,
  10320495.3125,
  10320731.25,
  10320787.5,
  10321564.0625,
  10322139.0625,
  10322340.625,
  10323145.3125,
  10323759.375,
  10324007.8125,
  10324859.375,
  10324929.6875,
  10325692.1875,
  10325932.8125,
  10326228.125,
  10326500.0,
  10327004.6875,
  10327128.125,
  10327271.875,
  10327967.1875,
  10328482.8125,
  10328885.9375,
  10328943.75,
  10329042.1875,
  10329334.375,
  10329971.875,
  10329978.125,
  10331275.0,
  10331809.375,
  10332126.5625,
  10332468.75,
  10332837.5,
  10333168.75,
  10333214.0625,
  10333496.875,
  10334393.75,
  10334576.5625,
  10335740.625,
  10335789.0625,
  10335893.75,
  10336884.375,
  10337393.75,
  10337604.6875,
  10338650.0,
  10339279.6875,
  10339659.375,
  10339803.125,
  10339995.3125,
  10340067.1875,
  10340114.0625,
  10340175.0,
  10340826.5625,
  10341267.1875,
  10342056.25,
  10343078.125,
  10343800.0,
  10343871.875,
  10344368.75,
  10345551.5625,
  10345629.6875,
  10345679.6875,
  10345717.1875,
  10346439.0625,
  10346998.4375,
  10350035.9375,
  10350800.0,
  10351921.875,
  10352040.625,
  10352481.25,
  10352781.25,
  10353268.75,
  10355371.875,
  10355970.3125,
  10356046.875,
  10358014.0625,
  10359557.8125,
  10360114.0625,
  10360432.8125,
  10360439.0625,
  10361782.8125,
  10361954.6875,
  10362314.0625,
  10362467.1875,
  10369871.875,
  10373509.375,
  10373914.0625,
  10388879.6875,
  10477029.6875,
  10498712.5,
  10501132.8125,
  10518271.875,
  10518284.375,
  10524175.0,
  10525876.5625,
  10528479.6875,
  10537621.875,
  10542632.8125,
  10546545.3125,
  10550665.625,
  10551973.4375,
  10568545.3125,
  10570568.75,
  10571379.6875,
  10575551.5625,
  10576414.0625,
  10583564.0625,
  10583575.0,
  10584500.0,
  10584512.5,
  10586576.5625,
  10586842.1875,
  10589575.0,
  10591045.3125,
  10593678.125,
  10593968.75,
  10599493.75,
  10604378.125,
  10606801.5625,
  10610676.5625,
  10620978.125,
  10621756.25,
  10625560.9375,
  10632734.375,
  10634217.1875,
  10637003.125,
  10641867.1875,
  10647748.4375,
  10650081.25,
  10650137.5,
  10651003.125,
  10651812.5,
  10652564.0625,
  10652648.4375,
  10652709.375,
  10653234.375,
  10653278.125,
  10660632.8125,
  10661256.25,
  10664578.125,
  10665196.875,
  10665653.125,
  10665859.375,
  10666350.0,
  10669685.9375,
  10672007.8125,
  10672059.375,
  10672292.1875,
  10674923.4375,
  10675757.8125,
  10677909.375,
  10677921.875,
  10677967.1875,
  10678023.4375,
  10680100.0,
  10683654.6875,
  10683773.4375,
  10687082.8125,
  10688095.3125,
  10689064.0625,
  10689368.75,
  10689915.625,
  10689975.0,
  10691718.75,
  10692067.1875,
  10694081.25,
  10694150.0,
  10694625.0,
  10695029.6875,
  10695581.25,
  10696026.5625,
  10697865.625,
  10698284.375,
  10698354.6875,
  10698550.0,
  10702175.0,
  10702253.125,
  10702845.3125,
  10703668.75,
  10705820.3125,
  10707839.0625,
  10709551.5625,
  10709714.0625,
  10709856.25,
  10711900.0,
  10713340.625,
  10714584.375,
  10715073.4375,
  10718389.0625,
  10722868.75,
  10801353.125,
  10801359.375,
  10805301.5625,
  10806518.75,
  10806573.4375,
  10806756.25,
  10807281.25,
  10810979.6875,
  10812654.6875,
  10821601.5625,
  10832964.0625,
  10836409.375,
  10838545.3125,
  10839285.9375,
  10839290.625,
  10840954.6875,
  10841642.1875,
  10841798.4375,
  10842746.875,
  10843996.875,
  10845215.625,
  10845837.5,
  10847076.5625,
  10847596.875,
  10847832.8125,
  10848335.9375,
  10848575.0,
  10849181.25,
  10849206.25,
  10849225.0,
  10850168.75,
  10850337.5,
  10851575.0,
  10851654.6875,
  10852021.875,
  10852251.5625,
  10852435.9375,
  10852439.0625,
  10852857.8125,
  10852867.1875,
  10853134.375,
  10853429.6875,
  10853762.5,
  10854515.625,
  10854521.875,
  10854643.75,
  10855050.0,
  10855078.125,
  10855346.875,
  10855401.5625,
  10859253.125,
  10859296.875,
  10859881.25,
  10860401.5625,
  10861117.1875,
  10861782.8125,
  10864178.125,
  10864331.25,
  10864904.6875,
  10865728.125,
  10865762.5,
  10866118.75,
  10866176.5625,
  10867034.375,
  10867429.6875,
  10867537.5,
  10867551.5625,
  10869642.1875,
  10869682.8125,
  10870432.8125,
  10871084.375,
  10875185.9375,
  10878128.125,
  10879346.875,
  10879662.5,
  10881098.4375,
  10881107.8125,
  10881153.125,
  10881179.6875,
  10881181.25,
  10881235.9375,
  10881264.0625,
  10881396.875,
  10881445.3125,
  10881451.5625,
  10881667.1875,
  10881771.875,
  10881839.0625,
  10881996.875,
  10882318.75,
  10882718.75,
  10882768.75,
  10882801.5625,
  10882878.125,
  10882926.5625,
  10882973.4375,
  10883081.25,
  10883118.75,
  10883125.0,
  10883190.625,
  10883379.6875,
  10883543.75,
  10883570.3125,
  10883667.1875,
  10883842.1875,
  10883878.125,
  10883901.5625,
  10883932.8125,
  10884385.9375,
  10884482.8125,
  10884726.5625,
  10885081.25,
  10885246.875,
  10885251.5625,
  10885273.4375,
  10885353.125,
  10885370.3125,
  10885437.5,
  10885493.75,
  10885643.75,
  10885656.25,
  10885673.4375,
  10885810.9375,
  10885904.6875,
  10886064.0625,
  10886087.5,
  10886093.75,
  10886135.9375,
  10886206.25,
  10886273.4375,
  10886292.1875,
  10886309.375,
  10886403.125,
  10886456.25,
  10886573.4375,
  10886714.0625,
  10886735.9375,
  10886764.0625,
  10886889.0625,
  10886921.875,
  10886973.4375,
  10886995.3125,
  10887090.625,
  10887093.75,
  10887131.25,
  10887260.9375,
  10887359.375,
  10887360.9375,
  10887437.5,
  10887570.3125,
  10887584.375,
  10887681.25,
  10887714.0625,
  10887718.75,
  10887746.875,
  10887750.0,
  10887770.3125,
  10887784.375,
  10887796.875,
  10887826.5625,
  10887960.9375,
  10888059.375,
  10888104.6875,
  10888156.25,
  10888250.0,
  10888350.0,
  10888368.75,
  10888415.625,
  10888465.625,
  10888537.5,
  10888548.4375,
  10888598.4375,
  10888628.125,
  10888739.0625,
  10888742.1875,
  10888746.875,
  10888860.9375,
  10888881.25,
  10888940.625,
  10888989.0625,
  10889076.5625,
  10889190.625,
  10889231.25,
  10889318.75,
  10889417.1875,
  10889478.125,
  10889490.625,
  10889529.6875,
  10889539.0625,
  10889576.5625,
  10889642.1875,
  10889710.9375,
  10889729.6875,
  10889800.0,
  10889804.6875,
  10889917.1875,
  10890192.1875,
  10890256.25,
  10890264.0625,
  10890296.875,
  10890343.75,
  10890446.875,
  10890506.25,
  10890820.3125,
  10891071.875,
  10891073.4375,
  10891145.3125,
  10891242.1875,
  10891323.4375,
  10891504.6875,
  10891531.25,
  10891567.1875,
  10891675.0,
  10891884.375,
  10892540.625,
  10893507.8125,
  10894684.375,
  10894834.375,
  10895203.125,
  10895268.75,
  10897412.5,
  10899225.0,
  10899317.1875,
  10902109.375,
  10906518.75,
  10911089.0625,
  10912175.0,
  10913253.125,
  10914320.3125,
  10914378.125,
  10914632.8125,
  10914775.0,
  10915095.3125,
  10915164.0625,
  10915228.125,
  10916170.3125,
  10916323.4375,
  10916429.6875,
  10916771.875,
  10917215.625,
  10917495.3125,
  10917556.25,
  10917664.0625,
  10918315.625,
  10918470.3125,
  10918898.4375,
  10919525.0,
  10919698.4375,
  10919989.0625,
  10920093.75,
  10920318.75,
  10920575.0,
  10920596.875,
  10920664.0625,
  10921070.3125,
  10921278.125,
  10921353.125,
  10921356.25,
  10921434.375,
  10921645.3125,
  10922104.6875,
  10922109.375,
  10922129.6875,
  10922198.4375,
  10922257.8125,
  10922264.0625,
  10922293.75,
  10922365.625,
  10922506.25,
  10922531.25,
  10922537.5,
  10922801.5625,
  10922806.25,
  10922835.9375,
  10922892.1875,
  10923278.125,
  10923501.5625,
  10923575.0,
  10923868.75,
  10923948.4375,
  10924093.75,
  10924262.5,
  10924364.0625,
  10924453.125,
  10924465.625,
  10925006.25,
  10925051.5625,
  10925153.125,
  10925348.4375,
  10925429.6875,
  10925596.875,
  10925925.0,
  10926098.4375,
  10926128.125,
  10926151.5625,
  10926454.6875,
  10926654.6875,
  10926770.3125,
  10926889.0625,
  10926915.625,
  10927321.875,
  10927596.875,
  10927726.5625,
  10927981.25,
  10928145.3125,
  10928543.75,
  10928789.0625,
  10929031.25,
  10929031.25,
  10929532.8125,
  10929587.5,
  10930218.75,
  10930218.75,
  10930232.8125,
  10930420.3125,
  10930454.6875,
  10930664.0625,
  10930860.9375,
  10930870.3125,
  10930900.0,
  10930928.125,
  10931104.6875,
  10931218.75,
  10931448.4375,
  10931451.5625,
  10931681.25,
  10931843.75,
  10932384.375,
  10932676.5625,
  10932721.875,
  10932757.8125,
  10932917.1875,
  10934662.5,
  10934687.5,
  10935039.0625,
  10935148.4375,
  10935503.125,
  10936704.6875,
  10937028.125,
  10937239.0625,
  10937546.875,
  10937560.9375,
  10937981.25,
  10938842.1875,
  10939871.875,
  10940214.0625,
  10940606.25,
  10940781.25,
  10940810.9375,
  10940960.9375,
  10941887.5,
  10942423.4375,
  10943606.25,
  10943848.4375,
  10944360.9375,
  10945128.125,
  10945957.8125,
  10946528.125,
  10946817.1875,
  10947159.375,
  10947937.5,
  10949028.125,
  10949226.5625,
  10949809.375,
  10950296.875,
  10950470.3125,
  10950673.4375,
  10951807.8125,
  10952471.875,
  10952725.0,
  10952778.125,
  10952784.375,
  10952873.4375,
  10952900.0,
  10953121.875,
  10954700.0,
  10955804.6875,
  10955865.625,
  10956117.1875,
  10956343.75,
  10956437.5,
  10958182.8125,
  10960417.1875,
  10961435.9375,
  10961564.0625,
  10961642.1875,
  10961832.8125,
  10961914.0625,
  10961915.625,
  10963120.3125,
  10963295.3125,
  10963384.375,
  10963650.0,
  10963768.75,
  10963793.75,
  10964090.625,
  10965314.0625,
  10965384.375,
  10965462.5,
  10966151.5625,
  10967268.75,
  10967918.75,
  10969079.6875,
  10969656.25,
  10969787.5,
  10971337.5,
  10972368.75,
  10972854.6875,
  10973035.9375,
  10973556.25,
  10973951.5625,
  10974559.375,
  10974835.9375,
  10975565.625,
  10975734.375,
  10977785.9375,
  10978518.75,
  10978632.8125,
  10978659.375,
  10978737.5,
  10979903.125,
  10980192.1875,
  10980512.5,
  10981354.6875,
  10982539.0625,
  10982834.375,
  10982890.625,
  10983115.625,
  10983232.8125,
  10983560.9375,
  10983590.625,
  10984081.25,
  10984831.25,
  10984979.6875,
  10985089.0625,
  10985362.5,
  10985415.625,
  10985467.1875,
  10985660.9375,
  10985714.0625,
  10985831.25,
  10986034.375,
  10986057.8125,
  10986114.0625,
  10986234.375,
  10986971.875,
  10987185.9375,
  10987251.5625,
  10987528.125,
  10987803.125,
  10987950.0,
  10988021.875,
  10988181.25,
  10988284.375,
  10988796.875,
  10988881.25,
  10989420.3125,
  10990409.375,
  10990484.375,
  10990981.25,
  10990982.8125,
  10991706.25,
  10992010.9375,
  10992103.125,
  10992479.6875,
  10992546.875,
  10992654.6875,
  10992743.75,
  10992884.375,
  10993004.6875,
  10993120.3125,
  10993168.75,
  10993217.1875,
  10993271.875,
  10993303.125,
  10993303.125,
  10993315.625,
  10993532.8125,
  10994193.75,
  10994250.0,
  10994668.75,
  10994884.375,
  10995017.1875,
  10995057.8125,
  10995295.3125,
  10995485.9375,
  10995782.8125,
  10995939.0625,
  10997456.25,
  10998806.25,
  11000315.625,
  11001193.75,
  11001395.3125,
  11001623.4375,
  11002075.0,
  11002526.5625,
  11002592.1875,
  11002617.1875,
  11002956.25,
  11003267.1875,
  11006328.125,
  11006500.0,
  11006693.75,
  11007612.5,
  11008003.125,
  11008017.1875,
  11008315.625,
  11008356.25,
  11009971.875],
 [22.021009242355184,
  38.4890783953468,
  47.02068742803747,
  28.392128982770146,
  66.74958698675533,
  55.22233102952246,
  5.767541976575615,
  85.98545176225845,
  76.9788880898114,
  25.313246423077054,
  26.001943124255078,
  92.83992457127137,
  22.140125767721898,
  6.775345668204583,
  23.317082582744632,
  29.846308908107645,
  60.91451679256374,
  7.328955250217951,
  10.965439446176706,
  54.96693523938706,
  5.0788987163568,
  49.9400207927214,
  8.494233821597778,
  17.79456726236055,
  26.614098956838923,
  24.077513630222768,
  7.881532955875548,
  51.47549429520111,
  38.09964661220158,
  108.72243320993135,
  48.32018342004417,
  12.352624450128737,
  15.216745508036997,
  47.56984033199902,
  22.11036563746642,
  14.560200383213502,
  13.21185426504061,
  17.00803663695596,
  77.58990132771481,
  18.928524810622605,
  18.26483942357279,
  5.528073577416909,
  27.787439662382923,
  15.147994241430474,
  77.08979221999557,
  44.22948272239949,
  5.82832345829075,
  41.997380502800254,
  121.56549943227482,
  89.5261471073743,
  16.279487112996037,
  18.588269635683083,
  31.178240603922152,
  133.54967171658197,
  46.77596037010173,
  26.653184434819728,
  79.97654101549485,
  65.43011049508247,
  17.955891771497473,
  61.90826001779653,
  20.868631608281472,
  30.946566757831633,
  61.02832664774631,
  165.8110614041866,
  53.76319944199383,
  20.34710226503025,
  26.1027783369538,
  18.50971128519458,
  12.056668303282246,
  17.21388380696328,
  89.09662348260042,
  11.821077989476422,
  200.72604604634267,
  5.4336971123951825,
  16.10064007640934,
  15.351378248519651,
  66.85188253566955,
  94.70977099622877,
  38.94508089319679,
  7.364347189617007,
  56.823631086798514,
  6.9515240150768,
  23.655785158591573,
  109.94492113787433,
  22.57318643421639,
  28.80568689671915,
  36.86961944261854,
  7.283532726832849,
  25.09127887004715,
  46.893205239032824,
  29.074893782244324,
  137.60778568136072,
  45.539803689582556,
  34.73445517965499,
  98.36953253093536,
  5.276151751162502,
  16.73576318457908,
  20.490710855333738,
  11.666800405834781,
  40.666781069542665,
  95.499035707235,
  60.32844883773057,
  53.09497305128083,
  11.138284956021902,
  34.41457539434679,
  18.538386078888358,
  47.81573168259413,
  13.823627868534967,
  88.07231232320218,
  16.841576235918975,
  6.313443538055007,
  46.01188192946974,
  35.4555126037289,
  62.31103758371546,
  18.260874149808338,
  70.50340101649549,
  16.998355474105303,
  32.349993325646366,
  36.35025351582656,
  9.905739734247941,
  6.272346554264856,
  16.548075951133608,
  173.4780343651496,
  14.66564857523588,
  18.176815182647395,
  25.105044093816872,
  16.020130987277902,
  19.486453323053702,
  20.430833529285074,
  27.70935736374711,
  36.63704332839521,
  27.408499881318978,
  19.671181140231766,
  11.382255464408123,
  16.694242426628662,
  72.4963068669911,
  16.224052722840963,
  32.937620370474534,
  13.054374663405506,
  53.593902462809424,
  37.660485575055304,
  124.6842273463259,
  11.973016827946285,
  20.17553459537155,
  24.393674603768687,
  5.041678731527717,
  23.24902946044456,
  76.2990777405718,
  18.977857621516563,
  62.036722606662316,
  7.0193588417710515,
  10.981624497878927,
  8.63548794277709,
  19.284447395602363,
  91.59509239318692,
  5.239379273968614,
  18.311129497436117,
  83.30861089082661,
  68.88100814463934,
  11.500091722865413,
  81.27556251771465,
  38.88334790456944,
  17.95606701292393,
  34.5942353457635,
  11.094001898801896,
  30.218668317589874,
  28.793985019134244,
  35.00334716868556,
  89.27003577678465,
  62.178988330270435,
  5.786125626670979,
  15.4954564536649,
  14.245958827256203,
  12.20793513230758,
  17.964454316839333,
  98.4255103702047,
  17.351265301932973,
  7.86100509754009,
  22.050736113141603,
  26.438877054565015,
  88.08701175903134,
  13.770208485666899,
  88.72686862228527,
  163.30563692716265,
  52.090468440495556,
  21.586448490358404,
  9.474736021000595,
  36.213764981720736,
  85.08699727616452,
  214.53261368846267,
  10.740221575344053,
  35.77077205935591,
  91.18864709509413,
  16.229036704589458,
  10.799802374587344,
  94.71323087772357,
  39.0077658431792,
  5.621626020592023,
  21.458026174122292,
  39.61612410696293,
  23.427202472568087,
  30.099284242926736,
  157.64861930971978,
  30.809233968790622,
  5.8023261510349915,
  13.519787243402202,
  13.73244761866103,
  14.153340818472302,
  17.512550123563887,
  5.799873310174771,
  65.81773487747353,
  31.09619083558779,
  24.440747308417343,
  144.62507467580127,
  65.51660079114662,
  33.451659394009475,
  8.00529218687174,
  59.24554381775309,
  8.748882696933995,
  55.41716628969277,
  5.831643974250857,
  223.5296764355819,
  30.710301430118875,
  8.697083915282567,
  30.55764708558197,
  29.11562981964126,
  7.71821965787399,
  5.332107034424148,
  8.953723696357054,
  12.740229987697337,
  14.068827868213816,
  36.27327213427412,
  32.020303071632846,
  44.60634230266353,
  16.57414431456052,
  6.434576193701606,
  14.907339095605545,
  22.890187732247284,
  5.24721410341351,
  25.26679854568272,
  15.271712538456171,
  36.6513410796253,
  5.383907988892966,
  24.901555215662682,
  25.036365624854696,
  139.8300860633072,
  9.680220528056239,
  6.11022503291402,
  15.076892778596527,
  31.1540587976988,
  22.547681036207095,
  11.572712877285262,
  65.87077213573423,
  37.741516255717784,
  98.86840405537988,
  8.814990455582633,
  41.37817426786875,
  5.204984731342885,
  41.52914443715207,
  106.0944401293016,
  15.021078390882824,
  22.82330875321483,
  16.66503134240758,
  64.0284232252087,
  43.765999967727694,
  56.41239548588396,
  31.73265188069879,
  31.92621215244271,
  39.5202020939635,
  12.743189561458843,
  19.349836969610667,
  17.93846462401773,
  49.29586602947651,
  15.279825888494805,
  34.50729847707833,
  77.14075514074692,
  37.47004128373334,
  10.928065331906891,
  22.9258506804753,
  151.8548261866377,
  50.36576040721054,
  57.160279098744134,
  5.1463006665050255,
  16.246580314279104,
  37.54345238010842,
  36.79520231005443,
  26.596550451061262,
  30.02696603035968,
  30.29493832764156,
  66.75335578822245,
  22.47519437477357,
  20.362495308891447,
  13.782695860583036,
  5.105228635646729,
  11.465829376306463,
  24.689489332707915,
  22.468568125350938,
  18.251894050617672,
  10.41531876261297,
  13.710940199974536,
  52.34867932798088,
  42.540561392271215,
  9.865211756656256,
  19.86159492684419,
  13.716269438775196,
  8.12563030131686,
  6.079465494844198,
  24.220883659054937,
  7.680922032547267,
  30.718974200800645,
  7.282300362708002,
  67.03591453283133,
  44.46697470882711,
  9.25112182635102,
  35.64006940966409,
  7.054792151097475,
  36.9233388624636,
  38.408505668803954,
  31.26842943641544,
  44.97547688189459,
  72.49704185320105,
  18.03348940304471,
  128.7066757451429,
  134.79254782541915,
  9.099341105235155,
  28.009804625469606,
  6.209783846193417,
  65.0208204057012,
  13.855440865426292,
  115.08518137467516,
  37.232730994110085,
  10.031061446947488,
  59.63945759145744,
  76.92252500768866,
  41.165820552689155,
  95.07188602774748,
  24.657813383502127,
  11.137720828219305,
  10.496877095393701,
  25.47161697023405,
  55.64807023954094,
  31.104682557426386,
  13.26363570340616,
  43.488304977597295,
  13.233473167898177,
  15.360079603334214,
  22.848670904501805,
  42.559909705301784,
  34.67602754924622,
  37.97093729504877,
  89.55418966788866,
  23.655749886430065,
  25.256919586640304,
  57.95666255489811,
  26.289583863198995,
  69.03233943083873,
  87.2797337057543,
  6.682150561842103,
  24.833240771015063,
  47.868893005645724,
  31.129884345645216,
  68.4942083232207,
  5.181325306025436,
  37.8727256560143,
  69.26176883543641,
  43.32495886566604,
  43.07275486455848,
  15.867927464700793,
  56.842124237070955,
  11.026938422475894,
  29.742461821466048,
  35.477831521480155,
  109.3976900512591,
  22.06295365220376,
  34.197541967758205,
  16.91099457022679,
  23.679096116763116,
  56.24427229374366,
  13.056530399021842,
  5.245650491351172,
  16.585874732501,
  9.519687281671542,
  5.3110270733573826,
  49.98210644744865,
  10.533285023869338,
  14.03520230293683,
  80.53185550481206,
  142.20907185489466,
  13.621867184182264,
  10.496795945218656,
  52.794487016302504,
  5.81852381302713,
  5.11988686230411,
  13.010602722255353,
  72.03845711836274,
  19.96349716340491,
  7.4510984617700595,
  6.519807299136634,
  76.72752828676161,
  14.955102447042108,
  17.488279106598295,
  13.2240439570794,
  77.17335669046227,
  33.601244916449915,
  43.65337426395233,
  10.535746741680837,
  21.54555258334053,
  20.873519761088723,
  28.238810736462195,
  5.778985931681963,
  7.711172389454429,
  52.85776382538101,
  81.54478801636546,
  9.975727577410039,
  5.156266483097807,
  88.13017486023823,
  17.535101307001124,
  37.6776960318424,
  62.6841803427774,
  14.224764176988328,
  18.411617199313568,
  41.170161100909084,
  15.974723954580107,
  21.634210021222763,
  35.126478197135846,
  29.38502145942212,
  12.113958991867475,
  91.67982697504785,
  38.62457881105571,
  31.782638781353505,
  64.6841965692721,
  11.417150361034581,
  73.87486338456806,
  8.806014972258357,
  18.95679500699656,
  39.14002171153467,
  22.41591445096192,
  65.99529546426692,
  28.67072992745748,
  8.805729919763305,
  65.99059132796356,
  28.11468680341484,
  76.44349714650832,
  66.16629340087987,
  10.152127676087854,
  6.357698827331542,
  84.66973836414083,
  10.592715304784974,
  65.64720746561785,
  47.66557987703953,
  23.33183279362798,
  44.579767139026444,
  16.380984123555233,
  7.905809331005398,
  65.33143082728952,
  38.16439296293021,
  61.836374435680966,
  37.69603815554154,
  22.733831852294873,
  10.911905571870372,
  74.65633083934625,
  12.508210428003846,
  9.725340240774221,
  48.09394253109295,
  97.45586146597684,
  84.18692986328307,
  7.898691548649248,
  10.68649086219551,
  37.46424444986359,
  24.73213961121516,
  41.872491913159266,
  40.62338597392003,
  17.05929120717752,
  12.748262745258055,
  9.773554885743303,
  28.227988516823736,
  5.433260964862621,
  17.68868122704911,
  42.4351797973745,
  48.472047921834985,
  9.856016171300812,
  9.06385286707796,
  18.72618452329767,
  17.19172638298216,
  77.62683033630735,
  5.350408738443339,
  81.43577629676666,
  57.4512005163295,
  65.90236095603541,
  28.558262024427805,
  36.25644384351509,
  114.36925168269248,
  16.345061629878007,
  78.91918042756956,
  5.3961460087826705,
  184.18845408512198,
  5.44601500182603,
  18.178213819118735,
  38.768675227088394,
  57.638076333475645,
  34.11240184579252,
  42.55537570015323,
  14.561482675649726,
  40.89795889076275,
  10.857965566203532,
  87.63891674092557,
  73.51387920871596,
  17.6120719174898,
  5.249086831333948,
  34.32455983600745,
  5.171124406431231,
  9.094198901350126,
  27.623466950636892,
  45.21400875983571,
  28.284419286661436,
  111.75056451168965,
  33.67624455001214,
  6.813122942536255,
  6.112532485850941,
  83.34096537565901,
  97.63111396509836,
  38.36411684364267,
  79.12613931639416,
  19.856801445326067,
  63.479904976588095,
  99.23966045508843,
  5.716553670554265,
  203.4042285228039,
  36.85084191552412,
  13.847658985820932,
  42.40245615647174,
  7.7619145515877115,
  10.125094502875928,
  22.59347082126714,
  35.5995571271542,
  49.83339886896031,
  13.023547661770873,
  86.32743717240116,
  5.365018795285513,
  30.15722298452897,
  53.3130026198923,
  40.11615393095157,
  104.90759889866919,
  22.038872957895677,
  111.06000894916396,
  86.21779869376928,
  16.648411062053647,
  15.122188221791173,
  119.67093414439576,
  56.10826696282055,
  30.4619104345695,
  85.57323752579141,
  26.64526035465237,
  5.512210331251039,
  16.168188756214402,
  50.16706543608371,
  7.862688686332413,
  21.174991404603016,
  29.770891245999692,
  53.83780875022511,
  45.27308263126477,
  18.849458107855693,
  22.682093444403883,
  9.267214533984296,
  84.65023467444472,
  19.710697184538045,
  22.2028086579542,
  29.509347136265312,
  6.933488576430554,
  22.683479401268528,
  5.271866949848501,
  48.45177151309455,
  94.34272152631485,
  50.74603832448135,
  65.8893793756094,
  59.18668021183142,
  6.411250187473843,
  34.330131028621736,
  40.6576972695079,
  11.104617550258228,
  7.887705892111198,
  5.18005252761478,
  19.746549438652668,
  92.05133031035159,
  35.05692637374621,
  45.4773152326746,
  23.73527356299942,
  22.16014330001128,
  7.680815981306755,
  28.249317472032615,
  94.17449764976114,
  12.395297198884663,
  65.55968764439035,
  8.188548098923544,
  19.76390527371317,
  13.782718952804494,
  27.463261022589784,
  7.253873108500993,
  7.098247328973293,
  87.34684372462829,
  66.49412408220236,
  7.283875763722794,
  17.330554785575572,
  36.52366381780125,
  8.01223805647984,
  56.429363548751994,
  54.06117003741237,
  14.770592630022241,
  26.42326911590387,
  59.868183019168306,
  43.23805501909495,
  10.214640248906216,
  59.72851590158683,
  20.659875269530673,
  23.493264175909378,
  31.204148809644746,
  25.198365358576496,
  10.903033004281209,
  51.03699798151449,
  26.935533597535574,
  52.816519902510024,
  28.4479667681646,
  13.573508083331795,
  13.31283556640876,
  59.11371970603487,
  27.41232303341799,
  63.92876644192435,
  16.3301677067414,
  30.17306679133136,
  73.05017821255817,
  43.76804684322735,
  9.377318534420954,
  101.34424601873934,
  41.035859801503136,
  24.265292619418787,
  8.277266186050756,
  34.45577091187768,
  63.05802953976904,
  11.163152759147,
  15.490200850043493,
  7.000659673439083,
  18.65522636836891,
  8.611027181783722,
  15.7071408820859,
  15.671686514707517,
  15.551000979005742,
  55.86499579745119,
  21.053437689644273,
  73.85920695239307,
  21.616976840544154,
  113.0408594774431,
  15.286231190939306,
  24.26820016859054,
  35.35622195433799,
  65.41433752029411,
  59.08459494152116,
  23.314415555225686,
  5.617862029888876,
  12.961658684330763,
  14.75647783470239,
  8.73183091086844,
  29.614177463402136,
  20.342343756181247,
  30.213886107224802,
  5.058864945783083,
  17.524693815422285,
  25.07481566810459,
  9.215964805219468,
  15.434852025089148,
  10.798638114703666,
  24.564514480703515,
  25.390587529801348,
  18.9742667465705,
  18.908219815034375,
  30.34449822091364,
  14.624239766767865,
  46.16688792442261,
  20.307999884402832,
  12.946293858889382,
  48.208410674804554,
  5.062074776884923,
  16.95465162004606,
  16.6640376812401,
  47.45443784378364,
  19.370594172227264,
  29.58239319234457,
  65.03935836284646,
  5.339394518650768,
  5.874130199960952,
  22.421037076783215,
  46.58539084372967,
  7.990716144309684,
  22.979159192473006,
  7.193342301813022,
  56.42877660935584,
  38.63725540578667,
  21.01436858523974,
  14.171886003128202,
  63.64548751945541,
  29.664315790507796,
  74.3046080298979,
  18.35924808820434,
  14.472804136288927,
  60.94876527352899,
  14.602847740925831,
  116.40782153342616,
  8.992212337913951,
  5.9875307768331725,
  78.45442709647688,
  74.62526212553288,
  46.40240576676584,
  12.054653876271365,
  5.285194301886191,
  7.839299582827493,
  63.392201058440435,
  54.77349197330699,
  12.847227483279545,
  71.80080956859267,
  14.586885534611513,
  11.755643620206742,
  16.436451660646064,
  10.46448021601351,
  16.851498424285793,
  8.626307341908115,
  35.91879378318802,
  46.550582331435805,
  8.046991122943412,
  55.17416982013076,
  46.85609323020705,
  123.32736537632918,
  10.134955948997298,
  29.108489257247527,
  21.27870430558702,
  25.34670583857635,
  49.5673741311398,
  71.76938397849558,
  23.111600725871835,
  50.19117509704823,
  8.680596800741602,
  55.06642405893906,
  10.329415999954492,
  90.42979675993865,
  29.71168000827199,
  23.82351004444102,
  35.71374341061582,
  5.510616214255254,
  55.83269388527026,
  5.203303906647725,
  7.079532911234448,
  10.34548884351672,
  23.688498973806844,
  44.93005624181225,
  47.62743080132242,
  19.80910117496984,
  15.33786883264623,
  61.1266382591012,
  43.013515943568244,
  124.00386625107237,
  52.01179770924548,
  81.10372725906613,
  15.249657599687554,
  12.384611721622107,
  17.71189451480888,
  34.728231331536676,
  20.475720516204397,
  65.64446515560198,
  9.566740644794663,
  11.445634567423177,
  68.9017362685597,
  29.98705052870775,
  11.139339005510818,
  13.30437144220927,
  27.846160543068166,
  5.674755432438352,
  12.343400811292362,
  16.94600758369818,
  20.42515646291394,
  29.364399613952926,
  61.024032258373,
  7.105237067849262,
  11.936150474254427,
  20.701798610012098,
  27.923475131570207,
  21.169829402338863,
  10.678265085728777,
  103.51254346744452,
  72.44795854086237,
  10.915504204385975,
  6.776016059219304,
  17.290112236251346,
  14.548129569537364,
  21.262300580920257,
  14.650993960765222,
  19.795831747331114,
  5.928722362066054,
  5.715891480488105,
  20.707832994601986,
  17.599363625216174,
  13.4017068908327,
  14.668583908856826,
  25.27435717097639,
  20.64032799362026,
  81.86736506133546,
  15.225643744904625,
  87.47437020994107,
  29.30794757172552,
  18.102969558140508,
  8.438859460490141,
  56.3849649656277,
  34.13813401389431,
  36.0818578019151,
  86.66110046385121,
  20.708465544551917,
  77.10642400050229,
  44.34881107527772,
  14.713167720706386,
  71.51845277837111,
  22.375394475295273,
  16.28818993360847,
  48.57040585597619,
  41.96772091864355,
  22.283426862356265,
  71.72591491862836,
  20.626442323309714,
  84.91838809504783,
  5.292777944683193,
  13.336394033125114,
  44.48343760206406,
  57.336228013175926,
  5.410824069280008,
  32.36725898468573,
  85.74413940363755,
  13.96525748953117,
  36.63488752594296,
  22.31510641626128,
  32.39794798722611,
  17.56963513315327,
  73.29718922332276,
  28.390550148637203,
  78.737325721106,
  40.13258517613334,
  24.2569644572232,
  5.73677101680262,
  17.007721219121464,
  27.546531131391724,
  17.88566874786028,
  29.04182172084941,
  5.098931286238759,
  5.362676915658964,
  60.05818336903028,
  21.74389677485656,
  37.48173029369633,
  7.21879373844304,
  44.515530418320196,
  142.56983849035865,
  18.138350794524957,
  7.616096618973378,
  14.289807409547906,
  35.20519866100101,
  27.285240800385132,
  14.368589102463558,
  51.62647468520487,
  25.991859649822743,
  9.074751637838316,
  9.046592571905157,
  14.584119126960067,
  5.937328771401298,
  10.045440512891627,
  23.313660511017435,
  17.107995963080118,
  5.818522107576128,
  71.08696593997973,
  25.442802520418088,
  41.45241381870106,
  66.22400949133461,
  146.57638690000948,
  105.88807152064071,
  5.960893475895645,
  6.148202775042389,
  21.192313707992646,
  44.09465620369771,
  25.457030630336206,
  12.425603059532982,
  24.6471178701932,
  35.382632947966556,
  17.186985693617892,
  5.161512798362923,
  11.884989819364487,
  5.0549004826156985,
  6.386668212430329,
  15.697027261721002,
  16.242766529518974,
  10.170805794759305,
  17.497511785217142,
  78.80800858733525,
  88.79806092413611,
  23.144164809664098,
  13.492445974154005,
  50.9021917385543,
  52.65022817104952,
  121.10695939461881,
  44.3338198699621,
  15.139243419728508,
  46.04415661762758,
  7.385940347398921,
  7.797147971087979,
  79.39219995992596,
  21.09274255073534,
  25.550954959178696,
  8.575364572269729,
  69.00936254497469,
  82.56555452336234,
  11.766986380727856,
  67.68651563337193,
  69.35916430244502,
  79.07667980291255,
  25.504299287895186,
  53.09170737629367,
  12.84771979132092,
  64.6783645165732,
  30.29273264034638,
  20.994389289618372,
  101.47276786301427,
  6.200436275231995,
  27.108089089319712,
  22.83534139160889,
  6.223195707648995,
  14.126174705444814,
  9.093728426326399,
  7.587944510358162,
  66.22406423386317,
  71.55801092103763])

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

([3095096.875,
  3130134.375,
  3173053.125,
  3184339.0625,
  3236659.375,
  3339829.6875,
  3361368.75,
  3394290.625,
  3478837.5,
  3623385.9375,
  3630409.375,
  3638814.0625,
  3644310.9375,
  3749187.5,
  3800285.9375,
  3825940.625,
  3833795.3125,
  3835598.4375,
  3839553.125,
  3881589.0625,
  3888871.875,
  3936321.875,
  3976903.125,
  4060014.0625,
  4067768.75,
  4092753.125,
  4106103.125,
  4359567.1875,
  4518362.5,
  5543140.625,
  5544182.8125,
  5604010.9375,
  5605882.8125,
  5606962.5,
  5812875.0,
  5815982.8125,
  5816301.5625,
  5817520.3125,
  5822973.4375,
  5824389.0625,
  5827267.1875,
  5828703.125,
  5879342.1875,
  5885393.75,
  5895771.875,
  5913675.0,
  5924465.625,
  5924928.125,
  5952151.5625,
  5955779.6875,
  5969759.375,
  5974773.4375,
  5975303.125,
  5987084.375,
  6001139.0625,
  6002212.5,
  6003457.8125,
  6004635.9375,
  6005606.25,
  6010192.1875,
  6010570.3125,
  6038145.3125,
  6044790.625,
  6054506.25,
  6057182.8125,
  6060404.6875,
  6061629.6875,
  6061817.1875,
  6061825.0,
  6061843.75,
  6062412.5,
  6063631.25,
  6067809.375,
  6069576.5625,
  6139709.375,
  6140701.5625,
  6146695.3125,
  6148635.9375,
  6154395.3125,
  6226862.5,
  6230390.625,
  6234929.6875,
  6271028.125,
  6281798.4375,
  6281946.875,
  6312445.3125,
  6349659.375,
  6354142.1875,
  6354982.8125,
  6360348.4375,
  6366464.0625,
  6429153.125,
  6430348.4375,
  6430471.875,
  6430485.9375,
  6430554.6875,
  6476718.75,
  6481134.375,
  6491523.4375,
  6543426.5625,
  6549515.625,
  6562106.25,
  6574489.0625,
  6574529.6875,
  6585251.5625,
  6594785.9375,
  6597854.6875,
  6652389.0625,
  6660606.25,
  6667795.3125,
  6686635.9375,
  6709715.625,
  6909517.1875,
  6920348.4375,
  6920362.5,
  7085657.8125,
  7090815.625,
  7262906.25,
  7471714.0625,
  7512268.75,
  7623790.625,
  7670892.1875,
  7702435.9375,
  7746726.5625,
  7747443.75,
  7747725.0,
  7821670.3125,
  8293815.625,
  8407873.4375,
  8549767.1875,
  8581109.375,
  8673853.125,
  8674500.0,
  8844775.0,
  8883093.75,
  8958075.0,
  9048804.6875,
  9134260.9375,
  9166535.9375,
  9167731.25,
  9185896.875,
  9252093.75,
  9264071.875,
  9267556.25,
  9276921.875,
  9277060.9375,
  9287301.5625,
  9293276.5625,
  9295248.4375,
  9298473.4375,
  9298865.625,
  9299448.4375,
  9301309.375,
  9301535.9375,
  9302603.125,
  9302654.6875,
  9307309.375,
  9308431.25,
  9308510.9375,
  9322225.0,
  9325959.375,
  9327304.6875,
  9327531.25,
  9328395.3125,
  9330585.9375,
  9331032.8125,
  9331407.8125,
  9332526.5625,
  9332551.5625,
  9338215.625,
  9338270.3125,
  9339307.8125,
  9342260.9375,
  9344514.0625,
  9346723.4375,
  9347164.0625,
  9348017.1875,
  9351326.5625,
  9352309.375,
  9352423.4375,
  9358867.1875,
  9380614.0625,
  9384768.75,
  9385475.0,
  9574679.6875,
  9588153.125,
  9589546.875,
  9596085.9375,
  9597867.1875,
  9598028.125,
  9610746.875,
  9645923.4375,
  9743623.4375,
  9774156.25,
  9807400.0,
  9915476.5625,
  10069742.1875,
  10158950.0,
  10280985.9375,
  10289282.8125,
  10295045.3125,
  10295084.375,
  10301118.75,
  10307501.5625,
  10308054.6875,
  10308654.6875,
  10309025.0,
  10309135.9375,
  10309784.375,
  10309954.6875,
  10310046.875,
  10310456.25,
  10310623.4375,
  10311826.5625,
  10312025.0,
  10312292.1875,
  10312490.625,
  10312698.4375,
  10314018.75,
  10314235.9375,
  10314509.375,
  10314610.9375,
  10314942.1875,
  10317314.0625,
  10317371.875,
  10317496.875,
  10317634.375,
  10317743.75,
  10317978.125,
  10318615.625,
  10319348.4375,
  10319507.8125,
  10320118.75,
  10320273.4375,
  10320300.0,
  10320495.3125,
  10320731.25,
  10320787.5,
  10321564.0625,
  10322139.0625,
  10322340.625,
  10323145.3125,
  10323759.375,
  10324007.8125,
  10324859.375,
  10324929.6875,
  10325692.1875,
  10325932.8125,
  10326228.125,
  10326500.0,
  10327004.6875,
  10327128.125,
  10327271.875,
  10327967.1875,
  10328482.8125,
  10328885.9375,
  10328943.75,
  10329042.1875,
  10329334.375,
  10329971.875,
  10329978.125,
  10331275.0,
  10331809.375,
  10332126.5625,
  10332468.75,
  10332837.5,
  10333168.75,
  10333214.0625,
  10333496.875,
  10334393.75,
  10334576.5625,
  10335740.625,
  10335789.0625,
  10335893.75,
  10336884.375,
  10337393.75,
  10337604.6875,
  10338650.0,
  10339279.6875,
  10339659.375,
  10339803.125,
  10339995.3125,
  10340067.1875,
  10340114.0625,
  10340175.0,
  10340826.5625,
  10341267.1875,
  10342056.25,
  10343078.125,
  10343800.0,
  10343871.875,
  10344368.75,
  10345551.5625,
  10345629.6875,
  10345679.6875,
  10345717.1875,
  10346439.0625,
  10346998.4375,
  10350035.9375,
  10350800.0,
  10351921.875,
  10352040.625,
  10352481.25,
  10352781.25,
  10353268.75,
  10355371.875,
  10355970.3125,
  10356046.875,
  10358014.0625,
  10359557.8125,
  10360114.0625,
  10360432.8125,
  10360439.0625,
  10361782.8125,
  10361954.6875,
  10362314.0625,
  10362467.1875,
  10369871.875,
  10373509.375,
  10373914.0625,
  10388879.6875,
  10477029.6875,
  10498712.5,
  10501132.8125,
  10518271.875,
  10518284.375,
  10524175.0,
  10525876.5625,
  10528479.6875,
  10537621.875,
  10542632.8125,
  10546545.3125,
  10550665.625,
  10551973.4375,
  10568545.3125,
  10570568.75,
  10571379.6875,
  10575551.5625,
  10576414.0625,
  10583564.0625,
  10583575.0,
  10584500.0,
  10584512.5,
  10586576.5625,
  10586842.1875,
  10589575.0,
  10591045.3125,
  10593678.125,
  10593968.75,
  10599493.75,
  10604378.125,
  10606801.5625,
  10610676.5625,
  10620978.125,
  10621756.25,
  10625560.9375,
  10632734.375,
  10634217.1875,
  10637003.125,
  10641867.1875,
  10647748.4375,
  10650081.25,
  10650137.5,
  10651003.125,
  10651812.5,
  10652564.0625,
  10652648.4375,
  10652709.375,
  10653234.375,
  10653278.125,
  10660632.8125,
  10661256.25,
  10664578.125,
  10665196.875,
  10665653.125,
  10665859.375,
  10666350.0,
  10669685.9375,
  10672007.8125,
  10672059.375,
  10672292.1875,
  10674923.4375,
  10675757.8125,
  10677909.375,
  10677921.875,
  10677967.1875,
  10678023.4375,
  10680100.0,
  10683654.6875,
  10683773.4375,
  10687082.8125,
  10688095.3125,
  10689064.0625,
  10689368.75,
  10689915.625,
  10689975.0,
  10691718.75,
  10692067.1875,
  10694081.25,
  10694150.0,
  10694625.0,
  10695029.6875,
  10695581.25,
  10696026.5625,
  10697865.625,
  10698284.375,
  10698354.6875,
  10698550.0,
  10702175.0,
  10702253.125,
  10702845.3125,
  10703668.75,
  10705820.3125,
  10707839.0625,
  10709551.5625,
  10709714.0625,
  10709856.25,
  10711900.0,
  10713340.625,
  10714584.375,
  10715073.4375,
  10718389.0625,
  10722868.75,
  10801353.125,
  10801359.375,
  10805301.5625,
  10806518.75,
  10806573.4375,
  10806756.25,
  10807281.25,
  10810979.6875,
  10812654.6875,
  10821601.5625,
  10832964.0625,
  10836409.375,
  10838545.3125,
  10839285.9375,
  10839290.625,
  10840954.6875,
  10841642.1875,
  10841798.4375,
  10842746.875,
  10843996.875,
  10845215.625,
  10845837.5,
  10847076.5625,
  10847596.875,
  10847832.8125,
  10848335.9375,
  10848575.0,
  10849181.25,
  10849206.25,
  10849225.0,
  10850168.75,
  10850337.5,
  10851575.0,
  10851654.6875,
  10852021.875,
  10852251.5625,
  10852435.9375,
  10852439.0625,
  10852857.8125,
  10852867.1875,
  10853134.375,
  10853429.6875,
  10853762.5,
  10854515.625,
  10854521.875,
  10854643.75,
  10855050.0,
  10855078.125,
  10855346.875,
  10855401.5625,
  10859253.125,
  10859296.875,
  10859881.25,
  10860401.5625,
  10861117.1875,
  10861782.8125,
  10864178.125,
  10864331.25,
  10864904.6875,
  10865728.125,
  10865762.5,
  10866118.75,
  10866176.5625,
  10867034.375,
  10867429.6875,
  10867537.5,
  10867551.5625,
  10869642.1875,
  10869682.8125,
  10870432.8125,
  10871084.375,
  10875185.9375,
  10878128.125,
  10879346.875,
  10879662.5,
  10881098.4375,
  10881107.8125,
  10881153.125,
  10881179.6875,
  10881181.25,
  10881235.9375,
  10881264.0625,
  10881396.875,
  10881445.3125,
  10881451.5625,
  10881667.1875,
  10881771.875,
  10881839.0625,
  10881996.875,
  10882318.75,
  10882718.75,
  10882768.75,
  10882801.5625,
  10882878.125,
  10882926.5625,
  10882973.4375,
  10883081.25,
  10883118.75,
  10883125.0,
  10883190.625,
  10883379.6875,
  10883543.75,
  10883570.3125,
  10883667.1875,
  10883842.1875,
  10883878.125,
  10883901.5625,
  10883932.8125,
  10884385.9375,
  10884482.8125,
  10884726.5625,
  10885081.25,
  10885246.875,
  10885251.5625,
  10885273.4375,
  10885353.125,
  10885370.3125,
  10885437.5,
  10885493.75,
  10885643.75,
  10885656.25,
  10885673.4375,
  10885810.9375,
  10885904.6875,
  10886064.0625,
  10886087.5,
  10886093.75,
  10886135.9375,
  10886206.25,
  10886273.4375,
  10886292.1875,
  10886309.375,
  10886403.125,
  10886456.25,
  10886573.4375,
  10886714.0625,
  10886735.9375,
  10886764.0625,
  10886889.0625,
  10886921.875,
  10886973.4375,
  10886995.3125,
  10887090.625,
  10887093.75,
  10887131.25,
  10887260.9375,
  10887359.375,
  10887360.9375,
  10887437.5,
  10887570.3125,
  10887584.375,
  10887681.25,
  10887714.0625,
  10887718.75,
  10887746.875,
  10887750.0,
  10887770.3125,
  10887784.375,
  10887796.875,
  10887826.5625,
  10887960.9375,
  10888059.375,
  10888104.6875,
  10888156.25,
  10888250.0,
  10888350.0,
  10888368.75,
  10888415.625,
  10888465.625,
  10888537.5,
  10888548.4375,
  10888598.4375,
  10888628.125,
  10888739.0625,
  10888742.1875,
  10888746.875,
  10888860.9375,
  10888881.25,
  10888940.625,
  10888989.0625,
  10889076.5625,
  10889190.625,
  10889231.25,
  10889318.75,
  10889417.1875,
  10889478.125,
  10889490.625,
  10889529.6875,
  10889539.0625,
  10889576.5625,
  10889642.1875,
  10889710.9375,
  10889729.6875,
  10889800.0,
  10889804.6875,
  10889917.1875,
  10890192.1875,
  10890256.25,
  10890264.0625,
  10890296.875,
  10890343.75,
  10890446.875,
  10890506.25,
  10890820.3125,
  10891071.875,
  10891073.4375,
  10891145.3125,
  10891242.1875,
  10891323.4375,
  10891504.6875,
  10891531.25,
  10891567.1875,
  10891675.0,
  10891884.375,
  10892540.625,
  10893507.8125,
  10894684.375,
  10894834.375,
  10895203.125,
  10895268.75,
  10897412.5,
  10899225.0,
  10899317.1875,
  10902109.375,
  10906518.75,
  10911089.0625,
  10912175.0,
  10913253.125,
  10914320.3125,
  10914378.125,
  10914632.8125,
  10914775.0,
  10915095.3125,
  10915164.0625,
  10915228.125,
  10916170.3125,
  10916323.4375,
  10916429.6875,
  10916771.875,
  10917215.625,
  10917495.3125,
  10917556.25,
  10917664.0625,
  10918315.625,
  10918470.3125,
  10918898.4375,
  10919525.0,
  10919698.4375,
  10919989.0625,
  10920093.75,
  10920318.75,
  10920575.0,
  10920596.875,
  10920664.0625,
  10921070.3125,
  10921278.125,
  10921353.125,
  10921356.25,
  10921434.375,
  10921645.3125,
  10922104.6875,
  10922109.375,
  10922129.6875,
  10922198.4375,
  10922257.8125,
  10922264.0625,
  10922293.75,
  10922365.625,
  10922506.25,
  10922531.25,
  10922537.5,
  10922801.5625,
  10922806.25,
  10922835.9375,
  10922892.1875,
  10923278.125,
  10923501.5625,
  10923575.0,
  10923868.75,
  10923948.4375,
  10924093.75,
  10924262.5,
  10924364.0625,
  10924453.125,
  10924465.625,
  10925006.25,
  10925051.5625,
  10925153.125,
  10925348.4375,
  10925429.6875,
  10925596.875,
  10925925.0,
  10926098.4375,
  10926128.125,
  10926151.5625,
  10926454.6875,
  10926654.6875,
  10926770.3125,
  10926889.0625,
  10926915.625,
  10927321.875,
  10927596.875,
  10927726.5625,
  10927981.25,
  10928145.3125,
  10928543.75,
  10928789.0625,
  10929031.25,
  10929031.25,
  10929532.8125,
  10929587.5,
  10930218.75,
  10930218.75,
  10930232.8125,
  10930420.3125,
  10930454.6875,
  10930664.0625,
  10930860.9375,
  10930870.3125,
  10930900.0,
  10930928.125,
  10931104.6875,
  10931218.75,
  10931448.4375,
  10931451.5625,
  10931681.25,
  10931843.75,
  10932384.375,
  10932676.5625,
  10932721.875,
  10932757.8125,
  10932917.1875,
  10934662.5,
  10934687.5,
  10935039.0625,
  10935148.4375,
  10935503.125,
  10936704.6875,
  10937028.125,
  10937239.0625,
  10937546.875,
  10937560.9375,
  10937981.25,
  10938842.1875,
  10939871.875,
  10940214.0625,
  10940606.25,
  10940781.25,
  10940810.9375,
  10940960.9375,
  10941887.5,
  10942423.4375,
  10943606.25,
  10943848.4375,
  10944360.9375,
  10945128.125,
  10945957.8125,
  10946528.125,
  10946817.1875,
  10947159.375,
  10947937.5,
  10949028.125,
  10949226.5625,
  10949809.375,
  10950296.875,
  10950470.3125,
  10950673.4375,
  10951807.8125,
  10952471.875,
  10952725.0,
  10952778.125,
  10952784.375,
  10952873.4375,
  10952900.0,
  10953121.875,
  10954700.0,
  10955804.6875,
  10955865.625,
  10956117.1875,
  10956343.75,
  10956437.5,
  10958182.8125,
  10960417.1875,
  10961435.9375,
  10961564.0625,
  10961642.1875,
  10961832.8125,
  10961914.0625,
  10961915.625,
  10963120.3125,
  10963295.3125,
  10963384.375,
  10963650.0,
  10963768.75,
  10963793.75,
  10964090.625,
  10965314.0625,
  10965384.375,
  10965462.5,
  10966151.5625,
  10967268.75,
  10967918.75,
  10969079.6875,
  10969656.25,
  10969787.5,
  10971337.5,
  10972368.75,
  10972854.6875,
  10973035.9375,
  10973556.25,
  10973951.5625,
  10974559.375,
  10974835.9375,
  10975565.625,
  10975734.375,
  10977785.9375,
  10978518.75,
  10978632.8125,
  10978659.375,
  10978737.5,
  10979903.125,
  10980192.1875,
  10980512.5,
  10981354.6875,
  10982539.0625,
  10982834.375,
  10982890.625,
  10983115.625,
  10983232.8125,
  10983560.9375,
  10983590.625,
  10984081.25,
  10984831.25,
  10984979.6875,
  10985089.0625,
  10985362.5,
  10985415.625,
  10985467.1875,
  10985660.9375,
  10985714.0625,
  10985831.25,
  10986034.375,
  10986057.8125,
  10986114.0625,
  10986234.375,
  10986971.875,
  10987185.9375,
  10987251.5625,
  10987528.125,
  10987803.125,
  10987950.0,
  10988021.875,
  10988181.25,
  10988284.375,
  10988796.875,
  10988881.25,
  10989420.3125,
  10990409.375,
  10990484.375,
  10990981.25,
  10990982.8125,
  10991706.25,
  10992010.9375,
  10992103.125,
  10992479.6875,
  10992546.875,
  10992654.6875,
  10992743.75,
  10992884.375,
  10993004.6875,
  10993120.3125,
  10993168.75,
  10993217.1875,
  10993271.875,
  10993303.125,
  10993303.125,
  10993315.625,
  10993532.8125,
  10994193.75,
  10994250.0,
  10994668.75,
  10994884.375,
  10995017.1875,
  10995057.8125,
  10995295.3125,
  10995485.9375,
  10995782.8125,
  10995939.0625,
  10997456.25,
  10998806.25,
  11000315.625,
  11001193.75,
  11001395.3125,
  11001623.4375,
  11002075.0,
  11002526.5625,
  11002592.1875,
  11002617.1875,
  11002956.25,
  11003267.1875,
  11006328.125,
  11006500.0,
  11006693.75,
  11007612.5,
  11008003.125,
  11008017.1875,
  11008315.625,
  11008356.25,
  11009971.875],
 [22.021009242355184,
  38.4890783953468,
  47.02068742803747,
  28.392128982770146,
  66.74958698675533,
  55.22233102952246,
  5.767541976575615,
  85.98545176225845,
  76.9788880898114,
  25.313246423077054,
  26.001943124255078,
  92.83992457127137,
  22.140125767721898,
  6.775345668204583,
  23.317082582744632,
  29.846308908107645,
  60.91451679256374,
  7.328955250217951,
  10.965439446176706,
  54.96693523938706,
  5.0788987163568,
  49.9400207927214,
  8.494233821597778,
  17.79456726236055,
  26.614098956838923,
  24.077513630222768,
  7.881532955875548,
  51.47549429520111,
  38.09964661220158,
  108.72243320993135,
  48.32018342004417,
  12.352624450128737,
  15.216745508036997,
  47.56984033199902,
  22.11036563746642,
  14.560200383213502,
  13.21185426504061,
  17.00803663695596,
  77.58990132771481,
  18.928524810622605,
  18.26483942357279,
  5.528073577416909,
  27.787439662382923,
  15.147994241430474,
  77.08979221999557,
  44.22948272239949,
  5.82832345829075,
  41.997380502800254,
  121.56549943227482,
  89.5261471073743,
  16.279487112996037,
  18.588269635683083,
  31.178240603922152,
  133.54967171658197,
  46.77596037010173,
  26.653184434819728,
  79.97654101549485,
  65.43011049508247,
  17.955891771497473,
  61.90826001779653,
  20.868631608281472,
  30.946566757831633,
  61.02832664774631,
  165.8110614041866,
  53.76319944199383,
  20.34710226503025,
  26.1027783369538,
  18.50971128519458,
  12.056668303282246,
  17.21388380696328,
  89.09662348260042,
  11.821077989476422,
  200.72604604634267,
  5.4336971123951825,
  16.10064007640934,
  15.351378248519651,
  66.85188253566955,
  94.70977099622877,
  38.94508089319679,
  7.364347189617007,
  56.823631086798514,
  6.9515240150768,
  23.655785158591573,
  109.94492113787433,
  22.57318643421639,
  28.80568689671915,
  36.86961944261854,
  7.283532726832849,
  25.09127887004715,
  46.893205239032824,
  29.074893782244324,
  137.60778568136072,
  45.539803689582556,
  34.73445517965499,
  98.36953253093536,
  5.276151751162502,
  16.73576318457908,
  20.490710855333738,
  11.666800405834781,
  40.666781069542665,
  95.499035707235,
  60.32844883773057,
  53.09497305128083,
  11.138284956021902,
  34.41457539434679,
  18.538386078888358,
  47.81573168259413,
  13.823627868534967,
  88.07231232320218,
  16.841576235918975,
  6.313443538055007,
  46.01188192946974,
  35.4555126037289,
  62.31103758371546,
  18.260874149808338,
  70.50340101649549,
  16.998355474105303,
  32.349993325646366,
  36.35025351582656,
  9.905739734247941,
  6.272346554264856,
  16.548075951133608,
  173.4780343651496,
  14.66564857523588,
  18.176815182647395,
  25.105044093816872,
  16.020130987277902,
  19.486453323053702,
  20.430833529285074,
  27.70935736374711,
  36.63704332839521,
  27.408499881318978,
  19.671181140231766,
  11.382255464408123,
  16.694242426628662,
  72.4963068669911,
  16.224052722840963,
  32.937620370474534,
  13.054374663405506,
  53.593902462809424,
  37.660485575055304,
  124.6842273463259,
  11.973016827946285,
  20.17553459537155,
  24.393674603768687,
  5.041678731527717,
  23.24902946044456,
  76.2990777405718,
  18.977857621516563,
  62.036722606662316,
  7.0193588417710515,
  10.981624497878927,
  8.63548794277709,
  19.284447395602363,
  91.59509239318692,
  5.239379273968614,
  18.311129497436117,
  83.30861089082661,
  68.88100814463934,
  11.500091722865413,
  81.27556251771465,
  38.88334790456944,
  17.95606701292393,
  34.5942353457635,
  11.094001898801896,
  30.218668317589874,
  28.793985019134244,
  35.00334716868556,
  89.27003577678465,
  62.178988330270435,
  5.786125626670979,
  15.4954564536649,
  14.245958827256203,
  12.20793513230758,
  17.964454316839333,
  98.4255103702047,
  17.351265301932973,
  7.86100509754009,
  22.050736113141603,
  26.438877054565015,
  88.08701175903134,
  13.770208485666899,
  88.72686862228527,
  163.30563692716265,
  52.090468440495556,
  21.586448490358404,
  9.474736021000595,
  36.213764981720736,
  85.08699727616452,
  214.53261368846267,
  10.740221575344053,
  35.77077205935591,
  91.18864709509413,
  16.229036704589458,
  10.799802374587344,
  94.71323087772357,
  39.0077658431792,
  5.621626020592023,
  21.458026174122292,
  39.61612410696293,
  23.427202472568087,
  30.099284242926736,
  157.64861930971978,
  30.809233968790622,
  5.8023261510349915,
  13.519787243402202,
  13.73244761866103,
  14.153340818472302,
  17.512550123563887,
  5.799873310174771,
  65.81773487747353,
  31.09619083558779,
  24.440747308417343,
  144.62507467580127,
  65.51660079114662,
  33.451659394009475,
  8.00529218687174,
  59.24554381775309,
  8.748882696933995,
  55.41716628969277,
  5.831643974250857,
  223.5296764355819,
  30.710301430118875,
  8.697083915282567,
  30.55764708558197,
  29.11562981964126,
  7.71821965787399,
  5.332107034424148,
  8.953723696357054,
  12.740229987697337,
  14.068827868213816,
  36.27327213427412,
  32.020303071632846,
  44.60634230266353,
  16.57414431456052,
  6.434576193701606,
  14.907339095605545,
  22.890187732247284,
  5.24721410341351,
  25.26679854568272,
  15.271712538456171,
  36.6513410796253,
  5.383907988892966,
  24.901555215662682,
  25.036365624854696,
  139.8300860633072,
  9.680220528056239,
  6.11022503291402,
  15.076892778596527,
  31.1540587976988,
  22.547681036207095,
  11.572712877285262,
  65.87077213573423,
  37.741516255717784,
  98.86840405537988,
  8.814990455582633,
  41.37817426786875,
  5.204984731342885,
  41.52914443715207,
  106.0944401293016,
  15.021078390882824,
  22.82330875321483,
  16.66503134240758,
  64.0284232252087,
  43.765999967727694,
  56.41239548588396,
  31.73265188069879,
  31.92621215244271,
  39.5202020939635,
  12.743189561458843,
  19.349836969610667,
  17.93846462401773,
  49.29586602947651,
  15.279825888494805,
  34.50729847707833,
  77.14075514074692,
  37.47004128373334,
  10.928065331906891,
  22.9258506804753,
  151.8548261866377,
  50.36576040721054,
  57.160279098744134,
  5.1463006665050255,
  16.246580314279104,
  37.54345238010842,
  36.79520231005443,
  26.596550451061262,
  30.02696603035968,
  30.29493832764156,
  66.75335578822245,
  22.47519437477357,
  20.362495308891447,
  13.782695860583036,
  5.105228635646729,
  11.465829376306463,
  24.689489332707915,
  22.468568125350938,
  18.251894050617672,
  10.41531876261297,
  13.710940199974536,
  52.34867932798088,
  42.540561392271215,
  9.865211756656256,
  19.86159492684419,
  13.716269438775196,
  8.12563030131686,
  6.079465494844198,
  24.220883659054937,
  7.680922032547267,
  30.718974200800645,
  7.282300362708002,
  67.03591453283133,
  44.46697470882711,
  9.25112182635102,
  35.64006940966409,
  7.054792151097475,
  36.9233388624636,
  38.408505668803954,
  31.26842943641544,
  44.97547688189459,
  72.49704185320105,
  18.03348940304471,
  128.7066757451429,
  134.79254782541915,
  9.099341105235155,
  28.009804625469606,
  6.209783846193417,
  65.0208204057012,
  13.855440865426292,
  115.08518137467516,
  37.232730994110085,
  10.031061446947488,
  59.63945759145744,
  76.92252500768866,
  41.165820552689155,
  95.07188602774748,
  24.657813383502127,
  11.137720828219305,
  10.496877095393701,
  25.47161697023405,
  55.64807023954094,
  31.104682557426386,
  13.26363570340616,
  43.488304977597295,
  13.233473167898177,
  15.360079603334214,
  22.848670904501805,
  42.559909705301784,
  34.67602754924622,
  37.97093729504877,
  89.55418966788866,
  23.655749886430065,
  25.256919586640304,
  57.95666255489811,
  26.289583863198995,
  69.03233943083873,
  87.2797337057543,
  6.682150561842103,
  24.833240771015063,
  47.868893005645724,
  31.129884345645216,
  68.4942083232207,
  5.181325306025436,
  37.8727256560143,
  69.26176883543641,
  43.32495886566604,
  43.07275486455848,
  15.867927464700793,
  56.842124237070955,
  11.026938422475894,
  29.742461821466048,
  35.477831521480155,
  109.3976900512591,
  22.06295365220376,
  34.197541967758205,
  16.91099457022679,
  23.679096116763116,
  56.24427229374366,
  13.056530399021842,
  5.245650491351172,
  16.585874732501,
  9.519687281671542,
  5.3110270733573826,
  49.98210644744865,
  10.533285023869338,
  14.03520230293683,
  80.53185550481206,
  142.20907185489466,
  13.621867184182264,
  10.496795945218656,
  52.794487016302504,
  5.81852381302713,
  5.11988686230411,
  13.010602722255353,
  72.03845711836274,
  19.96349716340491,
  7.4510984617700595,
  6.519807299136634,
  76.72752828676161,
  14.955102447042108,
  17.488279106598295,
  13.2240439570794,
  77.17335669046227,
  33.601244916449915,
  43.65337426395233,
  10.535746741680837,
  21.54555258334053,
  20.873519761088723,
  28.238810736462195,
  5.778985931681963,
  7.711172389454429,
  52.85776382538101,
  81.54478801636546,
  9.975727577410039,
  5.156266483097807,
  88.13017486023823,
  17.535101307001124,
  37.6776960318424,
  62.6841803427774,
  14.224764176988328,
  18.411617199313568,
  41.170161100909084,
  15.974723954580107,
  21.634210021222763,
  35.126478197135846,
  29.38502145942212,
  12.113958991867475,
  91.67982697504785,
  38.62457881105571,
  31.782638781353505,
  64.6841965692721,
  11.417150361034581,
  73.87486338456806,
  8.806014972258357,
  18.95679500699656,
  39.14002171153467,
  22.41591445096192,
  65.99529546426692,
  28.67072992745748,
  8.805729919763305,
  65.99059132796356,
  28.11468680341484,
  76.44349714650832,
  66.16629340087987,
  10.152127676087854,
  6.357698827331542,
  84.66973836414083,
  10.592715304784974,
  65.64720746561785,
  47.66557987703953,
  23.33183279362798,
  44.579767139026444,
  16.380984123555233,
  7.905809331005398,
  65.33143082728952,
  38.16439296293021,
  61.836374435680966,
  37.69603815554154,
  22.733831852294873,
  10.911905571870372,
  74.65633083934625,
  12.508210428003846,
  9.725340240774221,
  48.09394253109295,
  97.45586146597684,
  84.18692986328307,
  7.898691548649248,
  10.68649086219551,
  37.46424444986359,
  24.73213961121516,
  41.872491913159266,
  40.62338597392003,
  17.05929120717752,
  12.748262745258055,
  9.773554885743303,
  28.227988516823736,
  5.433260964862621,
  17.68868122704911,
  42.4351797973745,
  48.472047921834985,
  9.856016171300812,
  9.06385286707796,
  18.72618452329767,
  17.19172638298216,
  77.62683033630735,
  5.350408738443339,
  81.43577629676666,
  57.4512005163295,
  65.90236095603541,
  28.558262024427805,
  36.25644384351509,
  114.36925168269248,
  16.345061629878007,
  78.91918042756956,
  5.3961460087826705,
  184.18845408512198,
  5.44601500182603,
  18.178213819118735,
  38.768675227088394,
  57.638076333475645,
  34.11240184579252,
  42.55537570015323,
  14.561482675649726,
  40.89795889076275,
  10.857965566203532,
  87.63891674092557,
  73.51387920871596,
  17.6120719174898,
  5.249086831333948,
  34.32455983600745,
  5.171124406431231,
  9.094198901350126,
  27.623466950636892,
  45.21400875983571,
  28.284419286661436,
  111.75056451168965,
  33.67624455001214,
  6.813122942536255,
  6.112532485850941,
  83.34096537565901,
  97.63111396509836,
  38.36411684364267,
  79.12613931639416,
  19.856801445326067,
  63.479904976588095,
  99.23966045508843,
  5.716553670554265,
  203.4042285228039,
  36.85084191552412,
  13.847658985820932,
  42.40245615647174,
  7.7619145515877115,
  10.125094502875928,
  22.59347082126714,
  35.5995571271542,
  49.83339886896031,
  13.023547661770873,
  86.32743717240116,
  5.365018795285513,
  30.15722298452897,
  53.3130026198923,
  40.11615393095157,
  104.90759889866919,
  22.038872957895677,
  111.06000894916396,
  86.21779869376928,
  16.648411062053647,
  15.122188221791173,
  119.67093414439576,
  56.10826696282055,
  30.4619104345695,
  85.57323752579141,
  26.64526035465237,
  5.512210331251039,
  16.168188756214402,
  50.16706543608371,
  7.862688686332413,
  21.174991404603016,
  29.770891245999692,
  53.83780875022511,
  45.27308263126477,
  18.849458107855693,
  22.682093444403883,
  9.267214533984296,
  84.65023467444472,
  19.710697184538045,
  22.2028086579542,
  29.509347136265312,
  6.933488576430554,
  22.683479401268528,
  5.271866949848501,
  48.45177151309455,
  94.34272152631485,
  50.74603832448135,
  65.8893793756094,
  59.18668021183142,
  6.411250187473843,
  34.330131028621736,
  40.6576972695079,
  11.104617550258228,
  7.887705892111198,
  5.18005252761478,
  19.746549438652668,
  92.05133031035159,
  35.05692637374621,
  45.4773152326746,
  23.73527356299942,
  22.16014330001128,
  7.680815981306755,
  28.249317472032615,
  94.17449764976114,
  12.395297198884663,
  65.55968764439035,
  8.188548098923544,
  19.76390527371317,
  13.782718952804494,
  27.463261022589784,
  7.253873108500993,
  7.098247328973293,
  87.34684372462829,
  66.49412408220236,
  7.283875763722794,
  17.330554785575572,
  36.52366381780125,
  8.01223805647984,
  56.429363548751994,
  54.06117003741237,
  14.770592630022241,
  26.42326911590387,
  59.868183019168306,
  43.23805501909495,
  10.214640248906216,
  59.72851590158683,
  20.659875269530673,
  23.493264175909378,
  31.204148809644746,
  25.198365358576496,
  10.903033004281209,
  51.03699798151449,
  26.935533597535574,
  52.816519902510024,
  28.4479667681646,
  13.573508083331795,
  13.31283556640876,
  59.11371970603487,
  27.41232303341799,
  63.92876644192435,
  16.3301677067414,
  30.17306679133136,
  73.05017821255817,
  43.76804684322735,
  9.377318534420954,
  101.34424601873934,
  41.035859801503136,
  24.265292619418787,
  8.277266186050756,
  34.45577091187768,
  63.05802953976904,
  11.163152759147,
  15.490200850043493,
  7.000659673439083,
  18.65522636836891,
  8.611027181783722,
  15.7071408820859,
  15.671686514707517,
  15.551000979005742,
  55.86499579745119,
  21.053437689644273,
  73.85920695239307,
  21.616976840544154,
  113.0408594774431,
  15.286231190939306,
  24.26820016859054,
  35.35622195433799,
  65.41433752029411,
  59.08459494152116,
  23.314415555225686,
  5.617862029888876,
  12.961658684330763,
  14.75647783470239,
  8.73183091086844,
  29.614177463402136,
  20.342343756181247,
  30.213886107224802,
  5.058864945783083,
  17.524693815422285,
  25.07481566810459,
  9.215964805219468,
  15.434852025089148,
  10.798638114703666,
  24.564514480703515,
  25.390587529801348,
  18.9742667465705,
  18.908219815034375,
  30.34449822091364,
  14.624239766767865,
  46.16688792442261,
  20.307999884402832,
  12.946293858889382,
  48.208410674804554,
  5.062074776884923,
  16.95465162004606,
  16.6640376812401,
  47.45443784378364,
  19.370594172227264,
  29.58239319234457,
  65.03935836284646,
  5.339394518650768,
  5.874130199960952,
  22.421037076783215,
  46.58539084372967,
  7.990716144309684,
  22.979159192473006,
  7.193342301813022,
  56.42877660935584,
  38.63725540578667,
  21.01436858523974,
  14.171886003128202,
  63.64548751945541,
  29.664315790507796,
  74.3046080298979,
  18.35924808820434,
  14.472804136288927,
  60.94876527352899,
  14.602847740925831,
  116.40782153342616,
  8.992212337913951,
  5.9875307768331725,
  78.45442709647688,
  74.62526212553288,
  46.40240576676584,
  12.054653876271365,
  5.285194301886191,
  7.839299582827493,
  63.392201058440435,
  54.77349197330699,
  12.847227483279545,
  71.80080956859267,
  14.586885534611513,
  11.755643620206742,
  16.436451660646064,
  10.46448021601351,
  16.851498424285793,
  8.626307341908115,
  35.91879378318802,
  46.550582331435805,
  8.046991122943412,
  55.17416982013076,
  46.85609323020705,
  123.32736537632918,
  10.134955948997298,
  29.108489257247527,
  21.27870430558702,
  25.34670583857635,
  49.5673741311398,
  71.76938397849558,
  23.111600725871835,
  50.19117509704823,
  8.680596800741602,
  55.06642405893906,
  10.329415999954492,
  90.42979675993865,
  29.71168000827199,
  23.82351004444102,
  35.71374341061582,
  5.510616214255254,
  55.83269388527026,
  5.203303906647725,
  7.079532911234448,
  10.34548884351672,
  23.688498973806844,
  44.93005624181225,
  47.62743080132242,
  19.80910117496984,
  15.33786883264623,
  61.1266382591012,
  43.013515943568244,
  124.00386625107237,
  52.01179770924548,
  81.10372725906613,
  15.249657599687554,
  12.384611721622107,
  17.71189451480888,
  34.728231331536676,
  20.475720516204397,
  65.64446515560198,
  9.566740644794663,
  11.445634567423177,
  68.9017362685597,
  29.98705052870775,
  11.139339005510818,
  13.30437144220927,
  27.846160543068166,
  5.674755432438352,
  12.343400811292362,
  16.94600758369818,
  20.42515646291394,
  29.364399613952926,
  61.024032258373,
  7.105237067849262,
  11.936150474254427,
  20.701798610012098,
  27.923475131570207,
  21.169829402338863,
  10.678265085728777,
  103.51254346744452,
  72.44795854086237,
  10.915504204385975,
  6.776016059219304,
  17.290112236251346,
  14.548129569537364,
  21.262300580920257,
  14.650993960765222,
  19.795831747331114,
  5.928722362066054,
  5.715891480488105,
  20.707832994601986,
  17.599363625216174,
  13.4017068908327,
  14.668583908856826,
  25.27435717097639,
  20.64032799362026,
  81.86736506133546,
  15.225643744904625,
  87.47437020994107,
  29.30794757172552,
  18.102969558140508,
  8.438859460490141,
  56.3849649656277,
  34.13813401389431,
  36.0818578019151,
  86.66110046385121,
  20.708465544551917,
  77.10642400050229,
  44.34881107527772,
  14.713167720706386,
  71.51845277837111,
  22.375394475295273,
  16.28818993360847,
  48.57040585597619,
  41.96772091864355,
  22.283426862356265,
  71.72591491862836,
  20.626442323309714,
  84.91838809504783,
  5.292777944683193,
  13.336394033125114,
  44.48343760206406,
  57.336228013175926,
  5.410824069280008,
  32.36725898468573,
  85.74413940363755,
  13.96525748953117,
  36.63488752594296,
  22.31510641626128,
  32.39794798722611,
  17.56963513315327,
  73.29718922332276,
  28.390550148637203,
  78.737325721106,
  40.13258517613334,
  24.2569644572232,
  5.73677101680262,
  17.007721219121464,
  27.546531131391724,
  17.88566874786028,
  29.04182172084941,
  5.098931286238759,
  5.362676915658964,
  60.05818336903028,
  21.74389677485656,
  37.48173029369633,
  7.21879373844304,
  44.515530418320196,
  142.56983849035865,
  18.138350794524957,
  7.616096618973378,
  14.289807409547906,
  35.20519866100101,
  27.285240800385132,
  14.368589102463558,
  51.62647468520487,
  25.991859649822743,
  9.074751637838316,
  9.046592571905157,
  14.584119126960067,
  5.937328771401298,
  10.045440512891627,
  23.313660511017435,
  17.107995963080118,
  5.818522107576128,
  71.08696593997973,
  25.442802520418088,
  41.45241381870106,
  66.22400949133461,
  146.57638690000948,
  105.88807152064071,
  5.960893475895645,
  6.148202775042389,
  21.192313707992646,
  44.09465620369771,
  25.457030630336206,
  12.425603059532982,
  24.6471178701932,
  35.382632947966556,
  17.186985693617892,
  5.161512798362923,
  11.884989819364487,
  5.0549004826156985,
  6.386668212430329,
  15.697027261721002,
  16.242766529518974,
  10.170805794759305,
  17.497511785217142,
  78.80800858733525,
  88.79806092413611,
  23.144164809664098,
  13.492445974154005,
  50.9021917385543,
  52.65022817104952,
  121.10695939461881,
  44.3338198699621,
  15.139243419728508,
  46.04415661762758,
  7.385940347398921,
  7.797147971087979,
  79.39219995992596,
  21.09274255073534,
  25.550954959178696,
  8.575364572269729,
  69.00936254497469,
  82.56555452336234,
  11.766986380727856,
  67.68651563337193,
  69.35916430244502,
  79.07667980291255,
  25.504299287895186,
  53.09170737629367,
  12.84771979132092,
  64.6783645165732,
  30.29273264034638,
  20.994389289618372,
  101.47276786301427,
  6.200436275231995,
  27.108089089319712,
  22.83534139160889,
  6.223195707648995,
  14.126174705444814,
  9.093728426326399,
  7.587944510358162,
  66.22406423386317,
  71.55801092103763])

Graph for home page

Detected energies during the discharge + Energy spectrum

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