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

([3364696.875,
  3587743.75,
  3917159.375,
  4005256.25,
  4062275.0,
  4279912.5,
  4419003.125,
  5943579.6875,
  5946803.125,
  5976665.625,
  5977243.75,
  6036115.625,
  6060587.5,
  6162345.3125,
  6166395.3125,
  6167187.5,
  6167535.9375,
  6169118.75,
  6173065.625,
  6173517.1875,
  6174679.6875,
  6189114.0625,
  6189367.1875,
  6194743.75,
  6200681.25,
  6207485.9375,
  6229645.3125,
  6243259.375,
  6279590.625,
  6304725.0,
  6326243.75,
  6336606.25,
  6359184.375,
  6359246.875,
  6362032.8125,
  6364492.1875,
  6364865.625,
  6364917.1875,
  6367781.25,
  6368956.25,
  6370382.8125,
  6371445.3125,
  6375167.1875,
  6375825.0,
  6390464.0625,
  6402471.875,
  6406246.875,
  6408206.25,
  6428725.0,
  6437965.625,
  6438481.25,
  6438542.1875,
  6440595.3125,
  6445176.5625,
  6446743.75,
  6447146.875,
  6450457.8125,
  6459673.4375,
  6464581.25,
  6465748.4375,
  6466065.625,
  6466120.3125,
  6467882.8125,
  6468570.3125,
  6468723.4375,
  6469321.875,
  6470251.5625,
  6470309.375,
  6471634.375,
  6471650.0,
  6472171.875,
  6472173.4375,
  6473368.75,
  6473637.5,
  6474581.25,
  6474662.5,
  6474790.625,
  6475096.875,
  6475337.5,
  6475690.625,
  6475828.125,
  6485023.4375,
  6487828.125,
  6487896.875,
  6492175.0,
  6497145.3125,
  6498739.0625,
  6501164.0625,
  6501520.3125,
  6501609.375,
  6502396.875,
  6502448.4375,
  6503942.1875,
  6505600.0,
  6506640.625,
  6509953.125,
  6510768.75,
  6511710.9375,
  6512196.875,
  6512881.25,
  6513687.5,
  6513792.1875,
  6515100.0,
  6515623.4375,
  6515870.3125,
  6515882.8125,
  6515904.6875,
  6516553.125,
  6516596.875,
  6516915.625,
  6516987.5,
  6517542.1875,
  6517921.875,
  6518796.875,
  6518921.875,
  6519176.5625,
  6520000.0,
  6521760.9375,
  6523034.375,
  6523367.1875,
  6523532.8125,
  6523548.4375,
  6525879.6875,
  6525935.9375,
  6532328.125,
  6540443.75,
  6544451.5625,
  6546923.4375,
  6547420.3125,
  6549707.8125,
  6553698.4375,
  6557842.1875,
  6559315.625,
  6559487.5,
  6560742.1875,
  6560928.125,
  6561950.0,
  6561973.4375,
  6562503.125,
  6564526.5625,
  6564675.0,
  6564754.6875,
  6565093.75,
  6565120.3125,
  6565210.9375,
  6565664.0625,
  6565890.625,
  6566762.5,
  6567170.3125,
  6568292.1875,
  6568648.4375,
  6568820.3125,
  6571121.875,
  6571251.5625,
  6571876.5625,
  6572320.3125,
  6573785.9375,
  6574362.5,
  6574495.3125,
  6574496.875,
  6577071.875,
  6577881.25,
  6581834.375,
  6582689.0625,
  6583628.125,
  6587057.8125,
  6590504.6875,
  6592665.625,
  6592682.8125,
  6593306.25,
  6593353.125,
  6593975.0,
  6608518.75,
  6610675.0,
  6610989.0625,
  6611581.25,
  6613792.1875,
  6613826.5625,
  6614840.625,
  6619557.8125,
  6622720.3125,
  6626937.5,
  6628729.6875,
  6631796.875,
  6637106.25,
  6637487.5,
  6639490.625,
  6642393.75,
  6645617.1875,
  6648360.9375,
  6648451.5625,
  6649357.8125,
  6651246.875,
  6653803.125,
  6654814.0625,
  6655835.9375,
  6660771.875,
  6660776.5625,
  6660985.9375,
  6661476.5625,
  6662220.3125,
  6662435.9375,
  6662737.5,
  6663009.375,
  6663775.0,
  6664293.75,
  6664698.4375,
  6664707.8125,
  6664712.5,
  6664732.8125,
  6665246.875,
  6665481.25,
  6665907.8125,
  6666179.6875,
  6667107.8125,
  6667243.75,
  6667243.75,
  6667531.25,
  6667659.375,
  6667659.375,
  6667675.0,
  6667693.75,
  6667948.4375,
  6669767.1875,
  6669856.25,
  6669909.375,
  6671204.6875,
  6671242.1875,
  6672765.625,
  6672789.0625,
  6673117.1875,
  6673551.5625,
  6673593.75,
  6683201.5625,
  6683554.6875,
  6684132.8125,
  6684712.5,
  6684729.6875,
  6685514.0625,
  6686403.125,
  6686450.0,
  6686914.0625,
  6686942.1875,
  6686967.1875,
  6687146.875,
  6687534.375,
  6687709.375,
  6687842.1875,
  6687906.25,
  6687987.5,
  6688120.3125,
  6688150.0,
  6688173.4375,
  6688479.6875,
  6688576.5625,
  6688626.5625,
  6688725.0,
  6688809.375,
  6688832.8125,
  6688839.0625,
  6688839.0625,
  6688907.8125,
  6689046.875,
  6689062.5,
  6689137.5,
  6689223.4375,
  6689251.5625,
  6689256.25,
  6689273.4375,
  6689421.875,
  6689453.125,
  6689568.75,
  6689720.3125,
  6689850.0,
  6689853.125,
  6689895.3125,
  6689925.0,
  6689937.5,
  6690018.75,
  6690087.5,
  6690129.6875,
  6690229.6875,
  6690407.8125,
  6690440.625,
  6690625.0,
  6690628.125,
  6690648.4375,
  6690673.4375,
  6690689.0625,
  6690709.375,
  6690731.25,
  6690850.0,
  6690934.375,
  6690998.4375,
  6691109.375,
  6691165.625,
  6691282.8125,
  6691296.875,
  6691367.1875,
  6691382.8125,
  6691387.5,
  6691398.4375,
  6691412.5,
  6691457.8125,
  6691578.125,
  6691617.1875,
  6691618.75,
  6691746.875,
  6691756.25,
  6691771.875,
  6691870.3125,
  6691995.3125,
  6692050.0,
  6692085.9375,
  6692139.0625,
  6692145.3125,
  6692290.625,
  6692312.5,
  6692314.0625,
  6692348.4375,
  6692354.6875,
  6692365.625,
  6692465.625,
  6692543.75,
  6692548.4375,
  6692664.0625,
  6692703.125,
  6692753.125,
  6693200.0,
  6693221.875,
  6693418.75,
  6693481.25,
  6693589.0625,
  6693678.125,
  6693714.0625,
  6693890.625,
  6693945.3125,
  6694150.0,
  6694182.8125,
  6694321.875,
  6694412.5,
  6694496.875,
  6694518.75,
  6694560.9375,
  6694573.4375,
  6694604.6875,
  6694606.25,
  6694628.125,
  6694867.1875,
  6694929.6875,
  6694965.625,
  6695107.8125,
  6695304.6875,
  6695348.4375,
  6695868.75,
  6695875.0,
  6695934.375,
  6696343.75,
  6696643.75,
  6698365.625,
  6698726.5625,
  6699153.125,
  6699443.75,
  6700907.8125,
  6701806.25,
  6703776.5625,
  6704851.5625,
  6704968.75,
  6705842.1875,
  6705932.8125,
  6706200.0,
  6706468.75,
  6706506.25,
  6706637.5,
  6706678.125,
  6706687.5,
  6707142.1875,
  6708037.5,
  6708206.25,
  6713785.9375,
  6715167.1875,
  6715550.0,
  6716701.5625,
  6717729.6875,
  6719778.125,
  6719892.1875,
  6720403.125,
  6720643.75,
  6721895.3125,
  6722748.4375,
  6722806.25,
  6722854.6875,
  6722859.375,
  6723264.0625,
  6723445.3125,
  6723500.0,
  6723537.5,
  6723787.5,
  6724260.9375,
  6724289.0625,
  6725371.875,
  6725581.25,
  6725668.75,
  6725821.875,
  6725842.1875,
  6725887.5,
  6725910.9375,
  6726640.625,
  6726971.875,
  6727012.5,
  6727056.25,
  6727259.375,
  6727626.5625,
  6728512.5,
  6728743.75,
  6728762.5,
  6731779.6875,
  6732260.9375,
  6732432.8125,
  6732892.1875,
  6733714.0625,
  6734140.625,
  6734910.9375,
  6735528.125,
  6738829.6875,
  6738943.75,
  6739029.6875,
  6739410.9375,
  6739643.75,
  6740900.0,
  6742184.375,
  6742465.625,
  6743915.625,
  6744950.0,
  6745137.5,
  6746487.5,
  6746514.0625,
  6747059.375,
  6750815.625,
  6752242.1875,
  6752393.75,
  6752410.9375,
  6752435.9375,
  6753500.0,
  6754150.0,
  6758209.375,
  6759045.3125,
  6759801.5625,
  6760979.6875,
  6761115.625,
  6775076.5625,
  6775379.6875,
  6779660.9375,
  6782918.75,
  6788901.5625,
  6790290.625,
  6795151.5625,
  6805556.25,
  6807815.625,
  6822965.625,
  6823621.875,
  6823632.8125,
  6828728.125,
  6829656.25,
  6841743.75,
  6841768.75,
  6845446.875,
  6845564.0625,
  6846368.75,
  6846379.6875,
  6846431.25,
  6846573.4375,
  6846617.1875,
  6846870.3125,
  6847267.1875,
  6849175.0,
  6849739.0625,
  6850175.0,
  6850481.25,
  6852170.3125,
  6852385.9375,
  6868739.0625,
  6868917.1875,
  6882912.5,
  6888640.625,
  6889214.0625,
  6889606.25,
  6889843.75,
  6889950.0,
  6890810.9375,
  6891292.1875,
  6891485.9375,
  6892315.625,
  6892803.125,
  6892867.1875,
  6893407.8125,
  6895385.9375,
  6895560.9375,
  6896421.875,
  6897648.4375,
  6897671.875,
  6897778.125,
  6897896.875,
  6898139.0625,
  6898165.625,
  6898393.75,
  6898721.875,
  6898735.9375,
  6898770.3125,
  6898815.625,
  6898945.3125,
  6899265.625,
  6899350.0,
  6899509.375,
  6899554.6875,
  6899812.5,
  6900217.1875,
  6900429.6875,
  6900546.875,
  6902576.5625,
  6902656.25,
  6906939.0625,
  6913245.3125,
  6913442.1875,
  6913759.375,
  6913784.375,
  6914006.25,
  6923751.5625,
  6925965.625,
  6927595.3125,
  6932939.0625,
  6938412.5,
  6938500.0,
  6943160.9375,
  6945146.875,
  6950165.625,
  6951968.75,
  6957435.9375,
  6960814.0625,
  6964254.6875,
  6964367.1875,
  6967328.125,
  6967510.9375,
  6969062.5,
  6969065.625,
  6970078.125,
  6970951.5625,
  6970954.6875,
  6972007.8125,
  6972178.125,
  6972671.875,
  6973873.4375,
  6975665.625,
  6975689.0625,
  6975725.0,
  6976389.0625,
  6976860.9375,
  6977046.875,
  6977210.9375,
  6977351.5625,
  6978246.875,
  6978254.6875,
  6978870.3125,
  6979339.0625,
  6980842.1875,
  6980962.5,
  6983843.75,
  6984129.6875,
  6984187.5,
  6985410.9375,
  6985442.1875,
  6985471.875,
  6985501.5625,
  6985843.75,
  6986506.25,
  6987567.1875,
  6988475.0,
  6988809.375,
  6988904.6875,
  6991437.5,
  6991546.875,
  6992840.625,
  6993178.125,
  6993340.625,
  6994176.5625,
  6994184.375,
  6994189.0625,
  6994957.8125,
  6995317.1875,
  6995326.5625,
  6995410.9375,
  6995951.5625,
  6995968.75,
  6998239.0625,
  6998546.875,
  7000923.4375,
  7001118.75,
  7001329.6875,
  7002537.5,
  7002551.5625,
  7003082.8125,
  7003625.0,
  7003793.75,
  7003798.4375,
  7003840.625,
  7003860.9375,
  7004200.0,
  7005406.25,
  7007182.8125,
  7011151.5625,
  7011412.5,
  7011417.1875,
  7011759.375,
  7013165.625,
  7014240.625,
  7015301.5625,
  7018817.1875,
  7019067.1875,
  7021000.0,
  7022381.25,
  7022514.0625,
  7022531.25,
  7022943.75,
  7023767.1875,
  7024992.1875,
  7025025.0,
  7025112.5,
  7025460.9375,
  7025710.9375,
  7025939.0625,
  7025998.4375,
  7026079.6875,
  7026495.3125,
  7027435.9375,
  7027623.4375,
  7028498.4375,
  7028890.625,
  7029484.375,
  7030900.0,
  7031848.4375,
  7032140.625,
  7032292.1875,
  7032504.6875,
  7033356.25,
  7034859.375,
  7034962.5,
  7034968.75,
  7035606.25,
  7035607.8125,
  7036428.125,
  7037121.875,
  7037123.4375,
  7037503.125,
  7038662.5,
  7038715.625,
  7039207.8125,
  7039254.6875,
  7039350.0,
  7039354.6875,
  7039662.5,
  7039820.3125,
  7039956.25,
  7040417.1875,
  7040514.0625,
  7040625.0,
  7040729.6875,
  7041090.625,
  7041131.25,
  7041134.375,
  7041375.0,
  7041817.1875,
  7041831.25,
  7041843.75,
  7041859.375,
  7041875.0,
  7041893.75,
  7041990.625,
  7042067.1875,
  7042160.9375,
  7042171.875,
  7042221.875,
  7042579.6875,
  7042771.875,
  7042873.4375,
  7043117.1875,
  7043200.0,
  7043381.25,
  7043401.5625,
  7043826.5625,
  7044237.5,
  7044504.6875,
  7057339.0625,
  7058273.4375,
  7060925.0,
  7061534.375,
  7062500.0,
  7063012.5,
  7063031.25,
  7063771.875,
  7065518.75,
  7065803.125,
  7065828.125,
  7066556.25,
  7067959.375,
  7068065.625,
  7068767.1875,
  7070281.25,
  7074481.25,
  7074498.4375,
  7074615.625,
  7074687.5,
  7075543.75,
  7075842.1875,
  7075893.75,
  7076270.3125,
  7076410.9375,
  7076723.4375,
  7078187.5,
  7079053.125,
  7079450.0,
  7080343.75,
  7085562.5,
  7087193.75,
  7087200.0,
  7087267.1875,
  7087456.25,
  7091034.375,
  7097928.125,
  7098204.6875,
  7098815.625,
  7101859.375,
  7103120.3125,
  7108912.5,
  7108915.625,
  7110834.375,
  7117862.5,
  7121410.9375,
  7122609.375,
  7123043.75,
  7123385.9375,
  7123467.1875,
  7123743.75,
  7124987.5,
  7126746.875,
  7127312.5,
  7128600.0,
  7129276.5625,
  7129406.25,
  7130171.875,
  7130509.375,
  7131173.4375,
  7132093.75,
  7132645.3125,
  7134546.875,
  7137681.25,
  7138637.5,
  7142192.1875,
  7147798.4375,
  7148698.4375,
  7148828.125,
  7149118.75,
  7149551.5625,
  7151535.9375,
  7151878.125,
  7152326.5625,
  7152554.6875,
  7153806.25,
  7153823.4375,
  7155406.25,
  7156182.8125,
  7156943.75,
  7157062.5,
  7157703.125,
  7158659.375,
  7159415.625,
  7162153.125,
  7163025.0,
  7163293.75,
  7165400.0,
  7165432.8125,
  7165526.5625,
  7168846.875,
  7172923.4375,
  7174017.1875,
  7176732.8125,
  7178112.5,
  7181457.8125,
  7182885.9375,
  7184625.0,
  7186610.9375,
  7186668.75,
  7188978.125,
  7189812.5,
  7191117.1875,
  7193868.75,
  7195989.0625,
  7199860.9375,
  7201225.0,
  7202990.625,
  7204723.4375,
  7208831.25,
  7209640.625,
  7209660.9375,
  7210468.75,
  7210945.3125,
  7211884.375,
  7212140.625,
  7212409.375,
  7212729.6875,
  7213804.6875,
  7214389.0625,
  7214554.6875,
  7214753.125,
  7214789.0625,
  7214909.375,
  7215381.25,
  7215393.75,
  7215409.375,
  7215467.1875,
  7215560.9375,
  7215643.75,
  7215696.875,
  7215767.1875,
  7215785.9375,
  7215818.75,
  7215867.1875,
  7215906.25,
  7215917.1875,
  7215950.0,
  7216175.0,
  7216357.8125,
  7216450.0,
  7216520.3125,
  7216676.5625,
  7216940.625,
  7216970.3125,
  7217082.8125,
  7217226.5625,
  7217642.1875,
  7217714.0625,
  7218464.0625,
  7219112.5,
  7219385.9375,
  7221445.3125,
  7222723.4375,
  7231503.125,
  7243095.3125,
  7245731.25,
  7245950.0,
  7246095.3125,
  7248907.8125,
  7250520.3125,
  7258990.625,
  7260260.9375,
  7260695.3125,
  7261003.125,
  7261865.625,
  7266873.4375,
  7282703.125,
  7283145.3125,
  7284800.0,
  7286757.8125,
  7286831.25,
  7291665.625,
  7294006.25,
  7309860.9375,
  7312632.8125,
  7313651.5625,
  7314429.6875,
  7315393.75,
  7320562.5,
  7322117.1875,
  7328507.8125,
  7354607.8125,
  7354776.5625,
  7355726.5625,
  7356200.0,
  7356939.0625,
  7357431.25,
  7394865.625,
  7434064.0625,
  7576507.8125,
  7585914.0625,
  7620604.6875,
  7620907.8125,
  7624446.875,
  7653101.5625,
  7683304.6875,
  7684326.5625,
  7684378.125,
  7714845.3125,
  7736729.6875,
  7740359.375,
  7752598.4375,
  7753520.3125,
  7753707.8125,
  7778196.875,
  7781278.125,
  7781287.5,
  7781371.875,
  7879235.9375,
  8104754.6875,
  8173246.875,
  8539956.25,
  8551959.375,
  8552020.3125,
  8721107.8125,
  8761750.0,
  8872392.1875,
  8907807.8125,
  8908339.0625,
  8910898.4375,
  8916675.0,
  8920428.125,
  8931807.8125,
  9182890.625,
  9193826.5625,
  9195637.5,
  9261812.5,
  9275353.125,
  9286118.75,
  9290906.25,
  9294509.375,
  9318065.625,
  9352184.375,
  9356853.125,
  9359085.9375,
  9368718.75,
  9376306.25,
  9378681.25,
  9411829.6875,
  9454026.5625,
  9456243.75,
  9456651.5625,
  9456993.75,
  9457723.4375,
  9461206.25,
  9462717.1875,
  9463818.75,
  9474403.125,
  9482407.8125,
  9486654.6875,
  9487507.8125,
  9490128.125,
  9490478.125,
  9506351.5625,
  9507842.1875,
  9524360.9375,
  9554387.5,
  9581104.6875,
  9622239.0625,
  9628129.6875,
  9642587.5,
  9649184.375,
  9666600.0,
  9669510.9375,
  9669579.6875,
  9671387.5,
  9673656.25,
  9688442.1875,
  9711200.0,
  9711279.6875,
  9827278.125,
  9841485.9375,
  9889659.375,
  9893678.125,
  9894756.25,
  9894798.4375,
  9895709.375,
  9895735.9375,
  9897373.4375,
  9897423.4375,
  9936626.5625,
  9962996.875,
  9995064.0625,
  10013256.25,
  10026162.5,
  10026946.875,
  10029226.5625,
  10036209.375,
  10036742.1875,
  10042568.75,
  10042981.25,
  10046739.0625,
  10046743.75,
  10049373.4375,
  10087198.4375,
  10090228.125,
  10116943.75,
  10263607.8125,
  11257718.75,
  11267890.625,
  11267957.8125,
  11280456.25],
 [36.39723457256683,
  61.42011606542019,
  57.63331380236617,
  8.929290614776288,
  5.867756950569329,
  103.15167346425969,
  13.632253326347488,
  61.988275601110445,
  15.782696782507756,
  5.227109715887212,
  43.40145628763871,
  102.45907594437872,
  70.63837596907882,
  128.59382627377136,
  17.587650568578034,
  15.793167642101507,
  7.671813627472112,
  9.006379654396119,
  28.70609851781273,
  100.29981138667708,
  78.80329657045203,
  23.225821689737273,
  88.18321797392157,
  38.34357814715931,
  23.47776991851982,
  5.225485454867674,
  5.659364363078715,
  76.76677326861977,
  8.024261796596173,
  24.275962905745054,
  25.377084321643334,
  83.89258329860951,
  23.946462742652194,
  50.93071164949142,
  58.51405958943931,
  67.06249961993483,
  55.847138471574354,
  14.149140945397248,
  20.544470133004168,
  10.0180757997359,
  68.33173525422447,
  108.77033466229885,
  26.194573928714608,
  8.360969807334461,
  11.591584576740122,
  15.758665732968948,
  34.030938217134114,
  5.481389065289376,
  24.376277770638143,
  42.175379433676426,
  53.61122422065942,
  5.653454013029906,
  15.857975173415213,
  65.9539246648907,
  8.310939051713781,
  43.64728620668559,
  192.57406810768012,
  50.15962257049653,
  34.27502430628671,
  85.83646132075555,
  34.622424313148166,
  131.8763932836117,
  70.45556578208915,
  24.117761348660785,
  19.78736390407145,
  19.24625273551636,
  259.5303046208734,
  9.267275539627011,
  94.52168757313635,
  25.67429902289775,
  17.247923347038633,
  16.99510366372298,
  18.338551681601587,
  13.32147694962425,
  63.72091543156779,
  56.37430669301801,
  50.916425410211986,
  125.0897587093759,
  58.48166225709635,
  16.703579595561827,
  19.11011430135429,
  25.51252779403324,
  32.804139755043394,
  10.697281014364044,
  54.34203430499464,
  19.291934742150524,
  124.8974664064839,
  15.609415701052098,
  5.432567156193797,
  73.14207159300558,
  19.598637027356833,
  7.306629199403174,
  11.296177152251182,
  47.552568072581444,
  36.96907227087646,
  54.25984614013294,
  23.25976152933508,
  45.334650100167025,
  176.35599253956264,
  17.607075955173745,
  82.81061103496214,
  56.31303715647137,
  8.231195933068582,
  93.65401967694291,
  157.1143954299284,
  62.03569970000155,
  7.316580896670568,
  7.052124010949603,
  35.7759027797182,
  30.82360542216003,
  9.798540488645234,
  100.03988013196754,
  11.38314142977542,
  7.258916934182001,
  39.52951734617878,
  52.82393617615527,
  34.045747539411394,
  63.760109539502466,
  12.6494988393673,
  29.320364298043064,
  21.823657080944493,
  7.000983790584903,
  26.192950922073393,
  5.0758953453095605,
  14.418031387502166,
  33.6527390632858,
  76.86622242713834,
  13.311024654878873,
  29.35449753039608,
  65.95837123281866,
  12.413601516380288,
  29.296784189244242,
  97.81449871045096,
  99.83564237066203,
  40.2476894423131,
  51.4244131011592,
  79.48771264409444,
  43.938712458056315,
  13.011741221320914,
  12.867665300025267,
  44.411861212108604,
  5.698327666029522,
  23.220022487333814,
  7.677072247957201,
  123.11936950859335,
  77.70092946670269,
  51.24469015790008,
  206.5094388591854,
  22.04565298695637,
  48.40618971229868,
  29.481152250632725,
  75.84270141599688,
  107.3155784615033,
  20.330654925873958,
  101.113585193938,
  68.58657363330593,
  19.050999659011822,
  15.296627976197115,
  32.56451285716979,
  21.81137862581718,
  59.30814293658601,
  60.708985931199436,
  137.42236923465742,
  129.95393069598427,
  49.771180646118474,
  49.45136136964324,
  24.937731699140002,
  118.24452935403079,
  20.073764992819875,
  153.86990491667837,
  107.63172441621319,
  15.016307833826497,
  8.408690366765153,
  30.312961874645012,
  5.997958863078135,
  8.452608138764365,
  5.5570975373728,
  8.001892567757281,
  73.79057910086591,
  21.111881574438257,
  6.066153738741113,
  47.08419560585291,
  56.67510204822391,
  43.71788249107159,
  82.39327886377018,
  57.65663076026064,
  13.303113273263156,
  10.330633974957554,
  5.755159609170377,
  19.040784110635492,
  5.621669378083233,
  69.64824848732714,
  33.21099769977836,
  62.609628259637134,
  69.44656875655085,
  19.917878241918054,
  73.48391273876348,
  15.59866757654158,
  28.525795935709684,
  16.450191714611123,
  19.085978566616227,
  9.847199642537985,
  13.939621181470327,
  17.2774418111779,
  10.45611498110723,
  8.539391555424086,
  15.47772590263954,
  51.50271746735625,
  27.51102222060111,
  39.995722741166865,
  10.561223311222646,
  5.49784137905199,
  32.066418116179484,
  96.33037548367857,
  39.25960964222006,
  17.682133427314476,
  96.80524359126099,
  111.17831515119403,
  11.504910229982594,
  16.707816525480926,
  71.48760686188146,
  27.18776603881378,
  38.39960078758963,
  16.4254365726055,
  60.08746596578121,
  45.83385414665403,
  23.353312654207613,
  58.85832501460436,
  41.66846193034802,
  25.889133080981267,
  27.0149215014014,
  112.80447085716244,
  10.512909891136907,
  13.662628979492569,
  70.4005236572396,
  74.16656977228209,
  35.238825696528004,
  28.48184039740719,
  10.924899051597857,
  83.17739913908383,
  16.62831994178219,
  60.086958956855675,
  57.91383276906872,
  20.148785539171563,
  18.427659911220164,
  65.26106321836986,
  7.572559540440718,
  115.73928889501224,
  31.863727267981485,
  82.51041484203607,
  35.361293566290534,
  23.648882747165857,
  56.37135508295516,
  101.0228735666784,
  36.995774727011955,
  239.95573348165058,
  8.574192640784469,
  37.96156166995581,
  89.42040732347566,
  20.22137316037128,
  49.65146830305179,
  20.658802593847145,
  27.283433359698986,
  109.55191884710541,
  23.196824737527788,
  146.457331866237,
  10.2022433282344,
  37.952321428908974,
  13.13483637306454,
  9.545625899412434,
  13.385410919813104,
  38.842488194322364,
  17.783996143897205,
  14.750502476257601,
  6.112562033292657,
  18.771899827565584,
  24.80375480532396,
  78.04036275641678,
  46.72856783120531,
  5.295269109919355,
  93.29325981912626,
  99.98597977415547,
  5.966135850601805,
  14.261884710409893,
  16.048538295470127,
  74.57563138024784,
  5.3234613789131,
  34.57341420045836,
  5.3586324095810705,
  44.485182068559666,
  45.081714863343834,
  16.875670462561438,
  88.20729215095666,
  12.437618944272808,
  19.696139110195446,
  71.00349229313728,
  23.773735673476274,
  105.41442588999934,
  28.946077332056817,
  47.253734933968076,
  5.462575582795934,
  30.145835996026317,
  40.09501662518248,
  11.427284583854847,
  5.55686999448816,
  12.318646529880045,
  17.88578485995733,
  107.07422046782936,
  18.249795060960555,
  32.26400018416813,
  7.930158419839312,
  110.67809866341692,
  100.21268052669818,
  42.514685206195864,
  5.095016070893412,
  61.636161580364984,
  70.08682666651234,
  11.157789722497984,
  29.222816861491555,
  12.971300139419613,
  27.543324638141456,
  7.573531274032749,
  12.324898001981433,
  29.7710779204007,
  9.14804439315284,
  18.360638645955152,
  112.90374119668057,
  132.67254040535673,
  7.162134145727406,
  36.49372334823385,
  6.622016997883784,
  80.68135832342915,
  13.82305552678966,
  77.48743691819745,
  8.92290458525,
  14.580672675474293,
  45.27371879263896,
  7.602347372635561,
  46.92059890289319,
  23.19730413326542,
  7.855230077894405,
  15.048063946899774,
  32.50167566475677,
  16.400917327583734,
  18.337144223995534,
  40.4265125267047,
  6.59677491956983,
  39.655220661190405,
  40.04895538926286,
  56.05252352265611,
  15.331484932504182,
  13.099252738344084,
  22.280172788529264,
  70.29092211381443,
  47.92602721364865,
  195.334377726145,
  5.323626736826369,
  20.677837781591993,
  88.1626383636549,
  10.108983179239246,
  28.152152489077626,
  19.25695419449508,
  54.61985711274021,
  26.572763025421242,
  6.629374543258719,
  25.05963433413907,
  41.21229926008208,
  77.93295679007664,
  250.77590879265006,
  9.035192382103938,
  54.87175112220168,
  66.64111531820397,
  25.91709629781178,
  76.67381992941066,
  53.40916499705752,
  16.340801708160313,
  26.298788136451915,
  214.9276806945125,
  66.19947733960495,
  68.80751329046113,
  63.46934616769979,
  13.537669646613109,
  89.7327046603615,
  102.63149776891655,
  46.08703091630864,
  37.1145750733646,
  21.071852348393904,
  42.392824316576274,
  17.86161228472567,
  42.425304714793135,
  45.48455474609187,
  50.216542413045,
  5.360485862196153,
  19.175318322162177,
  11.012932792400386,
  72.29662011971712,
  135.01674338159373,
  61.074610629654856,
  60.62720243901956,
  15.357241098856736,
  7.759515874011704,
  5.989970420166393,
  5.024611979714108,
  9.060548530053877,
  38.32735034697739,
  62.119420008907156,
  24.17980582768046,
  13.694990792726543,
  57.4125608820768,
  9.893305914485136,
  81.69622814528651,
  21.74521986225286,
  8.750634094959812,
  29.14782057690727,
  63.97159745895686,
  7.130111006729035,
  59.81747671424412,
  12.539494478704636,
  7.291013191795108,
  26.3317224034226,
  72.43517609089575,
  116.95774507794263,
  91.63969363108122,
  15.227370288491633,
  35.72033358393782,
  27.30521793086246,
  5.210992339876186,
  22.58353781299215,
  64.3951221857838,
  13.805659203132972,
  22.311537867308818,
  8.16569461317244,
  107.99538676849032,
  17.077027999769637,
  47.66027591322897,
  41.37640024765623,
  42.07896026079423,
  39.4994745302631,
  69.97118573483006,
  55.399735251854224,
  15.92276263285325,
  17.585232347653264,
  51.09300150269347,
  141.83538047771194,
  25.860473614439737,
  51.70004967125205,
  138.24599936323148,
  6.296551660553312,
  27.066458099730465,
  66.53279747520995,
  81.67574145991111,
  115.38709229549974,
  18.80387966381403,
  43.13626442991365,
  117.90100940387464,
  193.27713109874847,
  106.92378733288527,
  28.925922881684542,
  42.98317903938009,
  104.82569972617829,
  18.09353498595276,
  32.274513881869055,
  45.67891945659834,
  94.57165448160939,
  47.54650021282325,
  14.895962407906907,
  26.768832498528994,
  17.551074736142564,
  30.699116105173232,
  6.935831046905322,
  115.22257517935395,
  7.469597110313175,
  87.76036775135202,
  12.25074032922134,
  41.11555506589667,
  23.83262040248162,
  18.002465826479334,
  23.81477842762854,
  8.838077822794379,
  21.68258229385632,
  28.00852373244167,
  90.11889493377922,
  163.07121689626797,
  20.540108270706625,
  65.80541814133274,
  47.786993726294426,
  29.335042397402365,
  93.86532569676896,
  99.56934917062974,
  17.757248337781288,
  55.39859283739911,
  66.25843115182607,
  118.28008005675935,
  30.671089556711923,
  5.051039471698173,
  70.7305224709122,
  119.90700956452221,
  230.26580581163404,
  25.971573098422343,
  12.149888203213914,
  5.703698023666907,
  152.02059815483406,
  66.6988898443032,
  52.15059682527077,
  28.950333905865026,
  7.908227402268309,
  184.35839680302357,
  117.90609263310223,
  46.24631401665481,
  14.498498191791379,
  54.689428424699976,
  20.35447639425396,
  61.38881776133482,
  38.96333242894729,
  25.431385167612866,
  84.12197391091124,
  12.184693790405046,
  6.327233705329153,
  31.548206682192713,
  15.271782150545015,
  5.238278021608867,
  57.11939878703496,
  97.18729218157577,
  57.077332620175994,
  15.686016685863741,
  17.702169318602504,
  66.30863660574207,
  27.550598380821416,
  163.15543480322236,
  150.88784637190432,
  30.871435262025717,
  91.10225921377655,
  25.816027637756207,
  78.26760028282013,
  26.755804301441145,
  12.876748788353112,
  25.8584221967763,
  24.449481349288803,
  30.709596758514653,
  231.14298469256718,
  18.970093123992815,
  48.1997996113733,
  25.00684685505102,
  29.385623027793265,
  44.436833822308195,
  14.250795298761048,
  207.57083810619093,
  130.15807504395414,
  119.25182099085103,
  24.66519540315926,
  101.99071395184572,
  12.308830987562384,
  57.82716184538903,
  9.731976670841103,
  129.23115679400652,
  16.72269343241506,
  5.423619487829046,
  70.55877603142835,
  19.425067156256375,
  25.788211884679193,
  21.750981785451568,
  172.121249594233,
  68.34847770250651,
  56.21839414178929,
  25.765890989278137,
  78.6671097282135,
  19.638280989855588,
  23.331589709919193,
  30.161396048893327,
  7.708445399624251,
  46.14224247816008,
  5.540215973232311,
  40.7158900641637,
  41.11487049581796,
  40.21710777591162,
  134.14653994173736,
  17.504486033334462,
  52.23380160910488,
  80.08579439489925,
  5.922549508149697,
  19.69870709776275,
  10.377037135534696,
  24.217721677720853,
  70.85770556994306,
  63.97574868434822,
  53.62949064490216,
  15.11559126517665,
  51.747009104773994,
  81.054476769162,
  43.86776692651202,
  202.6870560424551,
  18.108006834507794,
  7.7722192158968655,
  6.105094811494561,
  5.373087842024246,
  14.050546904503438,
  32.859317000322775,
  18.099902704058707,
  121.00330014735792,
  45.78391145108671,
  99.5478504324995,
  102.88623227747873,
  116.5460079068383,
  27.886941051564182,
  66.91752860174782,
  5.080617267890299,
  108.49749951238559,
  87.86397757255472,
  7.114315597064646,
  21.569724714252455,
  23.42700522992944,
  26.796810870875007,
  25.071296640276515,
  132.63278955868236,
  8.412439098701245,
  28.48464565877721,
  36.88131302685445,
  67.218326194331,
  5.5182498282627535,
  43.827979250495424,
  59.32039552962836,
  71.31527646754795,
  35.93046759540154,
  12.59109847741283,
  24.788957921633937,
  129.3571454260518,
  11.887722194011053,
  40.99694814611819,
  94.71308968063315,
  36.73503457778565,
  13.97215306463258,
  5.14119302486334,
  163.0448576896018,
  86.75943414643599,
  16.038040980769715,
  15.664011593266858,
  33.00629934178887,
  148.83283057136592,
  110.98418488934249,
  100.3834860382854,
  15.75717528846002,
  12.307760569011332,
  53.475879372101176,
  16.568942792146007,
  99.58280660718562,
  185.49261717894046,
  76.75495121468843,
  8.804101245994955,
  12.072985631697277,
  41.06437663797372,
  15.778021398277152,
  24.57910530495824,
  9.708315924291124,
  21.649317940346812,
  11.242234913006676,
  7.466023442864695,
  90.37652339293187,
  29.41549523078106,
  35.862949981336186,
  54.91764088297558,
  60.21665190189594,
  12.016740812354248,
  36.34313114773403,
  22.91076117181387,
  9.464872222184768,
  83.44011472879953,
  118.60474671374632,
  16.90482374613588,
  104.61256431685752,
  72.95280664984845,
  6.268585572900139,
  20.713233602278045,
  80.39620366502926,
  52.181210928792396,
  6.938050559871266,
  20.648118035248583,
  30.99681502206858,
  49.628281036959365,
  9.584193167767525,
  16.509802291189768,
  125.15645710139621,
  25.344732339623746,
  15.977223526087242,
  14.286513788791183,
  96.5675655045218,
  58.23532270886482,
  47.68076013957206,
  16.018604745092354,
  26.054669212727973,
  6.652869136204719,
  54.29843236753274,
  14.414812044652408,
  16.639164896276604,
  71.41174692229937,
  20.833132629259005,
  68.45718923072876,
  33.26616913682368,
  20.047345746189144,
  27.934179038482363,
  19.299217675171178,
  107.58302366828818,
  75.6500306749686,
  6.651848240630262,
  31.944831622204124,
  73.3765233425237,
  39.67320774469785,
  17.440999100494246,
  34.44164294826572,
  126.93734537630989,
  16.620727697732978,
  69.54665555127444,
  60.45629049673758,
  5.63105086912868,
  19.522051187072236,
  62.27096149541468,
  60.67195590491417,
  86.52526589928303,
  38.61207726252993,
  7.532671805134615,
  50.48846472674089,
  19.07056108463937,
  28.420606500143194,
  36.757386983235634,
  10.923075603192402,
  5.3225035018966596,
  33.21727700725195,
  54.134521734610104,
  12.75159216306036,
  128.43322641248153,
  54.06743389099918,
  101.81941678298769,
  5.101439543090858,
  16.94303162038046,
  27.3345310750813,
  10.484055949992534,
  29.616859538365,
  148.77454016864962,
  62.74704913622194,
  97.1840066969832,
  44.560913847901304,
  5.589348195126645,
  16.590903496044287,
  30.462369070591823,
  13.418757577771217,
  30.75327493158142,
  15.888223149917332,
  46.660505016633586,
  45.19916506197506,
  22.215699891678508,
  117.4328234780254,
  30.18347214809416,
  104.71470924944434,
  5.984992075253884,
  48.35175979250335,
  46.84055650898214,
  58.677855926794905,
  38.98927602372218,
  93.2079691839044,
  85.63224558858346,
  5.888573005613501,
  12.593361415831872,
  87.74769253283105,
  19.710567312575776,
  26.34004550831608,
  13.495745703797754,
  27.70258630395062,
  13.568252275120846,
  158.12285374798728,
  39.73804158175106,
  29.626687211705637,
  5.914663711753748,
  75.04314916062881,
  9.315255833188704,
  43.87834060130922,
  31.44055115826913,
  12.265292996886549,
  22.490785312920128,
  43.63774794211784,
  11.71939375103589,
  8.262439396794834,
  9.105780844777515,
  46.50149903319378,
  15.887104693677763,
  63.99154882662474,
  101.09270443189122,
  29.939639858062307,
  19.507409018022656,
  74.69168717780781,
  28.089010654182335,
  23.25564009212045,
  5.085162114670139,
  81.78955290350748,
  42.66747885033354,
  106.00899290103145,
  29.715198934034966,
  67.27489565399378,
  48.80267760135854,
  34.16138395616179,
  18.007079459133735,
  21.84921166905134,
  38.96767810201118,
  80.75319523370256,
  13.347909526829145,
  24.594326419619406,
  10.561663315977656,
  6.624031667018189,
  75.94491658220176,
  18.79511353010632,
  39.79939997999378,
  33.36096152171195,
  101.36317860277666,
  113.76091662574385,
  24.640488746882934,
  86.66010868403646,
  5.706687226396635,
  44.32853366457902,
  5.51405171921014,
  108.70895805014075,
  5.5310254389849804,
  70.41057491272504,
  47.12389223818455,
  108.6989571892436,
  182.4737320633952,
  29.998573992380294,
  49.69432611728146,
  5.029258328330005,
  32.304414962535944,
  17.36997124214708,
  121.19064655191032,
  101.87795678139386,
  45.04779284734584,
  48.593768426363894,
  35.64057058169847,
  107.42594130365052,
  31.45399264785776,
  21.200279346067923,
  8.461430234884595,
  24.283934472678013,
  82.48528117232584,
  32.821648162731165,
  57.605965967311306,
  43.30185708596995,
  45.644495911362725,
  17.320360771068028,
  183.16915853507663,
  8.520177703519378,
  18.312309167844806,
  121.68491974258711,
  42.862973821472316,
  47.50439172519498,
  99.57952092077544,
  54.71425161718154,
  34.623437296487715,
  20.865040033802842,
  11.440674212609661,
  17.025758523183047,
  30.839824659911425,
  90.49740866544097,
  27.157401723102364,
  86.66107886847195,
  33.56837965661565,
  54.88772738836616,
  50.96950534021002,
  24.65996221703349,
  37.96578564009686,
  63.34338405086122,
  43.34534683533772,
  13.98648715381339,
  32.32862240491876,
  94.97892029322108,
  16.403567208079387,
  85.04221522791974,
  77.79832017093501,
  6.5854865470334225,
  10.807054428923905,
  81.93777765066649,
  6.976146958220431,
  46.228462547775315,
  38.6641998882411,
  17.736808397770282,
  101.74583667175969,
  59.2889039077969,
  45.20372525382959,
  127.32025166287379,
  63.21859304486967,
  31.281834140319642,
  12.078817022313709,
  19.263869615970876,
  32.991081530648486,
  49.64258354205053,
  8.374601105832745,
  7.029297240494546,
  16.558432035663078,
  53.82449793540092,
  43.566399940304734,
  91.67528624738806,
  9.157749983205344,
  6.579909981910118,
  17.56094362280539,
  40.265505747996045,
  56.573520100534665,
  17.413538214154975,
  12.7797434429684,
  32.35867425774145,
  77.82026538818306,
  5.121533518980227,
  101.7578709870997,
  77.48368832063628,
  55.202222263927275,
  18.626238278642226,
  17.76215976406651,
  6.34910612362186,
  100.8895039327821,
  25.926771694425145,
  17.896138230928678,
  32.41485772983219,
  25.21402044751167,
  19.29398340171033,
  10.248217214793396,
  9.791959674920195,
  7.343725845564252,
  52.15446129367068,
  24.959982289937354,
  84.07555915953853,
  17.0113708133652,
  9.03658189124093,
  61.84057891527404,
  9.538604643810766,
  19.05735813538856,
  18.052449654779608,
  134.6959197010295,
  36.04259739875368,
  5.615267558339154,
  73.24430616310883,
  52.63145940532872,
  53.47126018274981,
  32.51113391013275,
  32.22240102863999,
  86.10351476706427,
  25.389287956549097,
  142.8491652309866,
  53.103076985387894,
  37.22687275226063,
  31.21248237943917,
  19.83284630745572,
  53.417590569020206,
  12.524809533106545,
  173.8357079722619,
  25.77554049379744,
  110.00950238043289,
  166.06392545198864,
  5.288620049276263,
  35.217247713507675,
  10.946266471633303,
  8.534230151249488,
  122.96767232576417,
  11.847801286111878,
  108.28968837739178,
  67.93743992899132,
  44.700093264336886,
  18.937412132932778,
  5.242587642431672,
  84.0679105394437,
  31.62858994253959,
  19.460901129165595,
  18.946922394244822,
  190.24575871275707,
  76.78417289047421,
  59.02407556178068,
  32.69123832598587,
  69.3843738199064,
  7.305782991260869,
  14.679698685690322,
  51.32037673368312,
  8.769510196013407,
  6.0803622344308685,
  41.08523329336502,
  6.117089677843529,
  45.86859084740329,
  181.14611119106755,
  29.55674198479838,
  37.99132297661872,
  25.980361459516537,
  12.50473661248837,
  52.30442459901546,
  7.494578440218566,
  74.4312715877626,
  35.26890893853241,
  54.144150710916286,
  6.435521839743759,
  18.174198918727328,
  54.50132881917325,
  6.512541534962708,
  73.1355307334673])

