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

([6370415.625,
  7493778.0300543625,
  7555534.375,
  7962801.129695198,
  7964143.75,
  7967822.601024829,
  7968253.340446217,
  8037436.8105644435,
  8042098.319185842,
  8042109.910255405,
  8075457.405665217,
  8079616.903053936,
  8079639.0625,
  8079639.0625,
  8082554.6875,
  8083149.985643593,
  8083509.375,
  8136437.5,
  8153898.4375,
  8154428.125,
  8174020.418735455,
  8175112.4367218865,
  8175304.6875,
  8176729.6875,
  8180535.107012575,
  8185240.625,
  8186591.116974672,
  8186642.596790322,
  8192020.3125,
  8192148.420630831,
  8192218.933158315,
  8192246.128306476,
  8192264.827273551,
  8192624.537612015,
  8193069.77111207,
  8193671.875,
  8194315.650395631,
  8194335.179351896,
  8194665.625,
  8198375.904300488,
  8199085.9329081075,
  8199108.879339639,
  8199257.788063141,
  8199354.130501113,
  8199365.625,
  8199607.8125,
  8200785.855714465,
  8200825.710530285,
  8200942.4443801725,
  8202126.5625,
  8203862.9426920945,
  8203944.91640513,
  8204364.0625,
  8209717.541984143,
  8212478.143119827,
  8213118.827638851,
  8213290.625,
  8216850.661793346,
  8219637.5,
  8220700.0,
  8222567.1875,
  8225080.360513974,
  8225993.288313369,
  8227054.6875,
  8270126.197484022,
  8271447.153493197,
  8273425.394544437,
  8276015.348340758,
  8276460.626891639,
  8281905.42122556,
  8282150.221704014,
  8282841.598249388,
  8283227.390570716,
  8283833.745651554,
  8283956.573109734,
  8284219.662687637,
  8284299.548969384,
  8284312.5,
  8284615.128322482,
  8284924.827108234,
  8285814.0625,
  8286107.154223692,
  8286499.360096596,
  8287085.9375,
  8288274.273612795,
  8290093.75,
  8405582.8125,
  8408315.40149184,
  8410287.5,
  8412654.6875,
  8414519.486750592,
  8417281.25,
  8417697.575686883,
  8420295.777292263,
  8423079.033895,
  8423748.898787437,
  8426082.188130599,
  8426984.27525257,
  8427049.827288557,
  8429182.921553101,
  8429360.9375,
  8436708.989427311,
  8438760.03815258,
  8441468.41491777,
  8443260.37791533,
  8498326.495300852,
  8503645.84830694,
  8504191.482654564,
  8504353.125,
  8505589.0625,
  8505802.686041163,
  8506169.95587018,
  8506878.572144937,
  8507337.254467461,
  8510127.244231887,
  8510898.080946347,
  8514135.459702633,
  8514590.104032371,
  8514784.375,
  8516200.0,
  8516287.9243066,
  8558854.615864597,
  8558854.728845127,
  8560740.886496548,
  8599375.0,
  8620446.521481171,
  8621854.718896056,
  8621908.004352933,
  8622296.875,
  8622376.5625,
  8623254.6875,
  8623889.16336461,
  8624125.615159553,
  8624963.867425408,
  8625138.678546667,
  8625142.1875,
  8625294.355791416,
  8625304.304217637,
  8625317.418941231,
  8625453.70795514,
  8625489.014595686,
  8625871.875,
  8626309.420773359,
  8626314.34711579,
  8626330.180160021,
  8626500.214021055,
  8626987.82077895,
  8627556.25,
  8627602.648604143,
  8627869.222985258,
  8627925.610292373,
  8628051.477168456,
  8628068.977281945,
  8628273.4375,
  8628335.361533169,
  8628447.67354922,
  8628890.625,
  8628948.167414863,
  8629007.646143977,
  8629078.125,
  8629108.293836331,
  8629115.481515514,
  8629471.875,
  8629586.199304296,
  8629597.238134941,
  8629612.5,
  8629784.309322458,
  8629855.481487747,
  8629856.738010824,
  8629859.624949783,
  8629945.380105993,
  8630056.25,
  8630219.832351342,
  8630350.0,
  8630360.9375,
  8630387.5,
  8630471.875,
  8630526.283756888,
  8630539.183062648,
  8630699.555421626,
  8630931.89094182,
  8630947.674384046,
  8631130.274577262,
  8631195.3125,
  8631726.5625,
  8631736.264508057,
  8631911.756756522,
  8631978.125,
  8632055.865419207,
  8632080.724470418,
  8632136.883974709,
  8632143.75,
  8632181.658875622,
  8632219.08711897,
  8632282.8125,
  8632531.04064556,
  8632536.625366537,
  8632645.24249484,
  8632884.037684638,
  8633043.196403252,
  8633146.695263699,
  8633228.657263376,
  8633828.125,
  8634537.346366692,
  8635047.095092937,
  8635993.75,
  8636199.481363596,
  8636318.043321159,
  8636370.3125,
  8637078.125,
  8637095.709937243,
  8637131.25,
  8637354.6875,
  8637481.070846988,
  8637659.520163998,
  8637735.945668709,
  8637758.100223372,
  8638475.295057032,
  8649636.234736966,
  8709416.469863476,
  8718930.804637767,
  8720954.6875,
  8720991.766065644,
  8721909.246923601,
  8722744.140997415,
  8723401.5625,
  8727917.607537014,
  8728189.743827332,
  8728550.810509421,
  8728790.220477886,
  8728839.807495305,
  8729002.999835046,
  8729168.75,
  8729222.352610271,
  8729265.625,
  8729268.681606099,
  8729285.561634906,
  8729371.270680608,
  8729645.3125,
  8729757.8125,
  8729776.5625,
  8729881.172019469,
  8730502.842709733,
  8730694.568198716,
  8730964.589418175,
  8731267.1875,
  8731302.535127826,
  8731329.340709189,
  8731512.384490192,
  8731532.8125,
  8731535.980263544,
  8731619.996473566,
  8731639.0625,
  8731662.72758086,
  8731764.082063837,
  8731782.8125,
  8731807.677842692,
  8731807.794366304,
  8731813.80389695,
  8731840.504328666,
  8731914.32290047,
  8731964.0625,
  8731992.769423552,
  8732096.875,
  8732143.75,
  8732151.037202025,
  8732214.063766768,
  8732340.598849488,
  8732457.863855394,
  8732539.747332722,
  8732703.125,
  8732704.665772485,
  8732723.4375,
  8732745.20205004,
  8732810.9375,
  8732828.856160294,
  8732873.366202427,
  8732968.99372651,
  8732993.391307807,
  8732995.592692751,
  8733025.646593425,
  8733176.051991904,
  8733509.721739609,
  8733581.25,
  8733712.297023658,
  8733740.582411997,
  8733807.315393174,
  8734290.625,
  8734438.938187826,
  8734583.68221809,
  8734821.029561054,
  8734943.733371412,
  8735001.5625,
  8735223.4375,
  8735657.8125,
  8736065.477876388,
  8736447.690357901,
  8737406.552472765,
  8737790.625,
  8737857.08753616,
  8737931.617944706,
  8738483.967364809,
  8739176.65991593,
  8740860.925612165,
  8742864.0625,
  8774249.843960293,
  8776625.369678203,
  8777550.0,
  8780571.875,
  8782018.75,
  8783288.906315159,
  8939999.679760123,
  8940659.135945218,
  8943594.460837632,
  8943672.788193561,
  8944380.940364253,
  8945540.24073969,
  8946936.601113407,
  8947242.853862116,
  8947565.453072768,
  8947692.965163244,
  8948229.6875,
  8948246.36259576,
  8948370.3125,
  8948792.1875,
  8948849.573628077,
  8949126.031295095,
  8949397.581026677,
  8949678.125,
  8949825.450688124,
  8949933.931057876,
  8950042.961132133,
  8950165.870370278,
  8950512.5,
  8951655.907223774,
  8951675.595838109,
  8952054.458046323,
  8955907.598836772,
  8955968.75,
  8965944.842796069,
  8965960.221309869,
  8967595.721422737,
  9014307.874400143,
  9015828.214325279,
  9018188.602514023,
  9019171.344940957,
  9023698.782810353,
  9142423.4375,
  9149682.8125,
  9149684.375,
  9160189.0625,
  9160331.955291709,
  9160423.710363453,
  9161468.432930669,
  9189379.00459956,
  9190366.260583293,
  9192260.354203666,
  9192628.125,
  9192715.625,
  9192741.383149007,
  9196253.125,
  9197878.725043751,
  9198312.5,
  9202979.6875,
  9202982.688012155,
  9223553.073441105,
  9241754.6875,
  9252871.91279967,
  9257876.55002248,
  9259709.375,
  9261424.756808786,
  9262584.375,
  9262589.686682368,
  9262654.297321755,
  9263220.755100373,
  9263274.411818977,
  9263409.206069665,
  9264044.421800464,
  9264051.5625,
  9266150.0,
  9266654.646965856,
  9266663.564005904,
  9267337.222910373,
  9267752.861241207,
  9267853.366124198,
  9268108.020947196,
  9268181.25,
  9268620.3125,
  9268640.294603912,
  9268651.42046744,
  9268667.1875,
  9268970.405585008,
  9269066.05474343,
  9269464.0625,
  9269828.862440998,
  9270009.375,
  9270177.310870457,
  9270337.288028244,
  9270753.125,
  9272040.625,
  9272367.159734856,
  9272393.317499336,
  9272560.9375,
  9272563.762278028,
  9273904.403491234,
  9274265.625,
  9274285.9375,
  9276315.625,
  9276516.907272683,
  9276638.694584612,
  9277175.790397238,
  9278025.055010123,
  9278710.607178574,
  9278831.034385173,
  9283077.211784298,
  9284546.116646703,
  9285478.943681853,
  9287675.0,
  9295871.512814384,
  9377578.8490112,
  9394276.331778187,
  9402774.743281849,
  9459024.175972257,
  9471436.258415814,
  9475014.0625,
  9478514.0625,
  9479576.115839161,
  9481256.819753585,
  9482750.0,
  9484232.146311082,
  9484928.125,
  9485565.185714068,
  9498817.157990567,
  9509428.125,
  9510382.077089146,
  9512260.637481933,
  9514498.4375,
  9516859.960155206,
  9518134.375,
  9518170.909750745,
  9521739.455142671,
  9526175.0,
  9526532.8125,
  9548804.08259865,
  9549356.25,
  9549737.5,
  9550221.550571797,
  9550562.482254831,
  9550642.97407402,
  9550743.75,
  9551079.6875,
  9551357.00724221,
  9551629.591670955,
  9551734.375,
  9552161.36489647,
  9552571.875,
  9553975.037426857,
  9555084.787061213,
  9559444.79496927,
  9565337.246534094,
  9573908.285895098,
  9577095.248650357,
  9579810.9375,
  9602998.4375,
  9631200.0,
  9631282.8125,
  9644339.681293476,
  9650043.050564112,
  9650177.913190836,
  9653409.157319397,
  9659559.375,
  9662420.490310647,
  9665681.191263411,
  9696797.25132596,
  9699216.769359216,
  9712976.5625,
  9727234.236797117,
  9731064.99977752,
  9731506.887808405,
  9737685.660092864,
  9747703.125,
  9749025.375544293,
  9749044.683595624,
  9756227.54930186,
  9770949.52215201,
  9774844.066328354,
  9775657.081608614,
  9775977.867568113,
  9776004.73781997,
  9776104.313250314,
  9777844.575861452,
  9777910.703827191,
  9778204.005161757,
  9778858.851219114,
  9780576.5625,
  9780636.776243567,
  9781240.625,
  9781306.25,
  9783660.240544371,
  9784708.035561578,
  9784973.651557056,
  9786601.34170449,
  9794216.488369154,
  9799759.375,
  9799789.696874274,
  9801318.75,
  9804746.875,
  9804774.68786805,
  9814371.73781952,
  9816131.25,
  9817379.6875,
  9821257.471137635,
  9821283.287336651,
  9830575.653223215,
  9830751.5625,
  9831013.215233142,
  9832269.387916176,
  9833201.5625,
  9833762.5,
  9835090.477278536,
  9835564.372364897,
  9836751.5625,
  9838115.678657932,
  9840328.125,
  9845176.5625,
  9846492.062683504,
  9849104.298607124,
  9852218.75,
  9853100.0,
  9855243.75,
  9868110.11708529,
  9874199.011829972,
  9891045.3125,
  9891351.34209926,
  9891484.643841166,
  9892084.24516413,
  9892211.560747504,
  9892398.28999623,
  9892785.347043065,
  9894764.118910417,
  9895080.232217988,
  9898143.75,
  9912034.181617172,
  9918207.589453574,
  9921815.468362415,
  9931720.608825784,
  9933442.071767509,
  9934359.24134096,
  9935242.1875,
  9937491.231255328,
  9939832.198448388,
  9940006.25,
  9940724.694115946,
  9941653.675227508,
  9941662.928390028,
  9956947.911210822,
  9984990.625,
  9986717.1875,
  10008379.55938739,
  10023350.0,
  10024503.066090113,
  10025137.338472044,
  10025656.429136688,
  10025744.408028519,
  10025777.28729333,
  10026291.864872677,
  10026648.4375,
  10026676.168974485,
  10026718.75,
  10026764.945529215,
  10026842.124798084,
  10026981.461079579,
  10027114.497558078,
  10029321.202449482,
  10029938.862153422,
  10030294.925983619,
  10031725.756052842,
  10031845.3125,
  10032181.25,
  10036345.618876524,
  10050852.778409323,
  10058517.014595399,
  10059911.835489884,
  10062943.75,
  10065045.913777815,
  10087254.601398285,
  10088404.580441693,
  10091015.509268532,
  10092418.520629097,
  10094149.513584454,
  10097903.57090108,
  10105518.874115529,
  10107595.081593383,
  10108841.281212363,
  10112163.5628987,
  10112171.845127827,
  10113179.6875,
  10114915.625,
  10115576.395564273,
  10115631.25,
  10115864.0625,
  10119639.0625,
  10120343.293410826,
  10120860.281938678,
  10139465.901280563,
  10142011.645172808,
  10146874.416979853,
  10147400.791614681,
  10148507.115401408,
  10148513.81805288,
  10149544.461560847,
  10153609.826628642,
  10154984.947824977,
  10176128.125,
  10176130.709966196,
  10176132.366017798,
  10193949.671064772,
  10208153.65062125,
  10267861.673072387,
  10275691.645052992,
  10286308.9560147,
  10287534.889953243,
  10603706.693600286,
  12588659.672908721,
  16404538.028302733,
  16408593.75,
  16411339.0625],
 [46.70187529791809,
  26.778519210483203,
  34.929810695380695,
  30.82141510413724,
  44.60358164593947,
  38.01323258785923,
  64.66884217258864,
  5.881532903219422,
  14.508945045254821,
  8.591476955280802,
  22.475700000774964,
  6.059010527876697,
  30.611034282837085,
  70.14888870175085,
  72.27188305490724,
  116.8866708191389,
  94.37312122572042,
  158.38593355274404,
  30.260793578755177,
  66.65105895372297,
  15.36250340168615,
  5.243377712831362,
  57.78847746656303,
  44.18555253348138,
  6.005793669874499,
  32.29867818125466,
  20.63718875758306,
  11.473883023585323,
  34.08157005468063,
  70.77641853027978,
  15.688816834759308,
  5.148328882203151,
  18.961382820747843,
  21.87668884328065,
  13.214105709573479,
  88.47133350724329,
  24.182897828988576,
  90.78455605728969,
  58.34272296964344,
  120.48726389314943,
  28.43452401417073,
  55.49341841196542,
  5.937280238167312,
  5.5839553042829735,
  39.4510554064264,
  113.46620777291861,
  27.005413788425766,
  39.7904998613743,
  16.020205783754733,
  65.63201083996525,
  5.639013552353593,
  33.076494537485466,
  63.73504673112545,
  100.01097137115934,
  15.011026434197817,
  6.499331063281637,
  55.082422375903064,
  12.348002524544372,
  65.57422597273087,
  73.7982607487158,
  105.24860805165144,
  18.325180253854814,
  13.432216601529163,
  122.94125135325994,
  5.329567590493867,
  6.8204090584222765,
  78.42636609359475,
  10.721157549866126,
  22.61850212526059,
  18.73973637693411,
  15.859704370756418,
  6.398524850085989,
  79.70370378671973,
  28.003184591389285,
  11.67332319821339,
  86.20356685221336,
  8.977083751780745,
  35.50449433951619,
  8.040301039044142,
  37.792583808466766,
  64.84517646200891,
  31.562885132014095,
  8.843211029024912,
  73.26754355115183,
  19.497453424705284,
  55.803908392947505,
  52.04814588801686,
  24.779598822289124,
  84.04060103225959,
  58.28971255007087,
  8.718814705991438,
  104.14716408288297,
  28.525893216885997,
  25.8808316667252,
  23.892057726192377,
  5.736758262867294,
  100.05019969162852,
  61.962023270583465,
  19.873583568623783,
  15.374124769516019,
  71.685059555052,
  13.650755213757346,
  6.596326083527557,
  87.59233808391186,
  93.75701709674506,
  52.63628011309282,
  125.9440198025944,
  19.667138072393843,
  60.569232496395806,
  83.19947276389935,
  13.496391256413157,
  23.053358246917433,
  11.906802185589642,
  22.288862132045832,
  59.54409555279975,
  19.103556687794065,
  10.454869328707407,
  21.286456227818473,
  60.897953504704304,
  64.29223309561432,
  11.862702552214545,
  5.808698804599091,
  83.96516192452341,
  6.816859917096688,
  53.12146184018825,
  23.09036625445928,
  24.957494487293186,
  142.2189622687855,
  176.9448703939842,
  76.62692829352429,
  61.00728137972065,
  235.04530665846653,
  12.247830850512686,
  9.171729376032681,
  29.018075632466125,
  51.464352574582975,
  5.993956786459045,
  6.689325892679059,
  76.35019276971984,
  8.620948308983042,
  215.01428540330227,
  88.55712331610883,
  99.16409201726042,
  19.81545907813914,
  54.3517845707529,
  15.825233633362279,
  105.35225873348286,
  63.42805987308003,
  21.731840707425633,
  25.34427159650101,
  9.906893852012285,
  26.943826383616134,
  15.033486848344133,
  94.79963746730922,
  33.26206419656441,
  19.2103804916421,
  121.47075273348338,
  40.648977610270634,
  83.5754857525676,
  76.05734160265169,
  48.00345733024239,
  36.54351281191254,
  183.85584201330803,
  7.192015712466861,
  6.314194553253625,
  80.49713227669442,
  27.28945603623275,
  7.918086126010027,
  6.1403621485599285,
  6.299302799259975,
  11.235049560900949,
  32.76524849824067,
  16.119081731880655,
  98.94558817675104,
  107.63240248617849,
  66.61378619660601,
  51.205806035019336,
  9.106330333642093,
  15.421851353338477,
  122.4022391184865,
  24.641696334303152,
  28.393866846438545,
  331.83909665368395,
  51.550164942476414,
  64.67215489166706,
  85.14722149533038,
  19.367273532339695,
  65.2681552959928,
  187.71955306086326,
  12.01059252865771,
  132.19415426244564,
  63.90082525530318,
  60.24878488034214,
  8.473277322168736,
  80.46228187453599,
  15.365893163035754,
  18.951752350717324,
  11.02001930919038,
  16.636538239386148,
  7.4145744059897405,
  158.1707401026462,
  118.36548762897249,
  57.078161900054596,
  5.336684648969083,
  15.523664315186855,
  79.3768881782577,
  6.664201222387136,
  32.830191941520106,
  165.0290558258976,
  42.889092279693536,
  58.16852990670368,
  58.75245768112613,
  47.202129791859754,
  9.184378021602278,
  11.362450332052997,
  40.926474840981655,
  22.77027582241642,
  263.9694528166787,
  5.7259397478116965,
  50.46132235134195,
  24.560045009209695,
  30.746574037063635,
  13.546495420504751,
  184.76117478753852,
  7.228171409420903,
  162.03123877781854,
  27.248886117258103,
  51.91949639004735,
  17.12707740789666,
  9.011278227691363,
  12.533170314865856,
  14.484976310635073,
  40.320256804183046,
  229.2243925645167,
  58.020873599503744,
  68.23593521651941,
  9.032599796589992,
  13.049140296246122,
  75.68453223371547,
  88.0263414247799,
  72.72855669296845,
  9.26632529624093,
  31.95568541010842,
  18.565288204072786,
  6.572989215706378,
  33.668481949351694,
  166.2288172988735,
  31.49381332153061,
  94.78094980666046,
  69.95620963513767,
  5.817151192690965,
  13.859823949122577,
  69.69210932208712,
  15.88621166047347,
  128.88385834261578,
  100.40517459160081,
  154.29904193148087,
  28.170930516091474,
  24.29433624785896,
  19.17730474305895,
  7.191627805762757,
  45.882471698025086,
  17.746097797005703,
  116.96305489355336,
  72.3706972161981,
  59.91569511576114,
  5.130529264647713,
  40.30946553961717,
  8.313922916278605,
  12.39843432757452,
  149.82636603202636,
  30.307824970585465,
  72.89444119752288,
  92.46323927551754,
  41.36082419237121,
  5.027807221132128,
  27.189549548308737,
  63.848055796521415,
  45.93119879337665,
  32.749421257107585,
  18.152133104070142,
  33.48537727923696,
  29.754083973826617,
  63.55459163992968,
  8.182079526561491,
  27.709249610686122,
  21.52200982709353,
  61.48813416578332,
  7.558886495050455,
  290.92460191285636,
  12.470054926683366,
  28.199942529957777,
  48.91223674312512,
  65.76840818404912,
  83.0356230666003,
  19.544457521226143,
  87.74474900680978,
  11.63586723039669,
  68.29038163460551,
  18.162090793620557,
  7.221652973018711,
  56.16122131601681,
  6.783111999289823,
  39.407042157451386,
  63.8205451197921,
  25.782035459965066,
  9.66485399243889,
  78.15853155327203,
  54.8818379950564,
  168.5959610840994,
  29.910783078630168,
  13.846890850248322,
  68.98559104847178,
  18.580597483642176,
  7.970672699358682,
  6.703225490321401,
  22.73096793591837,
  6.596405849687119,
  64.42642825755996,
  93.89526079896632,
  29.868655572523842,
  110.28123368116081,
  13.292201638715234,
  113.70526911240576,
  42.91729359508638,
  13.532323066322158,
  124.5544471192682,
  18.54759331897587,
  140.37946363020967,
  11.913690245298573,
  177.84964793066007,
  13.834385703313497,
  15.967198658668417,
  61.72913782223837,
  23.21929200300011,
  45.391326773925826,
  129.13186390674142,
  289.8578930967841,
  47.57153492728805,
  121.00351354903638,
  7.94097897408707,
  61.04952149451691,
  134.0361868713172,
  214.17243562702177,
  62.01118219847119,
  13.244468452883574,
  87.72355187342083,
  82.8408176591191,
  55.46424805223416,
  68.08558901313077,
  50.43535967618401,
  18.542564749561496,
  42.87885858443832,
  9.07695999857595,
  19.84394954097275,
  23.95951748844324,
  20.700558833604617,
  123.57815567731637,
  58.753197064301574,
  10.066847333883919,
  49.38895412278483,
  72.32167806734807,
  97.90205339136683,
  52.058015177071034,
  26.2853577981563,
  27.544337630720243,
  77.69744820677403,
  8.306668086960862,
  34.16270328747581,
  53.502427345286506,
  193.24332864560648,
  60.7331952820079,
  6.589600930168098,
  65.90143558178514,
  170.17334634892867,
  172.87848605710977,
  23.675753189865198,
  140.81252555868986,
  39.95165549837951,
  46.02447052032211,
  27.7473359302854,
  13.330172246213102,
  36.332422048138405,
  18.324056315770086,
  6.525565740453631,
  117.74747273904626,
  52.71835797388259,
  38.49040964502407,
  5.91113998552584,
  7.084609515230091,
  38.054669631797026,
  23.112362852906486,
  16.108614952900894,
  32.17757186603942,
  12.51600548626608,
  60.48537328649722,
  12.540862610173393,
  67.09730716201294,
  102.20363520730186,
  44.1482154547918,
  27.986934925557147,
  18.85263328166431,
  51.16044097513598,
  23.74790668158524,
  39.76571369114431,
  69.98954737835284,
  38.030985321506456,
  65.0733248508346,
  24.00663937459156,
  5.384592801821165,
  25.751840925971877,
  15.155411279226223,
  7.487981088188855,
  34.924594258909224,
  18.169561500631385,
  16.243358751616142,
  32.60861642361801,
  60.950109820315454,
  22.98762916282685,
  12.48566726451325,
  34.3570237280685,
  14.045674840285566,
  7.894111881548978,
  24.78786671620434,
  47.804558658359376,
  199.75240160884076,
  16.71756566968356,
  60.29502325428916,
  46.089573923116326,
  83.37064609866479,
  43.71183619906472,
  34.59568307414332,
  82.52801450426367,
  104.62665395840264,
  18.910165928715287,
  146.6039879194314,
  80.670907496638,
  8.622317475166557,
  31.226910651331885,
  28.261438151174957,
  61.12100158293882,
  73.68523527326911,
  101.43427681237512,
  5.436388966208519,
  31.15810270151968,
  67.81767839694805,
  36.088764046300724,
  68.36630109803939,
  5.083484642324268,
  35.49051144008356,
  195.52031483802747,
  44.135822555962626,
  87.54202243917635,
  35.51655732031087,
  41.84939694884301,
  78.7421538584909,
  40.89916675904822,
  47.13549197419221,
  121.40182023529263,
  24.36392402399269,
  8.553854977142008,
  6.7508260798501185,
  85.30027969648648,
  37.683851070644536,
  56.65996423617711,
  199.88335181180423,
  28.24942158257516,
  76.93091867108134,
  61.19748395650218,
  22.89215591677996,
  62.86926419825133,
  42.49217559907337,
  27.414056493984745,
  22.263883733765805,
  12.157407680701521,
  54.020758036195005,
  26.92548536086302,
  20.328716638784314,
  8.65534820567422,
  82.32073003455208,
  41.70842855353828,
  53.309874547430674,
  5.227877997793191,
  6.400479801491083,
  23.618034586963198,
  9.614009660866255,
  8.781437799388407,
  85.08405515989945,
  417.86527539383354,
  5.683346668273212,
  67.7976967204342,
  6.064131779144126,
  8.814189052367647,
  10.39708143194343,
  33.3499571716186,
  5.3715080002456705,
  69.35963564089003,
  56.33067338853214,
  8.8042783110149,
  28.68439347473634,
  248.7313731906599,
  96.36491044932919,
  6.949047935340456,
  46.55847659608011,
  12.288800014797953,
  65.54357471321505,
  50.09034354830554,
  5.3866740948357155,
  91.12887536554348,
  79.46717313107594,
  47.00454334581888,
  23.2363923581846,
  11.960971598207031,
  18.19515273658992,
  43.06355929057859,
  18.61088588298651,
  6.897914152057769,
  91.13953798748346,
  75.35457152362447,
  25.912475656617712,
  84.14866642330719,
  34.842502215644636,
  35.05145853207716,
  72.6123756367027,
  123.0584301955171,
  41.40845733110942,
  33.035603850593006,
  106.52153757221103,
  55.42619301913978,
  59.45955314783746,
  120.39101402378411,
  17.699981204886807,
  53.6524276291372,
  24.823705272529708,
  25.5723615970677,
  26.198388621020243,
  36.254440090333176,
  25.915910648464575,
  59.854136634367805,
  9.379561233948177,
  17.52313468120081,
  39.21106639586953,
  25.2153696824332,
  24.786140624698064,
  6.46281408689269,
  27.558196781412168,
  41.04480370114539,
  42.363107194378486,
  94.39277672362118,
  6.346958450260396,
  29.56123711211159,
  74.8117128443656,
  115.5617979284524,
  17.555682944633705,
  24.605936950663086,
  141.55727288859697,
  58.54508817669632,
  28.739716613864232,
  26.226336789426146,
  49.147688160152526,
  34.81408404971341,
  74.77143674226606,
  517.048739068739,
  18.226501409792487,
  8.711045803339946,
  46.97639509392978,
  214.27731171191277,
  30.349908951550802,
  29.469963141960495,
  24.802438125808365,
  124.94864884208872,
  52.409038775570245,
  37.219864730453786,
  105.9241515787472,
  82.91094895760472,
  170.55872320965108,
  101.07256238260877,
  185.1427970797945,
  40.43391062291363,
  16.253770377886774,
  10.625985547298702,
  25.522355193560234,
  6.021618595032555,
  49.49888562316846,
  29.957535500909582,
  173.88775941692262,
  203.56731118014963,
  10.951504453771573,
  24.696898632965173,
  13.363405267913505,
  6.855916618696353,
  28.812161492368105,
  31.664977251124427,
  12.33588087671406,
  55.6313817056442,
  45.469855298211236,
  72.02802761282098,
  78.30722344921077,
  9.194128811396613,
  40.65508410729425,
  64.99007470630723,
  65.56742580181779,
  8.973030365894747,
  45.04201811832388,
  55.71597256887867,
  87.86953020720199,
  132.5650753261358,
  19.158569090946646,
  72.5907387068944,
  49.393722498481225,
  18.585572143673335,
  28.667635431854634,
  100.06111359732867,
  30.905561285705293,
  10.37691647226926,
  10.494858268227146,
  13.820311339199439,
  17.4115480964118,
  6.91981521429716,
  21.079867375164344,
  104.97289321266481,
  170.21304565625908,
  7.771967007667543,
  49.66112050609631,
  60.92432429660513,
  68.1286865792631,
  91.45217897216054])

