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

([3516563.659866176,
  3522453.125,
  3616325.0,
  3650026.5625,
  3723081.982589908,
  3761592.036803253,
  3849862.938197821,
  3887018.75,
  3895548.4375,
  3895574.7294298923,
  3898183.742516835,
  3900555.2622644934,
  3901127.8228500597,
  3905273.4146939702,
  3913945.3125,
  3914570.7848742297,
  3918353.2176982756,
  3958741.5138268527,
  3970011.4385161516,
  3979799.652759667,
  3985645.7356110835,
  3994223.9795589126,
  4001441.8845774713,
  4033557.597831782,
  4058732.589787411,
  4067925.2378966655,
  4069324.8724835906,
  4086329.609696547,
  4109671.081690553,
  4114940.2682108535,
  4117817.0639380845,
  4125131.833253878,
  4125265.441827439,
  4131237.725549722,
  4131240.625,
  4150953.2070430163,
  4165110.9375,
  4189198.9001408773,
  4189249.166159348,
  4189553.8027579724,
  4189646.844598682,
  4192839.0507922294,
  4196190.625,
  4212550.0,
  4214835.9375,
  4218457.50684519,
  4221790.060107657,
  4223167.049386633,
  4223593.65076936,
  4236418.788865045,
  4241158.793132736,
  4246851.921913718,
  4260638.556822054,
  4262944.518657322,
  4263064.848772324,
  4263129.9123975225,
  4264238.760884773,
  4264580.94333534,
  4269443.575760565,
  4287044.235263441,
  4320328.393560219,
  4320855.548561205,
  4322767.342182601,
  4323065.625,
  4323229.6875,
  4323620.57602832,
  4329303.125,
  4333546.875,
  4336798.1013282575,
  4336802.888449427,
  4340942.448067964,
  4378198.25385809,
  4379257.716371073,
  4439030.943377524,
  4477696.875,
  4478680.133732051,
  4478698.647567996,
  4479185.133383408,
  4484549.066246027,
  4485729.642089781,
  4485742.44386866,
  4534054.199752554,
  4551546.875,
  4559218.112855378,
  4563014.0625,
  4564523.4375,
  4564706.823085847,
  4565094.13595249,
  4566114.0625,
  4566659.11163613,
  4581139.056734097,
  4593009.375,
  4594487.5,
  4633616.59770668,
  4635914.0625,
  4638708.995465376,
  4667414.0625,
  4669028.273503867,
  4670405.79649102,
  4683790.372460848,
  4753434.255329365,
  4754104.859598027,
  4767010.798104672,
  4786629.6875,
  4787223.169701784,
  4787296.81704278,
  4789484.375,
  4805927.3106456995,
  4835882.8125,
  4852521.875,
  4853962.457951159,
  4901526.5625,
  4901549.697415754,
  4910900.0,
  4911341.937839819,
  4930861.602936756,
  4961820.3125,
  4964982.8125,
  4966923.019864093,
  5055995.609911848,
  5056321.406174996,
  5057182.8125,
  5057650.0,
  5058361.203861539,
  5059817.890399953,
  5060273.4375,
  5060786.064762945,
  5072169.356117767,
  5107959.710756048,
  5161585.9375,
  5175995.354599998,
  5177107.376896036,
  5178205.430695763,
  5178362.5,
  5180024.552726621,
  5190468.888403881,
  5192246.14109956,
  5239840.625,
  5248677.229435093,
  5248686.407838635,
  5248691.690909987,
  5287351.5625,
  5294191.06705885,
  5294198.948810237,
  5294745.3125,
  5328335.50160979,
  5350006.521177989,
  5361277.704009917,
  5361306.943817168,
  5431927.875656571,
  5450353.56765452,
  5450434.4694744935,
  5465775.885359852,
  5477467.892138546,
  5477475.232293967,
  5480420.3125,
  5507400.907576567,
  5539351.247810291,
  5540832.8125,
  5569589.556552177,
  5637797.152746597,
  5674886.132797357,
  5676977.644666485,
  5709127.631710819,
  5839211.621115992,
  6506453.106326316,
  6514916.316738477,
  6527396.51166487,
  6528086.061037031,
  6528271.067004249,
  6528316.379315424,
  6528324.667390298,
  6528546.182630836,
  6528578.102955911,
  6528600.149832682,
  6671662.5,
  6805143.75,
  6805822.120942643,
  6819506.25,
  6820201.955834646,
  6820232.625928741,
  6822763.93460819,
  6908504.296530955,
  6928553.125,
  6929986.028857499,
  6932011.688049375,
  7070934.375,
  7070946.549266327,
  7184259.200758169,
  7185426.502731286,
  7185976.5625,
  7193747.411353213,
  7328268.812119185,
  7328951.200370186,
  7331853.125,
  7334797.1820328785,
  7337451.5625,
  7346221.875,
  7376234.375,
  7376248.497398315,
  7385063.05969682,
  7388484.587637318,
  7430701.5625,
  7457791.727214021,
  7457888.556073201,
  7459475.0,
  7462050.0,
  7462079.330295611,
  7466277.458640598,
  7479812.402334075,
  7480492.1875,
  7484383.5667245835,
  7487419.249669934,
  7489599.98710439,
  7492492.093171464,
  7493315.7290560575,
  7499106.96147352,
  7500438.968430275,
  7500525.883778526,
  7500663.9830243075,
  7500664.092543603,
  7501518.75,
  7501529.964638061,
  7502032.345091607,
  7502081.157814215,
  7502651.877024314,
  7502660.37610205,
  7502665.625,
  7502792.1875,
  7503533.725369856,
  7503871.007596663,
  7504671.739338248,
  7504673.096900943,
  7504750.966245916,
  7504963.172085897,
  7505013.184626922,
  7505110.59045512,
  7505275.632001752,
  7505277.89805831,
  7505306.25,
  7505354.6875,
  7505392.1875,
  7505503.125,
  7505539.081645869,
  7505595.49870544,
  7505752.272668069,
  7505762.358975728,
  7505903.1013427805,
  7505973.4375,
  7506116.900846033,
  7506243.558663944,
  7506325.0,
  7506437.5,
  7506442.887391759,
  7506475.840766794,
  7506678.713185535,
  7506785.814736297,
  7506843.722164545,
  7507137.5,
  7507309.005697363,
  7507353.125,
  7507767.310103823,
  7508012.48385256,
  7508507.640571699,
  7508546.430646867,
  7508647.820744732,
  7509394.974935964,
  7510080.158732645,
  7510089.650962341,
  7510096.845811033,
  7510366.4454933675,
  7511011.988574928,
  7511096.875,
  7511950.947257247,
  7513034.781752089,
  7513100.0,
  7513442.1875,
  7513553.369248716,
  7513650.0,
  7514315.16715968,
  7515747.343673917,
  7516033.86141414,
  7516110.248588808,
  7516182.732826867,
  7525296.51479232,
  7526988.996935411,
  7527384.375,
  7527796.875,
  7558241.802918784,
  7558776.248698344,
  7561118.276115736,
  7567262.806755644,
  7576693.75,
  7578590.5847103195,
  7621339.41884141,
  7621541.846642729,
  7621578.827208772,
  7622125.0,
  7622158.158586307,
  7623459.388809071,
  7623473.4375,
  7623816.769938677,
  7623865.28598947,
  7625083.38450811,
  7626881.25,
  7627787.5,
  7629726.325955533,
  7669481.028723511,
  7669493.083451879,
  7714184.61886304,
  7715714.500899253,
  7721634.375,
  7725970.871652735,
  7725998.7184672905,
  7727196.875,
  7728661.949523789,
  7731865.492599414,
  7739229.817074609,
  7754009.056972649,
  7758397.618610811,
  7758751.5625,
  7776972.935246239,
  7801048.249968431,
  7830308.459139842,
  7843560.9375,
  7845811.575901804,
  7850957.8125,
  7850977.587026083,
  7851108.426781962,
  7851614.0625,
  7857219.566250149,
  7858342.014396781,
  7861776.5625,
  7871777.579991633,
  7872139.766202867,
  7874850.0,
  7878893.75,
  7878911.689775805,
  7880784.344459365,
  7880838.023097858,
  7885476.5625,
  7885601.5625,
  7887880.148248193,
  7889028.1478606,
  7891131.712528821,
  7900521.939498802,
  7907994.174336239,
  7908909.375,
  7909715.581591999,
  7914585.9375,
  7917464.424186572,
  7919542.520877404,
  7919951.718326982,
  7920125.0,
  7920556.129278812,
  7920565.249017154,
  7930822.401616633,
  7931471.218914138,
  7935716.83438859,
  7935820.673969362,
  7936287.5,
  7937709.375,
  7938704.6875,
  7939506.675765278,
  7941996.636586881,
  7942496.597864096,
  7942515.4397056755,
  7945003.847510337,
  7950767.163955314,
  7951167.1875,
  7953253.297134905,
  7954354.248276682,
  7954564.61794186,
  7954773.3243605485,
  7954909.783049829,
  7954950.545597302,
  7955018.203875645,
  7955091.117820228,
  7955109.375,
  7955382.8125,
  7955549.207529363,
  7956063.255779153,
  7956654.938679654,
  7956731.9267000845,
  7956804.3683185885,
  7956948.598171674,
  7957244.0949284285,
  7957581.588556634,
  7957622.185436127,
  7957793.75,
  7958537.107324194,
  7958808.211149185,
  7959150.502406015,
  7959260.569680303,
  7959454.274185324,
  7959469.582145377,
  7959478.40863012,
  7960003.122486739,
  7960014.0625,
  7960076.5625,
  7961015.625,
  7961285.766834207,
  7963045.3125,
  7963910.9375,
  7968078.631979903,
  7968616.648323349,
  7978218.627613123,
  7978846.875,
  7979496.875,
  7980857.8125,
  7981467.502001785,
  7982330.224749004,
  7983459.304819967,
  7984481.25,
  7984924.764752028,
  7985193.696799739,
  7985203.125,
  7985332.8125,
  7985451.999197627,
  7985938.309492799,
  7986282.8125,
  7986432.508787645,
  7986596.0434040325,
  7986713.6566671375,
  7986745.3125,
  7987198.287837345,
  7987228.125,
  7987928.768118402,
  7988132.289831875,
  7988139.0625,
  7988231.25,
  7988332.8125,
  7988335.809714023,
  7988348.230784366,
  7988459.375,
  7988854.6875,
  7988875.558616115,
  7988876.314511786,
  7990432.574515623,
  7991405.464900492,
  7991929.237181021,
  7992604.6875,
  7992660.604839769,
  7992774.346633762,
  7992804.6875,
  7992811.072266713,
  7992923.153195937,
  7992981.090894862,
  7993001.41287108,
  7993037.117982465,
  7993169.199880819,
  7993277.526542489,
  7993702.261217761,
  7993907.4795841705,
  7994029.21684611,
  7994227.833224316,
  7994539.543126533,
  7995877.317126113,
  7995916.748147555,
  8000692.1875,
  8011682.8125,
  8014434.375,
  8015899.054810305,
  8025685.9375,
  8032218.1172520835,
  8047057.080455237,
  8047248.653650699,
  8050101.520064619,
  8050116.534531348,
  8053466.341432249,
  8054517.039802403,
  8059988.231678637,
  8060043.700539757,
  8061199.668521936,
  8061385.217663492,
  8061429.6875,
  8063539.024374357,
  8064121.888140827,
  8065523.4375,
  8067077.707971298,
  8067296.653043272,
  8068217.1875,
  8075822.911583884,
  8075920.3125,
  8077877.847967238,
  8082654.6875,
  8085544.452148123,
  8113979.242727022,
  8114130.202462289,
  8114766.771018609,
  8117300.3447585935,
  8117887.355749825,
  8118154.107213509,
  8118401.552081097,
  8119190.625,
  8119407.8125,
  8119784.286486917,
  8120111.389292315,
  8120337.5,
  8121123.741568458,
  8121775.46166475,
  8122257.300305839,
  8122411.211256543,
  8123620.902497167,
  8124050.235486784,
  8124765.245954152,
  8125026.876584971,
  8125293.129518989,
  8126296.875,
  8126345.3125,
  8126849.679525267,
  8126854.570839319,
  8127434.257318617,
  8127610.301968664,
  8128060.806240722,
  8128788.326300505,
  8130394.441011572,
  8148213.394377924,
  8154723.1846909765,
  8164401.5625,
  8166102.663778551,
  8192096.875,
  8201456.25,
  8201479.758199686,
  8201501.5625,
  8210608.064887031,
  8220403.651414304,
  8237193.75,
  8279245.3125,
  8279327.193462723,
  8281749.511716898,
  8282440.6902050655,
  8282845.3125,
  8282847.637554568,
  8282851.5625,
  8282851.682549346,
  8283321.875,
  8283490.302014268,
  8293037.5,
  8431257.6909584,
  8444651.574392842,
  8518851.419710709,
  8543912.817068344,
  8659973.68822171,
  8671880.570053363,
  10725091.041000225,
  10814141.783209614,
  10814938.866117498,
  10816056.25,
  10816510.66763128,
  10817462.353888037,
  10822502.389516778,
  10823751.5625,
  10844075.80897423,
  10845218.187316788,
  10845837.407292234,
  10846184.345855124,
  10847018.75,
  11805764.029285511,
  11921743.75,
  11949632.435519068,
  11950119.943488978,
  11995270.461129945,
  12021439.0625,
  12021790.213106142,
  12021903.76831068,
  12023928.058359113,
  12029901.5625,
  12032356.25,
  12032566.330043258,
  12033127.680007346,
  12033213.651444098,
  12033272.094065566,
  12033376.625534093,
  12033396.719807995,
  12033728.125,
  12033768.75,
  12033795.942995561,
  12033824.607349725,
  12034171.185366008,
  12034991.840709705,
  12035120.3125,
  12035567.1875,
  12036231.607004987,
  12037243.75,
  12045588.584479712,
  12048837.5,
  12050963.139607713,
  12058795.3125,
  12762587.5,
  12956062.554147929,
  13053614.281567069,
  13068352.457186054,
  13075096.875,
  13075639.945638333,
  13077987.5,
  13086819.859345548,
  13086821.652011884,
  13089836.260993473,
  13091025.479852986,
  13103281.93071233,
  13138178.125,
  13141278.125,
  13170725.91664031,
  13595979.993300736,
  13882415.300789392,
  16412244.98934245,
  16414216.73666595,
  16415442.1875,
  16416021.875,
  16416364.017114962,
  16416374.618944822,
  16416996.875,
  16417047.034869852,
  16417461.81319017,
  16420471.109906167,
  16420676.16772352,
  16421679.301742032,
  16421807.8125,
  16423380.131742321,
  16423774.339786593,
  16424723.008657621,
  16428105.59359269,
  16428929.6875,
  16432030.922882574,
  16437589.0625,
  16438021.201651422,
  16438662.37073636,
  16438719.414961193,
  16442312.300715653,
  16442533.38857625,
  16443037.5,
  16446536.208341783,
  16446818.70456963,
  16447123.4375,
  16448395.977873702,
  16454951.51280248],
 [7.01948642501036,
  66.5344752909091,
  60.31015088691055,
  58.807249870126746,
  44.03904779972045,
  59.37580299397437,
  11.889464342802523,
  37.13071448641459,
  45.1374148155577,
  5.599248888212579,
  12.977282600400674,
  58.27258785089079,
  34.82352582445307,
  21.62640621979545,
  46.920457061113126,
  27.71995546625717,
  27.423223830130336,
  19.919517997404416,
  94.05985280775124,
  6.434001922751892,
  7.763790653072909,
  190.74380114684226,
  46.35462406180823,
  24.905573582526053,
  51.52242626319719,
  11.521093839300718,
  8.220276325770298,
  27.074976173999428,
  5.7575248619425,
  13.736160203447916,
  7.0893454247500065,
  12.180726844509678,
  9.181178514383564,
  11.499558965669578,
  93.11351904384259,
  15.263719875607402,
  58.31284022111421,
  11.937019493941452,
  5.082390729549531,
  45.39026707473275,
  59.635129428750794,
  46.19213641925138,
  63.528374628132724,
  72.74818997139262,
  51.54304523246065,
  10.681349172057688,
  8.895443407269585,
  7.0856130415657494,
  26.70561811891589,
  15.091743275262589,
  93.65821982617186,
  5.354390188730218,
  92.70727539858609,
  40.220045566177134,
  7.9147349230221895,
  29.786147033160134,
  23.73013127136719,
  45.14792189579634,
  54.48486962159565,
  17.40342296553851,
  6.116479587207561,
  10.183711728245376,
  15.565793254460981,
  52.715997787074635,
  37.79805709352883,
  5.723737457916772,
  35.363379082082126,
  67.91281922761925,
  22.041392775048468,
  24.596528704896908,
  11.560974923742211,
  14.28489046241263,
  19.828347753205705,
  5.913140505281145,
  30.84927515445939,
  24.26478682847881,
  15.80757787538496,
  30.975052760046104,
  7.847658584364835,
  27.65711666150794,
  88.94033302056516,
  13.35972862675263,
  29.254265501214082,
  18.806178073471457,
  103.21065893487528,
  88.21821869322318,
  27.494018652988636,
  69.08801818417956,
  98.64406292105221,
  5.450575079151153,
  6.762355996254278,
  59.34291309079271,
  66.99383363677563,
  5.505124699710122,
  45.204167937297555,
  38.001460838548276,
  67.32652069615425,
  15.539450498463784,
  47.39848625237699,
  20.75853977662154,
  26.36428938027978,
  6.514411803077918,
  77.26025906445722,
  31.153081766231978,
  14.010404894940924,
  81.78434290161276,
  90.36293292109438,
  17.587385091516317,
  90.29259681613227,
  38.67964632334759,
  81.81294045640266,
  78.65935000452662,
  111.1872621520416,
  74.43634351252643,
  14.06827609612685,
  42.99775472276282,
  47.129850063882856,
  32.37198239892718,
  8.077418413775064,
  77.92308264200847,
  35.235842860851704,
  74.71835895380818,
  35.10502150440469,
  66.98134535188979,
  12.438428781658171,
  57.49066791707384,
  22.80013103028026,
  17.52453073898031,
  29.910045229436932,
  119.65984004010689,
  5.166343226460405,
  22.16047291462894,
  18.806627203819815,
  56.37830779705009,
  126.07292485672849,
  33.40471336476145,
  93.43805855070607,
  151.4559803758041,
  29.199491004698487,
  11.952120088919054,
  50.67350034480911,
  59.15314927195181,
  16.944692324623,
  12.03351395896264,
  32.560997706221364,
  22.157391524439905,
  9.571731739585307,
  22.319217459656556,
  23.644373979877276,
  7.515477096229745,
  39.81666046625503,
  6.097979143852031,
  8.81772402111189,
  10.007909185576217,
  6.524132495346957,
  136.0984072982323,
  20.302643101868046,
  97.37256540234979,
  39.9675556136904,
  5.738939188561282,
  9.577840303031515,
  5.470211823725413,
  16.271294710365076,
  21.560306464015177,
  27.033278559497933,
  21.513171879898294,
  62.996352305148434,
  77.2052304593661,
  11.326454062750171,
  118.43508598433738,
  7.3349376306266425,
  34.32894286558482,
  16.587694927966236,
  35.48357472115026,
  6.103820333534942,
  140.50183262448172,
  113.1949923592802,
  15.969823229073276,
  95.83185575224458,
  108.79473019544928,
  6.243910737215056,
  149.509639765815,
  95.81980058999328,
  71.34207437810113,
  56.45603567386273,
  39.82854734342361,
  55.34647400741533,
  8.121521393567496,
  111.53684735254822,
  8.255221518407486,
  51.023347804451,
  52.393514509451975,
  15.183685461573914,
  13.720195522184644,
  120.34628517000759,
  86.23788257531342,
  62.6340797745629,
  50.57864844152558,
  88.37294311637322,
  7.625046506672427,
  7.819042973771666,
  11.477175512841768,
  39.696999397822424,
  21.898632992768597,
  10.41876176278347,
  55.897453797976034,
  106.6211841686959,
  28.30799429673197,
  19.979940391073555,
  10.978404753405124,
  72.33804379425865,
  19.440692424021194,
  17.2629011687243,
  14.892584152340246,
  67.9990621148876,
  25.93454057980983,
  94.91739921347983,
  6.744870477656062,
  47.88189341391646,
  74.44776034856416,
  8.302529041937728,
  40.73971303264943,
  224.0079844782543,
  47.78312194934922,
  14.601408600712041,
  30.581149759162695,
  30.389216439980547,
  48.2979307383526,
  50.487077091111054,
  51.16086882617536,
  7.875602407220497,
  31.201754631804178,
  13.784890367903143,
  10.308187179208561,
  18.319028179826905,
  36.40186767606083,
  23.167850728245455,
  18.05838450240185,
  65.59209701472635,
  30.049248143749118,
  34.49003119887807,
  94.78521351051369,
  58.596978421854814,
  15.014990136100439,
  11.431789746274037,
  17.413647720303224,
  209.50872831091377,
  28.06506262159766,
  94.40611345419259,
  48.228949350596004,
  14.25914984041219,
  88.15775446547701,
  44.489907098302666,
  36.52467891332402,
  6.379607413714642,
  15.975034794336581,
  57.3329413597779,
  11.084475160518565,
  49.578195102523566,
  22.904101107776782,
  40.78986831271144,
  29.95084325543879,
  115.0640568074223,
  65.47489978498575,
  22.066768495516133,
  20.402444136158262,
  10.638110414522798,
  17.103848046619827,
  6.88696432167989,
  14.831409765521775,
  19.37674956878004,
  61.380033779791745,
  53.60163648888627,
  20.41568583529521,
  65.5778558280526,
  49.78491419066978,
  76.25758538163109,
  168.3829980692553,
  52.1180371437178,
  28.11794278661078,
  17.089728502558078,
  21.358636212983704,
  155.1979471989042,
  27.04226193325171,
  19.042624472018247,
  27.29147294117218,
  49.381286963217555,
  160.245494998116,
  30.895713893257824,
  10.670253439284917,
  20.88749173370106,
  37.4522393876569,
  44.70160904055175,
  27.751877502028695,
  51.684699814098614,
  13.784109211911154,
  136.3354759770187,
  106.00913812970329,
  19.439721314241293,
  9.34310688158804,
  53.073200540790225,
  10.532256861187141,
  23.264906042686434,
  17.685943543814037,
  64.92706328246186,
  87.68537595769766,
  50.49535000610757,
  24.811255717056913,
  19.978375402761404,
  15.960292882849512,
  17.98589602723802,
  83.82272723815153,
  6.148242591149478,
  54.74640876897011,
  30.174629157741613,
  6.983895201946333,
  26.156925450808387,
  15.4595102923306,
  40.72108683395517,
  12.530019091545823,
  93.22085259295635,
  90.62572792850523,
  9.177713572442869,
  99.69467097176545,
  91.6322453520697,
  21.30538608761715,
  65.57212150464937,
  24.936651580767368,
  17.94623503277572,
  62.033853505650875,
  5.224728526227363,
  83.36973492141257,
  76.3338083851332,
  36.668112067345135,
  17.98113415751144,
  70.53967302145612,
  71.45082073703959,
  14.114306917976986,
  65.76645680316268,
  12.057254425441004,
  68.03078799304134,
  28.619368815857676,
  32.2342117444501,
  86.62283885545854,
  17.05594592248561,
  32.88571888912515,
  9.717971190250434,
  39.884790934033774,
  87.04908712868436,
  58.060989958416485,
  63.14092819914653,
  34.27977224844754,
  5.716747544554541,
  65.56672993227527,
  38.32799547832274,
  22.826377473400147,
  17.41731766624705,
  82.12040332977067,
  48.07284305482992,
  18.433259980650668,
  133.33811851050376,
  103.93038694304727,
  46.30537083409534,
  35.73305287643594,
  66.25242789247426,
  79.29784991146443,
  10.849926253103032,
  16.90234586875447,
  60.633656655178385,
  28.752072152929404,
  24.498565015031357,
  103.67552141114194,
  30.150690250141423,
  21.750404567447887,
  6.847799763690607,
  17.528318592644773,
  31.034685844103823,
  11.996567471438304,
  45.498316296673,
  153.25420061452704,
  144.8597492405272,
  10.064211187533502,
  15.993892629083188,
  9.977616029560643,
  160.6320917151447,
  15.591445904642253,
  141.09952847124708,
  28.308173776293337,
  16.27334196853476,
  76.03808418491538,
  9.020068670224314,
  16.715084763437552,
  28.294177136650532,
  168.04138273387548,
  27.600518401003228,
  19.466610822228752,
  26.042167706909794,
  28.475501963253397,
  48.67486279225878,
  40.5105835016637,
  35.44033212604757,
  47.25530165951515,
  32.358615543203015,
  43.616184286091375,
  17.304429020167554,
  11.43838791393536,
  8.222900386474333,
  99.22619737560433,
  31.311807509626046,
  63.68146441032566,
  41.063377500850876,
  17.47926567105375,
  9.272758695425004,
  82.96550919229492,
  80.01576952123023,
  27.514415917045337,
  84.6703510835849,
  41.24248792088897,
  9.730131270595315,
  8.766677985748457,
  111.40842285845807,
  63.75469372622584,
  25.154042451332636,
  13.591565044220127,
  39.310501605966834,
  14.400183113800788,
  60.33421256013218,
  342.2368594861819,
  21.269915192104396,
  28.714816815368398,
  48.673591491272134,
  29.500433849373213,
  68.43004407218872,
  7.068152323284567,
  38.969453627299565,
  84.64689377389168,
  12.12968255252203,
  68.42734821563062,
  24.576602932665264,
  7.910889176187939,
  10.489930692405192,
  75.07853181726583,
  6.698862382409513,
  131.60487806034422,
  59.605057640178224,
  23.387563555654168,
  5.688899929952562,
  19.841634402843326,
  38.624652918136434,
  73.52091928047172,
  49.79276791155743,
  20.564115409617752,
  18.49774427546654,
  157.19093796273404,
  17.415993167074536,
  32.07119962571495,
  8.558218589611204,
  5.756528691766511,
  33.509400962997184,
  55.752147825836175,
  36.84681638070083,
  135.6936258451777,
  45.83424469415289,
  165.33858048887268,
  29.033364060284192,
  30.543214601143298,
  11.483252039375941,
  57.70960286832259,
  151.1827831237001,
  160.21111183262104,
  147.96336756459135,
  12.50296055739069,
  27.582708274366365,
  13.812566580984075,
  19.548654512312712,
  115.13868251898653,
  14.798156623344846,
  47.02515780789199,
  59.87398114485782,
  20.176143891995384,
  5.920306493509064,
  41.67058821352796,
  71.9424026223688,
  50.502555018418754,
  89.20243131304137,
  59.73242057507242,
  41.268084508405074,
  13.479888285782101,
  9.8084152873118,
  10.533662159678805,
  33.810940246315376,
  173.6890138636346,
  19.653637438452986,
  14.901945593704594,
  178.4985370515093,
  31.006518402273013,
  185.0410946524026,
  23.166356589444973,
  128.94315369225816,
  164.98825986850855,
  488.81047672436205,
  13.292772709478893,
  5.834517441355705,
  12.194814287675038,
  12.848638265611687,
  33.59336376401455,
  33.18862946613102,
  27.211577708324278,
  33.9285353726062,
  72.89399692287354,
  9.07436784130068,
  5.287497241740275,
  9.233840478109595,
  56.96200469664292,
  8.218363860710046,
  10.143708939945652,
  71.31742272978735,
  91.26734705277974,
  24.372857176076362,
  39.7621301569106,
  10.475971581836353,
  83.1366952685337,
  30.678923665187845,
  15.218009317075259,
  59.70791589325795,
  88.37888095748049,
  31.270848558495462,
  44.27785002988454,
  106.35383576934544,
  54.70533148395681,
  82.49152650932639,
  46.46836632088982,
  49.95622134319276,
  26.33377504486552,
  31.13012169604691,
  7.647242348643645,
  87.70916733128635,
  244.3491998399123,
  44.18998447137148,
  27.031806457671486,
  11.14634710319366,
  71.22879186649756,
  11.662318379562311,
  15.991783934276732,
  19.868060324019236,
  137.22063823086242,
  42.38645826180944,
  14.242339261180264,
  39.43918152336181,
  8.148842273041652,
  10.90675470673964,
  21.337718957501526,
  72.14724480940575,
  105.48717427939776,
  20.88979017664333,
  26.816781512226054,
  35.31175626525021,
  72.18762454173917,
  32.77324708509002,
  43.426352111209674,
  6.220716156886273,
  60.09042073058847,
  8.366770954639204,
  29.119503034392448,
  22.419461717742248,
  70.90023495316017,
  12.543898943657375,
  51.84963440419426,
  36.73649390827758,
  30.99849900974961,
  5.442985697440052,
  14.24141571826322,
  15.847843284458689,
  8.320194482960225,
  56.50965716339728,
  37.40800654950942,
  38.56510713587179,
  53.134957484725234,
  34.35026608747548,
  10.197422405449903,
  23.170910854904044,
  63.732374890056654,
  64.43720126789721,
  24.477392453311655,
  65.25096068944183,
  74.38590144420631,
  108.6739682549958,
  18.107091980841762,
  102.8265977671968,
  58.39070447933071,
  42.01451500795805,
  9.523770418462753,
  6.95631268600118,
  31.62766201396243,
  25.244079148727934,
  74.71889691396936,
  31.236141190981193,
  24.786967037988923,
  6.122398380105475,
  5.417140765627537,
  14.461384871439568,
  78.53648675725556,
  42.77447559141699,
  28.86931685524332,
  60.99998310461749,
  6.227087633799472,
  107.47634152240818,
  14.095113623373695,
  88.59071709566189,
  41.840726521834654,
  27.657580766254853,
  22.767794584542173,
  73.48899211355393,
  12.114093641176279,
  33.7663723096822,
  26.54422324200455,
  22.61116499362535,
  22.713810400142393,
  139.77593730435552,
  9.737580140139878,
  20.031078937212314,
  126.84831358045426,
  20.06291603638016,
  39.15799931658964,
  9.069319639681057,
  48.26116806569466,
  90.16281308540533,
  46.14274731649472,
  33.22964088097305,
  14.23269848518827,
  191.169578272674,
  108.14504140695574,
  7.704224694560673,
  14.771132032787731,
  237.37043882870378,
  5.75045442817992,
  11.050904857933716])

