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

([6181285.9375,
  6197781.25,
  6239173.4375,
  6247750.0,
  6248804.6875,
  6262335.9375,
  6272785.9375,
  6279440.625,
  6280354.6875,
  6286471.875,
  6286995.3125,
  6293664.0625,
  6295378.125,
  6295903.125,
  6314421.875,
  6353832.8125,
  6392850.0,
  6421156.25,
  6507553.125,
  6508446.875,
  6520773.4375,
  6536207.8125,
  6545845.3125,
  6552589.0625,
  6576779.6875,
  6580757.8125,
  6580773.4375,
  6595562.5,
  6629950.0,
  6631928.125,
  6666404.6875,
  6668471.875,
  6668517.1875,
  6669001.5625,
  6669265.625,
  6673604.6875,
  6675820.3125,
  6680060.9375,
  6691098.4375,
  6718495.3125,
  6723492.1875,
  6731370.3125,
  6733576.5625,
  6759671.875,
  6761328.125,
  6780406.25,
  6786860.9375,
  6822914.0625,
  6825746.875,
  6828890.625,
  6850154.6875,
  6850968.75,
  6850970.3125,
  6855337.5,
  6870803.125,
  6878707.8125,
  6882725.0,
  6895995.3125,
  6929346.875,
  6931592.1875,
  6934289.0625,
  6935462.5,
  6938654.6875,
  6939896.875,
  6960801.5625,
  6967254.6875,
  6993356.25,
  6994426.5625,
  6994682.8125,
  6998009.375,
  7004154.6875,
  7032918.75,
  7063692.1875,
  7097726.5625,
  7100550.0,
  7107506.25,
  7116837.5,
  7130043.75,
  7138718.75,
  7142651.5625,
  7149993.75,
  7169731.25,
  7170189.0625,
  7174621.875,
  7175293.75,
  7175346.875,
  7177920.3125,
  7177967.1875,
  7178021.875,
  7178173.4375,
  7179323.4375,
  7179867.1875,
  7180895.3125,
  7182092.1875,
  7212487.5,
  7215584.375,
  7217818.75,
  7229596.875,
  7229642.1875,
  7231545.3125,
  7242971.875,
  7263040.625,
  7263084.375,
  7268900.0,
  7270918.75,
  7289931.25,
  7291128.125,
  7292212.5,
  7292292.1875,
  7294689.0625,
  7295395.3125,
  7295718.75,
  7298059.375,
  7304248.4375,
  7309814.0625,
  7314185.9375,
  7319234.375,
  7327867.1875,
  7330954.6875,
  7331865.625,
  7333028.125,
  7335623.4375,
  7336623.4375,
  7336845.3125,
  7337079.6875,
  7337087.5,
  7337167.1875,
  7337795.3125,
  7339092.1875,
  7353020.3125,
  7356026.5625,
  7369845.3125,
  7371171.875,
  7373626.5625,
  7373648.4375,
  7377445.3125,
  7416301.5625,
  7433556.25,
  7440907.8125,
  7450426.5625,
  7451001.5625,
  7452140.625,
  7454321.875,
  7456137.5,
  7467356.25,
  7467365.625,
  7489384.375,
  7490357.8125,
  7490390.625,
  7503518.75,
  7509009.375,
  7509067.1875,
  7530709.375,
  7539157.8125,
  7548681.25,
  7581431.25,
  7587604.6875,
  7587640.625,
  7595981.25,
  7644993.75,
  7645037.5,
  7645812.5,
  7696973.4375,
  7725545.3125,
  7728612.5,
  7729962.5,
  7776073.4375,
  7804296.875,
  8175401.5625,
  8280001.5625,
  8305990.625,
  8334287.5,
  8393864.0625,
  8404465.625,
  8409240.625,
  8619264.0625,
  8683720.3125,
  8760975.0,
  8819106.25,
  8841190.625,
  8843214.0625,
  8918756.25,
  8942096.875,
  9151096.875,
  9233778.125,
  9249223.4375,
  9249692.1875,
  9252159.375,
  9253204.6875,
  9261684.375,
  9265754.6875,
  9273990.625,
  9284509.375,
  9285464.0625,
  9285531.25,
  9286167.1875,
  9286200.0,
  9287225.0,
  9287421.875,
  9291946.875,
  9294339.0625,
  9304896.875,
  9310596.875,
  9316385.9375,
  9324654.6875,
  9347229.6875,
  9351431.25,
  9354643.75,
  9391928.125,
  9391979.6875,
  9401175.0,
  9433898.4375,
  9435915.625,
  9485206.25,
  9486193.75,
  9607493.75,
  9645565.625,
  9647695.3125,
  9679243.75,
  9686637.5,
  9698545.3125,
  9698628.125,
  9766025.0,
  9797554.6875,
  9858392.1875,
  9975434.375,
  9989701.5625,
  10006487.5,
  10033015.625,
  10033667.1875,
  10034017.1875,
  10037151.5625,
  10040265.625,
  10040507.8125,
  10043556.25,
  10045089.0625,
  10054309.375,
  10054693.75,
  10058800.0,
  10058854.6875,
  10083485.9375,
  10084515.625,
  10085392.1875,
  10086518.75,
  10086662.5,
  10087318.75,
  10087570.3125,
  10087718.75,
  10088165.625,
  10089090.625,
  10091775.0,
  10092259.375,
  10093300.0,
  10094692.1875,
  10096720.3125,
  10098984.375,
  10099103.125,
  10101800.0,
  10102284.375,
  10168135.9375,
  10186507.8125,
  10190696.875,
  10233775.0,
  10233779.6875,
  10253465.625,
  10301259.375,
  10335398.4375,
  10376289.0625,
  10379604.6875,
  10380329.6875,
  10381626.5625,
  10382620.3125,
  10382701.5625,
  10382726.5625,
  10383701.5625,
  10426812.5,
  10456092.1875,
  10483621.875,
  10491740.625,
  10491762.5,
  10492215.625,
  10492378.125,
  10495167.1875,
  10496246.875,
  10497623.4375,
  10498098.4375,
  10499415.625,
  10505012.5,
  10507031.25,
  10508006.25,
  10508201.5625,
  10508820.3125,
  10509150.0,
  10509257.8125,
  10509415.625,
  10509523.4375,
  10509535.9375,
  10509587.5,
  10509787.5,
  10509801.5625,
  10509968.75,
  10509970.3125,
  10510062.5,
  10510109.375,
  10510214.0625,
  10510256.25,
  10510417.1875,
  10510421.875,
  10510587.5,
  10510615.625,
  10510857.8125,
  10510912.5,
  10511028.125,
  10511320.3125,
  10511385.9375,
  10511423.4375,
  10511593.75,
  10512109.375,
  10512242.1875,
  10513353.125,
  10513729.6875,
  10515221.875,
  10520542.1875,
  10532196.875,
  10532740.625,
  10573425.0,
  10628867.1875,
  10631404.6875,
  10649678.125,
  10710756.25,
  10715351.5625,
  10726079.6875,
  10730131.25,
  10734125.0,
  10734398.4375,
  10741070.3125,
  10743126.5625,
  10743710.9375,
  10747173.4375,
  10758114.0625,
  10760028.125,
  10769201.5625,
  10807843.75,
  10814521.875,
  10821571.875,
  10828520.3125,
  10832264.0625,
  10856915.625,
  10864357.8125,
  10918748.4375,
  10918775.0,
  10954748.4375,
  10961151.5625,
  10965453.125,
  10966560.9375,
  10970059.375,
  10970348.4375,
  10970737.5,
  10971503.125,
  10972400.0,
  10972715.625,
  10972768.75,
  10972875.0,
  10974218.75,
  10975379.6875,
  10976545.3125,
  10977331.25,
  10980662.5,
  10984535.9375,
  10988976.5625,
  10990892.1875,
  10991371.875,
  10991457.8125,
  10991687.5,
  10991867.1875,
  10992367.1875,
  10992642.1875,
  10992817.1875,
  10992898.4375,
  10992950.0,
  10993178.125,
  10993254.6875,
  10993287.5,
  10993353.125,
  10993793.75,
  10993875.0,
  10994062.5,
  10994084.375,
  10994106.25,
  10994137.5,
  10994182.8125,
  10994262.5,
  10994292.1875,
  10994300.0,
  10994303.125,
  10994307.8125,
  10994329.6875,
  10994479.6875,
  10994479.6875,
  10994523.4375,
  10994528.125,
  10994540.625,
  10994584.375,
  10994607.8125,
  10994617.1875,
  10994671.875,
  10994725.0,
  10994729.6875,
  10994820.3125,
  10994821.875,
  10994828.125,
  10994843.75,
  10994854.6875,
  10994867.1875,
  10994925.0,
  10994939.0625,
  10994954.6875,
  10995060.9375,
  10995089.0625,
  10995154.6875,
  10995164.0625,
  10995184.375,
  10995184.375,
  10995187.5,
  10995190.625,
  10995204.6875,
  10995214.0625,
  10995285.9375,
  10995293.75,
  10995334.375,
  10995342.1875,
  10995360.9375,
  10995381.25,
  10995381.25,
  10995385.9375,
  10995385.9375,
  10995392.1875,
  10995418.75,
  10995421.875,
  10995421.875,
  10995437.5,
  10995440.625,
  10995451.5625,
  10995468.75,
  10995493.75,
  10995493.75,
  10995506.25,
  10995528.125,
  10995539.0625,
  10995643.75,
  10995671.875,
  10995689.0625,
  10995709.375,
  10995731.25,
  10995734.375,
  10995745.3125,
  10995746.875,
  10995793.75,
  10995815.625,
  10995818.75,
  10995832.8125,
  10995864.0625,
  10995907.8125,
  10995912.5,
  10995957.8125,
  10995992.1875,
  10996059.375,
  10996064.0625,
  10996065.625,
  10996078.125,
  10996081.25,
  10996103.125,
  10996109.375,
  10996242.1875,
  10996250.0,
  10996276.5625,
  10996434.375,
  10996500.0,
  10996567.1875,
  10996648.4375,
  10996657.8125,
  10996684.375,
  10996748.4375,
  10996760.9375,
  10996768.75,
  10996817.1875,
  10997046.875,
  10997157.8125,
  10997289.0625,
  10997295.3125,
  10997431.25,
  10997637.5,
  10997876.5625,
  10998407.8125,
  10998465.625,
  10998496.875,
  10998604.6875,
  10998659.375,
  10998684.375,
  10999493.75,
  10999617.1875,
  11000021.875,
  11000035.9375,
  11000095.3125,
  11000306.25,
  11000457.8125,
  11000640.625,
  11000695.3125,
  11001281.25,
  11001317.1875,
  11001328.125,
  11001356.25,
  11001365.625,
  11001407.8125,
  11001440.625,
  11001450.0,
  11001453.125,
  11001468.75,
  11001473.4375,
  11001492.1875,
  11001523.4375,
  11001534.375,
  11001579.6875,
  11001651.5625,
  11001657.8125,
  11001762.5,
  11001943.75,
  11001946.875,
  11001962.5,
  11001987.5,
  11002057.8125,
  11002075.0,
  11002089.0625,
  11002118.75,
  11002140.625,
  11002153.125,
  11002257.8125,
  11002426.5625,
  11002478.125,
  11002529.6875,
  11002589.0625,
  11002701.5625,
  11002831.25,
  11002884.375,
  11003018.75,
  11003201.5625,
  11003214.0625,
  11003342.1875,
  11003464.0625,
  11003537.5,
  11003557.8125,
  11003579.6875,
  11003767.1875,
  11003779.6875,
  11003921.875,
  11003951.5625,
  11003967.1875,
  11004085.9375,
  11004234.375,
  11004459.375,
  11004590.625,
  11004685.9375,
  11004693.75,
  11004867.1875,
  11005126.5625,
  11005179.6875,
  11005812.5,
  11006681.25,
  11006735.9375,
  11006989.0625,
  11007001.5625,
  11007140.625,
  11007159.375,
  11007373.4375,
  11007415.625,
  11007446.875,
  11007482.8125,
  11007506.25,
  11007707.8125,
  11007728.125,
  11007759.375,
  11007814.0625,
  11007818.75,
  11007857.8125,
  11007890.625,
  11007918.75,
  11007929.6875,
  11007942.1875,
  11007946.875,
  11007959.375,
  11007959.375,
  11007971.875,
  11007998.4375,
  11008078.125,
  11008095.3125,
  11008101.5625,
  11008154.6875,
  11008204.6875,
  11008214.0625,
  11008257.8125,
  11008290.625,
  11008360.9375,
  11008364.0625,
  11008420.3125,
  11008464.0625,
  11008528.125,
  11008581.25,
  11008656.25,
  11008714.0625,
  11008756.25,
  11008767.1875,
  11008804.6875,
  11008807.8125,
  11008826.5625,
  11008834.375,
  11008862.5,
  11008948.4375,
  11009015.625,
  11009023.4375,
  11009110.9375,
  11009171.875,
  11009203.125,
  11009260.9375,
  11009473.4375,
  11009476.5625,
  11012743.75,
  11014612.5,
  11015421.875,
  11015943.75,
  11015976.5625,
  11016662.5,
  11016925.0,
  11016929.6875,
  11017753.125,
  11018825.0,
  11020676.5625,
  11027318.75,
  11027320.3125,
  11028362.5,
  11028370.3125,
  11037832.8125,
  11038293.75,
  11046914.0625,
  11051267.1875,
  11051351.5625,
  11051782.8125,
  11052800.0,
  11052903.125,
  11055126.5625,
  11055343.75,
  11055925.0,
  11056851.5625,
  11057420.3125,
  11058104.6875,
  11058665.625,
  11060215.625,
  11062082.8125,
  11062314.0625,
  11063559.375,
  11070346.875,
  11070935.9375,
  11070993.75,
  11072857.8125,
  11074231.25,
  11081189.0625,
  11082985.9375,
  11083876.5625,
  151831164.0625],
 [14.522412378749344,
  7.379992778300393,
  20.64045855497977,
  5.332989456358101,
  6.602010812981603,
  9.893424004460336,
  5.589378615973168,
  110.90538478687675,
  11.622872544491031,
  239.58983928274074,
  38.329248057040125,
  30.156670690525814,
  18.312726415391275,
  13.213046940069017,
  55.499476705499426,
  12.635008758636628,
  11.823825567831367,
  30.870911107769977,
  18.281691392176455,
  10.706446490391627,
  31.559052674672955,
  15.895855634874952,
  98.36152771154717,
  12.037312005452641,
  29.274761241454954,
  13.630214704948756,
  21.219575188132197,
  6.892693108972041,
  21.781125514213358,
  13.931343136319423,
  11.954458275133598,
  66.73275424795852,
  5.568376154226905,
  23.72959956898433,
  41.37533968905696,
  5.931186402321224,
  12.167238277426103,
  25.025206949742767,
  17.938851569538727,
  24.966618836864853,
  6.212778421386231,
  11.368993598320168,
  6.399029897659369,
  41.36418256274479,
  36.82693051703742,
  144.3227768443324,
  16.808172091539692,
  5.194234037820371,
  15.826552768804975,
  57.989383274771626,
  8.818435944969853,
  65.68688413196178,
  44.64478270348457,
  26.037015799208202,
  11.706944921685144,
  49.79007227663506,
  6.111296736451919,
  19.975285437579615,
  38.627272412796614,
  20.775486773095025,
  62.67762721634249,
  101.35731056713595,
  34.74045931880837,
  81.32856589187013,
  8.043412303762292,
  33.54386198889696,
  51.51636391299296,
  13.607063649881042,
  15.14491516945016,
  72.2685194215291,
  38.248752075632794,
  20.365521491486042,
  22.43134182092655,
  16.904271993482574,
  5.138220103948399,
  6.277932633580717,
  53.86370259348103,
  30.056699337048926,
  45.11949815537223,
  53.09602526003558,
  7.824758730285324,
  25.136474572811654,
  56.965452627250045,
  8.687201982735573,
  6.938473423169307,
  5.061121097064851,
  133.32872733439592,
  17.765340732483008,
  5.4015125207279135,
  14.192864833554385,
  14.53935567800184,
  34.63345622001923,
  60.79644814664843,
  50.80526275469903,
  218.1899287082013,
  56.883072513345894,
  79.25660573943448,
  40.87695836720245,
  6.615884398178277,
  81.50185818548036,
  102.97346715112994,
  20.027907061482868,
  8.699758766044383,
  8.653194328714939,
  5.11359448265291,
  32.62721118079563,
  35.0449669205988,
  55.33110389141369,
  5.110818410734519,
  93.43291721065732,
  52.821372728578766,
  23.473756405170686,
  77.05490659116731,
  20.26188446311911,
  5.183666374246042,
  8.414982115015128,
  20.462390948402188,
  5.129229862486688,
  48.21827717754643,
  40.19935426473932,
  43.857260625794666,
  42.096901830625924,
  106.04273309945751,
  21.24207218911864,
  18.085464898067954,
  36.84685085045022,
  71.4040101997111,
  6.0098415035890405,
  21.322507728958445,
  70.88419239166693,
  21.676578034426587,
  5.943412747191946,
  35.61649819708623,
  38.9113423970728,
  91.44714429286708,
  25.58894565757637,
  22.828962666450703,
  14.075128963828705,
  45.87313652808515,
  28.78495013587629,
  41.77044657502613,
  19.618697033701633,
  88.97796121291506,
  21.760966990165194,
  11.21048337112779,
  15.358939869131913,
  12.189919634745847,
  78.24750960999646,
  26.712369809335325,
  15.944981423681638,
  242.83830116691715,
  7.529317333603413,
  105.02939197105052,
  133.92301032128358,
  25.295187320693053,
  23.365490163923603,
  43.029147650747824,
  131.39301469355243,
  97.31593549434014,
  87.521590483323,
  61.39647622025864,
  20.490722849759596,
  9.604603074561172,
  60.61784056738845,
  18.211658646080693,
  20.698077379687824,
  6.137167577835072,
  49.77853747039155,
  10.207888441844103,
  11.246917884478261,
  47.25324192104496,
  65.1293254388128,
  24.19459803784773,
  90.26254115307393,
  6.368356953115741,
  41.99809989349299,
  69.9163799693189,
  20.103812425098596,
  37.73681860791247,
  35.39716302438328,
  67.22827888123679,
  8.150677163539159,
  15.261292399905834,
  35.62851392946617,
  37.95106830557609,
  88.34180028037551,
  38.56939648758064,
  18.58338911410286,
  34.78728885198922,
  71.95705252916208,
  12.201135906447666,
  31.560836051382026,
  65.68865727732911,
  14.157884227898197,
  29.045566711075583,
  78.57121345446706,
  16.51805070585848,
  37.07751422170219,
  39.811363934069014,
  66.91778807051035,
  39.97472551380958,
  113.47137292174581,
  54.68893151273588,
  73.52122735313056,
  143.92317869111557,
  32.30473596794936,
  20.765415229711678,
  17.19495925591494,
  34.37018966297808,
  5.17054730848159,
  31.298277229468493,
  8.813936063854511,
  90.3984176820932,
  53.923879350752244,
  30.9926977350519,
  42.35809036129966,
  92.02353793919525,
  12.943498850171082,
  8.49106700204624,
  16.261490936825016,
  23.145593562653506,
  7.103240281158632,
  33.597323951427995,
  47.61768654145746,
  25.227027592436468,
  11.770444260630278,
  18.887715137995105,
  29.1570302577205,
  8.975879372478412,
  14.633814264315477,
  57.21331724648192,
  71.59894134821812,
  9.524298882366455,
  42.37782901406701,
  18.515633035913623,
  30.293176114957753,
  14.024719768557528,
  51.103840693512296,
  32.55384439271148,
  5.129229784348003,
  22.610802466056544,
  58.84651464631091,
  37.923204299143606,
  51.98735238355061,
  48.59291253744345,
  31.419781165825405,
  45.019563310617,
  54.23399907653448,
  50.17771796853461,
  8.379946136630355,
  29.252189393661375,
  24.86740550440888,
  15.949126032680825,
  97.33661159595455,
  15.561169269121493,
  12.384890255553714,
  13.8193893929785,
  49.83061648020405,
  32.56327712838208,
  58.8399828527068,
  237.85023961168335,
  7.7736183744922736,
  73.43005564031378,
  10.87298271981529,
  32.91130203303636,
  20.128954653718885,
  57.489619366928174,
  18.654524622780475,
  57.803252877912925,
  34.52633665409199,
  10.908016378593254,
  90.55986605844916,
  39.87954447189963,
  7.361030551778287,
  44.30483807294905,
  11.604315674162388,
  9.778760016575207,
  41.88627705406574,
  141.11519655060658,
  42.29486581924668,
  11.205805691542148,
  51.04132662734624,
  72.65724673866539,
  94.95910394293396,
  128.22150229315847,
  42.96803135892449,
  39.70589819102901,
  62.97577086631011,
  56.040448701044035,
  26.89076669720656,
  28.70643271042061,
  23.057130051898277,
  24.656002302565632,
  175.58727149813038,
  96.56489921977861,
  73.81481916990167,
  16.73537135017913,
  16.32927302516183,
  31.48858583932537,
  10.071401486318006,
  9.215875831464261,
  108.5710705269465,
  5.2549536851257725,
  32.248156290893974,
  33.957869202641014,
  32.964145367463125,
  40.573141224991396,
  10.213821364515116,
  75.20565362000208,
  20.199526191126964,
  8.340869583580611,
  11.707933262941273,
  47.18005474009257,
  22.87487543238023,
  70.92953956136007,
  13.495330328738577,
  12.48425162350785,
  34.22227895519063,
  102.38712895673116,
  43.83356990840067,
  88.34045937630535,
  25.395164809745893,
  18.7900577058126,
  74.24061031034293,
  30.916743372846877,
  8.226542878960796,
  5.598321495504657,
  26.368431469700823,
  105.39909001328233,
  64.54209124714222,
  13.068654335637554,
  64.25442580482877,
  11.154716673104005,
  5.092807836105357,
  6.484505066144988,
  49.718081750008245,
  65.17715470391873,
  47.25530126913611,
  43.63685743542191,
  7.045660350381905,
  43.07968407081705,
  5.532713800152551,
  86.55904776009969,
  35.797390368598755,
  11.291655863136437,
  33.31790632291547,
  6.580314312994005,
  69.01943688051287,
  12.552490670146248,
  21.30785772606113,
  14.446971888551234,
  22.60316753586698,
  5.727697313815782,
  51.09028267891261,
  22.94967542553259,
  87.79325523340307,
  6.14436683665391,
  16.856345530811385,
  5.413100848723027,
  56.06986456302683,
  63.13061269126626,
  44.49310387576003,
  50.28868702987505,
  59.08755163070293,
  10.63697377265278,
  7.139064264230046,
  45.19612244852938,
  50.52060738802095,
  72.70241379232286,
  6.41688728005241,
  73.87571274690183,
  83.42023447666307,
  79.86328094613414,
  24.26078634954178,
  11.271004173151912,
  6.6723420695943965,
  43.81958941820018,
  76.5907632961646,
  10.742845244445874,
  67.45521577952377,
  47.74656676021197,
  24.90214744210141,
  89.79048978290905,
  73.29307364989036,
  41.27020773928539,
  91.32238095999966,
  54.006891372899105,
  45.166022761317116,
  53.72552244495233,
  19.00416897676892,
  10.234956892287236,
  31.150256449074554,
  11.375835129324832,
  9.362030025921142,
  12.54604748665758,
  5.93530311637763,
  66.94272458132205,
  43.20429333193299,
  40.26482319304364,
  46.15367575597653,
  50.3842935491331,
  50.538533223905105,
  15.141073759084229,
  82.27864007528866,
  79.40186545318585,
  40.61956623732796,
  75.02755334530434,
  61.61021772617565,
  53.00910946636314,
  12.266750407379732,
  40.28557965374314,
  53.02860554010428,
  79.78618347284286,
  6.055262311275415,
  5.442949892378668,
  7.605455311997575,
  18.206643181895934,
  74.89510468217578,
  6.143096656032979,
  24.005593130002406,
  90.69291186560322,
  10.18503435706545,
  9.124724430604719,
  26.91195178047335,
  7.789713694249968,
  47.305292758201816,
  5.578261085101192,
  23.779627651714872,
  140.5131643852066,
  57.684767173655224,
  71.74922248413849,
  33.15557240230683,
  16.645703464523116,
  37.09978429089474,
  7.894001915171967,
  57.65397880195902,
  17.908313813093184,
  23.728435632052374,
  46.471619302091355,
  18.873742977152524,
  7.402586579082309,
  79.60676749710838,
  15.644307866849804,
  5.104546910250049,
  46.156325523124764,
  54.78105417115934,
  24.350339269192034,
  11.270600071092053,
  5.970962042509847,
  29.729815414153855,
  55.42787021239333,
  11.511741688248748,
  22.317023808498988,
  6.353873037019125,
  47.78721957131938,
  71.83927135959942,
  52.812010606745616,
  11.73396454901019,
  6.128699943164398,
  72.72429150818456,
  6.2695976683471235,
  21.574502587042417,
  58.62398073711046,
  8.91774185806324,
  10.527213957027836,
  35.63637968256597,
  13.006443075551829,
  9.973367310432584,
  33.96371940830099,
  9.014396067181705,
  28.32770657305099,
  26.724477154180327,
  27.413789320210032,
  5.51665082833289,
  25.807665048402825,
  20.628770656301246,
  41.40895835186115,
  33.669607133354376,
  18.839680007110477,
  10.146518890411087,
  100.24903500901154,
  8.97526399302785,
  46.12767367685594,
  44.8528162546497,
  9.678039386429303,
  74.69354912229892,
  23.688487664224038,
  17.783860431202612,
  48.15473603194975,
  112.6960902922857,
  42.28987689297926,
  37.83654284102215,
  47.870869933050024,
  51.77640768792982,
  20.551599721648564,
  25.732533245162077,
  22.586433402874867,
  24.771827004221464,
  27.0755747299956,
  132.04403199094142,
  10.831232520090527,
  17.74755825429985,
  59.12946582589614,
  28.91359278764569,
  7.774654378966613,
  7.51386631583635,
  43.02177169499613,
  23.293211590279945,
  56.976170899326235,
  123.80055689238287,
  46.023159613972936,
  37.01139820308629,
  22.55573510783127,
  39.56422824316637,
  50.54096740408959,
  29.83191350763388,
  7.760900921146966,
  10.161441222067927,
  20.659340900155392,
  42.49410864870074,
  61.48666526072617,
  6.533091924921097,
  19.567501002232024,
  14.841877657825389,
  57.217800581628865,
  36.16787675286571,
  14.163766541589036,
  70.89983490439415,
  5.845871744444563,
  8.006113356538455,
  62.8551049545252,
  55.018646640655724,
  40.15736045756692,
  51.06734032581019,
  15.947039948299961,
  8.747152424439077,
  20.40536033958457,
  5.968212555372477,
  88.80924461630083,
  7.638357147601251,
  28.259044901676248,
  18.56045379661847,
  58.95708731326854,
  46.13927630527151,
  75.56172059662609,
  35.3242777282921,
  20.570639023681895,
  20.7570338432579,
  7.14513306217448,
  6.459396846820025,
  49.984578783008516,
  18.044423213754854,
  14.291447011305111,
  81.06273308354982,
  33.29065497448798,
  12.846827789713984,
  69.14397568365047,
  12.16431403940248,
  11.578217509744455,
  12.299540943370044,
  15.082951351340133,
  5.8731566575797,
  71.67585563358021,
  12.079145396137118,
  5.82930794838644,
  8.357428911293017,
  45.287975315295256,
  14.34790704387889,
  9.626000493755786,
  65.64604644327783,
  70.29590801542925,
  8.488143072221453,
  18.414947886021597,
  52.61199662979216,
  19.70833541476396,
  23.972217143270708,
  50.165246691928864,
  11.696039049206354,
  9.90131697830581,
  109.82602011050268,
  52.06545966448377,
  8.866436019682371,
  17.827981300158722,
  66.59581942026331,
  14.983261828746626,
  14.567153196507343,
  5.552611535652301,
  7.056759238299655,
  99.52177236852084,
  5.194565927259784,
  61.09886400546656,
  61.431918301097504,
  37.2401508968075,
  30.250622325670115,
  40.12151363885099,
  9.592989421076886,
  5.2557981885610765,
  14.18399562882538,
  50.59566804869893,
  18.625192550667332,
  48.490680927867345,
  22.87134627150619,
  109.37044623412048,
  106.34108044245549,
  6.464987511644463,
  14.570446903717254,
  66.03175718868997,
  6.387340647990015,
  51.078559952495766,
  20.18824715923365,
  8.222716356372963,
  16.52572507429681,
  90.60396527912548,
  26.42656241481182,
  58.496289283191125,
  8.504836798964465,
  17.214170734135003,
  21.910773970132258,
  8.066199133290901,
  16.426475372947156,
  31.55928003978359,
  9.643594835024379,
  33.109470026900816,
  47.76880120063448,
  50.35714206248249,
  7.822592524870712,
  11.65826066479001,
  66.16738335509929,
  17.6019875375003,
  95.67159483511644,
  18.740950540638618,
  67.93370128213776,
  10.591953068950675,
  39.55022381204277,
  16.680459756820152,
  55.75997547727652,
  12.082275523474362,
  16.024726709903177,
  13.41461577722364,
  66.3707870939449,
  72.02886249385459,
  16.255262604737894,
  12.137290851081648,
  8.96625395430296,
  5.0437588969339755,
  16.08265457403393,
  13.951048604127413,
  11.973080085858133,
  56.65635095026821,
  5.975008985054715,
  66.00135414520204,
  70.58633068872265,
  68.45497708632652,
  38.88817008071914,
  20.903892418355476,
  73.420858115235,
  27.77954518009039,
  13.20145840645084,
  10.306507253234834,
  5.353568382699863,
  7.116308624888159,
  30.497547025942197,
  14.986564311815544,
  18.141442000641227,
  5.365899879173057,
  17.799242445545108,
  18.197223218896994,
  42.14215427227429])

