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

([3497714.0379049713,
  3550652.8231579107,
  3597623.4375,
  3646509.905872267,
  3647718.565175181,
  3647873.893272848,
  3662370.3125,
  3712370.6367102996,
  3753201.500114644,
  3759410.5693005226,
  3778331.047855071,
  3778551.501607855,
  3885142.931882811,
  3886921.875,
  3895424.4876268203,
  3897929.6875,
  3999016.6810882944,
  4022796.8401280846,
  4026175.853823638,
  4104931.2364454693,
  4118038.843785412,
  4143906.25,
  4157779.1259316937,
  4237011.1149124345,
  4243728.092050614,
  4250592.1875,
  4272939.0625,
  4285544.280725929,
  4291148.4375,
  4301803.125,
  4318915.512789876,
  4319134.699317267,
  4341363.051125232,
  4358749.144718171,
  4359476.5625,
  4374908.634141004,
  4384786.82718285,
  4573999.094773447,
  4609680.879308154,
  4671587.5,
  4689415.012096634,
  4701065.625,
  4706709.375,
  4723603.372659559,
  4728463.560068321,
  4728487.877198473,
  4729141.411121211,
  4740712.5,
  4748134.236047245,
  4751720.753708192,
  4759046.875,
  4759460.648647059,
  4760292.1875,
  4767408.836273323,
  4778680.46903719,
  4789343.75,
  4812336.541152843,
  4815291.898985861,
  4843171.430583145,
  4845721.736798774,
  4851876.332625404,
  4875370.3125,
  4882090.625,
  4898103.18465866,
  4899745.266350453,
  4912727.28557971,
  4915565.625,
  4919764.18165474,
  4925560.9375,
  4933957.817040829,
  4962202.643552519,
  5007348.657002478,
  5054886.353251419,
  5154244.291191751,
  5552987.794967575,
  5576033.3656686945,
  5667795.779365314,
  5707228.125,
  5711295.008695379,
  5713535.943563967,
  5748974.027173245,
  5759927.501234855,
  5772332.820944044,
  5773675.0,
  5780288.018863097,
  5781492.843685747,
  5789097.569837585,
  5812277.270626946,
  5855300.946097331,
  5881275.470204075,
  5898355.417147607,
  5936006.25,
  5970291.935350806,
  5981324.012971827,
  6002224.595829856,
  6034448.42842006,
  6060244.896839247,
  6074544.324333707,
  6089863.381991046,
  6103499.3681991175,
  6147611.6818406815,
  6151015.625,
  6174600.64758485,
  6176357.690619808,
  6185452.579336756,
  6192746.875,
  6200629.790214532,
  6216687.5,
  6236551.427525158,
  6254308.061782735,
  6285741.730673833,
  6289832.533062202,
  6323686.032446726,
  6325717.236582025,
  6332524.873985,
  6336617.234567944,
  6336634.375,
  6365217.708124221,
  6381674.0505271135,
  6386181.379605349,
  6403343.095038348,
  6407563.196922803,
  6410983.086725494,
  6413475.364538637,
  6429345.3125,
  6444388.899515686,
  6480957.913874305,
  6491552.758400581,
  6517292.1875,
  6518781.679998049,
  6539635.9375,
  6542820.5052177245,
  6549210.258950732,
  6582260.657392049,
  6595543.591161964,
  6603073.510230481,
  6607812.393037335,
  6613220.352184222,
  6646252.243388383,
  6647116.049445571,
  6647126.276573213,
  6655516.381661021,
  6668478.283208494,
  6683402.641582886,
  6686925.879122401,
  6716466.948577763,
  6721753.125,
  6732743.297955024,
  6735214.759680619,
  6748725.199907782,
  6752338.069254645,
  6754197.583675046,
  6766492.1875,
  6768085.9375,
  6768144.708935806,
  6770055.75081373,
  6775074.9498236235,
  6776185.9375,
  6776710.9375,
  6778199.690927961,
  6778569.971350141,
  6780269.521893895,
  6785831.725897403,
  6785851.564869575,
  6786318.75,
  6786324.64460524,
  6787388.730570981,
  6813995.023475486,
  6818878.125,
  6829907.8125,
  6831415.254066735,
  6835192.790410658,
  6842255.819024344,
  6842885.9375,
  6844991.618727271,
  6846406.25,
  6882015.514014583,
  6886733.696111455,
  6888850.725431256,
  6909095.184679814,
  6909667.607798219,
  6910693.420281266,
  6912022.368144034,
  6912724.015210312,
  6913764.0625,
  6916705.271456131,
  6923899.951813152,
  6928611.065387513,
  6960775.0,
  6961187.347831112,
  6963243.75,
  6963678.675999479,
  6964652.820648306,
  6964930.785821577,
  6965654.6875,
  6969196.875,
  6969388.387159707,
  6975118.75,
  6977514.499256574,
  6982445.3125,
  6985221.679848768,
  6990301.500329422,
  6990960.651806744,
  6993072.2127529215,
  6993353.125,
  6999508.21965047,
  7000230.894740593,
  7010198.4375,
  7016158.6697281515,
  7019546.751769611,
  7021222.548676711,
  7022474.003446295,
  7027130.473610328,
  7028156.866133775,
  7030942.402873203,
  7031249.591404102,
  7033837.5,
  7047135.421266695,
  7050135.308272373,
  7052746.875,
  7068312.037819184,
  7076899.296273707,
  7076900.494830614,
  7087256.006900068,
  7089148.970594312,
  7090970.021301946,
  7091109.0149647035,
  7091255.533601181,
  7092292.068421058,
  7092456.25,
  7094270.493798151,
  7103714.029539407,
  7114850.215585469,
  7117503.0835208045,
  7117551.5625,
  7123250.964185511,
  7129998.873718578,
  7133856.924057238,
  7142151.5625,
  7142738.807376845,
  7143439.692075878,
  7146308.1694153175,
  7147660.809839507,
  7147671.445786514,
  7159628.125,
  7159734.705158518,
  7159788.814473894,
  7162144.475014605,
  7162150.0,
  7169017.318519995,
  7169230.019261797,
  7169935.753913544,
  7172920.902465539,
  7176882.652722928,
  7178188.578068851,
  7179061.876614159,
  7179324.403271313,
  7179490.466087875,
  7182333.517403187,
  7187177.27434857,
  7188528.1159116095,
  7191725.0,
  7191895.154397505,
  7204130.180128741,
  7204169.672283177,
  7209006.853022316,
  7217865.042091088,
  7219551.4300857885,
  7225241.972100578,
  7226648.4375,
  7227734.8177217785,
  7229245.3125,
  7229829.846345175,
  7230193.53695676,
  7231842.340945155,
  7231888.078734148,
  7232364.709200862,
  7232372.301829408,
  7232521.109969386,
  7232741.954768187,
  7233914.0625,
  7234162.5,
  7235432.499786478,
  7239157.1635331595,
  7244325.494834697,
  7244685.20464424,
  7245105.77201253,
  7255570.776843729,
  7257675.386293524,
  7258945.510961343,
  7259085.837440219,
  7265696.147384341,
  7266925.921563234,
  7269174.253436861,
  7272264.536199407,
  7272599.534228866,
  7276011.67404326,
  7284294.378563072,
  7287446.299962416,
  7288511.119753847,
  7290470.801825992,
  7306533.569651871,
  7310043.7140820995,
  7324809.375,
  7327637.5,
  7327841.897956264,
  7338386.692230745,
  7340909.6443142975,
  7347461.788161518,
  7354156.42569781,
  7356116.348642149,
  7356903.125,
  7358012.5,
  7360087.0025309045,
  7360179.6875,
  7360436.249767651,
  7361979.072640978,
  7363511.342652123,
  7365266.978793685,
  7372615.454972976,
  7373492.689963688,
  7373678.2771199625,
  7377810.864212106,
  7378082.156033357,
  7378364.0625,
  7381698.10140834,
  7382517.1875,
  7386766.038602638,
  7387265.698911208,
  7389809.375,
  7389899.238232341,
  7392726.5625,
  7393601.884181978,
  7396481.10109528,
  7397451.134688575,
  7400292.110994169,
  7414464.790526245,
  7421528.96751886,
  7422808.113973377,
  7427717.572494677,
  7427942.969064367,
  7429281.567143573,
  7430281.25,
  7439324.159790136,
  7442349.701613312,
  7444656.225462097,
  7453937.14808503,
  7454154.95430588,
  7454242.1875,
  7455711.802736418,
  7455733.876917875,
  7456356.651171738,
  7457641.5039581135,
  7466146.20728276,
  7484371.34212463,
  7497898.4375,
  7521130.884006161,
  7521133.164201415,
  7523756.166589465,
  7528310.904175868,
  7528841.000604379,
  7529612.788504662,
  7553704.663623589,
  7555316.991390835,
  7556739.0625,
  7556889.0625,
  7558663.637428357,
  7558677.1696301,
  7560010.576779641,
  7564865.617537673,
  7577990.814723356,
  7581129.399162411,
  7583572.602742677,
  7604113.274339929,
  7614928.125,
  7635136.702148178,
  7635627.514356926,
  7650650.993156718,
  7675844.6200247295,
  7677254.996354495,
  7719990.113402434,
  7785012.5,
  7799857.258461481,
  7816915.831979305,
  7839362.462219882,
  7840793.03806642,
  7843795.748948659,
  7844828.306650433,
  7848439.0625,
  7879345.8669913,
  7938216.936949867,
  7941979.6875,
  7948668.1748738,
  7973174.739725015,
  7991358.136018462,
  7992326.418546921,
  7995151.802341818,
  7997514.0625,
  8000912.506359556,
  8007521.716529981,
  8055546.736555785,
  8116830.846602414,
  8172903.262191973,
  8218063.744408554,
  8242548.490278894,
  8366773.40859431,
  8431325.0,
  8451342.619312404,
  8523376.450295145,
  8643443.214097332,
  8677638.948154803,
  8691441.30937179,
  8694154.6875,
  8718578.960704088,
  8757082.131675767,
  8757763.808451075,
  8772224.543148603,
  8835580.384768358,
  8843075.570937566,
  8920484.780158518,
  8940601.811470691,
  9273605.308321143,
  9275081.547567226,
  9277136.699257785,
  9277606.25,
  9289296.243383275],
 [28.047181200786373,
  52.81288382662554,
  59.31369258846639,
  40.69200492058198,
  10.850601384498821,
  7.24664611680004,
  31.970637247556404,
  28.043412397334862,
  65.03719527799717,
  31.362228493773916,
  25.086719456168513,
  43.67344923578488,
  32.69575665247008,
  42.470088496993576,
  52.45252826415969,
  36.729982506238834,
  44.42469951243292,
  56.75478974688252,
  39.76111884056752,
  28.259417068231965,
  60.856113731663626,
  61.134388504451756,
  8.897786059095294,
  72.95307006379514,
  16.43748203353476,
  73.79465463076647,
  49.3555012036899,
  5.307767432304067,
  72.60226357392536,
  37.54013743438695,
  7.5631511946967755,
  5.968155570067807,
  12.114294441841926,
  18.808228627338085,
  30.341201040163753,
  10.13838832581586,
  19.92656567920134,
  9.956742120988913,
  10.59383953246419,
  37.638482925576454,
  13.027112649553304,
  104.5680754853628,
  31.192497215454427,
  11.538213640744566,
  7.431092406943484,
  7.224266055269349,
  7.344660724805491,
  72.80697991499832,
  22.06348400138913,
  65.29178221407088,
  42.02122681229442,
  23.07323238386012,
  31.534361727249127,
  27.221789489270407,
  58.45971548475996,
  105.6643194717448,
  99.00583636601502,
  5.792936714093015,
  10.497514063818903,
  7.554005935898385,
  24.689824705063234,
  56.055550161433786,
  70.64195460670825,
  17.837358530517662,
  14.768478906956064,
  19.606577097963736,
  70.5026602558451,
  78.482036985448,
  49.32013134287927,
  8.288980052812715,
  47.98406781095085,
  15.849805727303025,
  24.27743116339196,
  12.09398211413425,
  51.547308986273045,
  9.847462457567982,
  44.38159916805248,
  50.65485108328809,
  23.702264249775965,
  8.28978665800843,
  66.17181027100118,
  91.4198323045596,
  14.9737921348668,
  51.310684655711526,
  17.372509733490787,
  15.257643690153508,
  8.691703287085197,
  16.402029900692686,
  22.62094137742749,
  18.012128151526753,
  49.17384243192064,
  37.696199885455485,
  24.381838420155493,
  5.85782369328562,
  6.217450071357888,
  14.907013515859788,
  33.607056307389094,
  9.869334034060586,
  7.374353577077692,
  74.43320274882579,
  49.04741858392995,
  46.705623942255315,
  46.59001618389503,
  15.338732477102035,
  69.91282537885836,
  101.19001327121701,
  8.341651787255232,
  33.01538228106095,
  10.923088953416391,
  11.541067886493234,
  21.934871954568848,
  24.01699744146813,
  34.13888208575462,
  15.13194648092543,
  14.482010151046797,
  12.051969831645707,
  59.34104475100405,
  9.814166855407752,
  25.576909109825156,
  11.336516989057625,
  39.6921771348044,
  47.941206998740896,
  53.56419613694311,
  6.545795633637129,
  52.87246571316592,
  25.674337526850373,
  31.565633950362958,
  57.35890132679978,
  42.2604644315674,
  52.238804781727524,
  75.1623870146333,
  8.391188575556326,
  7.957630062714162,
  8.14380527455751,
  14.368788222277708,
  19.086254641169877,
  60.9260738949855,
  9.3651853639999,
  11.522857634870155,
  14.411636771922012,
  6.058861630115405,
  84.15005757846126,
  7.661479135078816,
  21.66013388368943,
  97.79679500523503,
  5.521596170587194,
  55.564335532045135,
  33.73544726902156,
  18.487304981832192,
  9.506356982815628,
  9.864069599067932,
  44.594394375872525,
  49.30671294942751,
  33.48682953668613,
  8.86836478140469,
  15.513177138990631,
  27.56960062971737,
  70.34678586410934,
  86.66466858117408,
  89.74481291452517,
  13.783224221923044,
  101.30558214444994,
  71.00549986320404,
  5.706908830121373,
  48.74464751326761,
  7.4795636354752,
  15.835650655124148,
  39.558171187366064,
  48.36281442949397,
  42.40559492575114,
  15.202149361401078,
  76.60791105328192,
  66.93564275448935,
  53.76901450183988,
  20.83343363379231,
  63.383141669967394,
  26.508043366522575,
  7.374860489263898,
  111.60627292863768,
  26.231086484484873,
  26.63811785543406,
  101.29710657659585,
  29.99626990983278,
  23.432533025685906,
  41.539127978977405,
  17.7584959121628,
  11.053213824559645,
  15.452440843790876,
  87.83925683016827,
  25.842537887786637,
  43.092185394986664,
  47.51679330241817,
  36.305240942867144,
  71.69897097833267,
  29.282197138825836,
  39.72349837741832,
  12.87018189772502,
  41.929516214833896,
  27.803152458568256,
  33.19605701857011,
  24.45584729317172,
  24.984550494237013,
  9.100977281345864,
  9.634311008054024,
  64.26195053602083,
  43.57758816358024,
  13.740743016965135,
  34.88416294707728,
  27.465612757436936,
  14.4213415140087,
  5.652433363372423,
  143.61430511806967,
  136.6307351503348,
  60.64331279649914,
  15.831285287807232,
  85.14909054097939,
  114.29090349488894,
  32.34129623868612,
  20.29407741279457,
  46.07064296103692,
  21.878849170678983,
  5.763682986339976,
  17.23555506613659,
  7.517617005345822,
  6.336985359769647,
  23.863829237162214,
  23.013014157924616,
  35.50187985093534,
  69.75234977959825,
  73.93398631485508,
  40.56709089800573,
  8.269223785616477,
  11.48227375436185,
  6.26197195591292,
  35.64957613783387,
  25.920781978934958,
  7.769016795848976,
  7.31067968336355,
  72.1570257568756,
  9.12461284236896,
  18.04292194651764,
  21.6585464804233,
  119.20102109438125,
  22.229543880937424,
  68.20964891617646,
  9.627096908780667,
  72.29075737943148,
  5.367394632394701,
  60.16654452531593,
  15.502689824974281,
  16.37721931217523,
  7.538160316790053,
  34.69434858224939,
  25.151087848858246,
  5.898300123268472,
  26.76445812238596,
  32.5932685831057,
  25.73735804115053,
  27.198229246606548,
  21.45448653916853,
  40.07663218957094,
  60.43544261210587,
  25.74993213709288,
  17.22315094696406,
  20.199801535100832,
  19.236952466723956,
  67.4887344212514,
  14.459736181851824,
  41.79516367053969,
  61.61914865387153,
  26.08822337467371,
  69.00687877303376,
  56.47024004063678,
  24.928899665829253,
  5.622589873109586,
  67.81262017927335,
  71.34666466748479,
  18.410410968750718,
  23.40205004697094,
  9.142074244665693,
  61.998464340584185,
  81.11267108826205,
  97.46048904411307,
  34.315878195827,
  104.81459609577361,
  109.70204905321027,
  69.0909648679341,
  9.757527046146711,
  16.651252730563275,
  41.163571636702414,
  122.48669921452844,
  19.487890020845242,
  7.974422980307862,
  50.96902686071468,
  8.554038422124883,
  74.8051497672068,
  53.50039223560705,
  20.48430259352304,
  16.609040969886465,
  11.42505405219786,
  12.893670186861918,
  33.880880470375644,
  43.71588814504857,
  46.871822743386595,
  103.67199847158741,
  23.885231191208895,
  10.06304574844553,
  16.078014916554945,
  12.78104370778732,
  18.622448303905944,
  18.060756902927515,
  73.43940217994044,
  66.28429861433992,
  33.7469514556988,
  96.14622532875515,
  7.20603014320438,
  23.671345658718508,
  127.63036152868816,
  24.99141476271382,
  19.162497242442864,
  17.27874101672062,
  104.23407114531972,
  39.51113176772607,
  15.531284575352228,
  56.656045073872065,
  23.300414340189814,
  56.229104751129526,
  6.3209076637175645,
  7.630190659419213,
  101.75724331064139,
  8.760858435208862,
  64.95998502404379,
  42.63566003325964,
  8.209373014591698,
  5.382367230461233,
  11.010192493950838,
  31.730129881206008,
  59.7961105572717,
  61.00890953043228,
  11.787984679318953,
  26.08162639135969,
  6.828517735875828,
  104.77758392222887,
  59.758225657566946,
  23.771371378619747,
  63.31902433000927,
  23.109502804062043,
  16.066306686165007,
  66.76119973011633,
  7.949332276823246,
  20.278865788275453,
  73.72465826588441,
  19.838723548503427,
  44.15800448610915,
  19.538375771980746,
  33.89569259552787,
  62.85564904220214,
  16.47599408906276,
  45.546034296720975,
  96.58589807900131,
  16.59655903058439,
  86.07197862003892,
  8.273991310042685,
  14.243248335825816,
  49.42531714524736,
  63.163418310211085,
  7.456406526752239,
  17.901488081462286,
  9.043946617135283,
  70.17442404541902,
  34.32688785156431,
  72.48375772278058,
  7.326862348844956,
  7.909564405435314,
  36.61522691362226,
  60.86064750425102,
  8.863442832315798,
  13.140060231534514,
  42.57211126826223,
  30.484336342175855,
  25.879142172184977,
  43.859931570547474,
  5.774092782442045,
  6.805544068295682,
  16.551502797801795,
  19.61088697561319,
  16.914540103161247,
  33.79618180448718,
  115.23714275188097,
  9.848824150546228,
  71.02888533464782,
  34.568652919955404,
  13.124896907032698,
  10.74387336226806,
  45.71739073710705,
  25.972092200177656,
  15.941908762134132,
  41.86623601512538,
  7.60553944655945,
  10.888724338976685,
  5.526398921887358,
  18.613346748510104,
  7.15798836339459,
  23.522366898924822,
  6.774089408358032,
  8.271350436697809,
  50.38086160128853,
  9.725319042207698,
  7.563153053957649,
  69.30905633930662,
  14.51345069047326,
  55.61935060181155,
  48.584243737369604,
  10.154476525676786,
  19.86088900962037,
  48.729104396324885,
  5.780933189993433,
  60.90510569030765,
  6.342131654471429,
  16.748980263945683,
  11.540519002835799,
  40.28297302471581,
  42.27600040105143,
  42.19140810125096,
  36.22422207889886,
  29.082205074636683])

