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

([4865392.1875,
  4866430.372746798,
  4902988.877834031,
  4931209.506342441,
  4937882.430525896,
  4989789.699282084,
  5073603.342721403,
  5079138.615780119,
  5093085.45832693,
  5104781.20245201,
  5118369.247009864,
  5124782.27535554,
  5140497.130575997,
  5146756.0449977415,
  5150816.348784119,
  5153510.32937531,
  5167977.597559522,
  5172818.656928344,
  5193554.6875,
  5199787.575826112,
  5200972.652611837,
  5203475.506271819,
  5220719.785682488,
  5231376.29162031,
  5237002.271631566,
  5245253.370920323,
  5245967.1875,
  5247632.7693357095,
  5247760.907177379,
  5251951.340233685,
  5255084.084850079,
  5255424.468436674,
  5256010.9375,
  5263303.125,
  5265394.059089717,
  5275656.25,
  5276981.663172591,
  5279809.096651305,
  5280143.868922852,
  5280217.973431971,
  5282740.7915373845,
  5290157.5582990125,
  5291563.740521363,
  5296328.34780919,
  5298732.450366807,
  5301565.493446072,
  5301655.840839303,
  5305144.112784824,
  5317597.998484296,
  5317853.277808623,
  5319240.625,
  5325339.72645843,
  5327757.3651289325,
  5327767.1875,
  5327929.6875,
  5337276.5625,
  5338517.159864835,
  5342450.874912309,
  5349233.136463841,
  5360171.875,
  5374894.392739861,
  5377357.8125,
  5391899.3522231085,
  5393501.44976538,
  5393943.444008141,
  5396660.674678239,
  5406919.580455997,
  5407604.6875,
  5434826.5625,
  5441242.1875,
  5444014.977607985,
  5444230.828413099,
  5446717.084932526,
  5446754.99661403,
  5449906.246461507,
  5449924.498939607,
  5454240.625,
  5462168.74143095,
  5476964.340195499,
  5479483.608329411,
  5506009.027312048,
  5511290.801604016,
  5511955.605892921,
  5527270.8912632475,
  5528318.75,
  5564456.068549881,
  5575449.680018144,
  5581357.50364499,
  5595592.550203723,
  5598322.042129034,
  5617149.081652434,
  5624150.0,
  5648496.875,
  5649247.355164573,
  5649866.754714486,
  5653337.088260016,
  5654041.669712528,
  5657492.1875,
  5658982.774321432,
  5677289.376096702,
  5680786.293636969,
  5684165.625,
  5686352.953112831,
  5689164.601985192,
  5695773.73125102,
  5695848.7187655205,
  5706056.673345632,
  5714988.146156091,
  5722040.125979857,
  5722330.90480796,
  5743976.335496414,
  5751294.984614641,
  5751347.838897046,
  5760979.496384642,
  5774444.461648831,
  5775406.25,
  5796264.774612336,
  5815738.628183398,
  5825473.372968249,
  5840247.504002235,
  5878079.509601273,
  5890681.25,
  5910840.968365635,
  5912360.9375,
  5936604.598761982,
  5953924.510853552,
  5973079.6875,
  5983315.625,
  5997006.25,
  6001074.010843979,
  6002660.931751104,
  6006582.404836634,
  6011543.148967736,
  6014060.780350795,
  6027293.106376829,
  6029378.316641708,
  6063359.582652451,
  6082106.761413057,
  6082450.317008664,
  6087117.098017578,
  6097646.694387269,
  6105735.9375,
  6119614.9334454425,
  6127973.469032305,
  6128865.075723356,
  6141499.462182029,
  6149806.49389333,
  6154144.799703064,
  6193469.929502385,
  6202392.1875,
  6214989.698152065,
  6220368.830306384,
  6224790.625,
  6232175.303909209,
  6240180.36032181,
  6257892.221069465,
  6281338.799058635,
  6287671.875,
  6296088.45071476,
  6300289.367994876,
  6307860.05450089,
  6310993.857144763,
  6316745.190537277,
  6319001.440458031,
  6319027.418868509,
  6325835.83252889,
  6326250.503506313,
  6326490.625,
  6330628.724554141,
  6356588.320240547,
  6360415.335313595,
  6363562.221853637,
  6371274.374273011,
  6386523.422946879,
  6394667.1875,
  6394808.378273395,
  6394951.731885219,
  6403358.2815756025,
  6421995.7453817045,
  6437765.826745795,
  6438092.1875,
  6447632.190745184,
  6447646.992472689,
  6455367.293921407,
  6459441.354706894,
  6464406.181286556,
  6472130.82115724,
  6475175.0,
  6475564.552828369,
  6475599.2993196,
  6477937.928923169,
  6535332.618258111,
  6535743.835697867,
  6547525.0,
  6555413.456796719,
  6562930.399962301,
  6603127.438495427,
  6615900.671521435,
  6635863.955488209,
  6636370.300816379,
  6642985.021600708,
  6643038.471486368,
  6645496.875,
  6677741.959400974,
  6678039.534533036,
  6678332.126183774,
  6679654.6875,
  6680502.849347608,
  6680672.247446646,
  6685620.630531663,
  6688166.228500593,
  6695758.554384216,
  6713062.118196205,
  6713070.088983797,
  6725877.498739352,
  6733785.048467853,
  6736120.832033469,
  6756957.662575316,
  6758993.164265754,
  6759862.5,
  6767173.048328394,
  6769968.1656420585,
  6772373.67911551,
  6773278.756659159,
  6774296.8185508605,
  6774812.5,
  6776609.226339851,
  6822439.211374638,
  6822690.625,
  6827356.25,
  6827370.032493986,
  6838047.811853215,
  6864329.6875,
  6892826.139903439,
  6905046.199146711,
  6935284.180576586,
  6936528.697655473,
  6936781.25,
  6942010.130280411,
  6945760.9375,
  6949218.227563352,
  6951908.091333866,
  7000270.03845738,
  7000946.124665435,
  7008157.8125,
  7010532.065798446,
  7013267.952952388,
  7013883.764522348,
  7015869.7056191275,
  7019588.923166189,
  7053725.0,
  7063185.739657797,
  7067185.264894794,
  7072989.0625,
  7076690.728302234,
  7087073.196209238,
  7089940.096461754,
  7090184.375,
  7093756.193390366,
  7095231.760088102,
  7096520.08635461,
  7104253.53146217,
  7151514.0625,
  7159301.883181753,
  7296915.296155198,
  7398592.1875,
  7417698.280867868,
  7610501.410732309,
  7613122.4067467395,
  7633831.993895526,
  7694799.539224773,
  7698179.8815623205,
  7701135.712461473,
  7701135.9375,
  7710684.47880043,
  7716047.632769352,
  7727164.0625,
  7729813.215483495,
  7739449.389936655,
  7744217.940267222,
  7749068.75,
  7825326.5625,
  7847084.851139585,
  7851373.718111622,
  7854680.821309318,
  7917642.167052413,
  7918303.125,
  7919589.165673966,
  7921003.220929531,
  7921179.7965270085,
  7921527.851020623,
  7921530.550743104,
  7927358.943530014,
  7939742.105494164,
  7954414.0625,
  8007575.436703201,
  8007686.043733773,
  8022243.8229530845,
  8022659.496952706,
  8025723.957001001,
  8045795.318675183,
  8048839.353779782,
  8117460.786861867,
  8119758.507532099,
  8119759.972393939,
  8120473.281012251,
  8138021.6670857165,
  8166320.73808537,
  8201688.426360736,
  8206147.3120294465,
  8208581.975229668,
  8216937.580248538,
  8223488.295775776,
  8240317.1875,
  8294933.041968273,
  8320040.635791321,
  8326085.875880561,
  8328465.651843779,
  8369207.567085598,
  8501682.92761445,
  8501692.27164109,
  8501721.226743998,
  8518421.617780104,
  8549568.690662805,
  8552057.838524133,
  8553351.275278158,
  8562338.149069987,
  8568493.75,
  8597172.078720665,
  8608918.75,
  8656933.469497064,
  8659381.923756963,
  8681962.5,
  8745687.083034104,
  8782968.75,
  8783012.973672677,
  8786482.8125,
  8797340.17574899,
  8802787.551015569,
  8833018.512568837,
  8834316.016556965,
  8836015.625,
  8836911.6536696,
  8838725.0,
  8843837.602300921,
  8849181.386757499,
  8853117.845908893,
  8853892.326921102,
  8888759.375,
  8897020.3125,
  8910762.5,
  8911898.208020471,
  8916158.068736866,
  8917901.5625,
  8930655.663681623,
  8936133.845905786,
  8936922.108446503,
  8938055.517940393,
  8938203.125,
  8938923.228131622,
  8939000.355089176,
  8941153.067436485,
  8941657.231177805,
  8941778.020609286,
  8945236.18662382,
  8949599.838124905,
  8952072.544086674,
  8954128.70838222,
  8954268.40539949,
  8957745.79054372,
  8963076.382611705,
  8963167.81675824,
  8964237.5,
  8964747.823468974,
  8966033.80114346,
  8967071.0880002,
  8967099.415342307,
  8967130.264751952,
  8967186.834550261,
  8967311.329356162,
  8967542.1875,
  8968271.623086121,
  8980504.366962962,
  8984111.785309391,
  8987104.404950673,
  8991674.51527097,
  9054337.5,
  9106188.932355177,
  9125989.137378976,
  9130365.234000752,
  9132952.931051036,
  9150770.510395326,
  9155765.642201256,
  9168907.542526789,
  9178734.982517647,
  9233720.775440203,
  9277111.54348997,
  9472780.894757938,
  9900853.125,
  9914369.428507341,
  9915280.378040625,
  9925249.240184313,
  9928495.228508929,
  9936715.402180506,
  9937158.216579778,
  9941225.07845454,
  9944835.293526413,
  9945089.414011562,
  9945108.780579712,
  9947621.875,
  9948051.21522697,
  9951507.778116597,
  9956662.717577567,
  9963130.908218045,
  9963538.264386438,
  9963541.891482808,
  9967531.291752014,
  9967820.653793678,
  9969878.916244522,
  9970260.21780178,
  9971392.061634675,
  9973409.375,
  9974470.880837968,
  9977762.708106669,
  9980772.483151233,
  9985644.409075178,
  9986233.12707635,
  9988371.875,
  9995728.125,
  10001225.191173948,
  10006386.880484438,
  10010594.159571325,
  10013183.349512633,
  10014871.659683162,
  10015761.451107955,
  10021789.0625,
  10049328.194706488,
  10053358.641931511,
  10066326.120938992],
 [41.92171833418734,
  18.406778570380563,
  17.182598204863623,
  16.04072393902368,
  23.366057742176505,
  111.73035750812565,
  32.17263621088415,
  13.474375444383824,
  13.383450326542999,
  27.61778155813478,
  6.866960046466185,
  5.507567019094055,
  16.014175336720857,
  22.542489511418584,
  7.887822436325474,
  34.24940894886531,
  21.223665630181195,
  12.521139068320783,
  74.32773542709742,
  11.248333040615185,
  35.725563865620295,
  102.52744007144821,
  89.77889144838295,
  7.508117754744114,
  18.935935858196718,
  26.574759864928197,
  57.314746968540994,
  27.69858528190561,
  61.12364519878206,
  32.51491457194023,
  112.73840521746234,
  16.095115607142255,
  30.055594216079108,
  60.18235688516705,
  11.64783523824972,
  34.29117961954461,
  26.543888820978193,
  20.35007143511837,
  17.048044030941227,
  5.758671659712053,
  154.964331301543,
  59.17653832191581,
  13.841563984915918,
  9.527185249402326,
  22.988247212687696,
  26.170595781605172,
  9.007845191825291,
  81.60102991448929,
  96.72964398215204,
  62.29544449035191,
  63.12351423519161,
  42.581263059289,
  22.03467426212426,
  29.142767456464163,
  54.01520478968915,
  32.3446403650468,
  63.35122254770762,
  82.552795834706,
  16.33902645932599,
  56.66032586978959,
  42.6712759946148,
  55.47676703708981,
  52.184099426065124,
  14.5285498675742,
  23.674504510823443,
  5.388520763729763,
  19.453515682549053,
  38.70423609105755,
  51.03558424211106,
  30.29836203756138,
  40.99372581897157,
  23.876864381750938,
  14.564469163866981,
  16.266961002626093,
  14.928039472812234,
  21.482528768601814,
  91.87485407327918,
  22.286208876610225,
  37.6336973101285,
  12.650575376484758,
  156.98218384316613,
  9.485333663940402,
  20.16669359586452,
  162.7375945311784,
  43.91404499038363,
  14.026553857262552,
  29.715968347819018,
  64.40919262298496,
  127.8340116175549,
  66.24357607086449,
  7.854107644685509,
  43.8442993165237,
  90.08803383478521,
  5.73801395639558,
  24.078380117812017,
  22.42116708322602,
  13.277592907349314,
  39.832070479177574,
  27.79117039213128,
  72.32562738142042,
  49.30422006886399,
  207.79068931228366,
  18.323258252264726,
  12.090503156031742,
  8.448276550897713,
  5.046315850607294,
  72.31670914268936,
  6.899747738586512,
  24.11331186132422,
  26.295001169551043,
  24.73021146824968,
  221.4464167445425,
  7.993511107139844,
  20.196898697076485,
  8.702537247997737,
  32.423071887914034,
  78.61875089229471,
  90.8565925859285,
  7.104590919145231,
  17.99112147107802,
  85.31308960311398,
  36.72267188223688,
  35.65344189825217,
  69.45610045536584,
  26.88497676711342,
  10.440504649978429,
  68.29918668074761,
  155.24336803635688,
  71.9949963402927,
  25.694627540462186,
  22.36170146782888,
  23.97751578058194,
  32.65473515789066,
  53.86486148178036,
  30.93897808313987,
  17.06299458372032,
  14.310452392930296,
  152.61754044891873,
  11.662209833367863,
  18.733156659837153,
  82.89903039663994,
  91.8872044592244,
  17.205809415186046,
  66.65256591826027,
  22.53953929023367,
  21.123962669226557,
  40.00982907570861,
  10.410785603410375,
  25.768047455859197,
  103.30971703624303,
  34.03358804937414,
  11.255594593336562,
  80.02301407409479,
  18.054507105333013,
  44.41669704364794,
  20.110729860372608,
  62.62807676544802,
  95.52481211777905,
  20.44603768204785,
  40.095057170367056,
  65.83613246124716,
  106.67755094758725,
  16.513157022837653,
  40.240195836001504,
  17.6362914187126,
  5.528007558688627,
  6.332980053958412,
  71.11715682504537,
  7.8355036025151525,
  44.99224835778598,
  41.09603295048318,
  14.324079872713874,
  85.17457039982077,
  39.56788363755811,
  54.54017601491058,
  213.41744226137249,
  76.98631210072182,
  115.94930936287884,
  53.02697181438315,
  119.61291458060569,
  68.63845579837657,
  20.358829046692907,
  76.17867945682896,
  57.84908793124893,
  28.259205181569705,
  61.05173265173711,
  47.52541668551477,
  50.98278050150255,
  17.21021739233419,
  13.350754869949611,
  47.93421246628941,
  97.68083294676879,
  22.042131247716846,
  47.86036492203617,
  20.499707251894527,
  7.882989602247322,
  10.201067856043407,
  18.31528909845606,
  34.024119705998935,
  28.295759354659463,
  8.66091351814023,
  5.4369587574684255,
  31.45163209026017,
  9.1457029356276,
  17.108271175313337,
  19.816160940109455,
  68.87824639675871,
  10.722548843222299,
  6.315429758225304,
  19.648170141427272,
  101.5891540340275,
  5.222744089946213,
  67.88285399545867,
  9.149297886930816,
  6.966014626835337,
  141.2117462592151,
  9.813057895638435,
  25.660286854180207,
  30.263023988025257,
  44.956186737419365,
  10.56941302518406,
  62.18650836042743,
  6.298215119365958,
  12.282985189294758,
  27.455416535526005,
  58.224677592197324,
  38.990191625499115,
  15.541027777452356,
  43.193792243072224,
  72.91008299765167,
  24.009541643040983,
  93.93204348292733,
  40.79289963785829,
  72.89029492208391,
  18.59355293415175,
  29.1037590842152,
  8.614550137046685,
  51.88334600443296,
  52.31753879094182,
  107.73708786357828,
  85.31519926752722,
  6.303081720834172,
  46.787323055736074,
  84.70630547292899,
  143.85975401030373,
  7.3537901560015015,
  73.27950666612377,
  330.54533875957173,
  18.425495912807804,
  26.045620856974775,
  70.21511869375018,
  17.651154762992643,
  7.960273596734069,
  61.14223357952824,
  39.507815622332366,
  10.770396858787146,
  21.213062720514234,
  29.63253395404742,
  27.452509748097242,
  66.60691955384827,
  128.90192917860517,
  16.755790986489913,
  46.44826570511603,
  59.315598646901286,
  79.46153529905837,
  40.25577085773924,
  108.23979642176448,
  22.516456183692274,
  5.741456677168784,
  18.811594495253676,
  21.893520529197502,
  36.54576895803149,
  33.688813863701284,
  68.95964196801158,
  6.503486900934444,
  5.147034275536632,
  30.41197120249592,
  6.605420319494377,
  6.647531825946573,
  18.874867498320626,
  52.15835174953941,
  181.7403717657994,
  24.078806475310586,
  177.16699665156847,
  25.241317152568293,
  34.29508148292601,
  53.518070480521416,
  27.579970899651215,
  32.96840703842789,
  16.896268257507106,
  40.85194899710976,
  50.26677166194595,
  19.010876155511188,
  22.129783319390395,
  46.776678077127485,
  52.56579863462705,
  132.47698821146398,
  7.62983841837847,
  9.436694721913003,
  7.802598655730049,
  65.90542570284794,
  16.181547762364193,
  22.96191546722195,
  62.73762266520167,
  17.841042814899854,
  47.374685999798096,
  81.18315144530489,
  30.074014107145796,
  20.2346524867879,
  11.887206111573043,
  26.827274118708146,
  52.43006559295528,
  31.849280747464533,
  76.82246681915748,
  89.13051914494739,
  5.6146414211399716,
  14.711733662001759,
  13.816366731903393,
  75.86374802445212,
  18.57259661362213,
  11.261821955603915,
  5.502424591714042,
  23.936023638749308,
  35.73385136482164,
  24.327824452358765,
  23.915348838703206,
  6.158301841289045,
  75.8947431426148,
  15.779341696093134,
  63.94093662562174,
  154.88474496179057,
  29.345047728999045,
  30.920439227951356,
  33.1444244211425,
  79.69991577203174,
  11.958679077069876,
  121.10756591618468,
  10.491301494977966,
  64.07244632242285,
  52.28576675345078,
  16.92258996968778,
  54.27891163910192,
  44.26415586655249,
  52.96604387235075,
  11.291465444685056,
  7.157871142123615,
  15.645953312411521,
  9.452141141480972,
  38.6407084353018,
  53.76822579035287,
  57.02623948424671,
  62.05480028849644,
  8.426892858810728,
  126.46106618438952,
  20.67317908359952,
  30.690167757589336,
  14.438998032388177,
  69.54859315438307,
  97.70424466107696,
  24.981846845716742,
  5.729735752876958,
  29.38332104303501,
  61.970485907931206,
  24.106062675946326,
  11.540788343131759,
  98.10756601677218,
  6.16055637294852,
  13.735693101726685,
  29.73198078921328,
  6.564734087550668,
  98.00102546989545,
  13.434529785946006,
  29.024543272056835,
  56.89878424571103,
  27.6709659937862,
  100.50044738367049,
  10.322235437136158,
  5.362098907062477,
  19.987488349591487,
  13.857681181322183,
  132.16324307238256,
  24.385043352380485,
  23.491906890427916,
  44.66388570348595,
  81.31535588261062,
  59.61787451841501,
  46.549146042846495,
  34.82949706502397,
  5.247245845538025,
  9.021315415517376,
  10.837506222838575,
  15.753549652327102,
  15.007480400229264,
  85.83560456214487,
  12.231654049718982,
  5.640178978991911,
  139.5932121129734,
  13.740795005930728,
  56.36472326793371,
  18.3617060468938,
  64.02033677431709,
  6.935084607671953,
  7.099546612049477,
  24.789357288707183,
  16.743350206751078,
  25.814383584344597,
  28.810754563221913,
  16.47504598038929,
  95.89410128100813,
  60.57847072199497,
  6.433997047264384,
  72.87084656078866,
  11.50553350795801,
  10.632438960422588,
  10.073797253697192,
  10.694525586975601,
  6.272487883977548,
  11.706646961265921,
  7.916893884008875,
  10.162630226257125,
  68.68974515674206,
  68.78178402411277,
  5.7439091374785605,
  42.00238064554881,
  12.232992744449751,
  18.233348392344055,
  49.40488233255257,
  45.602511401472846,
  59.80945190365659,
  15.723934447434015,
  13.010220646899924,
  45.208060073584406,
  14.094553769932183,
  10.80704366956071,
  47.66663269538455,
  73.07094504566307,
  11.238439045063636,
  8.779980466783517,
  30.06709297597874])