Energy spectrum

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

([6181285.9375,
  6197781.25,
  6239173.4375,
  6247750.0,
  6248804.6875,
  6262335.9375,
  6272785.9375,
  6279440.625,
  6280354.6875,
  6286471.875,
  6286995.3125,
  6293664.0625,
  6295378.125,
  6295903.125,
  6314421.875,
  6353832.8125,
  6392850.0,
  6421156.25,
  6507553.125,
  6508446.875,
  6520773.4375,
  6536207.8125,
  6545845.3125,
  6552589.0625,
  6576779.6875,
  6580757.8125,
  6580773.4375,
  6595562.5,
  6629950.0,
  6631928.125,
  6666404.6875,
  6668471.875,
  6668517.1875,
  6669001.5625,
  6669265.625,
  6673604.6875,
  6675820.3125,
  6680060.9375,
  6691098.4375,
  6718495.3125,
  6723492.1875,
  6731370.3125,
  6733576.5625,
  6759671.875,
  6761328.125,
  6780406.25,
  6786860.9375,
  6822914.0625,
  6825746.875,
  6828890.625,
  6850154.6875,
  6850968.75,
  6850970.3125,
  6855337.5,
  6870803.125,
  6878707.8125,
  6882725.0,
  6895995.3125,
  6929346.875,
  6931592.1875,
  6934289.0625,
  6935462.5,
  6938654.6875,
  6939896.875,
  6960801.5625,
  6967254.6875,
  6993356.25,
  6994426.5625,
  6994682.8125,
  6998009.375,
  7004154.6875,
  7032918.75,
  7063692.1875,
  7097726.5625,
  7100550.0,
  7107506.25,
  7116837.5,
  7130043.75,
  7138718.75,
  7142651.5625,
  7149993.75,
  7169731.25,
  7170189.0625,
  7174621.875,
  7175293.75,
  7175346.875,
  7177920.3125,
  7177967.1875,
  7178021.875,
  7178173.4375,
  7179323.4375,
  7179867.1875,
  7180895.3125,
  7182092.1875,
  7212487.5,
  7215584.375,
  7217818.75,
  7229596.875,
  7229642.1875,
  7231545.3125,
  7242971.875,
  7263040.625,
  7263084.375,
  7268900.0,
  7270918.75,
  7289931.25,
  7291128.125,
  7292212.5,
  7292292.1875,
  7294689.0625,
  7295395.3125,
  7295718.75,
  7298059.375,
  7304248.4375,
  7309814.0625,
  7314185.9375,
  7319234.375,
  7327867.1875,
  7330954.6875,
  7331865.625,
  7333028.125,
  7335623.4375,
  7336623.4375,
  7336845.3125,
  7337079.6875,
  7337087.5,
  7337167.1875,
  7337795.3125,
  7339092.1875,
  7353020.3125,
  7356026.5625,
  7369845.3125,
  7371171.875,
  7373626.5625,
  7373648.4375,
  7377445.3125,
  7416301.5625,
  7433556.25,
  7440907.8125,
  7450426.5625,
  7451001.5625,
  7452140.625,
  7454321.875,
  7456137.5,
  7467356.25,
  7467365.625,
  7489384.375,
  7490357.8125,
  7490390.625,
  7503518.75,
  7509009.375,
  7509067.1875,
  7530709.375,
  7539157.8125,
  7548681.25,
  7581431.25,
  7587604.6875,
  7587640.625,
  7595981.25,
  7644993.75,
  7645037.5,
  7645812.5,
  7696973.4375,
  7725545.3125,
  7728612.5,
  7729962.5,
  7776073.4375,
  7804296.875,
  8175401.5625,
  8280001.5625,
  8305990.625,
  8334287.5,
  8393864.0625,
  8404465.625,
  8409240.625,
  8619264.0625,
  8683720.3125,
  8760975.0,
  8819106.25,
  8841190.625,
  8843214.0625,
  8918756.25,
  8942096.875,
  9151096.875,
  9233778.125,
  9249223.4375,
  9249692.1875,
  9252159.375,
  9253204.6875,
  9261684.375,
  9265754.6875,
  9273990.625,
  9284509.375,
  9285464.0625,
  9285531.25,
  9286167.1875,
  9286200.0,
  9287225.0,
  9287421.875,
  9291946.875,
  9294339.0625,
  9304896.875,
  9310596.875,
  9316385.9375,
  9324654.6875,
  9347229.6875,
  9351431.25,
  9354643.75,
  9391928.125,
  9391979.6875,
  9401175.0,
  9433898.4375,
  9435915.625,
  9485206.25,
  9486193.75,
  9607493.75,
  9645565.625,
  9647695.3125,
  9679243.75,
  9686637.5,
  9698545.3125,
  9698628.125,
  9766025.0,
  9797554.6875,
  9858392.1875,
  9975434.375,
  9989701.5625,
  10006487.5,
  10033015.625,
  10033667.1875,
  10034017.1875,
  10037151.5625,
  10040265.625,
  10040507.8125,
  10043556.25,
  10045089.0625,
  10054309.375,
  10054693.75,
  10058800.0,
  10058854.6875,
  10083485.9375,
  10084515.625,
  10085392.1875,
  10086518.75,
  10086662.5,
  10087318.75,
  10087570.3125,
  10087718.75,
  10088165.625,
  10089090.625,
  10091775.0,
  10092259.375,
  10093300.0,
  10094692.1875,
  10096720.3125,
  10098984.375,
  10099103.125,
  10101800.0,
  10102284.375,
  10168135.9375,
  10186507.8125,
  10190696.875,
  10233775.0,
  10233779.6875,
  10253465.625,
  10301259.375,
  10335398.4375,
  10376289.0625,
  10379604.6875,
  10380329.6875,
  10381626.5625,
  10382620.3125,
  10382701.5625,
  10382726.5625,
  10383701.5625,
  10426812.5,
  10456092.1875,
  10483621.875,
  10491740.625,
  10491762.5,
  10492215.625,
  10492378.125,
  10495167.1875,
  10496246.875,
  10497623.4375,
  10498098.4375,
  10499415.625,
  10505012.5,
  10507031.25,
  10508006.25,
  10508201.5625,
  10508820.3125,
  10509150.0,
  10509257.8125,
  10509415.625,
  10509523.4375,
  10509535.9375,
  10509587.5,
  10509787.5,
  10509801.5625,
  10509968.75,
  10509970.3125,
  10510062.5,
  10510109.375,
  10510214.0625,
  10510256.25,
  10510417.1875,
  10510421.875,
  10510587.5,
  10510615.625,
  10510857.8125,
  10510912.5,
  10511028.125,
  10511320.3125,
  10511385.9375,
  10511423.4375,
  10511593.75,
  10512109.375,
  10512242.1875,
  10513353.125,
  10513729.6875,
  10515221.875,
  10520542.1875,
  10532196.875,
  10532740.625,
  10573425.0,
  10628867.1875,
  10631404.6875,
  10649678.125,
  10710756.25,
  10715351.5625,
  10726079.6875,
  10730131.25,
  10734125.0,
  10734398.4375,
  10741070.3125,
  10743126.5625,
  10743710.9375,
  10747173.4375,
  10758114.0625,
  10760028.125,
  10769201.5625,
  10807843.75,
  10814521.875,
  10821571.875,
  10828520.3125,
  10832264.0625,
  10856915.625,
  10864357.8125,
  10918748.4375,
  10918775.0,
  10954748.4375,
  10961151.5625,
  10965453.125,
  10966560.9375,
  10970059.375,
  10970348.4375,
  10970737.5,
  10971503.125,
  10972400.0,
  10972715.625,
  10972768.75,
  10972875.0,
  10974218.75,
  10975379.6875,
  10976545.3125,
  10977331.25,
  10980662.5,
  10984535.9375,
  10988976.5625,
  10990892.1875,
  10991371.875,
  10991457.8125,
  10991687.5,
  10991867.1875,
  10992367.1875,
  10992642.1875,
  10992817.1875,
  10992898.4375,
  10992950.0,
  10993178.125,
  10993254.6875,
  10993287.5,
  10993353.125,
  10993793.75,
  10993875.0,
  10994062.5,
  10994084.375,
  10994106.25,
  10994137.5,
  10994182.8125,
  10994262.5,
  10994292.1875,
  10994300.0,
  10994303.125,
  10994307.8125,
  10994329.6875,
  10994479.6875,
  10994479.6875,
  10994523.4375,
  10994528.125,
  10994540.625,
  10994584.375,
  10994607.8125,
  10994617.1875,
  10994671.875,
  10994725.0,
  10994729.6875,
  10994820.3125,
  10994821.875,
  10994828.125,
  10994843.75,
  10994854.6875,
  10994867.1875,
  10994925.0,
  10994939.0625,
  10994954.6875,
  10995060.9375,
  10995089.0625,
  10995154.6875,
  10995164.0625,
  10995184.375,
  10995184.375,
  10995187.5,
  10995190.625,
  10995204.6875,
  10995214.0625,
  10995285.9375,
  10995293.75,
  10995334.375,
  10995342.1875,
  10995360.9375,
  10995381.25,
  10995381.25,
  10995385.9375,
  10995385.9375,
  10995392.1875,
  10995418.75,
  10995421.875,
  10995421.875,
  10995437.5,
  10995440.625,
  10995451.5625,
  10995468.75,
  10995493.75,
  10995493.75,
  10995506.25,
  10995528.125,
  10995539.0625,
  10995643.75,
  10995671.875,
  10995689.0625,
  10995709.375,
  10995731.25,
  10995734.375,
  10995745.3125,
  10995746.875,
  10995793.75,
  10995815.625,
  10995818.75,
  10995832.8125,
  10995864.0625,
  10995907.8125,
  10995912.5,
  10995957.8125,
  10995992.1875,
  10996059.375,
  10996064.0625,
  10996065.625,
  10996078.125,
  10996081.25,
  10996103.125,
  10996109.375,
  10996242.1875,
  10996250.0,
  10996276.5625,
  10996434.375,
  10996500.0,
  10996567.1875,
  10996648.4375,
  10996657.8125,
  10996684.375,
  10996748.4375,
  10996760.9375,
  10996768.75,
  10996817.1875,
  10997046.875,
  10997157.8125,
  10997289.0625,
  10997295.3125,
  10997431.25,
  10997637.5,
  10997876.5625,
  10998407.8125,
  10998465.625,
  10998496.875,
  10998604.6875,
  10998659.375,
  10998684.375,
  10999493.75,
  10999617.1875,
  11000021.875,
  11000035.9375,
  11000095.3125,
  11000306.25,
  11000457.8125,
  11000640.625,
  11000695.3125,
  11001281.25,
  11001317.1875,
  11001328.125,
  11001356.25,
  11001365.625,
  11001407.8125,
  11001440.625,
  11001450.0,
  11001453.125,
  11001468.75,
  11001473.4375,
  11001492.1875,
  11001523.4375,
  11001534.375,
  11001579.6875,
  11001651.5625,
  11001657.8125,
  11001762.5,
  11001943.75,
  11001946.875,
  11001962.5,
  11001987.5,
  11002057.8125,
  11002075.0,
  11002089.0625,
  11002118.75,
  11002140.625,
  11002153.125,
  11002257.8125,
  11002426.5625,
  11002478.125,
  11002529.6875,
  11002589.0625,
  11002701.5625,
  11002831.25,
  11002884.375,
  11003018.75,
  11003201.5625,
  11003214.0625,
  11003342.1875,
  11003464.0625,
  11003537.5,
  11003557.8125,
  11003579.6875,
  11003767.1875,
  11003779.6875,
  11003921.875,
  11003951.5625,
  11003967.1875,
  11004085.9375,
  11004234.375,
  11004459.375,
  11004590.625,
  11004685.9375,
  11004693.75,
  11004867.1875,
  11005126.5625,
  11005179.6875,
  11005812.5,
  11006681.25,
  11006735.9375,
  11006989.0625,
  11007001.5625,
  11007140.625,
  11007159.375,
  11007373.4375,
  11007415.625,
  11007446.875,
  11007482.8125,
  11007506.25,
  11007707.8125,
  11007728.125,
  11007759.375,
  11007814.0625,
  11007818.75,
  11007857.8125,
  11007890.625,
  11007918.75,
  11007929.6875,
  11007942.1875,
  11007946.875,
  11007959.375,
  11007959.375,
  11007971.875,
  11007998.4375,
  11008078.125,
  11008095.3125,
  11008101.5625,
  11008154.6875,
  11008204.6875,
  11008214.0625,
  11008257.8125,
  11008290.625,
  11008360.9375,
  11008364.0625,
  11008420.3125,
  11008464.0625,
  11008528.125,
  11008581.25,
  11008656.25,
  11008714.0625,
  11008756.25,
  11008767.1875,
  11008804.6875,
  11008807.8125,
  11008826.5625,
  11008834.375,
  11008862.5,
  11008948.4375,
  11009015.625,
  11009023.4375,
  11009110.9375,
  11009171.875,
  11009203.125,
  11009260.9375,
  11009473.4375,
  11009476.5625,
  11012743.75,
  11014612.5,
  11015421.875,
  11015943.75,
  11015976.5625,
  11016662.5,
  11016925.0,
  11016929.6875,
  11017753.125,
  11018825.0,
  11020676.5625,
  11027318.75,
  11027320.3125,
  11028362.5,
  11028370.3125,
  11037832.8125,
  11038293.75,
  11046914.0625,
  11051267.1875,
  11051351.5625,
  11051782.8125,
  11052800.0,
  11052903.125,
  11055126.5625,
  11055343.75,
  11055925.0,
  11056851.5625,
  11057420.3125,
  11058104.6875,
  11058665.625,
  11060215.625,
  11062082.8125,
  11062314.0625,
  11063559.375,
  11070346.875,
  11070935.9375,
  11070993.75,
  11072857.8125,
  11074231.25,
  11081189.0625,
  11082985.9375,
  11083876.5625,
  151831164.0625],
 [14.522412378749344,
  7.379992778300393,
  20.64045855497977,
  5.332989456358101,
  6.602010812981603,
  9.893424004460336,
  5.589378615973168,
  110.90538478687675,
  11.622872544491031,
  239.58983928274074,
  38.329248057040125,
  30.156670690525814,
  18.312726415391275,
  13.213046940069017,
  55.499476705499426,
  12.635008758636628,
  11.823825567831367,
  30.870911107769977,
  18.281691392176455,
  10.706446490391627,
  31.559052674672955,
  15.895855634874952,
  98.36152771154717,
  12.037312005452641,
  29.274761241454954,
  13.630214704948756,
  21.219575188132197,
  6.892693108972041,
  21.781125514213358,
  13.931343136319423,
  11.954458275133598,
  66.73275424795852,
  5.568376154226905,
  23.72959956898433,
  41.37533968905696,
  5.931186402321224,
  12.167238277426103,
  25.025206949742767,
  17.938851569538727,
  24.966618836864853,
  6.212778421386231,
  11.368993598320168,
  6.399029897659369,
  41.36418256274479,
  36.82693051703742,
  144.3227768443324,
  16.808172091539692,
  5.194234037820371,
  15.826552768804975,
  57.989383274771626,
  8.818435944969853,
  65.68688413196178,
  44.64478270348457,
  26.037015799208202,
  11.706944921685144,
  49.79007227663506,
  6.111296736451919,
  19.975285437579615,
  38.627272412796614,
  20.775486773095025,
  62.67762721634249,
  101.35731056713595,
  34.74045931880837,
  81.32856589187013,
  8.043412303762292,
  33.54386198889696,
  51.51636391299296,
  13.607063649881042,
  15.14491516945016,
  72.2685194215291,
  38.248752075632794,
  20.365521491486042,
  22.43134182092655,
  16.904271993482574,
  5.138220103948399,
  6.277932633580717,
  53.86370259348103,
  30.056699337048926,
  45.11949815537223,
  53.09602526003558,
  7.824758730285324,
  25.136474572811654,
  56.965452627250045,
  8.687201982735573,
  6.938473423169307,
  5.061121097064851,
  133.32872733439592,
  17.765340732483008,
  5.4015125207279135,
  14.192864833554385,
  14.53935567800184,
  34.63345622001923,
  60.79644814664843,
  50.80526275469903,
  218.1899287082013,
  56.883072513345894,
  79.25660573943448,
  40.87695836720245,
  6.615884398178277,
  81.50185818548036,
  102.97346715112994,
  20.027907061482868,
  8.699758766044383,
  8.653194328714939,
  5.11359448265291,
  32.62721118079563,
  35.0449669205988,
  55.33110389141369,
  5.110818410734519,
  93.43291721065732,
  52.821372728578766,
  23.473756405170686,
  77.05490659116731,
  20.26188446311911,
  5.183666374246042,
  8.414982115015128,
  20.462390948402188,
  5.129229862486688,
  48.21827717754643,
  40.19935426473932,
  43.857260625794666,
  42.096901830625924,
  106.04273309945751,
  21.24207218911864,
  18.085464898067954,
  36.84685085045022,
  71.4040101997111,
  6.0098415035890405,
  21.322507728958445,
  70.88419239166693,
  21.676578034426587,
  5.943412747191946,
  35.61649819708623,
  38.9113423970728,
  91.44714429286708,
  25.58894565757637,
  22.828962666450703,
  14.075128963828705,
  45.87313652808515,
  28.78495013587629,
  41.77044657502613,
  19.618697033701633,
  88.97796121291506,
  21.760966990165194,
  11.21048337112779,
  15.358939869131913,
  12.189919634745847,
  78.24750960999646,
  26.712369809335325,
  15.944981423681638,
  242.83830116691715,
  7.529317333603413,
  105.02939197105052,
  133.92301032128358,
  25.295187320693053,
  23.365490163923603,
  43.029147650747824,
  131.39301469355243,
  97.31593549434014,
  87.521590483323,
  61.39647622025864,
  20.490722849759596,
  9.604603074561172,
  60.61784056738845,
  18.211658646080693,
  20.698077379687824,
  6.137167577835072,
  49.77853747039155,
  10.207888441844103,
  11.246917884478261,
  47.25324192104496,
  65.1293254388128,
  24.19459803784773,
  90.26254115307393,
  6.368356953115741,
  41.99809989349299,
  69.9163799693189,
  20.103812425098596,
  37.73681860791247,
  35.39716302438328,
  67.22827888123679,
  8.150677163539159,
  15.261292399905834,
  35.62851392946617,
  37.95106830557609,
  88.34180028037551,
  38.56939648758064,
  18.58338911410286,
  34.78728885198922,
  71.95705252916208,
  12.201135906447666,
  31.560836051382026,
  65.68865727732911,
  14.157884227898197,
  29.045566711075583,
  78.57121345446706,
  16.51805070585848,
  37.07751422170219,
  39.811363934069014,
  66.91778807051035,
  39.97472551380958,
  113.47137292174581,
  54.68893151273588,
  73.52122735313056,
  143.92317869111557,
  32.30473596794936,
  20.765415229711678,
  17.19495925591494,
  34.37018966297808,
  5.17054730848159,
  31.298277229468493,
  8.813936063854511,
  90.3984176820932,
  53.923879350752244,
  30.9926977350519,
  42.35809036129966,
  92.02353793919525,
  12.943498850171082,
  8.49106700204624,
  16.261490936825016,
  23.145593562653506,
  7.103240281158632,
  33.597323951427995,
  47.61768654145746,
  25.227027592436468,
  11.770444260630278,
  18.887715137995105,
  29.1570302577205,
  8.975879372478412,
  14.633814264315477,
  57.21331724648192,
  71.59894134821812,
  9.524298882366455,
  42.37782901406701,
  18.515633035913623,
  30.293176114957753,
  14.024719768557528,
  51.103840693512296,
  32.55384439271148,
  5.129229784348003,
  22.610802466056544,
  58.84651464631091,
  37.923204299143606,
  51.98735238355061,
  48.59291253744345,
  31.419781165825405,
  45.019563310617,
  54.23399907653448,
  50.17771796853461,
  8.379946136630355,
  29.252189393661375,
  24.86740550440888,
  15.949126032680825,
  97.33661159595455,
  15.561169269121493,
  12.384890255553714,
  13.8193893929785,
  49.83061648020405,
  32.56327712838208,
  58.8399828527068,
  237.85023961168335,
  7.7736183744922736,
  73.43005564031378,
  10.87298271981529,
  32.91130203303636,
  20.128954653718885,
  57.489619366928174,
  18.654524622780475,
  57.803252877912925,
  34.52633665409199,
  10.908016378593254,
  90.55986605844916,
  39.87954447189963,
  7.361030551778287,
  44.30483807294905,
  11.604315674162388,
  9.778760016575207,
  41.88627705406574,
  141.11519655060658,
  42.29486581924668,
  11.205805691542148,
  51.04132662734624,
  72.65724673866539,
  94.95910394293396,
  128.22150229315847,
  42.96803135892449,
  39.70589819102901,
  62.97577086631011,
  56.040448701044035,
  26.89076669720656,
  28.70643271042061,
  23.057130051898277,
  24.656002302565632,
  175.58727149813038,
  96.56489921977861,
  73.81481916990167,
  16.73537135017913,
  16.32927302516183,
  31.48858583932537,
  10.071401486318006,
  9.215875831464261,
  108.5710705269465,
  5.2549536851257725,
  32.248156290893974,
  33.957869202641014,
  32.964145367463125,
  40.573141224991396,
  10.213821364515116,
  75.20565362000208,
  20.199526191126964,
  8.340869583580611,
  11.707933262941273,
  47.18005474009257,
  22.87487543238023,
  70.92953956136007,
  13.495330328738577,
  12.48425162350785,
  34.22227895519063,
  102.38712895673116,
  43.83356990840067,
  88.34045937630535,
  25.395164809745893,
  18.7900577058126,
  74.24061031034293,
  30.916743372846877,
  8.226542878960796,
  5.598321495504657,
  26.368431469700823,
  105.39909001328233,
  64.54209124714222,
  13.068654335637554,
  64.25442580482877,
  11.154716673104005,
  5.092807836105357,
  6.484505066144988,
  49.718081750008245,
  65.17715470391873,
  47.25530126913611,
  43.63685743542191,
  7.045660350381905,
  43.07968407081705,
  5.532713800152551,
  86.55904776009969,
  35.797390368598755,
  11.291655863136437,
  33.31790632291547,
  6.580314312994005,
  69.01943688051287,
  12.552490670146248,
  21.30785772606113,
  14.446971888551234,
  22.60316753586698,
  5.727697313815782,
  51.09028267891261,
  22.94967542553259,
  87.79325523340307,
  6.14436683665391,
  16.856345530811385,
  5.413100848723027,
  56.06986456302683,
  63.13061269126626,
  44.49310387576003,
  50.28868702987505,
  59.08755163070293,
  10.63697377265278,
  7.139064264230046,
  45.19612244852938,
  50.52060738802095,
  72.70241379232286,
  6.41688728005241,
  73.87571274690183,
  83.42023447666307,
  79.86328094613414,
  24.26078634954178,
  11.271004173151912,
  6.6723420695943965,
  43.81958941820018,
  76.5907632961646,
  10.742845244445874,
  67.45521577952377,
  47.74656676021197,
  24.90214744210141,
  89.79048978290905,
  73.29307364989036,
  41.27020773928539,
  91.32238095999966,
  54.006891372899105,
  45.166022761317116,
  53.72552244495233,
  19.00416897676892,
  10.234956892287236,
  31.150256449074554,
  11.375835129324832,
  9.362030025921142,
  12.54604748665758,
  5.93530311637763,
  66.94272458132205,
  43.20429333193299,
  40.26482319304364,
  46.15367575597653,
  50.3842935491331,
  50.538533223905105,
  15.141073759084229,
  82.27864007528866,
  79.40186545318585,
  40.61956623732796,
  75.02755334530434,
  61.61021772617565,
  53.00910946636314,
  12.266750407379732,
  40.28557965374314,
  53.02860554010428,
  79.78618347284286,
  6.055262311275415,
  5.442949892378668,
  7.605455311997575,
  18.206643181895934,
  74.89510468217578,
  6.143096656032979,
  24.005593130002406,
  90.69291186560322,
  10.18503435706545,
  9.124724430604719,
  26.91195178047335,
  7.789713694249968,
  47.305292758201816,
  5.578261085101192,
  23.779627651714872,
  140.5131643852066,
  57.684767173655224,
  71.74922248413849,
  33.15557240230683,
  16.645703464523116,
  37.09978429089474,
  7.894001915171967,
  57.65397880195902,
  17.908313813093184,
  23.728435632052374,
  46.471619302091355,
  18.873742977152524,
  7.402586579082309,
  79.60676749710838,
  15.644307866849804,
  5.104546910250049,
  46.156325523124764,
  54.78105417115934,
  24.350339269192034,
  11.270600071092053,
  5.970962042509847,
  29.729815414153855,
  55.42787021239333,
  11.511741688248748,
  22.317023808498988,
  6.353873037019125,
  47.78721957131938,
  71.83927135959942,
  52.812010606745616,
  11.73396454901019,
  6.128699943164398,
  72.72429150818456,
  6.2695976683471235,
  21.574502587042417,
  58.62398073711046,
  8.91774185806324,
  10.527213957027836,
  35.63637968256597,
  13.006443075551829,
  9.973367310432584,
  33.96371940830099,
  9.014396067181705,
  28.32770657305099,
  26.724477154180327,
  27.413789320210032,
  5.51665082833289,
  25.807665048402825,
  20.628770656301246,
  41.40895835186115,
  33.669607133354376,
  18.839680007110477,
  10.146518890411087,
  100.24903500901154,
  8.97526399302785,
  46.12767367685594,
  44.8528162546497,
  9.678039386429303,
  74.69354912229892,
  23.688487664224038,
  17.783860431202612,
  48.15473603194975,
  112.6960902922857,
  42.28987689297926,
  37.83654284102215,
  47.870869933050024,
  51.77640768792982,
  20.551599721648564,
  25.732533245162077,
  22.586433402874867,
  24.771827004221464,
  27.0755747299956,
  132.04403199094142,
  10.831232520090527,
  17.74755825429985,
  59.12946582589614,
  28.91359278764569,
  7.774654378966613,
  7.51386631583635,
  43.02177169499613,
  23.293211590279945,
  56.976170899326235,
  123.80055689238287,
  46.023159613972936,
  37.01139820308629,
  22.55573510783127,
  39.56422824316637,
  50.54096740408959,
  29.83191350763388,
  7.760900921146966,
  10.161441222067927,
  20.659340900155392,
  42.49410864870074,
  61.48666526072617,
  6.533091924921097,
  19.567501002232024,
  14.841877657825389,
  57.217800581628865,
  36.16787675286571,
  14.163766541589036,
  70.89983490439415,
  5.845871744444563,
  8.006113356538455,
  62.8551049545252,
  55.018646640655724,
  40.15736045756692,
  51.06734032581019,
  15.947039948299961,
  8.747152424439077,
  20.40536033958457,
  5.968212555372477,
  88.80924461630083,
  7.638357147601251,
  28.259044901676248,
  18.56045379661847,
  58.95708731326854,
  46.13927630527151,
  75.56172059662609,
  35.3242777282921,
  20.570639023681895,
  20.7570338432579,
  7.14513306217448,
  6.459396846820025,
  49.984578783008516,
  18.044423213754854,
  14.291447011305111,
  81.06273308354982,
  33.29065497448798,
  12.846827789713984,
  69.14397568365047,
  12.16431403940248,
  11.578217509744455,
  12.299540943370044,
  15.082951351340133,
  5.8731566575797,
  71.67585563358021,
  12.079145396137118,
  5.82930794838644,
  8.357428911293017,
  45.287975315295256,
  14.34790704387889,
  9.626000493755786,
  65.64604644327783,
  70.29590801542925,
  8.488143072221453,
  18.414947886021597,
  52.61199662979216,
  19.70833541476396,
  23.972217143270708,
  50.165246691928864,
  11.696039049206354,
  9.90131697830581,
  109.82602011050268,
  52.06545966448377,
  8.866436019682371,
  17.827981300158722,
  66.59581942026331,
  14.983261828746626,
  14.567153196507343,
  5.552611535652301,
  7.056759238299655,
  99.52177236852084,
  5.194565927259784,
  61.09886400546656,
  61.431918301097504,
  37.2401508968075,
  30.250622325670115,
  40.12151363885099,
  9.592989421076886,
  5.2557981885610765,
  14.18399562882538,
  50.59566804869893,
  18.625192550667332,
  48.490680927867345,
  22.87134627150619,
  109.37044623412048,
  106.34108044245549,
  6.464987511644463,
  14.570446903717254,
  66.03175718868997,
  6.387340647990015,
  51.078559952495766,
  20.18824715923365,
  8.222716356372963,
  16.52572507429681,
  90.60396527912548,
  26.42656241481182,
  58.496289283191125,
  8.504836798964465,
  17.214170734135003,
  21.910773970132258,
  8.066199133290901,
  16.426475372947156,
  31.55928003978359,
  9.643594835024379,
  33.109470026900816,
  47.76880120063448,
  50.35714206248249,
  7.822592524870712,
  11.65826066479001,
  66.16738335509929,
  17.6019875375003,
  95.67159483511644,
  18.740950540638618,
  67.93370128213776,
  10.591953068950675,
  39.55022381204277,
  16.680459756820152,
  55.75997547727652,
  12.082275523474362,
  16.024726709903177,
  13.41461577722364,
  66.3707870939449,
  72.02886249385459,
  16.255262604737894,
  12.137290851081648,
  8.96625395430296,
  5.0437588969339755,
  16.08265457403393,
  13.951048604127413,
  11.973080085858133,
  56.65635095026821,
  5.975008985054715,
  66.00135414520204,
  70.58633068872265,
  68.45497708632652,
  38.88817008071914,
  20.903892418355476,
  73.420858115235,
  27.77954518009039,
  13.20145840645084,
  10.306507253234834,
  5.353568382699863,
  7.116308624888159,
  30.497547025942197,
  14.986564311815544,
  18.141442000641227,
  5.365899879173057,
  17.799242445545108,
  18.197223218896994,
  42.14215427227429])

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

