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

import matplotlib.cm as cm
#from matplotlib.mlab import griddata

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" t3pa2cls_XV_pc - zkousim vyzobat vysoke energie (jednotlive interakce) - funkce single_interaction. Dale delam prumernou velikost stopy stopy interakce pro danou energii - funkce size_of_interactions_average. Dale delam spektra 2 casti vyboje podle zadaneho casu - primarne pro double breakdown, tj. funkce energy_spectra_doublebreakdown

#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 = 44395  #test discharge for which the notebook will definitely work
shot_no = 45362
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
                calib[j][i]=float(word[j]) #j = radek, i = sloupec0 - pouze pro stavajici kalibraci - verze XV
    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))
        ClmNo.append(int(int(matrix_index[i]))//int(256)) #ver XV - bude pro novou kalibraci - nyn9 to bere 3patnou kalibraci jednotlivych pixelu (cislovani radku a sloupcu stejne jako v pixetu)
        RowNo.append(int(int(matrix_index[i]))%int(256)) ##ver XV - bude super pro novou kalibraci - odpovida radkum a sloupcum v pixetu. Nyni bere spatne kalibrace pixelu (viz fce load_calib)
    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 Timewalk_parameters_Si():
    #SI - korekce na TimeWalk - parametry
    A = -1.21988
    B = 4.33638
    C = 29.5075
    D = 1
    sigma_A=0.7013
    sigma_B=0.1366
    sigma_C=4.753
    sigma_D=0
    return A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D

def Timewalk(E,T):
    i=0
    A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
    while i < len(T):
        timewalk=float(0)
        timewalk = (C / (E[i] - B)) + A 
        if(timewalk<0):
            timewalk=0
        T[i]=T[i]-timewalk
        i=i+1
        
    return E,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=[]
        N_tot=[]
        eventNo=0 #for high_energy_event function
        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))
            N_tot.append(len(TCl)) #new in ver. XV
            
            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') 
            #eventNo=single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path)
            i=i+1
    t3pa_cls.close
    return T_first, E_tot, N_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, drawstyle='steps-mid')
    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.plot(T_hit, count, drawstyle='steps-mid')
    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.fill_between(xmean, count,step='mid') #pre mid post
    plt.plot(xmean, count, drawstyle='steps-mid')
    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, drawstyle='steps-mid')
    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].fill_between(x2, y2,step='mid') #pre mid post
    ax[1].plot(x2, y2, drawstyle='steps-mid')
    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

def single_interaction(eventNo,Etot, Tfirst, ECl, TCl, RowNoCl, ClmNoCl,path): # ver XV
    #tuto funkci volam ve funkci file_t3pa_cls_new(file_t3pa_cls,T):
    #This function is calling in function file_t3pa_cls_new(file_t3pa_cls,T)
    
    high_E = 1500 #[keV] - zadani, jaka energie je pro me vysoka

    if(Etot>=high_E or len(ECl)>15):
    #if (len(ECl)>15): #len(ECl) = pocet pixelu dane interakce
        print("Total energy of track: "+str(Etot)+" keV\n")
        A, B, C, D, sigma_A, sigma_B, sigma_C, sigma_D = Timewalk_parameters_Si()
        
        ######################
        #SI - rekonstrukce z(t) - parametry
        a=0.0193947
        b=-0.0273934
        sigma_a=0.0004163
        sigma_b=0.009832
        #######################
        
        time=[]
        sigma_time=[]
        time_min = []
        time_max=[]
        z_rec=[]
        z_rec_min=[]
        z_rec_max=[]
        sigma_z_rec=[]
        
        for i in range(0, len(TCl)):
            time.append(TCl[i] - Tfirst)
            #chyba urceni casu = chyba timewalk ->
            sigma_time.append(np.sqrt(((sigma_C / (pow((ECl[i] - B), D)))*(sigma_C / (pow((ECl[i] - B), D)))) + (((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))*((sigma_B*C*D) / (pow((ECl[i] - B), (D + 1))))) + (((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))* ((sigma_D*C*math.log(ECl[i] - B)) / (pow((ECl[i] - B), D)))) + (sigma_A*sigma_A)))
            time_min.append(time[i] - sigma_time[i])
            time_max.append(time[i] + sigma_time[i])

            #//
            #//rekonstrukce hloubky interakce
            z_rec.append(a * time[i] + (b))
            sigma_z_rec.append(np.sqrt(((time[i]*time[i])*(sigma_a*sigma_a)) + (sigma_b*sigma_b) + (a*a)*(sigma_time[i]*sigma_time[i])))
            z_rec_min.append(z_rec[i] - sigma_z_rec[i])
            z_rec_max.append(z_rec[i] + sigma_z_rec[i])
            #//
            
        single_interaction_file=str(path)+'event'+str(eventNo)+'.txt'
        single_interaction_2D_fig=str(path)+'2D_event'+str(eventNo)+'.png'
        single_interaction_3D_fig=str(path)+'3D_event'+str(eventNo)+'.png'
        
        with open(single_interaction_file, "w", encoding="utf-8") as f:
            f.write('#\n')
            f.write('# x_mean   y_mean   z_rec   z_rec_min   z_rec_max   sigma_z_rec\n')
            f.write('#E_tot='+str(Etot)+'\n')
            f.write('\n')
            for m in range(0,len(TCl)):
                    #print(str(RowNoCl[m])+'     '+str(ClmNoCl[m])+'     '+str(z_rec[m])+'     '+str(z_rec_min[m])+'     '+str(z_rec_max[m])+'     '+str(sigma_z_rec[m])+'     '+str(ECl[m])+'\n')
                    f.write(str(RowNoCl[m])+'     '+str(ClmNoCl[m])+'     '+str(z_rec[m])+'     '+str(z_rec_min[m])+'     '+str(z_rec_max[m])+'     '+str(sigma_z_rec[m])+'     '+str(ECl[m])+'\n')

        f.close
        
        ############################### plot
        #binwidth=1
        if ((max(RowNoCl)-min(RowNoCl))>=(max(ClmNoCl)-min(ClmNoCl))):
            rang=max(RowNoCl)-min(RowNoCl)
            #bins=(range((min(RowNoCl)-1), (max(RowNoCl)+3), binwidth), range((min(ClmNoCl)-1), (min(ClmNoCl)+rang+3), binwidth))
            x_min=min(RowNoCl)-1
            x_max=max(RowNoCl)+2
            y_min=min(ClmNoCl)-1
            y_max=min(ClmNoCl)+rang+2
        else:
            rang = max(ClmNoCl)-min(ClmNoCl)
            #bins=(range((min(RowNoCl)-1), (min(RowNoCl)+rang+3), binwidth), range((min(ClmNoCl)-1), (max(ClmNoCl)+3), binwidth))
            x_min=min(RowNoCl)-1
            x_max=min(RowNoCl)+rang+2
            y_min=min(ClmNoCl)-1
            y_max=max(ClmNoCl)+2
            
        ###############$$$$$$$$$$$$$$$$$$$$
        rows, cols = (x_max-x_min, x_max-x_min)
        #arrECl = [[0]*cols]*rows
        arrECl= [[0] * cols for _ in range(rows)]
        for i in range (x_min,x_max):
            for j in range (y_min,y_max):
                for k in range(0,len(RowNoCl)):
                        if ((i==RowNoCl[k]) and (j==ClmNoCl[k])):
                            #print(str(i)+"   "+str(j)+"    "+str(RowNoCl[k])+"    "+str(ClmNoCl[k])+"\n")
                            arrECl[i-x_min][j-y_min]=ECl[k]
        
        arrECl=np.transpose(arrECl)
        arrECl[arrECl == 0] = np.nan
        
        #2D plot
        plt.imshow(arrECl,cmap='jet',  interpolation='none', origin=[(x_min-0.2),(y_min-0.2)], aspect='equal', extent=[x_min-0.2,x_max,y_min-0.2,y_max])
       
        ###############$$$$$$$$$$$$$$$$$$$$
        cb=plt.colorbar(pad=0.005)
        plt.tick_params(axis="x",direction="in")
        plt.tick_params(axis="y",direction="in")
        cb.ax.tick_params(axis='y', direction='in')
        cb.set_label('Energy in pixel [keV]')
        plt.xlabel('x [pixel]')
        plt.ylabel('y [pixel]')
        plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
        plt.savefig(single_interaction_2D_fig, dpi = 1000)
        plt.show()
        plt.close()
        ############################### plot
        
        
        ############## 3D plot
        # Creating figure
        fig = plt.figure(figsize = (16, 9))
        ax = plt.axes(projection ="3d")

        # Add x, y gridlines 
        ax.grid(b = True, color ='grey', # 
                linestyle ='-.', linewidth = 0.3, 
                alpha = 0.2) 

        # Creating plot
        sctt = ax.scatter3D(RowNoCl, ClmNoCl, z_rec, 
                            alpha = 0.8,
                            c = ECl, 
                            cmap = 'jet', 
                            s = 4000, #size of bubbles in plot
                            marker ='o')
        plt.title('$E_{tot}$ = '+str(round(Etot, 2))+' keV, $N_{tot}$ = '+str(len(ECl)))
        ax.set_xlabel('x [pixel]', fontweight ='bold') 
        ax.set_ylabel('y [pixel]', fontweight ='bold') 
        ax.set_zlabel('z [mm]', fontweight ='bold')
        cb=fig.colorbar(sctt, ax = ax, shrink = 0.5,aspect = 10,pad=0.00005) #
        cb.ax.tick_params(axis='y', direction='in')
        cb.set_label('Energy in pixel [keV]', fontweight ='bold')
        plt.savefig(single_interaction_3D_fig, dpi = 1000)
        # show plot
        plt.show()
        plt.close()
        ############## 3D plot
        
        eventNo=eventNo+1
    return eventNo

