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

([3587735.9375,
  6802431.25,
  6842682.8125,
  6920104.6875,
  6949962.5,
  6950360.9375,
  7015684.375,
  7017779.6875,
  7106253.125,
  7188793.75,
  7195482.8125,
  7195539.0625,
  7196743.75,
  7197140.625,
  7275157.8125,
  7382965.625,
  7386604.6875,
  7390162.5,
  7391800.0,
  7391832.8125,
  7415934.375,
  7424310.9375,
  7468432.8125,
  7468512.5,
  7469973.4375,
  7470921.875,
  7473642.1875,
  7473678.125,
  7475582.8125,
  7479473.4375,
  7498467.1875,
  7503434.375,
  7506903.125,
  7508653.125,
  7508834.375,
  7513351.5625,
  7513446.875,
  7514012.5,
  7516193.75,
  7516782.8125,
  7517039.0625,
  7517884.375,
  7526164.0625,
  7531517.1875,
  7532157.8125,
  7533625.0,
  7550404.6875,
  7550453.125,
  7552068.75,
  7553210.9375,
  7553343.75,
  7553731.25,
  7591914.0625,
  7592870.3125,
  7599943.75,
  7600496.875,
  7612810.9375,
  7623751.5625,
  7626632.8125,
  7626650.0,
  7628281.25,
  7639168.75,
  7641737.5,
  7642992.1875,
  7643592.1875,
  7644710.9375,
  7661100.0,
  7668985.9375,
  7669834.375,
  7683825.0,
  7694121.875,
  7704985.9375,
  7706587.5,
  7712809.375,
  7712868.75,
  7718031.25,
  7718665.625,
  7723415.625,
  7729760.9375,
  7739356.25,
  7745181.25,
  7746237.5,
  7748459.375,
  7752043.75,
  7760135.9375,
  7766709.375,
  7784137.5,
  7813943.75,
  7817521.875,
  7822085.9375,
  7855440.625,
  7857710.9375,
  7860057.8125,
  7860148.4375,
  7861342.1875,
  7861828.125,
  7861918.75,
  7863017.1875,
  7863412.5,
  7863564.0625,
  7863617.1875,
  7864025.0,
  7864389.0625,
  7868554.6875,
  7871267.1875,
  7876717.1875,
  7877762.5,
  7878728.125,
  7880871.875,
  7889887.5,
  7917645.3125,
  7934773.4375,
  7934795.3125,
  7943190.625,
  7967570.3125,
  7969285.9375,
  7971932.8125,
  7974818.75,
  7987201.5625,
  7989114.0625,
  7990282.8125,
  7991693.75,
  7993553.125,
  7995120.3125,
  7995895.3125,
  8009334.375,
  8013240.625,
  8020110.9375,
  8020296.875,
  8021992.1875,
  8040154.6875,
  8042485.9375,
  8044375.0,
  8045057.8125,
  8045390.625,
  8046200.0,
  8046748.4375,
  8046976.5625,
  8047296.875,
  8047625.0,
  8047684.375,
  8047728.125,
  8047748.4375,
  8048518.75,
  8048593.75,
  8048803.125,
  8048856.25,
  8049306.25,
  8049317.1875,
  8049546.875,
  8050032.8125,
  8050335.9375,
  8050664.0625,
  8050670.3125,
  8051039.0625,
  8051062.5,
  8051306.25,
  8051645.3125,
  8051739.0625,
  8052314.0625,
  8054231.25,
  8056465.625,
  8057725.0,
  8058593.75,
  8059387.5,
  8073864.0625,
  8079104.6875,
  8080493.75,
  8083139.0625,
  8083304.6875,
  8087073.4375,
  8087195.3125,
  8087335.9375,
  8087360.9375,
  8087509.375,
  8088851.5625,
  8091026.5625,
  8092528.125,
  8093196.875,
  8093304.6875,
  8094446.875,
  8118250.0,
  8118320.3125,
  8119912.5,
  8121423.4375,
  8121706.25,
  8121717.1875,
  8126790.625,
  8147593.75,
  8150389.0625,
  8153126.5625,
  8159525.0,
  8161943.75,
  8162251.5625,
  8164865.625,
  8167559.375,
  8168656.25,
  8176117.1875,
  8180609.375,
  8180653.125,
  8192901.5625,
  8203507.8125,
  8217309.375,
  8229478.125,
  8229889.0625,
  8230914.0625,
  8234614.0625,
  8235764.0625,
  8237190.625,
  8237760.9375,
  8245637.5,
  8266554.6875,
  8269339.0625,
  8272037.5,
  8275417.1875,
  8276156.25,
  8277659.375,
  8279089.0625,
  8280885.9375,
  8299728.125,
  8302331.25,
  8302506.25,
  8304726.5625,
  8305854.6875,
  8305926.5625,
  8308570.3125,
  8308771.875,
  8315718.75,
  8315821.875,
  8316964.0625,
  8338843.75,
  8345293.75,
  8347054.6875,
  8351057.8125,
  8377525.0,
  8380632.8125,
  8384343.75,
  8384429.6875,
  8385178.125,
  8388120.3125,
  8389101.5625,
  8393829.6875,
  8394310.9375,
  8395693.75,
  8399351.5625,
  8401098.4375,
  8409232.8125,
  8410437.5,
  8411906.25,
  8412567.1875,
  8412601.5625,
  8412737.5,
  8412912.5,
  8413031.25,
  8415842.1875,
  8417812.5,
  8418051.5625,
  8418064.0625,
  8418073.4375,
  8418089.0625,
  8418221.875,
  8422181.25,
  8428195.3125,
  8433054.6875,
  8434867.1875,
  8439981.25,
  8451942.1875,
  8455096.875,
  8455562.5,
  8462157.8125,
  8464734.375,
  8465123.4375,
  8468764.0625,
  8468856.25,
  8469284.375,
  8471246.875,
  8471862.5,
  8474565.625,
  8483171.875,
  8485726.5625,
  8485739.0625,
  8491293.75,
  8491321.875,
  8501990.625,
  8504814.0625,
  8506657.8125,
  8507259.375,
  8507540.625,
  8508042.1875,
  8508593.75,
  8509446.875,
  8510001.5625,
  8510521.875,
  8510620.3125,
  8510718.75,
  8510846.875,
  8513764.0625,
  8514793.75,
  8514920.3125,
  8515700.0,
  8516482.8125,
  8517231.25,
  8519948.4375,
  8529168.75,
  8529276.5625,
  8537926.5625,
  8541981.25,
  8548757.8125,
  8622457.8125,
  8636645.3125,
  8640654.6875,
  8643534.375,
  8644235.9375,
  8645557.8125,
  8645826.5625,
  8649232.8125,
  8650785.9375,
  8651528.125,
  8662882.8125,
  8662923.4375,
  8667153.125,
  8669218.75,
  8669221.875,
  8670859.375,
  8674385.9375,
  8679245.3125,
  8680443.75,
  8709396.875,
  8727407.8125,
  8737837.5,
  8739473.4375,
  8740848.4375,
  8747032.8125,
  8792084.375,
  8838993.75,
  8861401.5625,
  8864907.8125,
  8875876.5625,
  8876231.25,
  8877085.9375,
  8898010.9375,
  8907914.0625,
  8912482.8125,
  8912593.75,
  8912981.25,
  8913306.25,
  8913618.75,
  8914153.125,
  8914568.75,
  8914834.375,
  8914871.875,
  8914982.8125,
  8915453.125,
  8916064.0625,
  8916609.375,
  8918048.4375,
  8918089.0625,
  8920103.125,
  8928975.0,
  8945767.1875,
  8954939.0625,
  8956365.625,
  8958364.0625,
  8958376.5625,
  8960845.3125,
  8962526.5625,
  9010131.25,
  9010667.1875,
  9011198.4375,
  9011273.4375,
  9012182.8125,
  9014937.5,
  9014951.5625,
  9015003.125,
  9018128.125,
  9023029.6875,
  9044710.9375,
  9045473.4375,
  9045540.625,
  9048034.375,
  9094550.0,
  9103175.0,
  9103375.0,
  9107237.5,
  9114895.3125,
  9115864.0625,
  9115937.5,
  9116026.5625,
  9116031.25,
  9117192.1875,
  9117276.5625,
  9117504.6875,
  9118426.5625,
  9119031.25,
  9119481.25,
  9119520.3125,
  9120359.375,
  9121345.3125,
  9125962.5,
  9139135.9375,
  9139292.1875,
  9139881.25,
  9140506.25,
  9140873.4375,
  9144987.5,
  9145659.375,
  9155585.9375,
  9198287.5,
  9228348.4375,
  9249368.75,
  9250046.875,
  9251948.4375,
  9252273.4375,
  9257131.25,
  9268209.375,
  9280120.3125,
  9282406.25,
  9288718.75,
  9297796.875,
  9349317.1875,
  9370676.5625,
  9391842.1875,
  9402478.125,
  9402481.25,
  9415321.875,
  9421348.4375,
  9421409.375,
  9421640.625,
  9421818.75,
  9422517.1875,
  9422782.8125,
  9422992.1875,
  9454517.1875,
  9457410.9375,
  9457925.0,
  9458476.5625,
  9458532.8125,
  9458703.125,
  9459209.375,
  9550031.25,
  9558137.5,
  9562696.875,
  9562915.625,
  9563056.25,
  9563235.9375,
  9563315.625,
  9563631.25,
  9563784.375,
  9565028.125,
  9565415.625,
  9596843.75,
  9599223.4375,
  9600146.875,
  9601960.9375,
  9605525.0,
  9606675.0,
  9634103.125,
  9634468.75,
  9634960.9375,
  9635785.9375,
  9635796.875,
  9636626.5625,
  9636768.75,
  9637006.25,
  9637120.3125,
  9637945.3125,
  9641653.125,
  9643090.625,
  9643400.0,
  9661076.5625,
  9667389.0625,
  9673725.0,
  9697450.0,
  9739523.4375,
  9739548.4375,
  9770415.625,
  9788409.375,
  10077876.5625,
  10325946.875,
  10330818.75,
  10331992.1875,
  10371321.875,
  10398273.4375,
  10410256.25,
  10418064.0625,
  10418176.5625,
  10419162.5,
  10420442.1875,
  10425970.3125,
  10427573.4375,
  10432634.375,
  10447400.0,
  10447464.0625,
  10448446.875,
  10448559.375,
  10450784.375,
  10463826.5625,
  10464318.75,
  10470898.4375,
  10514060.9375,
  10515526.5625,
  10516678.125,
  10516709.375,
  10522045.3125,
  10639671.875,
  10731839.0625,
  10740009.375,
  10740075.0,
  10768415.625,
  10775800.0,
  10795846.875,
  10827765.625,
  10827873.4375,
  10833659.375,
  11007678.125,
  11009654.6875,
  11011651.5625,
  11027709.375,
  11029242.1875,
  11031825.0,
  11037675.0,
  11040745.3125,
  11041187.5,
  11041401.5625,
  11042421.875,
  11042721.875,
  11048273.4375,
  11049246.875,
  11051325.0,
  11051628.125,
  11052604.6875,
  11052657.8125,
  11052859.375,
  11052892.1875,
  11053489.0625,
  11060757.8125,
  11062001.5625,
  11066259.375,
  11067946.875,
  11069446.875,
  11070087.5,
  11070993.75,
  11071000.0,
  11148853.125,
  11155615.625,
  11158576.5625,
  11158876.5625,
  11163339.0625,
  11163839.0625,
  11165573.4375,
  11168075.0,
  11168362.5,
  11168370.3125,
  11170203.125,
  11178767.1875,
  11182475.0,
  11185226.5625,
  11185640.625,
  11191065.625,
  11196364.0625,
  11205457.8125,
  11213829.6875,
  11220510.9375,
  11223576.5625,
  11229248.4375,
  11229265.625,
  11240378.125,
  11257943.75,
  11278460.9375,
  11282954.6875,
  11286684.375,
  11296714.0625,
  11297750.0,
  11299289.0625,
  11300670.3125,
  11301668.75,
  11301767.1875,
  11302646.875,
  11303221.875,
  11308303.125,
  11308631.25,
  11308729.6875,
  11309160.9375,
  11310640.625,
  11311140.625,
  11312531.25,
  11313229.6875,
  11313371.875,
  11317139.0625,
  11317221.875,
  11318100.0,
  11318531.25,
  11318557.8125,
  11318907.8125,
  11318992.1875,
  11319001.5625,
  11319034.375,
  11319064.0625,
  11319132.8125,
  11319290.625,
  11319296.875,
  11319523.4375,
  11320175.0,
  11320187.5,
  11320196.875,
  11320200.0,
  11320207.8125,
  11320307.8125,
  11320379.6875,
  11320445.3125,
  11320484.375,
  11320534.375,
  11320609.375,
  11320632.8125,
  11320685.9375,
  11320806.25,
  11320835.9375,
  11320859.375,
  11320948.4375,
  11321134.375,
  11321195.3125,
  11321317.1875,
  11321325.0,
  11321334.375,
  11321375.0,
  11321445.3125,
  11321451.5625,
  11321548.4375,
  11321631.25,
  11321645.3125,
  11321650.0,
  11321651.5625,
  11321670.3125,
  11321673.4375,
  11321692.1875,
  11321706.25,
  11321707.8125,
  11321787.5,
  11321867.1875,
  11322164.0625,
  11322382.8125,
  11322437.5,
  11322495.3125,
  11324268.75,
  11324734.375,
  11325335.9375,
  11326656.25,
  11326773.4375,
  11327664.0625,
  11328546.875,
  11328968.75,
  11329412.5,
  11330157.8125,
  11330164.0625,
  11335054.6875,
  11336000.0,
  11339350.0,
  11348604.6875],
 [43.75392486847537,
  53.30134253245989,
  20.618502087399357,
  21.924406455376186,
  5.998002358215999,
  14.786688585495128,
  22.093196635757344,
  23.68024596119375,
  49.880642062778556,
  28.545659552558575,
  49.97113678068086,
  5.621115598902164,
  25.40780921453212,
  54.12124961084828,
  115.06108824214039,
  38.02860544665343,
  43.207414702213654,
  128.4824300140363,
  21.715094837566305,
  5.3886399917119405,
  42.922469130228535,
  26.738895845935634,
  25.530083219263126,
  26.733470676399556,
  57.11116651942968,
  17.387852869645496,
  16.529281526471504,
  5.782203723457968,
  35.539905034706365,
  110.89658736604568,
  5.2227790612697245,
  61.0837397063755,
  101.81588581564756,
  5.859871244315133,
  53.88475944792805,
  40.75024709223862,
  52.34716024748916,
  28.739353136881753,
  23.01820822849388,
  25.00720641955666,
  8.78088466984392,
  7.6218485553600335,
  7.174991517234949,
  16.113131589403842,
  51.131353566608624,
  61.233595888689834,
  25.45036325466539,
  8.072804957241642,
  79.32461414853991,
  8.49243355320877,
  12.96251965794333,
  7.078361901396194,
  23.605325737183904,
  46.74114375529501,
  77.12116530033984,
  29.40745167985136,
  158.35104411458886,
  8.418614803314128,
  26.015051982204774,
  55.56095569895757,
  67.79174536866789,
  7.282727483688339,
  5.306245525254198,
  39.616274915319664,
  38.365673509787314,
  15.168393593885424,
  103.00719759890382,
  30.59400773529808,
  60.46513376590568,
  13.59523915647309,
  13.98442482326472,
  18.83340320417493,
  92.45700577003812,
  164.7597848344842,
  6.2216284263079205,
  21.862049316945615,
  47.81199286191655,
  26.235870242488485,
  21.918383834820453,
  109.50129942928082,
  56.85705873108824,
  7.266005216878156,
  74.53906579978519,
  165.2228845648882,
  170.1498074003878,
  119.40422156835832,
  31.404124426836177,
  7.804210999045414,
  280.113029069459,
  22.614670440139495,
  67.94283851166814,
  64.37342961638569,
  21.266848516377248,
  31.599589344211594,
  20.59337399075316,
  12.87531967679679,
  5.061777665303518,
  16.1301893243049,
  5.216226737409376,
  27.132838893759356,
  11.593850624895591,
  81.83205461453588,
  9.913849821767581,
  41.86246897757699,
  73.87737571503175,
  19.69082384754751,
  127.76211110181643,
  5.8189128719845415,
  37.75297030915435,
  5.140248267368501,
  57.13899022260139,
  42.2129512169315,
  21.936811292404215,
  6.995837870222363,
  8.3438645572348,
  21.25739918940065,
  15.872419748789111,
  40.23707435461132,
  5.574288905078166,
  5.945044658998763,
  14.085976750598574,
  67.19019410745126,
  17.97092346727778,
  10.991760657873488,
  24.460192122423337,
  12.84821080003411,
  32.016047783948856,
  7.850894185283803,
  5.160396058089147,
  31.941227597115116,
  95.43425153435207,
  42.631870892468044,
  32.91635765563572,
  21.66045169580405,
  10.508535120488729,
  5.28306949989141,
  54.145080976936484,
  47.33039298365725,
  11.86945278746586,
  16.54348374789551,
  125.31657824084525,
  22.186633516588397,
  5.091154530554696,
  64.75589985183423,
  22.440970856910052,
  45.305497934321124,
  90.15879114157494,
  25.134685854975107,
  5.1347737761070285,
  13.687082883093103,
  96.36901726754965,
  20.16411602819774,
  56.94311118831999,
  11.492030932702049,
  114.24098958373553,
  33.54265949693855,
  9.931351135515833,
  5.887329859625653,
  36.7534253499761,
  15.735713464661806,
  40.83553411376458,
  24.774056696480542,
  53.564450313534635,
  107.53647075017517,
  29.949173515979236,
  5.2062649485443755,
  66.82711371818307,
  34.971916601410264,
  25.569493017896214,
  42.980025219009576,
  30.90573732860115,
  20.396730799438295,
  13.547814446505033,
  5.915621247658618,
  36.08713654857849,
  13.06910721637352,
  7.505004826910014,
  244.05993249000988,
  49.8579705444769,
  45.25762938516009,
  65.40244533665721,
  32.9014808539647,
  84.9042385848715,
  11.661303724722579,
  14.590991437742778,
  84.3619086805655,
  22.92885562569729,
  97.98584093477328,
  74.4383668533196,
  76.68229294493757,
  11.539071769791505,
  8.645717856359637,
  19.895290604577436,
  23.723269455149136,
  27.112570679417374,
  11.879210610849176,
  11.537239648147086,
  36.91157897781484,
  25.16423341634477,
  35.88155545615587,
  25.513483900735274,
  12.748184764742398,
  13.394873733489309,
  69.29429794685359,
  26.28281426419346,
  187.00466972933881,
  7.9704995115765795,
  173.67092396107054,
  76.10002814779887,
  91.88516735600359,
  22.223260547878148,
  111.47409997418461,
  24.043690545754973,
  16.928984469290558,
  33.31780430922666,
  22.369811688294686,
  120.62874672209836,
  54.21027400403291,
  16.591875927195197,
  40.565105739391385,
  36.855526511018034,
  7.333567135139411,
  54.87351521662981,
  5.266270017674993,
  8.432971762130554,
  36.96469894214835,
  104.52967296577759,
  107.82788024407577,
  38.73531221892093,
  36.51860141397936,
  12.18437043220287,
  28.839289086033226,
  15.758084107941405,
  20.700839828265995,
  62.47732521400617,
  48.969913191618375,
  75.01117140744776,
  66.23846180165555,
  85.55648562731454,
  65.1686255242071,
  46.754776157216945,
  6.11904876995819,
  9.799004156789668,
  7.179359480964491,
  31.726751026044646,
  11.892663500080214,
  8.602995999892196,
  25.884921803466998,
  212.8452123476046,
  50.89203883811902,
  5.253993486741754,
  6.552720861958359,
  5.056294872556809,
  86.27666082539653,
  92.0682438066245,
  5.6980376374954504,
  20.141043453919497,
  18.985395531791937,
  70.13011038949608,
  46.851635287902795,
  29.890694816886814,
  11.06396829750044,
  21.244447579686643,
  18.794297835226832,
  81.8384197889961,
  17.11567347736093,
  7.4347814435829465,
  17.11543832697124,
  25.765596526440994,
  171.0189770425673,
  166.62417570530727,
  47.617306423538224,
  87.52555138620906,
  5.099150847659307,
  6.160764826158617,
  33.207044223580425,
  6.225943962290561,
  82.44689715444312,
  130.23094394657033,
  110.29107486064109,
  53.88273073600563,
  44.84042821853197,
  29.59419248412873,
  34.10380275093801,
  86.78925581246611,
  17.36267224747698,
  31.718896982946124,
  12.230208666697209,
  17.204902799252732,
  37.53597700431143,
  92.54196077707205,
  57.055982939784954,
  15.84891231011288,
  19.56022588511453,
  144.26705097138995,
  22.964994816741147,
  29.16119566437618,
  20.622435762921583,
  17.811022479678297,
  97.6487166632215,
  6.511078623653669,
  49.47645177092092,
  28.57188099561196,
  72.048764994024,
  144.93386168548756,
  6.323083148656236,
  15.111041166099659,
  42.031232681637526,
  12.223654060772887,
  5.526743297279815,
  53.08652689302342,
  8.433700076865025,
  110.42941834796935,
  26.322587517522503,
  24.174092846941424,
  78.4021757391126,
  5.74261000245582,
  8.376521969029175,
  17.43699187361873,
  23.101694817855012,
  90.07423512268645,
  48.59528686526562,
  43.772488222419405,
  43.365409914294055,
  21.38342085023234,
  46.31468531162235,
  46.18205156078226,
  84.8699528777467,
  46.05404738548869,
  43.229078007261236,
  30.87182995664722,
  8.436366247547305,
  80.71462110275678,
  10.700143446109799,
  44.734507441241824,
  12.692501109911207,
  6.260513558007149,
  10.122689190026074,
  10.618749267772554,
  40.06237374100719,
  45.14372333344614,
  9.479292995380307,
  88.99744685885068,
  13.781359629670543,
  98.34673049687999,
  60.388927354927034,
  50.22720107871789,
  53.5090445784896,
  83.06641751766367,
  62.60647755524819,
  85.55358745679399,
  30.56143379005863,
  66.2142546385347,
  33.292564619169084,
  5.330715179626472,
  102.88503052828759,
  13.001123588736194,
  23.680216149907416,
  22.350929885657884,
  6.286787999733525,
  52.28705619053464,
  27.914073367470174,
  128.19404468544138,
  48.235538759442974,
  6.730400570813883,
  99.04659770009894,
  18.566916446525155,
  21.590357604646172,
  8.364921840770798,
  5.97513148222563,
  99.36807730698254,
  26.94518336903817,
  31.953018426734953,
  16.47842063191503,
  62.33788944563073,
  53.75135712823922,
  9.585770462043971,
  111.38256205787181,
  5.0856831092483405,
  17.669869388756783,
  18.632841150474132,
  16.592382657074097,
  78.0692523966197,
  5.516222382514227,
  501.37446045827375,
  17.87286395669639,
  33.24100517722623,
  22.572729366525067,
  118.9539994110268,
  17.383326089971284,
  6.177673255921603,
  30.26536075125395,
  34.89182847919665,
  12.935428667451065,
  68.13655334229341,
  69.14336503373974,
  68.68534139384671,
  29.47256093982479,
  63.69146312283162,
  5.9153370570827,
  48.29949717437159,
  66.67555906606268,
  15.204573859599215,
  33.052763747043606,
  52.014793102059656,
  32.75085541957406,
  33.68430278458382,
  17.164569270988352,
  14.505985158509112,
  70.94980489885096,
  155.25802732296813,
  66.06698038312743,
  197.4950889017045,
  52.85158449322144,
  177.9833347556298,
  54.42913933269631,
  37.66726863017692,
  31.321499521215586,
  39.856671284012265,
  20.896251788192018,
  12.741620673160543,
  78.81400404884599,
  20.495216122662654,
  22.719183493831196,
  21.484453049263323,
  55.08032137088479,
  95.71536992761766,
  14.675238278772825,
  27.04244470207559,
  5.677331531976479,
  10.675411814870174,
  43.485536273772304,
  140.01531243048245,
  6.547540752207912,
  6.709089619481133,
  22.07409890641973,
  5.152222989628288,
  25.070068647948318,
  36.17991549075103,
  54.22198599884703,
  75.15239261103841,
  47.02716273158394,
  17.91447372116932,
  31.28091963514686,
  43.05488593947,
  72.03202750997454,
  92.81886229976021,
  156.09729156777183,
  25.655530301825515,
  137.32711929510913,
  9.952254347881473,
  27.025776806655237,
  70.6236706877964,
  87.92116177428154,
  37.18081653779992,
  26.71413628601404,
  9.983651149520439,
  24.4262667547115,
  11.210259423102478,
  42.257604953995205,
  76.56037085296416,
  7.0024741209140124,
  18.612094973680065,
  19.546798175785707,
  29.715450231436222,
  8.692935143158158,
  73.99708871565566,
  21.859953673229406,
  151.0012550239196,
  15.306668697106522,
  26.097428890835417,
  98.34282794889336,
  5.132050572252867,
  11.239138031428249,
  20.125041169697816,
  380.74803274445173,
  33.401300993992955,
  35.594815898855565,
  5.6208934925531455,
  50.378191367867316,
  8.217028233409224,
  10.358428794141663,
  58.6030095307516,
  56.281365905722204,
  38.33917593442391,
  28.269961222263774,
  20.39754991860405,
  9.484820167183592,
  16.070756741727923,
  34.98034773276001,
  65.88387423109285,
  171.43762260308569,
  43.700894668918345,
  31.577902806041813,
  99.68172293680057,
  9.418639423545025,
  25.109288934753867,
  19.33079836537231,
  77.03104227277372,
  23.847740967838693,
  65.22663213820249,
  40.28309944929495,
  153.0397994298229,
  17.740475801007346,
  13.361176381119877,
  76.22616542834456,
  239.31482866215651,
  5.07025553405524,
  8.619375273750922,
  46.861527377546466,
  57.48950066574868,
  44.80035046353777,
  5.468272250057785,
  6.022209334269816,
  9.742251182113371,
  23.346196272156877,
  68.11083438952073,
  93.66070164549868,
  18.887996635082647,
  17.65389205380559,
  7.4263953058941405,
  53.245769431770015,
  127.62136587696972,
  28.086955784993386,
  62.91661044120437,
  32.87485273663541,
  12.093148533363456,
  11.724106924459827,
  17.21892141179511,
  46.35176743779546,
  49.13404485571881,
  5.181754878040886,
  21.038563658034015,
  18.146352252419614,
  27.96253172721343,
  32.4923063076619,
  46.4530751552082,
  20.422302890360182,
  16.966380366088654,
  38.345068188325044,
  27.789861771085373,
  70.16112336982464,
  23.704345732819228,
  53.491457619753724,
  106.05260753675385,
  373.34633863600675,
  14.503795962194422,
  65.94213909246844,
  11.268594675892013,
  21.130121465968895,
  12.50130017743082,
  19.36663731199885,
  59.41211713093734,
  20.98000624233178,
  55.20440709162843,
  20.2643863777146,
  15.56041246628888,
  10.441582803986975,
  69.56215342613369,
  7.0968916243289,
  23.554188031337553,
  49.534255126767874,
  50.69777400029342,
  44.94298789934566,
  28.878549813415905,
  103.96705180896483,
  22.65670289940123,
  77.97790652855814,
  28.599325195318094,
  16.628217858213844,
  5.6061062657207845,
  12.919991332486576,
  10.495335778690041,
  17.053526752436763,
  16.342034641945663,
  6.9474412753616575,
  19.78051052151996,
  20.698785774917553,
  42.92041201862068,
  176.96140839714167,
  7.533339199150585,
  62.956313733053065,
  74.45810518930293,
  20.747312610269052,
  15.651905448000619,
  115.27646964125759,
  71.35468514857598,
  58.42951223027682,
  49.98773070886058,
  50.88248829304828,
  5.992713265028529,
  80.22818625145914,
  57.20046222754387,
  21.803440267492345,
  16.132821937515907,
  72.54960544831206,
  7.625923488230528,
  30.638369917801967,
  10.857609179556457,
  15.918213539817387,
  20.431950772799226,
  6.378697211437925,
  114.58964839656181,
  7.544904728209705,
  27.297615513093156,
  54.638669822191154,
  23.62640971925502,
  37.93452452601168,
  41.41350395965277,
  49.979079118055594,
  14.770197475017461,
  19.47726597568899,
  48.82832267206836,
  68.87290089420236,
  62.40582323431573,
  30.597380181581805,
  18.503564931578502,
  64.67488627528273,
  40.31119269849788,
  26.020123432969818,
  18.897969795126812,
  28.49062108291132,
  76.82417461303996,
  49.864464652641814,
  43.86366796467256,
  38.08146936614241,
  5.273975770415842,
  27.998411736108974,
  12.289328547466589,
  33.91695098318348,
  16.864499712731543,
  12.465866299430061,
  37.470479294181175,
  54.15075259243452,
  5.078880195492166,
  103.72940872669705,
  6.645027463501742,
  7.420427532244959,
  17.407716903590135,
  10.861476476316085,
  8.52815305330454,
  5.783638738817864,
  26.57606061557894,
  95.45959219224747,
  80.39329081427108,
  122.72202930155345,
  18.150249471841473,
  5.716106541369781,
  17.43718729812661,
  41.30287211859103,
  75.30350851567653,
  8.267309959365393,
  10.312053380372127,
  5.144619479335722,
  65.9520710769438,
  9.02487398823923,
  57.56725091644046,
  103.4028634400867])