([6181285.9375,
  6197781.25,
  6239173.4375,
  6247750.0,
  6248804.6875,
  6262335.9375,
  6272785.9375,
  6279440.625,
  6280354.6875,
  6286471.875,
  6286995.3125,
  6293664.0625,
  6295378.125,
  6295903.125,
  6314421.875,
  6353832.8125,
  6392850.0,
  6421156.25,
  6507553.125,
  6508446.875,
  6520773.4375,
  6536207.8125,
  6545845.3125,
  6552589.0625,
  6576779.6875,
  6580757.8125,
  6580773.4375,
  6595562.5,
  6629950.0,
  6631928.125,
  6666404.6875,
  6668471.875,
  6668517.1875,
  6669001.5625,
  6669265.625,
  6673604.6875,
  6675820.3125,
  6680060.9375,
  6691098.4375,
  6718495.3125,
  6723492.1875,
  6731370.3125,
  6733576.5625,
  6759671.875,
  6761328.125,
  6780406.25,
  6786860.9375,
  6822914.0625,
  6825746.875,
  6828890.625,
  6850154.6875,
  6850968.75,
  6850970.3125,
  6855337.5,
  6870803.125,
  6878707.8125,
  6882725.0,
  6895995.3125,
  6929346.875,
  6931592.1875,
  6934289.0625,
  6935462.5,
  6938654.6875,
  6939896.875,
  6960801.5625,
  6967254.6875,
  6993356.25,
  6994426.5625,
  6994682.8125,
  6998009.375,
  7004154.6875,
  7032918.75,
  7063692.1875,
  7097726.5625,
  7100550.0,
  7107506.25,
  7116837.5,
  7130043.75,
  7138718.75,
  7142651.5625,
  7149993.75,
  7169731.25,
  7170189.0625,
  7174621.875,
  7175293.75,
  7175346.875,
  7177920.3125,
  7177967.1875,
  7178021.875,
  7178173.4375,
  7179323.4375,
  7179867.1875,
  7180895.3125,
  7182092.1875,
  7212487.5,
  7215584.375,
  7217818.75,
  7229596.875,
  7229642.1875,
  7231545.3125,
  7242971.875,
  7263040.625,
  7263084.375,
  7268900.0,
  7270918.75,
  7289931.25,
  7291128.125,
  7292212.5,
  7292292.1875,
  7294689.0625,
  7295395.3125,
  7295718.75,
  7298059.375,
  7304248.4375,
  7309814.0625,
  7314185.9375,
  7319234.375,
  7327867.1875,
  7330954.6875,
  7331865.625,
  7333028.125,
  7335623.4375,
  7336623.4375,
  7336845.3125,
  7337079.6875,
  7337087.5,
  7337167.1875,
  7337795.3125,
  7339092.1875,
  7353020.3125,
  7356026.5625,
  7369845.3125,
  7371171.875,
  7373626.5625,
  7373648.4375,
  7377445.3125,
  7416301.5625,
  7433556.25,
  7440907.8125,
  7450426.5625,
  7451001.5625,
  7452140.625,
  7454321.875,
  7456137.5,
  7467356.25,
  7467365.625,
  7489384.375,
  7490357.8125,
  7490390.625,
  7503518.75,
  7509009.375,
  7509067.1875,
  7530709.375,
  7539157.8125,
  7548681.25,
  7581431.25,
  7587604.6875,
  7587640.625,
  7595981.25,
  7644993.75,
  7645037.5,
  7645812.5,
  7696973.4375,
  7725545.3125,
  7728612.5,
  7729962.5,
  7776073.4375,
  7804296.875,
  8175401.5625,
  8280001.5625,
  8305990.625,
  8334287.5,
  8393864.0625,
  8404465.625,
  8409240.625,
  8619264.0625,
  8683720.3125,
  8760975.0,
  8819106.25,
  8841190.625,
  8843214.0625,
  8918756.25,
  8942096.875,
  9151096.875,
  9233778.125,
  9249223.4375,
  9249692.1875,
  9252159.375,
  9253204.6875,
  9261684.375,
  9265754.6875,
  9273990.625,
  9284509.375,
  9285464.0625,
  9285531.25,
  9286167.1875,
  9286200.0,
  9287225.0,
  9287421.875,
  9291946.875,
  9294339.0625,
  9304896.875,
  9310596.875,
  9316385.9375,
  9324654.6875,
  9347229.6875,
  9351431.25,
  9354643.75,
  9391928.125,
  9391979.6875,
  9401175.0,
  9433898.4375,
  9435915.625,
  9485206.25,
  9486193.75,
  9607493.75,
  9645565.625,
  9647695.3125,
  9679243.75,
  9686637.5,
  9698545.3125,
  9698628.125,
  9766025.0,
  9797554.6875,
  9858392.1875,
  9975434.375,
  9989701.5625,
  10006487.5,
  10033015.625,
  10033667.1875,
  10034017.1875,
  10037151.5625,
  10040265.625,
  10040507.8125,
  10043556.25,
  10045089.0625,
  10054309.375,
  10054693.75,
  10058800.0,
  10058854.6875,
  10083485.9375,
  10084515.625,
  10085392.1875,
  10086518.75,
  10086662.5,
  10087318.75,
  10087570.3125,
  10087718.75,
  10088165.625,
  10089090.625,
  10091775.0,
  10092259.375,
  10093300.0,
  10094692.1875,
  10096720.3125,
  10098984.375,
  10099103.125,
  10101800.0,
  10102284.375,
  10168135.9375,
  10186507.8125,
  10190696.875,
  10233775.0,
  10233779.6875,
  10253465.625,
  10301259.375,
  10335398.4375,
  10376289.0625,
  10379604.6875,
  10380329.6875,
  10381626.5625,
  10382620.3125,
  10382701.5625,
  10382726.5625,
  10383701.5625,
  10426812.5,
  10456092.1875,
  10483621.875,
  10491740.625,
  10491762.5,
  10492215.625,
  10492378.125,
  10495167.1875,
  10496246.875,
  10497623.4375,
  10498098.4375,
  10499415.625,
  10505012.5,
  10507031.25,
  10508006.25,
  10508201.5625,
  10508820.3125,
  10509150.0,
  10509257.8125,
  10509415.625,
  10509523.4375,
  10509535.9375,
  10509587.5,
  10509787.5,
  10509801.5625,
  10509968.75,
  10509970.3125,
  10510062.5,
  10510109.375,
  10510214.0625,
  10510256.25,
  10510417.1875,
  10510421.875,
  10510587.5,
  10510615.625,
  10510857.8125,
  10510912.5,
  10511028.125,
  10511320.3125,
  10511385.9375,
  10511423.4375,
  10511593.75,
  10512109.375,
  10512242.1875,
  10513353.125,
  10513729.6875,
  10515221.875,
  10520542.1875,
  10532196.875,
  10532740.625,
  10573425.0,
  10628867.1875,
  10631404.6875,
  10649678.125,
  10710756.25,
  10715351.5625,
  10726079.6875,
  10730131.25,
  10734125.0,
  10734398.4375,
  10741070.3125,
  10743126.5625,
  10743710.9375,
  10747173.4375,
  10758114.0625,
  10760028.125,
  10769201.5625,
  10807843.75,
  10814521.875,
  10821571.875,
  10828520.3125,
  10832264.0625,
  10856915.625,
  10864357.8125,
  10918748.4375,
  10918775.0,
  10954748.4375,
  10961151.5625,
  10965453.125,
  10966560.9375,
  10970059.375,
  10970348.4375,
  10970737.5,
  10971503.125,
  10972400.0,
  10972715.625,
  10972768.75,
  10972875.0,
  10974218.75,
  10975379.6875,
  10976545.3125,
  10977331.25,
  10980662.5,
  10984535.9375,
  10988976.5625,
  10990892.1875,
  10991371.875,
  10991457.8125,
  10991687.5,
  10991867.1875,
  10992367.1875,
  10992642.1875,
  10992817.1875,
  10992898.4375,
  10992950.0,
  10993178.125,
  10993254.6875,
  10993287.5,
  10993353.125,
  10993793.75,
  10993875.0,
  10994062.5,
  10994084.375,
  10994106.25,
  10994137.5,
  10994182.8125,
  10994262.5,
  10994292.1875,
  10994300.0,
  10994303.125,
  10994307.8125,
  10994329.6875,
  10994479.6875,
  10994479.6875,
  10994523.4375,
  10994528.125,
  10994540.625,
  10994584.375,
  10994607.8125,
  10994617.1875,
  10994671.875,
  10994725.0,
  10994729.6875,
  10994820.3125,
  10994821.875,
  10994828.125,
  10994843.75,
  10994854.6875,
  10994867.1875,
  10994925.0,
  10994939.0625,
  10994954.6875,
  10995060.9375,
  10995089.0625,
  10995154.6875,
  10995164.0625,
  10995184.375,
  10995184.375,
  10995187.5,
  10995190.625,
  10995204.6875,
  10995214.0625,
  10995285.9375,
  10995293.75,
  10995334.375,
  10995342.1875,
  10995360.9375,
  10995381.25,
  10995381.25,
  10995385.9375,
  10995385.9375,
  10995392.1875,
  10995418.75,
  10995421.875,
  10995421.875,
  10995437.5,
  10995440.625,
  10995451.5625,
  10995468.75,
  10995493.75,
  10995493.75,
  10995506.25,
  10995528.125,
  10995539.0625,
  10995643.75,
  10995671.875,
  10995689.0625,
  10995709.375,
  10995731.25,
  10995734.375,
  10995745.3125,
  10995746.875,
  10995793.75,
  10995815.625,
  10995818.75,
  10995832.8125,
  10995864.0625,
  10995907.8125,
  10995912.5,
  10995957.8125,
  10995992.1875,
  10996059.375,
  10996064.0625,
  10996065.625,
  10996078.125,
  10996081.25,
  10996103.125,
  10996109.375,
  10996242.1875,
  10996250.0,
  10996276.5625,
  10996434.375,
  10996500.0,
  10996567.1875,
  10996648.4375,
  10996657.8125,
  10996684.375,
  10996748.4375,
  10996760.9375,
  10996768.75,
  10996817.1875,
  10997046.875,
  10997157.8125,
  10997289.0625,
  10997295.3125,
  10997431.25,
  10997637.5,
  10997876.5625,
  10998407.8125,
  10998465.625,
  10998496.875,
  10998604.6875,
  10998659.375,
  10998684.375,
  10999493.75,
  10999617.1875,
  11000021.875,
  11000035.9375,
  11000095.3125,
  11000306.25,
  11000457.8125,
  11000640.625,
  11000695.3125,
  11001281.25,
  11001317.1875,
  11001328.125,
  11001356.25,
  11001365.625,
  11001407.8125,
  11001440.625,
  11001450.0,
  11001453.125,
  11001468.75,
  11001473.4375,
  11001492.1875,
  11001523.4375,
  11001534.375,
  11001579.6875,
  11001651.5625,
  11001657.8125,
  11001762.5,
  11001943.75,
  11001946.875,
  11001962.5,
  11001987.5,
  11002057.8125,
  11002075.0,
  11002089.0625,
  11002118.75,
  11002140.625,
  11002153.125,
  11002257.8125,
  11002426.5625,
  11002478.125,
  11002529.6875,
  11002589.0625,
  11002701.5625,
  11002831.25,
  11002884.375,
  11003018.75,
  11003201.5625,
  11003214.0625,
  11003342.1875,
  11003464.0625,
  11003537.5,
  11003557.8125,
  11003579.6875,
  11003767.1875,
  11003779.6875,
  11003921.875,
  11003951.5625,
  11003967.1875,
  11004085.9375,
  11004234.375,
  11004459.375,
  11004590.625,
  11004685.9375,
  11004693.75,
  11004867.1875,
  11005126.5625,
  11005179.6875,
  11005812.5,
  11006681.25,
  11006735.9375,
  11006989.0625,
  11007001.5625,
  11007140.625,
  11007159.375,
  11007373.4375,
  11007415.625,
  11007446.875,
  11007482.8125,
  11007506.25,
  11007707.8125,
  11007728.125,
  11007759.375,
  11007814.0625,
  11007818.75,
  11007857.8125,
  11007890.625,
  11007918.75,
  11007929.6875,
  11007942.1875,
  11007946.875,
  11007959.375,
  11007959.375,
  11007971.875,
  11007998.4375,
  11008078.125,
  11008095.3125,
  11008101.5625,
  11008154.6875,
  11008204.6875,
  11008214.0625,
  11008257.8125,
  11008290.625,
  11008360.9375,
  11008364.0625,
  11008420.3125,
  11008464.0625,
  11008528.125,
  11008581.25,
  11008656.25,
  11008714.0625,
  11008756.25,
  11008767.1875,
  11008804.6875,
  11008807.8125,
  11008826.5625,
  11008834.375,
  11008862.5,
  11008948.4375,
  11009015.625,
  11009023.4375,
  11009110.9375,
  11009171.875,
  11009203.125,
  11009260.9375,
  11009473.4375,
  11009476.5625,
  11012743.75,
  11014612.5,
  11015421.875,
  11015943.75,
  11015976.5625,
  11016662.5,
  11016925.0,
  11016929.6875,
  11017753.125,
  11018825.0,
  11020676.5625,
  11027318.75,
  11027320.3125,
  11028362.5,
  11028370.3125,
  11037832.8125,
  11038293.75,
  11046914.0625,
  11051267.1875,
  11051351.5625,
  11051782.8125,
  11052800.0,
  11052903.125,
  11055126.5625,
  11055343.75,
  11055925.0,
  11056851.5625,
  11057420.3125,
  11058104.6875,
  11058665.625,
  11060215.625,
  11062082.8125,
  11062314.0625,
  11063559.375,
  11070346.875,
  11070935.9375,
  11070993.75,
  11072857.8125,
  11074231.25,
  11081189.0625,
  11082985.9375,
  11083876.5625,
  151831164.0625],
 [14.522412378749344,
  7.379992778300393,
  20.64045855497977,
  5.332989456358101,
  6.602010812981603,
  9.893424004460336,
  5.589378615973168,
  110.90538478687675,
  11.622872544491031,
  239.58983928274074,
  38.329248057040125,
  30.156670690525814,
  18.312726415391275,
  13.213046940069017,
  55.499476705499426,
  12.635008758636628,
  11.823825567831367,
  30.870911107769977,
  18.281691392176455,
  10.706446490391627,
  31.559052674672955,
  15.895855634874952,
  98.36152771154717,
  12.037312005452641,
  29.274761241454954,
  13.630214704948756,
  21.219575188132197,
  6.892693108972041,
  21.781125514213358,
  13.931343136319423,
  11.954458275133598,
  66.73275424795852,
  5.568376154226905,
  23.72959956898433,
  41.37533968905696,
  5.931186402321224,
  12.167238277426103,
  25.025206949742767,
  17.938851569538727,
  24.966618836864853,
  6.212778421386231,
  11.368993598320168,
  6.399029897659369,
  41.36418256274479,
  36.82693051703742,
  144.3227768443324,
  16.808172091539692,
  5.194234037820371,
  15.826552768804975,
  57.989383274771626,
  8.818435944969853,
  65.68688413196178,
  44.64478270348457,
  26.037015799208202,
  11.706944921685144,
  49.79007227663506,
  6.111296736451919,
  19.975285437579615,
  38.627272412796614,
  20.775486773095025,
  62.67762721634249,
  101.35731056713595,
  34.74045931880837,
  81.32856589187013,
  8.043412303762292,
  33.54386198889696,
  51.51636391299296,
  13.607063649881042,
  15.14491516945016,
  72.2685194215291,
  38.248752075632794,
  20.365521491486042,
  22.43134182092655,
  16.904271993482574,
  5.138220103948399,
  6.277932633580717,
  53.86370259348103,
  30.056699337048926,
  45.11949815537223,
  53.09602526003558,
  7.824758730285324,
  25.136474572811654,
  56.965452627250045,
  8.687201982735573,
  6.938473423169307,
  5.061121097064851,
  133.32872733439592,
  17.765340732483008,
  5.4015125207279135,
  14.192864833554385,
  14.53935567800184,
  34.63345622001923,
  60.79644814664843,
  50.80526275469903,
  218.1899287082013,
  56.883072513345894,
  79.25660573943448,
  40.87695836720245,
  6.615884398178277,
  81.50185818548036,
  102.97346715112994,
  20.027907061482868,
  8.699758766044383,
  8.653194328714939,
  5.11359448265291,
  32.62721118079563,
  35.0449669205988,
  55.33110389141369,
  5.110818410734519,
  93.43291721065732,
  52.821372728578766,
  23.473756405170686,
  77.05490659116731,
  20.26188446311911,
  5.183666374246042,
  8.414982115015128,
  20.462390948402188,
  5.129229862486688,
  48.21827717754643,
  40.19935426473932,
  43.857260625794666,
  42.096901830625924,
  106.04273309945751,
  21.24207218911864,
  18.085464898067954,
  36.84685085045022,
  71.4040101997111,
  6.0098415035890405,
  21.322507728958445,
  70.88419239166693,
  21.676578034426587,
  5.943412747191946,
  35.61649819708623,
  38.9113423970728,
  91.44714429286708,
  25.58894565757637,
  22.828962666450703,
  14.075128963828705,
  45.87313652808515,
  28.78495013587629,
  41.77044657502613,
  19.618697033701633,
  88.97796121291506,
  21.760966990165194,
  11.21048337112779,
  15.358939869131913,
  12.189919634745847,
  78.24750960999646,
  26.712369809335325,
  15.944981423681638,
  242.83830116691715,
  7.529317333603413,
  105.02939197105052,
  133.92301032128358,
  25.295187320693053,
  23.365490163923603,
  43.029147650747824,
  131.39301469355243,
  97.31593549434014,
  87.521590483323,
  61.39647622025864,
  20.490722849759596,
  9.604603074561172,
  60.61784056738845,
  18.211658646080693,
  20.698077379687824,
  6.137167577835072,
  49.77853747039155,
  10.207888441844103,
  11.246917884478261,
  47.25324192104496,
  65.1293254388128,
  24.19459803784773,
  90.26254115307393,
  6.368356953115741,
  41.99809989349299,
  69.9163799693189,
  20.103812425098596,
  37.73681860791247,
  35.39716302438328,
  67.22827888123679,
  8.150677163539159,
  15.261292399905834,
  35.62851392946617,
  37.95106830557609,
  88.34180028037551,
  38.56939648758064,
  18.58338911410286,
  34.78728885198922,
  71.95705252916208,
  12.201135906447666,
  31.560836051382026,
  65.68865727732911,
  14.157884227898197,
  29.045566711075583,
  78.57121345446706,
  16.51805070585848,
  37.07751422170219,
  39.811363934069014,
  66.91778807051035,
  39.97472551380958,
  113.47137292174581,
  54.68893151273588,
  73.52122735313056,
  143.92317869111557,
  32.30473596794936,
  20.765415229711678,
  17.19495925591494,
  34.37018966297808,
  5.17054730848159,
  31.298277229468493,
  8.813936063854511,
  90.3984176820932,
  53.923879350752244,
  30.9926977350519,
  42.35809036129966,
  92.02353793919525,
  12.943498850171082,
  8.49106700204624,
  16.261490936825016,
  23.145593562653506,
  7.103240281158632,
  33.597323951427995,
  47.61768654145746,
  25.227027592436468,
  11.770444260630278,
  18.887715137995105,
  29.1570302577205,
  8.975879372478412,
  14.633814264315477,
  57.21331724648192,
  71.59894134821812,
  9.524298882366455,
  42.37782901406701,
  18.515633035913623,
  30.293176114957753,
  14.024719768557528,
  51.103840693512296,
  32.55384439271148,
  5.129229784348003,
  22.610802466056544,
  58.84651464631091,
  37.923204299143606,
  51.98735238355061,
  48.59291253744345,
  31.419781165825405,
  45.019563310617,
  54.23399907653448,
  50.17771796853461,
  8.379946136630355,
  29.252189393661375,
  24.86740550440888,
  15.949126032680825,
  97.33661159595455,
  15.561169269121493,
  12.384890255553714,
  13.8193893929785,
  49.83061648020405,
  32.56327712838208,
  58.8399828527068,
  237.85023961168335,
  7.7736183744922736,
  73.43005564031378,
  10.87298271981529,
  32.91130203303636,
  20.128954653718885,
  57.489619366928174,
  18.654524622780475,
  57.803252877912925,
  34.52633665409199,
  10.908016378593254,
  90.55986605844916,
  39.87954447189963,
  7.361030551778287,
  44.30483807294905,
  11.604315674162388,
  9.778760016575207,
  41.88627705406574,
  141.11519655060658,
  42.29486581924668,
  11.205805691542148,
  51.04132662734624,
  72.65724673866539,
  94.95910394293396,
  128.22150229315847,
  42.96803135892449,
  39.70589819102901,
  62.97577086631011,
  56.040448701044035,
  26.89076669720656,
  28.70643271042061,
  23.057130051898277,
  24.656002302565632,
  175.58727149813038,
  96.56489921977861,
  73.81481916990167,
  16.73537135017913,
  16.32927302516183,
  31.48858583932537,
  10.071401486318006,
  9.215875831464261,
  108.5710705269465,
  5.2549536851257725,
  32.248156290893974,
  33.957869202641014,
  32.964145367463125,
  40.573141224991396,
  10.213821364515116,
  75.20565362000208,
  20.199526191126964,
  8.340869583580611,
  11.707933262941273,
  47.18005474009257,
  22.87487543238023,
  70.92953956136007,
  13.495330328738577,
  12.48425162350785,
  34.22227895519063,
  102.38712895673116,
  43.83356990840067,
  88.34045937630535,
  25.395164809745893,
  18.7900577058126,
  74.24061031034293,
  30.916743372846877,
  8.226542878960796,
  5.598321495504657,
  26.368431469700823,
  105.39909001328233,
  64.54209124714222,
  13.068654335637554,
  64.25442580482877,
  11.154716673104005,
  5.092807836105357,
  6.484505066144988,
  49.718081750008245,
  65.17715470391873,
  47.25530126913611,
  43.63685743542191,
  7.045660350381905,
  43.07968407081705,
  5.532713800152551,
  86.55904776009969,
  35.797390368598755,
  11.291655863136437,
  33.31790632291547,
  6.580314312994005,
  69.01943688051287,
  12.552490670146248,
  21.30785772606113,
  14.446971888551234,
  22.60316753586698,
  5.727697313815782,
  51.09028267891261,
  22.94967542553259,
  87.79325523340307,
  6.14436683665391,
  16.856345530811385,
  5.413100848723027,
  56.06986456302683,
  63.13061269126626,
  44.49310387576003,
  50.28868702987505,
  59.08755163070293,
  10.63697377265278,
  7.139064264230046,
  45.19612244852938,
  50.52060738802095,
  72.70241379232286,
  6.41688728005241,
  73.87571274690183,
  83.42023447666307,
  79.86328094613414,
  24.26078634954178,
  11.271004173151912,
  6.6723420695943965,
  43.81958941820018,
  76.5907632961646,
  10.742845244445874,
  67.45521577952377,
  47.74656676021197,
  24.90214744210141,
  89.79048978290905,
  73.29307364989036,
  41.27020773928539,
  91.32238095999966,
  54.006891372899105,
  45.166022761317116,
  53.72552244495233,
  19.00416897676892,
  10.234956892287236,
  31.150256449074554,
  11.375835129324832,
  9.362030025921142,
  12.54604748665758,
  5.93530311637763,
  66.94272458132205,
  43.20429333193299,
  40.26482319304364,
  46.15367575597653,
  50.3842935491331,
  50.538533223905105,
  15.141073759084229,
  82.27864007528866,
  79.40186545318585,
  40.61956623732796,
  75.02755334530434,
  61.61021772617565,
  53.00910946636314,
  12.266750407379732,
  40.28557965374314,
  53.02860554010428,
  79.78618347284286,
  6.055262311275415,
  5.442949892378668,
  7.605455311997575,
  18.206643181895934,
  74.89510468217578,
  6.143096656032979,
  24.005593130002406,
  90.69291186560322,
  10.18503435706545,
  9.124724430604719,
  26.91195178047335,
  7.789713694249968,
  47.305292758201816,
  5.578261085101192,
  23.779627651714872,
  140.5131643852066,
  57.684767173655224,
  71.74922248413849,
  33.15557240230683,
  16.645703464523116,
  37.09978429089474,
  7.894001915171967,
  57.65397880195902,
  17.908313813093184,
  23.728435632052374,
  46.471619302091355,
  18.873742977152524,
  7.402586579082309,
  79.60676749710838,
  15.644307866849804,
  5.104546910250049,
  46.156325523124764,
  54.78105417115934,
  24.350339269192034,
  11.270600071092053,
  5.970962042509847,
  29.729815414153855,
  55.42787021239333,
  11.511741688248748,
  22.317023808498988,
  6.353873037019125,
  47.78721957131938,
  71.83927135959942,
  52.812010606745616,
  11.73396454901019,
  6.128699943164398,
  72.72429150818456,
  6.2695976683471235,
  21.574502587042417,
  58.62398073711046,
  8.91774185806324,
  10.527213957027836,
  35.63637968256597,
  13.006443075551829,
  9.973367310432584,
  33.96371940830099,
  9.014396067181705,
  28.32770657305099,
  26.724477154180327,
  27.413789320210032,
  5.51665082833289,
  25.807665048402825,
  20.628770656301246,
  41.40895835186115,
  33.669607133354376,
  18.839680007110477,
  10.146518890411087,
  100.24903500901154,
  8.97526399302785,
  46.12767367685594,
  44.8528162546497,
  9.678039386429303,
  74.69354912229892,
  23.688487664224038,
  17.783860431202612,
  48.15473603194975,
  112.6960902922857,
  42.28987689297926,
  37.83654284102215,
  47.870869933050024,
  51.77640768792982,
  20.551599721648564,
  25.732533245162077,
  22.586433402874867,
  24.771827004221464,
  27.0755747299956,
  132.04403199094142,
  10.831232520090527,
  17.74755825429985,
  59.12946582589614,
  28.91359278764569,
  7.774654378966613,
  7.51386631583635,
  43.02177169499613,
  23.293211590279945,
  56.976170899326235,
  123.80055689238287,
  46.023159613972936,
  37.01139820308629,
  22.55573510783127,
  39.56422824316637,
  50.54096740408959,
  29.83191350763388,
  7.760900921146966,
  10.161441222067927,
  20.659340900155392,
  42.49410864870074,
  61.48666526072617,
  6.533091924921097,
  19.567501002232024,
  14.841877657825389,
  57.217800581628865,
  36.16787675286571,
  14.163766541589036,
  70.89983490439415,
  5.845871744444563,
  8.006113356538455,
  62.8551049545252,
  55.018646640655724,
  40.15736045756692,
  51.06734032581019,
  15.947039948299961,
  8.747152424439077,
  20.40536033958457,
  5.968212555372477,
  88.80924461630083,
  7.638357147601251,
  28.259044901676248,
  18.56045379661847,
  58.95708731326854,
  46.13927630527151,
  75.56172059662609,
  35.3242777282921,
  20.570639023681895,
  20.7570338432579,
  7.14513306217448,
  6.459396846820025,
  49.984578783008516,
  18.044423213754854,
  14.291447011305111,
  81.06273308354982,
  33.29065497448798,
  12.846827789713984,
  69.14397568365047,
  12.16431403940248,
  11.578217509744455,
  12.299540943370044,
  15.082951351340133,
  5.8731566575797,
  71.67585563358021,
  12.079145396137118,
  5.82930794838644,
  8.357428911293017,
  45.287975315295256,
  14.34790704387889,
  9.626000493755786,
  65.64604644327783,
  70.29590801542925,
  8.488143072221453,
  18.414947886021597,
  52.61199662979216,
  19.70833541476396,
  23.972217143270708,
  50.165246691928864,
  11.696039049206354,
  9.90131697830581,
  109.82602011050268,
  52.06545966448377,
  8.866436019682371,
  17.827981300158722,
  66.59581942026331,
  14.983261828746626,
  14.567153196507343,
  5.552611535652301,
  7.056759238299655,
  99.52177236852084,
  5.194565927259784,
  61.09886400546656,
  61.431918301097504,
  37.2401508968075,
  30.250622325670115,
  40.12151363885099,
  9.592989421076886,
  5.2557981885610765,
  14.18399562882538,
  50.59566804869893,
  18.625192550667332,
  48.490680927867345,
  22.87134627150619,
  109.37044623412048,
  106.34108044245549,
  6.464987511644463,
  14.570446903717254,
  66.03175718868997,
  6.387340647990015,
  51.078559952495766,
  20.18824715923365,
  8.222716356372963,
  16.52572507429681,
  90.60396527912548,
  26.42656241481182,
  58.496289283191125,
  8.504836798964465,
  17.214170734135003,
  21.910773970132258,
  8.066199133290901,
  16.426475372947156,
  31.55928003978359,
  9.643594835024379,
  33.109470026900816,
  47.76880120063448,
  50.35714206248249,
  7.822592524870712,
  11.65826066479001,
  66.16738335509929,
  17.6019875375003,
  95.67159483511644,
  18.740950540638618,
  67.93370128213776,
  10.591953068950675,
  39.55022381204277,
  16.680459756820152,
  55.75997547727652,
  12.082275523474362,
  16.024726709903177,
  13.41461577722364,
  66.3707870939449,
  72.02886249385459,
  16.255262604737894,
  12.137290851081648,
  8.96625395430296,
  5.0437588969339755,
  16.08265457403393,
  13.951048604127413,
  11.973080085858133,
  56.65635095026821,
  5.975008985054715,
  66.00135414520204,
  70.58633068872265,
  68.45497708632652,
  38.88817008071914,
  20.903892418355476,
  73.420858115235,
  27.77954518009039,
  13.20145840645084,
  10.306507253234834,
  5.353568382699863,
  7.116308624888159,
  30.497547025942197,
  14.986564311815544,
  18.141442000641227,
  5.365899879173057,
  17.799242445545108,
  18.197223218896994,
  42.14215427227429])

Graph for home page

Detected energies during the discharge + Energy spectrum

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