Energy spectrum

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

([3497714.0379049713,
  3550652.8231579107,
  3597623.4375,
  3646509.905872267,
  3647718.565175181,
  3647873.893272848,
  3662370.3125,
  3712370.6367102996,
  3753201.500114644,
  3759410.5693005226,
  3778331.047855071,
  3778551.501607855,
  3885142.931882811,
  3886921.875,
  3895424.4876268203,
  3897929.6875,
  3999016.6810882944,
  4022796.8401280846,
  4026175.853823638,
  4104931.2364454693,
  4118038.843785412,
  4143906.25,
  4157779.1259316937,
  4237011.1149124345,
  4243728.092050614,
  4250592.1875,
  4272939.0625,
  4285544.280725929,
  4291148.4375,
  4301803.125,
  4318915.512789876,
  4319134.699317267,
  4341363.051125232,
  4358749.144718171,
  4359476.5625,
  4374908.634141004,
  4384786.82718285,
  4573999.094773447,
  4609680.879308154,
  4671587.5,
  4689415.012096634,
  4701065.625,
  4706709.375,
  4723603.372659559,
  4728463.560068321,
  4728487.877198473,
  4729141.411121211,
  4740712.5,
  4748134.236047245,
  4751720.753708192,
  4759046.875,
  4759460.648647059,
  4760292.1875,
  4767408.836273323,
  4778680.46903719,
  4789343.75,
  4812336.541152843,
  4815291.898985861,
  4843171.430583145,
  4845721.736798774,
  4851876.332625404,
  4875370.3125,
  4882090.625,
  4898103.18465866,
  4899745.266350453,
  4912727.28557971,
  4915565.625,
  4919764.18165474,
  4925560.9375,
  4933957.817040829,
  4962202.643552519,
  5007348.657002478,
  5054886.353251419,
  5154244.291191751,
  5552987.794967575,
  5576033.3656686945,
  5667795.779365314,
  5707228.125,
  5711295.008695379,
  5713535.943563967,
  5748974.027173245,
  5759927.501234855,
  5772332.820944044,
  5773675.0,
  5780288.018863097,
  5781492.843685747,
  5789097.569837585,
  5812277.270626946,
  5855300.946097331,
  5881275.470204075,
  5898355.417147607,
  5936006.25,
  5970291.935350806,
  5981324.012971827,
  6002224.595829856,
  6034448.42842006,
  6060244.896839247,
  6074544.324333707,
  6089863.381991046,
  6103499.3681991175,
  6147611.6818406815,
  6151015.625,
  6174600.64758485,
  6176357.690619808,
  6185452.579336756,
  6192746.875,
  6200629.790214532,
  6216687.5,
  6236551.427525158,
  6254308.061782735,
  6285741.730673833,
  6289832.533062202,
  6323686.032446726,
  6325717.236582025,
  6332524.873985,
  6336617.234567944,
  6336634.375,
  6365217.708124221,
  6381674.0505271135,
  6386181.379605349,
  6403343.095038348,
  6407563.196922803,
  6410983.086725494,
  6413475.364538637,
  6429345.3125,
  6444388.899515686,
  6480957.913874305,
  6491552.758400581,
  6517292.1875,
  6518781.679998049,
  6539635.9375,
  6542820.5052177245,
  6549210.258950732,
  6582260.657392049,
  6595543.591161964,
  6603073.510230481,
  6607812.393037335,
  6613220.352184222,
  6646252.243388383,
  6647116.049445571,
  6647126.276573213,
  6655516.381661021,
  6668478.283208494,
  6683402.641582886,
  6686925.879122401,
  6716466.948577763,
  6721753.125,
  6732743.297955024,
  6735214.759680619,
  6748725.199907782,
  6752338.069254645,
  6754197.583675046,
  6766492.1875,
  6768085.9375,
  6768144.708935806,
  6770055.75081373,
  6775074.9498236235,
  6776185.9375,
  6776710.9375,
  6778199.690927961,
  6778569.971350141,
  6780269.521893895,
  6785831.725897403,
  6785851.564869575,
  6786318.75,
  6786324.64460524,
  6787388.730570981,
  6813995.023475486,
  6818878.125,
  6829907.8125,
  6831415.254066735,
  6835192.790410658,
  6842255.819024344,
  6842885.9375,
  6844991.618727271,
  6846406.25,
  6882015.514014583,
  6886733.696111455,
  6888850.725431256,
  6909095.184679814,
  6909667.607798219,
  6910693.420281266,
  6912022.368144034,
  6912724.015210312,
  6913764.0625,
  6916705.271456131,
  6923899.951813152,
  6928611.065387513,
  6960775.0,
  6961187.347831112,
  6963243.75,
  6963678.675999479,
  6964652.820648306,
  6964930.785821577,
  6965654.6875,
  6969196.875,
  6969388.387159707,
  6975118.75,
  6977514.499256574,
  6982445.3125,
  6985221.679848768,
  6990301.500329422,
  6990960.651806744,
  6993072.2127529215,
  6993353.125,
  6999508.21965047,
  7000230.894740593,
  7010198.4375,
  7016158.6697281515,
  7019546.751769611,
  7021222.548676711,
  7022474.003446295,
  7027130.473610328,
  7028156.866133775,
  7030942.402873203,
  7031249.591404102,
  7033837.5,
  7047135.421266695,
  7050135.308272373,
  7052746.875,
  7068312.037819184,
  7076899.296273707,
  7076900.494830614,
  7087256.006900068,
  7089148.970594312,
  7090970.021301946,
  7091109.0149647035,
  7091255.533601181,
  7092292.068421058,
  7092456.25,
  7094270.493798151,
  7103714.029539407,
  7114850.215585469,
  7117503.0835208045,
  7117551.5625,
  7123250.964185511,
  7129998.873718578,
  7133856.924057238,
  7142151.5625,
  7142738.807376845,
  7143439.692075878,
  7146308.1694153175,
  7147660.809839507,
  7147671.445786514,
  7159628.125,
  7159734.705158518,
  7159788.814473894,
  7162144.475014605,
  7162150.0,
  7169017.318519995,
  7169230.019261797,
  7169935.753913544,
  7172920.902465539,
  7176882.652722928,
  7178188.578068851,
  7179061.876614159,
  7179324.403271313,
  7179490.466087875,
  7182333.517403187,
  7187177.27434857,
  7188528.1159116095,
  7191725.0,
  7191895.154397505,
  7204130.180128741,
  7204169.672283177,
  7209006.853022316,
  7217865.042091088,
  7219551.4300857885,
  7225241.972100578,
  7226648.4375,
  7227734.8177217785,
  7229245.3125,
  7229829.846345175,
  7230193.53695676,
  7231842.340945155,
  7231888.078734148,
  7232364.709200862,
  7232372.301829408,
  7232521.109969386,
  7232741.954768187,
  7233914.0625,
  7234162.5,
  7235432.499786478,
  7239157.1635331595,
  7244325.494834697,
  7244685.20464424,
  7245105.77201253,
  7255570.776843729,
  7257675.386293524,
  7258945.510961343,
  7259085.837440219,
  7265696.147384341,
  7266925.921563234,
  7269174.253436861,
  7272264.536199407,
  7272599.534228866,
  7276011.67404326,
  7284294.378563072,
  7287446.299962416,
  7288511.119753847,
  7290470.801825992,
  7306533.569651871,
  7310043.7140820995,
  7324809.375,
  7327637.5,
  7327841.897956264,
  7338386.692230745,
  7340909.6443142975,
  7347461.788161518,
  7354156.42569781,
  7356116.348642149,
  7356903.125,
  7358012.5,
  7360087.0025309045,
  7360179.6875,
  7360436.249767651,
  7361979.072640978,
  7363511.342652123,
  7365266.978793685,
  7372615.454972976,
  7373492.689963688,
  7373678.2771199625,
  7377810.864212106,
  7378082.156033357,
  7378364.0625,
  7381698.10140834,
  7382517.1875,
  7386766.038602638,
  7387265.698911208,
  7389809.375,
  7389899.238232341,
  7392726.5625,
  7393601.884181978,
  7396481.10109528,
  7397451.134688575,
  7400292.110994169,
  7414464.790526245,
  7421528.96751886,
  7422808.113973377,
  7427717.572494677,
  7427942.969064367,
  7429281.567143573,
  7430281.25,
  7439324.159790136,
  7442349.701613312,
  7444656.225462097,
  7453937.14808503,
  7454154.95430588,
  7454242.1875,
  7455711.802736418,
  7455733.876917875,
  7456356.651171738,
  7457641.5039581135,
  7466146.20728276,
  7484371.34212463,
  7497898.4375,
  7521130.884006161,
  7521133.164201415,
  7523756.166589465,
  7528310.904175868,
  7528841.000604379,
  7529612.788504662,
  7553704.663623589,
  7555316.991390835,
  7556739.0625,
  7556889.0625,
  7558663.637428357,
  7558677.1696301,
  7560010.576779641,
  7564865.617537673,
  7577990.814723356,
  7581129.399162411,
  7583572.602742677,
  7604113.274339929,
  7614928.125,
  7635136.702148178,
  7635627.514356926,
  7650650.993156718,
  7675844.6200247295,
  7677254.996354495,
  7719990.113402434,
  7785012.5,
  7799857.258461481,
  7816915.831979305,
  7839362.462219882,
  7840793.03806642,
  7843795.748948659,
  7844828.306650433,
  7848439.0625,
  7879345.8669913,
  7938216.936949867,
  7941979.6875,
  7948668.1748738,
  7973174.739725015,
  7991358.136018462,
  7992326.418546921,
  7995151.802341818,
  7997514.0625,
  8000912.506359556,
  8007521.716529981,
  8055546.736555785,
  8116830.846602414,
  8172903.262191973,
  8218063.744408554,
  8242548.490278894,
  8366773.40859431,
  8431325.0,
  8451342.619312404,
  8523376.450295145,
  8643443.214097332,
  8677638.948154803,
  8691441.30937179,
  8694154.6875,
  8718578.960704088,
  8757082.131675767,
  8757763.808451075,
  8772224.543148603,
  8835580.384768358,
  8843075.570937566,
  8920484.780158518,
  8940601.811470691,
  9273605.308321143,
  9275081.547567226,
  9277136.699257785,
  9277606.25,
  9289296.243383275],
 [28.047181200786373,
  52.81288382662554,
  59.31369258846639,
  40.69200492058198,
  10.850601384498821,
  7.24664611680004,
  31.970637247556404,
  28.043412397334862,
  65.03719527799717,
  31.362228493773916,
  25.086719456168513,
  43.67344923578488,
  32.69575665247008,
  42.470088496993576,
  52.45252826415969,
  36.729982506238834,
  44.42469951243292,
  56.75478974688252,
  39.76111884056752,
  28.259417068231965,
  60.856113731663626,
  61.134388504451756,
  8.897786059095294,
  72.95307006379514,
  16.43748203353476,
  73.79465463076647,
  49.3555012036899,
  5.307767432304067,
  72.60226357392536,
  37.54013743438695,
  7.5631511946967755,
  5.968155570067807,
  12.114294441841926,
  18.808228627338085,
  30.341201040163753,
  10.13838832581586,
  19.92656567920134,
  9.956742120988913,
  10.59383953246419,
  37.638482925576454,
  13.027112649553304,
  104.5680754853628,
  31.192497215454427,
  11.538213640744566,
  7.431092406943484,
  7.224266055269349,
  7.344660724805491,
  72.80697991499832,
  22.06348400138913,
  65.29178221407088,
  42.02122681229442,
  23.07323238386012,
  31.534361727249127,
  27.221789489270407,
  58.45971548475996,
  105.6643194717448,
  99.00583636601502,
  5.792936714093015,
  10.497514063818903,
  7.554005935898385,
  24.689824705063234,
  56.055550161433786,
  70.64195460670825,
  17.837358530517662,
  14.768478906956064,
  19.606577097963736,
  70.5026602558451,
  78.482036985448,
  49.32013134287927,
  8.288980052812715,
  47.98406781095085,
  15.849805727303025,
  24.27743116339196,
  12.09398211413425,
  51.547308986273045,
  9.847462457567982,
  44.38159916805248,
  50.65485108328809,
  23.702264249775965,
  8.28978665800843,
  66.17181027100118,
  91.4198323045596,
  14.9737921348668,
  51.310684655711526,
  17.372509733490787,
  15.257643690153508,
  8.691703287085197,
  16.402029900692686,
  22.62094137742749,
  18.012128151526753,
  49.17384243192064,
  37.696199885455485,
  24.381838420155493,
  5.85782369328562,
  6.217450071357888,
  14.907013515859788,
  33.607056307389094,
  9.869334034060586,
  7.374353577077692,
  74.43320274882579,
  49.04741858392995,
  46.705623942255315,
  46.59001618389503,
  15.338732477102035,
  69.91282537885836,
  101.19001327121701,
  8.341651787255232,
  33.01538228106095,
  10.923088953416391,
  11.541067886493234,
  21.934871954568848,
  24.01699744146813,
  34.13888208575462,
  15.13194648092543,
  14.482010151046797,
  12.051969831645707,
  59.34104475100405,
  9.814166855407752,
  25.576909109825156,
  11.336516989057625,
  39.6921771348044,
  47.941206998740896,
  53.56419613694311,
  6.545795633637129,
  52.87246571316592,
  25.674337526850373,
  31.565633950362958,
  57.35890132679978,
  42.2604644315674,
  52.238804781727524,
  75.1623870146333,
  8.391188575556326,
  7.957630062714162,
  8.14380527455751,
  14.368788222277708,
  19.086254641169877,
  60.9260738949855,
  9.3651853639999,
  11.522857634870155,
  14.411636771922012,
  6.058861630115405,
  84.15005757846126,
  7.661479135078816,
  21.66013388368943,
  97.79679500523503,
  5.521596170587194,
  55.564335532045135,
  33.73544726902156,
  18.487304981832192,
  9.506356982815628,
  9.864069599067932,
  44.594394375872525,
  49.30671294942751,
  33.48682953668613,
  8.86836478140469,
  15.513177138990631,
  27.56960062971737,
  70.34678586410934,
  86.66466858117408,
  89.74481291452517,
  13.783224221923044,
  101.30558214444994,
  71.00549986320404,
  5.706908830121373,
  48.74464751326761,
  7.4795636354752,
  15.835650655124148,
  39.558171187366064,
  48.36281442949397,
  42.40559492575114,
  15.202149361401078,
  76.60791105328192,
  66.93564275448935,
  53.76901450183988,
  20.83343363379231,
  63.383141669967394,
  26.508043366522575,
  7.374860489263898,
  111.60627292863768,
  26.231086484484873,
  26.63811785543406,
  101.29710657659585,
  29.99626990983278,
  23.432533025685906,
  41.539127978977405,
  17.7584959121628,
  11.053213824559645,
  15.452440843790876,
  87.83925683016827,
  25.842537887786637,
  43.092185394986664,
  47.51679330241817,
  36.305240942867144,
  71.69897097833267,
  29.282197138825836,
  39.72349837741832,
  12.87018189772502,
  41.929516214833896,
  27.803152458568256,
  33.19605701857011,
  24.45584729317172,
  24.984550494237013,
  9.100977281345864,
  9.634311008054024,
  64.26195053602083,
  43.57758816358024,
  13.740743016965135,
  34.88416294707728,
  27.465612757436936,
  14.4213415140087,
  5.652433363372423,
  143.61430511806967,
  136.6307351503348,
  60.64331279649914,
  15.831285287807232,
  85.14909054097939,
  114.29090349488894,
  32.34129623868612,
  20.29407741279457,
  46.07064296103692,
  21.878849170678983,
  5.763682986339976,
  17.23555506613659,
  7.517617005345822,
  6.336985359769647,
  23.863829237162214,
  23.013014157924616,
  35.50187985093534,
  69.75234977959825,
  73.93398631485508,
  40.56709089800573,
  8.269223785616477,
  11.48227375436185,
  6.26197195591292,
  35.64957613783387,
  25.920781978934958,
  7.769016795848976,
  7.31067968336355,
  72.1570257568756,
  9.12461284236896,
  18.04292194651764,
  21.6585464804233,
  119.20102109438125,
  22.229543880937424,
  68.20964891617646,
  9.627096908780667,
  72.29075737943148,
  5.367394632394701,
  60.16654452531593,
  15.502689824974281,
  16.37721931217523,
  7.538160316790053,
  34.69434858224939,
  25.151087848858246,
  5.898300123268472,
  26.76445812238596,
  32.5932685831057,
  25.73735804115053,
  27.198229246606548,
  21.45448653916853,
  40.07663218957094,
  60.43544261210587,
  25.74993213709288,
  17.22315094696406,
  20.199801535100832,
  19.236952466723956,
  67.4887344212514,
  14.459736181851824,
  41.79516367053969,
  61.61914865387153,
  26.08822337467371,
  69.00687877303376,
  56.47024004063678,
  24.928899665829253,
  5.622589873109586,
  67.81262017927335,
  71.34666466748479,
  18.410410968750718,
  23.40205004697094,
  9.142074244665693,
  61.998464340584185,
  81.11267108826205,
  97.46048904411307,
  34.315878195827,
  104.81459609577361,
  109.70204905321027,
  69.0909648679341,
  9.757527046146711,
  16.651252730563275,
  41.163571636702414,
  122.48669921452844,
  19.487890020845242,
  7.974422980307862,
  50.96902686071468,
  8.554038422124883,
  74.8051497672068,
  53.50039223560705,
  20.48430259352304,
  16.609040969886465,
  11.42505405219786,
  12.893670186861918,
  33.880880470375644,
  43.71588814504857,
  46.871822743386595,
  103.67199847158741,
  23.885231191208895,
  10.06304574844553,
  16.078014916554945,
  12.78104370778732,
  18.622448303905944,
  18.060756902927515,
  73.43940217994044,
  66.28429861433992,
  33.7469514556988,
  96.14622532875515,
  7.20603014320438,
  23.671345658718508,
  127.63036152868816,
  24.99141476271382,
  19.162497242442864,
  17.27874101672062,
  104.23407114531972,
  39.51113176772607,
  15.531284575352228,
  56.656045073872065,
  23.300414340189814,
  56.229104751129526,
  6.3209076637175645,
  7.630190659419213,
  101.75724331064139,
  8.760858435208862,
  64.95998502404379,
  42.63566003325964,
  8.209373014591698,
  5.382367230461233,
  11.010192493950838,
  31.730129881206008,
  59.7961105572717,
  61.00890953043228,
  11.787984679318953,
  26.08162639135969,
  6.828517735875828,
  104.77758392222887,
  59.758225657566946,
  23.771371378619747,
  63.31902433000927,
  23.109502804062043,
  16.066306686165007,
  66.76119973011633,
  7.949332276823246,
  20.278865788275453,
  73.72465826588441,
  19.838723548503427,
  44.15800448610915,
  19.538375771980746,
  33.89569259552787,
  62.85564904220214,
  16.47599408906276,
  45.546034296720975,
  96.58589807900131,
  16.59655903058439,
  86.07197862003892,
  8.273991310042685,
  14.243248335825816,
  49.42531714524736,
  63.163418310211085,
  7.456406526752239,
  17.901488081462286,
  9.043946617135283,
  70.17442404541902,
  34.32688785156431,
  72.48375772278058,
  7.326862348844956,
  7.909564405435314,
  36.61522691362226,
  60.86064750425102,
  8.863442832315798,
  13.140060231534514,
  42.57211126826223,
  30.484336342175855,
  25.879142172184977,
  43.859931570547474,
  5.774092782442045,
  6.805544068295682,
  16.551502797801795,
  19.61088697561319,
  16.914540103161247,
  33.79618180448718,
  115.23714275188097,
  9.848824150546228,
  71.02888533464782,
  34.568652919955404,
  13.124896907032698,
  10.74387336226806,
  45.71739073710705,
  25.972092200177656,
  15.941908762134132,
  41.86623601512538,
  7.60553944655945,
  10.888724338976685,
  5.526398921887358,
  18.613346748510104,
  7.15798836339459,
  23.522366898924822,
  6.774089408358032,
  8.271350436697809,
  50.38086160128853,
  9.725319042207698,
  7.563153053957649,
  69.30905633930662,
  14.51345069047326,
  55.61935060181155,
  48.584243737369604,
  10.154476525676786,
  19.86088900962037,
  48.729104396324885,
  5.780933189993433,
  60.90510569030765,
  6.342131654471429,
  16.748980263945683,
  11.540519002835799,
  40.28297302471581,
  42.27600040105143,
  42.19140810125096,
  36.22422207889886,
  29.082205074636683])

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