Energy spectrum

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

([3587735.9375,
  6802431.25,
  6842682.8125,
  6920104.6875,
  6949962.5,
  6950360.9375,
  7015684.375,
  7017779.6875,
  7106253.125,
  7188793.75,
  7195482.8125,
  7195539.0625,
  7196743.75,
  7197140.625,
  7275157.8125,
  7382965.625,
  7386604.6875,
  7390162.5,
  7391800.0,
  7391832.8125,
  7415934.375,
  7424310.9375,
  7468432.8125,
  7468512.5,
  7469973.4375,
  7470921.875,
  7473642.1875,
  7473678.125,
  7475582.8125,
  7479473.4375,
  7498467.1875,
  7503434.375,
  7506903.125,
  7508653.125,
  7508834.375,
  7513351.5625,
  7513446.875,
  7514012.5,
  7516193.75,
  7516782.8125,
  7517039.0625,
  7517884.375,
  7526164.0625,
  7531517.1875,
  7532157.8125,
  7533625.0,
  7550404.6875,
  7550453.125,
  7552068.75,
  7553210.9375,
  7553343.75,
  7553731.25,
  7591914.0625,
  7592870.3125,
  7599943.75,
  7600496.875,
  7612810.9375,
  7623751.5625,
  7626632.8125,
  7626650.0,
  7628281.25,
  7639168.75,
  7641737.5,
  7642992.1875,
  7643592.1875,
  7644710.9375,
  7661100.0,
  7668985.9375,
  7669834.375,
  7683825.0,
  7694121.875,
  7704985.9375,
  7706587.5,
  7712809.375,
  7712868.75,
  7718031.25,
  7718665.625,
  7723415.625,
  7729760.9375,
  7739356.25,
  7745181.25,
  7746237.5,
  7748459.375,
  7752043.75,
  7760135.9375,
  7766709.375,
  7784137.5,
  7813943.75,
  7817521.875,
  7822085.9375,
  7855440.625,
  7857710.9375,
  7860057.8125,
  7860148.4375,
  7861342.1875,
  7861828.125,
  7861918.75,
  7863017.1875,
  7863412.5,
  7863564.0625,
  7863617.1875,
  7864025.0,
  7864389.0625,
  7868554.6875,
  7871267.1875,
  7876717.1875,
  7877762.5,
  7878728.125,
  7880871.875,
  7889887.5,
  7917645.3125,
  7934773.4375,
  7934795.3125,
  7943190.625,
  7967570.3125,
  7969285.9375,
  7971932.8125,
  7974818.75,
  7987201.5625,
  7989114.0625,
  7990282.8125,
  7991693.75,
  7993553.125,
  7995120.3125,
  7995895.3125,
  8009334.375,
  8013240.625,
  8020110.9375,
  8020296.875,
  8021992.1875,
  8040154.6875,
  8042485.9375,
  8044375.0,
  8045057.8125,
  8045390.625,
  8046200.0,
  8046748.4375,
  8046976.5625,
  8047296.875,
  8047625.0,
  8047684.375,
  8047728.125,
  8047748.4375,
  8048518.75,
  8048593.75,
  8048803.125,
  8048856.25,
  8049306.25,
  8049317.1875,
  8049546.875,
  8050032.8125,
  8050335.9375,
  8050664.0625,
  8050670.3125,
  8051039.0625,
  8051062.5,
  8051306.25,
  8051645.3125,
  8051739.0625,
  8052314.0625,
  8054231.25,
  8056465.625,
  8057725.0,
  8058593.75,
  8059387.5,
  8073864.0625,
  8079104.6875,
  8080493.75,
  8083139.0625,
  8083304.6875,
  8087073.4375,
  8087195.3125,
  8087335.9375,
  8087360.9375,
  8087509.375,
  8088851.5625,
  8091026.5625,
  8092528.125,
  8093196.875,
  8093304.6875,
  8094446.875,
  8118250.0,
  8118320.3125,
  8119912.5,
  8121423.4375,
  8121706.25,
  8121717.1875,
  8126790.625,
  8147593.75,
  8150389.0625,
  8153126.5625,
  8159525.0,
  8161943.75,
  8162251.5625,
  8164865.625,
  8167559.375,
  8168656.25,
  8176117.1875,
  8180609.375,
  8180653.125,
  8192901.5625,
  8203507.8125,
  8217309.375,
  8229478.125,
  8229889.0625,
  8230914.0625,
  8234614.0625,
  8235764.0625,
  8237190.625,
  8237760.9375,
  8245637.5,
  8266554.6875,
  8269339.0625,
  8272037.5,
  8275417.1875,
  8276156.25,
  8277659.375,
  8279089.0625,
  8280885.9375,
  8299728.125,
  8302331.25,
  8302506.25,
  8304726.5625,
  8305854.6875,
  8305926.5625,
  8308570.3125,
  8308771.875,
  8315718.75,
  8315821.875,
  8316964.0625,
  8338843.75,
  8345293.75,
  8347054.6875,
  8351057.8125,
  8377525.0,
  8380632.8125,
  8384343.75,
  8384429.6875,
  8385178.125,
  8388120.3125,
  8389101.5625,
  8393829.6875,
  8394310.9375,
  8395693.75,
  8399351.5625,
  8401098.4375,
  8409232.8125,
  8410437.5,
  8411906.25,
  8412567.1875,
  8412601.5625,
  8412737.5,
  8412912.5,
  8413031.25,
  8415842.1875,
  8417812.5,
  8418051.5625,
  8418064.0625,
  8418073.4375,
  8418089.0625,
  8418221.875,
  8422181.25,
  8428195.3125,
  8433054.6875,
  8434867.1875,
  8439981.25,
  8451942.1875,
  8455096.875,
  8455562.5,
  8462157.8125,
  8464734.375,
  8465123.4375,
  8468764.0625,
  8468856.25,
  8469284.375,
  8471246.875,
  8471862.5,
  8474565.625,
  8483171.875,
  8485726.5625,
  8485739.0625,
  8491293.75,
  8491321.875,
  8501990.625,
  8504814.0625,
  8506657.8125,
  8507259.375,
  8507540.625,
  8508042.1875,
  8508593.75,
  8509446.875,
  8510001.5625,
  8510521.875,
  8510620.3125,
  8510718.75,
  8510846.875,
  8513764.0625,
  8514793.75,
  8514920.3125,
  8515700.0,
  8516482.8125,
  8517231.25,
  8519948.4375,
  8529168.75,
  8529276.5625,
  8537926.5625,
  8541981.25,
  8548757.8125,
  8622457.8125,
  8636645.3125,
  8640654.6875,
  8643534.375,
  8644235.9375,
  8645557.8125,
  8645826.5625,
  8649232.8125,
  8650785.9375,
  8651528.125,
  8662882.8125,
  8662923.4375,
  8667153.125,
  8669218.75,
  8669221.875,
  8670859.375,
  8674385.9375,
  8679245.3125,
  8680443.75,
  8709396.875,
  8727407.8125,
  8737837.5,
  8739473.4375,
  8740848.4375,
  8747032.8125,
  8792084.375,
  8838993.75,
  8861401.5625,
  8864907.8125,
  8875876.5625,
  8876231.25,
  8877085.9375,
  8898010.9375,
  8907914.0625,
  8912482.8125,
  8912593.75,
  8912981.25,
  8913306.25,
  8913618.75,
  8914153.125,
  8914568.75,
  8914834.375,
  8914871.875,
  8914982.8125,
  8915453.125,
  8916064.0625,
  8916609.375,
  8918048.4375,
  8918089.0625,
  8920103.125,
  8928975.0,
  8945767.1875,
  8954939.0625,
  8956365.625,
  8958364.0625,
  8958376.5625,
  8960845.3125,
  8962526.5625,
  9010131.25,
  9010667.1875,
  9011198.4375,
  9011273.4375,
  9012182.8125,
  9014937.5,
  9014951.5625,
  9015003.125,
  9018128.125,
  9023029.6875,
  9044710.9375,
  9045473.4375,
  9045540.625,
  9048034.375,
  9094550.0,
  9103175.0,
  9103375.0,
  9107237.5,
  9114895.3125,
  9115864.0625,
  9115937.5,
  9116026.5625,
  9116031.25,
  9117192.1875,
  9117276.5625,
  9117504.6875,
  9118426.5625,
  9119031.25,
  9119481.25,
  9119520.3125,
  9120359.375,
  9121345.3125,
  9125962.5,
  9139135.9375,
  9139292.1875,
  9139881.25,
  9140506.25,
  9140873.4375,
  9144987.5,
  9145659.375,
  9155585.9375,
  9198287.5,
  9228348.4375,
  9249368.75,
  9250046.875,
  9251948.4375,
  9252273.4375,
  9257131.25,
  9268209.375,
  9280120.3125,
  9282406.25,
  9288718.75,
  9297796.875,
  9349317.1875,
  9370676.5625,
  9391842.1875,
  9402478.125,
  9402481.25,
  9415321.875,
  9421348.4375,
  9421409.375,
  9421640.625,
  9421818.75,
  9422517.1875,
  9422782.8125,
  9422992.1875,
  9454517.1875,
  9457410.9375,
  9457925.0,
  9458476.5625,
  9458532.8125,
  9458703.125,
  9459209.375,
  9550031.25,
  9558137.5,
  9562696.875,
  9562915.625,
  9563056.25,
  9563235.9375,
  9563315.625,
  9563631.25,
  9563784.375,
  9565028.125,
  9565415.625,
  9596843.75,
  9599223.4375,
  9600146.875,
  9601960.9375,
  9605525.0,
  9606675.0,
  9634103.125,
  9634468.75,
  9634960.9375,
  9635785.9375,
  9635796.875,
  9636626.5625,
  9636768.75,
  9637006.25,
  9637120.3125,
  9637945.3125,
  9641653.125,
  9643090.625,
  9643400.0,
  9661076.5625,
  9667389.0625,
  9673725.0,
  9697450.0,
  9739523.4375,
  9739548.4375,
  9770415.625,
  9788409.375,
  10077876.5625,
  10325946.875,
  10330818.75,
  10331992.1875,
  10371321.875,
  10398273.4375,
  10410256.25,
  10418064.0625,
  10418176.5625,
  10419162.5,
  10420442.1875,
  10425970.3125,
  10427573.4375,
  10432634.375,
  10447400.0,
  10447464.0625,
  10448446.875,
  10448559.375,
  10450784.375,
  10463826.5625,
  10464318.75,
  10470898.4375,
  10514060.9375,
  10515526.5625,
  10516678.125,
  10516709.375,
  10522045.3125,
  10639671.875,
  10731839.0625,
  10740009.375,
  10740075.0,
  10768415.625,
  10775800.0,
  10795846.875,
  10827765.625,
  10827873.4375,
  10833659.375,
  11007678.125,
  11009654.6875,
  11011651.5625,
  11027709.375,
  11029242.1875,
  11031825.0,
  11037675.0,
  11040745.3125,
  11041187.5,
  11041401.5625,
  11042421.875,
  11042721.875,
  11048273.4375,
  11049246.875,
  11051325.0,
  11051628.125,
  11052604.6875,
  11052657.8125,
  11052859.375,
  11052892.1875,
  11053489.0625,
  11060757.8125,
  11062001.5625,
  11066259.375,
  11067946.875,
  11069446.875,
  11070087.5,
  11070993.75,
  11071000.0,
  11148853.125,
  11155615.625,
  11158576.5625,
  11158876.5625,
  11163339.0625,
  11163839.0625,
  11165573.4375,
  11168075.0,
  11168362.5,
  11168370.3125,
  11170203.125,
  11178767.1875,
  11182475.0,
  11185226.5625,
  11185640.625,
  11191065.625,
  11196364.0625,
  11205457.8125,
  11213829.6875,
  11220510.9375,
  11223576.5625,
  11229248.4375,
  11229265.625,
  11240378.125,
  11257943.75,
  11278460.9375,
  11282954.6875,
  11286684.375,
  11296714.0625,
  11297750.0,
  11299289.0625,
  11300670.3125,
  11301668.75,
  11301767.1875,
  11302646.875,
  11303221.875,
  11308303.125,
  11308631.25,
  11308729.6875,
  11309160.9375,
  11310640.625,
  11311140.625,
  11312531.25,
  11313229.6875,
  11313371.875,
  11317139.0625,
  11317221.875,
  11318100.0,
  11318531.25,
  11318557.8125,
  11318907.8125,
  11318992.1875,
  11319001.5625,
  11319034.375,
  11319064.0625,
  11319132.8125,
  11319290.625,
  11319296.875,
  11319523.4375,
  11320175.0,
  11320187.5,
  11320196.875,
  11320200.0,
  11320207.8125,
  11320307.8125,
  11320379.6875,
  11320445.3125,
  11320484.375,
  11320534.375,
  11320609.375,
  11320632.8125,
  11320685.9375,
  11320806.25,
  11320835.9375,
  11320859.375,
  11320948.4375,
  11321134.375,
  11321195.3125,
  11321317.1875,
  11321325.0,
  11321334.375,
  11321375.0,
  11321445.3125,
  11321451.5625,
  11321548.4375,
  11321631.25,
  11321645.3125,
  11321650.0,
  11321651.5625,
  11321670.3125,
  11321673.4375,
  11321692.1875,
  11321706.25,
  11321707.8125,
  11321787.5,
  11321867.1875,
  11322164.0625,
  11322382.8125,
  11322437.5,
  11322495.3125,
  11324268.75,
  11324734.375,
  11325335.9375,
  11326656.25,
  11326773.4375,
  11327664.0625,
  11328546.875,
  11328968.75,
  11329412.5,
  11330157.8125,
  11330164.0625,
  11335054.6875,
  11336000.0,
  11339350.0,
  11348604.6875],
 [43.75392486847537,
  53.30134253245989,
  20.618502087399357,
  21.924406455376186,
  5.998002358215999,
  14.786688585495128,
  22.093196635757344,
  23.68024596119375,
  49.880642062778556,
  28.545659552558575,
  49.97113678068086,
  5.621115598902164,
  25.40780921453212,
  54.12124961084828,
  115.06108824214039,
  38.02860544665343,
  43.207414702213654,
  128.4824300140363,
  21.715094837566305,
  5.3886399917119405,
  42.922469130228535,
  26.738895845935634,
  25.530083219263126,
  26.733470676399556,
  57.11116651942968,
  17.387852869645496,
  16.529281526471504,
  5.782203723457968,
  35.539905034706365,
  110.89658736604568,
  5.2227790612697245,
  61.0837397063755,
  101.81588581564756,
  5.859871244315133,
  53.88475944792805,
  40.75024709223862,
  52.34716024748916,
  28.739353136881753,
  23.01820822849388,
  25.00720641955666,
  8.78088466984392,
  7.6218485553600335,
  7.174991517234949,
  16.113131589403842,
  51.131353566608624,
  61.233595888689834,
  25.45036325466539,
  8.072804957241642,
  79.32461414853991,
  8.49243355320877,
  12.96251965794333,
  7.078361901396194,
  23.605325737183904,
  46.74114375529501,
  77.12116530033984,
  29.40745167985136,
  158.35104411458886,
  8.418614803314128,
  26.015051982204774,
  55.56095569895757,
  67.79174536866789,
  7.282727483688339,
  5.306245525254198,
  39.616274915319664,
  38.365673509787314,
  15.168393593885424,
  103.00719759890382,
  30.59400773529808,
  60.46513376590568,
  13.59523915647309,
  13.98442482326472,
  18.83340320417493,
  92.45700577003812,
  164.7597848344842,
  6.2216284263079205,
  21.862049316945615,
  47.81199286191655,
  26.235870242488485,
  21.918383834820453,
  109.50129942928082,
  56.85705873108824,
  7.266005216878156,
  74.53906579978519,
  165.2228845648882,
  170.1498074003878,
  119.40422156835832,
  31.404124426836177,
  7.804210999045414,
  280.113029069459,
  22.614670440139495,
  67.94283851166814,
  64.37342961638569,
  21.266848516377248,
  31.599589344211594,
  20.59337399075316,
  12.87531967679679,
  5.061777665303518,
  16.1301893243049,
  5.216226737409376,
  27.132838893759356,
  11.593850624895591,
  81.83205461453588,
  9.913849821767581,
  41.86246897757699,
  73.87737571503175,
  19.69082384754751,
  127.76211110181643,
  5.8189128719845415,
  37.75297030915435,
  5.140248267368501,
  57.13899022260139,
  42.2129512169315,
  21.936811292404215,
  6.995837870222363,
  8.3438645572348,
  21.25739918940065,
  15.872419748789111,
  40.23707435461132,
  5.574288905078166,
  5.945044658998763,
  14.085976750598574,
  67.19019410745126,
  17.97092346727778,
  10.991760657873488,
  24.460192122423337,
  12.84821080003411,
  32.016047783948856,
  7.850894185283803,
  5.160396058089147,
  31.941227597115116,
  95.43425153435207,
  42.631870892468044,
  32.91635765563572,
  21.66045169580405,
  10.508535120488729,
  5.28306949989141,
  54.145080976936484,
  47.33039298365725,
  11.86945278746586,
  16.54348374789551,
  125.31657824084525,
  22.186633516588397,
  5.091154530554696,
  64.75589985183423,
  22.440970856910052,
  45.305497934321124,
  90.15879114157494,
  25.134685854975107,
  5.1347737761070285,
  13.687082883093103,
  96.36901726754965,
  20.16411602819774,
  56.94311118831999,
  11.492030932702049,
  114.24098958373553,
  33.54265949693855,
  9.931351135515833,
  5.887329859625653,
  36.7534253499761,
  15.735713464661806,
  40.83553411376458,
  24.774056696480542,
  53.564450313534635,
  107.53647075017517,
  29.949173515979236,
  5.2062649485443755,
  66.82711371818307,
  34.971916601410264,
  25.569493017896214,
  42.980025219009576,
  30.90573732860115,
  20.396730799438295,
  13.547814446505033,
  5.915621247658618,
  36.08713654857849,
  13.06910721637352,
  7.505004826910014,
  244.05993249000988,
  49.8579705444769,
  45.25762938516009,
  65.40244533665721,
  32.9014808539647,
  84.9042385848715,
  11.661303724722579,
  14.590991437742778,
  84.3619086805655,
  22.92885562569729,
  97.98584093477328,
  74.4383668533196,
  76.68229294493757,
  11.539071769791505,
  8.645717856359637,
  19.895290604577436,
  23.723269455149136,
  27.112570679417374,
  11.879210610849176,
  11.537239648147086,
  36.91157897781484,
  25.16423341634477,
  35.88155545615587,
  25.513483900735274,
  12.748184764742398,
  13.394873733489309,
  69.29429794685359,
  26.28281426419346,
  187.00466972933881,
  7.9704995115765795,
  173.67092396107054,
  76.10002814779887,
  91.88516735600359,
  22.223260547878148,
  111.47409997418461,
  24.043690545754973,
  16.928984469290558,
  33.31780430922666,
  22.369811688294686,
  120.62874672209836,
  54.21027400403291,
  16.591875927195197,
  40.565105739391385,
  36.855526511018034,
  7.333567135139411,
  54.87351521662981,
  5.266270017674993,
  8.432971762130554,
  36.96469894214835,
  104.52967296577759,
  107.82788024407577,
  38.73531221892093,
  36.51860141397936,
  12.18437043220287,
  28.839289086033226,
  15.758084107941405,
  20.700839828265995,
  62.47732521400617,
  48.969913191618375,
  75.01117140744776,
  66.23846180165555,
  85.55648562731454,
  65.1686255242071,
  46.754776157216945,
  6.11904876995819,
  9.799004156789668,
  7.179359480964491,
  31.726751026044646,
  11.892663500080214,
  8.602995999892196,
  25.884921803466998,
  212.8452123476046,
  50.89203883811902,
  5.253993486741754,
  6.552720861958359,
  5.056294872556809,
  86.27666082539653,
  92.0682438066245,
  5.6980376374954504,
  20.141043453919497,
  18.985395531791937,
  70.13011038949608,
  46.851635287902795,
  29.890694816886814,
  11.06396829750044,
  21.244447579686643,
  18.794297835226832,
  81.8384197889961,
  17.11567347736093,
  7.4347814435829465,
  17.11543832697124,
  25.765596526440994,
  171.0189770425673,
  166.62417570530727,
  47.617306423538224,
  87.52555138620906,
  5.099150847659307,
  6.160764826158617,
  33.207044223580425,
  6.225943962290561,
  82.44689715444312,
  130.23094394657033,
  110.29107486064109,
  53.88273073600563,
  44.84042821853197,
  29.59419248412873,
  34.10380275093801,
  86.78925581246611,
  17.36267224747698,
  31.718896982946124,
  12.230208666697209,
  17.204902799252732,
  37.53597700431143,
  92.54196077707205,
  57.055982939784954,
  15.84891231011288,
  19.56022588511453,
  144.26705097138995,
  22.964994816741147,
  29.16119566437618,
  20.622435762921583,
  17.811022479678297,
  97.6487166632215,
  6.511078623653669,
  49.47645177092092,
  28.57188099561196,
  72.048764994024,
  144.93386168548756,
  6.323083148656236,
  15.111041166099659,
  42.031232681637526,
  12.223654060772887,
  5.526743297279815,
  53.08652689302342,
  8.433700076865025,
  110.42941834796935,
  26.322587517522503,
  24.174092846941424,
  78.4021757391126,
  5.74261000245582,
  8.376521969029175,
  17.43699187361873,
  23.101694817855012,
  90.07423512268645,
  48.59528686526562,
  43.772488222419405,
  43.365409914294055,
  21.38342085023234,
  46.31468531162235,
  46.18205156078226,
  84.8699528777467,
  46.05404738548869,
  43.229078007261236,
  30.87182995664722,
  8.436366247547305,
  80.71462110275678,
  10.700143446109799,
  44.734507441241824,
  12.692501109911207,
  6.260513558007149,
  10.122689190026074,
  10.618749267772554,
  40.06237374100719,
  45.14372333344614,
  9.479292995380307,
  88.99744685885068,
  13.781359629670543,
  98.34673049687999,
  60.388927354927034,
  50.22720107871789,
  53.5090445784896,
  83.06641751766367,
  62.60647755524819,
  85.55358745679399,
  30.56143379005863,
  66.2142546385347,
  33.292564619169084,
  5.330715179626472,
  102.88503052828759,
  13.001123588736194,
  23.680216149907416,
  22.350929885657884,
  6.286787999733525,
  52.28705619053464,
  27.914073367470174,
  128.19404468544138,
  48.235538759442974,
  6.730400570813883,
  99.04659770009894,
  18.566916446525155,
  21.590357604646172,
  8.364921840770798,
  5.97513148222563,
  99.36807730698254,
  26.94518336903817,
  31.953018426734953,
  16.47842063191503,
  62.33788944563073,
  53.75135712823922,
  9.585770462043971,
  111.38256205787181,
  5.0856831092483405,
  17.669869388756783,
  18.632841150474132,
  16.592382657074097,
  78.0692523966197,
  5.516222382514227,
  501.37446045827375,
  17.87286395669639,
  33.24100517722623,
  22.572729366525067,
  118.9539994110268,
  17.383326089971284,
  6.177673255921603,
  30.26536075125395,
  34.89182847919665,
  12.935428667451065,
  68.13655334229341,
  69.14336503373974,
  68.68534139384671,
  29.47256093982479,
  63.69146312283162,
  5.9153370570827,
  48.29949717437159,
  66.67555906606268,
  15.204573859599215,
  33.052763747043606,
  52.014793102059656,
  32.75085541957406,
  33.68430278458382,
  17.164569270988352,
  14.505985158509112,
  70.94980489885096,
  155.25802732296813,
  66.06698038312743,
  197.4950889017045,
  52.85158449322144,
  177.9833347556298,
  54.42913933269631,
  37.66726863017692,
  31.321499521215586,
  39.856671284012265,
  20.896251788192018,
  12.741620673160543,
  78.81400404884599,
  20.495216122662654,
  22.719183493831196,
  21.484453049263323,
  55.08032137088479,
  95.71536992761766,
  14.675238278772825,
  27.04244470207559,
  5.677331531976479,
  10.675411814870174,
  43.485536273772304,
  140.01531243048245,
  6.547540752207912,
  6.709089619481133,
  22.07409890641973,
  5.152222989628288,
  25.070068647948318,
  36.17991549075103,
  54.22198599884703,
  75.15239261103841,
  47.02716273158394,
  17.91447372116932,
  31.28091963514686,
  43.05488593947,
  72.03202750997454,
  92.81886229976021,
  156.09729156777183,
  25.655530301825515,
  137.32711929510913,
  9.952254347881473,
  27.025776806655237,
  70.6236706877964,
  87.92116177428154,
  37.18081653779992,
  26.71413628601404,
  9.983651149520439,
  24.4262667547115,
  11.210259423102478,
  42.257604953995205,
  76.56037085296416,
  7.0024741209140124,
  18.612094973680065,
  19.546798175785707,
  29.715450231436222,
  8.692935143158158,
  73.99708871565566,
  21.859953673229406,
  151.0012550239196,
  15.306668697106522,
  26.097428890835417,
  98.34282794889336,
  5.132050572252867,
  11.239138031428249,
  20.125041169697816,
  380.74803274445173,
  33.401300993992955,
  35.594815898855565,
  5.6208934925531455,
  50.378191367867316,
  8.217028233409224,
  10.358428794141663,
  58.6030095307516,
  56.281365905722204,
  38.33917593442391,
  28.269961222263774,
  20.39754991860405,
  9.484820167183592,
  16.070756741727923,
  34.98034773276001,
  65.88387423109285,
  171.43762260308569,
  43.700894668918345,
  31.577902806041813,
  99.68172293680057,
  9.418639423545025,
  25.109288934753867,
  19.33079836537231,
  77.03104227277372,
  23.847740967838693,
  65.22663213820249,
  40.28309944929495,
  153.0397994298229,
  17.740475801007346,
  13.361176381119877,
  76.22616542834456,
  239.31482866215651,
  5.07025553405524,
  8.619375273750922,
  46.861527377546466,
  57.48950066574868,
  44.80035046353777,
  5.468272250057785,
  6.022209334269816,
  9.742251182113371,
  23.346196272156877,
  68.11083438952073,
  93.66070164549868,
  18.887996635082647,
  17.65389205380559,
  7.4263953058941405,
  53.245769431770015,
  127.62136587696972,
  28.086955784993386,
  62.91661044120437,
  32.87485273663541,
  12.093148533363456,
  11.724106924459827,
  17.21892141179511,
  46.35176743779546,
  49.13404485571881,
  5.181754878040886,
  21.038563658034015,
  18.146352252419614,
  27.96253172721343,
  32.4923063076619,
  46.4530751552082,
  20.422302890360182,
  16.966380366088654,
  38.345068188325044,
  27.789861771085373,
  70.16112336982464,
  23.704345732819228,
  53.491457619753724,
  106.05260753675385,
  373.34633863600675,
  14.503795962194422,
  65.94213909246844,
  11.268594675892013,
  21.130121465968895,
  12.50130017743082,
  19.36663731199885,
  59.41211713093734,
  20.98000624233178,
  55.20440709162843,
  20.2643863777146,
  15.56041246628888,
  10.441582803986975,
  69.56215342613369,
  7.0968916243289,
  23.554188031337553,
  49.534255126767874,
  50.69777400029342,
  44.94298789934566,
  28.878549813415905,
  103.96705180896483,
  22.65670289940123,
  77.97790652855814,
  28.599325195318094,
  16.628217858213844,
  5.6061062657207845,
  12.919991332486576,
  10.495335778690041,
  17.053526752436763,
  16.342034641945663,
  6.9474412753616575,
  19.78051052151996,
  20.698785774917553,
  42.92041201862068,
  176.96140839714167,
  7.533339199150585,
  62.956313733053065,
  74.45810518930293,
  20.747312610269052,
  15.651905448000619,
  115.27646964125759,
  71.35468514857598,
  58.42951223027682,
  49.98773070886058,
  50.88248829304828,
  5.992713265028529,
  80.22818625145914,
  57.20046222754387,
  21.803440267492345,
  16.132821937515907,
  72.54960544831206,
  7.625923488230528,
  30.638369917801967,
  10.857609179556457,
  15.918213539817387,
  20.431950772799226,
  6.378697211437925,
  114.58964839656181,
  7.544904728209705,
  27.297615513093156,
  54.638669822191154,
  23.62640971925502,
  37.93452452601168,
  41.41350395965277,
  49.979079118055594,
  14.770197475017461,
  19.47726597568899,
  48.82832267206836,
  68.87290089420236,
  62.40582323431573,
  30.597380181581805,
  18.503564931578502,
  64.67488627528273,
  40.31119269849788,
  26.020123432969818,
  18.897969795126812,
  28.49062108291132,
  76.82417461303996,
  49.864464652641814,
  43.86366796467256,
  38.08146936614241,
  5.273975770415842,
  27.998411736108974,
  12.289328547466589,
  33.91695098318348,
  16.864499712731543,
  12.465866299430061,
  37.470479294181175,
  54.15075259243452,
  5.078880195492166,
  103.72940872669705,
  6.645027463501742,
  7.420427532244959,
  17.407716903590135,
  10.861476476316085,
  8.52815305330454,
  5.783638738817864,
  26.57606061557894,
  95.45959219224747,
  80.39329081427108,
  122.72202930155345,
  18.150249471841473,
  5.716106541369781,
  17.43718729812661,
  41.30287211859103,
  75.30350851567653,
  8.267309959365393,
  10.312053380372127,
  5.144619479335722,
  65.9520710769438,
  9.02487398823923,
  57.56725091644046,
  103.4028634400867])

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