Energy spectrum

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

([6370415.625,
  7493778.0300543625,
  7555534.375,
  7962801.129695198,
  7964143.75,
  7967822.601024829,
  7968253.340446217,
  8037436.8105644435,
  8042098.319185842,
  8042109.910255405,
  8075457.405665217,
  8079616.903053936,
  8079639.0625,
  8079639.0625,
  8082554.6875,
  8083149.985643593,
  8083509.375,
  8136437.5,
  8153898.4375,
  8154428.125,
  8174020.418735455,
  8175112.4367218865,
  8175304.6875,
  8176729.6875,
  8180535.107012575,
  8185240.625,
  8186591.116974672,
  8186642.596790322,
  8192020.3125,
  8192148.420630831,
  8192218.933158315,
  8192246.128306476,
  8192264.827273551,
  8192624.537612015,
  8193069.77111207,
  8193671.875,
  8194315.650395631,
  8194335.179351896,
  8194665.625,
  8198375.904300488,
  8199085.9329081075,
  8199108.879339639,
  8199257.788063141,
  8199354.130501113,
  8199365.625,
  8199607.8125,
  8200785.855714465,
  8200825.710530285,
  8200942.4443801725,
  8202126.5625,
  8203862.9426920945,
  8203944.91640513,
  8204364.0625,
  8209717.541984143,
  8212478.143119827,
  8213118.827638851,
  8213290.625,
  8216850.661793346,
  8219637.5,
  8220700.0,
  8222567.1875,
  8225080.360513974,
  8225993.288313369,
  8227054.6875,
  8270126.197484022,
  8271447.153493197,
  8273425.394544437,
  8276015.348340758,
  8276460.626891639,
  8281905.42122556,
  8282150.221704014,
  8282841.598249388,
  8283227.390570716,
  8283833.745651554,
  8283956.573109734,
  8284219.662687637,
  8284299.548969384,
  8284312.5,
  8284615.128322482,
  8284924.827108234,
  8285814.0625,
  8286107.154223692,
  8286499.360096596,
  8287085.9375,
  8288274.273612795,
  8290093.75,
  8405582.8125,
  8408315.40149184,
  8410287.5,
  8412654.6875,
  8414519.486750592,
  8417281.25,
  8417697.575686883,
  8420295.777292263,
  8423079.033895,
  8423748.898787437,
  8426082.188130599,
  8426984.27525257,
  8427049.827288557,
  8429182.921553101,
  8429360.9375,
  8436708.989427311,
  8438760.03815258,
  8441468.41491777,
  8443260.37791533,
  8498326.495300852,
  8503645.84830694,
  8504191.482654564,
  8504353.125,
  8505589.0625,
  8505802.686041163,
  8506169.95587018,
  8506878.572144937,
  8507337.254467461,
  8510127.244231887,
  8510898.080946347,
  8514135.459702633,
  8514590.104032371,
  8514784.375,
  8516200.0,
  8516287.9243066,
  8558854.615864597,
  8558854.728845127,
  8560740.886496548,
  8599375.0,
  8620446.521481171,
  8621854.718896056,
  8621908.004352933,
  8622296.875,
  8622376.5625,
  8623254.6875,
  8623889.16336461,
  8624125.615159553,
  8624963.867425408,
  8625138.678546667,
  8625142.1875,
  8625294.355791416,
  8625304.304217637,
  8625317.418941231,
  8625453.70795514,
  8625489.014595686,
  8625871.875,
  8626309.420773359,
  8626314.34711579,
  8626330.180160021,
  8626500.214021055,
  8626987.82077895,
  8627556.25,
  8627602.648604143,
  8627869.222985258,
  8627925.610292373,
  8628051.477168456,
  8628068.977281945,
  8628273.4375,
  8628335.361533169,
  8628447.67354922,
  8628890.625,
  8628948.167414863,
  8629007.646143977,
  8629078.125,
  8629108.293836331,
  8629115.481515514,
  8629471.875,
  8629586.199304296,
  8629597.238134941,
  8629612.5,
  8629784.309322458,
  8629855.481487747,
  8629856.738010824,
  8629859.624949783,
  8629945.380105993,
  8630056.25,
  8630219.832351342,
  8630350.0,
  8630360.9375,
  8630387.5,
  8630471.875,
  8630526.283756888,
  8630539.183062648,
  8630699.555421626,
  8630931.89094182,
  8630947.674384046,
  8631130.274577262,
  8631195.3125,
  8631726.5625,
  8631736.264508057,
  8631911.756756522,
  8631978.125,
  8632055.865419207,
  8632080.724470418,
  8632136.883974709,
  8632143.75,
  8632181.658875622,
  8632219.08711897,
  8632282.8125,
  8632531.04064556,
  8632536.625366537,
  8632645.24249484,
  8632884.037684638,
  8633043.196403252,
  8633146.695263699,
  8633228.657263376,
  8633828.125,
  8634537.346366692,
  8635047.095092937,
  8635993.75,
  8636199.481363596,
  8636318.043321159,
  8636370.3125,
  8637078.125,
  8637095.709937243,
  8637131.25,
  8637354.6875,
  8637481.070846988,
  8637659.520163998,
  8637735.945668709,
  8637758.100223372,
  8638475.295057032,
  8649636.234736966,
  8709416.469863476,
  8718930.804637767,
  8720954.6875,
  8720991.766065644,
  8721909.246923601,
  8722744.140997415,
  8723401.5625,
  8727917.607537014,
  8728189.743827332,
  8728550.810509421,
  8728790.220477886,
  8728839.807495305,
  8729002.999835046,
  8729168.75,
  8729222.352610271,
  8729265.625,
  8729268.681606099,
  8729285.561634906,
  8729371.270680608,
  8729645.3125,
  8729757.8125,
  8729776.5625,
  8729881.172019469,
  8730502.842709733,
  8730694.568198716,
  8730964.589418175,
  8731267.1875,
  8731302.535127826,
  8731329.340709189,
  8731512.384490192,
  8731532.8125,
  8731535.980263544,
  8731619.996473566,
  8731639.0625,
  8731662.72758086,
  8731764.082063837,
  8731782.8125,
  8731807.677842692,
  8731807.794366304,
  8731813.80389695,
  8731840.504328666,
  8731914.32290047,
  8731964.0625,
  8731992.769423552,
  8732096.875,
  8732143.75,
  8732151.037202025,
  8732214.063766768,
  8732340.598849488,
  8732457.863855394,
  8732539.747332722,
  8732703.125,
  8732704.665772485,
  8732723.4375,
  8732745.20205004,
  8732810.9375,
  8732828.856160294,
  8732873.366202427,
  8732968.99372651,
  8732993.391307807,
  8732995.592692751,
  8733025.646593425,
  8733176.051991904,
  8733509.721739609,
  8733581.25,
  8733712.297023658,
  8733740.582411997,
  8733807.315393174,
  8734290.625,
  8734438.938187826,
  8734583.68221809,
  8734821.029561054,
  8734943.733371412,
  8735001.5625,
  8735223.4375,
  8735657.8125,
  8736065.477876388,
  8736447.690357901,
  8737406.552472765,
  8737790.625,
  8737857.08753616,
  8737931.617944706,
  8738483.967364809,
  8739176.65991593,
  8740860.925612165,
  8742864.0625,
  8774249.843960293,
  8776625.369678203,
  8777550.0,
  8780571.875,
  8782018.75,
  8783288.906315159,
  8939999.679760123,
  8940659.135945218,
  8943594.460837632,
  8943672.788193561,
  8944380.940364253,
  8945540.24073969,
  8946936.601113407,
  8947242.853862116,
  8947565.453072768,
  8947692.965163244,
  8948229.6875,
  8948246.36259576,
  8948370.3125,
  8948792.1875,
  8948849.573628077,
  8949126.031295095,
  8949397.581026677,
  8949678.125,
  8949825.450688124,
  8949933.931057876,
  8950042.961132133,
  8950165.870370278,
  8950512.5,
  8951655.907223774,
  8951675.595838109,
  8952054.458046323,
  8955907.598836772,
  8955968.75,
  8965944.842796069,
  8965960.221309869,
  8967595.721422737,
  9014307.874400143,
  9015828.214325279,
  9018188.602514023,
  9019171.344940957,
  9023698.782810353,
  9142423.4375,
  9149682.8125,
  9149684.375,
  9160189.0625,
  9160331.955291709,
  9160423.710363453,
  9161468.432930669,
  9189379.00459956,
  9190366.260583293,
  9192260.354203666,
  9192628.125,
  9192715.625,
  9192741.383149007,
  9196253.125,
  9197878.725043751,
  9198312.5,
  9202979.6875,
  9202982.688012155,
  9223553.073441105,
  9241754.6875,
  9252871.91279967,
  9257876.55002248,
  9259709.375,
  9261424.756808786,
  9262584.375,
  9262589.686682368,
  9262654.297321755,
  9263220.755100373,
  9263274.411818977,
  9263409.206069665,
  9264044.421800464,
  9264051.5625,
  9266150.0,
  9266654.646965856,
  9266663.564005904,
  9267337.222910373,
  9267752.861241207,
  9267853.366124198,
  9268108.020947196,
  9268181.25,
  9268620.3125,
  9268640.294603912,
  9268651.42046744,
  9268667.1875,
  9268970.405585008,
  9269066.05474343,
  9269464.0625,
  9269828.862440998,
  9270009.375,
  9270177.310870457,
  9270337.288028244,
  9270753.125,
  9272040.625,
  9272367.159734856,
  9272393.317499336,
  9272560.9375,
  9272563.762278028,
  9273904.403491234,
  9274265.625,
  9274285.9375,
  9276315.625,
  9276516.907272683,
  9276638.694584612,
  9277175.790397238,
  9278025.055010123,
  9278710.607178574,
  9278831.034385173,
  9283077.211784298,
  9284546.116646703,
  9285478.943681853,
  9287675.0,
  9295871.512814384,
  9377578.8490112,
  9394276.331778187,
  9402774.743281849,
  9459024.175972257,
  9471436.258415814,
  9475014.0625,
  9478514.0625,
  9479576.115839161,
  9481256.819753585,
  9482750.0,
  9484232.146311082,
  9484928.125,
  9485565.185714068,
  9498817.157990567,
  9509428.125,
  9510382.077089146,
  9512260.637481933,
  9514498.4375,
  9516859.960155206,
  9518134.375,
  9518170.909750745,
  9521739.455142671,
  9526175.0,
  9526532.8125,
  9548804.08259865,
  9549356.25,
  9549737.5,
  9550221.550571797,
  9550562.482254831,
  9550642.97407402,
  9550743.75,
  9551079.6875,
  9551357.00724221,
  9551629.591670955,
  9551734.375,
  9552161.36489647,
  9552571.875,
  9553975.037426857,
  9555084.787061213,
  9559444.79496927,
  9565337.246534094,
  9573908.285895098,
  9577095.248650357,
  9579810.9375,
  9602998.4375,
  9631200.0,
  9631282.8125,
  9644339.681293476,
  9650043.050564112,
  9650177.913190836,
  9653409.157319397,
  9659559.375,
  9662420.490310647,
  9665681.191263411,
  9696797.25132596,
  9699216.769359216,
  9712976.5625,
  9727234.236797117,
  9731064.99977752,
  9731506.887808405,
  9737685.660092864,
  9747703.125,
  9749025.375544293,
  9749044.683595624,
  9756227.54930186,
  9770949.52215201,
  9774844.066328354,
  9775657.081608614,
  9775977.867568113,
  9776004.73781997,
  9776104.313250314,
  9777844.575861452,
  9777910.703827191,
  9778204.005161757,
  9778858.851219114,
  9780576.5625,
  9780636.776243567,
  9781240.625,
  9781306.25,
  9783660.240544371,
  9784708.035561578,
  9784973.651557056,
  9786601.34170449,
  9794216.488369154,
  9799759.375,
  9799789.696874274,
  9801318.75,
  9804746.875,
  9804774.68786805,
  9814371.73781952,
  9816131.25,
  9817379.6875,
  9821257.471137635,
  9821283.287336651,
  9830575.653223215,
  9830751.5625,
  9831013.215233142,
  9832269.387916176,
  9833201.5625,
  9833762.5,
  9835090.477278536,
  9835564.372364897,
  9836751.5625,
  9838115.678657932,
  9840328.125,
  9845176.5625,
  9846492.062683504,
  9849104.298607124,
  9852218.75,
  9853100.0,
  9855243.75,
  9868110.11708529,
  9874199.011829972,
  9891045.3125,
  9891351.34209926,
  9891484.643841166,
  9892084.24516413,
  9892211.560747504,
  9892398.28999623,
  9892785.347043065,
  9894764.118910417,
  9895080.232217988,
  9898143.75,
  9912034.181617172,
  9918207.589453574,
  9921815.468362415,
  9931720.608825784,
  9933442.071767509,
  9934359.24134096,
  9935242.1875,
  9937491.231255328,
  9939832.198448388,
  9940006.25,
  9940724.694115946,
  9941653.675227508,
  9941662.928390028,
  9956947.911210822,
  9984990.625,
  9986717.1875,
  10008379.55938739,
  10023350.0,
  10024503.066090113,
  10025137.338472044,
  10025656.429136688,
  10025744.408028519,
  10025777.28729333,
  10026291.864872677,
  10026648.4375,
  10026676.168974485,
  10026718.75,
  10026764.945529215,
  10026842.124798084,
  10026981.461079579,
  10027114.497558078,
  10029321.202449482,
  10029938.862153422,
  10030294.925983619,
  10031725.756052842,
  10031845.3125,
  10032181.25,
  10036345.618876524,
  10050852.778409323,
  10058517.014595399,
  10059911.835489884,
  10062943.75,
  10065045.913777815,
  10087254.601398285,
  10088404.580441693,
  10091015.509268532,
  10092418.520629097,
  10094149.513584454,
  10097903.57090108,
  10105518.874115529,
  10107595.081593383,
  10108841.281212363,
  10112163.5628987,
  10112171.845127827,
  10113179.6875,
  10114915.625,
  10115576.395564273,
  10115631.25,
  10115864.0625,
  10119639.0625,
  10120343.293410826,
  10120860.281938678,
  10139465.901280563,
  10142011.645172808,
  10146874.416979853,
  10147400.791614681,
  10148507.115401408,
  10148513.81805288,
  10149544.461560847,
  10153609.826628642,
  10154984.947824977,
  10176128.125,
  10176130.709966196,
  10176132.366017798,
  10193949.671064772,
  10208153.65062125,
  10267861.673072387,
  10275691.645052992,
  10286308.9560147,
  10287534.889953243,
  10603706.693600286,
  12588659.672908721,
  16404538.028302733,
  16408593.75,
  16411339.0625],
 [46.70187529791809,
  26.778519210483203,
  34.929810695380695,
  30.82141510413724,
  44.60358164593947,
  38.01323258785923,
  64.66884217258864,
  5.881532903219422,
  14.508945045254821,
  8.591476955280802,
  22.475700000774964,
  6.059010527876697,
  30.611034282837085,
  70.14888870175085,
  72.27188305490724,
  116.8866708191389,
  94.37312122572042,
  158.38593355274404,
  30.260793578755177,
  66.65105895372297,
  15.36250340168615,
  5.243377712831362,
  57.78847746656303,
  44.18555253348138,
  6.005793669874499,
  32.29867818125466,
  20.63718875758306,
  11.473883023585323,
  34.08157005468063,
  70.77641853027978,
  15.688816834759308,
  5.148328882203151,
  18.961382820747843,
  21.87668884328065,
  13.214105709573479,
  88.47133350724329,
  24.182897828988576,
  90.78455605728969,
  58.34272296964344,
  120.48726389314943,
  28.43452401417073,
  55.49341841196542,
  5.937280238167312,
  5.5839553042829735,
  39.4510554064264,
  113.46620777291861,
  27.005413788425766,
  39.7904998613743,
  16.020205783754733,
  65.63201083996525,
  5.639013552353593,
  33.076494537485466,
  63.73504673112545,
  100.01097137115934,
  15.011026434197817,
  6.499331063281637,
  55.082422375903064,
  12.348002524544372,
  65.57422597273087,
  73.7982607487158,
  105.24860805165144,
  18.325180253854814,
  13.432216601529163,
  122.94125135325994,
  5.329567590493867,
  6.8204090584222765,
  78.42636609359475,
  10.721157549866126,
  22.61850212526059,
  18.73973637693411,
  15.859704370756418,
  6.398524850085989,
  79.70370378671973,
  28.003184591389285,
  11.67332319821339,
  86.20356685221336,
  8.977083751780745,
  35.50449433951619,
  8.040301039044142,
  37.792583808466766,
  64.84517646200891,
  31.562885132014095,
  8.843211029024912,
  73.26754355115183,
  19.497453424705284,
  55.803908392947505,
  52.04814588801686,
  24.779598822289124,
  84.04060103225959,
  58.28971255007087,
  8.718814705991438,
  104.14716408288297,
  28.525893216885997,
  25.8808316667252,
  23.892057726192377,
  5.736758262867294,
  100.05019969162852,
  61.962023270583465,
  19.873583568623783,
  15.374124769516019,
  71.685059555052,
  13.650755213757346,
  6.596326083527557,
  87.59233808391186,
  93.75701709674506,
  52.63628011309282,
  125.9440198025944,
  19.667138072393843,
  60.569232496395806,
  83.19947276389935,
  13.496391256413157,
  23.053358246917433,
  11.906802185589642,
  22.288862132045832,
  59.54409555279975,
  19.103556687794065,
  10.454869328707407,
  21.286456227818473,
  60.897953504704304,
  64.29223309561432,
  11.862702552214545,
  5.808698804599091,
  83.96516192452341,
  6.816859917096688,
  53.12146184018825,
  23.09036625445928,
  24.957494487293186,
  142.2189622687855,
  176.9448703939842,
  76.62692829352429,
  61.00728137972065,
  235.04530665846653,
  12.247830850512686,
  9.171729376032681,
  29.018075632466125,
  51.464352574582975,
  5.993956786459045,
  6.689325892679059,
  76.35019276971984,
  8.620948308983042,
  215.01428540330227,
  88.55712331610883,
  99.16409201726042,
  19.81545907813914,
  54.3517845707529,
  15.825233633362279,
  105.35225873348286,
  63.42805987308003,
  21.731840707425633,
  25.34427159650101,
  9.906893852012285,
  26.943826383616134,
  15.033486848344133,
  94.79963746730922,
  33.26206419656441,
  19.2103804916421,
  121.47075273348338,
  40.648977610270634,
  83.5754857525676,
  76.05734160265169,
  48.00345733024239,
  36.54351281191254,
  183.85584201330803,
  7.192015712466861,
  6.314194553253625,
  80.49713227669442,
  27.28945603623275,
  7.918086126010027,
  6.1403621485599285,
  6.299302799259975,
  11.235049560900949,
  32.76524849824067,
  16.119081731880655,
  98.94558817675104,
  107.63240248617849,
  66.61378619660601,
  51.205806035019336,
  9.106330333642093,
  15.421851353338477,
  122.4022391184865,
  24.641696334303152,
  28.393866846438545,
  331.83909665368395,
  51.550164942476414,
  64.67215489166706,
  85.14722149533038,
  19.367273532339695,
  65.2681552959928,
  187.71955306086326,
  12.01059252865771,
  132.19415426244564,
  63.90082525530318,
  60.24878488034214,
  8.473277322168736,
  80.46228187453599,
  15.365893163035754,
  18.951752350717324,
  11.02001930919038,
  16.636538239386148,
  7.4145744059897405,
  158.1707401026462,
  118.36548762897249,
  57.078161900054596,
  5.336684648969083,
  15.523664315186855,
  79.3768881782577,
  6.664201222387136,
  32.830191941520106,
  165.0290558258976,
  42.889092279693536,
  58.16852990670368,
  58.75245768112613,
  47.202129791859754,
  9.184378021602278,
  11.362450332052997,
  40.926474840981655,
  22.77027582241642,
  263.9694528166787,
  5.7259397478116965,
  50.46132235134195,
  24.560045009209695,
  30.746574037063635,
  13.546495420504751,
  184.76117478753852,
  7.228171409420903,
  162.03123877781854,
  27.248886117258103,
  51.91949639004735,
  17.12707740789666,
  9.011278227691363,
  12.533170314865856,
  14.484976310635073,
  40.320256804183046,
  229.2243925645167,
  58.020873599503744,
  68.23593521651941,
  9.032599796589992,
  13.049140296246122,
  75.68453223371547,
  88.0263414247799,
  72.72855669296845,
  9.26632529624093,
  31.95568541010842,
  18.565288204072786,
  6.572989215706378,
  33.668481949351694,
  166.2288172988735,
  31.49381332153061,
  94.78094980666046,
  69.95620963513767,
  5.817151192690965,
  13.859823949122577,
  69.69210932208712,
  15.88621166047347,
  128.88385834261578,
  100.40517459160081,
  154.29904193148087,
  28.170930516091474,
  24.29433624785896,
  19.17730474305895,
  7.191627805762757,
  45.882471698025086,
  17.746097797005703,
  116.96305489355336,
  72.3706972161981,
  59.91569511576114,
  5.130529264647713,
  40.30946553961717,
  8.313922916278605,
  12.39843432757452,
  149.82636603202636,
  30.307824970585465,
  72.89444119752288,
  92.46323927551754,
  41.36082419237121,
  5.027807221132128,
  27.189549548308737,
  63.848055796521415,
  45.93119879337665,
  32.749421257107585,
  18.152133104070142,
  33.48537727923696,
  29.754083973826617,
  63.55459163992968,
  8.182079526561491,
  27.709249610686122,
  21.52200982709353,
  61.48813416578332,
  7.558886495050455,
  290.92460191285636,
  12.470054926683366,
  28.199942529957777,
  48.91223674312512,
  65.76840818404912,
  83.0356230666003,
  19.544457521226143,
  87.74474900680978,
  11.63586723039669,
  68.29038163460551,
  18.162090793620557,
  7.221652973018711,
  56.16122131601681,
  6.783111999289823,
  39.407042157451386,
  63.8205451197921,
  25.782035459965066,
  9.66485399243889,
  78.15853155327203,
  54.8818379950564,
  168.5959610840994,
  29.910783078630168,
  13.846890850248322,
  68.98559104847178,
  18.580597483642176,
  7.970672699358682,
  6.703225490321401,
  22.73096793591837,
  6.596405849687119,
  64.42642825755996,
  93.89526079896632,
  29.868655572523842,
  110.28123368116081,
  13.292201638715234,
  113.70526911240576,
  42.91729359508638,
  13.532323066322158,
  124.5544471192682,
  18.54759331897587,
  140.37946363020967,
  11.913690245298573,
  177.84964793066007,
  13.834385703313497,
  15.967198658668417,
  61.72913782223837,
  23.21929200300011,
  45.391326773925826,
  129.13186390674142,
  289.8578930967841,
  47.57153492728805,
  121.00351354903638,
  7.94097897408707,
  61.04952149451691,
  134.0361868713172,
  214.17243562702177,
  62.01118219847119,
  13.244468452883574,
  87.72355187342083,
  82.8408176591191,
  55.46424805223416,
  68.08558901313077,
  50.43535967618401,
  18.542564749561496,
  42.87885858443832,
  9.07695999857595,
  19.84394954097275,
  23.95951748844324,
  20.700558833604617,
  123.57815567731637,
  58.753197064301574,
  10.066847333883919,
  49.38895412278483,
  72.32167806734807,
  97.90205339136683,
  52.058015177071034,
  26.2853577981563,
  27.544337630720243,
  77.69744820677403,
  8.306668086960862,
  34.16270328747581,
  53.502427345286506,
  193.24332864560648,
  60.7331952820079,
  6.589600930168098,
  65.90143558178514,
  170.17334634892867,
  172.87848605710977,
  23.675753189865198,
  140.81252555868986,
  39.95165549837951,
  46.02447052032211,
  27.7473359302854,
  13.330172246213102,
  36.332422048138405,
  18.324056315770086,
  6.525565740453631,
  117.74747273904626,
  52.71835797388259,
  38.49040964502407,
  5.91113998552584,
  7.084609515230091,
  38.054669631797026,
  23.112362852906486,
  16.108614952900894,
  32.17757186603942,
  12.51600548626608,
  60.48537328649722,
  12.540862610173393,
  67.09730716201294,
  102.20363520730186,
  44.1482154547918,
  27.986934925557147,
  18.85263328166431,
  51.16044097513598,
  23.74790668158524,
  39.76571369114431,
  69.98954737835284,
  38.030985321506456,
  65.0733248508346,
  24.00663937459156,
  5.384592801821165,
  25.751840925971877,
  15.155411279226223,
  7.487981088188855,
  34.924594258909224,
  18.169561500631385,
  16.243358751616142,
  32.60861642361801,
  60.950109820315454,
  22.98762916282685,
  12.48566726451325,
  34.3570237280685,
  14.045674840285566,
  7.894111881548978,
  24.78786671620434,
  47.804558658359376,
  199.75240160884076,
  16.71756566968356,
  60.29502325428916,
  46.089573923116326,
  83.37064609866479,
  43.71183619906472,
  34.59568307414332,
  82.52801450426367,
  104.62665395840264,
  18.910165928715287,
  146.6039879194314,
  80.670907496638,
  8.622317475166557,
  31.226910651331885,
  28.261438151174957,
  61.12100158293882,
  73.68523527326911,
  101.43427681237512,
  5.436388966208519,
  31.15810270151968,
  67.81767839694805,
  36.088764046300724,
  68.36630109803939,
  5.083484642324268,
  35.49051144008356,
  195.52031483802747,
  44.135822555962626,
  87.54202243917635,
  35.51655732031087,
  41.84939694884301,
  78.7421538584909,
  40.89916675904822,
  47.13549197419221,
  121.40182023529263,
  24.36392402399269,
  8.553854977142008,
  6.7508260798501185,
  85.30027969648648,
  37.683851070644536,
  56.65996423617711,
  199.88335181180423,
  28.24942158257516,
  76.93091867108134,
  61.19748395650218,
  22.89215591677996,
  62.86926419825133,
  42.49217559907337,
  27.414056493984745,
  22.263883733765805,
  12.157407680701521,
  54.020758036195005,
  26.92548536086302,
  20.328716638784314,
  8.65534820567422,
  82.32073003455208,
  41.70842855353828,
  53.309874547430674,
  5.227877997793191,
  6.400479801491083,
  23.618034586963198,
  9.614009660866255,
  8.781437799388407,
  85.08405515989945,
  417.86527539383354,
  5.683346668273212,
  67.7976967204342,
  6.064131779144126,
  8.814189052367647,
  10.39708143194343,
  33.3499571716186,
  5.3715080002456705,
  69.35963564089003,
  56.33067338853214,
  8.8042783110149,
  28.68439347473634,
  248.7313731906599,
  96.36491044932919,
  6.949047935340456,
  46.55847659608011,
  12.288800014797953,
  65.54357471321505,
  50.09034354830554,
  5.3866740948357155,
  91.12887536554348,
  79.46717313107594,
  47.00454334581888,
  23.2363923581846,
  11.960971598207031,
  18.19515273658992,
  43.06355929057859,
  18.61088588298651,
  6.897914152057769,
  91.13953798748346,
  75.35457152362447,
  25.912475656617712,
  84.14866642330719,
  34.842502215644636,
  35.05145853207716,
  72.6123756367027,
  123.0584301955171,
  41.40845733110942,
  33.035603850593006,
  106.52153757221103,
  55.42619301913978,
  59.45955314783746,
  120.39101402378411,
  17.699981204886807,
  53.6524276291372,
  24.823705272529708,
  25.5723615970677,
  26.198388621020243,
  36.254440090333176,
  25.915910648464575,
  59.854136634367805,
  9.379561233948177,
  17.52313468120081,
  39.21106639586953,
  25.2153696824332,
  24.786140624698064,
  6.46281408689269,
  27.558196781412168,
  41.04480370114539,
  42.363107194378486,
  94.39277672362118,
  6.346958450260396,
  29.56123711211159,
  74.8117128443656,
  115.5617979284524,
  17.555682944633705,
  24.605936950663086,
  141.55727288859697,
  58.54508817669632,
  28.739716613864232,
  26.226336789426146,
  49.147688160152526,
  34.81408404971341,
  74.77143674226606,
  517.048739068739,
  18.226501409792487,
  8.711045803339946,
  46.97639509392978,
  214.27731171191277,
  30.349908951550802,
  29.469963141960495,
  24.802438125808365,
  124.94864884208872,
  52.409038775570245,
  37.219864730453786,
  105.9241515787472,
  82.91094895760472,
  170.55872320965108,
  101.07256238260877,
  185.1427970797945,
  40.43391062291363,
  16.253770377886774,
  10.625985547298702,
  25.522355193560234,
  6.021618595032555,
  49.49888562316846,
  29.957535500909582,
  173.88775941692262,
  203.56731118014963,
  10.951504453771573,
  24.696898632965173,
  13.363405267913505,
  6.855916618696353,
  28.812161492368105,
  31.664977251124427,
  12.33588087671406,
  55.6313817056442,
  45.469855298211236,
  72.02802761282098,
  78.30722344921077,
  9.194128811396613,
  40.65508410729425,
  64.99007470630723,
  65.56742580181779,
  8.973030365894747,
  45.04201811832388,
  55.71597256887867,
  87.86953020720199,
  132.5650753261358,
  19.158569090946646,
  72.5907387068944,
  49.393722498481225,
  18.585572143673335,
  28.667635431854634,
  100.06111359732867,
  30.905561285705293,
  10.37691647226926,
  10.494858268227146,
  13.820311339199439,
  17.4115480964118,
  6.91981521429716,
  21.079867375164344,
  104.97289321266481,
  170.21304565625908,
  7.771967007667543,
  49.66112050609631,
  60.92432429660513,
  68.1286865792631,
  91.45217897216054])

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