def size_of_interactions_average(Etot, Ntot,file_hist):  #ver XV
    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=[]
    count=[]
    nop=[] #number of pixels
    nop_mean=[] #mean number of pixels in dE
    sigma=[] #root mean square error of the mean - stredni kvadraticka chyba prumeru
    
    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.append(0)
        nop.append(0) 
        nop_mean.append(0) 
        sigma.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
            nop[E_index]=nop[E_index]+Ntot[i] #ver XV
    
    for l in range(0,pocet):
        if nop[l]!=0:
            nop_mean[l]=float(nop[l])/float(count[l])
    
    #####Computing sigma
    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])):
            if nop_mean[E_index]!=0:
                sigma[E_index]=sigma[E_index]+((float(Ntot[i])-float(nop_mean[E_index]))*(float(Ntot[i])-float(nop_mean[E_index])))
            
    for l in range(0,pocet):
        if nop[l]!=0 and count[l]>1:
            sigma[l]=np.sqrt(sigma[l]/(float(count[l])*(float(count[l])-1)))
        if nop[l]!=0 and count[l]==1:
            sigma[l]=10
    #####Computing sigma
    
    #####write to file
    with open(file_hist, "w", encoding="utf-8") as hist:
        hist.write('# x_mean [keV]     mean_size_of_interaction [pixel]     sigma_mean_size_of_interaction\n')
        for m in range(0,len(xmean)):
                hist.write(str(xmean[m])+'     '+str(nop_mean[m])+'     '+str(sigma[m])+'\n')
        
    hist.close
    #####write to file
    
    plt.rcParams.update({'font.size': 20})
    fig, ax = plt.subplots(figsize =(10, 7))
    #plt.fill_between(xmean, nop,step='mid',color='blue', alpha=0.3) #pre mid post
    #plt.fill_between(xmean, nop,step='mid') #pre mid post
    plt.plot(xmean, nop_mean, drawstyle='steps-mid')
    plt.errorbar(xmean, nop_mean, yerr = sigma,fmt ='.', linewidth=0.7)
    plt.title(detector+', #'+str(shot_no))
    plt.xlabel('Energy [keV]')
    plt.ylabel('Average size of interaction [pixel]')
    plt.xlim(0,500)
    plt.ylim(0,17)
    #ax.set_yscale('log') #log scale y
    caption = '# x = energy in keV, dE= '+str(dE)+' keV'
    plt.savefig(figure_size_interaction, dpi = 1000)
    return

def energy_spectra_doublebreakdown(Etot, Tfirst, file_hist):
    #this function separates discharge into two parts acording to time Tsep and creates energy spectrum for each part.
    #This is for discharges with double brakdown
    
    Tsep=14636000 # [ns] 
    
    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_first=[]
    count_second=[]
    for l in range(0,pocet):
        count_first.append(0)
        count_second.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])):
            if(Tfirst[i]<Tsep):
                count_first[E_index]=count_first[E_index]+1
            else:
                count_second[E_index]=count_second[E_index]+1
    
    caption = '# x = energy in keV, dE= '+str(dE)+' keV'
    with open(file_hist, "w", encoding="utf-8") as hist:
        hist.write('#\n')
        hist.write('#'+str(caption)+'\n')
        hist.write('# x_mean     count_before_2nd_breakdown     count_2nd_breakdown\n')
        hist.write('\n')
        for m in range(0,len(xmean)):
                hist.write(str(xmean[m])+'     '+str(count_first[m])+'     '+str(count_second[m])+'\n')
        
    hist.close
    
    plt.rcParams.update({'font.size': 20})
    fig, ax = plt.subplots(figsize =(10, 7))
    #plt.fill_between(xmean, count_first,step='mid') #pre mid post
    plt.plot(xmean, count_first, label='Before 2nd breakdown', drawstyle='steps-mid')
    plt.plot(xmean, count_second, label='2nd breakdown', drawstyle='steps-mid')
    plt.legend(loc='best')
    plt.title(detector+', #'+str(shot_no))
    plt.xlabel('Energy [keV]')
    plt.ylabel('Count')
    plt.xlim(0,500)
    ax.set_yscale('log') #log scale y
    plt.savefig(figure_energy_spectra_doublebreakdown, 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'
size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction.txt'
file_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown.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'
figure_size_interaction= 'H03-W0051_shot_'+str(shot)+'size_interaction'
figure_energy_spectra_doublebreakdown= 'H03-W0051_shot_'+str(shot)+'energy_spectra_doublebreakdown'

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)
E,T=Timewalk(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, Ntot = 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);

