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

([2785207.8125,
  3082431.25,
  3120943.75,
  3157382.8125,
  3202934.375,
  3272376.5625,
  3316748.4375,
  3388651.5625,
  3402350.0,
  3402814.0625,
  3462937.5,
  3538078.125,
  3541453.125,
  3546592.1875,
  3978398.4375,
  5193150.0,
  5193853.125,
  5193853.125,
  5288400.0,
  5290445.3125,
  5290582.8125,
  5291470.3125,
  5291560.9375,
  5518459.375,
  5521100.0,
  5522295.3125,
  5527396.875,
  5549023.4375,
  5589784.375,
  5674087.5,
  5734862.5,
  5735451.5625,
  5736109.375,
  5766712.5,
  5768162.5,
  5770014.0625,
  5770028.125,
  5771626.5625,
  5773115.625,
  5774379.6875,
  5775762.5,
  5777510.9375,
  5777618.75,
  5778564.0625,
  5781085.9375,
  5811298.4375,
  5815667.1875,
  5820326.5625,
  5832665.625,
  5836910.9375,
  5847814.0625,
  5853029.6875,
  5855212.5,
  5857195.3125,
  5861312.5,
  5896285.9375,
  5926950.0,
  5943704.6875,
  5944151.5625,
  5947159.375,
  5947201.5625,
  5948203.125,
  5951895.3125,
  5953640.625,
  5954620.3125,
  5957484.375,
  5984882.8125,
  6002948.4375,
  6007176.5625,
  6008689.0625,
  6009878.125,
  6013659.375,
  6016290.625,
  6026440.625,
  6052370.3125,
  6052632.8125,
  6054620.3125,
  6059753.125,
  6065460.9375,
  6067029.6875,
  6068957.8125,
  6082809.375,
  6082829.6875,
  6082945.3125,
  6088596.875,
  6107103.125,
  6109231.25,
  6109357.8125,
  6109453.125,
  6109592.1875,
  6109603.125,
  6109667.1875,
  6110323.4375,
  6112467.1875,
  6113878.125,
  6114006.25,
  6114971.875,
  6115895.3125,
  6120559.375,
  6121596.875,
  6122051.5625,
  6137234.375,
  6157212.5,
  6170460.9375,
  6184042.1875,
  6203723.4375,
  6226789.0625,
  6226837.5,
  6244317.1875,
  6245564.0625,
  6253967.1875,
  6254454.6875,
  6254785.9375,
  6256253.125,
  6256690.625,
  6257285.9375,
  6257748.4375,
  6258062.5,
  6258906.25,
  6259164.0625,
  6259389.0625,
  6259390.625,
  6259434.375,
  6259456.25,
  6259920.3125,
  6259973.4375,
  6261121.875,
  6261453.125,
  6262770.3125,
  6265626.5625,
  6266168.75,
  6266823.4375,
  6267245.3125,
  6268345.3125,
  6268496.875,
  6270157.8125,
  6270442.1875,
  6271295.3125,
  6272292.1875,
  6272906.25,
  6273165.625,
  6276104.6875,
  6276667.1875,
  6278050.0,
  6278053.125,
  6280064.0625,
  6280557.8125,
  6286060.9375,
  6287243.75,
  6287435.9375,
  6289389.0625,
  6290634.375,
  6296525.0,
  6310990.625,
  6321596.875,
  6334381.25,
  6338228.125,
  6341671.875,
  6342565.625,
  6352170.3125,
  6369259.375,
  6411226.5625,
  6413437.5,
  6413898.4375,
  6416093.75,
  6423415.625,
  6423535.9375,
  6424820.3125,
  6429246.875,
  6436790.625,
  6454834.375,
  6457175.0,
  6462757.8125,
  6462843.75,
  6466093.75,
  6467762.5,
  6475650.0,
  6519840.625,
  6553476.5625,
  6561409.375,
  6562390.625,
  6581778.125,
  6638431.25,
  6757337.5,
  7017410.9375,
  7256976.5625,
  8369315.625,
  8591932.8125,
  9111762.5,
  9183620.3125,
  9436548.4375,
  9449868.75,
  9467976.5625,
  9517912.5,
  9518856.25,
  9525717.1875,
  9561678.125,
  9580257.8125,
  9594206.25,
  9594498.4375,
  9600270.3125,
  9627248.4375,
  9629625.0,
  9629626.5625,
  9630092.1875,
  9638162.5,
  9783900.0,
  9784475.0,
  9787135.9375,
  9787645.3125,
  9819854.6875,
  9823010.9375,
  9825565.625,
  9825842.1875,
  9827115.625,
  9827421.875,
  9827426.5625,
  9836901.5625,
  9836909.375,
  9838504.6875,
  9839989.0625,
  9842937.5,
  9843367.1875,
  9844545.3125,
  9845146.875,
  9845664.0625,
  9845945.3125,
  9866448.4375,
  9877801.5625,
  9879606.25,
  9879992.1875,
  9880634.375,
  9880693.75,
  9888973.4375,
  9904965.625,
  9966787.5,
  10057748.4375,
  10064606.25,
  10065767.1875,
  10066778.125,
  10089332.8125,
  10090520.3125,
  10097212.5,
  10110789.0625,
  10203009.375,
  10223107.8125,
  10237484.375,
  10253071.875,
  10261001.5625,
  10545165.625,
  10644046.875,
  10653407.8125,
  10852626.5625,
  10896526.5625,
  10910004.6875,
  10910665.625,
  10913295.3125,
  10916296.875,
  10916959.375,
  10921948.4375,
  10934098.4375,
  10934223.4375,
  10937068.75,
  10938762.5,
  10961723.4375,
  11055554.6875,
  11061625.0,
  11061629.6875,
  11065234.375,
  11065939.0625,
  11068700.0,
  11068701.5625,
  11071746.875,
  11072135.9375,
  11073032.8125,
  11073306.25,
  11073357.8125,
  11073520.3125,
  11075271.875,
  11075721.875,
  11076321.875,
  11076446.875,
  11076973.4375,
  11077978.125,
  11078221.875,
  11078734.375,
  11079070.3125,
  11079393.75,
  11079415.625,
  11079607.8125,
  11080146.875,
  11080193.75,
  11080298.4375,
  11080310.9375,
  11080645.3125,
  11080929.6875,
  11080971.875,
  11081698.4375,
  11081718.75,
  11082050.0,
  11082281.25,
  11082307.8125,
  11082323.4375,
  11082353.125,
  11082643.75,
  11082701.5625,
  11082745.3125,
  11082853.125,
  11083339.0625,
  11083515.625,
  11083525.0,
  11083562.5,
  11083687.5,
  11083854.6875,
  11083906.25,
  11084318.75,
  11084345.3125,
  11084723.4375,
  11085229.6875,
  11085325.0,
  11085395.3125,
  11085893.75,
  11086029.6875,
  11086185.9375,
  11086209.375,
  11086237.5,
  11087143.75,
  11087195.3125,
  11087578.125,
  11087856.25,
  11087867.1875,
  11088007.8125,
  11088218.75,
  11088292.1875,
  11088304.6875,
  11088304.6875,
  11088334.375,
  11088343.75,
  11088398.4375,
  11088425.0,
  11088487.5,
  11088784.375,
  11088882.8125,
  11089237.5,
  11089368.75,
  11089403.125,
  11089643.75,
  11089903.125,
  11090220.3125,
  11090796.875,
  11090929.6875,
  11090968.75,
  11091234.375,
  11091460.9375,
  11091546.875,
  11091678.125,
  11092253.125,
  11092732.8125,
  11092785.9375,
  11093701.5625,
  11094317.1875,
  11094365.625,
  11094393.75,
  11095151.5625,
  11095529.6875,
  11097075.0,
  11097107.8125,
  11097926.5625,
  11098078.125,
  11099054.6875,
  11099410.9375,
  11099456.25,
  11099535.9375,
  11099950.0,
  11100032.8125,
  11100079.6875,
  11100256.25,
  11100276.5625,
  11101029.6875,
  11101259.375,
  11101271.875,
  11101281.25,
  11101484.375,
  11101603.125,
  11101653.125,
  11101803.125,
  11102021.875,
  11102029.6875,
  11102059.375,
  11102203.125,
  11102985.9375,
  11103226.5625,
  11103518.75,
  11103614.0625,
  11103675.0,
  11103842.1875,
  11104142.1875,
  11104189.0625,
  11104301.5625,
  11104323.4375,
  11104648.4375,
  11104776.5625,
  11104896.875,
  11105165.625,
  11105376.5625,
  11105423.4375,
  11105507.8125,
  11105529.6875,
  11105667.1875,
  11105717.1875,
  11105789.0625,
  11105839.0625,
  11105862.5,
  11105876.5625,
  11106178.125,
  11106242.1875,
  11106457.8125,
  11106540.625,
  11106589.0625,
  11106614.0625,
  11106667.1875,
  11106676.5625,
  11106762.5,
  11107221.875,
  11107334.375,
  11107454.6875,
  11107496.875,
  11107514.0625,
  11107515.625,
  11107596.875,
  11107668.75,
  11107745.3125,
  11107826.5625,
  11107850.0,
  11108053.125,
  11108179.6875,
  11108443.75,
  11108637.5,
  11108671.875,
  11108693.75,
  11108790.625,
  11108793.75,
  11108915.625,
  11108946.875,
  11108957.8125,
  11108987.5,
  11109128.125,
  11109173.4375,
  11109184.375,
  11109207.8125,
  11109217.1875,
  11109218.75,
  11109468.75,
  11109476.5625,
  11109526.5625,
  11109528.125,
  11109532.8125,
  11109546.875,
  11109589.0625,
  11109648.4375,
  11109682.8125,
  11109762.5,
  11109764.0625,
  11109815.625,
  11109823.4375,
  11109900.0,
  11109923.4375,
  11109978.125,
  11110109.375,
  11110126.5625,
  11110157.8125,
  11110182.8125,
  11110335.9375,
  11110478.125,
  11110481.25,
  11110584.375,
  11110587.5,
  11110662.5,
  11110773.4375,
  11110832.8125,
  11110976.5625,
  11110981.25,
  11111006.25,
  11111042.1875,
  11111175.0,
  11111257.8125,
  11111454.6875,
  11111690.625,
  11111814.0625,
  11111854.6875,
  11111932.8125,
  11112251.5625,
  11112401.5625,
  11112442.1875,
  11112443.75,
  11112501.5625,
  11112504.6875,
  11112654.6875,
  11112684.375,
  11112864.0625,
  11112985.9375,
  11112992.1875,
  11113348.4375,
  11113796.875,
  11114065.625,
  11114189.0625,
  11114834.375,
  11114839.0625,
  11115062.5,
  11115182.8125,
  11115473.4375,
  11115515.625,
  11115557.8125,
  11115837.5,
  11115998.4375,
  11116151.5625,
  11116334.375,
  11116417.1875,
  11116784.375,
  11117264.0625,
  11117395.3125,
  11117406.25,
  11117420.3125,
  11117450.0,
  11118028.125,
  11118385.9375,
  11118392.1875,
  11118506.25,
  11118564.0625,
  11118782.8125,
  11118900.0,
  11118932.8125,
  11119284.375,
  11119553.125,
  11119553.125,
  11119893.75,
  11120350.0,
  11120417.1875,
  11120482.8125,
  11120735.9375,
  11120801.5625,
  11120812.5,
  11120928.125,
  11120996.875,
  11121046.875,
  11121134.375,
  11121464.0625,
  11121756.25,
  11121914.0625,
  11122115.625,
  11122132.8125,
  11122139.0625,
  11122153.125,
  11122168.75,
  11122171.875,
  11122265.625,
  11122351.5625,
  11122634.375,
  11122896.875,
  11123429.6875,
  11123431.25,
  11123457.8125,
  11123589.0625,
  11123618.75,
  11123646.875,
  11123748.4375,
  11123879.6875,
  11123895.3125,
  11124029.6875,
  11124148.4375,
  11124975.0,
  11125065.625,
  11125103.125,
  11125273.4375,
  11125634.375,
  11125673.4375,
  11125709.375,
  11125756.25,
  11125792.1875,
  11126028.125,
  11126160.9375,
  11126193.75,
  11126301.5625,
  11126320.3125,
  11126664.0625,
  11126760.9375,
  11126865.625,
  11126873.4375,
  11127034.375,
  11127121.875,
  11127193.75,
  11127804.6875,
  11128012.5,
  11128984.375,
  11129028.125,
  11129656.25,
  11129875.0,
  11129925.0,
  11129929.6875,
  11129989.0625,
  11130218.75,
  11130428.125,
  11130431.25,
  11130459.375,
  11130676.5625,
  11130693.75,
  11130704.6875,
  11130859.375,
  11130937.5,
  11130945.3125,
  11131132.8125,
  11131207.8125,
  11131667.1875,
  11131887.5,
  11132007.8125,
  11132314.0625,
  11132329.6875,
  11132342.1875,
  11132429.6875,
  11132759.375,
  11132781.25,
  11132870.3125,
  11133593.75,
  11133628.125,
  11133837.5,
  11133950.0,
  11134250.0,
  11134343.75,
  11134365.625,
  11134392.1875,
  11134403.125,
  11134665.625,
  11134707.8125,
  11134935.9375,
  11135507.8125,
  11135540.625,
  11135845.3125,
  11136234.375,
  11136250.0,
  11136298.4375,
  11136492.1875,
  11136496.875,
  11136520.3125,
  11136692.1875,
  11136929.6875,
  11137028.125,
  11137075.0,
  11137279.6875,
  11137348.4375,
  11137368.75,
  11137387.5,
  11137432.8125,
  11137906.25,
  11137917.1875,
  11138284.375,
  11138446.875,
  11138593.75,
  11138715.625,
  11138945.3125,
  11139017.1875,
  11139164.0625,
  11139339.0625,
  11139967.1875,
  11140118.75,
  11140201.5625,
  11140876.5625,
  11140879.6875,
  11140982.8125,
  11141414.0625,
  11141537.5,
  11141678.125,
  11141781.25,
  11141795.3125,
  11141985.9375,
  11142148.4375,
  11142553.125,
  11143376.5625,
  11143431.25,
  11143962.5,
  11144178.125,
  11144229.6875,
  11144285.9375,
  11144418.75,
  11144917.1875,
  11145003.125,
  11145079.6875,
  11145237.5,
  11145259.375,
  11145731.25,
  11146504.6875,
  11147417.1875,
  11147837.5,
  11148417.1875,
  11148946.875,
  11149095.3125,
  11149917.1875,
  11150837.5,
  11152126.5625,
  11152293.75,
  11152412.5,
  11152615.625,
  11152701.5625,
  11153057.8125,
  11153790.625,
  11154282.8125,
  11154315.625,
  11154379.6875,
  11154762.5,
  11154815.625,
  11155117.1875,
  11155557.8125,
  11156307.8125,
  11156856.25,
  11156921.875,
  11157521.875,
  11157682.8125,
  11158201.5625,
  11158256.25,
  11158400.0,
  11158575.0,
  11158640.625,
  11158940.625,
  11159389.0625,
  11160835.9375,
  11162770.3125,
  11162832.8125,
  11162939.0625,
  11163557.8125,
  11163739.0625,
  11163743.75,
  11163989.0625,
  11164571.875,
  11165528.125,
  11165548.4375,
  11166264.0625,
  11167178.125,
  11167421.875,
  11168056.25,
  11168268.75,
  11168671.875,
  11169267.1875,
  11169307.8125,
  11169353.125,
  11169476.5625,
  11169884.375,
  11170520.3125,
  11170590.625,
  11170656.25,
  11172698.4375,
  11173348.4375,
  11174301.5625,
  11176385.9375,
  11176690.625,
  11179804.6875,
  11180020.3125,
  11180125.0,
  11181593.75,
  11190251.5625],
 [5.152539720628002,
  6.854280608202352,
  28.46977324915323,
  8.518592729869086,
  104.12549718709128,
  10.851741588919381,
  5.517981870959336,
  66.07082341469258,
  5.07205040950528,
  70.83343627793475,
  13.477841232271121,
  21.514784216537013,
  7.294620685653026,
  71.59663946091045,
  19.04167810178794,
  101.8615281647461,
  38.31463262727778,
  15.748185106962069,
  22.337371770480527,
  41.86135388760191,
  24.76487392530209,
  42.20275329711578,
  22.83905188675709,
  64.1082339267195,
  26.240131209688553,
  59.809587418109714,
  44.944165232861174,
  39.24467288238996,
  90.35129801094376,
  30.49577284348799,
  120.38109198228508,
  11.191893569629812,
  60.31832820016083,
  28.40224654850701,
  11.554734225267199,
  40.18066455224005,
  36.20438031671707,
  28.351095269272477,
  57.1538165009531,
  51.573978550881286,
  87.88166151492781,
  36.69209596588307,
  13.675773006736271,
  28.839588101807237,
  55.04759646551981,
  5.687725581577804,
  20.92132427970974,
  61.755145281881425,
  45.54590494190182,
  66.63260202559736,
  18.41951993227862,
  9.361202130713966,
  41.85733994915828,
  40.41770608476546,
  77.83344741274658,
  33.50197557675041,
  15.768944741324155,
  67.71139077171931,
  84.67955027596174,
  6.011356114618717,
  6.125334205525346,
  92.5483191641962,
  135.5897747452631,
  5.960620084817627,
  40.39701737275518,
  10.540728353701452,
  39.45172103609671,
  76.39916273657171,
  63.80271822176243,
  5.790528280843724,
  7.200659289731451,
  11.509572524544506,
  80.01810437947763,
  5.171117644277947,
  96.27683852930934,
  186.1417397460132,
  8.458468574236372,
  13.888748650162718,
  28.979197039357043,
  14.203136917788044,
  19.147186778991607,
  20.71329242183629,
  25.29131804774221,
  5.413523852409813,
  130.16779268193548,
  34.235723812894264,
  7.014783924191486,
  67.48605592880169,
  171.36108906068372,
  34.293081690683756,
  25.59189536940345,
  5.875564676840401,
  130.5400459512419,
  180.8139704964108,
  24.40088221057529,
  30.664071638417916,
  49.29916302104552,
  19.211923297463933,
  38.35488307977782,
  7.194868237328745,
  23.600278627113184,
  10.923813758005368,
  47.607297203554104,
  13.970180127193217,
  5.32309791551664,
  109.8117977837779,
  18.285130042481654,
  5.542866385732408,
  43.80287967358354,
  14.510289209141238,
  22.583907452686088,
  24.847664972428685,
  33.82653267657091,
  21.16403345986571,
  171.12548785117116,
  116.9795915399358,
  66.17091062104805,
  9.443723749831062,
  12.087585797711258,
  31.973466546309027,
  16.435148161313535,
  41.212228278081646,
  15.900425747520845,
  20.733557048683927,
  25.197709162988858,
  65.82600010526188,
  20.129005151250173,
  11.749626610644865,
  139.47939737716968,
  68.48859484153373,
  24.977419719971078,
  14.152626115950188,
  14.780748339516347,
  25.40293149882287,
  76.51991348766333,
  17.54669692670169,
  19.305919145468287,
  60.40740473952349,
  7.2109320187492,
  39.20246516494882,
  48.269934970371956,
  70.74524123177379,
  56.14594605387755,
  11.086447050420286,
  45.2014525582852,
  25.559896113786696,
  25.11342373677851,
  17.169847986045625,
  79.73208922562054,
  52.590421189295434,
  36.095362761543846,
  5.346624621020493,
  42.44837923419093,
  87.28161169490454,
  25.881576481000494,
  8.137587074235956,
  53.190537051731575,
  17.887233167459463,
  31.59946636889701,
  24.074466485893478,
  19.496495567196902,
  26.346880458735214,
  15.315357157168105,
  49.99310863578063,
  80.9253123742679,
  9.473031339365475,
  22.486791936312983,
  11.662643518001484,
  71.85230768789084,
  16.735752443184825,
  103.20065404595042,
  7.586015517187785,
  20.600063522329673,
  5.960674491227817,
  85.55176226011383,
  75.01433825157935,
  20.08003280395637,
  12.648045863603514,
  9.51613277849609,
  7.93866252653696,
  86.82205887681428,
  72.23923501461027,
  10.401945538605132,
  7.95739135679953,
  16.391244436923426,
  40.535221212822655,
  34.272901388006744,
  54.51319317104874,
  10.935181603713145,
  41.68904159387064,
  36.64664766442374,
  60.270084587438525,
  51.54528200840906,
  41.04865140503001,
  18.521115835466233,
  77.97960074148779,
  42.40383516276643,
  8.043868596005852,
  18.45246269732988,
  49.958918820606584,
  56.28467471045701,
  81.4668970979767,
  25.588125516294543,
  88.46723420294437,
  49.20938784171525,
  29.492253688002982,
  26.113114805826186,
  14.513556031106724,
  46.422392689496654,
  10.1746537804279,
  11.483746364289411,
  7.318606183786745,
  6.346254242713159,
  22.096643556250786,
  23.224134371512825,
  44.04413327415761,
  18.522602832307346,
  29.05964972547983,
  22.509880825783704,
  7.700183749296881,
  65.02246401930786,
  68.74819988327873,
  43.97472934340662,
  108.84548564690715,
  7.320219495069112,
  6.2761340964351415,
  13.336095984452928,
  55.226124766388004,
  45.175026213035906,
  114.73517871425935,
  15.995647993758702,
  6.930998450410038,
  5.031119491253183,
  7.396871343343425,
  6.767683151708588,
  155.31080594303745,
  27.818023408680585,
  21.12701594558999,
  29.66940350703561,
  5.137249344001191,
  50.33604985175072,
  71.5731356473262,
  18.301471469860463,
  20.948050266697503,
  100.21797413158507,
  31.743090192893078,
  5.67603872794253,
  54.764653430482234,
  70.81253963895017,
  28.084685354651697,
  86.83339496681926,
  20.704348569831264,
  106.73138931339427,
  8.623631646227464,
  57.041228718383074,
  11.20262373289785,
  43.64437065664903,
  82.90455393049896,
  25.242442820201994,
  16.884654122020006,
  84.29590427404828,
  15.943988677524754,
  123.91754786059644,
  20.105722669737105,
  23.761109782440226,
  5.153674506095576,
  34.220916806604734,
  12.41614136873189,
  5.473770256582496,
  10.890228747204613,
  28.15591359082218,
  8.062699944680636,
  20.74396740876828,
  128.9700307914263,
  58.59504639904654,
  48.804394379641764,
  21.95333066002647,
  11.1168484396999,
  15.1378607061571,
  15.434417533188464,
  40.7726903615389,
  11.46923288239749,
  41.89678193905154,
  122.78089239035424,
  19.92304984647383,
  15.66076906866832,
  27.883494730103067,
  102.75934387501144,
  36.012735343213365,
  22.691350507312702,
  27.252866080233144,
  15.277988413481557,
  20.665325164727093,
  59.37677793786111,
  11.966486557685885,
  8.122196823756076,
  5.041232477753647,
  10.780539126984296,
  16.895224306132164,
  20.48843869720723,
  7.578437925359917,
  19.297314369770156,
  31.88471947070518,
  35.19837571782373,
  54.87603439547455,
  25.230997811952303,
  71.91030404249504,
  13.597968080863915,
  21.952513914708014,
  19.449292463114418,
  8.341334518937044,
  59.999329074837476,
  11.38497874207974,
  6.961264139888034,
  13.00771208931076,
  7.1132072541006135,
  16.811420378171682,
  77.04137400714121,
  5.533212121326684,
  15.391321861575165,
  6.268804524413115,
  24.596142735302863,
  40.43658418948557,
  55.30297277765685,
  16.760382768928178,
  16.519302034079562,
  12.27867586412826,
  21.291185000968785,
  51.16257862654765,
  14.31632160919657,
  6.557504010903474,
  66.64742071867056,
  46.633192304423915,
  11.743559607178831,
  25.314954189756374,
  31.162060904972744,
  5.211515439271813,
  10.378713614348863,
  98.00176264788041,
  16.907295975214236,
  5.472559157478296,
  14.272425737459907,
  88.81167535663668,
  11.240867233831823,
  42.02708852416934,
  5.585498339324491,
  15.812385753677258,
  40.09483634410534,
  40.45859063594583,
  15.229290694036886,
  30.59559926462969,
  25.98196220755276,
  38.58029503520994,
  20.379248927128046,
  19.986745785314824,
  16.305593850579523,
  11.164446095152494,
  36.172450538249336,
  67.67295146877589,
  116.41649689985026,
  18.741962558021726,
  19.083042167773066,
  10.4630148627016,
  67.72030521534391,
  59.82808442854977,
  11.172606097271364,
  70.31555036349765,
  8.552947939283497,
  19.20767016588948,
  30.814969696915252,
  10.92591029708916,
  17.967127742497667,
  58.06875056584753,
  36.29990209987562,
  24.249659806568168,
  14.471739836965908,
  36.67635838868224,
  8.256167712041293,
  21.933269233382674,
  23.72130813587659,
  66.6259338899255,
  20.635910129127666,
  22.379783319646005,
  12.681214690936166,
  8.529439566140935,
  62.34803907660725,
  6.675531474630341,
  35.633463154170656,
  9.814290307496611,
  95.886496185331,
  15.590193431898967,
  65.90200675322568,
  20.4587117178921,
  25.627230975421504,
  13.064130728992723,
  43.20720915324519,
  31.278666802864652,
  45.33777963589865,
  46.146959600607595,
  18.956521288079514,
  7.156548773757226,
  8.571271260729043,
  46.92798890070423,
  5.682290111987992,
  18.40784042348468,
  12.131004937511317,
  6.373642381055459,
  21.931660196806877,
  20.065618757236813,
  63.510783073412284,
  22.920092270423453,
  112.33175443260717,
  13.27371108081708,
  5.716294736428011,
  12.81975667709802,
  70.82363526510603,
  11.28756711526393,
  23.183950803297233,
  33.732320090633976,
  63.94673550935585,
  14.183702210261345,
  5.2953254522313955,
  60.39336210447182,
  34.452494746811595,
  11.396605876286944,
  11.651889344045538,
  22.701832172064627,
  14.452035460802747,
  17.66267803526821,
  6.505809198579629,
  16.48959972276387,
  66.65388061844983,
  14.955776192994849,
  33.92185951219833,
  17.86996744696716,
  24.474558195116117,
  59.26694594642727,
  45.56915405824121,
  8.595740928366824,
  19.037648175269272,
  25.35777756700905,
  45.40384899156447,
  11.693675053578785,
  35.209035827186035,
  23.639042189965767,
  15.505203043197842,
  68.66356314370131,
  18.19565063851373,
  16.82124653958557,
  23.362626243755667,
  14.87575813333602,
  72.94167640088946,
  22.06683624877233,
  13.006993305602885,
  29.41842898033386,
  23.471531735608902,
  75.45367467594446,
  83.88596928035406,
  20.166283026243185,
  6.05005583064838,
  9.140306808718668,
  25.884140793399066,
  9.258400598827512,
  6.993925299398916,
  18.610203750433183,
  59.113875314393475,
  28.1253274360485,
  11.548869799180716,
  67.55106370241315,
  55.9881592055243,
  34.817506184456235,
  6.97554827512385,
  26.026504858549526,
  10.662979206593354,
  35.85399359675375,
  21.505246454104626,
  26.99636480989424,
  5.567337491295738,
  18.18712707321843,
  44.05514089594186,
  13.578008233895376,
  31.509836926867912,
  9.129941566361667,
  5.889746723795036,
  16.642400753730172,
  15.18501496634751,
  27.0943252053947,
  19.076704642784655,
  6.691163229226904,
  57.66642499346492,
  16.615214856301986,
  5.323664037217394,
  12.870639758178005,
  13.47855306367875,
  84.86008029315596,
  15.024038021536226,
  21.08832857150172,
  20.502913106869094,
  33.75538405899368,
  6.733575600312667,
  15.809364813691337,
  7.201286681865446,
  6.0816451145323285,
  79.26252707616494,
  13.726517039074029,
  47.781471066476385,
  22.24930516414614,
  60.7518547946961,
  31.221493428716915,
  24.966690284075703,
  14.347249701326785,
  21.599194252546575,
  38.3588353782252,
  15.564421283043233,
  13.597774954942903,
  39.99628468780775,
  22.88995648003838,
  40.48202396444261,
  35.67166497409103,
  134.39385721480718,
  5.300263663531363,
  85.30866242084073,
  69.29067386173162,
  41.64134077449923,
  43.43278645042357,
  59.87782170498564,
  55.78212196213634,
  6.899441414885073,
  14.38002948266162,
  38.737229854579226,
  19.287070746067666,
  19.365617318460707,
  7.858825854702455,
  6.779570110032547,
  8.623146387587262,
  9.429939411889002,
  5.3678137539094415,
  16.7341984319211,
  22.73427983825978,
  11.966264598200269,
  6.834493365791316,
  5.379952966447907,
  39.778012427051905,
  26.206673266390045,
  12.036368094820142,
  37.99691192933244,
  49.06396024947829,
  146.17828993538492,
  55.25260988315125,
  37.1440217974349,
  15.787236319903396,
  16.647127300143982,
  6.12673853908567,
  5.659348981629775,
  22.30105131438771,
  16.729821211280854,
  29.051625317992823,
  10.811494287935094,
  25.861532066629692,
  34.97851508179406,
  67.02990443086271,
  24.414695995324138,
  16.379813581841642,
  23.11032814379475,
  55.247479014106645,
  61.38022337818151,
  24.324547752369202,
  22.241070093804975,
  7.314593122274525,
  80.20173111609554,
  77.0908019147634,
  102.9154672109412,
  22.957495580830614,
  33.765342634374036,
  46.24134287797399,
  69.1674379915813,
  6.375937003992957,
  8.713720782090602,
  100.05871615653862,
  28.610017911759115,
  79.28598522153132,
  16.865336857350403,
  64.55196103296775,
  9.53284924231757,
  19.443259794927265,
  5.033671761238463,
  10.470431196475305,
  25.097019873028934,
  35.04574341438455,
  27.402774751803868,
  12.694683810796962,
  8.340559177073482,
  13.640831375927062,
  54.81265703640567,
  49.06493623588928,
  5.62065897658044,
  8.794972639428673,
  87.87589858414329,
  30.56781556105499,
  47.260721922823905,
  78.42667838875936,
  48.34745105241726,
  22.619691068853466,
  68.67742711842382,
  6.318123902707775,
  5.5245489791167115,
  5.7359263785661305,
  52.12040622413002,
  11.50279764375112,
  57.866386987602226,
  42.608079697869286,
  85.535365369179,
  10.253058998940517,
  13.472752065084217,
  42.71999881154682,
  24.339973270093914,
  6.913804574393271,
  23.523027679021425,
  79.90497877345776,
  6.785036547183423,
  30.94330828229974,
  56.959207546751344,
  104.11825331105041,
  20.765798431667406,
  6.6513033600003375,
  47.76379734947511,
  37.89704300608452,
  8.56893262957653,
  7.885937528121809,
  26.27196888419597,
  26.26207190317276,
  51.051474998370665,
  14.108970344869345,
  6.773956792795996,
  33.69043216476397,
  92.38500111882406,
  40.292837394265725,
  35.630836819228385,
  74.14495319564793,
  23.83770323817886,
  8.822136888904623,
  13.207788281550927,
  101.69144236874027,
  50.83514608014728,
  11.361083673417486,
  105.45912645222798,
  17.43275333135287,
  10.799588486890835,
  26.664916885708884,
  30.04018034377117,
  11.969509083779881,
  26.214602777262158,
  47.738302108751135,
  26.654833699605224,
  53.26300719531444,
  36.57322617394971,
  11.079982486593522,
  19.746098796914005,
  13.11481075609913,
  6.431811899239111,
  47.500103614519844,
  28.433849245706867,
  5.328613322419715,
  11.888137753680944,
  16.759864257868692,
  9.6620400973615,
  57.540275366133756,
  18.81247992716082,
  35.09621044822148,
  7.865318277315491,
  28.762790559811386,
  33.41290685122884,
  5.432056116735115,
  10.007605498857345,
  6.830449368245077,
  37.41825647814323,
  15.793492646300152,
  7.9181812397333475,
  60.219999080497615,
  51.71645802030075,
  12.279336418927917,
  5.339399799810742,
  11.160836850946943,
  13.716416215455652,
  51.61698277085947,
  9.043224983718092,
  88.66804603884162,
  45.37083132518059,
  21.57801995895946,
  77.2429299500754,
  62.76121402777798,
  7.603611602313181,
  35.33944535508393,
  30.609880227178472,
  12.177713182778449,
  18.675053405733717,
  7.654187800790934,
  54.36381559419603,
  27.994434279752134,
  12.665632173913874,
  45.660785740555816,
  18.92654886503231,
  32.17850497185195,
  18.233160117141868,
  40.79318550543543,
  54.72864586884015,
  32.364631241112306,
  12.836469344741218,
  18.740895444816292,
  5.629464643133254,
  12.523195639411641,
  17.26606467230877,
  22.11263852921472,
  6.407171551870401,
  49.31818741691774,
  7.333911131171782,
  81.60236692441008,
  9.725901949970657,
  11.261091657772775,
  11.530471496658208,
  6.702385151807667,
  14.80563544414959,
  29.65868461440448,
  7.804736296517852,
  11.551961936322066,
  99.99260820493305,
  24.514330279134825,
  17.07862801704698,
  42.52883438094541,
  5.978605888754163,
  75.74311853146145,
  97.90859183660349,
  56.948107276915316,
  67.48149726836735,
  33.42602915415409,
  33.50280221714644,
  21.360916216218968,
  41.48654089270738,
  5.896252938555228,
  81.94242824717769,
  77.20691742822233,
  10.152184116461445,
  13.017087975525696,
  23.577173353366046,
  109.27527899189722,
  114.18203835613429,
  50.30175731612823,
  5.146128300237676,
  12.865746802158746,
  33.40661388393106,
  69.20976579662678,
  53.32857712598188,
  81.99181316987061,
  91.74674516924516,
  141.83846404564252,
  15.66922277029354,
  41.07080556221892,
  69.69847771686334])