([6370415.625,
  7493778.0300543625,
  7555534.375,
  7962801.129695198,
  7964143.75,
  7967822.601024829,
  7968253.340446217,
  8037436.8105644435,
  8042098.319185842,
  8042109.910255405,
  8075457.405665217,
  8079616.903053936,
  8079639.0625,
  8079639.0625,
  8082554.6875,
  8083149.985643593,
  8083509.375,
  8136437.5,
  8153898.4375,
  8154428.125,
  8174020.418735455,
  8175112.4367218865,
  8175304.6875,
  8176729.6875,
  8180535.107012575,
  8185240.625,
  8186591.116974672,
  8186642.596790322,
  8192020.3125,
  8192148.420630831,
  8192218.933158315,
  8192246.128306476,
  8192264.827273551,
  8192624.537612015,
  8193069.77111207,
  8193671.875,
  8194315.650395631,
  8194335.179351896,
  8194665.625,
  8198375.904300488,
  8199085.9329081075,
  8199108.879339639,
  8199257.788063141,
  8199354.130501113,
  8199365.625,
  8199607.8125,
  8200785.855714465,
  8200825.710530285,
  8200942.4443801725,
  8202126.5625,
  8203862.9426920945,
  8203944.91640513,
  8204364.0625,
  8209717.541984143,
  8212478.143119827,
  8213118.827638851,
  8213290.625,
  8216850.661793346,
  8219637.5,
  8220700.0,
  8222567.1875,
  8225080.360513974,
  8225993.288313369,
  8227054.6875,
  8270126.197484022,
  8271447.153493197,
  8273425.394544437,
  8276015.348340758,
  8276460.626891639,
  8281905.42122556,
  8282150.221704014,
  8282841.598249388,
  8283227.390570716,
  8283833.745651554,
  8283956.573109734,
  8284219.662687637,
  8284299.548969384,
  8284312.5,
  8284615.128322482,
  8284924.827108234,
  8285814.0625,
  8286107.154223692,
  8286499.360096596,
  8287085.9375,
  8288274.273612795,
  8290093.75,
  8405582.8125,
  8408315.40149184,
  8410287.5,
  8412654.6875,
  8414519.486750592,
  8417281.25,
  8417697.575686883,
  8420295.777292263,
  8423079.033895,
  8423748.898787437,
  8426082.188130599,
  8426984.27525257,
  8427049.827288557,
  8429182.921553101,
  8429360.9375,
  8436708.989427311,
  8438760.03815258,
  8441468.41491777,
  8443260.37791533,
  8498326.495300852,
  8503645.84830694,
  8504191.482654564,
  8504353.125,
  8505589.0625,
  8505802.686041163,
  8506169.95587018,
  8506878.572144937,
  8507337.254467461,
  8510127.244231887,
  8510898.080946347,
  8514135.459702633,
  8514590.104032371,
  8514784.375,
  8516200.0,
  8516287.9243066,
  8558854.615864597,
  8558854.728845127,
  8560740.886496548,
  8599375.0,
  8620446.521481171,
  8621854.718896056,
  8621908.004352933,
  8622296.875,
  8622376.5625,
  8623254.6875,
  8623889.16336461,
  8624125.615159553,
  8624963.867425408,
  8625138.678546667,
  8625142.1875,
  8625294.355791416,
  8625304.304217637,
  8625317.418941231,
  8625453.70795514,
  8625489.014595686,
  8625871.875,
  8626309.420773359,
  8626314.34711579,
  8626330.180160021,
  8626500.214021055,
  8626987.82077895,
  8627556.25,
  8627602.648604143,
  8627869.222985258,
  8627925.610292373,
  8628051.477168456,
  8628068.977281945,
  8628273.4375,
  8628335.361533169,
  8628447.67354922,
  8628890.625,
  8628948.167414863,
  8629007.646143977,
  8629078.125,
  8629108.293836331,
  8629115.481515514,
  8629471.875,
  8629586.199304296,
  8629597.238134941,
  8629612.5,
  8629784.309322458,
  8629855.481487747,
  8629856.738010824,
  8629859.624949783,
  8629945.380105993,
  8630056.25,
  8630219.832351342,
  8630350.0,
  8630360.9375,
  8630387.5,
  8630471.875,
  8630526.283756888,
  8630539.183062648,
  8630699.555421626,
  8630931.89094182,
  8630947.674384046,
  8631130.274577262,
  8631195.3125,
  8631726.5625,
  8631736.264508057,
  8631911.756756522,
  8631978.125,
  8632055.865419207,
  8632080.724470418,
  8632136.883974709,
  8632143.75,
  8632181.658875622,
  8632219.08711897,
  8632282.8125,
  8632531.04064556,
  8632536.625366537,
  8632645.24249484,
  8632884.037684638,
  8633043.196403252,
  8633146.695263699,
  8633228.657263376,
  8633828.125,
  8634537.346366692,
  8635047.095092937,
  8635993.75,
  8636199.481363596,
  8636318.043321159,
  8636370.3125,
  8637078.125,
  8637095.709937243,
  8637131.25,
  8637354.6875,
  8637481.070846988,
  8637659.520163998,
  8637735.945668709,
  8637758.100223372,
  8638475.295057032,
  8649636.234736966,
  8709416.469863476,
  8718930.804637767,
  8720954.6875,
  8720991.766065644,
  8721909.246923601,
  8722744.140997415,
  8723401.5625,
  8727917.607537014,
  8728189.743827332,
  8728550.810509421,
  8728790.220477886,
  8728839.807495305,
  8729002.999835046,
  8729168.75,
  8729222.352610271,
  8729265.625,
  8729268.681606099,
  8729285.561634906,
  8729371.270680608,
  8729645.3125,
  8729757.8125,
  8729776.5625,
  8729881.172019469,
  8730502.842709733,
  8730694.568198716,
  8730964.589418175,
  8731267.1875,
  8731302.535127826,
  8731329.340709189,
  8731512.384490192,
  8731532.8125,
  8731535.980263544,
  8731619.996473566,
  8731639.0625,
  8731662.72758086,
  8731764.082063837,
  8731782.8125,
  8731807.677842692,
  8731807.794366304,
  8731813.80389695,
  8731840.504328666,
  8731914.32290047,
  8731964.0625,
  8731992.769423552,
  8732096.875,
  8732143.75,
  8732151.037202025,
  8732214.063766768,
  8732340.598849488,
  8732457.863855394,
  8732539.747332722,
  8732703.125,
  8732704.665772485,
  8732723.4375,
  8732745.20205004,
  8732810.9375,
  8732828.856160294,
  8732873.366202427,
  8732968.99372651,
  8732993.391307807,
  8732995.592692751,
  8733025.646593425,
  8733176.051991904,
  8733509.721739609,
  8733581.25,
  8733712.297023658,
  8733740.582411997,
  8733807.315393174,
  8734290.625,
  8734438.938187826,
  8734583.68221809,
  8734821.029561054,
  8734943.733371412,
  8735001.5625,
  8735223.4375,
  8735657.8125,
  8736065.477876388,
  8736447.690357901,
  8737406.552472765,
  8737790.625,
  8737857.08753616,
  8737931.617944706,
  8738483.967364809,
  8739176.65991593,
  8740860.925612165,
  8742864.0625,
  8774249.843960293,
  8776625.369678203,
  8777550.0,
  8780571.875,
  8782018.75,
  8783288.906315159,
  8939999.679760123,
  8940659.135945218,
  8943594.460837632,
  8943672.788193561,
  8944380.940364253,
  8945540.24073969,
  8946936.601113407,
  8947242.853862116,
  8947565.453072768,
  8947692.965163244,
  8948229.6875,
  8948246.36259576,
  8948370.3125,
  8948792.1875,
  8948849.573628077,
  8949126.031295095,
  8949397.581026677,
  8949678.125,
  8949825.450688124,
  8949933.931057876,
  8950042.961132133,
  8950165.870370278,
  8950512.5,
  8951655.907223774,
  8951675.595838109,
  8952054.458046323,
  8955907.598836772,
  8955968.75,
  8965944.842796069,
  8965960.221309869,
  8967595.721422737,
  9014307.874400143,
  9015828.214325279,
  9018188.602514023,
  9019171.344940957,
  9023698.782810353,
  9142423.4375,
  9149682.8125,
  9149684.375,
  9160189.0625,
  9160331.955291709,
  9160423.710363453,
  9161468.432930669,
  9189379.00459956,
  9190366.260583293,
  9192260.354203666,
  9192628.125,
  9192715.625,
  9192741.383149007,
  9196253.125,
  9197878.725043751,
  9198312.5,
  9202979.6875,
  9202982.688012155,
  9223553.073441105,
  9241754.6875,
  9252871.91279967,
  9257876.55002248,
  9259709.375,
  9261424.756808786,
  9262584.375,
  9262589.686682368,
  9262654.297321755,
  9263220.755100373,
  9263274.411818977,
  9263409.206069665,
  9264044.421800464,
  9264051.5625,
  9266150.0,
  9266654.646965856,
  9266663.564005904,
  9267337.222910373,
  9267752.861241207,
  9267853.366124198,
  9268108.020947196,
  9268181.25,
  9268620.3125,
  9268640.294603912,
  9268651.42046744,
  9268667.1875,
  9268970.405585008,
  9269066.05474343,
  9269464.0625,
  9269828.862440998,
  9270009.375,
  9270177.310870457,
  9270337.288028244,
  9270753.125,
  9272040.625,
  9272367.159734856,
  9272393.317499336,
  9272560.9375,
  9272563.762278028,
  9273904.403491234,
  9274265.625,
  9274285.9375,
  9276315.625,
  9276516.907272683,
  9276638.694584612,
  9277175.790397238,
  9278025.055010123,
  9278710.607178574,
  9278831.034385173,
  9283077.211784298,
  9284546.116646703,
  9285478.943681853,
  9287675.0,
  9295871.512814384,
  9377578.8490112,
  9394276.331778187,
  9402774.743281849,
  9459024.175972257,
  9471436.258415814,
  9475014.0625,
  9478514.0625,
  9479576.115839161,
  9481256.819753585,
  9482750.0,
  9484232.146311082,
  9484928.125,
  9485565.185714068,
  9498817.157990567,
  9509428.125,
  9510382.077089146,
  9512260.637481933,
  9514498.4375,
  9516859.960155206,
  9518134.375,
  9518170.909750745,
  9521739.455142671,
  9526175.0,
  9526532.8125,
  9548804.08259865,
  9549356.25,
  9549737.5,
  9550221.550571797,
  9550562.482254831,
  9550642.97407402,
  9550743.75,
  9551079.6875,
  9551357.00724221,
  9551629.591670955,
  9551734.375,
  9552161.36489647,
  9552571.875,
  9553975.037426857,
  9555084.787061213,
  9559444.79496927,
  9565337.246534094,
  9573908.285895098,
  9577095.248650357,
  9579810.9375,
  9602998.4375,
  9631200.0,
  9631282.8125,
  9644339.681293476,
  9650043.050564112,
  9650177.913190836,
  9653409.157319397,
  9659559.375,
  9662420.490310647,
  9665681.191263411,
  9696797.25132596,
  9699216.769359216,
  9712976.5625,
  9727234.236797117,
  9731064.99977752,
  9731506.887808405,
  9737685.660092864,
  9747703.125,
  9749025.375544293,
  9749044.683595624,
  9756227.54930186,
  9770949.52215201,
  9774844.066328354,
  9775657.081608614,
  9775977.867568113,
  9776004.73781997,
  9776104.313250314,
  9777844.575861452,
  9777910.703827191,
  9778204.005161757,
  9778858.851219114,
  9780576.5625,
  9780636.776243567,
  9781240.625,
  9781306.25,
  9783660.240544371,
  9784708.035561578,
  9784973.651557056,
  9786601.34170449,
  9794216.488369154,
  9799759.375,
  9799789.696874274,
  9801318.75,
  9804746.875,
  9804774.68786805,
  9814371.73781952,
  9816131.25,
  9817379.6875,
  9821257.471137635,
  9821283.287336651,
  9830575.653223215,
  9830751.5625,
  9831013.215233142,
  9832269.387916176,
  9833201.5625,
  9833762.5,
  9835090.477278536,
  9835564.372364897,
  9836751.5625,
  9838115.678657932,
  9840328.125,
  9845176.5625,
  9846492.062683504,
  9849104.298607124,
  9852218.75,
  9853100.0,
  9855243.75,
  9868110.11708529,
  9874199.011829972,
  9891045.3125,
  9891351.34209926,
  9891484.643841166,
  9892084.24516413,
  9892211.560747504,
  9892398.28999623,
  9892785.347043065,
  9894764.118910417,
  9895080.232217988,
  9898143.75,
  9912034.181617172,
  9918207.589453574,
  9921815.468362415,
  9931720.608825784,
  9933442.071767509,
  9934359.24134096,
  9935242.1875,
  9937491.231255328,
  9939832.198448388,
  9940006.25,
  9940724.694115946,
  9941653.675227508,
  9941662.928390028,
  9956947.911210822,
  9984990.625,
  9986717.1875,
  10008379.55938739,
  10023350.0,
  10024503.066090113,
  10025137.338472044,
  10025656.429136688,
  10025744.408028519,
  10025777.28729333,
  10026291.864872677,
  10026648.4375,
  10026676.168974485,
  10026718.75,
  10026764.945529215,
  10026842.124798084,
  10026981.461079579,
  10027114.497558078,
  10029321.202449482,
  10029938.862153422,
  10030294.925983619,
  10031725.756052842,
  10031845.3125,
  10032181.25,
  10036345.618876524,
  10050852.778409323,
  10058517.014595399,
  10059911.835489884,
  10062943.75,
  10065045.913777815,
  10087254.601398285,
  10088404.580441693,
  10091015.509268532,
  10092418.520629097,
  10094149.513584454,
  10097903.57090108,
  10105518.874115529,
  10107595.081593383,
  10108841.281212363,
  10112163.5628987,
  10112171.845127827,
  10113179.6875,
  10114915.625,
  10115576.395564273,
  10115631.25,
  10115864.0625,
  10119639.0625,
  10120343.293410826,
  10120860.281938678,
  10139465.901280563,
  10142011.645172808,
  10146874.416979853,
  10147400.791614681,
  10148507.115401408,
  10148513.81805288,
  10149544.461560847,
  10153609.826628642,
  10154984.947824977,
  10176128.125,
  10176130.709966196,
  10176132.366017798,
  10193949.671064772,
  10208153.65062125,
  10267861.673072387,
  10275691.645052992,
  10286308.9560147,
  10287534.889953243,
  10603706.693600286,
  12588659.672908721,
  16404538.028302733,
  16408593.75,
  16411339.0625],
 [46.70187529791809,
  26.778519210483203,
  34.929810695380695,
  30.82141510413724,
  44.60358164593947,
  38.01323258785923,
  64.66884217258864,
  5.881532903219422,
  14.508945045254821,
  8.591476955280802,
  22.475700000774964,
  6.059010527876697,
  30.611034282837085,
  70.14888870175085,
  72.27188305490724,
  116.8866708191389,
  94.37312122572042,
  158.38593355274404,
  30.260793578755177,
  66.65105895372297,
  15.36250340168615,
  5.243377712831362,
  57.78847746656303,
  44.18555253348138,
  6.005793669874499,
  32.29867818125466,
  20.63718875758306,
  11.473883023585323,
  34.08157005468063,
  70.77641853027978,
  15.688816834759308,
  5.148328882203151,
  18.961382820747843,
  21.87668884328065,
  13.214105709573479,
  88.47133350724329,
  24.182897828988576,
  90.78455605728969,
  58.34272296964344,
  120.48726389314943,
  28.43452401417073,
  55.49341841196542,
  5.937280238167312,
  5.5839553042829735,
  39.4510554064264,
  113.46620777291861,
  27.005413788425766,
  39.7904998613743,
  16.020205783754733,
  65.63201083996525,
  5.639013552353593,
  33.076494537485466,
  63.73504673112545,
  100.01097137115934,
  15.011026434197817,
  6.499331063281637,
  55.082422375903064,
  12.348002524544372,
  65.57422597273087,
  73.7982607487158,
  105.24860805165144,
  18.325180253854814,
  13.432216601529163,
  122.94125135325994,
  5.329567590493867,
  6.8204090584222765,
  78.42636609359475,
  10.721157549866126,
  22.61850212526059,
  18.73973637693411,
  15.859704370756418,
  6.398524850085989,
  79.70370378671973,
  28.003184591389285,
  11.67332319821339,
  86.20356685221336,
  8.977083751780745,
  35.50449433951619,
  8.040301039044142,
  37.792583808466766,
  64.84517646200891,
  31.562885132014095,
  8.843211029024912,
  73.26754355115183,
  19.497453424705284,
  55.803908392947505,
  52.04814588801686,
  24.779598822289124,
  84.04060103225959,
  58.28971255007087,
  8.718814705991438,
  104.14716408288297,
  28.525893216885997,
  25.8808316667252,
  23.892057726192377,
  5.736758262867294,
  100.05019969162852,
  61.962023270583465,
  19.873583568623783,
  15.374124769516019,
  71.685059555052,
  13.650755213757346,
  6.596326083527557,
  87.59233808391186,
  93.75701709674506,
  52.63628011309282,
  125.9440198025944,
  19.667138072393843,
  60.569232496395806,
  83.19947276389935,
  13.496391256413157,
  23.053358246917433,
  11.906802185589642,
  22.288862132045832,
  59.54409555279975,
  19.103556687794065,
  10.454869328707407,
  21.286456227818473,
  60.897953504704304,
  64.29223309561432,
  11.862702552214545,
  5.808698804599091,
  83.96516192452341,
  6.816859917096688,
  53.12146184018825,
  23.09036625445928,
  24.957494487293186,
  142.2189622687855,
  176.9448703939842,
  76.62692829352429,
  61.00728137972065,
  235.04530665846653,
  12.247830850512686,
  9.171729376032681,
  29.018075632466125,
  51.464352574582975,
  5.993956786459045,
  6.689325892679059,
  76.35019276971984,
  8.620948308983042,
  215.01428540330227,
  88.55712331610883,
  99.16409201726042,
  19.81545907813914,
  54.3517845707529,
  15.825233633362279,
  105.35225873348286,
  63.42805987308003,
  21.731840707425633,
  25.34427159650101,
  9.906893852012285,
  26.943826383616134,
  15.033486848344133,
  94.79963746730922,
  33.26206419656441,
  19.2103804916421,
  121.47075273348338,
  40.648977610270634,
  83.5754857525676,
  76.05734160265169,
  48.00345733024239,
  36.54351281191254,
  183.85584201330803,
  7.192015712466861,
  6.314194553253625,
  80.49713227669442,
  27.28945603623275,
  7.918086126010027,
  6.1403621485599285,
  6.299302799259975,
  11.235049560900949,
  32.76524849824067,
  16.119081731880655,
  98.94558817675104,
  107.63240248617849,
  66.61378619660601,
  51.205806035019336,
  9.106330333642093,
  15.421851353338477,
  122.4022391184865,
  24.641696334303152,
  28.393866846438545,
  331.83909665368395,
  51.550164942476414,
  64.67215489166706,
  85.14722149533038,
  19.367273532339695,
  65.2681552959928,
  187.71955306086326,
  12.01059252865771,
  132.19415426244564,
  63.90082525530318,
  60.24878488034214,
  8.473277322168736,
  80.46228187453599,
  15.365893163035754,
  18.951752350717324,
  11.02001930919038,
  16.636538239386148,
  7.4145744059897405,
  158.1707401026462,
  118.36548762897249,
  57.078161900054596,
  5.336684648969083,
  15.523664315186855,
  79.3768881782577,
  6.664201222387136,
  32.830191941520106,
  165.0290558258976,
  42.889092279693536,
  58.16852990670368,
  58.75245768112613,
  47.202129791859754,
  9.184378021602278,
  11.362450332052997,
  40.926474840981655,
  22.77027582241642,
  263.9694528166787,
  5.7259397478116965,
  50.46132235134195,
  24.560045009209695,
  30.746574037063635,
  13.546495420504751,
  184.76117478753852,
  7.228171409420903,
  162.03123877781854,
  27.248886117258103,
  51.91949639004735,
  17.12707740789666,
  9.011278227691363,
  12.533170314865856,
  14.484976310635073,
  40.320256804183046,
  229.2243925645167,
  58.020873599503744,
  68.23593521651941,
  9.032599796589992,
  13.049140296246122,
  75.68453223371547,
  88.0263414247799,
  72.72855669296845,
  9.26632529624093,
  31.95568541010842,
  18.565288204072786,
  6.572989215706378,
  33.668481949351694,
  166.2288172988735,
  31.49381332153061,
  94.78094980666046,
  69.95620963513767,
  5.817151192690965,
  13.859823949122577,
  69.69210932208712,
  15.88621166047347,
  128.88385834261578,
  100.40517459160081,
  154.29904193148087,
  28.170930516091474,
  24.29433624785896,
  19.17730474305895,
  7.191627805762757,
  45.882471698025086,
  17.746097797005703,
  116.96305489355336,
  72.3706972161981,
  59.91569511576114,
  5.130529264647713,
  40.30946553961717,
  8.313922916278605,
  12.39843432757452,
  149.82636603202636,
  30.307824970585465,
  72.89444119752288,
  92.46323927551754,
  41.36082419237121,
  5.027807221132128,
  27.189549548308737,
  63.848055796521415,
  45.93119879337665,
  32.749421257107585,
  18.152133104070142,
  33.48537727923696,
  29.754083973826617,
  63.55459163992968,
  8.182079526561491,
  27.709249610686122,
  21.52200982709353,
  61.48813416578332,
  7.558886495050455,
  290.92460191285636,
  12.470054926683366,
  28.199942529957777,
  48.91223674312512,
  65.76840818404912,
  83.0356230666003,
  19.544457521226143,
  87.74474900680978,
  11.63586723039669,
  68.29038163460551,
  18.162090793620557,
  7.221652973018711,
  56.16122131601681,
  6.783111999289823,
  39.407042157451386,
  63.8205451197921,
  25.782035459965066,
  9.66485399243889,
  78.15853155327203,
  54.8818379950564,
  168.5959610840994,
  29.910783078630168,
  13.846890850248322,
  68.98559104847178,
  18.580597483642176,
  7.970672699358682,
  6.703225490321401,
  22.73096793591837,
  6.596405849687119,
  64.42642825755996,
  93.89526079896632,
  29.868655572523842,
  110.28123368116081,
  13.292201638715234,
  113.70526911240576,
  42.91729359508638,
  13.532323066322158,
  124.5544471192682,
  18.54759331897587,
  140.37946363020967,
  11.913690245298573,
  177.84964793066007,
  13.834385703313497,
  15.967198658668417,
  61.72913782223837,
  23.21929200300011,
  45.391326773925826,
  129.13186390674142,
  289.8578930967841,
  47.57153492728805,
  121.00351354903638,
  7.94097897408707,
  61.04952149451691,
  134.0361868713172,
  214.17243562702177,
  62.01118219847119,
  13.244468452883574,
  87.72355187342083,
  82.8408176591191,
  55.46424805223416,
  68.08558901313077,
  50.43535967618401,
  18.542564749561496,
  42.87885858443832,
  9.07695999857595,
  19.84394954097275,
  23.95951748844324,
  20.700558833604617,
  123.57815567731637,
  58.753197064301574,
  10.066847333883919,
  49.38895412278483,
  72.32167806734807,
  97.90205339136683,
  52.058015177071034,
  26.2853577981563,
  27.544337630720243,
  77.69744820677403,
  8.306668086960862,
  34.16270328747581,
  53.502427345286506,
  193.24332864560648,
  60.7331952820079,
  6.589600930168098,
  65.90143558178514,
  170.17334634892867,
  172.87848605710977,
  23.675753189865198,
  140.81252555868986,
  39.95165549837951,
  46.02447052032211,
  27.7473359302854,
  13.330172246213102,
  36.332422048138405,
  18.324056315770086,
  6.525565740453631,
  117.74747273904626,
  52.71835797388259,
  38.49040964502407,
  5.91113998552584,
  7.084609515230091,
  38.054669631797026,
  23.112362852906486,
  16.108614952900894,
  32.17757186603942,
  12.51600548626608,
  60.48537328649722,
  12.540862610173393,
  67.09730716201294,
  102.20363520730186,
  44.1482154547918,
  27.986934925557147,
  18.85263328166431,
  51.16044097513598,
  23.74790668158524,
  39.76571369114431,
  69.98954737835284,
  38.030985321506456,
  65.0733248508346,
  24.00663937459156,
  5.384592801821165,
  25.751840925971877,
  15.155411279226223,
  7.487981088188855,
  34.924594258909224,
  18.169561500631385,
  16.243358751616142,
  32.60861642361801,
  60.950109820315454,
  22.98762916282685,
  12.48566726451325,
  34.3570237280685,
  14.045674840285566,
  7.894111881548978,
  24.78786671620434,
  47.804558658359376,
  199.75240160884076,
  16.71756566968356,
  60.29502325428916,
  46.089573923116326,
  83.37064609866479,
  43.71183619906472,
  34.59568307414332,
  82.52801450426367,
  104.62665395840264,
  18.910165928715287,
  146.6039879194314,
  80.670907496638,
  8.622317475166557,
  31.226910651331885,
  28.261438151174957,
  61.12100158293882,
  73.68523527326911,
  101.43427681237512,
  5.436388966208519,
  31.15810270151968,
  67.81767839694805,
  36.088764046300724,
  68.36630109803939,
  5.083484642324268,
  35.49051144008356,
  195.52031483802747,
  44.135822555962626,
  87.54202243917635,
  35.51655732031087,
  41.84939694884301,
  78.7421538584909,
  40.89916675904822,
  47.13549197419221,
  121.40182023529263,
  24.36392402399269,
  8.553854977142008,
  6.7508260798501185,
  85.30027969648648,
  37.683851070644536,
  56.65996423617711,
  199.88335181180423,
  28.24942158257516,
  76.93091867108134,
  61.19748395650218,
  22.89215591677996,
  62.86926419825133,
  42.49217559907337,
  27.414056493984745,
  22.263883733765805,
  12.157407680701521,
  54.020758036195005,
  26.92548536086302,
  20.328716638784314,
  8.65534820567422,
  82.32073003455208,
  41.70842855353828,
  53.309874547430674,
  5.227877997793191,
  6.400479801491083,
  23.618034586963198,
  9.614009660866255,
  8.781437799388407,
  85.08405515989945,
  417.86527539383354,
  5.683346668273212,
  67.7976967204342,
  6.064131779144126,
  8.814189052367647,
  10.39708143194343,
  33.3499571716186,
  5.3715080002456705,
  69.35963564089003,
  56.33067338853214,
  8.8042783110149,
  28.68439347473634,
  248.7313731906599,
  96.36491044932919,
  6.949047935340456,
  46.55847659608011,
  12.288800014797953,
  65.54357471321505,
  50.09034354830554,
  5.3866740948357155,
  91.12887536554348,
  79.46717313107594,
  47.00454334581888,
  23.2363923581846,
  11.960971598207031,
  18.19515273658992,
  43.06355929057859,
  18.61088588298651,
  6.897914152057769,
  91.13953798748346,
  75.35457152362447,
  25.912475656617712,
  84.14866642330719,
  34.842502215644636,
  35.05145853207716,
  72.6123756367027,
  123.0584301955171,
  41.40845733110942,
  33.035603850593006,
  106.52153757221103,
  55.42619301913978,
  59.45955314783746,
  120.39101402378411,
  17.699981204886807,
  53.6524276291372,
  24.823705272529708,
  25.5723615970677,
  26.198388621020243,
  36.254440090333176,
  25.915910648464575,
  59.854136634367805,
  9.379561233948177,
  17.52313468120081,
  39.21106639586953,
  25.2153696824332,
  24.786140624698064,
  6.46281408689269,
  27.558196781412168,
  41.04480370114539,
  42.363107194378486,
  94.39277672362118,
  6.346958450260396,
  29.56123711211159,
  74.8117128443656,
  115.5617979284524,
  17.555682944633705,
  24.605936950663086,
  141.55727288859697,
  58.54508817669632,
  28.739716613864232,
  26.226336789426146,
  49.147688160152526,
  34.81408404971341,
  74.77143674226606,
  517.048739068739,
  18.226501409792487,
  8.711045803339946,
  46.97639509392978,
  214.27731171191277,
  30.349908951550802,
  29.469963141960495,
  24.802438125808365,
  124.94864884208872,
  52.409038775570245,
  37.219864730453786,
  105.9241515787472,
  82.91094895760472,
  170.55872320965108,
  101.07256238260877,
  185.1427970797945,
  40.43391062291363,
  16.253770377886774,
  10.625985547298702,
  25.522355193560234,
  6.021618595032555,
  49.49888562316846,
  29.957535500909582,
  173.88775941692262,
  203.56731118014963,
  10.951504453771573,
  24.696898632965173,
  13.363405267913505,
  6.855916618696353,
  28.812161492368105,
  31.664977251124427,
  12.33588087671406,
  55.6313817056442,
  45.469855298211236,
  72.02802761282098,
  78.30722344921077,
  9.194128811396613,
  40.65508410729425,
  64.99007470630723,
  65.56742580181779,
  8.973030365894747,
  45.04201811832388,
  55.71597256887867,
  87.86953020720199,
  132.5650753261358,
  19.158569090946646,
  72.5907387068944,
  49.393722498481225,
  18.585572143673335,
  28.667635431854634,
  100.06111359732867,
  30.905561285705293,
  10.37691647226926,
  10.494858268227146,
  13.820311339199439,
  17.4115480964118,
  6.91981521429716,
  21.079867375164344,
  104.97289321266481,
  170.21304565625908,
  7.771967007667543,
  49.66112050609631,
  60.92432429660513,
  68.1286865792631,
  91.45217897216054])

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)