Energy spectrum

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

([3364696.875,
  3587743.75,
  3917159.375,
  4005256.25,
  4062275.0,
  4279912.5,
  4419003.125,
  5943579.6875,
  5946803.125,
  5976665.625,
  5977243.75,
  6036115.625,
  6060587.5,
  6162345.3125,
  6166395.3125,
  6167187.5,
  6167535.9375,
  6169118.75,
  6173065.625,
  6173517.1875,
  6174679.6875,
  6189114.0625,
  6189367.1875,
  6194743.75,
  6200681.25,
  6207485.9375,
  6229645.3125,
  6243259.375,
  6279590.625,
  6304725.0,
  6326243.75,
  6336606.25,
  6359184.375,
  6359246.875,
  6362032.8125,
  6364492.1875,
  6364865.625,
  6364917.1875,
  6367781.25,
  6368956.25,
  6370382.8125,
  6371445.3125,
  6375167.1875,
  6375825.0,
  6390464.0625,
  6402471.875,
  6406246.875,
  6408206.25,
  6428725.0,
  6437965.625,
  6438481.25,
  6438542.1875,
  6440595.3125,
  6445176.5625,
  6446743.75,
  6447146.875,
  6450457.8125,
  6459673.4375,
  6464581.25,
  6465748.4375,
  6466065.625,
  6466120.3125,
  6467882.8125,
  6468570.3125,
  6468723.4375,
  6469321.875,
  6470251.5625,
  6470309.375,
  6471634.375,
  6471650.0,
  6472171.875,
  6472173.4375,
  6473368.75,
  6473637.5,
  6474581.25,
  6474662.5,
  6474790.625,
  6475096.875,
  6475337.5,
  6475690.625,
  6475828.125,
  6485023.4375,
  6487828.125,
  6487896.875,
  6492175.0,
  6497145.3125,
  6498739.0625,
  6501164.0625,
  6501520.3125,
  6501609.375,
  6502396.875,
  6502448.4375,
  6503942.1875,
  6505600.0,
  6506640.625,
  6509953.125,
  6510768.75,
  6511710.9375,
  6512196.875,
  6512881.25,
  6513687.5,
  6513792.1875,
  6515100.0,
  6515623.4375,
  6515870.3125,
  6515882.8125,
  6515904.6875,
  6516553.125,
  6516596.875,
  6516915.625,
  6516987.5,
  6517542.1875,
  6517921.875,
  6518796.875,
  6518921.875,
  6519176.5625,
  6520000.0,
  6521760.9375,
  6523034.375,
  6523367.1875,
  6523532.8125,
  6523548.4375,
  6525879.6875,
  6525935.9375,
  6532328.125,
  6540443.75,
  6544451.5625,
  6546923.4375,
  6547420.3125,
  6549707.8125,
  6553698.4375,
  6557842.1875,
  6559315.625,
  6559487.5,
  6560742.1875,
  6560928.125,
  6561950.0,
  6561973.4375,
  6562503.125,
  6564526.5625,
  6564675.0,
  6564754.6875,
  6565093.75,
  6565120.3125,
  6565210.9375,
  6565664.0625,
  6565890.625,
  6566762.5,
  6567170.3125,
  6568292.1875,
  6568648.4375,
  6568820.3125,
  6571121.875,
  6571251.5625,
  6571876.5625,
  6572320.3125,
  6573785.9375,
  6574362.5,
  6574495.3125,
  6574496.875,
  6577071.875,
  6577881.25,
  6581834.375,
  6582689.0625,
  6583628.125,
  6587057.8125,
  6590504.6875,
  6592665.625,
  6592682.8125,
  6593306.25,
  6593353.125,
  6593975.0,
  6608518.75,
  6610675.0,
  6610989.0625,
  6611581.25,
  6613792.1875,
  6613826.5625,
  6614840.625,
  6619557.8125,
  6622720.3125,
  6626937.5,
  6628729.6875,
  6631796.875,
  6637106.25,
  6637487.5,
  6639490.625,
  6642393.75,
  6645617.1875,
  6648360.9375,
  6648451.5625,
  6649357.8125,
  6651246.875,
  6653803.125,
  6654814.0625,
  6655835.9375,
  6660771.875,
  6660776.5625,
  6660985.9375,
  6661476.5625,
  6662220.3125,
  6662435.9375,
  6662737.5,
  6663009.375,
  6663775.0,
  6664293.75,
  6664698.4375,
  6664707.8125,
  6664712.5,
  6664732.8125,
  6665246.875,
  6665481.25,
  6665907.8125,
  6666179.6875,
  6667107.8125,
  6667243.75,
  6667243.75,
  6667531.25,
  6667659.375,
  6667659.375,
  6667675.0,
  6667693.75,
  6667948.4375,
  6669767.1875,
  6669856.25,
  6669909.375,
  6671204.6875,
  6671242.1875,
  6672765.625,
  6672789.0625,
  6673117.1875,
  6673551.5625,
  6673593.75,
  6683201.5625,
  6683554.6875,
  6684132.8125,
  6684712.5,
  6684729.6875,
  6685514.0625,
  6686403.125,
  6686450.0,
  6686914.0625,
  6686942.1875,
  6686967.1875,
  6687146.875,
  6687534.375,
  6687709.375,
  6687842.1875,
  6687906.25,
  6687987.5,
  6688120.3125,
  6688150.0,
  6688173.4375,
  6688479.6875,
  6688576.5625,
  6688626.5625,
  6688725.0,
  6688809.375,
  6688832.8125,
  6688839.0625,
  6688839.0625,
  6688907.8125,
  6689046.875,
  6689062.5,
  6689137.5,
  6689223.4375,
  6689251.5625,
  6689256.25,
  6689273.4375,
  6689421.875,
  6689453.125,
  6689568.75,
  6689720.3125,
  6689850.0,
  6689853.125,
  6689895.3125,
  6689925.0,
  6689937.5,
  6690018.75,
  6690087.5,
  6690129.6875,
  6690229.6875,
  6690407.8125,
  6690440.625,
  6690625.0,
  6690628.125,
  6690648.4375,
  6690673.4375,
  6690689.0625,
  6690709.375,
  6690731.25,
  6690850.0,
  6690934.375,
  6690998.4375,
  6691109.375,
  6691165.625,
  6691282.8125,
  6691296.875,
  6691367.1875,
  6691382.8125,
  6691387.5,
  6691398.4375,
  6691412.5,
  6691457.8125,
  6691578.125,
  6691617.1875,
  6691618.75,
  6691746.875,
  6691756.25,
  6691771.875,
  6691870.3125,
  6691995.3125,
  6692050.0,
  6692085.9375,
  6692139.0625,
  6692145.3125,
  6692290.625,
  6692312.5,
  6692314.0625,
  6692348.4375,
  6692354.6875,
  6692365.625,
  6692465.625,
  6692543.75,
  6692548.4375,
  6692664.0625,
  6692703.125,
  6692753.125,
  6693200.0,
  6693221.875,
  6693418.75,
  6693481.25,
  6693589.0625,
  6693678.125,
  6693714.0625,
  6693890.625,
  6693945.3125,
  6694150.0,
  6694182.8125,
  6694321.875,
  6694412.5,
  6694496.875,
  6694518.75,
  6694560.9375,
  6694573.4375,
  6694604.6875,
  6694606.25,
  6694628.125,
  6694867.1875,
  6694929.6875,
  6694965.625,
  6695107.8125,
  6695304.6875,
  6695348.4375,
  6695868.75,
  6695875.0,
  6695934.375,
  6696343.75,
  6696643.75,
  6698365.625,
  6698726.5625,
  6699153.125,
  6699443.75,
  6700907.8125,
  6701806.25,
  6703776.5625,
  6704851.5625,
  6704968.75,
  6705842.1875,
  6705932.8125,
  6706200.0,
  6706468.75,
  6706506.25,
  6706637.5,
  6706678.125,
  6706687.5,
  6707142.1875,
  6708037.5,
  6708206.25,
  6713785.9375,
  6715167.1875,
  6715550.0,
  6716701.5625,
  6717729.6875,
  6719778.125,
  6719892.1875,
  6720403.125,
  6720643.75,
  6721895.3125,
  6722748.4375,
  6722806.25,
  6722854.6875,
  6722859.375,
  6723264.0625,
  6723445.3125,
  6723500.0,
  6723537.5,
  6723787.5,
  6724260.9375,
  6724289.0625,
  6725371.875,
  6725581.25,
  6725668.75,
  6725821.875,
  6725842.1875,
  6725887.5,
  6725910.9375,
  6726640.625,
  6726971.875,
  6727012.5,
  6727056.25,
  6727259.375,
  6727626.5625,
  6728512.5,
  6728743.75,
  6728762.5,
  6731779.6875,
  6732260.9375,
  6732432.8125,
  6732892.1875,
  6733714.0625,
  6734140.625,
  6734910.9375,
  6735528.125,
  6738829.6875,
  6738943.75,
  6739029.6875,
  6739410.9375,
  6739643.75,
  6740900.0,
  6742184.375,
  6742465.625,
  6743915.625,
  6744950.0,
  6745137.5,
  6746487.5,
  6746514.0625,
  6747059.375,
  6750815.625,
  6752242.1875,
  6752393.75,
  6752410.9375,
  6752435.9375,
  6753500.0,
  6754150.0,
  6758209.375,
  6759045.3125,
  6759801.5625,
  6760979.6875,
  6761115.625,
  6775076.5625,
  6775379.6875,
  6779660.9375,
  6782918.75,
  6788901.5625,
  6790290.625,
  6795151.5625,
  6805556.25,
  6807815.625,
  6822965.625,
  6823621.875,
  6823632.8125,
  6828728.125,
  6829656.25,
  6841743.75,
  6841768.75,
  6845446.875,
  6845564.0625,
  6846368.75,
  6846379.6875,
  6846431.25,
  6846573.4375,
  6846617.1875,
  6846870.3125,
  6847267.1875,
  6849175.0,
  6849739.0625,
  6850175.0,
  6850481.25,
  6852170.3125,
  6852385.9375,
  6868739.0625,
  6868917.1875,
  6882912.5,
  6888640.625,
  6889214.0625,
  6889606.25,
  6889843.75,
  6889950.0,
  6890810.9375,
  6891292.1875,
  6891485.9375,
  6892315.625,
  6892803.125,
  6892867.1875,
  6893407.8125,
  6895385.9375,
  6895560.9375,
  6896421.875,
  6897648.4375,
  6897671.875,
  6897778.125,
  6897896.875,
  6898139.0625,
  6898165.625,
  6898393.75,
  6898721.875,
  6898735.9375,
  6898770.3125,
  6898815.625,
  6898945.3125,
  6899265.625,
  6899350.0,
  6899509.375,
  6899554.6875,
  6899812.5,
  6900217.1875,
  6900429.6875,
  6900546.875,
  6902576.5625,
  6902656.25,
  6906939.0625,
  6913245.3125,
  6913442.1875,
  6913759.375,
  6913784.375,
  6914006.25,
  6923751.5625,
  6925965.625,
  6927595.3125,
  6932939.0625,
  6938412.5,
  6938500.0,
  6943160.9375,
  6945146.875,
  6950165.625,
  6951968.75,
  6957435.9375,
  6960814.0625,
  6964254.6875,
  6964367.1875,
  6967328.125,
  6967510.9375,
  6969062.5,
  6969065.625,
  6970078.125,
  6970951.5625,
  6970954.6875,
  6972007.8125,
  6972178.125,
  6972671.875,
  6973873.4375,
  6975665.625,
  6975689.0625,
  6975725.0,
  6976389.0625,
  6976860.9375,
  6977046.875,
  6977210.9375,
  6977351.5625,
  6978246.875,
  6978254.6875,
  6978870.3125,
  6979339.0625,
  6980842.1875,
  6980962.5,
  6983843.75,
  6984129.6875,
  6984187.5,
  6985410.9375,
  6985442.1875,
  6985471.875,
  6985501.5625,
  6985843.75,
  6986506.25,
  6987567.1875,
  6988475.0,
  6988809.375,
  6988904.6875,
  6991437.5,
  6991546.875,
  6992840.625,
  6993178.125,
  6993340.625,
  6994176.5625,
  6994184.375,
  6994189.0625,
  6994957.8125,
  6995317.1875,
  6995326.5625,
  6995410.9375,
  6995951.5625,
  6995968.75,
  6998239.0625,
  6998546.875,
  7000923.4375,
  7001118.75,
  7001329.6875,
  7002537.5,
  7002551.5625,
  7003082.8125,
  7003625.0,
  7003793.75,
  7003798.4375,
  7003840.625,
  7003860.9375,
  7004200.0,
  7005406.25,
  7007182.8125,
  7011151.5625,
  7011412.5,
  7011417.1875,
  7011759.375,
  7013165.625,
  7014240.625,
  7015301.5625,
  7018817.1875,
  7019067.1875,
  7021000.0,
  7022381.25,
  7022514.0625,
  7022531.25,
  7022943.75,
  7023767.1875,
  7024992.1875,
  7025025.0,
  7025112.5,
  7025460.9375,
  7025710.9375,
  7025939.0625,
  7025998.4375,
  7026079.6875,
  7026495.3125,
  7027435.9375,
  7027623.4375,
  7028498.4375,
  7028890.625,
  7029484.375,
  7030900.0,
  7031848.4375,
  7032140.625,
  7032292.1875,
  7032504.6875,
  7033356.25,
  7034859.375,
  7034962.5,
  7034968.75,
  7035606.25,
  7035607.8125,
  7036428.125,
  7037121.875,
  7037123.4375,
  7037503.125,
  7038662.5,
  7038715.625,
  7039207.8125,
  7039254.6875,
  7039350.0,
  7039354.6875,
  7039662.5,
  7039820.3125,
  7039956.25,
  7040417.1875,
  7040514.0625,
  7040625.0,
  7040729.6875,
  7041090.625,
  7041131.25,
  7041134.375,
  7041375.0,
  7041817.1875,
  7041831.25,
  7041843.75,
  7041859.375,
  7041875.0,
  7041893.75,
  7041990.625,
  7042067.1875,
  7042160.9375,
  7042171.875,
  7042221.875,
  7042579.6875,
  7042771.875,
  7042873.4375,
  7043117.1875,
  7043200.0,
  7043381.25,
  7043401.5625,
  7043826.5625,
  7044237.5,
  7044504.6875,
  7057339.0625,
  7058273.4375,
  7060925.0,
  7061534.375,
  7062500.0,
  7063012.5,
  7063031.25,
  7063771.875,
  7065518.75,
  7065803.125,
  7065828.125,
  7066556.25,
  7067959.375,
  7068065.625,
  7068767.1875,
  7070281.25,
  7074481.25,
  7074498.4375,
  7074615.625,
  7074687.5,
  7075543.75,
  7075842.1875,
  7075893.75,
  7076270.3125,
  7076410.9375,
  7076723.4375,
  7078187.5,
  7079053.125,
  7079450.0,
  7080343.75,
  7085562.5,
  7087193.75,
  7087200.0,
  7087267.1875,
  7087456.25,
  7091034.375,
  7097928.125,
  7098204.6875,
  7098815.625,
  7101859.375,
  7103120.3125,
  7108912.5,
  7108915.625,
  7110834.375,
  7117862.5,
  7121410.9375,
  7122609.375,
  7123043.75,
  7123385.9375,
  7123467.1875,
  7123743.75,
  7124987.5,
  7126746.875,
  7127312.5,
  7128600.0,
  7129276.5625,
  7129406.25,
  7130171.875,
  7130509.375,
  7131173.4375,
  7132093.75,
  7132645.3125,
  7134546.875,
  7137681.25,
  7138637.5,
  7142192.1875,
  7147798.4375,
  7148698.4375,
  7148828.125,
  7149118.75,
  7149551.5625,
  7151535.9375,
  7151878.125,
  7152326.5625,
  7152554.6875,
  7153806.25,
  7153823.4375,
  7155406.25,
  7156182.8125,
  7156943.75,
  7157062.5,
  7157703.125,
  7158659.375,
  7159415.625,
  7162153.125,
  7163025.0,
  7163293.75,
  7165400.0,
  7165432.8125,
  7165526.5625,
  7168846.875,
  7172923.4375,
  7174017.1875,
  7176732.8125,
  7178112.5,
  7181457.8125,
  7182885.9375,
  7184625.0,
  7186610.9375,
  7186668.75,
  7188978.125,
  7189812.5,
  7191117.1875,
  7193868.75,
  7195989.0625,
  7199860.9375,
  7201225.0,
  7202990.625,
  7204723.4375,
  7208831.25,
  7209640.625,
  7209660.9375,
  7210468.75,
  7210945.3125,
  7211884.375,
  7212140.625,
  7212409.375,
  7212729.6875,
  7213804.6875,
  7214389.0625,
  7214554.6875,
  7214753.125,
  7214789.0625,
  7214909.375,
  7215381.25,
  7215393.75,
  7215409.375,
  7215467.1875,
  7215560.9375,
  7215643.75,
  7215696.875,
  7215767.1875,
  7215785.9375,
  7215818.75,
  7215867.1875,
  7215906.25,
  7215917.1875,
  7215950.0,
  7216175.0,
  7216357.8125,
  7216450.0,
  7216520.3125,
  7216676.5625,
  7216940.625,
  7216970.3125,
  7217082.8125,
  7217226.5625,
  7217642.1875,
  7217714.0625,
  7218464.0625,
  7219112.5,
  7219385.9375,
  7221445.3125,
  7222723.4375,
  7231503.125,
  7243095.3125,
  7245731.25,
  7245950.0,
  7246095.3125,
  7248907.8125,
  7250520.3125,
  7258990.625,
  7260260.9375,
  7260695.3125,
  7261003.125,
  7261865.625,
  7266873.4375,
  7282703.125,
  7283145.3125,
  7284800.0,
  7286757.8125,
  7286831.25,
  7291665.625,
  7294006.25,
  7309860.9375,
  7312632.8125,
  7313651.5625,
  7314429.6875,
  7315393.75,
  7320562.5,
  7322117.1875,
  7328507.8125,
  7354607.8125,
  7354776.5625,
  7355726.5625,
  7356200.0,
  7356939.0625,
  7357431.25,
  7394865.625,
  7434064.0625,
  7576507.8125,
  7585914.0625,
  7620604.6875,
  7620907.8125,
  7624446.875,
  7653101.5625,
  7683304.6875,
  7684326.5625,
  7684378.125,
  7714845.3125,
  7736729.6875,
  7740359.375,
  7752598.4375,
  7753520.3125,
  7753707.8125,
  7778196.875,
  7781278.125,
  7781287.5,
  7781371.875,
  7879235.9375,
  8104754.6875,
  8173246.875,
  8539956.25,
  8551959.375,
  8552020.3125,
  8721107.8125,
  8761750.0,
  8872392.1875,
  8907807.8125,
  8908339.0625,
  8910898.4375,
  8916675.0,
  8920428.125,
  8931807.8125,
  9182890.625,
  9193826.5625,
  9195637.5,
  9261812.5,
  9275353.125,
  9286118.75,
  9290906.25,
  9294509.375,
  9318065.625,
  9352184.375,
  9356853.125,
  9359085.9375,
  9368718.75,
  9376306.25,
  9378681.25,
  9411829.6875,
  9454026.5625,
  9456243.75,
  9456651.5625,
  9456993.75,
  9457723.4375,
  9461206.25,
  9462717.1875,
  9463818.75,
  9474403.125,
  9482407.8125,
  9486654.6875,
  9487507.8125,
  9490128.125,
  9490478.125,
  9506351.5625,
  9507842.1875,
  9524360.9375,
  9554387.5,
  9581104.6875,
  9622239.0625,
  9628129.6875,
  9642587.5,
  9649184.375,
  9666600.0,
  9669510.9375,
  9669579.6875,
  9671387.5,
  9673656.25,
  9688442.1875,
  9711200.0,
  9711279.6875,
  9827278.125,
  9841485.9375,
  9889659.375,
  9893678.125,
  9894756.25,
  9894798.4375,
  9895709.375,
  9895735.9375,
  9897373.4375,
  9897423.4375,
  9936626.5625,
  9962996.875,
  9995064.0625,
  10013256.25,
  10026162.5,
  10026946.875,
  10029226.5625,
  10036209.375,
  10036742.1875,
  10042568.75,
  10042981.25,
  10046739.0625,
  10046743.75,
  10049373.4375,
  10087198.4375,
  10090228.125,
  10116943.75,
  10263607.8125,
  11257718.75,
  11267890.625,
  11267957.8125,
  11280456.25],
 [36.39723457256683,
  61.42011606542019,
  57.63331380236617,
  8.929290614776288,
  5.867756950569329,
  103.15167346425969,
  13.632253326347488,
  61.988275601110445,
  15.782696782507756,
  5.227109715887212,
  43.40145628763871,
  102.45907594437872,
  70.63837596907882,
  128.59382627377136,
  17.587650568578034,
  15.793167642101507,
  7.671813627472112,
  9.006379654396119,
  28.70609851781273,
  100.29981138667708,
  78.80329657045203,
  23.225821689737273,
  88.18321797392157,
  38.34357814715931,
  23.47776991851982,
  5.225485454867674,
  5.659364363078715,
  76.76677326861977,
  8.024261796596173,
  24.275962905745054,
  25.377084321643334,
  83.89258329860951,
  23.946462742652194,
  50.93071164949142,
  58.51405958943931,
  67.06249961993483,
  55.847138471574354,
  14.149140945397248,
  20.544470133004168,
  10.0180757997359,
  68.33173525422447,
  108.77033466229885,
  26.194573928714608,
  8.360969807334461,
  11.591584576740122,
  15.758665732968948,
  34.030938217134114,
  5.481389065289376,
  24.376277770638143,
  42.175379433676426,
  53.61122422065942,
  5.653454013029906,
  15.857975173415213,
  65.9539246648907,
  8.310939051713781,
  43.64728620668559,
  192.57406810768012,
  50.15962257049653,
  34.27502430628671,
  85.83646132075555,
  34.622424313148166,
  131.8763932836117,
  70.45556578208915,
  24.117761348660785,
  19.78736390407145,
  19.24625273551636,
  259.5303046208734,
  9.267275539627011,
  94.52168757313635,
  25.67429902289775,
  17.247923347038633,
  16.99510366372298,
  18.338551681601587,
  13.32147694962425,
  63.72091543156779,
  56.37430669301801,
  50.916425410211986,
  125.0897587093759,
  58.48166225709635,
  16.703579595561827,
  19.11011430135429,
  25.51252779403324,
  32.804139755043394,
  10.697281014364044,
  54.34203430499464,
  19.291934742150524,
  124.8974664064839,
  15.609415701052098,
  5.432567156193797,
  73.14207159300558,
  19.598637027356833,
  7.306629199403174,
  11.296177152251182,
  47.552568072581444,
  36.96907227087646,
  54.25984614013294,
  23.25976152933508,
  45.334650100167025,
  176.35599253956264,
  17.607075955173745,
  82.81061103496214,
  56.31303715647137,
  8.231195933068582,
  93.65401967694291,
  157.1143954299284,
  62.03569970000155,
  7.316580896670568,
  7.052124010949603,
  35.7759027797182,
  30.82360542216003,
  9.798540488645234,
  100.03988013196754,
  11.38314142977542,
  7.258916934182001,
  39.52951734617878,
  52.82393617615527,
  34.045747539411394,
  63.760109539502466,
  12.6494988393673,
  29.320364298043064,
  21.823657080944493,
  7.000983790584903,
  26.192950922073393,
  5.0758953453095605,
  14.418031387502166,
  33.6527390632858,
  76.86622242713834,
  13.311024654878873,
  29.35449753039608,
  65.95837123281866,
  12.413601516380288,
  29.296784189244242,
  97.81449871045096,
  99.83564237066203,
  40.2476894423131,
  51.4244131011592,
  79.48771264409444,
  43.938712458056315,
  13.011741221320914,
  12.867665300025267,
  44.411861212108604,
  5.698327666029522,
  23.220022487333814,
  7.677072247957201,
  123.11936950859335,
  77.70092946670269,
  51.24469015790008,
  206.5094388591854,
  22.04565298695637,
  48.40618971229868,
  29.481152250632725,
  75.84270141599688,
  107.3155784615033,
  20.330654925873958,
  101.113585193938,
  68.58657363330593,
  19.050999659011822,
  15.296627976197115,
  32.56451285716979,
  21.81137862581718,
  59.30814293658601,
  60.708985931199436,
  137.42236923465742,
  129.95393069598427,
  49.771180646118474,
  49.45136136964324,
  24.937731699140002,
  118.24452935403079,
  20.073764992819875,
  153.86990491667837,
  107.63172441621319,
  15.016307833826497,
  8.408690366765153,
  30.312961874645012,
  5.997958863078135,
  8.452608138764365,
  5.5570975373728,
  8.001892567757281,
  73.79057910086591,
  21.111881574438257,
  6.066153738741113,
  47.08419560585291,
  56.67510204822391,
  43.71788249107159,
  82.39327886377018,
  57.65663076026064,
  13.303113273263156,
  10.330633974957554,
  5.755159609170377,
  19.040784110635492,
  5.621669378083233,
  69.64824848732714,
  33.21099769977836,
  62.609628259637134,
  69.44656875655085,
  19.917878241918054,
  73.48391273876348,
  15.59866757654158,
  28.525795935709684,
  16.450191714611123,
  19.085978566616227,
  9.847199642537985,
  13.939621181470327,
  17.2774418111779,
  10.45611498110723,
  8.539391555424086,
  15.47772590263954,
  51.50271746735625,
  27.51102222060111,
  39.995722741166865,
  10.561223311222646,
  5.49784137905199,
  32.066418116179484,
  96.33037548367857,
  39.25960964222006,
  17.682133427314476,
  96.80524359126099,
  111.17831515119403,
  11.504910229982594,
  16.707816525480926,
  71.48760686188146,
  27.18776603881378,
  38.39960078758963,
  16.4254365726055,
  60.08746596578121,
  45.83385414665403,
  23.353312654207613,
  58.85832501460436,
  41.66846193034802,
  25.889133080981267,
  27.0149215014014,
  112.80447085716244,
  10.512909891136907,
  13.662628979492569,
  70.4005236572396,
  74.16656977228209,
  35.238825696528004,
  28.48184039740719,
  10.924899051597857,
  83.17739913908383,
  16.62831994178219,
  60.086958956855675,
  57.91383276906872,
  20.148785539171563,
  18.427659911220164,
  65.26106321836986,
  7.572559540440718,
  115.73928889501224,
  31.863727267981485,
  82.51041484203607,
  35.361293566290534,
  23.648882747165857,
  56.37135508295516,
  101.0228735666784,
  36.995774727011955,
  239.95573348165058,
  8.574192640784469,
  37.96156166995581,
  89.42040732347566,
  20.22137316037128,
  49.65146830305179,
  20.658802593847145,
  27.283433359698986,
  109.55191884710541,
  23.196824737527788,
  146.457331866237,
  10.2022433282344,
  37.952321428908974,
  13.13483637306454,
  9.545625899412434,
  13.385410919813104,
  38.842488194322364,
  17.783996143897205,
  14.750502476257601,
  6.112562033292657,
  18.771899827565584,
  24.80375480532396,
  78.04036275641678,
  46.72856783120531,
  5.295269109919355,
  93.29325981912626,
  99.98597977415547,
  5.966135850601805,
  14.261884710409893,
  16.048538295470127,
  74.57563138024784,
  5.3234613789131,
  34.57341420045836,
  5.3586324095810705,
  44.485182068559666,
  45.081714863343834,
  16.875670462561438,
  88.20729215095666,
  12.437618944272808,
  19.696139110195446,
  71.00349229313728,
  23.773735673476274,
  105.41442588999934,
  28.946077332056817,
  47.253734933968076,
  5.462575582795934,
  30.145835996026317,
  40.09501662518248,
  11.427284583854847,
  5.55686999448816,
  12.318646529880045,
  17.88578485995733,
  107.07422046782936,
  18.249795060960555,
  32.26400018416813,
  7.930158419839312,
  110.67809866341692,
  100.21268052669818,
  42.514685206195864,
  5.095016070893412,
  61.636161580364984,
  70.08682666651234,
  11.157789722497984,
  29.222816861491555,
  12.971300139419613,
  27.543324638141456,
  7.573531274032749,
  12.324898001981433,
  29.7710779204007,
  9.14804439315284,
  18.360638645955152,
  112.90374119668057,
  132.67254040535673,
  7.162134145727406,
  36.49372334823385,
  6.622016997883784,
  80.68135832342915,
  13.82305552678966,
  77.48743691819745,
  8.92290458525,
  14.580672675474293,
  45.27371879263896,
  7.602347372635561,
  46.92059890289319,
  23.19730413326542,
  7.855230077894405,
  15.048063946899774,
  32.50167566475677,
  16.400917327583734,
  18.337144223995534,
  40.4265125267047,
  6.59677491956983,
  39.655220661190405,
  40.04895538926286,
  56.05252352265611,
  15.331484932504182,
  13.099252738344084,
  22.280172788529264,
  70.29092211381443,
  47.92602721364865,
  195.334377726145,
  5.323626736826369,
  20.677837781591993,
  88.1626383636549,
  10.108983179239246,
  28.152152489077626,
  19.25695419449508,
  54.61985711274021,
  26.572763025421242,
  6.629374543258719,
  25.05963433413907,
  41.21229926008208,
  77.93295679007664,
  250.77590879265006,
  9.035192382103938,
  54.87175112220168,
  66.64111531820397,
  25.91709629781178,
  76.67381992941066,
  53.40916499705752,
  16.340801708160313,
  26.298788136451915,
  214.9276806945125,
  66.19947733960495,
  68.80751329046113,
  63.46934616769979,
  13.537669646613109,
  89.7327046603615,
  102.63149776891655,
  46.08703091630864,
  37.1145750733646,
  21.071852348393904,
  42.392824316576274,
  17.86161228472567,
  42.425304714793135,
  45.48455474609187,
  50.216542413045,
  5.360485862196153,
  19.175318322162177,
  11.012932792400386,
  72.29662011971712,
  135.01674338159373,
  61.074610629654856,
  60.62720243901956,
  15.357241098856736,
  7.759515874011704,
  5.989970420166393,
  5.024611979714108,
  9.060548530053877,
  38.32735034697739,
  62.119420008907156,
  24.17980582768046,
  13.694990792726543,
  57.4125608820768,
  9.893305914485136,
  81.69622814528651,
  21.74521986225286,
  8.750634094959812,
  29.14782057690727,
  63.97159745895686,
  7.130111006729035,
  59.81747671424412,
  12.539494478704636,
  7.291013191795108,
  26.3317224034226,
  72.43517609089575,
  116.95774507794263,
  91.63969363108122,
  15.227370288491633,
  35.72033358393782,
  27.30521793086246,
  5.210992339876186,
  22.58353781299215,
  64.3951221857838,
  13.805659203132972,
  22.311537867308818,
  8.16569461317244,
  107.99538676849032,
  17.077027999769637,
  47.66027591322897,
  41.37640024765623,
  42.07896026079423,
  39.4994745302631,
  69.97118573483006,
  55.399735251854224,
  15.92276263285325,
  17.585232347653264,
  51.09300150269347,
  141.83538047771194,
  25.860473614439737,
  51.70004967125205,
  138.24599936323148,
  6.296551660553312,
  27.066458099730465,
  66.53279747520995,
  81.67574145991111,
  115.38709229549974,
  18.80387966381403,
  43.13626442991365,
  117.90100940387464,
  193.27713109874847,
  106.92378733288527,
  28.925922881684542,
  42.98317903938009,
  104.82569972617829,
  18.09353498595276,
  32.274513881869055,
  45.67891945659834,
  94.57165448160939,
  47.54650021282325,
  14.895962407906907,
  26.768832498528994,
  17.551074736142564,
  30.699116105173232,
  6.935831046905322,
  115.22257517935395,
  7.469597110313175,
  87.76036775135202,
  12.25074032922134,
  41.11555506589667,
  23.83262040248162,
  18.002465826479334,
  23.81477842762854,
  8.838077822794379,
  21.68258229385632,
  28.00852373244167,
  90.11889493377922,
  163.07121689626797,
  20.540108270706625,
  65.80541814133274,
  47.786993726294426,
  29.335042397402365,
  93.86532569676896,
  99.56934917062974,
  17.757248337781288,
  55.39859283739911,
  66.25843115182607,
  118.28008005675935,
  30.671089556711923,
  5.051039471698173,
  70.7305224709122,
  119.90700956452221,
  230.26580581163404,
  25.971573098422343,
  12.149888203213914,
  5.703698023666907,
  152.02059815483406,
  66.6988898443032,
  52.15059682527077,
  28.950333905865026,
  7.908227402268309,
  184.35839680302357,
  117.90609263310223,
  46.24631401665481,
  14.498498191791379,
  54.689428424699976,
  20.35447639425396,
  61.38881776133482,
  38.96333242894729,
  25.431385167612866,
  84.12197391091124,
  12.184693790405046,
  6.327233705329153,
  31.548206682192713,
  15.271782150545015,
  5.238278021608867,
  57.11939878703496,
  97.18729218157577,
  57.077332620175994,
  15.686016685863741,
  17.702169318602504,
  66.30863660574207,
  27.550598380821416,
  163.15543480322236,
  150.88784637190432,
  30.871435262025717,
  91.10225921377655,
  25.816027637756207,
  78.26760028282013,
  26.755804301441145,
  12.876748788353112,
  25.8584221967763,
  24.449481349288803,
  30.709596758514653,
  231.14298469256718,
  18.970093123992815,
  48.1997996113733,
  25.00684685505102,
  29.385623027793265,
  44.436833822308195,
  14.250795298761048,
  207.57083810619093,
  130.15807504395414,
  119.25182099085103,
  24.66519540315926,
  101.99071395184572,
  12.308830987562384,
  57.82716184538903,
  9.731976670841103,
  129.23115679400652,
  16.72269343241506,
  5.423619487829046,
  70.55877603142835,
  19.425067156256375,
  25.788211884679193,
  21.750981785451568,
  172.121249594233,
  68.34847770250651,
  56.21839414178929,
  25.765890989278137,
  78.6671097282135,
  19.638280989855588,
  23.331589709919193,
  30.161396048893327,
  7.708445399624251,
  46.14224247816008,
  5.540215973232311,
  40.7158900641637,
  41.11487049581796,
  40.21710777591162,
  134.14653994173736,
  17.504486033334462,
  52.23380160910488,
  80.08579439489925,
  5.922549508149697,
  19.69870709776275,
  10.377037135534696,
  24.217721677720853,
  70.85770556994306,
  63.97574868434822,
  53.62949064490216,
  15.11559126517665,
  51.747009104773994,
  81.054476769162,
  43.86776692651202,
  202.6870560424551,
  18.108006834507794,
  7.7722192158968655,
  6.105094811494561,
  5.373087842024246,
  14.050546904503438,
  32.859317000322775,
  18.099902704058707,
  121.00330014735792,
  45.78391145108671,
  99.5478504324995,
  102.88623227747873,
  116.5460079068383,
  27.886941051564182,
  66.91752860174782,
  5.080617267890299,
  108.49749951238559,
  87.86397757255472,
  7.114315597064646,
  21.569724714252455,
  23.42700522992944,
  26.796810870875007,
  25.071296640276515,
  132.63278955868236,
  8.412439098701245,
  28.48464565877721,
  36.88131302685445,
  67.218326194331,
  5.5182498282627535,
  43.827979250495424,
  59.32039552962836,
  71.31527646754795,
  35.93046759540154,
  12.59109847741283,
  24.788957921633937,
  129.3571454260518,
  11.887722194011053,
  40.99694814611819,
  94.71308968063315,
  36.73503457778565,
  13.97215306463258,
  5.14119302486334,
  163.0448576896018,
  86.75943414643599,
  16.038040980769715,
  15.664011593266858,
  33.00629934178887,
  148.83283057136592,
  110.98418488934249,
  100.3834860382854,
  15.75717528846002,
  12.307760569011332,
  53.475879372101176,
  16.568942792146007,
  99.58280660718562,
  185.49261717894046,
  76.75495121468843,
  8.804101245994955,
  12.072985631697277,
  41.06437663797372,
  15.778021398277152,
  24.57910530495824,
  9.708315924291124,
  21.649317940346812,
  11.242234913006676,
  7.466023442864695,
  90.37652339293187,
  29.41549523078106,
  35.862949981336186,
  54.91764088297558,
  60.21665190189594,
  12.016740812354248,
  36.34313114773403,
  22.91076117181387,
  9.464872222184768,
  83.44011472879953,
  118.60474671374632,
  16.90482374613588,
  104.61256431685752,
  72.95280664984845,
  6.268585572900139,
  20.713233602278045,
  80.39620366502926,
  52.181210928792396,
  6.938050559871266,
  20.648118035248583,
  30.99681502206858,
  49.628281036959365,
  9.584193167767525,
  16.509802291189768,
  125.15645710139621,
  25.344732339623746,
  15.977223526087242,
  14.286513788791183,
  96.5675655045218,
  58.23532270886482,
  47.68076013957206,
  16.018604745092354,
  26.054669212727973,
  6.652869136204719,
  54.29843236753274,
  14.414812044652408,
  16.639164896276604,
  71.41174692229937,
  20.833132629259005,
  68.45718923072876,
  33.26616913682368,
  20.047345746189144,
  27.934179038482363,
  19.299217675171178,
  107.58302366828818,
  75.6500306749686,
  6.651848240630262,
  31.944831622204124,
  73.3765233425237,
  39.67320774469785,
  17.440999100494246,
  34.44164294826572,
  126.93734537630989,
  16.620727697732978,
  69.54665555127444,
  60.45629049673758,
  5.63105086912868,
  19.522051187072236,
  62.27096149541468,
  60.67195590491417,
  86.52526589928303,
  38.61207726252993,
  7.532671805134615,
  50.48846472674089,
  19.07056108463937,
  28.420606500143194,
  36.757386983235634,
  10.923075603192402,
  5.3225035018966596,
  33.21727700725195,
  54.134521734610104,
  12.75159216306036,
  128.43322641248153,
  54.06743389099918,
  101.81941678298769,
  5.101439543090858,
  16.94303162038046,
  27.3345310750813,
  10.484055949992534,
  29.616859538365,
  148.77454016864962,
  62.74704913622194,
  97.1840066969832,
  44.560913847901304,
  5.589348195126645,
  16.590903496044287,
  30.462369070591823,
  13.418757577771217,
  30.75327493158142,
  15.888223149917332,
  46.660505016633586,
  45.19916506197506,
  22.215699891678508,
  117.4328234780254,
  30.18347214809416,
  104.71470924944434,
  5.984992075253884,
  48.35175979250335,
  46.84055650898214,
  58.677855926794905,
  38.98927602372218,
  93.2079691839044,
  85.63224558858346,
  5.888573005613501,
  12.593361415831872,
  87.74769253283105,
  19.710567312575776,
  26.34004550831608,
  13.495745703797754,
  27.70258630395062,
  13.568252275120846,
  158.12285374798728,
  39.73804158175106,
  29.626687211705637,
  5.914663711753748,
  75.04314916062881,
  9.315255833188704,
  43.87834060130922,
  31.44055115826913,
  12.265292996886549,
  22.490785312920128,
  43.63774794211784,
  11.71939375103589,
  8.262439396794834,
  9.105780844777515,
  46.50149903319378,
  15.887104693677763,
  63.99154882662474,
  101.09270443189122,
  29.939639858062307,
  19.507409018022656,
  74.69168717780781,
  28.089010654182335,
  23.25564009212045,
  5.085162114670139,
  81.78955290350748,
  42.66747885033354,
  106.00899290103145,
  29.715198934034966,
  67.27489565399378,
  48.80267760135854,
  34.16138395616179,
  18.007079459133735,
  21.84921166905134,
  38.96767810201118,
  80.75319523370256,
  13.347909526829145,
  24.594326419619406,
  10.561663315977656,
  6.624031667018189,
  75.94491658220176,
  18.79511353010632,
  39.79939997999378,
  33.36096152171195,
  101.36317860277666,
  113.76091662574385,
  24.640488746882934,
  86.66010868403646,
  5.706687226396635,
  44.32853366457902,
  5.51405171921014,
  108.70895805014075,
  5.5310254389849804,
  70.41057491272504,
  47.12389223818455,
  108.6989571892436,
  182.4737320633952,
  29.998573992380294,
  49.69432611728146,
  5.029258328330005,
  32.304414962535944,
  17.36997124214708,
  121.19064655191032,
  101.87795678139386,
  45.04779284734584,
  48.593768426363894,
  35.64057058169847,
  107.42594130365052,
  31.45399264785776,
  21.200279346067923,
  8.461430234884595,
  24.283934472678013,
  82.48528117232584,
  32.821648162731165,
  57.605965967311306,
  43.30185708596995,
  45.644495911362725,
  17.320360771068028,
  183.16915853507663,
  8.520177703519378,
  18.312309167844806,
  121.68491974258711,
  42.862973821472316,
  47.50439172519498,
  99.57952092077544,
  54.71425161718154,
  34.623437296487715,
  20.865040033802842,
  11.440674212609661,
  17.025758523183047,
  30.839824659911425,
  90.49740866544097,
  27.157401723102364,
  86.66107886847195,
  33.56837965661565,
  54.88772738836616,
  50.96950534021002,
  24.65996221703349,
  37.96578564009686,
  63.34338405086122,
  43.34534683533772,
  13.98648715381339,
  32.32862240491876,
  94.97892029322108,
  16.403567208079387,
  85.04221522791974,
  77.79832017093501,
  6.5854865470334225,
  10.807054428923905,
  81.93777765066649,
  6.976146958220431,
  46.228462547775315,
  38.6641998882411,
  17.736808397770282,
  101.74583667175969,
  59.2889039077969,
  45.20372525382959,
  127.32025166287379,
  63.21859304486967,
  31.281834140319642,
  12.078817022313709,
  19.263869615970876,
  32.991081530648486,
  49.64258354205053,
  8.374601105832745,
  7.029297240494546,
  16.558432035663078,
  53.82449793540092,
  43.566399940304734,
  91.67528624738806,
  9.157749983205344,
  6.579909981910118,
  17.56094362280539,
  40.265505747996045,
  56.573520100534665,
  17.413538214154975,
  12.7797434429684,
  32.35867425774145,
  77.82026538818306,
  5.121533518980227,
  101.7578709870997,
  77.48368832063628,
  55.202222263927275,
  18.626238278642226,
  17.76215976406651,
  6.34910612362186,
  100.8895039327821,
  25.926771694425145,
  17.896138230928678,
  32.41485772983219,
  25.21402044751167,
  19.29398340171033,
  10.248217214793396,
  9.791959674920195,
  7.343725845564252,
  52.15446129367068,
  24.959982289937354,
  84.07555915953853,
  17.0113708133652,
  9.03658189124093,
  61.84057891527404,
  9.538604643810766,
  19.05735813538856,
  18.052449654779608,
  134.6959197010295,
  36.04259739875368,
  5.615267558339154,
  73.24430616310883,
  52.63145940532872,
  53.47126018274981,
  32.51113391013275,
  32.22240102863999,
  86.10351476706427,
  25.389287956549097,
  142.8491652309866,
  53.103076985387894,
  37.22687275226063,
  31.21248237943917,
  19.83284630745572,
  53.417590569020206,
  12.524809533106545,
  173.8357079722619,
  25.77554049379744,
  110.00950238043289,
  166.06392545198864,
  5.288620049276263,
  35.217247713507675,
  10.946266471633303,
  8.534230151249488,
  122.96767232576417,
  11.847801286111878,
  108.28968837739178,
  67.93743992899132,
  44.700093264336886,
  18.937412132932778,
  5.242587642431672,
  84.0679105394437,
  31.62858994253959,
  19.460901129165595,
  18.946922394244822,
  190.24575871275707,
  76.78417289047421,
  59.02407556178068,
  32.69123832598587,
  69.3843738199064,
  7.305782991260869,
  14.679698685690322,
  51.32037673368312,
  8.769510196013407,
  6.0803622344308685,
  41.08523329336502,
  6.117089677843529,
  45.86859084740329,
  181.14611119106755,
  29.55674198479838,
  37.99132297661872,
  25.980361459516537,
  12.50473661248837,
  52.30442459901546,
  7.494578440218566,
  74.4312715877626,
  35.26890893853241,
  54.144150710916286,
  6.435521839743759,
  18.174198918727328,
  54.50132881917325,
  6.512541534962708,
  73.1355307334673])

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