Energy spectrum

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

([4865392.1875,
  4866430.372746798,
  4902988.877834031,
  4931209.506342441,
  4937882.430525896,
  4989789.699282084,
  5073603.342721403,
  5079138.615780119,
  5093085.45832693,
  5104781.20245201,
  5118369.247009864,
  5124782.27535554,
  5140497.130575997,
  5146756.0449977415,
  5150816.348784119,
  5153510.32937531,
  5167977.597559522,
  5172818.656928344,
  5193554.6875,
  5199787.575826112,
  5200972.652611837,
  5203475.506271819,
  5220719.785682488,
  5231376.29162031,
  5237002.271631566,
  5245253.370920323,
  5245967.1875,
  5247632.7693357095,
  5247760.907177379,
  5251951.340233685,
  5255084.084850079,
  5255424.468436674,
  5256010.9375,
  5263303.125,
  5265394.059089717,
  5275656.25,
  5276981.663172591,
  5279809.096651305,
  5280143.868922852,
  5280217.973431971,
  5282740.7915373845,
  5290157.5582990125,
  5291563.740521363,
  5296328.34780919,
  5298732.450366807,
  5301565.493446072,
  5301655.840839303,
  5305144.112784824,
  5317597.998484296,
  5317853.277808623,
  5319240.625,
  5325339.72645843,
  5327757.3651289325,
  5327767.1875,
  5327929.6875,
  5337276.5625,
  5338517.159864835,
  5342450.874912309,
  5349233.136463841,
  5360171.875,
  5374894.392739861,
  5377357.8125,
  5391899.3522231085,
  5393501.44976538,
  5393943.444008141,
  5396660.674678239,
  5406919.580455997,
  5407604.6875,
  5434826.5625,
  5441242.1875,
  5444014.977607985,
  5444230.828413099,
  5446717.084932526,
  5446754.99661403,
  5449906.246461507,
  5449924.498939607,
  5454240.625,
  5462168.74143095,
  5476964.340195499,
  5479483.608329411,
  5506009.027312048,
  5511290.801604016,
  5511955.605892921,
  5527270.8912632475,
  5528318.75,
  5564456.068549881,
  5575449.680018144,
  5581357.50364499,
  5595592.550203723,
  5598322.042129034,
  5617149.081652434,
  5624150.0,
  5648496.875,
  5649247.355164573,
  5649866.754714486,
  5653337.088260016,
  5654041.669712528,
  5657492.1875,
  5658982.774321432,
  5677289.376096702,
  5680786.293636969,
  5684165.625,
  5686352.953112831,
  5689164.601985192,
  5695773.73125102,
  5695848.7187655205,
  5706056.673345632,
  5714988.146156091,
  5722040.125979857,
  5722330.90480796,
  5743976.335496414,
  5751294.984614641,
  5751347.838897046,
  5760979.496384642,
  5774444.461648831,
  5775406.25,
  5796264.774612336,
  5815738.628183398,
  5825473.372968249,
  5840247.504002235,
  5878079.509601273,
  5890681.25,
  5910840.968365635,
  5912360.9375,
  5936604.598761982,
  5953924.510853552,
  5973079.6875,
  5983315.625,
  5997006.25,
  6001074.010843979,
  6002660.931751104,
  6006582.404836634,
  6011543.148967736,
  6014060.780350795,
  6027293.106376829,
  6029378.316641708,
  6063359.582652451,
  6082106.761413057,
  6082450.317008664,
  6087117.098017578,
  6097646.694387269,
  6105735.9375,
  6119614.9334454425,
  6127973.469032305,
  6128865.075723356,
  6141499.462182029,
  6149806.49389333,
  6154144.799703064,
  6193469.929502385,
  6202392.1875,
  6214989.698152065,
  6220368.830306384,
  6224790.625,
  6232175.303909209,
  6240180.36032181,
  6257892.221069465,
  6281338.799058635,
  6287671.875,
  6296088.45071476,
  6300289.367994876,
  6307860.05450089,
  6310993.857144763,
  6316745.190537277,
  6319001.440458031,
  6319027.418868509,
  6325835.83252889,
  6326250.503506313,
  6326490.625,
  6330628.724554141,
  6356588.320240547,
  6360415.335313595,
  6363562.221853637,
  6371274.374273011,
  6386523.422946879,
  6394667.1875,
  6394808.378273395,
  6394951.731885219,
  6403358.2815756025,
  6421995.7453817045,
  6437765.826745795,
  6438092.1875,
  6447632.190745184,
  6447646.992472689,
  6455367.293921407,
  6459441.354706894,
  6464406.181286556,
  6472130.82115724,
  6475175.0,
  6475564.552828369,
  6475599.2993196,
  6477937.928923169,
  6535332.618258111,
  6535743.835697867,
  6547525.0,
  6555413.456796719,
  6562930.399962301,
  6603127.438495427,
  6615900.671521435,
  6635863.955488209,
  6636370.300816379,
  6642985.021600708,
  6643038.471486368,
  6645496.875,
  6677741.959400974,
  6678039.534533036,
  6678332.126183774,
  6679654.6875,
  6680502.849347608,
  6680672.247446646,
  6685620.630531663,
  6688166.228500593,
  6695758.554384216,
  6713062.118196205,
  6713070.088983797,
  6725877.498739352,
  6733785.048467853,
  6736120.832033469,
  6756957.662575316,
  6758993.164265754,
  6759862.5,
  6767173.048328394,
  6769968.1656420585,
  6772373.67911551,
  6773278.756659159,
  6774296.8185508605,
  6774812.5,
  6776609.226339851,
  6822439.211374638,
  6822690.625,
  6827356.25,
  6827370.032493986,
  6838047.811853215,
  6864329.6875,
  6892826.139903439,
  6905046.199146711,
  6935284.180576586,
  6936528.697655473,
  6936781.25,
  6942010.130280411,
  6945760.9375,
  6949218.227563352,
  6951908.091333866,
  7000270.03845738,
  7000946.124665435,
  7008157.8125,
  7010532.065798446,
  7013267.952952388,
  7013883.764522348,
  7015869.7056191275,
  7019588.923166189,
  7053725.0,
  7063185.739657797,
  7067185.264894794,
  7072989.0625,
  7076690.728302234,
  7087073.196209238,
  7089940.096461754,
  7090184.375,
  7093756.193390366,
  7095231.760088102,
  7096520.08635461,
  7104253.53146217,
  7151514.0625,
  7159301.883181753,
  7296915.296155198,
  7398592.1875,
  7417698.280867868,
  7610501.410732309,
  7613122.4067467395,
  7633831.993895526,
  7694799.539224773,
  7698179.8815623205,
  7701135.712461473,
  7701135.9375,
  7710684.47880043,
  7716047.632769352,
  7727164.0625,
  7729813.215483495,
  7739449.389936655,
  7744217.940267222,
  7749068.75,
  7825326.5625,
  7847084.851139585,
  7851373.718111622,
  7854680.821309318,
  7917642.167052413,
  7918303.125,
  7919589.165673966,
  7921003.220929531,
  7921179.7965270085,
  7921527.851020623,
  7921530.550743104,
  7927358.943530014,
  7939742.105494164,
  7954414.0625,
  8007575.436703201,
  8007686.043733773,
  8022243.8229530845,
  8022659.496952706,
  8025723.957001001,
  8045795.318675183,
  8048839.353779782,
  8117460.786861867,
  8119758.507532099,
  8119759.972393939,
  8120473.281012251,
  8138021.6670857165,
  8166320.73808537,
  8201688.426360736,
  8206147.3120294465,
  8208581.975229668,
  8216937.580248538,
  8223488.295775776,
  8240317.1875,
  8294933.041968273,
  8320040.635791321,
  8326085.875880561,
  8328465.651843779,
  8369207.567085598,
  8501682.92761445,
  8501692.27164109,
  8501721.226743998,
  8518421.617780104,
  8549568.690662805,
  8552057.838524133,
  8553351.275278158,
  8562338.149069987,
  8568493.75,
  8597172.078720665,
  8608918.75,
  8656933.469497064,
  8659381.923756963,
  8681962.5,
  8745687.083034104,
  8782968.75,
  8783012.973672677,
  8786482.8125,
  8797340.17574899,
  8802787.551015569,
  8833018.512568837,
  8834316.016556965,
  8836015.625,
  8836911.6536696,
  8838725.0,
  8843837.602300921,
  8849181.386757499,
  8853117.845908893,
  8853892.326921102,
  8888759.375,
  8897020.3125,
  8910762.5,
  8911898.208020471,
  8916158.068736866,
  8917901.5625,
  8930655.663681623,
  8936133.845905786,
  8936922.108446503,
  8938055.517940393,
  8938203.125,
  8938923.228131622,
  8939000.355089176,
  8941153.067436485,
  8941657.231177805,
  8941778.020609286,
  8945236.18662382,
  8949599.838124905,
  8952072.544086674,
  8954128.70838222,
  8954268.40539949,
  8957745.79054372,
  8963076.382611705,
  8963167.81675824,
  8964237.5,
  8964747.823468974,
  8966033.80114346,
  8967071.0880002,
  8967099.415342307,
  8967130.264751952,
  8967186.834550261,
  8967311.329356162,
  8967542.1875,
  8968271.623086121,
  8980504.366962962,
  8984111.785309391,
  8987104.404950673,
  8991674.51527097,
  9054337.5,
  9106188.932355177,
  9125989.137378976,
  9130365.234000752,
  9132952.931051036,
  9150770.510395326,
  9155765.642201256,
  9168907.542526789,
  9178734.982517647,
  9233720.775440203,
  9277111.54348997,
  9472780.894757938,
  9900853.125,
  9914369.428507341,
  9915280.378040625,
  9925249.240184313,
  9928495.228508929,
  9936715.402180506,
  9937158.216579778,
  9941225.07845454,
  9944835.293526413,
  9945089.414011562,
  9945108.780579712,
  9947621.875,
  9948051.21522697,
  9951507.778116597,
  9956662.717577567,
  9963130.908218045,
  9963538.264386438,
  9963541.891482808,
  9967531.291752014,
  9967820.653793678,
  9969878.916244522,
  9970260.21780178,
  9971392.061634675,
  9973409.375,
  9974470.880837968,
  9977762.708106669,
  9980772.483151233,
  9985644.409075178,
  9986233.12707635,
  9988371.875,
  9995728.125,
  10001225.191173948,
  10006386.880484438,
  10010594.159571325,
  10013183.349512633,
  10014871.659683162,
  10015761.451107955,
  10021789.0625,
  10049328.194706488,
  10053358.641931511,
  10066326.120938992],
 [41.92171833418734,
  18.406778570380563,
  17.182598204863623,
  16.04072393902368,
  23.366057742176505,
  111.73035750812565,
  32.17263621088415,
  13.474375444383824,
  13.383450326542999,
  27.61778155813478,
  6.866960046466185,
  5.507567019094055,
  16.014175336720857,
  22.542489511418584,
  7.887822436325474,
  34.24940894886531,
  21.223665630181195,
  12.521139068320783,
  74.32773542709742,
  11.248333040615185,
  35.725563865620295,
  102.52744007144821,
  89.77889144838295,
  7.508117754744114,
  18.935935858196718,
  26.574759864928197,
  57.314746968540994,
  27.69858528190561,
  61.12364519878206,
  32.51491457194023,
  112.73840521746234,
  16.095115607142255,
  30.055594216079108,
  60.18235688516705,
  11.64783523824972,
  34.29117961954461,
  26.543888820978193,
  20.35007143511837,
  17.048044030941227,
  5.758671659712053,
  154.964331301543,
  59.17653832191581,
  13.841563984915918,
  9.527185249402326,
  22.988247212687696,
  26.170595781605172,
  9.007845191825291,
  81.60102991448929,
  96.72964398215204,
  62.29544449035191,
  63.12351423519161,
  42.581263059289,
  22.03467426212426,
  29.142767456464163,
  54.01520478968915,
  32.3446403650468,
  63.35122254770762,
  82.552795834706,
  16.33902645932599,
  56.66032586978959,
  42.6712759946148,
  55.47676703708981,
  52.184099426065124,
  14.5285498675742,
  23.674504510823443,
  5.388520763729763,
  19.453515682549053,
  38.70423609105755,
  51.03558424211106,
  30.29836203756138,
  40.99372581897157,
  23.876864381750938,
  14.564469163866981,
  16.266961002626093,
  14.928039472812234,
  21.482528768601814,
  91.87485407327918,
  22.286208876610225,
  37.6336973101285,
  12.650575376484758,
  156.98218384316613,
  9.485333663940402,
  20.16669359586452,
  162.7375945311784,
  43.91404499038363,
  14.026553857262552,
  29.715968347819018,
  64.40919262298496,
  127.8340116175549,
  66.24357607086449,
  7.854107644685509,
  43.8442993165237,
  90.08803383478521,
  5.73801395639558,
  24.078380117812017,
  22.42116708322602,
  13.277592907349314,
  39.832070479177574,
  27.79117039213128,
  72.32562738142042,
  49.30422006886399,
  207.79068931228366,
  18.323258252264726,
  12.090503156031742,
  8.448276550897713,
  5.046315850607294,
  72.31670914268936,
  6.899747738586512,
  24.11331186132422,
  26.295001169551043,
  24.73021146824968,
  221.4464167445425,
  7.993511107139844,
  20.196898697076485,
  8.702537247997737,
  32.423071887914034,
  78.61875089229471,
  90.8565925859285,
  7.104590919145231,
  17.99112147107802,
  85.31308960311398,
  36.72267188223688,
  35.65344189825217,
  69.45610045536584,
  26.88497676711342,
  10.440504649978429,
  68.29918668074761,
  155.24336803635688,
  71.9949963402927,
  25.694627540462186,
  22.36170146782888,
  23.97751578058194,
  32.65473515789066,
  53.86486148178036,
  30.93897808313987,
  17.06299458372032,
  14.310452392930296,
  152.61754044891873,
  11.662209833367863,
  18.733156659837153,
  82.89903039663994,
  91.8872044592244,
  17.205809415186046,
  66.65256591826027,
  22.53953929023367,
  21.123962669226557,
  40.00982907570861,
  10.410785603410375,
  25.768047455859197,
  103.30971703624303,
  34.03358804937414,
  11.255594593336562,
  80.02301407409479,
  18.054507105333013,
  44.41669704364794,
  20.110729860372608,
  62.62807676544802,
  95.52481211777905,
  20.44603768204785,
  40.095057170367056,
  65.83613246124716,
  106.67755094758725,
  16.513157022837653,
  40.240195836001504,
  17.6362914187126,
  5.528007558688627,
  6.332980053958412,
  71.11715682504537,
  7.8355036025151525,
  44.99224835778598,
  41.09603295048318,
  14.324079872713874,
  85.17457039982077,
  39.56788363755811,
  54.54017601491058,
  213.41744226137249,
  76.98631210072182,
  115.94930936287884,
  53.02697181438315,
  119.61291458060569,
  68.63845579837657,
  20.358829046692907,
  76.17867945682896,
  57.84908793124893,
  28.259205181569705,
  61.05173265173711,
  47.52541668551477,
  50.98278050150255,
  17.21021739233419,
  13.350754869949611,
  47.93421246628941,
  97.68083294676879,
  22.042131247716846,
  47.86036492203617,
  20.499707251894527,
  7.882989602247322,
  10.201067856043407,
  18.31528909845606,
  34.024119705998935,
  28.295759354659463,
  8.66091351814023,
  5.4369587574684255,
  31.45163209026017,
  9.1457029356276,
  17.108271175313337,
  19.816160940109455,
  68.87824639675871,
  10.722548843222299,
  6.315429758225304,
  19.648170141427272,
  101.5891540340275,
  5.222744089946213,
  67.88285399545867,
  9.149297886930816,
  6.966014626835337,
  141.2117462592151,
  9.813057895638435,
  25.660286854180207,
  30.263023988025257,
  44.956186737419365,
  10.56941302518406,
  62.18650836042743,
  6.298215119365958,
  12.282985189294758,
  27.455416535526005,
  58.224677592197324,
  38.990191625499115,
  15.541027777452356,
  43.193792243072224,
  72.91008299765167,
  24.009541643040983,
  93.93204348292733,
  40.79289963785829,
  72.89029492208391,
  18.59355293415175,
  29.1037590842152,
  8.614550137046685,
  51.88334600443296,
  52.31753879094182,
  107.73708786357828,
  85.31519926752722,
  6.303081720834172,
  46.787323055736074,
  84.70630547292899,
  143.85975401030373,
  7.3537901560015015,
  73.27950666612377,
  330.54533875957173,
  18.425495912807804,
  26.045620856974775,
  70.21511869375018,
  17.651154762992643,
  7.960273596734069,
  61.14223357952824,
  39.507815622332366,
  10.770396858787146,
  21.213062720514234,
  29.63253395404742,
  27.452509748097242,
  66.60691955384827,
  128.90192917860517,
  16.755790986489913,
  46.44826570511603,
  59.315598646901286,
  79.46153529905837,
  40.25577085773924,
  108.23979642176448,
  22.516456183692274,
  5.741456677168784,
  18.811594495253676,
  21.893520529197502,
  36.54576895803149,
  33.688813863701284,
  68.95964196801158,
  6.503486900934444,
  5.147034275536632,
  30.41197120249592,
  6.605420319494377,
  6.647531825946573,
  18.874867498320626,
  52.15835174953941,
  181.7403717657994,
  24.078806475310586,
  177.16699665156847,
  25.241317152568293,
  34.29508148292601,
  53.518070480521416,
  27.579970899651215,
  32.96840703842789,
  16.896268257507106,
  40.85194899710976,
  50.26677166194595,
  19.010876155511188,
  22.129783319390395,
  46.776678077127485,
  52.56579863462705,
  132.47698821146398,
  7.62983841837847,
  9.436694721913003,
  7.802598655730049,
  65.90542570284794,
  16.181547762364193,
  22.96191546722195,
  62.73762266520167,
  17.841042814899854,
  47.374685999798096,
  81.18315144530489,
  30.074014107145796,
  20.2346524867879,
  11.887206111573043,
  26.827274118708146,
  52.43006559295528,
  31.849280747464533,
  76.82246681915748,
  89.13051914494739,
  5.6146414211399716,
  14.711733662001759,
  13.816366731903393,
  75.86374802445212,
  18.57259661362213,
  11.261821955603915,
  5.502424591714042,
  23.936023638749308,
  35.73385136482164,
  24.327824452358765,
  23.915348838703206,
  6.158301841289045,
  75.8947431426148,
  15.779341696093134,
  63.94093662562174,
  154.88474496179057,
  29.345047728999045,
  30.920439227951356,
  33.1444244211425,
  79.69991577203174,
  11.958679077069876,
  121.10756591618468,
  10.491301494977966,
  64.07244632242285,
  52.28576675345078,
  16.92258996968778,
  54.27891163910192,
  44.26415586655249,
  52.96604387235075,
  11.291465444685056,
  7.157871142123615,
  15.645953312411521,
  9.452141141480972,
  38.6407084353018,
  53.76822579035287,
  57.02623948424671,
  62.05480028849644,
  8.426892858810728,
  126.46106618438952,
  20.67317908359952,
  30.690167757589336,
  14.438998032388177,
  69.54859315438307,
  97.70424466107696,
  24.981846845716742,
  5.729735752876958,
  29.38332104303501,
  61.970485907931206,
  24.106062675946326,
  11.540788343131759,
  98.10756601677218,
  6.16055637294852,
  13.735693101726685,
  29.73198078921328,
  6.564734087550668,
  98.00102546989545,
  13.434529785946006,
  29.024543272056835,
  56.89878424571103,
  27.6709659937862,
  100.50044738367049,
  10.322235437136158,
  5.362098907062477,
  19.987488349591487,
  13.857681181322183,
  132.16324307238256,
  24.385043352380485,
  23.491906890427916,
  44.66388570348595,
  81.31535588261062,
  59.61787451841501,
  46.549146042846495,
  34.82949706502397,
  5.247245845538025,
  9.021315415517376,
  10.837506222838575,
  15.753549652327102,
  15.007480400229264,
  85.83560456214487,
  12.231654049718982,
  5.640178978991911,
  139.5932121129734,
  13.740795005930728,
  56.36472326793371,
  18.3617060468938,
  64.02033677431709,
  6.935084607671953,
  7.099546612049477,
  24.789357288707183,
  16.743350206751078,
  25.814383584344597,
  28.810754563221913,
  16.47504598038929,
  95.89410128100813,
  60.57847072199497,
  6.433997047264384,
  72.87084656078866,
  11.50553350795801,
  10.632438960422588,
  10.073797253697192,
  10.694525586975601,
  6.272487883977548,
  11.706646961265921,
  7.916893884008875,
  10.162630226257125,
  68.68974515674206,
  68.78178402411277,
  5.7439091374785605,
  42.00238064554881,
  12.232992744449751,
  18.233348392344055,
  49.40488233255257,
  45.602511401472846,
  59.80945190365659,
  15.723934447434015,
  13.010220646899924,
  45.208060073584406,
  14.094553769932183,
  10.80704366956071,
  47.66663269538455,
  73.07094504566307,
  11.238439045063636,
  8.779980466783517,
  30.06709297597874])

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