([4694004.6875,
  4808577.372469842,
  4829393.75,
  4832160.4729947215,
  4835189.0625,
  4835197.5757367695,
  4838463.587216016,
  4850960.866667724,
  4851387.279298142,
  4853970.350343791,
  4854046.716112161,
  4854624.717582297,
  4860087.5,
  4865018.75,
  4865732.184104946,
  4866298.8950518835,
  4866820.3125,
  4891287.5,
  4895410.9375,
  4897640.609005003,
  4901851.880021661,
  4906698.095913895,
  4907011.771203193,
  4907448.401623274,
  4910232.661428186,
  4912437.5,
  4916085.417639641,
  4918602.300552371,
  4922876.293468087,
  4928815.518279881,
  4932015.625,
  4943148.4375,
  4955174.565488515,
  4983501.271457321,
  4992959.375,
  4995539.0625,
  5015819.558059051,
  5016340.575258462,
  5023478.125,
  5032326.636951393,
  5034235.9375,
  5036161.641974806,
  5039305.537404446,
  5041470.280810332,
  5066207.655221228,
  5066903.849012415,
  5074100.0,
  5078420.201330655,
  5080255.2745319875,
  5083355.3471233165,
  5092399.822676634,
  5092486.9018794475,
  5092812.5,
  5093047.821994288,
  5093134.201920701,
  5093407.056890622,
  5094117.1875,
  5094623.488912772,
  5094734.073307234,
  5094740.3410823615,
  5095188.582553806,
  5097135.9375,
  5097265.40053432,
  5098068.331329957,
  5098865.437252875,
  5099009.16006795,
  5099365.420171658,
  5099410.750296834,
  5099601.5625,
  5099937.374837869,
  5100758.796193155,
  5100768.069609331,
  5101492.1875,
  5103541.296586633,
  5104068.611881942,
  5104094.034230023,
  5104752.890088723,
  5105220.570571521,
  5105323.4375,
  5105490.720660989,
  5105508.521919656,
  5105670.3125,
  5105671.654121327,
  5106156.772065578,
  5106372.944268811,
  5106595.3125,
  5107134.549287804,
  5107575.198086805,
  5107876.683707601,
  5107878.872365402,
  5107975.0,
  5108042.8833440915,
  5108192.2144353185,
  5108380.900189651,
  5108586.00717675,
  5109306.25,
  5109388.417477825,
  5109615.625,
  5109736.163102475,
  5110691.666187518,
  5111401.966758304,
  5111490.625,
  5112127.963118512,
  5112895.888247758,
  5112950.026901975,
  5113793.384749336,
  5114179.6875,
  5114364.852440559,
  5114462.5,
  5114473.4375,
  5115124.227090571,
  5115191.536472652,
  5115344.204021185,
  5115502.417075702,
  5115741.348361351,
  5116485.328295426,
  5116562.5,
  5116593.489144347,
  5116613.610862971,
  5116779.1328282235,
  5117247.056000437,
  5118132.679415019,
  5118261.128204859,
  5118478.937814005,
  5119128.125,
  5119217.1875,
  5119785.825764455,
  5124015.963600086,
  5124304.707913356,
  5124370.098579549,
  5124881.092059543,
  5124911.775081755,
  5126035.781837893,
  5127626.681530933,
  5129140.625,
  5129414.498736343,
  5129651.5625,
  5129696.547701512,
  5130052.235399146,
  5130768.75,
  5131029.6875,
  5131199.604959133,
  5131902.723537738,
  5131969.08454612,
  5132062.88811622,
  5134320.331568867,
  5134728.871460731,
  5136223.4375,
  5140659.941584441,
  5141100.326357179,
  5141117.990946873,
  5141498.206717823,
  5142395.3125,
  5142449.96567222,
  5142633.769782632,
  5142744.78360987,
  5142804.6875,
  5143046.584066203,
  5143387.347854304,
  5143638.822733251,
  5143854.660056851,
  5144356.25,
  5144444.642546348,
  5144589.0625,
  5144602.817701543,
  5144605.477248827,
  5144990.878646369,
  5145059.420718178,
  5145086.2114315815,
  5145370.3125,
  5145383.284397847,
  5145655.891391298,
  5145755.686767004,
  5145761.781111819,
  5145801.1513136625,
  5146231.099199403,
  5146361.458257177,
  5146466.736728932,
  5146586.37906749,
  5146621.82611822,
  5146706.984524036,
  5146741.718585128,
  5146804.4247723725,
  5146855.479355085,
  5146897.620963,
  5147001.877914878,
  5147096.737451783,
  5147191.88517687,
  5147254.955318017,
  5147313.762920726,
  5147448.499477593,
  5147700.235121299,
  5147839.200961352,
  5147848.448028363,
  5147914.7544843005,
  5148065.781773489,
  5148207.461630596,
  5148429.6875,
  5148530.879588512,
  5148617.980795272,
  5148625.470771401,
  5148792.743337879,
  5148895.988933192,
  5148960.9375,
  5149676.867155789,
  5149759.319703105,
  5151380.520745591,
  5151446.105416869,
  5152190.504844647,
  5152366.62584822,
  5152744.703942179,
  5153000.766416478,
  5153104.4182126755,
  5153393.75,
  5153457.600290114,
  5154430.814618393,
  5155607.8125,
  5155617.56739411,
  5155625.0,
  5156575.0,
  5158364.786376255,
  5158650.278738839,
  5158675.961852983,
  5160023.575736752,
  5160367.6540790135,
  5162783.560865514,
  5165412.055243464,
  5167618.54034967,
  5176725.0,
  5180611.517740609,
  5181646.733742297,
  5192275.0,
  5202403.125,
  5207536.565125276,
  5209086.084071497,
  5209101.659606263,
  5229930.358088509,
  5231028.125,
  5236405.181011105,
  5236820.9681162825,
  5238782.884217574,
  5238928.150866442,
  5241415.445744058,
  5248083.826437938,
  5249844.8623016765,
  5250024.358008581,
  5253009.422054304,
  5253775.0,
  5261548.4375,
  5265316.266870025,
  5267368.591254162,
  5267375.206075867,
  5268998.249328477,
  5269252.355843868,
  5269370.3125,
  5269401.5625,
  5269443.630010111,
  5269774.351209665,
  5270050.0,
  5271162.998748115,
  5272202.392863058,
  5272243.9340468235,
  5272246.875,
  5273335.172039438,
  5274944.39831536,
  5283174.041905419,
  5286038.1293557575,
  5290095.992354738,
  5296206.195201195,
  5296207.713871118,
  5296422.0488007,
  5297942.1875,
  5297959.086511491,
  5298113.7839112785,
  5298234.203882818,
  5299062.5,
  5300073.318449883,
  5301200.037387405,
  5301501.5625,
  5302742.1875,
  5302985.292294794,
  5320713.121855185,
  5322578.125,
  5328812.494700109,
  5331696.197855707,
  5332261.657305221,
  5337693.027363412,
  5337718.614399251,
  5339570.760605811,
  5355314.63303687,
  5365040.289862701,
  5376591.261865186,
  5386071.26815407,
  5390136.1325436365,
  5391610.9375,
  5405304.62692363,
  5406559.375,
  5416752.587103924,
  5435344.377887121,
  5438416.787787994,
  5442623.028988975,
  5453554.23859176,
  5456244.805857645,
  5465288.462945165,
  5469689.0625,
  5471504.0964952065,
  5485812.792433178,
  5493850.565478801,
  5512331.191068753,
  5517009.228468536,
  5521196.789962279,
  5522676.168382783,
  5524935.9375,
  5527273.9210504,
  5528385.9375,
  5529890.235479234,
  5530542.993582652,
  5531178.516369288,
  5539693.753242967,
  5547794.814748696,
  5561018.011655362,
  5563015.190728106,
  5565985.569669171,
  5575432.648347869,
  5575469.374932138,
  5580503.3122413205,
  5597721.07133487,
  5612261.95665256,
  5641439.011720579,
  5650640.475272585,
  5655115.520086922,
  5660510.809538817,
  5661272.090823038,
  5661665.469593941,
  5662285.9375,
  5662943.9259487055,
  5662973.009469303,
  5663434.404433025,
  5671745.3125,
  5675628.151393408,
  5681021.875,
  5681710.415618129,
  5684340.295648728,
  5686427.870681439,
  5687791.214183656,
  5690356.595456704,
  5691048.4375,
  5693648.587101451,
  5694278.544786245,
  5695262.5,
  5701011.746802646,
  5706006.801610432,
  5709016.010009828,
  5719288.640826281,
  5724500.0,
  5727643.60315394,
  5732699.495307683,
  5747499.749950118,
  5747500.300883567,
  5747596.62360939,
  5754865.231153525,
  5755421.875,
  5758662.332813297,
  5758682.76462499,
  5759554.6875,
  5760360.9375,
  5760698.5331816785,
  5771644.772551745,
  5772427.634172109,
  5792246.875,
  5795078.032519066,
  5795747.320811018,
  5803335.292934461,
  5804651.410439777,
  5815487.400264339,
  5817965.625,
  5818124.43470538,
  5818200.0,
  5831084.375,
  5838781.601757914,
  5840365.747804157,
  5842970.3125,
  5845291.933640746,
  5848140.625,
  5852334.3672083225,
  5859795.489483133,
  5860608.973634419,
  5865496.991196872,
  5866614.5174266165,
  5866670.6561953975,
  5869540.305823777,
  5869549.838206735,
  5869892.1875,
  5871568.093375797,
  5872000.0,
  5874910.819270179,
  5874958.745822031,
  5874968.916049484,
  5883385.765179415,
  5884901.5625,
  5890356.146796166,
  5900065.958948857,
  5900085.770721988,
  5926671.662790114,
  5947995.261029914,
  5974207.313325031,
  5981109.479358432,
  5982521.349250061,
  5987249.4706328,
  5987253.608432798,
  6037691.172186762,
  6043754.6875,
  6043765.625,
  6047271.875,
  6057065.625,
  6083463.233678017,
  6088084.90542871,
  6096451.047523742,
  6105148.960883869,
  6105162.7231516475,
  6107698.4375,
  6108743.321860101,
  6156955.0551039865,
  6185094.879078507,
  6192882.8125,
  6222450.0,
  6234778.164475562,
  6263332.321396833,
  6276848.125982724,
  6310924.397318744,
  6333128.125,
  6347029.6875,
  6361256.25,
  6361302.620720931,
  6363876.095776006,
  6383368.140021875,
  6596945.277879903,
  6792235.708640398,
  7180392.237749754,
  7190928.125,
  7190937.5,
  7249229.6875,
  7369176.5625,
  7435844.281191614,
  8259487.585124541,
  8295015.804396403,
  8566520.3125,
  8612320.3125,
  8638875.0,
  8638907.198572712,
  9010590.625,
  105461715.15321432,
  409286308.8955168],
 [44.42477307749312,
  49.96609244976557,
  40.87389347834934,
  6.416660493022236,
  53.99412772129534,
  18.51147953711033,
  49.29685290640293,
  8.24947188315964,
  14.759983092403104,
  15.08774115961008,
  5.9256870301143545,
  23.977959645296348,
  52.5315881816371,
  75.29230807725236,
  44.72466138105205,
  11.927069315553487,
  47.71099851526634,
  47.02324246901132,
  156.84240298813347,
  7.117253397298565,
  71.11684278789919,
  6.69717526522069,
  39.88360040306763,
  5.531341860525387,
  154.41072202614725,
  56.397085184607306,
  48.319044820227134,
  18.770221024742526,
  67.16164093873091,
  14.549767833479011,
  43.69399748449416,
  85.25761533493149,
  36.76558753728017,
  27.968108377929276,
  56.437519819866694,
  120.5752728293347,
  71.13826336146273,
  8.260447168588755,
  62.08046277269942,
  16.613998095810693,
  59.13776635557357,
  7.319827029866923,
  7.364350653734773,
  69.15973056167854,
  47.413872665968285,
  92.58409865493707,
  101.60660372620617,
  26.50497962053751,
  6.151349605940142,
  7.3063586266766585,
  163.89181835138777,
  13.065116925778804,
  32.33130019745596,
  34.455816493414886,
  9.189220641010868,
  8.764947767926815,
  79.43192316466265,
  72.62886529136904,
  13.904085510985677,
  5.915055840004793,
  83.18112934327156,
  45.81474067159499,
  24.76604611075996,
  105.53179415813676,
  8.189727717088006,
  53.3582493801643,
  14.214331793834546,
  44.23880686128197,
  89.385686159402,
  62.152328680400814,
  8.885663069648574,
  33.998705913139716,
  58.224245106142725,
  21.390524878587676,
  26.0650561336442,
  152.1060258518192,
  37.02639810704824,
  8.427933468919303,
  69.05314730332555,
  6.502192218445554,
  18.570853939544524,
  33.95059649285611,
  24.816909159076705,
  36.11464670441421,
  21.56088952275647,
  56.34380564508474,
  23.725841937296806,
  18.733550573794172,
  7.153682172343619,
  47.30632576217731,
  72.85410690034762,
  5.653735805873015,
  8.300872719440944,
  26.90854572086416,
  17.746468750421595,
  37.522337956694905,
  92.21999539938129,
  63.55272800046399,
  8.409595016818066,
  13.268052716350402,
  16.74428160044593,
  33.40799817094388,
  82.80651568416958,
  17.708567465364606,
  5.708472262549371,
  31.162506082363564,
  66.79079780214342,
  10.102447256865965,
  62.83779386096224,
  29.704808097780376,
  16.364790848191966,
  41.87933344021094,
  11.920181263548098,
  7.944622374560327,
  69.92587012299553,
  73.64431261133691,
  30.536156037380355,
  28.285278047233597,
  68.34973400802046,
  92.75471294896701,
  17.12784852125473,
  15.683252458747484,
  6.802186053313718,
  10.128335591211648,
  73.57093761112411,
  52.4857470051361,
  14.532068300701411,
  16.410912501860707,
  71.17734772356275,
  18.888498588073,
  14.878988275433365,
  125.842377086918,
  14.379632929063094,
  11.319000423430362,
  74.10033836997975,
  25.71754013625781,
  71.49050589577294,
  117.34423604334512,
  28.37038220070245,
  81.84163068087084,
  71.0766851747359,
  22.608173133296752,
  10.553272582885974,
  64.2517514590782,
  44.432395966905965,
  7.126473320351014,
  116.14751379342539,
  75.93760672031826,
  17.653267416994023,
  7.7257003141274,
  25.05854467101308,
  57.57002974650062,
  82.76895298427584,
  7.590868918783576,
  8.867214371107352,
  48.17216493073142,
  28.657355722719522,
  49.60671178415252,
  14.809124885747092,
  22.079503586088233,
  39.26427399234346,
  53.26045898789833,
  5.226773727087648,
  88.84364277310306,
  33.72349433589182,
  19.287053601784734,
  9.55549747610324,
  15.118676797729096,
  50.7703106398665,
  116.59436326066961,
  25.233652324344742,
  23.02989358557103,
  10.348364652401134,
  8.789489833793379,
  6.684102023871341,
  20.92533071832475,
  12.052528607734255,
  6.418676703138052,
  77.17715622439955,
  27.593306922475765,
  65.25675256492187,
  15.689895000357028,
  24.238813694852315,
  19.160359229874775,
  65.88559311981061,
  11.659312226482983,
  16.85494050797344,
  10.68589801610649,
  25.176541697049075,
  39.31630482109196,
  15.18312088131779,
  15.92040167073918,
  25.07080375845274,
  31.50864675235237,
  31.01999465649287,
  6.795214057494854,
  23.121999092390322,
  57.24715739192866,
  22.891154531330407,
  20.16403729787228,
  17.101300505381584,
  25.019278782390913,
  6.906507750933781,
  31.02256975274764,
  9.603014389412401,
  105.66071307436671,
  16.137206728772213,
  25.435562286269658,
  63.16629552267574,
  13.16030929712963,
  6.9701697769207165,
  122.67411940072881,
  8.149128680072858,
  64.14872004564745,
  45.129760620554784,
  12.499743086677732,
  86.45970410041727,
  11.77839891782654,
  60.19460400417759,
  42.03534379556447,
  12.48536356238823,
  6.8204604263486255,
  13.058598253132624,
  79.77274773816802,
  7.249751160571848,
  31.164652004782994,
  5.901587208411874,
  14.198412381846978,
  54.51179804631836,
  34.81278273675066,
  8.213264264961525,
  31.652169043108945,
  57.45057773315235,
  26.027299144740525,
  6.793125521662823,
  68.09653692684053,
  17.32280264686723,
  71.5788572363732,
  11.997941591951136,
  73.8548765044674,
  70.19616970751521,
  23.189001907835646,
  31.851842694091072,
  42.82807823049187,
  29.459411339099248,
  28.087197356946508,
  11.20206094750102,
  70.10433247809218,
  69.19580145186355,
  40.1702134012693,
  21.838704273804154,
  85.96499394424217,
  78.67133150388811,
  8.755962162840033,
  82.28764158208897,
  47.599723993880765,
  14.503071748369363,
  10.24536037278826,
  39.29326322839205,
  17.25768320880104,
  19.452795927391264,
  7.671188874461247,
  68.84857946513178,
  16.272385048053163,
  6.593760873153204,
  17.88451791532223,
  57.74382893486751,
  18.3706941299025,
  51.76665355504497,
  6.254817225170997,
  11.410688082498234,
  32.73657344433986,
  22.736245580257837,
  19.246571118349145,
  8.198114181846222,
  31.53422323672576,
  100.59066231016293,
  7.131314472582982,
  61.88055534954937,
  72.9641875114753,
  6.1889834937319455,
  6.894347710505344,
  32.56282622804951,
  96.56951937424725,
  27.05964304737305,
  47.44217314407133,
  6.180020700919893,
  5.337296513239306,
  37.79389966033233,
  9.86539767533493,
  7.486360838910614,
  5.364228849553785,
  47.27854356067586,
  15.74096597368917,
  30.656465737352296,
  63.905701320320304,
  59.34748791447035,
  111.72761188401005,
  62.60250734860647,
  22.55547462399196,
  10.544053604547504,
  35.68140491290819,
  80.53030142649618,
  49.82146271426399,
  66.0937694663576,
  6.650986946553472,
  86.53735264185588,
  6.14894666873077,
  107.89280121863177,
  25.931266151148222,
  26.94891679642271,
  75.56408012183371,
  68.85333998443237,
  8.559985482282832,
  63.691029646788635,
  13.639161572916619,
  26.009151477429022,
  16.677395392992437,
  14.953884516854211,
  10.429658250098077,
  138.00425059820108,
  18.556338710559476,
  27.942588659077394,
  25.656333005160427,
  12.033347385784475,
  98.19493626563742,
  33.53547893988432,
  12.90442083450786,
  5.040689238908884,
  62.723149496761295,
  9.244241212944024,
  22.64771036369743,
  21.95062387581726,
  25.294053364567645,
  40.4874824803596,
  34.20575784242585,
  22.242386712972735,
  6.270924193246074,
  50.05542576202717,
  30.49320250073631,
  147.36930740630018,
  85.00285230322807,
  32.92653292115592,
  24.352340380366304,
  5.309640403255765,
  15.915199053699894,
  51.71643672254985,
  8.367305845709916,
  63.98367415685296,
  51.045566192689165,
  31.20501213035617,
  5.419988743630674,
  11.788013189525348,
  97.27917222395683,
  66.23909454626711,
  23.53559935309491,
  30.240776992380425,
  63.52051520997735,
  11.632998700676932,
  10.756258746089319,
  70.49394540294152,
  51.62258860312241,
  25.729364095207302,
  5.7037178737320335,
  97.78418139390564,
  103.2130468329833,
  6.502195507358726,
  21.103641194300113,
  81.99109528255512,
  29.88911322099194,
  90.40398765404383,
  30.426762543990836,
  8.840004125316065,
  5.697148352625648,
  23.220114972545407,
  49.03590769431175,
  66.45174410289283,
  67.75785875415326,
  117.73665211892482,
  16.644000842663353,
  72.39957495838564,
  85.960267272527,
  6.062092028385488,
  71.32858963515577,
  14.91189413095649,
  11.416089729073404,
  13.60455148167163,
  11.314320674529938,
  37.50809757008781,
  78.28151685532504,
  23.50884397843321,
  7.54574958135721,
  70.5350554084501,
  8.831730565031169,
  87.13604870248486,
  5.527369857055622,
  35.80491186809685,
  15.61458045052274,
  14.323009493522108,
  125.72841311383547,
  61.25612427969707,
  8.4714395189329,
  9.194254941847579,
  12.366277002578592,
  27.545958796040416,
  19.25230820209658,
  12.909903731668056,
  6.66289562686799,
  21.205064994367977,
  23.100781330576016,
  57.98689882701888,
  70.36522070258081,
  91.05176625136534,
  48.836232998085286,
  31.6006457136473,
  18.73940213124516,
  12.072092887537075,
  13.28521588157706,
  60.990045209094816,
  56.22617848382819,
  42.67962376222278,
  32.11178450382196,
  16.555939843393517,
  22.18400194651775,
  82.55157499855349,
  63.049798928917596,
  33.11012574999999,
  97.89674012493117,
  23.604730846246586,
  20.526507986082866,
  64.43449775942518,
  34.01014003296156,
  78.19815881172786,
  7.4304929158720965,
  30.83426765823061,
  6.969835923341012,
  27.857701414143335,
  33.982842434147976,
  13.954149107583603,
  49.079138589792635,
  49.62162157954233,
  91.27539384673928,
  29.697659926888676,
  5.078684851720425,
  8.332111499103203,
  185.54946407150334,
  74.85303069547231,
  49.010656213890314,
  86.84510455858383,
  20.427220179174753,
  73.23311408252631,
  188.93552586386815,
  10.452731280726248])