([3497714.0379049713,
  3550652.8231579107,
  3597623.4375,
  3646509.905872267,
  3647718.565175181,
  3647873.893272848,
  3662370.3125,
  3712370.6367102996,
  3753201.500114644,
  3759410.5693005226,
  3778331.047855071,
  3778551.501607855,
  3885142.931882811,
  3886921.875,
  3895424.4876268203,
  3897929.6875,
  3999016.6810882944,
  4022796.8401280846,
  4026175.853823638,
  4104931.2364454693,
  4118038.843785412,
  4143906.25,
  4157779.1259316937,
  4237011.1149124345,
  4243728.092050614,
  4250592.1875,
  4272939.0625,
  4285544.280725929,
  4291148.4375,
  4301803.125,
  4318915.512789876,
  4319134.699317267,
  4341363.051125232,
  4358749.144718171,
  4359476.5625,
  4374908.634141004,
  4384786.82718285,
  4573999.094773447,
  4609680.879308154,
  4671587.5,
  4689415.012096634,
  4701065.625,
  4706709.375,
  4723603.372659559,
  4728463.560068321,
  4728487.877198473,
  4729141.411121211,
  4740712.5,
  4748134.236047245,
  4751720.753708192,
  4759046.875,
  4759460.648647059,
  4760292.1875,
  4767408.836273323,
  4778680.46903719,
  4789343.75,
  4812336.541152843,
  4815291.898985861,
  4843171.430583145,
  4845721.736798774,
  4851876.332625404,
  4875370.3125,
  4882090.625,
  4898103.18465866,
  4899745.266350453,
  4912727.28557971,
  4915565.625,
  4919764.18165474,
  4925560.9375,
  4933957.817040829,
  4962202.643552519,
  5007348.657002478,
  5054886.353251419,
  5154244.291191751,
  5552987.794967575,
  5576033.3656686945,
  5667795.779365314,
  5707228.125,
  5711295.008695379,
  5713535.943563967,
  5748974.027173245,
  5759927.501234855,
  5772332.820944044,
  5773675.0,
  5780288.018863097,
  5781492.843685747,
  5789097.569837585,
  5812277.270626946,
  5855300.946097331,
  5881275.470204075,
  5898355.417147607,
  5936006.25,
  5970291.935350806,
  5981324.012971827,
  6002224.595829856,
  6034448.42842006,
  6060244.896839247,
  6074544.324333707,
  6089863.381991046,
  6103499.3681991175,
  6147611.6818406815,
  6151015.625,
  6174600.64758485,
  6176357.690619808,
  6185452.579336756,
  6192746.875,
  6200629.790214532,
  6216687.5,
  6236551.427525158,
  6254308.061782735,
  6285741.730673833,
  6289832.533062202,
  6323686.032446726,
  6325717.236582025,
  6332524.873985,
  6336617.234567944,
  6336634.375,
  6365217.708124221,
  6381674.0505271135,
  6386181.379605349,
  6403343.095038348,
  6407563.196922803,
  6410983.086725494,
  6413475.364538637,
  6429345.3125,
  6444388.899515686,
  6480957.913874305,
  6491552.758400581,
  6517292.1875,
  6518781.679998049,
  6539635.9375,
  6542820.5052177245,
  6549210.258950732,
  6582260.657392049,
  6595543.591161964,
  6603073.510230481,
  6607812.393037335,
  6613220.352184222,
  6646252.243388383,
  6647116.049445571,
  6647126.276573213,
  6655516.381661021,
  6668478.283208494,
  6683402.641582886,
  6686925.879122401,
  6716466.948577763,
  6721753.125,
  6732743.297955024,
  6735214.759680619,
  6748725.199907782,
  6752338.069254645,
  6754197.583675046,
  6766492.1875,
  6768085.9375,
  6768144.708935806,
  6770055.75081373,
  6775074.9498236235,
  6776185.9375,
  6776710.9375,
  6778199.690927961,
  6778569.971350141,
  6780269.521893895,
  6785831.725897403,
  6785851.564869575,
  6786318.75,
  6786324.64460524,
  6787388.730570981,
  6813995.023475486,
  6818878.125,
  6829907.8125,
  6831415.254066735,
  6835192.790410658,
  6842255.819024344,
  6842885.9375,
  6844991.618727271,
  6846406.25,
  6882015.514014583,
  6886733.696111455,
  6888850.725431256,
  6909095.184679814,
  6909667.607798219,
  6910693.420281266,
  6912022.368144034,
  6912724.015210312,
  6913764.0625,
  6916705.271456131,
  6923899.951813152,
  6928611.065387513,
  6960775.0,
  6961187.347831112,
  6963243.75,
  6963678.675999479,
  6964652.820648306,
  6964930.785821577,
  6965654.6875,
  6969196.875,
  6969388.387159707,
  6975118.75,
  6977514.499256574,
  6982445.3125,
  6985221.679848768,
  6990301.500329422,
  6990960.651806744,
  6993072.2127529215,
  6993353.125,
  6999508.21965047,
  7000230.894740593,
  7010198.4375,
  7016158.6697281515,
  7019546.751769611,
  7021222.548676711,
  7022474.003446295,
  7027130.473610328,
  7028156.866133775,
  7030942.402873203,
  7031249.591404102,
  7033837.5,
  7047135.421266695,
  7050135.308272373,
  7052746.875,
  7068312.037819184,
  7076899.296273707,
  7076900.494830614,
  7087256.006900068,
  7089148.970594312,
  7090970.021301946,
  7091109.0149647035,
  7091255.533601181,
  7092292.068421058,
  7092456.25,
  7094270.493798151,
  7103714.029539407,
  7114850.215585469,
  7117503.0835208045,
  7117551.5625,
  7123250.964185511,
  7129998.873718578,
  7133856.924057238,
  7142151.5625,
  7142738.807376845,
  7143439.692075878,
  7146308.1694153175,
  7147660.809839507,
  7147671.445786514,
  7159628.125,
  7159734.705158518,
  7159788.814473894,
  7162144.475014605,
  7162150.0,
  7169017.318519995,
  7169230.019261797,
  7169935.753913544,
  7172920.902465539,
  7176882.652722928,
  7178188.578068851,
  7179061.876614159,
  7179324.403271313,
  7179490.466087875,
  7182333.517403187,
  7187177.27434857,
  7188528.1159116095,
  7191725.0,
  7191895.154397505,
  7204130.180128741,
  7204169.672283177,
  7209006.853022316,
  7217865.042091088,
  7219551.4300857885,
  7225241.972100578,
  7226648.4375,
  7227734.8177217785,
  7229245.3125,
  7229829.846345175,
  7230193.53695676,
  7231842.340945155,
  7231888.078734148,
  7232364.709200862,
  7232372.301829408,
  7232521.109969386,
  7232741.954768187,
  7233914.0625,
  7234162.5,
  7235432.499786478,
  7239157.1635331595,
  7244325.494834697,
  7244685.20464424,
  7245105.77201253,
  7255570.776843729,
  7257675.386293524,
  7258945.510961343,
  7259085.837440219,
  7265696.147384341,
  7266925.921563234,
  7269174.253436861,
  7272264.536199407,
  7272599.534228866,
  7276011.67404326,
  7284294.378563072,
  7287446.299962416,
  7288511.119753847,
  7290470.801825992,
  7306533.569651871,
  7310043.7140820995,
  7324809.375,
  7327637.5,
  7327841.897956264,
  7338386.692230745,
  7340909.6443142975,
  7347461.788161518,
  7354156.42569781,
  7356116.348642149,
  7356903.125,
  7358012.5,
  7360087.0025309045,
  7360179.6875,
  7360436.249767651,
  7361979.072640978,
  7363511.342652123,
  7365266.978793685,
  7372615.454972976,
  7373492.689963688,
  7373678.2771199625,
  7377810.864212106,
  7378082.156033357,
  7378364.0625,
  7381698.10140834,
  7382517.1875,
  7386766.038602638,
  7387265.698911208,
  7389809.375,
  7389899.238232341,
  7392726.5625,
  7393601.884181978,
  7396481.10109528,
  7397451.134688575,
  7400292.110994169,
  7414464.790526245,
  7421528.96751886,
  7422808.113973377,
  7427717.572494677,
  7427942.969064367,
  7429281.567143573,
  7430281.25,
  7439324.159790136,
  7442349.701613312,
  7444656.225462097,
  7453937.14808503,
  7454154.95430588,
  7454242.1875,
  7455711.802736418,
  7455733.876917875,
  7456356.651171738,
  7457641.5039581135,
  7466146.20728276,
  7484371.34212463,
  7497898.4375,
  7521130.884006161,
  7521133.164201415,
  7523756.166589465,
  7528310.904175868,
  7528841.000604379,
  7529612.788504662,
  7553704.663623589,
  7555316.991390835,
  7556739.0625,
  7556889.0625,
  7558663.637428357,
  7558677.1696301,
  7560010.576779641,
  7564865.617537673,
  7577990.814723356,
  7581129.399162411,
  7583572.602742677,
  7604113.274339929,
  7614928.125,
  7635136.702148178,
  7635627.514356926,
  7650650.993156718,
  7675844.6200247295,
  7677254.996354495,
  7719990.113402434,
  7785012.5,
  7799857.258461481,
  7816915.831979305,
  7839362.462219882,
  7840793.03806642,
  7843795.748948659,
  7844828.306650433,
  7848439.0625,
  7879345.8669913,
  7938216.936949867,
  7941979.6875,
  7948668.1748738,
  7973174.739725015,
  7991358.136018462,
  7992326.418546921,
  7995151.802341818,
  7997514.0625,
  8000912.506359556,
  8007521.716529981,
  8055546.736555785,
  8116830.846602414,
  8172903.262191973,
  8218063.744408554,
  8242548.490278894,
  8366773.40859431,
  8431325.0,
  8451342.619312404,
  8523376.450295145,
  8643443.214097332,
  8677638.948154803,
  8691441.30937179,
  8694154.6875,
  8718578.960704088,
  8757082.131675767,
  8757763.808451075,
  8772224.543148603,
  8835580.384768358,
  8843075.570937566,
  8920484.780158518,
  8940601.811470691,
  9273605.308321143,
  9275081.547567226,
  9277136.699257785,
  9277606.25,
  9289296.243383275],
 [28.047181200786373,
  52.81288382662554,
  59.31369258846639,
  40.69200492058198,
  10.850601384498821,
  7.24664611680004,
  31.970637247556404,
  28.043412397334862,
  65.03719527799717,
  31.362228493773916,
  25.086719456168513,
  43.67344923578488,
  32.69575665247008,
  42.470088496993576,
  52.45252826415969,
  36.729982506238834,
  44.42469951243292,
  56.75478974688252,
  39.76111884056752,
  28.259417068231965,
  60.856113731663626,
  61.134388504451756,
  8.897786059095294,
  72.95307006379514,
  16.43748203353476,
  73.79465463076647,
  49.3555012036899,
  5.307767432304067,
  72.60226357392536,
  37.54013743438695,
  7.5631511946967755,
  5.968155570067807,
  12.114294441841926,
  18.808228627338085,
  30.341201040163753,
  10.13838832581586,
  19.92656567920134,
  9.956742120988913,
  10.59383953246419,
  37.638482925576454,
  13.027112649553304,
  104.5680754853628,
  31.192497215454427,
  11.538213640744566,
  7.431092406943484,
  7.224266055269349,
  7.344660724805491,
  72.80697991499832,
  22.06348400138913,
  65.29178221407088,
  42.02122681229442,
  23.07323238386012,
  31.534361727249127,
  27.221789489270407,
  58.45971548475996,
  105.6643194717448,
  99.00583636601502,
  5.792936714093015,
  10.497514063818903,
  7.554005935898385,
  24.689824705063234,
  56.055550161433786,
  70.64195460670825,
  17.837358530517662,
  14.768478906956064,
  19.606577097963736,
  70.5026602558451,
  78.482036985448,
  49.32013134287927,
  8.288980052812715,
  47.98406781095085,
  15.849805727303025,
  24.27743116339196,
  12.09398211413425,
  51.547308986273045,
  9.847462457567982,
  44.38159916805248,
  50.65485108328809,
  23.702264249775965,
  8.28978665800843,
  66.17181027100118,
  91.4198323045596,
  14.9737921348668,
  51.310684655711526,
  17.372509733490787,
  15.257643690153508,
  8.691703287085197,
  16.402029900692686,
  22.62094137742749,
  18.012128151526753,
  49.17384243192064,
  37.696199885455485,
  24.381838420155493,
  5.85782369328562,
  6.217450071357888,
  14.907013515859788,
  33.607056307389094,
  9.869334034060586,
  7.374353577077692,
  74.43320274882579,
  49.04741858392995,
  46.705623942255315,
  46.59001618389503,
  15.338732477102035,
  69.91282537885836,
  101.19001327121701,
  8.341651787255232,
  33.01538228106095,
  10.923088953416391,
  11.541067886493234,
  21.934871954568848,
  24.01699744146813,
  34.13888208575462,
  15.13194648092543,
  14.482010151046797,
  12.051969831645707,
  59.34104475100405,
  9.814166855407752,
  25.576909109825156,
  11.336516989057625,
  39.6921771348044,
  47.941206998740896,
  53.56419613694311,
  6.545795633637129,
  52.87246571316592,
  25.674337526850373,
  31.565633950362958,
  57.35890132679978,
  42.2604644315674,
  52.238804781727524,
  75.1623870146333,
  8.391188575556326,
  7.957630062714162,
  8.14380527455751,
  14.368788222277708,
  19.086254641169877,
  60.9260738949855,
  9.3651853639999,
  11.522857634870155,
  14.411636771922012,
  6.058861630115405,
  84.15005757846126,
  7.661479135078816,
  21.66013388368943,
  97.79679500523503,
  5.521596170587194,
  55.564335532045135,
  33.73544726902156,
  18.487304981832192,
  9.506356982815628,
  9.864069599067932,
  44.594394375872525,
  49.30671294942751,
  33.48682953668613,
  8.86836478140469,
  15.513177138990631,
  27.56960062971737,
  70.34678586410934,
  86.66466858117408,
  89.74481291452517,
  13.783224221923044,
  101.30558214444994,
  71.00549986320404,
  5.706908830121373,
  48.74464751326761,
  7.4795636354752,
  15.835650655124148,
  39.558171187366064,
  48.36281442949397,
  42.40559492575114,
  15.202149361401078,
  76.60791105328192,
  66.93564275448935,
  53.76901450183988,
  20.83343363379231,
  63.383141669967394,
  26.508043366522575,
  7.374860489263898,
  111.60627292863768,
  26.231086484484873,
  26.63811785543406,
  101.29710657659585,
  29.99626990983278,
  23.432533025685906,
  41.539127978977405,
  17.7584959121628,
  11.053213824559645,
  15.452440843790876,
  87.83925683016827,
  25.842537887786637,
  43.092185394986664,
  47.51679330241817,
  36.305240942867144,
  71.69897097833267,
  29.282197138825836,
  39.72349837741832,
  12.87018189772502,
  41.929516214833896,
  27.803152458568256,
  33.19605701857011,
  24.45584729317172,
  24.984550494237013,
  9.100977281345864,
  9.634311008054024,
  64.26195053602083,
  43.57758816358024,
  13.740743016965135,
  34.88416294707728,
  27.465612757436936,
  14.4213415140087,
  5.652433363372423,
  143.61430511806967,
  136.6307351503348,
  60.64331279649914,
  15.831285287807232,
  85.14909054097939,
  114.29090349488894,
  32.34129623868612,
  20.29407741279457,
  46.07064296103692,
  21.878849170678983,
  5.763682986339976,
  17.23555506613659,
  7.517617005345822,
  6.336985359769647,
  23.863829237162214,
  23.013014157924616,
  35.50187985093534,
  69.75234977959825,
  73.93398631485508,
  40.56709089800573,
  8.269223785616477,
  11.48227375436185,
  6.26197195591292,
  35.64957613783387,
  25.920781978934958,
  7.769016795848976,
  7.31067968336355,
  72.1570257568756,
  9.12461284236896,
  18.04292194651764,
  21.6585464804233,
  119.20102109438125,
  22.229543880937424,
  68.20964891617646,
  9.627096908780667,
  72.29075737943148,
  5.367394632394701,
  60.16654452531593,
  15.502689824974281,
  16.37721931217523,
  7.538160316790053,
  34.69434858224939,
  25.151087848858246,
  5.898300123268472,
  26.76445812238596,
  32.5932685831057,
  25.73735804115053,
  27.198229246606548,
  21.45448653916853,
  40.07663218957094,
  60.43544261210587,
  25.74993213709288,
  17.22315094696406,
  20.199801535100832,
  19.236952466723956,
  67.4887344212514,
  14.459736181851824,
  41.79516367053969,
  61.61914865387153,
  26.08822337467371,
  69.00687877303376,
  56.47024004063678,
  24.928899665829253,
  5.622589873109586,
  67.81262017927335,
  71.34666466748479,
  18.410410968750718,
  23.40205004697094,
  9.142074244665693,
  61.998464340584185,
  81.11267108826205,
  97.46048904411307,
  34.315878195827,
  104.81459609577361,
  109.70204905321027,
  69.0909648679341,
  9.757527046146711,
  16.651252730563275,
  41.163571636702414,
  122.48669921452844,
  19.487890020845242,
  7.974422980307862,
  50.96902686071468,
  8.554038422124883,
  74.8051497672068,
  53.50039223560705,
  20.48430259352304,
  16.609040969886465,
  11.42505405219786,
  12.893670186861918,
  33.880880470375644,
  43.71588814504857,
  46.871822743386595,
  103.67199847158741,
  23.885231191208895,
  10.06304574844553,
  16.078014916554945,
  12.78104370778732,
  18.622448303905944,
  18.060756902927515,
  73.43940217994044,
  66.28429861433992,
  33.7469514556988,
  96.14622532875515,
  7.20603014320438,
  23.671345658718508,
  127.63036152868816,
  24.99141476271382,
  19.162497242442864,
  17.27874101672062,
  104.23407114531972,
  39.51113176772607,
  15.531284575352228,
  56.656045073872065,
  23.300414340189814,
  56.229104751129526,
  6.3209076637175645,
  7.630190659419213,
  101.75724331064139,
  8.760858435208862,
  64.95998502404379,
  42.63566003325964,
  8.209373014591698,
  5.382367230461233,
  11.010192493950838,
  31.730129881206008,
  59.7961105572717,
  61.00890953043228,
  11.787984679318953,
  26.08162639135969,
  6.828517735875828,
  104.77758392222887,
  59.758225657566946,
  23.771371378619747,
  63.31902433000927,
  23.109502804062043,
  16.066306686165007,
  66.76119973011633,
  7.949332276823246,
  20.278865788275453,
  73.72465826588441,
  19.838723548503427,
  44.15800448610915,
  19.538375771980746,
  33.89569259552787,
  62.85564904220214,
  16.47599408906276,
  45.546034296720975,
  96.58589807900131,
  16.59655903058439,
  86.07197862003892,
  8.273991310042685,
  14.243248335825816,
  49.42531714524736,
  63.163418310211085,
  7.456406526752239,
  17.901488081462286,
  9.043946617135283,
  70.17442404541902,
  34.32688785156431,
  72.48375772278058,
  7.326862348844956,
  7.909564405435314,
  36.61522691362226,
  60.86064750425102,
  8.863442832315798,
  13.140060231534514,
  42.57211126826223,
  30.484336342175855,
  25.879142172184977,
  43.859931570547474,
  5.774092782442045,
  6.805544068295682,
  16.551502797801795,
  19.61088697561319,
  16.914540103161247,
  33.79618180448718,
  115.23714275188097,
  9.848824150546228,
  71.02888533464782,
  34.568652919955404,
  13.124896907032698,
  10.74387336226806,
  45.71739073710705,
  25.972092200177656,
  15.941908762134132,
  41.86623601512538,
  7.60553944655945,
  10.888724338976685,
  5.526398921887358,
  18.613346748510104,
  7.15798836339459,
  23.522366898924822,
  6.774089408358032,
  8.271350436697809,
  50.38086160128853,
  9.725319042207698,
  7.563153053957649,
  69.30905633930662,
  14.51345069047326,
  55.61935060181155,
  48.584243737369604,
  10.154476525676786,
  19.86088900962037,
  48.729104396324885,
  5.780933189993433,
  60.90510569030765,
  6.342131654471429,
  16.748980263945683,
  11.540519002835799,
  40.28297302471581,
  42.27600040105143,
  42.19140810125096,
  36.22422207889886,
  29.082205074636683])

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)