([3364696.875,
  3587743.75,
  3917159.375,
  4005256.25,
  4062275.0,
  4279912.5,
  4419003.125,
  5943579.6875,
  5946803.125,
  5976665.625,
  5977243.75,
  6036115.625,
  6060587.5,
  6162345.3125,
  6166395.3125,
  6167187.5,
  6167535.9375,
  6169118.75,
  6173065.625,
  6173517.1875,
  6174679.6875,
  6189114.0625,
  6189367.1875,
  6194743.75,
  6200681.25,
  6207485.9375,
  6229645.3125,
  6243259.375,
  6279590.625,
  6304725.0,
  6326243.75,
  6336606.25,
  6359184.375,
  6359246.875,
  6362032.8125,
  6364492.1875,
  6364865.625,
  6364917.1875,
  6367781.25,
  6368956.25,
  6370382.8125,
  6371445.3125,
  6375167.1875,
  6375825.0,
  6390464.0625,
  6402471.875,
  6406246.875,
  6408206.25,
  6428725.0,
  6437965.625,
  6438481.25,
  6438542.1875,
  6440595.3125,
  6445176.5625,
  6446743.75,
  6447146.875,
  6450457.8125,
  6459673.4375,
  6464581.25,
  6465748.4375,
  6466065.625,
  6466120.3125,
  6467882.8125,
  6468570.3125,
  6468723.4375,
  6469321.875,
  6470251.5625,
  6470309.375,
  6471634.375,
  6471650.0,
  6472171.875,
  6472173.4375,
  6473368.75,
  6473637.5,
  6474581.25,
  6474662.5,
  6474790.625,
  6475096.875,
  6475337.5,
  6475690.625,
  6475828.125,
  6485023.4375,
  6487828.125,
  6487896.875,
  6492175.0,
  6497145.3125,
  6498739.0625,
  6501164.0625,
  6501520.3125,
  6501609.375,
  6502396.875,
  6502448.4375,
  6503942.1875,
  6505600.0,
  6506640.625,
  6509953.125,
  6510768.75,
  6511710.9375,
  6512196.875,
  6512881.25,
  6513687.5,
  6513792.1875,
  6515100.0,
  6515623.4375,
  6515870.3125,
  6515882.8125,
  6515904.6875,
  6516553.125,
  6516596.875,
  6516915.625,
  6516987.5,
  6517542.1875,
  6517921.875,
  6518796.875,
  6518921.875,
  6519176.5625,
  6520000.0,
  6521760.9375,
  6523034.375,
  6523367.1875,
  6523532.8125,
  6523548.4375,
  6525879.6875,
  6525935.9375,
  6532328.125,
  6540443.75,
  6544451.5625,
  6546923.4375,
  6547420.3125,
  6549707.8125,
  6553698.4375,
  6557842.1875,
  6559315.625,
  6559487.5,
  6560742.1875,
  6560928.125,
  6561950.0,
  6561973.4375,
  6562503.125,
  6564526.5625,
  6564675.0,
  6564754.6875,
  6565093.75,
  6565120.3125,
  6565210.9375,
  6565664.0625,
  6565890.625,
  6566762.5,
  6567170.3125,
  6568292.1875,
  6568648.4375,
  6568820.3125,
  6571121.875,
  6571251.5625,
  6571876.5625,
  6572320.3125,
  6573785.9375,
  6574362.5,
  6574495.3125,
  6574496.875,
  6577071.875,
  6577881.25,
  6581834.375,
  6582689.0625,
  6583628.125,
  6587057.8125,
  6590504.6875,
  6592665.625,
  6592682.8125,
  6593306.25,
  6593353.125,
  6593975.0,
  6608518.75,
  6610675.0,
  6610989.0625,
  6611581.25,
  6613792.1875,
  6613826.5625,
  6614840.625,
  6619557.8125,
  6622720.3125,
  6626937.5,
  6628729.6875,
  6631796.875,
  6637106.25,
  6637487.5,
  6639490.625,
  6642393.75,
  6645617.1875,
  6648360.9375,
  6648451.5625,
  6649357.8125,
  6651246.875,
  6653803.125,
  6654814.0625,
  6655835.9375,
  6660771.875,
  6660776.5625,
  6660985.9375,
  6661476.5625,
  6662220.3125,
  6662435.9375,
  6662737.5,
  6663009.375,
  6663775.0,
  6664293.75,
  6664698.4375,
  6664707.8125,
  6664712.5,
  6664732.8125,
  6665246.875,
  6665481.25,
  6665907.8125,
  6666179.6875,
  6667107.8125,
  6667243.75,
  6667243.75,
  6667531.25,
  6667659.375,
  6667659.375,
  6667675.0,
  6667693.75,
  6667948.4375,
  6669767.1875,
  6669856.25,
  6669909.375,
  6671204.6875,
  6671242.1875,
  6672765.625,
  6672789.0625,
  6673117.1875,
  6673551.5625,
  6673593.75,
  6683201.5625,
  6683554.6875,
  6684132.8125,
  6684712.5,
  6684729.6875,
  6685514.0625,
  6686403.125,
  6686450.0,
  6686914.0625,
  6686942.1875,
  6686967.1875,
  6687146.875,
  6687534.375,
  6687709.375,
  6687842.1875,
  6687906.25,
  6687987.5,
  6688120.3125,
  6688150.0,
  6688173.4375,
  6688479.6875,
  6688576.5625,
  6688626.5625,
  6688725.0,
  6688809.375,
  6688832.8125,
  6688839.0625,
  6688839.0625,
  6688907.8125,
  6689046.875,
  6689062.5,
  6689137.5,
  6689223.4375,
  6689251.5625,
  6689256.25,
  6689273.4375,
  6689421.875,
  6689453.125,
  6689568.75,
  6689720.3125,
  6689850.0,
  6689853.125,
  6689895.3125,
  6689925.0,
  6689937.5,
  6690018.75,
  6690087.5,
  6690129.6875,
  6690229.6875,
  6690407.8125,
  6690440.625,
  6690625.0,
  6690628.125,
  6690648.4375,
  6690673.4375,
  6690689.0625,
  6690709.375,
  6690731.25,
  6690850.0,
  6690934.375,
  6690998.4375,
  6691109.375,
  6691165.625,
  6691282.8125,
  6691296.875,
  6691367.1875,
  6691382.8125,
  6691387.5,
  6691398.4375,
  6691412.5,
  6691457.8125,
  6691578.125,
  6691617.1875,
  6691618.75,
  6691746.875,
  6691756.25,
  6691771.875,
  6691870.3125,
  6691995.3125,
  6692050.0,
  6692085.9375,
  6692139.0625,
  6692145.3125,
  6692290.625,
  6692312.5,
  6692314.0625,
  6692348.4375,
  6692354.6875,
  6692365.625,
  6692465.625,
  6692543.75,
  6692548.4375,
  6692664.0625,
  6692703.125,
  6692753.125,
  6693200.0,
  6693221.875,
  6693418.75,
  6693481.25,
  6693589.0625,
  6693678.125,
  6693714.0625,
  6693890.625,
  6693945.3125,
  6694150.0,
  6694182.8125,
  6694321.875,
  6694412.5,
  6694496.875,
  6694518.75,
  6694560.9375,
  6694573.4375,
  6694604.6875,
  6694606.25,
  6694628.125,
  6694867.1875,
  6694929.6875,
  6694965.625,
  6695107.8125,
  6695304.6875,
  6695348.4375,
  6695868.75,
  6695875.0,
  6695934.375,
  6696343.75,
  6696643.75,
  6698365.625,
  6698726.5625,
  6699153.125,
  6699443.75,
  6700907.8125,
  6701806.25,
  6703776.5625,
  6704851.5625,
  6704968.75,
  6705842.1875,
  6705932.8125,
  6706200.0,
  6706468.75,
  6706506.25,
  6706637.5,
  6706678.125,
  6706687.5,
  6707142.1875,
  6708037.5,
  6708206.25,
  6713785.9375,
  6715167.1875,
  6715550.0,
  6716701.5625,
  6717729.6875,
  6719778.125,
  6719892.1875,
  6720403.125,
  6720643.75,
  6721895.3125,
  6722748.4375,
  6722806.25,
  6722854.6875,
  6722859.375,
  6723264.0625,
  6723445.3125,
  6723500.0,
  6723537.5,
  6723787.5,
  6724260.9375,
  6724289.0625,
  6725371.875,
  6725581.25,
  6725668.75,
  6725821.875,
  6725842.1875,
  6725887.5,
  6725910.9375,
  6726640.625,
  6726971.875,
  6727012.5,
  6727056.25,
  6727259.375,
  6727626.5625,
  6728512.5,
  6728743.75,
  6728762.5,
  6731779.6875,
  6732260.9375,
  6732432.8125,
  6732892.1875,
  6733714.0625,
  6734140.625,
  6734910.9375,
  6735528.125,
  6738829.6875,
  6738943.75,
  6739029.6875,
  6739410.9375,
  6739643.75,
  6740900.0,
  6742184.375,
  6742465.625,
  6743915.625,
  6744950.0,
  6745137.5,
  6746487.5,
  6746514.0625,
  6747059.375,
  6750815.625,
  6752242.1875,
  6752393.75,
  6752410.9375,
  6752435.9375,
  6753500.0,
  6754150.0,
  6758209.375,
  6759045.3125,
  6759801.5625,
  6760979.6875,
  6761115.625,
  6775076.5625,
  6775379.6875,
  6779660.9375,
  6782918.75,
  6788901.5625,
  6790290.625,
  6795151.5625,
  6805556.25,
  6807815.625,
  6822965.625,
  6823621.875,
  6823632.8125,
  6828728.125,
  6829656.25,
  6841743.75,
  6841768.75,
  6845446.875,
  6845564.0625,
  6846368.75,
  6846379.6875,
  6846431.25,
  6846573.4375,
  6846617.1875,
  6846870.3125,
  6847267.1875,
  6849175.0,
  6849739.0625,
  6850175.0,
  6850481.25,
  6852170.3125,
  6852385.9375,
  6868739.0625,
  6868917.1875,
  6882912.5,
  6888640.625,
  6889214.0625,
  6889606.25,
  6889843.75,
  6889950.0,
  6890810.9375,
  6891292.1875,
  6891485.9375,
  6892315.625,
  6892803.125,
  6892867.1875,
  6893407.8125,
  6895385.9375,
  6895560.9375,
  6896421.875,
  6897648.4375,
  6897671.875,
  6897778.125,
  6897896.875,
  6898139.0625,
  6898165.625,
  6898393.75,
  6898721.875,
  6898735.9375,
  6898770.3125,
  6898815.625,
  6898945.3125,
  6899265.625,
  6899350.0,
  6899509.375,
  6899554.6875,
  6899812.5,
  6900217.1875,
  6900429.6875,
  6900546.875,
  6902576.5625,
  6902656.25,
  6906939.0625,
  6913245.3125,
  6913442.1875,
  6913759.375,
  6913784.375,
  6914006.25,
  6923751.5625,
  6925965.625,
  6927595.3125,
  6932939.0625,
  6938412.5,
  6938500.0,
  6943160.9375,
  6945146.875,
  6950165.625,
  6951968.75,
  6957435.9375,
  6960814.0625,
  6964254.6875,
  6964367.1875,
  6967328.125,
  6967510.9375,
  6969062.5,
  6969065.625,
  6970078.125,
  6970951.5625,
  6970954.6875,
  6972007.8125,
  6972178.125,
  6972671.875,
  6973873.4375,
  6975665.625,
  6975689.0625,
  6975725.0,
  6976389.0625,
  6976860.9375,
  6977046.875,
  6977210.9375,
  6977351.5625,
  6978246.875,
  6978254.6875,
  6978870.3125,
  6979339.0625,
  6980842.1875,
  6980962.5,
  6983843.75,
  6984129.6875,
  6984187.5,
  6985410.9375,
  6985442.1875,
  6985471.875,
  6985501.5625,
  6985843.75,
  6986506.25,
  6987567.1875,
  6988475.0,
  6988809.375,
  6988904.6875,
  6991437.5,
  6991546.875,
  6992840.625,
  6993178.125,
  6993340.625,
  6994176.5625,
  6994184.375,
  6994189.0625,
  6994957.8125,
  6995317.1875,
  6995326.5625,
  6995410.9375,
  6995951.5625,
  6995968.75,
  6998239.0625,
  6998546.875,
  7000923.4375,
  7001118.75,
  7001329.6875,
  7002537.5,
  7002551.5625,
  7003082.8125,
  7003625.0,
  7003793.75,
  7003798.4375,
  7003840.625,
  7003860.9375,
  7004200.0,
  7005406.25,
  7007182.8125,
  7011151.5625,
  7011412.5,
  7011417.1875,
  7011759.375,
  7013165.625,
  7014240.625,
  7015301.5625,
  7018817.1875,
  7019067.1875,
  7021000.0,
  7022381.25,
  7022514.0625,
  7022531.25,
  7022943.75,
  7023767.1875,
  7024992.1875,
  7025025.0,
  7025112.5,
  7025460.9375,
  7025710.9375,
  7025939.0625,
  7025998.4375,
  7026079.6875,
  7026495.3125,
  7027435.9375,
  7027623.4375,
  7028498.4375,
  7028890.625,
  7029484.375,
  7030900.0,
  7031848.4375,
  7032140.625,
  7032292.1875,
  7032504.6875,
  7033356.25,
  7034859.375,
  7034962.5,
  7034968.75,
  7035606.25,
  7035607.8125,
  7036428.125,
  7037121.875,
  7037123.4375,
  7037503.125,
  7038662.5,
  7038715.625,
  7039207.8125,
  7039254.6875,
  7039350.0,
  7039354.6875,
  7039662.5,
  7039820.3125,
  7039956.25,
  7040417.1875,
  7040514.0625,
  7040625.0,
  7040729.6875,
  7041090.625,
  7041131.25,
  7041134.375,
  7041375.0,
  7041817.1875,
  7041831.25,
  7041843.75,
  7041859.375,
  7041875.0,
  7041893.75,
  7041990.625,
  7042067.1875,
  7042160.9375,
  7042171.875,
  7042221.875,
  7042579.6875,
  7042771.875,
  7042873.4375,
  7043117.1875,
  7043200.0,
  7043381.25,
  7043401.5625,
  7043826.5625,
  7044237.5,
  7044504.6875,
  7057339.0625,
  7058273.4375,
  7060925.0,
  7061534.375,
  7062500.0,
  7063012.5,
  7063031.25,
  7063771.875,
  7065518.75,
  7065803.125,
  7065828.125,
  7066556.25,
  7067959.375,
  7068065.625,
  7068767.1875,
  7070281.25,
  7074481.25,
  7074498.4375,
  7074615.625,
  7074687.5,
  7075543.75,
  7075842.1875,
  7075893.75,
  7076270.3125,
  7076410.9375,
  7076723.4375,
  7078187.5,
  7079053.125,
  7079450.0,
  7080343.75,
  7085562.5,
  7087193.75,
  7087200.0,
  7087267.1875,
  7087456.25,
  7091034.375,
  7097928.125,
  7098204.6875,
  7098815.625,
  7101859.375,
  7103120.3125,
  7108912.5,
  7108915.625,
  7110834.375,
  7117862.5,
  7121410.9375,
  7122609.375,
  7123043.75,
  7123385.9375,
  7123467.1875,
  7123743.75,
  7124987.5,
  7126746.875,
  7127312.5,
  7128600.0,
  7129276.5625,
  7129406.25,
  7130171.875,
  7130509.375,
  7131173.4375,
  7132093.75,
  7132645.3125,
  7134546.875,
  7137681.25,
  7138637.5,
  7142192.1875,
  7147798.4375,
  7148698.4375,
  7148828.125,
  7149118.75,
  7149551.5625,
  7151535.9375,
  7151878.125,
  7152326.5625,
  7152554.6875,
  7153806.25,
  7153823.4375,
  7155406.25,
  7156182.8125,
  7156943.75,
  7157062.5,
  7157703.125,
  7158659.375,
  7159415.625,
  7162153.125,
  7163025.0,
  7163293.75,
  7165400.0,
  7165432.8125,
  7165526.5625,
  7168846.875,
  7172923.4375,
  7174017.1875,
  7176732.8125,
  7178112.5,
  7181457.8125,
  7182885.9375,
  7184625.0,
  7186610.9375,
  7186668.75,
  7188978.125,
  7189812.5,
  7191117.1875,
  7193868.75,
  7195989.0625,
  7199860.9375,
  7201225.0,
  7202990.625,
  7204723.4375,
  7208831.25,
  7209640.625,
  7209660.9375,
  7210468.75,
  7210945.3125,
  7211884.375,
  7212140.625,
  7212409.375,
  7212729.6875,
  7213804.6875,
  7214389.0625,
  7214554.6875,
  7214753.125,
  7214789.0625,
  7214909.375,
  7215381.25,
  7215393.75,
  7215409.375,
  7215467.1875,
  7215560.9375,
  7215643.75,
  7215696.875,
  7215767.1875,
  7215785.9375,
  7215818.75,
  7215867.1875,
  7215906.25,
  7215917.1875,
  7215950.0,
  7216175.0,
  7216357.8125,
  7216450.0,
  7216520.3125,
  7216676.5625,
  7216940.625,
  7216970.3125,
  7217082.8125,
  7217226.5625,
  7217642.1875,
  7217714.0625,
  7218464.0625,
  7219112.5,
  7219385.9375,
  7221445.3125,
  7222723.4375,
  7231503.125,
  7243095.3125,
  7245731.25,
  7245950.0,
  7246095.3125,
  7248907.8125,
  7250520.3125,
  7258990.625,
  7260260.9375,
  7260695.3125,
  7261003.125,
  7261865.625,
  7266873.4375,
  7282703.125,
  7283145.3125,
  7284800.0,
  7286757.8125,
  7286831.25,
  7291665.625,
  7294006.25,
  7309860.9375,
  7312632.8125,
  7313651.5625,
  7314429.6875,
  7315393.75,
  7320562.5,
  7322117.1875,
  7328507.8125,
  7354607.8125,
  7354776.5625,
  7355726.5625,
  7356200.0,
  7356939.0625,
  7357431.25,
  7394865.625,
  7434064.0625,
  7576507.8125,
  7585914.0625,
  7620604.6875,
  7620907.8125,
  7624446.875,
  7653101.5625,
  7683304.6875,
  7684326.5625,
  7684378.125,
  7714845.3125,
  7736729.6875,
  7740359.375,
  7752598.4375,
  7753520.3125,
  7753707.8125,
  7778196.875,
  7781278.125,
  7781287.5,
  7781371.875,
  7879235.9375,
  8104754.6875,
  8173246.875,
  8539956.25,
  8551959.375,
  8552020.3125,
  8721107.8125,
  8761750.0,
  8872392.1875,
  8907807.8125,
  8908339.0625,
  8910898.4375,
  8916675.0,
  8920428.125,
  8931807.8125,
  9182890.625,
  9193826.5625,
  9195637.5,
  9261812.5,
  9275353.125,
  9286118.75,
  9290906.25,
  9294509.375,
  9318065.625,
  9352184.375,
  9356853.125,
  9359085.9375,
  9368718.75,
  9376306.25,
  9378681.25,
  9411829.6875,
  9454026.5625,
  9456243.75,
  9456651.5625,
  9456993.75,
  9457723.4375,
  9461206.25,
  9462717.1875,
  9463818.75,
  9474403.125,
  9482407.8125,
  9486654.6875,
  9487507.8125,
  9490128.125,
  9490478.125,
  9506351.5625,
  9507842.1875,
  9524360.9375,
  9554387.5,
  9581104.6875,
  9622239.0625,
  9628129.6875,
  9642587.5,
  9649184.375,
  9666600.0,
  9669510.9375,
  9669579.6875,
  9671387.5,
  9673656.25,
  9688442.1875,
  9711200.0,
  9711279.6875,
  9827278.125,
  9841485.9375,
  9889659.375,
  9893678.125,
  9894756.25,
  9894798.4375,
  9895709.375,
  9895735.9375,
  9897373.4375,
  9897423.4375,
  9936626.5625,
  9962996.875,
  9995064.0625,
  10013256.25,
  10026162.5,
  10026946.875,
  10029226.5625,
  10036209.375,
  10036742.1875,
  10042568.75,
  10042981.25,
  10046739.0625,
  10046743.75,
  10049373.4375,
  10087198.4375,
  10090228.125,
  10116943.75,
  10263607.8125,
  11257718.75,
  11267890.625,
  11267957.8125,
  11280456.25],
 [36.39723457256683,
  61.42011606542019,
  57.63331380236617,
  8.929290614776288,
  5.867756950569329,
  103.15167346425969,
  13.632253326347488,
  61.988275601110445,
  15.782696782507756,
  5.227109715887212,
  43.40145628763871,
  102.45907594437872,
  70.63837596907882,
  128.59382627377136,
  17.587650568578034,
  15.793167642101507,
  7.671813627472112,
  9.006379654396119,
  28.70609851781273,
  100.29981138667708,
  78.80329657045203,
  23.225821689737273,
  88.18321797392157,
  38.34357814715931,
  23.47776991851982,
  5.225485454867674,
  5.659364363078715,
  76.76677326861977,
  8.024261796596173,
  24.275962905745054,
  25.377084321643334,
  83.89258329860951,
  23.946462742652194,
  50.93071164949142,
  58.51405958943931,
  67.06249961993483,
  55.847138471574354,
  14.149140945397248,
  20.544470133004168,
  10.0180757997359,
  68.33173525422447,
  108.77033466229885,
  26.194573928714608,
  8.360969807334461,
  11.591584576740122,
  15.758665732968948,
  34.030938217134114,
  5.481389065289376,
  24.376277770638143,
  42.175379433676426,
  53.61122422065942,
  5.653454013029906,
  15.857975173415213,
  65.9539246648907,
  8.310939051713781,
  43.64728620668559,
  192.57406810768012,
  50.15962257049653,
  34.27502430628671,
  85.83646132075555,
  34.622424313148166,
  131.8763932836117,
  70.45556578208915,
  24.117761348660785,
  19.78736390407145,
  19.24625273551636,
  259.5303046208734,
  9.267275539627011,
  94.52168757313635,
  25.67429902289775,
  17.247923347038633,
  16.99510366372298,
  18.338551681601587,
  13.32147694962425,
  63.72091543156779,
  56.37430669301801,
  50.916425410211986,
  125.0897587093759,
  58.48166225709635,
  16.703579595561827,
  19.11011430135429,
  25.51252779403324,
  32.804139755043394,
  10.697281014364044,
  54.34203430499464,
  19.291934742150524,
  124.8974664064839,
  15.609415701052098,
  5.432567156193797,
  73.14207159300558,
  19.598637027356833,
  7.306629199403174,
  11.296177152251182,
  47.552568072581444,
  36.96907227087646,
  54.25984614013294,
  23.25976152933508,
  45.334650100167025,
  176.35599253956264,
  17.607075955173745,
  82.81061103496214,
  56.31303715647137,
  8.231195933068582,
  93.65401967694291,
  157.1143954299284,
  62.03569970000155,
  7.316580896670568,
  7.052124010949603,
  35.7759027797182,
  30.82360542216003,
  9.798540488645234,
  100.03988013196754,
  11.38314142977542,
  7.258916934182001,
  39.52951734617878,
  52.82393617615527,
  34.045747539411394,
  63.760109539502466,
  12.6494988393673,
  29.320364298043064,
  21.823657080944493,
  7.000983790584903,
  26.192950922073393,
  5.0758953453095605,
  14.418031387502166,
  33.6527390632858,
  76.86622242713834,
  13.311024654878873,
  29.35449753039608,
  65.95837123281866,
  12.413601516380288,
  29.296784189244242,
  97.81449871045096,
  99.83564237066203,
  40.2476894423131,
  51.4244131011592,
  79.48771264409444,
  43.938712458056315,
  13.011741221320914,
  12.867665300025267,
  44.411861212108604,
  5.698327666029522,
  23.220022487333814,
  7.677072247957201,
  123.11936950859335,
  77.70092946670269,
  51.24469015790008,
  206.5094388591854,
  22.04565298695637,
  48.40618971229868,
  29.481152250632725,
  75.84270141599688,
  107.3155784615033,
  20.330654925873958,
  101.113585193938,
  68.58657363330593,
  19.050999659011822,
  15.296627976197115,
  32.56451285716979,
  21.81137862581718,
  59.30814293658601,
  60.708985931199436,
  137.42236923465742,
  129.95393069598427,
  49.771180646118474,
  49.45136136964324,
  24.937731699140002,
  118.24452935403079,
  20.073764992819875,
  153.86990491667837,
  107.63172441621319,
  15.016307833826497,
  8.408690366765153,
  30.312961874645012,
  5.997958863078135,
  8.452608138764365,
  5.5570975373728,
  8.001892567757281,
  73.79057910086591,
  21.111881574438257,
  6.066153738741113,
  47.08419560585291,
  56.67510204822391,
  43.71788249107159,
  82.39327886377018,
  57.65663076026064,
  13.303113273263156,
  10.330633974957554,
  5.755159609170377,
  19.040784110635492,
  5.621669378083233,
  69.64824848732714,
  33.21099769977836,
  62.609628259637134,
  69.44656875655085,
  19.917878241918054,
  73.48391273876348,
  15.59866757654158,
  28.525795935709684,
  16.450191714611123,
  19.085978566616227,
  9.847199642537985,
  13.939621181470327,
  17.2774418111779,
  10.45611498110723,
  8.539391555424086,
  15.47772590263954,
  51.50271746735625,
  27.51102222060111,
  39.995722741166865,
  10.561223311222646,
  5.49784137905199,
  32.066418116179484,
  96.33037548367857,
  39.25960964222006,
  17.682133427314476,
  96.80524359126099,
  111.17831515119403,
  11.504910229982594,
  16.707816525480926,
  71.48760686188146,
  27.18776603881378,
  38.39960078758963,
  16.4254365726055,
  60.08746596578121,
  45.83385414665403,
  23.353312654207613,
  58.85832501460436,
  41.66846193034802,
  25.889133080981267,
  27.0149215014014,
  112.80447085716244,
  10.512909891136907,
  13.662628979492569,
  70.4005236572396,
  74.16656977228209,
  35.238825696528004,
  28.48184039740719,
  10.924899051597857,
  83.17739913908383,
  16.62831994178219,
  60.086958956855675,
  57.91383276906872,
  20.148785539171563,
  18.427659911220164,
  65.26106321836986,
  7.572559540440718,
  115.73928889501224,
  31.863727267981485,
  82.51041484203607,
  35.361293566290534,
  23.648882747165857,
  56.37135508295516,
  101.0228735666784,
  36.995774727011955,
  239.95573348165058,
  8.574192640784469,
  37.96156166995581,
  89.42040732347566,
  20.22137316037128,
  49.65146830305179,
  20.658802593847145,
  27.283433359698986,
  109.55191884710541,
  23.196824737527788,
  146.457331866237,
  10.2022433282344,
  37.952321428908974,
  13.13483637306454,
  9.545625899412434,
  13.385410919813104,
  38.842488194322364,
  17.783996143897205,
  14.750502476257601,
  6.112562033292657,
  18.771899827565584,
  24.80375480532396,
  78.04036275641678,
  46.72856783120531,
  5.295269109919355,
  93.29325981912626,
  99.98597977415547,
  5.966135850601805,
  14.261884710409893,
  16.048538295470127,
  74.57563138024784,
  5.3234613789131,
  34.57341420045836,
  5.3586324095810705,
  44.485182068559666,
  45.081714863343834,
  16.875670462561438,
  88.20729215095666,
  12.437618944272808,
  19.696139110195446,
  71.00349229313728,
  23.773735673476274,
  105.41442588999934,
  28.946077332056817,
  47.253734933968076,
  5.462575582795934,
  30.145835996026317,
  40.09501662518248,
  11.427284583854847,
  5.55686999448816,
  12.318646529880045,
  17.88578485995733,
  107.07422046782936,
  18.249795060960555,
  32.26400018416813,
  7.930158419839312,
  110.67809866341692,
  100.21268052669818,
  42.514685206195864,
  5.095016070893412,
  61.636161580364984,
  70.08682666651234,
  11.157789722497984,
  29.222816861491555,
  12.971300139419613,
  27.543324638141456,
  7.573531274032749,
  12.324898001981433,
  29.7710779204007,
  9.14804439315284,
  18.360638645955152,
  112.90374119668057,
  132.67254040535673,
  7.162134145727406,
  36.49372334823385,
  6.622016997883784,
  80.68135832342915,
  13.82305552678966,
  77.48743691819745,
  8.92290458525,
  14.580672675474293,
  45.27371879263896,
  7.602347372635561,
  46.92059890289319,
  23.19730413326542,
  7.855230077894405,
  15.048063946899774,
  32.50167566475677,
  16.400917327583734,
  18.337144223995534,
  40.4265125267047,
  6.59677491956983,
  39.655220661190405,
  40.04895538926286,
  56.05252352265611,
  15.331484932504182,
  13.099252738344084,
  22.280172788529264,
  70.29092211381443,
  47.92602721364865,
  195.334377726145,
  5.323626736826369,
  20.677837781591993,
  88.1626383636549,
  10.108983179239246,
  28.152152489077626,
  19.25695419449508,
  54.61985711274021,
  26.572763025421242,
  6.629374543258719,
  25.05963433413907,
  41.21229926008208,
  77.93295679007664,
  250.77590879265006,
  9.035192382103938,
  54.87175112220168,
  66.64111531820397,
  25.91709629781178,
  76.67381992941066,
  53.40916499705752,
  16.340801708160313,
  26.298788136451915,
  214.9276806945125,
  66.19947733960495,
  68.80751329046113,
  63.46934616769979,
  13.537669646613109,
  89.7327046603615,
  102.63149776891655,
  46.08703091630864,
  37.1145750733646,
  21.071852348393904,
  42.392824316576274,
  17.86161228472567,
  42.425304714793135,
  45.48455474609187,
  50.216542413045,
  5.360485862196153,
  19.175318322162177,
  11.012932792400386,
  72.29662011971712,
  135.01674338159373,
  61.074610629654856,
  60.62720243901956,
  15.357241098856736,
  7.759515874011704,
  5.989970420166393,
  5.024611979714108,
  9.060548530053877,
  38.32735034697739,
  62.119420008907156,
  24.17980582768046,
  13.694990792726543,
  57.4125608820768,
  9.893305914485136,
  81.69622814528651,
  21.74521986225286,
  8.750634094959812,
  29.14782057690727,
  63.97159745895686,
  7.130111006729035,
  59.81747671424412,
  12.539494478704636,
  7.291013191795108,
  26.3317224034226,
  72.43517609089575,
  116.95774507794263,
  91.63969363108122,
  15.227370288491633,
  35.72033358393782,
  27.30521793086246,
  5.210992339876186,
  22.58353781299215,
  64.3951221857838,
  13.805659203132972,
  22.311537867308818,
  8.16569461317244,
  107.99538676849032,
  17.077027999769637,
  47.66027591322897,
  41.37640024765623,
  42.07896026079423,
  39.4994745302631,
  69.97118573483006,
  55.399735251854224,
  15.92276263285325,
  17.585232347653264,
  51.09300150269347,
  141.83538047771194,
  25.860473614439737,
  51.70004967125205,
  138.24599936323148,
  6.296551660553312,
  27.066458099730465,
  66.53279747520995,
  81.67574145991111,
  115.38709229549974,
  18.80387966381403,
  43.13626442991365,
  117.90100940387464,
  193.27713109874847,
  106.92378733288527,
  28.925922881684542,
  42.98317903938009,
  104.82569972617829,
  18.09353498595276,
  32.274513881869055,
  45.67891945659834,
  94.57165448160939,
  47.54650021282325,
  14.895962407906907,
  26.768832498528994,
  17.551074736142564,
  30.699116105173232,
  6.935831046905322,
  115.22257517935395,
  7.469597110313175,
  87.76036775135202,
  12.25074032922134,
  41.11555506589667,
  23.83262040248162,
  18.002465826479334,
  23.81477842762854,
  8.838077822794379,
  21.68258229385632,
  28.00852373244167,
  90.11889493377922,
  163.07121689626797,
  20.540108270706625,
  65.80541814133274,
  47.786993726294426,
  29.335042397402365,
  93.86532569676896,
  99.56934917062974,
  17.757248337781288,
  55.39859283739911,
  66.25843115182607,
  118.28008005675935,
  30.671089556711923,
  5.051039471698173,
  70.7305224709122,
  119.90700956452221,
  230.26580581163404,
  25.971573098422343,
  12.149888203213914,
  5.703698023666907,
  152.02059815483406,
  66.6988898443032,
  52.15059682527077,
  28.950333905865026,
  7.908227402268309,
  184.35839680302357,
  117.90609263310223,
  46.24631401665481,
  14.498498191791379,
  54.689428424699976,
  20.35447639425396,
  61.38881776133482,
  38.96333242894729,
  25.431385167612866,
  84.12197391091124,
  12.184693790405046,
  6.327233705329153,
  31.548206682192713,
  15.271782150545015,
  5.238278021608867,
  57.11939878703496,
  97.18729218157577,
  57.077332620175994,
  15.686016685863741,
  17.702169318602504,
  66.30863660574207,
  27.550598380821416,
  163.15543480322236,
  150.88784637190432,
  30.871435262025717,
  91.10225921377655,
  25.816027637756207,
  78.26760028282013,
  26.755804301441145,
  12.876748788353112,
  25.8584221967763,
  24.449481349288803,
  30.709596758514653,
  231.14298469256718,
  18.970093123992815,
  48.1997996113733,
  25.00684685505102,
  29.385623027793265,
  44.436833822308195,
  14.250795298761048,
  207.57083810619093,
  130.15807504395414,
  119.25182099085103,
  24.66519540315926,
  101.99071395184572,
  12.308830987562384,
  57.82716184538903,
  9.731976670841103,
  129.23115679400652,
  16.72269343241506,
  5.423619487829046,
  70.55877603142835,
  19.425067156256375,
  25.788211884679193,
  21.750981785451568,
  172.121249594233,
  68.34847770250651,
  56.21839414178929,
  25.765890989278137,
  78.6671097282135,
  19.638280989855588,
  23.331589709919193,
  30.161396048893327,
  7.708445399624251,
  46.14224247816008,
  5.540215973232311,
  40.7158900641637,
  41.11487049581796,
  40.21710777591162,
  134.14653994173736,
  17.504486033334462,
  52.23380160910488,
  80.08579439489925,
  5.922549508149697,
  19.69870709776275,
  10.377037135534696,
  24.217721677720853,
  70.85770556994306,
  63.97574868434822,
  53.62949064490216,
  15.11559126517665,
  51.747009104773994,
  81.054476769162,
  43.86776692651202,
  202.6870560424551,
  18.108006834507794,
  7.7722192158968655,
  6.105094811494561,
  5.373087842024246,
  14.050546904503438,
  32.859317000322775,
  18.099902704058707,
  121.00330014735792,
  45.78391145108671,
  99.5478504324995,
  102.88623227747873,
  116.5460079068383,
  27.886941051564182,
  66.91752860174782,
  5.080617267890299,
  108.49749951238559,
  87.86397757255472,
  7.114315597064646,
  21.569724714252455,
  23.42700522992944,
  26.796810870875007,
  25.071296640276515,
  132.63278955868236,
  8.412439098701245,
  28.48464565877721,
  36.88131302685445,
  67.218326194331,
  5.5182498282627535,
  43.827979250495424,
  59.32039552962836,
  71.31527646754795,
  35.93046759540154,
  12.59109847741283,
  24.788957921633937,
  129.3571454260518,
  11.887722194011053,
  40.99694814611819,
  94.71308968063315,
  36.73503457778565,
  13.97215306463258,
  5.14119302486334,
  163.0448576896018,
  86.75943414643599,
  16.038040980769715,
  15.664011593266858,
  33.00629934178887,
  148.83283057136592,
  110.98418488934249,
  100.3834860382854,
  15.75717528846002,
  12.307760569011332,
  53.475879372101176,
  16.568942792146007,
  99.58280660718562,
  185.49261717894046,
  76.75495121468843,
  8.804101245994955,
  12.072985631697277,
  41.06437663797372,
  15.778021398277152,
  24.57910530495824,
  9.708315924291124,
  21.649317940346812,
  11.242234913006676,
  7.466023442864695,
  90.37652339293187,
  29.41549523078106,
  35.862949981336186,
  54.91764088297558,
  60.21665190189594,
  12.016740812354248,
  36.34313114773403,
  22.91076117181387,
  9.464872222184768,
  83.44011472879953,
  118.60474671374632,
  16.90482374613588,
  104.61256431685752,
  72.95280664984845,
  6.268585572900139,
  20.713233602278045,
  80.39620366502926,
  52.181210928792396,
  6.938050559871266,
  20.648118035248583,
  30.99681502206858,
  49.628281036959365,
  9.584193167767525,
  16.509802291189768,
  125.15645710139621,
  25.344732339623746,
  15.977223526087242,
  14.286513788791183,
  96.5675655045218,
  58.23532270886482,
  47.68076013957206,
  16.018604745092354,
  26.054669212727973,
  6.652869136204719,
  54.29843236753274,
  14.414812044652408,
  16.639164896276604,
  71.41174692229937,
  20.833132629259005,
  68.45718923072876,
  33.26616913682368,
  20.047345746189144,
  27.934179038482363,
  19.299217675171178,
  107.58302366828818,
  75.6500306749686,
  6.651848240630262,
  31.944831622204124,
  73.3765233425237,
  39.67320774469785,
  17.440999100494246,
  34.44164294826572,
  126.93734537630989,
  16.620727697732978,
  69.54665555127444,
  60.45629049673758,
  5.63105086912868,
  19.522051187072236,
  62.27096149541468,
  60.67195590491417,
  86.52526589928303,
  38.61207726252993,
  7.532671805134615,
  50.48846472674089,
  19.07056108463937,
  28.420606500143194,
  36.757386983235634,
  10.923075603192402,
  5.3225035018966596,
  33.21727700725195,
  54.134521734610104,
  12.75159216306036,
  128.43322641248153,
  54.06743389099918,
  101.81941678298769,
  5.101439543090858,
  16.94303162038046,
  27.3345310750813,
  10.484055949992534,
  29.616859538365,
  148.77454016864962,
  62.74704913622194,
  97.1840066969832,
  44.560913847901304,
  5.589348195126645,
  16.590903496044287,
  30.462369070591823,
  13.418757577771217,
  30.75327493158142,
  15.888223149917332,
  46.660505016633586,
  45.19916506197506,
  22.215699891678508,
  117.4328234780254,
  30.18347214809416,
  104.71470924944434,
  5.984992075253884,
  48.35175979250335,
  46.84055650898214,
  58.677855926794905,
  38.98927602372218,
  93.2079691839044,
  85.63224558858346,
  5.888573005613501,
  12.593361415831872,
  87.74769253283105,
  19.710567312575776,
  26.34004550831608,
  13.495745703797754,
  27.70258630395062,
  13.568252275120846,
  158.12285374798728,
  39.73804158175106,
  29.626687211705637,
  5.914663711753748,
  75.04314916062881,
  9.315255833188704,
  43.87834060130922,
  31.44055115826913,
  12.265292996886549,
  22.490785312920128,
  43.63774794211784,
  11.71939375103589,
  8.262439396794834,
  9.105780844777515,
  46.50149903319378,
  15.887104693677763,
  63.99154882662474,
  101.09270443189122,
  29.939639858062307,
  19.507409018022656,
  74.69168717780781,
  28.089010654182335,
  23.25564009212045,
  5.085162114670139,
  81.78955290350748,
  42.66747885033354,
  106.00899290103145,
  29.715198934034966,
  67.27489565399378,
  48.80267760135854,
  34.16138395616179,
  18.007079459133735,
  21.84921166905134,
  38.96767810201118,
  80.75319523370256,
  13.347909526829145,
  24.594326419619406,
  10.561663315977656,
  6.624031667018189,
  75.94491658220176,
  18.79511353010632,
  39.79939997999378,
  33.36096152171195,
  101.36317860277666,
  113.76091662574385,
  24.640488746882934,
  86.66010868403646,
  5.706687226396635,
  44.32853366457902,
  5.51405171921014,
  108.70895805014075,
  5.5310254389849804,
  70.41057491272504,
  47.12389223818455,
  108.6989571892436,
  182.4737320633952,
  29.998573992380294,
  49.69432611728146,
  5.029258328330005,
  32.304414962535944,
  17.36997124214708,
  121.19064655191032,
  101.87795678139386,
  45.04779284734584,
  48.593768426363894,
  35.64057058169847,
  107.42594130365052,
  31.45399264785776,
  21.200279346067923,
  8.461430234884595,
  24.283934472678013,
  82.48528117232584,
  32.821648162731165,
  57.605965967311306,
  43.30185708596995,
  45.644495911362725,
  17.320360771068028,
  183.16915853507663,
  8.520177703519378,
  18.312309167844806,
  121.68491974258711,
  42.862973821472316,
  47.50439172519498,
  99.57952092077544,
  54.71425161718154,
  34.623437296487715,
  20.865040033802842,
  11.440674212609661,
  17.025758523183047,
  30.839824659911425,
  90.49740866544097,
  27.157401723102364,
  86.66107886847195,
  33.56837965661565,
  54.88772738836616,
  50.96950534021002,
  24.65996221703349,
  37.96578564009686,
  63.34338405086122,
  43.34534683533772,
  13.98648715381339,
  32.32862240491876,
  94.97892029322108,
  16.403567208079387,
  85.04221522791974,
  77.79832017093501,
  6.5854865470334225,
  10.807054428923905,
  81.93777765066649,
  6.976146958220431,
  46.228462547775315,
  38.6641998882411,
  17.736808397770282,
  101.74583667175969,
  59.2889039077969,
  45.20372525382959,
  127.32025166287379,
  63.21859304486967,
  31.281834140319642,
  12.078817022313709,
  19.263869615970876,
  32.991081530648486,
  49.64258354205053,
  8.374601105832745,
  7.029297240494546,
  16.558432035663078,
  53.82449793540092,
  43.566399940304734,
  91.67528624738806,
  9.157749983205344,
  6.579909981910118,
  17.56094362280539,
  40.265505747996045,
  56.573520100534665,
  17.413538214154975,
  12.7797434429684,
  32.35867425774145,
  77.82026538818306,
  5.121533518980227,
  101.7578709870997,
  77.48368832063628,
  55.202222263927275,
  18.626238278642226,
  17.76215976406651,
  6.34910612362186,
  100.8895039327821,
  25.926771694425145,
  17.896138230928678,
  32.41485772983219,
  25.21402044751167,
  19.29398340171033,
  10.248217214793396,
  9.791959674920195,
  7.343725845564252,
  52.15446129367068,
  24.959982289937354,
  84.07555915953853,
  17.0113708133652,
  9.03658189124093,
  61.84057891527404,
  9.538604643810766,
  19.05735813538856,
  18.052449654779608,
  134.6959197010295,
  36.04259739875368,
  5.615267558339154,
  73.24430616310883,
  52.63145940532872,
  53.47126018274981,
  32.51113391013275,
  32.22240102863999,
  86.10351476706427,
  25.389287956549097,
  142.8491652309866,
  53.103076985387894,
  37.22687275226063,
  31.21248237943917,
  19.83284630745572,
  53.417590569020206,
  12.524809533106545,
  173.8357079722619,
  25.77554049379744,
  110.00950238043289,
  166.06392545198864,
  5.288620049276263,
  35.217247713507675,
  10.946266471633303,
  8.534230151249488,
  122.96767232576417,
  11.847801286111878,
  108.28968837739178,
  67.93743992899132,
  44.700093264336886,
  18.937412132932778,
  5.242587642431672,
  84.0679105394437,
  31.62858994253959,
  19.460901129165595,
  18.946922394244822,
  190.24575871275707,
  76.78417289047421,
  59.02407556178068,
  32.69123832598587,
  69.3843738199064,
  7.305782991260869,
  14.679698685690322,
  51.32037673368312,
  8.769510196013407,
  6.0803622344308685,
  41.08523329336502,
  6.117089677843529,
  45.86859084740329,
  181.14611119106755,
  29.55674198479838,
  37.99132297661872,
  25.980361459516537,
  12.50473661248837,
  52.30442459901546,
  7.494578440218566,
  74.4312715877626,
  35.26890893853241,
  54.144150710916286,
  6.435521839743759,
  18.174198918727328,
  54.50132881917325,
  6.512541534962708,
  73.1355307334673])

Graph for home page

Detected energies during the discharge + Energy spectrum

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