Energy spectrum

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

([4694004.6875,
  4808577.372469842,
  4829393.75,
  4832160.4729947215,
  4835189.0625,
  4835197.5757367695,
  4838463.587216016,
  4850960.866667724,
  4851387.279298142,
  4853970.350343791,
  4854046.716112161,
  4854624.717582297,
  4860087.5,
  4865018.75,
  4865732.184104946,
  4866298.8950518835,
  4866820.3125,
  4891287.5,
  4895410.9375,
  4897640.609005003,
  4901851.880021661,
  4906698.095913895,
  4907011.771203193,
  4907448.401623274,
  4910232.661428186,
  4912437.5,
  4916085.417639641,
  4918602.300552371,
  4922876.293468087,
  4928815.518279881,
  4932015.625,
  4943148.4375,
  4955174.565488515,
  4983501.271457321,
  4992959.375,
  4995539.0625,
  5015819.558059051,
  5016340.575258462,
  5023478.125,
  5032326.636951393,
  5034235.9375,
  5036161.641974806,
  5039305.537404446,
  5041470.280810332,
  5066207.655221228,
  5066903.849012415,
  5074100.0,
  5078420.201330655,
  5080255.2745319875,
  5083355.3471233165,
  5092399.822676634,
  5092486.9018794475,
  5092812.5,
  5093047.821994288,
  5093134.201920701,
  5093407.056890622,
  5094117.1875,
  5094623.488912772,
  5094734.073307234,
  5094740.3410823615,
  5095188.582553806,
  5097135.9375,
  5097265.40053432,
  5098068.331329957,
  5098865.437252875,
  5099009.16006795,
  5099365.420171658,
  5099410.750296834,
  5099601.5625,
  5099937.374837869,
  5100758.796193155,
  5100768.069609331,
  5101492.1875,
  5103541.296586633,
  5104068.611881942,
  5104094.034230023,
  5104752.890088723,
  5105220.570571521,
  5105323.4375,
  5105490.720660989,
  5105508.521919656,
  5105670.3125,
  5105671.654121327,
  5106156.772065578,
  5106372.944268811,
  5106595.3125,
  5107134.549287804,
  5107575.198086805,
  5107876.683707601,
  5107878.872365402,
  5107975.0,
  5108042.8833440915,
  5108192.2144353185,
  5108380.900189651,
  5108586.00717675,
  5109306.25,
  5109388.417477825,
  5109615.625,
  5109736.163102475,
  5110691.666187518,
  5111401.966758304,
  5111490.625,
  5112127.963118512,
  5112895.888247758,
  5112950.026901975,
  5113793.384749336,
  5114179.6875,
  5114364.852440559,
  5114462.5,
  5114473.4375,
  5115124.227090571,
  5115191.536472652,
  5115344.204021185,
  5115502.417075702,
  5115741.348361351,
  5116485.328295426,
  5116562.5,
  5116593.489144347,
  5116613.610862971,
  5116779.1328282235,
  5117247.056000437,
  5118132.679415019,
  5118261.128204859,
  5118478.937814005,
  5119128.125,
  5119217.1875,
  5119785.825764455,
  5124015.963600086,
  5124304.707913356,
  5124370.098579549,
  5124881.092059543,
  5124911.775081755,
  5126035.781837893,
  5127626.681530933,
  5129140.625,
  5129414.498736343,
  5129651.5625,
  5129696.547701512,
  5130052.235399146,
  5130768.75,
  5131029.6875,
  5131199.604959133,
  5131902.723537738,
  5131969.08454612,
  5132062.88811622,
  5134320.331568867,
  5134728.871460731,
  5136223.4375,
  5140659.941584441,
  5141100.326357179,
  5141117.990946873,
  5141498.206717823,
  5142395.3125,
  5142449.96567222,
  5142633.769782632,
  5142744.78360987,
  5142804.6875,
  5143046.584066203,
  5143387.347854304,
  5143638.822733251,
  5143854.660056851,
  5144356.25,
  5144444.642546348,
  5144589.0625,
  5144602.817701543,
  5144605.477248827,
  5144990.878646369,
  5145059.420718178,
  5145086.2114315815,
  5145370.3125,
  5145383.284397847,
  5145655.891391298,
  5145755.686767004,
  5145761.781111819,
  5145801.1513136625,
  5146231.099199403,
  5146361.458257177,
  5146466.736728932,
  5146586.37906749,
  5146621.82611822,
  5146706.984524036,
  5146741.718585128,
  5146804.4247723725,
  5146855.479355085,
  5146897.620963,
  5147001.877914878,
  5147096.737451783,
  5147191.88517687,
  5147254.955318017,
  5147313.762920726,
  5147448.499477593,
  5147700.235121299,
  5147839.200961352,
  5147848.448028363,
  5147914.7544843005,
  5148065.781773489,
  5148207.461630596,
  5148429.6875,
  5148530.879588512,
  5148617.980795272,
  5148625.470771401,
  5148792.743337879,
  5148895.988933192,
  5148960.9375,
  5149676.867155789,
  5149759.319703105,
  5151380.520745591,
  5151446.105416869,
  5152190.504844647,
  5152366.62584822,
  5152744.703942179,
  5153000.766416478,
  5153104.4182126755,
  5153393.75,
  5153457.600290114,
  5154430.814618393,
  5155607.8125,
  5155617.56739411,
  5155625.0,
  5156575.0,
  5158364.786376255,
  5158650.278738839,
  5158675.961852983,
  5160023.575736752,
  5160367.6540790135,
  5162783.560865514,
  5165412.055243464,
  5167618.54034967,
  5176725.0,
  5180611.517740609,
  5181646.733742297,
  5192275.0,
  5202403.125,
  5207536.565125276,
  5209086.084071497,
  5209101.659606263,
  5229930.358088509,
  5231028.125,
  5236405.181011105,
  5236820.9681162825,
  5238782.884217574,
  5238928.150866442,
  5241415.445744058,
  5248083.826437938,
  5249844.8623016765,
  5250024.358008581,
  5253009.422054304,
  5253775.0,
  5261548.4375,
  5265316.266870025,
  5267368.591254162,
  5267375.206075867,
  5268998.249328477,
  5269252.355843868,
  5269370.3125,
  5269401.5625,
  5269443.630010111,
  5269774.351209665,
  5270050.0,
  5271162.998748115,
  5272202.392863058,
  5272243.9340468235,
  5272246.875,
  5273335.172039438,
  5274944.39831536,
  5283174.041905419,
  5286038.1293557575,
  5290095.992354738,
  5296206.195201195,
  5296207.713871118,
  5296422.0488007,
  5297942.1875,
  5297959.086511491,
  5298113.7839112785,
  5298234.203882818,
  5299062.5,
  5300073.318449883,
  5301200.037387405,
  5301501.5625,
  5302742.1875,
  5302985.292294794,
  5320713.121855185,
  5322578.125,
  5328812.494700109,
  5331696.197855707,
  5332261.657305221,
  5337693.027363412,
  5337718.614399251,
  5339570.760605811,
  5355314.63303687,
  5365040.289862701,
  5376591.261865186,
  5386071.26815407,
  5390136.1325436365,
  5391610.9375,
  5405304.62692363,
  5406559.375,
  5416752.587103924,
  5435344.377887121,
  5438416.787787994,
  5442623.028988975,
  5453554.23859176,
  5456244.805857645,
  5465288.462945165,
  5469689.0625,
  5471504.0964952065,
  5485812.792433178,
  5493850.565478801,
  5512331.191068753,
  5517009.228468536,
  5521196.789962279,
  5522676.168382783,
  5524935.9375,
  5527273.9210504,
  5528385.9375,
  5529890.235479234,
  5530542.993582652,
  5531178.516369288,
  5539693.753242967,
  5547794.814748696,
  5561018.011655362,
  5563015.190728106,
  5565985.569669171,
  5575432.648347869,
  5575469.374932138,
  5580503.3122413205,
  5597721.07133487,
  5612261.95665256,
  5641439.011720579,
  5650640.475272585,
  5655115.520086922,
  5660510.809538817,
  5661272.090823038,
  5661665.469593941,
  5662285.9375,
  5662943.9259487055,
  5662973.009469303,
  5663434.404433025,
  5671745.3125,
  5675628.151393408,
  5681021.875,
  5681710.415618129,
  5684340.295648728,
  5686427.870681439,
  5687791.214183656,
  5690356.595456704,
  5691048.4375,
  5693648.587101451,
  5694278.544786245,
  5695262.5,
  5701011.746802646,
  5706006.801610432,
  5709016.010009828,
  5719288.640826281,
  5724500.0,
  5727643.60315394,
  5732699.495307683,
  5747499.749950118,
  5747500.300883567,
  5747596.62360939,
  5754865.231153525,
  5755421.875,
  5758662.332813297,
  5758682.76462499,
  5759554.6875,
  5760360.9375,
  5760698.5331816785,
  5771644.772551745,
  5772427.634172109,
  5792246.875,
  5795078.032519066,
  5795747.320811018,
  5803335.292934461,
  5804651.410439777,
  5815487.400264339,
  5817965.625,
  5818124.43470538,
  5818200.0,
  5831084.375,
  5838781.601757914,
  5840365.747804157,
  5842970.3125,
  5845291.933640746,
  5848140.625,
  5852334.3672083225,
  5859795.489483133,
  5860608.973634419,
  5865496.991196872,
  5866614.5174266165,
  5866670.6561953975,
  5869540.305823777,
  5869549.838206735,
  5869892.1875,
  5871568.093375797,
  5872000.0,
  5874910.819270179,
  5874958.745822031,
  5874968.916049484,
  5883385.765179415,
  5884901.5625,
  5890356.146796166,
  5900065.958948857,
  5900085.770721988,
  5926671.662790114,
  5947995.261029914,
  5974207.313325031,
  5981109.479358432,
  5982521.349250061,
  5987249.4706328,
  5987253.608432798,
  6037691.172186762,
  6043754.6875,
  6043765.625,
  6047271.875,
  6057065.625,
  6083463.233678017,
  6088084.90542871,
  6096451.047523742,
  6105148.960883869,
  6105162.7231516475,
  6107698.4375,
  6108743.321860101,
  6156955.0551039865,
  6185094.879078507,
  6192882.8125,
  6222450.0,
  6234778.164475562,
  6263332.321396833,
  6276848.125982724,
  6310924.397318744,
  6333128.125,
  6347029.6875,
  6361256.25,
  6361302.620720931,
  6363876.095776006,
  6383368.140021875,
  6596945.277879903,
  6792235.708640398,
  7180392.237749754,
  7190928.125,
  7190937.5,
  7249229.6875,
  7369176.5625,
  7435844.281191614,
  8259487.585124541,
  8295015.804396403,
  8566520.3125,
  8612320.3125,
  8638875.0,
  8638907.198572712,
  9010590.625,
  105461715.15321432,
  409286308.8955168],
 [44.42477307749312,
  49.96609244976557,
  40.87389347834934,
  6.416660493022236,
  53.99412772129534,
  18.51147953711033,
  49.29685290640293,
  8.24947188315964,
  14.759983092403104,
  15.08774115961008,
  5.9256870301143545,
  23.977959645296348,
  52.5315881816371,
  75.29230807725236,
  44.72466138105205,
  11.927069315553487,
  47.71099851526634,
  47.02324246901132,
  156.84240298813347,
  7.117253397298565,
  71.11684278789919,
  6.69717526522069,
  39.88360040306763,
  5.531341860525387,
  154.41072202614725,
  56.397085184607306,
  48.319044820227134,
  18.770221024742526,
  67.16164093873091,
  14.549767833479011,
  43.69399748449416,
  85.25761533493149,
  36.76558753728017,
  27.968108377929276,
  56.437519819866694,
  120.5752728293347,
  71.13826336146273,
  8.260447168588755,
  62.08046277269942,
  16.613998095810693,
  59.13776635557357,
  7.319827029866923,
  7.364350653734773,
  69.15973056167854,
  47.413872665968285,
  92.58409865493707,
  101.60660372620617,
  26.50497962053751,
  6.151349605940142,
  7.3063586266766585,
  163.89181835138777,
  13.065116925778804,
  32.33130019745596,
  34.455816493414886,
  9.189220641010868,
  8.764947767926815,
  79.43192316466265,
  72.62886529136904,
  13.904085510985677,
  5.915055840004793,
  83.18112934327156,
  45.81474067159499,
  24.76604611075996,
  105.53179415813676,
  8.189727717088006,
  53.3582493801643,
  14.214331793834546,
  44.23880686128197,
  89.385686159402,
  62.152328680400814,
  8.885663069648574,
  33.998705913139716,
  58.224245106142725,
  21.390524878587676,
  26.0650561336442,
  152.1060258518192,
  37.02639810704824,
  8.427933468919303,
  69.05314730332555,
  6.502192218445554,
  18.570853939544524,
  33.95059649285611,
  24.816909159076705,
  36.11464670441421,
  21.56088952275647,
  56.34380564508474,
  23.725841937296806,
  18.733550573794172,
  7.153682172343619,
  47.30632576217731,
  72.85410690034762,
  5.653735805873015,
  8.300872719440944,
  26.90854572086416,
  17.746468750421595,
  37.522337956694905,
  92.21999539938129,
  63.55272800046399,
  8.409595016818066,
  13.268052716350402,
  16.74428160044593,
  33.40799817094388,
  82.80651568416958,
  17.708567465364606,
  5.708472262549371,
  31.162506082363564,
  66.79079780214342,
  10.102447256865965,
  62.83779386096224,
  29.704808097780376,
  16.364790848191966,
  41.87933344021094,
  11.920181263548098,
  7.944622374560327,
  69.92587012299553,
  73.64431261133691,
  30.536156037380355,
  28.285278047233597,
  68.34973400802046,
  92.75471294896701,
  17.12784852125473,
  15.683252458747484,
  6.802186053313718,
  10.128335591211648,
  73.57093761112411,
  52.4857470051361,
  14.532068300701411,
  16.410912501860707,
  71.17734772356275,
  18.888498588073,
  14.878988275433365,
  125.842377086918,
  14.379632929063094,
  11.319000423430362,
  74.10033836997975,
  25.71754013625781,
  71.49050589577294,
  117.34423604334512,
  28.37038220070245,
  81.84163068087084,
  71.0766851747359,
  22.608173133296752,
  10.553272582885974,
  64.2517514590782,
  44.432395966905965,
  7.126473320351014,
  116.14751379342539,
  75.93760672031826,
  17.653267416994023,
  7.7257003141274,
  25.05854467101308,
  57.57002974650062,
  82.76895298427584,
  7.590868918783576,
  8.867214371107352,
  48.17216493073142,
  28.657355722719522,
  49.60671178415252,
  14.809124885747092,
  22.079503586088233,
  39.26427399234346,
  53.26045898789833,
  5.226773727087648,
  88.84364277310306,
  33.72349433589182,
  19.287053601784734,
  9.55549747610324,
  15.118676797729096,
  50.7703106398665,
  116.59436326066961,
  25.233652324344742,
  23.02989358557103,
  10.348364652401134,
  8.789489833793379,
  6.684102023871341,
  20.92533071832475,
  12.052528607734255,
  6.418676703138052,
  77.17715622439955,
  27.593306922475765,
  65.25675256492187,
  15.689895000357028,
  24.238813694852315,
  19.160359229874775,
  65.88559311981061,
  11.659312226482983,
  16.85494050797344,
  10.68589801610649,
  25.176541697049075,
  39.31630482109196,
  15.18312088131779,
  15.92040167073918,
  25.07080375845274,
  31.50864675235237,
  31.01999465649287,
  6.795214057494854,
  23.121999092390322,
  57.24715739192866,
  22.891154531330407,
  20.16403729787228,
  17.101300505381584,
  25.019278782390913,
  6.906507750933781,
  31.02256975274764,
  9.603014389412401,
  105.66071307436671,
  16.137206728772213,
  25.435562286269658,
  63.16629552267574,
  13.16030929712963,
  6.9701697769207165,
  122.67411940072881,
  8.149128680072858,
  64.14872004564745,
  45.129760620554784,
  12.499743086677732,
  86.45970410041727,
  11.77839891782654,
  60.19460400417759,
  42.03534379556447,
  12.48536356238823,
  6.8204604263486255,
  13.058598253132624,
  79.77274773816802,
  7.249751160571848,
  31.164652004782994,
  5.901587208411874,
  14.198412381846978,
  54.51179804631836,
  34.81278273675066,
  8.213264264961525,
  31.652169043108945,
  57.45057773315235,
  26.027299144740525,
  6.793125521662823,
  68.09653692684053,
  17.32280264686723,
  71.5788572363732,
  11.997941591951136,
  73.8548765044674,
  70.19616970751521,
  23.189001907835646,
  31.851842694091072,
  42.82807823049187,
  29.459411339099248,
  28.087197356946508,
  11.20206094750102,
  70.10433247809218,
  69.19580145186355,
  40.1702134012693,
  21.838704273804154,
  85.96499394424217,
  78.67133150388811,
  8.755962162840033,
  82.28764158208897,
  47.599723993880765,
  14.503071748369363,
  10.24536037278826,
  39.29326322839205,
  17.25768320880104,
  19.452795927391264,
  7.671188874461247,
  68.84857946513178,
  16.272385048053163,
  6.593760873153204,
  17.88451791532223,
  57.74382893486751,
  18.3706941299025,
  51.76665355504497,
  6.254817225170997,
  11.410688082498234,
  32.73657344433986,
  22.736245580257837,
  19.246571118349145,
  8.198114181846222,
  31.53422323672576,
  100.59066231016293,
  7.131314472582982,
  61.88055534954937,
  72.9641875114753,
  6.1889834937319455,
  6.894347710505344,
  32.56282622804951,
  96.56951937424725,
  27.05964304737305,
  47.44217314407133,
  6.180020700919893,
  5.337296513239306,
  37.79389966033233,
  9.86539767533493,
  7.486360838910614,
  5.364228849553785,
  47.27854356067586,
  15.74096597368917,
  30.656465737352296,
  63.905701320320304,
  59.34748791447035,
  111.72761188401005,
  62.60250734860647,
  22.55547462399196,
  10.544053604547504,
  35.68140491290819,
  80.53030142649618,
  49.82146271426399,
  66.0937694663576,
  6.650986946553472,
  86.53735264185588,
  6.14894666873077,
  107.89280121863177,
  25.931266151148222,
  26.94891679642271,
  75.56408012183371,
  68.85333998443237,
  8.559985482282832,
  63.691029646788635,
  13.639161572916619,
  26.009151477429022,
  16.677395392992437,
  14.953884516854211,
  10.429658250098077,
  138.00425059820108,
  18.556338710559476,
  27.942588659077394,
  25.656333005160427,
  12.033347385784475,
  98.19493626563742,
  33.53547893988432,
  12.90442083450786,
  5.040689238908884,
  62.723149496761295,
  9.244241212944024,
  22.64771036369743,
  21.95062387581726,
  25.294053364567645,
  40.4874824803596,
  34.20575784242585,
  22.242386712972735,
  6.270924193246074,
  50.05542576202717,
  30.49320250073631,
  147.36930740630018,
  85.00285230322807,
  32.92653292115592,
  24.352340380366304,
  5.309640403255765,
  15.915199053699894,
  51.71643672254985,
  8.367305845709916,
  63.98367415685296,
  51.045566192689165,
  31.20501213035617,
  5.419988743630674,
  11.788013189525348,
  97.27917222395683,
  66.23909454626711,
  23.53559935309491,
  30.240776992380425,
  63.52051520997735,
  11.632998700676932,
  10.756258746089319,
  70.49394540294152,
  51.62258860312241,
  25.729364095207302,
  5.7037178737320335,
  97.78418139390564,
  103.2130468329833,
  6.502195507358726,
  21.103641194300113,
  81.99109528255512,
  29.88911322099194,
  90.40398765404383,
  30.426762543990836,
  8.840004125316065,
  5.697148352625648,
  23.220114972545407,
  49.03590769431175,
  66.45174410289283,
  67.75785875415326,
  117.73665211892482,
  16.644000842663353,
  72.39957495838564,
  85.960267272527,
  6.062092028385488,
  71.32858963515577,
  14.91189413095649,
  11.416089729073404,
  13.60455148167163,
  11.314320674529938,
  37.50809757008781,
  78.28151685532504,
  23.50884397843321,
  7.54574958135721,
  70.5350554084501,
  8.831730565031169,
  87.13604870248486,
  5.527369857055622,
  35.80491186809685,
  15.61458045052274,
  14.323009493522108,
  125.72841311383547,
  61.25612427969707,
  8.4714395189329,
  9.194254941847579,
  12.366277002578592,
  27.545958796040416,
  19.25230820209658,
  12.909903731668056,
  6.66289562686799,
  21.205064994367977,
  23.100781330576016,
  57.98689882701888,
  70.36522070258081,
  91.05176625136534,
  48.836232998085286,
  31.6006457136473,
  18.73940213124516,
  12.072092887537075,
  13.28521588157706,
  60.990045209094816,
  56.22617848382819,
  42.67962376222278,
  32.11178450382196,
  16.555939843393517,
  22.18400194651775,
  82.55157499855349,
  63.049798928917596,
  33.11012574999999,
  97.89674012493117,
  23.604730846246586,
  20.526507986082866,
  64.43449775942518,
  34.01014003296156,
  78.19815881172786,
  7.4304929158720965,
  30.83426765823061,
  6.969835923341012,
  27.857701414143335,
  33.982842434147976,
  13.954149107583603,
  49.079138589792635,
  49.62162157954233,
  91.27539384673928,
  29.697659926888676,
  5.078684851720425,
  8.332111499103203,
  185.54946407150334,
  74.85303069547231,
  49.010656213890314,
  86.84510455858383,
  20.427220179174753,
  73.23311408252631,
  188.93552586386815,
  10.452731280726248])