Energy spectrum

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

([3516563.659866176,
  3522453.125,
  3616325.0,
  3650026.5625,
  3723081.982589908,
  3761592.036803253,
  3849862.938197821,
  3887018.75,
  3895548.4375,
  3895574.7294298923,
  3898183.742516835,
  3900555.2622644934,
  3901127.8228500597,
  3905273.4146939702,
  3913945.3125,
  3914570.7848742297,
  3918353.2176982756,
  3958741.5138268527,
  3970011.4385161516,
  3979799.652759667,
  3985645.7356110835,
  3994223.9795589126,
  4001441.8845774713,
  4033557.597831782,
  4058732.589787411,
  4067925.2378966655,
  4069324.8724835906,
  4086329.609696547,
  4109671.081690553,
  4114940.2682108535,
  4117817.0639380845,
  4125131.833253878,
  4125265.441827439,
  4131237.725549722,
  4131240.625,
  4150953.2070430163,
  4165110.9375,
  4189198.9001408773,
  4189249.166159348,
  4189553.8027579724,
  4189646.844598682,
  4192839.0507922294,
  4196190.625,
  4212550.0,
  4214835.9375,
  4218457.50684519,
  4221790.060107657,
  4223167.049386633,
  4223593.65076936,
  4236418.788865045,
  4241158.793132736,
  4246851.921913718,
  4260638.556822054,
  4262944.518657322,
  4263064.848772324,
  4263129.9123975225,
  4264238.760884773,
  4264580.94333534,
  4269443.575760565,
  4287044.235263441,
  4320328.393560219,
  4320855.548561205,
  4322767.342182601,
  4323065.625,
  4323229.6875,
  4323620.57602832,
  4329303.125,
  4333546.875,
  4336798.1013282575,
  4336802.888449427,
  4340942.448067964,
  4378198.25385809,
  4379257.716371073,
  4439030.943377524,
  4477696.875,
  4478680.133732051,
  4478698.647567996,
  4479185.133383408,
  4484549.066246027,
  4485729.642089781,
  4485742.44386866,
  4534054.199752554,
  4551546.875,
  4559218.112855378,
  4563014.0625,
  4564523.4375,
  4564706.823085847,
  4565094.13595249,
  4566114.0625,
  4566659.11163613,
  4581139.056734097,
  4593009.375,
  4594487.5,
  4633616.59770668,
  4635914.0625,
  4638708.995465376,
  4667414.0625,
  4669028.273503867,
  4670405.79649102,
  4683790.372460848,
  4753434.255329365,
  4754104.859598027,
  4767010.798104672,
  4786629.6875,
  4787223.169701784,
  4787296.81704278,
  4789484.375,
  4805927.3106456995,
  4835882.8125,
  4852521.875,
  4853962.457951159,
  4901526.5625,
  4901549.697415754,
  4910900.0,
  4911341.937839819,
  4930861.602936756,
  4961820.3125,
  4964982.8125,
  4966923.019864093,
  5055995.609911848,
  5056321.406174996,
  5057182.8125,
  5057650.0,
  5058361.203861539,
  5059817.890399953,
  5060273.4375,
  5060786.064762945,
  5072169.356117767,
  5107959.710756048,
  5161585.9375,
  5175995.354599998,
  5177107.376896036,
  5178205.430695763,
  5178362.5,
  5180024.552726621,
  5190468.888403881,
  5192246.14109956,
  5239840.625,
  5248677.229435093,
  5248686.407838635,
  5248691.690909987,
  5287351.5625,
  5294191.06705885,
  5294198.948810237,
  5294745.3125,
  5328335.50160979,
  5350006.521177989,
  5361277.704009917,
  5361306.943817168,
  5431927.875656571,
  5450353.56765452,
  5450434.4694744935,
  5465775.885359852,
  5477467.892138546,
  5477475.232293967,
  5480420.3125,
  5507400.907576567,
  5539351.247810291,
  5540832.8125,
  5569589.556552177,
  5637797.152746597,
  5674886.132797357,
  5676977.644666485,
  5709127.631710819,
  5839211.621115992,
  6506453.106326316,
  6514916.316738477,
  6527396.51166487,
  6528086.061037031,
  6528271.067004249,
  6528316.379315424,
  6528324.667390298,
  6528546.182630836,
  6528578.102955911,
  6528600.149832682,
  6671662.5,
  6805143.75,
  6805822.120942643,
  6819506.25,
  6820201.955834646,
  6820232.625928741,
  6822763.93460819,
  6908504.296530955,
  6928553.125,
  6929986.028857499,
  6932011.688049375,
  7070934.375,
  7070946.549266327,
  7184259.200758169,
  7185426.502731286,
  7185976.5625,
  7193747.411353213,
  7328268.812119185,
  7328951.200370186,
  7331853.125,
  7334797.1820328785,
  7337451.5625,
  7346221.875,
  7376234.375,
  7376248.497398315,
  7385063.05969682,
  7388484.587637318,
  7430701.5625,
  7457791.727214021,
  7457888.556073201,
  7459475.0,
  7462050.0,
  7462079.330295611,
  7466277.458640598,
  7479812.402334075,
  7480492.1875,
  7484383.5667245835,
  7487419.249669934,
  7489599.98710439,
  7492492.093171464,
  7493315.7290560575,
  7499106.96147352,
  7500438.968430275,
  7500525.883778526,
  7500663.9830243075,
  7500664.092543603,
  7501518.75,
  7501529.964638061,
  7502032.345091607,
  7502081.157814215,
  7502651.877024314,
  7502660.37610205,
  7502665.625,
  7502792.1875,
  7503533.725369856,
  7503871.007596663,
  7504671.739338248,
  7504673.096900943,
  7504750.966245916,
  7504963.172085897,
  7505013.184626922,
  7505110.59045512,
  7505275.632001752,
  7505277.89805831,
  7505306.25,
  7505354.6875,
  7505392.1875,
  7505503.125,
  7505539.081645869,
  7505595.49870544,
  7505752.272668069,
  7505762.358975728,
  7505903.1013427805,
  7505973.4375,
  7506116.900846033,
  7506243.558663944,
  7506325.0,
  7506437.5,
  7506442.887391759,
  7506475.840766794,
  7506678.713185535,
  7506785.814736297,
  7506843.722164545,
  7507137.5,
  7507309.005697363,
  7507353.125,
  7507767.310103823,
  7508012.48385256,
  7508507.640571699,
  7508546.430646867,
  7508647.820744732,
  7509394.974935964,
  7510080.158732645,
  7510089.650962341,
  7510096.845811033,
  7510366.4454933675,
  7511011.988574928,
  7511096.875,
  7511950.947257247,
  7513034.781752089,
  7513100.0,
  7513442.1875,
  7513553.369248716,
  7513650.0,
  7514315.16715968,
  7515747.343673917,
  7516033.86141414,
  7516110.248588808,
  7516182.732826867,
  7525296.51479232,
  7526988.996935411,
  7527384.375,
  7527796.875,
  7558241.802918784,
  7558776.248698344,
  7561118.276115736,
  7567262.806755644,
  7576693.75,
  7578590.5847103195,
  7621339.41884141,
  7621541.846642729,
  7621578.827208772,
  7622125.0,
  7622158.158586307,
  7623459.388809071,
  7623473.4375,
  7623816.769938677,
  7623865.28598947,
  7625083.38450811,
  7626881.25,
  7627787.5,
  7629726.325955533,
  7669481.028723511,
  7669493.083451879,
  7714184.61886304,
  7715714.500899253,
  7721634.375,
  7725970.871652735,
  7725998.7184672905,
  7727196.875,
  7728661.949523789,
  7731865.492599414,
  7739229.817074609,
  7754009.056972649,
  7758397.618610811,
  7758751.5625,
  7776972.935246239,
  7801048.249968431,
  7830308.459139842,
  7843560.9375,
  7845811.575901804,
  7850957.8125,
  7850977.587026083,
  7851108.426781962,
  7851614.0625,
  7857219.566250149,
  7858342.014396781,
  7861776.5625,
  7871777.579991633,
  7872139.766202867,
  7874850.0,
  7878893.75,
  7878911.689775805,
  7880784.344459365,
  7880838.023097858,
  7885476.5625,
  7885601.5625,
  7887880.148248193,
  7889028.1478606,
  7891131.712528821,
  7900521.939498802,
  7907994.174336239,
  7908909.375,
  7909715.581591999,
  7914585.9375,
  7917464.424186572,
  7919542.520877404,
  7919951.718326982,
  7920125.0,
  7920556.129278812,
  7920565.249017154,
  7930822.401616633,
  7931471.218914138,
  7935716.83438859,
  7935820.673969362,
  7936287.5,
  7937709.375,
  7938704.6875,
  7939506.675765278,
  7941996.636586881,
  7942496.597864096,
  7942515.4397056755,
  7945003.847510337,
  7950767.163955314,
  7951167.1875,
  7953253.297134905,
  7954354.248276682,
  7954564.61794186,
  7954773.3243605485,
  7954909.783049829,
  7954950.545597302,
  7955018.203875645,
  7955091.117820228,
  7955109.375,
  7955382.8125,
  7955549.207529363,
  7956063.255779153,
  7956654.938679654,
  7956731.9267000845,
  7956804.3683185885,
  7956948.598171674,
  7957244.0949284285,
  7957581.588556634,
  7957622.185436127,
  7957793.75,
  7958537.107324194,
  7958808.211149185,
  7959150.502406015,
  7959260.569680303,
  7959454.274185324,
  7959469.582145377,
  7959478.40863012,
  7960003.122486739,
  7960014.0625,
  7960076.5625,
  7961015.625,
  7961285.766834207,
  7963045.3125,
  7963910.9375,
  7968078.631979903,
  7968616.648323349,
  7978218.627613123,
  7978846.875,
  7979496.875,
  7980857.8125,
  7981467.502001785,
  7982330.224749004,
  7983459.304819967,
  7984481.25,
  7984924.764752028,
  7985193.696799739,
  7985203.125,
  7985332.8125,
  7985451.999197627,
  7985938.309492799,
  7986282.8125,
  7986432.508787645,
  7986596.0434040325,
  7986713.6566671375,
  7986745.3125,
  7987198.287837345,
  7987228.125,
  7987928.768118402,
  7988132.289831875,
  7988139.0625,
  7988231.25,
  7988332.8125,
  7988335.809714023,
  7988348.230784366,
  7988459.375,
  7988854.6875,
  7988875.558616115,
  7988876.314511786,
  7990432.574515623,
  7991405.464900492,
  7991929.237181021,
  7992604.6875,
  7992660.604839769,
  7992774.346633762,
  7992804.6875,
  7992811.072266713,
  7992923.153195937,
  7992981.090894862,
  7993001.41287108,
  7993037.117982465,
  7993169.199880819,
  7993277.526542489,
  7993702.261217761,
  7993907.4795841705,
  7994029.21684611,
  7994227.833224316,
  7994539.543126533,
  7995877.317126113,
  7995916.748147555,
  8000692.1875,
  8011682.8125,
  8014434.375,
  8015899.054810305,
  8025685.9375,
  8032218.1172520835,
  8047057.080455237,
  8047248.653650699,
  8050101.520064619,
  8050116.534531348,
  8053466.341432249,
  8054517.039802403,
  8059988.231678637,
  8060043.700539757,
  8061199.668521936,
  8061385.217663492,
  8061429.6875,
  8063539.024374357,
  8064121.888140827,
  8065523.4375,
  8067077.707971298,
  8067296.653043272,
  8068217.1875,
  8075822.911583884,
  8075920.3125,
  8077877.847967238,
  8082654.6875,
  8085544.452148123,
  8113979.242727022,
  8114130.202462289,
  8114766.771018609,
  8117300.3447585935,
  8117887.355749825,
  8118154.107213509,
  8118401.552081097,
  8119190.625,
  8119407.8125,
  8119784.286486917,
  8120111.389292315,
  8120337.5,
  8121123.741568458,
  8121775.46166475,
  8122257.300305839,
  8122411.211256543,
  8123620.902497167,
  8124050.235486784,
  8124765.245954152,
  8125026.876584971,
  8125293.129518989,
  8126296.875,
  8126345.3125,
  8126849.679525267,
  8126854.570839319,
  8127434.257318617,
  8127610.301968664,
  8128060.806240722,
  8128788.326300505,
  8130394.441011572,
  8148213.394377924,
  8154723.1846909765,
  8164401.5625,
  8166102.663778551,
  8192096.875,
  8201456.25,
  8201479.758199686,
  8201501.5625,
  8210608.064887031,
  8220403.651414304,
  8237193.75,
  8279245.3125,
  8279327.193462723,
  8281749.511716898,
  8282440.6902050655,
  8282845.3125,
  8282847.637554568,
  8282851.5625,
  8282851.682549346,
  8283321.875,
  8283490.302014268,
  8293037.5,
  8431257.6909584,
  8444651.574392842,
  8518851.419710709,
  8543912.817068344,
  8659973.68822171,
  8671880.570053363,
  10725091.041000225,
  10814141.783209614,
  10814938.866117498,
  10816056.25,
  10816510.66763128,
  10817462.353888037,
  10822502.389516778,
  10823751.5625,
  10844075.80897423,
  10845218.187316788,
  10845837.407292234,
  10846184.345855124,
  10847018.75,
  11805764.029285511,
  11921743.75,
  11949632.435519068,
  11950119.943488978,
  11995270.461129945,
  12021439.0625,
  12021790.213106142,
  12021903.76831068,
  12023928.058359113,
  12029901.5625,
  12032356.25,
  12032566.330043258,
  12033127.680007346,
  12033213.651444098,
  12033272.094065566,
  12033376.625534093,
  12033396.719807995,
  12033728.125,
  12033768.75,
  12033795.942995561,
  12033824.607349725,
  12034171.185366008,
  12034991.840709705,
  12035120.3125,
  12035567.1875,
  12036231.607004987,
  12037243.75,
  12045588.584479712,
  12048837.5,
  12050963.139607713,
  12058795.3125,
  12762587.5,
  12956062.554147929,
  13053614.281567069,
  13068352.457186054,
  13075096.875,
  13075639.945638333,
  13077987.5,
  13086819.859345548,
  13086821.652011884,
  13089836.260993473,
  13091025.479852986,
  13103281.93071233,
  13138178.125,
  13141278.125,
  13170725.91664031,
  13595979.993300736,
  13882415.300789392,
  16412244.98934245,
  16414216.73666595,
  16415442.1875,
  16416021.875,
  16416364.017114962,
  16416374.618944822,
  16416996.875,
  16417047.034869852,
  16417461.81319017,
  16420471.109906167,
  16420676.16772352,
  16421679.301742032,
  16421807.8125,
  16423380.131742321,
  16423774.339786593,
  16424723.008657621,
  16428105.59359269,
  16428929.6875,
  16432030.922882574,
  16437589.0625,
  16438021.201651422,
  16438662.37073636,
  16438719.414961193,
  16442312.300715653,
  16442533.38857625,
  16443037.5,
  16446536.208341783,
  16446818.70456963,
  16447123.4375,
  16448395.977873702,
  16454951.51280248],
 [7.01948642501036,
  66.5344752909091,
  60.31015088691055,
  58.807249870126746,
  44.03904779972045,
  59.37580299397437,
  11.889464342802523,
  37.13071448641459,
  45.1374148155577,
  5.599248888212579,
  12.977282600400674,
  58.27258785089079,
  34.82352582445307,
  21.62640621979545,
  46.920457061113126,
  27.71995546625717,
  27.423223830130336,
  19.919517997404416,
  94.05985280775124,
  6.434001922751892,
  7.763790653072909,
  190.74380114684226,
  46.35462406180823,
  24.905573582526053,
  51.52242626319719,
  11.521093839300718,
  8.220276325770298,
  27.074976173999428,
  5.7575248619425,
  13.736160203447916,
  7.0893454247500065,
  12.180726844509678,
  9.181178514383564,
  11.499558965669578,
  93.11351904384259,
  15.263719875607402,
  58.31284022111421,
  11.937019493941452,
  5.082390729549531,
  45.39026707473275,
  59.635129428750794,
  46.19213641925138,
  63.528374628132724,
  72.74818997139262,
  51.54304523246065,
  10.681349172057688,
  8.895443407269585,
  7.0856130415657494,
  26.70561811891589,
  15.091743275262589,
  93.65821982617186,
  5.354390188730218,
  92.70727539858609,
  40.220045566177134,
  7.9147349230221895,
  29.786147033160134,
  23.73013127136719,
  45.14792189579634,
  54.48486962159565,
  17.40342296553851,
  6.116479587207561,
  10.183711728245376,
  15.565793254460981,
  52.715997787074635,
  37.79805709352883,
  5.723737457916772,
  35.363379082082126,
  67.91281922761925,
  22.041392775048468,
  24.596528704896908,
  11.560974923742211,
  14.28489046241263,
  19.828347753205705,
  5.913140505281145,
  30.84927515445939,
  24.26478682847881,
  15.80757787538496,
  30.975052760046104,
  7.847658584364835,
  27.65711666150794,
  88.94033302056516,
  13.35972862675263,
  29.254265501214082,
  18.806178073471457,
  103.21065893487528,
  88.21821869322318,
  27.494018652988636,
  69.08801818417956,
  98.64406292105221,
  5.450575079151153,
  6.762355996254278,
  59.34291309079271,
  66.99383363677563,
  5.505124699710122,
  45.204167937297555,
  38.001460838548276,
  67.32652069615425,
  15.539450498463784,
  47.39848625237699,
  20.75853977662154,
  26.36428938027978,
  6.514411803077918,
  77.26025906445722,
  31.153081766231978,
  14.010404894940924,
  81.78434290161276,
  90.36293292109438,
  17.587385091516317,
  90.29259681613227,
  38.67964632334759,
  81.81294045640266,
  78.65935000452662,
  111.1872621520416,
  74.43634351252643,
  14.06827609612685,
  42.99775472276282,
  47.129850063882856,
  32.37198239892718,
  8.077418413775064,
  77.92308264200847,
  35.235842860851704,
  74.71835895380818,
  35.10502150440469,
  66.98134535188979,
  12.438428781658171,
  57.49066791707384,
  22.80013103028026,
  17.52453073898031,
  29.910045229436932,
  119.65984004010689,
  5.166343226460405,
  22.16047291462894,
  18.806627203819815,
  56.37830779705009,
  126.07292485672849,
  33.40471336476145,
  93.43805855070607,
  151.4559803758041,
  29.199491004698487,
  11.952120088919054,
  50.67350034480911,
  59.15314927195181,
  16.944692324623,
  12.03351395896264,
  32.560997706221364,
  22.157391524439905,
  9.571731739585307,
  22.319217459656556,
  23.644373979877276,
  7.515477096229745,
  39.81666046625503,
  6.097979143852031,
  8.81772402111189,
  10.007909185576217,
  6.524132495346957,
  136.0984072982323,
  20.302643101868046,
  97.37256540234979,
  39.9675556136904,
  5.738939188561282,
  9.577840303031515,
  5.470211823725413,
  16.271294710365076,
  21.560306464015177,
  27.033278559497933,
  21.513171879898294,
  62.996352305148434,
  77.2052304593661,
  11.326454062750171,
  118.43508598433738,
  7.3349376306266425,
  34.32894286558482,
  16.587694927966236,
  35.48357472115026,
  6.103820333534942,
  140.50183262448172,
  113.1949923592802,
  15.969823229073276,
  95.83185575224458,
  108.79473019544928,
  6.243910737215056,
  149.509639765815,
  95.81980058999328,
  71.34207437810113,
  56.45603567386273,
  39.82854734342361,
  55.34647400741533,
  8.121521393567496,
  111.53684735254822,
  8.255221518407486,
  51.023347804451,
  52.393514509451975,
  15.183685461573914,
  13.720195522184644,
  120.34628517000759,
  86.23788257531342,
  62.6340797745629,
  50.57864844152558,
  88.37294311637322,
  7.625046506672427,
  7.819042973771666,
  11.477175512841768,
  39.696999397822424,
  21.898632992768597,
  10.41876176278347,
  55.897453797976034,
  106.6211841686959,
  28.30799429673197,
  19.979940391073555,
  10.978404753405124,
  72.33804379425865,
  19.440692424021194,
  17.2629011687243,
  14.892584152340246,
  67.9990621148876,
  25.93454057980983,
  94.91739921347983,
  6.744870477656062,
  47.88189341391646,
  74.44776034856416,
  8.302529041937728,
  40.73971303264943,
  224.0079844782543,
  47.78312194934922,
  14.601408600712041,
  30.581149759162695,
  30.389216439980547,
  48.2979307383526,
  50.487077091111054,
  51.16086882617536,
  7.875602407220497,
  31.201754631804178,
  13.784890367903143,
  10.308187179208561,
  18.319028179826905,
  36.40186767606083,
  23.167850728245455,
  18.05838450240185,
  65.59209701472635,
  30.049248143749118,
  34.49003119887807,
  94.78521351051369,
  58.596978421854814,
  15.014990136100439,
  11.431789746274037,
  17.413647720303224,
  209.50872831091377,
  28.06506262159766,
  94.40611345419259,
  48.228949350596004,
  14.25914984041219,
  88.15775446547701,
  44.489907098302666,
  36.52467891332402,
  6.379607413714642,
  15.975034794336581,
  57.3329413597779,
  11.084475160518565,
  49.578195102523566,
  22.904101107776782,
  40.78986831271144,
  29.95084325543879,
  115.0640568074223,
  65.47489978498575,
  22.066768495516133,
  20.402444136158262,
  10.638110414522798,
  17.103848046619827,
  6.88696432167989,
  14.831409765521775,
  19.37674956878004,
  61.380033779791745,
  53.60163648888627,
  20.41568583529521,
  65.5778558280526,
  49.78491419066978,
  76.25758538163109,
  168.3829980692553,
  52.1180371437178,
  28.11794278661078,
  17.089728502558078,
  21.358636212983704,
  155.1979471989042,
  27.04226193325171,
  19.042624472018247,
  27.29147294117218,
  49.381286963217555,
  160.245494998116,
  30.895713893257824,
  10.670253439284917,
  20.88749173370106,
  37.4522393876569,
  44.70160904055175,
  27.751877502028695,
  51.684699814098614,
  13.784109211911154,
  136.3354759770187,
  106.00913812970329,
  19.439721314241293,
  9.34310688158804,
  53.073200540790225,
  10.532256861187141,
  23.264906042686434,
  17.685943543814037,
  64.92706328246186,
  87.68537595769766,
  50.49535000610757,
  24.811255717056913,
  19.978375402761404,
  15.960292882849512,
  17.98589602723802,
  83.82272723815153,
  6.148242591149478,
  54.74640876897011,
  30.174629157741613,
  6.983895201946333,
  26.156925450808387,
  15.4595102923306,
  40.72108683395517,
  12.530019091545823,
  93.22085259295635,
  90.62572792850523,
  9.177713572442869,
  99.69467097176545,
  91.6322453520697,
  21.30538608761715,
  65.57212150464937,
  24.936651580767368,
  17.94623503277572,
  62.033853505650875,
  5.224728526227363,
  83.36973492141257,
  76.3338083851332,
  36.668112067345135,
  17.98113415751144,
  70.53967302145612,
  71.45082073703959,
  14.114306917976986,
  65.76645680316268,
  12.057254425441004,
  68.03078799304134,
  28.619368815857676,
  32.2342117444501,
  86.62283885545854,
  17.05594592248561,
  32.88571888912515,
  9.717971190250434,
  39.884790934033774,
  87.04908712868436,
  58.060989958416485,
  63.14092819914653,
  34.27977224844754,
  5.716747544554541,
  65.56672993227527,
  38.32799547832274,
  22.826377473400147,
  17.41731766624705,
  82.12040332977067,
  48.07284305482992,
  18.433259980650668,
  133.33811851050376,
  103.93038694304727,
  46.30537083409534,
  35.73305287643594,
  66.25242789247426,
  79.29784991146443,
  10.849926253103032,
  16.90234586875447,
  60.633656655178385,
  28.752072152929404,
  24.498565015031357,
  103.67552141114194,
  30.150690250141423,
  21.750404567447887,
  6.847799763690607,
  17.528318592644773,
  31.034685844103823,
  11.996567471438304,
  45.498316296673,
  153.25420061452704,
  144.8597492405272,
  10.064211187533502,
  15.993892629083188,
  9.977616029560643,
  160.6320917151447,
  15.591445904642253,
  141.09952847124708,
  28.308173776293337,
  16.27334196853476,
  76.03808418491538,
  9.020068670224314,
  16.715084763437552,
  28.294177136650532,
  168.04138273387548,
  27.600518401003228,
  19.466610822228752,
  26.042167706909794,
  28.475501963253397,
  48.67486279225878,
  40.5105835016637,
  35.44033212604757,
  47.25530165951515,
  32.358615543203015,
  43.616184286091375,
  17.304429020167554,
  11.43838791393536,
  8.222900386474333,
  99.22619737560433,
  31.311807509626046,
  63.68146441032566,
  41.063377500850876,
  17.47926567105375,
  9.272758695425004,
  82.96550919229492,
  80.01576952123023,
  27.514415917045337,
  84.6703510835849,
  41.24248792088897,
  9.730131270595315,
  8.766677985748457,
  111.40842285845807,
  63.75469372622584,
  25.154042451332636,
  13.591565044220127,
  39.310501605966834,
  14.400183113800788,
  60.33421256013218,
  342.2368594861819,
  21.269915192104396,
  28.714816815368398,
  48.673591491272134,
  29.500433849373213,
  68.43004407218872,
  7.068152323284567,
  38.969453627299565,
  84.64689377389168,
  12.12968255252203,
  68.42734821563062,
  24.576602932665264,
  7.910889176187939,
  10.489930692405192,
  75.07853181726583,
  6.698862382409513,
  131.60487806034422,
  59.605057640178224,
  23.387563555654168,
  5.688899929952562,
  19.841634402843326,
  38.624652918136434,
  73.52091928047172,
  49.79276791155743,
  20.564115409617752,
  18.49774427546654,
  157.19093796273404,
  17.415993167074536,
  32.07119962571495,
  8.558218589611204,
  5.756528691766511,
  33.509400962997184,
  55.752147825836175,
  36.84681638070083,
  135.6936258451777,
  45.83424469415289,
  165.33858048887268,
  29.033364060284192,
  30.543214601143298,
  11.483252039375941,
  57.70960286832259,
  151.1827831237001,
  160.21111183262104,
  147.96336756459135,
  12.50296055739069,
  27.582708274366365,
  13.812566580984075,
  19.548654512312712,
  115.13868251898653,
  14.798156623344846,
  47.02515780789199,
  59.87398114485782,
  20.176143891995384,
  5.920306493509064,
  41.67058821352796,
  71.9424026223688,
  50.502555018418754,
  89.20243131304137,
  59.73242057507242,
  41.268084508405074,
  13.479888285782101,
  9.8084152873118,
  10.533662159678805,
  33.810940246315376,
  173.6890138636346,
  19.653637438452986,
  14.901945593704594,
  178.4985370515093,
  31.006518402273013,
  185.0410946524026,
  23.166356589444973,
  128.94315369225816,
  164.98825986850855,
  488.81047672436205,
  13.292772709478893,
  5.834517441355705,
  12.194814287675038,
  12.848638265611687,
  33.59336376401455,
  33.18862946613102,
  27.211577708324278,
  33.9285353726062,
  72.89399692287354,
  9.07436784130068,
  5.287497241740275,
  9.233840478109595,
  56.96200469664292,
  8.218363860710046,
  10.143708939945652,
  71.31742272978735,
  91.26734705277974,
  24.372857176076362,
  39.7621301569106,
  10.475971581836353,
  83.1366952685337,
  30.678923665187845,
  15.218009317075259,
  59.70791589325795,
  88.37888095748049,
  31.270848558495462,
  44.27785002988454,
  106.35383576934544,
  54.70533148395681,
  82.49152650932639,
  46.46836632088982,
  49.95622134319276,
  26.33377504486552,
  31.13012169604691,
  7.647242348643645,
  87.70916733128635,
  244.3491998399123,
  44.18998447137148,
  27.031806457671486,
  11.14634710319366,
  71.22879186649756,
  11.662318379562311,
  15.991783934276732,
  19.868060324019236,
  137.22063823086242,
  42.38645826180944,
  14.242339261180264,
  39.43918152336181,
  8.148842273041652,
  10.90675470673964,
  21.337718957501526,
  72.14724480940575,
  105.48717427939776,
  20.88979017664333,
  26.816781512226054,
  35.31175626525021,
  72.18762454173917,
  32.77324708509002,
  43.426352111209674,
  6.220716156886273,
  60.09042073058847,
  8.366770954639204,
  29.119503034392448,
  22.419461717742248,
  70.90023495316017,
  12.543898943657375,
  51.84963440419426,
  36.73649390827758,
  30.99849900974961,
  5.442985697440052,
  14.24141571826322,
  15.847843284458689,
  8.320194482960225,
  56.50965716339728,
  37.40800654950942,
  38.56510713587179,
  53.134957484725234,
  34.35026608747548,
  10.197422405449903,
  23.170910854904044,
  63.732374890056654,
  64.43720126789721,
  24.477392453311655,
  65.25096068944183,
  74.38590144420631,
  108.6739682549958,
  18.107091980841762,
  102.8265977671968,
  58.39070447933071,
  42.01451500795805,
  9.523770418462753,
  6.95631268600118,
  31.62766201396243,
  25.244079148727934,
  74.71889691396936,
  31.236141190981193,
  24.786967037988923,
  6.122398380105475,
  5.417140765627537,
  14.461384871439568,
  78.53648675725556,
  42.77447559141699,
  28.86931685524332,
  60.99998310461749,
  6.227087633799472,
  107.47634152240818,
  14.095113623373695,
  88.59071709566189,
  41.840726521834654,
  27.657580766254853,
  22.767794584542173,
  73.48899211355393,
  12.114093641176279,
  33.7663723096822,
  26.54422324200455,
  22.61116499362535,
  22.713810400142393,
  139.77593730435552,
  9.737580140139878,
  20.031078937212314,
  126.84831358045426,
  20.06291603638016,
  39.15799931658964,
  9.069319639681057,
  48.26116806569466,
  90.16281308540533,
  46.14274731649472,
  33.22964088097305,
  14.23269848518827,
  191.169578272674,
  108.14504140695574,
  7.704224694560673,
  14.771132032787731,
  237.37043882870378,
  5.75045442817992,
  11.050904857933716])

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