([4865392.1875,
  4866430.372746798,
  4902988.877834031,
  4931209.506342441,
  4937882.430525896,
  4989789.699282084,
  5073603.342721403,
  5079138.615780119,
  5093085.45832693,
  5104781.20245201,
  5118369.247009864,
  5124782.27535554,
  5140497.130575997,
  5146756.0449977415,
  5150816.348784119,
  5153510.32937531,
  5167977.597559522,
  5172818.656928344,
  5193554.6875,
  5199787.575826112,
  5200972.652611837,
  5203475.506271819,
  5220719.785682488,
  5231376.29162031,
  5237002.271631566,
  5245253.370920323,
  5245967.1875,
  5247632.7693357095,
  5247760.907177379,
  5251951.340233685,
  5255084.084850079,
  5255424.468436674,
  5256010.9375,
  5263303.125,
  5265394.059089717,
  5275656.25,
  5276981.663172591,
  5279809.096651305,
  5280143.868922852,
  5280217.973431971,
  5282740.7915373845,
  5290157.5582990125,
  5291563.740521363,
  5296328.34780919,
  5298732.450366807,
  5301565.493446072,
  5301655.840839303,
  5305144.112784824,
  5317597.998484296,
  5317853.277808623,
  5319240.625,
  5325339.72645843,
  5327757.3651289325,
  5327767.1875,
  5327929.6875,
  5337276.5625,
  5338517.159864835,
  5342450.874912309,
  5349233.136463841,
  5360171.875,
  5374894.392739861,
  5377357.8125,
  5391899.3522231085,
  5393501.44976538,
  5393943.444008141,
  5396660.674678239,
  5406919.580455997,
  5407604.6875,
  5434826.5625,
  5441242.1875,
  5444014.977607985,
  5444230.828413099,
  5446717.084932526,
  5446754.99661403,
  5449906.246461507,
  5449924.498939607,
  5454240.625,
  5462168.74143095,
  5476964.340195499,
  5479483.608329411,
  5506009.027312048,
  5511290.801604016,
  5511955.605892921,
  5527270.8912632475,
  5528318.75,
  5564456.068549881,
  5575449.680018144,
  5581357.50364499,
  5595592.550203723,
  5598322.042129034,
  5617149.081652434,
  5624150.0,
  5648496.875,
  5649247.355164573,
  5649866.754714486,
  5653337.088260016,
  5654041.669712528,
  5657492.1875,
  5658982.774321432,
  5677289.376096702,
  5680786.293636969,
  5684165.625,
  5686352.953112831,
  5689164.601985192,
  5695773.73125102,
  5695848.7187655205,
  5706056.673345632,
  5714988.146156091,
  5722040.125979857,
  5722330.90480796,
  5743976.335496414,
  5751294.984614641,
  5751347.838897046,
  5760979.496384642,
  5774444.461648831,
  5775406.25,
  5796264.774612336,
  5815738.628183398,
  5825473.372968249,
  5840247.504002235,
  5878079.509601273,
  5890681.25,
  5910840.968365635,
  5912360.9375,
  5936604.598761982,
  5953924.510853552,
  5973079.6875,
  5983315.625,
  5997006.25,
  6001074.010843979,
  6002660.931751104,
  6006582.404836634,
  6011543.148967736,
  6014060.780350795,
  6027293.106376829,
  6029378.316641708,
  6063359.582652451,
  6082106.761413057,
  6082450.317008664,
  6087117.098017578,
  6097646.694387269,
  6105735.9375,
  6119614.9334454425,
  6127973.469032305,
  6128865.075723356,
  6141499.462182029,
  6149806.49389333,
  6154144.799703064,
  6193469.929502385,
  6202392.1875,
  6214989.698152065,
  6220368.830306384,
  6224790.625,
  6232175.303909209,
  6240180.36032181,
  6257892.221069465,
  6281338.799058635,
  6287671.875,
  6296088.45071476,
  6300289.367994876,
  6307860.05450089,
  6310993.857144763,
  6316745.190537277,
  6319001.440458031,
  6319027.418868509,
  6325835.83252889,
  6326250.503506313,
  6326490.625,
  6330628.724554141,
  6356588.320240547,
  6360415.335313595,
  6363562.221853637,
  6371274.374273011,
  6386523.422946879,
  6394667.1875,
  6394808.378273395,
  6394951.731885219,
  6403358.2815756025,
  6421995.7453817045,
  6437765.826745795,
  6438092.1875,
  6447632.190745184,
  6447646.992472689,
  6455367.293921407,
  6459441.354706894,
  6464406.181286556,
  6472130.82115724,
  6475175.0,
  6475564.552828369,
  6475599.2993196,
  6477937.928923169,
  6535332.618258111,
  6535743.835697867,
  6547525.0,
  6555413.456796719,
  6562930.399962301,
  6603127.438495427,
  6615900.671521435,
  6635863.955488209,
  6636370.300816379,
  6642985.021600708,
  6643038.471486368,
  6645496.875,
  6677741.959400974,
  6678039.534533036,
  6678332.126183774,
  6679654.6875,
  6680502.849347608,
  6680672.247446646,
  6685620.630531663,
  6688166.228500593,
  6695758.554384216,
  6713062.118196205,
  6713070.088983797,
  6725877.498739352,
  6733785.048467853,
  6736120.832033469,
  6756957.662575316,
  6758993.164265754,
  6759862.5,
  6767173.048328394,
  6769968.1656420585,
  6772373.67911551,
  6773278.756659159,
  6774296.8185508605,
  6774812.5,
  6776609.226339851,
  6822439.211374638,
  6822690.625,
  6827356.25,
  6827370.032493986,
  6838047.811853215,
  6864329.6875,
  6892826.139903439,
  6905046.199146711,
  6935284.180576586,
  6936528.697655473,
  6936781.25,
  6942010.130280411,
  6945760.9375,
  6949218.227563352,
  6951908.091333866,
  7000270.03845738,
  7000946.124665435,
  7008157.8125,
  7010532.065798446,
  7013267.952952388,
  7013883.764522348,
  7015869.7056191275,
  7019588.923166189,
  7053725.0,
  7063185.739657797,
  7067185.264894794,
  7072989.0625,
  7076690.728302234,
  7087073.196209238,
  7089940.096461754,
  7090184.375,
  7093756.193390366,
  7095231.760088102,
  7096520.08635461,
  7104253.53146217,
  7151514.0625,
  7159301.883181753,
  7296915.296155198,
  7398592.1875,
  7417698.280867868,
  7610501.410732309,
  7613122.4067467395,
  7633831.993895526,
  7694799.539224773,
  7698179.8815623205,
  7701135.712461473,
  7701135.9375,
  7710684.47880043,
  7716047.632769352,
  7727164.0625,
  7729813.215483495,
  7739449.389936655,
  7744217.940267222,
  7749068.75,
  7825326.5625,
  7847084.851139585,
  7851373.718111622,
  7854680.821309318,
  7917642.167052413,
  7918303.125,
  7919589.165673966,
  7921003.220929531,
  7921179.7965270085,
  7921527.851020623,
  7921530.550743104,
  7927358.943530014,
  7939742.105494164,
  7954414.0625,
  8007575.436703201,
  8007686.043733773,
  8022243.8229530845,
  8022659.496952706,
  8025723.957001001,
  8045795.318675183,
  8048839.353779782,
  8117460.786861867,
  8119758.507532099,
  8119759.972393939,
  8120473.281012251,
  8138021.6670857165,
  8166320.73808537,
  8201688.426360736,
  8206147.3120294465,
  8208581.975229668,
  8216937.580248538,
  8223488.295775776,
  8240317.1875,
  8294933.041968273,
  8320040.635791321,
  8326085.875880561,
  8328465.651843779,
  8369207.567085598,
  8501682.92761445,
  8501692.27164109,
  8501721.226743998,
  8518421.617780104,
  8549568.690662805,
  8552057.838524133,
  8553351.275278158,
  8562338.149069987,
  8568493.75,
  8597172.078720665,
  8608918.75,
  8656933.469497064,
  8659381.923756963,
  8681962.5,
  8745687.083034104,
  8782968.75,
  8783012.973672677,
  8786482.8125,
  8797340.17574899,
  8802787.551015569,
  8833018.512568837,
  8834316.016556965,
  8836015.625,
  8836911.6536696,
  8838725.0,
  8843837.602300921,
  8849181.386757499,
  8853117.845908893,
  8853892.326921102,
  8888759.375,
  8897020.3125,
  8910762.5,
  8911898.208020471,
  8916158.068736866,
  8917901.5625,
  8930655.663681623,
  8936133.845905786,
  8936922.108446503,
  8938055.517940393,
  8938203.125,
  8938923.228131622,
  8939000.355089176,
  8941153.067436485,
  8941657.231177805,
  8941778.020609286,
  8945236.18662382,
  8949599.838124905,
  8952072.544086674,
  8954128.70838222,
  8954268.40539949,
  8957745.79054372,
  8963076.382611705,
  8963167.81675824,
  8964237.5,
  8964747.823468974,
  8966033.80114346,
  8967071.0880002,
  8967099.415342307,
  8967130.264751952,
  8967186.834550261,
  8967311.329356162,
  8967542.1875,
  8968271.623086121,
  8980504.366962962,
  8984111.785309391,
  8987104.404950673,
  8991674.51527097,
  9054337.5,
  9106188.932355177,
  9125989.137378976,
  9130365.234000752,
  9132952.931051036,
  9150770.510395326,
  9155765.642201256,
  9168907.542526789,
  9178734.982517647,
  9233720.775440203,
  9277111.54348997,
  9472780.894757938,
  9900853.125,
  9914369.428507341,
  9915280.378040625,
  9925249.240184313,
  9928495.228508929,
  9936715.402180506,
  9937158.216579778,
  9941225.07845454,
  9944835.293526413,
  9945089.414011562,
  9945108.780579712,
  9947621.875,
  9948051.21522697,
  9951507.778116597,
  9956662.717577567,
  9963130.908218045,
  9963538.264386438,
  9963541.891482808,
  9967531.291752014,
  9967820.653793678,
  9969878.916244522,
  9970260.21780178,
  9971392.061634675,
  9973409.375,
  9974470.880837968,
  9977762.708106669,
  9980772.483151233,
  9985644.409075178,
  9986233.12707635,
  9988371.875,
  9995728.125,
  10001225.191173948,
  10006386.880484438,
  10010594.159571325,
  10013183.349512633,
  10014871.659683162,
  10015761.451107955,
  10021789.0625,
  10049328.194706488,
  10053358.641931511,
  10066326.120938992],
 [41.92171833418734,
  18.406778570380563,
  17.182598204863623,
  16.04072393902368,
  23.366057742176505,
  111.73035750812565,
  32.17263621088415,
  13.474375444383824,
  13.383450326542999,
  27.61778155813478,
  6.866960046466185,
  5.507567019094055,
  16.014175336720857,
  22.542489511418584,
  7.887822436325474,
  34.24940894886531,
  21.223665630181195,
  12.521139068320783,
  74.32773542709742,
  11.248333040615185,
  35.725563865620295,
  102.52744007144821,
  89.77889144838295,
  7.508117754744114,
  18.935935858196718,
  26.574759864928197,
  57.314746968540994,
  27.69858528190561,
  61.12364519878206,
  32.51491457194023,
  112.73840521746234,
  16.095115607142255,
  30.055594216079108,
  60.18235688516705,
  11.64783523824972,
  34.29117961954461,
  26.543888820978193,
  20.35007143511837,
  17.048044030941227,
  5.758671659712053,
  154.964331301543,
  59.17653832191581,
  13.841563984915918,
  9.527185249402326,
  22.988247212687696,
  26.170595781605172,
  9.007845191825291,
  81.60102991448929,
  96.72964398215204,
  62.29544449035191,
  63.12351423519161,
  42.581263059289,
  22.03467426212426,
  29.142767456464163,
  54.01520478968915,
  32.3446403650468,
  63.35122254770762,
  82.552795834706,
  16.33902645932599,
  56.66032586978959,
  42.6712759946148,
  55.47676703708981,
  52.184099426065124,
  14.5285498675742,
  23.674504510823443,
  5.388520763729763,
  19.453515682549053,
  38.70423609105755,
  51.03558424211106,
  30.29836203756138,
  40.99372581897157,
  23.876864381750938,
  14.564469163866981,
  16.266961002626093,
  14.928039472812234,
  21.482528768601814,
  91.87485407327918,
  22.286208876610225,
  37.6336973101285,
  12.650575376484758,
  156.98218384316613,
  9.485333663940402,
  20.16669359586452,
  162.7375945311784,
  43.91404499038363,
  14.026553857262552,
  29.715968347819018,
  64.40919262298496,
  127.8340116175549,
  66.24357607086449,
  7.854107644685509,
  43.8442993165237,
  90.08803383478521,
  5.73801395639558,
  24.078380117812017,
  22.42116708322602,
  13.277592907349314,
  39.832070479177574,
  27.79117039213128,
  72.32562738142042,
  49.30422006886399,
  207.79068931228366,
  18.323258252264726,
  12.090503156031742,
  8.448276550897713,
  5.046315850607294,
  72.31670914268936,
  6.899747738586512,
  24.11331186132422,
  26.295001169551043,
  24.73021146824968,
  221.4464167445425,
  7.993511107139844,
  20.196898697076485,
  8.702537247997737,
  32.423071887914034,
  78.61875089229471,
  90.8565925859285,
  7.104590919145231,
  17.99112147107802,
  85.31308960311398,
  36.72267188223688,
  35.65344189825217,
  69.45610045536584,
  26.88497676711342,
  10.440504649978429,
  68.29918668074761,
  155.24336803635688,
  71.9949963402927,
  25.694627540462186,
  22.36170146782888,
  23.97751578058194,
  32.65473515789066,
  53.86486148178036,
  30.93897808313987,
  17.06299458372032,
  14.310452392930296,
  152.61754044891873,
  11.662209833367863,
  18.733156659837153,
  82.89903039663994,
  91.8872044592244,
  17.205809415186046,
  66.65256591826027,
  22.53953929023367,
  21.123962669226557,
  40.00982907570861,
  10.410785603410375,
  25.768047455859197,
  103.30971703624303,
  34.03358804937414,
  11.255594593336562,
  80.02301407409479,
  18.054507105333013,
  44.41669704364794,
  20.110729860372608,
  62.62807676544802,
  95.52481211777905,
  20.44603768204785,
  40.095057170367056,
  65.83613246124716,
  106.67755094758725,
  16.513157022837653,
  40.240195836001504,
  17.6362914187126,
  5.528007558688627,
  6.332980053958412,
  71.11715682504537,
  7.8355036025151525,
  44.99224835778598,
  41.09603295048318,
  14.324079872713874,
  85.17457039982077,
  39.56788363755811,
  54.54017601491058,
  213.41744226137249,
  76.98631210072182,
  115.94930936287884,
  53.02697181438315,
  119.61291458060569,
  68.63845579837657,
  20.358829046692907,
  76.17867945682896,
  57.84908793124893,
  28.259205181569705,
  61.05173265173711,
  47.52541668551477,
  50.98278050150255,
  17.21021739233419,
  13.350754869949611,
  47.93421246628941,
  97.68083294676879,
  22.042131247716846,
  47.86036492203617,
  20.499707251894527,
  7.882989602247322,
  10.201067856043407,
  18.31528909845606,
  34.024119705998935,
  28.295759354659463,
  8.66091351814023,
  5.4369587574684255,
  31.45163209026017,
  9.1457029356276,
  17.108271175313337,
  19.816160940109455,
  68.87824639675871,
  10.722548843222299,
  6.315429758225304,
  19.648170141427272,
  101.5891540340275,
  5.222744089946213,
  67.88285399545867,
  9.149297886930816,
  6.966014626835337,
  141.2117462592151,
  9.813057895638435,
  25.660286854180207,
  30.263023988025257,
  44.956186737419365,
  10.56941302518406,
  62.18650836042743,
  6.298215119365958,
  12.282985189294758,
  27.455416535526005,
  58.224677592197324,
  38.990191625499115,
  15.541027777452356,
  43.193792243072224,
  72.91008299765167,
  24.009541643040983,
  93.93204348292733,
  40.79289963785829,
  72.89029492208391,
  18.59355293415175,
  29.1037590842152,
  8.614550137046685,
  51.88334600443296,
  52.31753879094182,
  107.73708786357828,
  85.31519926752722,
  6.303081720834172,
  46.787323055736074,
  84.70630547292899,
  143.85975401030373,
  7.3537901560015015,
  73.27950666612377,
  330.54533875957173,
  18.425495912807804,
  26.045620856974775,
  70.21511869375018,
  17.651154762992643,
  7.960273596734069,
  61.14223357952824,
  39.507815622332366,
  10.770396858787146,
  21.213062720514234,
  29.63253395404742,
  27.452509748097242,
  66.60691955384827,
  128.90192917860517,
  16.755790986489913,
  46.44826570511603,
  59.315598646901286,
  79.46153529905837,
  40.25577085773924,
  108.23979642176448,
  22.516456183692274,
  5.741456677168784,
  18.811594495253676,
  21.893520529197502,
  36.54576895803149,
  33.688813863701284,
  68.95964196801158,
  6.503486900934444,
  5.147034275536632,
  30.41197120249592,
  6.605420319494377,
  6.647531825946573,
  18.874867498320626,
  52.15835174953941,
  181.7403717657994,
  24.078806475310586,
  177.16699665156847,
  25.241317152568293,
  34.29508148292601,
  53.518070480521416,
  27.579970899651215,
  32.96840703842789,
  16.896268257507106,
  40.85194899710976,
  50.26677166194595,
  19.010876155511188,
  22.129783319390395,
  46.776678077127485,
  52.56579863462705,
  132.47698821146398,
  7.62983841837847,
  9.436694721913003,
  7.802598655730049,
  65.90542570284794,
  16.181547762364193,
  22.96191546722195,
  62.73762266520167,
  17.841042814899854,
  47.374685999798096,
  81.18315144530489,
  30.074014107145796,
  20.2346524867879,
  11.887206111573043,
  26.827274118708146,
  52.43006559295528,
  31.849280747464533,
  76.82246681915748,
  89.13051914494739,
  5.6146414211399716,
  14.711733662001759,
  13.816366731903393,
  75.86374802445212,
  18.57259661362213,
  11.261821955603915,
  5.502424591714042,
  23.936023638749308,
  35.73385136482164,
  24.327824452358765,
  23.915348838703206,
  6.158301841289045,
  75.8947431426148,
  15.779341696093134,
  63.94093662562174,
  154.88474496179057,
  29.345047728999045,
  30.920439227951356,
  33.1444244211425,
  79.69991577203174,
  11.958679077069876,
  121.10756591618468,
  10.491301494977966,
  64.07244632242285,
  52.28576675345078,
  16.92258996968778,
  54.27891163910192,
  44.26415586655249,
  52.96604387235075,
  11.291465444685056,
  7.157871142123615,
  15.645953312411521,
  9.452141141480972,
  38.6407084353018,
  53.76822579035287,
  57.02623948424671,
  62.05480028849644,
  8.426892858810728,
  126.46106618438952,
  20.67317908359952,
  30.690167757589336,
  14.438998032388177,
  69.54859315438307,
  97.70424466107696,
  24.981846845716742,
  5.729735752876958,
  29.38332104303501,
  61.970485907931206,
  24.106062675946326,
  11.540788343131759,
  98.10756601677218,
  6.16055637294852,
  13.735693101726685,
  29.73198078921328,
  6.564734087550668,
  98.00102546989545,
  13.434529785946006,
  29.024543272056835,
  56.89878424571103,
  27.6709659937862,
  100.50044738367049,
  10.322235437136158,
  5.362098907062477,
  19.987488349591487,
  13.857681181322183,
  132.16324307238256,
  24.385043352380485,
  23.491906890427916,
  44.66388570348595,
  81.31535588261062,
  59.61787451841501,
  46.549146042846495,
  34.82949706502397,
  5.247245845538025,
  9.021315415517376,
  10.837506222838575,
  15.753549652327102,
  15.007480400229264,
  85.83560456214487,
  12.231654049718982,
  5.640178978991911,
  139.5932121129734,
  13.740795005930728,
  56.36472326793371,
  18.3617060468938,
  64.02033677431709,
  6.935084607671953,
  7.099546612049477,
  24.789357288707183,
  16.743350206751078,
  25.814383584344597,
  28.810754563221913,
  16.47504598038929,
  95.89410128100813,
  60.57847072199497,
  6.433997047264384,
  72.87084656078866,
  11.50553350795801,
  10.632438960422588,
  10.073797253697192,
  10.694525586975601,
  6.272487883977548,
  11.706646961265921,
  7.916893884008875,
  10.162630226257125,
  68.68974515674206,
  68.78178402411277,
  5.7439091374785605,
  42.00238064554881,
  12.232992744449751,
  18.233348392344055,
  49.40488233255257,
  45.602511401472846,
  59.80945190365659,
  15.723934447434015,
  13.010220646899924,
  45.208060073584406,
  14.094553769932183,
  10.80704366956071,
  47.66663269538455,
  73.07094504566307,
  11.238439045063636,
  8.779980466783517,
  30.06709297597874])

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)