Count of hits during the discharge

caption, T_hit,count1 = 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, count1, caption);

([4694004.6875,
  4808577.372469842,
  4829393.75,
  4832160.4729947215,
  4835189.0625,
  4835197.5757367695,
  4838463.587216016,
  4850960.866667724,
  4851387.279298142,
  4853970.350343791,
  4854046.716112161,
  4854624.717582297,
  4860087.5,
  4865018.75,
  4865732.184104946,
  4866298.8950518835,
  4866820.3125,
  4891287.5,
  4895410.9375,
  4897640.609005003,
  4901851.880021661,
  4906698.095913895,
  4907011.771203193,
  4907448.401623274,
  4910232.661428186,
  4912437.5,
  4916085.417639641,
  4918602.300552371,
  4922876.293468087,
  4928815.518279881,
  4932015.625,
  4943148.4375,
  4955174.565488515,
  4983501.271457321,
  4992959.375,
  4995539.0625,
  5015819.558059051,
  5016340.575258462,
  5023478.125,
  5032326.636951393,
  5034235.9375,
  5036161.641974806,
  5039305.537404446,
  5041470.280810332,
  5066207.655221228,
  5066903.849012415,
  5074100.0,
  5078420.201330655,
  5080255.2745319875,
  5083355.3471233165,
  5092399.822676634,
  5092486.9018794475,
  5092812.5,
  5093047.821994288,
  5093134.201920701,
  5093407.056890622,
  5094117.1875,
  5094623.488912772,
  5094734.073307234,
  5094740.3410823615,
  5095188.582553806,
  5097135.9375,
  5097265.40053432,
  5098068.331329957,
  5098865.437252875,
  5099009.16006795,
  5099365.420171658,
  5099410.750296834,
  5099601.5625,
  5099937.374837869,
  5100758.796193155,
  5100768.069609331,
  5101492.1875,
  5103541.296586633,
  5104068.611881942,
  5104094.034230023,
  5104752.890088723,
  5105220.570571521,
  5105323.4375,
  5105490.720660989,
  5105508.521919656,
  5105670.3125,
  5105671.654121327,
  5106156.772065578,
  5106372.944268811,
  5106595.3125,
  5107134.549287804,
  5107575.198086805,
  5107876.683707601,
  5107878.872365402,
  5107975.0,
  5108042.8833440915,
  5108192.2144353185,
  5108380.900189651,
  5108586.00717675,
  5109306.25,
  5109388.417477825,
  5109615.625,
  5109736.163102475,
  5110691.666187518,
  5111401.966758304,
  5111490.625,
  5112127.963118512,
  5112895.888247758,
  5112950.026901975,
  5113793.384749336,
  5114179.6875,
  5114364.852440559,
  5114462.5,
  5114473.4375,
  5115124.227090571,
  5115191.536472652,
  5115344.204021185,
  5115502.417075702,
  5115741.348361351,
  5116485.328295426,
  5116562.5,
  5116593.489144347,
  5116613.610862971,
  5116779.1328282235,
  5117247.056000437,
  5118132.679415019,
  5118261.128204859,
  5118478.937814005,
  5119128.125,
  5119217.1875,
  5119785.825764455,
  5124015.963600086,
  5124304.707913356,
  5124370.098579549,
  5124881.092059543,
  5124911.775081755,
  5126035.781837893,
  5127626.681530933,
  5129140.625,
  5129414.498736343,
  5129651.5625,
  5129696.547701512,
  5130052.235399146,
  5130768.75,
  5131029.6875,
  5131199.604959133,
  5131902.723537738,
  5131969.08454612,
  5132062.88811622,
  5134320.331568867,
  5134728.871460731,
  5136223.4375,
  5140659.941584441,
  5141100.326357179,
  5141117.990946873,
  5141498.206717823,
  5142395.3125,
  5142449.96567222,
  5142633.769782632,
  5142744.78360987,
  5142804.6875,
  5143046.584066203,
  5143387.347854304,
  5143638.822733251,
  5143854.660056851,
  5144356.25,
  5144444.642546348,
  5144589.0625,
  5144602.817701543,
  5144605.477248827,
  5144990.878646369,
  5145059.420718178,
  5145086.2114315815,
  5145370.3125,
  5145383.284397847,
  5145655.891391298,
  5145755.686767004,
  5145761.781111819,
  5145801.1513136625,
  5146231.099199403,
  5146361.458257177,
  5146466.736728932,
  5146586.37906749,
  5146621.82611822,
  5146706.984524036,
  5146741.718585128,
  5146804.4247723725,
  5146855.479355085,
  5146897.620963,
  5147001.877914878,
  5147096.737451783,
  5147191.88517687,
  5147254.955318017,
  5147313.762920726,
  5147448.499477593,
  5147700.235121299,
  5147839.200961352,
  5147848.448028363,
  5147914.7544843005,
  5148065.781773489,
  5148207.461630596,
  5148429.6875,
  5148530.879588512,
  5148617.980795272,
  5148625.470771401,
  5148792.743337879,
  5148895.988933192,
  5148960.9375,
  5149676.867155789,
  5149759.319703105,
  5151380.520745591,
  5151446.105416869,
  5152190.504844647,
  5152366.62584822,
  5152744.703942179,
  5153000.766416478,
  5153104.4182126755,
  5153393.75,
  5153457.600290114,
  5154430.814618393,
  5155607.8125,
  5155617.56739411,
  5155625.0,
  5156575.0,
  5158364.786376255,
  5158650.278738839,
  5158675.961852983,
  5160023.575736752,
  5160367.6540790135,
  5162783.560865514,
  5165412.055243464,
  5167618.54034967,
  5176725.0,
  5180611.517740609,
  5181646.733742297,
  5192275.0,
  5202403.125,
  5207536.565125276,
  5209086.084071497,
  5209101.659606263,
  5229930.358088509,
  5231028.125,
  5236405.181011105,
  5236820.9681162825,
  5238782.884217574,
  5238928.150866442,
  5241415.445744058,
  5248083.826437938,
  5249844.8623016765,
  5250024.358008581,
  5253009.422054304,
  5253775.0,
  5261548.4375,
  5265316.266870025,
  5267368.591254162,
  5267375.206075867,
  5268998.249328477,
  5269252.355843868,
  5269370.3125,
  5269401.5625,
  5269443.630010111,
  5269774.351209665,
  5270050.0,
  5271162.998748115,
  5272202.392863058,
  5272243.9340468235,
  5272246.875,
  5273335.172039438,
  5274944.39831536,
  5283174.041905419,
  5286038.1293557575,
  5290095.992354738,
  5296206.195201195,
  5296207.713871118,
  5296422.0488007,
  5297942.1875,
  5297959.086511491,
  5298113.7839112785,
  5298234.203882818,
  5299062.5,
  5300073.318449883,
  5301200.037387405,
  5301501.5625,
  5302742.1875,
  5302985.292294794,
  5320713.121855185,
  5322578.125,
  5328812.494700109,
  5331696.197855707,
  5332261.657305221,
  5337693.027363412,
  5337718.614399251,
  5339570.760605811,
  5355314.63303687,
  5365040.289862701,
  5376591.261865186,
  5386071.26815407,
  5390136.1325436365,
  5391610.9375,
  5405304.62692363,
  5406559.375,
  5416752.587103924,
  5435344.377887121,
  5438416.787787994,
  5442623.028988975,
  5453554.23859176,
  5456244.805857645,
  5465288.462945165,
  5469689.0625,
  5471504.0964952065,
  5485812.792433178,
  5493850.565478801,
  5512331.191068753,
  5517009.228468536,
  5521196.789962279,
  5522676.168382783,
  5524935.9375,
  5527273.9210504,
  5528385.9375,
  5529890.235479234,
  5530542.993582652,
  5531178.516369288,
  5539693.753242967,
  5547794.814748696,
  5561018.011655362,
  5563015.190728106,
  5565985.569669171,
  5575432.648347869,
  5575469.374932138,
  5580503.3122413205,
  5597721.07133487,
  5612261.95665256,
  5641439.011720579,
  5650640.475272585,
  5655115.520086922,
  5660510.809538817,
  5661272.090823038,
  5661665.469593941,
  5662285.9375,
  5662943.9259487055,
  5662973.009469303,
  5663434.404433025,
  5671745.3125,
  5675628.151393408,
  5681021.875,
  5681710.415618129,
  5684340.295648728,
  5686427.870681439,
  5687791.214183656,
  5690356.595456704,
  5691048.4375,
  5693648.587101451,
  5694278.544786245,
  5695262.5,
  5701011.746802646,
  5706006.801610432,
  5709016.010009828,
  5719288.640826281,
  5724500.0,
  5727643.60315394,
  5732699.495307683,
  5747499.749950118,
  5747500.300883567,
  5747596.62360939,
  5754865.231153525,
  5755421.875,
  5758662.332813297,
  5758682.76462499,
  5759554.6875,
  5760360.9375,
  5760698.5331816785,
  5771644.772551745,
  5772427.634172109,
  5792246.875,
  5795078.032519066,
  5795747.320811018,
  5803335.292934461,
  5804651.410439777,
  5815487.400264339,
  5817965.625,
  5818124.43470538,
  5818200.0,
  5831084.375,
  5838781.601757914,
  5840365.747804157,
  5842970.3125,
  5845291.933640746,
  5848140.625,
  5852334.3672083225,
  5859795.489483133,
  5860608.973634419,
  5865496.991196872,
  5866614.5174266165,
  5866670.6561953975,
  5869540.305823777,
  5869549.838206735,
  5869892.1875,
  5871568.093375797,
  5872000.0,
  5874910.819270179,
  5874958.745822031,
  5874968.916049484,
  5883385.765179415,
  5884901.5625,
  5890356.146796166,
  5900065.958948857,
  5900085.770721988,
  5926671.662790114,
  5947995.261029914,
  5974207.313325031,
  5981109.479358432,
  5982521.349250061,
  5987249.4706328,
  5987253.608432798,
  6037691.172186762,
  6043754.6875,
  6043765.625,
  6047271.875,
  6057065.625,
  6083463.233678017,
  6088084.90542871,
  6096451.047523742,
  6105148.960883869,
  6105162.7231516475,
  6107698.4375,
  6108743.321860101,
  6156955.0551039865,
  6185094.879078507,
  6192882.8125,
  6222450.0,
  6234778.164475562,
  6263332.321396833,
  6276848.125982724,
  6310924.397318744,
  6333128.125,
  6347029.6875,
  6361256.25,
  6361302.620720931,
  6363876.095776006,
  6383368.140021875,
  6596945.277879903,
  6792235.708640398,
  7180392.237749754,
  7190928.125,
  7190937.5,
  7249229.6875,
  7369176.5625,
  7435844.281191614,
  8259487.585124541,
  8295015.804396403,
  8566520.3125,
  8612320.3125,
  8638875.0,
  8638907.198572712,
  9010590.625,
  105461715.15321432,
  409286308.8955168],
 [44.42477307749312,
  49.96609244976557,
  40.87389347834934,
  6.416660493022236,
  53.99412772129534,
  18.51147953711033,
  49.29685290640293,
  8.24947188315964,
  14.759983092403104,
  15.08774115961008,
  5.9256870301143545,
  23.977959645296348,
  52.5315881816371,
  75.29230807725236,
  44.72466138105205,
  11.927069315553487,
  47.71099851526634,
  47.02324246901132,
  156.84240298813347,
  7.117253397298565,
  71.11684278789919,
  6.69717526522069,
  39.88360040306763,
  5.531341860525387,
  154.41072202614725,
  56.397085184607306,
  48.319044820227134,
  18.770221024742526,
  67.16164093873091,
  14.549767833479011,
  43.69399748449416,
  85.25761533493149,
  36.76558753728017,
  27.968108377929276,
  56.437519819866694,
  120.5752728293347,
  71.13826336146273,
  8.260447168588755,
  62.08046277269942,
  16.613998095810693,
  59.13776635557357,
  7.319827029866923,
  7.364350653734773,
  69.15973056167854,
  47.413872665968285,
  92.58409865493707,
  101.60660372620617,
  26.50497962053751,
  6.151349605940142,
  7.3063586266766585,
  163.89181835138777,
  13.065116925778804,
  32.33130019745596,
  34.455816493414886,
  9.189220641010868,
  8.764947767926815,
  79.43192316466265,
  72.62886529136904,
  13.904085510985677,
  5.915055840004793,
  83.18112934327156,
  45.81474067159499,
  24.76604611075996,
  105.53179415813676,
  8.189727717088006,
  53.3582493801643,
  14.214331793834546,
  44.23880686128197,
  89.385686159402,
  62.152328680400814,
  8.885663069648574,
  33.998705913139716,
  58.224245106142725,
  21.390524878587676,
  26.0650561336442,
  152.1060258518192,
  37.02639810704824,
  8.427933468919303,
  69.05314730332555,
  6.502192218445554,
  18.570853939544524,
  33.95059649285611,
  24.816909159076705,
  36.11464670441421,
  21.56088952275647,
  56.34380564508474,
  23.725841937296806,
  18.733550573794172,
  7.153682172343619,
  47.30632576217731,
  72.85410690034762,
  5.653735805873015,
  8.300872719440944,
  26.90854572086416,
  17.746468750421595,
  37.522337956694905,
  92.21999539938129,
  63.55272800046399,
  8.409595016818066,
  13.268052716350402,
  16.74428160044593,
  33.40799817094388,
  82.80651568416958,
  17.708567465364606,
  5.708472262549371,
  31.162506082363564,
  66.79079780214342,
  10.102447256865965,
  62.83779386096224,
  29.704808097780376,
  16.364790848191966,
  41.87933344021094,
  11.920181263548098,
  7.944622374560327,
  69.92587012299553,
  73.64431261133691,
  30.536156037380355,
  28.285278047233597,
  68.34973400802046,
  92.75471294896701,
  17.12784852125473,
  15.683252458747484,
  6.802186053313718,
  10.128335591211648,
  73.57093761112411,
  52.4857470051361,
  14.532068300701411,
  16.410912501860707,
  71.17734772356275,
  18.888498588073,
  14.878988275433365,
  125.842377086918,
  14.379632929063094,
  11.319000423430362,
  74.10033836997975,
  25.71754013625781,
  71.49050589577294,
  117.34423604334512,
  28.37038220070245,
  81.84163068087084,
  71.0766851747359,
  22.608173133296752,
  10.553272582885974,
  64.2517514590782,
  44.432395966905965,
  7.126473320351014,
  116.14751379342539,
  75.93760672031826,
  17.653267416994023,
  7.7257003141274,
  25.05854467101308,
  57.57002974650062,
  82.76895298427584,
  7.590868918783576,
  8.867214371107352,
  48.17216493073142,
  28.657355722719522,
  49.60671178415252,
  14.809124885747092,
  22.079503586088233,
  39.26427399234346,
  53.26045898789833,
  5.226773727087648,
  88.84364277310306,
  33.72349433589182,
  19.287053601784734,
  9.55549747610324,
  15.118676797729096,
  50.7703106398665,
  116.59436326066961,
  25.233652324344742,
  23.02989358557103,
  10.348364652401134,
  8.789489833793379,
  6.684102023871341,
  20.92533071832475,
  12.052528607734255,
  6.418676703138052,
  77.17715622439955,
  27.593306922475765,
  65.25675256492187,
  15.689895000357028,
  24.238813694852315,
  19.160359229874775,
  65.88559311981061,
  11.659312226482983,
  16.85494050797344,
  10.68589801610649,
  25.176541697049075,
  39.31630482109196,
  15.18312088131779,
  15.92040167073918,
  25.07080375845274,
  31.50864675235237,
  31.01999465649287,
  6.795214057494854,
  23.121999092390322,
  57.24715739192866,
  22.891154531330407,
  20.16403729787228,
  17.101300505381584,
  25.019278782390913,
  6.906507750933781,
  31.02256975274764,
  9.603014389412401,
  105.66071307436671,
  16.137206728772213,
  25.435562286269658,
  63.16629552267574,
  13.16030929712963,
  6.9701697769207165,
  122.67411940072881,
  8.149128680072858,
  64.14872004564745,
  45.129760620554784,
  12.499743086677732,
  86.45970410041727,
  11.77839891782654,
  60.19460400417759,
  42.03534379556447,
  12.48536356238823,
  6.8204604263486255,
  13.058598253132624,
  79.77274773816802,
  7.249751160571848,
  31.164652004782994,
  5.901587208411874,
  14.198412381846978,
  54.51179804631836,
  34.81278273675066,
  8.213264264961525,
  31.652169043108945,
  57.45057773315235,
  26.027299144740525,
  6.793125521662823,
  68.09653692684053,
  17.32280264686723,
  71.5788572363732,
  11.997941591951136,
  73.8548765044674,
  70.19616970751521,
  23.189001907835646,
  31.851842694091072,
  42.82807823049187,
  29.459411339099248,
  28.087197356946508,
  11.20206094750102,
  70.10433247809218,
  69.19580145186355,
  40.1702134012693,
  21.838704273804154,
  85.96499394424217,
  78.67133150388811,
  8.755962162840033,
  82.28764158208897,
  47.599723993880765,
  14.503071748369363,
  10.24536037278826,
  39.29326322839205,
  17.25768320880104,
  19.452795927391264,
  7.671188874461247,
  68.84857946513178,
  16.272385048053163,
  6.593760873153204,
  17.88451791532223,
  57.74382893486751,
  18.3706941299025,
  51.76665355504497,
  6.254817225170997,
  11.410688082498234,
  32.73657344433986,
  22.736245580257837,
  19.246571118349145,
  8.198114181846222,
  31.53422323672576,
  100.59066231016293,
  7.131314472582982,
  61.88055534954937,
  72.9641875114753,
  6.1889834937319455,
  6.894347710505344,
  32.56282622804951,
  96.56951937424725,
  27.05964304737305,
  47.44217314407133,
  6.180020700919893,
  5.337296513239306,
  37.79389966033233,
  9.86539767533493,
  7.486360838910614,
  5.364228849553785,
  47.27854356067586,
  15.74096597368917,
  30.656465737352296,
  63.905701320320304,
  59.34748791447035,
  111.72761188401005,
  62.60250734860647,
  22.55547462399196,
  10.544053604547504,
  35.68140491290819,
  80.53030142649618,
  49.82146271426399,
  66.0937694663576,
  6.650986946553472,
  86.53735264185588,
  6.14894666873077,
  107.89280121863177,
  25.931266151148222,
  26.94891679642271,
  75.56408012183371,
  68.85333998443237,
  8.559985482282832,
  63.691029646788635,
  13.639161572916619,
  26.009151477429022,
  16.677395392992437,
  14.953884516854211,
  10.429658250098077,
  138.00425059820108,
  18.556338710559476,
  27.942588659077394,
  25.656333005160427,
  12.033347385784475,
  98.19493626563742,
  33.53547893988432,
  12.90442083450786,
  5.040689238908884,
  62.723149496761295,
  9.244241212944024,
  22.64771036369743,
  21.95062387581726,
  25.294053364567645,
  40.4874824803596,
  34.20575784242585,
  22.242386712972735,
  6.270924193246074,
  50.05542576202717,
  30.49320250073631,
  147.36930740630018,
  85.00285230322807,
  32.92653292115592,
  24.352340380366304,
  5.309640403255765,
  15.915199053699894,
  51.71643672254985,
  8.367305845709916,
  63.98367415685296,
  51.045566192689165,
  31.20501213035617,
  5.419988743630674,
  11.788013189525348,
  97.27917222395683,
  66.23909454626711,
  23.53559935309491,
  30.240776992380425,
  63.52051520997735,
  11.632998700676932,
  10.756258746089319,
  70.49394540294152,
  51.62258860312241,
  25.729364095207302,
  5.7037178737320335,
  97.78418139390564,
  103.2130468329833,
  6.502195507358726,
  21.103641194300113,
  81.99109528255512,
  29.88911322099194,
  90.40398765404383,
  30.426762543990836,
  8.840004125316065,
  5.697148352625648,
  23.220114972545407,
  49.03590769431175,
  66.45174410289283,
  67.75785875415326,
  117.73665211892482,
  16.644000842663353,
  72.39957495838564,
  85.960267272527,
  6.062092028385488,
  71.32858963515577,
  14.91189413095649,
  11.416089729073404,
  13.60455148167163,
  11.314320674529938,
  37.50809757008781,
  78.28151685532504,
  23.50884397843321,
  7.54574958135721,
  70.5350554084501,
  8.831730565031169,
  87.13604870248486,
  5.527369857055622,
  35.80491186809685,
  15.61458045052274,
  14.323009493522108,
  125.72841311383547,
  61.25612427969707,
  8.4714395189329,
  9.194254941847579,
  12.366277002578592,
  27.545958796040416,
  19.25230820209658,
  12.909903731668056,
  6.66289562686799,
  21.205064994367977,
  23.100781330576016,
  57.98689882701888,
  70.36522070258081,
  91.05176625136534,
  48.836232998085286,
  31.6006457136473,
  18.73940213124516,
  12.072092887537075,
  13.28521588157706,
  60.990045209094816,
  56.22617848382819,
  42.67962376222278,
  32.11178450382196,
  16.555939843393517,
  22.18400194651775,
  82.55157499855349,
  63.049798928917596,
  33.11012574999999,
  97.89674012493117,
  23.604730846246586,
  20.526507986082866,
  64.43449775942518,
  34.01014003296156,
  78.19815881172786,
  7.4304929158720965,
  30.83426765823061,
  6.969835923341012,
  27.857701414143335,
  33.982842434147976,
  13.954149107583603,
  49.079138589792635,
  49.62162157954233,
  91.27539384673928,
  29.697659926888676,
  5.078684851720425,
  8.332111499103203,
  185.54946407150334,
  74.85303069547231,
  49.010656213890314,
  86.84510455858383,
  20.427220179174753,
  73.23311408252631,
  188.93552586386815,
  10.452731280726248])

Graph for home page

Detected energies during the discharge + Energy spectrum

multiplot(icon_fig, T_int_first,E,xmean,count)

Average size of interactions

size_of_interactions_average(Etot, Ntot,size_interaction)

Energy spectra for parts of doublebreadown discharges

#energy_spectra_doublebreakdown(Etot, T_first, file_energy_spectra_doublebreakdown)