([3516563.659866176,
  3522453.125,
  3616325.0,
  3650026.5625,
  3723081.982589908,
  3761592.036803253,
  3849862.938197821,
  3887018.75,
  3895548.4375,
  3895574.7294298923,
  3898183.742516835,
  3900555.2622644934,
  3901127.8228500597,
  3905273.4146939702,
  3913945.3125,
  3914570.7848742297,
  3918353.2176982756,
  3958741.5138268527,
  3970011.4385161516,
  3979799.652759667,
  3985645.7356110835,
  3994223.9795589126,
  4001441.8845774713,
  4033557.597831782,
  4058732.589787411,
  4067925.2378966655,
  4069324.8724835906,
  4086329.609696547,
  4109671.081690553,
  4114940.2682108535,
  4117817.0639380845,
  4125131.833253878,
  4125265.441827439,
  4131237.725549722,
  4131240.625,
  4150953.2070430163,
  4165110.9375,
  4189198.9001408773,
  4189249.166159348,
  4189553.8027579724,
  4189646.844598682,
  4192839.0507922294,
  4196190.625,
  4212550.0,
  4214835.9375,
  4218457.50684519,
  4221790.060107657,
  4223167.049386633,
  4223593.65076936,
  4236418.788865045,
  4241158.793132736,
  4246851.921913718,
  4260638.556822054,
  4262944.518657322,
  4263064.848772324,
  4263129.9123975225,
  4264238.760884773,
  4264580.94333534,
  4269443.575760565,
  4287044.235263441,
  4320328.393560219,
  4320855.548561205,
  4322767.342182601,
  4323065.625,
  4323229.6875,
  4323620.57602832,
  4329303.125,
  4333546.875,
  4336798.1013282575,
  4336802.888449427,
  4340942.448067964,
  4378198.25385809,
  4379257.716371073,
  4439030.943377524,
  4477696.875,
  4478680.133732051,
  4478698.647567996,
  4479185.133383408,
  4484549.066246027,
  4485729.642089781,
  4485742.44386866,
  4534054.199752554,
  4551546.875,
  4559218.112855378,
  4563014.0625,
  4564523.4375,
  4564706.823085847,
  4565094.13595249,
  4566114.0625,
  4566659.11163613,
  4581139.056734097,
  4593009.375,
  4594487.5,
  4633616.59770668,
  4635914.0625,
  4638708.995465376,
  4667414.0625,
  4669028.273503867,
  4670405.79649102,
  4683790.372460848,
  4753434.255329365,
  4754104.859598027,
  4767010.798104672,
  4786629.6875,
  4787223.169701784,
  4787296.81704278,
  4789484.375,
  4805927.3106456995,
  4835882.8125,
  4852521.875,
  4853962.457951159,
  4901526.5625,
  4901549.697415754,
  4910900.0,
  4911341.937839819,
  4930861.602936756,
  4961820.3125,
  4964982.8125,
  4966923.019864093,
  5055995.609911848,
  5056321.406174996,
  5057182.8125,
  5057650.0,
  5058361.203861539,
  5059817.890399953,
  5060273.4375,
  5060786.064762945,
  5072169.356117767,
  5107959.710756048,
  5161585.9375,
  5175995.354599998,
  5177107.376896036,
  5178205.430695763,
  5178362.5,
  5180024.552726621,
  5190468.888403881,
  5192246.14109956,
  5239840.625,
  5248677.229435093,
  5248686.407838635,
  5248691.690909987,
  5287351.5625,
  5294191.06705885,
  5294198.948810237,
  5294745.3125,
  5328335.50160979,
  5350006.521177989,
  5361277.704009917,
  5361306.943817168,
  5431927.875656571,
  5450353.56765452,
  5450434.4694744935,
  5465775.885359852,
  5477467.892138546,
  5477475.232293967,
  5480420.3125,
  5507400.907576567,
  5539351.247810291,
  5540832.8125,
  5569589.556552177,
  5637797.152746597,
  5674886.132797357,
  5676977.644666485,
  5709127.631710819,
  5839211.621115992,
  6506453.106326316,
  6514916.316738477,
  6527396.51166487,
  6528086.061037031,
  6528271.067004249,
  6528316.379315424,
  6528324.667390298,
  6528546.182630836,
  6528578.102955911,
  6528600.149832682,
  6671662.5,
  6805143.75,
  6805822.120942643,
  6819506.25,
  6820201.955834646,
  6820232.625928741,
  6822763.93460819,
  6908504.296530955,
  6928553.125,
  6929986.028857499,
  6932011.688049375,
  7070934.375,
  7070946.549266327,
  7184259.200758169,
  7185426.502731286,
  7185976.5625,
  7193747.411353213,
  7328268.812119185,
  7328951.200370186,
  7331853.125,
  7334797.1820328785,
  7337451.5625,
  7346221.875,
  7376234.375,
  7376248.497398315,
  7385063.05969682,
  7388484.587637318,
  7430701.5625,
  7457791.727214021,
  7457888.556073201,
  7459475.0,
  7462050.0,
  7462079.330295611,
  7466277.458640598,
  7479812.402334075,
  7480492.1875,
  7484383.5667245835,
  7487419.249669934,
  7489599.98710439,
  7492492.093171464,
  7493315.7290560575,
  7499106.96147352,
  7500438.968430275,
  7500525.883778526,
  7500663.9830243075,
  7500664.092543603,
  7501518.75,
  7501529.964638061,
  7502032.345091607,
  7502081.157814215,
  7502651.877024314,
  7502660.37610205,
  7502665.625,
  7502792.1875,
  7503533.725369856,
  7503871.007596663,
  7504671.739338248,
  7504673.096900943,
  7504750.966245916,
  7504963.172085897,
  7505013.184626922,
  7505110.59045512,
  7505275.632001752,
  7505277.89805831,
  7505306.25,
  7505354.6875,
  7505392.1875,
  7505503.125,
  7505539.081645869,
  7505595.49870544,
  7505752.272668069,
  7505762.358975728,
  7505903.1013427805,
  7505973.4375,
  7506116.900846033,
  7506243.558663944,
  7506325.0,
  7506437.5,
  7506442.887391759,
  7506475.840766794,
  7506678.713185535,
  7506785.814736297,
  7506843.722164545,
  7507137.5,
  7507309.005697363,
  7507353.125,
  7507767.310103823,
  7508012.48385256,
  7508507.640571699,
  7508546.430646867,
  7508647.820744732,
  7509394.974935964,
  7510080.158732645,
  7510089.650962341,
  7510096.845811033,
  7510366.4454933675,
  7511011.988574928,
  7511096.875,
  7511950.947257247,
  7513034.781752089,
  7513100.0,
  7513442.1875,
  7513553.369248716,
  7513650.0,
  7514315.16715968,
  7515747.343673917,
  7516033.86141414,
  7516110.248588808,
  7516182.732826867,
  7525296.51479232,
  7526988.996935411,
  7527384.375,
  7527796.875,
  7558241.802918784,
  7558776.248698344,
  7561118.276115736,
  7567262.806755644,
  7576693.75,
  7578590.5847103195,
  7621339.41884141,
  7621541.846642729,
  7621578.827208772,
  7622125.0,
  7622158.158586307,
  7623459.388809071,
  7623473.4375,
  7623816.769938677,
  7623865.28598947,
  7625083.38450811,
  7626881.25,
  7627787.5,
  7629726.325955533,
  7669481.028723511,
  7669493.083451879,
  7714184.61886304,
  7715714.500899253,
  7721634.375,
  7725970.871652735,
  7725998.7184672905,
  7727196.875,
  7728661.949523789,
  7731865.492599414,
  7739229.817074609,
  7754009.056972649,
  7758397.618610811,
  7758751.5625,
  7776972.935246239,
  7801048.249968431,
  7830308.459139842,
  7843560.9375,
  7845811.575901804,
  7850957.8125,
  7850977.587026083,
  7851108.426781962,
  7851614.0625,
  7857219.566250149,
  7858342.014396781,
  7861776.5625,
  7871777.579991633,
  7872139.766202867,
  7874850.0,
  7878893.75,
  7878911.689775805,
  7880784.344459365,
  7880838.023097858,
  7885476.5625,
  7885601.5625,
  7887880.148248193,
  7889028.1478606,
  7891131.712528821,
  7900521.939498802,
  7907994.174336239,
  7908909.375,
  7909715.581591999,
  7914585.9375,
  7917464.424186572,
  7919542.520877404,
  7919951.718326982,
  7920125.0,
  7920556.129278812,
  7920565.249017154,
  7930822.401616633,
  7931471.218914138,
  7935716.83438859,
  7935820.673969362,
  7936287.5,
  7937709.375,
  7938704.6875,
  7939506.675765278,
  7941996.636586881,
  7942496.597864096,
  7942515.4397056755,
  7945003.847510337,
  7950767.163955314,
  7951167.1875,
  7953253.297134905,
  7954354.248276682,
  7954564.61794186,
  7954773.3243605485,
  7954909.783049829,
  7954950.545597302,
  7955018.203875645,
  7955091.117820228,
  7955109.375,
  7955382.8125,
  7955549.207529363,
  7956063.255779153,
  7956654.938679654,
  7956731.9267000845,
  7956804.3683185885,
  7956948.598171674,
  7957244.0949284285,
  7957581.588556634,
  7957622.185436127,
  7957793.75,
  7958537.107324194,
  7958808.211149185,
  7959150.502406015,
  7959260.569680303,
  7959454.274185324,
  7959469.582145377,
  7959478.40863012,
  7960003.122486739,
  7960014.0625,
  7960076.5625,
  7961015.625,
  7961285.766834207,
  7963045.3125,
  7963910.9375,
  7968078.631979903,
  7968616.648323349,
  7978218.627613123,
  7978846.875,
  7979496.875,
  7980857.8125,
  7981467.502001785,
  7982330.224749004,
  7983459.304819967,
  7984481.25,
  7984924.764752028,
  7985193.696799739,
  7985203.125,
  7985332.8125,
  7985451.999197627,
  7985938.309492799,
  7986282.8125,
  7986432.508787645,
  7986596.0434040325,
  7986713.6566671375,
  7986745.3125,
  7987198.287837345,
  7987228.125,
  7987928.768118402,
  7988132.289831875,
  7988139.0625,
  7988231.25,
  7988332.8125,
  7988335.809714023,
  7988348.230784366,
  7988459.375,
  7988854.6875,
  7988875.558616115,
  7988876.314511786,
  7990432.574515623,
  7991405.464900492,
  7991929.237181021,
  7992604.6875,
  7992660.604839769,
  7992774.346633762,
  7992804.6875,
  7992811.072266713,
  7992923.153195937,
  7992981.090894862,
  7993001.41287108,
  7993037.117982465,
  7993169.199880819,
  7993277.526542489,
  7993702.261217761,
  7993907.4795841705,
  7994029.21684611,
  7994227.833224316,
  7994539.543126533,
  7995877.317126113,
  7995916.748147555,
  8000692.1875,
  8011682.8125,
  8014434.375,
  8015899.054810305,
  8025685.9375,
  8032218.1172520835,
  8047057.080455237,
  8047248.653650699,
  8050101.520064619,
  8050116.534531348,
  8053466.341432249,
  8054517.039802403,
  8059988.231678637,
  8060043.700539757,
  8061199.668521936,
  8061385.217663492,
  8061429.6875,
  8063539.024374357,
  8064121.888140827,
  8065523.4375,
  8067077.707971298,
  8067296.653043272,
  8068217.1875,
  8075822.911583884,
  8075920.3125,
  8077877.847967238,
  8082654.6875,
  8085544.452148123,
  8113979.242727022,
  8114130.202462289,
  8114766.771018609,
  8117300.3447585935,
  8117887.355749825,
  8118154.107213509,
  8118401.552081097,
  8119190.625,
  8119407.8125,
  8119784.286486917,
  8120111.389292315,
  8120337.5,
  8121123.741568458,
  8121775.46166475,
  8122257.300305839,
  8122411.211256543,
  8123620.902497167,
  8124050.235486784,
  8124765.245954152,
  8125026.876584971,
  8125293.129518989,
  8126296.875,
  8126345.3125,
  8126849.679525267,
  8126854.570839319,
  8127434.257318617,
  8127610.301968664,
  8128060.806240722,
  8128788.326300505,
  8130394.441011572,
  8148213.394377924,
  8154723.1846909765,
  8164401.5625,
  8166102.663778551,
  8192096.875,
  8201456.25,
  8201479.758199686,
  8201501.5625,
  8210608.064887031,
  8220403.651414304,
  8237193.75,
  8279245.3125,
  8279327.193462723,
  8281749.511716898,
  8282440.6902050655,
  8282845.3125,
  8282847.637554568,
  8282851.5625,
  8282851.682549346,
  8283321.875,
  8283490.302014268,
  8293037.5,
  8431257.6909584,
  8444651.574392842,
  8518851.419710709,
  8543912.817068344,
  8659973.68822171,
  8671880.570053363,
  10725091.041000225,
  10814141.783209614,
  10814938.866117498,
  10816056.25,
  10816510.66763128,
  10817462.353888037,
  10822502.389516778,
  10823751.5625,
  10844075.80897423,
  10845218.187316788,
  10845837.407292234,
  10846184.345855124,
  10847018.75,
  11805764.029285511,
  11921743.75,
  11949632.435519068,
  11950119.943488978,
  11995270.461129945,
  12021439.0625,
  12021790.213106142,
  12021903.76831068,
  12023928.058359113,
  12029901.5625,
  12032356.25,
  12032566.330043258,
  12033127.680007346,
  12033213.651444098,
  12033272.094065566,
  12033376.625534093,
  12033396.719807995,
  12033728.125,
  12033768.75,
  12033795.942995561,
  12033824.607349725,
  12034171.185366008,
  12034991.840709705,
  12035120.3125,
  12035567.1875,
  12036231.607004987,
  12037243.75,
  12045588.584479712,
  12048837.5,
  12050963.139607713,
  12058795.3125,
  12762587.5,
  12956062.554147929,
  13053614.281567069,
  13068352.457186054,
  13075096.875,
  13075639.945638333,
  13077987.5,
  13086819.859345548,
  13086821.652011884,
  13089836.260993473,
  13091025.479852986,
  13103281.93071233,
  13138178.125,
  13141278.125,
  13170725.91664031,
  13595979.993300736,
  13882415.300789392,
  16412244.98934245,
  16414216.73666595,
  16415442.1875,
  16416021.875,
  16416364.017114962,
  16416374.618944822,
  16416996.875,
  16417047.034869852,
  16417461.81319017,
  16420471.109906167,
  16420676.16772352,
  16421679.301742032,
  16421807.8125,
  16423380.131742321,
  16423774.339786593,
  16424723.008657621,
  16428105.59359269,
  16428929.6875,
  16432030.922882574,
  16437589.0625,
  16438021.201651422,
  16438662.37073636,
  16438719.414961193,
  16442312.300715653,
  16442533.38857625,
  16443037.5,
  16446536.208341783,
  16446818.70456963,
  16447123.4375,
  16448395.977873702,
  16454951.51280248],
 [7.01948642501036,
  66.5344752909091,
  60.31015088691055,
  58.807249870126746,
  44.03904779972045,
  59.37580299397437,
  11.889464342802523,
  37.13071448641459,
  45.1374148155577,
  5.599248888212579,
  12.977282600400674,
  58.27258785089079,
  34.82352582445307,
  21.62640621979545,
  46.920457061113126,
  27.71995546625717,
  27.423223830130336,
  19.919517997404416,
  94.05985280775124,
  6.434001922751892,
  7.763790653072909,
  190.74380114684226,
  46.35462406180823,
  24.905573582526053,
  51.52242626319719,
  11.521093839300718,
  8.220276325770298,
  27.074976173999428,
  5.7575248619425,
  13.736160203447916,
  7.0893454247500065,
  12.180726844509678,
  9.181178514383564,
  11.499558965669578,
  93.11351904384259,
  15.263719875607402,
  58.31284022111421,
  11.937019493941452,
  5.082390729549531,
  45.39026707473275,
  59.635129428750794,
  46.19213641925138,
  63.528374628132724,
  72.74818997139262,
  51.54304523246065,
  10.681349172057688,
  8.895443407269585,
  7.0856130415657494,
  26.70561811891589,
  15.091743275262589,
  93.65821982617186,
  5.354390188730218,
  92.70727539858609,
  40.220045566177134,
  7.9147349230221895,
  29.786147033160134,
  23.73013127136719,
  45.14792189579634,
  54.48486962159565,
  17.40342296553851,
  6.116479587207561,
  10.183711728245376,
  15.565793254460981,
  52.715997787074635,
  37.79805709352883,
  5.723737457916772,
  35.363379082082126,
  67.91281922761925,
  22.041392775048468,
  24.596528704896908,
  11.560974923742211,
  14.28489046241263,
  19.828347753205705,
  5.913140505281145,
  30.84927515445939,
  24.26478682847881,
  15.80757787538496,
  30.975052760046104,
  7.847658584364835,
  27.65711666150794,
  88.94033302056516,
  13.35972862675263,
  29.254265501214082,
  18.806178073471457,
  103.21065893487528,
  88.21821869322318,
  27.494018652988636,
  69.08801818417956,
  98.64406292105221,
  5.450575079151153,
  6.762355996254278,
  59.34291309079271,
  66.99383363677563,
  5.505124699710122,
  45.204167937297555,
  38.001460838548276,
  67.32652069615425,
  15.539450498463784,
  47.39848625237699,
  20.75853977662154,
  26.36428938027978,
  6.514411803077918,
  77.26025906445722,
  31.153081766231978,
  14.010404894940924,
  81.78434290161276,
  90.36293292109438,
  17.587385091516317,
  90.29259681613227,
  38.67964632334759,
  81.81294045640266,
  78.65935000452662,
  111.1872621520416,
  74.43634351252643,
  14.06827609612685,
  42.99775472276282,
  47.129850063882856,
  32.37198239892718,
  8.077418413775064,
  77.92308264200847,
  35.235842860851704,
  74.71835895380818,
  35.10502150440469,
  66.98134535188979,
  12.438428781658171,
  57.49066791707384,
  22.80013103028026,
  17.52453073898031,
  29.910045229436932,
  119.65984004010689,
  5.166343226460405,
  22.16047291462894,
  18.806627203819815,
  56.37830779705009,
  126.07292485672849,
  33.40471336476145,
  93.43805855070607,
  151.4559803758041,
  29.199491004698487,
  11.952120088919054,
  50.67350034480911,
  59.15314927195181,
  16.944692324623,
  12.03351395896264,
  32.560997706221364,
  22.157391524439905,
  9.571731739585307,
  22.319217459656556,
  23.644373979877276,
  7.515477096229745,
  39.81666046625503,
  6.097979143852031,
  8.81772402111189,
  10.007909185576217,
  6.524132495346957,
  136.0984072982323,
  20.302643101868046,
  97.37256540234979,
  39.9675556136904,
  5.738939188561282,
  9.577840303031515,
  5.470211823725413,
  16.271294710365076,
  21.560306464015177,
  27.033278559497933,
  21.513171879898294,
  62.996352305148434,
  77.2052304593661,
  11.326454062750171,
  118.43508598433738,
  7.3349376306266425,
  34.32894286558482,
  16.587694927966236,
  35.48357472115026,
  6.103820333534942,
  140.50183262448172,
  113.1949923592802,
  15.969823229073276,
  95.83185575224458,
  108.79473019544928,
  6.243910737215056,
  149.509639765815,
  95.81980058999328,
  71.34207437810113,
  56.45603567386273,
  39.82854734342361,
  55.34647400741533,
  8.121521393567496,
  111.53684735254822,
  8.255221518407486,
  51.023347804451,
  52.393514509451975,
  15.183685461573914,
  13.720195522184644,
  120.34628517000759,
  86.23788257531342,
  62.6340797745629,
  50.57864844152558,
  88.37294311637322,
  7.625046506672427,
  7.819042973771666,
  11.477175512841768,
  39.696999397822424,
  21.898632992768597,
  10.41876176278347,
  55.897453797976034,
  106.6211841686959,
  28.30799429673197,
  19.979940391073555,
  10.978404753405124,
  72.33804379425865,
  19.440692424021194,
  17.2629011687243,
  14.892584152340246,
  67.9990621148876,
  25.93454057980983,
  94.91739921347983,
  6.744870477656062,
  47.88189341391646,
  74.44776034856416,
  8.302529041937728,
  40.73971303264943,
  224.0079844782543,
  47.78312194934922,
  14.601408600712041,
  30.581149759162695,
  30.389216439980547,
  48.2979307383526,
  50.487077091111054,
  51.16086882617536,
  7.875602407220497,
  31.201754631804178,
  13.784890367903143,
  10.308187179208561,
  18.319028179826905,
  36.40186767606083,
  23.167850728245455,
  18.05838450240185,
  65.59209701472635,
  30.049248143749118,
  34.49003119887807,
  94.78521351051369,
  58.596978421854814,
  15.014990136100439,
  11.431789746274037,
  17.413647720303224,
  209.50872831091377,
  28.06506262159766,
  94.40611345419259,
  48.228949350596004,
  14.25914984041219,
  88.15775446547701,
  44.489907098302666,
  36.52467891332402,
  6.379607413714642,
  15.975034794336581,
  57.3329413597779,
  11.084475160518565,
  49.578195102523566,
  22.904101107776782,
  40.78986831271144,
  29.95084325543879,
  115.0640568074223,
  65.47489978498575,
  22.066768495516133,
  20.402444136158262,
  10.638110414522798,
  17.103848046619827,
  6.88696432167989,
  14.831409765521775,
  19.37674956878004,
  61.380033779791745,
  53.60163648888627,
  20.41568583529521,
  65.5778558280526,
  49.78491419066978,
  76.25758538163109,
  168.3829980692553,
  52.1180371437178,
  28.11794278661078,
  17.089728502558078,
  21.358636212983704,
  155.1979471989042,
  27.04226193325171,
  19.042624472018247,
  27.29147294117218,
  49.381286963217555,
  160.245494998116,
  30.895713893257824,
  10.670253439284917,
  20.88749173370106,
  37.4522393876569,
  44.70160904055175,
  27.751877502028695,
  51.684699814098614,
  13.784109211911154,
  136.3354759770187,
  106.00913812970329,
  19.439721314241293,
  9.34310688158804,
  53.073200540790225,
  10.532256861187141,
  23.264906042686434,
  17.685943543814037,
  64.92706328246186,
  87.68537595769766,
  50.49535000610757,
  24.811255717056913,
  19.978375402761404,
  15.960292882849512,
  17.98589602723802,
  83.82272723815153,
  6.148242591149478,
  54.74640876897011,
  30.174629157741613,
  6.983895201946333,
  26.156925450808387,
  15.4595102923306,
  40.72108683395517,
  12.530019091545823,
  93.22085259295635,
  90.62572792850523,
  9.177713572442869,
  99.69467097176545,
  91.6322453520697,
  21.30538608761715,
  65.57212150464937,
  24.936651580767368,
  17.94623503277572,
  62.033853505650875,
  5.224728526227363,
  83.36973492141257,
  76.3338083851332,
  36.668112067345135,
  17.98113415751144,
  70.53967302145612,
  71.45082073703959,
  14.114306917976986,
  65.76645680316268,
  12.057254425441004,
  68.03078799304134,
  28.619368815857676,
  32.2342117444501,
  86.62283885545854,
  17.05594592248561,
  32.88571888912515,
  9.717971190250434,
  39.884790934033774,
  87.04908712868436,
  58.060989958416485,
  63.14092819914653,
  34.27977224844754,
  5.716747544554541,
  65.56672993227527,
  38.32799547832274,
  22.826377473400147,
  17.41731766624705,
  82.12040332977067,
  48.07284305482992,
  18.433259980650668,
  133.33811851050376,
  103.93038694304727,
  46.30537083409534,
  35.73305287643594,
  66.25242789247426,
  79.29784991146443,
  10.849926253103032,
  16.90234586875447,
  60.633656655178385,
  28.752072152929404,
  24.498565015031357,
  103.67552141114194,
  30.150690250141423,
  21.750404567447887,
  6.847799763690607,
  17.528318592644773,
  31.034685844103823,
  11.996567471438304,
  45.498316296673,
  153.25420061452704,
  144.8597492405272,
  10.064211187533502,
  15.993892629083188,
  9.977616029560643,
  160.6320917151447,
  15.591445904642253,
  141.09952847124708,
  28.308173776293337,
  16.27334196853476,
  76.03808418491538,
  9.020068670224314,
  16.715084763437552,
  28.294177136650532,
  168.04138273387548,
  27.600518401003228,
  19.466610822228752,
  26.042167706909794,
  28.475501963253397,
  48.67486279225878,
  40.5105835016637,
  35.44033212604757,
  47.25530165951515,
  32.358615543203015,
  43.616184286091375,
  17.304429020167554,
  11.43838791393536,
  8.222900386474333,
  99.22619737560433,
  31.311807509626046,
  63.68146441032566,
  41.063377500850876,
  17.47926567105375,
  9.272758695425004,
  82.96550919229492,
  80.01576952123023,
  27.514415917045337,
  84.6703510835849,
  41.24248792088897,
  9.730131270595315,
  8.766677985748457,
  111.40842285845807,
  63.75469372622584,
  25.154042451332636,
  13.591565044220127,
  39.310501605966834,
  14.400183113800788,
  60.33421256013218,
  342.2368594861819,
  21.269915192104396,
  28.714816815368398,
  48.673591491272134,
  29.500433849373213,
  68.43004407218872,
  7.068152323284567,
  38.969453627299565,
  84.64689377389168,
  12.12968255252203,
  68.42734821563062,
  24.576602932665264,
  7.910889176187939,
  10.489930692405192,
  75.07853181726583,
  6.698862382409513,
  131.60487806034422,
  59.605057640178224,
  23.387563555654168,
  5.688899929952562,
  19.841634402843326,
  38.624652918136434,
  73.52091928047172,
  49.79276791155743,
  20.564115409617752,
  18.49774427546654,
  157.19093796273404,
  17.415993167074536,
  32.07119962571495,
  8.558218589611204,
  5.756528691766511,
  33.509400962997184,
  55.752147825836175,
  36.84681638070083,
  135.6936258451777,
  45.83424469415289,
  165.33858048887268,
  29.033364060284192,
  30.543214601143298,
  11.483252039375941,
  57.70960286832259,
  151.1827831237001,
  160.21111183262104,
  147.96336756459135,
  12.50296055739069,
  27.582708274366365,
  13.812566580984075,
  19.548654512312712,
  115.13868251898653,
  14.798156623344846,
  47.02515780789199,
  59.87398114485782,
  20.176143891995384,
  5.920306493509064,
  41.67058821352796,
  71.9424026223688,
  50.502555018418754,
  89.20243131304137,
  59.73242057507242,
  41.268084508405074,
  13.479888285782101,
  9.8084152873118,
  10.533662159678805,
  33.810940246315376,
  173.6890138636346,
  19.653637438452986,
  14.901945593704594,
  178.4985370515093,
  31.006518402273013,
  185.0410946524026,
  23.166356589444973,
  128.94315369225816,
  164.98825986850855,
  488.81047672436205,
  13.292772709478893,
  5.834517441355705,
  12.194814287675038,
  12.848638265611687,
  33.59336376401455,
  33.18862946613102,
  27.211577708324278,
  33.9285353726062,
  72.89399692287354,
  9.07436784130068,
  5.287497241740275,
  9.233840478109595,
  56.96200469664292,
  8.218363860710046,
  10.143708939945652,
  71.31742272978735,
  91.26734705277974,
  24.372857176076362,
  39.7621301569106,
  10.475971581836353,
  83.1366952685337,
  30.678923665187845,
  15.218009317075259,
  59.70791589325795,
  88.37888095748049,
  31.270848558495462,
  44.27785002988454,
  106.35383576934544,
  54.70533148395681,
  82.49152650932639,
  46.46836632088982,
  49.95622134319276,
  26.33377504486552,
  31.13012169604691,
  7.647242348643645,
  87.70916733128635,
  244.3491998399123,
  44.18998447137148,
  27.031806457671486,
  11.14634710319366,
  71.22879186649756,
  11.662318379562311,
  15.991783934276732,
  19.868060324019236,
  137.22063823086242,
  42.38645826180944,
  14.242339261180264,
  39.43918152336181,
  8.148842273041652,
  10.90675470673964,
  21.337718957501526,
  72.14724480940575,
  105.48717427939776,
  20.88979017664333,
  26.816781512226054,
  35.31175626525021,
  72.18762454173917,
  32.77324708509002,
  43.426352111209674,
  6.220716156886273,
  60.09042073058847,
  8.366770954639204,
  29.119503034392448,
  22.419461717742248,
  70.90023495316017,
  12.543898943657375,
  51.84963440419426,
  36.73649390827758,
  30.99849900974961,
  5.442985697440052,
  14.24141571826322,
  15.847843284458689,
  8.320194482960225,
  56.50965716339728,
  37.40800654950942,
  38.56510713587179,
  53.134957484725234,
  34.35026608747548,
  10.197422405449903,
  23.170910854904044,
  63.732374890056654,
  64.43720126789721,
  24.477392453311655,
  65.25096068944183,
  74.38590144420631,
  108.6739682549958,
  18.107091980841762,
  102.8265977671968,
  58.39070447933071,
  42.01451500795805,
  9.523770418462753,
  6.95631268600118,
  31.62766201396243,
  25.244079148727934,
  74.71889691396936,
  31.236141190981193,
  24.786967037988923,
  6.122398380105475,
  5.417140765627537,
  14.461384871439568,
  78.53648675725556,
  42.77447559141699,
  28.86931685524332,
  60.99998310461749,
  6.227087633799472,
  107.47634152240818,
  14.095113623373695,
  88.59071709566189,
  41.840726521834654,
  27.657580766254853,
  22.767794584542173,
  73.48899211355393,
  12.114093641176279,
  33.7663723096822,
  26.54422324200455,
  22.61116499362535,
  22.713810400142393,
  139.77593730435552,
  9.737580140139878,
  20.031078937212314,
  126.84831358045426,
  20.06291603638016,
  39.15799931658964,
  9.069319639681057,
  48.26116806569466,
  90.16281308540533,
  46.14274731649472,
  33.22964088097305,
  14.23269848518827,
  191.169578272674,
  108.14504140695574,
  7.704224694560673,
  14.771132032787731,
  237.37043882870378,
  5.75045442817992,
  11.050904857933716])

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)