Energy spectrum

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

([2785207.8125,
  3082431.25,
  3120943.75,
  3157382.8125,
  3202934.375,
  3272376.5625,
  3316748.4375,
  3388651.5625,
  3402350.0,
  3402814.0625,
  3462937.5,
  3538078.125,
  3541453.125,
  3546592.1875,
  3978398.4375,
  5193150.0,
  5193853.125,
  5193853.125,
  5288400.0,
  5290445.3125,
  5290582.8125,
  5291470.3125,
  5291560.9375,
  5518459.375,
  5521100.0,
  5522295.3125,
  5527396.875,
  5549023.4375,
  5589784.375,
  5674087.5,
  5734862.5,
  5735451.5625,
  5736109.375,
  5766712.5,
  5768162.5,
  5770014.0625,
  5770028.125,
  5771626.5625,
  5773115.625,
  5774379.6875,
  5775762.5,
  5777510.9375,
  5777618.75,
  5778564.0625,
  5781085.9375,
  5811298.4375,
  5815667.1875,
  5820326.5625,
  5832665.625,
  5836910.9375,
  5847814.0625,
  5853029.6875,
  5855212.5,
  5857195.3125,
  5861312.5,
  5896285.9375,
  5926950.0,
  5943704.6875,
  5944151.5625,
  5947159.375,
  5947201.5625,
  5948203.125,
  5951895.3125,
  5953640.625,
  5954620.3125,
  5957484.375,
  5984882.8125,
  6002948.4375,
  6007176.5625,
  6008689.0625,
  6009878.125,
  6013659.375,
  6016290.625,
  6026440.625,
  6052370.3125,
  6052632.8125,
  6054620.3125,
  6059753.125,
  6065460.9375,
  6067029.6875,
  6068957.8125,
  6082809.375,
  6082829.6875,
  6082945.3125,
  6088596.875,
  6107103.125,
  6109231.25,
  6109357.8125,
  6109453.125,
  6109592.1875,
  6109603.125,
  6109667.1875,
  6110323.4375,
  6112467.1875,
  6113878.125,
  6114006.25,
  6114971.875,
  6115895.3125,
  6120559.375,
  6121596.875,
  6122051.5625,
  6137234.375,
  6157212.5,
  6170460.9375,
  6184042.1875,
  6203723.4375,
  6226789.0625,
  6226837.5,
  6244317.1875,
  6245564.0625,
  6253967.1875,
  6254454.6875,
  6254785.9375,
  6256253.125,
  6256690.625,
  6257285.9375,
  6257748.4375,
  6258062.5,
  6258906.25,
  6259164.0625,
  6259389.0625,
  6259390.625,
  6259434.375,
  6259456.25,
  6259920.3125,
  6259973.4375,
  6261121.875,
  6261453.125,
  6262770.3125,
  6265626.5625,
  6266168.75,
  6266823.4375,
  6267245.3125,
  6268345.3125,
  6268496.875,
  6270157.8125,
  6270442.1875,
  6271295.3125,
  6272292.1875,
  6272906.25,
  6273165.625,
  6276104.6875,
  6276667.1875,
  6278050.0,
  6278053.125,
  6280064.0625,
  6280557.8125,
  6286060.9375,
  6287243.75,
  6287435.9375,
  6289389.0625,
  6290634.375,
  6296525.0,
  6310990.625,
  6321596.875,
  6334381.25,
  6338228.125,
  6341671.875,
  6342565.625,
  6352170.3125,
  6369259.375,
  6411226.5625,
  6413437.5,
  6413898.4375,
  6416093.75,
  6423415.625,
  6423535.9375,
  6424820.3125,
  6429246.875,
  6436790.625,
  6454834.375,
  6457175.0,
  6462757.8125,
  6462843.75,
  6466093.75,
  6467762.5,
  6475650.0,
  6519840.625,
  6553476.5625,
  6561409.375,
  6562390.625,
  6581778.125,
  6638431.25,
  6757337.5,
  7017410.9375,
  7256976.5625,
  8369315.625,
  8591932.8125,
  9111762.5,
  9183620.3125,
  9436548.4375,
  9449868.75,
  9467976.5625,
  9517912.5,
  9518856.25,
  9525717.1875,
  9561678.125,
  9580257.8125,
  9594206.25,
  9594498.4375,
  9600270.3125,
  9627248.4375,
  9629625.0,
  9629626.5625,
  9630092.1875,
  9638162.5,
  9783900.0,
  9784475.0,
  9787135.9375,
  9787645.3125,
  9819854.6875,
  9823010.9375,
  9825565.625,
  9825842.1875,
  9827115.625,
  9827421.875,
  9827426.5625,
  9836901.5625,
  9836909.375,
  9838504.6875,
  9839989.0625,
  9842937.5,
  9843367.1875,
  9844545.3125,
  9845146.875,
  9845664.0625,
  9845945.3125,
  9866448.4375,
  9877801.5625,
  9879606.25,
  9879992.1875,
  9880634.375,
  9880693.75,
  9888973.4375,
  9904965.625,
  9966787.5,
  10057748.4375,
  10064606.25,
  10065767.1875,
  10066778.125,
  10089332.8125,
  10090520.3125,
  10097212.5,
  10110789.0625,
  10203009.375,
  10223107.8125,
  10237484.375,
  10253071.875,
  10261001.5625,
  10545165.625,
  10644046.875,
  10653407.8125,
  10852626.5625,
  10896526.5625,
  10910004.6875,
  10910665.625,
  10913295.3125,
  10916296.875,
  10916959.375,
  10921948.4375,
  10934098.4375,
  10934223.4375,
  10937068.75,
  10938762.5,
  10961723.4375,
  11055554.6875,
  11061625.0,
  11061629.6875,
  11065234.375,
  11065939.0625,
  11068700.0,
  11068701.5625,
  11071746.875,
  11072135.9375,
  11073032.8125,
  11073306.25,
  11073357.8125,
  11073520.3125,
  11075271.875,
  11075721.875,
  11076321.875,
  11076446.875,
  11076973.4375,
  11077978.125,
  11078221.875,
  11078734.375,
  11079070.3125,
  11079393.75,
  11079415.625,
  11079607.8125,
  11080146.875,
  11080193.75,
  11080298.4375,
  11080310.9375,
  11080645.3125,
  11080929.6875,
  11080971.875,
  11081698.4375,
  11081718.75,
  11082050.0,
  11082281.25,
  11082307.8125,
  11082323.4375,
  11082353.125,
  11082643.75,
  11082701.5625,
  11082745.3125,
  11082853.125,
  11083339.0625,
  11083515.625,
  11083525.0,
  11083562.5,
  11083687.5,
  11083854.6875,
  11083906.25,
  11084318.75,
  11084345.3125,
  11084723.4375,
  11085229.6875,
  11085325.0,
  11085395.3125,
  11085893.75,
  11086029.6875,
  11086185.9375,
  11086209.375,
  11086237.5,
  11087143.75,
  11087195.3125,
  11087578.125,
  11087856.25,
  11087867.1875,
  11088007.8125,
  11088218.75,
  11088292.1875,
  11088304.6875,
  11088304.6875,
  11088334.375,
  11088343.75,
  11088398.4375,
  11088425.0,
  11088487.5,
  11088784.375,
  11088882.8125,
  11089237.5,
  11089368.75,
  11089403.125,
  11089643.75,
  11089903.125,
  11090220.3125,
  11090796.875,
  11090929.6875,
  11090968.75,
  11091234.375,
  11091460.9375,
  11091546.875,
  11091678.125,
  11092253.125,
  11092732.8125,
  11092785.9375,
  11093701.5625,
  11094317.1875,
  11094365.625,
  11094393.75,
  11095151.5625,
  11095529.6875,
  11097075.0,
  11097107.8125,
  11097926.5625,
  11098078.125,
  11099054.6875,
  11099410.9375,
  11099456.25,
  11099535.9375,
  11099950.0,
  11100032.8125,
  11100079.6875,
  11100256.25,
  11100276.5625,
  11101029.6875,
  11101259.375,
  11101271.875,
  11101281.25,
  11101484.375,
  11101603.125,
  11101653.125,
  11101803.125,
  11102021.875,
  11102029.6875,
  11102059.375,
  11102203.125,
  11102985.9375,
  11103226.5625,
  11103518.75,
  11103614.0625,
  11103675.0,
  11103842.1875,
  11104142.1875,
  11104189.0625,
  11104301.5625,
  11104323.4375,
  11104648.4375,
  11104776.5625,
  11104896.875,
  11105165.625,
  11105376.5625,
  11105423.4375,
  11105507.8125,
  11105529.6875,
  11105667.1875,
  11105717.1875,
  11105789.0625,
  11105839.0625,
  11105862.5,
  11105876.5625,
  11106178.125,
  11106242.1875,
  11106457.8125,
  11106540.625,
  11106589.0625,
  11106614.0625,
  11106667.1875,
  11106676.5625,
  11106762.5,
  11107221.875,
  11107334.375,
  11107454.6875,
  11107496.875,
  11107514.0625,
  11107515.625,
  11107596.875,
  11107668.75,
  11107745.3125,
  11107826.5625,
  11107850.0,
  11108053.125,
  11108179.6875,
  11108443.75,
  11108637.5,
  11108671.875,
  11108693.75,
  11108790.625,
  11108793.75,
  11108915.625,
  11108946.875,
  11108957.8125,
  11108987.5,
  11109128.125,
  11109173.4375,
  11109184.375,
  11109207.8125,
  11109217.1875,
  11109218.75,
  11109468.75,
  11109476.5625,
  11109526.5625,
  11109528.125,
  11109532.8125,
  11109546.875,
  11109589.0625,
  11109648.4375,
  11109682.8125,
  11109762.5,
  11109764.0625,
  11109815.625,
  11109823.4375,
  11109900.0,
  11109923.4375,
  11109978.125,
  11110109.375,
  11110126.5625,
  11110157.8125,
  11110182.8125,
  11110335.9375,
  11110478.125,
  11110481.25,
  11110584.375,
  11110587.5,
  11110662.5,
  11110773.4375,
  11110832.8125,
  11110976.5625,
  11110981.25,
  11111006.25,
  11111042.1875,
  11111175.0,
  11111257.8125,
  11111454.6875,
  11111690.625,
  11111814.0625,
  11111854.6875,
  11111932.8125,
  11112251.5625,
  11112401.5625,
  11112442.1875,
  11112443.75,
  11112501.5625,
  11112504.6875,
  11112654.6875,
  11112684.375,
  11112864.0625,
  11112985.9375,
  11112992.1875,
  11113348.4375,
  11113796.875,
  11114065.625,
  11114189.0625,
  11114834.375,
  11114839.0625,
  11115062.5,
  11115182.8125,
  11115473.4375,
  11115515.625,
  11115557.8125,
  11115837.5,
  11115998.4375,
  11116151.5625,
  11116334.375,
  11116417.1875,
  11116784.375,
  11117264.0625,
  11117395.3125,
  11117406.25,
  11117420.3125,
  11117450.0,
  11118028.125,
  11118385.9375,
  11118392.1875,
  11118506.25,
  11118564.0625,
  11118782.8125,
  11118900.0,
  11118932.8125,
  11119284.375,
  11119553.125,
  11119553.125,
  11119893.75,
  11120350.0,
  11120417.1875,
  11120482.8125,
  11120735.9375,
  11120801.5625,
  11120812.5,
  11120928.125,
  11120996.875,
  11121046.875,
  11121134.375,
  11121464.0625,
  11121756.25,
  11121914.0625,
  11122115.625,
  11122132.8125,
  11122139.0625,
  11122153.125,
  11122168.75,
  11122171.875,
  11122265.625,
  11122351.5625,
  11122634.375,
  11122896.875,
  11123429.6875,
  11123431.25,
  11123457.8125,
  11123589.0625,
  11123618.75,
  11123646.875,
  11123748.4375,
  11123879.6875,
  11123895.3125,
  11124029.6875,
  11124148.4375,
  11124975.0,
  11125065.625,
  11125103.125,
  11125273.4375,
  11125634.375,
  11125673.4375,
  11125709.375,
  11125756.25,
  11125792.1875,
  11126028.125,
  11126160.9375,
  11126193.75,
  11126301.5625,
  11126320.3125,
  11126664.0625,
  11126760.9375,
  11126865.625,
  11126873.4375,
  11127034.375,
  11127121.875,
  11127193.75,
  11127804.6875,
  11128012.5,
  11128984.375,
  11129028.125,
  11129656.25,
  11129875.0,
  11129925.0,
  11129929.6875,
  11129989.0625,
  11130218.75,
  11130428.125,
  11130431.25,
  11130459.375,
  11130676.5625,
  11130693.75,
  11130704.6875,
  11130859.375,
  11130937.5,
  11130945.3125,
  11131132.8125,
  11131207.8125,
  11131667.1875,
  11131887.5,
  11132007.8125,
  11132314.0625,
  11132329.6875,
  11132342.1875,
  11132429.6875,
  11132759.375,
  11132781.25,
  11132870.3125,
  11133593.75,
  11133628.125,
  11133837.5,
  11133950.0,
  11134250.0,
  11134343.75,
  11134365.625,
  11134392.1875,
  11134403.125,
  11134665.625,
  11134707.8125,
  11134935.9375,
  11135507.8125,
  11135540.625,
  11135845.3125,
  11136234.375,
  11136250.0,
  11136298.4375,
  11136492.1875,
  11136496.875,
  11136520.3125,
  11136692.1875,
  11136929.6875,
  11137028.125,
  11137075.0,
  11137279.6875,
  11137348.4375,
  11137368.75,
  11137387.5,
  11137432.8125,
  11137906.25,
  11137917.1875,
  11138284.375,
  11138446.875,
  11138593.75,
  11138715.625,
  11138945.3125,
  11139017.1875,
  11139164.0625,
  11139339.0625,
  11139967.1875,
  11140118.75,
  11140201.5625,
  11140876.5625,
  11140879.6875,
  11140982.8125,
  11141414.0625,
  11141537.5,
  11141678.125,
  11141781.25,
  11141795.3125,
  11141985.9375,
  11142148.4375,
  11142553.125,
  11143376.5625,
  11143431.25,
  11143962.5,
  11144178.125,
  11144229.6875,
  11144285.9375,
  11144418.75,
  11144917.1875,
  11145003.125,
  11145079.6875,
  11145237.5,
  11145259.375,
  11145731.25,
  11146504.6875,
  11147417.1875,
  11147837.5,
  11148417.1875,
  11148946.875,
  11149095.3125,
  11149917.1875,
  11150837.5,
  11152126.5625,
  11152293.75,
  11152412.5,
  11152615.625,
  11152701.5625,
  11153057.8125,
  11153790.625,
  11154282.8125,
  11154315.625,
  11154379.6875,
  11154762.5,
  11154815.625,
  11155117.1875,
  11155557.8125,
  11156307.8125,
  11156856.25,
  11156921.875,
  11157521.875,
  11157682.8125,
  11158201.5625,
  11158256.25,
  11158400.0,
  11158575.0,
  11158640.625,
  11158940.625,
  11159389.0625,
  11160835.9375,
  11162770.3125,
  11162832.8125,
  11162939.0625,
  11163557.8125,
  11163739.0625,
  11163743.75,
  11163989.0625,
  11164571.875,
  11165528.125,
  11165548.4375,
  11166264.0625,
  11167178.125,
  11167421.875,
  11168056.25,
  11168268.75,
  11168671.875,
  11169267.1875,
  11169307.8125,
  11169353.125,
  11169476.5625,
  11169884.375,
  11170520.3125,
  11170590.625,
  11170656.25,
  11172698.4375,
  11173348.4375,
  11174301.5625,
  11176385.9375,
  11176690.625,
  11179804.6875,
  11180020.3125,
  11180125.0,
  11181593.75,
  11190251.5625],
 [5.152539720628002,
  6.854280608202352,
  28.46977324915323,
  8.518592729869086,
  104.12549718709128,
  10.851741588919381,
  5.517981870959336,
  66.07082341469258,
  5.07205040950528,
  70.83343627793475,
  13.477841232271121,
  21.514784216537013,
  7.294620685653026,
  71.59663946091045,
  19.04167810178794,
  101.8615281647461,
  38.31463262727778,
  15.748185106962069,
  22.337371770480527,
  41.86135388760191,
  24.76487392530209,
  42.20275329711578,
  22.83905188675709,
  64.1082339267195,
  26.240131209688553,
  59.809587418109714,
  44.944165232861174,
  39.24467288238996,
  90.35129801094376,
  30.49577284348799,
  120.38109198228508,
  11.191893569629812,
  60.31832820016083,
  28.40224654850701,
  11.554734225267199,
  40.18066455224005,
  36.20438031671707,
  28.351095269272477,
  57.1538165009531,
  51.573978550881286,
  87.88166151492781,
  36.69209596588307,
  13.675773006736271,
  28.839588101807237,
  55.04759646551981,
  5.687725581577804,
  20.92132427970974,
  61.755145281881425,
  45.54590494190182,
  66.63260202559736,
  18.41951993227862,
  9.361202130713966,
  41.85733994915828,
  40.41770608476546,
  77.83344741274658,
  33.50197557675041,
  15.768944741324155,
  67.71139077171931,
  84.67955027596174,
  6.011356114618717,
  6.125334205525346,
  92.5483191641962,
  135.5897747452631,
  5.960620084817627,
  40.39701737275518,
  10.540728353701452,
  39.45172103609671,
  76.39916273657171,
  63.80271822176243,
  5.790528280843724,
  7.200659289731451,
  11.509572524544506,
  80.01810437947763,
  5.171117644277947,
  96.27683852930934,
  186.1417397460132,
  8.458468574236372,
  13.888748650162718,
  28.979197039357043,
  14.203136917788044,
  19.147186778991607,
  20.71329242183629,
  25.29131804774221,
  5.413523852409813,
  130.16779268193548,
  34.235723812894264,
  7.014783924191486,
  67.48605592880169,
  171.36108906068372,
  34.293081690683756,
  25.59189536940345,
  5.875564676840401,
  130.5400459512419,
  180.8139704964108,
  24.40088221057529,
  30.664071638417916,
  49.29916302104552,
  19.211923297463933,
  38.35488307977782,
  7.194868237328745,
  23.600278627113184,
  10.923813758005368,
  47.607297203554104,
  13.970180127193217,
  5.32309791551664,
  109.8117977837779,
  18.285130042481654,
  5.542866385732408,
  43.80287967358354,
  14.510289209141238,
  22.583907452686088,
  24.847664972428685,
  33.82653267657091,
  21.16403345986571,
  171.12548785117116,
  116.9795915399358,
  66.17091062104805,
  9.443723749831062,
  12.087585797711258,
  31.973466546309027,
  16.435148161313535,
  41.212228278081646,
  15.900425747520845,
  20.733557048683927,
  25.197709162988858,
  65.82600010526188,
  20.129005151250173,
  11.749626610644865,
  139.47939737716968,
  68.48859484153373,
  24.977419719971078,
  14.152626115950188,
  14.780748339516347,
  25.40293149882287,
  76.51991348766333,
  17.54669692670169,
  19.305919145468287,
  60.40740473952349,
  7.2109320187492,
  39.20246516494882,
  48.269934970371956,
  70.74524123177379,
  56.14594605387755,
  11.086447050420286,
  45.2014525582852,
  25.559896113786696,
  25.11342373677851,
  17.169847986045625,
  79.73208922562054,
  52.590421189295434,
  36.095362761543846,
  5.346624621020493,
  42.44837923419093,
  87.28161169490454,
  25.881576481000494,
  8.137587074235956,
  53.190537051731575,
  17.887233167459463,
  31.59946636889701,
  24.074466485893478,
  19.496495567196902,
  26.346880458735214,
  15.315357157168105,
  49.99310863578063,
  80.9253123742679,
  9.473031339365475,
  22.486791936312983,
  11.662643518001484,
  71.85230768789084,
  16.735752443184825,
  103.20065404595042,
  7.586015517187785,
  20.600063522329673,
  5.960674491227817,
  85.55176226011383,
  75.01433825157935,
  20.08003280395637,
  12.648045863603514,
  9.51613277849609,
  7.93866252653696,
  86.82205887681428,
  72.23923501461027,
  10.401945538605132,
  7.95739135679953,
  16.391244436923426,
  40.535221212822655,
  34.272901388006744,
  54.51319317104874,
  10.935181603713145,
  41.68904159387064,
  36.64664766442374,
  60.270084587438525,
  51.54528200840906,
  41.04865140503001,
  18.521115835466233,
  77.97960074148779,
  42.40383516276643,
  8.043868596005852,
  18.45246269732988,
  49.958918820606584,
  56.28467471045701,
  81.4668970979767,
  25.588125516294543,
  88.46723420294437,
  49.20938784171525,
  29.492253688002982,
  26.113114805826186,
  14.513556031106724,
  46.422392689496654,
  10.1746537804279,
  11.483746364289411,
  7.318606183786745,
  6.346254242713159,
  22.096643556250786,
  23.224134371512825,
  44.04413327415761,
  18.522602832307346,
  29.05964972547983,
  22.509880825783704,
  7.700183749296881,
  65.02246401930786,
  68.74819988327873,
  43.97472934340662,
  108.84548564690715,
  7.320219495069112,
  6.2761340964351415,
  13.336095984452928,
  55.226124766388004,
  45.175026213035906,
  114.73517871425935,
  15.995647993758702,
  6.930998450410038,
  5.031119491253183,
  7.396871343343425,
  6.767683151708588,
  155.31080594303745,
  27.818023408680585,
  21.12701594558999,
  29.66940350703561,
  5.137249344001191,
  50.33604985175072,
  71.5731356473262,
  18.301471469860463,
  20.948050266697503,
  100.21797413158507,
  31.743090192893078,
  5.67603872794253,
  54.764653430482234,
  70.81253963895017,
  28.084685354651697,
  86.83339496681926,
  20.704348569831264,
  106.73138931339427,
  8.623631646227464,
  57.041228718383074,
  11.20262373289785,
  43.64437065664903,
  82.90455393049896,
  25.242442820201994,
  16.884654122020006,
  84.29590427404828,
  15.943988677524754,
  123.91754786059644,
  20.105722669737105,
  23.761109782440226,
  5.153674506095576,
  34.220916806604734,
  12.41614136873189,
  5.473770256582496,
  10.890228747204613,
  28.15591359082218,
  8.062699944680636,
  20.74396740876828,
  128.9700307914263,
  58.59504639904654,
  48.804394379641764,
  21.95333066002647,
  11.1168484396999,
  15.1378607061571,
  15.434417533188464,
  40.7726903615389,
  11.46923288239749,
  41.89678193905154,
  122.78089239035424,
  19.92304984647383,
  15.66076906866832,
  27.883494730103067,
  102.75934387501144,
  36.012735343213365,
  22.691350507312702,
  27.252866080233144,
  15.277988413481557,
  20.665325164727093,
  59.37677793786111,
  11.966486557685885,
  8.122196823756076,
  5.041232477753647,
  10.780539126984296,
  16.895224306132164,
  20.48843869720723,
  7.578437925359917,
  19.297314369770156,
  31.88471947070518,
  35.19837571782373,
  54.87603439547455,
  25.230997811952303,
  71.91030404249504,
  13.597968080863915,
  21.952513914708014,
  19.449292463114418,
  8.341334518937044,
  59.999329074837476,
  11.38497874207974,
  6.961264139888034,
  13.00771208931076,
  7.1132072541006135,
  16.811420378171682,
  77.04137400714121,
  5.533212121326684,
  15.391321861575165,
  6.268804524413115,
  24.596142735302863,
  40.43658418948557,
  55.30297277765685,
  16.760382768928178,
  16.519302034079562,
  12.27867586412826,
  21.291185000968785,
  51.16257862654765,
  14.31632160919657,
  6.557504010903474,
  66.64742071867056,
  46.633192304423915,
  11.743559607178831,
  25.314954189756374,
  31.162060904972744,
  5.211515439271813,
  10.378713614348863,
  98.00176264788041,
  16.907295975214236,
  5.472559157478296,
  14.272425737459907,
  88.81167535663668,
  11.240867233831823,
  42.02708852416934,
  5.585498339324491,
  15.812385753677258,
  40.09483634410534,
  40.45859063594583,
  15.229290694036886,
  30.59559926462969,
  25.98196220755276,
  38.58029503520994,
  20.379248927128046,
  19.986745785314824,
  16.305593850579523,
  11.164446095152494,
  36.172450538249336,
  67.67295146877589,
  116.41649689985026,
  18.741962558021726,
  19.083042167773066,
  10.4630148627016,
  67.72030521534391,
  59.82808442854977,
  11.172606097271364,
  70.31555036349765,
  8.552947939283497,
  19.20767016588948,
  30.814969696915252,
  10.92591029708916,
  17.967127742497667,
  58.06875056584753,
  36.29990209987562,
  24.249659806568168,
  14.471739836965908,
  36.67635838868224,
  8.256167712041293,
  21.933269233382674,
  23.72130813587659,
  66.6259338899255,
  20.635910129127666,
  22.379783319646005,
  12.681214690936166,
  8.529439566140935,
  62.34803907660725,
  6.675531474630341,
  35.633463154170656,
  9.814290307496611,
  95.886496185331,
  15.590193431898967,
  65.90200675322568,
  20.4587117178921,
  25.627230975421504,
  13.064130728992723,
  43.20720915324519,
  31.278666802864652,
  45.33777963589865,
  46.146959600607595,
  18.956521288079514,
  7.156548773757226,
  8.571271260729043,
  46.92798890070423,
  5.682290111987992,
  18.40784042348468,
  12.131004937511317,
  6.373642381055459,
  21.931660196806877,
  20.065618757236813,
  63.510783073412284,
  22.920092270423453,
  112.33175443260717,
  13.27371108081708,
  5.716294736428011,
  12.81975667709802,
  70.82363526510603,
  11.28756711526393,
  23.183950803297233,
  33.732320090633976,
  63.94673550935585,
  14.183702210261345,
  5.2953254522313955,
  60.39336210447182,
  34.452494746811595,
  11.396605876286944,
  11.651889344045538,
  22.701832172064627,
  14.452035460802747,
  17.66267803526821,
  6.505809198579629,
  16.48959972276387,
  66.65388061844983,
  14.955776192994849,
  33.92185951219833,
  17.86996744696716,
  24.474558195116117,
  59.26694594642727,
  45.56915405824121,
  8.595740928366824,
  19.037648175269272,
  25.35777756700905,
  45.40384899156447,
  11.693675053578785,
  35.209035827186035,
  23.639042189965767,
  15.505203043197842,
  68.66356314370131,
  18.19565063851373,
  16.82124653958557,
  23.362626243755667,
  14.87575813333602,
  72.94167640088946,
  22.06683624877233,
  13.006993305602885,
  29.41842898033386,
  23.471531735608902,
  75.45367467594446,
  83.88596928035406,
  20.166283026243185,
  6.05005583064838,
  9.140306808718668,
  25.884140793399066,
  9.258400598827512,
  6.993925299398916,
  18.610203750433183,
  59.113875314393475,
  28.1253274360485,
  11.548869799180716,
  67.55106370241315,
  55.9881592055243,
  34.817506184456235,
  6.97554827512385,
  26.026504858549526,
  10.662979206593354,
  35.85399359675375,
  21.505246454104626,
  26.99636480989424,
  5.567337491295738,
  18.18712707321843,
  44.05514089594186,
  13.578008233895376,
  31.509836926867912,
  9.129941566361667,
  5.889746723795036,
  16.642400753730172,
  15.18501496634751,
  27.0943252053947,
  19.076704642784655,
  6.691163229226904,
  57.66642499346492,
  16.615214856301986,
  5.323664037217394,
  12.870639758178005,
  13.47855306367875,
  84.86008029315596,
  15.024038021536226,
  21.08832857150172,
  20.502913106869094,
  33.75538405899368,
  6.733575600312667,
  15.809364813691337,
  7.201286681865446,
  6.0816451145323285,
  79.26252707616494,
  13.726517039074029,
  47.781471066476385,
  22.24930516414614,
  60.7518547946961,
  31.221493428716915,
  24.966690284075703,
  14.347249701326785,
  21.599194252546575,
  38.3588353782252,
  15.564421283043233,
  13.597774954942903,
  39.99628468780775,
  22.88995648003838,
  40.48202396444261,
  35.67166497409103,
  134.39385721480718,
  5.300263663531363,
  85.30866242084073,
  69.29067386173162,
  41.64134077449923,
  43.43278645042357,
  59.87782170498564,
  55.78212196213634,
  6.899441414885073,
  14.38002948266162,
  38.737229854579226,
  19.287070746067666,
  19.365617318460707,
  7.858825854702455,
  6.779570110032547,
  8.623146387587262,
  9.429939411889002,
  5.3678137539094415,
  16.7341984319211,
  22.73427983825978,
  11.966264598200269,
  6.834493365791316,
  5.379952966447907,
  39.778012427051905,
  26.206673266390045,
  12.036368094820142,
  37.99691192933244,
  49.06396024947829,
  146.17828993538492,
  55.25260988315125,
  37.1440217974349,
  15.787236319903396,
  16.647127300143982,
  6.12673853908567,
  5.659348981629775,
  22.30105131438771,
  16.729821211280854,
  29.051625317992823,
  10.811494287935094,
  25.861532066629692,
  34.97851508179406,
  67.02990443086271,
  24.414695995324138,
  16.379813581841642,
  23.11032814379475,
  55.247479014106645,
  61.38022337818151,
  24.324547752369202,
  22.241070093804975,
  7.314593122274525,
  80.20173111609554,
  77.0908019147634,
  102.9154672109412,
  22.957495580830614,
  33.765342634374036,
  46.24134287797399,
  69.1674379915813,
  6.375937003992957,
  8.713720782090602,
  100.05871615653862,
  28.610017911759115,
  79.28598522153132,
  16.865336857350403,
  64.55196103296775,
  9.53284924231757,
  19.443259794927265,
  5.033671761238463,
  10.470431196475305,
  25.097019873028934,
  35.04574341438455,
  27.402774751803868,
  12.694683810796962,
  8.340559177073482,
  13.640831375927062,
  54.81265703640567,
  49.06493623588928,
  5.62065897658044,
  8.794972639428673,
  87.87589858414329,
  30.56781556105499,
  47.260721922823905,
  78.42667838875936,
  48.34745105241726,
  22.619691068853466,
  68.67742711842382,
  6.318123902707775,
  5.5245489791167115,
  5.7359263785661305,
  52.12040622413002,
  11.50279764375112,
  57.866386987602226,
  42.608079697869286,
  85.535365369179,
  10.253058998940517,
  13.472752065084217,
  42.71999881154682,
  24.339973270093914,
  6.913804574393271,
  23.523027679021425,
  79.90497877345776,
  6.785036547183423,
  30.94330828229974,
  56.959207546751344,
  104.11825331105041,
  20.765798431667406,
  6.6513033600003375,
  47.76379734947511,
  37.89704300608452,
  8.56893262957653,
  7.885937528121809,
  26.27196888419597,
  26.26207190317276,
  51.051474998370665,
  14.108970344869345,
  6.773956792795996,
  33.69043216476397,
  92.38500111882406,
  40.292837394265725,
  35.630836819228385,
  74.14495319564793,
  23.83770323817886,
  8.822136888904623,
  13.207788281550927,
  101.69144236874027,
  50.83514608014728,
  11.361083673417486,
  105.45912645222798,
  17.43275333135287,
  10.799588486890835,
  26.664916885708884,
  30.04018034377117,
  11.969509083779881,
  26.214602777262158,
  47.738302108751135,
  26.654833699605224,
  53.26300719531444,
  36.57322617394971,
  11.079982486593522,
  19.746098796914005,
  13.11481075609913,
  6.431811899239111,
  47.500103614519844,
  28.433849245706867,
  5.328613322419715,
  11.888137753680944,
  16.759864257868692,
  9.6620400973615,
  57.540275366133756,
  18.81247992716082,
  35.09621044822148,
  7.865318277315491,
  28.762790559811386,
  33.41290685122884,
  5.432056116735115,
  10.007605498857345,
  6.830449368245077,
  37.41825647814323,
  15.793492646300152,
  7.9181812397333475,
  60.219999080497615,
  51.71645802030075,
  12.279336418927917,
  5.339399799810742,
  11.160836850946943,
  13.716416215455652,
  51.61698277085947,
  9.043224983718092,
  88.66804603884162,
  45.37083132518059,
  21.57801995895946,
  77.2429299500754,
  62.76121402777798,
  7.603611602313181,
  35.33944535508393,
  30.609880227178472,
  12.177713182778449,
  18.675053405733717,
  7.654187800790934,
  54.36381559419603,
  27.994434279752134,
  12.665632173913874,
  45.660785740555816,
  18.92654886503231,
  32.17850497185195,
  18.233160117141868,
  40.79318550543543,
  54.72864586884015,
  32.364631241112306,
  12.836469344741218,
  18.740895444816292,
  5.629464643133254,
  12.523195639411641,
  17.26606467230877,
  22.11263852921472,
  6.407171551870401,
  49.31818741691774,
  7.333911131171782,
  81.60236692441008,
  9.725901949970657,
  11.261091657772775,
  11.530471496658208,
  6.702385151807667,
  14.80563544414959,
  29.65868461440448,
  7.804736296517852,
  11.551961936322066,
  99.99260820493305,
  24.514330279134825,
  17.07862801704698,
  42.52883438094541,
  5.978605888754163,
  75.74311853146145,
  97.90859183660349,
  56.948107276915316,
  67.48149726836735,
  33.42602915415409,
  33.50280221714644,
  21.360916216218968,
  41.48654089270738,
  5.896252938555228,
  81.94242824717769,
  77.20691742822233,
  10.152184116461445,
  13.017087975525696,
  23.577173353366046,
  109.27527899189722,
  114.18203835613429,
  50.30175731612823,
  5.146128300237676,
  12.865746802158746,
  33.40661388393106,
  69.20976579662678,
  53.32857712598188,
  81.99181316987061,
  91.74674516924516,
  141.83846404564252,
  15.66922277029354,
  41.07080556221892,
  69.69847771686334])

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