([3587735.9375,
  6802431.25,
  6842682.8125,
  6920104.6875,
  6949962.5,
  6950360.9375,
  7015684.375,
  7017779.6875,
  7106253.125,
  7188793.75,
  7195482.8125,
  7195539.0625,
  7196743.75,
  7197140.625,
  7275157.8125,
  7382965.625,
  7386604.6875,
  7390162.5,
  7391800.0,
  7391832.8125,
  7415934.375,
  7424310.9375,
  7468432.8125,
  7468512.5,
  7469973.4375,
  7470921.875,
  7473642.1875,
  7473678.125,
  7475582.8125,
  7479473.4375,
  7498467.1875,
  7503434.375,
  7506903.125,
  7508653.125,
  7508834.375,
  7513351.5625,
  7513446.875,
  7514012.5,
  7516193.75,
  7516782.8125,
  7517039.0625,
  7517884.375,
  7526164.0625,
  7531517.1875,
  7532157.8125,
  7533625.0,
  7550404.6875,
  7550453.125,
  7552068.75,
  7553210.9375,
  7553343.75,
  7553731.25,
  7591914.0625,
  7592870.3125,
  7599943.75,
  7600496.875,
  7612810.9375,
  7623751.5625,
  7626632.8125,
  7626650.0,
  7628281.25,
  7639168.75,
  7641737.5,
  7642992.1875,
  7643592.1875,
  7644710.9375,
  7661100.0,
  7668985.9375,
  7669834.375,
  7683825.0,
  7694121.875,
  7704985.9375,
  7706587.5,
  7712809.375,
  7712868.75,
  7718031.25,
  7718665.625,
  7723415.625,
  7729760.9375,
  7739356.25,
  7745181.25,
  7746237.5,
  7748459.375,
  7752043.75,
  7760135.9375,
  7766709.375,
  7784137.5,
  7813943.75,
  7817521.875,
  7822085.9375,
  7855440.625,
  7857710.9375,
  7860057.8125,
  7860148.4375,
  7861342.1875,
  7861828.125,
  7861918.75,
  7863017.1875,
  7863412.5,
  7863564.0625,
  7863617.1875,
  7864025.0,
  7864389.0625,
  7868554.6875,
  7871267.1875,
  7876717.1875,
  7877762.5,
  7878728.125,
  7880871.875,
  7889887.5,
  7917645.3125,
  7934773.4375,
  7934795.3125,
  7943190.625,
  7967570.3125,
  7969285.9375,
  7971932.8125,
  7974818.75,
  7987201.5625,
  7989114.0625,
  7990282.8125,
  7991693.75,
  7993553.125,
  7995120.3125,
  7995895.3125,
  8009334.375,
  8013240.625,
  8020110.9375,
  8020296.875,
  8021992.1875,
  8040154.6875,
  8042485.9375,
  8044375.0,
  8045057.8125,
  8045390.625,
  8046200.0,
  8046748.4375,
  8046976.5625,
  8047296.875,
  8047625.0,
  8047684.375,
  8047728.125,
  8047748.4375,
  8048518.75,
  8048593.75,
  8048803.125,
  8048856.25,
  8049306.25,
  8049317.1875,
  8049546.875,
  8050032.8125,
  8050335.9375,
  8050664.0625,
  8050670.3125,
  8051039.0625,
  8051062.5,
  8051306.25,
  8051645.3125,
  8051739.0625,
  8052314.0625,
  8054231.25,
  8056465.625,
  8057725.0,
  8058593.75,
  8059387.5,
  8073864.0625,
  8079104.6875,
  8080493.75,
  8083139.0625,
  8083304.6875,
  8087073.4375,
  8087195.3125,
  8087335.9375,
  8087360.9375,
  8087509.375,
  8088851.5625,
  8091026.5625,
  8092528.125,
  8093196.875,
  8093304.6875,
  8094446.875,
  8118250.0,
  8118320.3125,
  8119912.5,
  8121423.4375,
  8121706.25,
  8121717.1875,
  8126790.625,
  8147593.75,
  8150389.0625,
  8153126.5625,
  8159525.0,
  8161943.75,
  8162251.5625,
  8164865.625,
  8167559.375,
  8168656.25,
  8176117.1875,
  8180609.375,
  8180653.125,
  8192901.5625,
  8203507.8125,
  8217309.375,
  8229478.125,
  8229889.0625,
  8230914.0625,
  8234614.0625,
  8235764.0625,
  8237190.625,
  8237760.9375,
  8245637.5,
  8266554.6875,
  8269339.0625,
  8272037.5,
  8275417.1875,
  8276156.25,
  8277659.375,
  8279089.0625,
  8280885.9375,
  8299728.125,
  8302331.25,
  8302506.25,
  8304726.5625,
  8305854.6875,
  8305926.5625,
  8308570.3125,
  8308771.875,
  8315718.75,
  8315821.875,
  8316964.0625,
  8338843.75,
  8345293.75,
  8347054.6875,
  8351057.8125,
  8377525.0,
  8380632.8125,
  8384343.75,
  8384429.6875,
  8385178.125,
  8388120.3125,
  8389101.5625,
  8393829.6875,
  8394310.9375,
  8395693.75,
  8399351.5625,
  8401098.4375,
  8409232.8125,
  8410437.5,
  8411906.25,
  8412567.1875,
  8412601.5625,
  8412737.5,
  8412912.5,
  8413031.25,
  8415842.1875,
  8417812.5,
  8418051.5625,
  8418064.0625,
  8418073.4375,
  8418089.0625,
  8418221.875,
  8422181.25,
  8428195.3125,
  8433054.6875,
  8434867.1875,
  8439981.25,
  8451942.1875,
  8455096.875,
  8455562.5,
  8462157.8125,
  8464734.375,
  8465123.4375,
  8468764.0625,
  8468856.25,
  8469284.375,
  8471246.875,
  8471862.5,
  8474565.625,
  8483171.875,
  8485726.5625,
  8485739.0625,
  8491293.75,
  8491321.875,
  8501990.625,
  8504814.0625,
  8506657.8125,
  8507259.375,
  8507540.625,
  8508042.1875,
  8508593.75,
  8509446.875,
  8510001.5625,
  8510521.875,
  8510620.3125,
  8510718.75,
  8510846.875,
  8513764.0625,
  8514793.75,
  8514920.3125,
  8515700.0,
  8516482.8125,
  8517231.25,
  8519948.4375,
  8529168.75,
  8529276.5625,
  8537926.5625,
  8541981.25,
  8548757.8125,
  8622457.8125,
  8636645.3125,
  8640654.6875,
  8643534.375,
  8644235.9375,
  8645557.8125,
  8645826.5625,
  8649232.8125,
  8650785.9375,
  8651528.125,
  8662882.8125,
  8662923.4375,
  8667153.125,
  8669218.75,
  8669221.875,
  8670859.375,
  8674385.9375,
  8679245.3125,
  8680443.75,
  8709396.875,
  8727407.8125,
  8737837.5,
  8739473.4375,
  8740848.4375,
  8747032.8125,
  8792084.375,
  8838993.75,
  8861401.5625,
  8864907.8125,
  8875876.5625,
  8876231.25,
  8877085.9375,
  8898010.9375,
  8907914.0625,
  8912482.8125,
  8912593.75,
  8912981.25,
  8913306.25,
  8913618.75,
  8914153.125,
  8914568.75,
  8914834.375,
  8914871.875,
  8914982.8125,
  8915453.125,
  8916064.0625,
  8916609.375,
  8918048.4375,
  8918089.0625,
  8920103.125,
  8928975.0,
  8945767.1875,
  8954939.0625,
  8956365.625,
  8958364.0625,
  8958376.5625,
  8960845.3125,
  8962526.5625,
  9010131.25,
  9010667.1875,
  9011198.4375,
  9011273.4375,
  9012182.8125,
  9014937.5,
  9014951.5625,
  9015003.125,
  9018128.125,
  9023029.6875,
  9044710.9375,
  9045473.4375,
  9045540.625,
  9048034.375,
  9094550.0,
  9103175.0,
  9103375.0,
  9107237.5,
  9114895.3125,
  9115864.0625,
  9115937.5,
  9116026.5625,
  9116031.25,
  9117192.1875,
  9117276.5625,
  9117504.6875,
  9118426.5625,
  9119031.25,
  9119481.25,
  9119520.3125,
  9120359.375,
  9121345.3125,
  9125962.5,
  9139135.9375,
  9139292.1875,
  9139881.25,
  9140506.25,
  9140873.4375,
  9144987.5,
  9145659.375,
  9155585.9375,
  9198287.5,
  9228348.4375,
  9249368.75,
  9250046.875,
  9251948.4375,
  9252273.4375,
  9257131.25,
  9268209.375,
  9280120.3125,
  9282406.25,
  9288718.75,
  9297796.875,
  9349317.1875,
  9370676.5625,
  9391842.1875,
  9402478.125,
  9402481.25,
  9415321.875,
  9421348.4375,
  9421409.375,
  9421640.625,
  9421818.75,
  9422517.1875,
  9422782.8125,
  9422992.1875,
  9454517.1875,
  9457410.9375,
  9457925.0,
  9458476.5625,
  9458532.8125,
  9458703.125,
  9459209.375,
  9550031.25,
  9558137.5,
  9562696.875,
  9562915.625,
  9563056.25,
  9563235.9375,
  9563315.625,
  9563631.25,
  9563784.375,
  9565028.125,
  9565415.625,
  9596843.75,
  9599223.4375,
  9600146.875,
  9601960.9375,
  9605525.0,
  9606675.0,
  9634103.125,
  9634468.75,
  9634960.9375,
  9635785.9375,
  9635796.875,
  9636626.5625,
  9636768.75,
  9637006.25,
  9637120.3125,
  9637945.3125,
  9641653.125,
  9643090.625,
  9643400.0,
  9661076.5625,
  9667389.0625,
  9673725.0,
  9697450.0,
  9739523.4375,
  9739548.4375,
  9770415.625,
  9788409.375,
  10077876.5625,
  10325946.875,
  10330818.75,
  10331992.1875,
  10371321.875,
  10398273.4375,
  10410256.25,
  10418064.0625,
  10418176.5625,
  10419162.5,
  10420442.1875,
  10425970.3125,
  10427573.4375,
  10432634.375,
  10447400.0,
  10447464.0625,
  10448446.875,
  10448559.375,
  10450784.375,
  10463826.5625,
  10464318.75,
  10470898.4375,
  10514060.9375,
  10515526.5625,
  10516678.125,
  10516709.375,
  10522045.3125,
  10639671.875,
  10731839.0625,
  10740009.375,
  10740075.0,
  10768415.625,
  10775800.0,
  10795846.875,
  10827765.625,
  10827873.4375,
  10833659.375,
  11007678.125,
  11009654.6875,
  11011651.5625,
  11027709.375,
  11029242.1875,
  11031825.0,
  11037675.0,
  11040745.3125,
  11041187.5,
  11041401.5625,
  11042421.875,
  11042721.875,
  11048273.4375,
  11049246.875,
  11051325.0,
  11051628.125,
  11052604.6875,
  11052657.8125,
  11052859.375,
  11052892.1875,
  11053489.0625,
  11060757.8125,
  11062001.5625,
  11066259.375,
  11067946.875,
  11069446.875,
  11070087.5,
  11070993.75,
  11071000.0,
  11148853.125,
  11155615.625,
  11158576.5625,
  11158876.5625,
  11163339.0625,
  11163839.0625,
  11165573.4375,
  11168075.0,
  11168362.5,
  11168370.3125,
  11170203.125,
  11178767.1875,
  11182475.0,
  11185226.5625,
  11185640.625,
  11191065.625,
  11196364.0625,
  11205457.8125,
  11213829.6875,
  11220510.9375,
  11223576.5625,
  11229248.4375,
  11229265.625,
  11240378.125,
  11257943.75,
  11278460.9375,
  11282954.6875,
  11286684.375,
  11296714.0625,
  11297750.0,
  11299289.0625,
  11300670.3125,
  11301668.75,
  11301767.1875,
  11302646.875,
  11303221.875,
  11308303.125,
  11308631.25,
  11308729.6875,
  11309160.9375,
  11310640.625,
  11311140.625,
  11312531.25,
  11313229.6875,
  11313371.875,
  11317139.0625,
  11317221.875,
  11318100.0,
  11318531.25,
  11318557.8125,
  11318907.8125,
  11318992.1875,
  11319001.5625,
  11319034.375,
  11319064.0625,
  11319132.8125,
  11319290.625,
  11319296.875,
  11319523.4375,
  11320175.0,
  11320187.5,
  11320196.875,
  11320200.0,
  11320207.8125,
  11320307.8125,
  11320379.6875,
  11320445.3125,
  11320484.375,
  11320534.375,
  11320609.375,
  11320632.8125,
  11320685.9375,
  11320806.25,
  11320835.9375,
  11320859.375,
  11320948.4375,
  11321134.375,
  11321195.3125,
  11321317.1875,
  11321325.0,
  11321334.375,
  11321375.0,
  11321445.3125,
  11321451.5625,
  11321548.4375,
  11321631.25,
  11321645.3125,
  11321650.0,
  11321651.5625,
  11321670.3125,
  11321673.4375,
  11321692.1875,
  11321706.25,
  11321707.8125,
  11321787.5,
  11321867.1875,
  11322164.0625,
  11322382.8125,
  11322437.5,
  11322495.3125,
  11324268.75,
  11324734.375,
  11325335.9375,
  11326656.25,
  11326773.4375,
  11327664.0625,
  11328546.875,
  11328968.75,
  11329412.5,
  11330157.8125,
  11330164.0625,
  11335054.6875,
  11336000.0,
  11339350.0,
  11348604.6875],
 [43.75392486847537,
  53.30134253245989,
  20.618502087399357,
  21.924406455376186,
  5.998002358215999,
  14.786688585495128,
  22.093196635757344,
  23.68024596119375,
  49.880642062778556,
  28.545659552558575,
  49.97113678068086,
  5.621115598902164,
  25.40780921453212,
  54.12124961084828,
  115.06108824214039,
  38.02860544665343,
  43.207414702213654,
  128.4824300140363,
  21.715094837566305,
  5.3886399917119405,
  42.922469130228535,
  26.738895845935634,
  25.530083219263126,
  26.733470676399556,
  57.11116651942968,
  17.387852869645496,
  16.529281526471504,
  5.782203723457968,
  35.539905034706365,
  110.89658736604568,
  5.2227790612697245,
  61.0837397063755,
  101.81588581564756,
  5.859871244315133,
  53.88475944792805,
  40.75024709223862,
  52.34716024748916,
  28.739353136881753,
  23.01820822849388,
  25.00720641955666,
  8.78088466984392,
  7.6218485553600335,
  7.174991517234949,
  16.113131589403842,
  51.131353566608624,
  61.233595888689834,
  25.45036325466539,
  8.072804957241642,
  79.32461414853991,
  8.49243355320877,
  12.96251965794333,
  7.078361901396194,
  23.605325737183904,
  46.74114375529501,
  77.12116530033984,
  29.40745167985136,
  158.35104411458886,
  8.418614803314128,
  26.015051982204774,
  55.56095569895757,
  67.79174536866789,
  7.282727483688339,
  5.306245525254198,
  39.616274915319664,
  38.365673509787314,
  15.168393593885424,
  103.00719759890382,
  30.59400773529808,
  60.46513376590568,
  13.59523915647309,
  13.98442482326472,
  18.83340320417493,
  92.45700577003812,
  164.7597848344842,
  6.2216284263079205,
  21.862049316945615,
  47.81199286191655,
  26.235870242488485,
  21.918383834820453,
  109.50129942928082,
  56.85705873108824,
  7.266005216878156,
  74.53906579978519,
  165.2228845648882,
  170.1498074003878,
  119.40422156835832,
  31.404124426836177,
  7.804210999045414,
  280.113029069459,
  22.614670440139495,
  67.94283851166814,
  64.37342961638569,
  21.266848516377248,
  31.599589344211594,
  20.59337399075316,
  12.87531967679679,
  5.061777665303518,
  16.1301893243049,
  5.216226737409376,
  27.132838893759356,
  11.593850624895591,
  81.83205461453588,
  9.913849821767581,
  41.86246897757699,
  73.87737571503175,
  19.69082384754751,
  127.76211110181643,
  5.8189128719845415,
  37.75297030915435,
  5.140248267368501,
  57.13899022260139,
  42.2129512169315,
  21.936811292404215,
  6.995837870222363,
  8.3438645572348,
  21.25739918940065,
  15.872419748789111,
  40.23707435461132,
  5.574288905078166,
  5.945044658998763,
  14.085976750598574,
  67.19019410745126,
  17.97092346727778,
  10.991760657873488,
  24.460192122423337,
  12.84821080003411,
  32.016047783948856,
  7.850894185283803,
  5.160396058089147,
  31.941227597115116,
  95.43425153435207,
  42.631870892468044,
  32.91635765563572,
  21.66045169580405,
  10.508535120488729,
  5.28306949989141,
  54.145080976936484,
  47.33039298365725,
  11.86945278746586,
  16.54348374789551,
  125.31657824084525,
  22.186633516588397,
  5.091154530554696,
  64.75589985183423,
  22.440970856910052,
  45.305497934321124,
  90.15879114157494,
  25.134685854975107,
  5.1347737761070285,
  13.687082883093103,
  96.36901726754965,
  20.16411602819774,
  56.94311118831999,
  11.492030932702049,
  114.24098958373553,
  33.54265949693855,
  9.931351135515833,
  5.887329859625653,
  36.7534253499761,
  15.735713464661806,
  40.83553411376458,
  24.774056696480542,
  53.564450313534635,
  107.53647075017517,
  29.949173515979236,
  5.2062649485443755,
  66.82711371818307,
  34.971916601410264,
  25.569493017896214,
  42.980025219009576,
  30.90573732860115,
  20.396730799438295,
  13.547814446505033,
  5.915621247658618,
  36.08713654857849,
  13.06910721637352,
  7.505004826910014,
  244.05993249000988,
  49.8579705444769,
  45.25762938516009,
  65.40244533665721,
  32.9014808539647,
  84.9042385848715,
  11.661303724722579,
  14.590991437742778,
  84.3619086805655,
  22.92885562569729,
  97.98584093477328,
  74.4383668533196,
  76.68229294493757,
  11.539071769791505,
  8.645717856359637,
  19.895290604577436,
  23.723269455149136,
  27.112570679417374,
  11.879210610849176,
  11.537239648147086,
  36.91157897781484,
  25.16423341634477,
  35.88155545615587,
  25.513483900735274,
  12.748184764742398,
  13.394873733489309,
  69.29429794685359,
  26.28281426419346,
  187.00466972933881,
  7.9704995115765795,
  173.67092396107054,
  76.10002814779887,
  91.88516735600359,
  22.223260547878148,
  111.47409997418461,
  24.043690545754973,
  16.928984469290558,
  33.31780430922666,
  22.369811688294686,
  120.62874672209836,
  54.21027400403291,
  16.591875927195197,
  40.565105739391385,
  36.855526511018034,
  7.333567135139411,
  54.87351521662981,
  5.266270017674993,
  8.432971762130554,
  36.96469894214835,
  104.52967296577759,
  107.82788024407577,
  38.73531221892093,
  36.51860141397936,
  12.18437043220287,
  28.839289086033226,
  15.758084107941405,
  20.700839828265995,
  62.47732521400617,
  48.969913191618375,
  75.01117140744776,
  66.23846180165555,
  85.55648562731454,
  65.1686255242071,
  46.754776157216945,
  6.11904876995819,
  9.799004156789668,
  7.179359480964491,
  31.726751026044646,
  11.892663500080214,
  8.602995999892196,
  25.884921803466998,
  212.8452123476046,
  50.89203883811902,
  5.253993486741754,
  6.552720861958359,
  5.056294872556809,
  86.27666082539653,
  92.0682438066245,
  5.6980376374954504,
  20.141043453919497,
  18.985395531791937,
  70.13011038949608,
  46.851635287902795,
  29.890694816886814,
  11.06396829750044,
  21.244447579686643,
  18.794297835226832,
  81.8384197889961,
  17.11567347736093,
  7.4347814435829465,
  17.11543832697124,
  25.765596526440994,
  171.0189770425673,
  166.62417570530727,
  47.617306423538224,
  87.52555138620906,
  5.099150847659307,
  6.160764826158617,
  33.207044223580425,
  6.225943962290561,
  82.44689715444312,
  130.23094394657033,
  110.29107486064109,
  53.88273073600563,
  44.84042821853197,
  29.59419248412873,
  34.10380275093801,
  86.78925581246611,
  17.36267224747698,
  31.718896982946124,
  12.230208666697209,
  17.204902799252732,
  37.53597700431143,
  92.54196077707205,
  57.055982939784954,
  15.84891231011288,
  19.56022588511453,
  144.26705097138995,
  22.964994816741147,
  29.16119566437618,
  20.622435762921583,
  17.811022479678297,
  97.6487166632215,
  6.511078623653669,
  49.47645177092092,
  28.57188099561196,
  72.048764994024,
  144.93386168548756,
  6.323083148656236,
  15.111041166099659,
  42.031232681637526,
  12.223654060772887,
  5.526743297279815,
  53.08652689302342,
  8.433700076865025,
  110.42941834796935,
  26.322587517522503,
  24.174092846941424,
  78.4021757391126,
  5.74261000245582,
  8.376521969029175,
  17.43699187361873,
  23.101694817855012,
  90.07423512268645,
  48.59528686526562,
  43.772488222419405,
  43.365409914294055,
  21.38342085023234,
  46.31468531162235,
  46.18205156078226,
  84.8699528777467,
  46.05404738548869,
  43.229078007261236,
  30.87182995664722,
  8.436366247547305,
  80.71462110275678,
  10.700143446109799,
  44.734507441241824,
  12.692501109911207,
  6.260513558007149,
  10.122689190026074,
  10.618749267772554,
  40.06237374100719,
  45.14372333344614,
  9.479292995380307,
  88.99744685885068,
  13.781359629670543,
  98.34673049687999,
  60.388927354927034,
  50.22720107871789,
  53.5090445784896,
  83.06641751766367,
  62.60647755524819,
  85.55358745679399,
  30.56143379005863,
  66.2142546385347,
  33.292564619169084,
  5.330715179626472,
  102.88503052828759,
  13.001123588736194,
  23.680216149907416,
  22.350929885657884,
  6.286787999733525,
  52.28705619053464,
  27.914073367470174,
  128.19404468544138,
  48.235538759442974,
  6.730400570813883,
  99.04659770009894,
  18.566916446525155,
  21.590357604646172,
  8.364921840770798,
  5.97513148222563,
  99.36807730698254,
  26.94518336903817,
  31.953018426734953,
  16.47842063191503,
  62.33788944563073,
  53.75135712823922,
  9.585770462043971,
  111.38256205787181,
  5.0856831092483405,
  17.669869388756783,
  18.632841150474132,
  16.592382657074097,
  78.0692523966197,
  5.516222382514227,
  501.37446045827375,
  17.87286395669639,
  33.24100517722623,
  22.572729366525067,
  118.9539994110268,
  17.383326089971284,
  6.177673255921603,
  30.26536075125395,
  34.89182847919665,
  12.935428667451065,
  68.13655334229341,
  69.14336503373974,
  68.68534139384671,
  29.47256093982479,
  63.69146312283162,
  5.9153370570827,
  48.29949717437159,
  66.67555906606268,
  15.204573859599215,
  33.052763747043606,
  52.014793102059656,
  32.75085541957406,
  33.68430278458382,
  17.164569270988352,
  14.505985158509112,
  70.94980489885096,
  155.25802732296813,
  66.06698038312743,
  197.4950889017045,
  52.85158449322144,
  177.9833347556298,
  54.42913933269631,
  37.66726863017692,
  31.321499521215586,
  39.856671284012265,
  20.896251788192018,
  12.741620673160543,
  78.81400404884599,
  20.495216122662654,
  22.719183493831196,
  21.484453049263323,
  55.08032137088479,
  95.71536992761766,
  14.675238278772825,
  27.04244470207559,
  5.677331531976479,
  10.675411814870174,
  43.485536273772304,
  140.01531243048245,
  6.547540752207912,
  6.709089619481133,
  22.07409890641973,
  5.152222989628288,
  25.070068647948318,
  36.17991549075103,
  54.22198599884703,
  75.15239261103841,
  47.02716273158394,
  17.91447372116932,
  31.28091963514686,
  43.05488593947,
  72.03202750997454,
  92.81886229976021,
  156.09729156777183,
  25.655530301825515,
  137.32711929510913,
  9.952254347881473,
  27.025776806655237,
  70.6236706877964,
  87.92116177428154,
  37.18081653779992,
  26.71413628601404,
  9.983651149520439,
  24.4262667547115,
  11.210259423102478,
  42.257604953995205,
  76.56037085296416,
  7.0024741209140124,
  18.612094973680065,
  19.546798175785707,
  29.715450231436222,
  8.692935143158158,
  73.99708871565566,
  21.859953673229406,
  151.0012550239196,
  15.306668697106522,
  26.097428890835417,
  98.34282794889336,
  5.132050572252867,
  11.239138031428249,
  20.125041169697816,
  380.74803274445173,
  33.401300993992955,
  35.594815898855565,
  5.6208934925531455,
  50.378191367867316,
  8.217028233409224,
  10.358428794141663,
  58.6030095307516,
  56.281365905722204,
  38.33917593442391,
  28.269961222263774,
  20.39754991860405,
  9.484820167183592,
  16.070756741727923,
  34.98034773276001,
  65.88387423109285,
  171.43762260308569,
  43.700894668918345,
  31.577902806041813,
  99.68172293680057,
  9.418639423545025,
  25.109288934753867,
  19.33079836537231,
  77.03104227277372,
  23.847740967838693,
  65.22663213820249,
  40.28309944929495,
  153.0397994298229,
  17.740475801007346,
  13.361176381119877,
  76.22616542834456,
  239.31482866215651,
  5.07025553405524,
  8.619375273750922,
  46.861527377546466,
  57.48950066574868,
  44.80035046353777,
  5.468272250057785,
  6.022209334269816,
  9.742251182113371,
  23.346196272156877,
  68.11083438952073,
  93.66070164549868,
  18.887996635082647,
  17.65389205380559,
  7.4263953058941405,
  53.245769431770015,
  127.62136587696972,
  28.086955784993386,
  62.91661044120437,
  32.87485273663541,
  12.093148533363456,
  11.724106924459827,
  17.21892141179511,
  46.35176743779546,
  49.13404485571881,
  5.181754878040886,
  21.038563658034015,
  18.146352252419614,
  27.96253172721343,
  32.4923063076619,
  46.4530751552082,
  20.422302890360182,
  16.966380366088654,
  38.345068188325044,
  27.789861771085373,
  70.16112336982464,
  23.704345732819228,
  53.491457619753724,
  106.05260753675385,
  373.34633863600675,
  14.503795962194422,
  65.94213909246844,
  11.268594675892013,
  21.130121465968895,
  12.50130017743082,
  19.36663731199885,
  59.41211713093734,
  20.98000624233178,
  55.20440709162843,
  20.2643863777146,
  15.56041246628888,
  10.441582803986975,
  69.56215342613369,
  7.0968916243289,
  23.554188031337553,
  49.534255126767874,
  50.69777400029342,
  44.94298789934566,
  28.878549813415905,
  103.96705180896483,
  22.65670289940123,
  77.97790652855814,
  28.599325195318094,
  16.628217858213844,
  5.6061062657207845,
  12.919991332486576,
  10.495335778690041,
  17.053526752436763,
  16.342034641945663,
  6.9474412753616575,
  19.78051052151996,
  20.698785774917553,
  42.92041201862068,
  176.96140839714167,
  7.533339199150585,
  62.956313733053065,
  74.45810518930293,
  20.747312610269052,
  15.651905448000619,
  115.27646964125759,
  71.35468514857598,
  58.42951223027682,
  49.98773070886058,
  50.88248829304828,
  5.992713265028529,
  80.22818625145914,
  57.20046222754387,
  21.803440267492345,
  16.132821937515907,
  72.54960544831206,
  7.625923488230528,
  30.638369917801967,
  10.857609179556457,
  15.918213539817387,
  20.431950772799226,
  6.378697211437925,
  114.58964839656181,
  7.544904728209705,
  27.297615513093156,
  54.638669822191154,
  23.62640971925502,
  37.93452452601168,
  41.41350395965277,
  49.979079118055594,
  14.770197475017461,
  19.47726597568899,
  48.82832267206836,
  68.87290089420236,
  62.40582323431573,
  30.597380181581805,
  18.503564931578502,
  64.67488627528273,
  40.31119269849788,
  26.020123432969818,
  18.897969795126812,
  28.49062108291132,
  76.82417461303996,
  49.864464652641814,
  43.86366796467256,
  38.08146936614241,
  5.273975770415842,
  27.998411736108974,
  12.289328547466589,
  33.91695098318348,
  16.864499712731543,
  12.465866299430061,
  37.470479294181175,
  54.15075259243452,
  5.078880195492166,
  103.72940872669705,
  6.645027463501742,
  7.420427532244959,
  17.407716903590135,
  10.861476476316085,
  8.52815305330454,
  5.783638738817864,
  26.57606061557894,
  95.45959219224747,
  80.39329081427108,
  122.72202930155345,
  18.150249471841473,
  5.716106541369781,
  17.43718729812661,
  41.30287211859103,
  75.30350851567653,
  8.267309959365393,
  10.312053380372127,
  5.144619479335722,
  65.9520710769438,
  9.02487398823923,
  57.56725091644046,
  103.4028634400867])

Graph for home page

Detected energies during the discharge + Energy spectrum

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