([2785207.8125,
  3082431.25,
  3120943.75,
  3157382.8125,
  3202934.375,
  3272376.5625,
  3316748.4375,
  3388651.5625,
  3402350.0,
  3402814.0625,
  3462937.5,
  3538078.125,
  3541453.125,
  3546592.1875,
  3978398.4375,
  5193150.0,
  5193853.125,
  5193853.125,
  5288400.0,
  5290445.3125,
  5290582.8125,
  5291470.3125,
  5291560.9375,
  5518459.375,
  5521100.0,
  5522295.3125,
  5527396.875,
  5549023.4375,
  5589784.375,
  5674087.5,
  5734862.5,
  5735451.5625,
  5736109.375,
  5766712.5,
  5768162.5,
  5770014.0625,
  5770028.125,
  5771626.5625,
  5773115.625,
  5774379.6875,
  5775762.5,
  5777510.9375,
  5777618.75,
  5778564.0625,
  5781085.9375,
  5811298.4375,
  5815667.1875,
  5820326.5625,
  5832665.625,
  5836910.9375,
  5847814.0625,
  5853029.6875,
  5855212.5,
  5857195.3125,
  5861312.5,
  5896285.9375,
  5926950.0,
  5943704.6875,
  5944151.5625,
  5947159.375,
  5947201.5625,
  5948203.125,
  5951895.3125,
  5953640.625,
  5954620.3125,
  5957484.375,
  5984882.8125,
  6002948.4375,
  6007176.5625,
  6008689.0625,
  6009878.125,
  6013659.375,
  6016290.625,
  6026440.625,
  6052370.3125,
  6052632.8125,
  6054620.3125,
  6059753.125,
  6065460.9375,
  6067029.6875,
  6068957.8125,
  6082809.375,
  6082829.6875,
  6082945.3125,
  6088596.875,
  6107103.125,
  6109231.25,
  6109357.8125,
  6109453.125,
  6109592.1875,
  6109603.125,
  6109667.1875,
  6110323.4375,
  6112467.1875,
  6113878.125,
  6114006.25,
  6114971.875,
  6115895.3125,
  6120559.375,
  6121596.875,
  6122051.5625,
  6137234.375,
  6157212.5,
  6170460.9375,
  6184042.1875,
  6203723.4375,
  6226789.0625,
  6226837.5,
  6244317.1875,
  6245564.0625,
  6253967.1875,
  6254454.6875,
  6254785.9375,
  6256253.125,
  6256690.625,
  6257285.9375,
  6257748.4375,
  6258062.5,
  6258906.25,
  6259164.0625,
  6259389.0625,
  6259390.625,
  6259434.375,
  6259456.25,
  6259920.3125,
  6259973.4375,
  6261121.875,
  6261453.125,
  6262770.3125,
  6265626.5625,
  6266168.75,
  6266823.4375,
  6267245.3125,
  6268345.3125,
  6268496.875,
  6270157.8125,
  6270442.1875,
  6271295.3125,
  6272292.1875,
  6272906.25,
  6273165.625,
  6276104.6875,
  6276667.1875,
  6278050.0,
  6278053.125,
  6280064.0625,
  6280557.8125,
  6286060.9375,
  6287243.75,
  6287435.9375,
  6289389.0625,
  6290634.375,
  6296525.0,
  6310990.625,
  6321596.875,
  6334381.25,
  6338228.125,
  6341671.875,
  6342565.625,
  6352170.3125,
  6369259.375,
  6411226.5625,
  6413437.5,
  6413898.4375,
  6416093.75,
  6423415.625,
  6423535.9375,
  6424820.3125,
  6429246.875,
  6436790.625,
  6454834.375,
  6457175.0,
  6462757.8125,
  6462843.75,
  6466093.75,
  6467762.5,
  6475650.0,
  6519840.625,
  6553476.5625,
  6561409.375,
  6562390.625,
  6581778.125,
  6638431.25,
  6757337.5,
  7017410.9375,
  7256976.5625,
  8369315.625,
  8591932.8125,
  9111762.5,
  9183620.3125,
  9436548.4375,
  9449868.75,
  9467976.5625,
  9517912.5,
  9518856.25,
  9525717.1875,
  9561678.125,
  9580257.8125,
  9594206.25,
  9594498.4375,
  9600270.3125,
  9627248.4375,
  9629625.0,
  9629626.5625,
  9630092.1875,
  9638162.5,
  9783900.0,
  9784475.0,
  9787135.9375,
  9787645.3125,
  9819854.6875,
  9823010.9375,
  9825565.625,
  9825842.1875,
  9827115.625,
  9827421.875,
  9827426.5625,
  9836901.5625,
  9836909.375,
  9838504.6875,
  9839989.0625,
  9842937.5,
  9843367.1875,
  9844545.3125,
  9845146.875,
  9845664.0625,
  9845945.3125,
  9866448.4375,
  9877801.5625,
  9879606.25,
  9879992.1875,
  9880634.375,
  9880693.75,
  9888973.4375,
  9904965.625,
  9966787.5,
  10057748.4375,
  10064606.25,
  10065767.1875,
  10066778.125,
  10089332.8125,
  10090520.3125,
  10097212.5,
  10110789.0625,
  10203009.375,
  10223107.8125,
  10237484.375,
  10253071.875,
  10261001.5625,
  10545165.625,
  10644046.875,
  10653407.8125,
  10852626.5625,
  10896526.5625,
  10910004.6875,
  10910665.625,
  10913295.3125,
  10916296.875,
  10916959.375,
  10921948.4375,
  10934098.4375,
  10934223.4375,
  10937068.75,
  10938762.5,
  10961723.4375,
  11055554.6875,
  11061625.0,
  11061629.6875,
  11065234.375,
  11065939.0625,
  11068700.0,
  11068701.5625,
  11071746.875,
  11072135.9375,
  11073032.8125,
  11073306.25,
  11073357.8125,
  11073520.3125,
  11075271.875,
  11075721.875,
  11076321.875,
  11076446.875,
  11076973.4375,
  11077978.125,
  11078221.875,
  11078734.375,
  11079070.3125,
  11079393.75,
  11079415.625,
  11079607.8125,
  11080146.875,
  11080193.75,
  11080298.4375,
  11080310.9375,
  11080645.3125,
  11080929.6875,
  11080971.875,
  11081698.4375,
  11081718.75,
  11082050.0,
  11082281.25,
  11082307.8125,
  11082323.4375,
  11082353.125,
  11082643.75,
  11082701.5625,
  11082745.3125,
  11082853.125,
  11083339.0625,
  11083515.625,
  11083525.0,
  11083562.5,
  11083687.5,
  11083854.6875,
  11083906.25,
  11084318.75,
  11084345.3125,
  11084723.4375,
  11085229.6875,
  11085325.0,
  11085395.3125,
  11085893.75,
  11086029.6875,
  11086185.9375,
  11086209.375,
  11086237.5,
  11087143.75,
  11087195.3125,
  11087578.125,
  11087856.25,
  11087867.1875,
  11088007.8125,
  11088218.75,
  11088292.1875,
  11088304.6875,
  11088304.6875,
  11088334.375,
  11088343.75,
  11088398.4375,
  11088425.0,
  11088487.5,
  11088784.375,
  11088882.8125,
  11089237.5,
  11089368.75,
  11089403.125,
  11089643.75,
  11089903.125,
  11090220.3125,
  11090796.875,
  11090929.6875,
  11090968.75,
  11091234.375,
  11091460.9375,
  11091546.875,
  11091678.125,
  11092253.125,
  11092732.8125,
  11092785.9375,
  11093701.5625,
  11094317.1875,
  11094365.625,
  11094393.75,
  11095151.5625,
  11095529.6875,
  11097075.0,
  11097107.8125,
  11097926.5625,
  11098078.125,
  11099054.6875,
  11099410.9375,
  11099456.25,
  11099535.9375,
  11099950.0,
  11100032.8125,
  11100079.6875,
  11100256.25,
  11100276.5625,
  11101029.6875,
  11101259.375,
  11101271.875,
  11101281.25,
  11101484.375,
  11101603.125,
  11101653.125,
  11101803.125,
  11102021.875,
  11102029.6875,
  11102059.375,
  11102203.125,
  11102985.9375,
  11103226.5625,
  11103518.75,
  11103614.0625,
  11103675.0,
  11103842.1875,
  11104142.1875,
  11104189.0625,
  11104301.5625,
  11104323.4375,
  11104648.4375,
  11104776.5625,
  11104896.875,
  11105165.625,
  11105376.5625,
  11105423.4375,
  11105507.8125,
  11105529.6875,
  11105667.1875,
  11105717.1875,
  11105789.0625,
  11105839.0625,
  11105862.5,
  11105876.5625,
  11106178.125,
  11106242.1875,
  11106457.8125,
  11106540.625,
  11106589.0625,
  11106614.0625,
  11106667.1875,
  11106676.5625,
  11106762.5,
  11107221.875,
  11107334.375,
  11107454.6875,
  11107496.875,
  11107514.0625,
  11107515.625,
  11107596.875,
  11107668.75,
  11107745.3125,
  11107826.5625,
  11107850.0,
  11108053.125,
  11108179.6875,
  11108443.75,
  11108637.5,
  11108671.875,
  11108693.75,
  11108790.625,
  11108793.75,
  11108915.625,
  11108946.875,
  11108957.8125,
  11108987.5,
  11109128.125,
  11109173.4375,
  11109184.375,
  11109207.8125,
  11109217.1875,
  11109218.75,
  11109468.75,
  11109476.5625,
  11109526.5625,
  11109528.125,
  11109532.8125,
  11109546.875,
  11109589.0625,
  11109648.4375,
  11109682.8125,
  11109762.5,
  11109764.0625,
  11109815.625,
  11109823.4375,
  11109900.0,
  11109923.4375,
  11109978.125,
  11110109.375,
  11110126.5625,
  11110157.8125,
  11110182.8125,
  11110335.9375,
  11110478.125,
  11110481.25,
  11110584.375,
  11110587.5,
  11110662.5,
  11110773.4375,
  11110832.8125,
  11110976.5625,
  11110981.25,
  11111006.25,
  11111042.1875,
  11111175.0,
  11111257.8125,
  11111454.6875,
  11111690.625,
  11111814.0625,
  11111854.6875,
  11111932.8125,
  11112251.5625,
  11112401.5625,
  11112442.1875,
  11112443.75,
  11112501.5625,
  11112504.6875,
  11112654.6875,
  11112684.375,
  11112864.0625,
  11112985.9375,
  11112992.1875,
  11113348.4375,
  11113796.875,
  11114065.625,
  11114189.0625,
  11114834.375,
  11114839.0625,
  11115062.5,
  11115182.8125,
  11115473.4375,
  11115515.625,
  11115557.8125,
  11115837.5,
  11115998.4375,
  11116151.5625,
  11116334.375,
  11116417.1875,
  11116784.375,
  11117264.0625,
  11117395.3125,
  11117406.25,
  11117420.3125,
  11117450.0,
  11118028.125,
  11118385.9375,
  11118392.1875,
  11118506.25,
  11118564.0625,
  11118782.8125,
  11118900.0,
  11118932.8125,
  11119284.375,
  11119553.125,
  11119553.125,
  11119893.75,
  11120350.0,
  11120417.1875,
  11120482.8125,
  11120735.9375,
  11120801.5625,
  11120812.5,
  11120928.125,
  11120996.875,
  11121046.875,
  11121134.375,
  11121464.0625,
  11121756.25,
  11121914.0625,
  11122115.625,
  11122132.8125,
  11122139.0625,
  11122153.125,
  11122168.75,
  11122171.875,
  11122265.625,
  11122351.5625,
  11122634.375,
  11122896.875,
  11123429.6875,
  11123431.25,
  11123457.8125,
  11123589.0625,
  11123618.75,
  11123646.875,
  11123748.4375,
  11123879.6875,
  11123895.3125,
  11124029.6875,
  11124148.4375,
  11124975.0,
  11125065.625,
  11125103.125,
  11125273.4375,
  11125634.375,
  11125673.4375,
  11125709.375,
  11125756.25,
  11125792.1875,
  11126028.125,
  11126160.9375,
  11126193.75,
  11126301.5625,
  11126320.3125,
  11126664.0625,
  11126760.9375,
  11126865.625,
  11126873.4375,
  11127034.375,
  11127121.875,
  11127193.75,
  11127804.6875,
  11128012.5,
  11128984.375,
  11129028.125,
  11129656.25,
  11129875.0,
  11129925.0,
  11129929.6875,
  11129989.0625,
  11130218.75,
  11130428.125,
  11130431.25,
  11130459.375,
  11130676.5625,
  11130693.75,
  11130704.6875,
  11130859.375,
  11130937.5,
  11130945.3125,
  11131132.8125,
  11131207.8125,
  11131667.1875,
  11131887.5,
  11132007.8125,
  11132314.0625,
  11132329.6875,
  11132342.1875,
  11132429.6875,
  11132759.375,
  11132781.25,
  11132870.3125,
  11133593.75,
  11133628.125,
  11133837.5,
  11133950.0,
  11134250.0,
  11134343.75,
  11134365.625,
  11134392.1875,
  11134403.125,
  11134665.625,
  11134707.8125,
  11134935.9375,
  11135507.8125,
  11135540.625,
  11135845.3125,
  11136234.375,
  11136250.0,
  11136298.4375,
  11136492.1875,
  11136496.875,
  11136520.3125,
  11136692.1875,
  11136929.6875,
  11137028.125,
  11137075.0,
  11137279.6875,
  11137348.4375,
  11137368.75,
  11137387.5,
  11137432.8125,
  11137906.25,
  11137917.1875,
  11138284.375,
  11138446.875,
  11138593.75,
  11138715.625,
  11138945.3125,
  11139017.1875,
  11139164.0625,
  11139339.0625,
  11139967.1875,
  11140118.75,
  11140201.5625,
  11140876.5625,
  11140879.6875,
  11140982.8125,
  11141414.0625,
  11141537.5,
  11141678.125,
  11141781.25,
  11141795.3125,
  11141985.9375,
  11142148.4375,
  11142553.125,
  11143376.5625,
  11143431.25,
  11143962.5,
  11144178.125,
  11144229.6875,
  11144285.9375,
  11144418.75,
  11144917.1875,
  11145003.125,
  11145079.6875,
  11145237.5,
  11145259.375,
  11145731.25,
  11146504.6875,
  11147417.1875,
  11147837.5,
  11148417.1875,
  11148946.875,
  11149095.3125,
  11149917.1875,
  11150837.5,
  11152126.5625,
  11152293.75,
  11152412.5,
  11152615.625,
  11152701.5625,
  11153057.8125,
  11153790.625,
  11154282.8125,
  11154315.625,
  11154379.6875,
  11154762.5,
  11154815.625,
  11155117.1875,
  11155557.8125,
  11156307.8125,
  11156856.25,
  11156921.875,
  11157521.875,
  11157682.8125,
  11158201.5625,
  11158256.25,
  11158400.0,
  11158575.0,
  11158640.625,
  11158940.625,
  11159389.0625,
  11160835.9375,
  11162770.3125,
  11162832.8125,
  11162939.0625,
  11163557.8125,
  11163739.0625,
  11163743.75,
  11163989.0625,
  11164571.875,
  11165528.125,
  11165548.4375,
  11166264.0625,
  11167178.125,
  11167421.875,
  11168056.25,
  11168268.75,
  11168671.875,
  11169267.1875,
  11169307.8125,
  11169353.125,
  11169476.5625,
  11169884.375,
  11170520.3125,
  11170590.625,
  11170656.25,
  11172698.4375,
  11173348.4375,
  11174301.5625,
  11176385.9375,
  11176690.625,
  11179804.6875,
  11180020.3125,
  11180125.0,
  11181593.75,
  11190251.5625],
 [5.152539720628002,
  6.854280608202352,
  28.46977324915323,
  8.518592729869086,
  104.12549718709128,
  10.851741588919381,
  5.517981870959336,
  66.07082341469258,
  5.07205040950528,
  70.83343627793475,
  13.477841232271121,
  21.514784216537013,
  7.294620685653026,
  71.59663946091045,
  19.04167810178794,
  101.8615281647461,
  38.31463262727778,
  15.748185106962069,
  22.337371770480527,
  41.86135388760191,
  24.76487392530209,
  42.20275329711578,
  22.83905188675709,
  64.1082339267195,
  26.240131209688553,
  59.809587418109714,
  44.944165232861174,
  39.24467288238996,
  90.35129801094376,
  30.49577284348799,
  120.38109198228508,
  11.191893569629812,
  60.31832820016083,
  28.40224654850701,
  11.554734225267199,
  40.18066455224005,
  36.20438031671707,
  28.351095269272477,
  57.1538165009531,
  51.573978550881286,
  87.88166151492781,
  36.69209596588307,
  13.675773006736271,
  28.839588101807237,
  55.04759646551981,
  5.687725581577804,
  20.92132427970974,
  61.755145281881425,
  45.54590494190182,
  66.63260202559736,
  18.41951993227862,
  9.361202130713966,
  41.85733994915828,
  40.41770608476546,
  77.83344741274658,
  33.50197557675041,
  15.768944741324155,
  67.71139077171931,
  84.67955027596174,
  6.011356114618717,
  6.125334205525346,
  92.5483191641962,
  135.5897747452631,
  5.960620084817627,
  40.39701737275518,
  10.540728353701452,
  39.45172103609671,
  76.39916273657171,
  63.80271822176243,
  5.790528280843724,
  7.200659289731451,
  11.509572524544506,
  80.01810437947763,
  5.171117644277947,
  96.27683852930934,
  186.1417397460132,
  8.458468574236372,
  13.888748650162718,
  28.979197039357043,
  14.203136917788044,
  19.147186778991607,
  20.71329242183629,
  25.29131804774221,
  5.413523852409813,
  130.16779268193548,
  34.235723812894264,
  7.014783924191486,
  67.48605592880169,
  171.36108906068372,
  34.293081690683756,
  25.59189536940345,
  5.875564676840401,
  130.5400459512419,
  180.8139704964108,
  24.40088221057529,
  30.664071638417916,
  49.29916302104552,
  19.211923297463933,
  38.35488307977782,
  7.194868237328745,
  23.600278627113184,
  10.923813758005368,
  47.607297203554104,
  13.970180127193217,
  5.32309791551664,
  109.8117977837779,
  18.285130042481654,
  5.542866385732408,
  43.80287967358354,
  14.510289209141238,
  22.583907452686088,
  24.847664972428685,
  33.82653267657091,
  21.16403345986571,
  171.12548785117116,
  116.9795915399358,
  66.17091062104805,
  9.443723749831062,
  12.087585797711258,
  31.973466546309027,
  16.435148161313535,
  41.212228278081646,
  15.900425747520845,
  20.733557048683927,
  25.197709162988858,
  65.82600010526188,
  20.129005151250173,
  11.749626610644865,
  139.47939737716968,
  68.48859484153373,
  24.977419719971078,
  14.152626115950188,
  14.780748339516347,
  25.40293149882287,
  76.51991348766333,
  17.54669692670169,
  19.305919145468287,
  60.40740473952349,
  7.2109320187492,
  39.20246516494882,
  48.269934970371956,
  70.74524123177379,
  56.14594605387755,
  11.086447050420286,
  45.2014525582852,
  25.559896113786696,
  25.11342373677851,
  17.169847986045625,
  79.73208922562054,
  52.590421189295434,
  36.095362761543846,
  5.346624621020493,
  42.44837923419093,
  87.28161169490454,
  25.881576481000494,
  8.137587074235956,
  53.190537051731575,
  17.887233167459463,
  31.59946636889701,
  24.074466485893478,
  19.496495567196902,
  26.346880458735214,
  15.315357157168105,
  49.99310863578063,
  80.9253123742679,
  9.473031339365475,
  22.486791936312983,
  11.662643518001484,
  71.85230768789084,
  16.735752443184825,
  103.20065404595042,
  7.586015517187785,
  20.600063522329673,
  5.960674491227817,
  85.55176226011383,
  75.01433825157935,
  20.08003280395637,
  12.648045863603514,
  9.51613277849609,
  7.93866252653696,
  86.82205887681428,
  72.23923501461027,
  10.401945538605132,
  7.95739135679953,
  16.391244436923426,
  40.535221212822655,
  34.272901388006744,
  54.51319317104874,
  10.935181603713145,
  41.68904159387064,
  36.64664766442374,
  60.270084587438525,
  51.54528200840906,
  41.04865140503001,
  18.521115835466233,
  77.97960074148779,
  42.40383516276643,
  8.043868596005852,
  18.45246269732988,
  49.958918820606584,
  56.28467471045701,
  81.4668970979767,
  25.588125516294543,
  88.46723420294437,
  49.20938784171525,
  29.492253688002982,
  26.113114805826186,
  14.513556031106724,
  46.422392689496654,
  10.1746537804279,
  11.483746364289411,
  7.318606183786745,
  6.346254242713159,
  22.096643556250786,
  23.224134371512825,
  44.04413327415761,
  18.522602832307346,
  29.05964972547983,
  22.509880825783704,
  7.700183749296881,
  65.02246401930786,
  68.74819988327873,
  43.97472934340662,
  108.84548564690715,
  7.320219495069112,
  6.2761340964351415,
  13.336095984452928,
  55.226124766388004,
  45.175026213035906,
  114.73517871425935,
  15.995647993758702,
  6.930998450410038,
  5.031119491253183,
  7.396871343343425,
  6.767683151708588,
  155.31080594303745,
  27.818023408680585,
  21.12701594558999,
  29.66940350703561,
  5.137249344001191,
  50.33604985175072,
  71.5731356473262,
  18.301471469860463,
  20.948050266697503,
  100.21797413158507,
  31.743090192893078,
  5.67603872794253,
  54.764653430482234,
  70.81253963895017,
  28.084685354651697,
  86.83339496681926,
  20.704348569831264,
  106.73138931339427,
  8.623631646227464,
  57.041228718383074,
  11.20262373289785,
  43.64437065664903,
  82.90455393049896,
  25.242442820201994,
  16.884654122020006,
  84.29590427404828,
  15.943988677524754,
  123.91754786059644,
  20.105722669737105,
  23.761109782440226,
  5.153674506095576,
  34.220916806604734,
  12.41614136873189,
  5.473770256582496,
  10.890228747204613,
  28.15591359082218,
  8.062699944680636,
  20.74396740876828,
  128.9700307914263,
  58.59504639904654,
  48.804394379641764,
  21.95333066002647,
  11.1168484396999,
  15.1378607061571,
  15.434417533188464,
  40.7726903615389,
  11.46923288239749,
  41.89678193905154,
  122.78089239035424,
  19.92304984647383,
  15.66076906866832,
  27.883494730103067,
  102.75934387501144,
  36.012735343213365,
  22.691350507312702,
  27.252866080233144,
  15.277988413481557,
  20.665325164727093,
  59.37677793786111,
  11.966486557685885,
  8.122196823756076,
  5.041232477753647,
  10.780539126984296,
  16.895224306132164,
  20.48843869720723,
  7.578437925359917,
  19.297314369770156,
  31.88471947070518,
  35.19837571782373,
  54.87603439547455,
  25.230997811952303,
  71.91030404249504,
  13.597968080863915,
  21.952513914708014,
  19.449292463114418,
  8.341334518937044,
  59.999329074837476,
  11.38497874207974,
  6.961264139888034,
  13.00771208931076,
  7.1132072541006135,
  16.811420378171682,
  77.04137400714121,
  5.533212121326684,
  15.391321861575165,
  6.268804524413115,
  24.596142735302863,
  40.43658418948557,
  55.30297277765685,
  16.760382768928178,
  16.519302034079562,
  12.27867586412826,
  21.291185000968785,
  51.16257862654765,
  14.31632160919657,
  6.557504010903474,
  66.64742071867056,
  46.633192304423915,
  11.743559607178831,
  25.314954189756374,
  31.162060904972744,
  5.211515439271813,
  10.378713614348863,
  98.00176264788041,
  16.907295975214236,
  5.472559157478296,
  14.272425737459907,
  88.81167535663668,
  11.240867233831823,
  42.02708852416934,
  5.585498339324491,
  15.812385753677258,
  40.09483634410534,
  40.45859063594583,
  15.229290694036886,
  30.59559926462969,
  25.98196220755276,
  38.58029503520994,
  20.379248927128046,
  19.986745785314824,
  16.305593850579523,
  11.164446095152494,
  36.172450538249336,
  67.67295146877589,
  116.41649689985026,
  18.741962558021726,
  19.083042167773066,
  10.4630148627016,
  67.72030521534391,
  59.82808442854977,
  11.172606097271364,
  70.31555036349765,
  8.552947939283497,
  19.20767016588948,
  30.814969696915252,
  10.92591029708916,
  17.967127742497667,
  58.06875056584753,
  36.29990209987562,
  24.249659806568168,
  14.471739836965908,
  36.67635838868224,
  8.256167712041293,
  21.933269233382674,
  23.72130813587659,
  66.6259338899255,
  20.635910129127666,
  22.379783319646005,
  12.681214690936166,
  8.529439566140935,
  62.34803907660725,
  6.675531474630341,
  35.633463154170656,
  9.814290307496611,
  95.886496185331,
  15.590193431898967,
  65.90200675322568,
  20.4587117178921,
  25.627230975421504,
  13.064130728992723,
  43.20720915324519,
  31.278666802864652,
  45.33777963589865,
  46.146959600607595,
  18.956521288079514,
  7.156548773757226,
  8.571271260729043,
  46.92798890070423,
  5.682290111987992,
  18.40784042348468,
  12.131004937511317,
  6.373642381055459,
  21.931660196806877,
  20.065618757236813,
  63.510783073412284,
  22.920092270423453,
  112.33175443260717,
  13.27371108081708,
  5.716294736428011,
  12.81975667709802,
  70.82363526510603,
  11.28756711526393,
  23.183950803297233,
  33.732320090633976,
  63.94673550935585,
  14.183702210261345,
  5.2953254522313955,
  60.39336210447182,
  34.452494746811595,
  11.396605876286944,
  11.651889344045538,
  22.701832172064627,
  14.452035460802747,
  17.66267803526821,
  6.505809198579629,
  16.48959972276387,
  66.65388061844983,
  14.955776192994849,
  33.92185951219833,
  17.86996744696716,
  24.474558195116117,
  59.26694594642727,
  45.56915405824121,
  8.595740928366824,
  19.037648175269272,
  25.35777756700905,
  45.40384899156447,
  11.693675053578785,
  35.209035827186035,
  23.639042189965767,
  15.505203043197842,
  68.66356314370131,
  18.19565063851373,
  16.82124653958557,
  23.362626243755667,
  14.87575813333602,
  72.94167640088946,
  22.06683624877233,
  13.006993305602885,
  29.41842898033386,
  23.471531735608902,
  75.45367467594446,
  83.88596928035406,
  20.166283026243185,
  6.05005583064838,
  9.140306808718668,
  25.884140793399066,
  9.258400598827512,
  6.993925299398916,
  18.610203750433183,
  59.113875314393475,
  28.1253274360485,
  11.548869799180716,
  67.55106370241315,
  55.9881592055243,
  34.817506184456235,
  6.97554827512385,
  26.026504858549526,
  10.662979206593354,
  35.85399359675375,
  21.505246454104626,
  26.99636480989424,
  5.567337491295738,
  18.18712707321843,
  44.05514089594186,
  13.578008233895376,
  31.509836926867912,
  9.129941566361667,
  5.889746723795036,
  16.642400753730172,
  15.18501496634751,
  27.0943252053947,
  19.076704642784655,
  6.691163229226904,
  57.66642499346492,
  16.615214856301986,
  5.323664037217394,
  12.870639758178005,
  13.47855306367875,
  84.86008029315596,
  15.024038021536226,
  21.08832857150172,
  20.502913106869094,
  33.75538405899368,
  6.733575600312667,
  15.809364813691337,
  7.201286681865446,
  6.0816451145323285,
  79.26252707616494,
  13.726517039074029,
  47.781471066476385,
  22.24930516414614,
  60.7518547946961,
  31.221493428716915,
  24.966690284075703,
  14.347249701326785,
  21.599194252546575,
  38.3588353782252,
  15.564421283043233,
  13.597774954942903,
  39.99628468780775,
  22.88995648003838,
  40.48202396444261,
  35.67166497409103,
  134.39385721480718,
  5.300263663531363,
  85.30866242084073,
  69.29067386173162,
  41.64134077449923,
  43.43278645042357,
  59.87782170498564,
  55.78212196213634,
  6.899441414885073,
  14.38002948266162,
  38.737229854579226,
  19.287070746067666,
  19.365617318460707,
  7.858825854702455,
  6.779570110032547,
  8.623146387587262,
  9.429939411889002,
  5.3678137539094415,
  16.7341984319211,
  22.73427983825978,
  11.966264598200269,
  6.834493365791316,
  5.379952966447907,
  39.778012427051905,
  26.206673266390045,
  12.036368094820142,
  37.99691192933244,
  49.06396024947829,
  146.17828993538492,
  55.25260988315125,
  37.1440217974349,
  15.787236319903396,
  16.647127300143982,
  6.12673853908567,
  5.659348981629775,
  22.30105131438771,
  16.729821211280854,
  29.051625317992823,
  10.811494287935094,
  25.861532066629692,
  34.97851508179406,
  67.02990443086271,
  24.414695995324138,
  16.379813581841642,
  23.11032814379475,
  55.247479014106645,
  61.38022337818151,
  24.324547752369202,
  22.241070093804975,
  7.314593122274525,
  80.20173111609554,
  77.0908019147634,
  102.9154672109412,
  22.957495580830614,
  33.765342634374036,
  46.24134287797399,
  69.1674379915813,
  6.375937003992957,
  8.713720782090602,
  100.05871615653862,
  28.610017911759115,
  79.28598522153132,
  16.865336857350403,
  64.55196103296775,
  9.53284924231757,
  19.443259794927265,
  5.033671761238463,
  10.470431196475305,
  25.097019873028934,
  35.04574341438455,
  27.402774751803868,
  12.694683810796962,
  8.340559177073482,
  13.640831375927062,
  54.81265703640567,
  49.06493623588928,
  5.62065897658044,
  8.794972639428673,
  87.87589858414329,
  30.56781556105499,
  47.260721922823905,
  78.42667838875936,
  48.34745105241726,
  22.619691068853466,
  68.67742711842382,
  6.318123902707775,
  5.5245489791167115,
  5.7359263785661305,
  52.12040622413002,
  11.50279764375112,
  57.866386987602226,
  42.608079697869286,
  85.535365369179,
  10.253058998940517,
  13.472752065084217,
  42.71999881154682,
  24.339973270093914,
  6.913804574393271,
  23.523027679021425,
  79.90497877345776,
  6.785036547183423,
  30.94330828229974,
  56.959207546751344,
  104.11825331105041,
  20.765798431667406,
  6.6513033600003375,
  47.76379734947511,
  37.89704300608452,
  8.56893262957653,
  7.885937528121809,
  26.27196888419597,
  26.26207190317276,
  51.051474998370665,
  14.108970344869345,
  6.773956792795996,
  33.69043216476397,
  92.38500111882406,
  40.292837394265725,
  35.630836819228385,
  74.14495319564793,
  23.83770323817886,
  8.822136888904623,
  13.207788281550927,
  101.69144236874027,
  50.83514608014728,
  11.361083673417486,
  105.45912645222798,
  17.43275333135287,
  10.799588486890835,
  26.664916885708884,
  30.04018034377117,
  11.969509083779881,
  26.214602777262158,
  47.738302108751135,
  26.654833699605224,
  53.26300719531444,
  36.57322617394971,
  11.079982486593522,
  19.746098796914005,
  13.11481075609913,
  6.431811899239111,
  47.500103614519844,
  28.433849245706867,
  5.328613322419715,
  11.888137753680944,
  16.759864257868692,
  9.6620400973615,
  57.540275366133756,
  18.81247992716082,
  35.09621044822148,
  7.865318277315491,
  28.762790559811386,
  33.41290685122884,
  5.432056116735115,
  10.007605498857345,
  6.830449368245077,
  37.41825647814323,
  15.793492646300152,
  7.9181812397333475,
  60.219999080497615,
  51.71645802030075,
  12.279336418927917,
  5.339399799810742,
  11.160836850946943,
  13.716416215455652,
  51.61698277085947,
  9.043224983718092,
  88.66804603884162,
  45.37083132518059,
  21.57801995895946,
  77.2429299500754,
  62.76121402777798,
  7.603611602313181,
  35.33944535508393,
  30.609880227178472,
  12.177713182778449,
  18.675053405733717,
  7.654187800790934,
  54.36381559419603,
  27.994434279752134,
  12.665632173913874,
  45.660785740555816,
  18.92654886503231,
  32.17850497185195,
  18.233160117141868,
  40.79318550543543,
  54.72864586884015,
  32.364631241112306,
  12.836469344741218,
  18.740895444816292,
  5.629464643133254,
  12.523195639411641,
  17.26606467230877,
  22.11263852921472,
  6.407171551870401,
  49.31818741691774,
  7.333911131171782,
  81.60236692441008,
  9.725901949970657,
  11.261091657772775,
  11.530471496658208,
  6.702385151807667,
  14.80563544414959,
  29.65868461440448,
  7.804736296517852,
  11.551961936322066,
  99.99260820493305,
  24.514330279134825,
  17.07862801704698,
  42.52883438094541,
  5.978605888754163,
  75.74311853146145,
  97.90859183660349,
  56.948107276915316,
  67.48149726836735,
  33.42602915415409,
  33.50280221714644,
  21.360916216218968,
  41.48654089270738,
  5.896252938555228,
  81.94242824717769,
  77.20691742822233,
  10.152184116461445,
  13.017087975525696,
  23.577173353366046,
  109.27527899189722,
  114.18203835613429,
  50.30175731612823,
  5.146128300237676,
  12.865746802158746,
  33.40661388393106,
  69.20976579662678,
  53.32857712598188,
  81.99181316987061,
  91.74674516924516,
  141.83846404564252,
  15.66922277029354,
  41.07080556221892,
  69.69847771686334])

Graph for home page

Detected energies during the discharge + Energy spectrum

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